Diversity, biogeography and assembly mechanisms of groundwater microbiomes
Explore the source record for details and available documents.
Wu, Zongzhi (ORCID:0000000251535246)↗
Engineering topics
Arkin, Adam P.
Publications and source records attributed to Arkin, Adam P..
Fitness factors impacting survival of a subsurface bacterium in contaminated groundwater
Many factors contribute to the ability of a microbial species to persist when encountering complexly contaminated environments, including time of exposure, the nature and concentration of contaminants, availability of nutritional resources, and possession of a combination of appropriate molecular mechanisms needed for survival. Herein we sought to identify genes that are most important for survival of Gram-negative Enterobacteriaceae in contaminated groundwater environments containing high concentrations of nitrate and metals using the metal-tolerant Oak Ridge Reservation isolate, Pantoea sp. MT58 (MT58). Survival fitness experiments in which a randomly barcoded transposon insertion (RB-TnSeq) library of MT58 was exposed directly to contaminated Oak Ridge Reservation groundwater samples from across a nitrate and mixed metal contamination plume were used to identify genes important for survival with increasing exposure times and concentrations of contaminants, and availability of a carbon source. Genes involved in controlling and using carbon, encoding transcriptional regulators, and related to Gram-negative outer membrane processes were among those found to be important for survival in contaminated Oak Ridge Reservation groundwater. A comparative genomics analysis of 75 Pantoea genus strains allowed us to further separate the survival determinants into core and non-core genes in the Pantoea pangenome, revealing insights into the survival of subsurface microorganisms during contaminant plume intrusion.
59 BASIC BIOLOGICAL SCIENCES↗
Interactive tools for functional annotation of bacterial genomes
Automated annotations of protein functions are error-prone because of our lack of knowledge of protein functions. For example, it is often impossible to predict the correct substrate for an enzyme or a transporter. Furthermore, much of the knowledge that we do have about the functions of proteins is missing from the underlying databases. We discuss how to use interactive tools to quickly find different kinds of information relevant to a protein’s function. Many of these tools are available via PaperBLAST (http://papers.genomics.lbl.gov). Combining these tools often allows us to infer a protein’s function. Ideally, accurate annotations would allow us to predict a bacterium’s capabilities from its genome sequence, but in practice, this remains challenging. We describe interactive tools that infer potential capabilities from a genome sequence or that search a genome to find proteins that might perform a specific function of interest.
59 BASIC BIOLOGICAL SCIENCES↗
High-throughput genetics enables identification of nutrient utilization and accessory energy metabolism genes in a model methanogen
Archaea are widespread in the environment and play fundamental roles in diverse ecosystems; however, characterization of their unique biology requires advanced tools. This is particularly challenging when characterizing gene function. Here, we generate randomly barcoded transposon libraries in the model methanogenic archaeon Methanococcus maripaludis and use high-throughput growth methods to conduct fitness assays (RB-TnSeq) across over 100 unique growth conditions. Using our approach, we identified new genes involved in nutrient utilization and response to oxidative stress. We identified novel genes for the usage of diverse nitrogen sources in M. maripaludis including a putative regulator of alanine deamination and molybdate transporters important for nitrogen fixation. Furthermore, leveraging the fitness data, we inferred that M. maripaludis can utilize additional nitrogen sources including $\tiny{L}$-glutamine, $\tiny{D}$-glucuronamide, and adenosine. Under autotrophic growth conditions, we identified a gene encoding a domain of unknown function (DUF166) that is important for fitness and hypothesize that it has an accessory role in carbon dioxide assimilation. Finally, comparing fitness costs of oxygen versus sulfite stress, we identified a previously uncharacterized class of dissimilatory sulfite reductase-like proteins (Dsr-LP; group IIId) that is important during growth in the presence of sulfite. When overexpressed, Dsr-LP conferred sulfite resistance and enabled use of sulfite as the sole sulfur source. The high-throughput approach employed here allowed for generation of a large-scale data set that can be used as a resource to further understand gene function and metabolism in the archaeal domain.
59 BASIC BIOLOGICAL SCIENCES↗
Molecular mechanisms and environmental adaptations of flagellar loss and biofilm growth of Rhodanobacter under environmental stress
Biofilms aid bacterial adhesion to surfaces via direct and indirect mechanisms, and formation of biofilms is considered as an important strategy for adaptation and survival in suboptimal environmental conditions. However, the molecular underpinnings of biofilm formation in subsurface sediment/groundwater ecosystems where microorganisms often experience fluctuations in nutrient input, pH, and nitrate or metal concentrations are underexplored. Here, we examined biofilm formation under different nutrient, pH, metal, and nitrate regimens of 16 Rhodanobacter strains isolated from subsurface groundwater wells spanning diverse levels of pH (3.5 to 5) and nitrates (13.7 to 146 mM). Eight Rhodanobacter strains demonstrated significant biofilm growth under low pH, suggesting adaptations for survival and growth at low pH. Biofilms were intensified under aluminum stress, particularly in strains possessing fewer genetic traits associated with biofilm formation, findings warranting further investigation. Through random barcode transposon-site sequencing (RB-TnSeq), proteomics, use of specific mutants, and transmission electron microscopy analysis, we discovered flagellar loss under aluminum stress, indicating a potential relationship between motility, metal tolerance, and biofilm growth. Comparative genomic analyses revealed the absence of flagella and chemotaxis genes and the presence of a putative type VI secretion system in the highly biofilm-forming strain FW021-MT20. In this study we identified genetic determinants associated with biofilm growth under metal stress in a predominant environmental genus, Rhodanobacter, and identified traits aiding survival and adaptation to contaminated subsurface environments.
59 BASIC BIOLOGICAL SCIENCES↗
Multiplex characterization of microbial traits using dual barcoded nucleic acid fragment expression library
Disclosed herein are barcoded expression libraries comprising a plurality of expression vectors, wherein each expression vector comprises a nucleic acid fragment flanked by a first barcode and a second barcode. Further disclosed herein are methods of making the barcoded expression libraries and methods of conducting functional analysis using the barcoded expression libraries.
Mutalik, Vivek K.↗
Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes
Abstract 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.
59 BASIC BIOLOGICAL SCIENCES↗
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.
59 BASIC BIOLOGICAL SCIENCES↗
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. Further, 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.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 .
59 BASIC BIOLOGICAL SCIENCES↗
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. As a result, our findings highlight the importance of microbial dispersal limitation and environmental heterogeneity in ecosystem restoration and management.
54 ENVIRONMENTAL SCIENCES↗
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.
97 MATHEMATICS AND COMPUTING↗
Evaluating E. coli genome‐scale metabolic model accuracy with high‐throughput mutant fitness data
Abstract 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.
59 BASIC BIOLOGICAL SCIENCES↗
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. Here, 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.
54 ENVIRONMENTAL SCIENCES↗
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.
59 BASIC BIOLOGICAL SCIENCES↗
Engineering osmolysis susceptibility in Cupriavidus necator and Escherichia coli for recovery of intracellular products
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. 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. 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.
59 BASIC BIOLOGICAL SCIENCES↗
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.
59 BASIC BIOLOGICAL SCIENCES↗
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.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗