Next Article in Journal
Incorporating the Soil Gas Gradient Method and Functional Genes to Assess the Natural Source Zone Depletion at a Petroleum-Hydrocarbon-Contaminated Site of a Purification Plant in Northwest China
Previous Article in Journal
Predicting Factors for Seizures after Cerebral Venous Thrombosis: A Retrospective Single Center Cohort Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 13
Issue 1
10.3390/life13010113
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Snow Surface Microbial Diversity at the Detection Limit within the Vicinity of the Concordia Station, Antarctica

by
Alessandro Napoli
1,2,†,
Claudia Coleine
3,†,
Nikea J. Ulrich
4,
Ralf Moeller
5,6,*,
Daniela Billi
1,* and
Laura Selbmann
3,7
1
Department of Biology, University of Rome Tor Vergata, 00133 Rome, Italy
2
Ph.D. Program in Cellular and Molecular Biology, Department of Biology, University of Rome Tor Vergata, 00133 Rome, Italy
3
Department of Ecological and Biological Sciences, University of Tuscia, 01100 Viterbo, Italy
4
Division of Biological Sciences, University of Montana, Missoula, MT 59812, USA
5
Aerospace Microbiology Research Group, Radiation Biology Department, Institute of Aerospace Medicine, German Aerospace Center (DLR), 28359 Cologne, Germany
6
Department of Natural Sciences, University of Applied Sciences Bonn-Rhein-Sieg (BRSU), 53359 Rheinbach, Germany
7
Mycological Section, Italian Antarctic National Museum (MNA), 16128 Genoa, Italy
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Life 2023, 13(1), 113; https://doi.org/10.3390/life13010113
Submission received: 4 November 2022 / Revised: 17 December 2022 / Accepted: 23 December 2022 / Published: 30 December 2022
(This article belongs to the Section Astrobiology)
Browse Figures
Versions Notes

Abstract

:
The Concordia Research Station provides a unique location for preparatory activities for future human journey to Mars, to explore microbial diversity at subzero temperatures, and monitor the dissemination of human-associated microorganisms within the pristine surrounding environment. Amplicon sequencing was leveraged to investigate the microbial diversity of surface snow samples collected monthly over a two-year period, at three distances from the Station (10, 500, and 1000 m). Even when the extracted total DNA was below the detection limit, 16S rRNA gene sequencing was successfully performed on all samples, while 18S rRNA was amplified on 19 samples out of 51. No significant relationships were observed between microbial diversity and seasonality (summer or winter) or distance from the Concordia base. This suggested that if present, the anthropogenic impact should have been below the detectable limit. While harboring low microbial diversity, the surface snow samples were characterized by heterogeneous microbiomes. Ultimately, our study corroborated the use of DNA sequencing-based techniques for revealing microbial presence in remote and hostile environments, with implications for Planetary Protection during space missions and for life-detection in astrobiology relevant targets.

1. Introduction

Extraterrestrial analogue sites on Earth are widely exploited to support planetary exploration by testing instruments and protocols for sampling and analysis procedures for robotic explorations as well as to train astronauts for future missions on Mars [1,2]. The French–Italian Concordia Research Station located on a plateau 3233 m above sea level is a confined, isolated habitat that provides a unique location for preparatory activities for future human journey to Mars, but also to explore microbial diversity in surface snow. With its permanent inhabitation, the Concordia Station is an ideal test site to monitor the dissemination of human-associated microorganisms within the natural (pristine) surrounding habitat and its indigenous microbial diversity. Unraveling the permanence of terrestrial microorganisms in near-surface frozen environments is also of importance to protect target planets from forward contamination by human and/or robotic exploration [3,4]. Indeed, planetary protection is a main concern for space agencies and governs for a responsible exploration of other worlds [5]. Conversely, investigating the limit of microbial life detectability in frozen habitats is relevant for future exploration missions of icy bodies in the Solar System such as Europa and Enceladus [6]. Currently, life-detection technologies are based on miniaturized mass spectrometry and Raman spectroscopy [7], although, other promising techniques are under development based on nanopore sequencing for detecting nucleic acid-based life [8], and an antibody-based technique [9].
A fundamental challenge to the investigation of microbial diversity in glaciers, namely snow, glacier surface, glacier ice, and basal ice, is given by cell abundancy that increases from snow to basal ice (for a review see, [10]). A bacterial content as low as 0.2–5 × 103 cells/mL was detected in the snowmelt of surface snow from the South Pole [11], while over 6 × 107 cells/mL occurred in the basal ice core from Greenland [12]. The lack of cultivation technologies that mimic natural conditions is the major obstacle to culturomics and, so far, only a few hundred bacterial strains have been isolated from glacial ice from Antarctica, Greenland, Canada, Tibet, and Alaska (for a review see, [10]). For instance, more bacteria were recovered from Late Holocene ice from the Taylor Dome region in Antarctica than from ice from the Antarctic peninsula or from Greenland [12,13], while viable microorganisms were identified after a gradual melting of Vostok ice samples [14,15].
The low biomass in glacier ice samples also impairs molecular approaches to unravel microbial diversity. A few studies based on rRNA amplicon sequencing were performed, like the one on Lake Vostok accretion ice [12,16], although recently, metagenome-assembled genomes were reported for Tibetan glaciers [17] and high-throughput sequencing was used for subglacial ice samples from Svalbard [18].
The present study was performed in the framework of the BacFinder project (European Space Agency, ESA AO-13-Concordia) with the aim to unravel the environmental- and human-associated microbial diversity in the surroundings of the Concordia Station. This is the first intensive and extensive surface snow sampling performed monthly over a two-year period at three distances (10, 500, and 1000 m) from the Concordia Station, and investigated by a high-throughput sequencing approach. Emphasis was laid on the relation between microbial presence and both seasonality and distance from the Base. The challenges of monthly sampling of Antarctic snow with low microbial biomass over a 2-year period generated certain study constraints; however, this study provides a more comprehensive view of microbial communities in surface snow of the Antarctic Polar Plateau. Our intensive study corroborated the use of DNA sequencing for revealing microbial presence in remote and hostile environments, with implications for Planetary Protections and for life-detection in astrobiology relevant targets.

2. Materials and Methods

2.1. Study Area

Ice samples were collected in the sampling area located at the Concordia Research Station (75°05′59″ S 123°19′56″ E) (Figure 1), which opened in 2005 as a French–Italian research facility that was built in a location called Dome C on the Antarctic Plateau, 3233 m above sea level. There, temperatures can drop to −80 °C in the winter, with a yearly average temperature of −50°C. Dome C does not experience the katabatic winds typical for the coastal regions of Antarctica because of its elevated location and its relative distance from the edges of the Antarctic Plateau. Typical wind speed in winter is 2.8 m/s. The nearest human beings are stationed some 600 km away at the Russian Vostok base, making Concordia more remote than the International Space Station. Concordia Base hosts up to 80 people during the Austral summer, including technicians, logistics, and researchers. During winter, a small group of 16 persons, named “Winter Over”, remain confined for at least 9 months.
During January–December 2015 and January–December 2016, two sampling campaigns were conducted monthly in the proximity of the Concordia station (10 m, L1), in the outer region of the research area (500 m, L2), and in the relative pristine-like area (in relative distance) to the Concordia area (1000 m, L3) as shown in Figure 1. Summer samples were collected in the period September–April; Winter samples were collected in the period May–August (see Supplementary Table S1 for more details on sampling). For each sampling site, ten replicates of snow in 50 mL-Falcons were sampled using clean, disposable nitrile gloves. Samples were preserved at −80 °C immediately upon collection and finally transported and stored at −80 °C at the German Aerospace Center (DLR, Cologne), until processing. During the processing, a total of 60 samples were yielded by pulling together 10 replicates for each site that yielded a total of 500 mL of snow for each sampling site.

2.2. Total DNA Extraction and Amplicon Sequencing

A total of 51 samples were processed for amplicon sequencing, while 9 (out of 60) samples were stored for culturomics experiments. Replicates were pooled resulting in a final volume of 500 mL for each sample. All samples were then melted in sterile conditions at 4 °C and filtered using Millipore filters (0.20 mm pore size). Total DNA was isolated from filters using the MOBIO Power Soil DNA Extraction kit (MOBIO Laboratories, Carlsbad, CA, USA) following the manufacturer’s instructions, and finally eluted in 50 µL. To account for possible contamination during the filtering and extraction steps, PCR-grade H2O was filtered through a Millipore filter and extracted following the same procedure as the snow samples.
The V4 region of the 16S rRNA and the 18S rRNA were targeted, respectively. The V4 region was amplified using the new developed barcoded F515 (GTGCCAGCMGCCGCGGTAA)/R806 (GGACTACHVGGGTWTCTAAT) primers as described by [19]; Euk_1391f (GTACACACCGCCCGTC) and EukBr (TGATCCTTCTGCAGGTTCACCTAC) primers were used to amplify the 18S rRNA region.
PCR was carried out with a total volume of 25 μL, containing 1 μL of each primer, 12.5 μL of BioMix (BioLine GmbH, Luckenwalde, Germany), 9.5 μL of nuclease-free water (Sigma–Aldrich, Gillingham, UK), and 5 ng of DNA template using an automated thermal cycler (Bio-Rad, Hercules, CA, USA). A major concern for low microbial biomass eDNA samples is the presence of contamination [20]. In addition to PCR controls, we performed PCR on extraction controls to identify any contaminants from the extraction kit. For bacteria, amplification for the V4 region followed a protocol of an initial denaturation at 94 °C for 3 min, 35 cycles of denaturation at 94 °C for 45 s, annealing at 50 °C for 1 min, extension at 72 °C for 90 s, followed by a final extension at 72 °C for 10 min. For eukaryotes, the 18S rRNA regions were amplified following initial denaturation at 94 °C for 1 min, 35 cycles of denaturation at 94 °C for 30 s, annealing at 52 °C for 30 s, extension at 68 °C for 90 s, followed by a final extension at 68 °C for 7 min. Amplicons, purified and quantified by Qubit dsDNA HS Assay Kit (Life Technologies, Carlsbad, CA, USA), were then barcoded and pooled, including the extraction and 16S rRNA amplification of PCR controls. Paired-end sequencing (2 × 300 bp) was carried out on the Illumina MiSeq platform at the Vincent J. Coates Genomics Sequencing Laboratory at the University of California (Berkeley, CA, USA), at the Edmund Mach Foundation (Trento, Italy), and at the Medical University of Graz, Center for Microbial Research (Graz, Austria).

2.3. Bioinformatics

The 18S rRNA and 16S rRNA amplicon sequencing datasets were processed with AMPtk: Amplicon ToolKit for Next Generation Sequence data (formally UFITS) (https://github.com/nextgenusfs/amptk accessed on 10 May 2022) v.1.0.0 [21], removing barcodes and primers from raw data. Briefly, reads were subjected to quality trimming to a maximum of 250 bp, discarding reads less than 100 bp in length, and chimera removal was performed utilizing USEARCH with default parameters v. 9.1.13 [22]. Sequence quality filtering was performed with the expected error parameter of 0.9 [23] and the cleaned dataset was clustered with UPARSE using a 97% percent identity parameter to generate the Operational Taxonomic Units (OTUs). Global singletons and rare taxa (<5 reads) were eliminated as likely false positives due to sequencing errors as described [24]. Finally, taxonomic identification was performed with hybrid database SINTAX/UTAX [22] as implemented in AMPtk. Negative controls (n = 3, average 718 reads) had substantially fewer reads than snow samples and these reads were taxonomically assigned to a common contaminant in extraction kits Ralstonia. Contaminating reads were removed from the 16S rRNA dataset before downstream analyses. Raw sequencing data have been archived in the NCBI SRA database linked to BioProject accession number PRJNA634088.

2.4. Downstream Analysis

Prior to further analysis, read counts were normalized; the number of reads in each sample using median sequencing depth (57,719 for 18S rRNA and 65,897 for 16S rRNA) and analyses were conducted on the rarefied dataset. The OTUs were further investigated by displaying bar plots of the relative abundance.
Statistical analyses were performed in R 2.15.1 (The R Foundation for Statistical Computing, Vienna, Austria) environment using the package phyloseq v.1.24.2 [25]. Several biodiversity indices such as richness in species, Shannon’s diversity [26], and Simpson’s dominance [27] were calculated. A pairwise comparison using the Wilcoxon rank sum test was performed to assess whether biodiversity indices differed (p < 0.05) among samples according to season sampling and distance from the research station.
Changes in alpha diversity and taxonomic composition were compared using one-way ANOVA and pairwise multiple comparison (Tukey test) using PAST v.2.17 software (Paleontological Statistics) [28]. A small probability p-value (<0.05) indicated a significant difference of community composition among all samples.
The effect of season sampling and distance was also tested using the permutational ANOVA (PERMANOVA) with 999 numbers of permutations applied with the “adonis” function in the “vegan” v.2.5-4 [29] package. Changes in community composition were displayed with PCoA ordination plots using the unweighted UniFrac distance.
Venn diagrams were constructed to show the number of 16S rRNA and 18S rRNA shared and unique OTUs, based on OTU tables normalized, using the VENNY tool [30].

3. Results

3.1. Raw Data Processing

DNA concentration was below the detection limit (0.005 ng/μL) for most of the samples. The 16S rRNA gene amplicon sequencing from 51 samples produced a total of 6,257,600 paired-end reads. Clustering of Operational Taxonomic Units (OTUs) was performed at 97% identity threshold, resulting in 5,187,083 reads mapped to OTUs. Sequences were grouped into 7171 OTUs (Supplementary Table S2). After singletons and rare taxa removal (9 out of 4536 OTUs), and normalization using median sequencing depth, a total of 2635 quality-filtered bacterial OTUs were obtained.
For eukaryotes, the 18S rRNA sequence reads were obtained from 19 samples (out of 54) and resulted in 1,305,012 sequences and upon filtering, 1,281,800 reads were mapped into 105 quality-filtered OTUs (Supplementary Table S3). The OTUs table was then normalized among samples using median sequencing depth.
In the 16S rRNA dataset, the rarefaction curve of total OTUs per sample reached a plateau only for a few samples, while the rarefaction curves generated from 18S rRNA datasets never reached a plateau (Supplementary Figure S1).

3.2. Taxonomy Structure and Composition

The majority of the identified 16S sequences recovered among all samples belonged to Proteobacteria (57%), followed by Cyanobacteria (15%), Actinobacteria (11%), Firmicutes (9%), Bacteroides (4%), and Acidobacteria (1%) as shown in Figure 2a. Members belonging to Planctomycetes, Verrucomicrobia, Gemmatimonadetes, Nitrospirae, Deinococcus-Thermus, Chloroflexi, Armatimonadetes, Fusobacteria, Parcubacteria, candidate division WPS-1, Chlamydiae, BRC1, Microgenomates, Spirochaetes, and Diapherotrites have been retrieved in the lowest percentage (<1%). Additionally, several archaeal taxa such as members belonging to Woesearchaeota, Thaumarchaeota, Pacearchaeota, and Euryarchaeota were also found, even if at lowest abundance (<1%).
As shown in Figure 2b, among the classes, Betaproteobacteria (30%) was the most abundant, followed by Gammaproteobacteria (18%), Actinobacteria (11%), Alphaproteobacteria (9%), Bacilli (9%), and Sphingobacteriia (2%). Members belonging to other classes were less frequent across the dataset (<1%), as follows: Acidobacteria (0.80%), Clostridia (0.69%), Flavobacteriia (0.61%), Acidobacteria Gp16 (0.50%), Cytophagia (0.41%), Cyanobacteria (0.40%), Planctomycetia (0.35%), Armatimonadia (0.26%), Fusobacteria (0.22%), Bacteroidia (0.19%), Gemmatimonadetes (0.14%), Acidobacteria Gp6 (0.13%), Spartobacteria (0.11%), and Deinococci (0.10%).
Instead, among eukaryotes, most of the identified sequences belonged to fungi, in particular the phylum Ascomycota that represents the most abundant phylum (66%) for fungi and occurred in all analyzed samples, followed by Basidiomycota (5%), while members belonging to Glomeromycota and Rozellomyconta were present in the lowest percentage (<1%). Algae belonging to Chlorophyta were represented by 28% (Figure 2c).
Among the fungal classes, the most representative were Lecanoromycetes (12.2%), Xylonomycetes (12.7%), Agaricomycetes (8.37%), Pezizomycetes (4.9%), Saccharomycetes (3.7%), Eurotiomycetes (3.1%), Dothideomycetes (2.9%), and Sordariomycetes (2.2%), while algae were represented by Trebouxiophyceae (49.5%) and Chlorophyceae (0.6%), respectively (Figure 2d).
We further examined the distribution of both 16S rRNA and 18S rRNA taxa, showing that their relative abundance, in general, was highly variable across all samples; indeed, most taxa were present intermittently among sampled sites (Supplementary Figures S2 and S3).
As OTUS were unidentified at a finer level, we proceeded with a comparison of taxonomic composition by distance and seasonality at phylum and class level. We found that changes in bacterial community composition at phylum level was not observed. Instead, at class level, Gammaproteobacteria and Alphaproteobacteria significantly (p < 0.05) varied according to season; in particular, Gammaproteobacteria were most abundant during winter, while the latter was in the summer (Figure 3a); moreover, Gammaproteobacteria were enriched (p < 0.05) in the proximity of the research station (L1, 10 m) (Figure 3b).
When 18S rRNA dataset was examined, we found that several classes significantly varied (p < 0.05) according to the season: in particular, Eurotiomycetes were particularly abundant in the winter as well as Xylonomycetes, Agaricomycetes, Chlorophyceae, and Dothideomycetes (Figure 4a). Instead, when comparing samples collected at different distances from the research station, we reported that Xylonomycetes were particularly abundant in L2 and L3, contrary to Lecanoromycetes that decreased in those two groups. In contrast, Dothideomycetes was predominant in L2, while Chlorophyceae were present mainly in L1 (p < 0.05) (Figure 4b).

3.3. Preliminary Data of Anthropogenic Effect on Alpha and Beta Diversity

The number of OTUs and biodiversity indices such as Shannon’s and Simpson’s metrics were calculated for each sample, as well as for 16S rRNA and 18S rRNA components. For prokaryotes, the number of observed OTUs ranged from 37 (sample S6) to 1507 (S12), while biodiversity varied from 1.95 (S18) to 4.93 (S16) and from 0.72 (S8) to 0.96 (S16 and S1) for Shannon’s and Simpson’s indices, respectively (Figure 5a,b; Supplementary Table S4).
Overall, the 18S rRNA biodiversity was lower than the 16S rRNA counterpart (Figure 5a,b). The richness ranged from 8 (sample 18S_S9) to 98 (18S_S12 and 18S_S14), Shannon’s index from 0.89 (18S_S17) to 2.49 (18S_S2), and Simpson’s index from 0.49 (18S_S17) to 0.86 (18S_S2) (Figure 5c,d; Supplementary Table S5).
We did not observe any significant change (Wilcoxon-Mann-Whitney) in 16S rRNA and 18S rRNA biodiversity according to site and sampling.
In both 16S rRNA and 18S rRNA data sets, the beta diversity was assessed from unweighted UniFrac distance and presented with PCoA plots and PERMANOVA statistical analysis. The effect of season sampling and distance were not found to be significant in 16S rRNA (p-value > 0.05), while the 18S rRNA samples were found to be significant (p-value = 0.028 < 0.05) for seasons only (Figure 6).

3.4. Unique and Shared Taxa According to Distances and Season

Venn diagrams were used to explore the percentage of shared and unique identified taxa at class level among all samples. As shown in Figure 7A, a substantial fraction of prokaryotic OTUs (1041 OTUs, 39.5%) were shared among the three sampling distances s, while 625 (23.7%), 181 (6.9%), and 114 (4.3%) OTUs were unique for L1, L2, and L3, respectively. When seasonality is considered, we found that almost 60% of prokaryotic taxa (1536 OTUs) were shared between winter and summer; further, a lower percentage of unique prokaryotes was found in the two seasons (400 OTUs, 15.2% in winter and 699 OTUs, 26.5% in summer, respectively) (Figure 7B).
For eukaryotes, Venn diagram showed that almost half (49.2%) of the total OTUs were shared among the three distances, while only 4 taxa (3%) were unique to L1 and L3 distances, including members of Sordariomycetes (Trichoderma sp.) and Chlorophyta, while 19 (14.9%) taxa were found exclusively in L2 (e.g., members belonging to Chytridiomycota, Pezizales, Aspergillus sp., and Trebouxiophyceae) (Figure 7C). Figure 7D showed that a few taxa (10 OTUs, 7.6%) were found in summer, exclusively; while a greater percentage of taxa (31 OTUs, 23.5%) were unique for the winter. The highest fraction of OTUs (91, 68.9%) was shared between the two seasons.

4. Discussion

Glacial habitats are among the most challenging natural environments for life. In particular, the ice sheet of the Antarctic Polar Plateau, the ancient, most isolated, stable and coldest icy environment on Earth, represents a peculiar case where liquid water essential for active life is virtually absent. This work contributes to deepening our knowledge on prokaryotic (V4 region of 16S rRNA gene) and eukaryotic (V4 region of 18S rRNA gene) microbial assemblages inhabiting the Antarctic ice shelf of the Polar Plateau, accounted as the most remote base on Earth, during a 2-year monthly sampling and over a distance of 1 km from the Concordia Research Station. The only microbiological study to date that is available on this habitat was performed several years ago and was based on a single snow sampling at 2 km distance from the Concordia Research Station and investigated only the bacterial 16S rRNA component. In that study, it was reported that the snow melt samples may contain approximately 3.4 × 102 cells mL−1 of autotrophic and heterotrophic bacteria and a negligible amount of eukaryotes [31].
In the present study, the PCR amplification and subsequent high throughput sequencing was performed on 51 samples for the 16S rRNA target, while only 19 samples were successfully targeted for the 18S rRNA amplification. Indeed, despite organisms occurring in an environment at a very low abundance, in the present study, microbial diversity could be revealed by PCR for Illumina sequencing of the almost undetectable (0.005 ng/µL) amounts of DNA that were obtained from each sampling site (500 mL of surface snow/site).
Among prokaryotes, members of phylum Proteobacteria resulted largely as the most abundant followed by Actinobacteria and Firmicutes. The first two, in particular, are reported among the most abundant prokaryotic representatives both in Earth’s major habitat types [32] and airborne bacteria in atmospheric bioaerosols [33,34]. Proteobacteria were also found to be dominant by Michaud and collaborators (2014) in a snow sample near Concordia Station and in other extreme cold environments [35,36,37,38,39]. Conversely, Cyanobacteria, despite being dominant in Antarctic terrestrial environments and having the capability to directly influence components of the cryosphere [40], were found only in a single snow sample (S7), at a 500-m distance from the Concordia Station. This very low incidence of autotrophic prokaryotes does not support the hypothesis of their role as primary producers and ecosystem engineers in an active microbial community in the Antarctic snow of the Polar Plateau as previously suggested [31].
Eukaryotes were found at very low abundance and were mostly represented by fungi. This is not surprising because most of these microorganisms produce abundant airborne propagules that are easily transported, even throughout remote environments. The majority belonged to the most abundant fungal phylum Ascomycota, which is by far the largest of the fungal Kingdom, hosting numerous abundantly sporulating fungi. Among the most abundant fungal classes found, Eurotiomycetes, whose members are of the Aspergillus/Penicillium group, are known to be globally distributed and copiously producing airborne propagules. The class Dothideomycetes is the widest in Ascomycota producing bitunicate asci; frequently encountered taxa in this class, such as Aureobasidium, Cladosporium, or Hortaea have common characteristics of adaptation or resistance to low temperature and high osmotic pressure, which may be essential to survival under the icy conditions of the Polar Plateau. Other common taxa in this class were, instead, typically plant associated as for the genus Phoma.
Surprisingly, we found members of the class Xylonomycetes, only relatively recently discovered and described, which are exclusively found associated as endophytes with plants of remote forests in Peru [41]. Despite the fact that they also are sporulating and producing picnidia even under lab conditions, they are not widespread fungi and their random presence in ice samples of the Polar Plateau is difficult to explain. Their presence may be more probably ascribed to their capacity to spread globally since it is becoming evident that much of the microbial life within remote environments, including Antarctica, is transported by air currents [42]. This is true even for Chlorophyceae algae, here found in a relatively high abundance, which are reported to be easily transported overseas among continents [43].
Our interpretation of the coincidental occurrence of microbial eukaryotes found was also corroborated by the almost undetectable presence and lack of any correlation with distance to the Concordia Station. Even their apparently higher presence in the winter season, when the anthropogenic impact is surely much lower, is difficult to explain, if not as an accidental datum due to the high variability among samples and the very insubstantial presence of DNA.
OTUs richness and other biodiversity indices such as Shannon’s and Simpson’s metrics were calculated for each sample for comparative purposes to test whether different conditions (i.e., distance from the station and seasonality) significantly impact community composition. Overall, this work suggested that while harboring low microbial diversity, the investigated snow samples were characterized by heterogeneous microbiomes. Indeed, a considerable amount of heterogeneity in alpha-diversity has been observed between samples, both for prokaryotes and eukaryotes. Specifically, in our survey, richness and Shannon’s index values showed relatively low diversity; they significantly varied among all samples: number of observed species ranged from 37 to 1507 and 8 to 98 in fungi and bacteria, respectively, while Shannon’s diversity values ranged from 1.95 to 4.93 and 0.89 to 2.49 in the two assemblages, respectively. However, this high variability was correlated to neither distance from the station nor sampling seasons (p > 0.05). Interestingly, snow samples were found to have a core (i.e., OTUs present in at least 75% of the samples) microbiome which persisted regardless of the environmental factors and level of human activity.
Taken together, this preliminary study provided a high-throughput sequencing of the prokaryotic and eukaryotic icy communities occurring in the Antarctic plateau in the vicinity of the Concordia Research Station. In order to provide new insights into the anthropogenic impact on these fragile ecosystems, the taxonomy and community composition of the surface snow samples were evaluated by considering the distance from the Concordia Research Station and the seasonality of the sampling.
A further and more extensive sampling is necessary to better investigate the microbial diversity in superficial snow over a prolonged period of time. Therefore, in order to increase the biomass to obtain adequate amounts of DNA to enable a comprehensive meta-community profiling, the crew members of the Concordia Research Station collected during the 2018–2019 campaign 15 samples of 50 mL-Falcon, instead of 10 samples of 50 mL-Falcon as in the 2015–2016 campaign.
Data suggested that if present, the anthropogenic impact should have been below the detectable limit. Snow samples had variable microbial assemblages that did not associate with changes in seasonality or distance from station, suggesting a transient existence for resident microbes that might not be significantly influenced by environmental variables. Furthermore, RNA analysis associated with a metagenomic approach could also provide insights into the metabolic status of the microbial communities.
Finally, our study corroborated the use of DNA sequencing-based techniques for revealing microbial presence in remote and hostile environments, with implications for Planetary Protection during space missions and for the development of life–life-detection instrumentations for astrobiological relevant targets such as Mars and the icy moons Europa and Enceladus.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life13010113/s1, Figure S1: Rarefaction curve of total OTUs, 16S rRNA (A) and 18S rRNA (B) datasets; Figure S2: Barplot of Class in 16S rRNA dataset, Barplot based on season, as shown in Supplementary Table S1: Summer and Winter (A), Barplot based on distance (L1 = 10 m, L2 = 500 m and L3 = 1000 m) from Concordia Research Station (B); Figure S3: Barplot of Class 18S rRNA dataset. Barplot based on season, as shown in Supplementary Table S1: Summer and Winter (A), Barplot based on distance (L1 = 10 m, L2 = 500 m and L3 = 1000 m) from Concordia Research Station (B); Table S1: Distances at Concordia Station and metadata. Data include samples correspondence, distance from research station, coordinates, season and sampling date are reported. The samples for which 16S rRNA and 18S rRNA amplicon sequencing was successful are reported in green, while the unsuccessfully ones are in in red; Table S2: 16S rRNA taxonomy. Taxonomy assignment (results at 97% of identity) was performed with SINTAX/UTAX databases (Edgar, 2010); Table S3: 18S rRNA taxonomy. Taxonomy assignment (results at 97% of identity) was performed with SINTAX/UTAX databases (Edgar, 2010); Table S4: 16S rRNA biodiversity. Observed OTUs Shannon’s and Simpson’s indices for each sample; Table S5: 18S rRNA biodiversity. Observed OTUs Shannon’s and Simpson’s indices for each sample.

Author Contributions

Conceptualization, L.S., D.B., and R.M.; methodology, A.N., C.C., and N.J.U.; formal analysis, A.N., C.C., and N.J.U.; funding acquisition, D.B., L.S., and R.M.; writing—original draft, C.C., L.S., and D.B.; writing—review and editing, R.M., L.S., D.B., A.N., C.C., and N.J.U. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Antarctic Research Program (MIcrobial Diversity within the Vicinity of the Concordia Antarctic Station _MIDAS project grant PNRA 2016_00101 to DB). N.J.U. and R.M were supported by the DLR grant FuE-Projekt “ISS LIFE” (Programm RF-FuW, TP 475). The BacFinder project was supported by a Europlanet 2020 RI grant from the European Union’s Horizon 2020 research and innovation program under grant agreement No. 654208.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

It does not involve humans.

Data Availability Statement

Data available in a publicly accessible repository. NCBI SRA database linked to BioProject accession number PRJNA634088.

Acknowledgments

We extend thanks to the ESA scientists for sample collection. The authors express their gratitude to Christine Moissl-Eichinger for help with 16S rRNA sequencing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Foucher, F.; Hickman-Lewis, K.; Hutzler, A.; Joy, K.H.; Folco, L.; Bridges, J.C.; Wozniakiewicz, P.; Martínez-Frías, J.; Debaille, V.; Zolensky, M.; et al. Definition and use of functional analogues in planetary exploration. Planet. Space Sci. 2021, 197, 105162. [Google Scholar] [CrossRef]
  2. Preston, L.J.; Dartnell, L.R. Planetary habitability: Lessons learned from terrestrial analogues. Int. J. Astrobiol. 2014, 13, 81–98. [Google Scholar] [CrossRef] [Green Version]
  3. Marion, G.M.; Fritsen, C.H.; Eicken, H.; Payne, M.C. The search for life on Europa: Limiting environmental factors, potential habitats, and Earth analogues. Astrobiology 2003, 3, 785–811. [Google Scholar] [CrossRef] [PubMed]
  4. Nicholson, W.L.; Schuerger, A.C.; Race, M.S. Migrating microbes and planetary protection. Trends Microbiol. 2009, 17, 389–392. [Google Scholar] [CrossRef] [PubMed]
  5. Cheney, T.; Newman, C.; Olsson-Francis, K.; Steele, S.; Pearson, V.; Lee, S. Planetary protection in the new space era: Science and Governance. Front. Astron. Space Sci. 2020, 7, 90. [Google Scholar] [CrossRef]
  6. Martins, Z.; Cottin, H.; Kotler, J.M.; Carrasco, N.; Cockell, C.S.; de la Torre Noetzel, R.; Demets, R.; de Vera, J.P.; d’Hendecourt, L.; Ehrenfreund, P.; et al. Earth as a tool for astrobiology—A European perspective. Space Sci. Rev. 2017, 209, 43–81. [Google Scholar] [CrossRef] [Green Version]
  7. Marshall, C.P.; Edwards, H.G.; Jehlicka, J. Understanding the application of Raman spectroscopy to the detection of traces of life. Astrobiology 2010, 10, 229–243. [Google Scholar] [CrossRef] [Green Version]
  8. Arevalo, R., Jr.; Ni, Z.; Danell, R.M. Mass spectrometry and planetary exploration: A brief review and future projection. J. Mass. Spectrom. 2019, 55, e4454. [Google Scholar]
  9. Carr, C.E.; Bryan, N.C.; Saboda, K.N.; Bhattaru, S.A.; Ruvkun, G.; Zuber, M.T. Nanopore sequencing at Mars, Europa, and microgravity conditions. NPJ Micro. 2020, 6, 24. [Google Scholar] [CrossRef]
  10. Fairén, A.G.; Gómez-Elvira, J.; Briones, C.; Prieto-Ballesteros, O.; Rodríguez-Manfredi, J.A.; López Heredero, R.; Belenguer, T.; Moral, A.G.; Moreno-Paz, M.; Parro, V. The complex molecules detector (CMOLD): A fluidic-based instrument suite to search for (bio)chemical complexity on Mars and icy moons. Astrobiology 2020, 20, 1076–1096. [Google Scholar] [CrossRef]
  11. Kim, S.; Lee, H.; Hur, S.D.; Sul, W.J.; Kim, O.S. Glaciers as microbial habitats: Current knowledge and implication. J. Microbiol. 2022, 60, 67–779. [Google Scholar] [CrossRef] [PubMed]
  12. Carpenter, E.J.; Lin, S.; Capone, D.G. Bacterial activity in South Pole snow. Appl. Environ. Microbiol. 2000, 66, 4514–4517. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sheridan, P.P.; Miteva, V.I.; Brenchley, J.E. Phylogenetic analysis of anaerobic psychrophilic enrichment cultures obtained from a Greenland glacier ice core. Appl. Environ. Microbiol. 2003, 69, 2153–2160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Christner, B.C.; Mosley-Thompson, E.; Thompson, L.G.; Zagorodnov, V.; Sandman, K.; Reeve, J.N. Recovery and identification of viable bacteria immured in glacial ice. Icarus 2000, 144, 479–485. [Google Scholar] [CrossRef]
  15. Karl, D.M.; Bird, D.F.; Björkman, K.; Houlihan, T.; Shackelford, R.; Tupas, L. Microorganisms in the accreted ice of Lake Vostok. Science 1999, 286, 2144–2147. [Google Scholar] [CrossRef] [Green Version]
  16. Abyzov, S.S.; Mitskevich, I.N.; Poglazova, M.N.; Barkov, N.I.; Lipenkov, V.Y.; Bobin, N.E.; Koudryashov, B.B.; Pashkevich, V.M.; Ivanov, M.V. Microflora in the basal strata at the Antarctic ice core above the Vostok Lake. Adv. Space Res. 2001, 28, 701–706. [Google Scholar] [CrossRef]
  17. Liu, Y.; Ji, M.; Yu, T.; Zaugg, J.; Anesio, A.M.; Zhang, Z.; Hu, S.; Hugenholtz, P.; Liu, K.; Liu, P.; et al. A genome and gene catalog of glacier microbiomes. Nat. Biotechnol. 2022, 40, 1341–1348. [Google Scholar] [CrossRef]
  18. Perini, L.; Gostinčar, C.; Gunde-Cimerman, N. Fungal and bacterial diversity of Svalbard subglacial ice. Sci. Rep. 2019, 9, 20230. [Google Scholar] [CrossRef] [Green Version]
  19. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Huntley, J.; Fierer, N.; Owens, S.M.; Betley, J.; Fraswer, L.; Bauer, M.; et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012, 6, 1621–1624. [Google Scholar] [CrossRef] [Green Version]
  20. Hornung, B.V.H.; Zwittink, R.D.; Kuijper, E.J. Issues and current standards of controls in microbiome research. FEMS Microbiol. Ecol. 2019, 95, fiz045. [Google Scholar] [CrossRef] [Green Version]
  21. Palmer, J.M.; Jusino, M.A.; Banik, M.T.; Lindner, D.L. Non-biological synthetic spike-in controls and the AMPtk software pipeline improve mycobiome data. PeerJ 2018, 6, e4925. [Google Scholar] [CrossRef] [PubMed]
  22. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Edgar, R.C.; Flyvbjerg, H. Error filtering, pair assembly and error correction for next-generation sequencing reads. Bioinformatics 2015, 31, 3476–3482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lindahl, B.D.; Nilsson, R.H.; Tedersoo, L.; Abarenkov, K.; Carlsen, T.; Kjøller, R.; Kõljalg, U.; Pennanen, T.; Rosendahl, S.; Stenlid, J.; et al. Fungal community analysis by high-throughput sequencing of amplified markers—A user’s guide. New Phytol. 2013, 199, 288–299. [Google Scholar] [CrossRef] [Green Version]
  25. McMurdie, P.J.; Holmes, S. Phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [Green Version]
  26. Shannon, C.E.; Weaver, W. The Mathematical Theory of Communication; University of Illinois Press: Champaign, IL, USA, 1963. [Google Scholar]
  27. Simpson, E.H. Measurement of diversity. Nature 1949, 163, 688. [Google Scholar] [CrossRef]
  28. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST-palaeontological statistics. Software package for education and data analysis. Palaeontol. Electr. 2001, 4, 1–9. [Google Scholar]
  29. Dixon, P. VEGAN, a package of R functions for community ecology. J. Veg. Sci. 2003, 14, 927–930. [Google Scholar] [CrossRef]
  30. Oliveros, J.C. VENNY. An Interactive Tool for Comparing Lists with Venn Diagrams. 2007. Available online: http://bioinfogp.cnb.csic.es/tools/venny/index.html (accessed on 10 May 2022).
  31. Michaud, L.; Giudice, A.L.; Mysara, M.; Monsieurs, P.; Raffa, C.; Leys, N.; Amalfitano, S.; Van Houdt, R. Snow surface microbiome on the high Antarctic plateau (DOME C). PLoS ONE 2014, 9, e104505. [Google Scholar] [CrossRef]
  32. Nemergut, D.R.; Costello, E.K.; Hamady, M.; Lozupone, C.; Jiang, L.; Schmidt, S.K.; Fierer, N.; Townsend, A.R.; Cleveland, C.C.; Stanish, L.; et al. Global patterns in the biogeography of bacterial taxa. Environ. Microbiol. 2011, 13, 135–144. [Google Scholar] [CrossRef]
  33. Tang, K.; Huang, Z.; Huang, J.; Maki, T.; Zhang, S.; Shimizu, A.; Ma, X.; Shi, J.; Bi, J.; Zhou, T.; et al. Characterization of atmospheric bioaerosols along the transport pathway of Asian dust during the dust-bioaerosol 2016 Campaign. Atmos. Chem. Phys. 2018, 18, 7131–7148. [Google Scholar] [CrossRef] [Green Version]
  34. Fang, J.; Dong, Q.; Shen, W.; Liu, X.; Dou, N.; Xian, L.; Chen, H. Variation in near-surface airborne bacterial communities among five forest types. Forests 2020, 11, 561. [Google Scholar] [CrossRef]
  35. Shivaji, S.; Begum, Z.; Shiva Nageswara Rao, S.S.; Reddy, P.V.V.V.; Manasa, P.; Sailaja, B.; Prathiba, M.S.; Thamban, M.; Krishnan, K.P.; Singh, S.M.; et al. Antarctic ice core samples: Culturable bacterial diversity. Res. Microbiol. 2013, 164, 70–82. [Google Scholar] [CrossRef] [PubMed]
  36. Canini, F.; Zucconi, L.; Coleine, C.; D′Alò, F.; Onofri, S.; Geml, J. Expansion of shrubs could result in local loss of soil bacterial richness in Western Greenland. FEMS Microbiol. Ecol. 2020, 96, fiaa089. [Google Scholar] [CrossRef] [PubMed]
  37. Ertekin, E.; Meslier, V.; Browning, A.; Treadgold, J.; DiRuggiero, J. Rock structure drives the taxonomic and functional diversity of endolithic microbial communities in extreme environments. Environ. Microbiol. 2021, 23, 3937–3956. [Google Scholar] [CrossRef] [PubMed]
  38. Albanese, D.; Coleine, C.; Rota-Stabelli, O.; Onofri, S.; Tringe, S.G.; Stajich, J.E.; Selbmann, L.; Donati, C. Pre-Cambrian roots of novel Antarctic cryptoendolithic bacterial lineages. Microbiome 2021, 9, 63. [Google Scholar] [CrossRef] [PubMed]
  39. Coleine, C.; Biagioli, F.; de Vera, J.P.; Onofri, S.; Selbmann, L. Endolithic microbial composition in Helliwell Hills, a newly investigated Mars-like area in Antarctica. Environ. Microbiol. 2021, 23, 4002–4016. [Google Scholar] [CrossRef]
  40. Van Goethem, M.W.; Cowan, D.A. Role of cyanobacteria in the ecology of polar environments. In The Ecological Role of Micro-Organisms in the Antarctic Environment; Castro-Sowinski, S., Ed.; Springer: Cham, Switzerland, 2019; pp. 3–23. [Google Scholar]
  41. Gazis, R.; Miadlikowska, J.; Lutzoni, F.; Arnold, A.E.; Chaverri, P. Culture-based study of endophytes associated with rubber trees in Peru reveals a new class of Pezizomycotina: Xylonomycetes. Mol. Phylogenet. Evol. 2012, 65, 294–304. [Google Scholar] [CrossRef]
  42. Pearce, D.A.; Bridge, P.D.; Hughes, K.A.; Sattler, B.; Psenner, R.; Russell, N.J. Microorganisms in the atmosphere over Antarctica. FEMS Microbiol. Ecol. 2009, 69, 143–157. [Google Scholar] [CrossRef]
  43. Elster, J. Algal versatility in various extreme environments. In Enigmatic Microorganisms and Life in Extreme Environments; Seckbach, J., Ed.; Springer: Dordrecht, The Netherlands, 1999; pp. 215–227. [Google Scholar]
Life 13 00113 g001 550
Figure 1. Concordia Research Station and sampling sites. (A) Map of Antarctica showing the Concordia Research Station at Dome C (B); sampling site L1 at 10 m (C); sampling site L2 at 500 m (D); and sampling site L3 at 1000 m (E). (A–E) photo credit: European Space Agency.
Figure 1. Concordia Research Station and sampling sites. (A) Map of Antarctica showing the Concordia Research Station at Dome C (B); sampling site L1 at 10 m (C); sampling site L2 at 500 m (D); and sampling site L3 at 1000 m (E). (A–E) photo credit: European Space Agency.
Life 13 00113 g001
Life 13 00113 g002 550
Figure 2. Microbial community composition. Relative abundance of top 10 identified Phyla and Classes in 16S (a,b) and 18S (c,d) datasets.
Figure 2. Microbial community composition. Relative abundance of top 10 identified Phyla and Classes in 16S (a,b) and 18S (c,d) datasets.
Life 13 00113 g002
Life 13 00113 g003 550
Figure 3. Barplot of Class relative abundances in 16S rRNA dataset. Barplot based on the sampling season as shown in supplementary Table S1: Summer and Winter (a). Barplot based on sampling distances from Concordia Research Station: L1 = 10 m, L2 = 500 m, and L3 = 1.000 m (b).
Figure 3. Barplot of Class relative abundances in 16S rRNA dataset. Barplot based on the sampling season as shown in supplementary Table S1: Summer and Winter (a). Barplot based on sampling distances from Concordia Research Station: L1 = 10 m, L2 = 500 m, and L3 = 1.000 m (b).
Life 13 00113 g003
Life 13 00113 g004 550
Figure 4. Barplot of Classes relative abundances in 18S rRNA dataset. Barplot based on the sampling season, as shown in Table S1: Summer and Winter (a). Barplot based on sampling distances from Concordia Research Station: L1 = 10 m, L2 = 500 m, and L3 = 1.000 m (b).
Figure 4. Barplot of Classes relative abundances in 18S rRNA dataset. Barplot based on the sampling season, as shown in Table S1: Summer and Winter (a). Barplot based on sampling distances from Concordia Research Station: L1 = 10 m, L2 = 500 m, and L3 = 1.000 m (b).
Life 13 00113 g004
Life 13 00113 g005 550
Figure 5. Biodiversity indices. Observed OTUs, Shannon’s and Simpson’s indices for 16S, based on season (a), distances (b), and 18rRNA datasets (c,d). Boxplots show 25th and 75th percentile, while error bars 1st and 99th percentile.
Figure 5. Biodiversity indices. Observed OTUs, Shannon’s and Simpson’s indices for 16S, based on season (a), distances (b), and 18rRNA datasets (c,d). Boxplots show 25th and 75th percentile, while error bars 1st and 99th percentile.
Life 13 00113 g005
Life 13 00113 g006 550
Figure 6. Principal Coordinate Analysis (PCoA) plot of the beta diversity distances obtained from the unweighted UNIFRAC analysis of Observation, Shannon index, and Simpson index. PCoA of 16S rRNA sequences based on season (a), distances (b), and PCoA of 18S rRNA sequences based on season (c) and distances (d).
Figure 6. Principal Coordinate Analysis (PCoA) plot of the beta diversity distances obtained from the unweighted UNIFRAC analysis of Observation, Shannon index, and Simpson index. PCoA of 16S rRNA sequences based on season (a), distances (b), and PCoA of 18S rRNA sequences based on season (c) and distances (d).
Life 13 00113 g006
Life 13 00113 g007 550
Figure 7. Venn diagram showing the distribution of bacterial taxa (16S rRNA dataset) relative to distances (A) and season (B) data. Venn diagrams showing the distribution of eukaryotic taxa (18S rRNA taxa dataset) on distances and season data, respectively (C,D). Percentages of both shared and unique OTUs are shown in parentheses.
Figure 7. Venn diagram showing the distribution of bacterial taxa (16S rRNA dataset) relative to distances (A) and season (B) data. Venn diagrams showing the distribution of eukaryotic taxa (18S rRNA taxa dataset) on distances and season data, respectively (C,D). Percentages of both shared and unique OTUs are shown in parentheses.
Life 13 00113 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Napoli, A.; Coleine, C.; Ulrich, N.J.; Moeller, R.; Billi, D.; Selbmann, L. Snow Surface Microbial Diversity at the Detection Limit within the Vicinity of the Concordia Station, Antarctica. Life 2023, 13, 113. https://doi.org/10.3390/life13010113

AMA Style

Napoli A, Coleine C, Ulrich NJ, Moeller R, Billi D, Selbmann L. Snow Surface Microbial Diversity at the Detection Limit within the Vicinity of the Concordia Station, Antarctica. Life. 2023; 13(1):113. https://doi.org/10.3390/life13010113

Chicago/Turabian Style

Napoli, Alessandro, Claudia Coleine, Nikea J. Ulrich, Ralf Moeller, Daniela Billi, and Laura Selbmann. 2023. "Snow Surface Microbial Diversity at the Detection Limit within the Vicinity of the Concordia Station, Antarctica" Life 13, no. 1: 113. https://doi.org/10.3390/life13010113

APA Style

Napoli, A., Coleine, C., Ulrich, N. J., Moeller, R., Billi, D., & Selbmann, L. (2023). Snow Surface Microbial Diversity at the Detection Limit within the Vicinity of the Concordia Station, Antarctica. Life, 13(1), 113. https://doi.org/10.3390/life13010113

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop