Genetic manipulation of gut microbes enables single-gene interrogation in a complex microbiome.

Hundreds of microbiota genes are associated with host biology/disease. Unraveling the causal contribution of a microbiota gene to host biology remains difficult because many are encoded by nonmodel gut commensals and not genetically targetable. A general approach to identify their gene transfer methodology and build their gene manipulation tools would enable mechanistic dissections of their impact on host physiology. We developed a pipeline that identifies the gene transfer methods for multiple nonmodel microbes spanning five phyla, and we demonstrated the utility of their genetic tools by modulating microbiome-derived short-chain fatty acids and bile acids in vitro and in the host. In a proof-of-principle study, by deleting a commensal gene for bile acid synthesis in a complex microbiome, we discovered an intriguing role of this gene in regulating colon inflammation. This technology will enable genetically engineering the nonmodel gut microbiome and facilitate mechanistic dissection of microbiota-host interactions.

Eukaryotic Acquisition of a Bacterial Operon.

Operons are a hallmark of bacterial genomes, where they allow concerted expression of functionally related genes as single polycistronic transcripts. They are rare in eukaryotes, where each gene usually drives expression of its own independent messenger RNAs. Here, we report the horizontal operon transfer of a siderophore biosynthesis pathway from relatives of Escherichia coli into a group of budding yeast taxa. We further show that the co-linearly arranged secondary metabolism genes are expressed, exhibit eukaryotic transcriptional features, and enable the sequestration and uptake of iron. After transfer, several genetic changes occurred during subsequent evolution, including the gain of new transcription start sites that were sometimes within protein-coding sequences, acquisition of polyadenylation sites, structural rearrangements, and integration of eukaryotic genes into the cluster. We conclude that the genes were likely acquired as a unit, modified for eukaryotic gene expression, and maintained by selection to adapt to the highly competitive, iron-limited environment.

Peptidoglycan Production by an Insect-Bacterial Mosaic.

Peptidoglycan (PG) is a defining feature of bacteria, involved in cell division, shape, and integrity. We previously reported that several genes related to PG biosynthesis were horizontally transferred from bacteria to the nuclear genome of mealybugs. Mealybugs are notable for containing a nested bacteria-within-bacterium endosymbiotic structure in specialized insect cells, where one bacterium, Moranella, lives in the cytoplasm of another bacterium, Tremblaya. Here we show that horizontally transferred genes on the mealybug genome work together with genes retained on the Moranella genome to produce a PG layer exclusively at the Moranella cell periphery. Furthermore, we show that an insect protein encoded by a horizontally transferred gene of bacterial origin is transported into the Moranella cytoplasm. These results provide a striking parallel to the genetic and biochemical mosaicism found in organelles, and prove that multiple horizontally transferred genes can become integrated into a functional pathway distributed between animal and bacterial endosymbiont genomes.

A Comprehensive, CRISPR-based Functional Analysis of Essential Genes in Bacteria.

Essential gene functions underpin the core reactions required for cell viability, but their contributions and relationships are poorly studied in vivo. Using CRISPR interference, we created knockdowns of every essential gene in Bacillus subtilis and probed their phenotypes. Our high-confidence essential gene network, established using chemical genomics, showed extensive interconnections among distantly related processes and identified modes of action for uncharacterized antibiotics. Importantly, mild knockdown of essential gene functions significantly reduced stationary-phase survival without affecting maximal growth rate, suggesting that essential protein levels are set to maximize outgrowth from stationary phase. Finally, high-throughput microscopy indicated that cell morphology is relatively insensitive to mild knockdown but profoundly affected by depletion of gene function, revealing intimate connections between cell growth and shape. Our results provide a framework for systematic investigation of essential gene functions in vivo broadly applicable to diverse microorganisms and amenable to comparative analysis.

Intragenic DNA inversions expand bacterial coding capacity.

Bacterial populations that originate from a single bacterium are not strictly clonal and often contain subgroups with distinct phenotypes1. Bacteria can generate heterogeneity through phase variation-a preprogrammed, reversible mechanism that alters gene expression levels across a population1. One well-studied type of phase variation involves enzyme-mediated inversion of specific regions of genomic DNA2. Frequently, these DNA inversions flip the orientation of promoters, turning transcription of adjacent coding regions on or off2. Through this mechanism, inversion can affect fitness, survival or group dynamics3,4. Here, we describe the development of PhaVa, a computational tool that identifies DNA inversions using long-read datasets. We also identify 372 'intragenic invertons', a novel class of DNA inversions found entirely within genes, in genomes of bacterial and archaeal isolates. Intragenic invertons allow a gene to encode two or more versions of a protein by flipping a DNA sequence within the coding region, thereby increasing coding capacity without increasing genome size. We validate ten intragenic invertons in the gut commensal Bacteroides thetaiotaomicron, and experimentally characterize an intragenic inverton in the thiamine biosynthesis gene thiC.

Functional and evolutionary significance of unknown genes from uncultivated taxa.

Many of the Earth's microbes remain uncultured and understudied, limiting our understanding of the functional and evolutionary aspects of their genetic material, which remain largely overlooked in most metagenomic studies1. Here we analysed 149,842 environmental genomes from multiple habitats2-6 and compiled a curated catalogue of 404,085 functionally and evolutionarily significant novel (FESNov) gene families exclusive to uncultivated prokaryotic taxa. All FESNov families span multiple species, exhibit strong signals of purifying selection and qualify as new orthologous groups, thus nearly tripling the number of bacterial and archaeal gene families described to date. The FESNov catalogue is enriched in clade-specific traits, including 1,034 novel families that can distinguish entire uncultivated phyla, classes and orders, probably representing synapomorphies that facilitated their evolutionary divergence. Using genomic context analysis and structural alignments we predicted functional associations for 32.4% of FESNov families, including 4,349 high-confidence associations with important biological processes. These predictions provide a valuable hypothesis-driven framework that we used for experimental validatation of a new gene family involved in cell motility and a novel set of antimicrobial peptides. We also demonstrate that the relative abundance profiles of novel families can discriminate between environments and clinical conditions, leading to the discovery of potentially new biomarkers associated with colorectal cancer. We expect this work to enhance future metagenomics studies and expand our knowledge of the genetic repertory of uncultivated organisms.

In situ targeted base editing of bacteria in the mouse gut.

Microbiome research is now demonstrating a growing number of bacterial strains and genes that affect our health1. Although CRISPR-derived tools have shown great success in editing disease-driving genes in human cells2, we currently lack the tools to achieve comparable success for bacterial targets in situ. Here we engineer a phage-derived particle to deliver a base editor and modify Escherichia coli colonizing the mouse gut. Editing of a ß-lactamase gene in a model E. coli strain resulted in a median editing efficiency of 93% of the target bacterial population with a single dose. Edited bacteria were stably maintained in the mouse gut for at least 42 days following treatment. This was achieved using a non-replicative DNA vector, preventing maintenance and dissemination of the payload. We then leveraged this approach to edit several genes of therapeutic relevance in E. coli and Klebsiella pneumoniae strains in vitro and demonstrate in situ editing of a gene involved in the production of curli in a pathogenic E. coli strain. Our work demonstrates the feasibility of modifying bacteria directly in the gut, offering a new avenue to investigate the function of bacterial genes and opening the door to the design of new microbiome-targeted therapies.

A Pseudomonas aeruginosa small RNA regulates chronic and acute infection.

The ability to switch between different lifestyles allows bacterial pathogens to thrive in diverse ecological niches1,2. However, a molecular understanding of their lifestyle changes within the human host is lacking. Here, by directly examining bacterial gene expression in human-derived samples, we discover a gene that orchestrates the transition between chronic and acute infection in the opportunistic pathogen Pseudomonas aeruginosa. The expression level of this gene, here named sicX, is the highest of the P. aeruginosa genes expressed in human chronic wound and cystic fibrosis infections, but it is expressed at extremely low levels during standard laboratory growth. We show that sicX encodes a small RNA that is strongly induced by low-oxygen conditions and post-transcriptionally regulates anaerobic ubiquinone biosynthesis. Deletion of sicX causes P. aeruginosa to switch from a chronic to an acute lifestyle in multiple mammalian models of infection. Notably, sicX is also a biomarker for this chronic-to-acute transition, as it is the most downregulated gene when a chronic infection is dispersed to cause acute septicaemia. This work solves a decades-old question regarding the molecular basis underlying the chronic-to-acute switch in P. aeruginosa and suggests oxygen as a primary environmental driver of acute lethality.

A naturally inspired antibiotic to target multidrug-resistant pathogens.

Gram-negative bacteria are responsible for an increasing number of deaths caused by antibiotic-resistant infections1,2. The bacterial natural product colistin is considered the last line of defence against a number of Gram-negative pathogens. The recent global spread of the plasmid-borne mobilized colistin-resistance gene mcr-1 (phosphoethanolamine transferase) threatens the usefulness of colistin3. Bacteria-derived antibiotics often appear in nature as collections of similar structures that are encoded by evolutionarily related biosynthetic gene clusters. This structural diversity is, at least in part, expected to be a response to the development of natural resistance, which often mechanistically mimics clinical resistance. Here we propose that a solution to mcr-1-mediated resistance might have evolved among naturally occurring colistin congeners. Bioinformatic analysis of sequenced bacterial genomes identified a biosynthetic gene cluster that was predicted to encode a structurally divergent colistin congener. Chemical synthesis of this structure produced macolacin, which is active against Gram-negative pathogens expressing mcr-1 and intrinsically resistant pathogens with chromosomally encoded phosphoethanolamine transferase genes. These Gram-negative bacteria include extensively drug-resistant Acinetobacter baumannii and intrinsically colistin-resistant Neisseria gonorrhoeae, which, owing to a lack of effective treatment options, are considered among the highest level threat pathogens4. In a mouse neutropenic infection model, a biphenyl analogue of macolacin proved to be effective against extensively drug-resistant A. baumannii with colistin-resistance, thus providing a naturally inspired and easily produced therapeutic lead for overcoming colistin-resistant pathogens.

Promoter recruitment drives the emergence of proto-genes in a long-term evolution experiment with Escherichia coli.

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli long-term evolution experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, with levels of transcription across low-expressed regions increasing in later generations of the experiment. Proto-genes formed downstream of new mutations result either from insertion element activity or chromosomal translocations that fused preexisting regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter, although such cases were rare compared to those caused by recruitment of preexisting promoters. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, can persist stably, and can serve as potential substrates for new gene formation.

Mycobacteria that cause tuberculosis have retained ancestrally acquired genes for the biosynthesis of chemically diverse terpene nucleosides.

Mycobacterium tuberculosis (Mtb) releases the unusual terpene nucleoside 1-tuberculosinyladenosine (1-TbAd) to block lysosomal function and promote survival in human macrophages. Using conventional approaches, we found that genes Rv3377c and Rv3378c, but not Rv3376, were necessary for 1-TbAd biosynthesis. Here, we introduce linear models for mass spectrometry (limms) software as a next-generation lipidomics tool to study the essential functions of lipid biosynthetic enzymes on a whole-cell basis. Using limms, whole-cell lipid profiles deepened the phenotypic landscape of comparative mass spectrometry experiments and identified a large family of approximately 100 terpene nucleoside metabolites downstream of Rv3378c. We validated the identity of previously unknown adenine-, adenosine-, and lipid-modified tuberculosinol-containing molecules using synthetic chemistry and collisional mass spectrometry, including comprehensive profiling of bacterial lipids that fragment to adenine. We tracked terpene nucleoside genotypes and lipid phenotypes among Mycobacterium tuberculosis complex (MTC) species that did or did not evolve to productively infect either human or nonhuman mammals. Although 1-TbAd biosynthesis genes were thought to be restricted to the MTC, we identified the locus in unexpected species outside the MTC. Sequence analysis of the locus showed nucleotide usage characteristic of plasmids from plant-associated bacteria, clarifying the origin and timing of horizontal gene transfer to a pre-MTC progenitor. The data demonstrated correlation between high level terpene nucleoside biosynthesis and mycobacterial competence for human infection, and 2 mechanisms of 1-TbAd biosynthesis loss. Overall, the selective gain and evolutionary retention of tuberculosinyl metabolites in modern species that cause human TB suggest a role in human TB disease, and the newly discovered molecules represent candidate disease-specific biomarkers.

A metabolomics pipeline for the mechanistic interrogation of the gut microbiome.

Gut microorganisms modulate host phenotypes and are associated with numerous health effects in humans, ranging from host responses to cancer immunotherapy to metabolic disease and obesity. However, difficulty in accurate and high-throughput functional analysis of human gut microorganisms has hindered efforts to define mechanistic connections between individual microbial strains and host phenotypes. One key way in which the gut microbiome influences host physiology is through the production of small molecules1-3, yet progress in elucidating this chemical interplay has been hindered by limited tools calibrated to detect the products of anaerobic biochemistry in the gut. Here we construct a microbiome-focused, integrated mass-spectrometry pipeline to accelerate the identification of microbiota-dependent metabolites in diverse sample types. We report the metabolic profiles of 178 gut microorganism strains using our library of 833 metabolites. Using this metabolomics resource, we establish deviations in the relationships between phylogeny and metabolism, use machine learning to discover a previously undescribed type of metabolism in Bacteroides, and reveal candidate biochemical pathways using comparative genomics. Microbiota-dependent metabolites can be detected in diverse biological fluids from gnotobiotic and conventionally colonized mice and traced back to the corresponding metabolomic profiles of cultured bacteria. Collectively, our microbiome-focused metabolomics pipeline and interactive metabolomics profile explorer are a powerful tool for characterizing microorganisms and interactions between microorganisms and their host.

Transcriptional Pausing as a Mediator of Bacterial Gene Regulation.

Cellular life depends on transcription of DNA by RNA polymerase to express genetic information. RNA polymerase has evolved not just to read information from DNA and write it to RNA but also to sense and process information from the cellular and extracellular environments. Much of this information processing occurs during transcript elongation, when transcriptional pausing enables regulatory decisions. Transcriptional pauses halt RNA polymerase in response to DNA and RNA sequences and structures at locations and times that help coordinate interactions with small molecules and transcription factors important for regulation. Four classes of transcriptional pause signals are now evident after decades of study: elemental pauses, backtrack pauses, hairpin-stabilized pauses, and regulator-stabilized pauses. In this review, I describe current understanding of the molecular mechanisms of these four classes of pause signals, remaining questions about how RNA polymerase responds to pause signals, and the many exciting directions now open to understand pausing and the regulation of transcript elongation on a genome-wide scale.

Action of a minimal contractile bactericidal nanomachine.

R-type bacteriocins are minimal contractile nanomachines that hold promise as precision antibiotics1-4. Each bactericidal complex uses a collar to bridge a hollow tube with a contractile sheath loaded in a metastable state by a baseplate scaffold1,2. Fine-tuning of such nucleic acid-free protein machines for precision medicine calls for an atomic description of the entire complex and contraction mechanism, which is not available from baseplate structures of the (DNA-containing) T4 bacteriophage5. Here we report the atomic model of the complete R2 pyocin in its pre-contraction and post-contraction states, each containing 384 subunits of 11 unique atomic models of 10 gene products. Comparison of these structures suggests the following sequence of events during pyocin contraction: tail fibres trigger lateral dissociation of baseplate triplexes; the dissociation then initiates a cascade of events leading to sheath contraction; and this contraction converts chemical energy into mechanical force to drive the iron-tipped tube across the bacterial cell surface, killing the bacterium.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25.

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

CRISPR-Cas immunity is repressed by the LysR-type transcriptional regulator PigU.

Bacteria protect themselves from infection by bacteriophages (phages) using different defence systems, such as CRISPR-Cas. Although CRISPR-Cas provides phage resistance, fitness costs are incurred, such as through autoimmunity. CRISPR-Cas regulation can optimise defence and minimise these costs. We recently developed a genome-wide functional genomics approach (SorTn-seq) for high-throughput discovery of regulators of bacterial gene expression. Here, we applied SorTn-seq to identify loci influencing expression of the two type III-A Serratia CRISPR arrays. Multiple genes affected CRISPR expression, including those involved in outer membrane and lipopolysaccharide synthesis. By comparing loci affecting type III CRISPR arrays and cas operon expression, we identified PigU (LrhA) as a repressor that co-ordinately controls both arrays and cas genes. By repressing type III-A CRISPR-Cas expression, PigU shuts off CRISPR-Cas interference against plasmids and phages. PigU also represses interference and CRISPR adaptation by the type I-F system, which is also present in Serratia. RNA sequencing demonstrated that PigU is a global regulator that controls secondary metabolite production and motility, in addition to CRISPR-Cas immunity. Increased PigU also resulted in elevated expression of three Serratia prophages, indicating their likely induction upon sensing PigU-induced cellular changes. In summary, PigU is a major regulator of CRISPR-Cas immunity in Serratia.

Long range segmentation of prokaryotic genomes by gene age and functionality.

Bacterial and archaeal genomes encompass numerous operons that typically consist of two to five genes. On larger scales, however, gene order is poorly conserved through the evolution of prokaryotes. Nevertheless, non-random localization of different classes of genes on prokaryotic chromosomes could reflect important functional and evolutionary constraints. We explored the patterns of genomic localization of evolutionarily conserved (ancient) and variable (young) genes across the diversity of bacteria and archaea. Nearly all bacterial and archaeal chromosomes were found to encompass large segments of 100-300 kb that were significantly enriched in either ancient or young genes. Similar clustering of genes with lethal knockout phenotype (essential genes) was observed as well. Mathematical modeling of genome evolution suggests that this long-range gene clustering in prokaryotic chromosomes reflects perpetual genome rearrangement driven by a combination of selective and neutral processes rather than evolutionary conservation.

Rapid functional activation of horizontally transferred eukaryotic intron-containing genes in the bacterial recipient.

Horizontal gene transfer has occurred across all domains of life and contributed substantially to the evolution of both prokaryotes and eukaryotes. Previous studies suggest that many horizontally transferred eukaryotic genes conferred selective advantages to bacterial recipients, but how these eukaryotic genes evolved into functional bacterial genes remained unclear, particularly how bacteria overcome the expressional barrier posed by eukaryotic introns. Here, we first confirmed that the presence of intron would inactivate the horizontally transferred gene in Escherichia coli even if this gene could be efficiently transcribed. Subsequent large-scale genetic screens for activation of gene function revealed that activation events could rapidly occur within several days of selective cultivation. Molecular analysis of activation events uncovered two distinct mechanisms how bacteria overcome the intron barrier: (i) intron was partially deleted and the resulting stop codon-removed mutation led to one intact foreign protein or (ii) intron was intactly retained but it mediated the translation initiation and the interaction of two split small proteins (derived from coding sequences up- and downstream of intron, respectively) to restore gene function. Our findings underscore the likelihood that horizontally transferred eukaryotic intron-containing genes could rapidly acquire functionality if they confer a selective advantage to the prokaryotic recipient.

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

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

Bacterial origin of a key innovation in the evolution of the vertebrate eye.

The vertebrate eye was described by Charles Darwin as one of the greatest potential challenges to a theory of natural selection by stepwise evolutionary processes. While numerous evolutionary transitions that led to the vertebrate eye have been explained, some aspects appear to be vertebrate specific with no obvious metazoan precursor. One critical difference between vertebrate and invertebrate vision hinges on interphotoreceptor retinoid-binding protein (IRBP, also known as retinol-binding protein, RBP3), which enables the physical separation and specialization of cells in the vertebrate visual cycle by promoting retinoid shuttling between cell types. While IRBP has been functionally described, its evolutionary origin has remained elusive. Here, we show that IRBP arose via acquisition of novel genetic material from bacteria by interdomain horizontal gene transfer (iHGT). We demonstrate that a gene encoding a bacterial peptidase was acquired prior to the radiation of extant vertebrates >500 Mya and underwent subsequent domain duplication and neofunctionalization to give rise to vertebrate IRBP. Our phylogenomic analyses on >900 high-quality genomes across the tree of life provided the resolution to distinguish contamination in genome assemblies from true instances of horizontal acquisition of IRBP and led us to discover additional independent transfers of the same bacterial peptidase gene family into distinct eukaryotic lineages. Importantly, this work illustrates the evolutionary basis of a key transition that led to the vertebrate visual cycle and highlights the striking impact that acquisition of bacterial genes has had on vertebrate evolution.