The cryosphere is an integral part of the global climate system and one of the major habitable ecosystems of Earth's biosphere. These permanently frozen environments harbor diverse, viable and metabolically active microbial populations that represent almost all the major phylogenetic groups. In this study, we investigated the microbial diversity in the surface snow surrounding the Concordia Research Station on the High Antarctic Plateau through a polyphasic approach, including direct prokaryotic quantification by flow cytometry and catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH), and phylogenetic identification by 16S RNA gene clone library sequencing and 454 16S amplicon pyrosequencing. Although the microbial abundance was low (<103 cells/ml of snowmelt), concordant results were obtained with the different techniques. The microbial community was mainly composed of members of the Alpha-proteobacteria class (e.g. Kiloniellaceae and Rhodobacteraceae), which is one of the most well-represented bacterial groups in marine habitats, Bacteroidetes (e.g. Cryomorphaceae and Flavobacteriaceae) and Cyanobacteria. Based on our results, polar microorganisms could not only be considered as deposited airborne particles, but as an active component of the snowpack ecology of the High Antarctic Plateau.
Citation: Michaud L, Lo Giudice A, Mysara M, Monsieurs P, Raffa C, Leys N, et al. (2014) Snow Surface Microbiome on the High Antarctic Plateau (DOME C). PLoS ONE 9(8): e104505. https://doi.org/10.1371/journal.pone.0104505
Editor: Mitchell Ho, National Cancer Institute, NIH, United States of America
Received: January 28, 2014; Accepted: July 14, 2014; Published: August 7, 2014
Copyright: © 2014 Michaud et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the European Space Agency (ESA-PRODEX) and the Belgian Science Policy (Belspo) through the EXANAM project (C90359). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
The cryosphere, the portion of Earth's biosphere where water is in solid form as snow or ice, includes vast areas of sea ice, freshwater ice, glaciers, ice sheets, snow cover and permafrost . The extremely harsh climatic conditions act as a severe ecological filter for all immigrant and resident organisms and represent the best analogue on Earth for the search of extraterrestrial life , . Cold/frozen environments harbour abundant and diverse microorganisms, which are not merely repositories for wind-transported microorganisms . In particular snow ecological systems, which cover about 35% of the Earth's surface permanently or for varying times during the year , constitute dynamic reservoirs of nutrients and microorganisms .
Studies on the cryosphere revealed the presence of a number of bacteria (mainly belonging to the Bacteroidetes, Actinobacteria, Firmicutes, and Proteobacteria) , which must cope with a combination of severe environmental stressors including (ultra)oligotrophic conditions, high solar UV radiation, freeze/thaw cycles, and limited liquid water . Carpenter et al.  retrieved from the South Pole snow sequences related to Thermus-Deinococcus-like organisms, and reported low rates of bacterial DNA and protein synthesis. Junge et al.  reported that bacterial activity can occur at subzero temperatures. Recently, Lopatina et al.  retrieved from the snow surrounding Antarctic coastal Russian Stations some bacterial genera (Variovorax, Janthinobacterium, Pseudomonas, and Sphingomonas) that may be considered as endogenous Antarctic snow inhabitants.
Several physiological characteristics allow prokaryotes to survive in these harsh environments such as spore formation, pigmentation, increased membrane fluidity and production of enzymes active at low temperature .
Interest in the survival mechanisms of microorganisms in cold environments is increasing because of the possibility to discover new sources of biotechnologically interesting molecules and of the desire to understand the boundaries of life on Earth or elsewhere in our solar system. The evidence of ice on Mars and Europa has made the isolation and study of psychrophiles particularly important for determining the types of organisms that can survive in frozen environments .
The Dome C, situated on top of the Antarctic Plateau (3 233 m above sea level), is one of the coldest and driest places on Earth with temperatures that hardly rise above −25°C in summer (annual average air temperature is −54.5°C). Only a few reports on the diversity of Antarctic Plateau-associated microorganisms exist, although the diversity of Antarctic prokaryotes has been extensively studied , .
The aim of the present work was to explore the prokaryotic diversity and source of microorganisms in the surface snow (corresponding to the annual deposition) of a pristine site on the Dome C, near the Italian-French Research base Concordia. A polyphasic approach, including direct prokaryotic identification by epifluorescence, catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH), flow-cytometry, 16S rRNA gene clone library sequencing and high-throughput 454 16S amplicon pyrosequencing was adopted.
Materials and Methods
Sampling and preliminary treatment of samples
Snow surface sample (top 10 cm corresponding to the annual deposition) was collected in triplicate from a “Clean Area” 2 km from the Research Base “Concordia” (75°06′S–123°20′E) (Fig. 1). No specific permissions were required for sampling activities. Field studies did not involve endangered or protected species.
Snow temperature was −31°C, while the air temperature was −29°C with a wind chill of −45.6°C. The snow was very compact with the presence of extended wavelike sharp and irregular ridges of hard snow (termed “sastrugi”) .
Sampling was performed by using polyethylene boxes pre-treated with 1M hydrogen chloride and hydrogen peroxide, as previously described . Sterile gloves and suit, and an ethanol flame-sterilized shovel were used. Collected samples were allowed to thaw at 4°C for 24–48 h in the laboratory facility of the Base (one hundred litres of packed snow per sample resulted in approximately 20 litres of snowmelt). Autoclave-sterilized Milli-Q water was treated in tandem with the snow samples as a negative-control field blank.
Melted snow samples were processed as follows: a) for total counting and CARD-FISH, 500 ml was fixed with formalin, incubated at 4°C for 1 h and stored at −20°C; b) for flow cytometry analysis, 5 ml was fixed with a freshly prepared solution of 1% paraformaldehyde and 0.05% glutaraldehyde (final concentrations) and then stored at −20°C; c) for DNA extraction, 15 l was filtered through a 0.2-µm-pore-size Sterivex filter unit (Millipore). The filters were stored at −20°C in lysis buffer (50 mM tris, 40 mM EDTA, and 750 mM sucrose).
Samples were packed for shipment and maintained refrigerated at −20°C during the trip back to Europe. Sample temperature was constantly monitored and registered by c\temp™ single-use data logger (Marathon Products Inc.).
Total microbial cell counts by fluorescence microscopy and flow cytometry
For microscopic observation, samples were filtered through 0.2-µm black 25-mm-diameter polycarbonate filters (Isopore) and stained with 4′,6-diamidino-2-phenylindole (DAPI) (1 µg ml−1, final concentration; Sigma) for 10 min. Filters were mounted on slides with non-fluorescent mounting oil (Immersol 518 M). Subsequently, cell enumeration by epifluorescence microscopy was carried out at magnification 1000× by using an Axioplan (Zeiss) fluorescence microscope .
Microbial cells in the snowmelt samples were quantified and characterized by the Flow Cytometer A50-micro (Apogee Flow System, Hertfordshire, England) equipped with a solid state laser set at 20 mV and tuned to an excitation wave length of 488 nm. The volumetric absolute counting was carried out by staining with SYBR Green I (1∶10000 final concentration; Molecular Probes, Invitrogen), while the red autofluorescent cell signals, due to photosynthetic pigments, were used to identify and quantify photoautotrophic microorganisms. Thresholding was carried out using the green or red channel, respectively. The light scattering signals (forward and side light scatter), green fluorescence (535/35 nm), orange fluorescence (590/35 nm) and red fluorescence (>610 nm) were acquired for cell group discrimination. Samples were run at a flow rate of 38.9 µl min−1, in order to keep the number of total events below 1000 event/s. An exclusion gate was applied to avoid visualization of abiotic particles characterized by low green and high red fluorescence. Since cell numbers were expected to be around the detection limit of the instrument, repeated runs were computed together to reach at least 300 positive cells in the counting gate. Data were analyzed using the Apogee Histogram Software (v2.05).
The relative abundance of bacterial groups was examined by performing CARD-FISH reactions with bacteria-specific probes (EUB338 mix), an Archaea-specific probe (ARCH915) and 14 different bacterial group-specific probes (Table 1) –. The probe NON338 was used as a negative control . Samples (100 ml) were filtered at low pressure (200 mbar) onto 0.22 µm white polycarbonate filters (type GTTP; diameter, 47 mm, Millipore). Filters were cut in sections and dipped in 0.2% (w/v) low-gelling point agarose (gel strength 300 g cm−2, Biozym, USA), and then treated as described in Pernthaler et al. ,  and Teira et al. . Filter sections were counterstained with DAPI at a final concentration of 1 µg ml−1 and mounted on a glass slide with Citifluor and Vectashield in a 4∶1 proportion. Slides were stored at −20°C. CARD-FISH filter sections were visualized under a Zeiss epifluorescence microscope. Cell count was done by enumerating at least 500 DAPI-stained cells in 20 fields covering an area of 100 µm×100 µm each and numbers were compared to probe-stained cells in the same field.
Total genomic DNA from melted snow was extracted in triplicate using the phenol-chloroform method according to Zhou et al. , precipitated by adding 0.7 volumes of 100% isopropanol followed by a wash with ice-cold 70% ethanol, and after air-drying resuspended in 50 µl of deionizated sterile water. Quantity and quality of extracted DNA was checked by nanodrop ND-1000 device and the Quant-iT PicoGreen dsDNA reagent and kit (Life Tech, Carlsbad, USA) following the manufacturer's instructions.
454 16S amplicon pyrosequencing
PCR of a bacterial 16S rRNA gene fragment (V1–V3 region, 507 bp) and subsequent tag-encoded pyrosequencing were performed at DNAVision (Charleroi, Belgium). The 16S rRNA genes were amplified using the two universal primers 8F (5′- AGAGTTTGATCCTGGCTCAG -3′) and 518R (5′- ATTACCGCGGCTGCTGG -3′). The forward primer contained the sequence of the Titanium A adaptor (5′-CCATCTCATCCCTGCGTGTCTCCGACTCAG-3′) and a barcode sequence. For each sample, a PCR mix of 100 µl was prepared containing 1×PCR buffer, 2U of KAPA HiFi Hotstart polymerase blend and dNTPs (Kapabiosystems), 300 nM primers (Eurogentec, Liege, Belgium), and 60 ng gDNA. Thermal cycling consisted of initial denaturation at 95°C for 5 min, followed by 25 cycles of denaturation at 98°C for 20 s, annealing at 56°C for 40 s, and extension at 72°C for 20 s, with a final extension of 5 min at 72°C. 3 µl of PCR product were added to a new PCR mix (identical as first round of PCR) for the nested PCR of 15 cycles. Amplicons were visualized on 1% agarose gels using GelGreen Nucleic Acid gel stain in 1× TAE (Biotium) and were cleaned using the Wizard SV Gel and PCR Clean-up System (Promega) according to the manufacturer's instructions. Pyrosequencing was carried out using the forward primer on a 454 Life Sciences Genome Sequencer FLX instrument (Roche) following titanium chemistry.
These pyrosequencing data have been submitted to the Sequence Read Archive of NCBI under the accession number SRX518097, BioProjectID number PRJNA244715.
16S rRNA gene clone library
Extracted DNA from triplicate samples was pooled and firstly used for clone libraries for Bacteria and Archaea. All PCR amplifications were performed with an ABI9600 thermocycler (PE; Applied Biosystems).
For amplification of the bacterial 16S rRNA gene the forward primer 8f (5′-AGAGTTTGATCCTGGCTCAG-3′) and the reverse primer 907r (5′-CCGTCAATTCCTTTRAGTTT-3′) were used. The PCR program and reaction conditions are described in Michaud et al. .
For archaeal 16S rRNA gene amplification two different protocols were tested: a) primers 340f (5′-CCCTAYGGGGYGCASCAG-3′) and 1000r (5′-GAGARGWRGTGCATGGCC-3′)  with GoTaq Green Master Mix (Promega); b) primers Arc20f (5′-GTTTCCGGTTGATCCYGCCRG-3′) and Arc958r (5′-GTTTYCCGGCGTTGAMTCCAATT-3′)  with Takara's Ex Taq DNA Polymerase. No amplification of archaeal 16S rRNA genes was detected.
The bacterial 16S rRNA gene library was constructed using the products of five PCR reactions that were pooled and ligated into the pGEM Easy Vector System (Promega) according to the manufacturer's instructions. The resulting ligation products were used to transform Escherichia coli ElectroMAX DH10B cells (Invitrogen). Inserts were subsequently PCR amplified from lysed white colonies using vector-specific primers, M13f (5′-GTAAAACGACGGCCAGT-3′) and M13r (5′-CAGGAAACAGCTATGACC-3′), under the same previously described PCR conditions . Results were analyzed by agarose gel electrophoresis as described previously . The amplified 16S rRNA gene fragments were directly purified from PCR reaction mixture using the QIAquick PCR Purification Kit (Qiagen) according to the supplier's instructions. Automated sequencing of the 16S rRNA gene clone library was carried out by Sanger sequencing using the dye-terminator method (Macrogen, Korea). Sequence data have been submitted to the GenBank database under accession Nos KJ870937-KJ870986.
Analysis of 16S rRNA gene clone library and 454 amplicon pyrosequencing data
Reads produced by Sanger sequencing and 454 pyrosequencing were trimmed using the mothur software  as described below. 454 generated sequences containing more than two errors in the primer sequence, or more than one error in the barcode sequence were removed. Pyrosequencing reads and Sanger reads were first subjected to basic trimming as recommended  via the removal of reads with read length lower than 200 nts, containing one or more ambiguous bases, or having a homopolymer longer than 8 (only applied for pyrosequencing data). Next, a quality trimming approach was applied using a sliding window of 100 nt length with a cut-off on the Phred score of 30 and 20 for pyrosequencing and Sanger reads, respectively. For Sanger, we used the TraceTuner software (3.0.6; http://sourceforge.net/projects/tracetuner/) in order to convert the quality files in Phred score. All reads were aligned to the Silva reference database (http://www.mothur.org/wiki/Silva_reference_alignment), thereby removing reads starting extremely late, or ending extremely early in the alignment (belong in the 2% outliers of the start and end) using the mothur filter and screen commands. Pyrosequencing reads were denoised using NoDe (Mysara et al., unpublished data; Noise Detector (NoDe) is an artificial intelligence-based classifier that first identifies positions likely to be erroneous and subsequently clusters those error-prone reads with correct reads resulting in error-free consensus reads.), while the Sanger sequencing data were denoised using the SLP algorithm  as implemented in mothur, using a maximum of 2% differences. Chimera removal and operational taxonomic unit (OTU) clustering of the pyrosequencing reads was performed using UPARSE . For the Sanger sequencing data we used the UCHIME chimera detection tool  together with nearest neighbour clustering algorithm implemented in mothur (with 0.03 distance). Phylotyping was performed using the RDP classifier  retrained with the greengenes dataset (using 80% confidence threshold as cut off), as recommended in .
Phylogenetic tree construction
For each phylotype the most representative read was selected by calculating the pairwise distance between all reads in a single Phylotype and selecting the read with the smallest maximum distance to all other reads (mothur Get.Oturep command), which led to 23 and 178 reads for the clone library and pyrosequencing data, respectively. For the pyrosequencing reads, only the 15 most abundant phylotypes were selected in order to improve readability of the tree. The closest relatives of these phylotypes were obtained using the RDP classifier  with the RDP cultured and the greengenes dataset, leading to 33 and 34 unique hits, respectively. The full set of sequences was aligned against the SILVA database and the resulting alignments were cleaned (using the align.ses and filter.seqs commands as available in mothur ). The distance was calculated and parsed for tree creation using mothur (dist.seqs and clearcut commands) . Visualization of the phylogenetic tree was done using FigTree [http://tree.bio.ed.ac.uk/software/figtree/].
Scanning electron microscopy
For scanning electron microscopy (SEM) samples were gently thawed at 4°C and filtered on a 0.02 µm (25 mm diameters) inorganic membrane (Anodisc; Whatman). Filters were fixed with glutaraldehyde (4%) in 0.1% cacodylate buffer at pH 7.4 for 3 h and then dehydrated with 25, 50, 70, 85, and 95% ethanol for 5 min. Subsequent process manipulation and photography were carried out using a SEM JEOL JSM-6500F (JEOL-USA Inc., Peabody, MA USA).
Prokaryotic and pico-eukaryotic abundances
Direct counts by epifluorescence microscopy revealed population densities of 4.3±2.2×102 cells ml−1 of snowmelt (mean ± standard deviation). In turn, 1.7±0.6×102 autofluorescent cells ml−1 of snowmelt were detected.
Flow cytometry showed a comparable number of non-autofluorescent prokaryotes in the snowmelt samples (2.5±0.5×102 cells ml−1) (Fig. 2a). An abundance of 0.9±0.3×102 cyanobacterial cells ml−1 (26.5% of the total counted prokaryotes) was determined based on the autofluorescent signals in an orange-vs-red cytograms (Fig. 2b). The autofluorescent signals were coherent with those of cells containing phycocyanin and phycoerythrin. A further cytometric population of photoautotrophs showed red fluorescence and forward scatter signals similar to those of picoeukaryotes, but the exact number was not properly determinable by flow cytometry. However, a qualitative microscopic control confirmed the presence of a few picoeukaryotic cells in the snowmelt sample (Fig. 2c).
Total prokaryotes (a), Cyanobacteria and picoeukaryotes (b) retrieved in the snowmelt samples by SYBR Green I staining and autofluorescence signals, respectively. The microscopic control confirmed their presence according to cell fluorescence and size (c).
Around 81.9±7.6% of DAPI-stained cells could be visualized by CARD-FISH with HRP-labeled oligonucleotide probe EUB338. By the applied protocol and stringency conditions, only very few cells appeared as positively hybridized by the probe ARCH915. The set of probes used to detect the major divisions within the Bacteria domain associated 98% of the EUB338 hybridized cells to known bacterial groups (Fig. 3). The negative control probe detected 0 to 2% of the DAPI-stained cells.
Only phyla with abundance above 1% are shown.
Among the Proteobacteria (32.5±10.5% of DAPI- stained cells), Alpha-, Beta- and Gammaproteobacteria accounted for 11.8±3.6, 5.8±2.0, and 12.8±2.8%, respectively. Only 1.2±0.9 and 0.9±1% of EUB338 stained cells were detected by probes for the Delta- and Epsilonproteobacteria. The LGC, HGC and CF319 stained cells (Firmicutes, Actinobacteria and Cytophaga-Flavobacterium-Bacteroidetes group, respectively) accounted for 7.5±2.9, 1.8±0.9 and 20.3±4.2%, respectively. Autotrophic pico-cyanobacteria were particularly abundant and represented 35.9±10.0% of DAPI-stained cells (Fig. 3).
454 16S amplicon pyrosequencing
Extracted DNA from snow sample amounted to 8.9 ng µl−1. From the initial data set containing 16914 reads (mean read length 453 bp), 16674 reads (mean read length 419 bp) were retained after filtering based on the quality of the barcode and primer sequence (barcode and primer sequence were removed in this step), 14657 reads (mean read length 404 bp) after basic and quality trimming, 14142 reads (mean read length 408 bp) after screening based on the alignment of the reads to the SILVA database. After denoising, those 14142 reads can be represented by 3772 unique reads. Applying UPARSE (combining chimera detection and OTU clustering) resulted in a total number of 820 OTUs representing 14043 reads with an average length of 420 bp (Table S1 in File S1). Next, the RDP classifier assigned those 820 OTUs to 303 different phylotypes (Table S2 in File S1). The obtained data covered a broad spectrum of known bacterial phyla, with the predominance of Proteobacteria (40.4%; 140 phylotypes) and Bacteroidetes (23.9%; 47 phylotypes), followed by Cyanobacteria (15.6%; 14 phylotypes), Actinobacteria (2.3%; 14 phylotypes) and Verrucomicrobia (1.5%; 15 phylotypes) (Fig. 3). The Fusobacteria, Planctomycetes, Lentisphaerae, Firmicutes, Chlorobi, Tenericutes, Acidobacteria, Spirochaetes and Chloroflexi were less abundant and were designated as “others” in Fig. 3. Most of the sequences that remained unclassified at the phylum level (11.1%) belonged to Phylotype_P001 (Table S3 in File S1).
Among the Proteobacteria, members of the Alpha- (30.0%; 52 phylotypes), Beta- (2.8%; 15 phylotypes), Gamma- (3.5%; 40 phylotypes), Delta- (2.1%; 24 phylotypes), and Epsilonproteobacteria (1.9%; 7 phylotypes) were detected (Fig. 4, Tables S3 and S4 in File S1). The classified Alphaproteobacteria were mainly represented by Kiloniellaceae (7.5%; 3 phylotypes), Rhodobacteraceae (5.7%; 10 phylotypes) and Rhodospirillaceae (1.7%; 6 phylotypes). Almost half of the Alphaproteobacteria (40.7%) remained unclassified and most belonged to Phylotype_P002 (Table S3 in File S1). Also most of the Gammaproteobacteria remained unclassified, although 22.3% were SUP05-related (sulfur oxidizers). The classified Betaproteobacteria were mostly assigned to Methylophilales (Methylophilaceae), Burkholderiales (Comamonadaceae) and Rhodocyclales (Rhodocyclaceae) (Tables S3 and S4 in File S1). The order Sva085, Desulfobacterales (Desulfobacteraceae and Desulfobulbaceae families) were most abundant among Deltaproteobacteria. These families include sulfate-reducing organisms such as those belonging to the Desulfotalea genus (Phylotype_P088). Most Epsilonproteobacteria (96.6%) were affiliated to the Campylobacterales with the predominance of members of the genera Arcobacter and Sulfurimonas (Tables S3 and S4 in File S1).
Only classes with abundance above 1% are shown.
Classified Bacteroidetes were predominantly affiliated to the Flavobacteriaceae. Approximately 70% of the Cyanobacteria/Chloroplast sequences were related to algal chloroplasts, including genera from Chlorophyta, Cryptophyta, Stramenopiles, Haptophyceae and Streptophyta. In addition, 27.2% of the Cyanobacteria/Chloroplast sequences were related to Synechococcus. Among Actinobacteria, Acidimicrobiales (OCS155 and C111) and Microbacteriaceae were most abundant. Among Verrucomicrobia, Opitutae (unclassified and Puniceicoccaceae) and Verrucomicrobiaceae were most abundant (Tables S3 and S4 in File S1).
16S rRNA gene clone library
Amplification of bacterial 16S rRNA genes was positive. However, no amplification of archaeal 16S rRNA genes was detected, thus supporting the outcome of the CARD-FISH analysis. The construction and screening of a 16S rRNA gene clone library (95 clones) yielded an initial data set containing 72 reads (mean read length 1078 bp). From this set, 63 reads (mean read length 924 bp) were retained after basic and quality trimming, and screening based on the alignment of the reads to the SILVA database (Table S1 in File S1). After denoising, those 63 reads can be represented by 59 unique reads. Applying chimera detection and OTU clustering resulted in a total number of 50 OTUs representing 63 reads with an average length of 923 bp (Table S1 in File S1). Next, the RDP classifier assigned those 50 OTUs to 25 different phylotypes (Tables S5 and S7 in File S1). The obtained data covered four bacterial phyla, with the predominance of Proteobacteria (63.5%; 15 phylotypes) and Cyanobacteria/Chloroplast (22.2%; 4 phylotypes), followed by Bacteroidetes (11.1%; 4 phylotypes) and Actinobacteria (1.6%; 1 phylotype) (Fig. 3). The classified Alphaproteobacteria were mainly represented by Pelagibacteraceae (25.4%), and Rhodobacteraceae (6.4%). The classified Gammaproteobacteria included Aeromonas. The Bacteroidetes were affiliated to the Cryomorphaceae and Flavobacteriaceae families. Fifty per cent of the Cyanobacteria/chloroplast sequences were related to algal chloroplasts (Table S7 in File S1).
Comparison of different approaches
Overall a reasonable agreement between the composition of bacterial phyla were found with the three techniques (Fig. 3). Also at class level, a good agreement was found with both sequencing techniques (Figs. 4 and 5). However, some differences were found between the contribution of specific bacterial groups to the total community for instance for Flavobacteria and Cyanobacteria. Nevertheless, only the increased contribution of Cyanobacteria in the clone library was found to be statistically deviating from the results based on the pyrosequencing reads (p-value <0.05 using the hypergeometric distribution). As the amount of extracted DNA was limited, our outcome could be biased by contaminating bacterial DNA in reagents , . However, the comparable outcome of both sequencing techniques, which used different types of polymerase and reagents, indicated that the overall community composition is unlikely to be attributed to contamination.
This study represents the first report on the microbiome of the surface snow surrounding the Research Base “Concordia” on the high Antarctic Plateau (3233 m above the sea level at Dome C on the Eastern Antarctic plateau). This unique environment is used (next to glaciology, atmospheric sciences, astronomy, etc.) as a test bed for long-duration spaceflights to study the physiologic and psychological adaptation to isolated environments . In addition, Antarctica is seen as an analogue of extraterrestrial environments .
Bacterial densities estimated by direct epifluorescence microscopy and flow cytometry are in line with previous reports for South Pole snow (200 to 5000 cells ml−1) , but were two orders of magnitude lower than in snow at Ny-Alesund (Spitzberg) (6×104 cells ml−1) , on the Tibetan Plateau (0.68 to 720×103 cells ml−1) , or in Alpine sites (1.1×104 cells ml−1) . Moreover, about 30% of the total counted prokaryotes showed the typical autofluorescence signature of phycocyanin- and phycoerythrin-containing Cyanobacteria, as estimated by either flow cytometry (26.5%) or epifluorescence microscopy (35.9%). This finding was supported by the identification of Cyanobacteria by both sequencing techniques.
Overall, the bacterial community composition in surface snow surrounding the Concordia Base at Dome C was similar to that reported in glaciers, snow, lake ice, sea ice, and atmospheric clouds . Although the site is remote from the sea, a strong influence from the marine environment seems to exist with Proteobacteria being the dominant members, followed by Bacteroidetes and Cyanobacteria . Allochthonous input of marine microorganisms and other particulate organic matter are introduced by atmospheric circulation and deposition , , , . Free-living microorganisms have high dispersal rates and can be transported over very long distances because of their microscopic size and resistance to environmental stressors , . However, we cannot a priori exclude the possibility that a minor part of microorganisms could come from human activity even when the chosen sampling site was two kilometers from Concordia Station in the restricted access “Clean Area” zone.
Among the Alphaproteobacteria, Kiloniellaceae and Rhodobacteraceae were the most abundant groups. Many were phylogenetically related to sequences retrieved from marine origin (Fig. 6). Such a family exhibits a variety of growth modes like photo- and chemo-heterotrophic growth under light anoxic conditions and dark micro-oxic to oxic conditions, respectively, as well as photolithoautotrophic growth. Different members of the Rhodobacteraceae (Amaricoccus, Anaerospora, Loktanella, Octadecabacter, Paracoccus, Phaeobacter and Thalassobius) were detected with 454 16S amplicon pyrosequencing. Nevertheless, almost all of the Kiloniellaceae and 75.5% of the Rhodobacteraceae reads remained unaffiliated.
The closest relatives to the phylotypes from the 16S rRNA gene clone library (green; Phylotype S) and 454 16S amplicon pyrosequencing (15 most abundant) (red; Phylotype P) were retrieved from the RDP cultured and the greengenes dataset and compiled to calculate an unrooted phylogenetic tree via mothur. All relatives isolated from marine origin are shown in blue.
The Gamma-, Beta-, Epsilon- and Deltaproteobacteria were detected at lower levels (454 16S amplicon pyrosequencing). Among the Epsilonproteobacteria, the most represented genera were Arcobacter and Sulfurimonas. Arcobacter is a unique genus harboring species that can grow microaerobically, occupying various ecological niches, e.g. in association with animals (and humans) and other species that are free-living and found in seawater . Among the classified Gamma- and Deltaproteobacteria, members related to the sulfur and nitrogen cycles were retrieved.
Members of the Bacteroidetes, predominantly Flavobacteriaceae, represented a significant fraction of the bacterial community in the surface snow at Dome C. Members of Flavobacteriaceae, originate from marine environments and are important for the mineralization of organic matter in marine ecosystems  (Fig. 6).
Another main group identified was the Cyanobacteria class (dominated by the Synechococcus genus). Marine synechococci are ubiquitously distributed throughout the oceans and can be found in any ecosystem type up to the polar circle  including extreme regions like the Mount Everest (in the snow at 6 500 m altitude)  and in extensive snow cover on the Byers Peninsula (Livingston Island, maritime Antarctica) . The Cyanobacteria are organisms that occur in most sun-exposed ecosystems on Earth and are frequently found in ice as well as cold desert habitats . They are relatively rare in non-stable cryo-ecosystems as melting snow, but have been recorded in glacier snow from several areas of the world , .
The phylogenetic affiliation also revealed the presence of eukaryotic (chloroplasts) taxa. The closest relatives of many phylotypes were related to marine environments. In particular, Chlorophyta and Cryptophyta were identified in both the cloning and the pyrosequencing dataset. This data is supported by the cytometric and microscopic observation of picoeukaryotes (Fig. 2).
The more than 300 phylotypes identified in this work represent a significant diversity taking into account the environmental conditions and remoteness of the site. The dataset covers the majority (98.2% coverage) of the diversity and is in line with previous reports, which include heterotrophic bacteria, cyanobacteria and eukaryotes, with many of them that are related to known psychrophilic and psychrotolerant species .
Although our results provide no direct evidence for possible activity of the depicted bacterial community, some may be indicative for activity. The community could be organized in a patchy biofilm structure suggested by the presence of exopolysaccharide-like debris on the DAPI-stained filters (not shown) and by scanning electron microscopy (Fig. 7). Furthermore, exopolysaccharides are thought to be a primary protection mechanism by reducing water loss and by restricting ice crystal formation to sites outside the cells . In addition, Cyanobacteria, detected as auto-fluorescent cells via flow cytometry, as well as bacteria involved in carbon, nitrogen and sulfur cycles, were identified.
Pictures are at increasing magnification from a) to d).
It has been reported that microbial abundance in the snowpack is usually positively correlated with Ca2+ concentrations, serving as a proxy for dust . Moreover, microorganisms in dry polar snow were shown to be involved in active exchange of reactive nitrogen species with the atmosphere and were thus contributing to biogeochemical cycling at low temperature . Despite the low temperatures and limited liquid water supply, microbes inhabiting the snow cover seemed to constitute a dynamic community, as it was recently suggested by the observation of a link between microbial community structure and snow chemistry . Survival of prokaryotes and unicellular eukaryotes has been demonstrated down to −18°C in unfrozen hypersaline polar lakes, as well as metabolism of bacteria from snow at −17°C and from permafrost at −20°C , . However, the lower temperature limit for cell growth and division remains unclear. At least for nitrification, some reports exist on activities at temperatures close to the Concordia surface snow at sampling time. Miteva and colleagues  observed experimentally nitrification by Nitrosomonas cryotolerans in solutions frozen at −32°C. Price and Sowers  reported evidence for metabolically active bacteria and fungi at temperatures down to −40°C. Finally, Sowers , who studied the N2O in Antarctic ice cores, suggested that the recorded artifacts correlated with high concentrations of dust and high bacterial counts could be explained by the presence of metabolically active bacteria, presumably brought to the ice caps with the dust, that had produced sufficient N2O over the years to alter the composition of the archived air.
This study represents the first investigation focused on the microbiome of the surface snow surrounding the Research Base “Concordia” on the high Antarctic Plateau. In such remote and extreme environment, a strong influence from the coastal areas seems to exist with marine Proteobacteria (i.e. Kiloniellaceae and Rhodobacteraceae) representing the dominant members within the microbial community, followed by Bacteroidetes and Cyanobacteria. Thus, our findings support the evidence that polar microorganisms should be considered not only as deposited airborne particles, but also as an active component of the snowpack ecology.
Dedicated to the memory of Luigi Michaud, passionate researcher, who died in Antarctica during sampling activities.
Authors wish to thank all of the staff at “Concordia” Station, for the logistic help and support, the “Institute Polaire Paul Emile Victor” (IPEV), the Italian “Programma Nazionale di Richerche in Antartide” (PNRA) and the MNA (Museo Nazionale dell'Antartide), which made the expedition possible.
Conceived and designed the experiments: LM. Performed the experiments: LM PM CR SA RVH. Analyzed the data: LM ALG MM PM NL SA RVH. Wrote the paper: LM ALG MM PM SA RVH.
- 1. Miteva V (2008) Bacteria in snow and glacier ice. In: Margesin R, Schinner F, Marx JC, Gerday C, editors. Psychrophiles: from biodiversity to biotechnology. Springer: Berlin. pp. 31–50.
- 2. Lugg D, Shepanek M (1999) Space analogue studies in Antarctica. Acta Astronaut 44: 693–699.
- 3. Pyne SJ (2007) The extraterrestrial Earth: Antarctica as analogue for space exploration. Space Policy 23: 147–149.
- 4. Harding T, Jungblut AD, Lovejoy C, Vincent WF (2011) Microbes in high Arctic snow and implications for the cold biosphere. Appl Environ Microbiol 77: 3234–3243.
- 5. Liu Y, Yao T, Jiao N, Kang S, Xu B, et al. (2009) Bacterial diversity in the snow over Tibetan Plateau Glaciers. Extremophiles 13: 411–423.
- 6. Yan F, Hou S, Chen T, Ma X, Zhang S (2012) Culturable bacteria isolated from snow cores along the 1300 km traverse from Zhongshan Station to Dome A, East Antarctica. Extremophiles 16: 345–354.
- 7. Carpenter EJ, Lin S, Capone DG (2000) Bacterial activity in South Pole snow. Appl Environ Microbiol 66: 4514–4517.
- 8. Junge K, Eicken H, Deming JW (2004) Bacterial activity at −2 to −20 degrees C in Arctic wintertime sea ice. Appl Environ Microbiol 70: 550–557.
- 9. Lopatina A, Krylenkov V, Severinov K (2013) Activity and bacterial diversity of snow around Russian Antarctic stations. Res Microbiol 164: 949–958.
- 10. Amato P, Hennebelle R, Magand O, Sancelme M, Delort AM, et al. (2007) Bacterial characterization of the snow cover at Spitzberg, Svalbard. FEMS Microbiol Ecol 59: 255–264.
- 11. Miteva VI, Sheridan PP, Brenchley JE (2004) Phylogenetic and physiological diversity of microorganisms isolated from a deep Greenland glacier ice core. Appl Environ Microbiol 70: 202–213.
- 12. Tindall BJ (2004) Prokaryotic diversity in the Antarctic: the tip of the iceberg. Microb Ecol 47: 271–283.
- 13. Michaud L, Caruso C, Mangano S, Interdonato F, Bruni V, et al. (2012) Predominance of Flavobacterium, Pseudomonas, and Polaromonas within the prokaryotic community of freshwater shallow lakes in the northern Victoria Land, East Antarctica. FEMS Microbiol Ecol 82: 391–404.
- 14. Van Houdt R, Deghorain M, Vermeersch M, Proovost A, Lo Giudice A, et al. (2013) Characterization of culturable Paenibacillus spp. from the snow surface on the high Antarctic Plateau (DOME C) and their dissemination in the Concordia research station. Extremophiles 17: 565–573.
- 15. Porter KG, Feig YS (1980) The use of DAPI for identifying and counting aquatic microflora. Limnol Oceanogr 25: 943–948.
- 16. Amann RI, Binder BJ, Olson RJ, Chisholm SW, Devereux R, et al. (1990) Combination of 16S rRNA- targeted oligonucleotide probes with flow cytometry for analyzing mixed microbial populations. Appl Environ Microbiol 56: 1919–1925.
- 17. Daims H, Brühl A, Amann R, Schleifer K-H, Wagner M (1999) The domain-specific probe EUB338 is insufficient for the detection of all Bacteria: development and evaluation of a more comprehensive probe set. Syst Appl Microbiol 22: 434–444.
- 18. Neef A (1997) Application of in situ identification of bacteria to population analysis in complex microbial communities. PhD Dissertation, Munich.
- 19. Manz W, Amann R, Ludwig W, Wagner M, Schaefer KH (1992) Phylogenetic oligodeoxynucleotide probes for the major subclasses of proteobacteria: problems and solutions. Syst Appl Microbiol 15: 593–600.
- 20. Loy A, Lehner A, Lee N, Adamczyk J, Meier H, et al. (2002) Oligonucleotide microarray for 16S rRNA gene based detection of all recognized lineages of sulfate reducing prokaryotes in the environment. Appl Environ Microbiol 68: 5064–5081.
- 21. Loy A, Maixner F, Wagner M, Horn M (2007) probeBase: an online resource for rRNA-targeted oligonucleotide probes. Nucleic Acids Res 35: 800–804.
- 22. Manz W, Amann R, Ludwig W, Wagner M, Schaefer KH (1996) Application of a suite of 16S rRNA specific oligonucleotide probes designed to investigate bacteria of the phylum Cytophaga-Flavobacter-Bacteroides in the natural environment. Microbiology 142: 1097–1106.
- 23. Roller C, Wagner M, Amann R, Ludwig W, Schleifer KH (1994) In situ probing of Gram-positive bacteria with high DNA G+C content using 23S rRNA-targeted oligonucleotides. Microbiology 140: 2849–2858.
- 24. Meier H, Amann R, Ludwig W, Schleifer KH (1999) Specific oligonucleotide probes for in situ detection of a major group of Gram-positive bacteria with low DNA G+C content. Syst Appl Microbiol 22: 186–196.
- 25. Schonhuber W, Zarda B, Eix S, Rippka R, Herdman M, et al. (1999) In situ identification of cyanobacteria with horseradish peroxidase-labeled, rRNA-targeted oligonucleotide probes. Appl Environ Microbiol 65: 1259–1267.
- 26. Stahl DA, Amann R (1991) Development and application of nucleic acid probes. In: Stackebrand E, Goodfellow M, editors. Nucleic acid techniques in bacterial systematics. John Wiley and Sons Ltd.: Chichester, UK. Pp. 205–248.
- 27. Wallner G, Amann R, Beisker W (1993) Optimizing fluorescent in situ hybridization of suspended cells with rRNA-targeted oligonucleotide probes for the flow cytometric identification of microorganisms. Cytometry 14: 136–143.
- 28. Pernthaler A, Pernthaler J, Amann R (2002) Fluorescence in situ hybridization and catalyzed reporter deposition for the identification of marine bacteria. Appl Environ Microbiol 68: 3094–3101.
- 29. Pernthaler A, Preston CM, Pernthaler J, DeLong EF, Amann R (2002) Comparison of fluorescently labelled oligonucleotide and polynucleotide probes for the detection of pelagic marine bacteria and Archaea. Appl Environ Microbiol 68: 661–667.
- 30. Teira E, Reinthaler T, Pernthaler A, Pernthaler J, Herndl GJ (2004) Combining catalyzed reporter deposition-fluorescence in situ hybridization and microautoradiography to detect substrate utilization by bacteria and Archaea in the deep ocean. Appl Environ Microbiol 70: 4411–4414.
- 31. Zhou J, Bruns MA, Tiedje JM (1996) DNA recovery from soils of diverse composition. Appl Environ Microbiol 62: 316–322.
- 32. Gantner S, Andersson AF, Alonso-Sáez L, Bertilsson S (2011) Novel primers for 16S rRNA-based archaeal community analyses in environmental samples. J Microbiol Met 84: 12–18.
- 33. La Cono V, La Spada G, Arcadi E, Placenti F, Smedile F, et al. (2013) Partaking of Archaea to biogeochemical cycling in oxygen-deficient zones of meromictic saline Lake Faro (Messina, Italy). Environ Microbiol 15: 1717–1733.
- 34. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, et al. (2009) Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75: 7537–7541.
- 35. Schloss PD, Gevers D, Westcott SL (2011) Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS One 6(12): e27310.
- 36. Huse SM, Welch DM, Morrison HG, Sogin ML (2010) Ironing out the wrinkles in the rare biosphere through improved OTU clustering. Environ Microbiol 12: 1889–1898.
- 37. Edgar RC (2013) UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods 10: 996–998.
- 38. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R (2011) UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27: 2194–2200.
- 39. Wang Q, Garrity GM, Tiedje JM, Cole JR (2007) Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol 16: 5261–5267.
- 40. Böttger EC (1990) Frequent contamination of Taq polymerase with DNA. Clin Chem 36: 1258–1259.
- 41. Hughes MS, Beck LA, Skuce RA (1994) Identification and elimination of DNA sequences in Taq DNA polymerase. J Clin Microbiol 32: 2007–2008.
- 42. Sattler B, Puxbaum H, Psenner R (2001) Bacterial growth in supercooled cloud droplets. Geophys Res Lett 28: 239–242.
- 43. Giovannoni SJ, Stingl U (2005) Molecular diversity and ecology of microbial plankton. Nature 437: 343–348.
- 44. Cowan DA, Chown SL, Convey P, Tuffin M, Hughes K, et al. (2011) Non-indigenous microorganisms in the Antarctic: assessing the risks. Trends Microbiol 19: 540–548.
- 45. Vincent WF (2000) Evolutionary origins of Antarctic microbiota: invasion, selection and endemism. Antarct Sci 12: 374–385.
- 46. Aller JY, Kuznetsova MR, Jahns CJ, Kemp PF (2005) The sea surface microlayer as a source of viral and bacterial enrichment in marine aerosols. J Aerosol Sci 6: 801–812.
- 47. Donachie SP, Bowman JP, On SLW, Alam M (2005) Arcobacter halophilus sp. nov., the first obligate halophile in the genus Arcobacter. Int J Syst Evol Microbiol 55: 1271–1277.
- 48. Cottrell MT, Kirchman DL (2000) Natural assemblages of marine proteobacteria and members of the Cytophaga-Flavobacter cluster consuming low- and high-molecular-weight dissolved organic matter. Appl Environ Microbiol 66: 1692–1697.
- 49. Scanlan DJ, Ostrowski M, Mazard S, Dufresne A, Garczarek L, et al. (2009) Ecological genomics of marine picocyanobacteria. Microbiol Mol Biol Rev 73: 249–299.
- 50. Liu Y, Yao T, Kang S, Jiao N, Zeng YH, et al. (2006) Seasonal variation of snow microbial community structure in the East Rogbuk glacier, Mt. Everest. Chin Sci Bull 51: 1476–1486.
- 51. Velázquez D, Rochera C, Camacho A, Quesada A (2011) Temperature effects on Antarctic non-marine phototrophic communities. Polar Biol 34: 1045–1055.
- 52. Vincent WF (2007) Cold tolerance in cyanobacteria and life in the cryosphere. In: Seckbach J, editor. Algae and cyanobacteria in extreme environments. Springer: Heidelberg. Pp. 289–304.
- 53. Quesada A, Vincent V (2012) Cyanobacteria in the cryosphere: snow, ice and extreme cold. In: Whitton BA, editor. Ecology of Cyanobacteria II: their diversity in space and time. Springer. pp. 387–399.
- 54. Margesin R, Miteva V (2011) Diversity and ecology of psychrophilic microorganisms. Res Microbiol 162: 346–361.
- 55. Amoroso A, Domine F, Esposito G, Morin S, Savarino J, et al. (2010) Microorganisms in dry polar snow are involved in the exchanges of reactive nitrogen species with the atmosphere. Environ Sci Technol 44: 714–719.
- 56. Larose C, Dommergue A, Vogel TM (2013) Microbial nitrogen cycling in Arctic snowpacks. Environ Res Lett 8: 035004.
- 57. Clarke A, Morris GJ, Fonseca F, Murray BJ, Acton E, et al. (2013) A low temperature limit for life on earth. Plos One 8: 1–11.
- 58. Miteva V, Sowers T, Brenchley J (2007) Production of N2O by ammonia oxidizing bacteria at subfreezing temperatures as a model for assessing the N2O anomalies in the Vostok ice core. Geomicrobiol J 24: 451–459.
- 59. Price PB, Sowers T (2004) Temperature dependence of metabolic rates for microbial growth, maintenance, and survival. Proc Natl Acad Sci USA 101: 4631–4636.
- 60. Sowers TN (2001) N2O record spanning the penultimate deglaciation from the Vostok ice core. J Geophys Res -Atmos 106: 31903–31914.