Next Article in Journal
Bifidobacterium bifidum: A Key Member of the Early Human Gut Microbiota
Next Article in Special Issue
Effects of a Simulated Acute Oil Spillage on Bacterial Communities from Arctic and Antarctic Marine Sediments
Previous Article in Journal
First Proteomic Approach to Identify Cell Death Biomarkers in Wine Yeasts during Sparkling Wine Production
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 7
Issue 11
10.3390/microorganisms7110543
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Nitrogen Content on Archaeal Diversity in an Arctic Lake Region

by
Jinjiang Lv
1,†,
Feng Liu
2,†,
Wenbing Han
1,
Yu Wang
1,
Qian Zhu
1,
Jiaye Zang
2,
Shuang Wang
3,
Botao Zhang
1 and
Nengfei Wang
2,*
1
College of Chemistry and Chemical Engineering, Qingdao University, Qingdao 266071, China
2
Key Lab of Marine Bioactive Substances, First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China
3
Department of Biochemistry and Molecular biology, School of Basic Medicine, Qingdao University, Qingdao 266071, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Microorganisms 2019, 7(11), 543; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms7110543
Submission received: 15 October 2019 / Revised: 4 November 2019 / Accepted: 7 November 2019 / Published: 8 November 2019
(This article belongs to the Special Issue Polar Microbial Ecology: The Role of Microbes in the Functioning of Extremely Cold Ecosystems)
Browse Figures
Review Reports Versions Notes

Abstract

:
The function of Arctic soil ecosystems is crucially important for the global climate, and nitrogen (N) is the major limiting nutrient in these environments. This study assessed the effects of changes in nitrogen content on archaeal community diversity and composition in the Arctic lake area (London Island, Svalbard). A total of 16S rRNA genes were sequenced to investigate archaeal community composition. First, the soil samples and sediment samples were significantly different for the geochemical properties and archaeal community composition. Thaumarchaeota was an abundant phylum in the nine soil samples. Moreover, Euryarchaeota, Woesearchaeota, and Bathyarchaeota were significantly abundant phyla in the three sediment samples. Second, it was found that the surface runoff caused by the thawing of frozen soil and snow changed the geochemical properties of soils. Then, changes in geochemical properties affected the archaeal community composition in the soils. Moreover, a distance-based redundancy analysis revealed that NH4+–N (p < 0.05) and water content were the most significant factors that correlated with the archaeal community composition. Our study suggests that nitrogen content plays an important role in soil archaeal communities. Moreover, archaea play an important role in the carbon and nitrogen cycle in the Arctic lake area.

1. Introduction

Because of the amplification effect of global warming in the Arctic, Arctic temperatures are rising faster than those in low latitudes [1], and the latest research shows that Arctic warming has intensified over the past decade [2]. Climate warming may increase the thawing of permafrost and lead to changes in the distribution and abundance of thaw lakes and ponds [3,4].
In addition to the changes in the landscape, the geochemical properties of the Arctic ecosystem could also change. These properties include the pH, transitions between aerobic and anaerobic conditions, and the content of carbon and nitrogen [5], which will significantly change the biota and the function of ecosystems and may reshape the feedback process in climate change [6]. However, how climate changes translate into alterations of ecosystem dynamics is not well understood [7,8]. Microbial communities driving biogeochemical processes underlying nutrient fluxes may play an important role in the feedback effect of climate warming [9,10]. Therefore, understanding the structure and function of Arctic microbial communities is of great significance in predicting the response of Arctic ecosystems to climate warming [11].
Moreover, archaeal populations are rich in cold and temperate environments [12,13], and archaea not only exist in extremely high-temperature and halophilic environments but also in soils, seawater, fresh water, lake sediments, and other environments [14]. Archaea play an important role in the nitrogen cycle [15]. Ammonia-oxidizing archaea (AOA) mainly belong to Thaumarchaeota, which are ubiquitous and abundant in oligotrophic marine waters, estuaries, and sediments [16,17]. AOA participate in the rate-limiting reactions of nitrification processes and have recently been found to be potentially important in nitrogen cycling in a variety of environments [15]. Methanogenic archaea belong to Euryarchaeota. Some methanogenic archaea contain nitrogen-fixing genes [18]. Moreover, Methanobacterium ivanovi has been demonstrated to be able to fix molecular nitrogen [19,20].
There is abundant evidence indicating that nitrogen is the major limiting nutrient in Arctic soils [21,22]. Over geological timescales, the N cycle is thought to have an effect on the global biogeochemical cycle of carbon and change the content of atmospheric CO2 [23].
The thawing of permafrost and subsequent hydrological changes have also been shown to increase the availability of dissolved organic and inorganic nitrogen [24,25]. Surface runoff is formed by the thawing of frozen soil, and glaciers can erode ammonium nitrogen and nitrate nitrogen in soil, thereby enriching nitrogen during water erosion [26,27]. These changes raised the net primary productivity via more rapid N mineralization and uptake [28], which could affect carbon dynamics by ultimately changing the metabolic efficiency of the microbiome [5].
In this study, we selected the London Island of Konsfjorden as a sample point, as this island has typical slope features. Previous studies mainly focused on the correlation between bacterial communities and environmental factors [29]. The aim of this research is mainly to study the influence of nitrogen content changes caused by running water erosion on archaeal communities.

2. Materials and Methods

2.1. Sampling Site Description and Sample Collection

The study area is located on the London Island of Konsfjorden, which is on the west coast of Spitsbergen in the Svalbard archipelago (Figure 1a). Surface runoff is formed by the thawing of frozen soil and accumulated snow flows down the hillside into the lake during the summer. Samples were collected at four sites in July 2016 during China’s Arctic expedition. Three soil sites and one sediment site (Hill, Up, Down, and Sedi; Figure 1b) were sampled from the surface (5 cm) near each other (about 1 m apart) in triplicate. Fifty grams of samples were collected using a sterile shovel and put directly into TWIRL’EM sterile sampling bags (Labplas Inc., Sainte-Julie, QC, Canada). A total of nine soil samples and three lake sediment samples were obtained. After collection, samples were placed in centrifuge tubes at −20 °C in the Yellow River Station (Ny-Aalesund, Norway) for about 20 days and then taken to the home laboratory in China by plane. During the flights, the samples were conserved in an incubator with frozen ice bags. In the home laboratory, soil samples were frozen at −80 °C until nucleic acid extraction.

2.2. Geochemical Properties of Soils and Lake Sediments

A total of seven soil geochemical properties were measured, including pH, water content, organic carbon (OrC), and four soluble nutrients (NH4+–N, SiO42−–Si, NO3−–N and NO2−–N; Table 1). Soil pH was measured by adding 10 mL of distilled water to four grams of soil for pH measurement using a pH meter (PHS-3C, Shanghai REX Instrument Factory, Shanghai, China). Water content (10 g of each sample) was determined as the gravimetric weight loss after drying the soil/sediment at 105 ℃ until reaching a constant weight. OrC was measured according to the following procedure. After freeze-drying, the soil sample was ground into powder and reacted with 5 mL 10% hydrochloric acid for 12 h, subsequently rinsed with Milli-Q water for complete acid removal, and dried overnight at 50 °C to be analyzed in an element analyzer (EA3000, Euro Vector SpA, Milan, Italy). The soils used to determine the nutrients were also freeze-dried, ground using an Agate mortar, and then water was added at a ratio of 1:10 (g/mL). After shaking once every 4 h for 48 h, a nutrient auto-analyzer (QuAAtro, SEAL, Germany) was used to determine other physical and chemical properties at a relative standard deviation of <5%.

2.3. DNA Extraction, Polymerase Chain Reaction (PCR) Amplification, and Sequencing

2.3.1. DNA Extraction and PCR Amplification

The total genome DNA from the samples was extracted using the CTAB/SDS method. The DNA was extracted from 0.25 g soil from each sample using a Power Soil DNA Isolation Kit (MOBIO Laboratories, San Diego, CA, USA) according to the manufacturer’s instructions. DNA concentration and purity were monitored on 1% agarose gels. According to the concentration, DNA was diluted to 1 ng/μL using sterile water. The V4 and V5 hypervariable regions of the archaea 16S ribosomal RNA gene were amplified by polymerase chain reaction (PCR) using the primers Arch519F (5′-CAGCCGCCGCGGTAA-3′) and Arch915R (5′-GTGCTCCCCCGCCAATTCCT-3′). All PCR reactions were carried out in 30 μL reactions, including 15 μL of Phusion® High-Fidelity PCR Master Mix with a GC Buffer (New England Biolabs, Ipswich, MA, United States), 0.2 μM of forward and reverse primers, and 10 ng template DNA. The PCR amplification cycle was as follows: Initial denaturation at 94 °C for 4 min, followed by 30 cycles of denaturation at 94 °C for 15 s, annealing at 56 °C for 30 s, and elongation at 68 °C for 80 s, with a final extension of 72 °C for 10 min.

2.3.2. PCR Product Quantification, Qualification, and Purification

PCR products were mixed with equal volume of 1X loading buffer (containing SYB green) and loaded onto 2% agarose gel for detection. Samples with a bright main strip between 400–450 bp were chosen for further experiments. PCR products were mixed in equidense ratios. Then, a mixture of PCR products was purified with a Qiagen Gel Extraction Kit (Qiagen, Hilden, Germany).

2.3.3. Library Preparation and Sequencing

Sequencing libraries were generated using the TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, San Diego, CA, USA) following the manufacturer’s recommendations, and index codes were added. The library quality was assessed on a Qubit 2.0 Fluorometer (Thermo Scientific, Wilmington, DE, USA) and Agilent Bioanalyzer 2100 system (Agilent Technologies, Palo Alto, CA, USA). Lastly, the library was sequenced on an Illumina HiSeq2500 platform (San Diego, CA, USA), and 250 bp paired-end reads were generated.

2.3.4. Processing of Sequencing Data

Paired-end reads were assigned to samples based on their unique barcode and truncated by cutting off the barcode and primer sequence, using FLASH (V1.2.7) to merge the paired-end reads to get the Raw Tags. Quality filtering on the raw tags was performed under specific filtering conditions to obtain high-quality clean tags according to QIIME (V1.7.0); then, chimeras were detected by the UCHIME algorithm and against the Gold database. Finally, the chimeric sequences and lllumina adapters were removed. The resulting effective reads of all samples were clustered using the Uparse software (V7.0, http://drive5.com/uparse/), and the sequences were clustered into OTUs (operational taxonomic units) with a 97% identity. Meanwhile, the most frequent sequence for an OTU was selected as the representative sequence of the OTU. The species were annotated and analyzed with a representative sequence of OTUs using QIIME and the SSU rRNA database SILVA (V 128) to obtain taxonomic information and calculate the abundance at each classification level in all the samples. The raw reads were deposited into the NCBI Sequence Read Archive (SRA) database (accession number: SUB6251077).

2.4. Statistical Analyses

The OTU abundance information was normalized using a standard sequence number corresponding to the sample with the least sequences. A subsequent analysis of the alpha diversity and beta diversity were all performed based on the output normalized data. Alpha diversity was applied in analyzing the complexity of species diversity for a sample based on 6 indices: Chao1, Shannon’s index (H′), Pielou, Simpson, ACE, and Good’s coverage. All indices in our samples were performed on QIIME (V 1.7.0) and displayed with R software (V 3.6.0). A one-way analysis of variance (ANOVA) followed by Tukey’s HSD (honest significant difference) test was performed for the geochemical properties and diversity parameters of the samples to determine the level of significance using the Statistical Package for the Social Sciences (SPSS) software (V.17.0). The relationships among the archaeal communities in the 12 samples were analyzed by a hierarchical clustering analysis using the R (V.3.6.0) statistical software. An analysis of similarities (ANOSIM) test was performed to determine whether the four sampling sites had statistically significantly different archaeal communities by using QIIME (V.1.7.0) software. The abundance-based Bray–Curtis similarity coefficient was used to examine the dissimilarity of the different samples. The relevance of environmental factors in explaining the distribution patterns of archaeal communities in different samples was analyzed by a Bray–Curtis distance-based redundancy analysis (db-RDA) using the R (V.3.6.0) statistical software. The input data can be standardized by taking the logarithm and other methods to ensure that continuous variables are sufficiently normally distributed. A Monte Carlo permutation test was also performed to examine the relationship between the 7 geochemical properties and the archaeal community composition in this Arctic lake area. In order to study the phylogenetic relationship of the different OTUs and the different dominant species in the different samples, a cluster analysis (heat map) of the archaeal genera at different sampling sites was performed with the R (V.3.6.0) statistical software. The species accumulation box-plot and Venn diagram were also developed using the by R (V.3.6.0) statistical software A linear discriminant analysis effect size (LEfSe) method was used to identify the significantly different archaeal groups in the different sampling sites.

3. Results

3.1. Geochemical Properties of Soil and Sediment Samples

The experimental results showed that there were some differences in geochemical properties among the four sites. The sediment samples were quite different between the soil samples. Moreover, there was a particular regularity in the change of the geochemical properties along the path of the surface runoff. Firstly, the contents of organic carbon, NO3−–N, and NO2−–N gradually decreased, while the contents of NH4+–N and SiO42−–Si showed an increasing trend. The highest values of NH4+–N and SiO42−–Si were observed from the sediment samples (site Sedi), where the lowest values of NO3−–N, NO2−–N, and organic carbon were also detected. Interestingly, the highest values of NO3−–N and NO2−–N and the lowest values of SiO42−–Si were detected at the Up site, where there are some lichens. Furthermore, the water content of the three soil samples was basically the same (from 8% to 16%). However, the water content of the sediment samples was high at 22%. The pH values of all samples ranged from 7.63 to 8.22, which means that the soils were weakly alkaline. On the other hand, the pH value of the Down site was relatively high. The SiO42−–Si of the Down and Sedi sites was higher than that of the Hill and Up sites, while the organic carbon content was lower.

3.2. Archaeal Diversity and Community Composition

In the 12 samples, a total of 716,728 raw sequences and 669,370 clean reads were retained after a series of quality control measures, and 2976 OTUs (at a 3% evolutionary distance) were identified in this study. Then, we eliminated the bacteria OTUs through the reference database and obtained 344 archaea OTUs (Figure 3). Down.1 contains the fewest OTUs (42), while Sedi.2 contains the most OTUs (276). The Good’s coverage estimator of the OTUs in the samples ranged from 0.997 to 1 (Table 2), indicating that the sequences sufficiently covered most of the archaeal populations in all samples. The Chao and AEC values of the three soil sample sites varied considerably. The lowest value appeared at the Down site, which indicates that the Down site had low species richness. In contrast, the sediment samples had high species richness. The changing trend of Shannon and Simpson is the same; the value of the sediment sample was higher than that of the soil sample, and the three soil sample sites were basically the same, indicating that the sediment had high species diversity. The Pielou of the soil sample ranged from 0.309 to 0.450, which is lower than the Pielou of the sediment samples (0.529–0.705), indicating that the archaeal community evenness of the sediment samples was high. The Shannon (4.116), Simpson (0.851), Chao1 (280.41), and AEC values (284.337) of Sedi.2 were the highest value among the 12 samples. The OTU accumulation box-plot (Figure 2) is saturated with all 12 samples, indicating that the OTUs of the study sites are well represented by these samples.
A Venn diagram (Figure 3) demonstrates that the OTUs differed among the four sites. The number of site-specific OTUs ranged from 5 (Site Down) to 190 (Site Sedi). Only 38 in 344 OTUs were shared by all four sampling sites.
In order to study the effects of nitrogen content change on archaeal communities, the general archaeal community structure was first studied. At different taxonomic ranks, the soil sample sites can be significantly distinguished from the sediment sample sites. For example, at the phylum rank (Figure 4), all OTUs were mapped to 9 phyla; members of Euryarchaeota (44–63%), Woesearchaeota (3–10%), Bathyarchaeota (MCG, Miscellaneous Crenarchaeota Group) (2–8%), and Sm1K20 (1–6%) were significantly abundant in the sediment samples. Thaumarchaeota was distributed in all 12 samples and was the dominate phylum, with a content of 97–99% in the three soil sample sites (Hill, Up, and Down). At the genus rank (Figure 5), the community structures and dominant genera of the Archaea at different sampling sites were quite different. Diapherotrites were unidentified, Parvarchaeota were unidentified, and Woesearchaeota were unidentified; they are close to each other in the phylogenetic tree and are mainly distributed in sediment samples. Indeed, they all belong to the superphyla DPANN [30]. Candidatus_Methanoperedens, Methanobacterium, Methanoregula, Methanosaeta, Methanobacteriaceae unidentified, Methanosarcina, and Methanomassiliicoccus were mainly distributed in sediment samples; except for Candidatus_Methanoperedens, all of them were methanogenic archaea. The abundance of Candidatus_Methanoperedens and Methanobacterium were also relatively high in the sediment samples. Some groups with few studies, such as SM1K20 (unidentified), MCG (Bathyarchaeota) (unidentified), MCG (unidentified), Candidatus lainarchaeum were also mainly distributed in sediments. Candidatus Nitrososphaera were distributed in all sampling sites and were more abundant in the soil samples. Moreover, they were far away from other archaea in the phylogenetic tree.
Besides the community composition, we also focused on archaea with large differences in archaeal abundance. Based on the LEfSe results, 17 taxa showed a linear discriminant analysis (LDA) score greater than 4 (the cutoff for significance test) in the 12 samples (Figure 6). Among all the samples, we found that Thaumarchaeota had the highest LDA (5.55) score, which belonged to the Up site, followed by Euryarchaeota (LDA 5.43). In these 12 samples, there were three phyla with significant abundance differences, including Thaumarchaeota, Euryarchaeota, and Woesearchaeota (LDA score 4.48). The Sedi site had 13 taxa that were more abundant than those of other sites. Methanobacteria and Methanomicrobia were two classes with significant abundance differences among the sediments.

3.3. Correlations between Environmental Variables and Archaeal Community Structure

A distance-based redundancy analysis (db-rda; Figure 7) and Monte Carlo permutation test (Table 3) were performed to examine the relationship between the seven geochemical properties and archaeal community composition in the Arctic area. This analysis also revealed the top 4 genera in the 12 sites (Table S1): Thaumarchaeota (unclassified), Candidatus Nitrososphaera, Methanobacterium, and Candidatus Methanoperedens. A combination of the seven environment factors showed a significant correlation with the archaeal community composition, which explained 69.92% of the archaeal community variation. Among the seven geochemical factors, water content (r2 = 0.6685, p < 0.05) and NH4+–N (r2 = 0.6897, p < 0.05) were the most important factors associated with archaeal community composition in the area. Candidatus Methanoperedens and Methanobacterium were positively correlated with NH4+–N, SiO42−–Si, and Water content, while Thaumarchaeota (unclassified) and Candidatus Nitrososphaera were positively correlated with NO3−–N, NO2−–N, and organic carbon.

4. Discussion

Our research has illustrated how the composition and diversity of archaeal communities in a lake area of the London Island changes with geochemical factors. In the present study, we only examined 344 archaea OTUs (Figure 3). The diversity of archaeal communities is much less than that of bacterial communities [29]. However, we have found that archaea play an important role in the carbon and nitrogen cycle in the Arctic region, by analyzing archaeal communities. Thaumarchaeota was significantly abundant in the three soil-sample sites and also distributed in the sediment samples. Euryarchaeota, Woesearchaeota, and Bathyarchaeota were mainly distributed in the sediment samples (Figure 4).
On London Island, geochemical properties play an important role in the changes of archaeal community diversity in the soil. The combination of the seven geochemical factors showed a significant correlation with archaeal community composition. In this study, NH4+–N (r2 = 0.6897, p < 0.05) and water content (r2 = 0.6685, p < 0.05) showed the highest correlation with archaeal community composition in all samples (Table 3), and other geochemical factors had different effects on different archaeal communities except for the pH factor (Figure 7). The dominant genera of archaea on London Island were Candidatus Nitrososphaera, Methanobacterium, Candidatus Methanoperedens, and Thaumarchaeota (unclassified) (Figure 7). Thaumarchaeota (unclassified) mainly belongs to the SCG (Soil Crenarchaeotic Group), and Candidatus Nitrososphaera belongs to AOA, which can oxidize NH4+ into NO2- and fix carbon dioxide [31,32]. Furthermore, NO2- can be reduced to N2O, which is a greenhouse gas [33]. Oxygen is needed in the process of ammonia-nitrogen oxidation so that the number of AOA is absolutely dominant in surface soil samples and less dominant in sediment samples [34]. All of these observations are consistent with the experimental results. Candidatus Nitrososphaera is negatively correlated with water content and NH4+–N but positively correlated with NO3−–N and NO2−–N (Figure 7). In general, soil ammonia-oxidizing archaea play a more important role in low total ammonia concentration environments where they were below 15 μg NH4+–N (g dw. soil)−1 [13,35]. The total ammonia concentrations of London Island were relatively low, ranging from 0.5 to 4.1 NH4+–N (g dw. soil)−1 (Table 1). Therefore, the AOA on London Island may play an important role in the nitrogen cycle and may be affected by geochemical properties.
Methanobacteria and Methanomicrobia were two significant classes of abundance differences for methanogenic archaea in the sediments (Figure 6). The dominant genus, Methanobacterium, which belongs to Methanobacteria, is a kind of methanogenic archaea and is strictly anaerobic. The formation of an anaerobic environment in flooded soil could explain the positive correlation between methanogenic archaea and water content [36]. However, the current study cannot disclose whether the observed effects of the water regime on methanogenic archaea are predominantly caused by redox effects, substrate availability, or a combination of these two factors. Further studies are needed [36]. In this study, Methanobacterium also had a positive correlation with NH4+–N; however, Methanogenic archaea are sensitive to ammonia nitrogen, and high concentrations of ammonia nitrogen inhibit methane production [37,38]. We suggest that the low concentration of ammonia nitrogen is not enough for inhibition, and the symbiotic relationship of the other microorganisms leads to this result. For example, the latest research results show a potential syntrophic relationship between Woesearchaeota and methanogens [39]. More research is needed in the future.
Candidatus Methanoperedens is a kind of Anaerobic Methanotrophic Archaea (ANME). In this study, Candidatus Methanoperedens belongs to Methanomicrobia, which is classified as ANME-1. ANME can denitrify with CH4 as an electron donor and NO3- as an electron acceptor. This process is named nitrate-/nitrite-dependent anaerobic methane oxidation (N-DAMO) [40,41]. Without this process, there would be an additional 10–60% of CH4 in the atmosphere [42]. In this study, Candidatus Methanoperedens is positively correlated with water content and NH4+–N but negatively correlated with Orc, NO3−–N, and NO2−–N (Figure 7). All of these observations are supported by previous studies. Candidatus Methanoperedens plays an important role in the carbon and nitrogen cycle on London Island.
On London Island, we focused on the influence of nitrogen content changes caused by running water erosion on archaeal communities. Because the sample sites were close to each other, and the soil properties were basically the same as those of the surrounding environment, we believe that the change in the trend of the geochemical factors was due to the erosion of melting water caused by the melting of permafrost and snow. The close proximity of these sample sites to each other could also lead to some limitations. Moreover, the generality of the study is partially limited by the single locale where the samples were taken. Although our use of only 12 sampling sites means that there are not enough data to determine what is happening widely in the Arctic during the process of climate warming, our suite of examined variables and statistical analyses provide an interesting case study of likely ecological changes that are occurring more broadly across portions of the Arctic. In this study, the contents of NO3−–N and NO2−–N gradually decreased, while the contents of NH4+−N and SiO42−–Si showed an increasing trend. This trend was the same as that of some other research results [43]. A portion of NO3−–N comes from the enrichment of snow to the atmosphere [44]. Moreover, NO3−–N can accumulate more significantly in dry soil sites because of the nitrification of AOA. By contrast, NO3−–N is consumed as N-DAMO increases in the sediments. NO2−–N, as an intermediate product of the nitrification and denitrification process, has a limited existence time [45]. Diverse biotopes, metabolic complementation, and the effects of the seven geochemical properties lead us to conclude that NH4+–N significantly affects the soil archaeal communities in this study. NO3−–N and NO2−–N also indicate the distribution and composition of some archaeal communities. The dominant genera of archaea like Candidatus Nitrososphaera, Methanobacterium, and Candidatus Methanoperedens all play an important role in carbon and nitrogen cycles. Through our analysis of archaeal communities, we infer that the main function of water content is to form an anoxic environment; we will undertake further research by determining oxygen content. We will also measure methane or CO2 to study the actual effects of archaea on carbon and nitrogen cycles in the Arctic region. Small ponds are not usually considered in large-scale greenhouse gas (GHG) and global carbon-cycling studies because they cannot be seen with remote sensing tools [46,47]. There are also some other factors that should be analyzed in future studies. Some studies suggest that soil depth and salinity may influence the distribution and composition of archaeal communities [4,48]. The results of this study may also show that, through runoff, terrestrial archaea may be transported into the seawater and sediments along the flow path. Considering the effects of global warming, our study provides reliable data for lakes formed by melting snow and ice in parts of the Arctic region.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/2076-2607/7/11/543/s1.

Author Contributions

Conceptualization, F.L., S.W., J.Z., B.Z. and N.W.; data curation, J.L. and F.L.; formal analysis, J.L.; funding acquisition, N.W.; investigation, N.W.; methodology, J.L., Y.W. and Q.Z.; project administration, F.L. and N.W.; resources, F.L. and N.W.; software, J.L.; supervision, F.L. and S.W. and B.Z.; validation, F.L. and J.Z. and N.W.; visualization, J.L.; writing—original draft, J.L.; writing—review and editing, F.L., J.L. and W.H.

Funding

This research was funded by the National Natural Science Foundation of China (Nos. 41776198), the Basic Scientific Fund for National Public Research Institutes of China (Nos.GY0219Q10), and the Open fund of the Key Lab of Marine Bioactive Substances of The First institute of Oceanography, MNR, (No. MBSMAT-2017-01).

Acknowledgments

Special thanks go to Jiaye Zang, Nengfei Wang, and Botao Zhang for their constructive suggestions on the revision of this article and the other authors for their excellent technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kirtman, B.; Power, S.; Adedoyin, A.; Boer, G.; Bojariu, R.; Camilloni, I.; Doblas-Reyes, F.; Fiore, A.; Kimoto, M.; Meehl, G. Near-term climate change: Projections and predictability. In Climate Change 2013: The Physical Science Basis: Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013. [Google Scholar]
  2. Huang, J.; Zhang, X.; Zhang, Q.; Lin, Y.; Hao, M.; Luo, Y.; Zhao, Z.; Yao, Y.; Chen, X.; Wang, L.; et al. Recently amplified arctic warming has contributed to a continual global warming trend. Nat. Clim. Chang. 2017, 7, 875–879. [Google Scholar] [CrossRef]
  3. Vincent, W.F.; Lemay, M.; Allard, M.; Wolfe, B.B. Adapting to permafrost change: A science framework. Eos Trans. Am. Geophys. Union 2013, 94, 373–375. [Google Scholar] [CrossRef]
  4. Tripathi, B.M.; Kim, M.; Kim, Y.; Byun, E.; Yang, J.W.; Ahn, J.; Lee, Y.K. Variations in bacterial and archaeal communities along depth profiles of alaskan soil cores. Sci. Rep. 2018, 8, 504. [Google Scholar] [CrossRef] [PubMed]
  5. Kwon, M.J.; Jung, J.Y.; Tripathi, B.M.; Gockede, M.; Lee, Y.K.; Kim, M. Dynamics of microbial communities and CO2 and CH4 fluxes in the tundra ecosystems of the changing arctic. J. Microbiol. 2019, 57, 325–336. [Google Scholar] [CrossRef]
  6. Schuur, E.A.; Mack, M.C. Ecological response to permafrost thaw and consequences for local and global ecosystem services. Annu. Rev. Ecol. Evol. Syst. 2018, 49, 279–301. [Google Scholar] [CrossRef]
  7. Doney, S.C.; Ruckelshaus, M.; Duffy, J.E.; Barry, J.P.; Chan, F.; English, C.A.; Galindo, H.M.; Grebmeier, J.M.; Hollowed, A.B.; Knowlton, N. Climate change impacts on marine ecosystems. Ann. Rev. Mar. Sci. 2012, 4, 11–37. [Google Scholar] [CrossRef]
  8. Bellard, C.; Bertelsmeier, C.; Leadley, P.; Thuiller, W.; Courchamp, F. Impacts of climate change on the future of biodiversity. Ecol. Lett. 2012, 15, 365–377. [Google Scholar] [CrossRef] [Green Version]
  9. Singh, B.K.; Bardgett, R.D.; Smith, P.; Reay, D.S. Microorganisms and climate change: Terrestrial feedbacks and mitigation options. Nat. Rev. Microbiol. 2010, 8, 779–790. [Google Scholar] [CrossRef]
  10. Chonggang, X.; Chao, L.; Stan, W.; Cathy, W.; Nathan, M.D. Importance of feedback loops between soil inorganic nitrogen and microbial communities in the heterotrophic soil respiration response to global warming. Nat. Rev. Microbiol. 2011, 9, 222. [Google Scholar]
  11. Blaud, A.; Lerch, T.Z.; Phoenix, G.K.; Osborn, A.M. Arctic soil microbial diversity in a changing world. Res. Microbiol. 2015, 166, 796–813. [Google Scholar] [CrossRef]
  12. Høj, L.; Olsen, R.A.; Torsvik, V.L. Effects of temperature on the diversity and community structure of known methanogenic groups and other archaea in high arctic peat. ISME J. 2008, 2, 37–48. [Google Scholar] [CrossRef] [PubMed]
  13. Schleper, C.; Nicol, G.W. Ammonia-oxidising archaea—Physiology, ecology and evolution. Adv. Microb. Physiol. 2010, 57, 1–41. [Google Scholar] [PubMed]
  14. Francis, C.A.; Roberts, K.J.; Beman, J.M.; Santoro, A.E.; Oakley, B.B. Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc. Natl. Acad. Sci. USA 2005, 102, 14683–14688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Hongyue, D.; Xiaoxia, Z.; Jin, S.; Tiegang, L.; Zhinan, Z.; Guanpin, Y. Diversity and spatial distribution of sediment ammonia-oxidizing crenarchaeota in response to estuarine and environmental gradients in the changjiang estuary and east china sea. Microbiology 2008, 154, 2084–2095. [Google Scholar]
  16. Vetriani, C.; Jannasch, H.W.; Macgregor, B.J.; Stahl, D.A.; Reysenbach, A.L. Population structure and phylogenetic characterization of marine benthic archaea in deep-sea sediments. Appl. Environ. Microbiol. 1999, 65, 4375–4384. [Google Scholar]
  17. Brochier-Armanet, C.; Boussau, B.; Gribaldo, S.; Forterre, P. Mesophilic crenarchaeota: Proposal for a third archaeal phylum, the thaumarchaeota. Nat. Rev. Microbiol. 2008, 6, 245. [Google Scholar] [CrossRef]
  18. Kessler, P.; Blank, C.; Leigh, J.A. The nif gene operon of the methanogenic archaeon methanococcus maripaludis. J. Bacteriol. 1998, 180, 1504. [Google Scholar]
  19. Magingo, F.S.S.; Stumm, C.K. Nitrogen fixation by methanobacterium formicicum. FEMS Microbiol. Lett. 1991, 81, 273–277. [Google Scholar] [CrossRef]
  20. Leigh, J.A. Nitrogen fixation in methanogens: The archaeal perspective. Curr. Issues Mol. Biol. 2000, 2, 125–131. [Google Scholar]
  21. Nordin, A.; Schmidt, I.K.; Shaver, G.R. Nitrogen uptake by arctic soil microbes and plants in relation to soil nitrogen supply. Ecology 2004, 85, 955–962. [Google Scholar] [CrossRef]
  22. Shaver, G.R.; Chapin, F.S., III. Response to fertilization by various plant growth forms in an alaskan tundra: Nutrient accumulation and growth. Ecology 1980, 61, 662–675. [Google Scholar] [CrossRef]
  23. Falkowski, P.G. Evolution of the nitrogen cycle and its influence on the biological sequestration of co2 in the ocean. Nature 1997, 387, 272–275. [Google Scholar] [CrossRef]
  24. Chen, L.; Liu, L.; Mao, C.; Qin, S.; Wang, J.; Liu, F.; Blagodatsky, S.; Yang, G.; Zhang, Q.; Zhang, D. Nitrogen availability regulates topsoil carbon dynamics after permafrost thaw by altering microbial metabolic efficiency. Nat. Commun. 2018, 9, 3951. [Google Scholar] [CrossRef] [PubMed]
  25. Finger, R.A.; Turetsky, M.R.; Kielland, K.; Ruess, R.W.; Mack, M.C.; Euskirchen, E.S. Effects of permafrost thaw on nitrogen availability and plant-soil interactions in a boreal alaskan lowland. J. Ecol. 2016, 104, 1542–1554. [Google Scholar] [CrossRef]
  26. Palis, R.G.; Okwach, G.; Rose, C.W.; Saffigna, P.G. Soil erosion processes and nutrient loss. I. The interpretation of enrichment ratio and nitrogen loss in runoff sediment. Aust. J. Soil Res. 1990, 28, 623–639. [Google Scholar] [CrossRef]
  27. Douglas, C.L.; King, K.A.; Zuzel, J.F. Nitrogen and phosphorus in surface runoff and sediment from a wheat-pea rotation in northeastern oregon. J. Environ. Qual. 1998, 27, 1170–1177. [Google Scholar] [CrossRef]
  28. Grant, R.F. Ecosystem co2 and ch4 exchange in a mixed tundra and a fen within a hydrologically diverse arctic landscape: 2. Modeled impacts of climate change. J. Geophys. Res. Biogeosci. 2015, 120, 1388–1406. [Google Scholar] [CrossRef]
  29. Wang, N.F.; Zhang, T.; Yang, X.; Wang, S.; Yu, Y.; Dong, L.L.; Guo, Y.D.; Ma, Y.X.; Zang, J.Y. Diversity and composition of bacterial community in soils and lake sediments from an arctic lake area. Front. Microbiol. 2016, 7, 1170. [Google Scholar] [CrossRef]
  30. Dombrowski, N.; Lee, J.-H.; Williams, T.A.; Offre, P.; Spang, A. Genomic diversity, lifestyles and evolutionary origins of dpann archaea. FEMS Microbiol. Lett. 2019, 366, fnz008. [Google Scholar] [CrossRef]
  31. Roland, H.; Lebedeva, E.V.; Eva, S.; Kilian, S.; Andreas, R.; Holger, D.; Michael, W. A moderately thermophilic ammonia-oxidizing crenarchaeote from a hot spring. Proc. Natl. Acad. Sci. USA 2008, 105, 2134–2139. [Google Scholar] [Green Version]
  32. Li-Mei, Z.; Offre, P.R.; Ji-Zheng, H.; Verhamme, D.T.; Nicol, G.W.; Prosser, J.I. Autotrophic ammonia oxidation by soil thaumarchaea. Proc. Natl. Acad. Sci. USA 2010, 107, 17240–17245. [Google Scholar] [Green Version]
  33. Miranda, K.M.; Ridnour, L.; Esprey, M.; Citrin, D.; Thomas, D.; Mancardi, D.; Donzelli, S.; Wink, D.A.; Katori, T.; Tocchetti, C.G. Comparison of the chemical biology of NO and HNO: An inorganic perspective. Prog. Inorg. Chem. 2005, 54, 349–384. [Google Scholar]
  34. Michael, P.; Christa, S.; Michael, W. The thaumarchaeota: An emerging view of their phylogeny and ecophysiology. Curr. Opin. Microbiol. 2011, 14, 300–306. [Google Scholar]
  35. Jennifer, P.; Dumont, M.G.; Ralf, C. Ammonia oxidation coupled to co2 fixation by archaea and bacteria in an agricultural soil. Proc. Natl. Acad. Sci. USA 2011, 108, 4170–4175. [Google Scholar]
  36. Lone, H.; Marte, R.; Lars Egil, H.; Olsen, R.A.; Torsvik, V.L. Effects of water regime on archaeal community composition in arctic soils. Environ. Microbiol. 2010, 8, 984–996. [Google Scholar]
  37. Kayhanian, M. Performance of a high--solids anaerobic digestion process under various ammonia concentrations. J. Chem. Technol. Biotechnol. 2010, 59, 349–352. [Google Scholar] [CrossRef]
  38. Jin, R.; Huang, G.; Ma, C.; Yang, G. Ammonia inhibition of anaerobic digestion. Ind. Water Treat. 2010, 30, 9–12. [Google Scholar]
  39. Liu, X.; Li, M.; Castelle, C.J.; Probst, A.J.; Zhou, Z.; Pan, J.; Liu, Y.; Banfield, J.F.; Gu, J.D. Insights into the ecology, evolution, and metabolism of the widespread woesearchaeotal lineages. Microbiome 2018, 6, 102. [Google Scholar] [CrossRef]
  40. Cui, M.; Ma, A.; Qi, H.; Zhuang, X.; Zhuang, G. Anaerobic oxidation of methane: An “active” microbial process. Microbiologyopen 2015, 4, 1–11. [Google Scholar] [CrossRef]
  41. Vaksmaa, A.; Guerrero-Cruz, S.; Alen, T.A.V.; Cremers, G.; Ettwig, K.F.; Lüke, C.; Jetten, M.S.M. Enrichment of anaerobic nitrate-dependent methanotrophic ‘candidatus methanoperedens nitroreducens’ archaea from an italian paddy field soil. Appl. Microbiol. Biotechnol. 2017, 101, 7075–7084. [Google Scholar] [CrossRef]
  42. Conrad, R. The global methane cycle: Recent advances in understanding the microbial processes involved. Environ. Microbiol. Rep. 2010, 1, 285–292. [Google Scholar] [CrossRef] [PubMed]
  43. Molénat, J.; Durand, P.; Gascuel-Odoux, C.; Davy, P.; Gruau, G. Mechanisms of nitrate transfer from soil to stream in an agricultural watershed of french brittany. Water Air Soil Pollut. 2002, 133, 161–183. [Google Scholar] [CrossRef]
  44. Sharma, B.M.; Tayal, S.; Chakraborty, P.; Bharat, G.K. Chemical characterization of meltwater from east rathong glacier vis-à-vis Western Himalayan Glaciers. In Dynamics of Climate Change and Water Resources of Northwestern Himalaya; Springer: Cham, Switzerland, 2015. [Google Scholar]
  45. Cao, L.K.; Chen, G.J.; Lu, Y.T. Nitrogen leaching in vegetable fields in the suburbs of shanghai. Pedosphere 2005, 15, 641–645. [Google Scholar]
  46. Tranvik, L.J.; Downing, J.A.; Cotner, J.B.; Loiselle, S.A.; Striegl, R.G.; Ballatore, T.J.; Dillon, P.; Finlay, K.; Fortino, K.; Knoll, L.B. Lakes and reservoirs as regulators of carbon cycling and climate. Limnol. Oceanogr. 2009, 54, 2298–2314. [Google Scholar] [CrossRef] [Green Version]
  47. Muster, S.; Langer, M.; Heim, B.; Westermann, S.; Boike, J. Subpixel heterogeneity of ice-wedge polygonal tundra: A multi-scale analysis of land cover and evapotranspiration in the lena river delta, siberia. Tellus Ser. B Chem. Phys. Meteorol. 2012, 64, 1–14. [Google Scholar] [CrossRef]
  48. Edmonds, J.W.; Weston, N.B.; Joye, S.B.; Mou, X.; Moran, M.A. Microbial community response to seawater amendment in low-salinity tidal sediments. Microb. Ecol. 2009, 58, 558–568. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 07 00543 g001 550
Figure 1. (a) Location; (b) the sampling sites from an Arctic lake area in the present study.
Figure 1. (a) Location; (b) the sampling sites from an Arctic lake area in the present study.
Microorganisms 07 00543 g001
Microorganisms 07 00543 g002 550
Figure 2. The species accumulation box-plot of the 12 samples.
Figure 2. The species accumulation box-plot of the 12 samples.
Microorganisms 07 00543 g002
Microorganisms 07 00543 g003 550
Figure 3. A Venn diagram displaying the degree of overlap of the archaea OTUs (at a 3% evolutionary distance) among the four sampling sites.
Figure 3. A Venn diagram displaying the degree of overlap of the archaea OTUs (at a 3% evolutionary distance) among the four sampling sites.
Microorganisms 07 00543 g003
Microorganisms 07 00543 g004 550
Figure 4. The top 10 abundance of the different phyla in the 12 soil samples of the present study.
Figure 4. The top 10 abundance of the different phyla in the 12 soil samples of the present study.
Microorganisms 07 00543 g004
Microorganisms 07 00543 g005 550
Figure 5. A heat map of the archaeal genus at different sampling sites in the Arctic lake area. A phylogenetic tree is on the left. The changes of the color gradient indicate the richness of the archaeal genus.
Figure 5. A heat map of the archaeal genus at different sampling sites in the Arctic lake area. A phylogenetic tree is on the left. The changes of the color gradient indicate the richness of the archaeal genus.
Microorganisms 07 00543 g005
Microorganisms 07 00543 g006 550
Figure 6. The taxa showing different relative abundances among the sites using the linear discriminant analysis (LDA).
Figure 6. The taxa showing different relative abundances among the sites using the linear discriminant analysis (LDA).
Microorganisms 07 00543 g006
Microorganisms 07 00543 g007 550
Figure 7. A distance-based redundancy analysis to show the correlations between the bacterial communities and environmental factors of the 12 samples from the four sampling sites.
Figure 7. A distance-based redundancy analysis to show the correlations between the bacterial communities and environmental factors of the 12 samples from the four sampling sites.
Microorganisms 07 00543 g007
Table
Table 1. Geochemical properties of the 12 samples investigated in the present study.
Table 1. Geochemical properties of the 12 samples investigated in the present study.
SiteSampleWater ContentpHOrganic Carbon (OrC) %NH4+–N (μg/g)SiO42−–Si (μg/g)NO3−–N (μg/g)NO2−–N (μg/g)
HillHill.10.147.631.180.942.250.060.04
Hill.20.137.971.001.892.320.300.22
Hill.30.148.031.060.501.850.620.04
average0.14 ± 0.017.88 ± 0.221.08 ± 0.091.11 ± 0.712.14 ± 0.250.32 ± 0.280.10 ± 0.10
UpUp.10.167.621.532.851.430.620.23
Up.20.117.681.161.851.750.520.15
Up.30.107.850.120.752.190.030.04
average0.12 ± 0.037.72 ± 120.94 ± 0.731.82 ± 1.051.79 ± 0.380.39 ± 0.320.14 ± 0.10
DownDown.10.118.130.251.154.360.200.03
Down.20.118.220.210.933.110.200.03
Down.30.088.110.161.394.500.440.06
average0.10 ± 0.028.15 ± 0.060.21 ± 0.051.16 ± 0.233.99 ± 0.770.28 ± 0.140.04 ± 0.02
SedimentSedi.10.187.910.113.273.110.050.01
Sedi.20.237.840.514.084.290.010.01
Sedi.30.197.900.313.443.060.010.01
average0.20 ± 0.037.88 ± 0.040.31 ± 0.23.60 ± 0.423.49 ± 0.700.02 ± 0.020.01 ± 0
Table
Table 2. Summary data for the Miseq sequencing data from the 12 samples in the present study.
Table 2. Summary data for the Miseq sequencing data from the 12 samples in the present study.
Sample NameRaw TagEffectiveOTUsShannon
(H’)		Pielou
(J’)		Simpson
(1-D)		chao1ACEGood’s Coverage
Hill.12921723991752.4910.426 0.7738390.2540.999
Hill.22155817420742.6290.450 0.79665.88271.3680.999
Hill.32560818117872.2480.385 0.707122.667123.3250.998
Up.12685423756562.5640.439 0.7876674.7210.999
Up.22650322461602.3920.410 0.77554.66756.9630.999
Up.32363520012692.2130.379 0.75662.54566.6660.999
Down.12857525156421.8030.309 0.63236.2536.2961
Down.22794422992502.1870.374 0.71540.90943.9271
Down.32660723313652.1140.362 0.64353.07757.3510.999
Sedi.127538205141963.090.529 0.789191.091195.7250.998
Sedi.231650235422764.1160.705 0.851280.41284.3370.997
Sedi.319396147102063.4140.585 0.757190.634194.9810.999
Table
Table 3. A Monte Carlo permutation test of the relationship between environmental factors and archaeal community composition.
Table 3. A Monte Carlo permutation test of the relationship between environmental factors and archaeal community composition.
RDA1RDA2r2Pr (>r)
Water content−0.9634780.2677880.66850.008**
pH0.123615−0.9923300.02090.908
OrC0.8977150.4405760.33450.149
NH4+–N−0.9958610.0908890.68970.008**
SiO42−–Si−0.985292−0.1708790.16670.431
NO3−–N0.9898570.1420670.38760.108
NO2−–N0.999554−0.0298760.20800.362
**correlation is significant at the 0.01 level.

Share and Cite

MDPI and ACS Style

Lv, J.; Liu, F.; Han, W.; Wang, Y.; Zhu, Q.; Zang, J.; Wang, S.; Zhang, B.; Wang, N. The Effect of Nitrogen Content on Archaeal Diversity in an Arctic Lake Region. Microorganisms 2019, 7, 543. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms7110543

AMA Style

Lv J, Liu F, Han W, Wang Y, Zhu Q, Zang J, Wang S, Zhang B, Wang N. The Effect of Nitrogen Content on Archaeal Diversity in an Arctic Lake Region. Microorganisms. 2019; 7(11):543. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms7110543

Chicago/Turabian Style

Lv, Jinjiang, Feng Liu, Wenbing Han, Yu Wang, Qian Zhu, Jiaye Zang, Shuang Wang, Botao Zhang, and Nengfei Wang. 2019. "The Effect of Nitrogen Content on Archaeal Diversity in an Arctic Lake Region" Microorganisms 7, no. 11: 543. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms7110543

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop