The Snowmelt Niche Differentiates Three Microbial Life Strategies That Influence Soil Nitrogen Availability During and After Winter

Sorensen, Patrick O.; Beller, Harry R.; Bill, Markus; Bouskill, Nicholas J.; Hubbard, Susan S.; Karaoz, Ulas; Polussa, Alexander; Steltzer, Heidi; Wang, Shi; Williams, Kenneth H.; Wu, Yuxin; Brodie, Eoin L.

doi:10.3389/fmicb.2020.00871

ORIGINAL RESEARCH article

Front. Microbiol., 15 May 2020

Sec. Terrestrial Microbiology

Volume 11 - 2020 | https://doi.org/10.3389/fmicb.2020.00871

The Snowmelt Niche Differentiates Three Microbial Life Strategies That Influence Soil Nitrogen Availability During and After Winter

$\r\nPatrick O. Sorensen*$ Patrick O. Sorensen^1*

Harry R. Beller¹

Markus Bill¹

Nicholas J. Bouskill¹

Susan S. Hubbard¹

Ulas Karaoz¹

Alexander Polussa^1,2

Heidi Steltzer^3,4

Shi Wang¹

Kenneth H. Williams^1,4

Yuxin Wu¹

Eoin L. Brodie^1,5*

¹Earth and Environmental Sciences, Lawrence Berkeley National Laboratory, Berkeley, CA, United States
²School of Forestry and Environmental Studies, Yale University, New Haven, CT, United States
³Fort Lewis College, Durango, CO, United States
⁴Rocky Mountain Biological Laboratory, Gothic, CO, United States
⁵Department of Environmental Science, Policy and Management, University of California, Berkeley, Berkeley, CA, United States

Soil microbial biomass can reach its annual maximum pool size beneath the winter snowpack and is known to decline abruptly following snowmelt in seasonally snow-covered ecosystems. Observed differences in winter versus summer microbial taxonomic composition also suggests that phylogenetically conserved traits may permit winter- versus summer-adapted microorganisms to occupy distinct niches. In this study, we sought to identify archaea, bacteria, and fungi that are associated with the soil microbial bloom overwinter and the subsequent biomass collapse following snowmelt at a high-altitude watershed in central Colorado, United States. Archaea, bacteria, and fungi were categorized into three life strategies (Winter-Adapted, Snowmelt-Specialist, Spring-Adapted) based upon changes in abundance during winter, the snowmelt period, and after snowmelt in spring. We calculated indices of phylogenetic relatedness (archaea and bacteria) or assigned functional attributes (fungi) to organisms within life strategies to infer whether phylogenetically conserved traits differentiate Winter-Adapted, Snowmelt-Specialist, and Spring-Adapted groups. We observed that the soil microbial bloom was correlated in time with a pulse of snowmelt infiltration, which commenced 65 days prior to soils becoming snow-free. A pulse of nitrogen (N, as nitrate) occurred after snowmelt, along with a collapse in the microbial biomass pool size, and an increased abundance of nitrifying archaea and bacteria (e.g., Thaumarchaeota, Nitrospirae). Winter- and Spring-Adapted archaea and bacteria were phylogenetically clustered, suggesting that phylogenetically conserved traits allow Winter- and Spring-Adapted archaea and bacteria to occupy distinct niches. In contrast, Snowmelt-Specialist archaea and bacteria were phylogenetically overdispersed, suggesting that the key mechanism(s) of the microbial biomass crash are likely to be density-dependent (e.g., trophic interactions, competitive exclusion) and affect organisms across a broad phylogenetic spectrum. Saprotrophic fungi were the dominant functional group across fungal life strategies, however, ectomycorrhizal fungi experienced a large increase in abundance in spring. If well-coupled plant-mycorrhizal phenology currently buffers ecosystem N losses in spring, then changes in snowmelt timing may alter ecosystem N retention potential. Overall, we observed that snowmelt separates three distinct soil niches that are occupied by ecologically distinct groups of microorganisms. This ecological differentiation is of biogeochemical importance, particularly with respect to the mobilization of nitrogen during winter, before and after snowmelt.

Introduction

Snowmelt is a hydrologic event that exerts significant influence on annual water and nitrogen (N) export in seasonally snow-covered, mountainous catchments (Bales et al., 2006). Historically, the timing of snowmelt has coincided with a suite of environmental conditions that characterize the transition from winter to spring, including rising air temperature, longer days having higher photosynthetically active radiation, and greater soil moisture availability (Marks and Dozier, 1992; Harpold and Molotch, 2015; Musselman et al., 2017). These abiotic cues initiate phenological transitions between winter and spring metabolic states for both plants and soil microbial communities (Richardson et al., 2006; Inouye, 2008; Miller-Rushing et al., 2010; Contosta et al., 2016; Ladwig et al., 2016; Thackeray et al., 2016). As a consequence, snowmelt is found to be a critical period that influences nutrient mobilization, assimilation, and retention in terrestrial ecosystems (Brooks et al., 1998; Brooks and Williams, 1999; Grogan and Jonasson, 2003; Kielland et al., 2006; Campbell et al., 2007). Rising global air temperature has led to reductions in winter snowpack extent, earlier snowmelt dates, and altered growing season length in many mountainous catchments (Mote et al., 2005; Steltzer and Post, 2009; Harte et al., 2015; Sloat et al., 2015; Bormann et al., 2018; Hock et al., 2019; Prevéy et al., 2019). The ecological consequences of such environmental changes, however, are not well understood (Ernakovich et al., 2014).

Soil microbial biomass production beneath the winter snowpack can be significant, and in some ecosystems, microbial biomass may reach its annual maximum pool size beneath the winter snowpack (Brooks et al., 1998; Lipson et al., 1999; Grogan and Jonasson, 2003; Edwards et al., 2006; Larsen et al., 2007; Buckeridge and Grogan, 2010; Buckeridge et al., 2013). The microbial bloom in soil during winter immobilizes N within microbial biomass (Brooks et al., 1998). An abrupt collapse in the microbial biomass pool size after snowmelt results in the subsequent release of N from the microbial biomass pool (Lipson et al., 1999; Buckeridge et al., 2010; Isobe et al., 2018). Numerous mechanisms have been hypothesized to induce the microbial biomass crash, including substrate depletion and starvation of the winter community, cell lysis due to soil freeze-thaw events, sudden changes in the osmotic environment with snowmelt, grazing by protozoa and mesofauna, and mortality and replacement of winter-adapted psychrophiles by summer-adapted mesophiles (Jefferies et al., 2010). A better understanding of the factors that promote microbial biomass production under-snow and the crash after snowmelt is necessary because the flux of N released from microbial biomass after snowmelt can be the largest annual pulse of soil N in some ecosystems (Grogan and Jonasson, 2003; Schmidt et al., 2007).

In addition to differences in the size of the microbial biomass pool across seasons, bacterial and fungal community composition are also known to differ taxonomically in winter compared to summer (Schadt et al., 2003; Lipson and Schmidt, 2004; Aanderud et al., 2013; Isobe et al., 2018), although this has typically been reported at coarse (e.g., phylum) levels of taxonomic resolution. For example, Acidobacteria, Verrucomicrobia, and Bacteriodetes co-exist in alpine soils, but strong successional differences in relative abundance are observed among these phyla during winter and summer (Schmidt et al., 2007). Because two species competing for the same resources cannot persist together indefinitely, such co-occurrence and succession implies that these phyla do not occupy the same niche or partition resources in space or time (Hutchinson, 1957). Differences in growth rates and biomass yields among winter and summer bacterial isolates also indicate that phenotypic traits with physiological trade-offs may differentiate winter- and summer-adapted microorganisms (Lipson et al., 2009).

More broadly, an organism’s fitness across varying environmental gradients should be reflected by its abundance in the environment and constrained by a suite of underlying physiological traits, collectively referred to as an organism’s life strategy (Schimel et al., 2007; Placella et al., 2012; Evans et al., 2014; Ho et al., 2017). If such traits are adaptive and phylogenetically conserved (Webb et al., 2002), then quantifying the phylogenetic relatedness of organisms grouped into life strategies may provide insights into the community assembly processes (e.g., niche partitioning) that underlie the microbial bloom beneath the winter snowpack and subsequent biomass collapse after snowmelt.

Within this framework we addressed the following questions: (1) are closely related microorganisms responsible for the overwinter microbial biomass bloom, or is the phenomenon widespread and characteristic among distantly related taxa? (2) what is the mechanism responsible for the microbial biomass crash following snowmelt and how widespread is the snowmelt decline among distantly related bacterial, archaeal, and fungal taxa? (3) are there identifiable temporal abundance patterns that differentiate microorganisms into life strategies that can be used to infer an organism’s snowmelt niche?

To address these questions, we sampled soils along a 200-m upland hillslope-to-riparian floodplain transect in a mountainous catchment at the East River Watershed, Colorado, United States (Hubbard et al., 2018). Soil samples were collected over a time-course that began after plant senescence in autumn, through snow accumulation during winter, during snowmelt, and through plant green-up in spring and early summer. We monitored soil temperature and moisture and quantified the dynamics of the soil microbial biomass pool as well as extractable N in soil that was measured as extractable soil nitrate ( ${NO}_{3}^{-}$ ). Bacterial, archaeal, and fungal species relative abundances were also measured over this seasonal time-course. Here, we show that soil temperature and soil moisture during winter are controlled by winter snowpack dynamics and that snowmelt induces rapid changes in the soil abiotic environment. In addition, snowmelt triggered marked increases in the size of the microbial biomass pool, which also reduced soil N availability. Finally, we show that specific bacterial, archaeal, and fungal taxa are associated with the rise and fall of microbial biomass overwinter and that the abundance responses can be used to differentiate taxa into three distinct life strategies separated in time by the snowmelt event.

Materials and Methods

The Upper East River Watershed is located in Gunnison County, CO, United States in the West Elk Mountains near the towns of Crested Butte and Gothic (38°57.5′ N, 106°59.3′ W) and is home to the Rocky Mountain Biological Laboratory. Elevation at the East River watershed spans from 2750 to 4000 m. The climate at East River is continental and characterized by periods of persistent snow cover in winter lasting 4 to 6 months (e.g., November through May) followed by an arid summer with intermittent, monsoonal rain events that occur from July through September. The lowest daily air temperatures typically occur in January (−14 ± 3°C), whereas the highest daily air temperatures typically occur in July (23.5 ± 3°C). The average summer air temperature has increased by 0.5 ± 0.1°C per decade since 1974 (CaraDonna and Inouye, 2015). The average annual precipitation is approximately 1200 cm, with the majority (>80%) falling as snowfall during winter (Harte and Shaw, 1995; Hubbard et al., 2018). The maximum annual snow depth between 1974 and 2017 was 465 cm and the date of snowmelt in spring has advanced by 3.5 ± 2 days per decade since 1974 (Iler et al., 2013).

We sampled a 200-m transect that transitioned from an upland hillslope (hereafter referred to as Hillslope) to a riparian floodplain (hereafter, Floodplain) adjacent to the main stem of the East River (elevation ∼2760 to 2775 m). We chose this transect because these ecosystem types are representative of dominant land cover adjacent to the East River. Vegetation in the Hillslope is a mix of perennial bunchgrasses (e.g., Festuca arizonica), forbs (e.g., Lupinus spp., Potentilla gracilis, Veratrum californicum), and shrubs (Artemisia tridentata), whereas plant communities in the Floodplain are dominated by dwarf birch (Betula grandulosa) and mountain willow (Salix monticola, Falco et al., 2019).

Soil temperature, soil volumetric water content, and water potential were measured continuously starting in October 2016 at three locations on the Hillslope and at three locations on the Floodplain. Soil temperature was measured continuously at 6 and 9 cm below the soil surface (sensor model RT-1, 5TE and MPS6, METER). Soil volumetric water content was measured at 9 cm using a time-domain reflectometry probe (model 5TE, METER Group). Soil water potential was measured at 17 cm using MPS6 from METER. We measured snow depth during winter with a marked snow pole or meter stick on the dates of soil sampling that had snow cover at the site.

Soil Sampling and Field Processing

Soil samples were collected on four dates to characterize the soil microbial biomass carbon (C) pool, extractable soil ${NO}_{3}^{-}$ , and microbial community composition, starting first in autumn after plant senescence (September 12, 2016), at peak winter snow depth (March 7, 2017), during the snowmelt period (May 9, 2017), and following the complete loss of snow and the start of the plant growing season (June 9, 2017). During snow-free times of the year, soils were collected using a 4 cm diameter soil bulk density corer at 12 plots at the Hillslope and at 3 plots at the Floodplain. The six plots on the Hillslope were encompassed by a representative 80 m × 100 m area that was approximately 100 m upslope from the edge of the Floodplain. Each plot on the Hillslope was 20 m × 20 m and separated by at least 2.5 m. The Hillslope plots were arranged in three pairs that ran parallel to the river. The paired plots will serve as experimental blocks for a future field experiment. The three plots in the Floodplain were similarly 20 m × 20 m and arranged in a 100 m × 30 m transect that was perpendicular to the river and spanned the area from the toeslope-floodplain transition to near the edge of the river. Soil cores were subsampled and split into three discrete depth increments; 0 to 5 cm, 5 to 15 cm, and 15 cm + below the soil surface. These relatively shallow soils are best described as A-horizons at the Hillslope and Floodplain locations. Soil horizon development and differentiation is greatly reduced because of substantial bioturbation from burrowing mammals on the Hillslope and periodic sediment deposition during spring flooding events in the Floodplain.

During periods of winter snow cover (i.e., March and May 2017), snow-pits were dug down to the soil surface at three locations on the Hillslope and one location in the Floodplain in order to sample soils from beneath the snowpack. In each snow-pit, soils were sampled at two adjacent locations separated by more than 1 m using the soil coring method described above. New snow pits were dug in different locations on the Hillslope and Floodplain during snowmelt in May 2017 to avoid repeatedly sampling the same locations that were sampled at peak winter snow depth in March 2017. Thus, during the snow-free time of the year, we sampled and analyzed 180 total soil samples (2 ecosystem types × 15 Plots × 2 Time Points × 3 Depths) and 48 soils total during winter (4 Snow-pits × 2 Time Points × 3 Depths × 2 replicate cores). A ∼10 g subsample from each soil core at each depth was placed immediately on dry ice in the field, frozen, and archived for archaeal, bacterial and fungal community analysis. The remainder of the soil core was allocated to physical and chemical characterization (described below) and was stored at 4°C until further analysis.

Soil Physical and Chemical Properties

Gravimetric soil moisture content was measured for each sample on each sampling date by determining the mass of water lost during a 48-h incubation at 60°C. Soil pH was determined in a 1:1 w/v slurry in reagent water (18 MΩ resistance; GenPure UltraPure Water System; Thermo Scientific) using an Orion soil pH probe (Thermo Scientific). Total soil organic C and N concentrations and soil δ¹³C and δ¹⁵N were measured using a Costech ECS 4010 elemental analyzer coupled to a Thermo Delta V Plus isotope ratio mass spectrometer at the Center for Isotope Geochemistry at Lawrence Berkeley National Laboratory (Berkeley, CA, United States).

Extractable nitrate ( ${NO}_{3}^{-}$ ) as well as microbial biomass carbon (C) were quantified in 0.5M K₂SO₄ extracts. A 5 g subsample of each soil (field-moist) was extracted in 25 mL of 0.5M K₂SO₄ on an orbital shaker table for 60 min, then gravity filtered through pre-leached #42 Whatman filter paper and frozen until further analysis. ${NO}_{3}^{-}$ concentration in the 0.5M K₂SO₄ extracts was measured colorimetrically using a Versamax microplate spectrophotometer (Molecular Devices) using a modified Greiss reaction protocol (Doane and Horwath, 2003). Microbial biomass carbon (C) was measured using the chloroform-fumigation extraction method (Brookes et al., 1985). A 5-g subsample was fumigated with ethanol-free chloroform for 7 days and then extracted as stated above. Microbial biomass C was estimated as organic C measured in the fumigated extract minus organic C measured in the non-fumigated sample (Brookes et al., 1985). We did not apply an extraction correction to account for incomplete microbial biomass lysis during the fumigation. The concentration of dissolved organic C in the fumigated and non-fumigated samples was quantified using the Mn(III) pyrophosphate oxidation method (Bartlett and Ross, 1988). We applied a universal correction factor to the measured organic C concentrations to account for differences in reaction efficiency across soil types (Giasson et al., 2014). All concentrations were corrected for the field-moist water content of the soil and values are reported based on the oven-dried soil mass used for the extractions.

Illumina Library Preparation and Bioinformatics

Total genomic DNA was extracted in triplicate from each soil sample collected on each date using the DNeasy PowerSoil DNA Extraction Kit (QIAGEN) with the protocol modified to include a 5-min incubation at 65°C prior to bead-beating to increase lysis efficiency. Replicate soil DNA extracts were combined prior to PCR amplification. PCR reaction conditions are summarized in the Supplementary Material (see Supplementary Datasheet S1). Purified PCR products were pooled in equimolar concentrations and sequenced on a single lane for 300-bp paired-end v3 Illumina MiSeq sequencing conducted at the Vincent J. Coates Genomics Sequencing Laboratory at the University of California, Berkeley.

Forward and reverse reads were aligned and paired using usearch [v8.1.1861 (Edgar, 2010)] fastq_mergepairs command (maximum diff = 3). The aligned reads were quality filtered (command fastq_filter with -fastq_trunclen = 230, -fastq_maxee = 0.1), concatenated into a single fasta file, and singletons were removed (command sortbysize with minsize = 2). These filtered sequences were used for operational taxonomic unit (OTU) clustering with the uparse pipeline (Edgar, 2013) setting the OTU cut-off threshold to 97%. Chimeric sequences were filtered with uchime (Edgar et al., 2011) with reference to the ChimeraSlayer database downloaded from http://drive5.com/uchime/gold.fa. OTU abundances across individual samples were calculated by mapping chimera-filtered OTUs against the quality-filtered reads (command usearch_global with -strand plus -id 0.97). Taxonomy was assigned to each OTU by a Naïve Bayes classifier using the assign_taxonomy. py command in QIIME (Wang et al., 2007) with reference to the SILVA database accessed from mothur (Schloss et al., 2009) release 119: https://mothur.org/wiki/Silva_reference_files#Release_119. For phylogenetic inference of bacterial and archaeal OTUs, representative sequences for each bacterial OTU were aligned to a SILVA SEED sliced alignment using the PyNAST algorithm (Caporaso et al., 2010) and archaeal and bacterial phylogeny was inferred using FastTree (Price et al., 2010).

This workflow resulted in 33 archaeal, 18,683 bacterial, and 4,194 fungal OTUs. OTUs were further filtered to include only OTUs that occurred in at least 25% of samples across sampling dates (25 archaeal, 10,555 bacterial, and 3,157 fungal OTUs after filtering). This prevalence filter was applied to decrease the effect of spatial variation during field sampling and to decrease the potential for sequencing artifacts that might affect our temporal analysis (see section “Statistical Analysis”). Fungal OTUs were assigned to functional guilds by referencing annotated databases using the open-source software FUNGuild (Nguyen et al., 2016). FUNGuild assigned functional attributes (e.g., trophic mode, guild) to 1,528 fungal OTUs. Fungal OTUs with ‘unknown’ FUNGuild functional annotations (1,629 OTUs) were also included in subsequent downstream analyses. Raw de-multiplexed sequences have been archived in the NCBI Bioproject database and are publicly available under accession no. PRJNA587134.

Statistical Analysis

All statistical analyses were completed using R v 3.5.2 (R Development Core Team, 2010). Differences in soil chemical properties between locations (Hillslope versus Floodplain) were tested using linear-mixed effect models using the R package ‘nlme’ (Pinheiro et al., 2016). Location was the fixed-effect and Plot was the random-effect in the model. The effect of sampling date on the size of the microbial biomass C pool, extractable soil ${NO}_{3}^{-}$ , and gravimetric water content were also assessed using linear-mixed effect models with Time of Sampling as the fixed-effect and Plot as the random-effect in the model. Differences in bacterial and fungal community composition across sampling dates was determined by permutational multivariate analysis of variance (perMANOVA, permutations n = 999) using the R package ‘vegan’ (Oksanen et al., 2011). Dissimilarity between samples was calculated as Bray–Curtis distances for fungi or weighted-unifrac distances for bacteria and archaea (Bray and Curtis, 1957; Lozupone et al., 2007). Pairwise comparisons of community dissimilarity was tested using the R package ‘mctools’ (Leff, 2017).

We identified OTUs that had a statistically significant change in abundance between at least one of three paired time periods [September to March (i.e., ‘Winter’), March to May (i.e., ‘Snowmelt’), May to June (i.e., ‘Spring’)] by calculating log2foldchanges in relative abundance between time-points using the R package ‘gtools’ (Warnes et al., 2018). A 95% confidence interval for the log2foldchange was derived by applying a formula for error propagation for the product of a ratio,

95 % C I_{f o l d c h a n g e} =

| F o l d C h a n g e | \times \sqrt{{(\frac{95 % C I_{t 1}}{m e a n_{t 1}})}^{2} + {(\frac{95 % C I_{t 2}}{m e a n_{t 2}})}^{2}} (1)

where mean_t₁_,t₂ was the mean relative abundance of each OTU across sampling locations on each date and 95% CI _t₁_,t₂ was the 95% confidence interval of each OTU on each date.

P-values were calculated from the confidence intervals (Altman and Bland, 2011) and were adjusted for multiple comparisons by the Benjimini and Hochberg false discovery rate (‘fdr’) method using the R package ‘stats.’ Only OTUs with at least one significant (e.g., FDR-adjusted P < 0.05) log2foldchange across time points were included in further analysis (8 archaeal, 4,583 bacterial, and 1,110 fungal OTUs).

Next, we applied agglomerative hierarchical clustering to group OTUs according to differences in relative abundance patterns across sampling time points (Placella et al., 2012; Evans et al., 2014). Distances between clusters were calculated by the average-linkage method. Dendrograms were visualized using ‘ggdendro’ (de Vries and Ripley, 2016) and cut at even height to categorize OTUs into one of three life strategies for archaea and bacteria. Here we define life strategy to be an assortment of physiological traits, selected for by abiotic and biotic factors, that determine an organism’s acquisition of resources, growth and reproduction, and responses to varying environmental gradients (Chagnon et al., 2013; Ho et al., 2017; Malik et al., 2018). In our study, we assume that changes in relative abundance underlie distinct life strategies related to winter snowpack dynamics and reflect each organism’s degree of adaptation and fitness during winter, snowmelt, and spring. For fungi, we observed a fourth life strategy that had high abundance in fall after plant senescence, but did not analyze this group further. There was substantial taxonomic overlap among strategies at all taxonomic ranks in both the Hillslope and Floodplain. Therefore, we interpreted the ecological significance of each life strategy based on the collective characteristics of a subset of OTUs that accounted for a significant portion of the temporal change in abundance or had large differential abundance across sampling dates.

Lastly, we calculated two measures of phylogenetic community relatedness within life strategies for bacteria and archaea, the net relatedness index (NRI) and the nearest taxon index (NTI), using the R package ‘picante’ with a null model based on random taxa shuffling and 1000 permutations (Webb et al., 2002; Kembel et al., 2010). Phylogenetic clustering (positive NRI/NTI) indicates that taxa in a group are more phylogenetically related than expected, compared to a random sampling from the regional species pool. Phylogenetic clustering suggests a group of taxa that are ecologically similar and share a common niche with traits that have been retained through speciation events (i.e., phylogenetic niche conservatism, Crisp and Cook, 2012). On the other hand, phylogenetic overdispersion (negative NRI/NTI) indicates a group of organisms less phylogenetically related than expected, as compared to the regional species pool. Phylogenetic overdispersion can arise from community assembly processes (e.g., trophic interactions, competitive exclusion) that result in a group of ecologically dissimilar taxa with non-overlapping niches (i.e., niche partitioning) or adaptive traits with a broad phylogenetic distribution.

Results

Soil Characteristics and Winter Snowpack Control Over Soil Temperature and Moisture

Average gravimetric soil moisture content and total soil C were greater on the Floodplain compared to the Hillslope (Table 1). Microbial biomass C in the top-soil (0 to 5 cm depth) was also approximately three times greater on the Floodplain compared to the Hillslope, but was comparable between the Hillslope and Floodplain in soils deeper than 5 cm (Table 1). Across all sampling depths, soils on the Hillslope had consistently greater δ¹³C and δ¹⁵N values than soils in the Floodplain (Table 1).

TABLE 1

Table 1. Soil properties measured in the Hillslope and Floodplain watershed locations.

The onset of winter snowpack accumulation, persistence of snow cover, and complete loss of snow (approximately May 15, 2017) resulted in several different soil temperature and moisture regimes between November 2016 and July 2017 (Table 2 and Figure 1). Soil temperature at 6 cm depth was less than 0°C on the Hillslope in early December prior to the development of a persistent snowpack (Figure 1A). Snow accumulation and persistence coincided with a gradual increase in soil temperature to slightly above 0°C from January throughout the remainder of winter. Complete loss of snow led to higher soil temperatures at both the Hillslope and Floodplain (Figure 1A and Table 2) and trends in soil temperature in spring generally tracked trends in air temperature after snowmelt. We did not observe soil freeze-thaw cycles during or after snowmelt, but did observe soil freeze-thaw cycles in the winter time-period (Supplementary Table S1).

TABLE 2

Table 2. Snowpack effects on soil temperature and moisture.

FIGURE 1

Figure 1. Soil temperature (A) at 6 cm and 17 cm below the soil surface remained above 0°C when soils were covered with snow during winter. Loss of snow cover in May 2017 resulted in a rapid increase in soil temperature. The onset of snowmelt in March triggered a large increase in soil volumetric water content (B), as well as soil water potential (C), which lasted through early June 2017. Volumetric water content was measured at 9-cm depth and water potential was measured at 17-cm depth below the soil surface. Arrows indicate the dates of soil sampling.

Soils in winter generally had lower volumetric water content and lower soil water potential compared to soils collected during snowmelt or spring (Table 2). Because over 90% of the winter snowpack was melted between early March and early May 2017, we refer to that time-period as the snowmelt period (Table 2). Snowmelt resulted in an increase in soil volumetric water content and soil water potential (Figures 1B,C and Table 2) with an extended period of soil saturation beginning more than 60 days prior to the complete loss of winter snow. At the Hillslope, soil moisture declined rapidly after snowmelt as plants broke winter dormancy in spring (Figures 1B,C). In contrast, soils on the Floodplain remained saturated long after snowmelt (Figure 1C and Table 2). In summary, the winter period was characterized by deep snow with generally cold and dry soils, the snowmelt period by the rapid loss of the winter snowpack with cold and wet soils, and the spring soil environment was characterized by rapid warming and soil drying (Table 2 and Figures 1A–C).

Microbial Biomass Bloom and Crash With Subsequent Pulse of Soil N

Soil microbial biomass and the extractable soil nitrate ( ${NO}_{3}^{-}$ ) pool responded strongly to variations in snowpack depth, soil temperature, and moisture (Figures 1, 2A–C). For example, soil microbial biomass C increased 2- to 5-fold at all three soil depths (0 to 5 cm, 5 to 15 cm, 15 cm +) during the snowmelt period on the Hillslope (e.g., March to May, Figure 2B) and trends in microbial biomass pool size were well-correlated with trends in soil water content (Figure 2A). Similar to observations for the Hillslope, microbial biomass C in the Floodplain also increased during snowmelt in the shallowest soils (e.g., 0 to 5 cm, Figure 2B), however, this was not observed in soils sampled more than 5 cm below the soil surface. Soil microbial biomass decreased more than twofold across all depths after snowmelt (between May and June) on the Hillslope (Figure 2B), which coincided with a 3- to 5-fold increase of extractable soil ${NO}_{3}^{-}$ (Figure 2C). Although microbial biomass also collapsed on the Floodplain in topsoils after snowmelt (Figure 2B), in contrast to the Hillslope, we did not observe an increase in the soil ${NO}_{3}^{-}$ pool size after snowmelt in the Floodplain (Figure 2C).

FIGURE 2

Figure 2. Soil water content (A) and soil microbial biomass (B) in the Hillslope and Floodplain. A pulse of extractable soil nitrate (C) was observed in the Hillslope after snowmelt in June 2017. P-values are the outcome of mixed-models testing for the effect of time of sampling on soil water content, microbial biomass, and extractable nitrate.

Similar to changes in the size of soil microbial biomass pool, archaeal and bacterial (Figures 3A,B) as well as fungal community structure (Figures 3C,D) varied significantly during winter, snowmelt, and following the loss of snow in spring. Date of sampling was a significant factor explaining archaeal and bacterial community structure in the Hillslope (perMANOVA pseudo-R² = 0.11, P ≤ 0.001, see Supplementary Table S2 for pair-wise comparisons) and Floodplain (perMANOVA pseudo-R² = 0.33, P ≤ 0.001) (Figures 3A,B). Fungal community structure was likewise significantly affected by the date of soil sampling in both the Hillslope (pseudo-R² = 0.05, P ≤ 0.001) and Floodplain (perMANOVA pseudo-R² = 0.22, P < 0.001, Figures 3C,D). Depth significantly affected archaeal and bacterial (perMANOVA pseudo-R² = 0.33, P < 0.001) as well as fungal community structure (perMANOVA pseudo-R² = 0.05, P < 0.001) in the Hillslope, but not in the Floodplain. We did not observe a significant Date × Depth of sampling interaction for either archaea and bacteria or for fungi in either the Hillslope or Floodplain (perMANOVA P ≥ 0.05 in all cases).

FIGURE 3

Figure 3. Archaeal and bacterial community structure represented by non-metric dimensional scaling (NMDS) in the Hillslope (A) and Floodplain (B) as well as fungal community structure in the Hillslope (C) and Floodplain (D).

Snowmelt Selected for Phylogenetically Clustered Bacterial Life Strategies

Bacterial and archaeal species were grouped together into one of three life strategies based on changes in relative abundance in response to snow accumulation, snowmelt, and the onset of spring (Figures 4A–H). As a group, Winter-Adapted archaea and bacteria had the highest group relative abundance (i.e., all Winter-Adapted OTUs in group summed together, separated by depth) in March (Figures 4G,H), a time when winter snow depth was greatest and soils were relatively cool and dry (Table 2). The relative abundance of Snowmelt-Specialist archaea and bacteria increased 1.8 to 2.4-fold on the Hillslope (Figure 4G) and 2- to 6-fold in the Floodplain (Figure 4H) during the snowmelt period (i.e., March–May). The majority of bacterial and archaeal taxa across life strategies were Spring-Adapted (Table 3). Spring-Adapted archaea and bacteria group relative abundance increased after snowmelt and ranged from 20 to 45% of the total community in the Hillslope or 5 to 8% of the total community in the Floodplain (Figures 4G,H).

FIGURE 4

Figure 4. Archaeal and bacterial OTUs that had a significant change in abundance between any two sampling time points (September, March, May, June) were grouped by hierarchical clustering in the Hillslope (A,C,E) and Floodplain (B,D,F). Three life-history strategies related to winter snow cover, snowmelt, and loss of snow cover were identified and the group relative abundance (i.e., sum of all OTUs in a group at each depth) patterns of each strategy is shown (G,H).

TABLE 3

Table 3. Phylogenetic relatedness of bacteria and archaea grouped by life-history strategy.

Generally speaking, the direction of the phylogenetic relatedness patterns within archaeal and bacterial life strategies were most often toward phylogenetic clustering (positive NRI/NTI values), however, the responses were not always observed at every sampling depth in the Hillslope or Floodplain (Table 3). For example, Spring-Adapted archaea and bacteria were phylogenetically clustered across all soil depths in the Hillslope based on the Net Relatedness Index (NRI, Table 3). In the Floodplain, Snowmelt-Specialist archaea and bacteria were similarly phylogenetically clustered across all sampling depths. In contrast to observations for the Floodplain, Snowmelt-Specialist archaea and bacteria sampled 0 to 5 cm below the soil surface were phylogenetically overdispersed on the Hillslope (NRI = −2.15, Table 3).

Bradyrhizobium (α-Proteobacteria), Tardiphaga (α-Proteobacteria), Sphingomonas (α-Proteobacteria), and Massilia (β-Proteobacteria) collectively accounted for 73% of the total increase in group relative abundance of Winter-Adapted bacteria at the Hillslope during winter (i.e., September 2016 to March 2017; Figure 5A). In the Floodplain, the DA101 soil group (Verrucomicrobia), Thaumarchaeota spp., and Pedobacter (Bacteroidetes) likewise contributed substantially to the increase in relative abundance of Winter-Adapted bacteria (Figure 5B). Streptomyces (Actinobacteria) and Candidatus Nitrososphaera (Thaumarchaeota) together accounted for 11% of the total increase in relative abundance of Snowmelt-Specialist archaea and bacteria from March to May on the Hillslope (Figure 5C), whereas Bacteriovoraceae (δ-Proteobacteria; 44% contribution alone) dominated the response of Snowmelt-Specialists in the Floodplain (Figure 5D). Various phylotypes of Acidobacteria [e.g., Acidobacteriaceae (Subgroup 1) spp., Subgroup 2 spp., Candidatus Solibacter, RB41 spp.] composed the largest fraction of Spring-Adapted bacteria in the Hillslope (Figure 5E). Nitrifying taxa (e.g., Nitrospira spp. and Thaumarchaeota spp.) also contributed 7% to the increase of Spring-Adapted archaea and bacteria from May to June on the Hillslope (Figure 5E). Similar to the Hillslope, other potentially nitrifying organisms belonging to the phylum Nitrospirae (e.g., 4-29 spp.) were major contributors to the increase of Spring-Adapted bacteria after snowmelt and during plant-green up in the Floodplain (Figure 5F). All archaeal and bacterial OTUs with life strategies assigned are listed in Supplementary Table S3.

FIGURE 5

Figure 5. The top 10 archaeal and bacterial phylotypes within each life-history strategy in the Hillslope (A,C,E) and Floodplain (B,D,F). Phylotypes were ranked by their percent change in abundance between the specified sampling time-points (x-axis). The bars represent log2fold changes and the percent contributions to the change in group relative abundance are given at the side of the bar. Abbreviations for taxonomy (phylum or subphylum level) are Acido – Acidobacteria; Actino – Actinobacteria; Bactero – Bacteriodetes; Gemmo – Gemmatimonadetes; Plancto – Planctomycetes; α-Proteo – Alphaproteobacteria; β-Proteo – Betaproteobacteria; δ-Proteo – Deltaproteobacteria; γ-Proteo – Gammproteobacteria; Verruco – Verrucomicrobia.

Fungal Taxa and Functional Groups Were Also Differentiated by Snowmelt Niche

Winter-Adapted, Snowmelt-Specialist, and Spring-Adapted fungi were also observed across all three soil depths in both the Hillslope and Floodplain (Figures 6A–F). Similar to the archaeal and bacterial communities, Winter-Adapted fungi had the highest group abundance (∼8% of total fungal community) at peak winter snow depth (Figure 6G). Snowmelt-Specialist fungi increased in abundance during the snowmelt period (i.e., March to May), and Spring-Adapted fungi had the highest group abundance during plant green-up in June (Figures 6G,H). The majority of fungal taxa in the Floodplain were Spring-Adapted, with far fewer fungal species categorized as either Winter-Adapted or Snowmelt-Specialists (Figures 6B,D,F and Table 4). For example, Spring-Adapted fungi accounted for 27 to 43% of the total fungal community in the Floodplain in June (Figure 6H), whereas Winter-Adapted and Snowmelt-Specialist fungi were rarely more than 2% of the total fungal community at any time.

FIGURE 6

Figure 6. Fungal OTUs that had a significant change in abundance between any two sampling time points were grouped by hierarchical clustering in the Hillslope (A,C,E) and Floodplain (B,D,F). Similar to bacterial OTUs, we observed Winter-Adapted, Snowmelt-Specialists, and Spring-Adapted fungi across all depths and in both the Hillslope (G) and Floodplain (H). A fourth group with high abundance in autumn was also observed, but this group was not included in further analyses.

TABLE 4

Table 4. Composition of functional guilds for fungal life-history strategies related to snowmelt.

Saprotrophy was the dominant fungal trophic mode at both the Hillslope and Floodplain (Table 4). For example, across all soil depths, 80 to 100% of Winter-Adapted fungal taxa were classified as saprotrophic. Symbiotic fungi that associate with plant roots (e.g., arbuscular mycorrhizae, ectomycorrhizae, and root-endophytes) had higher prevalence within the Spring-Adapted life strategy compared to Winter-Adapted or Snowmelt-Specialist fungi (Table 4). Across soil depths, root-associated fungi constituted 16 to 22% of all Spring-Adapted fungi on the Hillslope and 24 to 38% of Spring-Adapted fungi in the Floodplain. Ectomycorrhizal fungi were the predominant root-associated functional group, especially in the Floodplain, where arbuscular mycorrhizal fungi were not observed among any of the life strategies (Table 4).

Genabea (Pezizomycetes), Psilocybe (Agaricomycetes), and Crepidotus (Agaricomycetes) collectively accounted for 77% of the increase in group relative abundance observed for Winter-Adapted fungi occurring from September to March on the Hillslope (Figure 7A). A few select fungal genera also contributed disproportionately to changes in group abundance for Snowmelt-Specialists and Spring-Adapted fungi during and after snowmelt (Figures 7C–F). For example, Pterula (Agaricomycetes) alone accounted for 53% of the increase in relative abundance for Snowmelt-Specialist fungi during the snowmelt period (i.e., March to May) and Tricholoma (Agaricomyetes), Botrytis (Leotiomycetes), and Cuphophyllus (Agaricomycetes) together accounted for 48% of the increase in group abundance for Spring-Adapted fungi from May to June on the Hillslope (Figures 7C,E). Gliomastix (Sordariomycetes), Subplenodomus (Dothideomycetes), Knufia (Ascomycetes familia incertae), Articulospora (Leotiomycetes), and Ganoderma (Agaricomycetes) accounted for 78% of the increase in group relative abundance of Winter-Adapted fungi from September to March in the Floodplain (Figure 7B). Thelephora (Agaricomycetes), Hebeloma (Agaricomycetes), Archaeorhizomyces (Archaeorhizomycetes), and Tetracladium (Leotiomycetes) were Spring-Adapted fungi with substantial increases in relative abundance in the Floodplain after snowmelt (Figure 7F). Very few fungal OTUs employed the Snowmelt-Specialist life strategy in the Floodplain (Table 4), which was composed of only two fungal genera [i.e., Minutisphaera (Dothideomycetes) and Endosporium (Dothideomycetes); Figure 7D]. All fungal OTUs with life strategies assigned are listed in Supplementary Table S4.

FIGURE 7

Figure 7. The top 10 fungal phylotypes within each life-history strategy in the Hillslope (A,C,E) and Floodplain (B,D,F). Phylotypes were ranked by their percent change in abundance between the specified sampling time-points. The bars represent log2fold changes and the percent contributions to the change in group relative abundance are given at the side of the bar. Only two fungal genera contributed to the significant change in fungal Snowmelt-Specialists from March to May in the Floodplain location (D). Abbreviations for taxonomy (class-level) are Agarico – Agaricomycetes; Archaeo – Archaeorhizomycetes; Asco – Ascomycetes; Dothideo – Dothideomycetes; Eurotio – Eurotiomycetes; Leotio – Leotiomycetes; Pezizo – Pezizomycetes; Sordario – Sordariomycetes. Abbreviations for functional groups are Ecto – ectomycorrhizae; Endophyte – root endophyte; Mixo – mixotrophic; Patho – pathotrophic; Sapro – saprotrophic.

Discussion

Soil microbial biomass is commonly observed to bloom beneath the winter snowpack and subsequently decline following snowmelt, often resulting in a significant pulse of soil N in spring (Grogan and Jonasson, 2003; Edwards et al., 2006; Schmidt et al., 2007). In this study, we primarily sought to identify bacteria, archaea, and fungi that were associated with the microbial biomass bloom during winter and biomass crash following snowmelt. We also sought to infer whether the traits that govern microbial community assembly during and after snowmelt were phylogenetically conserved. This was accomplished by measuring the phylogenetic relatedness of bacteria and archaea, or by comparing guilds of soil fungi, that were grouped into life strategies based upon relative abundance patterns before, during, and after snowmelt.

Overall, we found that (1) increases in soil microbial biomass production beneath the winter snowpack were observed at both locations in our study (Hillslope and Floodplain) and appeared to be induced by a pulse of snowmelt infiltration (Figure 2), (2) the microbial biomass collapse was associated with a significant pulse of N measured as extractable soil ${NO}_{3}^{-}$ (Figure 2), (3) three microbial life strategies (Winter-Adapted, Snowmelt-Specialist, and Spring-Adapted) were identified at each soil depth at both locations (Figures 4, 6), (4) life strategies were most often phylogenetically clustered (bacteria and archaea) or shared similar trophic modes (e.g., saprotrophic fungi; Tables 3, 4), and finally that (5) a few select taxa within each life strategy contributed disproportionately to the abundance responses (Figures 5, 7). Thus, we have shown that bacteria, archaea, and fungi with similar responses to snow accumulation and snowmelt are likely to share adaptive traits and we conclude that this framework may be useful in understanding an organism’s snowmelt niche and response to changing winter snowpack conditions.

Role of Moisture and Substrate Utilization in Structuring the Winter-Adapted Niche

Phylogenetic clustering of Winter-Adapted archaea and bacteria may partly be explained by poor connectivity of soil pore spaces due to dry soil conditions that also lead to specific substrate utilization patterns during the winter. Because plant detritus is thought to be the primary substrate utilized by bacteria and fungi during winter (Taylor and Jones, 1990; Uchida et al., 2005; Hooker and Stark, 2012; Isobe et al., 2018), low substrate diffusion rates, combined with high spatial heterogeneity and low resource availability, probably selects for bacteria and fungi with specialized abilities to degrade complex soil organic matter and plant root or leaf litter. A predominance of α-Proteobacteria (e.g., Tardiphaga, Massilia, Phenylobacteria, Rhodoplanes) known to be early- to mid-stage colonizers of decomposing root and leaf-litter (Purahong et al., 2016; Bonanomi et al., 2019) supports this hypothesis (Figure 5). Similar to the bacterial community, we observed that Winter-Adapted fungi were predominantly saprotrophic, with a few select taxa that are known wood or detritus decomposers (e.g., Gliomastix, Psilocybe, Ganoderma) contributing significantly to the increase in Winter-Adapted fungi (Figure 7). Additional fungal taxa commonly found during early- to mid-stages of decomposition of leaf or root litter (e.g., Acremonium, Mycosphaerella, Saccharicola, Tetracladium; Kuhnert et al., 2012) also had significant positive responses to winter snow accumulation at our study site.

Lipson et al. (2009) previously found that bacteria and fungi isolated from beneath the snowpack in winter had a propensity for high growth rates, high mass-specific respiration, and low growth yields. If those observations are generalizable to the winter-adapted taxa observed here, we hypothesize that such physiological traits would select for a winter-adapted phenotype characterized by enhanced investment in extracellular enzyme production to optimize nutrient acquisition via depolymerization of soil organic matter at the expense of traits that optimize stress tolerance and growth yield (Schimel et al., 2007; Malik et al., 2018; Zhalnina et al., 2018).

We also observed two potentially N-fixing genera (e.g., Rhizobium and Bradyrhizobium) that contributed significantly to the overwinter increase in Winter-Adapted archaea and bacteria on the Hillslope, but these genera were not observed in the Floodplain (Figure 5). Recent studies have shown high abundance of N-fixing Bradyrhizobium and Rhizobium in mid- to late-stages of litter decomposition and their presence is thought to alleviate N-limitation that occurs during late stages of decomposition, which thereby facilitates the activity of decomposers that degrade more complex biopolymers (Purahong et al., 2016; Lladó et al., 2017; Bonanomi et al., 2019). In seasonally snow-covered ecosystems, like the East River study site, a similar effect might occur as winter progresses due to the prolonged absence of fresh C inputs from plants aboveground and the potential for progressive N-limitation as the decomposer community propagates. In contrast to the Hillslope, we observed greater total soil N availability in the Floodplain as well as autotrophic, nitrifying archaea (e.g., Thaumarchaeota) during winter. These observations, along with a lack of a significant contribution from N-fixing genera (e.g., Rhizobium and Bradyrhizobium), might indicate that the N-fixing niche is not present beneath the snowpack in the Floodplain.

Snowmelt Infiltration Triggers Microbial Biomass Bloom

The microbial biomass bloom occurred during a period of extended soil saturation that lasted for approximately 65 days and that coincided with net snowpack loss and snowmelt infiltration (Figures 1, 2). Snowmelt-triggered biomass production stands in contrast to previous hypotheses predicting that a sudden decrease in osmolarity due to snowmelt infiltration should cause lysis of microbial cells (Schimel et al., 2007; Jefferies et al., 2010). Our results are similar to those of other studies that have observed an increase in soil microbial biomass production and peak microbial biomass during snowmelt, when soil temperatures were near or above 0°C, followed by a collapse in the microbial biomass pool size when soils became snow-free and were relatively drier in spring (Brooks et al., 1998; Lipson et al., 1999; Edwards et al., 2006; Buckeridge and Grogan, 2010; Isobe et al., 2018). However, we acknowledge that the time intervals of our sampling do not allow us to fully resolve the temporal trends of the microbial biomass bloom and crash or to identify the exact date of peak microbial biomass.

We hypothesize that alleviation from moisture and substrate limitation may have triggered the microbial biomass bloom that we observed under the winter snowpack. While microbial metabolic activity has been shown to occur in soils with temperatures well below 0°C (McMahon et al., 2009), our findings suggest that soil and air temperatures need to be warm enough to induce net snowmelt infiltration and for soil water to be in the liquid phase in order to induce the soil microbial bloom. In addition, the microbial bloom did not occur in Floodplain soils deeper than 5 cm, where higher soil moisture availability existed throughout the winter.

We also observed phylogenetic overdispersion of Snowmelt-Specialist bacteria and archaea on the Hillslope in our study (Table 3), whereas Snowmelt-Specialists were phylogenetically clustered in the Floodplain. For the Hillslope, these results imply that ecologically dissimilar organisms with a broad phylogenetic distribution contribute to the biomass bloom during snowmelt and biomass collapse thereafter. Thus, antagonistic biotic interactions, such as competitive exclusion and density-dependent regulation of population size (e.g., trophic interactions), are likely to be the community assembly mechanisms underlying the microbial biomass bloom, biomass crash, and taxonomic succession that we observed before and after snowmelt. For example, antagonism resulting from toxin production may promote biomass production by reducing competition during a time of high resource availability (Moons et al., 2009; Hibbing et al., 2010). After snowmelt, the population decline on the Hillslope also indicates far greater mortality compared to reproduction rates, which may result from resource depletion or predation by viruses, microfauna, and mesofauna (Moons et al., 2009; Buckeridge et al., 2013; Phaiboun et al., 2015; Georgiou et al., 2017).

Alternatively, competitive exclusion of Winter-Adapted organisms may occur during snowmelt if Snowmelt-Specialist archaea and bacteria have higher growth efficiencies compared to Winter-Adapted taxa (Georgiou et al., 2017). Dissolved organic substrate quality and quantity is known to increase during snowmelt (Campbell et al., 2014; Patel et al., 2018), thus organisms that are characteristic of the biomass bloom (i.e., Snowmelt-Specialists) should be able to re-allocate energy resources away from extracellular enzyme production and nutrient uptake toward protein/fatty acid/DNA synthesis to maximize growth efficiency (Schimel et al., 2007; Malik et al., 2018). Furthermore, it is unlikely that organisms with high growth efficiencies or whole communities with high rates of biomass production would be closely related phylogenetically, but rather could possess genomic features and adaptive traits that are shared due to convergent evolutionary trajectories (Roller et al., 2016; Muscarella and Lennon, 2018).

Spring Nitrogen Dynamics and Spring-Adapted Bacteria and Fungi

We mainly attribute the pulse of extractable soil ${NO}_{3}^{-}$ observed in spring to elevated N mineralization and nitrification rates resulting from a flush of labile C and N that was released from microbial biomass after snowmelt. Lipson et al. (1999) similarly showed a pulse of soil protein that was attributed to the lysis of soil microbial biomass after snowmelt. Furthermore, four ammonia- or nitrite-oxidizing groups [e.g., 0319-6A21 (Nitrospirae), unassigned Nitrospirae spp., and unassigned Thaumarchaeota spp.] contributed substantially to the increase of Spring-Adapted archaea and bacteria on the Hillslope (Figure 5), which also coincided with the pulse of soil ${NO}_{3}^{-}$ in spring (Figure 2) and suggests substantial autotrophic nitrification. It should be noted that nitrogen in the snowpack also likely contributed to the pulse of soil ${NO}_{3}^{-}$ that we observed after snowmelt, however, previous studies have shown that snow is usually a minor contributor compared to nitrogen that has passed through, and been released from, the soil microbial pool (Williams et al., 1996; Williard et al., 2001; Burns and Kendall, 2002; Grogan and Jonasson, 2003; Campbell et al., 2007). In contrast to the Hillslope, we did not observe ${NO}_{3}^{-}$ accumulation in the Floodplain, despite a twofold decline in microbial biomass in the upper soil profile and increases in the abundance of some nitrifying organisms. Lower nitrate accumulation could arise from any combination of the following: lower mineralization and nitrification rates, higher N uptake rates (plant and microbial) with higher N immobilization, or higher rates of ${NO}_{3}^{-}$ assimilation and denitrification at the Floodplain compared to the Hillslope.

In addition to the prevalence of ammonia-oxidizing genera, we also found that four groups of Acidobacteria (i.e., Subgroup 1, Subgroup 2, Candidatus Solibacter, and RB41) made large contributions to the increase in Spring-Adapted bacteria on the Hillslope (Figure 5E). Acidobacteria are traditionally considered to be oligotrophs with abundances negatively correlated to soil organic C (Fierer et al., 2007). However, an increasing number of studies have shown that Acidobacteria responses to gradients in organic C availability can be either positive or negative, challenging their status as strictly oligotrophic (Kielak et al., 2016; Lladó et al., 2017). In addition, some Acidobacteria can be isolated from the environment under anoxic conditions using microbial necromass (i.e., gellan gum) as a growth substrate and Ca. Solibacter has 4- to 6-fold higher nutrient-transporting genes in its genome compared to most other bacteria (Kielak et al., 2016). Therefore, some Acidobacteria phylotypes may occupy the Spring-Adapted niche by capitalizing on the necromass from the biomass crash along with high investment in nutrient transporters after snowmelt.

Mycorrhizal fungi made significant contributions to the increase in Spring-Adapted fungi in both the Hillslope and Floodplain (Figure 7); this pattern was driven primarily by an increase in ectomycorrhizal fungi. For example, in the Floodplain, 6 of the top 10 Spring-Adapted fungi were ectomycorrhizae. If ecosystem N retention in spring is enhanced by well-coupled plant-mycorrhizae phenology, then the implications of our results are noteworthy with respect to changing winter conditions. For example, although earlier snowmelt could result in an earlier pulse of N released from the crash in microbial biomass, floodplain ecosystems may buffer ecosystem N losses via well-coupled plant-mycorrhizal N uptake in spring. It may also be possible to predict the capacity of other watershed locations to retain N in response to earlier snowmelt based upon plant distributions and plant-mycorrhizae associations, given that plant species assemblages and plant-mycorrhizal associations can be mapped and inferred at high spatial resolution using remote-sensing methods (Fisher et al., 2016; Falco et al., 2019).

Conclusion

Because snowmelt rates have been declining in mountainous regions in recent decades (Musselman et al., 2017; Harpold and Brooks, 2018), the length of time during which soils are saturated beneath a melting snowpack may increase in the future and thereby promote a larger snowmelt microbial biomass bloom and crash. Alternatively, a future with warmer air temperatures and a shallower winter snowpack would be expected to result in an increase in the frequency of sustained soil freezing or more frequent soil freeze-thaw events during winter (Hardy et al., 2001; Brown and DeGaetano, 2011; Kreyling and Henry, 2011), which would most likely reduce the magnitude of the microbial biomass bloom, biomass collapse, and pulse of N after snowmelt (Brooks et al., 1998; Buckeridge et al., 2010; Ueda et al., 2013). Whether future environmental conditions will sustain an overwinter microbial bloom warrants further study, because the nutrient flush from microbial biomass following snowmelt can be one of the largest annual soil N pulses in some high-latitude or high-altitude ecosystems (Grogan and Jonasson, 2003; Schmidt et al., 2007).

Here, we have provided novel evidence that the snowmelt period is an environmental filter that differentiates soil bacteria, archaea, and fungi into three distinct life strategies. The phylogenetic relatedness of bacteria and archaea and enrichment of fungal functional groups within life strategies suggests that selective, deterministic processes allow Winter-Adapted, Snowmelt-Specialist, and Spring-Adapted organisms to occupy distinct niches that are related to the winter snowpack and snowmelt. Although our main goal was not to compare and contrast two distinct watershed locations (Hillslope and Floodplain), we did observe differences in the magnitude of the microbial biomass bloom and collapse among locations and soil depths, patterns of ${NO}_{3}^{-}$ accumulation in spring, and differences in contributions from various taxa to each life strategy. These results indicate that local factors, in part, shape the response of these watershed locations to snowmelt. We contend, however, that common microbial traits allow for ecologically similar organisms to occupy the same life strategy, in spite of watershed location. This hypothesis can be tested using high-resolution, culture-independent methods (e.g., genome-resolved metagenomics and metatranscriptomics) and by quantifying traits in controlled laboratory experiments and with isolates to further explore the microbial snowmelt niche.

Data Availability Statement

Sequencing reads can be found in the NCBI Bioproject database accession no. PRJNA587134. Data for this manuscript were downloaded from the Watershed Function SFA data portal (http://wfsfa-data.lbl.gov) and are archived at the U.S. Department of Energy’s Environmental System Science Data Infrastructure for a Virtual Ecosystem (ESS-DIVE, Sorensen et al., 2019).

Author Contributions

PS, HB, and EB designed the study, which is a contribution to the Watershed Function Scientific Focus Area lead by SH and KW. Data collection and analysis was performed by PS (primarily) and SW, AP, MB, UK, NB, HB, and EB. PS, HB, and EB wrote the manuscript. All authors contributed to revising the manuscript.

Funding

RMBL lab spaces and equipment are funded in part by the RMBL Equipment (Understanding Genetic Mechanisms) Grant, DBI-1315705. This material is based upon work supported as part of the Watershed Function Scientific Focus Area funded by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research under Award Number DE-AC02-05CH11231 to Lawrence Berkeley National Laboratory.

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.

Acknowledgments

We would like to thank Erica Dorr, Biz Whitney, Hans Wu Singh, and Chelsea Wilmer for their assistance in the lab and in the field. Mark Conrad, Rosemary Carroll, and Wendy Brown were extremely helpful digging snow pits during winter. We also wish to thank Jenny Reithel and Shannon Sprott at Rocky Mountain Biological Laboratory (RMBL) for their help at the field site. Lastly, we thank Zarine Kakalia and Joan Damerow for their help to archive data on ESS-Dive.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.00871/full#supplementary-material

TABLE S1 | Daily soil freeze-thaw cycles measured in the Hillslope and Floodplain during winter and before and after snowmelt in 2016/17. We considered a freeze-thaw cycle to have occurred if the daily maximum soil temperature was > 0°C and the daily minimum soil temperature was < 0°C.

TABLE S2 | Pair-wise comparisons of permANOVAs for bacteria and fungi.

TABLE S3 | Life strategies assigned to archaeal and bacterial OTUs.

TABLE S4 | Life strategies assigned to fungal OTUs.

References

Aanderud, Z. T., Jones, S. E., Schoolmaster, D. R. Jr., Fierer, N., and Lennon, J. T. (2013). Sensitivity of soil respiration and microbial communities to altered snowfall. Soil Biol. Biochem. 57, 217–227.