Introduction
Ice cores are long-term, chronological repositories of microorganisms. The assemblages of airborne microorganisms immured chronologically within glacial ice are species that were atmospheric constituents originating from a variety of ecological sources at defined ages in history (Reference Priscu, Christner and BullPriscu and Christner, 2004).
With the progress of biological techniques and the increase of interest in the diversity and survival mechanism of microorganisms in frozen samples, studies of diversity, viability and physiology of the organisms in glacial ice became feasible. Although viruses, bacteria, fungi, pollen grains and seeds have been identified in glacial ice (Reference Abyzov and FriedmannAbyzov, 1993; Reference Dancer, Shears and PlattDancer and others, 1997; Reference Abyzov, Mitskevich and PoglazovaAbyzov and others, 1998; Reference CastelloCastello and others, 1999; Reference Christner, Mosley-Thompson, Thompson, Zagorodnov, Sandman and ReeveChristner and others, 2000; Reference Zhang, Yao, Ma and WangZhang and others, 2001), the detection, isolation and characterization of microorganisms from glaciers are still relatively novel enterprises.
A number of studies have focused on DNA diversity, bacterial recovery or longevity from glaciers (Reference Christner, Mosley-Thompson, Thompson and ReeveChristner and others, 2001; Reference Zhang, Yao, Ma and WangZhang and others, 2001; Reference Sheridan, Miteva and BrenchleySheridan and others, 2003), but few attempts have been made to determine the structure of bacterial populations preserved in ice deposited under diverse environments, and there are no documented reports of the changes of the bacterial population profile at different depths of the ice core. However, how the structure of microbial populations entrapped in ice cores at defined climatic and environmental boundaries changed and how these changes are related to the corresponding environmental conditions are interesting questions (Reference Castello and RogersCastello and Rogers, 2005).
The denaturing gradient gel electrophoresis (DGGE) DNA fingerprinting technique offers a more rapid and comprehensive means to simultaneously survey various bacterial populations, allowing the comparison of total bacterial communities of different samples through pattern analysis, and at the same time describing the community structure intuitively. Since DGGE was introduced into microbial ecology by Reference Muyzer, de Waal and UitterlindenMuyzer and others (1993), many studies have been carried out to obtain outlines of the bacterial communities associated with environmental perturbations or seasonal, spatial and geographical variability (Reference Ferris, Muyzer and WardFerris and others, 1996; Reference IwamotoIwamoto and others, 2000; Reference Riemann, Steward and AzamRiemann and others, 2000). Here, we examine the bacterial DNA profile of 13 selected Puruogangri (Tibetan Plateau) ice-core samples along its 83.45 m length by DGGE with the following three objectives: (1) to investigate the abundance and diversity of organisms in this sample and make comparisons with studies of geographically different ice cores and different habitats; (2) to evaluate changes in the microbial community structure along its length, and relations between the changes and environmental fluctuations; and (3) to describe what are the dominant indices in this glacier ice and how the organisms exist in situ.
Materials and Methods
Site feature
The Puruogangri ice field (33˚44’ –34˚04’ N, 89˚00’–89˚20’ E; 5970 ma.s.l.) is the largest glacier on the Tibetan Plateau, with a total area of 422.58 km2, and is located in the middle of the plateau on the western end of the Tanggula mountain range (Fig. 1). The ice field is surrounded by lakes and sand dunes that lap onto its margins, and the sand is most mobile in the spring. Its vertical relief facilitates the transportation of fragments of the sparse vegetation in the area to the top of the ice field (Reference YaoYao, 2000). In October 2000, the average temperature on the glacier surface (6000ma.s.l.) was −14.6˚C and the highest temperature was −3.6˚C, whereas the temperature at 80 m in the borehole was −7˚C.
Sampling
In October 2000, an ice core was drilled from the flat part (33˚5’057"N, 89˚06’59" E; 6000ma.s.l.) of the Puruogangri ice field, extending to 83.45 m below the surface. Without thawing, sections of the ice core were sealed in plastic tubes which were then placed inside aluminum-covered cardboard tubes. These were immediately transported to the Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, and kept in an ice-core storeroom (−20˚C) prior to sampling. Thirteen samples with a vertical length of about 15 cm were selected at different depths, as shown in Table 1.
To aseptically sample the ice core for DNA extraction, the samples were manipulated in an ultraviolet (UV)-sterilized positive-pressure laminar-flow hood and placed in a biologically clean room according to Reference Christner, Mikucki, Foreman, Denson and PriscuChristner and others (2005). The 5 mm outermost layer of ice was scraped with a newly autoclaved sawtoothed knife, and the inner 5 mm was sliced with another autoclaved sawtoothed knife. Sections of the scraped ice cores were held with sterile forceps within the laminar-flow hood and completely rinsed with −20˚C 0.2 μm filtered 95% ethanol, and then rinsed with 500 mL of 4˚C sterile, deionized water (Milli-Q Element System, Millipore, USA). The remaining ice was placed in sterile containers with a lid and allowed to thaw completely at 4˚C. The control core was processed and melted following an identical protocol to those of the samples. Ten mL of the newly melted subsamples prepared as above were used for the analysis of ions using an ion chromatography system (DX 600, Dionex Corp., USA). A control core constructed with Milli-Q-purified redistilled water was processed by methods identical with those of the sample.
Extraction of genomic DNA
The melted ice samples (∼200 mL) were filtered with a 0.2 μm pore-size autoclaved filter (diameter 25 mm; Millipore). Total DNA was extracted from bacterial cells trapped on the filter, which was suspended in 1 mL of TE buffer (pH 8.5), supplemented with 0.13 g of lysozyme and 20 μL proteinase K (20 mg mL-1), and incubated at 37˚C for 30min. After that, the resulting suspensions were mixed with 80 μL SDS (20%) and incubated at 65˚C for 2 hours. DNA was extracted with equal volumes of phenol-chloro-form-isoamyl alcohol (25:24:1, vol/vol/vol) and chloro-form-isoamyl alcohol (24:1, vol/vol). The supernatant fluid was transferred to another Eppendorf tube and precipitated with 0.7 volume of cold isopropanol and 0.1 volume of 3 M sodium acetate (pH 5.8) at −20˚C overnight, and centrifuged at >16 000g for 30min at 4˚C. The DNA pellets were washed with 70% cold ethanol and dried. Finally, the dried pellets were resuspended in a TE buffer and pooled to give a final volume of 30 μL.
PCR and DGGE
DGGE was performed as described by Reference Muyzer, de Waal and UitterlindenMuyzer and others (1993). PCR amplification was performed in a thermocycler (Icycler, Bio-RAD) in a 25 μL reaction mixture. The V3 hypervariable region of bacterial 16S rDNA fragments was amplified by using primers 357F-GC (5’CGCCCGCCGCGC-GCGGCGGGCGGGGCGGGGGCACGGGGGGCC-TACGGGAGGCAGCAG-3’) and 518R (5’-GTATTACCGCG-GCTGCTGG-3’) as the PCR primers (Reference Muyzer, de Waal and UitterlindenMuyzer and others, 1993). Each amplification reaction mixture (25 μL) consisted of 0.5 U of Ex Taq DNA polymerase (Takara Shuzo, Otsu, Japan), 50 μg of genomic DNA solution, 2.5 μL of 10 x PCR buffer,each primer at a concentration of 0.25 μM, and a mixture containing deoxynucleoside triphosphate at a concentration of 100 μM in a tube. A touchdown program (Reference Muyzer, de Waal and UitterlindenMuyzer and others 1993; Reference Ogino, Koshikawa, Nakahara and UchiyamaOgino and others, 2001) was implemented as follows: firstly there was an initial denaturation at 94˚C for 5 min followed by 30 cycles at 94°C for 1 min, an annealing temperature for 1 min, and 72°C for 1 min, and finally the reaction mixture was kept at 72°C for 7 min. During the reaction cycle, the annealing temperature was decreased by 1°C from 65 to 56°C every second cycle in the first 20 cycles. The annealing temperature was 55°C in the last 10 cycles. The amplified product was analyzed by electrophoresis on 1.5% agarose gels and quantified by UV spectrophotometer (UV 7501, Techcomp) for DGGE; 250 ng of purified amplicons was used. DGGE was performed using a D-Code system (Bio-Rad Laboratories, Inc., Hercules, CA).
PCR products were separated on a 1.5 mm thick vertical gel containing 8% (wt/vol.) polyacrylamide and a linear gradient of the denaturants urea and formamide, increasing from 30% at the top of the gel to 60% at the bottom. Here, 100% denaturant is defined as 7M urea and 40% (vol/vol) formamide. DGGE was conducted with 0.5 x TAE buffer (1 x TAE buffer is 0.04M Tris base, 0.02 M sodium acetate and 10 mm EDTA (pH 7.4)) at a constant voltage of 35 V at 60°C for 18 hours. The gel was stained with SYBR Green-I (Molecular Probes, Eugene, OR, USA) and photographed on a UV transilluminator (UVP, Bioimaging system).
Statistical analysis of DGGE banding patterns
To determine the information content of the banding patterns in terms of their structural diversity, the DGGE profile was analyzed in the following two ways:
-
1. The similarities of all possible pairs of gel tracks were calculated, and then a cluster analysis of the matrix of similarity values and visualization in a dendrogram were performed. Bacterial population similarity between samples was compared using the Dice-Sorensen similarity index (Reference PielouPielou, 1966) with MultiVariate Statistical Package (MVSP) version 3.12h (GeoMem, Blairgowrie, UK). The Dice coefficient SD = 2j/(a+b), where j is the number of bands common to both samples, a is the number of bands in sample A, and b is the number of bands in sample B. This number is then multiplied by 100 to obtain the percentage similarity. A value of 0 indicates that the samples are completely different, whereas a value of 100 indicates that they are identical.
-
2. Bacterial diversity was estimated with the Shannon-Weaver index (H) from the DGGE banding pattern of the samples, which made it possible to determine changes in species diversity. Then the Shannon-Weaver index was performed with the Labworks 4.0 software package (UVP Bioimaging system) based on the number and relative intensities of bands on a gel strip. For banding profile comparisons, it was assumed that single DGGE bands contained a single taxon. The Shannon-Weaver index of general diversity, H (Reference Shannon and WeaverShannon and Weaver, 1963), was calculated using the following function:
where P is the importance probability of the bands in a track. H was calculated on the basis of the bands on the gel tracks using the intensity of the bands as judged by peak heights in the densitometric curves. The importance probability Pi was calculated as:
where ni is the height of a peak and N is the sum of all the peak heights of the bands in the densitometry profile.
The Simpson index of dominance, S (Reference SimpsonSimpson, 1949), was calculated using the following function:
Results
The DGGE profiles of the bacterial DNA from ice-core samples taken at different depths are shown in Figure 2. The number, precise position and intensity of the bands reflect the number, species and relative abundance of dominant rDNA types in the samples and thus allow a comparison of bacterial DNA structures with each other. Each lane of the DGGE profile contained many bands, numbering approximately 17–29 (Fig. 2), which indicated that a great diversity of bacterial species has evolved in the Puruogangri ice core, but the patterns of diversity did not seem to be related to sample depth. The band patterns appeared to be different for the different depths. However, some bands appear to be common to all samples of different depths, and several intensive bands were clearly visible in some samples, which demonstrated that the dominant species are not limited to one or two.
Dendrograms were constructed from the matrix of similarity values of the DGGE banding patterns based on band existence (Fig. 3). Differences between the samples along the length were higher, as indicated by the low similarity indices (0–43%). Three distinct clusters were recognized, which corresponded to the microbial DNA structure in each of the three periods when snow was deposited under different environmental conditions. The DGGE profile of samples 1–6 clustered together; samples 7, 8, 10, 1 1 , 12 and 13 clustered into another major group; and sample 9 clustered separately. The closer sections of the ice core show higher similarity of bacterial DNA structure, and the highest similarity was found between samples 4 and 6 (43%).
The concentration of dominance (S) and the general index of biodiversity (H) were calculated on the basis of the number and relative intensities of bands on a gel track (Table 1). Biodiversity as well as concentration of dominance varied. The Shannon index (H) fluctuated within the range 2.46−3.17 (mean 2.91; SD 0.25; n=13). These results indicated bacterial species diversity in the ice core was high and fluctuated periodically, with three peaks in values along the whole length (Fig. 4). The concentration of dominance (mean 0.062; SD 0.022; n = 13) fluctuated within the range 0.047−0.112. A negative correlation (R = 0.91; P<0.01) was found between diversity index and concentration of dominance (Fig. 4).
Concentrations of Ca2+, Na+, Mg2+ and K+ at 13 increasing depth samples were investigated, and the variation in trends of Ca2+ concentration (mean 1.7 x 103; SD 0.9 x 103; n = 13) coincides with those of Na+, Mg2+ and K+ (data not shown). A significant, positive correlation (R = 0.71; P<0.01) was found between H and Ca2+ concentration (Fig. 4).
Discussion and Conclusions
There have been many studies of bacterial diversity in ice cores using both culture-dependent and culture-independent methods. However, this is the first attempt to apply a DGGE DNA fingerprinting technique to study changes of bacterial DNA structure in an ice core from a non-polar ice core.
DGGE was introduced into microbial ecology to obtain an overview of the structural diversity of microbial communities (Reference Muyzer, de Waal and UitterlindenMuyzer and others, 1993). It is based on the separation of ribosomal gene sequences directly amplified from community DNA by using conserved primers on a denaturing gel according to their melting point. DGGE analysis with similarity and diversity algorithms proved to be a very powerful tool to monitor changes in the structure of the microbial community (Reference Ferris and WardFerris and Ward, 1997; Reference Eichner, Erb, Timmis and Wagner-DoblerEichner and others, 1999; Reference IwamotoIwamoto and others, 2000). This procedure allowed one for the first time to directly identify the presence and relative abundance of different species and to profile microbial populations in both a qualitative and a semi-quantitative way. In this study, we investigated the changes in the bacterial DNA structure of the Puruogangri ice core by DGGE, and the relationship between bacterial diversity and Ca2+ concentration was also explored.
DGGE profiles indicated that bacterial species diversity is high in the Puruogangri ice core. The same result was verified by genomic DNA amplification and the recovery of viable bacteria in the Tibetan Plateau glacier (Reference Zhang, Yao and MaZhang and others, 2002, Reference Zhang, Ma, Yao and Zhang2003; Reference XiangXiang and others, 2004). In general, ice-core samples from non-polar, low-latitude, high-altitude glaciers contain a greater number, as well as more diversity, of bacterial species than those from polar ice. This can most probably be attributed to the proximity of major biological ecosystems (Reference Abyzov, Mitskevich and PoglazovaAbyzov and others, 1998; Reference PriscuPriscu and others, 1999; Reference Christner, Mosley-Thompson, Thompson, Zagorodnov, Sandman and ReeveChristner and others, 2000; Reference Zhang, Ma, Yao and ZhangZhang and others, 2003). It can be hypothesized that the monsoon, local and global environmental conditions and surrounding ecosystems contributing biological particles to a particular air mass influence the concentration and diversity of airborne microorganisms (Reference Christner, Mosley-Thompson, Thompson and ReeveChristner and others, 2003; Reference Segawa, Miyamoto, Ushida, Agata, Okada and KohshimaSegawa and others, 2005; Reference Xiang, Yao, An, Xu and WangXiang and others, 2005). Puruogangri glacier is surrounded by lakes and sand dunes, which might contribute to the higher bacterial diversity in this glacier. Increased microbial deposition occurred in parts of the glacier that were contiguous to environments that supply airborne rock grains and soils (Reference Christner, Mosley-Thompson, Thompson, Zagorodnov, Sandman and ReeveChristner and others, 2000), which served to transport and protect attached microorganisms. In addition, no relationship was perceived between DGGE pattern and depth, which probably implied that the microorganisms were deposited in the glacier episodically under diverse climatic conditions and at different historical times. This testified that perennial ice was a long-term, chronological repository of microbial assemblages corresponding to prevailing climate and environment conditions (Reference Priscu, Christner and BullPriscu and Christner, 2004; Reference Castello and RogersCastello and Rogers, 2005).
The commonality of some intensive bands in all samples (Fig. 2) may suggest that the corresponding bacterium was an endemic species and that it grew on the glacial surface. Reference Christner, Mosley-Thompson, Thompson and ReeveChristner and others (2003) found cryoconite holes might serve as biological refuges that, on glacial melting, can repopulate the local environments. However, the fact that these bacteria were only detected in a specific sample also suggests that they were occasionally supplied from the atmosphere and were not indigenous members of the snow microbiota in this region. In addition, several intensive bands were observable in some lanes (Fig. 2), which probably implied that the psychrophilic bacterial species having adapted to the unique habitat were not limited to one or two.
Community diversity is a key concept in ecology, and its quantification is fundamental for analyzing such phenomena as succession, colonization and response to disturbances. Changes in bacterial diversity of ancient organisms from an ice core can be used to formulate correlations between microbial distributions and environmental characteristics (climatic conditions, volcanism, etc.) at each of the time points (Reference Castello and RogersCastello and Rogers, 2005). Here we applied the Shannon index (Reference Shannon and WeaverShannon and Weaver, 1963), based on the number and intensity of DGGE bands, to reflect the changes in bacterial diversity within the microbial assemblages in the ice core with no need for cultivation. The biodiversity index H’ and concentration of dominance combines the relative abundance of species and the total specific richness, thus reflecting the distribution of the dominance amongst the major species present in the samples studied. High H’ values (mean = 2.91) revealed that the bacterial diversity in the Puruogangri ice core was high but the variability was not related to the sample depth (Fig. 4). The lower index of dominance (Table 1) compared to another active community (Reference Ampe and MiambiAmpe and Miambi, 2000) is further evidence that microbial assemblages in glacial ice are a relatively silent world. Interestingly, a significantly negative correlation between H’ and S values was found (Fig. 4). This reflects the inherent dynamics of the bacterial community.
Noticeably, the Shannon index, H’, displayed a positive correlation with the Ca2+ concentration (Fig. 4), which is a good proxy of the dust in the ice core from the Tibetan Plateau (Reference Sheng, Yao and LiSheng and others, 1999; Reference Wu, Yao and ThompsonWu and others, 2004). Variations of Ca2+ concentration are generally related to climatic changes: high Ca2+ concentration periods coincided with cold periods, and low Ca2+ concentration periods with warm periods (Reference Sheng, Yao and LiSheng and others, 1999; Reference Wu, Yao and ThompsonWu and others, 2004; Reference Wang, Yao, Thompson and DavisWang and others, 2006). Thus we conclude that more bacterial species are deposited in the glacier during cold periods, and fewer during warm periods. Logically, the intensified atmospheric circulation during cold climatic periods helped increase the concentration of airborne dust particles. Therefore, more organisms attached to the macroparticles were moved along the glacier surface during cold periods. Such conditions also decreased desiccation and increased the likelihood of bacterial survival. Moreover, dust particles protected microorganisms and influenced the amount as well as the physiological state of the microbes in ice (Reference PriscuPriscu and others, 1998; Reference Christner, Mosley-Thompson, Thompson, Zagorodnov, Sandman and ReeveChristner and others, 2000). For example, dust may form a cryoconite hole on a glacier, as it absorbs solar radiation and promotes melting of the ice beneath it, where the activities of algae and bacteria might be supported by photosynthesis and bacterial decomposition of old algae and other organic particles (Reference Takeuchi, Kohshima, Shiraiwa and KubotaTakeuchi and others, 2001, Reference Takeuchi, Kohshima and Segawa2003). The results of this investigation indicated that the dust content might be more vital than other factors in determining bacterial species diversity in the Puruogangri glacial ice.
Studies indicated the bacterial diversity index (H’) in other habitats, the potential originating sources of microbes entrapped in the Puruogangri ice core, was varied. The Shannon index ranged from 3.72–4.44 in horticultural soil of southern Italy to 2.13 in an Italian agricultural soil and 2.87 in a Belgian soil (Reference Andreoni, Cavalca and RaoAndreoni and others, 2004; Reference Crecchio, Gelsomino and AmbrosoliCrecchio and others, 2004). It ranged from 1.20 to 1.34 in Biwa lake (Japan) sediment (Reference Koizumi, Kojima and FukuiKoizumi and others, 2003) and from 0.887 to 1.271 in a North Sea marine environment (Reference DavidDavid, 2004). Shannon index analysis of the DGGE pattern revealed that the bacterial DNA diversity in the Puruogangri ice core was higher, ranging from 2.46 to 3.17 (mean 2.91; SD 0.25; n = 13). The index value in the Puruogangri ice core is comparative to that of soil. This may imply that the microorganisms in the Puruogangri glacier mainly originated from a soil ecosystem, which is consistent with the result derived by Ca2+ analysis.
Another notable point is that periodic fluctuation of the Shannon index of bacterial DNA structure was observed. This probably reflected a cycle of climatic changes (Fig. 4). Unfortunately, due to the lack of climate analysis and dating data, we cannot elicit a precise relationship between bacterial DNA diversity and climatic factors or geologic events in a specific time frame in the Puruogangri ice core. But following an interdisciplinary perspective could be promising to verify the changes of climate.
Clustering analysis performed here demonstrated the percentage similarity between bacterial assemblages in the ice core. The closer sections of ice core showed higher resemblance in bacterial DNA structure (Fig. 3), which may be attributed to the fact that the variation of climate and corresponding environment was a gradual process. Parallel climatic conditions support similar vegetation. Therefore, the microorganisms transported to the glacier from distant sources by a defined monsoon were relatively similar at longer timescales (Reference PriscuPriscu and others, 1998). Additionally, with the change of climatic and environmental conditions, the pattern of indigenous bacteriagrowinginthe unique habitat of a glacial surface could shift gradually (Reference Segawa, Miyamoto, Ushida, Agata, Okada and KohshimaSegawa and others, 2005). This might also account for the higher resemblance of bacterial population in closer sections of the ice core. Furthermore, there was less similarity found between the DGGE patterns of the bacterial structure (Fig. 3). The limited percentage similarity (0–15%) between the bacterial community of sea-ice samples was reported as well (Reference Brown and BowmanBrown and Bowman, 2001). However, the similarity between the bacterial community is high in other habitats where the bacterial community is active and shifting (Reference MurrayMurray and others, 1998; Reference Ampe and MiambiAmpe and Miambi, 2000). As discussed below, the low similarity might reflect the episodic deposition of airborne microorganisms and low activity of bacteria in glacial ice.
We performed DGGE to show the bacterial DNA structure in a Puruogangri ice core. Using the Shannon– Weaver index of diversity in combination with cluster analysis of the DGGE banding patterns, we were able to describe the bacterial structure and monitor the responses of bacterial structure to environmental changes. Bacterial species diversity is high in Puruogangri glacier ice, and the dust content is one key determining factor. Furthermore, the Shannon index based on the DGGE pattern of bacterial assemblages in glacial ice is a bio-indicator of climatic changes. More cultivation work and molecular analyses are required for high-altitude glacier habitats in which the microbial diversity still remains largely unstudied.
Acknowledgements
We thank S. Rogers, J. Priscu and E. Mosley-Thompson for reviewing the manuscript and providing constructive comments, and X.J. Zhang for valuable assistance with DGGE analysis. We acknowledge the editorial assistance of P. Mayes. We thank members of the 2000 Puruogangri ice field expedition for their hard work in the field. This study was funded jointly by the National Science Foundation of China (40121101), the Knowledge Innovation Project of the Chinese Academy of Sciences (KZCX1-1-02), National Basic Research program of China (grant No. 2005CB422004) and the Science and Technology Department Foundation of China (2001CCB00300).