Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes.

Mobile genetic elements (MGEs) like plasmids, viruses, and transposable elements can provide fitness benefits to their hosts for survival in the presence of environmental stressors. Heavy metal resistance genes (HMRGs) are frequently observed on MGEs, suggesting that MGEs may be an important driver of adaptive evolution in environments contaminated with heavy metals. Here, we report the meta-mobilome of the heavy metal-contaminated regions of the Oak Ridge Reservation subsurface. This meta-mobilome was compared with one derived from samples collected from unimpacted regions of the Oak Ridge Reservation subsurface. We assembled 1615 unique circularized DNA elements that we propose to be MGEs. The circular elements from the highly contaminated subsurface were enriched in HMRG clusters relative to those from the nearby unimpacted regions. Additionally, we found that these HMRGs were associated with Gamma and Betaproteobacteria hosts in the contaminated subsurface and potentially facilitate the persistence and dominance of these taxa in this region. Finally, the HMRGs were associated with conjugative elements, suggesting their potential for future lateral transfer. We demonstrate how our understanding of MGE ecology, evolution, and function can be enhanced through the genomic context provided by completed MGE assemblies.

Genomic and environmental controls on Castellaniella biogeography in an anthropogenically disturbed subsurface.

Castellaniella species have been isolated from a variety of mixed-waste environments including the nitrate and multiple metal-contaminated subsurface at the Oak Ridge Reservation (ORR). Previous studies examining microbial community composition and nitrate removal at ORR during biostimulation efforts reported increased abundances of members of the Castellaniella genus concurrent with increased denitrification rates. Thus, we asked how genomic and abiotic factors control the Castellaniella biogeography at the site to understand how these factors may influence nitrate transformation in an anthropogenically impacted setting. We report the isolation and characterization of several Castellaniella strains from the ORR subsurface. Five of these isolates match at 100% identity (at the 16S rRNA gene V4 region) to two Castellaniella amplicon sequence variants (ASVs), ASV1 and ASV2, that have persisted in the ORR subsurface for at least 2 decades. However, ASV2 has consistently higher relative abundance in samples taken from the site and was also the dominant blooming denitrifier population during a prior biostimulation effort. We found that the ASV2 representative strain has greater resistance to mixed metal stress than the ASV1 representative strains. We attribute this resistance, in part, to the large number of unique heavy metal resistance genes identified on a genomic island in the ASV2 representative genome. Additionally, we suggest that the relatively lower fitness of ASV1 may be connected to the loss of the nitrous oxide reductase (nos) operon (and associated nitrous oxide reductase activity) due to the insertion at this genomic locus of a mobile genetic element carrying copper resistance genes. This study demonstrates the value of integrating genomic, environmental, and phenotypic data to characterize the biogeography of key microorganisms in contaminated sites.

Contribution of Microorganisms with the Clade II Nitrous Oxide Reductase to Suppression of Surface Emissions of Nitrous Oxide.

The sources and sinks of nitrous oxide, as control emissions to the atmosphere, are generally poorly constrained for most environmental systems. Initial depth-resolved analysis of nitrous oxide flux from observation wells and the proximal surface within a nitrate contaminated aquifer system revealed high subsurface production but little escape from the surface. To better understand the environmental controls of production and emission at this site, we used a combination of isotopic, geochemical, and molecular analyses to show that chemodenitrification and bacterial denitrification are major sources of nitrous oxide in this subsurface, where low DO, low pH, and high nitrate are correlated with significant nitrous oxide production. Depth-resolved metagenomes showed that consumption of nitrous oxide near the surface was correlated with an enrichment of Clade II nitrous oxide reducers, consistent with a growing appreciation of their importance in controlling release of nitrous oxide to the atmosphere. Our work also provides evidence for the reduction of nitrous oxide at a pH of 4, well below the generally accepted limit of pH 5.

A fast comparative genome browser for diverse bacteria and archaea.

Genome sequencing has revealed an incredible diversity of bacteria and archaea, but there are no fast and convenient tools for browsing across these genomes. It is cumbersome to view the prevalence of homologs for a protein of interest, or the gene neighborhoods of those homologs, across the diversity of the prokaryotes. We developed a web-based tool, fast.genomics, that uses two strategies to support fast browsing across the diversity of prokaryotes. First, the database of genomes is split up. The main database contains one representative from each of the 6,377 genera that have a high-quality genome, and additional databases for each taxonomic order contain up to 10 representatives of each species. Second, homologs of proteins of interest are identified quickly by using accelerated searches, usually in a few seconds. Once homologs are identified, fast.genomics can quickly show their prevalence across taxa, view their neighboring genes, or compare the prevalence of two different proteins. Fast.genomics is available at https://fast.genomics.lbl.gov.

Environmental stress mediates groundwater microbial community assembly.

Community assembly describes how different ecological processes shape microbial community composition and structure. How environmental factors impact community assembly remains elusive. Here we sampled microbial communities and >200 biogeochemical variables in groundwater at the Oak Ridge Field Research Center, a former nuclear waste disposal site, and developed a theoretical framework to conceptualize the relationships between community assembly processes and environmental stresses. We found that stochastic assembly processes were critical (>60% on average) in shaping community structure, but their relative importance decreased as stress increased. Dispersal limitation and 'drift' related to random birth and death had negative correlations with stresses, whereas the selection processes leading to dissimilar communities increased with stresses, primarily related to pH, cobalt and molybdenum. Assembly mechanisms also varied greatly among different phylogenetic groups. Our findings highlight the importance of microbial dispersal limitation and environmental heterogeneity in ecosystem restoration and management.

Eliminating genes for a two-component system increases PHB productivity in Cupriavidus basilensis 4G11 under PHB suppressing, nonstress conditions.

Species of bacteria from the genus Cupriavidus are known, in part, for their ability to produce high amounts of poly-hydroxybutyrate (PHB) making them attractive candidates for bioplastic production. The native synthesis of PHB occurs during periods of metabolic stress, and the process regulating the initiation of PHB accumulation in these organisms is not fully understood. Screening an RB-TnSeq transposon library of Cupriavidus basilensis 4G11 allowed us to identify two genes of an apparent, uncharacterized two-component system, which when omitted from the genome enable increased PHB productivity in balanced, nonstress growth conditions. We observe average increases in PHB productivity of 56% and 41% relative to the wildtype parent strain upon deleting each gene individually from the genome. The increased PHB phenotype disappears, however, in nitrogen-free unbalanced growth conditions suggesting the phenotype is specific to fast-growing, replete, nonstress growth. Bioproduction modeling suggests this phenotype could be due to a decreased reliance on metabolic stress induced by nitrogen limitation to initiate PHB production in the mutant strains. Due to uncertainty in the two-component system's input signal and regulon, the mechanism by which these genes impart this phenotype remains unclear. Such strains may allow for the use of single-stage, continuous bioreactor systems, which are far simpler than many PHB bioproduction schemes used previously, given a similar product yield to batch systems in such a configuration. Bioproductivity modeling suggests that omitting this regulation in the cells may increase PHB productivity up to 24% relative to the wildtype organism when using single-stage continuous systems. This work expands our understanding of the regulation of PHB accumulation in Cupriavidus, in particular the initiation of this process upon transition into unbalanced growth regimes.

A bacterial sensor taxonomy across earth ecosystems for machine learning applications.

Microbial communities have evolved to colonize all ecosystems of the planet, from the deep sea to the human gut. Microbes survive by sensing, responding, and adapting to immediate environmental cues. This process is driven by signal transduction proteins such as histidine kinases, which use their sensing domains to bind or otherwise detect environmental cues and "transduce" signals to adjust internal processes. We hypothesized that an ecosystem's unique stimuli leave a sensor "fingerprint," able to identify and shed insight on ecosystem conditions. To test this, we collected 20,712 publicly available metagenomes from Host-associated, Environmental, and Engineered ecosystems across the globe. We extracted and clustered the collection's nearly 18M unique sensory domains into 113,712 similar groupings with MMseqs2. We built gradient-boosted decision tree machine learning models and found we could classify the ecosystem type (accuracy: 87%) and predict the levels of different physical parameters (R2 score: 83%) using the sensor cluster abundance as features. Feature importance enables identification of the most predictive sensors to differentiate between ecosystems which can lead to mechanistic interpretations if the sensor domains are well annotated. To demonstrate this, a machine learning model was trained to predict patient's disease state and used to identify domains related to oxygen sensing present in a healthy gut but missing in patients with abnormal conditions. Moreover, since 98.7% of identified sensor domains are uncharacterized, importance ranking can be used to prioritize sensors to determine what ecosystem function they may be sensing. Furthermore, these new predictive sensors can function as targets for novel sensor engineering with applications in biotechnology, ecosystem maintenance, and medicine.IMPORTANCEMicrobes infect, colonize, and proliferate due to their ability to sense and respond quickly to their surroundings. In this research, we extract the sensory proteins from a diverse range of environmental, engineered, and host-associated metagenomes. We trained machine learning classifiers using sensors as features such that it is possible to predict the ecosystem for a metagenome from its sensor profile. We use the optimized model's feature importance to identify the most impactful and predictive sensors in different environments. We next use the sensor profile from human gut metagenomes to classify their disease states and explore which sensors can explain differences between diseases. The sensors most predictive of environmental labels here, most of which correspond to uncharacterized proteins, are a useful starting point for the discovery of important environment signals and the development of possible diagnostic interventions.

Systematic and scalable genome-wide essentiality mapping to identify nonessential genes in phages.

Phages are one of the key ecological drivers of microbial community dynamics, function, and evolution. Despite their importance in bacterial ecology and evolutionary processes, phage genes are poorly characterized, hampering their usage in a variety of biotechnological applications. Methods to characterize such genes, even those critical to the phage life cycle, are labor intensive and are generally phage specific. Here, we develop a systematic gene essentiality mapping method scalable to new phage-host combinations that facilitate the identification of nonessential genes. As a proof of concept, we use an arrayed genome-wide CRISPR interference (CRISPRi) assay to map gene essentiality landscape in the canonical coliphages λ and P1. Results from a single panel of CRISPRi probes largely recapitulate the essential gene roster determined from decades of genetic analysis for lambda and provide new insights into essential and nonessential loci in P1. We present evidence of how CRISPRi polarity can lead to false positive gene essentiality assignments and recommend caution towards interpreting CRISPRi data on gene essentiality when applied to less studied phages. Finally, we show that we can engineer phages by inserting DNA barcodes into newly identified inessential regions, which will empower processes of identification, quantification, and tracking of phages in diverse applications.

Evaluating E. coli genome-scale metabolic model accuracy with high-throughput mutant fitness data.

The Escherichia coli genome-scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision-recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene-protein-reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high-throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.

Mixed nitrate and metal contamination influences operational speciation of toxic and essential elements.

Environmental contamination constrains microbial communities impacting diversity and total metabolic activity. The former S-3 Ponds contamination site at Oak Ridge Reservation (ORR), TN, has elevated concentrations of nitric acid and multiple metals from decades of processing nuclear material. To determine the nature of the metal contamination in the sediment, a three-step sequential chemical extraction (BCR) was performed on sediment segments from a core located upgradient (EB271, non-contaminated) and one downgradient (EB106, contaminated) of the S-3 Ponds. The resulting exchangeable, reducing, and oxidizing fractions were analyzed for 18 different elements. Comparison of the two cores revealed changes in operational speciation for several elements caused by the contamination. Those present from the S-3 Ponds, including Al, U, Co, Cu, Ni, and Cd, were not only elevated in concentration in the EB106 core but were also operationally more available with increased mobility in the acidic environment. Other elements, including Mg, Ca, P, V, As, and Mo, were less operationally available in EB106 having decreased concentrations in the exchangeable fraction. The bioavailability of essential macro nutrients Mg, Ca, and P from the two types of sediment was determined using three metal-tolerant bacteria previously isolated from ORR. Mg and Ca were available from both sediments for all three strains; however, P was not bioavailable from either sediment for any strain. The decreased operational speciation of P in contaminated ORR sediment may increase the dependence of the microbial community on other pools of P or select for microorganisms with increased P scavenging capabilities. Hence, the microbial community at the former S-3 Ponds contamination site may be constrained not only by increased toxic metal concentrations but also by the availability of essential elements, including P.

Geochemical constraints on bacteriophage infectivity in terrestrial environments.

Lytic phages can be potent and selective inhibitors of microbial growth and can have profound impacts on microbiome composition and function. However, there is uncertainty about the biogeochemical conditions under which phage predation modulates microbial ecosystem function, particularly in terrestrial systems. Ionic strength is critical for infection of bacteria by many phages, but quantitative data is limited on the ion thresholds for phage infection that can be compared with environmental ion concentrations. Similarly, while carbon composition varies in the environment, we do not know how this variability influences the impact of phage predation on microbiome function. Here, we measured the half-maximal effective concentrations (EC50) of 80 different inorganic ions for the infection of E. coli with two canonical dsDNA and ssRNA phages, T4 and MS2, respectively. Many alkaline earth metals and alkali metals enabled lytic infection but the ionic strength thresholds varied for different ions between phages. Additionally, using a freshwater nitrate-reducing microbiome, we found that the ability of lytic phages to influence nitrate reduction end-products depended upon the carbon source as well as ionic strength. For all phage:host pairs, the ion EC50s for phage infection exceeded the ion concentrations found in many terrestrial freshwater systems. Thus, our findings support a model where phages most influence terrestrial microbial functional ecology in hot spots and hot moments such as metazoan guts, drought influenced soils, or biofilms where ion concentration is locally or transiently elevated and nutrients are available to support the growth of specific phage hosts.

Distinct Depth-Discrete Profiles of Microbial Communities and Geochemical Insights in the Subsurface Critical Zone.

Microbial assembly and metabolic potential in the subsurface critical zone (SCZ) are substantially impacted by subsurface geochemistry and hydrogeology, selecting for microbes distinct from those in surficial soils. In this study, we integrated metagenomics and geochemistry to elucidate how microbial composition and metabolic potential are shaped and impacted by vertical variations in geochemistry and hydrogeology in terrestrial subsurface sediment. A sediment core from an uncontaminated, pristine well at Oak Ridge Field Research Center in Oak Ridge, Tennessee, including the shallow subsurface, vadose zone, capillary fringe, and saturated zone, was used in this study. Our results showed that subsurface microbes were highly localized and that communities were rarely interconnected. Microbial community composition as well as metabolic potential in carbon and nitrogen cycling varied even over short vertical distances. Further analyses indicated a strong depth-related covariation of community composition with a subset of 12 environmental variables. An analysis of dissolved organic carbon (DOC) quality via ultrahigh resolution mass spectrometry suggested that the SCZ was generally a low-carbon environment, with the relative portion of labile DOC decreasing and that of recalcitrant DOC increasing along the depth, selecting microbes from copiotrophs to oligotrophs and also impacting the microbial metabolic potential in the carbon cycle. Our study demonstrates that sediment geochemistry and hydrogeology are vital in the selection of distinct microbial populations and metabolism in the SCZ. IMPORTANCE In this study, we explored the links between geochemical parameters, microbial community structure and metabolic potential across the depth of sediment, including the shallow subsurface, vadose zone, capillary fringe, and saturated zone. Our results revealed that microbes in the terrestrial subsurface can be highly localized, with communities rarely being interconnected along the depth. Overall, our research demonstrates that sediment geochemistry and hydrogeology are vital in the selection of distinct microbial populations and metabolic potential in different depths of subsurface terrestrial sediment. Such studies correlating microbial community analyses and geochemistry analyses, including high resolution mass spectrometry analyses of natural organic carbon, will further the fundamental understanding of microbial ecology and biogeochemistry in subsurface terrestrial ecosystems and will benefit the future development of predictive models on nutrient turnover in these environments.

Large-scale genetic characterization of the model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough.

Sulfate-reducing bacteria (SRB) are obligate anaerobes that can couple their growth to the reduction of sulfate. Despite the importance of SRB to global nutrient cycles and their damage to the petroleum industry, our molecular understanding of their physiology remains limited. To systematically provide new insights into SRB biology, we generated a randomly barcoded transposon mutant library in the model SRB Desulfovibrio vulgaris Hildenborough (DvH) and used this genome-wide resource to assay the importance of its genes under a range of metabolic and stress conditions. In addition to defining the essential gene set of DvH, we identified a conditional phenotype for 1,137 non-essential genes. Through examination of these conditional phenotypes, we were able to make a number of novel insights into our molecular understanding of DvH, including how this bacterium synthesizes vitamins. For example, we identified DVU0867 as an atypical L-aspartate decarboxylase required for the synthesis of pantothenic acid, provided the first experimental evidence that biotin synthesis in DvH occurs via a specialized acyl carrier protein and without methyl esters, and demonstrated that the uncharacterized dehydrogenase DVU0826:DVU0827 is necessary for the synthesis of pyridoxal phosphate. In addition, we used the mutant fitness data to identify genes involved in the assimilation of diverse nitrogen sources and gained insights into the mechanism of inhibition of chlorate and molybdate. Our large-scale fitness dataset and RB-TnSeq mutant library are community-wide resources that can be used to generate further testable hypotheses into the gene functions of this environmentally and industrially important group of bacteria.

Engineering osmolysis susceptibility in Cupriavidus necator and Escherichia coli for recovery of intracellular products.

BACKGROUND: Intracellular biomacromolecules, such as industrial enzymes and biopolymers, represent an important class of bio-derived products obtained from bacterial hosts. A common key step in the downstream separation of these biomolecules is lysis of the bacterial cell wall to effect release of cytoplasmic contents. Cell lysis is typically achieved either through mechanical disruption or reagent-based methods, which introduce issues of energy demand, material needs, high costs, and scaling problems. Osmolysis, a cell lysis method that relies on hypoosmotic downshock upon resuspension of cells in distilled water, has been applied for bioseparation of intracellular products from extreme halophiles and mammalian cells. However, most industrial bacterial strains are non-halotolerant and relatively resistant to hypoosmotic cell lysis. RESULTS: To overcome this limitation, we developed two strategies to increase the susceptibility of non-halotolerant hosts to osmolysis using Cupriavidus necator, a strain often used in electromicrobial production, as a prototypical strain. In one strategy, C. necator was evolved to increase its halotolerance from 1.5% to 3.25% (w/v) NaCl through adaptive laboratory evolution, and genes potentially responsible for this phenotypic change were identified by whole genome sequencing. The evolved halotolerant strain experienced an osmolytic efficiency of 47% in distilled water following growth in 3% (w/v) NaCl. In a second strategy, the cells were made susceptible to osmolysis by knocking out the large-conductance mechanosensitive channel (mscL) gene in C. necator. When these strategies were combined by knocking out the mscL gene from the evolved halotolerant strain, greater than 90% osmolytic efficiency was observed upon osmotic downshock. A modified version of this strategy was applied to E. coli BL21 by deleting the mscL and mscS (small-conductance mechanosensitive channel) genes. When grown in medium with 4% NaCl and subsequently resuspended in distilled water, this engineered strain experienced 75% cell lysis, although decreases in cell growth rate due to higher salt concentrations were observed. CONCLUSIONS: Our strategy is shown to be a simple and effective way to lyse cells for the purification of intracellular biomacromolecules and may be applicable in many bacteria used for bioproduction.

Microbial biomanufacturing for space-exploration-what to take and when to make.

As renewed interest in human space-exploration intensifies, a coherent and modernized strategy for mission design and planning has become increasingly crucial. Biotechnology has emerged as a promising approach to increase resilience, flexibility, and efficiency of missions, by virtue of its ability to effectively utilize in situ resources and reclaim resources from waste streams. Here we outline four primary mission-classes on Moon and Mars that drive a staged and accretive biomanufacturing strategy. Each class requires a unique approach to integrate biomanufacturing into the existing mission-architecture and so faces unique challenges in technology development. These challenges stem directly from the resources available in a given mission-class-the degree to which feedstocks are derived from cargo and in situ resources-and the degree to which loop-closure is necessary. As mission duration and distance from Earth increase, the benefits of specialized, sustainable biomanufacturing processes also increase. Consequentially, we define specific design-scenarios and quantify the usefulness of in-space biomanufacturing, to guide techno-economics of space-missions. Especially materials emerged as a potentially pivotal target for biomanufacturing with large impact on up-mass cost. Subsequently, we outline the processes needed for development, testing, and deployment of requisite technologies. As space-related technology development often does, these advancements are likely to have profound implications for the creation of a resilient circular bioeconomy on Earth.

Multicopy suppressor screens reveal convergent evolution of single-gene lysis proteins.

Single-strand RNA (ssRNA) Fiersviridae phages cause host lysis with a product of single gene (sgl for single-gene lysis; product Sgl) that induces autolysis. Many different Sgls have been discovered, but the molecular targets of only a few have been identified. In this study, we used a high-throughput genetic screen to uncover genome-wide host suppressors of diverse Sgls. In addition to validating known molecular mechanisms, we discovered that the Sgl of PP7, an ssRNA phage of Pseudomonas aeruginosa, targets MurJ, the flippase responsible for lipid II export, previously shown to be the target of the Sgl of coliphage M. These two Sgls, which are unrelated and predicted to have opposite membrane topology, thus represent a case of convergent evolution. We extended the genetic screens to other uncharacterized Sgls and uncovered a common set of multicopy suppressors, suggesting that these Sgls act by the same or similar mechanism.

Quantifying the local adaptive landscape of a nascent bacterial community.

The fitness effects of all possible mutations available to an organism largely shape the dynamics of evolutionary adaptation. Yet, whether and how this adaptive landscape changes over evolutionary times, especially upon ecological diversification and changes in community composition, remains poorly understood. We sought to fill this gap by analyzing a stable community of two closely related ecotypes ("L" and "S") shortly after they emerged within the E. coli Long-Term Evolution Experiment (LTEE). We engineered genome-wide barcoded transposon libraries to measure the invasion fitness effects of all possible gene knockouts in the coexisting strains as well as their ancestor, for many different, ecologically relevant conditions. We find consistent statistical patterns of fitness effect variation across both genetic background and community composition, despite the idiosyncratic behavior of individual knockouts. Additionally, fitness effects are correlated with evolutionary outcomes for a number of conditions, possibly revealing shifting patterns of adaptation. Together, our results reveal how ecological and epistatic effects combine to shape the adaptive landscape in a nascent ecological community.

Data-driven flow-map models for data-efficient discovery of dynamics and fast uncertainty quantification of biological and biochemical systems.

Computational models are increasingly used to investigate and predict the complex dynamics of biological and biochemical systems. Nevertheless, governing equations of a biochemical system may not be (fully) known, which would necessitate learning the system dynamics directly from, often limited and noisy, observed data. On the other hand, when expensive models are available, systematic and efficient quantification of the effects of model uncertainties on quantities of interest can be an arduous task. This paper leverages the notion of flow-map (de)compositions to present a framework that can address both of these challenges via learning data-driven models useful for capturing the dynamical behavior of biochemical systems. Data-driven flow-map models seek to directly learn the integration operators of the governing differential equations in a black-box manner, irrespective of structure of the underlying equations. As such, they can serve as a flexible approach for deriving fast-to-evaluate surrogates for expensive computational models of system dynamics, or, alternatively, for reconstructing the long-term system dynamics via experimental observations. We present a data-efficient approach to data-driven flow-map modeling based on polynomial chaos Kriging. The approach is demonstrated for discovery of the dynamics of various benchmark systems and a coculture bioreactor subject to external forcing, as well as for uncertainty quantification of a microbial electrosynthesis reactor. Such data-driven models and analyses of dynamical systems can be paramount in the design and optimization of bioprocesses and integrated biomanufacturing systems.