Multi-omics of permafrost, active layer and thermokarst bog soil microbiomes.

Over 20% of Earth's terrestrial surface is underlain by permafrost with vast stores of carbon that, once thawed, may represent the largest future transfer of carbon from the biosphere to the atmosphere. This process is largely dependent on microbial responses, but we know little about microbial activity in intact, let alone in thawing, permafrost. Molecular approaches have recently revealed the identities and functional gene composition of microorganisms in some permafrost soils and a rapid shift in functional gene composition during short-term thaw experiments. However, the fate of permafrost carbon depends on climatic, hydrological and microbial responses to thaw at decadal scales. Here we use the combination of several molecular 'omics' approaches to determine the phylogenetic composition of the microbial communities, including several draft genomes of novel species, their functional potential and activity in soils representing different states of thaw: intact permafrost, seasonally thawed active layer and thermokarst bog. The multi-omics strategy reveals a good correlation of process rates to omics data for dominant processes, such as methanogenesis in the bog, as well as novel survival strategies for potentially active microbes in permafrost.

Permafrost microbial communities and functional genes are structured by latitudinal and soil geochemical gradients.

Permafrost underlies approximately one quarter of Northern Hemisphere terrestrial surfaces and contains 25-50% of the global soil carbon (C) pool. Permafrost soils and the C stocks within are vulnerable to ongoing and future projected climate warming. The biogeography of microbial communities inhabiting permafrost has not been examined beyond a small number of sites focused on local-scale variation. Permafrost is different from other soils. Perennially frozen conditions in permafrost dictate that microbial communities do not turn over quickly, thus possibly providing strong linkages to past environments. Thus, the factors structuring the composition and function of microbial communities may differ from patterns observed in other terrestrial environments. Here, we analyzed 133 permafrost metagenomes from North America, Europe, and Asia. Permafrost biodiversity and taxonomic distribution varied in relation to pH, latitude and soil depth. The distribution of genes differed by latitude, soil depth, age, and pH. Genes that were the most highly variable across all sites were associated with energy metabolism and C-assimilation. Specifically, methanogenesis, fermentation, nitrate reduction, and replenishment of citric acid cycle intermediates. This suggests that adaptations to energy acquisition and substrate availability are among some of the strongest selective pressures shaping permafrost microbial communities. The spatial variation in metabolic potential has primed communities for specific biogeochemical processes as soils thaw due to climate change, which could cause regional- to global- scale variation in C and nitrogen processing and greenhouse gas emissions.

Phylogenetic and ecological analyses of soil and sporocarp DNA sequences reveal high diversity and strong habitat partitioning in the boreal ectomycorrhizal genus Russula (Russulales; Basidiomycota).

SUMMARY: *Although critical for the functioning of ecosystems, fungi are poorly known in high-latitude regions. Here, we provide the first genetic diversity assessment of one of the most diverse and abundant ectomycorrhizal genera in Alaska: Russula. *We analyzed internal transcribed spacer rDNA sequences from sporocarps and soil samples using phylogenetic methods, operational taxonomic unit (OTU) delimitations and ordinations to compare species composition in various types of boreal forest. *The genus Russula is highly diverse in Alaska, with at least 42 nonsingleton OTUs (soil) and 50 phylogroups (soil + sporocarp). Russula taxa showed strong habitat preference to one of the two major forest types in the sampled regions (black spruce and birch-aspen-white spruce), and some preference for soil horizon. *Our results show that the vast majority of Russula species are present in the soil samples, although some additional taxa are expected to be found with extended sampling. OTU diversity in black spruce forests was only one-third of the diversity observed in mixed upland forests. Our findings suggest that some of the diversity is niche based, especially along host and successional axes, because most OTUs predictably occurred in specific habitats, regardless of geographical location.

Life at the Frozen Limit: Microbial Carbon Metabolism Across a Late Pleistocene Permafrost Chronosequence.

Permafrost is an extreme habitat yet it hosts microbial populations that remain active over millennia. Using permafrost collected from a Pleistocene chronosequence (19 to 33 ka), we hypothesized that the functional genetic potential of microbial communities in permafrost would reflect microbial strategies to metabolize permafrost soluble organic matter (OM) in situ over geologic time. We also hypothesized that changes in the metagenome across the chronosequence would correlate with shifts in carbon chemistry, permafrost age, and paleoclimate at the time of permafrost formation. We combined high-resolution characterization of water-soluble OM by Fourier-transform ion-cyclotron-resonance mass spectrometry (FT-ICR MS), quantification of organic anions in permafrost water extracts, and metagenomic sequencing to better understand the relationships between the molecular-level composition of potentially bioavailable OM, the microbial community, and permafrost age. Both age and paleoclimate had marked effects on both the molecular composition of dissolved OM and the microbial community. The relative abundance of genes associated with hydrogenotrophic methanogenesis, carbohydrate active enzyme families, nominal oxidation state of carbon (NOSC), and number of identifiable molecular formulae significantly decreased with increasing age. In contrast, genes associated with fermentation of short chain fatty acids (SCFAs), the concentration of SCFAs and ammonium all significantly increased with age. We present a conceptual model of microbial metabolism in permafrost based on fermentation of OM and the buildup of organic acids that helps to explain the unique chemistry of ancient permafrost soils. These findings imply long-term in situ microbial turnover of ancient permafrost OM and that this pooled biolabile OM could prime ancient permafrost soils for a larger and more rapid microbial response to thaw compared to younger permafrost soils.

Molecular phylogenetic biodiversity assessment of arctic and boreal ectomycorrhizal Lactarius Pers. (Russulales; Basidiomycota) in Alaska, based on soil and sporocarp DNA.

Despite the critical roles fungi play in the functioning of ecosystems, especially as symbionts of plants and recyclers of organic matter, their biodiversity is poorly known in high-latitude regions. In this paper, we discuss the molecular diversity of one of the most diverse and abundant groups of ectomycorrhizal fungi: the genus Lactarius Pers. We analysed internal transcribed spacer rDNA sequences from both curated sporocarp collections and soil polymerase chain reaction clone libraries sampled in the arctic tundra and boreal forests of Alaska. Our genetic diversity assessment, based on various phylogenetic methods and operational taxonomic unit (OTU) delimitations, suggests that the genus Lactarius is diverse in Alaska, with at least 43 putative phylogroups, and 24 and 38 distinct OTUs based on 95% and 97% internal transcribed spacer sequence similarity, respectively. Some OTUs were identified to known species, while others were novel, previously unsequenced groups. Non-asymptotic species accumulation curves, the disparity between observed and estimated richness, and the high number of singleton OTUs indicated that many Lactarius species remain to be found and identified in Alaska. Many Lactarius taxa show strong habitat preference to one of the three major vegetation types in the sampled regions (arctic tundra, black spruce forests, and mixed birch-aspen-white spruce forests), as supported by statistical tests of UniFrac distances and principal coordinates analyses (PCoA). Together, our data robustly demonstrate great diversity and nonrandom ecological partitioning in an important boreal ectomycorrhizal genus within a relatively small geographical region. The observed diversity of Lactarius was much higher in either type of boreal forest than in the arctic tundra, supporting the widely recognized pattern of decreasing species richness with increasing latitude.

Large loss of CO2 in winter observed across the northern permafrost region.

Recent warming in the Arctic, which has been amplified during the winter1-3, greatly enhances microbial decomposition of soil organic matter and subsequent release of carbon dioxide (CO2)4. However, the amount of CO2 released in winter is highly uncertain and has not been well represented by ecosystem models or by empirically-based estimates5,6. Here we synthesize regional in situ observations of CO2 flux from arctic and boreal soils to assess current and future winter carbon losses from the northern permafrost domain. We estimate a contemporary loss of 1662 Tg C yr-1 from the permafrost region during the winter season (October through April). This loss is greater than the average growing season carbon uptake for this region estimated from process models (-1032 Tg C yr-1). Extending model predictions to warmer conditions in 2100 indicates that winter CO2 emissions will increase 17% under a moderate mitigation scenario-Representative Concentration Pathway (RCP) 4.5-and 41% under business-as-usual emissions scenario-RCP 8.5. Our results provide a new baseline for winter CO2 emissions from northern terrestrial regions and indicate that enhanced soil CO2 loss due to winter warming may offset growing season carbon uptake under future climatic conditions.

The evolutionary history of mycorrhizal specificity among lady's slipper orchids.

Although coevolution is acknowledged to occur in nature, coevolutionary patterns in symbioses not involving species-to-species relationships are poorly understood. Mycorrhizal plants are thought to be too generalist to coevolve with their symbiotic fungi; yet some plants, including some orchids, exhibit strikingly narrow mycorrhizal specificity. Here, we assess the evolutionary history of mycorrhizal specificity in the lady's slipper orchid genus, Cypripedium. We sampled 90 populations of 15 taxa across three continents, using DNA methods to identify fungal symbionts and quantify mycorrhizal specificity. We assessed phylogenetic relationships among sampled Cypripedium taxa, onto which we mapped mycorrhizal specificity. Cypripedium taxa associated almost exclusively with fungi within family Tulasnellaceae. Ancestral specificity appears to have been narrow, followed by a broadening after the divergence of C. debile. Specificity then narrowed, resulting in strikingly narrow specificity in most of the taxa in this study, with no taxon rewidening to the same extant as basal members of the genus. Sympatric taxa generally associated with different sets of fungi, and most clades of Cypripedium-mycorrhizal fungi were found throughout much of the northern hemisphere, suggesting that these evolutionary patterns in specificity are not the result of biogeographic lack of opportunity to associate with potential partners. Mycorrhizal specificity in genus Cypripedium appears to be an evolvable trait, and associations with particular fungi are phylogenetically conserved.

Patterns in wetland microbial community composition and functional gene repertoire associated with methane emissions.

UNLABELLED: Wetland restoration on peat islands previously drained for agriculture has potential to reverse land subsidence and sequester atmospheric carbon dioxide as peat accretes. However, the emission of methane could potentially offset the greenhouse gas benefits of captured carbon. As microbial communities play a key role in governing wetland greenhouse gas fluxes, we are interested in how microbial community composition and functions are associated with wetland hydrology, biogeochemistry, and methane emission, which is critical to modeling the microbial component in wetland methane fluxes and to managing restoration projects for maximal carbon sequestration. Here, we couple sequence-based methods with biogeochemical and greenhouse gas measurements to interrogate microbial communities from a pilot-scale restored wetland in the Sacramento-San Joaquin Delta of California, revealing considerable spatial heterogeneity even within this relatively small site. A number of microbial populations and functions showed strong correlations with electron acceptor availability and methane production; some also showed a preference for association with plant roots. Marker gene phylogenies revealed a diversity of major methane-producing and -consuming populations and suggested novel diversity within methanotrophs. Methanogenic archaea were observed in all samples, as were nitrate-, sulfate-, and metal-reducing bacteria, indicating that no single terminal electron acceptor was preferred despite differences in energetic favorability and suggesting spatial microheterogeneity and microniches. Notably, methanogens were negatively correlated with nitrate-, sulfate-, and metal-reducing bacteria and were most abundant at sampling sites with high peat accretion and low electron acceptor availability, where methane production was highest. IMPORTANCE: Wetlands are the largest nonanthropogenic source of atmospheric methane but also a key global carbon reservoir. Characterizing belowground microbial communities that mediate carbon cycling in wetlands is critical to accurately predicting their responses to changes in land management and climate. Here, we studied a restored wetland and revealed substantial spatial heterogeneity in biogeochemistry, methane production, and microbial communities, largely associated with the wetland hydraulic design. We observed patterns in microbial community composition and functions correlated with biogeochemistry and methane production, including diverse microorganisms involved in methane production and consumption. We found that methanogenesis gene abundance is inversely correlated with genes from pathways exploiting other electron acceptors, yet the ubiquitous presence of genes from all these pathways suggests that diverse electron acceptors contribute to the energetic balance of the ecosystem. These investigations represent an important step toward effective management of wetlands to reduce methane flux to the atmosphere and enhance belowground carbon storage.

Impact of fire on active layer and permafrost microbial communities and metagenomes in an upland Alaskan boreal forest.

Permafrost soils are large reservoirs of potentially labile carbon (C). Understanding the dynamics of C release from these soils requires us to account for the impact of wildfires, which are increasing in frequency as the climate changes. Boreal wildfires contribute to global emission of greenhouse gases (GHG-CO2, CH4 and N2O) and indirectly result in the thawing of near-surface permafrost. In this study, we aimed to define the impact of fire on soil microbial communities and metabolic potential for GHG fluxes in samples collected up to 1 m depth from an upland black spruce forest near Nome Creek, Alaska. We measured geochemistry, GHG fluxes, potential soil enzyme activities and microbial community structure via 16SrRNA gene and metagenome sequencing. We found that soil moisture, C content and the potential for respiration were reduced by fire, as were microbial community diversity and metabolic potential. There were shifts in dominance of several microbial community members, including a higher abundance of candidate phylum AD3 after fire. The metagenome data showed that fire had a pervasive impact on genes involved in carbohydrate metabolism, methanogenesis and the nitrogen cycle. Although fire resulted in an immediate release of CO2 from surface soils, our results suggest that the potential for emission of GHG was ultimately reduced at all soil depths over the longer term. Because of the size of the permafrost C reservoir, these results are crucial for understanding whether fire produces a positive or negative feedback loop contributing to the global C cycle.

Increasing ecological inference from high throughput sequencing of fungi in the environment through a tagging approach.

High throughput sequencing methods are widely used in analyses of microbial diversity, but are generally applied to small numbers of samples, which precludes characterization of patterns of microbial diversity across space and time. We have designed a primer-tagging approach that allows pooling and subsequent sorting of numerous samples, which is directed to amplification of a region spanning the nuclear ribosomal internal transcribed spacers and partial large subunit from fungi in environmental samples. To test the method for phylogenetic biases, we constructed a controlled mixture of four taxa representing the Chytridiomycota, Zygomycota, Ascomycota and Basidiomycota. Following cloning and colony restriction fragment length polymorphism analysis, we found no significant difference in representation in 19 of the 23 tested primers. We also generated a clone library from two soil DNA extracts using two primers for each extract and compared 456 clone sequences. Community diversity statistics and contingency table tests applied to counts of operational taxonomic units revealed that the two DNA extracts differed significantly, while the pairs of tagged primers from each extract were indistinguishable. Similar results were obtained using UniFrac phylogenetic comparisons. Together, these results suggest that the pig-tagged primers can be used to increase ecological inference in high throughput sequencing projects on fungi.