Microbial Community Structure of Subalpine Snow in the Sierra Nevada, California

ABSTRCT

Mounting evidence suggests that Earth’s cryosphere harbors diverse and active microbial communities. However, our understanding of microbial composition and diversity in seasonal snowpack of montane ecosystems remains limited. We sequenced the 16S rRNA gene to determine microbial structure (composition and diversity) of snow at two depths (0–15 and 15–30 cm) of a subalpine site in the Southern Sierra Critical Zone Observatory, California, U.S.A. Proteobacteria dominated both depths (~72% of sequences), and this phylum was composed mostly of bacteria within the Rhodospirillales order. Cyanobacteria were almost exclusively present in the upper snow layer, while Actinobacteria and Firmicutes were more abundant in the deep snow layer. Many of the most abundant phylotypes were Acetobacteraceae. Phylotype NCR4874, which comprised 22%–32% of the sequences, was most closely related to the N₂-fixing bacteria Asaia siamensis, suggesting that N₂ fixation may be an important process within the Sierra snowpack. In addition, just under half (45%) of the numerically dominant phylotypes shared >98% similarity with sequences recovered from other cold environments. Our results suggest that microbial communities of subalpine Sierra Nevada snowpack are diverse, with both snow depths harboring distinct but overlapping communities consisting largely of cold-adapted bacteria.

Introduction

The cryosphere is a vital, but still incompletely understood, portion of Earth that includes sea ice, glaciers, permafrost, and seasonally snow-covered land (Miteva, 2008). Observations suggest that seasonal snowpack, the largest component of the cryosphere covering over a third of the land surface (Hinkler et al., 2008), can influence the local climate (Cohen and Entekhabi, 2001; Déry and Brown, 2007; Qian et al., 2014), underlying soil (Buckeridge and Grogan, 2008; Schimel et al., 2004), and adjacent ecosystems such as streams (Brown et al., 2007). For example, by regulating freeze-thaw events, the extent and duration of snow cover can affect soil microbial community composition (Zinger et al., 2009), microbial-mediated soil nitrogen (N) cycling (Edwards et al., 2007), and greenhouse gas exchange with the atmosphere (Jones, 1999). However, the degree to which snow acts as an ecosystem itself has only recently been studied, and much of the biology and ecology of snow remains unknown (Larose et al., 2013a).

Despite conditions that include desiccation, freeze-thaw cycling, high solar irradiation, and low nutrient concentrations (Miteva, 2008), snow can harbor diverse and active microbial communities (Lopatina et al., 2013; Margesin and Miteva, 2011; Møller et al., 2013) that are deposited from both local and distant sources (Chuvochina et al., 2011a; Harding et al., 2011). For example, Chuvochina et al. (2011a) found that microorganisms associated with Saharan dust collect on and likely colonize snow in the Alps. In addition, using a differential staining method, Lopatina et al. (2013) revealed that up to 45.8% of bacteria were viable in snow samples collected from Antarctica. Many of these microorganisms can serve as biological ice nuclei, catalyzing ice formation and snowfall events (Christner et al., 2008)—and are able to cope with the potentially harsh conditions associated with snow by increasing membrane fluidity, synthesizing “antifreeze” proteins, and producing detoxifying enzymes (Casanueva et al., 2010; De Maayer et al., 2014; Maccario et al., 2014; Rothschild and Mancinelli, 2001).

Within snowpack, changes in physicochemical properties (e.g., ultraviolet [UV] light and snow density) with depth may stimulate postdepositional selection and subsequent vertical stratification of microbial communities (Møller et al., 2013; Seok et al., 2009; Xiang et al., 2009a, 2009b). For example, Møller et al. (2013) found that Cyanobacteria in High Arctic snow were confined to the top 90 cm of snow and speculated that the absence of this phylum from deeper snow layers was due to reduced light availability. Xiang et al. (2009a) reported similar findings from glacial snow in the Tianshan Mountains. However, the question of how microorganisms respond to burial in subalpine montane snowpack remains relatively unexplored, as much of our current understanding of the abundance, viability, and composition of microbial communities within snow comes from studies of polar environments (Harding et al., 2011; Hell et al., 2013; Larose et al., 2010; Lopatina et al., 2013; Maccario et al., 2014; Møller et al., 2013).

Montane snowpack differs significantly from high latitude and glacial snowpack because (sub)alpine environments experience a comparatively milder climate, receive more direct solar radiation, and the snowpack itself is more ephemeral. In addition, montane environments likely receive microbial inputs from different source locations. The few available studies of montane snow-associated microbial communities have revealed microbial densities ranging from 680 to 720,000 cells mL^-1 snow (Bauer et al., 2002; Lazzaro et al., 2015; Liu et al., 2009), with a predominance of bacteria in the Proteobacteria and Actinobacteria phyla (Chuvochina et al., 2011b; Liu et al., 2007; Meola et al., 2015). In the Sierra Nevada, California, bacterial abundances were found to range from 9.8 × 10⁴ to 5.5 × 10⁵ cells mL^-1 in white (no algae) and red (algae present) snow areas, respectively (Thomas and Duval, 1995); however, bacterial community composition was not assessed. Segawa et al. (2005) further demonstrated the ability of psychrophilic and psychrotrophic bacterial species (e.g., Clostridium psychrophilum) to grow in a montane snow environment, confirming that snow can be a viable habitat for microbial life, and not simply a reservoir of deposited, dormant cells.

Because microbial communities within the cryosphere may alter surface albedo (Thomas and Duval, 1995; Yallop et al., 2012) and regulate nutrient cycling (Amoroso et al., 2010; Hodson et al., 2008), characterizing snow-associated microbial communities of subalpine systems is important for understanding current and future dynamics of snow in these montane ecosystems. In particular, measuring microbial richness and composition of subalpine snowpack can provide information on the role of snow as a vector for microbial dispersal, the role of ephemeral snow in the maintenance of soil microbial diversity, and the potential for microbial-mediated nutrient transformations to occur in situ (Maccario et al., 2014).

The overall objective of this exploratory study was to assess microbial structure (composition and diversity) in different layers (0–15 and 15–30 cm depth) of subalpine snow in the Sierra Nevada, California, U.S.A. Using a next-generation sequencing approach (Illumina MiSeq), we expected to find a diverse microbial community composed, in part, of stress-tolerant and putative cold-adapted taxa. We further expected to find shifts in composition by snow depth that would be indicative of differences in predepositional microbial dispersal or postdepositional selection.

Methods

Site Characteristics and Sample Collection

Samples were collected from Shorthair Creek (37.06698°N, 118.98711°W), the highest elevation site of the Southern Sierra Critical Zone Observatory (Fig. 1; http://criticalzone.org/sierra/). This subalpine site is located at 2700 m elevation along the western slope of the Sierra Nevada and is composed primarily of lodgepole pine (Pinus contorta) and western white pine (Pinus monticola). Annual temperature ranges from -1.9 °C (mean minimum) to 10.2 °C (mean maximum), and mean annual precipitation is 1078 mm y^-1 (Goulden et al., 2012).

On 12 October 2014, prior to snowfall, a 3 × 2.7 m Tyvek tarp (DuPont HomeWrap, E. I. du Pont de Nemours and Company, Wilmington, Delaware, U.S.A.) was secured to the ground with polyvinyl chloride (PVC) pipes and sandbags, and marked with t-posts in order to provide an area where snow-deposited microbial communities would not be influenced by local soil. On 2 February 2015, we collected snow samples from the snowpack that had accumulated on the tarp-covered area. The snowpack was approximately 30 cm deep and had accumulated primarily from a series of small snow events in late 2014 followed by a period of interstorm snow ripening (Fig. 2). The only snow event within the month preceding our sampling campaign resulted in less than 1 cm of precipitation. We therefore sampled at two different depths that largely reflected different postdeposition ripening rather than different snow events: A single sample (13 × 13 cm in area) was taken from the upper layer (0–15 cm depth) and underlying deeper layer of snow (15–30 cm depth; total snow depth was 30 cm). After digging a pit adjacent to the collection site, snow was removed using a saw and placed in a 5.4 L sterile Whirlpak bag (Nasco Inc., Salida, California, U.S.A.). To minimize the potential for contamination, we sterilized all tools with 10% bleach followed by 70% ethanol and wore Tyvek body suits (DuPont Coveralls) and latex gloves during sampling. The snow samples were placed on dry ice while transported back to the laboratory, where they were stored at -20 °C until analysis.

Sample Processing

We left the two samples to melt in a cold room at 4 °C prior to filtering them through sterile 0.2 μm filters (47-mm diameter; Pall Supor 200 Sterile Grid filters, Pall Corporation, Port Washington, New York, U.S.A.) in an acid-wash sterilized vacuum filtration unit. Once filtration was complete, filters were stored in sterile Whirlpak bags at -20 °C until further analysis. To test for contamination, we created a negative control by filtering 1 L sterile Milli-Q water through a blank filter, and processed that control alongside the snow samples. After the samples were extracted and amplified, the blank filter was free of DNA, indicating that contamination did not occur during processing.

We extracted microbial DNA from the filters using a MO BIO PowerWater DNA Isolation Kit (MO BIO Laboratories, Carlsbad, California, U.S.A.) following the manufacturer's instructions. The extracted DNA was then quantified using a NanoDrop 2000 (Thermo Fisher Scientific, Wilmington, Delaware, U.S.A.) and amplified in duplicate using primers targeting the V3-V4 region of the 16S rRNA gene (S-D-Bact-0341-b-S-17 and S-D-Bact-0785-a-A-21; Klindworth et al., 2013). We conducted polymerase chain reaction (PCR) by combining 2.5 μL DNA template, 5 μL each of 1 μM forward and reverse primers, and 12.5 μL KAPA HiFi HotStart ReadyMix (KAPA Biosystems, Wilmington, Massachusetts, U.S.A.), totaling a 25 μL reaction. Thermocycler conditions were as follows: 95 °C for 3 minutes, followed by 25 cycles of 95 °C for 30 seconds, 55 °C for 30 seconds, 72 °C for 30 seconds, followed by an extension step for 5 minutes at 72 °C. After amplification, we combined and purified the duplicate PCR products using Agencourt AMPure XP Beads (Beckman Coulter Genomics, Danvers, Massachusetts, U.S.A.). A second round of PCR was subsequently conducted to attach dual indices and sequencing adapters to each sample using the Nextera XT Index Kit (Illumina, San Diego, California, U.S.A.). Briefly, 5 μL DNA, 5 μL each of 1 μM forward and reverse index primers, 25 μL KAPA HiFi HotStart ReadyMix, and 10 μL PCR grade water were combined to create a 50 μL mixture. Thermocycler conditions were as follows: 95 °C for 3 minutes, followed by 8 cycles of 95 °C for 30 seconds, 55 °C for 30 seconds, 72 °C for 30 seconds, followed by an extension step for 5 minutes at 72 °C. We then conducted a second purification step (same as described above) on the indexed amplicons and quantified the products using the Quant-iT PicoGreen dsDNA assay kit (Life Technologies, Grand Island, New York, U.S.A.). As a final step, we pooled the samples together in equimolar concentrations with samples from another project and sequenced them using an Illumina MiSeq instrument at the University of California Riverside (UCR) Genomics Core Facility.

FIGURE 1. Map showing location of Shorthair Creek study site (37.06698°N, 118.98711°W) along an elevation transect through the Sierra Nevada, California. (left) The elevation of the site, at ~2.7 km, puts it well within the limits of seasonal snowpack, and thus provides an opportunity to investigate the composition and diversity of microorganisms deposited on the ecosystem by snow.

FIGURE 2. Maximum, minimum, and average daily air temperature (°C), and (bottom) daily precipitation and snow water equivalent (SWE) in centimeters from 1 October 2014 through 28 February 2015. Data were obtained from the Department of Water Resources California Data Exchange Center from two locations near Shorthair Creek along the western slope of the Sierra Nevada. Temperature and SWE data were obtained from West Woodchuck Meadow (2773 m elevation; 37.030°N, 118.918°W); precipitation data were obtained from the Dinkey Forest Service site (1728 m elevation; 37.066°N, 119.039°W). Samples were collected on 2 February 2015, which is highlighted on the bottom portion of the figure with an arrow.

Sequence Analysis

We obtained the sequences already demultiplexed from the UCR Genomics Core Facility and subsequently processed them using Quantitative Insights into Microbial Ecology (QIIME; Caporaso et al., 2010). After we joined the forward and reverse reads (allowing for 20% maximum difference within the region of overlap), we used default parameters to conduct quality control. Specifically, reads were excluded if less than 75% of the read length was consecutive high-quality base calls, if there were more than three consecutive low-quality base calls, if a Phred score was below three, if the length was less than 75 bases, or if one or more ambiguous calls were present (Bokulich et al., 2013). We picked operational taxonomic units (OTU or phylotype) at 97% similarity using open reference UCLUST against the 13_8 release of the Greengenes database. Reads that did not match any sequences in the database were clustered de novo. Raw sequences have been deposited to the NCBI Sequence Read Archive under the study accession number SRP068775.

In QIIME, we removed chloroplast and unassigned sequences prior to subsampling to an even depth of 24,678 reads. Singletons were not removed. Using these rarefied OTU tables, we then estimated relative taxon abundances and ran alpha diversity analyses for each depth. Alpha diversity was measured using a variety of complementary methods, including Chao1 richness, Shannon diversity index, Simpson's diversity, and Faith's phylogenetic diversity (Faith, 1992).

Organisms closely related to numerically dominant phylotypes (>200 sequences) were identified using the NCBI Basic Local Alignment Search Tool (BLASTn) and EzTaxon (Johnson et al., 2008; Kim et al., 2012a). The EzTaxon server uses BLASTn and a pair-wise sequence alignment algorithm to identify both validly published prokaryotic names and representatives of uncultured phylotypes. In our study, we used the EzTaxon server to filter for closely related, validly published (isolated) prokaryotic names, providing information on the closest related microorganism. We used BLASTn to identify closely related uncultured/environmental sample sequences, providing insight into the environments that harbor similar phylotypes.

Results

Sequencing resulted in 24,678 reads for the shallow snow layer (0–15 cm) and 49,336 reads for the deep snow layer (15–30 cm). At an even sequencing depth of 24,678 reads, our efforts yielded 1,173 and 1,861 OTUs (phylotypes) from the shallow and deep snow layers, respectively (Table 1). These phylotypes, 674 of which were shared between the two snow depths (i.e., recovered from both depths; data not shown), belonged to 22 phyla, 67 classes, 123 orders, 223 families, and 373 genera. Both samples were dominated by Proteobacteria, which accounted for approximately 72% of the sequences (Fig. 3; Supplementary Table 1). Acidobacteria were less prevalent, but the relative abundance of this phylum was also similar across both samples (~7%). In contrast, the relative abundances of Cyanobacteria, Chloroflexi, Firmicutes, and Actinobacteria varied by snow depth (Fig. 3). Cyanobacteria accounted for 7.0% of sequences in the shallow snow layer, while only 0.01% of sequences were Cyanobacteria in the deep snow layer. We observed a similar trend with Chloroflexi (5.5% and 1.1% in shallow and deep snow layers, respectively). In contrast, the relative abundance of Actinobacteria was greater in the deep snow layer (13.5%) than the shallow snow layer (2.3%), as was the relative abundance of Firmicutes (0.3% and 3.6% in shallow and deep snow layers, respectively).

TABLE 1 Values of alpha diversity by snow depth (0–15 cm and 15–30 cm) from the Shorthair Creek study site.

FIGURE 3. Relative abundances of the 22 phyla recovered between both snow depths (0–15 and 15–30 cm) from the Shorthair Creek study site. Both graph and legend share the same order, sequential from top to bottom. Proteobacteria dominated both samples, while Cyanobacteria were more abundant in the upper snow layer and Actinobacteria and Firmicutes were more abundant in the deeper snow layer.

Within Proteobacteria, the order Rhodospirillales was the most prevalent (Fig. 4; 70.3%–80.9% of sequences). Of the Rhodospirillales, 99% were Acetobacteraceae. The order Myxococcales was the next most prevalent, accounting for 5.2% (shallow snow) and 5.9% (deep snow) of sequences within this phylum. While Burkholderiales made up 18.0% of Proteobacteria sequences in the upper snow layer, this order made up only 2.5% in the deep snow layer. Twenty-six other orders of the Proteobacteria were also present, including Legionellales, Sphingomaonadales, Rhizobiales, and Pseudomonadales (Fig. 4; Supplementary Table 2).

FIGURE 4. Relative abundances of each order within the Proteobacteria phylum recovered from each snow depth (0–15 and 15–30 cm) at the Shorthair Creek site. Both graph and legend share the same order, sequential from top to bottom. Rhodospirillales accounted for the majority of Proteobacteria sequences in both samples. “Other” indicates the combined relative sequence abundance of the 17 additional, rare taxonomic Proteobacteria orders.

We detected a number of genera that are known to tolerate and thrive in conditions within the cryosphere. These genera included the following: Methylobacterium (0.05% and 0.32% in shallow and deep layers, respectively), Geodermatophilus (0% and 2.24% in shallow and deep layers, respectively), and Hymenobacter (2.03% and 0.55% in shallow and deep layers, respectively). Janthinobacterium (Kim et al., 2012b; Segawa et al., 2005), a genus that contains psychrotolerant bacteria, was also present in both samples (0.63% and 0.15% in shallow and deep layers, respectively).

Of the most abundant phylotypes (29 phylotypes identified as those that occurred more than 200 times total across both samples), 11 were from the Acetobacteraceae family and 25 were shared between the shallow and deep snow layers (Fig. 5). The most predominant phylotype across both snow depths (phylotype NCR4874, 22%–32% of sequences) was most closely related (96.8%) to Asaia siamensis, a dinitrogen (N2) fixing Acetobacteraceae. The top uncultured BLASTn hit for this phylotype was with a human-associated bacterium recovered by Kong et al. (2012), but it also shared 98% similarity with sequences recovered from cold environments such as snow in the Alps, glacial ice worms, and Arctic lichens (data not shown). While it is therefore possible that this phylotype (NCR4874) was a contaminant, this possibility is unlikely given that the blank control consistently tested negative for contaminant DNA. Thirteen additional abundant phylotypes matched at ≥⃒97% to sequences found in other cold environments. For example, phylotype NCR5142, which was common in both snow layers, shared 99% similarity with a sequence previously identified from peat soil in the Arctic. However, not all dominant phylotypes occurred in both snow layers (Fig. 5). Two phylotypes that were common in the shallow snow layer but were not identified in the deep snow layer belonged to Cyanobacteria (phylotypes 818835 and 558201) and one belonged to Chloroflexi (phylotype NCR6630). Many phylotypes that occurred more frequently in the bottom snow layer than the top were Actinobacteria, one of which shared >99% similarity with Geodermatophilus siccatus (phylotype 453616) and another with Galbitalea soli (phylotype 279325).

In addition to harboring differences in microbial community composition, the snow samples also differed in alpha diversity. The shallow snow layer had fewer observed OTUs than the deep snow layer (1173 and 1839 OTUs, respectively). The same was true for Chao1 richness estimator, Faith's phylogenetic diversity, and Shannon diversity, while Simpson's diversity was more equitable across snow depth (Table 1).

Discussion

Mounting evidence suggests that the cryosphere harbors diverse and active microbial communities (Margesin and Miteva, 2011). However, our understanding of how microorganisms vary by snow depth in seasonally dynamic montane ecosystems remains limited (Møller et al., 2013). We used high throughput sequencing to determine microbial community composition from two depths in a Sierra Nevada subalpine site (California, U.S.A.) and found that community membership between the shallow and deep snow layers was distinct but overlapping. While lack of replication limits our ability to generalize this finding, it suggests that vertical spatial variation in microbial community composition may occur within ephemeral snowpack just as it does within other terrestrial and aquatic environments (Eilers et al., 2012; Humayoun et al., 2003).

Proteobacteria dominated both snow depths, which was expected given that this phylum is widespread in cold environments and has been found to predominate snowpack in Greenland (Møller et al., 2013), Svalbard (Amato et al., 2007), Antarctica (Lopatina et al., 2013), Alaska (Sheik et al., 2015), and Mount Everest (Liu et al., 2007). In our study, the Rhodospirillales order (α-Proteobacteria) comprised more than 70% of the Proteobacteria sequences, and 11 of the most abundant phylotypes were Acetobacteraceae (acetic acid bacteria, a family of Rhodospirillales). Many Acetobacteraceae are common symbionts of insects (Crotti et al., 2010), and studies have detected this family in eolian dust samples (Bowers et al., 2012); these two dispersal mechanisms could have contributed to the high rates of Acetobacteraceae colonization at our site (Gilbert, 1980). Another possibility is that these airborne microorganisms acted as ice nuclei for the snow itself (Christner et al., 2008).

The high relative abundance of Acetobacteraceae in both snow depths not only illustrates their ability to withstand the relatively harsh conditions of snow (i.e., low temperatures, physical stress from freeze-thawing cycles, etc.), but also suggests that they are likely ecologically and functionally significant within this environment. Species in this family perform a number of critical functions, including N₂ fixation and phosphate solubilization (Loganathan, 2004), and many are important endophytes (Saravanan et al., 2008). The most abundant phylotype (phylotype NCR4874) within both samples shared 96.8% similarity with A. siamensis, an aerobic Acetobacteraceae that performs N₂ fixation (Katsura et al., 2001; Samaddar et al., 2011). While N₂ fixation appears to be a relatively unconserved microbial trait (genus to species level; Martiny et al., 2015), the prevalence of this phylotype and its similarity with A. siamensis suggests that there may be the genetic potential for N₂ fixation in this environment.

FIGURE 5. Heatmap showing the numerically dominant phylotypes (>200 sequences) across both snow depths (0–15 and 15–30 cm) and their relative abundances as percentages of all 16S rRNA sequences recovered within their respective depth. The color spectrum ranges from dark blue (the lowest relative abundance) to dark red (the highest relative abundance). OTU ID numbers are also listed (NCR = new.cleanup.referenceOTU as designated during the OTU picking process). Also included are the (1) taxonomic assignments to the family level, (2) the closest related NCBI BLAST hit and the environment from which it was recovered; and (3) the closest related described (isolated) microorganism, environment from which it was isolated, and the associated reference. OTU NRC4874, an Acetobacteraceae, was the most abundant phylotype across both samples, accounting for 22.61% and 31.78% of sequences in the upper and deeper snow layers, respectively. % sim = % sequence similarity, Proteo = Proteobacteria; Actino = Actinobacteria; Acido = Acidobacteria. Alpha, Beta, Delta, Gamma denote classes of Proteobacteria. Full citations can be found in the Appendix.

Previous studies suggested that snow-associated microbial communities are largely derived from eolian movement and deposition of microorganisms (Bowers et al., 2012; Chuvochina et al., 2011b). Windborne dust is one primary vector of such movement, although aerosolized microorganisms do not necessarily require association with dust particles. The composition of eolian microbial communities is influenced by provenance (source location) and temporal trends in microbial composition of the source material (Bowers et al., 2012, 2011). In the Sierra Nevada, westerly winds move air from Asia across the Pacific Ocean; accordingly, dust inputs to the Sierra Nevada originate from a combination of local (e.g., Central Valley, California, U.S.A.) and distant (e.g., Asia) sources (Kellogg and Griffin, 2006; Vicars and Sickman, 2011), depending on the elevation and time of year (Aciego et al., unpublished). While it is possible that the microorganisms observed in our study are cosmopolitan (i.e., they are found everywhere), mounting evidence suggests that many microorganisms display biogeographic patterns (Martiny et al., 2006; Papke and Ward, 2004). Hence, the presence of particular phylotypes within the snow can provide insights into which regions or sources are contributing to microbial colonization over time (Bowers et al., 2012). For example, at the Shorthair Creek site, one of the most abundant phylotypes was closely related (>99%) to a sequence previously recovered from the Sahara (del Carmen Montero-Calasanz et al., 2013), a desert that is one of the greatest sources of dust globally (Goudie and Middleton, 2001). Chuvochina et al. (2011a) also demonstrated that Saharan dust significantly influenced microbial composition of snowpack in the Alps. Other abundant phylotypes shared >98.5% similarity with microorganisms isolated from the air (phylotype 539915) and soil (phylotype 279325) in Korea, which is a known source of dust for the western United States (Kellogg and Griffin, 2006). The presence of these phylotypes suggests that part of the microbial community in the subalpine Sierra Nevada snowpack is derived from distant desert sources (Griffin, 2007).

While differences in the relative abundance of microbial groups between snow depths may reflect changes in the place of origin, temporal trends in microbial composition of the source material, or redistribution of microbial communities within snowpack (Lazzaro et al., 2015), they may also result from postdepositional environmental selection (Xiang et al., 2009b). In our study, both the upper and lower snow layers likely were deposited from similar sources (Pacific storm events from Asia), but postdepositional time differed; therefore, it is conceivable that postdepositional selection and/or redistribution within the snowpack played a larger role than place of origin in developing vertical stratification. For instance, the elevated relative abundance of phylotypes 1110021 and 453616 within the deeper layer of snow may reflect the ability of Geodermaphilaceae to tolerate stressful conditions (Ivanova et al., 2010). Indeed, phylotype 1110021 shared 99.3% similarity with Modestobacter marinus, an isolate that has often been recovered from harsh environments such as deserts and is able to withstand high gamma radiation, UV levels, and desiccation typical of snowpack conditions (Gtari et al., 2012). Similar tolerance for extended periods of extreme desiccation may explain the higher abundance of Actinobacteria and Firmicutes in the deeper snow layer (Connon et al., 2007), and may also account for the commonness of these phyla in other snow-dominated environments (Liu et al., 2007). It is also possible that environmental conditions selected for Cyanobacteria (a phylum of photosynthetic and N2-fixing bacteria) in the shallow snow layer but not the deeper layer, thereby shaping the vertical distribution of this phylum within the snowpack. A similar inverse relationship between the abundance of Cyanobacteria and snow depth was reported by Møller et al. (2013), which led those authors to conclude that light intensity controls the spatial pattern of these photosynthetic bacteria.

Thirteen of the dominant phylotypes matched with >98% similarity to sequences previously recovered from cold environments such as glaciers, polar deserts, and the arctic tundra. Phylotype NCR652, for example, was closely related (99%) to a lichen-associated microorganism found within the polar regions. In addition, we recovered four lesser abundant genera (Methylobacterium, Geodermatophilus, Hymenobacter, and Janthinobacterium) that are known to survive in the cryosphere (Chuvochina et al., 2011a; Kim et al., 2012b; Rainey et al., 2005). The first three genera can tolerate high doses of ionizing radiation (Chuvochina et al., 2011b; Rainey et al., 2005), while the fourth is a genus that contains psychrotolerant bacteria (Segawa et al., 2005). Burkholderia, a genus comprised of bacteria found in polar environments (Lee et al., 2014) and that dominate some subalpine soils (Lipson et al., 2008), was also found in the snowpack, albeit at a low relative abundance (0.15%–0.45% of sequences in upper and lower layers of snow, respectively). The detection of these sequences within our study and others confirms prior findings that cold-adapted microorganisms can disperse globally, and that many are able to survive and persist in geographically diverse cold environments.

With only two samples collected (and only one per depth), lack of replication is a clear limitation of this study. However, our exploratory approach and results can be used to develop testable hypotheses and highlight several research questions that warrant further investigation in the California Sierra Nevada and other (sub)alpine ecosystems. Some notable questions include: Is there genetic potential for N₂ fixation in ephemeral montane snowpack, similar to what has been found in the Arctic (Telling et al., 2011; Larose et al., 2013b)? And, if there is genetic potential, does this result in ecologically meaningful rates of in situ N₂ fixation? Addressing these and other questions regarding the processes that govern vertical variability in microbial composition will provide critical information on microbial structure and function in this potentially important, but understudied, component of montane ecosystems.

Appendix

TABLE A1 (A) Relative abundances at the phyla level for each snow depth. (B) Richness (number of operational taxonomic units [OTUs]) of the top 10 most abundant phyla.

TABLE A2 Relative abundances of each order comprising the Proteobacteria for each snow depth

TABLE A3 References for Figure 5.

Table

Acknowledgments

We would like to thank Erin Stacy and Patrick Womble for their help with sample collection. We would also like to thank two anonymous reviewers for providing constructive comments. This research is part of the Southern Sierra Critical Zone Observatory, which is funded by the National Science Foundation (EAR 1331939). Sample collection was funded in part by the National Science Foundation Rapid Response Research (RAPID) program (1449197), and partial funding was provided by a Packard Fellowship to SMA.