- Research
- Open access
- Published:
Differences in archaeal diversity and potential ecological functions between saline and hypersaline lakes on Qinghai-Tibet Plateau were driven by multiple environmental and non-environmental factors beyond the salinity
BMC Microbiology volume 24, Article number: 153 (2024)
Abstract
Background
Saline lakes are home to various archaea that play special and crucial roles in the global biogeochemical cycle. The Qinghai-Tibet Plateau hosts a large number of lakes with diverse salinity ranging from 0.1 to over 400Â g/L, harboring complex and diverse archaea. To the best of our knowledge, the formation mechanisms and potential ecological roles of archaea in Qinghai-Tibetan Plateau saline lakes remain largely unknown.
Results
Using High-throughput Illumina sequencing, we uncovered the vastly distinct archaea communities between two typical saline lakes with significant salinity differences on the Qinghai Tibet Plateau (Qinghai saline lake and Chaka hypersaline lake) and suggested archaea played different important roles in methanogenesis-related and nitrate reduction-related functions of these two lakes, respectively. Rather than the individual effect of salinity, the composite effect of salinity with diverse environmental parameters (e.g., temperature, chlorophyll a, total nitrogen, and total phosphorus) dominated the explanation of the variations in archaeal community structure in different habitats. Based on the network analysis, we further found the correlations between dominant archaeal OTUs were tight but significantly different between the two habitats, implying that archaeal interactions may also largely determine the shape of archaeal communities.
Conclusion
The present study improved our understanding of the structure and function of archaea in different saline lakes on the Qinghai-Tibet Plateau and provided a new perspective on the mechanisms underlying shaping their communities.
Background
Saline lakes are found all over the world and play important roles in global biogeochemical cycling [1, 2]. Although all three domains of life (Archaea, Bacteria, and Eukarya) could be found in saline lakes, archaea were thought to represent the unneglectable important component of the biological community in saline lakes, particularly in some hypersaline lakes (salinity > 50.0 g/L) [3,4,5,6]. For example, in the high-salinity sediments of hypersaline lake Chaka, the abundance of archaea was found two times higher than bacteria [5]. In the Great Salt Lake (salinity > 240.0 g/L) waters, more than 500 archaea sequences were detected using the clone-library method, which were attributed to 4 known genera (Haloplanus, Halorubrum, Natronococcus, Haloquadratum) and 5 unknown groups [6]. Comparatively, only a few bacteria sequences were detected in the Great Salt Lake at the same time, and almost all of them belonged to the bacterial genus Salinibacter. Furthermore, some previous studies discovered different taxa of archaea presented different metabolic characteristics, implying that they played varied ecological functions in the environment [7,8,9]. Thus, it is important to generate information on the archaea communities and ecological functions present in different saline lakes, as well as how they are influenced by their surroundings.
Salinity has frequently been assumed to be the primary important environmental factor influencing the composition of microorganisms in different saline lakes [5, 6]. For the archaea, a few studies have also revealed that salinity was a significant factor in the formation of their community in saline lakes, such as Gahai, Xiaochaidan, and Charhan lake [10, 11]. However, more and more recent evidences suggested the explanation of salinity on microorganisms’ community variations in different saline lakes was still very limited. For example, in saline lakes with different salinity on the Qinghai Tibet Plateau, it has been reported that salinity can only explain 14–33% of the bacterial community variations [12]. On one hand, these unexplained variations can be attributed to the unmeasured environmental factors [13, 14]. On the other hand, microbial interactions have recently been proposed to be an alternative explanation [15, 16]. However, up to our knowledge, how the environmental factors and microbial interactions influence the archaeal communities in different saline lakes still remains largely unexplored.
The Qinghai-Tibet Plateau hosts a large number of lakes with diverse salinity ranging from 0.1 to over 400 g/L, making it a natural laboratory and resource treasure house for studying archaea. However, although a few studies were used to study archaea diversity in different Qinghai-Tibet Plateau saline lakes [5, 11], the formation mechanisms and potential ecological roles of archaea in these lakes remain largely unknown. In this study, we performed 16 S rRNA gene phylogenetic analyses in two typical salinity lakes on the Qinghai-Tibet Plateau: Qinghai saline lake and Chaka hypersaline lake, which were nearly in graphic and altitude but significantly different in salinity. The ultimate goal of this study was to improve our understanding of the structure and potential ecological function of archaea in different plateau saline lakes, as well as to identify the dominant driving factors of the archaea communities.
Methods
Sample collection
Water samples were collected from two lakes in July 2022. Qinghai saline lake samples were marked as QW1, QW2, QW3, QW4, QW5, and QW6, and Chaka hyper-saline lake samples were marked as CW1, CW2, CW3, CW4, and CW5 (Table S1). Two liters of water were collected from each sampling site and stored in sterile containers to prevent contamination. One liter of water in individual stations was filtered through 0.22 μm polycarbonate Ispore membrane filters (Millipore, USA) and stored at -80 °C for sequencing analysis. The remaining water was analyzed for environmental variables within 24 h of collection.
Analysis of environmental variables
Dissolved oxygen (DO), pH, and salinity data were monitored on-site using YSI Pro Plus sensors (YSI, Yellow Springs, OH, USA). Ammonium (NH4+), nitrate (NO3−), and phosphate (PO43−) levels were measured using a QuAAtro39 Continuous Segment-ed Flow Analyzer (SEAL Analytical, Inc., Mequon, WI, USA). Chlorophyll a (Chl_a), total nitrogen (TN), and total phosphate (TP) were measured following the methods described by He [17]. The results are summarized in Table S1.
DNA extraction and amplicon sequencing
Total genomic DNA of environmental samples was extracted from 0.22-µm filters stored at -80 °C using an E.N.Z.A.™ Water DNA Kit (Omega Biotek, Inc., Norcross, GA, USA). The isolated DNA pellet was resuspended in 100 µL of Milli-Q water (Merck Millipore, Burlington, MA, USA), quantified using a Nanodrop Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), and stored at -20 °C until further use. To increase the recovery rate and specificity of archaeal sequences during high-throughput sequencing, two rounds of nested PCR were used in this study. For the first round, the archaeal-specific primers GU1ST-340 F (5′-CCCTAYGGGGYGCASCAG-3′) and GU1ST-1000R (5′-GGCCATGCACYWCYTCTC-3′) [18] were used within the following program: a denaturation step at 95 °C for 3 min, followed by five cycles consisting of 30 s at 94 °C, 20 s at 45 °C, 30 s at 65 °C; 20 cycles consisting of 20 s at 94 °C, 20 s at 55 °C, 30 s at 72 °C, a final extension at 72 °C for 5 min. 20–30 ng of the PCR product from the first round was used as template DNA for the second PCR, which was conducted using the same PCR conditions and the general V3–V4 primer sets 349 F (5′-GYGCASCAGKCGMGAAW-3′) and 806R (5′-GG′ACTACVSGGGTATCTAAT-3′) [19]. PCR products were examined by 2% (w/v) agarose gel electrophoresis in (Tris, Boric acid, TBE EDTA) buffer solution. Subsequently, the amplicons were loaded onto an Illumina HiSeq platform (Illumina Inc., San Diego, CA, USA) according to the manufacturer’s guidelines.
Bioinformatics analysis
The fastq file that resulted from the interpreted flowgrams was deposited in the NCBI Sequence Read Archive (BioProject number PRJNA983497). After sequencing, Cutadapt (v.1.18) [20] was used to remove the primer adapter and the read pairs were sutured (merged) into sequences using PEAR (v.0.9.8) [21], based on the overlapping relationship between PE read segments. FASTQ files were processed to generate separate FASTA and QUAL files, which were then analyzed using standard methods. PRINSEQ (v.0.20.4) [22] was used to excise bases with a mass value below 20 at the tail of the reads, and a window of 10 bp was set if the mean mass value within the window was below 20. The back-end bases were cut off from the window, sequences containing N and short sequences were filtered after quality control, and sequences with low complexity were filtered out. The FASTA files were de-replicated, abundance sorted, and singleton sequences were removed. The OTUs (using a 97% similarity cut-off) of archaea were clustered de novo using USEARCH (v.11.0.667) [23]. Mothur (1.43.0) [24] was used to determine the alpha diversity index. The Redundancy, Venn, Mantel-test, and network analysis of different samples were performed using R (v.4.3.0) software. Archaeal metabolism functional prediction was conducted using PICRUSt (v.1.1.4) [25, 26]. The ecological functional prediction of the archaeal communities was conducted using FAPROTAX (v1.2.1) [27].
Results
Differences in archaeal diversity, taxonomic compositions, and potential ecological functions
The overall number of archaeal sequences in Qinghai saline lake and Chaka hypersaline lake was 233,097 and 179,359, which were clustered into 791 and 428 operational taxonomic units (OTUs, at a 97% cut-off), respectively. As Table 1 shown, the richness (Chao index) of archaea in Qinghai saline lake (in an average of 532) was significantly higher (Wilcox.test, p = 0.004 < 0.05) than that in Chaka hypersaline lake (in an average of 352), suggesting more archaeal species were adaptable to the saline lake, but sensitive to the hypersaline environments. However, the alpha diversity (both Shannon and Simpson Index) in Chaka hypersaline lake showed no significant difference (Wilcox.test, p = 0.823 and 0.123 > 0.05) with that of Qinghai saline lake, demonstrating a stable and highly adapted archaea community existed in the hypersaline lake. Using Principal coordinate analysis (PCoA) with Bray-Curtis distance, we compared the beta diversity of archaeal communities in different lakes. Based on our results, the first two principal components explained 53.9% and 11.5% of the total variation, respectively. Furthermore, the PERMANOVA (p = 0.002 < 0.005) results and PCoA ordination plots indicated that there were significant differences in the composition of archaea between the two lakes (Fig. S1).
At the phyla level, the relative abundance of archaeal sequences in Qinghai saline lake was dominated by Woesearchaeota (43.94%), followed by Euryarchaeota (29.16%) and Thaumarchaeota (4.8%) (Fig. 1a). By contrast, nearly all the sequences (95.74%) in Chaka hypersaline lake were annotated on Euryarchaeota. Except for Euryarchaeota, a few archaeal sequences in lake Chaka were also annotated on Nanohaloarchaeota (2.51%), Woesearchaeota (0.4%), and Thaumarchaeota (0.06%) (Fig. 1a). Notable, although phylum Euryarchaeota was dominant in both two lakes, at the genus level, the Euryarchaeota sequences in lake Qinghai and Chaka were mainly annotated on Methanothrix/ Methanocalculus/ Methanomassiliicoccus and Halorubrum/ Halohasta/ Halonotius/ Natronomonas, respectively (Fig. 1b). Thus, these evidences revealed archaea with markedly different taxonomic compositions in Qinghai saline lake and Chaka hypersaline lake.
Using the PICRUSt software, we predicted the metabolic functions of archaea based on the COG database. In the Chaka hypersaline lake, archaeal metabolism was advantaged by ‘replication, recombination and repair’ and ‘inorganic ion transport metabolism’, which were some functions to tackle the extremely high-salinity environment. Comparatively, archaea in the Qinghai saline lake highly in ‘energy production and conversion’ and ‘translation, ribosomal structure and biogenesis’ (Fig. 2a). Furthermore, we utilized FAPROTAX to predict the ecological functions of the archaea communities in different habitats. Archaea in Qinghai saline lake showed a higher abundance of ‘methanogenesis’, ‘hydrogenotrophic methanogenesis’, and ‘methanogenesis by CO2 reduction with H2’ functions. In contrast, the archaea predominantly displayed ‘nitrate reduction’, ‘chemoheterotrophy’, and ‘aerobic chemoheterotrophy’ functionality in Chaka hypersaline lake (Fig. 2b).
Impact of the environmental factors on archaeal communities
Following the exclusion with high inflation factors (> 10) (such as pH, NO3− and DO), redundancy analysis (RDA) was used to understand the difference in the archaeal community structure of the two habitats and identify the significant environmental factors regulating their community structure. The RDA plot clearly demonstrated that the archaeal communities in the saline lake and hypersaline lake were substantially distinct (Fig. 3a). Salinity, Chlorophyll a (Chl_a), total nitrogen (TN), and temperature were identified as the most critical elements influencing the development of archaeal communities. Further variance partitioning analyses (VPA) confirmed that salinity largely impacted the LP community structure variation (Fig. 3b). But the contributions of Chl_a, temperature, and nutrient (TP + TN + PO43−) to the archaeal communities’ variation were also unneglectable. In fact, archaeal variation explained by salinity combined temperature (representing 34% of the variation), Chl_a (representing 26% of the variation), or nutrient (representing 21% of the variation) was largely higher than the variation single explained by the salinity (represent 13% of the variation).
The Mantel Test was used to examine the relationships between environmental conditions and the archaeal phyla. Consistent with prior findings, although salinity was found to significantly affect a lot of phyla, such as Euryarchaeota and Crenarchaeota, but it was not the only important factor affecting these phyla. Especially for the Euryarchaeota, which was one of the most important phyla to contribute to archaeal communities in both Chaka hypersaline and Qinghai saline lakes, DO, pH, and PO4− were also found to have a large effect on them (Fig. 4).
Archaeal interactions in different habitats and their response to environmental factors
Co-occurrence network analysis of the dominant archaeal OTUs (Top 50) was used to infer their interactions in different habitats. In Qinghai saline lake, the total nodes and correlation number of the network were 51 and 118, respectively (Fig. 5a). Surprisingly, the correlation between archaeal dominant OTUs in Qinghai saline lake was entirely positive, indicating the wide co-presence of the dominant archaea species. The number of the total nodes and correlation number of Chaka hypersaline lake were 53 and 120, respectively, which was not statistically different from those of Qinghai saline lake (Fig. 5b). However, some negative correlations between different dominant archaeal OTUs were found in Chaka hypersaline lake (∼ 8% of the total dominant archaeal OTUs’ correlations), representing the mutual exclusion existing in dominant archaeal OTUs under the extremely high salt environment.
Among the correlation relationships identified between dominant archaeal OTUs and environmental factors, there were also significant differences between the two habitats (Fig. 5). In the Qinghai saline lake, dominant archaeal OTUs were mainly correlated with salinity, total phosphorus (TP), and Chlorophyll a, but in the Chaka hypersaline lake, dominant archaeal OTUs were associated with DO, TN, pH, and PO43−. Notable, although salinity was observed with correlation with dominant archaeal OTUs in both two habitats, their correlations’ positive/negative relationship was different. Specifically, the relationships between salinity and dominant archaeal OTUs were all positive in the network of Qinghai saline lake, but all negative in the network of Chaka hyper-saline lake. Moreover, in each habitat, the importance of some nutrient factors was higher than that of salinity. Especially, the number of correlations between dominant archaeal OTUs and phosphorus (TP/PO43−) was equal to and exceeded that of salinity in the Qinghai and Chaka lakes, respectively.
Discussion
In this study, the results of high-throughput sequencing revealed that the predominant phylum within the archaeal communities in Chaka hypersaline lake was Euryarchaeota (Fig. 1a). This observation was consistent with those of previous studies on the dominant archaeal phylum in the Great Salt Lake [28] and the Dead Sea [29]. Genus-level annotation further revealed Euryarchaeota in Chaka hypersaline lake was dominated by Halorubrum, Halohasta, Halonotius and Natronomonas, which was consistent with the archaeal community composition in high-altitude Andean lakes [30]. Based on culturable and genomic analysis, previous studies uncovered many Natronomonas and Halorubrum strains (e.g. Natronomonas pharaonis DSM 2160, and Halorubrum lacusprofundi ATCC 49,239) with the ability to reduce NO3− and NO2− assimilatorily [31, 32]. Interestingly, nitrate-reduction function was also found abundant in the hypersaline lake (Fig. 2b). Thus, these lines of evidence suggested archaea played important role in the nitrogen cycling of Chaka hypersaline lake. In the Qinghai saline lake, a previous culturable-dependent study supposed Euryarchaeota and Woesearchaeota were the first and second dominant phyla [33], respectively, which was different from the results of this study. Clearly, the most likely explanation for that variance was the difference in study methods. As a relatively new member of the superphylum DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaea), Woesearchaeota was surprisingly diverse and abundant in a wide range of extreme environments, such as deep oil reservoir, oligo-trophic lakes and indicating a high diversity of their roles in global biogeochemical cycles [34, 35]. Recently, multivariate regression analysis further revealed that Woesearchaeota might function in consortium with methanogens in the cycling of carbon [36]. Interestingly, our results also detected the wide co-presence of Woesearchaeota and methanogenesis-functional Euryarchaeota (Fig. 5a), which supported the previous observation. Furthermore, in the Qinghai saline lake, methanogenesis-related functions and taxa were distinctive characteristics of the archaeal annotation. Especially, the dominant genera in Qinghai saline lake were Methanothrix (average relative abundance > 8.30%) (Fig. 1b), which was widely distributed in both natural and artificial anoxic environments and played a major role in global methane production [37, 38]. Previous research has shown that some methanogenic archaea are halotolerant and can function with a salinity up to 20 g/L [39]. Thus, the appropriate salinity level (in an average of 16.5) and other environmental factors of Qinghai saline lake may suitable for the growth of these methanogenic archaea and allow them to gain some advantage in the communities. Considering the importance of methane in the greenhouse effect [40], the contribution of Qinghai saline lake to the greenhouse effect is worth further evaluating.
In this study, both network, RDA and Mantel-test results suggested the single influence of salinity on the structuring of archaeal communities was limited, which distinguished from prior study results that focused on the variation of bacterioplankton communities in different saline lakes [12]. This could be owing to the fact that many archaea have the ability to ‘live with salt’ [9]. Instead, the composite effect of salinity with diverse environmental parameters (e.g., temperature, chlorophyll a, total nitrogen, and total phosphorus) dominated the explanation of the variations in archaeal community structure. In oceanic habitats, several previous studies have suggested the structure of archaeal communities was significantly influenced by temperature, nitrogen, and chlorophyll-a concentrations [41, 42]. While in some freshwater lakes, a growing body of evidence suggested that phosphate or total phosphorus can essentially alter the structure of archaeal communities [43, 44]. The influence of environmental factors on archaea in saline lakes seems to be a combination of the two habitats mentioned above, but the specific regulatory effects of these factors and their mechanisms still need further analysis, especially based on culture-dependent studies and in situ investigations.
Although a lot of environmental factors have been considered, nearly 40% of the constrained variance of archaeal communities in different habitats was still unexplained, highlighting the limitations of employing environmental factors to explain community shifts. Currently, it has been reported that community structure predictions considering biotic information are more accurate than those based on only environmental factors [15, 45]. In this study, more than 80% of the dominant OTUs in each habitat correlated with others and the average correlation number of each note was greater than 4, indicating the strong internal relationships between archaea in each habitat. However, quantifying the contribution of archaeal interactions to their community shaping still remains a big challenge. Moreover, the proportion of positive/negative correlations of dominant archaeal OTUs was significantly different between the two habitats. In particular, negatively correlated relationships of dominant OTUs can only be detected in Chaka hypersaline lake. In Tibetan Plateau lakes, previous studies also uncovered the proportion of bacterial negative correlations increased from 21% in saline lakes to 43% in hypersaline lakes and postulated that high nutrient fluxes in well-connected habitats support the expression of competitive advantage and lead to the exclusion of inferior species [12]. According to our findings, the above theory may also apply to archaea communities.
Conclusions
Using high-throughput sequencing, this study revealed the significantly different communities and potential ecological functions of archaea between Qinghai saline lake and Chaka hypersaline lake. Based on RDA, VPA, Mantel-test, and network analysis, we further suggested the differences of archaeal communities between saline and hypersaline lakes were driven by the combination of multiple environmental and non-environmental factors (archaeal interactions). Overall, this research improved our understanding of the structure and ecological role of archaea in saline lakes and provided an update on the mechanisms that shape their communities.
Data availability
Sequence data that support the findings of this study was deposited in the NCBI Sequence Read Archive (BioProject number PRJNA983497).
References
Saccò M, White NE, Harrod C, Salazar G, Aguilar P, Cubillos CF, Meredith K, Baxter BK, Oren A, Anufriieva E, et al. Salt to conserve: a review on the ecology and preservation of hypersaline ecosystems. Biol Rev. 2021;96(6):2828–50.
Liu Y, Priscu JC, Xiong J, Conrad R, Vick-Majors T, Chu H, Hou J. Salinity drives archaeal distribution patterns in high altitude lake sediments on the Tibetan Plateau. FEMS Microbiol Ecol. 2016;92(3):fiw033.
Oren A. Diversity of halophilic microorganisms: environments, phylogeny, physiology, and applications. J Ind Microbiol Biotechnol. 2002;28(1):56–63.
Margesin R, Schinner F. Potential of halotolerant and halophilic microorganisms for biotechnology. Extremophiles. 2001;5(2):73–83.
Jiang H, Dong H, Yu B, Liu X, Li Y, Ji S, Zhang CL. Microbial response to salinity change in Lake Chaka, a hypersaline lake on tibetan plateau. Environ Microbiol. 2007;9(10):2603–21.
Almeida-Dalmet S, Sikaroodi M, Gillevet PM, Litchfield CD, Baxter BK. Temporal study of the microbial diversity of the north arm of Great Salt Lake, Utah, U.S. Microorganisms. 2015; 3(3):310 – 26.
Baker BJ, De Anda V, Seitz KW, Dombrowski N, Santoro AE, Lloyd KG. Diversity, ecology and evolution of Archaea. Nat Microbiol. 2020;5(7):887–900.
Spang A, Caceres EF, Ettema TJG. Genomic exploration of the diversity, ecology, and evolution of the archaeal domain of life. Science. 2017;357(6351):eaaf3883.
Andrei AŞ, Banciu HL, Oren A. Living with salt: metabolic and phylogenetic diversity of archaea inhabiting saline ecosystems. FEMS Microbiol Lett. 2012;330(1):1–9.
Jiang H, Dong H, Deng S, Yu B, Huang Q, Wu Q. Response of archaeal community structure to environmental changes in lakes on the Tibetan Plateau, Northwestern China. Geomicrobiol J. 2009;26(4):289–97.
Zhong ZP, Liu Y, Miao LL, Wang F, Chu LM, Wang JL, Liu ZP. Prokaryotic community structure driven by salinity and ionic concentrations in plateau lakes of the Tibetan Plateau. Appl Environ Microb. 2016;82(6):1846–58.
Ji M, Kong W, Yue L, Wang J, Deng Y, Zhu L. Salinity reduces bacterial diversity, but increases network complexity in tibetan Plateau lakes. FEMS Microbiol Ecol. 2019;95(12):fiz190.
Stegen JC, Lin X, Fredrickson JK, Chen X, Kennedy DW, Murray CJ, Rockhold ML, Konopka A. Quantifying community assembly processes and identifying features that impose them. ISME J. 2013;7(11):2069–79.
Cremer J, Melbinger A, Wienand K, Henriquez T, Jung H, Frey E. Cooperation in microbial populations: theory and experimental model systems. J Mol Biol. 2019;431(23):4599–644.
Lima-Mendez G, Faust K, Henry N, Decelle J, Colin S, Carcillo F, Chaffron S, Ignacio-Espinosa JC, Roux S, Vincent F, et al. Determinants of community structure in the global plankton interactome. Science. 2015;348(6237):1262073.
Zhou J, Song X, Zhang C-Y, Chen G-F, Lao Y-M, Jin H, Cai Z-H. Distribution patterns of microbial community structure along a 7000-mile latitudinal transect from the Mediterranean Sea across the Atlantic Ocean to the Brazilian Coastal Sea. Microb Ecol. 2018;76(3):592–609.
He Y, Sen B, Shang J, He Y, Xie N, Zhang Y, Zhang J, Johnson ZI, Wang G. Seasonal influence of scallop culture on nutrient flux, bacterial pathogens and bacterioplankton diversity across estuaries off the Bohai Sea Coast of Northern China. Mar Pollut Bull. 2017;124(1):411–20.
He C, Lin W, Zheng X, Wang C, Hu Z, Wang W. Synergistic effect of magnetite and zero-valent iron on anaerobic degradation and methanogenesis of phenol. Bioresour Technol. 2019;291:121874.
Wang G, Li Q, Gao X, Wang XC. Synergetic promotion of syntrophic methane production from anaerobic digestion of complex organic wastes by biochar: performance and associated mechanisms. Bioresour Technol. 2018;250:812–20.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2.
Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina paired-end reAd mergeR. Bioinformatics. 2014;30(5):614–20.
Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27(6):863–4.
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microb. 2009;75(23):7537–41.
Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21.
Tatusov RL, Koonin EV, Lipman DJ. A genomic perspective on protein families. Science. 1997;278(5338):631–7.
Louca S, Parfrey LW, Doebeli M. Decoupling function and taxonomy in the global ocean microbiome. Science. 2016;353(6305):1272–77.
Tazi L, Breakwell DP, Harker AR, Crandall KA. Life in extreme environments: microbial diversity in Great Salt Lake. Utah Extremophiles. 2014;18(3):525–35.
Jacob JH, Hussein EI, Shakhatreh MAK, Cornelison CT. Microbial community analysis of the hypersaline water of the Dead Sea using high-throughput amplicon sequencing. MicrobiologyOpen. 2017;6(5):e00500.
Maldonado MJ, AlbarracÃn VH, Lara JA, Ferrero MA, FarÃas ME. Culture-dependent and -independent methods reveal dominance of halophilic Euryarchaeota in high-altitude Andean lakes. Aquat Microb Ecol. 2018;81(2):171–88.
Torregrosa-Crespo J, Bergaust L, Pire C, MartÃnez-Espinosa RM. Denitrifying haloarchaea: sources and sinks of nitrogenous gases. FEMS Microbiol Lett. 2017;365(3):fnx270.
Alcántara-Hernández RJ, Cs V-E, Zavala-DÃÂaz, de la Serna FJ, Rodriguez-Revilla J, Dendooven L, Marsch R. Haloarchaeal assimilatory nitrate-reducing communities from a saline alkaline soil. FEMS Microbiol Lett. 2009; 298(1):56–66.
Jiang H, Huang J, Yang J. Halotolerant and halophilic microbes and their environmental implications in saline and hypersaline lakes in Qinghai Province, China. In: Extremophiles in eurasian ecosystems: ecology, diversity, and applications. vol. 8; 2018: 299–316.
Liu X, Li M, Castelle CJ, Probst AJ, Zhou Z, Pan J, Liu Y, Banfield JF, Gu J-D. Insights into the ecology, evolution, and metabolism of the widespread woesearchaeotal lineages. Microbiome. 2018;6(1):1–16.
Ortiz-Alvarez R, Casamayor EO. High occurrence of Pacearchaeota and Woesearchaeota (Archaea Superphylum DPANN) in the surface waters of oligotrophic high-altitude lakes. Environ Microbiol Rep. 2016;8(2):210–7.
Liu X, Wang Y, Gu JD. Ecological distribution and potential roles of Woesearchaeota in anaerobic biogeochemical cycling unveiled by genomic analysis. Comput Struct Biotechnol J. 2021;19:794–800.
Zhou J, Smith JA, Li M, Holmes DE. Methane production by Methanothrix thermoacetophila via direct interspecies electron transfer with Geobacter metallireducens. mBio. 2023;14(4):e0036023.
Vigderovich H, Eckert W, Elvert M, Gafni A, Rubin-Blum M, Bergman O, Sivan O. Aerobic methanotrophy increases the net iron reduction in methanogenic lake sediments. Front Microbiol. 2023;14:1206414.
Gagliano MC, Sampara P, Plugge CM, Temmink H, Sudmalis D, Ziels RM, Atomi H. Functional insights of salinity stress-related pathways in metagenome-resolved Methanothrix genomes. Appl Environ Microb. 2022;88(10):e02449–21.
Camacho A, Picazo A, Rochera C, Santamans A, Morant D, Miralles-Lorenzo J, Castillo-Escrivà A. Methane emissions in Spanish saline lakes: current rates, temperature and salinity responses, and evolution under different climate change scenarios. Water. 2017;9(9):659.
Parada AE, Fuhrman JA. Marine archaeal dynamics and interactions with the microbial community over 5 years from surface to seafloor. ISME J. 2017;11(11):2510–25.
Quan Z-X, Kim J-G, Gwak J-H, Jung M-Y, An S-U, Hyun J-H, Kang S, Rhee S-K. Distinct temporal dynamics of planktonic archaeal and bacterial assemblages in the bays of the Yellow Sea. PLoS ONE. 2019;14(8):e0221408.
Yuan H, Zhang R, Li Q, Han Q, Lu Q, Wu J. Unveiling the ecological significance of phosphorus fractions in shaping bacterial and archaeal beta diversity in mesotrophic lakes. Front Microbiol. 2023;14:1279751.
Zhang W, Chen R, Meng F, Yuan H, Geng M, Cheng L, Yin H, Xue B, Wang J. Ecosystem functioning is linked to microbial evenness and community composition along depth gradient in a semiarid lake. Ecol Indic. 2021; 132(2021):108314.
Weiss S, Van Treuren W, Lozupone C, Faust K, Friedman J, Deng Y, Xia LC, Xu ZZ, Ursell L, Alm EJ, et al. Correlation detection strategies in microbial data sets vary widely in sensitivity and precision. ISME J. 2016;10(7):1669–81.
Acknowledgements
Not applicable.
Funding
This work was supported by the Research Program of Application Foundation of Qinghai Province (2024-ZJ-749), 2022 Planning Project of Qinghai Minzu University (2022GH01), The National Natural Science Foundation of China (NSFC, grant number 32170063), Innovation Program of Chinese Academy of Agricultural Sciences (2024-AEPI, CAAS-CSGLCA-202302), The Central Public-interest Scientific Institution Basal Research Fund (Agro-Environmental Protection Institute, Ministry of Agriculture and Rural Affairs).
Author information
Authors and Affiliations
Contributions
Conceptualization, M.B. and H.Y.; methodology, Y.W. and W.L.; software, Y.W.; validation, Y.W., M.B. and H.Y.; formal analysis, W.L. and G.B.; investigation, W.L. and G.B.; resources, Y.W. and W.L.; data curation, M.B. and H.Y.; writing, Y.W., W.L. and M.B.; funding acquisition, Y.W. and M.B. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not Applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wang, Y., Li, W., Bao, G. et al. Differences in archaeal diversity and potential ecological functions between saline and hypersaline lakes on Qinghai-Tibet Plateau were driven by multiple environmental and non-environmental factors beyond the salinity. BMC Microbiol 24, 153 (2024). https://doi.org/10.1186/s12866-024-03307-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12866-024-03307-3