Design, construction, and in vivo augmentation of a complex gut microbiome.

Efforts to model the human gut microbiome in mice have led to important insights into the mechanisms of host-microbe interactions. However, the model communities studied to date have been defined or complex, but not both, limiting their utility. Here, we construct and characterize in vitro a defined community of 104 bacterial species composed of the most common taxa from the human gut microbiota (hCom1). We then used an iterative experimental process to fill open niches: germ-free mice were colonized with hCom1 and then challenged with a human fecal sample. We identified new species that engrafted following fecal challenge and added them to hCom1, yielding hCom2. In gnotobiotic mice, hCom2 exhibited increased stability to fecal challenge and robust colonization resistance against pathogenic Escherichia coli. Mice colonized by either hCom2 or a human fecal community are phenotypically similar, suggesting that this consortium will enable a mechanistic interrogation of species and genes on microbiome-associated phenotypes.

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.

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.

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.

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.

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

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

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.

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.

4'-Deoxypyridoxine disrupts vitamin B6 homeostasis in Escherichia coli K12 through combined inhibition of cumulative B6 uptake and PLP-dependent enzyme activity.

Pyridoxal 5'-phosphate (PLP) is the active form of vitamin B6 and a cofactor for many essential metabolic processes such as amino acid biosynthesis and one carbon metabolism. 4'-deoxypyridoxine (4dPN) is a long known B6 antimetabolite but its mechanism of action was not totally clear. By exploring different conditions in which PLP metabolism is affected in the model organism Escherichia coli K12, we showed that 4dPN cannot be used as a source of vitamin B6 as previously claimed and that it is toxic in several conditions where vitamin B6 homeostasis is affected, such as in a B6 auxotroph or in a mutant lacking the recently discovered PLP homeostasis gene, yggS. In addition, we found that 4dPN sensitivity is likely the result of multiple modes of toxicity, including inhibition of PLP-dependent enzyme activity by 4'-deoxypyridoxine phosphate (4dPNP) and inhibition of cumulative pyridoxine (PN) uptake. These toxicities are largely dependent on the phosphorylation of 4dPN by pyridoxal kinase (PdxK).

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.

Construction and characterization of a genome-scale ordered mutant collection of Bacteroides thetaiotaomicron.

BACKGROUND: Ordered transposon-insertion collections, in which specific transposon-insertion mutants are stored as monocultures in a genome-scale collection, represent a promising tool for genetic dissection of human gut microbiota members. However, publicly available collections are scarce and the construction methodology remains in early stages of development. RESULTS: Here, we describe the assembly of a genome-scale ordered collection of transposon-insertion mutants in the model gut anaerobe Bacteroides thetaiotaomicron VPI-5482 that we created as a resource for the research community. We used flow cytometry to sort single cells from a pooled library, located mutants within this initial progenitor collection by applying a pooling strategy with barcode sequencing, and re-arrayed specific mutants to create a condensed collection with single-insertion strains covering >2500 genes. To demonstrate the potential of the condensed collection for phenotypic screening, we analyzed growth dynamics and cell morphology. We identified both growth defects and altered cell shape in mutants disrupting sphingolipid synthesis and thiamine scavenging. Finally, we analyzed the process of assembling the B. theta condensed collection to identify inefficiencies that limited coverage. We demonstrate as part of this analysis that the process of assembling an ordered collection can be accurately modeled using barcode sequencing data. CONCLUSION: We expect that utilization of this ordered collection will accelerate research into B. theta physiology and that lessons learned while assembling the collection will inform future efforts to assemble ordered mutant collections for an increasing number of gut microbiota members.

Study of the Production of Poly(Hydroxybutyrate-co-Hydroxyhexanoate) and Poly(Hydroxybutyrate-co-Hydroxyvalerate-co-Hydroxyhexanoate) in Rhodospirillum rubrum.

Poly(hydroxybutyrate-co-hydroxyhexanoate) [P(HB-co-HHx)] and poly(hydroxybutyrate-co-hydroxyvalerate-co-hydroxyhexanoate) [P(HB-co-HV-co-HHx)] demonstrate interesting mechanical and thermal properties as well as excellent biocompatibility, making them suitable for multiple applications and notably biomedical purposes. The production of such polymers was described in Rhodospirillum rubrum, a purple nonsulfur bacteria in a nutrient-lacking environment where the HHx synthesis is triggered by the presence of hexanoate in the medium. However, the production of P(HB-co-HHx) under nutrient-balanced growth conditions in R. rubrum has not been described so far, and the assimilation of hexanoate is poorly documented. In this study, we used proteomic analysis and a mutant fitness assay to demonstrate that hexanoate assimilation involve ß-oxidation and the ethylmalonyl-coenzyme A (CoA) (EMC) and methylbutanoyl-CoA (MBC) pathways, both being anaplerotic pathways already described in R. rubrum. Polyhydroxyalkanoate (PHA) production is likely to involve the de novo fatty acid synthesis pathway. Concerning the polymer composition, HB is the main component of the polymer, probably as acetyl-CoA and butyryl-CoA are intermediates of hexanoate assimilation pathways. When no essential nutrient is lacking in the medium, the synthesis of PHA seems to help maintain the redox balance of the cell. In this framework, we showed that the fixation of CO2 is required to sustain the growth. An increase in the proportion of HHx in the polymer was observed when redox stress was engendered in the cell under bicarbonate-limiting growth conditions. The addition of isoleucine or valerate in the medium also increased the HHx content of the polymer and allowed the production of a terpolymer of P(HB-co-HV-co-HHx). IMPORTANCE The use of purple bacteria, which can assimilate volatile fatty acids, for biotechnological applications has increased, since they reduce the production costs of added-value compounds such as PHA. P(HB-co-HHx) and P(HB-co-HV-co-HHx) have demonstrated interesting properties, notably for biomedical applications. In a nutrient-lacking environment, R. rubrum is known to synthesize such polymers when hexanoate is used as the carbon source. However, their production in R. rubrum in non-nutrient-lacking growth conditions has not been described so far, and the assimilation of hexanoate is poorly documented. As the carbon source and its assimilation directly impact the polymer composition, we studied under non-nutrient-lacking growth conditions the assimilation pathway of hexanoate and PHA production in R. rubrum. Proteomic analysis and mutant fitness assays allowed us to explain PHA production and composition. An increase in HHx content of the polymer and production of P(HB-co-HV-co-HHx) was possible using the knowledge gained on metabolism under hexanoate growth conditions.

Development of a Markerless Deletion Mutagenesis System in Nitrate-Reducing Bacterium Rhodanobacter denitrificans.

Rhodanobacter has been found as the dominant genus in aquifers contaminated with high concentrations of nitrate and uranium in Oak Ridge, TN, USA. The in situ stimulation of denitrification has been proposed as a potential method to remediate nitrate and uranium contamination. Among the Rhodanobacter species, Rhodanobacter denitrificans strains have been reported to be capable of denitrification and contain abundant metal resistance genes. However, due to the lack of a mutagenesis system in these strains, our understanding of the mechanisms underlying low-pH resistance and the ability to dominate in the contaminated environment remains limited. Here, we developed an in-frame markerless deletion system in two R. denitrificans strains. First, we optimized the growth conditions, tested antibiotic resistance, and determined appropriate transformation parameters in 10 Rhodanobacter strains. We then deleted the upp gene, which encodes uracil phosphoribosyltransferase, in R. denitrificans strains FW104-R3 and FW104-R5. The resulting strains were designated R3_Δupp and R5_Δupp and used as host strains for mutagenesis with 5-fluorouracil (5-FU) resistance as the counterselection marker to generate markerless deletion mutants. To test the developed protocol, the narG gene encoding nitrate reductase was knocked out in the R3_Δupp and R5_Δupp host strains. As expected, the narG mutants could not grow in anoxic medium with nitrate as the electron acceptor. Overall, these results show that the in-frame markerless deletion system is effective in two R. denitrificans strains, which will allow for future functional genomic studies in these strains furthering our understanding of the metabolic and resistance mechanisms present in Rhodanobacter species. IMPORTANCE Rhodanobacter denitrificans is capable of denitrification and is also resistant to toxic heavy metals and low pH. Accordingly, the presence of Rhodanobacter species at a particular environmental site is considered an indicator of nitrate and uranium contamination. These characteristics suggest its future potential application in bioremediation of nitrate or concurrent nitrate and uranium contamination in groundwater ecosystems. Due to the lack of genetic tools in this organism, the mechanisms of low-pH and heavy metal resistance in R. denitrificans strains remain elusive, which impedes its use in bioremediation strategies. Here, we developed a genome editing method in two R. denitrificans strains. This work marks a crucial step in developing Rhodanobacter as a model for studying the diverse mechanisms of low-pH and heavy metal resistance associated with denitrification.

Nitrogen Metabolism in Pseudomonas putida: Functional Analysis Using Random Barcode Transposon Sequencing.

Pseudomonas putida KT2440 has long been studied for its diverse and robust metabolisms, yet many genes and proteins imparting these growth capacities remain uncharacterized. Using pooled mutant fitness assays, we identified genes and proteins involved in the assimilation of 52 different nitrogen containing compounds. To assay amino acid biosynthesis, 19 amino acid drop-out conditions were also tested. From these 71 conditions, significant fitness phenotypes were elicited in 672 different genes including 100 transcriptional regulators and 112 transport-related proteins. We divide these conditions into 6 classes, and propose assimilatory pathways for the compounds based on this wealth of genetic data. To complement these data, we characterize the substrate range of three promiscuous aminotransferases relevant to metabolic engineering efforts in vitro. Furthermore, we examine the specificity of five transcriptional regulators, explaining some fitness data results and exploring their potential to be developed into useful synthetic biology tools. In addition, we use manifold learning to create an interactive visualization tool for interpreting our BarSeq data, which will improve the accessibility and utility of this work to other researchers. IMPORTANCE Understanding the genetic basis of P. putida's diverse metabolism is imperative for us to reach its full potential as a host for metabolic engineering. Many target molecules of the bioeconomy and their precursors contain nitrogen. This study provides functional evidence linking hundreds of genes to their roles in the metabolism of nitrogenous compounds, and provides an interactive tool for visualizing these data. We further characterize several aminotransferases, lactamases, and regulators, which are of particular interest for metabolic engineering.

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.

Genome-wide identification of Pseudomonas syringae genes required for fitness during colonization of the leaf surface and apoplast.

The foliar plant pathogen Pseudomonas syringae can establish large epiphytic populations on leaf surfaces before apoplastic colonization. However, the bacterial genes that contribute to these lifestyles have not been completely defined. The fitness contributions of 4,296 genes in P. syringae pv. syringae B728a were determined by genome-wide fitness profiling with a randomly barcoded transposon mutant library that was grown on the leaf surface and in the apoplast of the susceptible plant Phaseolus vulgaris Genes within the functional categories of amino acid and polysaccharide (including alginate) biosynthesis contributed most to fitness both on the leaf surface (epiphytic) and in the leaf interior (apoplast), while genes involved in type III secretion system and syringomycin synthesis were primarily important in the apoplast. Numerous other genes that had not been previously associated with in planta growth were also required for maximum epiphytic or apoplastic fitness. Fourteen hypothetical proteins and uncategorized glycosyltransferases were also required for maximum competitive fitness in and on leaves. For most genes, no relationship was seen between fitness in planta and either the magnitude of their expression in planta or degree of induction in planta compared to in vitro conditions measured in other studies. A lack of association of gene expression and fitness has important implications for the interpretation of transcriptional information and our broad understanding of plant-microbe interactions.

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

[This corrects the article DOI: 10.1371/journal.pgen.1007147.].

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.

Engineering Pseudomonas putida for efficient aromatic conversion to bioproduct using high throughput screening in a bioreactor.

Pseudomonas putida KT2440 is an emerging biomanufacturing host amenable for use with renewable carbon streams including aromatics such as para-coumarate. We used a pooled transposon library disrupting nearly all (4,778) non-essential genes to characterize this microbe under common stirred-tank bioreactor parameters with quantitative fitness assays. Assessing differential fitness values by monitoring changes in mutant strain abundance identified 33 gene mutants with improved fitness across multiple stirred-tank bioreactor formats. Twenty-one deletion strains from this subset were reconstructed, including GacA, a regulator, TtgB, an ABC transporter, and PP_0063, a lipid A acyltransferase. Thirteen deletion strains with roles in varying cellular functions were evaluated for conversion of para-coumarate, to a heterologous bioproduct, indigoidine. Several mutants, such as the ΔgacA strain improved fitness in a bioreactor by 35 fold and showed an 8-fold improvement in indigoidine production (4.5 g/L, 0.29 g/g, 23% of maximum theoretical yield) from para-coumarate as the carbon source.

Deciphering Microbial Metal Toxicity Responses via Random Bar Code Transposon Site Sequencing and Activity-Based Metabolomics.

To uncover metal toxicity targets and defense mechanisms of the facultative anaerobe Pantoea sp. strain MT58 (MT58), we used a multiomic strategy combining two global techniques, random bar code transposon site sequencing (RB-TnSeq) and activity-based metabolomics. MT58 is a metal-tolerant Oak Ridge Reservation (ORR) environmental isolate that was enriched in the presence of metals at concentrations measured in contaminated groundwater at an ORR nuclear waste site. The effects of three chemically different metals found at elevated concentrations in the ORR contaminated environment were investigated: the cation Al3+, the oxyanion CrO42-, and the oxycation UO22+. Both global techniques were applied using all three metals under both aerobic and anaerobic conditions to elucidate metal interactions mediated through the activity of metabolites and key genes/proteins. These revealed that Al3+ binds intracellular arginine, CrO42- enters the cell through sulfate transporters and oxidizes intracellular reduced thiols, and membrane-bound lipopolysaccharides protect the cell from UO22+ toxicity. In addition, the Tol outer membrane system contributed to the protection of cellular integrity from the toxic effects of all three metals. Likewise, we found evidence of regulation of lipid content in membranes under metal stress. Individually, RB-TnSeq and metabolomics are powerful tools to explore the impact various stresses have on biological systems. Here, we show that together they can be used synergistically to identify the molecular actors and mechanisms of these pertubations to an organism, furthering our understanding of how living systems interact with their environment. IMPORTANCE Studying microbial interactions with their environment can lead to a deeper understanding of biological molecular mechanisms. In this study, two global techniques, RB-TnSeq and activity metabolomics, were successfully used to probe the interactions between a metal-resistant microorganism, Pantoea sp. strain MT58, and metals contaminating a site where the organism can be located. A number of novel metal-microbe interactions were uncovered, including Al3+ toxicity targeting arginine synthesis, which could lead to a deeper understanding of the impact Al3+ contamination has on microbial communities as well as its impact on higher-level organisms, including plants for whom Al3+ contamination is an issue. Using multiomic approaches like the one described here is a way to further our understanding of microbial interactions and their impacts on the environment overall.