Search

Abstract A mechanistic understanding of host-microbe interactions in the gut microbiome is hindered by poorly annotated bacterial genomes. While functional genomics can generate large gene-to-phenotype datasets to accelerate functional discovery, their applications to study gut anaerobes have been limited. For instance, most gain-of-function screens of gut-derived genes have been performed in Escherichia coli and assayed in a small number of conditions. To address these challenges, we develop Barcoded Overexpression BActerial shotgun library sequencing (Boba-seq). We demonstrate the power of this approach by assaying genes from diverse gut Bacteroidales overexpressed in Bacteroides thetaiotaomicron. From hundreds of experiments, we identify new functions and phenotypes for 29 genes important for carbohydrate metabolism or tolerance to antibiotics or bile salts. Highlights include the discovery of a d-glucosamine kinase, a raffinose transporter, and several routes that increase tolerance to ceftriaxone and bile salts through lipid biosynthesis. This approach can be readily applied to develop screens in other strains and additional phenotypic assays.

ABSTRACT Archaea are widespread in the environment and play fundamental roles in diverse ecosystems; however, characterization of their unique biology requires advanced tools. This is particularly challenging when characterizing gene function. Here, we generate randomly barcoded transposon libraries in the model methanogenic archaeon Methanococcus maripaludis and use high-throughput growth methods to conduct fitness assays (RB-TnSeq) across over 100 unique growth conditions. Using our approach, we identified new genes involved in nutrient utilization and response to oxidative stress. We identified novel genes for the usage of diverse nitrogen sources in M. maripaludis including a putative regulator of alanine deamination and molybdate transporters important for nitrogen fixation. Furthermore, leveraging the fitness data, we inferred that M. maripaludis can utilize additional nitrogen sources including ʟ-glutamine, ᴅ-glucuronamide, and adenosine. Under autotrophic growth conditions, we identified a gene encoding a domain of unknown function (DUF166) that is important for fitness and hypothesize that it has an accessory role in carbon dioxide assimilation. Finally, comparing fitness costs of oxygen versus sulfite stress, we identified a previously uncharacterized class of dissimilatory sulfite reductase-like proteins (Dsr-LP; group IIId) that is important during growth in the presence of sulfite. When overexpressed, Dsr-LP conferred sulfite resistance and enabled use of sulfite as the sole sulfur source. The high-throughput approach employed here allowed for generation of a large-scale data set that can be used as a resource to further understand gene function and metabolism in the archaeal domain.IMPORTANCEArchaea are widespread in the environment, yet basic aspects of their biology remain underexplored. To address this, we apply randomly barcoded transposon libraries (RB-TnSeq) to the model archaeon Methanococcus maripaludis. RB-TnSeq coupled with high-throughput growth assays across over 100 unique conditions identified roles for previously uncharacterized genes, including several encoding proteins with domains of unknown function (DUFs). We also expand on our understanding of carbon and nitrogen metabolism and characterize a group IIId dissimilatory sulfite reductase-like protein as a functional sulfite reductase. This data set encompasses a wide range of additional conditions including stress, nitrogen fixation, amino acid supplementation, and autotrophy, thus providing an extensive data set for the archaeal community to mine for characterizing additional genes of unknown function.

Studying plant–microbe–soil interactions is challenging due to their high complexity and variability in natural ecosystems. While fabricated ecosystems provide opportunities to recapitulate aspects of these systems in reduced complexity and controlled environments, inoculation can be a significant source of variation. To tackle this, we evaluated how different bacterial inoculation practices and plant harvesting time points affect the reproducibility of a synthetic microbial community (SynCom) in association with the model grass Brachypodium distachyon. We tested three microbial inoculation practices, seed inoculation, transplant inoculation, and seedling inoculation, and two harvesting points, early (14-day-old plants) and late (21 days postinoculation). We grew our plants and bacterial strains in sterile devices (EcoFABs) and characterized the microbial community from root, rhizosphere, and sand using 16S ribosomal RNA gene sequencing. The results showed that inoculation practices significantly affected the rhizosphere microbial community when harvesting at an early time point but not at the late time point. As the SynCom showed a persistent association with B. distachyon at 21 days postinoculation regardless of inoculation practices, we assessed the reproducibility of each inoculation method and found that transplant inoculation showed the highest reproducibility. Moreover, plant biomass was not adversely affected by transplant inoculation treatment. We concluded that bacterial inoculation while transplanting coupled with a later harvesting time point gives the most reproducible microbial community in the EcoFAB−B. distachyon−SynCom fabricated ecosystem. We recommend this method as a standardized protocol for use with fabricated ecosystem experimental systems. [Graphic: see text] This manuscript has been authored by an author at Lawrence Berkeley National Laboratory under Contract number DE-AC02-05CH11231 with the U.S. Department of Energy. The U.S. Government retains, and the publisher, by accepting the article for publication, acknowledges, that the U.S. Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for U.S. Government purposes. This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Abstract 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.

Summary: 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.

Summary: 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.

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.

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.

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. Author summary For many microbes, we know little about them beyond their genome sequences. In principle, we could use genome sequences to predict microbes’ traits, such as which carbon sources they can eat, but first we need to identify more of the genes involved. We built an automated tool, GapMind, to annotate the transporters and enzymes for utilizing 62 common carbon sources, and used GapMind to identify gaps: transporters or enzymes that should be present, to explain how a bacterium uses a carbon source, but could not be found in the genome. By comparing these gaps to large-scale genetic data for 29 bacteria, we identified hundreds of novel transporters and enzymes, and a new metabolic pathway for consuming glucosamine. When we added these novel genes to GapMind, its results for diverse bacteria and archaea improved significantly.

ABSTRACT Rhodanobacter species dominate in the Oak Ridge Reservation (ORR) subsurface environments contaminated with acids, nitrate, metal radionuclides, and other heavy metals. To uncover the genomic features underlying adaptations to these mixed-waste environments and to guide genetic tool development, we sequenced the whole genomes of eight Rhodanobacter strains isolated from the ORR site. The genome sizes ranged from 3.9 to 4.2 Mb harboring 3,695 to 4,035 protein-coding genes and GC contents approximately 67%. Seven strains were classified as R. denitrificans and one strain, FW510-R12, as R. thiooxydans based on full length 16S rRNA sequences. According to gene annotation, the top two Cluster of Orthologous Groups (COGs) with high pan-genome expansion rates (Pan/Core gene ratio) were “replication, recombination and repair” and “defense mechanisms.” The denitrifying genes had high DNA homologies except the predicted protein structure variances in NosZ. In contrast, heavy metal resistance genes were diverse with between 7 to 34% of them were located in genomic islands, and these results suggested origins from horizontal gene transfer. Analysis of the methylation patterns in four strains revealed the unique 5mC methylation motifs. Most orthologs (78%) had ratios of nonsynonymous to synonymous substitutions (dN/dS) less than one when compared to the type strain 2APBS1, suggesting the prevalence of negative selection. Overall, the results provide evidence for the important roles of horizontal gene transfer and negative selection in genomic adaptation at the contaminated field site. The complex restriction-modification system genes and the unique methylation motifs in Rhodanobacter strains suggest the potential recalcitrance to genetic manipulation. IMPORTANCE Despite the dominance of Rhodanobacter species in the subsurface of the contaminated Oak Ridge Reservation (ORR) site, very little is known about the mechanisms underlying their adaptions to the various stressors present at ORR. Recently, multiple Rhodanobacter strains have been isolated from the ORR groundwater samples from several wells with varying geochemical properties. Using Illumina, PacBio, and Oxford Nanopore sequencing platforms, we obtained the whole genome sequences of eight Rhodanobacter strains. Comparison of the whole genomes demonstrated the genetic diversity, and analysis of the long nanopore reads revealed the heterogeneity of methylation patterns in strains isolated from the same well. Although all strains contained a complete set of denitrifying genes, the predicted tertiary structures of NosZ differed. The sequence comparison results demonstrate the important roles of horizontal gene transfer and negative selection in adaptation. In addition, these strains may be recalcitrant to genetic manipulation due to the complex restriction-modification systems and methylations.

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.

ABSTRACT Plant-pathogenic Ralstonia spp. colonize plant xylem and cause wilt diseases on a broad range of host plants. To identify genes that promote growth of diverse Ralstonia strains in xylem sap from tomato plants, we performed genome-scale genetic screens (random barcoded transposon mutant sequencing screens [RB-TnSeq]) in three strains spanning the genetic, geographical, and physiological range of plant-pathogenic Ralstonia: Ralstonia solanacearum IBSBF1503, Ralstonia pseudosolanacearum GMI1000, and Ralstonia syzygii PSI07. Contrasting mutant fitness phenotypes in culture media versus in xylem sap suggest that Ralstonia strains are adapted to ex vivo xylem sap and that culture media impose foreign selective pressures. Although wild-type Ralstonia grew in sap and in rich medium with similar doubling times and to a similar carrying capacity, more genes were essential for growth in sap than in rich medium. Each strain required many genes associated with envelope remodeling and repair processes for full fitness in xylem sap. These genes were associated with peptidoglycan peptide formation (murI), secretion of periplasmic proteins (tatC), periplasmic protein folding (dsbA), synthesis of osmoregulated periplasmic glucans (mdoGH), and lipopolysaccharide (LPS) biosynthesis. Mutant strains with mutations in four genes had strong, sap-specific fitness defects in all strain backgrounds: murI, thiC, purU, and a lipoprotein (RSc2007). Many amino acid biosynthesis genes were required for fitness in both minimal medium and xylem sap. Multiple mutants with insertions in virulence regulators had gains of fitness in culture media and neutral fitness in sap. Our genome-scale genetic screen identified Ralstonia fitness factors that promote growth in xylem sap, an ecologically relevant condition. IMPORTANCE Traditional transposon mutagenesis genetic screens pioneered molecular plant pathology and identified core virulence traits like the type III secretion system. TnSeq approaches that leverage next-generation sequencing to rapidly quantify transposon mutant phenotypes are ushering in a new wave of biological discovery. Here, we have adapted a genome-scale approach, random barcoded transposon mutant sequencing (RB-TnSeq), to discover fitness factors that promote growth of three related bacterial strains in a common niche, tomato xylem sap. Fitness of the wild type and mutants show that Ralstonia spp. are adapted to grow well in xylem sap from their natural host plant, tomato. Our screen identified multiple sap-specific fitness factors with roles in maintaining the bacterial envelope. These factors include putative adaptations to resist plant defenses that may include antimicrobial proteins and specialized metabolites that damage bacterial membranes.

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.

Gain of function methods based on gene overexpression are not easily applied to high-throughput screening across different experimental conditions. Here, the authors present Dub-seq, which separates overexpression library characterization from functional screening and uses random DNA barcodes to increase the experimental throughput.

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.

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 (

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.

ABSTRACT 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.

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.

ABSTRACT Plasmid-mediated horizontal gene transfer (HGT) is a major driver of genetic diversity in bacteria. We experimentally validated the function of a putative mercury resistance operon present on an abundant 8-kbp native plasmid found in groundwater samples without detectable levels of mercury. Phylogenetic analyses of the plasmid-encoded mercury reductases from the studied groundwater site show them to be distinct from those reported in proximal metal-contaminated sites. We synthesized the entire native plasmid and demonstrated that the plasmid was sufficient to confer functional mercury resistance in Escherichia coli. Given the possibility that natural transformation is a prevalent HGT mechanism in the low-cell-density environments of groundwaters, we also assayed bacterial strains from this environment for competence. We used the native plasmid-encoded metal resistance to design a screen and identified 17 strains positive for natural transformation. We selected 2 of the positive strains along with a model bacterium to fully confirm HGT via natural transformation. From an ecological perspective, the role of the native plasmid population in providing advantageous traits combined with the microbiome’s capacity to take up environmental DNA enables rapid adaptation to environmental stresses. IMPORTANCE Horizontal transfer of mobile genetic elements via natural transformation has been poorly understood in environmental microbes. Here, we confirm the functionality of a native plasmid-encoded mercury resistance operon in a model microbe and then query for the dissemination of this resistance trait via natural transformation into environmental bacterial isolates. We identified 17 strains including Gram-positive and Gram-negative bacteria to be naturally competent. These strains were able to successfully take up the plasmid DNA and obtain a clear growth advantage in the presence of mercury. Our study provides important insights into gene dissemination via natural transformation enabling rapid adaptation to dynamic stresses in groundwater environments.

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.

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.

ABSTRACT 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.

ABSTRACT Despite intensive study for 50 years, the biochemical and genetic links between lysine metabolism and central metabolism in Pseudomonas putida remain unresolved. To establish these biochemical links, we leveraged random barcode transposon sequencing (RB-TnSeq), a genome-wide assay measuring the fitness of thousands of genes in parallel, to identify multiple novel enzymes in both l- and d-lysine metabolism. We first describe three pathway enzymes that catabolize l-2-aminoadipate (l-2AA) to 2-ketoglutarate (2KG), connecting d-lysine to the TCA cycle. One of these enzymes, P. putida 5260 (PP_5260), contains a DUF1338 domain, representing a family with no previously described biological function. Our work also identified the recently described coenzyme A (CoA)-independent route of l-lysine degradation that results in metabolization to succinate. We expanded on previous findings by demonstrating that glutarate hydroxylase CsiD is promiscuous in its 2-oxoacid selectivity. Proteomics of selected pathway enzymes revealed that expression of catabolic genes is highly sensitive to the presence of particular pathway metabolites, implying intensive local and global regulation. This work demonstrated the utility of RB-TnSeq for discovering novel metabolic pathways in even well-studied bacteria, as well as its utility a powerful tool for validating previous research. IMPORTANCE P. putida lysine metabolism can produce multiple commodity chemicals, conferring great biotechnological value. Despite much research, the connection of lysine catabolism to central metabolism in P. putida remained undefined. Here, we used random barcode transposon sequencing to fill the gaps of lysine metabolism in P. putida. We describe a route of 2-oxoadipate (2OA) catabolism, which utilizes DUF1338-containing protein P. putida 5260 (PP_5260) in bacteria. Despite its prevalence in many domains of life, DUF1338-containing proteins have had no known biochemical function. We demonstrate that PP_5260 is a metalloenzyme which catalyzes an unusual route of decarboxylation of 2OA to d-2-hydroxyglutarate (d-2HG). Our screen also identified a recently described novel glutarate metabolic pathway. We validate previous results and expand the understanding of glutarate hydroxylase CsiD by showing that can it use either 2OA or 2KG as a cosubstrate. Our work demonstrated that biological novelty can be rapidly identified using unbiased experimental genetics and that RB-TnSeq can be used to rapidly validate previous results.

ABSTRACT Using genome-wide mutant fitness assays in diverse bacteria, we identified novel oxidative pathways for the catabolism of 2-deoxy-d-ribose and 2-deoxy-d-ribonate. We propose that deoxyribose is oxidized to deoxyribonate, oxidized to ketodeoxyribonate, and cleaved to acetyl coenzyme A (acetyl-CoA) and glyceryl-CoA. We have genetic evidence for this pathway in three genera of bacteria, and we confirmed the oxidation of deoxyribose to ketodeoxyribonate in vitro. In Pseudomonas simiae, the expression of enzymes in the pathway is induced by deoxyribose or deoxyribonate, while in Paraburkholderia bryophila and in Burkholderia phytofirmans, the pathway proceeds in parallel with the known deoxyribose 5-phosphate aldolase pathway. We identified another oxidative pathway for the catabolism of deoxyribonate, with acyl-CoA intermediates, in Klebsiella michiganensis. Of these four bacteria, only P. simiae relies entirely on an oxidative pathway to consume deoxyribose. The deoxyribose dehydrogenase of P. simiae is either nonspecific or evolved recently, as this enzyme is very similar to a novel vanillin dehydrogenase from Pseudomonas putida that we identified. So, we propose that these oxidative pathways evolved primarily to consume deoxyribonate, which is a waste product of metabolism. IMPORTANCE Deoxyribose is one of the building blocks of DNA and is released when cells die and their DNA degrades. We identified a bacterium that can grow with deoxyribose as its sole source of carbon even though its genome does not contain any of the known genes for breaking down deoxyribose. By growing many mutants of this bacterium together on deoxyribose and using DNA sequencing to measure the change in the mutants’ abundance, we identified multiple protein-coding genes that are required for growth on deoxyribose. Based on the similarity of these proteins to enzymes of known function, we propose a 6-step pathway in which deoxyribose is oxidized and then cleaved. Diverse bacteria use a portion of this pathway to break down a related compound, deoxyribonate, which is a waste product of metabolism. Our study illustrates the utility of large-scale bacterial genetics to identify previously unknown metabolic pathways.

ABSTRACT Transposon mutagenesis coupled to next-generation sequencing (TnSeq) is a powerful approach for discovering the functions of bacterial genes. However, the development of a suitable TnSeq strategy for a given bacterium can be costly and time-consuming. To meet this challenge, we describe a part-based strategy for constructing libraries of hundreds of transposon delivery vectors, which we term “magic pools.” Within a magic pool, each transposon vector has a different combination of upstream sequences (promoters and ribosome binding sites) and antibiotic resistance markers as well as a random DNA barcode sequence, which allows the tracking of each vector during mutagenesis experiments. To identify an efficient vector for a given bacterium, we mutagenize it with a magic pool and sequence the resulting insertions; we then use this efficient vector to generate a large mutant library. We used the magic pool strategy to construct transposon mutant libraries in five genera of bacteria, including three genera of the phylum Bacteroidetes. IMPORTANCE Molecular genetics is indispensable for interrogating the physiology of bacteria. However, the development of a functional genetic system for any given bacterium can be time-consuming. Here, we present a streamlined approach for identifying an effective transposon mutagenesis system for a new bacterium. Our strategy first involves the construction of hundreds of different transposon vector variants, which we term a “magic pool.” The efficacy of each vector in a magic pool is monitored in parallel using a unique DNA barcode that is introduced into each vector design. Using archived DNA “parts,” we next reassemble an effective vector for making a whole-genome transposon mutant library that is suitable for large-scale interrogation of gene function using competitive growth assays. Here, we demonstrate the utility of the magic pool system to make mutant libraries in five genera of bacteria.

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.

ABSTRACT Unraveling the drivers controlling the response and adaptation of biological communities to environmental change, especially anthropogenic activities, is a central but poorly understood issue in ecology and evolution. Comparative genomics studies suggest that lateral gene transfer (LGT) is a major force driving microbial genome evolution, but its role in the evolution of microbial communities remains elusive. To delineate the importance of LGT in mediating the response of a groundwater microbial community to heavy metal contamination, representative Rhodanobacter reference genomes were sequenced and compared to shotgun metagenome sequences. 16S rRNA gene-based amplicon sequence analysis indicated that Rhodanobacter populations were highly abundant in contaminated wells with low pHs and high levels of nitrate and heavy metals but remained rare in the uncontaminated wells. Sequence comparisons revealed that multiple geochemically important genes, including genes encoding Fe2+/Pb2+ permeases, most denitrification enzymes, and cytochrome c553, were native to Rhodanobacter and not subjected to LGT. In contrast, the Rhodanobacter pangenome contained a recombinational hot spot in which numerous metal resistance genes were subjected to LGT and/or duplication. In particular, Co2+/Zn2+/Cd2+ efflux and mercuric resistance operon genes appeared to be highly mobile within Rhodanobacter populations. Evidence of multiple duplications of a mercuric resistance operon common to most Rhodanobacter strains was also observed. Collectively, our analyses indicated the importance of LGT during the evolution of groundwater microbial communities in response to heavy metal contamination, and a conceptual model was developed to display such adaptive evolutionary processes for explaining the extreme dominance of Rhodanobacter populations in the contaminated groundwater microbiome. IMPORTANCE Lateral gene transfer (LGT), along with positive selection and gene duplication, are the three main mechanisms that drive adaptive evolution of microbial genomes and communities, but their relative importance is unclear. Some recent studies suggested that LGT is a major adaptive mechanism for microbial populations in response to changing environments, and hence, it could also be critical in shaping microbial community structure. However, direct evidence of LGT and its rates in extant natural microbial communities in response to changing environments is still lacking. Our results presented in this study provide explicit evidence that LGT played a crucial role in driving the evolution of a groundwater microbial community in response to extreme heavy metal contamination. It appears that acquisition of genes critical for survival, growth, and reproduction via LGT is the most rapid and effective way to enable microorganisms and associated microbial communities to quickly adapt to abrupt harsh environmental stresses.

ABSTRACT Reactive chlorine species (RCS) defense mechanisms are important for bacterial fitness in diverse environments. In addition to the anthropogenic use of RCS in the form of bleach, these compounds are also produced naturally through photochemical reactions of natural organic matter and in vivo by the mammalian immune system in response to invading microorganisms. To gain insight into bacterial RCS defense mechanisms, we investigated Azospira suillum strain PS, which produces periplasmic RCS as an intermediate of perchlorate respiration. Our studies identified an RCS response involving an RCS stress-sensing sigma/anti-sigma factor system (SigF/NrsF), a soluble hypochlorite-scavenging methionine-rich periplasmic protein (MrpX), and a putative periplasmic methionine sulfoxide reductase (YedY1). We investigated the underlying mechanism by phenotypic characterization of appropriate gene deletions, chemogenomic profiling of barcoded transposon pools, transcriptome sequencing, and biochemical assessment of methionine oxidation. Our results demonstrated that SigF was specifically activated by RCS and initiated the transcription of a small regulon centering around yedY1 and mrpX. A yedY1 paralog (yedY2) was found to have a similar fitness to yedY1 despite not being regulated by SigF. Markerless deletions of yedY2 confirmed its synergy with the SigF regulon. MrpX was strongly induced and rapidly oxidized by RCS, especially hypochlorite. Our results suggest a mechanism involving hypochlorite scavenging by sacrificial oxidation of the MrpX in the periplasm. Reduced MrpX is regenerated by the YedY methionine sulfoxide reductase activity. The phylogenomic distribution of this system revealed conservation in several Proteobacteria of clinical importance, including uropathogenic Escherichia coli and Brucella spp., implying a putative role in immune response evasion in vivo. IMPORTANCE Bacteria are often stressed in the environment by reactive chlorine species (RCS) of either anthropogenic or natural origin, but little is known of the defense mechanisms they have evolved. Using a microorganism that generates RCS internally as part of its respiratory process allowed us to uncover a novel defense mechanism based on RCS scavenging by reductive reaction with a sacrificial methionine-rich peptide and redox recycling through a methionine sulfoxide reductase. This system is conserved in a broad diversity of organisms, including some of clinical importance, invoking a possible important role in innate immune system evasion.

Sulfate-reducing bacteria play major roles in the global carbon and sulfur cycles, but it remains unclear how reducing sulfate yields energy. To determine the genetic basis of energy conservation, we measured the fitness of thousands of pooled mutants of Desulfovibrio alaskensis G20 during growth in 12 different combinations of electron donors and acceptors. We show that ion pumping by the ferredoxin:NADH oxidoreductase Rnf is required whenever substrate-level phosphorylation is not possible. The uncharacterized complex Hdr/flox-1 (Dde_1207:13) is sometimes important alongside Rnf and may perform an electron bifurcation to generate more reduced ferredoxin from NADH to allow further ion pumping. Similarly, during the oxidation of malate or fumarate, the electron-bifurcating transhydrogenase NfnAB-2 (Dde_1250:1) is important and may generate reduced ferredoxin to allow additional ion pumping by Rnf. During formate oxidation, the periplasmic [NiFeSe] hydrogenase HysAB is required, which suggests that hydrogen forms in the periplasm, diffuses to the cytoplasm, and is used to reduce ferredoxin, thus providing a substrate for Rnf. During hydrogen utilization, the transmembrane electron transport complex Tmc is important and may move electrons from the periplasm into the cytoplasmic sulfite reduction pathway. Finally, mutants of many other putative electron carriers have no clear phenotype, which suggests that they are not important under our growth conditions, although we cannot rule out genetic redundancy.

ABSTRACT Transcription start sites (TSSs) lying inside annotated genes, on the same or opposite strand, have been observed in diverse bacteria, but the function of these unexpected transcripts is unclear. Here, we use the metal-reducing bacterium Shewanella oneidensis MR-1 and its relatives to study the evolutionary conservation of unexpected TSSs. Using high-resolution tiling microarrays and 5′-end RNA sequencing, we identified 2,531 TSSs in S. oneidensis MR-1, of which 18% were located inside coding sequences (CDSs). Comparative transcriptome analysis with seven additional Shewanella species revealed that the majority (76%) of the TSSs within the upstream regions of annotated genes (gTSSs) were conserved. Thirty percent of the TSSs that were inside genes and on the sense strand (iTSSs) were also conserved. Sequence analysis around these iTSSs showed conserved promoter motifs, suggesting that many iTSS are under purifying selection. Furthermore, conserved iTSSs are enriched for regulatory motifs, suggesting that they are regulated, and they tend to eliminate polar effects, which confirms that they are functional. In contrast, the transcription of antisense TSSs located inside CDSs (aTSSs) was significantly less likely to be conserved (22%). However, aTSSs whose transcription was conserved often have conserved promoter motifs and drive the expression of nearby genes. Overall, our findings demonstrate that some internal TSSs are conserved and drive protein expression despite their unusual locations, but the majority are not conserved and may reflect noisy initiation of transcription rather than a biological function. IMPORTANCE The first step of gene expression is the initiation of transcription from promoters, which have been traditionally thought to be located upstream of genes. Recently, studies showed that in diverse bacteria, promoters are often located inside genes. It has not been clear if these unexpected promoters are important to the organism or if they result from transcriptional noise. Here, we identify and examine promoters in eight related bacterial species. Promoters that lie within genes on the sense strand are often conserved as locations and in their sequences. Furthermore, these promoters often affect the bacterium’s growth. Thus, many of these unexpected promoters are likely functional. Fewer promoters that lie within genes on the antisense strand are conserved, but the conserved ones seem to drive the expression of nearby genes.