Next Article in Journal
The Mechanism by Which Umbrella-Shaped Ratchet Trichomes on the Elaeagnus angustifolia Leaf Surface Collect Water and Reflect Light
Next Article in Special Issue
Harnessing Diesel-Degrading Potential of an Antarctic Microalga from Greenwich Island and Its Physiological Adaptation
Previous Article in Journal
The Spatial Precision of Contextual Feedback Signals in Human V1
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 12
Issue 7
10.3390/biology12071023
biology-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

Study of Archaeal Diversity in the Arctic Meltwater Lake Region

by
Yiling Qin
1,†,
Nengfei Wang
2,*,†,
Li Zheng
1,
Qinxin Li
3,
Long Wang
4,
Xiaoyu Xu
2 and
Xiaofei Yin
1
1
First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China
2
School of Chemistry and Chemical Engineering, Linyi University, Linyi 276000, China
3
College of Chemistry and Chemical Engineering, Qingdao University, Qingdao 266071, China
4
Department of Bioengineering, College of Marine Sciences and Biological Engineering, Qingdao University of Science & Technology, Qingdao 266042, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Biology 2023, 12(7), 1023; https://doi.org/10.3390/biology12071023
Submission received: 9 May 2023 / Revised: 4 July 2023 / Accepted: 16 July 2023 / Published: 20 July 2023
(This article belongs to the Special Issue Microbial Diversity and Environmental Adaptability of Microorganisms in Extreme Environments)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

The Arctic is experiencing a significant temperature increase under the global warming trend. As a result, the Arctic permafrost is thawing, glacial meltwater is gathering, and the depressions are gradually forming lakes, affecting the composition and material cycles of the terrestrial ecosystems. Two meltwater lakes with different landscapes in the Ny-Ålesund region of the Arctic were taken as study objects. The structure of the archaeal community and the influence of soil physiochemical factors on the archaeal community were investigated. The differences in the archaeal community structure between the intertidal and subtidal zones of the two lakes were compared, and the reasons for the differences were analyzed and discussed. A redundancy analysis identified NH4+, SiO32−, MC, NO3−, and NO2− as key soil physicochemical factors that have significantly influenced the structure of the archaeal community. The hub archaea in the archaeal community were identified by weighted gene co-expression network analysis (WGCNA). The use of WGCNA in this study provides new ideas for future research on the microbial community structure. In the context of global warming, this study contributes to research on archaeal communities in Arctic meltwater lakes in response to climate change.

Abstract

Two typical lakes formed from meltwater in the Ny-Ålesund area were taken as the study subjects in 2018. To investigate the archaeal community compositions of the two lakes, 16S rRNA genes from soil samples from the intertidal and subtidal zones of the two lakes were sequenced with high throughput. At the phylum level, the intertidal zone was dominated by Crenarchaeota and the subtidal zone was dominated by Halobacter; at the genus level, the intertidal zone was dominated by Nitrososphaeraceae_unclassified and Candidatus_Nitrocosmicus, while the subtidal zone was dominated by Methanoregula. The soil physicochemical factors pH, moisture content (MC), total organic carbon (TOC), total organic nitrogen (TON), nitrite nitrogen (NO2−-N), and nitrate nitrogen (NO3−-N) were significantly different in the intertidal and subtidal zones of the lake. By redundancy analysis, the results indicated that NH4+-N, SiO32−-Si, MC, NO3−-N, and NO2−-N have had highly significant effects on the archaeal diversity and distribution. A weighted gene co-expression network analysis (WGCNA) was used to search for hub archaea associated with physicochemical factors. The results suggested that these physicochemical factors play important roles in the diversity and structure of the archaeal community at different sites by altering the abundance of certain hub archaea. In addition, Woesearchaeales was found to be the hub archaea genus at every site.

1. Introduction

Under the global warming trend, the Arctic region is one of the most obvious regions affected by the temperature rise in the last hundred years [1]. The surface temperature in the Arctic is changing more than twice as fast as the global warming average, which means that there is an “Arctic amplification effect” [2]. As Arctic temperatures continue to rise, depressions created by thawing permafrost and then the melting of glaciers gradually form lakes, which leads to changes in terrestrial ecosystems and material cycles [3,4]. Lakes are an important part of the Arctic environment, and climate change has had a significant impact on Arctic lakes, with repeated freezing and thawing directly affecting their abundance and distribution [5]. The Svalbard lake survey showed that meltwater lakes are one of the major lake types, accounting for 30.4% of the lake area [6].
Changes in the abundance and distribution of lakes in the Arctic cause changes in the soil water content, nutrients, microbial diversity, and ecosystem function [5]. In Arctic soils, the microbial diversity is rich and dominates as decomposers in the terrestrial environment of the polar regions [7]. Changes in some key elements of the soil can cause changes in the original biogeochemical cycle, which can significantly alter the structure of biological communities as well as the ecosystem function [8]. Archaea are important microorganisms that play important roles in the material cycling processes of geochemistry, such as carbon cycling and nitrogen cycling processes [9]. The community structure and distribution of archaea are influenced by various physicochemical factors, such as the temperature [10], pH [11], and nitrogen [12], phosphorus [13], oxygen [14], soil moisture, and organic matter contents [15].
Microbial studies of Arctic lakes have mainly focused on bacteria with few studies conducted on archaea. The research on archaea in the Arctic Svalbard region has focused on marine sediments [16], glaciers [17], tundra [18], and Arctic fauna [19]. In the few studies of archaea in Svalbard lakes, some researchers have investigated the vertical structure of lake archaeal communities [20], but there have been no relevant studies on the structural differences between intertidal and subtidal archaeal communities. It is clear that the archaea in the lakes in the Ny-Ålesund area have been poorly and incompletely studied. Not only the structure and differences in the archaeal community composition in the lake sediments but also the correlation between the archaeal communities and relevant environmental factors have not been investigated. In addition, there are complex interrelationships among microorganisms, such as competition, symbiosis, and parasitism [21], and the specific forms are competitive substrates [22], direct electron transfer between species [23], and complementary metabolites [24].
In this study, intertidal and subtidal soil samples were collected from two lakes with different characteristics in the Ny-Ålesund area to reveal the structure and differences in soil archaeal communities and to explore the role and influence of environmental physicochemical factors on the structure and diversity of the archaeal communities. The archaeal diversity composition of the soil samples was analyzed by high-throughput 16S rRNA sequencing. At the same time, the physicochemical properties of the soil samples were determined by conventional instrumental analysis to further study their effects on the structure of soil archaeal communities and to identify the most important influencing factors. WGCNA was then used to find the hub that influences the archaeal community’s composition.

2. Materials and Methods

2.1. Study Site Description and Sample Collection

The study site is the Ny-Ålesund area, which is located in the Svalbard archipelago of the Arctic. Svalbard is located at the intersection of the northern North Atlantic Ocean and the Arctic Ocean basin, and therefore, the climate there differs significantly from the rest of the Arctic [25,26]. Ny-Ålesund is located on the western coast of Spitsbergen, which is the largest island in Svalbard. Due to its proximity to, and influence from, the warm Western Spitsbergen Current, the climate of Ny-Ålesund is generally warmer and wetter than that of other land masses at the same latitude [27]. Samples were collected in July 2018 from 12 sites in the Ny-Ålesund region of the Arctic (Table 1) with three parallel samples set at each site, giving a total of 36 soil individual samples. Figure 1 shows the locations of the two lakes on the map. The intertidal zone refers to the area that used to be above water, but which became submerged as meltwater increased due to warmer temperatures. The subtidal zone refers to the area that is submerged throughout the year. Approximately 50 g of soil was dug up with a sterile spoon into sterile bags at a depth of approximately 5 cm. At the end of sampling, the samples were stored at −20 °C immediately and promptly transferred to an ultra-low temperature refrigerator (−80 °C).

2.2. Physicochemical Properties of Soils

Nine physicochemical factors of the soil samples were measured, among which were the pH, moisture content (MC), total organic carbon (TOC), total organic nitrogen (TON), and the concentrations of five soluble nutrients: nitrite nitrogen (NO2−-N), nitrate nitrogen (NO3−-N), ammonium nitrogen (NH4+-N), silicate silicon (SiO32−-Si), and phosphate phosphorus (PO43−-P). The pH of the soil was measured with a pH meter (PHS-3C, Shanghai REX Instrument Factory, Shanghai, China) by taking 4 g of the soil sample and adding 10 mL of distilled water. The MC was measured by weighing 10 g of the soil samples, drying at 105 °C to a constant weight, and calculating the weight loss. Soil samples were freeze-dried immediately after removal from the −80 °C refrigerator for the determination of other physicochemical factors. The freeze-dried samples were decalcified with hydrochloric acid, rinsed with deionized water, and dried before being analyzed for TOC and TON using an element analyzer (EA30000, Euro Vector SpA, Milan, Italy) [28]. The contents of the five soluble nutrients were determined using freeze-dried samples treated with deionized water at a ratio of 1:10 (g∙mL−1) and analyzed by an automated nutrient analyzer (QuAAtro, SEAL, Germany) with a relative standard deviation of less than 0.5%. One-way ANOVA and Duncan’s test were performed on the obtained physicochemical factors’ data using Origin 2021.

2.3. DNA Extraction, PCR Amplification and Sequencing

Total DNA extraction from the soil was performed according to the Mo Bio Power SOIL DNA isolation kit manufacturer’s instructions, and the DNA obtained was tested for DNA purity and integrity in agarose gel. DNA at concentrations of 10–30 ng/μL−1 was used for PCR amplification with Arch519F (5′-CAGCCGCCGCGGTAA-3′) and Arch915R (5′-GTGCTCCCCCCCGCCAATTCCT-3′) [29] as primers. The V4-V5 region of the 16S rRNA gene was amplified according to the method described by He et al. [15], and PCR products of 400–450 bp were recovered with a gel extraction kit and sent to a commercial sequencing company for sequencing.

2.4. High-Throughput Sequencing and Statistical Analysis

The 16S rRNA gene V4-V5 region was sequenced with high throughput on the Illumina Miseq platform [30], and low-quality sequences were removed by quality control and cascading in order to obtain clean data. Original read segments were submitted to the NCBI Sequence Read Archive database (accession number: PRJNA956435). We performed quality control on the raw data, which included data splicing, filtering, and chimera removal processing. The DADA2 method in the QIIME 2.0 system was used for noise reduction to perform dereplication or the equivalent to form clustering at 100% similarity [31]. The ASVs (Amplicon Sequence Variants) were each de-duplicated sequence (corresponding to OTU representative sequences) obtained after noise reduction. Species annotations were made for each ASV by the classify-sklearn algorithm in QIIME 2.0 using Silva 138.1 [32]. Statistical analysis of the diversity indices of ASVs was performed using the QIIME 2.0 system. The QIIME2 system was used to calculate the observed_otus, shannon, simpson, chao1, goods_coverage, dominance, and pielou_e indices. A significantly different species analysis between groups was performed using LEFSE in the QIIME2 system. A redundancy analysis (RDA) was conducted using R software (version 4.0.5) to explore the correlation between soil physicochemical properties and the distribution of the archaeal community. The RDA analysis started with a multiple regression of the matrix of the species composition and the matrix of the environmental variables, and then the fitted value matrix was obtained. PCA was performed on the fitted value matrix to obtain the canonical eigenvector matrix, and then the image was plotted. The weighted gene co-expression network analysis (WGCNA) software package was used to differentiate the modules [21]. Then, the heatmap of modules and physicochemical factors and the network map of the core archaea were drawn. The results of the WGCNA were combined with the results of the RDA to analyze the influence of soil physicochemical factors on the structure of the archaeal community.

3. Results

3.1. Physicochemical Properties of Soil Samples

Nine physicochemical properties, including the pH, MC, TOC, TON, NO2−-N, NO3−-N, NH4+-N, SiO32−-Si, and PO43−-P, were measured in 36 soil samples from 12 sites. As can be seen from Table 2, except for NH4+-N, SiO32−-Si, and PO43−-P, other physicochemical factors were significantly different in the intertidal and subtidal zone of the lake. The contents of NO2−-N and NO3−-N in NHS were significantly higher than that of NHX, the content of NH4+-N was slightly higher, while PO43−-P was the opposite. The pH of XHS was significantly higher than that of XHX, while the TOC was reversed. The contents of NO2−-N and NO3−-N in XHX2 and XHS3 were significantly higher than those in other places. Among them, XHX2 had the highest content of SiO32−-Si, and XHX3 had the highest contents of NH4+-N and PO43−-P.

3.2. Diversity and Structure Analysis of Archaeal Community

A total of 36 samples from two lakes in Ny-Ålesund, Arctic, were sequenced with high throughput, and 1,960,571 archaeal sequences were obtained, with an average of 54,460 archaeal sequences per sample. Clustering by the DADA2 method yielded 1906 ASVs for subsequent analysis. After the alpha diversity analysis, the values of Good’s coverage in all sites were greater than 0.9999, indicating that the coverage of sequencing was high enough, and the vast majority of sequence regions were obtained by sequencing.
As can be seen from Figure 2a, except for the unclassified ones, all samples of ASVs were clustered into seven phyla. Among them, the highest relative abundance was found for Halobacterota, followed by Crenarchaeota, and these two phyla accounted for a much higher proportion of all sites than other phyla. However, in general, the distributions of the two phyla in intertidal and subtidal soils differed significantly, with Crenarchaeota predominating in the intertidal zones and Halobacterota predominating in the subtidal zones. Additionally, the relative abundance of Euryarchaeota in the intertidal and subtidal samples also differed significantly, as it was more abundant in subtidal zones than in intertidal zones.
In Figure 2b, it can be seen that the relative abundance of each genus at each sampling site differed significantly. In general, the distributions of archaea genera in two lakes were significantly different, and the dominant genera in the lake with birds were less significant than in the small lake. For example, Nitrososphaeraceae_unclassified, Methanoregula, and Candidatus_Nitrocosmicus had high relative abundances in the small lake. There were also differences between the intertidal and subtidal zones with high relative abundances of Nitrososphaeraceae_unclassified and Candidatus_Nitrocosmicus in the intertidal zones and Methanoregula and Methanosaeta in the subtidal zones.

3.3. Correlation between Soil Physicochemical Factors and Archaeal Community Structure

RDA was used to explore the relationships among nine soil physicochemical factors and archaeal communities. As shown in Figure 3, the first two axes explained 40.42% of the total variation in the structure of the archaeal community. The results showed that NH4+, SiO32−, MC, NO3−, and NO2− had highly significant effects on the archaeal diversity and distribution. Among all sampling sites, NH4+, SiO32−, and MC had the greatest effects on XHX2, while NO3− and NO2− had the greatest effects onXHS3. Then, the Monte Carlo permutation test was applied, and the results were as follows (Table 3): NH4+ (r2 = 0.4399, Pr = 0.001), SiO32− (r2 = 0.4911, Pr = 0.001), MC (r2 = 0.4264, Pr = 0.001), NO3− (r2 = 0.4254, Pr = 0.001), and NO2− (r2 = 0.5296, Pr = 0.001). The results further verified the significance of the effects of these five soil physicochemical factors on the archaeal community. In addition, TOC and pH had significant effects on the diversity and distribution of archaeal communities.
We not only studied the physicochemical factors but also applied LEFSE to focus on archaea with significant large differences in the relative abundance at each site. The results are shown in Figure 4a. There were significant differences in archaea at nine of the twelve studied sites. For genera with significant differences, NHX2 had the most differential archaea, followed by NHX3, and then XHS3 and NHX3. The LEFSE results for the two lakes are shown in Figure 4b with five distinctly different archaeal genera in NH (lake with birds) and XH (small lake) identified.

3.4. The Weighted Gene Co-Expression Network Analysis

The weighted gene co-expression network analysis (WGCNA) method was used to analyze the relationship between soil physicochemical factors and all ASVs annotated clearly as archaea, and the results are presented as a heatmap of the modules and physicochemical factors (Figure 5). The turquoise module showed a highly significant positive correlation with PO43− (r = 0.85, p < 0.01), NH4+ (r = 0.68, p < 0.01), and TOC (r = 0.57, p < 0.01). The gray module also showed a significant positive correlation with NH4+ (r = 0.61, p < 0.01) and significant negative correlations with NO3− (r = −0.48, p < 0.01), NO2− (r = −0.46, p < 0.01), and TON (r = −0.51, p < 0.01). The results of the turquoise module were imported into Cystoscope and plotted as a network diagram (Figure 6). As can be seen from the figure, there were 19 hub archaea groups, including three genera, Woesearchaeales, Candidatus_Nitrososphaera, and Methanobacterium, which were mainly Woesearchaeales. The network diagram of the gray module can be seen in Figure 7. The results show that the hub archaea differed among sites, but the hub genera all included Woesearchaeales. The hub genera in NHX2 were Woesearchaeales and Methanosaeta, the hub genera in XHS were Woesearchaeales and Candidatus_Nitrososphaera, and the hub genera in XHX3.1 and XHX3.2 were Woesearchaeales and Methanoregula, while Woesearchaeales, Candidatus_Nitrososphaera and Methanobacterium were identified as the hub genera in XHX3.3.

4. Discussion

This study revealed the composition and diversity of archaeal communities in a typical Arctic lake area and analyzed the differences between intertidal and subtidal areas. The results indicate that soil physicochemical factors influence the archaeal community composition. In this study, 1906 archaeal ASVs were identified with Good’s coverage values exceeding 99.99% at each sampling site. Although not as diverse as the bacterial communities in similar meltwater areas [33], richly diverse archaeal communities were found when compared to previous archaeal studies in the Svalbard region [34]. The 16S rRNA high-throughput sequencing data showed that the phyla Halobacterota and Crenarchaeota were absolutely dominant, and there were significant differences in the intertidal and subtidal zones. Halobacterota was reported as the major archaeal phylum in extreme environments, such as mud volcanoes [35], Antarctic salt cones [36], and deserts [37], as it can mostly tolerate saline environments but is less tolerant to oxygen and desiccation [38]. Oxygen levels and aridity are higher in the intertidal zone than in the subtidal zone, so the abundance of Halobacterota archaea is higher in the subtidal zone than in the intertidal zone. It has been shown that Halobacterota has sulfate allosteric reduction genes involved in the sulfur cycle, and some genes participate in the carbon cycle, such as those of the order Methanomicrobiales [39,40]. Crenarchaeota has been reported as the major archaeal phylum in both Antarctic and Arctic archaeal studies [41,42], and it was also the major archaeal phylum in previous Svalbard archaeal studies [34]. This is a diverse and widespread phylum that is mainly characterized by acidophilic, thermophilic, and anaerobic properties [43]. It also contains ammonia-oxidizing archaea, which are important players in the natural nitrogen cycle [44]. In the genus classification, the dominant genera were Nitrososphaeraceae_unclassified, Methanoregula, and Candidatus_Nitrocosmicus. The archaeal community structures in the intertidal and subtidal zones were significantly different, and those in the two lakes were also significantly different.
In the Arctic, soil physicochemical factors influence the diversity of archaeal communities and show significant correlations with the archaeal community composition. In this study, there were differences in the physicochemical properties of the soils that had different effects on the archaeal communities at different sampling sites. According to the results of the RDA, NH4+, SiO32−, MC, NO3−, and NO2− were highly correlated with the archaeal community composition of all samples. In general, these five physicochemical factors had greater influences on the small lake (XH) than on the lake with birds (NH). However, the different physicochemical factors affected the intertidal and subtidal zones to different extents. MC was significantly higher in the subtidal than in the intertidal zone, while the opposite was true for NO3− and NO2−. pH and TOC differed significantly only between XHS and XHX, while NH4+, SiO32−, and PO43− were high at some sites but not at most. These differences may result in differences in the composition of archaeal communities, while at the same time, the archaeal communities may also counteract the physicochemical properties of the soil [45]. In addition, global warming has affected the abundance and distribution of Arctic meltwater lakes which, in turn, has led to changes in the soil water content, nutrients, and microbial diversity [5]. In this study, the water content had a significant effect on the archaeal community. Therefore, the archaeal communities of the lakes might be affected by changes in the meltwater volume due to climate change. There was a study that showed that Arctic peat archaeal communities were dominated by nonmethanogenic archaea at lower temperatures, and the diversity and relative abundance of methanogenic archaea increased at higher temperatures [10]. Research has shown that glaciers have unique archaeal community physicochemical factors [46], so the glacial meltwater input may affect the local archaeal communities of lakes. These changes may show the response of the archaeal community to environmental changes.
The dominant genera in both lakes were mainly Nitrososphaeraceae (unclassified) and Candidatus_Nitrocosmicus in the intertidal zones and Methanoregula and Methanosaeta in the subtidal zones. The difference in dominant genera between the two regions may be related to the significantly higher NO3− and NO2− concentrations in the intertidal zones than in the subtidal zones. Nitrososphaeraceae belongs to the phylum Crenarchaeota and produces energy by oxidizing ammonia to nitrite under aerobic conditions [47]. Nitrososphaeraceae species are involved in the soil nitrogen cycle and play an important role [48], which explains the higher abundance of Nitrososphaeraceae in the intertidal than in the subtidal zones. Candidatus_Nitrocosmicus also belongs to the phylum Crenarchaeota and is a genus in the family Nitrososphaeraceae. It has a strong aerobic ammonia oxidation biofunction and can oxidize high concentrations of ammonia [49,50]. As observed in both Figure 2b and Figure 4a, the relative abundance of the genus Candidatus_Nitrocosmicus at XHS3 was significantly higher than that at the other sites, corresponding to the apparently high NO3− and NO2− concentrations at this site. Methanoregula and Methanosaeta are both members of the phylum Halobacterota, but in different classes. Methanoregula belongs to the class Methanomicrobia and is a strictly anaerobic, thermophilic organism that produces methane from H2, CO2 or formate [51,52]. Global distribution studies of sequences from this genus indicate that this lineage survives in a variety of environments [53], and some species tend to grow in nitrogen-rich environments [54]. This partly explains the significantly higher relative abundance of Methanoregula at XHX2 and XHX3, which have significantly higher NH4+ contents, than at other sites. Methanosaeta belongs to the class Methanosarcinia and is also an anaerobic methanogenic archaeon. It uses CO2 and methyl acetate as substrates for methane production and is more competitive than the genus Methanosarcina, which also consumes methyl acetate, regardless of the concentration of acetic acid [22,55]. It has been shown that 60% of the methane produced by organisms and released into the atmosphere comes from methyl acetate [56], so it is thought that Methanosaeta is the main methane producer on the Earth [55]. Methanosarcina spp. are acetic acid trophic methanogenic bacteria that compete with Methanosaeta spp. for acetic acid while also producing methane and growing on formate, methanol, and hydrogen as substrates [57]. Methanosaeta and Methanosarcinia frequently combine [58], which may be one reason why both genera simultaneously have a significantly higher relative abundance of NHX2 than other sites. The relative abundance of Halobacterales was significantly higher in XHS2. It is a salt-loving archaea, generally considered to have evolved from salt-loving methanogenic archaea [59], and is commonly found in sea salt [60].
The comparison of LEFSE from the two lakes showed that NH (lake with birds) had five archaeal genera whose relative abundances were significantly different from those in XH (small lake). They are mainly concentrated in classifications under the phylum Halobacterota, such as the genera Methanosarcina, Methanosphaerula, and Methanospirillum. The latter two are both strictly anaerobic archaea that use hydrogen, carbon dioxide, or formate to produce methane [61,62]. Combined with Figure 2a, it can be seen that the relative abundance of the phylum Euryarchaeota of NH is significantly higher than that of XH, and the relative abundance of the genus Methanobacterium in this phylum varies significantly. The genus Methanobacterium uses H2/CO2 and formate as substrates to synthesize methane and also to syntrophically degrade formate [63]. SiO32−, NO3−, NO2−, TOC, and MC were significantly different in the two lakes, and the RDA results show that TOC was correlated, and several other physicochemical factors were significantly correlated. However, the archaea with different abundances in the two lakes were concentrated in the methanogenic archaea type, so whether and how these physicochemical factors caused the differences in archaea remain to be investigated in a follow-up study. The two lakes have very different landscapes, with NH being the lake with a lot of bird habitats. Previous studies have shown that the archaeal diversity index increases with an increasing bird density [64], which is consistent with the results of this study. The results of the alpha diversity analysis in this study (Table S2) also show that the archaea diversity in NH is higher than that in XH. The physiological activities of birds affect the soil environment [65] and thus the structure of archaeal communities, but there is a lack of studies in this area.
WGCNA is an R package, and its basic process is to first cluster highly correlated genes and form a module, summarize these clusters with hub genes, correlate the modules with each other and with out-of-sample features and calculate the correlation, and finally find the hub genes [66]. It is widely used in the biomedical field to analyze hub genes for diseases, such as aortic dissection [67], colorectal cancer [68], and endometrial carcinoma [69]. In this study, we used WGCNA to explore the correlations between highly correlated archaeal modules and physicochemical factors to identify significantly influential factors and hub archaea. The results showed that the turquoise and gray modules correlated strongly with a number of physicochemical factors, and the correlations of NH4+, SiO32−, NO3−, and NO2− could be corroborated with the RDA results. There were differences in the hub archaea genera compositions in different environmental sites, but Woesearchaeales was the hub genus for each site. For example, the hub archaea were AR5 in NHS2, Methanosaeta in NHX2, and Methanospirillum in NHX3. The hub archaea of XHS were Candidatus_Nitrososphaera and Methanoregula of XHX, while the hub archaea of XHX3.3 had the same Candidatus_Nitrososphaera as XHS and Methanobacterium. Candidatus_Nitrososphaera is an ammonia-oxidizing archaea that uses NH3 and O2 for energy and has been shown to have the potential to use urea as a source of NH3 [70]. Methanobacterium is a hydrogenotrophic methanogenic archaea that uses CO2 and can also fix nitrogen [71,72]. A previous showed a significantly positive correlation between Methanobacterium and NH4+ [34]. Therefore, the differences in hub archaea in XHX3.3 may be caused by the significantly high NH4+ concentration at this site. Woesearchaeales are anaerobic archaea from the phylum Nanoarchaeota and are widely distributed in the environment, such as in Antarctic soils [41], trench sediments [73], sewage [74], and cows [75]. It has been shown that Woesearchaeales has positive correlations with Methanoregula and Methanosarcina, which have symbiotic or parasitic relationships with methanogens [24,74]. This may be an important reason why Woesearchaeales was identified as a hub archaeon.
In bacterial culture, a ‘sandwich agar plate method’ has been used, in which ‘helper’ bacteria are used as a medium sandwich to co-culture other marine bacteria [76]. Subsequent experiments have also shown the potential of the sandwich agar plate method for isolating uncultured bacteria [77]. Generally speaking, archaea are difficult to cultivate. Using bacteria as a reference, we would be able to co-cultivate more archaea with ‘helper’ archaea. However, the method of finding ‘helper’ bacteria through experimental testing is not applicable to archaea due to their noncultivability. Rapidly developing genomic technologies provide tremendous support for archaeal research [78]. However, many hypotheses about the evolution, physiology and diversity of archaea proposed by culture-independent methods have yet to be confirmed by culture experiments [79]. Improvements in traditional culture methods are therefore urgently needed. Therefore, if the hub archaea can be found through WGCNA, it may be possible to improve the traditional medium by virtue of its association with other archaea, so as to achieve the purpose of enriching culturable archaea. WGCNA was well used in our research to identify hub archaea.

5. Conclusions

It was found that significant differences in the archaea community structure exist not only between the intertidal zone and the subtidal zone in typical Arctic meltwater lakes, but also in lakes with different characteristics. According to the RDA analysis, NH4+, SiO32−, MC, NO3−, and NO2− were significantly and strongly correlated with the composition of the archaea community, as TOC and pH also had certain effects. Furthermore, WGCNA was used to explore the correlations among archaea modules and physicochemical factors. Then, the hub archaea were obtained. The sites with different environmental characteristics had different hub archaeal structures, but Woesearchaeales was the hub archaeon at each site. This study provides new ideas for future microbial diversity research and archaea culturability.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology12071023/s1, Table S1: Number of sample sequences; Table S2: Alpha index.

Author Contributions

Conceptualization, N.W., X.Y. and L.Z.; methodology, N.W., X.Y. and L.Z.; software, Y.Q., Q.L. and L.W.; validation, X.Y.; formal analysis, N.W. and X.Y.; resources, N.W.; data curation, Y.Q., Q.L., L.W. and X.X.; writing—original draft preparation, Y.Q. and X.X.; writing—review and editing, N.W., X.Y. and L.Z.; visualization, Y.Q., Q.L. and L.W.; supervision, N.W., X.Y. and L.Z.; project administration, N.W.; funding acquisition, N.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Foundation of Shandong Province National (No. ZR2020KC036), the National Natural Science Foundation of China (No. 41776198), and the Key R&D Program of China (2018YFC1406700).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. 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]
  2. Ding, M.; Wang, S.; Sun, W. Decadal Climate Change in Ny-Ålesund, Svalbard, A Representative Area of the Arctic. Condens. Matter 2018, 3, 12. [Google Scholar] [CrossRef] [Green Version]
  3. Christensen, T.R. It’s a gas. Nat. Geosci. 2016, 9, 647–648. [Google Scholar] [CrossRef]
  4. Kramshoj, M.; Albers, C.N.; Svendsen, S.H.; Bjorkman, M.P.; Lindwall, F.; Bjork, R.G.; Rinnan, R. Volatile emissions from thawing permafrost soils are influenced by meltwater drainage conditions. Glob. Chang. Biol. 2019, 25, 1704–1716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Schuur, E.A.G.; 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]
  6. Urbański, J.A. Monitoring and classification of high Arctic lakes in the Svalbard Islands using remote sensing. Int. J. Appl. Earth Obs. Geoinf. 2022, 112, 102911. [Google Scholar] [CrossRef]
  7. Malard, L.A.; Pearce, D.A. Microbial diversity and biogeography in Arctic soils. Environ. Microbiol. Rep. 2018, 10, 611–625. [Google Scholar] [CrossRef]
  8. Sjursen, H.; Michelsen, A.; Holmstrup, M. Effects of freeze–thaw cycles on microarthropods and nutrient availability in a sub-Arctic soil. Appl. Soil Ecol. 2005, 28, 79–93. [Google Scholar] [CrossRef]
  9. Dang, H.; Zhang, X.; Sun, J.; Li, T.; Zhang, Z.; Yang, G. 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] [CrossRef] [Green Version]
  10. 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. 2007, 2, 37–48. [Google Scholar] [CrossRef]
  11. Gubry-Rangin, C.; Novotnik, B.; Mandič-Mulec, I.; Nicol, G.W.; Prosser, J.I. Temperature responses of soil ammonia-oxidising archaea depend on pH. Soil Biol. Biochem. 2017, 106, 61–68. [Google Scholar] [CrossRef]
  12. Peay, K.G.; von Sperber, C.; Cardarelli, E.; Toju, H.; Francis, C.A.; Chadwick, O.A.; Vitousek, P.M. Convergence and contrast in the community structure of Bacteria, Fungi and Archaea along a tropical elevation-climate gradient. FEMS Microbiol. Ecol. 2017, 93, fix045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Paula, F.S.; Chin, J.P.; Schnurer, A.; Muller, B.; Manesiotis, P.; Waters, N.; Macintosh, K.A.; Quinn, J.P.; Connolly, J.; Abram, F.; et al. The potential for polyphosphate metabolism in Archaea and anaerobic polyphosphate formation in Methanosarcina mazei. Sci. Rep. 2019, 9, 17101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Qin, W.; Meinhardt, K.A.; Moffett, J.W.; Devol, A.H.; Virginia Armbrust, E.; Ingalls, A.E.; Stahl, D.A. Influence of oxygen availability on the activities of ammonia-oxidizing archaea. Environ. Microbiol. Rep. 2017, 9, 250–256. [Google Scholar] [CrossRef]
  15. He, S.; Tan, J.; Hu, W.; Mo, C. Diversity of Archaea and Its Correlation with Environmental Factors in the Ebinur Lake Wetland. Curr. Microbiol. 2019, 76, 1417–1424. [Google Scholar] [CrossRef] [PubMed]
  16. Treude, T.; Krause, S.; Steinle, L.; Burwicz, E.; Hamdan, L.J.; Niemann, H.; Feseker, T.; Liebetrau, V.; Krastel, S.; Berndt, C. Biogeochemical Consequences of Nonvertical Methane Transport in Sediment Offshore Northwestern Svalbard. J. Geophys. Res. Biogeosciences 2020, 125, e2019JG005371. [Google Scholar] [CrossRef]
  17. Zarsky, J.D.; Stibal, M.; Hodson, A.; Sattler, B.; Schostag, M.; Hansen, L.H.; Jacobsen, C.S.; Psenner, R. Large cryoconite aggregates on a Svalbard glacier support a diverse microbial community including ammonia-oxidizing archaea. Environ. Res. Lett. 2013, 8, 035044. [Google Scholar] [CrossRef]
  18. Lee, S.H.; Jang, I.; Chae, N.; Choi, T.; Kang, H. Organic layer serves as a hotspot of microbial activity and abundance in Arctic tundra soils. Microb. Ecol. 2013, 65, 405–414. [Google Scholar] [CrossRef]
  19. Sundset, M.A.; Edwards, J.E.; Cheng, Y.F.; Senosiain, R.S.; Fraile, M.N.; Northwood, K.S.; Praesteng, K.E.; Glad, T.; Mathiesen, S.D.; Wright, A.D. Rumen microbial diversity in Svalbard reindeer, with particular emphasis on methanogenic archaea. FEMS Microbiol. Ecol. 2009, 70, 553–562. [Google Scholar] [CrossRef] [Green Version]
  20. Pouliot, J.; Galand, P.E.; Lovejoy, C.; Vincent, W.F. Vertical structure of archaeal communities and the distribution of ammonia monooxygenase A gene variants in two meromictic High Arctic lakes. Environ. Microbiol. 2009, 11, 687–699. [Google Scholar] [CrossRef]
  21. Zuo, J.; Zu, M.; Liu, L.; Song, X.; Yuan, Y. Composition and diversity of bacterial communities in the rhizosphere of the Chinese medicinal herb Dendrobium. BMC Plant Biol. 2021, 21, 127. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, S.; Cheng, H.; Liu, J.; Hazen, T.C.; Huang, V.; He, Q. Unexpected competitiveness of Methanosaeta populations at elevated acetate concentrations in methanogenic treatment of animal wastewater. Appl. Microbiol. Biotechnol. 2017, 101, 1729–1738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Zheng, S.; Li, M.; Liu, Y.; Liu, F. Desulfovibrio feeding Methanobacterium with electrons in conductive methanogenic aggregates from coastal zones. Water Res. 2021, 202, 117490. [Google Scholar] [CrossRef] [PubMed]
  24. 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]
  25. van der Bilt, W.G.M.; Bakke, J.; Vasskog, K.; D’Andrea, W.J.; Bradley, R.S.; Ólafsdóttir, S. Reconstruction of glacier variability from lake sediments reveals dynamic Holocene climate in Svalbard. Quat. Sci. Rev. 2015. [Google Scholar] [CrossRef] [Green Version]
  26. de Wet, G.A.; Balascio, N.L.; D’Andrea, W.J.; Bakke, J.; Bradley, R.S.; Perren, B. Holocene glacier activity reconstructed from proglacial lake Gjøavatnet on Amsterdamøya, NW Svalbard. Quat. Sci. Rev. 2018, 183, 188–203. [Google Scholar] [CrossRef]
  27. Pedersen, Å.Ø.; Convey, P.; Newsham, K.K.; Mosbacher, J.B.; Fuglei, E.; Ravolainen, V.; Hansen, B.B.; Jensen, T.C.; Augusti, A.; Biersma, E.M.; et al. Five decades of terrestrial and freshwater research at Ny-Ålesund, Svalbard. Polar Res. 2022, 41. [Google Scholar] [CrossRef]
  28. Hu, L.; Shi, X.; Yu, Z.; Lin, T.; Wang, H.; Ma, D.; Guo, Z.; Yang, Z. Distribution of sedimentary organic matter in estuarine–inner shelf regions of the East China Sea: Implications for hydrodynamic forces and anthropogenic impact. Mar. Chem. 2012, 142–144, 29–40. [Google Scholar] [CrossRef]
  29. Wei, S.; Cui, H.; Zhang, Y.; Su, X.; Dong, H.; Chen, F.; Zhu, Y. Comparative evaluation of three archaeal primer pairs for exploring archaeal communities in deep-sea sediments and permafrost soils. Extremophiles 2019, 23, 747–757. [Google Scholar] [CrossRef]
  30. Li, Y.; Zhu, X.; Zhang, W.; Zhu, D.; Zhou, X.; Zhang, L. Archaeal communities in the deep-sea sediments of the South China Sea revealed by Illumina high-throughput sequencing. Ann. Microbiol. 2019, 69, 839–848. [Google Scholar] [CrossRef]
  31. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef] [PubMed]
  33. 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] [Green Version]
  34. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Megyes, M.; Móga, J.; Strat, D.; Borsodi, A.K. Bacterial and Archaeal Taxonomic Diversity of Mud Volcanoes (Beciu, Romania) via Metagenomic Approach. Geomicrobiol. J. 2021, 38, 532–539. [Google Scholar] [CrossRef]
  36. Azzaro, M.; Papale, M.; Rizzo, C.; Forte, E.; Lenaz, D.; Guglielmin, M.; Lo Giudice, A. Antarctic Salt-Cones: An Oasis of Microbial Life? The Example of Boulder Clay Glacier (Northern Victoria Land). Microorganisms 2022, 10, 1753. [Google Scholar] [CrossRef]
  37. Fuentes, B.; Choque, A.; Gomez, F.; Alarcon, J.; Castro-Nallar, E.; Arenas, F.; Contreras, D.; Morchen, R.; Amelung, W.; Knief, C.; et al. Influence of Physical-Chemical Soil Parameters on Microbiota Composition and Diversity in a Deep Hyperarid Core of the Atacama Desert. Front. Microbiol. 2021, 12, 794743. [Google Scholar] [CrossRef]
  38. Conrad, R. Methane Production in Soil Environments-Anaerobic Biogeochemistry and Microbial Life between Flooding and Desiccation. Microorganisms 2020, 8, 881. [Google Scholar] [CrossRef]
  39. Mo, S.; Li, J.; Li, B.; Yu, R.; Nie, S.; Zhang, Z.; Liao, J.; Jiang, Q.; Yan, B.; Jiang, C. Impacts of Desulfobacterales and Chromatiales on sulfate reduction in the subtropical mangrove ecosystem as revealed by SMDB analysis. bioRxiv 2020. [CrossRef]
  40. Anderson, I.; Ulrich, L.E.; Lupa, B.; Susanti, D.; Porat, I.; Hooper, S.D.; Lykidis, A.; Sieprawska-Lupa, M.; Dharmarajan, L.; Goltsman, E.; et al. Genomic characterization of methanomicrobiales reveals three classes of methanogens. PLoS ONE 2009, 4, e5797. [Google Scholar] [CrossRef] [Green Version]
  41. Li, Q.; Wang, N.; Han, W.; Zhang, B.; Zang, J.; Qin, Y.; Wang, L.; Liu, J.; Zhang, T. Soil Geochemical Properties Influencing the Diversity of Bacteria and Archaea in Soils of the Kitezh Lake Area, Antarctica. Biology 2022, 11, 1855. [Google Scholar] [CrossRef] [PubMed]
  42. Bano, N.; Ruffin, S.; Ransom, B.; Hollibaugh, J.T. Phylogenetic Composition of Arctic Ocean Archaeal Assemblages and Comparison with Antarctic Assemblages. Appl. Environ. Microbiol. 2004, 70, 781–789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Santos, A.; Bruna, P.; Martinez-Urtaza, J.; Solis, F.; Valenzuela, B.; Zamorano, P.; Barrientos, L. Two Archaeal Metagenome-Assembled Genomes from El Tatio Provide New Insights into the Crenarchaeota Phylum. Genes 2021, 12, 391. [Google Scholar] [CrossRef] [PubMed]
  44. Nicol, G.W.; Schleper, C. Ammonia-oxidising Crenarchaeota: Important players in the nitrogen cycle? Trends Microbiol. 2006, 14, 207–212. [Google Scholar] [CrossRef] [PubMed]
  45. Lin, X.; Green, S.; Tfaily, M.M.; Prakash, O.; Konstantinidis, K.T.; Corbett, J.E.; Chanton, J.P.; Cooper, W.T.; Kostka, J.E. Microbial community structure and activity linked to contrasting biogeochemical gradients in bog and fen environments of the Glacial Lake Agassiz Peatland. Appl. Environ. Microbiol. 2012, 78, 7023–7031. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Anesio, A.M.; Lutz, S.; Chrismas, N.A.M.; Benning, L.G. The microbiome of glaciers and ice sheets. NPJ Biofilms Microbiomes 2017, 3, 10. [Google Scholar] [CrossRef] [Green Version]
  47. Stieglmeier, M.; Klingl, A.; Alves, R.J.E.; Rittmann, S.K.R.; Melcher, M.; Leisch, N.; Schleper, C. Nitrososphaera viennensis gen. nov., sp. nov., an aerobic and mesophilic, ammonia-oxidizing archaeon from soil and a member of the archaeal phylum Thaumarchaeota. Int. J. Syst. Evol. Microbiol. 2014, 64, 2738–2752. [Google Scholar] [CrossRef]
  48. Qiu, L.; Zhang, Q.; Zhu, H.; Reich, P.B.; Banerjee, S.; van der Heijden, M.G.A.; Sadowsky, M.J.; Ishii, S.; Jia, X.; Shao, M.; et al. Erosion reduces soil microbial diversity, network complexity and multifunctionality. ISME J. 2021, 15, 2474–2489. [Google Scholar] [CrossRef]
  49. Cheng, X.; Xiang, X.; Yun, Y.; Wang, W.; Wang, H.; Bodelier, P.L.E. Archaea and their interactions with bacteria in a karst ecosystem. Front. Microbiol. 2023, 14, 1068595. [Google Scholar] [CrossRef]
  50. Wu, F.; Zhang, Y.; He, D.; Gu, J.-D.; Guo, Q.; Liu, X.; Duan, Y.; Zhao, J.; Wang, W.; Feng, H. Community structures of bacteria and archaea associated with the biodeterioration of sandstone sculptures at the Beishiku Temple. Int. Biodeterior. Biodegrad. 2021, 164, 105290. [Google Scholar] [CrossRef]
  51. Brauer, S.L.; Cadillo-Quiroz, H.; Ward, R.J.; Yavitt, J.B.; Zinder, S.H. Methanoregula boonei gen. nov., sp. nov., an acidiphilic methanogen isolated from an acidic peat bog. Int. J. Syst. Evol. Microbiol. 2011, 61, 45–52. [Google Scholar] [CrossRef] [Green Version]
  52. Yashiro, Y.; Sakai, S.; Ehara, M.; Miyazaki, M.; Yamaguchi, T.; Imachi, H. Methanoregula formicica sp. nov., a methane-producing archaeon isolated from methanogenic sludge. Int. J. Syst. Evol. Microbiol. 2011, 61, 53–59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Yang, S.; Liebner, S.; Winkel, M.; Alawi, M.; Horn, F.; Dörfer, C.; Ollivier, J.; He, J.-S.; Jin, H.; Kühn, P.; et al. In-depth analysis of core methanogenic communities from high elevation permafrost-affected wetlands. Soil. Biol. Biochem. 2017, 111, 66–77. [Google Scholar] [CrossRef]
  54. Zinder, S.; Bräuer, S. Methanoregula. In Bergey’s Manual of Systematics of Archaea and Bacteria; Wiley: Hoboken, NJ, USA, 2016; pp. 1–8. [Google Scholar] [CrossRef]
  55. Smith, K.S.; Ingram-Smith, C. Methanosaeta, the forgotten methanogen? Trends Microbiol. 2007, 15, 150–155. [Google Scholar] [CrossRef]
  56. Bridgham, S.D.; Cadillo-Quiroz, H.; Keller, J.K.; Zhuang, Q. Methane emissions from wetlands: Biogeochemical, microbial, and modeling perspectives from local to global scales. Glob. Chang. Biol. 2013, 19, 1325–1346. [Google Scholar] [CrossRef] [PubMed]
  57. Janssen, P.H. Selective enrichment and purification of cultures of Methanosaeta spp. J. Microbiol. Methods 2003, 52, 239–244. [Google Scholar] [CrossRef]
  58. Falz, K.Z.; Holliger, C.; Grosskopf, R.; Liesack, W.; Hahn, D. Vertical Distribution of Methanogens in the Anoxic Sediment of Rotsee (Switzerland). Appl. Environ. Microbiol. 1999, 65, 2402–2408. [Google Scholar] [CrossRef]
  59. Andam, C.P.; Harlow, T.J.; Papke, R.T.; Gogarten, J.P. Ancient origin of the divergent forms of leucyl-tRNA synthetases in the Halobacteriales. BMC Evol. Biol. 2012, 12, 85. [Google Scholar] [CrossRef] [Green Version]
  60. Satari, L.; Guillen, A.; Latorre-Perez, A.; Porcar, M. Beyond Archaea: The Table Salt Bacteriome. Front. Microbiol. 2021, 12, 714110. [Google Scholar] [CrossRef]
  61. Iino, T.; Mori, K.; Suzuki, K.I. Methanospirillum lacunae sp. nov., a methane-producing archaeon isolated from a puddly soil, and emended descriptions of the genus Methanospirillum and Methanospirillum hungatei. Int. J. Syst. Evol. Microbiol. 2010, 60, 2563–2566. [Google Scholar] [CrossRef] [Green Version]
  62. Cadillo-Quiroz, H.; Yavitt, J.B.; Zinder, S.H. Methanosphaerula palustris gen. nov., sp. nov., a hydrogenotrophic methanogen isolated from a minerotrophic fen peatland. Int. J. Syst. Evol. Microbiol. 2009, 59, 928–935. [Google Scholar] [CrossRef] [PubMed]
  63. Xu, J.; Jia, Z.; Lin, X.; Feng, Y. DNA-based stable isotope probing identifies formate-metabolizing methanogenic archaea in paddy soil. Microbiol. Res. 2017, 202, 36–42. [Google Scholar] [CrossRef] [PubMed]
  64. Domínguez, M.T.; Gutiérrez, E.; González-Domínguez, B.; Román, M.; Ávila, J.M.; Ramo, C.; Gonzalez, J.M.; García, L.V. Impacts of protected colonial birds on soil microbial communities: When protection leads to degradation. Soil. Biol. Biochem. 2017, 105, 59–70. [Google Scholar] [CrossRef] [Green Version]
  65. Blonska, E.; Lasota, J.; Jankowiak, R.; Michalcewicz, J.; Wojas, T.; Zbyryt, A.; Ciach, M. Biological and physicochemical properties of the nests of White Stork Ciconia ciconia reveal soil entirely formed, modified and maintained by birds. Sci. Total Environ. 2021, 763, 143020. [Google Scholar] [CrossRef]
  66. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef] [Green Version]
  67. Zhang, H.; Chen, T.; Zhao, W.; Zhang, Y.; Tang, L.; Xu, B.; Chen, X. Identification of crucial genes in aortic dissection by WGCNA. Atherosclerosis 2021, 331, e244–e246. [Google Scholar] [CrossRef]
  68. Di, Z.; Zhou, S.; Xu, G.; Ren, L.; Li, C.; Ding, Z.; Huang, K.; Liang, L.; Yuan, Y. Single-cell and WGCNA uncover a prognostic model and potential oncogenes in colorectal cancer. Biol. Proced. Online 2022, 24, 13. [Google Scholar] [CrossRef]
  69. Liu, J.; Wu, Z.; Sun, R.; Nie, S.; Meng, H.; Zhong, Y.; Nie, X.; Cheng, W. Using mRNAsi to identify prognostic-related genes in endometrial carcinoma based on WGCNA. Life Sci. 2020, 258, 118231. [Google Scholar] [CrossRef]
  70. Zhalnina, K.V.; Dias, R.; Leonard, M.T.; Quadros, P.; Camargo, F.; Drew, J.C.; Farmerie, W.G.; Daroub, S.H.; Triplett, E.W. Genome Sequence of Candidatus Nitrososphaera evergladensis from Group I.1b Enriched from Everglades Soil Reveals Novel Genomic Features of the Ammonia-Oxidizing Archaea. PLoS ONE 2014, 9, e101648. [Google Scholar] [CrossRef] [Green Version]
  71. Bhatnagar, L.; Jain, M.K.; Aubert, J.P.; Zeikus, J.G. Comparison of Assimilatory Organic Nitrogen, Sulfur, and Carbon Sources for Growth of Methanobacterium Species. Appl. Environ. Microbiol. 1984, 48, 785–790. [Google Scholar] [CrossRef] [Green Version]
  72. Magingo, F.S.S.; Stumm, C.K. Nitrogen fixation by Methanobacterium formicicum. FEMS Microbiol. Lett. 1991, 81, 273–277. [Google Scholar] [CrossRef]
  73. Yang, N.; Tian, C.; Lv, Y.; Hou, J.; Yang, Z.; Xiao, X.; Zhang, Y. Novel primers for 16S rRNA gene-based archaeal and bacterial community analysis in oceanic trench sediments. Appl. Microbiol. Biotechnol. 2022, 106, 2795–2809. [Google Scholar] [CrossRef] [PubMed]
  74. Kuroda, K.; Narihiro, T.; Shinshima, F.; Yoshida, M.; Yamaguchi, H.; Kurashita, H.; Nakahara, N.; Nobu, M.K.; Noguchi, T.Q.P.; Yamauchi, M.; et al. High-rate cotreatment of purified terephthalate and dimethyl terephthalate manufacturing wastewater by a mesophilic upflow anaerobic sludge blanket reactor and the microbial ecology relevant to aromatic compound degradation. Water Res. 2022, 219, 118581. [Google Scholar] [CrossRef] [PubMed]
  75. Walker, M.; Holton, M.; Lourenco, J.; Fontes, P. 299 Differences in Microbial Community Composition between Uterine Horns Ipsilateral and Contralateral to the Corpus Luteum in Beef Cows on day 15 of the Estrous Cycle. J. Anim. Sci. 2022, 100, 142–143. [Google Scholar] [CrossRef]
  76. Liu, Z. Preliminary Application of Sandwich Plate Method in Isolation of Bacteria from Marine Sediments and Multiphase Classification of Two New Strain. Master’s Thesis, Shandong University, Ji’nan, China, 2020. [Google Scholar]
  77. Zhang, J. The Establishment of the Sandwich Plate Method and Its Application in the Mining of Bacterial Resources in Offshore Sediments. PhD Thesis, Shandong University, Ji’nan, China, 2021. [Google Scholar]
  78. Schleper, C.; Jurgens, G.; Jonuscheit, M. Genomic studies of uncultivated archaea. Nat. Rev. Microbiol. 2005, 3, 479–488. [Google Scholar] [CrossRef] [PubMed]
  79. Sun, Y.; Liu, Y.; Pan, J.; Wang, F.; Li, M. Perspectives on Cultivation Strategies of Archaea. Microb. Ecol. 2019, 79, 770–784. [Google Scholar] [CrossRef] [PubMed]
Biology 12 01023 g001 550
Figure 1. Map of sampling sites.
Figure 1. Map of sampling sites.
Biology 12 01023 g001
Biology 12 01023 g002 550
Figure 2. Bar chart of the relative species abundance of archaea at the phylum level (a) and genus level (b).
Figure 2. Bar chart of the relative species abundance of archaea at the phylum level (a) and genus level (b).
Biology 12 01023 g002
Biology 12 01023 g003 550
Figure 3. Redundancy analysis showing correlations between the composition of the archaeal communities and nine environmental factors for 36 samples from 12 sampling sites. The arrows represent the measured physicochemical factors.
Figure 3. Redundancy analysis showing correlations between the composition of the archaeal communities and nine environmental factors for 36 samples from 12 sampling sites. The arrows represent the measured physicochemical factors.
Biology 12 01023 g003
Biology 12 01023 g004 550
Figure 4. LDA scores showing archaeal taxa with significant different relative abundances. (a) shows the results among sampling sites and (b) shows the results between the two lakes.
Figure 4. LDA scores showing archaeal taxa with significant different relative abundances. (a) shows the results among sampling sites and (b) shows the results between the two lakes.
Biology 12 01023 g004
Biology 12 01023 g005 550
Figure 5. Relationships between archaeal modules and soil physicochemical factors. The horizontal coordinates are physicochemical factors and the vertical coordinates indicate modules. Red squares indicate positive correlations and green squares indicate negative correlations. The darker the color, the stronger the correlation.
Figure 5. Relationships between archaeal modules and soil physicochemical factors. The horizontal coordinates are physicochemical factors and the vertical coordinates indicate modules. Red squares indicate positive correlations and green squares indicate negative correlations. The darker the color, the stronger the correlation.
Biology 12 01023 g005
Biology 12 01023 g006 550
Figure 6. The network diagram of the turquoise module. The nodes represent ASVs, and the lines show the connections between them. The larger the size of a node, the more relevant it is to other nodes.
Figure 6. The network diagram of the turquoise module. The nodes represent ASVs, and the lines show the connections between them. The larger the size of a node, the more relevant it is to other nodes.
Biology 12 01023 g006
Biology 12 01023 g007 550
Figure 7. The network diagram of the gray module. The nodes represent ASVs, and the lines show the connections between them. The larger the size of a node, the more relevant it is to other nodes. The individual submodules are connected to each other by a number of ASVs with relatively low weight values.
Figure 7. The network diagram of the gray module. The nodes represent ASVs, and the lines show the connections between them. The larger the size of a node, the more relevant it is to other nodes. The individual submodules are connected to each other by a number of ASVs with relatively low weight values.
Biology 12 01023 g007
Table
Table 1. Sample collection sites.
Table 1. Sample collection sites.
SitesLocation Profile Coordinate
NHS1Intertidal soil of the lake with birds78.9656° N 12.0729° E
NHS278.9657° N 12.0701° E
NHS378.9662° N 12.0680° E
NHX1Subtidal soil of the lake with birds78.9656° N 12.0729° E
NHX278.9657° N 12.0701° E
NHX378.9662° N 12.0680° E
XHS1Intertidal soil of the small lake78.9631° N 12.0690° E
XHS278.9638° N 12.0678° E
XHS378.9630° N 12.0698° E
XHX1Subtidal soil of the small lake78.9631° N 12.0690°E
XHX278.9638° N 12.0678° E
XHX378.9630° N 12.0698° E
Table
Table 2. Physicochemical properties of soil samples.
Table 2. Physicochemical properties of soil samples.
SitespHMC(%)TOC (%)TON (%)NO2−-N (μg∙g−1)NO3−-N (μg∙g−1)NH4+-N (μg∙g−1)SiO32−-Si
(μg∙g−1)		PO43−-P (μg∙g−1)
NHS17.08 ± 0.22 fg26.44 ± 1.36 c6.74 ± 0.81 bcde1.479 ± 0.116 bc0.424 ± 0.138 bc2.016 ± 0.740 c3.249 ± 0.262 c17.655 ± 0.635 c0.049 ± 0.016 b
NHS27.68 ± 0.15 de27.14 ± 1.77 c4.23 ± 1.11 cde1.052 ± 0.487 cd0.587 ± 0.149 b2.159 ± 0.021 bc9.597 ± 1.901 c16.852 ± 1.179 c0.219 ± 0.034 b
NHS38.11 ± 0.15 ab17.36 ± 2.37 d1.08 ± 0.32 de0.231 ± 0.064 e0.173 ± 0.079 cde0.998 ± 0.235 de5.467 ± 4.379 c20.347 ± 3.880 c0.040 ± 0.023 b
NHX17.42 ± 0.15 ef25.62 ± 2.88 c8.10 ± 1.43 bcd1.881 ± 0.328 ab0.080 ± 0.043 de0.764 ± 0.242 e2.458 ± 1.417 c14.402 ± 1.366 c0.095 ± 0.060 b
NHX27.52 ± 0.14 de15.58 ± 0.91 d3.38 ± 1.29 cde0.755 ± 0.321 de0.039 ± 0.005 e0.332 ± 0.070 e1.423 ± 0.392 c14.046 ± 0.619 c0.074 ± 0.019 b
NHX38.10 ± 0.11 abc14.06 ± 1.35 d0.23 ± 0.10 e0.132 ± 0.183 e0.080 ± 0.038 de0.364 ± 0.058 e4.613 ± 2.351 c27.733 ± 7.466 c0.045 ± 0.048 b
XHS18.30 ± 0.15 a12.36 ± 0.47 d0.87 ± 0.48 de0.207 ± 0.206 e0.196 ± 0.078 cde0.693 ± 0.187 e2.318 ± 1.518 c34.942 ± 4.522 c0.028 ± 0.022 b
XHS27.76 ± 0.02 cde31.29 ± 2.00 c6.39 ± 0.52 bcde1.736 ± 0.052 ab0.590 ± 0.042 b1.596 ± 0.158 cd1.789 ± 0.379 c19.337 ± 0.853 c0.143 ± 0.228 b
XHS37.82 ± 0.08 bcd17.60 ± 1.68 d5.19 ± 3.13 bcde1.289 ± 0.628 bcd2.230 ± 0.430 a5.975 ± 0.472 a3.342 ± 0.239 c12.103 ± 0.746 c0.471 ± 0.189 b
XHX16.81 ± 0.07 gh40.25 ± 6.82 b12.27 ± 4.39 ab2.232 ± 0.407 a0.351 ± 0.101 bcde2.767 ± 0.488 b0.524 ± 0.257 c10.295 ± 0.510 c0.017 ± 0.008 b
XHX27.02 ± 0.36 gh58.56 ± 2.92 a9.64 ± 10.03 bc0.320 ± 0.397 e0.039 ± 0.022 e0.609 ± 0.128 e23.377 ± 18.074 b191.481 ± 152.678 a0.013 ± 0.010 b
XHX36.67 ± 0.23 h62.43 ± 7.10 a18.95 ± 4.56 a0.135 ± 0.161 e0.402 ± 0.259 bcd1.612 ± 0.560 cd65.510 ± 10.494 a168.284 ± 53.541 b18.857 ± 16.914 a
In a one-way ANOVA, significant differences between study sites were assumed at p < 0.05, and were followed by Duncan’s test. The letters a, b, c, d, e, f, g and h are used to show statistically significant differences.
Table
Table 3. Monte Carlo permutation test for the physicochemical factors and archaeal community.
Table 3. Monte Carlo permutation test for the physicochemical factors and archaeal community.
RDA1RDA2r2Pr (>r)
NH4+−0.886700.462340.43390.001***
SiO32−−0.882290.470710.49110.001***
NO3−0.766900.641760.42540.001***
PO43−−0.911350.411640.14560.054.
NO2−0.716690.697390.52960.001***
TON0.89069−0.454620.13220.100.
TOC−0.972680.232150.17730.039*
MC−0.919610.392840.42640.001***
pH0.957200.289430.25270.015*
* Correlation is significant at the 0.05 level; *** correlation is significant at the 0.001 level.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Qin, Y.; Wang, N.; Zheng, L.; Li, Q.; Wang, L.; Xu, X.; Yin, X. Study of Archaeal Diversity in the Arctic Meltwater Lake Region. Biology 2023, 12, 1023. https://doi.org/10.3390/biology12071023

AMA Style

Qin Y, Wang N, Zheng L, Li Q, Wang L, Xu X, Yin X. Study of Archaeal Diversity in the Arctic Meltwater Lake Region. Biology. 2023; 12(7):1023. https://doi.org/10.3390/biology12071023

Chicago/Turabian Style

Qin, Yiling, Nengfei Wang, Li Zheng, Qinxin Li, Long Wang, Xiaoyu Xu, and Xiaofei Yin. 2023. "Study of Archaeal Diversity in the Arctic Meltwater Lake Region" Biology 12, no. 7: 1023. https://doi.org/10.3390/biology12071023

APA Style

Qin, Y., Wang, N., Zheng, L., Li, Q., Wang, L., Xu, X., & Yin, X. (2023). Study of Archaeal Diversity in the Arctic Meltwater Lake Region. Biology, 12(7), 1023. https://doi.org/10.3390/biology12071023

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