Steven J. Blazewicz

Metagenomic analysis of a permafrost microbial community reveals a rapid response to thaw

Permafrost contains an estimated 1672 Pg carbon (C), an amount roughly equivalent to the total currently contained within land plants and the atmosphere. This reservoir of C is vulnerable to decomposition as rising global temperatures cause the permafrost to thaw. During thaw, trapped organic matter may become more accessible for microbial degradation and result in greenhouse gas emissions. Despite recent advances in the use of molecular tools to study permafrost microbial communities, their response to thaw remains unclear. Here we use deep metagenomic sequencing to determine the impact of thaw on microbial phylogenetic and functional genes, and relate these data to measurements of methane emissions. Metagenomics, the direct sequencing of DNA from the environment, allows the examination of whole biochemical pathways and associated processes, as opposed to individual pieces of the metabolic puzzle. Our metagenome analyses reveal that during transition from a frozen to a thawed state there are rapid shifts in many microbial, phylogenetic and functional gene abundances and pathways. After one week of incubation at 5 °C, permafrost metagenomes converge to be more similar to each other than while they are frozen. We find that multiple genes involved in cycling of C and nitrogen shift rapidly during thaw. We also construct the first draft genome from a complex soil metagenome, which corresponds to a novel methanogen. Methane previously accumulated in permafrost is released during thaw and subsequently consumed by methanotrophic bacteria. Together these data point towards the importance of rapid cycling of methane and nitrogen in thawing permafrost.

Evaluating rRNA as an indicator of microbial activity in environmental communities: limitations and uses

Microbes exist in a range of metabolic states (for example, dormant, active and growing) and analysis of ribosomal RNA (rRNA) is frequently employed to identify the ‘active’ fraction of microbes in environmental samples. While rRNA analyses are no longer commonly used to quantify a population’s growth rate in mixed communities, due to rRNA concentration not scaling linearly with growth rate uniformly across taxa, rRNA analyses are still frequently used toward the more conservative goal of identifying populations that are currently active in a mixed community. Yet, evidence indicates that the general use of rRNA as a reliable indicator of metabolic state in microbial assemblages has serious limitations. This report highlights the complex and often contradictory relationships between rRNA, growth and activity. Potential mechanisms for confounding rRNA patterns are discussed, including differences in life histories, life strategies and non-growth activities. Ways in which rRNA data can be used for useful characterization of microbial assemblages are presented, along with questions to be addressed in future studies.

Abundance of microbial genes associated with nitrogen cycling as indices of biogeochemical process rates across a vegetation gradient in Alaska

Nitrification and denitrification processes are crucial to plant nutrient availability, eutrophication and greenhouse gas production both locally and globally. Unravelling the major environmental predictors for nitrification and denitrification is thus pivotal in order to understand and model environmental nitrogen (N) cycling. Here, we sampled five plant community types characteristic of interior Alaska, including black spruce, bog birch, tussock grass and two fens. We assessed abundance of functional genes affiliated with nitrification (bacterial and archaeal amoA) and denitrification (nirK/S and nosZ) using qPCR, soil characteristics, potential nitrification and denitrification rates (PNR and PDR) and gross mineralization rates. The main chemical and biological predictors for PNR and PDR were assigned through path analysis. The potential N cycling rates varied dramatically between sites, from some of the highest (in fens) to some of the lowest (in black spruce) measured globally. Based on path analysis, functional gene abundances were the most important variables to predict potential rates. PNR was best explained by bacterial amoA gene abundance followed by ammonium content, whereas PDR was best explained directly by nosZ gene abundance and indirectly by nirK/S gene abundance and nitrate. Hence, functional gene abundance is a valuable index that integrates recent environmental history and recent process activity, and therefore is a good predictor of potential rates. The results of this study contribute to our understanding of the relative importance of different biological and chemical factors in driving the potential for nitrification and denitrification across terrestrial ecosystems.

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.

Growth and death of bacteria and fungi underlie rainfall-induced carbon dioxide pulses from seasonally dried soil

The rapid increase in microbial activity that occurs when a dry soil is rewetted has been well documented and is of great interest due to implications of changing precipitation patterns on soil C dynamics. Several studies have shown minor net changes in microbial population diversity or abundance following wet-up, but the gross population dynamics of bacteria and fungi resulting from soil wet-up are virtually unknown. Here we applied DNA stable isotope probing with H218O coupled with quantitative PCR to characterize new growth, survival, and mortality of bacteria and fungi following the rewetting of a seasonally dried California annual grassland soil. Microbial activity, as determined by CO2 production, increased significantly within three hours of wet-up, yet new growth was not detected until after three hours, suggesting a pulse of nongrowth activity immediately following wet-up, likely due to osmo-regulation and resuscitation from dormancy in response to the rapid change in water potential. Total microbial abundance revealed little change throughout the seven-day post-wet incubation, but there was substantial turnover of both bacterial and fungal populations (49% and 52%, respectively). New growth was linear between 24 and 168 hours for both bacteria and fungi, with average growth rates of 2.3 x 10(8) bacterial 16S rRNA gene copies x [g dry mass](-1) x h(-1) and 4.3 x 10(7) fungal ITS copies x [g dry mass](-1) x h(-1). While bacteria and fungi differed in their mortality and survival characteristics during the seven-day incubation, mortality that occurred within the first three hours was similar, with 25% and 27% of bacterial and fungal gene copies disappearing from the pre-wet community, respectively. The rapid disappearance of gene copies indicates that cell death, occurring either during the extreme dry down period (preceding five months) or during the rapid change in water potential due to wet-up, generates a significant pool of available C that likely contributes to the large pulse in CO2 associated with wet-up. A dynamic assemblage of growing and dying organisms controlled the CO2 pulse, but the balance between death and growth resulted in relatively stable total population abundances, even after a profound and sudden change in environment.

Modeling CH4 and CO2 cycling using porewater stable isotopes in a thermokarst bog in Interior Alaska: results from three conceptual reaction networks

Quantifying rates of microbial carbon transformation in peatlands is essential for gaining mechanistic understanding of the factors that influence methane emissions from these systems, and for predicting how emissions will respond to climate change and other disturbances. In this study, we used porewater stable isotopes collected from both the edge and center of a thermokarst bog in Interior Alaska to estimate in situ microbial reaction rates. We expected that near the edge of the thaw feature, actively thawing permafrost and greater abundance of sedges would increase carbon, oxygen and nutrient availability, enabling faster microbial rates relative to the center of the thaw feature. We developed three different conceptual reaction networks that explained the temporal change in porewater CO2, CH4, δ13C–CO2 and δ13C–CH4. All three reaction-network models included methane production, methane oxidation and CO2 production, and two of the models included homoacetogenesis—a reaction not previously included in isotope-based porewater models. All three models fit the data equally well, but rates resulting from the models differed. Most notably, inclusion of homoacetogenesis altered the modeled pathways of methane production when the reaction was directly coupled to methanogenesis, and it decreased gross methane production rates by up to a factor of five when it remained decoupled from methanogenesis. The ability of all three conceptual reaction networks to successfully match the measured data indicate that this technique for estimating in situ reaction rates requires other data and information from the site to confirm the considered set of microbial reactions. Despite these differences, all models indicated that, as expected, rates were greater at the edge than in the center of the thaw bog, that rates at the edge increased more during the growing season than did rates in the center, and that the ratio of acetoclastic to hydrogenotrophic methanogenesis was greater at the edge than in the center. In both locations, modeled rates (excluding methane oxidation) increased with depth. A puzzling outcome from the effort was that none of the models could fit the porewater dataset without generating “fugitive” carbon (i.e., methane or acetate generated by the models but not detected at the field site), indicating that either our conceptualization of the reactions occurring at the site remains incomplete or our site measurements are missing important carbon transformations and/or carbon fluxes. This model–data discrepancy will motivate and inform future research efforts focused on improving our understanding of carbon cycling in permafrost wetlands.

Stable isotope informed genome-resolved metagenomics reveals that Saccharibacteria utilize microbially-processed plant-derived carbon

BackgroundThe transformation of plant photosynthate into soil organic carbon and its recycling to CO2 by soil microorganisms is one of the central components of the terrestrial carbon cycle. There are currently large knowledge gaps related to which soil-associated microorganisms take up plant carbon in the rhizosphere and the fate of that carbon.ResultsWe conducted an experiment in which common wild oats (Avena fatua) were grown in a 13CO2 atmosphere and the rhizosphere and non-rhizosphere soil was sampled for genomic analyses. Density gradient centrifugation of DNA extracted from soil samples enabled distinction of microbes that did and did not incorporate the 13C into their DNA. A 1.45-Mbp genome of a Saccharibacteria (TM7) was identified and, despite the microbial complexity of rhizosphere soil, curated to completion. The genome lacks many biosynthetic pathways, including genes required to synthesize DNA de novo. Rather, it requires externally derived nucleotides for DNA and RNA synthesis. Given this, we conclude that rhizosphere-associated Saccharibacteria recycle DNA from bacteria that live off plant exudates and/or phage that acquired 13C because they preyed upon these bacteria and/or directly from the labeled plant DNA. Isotopic labeling indicates that the population was replicating during the 6-week period of plant growth. Interestingly, the genome is ~ 30% larger than other complete Saccharibacteria genomes from non-soil environments, largely due to more genes for complex carbon utilization and amino acid metabolism. Given the ability to degrade cellulose, hemicellulose, pectin, starch, and 1,3-β-glucan, we predict that this Saccharibacteria generates energy by fermentation of soil necromass and plant root exudates to acetate and lactate. The genome also encodes a linear electron transport chain featuring a terminal oxidase, suggesting that this Saccharibacteria may respire aerobically. The genome encodes a hydrolase that could breakdown salicylic acid, a plant defense signaling molecule, and genes to interconvert a variety of isoprenoids, including the plant hormone zeatin.ConclusionsRhizosphere Saccharibacteria likely depend on other bacteria for basic cellular building blocks. We propose that isotopically labeled CO2 is incorporated into plant-derived carbon and then into the DNA of rhizosphere organisms capable of nucleotide synthesis, and the nucleotides are recycled into Saccharibacterial genomes.

Inter-laboratory testing of the effect of DNA blocking reagent G2 on DNA extraction from low-biomass clay samples

Here we show that a commercial blocking reagent (G2) based on modified eukaryotic DNA significantly improved DNA extraction efficiency. We subjected G2 to an inter-laboratory testing, where DNA was extracted from the same clay subsoil using the same batch of kits. The inter-laboratory extraction campaign revealed large variation among the participating laboratories, but the reagent increased the number of PCR-amplified16S rRNA genes recovered from biomass naturally present in the soils by one log unit. An extensive sequencing approach demonstrated that the blocking reagent was free of contaminating DNA, and may therefore also be used in metagenomics studies that require direct sequencing.

Multiple element isotope probes, NanoSIMS, and the functional genomics of microbial carbon cycling in soils in response to chronic climatic change

Center for Ecosystem Science and Society, Northern Arizona University, Flagstaff, AZ, 86011, USA; Division of Plant and Soil Sciences, West Virginia University, Morgantown, WV, 26505; Department of Biological Sciences, Northern Arizona University, Flagstaff, AZ, 86011, USA; Center for Microbial Genetics and Genomics, Northern Arizona University, Flagstaff, AZ, 86011, USA; Center for Microbiomics and Human Health, Translational Genomics Research Institute, Flagstaff AZ, 86001, USA; Department of Environmental and Occupational Health, Milken Institute School of Public Health, George Washington University, Washington DC, 20037, USA; 7 Lawrence Livermore National Laboratory, Livermore, CA, 94550, USA