Intra-genomic variation in G + C content and its implications for DNA stable isotope probing

The effect of the 13C abundance of soil microbial DNA on identifying labelled fractions after ultracentrifugation

DNA-based stable isotope probing (DNA-SIP) technology has been widely employed to trace microbes assimilating target substrates. However, the fractions with labelled universal genes are sometimes difficult to distinguish when detected by quantitative real-time PCR. In this experiment, three paddy soils (AQ, CZ, and NB) were amended with 0.1% glucose containing 13C at six levels, and DNA was then extracted after a 7-day incubation and subjected to isopycnic gradient centrifugation. The results showed that the amount of labelled DNA was notably related to the 13C-glucose percentage, while the separation spans of 18S rRNA and 16S rRNA genes between labelled and unlabelled treatments became notably clearer when the δ13C values of the total DNA were 90.9, 61.6, and 38.9‰ and 256.2, 104.5 and 126.1‰ in the AQ, CZ, and NB soils, respectively. Moreover, fractionated DNA was also labelled by determining the δ13C values while adding only 5 atom% 13C-glucose to the soil. The results suggest that the optimal labelling fractions were not always those fractions with the maximal gene abundance, and detecting the δ13C values of the total and fractionated DNA was beneficial in estimating the results of DNA-SIP. KEY POINTS: • Appropriate 13C-DNA amount was needed for DNA-SIP. • Detecting the 13C ratio of fractionated DNA directly was an assistant method for identifying the labelled fractions. • Fractions with the maximal 18S or 16S rRNA gene abundance always were not labelled.

MISIP: a data standard for the reuse and reproducibility of any stable isotope probing-derived nucleic acid sequence and experiment

DNA/RNA-stable isotope probing (SIP) is a powerful tool to link in situ microbial activity to sequencing data. Every SIP dataset captures distinct information about microbial community metabolism, process rates, and population dynamics, offering valuable insights for a wide range of research questions. Data reuse maximizes the information derived from the labor and resource-intensive SIP approaches. Yet, a review of publicly available SIP sequencing metadata showed that critical information necessary for reproducibility and reuse was often missing. Here, we outline the Minimum Information for any Stable Isotope Probing Sequence (MISIP) according to the Minimum Information for any (x) Sequence (MIxS) framework and include examples of MISIP reporting for common SIP experiments. Our objectives are to expand the capacity of MIxS to accommodate SIP-specific metadata and guide SIP users in metadata collection when planning and reporting an experiment. The MISIP standard requires 5 metadata fields-isotope, isotopolog, isotopolog label, labeling approach, and gradient position-and recommends several fields that represent best practices in acquiring and reporting SIP sequencing data (e.g., gradient density and nucleic acid amount). The standard is intended to be used in concert with other MIxS checklists to comprehensively describe the origin of sequence data, such as for marker genes (MISIP-MIMARKS) or metagenomes (MISIP-MIMS), in combination with metadata required by an environmental extension (e.g., soil). The adoption of the proposed data standard will improve the reuse of any sequence derived from a SIP experiment and, by extension, deepen understanding of in situ biogeochemical processes and microbial ecology.

Multisubstrate DNA stable isotope probing reveals guild structure of bacteria that mediate soil carbon cycling

Soil microorganisms determine the fate of soil organic matter (SOM), and their activities compose a major component of the global carbon (C) cycle. We employed a multisubstrate, DNA-stable isotope probing experiment to track bacterial assimilation of C derived from distinct sources that varied in bioavailability. This approach allowed us to measure microbial contributions to SOM processing by measuring the C assimilation dynamics of diverse microorganisms as they interacted within soil. We identified and tracked 1,286 bacterial taxa that assimilated 13C in an agricultural soil over a period of 48 d. Overall 13C-assimilation dynamics of bacterial taxa, defined by the source and timing of the 13C they assimilated, exhibited low phylogenetic conservation. We identified bacterial guilds composed of taxa that had similar 13C assimilation dynamics. We show that C-source bioavailability explained significant variation in both C mineralization dynamics and guild structure, and that the growth dynamics of bacterial guilds differed significantly in response to C addition. We also demonstrate that the guild structure explains significant variation in the biogeographical distribution of bacteria at continental and global scales. These results suggest that an understanding of in situ growth dynamics is essential for understanding microbial contributions to soil C cycling. We interpret these findings in the context of bacterial life history strategies and their relationship to terrestrial C cycling.

Metagenomics and Quantitative Stable Isotope Probing Offer Insights into Metabolism of Polycyclic Aromatic Hydrocarbon Degraders in Chronically Polluted Seawater

Bacterial biodegradation is a significant contributor to remineralization of polycyclic aromatic hydrocarbons (PAHs)-toxic and recalcitrant components of crude oil as well as by-products of partial combustion chronically introduced into seawater via atmospheric deposition. The Deepwater Horizon oil spill demonstrated the speed at which a seed PAH-degrading community maintained by chronic inputs responds to acute pollution. We investigated the diversity and functional potential of a similar seed community in the chronically polluted Port of Los Angeles (POLA), using stable isotope probing with naphthalene, deep-sequenced metagenomes, and carbon incorporation rate measurements at the port and in two sites in the San Pedro Channel. We demonstrate the ability of the community of degraders at the POLA to incorporate carbon from naphthalene, leading to a quick shift in microbial community composition to be dominated by the normally rare Colwellia and Cycloclasticus We show that metagenome-assembled genomes (MAGs) belonged to these naphthalene degraders by matching their 16S-rRNA gene with experimental stable isotope probing data. Surprisingly, we did not find a full PAH degradation pathway in those genomes, even when combining genes from the entire microbial community, leading us to hypothesize that promiscuous dehydrogenases replace canonical naphthalene degradation enzymes in this site. We compared metabolic pathways identified in 29 genomes whose abundance increased in the presence of naphthalene to generate genomic-based recommendations for future optimization of PAH bioremediation at the POLA, e.g., ammonium as opposed to urea, heme or hemoproteins as an iron source, and polar amino acids.IMPORTANCE Oil spills in the marine environment have a devastating effect on marine life and biogeochemical cycles through bioaccumulation of toxic hydrocarbons and oxygen depletion by hydrocarbon-degrading bacteria. Oil-degrading bacteria occur naturally in the ocean, especially where they are supported by chronic inputs of oil or other organic carbon sources, and have a significant role in degradation of oil spills. Polycyclic aromatic hydrocarbons are the most persistent and toxic component of crude oil. Therefore, the bacteria that can break those molecules down are of particular importance. We identified such bacteria at the Port of Los Angeles (POLA), one of the busiest ports worldwide, and characterized their metabolic capabilities. We propose chemical targets based on those analyses to stimulate the activity of these bacteria in case of an oil spill in the Port POLA.

Competitive Exclusion and Metabolic Dependency among Microorganisms Structure the Cellulose Economy of an Agricultural Soil

Microorganisms that degrade cellulose utilize extracellular reactions that yield free by-products which can promote interactions with noncellulolytic organisms. We hypothesized that these interactions determine the ecological and physiological traits governing the fate of cellulosic carbon (C) in soil. We performed comparative genomics with genome bins from a shotgun metagenomic-stable isotope probing experiment to characterize the attributes of cellulolytic and noncellulolytic taxa accessing ¹³C from cellulose. We hypothesized that cellulolytic taxa would exhibit competitive traits that limit access, while noncellulolytic taxa would display greater metabolic dependency, such as signatures of adaptive gene loss. We tested our hypotheses by evaluating genomic traits indicative of competitive exclusion or metabolic dependency, such as antibiotic production, growth rate, surface attachment, biomass degrading potential, and auxotrophy. The most ¹³C-enriched taxa were cellulolytic Cellvibrio (Gammaproteobacteria) and Chaetomium (Ascomycota), which exhibited a strategy of self-sufficiency (prototrophy), rapid growth, and competitive exclusion via antibiotic production. Auxotrophy was more prevalent in cellulolytic Actinobacteria than in cellulolytic Proteobacteria, demonstrating differences in dependency among cellulose degraders. Noncellulolytic taxa that accessed ¹³C from cellulose (Planctomycetales, Verrucomicrobia, and Vampirovibrionales) were also more dependent, as indicated by patterns of auxotrophy and ¹³C labeling (i.e., partial labeling or labeling at later stages). Major ¹³C-labeled cellulolytic microbes (e.g., Sorangium, Actinomycetales, Rhizobiales, and Caulobacteraceae) possessed adaptations for surface colonization (e.g., gliding motility, hyphae, attachment structures) signifying the importance of surface ecology in decomposing particulate organic matter. Our results demonstrated that access to cellulosic C was accompanied by ecological trade-offs characterized by differing degrees of metabolic dependency and competitive exclusion.

Influence of reductive soil disinfestation or biochar amendment on bacterial communities and their utilization of plant-derived carbon in the rhizosphere of tomato

Root-associated microorganisms play an important role in plant nutrition and productivity. However, our understanding of how a plant-microbiome system responds to pre-planting soil management remains limited. Here, continuous labeling with ¹³CO₂ gas combined with stable isotope probing (SIP) was applied to explore bacterial utilization of plant-derived carbon (C) in the tomato rhizosphere as affected by biochar amendment or reductive soil disinfestation (RSD). Our results showed that RSD treatment strongly shaped the soil bacterial community composition, while biochar soil amendment had little impact on the community in the rhizosphere of tomato. We observed that the bacterial community in the RSD treatment, which actively utilized plant-derived C, belonged to various phyla (i.e., Proteobacteria, Cyanobacteria, Verrucomicrobia, and Acidobacteria), while the genus Streptomyces (phylum Actinobacteria) was the main bacterial taxa that actively utilized plant-derived C in the biochar and control treatments. This study provides evidence that biochar application or RSD pre-planting soil management practices induced distinct bacterial utilization of plant-derived C, which may in turn regulate plant productivity in agricultural systems. KEY POINTS: • Genus Streptomyces was the main bacterial group utilizing plant-derived carbon in both control and biochar treatments. • Reductive soil disinfestation altered bacterial utilization of plant-derived carbon. • Biochar did not alter the composition of the bacterial communities but had more labeled bacterial taxa utilizing plant-derived carbon.

DNA-Stable Isotope Probing Shotgun Metagenomics Reveals the Resilience of Active Microbial Communities to Biochar Amendment in Oxisol Soil

The functions and interactions of individual microbial populations and their genes in agricultural soils amended with biochar remain elusive but are crucial for a deeper understanding of nutrient cycling and carbon (C) sequestration. In this study, we coupled DNA stable isotope probing (SIP) with shotgun metagenomics in order to target the active community in microcosms which contained soil collected from biochar-amended and control plots under napiergrass cultivation. Our analyses revealed that the active community was composed of high-abundant and low-abundant populations, including Actinobacteria, Proteobacteria, Gemmatimonadetes, and Acidobacteria. Although biochar did not significantly shift the active taxonomic and functional communities, we found that the narG (nitrate reductase) gene was significantly more abundant in the control metagenomes. Interestingly, putative denitrifier genomes generally encoded one gene or a partial denitrification pathway, suggesting denitrification is typically carried out by an assembly of different populations within this Oxisol soil. Altogether, these findings indicate that the impact of biochar on the active soil microbial community are transient in nature. As such, the addition of biochar to soils appears to be a promising strategy for the long-term C sequestration in agricultural soils, does not impart lasting effects on the microbial functional community, and thus mitigates un-intended microbial community shifts that may lead to fertilizer loss through increased N cycling.

Optimizing ultracentrifugation conditions for DNA-based stable isotope probing (DNA-SIP)

DNA-SIP (DNA-based stable isotope probing) is increasingly being employed in soil microbial ecology to identify those microbes assimilating the ¹³C/¹⁵N labelled substrate. Isopycnic gradient centrifugation is the primary experimental process for conducting DNA-SIP. However, diverse centrifugal conditions have been used in various recent studies. In order to get the optimum conditions of centrifugation for DNA-SIP, centrifugation time (36, 42, 48, 60 h), speed (45,000, 55,000 rpm) and the initial buoyant density (1.69, 1.71, 1.725 g ml^-1), as were used extensively in related studies, were tested in this experiment with the Vti 65.2 rotor. DNA with either ¹³C-labelling or unlabelled was extracted from a paddy soil pre-incubated with either ¹³C-labelled or natural abundance glucose. After ultracentrifugation, the gene abundance of bacterial 16S rRNA, fungal 18S rRNA, bacterial and archaeal amoA within the fractioned DNA was detected. The results showed that centrifugation for 48 h was enough for the DNA to reach stabilization in the CsCl solution. The initial density of the mixed solution was best adjusted to 1.71 g ml^-1 to ensure that most of the genes were concentrated on the middle fractions of the density gradient. Increasing the centrifugation speed would increase the density gradient of fractions; therefore, 45,000 rpm (184,000 g) was recommended so as to obtain the more widespread pattern of DNA in the centrifugal tube. We hope these findings will assist future researchers to conduct optimum ultracentrifugation for DNA-SIP.

Biochar Amendment Stimulates Utilization of Plant-Derived Carbon by Soil Bacteria in an Intercropping System

Plant-derived carbon (C) is considered fundamental to understand the interaction between rhizosphere microbes and plants in terrestrial ecosystems. Biochar soil amendment may enhance plant performance via changing soil properties or microbial diversity in the rhizosphere. However, our knowledge of how plant-microbiome associations respond to biochar amendment remains rather limited. Herein, ¹³CO₂ steady-state labeling combined with DNA stable-isotope probing was used to characterize soil bacterial communities in the rhizosphere contributing to the utilization of plant-derived C. The diversity of bacteria active in the utilization of root exudates was determined after biochar amendment in a legume-based intercropping system (Vicia faba L., with Zea mays L.). The results showed the biochar application not only changed the bacterial community structure and diversity in the rhizosphere, but also altered bacterial members actively assimilating plant-derived C. There were more labeled species in the biochar-amended soils than the control soils. Compared with the control, the biochar amendment increased the relative abundances of Firmicutes and Bacteroidetes members (i.e., Bacillus, Clostridium, Sporomusa, Desulfosporosinus, and Alicyclobacillus) while decreasing the abundances of Proteobacteria members (e.g., Methylobacterium and Sphingomonas) utilizing plant-derived C. In contrast, slow-growing species of the phyla Acidobacteria, Planctomycetes, and Gemmatimonadetes were barely labeled. The bacteria found stimulated by the biochar amendment are known for their ability to fix nitrogen, solubilize phosphorus, or reduce iron and sulfur, which may potentially contribute to the “biochar effect” in the rhizosphere. This study is the first to provide empirical evidence that biochar amendment can alter the soil bacterial community assimilating plant-derived C; this may have consequences for nutrient cycling and improving plant performance in intercropping systems.

Experimental Setup and Data Analysis Considerations for DNA- and RNA-SIP Experiments in the Omics Era

Careful and thoughtful experimental design is crucial to the success of any SIP experiment. This chapter discusses the essential aspects of designing a SIP experiment, focusing primarily on DNA- and RNA-SIP. The design aspects discussed here begin with considerations for carrying out the incubation, such as, the effect of choosing different stable isotopes and target biomolecules, to what degree should a labeled substrate be enriched, what concentration to use, and how long should the incubation take. Then tips and pitfalls in the technical execution of SIP are listed, including how much nucleic acids should be loaded, how many fractions to collect, and what centrifuge rotor to use. Lastly, a brief overview of the current methods for analyzing SIP data is presented, focusing on high-throughput amplicon sequencing, together with a discussion on how the choice of analysis method might affect the experimental design.

Data Analysis for DNA Stable Isotope Probing Experiments Using Multiple Window High-Resolution SIP

DNA stable isotope probing (DNA-SIP) allows for the identification of microbes that assimilate isotopically labeled substrates into DNA. Here we describe the analysis of sequencing data using the multiple window high-resolution DNA-SIP method (MW-HR-SIP). MW-HR-SIP has improved accuracy over other methods and is easily implemented on the statistical platform R. We also discuss key experimental parameters to consider when designing DNA-SIP experiments and how these parameters affect accuracy of analysis.

Bacterial contributions to delignification and lignocellulose degradation in forest soils with metagenomic and quantitative stable isotope probing

Delignification, or lignin-modification, facilitates the decomposition of lignocellulose in woody plant biomass. The extant diversity of lignin-degrading bacteria and fungi is underestimated by culture-dependent methods, limiting our understanding of the functional and ecological traits of decomposers populations. Here, we describe the use of stable isotope probing (SIP) coupled with amplicon and shotgun metagenomics to identify and characterize the functional attributes of lignin, cellulose and hemicellulose-degrading fungi and bacteria in coniferous forest soils from across North America. We tested the extent to which catabolic genes partitioned among different decomposer taxa; the relative roles of bacteria and fungi, and whether taxa or catabolic genes correlated with variation in lignocellulolytic activity, measured as the total assimilation of ¹³C-label into DNA and phospholipid fatty acids. We found high overall bacterial degradation of our model lignin substrate, particularly by gram-negative bacteria (Comamonadaceae and Caulobacteraceae), while fungi were more prominent in cellulose-degradation. Very few taxa incorporated ¹³C-label from more than one lignocellulosic polymer, suggesting specialization among decomposers. Collectively, members of Caulobacteraceae could degrade all three lignocellulosic polymers, providing new evidence for their importance in lignocellulose degradation. Variation in lignin-degrading activity was better explained by microbial community properties, such as catabolic gene content and community structure, than cellulose-degrading activity. SIP significantly improved shotgun metagenome assembly resulting in the recovery of several high-quality draft metagenome-assembled genomes and over 7500 contigs containing unique clusters of carbohydrate-active genes. These results improve understanding of which organisms, conditions and corresponding functional genes contribute to lignocellulose decomposition.

Community dynamics and functional characteristics of naphthalene-degrading populations in contaminated surface sediments and hypoxic/anoxic groundwater

Earlier research on the biogeochemical factors affecting natural attenuation in coal-tar contaminated groundwater, at South Glens Falls, NY, revealed the importance of anaerobic metabolism and trophic interactions between degrader and bacterivore populations. Field-based characterizations of both phenomena have proven challenging, but advances in stable isotope probing (SIP), single-cell imaging and shotgun metagenomics now provide cultivation-independent tools for their study. We tracked carbon from ¹³ C-labelled naphthalene through microbial populations in contaminated surface sediments over 6 days using respiration assays, secondary ion mass spectrometry imaging and shotgun metagenomics to disentangle the contaminant-based trophic web. Contaminant-exposed communities in hypoxic/anoxic groundwater were contrasted with those from oxic surface sediments to identify putative features of anaerobic catabolism of naphthalene. In total, six bacteria were responsible for naphthalene degradation. Cupriavidus, Ralstonia and Sphingomonas predominated at the earliest stages of SIP incubations and were succeeded in later stages by Stenotrophomonas and Rhodococcus. Metagenome-assembled genomes provided evidence for the ecological and functional characteristics underlying these temporal shifts. Identical species of Stenotrophomonas and Rhodococcus were abundant in the most contaminated, anoxic groundwater. Apparent increases in bacterivorous protozoa were observed following exposure to naphthalene, though insignificant amounts of carbon were transferred between bacterial degraders and populations of secondary feeders.

SIPSim: A Modeling Toolkit to Predict Accuracy and Aid Design of DNA-SIP Experiments

DNA Stable isotope probing (DNA-SIP) is a powerful method that links identity to function within microbial communities. The combination of DNA-SIP with multiplexed high throughput DNA sequencing enables simultaneous mapping of in situ assimilation dynamics for thousands of microbial taxonomic units. Hence, high throughput sequencing enabled SIP has enormous potential to reveal patterns of carbon and nitrogen exchange within microbial food webs. There are several different methods for analyzing DNA-SIP data and despite the power of SIP experiments, it remains difficult to comprehensively evaluate method accuracy across a wide range of experimental parameters. We have developed a toolset (SIPSim) that simulates DNA-SIP data, and we use this toolset to systematically evaluate different methods for analyzing DNA-SIP data. Specifically, we employ SIPSim to evaluate the effects that key experimental parameters (e.g., level of isotopic enrichment, number of labeled taxa, relative abundance of labeled taxa, community richness, community evenness, and beta-diversity) have on the specificity, sensitivity, and balanced accuracy (defined as the product of specificity and sensitivity) of DNA-SIP analyses. Furthermore, SIPSim can predict analytical accuracy and power as a function of experimental design and community characteristics, and thus should be of great use in the design and interpretation of DNA-SIP experiments.

HTSSIP: An R package for analysis of high throughput sequencing data from nucleic acid stable isotope probing (SIP) experiments

Combining high throughput sequencing with stable isotope probing (HTS-SIP) is a powerful method for mapping in situ metabolic processes to thousands of microbial taxa. However, accurately mapping metabolic processes to taxa is complex and challenging. Multiple HTS-SIP data analysis methods have been developed, including high-resolution stable isotope probing (HR-SIP), multi-window high-resolution stable isotope probing (MW-HR-SIP), quantitative stable isotope probing (qSIP), and ΔBD. Currently, there is no publicly available software designed specifically for analyzing HTS-SIP data. To address this shortfall, we have developed the HTSSIP R package, an open-source, cross-platform toolset for conducting HTS-SIP analyses in a straightforward and easily reproducible manner. The HTSSIP package, along with full documentation and examples, is available from CRAN at https://cran.r-project.org/web/packages/HTSSIP/index.html and Github at https://github.com/buckleylab/HTSSIP.

Inhibition of anaerobic digestion processes: Applications of molecular tools

Inhibition of anaerobic digestion (AD) due to perturbation caused by substrate composition and/or operating conditions can significantly reduce performance. Such perturbations could be limited by elucidating microbial community response to inhibitors and devising strategies to increase community resilience. To this end, advanced molecular methods are increasingly being applied to study the AD microbiome, a diverse community of microbial populations with complex interactions. This literature review of AD inhibition studies indicates that inhibitory concentrations are highly variable, likely stemming from differences in community structure or activity profile and previous exposure to inhibitors. More recent molecular methods such as 'omics' tools, substrate mapping, and real-time sequencing are helping to unravel the complexity of AD inhibition by elucidating physiological and ecological significance of key microbial populations. The AD community must strive towards developing predictive abilities to avoid system failure (e.g., real-time tracking of an indicator species) to improve resilience of AD systems.

A metagenomic survey of forest soil microbial communities more than a decade after timber harvesting

The scarcity of long-term data on soil microbial communities in the decades following timber harvesting limits current understanding of the ecological problems associated with maintaining the productivity of managed forests. The high complexity of soil communities and the heterogeneity of forest and soil necessitates a comprehensive approach to understand the role of microbial processes in managed forest ecosystems. Here, we describe a curated collection of well replicated, multi-faceted data from eighteen reforested sites in six different North American ecozones within the Long-term Soil Productivity (LTSP) Study, without detailed analysis of results or discussion. The experiments were designed to contrast microbial community composition and function among forest soils from harvested treatment plots with varying intensities of organic matter removal. The collection includes 724 bacterial (16S) and 658 fungal (ITS2) amplicon libraries, 133 shotgun metagenomic libraries as well as stable isotope probing amplicon libraries capturing the effects of harvesting on hemicellulolytic and cellulolytic populations. This collection serves as a foundation for the LTSP Study and other studies of the ecology of forest soil and forest disturbance.

Non-cyanobacterial diazotrophs mediate dinitrogen fixation in biological soil crusts during early crust formation

Biological soil crusts (BSCs) are key components of ecosystem productivity in arid lands and they cover a substantial fraction of the terrestrial surface. In particular, BSC N2-fixation contributes significantly to the nitrogen (N) budget of arid land ecosystems. In mature crusts, N2-fixation is largely attributed to heterocystous cyanobacteria; however, early successional crusts possess few N2-fixing cyanobacteria and this suggests that microorganisms other than cyanobacteria mediate N2-fixation during the critical early stages of BSC development. DNA stable isotope probing with (15)N2 revealed that Clostridiaceae and Proteobacteria are the most common microorganisms that assimilate (15)N2 in early successional crusts. The Clostridiaceae identified are divergent from previously characterized isolates, though N2-fixation has previously been observed in this family. The Proteobacteria identified share >98.5% small subunit rRNA gene sequence identity with isolates from genera known to possess diazotrophs (for example, Pseudomonas, Klebsiella, Shigella and Ideonella). The low abundance of these heterotrophic diazotrophs in BSCs may explain why they have not been characterized previously. Diazotrophs have a critical role in BSC formation and characterization of these organisms represents a crucial step towards understanding how anthropogenic change will affect the formation and ecological function of BSCs in arid ecosystems.