Diversity and distribution of archaea community along a stratigraphic permafrost profile from Qinghai-Tibetan plateau, China.
Permafrost occupies approximately 25% of Earth's terrestrial surface, occurring most frequently in high latitudes regions, especially those of the Northern Hemisphere [1,2]. However, alpine permafrost usually exists at low-temperature, high-elevation sites in temperate latitudes. It is well known that permafrost stores massive amounts of carbon; a recent estimate indicates that 1672 Pg of organic carbon, an amount roughly equivalent to the total carbon contained within land plants and the atmosphere [3-5], may exist in the northern permafrost region, which accounts for approximately 50% of the estimated global belowground organic carbon pool . With global warming, the permafrost is starting to thaw, with estimates of as much as 90% of the permafrost being lost by 2100 , which raises the question about the fate of carbon in thawing permafrost. It is inferred that release of carbon from permafrost to the atmosphere occurs primarily through accelerated microbial decomposition of organic matter . Previous studies showed that a large variety of microorganisms inhabit permafrost environments [8, 9]. During the thawing, the organic matter becomes more accessible to microbial degradation and results in greenhouse gas emissions , which is thought to be one of the most significant feedbacks from terrestrial ecosystems to the atmosphere, thus potentially exacerbating the greenhouse effect and further risking significant climate change .
Earlier studies showed that permafrost harbors a diverse microbial community including bacteria and archaea [10-12]. Most studies addressing permafrost microbial community were limited to the sites in Siberian permafrost [13, 14]. A research of Siberian tundra revealed that the total number of bacterial cells from the top to the bottom of the active layer (the top layer of soil that thaws in the summer and refreezes in the winter) range from 2.3 x [10.sup.9] - 1.2 x [10.sup.8] cells per gram dry soil . Archaea may constitute between 1% and 12% of the total cells in Siberian active layer soils [15,16]. The community of archaea was composed of Euryarchaeota (61%) and Crenarchaeota (39%) in the perennially frozen sediments . Methanogens, affiliated to the phylum of Euryarchaeota, are a group of archaea that produce methane under anaerobic conditions, which are ubiquitous in the biosphere and are particularly found in a variety of ecosystems such as rice paddies, lakes, hydrothermal vents, and permafrost soil and sediments . Some of them were proven to be viable inhabitants in the high latitude permafrost [18,19]. However, the uncultured methanogen cluster, Zoige cluster I (ZC-I), affiliated to Methanosarcinales within Euryarchaeota, was found to be the dominant group in the alpine permafrost wetland of Tibet plateau . Phylogenetic analyses of the methanogenic Archaea community revealed a great diversity of methanogens in the permafrost including families of Methanobacteriaceae, Methanomicrobiaceae, Methanosarcinaceae, Methanosaetaceae [16,17,19,20], methanogen group Rice cluster I (RC-I) and uncultured Rice cluster II (RC II) within the phylogenetic radiation of Methanosarcinales and Methanomicrobiales [19, 21-23] and ZC-I affiliated to Methanosarcinales , permafrost cluster I affiliated to Methanosarcinaceae , and permafrost cluster II and III affiliated to Methanosarcinales . The nonmethanogenic archaea Group 1.1b/MCG affiliated to the uncultured Crenarchaeota was mostly found in the High Arctic wetland permafrost . Even though we have already accumulated a modest amount of knowledge about the permafrost's methanogenic archaea, they continue to attract the attention of researchers due to their production of methane and the implications for the global warming.
The Qinghai-Tibetan Plateau is the largest alpine permafrost area on Earth. Earlier studies estimated the annual methane emissions from cold wetlands in the Qinghai-Tibetan Plateau at about 0.7-0.9 Tg . With the permafrost thawing due to the greenhouse effect, the organic carbon stored in permafrost would become more accessible for microbes, whose activities will dictate whether permafrost environments will be a net source or sink of greenhouse gas . However, current knowledge of the microbial community in the high-altitude permafrost is poorly understood. In this study, we described the diversity and vertical distribution of archaeal community in the permafrost profile, attempting to elucidate the composition of archaeacolonizing both the active layer and permafrost and to further characterize the methanogen community in Qinghai-Tibetan Plateau.
2. Material and Methods
2.1. Site Description and Sampling. The sampling site is located in Qinghai-Tibetan Plateau (QTP, N 38[degrees]05,38.10,, and E 99[degrees]10,05.13") with an elevation of 4300 m above sea level, where summers are mild and winters are quite cold. The mean annual temperature is approximately -5.8[degrees]C, with January temperatures ranging from -18 to -7[degrees]C and July temperatures ranging from 15 to 21[degrees]C. Significant rainfall occurs mainly in summer, while precipitation is very low in winter and spring. The study area belongs to the natural ecosystem of alpine swamp meadows with a large quantity of soil organic carbon storage. The active layer of the permafrost thaws every summer, with the thickness varying from 20 to 60 cm, while the permafrost layer occurs throughout the area to a depth of 600-700 m.
Samples were collected in June, 2013, and were obtained by digging to a depth of 65 cm. Samples were then collected aseptically from the uppermost 3-5 cm of the active layer to 63-65 cm of the permafrost layer at 10 cm intervals. The soil was placed into sterilized plastic bags for chemical analyses and 50 mL centrifuge tubes for microbiological analyses. All the samples were kept frozen in insulated containers and transported immediately to laboratory for further analyses.
2.2. Characteristics of the Permafrost Soil. The soil temperature was measured using the in situ method by inserting the pH probe into each layer of the soil. The pH of soil was determined separately on each of the soil samples for each site using a fresh soil to water ratio of 1: 5 (pH meter, Sartorius PB-10). The water content was determined as the weight loss of fresh soil dried at 105[degrees]C for 24 h. Total soil carbon content for each soil sample was determined by combustion for 16 h at 375[degrees]C. For analysis of methane concentrations in the soil, 5 g of each sample was taken with a cut-off syringe and injected into a 25 mL bottle containing 10 mL of saturated NaCl solution. The bottles were sealed by a solid butyl rubber septum and shaken for 1 h before they were used. The headspace methane concentration was determined by gas chromatography using a Packard Model 438A fitted with a flame ionization detector.
2.3. DNA Isolation. The DNA of the soil samples was extracted using the Soil DNA Isolation Kit (MP) following the manufacturer's instruction. The extracted DNA was visualized on a 1% agarose gel using ethidium bromide staining and quantified by Nanodrop (Thermo Scientiific, NanoDrop ND-2000).
2.4. Denaturing Gradient Gel Electrophoresis (DGGE). A nested-PCR approach was used to amplify V3 region of archaeal 16S rRNA. The full-length archaeal 16S rDNA gene was amplified from genomic DNA using the primer set 8f (5'-CGGTTGATCCTGCCGGA-3') and 1492r (5'GGCTACCTTGTTACGACTT-3') . The first round PCR amplification was performed in 50 [micro]L of reaction mixture containing 0.25 [micro]M of each primer, 0.2 mM dNTP, 1.5 mM MgCl2, 5 [micro]L of Taq buffer, 5 U Taq DNA polymerase (Invitrogen, USA), and 10 ng template DNA. The PCR conditions were as follows: 95[degrees]C for 5 min, 35 cycles of 50 s at 95[degrees]C, 50 s at 45[degrees]C, and 1.5 min at 72[degrees]C, followed by a final extension for 10 min at 72[degrees]C. The PCR products were gel-purified and used as a template to amplify the V3 region by using the primer set 340f-GC (57-CCCTACGGGGYGCASCAG) and 519r (TTACCGCGGCKGCTG) . The second round PCR conditions were used as follows: 95[degrees]C for 5 min, 35 cycles of 50 s at 95[degrees]C, 50 s at 45[degrees]C, and 1.5 min at 72[degrees]C, followed by a final extension for 10 min at 72[degrees]C.
DGGE was performed as described by Lazar et al.  with modifications. PCR products were applied onto 8% (wt/vol) polyacrylamide gels in 1 x TAE buffer with denaturant gradient from 20 to 40% (100% denaturant consists of 7 M urea and 40% formamide). Electrophoresis was performed at 60[degrees]C and 200 V for 6h. After electrophoresis, the gels were stained with ethidium bromide for 30 min and photographed with a UV imager (AlphaImager Mini System). DGGE bands were gel purified and reamplified with primer set 340f and 519r. The purified PCR products were cloned into vector and transformed into Escherichia coli TOP 10 for sequencing.
2.5. Clone Library Construction. Archaeal 16S rRNA genes were amplified from soil community DNA using the primers 21f (5' TTCCGGTTGATCCYGCCGGA) and 958r (5' YCCGGCGTTGAMTCCAATT); PCR reactions volume was the same with amplification of full length 16S rDNA of archaea, and PCR reactions were carried out according to DeLong  with modifications as follows: 1 cycle of 95[degrees]C for 1min, 30 cycles of 94[degrees]C for 1min, 50[degrees]C for 1min, 72[degrees]C for 2 min, and a final extension at 72[degrees]C for 10 min. The PCR products were purified as described above and cloned following the manufacturer's instructions using the pMD18T vector system (TaKaRa) together with competent E. coli JM109 cells. Randomly selected clones were sequenced.
2.6. Sequencing and Phylogenetic Analyses. The sequences were determined on an ABI 3730 automated DNA sequencer (Applied Biosystems) using the universal primer M13-47 or archaeal 21f primer. Sequences of 16S rDNA from DGGE bands and clone libraries were analyzed by using NCBI BLASTN search program to identify their putative closest phylogenetic relatives. The sequences were aligned with their relatives using Clustal W, and phylogenetic trees were constructed by the neighbor-joining method using the maximum-parsimony algorithm in MEG 4 software with 1000 bootstrap replicates. Phylotypes or operational taxonomic units (OTUs) were defined as sequences showing >97% homology to each other. Nucleotide sequences have been deposited in the GenBank database and the accession numbers are as follows: DGGE bands, KM251579-KM251601; archaeal clones, KM251602-KM251632.
2.7. Statistical Analysis. Clone library coverage was calculated according to Good . A computer program DOTUR was used to assign sequences to OTUs, and then the rarefaction curves were constructed by using the partial 16S rRNA gene sequences . Another computer program MOTHUR was used to calculate the Shannon and Simpson diversity indices, theabundance-based coverage estimator (ACE), and the bias-corrected Chao1 .
3.1. Physical and Chemical Characterization of the Permafrost Soil. The samples were collected at the end of June, 2013. The temperature along the soil profile showed a distinct decreasing gradient from the top at 18.34[degrees]C to the bottom at -0.45[degrees]C. The pH was slightly acidic, ranging from 6.43 to 6.79. Granulometric analysis revealed that the soils in the stratigraphic profile were composed of sand (799%18.73%), silt (66.84%-75.56%), and clay (11.26%-20.49%) with the average grain size ranging from 5.68 to 6.46 [micro]m. Sorting coefficients of the sample varied from 1.73 to 1.99, which indicated that the sediments in the studied sites are well-sorted (Table 1). The water content was highest in the uppermost layer of the soil because of the frequent rainfall in this area during summer months and waterlogging was often observed in the swamp. The organic carbon (TOC) in dried sediment decreased with depth from 9.01% to 1.76%, of which TOC in the topmost layer soil was significantly higher than that of other samples, suggesting an abundant organic matter accumulation in the near surface soil. The methane content was low in the topmost few centimeters but rapidly increased to the maximum value at 188.3 nmol/g in the other deeper samples (Table 1).
3.2. Diversity of Archaeal 16S rRNA Gene Sequences in the Clone Libraries. The diversity and composition of archaea community was investigated by constructing clone libraries of archaeal 16S rRNA gene fragments. Seven clone libraries were constructed from the uppermost active soil layer to the lowermost permafrost layer at 10 cm intervals. In each library, the coverage ranged from 95.2% to 97.7% (Table 2). Thirty to fifty-five sequences in each library were retrieved, and a total of 282 sequences were combined for analysis of their diversity and composition in the community. Fifteen operational taxonomic units (OTUs) were identified at the species level ([greater than or equal to]97% sequence similarity) through all the depths. Twelve and ten of them were detected in the active layer and the permafrost layer, respectively, and 7 OTUs were presented in both layers. Relatively more OTUs were observed in the depths of 3-5 cm (7 OTUs), 23-25 cm (7 OTUs), and 63-65 cm (OTUs). The Shannon and Simpson diversity indices varied from 0.79 to 1.74 and from 0.18 to 0.58, respectively, and the species richness estimators ACE and Chao 1 varied from 0 to 17 and from 4 to 10, respectively. The Shannon and Simpson diversity indices showed no significant difference, suggesting that diversity of the archaeal 16S rRNA gene libraries was similar among the layers (Table 2).
3.3. Composition of Archaeal Communities. Phylogenetic composition of archaeal 16S rRNA gene sequences revealed that they were mainly affiliated to two phyla, Crenarchaeota and Euryarchaeota, which comprise 47.2%% and 52.8% of the total sequences, respectively. The retrieved crenarchaeal sequences were classified into only one lineage, Group 1.3b/MCG-A [31, 32], which accounted for 47.2% of archaeal clone sequences. The retrieved euryarchaeal sequences could be classified into four lineages, Methanomicrobiales, Methanosarcinaceae, Methanosaetaceae, and unclassified Euryarchaeota, which accounted for 0.3%, 1.1%, 1.4%, and 50.0% of archaeal clone sequences, respectively (Figure 1).
The relative abundance and distribution of different archaeal phylogenetic lineages with the depth were shown in Figure 2. Both active layer and permafrost consisted predominantly of Crenarchaeota and Euryarchaeota. The active layer contained a proportion of Crenarchaeota at 51.2% compared to 48.8% Euryarchaeota, whereas the permafrost had a higher proportion of Euryarchaeota (58.8%) than Crenarchaeota (41.2%). Vertical distribution of OTUs along the sediment profile showed that OTU1 and OTU11, belonging to Group 1.3b/MCG-A and the unclassified group within Euryarchaeota, respectively, were widely distributed in all sediment layers. The abundance of OTU1 was gradually increased with depths, whereas OTU11 gradually decreased with depths except for the topmost three soil layers. OTU5, another lineage belonging to Group 1.3b/MCG-A within Crenarchaeota, was predominantly distributed only in the top parts of the active layers. It was very interesting that the occurrence of OTU5 was reduced sharply with depths, especially when the temperature dropped below 0.4[degrees]C; we inferred that OTU5 was probably controlled by the low temperatures. Three OTUs (OTU2, OTU3, and OTU4) were only detected in the permafrost layer and were not found in the active layer. Conversely, distributions of OTU7, OTU8, OTU9, OTU12, and OTU14 were only detected in the active layers, which suggested that the communities of archaea varied over the different depths.
3.4. Denaturing Gradient Gel Electrophoresis Analysis of 16S rRNA of Methanogenic Archaea. Denaturing gradient gel electrophoresis (DGGE) fingerprinting was used to characterize the community composition of methanogenic archaea. DGGE band patterns indicated that most DNA bands were observed in the topmost layer of the profile (3-5 cm), and similar band patterns were observed at all the depths below, which indicated the occurrence of similar methanogenic archaea communities in the lower layers of the profile. Major bands with strong intensity were commonly observed in all depths, which suggested that dominant DGGE bands represented the abundant methanogenic archaea in the community. However, several DGGE bands detected in the different layers of soil were not present in other layers, indicating depth-specific methanogenic populations. A total of 23 bands were excised and sequenced. Phylogenetic analysis showed that twelve sequences fell within the euryarchaeotal lineage, comprising three groups of ZCI cluster (12 sequences), Methanosarcinales (1 sequence), and Methanomicrobiales (8 sequences). Only two sequences were retrieved in the Miscellaneous Crenarchaeal Group (MCG) under crenarchaeotal lineage. Sequences affiliated to ZCI cluster were distributed in all depths, whereas sequences associated with Methanomicrobiales were mainly distributed in the top and middle of the permafrost sediments, but the sequences affiliated to MCG were only observed in the topmost active layer of permafrost. It was very interesting that only two bands affiliated to ZC-I cluster were observed at the depth of23-25 cm, which showed the maximum methane concentration (Figures 3 and 4 and Table 1).
Permafrost is well-known for storing massive amounts of organic carbon and is thought to be the most important natural methane emission source [3, 33]. Previous studies estimated that a total of 33.52 Pg of organic carbon and 0.7-0.9Tg of methane emission were emitted from the Qinghai-Tibetan Plateau grassland soils [25, 34]. Recent studies indicated a continuous increase in air temperatures on the QTP over the last 40 years, making them a potential source of considerable greenhouse gas emissions . Our data showed that the organic carbon content ranged from 1.61% to 3.02% at depths of 13-65 cm, which is comparable with 1.6% in a high Arctic permafrost soil from Spitsbergen reported by Hansen et al.  and 1.2-3.0% in a permafrost active layer soil from Lena Delta reported by Liebner and Wagner . However, the organic carbon content reached 9.01% in the topmost soil (3-5 cm) layers of our QTP sites. The higher organic carbon was inferred to be related to the decayed grass roots in the wetland. The methane content was much lower in both the top and bottom sediments along the permafrost profile except the subsurface at depth of 23-25 cm, where the methane content could reach a maximum value at 188.3 nmol/g. Patterns of methane content did not exhibit any trends with depth nor correlate noticeably with the physical and chemical properties of the sediments, such as organic carbon concentration, sedimentary properties, or water content. This observation was consistent with that Rivkina et al.  reported in the north-eastern Arctic tundra. The higher methane content at depth of 23-25 cm was inferred to be related with the lower redox potential and suitable substrates in this layer, which needs to be further proved.
A total of 15 OTUs of archaea were observed from the Qinghai-Tibetan Plateau permafrost soil, which presented a lower archaeal diversity when compared with other environments. Steven et al.  revealed 43 OTUs of archaea in a permafrost/ground ice core from the Canadian High Arctic. The lower diversity of archaea in QTP permafrost soil may be related to the habitat environments, which did not vary greatly as evidenced by the similar physiochemical variables shown in Table 1. Frank-Fahle et al.  reported that microbial diversity was highest in the surface layers and decreased towards the permafrost layers. Our data showed that the diversity of the archaea was similar among the layers and lower overall OTU numbers in each layer were observed. Existence of overlapping OTUs of archaeal communities from different depths suggested that the dominant archaeal community was similar over the sampling range.
The communtiy composition revealed by 16S rRNA gene clone libraries showed that Crenarchaeota and Euryarchaeota were the dominant phyla in both active layer and permafrost. A similar abundance of Crenarchaeota (51.2%) and Euryarchaeota (48.8%) was observed in the active layer. In contrast, a little higher abundance of Euryarchaeota (58.8%) and a little lower abundance of Crenarchaeota (41.2%) were observed in the permafrost layer (Figure 2). This result is quite similar with the report by Steven et al.  from the Canadian High Arctic, in which the 16S rRNA gene sequences belonging to the Crenarchaeota dominated the active layer and permafrost table horizons, while Euryarchaeota were predominant in the permafrost. However, Wilhelm et al.  reported that both active layer and permafrost consisted predominantly of Crenarchaeota at 71% and 95%, respectively; the active layer had a greater proportion of Euryarchaeota (22%) compared with permafrost (4%). The abundance of Crenarchaeota was inferred to be an increasing trend with a decreasing temperature, which suggested that members of Crenarchaeota may be particularly adapted to cold conditions . Nevertheless, our data indicated that even though both OTU1 and OTU5 belong to Group 1.3b/MCG-A within Crenarchaeota; 16.2% of OTU1 and 26.8% of OTU5 were distributed in the active layer; in comparison, 34.2% of OTU1 and 3.5% of OTU5 were distributed in the permafrost (Figure 2), suggesting that OTU1 may be more adapted to the cold environment whereas OTU5 appears more cold sensitive.
Based on 16S rRNA gene sequence analysis, both Crenarchaeota and Euryarchaeota showed a limited diversity with only one group (Group 1.3b/MCG-A) underlying Crenarchaeota and four groups (Methanomicrobiales, Methanosarcinaceae, Methanosaetaceae, and an unclassified group) underlying Euryarchaeota. A total of 97.2% sequences belonged to both Group 1.3b/MCG-A (47.2%) and unclassified Euryarchaeota (50.0%), which were predominantly comprised of OTU11 (48.9%) and OTU1 (23.4%), respectively (Figure 2). Ochsenreiter et al.  reported that Group 1.3b, an non-cultured group of Crenarchaeota, was widely distributed in different environments such as freshwater, wastewater, and soil. Sequences of the 16S rRNA gene affiliated to unclassified group within noncultured Euryarchaeota were frequently present in freshwater environments. Our data showed that Group 1.3b/MCG-A and unclassified Euryarchaeota were predominant archaea distributed over all the depths. It was very interesting that OTU5 affiliated to Group 1.3b/MCGA was overwhelmingly distributed in the active layers, and sequences belonging to OTU1 tended to increase gradually with depth, while sequences belonging to OTU11, affiliated to unclassified Euryarchaeota, were primarily distributed in the layers with much lower tempertures (Figure 2). To our knowledge, those findings have never been reported previously.
Methanogenic archaea, also known as methanogens, are an important group of Euryarchaeota that produce methane under anaerobic conditions and have proven to be viable inhabitants of permafrost [36, 38]. Our DGGE sequences analyses revealed that members of Methanosarcinales and Methanomicrobiales constituted the majority of methanogens in the soils, and the uncultured methanogen of ZC-I cluster affiliated to Methanosarcinales showed the strongest intensity of bands in DGGE and was distributed in all the depths along the permafrost profile, which indicated a dominant group of methanogens occurring in cold ecosystems. Similar observation was reported in the Zoige wetland of the Tibetan Plateau . Methanosarcinales are widespread in diverse anaerobic habitats including freshwater and marine mud and sediments, rumens, and sewage sludge digestors . They are capable of using acetate as a substrate for methanogenesis, distinguishing them from Methanomicrobiales, which can use [H.sub.2] and C[O.sub.2] as a substrate for methanogenesis . Our data showed that members of Methanomicrobiales generally showed less intense bands in DGGE and were distributed in limited layers, which indicated that the acetotrophic methanogens were relatively more abundant than the hydrogenotrophic methnogens in the permafrost. Chin and Conrad  reported when a paddy soil shifts to a low temperature, it resulted in a transient accumulation of acetate. Another report from Wagner and Pfeiffer  showed that acetate usually serves as the substrate for methanogenesis at a lower temperature. The distribution pattern of methanogens in QTP permafrost is consistent with a report from a Siberian Arctic permafrost in which Methanosarcinales were the dominant methanogens in the low temperature environment . Using the DGGE primers for amplification of methanogenic archaea, we recovered two sequences affiliated to MCG within Crenarchaeota from the topmost layer of the active permafrost, where the organic carbon content is very high (9.01%). The MCG archaea, a group of heterotrophic anaerobes, have a much wider habitat range that includes terrestrial and marine, hot and cold, and surface and subsurface environments [43, 44]. Currently, the carbon and energy sources for the MCG are unknown. Previously researches indicated that they were usually dominant in nutrient-rich environments and may utilize complex organic substrates [44-47].
In this study, we have investigated the archaeal community composition in stratigraphic soils from an alpine permafrost of Qinghai-Tibetan Plateau. Diversity of archaea was similar among all the depths, and the community was predominantly constituted by Crenarchaeota and Euryarchaeota. Each group of those two phyla has a unique distribution pattern between the active layers and the permafrost layers. The majority of methanogens belonged to Methanosarcinales and Methanomicrobiales under the phyla of Euryarchaeota. The present study will help improve our understanding of the community structure of archaea in the permafrost environment.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors acknowledge James Hurley for making a critical reading and revision of this paper. This research was supported by the National Special Research Fund, China (Gas Hydrate Resource Exploration and Production Testing Project).
 D. E. Graham, M. D. Wallenstein, T. A. Vishnivetskaya et al., "Microbes in thawing permafrost: the unknown variable in the climate change equation," The ISME Journal, vol. 6, no. 4, pp. 709-712, 2012.
 T. Zhang, R. G. Barry, K. Knowles, J. A. Heginbottom, and J. Brown, "Statistics and characteristics of permafrost and ground-ice distribution in the Northern Hemisphere," Polar Geography, vol. 23, no. 2, pp. 132-154, 1999.
 E. A. G. Schuur, J. Bockheim, J. G. Canadell et al., "Vulnerability of permafrost carbon to climate change: implications for the global carbon cycle," BioScience, vol. 58, no. 8, pp. 701-714,2008.
 R. MacKelprang, M. P. Waldrop, K. M. Deangelis et al., "Metagenomic analysis of a permafrost microbial community reveals a rapid response to thaw," Nature, vol. 480, no. 7377, pp. 368-371, 2011.
 S. A. Zimov, E. A. G. Schuur, and F. S. Chapin, "Permafrost and the global carbon budget," Science, vol. 312, no. 5780, pp. 1612-1613, 2006.
 C. Tamocai, J. G. Canadell, E. A. G. Schuur, P. Kuhry, G. Mazhitova, and S. Zimov, "Soil organic carbon pools in the Northern circumpolar permafrost region," Global Biogeochemical Cycles, vol. 23, no. 2, Article ID GB2023, 2009.
 D. M. Lawrence and A. G. Slater, "A projection of severe near-surface permafrost degradation during the 21st century," Geophysical Research Letters, vol. 32, no. 24, pp. 1-5, 2005.
 A. A. Hansen, R. A. Herbert, K. Mikkelsen et al., "Viability, diversity and composition of the bacterial community in a high Arctic permafrost soil from Spitsbergen, Northern Norway," Environmental Microbiology, vol. 9, no. 11, pp. 2870-2884, 2007
 B. Steven, W. H. Pollard, C. W. Greer, and L. G. Whyte, "Microbial diversity and activity through a permafrost/ground ice core profile from the Canadian high Arctic," Environmental Microbiology, vol. 10, no. 12, pp. 3388-3403, 2008.
 B. Steven, R. Leveille, W. H. Pollard, and L. G. Whyte, "Microbial ecology and biodiversity in permafrost," Extremophiles, vol. 10, no. 4, pp. 259-267, 2006.
 B. Steven, G. Briggs, C. P. McKay, W. H. Pollard, C. W. Greer, and L. G. Whyte, "Characterization of the microbial diversity in a permafrost sample from the Canadian high Arctic using culture-dependent and culture-independent methods," FEMS Microbiology Ecology, vol. 59, no. 2, pp. 513-523, 2007
 D. Wagner, "Microbial communities and processes in Arctic permafrost environments," in Microbiology of Extreme Soils, P. Dion and C. S. Nautiyal, Eds., pp. 133-154, Springer, Berlin, Germany, 2008.
 S. Liebner and D. Wagner, "Abundance, distribution and potential activity of methane oxidizing bacteria in permafrost soils from the Lena Delta, Siberia," Environmental Microbiology, vol. 9, no. 1, pp. 107-117, 2007
 T Shi, R. H. Reeves, D. A. Gilichinsky, and E. I. Friedmann, "Characterization of viable bacteria from siberian permafrost by 16S rDNA sequencing," Microbial Ecology, vol. 33, no. 3, pp. 169-179, 1997.
 S. Kobabe, D. Wagner, and E.-M. Pfeiffer, "Characterisation of microbial community composition of a Siberian tundra soil by fluorescence in situ hybridisation," FEMS Microbiology Ecology, vol. 50, no. 1, pp. 13-23, 2004.
 L. Hoj, R. A. Olsen, and V. L. Torsvik, "Effects of temperature on the diversity and community structure of known methanogenic groups and other archaea in High Arctic peat," The ISME Journal, vol. 2, no. 1, pp. 37-48, 2008.
 D. Wagner and S. Liebner, "Methanogenesis in Arctic permafrost habitats," in Handbook of Hydrocarbon and Lipid Microbiology, K. N. Timmis, Ed., pp. 655-623, Springer, Berlin, Germany, 2010.
 E. Rivkina, D. Gilichinsky, S. Wagener, J. Tiedje, and J. Mcgrath, "Biogeochemical activity of anaerobic microorganisms from buried permafrost sediments," Geomicrobiology Journal, vol. 15, no. 3, pp. 187-193, 1998.
 L. Ganzert, G. Jurgens, U. Munster, and D. Wagner, "Methanogenic communities in permafrost-affected soils of the Laptev Sea coast, Siberian Arctic, characterized by 16S rRNA gene fingerprints," FEMS Microbiology Ecology, vol. 59, no. 2, pp. 476-488, 2007.
 G. Zhang, J. Tian, N. Jiang, X. Guo, Y. Wang, and X. Dong, "Methanogen community in Zoige wetland of Tibetan plateau and phenotypic characterization of a dominant uncultured methanogen cluster ZC-I," Environmental Microbiology, vol. 10, no. 7, pp. 1850-1860, 2008.
 C. Erkel, M. Kuba, R. Heinhardt, and W. Liesack, "Genome of rice cluster I Archaea--the key methane producers in the rice rhizosphere," Science, vol. 313, no. 5785, pp. 370-372, 2006.
 J. B. Yavitt, N. Basiliko, M. R. Turetsky, and A. G. Hay, "Methanogenesis and methanogen diversity in three peatland types of the discontinuous permafrost zone, Boreal Western Continental Canada," Geomicrobiology Journal, vol. 23, no. 8, pp. 641-651, 2006.
 R. Grofikopf, S. Stubner, and W. Liesack, "Novel euryarchaeotal lineages detected on rice roots and in the anoxic bulk soil of flooded rice microcosms," Applied and Environmental Microbiology, vol. 64, no. 12, pp. 4983-4989, 1998.
 R. C. Wilhelm, T. D. Niederberger, C. Greer, and L. G. Whyte, "Microbial diversity of active layer and permafrost in an acidic wetland from the Canadian high arctic," Canadian Journal of Microbiology, vol. 57, no. 4, pp. 303-315, 2011.
 H. Jin, J. Wu, G. Cheng, N. Tomoko, and G. Sun, "Methane emissions from wetlands on the Qinghai-Tibet Plateau," Chinese Science Bulletin, vol. 44, no. 24, pp. 2282-2286, 1999.
 C. S. Lazar, S. L'Haridon, P. Pignet, and L. Toffin, "Archaeal populations in hypersaline sediments underlying orange microbial mats in the napoli mud volcano," Applied and Environmental Microbiology, vol. 77, no. 9, pp. 3120-3131, 2011.
 E. F. DeLong, "Archaea in coastal marine environments," Proceedings of the National Academy of Sciences of the United States of America, vol. 89, no. 12, pp. 5685-5689,1992.
 I. J. Good, "The population frequencies of species and the estimation of population parameters," Biometrika, vol. 40, pp. 237-264, 1953.
 P. D. Schloss and J. Handelsman, "Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness," Applied and Environmental Microbiology, vol. 71, no. 3, pp. 1501-1506, 2005.
 P. D. Schloss, S. L. Westcott, T. Ryabin et al., "Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities," Applied and Environmental Microbiology, vol. 75, no. 23, pp. 7537-7541, 2009.
 T. Ochsenreiter, D. Selezi, A. Quaiser, L. Bonch-Osmolovskaya, and C. Schleper, "Diversity and abundance of Crenarchaeota in terrestrial habitats studied by 16S RNA surveys and real time PCR," Environmental Microbiology, vol. 5, no. 9, pp. 787-797, 2003.
 J. Meng, J. Xu, D. Qin, Y. He, X. Xiao, and F. Wang, "Genetic and functional properties of uncultivated MCG archaea assessed by metagenome and gene expression analyses," The ISME Journal, vol. 8, no. 3, pp. 650-659, 2014.
 R. Cavicchioli, "Cold-adapted Archaea," Nature Reviews Microbiology, vol. 4, no. 5, pp. 331-343, 2006.
 W. Genxu, Q. Ju, C. Guodong, and L. Yuanmin, "Soil organic carbon pool of grassland soils on the Qinghai-Tibetan Plateau and its global implication," Science of the Total Environment, vol. 291, no. 1-3, pp. 207-217, 2002.
 X. Kang, "The features of climate change in the Qinghai-Tibetan Plateau region in the past 40 years," Journal of Glaciology and Geocryology, vol. 18, supplement, pp. 281-288, 1996.
 E. Rivkina, V. Shcherbakova, K. Laurinavichius et al., "Biogeochemistry of methane and methanogenic archaea in permafrost," FEMS Microbiology Ecology, vol. 61, no. 1, pp. 1-15, 2007
 B. A. Frank-Fahle, E. Yergeau, C. W. Greer, H. Lantuit, and D. Wagner, "Microbial functional potential and community composition in permafrost-affected soils of the NW Canadian Arctic," PLoS ONE, vol. 9, no. 1, Article ID e84761, 2014.
 D. Wagner, A. Lipski, A. Embacher, and A. Gattinger, "Methane fluxes in permafrost habitats of the Lena Delta: effects of microbial community structure and organic matter quality," Environmental Microbiology, vol. 7, no. 10, pp. 1582-1592, 2005.
 M. M. Kendall and A. D. Boone, "The order Methanosarciales," Prokaryotes, vol. 3, pp. 244-256, 2006.
 J. L. Garcia, B. Ollivier, and W. B. Whitman, "The order Methanomicrobiales," Prokaryotes, vol. 3, pp. 208-230, 2006.
 K.-J. Chin and R. Conrad, "Intermediary metabolism in methanogenic paddy soil and the influence of temperature," FEMS Microbiology Ecology, vol. 18, no. 2, pp. 85-102, 1995.
 D. Wagner and E.-M. Pfeiffer, "Two temperature optima of methane production in a typical soil of the Elbe river marshland," FEMS Microbiology Ecology, vol. 22, no. 2, pp. 145-153, 1997.
 A. P. Teske, "Microbial communities of deep marine subsurface sediments: molecular and cultivation surveys," Handbook of Environmental Chemistry, vol. 23, no. 6, pp. 357-368, 2006.
 A. Teske and K. B. S0rensen, "Uncultured archaea in deep marine subsurface sediments: have we caught them all?" The ISME Journal, vol. 2, no. 1, pp. 3-18, 2008.
 F. Inagaki, M. Suzuki, K. Takai et al., "Microbial communities associated with geological horizons in coastal subseafloor sediments from the sea of Okhotsk," Applied and Environmental Microbiology, vol. 69, no. 12, pp. 7224-7235, 2003.
 J. F Biddle, J. S. Lipp, M. A. Lever et al., "Heterotrophic Archaea dominate sedimentary subsurface ecosystems off Peru," Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 10, pp. 3846-3851, 2006.
 M. J. L. Coolen, H. Cypionka, A. M. Sass, H. Sass, and J. Overmann, "Ongoing modification of Mediterranean pleistocene sapropels mediated by prokaryotes," Science, vol. 296, no. 5577, pp. 2407-2410, 2002.
Shiping Wei, (1) Hongpeng Cui, (1) Hao He, (1) Fei Hu, (1) Xin Su, (1) and Youhai Zhu (2)
(1) School of Marine Sciences, China University of Geosciences, Beijing 100083, China
(2) Oil and Gas Survey, China Geological Survey, Beijing 100029, China
Correspondence should be addressed to Shiping Wei; firstname.lastname@example.org and Xin Su; email@example.com
Received 25 August 2014; Revised 22 October 2014; Accepted 3 November 2014; Published 25 November 2014
Academic Editor: Chuanlun Zhang
Table 1: Physiochemical properties of soil samples at different depths. Soil depth Temperature pH Sand Silt (cm) ([degrees]C) (%) (%) 3-5 18.34 6.79 16.33 72.40 13-15 3.68 6.70 7.99 75.56 23-25 1.85 6.74 13.19 69.81 33-35 0.44 6.60 8.06 71.45 43-45 -0.12 6.44 10.85 69.79 53-55 -0.24 6.57 18.73 66.84 63-65 -0.45 6.43 15.00 68.96 Standard Soil depth Clay Average grain deviation (cm) (%) size ([micro]m) ([sigma]) 3-5 11.26 5.68 1.80 13-15 16.45 6.25 1.73 23-25 17.00 6.18 1.91 33-35 20.49 6.46 1.82 43-45 19.36 6.30 1.90 53-55 14.43 5.84 1.99 63-65 16.05 6.06 1.94 Soil depth Water TOC C[H.sub.4] (cm) content (%) (%) (nmol/g) 3-5 46.52 9.01 4.1 13-15 25.56 3.02 19.9 23-25 21.35 2.43 188.3 33-35 17.18 2.43 79.3 43-45 19.23 1.76 6.7 53-55 20.91 1.61 22.5 63-65 25.54 1.96 16.1 Table 2: Sequencing information and statistical analyses of archaeal 16S rRNA gene clone libraries. Soil depth NC NO Coverage (%) ACE 3-5 cm 50 7 96.0 9.24 (73-24.6) 13-15 cm 42 5 95.2 6.94 (5.2-24.2) 23-25 cm 33 7 97.0 7.38 (70-12.2) 33-35 cm 43 4 97.7 4.78 (4.1-14.0) 43-45 cm 45 5 95.6 8 (5.6-20.8) 53-55 cm 39 5 97.4 0 (0-0) 63-65 cm 30 7 96.7 17 (8.1-98.3) Soil depth Chao 1 Shannon Simpson 3-5 cm 7.5 (7.0-15.3) 1.42 (1.2-1.7) 0.30 (0.21-0.40) 13-15 cm 6 (5.1-18.5) 1.07 (0.8-1.3) 0.43 (0.29-0.56) 23-25 cm 7 (0-7) 1.74 (1.5-2.0) 0.18 (0.11-0.25) 33-35 cm 4 (0-4) 0.79 (0.5-1.1) 0.58 (0.41-0.74) 43-45 cm 5.5 (5.0-13.3) 0.94 (0.7-1.2) 0.48 (0.35-0.60) 53-55 cm 8 (5.4-29.4) 0.99 (0.7-1.3) 0.42 (0.33-0.51) 63-65 cm 10 (7.4-30.0) 1.35 (1.0-1.7) 0.31 (0.21-0.41) Notes: NC is the number of clones in each library. NO is the number of Operational Taxonomic Unit (OTU) based on 97% nucleotide identity. Figure 2: Frequency of individual types (OTU) and relative depth-dependent distribution of groups of the archaea from the stratigraphic permafrost profile based on the gene clone libraries. Green segments indicate OTUs belonging to unclassified Euryarchaeota, blue segments indicate OTUs belonging to Group 1.3b/MCG-A, yellow segments indicate OTUs belonging to Methanosaetaceae, red segments indicate OTUs belonging to Methanosarcinaceae, and pink segments indicate OTUs belonging to Methanomicrobiales. OTU1 16.1% OTU2 0.0% OTU3 0.0% OTU4 0.0% OTU5 26.8% OTU6 4.2% OTU7 1.8% OTU8 1.8% OTU9 0.6% OTU10 0.6% OTU11 45.2% OTU12 0.6% OTU13 1.2% OTU14 0.6% OTU15 0.6% Active layer OTU1 14% OTU5 48% OTU6 4% OTU8 6% OTU9 2% OTU10 2% OTU11 24% 3-5cm OTU1 16.7% OTU5 16.7% OTU11 61.9% OTU12 2.4% OTU15 2.4% 13-15cm OTU1 18.2% OTU5 33.3% OTU6 12.1% OTU7 9.1% OTU11 18.2% OTU13 6.1% OTU14 3.0% 23-25cm OTU1 16.3% OTU5 7.0% OTU6 2.3% OTU11 74.4% 33-35cm OTU1 34.2% OTU2 0.9% OTU3 0.9% OTU4 0.9% OTU5 2.6% OTU6 1.8% OTU7 0.0% OTU8 0.0% OTU9 0.0% OTU10 0.9% OTU11 54.4% OTU12 0.0% OTU13 1.8% OTU14 0.0% OTU15 1.8% Permafrost OTU1 26.7% OTU5 4.4% OTU6 2.2% OTU11 64.4% OTU15 2.2% 43-45cm OTU1 39.5% OTU5 2.6% OTU6 2.6% OTU11 52.6% OTU15 2.6% 53-55cm OTU1 38.7% OTU2 3.2% OTU3 3.2% OTU4 3.2% OTU10 3.2% OTU11 41.9% OTU13 6.5% 63-65cm Note: Table made from pie chart.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Wei, Shiping; Cui, Hongpeng; He, Hao; Hu, Fei; Su, Xin; Zhu, Youhai|
|Date:||Jan 1, 2014|
|Previous Article:||Archaeal genome guardians give insights into eukaryotic DNA replication and damage response proteins.|
|Next Article:||The alternative route to heme in the methanogenic archaeon Methanosarcina barkeri.|