Rediversification following ecotype isolation reveals hidden adaptive potential.

Microbial communities play a critical role in ecological processes, and their diversity is key to their functioning. However, little is known about whether communities can regenerate ecological diversity following ecotype removal or extinction and how the rediversified communities would compare to the original ones. Here, we show that simple two-ecotype communities from the E. coli long-term evolution experiment (LTEE) consistently rediversified into two ecotypes following the isolation of one of the ecotypes, coexisting via negative frequency-dependent selection. Communities separated by more than 30,000 generations of evolutionary time rediversify in similar ways. The rediversified ecotype appears to share a number of growth traits with the ecotype it replaces. However, the rediversified community is also different from the original community in ways relevant to the mechanism of ecotype coexistence-for example, in stationary phase response and survival. We found substantial variation in the transcriptional states between the two original ecotypes, whereas the differences within the rediversified community were comparatively smaller, although the rediversified community showed unique patterns of differential expression. Our results suggest that evolution may leave room for alternative diversification processes even in a maximally reduced community of only two strains. We hypothesize that the presence of alternative evolutionary pathways may be even more pronounced in communities of many species where there are even more potential niches, highlighting an important role for perturbations, such as species removal, in evolving ecological communities.

Rediversification Following Ecotype Isolation Reveals Hidden Adaptive Potential.

Microbial communities play a critical role in ecological processes, and their diversity is key to their functioning. However, little is known about if communities can regenerate ecological diversity following species removal or extinction, and how the rediversified communities would compare to the original ones. Here we show that simple two-ecotype communities from the E. coli Long Term Evolution Experiment (LTEE) consistently rediversified into two ecotypes following the isolation of one of the ecotypes, coexisting via negative frequency-dependent selection. Communities separated by more than 30,000 generations of evolutionary time rediversify in similar ways. The rediversified ecotype appears to share a number of growth traits with the ecotype it replaces. However, the rediversified community is also different compared to the original community in ways relevant to the mechanism of ecotype coexistence, for example in stationary phase response and survival. We found substantial variation in the transcriptional states between the two original ecotypes, whereas the differences within the rediversified community were comparatively smaller, but with unique patterns of differential expression. Our results suggest that evolution may leave room for alternative diversification processes even in a maximally reduced community of only two strains. We hypothesize that the presence of alternative evolutionary pathways may be even more pronounced in communities of many species, highlighting an important role for perturbations, such as species removal, in evolving ecological communities.

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.

Quantifying the local adaptive landscape of a nascent bacterial community.

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

Global Analysis of Biomineralization Genes in Magnetospirillum magneticum AMB-1.

Magnetotactic bacteria (MTB) are a phylogenetically diverse group of bacteria remarkable for their ability to biomineralize magnetite (Fe3O4) or greigite (Fe3S4) in organelles called magnetosomes. The majority of genes required for magnetosome formation are encoded by a magnetosome gene island (MAI). Most previous genetic studies of MTB have focused on the MAI, using screens to identify key MAI genes or targeted genetics to isolate specific genes and their function in one specific growth condition. This is the first study that has taken an unbiased approach to look at many different growth conditions to reveal key genes both inside and outside the MAI. Here, we conducted random barcoded transposon mutagenesis (RB-TnSeq) in Magnetospirillum magneticum AMB-1. We generated a library of 184,710 unique strains in a wild-type background, generating ∼34 mutant strains for each gene. RB-TnSeq also allowed us to determine the essential gene set of AMB-1 under standard laboratory growth conditions. To pinpoint novel genes that are important for magnetosome formation, we subjected the library to magnetic selection screens under varied growth conditions. We compared biomineralization under standard growth conditions to biomineralization under high-iron and anaerobic conditions, respectively. Strains with transposon insertions in the MAI gene mamT had an exacerbated biomineralization defect under both high-iron and anaerobic conditions compared to standard conditions, adding to our knowledge of the role of MamT in magnetosome formation. Mutants in an ex-MAI gene, amb4151, are more magnetic than wild-type cells under anaerobic conditions. All three of these phenotypes were validated by creating a markerless deletion strain of the gene and evaluating with TEM imaging. Overall, our results indicate that growth conditions affect which genes are required for biomineralization and that some MAI genes may have more nuanced functions than was previously understood. IMPORTANCE Magnetotactic bacteria (MTB) are a group of bacteria that can form nano-sized crystals of magnetic minerals. MTB are likely an important part of their ecosystems, because they can account for up to a third of the microbial biomass in an aquatic habitat and consume large amounts of iron, potentially impacting the iron cycle. The ecology of MTB is relatively understudied; however, the cell biology and genetics of MTB have been studied for decades. Here, we leverage genetic studies of MTB to inform environmental studies. We expand the genetic toolset for studying MTB in the lab and identify novel genes, or functions of genes, that have an impact on biomineralization.

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.

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.

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.

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][Formula: see text]) or uranyl oxycation (U[VI][Formula: see text]). 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.

Mechanism of H2S Oxidation by the Dissimilatory Perchlorate-Reducing Microorganism Azospira suillum PS.

The genetic and biochemical basis of perchlorate-dependent H2S oxidation (PSOX) was investigated in the dissimilatory perchlorate-reducing microorganism (DPRM) Azospira suillum PS (PS). Previously, it was shown that all known DPRMs innately oxidize H2S, producing elemental sulfur (So). Although the process involving PSOX is thermodynamically favorable (ΔG°' = -206 kJ â‹… mol-1 H2S), the underlying biochemical and genetic mechanisms are currently unknown. Interestingly, H2S is preferentially utilized over physiological electron donors such as lactate or acetate although no growth benefit is obtained from the metabolism. Here, we determined that PSOX is due to a combination of enzymatic and abiotic interactions involving reactive intermediates of perchlorate respiration. Using various approaches, including barcode analysis by sequencing (Bar-seq), transcriptome sequencing (RNA-seq), and proteomics, along with targeted mutagenesis and biochemical characterization, we identified all facets of PSOX in PS. In support of our proposed model, deletion of identified upregulated PS genes traditionally known to be involved in sulfur redox cycling (e.g., Sox, sulfide:quinone reductase [SQR]) showed no defect in PSOX activity. Proteomic analysis revealed differential abundances of a variety of stress response metal efflux pumps and divalent heavy-metal transporter proteins, suggesting a general toxicity response. Furthermore, in vitro biochemical studies demonstrated direct PSOX mediated by purified perchlorate reductase (PcrAB) in the absence of other electron transfer proteins. The results of these studies support a model in which H2S oxidation is mediated by electron transport chain short-circuiting in the periplasmic space where the PcrAB directly oxidizes H2S to So The biogenically formed reactive intermediates (ClO2- and O2) subsequently react with additional H2S, producing polysulfide and So as end products.IMPORTANCE Inorganic sulfur compounds are widespread in nature, and microorganisms are central to their transformation, thereby playing a key role in the global sulfur cycle. Sulfur oxidation is mediated by a broad phylogenetic diversity of microorganisms, including anoxygenic phototrophs and either aerobic or anaerobic chemotrophs coupled to oxygen or nitrate respiration, respectively. Recently, perchlorate-respiring microorganisms were demonstrated to be innately capable of sulfur oxidation regardless of their phylogenetic affiliation. As recognition of the prevalence of these organisms intensifies, their role in global geochemical cycles is being queried. This is further highlighted by the recently recognized environmental pervasiveness of perchlorate not only across Earth but also throughout our solar system. The inferred importance of this metabolism not only is that it is a novel and previously unrecognized component of the global sulfur redox cycle but also is because of the recently demonstrated applicability of perchlorate respiration in the control of biogenic sulfide production in engineered environments such as oil reservoirs and wastewater treatment facilities, where excess H2S represents a significant environmental, process, and health risk, with associated costs approximating $90 billion annually.

A Comparison of the Costs and Benefits of Bacterial Gene Expression.

To study how a bacterium allocates its resources, we compared the costs and benefits of most (86%) of the proteins in Escherichia coli K-12 during growth in minimal glucose medium. The cost or investment in each protein was estimated from ribosomal profiling data, and the benefit of each protein was measured by assaying a library of transposon mutants. We found that proteins that are important for fitness are usually highly expressed, and 95% of these proteins are expressed at above 13 parts per million (ppm). Conversely, proteins that do not measurably benefit the host (with a benefit of less than 5% per generation) tend to be weakly expressed, with a median expression of 13 ppm. In aggregate, genes with no detectable benefit account for 31% of protein production, or about 22% if we correct for genetic redundancy. Although some of the apparently unnecessary expression could have subtle benefits in minimal glucose medium, the majority of the burden is due to genes that are important in other conditions. We propose that at least 13% of the cell's protein is "on standby" in case conditions change.

Novel Metal Cation Resistance Systems from Mutant Fitness Analysis of Denitrifying Pseudomonas stutzeri.

UNLABELLED: Metal ion transport systems have been studied extensively, but the specificity of a given transporter is often unclear from amino acid sequence data alone. In this study, predicted Cu(2+) and Zn(2+) resistance systems in Pseudomonas stutzeri strain RCH2 are compared with those experimentally implicated in Cu(2+) and Zn(2+) resistance, as determined by using a DNA-barcoded transposon mutant library. Mutant fitness data obtained under denitrifying conditions are combined with regulon predictions to yield a much more comprehensive picture of Cu(2+) and Zn(2+) resistance in strain RCH2. The results not only considerably expand what is known about well-established metal ion exporters (CzcCBA, CzcD, and CusCBA) and their accessory proteins (CzcI and CusF), they also reveal that isolates with mutations in some predicted Cu(2+) resistance systems do not show decreased fitness relative to the wild type when exposed to Cu(2+) In addition, new genes are identified that have no known connection to Zn(2+) (corB, corC, Psest_3226, Psest_3322, and Psest_0618) or Cu(2+) resistance (Mrp antiporter subunit gene, Psest_2850, and Psest_0584) but are crucial for resistance to these metal cations. Growth of individual deletion mutants lacking corB, corC, Psest_3226, or Psest_3322 confirmed the observed Zn-dependent phenotypes. Notably, to our knowledge, this is the first time a bacterial homolog of TMEM165, a human gene responsible for a congenital glycosylation disorder, has been deleted and the resulting strain characterized. Finally, the fitness values indicate Cu(2+)- and Zn(2+)-based inhibition of nitrite reductase and interference with molybdenum cofactor biosynthesis for nitrate reductase. These results extend the current understanding of Cu(2+) and Zn(2+) efflux and resistance and their effects on denitrifying metabolism. IMPORTANCE: In this study, genome-wide mutant fitness data in P. stutzeri RCH2 combined with regulon predictions identify several proteins of unknown function that are involved in resisting zinc and copper toxicity. For zinc, these include a member of the UPF0016 protein family that was previously implicated in Ca(2+)/H(+) antiport and a human congenital glycosylation disorder, CorB and CorC, which were previously linked to Mg(2+) transport, and Psest_3322 and Psest_0618, two proteins with no characterized homologs. Experiments using mutants lacking Psest_3226, Psest_3322, corB, corC, or czcI verified their proposed functions, which will enable future studies of these little-characterized zinc resistance determinants. Likewise, Psest_2850, annotated as an ion antiporter subunit, and the conserved hypothetical protein Psest_0584 are implicated in copper resistance. Physiological connections between previous studies and phenotypes presented here are discussed. Functional and mechanistic understanding of transport proteins improves the understanding of systems in which members of the same protein family, including those in humans, can have different functions.

Genetic dissection of chlorate respiration in Pseudomonas stutzeri†PDA reveals syntrophic (per)chlorate reduction.

Genes important for growth of Pseudomonas stutzeri PDA on chlorate were identified using a randomly DNA bar-coded transposon mutant library. During chlorate reduction, mutations in genes encoding the chlorate reductase clrABC, predicted molybdopterin cofactor chaperon clrD, molybdopterin biosynthesis and two genes of unknown function (clrE, clrF) had fitness defects in pooled mutant assays (Bar-seq). Markerless in-frame deletions confirmed that clrA, clrB and clrC were essential for chlorate reduction, while clrD, clrE and clrF had less severe growth defects. Interestingly, the key detoxification gene cld was essential for chlorate reduction in isogenic pure culture experiments, but showed only minor fitness defects in Bar-seq experiments. We hypothesized this was enabled through chlorite dismutation by the community, as most strains in the Bar-seq library contained an intact cld. In support of this, Δcld grew with wild-type PDA or ΔclrA, and purified Cld also restored growth to the Δcld mutant. Expanding on this, wild-type PDA and a Δcld mutant of the perchlorate reducer Azospira suillum PS grew on perchlorate in co-culture, but not individually. These results demonstrate that co-occurrence of cld and a chloroxyanion reductase within a single organism is not necessary and raises the possibility of syntrophic (per)chlorate respiration in the environment.

Determining Roles of Accessory Genes in Denitrification by Mutant Fitness Analyses.

Enzymes of the denitrification pathway play an important role in the global nitrogen cycle, including release of nitrous oxide, an ozone-depleting greenhouse gas. In addition, nitric oxide reductase, maturation factors, and proteins associated with nitric oxide detoxification are used by pathogens to combat nitric oxide release by host immune systems. While the core reductases that catalyze the conversion of nitrate to dinitrogen are well understood at a mechanistic level, there are many peripheral proteins required for denitrification whose basic function is unclear. A bar-coded transposon DNA library from Pseudomonas stutzeri strain RCH2 was grown under denitrifying conditions, using nitrate or nitrite as an electron acceptor, and also under molybdenum limitation conditions, with nitrate as the electron acceptor. Analysis of sequencing results from these growths yielded gene fitness data for 3,307 of the 4,265 protein-encoding genes present in strain RCH2. The insights presented here contribute to our understanding of how peripheral proteins contribute to a fully functioning denitrification pathway. We propose a new low-affinity molybdate transporter, OatABC, and show that differential regulation is observed for two MoaA homologs involved in molybdenum cofactor biosynthesis. We also propose that NnrS may function as a membrane-bound NO sensor. The dominant HemN paralog involved in heme biosynthesis is identified, and a CheR homolog is proposed to function in nitrate chemotaxis. In addition, new insights are provided into nitrite reductase redundancy, nitric oxide reductase maturation, nitrous oxide reductase maturation, and regulation.

The essential gene set of a photosynthetic organism.

Synechococcus elongatus PCC 7942 is a model organism used for studying photosynthesis and the circadian clock, and it is being developed for the production of fuel, industrial chemicals, and pharmaceuticals. To identify a comprehensive set of genes and intergenic regions that impacts fitness in S. elongatus, we created a pooled library of ∼ 250,000 transposon mutants and used sequencing to identify the insertion locations. By analyzing the distribution and survival of these mutants, we identified 718 of the organism's 2,723 genes as essential for survival under laboratory conditions. The validity of the essential gene set is supported by its tight overlap with well-conserved genes and its enrichment for core biological processes. The differences noted between our dataset and these predictors of essentiality, however, have led to surprising biological insights. One such finding is that genes in a large portion of the TCA cycle are dispensable, suggesting that S. elongatus does not require a cyclic TCA process. Furthermore, the density of the transposon mutant library enabled individual and global statements about the essentiality of noncoding RNAs, regulatory elements, and other intergenic regions. In this way, a group I intron located in tRNA(Leu), which has been used extensively for phylogenetic studies, was shown here to be essential for the survival of S. elongatus. Our survey of essentiality for every locus in the S. elongatus genome serves as a powerful resource for understanding the organism's physiology and defines the essential gene set required for the growth of a photosynthetic organism.

Novel mechanism for scavenging of hypochlorite involving a periplasmic methionine-rich Peptide and methionine sulfoxide reductase.

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

Rapid quantification of mutant fitness in diverse bacteria by sequencing randomly bar-coded transposons.

UNLABELLED: Transposon mutagenesis with next-generation sequencing (TnSeq) is a powerful approach to annotate gene function in bacteria, but existing protocols for TnSeq require laborious preparation of every sample before sequencing. Thus, the existing protocols are not amenable to the throughput necessary to identify phenotypes and functions for the majority of genes in diverse bacteria. Here, we present a method, random bar code transposon-site sequencing (RB-TnSeq), which increases the throughput of mutant fitness profiling by incorporating random DNA bar codes into Tn5 and mariner transposons and by using bar code sequencing (BarSeq) to assay mutant fitness. RB-TnSeq can be used with any transposon, and TnSeq is performed once per organism instead of once per sample. Each BarSeq assay requires only a simple PCR, and 48 to 96 samples can be sequenced on one lane of an Illumina HiSeq system. We demonstrate the reproducibility and biological significance of RB-TnSeq with Escherichia coli, Phaeobacter inhibens, Pseudomonas stutzeri, Shewanella amazonensis, and Shewanella oneidensis. To demonstrate the increased throughput of RB-TnSeq, we performed 387 successful genome-wide mutant fitness assays representing 130 different bacterium-carbon source combinations and identified 5,196 genes with significant phenotypes across the five bacteria. In P. inhibens, we used our mutant fitness data to identify genes important for the utilization of diverse carbon substrates, including a putative d-mannose isomerase that is required for mannitol catabolism. RB-TnSeq will enable the cost-effective functional annotation of diverse bacteria using mutant fitness profiling. IMPORTANCE: A large challenge in microbiology is the functional assessment of the millions of uncharacterized genes identified by genome sequencing. Transposon mutagenesis coupled to next-generation sequencing (TnSeq) is a powerful approach to assign phenotypes and functions to genes. However, the current strategies for TnSeq are too laborious to be applied to hundreds of experimental conditions across multiple bacteria. Here, we describe an approach, random bar code transposon-site sequencing (RB-TnSeq), which greatly simplifies the measurement of gene fitness by using bar code sequencing (BarSeq) to monitor the abundance of mutants. We performed 387 genome-wide fitness assays across five bacteria and identified phenotypes for over 5,000 genes. RB-TnSeq can be applied to diverse bacteria and is a powerful tool to annotate uncharacterized genes using phenotype data.

The genetic basis of energy conservation in the sulfate-reducing bacterium Desulfovibrio alaskensis G20.

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.