Matthias Winkel1*†‡
Christopher B. Trivedi1‡
Rey Mourot1
James A. Bradley1,2Andrea Vieth-Hillebrand1
Liane G. Benning1,3*- 1GFZ German Research Centre for Geosciences, Helmholtz Centre for Geosciences, Potsdam, Germany
- 2School of Geography, Queen Mary University of London, London, United Kingdom
- 3Department of Earth Sciences, Freie Universität Berlin, Berlin, Germany
Blooms of microalgae on glaciers and ice sheets are amplifying surface ice melting rates, which are already affected by climate change. Most studies on glacial microorganisms (including snow and glacier ice algae) have so far focused on the spring and summer melt season, leading to a temporal bias, and a knowledge gap in our understanding of the variations in microbial diversity, productivity, and physiology on glacier surfaces year-round. Here, we investigated the microbial communities from Icelandic glacier surface snow and bare ice habitats, with sampling spanning two consecutive years and carried out in both winter and two summer seasons. We evaluated the seasonal differences in microbial community composition using Illumina sequencing of the 16S rRNA, 18S rRNA, and ITS marker genes and correlating them with geochemical signals in the snow and ice. During summer, Chloromonas, Chlainomonas, Raphidonema, and Hydrurus dominated surface snow algal communities, while Ancylonema and Mesotaenium dominated the surface bare ice habitats. In winter, algae could not be detected, and the community composition was dominated by bacteria and fungi. The dominant bacterial taxa found in both winter and summer samples were Bacteriodetes, Actinobacteria, Alphaproteobacteria, and Gammaproteobacteria. The winter bacterial communities showed high similarities to airborne and fresh snow bacteria reported in other studies. This points toward the importance of dry and wet deposition as a wintertime source of microorganisms to the glacier surface. Winter samples were also richer in nutrients than summer samples, except for dissolved organic carbon—which was highest in summer snow and ice samples with blooming microalgae, suggesting that nutrients are accumulated during winter but primarily used by the microbial communities in the summer. Overall, our study shows that glacial snow and ice microbial communities are highly variable on a seasonal basis.
Introduction
Glaciers and ice sheets cover roughly 10% of Earth’s land surface and store approximately 69% of Earth’s freshwater (Meredith et al., 2019), thus serving as a critical natural resource. Moreover, glaciers and ice sheets are important components of Earth’s climate system—responding to global temperature changes, and are a major driver of sea level change (Meier et al., 2007). There has been increased surface melting of glaciers and ice sheets in recent decades (Nghiem et al., 2012; Lenaerts et al., 2013) resulting from both land surface air temperature increases as well as reductions in glacier surface albedo (Box et al., 2012; Ji et al., 2014; Cook et al., 2020; Tedstone et al., 2020)—a change to which the biodiversity of glacier and glacier adjacent habitats is sensitive (Anesio et al., 2017; Cauvy-Fraunié and Dangles, 2019; Stibal et al., 2020). Pigmented algae inhabiting glacier and ice sheet surfaces are part of the light absorbing particulates present on snow and ice surfaces. Such particulates have been shown to decrease surface albedo and increase surface melt rates of both snow and bare ice surfaces (Benning et al., 2014; Dumont et al., 2014; Lutz et al., 2016; Cook et al., 2020; Williams et al., 2020). Among these particulates, pigmented snow and glacier ice algal blooms play a major role in changing the albedo during the late spring to late summer melt season. In the dark zone on the western margin of the Greenland Ice Sheet, glacier ice algae contribute between 13% and 26% of the total annual surface melt (Cook et al., 2020) by lowering the albedo. As the main phototrophic primary producers on glacier surfaces, algae are also important for carbon and nutrient cycling (Anesio et al., 2017) as they produce dissolved organic carbon (DOC; Hood et al., 2009, 2015; Singer et al., 2012) that acts as a substrate for heterotrophic organisms (Stibal et al., 2012) and fertilizes downstream ecosystems when exported in meltwater (Lawson et al., 2014; Musilova et al., 2017; Irvine-Fynn et al., 2021).
Algae are found in numerous different habitats on glaciers, including cryoconite holes (Hamilton et al., 2013; Lutz et al., 2019), snow (Lutz et al., 2015a; Havig and Hamilton, 2019), bare ice (Yallop et al., 2012; Lutz et al., 2016, 2018; Williamson et al., 2018), biofilms (Lutz et al., 2017), and glacial streams (Smith et al., 2017). Such pigmented algae are often visible in distinct red or deep purple blooms on the snow and glacial ice surfaces (Lutz et al., 2016, 2018; Cook et al., 2020; Stibal et al., 2020; Williamson et al., 2020). The most commonly reported green and red snow algae are those of the phyla Chlorophyta and Ochrophyta, and the genera Chloromonas, Chlamydomonas, and Raphidonema (Lutz et al., 2015a), while Hydrurus has been reported as the dominant genus in yellow snow (Remias et al., 2013). In contrast, algal communities on bare ice surfaces are typically dominated by the phylum Streptophyta (Zygnematophyceae) consisting primarily of Mesotaenium (recently reclassified as Ancylonema alaskana; Procházková et al., 2021) and Ancylonema nordenskioeldii (Remias et al., 2009; Yallop et al., 2012; Lutz et al., 2018). Bacteria, archaea, fungi, and viruses are also found in glacial surface habitats (Larose et al., 2010; Maccario et al., 2014; Bellas et al., 2015; Lutz et al., 2015b; Perini et al., 2019). Most of the bacteria and fungi inhabiting glacier surfaces are heterotrophs that can feed on algal exudates (Lutz et al., 2015b), and they drive important biogeochemical processes including the fixation of atmospheric nitrogen (Telling et al., 2011; Larose et al., 2013).
Studies of microbial processes in supraglacial ecosystem have been limited to the spring and summer melt seasons, and an evaluation of variation in summer vs. winter dynamics is lacking. Challenges associated with sampling glacial and high-latitude environments during winter conditions include inaccessibility of sampling sites, adverse weather and light conditions, and low rates of biological activity (Carpenter et al., 2000). Furthermore, deep snow tends to cover glacier surface during winter, making the ice surface less accessible for sampling. Bacterial communities in snowpacks have been well characterized, even in early spring prior to melt (Larose et al., 2010; Maccario et al., 2014), however knowledge of the seasonal dynamics of bacterial, archaeal, fungal, viral, or algal communities is poor. Microbial communities have been analyzed seasonally in other cryospheric environments, including the active layers of permafrost (Schostag et al., 2015), thermokarst lake water columns (Vigneron et al., 2019), and alpine soils (Lazzaro et al., 2015), yet these systems differ significantly from glacial ice and snow and therefore offer little to no comparison of biological, chemical and physical processes.
Here, we assessed the microbial community composition and nutrient variations of glacier surface ice and snow habitats on various Icelandic glaciers, with sampling carried out during two consecutive summers and one winter time period. We analyzed the structure of glacial microbial communities using bacterial and eukaryotic marker genes and measured physico-chemical parameters including the concentrations of DOC, major ions, pH, and conductivity. For eukaryotic communities, we used two marker genes (18S rRNA and ITS) as previous studies (Lutz et al., 2018; Engstrom et al., 2020) showed these to give a more complete overview of eukaryotic diversity. Our data revealed distinct differences between summer and winter supraglacial microbial communities and we discuss their links to variable geochemical signatures of the surface snow and ice environments.
Materials and Methods
Field Sites, Sampling, and Sample Handling
Snow and ice samples were collected from four glaciers and one ice cap in Iceland (Figure 1) during the summers in 2018 (September) and 2019 (August), and in the winter 2019 (February; Table 1). Not all locations were sampled during each time point and snow/ice sample pairs could not always be collected at each site. In particular, during the 2019 winter sampling at Snæfellsjökull and Langjökull, only fresh snow could be collected. In summer, additional samples representing snow–ice interfaces were collected at Langjökull (2018) and Snæfellsjökull (2019), and pro-glacial waters were collected only at Langjökull (2018). The top few cm of 10 snow, seven ice, two snow–ice interface habitats, and one pro-glacial water sample were collected into sterile whirl pack bags or into 50 ml sterile Falcon® tubes using pre-sterilized plastic or metal trowels. The samples were used for marker gene sequencing, solution and particulate analysis, pH, and conductivity. For DOC analysis, separate samples were collected with sterile metal trowels into acid-washed (~10% HCl, 24 h) and ashed (550°C, 6 h) glass jars.
Figure 1. Sampling collection sites and habitats in Iceland. (A) Map of Iceland showing the geographical locations and time of season of the glaciers sampled for this study. (B) Images of representative habitats collected in the different seasons (colors; ice summer 2018 © J. Perez, snow winter 2018 © C. Trivedi, snow–ice interface, snow summer 2019 © J. Bradley, pro-glacial water summer 2018, ice winter 2019 © M. Winkel). The scale bar on all photographs represents 12 cm.
Table 1. Overview of samples, locations, habitats, coordinates, and field measurements.
All samples were returned to the field laboratory, thawed at room temperature, and immediately processed once fully melted. For dissolved inorganic ions, sample aliquots were filtered through 0.2 μm polycarbonate filters into either acid-washed Nalgene® bottles (for cations, acidified upon return to the home institution with ultra-pure nitric acid), or into 15 ml Falcon® tubes (for anions). Samples for DOC analysis were filtered using acid-cleaned glass syringes, 0.7 μm ashed glass fiber filters (GFF), and acid-washed metal filtration units, and the solutions were filtered directly into acid-washed, ashed 40 ml amber vials sealed with an acid-cleaned Teflon septa. For solid particulate bulk carbon, nitrogen, and isotope analysis, we filtered between 0.5 and ~20 L of melted snow and ice through 0.7 μm ashed GFF filters. Resulting solids were transferred with ethanol-sterilized flat forceps into ashed glass vials. Upon return to the home laboratory, these solids were freeze-dried and milled with a ball mill (Retsch MM2000). For DNA analysis, between 0.05 and ~12 L (exact volumes in Table 1) of melted snow or ice, depending on season and particle load, were filtered through single use, sterile, 0.2 μm cellulose nitrate filters (Thermo Scientific Nalgene). Filters were removed from the units, folded with sterile forceps, and transferred, into sterile 5 ml cryotubes that were returned to the home laboratory in a cryo-shipper (filled with liquid nitrogen) and subsequently stored at −80°C, until extraction and analysis. A separate 15 ml Falcon® tube was filled with a sample preserved with glutaraldehyde (final concentration 2%) and returned to the home laboratory for microscopic observations.
Dissolved and Particulates Analyses
The concentrations of dissolved inorganic major, minor, and trace elements in snow and ice samples were analyzed by inductively coupled plasma mass spectrometry (ICP-MS; Thermo Fisher iCAPQc) following methods described in McCutcheon et al. (2021). The precision of the analyses varied between 2% and 5%. Inorganic anions were analyzed by ion chromatography (IC) using conductivity detection (ICS 3000, Dionex). DOC was analyzed by LC-OCD (liquid chromatography—organic carbon detector) where organic carbon was quantified by IR detection of released CO2 after UV photo-oxidation (185 nm) in a Gräntzel thin-film reactor. Detailed descriptions of the latter two methods have been published elsewhere (Regenspurg et al., 2018).
Fine-milled particulate samples were analyzed for total carbon (TC), total nitrogen (TN), and carbon isotopic composition (δ13C) of TC using a Carlo Erba NC-2500 elemental analyzer coupled to a Finnigan DELTAplusXL isotope ratio mass spectrometer following approaches described previously in Lutz et al. (2016). Reproducibility was better than 0.1 wt % C, 0.01 wt % N, and <0.1‰ δ13C based on repeatedly measured standards.
Microscopy
Glutaraldehyde fixed sample aliquots were filtered through 0.2 μm polycarbonate filters that were subsequently stained with 4′,6′-diamidino-2-phenylindole (DAPI; 1 μg ml−1). Slides with filter sections were examined under a fluorescent microscope (Leica DM 2000) and images acquired with a DFC 420C camera and a FI/RH (Leica) filter system by following Winkel et al. (2018).
DNA Extraction, PCR Amplification, and Sequencing
DNA extractions were performed using the DNeasy PowerSoil Kit (QIAgen) according to the manufacturer’s protocol. DNA was eluted in 50 μl ultra-pure water. Extracted DNA concentration was measured on a Qubit 3.0 (Invitrogen) with the broad-range dsDNA kit (Invitrogen). Extracted DNA was amplified with bacterial primers for the 16S rRNA gene S-D-Bact-0341-b-S-17 (5′-CCTACGGGNGGCWGCAG-3′) and S-D-Bact-0785-a-A-21 (5′-GACTACHVGGGTATCTAATCC-3′; Herlemann et al., 2011); eukaryotic primers for the 18S rRNA gene 528F (5′-GCGGTAATTCCAGCTCCAA-3′) and 706R (5′-AATCCRAGAATTTCACCTCT-3′; Cheung et al., 2010), and the internal transcribed spacer 2 (ITS2) gene primers 5.8SbF (5′-GATGAAGAACGCAGCG-3′; Mikhailyuk et al., 2008) and ITS4R (5′-TCCTCCGCTTATTGATATGC-3′; White et al., 1990). Primers were tagged with Illumina adapter sequences. PCR was performed using KAPA HiFi HotStart ready mix (Roche) and the following PCR conditions: initial denaturation at 95°C for 3 min, 25 cycles (except samples collected in winter 2019 which were run for 35 cycles) of denaturation at 95°C for 30 s, annealing at 55°C for 30 s and elongation at 72°C for 30 s. Final elongation was run at 72°C for 5 min. All PCR runs were carried out in 25 μl reaction volumes. Amplified PCR products were sent to the Bristol Genomics Facility (United Kingdom), where they were barcoded with the Nextera XT Index kit. Pooled products were sequenced on the Illumina MiSeq platform using the V3 paired-end 2 × 300 bp chemistry.
Bioinformatics
Raw sequences were first quality checked with the FastQC program (Andrews, 2010) for quality and read length. All three targeted marker genes (16S rRNA, 18S rRNA, and ITS2) were analyzed with the DADA2 R package (Callahan et al., 2016). Documentation of the marker gene analysis pipelines is publicly available in the project’s GitHub repository.1 Taxonomy was assigned for 16S and 18S rRNA gene sequences using the SILVA database (Release 132; Quast et al., 2012), while ITS2 gene sequences were assigned using the UNITE database (version 8.0; UNITE Community, 2019). The most abundant ASVs were manually blasted against the NCBI database to cross-correlate results and obtain detailed taxonomic information. Detailed sequencing statistics including the number of sequences retained at each step as well as final sequences retrained for each sample can be found in Supplementary Table S1. Note that recently the glacier ice algae Mesotaenium has been renamed (now Ancylonema alaskana; Procházková et al., 2021) as have several bacterial phyla (e.g., Proteobacteria are now Pseudomonadota; Bacteriodetes are now Bacteriodota; Oren and Garrity, 2021). Nevertheless, in the current contribution, we used the original names, because at the time of data evaluations the notations in the databases used had not yet been updated.
In our low biomass winter season (February 2019) samples, we detected many known bacterial contaminants that we could show derived from our extraction kits and those were removed from the further analyses of amplicon sequence variants (ASVs; Callahan et al., 2017) as suggested by Salter et al. (2014); Glassing et al. (2016). We retained sequences that have been previously documented as dominant in snow and ice environments (Lutz et al., 2015a; Perini et al., 2019) including Polaromonas, Sediminibacterium, Pedobacter, Sphingomonas, Methylobacterium, Deinococcus, and Novosphingobium.
Statistical Analyses
Multivariate statistics were applied to environmental and sequencing data using with PAST v4 software (Hammer et al., 2001). Environmental parameters for principal component analysis (PCA) were log-transformed prior to statistical analyses, with the exception of pH according to Gobet et al. (2012). PCA was performed using an Euclidean distance matrix run with 1,000 iterations and PERMANOVA was used to test the difference between summer and winter samples. Sequencing data were Hellinger transformed (Legendre and Gallagher, 2001) prior to non-metric multidimensional scaling (NMDS) and hierarchical clustering. Diversity indices have been calculated on total ASVs with the implemented function of the PAST software.
Data Deposition
Marker gene sequences of the snow and ice samples are deposited in the NCBI Sequence Read Archive (SRA) under BioProject ID# PRJNA657180, under the accession numbers SRR12490806-SRR12490826 for 16S rRNA and 18S rRNA gene sequences, and SRR12589324-SRR12589336 for ITS2 gene sequences.
Results
Geochemical Characteristics
A total of 20 winter and summer snow and bare ice samples were collected in 2018 and 2019 from 5 different glaciers and ice sheets in Iceland (Figure 1; Table 1).
Overall, we found significant differences in geochemical signals and physical parameters (PERMANOVA, p = 0.0001) between summer and winter samples (Table 1; Figure 2). For example, conductivity values in winter (3.29–11.68 μs cm−1) were roughly double those in the summer samples (1.42–5.93 μs cm−1; Table 1), while pH values showed no seasonal trends (Figure 2). Comparison of the main dissolved ion concentrations revealed that snow had on average higher concentrations of sea-salt derived ions (Na, Cl, Mg, and K, Supplementary Figure S2 lower right quadrant), while ice contained higher concentrations of mineral-derived elements (e. g., Ca, Mn and Fe; Supplementary Figure S2, upper right quadrant). The major dissolved anions (Cl−, NO3−, and SO42−) varied dramatically seasonally (Table 2). Chloride was detected in all samples, although in the winter samples it was up to many orders of magnitude higher than in the summer samples. Nitrate and sulfate were below detection limits (50 ppb) in all samples collected during summer [with the exception of measurable nitrate in the Langjökull pro-glacial outflow water sample (IS18-6)], but were detectable in the winter samples. Notably, sulfate concentrations in the winter snow samples were up to 10 times higher than in the winter ice samples. DOC was variable in all surface snow and ice samples, but overall, the summer 2018 samples were characterized by up to 20 times higher concentrations of DOC compared to both the winter 2018 and the summer 2019 samples. Finally, trace elements Fe, Mo, Mn, and in part Zn, which are crucial for metabolic purposes, varied seasonally and by habitat, with higher concentrations in the winter snow samples compared to summer snow and winter and summer ice samples (Table 2). Most other elements showed no seasonal trends. Comparing these parameters by location and by habitat (snow vs. ice) showed that Skaftafellsjökull samples were separated from all other sampling locations (Supplementary Figure S1), but that the ice and snow from Skaftafellsjökull also plotted distantly from each other on a PCA plot (Supplementary Figures S1, S2). The other snow–ice sample (winter pairs from Sólheimajökull and summer pairs from Langjökull and Snæfellsjökull) all plotted closer to each other compared to Skaftafellsjökull. Also, of note is that sample IS19-12, which is a red snow sample from the snow–ice surface interface from Snæfellsjökull plotted apart from all other samples in all PCA plots (Supplementary Figures S1, S2; Figure 1 for the snow–ice interface image, bottom left).
Figure 2. Principal component analysis of data by season. Black vectors indicate environmental/geochemical parameters while differently colored dots represent season of sample collection: purple for summer 2018, light blue for winter 2019, and light red for summer 2019. A field blank of filtered distilled water is shown in grey.
Table 2. Geochemical measurements of organic and inorganic species in aqueous phase.
The difference between winter and summer samples was also mirrored in the total organic carbon (TOC %, Table 3) contents and isotope signature of the winter vs. the summer samples. The summer solid particulate samples contained far higher TC contents ranging from BDL to an exceptionally high value of 6.8% in the snow–ice interface sample IS19-12. In contrast, all winter samples were characterized by low TC values (max 0.41%, Table 3) that were on average one order of magnitude lower compared to the summer samples. In addition, all TN values in the winter samples were below detection limit of our analyses (<0.1%).
Table 3. Bulk measurements of carbon, nitrogen and carbon stable isotopes based on dry weight. C/N ratios calculated from the TC and TN values.
In terms of the carbon isotopic signatures, we found that the winter samples were characterized by heavier stable carbon isotope signals (summer average δ13C > ~ −27‰; and winter average δ13C < ~ −25‰). Differences were also noticeable between ice and snow pairs from the same glacier and sampled at the same time in the summer samples, with ice having heavier δ13C (Table 3). Winter samples showed similar values for both ice and snow. Because all winter sample had BDL TN contents, C/N vs. δ13C values could only be plotted for the summer samples (Supplementary Figure S5; Table 3).
Habitat and Seasonal Variation in Microbial Community Composition
Eukaryotic Communities
We targeted the 18S rRNA gene and ITS2 genes for amplification and sequencing of eukaryotes, but only had successful amplification of these genes in samples collected during summer. Overall, Fungi, Plantae, and Chromista dominated the eukaryotic community in both our snow and bare ice samples (Figure 3). The dominant fungi phyla Basidiomycota and Chytridiomycota made up between 19% and 50% of reads in all samples. The Chromista were dominated by Cercozoa (5%–17%), Ciliophora (1%–21%), and Ochrophyta (5%–35%). The latter group contained snow algae of the genus Hydrurus, commonly found in “yellow snow” (Remias et al., 2013). The Plantae consisted of the two phyla Chlorophyta (13%–65%) and Phragmoplastophyta (<1%–15%; Figure 3). The Chlorophyta include the typical snow algae of the genera Chloromonas, Chlainomonas, and Sanguina, which cause the often observed red and green snow blooms (Lutz et al., 2016; Remias et al., 2016; Engstrom et al., 2020). The dominant glacier ice algae in our samples were Phragmoplastophyta, which include the well-known glacier ice algae of the genera Mesotaenium and Ancylonema (Remias et al., 2009, 2012).
Figure 3. Stacked bar plots showing 18S rRNA gene relative abundance at the phylum level. Samples are ordered chronologically for the each glacier. Letters above bars indicate habitat with I = ice, S, snow; SI-I, snow–ice interface; and PW, pro-glacial water and numbers below bars indicate sample number (Table 1). Only summer samples are shown because 18S rRNA amplification of the winter 2019 samples did not yield a product.
We extracted algal ASVs and clustered them using Bray–Curtis dissimilarities, revealing that samples were mainly clustering according to differences in habitat—with the exception of the snow–ice interface sample IS19-12 (Figure 4A) and ice samples from Snæfellsjökull (IS18-1 and IS19-11). In sample IS19-12, more than 90% of the reads belonged to the genus Chlainomonas (Figure 4B), with its typical red xanthophyll rich pigmentation that causes the dark red color of this sample (Figure 4D). One cluster (IS18-6, IS19-14, IS18-3, IS18-8), which contained only ice samples and the pro-glacial water sample, had high numbers of the genera Mesotaenium and Ancylonema (Figure 4B). Summer ice samples from Snæfellsjökull (IS18-1 and IS19-11) contained negligible proportions of glacier ice algae and clustered with snow from the other summer samples (Figure 4A). Chloromonas and Hydrurus and the permafrost/soil algae Raphidonema were dominant in all snow samples and in the snow–ice interface sample IS18-4 (Figure 4B). Beta diversity analysis via NMDS of all samples showed a separation between snow and ice samples, with snow–ice interface samples (IS18-4 and IS19-12) plotting in between. The pro-glacial water sample clustered with the ice samples (Figure 4C).
Figure 4. Algae community structure based on 18S rRNA gene abundance. (A) Hierarchical clustering of 18S rRNA algae ASVs using Bray–Curtis and 1,000 iterations underlain by colored boxes that show season (red, purple), location (pink, orange, cyan), and habitat (light blue, green, gray, yellow). (B) Stacked barplots showing 18S rRNA gene ASVs relative abundance at the genus level with bars ordered according to hierarchical clustering. (C) NMDS plot of 18S rRNA gene ASVs using Bray–Curtis dissimilarities of Hellinger transformed ASVs with sampling sites color coded according to habitats as shown in legend in (A). (D) Close up image of field site IS19-12 showing a dominance of red snow with the 15 ml Falcon tube as size reference and the scale bar representing 12 cm.
In all summer samples, sequencing of the ITS2 gene showed a dominance of Chloromonas across all habitats and locations, with the exception of the snow–ice interface sample IS19-12, which was dominated by Chlainomonas (Supplementary Figure S3B). This result mirrors our 18S rRNA gene data (Figure 4), but the sequencing of the ITS2 gene missed the glacier ice algae of the phylum Phragmoplastophyta that were present in our 18S rRNA data (Figure 4 vs. Supplementary Figure S3A). Contradictory to the 18S rRNA gene sequences (Figure 4), we found that in the ITS our samples were dominated by Koliellaceae/Pseudochlorella (Supplementary Figure S3B), but other dominant taxa such as Hydrurus, seen in the 18S rRNA data, were missing (Figure 4B).
Bacterial Communities
Bacterial 16S rRNA genes were amplified and sequenced in all winter and summer samples and across all habitats. DNA concentrations were low or below Qubit detection levels (<0.01 μg ml−1) in all winter samples. Sequences were clustered using Bray–Curtis dissimilarities. We found that summer and winter bacterial communities were distinct from each other (Figure 5A). Bacterial communities in summer (2018 and 2019) were very similar to each other, only partially clustering by habitat type, while they clearly differed from winter communities (Figure 5).
Figure 5. Bacterial community structure based on 16S rRNA gene abundances. (A) Hierarchical clustering of 16S rRNA algae ASVs using Bray–Curtis and 1,000 iterations underlain by colored boxes that show season (red, purple, blue), locations (pink, orange, cyan, light purple, light yellow), and habitats (light blue, green, gray, yellow). (B) Stacked barplots showing 16S rRNA gene ASVs relative abundances at the class level. Bars have been ordered according to hierarchical clustering. (C) NMDS plot of 16S rRNA gene ASVs using Bray–Curtis dissimilarities of Hellinger transformed ASVs with sampling sites color coded per season as shown in legend in (A). (D) Microscopic images of DAPI-stained cells. Arrow points toward different cell morphologies: solid = rod, dotted = coccoid, and dashed = filamentous. Scale bar represents 10 μm. Direct imaging of DAPI-stained cells using epifluorescence microscopy revealed a range of microbial morphologies from cocci, and rod-shaped cells ((D), arrows in main image) to filaments ((D), inset).
Our summer samples were dominated by Bacteriodetes (22% and 49%), Gammaproteobacteria (2% and 45%), Alphaproteobacteria (12% and 30%), Cyanobacteria (4% and 26%), Actinobacteria (up to 8%), and Acidobacteria (up to 4%). The summer snow samples showed slightly lower proportions of Gammaproteobacteria, negligible Acidobacteria (Figure 5B), and had the lowest alpha diversity (Observed ASVs and Simpson index, Supplementary Figure S4). The summer ice samples clustered together with the snow–ice interface and the pro-glacial water samples and were characterized by low abundance taxa such as Planctomycetes, Gemmatimonadetes, Firmicutes, Deinococcus-Thermus, Chloroflexi, Verrucomicrobia, and Eremiobacteraeota (WPS-2). These were mostly absent in the summer snow samples cluster (Supplementary Figure S6). The winter samples generally had lower numbers of observed ASVs (~100), compared to the summer samples (between 50 and 500), and contained different taxa to the summer samples. Nevertheless, summer and winter snow and ice had similar Simpson indices (Supplementary Figure S4). In all winter samples, we saw an increased proportion of reads assigned to Actinobacteria (11%–38%), Firmicutes (up to 32%), Acidobacteria (1%–22%), Chloroflexi (1%–7%), Verrucomicrobia (up to 4%), and Armatimondetes (up to 2%); these phyla were almost exclusively found in our winter samples. Moreover, the proportion of reads for taxa that were dominant in the summer, such as Bacteriodetes, Gammaproteobacteria, and Alphaproteobacteria were on average 2.5-fold lower in the winter samples (Figure 5C). Some highly abundant taxa such as Cyanobacteria (more than 25% all ASVs in all samples with the exception of IS19-7 and IS19-3) showed similar abundances across both winter and summer. An NMDS plot using Bray–Curtis dissimilarities showed a separation between winter and summer samples: the summer samples clustered close to each other, while the winter samples were highly distinct (Figure 5C). Bacterial cells could be visualized by DAPI staining and showed a variety of different morphologies (Figure 5D).
Comparing winter and summer samples on an ASV level showed that they only share 30 ASVs, 7 of which belonging to chloroplast and mitochondria (Supplementary Table S2). The remaining 976 ASVs in our summer samples and 678 ASVs in our winter samples belonged to similar higher taxonomical levels (order to genus), but showed differences on ASV frequencies. For example, Bacteriodetes contained specific summer genera (Paludibacter, Prolixibacteraceae/BSV13, and WCHB1-32, Sedminibacterium, and Lentimicrobiaceae), whereas other genera (Alloprevotella, Terrimonas, and Flavisolibacter) only appeared in winter (Supplementary Figure S6). Similar patterns as for Bacteriodetes were found for Actinobacteria, Cyanobacteria, Alphaproteobacteria, and Gammaproteobacteria (Supplementary Table S2).
Discussion
Our data reveal clear seasonal shifts in the microbial community with a link to the geochemical characteristics of the glacier surface environment.
The three dominant eukaryotic phyla found in summer snow and ice samples (Chromista, Fungi, and Plantae) have been found previously on Icelandic glaciers (Lutz et al., 2015a) and elsewhere including snow from the Rocky Mountains (United States; Havig and Hamilton, 2019) and cryoconite holes in Queen Maud Land, Antarctica (Lutz et al., 2019), however in lower relative abundance than we detected here. The glacier ice algae Ancylonema, Mesotaenium, and Cylindrocystis (Figure 4) are also typical of glacier ice samples, as determined by conventional microscopic analysis (Takeuchi, 2001; Takeuchi and Kohshima, 2004; Uetake et al., 2010; Remias et al., 2012; Williamson et al., 2018). However, so far, only two studies have documented a large dominance of Ancylonema and Mesotaenium in glacier surface ice samples using molecular techniques (Lutz et al., 2018; McCutcheon et al., 2021). The clade Phragmoplastophyta, which was present in our 18S rRNA data but not in our ITS2 data (pointing toward a potential PCR primer bias) has previously been shown to comprise the majority of algae inhabiting glacier ice (Takeuchi, 2001; Takeuchi and Kohshima, 2004; Uetake et al., 2010; Remias et al., 2012; Lutz et al., 2018; Williamson et al., 2018; McCutcheon et al., 2021).
The summer snow algal community in our samples largely reflects findings from other similar (i.e., glacier surface snow) environments. In particular, Raphidonema and Hydrurus sp. have been shown to be common on Icelandic glaciers (Lutz et al., 2015a). Interestingly, the red-colored summer snow–ice interface sample from Snæfellsjökull (IS19-12) was dominated by Chlainomonas and not by the previously reported Chloromanas and Chlamydomonadaceae (Lutz et al., 2015a). Chlainomonas has rarely been documented among snow algal communities (Novis, 2002b), and has previously been reported in few locations (Hoham, 1975; Novis, 2002a, b; Remias et al., 2016; Engstrom et al., 2020). The abundance and diversity of algae found in our summer samples was not mirrored in our winter samples; they were fully absent from genomic data from our winter samples. Eukaryotic ice and snow algae require liquid water to bloom (Lutz et al., 2014; Anesio et al., 2017), and the scarcity of water during winter likely prohibits or slows algal metabolism. Moreover, in several of our winter sampling locations (i.e., Snæfellsjökull and Langjökull) glacier ice algae may have been present in deeper ice layers below the snow cover that we could not reach in our sampling efforts. The presence of such layers representing algal communities “frozen-in” at the end of the previous melt season is not well established and should be the target of future studies. Also, it remains unclear why glacial ice algae Ancylonema and Mesotaenium are so under-represented on Icelandic glaciers in comparison to the Greenland Ice Sheet or other Arctic settings where they are often found to dominate the glacier ice algae community (Lutz et al., 2017, 2018).
Bacteria were detected in all samples, but notably, the community composition differed significantly between summer and winter. The winter bacterial community was similar to air samples from Antarctica (Archer et al., 2019) and the Austrian Alps (Els et al., 2020), and also to fresh early spring snow from Svalbard, Norway (Larose et al., 2010; Maccario et al., 2014), based on previous studies that analyzed the same marker gene (16S rRNA) or the entire metagenome. The bacterial communities present in our summer samples shared similarities with those communities characterized in previous studies of summer-time glacial ecology (Hell et al., 2013; Lutz et al., 2015b; Hamilton and Havig, 2017; Havig and Hamilton, 2019).
Heterotrophic groups such as Bacteriodetes, Alphaproteobacteria, and Gammaproteobacteria dominated in all summer samples and belong to known groups that have been identified in marine systems to degrade algae-derived polysaccharides (Ruff et al., 2014; Teeling et al., 2016). These heterotrophs express carbohydrate-active enzymes known as CAZymes (Lombard et al., 2014), such as glycoside hydrolases (GH) that are most likely involved in the degradation of algae exudates (e.g., polysaccharides; Teeling et al., 2016). Cyanobacteria, while previously shown to inhabit cryoconite holes on glaciers and ice sheets (Stibal et al., 2008; Edwards et al., 2011), ice surfaces (Yallop et al., 2012), and some snow (Lutz et al., 2015a, 2016), often classify as Chloroplast in 16S rRNA gene sequencing. However, in this case the majority of sequences classified as Chloroplast are likely derived from the high abundance of glacier ice algae in our algal rich summer samples.
The variability in microbial communities between winter and summer was mirrored in the geochemical signals measured in the snow and ice samples. High concentrations of DOC were associated with summer algal blooms. During summer, glacial meltwater exports significant quantities of DOC (Lawson et al., 2014; Musilova et al., 2017), while during winter, freezing conditions on the glacier surface likely prevent supraglacial DOC export (Hood et al., 2009; Singer et al., 2012). It has recently been shown that primary production by supraglacial algae and other microorganisms over a summer melt season leads to an increase in the bioavailability of the meltwater-exported DOC (Kellerman et al., 2020). The high DOC concentrations measured in our summer 2018 snow samples (between 4 and 20 ppm, Table 2) are at the upper-most range of what has previously been measured for glacial ice (Hood et al., 2015), exceeding typical glacial snowpack and glacial outflow concentrations (Singer et al., 2012; Stubbins et al., 2012; Yan et al., 2016), and are largely in line with measured DOC concentrations in green snowpacks (i.e., snow containing abundant Chlorophyta) in Svalbard (Lutz et al., 2015b), and in ice samples from Antarctica (Barker et al., 2006). We suggest that the relatively high DOC concentrations in our samples compared to others (Singer et al., 2012; Stubbins et al., 2012; Hood et al., 2015; Yan et al., 2016) could arise from (i) our sampling at the end of the summer season (i.e., ~mid-September 2018) compared to other studies sampling between May and August (Singer et al., 2012; Stubbins et al., 2012; Wright et al., 2013; Lutz et al., 2015a,b; Yan et al., 2016), and (ii) the relatively high nutrient loading (and therefore high primary productivity) experienced by Icelandic glaciers due to volcanic input (Bagnato et al., 2013; Galeczka et al., 2017).
The combined organic δ13C signature and the C/N ratios for the summer samples provide an indication of the origin of the organic material (Lamb et al., 2006). Summer DOC (Table 2) is likely strongly influenced by microbially derived organic matter, evidenced by low C/N ratios (<13) and low δ13C-values (← ~ −25‰), typical of freshwater values (Lamb et al., 2006; Spijkerman et al., 2012). Compared to other glaciers in Svalbard and Sweden, the C/N and δ13C characteristics of our summer samples were similar to previously reported green snow and dirty ice samples (Lutz et al., 2017).
We found seasonal differences in the major macro- and micronutrients that are relevant for catabolic and anabolic cell processes (Goldman et al., 2012), for co-factors for enzyme activity, and that are used in anti-freeze protection mechanisms (Garnham et al., 2011). The concentrations of major dissolved ions in all our snow and ice samples (Table 2), fell within the range of concentrations reported in previous studies for Iceland snow and ice samples (Lutz et al., 2015a; Galeczka et al., 2017). In most cases, the ion concentrations were higher during winter than summer regardless of habitat (i.e., bare ice or snow; Table 2). Low summer nutrient and trace elements concentrations (e.g., Na, K, Mg, Ca, SO42−, NO3−, Mn, Mo, and Fe) correspond to the timing of algal blooms (i.e., summer). The availability of nitrogen and phosphorous are known to limit biological productivity of glacier surface ecosystems (Wolff, 2013; Williams et al., 2016; McCutcheon et al., 2021).
The high concentrations of NO32− detected in Snæfellsjökull, Skaftfellajökull, and Sólheimajökull winter snowpacks may arise from direct dry deposition of reactive inorganic nitrogen in snow (as observed in Svalbard; Kühnel et al., 2012, 2013; Björkman et al., 2013), biological nitrogen fixation (for example, by Cyano- and Proteobacteria; Telling et al., 2011), by cycling of organic nitrogen (Holland et al., 2019) or microbial necromass from glacial or pro-glacial environments (Bradley et al., 2016). While some chemical measurements point toward nitrification (Telling et al., 2014; Wadham et al., 2016), only one study directly identified putatively involved microorganisms (Segawa et al., 2014). Ammonia-oxidizing microorganisms such as Nitrosomonadaceae (Prosser et al., 2014) were present in our ice samples but were rarely found in our snow samples. Sulfate, in both our summer and winter samples, was likely derived from atmospheric deposition as sea-spray (e.g., Lutz et al., 2015a, 2016), or aerosol sulfate.
Conclusion and Outlook
Our data reveal seasonal characteristics of glacier surface ecosystems. We show a clear shift in the microbial community structure and corresponding geochemical environmental characteristics from the summer to the winter season.
Summer periods were characterized by a higher availability of microbially derived DOC and low nutrient concentrations, whereas winter was characterized by lower DOC and high nutrient concentrations.
The summer eukaryotic algae inhabiting the glacier surface ice and snow were distinct and were closely associated with a diverse community of heterotrophic bacteria likely relying on algal-delivered DOC. Winter microbial communities contained bacteria resembling communities typically found in air and snow and lacked eukaryotic algae.
The low abundance of glacial ice algae Ancylonema and Mesotaenium during summer in Iceland in comparison to the dominance in Greenlandic glaciers (Lutz et al., 2018; McCutcheon et al., 2021) is still unclear and should be further explored. Additionally, future work should identify the seasonal dynamics of glacier ice algae—specially to clarify whether their apparent absence in winter samples (as suggested in this study) is a consequence of sampling (i.e., algae may be present in deeper snow/ice interface layers) or is due to other ecological factors that drive algal distribution and abundance seasonally. In addition, summer and winter are punctuated by spring thaw and autumn freezing periods, which must also be studied to fully close the remaining gaps in understanding year-round characteristics of glacial ecosystems.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
MW and LGB designed the study. MW, CBT, and RM carried out laboratory work. MW and CBT performed bioinformatic analysis. MW wrote the first draft of the manuscript. All authors contributed to the sample collection in Iceland and contributed guidance and edits to the final manuscript preparation. All authors contributed to the article and approved the submitted version.
Funding
MW, RM, CBT, and LGB were funded through a Helmholtz Recruiting Initiative grant no. I-044-16-0. During the writing of this manuscript, LGB was also supported through an ERC Synergy grant (“DeepPurple” grant # 856416) from the European Research Council (ERC) funded under the European Union’s Horizon 2020 research and innovation program. JAB was part funded by an Alexander von Humboldt Foundation postdoctoral fellowship and a UK NERC grant (NE/T010967/1). The open access publication of this paper was supported within the funding program “Open Access Publikationskosten” Deutsche Forschungsgemeinschaft (DFG, German ResearchFoundation)—project number 491075472.
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.
Publisher’s Note
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.
Acknowledgments
We thank Kristin Günther, Silvia Pinkerneil, Stephen Reid, and Sven Peters from the GFZ for their laboratory support. We also thank Deirdre Clark and Sigurður Reynir Gíslason from the University of Iceland for logistical support during our field work in Iceland.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.876848/full#supplementary-material
Footnotes
References
Andrews, S. (2010). FastQC: A Quality Control Tool for High Throughput Sequence Data [Online]. Available at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (Accessed September 29, 2020).
Anesio, A. M., Lutz, S., Chrismas, N. A. M., and Benning, L. G. (2017). The microbiome of glaciers and ice sheets. NPJ Biofilms Microbiomes 3:10. doi: 10.1038/s41522-017-0019-0
Archer, S. D. J., Lee, K. C., Caruso, T., Maki, T., Lee, C. K., Cary, S. C., et al. (2019). Airborne microbial transport limitation to isolated antarctic soil habitats. Nat. Microbiol. 4, 925–932. doi: 10.1038/s41564-019-0370-4
Bagnato, E., Aiuppa, A., Bertagnini, A., Bonadonna, C., Cioni, R., Pistolesi, M., et al. (2013). Scavenging of sulphur, halogens and trace metals by volcanic ash: the 2010 Eyjafjallajökull eruption. Geochim. Cosmochim. Acta 103, 138–160. doi: 10.1016/j.gca.2012.10.048
Barker, J. D., Sharp, M. J., Fitzsimons, S. J., and Turner, R. J. (2006). Abundance and dynamics of dissolved organic carbon in glacier systems. Arct. Antarct. Alp. Res. 38, 163–172. doi: 10.1657/1523-0430(2006)38[163:AADODO]2.0.CO;2
Bellas, C. M., Anesio, A. M., and Barker, G. (2015). Analysis of virus genomes from glacial environments reveals novel virus groups with unusual host interactions. Front. Microbiol. 6:656. doi: 10.3389/fmicb.2015.00656
Benning, L. G., Anesio, A. M., Lutz, S., and Tranter, M. (2014). Biological impact on Greenland’s albedo. Nat. Geosci. 7:691. doi: 10.1038/ngeo2260
Björkman, M. P., Kühnel, R., Partridge, D. G., Roberts, T. J., Aas, W., Mazzola, M., et al. (2013). Nitrate dry deposition in Svalbard. Tellus B 65:19071. doi: 10.3402/tellusb.v65i0.19071
Box, J. E., Fettweis, X., Stroeve, J. C., Tedesco, M., Hall, D. K., and Steffen, K. (2012). Greenland ice sheet albedo feedback: thermodynamics and atmospheric drivers. Cryosphere 6, 821–839. doi: 10.5194/tc-6-821-2012
Bradley, J. A., Arndt, S., Šabacká, M., Benning, L. G., Barker, G. L., Blacker, J. J., et al. (2016). Microbial dynamics in a high arctic glacier forefield: a combined field, laboratory, and modelling approach. Biogeosciences 13, 5677–5696. doi: 10.5194/bg-13-5677-2016
Callahan, B. J., McMurdie, P. J., and Holmes, S. P. (2017). Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 11, 2639–2643. doi: 10.1038/ismej.2017.119
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Carpenter, E. J., Lin, S., and Capone, D. G. (2000). Bacterial activity in south pole snow. Appl. Environ. Microbiol. 66, 4514–4517. doi: 10.1128/AEM.66.10.4514-4517.2000
Cauvy-Fraunié, S., and Dangles, O. (2019). A global synthesis of biodiversity responses to glacier retreat. Nat. Ecol. Evol. 3, 1675–1685. doi: 10.1038/s41559-019-1042-8
Cheung, M. K., Au, C. H., Chu, K. H., Kwan, H. S., and Wong, C. K. (2010). Composition and genetic diversity of picoeukaryotes in subtropical coastal waters as revealed by 454 pyrosequencing. ISME J. 4, 1053–1059. doi: 10.1038/ismej.2010.26
Cook, J. M., Tedstone, A. J., Williamson, C., McCutcheon, J., Hodson, A. J., Dayal, A., et al. (2020). Glacier algae accelerate melt rates on the South-Western Greenland ice sheet. Cryosphere 14, 309–330. doi: 10.5194/tc-14-309-2020
Dumont, M., Brun, E., Picard, G., Michou, M., Libois, Q., Petit, J.-R., et al. (2014). Contribution of light-absorbing impurities in snow to Greenland’s darkening since 2009. Nat. Geosci. 7, 509–512. doi: 10.1038/ngeo2180
Edwards, A., Anesio, A. M., Rassner, S. M., Sattler, B., Hubbard, B., Perkins, W. T., et al. (2011). Possible interactions between bacterial diversity, microbial activity and supraglacial hydrology of cryoconite holes in Svalbard. ISME J. 5, 150–160. doi: 10.1038/ismej.2010.100
Els, N., Greilinger, M., Reisecker, M., Tignat-Perrier, R., Baumann-Stanzer, K., Kasper-Giebl, A., et al. (2020). Comparison of bacterial and fungal composition and their chemical interaction in free tropospheric air and snow Over an entire winter season at mount Sonnblick, Austria. Front. Microbiol. 11:980. doi: 10.3389/fmicb.2020.00980
Engstrom, C. B., Yakimovich, K. M., and Quarmby, L. M. (2020). Variation in snow algae blooms in the coast range of British Columbia. Front. Microbiol. 11:569. doi: 10.3389/fmicb.2020.00569
Galeczka, I., Eiriksdottir, E. S., Pálsson, F., Oelkers, E., Lutz, S., Benning, L. G., et al. (2017). Pollution from the 2014–15 Bárðarbunga eruption monitored by snow cores from the Vatnajökull glacier, Iceland. J. Volcanol. Geotherm. Res. 347, 371–396. doi: 10.1016/j.jvolgeores.2017.10.006
Garnham, C. P., Campbell, R. L., and Davies, P. L. (2011). Anchored clathrate waters bind antifreeze proteins to ice. Proc. Natl. Acad. Sci. U. S. A. 108, 7363–7367. doi: 10.1073/pnas.1100429108
Glassing, A., Dowd, S. E., Galandiuk, S., Davis, B., and Chiodini, R. J. (2016). Inherent bacterial DNA contamination of extraction and sequencing reagents may affect interpretation of microbiota in low bacterial biomass samples. Gut Pathog. 8:24. doi: 10.1186/s13099-016-0103-7
Gobet, A., Böer, S. I., Huse, S. M., van Beusekom, J. E. E., Quince, C., Sogin, M. L., et al. (2012). Diversity and dynamics of rare and of resident bacterial populations in coastal sands. ISME J. 6, 542–553. doi: 10.1038/ismej.2011.132
Goldman, A. D., Baross, J. A., and Samudrala, R. (2012). The enzymatic and metabolic capabilities of early life. PLoS One 7:e39912. doi: 10.1371/journal.pone.0039912
Hamilton, T. L., and Havig, J. (2017). Primary productivity of snow algae communities on stratovolcanoes of the Pacific northwest. Geobiology 15, 280–295. doi: 10.1111/gbi.12219
Hamilton, T. L., Peters, J. W., Skidmore, M. L., and Boyd, E. S. (2013). Molecular evidence for an active endogenous microbiome beneath glacial ice. ISME J. 7, 1402–1412. doi: 10.1038/ismej.2013.31
Hammer, Ø., Harper, D. A. T., and Ryan, P. D. (2001). Past: paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4.
Havig, J. R., and Hamilton, T. L. (2019). Snow algae drive productivity and weathering at volcanic rock-hosted glaciers. Geochim. Cosmochim. Acta 247, 220–242. doi: 10.1016/j.gca.2018.12.024
Hell, K., Edwards, A., Zarsky, J., Podmirseg, S. M., Girdwood, S., Pachebat, J. A., et al. (2013). The dynamic bacterial communities of a melting high Arctic glacier snowpack. ISME J. 7, 1814–1826. doi: 10.1038/ismej.2013.51
Herlemann, D. P., Labrenz, M., Jürgens, K., Bertilsson, S., Waniek, J. J., and Andersson, A. F. (2011). Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 5, 1571–1579. doi: 10.1038/ismej.2011.41
Hoham, R. W. (1975). Optimum temperatures and temperature ranges for growth of snow algae. Arct. Alp. Res. 7, 13–24. doi: 10.1080/00040851.1975.12003805
Holland, A. T., Williamson, C. J., Sgouridis, F., Tedstone, A. J., McCutcheon, J., Cook, J. M., et al. (2019). Dissolved organic nutrients dominate melting surface ice of the dark zone (Greenland ice sheet). Biogeosciences 16, 3283–3296. doi: 10.5194/bg-16-3283-2019
Hood, E., Battin, T. J., Fellman, J., O’Neel, S., and Spencer, R. G. M. (2015). Storage and release of organic carbon from glaciers and ice sheets. Nat. Geosci. 8, 91–96. doi: 10.1038/ngeo2331
Hood, E., Fellman, J., Spencer, R. G. M., Hernes, P. J., Edwards, R., D’Amore, D., et al. (2009). Glaciers as a source of ancient and labile organic matter to the marine environment. Nature 462, 1044–1047. doi: 10.1038/nature08580
Irvine-Fynn, T. D. L., Edwards, A., Stevens, I. T., Mitchell, A. C., Bunting, P., Box, J. E., et al. (2021). Storage and export of microbial biomass across the western Greenland ice sheet. Nat. Commun. 12:3960. doi: 10.1038/s41467-021-24040-9
Ji, F., Wu, Z., Huang, J., and Chassignet, E. P. (2014). Evolution of land surface air temperature trend. Nat. Clim. Chang. 4, 462–466. doi: 10.1038/nclimate2223
Kellerman, A. M., Hawkings, J. R., Wadham, J. L., Kohler, T. J., Stibal, M., Grater, E., et al. (2020). Glacier outflow dissolved organic matter as a window Into seasonally changing carbon sources: Leverett glacier, Greenland. J. Geophys. Res. Biogeosci. 125:e2019JG005161. doi: 10.1029/2019JG005161
Kühnel, R., Björkman, M. P., Vega, C. P., Hodson, A., Isaksson, E., and Ström, J. (2013). Reactive nitrogen and sulphate wet deposition at Zeppelin Station, Ny-Ålesund, Svalbard. Polar Res. 32:19136. doi: 10.3402/polar.v32i0.19136
Kühnel, R., Roberts, T. J., Björkman, M. P., Isaksson, E., Aas, W., Holmén, K., et al. (2012). 20-year climatology of and wet deposition at Ny-Ålesund, Svalbard. Adv. Meteorol. 2011:e406508. doi: 10.1155/2011/406508
Lamb, A. L., Wilson, G. P., and Leng, M. J. (2006). A review of coastal palaeoclimate and relative sea-level reconstructions using δ13C and C/N ratios in organic material. Earth-Sci. Rev. 75, 29–57. doi: 10.1016/j.earscirev.2005.10.003
Larose, C., Berger, S., Ferrari, C., Navarro, E., Dommergue, A., Schneider, D., et al. (2010). Microbial sequences retrieved from environmental samples from seasonal Arctic snow and meltwater from Svalbard, Norway. Extremophiles 14, 205–212. doi: 10.1007/s00792-009-0299-2
Larose, C., Dommergue, A., and Vogel, T. M. (2013). Microbial nitrogen cycling in Arctic snowpacks. Environ. Res. Lett. 8:35004. doi: 10.1088/1748-9326/8/3/035004
Lawson, E. C., Wadham, J. L., Tranter, M., Stibal, M., Lis, G. P., Butler, C. E. H., et al. (2014). Greenland ice sheet exports labile organic carbon to the Arctic oceans. Biogeosciences 11, 4015–4028. doi: 10.5194/bg-11-4015-2014
Lazzaro, A., Hilfiker, D., and Zeyer, J. (2015). Structures of microbial communities in alpine soils: seasonal and Elevational effects. Front. Microbiol. 6:1330. doi: 10.3389/fmicb.2015.01330
Legendre, P., and Gallagher, E. D. (2001). Ecologically meaningful transformations for ordination of species data. Oecologia 129, 271–280. doi: 10.1007/s004420100716
Lenaerts, J. T. M., Angelen, J. H., van Broeke, M. R., van den Gardner, A. S., Wouters, B., and van Meijgaard, E. (2013). Irreversible mass loss of Canadian Arctic archipelago glaciers. Geophys. Res. Lett. 40, 870–874. doi: 10.1002/grl.50214
Lombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M., and Henrissat, B. (2014). The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 42, D490–D495. doi: 10.1093/nar/gkt1178
Lutz, S., Anesio, A. M., Edwards, A., and Benning, L. G. (2015a). Microbial diversity on Icelandic glaciers and ice caps. Front. Microbiol. 6:307. doi: 10.3389/fmicb.2015.00307
Lutz, S., Anesio, A. M., Edwards, A., and Benning, L. G. (2017). Linking microbial diversity and functionality of arctic glacial surface habitats. Environ. Microbiol. 19, 551–565. doi: 10.1111/1462-2920.13494
Lutz, S., Anesio, A. M., Field, K., and Benning, L. G. (2015b). Integrated ‘Omics’, targeted metabolite and single-cell analyses of Arctic snow algae functionality and adaptability. Front. Microbiol. 6:1323. doi: 10.3389/fmicb.2015.01323
Lutz, S., Anesio, A. M., Raiswell, R., Edwards, A., Newton, R. J., Gill, F., et al. (2016). The biogeography of red snow microbiomes and their role in melting arctic glaciers. Nat. Commun. 7:11968. doi: 10.1038/ncomms11968
Lutz, S., Anesio, A. M., Villar, S. E. J., and Benning, L. G. (2014). Variations of algal communities cause darkening of a Greenland glacier. FEMS Microbiol. Ecol. 89, 402–414. doi: 10.1111/1574-6941.12351
Lutz, S., McCutcheon, J., McQuaid, J. B., and Benning, L. G. (2018). The diversity of ice algal communities on the Greenland ice sheet as revealed by oligotyping. Microb. Genom. 4. doi: 10.1099/mgen.0.000159
Lutz, S., Ziolkowski, L. A., and Benning, L. G. (2019). The biodiversity and geochemistry of Cryoconite holes in Queen Maud Land, East Antarctica. Microorganisms 7:160. doi: 10.3390/microorganisms7060160
Maccario, L., Vogel, T. M., and Larose, C. (2014). Potential drivers of microbial community structure and function in Arctic spring snow. Front. Microbiol. 5:413. doi: 10.3389/fmicb.2014.00413
McCutcheon, J., Lutz, S., Williamson, C., Cook, J. M., Tedstone, A. J., Vanderstraeten, A., et al. (2021). Mineral phosphorus drives glacier algal blooms on the Greenland ice sheet. Nat. Commun. 12:570. doi: 10.1038/s41467-020-20627-w
Meier, M. F., Dyurgerov, M. B., Rick, U. K., O’Neel, S., Pfeffer, W. T., Anderson, R. S., et al. (2007). Glaciers dominate Eustatic Sea-level rise in the 21st century. Science 317, 1064–1067. doi: 10.1126/science.1143906
Meredith, M., Sommerkorn, M., Cassotta, S., Derksen, C., Ekaykin, A., Hollowed, A., et al. (2019). “Polar Regions,” in IPCC Special Report on the Ocean and Cryosphere in a Changing Climate Available at: https://www.ipcc.ch/srocc/chapter/chapter-3-2/ (Accessed May 21, 2021).
Mikhailyuk, T. I., Sluiman, H. J., Massalski, A., Mudimu, O., Demchenko, E. M., Kondratyuk, S. Y., et al. (2008). New Streptophyte green algae from terrestrial habitats and an assessment of the genus nterfilum (Klebsormidiophyceae, Streptophyta)1. J. Phycol. 44, 1586–1603. doi: 10.1111/j.1529-8817.2008.00606.x
Musilova, M., Tranter, M., Wadham, J., Telling, J., Tedstone, A., and Anesio, A. M. (2017). Microbially driven export of labile organic carbon from the Greenland ice sheet. Nat. Geosci. 10, 360–365. doi: 10.1038/ngeo2920
Nghiem, S. V., Hall, D. K., Mote, T. L., Tedesco, M., Albert, M. R., Keegan, K., et al. (2012). The extreme melt across the Greenland ice sheet in 2012. Geophys. Res. Lett. 39. doi: 10.1029/2012GL053611
Novis, P. M. (2002a). Ecology of the snow alga Chlainomonas kolii (Chlamydomonadales, Chlorophyta) in New Zealand. Phycologia 41, 280–292. doi: 10.2216/i0031-8884-41-3-280.1
Novis, P. M. (2002b). New records of snow algae for New Zealand, from Mt philistine, Arthur’s pass National Park. N. Z. J. Bot. 40, 297–312. doi: 10.1080/0028825X.2002.9512789
Oren, A., and Garrity, G. M. Y. (2021). Valid publication of the names of forty-two phyla of prokaryotes. Int. J. Syst. Evol. Microbiol. 71:004645. doi: 10.1099/ijsem.0.005056
Perini, L., Gostinčar, C., and Gunde-Cimerman, N. (2019). Fungal and bacterial diversity of Svalbard subglacial ice. Sci. Rep. 9, 20230–20215. doi: 10.1038/s41598-019-56290-5
Procházková, L., Řezanka, T., Nedbalová, L., and Remias, D. (2021). Unicellular versus filamentous: The glacial alga Ancylonema alaskana comb. et stat. Nov. and its Ecophysiological relatedness to Ancylonema nordenskioeldii (Zygnematophyceae, Streptophyta). Microorganisms 9:1103. doi: 10.3390/microorganisms9051103
Prosser, J. I., Head, I. M., and Stein, L. Y. (2014). “The family Nitrosomonadaceae,” in The Prokaryotes: Alphaproteobacteria and Betaproteobacteria. eds. E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin, Heidelberg: Springer), 901–918.
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Regenspurg, S., Alawi, M., Blöcher, G., Börger, M., Kranz, S., Norden, B., et al. (2018). Impact of drilling mud on chemistry and microbiology of an upper Triassic groundwater after drilling and testing an exploration well for aquifer thermal energy storage in Berlin (Germany). Environ. Earth Sci. 77:516. doi: 10.1007/s12665-018-7696-8
Remias, D., Holzinger, A., Aigner, S., and Lütz, C. (2012). Ecophysiology and ultrastructure of Ancylonema nordenskiöldii (Zygnematales, Streptophyta), causing brown ice on glaciers in Svalbard (high arctic). Polar Biol. 35, 899–908. doi: 10.1007/s00300-011-1135-6
Remias, D., Holzinger, A., and Lütz, C. (2009). Physiology, ultrastructure and habitat of the ice alga Mesotaenium berggrenii (Zygnemaphyceae, Chlorophyta) from glaciers in the European Alps. Phycologia 48, 302–312. doi: 10.2216/08-13.1
Remias, D., Jost, S., Boenigk, J., Wastian, J., and Lütz, C. (2013). Hydrurus-related golden algae (Chrysophyceae) cause yellow snow in polar summer snowfields. Phycol. Res. 61, 277–285. doi: 10.1111/pre.12025
Remias, D., Pichrtová, M., Pangratz, M., Lütz, C., and Holzinger, A. (2016). Ecophysiology, secondary pigments and ultrastructure of Chlainomonas sp. (Chlorophyta) from the European Alps compared with Chlamydomonas nivalis forming red snow. FEMS Microbiol. Ecol. 92:fiw030. doi: 10.1093/femsec/fiw030
Ruff, E. S., Probandt, D., Zinkann, A.-C., Iversen, M. H., Klaas, C., Würzberg, L., et al. (2014). Indications for algae-degrading benthic microbial communities in deep-sea sediments along the Antarctic polar front. Deep-Sea Res. II: Top. Stud. Oceanogr. 108, 6–16. doi: 10.1016/j.dsr2.2014.05.011
Salter, S. J., Cox, M. J., Turek, E. M., Calus, S. T., Cookson, W. O., Moffatt, M. F., et al. (2014). Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 12:87. doi: 10.1186/s12915-014-0087-z
Schostag, M., Stibal, M., Jacobsen, C. S., Bælum, J., Taş, N., Elberling, B., et al. (2015). Distinct summer and winter bacterial communities in the active layer of Svalbard permafrost revealed by DNA- and RNA-based analyses. Front. Microbiol. 6:399. doi: 10.3389/fmicb.2015.00399
Segawa, T., Ishii, S., Ohte, N., Akiyoshi, A., Yamada, A., Maruyama, F., et al. (2014). The nitrogen cycle in cryoconites: naturally occurring nitrification-denitrification granules on a glacier. Environ. Microbiol. 16, 3250–3262. doi: 10.1111/1462-2920.12543
Singer, G. A., Fasching, C., Wilhelm, L., Niggemann, J., Steier, P., Dittmar, T., et al. (2012). Biogeochemically diverse organic matter in alpine glaciers and its downstream fate. Nat. Geosci. 5, 710–714. doi: 10.1038/ngeo1581
Smith, H. J., Foster, R. A., McKnight, D. M., Lisle, J. T., Littmann, S., Kuypers, M. M. M., et al. (2017). Microbial formation of labile organic carbon in Antarctic glacial environments. Nat. Geosci. 10, 356–359. doi: 10.1038/ngeo2925
Spijkerman, E., Wacker, A., Weithoff, G., and Leya, T. (2012). Elemental and fatty acid composition of snow algae in Arctic habitats. Front. Microbiol. 3:380. doi: 10.3389/fmicb.2012.00380
Stibal, M., Bradley, J. A., Edwards, A., Hotaling, S., Zawierucha, K., Rosvold, J., et al. (2020). Glacial ecosystems are essential to understanding biodiversity responses to glacier retreat. Nat. Ecol. Evol. 4, 686–687. doi: 10.1038/s41559-020-1163-0
Stibal, M., Šabacká, M., and Žárský, J. (2012). Biological processes on glacier and ice sheet surfaces. Nat. Geosci. 5, 771–774. doi: 10.1038/ngeo1611
Stibal, M., Tranter, M., Benning, L. G., and Rehak, J. (2008). Microbial primary production on an Arctic glacier is insignificant in comparison with allochthonous organic carbon input. Environ. Microbiol. 10, 2172–2178. doi: 10.1111/j.1462-2920.2008.01620.x
Stubbins, A., Hood, E., Raymond, P. A., Aiken, G. R., Sleighter, R. L., Hernes, P. J., et al. (2012). Anthropogenic aerosols as a source of ancient dissolved organic matter in glaciers. Nat. Geosci. 5, 198–201. doi: 10.1038/ngeo1403
Takeuchi, N. (2001). The altitudinal distribution of snow algae on an Alaska glacier (Gulkana glacier in the Alaska range). Hydrol. Process. 15, 3447–3459. doi: 10.1002/hyp.1040
Takeuchi, N., and Kohshima, S. (2004). A snow algal community on Tyndall glacier in the southern Patagonia Icefield, Chile. Arct. Antarct. Alp. Res. 36, 92–99. doi: 10.1657/1523-0430(2004)036[0092:ASACOT]2.0.CO;2
Tedstone, A. J., Cook, J. M., Williamson, C. J., Hofer, S., McCutcheon, J., Irvine-Fynn, T., et al. (2020). Algal growth and weathering crust state drive variability in western Greenland ice sheet ice albedo. Cryosphere 14, 521–538. doi: 10.5194/tc-14-521-2020
Teeling, H., Fuchs, B. M., Bennke, C. M., Krüger, K., Chafee, M., Kappelmann, L., et al. (2016). Recurring patterns in bacterioplankton dynamics during coastal spring algae blooms. elife 5:e11888. doi: 10.7554/ELIFE.11888
Telling, J., Anesio, A. M., Tranter, M., Fountain, A. G., Nylen, T., Hawkings, J., et al. (2014). Spring thaw ionic pulses boost nutrient availability and microbial growth in entombed Antarctic Dry Valley cryoconite holes. Front. Microbiol. 5:694. doi: 10.3389/fmicb.2014.00694
Telling, J., Anesio, A. M., Tranter, M., Irvine-Fynn, T., Hodson, A., Butler, C., et al. (2011). Nitrogen fixation on Arctic glaciers, Svalbard. J. Geophys. Res. Biogeosci. 116. doi: 10.1029/2010JG001632
Uetake, J., Naganuma, T., Hebsgaard, M. B., Kanda, H., and Kohshima, S. (2010). Communities of algae and cyanobacteria on glaciers in West Greenland. Policy. Sci. 4, 71–80. doi: 10.1016/j.polar.2010.03.002
UNITE Community (2019). UNITE QIIME release for eukaryotes. Version 18.11.2018. UNITE Community. Available at: https://doi.org/10.15156/BIO/786335
Vigneron, A., Lovejoy, C., Cruaud, P., Kalenitchenko, D., Culley, A., and Vincent, W. F. (2019). Contrasting winter versus summer microbial communities and metabolic functions in a permafrost thaw Lake. Front. Microbiol. 10:1656. doi: 10.3389/fmicb.2019.01656
Wadham, J. L., Hawkings, J., Telling, J., Chandler, D., Alcock, J., O’Donnell, E., et al. (2016). Sources, cycling and export of nitrogen on the Greenland ice sheet. Biogeosciences 13, 6339–6352. doi: 10.5194/bg-13-6339-2016
White, T. J., Bruns, T., Lee, S., and Taylor, J. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics,” in PCR Protocols. eds. Innis, M. A., Gelfand, D. H., Sninsky, J. J., and White, T. J. (Academic Press Inc.), 315–322.
Williams, J. J., Gourmelen, N., and Nienow, P. (2020). Dynamic response of the Greenland ice sheet to recent cooling. Sci. Rep. 10:1647. doi: 10.1038/s41598-020-58355-2
Williams, M. W., Wilson, A., Tshering, D., Thapa, P., and Kayastha, R. B. (2016). Using geochemical and isotopic chemistry to evaluate glacier melt contributions to the Chamkar Chhu (river), Bhutan. Ann. Glaciol. 57, 339–348. doi: 10.3189/2016AoG71A068
Williamson, C. J., Anesio, A. M., Cook, J., Tedstone, A., Poniecka, E., Holland, A., et al. (2018). Ice algal bloom development on the surface of the Greenland ice sheet. FEMS Microbiol. Ecol. 94:fiy025. doi: 10.1093/femsec/fiy025
Williamson, C. J., Cook, J., Tedstone, A., Yallop, M., McCutcheon, J., Poniecka, E., et al. (2020). Algal photophysiology drives darkening and melt of the Greenland ice sheet. Proc. Natl. Acad. Sci. U. S. A. 117, 5694–5705. doi: 10.1073/pnas.1918412117
Winkel, M., Mitzscherling, J., Overduin, P. P., Horn, F., Winterfeld, M., Rijkers, R., et al. (2018). Anaerobic methanotrophic communities thrive in deep submarine permafrost. Sci. Rep. 8:1291. doi: 10.1038/s41598-018-19505-9
Wolff, E. W. (2013). Ice sheets and nitrogen. Philos. Trans. R. Soc. B: Biol. Sci. 368:20130127. doi: 10.1098/rstb.2013.0127
Wright, K. E., Williamson, C., Grasby, S. E., Spear, J. R., and Templeton, A. S. (2013). Metagenomic evidence for sulfur lithotrophy by Epsilonproteobacteria as the major energy source for primary productivity in a sub-aerial arctic glacial deposit, Borup Fiord Pass. Front. Microbiol. 4:63. doi: 10.3389/fmicb.2013.00063
Yallop, M. L., Anesio, A. M., Perkins, R. G., Cook, J., Telling, J., Fagan, D., et al. (2012). Photophysiology and albedo-changing potential of the ice algal community on the surface of the Greenland ice sheet. ISME J. 6, 2302–2313. doi: 10.1038/ismej.2012.107
Keywords: glacier ice algae, snow algae, seasonality, microbial dynamics, cryosphere
Citation: Winkel M, Trivedi CB, Mourot R, Bradley JA, Vieth-Hillebrand A and Benning LG (2022) Seasonality of Glacial Snow and Ice Microbial Communities. Front. Microbiol. 13:876848. doi: 10.3389/fmicb.2022.876848
Edited by:
Pietro Buzzini, University of Perugia, ItalyReviewed by:
Maria Papale, Institute of Polar Sciences (CNR), ItalyIan Hawes, University of Waikato, New Zealand
Copyright © 2022 Winkel, Trivedi, Mourot, Bradley, Vieth-Hillebrand and Benning. 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: Matthias Winkel, bXdpbmtlbDE5ODJAZ21haWwuY29t; Liane G. Benning, YmVubmluZ0BnZnotcG90c2RhbS5kZQ==
†Present address: Matthias Winkel, Analytic and Toxikogenomics, BfR, Federal Institute for Risk Assessment, Berlin, Germany
‡These authors have contributed equally to this work