Abstract
The nitrogen (N) cycle is the foundation of the biogeochemistry on Earth and plays a crucial role in global climate stability. It is one of the most important nutrient cycles in high-altitude lakes. The biogeochemistry of nitrogen is almost entirely dependent on redox reactions mediated by microorganisms. However, the nitrogen cycling of microbial communities in the high-altitude saline lakes of the Qinghai-Tibet Plateau (QTP), the world’s “third pole” has not been investigated extensively. In this study, we used a metagenomic approach to investigate the microbial communities in four high-altitude pristine saline lakes in the Altun mountain on the QTP. We observed that Proteobacteria, Bacteroidota, and Actinobacteriota were dominant in these lakes. We reconstructed 1,593 bacterial MAGs and 8 archaeal MAGs, 1,060 of which were found to contain nitrogen cycle related genes. Our analysis revealed that nitrite reduction, nitrogen fixation, and assimilatory nitrate reduction processes might be active in the lakes. Denitrification might be a major mechanism driving the potential nitrogen loss, while nitrification might be inactive. A wide variety of microorganisms in the lake, dominated by Proteobacteria, participate together in the nitrogen cycle. The prevalence of the dominant taxon Yoonia in these lakes may be attributed to its well-established nitrogen functions and the coupled proton dynamics. This study is the first to systematically investigate the structure and nitrogen function of the microbial community in the high-altitude pristine saline lakes in the Altun mountain on the QTP. As such, it contributes to a better comprehension of biogeochemistry of high-altitude saline lakes.
1 Introduction
Nitrogen is a vital element for organisms, critical for nucleic acid and protein synthesis, and is often the key element limiting the primary productivity of ecosystems (). Nitrogen exists in the environment in many chemical compounds with oxidation states ranging from +5 to −3. The processes by which nitrogen converts between its different chemical forms are described in the nitrogen cycle, which is the foundation of the biogeochemistry on Earth (). The biogeochemistry of nitrogen is almost entirely dependent on redox reactions mediated by microorganisms (). The microorganisms are involved in seven main processes of nitrogen cycling in ecosystems: nitrogen fixation, ammonification, DNRA (dissimilatory nitrate reduction), assimilation, nitrification, denitrification, and anammox (; ). The structure and function of microbial communities are closely related to the ecosystem function (; Zhang et al., 2023).
The global nitrogen cycle has been severely disturbed since the 20th century by human activities including burning fossil fuels and applying artificial nitrogen fertilizers (). Numerous environmental issues, including eutrophication of terrestrial and aquatic systems, an increase in the greenhouse gas nitrous oxide, and ozone depletion in the stratosphere, have been linked to the significant nitrogen input (; ). Over the past few decades, advancements in ecological studies have provided profound insights into the intricate pathways and mechanisms of nitrogen cycling in various terrestrial and aquatic environments, particularly in the soil and marine environment (; ; ; ). From unraveling the complexities of microbial nitrogen transformations to assessing the effects of human activities on nitrogen fluxes, recent research has significantly expanded our understanding of the complex networks of interactions among nitrogen-transforming microorganisms and their implications for the global biogeochemical nitrogen cycling. In the last 15 years, four new reactions have been discovered: oxidation of hydroxylamine to NO (), dismutation of NO to N2 and O2 (), hydrazine synthesis () and oxidation of hydrazine to N2 (). Numerous new metabolic processes were also discovered, including complete ammonia oxidation (comammox) to nitrate (; ) and phototrophic nitrite oxidation (). In addition, novel microorganisms, such as symbiotic heterotrophic nitrogen-fixing cyanobacteria () and ammonia-oxidizing archaea () were also identified. Nitrogen cycling in microbial communities is carried out by complex networks of microorganisms with great diversity. Studies of microbial communities in marine, agricultural and soil environments have found cyanobacteria (), Nitrospina spp. (), “Candidatus. Methanoperedens” spp. () and other significant nitrogen-functioning microbial species. However, although nitrogen cycling in these environments has been studied in some depth, little research has been done on nitrogen cycling in high-altitude lakes, especially high-altitude saline lakes in the pristine areas of the QTP. Although some studies have analyzed the taxa composition, functions and distribution of microorganisms involved in specific nitrogen functions in lakes on the QTP (; Yang et al., 2012, 2013, 2015), the majority of the studies focused on a single functional taxon and lacked an understanding of the microbial community’s role in the overall nitrogen cycle. Furthermore, the investigation locations were situated in places with lower altitude and higher levels of human activity, rather than high-altitude pristine areas. It was not until 2023 that the first Shotgun metagenomics study for the microbial community of a lake above 4,000 meters above sea level (lake Bamucuo on the QTP) was reported (Wei et al., 2023). However, because this study only included data from one lake, it lacked systematization and did not comprehensively analyze the nitrogen cycling function of the lake’s microbial community. There is yet to be a comprehensive and systematic research of the microbial communities involved in nitrogen cycling in several lakes at an altitude of more than 4,000 m on the Altun Mountain of the Tibetan Plateau.
The “third pole,” the Qinghai-Tibetan Plateau (QTP), is the world’s highest plateau (average altitude of 4,500 m) and the largest plateau (2 × 106 km2) (Yang et al., 2013). It contains thousands of lakes, which account for more than half of China’s total lake area (). The plateau lake ecosystem is an important component of the plateau ecosystem. The QTP, along with the Arctic and Antarctic regions, is one of the most sensitive and vulnerable ecosystems in the world to global climate change, as well as a significant regulator of regional and global climate (). However, despite being known as the “third pole,” there is far less data available for the QTP than for the Arctic and Antarctic (). Although many high-altitude lakes are located in areas of low human activity and are typically oligotrophic environments, in the context of a substantial increase in global nitrogen fluxes and an intensification of the greenhouse effect, and influenced by factors such as increased snow and ice melt, increased erosive power of rainfall, and atmospheric deposition, high-altitude lakes have continued to receive an increased amount of nitrogen from the environment (Wolfe et al., 2001; ; ; Zhang et al., 2019). On the QTP, data from the Lake Qinghai reveal that nitrogen from human activities enters the lake through atmospheric deposition, potentially impacting the aquatic ecosystem of the lake (Zhang et al., 2019). Four alpine lakes in the Sanjiangyuan Region (Lake Eling, Lake Longbao, Lake Sea Star, and Lake Zhaling) have shown an increasing tendency in numerous nitrogen and phosphorus indicators year after year (). Over the last three decades, Lake Hurleg in the northeastern portion of the QTP has experienced an increase in nitrogen and phosphorus concentrations, although receiving very little direct human disturbance (Wu et al., 2021). Therefore, understanding the nitrogen cycling function of microbial communities in high-altitude lakes on the QTP in the context of increasing global nitrogen fluxes can contribute to our understanding of global biogeochemical cycling and microbial community responses to climate change.
In this context, the objectives of this study were (i) to investigate the structure of microbial communities in four high-altitude pristine saline lakes in the Altun Mountains on the QTP and (ii) to reveal the characteristics of the microbial nitrogen cycle and the microbial taxa involved in the nitrogen cycle in these lakes.
2 Materials and methods
2.1 Environmental samples collection
The study area is located in the Altun Shan National Nature Reserve on the northern edge of the QTP in China. From 2019 to 2021, we sampled four high-altitude pristine saline lakes in the Altun mountain (Lake Aqqikkol, Lake Ayakkum, Lake Jingyu and Lake Wusuxiao) multiple times (Supplementary Figures S1–S5). The sampling period covered the entire year, including April, July, October and November. Therefore, the samples were able to provide a comprehensive picture of the state of the lake microbial communities across time (Supplementary Table S1). 6–10 L of surface lake water were collected and promptly filtered using sterile filter membranes. The water was initially filtered with 10 μm membranes to remove bulk contaminants, then with 3 μm and 0.22 μm membranes, sequentially. The acquired filter membranes were folded inward twice into a 90-degree fan form, coated in sterile tin foil, and stored in plastic ziplock bags. The plastic ziplock bags were kept in dry ice and shipped to the company (Magigen, China) for DNA extraction and sequencing. Alkaline potassium persulfate digestion UV spectrophotometric method (HJ 636-2012) was used to determine the total nitrogen of the water samples (). Environmental factors were measured using a YSI multiparameter digital water quality meter (YSI, America).
2.2 DNA extraction and metagenome sequencing
DNA from 30 lake water samples was extracted with the Advanced Water DNA Kit (ALFA-SEQ). Integrity and purity of extracted DNA were monitored on 1% agarose gels. PE150 strategy was performed on the Illumina HiSeq 2,500 high-throughput sequencing platform. The quality of raw sequencing data was evaluated by Illumina filters and Trimmomatic v0.36 (), which removed adapters, primers, and low-quality data (Zhao et al., 2023). The cleaned data were assessed for quality and basic information statistics with the SeqKit v0.9.3 ().
2.3 Metagenomic assembly and binning
The cleaned data were de novo assembled using MEGAHIT v1.2.9 () with the parameter --presets meta-sensitive (--min-count 1 --k-list 21,29,39,49,…,129,141). Contigs shorter than 2Kb were removed from the assembled metagenomic sequences. The assembled and filtered sequences were analyzed for quality and basic information statistics using the QUAST v5.0.2 (). Metagenomic binning was performed using MetaBAT v2.12.1 () and CONCOCT v0.5.0 (). The binning results produced by these two software programs were later integrated using DAS Tool v1.1.1 () to build non-redundant MAGs (Metagenome-Assembled Genomes). The completeness and contamination of MAGs were assessed using CheckM v1.1.3 with “lineage_wf” workflow (). For further analysis, only high- and medium-quality MAGs meeting the MIMAG standard (completeness ≥50% and contamination <10%) () were retained. Taxonomic classification of MAGs was conducted using GTDB-Tk (v2.1.1, database r207) (). Because the samples were collected from special high-altitude pristine saline lake water habitats, the majority of MAGs could not be identified to the species level in the GTDB database. At the species level, ANI between MAGs was calculated using FastANI v1.32 () to distinguish between species with a 95% threshold ().
2.4 Data analysis
2.4.1 Gene annotation and abundance calculation
Genes in assembled contigs or MAGs were predicted using Prodigal v2.6.3 () with the -p meta parameter. KAAS (KEGG automatic annotation server) () was used to perform functional annotation of predicted genes using KEGG (Kyoto Encyclopedia of Genes and Genomes) database (). For the genomic features of a specific taxon, we used all the MAGs associated with that taxon. A gene was considered to be present in the taxon if it appeared in more than 25% of the MAGs or in high-quality MAGs with greater than 90% completeness ().
To evaluate the abundance of a given gene in the samples, we first performed a HMMER v3.3.1 () search across contigs for both the gene of interest and the marker gene rpsB (ribosomal protein S2), assumed to be present in every bacterial genome. Then, we divided the total depth of the contigs containing the gene of interest by the total depth of the contigs containing the marker gene. This allowed us to interpret the gene abundance as the average copy number of that gene in the sample. Specifically, we calculated the gene abundance as:
where RPKMi () is the count of reads per kilobase million for contig i, Q is the set of contigs that contain the gene of interest, and C is the set of contigs that contain the marker gene.
2.4.2 Genome-wide phylogenetic tree construction
Only MAGs with ≥70% genomic completeness and < 7% contamination were retained to construct the phylogenetic tree with the GTDB marker genes (120 bacterial marker genes or 122 archaeal marker genes). Sequence alignment was first performed using MAFFT v7.453 (), then again using MUSCLE v5.0.1278 (), and finally the aligned sequences were trimmed using trimAL v1.4 (). Maximum likelihood (ML) tree building was performed using RAxML v8.2.11 () with the parameter -m PROTGAMMAAUTO -N autoMRE to automatically set the number of replications of tree building and to calculate the bootstraps. Phylogenetic trees were visualized using the iTOL (Interactive Tree Of Life) online tool ().
2.4.3 MAGs abundance calculation
Reads from cleaned sequencing data for each sample were recruited to MAGs using the “mem” module of BWA v0.7.17 (). The result files were formatted and sorted using SAMtools v1.7 (). The depth of each contig within each MAG for every sample, i.e., the count of reads mapped to each genomic locus, was calculated using the “jgi_summarise_bam_contig_depths” program contained in MetaBAT v2.12.1 (). The aggregate depth of all contigs within a sample for each MAG, termed as “MAG Raw Depth,” was derived by summing the depths of individual contigs. This “MAG Raw Depth” was then normalized for the size of the MAG and the size of the sample sequencing reads data. The formula is as follows:
The resulting depth values can reflect the abundance of MAGs in the samples. Furthermore, these normalized results can be compared across different MAGs and samples.
2.4.4 Species abundance in microbial communities
Species-level abundance was obtained using mOTUs v3 () based on unassembled metagenomic sequencing data. This information was then utilized to assess the structural composition of the microbial communities in the samples.
2.4.5 Statistical analysis
All statistical analyses were performed in the R programming language version 4.2 (). The significance of the difference in the abundance of genes across samples was assessed using the wilcox.test() function. The plots were created using the ggplot2 package ().
3 Results
3.1 The microbially mediated nitrogen cycle
The high-altitude pristine saline lakes in the Altun mountain are an oligonitrogenous environment, with a total nitrogen content of 0.2–0.5 mg/L. According to the metagenomic data, the nirB, nifH, and nasA genes were abundant in high-altitude pristine saline lakes in the Altun mountain, with median values of 4.27, 2.12, and 1.12, respectively (p < 0.05) (Figure 1; Supplementary Figure S6). Additionally, the nirA, nirK, and nirS genes also showed relatively high abundance, with median values of 0.33, 0.30, and 0.26, respectively (p < 0.05 compared to other genes, except for these three) (Figure 1; Supplementary Figure S6). Therefore, the nitrite reduction, nitrogen fixation, and assimilatory nitrate reduction processes might be active in the lakes. Especially the nitrite reduction process. The abundance of genes for the processes of dissimilatory nitrite reduction (nirB), assimilatory nitrite reduction (nirA) and nitrite reduction to nitric oxide (nirK/nirS) were all relatively high, implying the presence of numerous active nitrite reduction in the lake. Only the [FeMo]-nitrogenase gene (nif) was found in lake water samples; neither the [FeV]-nitrogenase gene (vnf) nor the [FeFe]-nitrogenase gene (anf) were discovered. In the four lakes, the depth profiles of nitrogen related genes showed similar patterns (Supplementary Figure S7). Among the six genes with high abundance, the nirB, nifH, nirK, and nirS genes showed no significant difference among the different lakes. The depth of the nasA gene was lower in Lake Ayakkum compared to Lake Aqqikkol and Lake Jingyu. Similarly, the depth of the nirA gene was lower in Lake Ayakkum than in Lake Aqqikkol (Supplementary Figure S8). This might suggest that dissimilatory nitrate reduction is relatively inactive in Lake Ayakkum compared to the other three lakes. In different sampling months (April, July, October, November), the nirB, nifH, nirK, and nirS genes also showed no significant difference. However, the depth of the nasA gene and the nirA gene was higher in October than in the other months (Supplementary Figure S9). These results indicated that the nirB, nifH, nirK, and nirS genes had a consistent and stable distribution in the lakes, whereas the nasA and nirA genes had a variable distribution across different lakes and varied in abundance over time.
Figure 1
3.2 Genome reconstruction and microbial community composition
A total of 1,975.8 GB of Illumina paired-end shotgun sequencing data was obtained from 30 lake water samples. The total size of de novo assembly (after removing contigs shorter than 2 KB) was 15.4 GB with average N50 7,618 bp (Supplementary Table S2). Binning the metagenomic data resulted in the reconstruction of 1,601 MAGs with quality higher than the medium quality standard (completeness ≥50% and contamination <10%) (). Among these MAGs, 1,593 were bacterial genomes and 8 were archaeal genomes (Figure 2; Supplementary Figure S10; Supplementary Table S3).
Figure 2
Of the 1,593 bacterial MAGs recovered from the lake water, a total of 24 phyla, 54 classes, 130 orders, 208 families, and 314 genera could be clearly categorized using the GTDB database. The ANI between MAGs was calculated, and a total of 701 species were obtained by dividing them with a 95% threshold. One MAG was classified as an unknown phylum, five MAGs were classified as unknown orders, 38 MAGs were classified as unknown families (2.4%), 376 MAGs were classified as unknown genera (23.6%), and 1,440 MAGs were classified as unknown species (90.4%). The 8 archaeal MAGs obtained from the lake water belonged to three phyla, three classes, three orders, and five families. They were all categorized as unknown species. Among them, two MAGs of the phylum Thermoplasmatota were classified as unknown families, and both were of good quality, suggesting the possibility of new taxonomic groups. The highest number of MAGs recovered from lake water samples was from Proteobacteria (696, including 377 Alphaproteobacteria, 318 Gammaproteobacteria, and 1 Zetaproteobacteria), followed by Bacteroidota (461) and Actinobacteriota (148) (Figure 2; Supplementary Figure S10; Supplementary Table S3). Based on the depth of the MAGs, the three most abundant phyla in the lakes were the same: Bacteroidota (33.1%), Proteobacteria (30.3%), and Actinobacteriota (18.8%) (Supplementary Figure S10; Supplementary Table S3). At the family level, Rhodobacteraceae (Alphaproteobacteria) had the most MAGs (201), followed by Flavobacteriaceae (Bacteroidota) (101), and Microbacteriaceae (Actinobacteriota) (81) (Supplementary Table S3). Based on the depth of the MAGs, the three most abundant families in the lakes were the same: Flavobacteriaceae (13.7%), Microbacteriaceae (13.5%), and Rhodobacteraceae (8.4%) (Supplementary Figure S10; Supplementary Table S3). In all four lakes, the highest number of MAGs was found in these three phyla. At the family level, the distribution of MAGs varied among the lakes (Table 1). Notably, Flavobacteriaceae and Rhodobacteraceae consistently ranked within the top five families in terms of MAG counts across all four lakes. Furthermore, Balneolaceae, Burkholderiaceae, Cyclobacteriaceae, Microbacteriaceae, and Saprospiraceae featured among the top five MAG-rich families in two lakes. All of the 14 genera present in all the four sampled lakes (Supplementary Table S4) encoded complete pathways related to low nitrogen in the two-component system (i.e., the NtrC family of nitrogen-regulation systems GlnL, GlnG, GlnA or NtrY, NtrX, NifA), whereas pathways associated with low phosphorus or low potassium were missing. Considering this, it could be hypothesized that functions related to nitrogen cycling may be one of the reasons why these taxa were able to prevail in these lakes. The classification results of MAGs suggested that microorganisms in high-altitude pristine saline lakes in the Altun mountain had a high species diversity, and a huge number of species had not yet been discovered and studied. Besides, nitrogen related functions played a significant role in the microbial communities of these lakes.
Table 1
| Lakes | Lake Aqqikkol (AQK) | Lake Ayakkum (AYK) | Lake Jingyu (JYH) | Lake Wusuxiao (WSX) | ||||
|---|---|---|---|---|---|---|---|---|
| No. | Family | Counts of MAGs | Family | Counts of MAGs | Family | Counts of MAGs | Family | Counts of MAGs |
| 1 | f__ | 17 | Rhodobacteraceae | 173 | Burkholderiaceae | 12 | Rhodobacteraceae | 12 |
| 2 | Planctomycetaceae | 8 | Flavobacteriaceae | 77 | Microbacteriaceae | 11 | Flavobacteriaceae | 8 |
| 3 | Cyclobacteriaceae | 8 | Microbacteriaceae | 63 | Flavobacteriaceae | 10 | Alteromonadaceae | 3 |
| 4 | Saprospiraceae | 7 | Saprospiraceae | 61 | Cyclobacteriaceae | 9 | Balneolaceae | 3 |
| 5 | Rhodobacteraceae | 7 | Burkholderiaceae | 55 | Rhodobacteraceae | 9 | f__ | 3 |
| 6 | Flavobacteriaceae | 6 | Balneolaceae | 50 | Pseudohongiellaceae | 7 | Parvularculaceae | 3 |
The top 5 families in terms of counts of MAGs in the water of four lakes.
3.3 MAGs contain nitrogen cycle related genes and their classification
A total of 1,060 MAGs in 24 phyla and 157 families were found to contain nitrogen cycle related genes (Figure 2; Supplementary Table S5). In all four lakes, the most abundant phyla, both in terms of number and depth of MAGs, were Proteobacteria, Bacteroidota, and Actinobacteriota. The six most abundant nitrogen circle related genes were also mainly encoded by these phyla. Specifically, the nirB gene was mainly encoded by Bacteroidota and Proteobacteria, the nifH gene by Proteobacteria, the nasA gene by Proteobacteria, the nirA gene by Bacteroidota, the nirK gene by Proteobacteria, Bacteroidota, and Actinobacteriota, and the nirS gene by Proteobacteria (Supplementary Table S5). The major phyla encoding the six most abundant nitrogen cycle related genes were broadly similar across each lake, though some differences were observed. In Lake Aqqikkol, the nasA gene was mainly encoded by the phyla Proteobacteria, Bacteroidota, and Actinobacteriota, while the nirA gene was mainly encoded by Verrucomicrobiota. In Lake Jingyu, the nirA gene was mainly encoded by Bacteroidota and Proteobacteria, and the MAG encoding the nirS gene was not annotated. In Lake Wusuxiao, the nirA gene was mainly encoded by Cyanobacteria (Supplementary Table S6). Besides, we found two archaeal MAGs that encoded nitrogen circle related genes (Figure 2; Supplementary Table S5). One belonged to Halobacteriota and encoded seven types of nitrogen circle related genes: napA, narG, narH, nirA, nirK, norB, and nosZ. The other one belonged to Thermoplasmatota and encoded napA and nirK genes. (Figure 2; Supplementary Table S5).
The genus Yoonia of Alphaproteobacteria, Rhodobacterales, Rhodobacteraceae was the most prevalent taxon in the majority of lake samples, accounting for 17 to 80% of the microbial community (Figure 3B). This genus was more abundant in Lake Aqqikkol and Lake Jingyu, and less abundant in Lake Ayakkum and Lake Wusuxiao (Figure 3B). Yoonia encoded a variety of nitrogen cycle functional genes, including the [FeMo]-nitrogenase nifH gene, the nirK gene, the narGHI genes, and the nasAB genes, etc. (Figure 4A). It contributed an average of 23.9% of the depth of the nirH gene, 47.2% of the depth of the nasA gene, and 23.9% of the depth of the nirA gene in the microbial communities of the lakes (Supplementary Figure S11). Therefore, Yoonia might play an important role in the nitrogen cycle of high-altitude pristine saline lakes. Bacteroidota mainly encoded nirB, nosZ, and nirK genes in the lake water (Figure 3A), and Cyclobacterium (Bacteroidota, Bacteroidia, Cytophagales, Cyclobacteriacea) was a typical genus of Bacteroidota that encoded all these three genes simultaneously (Figure 3C).
Figure 3
Figure 4
3.4 Nitrogen cycle related genes in Yoonia and Cyclobacterium
Yoonia was the genus with the greatest abundance in lake water samples. There were five mOTUs within the genus Yoonia, including four species. In lake water samples, two mOTUs of Y. vestfoldensis, [ref_mOTU_v3_10885] and [ref_mOTU_v3_06935], accounted for 99% of the overall abundance of Yoonia. Y. vestfoldensis strains were recovered from microbial mats in Antarctic lakes, indicating its preference for pristine habitats (). It has been reported that Yoonia dominates the microbial community in the polar ocean with a proportion of approximately 11.1% (Zeng et al., 2022; ). However, the strong dominance (17–80%) of Yoonia in the lakes of the Altun mountain has not been observed in any other environment, thus making Yoonia a unique geospecific taxon in this setting. A total of 27 MAGs of Yoonia were recovered. Of these 27 MAGs, 11 were classified as invalid species name Y. sp018402175 and the remainders could not be categorized at the species level in the GTDB database.
Yoonia had the ability to convert inorganic nitrogen into amino acids. In addition, it encoded cyanase (CynS), urease (UreABC) and nitrilase (NIT1), which could use cyanate, urea and nitrile as nitrogen sources. This might favor its survival in low-nitrogen environments. As regards the nitrogen transport system, Yoonia encoded the Amt family ammonium transporter (AMT), the NNP family nitrate/nitrite transporter (NRT) and the nitrate/nitrite/Cyanate ABC transport system (NrtABCD). Yoonia had the ability to efficiently take up nitrogen from its surroundings (Figure 4A).
Yoonia encoded membrane-bound respiratory nitrate reductases (NarGHI), which is a membrane protein complex consisting of three different subunits. The NarI subunit anchors the NarG and NarH subunits, which are situated in the cytoplasm, to the membrane. NarG and NarH form the reductase module, receiving electrons from NarI, an integral membrane protein that binds two b-type hemes and draws electrons from the menaquinol pool to NarH. Menaquinol is oxidized at the periplasmic side of NarI, where electrons are transferred from the haem on the P side to the haem on the N side, ultimately reducing nitrate at the N side of the cytoplasmic membrane. The electrogenic (or proton) motive of NarGHI is caused by the transmembrane charge separation. The nitrate reduction by NarGHI is coupled to electron input from formate via the nitrate-inducible formate dehydrogenase (Fdh-N), forming the full proton motive redox loop ().
Yoonia encoded the [FeMo]-nitrogenase (NifH), which uses eight electrons and 16 ATP to reduce each N2 molecule to two NH3 molecules (). Electrons and ATP can be obtained via coupling with oxidative phosphorylation. However, a greater challenge for nitrogen fixation is the management of oxygen. Since nitrogenase is oxygen-sensitive, aerobic nitrogen-fixing bacteria must employ techniques to protect nitrogenase from oxygen (). Yoonia is exposed to lower oxygen pressures because high altitudes are low-oxygen habitats. Yoonia encoded Fe-Mn family superoxide dismutase (SOD2), which may help protect its nitrogenase. Additionally, it has been found that aquatic heterotrophic organisms can fix nitrogen as long as there is a sufficient amount of carbohydrate substrate (). Yoonia’s ability to use a wide variety of carbon substrates, as well as its extensive transport systems, might also potentially be helpful for its nitrogen-fixing function.
Yoonia encoded not only a protein system associated with bacterial chemotaxis, namely: methyl-accepting chemotaxis protein (MCP), chemotaxis protein (Che), flagellar motor switch protein (Fli), and chemotaxis protein (Mot), but also four transmembrane chemoreceptors homologous to Escherichia coli, namely: serine sensor receptor (Tsr), aspartate sensor receptor (Tar), ribose and galactose sensor (Trg), and peptide sensor receptor (Tap) (). Yoonia also encoded genes for flagellar assembly, therefore it might be motile. Yoonia might be more capable of sensing chemicals and chemotaxis, allowing it have a greater ability to thrive in the oligonitrogenous lake environment.
A total of 13 MAGs of Cyclobacterium were recovered. Using the GTDB database, all of these MAGs could not be categorized at the species level. Except for organic nitrogen synthase like GlnA, GltBD and GdhA and nitrogen-containing compound hydrolase NIT1, as a typical genus of Bacteroidota, Cyclobacterium also encoded nirB, nrfA, nosZ, and nirK genes (Figure 4).
The soluble cytoplasmic NADH-dependent nitrite reductase (Nir) catalyzes the reduction of nitrite to ammonium (). Nir does not generate a proton gradient, but instead specifically uses NADH as the electron donor. NirB is the catalytic subunit. Cytochrome c-552 nitrite reductase (Nrf) is a soluble periplasmic enzyme. Different species may use NrfH protein or a complex consisting of NrfB, NrfC and NrfD proteins as the electron donor for NrfA. These are two genetically distinct systems, but both use NrfA as the catalytic enzyme (). Both Nir and Nrf catalyze dissimilatory nitrite reduction to ammonium, which is the critical reaction in the dissimilatory nitrate reduction to ammonium (DNRA). Cyclobacterium may grow by coupling DNRA with the oxidation of electron donors such as organic substances and sulfide (). Encoding multiple kinds of nitrite reductases may provide strains with enhanced environmental adaptation. Escherichia coli encodes both Nir and Nrf. Their expression is complementary: in environments with low nitrate or nitrite concentrations, E. coli synthesizes Nrf; in environments with high nitrate or nitrite concentrations, E. coli synthesizes Nir; and in environments with moderate nitrate concentrations, both enzymes are produced (). The coding and expression of these enzymes gives E. coli the ability to cope with a wide range of nitrate and nitrite concentrations in various settings. A similar pattern may exist in Cyclobacterium. NirK is located in the periplasm and use a c-type cytochrome as its electron source (). NosZ is also located in the periplasm and use a c-type cytochrome or accessory proteins as its electron source (Wunsch and Zumft, 2005). Neither of them directly contributes to energy conservation through the proton motive force.
Yoonia encoded nitrate reduction but not nitrite reduction, and Cyclobacterium did not encode nitrate reduction but does encoded nitrite reduction. Cyclobacterium also encoded nosZ which was not encoded by Yoonia. Yoonia and Cyclobacterium is a good example to show that the classical process of nitrogen transformation is seldom accomplished by one single microorganism. In microbial communities, the nitrogen cycling is usually engaged by a wide variety of microorganisms.
4 Discussion
Nitrogenase (Nif) genes was found in all four lakes, with a total of 204 MAGs encoding nitrogenase genes belonging to four phyla: Proteobacteria, Cyanobacteria, Myxococcota and Desulfobacterota (Supplementary Table S5). Proteobacteria accounted for the majority, with 197 MAGs that encoded nitrogenase being extensively dispersed among the four lakes. High abundance of the nifH gene was observed in high-altitude pristine saline lakes in the Altun Mountains on the QTP, implying active nitrogen fixation. These lakes are oligonitrogenous environments. In such environments, microbial nitrogen fixation is an important source of nitrogen input to the ecosystem (). High diversity of the nifH gene has been observed in oligonitrogenous environments like Antarctica and the Arctic (; ). Additionally, the nifH gene of Alphaproteobacteria has been detected in the oligotrophic waters of the Atlantic and Pacific Oceans and alpine lake lakes (Zehr et al., 1998; Wang et al., 2020; ). This is consistent with our findings that Alphaproteobacteria might assume a nitrogen-fixing function in oligonitrogenous aquatic environments. Rhodobacteraceae (Yoonia’s family, Alphaproteobacteria) was the only family with nitrogenase (Nif) genes annotated in all four lakes. Furthermore, Rhodobacteraceae was one of the most abundant taxa in the four lakes. Therefore, it is reasonable to hypothesize that Rhodobacteraceae is the main bearer of the nitrogen fixation function that is prevalent in the high-altitude pristine saline lakes in the Altun mountain. In addition to this family, each lake has its own specific nitrogen-fixing taxa.
In the denitrification pathway, nitrite reductase (NirK/NirS) and nitrous oxide reductase (NosZ) genes were more abundant, but nitric oxide reductase (NorBC) genes were less abundant (Figure 1). Similarly, a total of 115 MAGs encoding both subunits of nitric oxide reductase (NorBC) belonging to 6 phyla and 24 families were found, which was lower in number and variety compared to the taxa encoding nitrite reductase (NirK/NirS) (310 MAGs belonging to 17 phyla and 77 families) and nitrous oxide reductase (NosZ) (239 MAGs, belonging to 10 phyla and 10 families) (Supplementary Table S5). The low abundance and limited distribution of nitric oxide reductase genes (norBC) implied that the microbial community was weak in reducing nitrous oxide to nitrous oxide, whereas the high abundance and wide distribution of nitrous oxide reductase gene (nosZ) implied that the microbial community possessed a relatively high potential for reducing nitrous oxide. Therefore, the microbial community in high-altitude pristine saline lakes in the Altun mountain may prefer to release the end product nitrogen gas after completing the reduction of nitrous oxide in the process of denitrification.
The hydrazine synthase gene (hzs) was found to be less abundant in the lake than the nitrite reductase genes (nirK/nirS) (Figure 1; Supplementary Figure S12), and only four MAGs of Planctomycetota were found to encode the hzs gene (Supplementary Table S5). Therefore, anammox may not be active in the lakes, however denitrification may be a major mechanism driving the potential nitrogen loss in the high-altitude pristine saline lakes in the Altun mountain, which is consistent with the findings from high-altitude rivers on the QTP (Zhang et al., 2021). The abundance of the nirK gene in the lake water was similar to that of the nirS gene (Figure 1; Supplementary Figure S12), indicating that the denitrifier communities in the lake water were the mixed nirK/nirS type. Usually, nirK- and nirS-type denitrifier communities have different species compositions and respond differently to different environmental factors, which leads to changes in the species structure of the communities, and consequently affects the nitrogen loss process in the environment (; ). The vast majority of MAGs encoding the nitrite reductase gene (nirK/nirS) (303 MAGs, 97.7%) encoded only one type of nitrite reductase gene, with only a few MAGs encoding both nirK and nirS genes (7 MAGs, 2.3%, 4 Gammaproteobacteria, 2 Deinococcota, and 1 Acidobacteriae). Recent research has found that the two types of nir genes are not mutually exclusive (), however the functions of two different nitrite reductase types in one organism still need to be confirmed by further studies (). Interestingly, 27 MAGs encoding nirS/nirK genes were found in Rhodobacteraceae, of which 23 (belonging to four genera Yoonia, Roseovarius, Flavimaricola and Cypionkella) also encoded the nitrogenase gene (nif). These strains had the potential to both increase and decrease nitrogen in the environment.
Microbial communities in the high-altitude pristine saline lakes in the Altun mountain encoded all three types of nitrate reductase common in prokaryotes. A total of 341 MAGs encoding at least one type of nitrate reductase were annotated from all lake water samples, accounting for 21.3% of the total number of MAGs (1,601) and belonging to 17 phyla and 88 families. The majority of them encoded only one type of nitrate reductase, and there were 47 MAGs (13.78%) with more than one kind of nitrate reductase. The encoding of more than one types of nitrate reductase was observed mainly in Proteobacteria, especially in Gammaproteobacteria. Species encoding each nitrate reductase had considerable diversity, with significant variations in species composition among lakes. This implied that, whereas nitrate reduction to nitrite was a prevalent pathway in the lakes, different microbial taxa might achieve it in different lakes. It is important to note that the boundary between assimilatory and dissimilatory nitrate reduction is not absolute. Nitrite produced by assimilatory nitrate reduction can be further reduced in the respiratory chain. Experiments have demonstrated that Mycobacterium tuberculosis can use the NarGHI complex for assimilatory nitrate reduction ().
Notably, although nitrate reductase (NarGHI) catalyzes nitrate reduction to nitrite, and nitrite oxidoreductase (NxrAB) catalyzes nitrite oxidation to nitrate, narGHI and nxrAB genes are homologs (). In the KEGG database, narG and nxrA genes share the same KO ID, as do narH and nxrB genes. As a result, they are difficult to distinguish from each other through annotations. Therefore, we did not single out the nxr genes, and the results annotated to each of these KOs were counted as narGHI genes.
Microbial communities in the high-altitude pristine saline lakes in the Altun mountain also encoded all three types of nitrite reductase common in prokaryotes. A total of 206 MAGs encoding one or more nitrite reductase genes were annotated from all lake water samples, accounting for 12.87% of the total number of MAGs (1,601) and belonging to 18 phyla and 61 families. The majority of them encoded only one type of nitrite reductase, and there were only18 MAGs (8.74%) with more than one kind of nitrite reductase, including Cyclobacterium, which encoded two kinds. Similar to the microbial taxa encoding nitrate reductase, species encoding each nitrite reductase had considerable diversity, with significant variations in species composition among lakes. Although encoding multiple types of nitrite reductases may provide strains with adaptations to cope with the wide range of nitrate and nitrite concentrations in the environment (), it may be that because nitrate and nitrite concentrations in high-altitude pristine saline lakes remain stable and low for long periods of time, microorganisms living in the water do not have the need to cope with drastic changes in nitrate and nitrite concentrations, and therefore do not require the complementary action of multiple enzymes.
In all lake water samples, a total of 424 MAGs were annotated with the potential to participate in the process of assimilatory or dissimilatory nitrate reduction, i.e., encoding nitrate reductase and/or nitrite reductase. There were 123 MAGs encoding both nitrate lyase and nitrite reductase, 218 MAGs encoding only nitrate reductase, and 83 MAGs encoding only nitrite reductase. Even taking into account missing genes due to incomplete MAGs, there are still many microorganisms that in fact encode only one functional enzyme and complete only one step of the catalytic reaction in the lakes. Because of the metabolic diversity of microorganisms and the complicated network of interactions in their communities, the classical process of nitrogen transformation is seldom accomplished by one single microorganism. Microbial communities are the units that carry out the nitrogen cycle in the environment, and communities are able to recycle nitrogen more efficiently than individual microbes. In microbial communities, the wide variety of microorganisms involved in nitrogen transformation processes generate complicated nitrogen cycling networks. To achieve optimal growth, each microbe has unique physiological demands and performs distinct nitrogen cycling functions ().
The number of MAGs encoding the napA gene was detected to be only half of the number of MAGs encoding the narGHI genes, but the number of families to which they belonged was similar to the number of families to which MAGs encoding the narGHI genes belonged, and the number of phyla to which they belonged was even greater than that of the phyla to which MAGs encoding the narGHI genes belonged. It suggested that the distribution pattern of microorganisms encoding the napA gene in high-altitude pristine saline lakes in the Altun mountain was dispersed at the family level and above, with the napA gene encoded in microorganisms of many families but only a few strains of microorganisms in each family rather than universally. This distribution pattern could suggest the potential for the transfer of the napA gene across high-level taxonomic units, or gene acquisition and loss events that occurred relatively early in the evolutionary process. The abundance of the amo gene in the lake water was low and only one MAG encoding amo gene was found, belonging to Nitrosomonas in Nitrosomonadaceae of Gammaproteobacteria. It was implied that nitrification was not abundant in the lake water microbial community of high-altitude pristine saline lakes in the Altun mountain, although it might exist.
We noticed that for some genes, such as nirB, nifH, and nasA, the depth of these genes and the number and depth of the MAGs encoding these genes were consistent. However, for urease (ureABC genes), the gene depth was low, yet we found many MAGs encoding these genes (Figures 1, 3A). This might be due to binning biases. While powerful, the binning process in metagenomics is not without its biases. found that repeat sequences and regions with variant nucleotide compositions were frequently not binned, leading to biases in the sequence characteristics and gene content of MAGs. They also discovered that genes located in the non-binned regions were strongly biased towards ribosomal RNAs, transfer RNAs, mobile element functions, and genes of unknown function. Conversely, coding sequences with known functions are more likely to be binned. While it is unclear if this preference reflects differences in specific genes, it is plausible to speculate that microorganisms containing particular genes (e.g., the ureABC gene) are more likely to be binned. Although this bias may exist during binning, it has a greater impact on the annotation and abundance analyses of non-coding sequences and unknown genes in the microbial communities, and a lesser impact on the analyses of the distribution of known genes in different taxa. Therefore, combining gene abundance and functional annotation of MAGs to jointly determine the abundance of functional genes in microbial communities and their distribution across different taxa is a more reliable method.
High abundance of the nirB gene was observed in the lakes, along with high abundance of the nasA gene (Figure 1), indicating a high overall abundance of the nitrate reduction pathways. However, genes related to nitrite reduction were more abundant than those related to nitrate reduction. The abundance of nitrogen related genes in the environment is influenced by various environmental factors. Under low oxygen conditions, nitrite concentrations are high and nitrate concentrations are low, leading to increased activity and abundance of nitrite reductase (). Salinity may also affect the abundance of nitrogen related genes in the environment. For example, the abundance of nasA genes was significantly negatively correlated with environmental salinity (). Thus, the relative abundance of nitrate reduction related genes and nitrite reduction related genes varies across different environments. In some settings, nitrate reduction related genes are more abundant, while in others, nitrite reduction related genes are more abundant (; ; Zhang and Lv, 2021; ).
No nitrate-reducing activity has been observed in the nine species of the genus Yoonia that have been isolated and cultured, according to the List of Prokaryotic names with Standing in Nomenclature (LPSN)1 (). However, Yoonia was initially isolated from Loktanella in 2018 (Wirth and Whitman, 2018). Notably, nitrate-reducing activity has been observed in species of the genus Loktanella that are taxonomically related to Yoonia (). The genus Yoonia has been poorly studied, and its nitrite-reducing or nitrogen-fixing enzyme activities are unknown. Yoonia has a distinctive distribution, being more abundant in water bodies that are less disturbed by human activities, cold, at high altitudes, or at high latitudes (Zeng et al., 2021, 2022; ). This unique distribution makes it worthy of further in-depth study. All strains of Cyclobacterium obtained from our samples were unknown species. However, previous isolation and culture studies have shown that C. plantarum within this genus has positive nitrate and nitrite reduction, with gas formation from nitrate (). Nevertheless, nitrate reduction, but not nitrite reduction, was observed in several species, including C. amurskyense; C. qasimii; C. marinum; C. lianum, and C. caenipelagi (). Therefore, it was possible that the MAGs we obtained had nitrite-reducing activities but not nitrate-reducing activities, and further experimental verification is needed.
In conclusion, we revealed the structure and the nitrogen cycle of microbial communities in four high-altitude pristine saline lakes in the Altun mountain on the QTP. We observed that Proteobacteria, Bacteroidota, and Actinobacteriota were dominant in these lakes. We recovered 1,601 MAGs, of which 1,440 belong to unknown species, and 1,060 contained nitrogen cycle related genes. This significantly expanded the microbial genomic data resources. We found that the nitrite reduction, nitrogen fixation, and assimilatory nitrate reduction processes may be active in the lakes. Denitrification may be a major mechanism driving the potential nitrogen loss. The dominance of the taxon Yoonia in the lake could be attributed to its well-established nitrogen functions. This study contributes to our understanding of microbial diversity and nitrogen function in high-altitude saline lakes on the QTP.
Statements
Data availability statement
The original contributions presented in the study are publicly available. The data can be found here: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1121972/.
Author contributions
ZZ: Data curation, Formal analysis, Writing – original draft. YZ: Funding acquisition, Conceptualization, Methodology, Writing – review & editing. FM: Visualization, Writing – review & editing. MX: Investigation, Writing – review & editing. HL: Investigation, Resources, Writing – review & editing. JX: Investigation, Resources, Writing – review & editing. BH: Conceptualization, Supervision, Writing – review & editing. MW: Funding acquisition, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Science & Technology Basic Resources Investigation Program of China (Grant No. 2017FY100300).
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.
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/fmicb.2024.1424368/full#supplementary-material
Footnotes
References
1
AlnebergJ.BjarnasonB. S.de BruijnI.SchirmerM.QuickJ.IjazU. Z.et al. (2014). Binning metagenomic contigs by coverage and composition. Nat. Methods11, 1144–1146. doi: 10.1038/nmeth.3103
2
Bentzon-TiliaM.SeverinI.Hansen LarsH.RiemannL. (2015). Genomics and ecophysiology of heterotrophic nitrogen-fixing Bacteria isolated from estuarine surface water. MBio:6. doi: 10.1128/mbio.00929-15
3
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics30, 2114–2120. doi: 10.1093/bioinformatics/btu170
4
BottomleyP. J.TaylorA. E.MyroldD. D. (2012). A consideration of the relative contributions of different microbial subpopulations to the soil N cycle. Front. Microbiol.3:373. doi: 10.3389/fmicb.2012.00373
5
BowersR. M.KyrpidesN. C.StepanauskasR.Harmon-SmithM.DoudD.ReddyT. B. K.et al. (2017). Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat. Biotechnol.35, 725–731. doi: 10.1038/nbt.3893
6
BrahneyJ.MahowaldN.WardD. S.BallantyneA. P.NeffJ. C. (2015). Is atmospheric phosphorus pollution altering global alpine Lake stoichiometry?Global Biogeochem. Cycles29, 1369–1383. doi: 10.1002/2015gb005137
7
BrunetR.Garcia-GilL. (1996). Sulfide-induced dissimilatory nitrate reduction to ammonia in anaerobic freshwater sediments. FEMS Microbiol. Ecol.21, 131–138. doi: 10.1111/j.1574-6941.1996.tb00340.x
8
CanfieldD. E.GlazerA. N.FalkowskiP. G. (2010). The evolution and future of Earth’s nitrogen cycle. Science330, 192–196. doi: 10.1126/science.1186120
9
Capella-GutierrezS.Silla-MartinezJ. M.GabaldonT. (2009). trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics25, 1972–1973. doi: 10.1093/bioinformatics/btp348
10
CarantoJ. D.LancasterK. M. (2017). Nitric oxide is an obligate bacterial nitrification intermediate produced by hydroxylamine oxidoreductase. Proc. Natl. Acad. Sci. USA114, 8217–8222. doi: 10.1073/pnas.1704504114
11
CasciottiK. L. (2016). Nitrogen and oxygen isotopic studies of the marine nitrogen cycle. Annu. Rev. Mar. Sci.8, 379–407. doi: 10.1146/annurev-marine-010213-135052
12
ChaumeilP. A.MussigA. J.HugenholtzP.ParksD. H. (2022). GTDB-Tk v2: memory friendly classification with the genome taxonomy database. Bioinformatics38, 5315–5316. doi: 10.1093/bioinformatics/btac672
13
ChenX.WeiH.ZhangJ. E. (2021). Nitrogen and sulfur additions improved the diversity of nirK- and nirS-type denitrifying bacterial communities of farmland soil. Biology (Basel)10:1191. doi: 10.3390/biology10111191
14
ChurchM. J.ShortC. M.JenkinsB. D.KarlD. M.ZehrJ. P. (2005). Temporal patterns of nitrogenase gene (nifH) expression in the oligotrophic North Pacific Ocean. Appl. Environ. Microbiol.71, 5362–5370. doi: 10.1128/AEM.71.9.5362-5370.2005
15
DaimsH.LebedevaE. V.PjevacP.HanP.HerboldC.AlbertsenM.et al. (2015). Complete nitrification by Nitrospira bacteria. Nature528, 504–509. doi: 10.1038/nature16461
16
DanecekP.BonfieldJ. K.LiddleJ.MarshallJ.OhanV.PollardM. O.et al. (2021). Twelve years of SAMtools and BCFtools. Gigascience10:giab008. doi: 10.1093/gigascience/giab008
17
DíezB.BergmanB.Pedrós-AlióC.AntóM.SnoeijsP. (2012). High cyanobacterial nifH gene diversity in Arctic seawater and sea ice brine. Environ. Microbiol. Rep.4, 360–366. doi: 10.1111/j.1758-2229.2012.00343.x
18
DingW.WangS.QinP.FanS.SuX.CaiP.et al. (2023). Anaerobic thiosulfate oxidation by the Roseobacter group is prevalent in marine biofilms. Nat. Commun.14:2033. doi: 10.1038/s41467-023-37759-4
19
EddyS. R. (2011). Accelerated profile HMM searches. PLoS Comput. Biol.7:e1002195. doi: 10.1371/journal.pcbi.1002195
20
EdgarR. C. (2004). MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics5, 1–19. doi: 10.1186/1471-2105-5-113
21
EttwigK. F.ButlerM. K.Le PaslierD.PelletierE.MangenotS.KuypersM. M.et al. (2010). Nitrite-driven anaerobic methane oxidation by oxygenic bacteria. Nature464, 543–548. doi: 10.1038/nature08883
22
FalkowskiP. G. (1997). Evolution of the nitrogen cycle and its influence on the biological sequestration of CO2 in the ocean. Nature387, 272–275. doi: 10.1038/387272a0
23
FengX.XingP. (2023). Genomics of Yoonia sp. isolates (family Roseobacteraceae) from Lake Zhangnai on the Tibetan plateau. Microorganisms11:2817. doi: 10.3390/microorganisms11112817
24
GallowayJ. N.AberJ. D.ErismanJ. W.SeitzingerS. P.HowarthR. W.CowlingE. B.et al. (2003). The nitrogen cascade. Bioscience53, 341–356. doi: 10.1641/0006-3568(2003)053[0341:tnc]2.0.co;2
25
GrafD. R.JonesC. M.HallinS. (2014). Intergenomic comparisons highlight modularity of the denitrification pathway and underpin the importance of community structure for N2O emissions. PLoS One9:e114118. doi: 10.1371/journal.pone.0114118
26
GriffinB. M.SchottJ.SchinkB. (2007). Nitrite, an electron donor for anoxygenic photosynthesis. Science316:1870. doi: 10.1126/science.1139478
27
GruberN. (2008). The marine nitrogen cycle: overview and challenges. Nitrogen Mar. Environ.2, 1–50.
28
GruberN.GallowayJ. N. (2008). An earth-system perspective of the global nitrogen cycle. Nature451, 293–296. doi: 10.1038/nature06592
29
GurevichA.SavelievV.VyahhiN.TeslerG. (2013). QUAST: quality assessment tool for genome assemblies. Bioinformatics29, 1072–1075. doi: 10.1093/bioinformatics/btt086
30
HelenD.KimH.TytgatB.AnneW. (2016). Highly diverse nirK genes comprise two major clades that harbour ammonium-producing denitrifiers. BMC Genomics17:155. doi: 10.1186/s12864-016-2465-0
31
HengL. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM: arXiv preprint.
32
HuangY.FanZ.ZhaoC.ChenG.HuangJ.ZhouZ.et al. (2023). Evaluating the impacts of biochemical processes on nitrogen dynamics in a tide gate-controlled river flowing into the South China Sea. Sci. Total Environ.881:163363. doi: 10.1016/j.scitotenv.2023.163363
33
HutchinsD. A.CaponeD. C. (2022). The marine nitrogen cycle: new developments and global change. Nat. Rev. Microbiol.20, 401–414. doi: 10.1038/s41579-022-00687-z
34
HyattD.ChenG. L.LoCascioP. F.LandM. L.LarimerF. W.HauserL. J. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics11:119. doi: 10.1186/1471-2105-11-119
35
InomuraK.BraggJ.FollowsM. J. (2017). A quantitative analysis of the direct and indirect costs of nitrogen fixation: a model based on Azotobacter vinelandii. ISME J.11, 166–175. doi: 10.1038/ismej.2016.97
36
JainC.Rodriguez-RL. M.PhillippyA. M.KonstantinidisK. T.AluruS. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun.9:5114. doi: 10.1038/s41467-018-07641-9
37
JiangH.DongH.YuB.LvG.DengS.BerzinsN.et al. (2009). Diversity and abundance of Ammonia-oxidizing Archaea and Bacteria in Qinghai Lake, northwestern China. Geomicrobiol J.26, 199–211. doi: 10.1080/01490450902744004
38
JungY. T.LeeJ. S.YoonJ. H. (2013). Cyclobacterium caenipelagi sp. nov., isolated from a tidal flat sediment, and emended description of the genus Cyclobacterium. Int. J. Syst. Evol. Microbiol.63, 3158–3163. doi: 10.1099/ijs.0.050161-0
39
JungblutA. D.NeilanB. A. (2010). NifH gene diversity and expression in a microbial mat community on the McMurdo ice shelf, Antarctica. Antarct. Sci.22, 117–122. doi: 10.1017/S0954102009990514
40
KanehisaM.GotoS. (2000). KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res.28, 27–30. doi: 10.1093/nar/28.1.27
41
KangD. W. D.LiF.KirtonE.ThomasA.EganR.AnH.et al. (2019). MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ7:e7359. doi: 10.7717/peerj.7359
42
KartalB.MaalckeW. J.De AlmeidaN. M.CirpusI.GloerichJ.GeertsW.et al. (2011). Molecular mechanism of anaerobic ammonium oxidation. Nature479, 127–130. doi: 10.1038/nature10453
43
KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol.30, 772–780. doi: 10.1093/molbev/mst010
44
KongH.LinJ.ZhangY.LiC.XuC.ShenL.et al. (2023). High natural nitric oxide emissions from lakes on Tibetan plateau under rapid warming. Nat. Geosci.16, 474–477. doi: 10.1038/s41561-023-01200-8
45
KonnekeM.BernhardA. E.de la TorreJ. R.WalkerC. B.WaterburyJ. B.StahlD. A. (2005). Isolation of an autotrophic ammonia-oxidizing marine archaeon. Nature437, 543–546. doi: 10.1038/nature03911
46
KuangX.JiaoJ. J. (2016). Review on climate change on the Tibetan plateau during the last half century. J. Geophys. Res. Atmos.121, 3979–4007. doi: 10.1002/2015JD024728
47
KuypersM. M. M.MarchantH. K.KartalB. (2018). The microbial nitrogen-cycling network. Nat. Rev. Microbiol.16, 263–276. doi: 10.1038/nrmicro.2018.9
48
LaneM. C.LloydA. L.MarkyvechT. A.HaganE. C.MobleyH. L. T. (2006). Uropathogenic Escherichia coli strains generally lack functional Trg and tap chemoreceptors found in the majority of E. coli strains strictly residing in the gut. J. Bacteriol.188, 5618–5625. doi: 10.1128/jb.00449-06
49
LeeS. D. (2012). Loktanella tamlensis sp. nov., isolated from seawater. Int. J. Syst. Evol. Microbiol.62, 586–590. doi: 10.1099/ijs.0.029462-0
50
LetunicI.BorkP. (2007). Interactive tree of life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics23, 127–128. doi: 10.1093/bioinformatics/btl529
51
LiD. H.LiuC. M.LuoR. B.SadakaneK.LamT. W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics31, 1674–1676. doi: 10.1093/bioinformatics/btv033
52
LiB.WangL.LiH.XueJ.LuoW.XingP.et al. (2024). Phosphorus-driven regime shift from heterotrophic to autotrophic diazotrophs in a deep alpine lake. Water Res.248:120848. doi: 10.1016/j.watres.2023.120848
53
LinG.LinX. (2022). Bait input altered microbial community structure and increased greenhouse gases production in coastal wetland sediment. Water Res.218:118520. doi: 10.1016/j.watres.2022.118520
54
LiuQ.YangJ.WangB.LiuW.HuaZ.JiangH. (2022). Influence of salinity on the diversity and composition of carbohydrate metabolism, nitrogen and sulfur cycling genes in lake surface sediments. Front. Microbiol.13:1019010. doi: 10.3389/fmicb.2022.1019010
55
LockwoodC.ButtJ. N.ClarkeT. A.RichardsonD. J. (2011). Molecular interactions between multihaem cytochromes: probing the protein-protein interactions between pentahaem cytochromes of a nitrite reductase complex. Biochem. Soc. Trans.39, 263–268. doi: 10.1042/bst0390263
56
LuS.-J.SiJ.-H.HouC.-Y.LiY.-S.WangM.-M.YanX.-X.et al. (2017). Spatiotemporal distribution of nitrogen and phosphorus in alpine lakes in the Sanjiangyuan region of the Tibetan plateau. Water Sci. Technol.76, 396–412. doi: 10.2166/wst.2017.091
57
MacdonaldH.PopeN. R.ColeJ. A. (1985). Isolation, characterization and complementation analysis of nirB mutants of Escherichia Coli deficient only in NADH-dependent nitrite reductase activity. J. Gen. Microbiol.131, 2771–2782. doi: 10.1099/00221287-131-10-2771
58
MalmS.TiffertY.MicklinghoffJ.SchultzeS.JoostI.WeberI.et al. (2009). The roles of the nitrate reductase NarGHJI, the nitrite reductase NirBD and the response regulator GlnR in nitrate assimilation of Mycobacterium tuberculosis. Microbiology (Reading)155, 1332–1339. doi: 10.1099/mic.0.023275-0
59
Martínez-EspinosaR. M.ColeJ. A.RichardsonD. J.WatmoughN. J. (2011). Enzymology and ecology of the nitrogen cycle. Biochem. Soc. Trans.39, 175–178. doi: 10.1042/BST0390175
60
MilaneseA.MendeD. R.PaoliL.SalazarG.RuscheweyhH. J.CuencaM.et al. (2019). Microbial abundance, activity and population genomic profiling with mOTUs2. Nat. Commun.10:1014. doi: 10.1038/s41467-019-08844-4
61
MoriyaY.ItohM.OkudaS.YoshizawaA. C.KanehisaM. (2007). KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res.35, W182–W185. doi: 10.1093/nar/gkm321
62
MortazaviA.WilliamsB. A.McCueK.SchaefferL.WoldB. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods5, 621–628. doi: 10.1038/nmeth.1226
63
MosleyO. E.GiosE.CloseM.WeaverL.DaughneyC.HandleyK. M. (2022). Nitrogen cycling and microbial cooperation in the terrestrial subsurface. ISME J.16, 2561–2573. doi: 10.1038/s41396-022-01300-0
64
NelsonM. B.BerlemontR.MartinyA. C.MartinyJ. B. H. (2015). Nitrogen cycling potential of a grassland litter microbial community. Appl. Environ. Microbiol.81, 7012–7022. doi: 10.1128/AEM.02222-15
65
NelsonM. B.MartinyA. C.MartinyJ. B. (2016). Global biogeography of microbial nitrogen-cycling traits in soil. Proc. Natl. Acad. Sci. USA113, 8033–8040. doi: 10.1073/pnas.1601070113
66
NelsonW. C.TullyB. J.MobberleyJ. M. (2020). Biases in genome reconstruction from metagenomic data. PeerJ8:e10119. doi: 10.7717/peerj.10119
67
PachiadakiM. G.SintesE.BergauerK.BrownJ. M.RecordN. R.SwanB. K.et al. (2017). Major role of nitrite-oxidizing bacteria in dark ocean carbon fixation. Science358, 1046–1051. doi: 10.1126/science.aan8260
68
PajaresS.RamosR. (2019). Processes and microorganisms involved in the marine nitrogen cycle: knowledge and gaps. Front. Mar. Sci.6:739. doi: 10.3389/fmars.2019.00739
69
ParksD. H.ImelfortM.SkennertonC. T.HugenholtzP.TysonG. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res.25, 1043–1055. doi: 10.1101/gr.186072.114
70
ParteA. C.CarbasseJ. S.Meier-KolthoffJ. P.ReimerL. C.GokerM. (2020). List of prokaryotic names with standing in nomenclature (LPSN) moves to the DSMZ. Int. J. Syst. Evol. Microbiol.70, 5607–5612. doi: 10.1099/ijsem.0.004332
71
QiuJ. (2008). China: the third pole. Nature454, 393–396. doi: 10.1038/454393a
72
R Core Team. (2013). R: a language and environment for statistical computing.
73
RenM.ZhangZ.WangX.ZhouZ.ChenD.ZengH.et al. (2018). Diversity and contributions to nitrogen cycling and carbon fixation of soil salinity shaped microbial communities in tarim basin. Front. Microbiol.9:431. doi: 10.3389/fmicb.2018.00431
74
RichterM.Rossello-MoraR. (2009). Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA106, 19126–19131. doi: 10.1073/pnas.0906412106
75
ShahinpeiA.AmoozegarM. A.MirfeiziL.NikouM. M.VentosaA.Sánchez-PorroC. (2020). Taxogenomics of the genus Cyclobacterium: Cyclobacterium xiamenense and Cyclobacterium halophilum as synonyms and description of Cyclobacterium plantarum sp. nov. Microorganisms8:610. doi: 10.3390/microorganisms8040610
76
ShenW.LeS.LiY.HuF. Q. (2016). SeqKit: a cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PLoS One11:e0163962. doi: 10.1371/journal.pone.0163962
77
SieberC. M. K.ProbstA. J.SharrarA.ThomasB. C.HessM.TringeS. G.et al. (2018). Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy. Nat. Microbiol.3, 836–843. doi: 10.1038/s41564-018-0171-1
78
SimonJ.KlotzM. G. (2013). Diversity and evolution of bioenergetic systems involved in microbial nitrogen compound transformations. Biochim. Biophys. Acta1827, 114–135. doi: 10.1016/j.bbabio.2012.07.005
79
SimonJ.van SpanningR. J. M.RichardsonD. J. (2008). The organisation of proton motive and non-proton motive redox loops in prokaryotic respiratory systems. Biochim. Biophys. Acta1777, 1480–1490. doi: 10.1016/j.bbabio.2008.09.008
80
SohmJ. A.WebbE. A.CaponeD. G. (2011). Emerging patterns of marine nitrogen fixation. Nat. Rev. Microbiol.9, 499–508. doi: 10.1038/nrmicro2594
81
StamatakisA. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics30, 1312–1313. doi: 10.1093/bioinformatics/btu033
82
ThompsonA. W.FosterR. A.KrupkeA.CarterB. J.MusatN.VaulotD.et al. (2012). Unicellular cyanobacterium symbiotic with a single-celled eukaryotic alga. Science337, 1546–1550. doi: 10.1126/science.1222700
83
TuQ.HeZ.WuL.XueK.XieG.ChainP.et al. (2017). Metagenomic reconstruction of nitrogen cycling pathways in a CO2-enriched grassland ecosystem. Soil Biol. Biochem.106, 99–108. doi: 10.1016/j.soilbio.2016.12.017
84
VaksmaaA.Guerrero-CruzS.van AlenT. A.CremersG.EttwigK. F.LükeC.et al. (2017). Enrichment of anaerobic nitrate-dependent methanotrophic “Candidatus Methanoperedens nitroreducens” archaea from an Italian paddy field soil. Appl. Microbiol. Biotechnol.101, 7075–7084. doi: 10.1007/s00253-017-8416-0
85
Van KesselM. A.SpethD. R.AlbertsenM.NielsenP. H.Op den CampH. J.KartalB.et al. (2015). Complete nitrification by a single microorganism. Nature528, 555–559. doi: 10.1038/nature16459
86
Van TrappenS.MergaertJ.SwingsJ. (2004). Loktanella salsilacus gen. Nov., sp. nov., Loktanella fryxellensis sp. nov and Loktanella vestfoldensis sp. nov., new members of the Rhodobacter group, isolated from microbial mats in Antarctic lakes. Int. J. Syst. Evol. Microbiol.54, 1263–1269. doi: 10.1099/ijs.0.03006-0
87
VillanuevaR. A. M.ChenZ. J. (2019). ggplot2: elegant graphics for data analysis. Meas.: Interdiscip. Res. Pers.17, 160–167. doi: 10.1080/15366367.2019.1565254
88
VrecaP.MuriG. (2006). Changes in accumulation of organic matter and stable carbon and nitrogen isotopes in sediments of two Slovenian mountain lakes (Lake Ledvica and Lake Planina), induced by eutrophication changes. Limnol. Oceanogr.51, 781–790. doi: 10.4319/lo.2006.51.1_part_2.0781
89
WangH. N.GunsalusR. P. (2000). The nrfA and nirB nitrite reductase operons in Escherichia coli are expressed differently in response to nitrate than to nitrite. J. Bacteriol.182, 5813–5822. doi: 10.1128/jb.182.20.5813-5822.2000
90
WangM.WuJ.ZhouT.LiangY.ZhengL. X.SunY. X. (2021). Effects of copper and florfenicol on nirS- and nirK-type denitrifier communities and related antibiotic resistance in vegetable soils. Ecotoxicol. Environ. Saf.213:112011. doi: 10.1016/j.ecoenv.2021.112011
91
WangL.XingP.LiH.ZhouL.WuQ. L. (2020). Distinct intra-lake heterogeneity of diazotrophs in a deep oligotrophic mountain lake. Microb. Ecol.79, 840–852. doi: 10.1007/s00248-019-01461-0
92
WeiC.SunD.YuanW. L.LiL.DaiC. X.ChenZ. Z.et al. (2023). Metagenomics revealing molecular profiles of microbial community structure and metabolic capacity in Bamucuo lake, Tibet. Environ. Res.217:114847. doi: 10.1016/j.envres.2022.114847
93
WirthJ. S.WhitmanW. B. (2018). Phylogenomic analyses of a clade within the roseobacter group suggest taxonomic reassignments of species of the genera Aestuariivita, Citreicella, Loktanella, Nautella, Pelagibaca, Ruegeria, Thalassobius, Thiobacimonas and Tropicibacter, and the proposal of six novel genera. Int. J. Syst. Evol. Microbiol.68, 2393–2411. doi: 10.1099/ijsem.0.002833
94
WolfeA. P.BaronJ. S.CornettR. J. (2001). Anthropogenic nitrogen deposition induces rapid ecological changes in alpine lakes of the Colorado front range (USA). J. Paleolimnol.25, 1–7. doi: 10.1023/a:1008129509322
95
WuY.WangS.NiZ.LiH.MayL.PuJ. (2021). Emerging water pollution in the world’s least disturbed lakes on Qinghai-Tibetan plateau. Environ. Pollut.272:116032. doi: 10.1016/j.envpol.2020.116032
96
WunschP.ZumftW. G. (2005). Functional domains of NosR, a novel transmembrane iron-sulfur flavoprotein necessary for nitrous oxide respiration. J. Bacteriol.187, 1992–2001. doi: 10.1128/jb.187.6.1992-2001.2005
97
YangJ.JiangH.DongH.HouW.LiG.WuG. (2015). Sedimentary archaeal amoA gene abundance reflects historic nutrient level and salinity fluctuations in Qinghai Lake, Tibetan plateau. Sci. Rep.5:18071. doi: 10.1038/srep18071
98
YangJ.JiangH.DongH.WangH.WuG.HouW.et al. (2013). amoA-encoding archaea and thaumarchaeol in the lakes on the northeastern Qinghai-Tibetan plateau, China. Front. Microbiol.4:329. doi: 10.3389/fmicb.2013.00329
99
YangJ.JiangH.WuG.HouW.SunY.LaiZ.et al. (2012). Co-occurrence of nitrite-dependent anaerobic methane oxidizing and anaerobic ammonia oxidizing bacteria in two Qinghai-Tibetan saline lakes. Front. Earth Sci.6, 383–391. doi: 10.1007/s11707-012-0336-9
100
ZehrJ. P.MellonM. T.ZaniS. (1998). New nitrogen-fixing microorganisms detected in oligotrophic oceans by amplification of Nitrogenase (nifH) genes. Appl. Environ. Microbiol.64, 3444–3450. doi: 10.1128/AEM.64.9.3444-3450.1998
101
ZengY. X.LiH. R.LuoW. (2022). Gene transfer agent g5 gene reveals bipolar and endemic distribution of Roseobacter clade members in polar coastal seawater. Diversity (Basel)14:392. doi: 10.3390/d14050392
102
ZengY. X.LuoW.LiH. R.YuY. (2021). High diversity of planktonic prokaryotes in Arctic Kongsfjorden seawaters in summer 2015. Polar Biol.44, 195–208. doi: 10.1007/s00300-020-02791-3
103
ZhangX.LinC.ZhouX.LeiK.GuoB.CaoY.et al. (2019). Concentrations, fluxes, and potential sources of nitrogen and phosphorus species in atmospheric wet deposition of the Lake Qinghai watershed, China. Sci. Total Environ.682, 523–531. doi: 10.1016/j.scitotenv.2019.05.224
104
ZhangL.LvJ. (2021). Land-use change from cropland to plantations affects the abundance of nitrogen cycle-related microorganisms and genes in the loess plateau of China. Appl. Soil Ecol.161:103873. doi: 10.1016/j.apsoil.2020.103873
105
ZhangS.QinW.BaiY.ZhangZ.WangJ.GaoH.et al. (2021). Linkages between anammox and denitrifying bacterial communities and nitrogen loss rates in high-elevation rivers. Limnol. Oceanogr.66, 765–778. doi: 10.1002/lno.11641
106
ZhangX.YaoC.ZhangB.TanW.GongJ.WangG.-Y.et al. (2023). Dynamics of benthic nitrate reduction pathways and associated microbial communities responding to the development of seasonal deoxygenation in a coastal mariculture zone. Environ. Sci. Technol.57, 15014–15025. doi: 10.1021/acs.est.3c03994
107
ZhaoY.LiuZ.ZhangB.CaiJ.YaoX.ZhangM.et al. (2023). Inter-bacterial mutualism promoted by public goods in a system characterized by deterministic temperature variation. Nat. Commun.14:5394. doi: 10.1038/s41467-023-41224-7
Summary
Keywords
metagenomics, high-altitude pristine saline lakes, nitrogen cycle, Qinghai-Tibetan plateau, metagenome-assembled genomes
Citation
Zhao Z, Zhao Y, Marotta F, Xamxidin M, Li H, Xu J, Hu B and Wu M (2024) The microbial community structure and nitrogen cycle of high-altitude pristine saline lakes on the Qinghai-Tibetan plateau. Front. Microbiol. 15:1424368. doi: 10.3389/fmicb.2024.1424368
Received
27 April 2024
Accepted
18 June 2024
Published
09 July 2024
Volume
15 - 2024
Edited by
Xianbiao Lin, Ocean University of China, China
Reviewed by
Lianhua Xu, Zhejiang Sci-Tech University, China
Xin Sun, Carnegie Institution for Science (CIS), United States
Updates
Copyright
© 2024 Zhao, Zhao, Marotta, Xamxidin, Li, Xu, Hu and Wu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Min Wu, wumin@zju.edu.cnBaolan Hu, blhu@zju.edu.cn
†These authors have contributed equally to this work and share first authorship
Disclaimer
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.