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.

An expanded transcriptome atlas for Bacteroides thetaiotaomicron reveals a small RNA that modulates tetracycline sensitivity

Plasticity in gene expression allows bacteria to adapt to diverse environments. This is particularly relevant in the dynamic niche of the human intestinal tract; however, transcriptional networks remain largely unknown for gut-resident bacteria. Here we apply differential RNA sequencing (RNA-seq) and conventional RNA-seq to the model gut bacterium Bacteroides thetaiotaomicron to map transcriptional units and profile their expression levels across 15 in vivo-relevant growth conditions. We infer stress- and carbon source-specific transcriptional regulons and expand the annotation of small RNAs (sRNAs). Integrating this expression atlas with published transposon mutant fitness data, we predict conditionally important sRNAs. These include MasB, which downregulates tetracycline tolerance. Using MS2 affinity purification and RNA-seq, we identify a putative MasB target and assess its role in the context of the MasB-associated phenotype. These data-publicly available through the Theta-Base web browser ( http://micromix.helmholtz-hiri.de/bacteroides/ )-constitute a valuable resource for the microbiome community.

A weaponized phage suppresses competitors in historical and modern metapopulations of pathogenic bacteria

Bacteriophages, the viruses of bacteria, are proposed to drive bacterial population dynamics, yet direct evidence of their impact on natural populations is limited. Here we identified viral sequences in a metapopulation of wild plant-associated Pseudomonas spp. genomes. We discovered that the most abundant viral cluster does not encode an intact phage but instead encodes a tailocin - a phage-derived element that bacteria use to kill competitors for interbacterial warfare. Each pathogenic Pseudomonas sp. strain carries one of a few distinct tailocin variants, which target variable polysaccharides in the outer membrane of co-occurring pathogenic strains. Analysis of historic herbarium samples from the last 170 years revealed that the same tailocin and receptor variants have persisted in the Pseudomonas populations for at least two centuries, suggesting the continued use of a defined set of tailocin haplotypes and receptors. These results indicate that tailocin genetic diversity can be mined to develop targeted "tailocin cocktails" for microbial control.

One-Sentence Summary

Bacterial pathogens in a host-associated metapopulation use a repurposed prophage to kill their competitors.

Randomly barcoded transposon mutant libraries for gut commensals II: Applying libraries for functional genetics

The critical role of the intestinal microbiota in human health and disease is well recognized. Nevertheless, there are still large gaps in our understanding of the functions and mechanisms encoded in the genomes of most members of the gut microbiota. Genome-scale libraries of transposon mutants are a powerful tool to help us address this gap. Recent advances in barcoded transposon mutagenesis have dramatically lowered the cost of mutant fitness determination in hundreds of in vitro and in vivo experimental conditions. In an accompanying review, we discuss recent advances and caveats for the construction of pooled and arrayed barcoded transposon mutant libraries in human gut commensals. In this review, we discuss how these libraries can be used across a wide range of applications, the technical aspects involved, and expectations for such screens.

Randomly barcoded transposon mutant libraries for gut commensals I: Strategies for efficient library construction

Randomly barcoded transposon mutant libraries are powerful tools for studying gene function and organization, assessing gene essentiality and pathways, discovering potential therapeutic targets, and understanding the physiology of gut bacteria and their interactions with the host. However, construction of high-quality libraries with uniform representation can be challenging. In this review, we survey various strategies for barcoded library construction, including transposition systems, methods of transposon delivery, optimal library size, and transconjugant selection schemes. We discuss the advantages and limitations of each approach, as well as factors to consider when selecting a strategy. In addition, we highlight experimental and computational advances in arraying condensed libraries from mutant pools. We focus on examples of successful library construction in gut bacteria and their application to gene function studies and drug discovery. Given the need for understanding gene function and organization in gut bacteria, we provide a comprehensive guide for researchers to construct randomly barcoded transposon mutant libraries.

A chemically-defined growth medium to support Lactobacillus-Acetobacter sp. community analysis

Lactobacilli and Acetobacter sp. are commercially important bacteria that often form communities in natural fermentations, including food preparations, spoilage, and in the digestive tract of the fruit fly Drosophila melanogaster. Communities of these bacteria are widespread and prolific, despite numerous strain-specific auxotrophies, suggesting they have evolved nutrient interdependencies that regulate their growth. The use of a chemically-defined medium (CDM) supporting the growth of both groups of bacteria would facilitate the identification of the molecular mechanisms for the metabolic interactions between them. While numerous CDMs have been developed that support specific strains of lactobacilli or Acetobacter, there has not been a medium formulated to support both genera. We developed such a medium, based on a previous CDM designed for growth of lactobacilli, by modifying the nutrient abundances to improve growth yield. We further simplified the medium by substituting casamino acids in place of individual amino acids and the standard Wolfe's vitamins and mineral stocks in place of individual vitamins and minerals, resulting in a reduction from 40 to 8 stock solutions. These stock solutions can be used to prepare several CDM formulations that support robust growth of numerous lactobacilli and Acetobacters. Here, we provide the composition and several examples of its use, which is important for tractability in dissecting the genetic and metabolic basis of natural bacterial species interactions.

Simvastatin induces human gut bacterial cell surface genes

Drugs intended to target mammalian cells can have broad off-target effects on the human gut microbiota with potential downstream consequences for drug efficacy and side effect profiles. Yet, despite a rich literature on antibiotic resistance, we still know very little about the mechanisms through which commensal bacteria evade non-antibiotic drugs. Here, we focus on statins, one of the most prescribed drug types in the world and an essential tool in the prevention and treatment of high circulating cholesterol levels. Prior work in humans, mice, and cell culture support an off-target effect of statins on human gut bacteria; however, the genetic determinants of statin sensitivity remain unknown. We confirmed that simvastatin inhibits the growth of diverse human gut bacterial strains grown in communities and in pure cultures. Drug sensitivity varied between phyla and was dose-dependent. We selected two representative simvastatin-sensitive species for more in-depth analysis: Eggerthella lenta (phylum: Actinobacteriota) and Bacteroides thetaiotaomicron (phylum: Bacteroidota). Transcriptomics revealed that both bacterial species upregulate genes in response to simvastatin that alter the cell membrane, including fatty acid biogenesis (E. lenta) and drug efflux systems (B. thetaiotaomicron). Transposon mutagenesis identified a key efflux system in B. thetaiotaomicron that enables growth in the presence of statins. Taken together, these results emphasize the importance of the bacterial cell membrane in countering the off-target effects of host-targeted drugs. Continued mechanistic dissection of the various mechanisms through which the human gut microbiota evades drugs will be essential to understand and predict the effects of drug administration in human cohorts and the potential downstream consequences for health and disease.

A mutant fitness compendium in Bifidobacteria reveals molecular determinants of colonization and host-microbe interactions

Bifidobacteria commonly represent a dominant constituent of human gut microbiomes during infancy, influencing nutrition, immune development, and resistance to infection. Despite interest as a probiotic therapy, predicting the nutritional requirements and health-promoting effects of Bifidobacteria is challenging due to major knowledge gaps. To overcome these deficiencies, we used large-scale genetics to create a compendium of mutant fitness in Bifidobacterium breve (Bb). We generated a high density, randomly barcoded transposon insertion pool in Bb, and used this pool to determine Bb fitness requirements during colonization of germ-free mice and chickens with multiple diets and in response to hundreds of in vitro perturbations. To enable mechanistic investigation, we constructed an ordered collection of insertion strains covering 1462 genes. We leveraged these tools to improve models of metabolic pathways, reveal unexpected host- and diet-specific requirements for colonization, and connect the production of immunomodulatory molecules to growth benefits. These resources will greatly reduce the barrier to future investigations of this important beneficial microbe.

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.

Genome-wide fitness profiling reveals molecular mechanisms that bacteria use to interact with Trichoderma atroviride exometabolites

Trichoderma spp. are ubiquitous rhizosphere fungi capable of producing several classes of secondary metabolites that can modify the dynamics of the plant-associated microbiome. However, the bacterial-fungal mechanisms that mediate these interactions have not been fully characterized. Here, a random barcode transposon-site sequencing (RB-TnSeq) approach was employed to identify bacterial genes important for fitness in the presence of Trichoderma atroviride exudates. We selected three rhizosphere bacteria with RB-TnSeq mutant libraries that can promote plant growth: the nitrogen fixers Klebsiella michiganensis M5aI and Herbaspirillum seropedicae SmR1, and Pseudomonas simiae WCS417. As a non-rhizosphere species, Pseudomonas putida KT2440 was also included. From the RB-TnSeq data, nitrogen-fixing bacteria competed mainly for iron and required the siderophore transport system TonB/ExbB for optimal fitness in the presence of T. atroviride exudates. In contrast, P. simiae and P. putida were highly dependent on mechanisms associated with membrane lipid modification that are required for resistance to cationic antimicrobial peptides (CAMPs). A mutant in the Hog1-MAP kinase (Δtmk3) gene of T. atroviride showed altered expression patterns of many nonribosomal peptide synthetase (NRPS) biosynthetic gene clusters with potential antibiotic activity. In contrast to exudates from wild-type T. atroviride, bacterial mutants containing lesions in genes associated with resistance to antibiotics did not show fitness defects when RB-TnSeq libraries were exposed to exudates from the Δtmk3 mutant. Unexpectedly, exudates from wild-type T. atroviride and the Δtmk3 mutant rescued purine auxotrophic mutants of H. seropedicae, K. michiganensis and P. simiae. Metabolomic analysis on exudates from wild-type T. atroviride and the Δtmk3 mutant showed that both strains excrete purines and complex metabolites; functional Tmk3 is required to produce some of these metabolites. This study highlights the complex interplay between Trichoderma-metabolites and soil bacteria, revealing both beneficial and antagonistic effects, and underscoring the intricate and multifaceted nature of this relationship.

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.

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.

An integrated transcriptomics-functional genomics approach reveals a small RNA that modulates Bacteroides thetaiotaomicron sensitivity to tetracyclines

Gene expression plasticity allows bacteria to adapt to diverse environments, tie their metabolism to available nutrients, and cope with stress. This is particularly relevant in a niche as dynamic and hostile as the human intestinal tract, yet transcriptional networks remain largely unknown in gut Bacteroides spp. Here, we map transcriptional units and profile their expression levels in Bacteroides thetaiotaomicron over a suite of 15 defined experimental conditions that are relevant in vivo , such as variation of temperature, pH, and oxygen tension, exposure to antibiotic stress, and growth on simple carbohydrates or on host mucin-derived glycans. Thereby, we infer stress- and carbon source-specific transcriptional regulons, including conditional expression of capsular polysaccharides and polysaccharide utilization loci, and expand the annotation of small regulatory RNAs (sRNAs) in this organism. Integrating this comprehensive expression atlas with transposon mutant fitness data, we identify conditionally important sRNAs. One example is MasB, whose inactivation led to increased bacterial tolerance of tetracyclines. Using MS2 affinity purification coupled with RNA sequencing, we predict targets of this sRNA and discuss their potential role in the context of the MasB-associated phenotype. Together, this transcriptomic compendium in combination with functional sRNA genomics-publicly available through a new iteration of the 'Theta-Base' web browser (www.helmholtz-hiri.de/en/datasets/bacteroides-v2)-constitutes a valuable resource for the microbiome and sRNA research communities alike.

Genome-wide identification of fitness determinants in the Xanthomonas campestris bacterial pathogen during early stages of plant infection

Plant diseases are an important threat to food production. While major pathogenicity determinants required for disease have been extensively studied, less is known on how pathogens thrive during host colonization, especially at early infection stages. Here, we used randomly barcoded-transposon insertion site sequencing (RB-TnSeq) to perform a genome-wide screen and identify key bacterial fitness determinants of the vascular pathogen Xanthomonas campestris pv campestris (Xcc) during infection of the cauliflower host plant (Brassica oleracea). This high-throughput analysis was conducted in hydathodes, the natural entry site of Xcc, in xylem sap and in synthetic media. Xcc did not face a strong bottleneck during hydathode infection. In total, 181 genes important for fitness were identified in plant-associated environments with functional enrichment in genes involved in metabolism but only few genes previously known to be involved in virulence. The biological relevance of 12 genes was independently confirmed by phenotyping single mutants. Notably, we show that XC_3388, a protein with no known function (DUF1631), plays a key role in the adaptation and virulence of Xcc possibly through c-di-GMP-mediated regulation. This study revealed yet unsuspected social behaviors adopted by Xcc individuals when confined inside hydathodes at early infection stages.

Caulobacter lipid A is conditionally dispensable in the absence of fur and in the presence of anionic sphingolipids

Lipid A, the membrane-anchored portion of lipopolysaccharide (LPS), is an essential component of the outer membrane (OM) of nearly all Gram-negative bacteria. Here we identify regulatory and structural factors that together render lipid A nonessential in Caulobacter crescentus. Mutations in the ferric uptake regulator fur allow Caulobacter to survive in the absence of either LpxC, which catalyzes an early step of lipid A synthesis, or CtpA, a tyrosine phosphatase homolog we find is needed for wild-type lipid A structure and abundance. Alterations in Fur-regulated processes, rather than iron status per se, underlie the ability to survive when lipid A synthesis is blocked. Fitness of lipid A-deficient Caulobacter requires an anionic sphingolipid, ceramide phosphoglycerate (CPG), which also mediates sensitivity to the antibiotic colistin. Our results demonstrate that, in an altered regulatory landscape, anionic sphingolipids can support the integrity of a lipid A-deficient OM.

Lipid A, the membrane-anchoring segment of lipopolysaccharide, is generally considered to be an essential component of the Gram-negative bacterial outer membrane. Zik et al. show that deletion of the transcriptional regulator fur and synthesis of the anionic sphingolipid ceramide phosphoglycerate enable Caulobacter crescentus to survive without lipid A.

A Defined Medium for Cultivation and Exometabolite Profiling of Soil Bacteria

Exometabolomics is an approach to assess how microorganisms alter, or react to their environments through the depletion and production of metabolites. It allows the examination of how soil microbes transform the small molecule metabolites within their environment, which can be used to study resource competition and cross-feeding. This approach is most powerful when used with defined media that enable tracking of all metabolites. However, microbial growth media have traditionally been developed for the isolation and growth of microorganisms but not metabolite utilization profiling through Liquid Chromatography Tandem Mass Spectrometry (LC-MS/MS). Here, we describe the construction of a defined medium, the Northen Lab Defined Medium (NLDM), that not only supports the growth of diverse soil bacteria but also is defined and therefore suited for exometabolomic experiments. Metabolites included in NLDM were selected based on their presence in R2A medium and soil, elemental stoichiometry requirements, as well as knowledge of metabolite usage by different bacteria. We found that NLDM supported the growth of 108 of the 110 phylogenetically diverse (spanning 36 different families) soil bacterial isolates tested and all of its metabolites were trackable through LC–MS/MS analysis. These results demonstrate the viability and utility of the constructed NLDM medium for growing and characterizing diverse microbial isolates and communities.

Filling gaps in bacterial catabolic pathways with computation and high-throughput genetics

To discover novel catabolic enzymes and transporters, we combined high-throughput genetic data from 29 bacteria with an automated tool to find gaps in their catabolic pathways. GapMind for carbon sources automatically annotates the uptake and catabolism of 62 compounds in bacterial and archaeal genomes. For the compounds that are utilized by the 29 bacteria, we systematically examined the gaps in GapMind's predicted pathways, and we used the mutant fitness data to find additional genes that were involved in their utilization. We identified novel pathways or enzymes for the utilization of glucosamine, citrulline, myo-inositol, lactose, and phenylacetate, and we annotated 299 diverged enzymes and transporters. We also curated 125 proteins from published reports. For the 29 bacteria with genetic data, GapMind finds high-confidence paths for 85% of utilized carbon sources. In diverse bacteria and archaea, 38% of utilized carbon sources have high-confidence paths, which was improved from 27% by incorporating the fitness-based annotations and our curation. GapMind for carbon sources is available as a web server (http://papers.genomics.lbl.gov/carbon) and takes just 30 seconds for the typical genome.

Species- and site-specific genome editing in complex bacterial communities

Understanding microbial gene functions relies on the application of experimental genetics in cultured microorganisms. However, the vast majority of bacteria and archaea remain uncultured, precluding the application of traditional genetic methods to these organisms and their interactions. Here, we characterize and validate a generalizable strategy for editing the genomes of specific organisms in microbial communities. We apply environmental transformation sequencing (ET-seq), in which nontargeted transposon insertions are mapped and quantified following delivery to a microbial community, to identify genetically tractable constituents. Next, DNA-editing all-in-one RNA-guided CRISPR-Cas transposase (DART) systems for targeted DNA insertion into organisms identified as tractable by ET-seq are used to enable organism- and locus-specific genetic manipulation in a community context. Using a combination of ET-seq and DART in soil and infant gut microbiota, we conduct species- and site-specific edits in several bacteria, measure gene fitness in a nonmodel bacterium and enrich targeted species. These tools enable editing of microbial communities for understanding and control.

The genetic basis of phage susceptibility, cross-resistance and host-range in Salmonella

Though bacteriophages (phages) are known to play a crucial role in bacterial fitness and virulence, our knowledge about the genetic basis of their interaction, cross-resistance and host-range is sparse. Here, we employed genome-wide screens in Salmonella enterica serovar Typhimurium to discover host determinants involved in resistance to eleven diverse lytic phages including four new phages isolated from a therapeutic phage cocktail. We uncovered 301 diverse host factors essential in phage infection, many of which are shared between multiple phages demonstrating potential cross-resistance mechanisms. We validate many of these novel findings and uncover the intricate interplay between RpoS, the virulence-associated general stress response sigma factor and RpoN, the nitrogen starvation sigma factor in phage cross-resistance. Finally, the infectivity pattern of eleven phages across a panel of 23 genome sequenced Salmonella strains indicates that additional constraints and interactions beyond the host factors uncovered here define the phage host range.

Biofilm Interaction Mapping and Analysis (BIMA) of Interspecific Interactions in Pseudomonas Co-culture Biofilms

Pseudomonas species are ubiquitous in nature and include numerous medically, agriculturally and technologically beneficial strains of which the interspecific interactions are of great interest for biotechnologies. Specifically, co-cultures containing Pseudomonas stutzeri have been used for bioremediation, biocontrol, aquaculture management and wastewater denitrification. Furthermore, the use of P. stutzeri biofilms, in combination with consortia-based approaches, may offer advantages for these processes. Understanding the interspecific interaction within biofilm co-cultures or consortia provides a means for improvement of current technologies. However, the investigation of biofilm-based consortia has been limited. We present an adaptable and scalable method for the analysis of macroscopic interactions (colony morphology, inhibition, and invasion) between colony-forming bacterial strains using an automated printing method followed by analysis of the genes and metabolites involved in the interactions. Using Biofilm Interaction Mapping and Analysis (BIMA), these interactions were investigated between P. stutzeri strain RCH2, a denitrifier isolated from chromium (VI) contaminated soil, and 13 other species of pseudomonas isolated from non-contaminated soil. One interaction partner, Pseudomonas fluorescens N1B4 was selected for mutant fitness profiling of a DNA-barcoded mutant library; with this approach four genes of importance were identified and the effects on interactions were evaluated with deletion mutants and mass spectrometry based metabolomics.

Systematic discovery of pseudomonad genetic factors involved in sensitivity to tailocins

Tailocins are bactericidal protein complexes produced by a wide variety of bacteria that kill closely related strains and may play a role in microbial community structure. Thanks to their high specificity, tailocins have been proposed as precision antibacterial agents for therapeutic applications. Compared to tailed phages, with whom they share an evolutionary and morphological relationship, bacterially produced tailocins kill their host upon production but producing strains display resistance to self-intoxication. Though lipopolysaccharide (LPS) has been shown to act as a receptor for tailocins, the breadth of factors involved in tailocin sensitivity, and the mechanisms behind resistance to self-intoxication, remain unclear. Here, we employed genome-wide screens in four non-model pseudomonads to identify mutants with altered fitness in the presence of tailocins produced by closely related pseudomonads. Our mutant screens identified O-antigen composition and display as most important in defining sensitivity to our tailocins. In addition, the screens suggest LPS thinning as a mechanism by which resistant strains can become more sensitive to tailocins. We validate many of these novel findings, and extend these observations of tailocin sensitivity to 130 genome-sequenced pseudomonads. This work offers insights into tailocin-bacteria interactions, informing the potential use of tailocins in microbiome manipulation and antibacterial therapy.

Microfabrication of a Chamber for High-Resolution, In Situ Imaging of the Whole Root for Plant-Microbe Interactions

Fabricated ecosystems (EcoFABs) offer an innovative approach to in situ examination of microbial establishment patterns around plant roots using nondestructive, high-resolution microscopy. Previously high-resolution imaging was challenging because the roots were not constrained to a fixed distance from the objective. Here, we describe a new 'Imaging EcoFAB' and the use of this device to image the entire root system of growing Brachypodium distachyon at high resolutions (20×, 40×) over a 3-week period. The device is capable of investigating root-microbe interactions of multimember communities. We examined nine strains of Pseudomonas simiae with different fluorescent constructs to B. distachyon and individual cells on root hairs were visible. Succession in the rhizosphere using two different strains of P. simiae was examined, where the second addition was shown to be able to establish in the root tissue. The device was suitable for imaging with different solid media at high magnification, allowing for the imaging of fungal establishment in the rhizosphere. Overall, the Imaging EcoFAB could improve our ability to investigate the spatiotemporal dynamics of the rhizosphere, including studies of fluorescently-tagged, multimember, synthetic communities.

Ecogenomics of Groundwater Phages Suggests Niche Differentiation Linked to Specific Environmental Tolerance

Viruses are ubiquitous microbiome components, shaping ecosystems via strain-specific predation, horizontal gene transfer and redistribution of nutrients through host lysis. Viral impacts are important in groundwater ecosystems, where microbes drive many nutrient fluxes and metabolic processes; however, little is known about the diversity of viruses in these environments. We analyzed four groundwater plasmidomes (the entire plasmid content of an environment) and identified 200 viral sequences, which clustered into 41 genus-level viral clusters (approximately equivalent to viral genera) including 9 known and 32 putative new genera. We used publicly available bacterial whole-genome sequences (WGS) and WGS from 261 bacterial isolates from this groundwater environment to identify potential viral hosts. We linked 76 of the 200 viral sequences to a range of bacterial phyla, the majority associated with Proteobacteria, followed by Firmicutes, Bacteroidetes, and Actinobacteria. The publicly available WGS enabled mapping bacterial hosts to several viral sequences. The WGS of groundwater isolates increased the depth of host prediction by allowing host identification at the strain level. The latter included 4 viruses that were almost entirely (>99% query coverage, >99% identity) identified as integrated in the genomes of Pseudomonas, Acidovorax, and Castellaniella strains, resulting in high-confidence host assignments. Lastly, 21 of these viruses carried putative auxiliary metabolite genes for metal and antibiotic resistance, which might drive their infection cycles and/or provide selective advantage to infected hosts. Exploring the groundwater virome provides a necessary foundation for integration of viruses into ecosystem models where they are key players in microbial adaption to environmental stress. IMPORTANCE To our knowledge, this is the first study to identify the bacteriophage distribution in a groundwater ecosystem shedding light on their prevalence and distribution across metal-contaminated and background sites. Our study is uniquely based on selective sequencing of solely the extrachromosomal elements of a microbiome followed by analysis for viral signatures, thus establishing a more focused approach for phage identifications. Using this method, we detected several novel phage genera along with those previously established. Our approach of using the whole-genome sequences of hundreds of bacterial isolates from the same site enabled us to make host assignments with high confidence, several at strain levels. Certain phage genes suggest that they provide an environment-specific selective advantage to their bacterial hosts. Our study lays the foundation for future research on directed phage isolations using specific bacterial host strains to further characterize groundwater phages, their life cycles, and their effects on groundwater microbiome and biogeochemistry.

Deletion Mutants, Archived Transposon Library, and Tagged Protein Constructs of the Model Sulfate-Reducing Bacterium Desulfovibrio vulgaris Hildenborough

The dissimilatory sulfate-reducing deltaproteobacterium Desulfovibrio vulgaris Hildenborough (ATCC 29579) was chosen by the research collaboration ENIGMA to explore tools and protocols for bringing this anaerobe to model status. Here, we describe a collection of genetic constructs generated by ENIGMA that are available to the research community.

Microbial Community Field Surveys Reveal Abundant Pseudomonas Population in Sorghum Rhizosphere Composed of Many Closely Related Phylotypes

While the root-associated microbiome is typically less diverse than the surrounding soil due to both plant selection and microbial competition for plant derived resources, it typically retains considerable complexity, harboring many hundreds of distinct bacterial species. Here, we report a time-dependent deviation from this trend in the rhizospheres of field grown sorghum. In this study, 16S rRNA amplicon sequencing was used to determine the impact of nitrogen fertilization on the development of the root-associated microbiomes of 10 sorghum genotypes grown in eastern Nebraska. We observed that early rhizosphere samples exhibit a significant reduction in overall diversity due to a high abundance of the bacterial genus Pseudomonas that occurred independent of host genotype in both high and low nitrogen fields and was not observed in the surrounding soil or associated root endosphere samples. When clustered at 97% identity, nearly all the Pseudomonas reads in this dataset were assigned to a single operational taxonomic unit (OTU); however, exact sequence variant (ESV)-level resolution demonstrated that this population comprised a large number of distinct Pseudomonas lineages. Furthermore, single-molecule long-read sequencing enabled high-resolution taxonomic profiling revealing further heterogeneity in the Pseudomonas lineages that was further confirmed using shotgun metagenomic sequencing. Finally, field soil enriched with specific carbon compounds recapitulated the increase in Pseudomonas, suggesting a possible connection between the enrichment of these Pseudomonas species and a plant-driven exudate profile.

Functional genetics of human gut commensal Bacteroides thetaiotaomicron reveals metabolic requirements for growth across environments

Harnessing the microbiota for beneficial outcomes is limited by our poor understanding of the constituent bacteria, as the functions of most of their genes are unknown. Here, we measure the growth of a barcoded transposon mutant library of the gut commensal Bacteroides thetaiotaomicron on 48 carbon sources, in the presence of 56 stress-inducing compounds, and during mono-colonization of gnotobiotic mice. We identify 516 genes with a specific phenotype under only one or a few conditions, enabling informed predictions of gene function. For example, we identify a glycoside hydrolase important for growth on type I rhamnogalacturonan, a DUF4861 protein for glycosaminoglycan utilization, a 3-keto-glucoside hydrolase for disaccharide utilization, and a tripartite multidrug resistance system specifically for bile salt tolerance. Furthermore, we show that B. thetaiotaomicron uses alternative enzymes for synthesizing nitrogen-containing metabolic precursors based on ammonium availability and that these enzymes are used differentially in vivo in a diet-dependent manner.

Four families of folate-independent methionine synthases

Although most organisms synthesize methionine from homocysteine and methyl folates, some have “core” methionine synthases that lack folate-binding domains and use other methyl donors. In vitro, the characterized core synthases use methylcobalamin as a methyl donor, but in vivo, they probably rely on corrinoid (vitamin B12-binding) proteins. We identified four families of core methionine synthases that are distantly related to each other (under 30% pairwise amino acid identity). From the characterized enzymes, we identified the families MesA, which is found in methanogens, and MesB, which is found in anaerobic bacteria and archaea with the Wood-Ljungdahl pathway. A third uncharacterized family, MesC, is found in anaerobic archaea that have the Wood-Ljungdahl pathway and lack known forms of methionine synthase. We predict that most members of the MesB and MesC families accept methyl groups from the iron-sulfur corrinoid protein of that pathway. The fourth family, MesD, is found only in aerobic bacteria. Using transposon mutants and complementation, we show that MesD does not require 5-methyltetrahydrofolate or cobalamin. Instead, MesD requires an uncharacterized protein family (DUF1852) and oxygen for activity.

Methionine is one of the amino acids that make up proteins, and the final step in methionine synthesis is the transfer of a methyl group. In most organisms, the methyl group is obtained from methyl folates, but some anaerobic bacteria and archaea are thought to use corrinoid (vitamin B12-binding) proteins instead. By analyzing the sequences of the potential methionine synthases across the genomes of diverse bacteria and archaea, we identified four families of folate-independent methionine synthases. For three of these families, we can use co-occurrence with corrinoid proteins to predict their likely partners. We show that the fourth family does not require vitamin B12; instead, it obtains methyl groups from an oxygen-dependent partner protein. Our results will help us understand the growth requirements of diverse bacteria and archaea.

Characterization of a Metal-Resistant Bacillus Strain With a High Molybdate Affinity ModA From Contaminated Sediments at the Oak Ridge Reservation

A nitrate- and metal-contaminated site at the Oak Ridge Reservation (ORR) was previously shown to contain the metal molybdenum (Mo) at picomolar concentrations. This potentially limits microbial nitrate reduction, as Mo is required by the enzyme nitrate reductase, which catalyzes the first step of nitrate removal. Enrichment for anaerobic nitrate-reducing microbes from contaminated sediment at the ORR yielded Bacillus strain EB106-08-02-XG196. This bacterium grows in the presence of multiple metals (Cd, Ni, Cu, Co, Mn, and U) but also exhibits better growth compared to control strains, including Pseudomonas fluorescens N2E2 isolated from a pristine ORR environment under low molybdate concentrations (<1 nM). Molybdate is taken up by the molybdate binding protein, ModA, of the molybdate ATP-binding cassette transporter. ModA of XG196 is phylogenetically distinct from those of other characterized ModA proteins. The genes encoding ModA from XG196, P. fluorescens N2E2 and Escherichia coli K12 were expressed in E. coli and the recombinant proteins were purified. Isothermal titration calorimetry analysis showed that XG196 ModA has a higher affinity for molybdate than other ModA proteins with a molybdate binding constant (K_D) of 2.2 nM, about one order of magnitude lower than those of P. fluorescens N2E2 (27.0 nM) and E. coli K12 (25.0 nM). XG196 ModA also showed a fivefold higher affinity for molybdate than for tungstate (11 nM), whereas the ModA proteins from P. fluorescens N2E2 [K_D (Mo) 27.0 nM, K_D (W) 26.7 nM] and E. coli K12[(K_D (Mo) 25.0 nM, K_D (W) 23.8 nM] had similar affinities for the two oxyanions. We propose that high molybdate affinity coupled with resistance to multiple metals gives strain XG196 a competitive advantage in Mo-limited environments contaminated with high concentrations of metals and nitrate, as found at ORR.

Fatty Acid and Alcohol Metabolism in Pseudomonas putida: Functional Analysis Using Random Barcode Transposon Sequencing

With its ability to catabolize a wide variety of carbon sources and a growing engineering toolkit, Pseudomonas putida KT2440 is emerging as an important chassis organism for metabolic engineering. Despite advances in our understanding of the organism, many gaps remain in our knowledge of the genetic basis of its metabolic capabilities. The gaps are particularly noticeable in our understanding of both fatty acid and alcohol catabolism, where many paralogs putatively coding for similar enzymes coexist, making biochemical assignment via sequence homology difficult. To rapidly assign function to the enzymes responsible for these metabolisms, we leveraged random barcode transposon sequencing (RB-Tn-Seq). Global fitness analyses of transposon libraries grown on 13 fatty acids and 10 alcohols produced strong phenotypes for hundreds of genes. Fitness data from mutant pools grown on fatty acids of varying chain lengths indicated specific enzyme substrate preferences and enabled us to hypothesize that DUF1302/DUF1329 family proteins potentially function as esterases. From the data, we also postulate catabolic routes for the two biogasoline molecules isoprenol and isopentanol, which are catabolized via leucine metabolism after initial oxidation and activation with coenzyme A (CoA). Because fatty acids and alcohols may serve as both feedstocks and final products of metabolic-engineering efforts, the fitness data presented here will help guide future genomic modifications toward higher titers, rates, and yields.IMPORTANCE To engineer novel metabolic pathways into P. putida, a comprehensive understanding of the genetic basis of its versatile metabolism is essential. Here, we provide functional evidence for the putative roles of hundreds of genes involved in the fatty acid and alcohol metabolism of the bacterium. These data provide a framework facilitating precise genetic changes to prevent product degradation and to channel the flux of specific pathway intermediates as desired.

High-throughput mapping of the phage resistance landscape in E. coli

Bacteriophages (phages) are critical players in the dynamics and function of microbial communities and drive processes as diverse as global biogeochemical cycles and human health. Phages tend to be predators finely tuned to attack specific hosts, even down to the strain level, which in turn defend themselves using an array of mechanisms. However, to date, efforts to rapidly and comprehensively identify bacterial host factors important in phage infection and resistance have yet to be fully realized. Here, we globally map the host genetic determinants involved in resistance to 14 phylogenetically diverse double-stranded DNA phages using two model Escherichia coli strains (K-12 and BL21) with known sequence divergence to demonstrate strain-specific differences. Using genome-wide loss-of-function and gain-of-function genetic technologies, we are able to confirm previously described phage receptors as well as uncover a number of previously unknown host factors that confer resistance to one or more of these phages. We uncover differences in resistance factors that strongly align with the susceptibility of K-12 and BL21 to specific phage. We also identify both phage-specific mechanisms, such as the unexpected role of cyclic-di-GMP in host sensitivity to phage N4, and more generic defenses, such as the overproduction of colanic acid capsular polysaccharide that defends against a wide array of phages. Our results indicate that host responses to phages can occur via diverse cellular mechanisms. Our systematic and high-throughput genetic workflow to characterize phage-host interaction determinants can be extended to diverse bacteria to generate datasets that allow predictive models of how phage-mediated selection will shape bacterial phenotype and evolution. The results of this study and future efforts to map the phage resistance landscape will lead to new insights into the coevolution of hosts and their phage, which can ultimately be used to design better phage therapeutic treatments and tools for precision microbiome engineering.

Distinctiveness of genes contributing to growth of Pseudomonas syringae in diverse host plant species

A variety of traits are necessary for bacterial colonization of the interior of plant hosts, including well-studied virulence effectors as well as other phenotypes contributing to bacterial growth and survival within the apoplast. High-throughput methods such as transposon sequencing (TnSeq) are powerful tools to identify such genes in bacterial pathogens. However, there is little information as to the distinctiveness of traits required for bacterial colonization of different hosts. Here, we utilize randomly barcoded TnSeq (RB-TnSeq) to identify the genes that contribute to the ability of Pseudomonas syringae strain B728a to grow within common bean (Phaseolus vulgaris), lima bean (Phaseolus lunatus), and pepper (Capsicum annuum); species representing two different plant families. The magnitude of contribution of most genes to apoplastic fitness in each of the plant hosts was similar. However, 50 genes significantly differed in their fitness contributions to growth within these species. These genes encoded proteins in various functional categories including polysaccharide synthesis and transport, amino acid metabolism and transport, cofactor metabolism, and phytotoxin synthesis and transport. Six genes that encoded unannotated, hypothetical proteins also contributed differentially to growth in these hosts. The genetic repertoire of a relatively promiscuous pathogen such as P. syringae may thus be shaped, at least in part, by the conditional contribution of some fitness determinants.

Selective carbon sources influence the end products of microbial nitrate respiration

Respiratory and catabolic genes are differentially distributed across microbial genomes. Thus, specific carbon sources may favor different respiratory processes. We profiled the influence of 94 carbon sources on the end products of nitrate respiration in microbial enrichment cultures from diverse terrestrial environments. We found that some carbon sources consistently favor dissimilatory nitrate reduction to ammonium (DNRA/nitrate ammonification) while other carbon sources favor nitrite accumulation or denitrification. For an enrichment culture from aquatic sediment, we sequenced the genomes of the most abundant strains, matched these genomes to 16S rDNA exact sequence variants (ESVs), and used 16S rDNA amplicon sequencing to track the differential enrichment of functionally distinct ESVs on different carbon sources. We found that changes in the abundances of strains with different genetic potentials for nitrite accumulation, DNRA or denitrification were correlated with the nitrite or ammonium concentrations in the enrichment cultures recovered on different carbon sources. Specifically, we found that either L-sorbose or D-cellobiose enriched for a Klebsiella nitrite accumulator, other sugars enriched for an Escherichia nitrate ammonifier, and citrate or formate enriched for a Pseudomonas denitrifier and a Sulfurospirillum nitrate ammonifier. Our results add important nuance to the current paradigm that higher concentrations of carbon will always favor DNRA over denitrification or nitrite accumulation, and we propose that, in some cases, carbon composition can be as important as carbon concentration in determining nitrate respiratory end products. Furthermore, our approach can be extended to other environments and metabolisms to characterize how selective parameters influence microbial community composition, gene content, and function.

Iron Supplementation Eliminates Antagonistic Interactions Between Root-Associated Bacteria

The rhizosphere microbiome (rhizobiome) plays a critical role in plant health and development. However, the processes by which the constituent microbes interact to form and maintain a community are not well understood. To investigate these molecular processes, we examined pairwise interactions between 11 different microbial isolates under select nutrient-rich and nutrient-limited conditions. We observed that when grown with media supplemented with 56 mM glucose, two microbial isolates were able to inhibit the growth of six other microbes. The interaction between microbes persisted even after the antagonistic microbe was removed, upon exposure to spent media. To probe the genetic basis for these antagonistic interactions, we used a barcoded transposon library in a proxy bacterium, Pseudomonas putida, to identify genes which showed enhanced sensitivity to the antagonistic factor(s) secreted by Acinetobacter sp. 02. Iron metabolism-related gene clusters in P. putida were implicated by this systems-level analysis. The supplementation of iron prevented the antagonistic interaction in the original microbial pair, supporting the hypothesis that iron limitation drives antagonistic microbial interactions between rhizobionts. We conclude that rhizobiome community composition is influenced by competition for limiting nutrients, with implications for growth and development of the plant.

GapMind: Automated Annotation of Amino Acid Biosynthesis

GapMind is a Web-based tool for annotating amino acid biosynthesis in bacteria and archaea (http://papers.genomics.lbl.gov/gaps). GapMind incorporates many variant pathways and 130 different reactions, and it analyzes a genome in just 15 s. To avoid error-prone transitive annotations, GapMind relies primarily on a database of experimentally characterized proteins. GapMind correctly handles fusion proteins and split proteins, which often cause errors for best-hit approaches. To improve GapMind's coverage, we examined genetic data from 35 bacteria that grow in defined media without amino acids, and we filled many gaps in amino acid biosynthesis pathways. For example, we identified additional genes for arginine synthesis with succinylated intermediates in Bacteroides thetaiotaomicron, and we propose that Dyella japonica synthesizes tyrosine from phenylalanine. Nevertheless, for many bacteria and archaea that grow in minimal media, genes for some steps still cannot be identified. To help interpret potential gaps, GapMind checks if they match known gaps in related microbes that can grow in minimal media. GapMind should aid the identification of microbial growth requirements.IMPORTANCE Many microbes can make all of the amino acids (the building blocks of proteins). In principle, we should be able to predict which amino acids a microbe can make, and which it requires as nutrients, by checking its genome sequence for all of the necessary genes. However, in practice, it is difficult to check for all of the alternative pathways. Furthermore, new pathways and enzymes are still being discovered. We built an automated tool, GapMind, to annotate amino acid biosynthesis in bacterial and archaeal genomes. We used GapMind to list gaps: cases where a microbe makes an amino acid but a complete pathway cannot be identified in its genome. We used these gaps, together with data from mutants, to identify new pathways and enzymes. However, for most bacteria and archaea, we still do not know how they can make all of the amino acids.

Leveraging host metabolism for bisdemethoxycurcumin production in Pseudomonas putida

Pseudomonas putida is a saprophytic bacterium with robust metabolisms and strong solvent tolerance making it an attractive host for metabolic engineering and bioremediation. Due to its diverse carbon metabolisms, its genome encodes an array of proteins and enzymes that can be readily applied to produce valuable products. In this work we sought to identify design principles and bottlenecks in the production of type III polyketide synthase (T3PKS)-derived compounds in P. putida. T3PKS products are widely used as nutraceuticals and medicines and often require aromatic starter units, such as coumaroyl-CoA, which is also an intermediate in the native coumarate catabolic pathway of P. putida. Using a randomly barcoded transposon mutant (RB-TnSeq) library, we assayed gene functions for a large portion of aromatic catabolism, confirmed known pathways, and proposed new annotations for two aromatic transporters. The 1,3,6,8-tetrahydroxynapthalene synthase of Streptomyces coelicolor (RppA), a microbial T3PKS, was then used to rapidly assay growth conditions for increased T3PKS product accumulation. The feruloyl/coumaroyl CoA synthetase (Fcs) of P. putida was used to supply coumaroyl-CoA for the curcuminoid synthase (CUS) of Oryza sativa, a plant T3PKS. We identified that accumulation of coumaroyl-CoA in this pathway results in extended growth lag times in P. putida. Deletion of the second step in coumarate catabolism, the enoyl-CoA hydratase-lyase (Ech), resulted in increased production of the type III polyketide bisdemethoxycurcumin.

Genome-Wide Transposon Screen of a Pseudomonas syringae mexB Mutant Reveals the Substrates of Efflux Transporters

Bacteria express numerous efflux transporters that confer resistance to diverse toxicants present in their environment. Due to a high level of functional redundancy of these transporters, it is difficult to identify those that are of most importance in conferring resistance to specific compounds. The resistance-nodulation-division (RND) protein family is one such example of redundant transporters that are widespread among Gram-negative bacteria. Within this family, the MexAB-OprM protein complex is highly expressed and conserved among Pseudomonas species. We exposed barcoded transposon mutant libraries in isogenic wild-type and ΔmexB backgrounds in P. syringae B728a to diverse toxic compounds in vitro to identify mutants with increased susceptibility to these compounds. Mutants with mutations in genes encoding both known and novel redundant transporters but with partially overlapping substrate specificities were observed in a ΔmexB background. Psyr_0228, an uncharacterized member of the major facilitator superfamily of transporters, preferentially contributes to tolerance of acridine orange and acriflavine. Another transporter located in the inner membrane, Psyr_0541, contributes to tolerance of acriflavine and berberine. The presence of multiple redundant, genomically encoded efflux transporters appears to enable bacterial strains to tolerate a diversity of environmental toxins. This genome-wide screen performed in a hypersusceptible mutant strain revealed numerous transporters that would otherwise be dispensable under these conditions. Bacterial strains such as P. syringae that likely encounter diverse toxins in their environment, such as in association with many different plant species, probably benefit from possessing multiple redundant transporters that enable versatility with respect to toleration of novel toxicants.IMPORTANCE Bacteria use protein pumps to remove toxic compounds from the cell interior, enabling survival in diverse environments. These protein pumps can be highly redundant, making their targeted examination difficult. In this study, we exposed mutant populations of Pseudomonas syringae to diverse toxicants to identify pumps that contributed to survival in those conditions. In parallel, we examined pump redundancy by testing mutants of a population lacking the primary efflux transporter responsible for toxin tolerance. We identified partial substrate overlap for redundant transporters, as well as several pumps that appeared more substrate specific. For bacteria that are found in diverse environments, having multiple, partially redundant efflux pumps likely allows flexibility in habitat colonization.

Omics-driven identification and elimination of valerolactam catabolism in Pseudomonas putida KT2440 for increased product titer

Pseudomonas putida is a promising bacterial chassis for metabolic engineering given its ability to metabolize a wide array of carbon sources, especially aromatic compounds derived from lignin. However, this omnivorous metabolism can also be a hindrance when it can naturally metabolize products produced from engineered pathways. Herein we show that P. putida is able to use valerolactam as a sole carbon source, as well as degrade caprolactam. Lactams represent important nylon precursors, and are produced in quantities exceeding one million tons per year (Zhang et al., 2017). To better understand this metabolism we use a combination of Random Barcode Transposon Sequencing (RB-TnSeq) and shotgun proteomics to identify the oplBA locus as the likely responsible amide hydrolase that initiates valerolactam catabolism. Deletion of the oplBA genes prevented P. putida from growing on valerolactam, prevented the degradation of valerolactam in rich media, and dramatically reduced caprolactam degradation under the same conditions. Deletion of oplBA, as well as pathways that compete for precursors L-lysine or 5-aminovalerate, increased the titer of valerolactam from undetectable after 48 h of production to ~90 mg/L. This work may serve as a template to rapidly eliminate undesirable metabolism in non-model hosts in future metabolic engineering efforts.

•

P. putida utilizes valerolactam as a sole carbon source and degrades caprolactam.

•

With RB-TnSeq and proteomics we identify the lactam hydrolytic enzyme OplBA.

•

Deleting oplBA prevents P. putida from growing on valerolactam.

•

We increased titers of valerolactam in P. putida from 0 mg/L to ~90 mg/L.

The selective pressures on the microbial community in a metal-contaminated aquifer

In many environments, toxic compounds restrict which microorganisms persist. However, in complex mixtures of inhibitory compounds, it is challenging to determine which specific compounds cause changes in abundance and prevent some microorganisms from growing. We focused on a contaminated aquifer in Oak Ridge, Tennessee, USA that has large gradients of pH and widely varying concentrations of uranium, nitrate, and many other inorganic ions. In the most contaminated wells, the microbial community is enriched in the Rhodanobacter genus. Rhodanobacter abundance is positively correlated with low pH and high concentrations of uranium and 13 other ions and we sought to determine which of these ions are selective pressures that favor the growth of Rhodanobacter over other taxa. Of these ions, low pH and high UO₂²⁺, Mn²⁺, Al³⁺, Cd²⁺, Zn²⁺, Co²⁺, and Ni²⁺ are both (a) selectively inhibitory of a Pseudomonas isolate from an uncontaminated well vs. a Rhodanobacter isolate from a contaminated well, and (b) reach toxic concentrations (for the Pseudomonas isolate) in the Rhodanobacter-dominated wells. We used mixtures of ions to simulate the groundwater conditions in the most contaminated wells and verified that few isolates aside from Rhodanobacter can tolerate these eight ions. These results clarify which ions are likely causal factors that impact the microbial community at this field site and are not merely correlated with taxonomic shifts. Furthermore, our general high-throughput approach can be applied to other environments, isolates, and conditions to systematically help identify selective pressures on microbial communities.

Principles of Systems Biology, No. 30

This month: two examples of door-opening, innovative microscopy (Garcia and also Benzinger et al.), expanding our functional knowledge of bacteria by over 10,000 genes (Deutschbauer), and probing how RNA structure dictates inclusion in liquid-like droplets in vivo (Langdon and Gladfelter).

Mutant phenotypes for thousands of bacterial genes of unknown function

One-third of all protein-coding genes from bacterial genomes cannot be annotated with a function. Here, to investigate the functions of these genes, we present genome-wide mutant fitness data from 32 diverse bacteria across dozens of growth conditions. We identified mutant phenotypes for 11,779 protein-coding genes that had not been annotated with a specific function. Many genes could be associated with a specific condition because the gene affected fitness only in that condition, or with another gene in the same bacterium because they had similar mutant phenotypes. Of the poorly annotated genes, 2,316 had associations that have high confidence because they are conserved in other bacteria. By combining these conserved associations with comparative genomics, we identified putative DNA repair proteins; in addition, we propose specific functions for poorly annotated enzymes and transporters and for uncharacterized protein families. Our study demonstrates the scalability of microbial genetics and its utility for improving gene annotations.

Filling gaps in bacterial amino acid biosynthesis pathways with high-throughput genetics

For many bacteria with sequenced genomes, we do not understand how they synthesize some amino acids. This makes it challenging to reconstruct their metabolism, and has led to speculation that bacteria might be cross-feeding amino acids. We studied heterotrophic bacteria from 10 different genera that grow without added amino acids even though an automated tool predicts that the bacteria have gaps in their amino acid synthesis pathways. Across these bacteria, there were 11 gaps in their amino acid biosynthesis pathways that we could not fill using current knowledge. Using genome-wide mutant fitness data, we identified novel enzymes that fill 9 of the 11 gaps and hence explain the biosynthesis of methionine, threonine, serine, or histidine by bacteria from six genera. We also found that the sulfate-reducing bacterium Desulfovibrio vulgaris synthesizes homocysteine (which is a precursor to methionine) by using DUF39, NIL/ferredoxin, and COG2122 proteins, and that homoserine is not an intermediate in this pathway. Our results suggest that most free-living bacteria can likely make all 20 amino acids and illustrate how high-throughput genetics can uncover previously-unknown amino acid biosynthesis genes.

For a few bacteria, it is well known how they can make all 20 of the standard amino acids (the building blocks of proteins). For many other bacteria, their genome sequence implies that there are gaps in these biosynthetic pathways, so that the bacteria cannot make all of the amino acids and would need to take up some of them from their environment instead. But many bacteria can grow in minimal media (without any amino acids) despite these apparent gaps. We studied 10 bacteria with predicted gaps in amino acid biosynthesis that nevertheless grow in minimal media. Most of these gaps were spurious, but 11 of the gaps were genuine and could not be explained by current knowledge. Using high-throughput genetics, we systematically identified genes that were required for growth in minimal media and identified the biosynthetic genes that fill 9 of the 11 gaps. We hope that this approach can be applied to many more bacteria and will eventually allow us to accurately predict the nutritional requirements of a bacterium from its genome sequence.

A metabolic pathway for catabolizing levulinic acid in bacteria

Microorganisms can catabolize a wide range of organic compounds and therefore have the potential to perform many industrially relevant bioconversions. One barrier to realizing the potential of biorefining strategies lies in our incomplete knowledge of metabolic pathways, including those that can be used to assimilate naturally abundant or easily generated feedstocks. For instance, levulinic acid (LA) is a carbon source that is readily obtainable as a dehydration product of lignocellulosic biomass and can serve as the sole carbon source for some bacteria. Yet, the genetics and structure of LA catabolism have remained unknown. Here, we report the identification and characterization of a seven-gene operon that enables LA catabolism in Pseudomonas putida KT2440. When the pathway was reconstituted with purified proteins, we observed the formation of four acyl-CoA intermediates, including a unique 4-phosphovaleryl-CoA and the previously observed 3-hydroxyvaleryl-CoA product. Using adaptive evolution, we obtained a mutant of Escherichia coli LS5218 with functional deletions of fadE and atoC that was capable of robust growth on LA when it expressed the five enzymes from the P. putida operon. This discovery will enable more efficient use of biomass hydrolysates and metabolic engineering to develop bioconversions using LA as a feedstock.

Genome-wide identification of bacterial plant colonization genes

Diverse soil-resident bacteria can contribute to plant growth and health, but the molecular mechanisms enabling them to effectively colonize their plant hosts remain poorly understood. We used randomly barcoded transposon mutagenesis sequencing (RB-TnSeq) in Pseudomonas simiae, a model root-colonizing bacterium, to establish a genome-wide map of bacterial genes required for colonization of the Arabidopsis thaliana root system. We identified 115 genes (2% of all P. simiae genes) with functions that are required for maximal competitive colonization of the root system. Among the genes we identified were some with obvious colonization-related roles in motility and carbon metabolism, as well as 44 other genes that had no or vague functional predictions. Independent validation assays of individual genes confirmed colonization functions for 20 of 22 (91%) cases tested. To further characterize genes identified by our screen, we compared the functional contributions of P. simiae genes to growth in 90 distinct in vitro conditions by RB-TnSeq, highlighting specific metabolic functions associated with root colonization genes. Our analysis of bacterial genes by sequence-driven saturation mutagenesis revealed a genome-wide map of the genetic determinants of plant root colonization and offers a starting point for targeted improvement of the colonization capabilities of plant-beneficial microbes.

Plants fix carbon to create an abundance of sugars and amino acids, thus providing an enticing environment for microorganisms that reside in soil. Once these microorganisms have colonized the root environment, they can dramatically influence plant growth and development. We set out to identify a comprehensive set of microbial genes that control or influence root colonization, using a genome-wide transposon mutagenesis approach (randomly barcoded transposon sequencing [RB-TnSeq]). By using this method, we identified several hundred genes that, when mutated, affect the ability of the bacterium P. simiae to competitively colonize the root system of the model plant A. thaliana. These included many genes purported to be involved in carbohydrate metabolism, cell wall biosynthesis, and motility, underscoring the notion that sugar metabolism, defense, and motility are all key features of a root-colonizing microbe. We also identified several amino acid transport and metabolism genes with mutations that confer a fitness advantage in root colonization. Lastly, we identified several genes with no known function that significantly alter root colonization ability when mutated. These findings suggest novel engineering strategies to improve biological product development, and will facilitate the mechanistic exploration of the root colonization process.

Mechanisms of Chromium and Uranium Toxicity in Pseudomonas stutzeri RCH2 Grown under Anaerobic Nitrate-Reducing Conditions

Chromium and uranium are highly toxic metals that contaminate many natural environments. We investigated their mechanisms of toxicity under anaerobic conditions using nitrate-reducing Pseudomonas stutzeri RCH2, which was originally isolated from a chromium-contaminated aquifer. A random barcode transposon site sequencing library of RCH2 was grown in the presence of the chromate oxyanion (Cr[VI]O42−) or uranyl oxycation (U[VI]O22+). Strains lacking genes required for a functional nitrate reductase had decreased fitness as both metals interacted with heme-containing enzymes required for the later steps in the denitrification pathway after nitrate is reduced to nitrite. Cr[VI]-resistance also required genes in the homologous recombination and nucleotide excision DNA repair pathways, showing that DNA is a target of Cr[VI] even under anaerobic conditions. The reduced thiol pool was also identified as a target of Cr[VI] toxicity and psest_2088, a gene of previously unknown function, was shown to have a role in the reduction of sulfite to sulfide. U[VI] resistance mechanisms involved exopolysaccharide synthesis and the universal stress protein UspA. As the first genome-wide fitness analysis of Cr[VI] and U[VI] toxicity under anaerobic conditions, this study provides new insight into the impact of Cr[VI] and U[VI] on an environmental isolate from a chromium contaminated site, as well as into the role of a ubiquitous protein, Psest_2088.

Draft Genome Sequences of Two Janthinobacteriumlividum Strains, Isolated from Pristine Groundwater Collected from the Oak Ridge Field Research Center

We present here the draft genome sequences of two Janthinobacterium lividum strains, GW456P and GW458P, isolated from groundwater samples collected from a background site at the Oak Ridge Field Research Center. Production of a purple pigment by these two strains was observed when grown on diluted (1/10) LB agar plates.

System-Wide Adaptations of Desulfovibrio alaskensis G20 to Phosphate-Limited Conditions

The prevalence of lipids devoid of phosphorus suggests that the availability of phosphorus limits microbial growth and activity in many anoxic, stratified environments. To better understand the response of anaerobic bacteria to phosphate limitation and starvation, this study combines microscopic and lipid analyses with the measurements of fitness of pooled barcoded transposon mutants of the model sulfate reducing bacterium Desulfovibrio alaskensis G20. Phosphate-limited G20 has lower growth rates and replaces more than 90% of its membrane phospholipids by a mixture of monoglycosyl diacylglycerol (MGDG), glycuronic acid diacylglycerol (GADG) and ornithine lipids, lacks polyphosphate granules, and synthesizes other cellular inclusions. Analyses of pooled and individual mutants reveal the importance of the high-affinity phosphate transport system (the Pst system), PhoR, and glycolipid and ornithine lipid synthases during phosphate limitation. The phosphate-dependent synthesis of MGDG in G20 and the widespread occurrence of the MGDG/GADG synthase among sulfate reducing ∂-Proteobacteria implicate these microbes in the production of abundant MGDG in anaerobic environments where the concentrations of phosphate are lower than 10 μM. Numerous predicted changes in the composition of the cell envelope and systems involved in transport, maintenance of cytoplasmic redox potential, central metabolism and regulatory pathways also suggest an impact of phosphate limitation on the susceptibility of sulfate reducing bacteria to other anthropogenic or environmental stresses.

Exometabolomics Assisted Design and Validation of Synthetic Obligate Mutualism

Synthetic microbial ecology has the potential to enhance the productivity and resiliency of biotechnology processes compared to approaches using single isolates. Engineering microbial consortia is challenging; however, one approach that has attracted significant attention is the creation of synthetic obligate mutualism using auxotrophic mutants that depend on each other for exchange or cross-feeding of metabolites. Here, we describe the integration of mutant library fitness profiling with mass spectrometry based exometabolomics as a method for constructing synthetic mutualism based on cross-feeding. Two industrially important species lacking known ecological interactions, Zymomonas mobilis and Escherichia coli, were selected as the test species. Amino acid exometabolites identified in the spent medium of Z. mobilis were used to select three corresponding E. coli auxotrophs (proA, pheA and IlvA), as potential E. coli counterparts for the coculture. A pooled mutant fitness assay with a Z. mobilis transposon mutant library was used to identify mutants with improved growth in the presence of E. coli. An auxotroph mutant in a gene (ZMO0748) with sequence similarity to cysteine synthase A (cysK), was selected as the Z. mobilis counterpart for the coculture. Exometabolomic analysis of spent E. coli medium identified glutathione related metabolites as potentially available for rescue of the Z. mobilis cysteine synthase mutant. Three sets of cocultures between the Z. mobilis auxotroph and each of the three E. coli auxotrophs were monitored by optical density for growth and analyzed by flow cytometry to confirm high cell counts for each species. Taken together, our methods provide a technological framework for creating synthetic mutualisms combining existing screening based methods and exometabolomics for both the selection of obligate mutualism partners and elucidation of metabolites involved in auxotroph rescue.

The Epigenomic Landscape of Prokaryotes

DNA methylation acts in concert with restriction enzymes to protect the integrity of prokaryotic genomes. Studies in a limited number of organisms suggest that methylation also contributes to prokaryotic genome regulation, but the prevalence and properties of such non-restriction-associated methylation systems remain poorly understood. Here, we used single molecule, real-time sequencing to map DNA modifications including m6A, m4C, and m5C across the genomes of 230 diverse bacterial and archaeal species. We observed DNA methylation in nearly all (93%) organisms examined, and identified a total of 834 distinct reproducibly methylated motifs. This data enabled annotation of the DNA binding specificities of 620 DNA Methyltransferases (MTases), doubling known specificities for previously hard to study Type I, IIG and III MTases, and revealing their extraordinary diversity. Strikingly, 48% of organisms harbor active Type II MTases with no apparent cognate restriction enzyme. These active ‘orphan’ MTases are present in diverse bacterial and archaeal phyla and show motif specificities and methylation patterns consistent with functions in gene regulation and DNA replication. Our results reveal the pervasive presence of DNA methylation throughout the prokaryotic kingdoms, as well as the diversity of sequence specificities and potential functions of DNA methylation systems.

DNA methylation is a chemical modification of DNA present in many prokaryotic genomes. The best-known role of DNA methylation is as a component of restriction-modification systems. In these systems, restriction enzymes target foreign DNA for cleavage, while DNA methylation protects the host genome from destruction. Studies in a handful of organisms show that DNA methylation may also act independently of restriction systems and function in genome regulation. However, a lack of technologies has limited the study of DNA methylation to a small number of organisms, and the broader patterns and functions of DNA methylation remain unknown. Here we use SMRT-sequencing to determine the genome wide DNA methylation patterns of more than 200 diverse bacteria and archaea. We show that DNA methylation is pervasive and present in more than 90% of studied organisms. Analysis of this data enabled annotation of the specific DNA binding sites of more than 600 restriction systems, revealing their extraordinary diversity. Strikingly, we observed widespread DNA methylation in the absence of restriction systems. Analyses of these patterns reveal that they are conserved through evolution, and likely function in genome regulation. Thus DNA methylation may play a far wider function in prokaryotic genome biology than was previously supposed.

Monofluorophosphate is a selective inhibitor of respiratory sulfate-reducing microorganisms

Despite the environmental and economic cost of microbial sulfidogenesis in industrial operations, few compounds are known as selective inhibitors of respiratory sulfate reducing microorganisms (SRM), and no study has systematically and quantitatively evaluated the selectivity and potency of SRM inhibitors. Using general, high-throughput assays to quantitatively evaluate inhibitor potency and selectivity in a model sulfate-reducing microbial ecosystem as well as inhibitor specificity for the sulfate reduction pathway in a model SRM, we screened a panel of inorganic oxyanions. We identified several SRM selective inhibitors including selenate, selenite, tellurate, tellurite, nitrate, nitrite, perchlorate, chlorate, monofluorophosphate, vanadate, molydate, and tungstate. Monofluorophosphate (MFP) was not known previously as a selective SRM inhibitor, but has promising characteristics including low toxicity to eukaryotic organisms, high stability at circumneutral pH, utility as an abiotic corrosion inhibitor, and low cost. MFP remains a potent inhibitor of SRM growing by fermentation, and MFP is tolerated by nitrate and perchlorate reducing microorganisms. For SRM inhibition, MFP is synergistic with nitrite and chlorite, and could enhance the efficacy of nitrate or perchlorate treatments. Finally, MFP inhibition is multifaceted. Both inhibition of the central sulfate reduction pathway and release of cytoplasmic fluoride ion are implicated in the mechanism of MFP toxicity.