Metazoan diversity in Chilean hypersaline lakes unveiled by environmental DNA

Saline and hypersaline wetlands are biodiversity hotspots for metazoans such as aquatic invertebrates and wading birds. However, the survival of these habitats and their biota is increasingly threatened by a combination of pressures from climate change and extractive processes. With the goal of improving conservation efforts in hypersaline ecosystems, this study tests the use of eDNA methods for metazoan biomonitoring. We employed a multi-assay approach utilizing three genetic markers—12S rRNA, 18S rRNA, and COI —to analyze biodiversity in sediment and water. Samples were collected from three hypersaline lakes in Northern Chile: Salar de Atacama (Laguna Puilar), Salar de Pujsa, and Salar de Tara. eDNA outputs were also compared with results generated from aquatic macroinvertebrate assessments using kick-nets. Our eDNA analyses revealed a total of 21 and 22 taxa across the three hypersaline lakes in sediment and water, respectively. The highest diversity was found in Salar de Tara (15 taxa within sediment and 13 taxa from water). Our multi-assay design detected a range of resident hypersaline taxa with different conservation status, spanning from rotifers to endangered snails, to amphipods and flamingos. Macroinvertebrate data via kick-net surveys further validated Salar de Tara as the most biodiverse system. Compared to net-based assessments, eDNA analysis allowed more refined taxonomic assignments for copepods and ostracods, while certain taxa such as Ephydridae or Hirudinea were not detected through molecular tests. Overall, this study provides evidence that eDNA is an effective tool to elucidate fine scale taxa assemblages and can refine conservation efforts in hypersaline lakes.

1 Introduction

Saline and hypersaline lakes worldwide provide immense ecological value (e.g., carbon sinks, microbial diversity hotspots) and host an array of metazoan communities, including charismatic birds such as flamingos and functionally pivotal groups of macroinvertebrate fauna (Balushkina et al., 2009; Hammer, 1986; Saccò et al., 2021b). For instance, the largest saltwater lake (4,400 km²) in the western hemisphere – the Great Salt Lake (Utah, USA), supports 338 migratory bird species which depend on the resident benthic invertebrates as an essential primary food source for their migrations along the Pacific Flyway (Baxter and Butler, 2020). In the southern hemisphere, the Salar de Atacama (3,000 km²), the third-largest salt-flat in the world after Salar de Uyuni in Bolivia (10,582 km²) and Salinas Grandes in Argentina (6,000 km²), hosts one vulnerable (Phoenicoparrus andinus) and two near threatened (Phoenicopterus chilensis and Phoenicoparrus jamesi) species of flamingos (Gajardo and Redón, 2019). Located in the pre-Andean depression of the Atacama Desert together with more than fifty other hypersaline lakes, Salar de Atacama has been the focus of extensive geohydrological studies (e.g., Finstad et al., 2016; Lowenstein et al., 2003; Marazuela et al., 2019), although with far less work being generated on ecological patterns.

In the Atacama region and worldwide, molecular-based studies on saline and hypersaline lakes have predominantly focused on microbes rather than metazoans. Microbial processes in these systems are essential to biogeochemical flows and ecosystem functioning, including nutrient cycling, carbon fixation, and the degradation of organic matter (Oren, 2015; Oren et al., 2009). This microbial focus has attracted significant industrial interest in regions with strong extractive drives due to the mineral resources found in these lakes (Garcés and Álvarez, 2020; Kidder et al., 2020). Consequently, there has been a diminished interest in studying macro-organismal groups. This trend is likely exacerbated by the lack of methodological and technological advancements, despite the important role metazoans play in energy transfer and habitat modification. As a result, our understanding of the functional dynamics regulating these ecosystems and the services they provide is limited. The generally low diversity and visibility of metazoans, combined with the significant logistical and sampling challenges of studying larger organisms, further diverts research attention away from metazoan taxa such as zooplankton, invertebrates, and birds.

Investigating metazoans in saline ecosystems is critical due to their functional role in the ecosystem, their capacity to act as bioindicators of environmental health and their unique biodiversity, which can reveal novel adaptations and evolutionary processes (Zadereev et al., 2020). Recent advancements in DNA sequencing technologies, bioinformatics, and sampling methodologies are facilitating a more comprehensive exploration of biodiversity, particularly within extreme or unusual environments such as hypersaline lakes or saltpans (Saccò et al., 2021b). One non-invasive technique gaining considerable momentum in aquatic biomonitoring studies is environmental DNA (eDNA) analysis – the detection of the total pool of DNA isolated from environmental samples (Takahashi et al., 2023 and references therein). Analysis of eDNA offers several logistical benefits over traditional methods (Blackman et al., 2024). Firstly, it enhances efficiency as eDNA sampling requires less time and effort compared to traditional biodiversity surveys, which often involve extensive fieldwork and specimen collection (Effendi et al., 2023). Secondly, eDNA analysis is cost-effective because its non-invasive nature reduces the need for specialized equipment and personnel, lowering the overall cost of biodiversity monitoring (Jo et al., 2022). Finally, eDNA analysis provides a broader scope by capturing a wide range of species, including those that are elusive or rare, thereby offering a more comprehensive understanding of the ecosystem’s biodiversity (Hansen et al., 2018). These three aspects have the potential to be highly beneficial for hypersaline lake assessments, understudied ecosystems often located in remote poly-extreme environments with potential risks for the health of researchers (due to high altitude and extreme temperature ranges, strong UV radiation, etc.) when conventional survey methods (e.g., trapping, counts) are utilized. However, the use of eDNA applications in studies of metazoan taxa in inland saline ecology has been limited to-date, mainly due to the technical challenges with high salt concentration and elevated UV radiation affecting the DNA extraction processes and eDNA persistence/quality in the environment, respectively (e.g., Saccò et al., 2021b; Tazi et al., 2014).

As demonstrated in other remote and extreme environments (such as deep oceans or groundwaters), eDNA has the potential to become a routinary tool for assessing biodiversity patterns due to recent analytical and technological progress (Hinz et al., 2022; Saccò et al., 2022; Takahashi et al., 2023). Building upon successful eDNA metabarcoding assessments of invertebrate community patterns in Australian hypersaline lakes (Campbell et al., 2023), this study aims to evaluate the applicability of eDNA metabarcoding techniques in delineating broader metazoan communities present in both water and sediment samples from hypersaline lakes in Chile. Specifically, we aim to (i) use an eDNA approach to describe the avian and invertebrate biodiversity present in three hypersaline lakes located at different altitudes in Northern Chile and to (ii) compare the macroinvertebrate biodiversity detected by eDNA against conventional net sampling and provide conservational guidelines. The results of this study bring new light on hypersaline diversity patterns and provide preliminary evidence on the validity and applicability of the eDNA technique in inland high salt conditions.

2 Materials and methods

2.1 Study area

Field work was carried out at three hypersaline lakes located in the Los Flamencos National Reserve, in Northern Chile. The reserve comprises a total of seven sectors from the Salar de Atacama to the salares and lagoons of the Altiplano in the Antofagasta region, covering an area of 740 km² (CONAF, 2008). Water, sediment and macroinvertebrate samples were collected in September 2021 from three sites: Salar de Tara (4322 masl), Salar de Pujsa (4514 masl) and Laguna Puilar (2306 masl). Samples were collected from three sampling points in each salar to encompass the environmental heterogeneity apparent in the sites (Figure 1). Geographical coordinates and physiochemical data for the three hypersaline lakes can be found in Supplementary Table 1.

Figure 1

Figure 1. Location of sampling points at Laguna Puilar (Salar de Atacama), Salar de Pujsa and Salar de Tara within the Los Flamencos National Reserve. The map depicts the difference in size of each salar and the three sampling points (in red) where samples were collected (Laguna Puilar displayed not in scale within the Salar the Atacama for representative reasons). At the bottom of the figure, some photos of the landscapes are shown: (A) Salar de Tara, (B) Salar de Pujsa and (C) Salar de Atacama (Laguna Puilar) (Credit: Gonzalo Salazar).

2.2 Kick-net sampling for macroinvertebrate diversity

Benthic invertebrate diversity was estimated using kick net sampling (mesh = 500 μm), a methodology that is frequently used for qualitative and semi-quantitative assessments in rivers and streams (Everall et al., 2017), but also in hypersaline lakes (Herbst, 1988). Three replicates per sampling point were collected, and each replicate involved passing the net through sediment/water interface across an area of approximate 10 m² and filtering 2000 L of water (totaling 18000 L of filtered water per lake). All macroinvertebrates found in the kick net were passed through a sieving system, consisting of four sieves of decreasing mesh size (500, 200, 100 and 44 μm) to remove excess sediment and eventually trap meiofauna (invertebrates smaller than 500 μm) from the sediments themselves. Samples were stored in plastic bags and then analyzed in the laboratory with a binocular stereo microscope. Taxa were identified to the lowest level possible by using taxonomic keys (e.g., Domínguez and Fernández, 2009; González, 2003; Roig-Juñent et al., 2014) and a table with the diversity data was generated.

2.3 eDNA sample collection

Water samples were collected using a 3 L sterilized bottle, and filtration was conducted in situ using a vacuum pump equipped with 0.4 and 0.22 µm mesh size filters. These filters were preserved in RNAlater within 2 mL tubes and stored at -20°C in darkness to maintain integrity. Although the initial goal was to obtain one 0.44 µm and one 0.22 µm filter for each sample, clogging issues due to the high salt concentrations in water occasionally necessitated the use of additional filters. For representative purposes, filters were combined to homogenize the volume analyzed per sampling point (for both the 0.44 and the 0.22 µm mesh size filters) which finally ranged between 200 to 237 ml (Supplementary Table 2). Additionally, one negative control for each sampling point, consisting of filtering 1 L of distilled water in the field, was included in the analysis. Sediment samples from the top 3 cm were collected in triplicate at each sampling point. These were stored in sterile 50 ml tubes containing analytical grade ethanol (99.9%) and kept frozen at –20°C until further analysis.

2.4 eDNA extraction and purification

All laboratory processes were conducted in dedicated laboratories within the Trace and Environmental DNA (TrEnD) Laboratory, Curtin University, Perth, Western Australia, and all controls were processed alongside samples through to sequencing. In total, 31 filters were used for the samples and 18 for the negative controls. DNA from the water was extracted using a modified Qiagen DNeasy blood and tissue kit protocol (Qiagen, Germany). The proportion of the filter papers from each sample replicate that provided the equivalent of approximately 225 mL of filtered water was digested (see Supplementary Table 2 for actual amounts). This quantity was chosen because it represented the minimum amount of water filtered at any one site. As all or a portion of some filters were digested for extraction, all of each of the sample control filters were digested alongside the sample filters. Two extraction controls were also included. Samples were digested using 810 µL Buffer ATL and 120 µL Proteinase K and were incubated overnight at 56°C with rotation. The following day 400 µL of the digest was combined with 400 µL Buffer AL and 400 µL 100% ethanol. The remainder of the protocol was as prescribed by Qiagen with the addition of a second wash with Buffer AW2 to remove any extra salts introduced by the increased ratio of Proteinase K. The final modification was elution with a double elution (2 x 50 µL) to maximize total DNA recovery for a combined extract volume of 100 µL in Buffer AE.

Sediment samples were extracted using a modified Power Soil Pro Kit protocol (Qiagen, Germany). Samples were defrosted and homogenized with two to four 10 mm silver balls which were added to each falcon tubes (depending on the volume) and placed in the TissueLyser (Qiagen, Germany) for 30 s at 30 Hz. Approximately 250 mg of the slurry was subsampled into Powerbead Tubes. Four small size, sterile glass beads were also added to all sample tubes and two extraction controls. After the addition of 800 µL of Buffer CD1 the samples were placed in the TissueLyser for two rounds at 25 Hz for 5 min (with inversion). Extractions were continued using a (QIAcube Qiagen, Germany) and following the Power Soil Pro IRT method. Extractions were eluted in a final volume of 100µL AE buffer.

2.5 eDNA metabarcoding

The eDNA in the samples was analyzed using three metabarcoding assays: 12S rRNA (L1091/H1985 - Cooper, 1994), 18S rRNA (18S_1F/18S_400R - Pochon et al., 2013), and COI (fwhF2/fwhR2n - Vamos et al., 2017) (Supplementary Table 3). To evaluate DNA extraction inhibition and to determine the required dilution for optimal amplification, PCR reactions were performed on each DNA extract by adding DNA template directly to the PCR master mix (neat), then performing serial dilutions (water to 1 in 10, and sediment to 1 in 10 and to 1 in 100). The PCRs were performed at a final volume of 25 µL where each reaction comprised of: 1× PCR Gold Buffer (Applied Biosystems), 0.25 mM dNTP mix (Astral Scientific, Australia), 2 mM MgCl₂ (Applied Biosystems), 1 U AmpliTaq Gold DNA polymerase (Applied Biosystems), 0.4 mg/mL bovine serum albumin (Fisher Biotec), 0.4 µM forward and reverse primers, 0.6 μl of a 1:10,000 solution of SYBR Green dye (Life Technologies), and 2 µL template DNA or dilution. PCRs were performed on StepOne Plus instruments (Applied Biosystems) with the following cycling conditions: 95°C for 5min, followed by 50 cycles of: 95°C for 30 s, Assay Tm °C (Supplementary Table 3) for 30 s, 72°C for 45 s, then a melt-curve analysis of: 95°C for 15 s, 60°C for 1 min, 95°C for 15 s, finishing with a final extension stage at 72°C for 10 min.

After selection and creation of the optimal dilution (neat, 1 in 10, or 1 in 100), PCRs were repeated in duplicate as described above with the addition of unique, single use combinations of 8 bp multiplex identifier-tagged (MID-tag) primers as described in Koziol et al. (2019) and van der Heyde et al. (2020). Master mixes were prepared using a QIAgility instrument (Qiagen) in an ultra-clean lab facility, with negative and positive PCR controls included on every plate. A sequencing library was created by combining samples into equimolar mini-pools based on the PCR amplification end points from each sample. The mini-pools were then analyzed using a QIAxcel fragment analyzer (Qiagen) to combine them into equimolar concentrations and form the final sequencing libraries. Each library was then size selected (Pippin Prep: Sage Sciences) with 2% dye-free cassettes, cleaned using a QIAquick PCR purification kit, quantified by Qubit (Thermo Fisher), and diluted to 2 nM. All libraries were sequenced on an Illumina MiSeq instrument using MiSeq V2 sequencing kits with custom sequencing primers.

2.6 Bioinformatics and data analysis

eDNA data provided a qualitative survey of the biodiversity hosted within the water and sediment targeted in the study. The diversity of the aquatic macroinvertebrate community was also assessed via a conventional semi-quantitative taxonomic approach, with the broader goal of comparing two distinctive analytical designs: eDNA (molecular) and kick-net (morphologically-based). Raw metabarcoding sequence data was filtered and analyzed using the eDNAFlow pipeline (Mousavi-Derazmahalleh et al., 2021) with custom filtering parameters applied (–minAlignLeng ‘12’, –minLen ‘(see Supplementary Table 3)’, –minsize ‘2’, perc_identity ‘95’ (‘90’ for 18S), –qcov ‘100’, –minMatch_lulu ‘97’). Resulting zero-radius operational taxonomic units (ZOTUs) were queried against the nucleotide database NCBI (GenBank) and assigned to the species level where possible or dropped back to the lowest common ancestor if multiple possible taxonomic assignments were given. Low abundant ZOTUs without at least one sample above 0.003% sequences assigned, were considered unreliable and excluded from the dataset (Vamos et al., 2017). Any ZOTUs from field, extraction and negative controls (a total of 1 field control per water sampling site, 4 extraction controls and 1 negative controls per assay) were removed from the dataset, and 3 different positive controls (one per assay) were also sequenced for testing purposes. To depict the diversity of taxa identified, a phylogenetic tree was constructed utilizing phyloT (Phylogenetic Tree Generator) relying on the taxonomy from the National Center for Biotechnology Information (NCBI), and subsequently rendered using the Interactive Tree of Life (iTOL) (Letunic and Bork, 2021).

3 Results

3.1 Biodiversity of the Chilean precordillera

Following read filtering, an average of 19,411 reads for 12S rRNA, 75,662 reads for 18S rRNA, and 4,858 reads for COI were achieved per sample. Across the three wetlands and two substrates sampled, 31 taxa were found via eDNA approaches. Salar de Tara exhibited the highest diversity in both sediment (15 taxa) and water (13 taxa), followed by Laguna Puilar (6 taxa in sediment and 11 taxa in water) and Salar de Pujsa (1 taxon in sediment and 4 taxa in water). The richness of taxa was comparable between water (n = 21) and sediment (n = 22) samples. For visual representation, see Figure 2 displaying the taxa distribution. Several groups were only recorded from sediments (Epirodrilus, Branchinectidae, Cyprideis, Limnocythere, Orthocladiinae, Trhypochthonius, Cyatholaimidae, Gieysztoria and Encentrum) or in water (Nais, Hyalella, Liposcelis, Smittia, Chroicocephalus, Scolopacidae, Spatula, Heleobia atacamensis, Macrostomum and Jakoba libera). The only taxon that was found in all the hypersaline lakes was the greater flamingo (Phoenicopterus). Arthropoda was the most diverse Phylum with 12 taxa within Hexanauplia (5 taxa), Ostracoda (4 taxa) and Insecta (3 taxa) followed by Chordata with Aves (6 taxa), while several other Phyla yielded one taxon (e.g., Tardigrada and Malacostraca).

Figure 2

Figure 2. iTOL graph showing the taxa found at the three hypersaline lakes (Salar de Tara, Laguna Puilar (salar de Atacama) and Salar de Pujsa) in water and sediment through the multi-assay (12S, 18S and COI) approach. On the left photos of some specimens from the field and sorting in the laboratory are shown (Credit: Gonzalo Salazar).

Metabarcoding analysis using the bird universal primer 12S yielded 6 bird species. The 18S dataset detected 18 taxa across a broader taxonomic spectrum and notably captured copepod diversity consistently across both water and sediment samples from each salar surveyed. Additionally, the 18S assay proved effective in detecting other microscopic taxa such as cestodes (Wardoides) and rotifers (Encentrum). Conversely, COI performed better than the other two assays at detecting gastropods (Heleobia atacamensis), dipterans (e.g., Smittia) and amphipods (Hyalella).

3.2 eDNA vs net sampling

Taxonomic assessment with eDNA and net sampling techniques revealed distinct patterns in the distribution and detection of macroinvertebrate taxa across the hypersaline lakes (Table 1). Within the Clitellata, Hirudinea was exclusively identified in Salar de Tara through net sampling. Oligochaeta, specifically the family Naididae, was detected in Laguna Puilar using net sampling methods. At a finer taxonomic resolution, the species Nais elinguis was discerned in Salar de Tara through metabarcoding techniques (COI). Interestingly, Naididae was exclusively detected via net sampling in Laguna Puilar, contrasting with the specific identification of Nais elinguis through eDNA metabarcoding in Salar de Tara. Among Gastropoda, Heleobia was identified in Laguna Puilar using both methodologies, with species-level identification achieved for Heleobia atacamensis through metabarcoding (COI). Crustaceans, notably Amphipoda such as Hyalella sp., were consistently found across all three ecosystems using both net sampling and metabarcoding approaches, although Hyalella sp. were not detected via eDNA analysis at Laguna Puilar. The copepod community included Calanoida (Boeckella sp.) observed in Salar de Tara via net sampling, Diaptomidae in Pujsa through metabarcoding, and Cyclopoida in Tara via metabarcoding. Harpacticoida species were observed in Tara and Puilar through net sampling, and in Salar de Pujsa through metabarcoding, with various copepod taxa (Nitokra sp., Amphiascoides sp., Mesochra sp., Itunella sp.) identified across multiple sites using metabarcoding techniques. Ostracoda species such as Podocopida sp. and Heterocypris sp. were detected in Salar de Pujsa and Salar de Atacama (Puilar) through metabarcoding, respectively, while Branchiopoda (Chydoridae and Ilyocryptidae) were detected in Salar de Tara through net sampling. Insecta included Hemiptera (Corixidae) found in Salar de Tara via net sampling, and Diptera families Ephydridae, Dolichopodidae, and Chironomidae identified using both sampling methods. Smittia sp. was detected in Laguna Puilar through metabarcoding.

Table 1

Table 1. Comparison of morphology-based taxonomic data obtained via kick-net surveying versus results through eDNA metabarcoding analysis for aquatic macroinvertebrates within the three hypersaline lakes targeted [Salar de Tara, Salar de Pujsa and Laguna Puilar (Salar de Atacama)].

4 Discussion

This study provides the first multi-assay eDNA analysis in water and sediment samples from Chilean high altitude hypersaline lakes. We successfully generated eDNA metabarcoding sequence data from both water and sediment samples, overcoming the inherent challenges associated with DNA preservation in these environments, which include high salinity, UV exposure, and difficulties in filtering saline water to collect eDNA. The findings from the eDNA assays revealed several endemic taxa (i.e., Heleobia atacamensis, Pheonicopterus) frequently reported in the hypersaline lakes of the Chilean Precordillera and Altiplano. These results extend recent studies indicating that eDNA approaches represent a promising and effective method for assessing the biodiversity of salt lakes (Campbell et al., 2023; Saccò et al., 2021a). All the taxa, apart from the ostracod Heterocypris, were only identified by a single assay, confirming the need for the implementation of multiple universal assays in eDNA studies attempting to examine metazoan community composition (Takahashi et al., 2023). At ecosystem level, our integrated eDNA-based findings demonstrate that the Salar de Tara hosts the highest metazoan diversity. These results underscore the necessity for enhanced biomonitoring activities (involving ideally extensive sampling campaigns over multiple seasons) to fully disclose the biodiversity patterns on site and maintain the ecohydrological integrity of the wetland (García-Sanz et al., 2021).

4.1 Avifaunal diversity

Six different waterbird taxa were identified in the eDNA data across three sites with Salar de Tara supporting all of them. Out of the three endemic species breeding in the Atacama Desert and Chilean Altiplano (Phoenicopterus chilensis, Phoenicoparrus jamesi and Phoenicoparrus andinus), only Chilean flamingos (Phoenicopterus sp.) - the most widespread and abundant overall (Gutiérrez et al., 2022 and references therein) - were identified in waters from all three sites. The reliable detection of Chilean flamingos using eDNA from water samples is important from a conservation perspective because flamingos are useful “flagship species” to promote the conservation of hypersaline lagoons which are under threat (Gajardo and Redón, 2019; Williams, 1998). In addition, Chilean flamingos are known for their complex movements and long migratory journeys during which they use different hypersaline habitats for feeding and breeding (Polla et al., 2018), and eDNA can help detect recent presence in the environment. Conversely, the smaller Phoenicoparrus specimens were abundant in all the three hypersaline lakes (personal observation Gonzalo Salazar, Figure 3), but the two species (Phoenicoparrus jamesi and Phoenicoparrus andinus) belonging to this genus were not detected via eDNA approaches. Stochastic sampling factors and reference sequence availability for Phoenicoparrus (e.g., no 12S sequences for this genus are currently available on GenBank) may have influenced the non-detection, stressing the need for development of species-specific assays that could even identify the sex of specimens as indicated by Chapman (2012).

Figure 3

Figure 3. Photos of Phoenicoparrus Jamesi individuals collected during fieldwork at Salar de Tara (A) and Salar de Pujsa (B). Credit: Gonzalo Salazar.

Another taxon of birds identified through aquatic eDNA is the South American crested duck (Lophonetta specularioides), whose distribution ranges in Chile include the Altiplano, a high-altitude plateau area situated among the high peaks of the Andes at over 4000 m a.s.l (Muñoz-Pedreros et al., 2015). Our molecular analyses revealed the presence of this species of native South American duck in waters of the Salar de Pujsa and Salar de Tara, the two wetlands located within the puna. Other avifaunal groups identified through aquatic eDNA and frequently reported in the region were gulls (Chroicocephalus) and sandpipers (Scalapacidae). Only two taxa were found in sediment (Phoenicopterus and Anas and both at Salar de Tara), suggesting that this substrate is not ideal for investigating the use of habitat of avifaunal species. Overall, these findings indicate that eDNA research has potential to build on traditional visual avifaunal surveys in hypersaline lakes, as already proved in waterholes wetlands (e.g., (Li et al., 2023) and for rare and elusive marsh bird species (e.g., Neice and Mcrae, 2021).

4.2 Invertebrate diversity

The majority of taxa identified in the waters and sediments of the three hypersaline lakes via eDNA approaches (n = 25) belonged to invertebrate groups. In sediment samples, most taxa were benthic, including nematodes (Cyatholaimidae), rotifers (Encentrum), oligochaetes (Epirodrilus), chironomids (Orthocladiinae), platyhelminthes (Gieysztoria) or potentially derived from dead tissue of species associated with the water column (ostracods Cyprodeis and Limnocythere). In water samples, our metabarcoding analyses identified the critically endangered gastropod species Heleobia atacamensis that is endemic to the Salar de Atacama. This snail was originally described in Tilopozo, South of the Salar de Atacama in 1860 (Collado et al., 2013). However, recent phylogenetic analyses and multi-locus research confirm that a much larger population exists (Collado et al., 2023) and reference therein), with specimens found at Peine, Tilomonte and Puilar wetlands, amongst others. Our results confirm the presence of the species at Laguna Puilar, providing evidence of the validity of eDNA approaches for hypersaline biomonitoring when reference sequence data is available. Our molecular-based design was also successful in bringing new light to the diversity of aquatic micro-animals (e.g., rotifers, tardigrades), amphipods (Hyalella), and certain meiofaunal groups such as copepods (harpacticoids at Laguna Puilar) and ostracods (Heterocypris at Tara and Puilar). Further, the eDNA data identified the cestode Wardoides that has not been reported in hypersaline lakes previously highlighting the potential capacity for eDNA methods to reveal otherwise overlooked biodiversity in these environments (Beng and Corlett, 2020). The only calanoid detected through molecular approaches was in Salar de Tara and was only able to be identified to the family Diaptomidae, a widespread and specious family comprising species that have adapted to a range of different aquatic environments including inland saline lakes (Bayly, 1993), but not reported in Northern Chile to date. Net surveys from the same salar identified the presence of the calanoid copepod genus Boeckella sp., a taxon that is known to occur in salt lakes from the Chilean Altiplano (Williams, 1998).

Overall, the eDNA approaches in our study generated greater taxonomic resolution compared to morphological identifications from net sampling in most faunal groups, particularly aquatic micro-animals (e.g., rotifers, tardigrades, copepods and ostracods) (Table 1). Nonetheless, the designed eDNA analysis was not able to detect some macroinvertebrate groups found with planktonic nets such as Ephydridae, Cladocera or Hirudinea, also with other groups being not recorded at sites when present (Hyalella at Laguna Puilar, Salazar’s personal communication). This aspect stresses the need for eDNA-net complementary approaches, and the development of reference DNA databases for hypersaline lakes to avoid false-negative results due to a lack of reference sequences or primer biases (Takahashi et al., 2023).

Indeed, the absence of the brine shrimps Artemia (Artemiidae) in this study is also noteworthy, as this taxon is typically considered a conspicuous component of southern South American salt lakes (Williams, 1998). Neither eDNA nor kick-net samples suggested the presence of Artemia in the three lakes studied, although an Anostracan of the family Branchinectidae was detected in sediment samples at Salar de Tara. These results could indicate that Artemia does not provide a major food source for flamingos in the three lakes targeted in this study, in line with Polla et al. (2018). Previous research also indicated that flamingos’ diet preferences are diverse and vary depending on seasonal food availability and species-specific trophic behaviors (Dorador et al., 2018; Gajardo and Redón, 2019), an eco-evolutionary pattern eventually contributing to the coexistence of sympatric species in similar waters (Polla et al., 2018). Nonetheless, given the relatively limited scope of this study, other stochastic factors such as sampling random effects or PCR primer binding issues might have played a role in shaping these findings presented. Another potential factor shaping our results relates to the fact that zooplankton and crustaceans often display dormancy patterns in hypersaline lakes (Moscatello and Belmonte, 2004; Shadrin et al., 2015) as a response of fluctuations in environmental conditions and resources availability (Shadrin and Anufriieva, 2020). However, our water and sediment samples were collected in early spring, which is when Artemia populations reach the maximum size in the Atacama region (Gajardo and Beardmore, 1993), indicating that a combination of factors is involved. Integrated water/sediment approaches are essential for elucidating seasonally-driven shifts in hypersaline (and freshwater) ecological community composition and dynamics (Schindler and Scheuerell, 2002). Moscatello et al. (2004) conducted a seasonal study along the Salento Peninsula in South East Italy, targeting the microplankton community both at active (via water column samples) and resting (via sediment samples) stages. The researchers identified more than 400 taxa - with 64 new records for the area - and also reported that per each period the most abundant taxa in water did not correspond to the dominant groups in sediment samples. Overall, these findings highlight the role of sediment as a reservoir of temporarily rare species (e.g., in form of cysts - see Moscatello and Belmonte, 2009), as this substrate provides an essential component for understanding the spatial and temporal dynamics shaping hypersaline biodiversity patterns. Our eDNA-based dataset reflects the uniqueness of sediment biodiversity at the targeted hypersaline lakes, as potentially resting groups of e.g., rotifers (Encentrum) or anostracans (Branchinectidae) were solely detected within this substrate. Recently, multi-substrate sampling designs have also been applied in hypersaline lakes through DNA metabarcoding techniques. Using a multi-marker approach, Zhang et al. (2024) targeted the broad biodiversity spectrum - from cyanobacteria (16S rRNA primer) to vertebrates (12S rRNA primer) - of Mang co Lake, located in the Tibetan Plateau. Insects dominated in abundance within water samples, while the invertebrate community in the sediment was more diverse and evenly distributed across several taxa, including Branchipoda, Flabellinia, and Insecta, among others. Our results - indicating the presence of branchiopods only in sediment samples - only partially matched the patterns observed. In addition, similarly to what our dataset indicated, several endemic waterbird species (e.g., the ruddy shelduck Tadorna ferruginea and the mallard Anas platyrhynchos) were detected at Mang co Lake, a finding that further confirms the potential of eDNA in revealing the use of hypersaline habitat by avifauna. Interestingly, Zhang et al. (2024) also reported the presence of Artemia urmiana (the prevalent arthropod in the Tibetan Plateau), confirming that the detection of endemic brine shrimp species via eDNA analyses in hypersaline lakes can be achieved. Hypersaline eukaryotic diversity was also investigated through DNA metabarcoding techniques (18S rRNA primer) in another Tibetan lake, Kyêbxang Co (He et al., 2022). Here, littoral, sublittoral, and profundal zones were sampled to comprehensively survey the resident benthic community. Out of a total of 42 species of metazoans, Chlorophyta, unclassified Eukaryota, Rotifera, Ciliophora, Arthropoda, Ascomycota, Haptista, Cercozoa, Nematoda, and Basidiomycota were the ten most abundant taxa. Significant differences in species composition among the three study zones were found, with littoral zones hosting higher eukaryotic diversity ranges than profundal and sublittoral zones. Given that salt lakes in the Atacama Altiplano are characterized by a high degree geo-hydrological heterogeneity among these three sub-environments (littoral, sublittoral and profundal) (Valero-Garcés et al., 2003), the molecularly-based sampling design applied by He et al. (2022) hosts great potential for fine scale ecological studies in the region. In another study employing 18S rRNA DNA metabarcoding, Gerasimova et al. (2023) focused on the protist community of 38 saline and hypersaline environments (located in South Ural and Crimea) characterized by a remarkable range of salinities (2–390‰). The differences in the taxonomic structure of the resident protist communities were driven by salinity, geographic location and pH, although the authors also acknowledged the presence of a large number of unclassified 18S rDNA sequences. Unfortunately, no protists were detected in our study, mainly due to the fact that - differently to Gerasimova et al. (2023) who employed a primer set previously tested and validated on protists (Gerasimova et al., 2020) - we used a more variable region of the 18s rRNA marker gene to target a broader inter-group diversity spectrum. Despite these promising preliminary results, further research on key variables – such as dormancy, seasonality, stratification patterns and fluctuations in salinity – is needed. In particular, the combined effect of these variables on the detectability of DNA fragments in the environment will be crucial to further validate the use eDNA analysis in unveiling the intricated biodiversity patterns of hypersaline lakes.

4.3 Conservational considerations

Our findings show that species of conservation concern (e.g., Chilean flamingos and H. atacamensis) are utilizing salt lakes in the Chilean precordillera and Altiplano. This is a crucial observation as saline lakes are increasingly subject to exploration to support mining of minerals such as lithium (Gajardo and Redón, 2019). A recent study targeting the Chilean portion of the so called “Lithium Triangle” (Salar de Atacama, Salar de Tara and Salar de Pujsa) by Gutiérrez et al. (2022) provided evidence of decreased flamingo diversity in correspondence of lakes subjected to increased lithium extraction over the last 30 years. This trend is likely to be indicative of overall ecological impoverishment of salars, an alarming dynamic that requires urgent targeted conservational management strategies to avoid environmental collapse (García-Sanz et al., 2021). Worldwide, the combination of uncontrolled resource extraction and climate change effects is threatening the long-term ecological integrity of hypersaline lakes. Inter-disciplinary approaches are the gateway to investigate such complex dynamics, an essential task in a world that is unequivocally going to face increased freshwater salinization due to anthropogenic impacts (Cunillera-Montcusí et al., 2022). Indeed, eDNA provides a powerful complementary platform for unveiling otherwise overlooked biodiversity patterns under high-salt conditions.

Data availability statement

All additional data are available in the Supplementary Material and the DNA-based data are archived in the Dryad repository (https://doi.org/10.5061/dryad.80gb5mm08).

Author contributions

MS: Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Validation, Writing – original draft, Writing – review & editing. MC: Data curation, Formal analysis, Investigation, Visualization, Writing – review & editing, Writing – original draft. PA: Investigation, Methodology, Visualization, Writing – review & editing. GS: Investigation, Writing – review & editing. TB: Data curation, Formal analysis, Investigation, Writing – review & editing. MH: Investigation, Writing – review & editing. AL: Writing – review & editing. NW: Writing – review & editing. CH: Investigation, Writing – review & editing. MA: Project administration, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research forms part of the eDNA for Global Biodiversity (eDGES) program.

Acknowledgments

This work was supported by resources provided by the Pawsey Supercomputing Centre with funding from the Australian Government and the Government of Western Australia. MS, GS, AL and NW acknowledge support from the BHP-Curtin alliance within the framework of the “eDNA for Global Environment Studies (eDGES)” program.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2025.1504666/full#supplementary-material

References

Balushkina E. V., Golubkov S. M., Golubkov M. S., Litvinchuk L. F., Shadrin N. V. (2009). Effect of abiotic and biotic factors on the structural and functional organization of the saline lake ecosystems in Crimea. Zh Obshch Biol. 70, 504–514.

PubMed Abstract | Google Scholar

Baxter B. K., Butler J. K. (2020). “Climate change and great salt lake,” in Great Salt Lake Biology: A Terminal Lake in a Time of Change. Eds. Baxter B. K., Butler J. K. (Springer International Publishing, Cham), 23–52. doi: 10.1007/978-3-030-40352-2_2

Crossref Full Text | Google Scholar

Bayly I. A. E. (1993). The fauna of athalassic saline waters in Australia and the Altiplano of South America: comparisons and historical perspectives. Hydrobiologia 267, 225–231. doi: 10.1007/BF00018804

Crossref Full Text | Google Scholar

Beng K. C., Corlett R. T. (2020). Applications of environmental DNA (eDNA) in ecology and conservation: opportunities, challenges and prospects. Biodivers Conserv. 29, 2089–2121. doi: 10.1007/s10531-020-01980-0

Crossref Full Text | Google Scholar

Blackman R., Couton M., Keck F., Kirschner D., Carraro L., Cereghetti E., et al. (2024). Environmental DNA: The next chapter. Mol. Ecol. 33, e17355. doi: 10.1111/mec.17355

PubMed Abstract | Crossref Full Text | Google Scholar

Campbell M. A., Laini A., White N. E., Allentoft M. E., Saccò M. (2023). When nets meet environmental DNA metabarcoding: integrative approach to unveil invertebrate community patterns of hypersaline lakes. J. Ocean. Limnol. 41, 1331–1340. doi: 10.1007/s00343-022-2151-9

Crossref Full Text | Google Scholar

Chapman A. (2012). Development of novel high-resolution melting (HRM) assays for gender identification of Caribbean Flamingo (Phoenicopterus ruber ruber) and other birds. (Corpus Christi, USA: The Texas A&M Texas University).

Google Scholar

Collado G. A., Torres-Díaz C., Vidal M. A., Valladares M. A. (2023). Genetic diversity, morphometric characterization, and conservation reassessment of the critically endangered freshwater snail, heleobia atacamensis, in the Atacama Saltpan, Northern Chile. Biology 12, 791. doi: 10.3390/biology12060791

PubMed Abstract | Crossref Full Text | Google Scholar

Collado G. A., Valladares M. A., Méndez M. A. (2013). Hidden diversity in spring snails from the Andean Altiplano, the second highest plateau on Earth, and the Atacama Desert, the driest place in the world. Zoological Stud. 52, 50. doi: 10.1186/1810-522X-52-50

Crossref Full Text | Google Scholar

CONAF (2008). Actualización Plan de Manejo Participativo Reserva Nacional Los Flamencos.

Google Scholar

Cooper A. (1994). DNA from museum specimens. In Ancient DNA: recovery and analysis of genetic material from paleontological,archaeological, museum, medical, and forensic specimens. (New York, NY: Springer New York). pp. 149–165.

Google Scholar

Cunillera-Montcusí D., Beklioğlu M., Cañedo-Argüelles M., Jeppesen E., Ptacnik R., Amorim C. A., et al. (2022). Freshwater salinisation: a research agenda for a saltier world. Trends Ecol. Evol. 37, 440–453. doi: 10.1016/j.tree.2021.12.005

PubMed Abstract | Crossref Full Text | Google Scholar

Domínguez E., Fernández H. R. (2009). “Macroinvertebrados bentónicos sudamericanos,” in Sistemática y biología (Fundación Miguel Lillo, Tucumán, Argentina), 656.

Google Scholar

Dorador C., Fink P., Hengst M., Icaza G., Villalobos A. S., Vejar D., et al. (2018). Microbial community composition and trophic role along a marked salinity gradient in Laguna Puilar, Salar de Atacama, Chile. Antonie Van Leeuwenhoek 111, 1361–1374. doi: 10.1007/s10482-018-1091-z

PubMed Abstract | Crossref Full Text | Google Scholar

Effendi H., Aprilia M., Hariyadi S., Permatasari P. A. (2023). Performance of environmental DNA metabarcoding to identify and quantify benthic organisms in river ecosystems. IOP Conf. Ser.: Earth Environ. Sci. 1266, 12079. doi: 10.1088/1755-1315/1266/1/012079

Crossref Full Text | Google Scholar

Everall N. C., Johnson M. F., Wood P., Farmer A., Wilby R. L., Measham N. (2017). Comparability of macroinvertebrate biomonitoring indices of river health derived from semi-quantitative and quantitative methodologies. Ecol. Indic. 78, 437–448. doi: 10.1016/j.ecolind.2017.03.040

Crossref Full Text | Google Scholar

Finstad K., Pfeiffer M., McNicol G., Barnes J., Demergasso C., Chong G., et al. (2016). Rates and geochemical processes of soil and salt crust formation in Salars of the Atacama Desert, Chile. Geoderma 284, 57–72. doi: 10.1016/j.geoderma.2016.08.020

Crossref Full Text | Google Scholar

Gajardo G. M., Beardmore J. A. (1993). Electrophoretic evidence suggests that the Artemia found in the Salar de Atacama, Chile, is A. franciscana Kellogg. Hydrobiologia 257, 65–71. doi: 10.1007/BF00005947

Crossref Full Text | Google Scholar

Gajardo G., Redón S. (2019). Andean hypersaline lakes in the Atacama Desert, northern Chile: Between lithium exploitation and unique biodiversity conservation. Conserv. Sci. Pract. 1, e94. doi: 10.1111/csp2.94

Crossref Full Text | Google Scholar

Garcés I., Álvarez G. (2020). Water mining and extractivism of the Salar de Atacama, Chile. WIT Trans. Ecol. Environ. 245, 189–199. doi: 10.2495/EID200181

Crossref Full Text | Google Scholar

García-Sanz I., Heine-Fuster I., Luque J. A., Pizarro H., Castillo R., Pailahual M., et al. (2021). Limnological response from high-altitude wetlands to the water supply in the Andean Altiplano. Sci. Rep. 11, 7681. doi: 10.1038/s41598-021-87162-6

PubMed Abstract | Crossref Full Text | Google Scholar

Gerasimova E. A., Balkin A. S., Filonchikova E. S., Mindolina Y. V., Zagumyonnyi D. G., Tikhonenkov D. V. (2023). Taxonomic structure of planktonic protist communities in saline and hypersaline continental waters revealed by metabarcoding. Water 15, 2008. doi: 10.3390/w15112008

Crossref Full Text | Google Scholar

Gerasimova E. A., Plotnikov A. O., Khlopko Y. A., Zlatogursky V. V. (2020). Multiple euryhaline lineages of centrohelids (Haptista: Centroplasthelida) in inland saline waters revealed with metabarcoding. J. Eukaryotic Microbiol. 67, 223–231. doi: 10.1111/jeu.12776

PubMed Abstract | Crossref Full Text | Google Scholar

González E. R. (2003). The freshwater amphipods Hyalella Smith 1874 in Chile (Crustacea: Amphipoda). Rev. Chil. Hist. Natural 76, 623–637. doi: 10.4067/S0716-078X2003000400007

Crossref Full Text | Google Scholar

Gutiérrez J. S., Moore J. N., Donnelly J. P., Dorador C., Navedo J. G., Senner N. R. (2022). Climate change and lithium mining influence flamingo abundance in the Lithium Triangle. Proc. R. Soc. B: Biol. Sci. 289, 20212388. doi: 10.1098/rspb.2021.2388

PubMed Abstract | Crossref Full Text | Google Scholar

Hammer U. T. (1986). Saline Lake Ecosystems of the World (Dordrecht, The Netherlands: Springer Science & Business Media).

Google Scholar

Hansen B. K., Bekkevold D., Clausen L. W., Nielsen E. E. (2018). The sceptical optimist: challenges and perspectives for the application of environmental DNA in marine fisheries. Fish Fisheries 19, 751–768. doi: 10.1111/faf.12286

Crossref Full Text | Google Scholar

He L., Wang Q., Wang Z., Wang F., Sun S., Liu X. (2022). Eukaryotic diversity based on high-throughput 18S rRNA sequencing and its relationship with environmental factors in a Salt Lake in tibet, China. Water 14, 2724. doi: 10.3390/w14172724

Crossref Full Text | Google Scholar

Herbst D. B. (1988). Comparative population ecology of Ephydra hians Say (Diptera: Ephydridae) at Mono Lake (California) and Abert Lake (Oregon). Hydrobiologia 158, 145–166. doi: 10.1007/BF00026273

Crossref Full Text | Google Scholar

Hinz S., Coston-Guarini J., Marnane M., Guarini J.-M. (2022). Evaluating eDNA for Use within Marine Environmental Impact Assessments. J. Mar. Sci. Eng. 10, 375. doi: 10.3390/jmse10030375

Crossref Full Text | Google Scholar

Jo T., Takao K., Minamoto T. (2022). Linking the state of environmental DNA to its application for biomonitoring and stock assessment: Targeting mitochondrial/nuclear genes, and different DNA fragment lengths and particle sizes. Environ. DNA 4, 271–283. doi: 10.1002/edn3.253

Crossref Full Text | Google Scholar

Kidder J. A., Leybourne M. I., Layton-Matthews D., Bowell R. J., Rissmann C. F. W. (2020). A review of hydrogeochemical mineral exploration in the Atacama Desert, Chile. Ore Geology Rev. 124, 103562. doi: 10.1016/j.oregeorev.2020.103562

Crossref Full Text | Google Scholar

Koziol A., Stat M., Simpson T., Jarman S., Di Battista J. D., Harvey E. S., et al. (2019). Environmental DNA metabarcoding studies are critically affected by substrate selection. Mol. Ecol. Resour. 19, 366–376. doi: 10.1111/1755-0998.12971

PubMed Abstract | Crossref Full Text | Google Scholar

Letunic I., Bork P. (2021). Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49, W293–W296. doi: 10.1093/nar/gkab301

PubMed Abstract | Crossref Full Text | Google Scholar

Li J., Seeber P., Axtner J., Crouthers R., Groenenberg M., Koehncke A., et al. (2023). Monitoring terrestrial wildlife by combining hybridization capture and metabarcoding data from waterhole environmental DNA. Biol. Conserv. 284, 110168. doi: 10.1016/j.biocon.2023.110168

Crossref Full Text | Google Scholar

Lowenstein T. K., Hein M. C., Bobst A. L., Jordan T. E., Ku T.-L., Luo S. (2003). An Assessment of stratigraphic completeness in climate-sensitive closed-basin lake sediments: Salar de Atacama, Chile. J. Sedimentary Res. 73, 91–104. doi: 10.1306/061002730091

Crossref Full Text | Google Scholar

Marazuela M. A., Vázquez-Suñé E., Ayora C., García-Gil A., Palma T. (2019). Hydrodynamics of salt flat basins: The Salar de Atacama example. Sci. Total Environ. 651, 668–683. doi: 10.1016/j.scitotenv.2018.09.190

PubMed Abstract | Crossref Full Text | Google Scholar

Moscatello S., Belmonte G. (2004). Active and resting stages of zooplankton and its seasonal evolution in a hypersaline temporary pond of the Mediterranean coast (the Vecchia Salina, SE Italy). Scientia Marina 68, 491–500. doi: 10.3989/scimar.2004.68n4491

Crossref Full Text | Google Scholar

Moscatello S., Belmonte G. (2009). Egg banks in hypersaline lakes of the South-East Europe. Saline Syst. 5, 1–7. doi: 10.1186/1746-1448-5-3

PubMed Abstract | Crossref Full Text | Google Scholar

Moscatello S., Rubino F., Saracino O. D., Fanelli G., Belmonte G., Boero F. (2004). Plankton biodiversity around the Salento Peninsula (South East Italy): an integrated water/sediment approach. Scientia Marina 68, 85–102. doi: 10.3989/scimar.2004.68s185

Crossref Full Text | Google Scholar

Mousavi-Derazmahalleh M., Stott A., Lines R., Peverley G., Nester G., Simpson T., et al. (2021). eDNAFlow, an automated, reproducible and scalable workflow for analysis of environmental DNA sequences exploiting Nextflow and Singularity. Mol. Ecol. Resour. 21, 1697–1704. doi: 10.1111/1755-0998.13356

PubMed Abstract | Crossref Full Text | Google Scholar

Muñoz-Pedreros A., De Los Ríos-Escalante P. R., Möller P. (2015). Zooplankton of the highland bogs of Putana, a desert wetland of the high puna, northern Chile. Crustaceana 88, 1235–1244. doi: 10.1163/15685403-00003482

Crossref Full Text | Google Scholar

Neice A., Mcrae S. (2021). An eDNA diagnostic test to detect a rare, secretive marsh bird. Global Ecol. Conserv. 27, e01529. doi: 10.1016/j.gecco.2021.e01529

Crossref Full Text | Google Scholar

Oren A. (2015). Halophilic microbial communities and their environments. Curr. Opin. Biotechnol. 33, 119–124. doi: 10.1016/j.copbio.2015.02.005

PubMed Abstract | Crossref Full Text | Google Scholar

Oren A., Baxter B., Weimer B. (2009). Microbial communities in salt lakes: Phylogenetic diversity, metabolic diversity, and in situ activities. Natural Resour. Environ. Issues 15, 51.

Google Scholar

Pochon X., Bott N. J., Smith K. F., Wood S. A. (2013). Evaluating detection limits of next-generation sequencing for the surveillance and monitoring of international marine pests. PLoS One 8, e73935. doi: 10.1371/journal.pone.0073935

PubMed Abstract | Crossref Full Text | Google Scholar

Polla W., Di Pasquale V. A., Rasuk M. C., Barberis I. M., Romano M., Manzo R. A., et al. (2018). Diet and feeding selectivity of the Andean Flamingo Phoenicoparrus andinus and Chilean Flamingo Phoenicopterus Chilensis in lowland wintering areas. Wildfowl 68, 2–29.

Google Scholar

Roig-Juñent S., Claps L. E., Morrone J. J. (2014). Biodiversidad de artrópodos argentinos (Mendoza: Editorial Sociedad Entomológica Argentina).

Google Scholar

Saccò M., Guzik M. T., van der Heyde M., Nevill P., Cooper S. J. B., Austin A. D., et al. (2022). eDNA in subterranean ecosystems: Applications, technical aspects, and future prospects. Sci. Total Environ. 820, 153223. doi: 10.1016/j.scitotenv.2022.153223

PubMed Abstract | Crossref Full Text | Google Scholar

Saccò M., White N. E., Campbell M., Allard S., Humphreys W. F., Pringle P., et al. (2021a). Metabarcoding under Brine: Microbial Ecology of Five Hypersaline Lakes at Rottnest Island (WA, Australia). Water 13, 1899. doi: 10.3390/w13141899

Crossref Full Text | Google Scholar

Saccò M., White N. E., Harrod C., Salazar G., Aguilar P., Cubillos C. F., et al. (2021b). Salt to conserve: a review on the ecology and preservation of hypersaline ecosystems. Biol. Rev. Camb Philos. Soc. 96, 2828–2850. doi: 10.1111/brv.12780

PubMed Abstract | Crossref Full Text | Google Scholar

Schindler D. E., Scheuerell M. D. (2002). Habitat coupling in lake ecosystems. Oikos 98, 177–189. doi: 10.1034/j.1600-0706.2002.980201.x

Crossref Full Text | Google Scholar

Shadrin N. V., Anufriieva E. V. (2020). Structure and trophic relations in hypersaline environments. Biol. Bull. Rev. 10, 48–56. doi: 10.1134/S2079086420010065

Crossref Full Text | Google Scholar

Shadrin N. V., Anufriieva E. V., Amat F., Eremin O. Y. (2015). Dormant stages of crustaceans as a mechanism of propagation in the extreme and unpredictable environment in the Crimean hypersaline lakes. Chin. J. Oceanology Limnology 33, 1362–1367. doi: 10.1007/s00343-015-4363-8

Crossref Full Text | Google Scholar

Takahashi M., Saccò M., Kestel J. H., Nester G., Campbell M. A., van der Heyde M., et al. (2023). Aquatic environmental DNA: A review of the macro-organismal biomonitoring revolution. Sci. Total Environ. 873, 162322. doi: 10.1016/j.scitotenv.2023.162322

PubMed Abstract | Crossref Full Text | Google Scholar

Tazi L., Breakwell D. P., Harker A. R., Crandall K. A. (2014). Life in extreme environments: microbial diversity in Great Salt Lake, Utah. Extremophiles 18, 525–535. doi: 10.1007/s00792-014-0637-x

PubMed Abstract | Crossref Full Text | Google Scholar

Valero-Garcés B. L., Delgado-Huertas A., Navas A., Edwards L., Schwalb A., Ratto N. (2003). Patterns of regional hydrological variability in central-southern Altiplano (18–26 S) lakes during the last 500 years. Palaeogeography Palaeoclimatology Palaeoecol. 194, 319–338. doi: 10.1016/S0031-0182(03)00284-0

Crossref Full Text | Google Scholar

Vamos E., Elbrecht V., Leese F. (2017). Short COI markers for freshwater macroinvertebrate metabarcoding. Metabarcoding Metagenomics 1, e14625. doi: 10.3897/mbmg.1.14625

Crossref Full Text | Google Scholar

van der Heyde M., Bunce M., Wardell-Johnson G., Fernandes K., White N. E., Nevill P. (2020). Testing multiple substrates for terrestrial biodiversity monitoring using environmental DNA metabarcoding. Mol. Ecol. Resour. 20, 732–745. doi: 10.1111/1755-0998.13148

PubMed Abstract | Crossref Full Text | Google Scholar

Williams W. D. (1998). Salinity as a determinant of the structure of biological communities in salt lakes. Hydrobiologia 381, 191–201. doi: 10.1023/A:1003287826503

Crossref Full Text | Google Scholar

Zadereev E., Lipka O., Karimov B., Krylenko M., Elias V., Pinto I. S., et al. (2020). Overview of past, current, and future ecosystem and biodiversity trends of inland saline lakes of Europe and Central Asia. Inland Waters 10, 438–452. doi: 10.1080/20442041.2020.1772034

Crossref Full Text | Google Scholar

Zhang S. Y., Yan Q., Zhao J., Liu Y., Yao M. (2024). Distinct multitrophic biodiversity composition and community organization in a freshwater lake and a hypersaline lake on the Tibetan Plateau. Iscience 27, 110124. doi: 10.1016/j.isci.2024.110124

PubMed Abstract | Crossref Full Text | Google Scholar