Introduction
The Tibetan Pateau, the world’s highest Plateau, contains the most glaciers within the high latitudes, with an area of 4.9×10
5 km
2 (
Shi, 2000). Glaciers and snow that cover the Plateau not only influence the global climate but also sensitively indicate global changes (
Barnett et al., 1989;
Thompson et al., 2000;
Duan et al., 2006;
Tian et al., 2006). Temporal and spatial characteristics of major ions, organic compounds, and dust particles in the snow have been studied in detail. Snow chemical compositions all exhibit strong seasonal variation and signify climatic change (
Xie et al., 2000;
Kang et al., 2002, 2004, 2007a;
Xiao et al., 2002a;
Xu et al., 2006;
Zhang et al., 2008). The snow chemistry spatial pattern decreasing from north to south implies differing atmospheric circulation (
Xiao et al., 2002a, b;
Xu et al., 2006). However, less is reported on bacterial features in snow on the plateau. Only bacteria in the snow at Xixiabangma Peak and Mt. Everest in the Himalayas have been reported (
Xiong et al., 1999;
Liu et al., 2006,
2007). Up to now, the snow bacterial response to climatic and environmental changes in different regions on the plateau has not been studied.
In recent years, bacteria in glacial ice cores over the Tibetan Plateau have received more and more attention (
Zhang et al., 2001;
Christner et al., 2003;
Xiang et al., 2005;
Zhang et al., 2007). It was found that microbial records from ice cores could potentially reflect past climatic and environmental changes (
Yao et al., 2006). A glacier is formed from snow as a result of gradual compression and transformation into ice; an understanding of bacteria in snow is important for the study of ice core records.
In this study, we investigated the bacterial abundance and diversity and the related chemical parameters of snow at the East Rongbuk Glacier and Guoqu Glacier. The two glaciers are located at the central and the southern Tibetan Plateau, respectively (Fig. 1). Our aim was to explore the seasonal and spatial variations of snow bacterial communities and their relationship with climate and the environment at the different glaciers in the Tibetan Plateau.
Sample and methods
Study area
Samples analyzed herein were collected from two snowpits (G and RBL) at the southern and the central Tibetan Plateau. Snowpit G (33°34′37.8″N, 91°10′35.3″E, 5720 m) was located at the Guoqu Glacier to the north of Mt. Geladaindong, the summit of the Tanggula Mountain, with the samples being taken on November, 2005. Snowpit RBL (27°59′N, 86°55′E, 6500 m) was situated at the East Rongbuk Glacier on the northern slope of Mt. Everest, the highest mountain in the world, with the samples being collected on April 2005.
The East Rongbuk Glacier at Mt. Everest in the Himalayas is located at the boundary region between the Indian monsoon-dominated wet and warm climate in the south and the westerly jet stream-dominated continental climate in the north. Snowfall at the glacier is caused by the moisture transported by the Indian monsoon and by local moisture from a short distance convective air mass in the summer, and by the westerly jet stream during other seasons.
The Guoqu Glacier on Mt. Geladaindong represents the northernmost extent of the summer Indian monsoon. Its summertime precipitation is derived from the moisture transported by the Indian monsoon and regional convection, with limited precipitation during the winter due to the occasional influx of westerly disturbances.
Sampling
Snowpit RBL was 1.7 m deep. Snow samples for bacterial analysis were collected continuously at 15 cm depth intervals, and at 5 cm intervals for the determination of δ18O and the concentrations of major dissolved ions.
Snowpit G was 0.86 m deep. Snow samples for bacterial abundance, δ18O and the concentrations of major dissolved ions were collected continuously at 5 cm depth intervals along the pit wall. Snow samples for bacterial diversity analysis were collected at 15 cm intervals.
Extreme care was taken at all times to ensure minimal contamination. At the beginning of sampling, the surface in contact with air was systematically scratched using a sterile spoon which was discarded. Non-particulating sterile suits, sterile gloves and masks were worn during the entire sampling process. Scoops for collecting snow were sterilized beforehand and used only once, for one sample. Snow samples forδ18O and major ions were put into Whirl Pack bags. Samples for bacterial analysis were placed in sterile Nalgene bottles. Snow samples were kept frozen until analysis in the laboratory.
Enumeration of bacterial cells
Cell abundance was analyzed using flow cytometry (Beckman Coulter, Epics Altra II). Snow meltwater (1 mL) was fixed with glutaraldehyde (final concentration: 1%) for flow cytometry analysis. SYBR Green I was applied as the nucleic acid stain (
Marie et al., 1997). Samples were run on an EPICS ALTRA II flow cytometer (Beckman Coulter), equipped with a 100-mW 488-nm water-cooled argon-ion laser (Cohenrent Inc., USA) and a standard filter set-up. An external sample injector (Harvard Apparatus PHD 2000) was employed for accurate quantification (
Jiao et al., 2005). 1 μm (diameter) fluorescent beads (Polyscience Inc.) were added to the sample for internal reference.
Measuring ofδ18O, concentration of major dissolved ions
The
δ18O values were determined using MAT-253. The concentrations of major dissolved ions (Ca
2+, Mg
2+, Na
+, Cl
1, SO
42-, NO
3-, and NH
4+) were measured using a Dionex Ion Chromatograph model 2010 (
Buck et al., 1992).
Community DNA extraction, PCR amplification and 16S rRNA gene clone library construction
Two samples from the Guoqu Glacier and six samples from the East Rongbuk Glacier were taken for molecular analysis. Snow samples were melted overnight at 4°C, and approximately 1 L meltwater was filtered through a 0.22 μm filter (Millipore). The membrane was immediately treated with 1 mL GTE buffer (25 mmol/L Tris, 10 mmol/L EDTA, 50 mmol/L Sucrose, 20 mg/mL, pH 8.0) at 37°C for 2 h. Then proteinase K, NaCl, and sodium dodecyl sulfate were added at final concentrations of 0.2 mg/mL, 0.7 mol/L, and 1%, respectively. The mixture was incubated at 65°C for 1.5 h. After extraction with phenol-chloroform-isoamyl alcohol (25∶24∶1) and chloroform-isoamyl alcohol (24∶1), DNA in the aqueous phase was precipitated with an equal volume of isopropanol overnight at –20°C. After being centrifuged and washed with 75% ethanol, the DNA was dissolved in TE buffer. To assess any contamination in the DNA extraction, a negative parallel control was established by filtering 1 L of autoclaved deionized water using the same method as described above. DNA preparations from both the samples and negative control were used as templates to amplify the bacterial 16S rRNA genes with the universal primer pair 27F (forward, 5′-AGA GTT TGA TCM TGG CTC AG-3′) and 1392R (reverse, 5′-ACG GGC GGT GTG TRC-3′) (
Brosius et al., 1978;
Christner, 2002) and La Taq polymerase (TaKaRa Co., Dalian, China). The reaction mixture (30 μL) consisted of 1 U of La Taq (TaKaRa Co., Dalian, China), 0.2 mM each dNTP, 3 μL of 10×Buffer, 0.15 mM of each primer and 1 μL (ca. 10 ng) DNA of template. The PCR program was as follows: initial incubation at 94°C for 5 min, followed by 30 cycles (94°C for 1 min, 56°C for 1 min and 72°C for 1.5 min), and then by a final extension at 72°C for 8 min.
The PCR products were purified using an agarose gel DNA purification kit (TaKaRa Co., Dalian, China), ligated into pGEM-T vector (TaKaRa Co., Dalian, China), and then transformed into E. coli DH5α. The presence of the inserts was checked by colony PCR. Clones from each library were picked randomly for re-amplification and restriction digestion by enzymes Hha I and Afa I. The digested fragments were visualized on a 3% agarose gel, and different clones were discriminated according to their RFLP patterns. One clone of each RFLP type was sequenced with 27F as the sequencing primer.
Phylogenetic analysis
All sequences obtained were checked for chimeric artifacts using the CHIMERA_CHECK program (
Maidak et al., 2001). Clones identified as potential chimeras were discarded. Sequences were calculated with DOTOUR (Schloss and Handelsman 2005) and similarities greater than 97% were grouped into one operational taxonomic unit (OTU). Sequences were assigned to genus level grouping with 80% confidence using the “Classifier” program of RDP (
Cole et al., 2005). The closest neighbors were retrieved from the NCBI (http://www.ncbi.nih.gov/ BLAST) through blasting. A phylogenetic tree including obtained OTUs and their closest relatives were constructed using MEGA software 3.1 (
Kumar et al., 2004). Neighbor-joining phylogenies were constructed from dissimilatory distances and pair-wise comparisons with the Jukes-Cantor distance model. Bootstrap analysis of 1000 replicates was performed.
Statistical analyses
Coverage of clone libraries was calculated using the equation
,
where
N is the number of clones that occurred only once (
Kemp and Aller, 2004). Diversity index (
Shannon) was calculated using the statistical program PAST (http://folk.uio.no/ohammer/past).
Nucleotide sequence accession numbers
The nucleotide sequences of partial 16S rRNA genes obtained in this study have been deposited in the GenBank database under accession numbers: DQ323081-DQ323115, EU152996-EU153042 .
Results
δ18O value and identification of seasonal snow
The
δ18O values of snow ranged from-19.95‰ to-10.68‰ in snowpit RBL, and from-23.68‰ to-9.65‰ in snowpit G. The
δ18O values of rainfall in the Tibetan Plateau are seasonally different due mainly to a precipitation “amount effect” (
Tian et al., 2001a, b,2003). In general,
δ18O is enriched during the dry non-monsoon season and is depleted during the wet monsoon season. Based on the
δ18O value, we identified the seasonal precipitation at the two snowpits. The 1.7 m-depth snowpit RBL spanned two years (2003 and 2004). Two intervals with high δ
18O values corresponded to snow deposition in the non-monsoon season, and other intervals with low
δ18O values were monsoon snow (Fig. 2(a)). The 0.86 m-depth snowpit G recorded one year (2005) of snow deposition. Snow from 0.25 to 0.55 m with the depleted
δ18O value deposit during the monsoon season (Fig. 2(b)).
Major ion concentrations
The vertical profiles of major ions Ca
2+, Mg
2+, Na
+, Cl
-, and SO
42- in snowpit G are characterised by their increased concentrations in non-monsoon snow, consistent with previous results (
Xiao et al., 2002a). Major ion concentrations are very high from the surface downward to a 0.15 m depth in the snow profile, corresponding to the non-monsoon period, suggesting the occurrence of a sandstorm event. However major ions do not show any obvious seasonal variation in snowpit RBL. The major ion concentrations are much higher in snowpit G than in RBL. The average concentration values of Ca
2+, Mg
2+, Na
+, Cl
-, and SO
42- in snowpit G are about 49, 20, 65, 13, and 22 times of those in RBL. The average content values of NO
3- and NH
4+ in snowpit G are about 4 and 3 times those in RBL.
Bacterial abundance
Bacterial abundance in snowpit RBL was from 5.7×103 to 2.3×104 cell/mL, with an obvious seasonal difference. Bacteria in the monsoon snow were more abundant than those in non-monsoon snow, with the average bacterial abundance of the monsoon and the non-monsoon seasons being 1.5×104 and 9.2×103 cell/mL, respectively.
Bacterial abundance in snowpit G ranged from 877×105 to 4.5×105 cell/mL. The layers with high Ca2+ concentrations from the post-monsoon season contained the most abundant bacteria among three seasons, with an average of 2.5×105 cell/mL. In the layers with similar Ca2+ concentrations, bacterial abundance was slightly higher during the monsoon season than during the pre-monsoon period, with the average abundance being 1.7×103 and 1.3×103 cell/mL, respectively.
16S rRNA gene clone library analysis
Four 16S rRNA gene clone libraries of the snow bacteria were contracted. RBL-m and RBL-n are the libraries representing the monsoon and the non-monsoon snow in the snowpit RBL, respectively. The RBL-m library contains clones from samples RBL-3, RBL-4, RBL-9, and RBL-10, and RBL-n contains clones from RBL-5 and RBL-7. G2 and G6 are the libraries corresponding to the monsoon and the non-monsoon seasons, respectively, in the snowpit G. In total, 349 clones were subjected to RFLP screening. The negative control did not yield any amplification products, supporting the credibility of the DNA extraction procedures. After RFLP screening, 81 unique patterns were obtained. One clone was chosen from each pattern and subjected to sequence analysis. The average sequence length was ca. 900 bp, covering V1-5 regions of 16S rRNA genes. Seventy-eight sequences were identified as normal bacterial 16S rRNA gene sequences after being checked using the CHIMERA_CHECK program. Coverage analyses showed that the sequenced clones represent 82%, 84%, 73% and 96% of the RBL-m, RBL-n, G2 and G6 libraries, indicating that the number of clones sequenced in this study provides a good inventory of the bacterial 16S rRNA gene sequences in snow from various seasons. All the sequences are associated with eight phylogenetic groups; they include the Alpha-, Beta-, and Gamma-proteobacteria, Actinobacteria, Cytophaga-Flavobacterium-Bacterioides (CFB), Firmicutes, Cyanobacteria, andTM7 candidate phylum. Of these, members of Gamma-proteobacteria dominate, and account for 46% of the total clones.
RBL-m and RBL-n libraries
Sequences in RBL-m and RBL-n are both affiliated with Alpha-, Beta-, and Gamma-proteobacteria, Actinobacteria, Firmicutes, CFB group, and Cyanobacteria. The Gamma-proteobacteria dominate in the two libraries (Fig. 2). Eight sequences, representing 77% of the total clones of the two libraries, are common in RBL-n and RBL-m. They are affiliated with genera Acidovorax, Acinetobacter, Enterobacter, Kocuria, Anoxybacillus, Chamaesiphon, and unclassified (Table 1). Among them, 35% of clones are closely related to soil bacteria, and 41% of clones are related to bacteria from the Arctic or Antarctic.
There are 14 sequences present only in RBL-m. They belong to ten genera, and four were unclassified (Table 1). Eight sequences are similar to bacteria from aquatic environments, e.g. marines, lakes, springs, and sediments in bays.
Ten sequences in different genera only occur in RBL-n. Their close relative is from Antarctica, marines, lakes, and sludges.
G2 and G6 libraries
Bacteria in the G2 library are different from those in the G6 library. Sequences in G2 are distributed in the Alpha-, Beta-, and Gamma-proteobacteria, Actinobacteria, Firmicutes, and the CFB groups (Fig. 2). Of these, members of Alpha-proteobacteria and Actinobacteria dominate, accounting for 35% and 31% of the total clones, respectively. Furthermore, they relate closely to sequences isolated from glaciers, permafrost, soil, high altitude plants, Antarctica sea water, ice and soil, sediment from the deep sea, lakes and estuaries, water from oceans, lakes and rivers, and animals (including humans) (Table 1).
The G6 library is composed of three (sub-)phyla, the CFB group, Beta-proteobacteria and Actinobacteria (Table 1). Sequences of the CFB group and Beta-proteobacteria are dominant, and account for 58% and 40% of the total clones, respectively. Sequences in G6 are associated with those from glacial and soil environments (Table 1). Seventy-six percent of the total clones’ nearest neighbors have been isolated from cold environments, i.e., high mountains, the Arctic, and Antarctica.
Discussion
Seasonal variation of bacterial abundance
Bacterial abundance at the two glaciers varies during the monsoon and the non-monsoon seasons. However, they exhibit different change trends at the East Rongbuk and Guoqu Glaciers. Bacterial abundance is correlated positively with Ca2+, Mg2+, Na+, Cl-, and SO42- concentrations (r=0.75, 0.73, 0.62, 0.72, 0.65 and 0.72, P<0.01) at the Guoqu Glacier. The bacterial abundance is higher in the non-monsoon season than in the monsoon season. In contrast, bacterial abundance does not show any obvious relationship with ion concentrations at the East Rongbuk Glacier. Meanwhile, the snow contains more bacteria in the monsoon season than in the non-monsoon season.
It has been reported that bacteria in snow are supplied by the atmosphere by mineral particles, and some psychrophilic or psychrotrophic bacteria grow rapidly after deposition (
Segawa et al., 2005). Major ions in the snow at the Guoqu Glacier originate mainly from the desert to the north. The ion concentrations indicate the amount of input from the atmosphere (
Xiao et al., 2002a). Ca
2+ concentrations in the snow provide an index of dust concentrations in the Tibetan Plateau (
Yao et al., 2004). The significant correlation between bacterial abundance and ion concentration indicates that snow bacteria in the Guoqu Glacier are strongly influenced by the dust input. Sandstorms often happen in the region during the non-monsoon season. A large number of dust resulting from sandstorm events transports a relatively great amount of bacteria into the glacier in the non-monsoon season. However, bacterial abundance is slightly lower in snow deposited in the non-monsoon season than in the monsoon season. It is suggested that bacteria grow weakly during the monsoon season. Therefore, both the bacterial amount supplied from the atmosphere and the growth conditions impact on the bacterial abundance in the Guoqu Glacier, with the supplied amount being the key role.
At the East Rongbuk Glacier, snow bacterial abundance in the monsoon season (average 15.2 cell/mL) is higher than in the non-monsoon season (average 9.2 cell/mL). The ion concentrations in the snow do not show any seasonal variation. It is suggested that the abundant bacteria in the monsoon season are not transported by dust. Bacterial growth after deposition is the possible reason causing such high bacterial abundance during the monsoon season. The automatic weather station at the sample site in our study recorded that the daily maximum temperature of most days in June and July 2005 was above 0°C (
Xie et al., 2007). The temperature makes bacterial growth possible. However, the bacterial growth conditions are currently not quite clear and await further investigation. Based on our data, it could be deduced that bacterial growth after deposition plays a more important role on the seasonal variation of bacterial abundance.
Seasonal variation of bacterial diversity
Snow bacteria at the Guoqu Glacier show significant seasonal differences. Firstly, bacterial diversity is greater during the monsoon season than during the non-monsoon season. Bacteria found in the samples from the monsoon season belong to 7 phyla and 15 genera, and 2 were unclassified. However, bacteria from the non-monsoon season belong to only 3 phyla and 3 genera, and 2 were unclassified. Secondly, the bacterial components of the two seasons are distinctly different; the CFB group dominates during the non-monsoon season, but contributes only a minor portion during the monsoon season. Beta-proteobacteria account for nearly half of the total in the non-monsoon season samples, but account for only 4% in the monsoon season samples. Alpha-proteobacter and Actinobacteria are the dominant components in the monsoon season. Thirdly, bacteria from different seasons were identified to originate from different environments. Bacteria in the non-monsoon season are characterized by having originated from two kinds of environments, cold environments and soil. Bacteria in the monsoon season are multifarious in origin, including marine rivers, lakes, plants, animals and humans, as well as cold environments and soil.
Snow bacteria at the East Rongbuk Glacier also exhibit seasonal variations. But the bacterial difference between the two seasons are not so distinct as that in the Guoqu Glacier. The bacterial components in the two seasons are similar (Fig. 2). The bacterial diversity of the Gamma-proteobacter subgroup decreases in the non-monsoon season, while the bacterial diversities of the CFB and Cyanobacteria groups increase. Both the monsoon and the non-monsoon snow contain their own genera, but the major clones are common in both seasons. The two seasons’ bacteria originate from similar environments, i.e. soil, marine lakes, the Arctic and Antarctic.
The seasonal difference in bacterial diversity is likely influenced by the atmospheric circulation and the temperature. Snow deposited in the monsoon and the non-monsoon seasons comes from different sources at the two glaciers. The summer Indian monsoon and the local convective water vapor contribute to the precipitation during the monsoon season, but the westerlies carry water vapor during the non-monsoon season (
Kang et al., 2007a, b). Differing atmospheric circulation carries various kinds of bacteria which are later loaded in the glacier. Otherwise, relatively high temperatures and wet climatic conditions during the monsoon season could result in the survival of more bacteria. Under the influence of atmospheric circulation and temperature, the seasonal variation in bacterial diversity at the East Rongbuk Glacier is different from that at the Guoqu Glacier. It has been reported that marine moisture is westerly transported to the middle Himalayas from the North Atlantic during winter (
Tian et al., 2005). Therefore, marine vapors could reach the East Rongbuk Glacier during both the monsoon and the non-monsoon seasons, but come to the Guoqu Glacier only during the monsoon season. The seasonal variations in temperature are also different at the two glaciers. The daily average temperature from the nearest meteorological station, Dingri Station, during the monsoon season is similar to that during the non-monsoon season, being at 18.1 and 10.7°C, respectively, while the temperatures measured from Anquo Station, which is the nearest station to the Guoqu Glacier, are vastly different during the two seasons, at 7.1 and-5.3°C. The distinct differences in atmospheric cycling and the temperatures between the monsoon and the non-monsoon seasons cause the significant difference in bacterial diversity at the Guoqu Glacier.
Spatial variation of snow bacteria
Three genera identified occur in both the East Rongbuk Glacier and the Guoqu Glacier; they include
Brevibacterium in
Actinobacteria,
Sphingomonas in
Alpha-proteobacteria, and
Polaromonas in
Gamma-proteobacteria. These genera are widely distributed in various glaciers on the Tibetan Plateau, such as the Manlan Glacier (
Xiang et al., 2004), the Guliya Glacier (
Christner et al., 2003), and the Puruogangri Glacier (
Zhang et al., 2006). They were also detected in Arctic and Antarctic glaciers (
Christner et al., 2000;
Miteva et al., 2004;
Miteva, 2007).
Except for the three common genera, the other bacteria in the two glaciers are affiliated with different genera and show different features. Sixteen of 46 sequences in the Guoqu Glacier are similar to soil bacteria. However, only 5 of 32 sequences in the East Rongbuk Glacier are relative to the bacteria in soil environments. Several genera, such as the
Stappia,
Roseobacter,
Ruegeria,
Rhodoferax at the East Rongbuk Glacier, are connected with those from marine environments (
Finneran et al., 2003;
Moran et al., 2003;
Kim et al., 2006). Bacterial diversity variation at the two glaciers indicates that snow bacteria are responsive to environmental conditions. The Guoqu Glacier is located at the central plateau and near the deserts, and thus receives relatively much continental dust. The East Rongbuk Glacier is located at the southern part and is strongly influenced by moisture from the Indian Ocean. The spatial variation of snow bacteria observed herein is in line with the snow chemistry. The spatial patterns of snow chemistry demonstrate that the continental dust originates mainly from the deserts to the north, while monsoon sources of the moisture exert greater influence on snow chemistry on the southern plateau (
Xiao et al., 2002a). Our data suggest that snow bacteria could be an index of climatic and environmental conditions, as good as other chemical parameters in the snow.
Conclusions
Abundant snow bacteria are present at the East Rongbuk Glacier and the Guoqu Glacier, ranging in abundance from 877×105 to 4.5×105 cell/mL. Both the bacterial abundance and the diversity show various patterns during the monsoon and the non-monsoon seasons at the two glaciers. At the Guoqu Glacier, more bacteria are transported in the non-monsoon season than in the monsoon season, with significant seasonal variation being observed in bacterial diversity. However, at the East Rongbuk Glacier, bacterial abundance is higher in the monsoon season than in the non-monsoon season due to preferential bacterial growth during the monsoon season. Bacterial differences between the two seasons are not so distinct as that seen in the Guoqu Glacier. The seasonal difference is likely influenced by various atmospheric circulations and temperature. Bacterial diversity at the two glaciers exhibits different features in response to various environments.
Higher Education Press and Springer-Verlag Berlin Heidelberg