Introduction

In recent years, researchers have proposed the concept of glacial ecology1. Three key glacial ecosystems have emerged: the supraglacial, subglacial, and englacial ecosystems, and most studies have focused on supraglacial ecosystems, such as cryoconite2, snow3, soil4, and glacier streams5. Moreover, most of these studies have been performed in one single habitat. However, a few studies have examined multiple habitats6,7,8, including glacial runoff and aquatic habitats. These studies have shown that different habitats harbour different microbial communities and biodiversity, albeit with some common groups, even within the same glacier.

Notably, these results with multiple habitats have focused on aquatic environments. In glaciers, when temperatures rise above zero, large volumes of melt water are produced; this melt water first enters glacial soil and then runs into moraine lakes or streams. Glacial soil may represent an important microbial source for the entire supraglacial system since it is connected with glacial snow and moraine lakes/streams. With the continuing glacial retreat, the interactions between bacterial communities in glacial snow and in glacial soil are becoming more significant. However, the bacterial diversity and interactions between communities in glacial snow and in glacial soil proximal to glacial snow are still unclear. Moreover, information regarding unclassified sequences and functional genera in Chongce Ice Cap is still needed. This type of research is necessary to explore changes in microbial resources due to accelerated glacier retreat.

Accordingly, in this study, we aimed to compare the community structures of glacial snow and glacial soil using 454 pyrosequencing of the 16S rRNA gene of microbial samples collected from the Chongce Ice Cap. We hypothesised that snow and soil would harbour different bacterial communities owing to their different environmental features. Moreover, we evaluated certain bacteria that, to the best of our knowledge, have not been reported in other glacial environments. Our study provides a basis for the understanding of microbial ecology in different supraglacial habitats, which is particularly important because of the unprecedented speed of glacier retreat.

Results and Discussion

A total of 126,728 raw reads and 96,012 operational taxonomic units (OTUs, Table 1) were obtained from the nine samples through 454 pyrosequencing analysis. The average length of the sequences was 478 bp. Each of the nine samples contained between 12,047 and 15,798 raw reads, with OTUs ranging from 7,769 to 13,434 (Table 1).

Table 1 Observed and estimated bacterial diversity from each sample.

The rarefaction curves, based on OTUs at 3% dissimilarity, showed bacterial diversity and richness of the samples (Fig. 1). Rarefaction curves were used to determine whether sufficient sequencing depth was achieved. The glacial snow displayed saturated trends, whereas the glacial soil showed upward trends. In addition, the variations observed for the different curves were supported by the coverage (Table 1) of different samples (95.7–99.7%), indicating that sufficient sequencing depth was obtained in our study.

Figure 1
figure 1

Rarefaction analysis of the nine samples.

Rank-abundance curves indicated that the majority of the reads in glacial snow showed more dominance and fewer OTUs, whereas the reads in glacial soil showed the opposite pattern with low dominance and higher richness (Fig. 2).

Figure 2
figure 2

Rank-abundance curve based on bacterial OTUs at a dissimilarity level of 3%.

Different bacterial diversities between glacial snow and glacial soil

The rarefaction analysis of the OTUs (Fig. 1) indicated that bacterial diversity and richness were higher in glacial soil than in glacial snow. Based on the results of Shannon index analysis, bacterial diversities were 1.175–2.081 in glacial snow (mean 1.568, n = 4) and 5.354–5.892 (mean 5.586, n = 5) in glacial soil (Table 1), consistent with the above rarefaction analysis and the rank-abundance curves. The higher bacterial diversity index in glacial soil than in glacial snow was similar to the study on cryoconite holes and other glacial habitats. Takeuchi et al.9 showed that cryoconite holes may contain the highest level of biodiversity relative to other glacial habitats. In the Ecology Glacier, Antarctica10, functional diversity in Biolog Ecoplates was low in the ice samples, ranging from 1 positive response to 15, and relatively high in cryoconite holes (from 13 to 26). The mean value of Shannon index with DGGE analysis in surface ice was 0.53 (n = 6). A higher bacterial diversity was found in cryoconite holes (mean = 0.68, n = 4). The above results were consistent to the distribution of TOC (3.36 mg/L in ice, 5.55 mg/L in cryoconite), DOC (1.54 mg/L in ice, 2.77 mg/L in cryoconite), POC (1.82 mg/L in ice, 2.79 mg/L in cryoconite), total nitrogen (0.26 mg/L in ice, 1.04 mg/L in cryoconite), total phosphorus (0.05 mg/L in ice, 0.15 mg/L in cryoconite), organic phosphorus (0.02 mg/L in ice, 0.14 mg/L in cryoconite). The results indicated that changing environmental factors influenced the diversity of the microbiocenosis in both habitat types10. However, our study is absent of environmental parameters, more studies are still needed to investigate why bacterial diversity in glacial snow is lower than in glacial soil.

The values of the Simpson index were reported as greater than 0 and less than 1. Values near 0 corresponded to highly diverse or heterogeneous ecosystems, and values near 1 corresponded to ecosystems that are more homogeneous. The Simpson index in our study showed a decreasing trend (Table 1) and reflected the increasing heterogeneity of ecosystems from glacial snow to glacial soil.

Comparison of the bacterial diversity indices with other study sites

The Shannon index in our study was 1.175–2.081 in glacial snow, with an altitude of 5,800–6,100 m. In the Yala glacier, the index was 2.7–3.06, with an altitude of 5,800–6,100 m. In Laohugou No. 12, the Shannon index was 3.83–5.5811, with an altitude of 4,271–4,960 m. For glacial soil samples, in our study of the Chongce glacier, the Shannon index was 5.354–5.783, with an altitude of 5,700 m. In Laohugou No. 12, the Shannon index was 6.55, with an altitude of 4,271 m11. In Tianshan No. 1, the Shannon index was 6.35–7.55, with an altitude of 3,740–4,480 m12. Thus, bacterial diversity index at different study sites was negatively correlated with altitude. Decreased temperature with increasing altitude was the most significant environmental gradient. Effect of temperature on the altitudinal distribution of microbes was also observed by Bahram et al.13 and Wang et al.14.

Different bacterial communities between glacial snow and glacial soil

To determine the distribution and biogeography of the bacterial community, the 454 pyrosequencing data were analysed in relation to sampling locations using PCoA analysis. As shown in Fig. 3, five soil samples were clustered together and were well separated from the four snow samples. Hierarchy cluster heatmap results (see Supplementary Fig. S1) also showed two clusters: one was from five soil samples, the other was from four snow samples, confirming the above PCoA results. Both PCoA and heatmap analysis highlighted the particular differences in the bacterial communities of glacial snow and glacial soil. The above results were similar to the findings of the Svalbard glacier8 and Yala glacier6. For the Svalbard glacier, bacterial communities in the snow libraries were significantly different from those of the meltwater libraries. For the Yala glacier, microbial communities in the moraine lakes and glacier streams were completely different from those in snow.

Figure 3
figure 3

Principal co-ordinate analysis (PCoA) of bacterial communities from the nine samples based on pyrosequencing of the 16S rRNA gene.

Next, 454 pyrosequencing was used to detect 26 phyla from the nine samples (Fig. 4). The proportion of dominant sequence must be above 1% with mean value within 4 snow samples or 5 soil samples. Actinobacteria and Proteobacteria were dominant in both glacial snow and glacial soil, ubiquitous, with extensive metabolic versatility. Additionally, Deinococcus-Thermus were dominant in glacial snow. Whereas, Acidobacteria, Armatimonadetes, Bacteroidetes, Chloroflexi, Cyanobacteria, Candidate division FBP, Gemmatimonadetes, Planctomycetes, Verrucomicrobia were dominant in glacial soil (see Supplementary Table S1).

Figure 4
figure 4

Bacterial community distribution at the phylum level in the nine samples.

The sequences that could be classified were assigned to 190 genera. Acinetobacter, Pseudomonas and Hymenobacter were dominant in both glacial snow and glacial soil. Five genera (Deinococcus, Methylobacterium, Polaromonas, Rhodococcus, and Sphingomonas) were dominant among the four snow samples (Fig. 5), whereas 11 genera (Flavisolibacter, Gemmata, Kaistobacter, Lysobacter, Methylibium, Phormidium, Rubellimicrobium, Spirosoma, and Thermomonas) were dominant among the five soil samples.

Figure 5
figure 5

Bacterial community distribution at the genus level in the nine samples.

Different abundances of bacterial phyla and, to some extent, genera between glacial snow and glacial soil were attributed to environmental conditions. Microorganisms in glacial snow were from aeolian deposition, in which the microbial load from aerosol, dust, and precipitation events directly determined the microbial compositions15. Seasonal variations in temperature, nutrient concentration, and solar radiation have the potential to change microbial community composition and to promote the occurrence of species that are well adapted to the glacial snow constraints16. In our study, Deinococcus-Thermus, which can endure radiation or high temperature17 dominated among the four snow samples. Consistent with this, glacial snow is a harsh but relatively constant environment harbouring a constrained microbial community, which may be particularly well adapted to its environment18. Glacier soil, however, is more dynamic in terms of ecosystem development and suffers with primary succession. Gemmatimonadetes is among the most common bacterial groups found in soil and sediment samples from other part of the globe19. Armatimonadetes and Planctomycetes are present in other cold habitats20,21. The dominance of Chloroflexi and Cyanobacteria in glacial soil in our study was consistent with a report demonstrating that photosynthetic autotrophs appeared to be common in all of the soil samples in the Puca glacier22. These distributions have been shown to be associated with communities in later stages of ecological succession. Despite these differences, glacial snow and glacial soil shared a few common bacterial phyla with high dominance. Proteobacteria and Actinobacteria, which are the major groups in most freshwater environments around the world, were the most dominant in both glacial snow (Proteobacteria averaged 83.33% and averaged 14.22%) and glacial soil (Proteobacteria averaged 25.75% and averaged 8.12%). Three genera of Acinetobacter (0–1.3%), Pseudomonas (0–1.1%), and Hymenobacter (0–2.2%) were dominant in both glacial snow and glacial soil. These findings were similar to the report of Larose et al.8, who showed that a common core group of microbial populations existed within snow libraries and meltwater libraries on the Svalbard glacier. In previous studies, Acinetobacter was present in East Rongbuk (ER) Glacier21, Greenland Ice Sheet Project (GISP 2)23, and Taylor Glacier, Antarctica24. Hymenobacter was present in Gulkana Glacier, Alaska25 and Victoria Upper Glacier, Antarctica26, and Pseudomonas was present in Damma glacier, Central Alps of Switzerland27, soil of the Kafni glacier4, and the GISP 2 ice core28. Thus, the common bacteria present in both glacial snow and glacial soil indicated their wide distribution in the glacial environment.

Shared and rare orders/families with studies of other sites

For data at the genus and phylum levels, we compared our pyrosequencing results with those of other studies, such as the foreland of Tianshan No. 1 glacier12, the High Arctic29, the rock-water interface in an East Antarctic freshwater ecosystem, Lake Tawani30, and Laohugou No. 1211. We also compared our results with those from culture-dependent and culture-independent clone libraries as well as denaturing gradient gel electrophoresis (DGGE), such as ER16,31, Malan32,33, Palong16, soil of the Kafni glacier4, Hailuogou31, GISP 223, and the water column at Terra Nova Bay, Antarctica34. This comparison indicated that among the 26 phyla, 22 were present in most of these habitats, while four (Candidate division GAL15, Candidate division GN02, Candidate division FBP, and Candidate division MVP-21) were only recovered in our study.

At the genus level, 65 genera from a total of 189 genera were not recovered from the above-mentioned habitats; these genera were as follows: Acaryochloris, Actinoplanes, Actinotalea, Afifella, Agrobacterium, Ammoniphilus, Anaerococcus, Ardenscatena, Brevibacterium, Caloramator, Candidatus Azobacteroides, Candidatus Solibacter, Capnocytophaga, Carnobacterium, Chroococcidiopsis, Chthonomonas, Cloacibacterium, Couchioplanes, Cupriavidus, Dermacoccus, Desulfococcus, Enhydrobacter, Erythromicrobium, Erwinia, Facklamia, Filifactor, Fimbriimonas, Hylemonella, Kaistobacter, Klebsiella, Kouleothrix, Lautropia, Leptonema, Macrococcus, Magnetospirillum, Marinobacterium, Microbispora, Micrococcus, Microlunatus, Modestobacter, Mycoplana, Myroides, Neisseria, Nitriliruptor, Oscillochloris, Parvibaculum, Planomicrobium, Pseudanabaena, Riemerella, Rubrivivax, Salinimicrobium, Sanguibacter, Selenomonas, Sporichthya, Sporocytophaga, Sulfuritalea, Tepidimonas, Thauera, Thermobaculum, Thermus, Vagococcus, Veillonella, Vibrio, Virgisporangium, and Zymomonas. However, most of the above genera were not dominant in the Chongce glacial environment, except for Kaistobacter and Pseudanabaena, which were dominant in glacial soil. These latter two genera were photosynthetic bacteria, permitting their survival in a nutrient-limited environment.

Resource exploitation in the Chongce Ice Cap

In our study, about 26.28–78.67% of the sequences could not be classified into any known genus. This result was consistent with reports from a glacier foreland of the High Arctic29 and forest soil35. Thus, many of the bacterial phylotypes in the Chongce glacier may be novel. Accordingly, further studies are needed in order to gain a better understanding of the full extent of prokaryotic diversity and its distribution on Earth29.

Some known functional bacterial genera were present in the Chongce glacier, including nitrifying bacteria of Nitrospira, nitrogen-fixing bacteria of Bradyrhizobium and Rhizobium, methane-oxidising bacteria of Methylobacterium and Methylopila, and sulphur- and sulphate-reducing bacteria of Desulfococcus (see Supplementary Table S2). Bradyrhizobium and Methylobacterium were represented by more sequences in glacial snow than in glacial soil, while Nitrospira exhibited the opposite distribution pattern. The higher abundance of Methylobacterium in glacial snow than in glacial soil was probably due to a prevalence of C1 (such as methanol or formaldehyde) metabolism in snow samples. C1 compound of Formaldehyde is reported to be produced in snowpack36 and is a central intermediate in methylotrophic growth. However, we have no direct evidence of the relationship between C1 compound and Methylobacterium. Further focused studies using the mchA probe are still needed to improve our understanding of the functional bacteria in the glacier. In summary, our data suggested that, despite the potential transfer of cells from ice-marginal habitats, bacterial communities in glacial soil were essentially distinct from those in glacial snow, with most phylotypes occupying distinct niches. Some rare genera and phyla were identified in the Chongce glacier. Although total DNA was obtained from three replicates and mixed to minimise spatial variation, 454 pyrosequencing data from each sampling site were not replicated. Thus, it is impossible to identify statistically different bacterial indicators specialized within different samples using PCoA, LEfSe, NMDS, and UPGMA, in order to improve our understanding of spatial succession from glacial snow to glacial soil. Further detailed studies, such as metagenomic analysis, are needed to determine potential drivers for the differences in bacterial diversity and communities between glacial snow and glacial soil based on the chemical properties of snow and soil samples. Metatranscriptomics is also suitable for measuring changes in functional gene expression and their regulation with respect to changing environmental conditions from glacial snow to glacial soil.

Methods

Study site and sampling

The Chongce Ice Cap is located on the southern slope of the West Kunlun Mountains, which lie along the northern edge of the Tibetan Plateau close to the Taklimakan Desert. The ice cap extends over a distance of 7 km and contains two prominent domes at 6,530 and 6,374 m37,38. The elevation of the terminus is about 5,800 m, with an equilibrium line of 5,930 m39. A comprehensive field investigation in association with over-year meteorological observations has shown that the Chongce Ice Cap is a summer-accumulation type inland glacier40. There is an annual precipitation of 400–450 mm, which is mainly in the form of snowfall41. The ice surface temperature ranges from −2 to −3 °C in the summer; however, the temperature has occasionally risen to above 0 °C during the day, providing a small amount of heat to melt the ice. In addition, strong solar radiation can produce a certain amount of ablation. At the local scale, depending on the geomorphology (e.g., slopes and slope breaks), meltwater from our study glacier flowed over bare moraine. Our samples were taken after the melting season in September–October 2013, when meltwater had already inflowed into the soil in the terminus zone.

Glaciers can be divided into two parts based on the equilibrium line: the accumulation zone above the equilibrium line, and the ablation zone below the equilibrium line. The accumulation zone is where material is obtained, whereas the ablation area is where material is expended. The frontal part of the glacier tongue is called the glacier terminus42. Snow samples were collected at the terminus zone, ablation zone, balance line, and accumulation zone. Soil samples were collected at the terminus zone and at 100, 200, 300, and 400 m from the terminus zone (Fig. 6). Detailed information regarding the latitude, longitude, and altitude of each sample is given in Supplementary Table S3. Surface samples collected at 0.01 m in depth were discarded, and the underlying snow and soil samples were placed into separate sterile Whirl-Pak bags (Nasco, Salida, CA, USA). For each site, three samples were collected. The methods for extreme care and sample transport were described previously11.

Figure 6
figure 6

Location of glacial snow and glacial soil sampled in the Chongce Ice Cap.

DNA extraction

In a class 100,000 clean room, each of the melted snow samples (300 mL) was filtered through hydrophilic polyethersulfone membranes (Pall; 0.22-μm pore size) with a Millipore Labscale Tangential Flow Filtration TFF System. The microorganisms on the membranes were eluted by agitation for 2 min with a whirlpool mixer and suspended in 2.0 mL phosphate-buffered saline. Soil samples (5 g) were placed in sterile ceramic mortars after freezing in liquid nitrogen and ground until reaching a powder-like consistency. Total DNA from soil and snow samples was isolated according to the methods described by Zhou et al.43.

16S rRNA gene pyrosequencing

Equal amounts of total DNA extracted from three replicates collected at each point were mixed to minimise spatial variation. Samples were then submitted for pyrosequencing at Personal Biotechnology Limited Company (Shanghai, China). The V1–V3 region of the 16S rRNA gene was selected to construct a community library through tag pyrosequencing. The forward primer of each sample was described in detail previously11. The reverse primer was 5′-CCATCTCATCCCTGCGTGTCTCCGACGACTNNNNNNNTTACCGCGGCTGCTGGCAC-3′ (the underlined sequence is 454 Life Sciences’ primer A, and the sequence in italics is the broad-range bacterial primer 533R). NNNNNNN designates the unique seven-base barcode used to tag each polymerase chain reaction (PCR) product (see Supplementary Table S3). PCR was carried out with the following conditions: 0.4 μM each primer, 3.2 ng template DNA, 1 × PCR buffer, 0.2 U Q5 Polymerase (M0491S; New England Biolabs [NEB]), and sterile deionised water to a final volume of 25 μL. PCR was performed according to the following thermal cycling conditions: 98 °C for 4 min; 27 cycles at 98 °C for 30 s, 55 °C for 45 s, and 72 °C for 60 s; and a final extension at 72 °C for 7 min.

Fluorescent quantification and emulsion PCR were described previously11. Amplicon pyrosequencing was performed on a Roche Genome Sequencer GS FLX+ System (Roche Applied Science, Indianapolis, IN, USA). The results of the raw data have been deposited into the NCBI Sequence Read Archive (SRA) under the accession number SRP058670.

Processing of pyrosequencing data

Data were processed using the Quantitative Insights into Microbial Ecology (QIIME) pipeline44. Briefly, sequences that were less than 200 bp or greater than 1,000 bp in length, contained incorrect primer sequences, or contained more than 1 ambiguous base were discarded. Chimeric sequences were checked and deleted using the “align” procedure of the UCHIME algorithm45.

Bioinformatics analysis

Sequence reads from each sample were clustered to give similarity-based OTUs using Cluster Database at High Identity with Tolerance (CD-HIT)46 with the minimum sequence identity set to 97%. OTU rarefaction curves, Good’s coverage, and all community richness and diversity indices (ACE, Chao1, Shannon, and Simpson) were generated with Mothur47.

A heatmap was constructed using R package vegan to compare the bacterial communities of the most abundant genera in each sample48. In addition, PCoA was performed based on the weighted UniFrac distance49 in R package vegan, to demonstrate the clustering of different samples.

Additional Information

How to cite this article: Yang, G. L. et al. Differences in Bacterial Diversity and Communities Between Glacial Snow and Glacial Soil on the Chongce Ice Cap, West Kunlun Mountains. Sci. Rep. 6, 36548; doi: 10.1038/srep36548 (2016).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.