Abstract
It was once a long-held view that the Antarctic was a pristine environment with low biomass, low biodiversity and low rates of microbial activity. However, as the intensity of scientific investigation has increased, so these views have started to change. In particular, the role and impact of human activity toward indigenous microbial communities has started to come under more intense scrutiny. During the Subglacial Lake Ellsworth exploration campaign in December 2012, a microbiological survey was conducted to determine the extent and likelihood of exogenous input into the subglacial lake system during the hot-water drilling process. Snow was collected from the surface to represent that used for melt water production for hot-water drilling. The results of this study showed that snow used to provide melt water differed in its microbiological composition from that of the surrounding area and raised the question of how the biogeography of snow-borne microorganisms might influence the potential outcome of scientific analyses. In this study, we investigated the biogeography of microorganisms in snow around a series of Antarctic logistic hubs, where human activity was clearly apparent, and from which scientific investigations have been undertaken. A change in microbial community structure with geographical location was apparent and, notably, a decrease in alpha diversity at more remote southern latitudes. Soil-related microorganisms dominated microbial assemblages suggesting terrestrial input, most likely from long-range aeolian transport into continental Antarctica. We also observed that relic DNA was not a major issue when assessing snow samples. Overall, our observations might have profound implications for future scientific activities in Antarctica, such as the need to establish “no-go” protected areas, the need for better characterization of field sites and improved protocols for sterilization and verification of ice drilling equipment.
Introduction
The cryosphere covers about 20% of Earth’s surface () and seasonally snow-covered areas cover up to 35% of the terrestrial surface (). Snow represents a climatically sensitive transient ecosystem, linking the atmosphere to the soil or ice below (); it is physically, chemically and biologically dynamic, has a strong impact on the hydrological cycle, and interacts with a large number of other ecosystems worldwide (Vavrus, 2007; ). Pristine snow is primarily seeded from the atmosphere and unique microbial communities thus develop within the snow over time (; ), which then actively participate in nutrient cycling (). Within the snow, microorganisms encounter very low temperatures, variable UV radiation (at the surface), limited water availability and oligotrophic nutrient conditions (). Yet, microorganisms living in such habitats have been identified worldwide, including within glacial ice (), cryoconite holes (), sea ice (Staley and Gosink, 1999) and snow (). Antarctic snow is an important habitat as over 99% of the continent is permanently ice covered () and microbial life dominates any other form of terrestrial life there (Vincent, 1988).
Studies of snow-borne microorganisms have mainly focused on Arctic regions, with fewer studies on Alpine and Antarctic snow communities. Arctic snow has been shown to harbor primarily Proteobacteria (Alpha- and Beta-), along with Bacteroidetes, Firmicutes, and Cyanobacteria (; ; ; ; ). The two large high throughput sequencing studies on Alpine snow microbes show a dominance of Alphaproteobacteria and Betaproteobacteria, with variable relative abundances of Firmicutes, Bacteroidetes, Acidobacteria, Actinobacteria and Cyanobacteria (Wunderlin et al., 2016; ). While more recent studies have investigated the physico-chemical properties of Antarctic snow (), patterns of accumulation and ablation () and optical properties (), only a limited number have investigated microbial diversity and distribution, and these report high abundances of Proteobacteria (Alpha-, Beta-) and variable levels of Firmicutes, Bacteroidetes and Cyanobacteria (; ; ). These studies are restricted to a limited number of specific locations leaving the vast majority of the continent unexplored. Hence, broad scale biogeographic patterns of Antarctic snow microorganisms are unknown. In addition, knowledge of colonization pathways, and whether Antarctic ecosystems are seeded primarily by local dispersal processes or long-range transport, is still open to question (Pearce et al., 2016).
There have been a number of Antarctic subglacial lake access projects in recent years such as those aiming to explore Lake Vostok in East Antarctica (), Lake Whillans at the edge of the West Antarctic Ice Sheet (Tulaczyk et al., 2014; Siegert et al., 2016), Lake Ellsworth in the center of West Antarctica (Siegert et al., 2014) and the recent ‘subglacial Antarctic lakes scientific access’ (SALSA) project focused on Lake Mercer and Lake Engelhardt on the Whillans Ice Plain. While these programs follow strict protocols to avoid microbial contamination of the subglacial ecosystems they intend to access (Priscu et al., 2013; ), they have exposed an issue that little attention is being given to indigenous microbial communities at the snow surface and what impact they might potentially have on the results obtained. For example, it is not yet known whether the snow surface communities enter the subglacial system directly via migration through the ice or indeed anthropogenically via ineffective sterilization of hot water drill fluid. Elsewhere, studies have shown the presence of human-associated bacteria in the vicinity of research stations () and field camps (Sjoling and Cowan, 2000), suggesting a potential impact of human activity, including in deep-field locations ().
Here, we investigate the spatial distribution of microbial communities in Antarctic surface snow in areas close to logistic hubs and areas of strong scientific research interest, from Signy Island to the Ellsworth Mountains. We also characterized the proportion of the bacterial community that was potentially viable using a propidium monoazide (PMA) assay, thus allowing us to assess the impact of relic DNA on Antarctic snow diversity studies. From these analyses, we identified patterns of microbial distribution and spatial variability within the Antarctic surface snowpack, and assessed bacterial viability. Our investigation underlines the importance of a full understanding of the microbiology of the field site prior to invasive field operations.
Materials and Methods
Sample Collection
Snow samples were collected between the South Orkney Islands and the Ellsworth Mountains between December 2012 and January 2014 (Figure 1 and Supplementary Table 1). On Signy Island, two sites were sampled; Gourlay Snowfield and Tuva Glacier. Signy Island is an important bird area and hosts a research station, maintained by the British Antarctic Survey (BAS). While the sampled areas are not directly impacted by the research station, there is human presence on the island and ongoing research activity in the vicinity. Livingston Island samples are described in the study by , who investigated the relationship between microbial communities and optical and physico-chemical properties of Antarctic snow. Livingston Island hosts five scientific bases, and occasional field camps. SkyBlu samples came from the vicinity of the BAS blue ice runway (a logistics hub for deep field operations). Pine Island Bay area is of research interest as the Pine Island Glacier and Thwaites Glacier are fast retreating and their collapse could lead to major sea level rise (Rignot et al., 2014). Ellsworth off-site samples were collected approximately one kilometer upwind of the Lake Ellsworth deep-field camp (Pearce D.A. et al., 2016). Camp samples were collected around the kitchen, generator and drilling areas. Snow was collected from the surface to represent collection during melt water production for hot-water drilling. At each sampling location, a 1 m snow pit was excavated and the top 30 cm of snow was sampled using an ethanol sterilized shovel and Whirl-Pak bags (Nasco, WI, United States). Samples were transported to the field laboratory, where they were melted and passed through 0.2 μm Sterivex filters (Merck, Darmstadt, Germany), before freezing at −20°C for transport and processing in the United Kingdom.
FIGURE 1
PMA Treatment and DNA Extractions
SkyBlu, Pine Island Bay and Ellsworth snow samples were processed in two ways; PMA-treated and non-PMA-treated in order to differentiate the potentially viable microbial community from relic DNA. Each 0.2 μm filter was cut in half using sterile and DNAase/ethanol treated razors. For each sample, one-half was processed with PMA as per and , using a 20 mM stock solution of PMA (Biotium, Hayward, CA, United States) in a 20% (v/v) aqueous solution of dimethyl sulfoxide (DMSO). Filters for PMA treatment were placed in a 6-well plate and a PMA solution at a final concentration of 100 μM was added. Cross-linking was initiated by 10 min incubation on ice, in the dark with occasional mixing, followed by 5 min of light exposure using a 650 W halogen lamp (FLASH 2000 L, DTS, Italy), at a 20 cm distance. The process was carried out in a laminar flow hood to avoid contamination of the samples. Non-PMA-treated samples were incubated in a 20% (v/v) aqueous solution of DMSO, and treated following the same protocol as PMA-treated samples. All samples were washed twice with sterile phosphate buffered saline (PBS) and all samples were then used for DNA extraction.
DNA from snow samples was extracted using the PowerWater kit from MoBio (Qiagen, Carlsbad, CA, United States). Each sample was PCR amplified using the primers 341F and 785R covering the V3-V4 regions of the 16S rRNA gene (), under the following conditions: initial denaturation at 95°C for 5 min then 25 cycles of 40 s denaturation at 95°C; primer annealing at 55°C for 2 min; and elongation at 72°C for 1 min then a final elongation at 72°C for 7 min (). DNA extraction kit controls were included alongside the snow derived DNA and sequenced. PCR amplicons were cleaned, normalized, quantified and supplemented with 5% PhiX before being loaded on Illumina MiSeq, as per the Illumina standard operating protocol ().
Illumina Sequencing and Data Processing
Forward and reverse reads longer than 420 bp were merged using FLASH (fast length adjustment of short reads) () for a total of 4.2 million reads. Vsearch (Rognes et al., 2016) was used for downstream analyses. Quality filtering by expected errors was run to keep only sequences with high Phred scores. Dereplication was used to identify unique sequences. Chimeras were removed using a two-step chimera detection method; first by aligning against ChimeraSlayer Gold database provided with SILVA (Pruesse et al., 2007), and secondly by using the de novo detection module in Vsearch. De novo operational taxonomic unit (OTU) calling was carried out at a 97% similarity level to generate operational taxonomical units (OTUs), and aligned using the Python Nearest Alignment Space Termination (PyNAST) () tool with a relaxed neighbor-joining tree built using FastTree (Price et al., 2010). The taxonomy was determined using the Classification Resources for Environmental Sequence Tags (CREST) () classifier with a confidence threshold of 0.80 against SILVA release 128 as a reference database. A total of 42 DNA samples, equal to 3 million reads (63,202 ± 60,525 reads/sample), were assigned against 6,793 OTUs using the USEARCH () algorithm for clustering with a threshold of 97 % similarity. Datasets are deposited on the European Nucleotide Archive under the accession numbers SMAN06368443 to SAMN06368447 and PRJEB29215.
Statistical Analysis
All statistical analyses were performed with a combination of Qiime1 V 1.90 () and R environment (R Core Team, 2013) using phyloseq (), vegan () and indicspecies () packages. We calculated the alpha diversity using matrices of richness (number of observed OTUs) and diversity (Shannon and Simpson indices) based on OTU tables rarefied to the lower number of reads. We tested differences in alpha diversity using ANOVA and Tukey honestly significant differences (HSD) test to define which location exhibited significant differences. Beta diversity using Bray-Curtis distance was then calculated by normalizing the OTU table using cumulative-sum scaling (CSS), in which raw counts are divided by the cumulative sum (Paulson et al., 2013). The dissimilarity matrix was plotted using principal coordinates analysis (PCoA). ANOSIM from vegan was used to analyze the similarities based on a Bray-Curtis dissimilarity matrix across locations with free permutations. Indicator species were determined by the Dufrene-Legendre indicator species analysis method () to identify OTUs that were specifically associated with the different geographic locations.
Results
The Surface Snow
We first investigated the bacterial communities excluding the Ellsworth field camp samples from the analysis as these were specifically chosen to maximize the chance of detecting some human impact. The alpha diversity analysis identified geographic differences in the number of observed OTUs (ANOVA F6,25 = 4.49, p = 0.0032) and Tukey’s HSD test suggested a latitudinal gradient with significant differences between Signy Island and the Ellsworth samples (Tukey HSD, p = 0.016). The Shannon and Simpson indices showed that samples from Livingston were more diverse and even that at any other location. Overall, Pine Island Bay, SkyBlu and Ellsworth harbored fewer OTUs than Livingston and Signy Island samples (Figure 2A–C).
FIGURE 2
Principal component analysis of bacterial communities identified three clusters (Figure 3A), corresponding to three locations: Livingston Island, Signy Island and continental Antarctica (encompassing the Lake Ellsworth area, Pine Island Bay and SkyBlu). A heatmap and dendrogram (Figure 3B) further illustrated the grouping into the three distinct clusters by geographic locations; each Antarctic island and continental Antarctica. Not only were each of the clusters different, samples within clusters differed in relative abundance and composition, as highlighted in the PCoA (Figure 3A) and heatmap (Figure 3B) (ANOSIM R = 0.782, p < 0.001).
FIGURE 3
Of 6,793 taxa identified, only 126 were characterized as abundant taxa, with relative abundances over 0.1%. Proteobacteria, mainly Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria dominated all samples collected. Bacteroidetes were in high abundance in Antarctic island samples but virtually absent from continental Antarctica samples. On the other hand, high abundances of Firmicutes were recorded in continental Antarctica (Figure 4A,B). Actinobacteria decreased in abundance at SkyBlu and Pine Island Bay but not in the Ellsworth samples. Bacteroidetes changed from Cytophaga, Flavobacteria, Sphingobacteriia, and Chitinophagia in northern sites, through unidentified at SkyBlu and Pine Island Bay, to Chitinophaga only at the southern sites.
FIGURE 4
Overall, at the genus level, distinctions between each site were clear (Supplementary Figure S1). A few instances of notable variations included the high proportion of Methylobacterium and Rhizobium in continental Antarctica (up to 20 and 13%, respectively), which are present in very low abundances in the Antarctic islands.
The core microbiome, as comprised of OTUs identified in over 80% of all samples analyzed, was represented by only one OTU, classified as Bradyrhizobiaceae, at a total relative abundance of 0.52%. The indicator species analysis identified OTUs specifically associated with geographic locations. Although over six thousand OTUs were detected in total, only 38 were associated with Signy Island, mainly in the Alphaproteobacteria, Actinobacteria and Bacteroidetes. 48 indicator species OTUs were identified among Livingston Island samples, mainly Betaproteobacteria, Bacteroidetes and Alphaproteobacteria. No indicator species were identified for continental Antarctic sites, further highlighting the difference in microbial communities between each location, regardless of geographic distances.
Cell Viability and Impact of Relic DNA on Snow Microbiological Diversity
To investigate the impact of relic DNA from dead cells, on microbial diversity of Antarctic snow, the SkyBlu, Pine Island Bay, Ellsworth pristine, camp and drill samples were treated with PMA before DNA extraction and sequencing. The number of OTUs, Shannon index and Simpson index were all lower in PMA-treated samples than in non-PMA-treated samples (Supplementary Figure S2), showing the presence of relic DNA in snow. While the number of OTUs identified decreased by up to 50% after PMA treatment in some cases (Supplementary Figure S2A), differences in diversity indices by treatment were not significant (ANOVA F1,18 = 1.1, p > 0.05) (Supplementary Figures S2B,C). Similarly, the community composition did not significantly change after PMA treatment (ANOSIM, R = 0.06, p = 0.125). Differences in diversity were mainly characterized by the decrease in Actinobacteria, an increase in Proteobacteria and variation in Firmicute abundance in PMA-treated samples (Figure 5).
FIGURE 5
The Lake Ellsworth Site
For this site, the Antarctic snow community was determined from an area of current scientific research interest. We sampled surface snow within the Lake Ellsworth deep-field camp in December 2012 (Pearce D.A. et al., 2016), close to the kitchen, generators and the deep hot-water drilling head, and away from the camp to use as a proxy for pristine snow communities. This deep field site was characterized by human presence, activity and disturbance of the natural environment. This is reflected by the isolation of the drilling head site in the PCoA (Figure 2A and Supplementary Figure 3), with an unique microbial assemblage. While both off-site and camp samples were highly similar in richness and diversity, the drilling site harbored a higher number of OTUs and appeared much more diverse than any other site sampled, including Signy and Livingston Islands (Figure 2). Differences in the microbial communities of the drilling site were apparent at the genus level (Supplementary Figure 1), and this is further illustrated by the indicator species analysis, which identified 49 drill-associated taxa, mainly classified as Firmicutes, Actinobacteria and Alphaproteobacteria. A BLAST search identified many of these taxa as soil-related, often with 99% identity match to worldwide soils, from Brazil to the Chinese mountains of Yanshan, to Moab National Park and Mt. Kilimanjaro.
Discussion
Our study of microbial communities in Antarctic snow identified three clusters: Livingston Island, Signy Island and continental Antarctica, which included Pine Island Bay, SkyBlu and the Lake Ellsworth region (Figure 3). The differences likely result from geographical distance and variation in environmental conditions from the Antarctic islands to continental Antarctica. Livingston Island has an annual mean temperature of −1°C by the coast and is characterized by seasonal snow cover and low altitude (). Similarly, Signy Island has mean annual temperatures around −3.5°C also with seasonal snow cover (). The sampling locations around continental Antarctica are at high altitude, with annual mean temperatures between −20 and −30°C () and permanent snow cover. Therefore, the difference in snow microbial communities between the Antarctic islands and continental Antarctica are likely due to snow type, density and physicochemical properties (; ). Differences observed between sites within the same cluster reflect natural variation of microbial communities, which have been observed down to the pore scale in soils (Ruamps et al., 2011).
Overall, Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria were abundant in all samples, which is in accordance with previous studies on the Antarctic snow microbial community (; ; ). However, compared to previous studies, Firmicutes were consistently identified in all samples and Cyanobacteria were generally below the 0.1% abundance threshold. Firmicutes are known to harbor a number of survival strategies such as biofilm formation or the formation of endospores. These strategies allow them to survive extreme environmental conditions such as physicochemical gradients, UV radiation and desiccation ().
In all samples, a large proportion of microorganisms were identified as soil-related. The analysis of extraction kit controls suggested that they were unlikely to be the result of kit contamination, as all samples, except SkyBlu and Blank 2, had more than 2000 reads. Beside SkyBlu and Blank2, PMA-treated samples and Blank1 had the smallest number of reads (below 10,000 reads). Furthermore, the blanks clustered away from other samples in the PCoA, illustrating the differences in microbial communities (Supplementary Figure 3). All other samples had over 12,000 reads; above the 2,000 minimum reads required to attain reproducible and reliable patterns of diversity ().
The study by identified soil microorganisms as a primary source of the snow microbial assemblage of the West Greenland Ice Sheet, while studies on airborne microbial communities have reported soil ecosystems as a primary source of microorganisms in the Arctic (; Šantl-Temkiv et al., 2018). Airborne microbial communities are consistently dominated by Firmicutes, Alpha-, Beta-, and Gamma- Proteobacteria () including in Antarctica (Pearce et al., 2010; ). Spore-forming genera (i.e., Bacillus spp.) are a major component of airborne microbial communities and are associated with long-distance dispersal (). Our study shows that snow and airborne microbial communities appear similar, especially in continental Antarctica, including a high abundance of Bacilli (Figure 4B) suggesting the importance of long-distance dispersal in seeding continental Antarctic snow ecosystems. However, we identified a small abundance of Bacteroidetes in the Ellsworth snow samples. This could reflect small terrestrial inputs from the Ellsworth Mountains, dominated by Bacteroidetes (Yergeau et al., 2007) and the closest terrestrial system, located approximately 70 km away from the sampling site (Pearce D.A. et al., 2016).
Antarctic island snow communities had higher abundances of Bacteroidetes and Actinobacteria and lower abundances of Firmicutes. Antarctic soil microbial communities from both Signy and Livingston Islands, harbored high abundances of Bacteroidetes, including Sphingobacteria and Chitinophagia (Yergeau et al., 2007; , ; ), also identified in the collected snow samples (Figure 4B). This may reflect local dispersal mechanisms, such as vector-based dispersal through the local wildlife (), or aerosolization of local soil bacteria dispersed over short-distances (). It should be noted that soil systems occur on Antarctic islands and the snow cover is seasonal, while in continental Antarctica, snow cover is permanent, soil ecosystems are sparse and, often, the closest source is the coastline or the Antarctic islands themselves. Overall, the observed microbial communities and clustering pattern support the hypothesis that soils are the primary source of snow microorganisms through local or long-range dispersal. While Antarctic islands snow microbial communities may be influenced by both pathways (Vincent, 2000), continental Antarctica appears to be seeded by air-borne microorganisms through long-range aeolian transport.
The indicator species analysis identified limited OTUs with geographic associations. The majority of indicator OTUs were associated with the Antarctic islands and far fewer with continental Antarctica. This analysis suggests that Signy and Livingston Islands have unique, geographically associated communities, but locally homogeneous microbial communities, likely influenced by surrounding environmental conditions. Continental sites with harsher environmental conditions, primarily seeded by air-borne microorganisms, had more variable microbial communities, even within the same location, likely dependant on weather and wind patterns.
We assessed the presence of relic DNA in snow samples collected from continental Antarctica. Relic DNA coming from dead cells, has been identified representing on average 40% of DNA extracted from soil samples (; ). While studies of microbial communities are essential to understand functional activities and ecological roles of microorganisms, relic DNA could lead to false estimations of putative functions and lead to misunderstanding of the role of microorganisms in geochemical cycling. Identifying relic DNA in environmental samples is relatively new and focused on soils (; ), where microbial biomass is high and concentrated. As samples were filtered prior to PMA treatment and extraction, free environmental DNA would have been discarded during the filtering process, unless attached to particles and, as such, relic DNA would only be present in low abundance in snow samples. Yet, in some cases, the number of reads and identified OTUs was reduced by 50% while, in others, differences were imperceptible. Despite variation in OTUs between non-PMA and PMA-treated samples, alpha diversity and community composition did not change significantly. While relic DNA can deeply affect some microbial community analyses, relic DNA is not a major issue in snow samples.
Finally, we did not identify directly human-related taxa in any of the samples. This is likely due to environmental conditions preventing human bacteria from colonizing Antarctic snow ecosystems rather than the lack of contamination of the surrounding environment (). However, strong winds and snow movements could have potentially transported human-related bacteria to distant sites and, thus, would remain undetected in this study. Furthermore, we used sites from distant locations as a proxy for the initial indigenous microbial populations present before human disturbance at the Ellsworth camp site. The high similarity of microbial assemblages in the Ellsworth off-site samples and camp samples suggests that human presence may not have a major impact on snow microbial communities. However, strong human-derived disturbances and activities, such as the import of equipment, snow digging, intense stepping, and sitting and manipulating instruments, may lead to changes in natural assemblages, as illustrated at the drilling site. While it is unlikely that the drilling site communities differ from the other samples purely by chance, we cannot exclude it completely as we did not sample the site prior to the research activity. In addition to human activities impacting the indigenous snow microbial community, the presence of exogenous microorganisms on camp and drilling material could also have an impact. This could be mitigated in future by careful and thorough cleaning and sterilization of equipment prior to landing in Antarctica and throughout the drilling process (Priscu et al., 2013; ). Due to the observed differences in microbial communities from site to site, we recommend indigenous surface microbial communities be considered in planning deep-field subsurface access operations and characterized prior to research activities.
Conclusion
This study, of spatial variation in microbial communities from Signy Island down toward the Ellsworth Mountains, covers the largest geographical extent of any investigation of Antarctic snow microbial communities. We identified three distinct geographic regions from the analysis, with each Antarctic island as individual clusters, and all sites of continental Antarctica clustered together. While each cluster harbored distinct microbial assemblages, within-cluster differences were also observed, illustrating the spatial variability of Antarctic snow bacterial communities. Overall, communities were primarily composed of soil-related microorganisms, suggesting colonization from both local and distant terrestrial ecosystems. The Lake Ellsworth case study suggests that while low-level human presence may not affect natural systems, intense human activities and ecosystem disturbances might. When undertaking research aimed at detecting life in extreme environments, such as in subglacial lakes, we believe it is important to consider how the indigenous microbial communities may be affected across the field site as well as with in subglacial targets themselves (Siegert and Kennicutt, 2018).
Statements
Author contributions
LM conducted the bioinformatics processing, statistical analyses, and wrote the manuscript. MS conducted the fieldwork, laboratory work of the Signy samples, and contributed to bioinformatic processing. IM and MM conducted the laboratory work from the Ellsworth, Pine Island Bay, and SkyBlu samples. MT participated in the Ellsworth sampling campaign. MJS was the Principal Investigator of the Lake Ellsworth exploration program. DP designed the study, co-wrote the manuscript and co-ordinated the sampling campaign. All authors have read, revised, and approved the manuscript.
Funding
This work was supported by grants from the United Kingdom Natural Environment Research Council grants G00465X/1, G00465X/2, G00465X/3, NE/H014446/1, and NE/H014802/1 and from the European Commission’s Marie Skłodowska-Curie Actions program under project number 675546.
Acknowledgments
We acknowledge the Lake Ellsworth Consortium, the British Antarctic Survey logistics and field support staff, air unit and base personnel at Rothera Research Station.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00461/full#supplementary-material
FIGURE S1Stacked bar chart of the bacterial diversity by sampling site.
FIGURE S2Plots of Alpha diversity measures by Location and by Treatment. (A) Observed OTUs. (B) Shannon Index. (C) Simpson Index.
FIGURE S3PCoA of microbial communities focused on the continental Antarctica cluster. It includes samples from the Ellsworth Camp, drilling site and extraction blanks clustered away.
TABLE S1Metadata table of sample sites and GPS locations.
References
1
AntonyR.SanyalA.KapseN.DhakephalkarP. K.ThambanM.NairS. (2016). Microbial communities associated with Antarctic snow pack and their biogeochemical implications.Microbiol. Res.192192–202. 10.1016/j.micres.2016.07.004
2
AzzoniR. S.TagliaferriI.FranzettiA.MayerC.LambrechtA.CompostellaC.et al (2018). Bacterial diversity in snow from mid-latitude mountain areas: alps, Eastern Anatolia, Karakoram and Himalaya.Anna. Glaciol.1–11.
3
BertlerN.MayewskiP. A.AristarainA.BarrettP.BecagliS.BernardoR.et al (2005). Snow chemistry across Antarctica.Anna. Glaciol.41167–179. 10.3189/172756405781813320
4
BoetiusA.AnesioA. M.DemingJ. W.MikuckiJ. A.RappJ. Z. (2015). Microbial ecology of the cryosphere: sea ice and glacial habitats.Nat. Rev. Microbiol.13677–690. 10.1038/nrmicro3522
5
BohlanderJ.ScambosT. (2001). THERMAP Antarctic Ice Sheet Temperature Data.Boulder, CO: National Snow and Ice Data Centre.
6
BottosE. M.WooA. C.Zawar-RezaP.PointingS. B.CaryS. C. (2014). Airborne bacterial populations above desert soils of the McMurdo Dry Valleys, Antarctica.Microb. Ecol.67120–128. 10.1007/s00248-013-0296-y
7
BulatS. A. (2016). Microbiology of the subglacial Lake Vostok: first results of borehole-frozen lake water analysis and prospects for searching for lake inhabitants.Philos. Trans. R. Soc. A374:20140292. 10.1098/rsta.2014.0292
8
CáceresM. D.LegendreP. (2009). Associations between species and groups of sites: indices and statistical inference.Ecology903566–3574. 10.1890/08-1823.1
9
CallaghanT. V.JohanssonM.BrownR. D.GroismanP. Y.LabbaN.RadionovV.et al (2011). Multiple effects of changes in Arctic snow cover.Ambio4032–45. 10.1007/s13280-011-0213-x
10
CameronK. A.HagedornB.DieserM.ChristnerB. C.ChoquetteK.SlettenR.et al (2015). Diversity and potential sources of microbiota associated with snow on western portions of the G reenland I ce S heet.Environ. Microbiol.17594–609. 10.1111/1462-2920.12446
11
CaporasoJ. G.BittingerK.BushmanF. D.DeSantisT. Z.AndersenG. L.KnightR. (2009). PyNAST: a flexible tool for aligning sequences to a template alignment.Bioinformatics26266–267. 10.1093/bioinformatics/btp636
12
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al (2010a). QIIME allows analysis of high-throughput community sequencing data.Nat. Methods7:335. 10.1038/nmeth.f.303
13
CaporasoJ. G.LauberC. L.WaltersW. A.Berg-LyonsD.LozuponeC. A.TurnbaughP. J.et al (2010b). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample.Proc. Natl. Acad. Sci. U.S.A.108(Suppl. 1), 4516–4522. 10.1073/pnas.1000080107
14
CariniP.MarsdenP. J.LeffJ. W.MorganE. E.StricklandM. S.FiererN. (2017). Relic DNA is abundant in soil and obscures estimates of soil microbial diversity.Nat. Microbiol.2:16242. 10.1038/nmicrobiol.2016.242
15
CarpenterE. J.LinS.CaponeD. G. (2000). Bacterial activity in South Pole snow.Appl. Environ. Microbiol.664514–4517. 10.1128/AEM.66.10.4514-4517.2000
16
ChongC. W.DunnM. J.ConveyP.TanG. A.WongR. C.TanI. K. (2009). Environmental influences on bacterial diversity of soils on Signy Island, maritime Antarctic.Polar Biol.321571–1582. 10.1007/s00300-009-0656-8
17
ChongC. W.PearceD. A.ConveyP.TanG. A.WongR. C.TanI. K. (2010). High levels of spatial heterogeneity in the biodiversity of soil prokaryotes on Signy Island, Antarctica.Soil Biol. Biochem.42601–610. 10.1016/j.soilbio.2009.12.009
18
ChristnerB. C.Mosley-ThompsonE.ThompsonL. G.ZagorodnovV.SandmanK.ReeveJ. N. (2000). Recovery and identification of viable bacteria immured in glacial ice.Icarus144479–485. 10.1006/icar.1999.6288
19
CowanD. A.ChownS. L.ConveyP.TuffinM.HughesK.PointingS.et al (2011). Non-indigenous microorganisms in the Antarctic: assessing the risks.Trends Microbiol.19540–548. 10.1016/j.tim.2011.07.008
20
CuthbertsonL.Amores-ArrochaH.MalardL. A.ElsN.SattlerB.PearceD. A. (2017). characterisation of arctic bacterial communities in the air above Svalbard.Biology6:29. 10.3390/biology6020029
21
DixonP. (2003). VEGAN, a package of R functions for community ecology.J. Veg. Sci.14927–930. 10.1111/j.1654-1103.2003.tb02228.x
22
EdgarR. C. (2010). Search and clustering orders of magnitude faster than BLAST.Bioinformatics262460–2461. 10.1093/bioinformatics/btq461
23
EvelandJ.GooseffM. N.LampkinD. J.BarrettJ.Takacs-VesbachC. (2013). Spatial and temporal patterns of snow accumulation and aerial ablation across the McMurdo Dry Valleys, Antarctica.Hydrol. Process.272864–2875.
24
FiererN. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome.Nat. Rev. Microbiol.15:579. 10.1038/nrmicro.2017.87
25
FiererN.LiuZ.Rodríguez-HernándezM.KnightR.HennM.HernandezM. T. (2008). Short-term temporal variability in airborne bacterial and fungal populations.Appl. Environ. Microbiol.74200–207. 10.1128/AEM.01467-07
26
FilippidouS.WunderlinT.JunierT.JeanneretN.DoradorC.MolinaV.et al (2016). A combination of extreme environmental conditions favor the prevalence of endospore-forming firmicutes.Front. Microbiol.7:1707. 10.3389/fmicb.2016.01707
27
FittipaldiM.NockerA.CodonyF. (2012). Progress in understanding preferential detection of live cells using viability dyes in combination with DNA amplification.J. Microbiol. Methods91276–289. 10.1016/j.mimet.2012.08.007
28
FranceJ.KingM.FreyM.ErblandJ.PicardG.PreunkertS.et al (2011). Snow optical properties at Dome C (Concordia), Antarctica; implications for snow emissions and snow chemistry of reactive nitrogen.Atmospher. Chem. Phys.119787–9801. 10.5194/acp-11-9787-2011
29
FretwellP.PritchardH. D.VaughanD. G.BamberJ.BarrandN.BellR.et al (2013). Bedmap2: improved ice bed, surface and thickness datasets for Antarctica.Cryosphere7375–393. 10.5194/tc-7-375-2013
30
GanzertL.BajerskiF.WagnerD. (2014). Bacterial community composition and diversity of five different permafrost-affected soils of Northeast Greenland.FEMS Microbiol. Ecol.89426–441. 10.1111/1574-6941.12352
31
GreenJ.BohannanB. J. (2006). Spatial scaling of microbial biodiversity.Trends .Ecol. Evol.21501–507. 10.1016/j.tree.2006.06.012
32
GuglielminM.EvansC. J. E.CannoneN. (2008). Active layer thermal regime under different vegetation conditions in permafrost areas. A case study at Signy Island (Maritime Antarctica).Geoderma14473–85. 10.1016/j.geoderma.2007.10.010
33
HardingT.JungblutA. D.LovejoyC.VincentW. F. (2011). Microbes in High Arctic snow and implications for the cold biosphere.Appl. Environ. Microbiol.773234–3243. 10.1128/AEM.02611-10
34
HellK.EdwardsA.ZarskyJ.PodmirsegS. M.GirdwoodS.PachebatJ. A.et al (2013). The dynamic bacterial communities of a melting High Arctic glacier snowpack.ISME J.7:1814. 10.1038/ismej.2013.51
35
HodsonA.NowakA.CookJ.SabackaM.WharfeE.PearceD.et al (2017). Microbes influence the biogeochemical and optical properties of maritime Antarctic snow.J. Geophys. Res. Biogeosci.1221456–1470. 10.1002/2016JG003694
36
HughesK. A.ThompsonA. (2004). Distribution of sewage pollution around a maritime Antarctic research station indicated by faecal coliforms, Clostridium perfringens and faecal sterol markers.Environ. Pollut.127315–321. 10.1016/j.envpol.2003.09.004
37
JonesH. (1999). The ecology of snow-covered systems: a brief overview of nutrient cycling and life in the cold.Hydrol. Process.132135–2147. 10.1002/(SICI)1099-1085(199910)13:14/15<2135::AID-HYP862>3.0.CO;2-Y
38
KelloggC. A.GriffinD. W.GarrisonV. H.PeakK. K.RoyallN.SmithR. R.et al (2004). Characterization of aerosolized bacteria and fungi from desert dust events in Mali, West Africa.Aerobiologia2099–110. 10.1023/B:AERO.0000032947.88335.bb
39
KlindworthA.PruesseE.SchweerT.PepliesJ.QuastC.HornM.et al (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies.Nucleic Acids Res.41:e1. 10.1093/nar/gks808
40
KozichJ.SchlossP.BaxterN.JeniorM.KoumpourasC. (2013). 16S rRNA Sequencing with the Illumina MiSeq: Library Generation, QC, & Sequencing. Available at: https://aem.asm.org/content/aem/suppl/2013/08/06/AEM.01043-13.DCSupplemental/zam999104626so1.pdf.
41
KuhnM. (2001). The nutrient cycle through snow and ice, a review.Aquat. Sci.63150–167. 10.1007/PL00001348
42
LanzénA.JørgensenS. L.HusonD. H.GorferM.GrindhaugS. H.JonassenI.et al (2012). CREST–classification resources for environmental sequence tags.PLoS One7:e49334. 10.1371/journal.pone.0049334
43
LaroseC.DommergueA.VogelT. M. (2013). The dynamic arctic snow pack: an unexplored environment for microbial diversity and activity.Biology2317–330. 10.3390/biology2010317
44
LennonJ. T.MuscarellaM. E.PlacellaS. A.LehmkuhlB. K. (2018). How, when, and where relic DNA affects microbial diversity.MBio9:e00637-18. 10.1128/mBio.00637-18
45
LopatinaA.MedvedevaS.ShmakovS.LogachevaM. D.KrylenkovV.SeverinovK. (2016). Metagenomic analysis of bacterial communities of Antarctic surface snow.Front. Microbiol.7:398. 10.3389/fmicb.2016.00398
46
MaccarioL.VogelT. M.LaroseC. (2014). Potential drivers of microbial community structure and function in Arctic spring snow.Front. Microbiol.5:413. 10.3389/fmicb.2014.00413
47
MagoèT.SalzbergS. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies.Bioinformatics272957–2963. 10.1093/bioinformatics/btr507
48
MakinsonK.PearceD.HodgsonD. A.BentleyM. J.SmithA. M.TranterM.et al (2016). Clean subglacial access: prospects for future deep hot-water drilling.Philos. Trans. R. Soc. A374:20140304. 10.1098/rsta.2014.0304
49
MargesinR.ZackeG.SchinnerF. (2002). Characterization of heterotrophic microorganisms in alpine glacier cryoconite.Arct. Antarct. Alpine Res.3488–93. 10.1080/15230430.2002.12003472
50
MartinyJ. B. H.BohannanB. J.BrownJ. H.ColwellR. K.FuhrmanJ. A.GreenJ. L.et al (2006). Microbial biogeography: putting microorganisms on the map.Nat. Rev. Microbiol.4102–112. 10.1038/nrmicro1341
51
Masson-DelmotteV.HouS.EkaykinA.JouzelJ.AristarainA.BernardoR.et al (2008). A review of Antarctic surface snow isotopic composition: observations, atmospheric circulation, and isotopic modeling.J. Clim.213359–3387. 10.1175/2007JCLI2139.1
52
McMurdieP. J.HolmesS. (2013). phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data.PLoS One8:e61217. 10.1371/journal.pone.0061217
53
MichaudL.GiudiceA. L.MysaraM.MonsieursP.RaffaC.LeysN.et al (2014). Snow surface microbiome on the High Antarctic Plateau (DOME C).PLoS One9:e104505. 10.1371/journal.pone.0104505
54
MitevaV. (2008). “Bacteria in snow and glacier ice,” in Psychrophiles: From Biodiversity to Biotechnology, edsSchinnerF.MargesinR.MarxJ. C. (Berlin: Springer), 31–50. 10.1007/978-3-540-74335-4_3
55
MøllerA. K.SøborgD. A.Abu Al-SoudW.SørensenS. J.KroerN. (2013). Bacterial community structure in High-Arctic snow and freshwater as revealed by pyrosequencing of 16S rRNA genes and cultivation.Polar Res.32:17390. 10.3402/polar.v32i0.17390
56
NockerA.CamperA. K. (2009). Novel approaches toward preferential detection of viable cells using nucleic acid amplification techniques.FEMS Microbiol. Lett.291137–142. 10.1111/j.1574-6968.2008.01429.x
57
PaulsonJ. N.StineO. C.BravoH. C.PopM. (2013). Differential abundance analysis for microbial marker-gene surveys.Nat. Methods10:1200. 10.1038/nmeth.2658
58
PearceD.MagiopoulosI.MowlemM.TranterM.HoltG.WoodwardJ.et al (2016). Microbiology: lessons from a first attempt at Lake Ellsworth.Philos. Trans. R. Soc. A.374:20140291. 10.1098/rsta.2014.0291
59
PearceD. A.AlekhinaI. A.TeraudsA.WilmotteA.QuesadaA.EdwardsA.et al (2016). Aerobiology over Antarctica–a new initiative for atmospheric ecology.Front. Microbiol.7:16. 10.3389/fmicb.2016.00016
60
PearceD. A.HughesK.Lachlan-CopeT.HarangozoS.JonesA. E. (2010). Biodiversity of air-borne microorganisms at Halley station, Antarctica.Extremophiles14145–159. 10.1007/s00792-009-0293-8
61
PriceM. N.DehalP. S.ArkinA. P. (2010). FastTree 2–approximately maximum-likelihood trees for large alignments.PLoS One5:e9490. 10.1371/journal.pone.0009490
62
PriscuJ. C.AchbergerA. M.CahoonJ. E.ChristnerB. C.EdwardsR. L.JonesW. L.et al (2013). A microbiologically clean strategy for access to the Whillans Ice Stream subglacial environment.Antarct. Sci.25637–647. 10.1017/S0954102013000035
63
PruesseE.QuastC.KnittelK.FuchsB. M.LudwigW.PepliesJ.et al (2007). SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB.Nucleic Acids Res.357188–7196. 10.1093/nar/gkm864
64
R Core Team (2013). R: A Language and Environment for Statistical Computing. Available at: https://www.gbif.org/tool/81287/r-a-language-and-environment-for-statistical-computing
65
RignotE.MouginotJ.MorlighemM.SeroussiH.ScheuchlB. (2014). Widespread, rapid grounding line retreat of Pine Island, Thwaites, Smith, and Kohler glaciers, West Antarctica, from 1992 to 2011.Geophys. Res. Lett.413502–3509. 10.1002/2014GL060140
66
RognesT.FlouriT.NicholsB.QuinceC.MahéF. (2016). VSEARCH: a versatile open source tool for metagenomics.PeerJ4:e2584. 10.7717/peerj.2584
67
RuampsL. S.NunanN.ChenuC. (2011). Microbial biogeography at the soil pore scale.Soil Biol. Biochem.43280–286. 10.1111/j.1574-6941.2008.00490.x
68
Šantl-TemkivT.GosewinkelU.StarnawskiP.LeverM.FinsterK. (2018). Aeolian dispersal of bacteria in southwest Greenland: their sources, abundance, diversity and physiological states.FEMS Microbiol. Ecol.94:fiy031. 10.1093/femsec/fiy031
69
SiegertM.KennicuttM. (2018). Governance of the exploration of subglacial antarctica.Front. Environ. Sci.6:103. 10.3389/fenvs.2018.00103
70
SiegertM. J.MakinsonK.BlakeD.MowlemM.RossN. (2014). An assessment of deep hot-water drilling as a means to undertake direct measurement and sampling of Antarctic subglacial lakes: experience and lessons learned from the Lake Ellsworth field season 2012/13.Ann. Glaciol.5559–73. 10.3189/2014AoG65A008
71
SiegertM. J.PriscuJ. C.AlekhinaI. A.WadhamJ. L.LyonsW. B. (2016). Antarctic Subglacial Lake Exploration: First Results and Future Plans.London: The Royal Society Publishing.
72
SjolingS.CowanD. (2000). Presence of human-specific enteric micro-organisms in current and historic field camp sites.Polar Biol.23644–650.
73
StaleyJ. T.GosinkJ. J. (1999). Poles apart: biodiversity and biogeography of sea ice bacteria.Annu. Rev. Microbiol.53189–215. 10.1146/annurev.micro.53.1.189
74
TulaczykS.MikuckiJ. A.SiegfriedM. R.PriscuJ. C.BarcheckC. G.BeemL. H.et al (2014). WISSARD at Subglacial Lake Whillans, West Antarctica: scientific operations and initial observations.Anna. Glaciol.5551–58. 10.3189/2014AoG65A009
75
VavrusS. (2007). The role of terrestrial snow cover in the climate system.Clim. Dyn.2973–88. 10.1371/journal.pone.0137275
76
VincentW. F. (1988). Microbial Ecosystems of Antarctica.Cambridge: Cambridge University Press.
77
VincentW. F. (2000). Evolutionary origins of Antarctic microbiota: invasion, selection and endemism.Antarct. Sci.12374–385. 10.1017/S0954102000000420
78
WunderlinT.FerrariB.PowerM. (2016). Global and local-scale variation in bacterial community structure of snow from the Swiss and Australian Alps.FEMS Microbiol. Ecol.92:fiw132. 10.1093/femsec/fiw132
79
YergeauE.NewshamK. K.PearceD. A.KowalchukG. A. (2007). Patterns of bacterial diversity across a range of Antarctic terrestrial habitats.Environ. Microbiol.92670–2682. 10.1111/j.1462-2920.2007.01379.x
Summary
Keywords
Antarctic, snow, biogeography, Ellsworth Lake, microbial diversity, relic DNA
Citation
Malard LA, Šabacká M, Magiopoulos I, Mowlem M, Hodson A, Tranter M, Siegert MJ and Pearce DA (2019) Spatial Variability of Antarctic Surface Snow Bacterial Communities. Front. Microbiol. 10:461. doi: 10.3389/fmicb.2019.00461
Received
19 October 2018
Accepted
21 February 2019
Published
26 March 2019
Volume
10 - 2019
Edited by
David Velazquez, Autonomous University of Madrid, Spain
Reviewed by
Charles K. Lee, University of Waikato, New Zealand; Connie Lovejoy, Laval University, Canada
Updates
Copyright
© 2019 Malard, Šabacká, Magiopoulos, Mowlem, Hodson, Tranter, Siegert and Pearce.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: David A. Pearce, david.pearce@northumbria.ac.uk
This article was submitted to Extreme Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.