The design and analysis of transposon insertion sequencing experiments

Identification of essential genes by transposon insertion sequencing and genome-scale metabolic model construction in Streptococcus suis

Bacterial essential genes are indispensable for the survival of bacteria and therefore are attractive targets for novel anti-microbial drugs. Identifying essential genes provides a roadmap for developing novel antibiotics and anti-microbial therapies. In this study, combining high-throughput transposon sequencing (Tn-seq) and genome-scale metabolic model (GEM) construction, essential genes of Streptococcus suis, an important emerging zoonotic bacterial pathogen, were analyzed. A highly efficient transposon (Tn) mutagenesis system was developed in S. suis. This system facilitated the construction of a high-density library containing over 160,000 Tn mutants. By sequencing the library and data analysis, more than 21,000 insertion sites and 150 essential genes for growth in the rich medium were identified. Subsequently, a GEM of S. suis SC19 strain was constructed, and 165 essential genes were predicted via flux balance analysis (FBA). A total of 244 essential genes were obtained by combining the results of Tn-seq, and FBA performed. Gene identity analysis revealed 101 essential genes as potential anti-bacterial drug targets. Among them, apart from many known antibiotic targets, some interesting essential genes were also identified, including those involved in capsule biosynthesis, aminoacyl-tRNA biosynthesis, lipid biosynthesis, cell division, and cell signaling. This work identified essential genes of S. suis at the whole-genome level, providing a reference for the mining of novel anti-microbial drug targets.

IMPORTANCE

Anti-microbial resistance (AMR) presents an escalating challenge, making anti-microbial drug development an urgent need. Bacterial essential genes represent promising targets for anti-microbial drugs. However, conventional approaches to identifying bacterial essential genes are time and labor intensive. Techniques such as Tn-seq and GEM construction offer a high-throughput approach for this identification. Streptococcus suis is an emerging zoonotic bacterial pathogen, posing a big threat to public health as well as the pig industry, and the levels of AMR are increasing. Our study has successfully identified essential genes in S. suis, providing crucial insights for the discovery of new anti-microbial drug targets.

Inducible transposon mutagenesis identifies bacterial fitness determinants during infection in mice

Transposon insertion sequencing (Tn-seq) is a powerful method for genome-scale forward genetics in bacteria. However, inefficient transposon delivery or stochastic loss of mutants due to population bottlenecks can limit its effectiveness. Here we have developed 'InducTn-seq', where an arabinose-inducible Tn5 transposase enables temporal control of mini-Tn5 transposition. InducTn-seq generated up to 1.2 million transposon mutants from a single colony of enterotoxigenic Escherichia coli, Salmonella typhimurium, Shigella flexneri and Citrobacter rodentium. This mutant diversity enabled more sensitive detection of subtle fitness defects and measurement of quantitative fitness effects for essential and non-essential genes. Applying InducTn-seq to C. rodentium in a mouse model of infectious colitis bypassed a highly restrictive host bottleneck, generating a diverse population of >5 × 105 unique transposon mutants compared to 10-102 recovered by traditional Tn-seq. This in vivo screen revealed that the C. rodentium type I-E CRISPR system is required to suppress a toxin otherwise activated during gut colonization. Our findings highlight the potential of InducTn-seq for genome-scale forward genetic screens in bacteria.

Functional genomics of Campylobacter -host interactions in an intestinal tissue model reveals a small lipoprotein essential for flagellar assembly

Campylobacter jejuni is currently the most common cause of bacterial gastroenteritis worldwide. However, its genome provides few clues about how it interacts with the host. Moreover, infection screens have often been limited to classical cell culture or animal models. To identify C. jejuni genes involved in host cell interactions, we applied transposon sequencing in a humanized 3D intestinal infection model based on tissue engineering. This revealed key proteins required for host cell adherence and/or internalization, including an Rrf2 family transcriptional regulator as well as three so far uncharacterized genes ( pflC / Cj1643 , pflD / Cj0892c , pflE / Cj0978c ), which we demonstrate to encode factors essential for motility. Deletion mutants of pflC / D / E are non-motile but retain intact, paralysed flagella filaments. We demonstrate that two of these newly identified motility proteins, PflC and PflD, are components of the C. jejuni 's periplasmic disk structures of the high torque motor. The third gene, pflE , encodes a small protein of only 57 aa. Using CryoET imaging we uncovered that the small protein has a striking effect on motor biogenesis, leading to a complete loss of the flagellar disk and motor structures upon its deletion. While PflE does not appear to be a structural component of the motor itself, our data suggests that it is a lipoprotein and supports localization of the main basal disk protein FlgP, which is the first assembly step of the flagellar disk structure. Despite being annotated as a lipoprotein, we find that C. jejuni FlgP instead relies on PflE for its association with the outer membrane. Overall, our genome-wide screen revealed novel C. jejuni host interaction factors including a transcriptional regulator as well as two structural components and a small protein crucial for biogenesis of the C. jejuni high torque flagella motor. Since the flagella machinery is a critical virulence determining factor for C. jejuni , our work demonstrates how such a small protein can, quite literally, bring a bacterial pathogen to a halt.

Insights into ecology, pathogenesis, and biofilm formation of Enterococcus faecalis from functional genomics

SUMMARYEnterococcus faecalis is a significant resident of the gastrointestinal tract of most animals, including humans. Although generally non-pathogenic in healthy hosts, this microbe is adept at the exploitation of compromises in host immune functions, resulting in life-threatening opportunistic infections whose treatments are complicated by a high degree of intrinsic and acquired resistance to antimicrobial chemotherapy. Historically, progress in enterococcal research was limited by a lack of experimental models that replicate natural infection pathways and the relevance of in vitro studies to the natural biology of the organism. In this review, we summarize the history of enterococcal research during the 20th and early 21st centuries and describe more recent genetic and genomic tools and screens developed to address challenges in the field. We also describe how the results of recent studies reveal the importance of previously uncharacterized enterococcal genes, and we provide examples of interesting determinants that have emerged as important contributors to enterococcal biology. These factors may also serve as targets for future vaccines and chemotherapeutic agents to combat life-threatening hospital infections.

Genome-wide mutagenesis identifies factors involved in MRSA vaginal colonization

Methicillin-resistant Staphylococcus aureus (MRSA) is an opportunistic pathogen that colonizes various body sites, including the nares, skin, and vagina. During pregnancy,colonization can lead to dysbiosis, adverse pregnancy outcomes, and invasive disease. To identify genes contributing to MRSA vaginal fitness, we performed transposon sequencing (Tn-seq) using a murine model of vaginal colonization, identifying over 250 conditionally essential genes. Five genes were validated in our murine model, including those encoding the aerobic respiration protein QoxB, bacillithiol biosynthesis component BshB2, sialic acid catabolism enzyme NanE, and staphylococcal regulator of respiration SrrAB. RNA sequencing and comparative analysis identified over 30 SrrAB-regulated genes potentially important for fitness in vaginal-like conditions, particularly under oxygen stress. These findings highlight pathways such as aerobic respiration, bacillithiol biosynthesis, sialic acid catabolism, and transcriptional regulation that support MRSA's competitive fitness in the vaginal tract.

Context matters: assessing the impacts of genomic background and ecology on microbial biosynthetic gene cluster evolution

Encoded within many microbial genomes, biosynthetic gene clusters (BGCs) underlie the synthesis of various secondary metabolites that often mediate ecologically important functions. Several studies and bioinformatics methods developed over the past decade have advanced our understanding of both microbial pangenomes and BGC evolution. In this minireview, we first highlight challenges in broad evolutionary analysis of BGCs, including delineation of BGC boundaries and clustering of BGCs across genomes. We further summarize key findings from microbial comparative genomics studies on BGC conservation across taxa and habitats and discuss the potential fitness effects of BGCs in different settings. Afterward, recent research showing the importance of genomic context on the production of secondary metabolites and the evolution of BGCs is highlighted. These studies draw parallels to recent, broader, investigations on gene-to-gene associations within microbial pangenomes. Finally, we describe mechanisms by which microbial pangenomes and BGCs evolve, ranging from the acquisition or origination of entire BGCs to micro-evolutionary trends of individual biosynthetic genes. An outlook on how expansions in the biosynthetic capabilities of some taxa might support theories that open pangenomes are the result of adaptive evolution is also discussed. We conclude with remarks about how future work leveraging longitudinal metagenomics across diverse ecosystems is likely to significantly improve our understanding on the evolution of microbial genomes and BGCs.

ArgR-dependent bacterial resistance to host lipid droplets in Edwardsiella piscicida

Lipid droplets (LDs), as innate immune hubs, function in the front line of antimicrobial defense involved in the host-pathogen arms race. Particularly for intracellular bacterial pathogens, the endowed capacity to resist host LDs can effectively facilitate pathogen in vivo colonization and evasion from the host's innate immune response. Here, to investigate the genetic mechanisms of intracellular bacteria response to host LDs, we utilized transposon insertion sequencing to dissect critical fitness determinants of Edwardsiella piscicida under the treatment of LDs isolated from its native host, turbot. Targeted metabolomics indicated that LD challenge resulted in the accumulation of intracellular arginine. The core arginine metabolism regulatory factor, ArgR, was found to play a pivotal role in combating LDs, emphasizing the importance of orchestrating intracellular arginine levels for bacterial LD adaptation. Specifically, ArgR enhanced the expressions of genes involved in arginine catabolism (speA/B and arcC) and diminished gene transcripts associated with arginine import (artP) and synthesis (argD/E/H). Furthermore, ArgR contributed to the pathogenesis of E. piscicida, promoting the proliferation in host cells and virulence in turbot. Collectively, our results shed light on the underlying mechanism of intracellular pathogens resisting LDs during bacterial infections and highlighting the crucial role of arginine in the host-pathogen interactions.

Saturated Transposon Analysis in Yeast as a one-step method to quantify the fitness effects of gene disruptions on a genome-wide scale

Transposon insertion site sequencing (TIS) is a powerful tool that has significantly advanced our knowledge of functional genomics. For example, TIS has been used to identify essential genes of Saccharomyces cerevisiae, screen for antibiotic resistance genes in Klebsiella pneumoniae and determine the set of genes required for virulence of Mycobacterium tuberculosis. While providing valuable insights, these applications of TIS focus on (conditional) gene essentiality and neglect possibly interesting but subtle differences in the importance of genes for fitness. Notably, it has been demonstrated that data obtained from TIS experiments can be used for fitness quantification and the construction of genetic interaction maps, but this potential is only sporadically exploited. Here, we present a method to quantify the fitness of gene disruption mutants using data obtained from a TIS screen developed for the yeast Saccharomyces cerevisiae called SATAY. We show that the mean read count per transposon insertion site provides a metric for fitness that is robust across biological and technical replicate experiments. Importantly, the ability to resolve differences between gene disruption mutants with low fitness depends crucially on the inclusion of insertion sites that are not observed in the sequencing data to estimate the mean. While our method provides reproducible results between replicate SATAY datasets, the obtained fitness distribution differs substantially from those obtained using other techniques. It is currently unclear whether these inconsistencies are due to biological or technical differences between the methods. We end with suggestions for modifications of the SATAY procedure that could improve the resolution of the fitness estimates. Our analysis indicates that increasing the sequencing depth does very little to reduce the uncertainty in the estimates, while replacing the PCR amplification with methods that avoid or reduce the number of amplification cycles will likely be most effective in reducing noise.

Discovering genetic modulators of the protein homeostasis system through multilevel analysis

Every protein progresses through a natural lifecycle from birth to maturation to death; this process is coordinated by the protein homeostasis system. Environmental or physiological conditions trigger pathways that maintain the homeostasis of the proteome. An open question is how these pathways are modulated to respond to the many stresses that an organism encounters during its lifetime. To address this question, we tested how the fitness landscape changes in response to environmental and genetic perturbations using directed and massively parallel transposon mutagenesis in Caulobacter crescentus. We developed a general computational pipeline for the analysis of gene-by-environment interactions in transposon mutagenesis experiments. This pipeline uses a combination of general linear models, statistical knockoffs, and a nonparametric Bayesian statistical model to identify essential genetic network components that are shared across environmental perturbations. This analysis allows us to quantify the similarity of proteotoxic environmental perturbations from the perspective of the fitness landscape. We find that essential genes vary more by genetic background than by environmental conditions, with limited overlap among mutant strains targeting different facets of the protein homeostasis system. We also identified 146 unique fitness determinants across different strains, with 19 genes common to at least two strains, showing varying resilience to proteotoxic stresses. Experiments exposing cells to a combination of genetic perturbations and dual environmental stressors show that perturbations that are quantitatively dissimilar from the perspective of the fitness landscape are likely to have a synergistic effect on the growth defect.

Transposon mutagenesis screen in Klebsiella pneumoniae identifies genetic determinants required for growth in human urine and serum

Klebsiella pneumoniae is a global public health concern due to the rising myriad of hypervirulent and multidrug-resistant clones both alarmingly associated with high mortality. The molecular mechanisms underpinning these recalcitrant K. pneumoniae infection, and how virulence is coupled with the emergence of lineages resistant to nearly all present-day clinically important antimicrobials, are unclear. In this study, we performed a genome-wide screen in K. pneumoniae ECL8, a member of the endemic K2-ST375 pathotype most often reported in Asia, to define genes essential for growth in a nutrient-rich laboratory medium (Luria-Bertani [LB] medium), human urine, and serum. Through transposon directed insertion-site sequencing (TraDIS), a total of 427 genes were identified as essential for growth on LB agar, whereas transposon insertions in 11 and 144 genes decreased fitness for growth in either urine or serum, respectively. These studies not only provide further knowledge on the genetics of this pathogen but also provide a strong impetus for discovering new antimicrobial targets to improve current therapeutic options for K. pneumoniae infections.

The uncharacterized PA3040-3042 operon is part of the cell envelope stress response and a tobramycin resistance determinant in a clinical isolate of Pseudomonas aeruginosa

Bacteriophages (hereafter "phages") are ubiquitous predators of bacteria in the natural world, but interest is growing in their development into antibacterial therapy as complement or replacement for antibiotics. However, bacteria have evolved a huge variety of antiphage defense systems allowing them to resist phage lysis to a greater or lesser extent. In addition to dedicated phage defense systems, some aspects of the general stress response also impact phage susceptibility, but the details of this are not well known. In order to elucidate these factors in the opportunistic pathogen Pseudomonas aeruginosa, we used the laboratory-conditioned strain PAO1 as host for phage infection experiments as it is naturally poor in dedicated phage defense systems. Screening by transposon insertion sequencing indicated that the uncharacterized operon PA3040-PA3042 was potentially associated with resistance to lytic phages. However, we found that its primary role appeared to be in regulating biofilm formation, particularly in a clinical isolate of P. aeruginosa in which it also altered tobramycin resistance. Its expression was highly growth-phase dependent and responsive to phage infection and cell envelope stress. Our results suggest that this operon may be a cryptic but important locus for P. aeruginosa stress tolerance.

IMPORTANCE

An important category of bacterial stress response systems is bacteriophage defense, where systems are triggered by bacteriophage infection and activate a response which may either destroy the phage genome or destroy the infected cell so that the rest of the population survives. In some bacteria, the cell envelope stress response is activated by bacteriophage infection, but it is unknown whether this contributes to the survival of the infection. We have found that a conserved uncharacterized operon (PA3040-PA3042) of the cell envelope stress regulon in Pseudomonas aeruginosa, which has very few dedicated phage defense systems, responds to phage infection and stationary phase as well as envelope stress and is important for growth and biofilm formation in a clinical isolate of P. aeruginosa, even in the absence of phages. As homologs of these genes are found in other bacteria, they may be a novel component of the general stress response.

Shedding light on bacteria-host interactions with the aid of TnSeq approaches

Bacteria are highly adaptable and grow in diverse niches, where they often interact with eukaryotic organisms. These interactions with different hosts span the entire spectrum from symbiosis to pathogenicity and thus determine the lifestyle of the bacterium. Knowledge of the genetic determinants involved in animal and plant host colonization by pathogenic and mutualistic bacteria is not only crucial to discover new drug targets for disease management but also for developing novel biostimulant strategies. In the last decades, significant progress in genome-wide high-throughput technologies such as transposon insertion sequencing has led to the identification of pathways that enable efficient host colonization. However, the extent to which similar genes play a role in this process in different bacteria is yet unclear. This review highlights the commonalities and specificities of bacterial determinants important for bacteria-host interaction.

Inducible transposon mutagenesis for genome-scale forward genetics

Transposon insertion sequencing (Tn-seq) is a powerful method for genome-scale functional genetics in bacteria. However, its effectiveness is often limited by a lack of mutant diversity, caused by either inefficient transposon delivery or stochastic loss of mutants due to population bottlenecks. Here, we introduce "InducTn-seq", which leverages inducible mutagenesis for temporal control of transposition. InducTn-seq generates millions of transposon mutants from a single colony, enabling the sensitive detection of subtle fitness defects and transforming binary classifications of gene essentiality into a quantitative fitness measurement across both essential and non-essential genes. Using a mouse model of infectious colitis, we show that InducTn-seq bypasses a highly restrictive host bottleneck to generate a diverse transposon mutant population from the few cells that initiate infection, revealing the role of oxygen-related metabolic plasticity in pathogenesis. Overall, InducTn-seq overcomes the limitations of traditional Tn-seq, unlocking new possibilities for genome-scale forward genetic screens in bacteria.

Biochemical and structural characterization reveals Rv3400 codes for β-phosphoglucomutase in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) adapt to various host environments and utilize a variety of sugars and lipids as carbon sources. Among these sugars, maltose and trehalose, also play crucial role in bacterial physiology and virulence. However, some key enzymes involved in trehalose and maltose metabolism in Mtb are not yet known. Here we structurally and functionally characterized a conserved hypothetical gene Rv3400. We determined the crystal structure of Rv3400 at 1.7 Å resolution. The crystal structure revealed that Rv3400 adopts Rossmann fold and shares high structural similarity with haloacid dehalogenase family of proteins. Our comparative structural analysis suggested that Rv3400 could perform either phosphatase or pyrophosphatase or β-phosphoglucomutase (β-PGM) activity. Using biochemical studies, we further confirmed that Rv3400 performs β-PGM activity and hence, Rv3400 encodes for β-PGM in Mtb. Our data also confirm that Mtb β-PGM is a metal dependent enzyme having broad specificity for divalent metal ions. β-PGM converts β-D-glucose-1-phosphate to β-D-glucose-6-phosphate which is required for the generation of ATP and NADPH through glycolysis and pentose phosphate pathway, respectively. Using site directed mutagenesis followed by biochemical studies, we show that two Asp residues in the highly conserved DxD motif, D29 and D31, are crucial for enzyme activity. While D29A, D31A, D29E, D31E and D29N mutants lost complete activity, D31N mutant retained about 30% activity. This study further helps in understanding the role of β-PGM in the physiology of Mtb.

Genome-wide analysis of acid tolerance genes of Enterococcus faecalis with RNA-seq and Tn-seq

Enterococcus faecalis, a formidable nosocomial and community-acquired opportunistic pathogen, can persist a wide range of extreme environments, including low pH and nutrient deficiency. Clarifying the survival mechanism of E. faecalis in low-pH conditions is the key to combating the infectious diseases caused by E. faecalis. In this study, we combined transcriptome profiling (RNA-seq) and transposon insertion sequencing (TIS) to comprehensively understand the genes that confer these features on E. faecalis. The metadata showed that genes whose products are involved in cation transportation and amino acid biosynthesis were predominantly differentially expressed under acid conditions. The products of genes such as opp1C and copY reduced the hydrion concentration in the cell, whereas those of gldA2, gnd2, ubiD, and ubiD2 mainly participated in amino metabolism, increasing matters to neutralize excess acid. These, together with the folE and hexB genes, which are involved in mismatch repair, form a network of E. faecalis genes necessary for its survival under acid conditions.

Discovering Genetic Modulators of the Protein Homeostasis System through Multilevel Analysis

Every protein progresses through a natural lifecycle from birth to maturation to death; this process is coordinated by the protein homeostasis system. Environmental or physiological conditions trigger pathways that maintain the homeostasis of the proteome. An open question is how these pathways are modulated to respond to the many stresses that an organism encounters during its lifetime. To address this question, we tested how the fitness landscape changes in response to environmental and genetic perturbations using directed and massively parallel transposon mutagenesis in Caulobacter crescentus. We developed a general computational pipeline for the analysis of gene-by-environment interactions in transposon mutagenesis experiments. This pipeline uses a combination of general linear models (GLMs), statistical knockoffs, and a nonparametric Bayesian statistical model to identify essential genetic network components that are shared across environmental perturbations. This analysis allows us to quantify the similarity of proteotoxic environmental perturbations from the perspective of the fitness landscape. We find that essential genes vary more by genetic background than by environmental conditions, with limited overlap among mutant strains targeting different facets of the protein homeostasis system. We also identified 146 unique fitness determinants across different strains, with 19 genes common to at least two strains, showing varying resilience to proteotoxic stresses. Experiments exposing cells to a combination of genetic perturbations and dual environmental stressors show that perturbations that are quantitatively dissimilar from the perspective of the fitness landscape are likely to have a synergistic effect on the growth defect.

CRISPR-based screening of small RNA modulators of bile susceptibility in Bacteroides thetaiotaomicron

Microbiota-centric interventions are limited by our incomplete understanding of the gene functions of many of its constituent species. This applies in particular to small RNAs (sRNAs), which are emerging as important regulators in microbiota species yet tend to be missed by traditional functional genomics approaches. Here, we establish CRISPR interference (CRISPRi) in the abundant microbiota member Bacteroides thetaiotaomicron for genome-wide sRNA screens. By assessing the abundance of different protospacer-adjacent motifs, we identify the Prevotella bryantii B14 Cas12a as a suitable nuclease for CRISPR screens in these bacteria and generate an inducible Cas12a expression system. Using a luciferase reporter strain, we infer guide design rules and use this knowledge to assemble a computational pipeline for automated gRNA design. By subjecting the resulting guide library to a phenotypic screen, we uncover the sRNA BatR to increase susceptibility to bile salts through the regulation of genes involved in Bacteroides cell surface structure. Our study lays the groundwork for unlocking the genetic potential of these major human gut mutualists and, more generally, for identifying hidden functions of bacterial sRNAs.

A peptidoglycan N-deacetylase specific for anhydroMurNAc chain termini in Agrobacterium tumefaciens

During growth, bacteria remodel and recycle their peptidoglycan (PG). A key family of PG-degrading enzymes is the lytic transglycosylases, which produce anhydromuropeptides, a modification that caps the PG chains and contributes to bacterial virulence. Previously, it was reported that the polar-growing Gram-negative plant pathogen Agrobacterium tumefaciens lacks anhydromuropeptides. Here, we report the identification of an enzyme, MdaA (MurNAc deacetylase A), which specifically removes the acetyl group from anhydromuropeptide chain termini in A. tumefaciens, resolving this apparent anomaly. A. tumefaciens lacking MdaA accumulates canonical anhydromuropeptides, whereas MdaA was able to deacetylate anhydro-N-acetyl muramic acid in purified sacculi that lack this modification. As for other PG deacetylases, MdaA belongs to the CE4 family of carbohydrate esterases but harbors an unusual Cys residue in its active site. MdaA is conserved in other polar-growing bacteria, suggesting a possible link between PG chain terminus deacetylation and polar growth.

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

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

CRISPR Interference-Based Functional Small RNA Genomics

Small RNAs (sRNAs) are versatile regulators universally present in species across the prokaryotic kingdom, yet their functional characterization remains a major bottleneck. Gene inactivation through random transposon insertion has proven extremely valuable in discovering hidden gene functions. However, this approach is biased toward long genes and usually results in the underrepresentation of sRNA mutants. In contrast, CRISPR interference (CRISPRi) harnesses guide RNAs to recruit cleavage-deficient Cas nucleases to specific DNA loci. The ensuing steric hindrance inhibits RNA polymerase assembly at-or migration along-predefined genes, allowing for targeted knockdown screens without major length bias. In this chapter, we provide a detailed protocol for CRISPRi-based functional screening of bacterial sRNAs. Using the abundant microbiota species Bacteroides thetaiotaomicron as a model, we describe the design and generation of a guide library targeting the full intergenic sRNA repertoire of this organism and its application to identify sRNA knockdown-associated fitness effects. Our protocol is generic and thus suitable for the systematic assessment of sRNA-associated phenotypes in a wide range of bacterial species and experimental conditions. We expect CRISPRi-based functional genomics to boost sRNA research in understudied bacterial taxa, for instance, members of the gut microbiota.

Recent Advances and Techniques for Identifying Novel Antibacterial Targets

BACKGROUND

With the emergence of drug-resistant bacteria, the development of new antibiotics is urgently required. Target-based drug discovery is the most frequently employed approach for the drug development process. However, traditional drug target identification techniques are costly and time-consuming. As research continues, innovative approaches for antibacterial target identification have been developed which enabled us to discover drug targets more easily and quickly.

METHODS

In this review, methods for finding drug targets from omics databases have been discussed in detail including principles, procedures, advantages, and potential limitations. The role of phage-driven and bacterial cytological profiling approaches is also discussed. Moreover, current article demonstrates the advancements being made in the establishment of computational tools, machine learning algorithms, and databases for antibacterial target identification.

RESULTS

Bacterial drug targets successfully identified by employing these aforementioned techniques are described as well.

CONCLUSION

The goal of this review is to attract the interest of synthetic chemists, biologists, and computational researchers to discuss and improve these methods for easier and quicker development of new drugs.

Analysis of key genes for the survival of Pantoea agglomerans under nutritional stress

The absolute amount of nutrients on plant leaves is usually low, and the growth of epiphytic bacteria is typically limited by nutrient content. Thus, is of great significance to study the survival mechanism of epiphytes under nutritional stress for plant disease control. In this paper, Pantoea agglomerans CHTF15 isolated from walnut leaves was used to detect the key genes for the survival of the bacterium under simulated nutrient stress in artificial medium. Genome sequencing was combined with transposon insertion sequencing (Tn-seq) for the detection technique. A total of 105 essential genes were screened from the whole genome. The genes were mainly related to the nucleotide metabolism, protein metabolism, biological oxidation and the gene repair of bacteria analyzed by gene ontology (GO) enrichment analysis. Volcano map analysis demonstrated that the functions of the 15 genes with the most significant differences were generally related to the synthesis of amino acids or proteins, the nucleotide metabolism and homologous recombination and repair. Competitive index analysis revealed that the deletion of the genes dksA and epmA regulating protein synthesis, the gene ribB involved in the nucleotide metabolism and the gene xerD involved in recombination repair induced a significant reduction in the survival ability of the corresponding mutants in the 0.10 % YEP medium and the walnut leaf surface. The results act as a foundation for further in-depth research on the infection process and the mechanisms of pathogenic bacteria.

Enhanced transformation efficiency in Treponema denticola enabled by SyngenicDNA-based plasmids lacking restriction-modification target motifs

Oral spirochetes are among a small group of keystone pathogens contributing to dysregulation of tissue homeostatic processes that leads to breakdown of the tissue and bone supporting the teeth in periodontal disease. Additionally, our group has recently demonstrated that Treponema are among the dominant microbial genera detected intracellularly in tumor specimens from patients with oral squamous cell carcinoma. While over 60 species and phylotypes of oral Treponema have been detected, T. denticola is one of the few that can be grown in culture and the only one in which genetic manipulation is regularly performed. Thus, T. denticola is a key model organism for studying spirochete metabolic processes, interactions with other microbes, and host cell and tissue responses relevant to oral diseases, as well as venereal and nonvenereal treponematoses whose agents lack workable genetic systems. We previously demonstrated improved transformation efficiency using an Escherichia coli-T. denticola shuttle plasmid and its utility for expression in T. denticola of an exogenous fluorescent protein that is active under anaerobic conditions. Here, we expand on this work by characterizing T. denticola Type I and Type II restriction-modification (R-M) systems and designing a high-efficiency R-M-silent "SyngenicDNA" shuttle plasmid resistant to all T. denticola ATCC 35405 R-M systems. Resequencing of the ATCC 33520 genome revealed an additional Type I R-M system consistent with the relatively low transformation efficiency of the shuttle plasmid in this strain. Using SyngenicDNA approaches, we optimized shuttle plasmid transformation efficiency in T. denticola and used it to complement a defined T. denticola ΔfhbB mutant strain. We further report the first high-efficiency transposon mutagenesis of T. denticola using an R-M-silent, codon-optimized, himarC9 transposase-based plasmid. Thus, use of SyngenicDNA-based strategies and tools can enable further mechanistic examinations of T. denticola physiology and behavior.

Genome-wide identification of genes critical for in vivo fitness of multi-drug resistant porcine extraintestinal pathogenic Escherichia coli by transposon-directed insertion site sequencing using a mouse infection model

Extraintestinal pathogenic Escherichia coli (ExPEC) is an important zoonotic pathogen. Recently, ExPEC has been reported to be an emerging problem in pig farming. However, the mechanism of pathogenicity of porcine ExPEC remains to be revealed. In this study, we constructed a transposon (Tn) mutagenesis library covering Tn insertion in over 72% of the chromosome-encoded genes of a virulent and multi-drug resistant porcine ExPEC strain PCN033. By using a mouse infection model, a transposon-directed insertion site sequencing (TraDIS) assay was performed to identify in vivo fitness factors. By comparing the Tn insertion frequencies between the input Tn library and the recovered library from different organs, 64 genes were identified to be involved in fitness during systemic infection. 15 genes were selected and individual gene deletion mutants were constructed. The in vivo fitness was evaluated by using a competitive infection assay. Among them, ΔfimG was significantly outcompeted by the WT strain in vivo and showed defective adhesion to host cells. rfa which was involved in lipopolysaccharide biosynthesis was shown to be critical for in vivo fitness which may have resulted from its role in the resistance to serum killing. In addition, several metabolic genes including fepB, sdhC, fepG, gltS, dcuA, ccmH, ddpD, narU, glpD, malM, and yabL and two regulatory genes metJ and baeS were shown as important determinants of in vivo fitness of porcine ExPEC. Collectively, this study performed a genome-wide screening for in vivo fitness factors which will be important for understanding the pathogenicity of porcine ExPEC.

Microbial Primer: Transposon directed insertion site sequencing (TraDIS): A high throughput method for linking genotype to phenotype

Genetic screens are a key tool for linking phenotype and genotype. Transposon mutagenesis was one of the first genetic methodologies to associate genetic loci with phenotypes. The advent of next-generation sequencing transformed the use of this technique allowing rapid interrogation of whole genomes for genes that correlate with phenotype. One method is transposon directed insertion-site sequencing (TraDIS). Here we describe the method, recent developments in technology, and the advantages and disadvantages of this method compared to other genetic screening tools.

acCRISPR: an activity-correction method for improving the accuracy of CRISPR screens

High throughput CRISPR screens are revolutionizing the way scientists unravel the genetic underpinnings of engineered and evolved phenotypes. One of the critical challenges in accurately assessing screening outcomes is accounting for the variability in sgRNA cutting efficiency. Poorly active guides targeting genes essential to screening conditions obscure the growth defects that are expected from disrupting them. Here, we develop acCRISPR, an end-to-end pipeline that identifies essential genes in pooled CRISPR screens using sgRNA read counts obtained from next-generation sequencing. acCRISPR uses experimentally determined cutting efficiencies for each guide in the library to provide an activity correction to the screening outcomes via calculation of an optimization metric, thus determining the fitness effect of disrupted genes. CRISPR-Cas9 and -Cas12a screens were carried out in the non-conventional oleaginous yeast Yarrowia lipolytica and acCRISPR was used to determine a high-confidence set of essential genes for growth under glucose, a common carbon source used for the industrial production of oleochemicals. acCRISPR was also used in screens quantifying relative cellular fitness under high salt conditions to identify genes that were related to salt tolerance. Collectively, this work presents an experimental-computational framework for CRISPR-based functional genomics studies that may be expanded to other non-conventional organisms of interest.

Engineering bacteriophages for enhanced host range and efficacy: insights from bacteriophage-bacteria interactions

Bacteriophages, the most abundant organisms on earth, have the potential to address the rise of multidrug-resistant bacteria resulting from the overuse of antibiotics. However, their high specificity and limited host range can hinder their effectiveness. Phage engineering, through the use of gene editing techniques, offers a means to enhance the host range of bacteria, improve phage efficacy, and facilitate efficient cell-free production of phage drugs. To engineer phages effectively, it is necessary to understand the interaction between phages and host bacteria. Understanding the interaction between the receptor recognition protein of bacteriophages and host receptors can serve as a valuable guide for modifying or replacing these proteins, thereby altering the receptor range of the bacteriophage. Research and development focused on the CRISPR-Cas bacterial immune system against bacteriophage nucleic acids can provide the necessary tools to promote recombination and counter-selection in engineered bacteriophage programs. Additionally, studying the transcription and assembly functions of bacteriophages in host bacteria can facilitate the engineered assembly of bacteriophage genomes in non-host environments. This review highlights a comprehensive summary of phage engineering methods, including in-host and out-of-host engineering, and the use of high-throughput methods to understand their role. The main aim of these techniques is to harness the intricate interactions between bacteriophages and hosts to inform and guide the engineering of bacteriophages, particularly in the context of studying and manipulating the host range of bacteriophages. By employing advanced high-throughput methods to identify specific bacteriophage receptor recognition genes, and subsequently introducing modifications or performing gene swapping through in-host recombination or out-of-host synthesis, it becomes possible to strategically alter the host range of bacteriophages. This capability holds immense significance for leveraging bacteriophages as a promising therapeutic approach against antibiotic-resistant bacteria.

Functional annotation and importance of marine bacterial transporters of plankton exometabolites

Metabolite exchange within marine microbial communities transfers carbon and other major elements through global cycles and forms the basis of microbial interactions. Yet lack of gene annotations and concern about the quality of existing ones remain major impediments to revealing currencies of carbon flux. We employed an arrayed mutant library of the marine bacterium Ruegeria pomeroyi DSS-3 to experimentally annotate substrates of organic compound transporter systems, using mutant growth and compound drawdown analyses to link transporters to their cognate substrates. Mutant experiments verified substrates for thirteen R. pomeroyi transporters. Four were previously hypothesized based on gene expression data (taurine, glucose/xylose, isethionate, and cadaverine/putrescine/spermidine); five were previously hypothesized based on homology to experimentally annotated transporters in other bacteria (citrate, glycerol, N-acetylglucosamine, fumarate/malate/succinate, and dimethylsulfoniopropionate); and four had no previous annotations (thymidine, carnitine, cysteate, and 3-hydroxybutyrate). These bring the total number of experimentally-verified organic carbon influx transporters to 18 of 126 in the R. pomeroyi genome. In a longitudinal study of a coastal phytoplankton bloom, expression patterns of the experimentally annotated transporters linked them to different stages of the bloom, and also led to the hypothesis that citrate and 3-hydroxybutyrate were among the most highly available bacterial substrates. Improved functional annotation of the gatekeepers of organic carbon uptake is critical for deciphering carbon flux and fate in microbial ecosystems.

Identification of determinants for entering into a viable but nonculturable state in Vibrio alginolyticus by Tn-seq

The viable but nonculturable (VBNC) state is a dormant state of nonsporulating bacteria that enhances survival in adverse environments. Systematic genome-wide research on the genetic basis of VBNC formation is warranted. In this study, we demonstrated that the marine bacterium Vibrio alginolyticus lost culturability but remained viable and entered into the VBNC state when exposed to low nutrient concentrations for prolonged periods of time. Using transposon-insertion sequencing (Tn-seq), we identified 635 determinants governing the formation of the VBNC state, including 322 genes with defective effects on VBNC formation and 313 genes contributing to entry into the VBNC state. Tn-seq analysis revealed that genes involved in various metabolic pathways were shown to have an inhibitory effect on VBNC formation, while genes related to chemotaxis or folate biosynthesis promoted entry into the VBNC state. Moreover, the effects of these genes on the formation of VBNC were validated with the growth of deletion mutants of eight selected genes under nutrient-limited conditions. Interestingly, fleQ and pyrI were identified as essential for entry into the VBNC state, and they affected the formation of the VBNC state independent of RpoE or ToxR regulation. Collectively, these results provide new insights into the mechanism of VBNC formation. KEY POINTS: • Vibrio alginolyticus has the ability to enter into the VBNC state under low nutrient conditions at low temperature. • The 635 determinants for entry into the VBNC state were systematically identified by transposon-insertion sequencing. • PyrI and FleQ were validated to play significant roles in the formation of the VBNC state.

The SMC-like RecN protein is at the crossroads of several genotoxic stress responses in Escherichia coli

Introduction

DNA damage repair (DDR) is an essential process for living organisms and contributes to genome maintenance and evolution. DDR involves different pathways including Homologous recombination (HR), Nucleotide Excision Repair (NER) and Base excision repair (BER) for example. The activity of each pathway is revealed with particular drug inducing lesions, but the repair of most DNA lesions depends on concomitant or subsequent action of the multiple pathways.

Methods

In the present study, we used two genotoxic antibiotics, mitomycin C (MMC) and Bleomycin (BLM), to decipher the interplays between these different pathways in E. coli. We combined genomic methods (TIS and Hi-SC2) and imaging assays with genetic dissections.

Results

We demonstrate that only a small set of DDR proteins are common to the repair of the lesions induced by these two drugs. Among them, RecN, an SMC-like protein, plays an important role by controlling sister chromatids dynamics and genome morphology at different steps of the repair processes. We further demonstrate that RecN influence on sister chromatids dynamics is not equivalent during the processing of the lesions induced by the two drugs. We observed that RecN activity and stability requires a pre-processing of the MMC-induced lesions by the NER but not for BLM-induced lesions.

Discussion

Those results show that RecN plays a major role in rescuing toxic intermediates generated by the BER pathway in addition to its well-known importance to the repair of double strand breaks by HR.

Transposon mutagenesis in oral streptococcus

Oral streptococci are gram-positive facultative anaerobic bacteria that are normal inhabitants of the human oral cavity and play an important role in maintaining oral microecological balance and pathogenesis. Transposon mutagenesis is an effective genetic manipulation strategy for studying the function of genomic features. In order to study cariogenic related genes and crucial biological element genes of oral Streptococcus, transposon mutagenesis was widely used to identify functional genes. With the advent of next-generation sequencing (NGS) technology and the development of transposon random mutation library construction methods, transposon insertion sequencing (TIS) came into being. Benefiting from high-throughput advances in NGS, TIS was able to evaluate the fitness contribution and essentiality of genetic features in the bacterial genome. The application of transposon mutagenesis, including TIS, to oral streptococci provided a massive amount of valuable detailed linkage data between genetic fitness and genetic backgrounds, further clarify the processes of colonization, virulence, and persistence and provides a more reliable basis for investigating relationships with host ecology and disease status. This review focuses on transposon mutagenesis, including TIS, and its applicability in oral streptococci.

Genome-wide transposon mutagenesis analysis of Burkholderia pseudomallei reveals essential genes for in vitro and in vivo survival

Introduction

Burkholderia pseudomallei, a soil-dwelling microbe that infects humans and animals is the cause of the fatal disease melioidosis. The molecular mechanisms that underlie B. pseudomallei's versatility to survive within a broad range of environments are still not well defined.

Methods

We used the genome-wide screening tool TraDIS (Transposon Directed Insertion-site Sequencing) to identify B. pseudomallei essential genes. Transposon-flanking regions were sequenced and gene essentiality was assessed based on the frequency of transposon insertions within each gene. Transposon mutants were grown in LB and M9 minimal medium to determine conditionally essential genes required for growth under laboratory conditions. The Caenorhabditis elegans infection model was used to assess genes associated with in vivo B. pseudomallei survival. Transposon mutants were fed to the worms, recovered from worm intestines, and sequenced. Two selected mutants were constructed and evaluated for the bacteria's ability to survive and proliferate in the nematode intestinal lumen.

Results

Approximately 500,000 transposon-insertion mutants of B. pseudomallei strain R15 were generated. A total of 848,811 unique transposon insertion sites were identified in the B. pseudomallei R15 genome and 492 genes carrying low insertion frequencies were predicted to be essential. A total of 96 genes specifically required to support growth under nutrient-depleted conditions were identified. Genes most likely to be involved in B. pseudomallei survival and adaptation in the C. elegans intestinal lumen, were identified. When compared to wild type B. pseudomallei, a Tn5 mutant of bpsl2988 exhibited reduced survival in the worm intestine, was attenuated in C. elegans killing and showed decreased colonization in the organs of infected mice.

Discussion

The B. pseudomallei conditional essential proteins should provide further insights into the bacteria's niche adaptation, pathogenesis, and virulence.

Elucidation of Essential Genes and Mutant Fitness during Adaptation toward Nitrogen Fixation Conditions in the Endophyte Azoarcus olearius BH72 Revealed by Tn-Seq

Azoarcus olearius BH72 is a diazotrophic model endophyte that contributes fixed nitrogen to its host plant, Kallar grass, and expresses nitrogenase genes endophytically. Despite extensive studies on biological nitrogen fixation (BNF) of diazotrophic endophytes, little is known about global genetic players involved in survival under respective physiological conditions. Here, we report a global genomic screen for putatively essential genes of A. olearius employing Tn5 transposon mutagenesis with a modified transposon combined with high-throughput sequencing (Tn-Seq). A large Tn5 master library of ~6 × 105 insertion mutants of strain BH72 was obtained. Next-generation sequencing identified 183,437 unique insertion sites into the 4,376,040-bp genome, displaying one insertion every 24 bp on average. Applying stringent criteria, we describe 616 genes as putatively essential for growth on rich medium. COG (Clusters of Orthologous Groups) assignment of the 564 identified protein-coding genes revealed enrichment of genes related to core cellular functions and cell viability. To mimic gradual adaptations toward BNF conditions, the Tn5 mutant library was grown aerobically in synthetic medium or microaerobically on either combined or atmospheric nitrogen. Enrichment and depletion analysis of Tn5 mutants not only demonstrated the role of BNF- and metabolism-related proteins but also revealed that, strikingly, many genes relevant for plant-microbe interactions decrease bacterial competitiveness in pure culture, such type IV pilus- and bacterial envelope-associated genes. IMPORTANCE A constantly growing world population and the daunting challenge of climate change demand new strategies in agricultural crop production. Intensive usage of chemical fertilizers, overloading the world's fields with organic input, threaten terrestrial and marine ecosystems as well as human health. Long overlooked, the beneficial interaction of endophytic bacteria and grasses has attracted ever-growing interest in research in the last decade. Capable of biological nitrogen fixation, diazotrophic endophytes not only provide a valuable source of combined nitrogen but also are known for diverse plant growth-promoting effects, thereby contributing to plant productivity. Elucidation of an essential gene set for a prominent model endophyte such as A. olearius BH72 provides us with powerful insights into its basic lifestyle. Knowledge about genes detrimental or advantageous under defined physiological conditions may point out a way of manipulating key steps in the bacterium's lifestyle and plant interaction toward a more sustainable agriculture.

How It All Begins: Bacterial Factors Mediating the Colonization of Invertebrate Hosts by Beneficial Symbionts

Beneficial associations with bacteria are widespread across animals, spanning a range of symbiont localizations, transmission routes, and functions. While some of these associations have evolved into obligate relationships with permanent symbiont localization within the host, the majority require colonization of every host generation from the environment or via maternal provisions. Across the broad diversity of host species and tissue types that beneficial bacteria can colonize, there are some highly specialized strategies for establishment yet also some common patterns in the molecular basis of colonization. This review focuses on the mechanisms underlying the early stage of beneficial bacterium-invertebrate associations, from initial contact to the establishment of the symbionts in a specific location of the host's body. We first reflect on general selective pressures that can drive the transition from a free-living to a host-associated lifestyle in bacteria. We then cover bacterial molecular factors for colonization in symbioses from both model and nonmodel invertebrate systems where these have been studied, including terrestrial and aquatic host taxa. Finally, we discuss how interactions between multiple colonizing bacteria and priority effects can influence colonization. Taking the bacterial perspective, we emphasize the importance of developing new experimentally tractable systems to derive general insights into the ecological factors and molecular adaptations underlying the origin and establishment of beneficial symbioses in animals.

Construction of high-density transposon mutant library of Staphylococcus aureus using bacteriophage ϕ11

Transposon mutant libraries are an important resource to study bacterial metabolism and pathogenesis. The fitness analysis of mutants in the libraries under various growth conditions provides important clues to study the physiology and biogenesis of structural components of a bacterial cell. A transposon library in conjunction with next-generation sequencing techniques, collectively named transposon sequencing (Tn-seq), enables high-throughput genome profiling and synthetic lethality analysis. Tn-seq has also been used to identify essential genes and to study the mode of action of antibacterials. To construct a high-density transposon mutant library, an efficient delivery system for transposition in a model bacterium is essential. Here, I describe a detailed protocol for generating a high-density phage-based transposon mutant library in a Staphylococcus aureus strain, and this protocol is readily applicable to other S. aureus strains including USA300 and MW2.

Surviving the host: Microbial metabolic genes required for growth of Pseudomonas aeruginosa in physiologically-relevant conditions

Pseudomonas aeruginosa, like other pathogens, adapts to the limiting nutritional environment of the host by altering patterns of gene expression and utilizing alternative pathways required for survival. Understanding the genes essential for survival in the host gives insight into pathways that this organism requires during infection and has the potential to identify better ways to treat infections. Here, we used a saturated transposon insertion mutant pool of P. aeruginosa strain PAO1 and transposon insertion sequencing (Tn-Seq), to identify genes conditionally important for survival under conditions mimicking the environment of a nosocomial infection. Conditions tested included tissue culture medium with and without human serum, a murine abscess model, and a human skin organoid model. Genes known to be upregulated during infections, as well as those involved in nucleotide metabolism, and cobalamin (vitamin B₁₂) biosynthesis, etc., were required for survival in vivo- and in host mimicking conditions, but not in nutrient rich lab medium, Mueller Hinton broth (MHB). Correspondingly, mutants in genes encoding proteins of nucleotide and cobalamin metabolism pathways were shown to have growth defects under physiologically-relevant media conditions, in vivo, and in vivo-like models, and were downregulated in expression under these conditions, when compared to MHB. This study provides evidence for the relevance of studying P. aeruginosa fitness in physiologically-relevant host mimicking conditions and identified metabolic pathways that represent potential novel targets for alternative therapies.

Transposon-Directed Insertion-Site Sequencing Reveals Glycolysis Gene gpmA as Part of the H2O2 Defense Mechanisms in Escherichia coli

Hydrogen peroxide (H₂O₂) is a common effector of defense mechanisms against pathogenic infections. However, bacterial factors involved in H₂O₂ tolerance remain unclear. Here we used transposon-directed insertion-site sequencing (TraDIS), a technique allowing the screening of the whole genome, to identify genes implicated in H₂O₂ tolerance in Escherichia coli. Our TraDIS analysis identified 10 mutants with fitness defect upon H₂O₂ exposure, among which previously H₂O₂-associated genes (oxyR, dps, dksA, rpoS, hfq and polA) and other genes with no known association with H₂O₂ tolerance in E. coli (corA, rbsR, nhaA and gpmA). This is the first description of the impact of gpmA, a gene involved in glycolysis, on the susceptibility of E. coli to H₂O₂. Indeed, confirmatory experiments showed that the deletion of gpmA led to a specific hypersensitivity to H₂O₂ comparable to the deletion of the major H₂O₂ scavenger gene katG. This hypersensitivity was not due to an alteration of catalase function and was independent of the carbon source or the presence of oxygen. Transcription of gpmA was upregulated under H₂O₂ exposure, highlighting its role under oxidative stress. In summary, our TraDIS approach identified gpmA as a member of the oxidative stress defense mechanism in E. coli.

Genome-wide analysis of Brucella melitensis genes required throughout intranasal infection in mice

Brucellae are facultative intracellular Gram-negative coccobacilli that chronically infect various mammals and cause brucellosis. Human brucellosis is among the most common bacterial zoonoses and the vast majority of cases are attributed to B. melitensis. Using transposon sequencing (Tn-seq) analysis, we showed that among 3369 predicted genes of the B. melitensis genome, 861 are required for optimal growth in rich medium and 186 additional genes appeared necessary for survival of B. melitensis in RAW 264.7 macrophages in vitro. As the mucosal immune system represents the first defense against Brucella infection, we investigated the early phase of pulmonary infection in mice. In situ analysis at the single cell level indicates a succession of killing and growth phases, followed by heterogenous proliferation of B. melitensis in alveolar macrophages during the first 48 hours of infection. Tn-seq analysis identified 94 additional genes that are required for survival in the lung at 48 hours post infection. Among them, 42 genes are common to RAW 264.7 macrophages and the lung conditions, including the T4SS and purine synthesis genes. But 52 genes are not identified in RAW 264.7 macrophages, including genes implicated in lipopolysaccharide (LPS) biosynthesis, methionine transport, tryptophan synthesis as well as fatty acid and carbohydrate metabolism. Interestingly, genes implicated in LPS synthesis and β oxidation of fatty acids are no longer required in Interleukin (IL)-17RA^-/- mice and asthmatic mice, respectively. This demonstrates that the immune status determines which genes are required for optimal survival and growth of B. melitensis in vivo.

Brucellosis is one of the most widespread bacterial zoonoses worldwide. Using transposon sequencing (Tn-seq) analysis, we showed that among 3369 predicted genes of the Brucella melitensis genome, 861 are required for optimal growth in rich medium and 186 additional genes appeared necessary for survival of B. melitensis in RAW 264.7 macrophages in vitro. We also investigated the early phase of pulmonary infection in mice and identified 94 additional genes that are required for survival in the lung at 48 hours post infection. Among them, 42 genes are common to RAW 264.7 macrophages and the lung conditions, including the T4SS and purine synthesis genes. But 52 genes are not identified in RAW 264.7 macrophages, including genes implicated in lipopolysaccharide (LPS) biosynthesis, methionine transport, tryptophan synthesis as well as fatty acid and carbohydrate metabolism. Interestingly, genes implicated in LPS synthesis and β oxidation of fatty acids are no longer required in Interleukin (IL)-17RA^-/- mice and asthmatic mice, respectively. Our work demonstrates that both the immune status and the nature of the infected cell type determines which genes are required for optimal survival and growth of B. melitensis in vivo.

Functional Genomics Identified Novel Genes Involved in Growth at Low Temperatures in Listeria monocytogenes

Listeria monocytogenes (Lm) is a foodborne pathogen that can cause severe human illness. Standard control measures for restricting bacterial growth, such as refrigeration, are often inadequate as Lm grows well at low temperatures. To identify genes involved in growth at low temperatures, a powerful functional genomics method Tn-seq was performed in this study. This genome-wide screening comprehensively identified the known and novel genetic determinants involved in low-temperature growth. A novel gene lmo1366, encoding rRNA methyltransferase, was identified to play an essential role in Lm growth at 16°C. In contrast, the inactivation of lmo2301, a gene encoding the terminase of phage A118, significantly enhanced the growth of Lm at 16°C. The deletion of lmo1366 or lmo2301 resulted in cell morphology alterations and impaired the growth rate in milk and other conditions at low temperatures. Transcriptomic analysis revealed that the Δlmo1366 and Δlmo2301 mutants exhibited altered transcriptional patterns compared to the wild-type strain at 16°C with significant differences in genes involved in ribosome structural stability and function, and membrane biogenesis, respectively. This work uncovered novel genetic determinants involved in Lm growth at 16°C, which could lead to a better understanding of how bacteria survive and multiply at low temperatures. Furthermore, these findings could potentially contribute to developing novel antibacterial strategies under low-temperature conditions.

IMPORTANCEListeria monocytogenes is a Gram-positive pathogen that contributes to foodborne outbreaks due to its ability to survive at low temperatures. However, the genetic determinants of Lm involved in growth at low temperatures have not been fully defined. Here, the genetic determinants involved in the low-temperature growth of Lm were comprehensively identified on a genome-wide scale by Tn-seq. The gene lmo1366, encoding rRNA methyltransferase, was identified essential for growth under low-temperature conditions. On the other hand, the gene lmo2301, encoding terminase of phage A118, plays a negative role in bacterial growth at low temperatures. The transcriptomic analysis revealed the potential mechanisms. These findings lead to a better understanding of how bacteria survive and multiply at low temperatures and could provide unique targets for novel antibacterial strategies under low-temperature conditions.

Experimental and Analytical Approaches for Improving the Resolution of Randomly Barcoded Transposon Insertion Sequencing (RB-TnSeq) Studies

Randomly barcoded transposon insertion sequencing (RB-TnSeq) is an efficient, multiplexed method to determine microbial gene function during growth under a selection condition of interest. This technique applies to growth, tolerance, and persistence studies in a variety of hosts, but the wealth of data generated can complicate the identification of the most critical gene targets. Experimental and analytical methods for improving the resolution of RB-TnSeq are proposed, using Pseudomonas putida KT2440 as an example organism. Several key parameters, such as baseline media selection, substantially influence the determination of gene fitness. We also present options to increase statistical confidence in gene fitness, including increasing the number of biological replicates and passaging the baseline culture in parallel with selection conditions. These considerations provide practitioners with several options to identify genes of importance in TnSeq data sets, thereby streamlining metabolic characterization.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

A Comprehensive Overview of the Genes and Functions Required for Lettuce Infection by the Hemibiotrophic Phytopathogen Xanthomonas hortorum pv. vitians

The successful infection of a host plant by a phytopathogenic bacterium depends on a finely tuned molecular cross talk between the two partners. Thanks to transposon insertion sequencing techniques (Tn-seq), whole genomes can now be assessed to determine which genes are important for the fitness of several plant-associated bacteria in planta. Despite its agricultural relevance, the dynamic molecular interaction established between the foliar hemibiotrophic phytopathogen Xanthomonas hortorum pv. vitians and its host, lettuce (Lactuca sativa), remains completely unknown. To decipher the genes and functions mobilized by the pathogen throughout the infection process, we conducted a Tn-seq experiment in lettuce leaves to mimic the selective pressure occurring during natural infection. This genome-wide screening identified 170 genes whose disruption caused serious fitness defects in lettuce. A thorough examination of these genes using comparative genomics and gene set enrichment analyses highlighted that several functions and pathways were highly critical for the pathogen's survival. Numerous genes involved in amino acid, nucleic acid, and exopolysaccharide biosynthesis were critical. The xps type II secretion system operon, a few TonB-dependent transporters involved in carbohydrate or siderophore scavenging, and multiple genes of the carbohydrate catabolism pathways were also critical, emphasizing the importance of nutrition systems in a nutrient-limited environment. Finally, several genes implied in camouflage from the plant immune system and resistance to immunity-induced oxidative stress were strongly involved in host colonization. As a whole, these results highlight some of the central metabolic pathways and cellular functions critical for Xanthomonas host adaptation and pathogenesis. IMPORTANCE Xanthomonas hortorum was recently the subject of renewed interest, as several studies highlighted that its members were responsible for diseases in a wide range of plant species, including crops of agricultural relevance (e.g., tomato and carrot). Among X. hortorum variants, X. hortorum pv. vitians is a reemerging foliar hemibiotrophic phytopathogen responsible for severe outbreaks of bacterial leaf spot of lettuce all around the world. Despite recent findings, sustainable and practical means of disease control remain to be developed. Understanding the host-pathogen interaction from a molecular perspective is crucial to support these efforts. The genes and functions mobilized by X. hortorum pv. vitians during its interaction with lettuce had never been investigated. Our study sheds light on these processes by screening the whole pathogen genome for genes critical for its fitness during the infection process, using transposon insertion sequencing and comparative genomics.

Exploiting transposons in the study of Staphylococcus aureus pathogenesis and virulence

The opportunistic pathogen Staphylococcus aureus has an extremely complex relationship with humans. While the bacteria can exist as a commensal in many, it can cause a wide range of diseases and infections when turned pathogenic. Its presence is a determinant of chronicity and poor prognosis in numerous diseases, and its genomic plasticity causes S. aureus antimicrobial resistance to be one of the most dire contemporary medical problems to solve. Genetic manipulation of S. aureus has led to numerous findings that are vital in the fight against its pathogenesis. The utilisation of transposon mutant libraries for the systematic inspection of the S. aureus genome led to many landmark discoveries pertaining to the bacteria's pathogenicity, antimicrobial resistance acquisition, and virulence regulation. In this review, we describe mutant libraries, and their significant contributions, from various S. aureus strains created with commonly used transposons. The general workflow for the construction of libraries will be presented, along with a discussion of the challenges of undertaking the task of large-scale library construction. As the accessibility of transposon mutant library construction, screening, and analysis increases, this genetic tool could be further exploited in the study of the S. aureus genome.

Using Structure-guided Fragment-Based Drug Discovery to Target Pseudomonas aeruginosa Infections in Cystic Fibrosis

Cystic fibrosis (CF) is progressive genetic disease that predisposes lungs and other organs to multiple long-lasting microbial infections. Pseudomonas aeruginosa is the most prevalent and deadly pathogen among these microbes. Lung function of CF patients worsens following chronic infections with P. aeruginosa and is associated with increased mortality and morbidity. Emergence of multidrug-resistant, extensively drug-resistant and pandrug-resistant strains of P. aeruginosa due to intrinsic and adaptive antibiotic resistance mechanisms has failed the current anti-pseudomonal antibiotics. Hence new antibacterials are urgently needed to treat P. aeruginosa infections. Structure-guided fragment-based drug discovery (FBDD) is a powerful approach in the field of drug development that has succeeded in delivering six FDA approved drugs over the past 20 years targeting a variety of biological molecules. However, FBDD has not been widely used in the development of anti-pseudomonal molecules. In this review, we first give a brief overview of our structure-guided FBDD pipeline and then give a detailed account of FBDD campaigns to combat P. aeruginosa infections by developing small molecules having either bactericidal or anti-virulence properties. We conclude with a brief overview of the FBDD efforts in our lab at the University of Cambridge towards targeting P. aeruginosa infections.

Required Gene Set for Autotrophic Growth of Clostridium autoethanogenum

The majority of the genes present in bacterial genomes remain poorly characterized, with up to one-third of those that are protein encoding having no definitive function. Transposon insertion sequencing represents a high-throughput technique that can help rectify this deficiency. The technology, however, can only be realistically applied to those species in which high rates of DNA transfer can be achieved. Here, we have developed a number of approaches that overcome this barrier in the autotrophic species Clostridium autoethanogenum by using a mariner-based transposon system. The inherent instability of such systems in the Escherichia coli conjugation donor due to transposition events was counteracted through the incorporation of a conditionally lethal codA marker on the plasmid backbone. Relatively low frequencies of transformation of the plasmid into C. autoethanogenum were circumvented through the use of a plasmid that is conditional for replication coupled with the routine implementation of an Illumina library preparation protocol that eliminates plasmid-based reads. A transposon library was then used to determine the essential genes needed for growth using carbon monoxide as the sole carbon and energy source.

IMPORTANCE Although microbial genome sequences are relatively easily determined, assigning gene function remains a bottleneck. Consequently, relatively few genes are well characterized, leaving the function of many as either hypothetical or entirely unknown. High-throughput transposon sequencing can help remedy this deficiency, but is generally only applicable to microbes with efficient DNA transfer procedures. These exclude many microorganisms of importance to humankind either as agents of disease or as industrial process organisms. Here, we developed approaches to facilitate transposon insertion sequencing in the acetogen Clostridium autoethanogenum, a chassis being exploited to convert single-carbon waste gases CO and CO₂ into chemicals and fuels at an industrial scale. This allowed the determination of gene essentiality under heterotrophic and autotrophic growth, providing insights into the utilization of CO as a sole carbon and energy source. The strategies implemented are translatable and will allow others to apply transposon insertion sequencing to other microbes where DNA transfer has until now represented a barrier to progress.

Long-read sequencing for identification of insertion sites in large transposon mutant libraries

Transposon insertion site sequencing (TIS) is a powerful method for associating genotype to phenotype. However, all TIS methods described to date use short nucleotide sequence reads which cannot uniquely determine the locations of transposon insertions within repeating genomic sequences where the repeat units are longer than the sequence read length. To overcome this limitation, we have developed a TIS method using Oxford Nanopore sequencing technology that generates and uses long nucleotide sequence reads; we have called this method LoRTIS (Long-Read Transposon Insertion-site Sequencing). LoRTIS enabled the unique localisation of transposon insertion sites within long repetitive genetic elements of E. coli, such as the transposase genes of insertion sequences and copies of the ~ 5 kb ribosomal RNA operon. We demonstrate that LoRTIS is reproducible, gives comparable results to short-read TIS methods for essential genes, and better resolution around repeat elements. The Oxford Nanopore sequencing device that we used is cost-effective, small and easily portable. Thus, LoRTIS is an efficient means of uniquely identifying transposon insertion sites within long repetitive genetic elements and can be easily transported to, and used in, laboratories that lack access to expensive DNA sequencing facilities.

Identification of putative essential protein domains from high-density transposon insertion sequencing

A first clue to gene function can be obtained by examining whether a gene is required for life in certain standard conditions, that is, whether a gene is essential. In bacteria, essential genes are usually identified by high-density transposon mutagenesis followed by sequencing of insertion sites (Tn-seq). These studies assign the term "essential" to whole genes rather than the protein domain sequences that encode the essential functions. However, genes can code for multiple protein domains that evolve their functions independently. Therefore, when essential genes code for more than one protein domain, only one of them could be essential. In this study, we defined this subset of genes as "essential domain-containing" (EDC) genes. Using a Tn-seq data set built-in Burkholderia cenocepacia K56-2, we developed an in silico pipeline to identify EDC genes and the essential protein domains they encode. We found forty candidate EDC genes and demonstrated growth defect phenotypes using CRISPR interference (CRISPRi). This analysis included two knockdowns of genes encoding the protein domains of unknown function DUF2213 and DUF4148. These putative essential domains are conserved in more than two hundred bacterial species, including human and plant pathogens. Together, our study suggests that essentiality should be assigned to individual protein domains rather than genes, contributing to a first functional characterization of protein domains of unknown function.

High-throughput identification of genes influencing the competitive ability to obtain nutrients and performance of biocontrol in Pseudomonas putida JBC17

Elucidating underlying mechanisms of biocontrol agents (BCAs) could aid in selecting potent BCAs and increasing their biocontrol efficacy. Nutrient competition is an important biocontrol mechanism; however, essential nutrient sources, and contributing genes for nutrient competition still remain to be explored. Pseudomonas putida JBC17 (JBC17WT) suppressed green mold in satsuma mandarins by inhibiting conidial germination of Penicillium digitatum via nutrient competition. To analyze genes essential for biocontrol performance of JBC17WT, we generated a transposon (Tn)-mediated mutant library and selected mutants with the ability to suppress conidial germination. Several mutants in the genes of flagella-formation, including fliR, fliH, and flgG, increased biocontrol performance and enhanced inhibition of conidial germination. They lost swimming motility, exhibited increased growth and rapid carbon and nitrogen utilization than the wild type under nutrient-poor conditions. The nutrient competition assay using polytetrafluoroethylene cylinders revealed that conidial germination was inhibited by nutrient absorption under nutrient-poor conditions. In addition, genes, including amidohydrolase (ytcJ), tonB-dependent receptor (cirA), argininosuccinate synthase (argG), D-3-phosphoglycerate dehydrogenase (serA), and chaperone protein (dnaJ), were involved in the inhibition of conidial germination. The results of this study indicate that rapid and continuous absorption of nutrients by JBC17WT restrict nutrient availability for conidial germination on nutrient-limited fruit surfaces, thereby decreasing the chances of fungal spores infecting fruits. The high-throughput analysis of Tn mutants of this study highlighted the importance of nutrient competition and the genes that influence biocontrol ability, which contributes to the development of biocontrol applications.

Analysis of Gene Essentiality from TnSeq Data Using Transit

TnSeq, or sequencing of transposon insertion libraries, has proven to be a valuable method for probing the functions of genes in a wide range of bacteria. TnSeq has found many applications for studying genes involved in core functions (such as cell division or metabolism), stress response, virulence, etc., as well as to identify potential drug targets. Two of the most commonly used transposons in practice are Himar1, which inserts randomly at TA dinucleotides, and Tn5, which can insert more broadly throughout the genome. These insertions cause putative gene function disruption, and clones with insertions in genes that cannot tolerate disruption (in a given condition) are eliminated from the population. Deep sequencing can be used to efficiently profile the surviving members, with insertions in genes that can be inferred to be non-essential. Data from TnSeq experiments (i.e. transposon insertion counts at specific genomic locations) is inherently noisy, making rigorous statistical analysis (e.g. quantifying significance) challenging. In this chapter, we describe Transit, a Python-based software package for analyzing TnSeq data that combines a variety of data processing tools, quality assessment methods, and analytical algorithms for identifying essential (or conditionally essential) genes.