ARTIST: high-resolution genome-wide assessment of fitness using transposon-insertion sequencing

High-Resolution Whole-Genome Analysis of Sister-Chromatid Contacts

Sister-chromatid cohesion describes the orderly association of newly replicated DNA molecules behind replication forks. It plays an essential role in the maintenance and faithful transmission of genetic information. Cohesion is created by DNA topological links and proteinaceous bridges, whose formation and deposition could be potentially affected by many processes. Current knowledge on cohesion has been mainly gained by fluorescence microscopy observation. However, the resolution limit of microscopy and the restricted number of genomic positions that can be simultaneously visualized considerably hampered progress. Here, we present a high-throughput methodology to monitor sister-chromatid contacts (Hi-SC2). Using the multi-chromosomal Vibrio cholerae bacterium as a model, we show that Hi-SC2 permits to monitor local variations in sister-chromatid cohesion at a high resolution over a whole genome.

A decade of advances in transposon-insertion sequencing

It has been 10 years since the introduction of modern transposon-insertion sequencing (TIS) methods, which combine genome-wide transposon mutagenesis with high-throughput sequencing to estimate the fitness contribution or essentiality of each genetic component in a bacterial genome. Four TIS variations were published in 2009: transposon sequencing (Tn-Seq), transposon-directed insertion site sequencing (TraDIS), insertion sequencing (INSeq) and high-throughput insertion tracking by deep sequencing (HITS). TIS has since become an important tool for molecular microbiologists, being one of the few genome-wide techniques that directly links phenotype to genotype and ultimately can assign gene function. In this Review, we discuss the recent applications of TIS to answer overarching biological questions. We explore emerging and multidisciplinary methods that build on TIS, with an eye towards future applications.

In this Review, several experts discuss progress in the decade since the development of transposon-based approaches for bacterial genetic screens. They describe how advances in both experimental technologies and analytical strategies are resulting in insights into diverse biological processes.

Tn-Seq reveals hidden complexity in the utilization of host-derived glutathione in Francisella tularensis

Host-derived glutathione (GSH) is an essential source of cysteine for the intracellular pathogen Francisella tularensis. In a comprehensive transposon insertion sequencing screen, we identified several F. tularensis genes that play central and previously unappreciated roles in the utilization of GSH during the growth of the bacterium in macrophages. We show that one of these, a gene we named dptA, encodes a proton-dependent oligopeptide transporter that enables growth of the organism on the dipeptide Cys-Gly, a key breakdown product of GSH generated by the enzyme γ-glutamyltranspeptidase (GGT). Although GGT was thought to be the principal enzyme involved in GSH breakdown in F. tularensis, our screen identified a second enzyme, referred to as ChaC, that is also involved in the utilization of exogenous GSH. However, unlike GGT and DptA, we show that the importance of ChaC in supporting intramacrophage growth extends beyond cysteine acquisition. Taken together, our findings provide a compendium of F. tularensis genes required for intracellular growth and identify new players in the metabolism of GSH that could be attractive targets for therapeutic intervention.

Selection or drift: The population biology underlying transposon insertion sequencing experiments

Transposon insertion sequencing methods such as Tn-seq revolutionized microbiology by allowing the identification of genomic loci that are critical for viability in a specific environment on a genome-wide scale. While powerful, transposon insertion sequencing suffers from limited reproducibility when different analysis methods are compared. From the perspective of population biology, this may be explained by changes in mutant frequency due to chance (drift) rather than differential fitness (selection). Here, we develop a mathematical model of the population biology of transposon insertion sequencing experiments, i.e. the changes in size and composition of the transposon-mutagenized population during the experiment. We use this model to investigate mutagenesis, the growth of the mutant library, and its passage through bottlenecks. Specifically, we study how these processes can lead to extinction of individual mutants depending on their fitness and the distribution of fitness effects (DFE) of the entire mutant population. We find that in typical in vitro experiments few mutants with high fitness go extinct. However, bottlenecks of a size that is common in animal infection models lead to so much random extinction that a large number of viable mutants would be misclassified. While mutants with low fitness are more likely to be lost during the experiment, mutants with intermediate fitness are expected to be much more abundant and can constitute a large proportion of detected hits, i.e. false positives. Thus, incorporating the DFEs of randomly generated mutations in the analysis may improve the reproducibility of transposon insertion experiments, especially when strong bottlenecks are encountered.

Heat-shock proteases promote survival of Pseudomonas aeruginosa during growth arrest

When nutrients in their environment are exhausted, bacterial cells become arrested for growth. During these periods, a primary challenge is maintaining cellular integrity with a reduced capacity for renewal or repair. Here, we show that the heat-shock protease FtsH is generally required for growth arrest survival of Pseudomonas aeruginosa, and that this requirement is independent of a role in regulating lipopolysaccharide synthesis, as has been suggested for Escherichia coli We find that ftsH interacts with diverse genes during growth and overlaps functionally with the other heat-shock protease-encoding genes hslVU, lon, and clpXP to promote survival during growth arrest. Systematic deletion of the heat-shock protease-encoding genes reveals that the proteases function hierarchically during growth arrest, with FtsH and ClpXP having primary, nonredundant roles, and HslVU and Lon deploying a secondary response to aging stress. This hierarchy is partially conserved during growth at high temperature and alkaline pH, suggesting that heat, pH, and growth arrest effectively impose a similar type of proteostatic stress at the cellular level. In support of this inference, heat and growth arrest act synergistically to kill cells, and protein aggregation appears to occur more rapidly in protease mutants during growth arrest and correlates with the onset of cell death. Our findings suggest that protein aggregation is a major driver of aging and cell death during growth arrest, and that coordinated activity of the heat-shock response is required to ensure ongoing protein quality control in the absence of growth.

An Oral Inoculation Infant Rabbit Model for Shigella Infection

Shigella species cause diarrheal disease globally. Shigellosis is typically characterized by bloody stools and colitis with mucosal damage and is the leading bacterial cause of diarrheal death worldwide. After the pathogen is orally ingested, it invades and replicates within the colonic epithelium through mechanisms that rely on its type III secretion system (T3SS). Currently, oral infection-based small animal models to study the pathogenesis of shigellosis are lacking. Here, we found that orogastric inoculation of infant rabbits with Shigella flexneri resulted in diarrhea and colonic pathology resembling that found in human shigellosis. Fasting animals prior to S. flexneri inoculation increased the frequency of disease. The pathogen colonized the colon, where both luminal and intraepithelial foci were observed. The intraepithelial foci likely arise through S. flexneri spreading from cell to cell. Robust S. flexneri intestinal colonization, invasion of the colonic epithelium, and epithelial sloughing all required the T3SS as well as IcsA, a factor required for bacterial spreading and adhesion in vitro Expression of the proinflammatory chemokine interleukin 8 (IL-8), detected with in situ mRNA labeling, was higher in animals infected with wild-type S. flexneri versus mutant strains deficient in icsA or T3SS, suggesting that epithelial invasion promotes expression of this chemokine. Collectively, our findings suggest that oral infection of infant rabbits offers a useful experimental model for studies of the pathogenesis of shigellosis and for testing of new therapeutics.IMPORTANCEShigella species are the leading bacterial cause of diarrheal death globally. The pathogen causes bacillary dysentery, a bloody diarrheal disease characterized by damage to the colonic mucosa and is usually spread through the fecal-oral route. Small animal models of shigellosis that rely on the oral route of infection are lacking. Here, we found that orogastric inoculation of infant rabbits with S. flexneri led to a diarrheal disease and colonic pathology reminiscent of human shigellosis. Diarrhea, intestinal colonization, and pathology in this model were dependent on the S. flexneri type III secretion system and IcsA, canonical Shigella virulence factors. Thus, oral infection of infant rabbits offers a feasible model to study the pathogenesis of shigellosis and to develop and test new therapeutics.

Statistical analysis of variability in TnSeq data across conditions using zero-inflated negative binomial regression

Background

Deep sequencing of transposon mutant libraries (or TnSeq) is a powerful method for probing essentiality of genomic loci under different environmental conditions. Various analytical methods have been described for identifying conditionally essential genes whose tolerance for insertions varies between two conditions. However, for large-scale experiments involving many conditions, a method is needed for identifying genes that exhibit significant variability in insertions across multiple conditions.

Results

In this paper, we introduce a novel statistical method for identifying genes with significant variability of insertion counts across multiple conditions based on Zero-Inflated Negative Binomial (ZINB) regression. Using likelihood ratio tests, we show that the ZINB distribution fits TnSeq data better than either ANOVA or a Negative Binomial (in a generalized linear model). We use ZINB regression to identify genes required for infection of M. tuberculosis H37Rv in C57BL/6 mice. We also use ZINB to perform a analysis of genes conditionally essential in H37Rv cultures exposed to multiple antibiotics.

Conclusions

Our results show that, not only does ZINB generally identify most of the genes found by pairwise resampling (and vastly out-performs ANOVA), but it also identifies additional genes where variability is detectable only when the magnitudes of insertion counts are treated separately from local differences in saturation, as in the ZINB model.

A live vaccine rapidly protects against cholera in an infant rabbit model

Outbreaks of cholera, a rapidly fatal diarrheal disease, often spread explosively. The efficacy of reactive vaccination campaigns-deploying Vibrio cholerae vaccines during epidemics-is partially limited by the time required for vaccine recipients to develop adaptive immunity. We created HaitiV, a live attenuated cholera vaccine candidate, by deleting diarrheagenic factors from a recent clinical isolate of V. cholerae and incorporating safeguards against vaccine reversion. We demonstrate that administration of HaitiV 24 hours before lethal challenge with wild-type V. cholerae reduced intestinal colonization by the wild-type strain, slowed disease progression, and reduced mortality in an infant rabbit model of cholera. HaitiV-mediated protection required viable vaccine, and rapid protection kinetics are not consistent with development of adaptive immunity. These features suggest that HaitiV mediates probiotic-like protection from cholera, a mechanism that is not known to be elicited by traditional vaccines. Mathematical modeling indicates that an intervention that works at the speed of HaitiV-mediated protection could improve the public health impact of reactive vaccination.

MarTrack: A versatile toolbox of mariner transposon derivatives used for functional genetic analysis of bacterial genomes

The mariner transposon family of Himar1 has been widely used for the random mutagenesis of bacteria to generate single insertions into the chromosome. Here, a versatile toolbox of mariner transposon derivatives was generated and applied to the functional genomics investigation of fish pathogen Edwardsiella piscicida. In this study, we combined the merits of the random mutagenesis of mariner transposon and common efficient reporter marker genes or regulatory elements, mCherry, gfp, luxAB, lacZ, sacBR, and P_BAD and antibiotic resistance cassettes to construct a series of derivative transposon vectors, pMmch, pMKGR, pMCGR, pMXKGR, pMLKGR, pMSGR, and pMPR, based on the initial transposon pMar2xT7. The function and effectiveness of the modified transposons were verified by introducing them into E. piscicida EIB202. Based on the toolbox, a transposon insertion mutant library containing approximately 3.0 × 10⁵ distinct mutants was constructed to explore the upstream regulators of esrB, the master regulator of the type III and type VI secretion systems (T3/T6SS) in E. piscicida. Following analysis by Con-ARTIST, ETAE_3474, annotated as fabR and involved in fatty acid metabolism, was screened out and identified as a novel regulator mediating T3SS and T6SS expression. In addition, the fabR mutants displayed critical virulence attenuation in turbot. Due to the broad-range host compatibility of mariner transposons, the newly built transposon toolbox can be applied for functional genomics studies in various bacteria.

Transposon-insertion sequencing screens unveil requirements for EHEC growth and intestinal colonization

Enterohemorrhagic Escherichia coli O157:H7 (EHEC) is an important food-borne pathogen that colonizes the colon. Transposon-insertion sequencing (TIS) was used to identify genes required for EHEC and E. coli K-12 growth in vitro and for EHEC growth in vivo in the infant rabbit colon. Surprisingly, many conserved loci contribute to EHEC's but not to K-12's growth in vitro. There was a restrictive bottleneck for EHEC colonization of the rabbit colon, which complicated identification of EHEC genes facilitating growth in vivo. Both a refined version of an existing analytic framework as well as PCA-based analysis were used to compensate for the effects of the infection bottleneck. These analyses confirmed that the EHEC LEE-encoded type III secretion apparatus is required for growth in vivo and revealed that only a few effectors are critical for in vivo fitness. Over 200 mutants not previously associated with EHEC survival/growth in vivo also appeared attenuated in vivo, and a subset of these putative in vivo fitness factors were validated. Some were found to contribute to efficient type-three secretion while others, including tatABC, oxyR, envC, acrAB, and cvpA, promote EHEC resistance to host-derived stresses. cvpA is also required for intestinal growth of several other enteric pathogens, and proved to be required for EHEC, Vibrio cholerae and Vibrio parahaemolyticus resistance to the bile salt deoxycholate, highlighting the important role of this previously uncharacterized protein in pathogen survival. Collectively, our findings provide a comprehensive framework for understanding EHEC growth in the intestine.

Methods for Tn-Seq Analysis in Acinetobacter baumannii

Transposon insertion sequencing (Tn-seq) is a powerful method for identifying genes required for virtually any growth or survival trait in bacteria. The technology employs next-generation DNA sequencing to identify and quantify the relative abundances of individual transposon mutants within complex pools of such mutants. When applied to pools of thousands to millions of random transposon mutants grown under selective pressure, the technique can rapidly identify, at genome scale, the mutants and corresponding genes negatively or positively selected. This chapter presents core protocols for Tn-seq analysis of Acinetobacter baumannii: generation of a high-saturation random transposon mutant pool, and isolation and sequencing of transposon-genome junctions from such a pool for identifying and quantifying the individual mutants. With these tools, the researcher can address diverse biological questions by carrying out selective growth of a mutant pool followed by Tn-seq analysis to identify genotype-phenotype associations.

A Vibrio cholerae BolA-Like Protein Is Required for Proper Cell Shape and Cell Envelope Integrity

BolA family proteins are conserved in Gram-negative bacteria and many eukaryotes. While diverse cellular phenotypes have been linked to this protein family, the molecular pathways through which these proteins mediate their effects are not well described. Here, we investigated the roles of BolA family proteins in Vibrio cholerae, the cholera pathogen. Like Escherichia coli, V. cholerae encodes two BolA proteins, BolA and IbaG. However, in marked contrast to E. coli, where bolA is linked to cell shape and ibaG is not, in V. cholerae, bolA mutants lack morphological defects, whereas ibaG proved critical for the generation and/or maintenance of the pathogen's morphology. Notably, the bizarre-shaped, multipolar, elongated, and wide cells that predominated in exponential-phase ΔibaG V. cholerae cultures were not observed in stationary-phase cultures. The V. cholerae ΔibaG mutant exhibited increased sensitivity to cell envelope stressors, including cell wall-acting antibiotics and bile, and was defective in intestinal colonization. ΔibaG V. cholerae had reduced peptidoglycan and lipid II and altered outer membrane lipids, likely contributing to the mutant's morphological defects and sensitivity to envelope stressors. Transposon insertion sequencing analysis of ibaG's genetic interactions suggested that ibaG is involved in several processes involved in the generation and homeostasis of the cell envelope. Furthermore, copurification studies revealed that IbaG interacts with proteins containing iron-sulfur clusters or involved in their assembly. Collectively, our findings suggest that V. cholerae IbaG controls cell morphology and cell envelope integrity through its role in biogenesis or trafficking of iron-sulfur cluster proteins.IMPORTANCE BolA-like proteins are conserved across prokaryotes and eukaryotes. These proteins have been linked to a variety of phenotypes, but the pathways and mechanisms through which they act have not been extensively characterized. Here, we unraveled the role of the BolA-like protein IbaG in the cholera pathogen Vibrio cholerae The absence of IbaG was associated with dramatic changes in cell morphology, sensitivity to envelope stressors, and intestinal colonization defects. IbaG was found to be required for biogenesis of several components of the V. cholerae cell envelope and to interact with numerous iron-sulfur cluster-containing proteins and factors involved in their assembly. Thus, our findings suggest that IbaG governs V. cholerae cell shape and cell envelope homeostasis through its effects on iron-sulfur proteins and associated pathways. The diversity of processes involving iron-sulfur-containing proteins is likely a factor underlying the range of phenotypes associated with BolA family proteins.

In Vivo Gene Essentiality and Metabolism in Bordetella pertussis

Bordetella pertussis is the causative agent of whooping cough, a serious respiratory illness affecting children and adults, associated with prolonged cough and potential mortality. Whooping cough has reemerged in recent years, emphasizing a need for increased knowledge of basic mechanisms of B. pertussis growth and pathogenicity. While previous studies have provided insight into in vitro gene essentiality of this organism, very little is known about in vivo gene essentiality, a critical gap in knowledge, since B. pertussis has no previously identified environmental reservoir and is isolated from human respiratory tract samples. We hypothesize that the metabolic capabilities of B. pertussis are especially tailored to the respiratory tract and that many of the genes involved in B. pertussis metabolism would be required to establish infection in vivo In this study, we generated a diverse library of transposon mutants and then used it to probe gene essentiality in vivo in a murine model of infection. Using the CON-ARTIST pipeline, 117 genes were identified as conditionally essential at 1 day postinfection, and 169 genes were identified as conditionally essential at 3 days postinfection. Most of the identified genes were associated with metabolism, and we utilized two existing genome-scale metabolic network reconstructions to probe the effects of individual essential genes on biomass synthesis. This analysis suggested a critical role for glucose metabolism and lipooligosaccharide biosynthesis in vivo This is the first genome-wide evaluation of in vivo gene essentiality in B. pertussis and provides tools for future exploration.IMPORTANCE Our study describes the first in vivo transposon sequencing (Tn-seq) analysis of B. pertussis and identifies genes predicted to be essential for in vivo growth in a murine model of intranasal infection, generating key resources for future investigations into B. pertussis pathogenesis and vaccine design.

Suppression of Alternative Lipooligosaccharide Glycosyltransferase Activity by UDP-Galactose Epimerase Enhances Murine Lung Infection and Evasion of Serum IgM

In pathogens that produce lipooligosaccharide (LOS), sugar residues within the surface-exposed LOS outer core mediate interactions with components of the host immune system, promoting bacterial infection. Many LOS structures are controlled by phase variation mediated by random slipped-strand base mispairing, which can reversibly switch gene expression on or off. Phase variation diversifies the LOS, however its adaptive role is not well-understood. Nontypeable Haemophilus influenzae (NTHi) is an important pathogen that causes a range of illnesses in the upper and lower respiratory tract. In NTHi a phase variable galactosyltransferase encoded by lic2A initiates galactose chain extension of the LOS outer core. The donor substrate for Lic2A, UDP-galactose, is generated from UDP-glucose by UDP-galactose epimerase encoded by galE. Our previous fitness profiling of H. influenzae mutants in a murine lung model showed that the galE mutant had a severe survival defect, while the lic2A mutant's defect was modest, leading us to postulate that unidentified factors act as suppressors of potential defects in a lic2A mutant. Herein we conducted a genome-wide genetic interaction screen to identify genes epistatic on lic2A for survival in the murine lung. An unexpected finding was that galE mutants exhibited restored virulence properties in a lic2A mutant background. We identified an alternative antibody epitope generated by Lic2A in the galE mutant that increased sensitivity to classical complement mediated killing in human serum. Deletion of lic2A or restoration of UDP-galactose synthesis alleviated the galE mutant's virulence defects. These studies indicate that when deprived of its galactosyl substrate, Lic2A acquires an alternative activity leading to increased recognition of NTHi by IgM and decreased survival in the lung model. Biofilm formation was increased by deletion of galE and by increased availability of UDP-GlcNAc precursors that can compete with UDP-galactose production. NTHi's ability to reversibly inactivate lic2A by phase-variation may influence survival in niches of infection in which UDP-Galactose levels are limiting.

Defining the core essential genome of Pseudomonas aeruginosa

Genomics offered the promise of transforming antibiotic discovery by revealing many new essential genes as good targets, but the results fell short of the promise. While numerous factors contributed to the disappointing yield, one factor was that essential genes for a bacterial species were often defined based on a single or limited number of strains grown under a single or limited number of in vitro laboratory conditions. In fact, the essentiality of a gene can depend on both the genetic background and growth condition. We thus developed a strategy for more rigorously defining the core essential genome of a bacterial species by studying many pathogen strains and growth conditions. We assessed how many strains must be examined to converge on a set of core essential genes for a species. We used transposon insertion sequencing (Tn-Seq) to define essential genes in nine strains of Pseudomonas aeruginosa on five different media and developed a statistical model, FiTnEss, to classify genes as essential versus nonessential across all strain-medium combinations. We defined a set of 321 core essential genes, representing 6.6% of the genome. We determined that analysis of four strains was typically sufficient in P. aeruginosa to converge on a set of core essential genes likely to be essential across the species across a wide range of conditions relevant to in vivo infection, and thus to represent attractive targets for novel drug discovery.

The biotroph Agrobacterium tumefaciens thrives in tumors by exploiting a wide spectrum of plant host metabolites

Agrobacterium tumefaciens is a niche-constructing biotroph that exploits host plant metabolites. We combined metabolomics, transposon-sequencing (Tn-seq), transcriptomics, and reverse genetics to characterize A. tumefaciens pathways involved in the exploitation of resources from the Solanum lycopersicum host plant. Metabolomics of healthy stems and plant tumors revealed the common (e.g. sucrose, glutamate) and enriched (e.g. opines, γ-aminobutyric acid (GABA), γ-hydroxybutyric acid (GHB), pyruvate) metabolites that A. tumefaciens could use as nutrients. Tn-seq and transcriptomics pinpointed the genes that are crucial and/or upregulated when the pathogen grew on either sucrose (pgi, kdgA, pycA, cisY) or GHB (blcAB, pckA, eno, gpsA) as a carbon source. While sucrose assimilation involved the Entner-Doudoroff and tricarboxylic acid (TCA) pathways, GHB degradation required the blc genes, TCA cycle, and gluconeogenesis. The tumor-enriched metabolite pyruvate is at the node connecting these pathways. Using reverse genetics, we showed that the blc, pckA, and pycA loci were important for aggressiveness (tumor weight), proliferation (bacterial charge), and/or fitness (competition between the constructed mutants and wild-type) of A. tumefaciens in plant tumors. This work highlighted how a biotroph mobilizes its central metabolism for exploiting a wide diversity of resources in a plant host. It further shows the complementarity of functional genome-wide scans by transcriptomics and Tn-seq to decipher the lifestyle of a plant pathogen.

Genome-Wide Screening for Enteric Colonization Factors in Carbapenem-Resistant ST258 Klebsiella pneumoniae

A diverse, antibiotic-naive microbiota prevents highly antibiotic-resistant microbes, including carbapenem-resistant Klebsiella pneumoniae (CR-Kp), from achieving dense colonization of the intestinal lumen. Antibiotic-mediated destruction of the microbiota leads to expansion of CR-Kp in the gut, markedly increasing the risk of bacteremia in vulnerable patients. While preventing dense colonization represents a rational approach to reduce intra- and interpatient dissemination of CR-Kp, little is known about pathogen-associated factors that enable dense growth and persistence in the intestinal lumen. To identify genetic factors essential for dense colonization of the gut by CR-Kp, we constructed a highly saturated transposon mutant library with >150,000 unique mutations in an ST258 strain of CR-Kp and screened for in vitro growth and in vivo intestinal colonization in antibiotic-treated mice. Stochastic and partially reversible fluctuations in the representation of different mutations during dense colonization revealed the dynamic nature of intestinal microbial populations. We identified genes that are crucial for early and late stages of dense gut colonization and confirmed their role by testing isogenic mutants in in vivo competition assays with wild-type CR-Kp Screening of the transposon library also identified mutations that enhanced in vivo CR-Kp growth. These newly identified colonization factors may provide novel therapeutic opportunities to reduce intestinal colonization by CR-KpIMPORTANCEKlebsiella pneumoniae is a common cause of bloodstream infections in immunocompromised and hospitalized patients, and over the last 2 decades, some strains have acquired resistance to nearly all available antibiotics, including broad-spectrum carbapenems. The U.S. Centers for Disease Control and Prevention has listed carbapenem-resistant K. pneumoniae (CR-Kp) as an urgent public health threat. Dense colonization of the intestine by CR-Kp and other antibiotic-resistant bacteria is associated with an increased risk of bacteremia. Reducing the density of gut colonization by CR-Kp is likely to reduce their transmission from patient to patient in health care facilities as well as systemic infections. How CR-Kp expands and persists in the gut lumen, however, is poorly understood. Herein, we generated a highly saturated mutant library in a multidrug-resistant K. pneumoniae strain and identified genetic factors that are associated with dense gut colonization by K. pneumoniae This study sheds light on host colonization by K. pneumoniae and identifies potential colonization factors that contribute to high-density persistence of K. pneumoniae in the intestine.

The bare necessities: Uncovering essential and condition-critical genes with transposon sequencing

Querying gene function in bacteria has been greatly accelerated by the advent of transposon sequencing (Tn-seq) technologies (related Tn-seq strategies are known as TraDIS, INSeq, RB-TnSeq, and HITS). Pooled populations of transposon mutants are cultured in an environment and next-generation sequencing tools are used to determine areas of the genome that are important for bacterial fitness. In this review we provide an overview of Tn-seq methodologies and discuss how Tn-seq has been applied, or could be applied, to the study of oral microbiology. These applications include studying the essential genome as a means to rationally design therapeutic agents. Tn-seq has also contributed to our understanding of well-studied biological processes in oral bacteria. Other important applications include in vivo pathogenesis studies and use of Tn-seq to probe the molecular basis of microbial interactions. We also highlight recent advancements in techniques that act in synergy with Tn-seq such as clustered regularly interspaced short palindromic repeats (CRISPR) interference and microfluidic chip platforms.

Unsupervised Learning Approach for Comparing Multiple Transposon Insertion Sequencing Studies

Transposon insertion sequencing (TIS) is a widely used technique for conducting genome-scale forward genetic screens in bacteria. However, few methods enable comparison of TIS data across multiple replicates of a screen or across independent screens, including screens performed in different organisms. Here, we introduce a post hoc analytic framework, comparative TIS (CompTIS), which utilizes unsupervised learning to enable meta-analysis of multiple TIS data sets. CompTIS first implements screen-level principal-component analysis (PCA) and clustering to identify variation between the TIS screens. This initial screen-level analysis facilitates the selection of related screens for additional analyses, reveals the relatedness of complex environments based on growth phenotypes measured by TIS, and provides a useful quality control step. Subsequently, PCA is performed on genes to identify loci whose corresponding mutants lead to concordant/discordant phenotypes across all or in a subset of screens. We used CompTIS to analyze published intestinal colonization TIS data sets from two vibrio species. Gene-level analyses identified both pan-vibrio genes required for intestinal colonization and conserved genes that displayed species-specific requirements. CompTIS is applicable to virtually any combination of TIS screens and can be implemented without regard to either the number of screens or the methods used for upstream data analysis.IMPORTANCE Forward genetic screens are powerful tools for functional genomics. The comparison of similar forward genetic screens performed in different organisms enables the identification of genes with similar or different phenotypes across organisms. Transposon insertion sequencing is a widely used method for conducting genome-scale forward genetic screens in bacteria, yet few bioinformatic approaches have been developed to compare the results of screen replicates and different screens conducted across species or strains. Here, we used principal-component analysis (PCA) and hierarchical clustering, two unsupervised learning approaches, to analyze the relatedness of multiple in vivo screens of pathogenic vibrios. This analytic framework reveals both shared pan-vibrio requirements for intestinal colonization and strain-specific dependencies. Our findings suggest that PCA-based analytics will be a straightforward widely applicable approach for comparing diverse transposon insertion sequencing screens.

MarTrack: A versatile toolbox of mariner transposon derivatives used for functional genetic analysis of bacterial genomes

The mariner transposon family of Himar1 has been widely used for the random mutagenesis of bacteria to generate single insertions into the chromosome. Here, a versatile toolbox of mariner transposon derivatives was generated and applied to the functional genomics investigation of fish pathogen Edwardsiella piscicida. In this study, we combined the merits of the random mutagenesis of mariner transposon and common efficient reporter marker genes or regulatory elements, mcherry, gfp, luxAB, lacZ, sacBR, and P_BAD and antibiotic resistance cassettes to construct a series of derivative transposon vectors, pMmch, pMKGR, pMCGR, pMXKGR, pMLKGR, pMSGR, and pMPR, based on the initial transposon pMar2xT7. The function and effectiveness of the modified transposons were verified by introducing them into E. piscicida EIB202. Based on the toolbox, a transposon insertion mutant library containing approximately 3.0 × 10⁵ separated mutants was constructed to explore the upstream regulators of esrB, the master regulator of the type III and type VI secretion systems (T3/T6SS) in E. piscicida. Following analysis by Con-ARTIST, ETAE_2184 (renamed as EsrR) was screened out and identified as a novel regulator mediating T3SS expression. In addition, the esrR mutants displayed critical virulence attenuation. Due to the broad-range host compatibility of mariner transposons, the newly built transposon toolbox can be broadly applied for functional genomics studies in various bacteria.

The RNA degradosome promotes tRNA quality control through clearance of hypomodified tRNA

The factors and mechanisms that govern tRNA stability in bacteria are not well understood. Here, we investigated the influence of posttranscriptional modification of bacterial tRNAs (tRNA modification) on tRNA stability. We focused on ThiI-generated 4-thiouridine (s⁴U), a modification found in bacterial and archaeal tRNAs. Comprehensive quantification of Vibrio cholerae tRNAs revealed that the abundance of some tRNAs is decreased in a ΔthiI strain in a stationary phase-specific manner. Multiple mechanisms, including rapid degradation of a subset of hypomodified tRNAs, account for the reduced abundance of tRNAs in the absence of thiI Through transposon insertion sequencing, we identified additional tRNA modifications that promote tRNA stability and bacterial viability. Genetic analysis of suppressor mutants as well as biochemical analyses revealed that rapid degradation of hypomodified tRNA is mediated by the RNA degradosome. Elongation factor Tu seems to compete with the RNA degradosome, protecting aminoacyl tRNAs from decay. Together, our observations describe a previously unrecognized bacterial tRNA quality control system in which hypomodification sensitizes tRNAs to decay mediated by the RNA degradosome.

Comprehensive identification of Vibrio vulnificus genes required for growth in human serum

Vibrio vulnificus can be a highly invasive pathogen capable of spreading from an infection site to the bloodstream, causing sepsis and death. To survive and proliferate in blood, the pathogen requires mechanisms to overcome the innate immune defenses and metabolic limitations of this host niche. We created a high-density transposon mutant library in YJ016, a strain representative of the most virulent V. vulnificus lineage (or phylogroup) and used transposon insertion sequencing (TIS) screens to identify loci that enable the pathogen to survive and proliferate in human serum. Initially, genes underrepresented for insertions were used to estimate the V. vulnificus essential gene set; comparisons of these genes with similar TIS-based classification of underrepresented genes in other vibrios enabled the compilation of a common Vibrio essential gene set. Analysis of the relative abundance of insertion mutants in the library after exposure to serum suggested that genes involved in capsule biogenesis are critical for YJ016 complement resistance. Notably, homologues of two genes required for YJ016 serum-resistance and capsule biogenesis were not previously linked to capsule biogenesis and are largely absent from other V. vulnificus strains. The relative abundance of mutants after exposure to heat inactivated serum was compared with the findings from the serum screen. These comparisons suggest that in both conditions the pathogen relies on its Na⁺ transporting NADH-ubiquinone reductase (NQR) complex and type II secretion system to survive/proliferate within the metabolic constraints of serum. Collectively, our findings reveal the potency of comparative TIS screens to provide knowledge of how a pathogen overcomes the diverse limitations to growth imposed by serum.

Genetic Determinants of Penicillin Tolerance in Vibrio cholerae

Many bacteria are resistant to killing (tolerant) by typically bactericidal antibiotics due to their ability to counteract drug-induced cell damage. Vibrio cholerae, the cholera agent, displays an unusually high tolerance to diverse inhibitors of cell wall synthesis. Exposure to these agents, which in other bacteria leads to lysis and death, results in a breakdown of the cell wall and subsequent sphere formation in V. cholerae Spheres readily recover to rod-shaped cells upon antibiotic removal, but the mechanisms mediating the recovery process are not well characterized. Here, we found that the mechanisms of recovery are dependent on environmental conditions. Interestingly, on agarose pads, spheres undergo characteristic stages during the restoration of rod shape. Drug inhibition and microscopy experiments suggest that class A penicillin binding proteins (aPBPs) play a more active role than the Rod system, especially early in sphere recovery. Transposon insertion sequencing (TnSeq) analyses revealed that lipopolysaccharide (LPS) and cell wall biogenesis genes, as well as the sigma E cell envelope stress response, were particularly critical for recovery. LPS core and O-antigen appear to be more critical for sphere formation/integrity and viability than lipid A modifications. Overall, our findings demonstrate that the outer membrane is a key contributor to beta lactam tolerance and suggest a role for aPBPs in cell wall biogenesis in the absence of rod-shape cues. Factors required for postantibiotic recovery could serve as targets for antibiotic adjuvants that enhance the efficacy of antibiotics that inhibit cell wall biogenesis.

Critical role for a promoter discriminator in RpoS control of virulence in Edwardsiella piscicida

Edwardsiella piscicida is a leading fish pathogen that causes significant economic loses in the aquaculture industry. The pathogen depends on type III and type VI secretion systems (T3/T6SS) for growth and virulence in fish and the expression of both systems is controlled by the EsrB transcription activator. Here, we performed a Tn-seq-based screen to uncover factors that govern esrB expression. Unexpectedly, we discovered that RpoS antagonizes esrB expression and thereby inhibits production of E. piscicida’s T3/T6SS. Using in vitro transcription assays, we showed that RpoS can block RpoD-mediated transcription of esrB. ChIP-seq- and RNA-seq-based profiling, as well as mutational and biochemical analyses revealed that RpoS-repressed promoters contain a -6G in their respective discriminator sequences; moreover, this -6G proved critical for RpoS to inhibit esrB expression. Mutation of the RpoS R99 residue, an amino acid that molecular modeling predicts interacts with -6G in the esrB discriminator, abolished RpoS’ capacity for repression. In a turbot model, an rpoS deletion mutant was attenuated early but not late in infection, whereas a mutant expressing RpoS^R99A exhibited elevated fitness throughout the infection period. Collectively, these findings deepen our understanding of how RpoS can inhibit gene expression and demonstrate the temporal variation in the requirement for this sigma factor during infection.

Edwardsiella piscicida, a major fish pathogen, relies on T3/T6SSs for virulence and the EsrB transcription activator promotes the expression of these secretion systems and many other genes that enable growth in fish. Here, we found that the alternative sigma factor RpoS inhibits expression of esrB thereby diminishing expression of virulence-associated genes. Transcriptome profiling revealed that, as in many other organisms, RpoS enables expression of hundreds of genes, many of which are linked to stress responses, suggesting that RpoS may mediate a trade-off between stress adaptation and virulence. Consistent with this idea, we found that an rpoS mutant was attenuated early, but not late in infection of turbot, whereas an esrB mutant was attenuated late and not early in infection. Molecular analyses demonstrated that RpoS inhibition of esrB expression involves a direct interaction between RpoS and the esrB promoter; in particular, interactions between RpoS residue R99 and the -6G nucleotide in the esrB promoter discriminator appear to be critical for repression of esrB expression. These findings provide new insight into how a sigma factor can impede transcription and demonstrate the temporal dynamics of the requirement for a sigma factor during the course of infection.

Genome-Wide Determination of Gene Essentiality by Transposon Insertion Sequencing in Yeast Pichia pastoris

In many prokaryotes but limited eukaryotic species, the combination of transposon mutagenesis and high-throughput sequencing has greatly accelerated the identification of essential genes. Here we successfully applied this technique to the methylotrophic yeast Pichia pastoris and classified its conditionally essential/non-essential gene sets. Firstly, we showed that two DNA transposons, TcBuster and Sleeping beauty, had high transposition activities in P. pastoris. By merging their insertion libraries and performing Tn-seq, we identified a total of 202,858 unique insertions under glucose supported growth condition. We then developed a machine learning method to classify the 5,040 annotated genes into putatively essential, putatively non-essential, ambig1 and ambig2 groups, and validated the accuracy of this classification model. Besides, Tn-seq was also performed under methanol supported growth condition and methanol specific essential genes were identified. The comparison of conditionally essential genes between glucose and methanol supported growth conditions helped to reveal potential novel targets involved in methanol metabolism and signaling. Our findings suggest that transposon mutagenesis and Tn-seq could be applied in the methylotrophic yeast Pichia pastoris to classify conditionally essential/non-essential gene sets. Our work also shows that determining gene essentiality under different culture conditions could help to screen for novel functional components specifically involved in methanol metabolism.

Genome-wide mutant profiling predicts the mechanism of a Lipid II binding antibiotic

Identifying targets of antibacterial compounds remains a challenging step in antibiotic development. We have developed a two-pronged functional genomics approach to predict mechanism of action that uses mutant fitness data from antibiotic-treated transposon libraries containing both upregulation and inactivation mutants. We treated a Staphylococcus aureus transposon library containing 690,000 unique insertions with 32 antibiotics. Upregulation signatures, identified from directional biases in insertions, revealed known molecular targets and resistance mechanisms for the majority of these. Because single gene upregulation does not always confer resistance, we used a complementary machine learning approach to predict mechanism from inactivation mutant fitness profiles. This approach suggested the cell wall precursor Lipid II as the molecular target of the lysocins, a mechanism we have confirmed. We conclude that docking to membrane-anchored Lipid II precedes the selective bacteriolysis that distinguishes these lytic natural products, showing the utility of our approach for nominating antibiotic mechanism of action.

The Streptococcus agalactiae Stringent Response Enhances Virulence and Persistence in Human Blood

Streptococcus agalactiae (group B Streptococcus [GBS]) causes serious infections in neonates. We previously reported a transposon sequencing (Tn-seq) system for performing genomewide assessment of gene fitness in GBS. In order to identify molecular mechanisms required for GBS to transition from a mucosal commensal lifestyle to bloodstream invasion, we performed Tn-seq on GBS strain A909 with human whole blood. Our analysis identified 16 genes conditionally essential for GBS survival in blood, of which 75% were members of the capsular polysaccharide (cps) operon. Among the non-cps genes identified as conditionally essential was relA, which encodes an enzyme whose activity is central to the bacterial stringent response—a conserved adaptation to environmental stress. We used blood coincubation studies of targeted knockout strains to confirm the expected growth defects of GBS deficient in capsule or stringent response activation. Unexpectedly, we found that the relA knockout strains demonstrated decreased expression of β-hemolysin/cytolysin, an important cytotoxin implicated in facilitating GBS invasion. Furthermore, chemical activation of the stringent response with serine hydroxamate increased β-hemolysin/cytolysin expression. To establish a mechanism by which the stringent response leads to increased cytotoxicity, we performed transcriptome sequencing (RNA-seq) on two GBS strains grown under stringent response or control conditions. This revealed a conserved decrease in the expression of genes in the arginine deiminase pathway during stringent response activation. Through coincubation with supplemental arginine and the arginine antagonist canavanine, we show that arginine availability is a determinant of GBS cytotoxicity and that the pathway between stringent response activation and increased virulence is arginine dependent.

A Comprehensive Assessment of the Genetic Determinants in Salmonella Typhimurium for Resistance to Hydrogen Peroxide Using Proteogenomics

Salmonella is an intracellular pathogen infecting a wide range of hosts and can survive in macrophages. An essential mechanism used by macrophages to eradicate Salmonella is production of reactive oxygen species. Here, we used proteogenomics to determine the candidate genes and proteins that have a role in resistance of S. Typhimurium to H₂O₂. For Tn-seq, a saturated Tn5 insertion library was grown in vitro under either 2.5 (H₂O₂L) or 3.5 mM H₂O₂ (H₂O₂H). We identified two sets of overlapping genes required for resistance of S. Typhimurium to H₂O₂L and H₂O₂H, and the results were validated via phenotypic evaluation of 50 selected mutants. The enriched pathways for H₂O₂ resistance included DNA repair, aromatic amino acid biosynthesis (aroBK), Fe-S cluster biosynthesis, iron homeostasis and a putative iron transporter system (ybbKLM), and H₂O₂ scavenging enzymes. Proteomics revealed that the majority of essential proteins, including ribosomal proteins, were downregulated upon exposure to H₂O₂. On the contrary, a subset of conditionally essential proteins identified by Tn-seq were analyzed by targeted proteomics, and 70% of them were upregulated by H₂O₂. The identified genes will deepen our understanding on S. Typhimurium survival mechanisms in macrophages, and can be exploited to develop new antimicrobial drugs.

Identification of Fitness Determinants during Energy-Limited Growth Arrest in Pseudomonas aeruginosa

Microbial growth arrest can be triggered by diverse factors, one of which is energy limitation due to scarcity of electron donors or acceptors. Genes that govern fitness during energy-limited growth arrest and the extent to which they overlap between different types of energy limitation are poorly defined. In this study, we exploited the fact that Pseudomonas aeruginosa can remain viable over several weeks when limited for organic carbon (pyruvate) as an electron donor or oxygen as an electron acceptor. ATP values were reduced under both types of limitation, yet more severely in the absence of oxygen. Using transposon-insertion sequencing (Tn-seq), we identified fitness determinants in these two energy-limited states. Multiple genes encoding general functions like transcriptional regulation and energy generation were required for fitness during carbon or oxygen limitation, yet many specific genes, and thus specific activities, differed in their relevance between these states. For instance, the global regulator RpoS was required during both types of energy limitation, while other global regulators such as DksA and LasR were required only during carbon or oxygen limitation, respectively. Similarly, certain ribosomal and tRNA modifications were specifically required during oxygen limitation. We validated fitness defects during energy limitation using independently generated mutants of genes detected in our screen. Mutants in distinct functional categories exhibited different fitness dynamics: regulatory genes generally manifested a phenotype early, whereas genes involved in cell wall metabolism were required later. Together, these results provide a new window into how P. aeruginosa survives growth arrest.

Growth-arrested bacteria are ubiquitous in nature and disease yet understudied at the molecular level. For example, growth-arrested cells constitute a major subpopulation of mature biofilms, serving as an antibiotic-tolerant reservoir in chronic infections. Identification of the genes required for survival of growth arrest (encompassing entry, maintenance, and exit) is an important first step toward understanding the physiology of bacteria in this state. Using Tn-seq, we identified and validated genes required for fitness of Pseudomonas aeruginosa when energy limited for organic carbon or oxygen, which represent two common causes of growth arrest for P. aeruginosa in diverse habitats. This unbiased, genome-wide survey is the first to reveal essential activities for a pathogen experiencing different types of energy limitation, finding both shared and divergent activities that are relevant at different survival stages. Future efforts can now be directed toward understanding how the biomolecules responsible for these activities contribute to fitness under these conditions.

A Comprehensive Overview of Online Resources to Identify and Predict Bacterial Essential Genes

Genes critical for the survival or reproduction of an organism in certain circumstances are classified as essential genes. Essential genes play a significant role in deciphering the survival mechanism of life. They may be greatly applied to pharmaceutics and synthetic biology. The continuous progress of experimental method for essential gene identification has accelerated the accumulation of gene essentiality data which facilitates the study of essential genes in silico. In this article, we present some available online resources related to gene essentiality, including bioinformatic software tools for transposon sequencing (Tn-seq) analysis, essential gene databases and online services to predict bacterial essential genes. We review several computational approaches that have been used to predict essential genes, and summarize the features used for gene essentiality prediction. In addition, we evaluate the available online bacterial essential gene prediction servers based on the experimentally validated essential gene sets of 30 bacteria from DEG. This article is intended to be a quick reference guide for the microbiologists interested in the essential genes.

Statistical analysis of genetic interactions in Tn-Seq data

Tn-Seq is an experimental method for probing the functions of genes through construction of complex random transposon insertion libraries and quantification of each mutant's abundance using next-generation sequencing. An important emerging application of Tn-Seq is for identifying genetic interactions, which involves comparing Tn mutant libraries generated in different genetic backgrounds (e.g. wild-type strain versus knockout strain). Several analytical methods have been proposed for analyzing Tn-Seq data to identify genetic interactions, including estimating relative fitness ratios and fitting a generalized linear model. However, these have limitations which necessitate an improved approach. We present a hierarchical Bayesian method for identifying genetic interactions through quantifying the statistical significance of changes in enrichment. The analysis involves a four-way comparison of insertion counts across datasets to identify transposon mutants that differentially affect bacterial fitness depending on genetic background. Our approach was applied to Tn-Seq libraries made in isogenic strains of Mycobacterium tuberculosis lacking three different genes of unknown function previously shown to be necessary for optimal fitness during infection. By analyzing the libraries subjected to selection in mice, we were able to distinguish several distinct classes of genetic interactions for each target gene that shed light on their functions and roles during infection.

Bacterial secretion of D-arginine controls environmental microbial biodiversity

Bacteria face tough competition in polymicrobial communities. To persist in a specific niche, many species produce toxic extracellular effectors to interfere with the growth of nearby microbes. These effectors include the recently reported non-canonical D-amino acids (NCDAAs). In Vibrio cholerae, the causative agent of cholera, NCDAAs control cell wall integrity in stationary phase. Here, an analysis of the composition of the extracellular medium of V. cholerae revealed the unprecedented presence of D-Arg. Compared with other D-amino acids, D-Arg displayed higher potency and broader toxicity in terms of the number of bacterial species affected. Tolerance to D-Arg was associated with mutations in the phosphate transport and chaperone systems, whereas D-Met lethality was suppressed by mutations in cell wall determinants. These observations suggest that NCDAAs target different cellular processes. Finally, even though virtually all Vibrio species are tolerant to D-Arg, only a few can produce this D-amino acid. Indeed, we demonstrate that D-Arg may function as part of a cooperative strategy in vibrio communities to protect non-producing members from competing bacteria. Because NCDAA production is widespread in bacteria, we anticipate that D-Arg is a relevant modulator of microbial subpopulations in diverse ecosystems.

Time-Resolved Transposon Insertion Sequencing Reveals Genome-Wide Fitness Dynamics during Infection

Transposon insertion sequencing (TIS) is a powerful high-throughput genetic technique that is transforming functional genomics in prokaryotes, because it enables genome-wide mapping of the determinants of fitness. However, current approaches for analyzing TIS data assume that selective pressures are constant over time and thus do not yield information regarding changes in the genetic requirements for growth in dynamic environments (e.g., during infection). Here, we describe structured analysis of TIS data collected as a time series, termed pattern analysis of conditional essentiality (PACE). From a temporal series of TIS data, PACE derives a quantitative assessment of each mutant’s fitness over the course of an experiment and identifies mutants with related fitness profiles. In so doing, PACE circumvents major limitations of existing methodologies, specifically the need for artificial effect size thresholds and enumeration of bacterial population expansion. We used PACE to analyze TIS samples of Edwardsiella piscicida (a fish pathogen) collected over a 2-week infection period from a natural host (the flatfish turbot). PACE uncovered more genes that affect E. piscicida’s fitness in vivo than were detected using a cutoff at a terminal sampling point, and it identified subpopulations of mutants with distinct fitness profiles, one of which informed the design of new live vaccine candidates. Overall, PACE enables efficient mining of time series TIS data and enhances the power and sensitivity of TIS-based analyses.

Transposon insertion sequencing (TIS) enables genome-wide mapping of the genetic determinants of fitness, typically based on observations at a single sampling point. Here, we move beyond analysis of endpoint TIS data to create a framework for analysis of time series TIS data, termed pattern analysis of conditional essentiality (PACE). We applied PACE to identify genes that contribute to colonization of a natural host by the fish pathogen Edwardsiella piscicida. PACE uncovered more genes that affect E. piscicida’s fitness in vivo than were detected using a terminal sampling point, and its clustering of mutants with related fitness profiles informed design of new live vaccine candidates. PACE yields insights into patterns of fitness dynamics and circumvents major limitations of existing methodologies. Finally, the PACE method should be applicable to additional “omic” time series data, including screens based on clustered regularly interspaced short palindromic repeats with Cas9 (CRISPR/Cas9).

Global Screening of Salmonella enterica Serovar Typhimurium Genes for Desiccation Survival

Salmonella spp., one of the most common foodborne bacterial pathogens, has the ability to survive under desiccation conditions in foods and food processing facilities for years. This raises the concerns of Salmonella infection in humans associated with low water activity foods. Salmonella responds to desiccation stress via complex pathways involving immediate physiological actions as well as coordinated genetic responses. However, the exact mechanisms of Salmonella to resist desiccation stress remain to be fully elucidated. In this study, we screened a genome-saturating transposon (Tn5) library of Salmonella Typhimurium (S. Typhimurium) 14028s under the in vitro desiccation stress using transposon sequencing (Tn-seq). We identified 61 genes and 6 intergenic regions required to overcome desiccation stress. Salmonella desiccation resistance genes were mostly related to energy production and conversion; cell wall/membrane/envelope biogenesis; inorganic ion transport and metabolism; regulation of biological process; DNA metabolic process; ABC transporters; and two component system. More than 20% of the Salmonella desiccation resistance genes encode either putative or hypothetical proteins. Phenotypic evaluation of 12 single gene knockout mutants showed 3 mutants (atpH, atpG, and corA) had significantly (p < 0.02) reduced survival as compared to the wild type during desiccation survival. Thus, our study provided new insights into the molecular mechanisms utilized by Salmonella for survival against desiccation stress. The findings might be further exploited to develop effective control strategies against Salmonella contamination in low water activity foods and food processing facilities.

Essential genome of Campylobacter jejuni

BACKGROUND

Campylobacter species are a leading cause of bacterial foodborne illness worldwide. Despite the global efforts to curb them, Campylobacter infections have increased continuously in both developed and developing countries. The development of effective strategies to control the infection by this pathogen is warranted. The essential genes of bacteria are the most prominent targets for this purpose. In this study, we used transposon sequencing (Tn-seq) of a genome-saturating library of Tn5 insertion mutants to define the essential genome of C. jejuni at a high resolution.

RESULT

We constructed a Tn5 mutant library of unprecedented complexity in C. jejuni NCTC 11168 with 95,929 unique insertions throughout the genome and used the genomic DNA of the library for the reconstruction of Tn5 libraries in the same (C. jejuni NCTC 11168) and different strain background (C. jejuni 81-176) through natural transformation. We identified 166 essential protein-coding genes and 20 essential transfer RNAs (tRNA) in C. jejuni NCTC 11168 which were intolerant to Tn5 insertions during in vitro growth. The reconstructed C. jejuni 81-176 library had 384 protein coding genes with no Tn5 insertions. Essential genes in both strain backgrounds were highly enriched in the cluster of orthologous group (COG) categories of 'Translation, ribosomal structure and biogenesis (J)', 'Energy production and conversion (C)', and 'Coenzyme transport and metabolism (H)'.

CONCLUSION

Comparative analysis among this and previous studies identified 50 core essential genes of C. jejuni, which can be further investigated for the development of novel strategies to control the spread of this notorious foodborne bacterial pathogen.

The essential genomic landscape of the commensal Bifidobacterium breve UCC2003

Bifidobacteria are common gut commensals with purported health-promoting effects. This has encouraged scientific research into bifidobacteria, though recalcitrance to genetic manipulation and scarcity of molecular tools has hampered our knowledge on the precise molecular determinants of their health-promoting attributes and gut adaptation. To overcome this problem and facilitate functional genomic analyses in bifidobacteria, we created a large Tn5 transposon mutant library of the commensal Bifidobacterium breve UCC2003 that was further characterized by means of a Transposon Directed Insertion Sequencing (TraDIS) approach. Statistical analysis of transposon insertion distribution revealed a set of 453 genes that are essential for or markedly contribute to growth of this strain under laboratory conditions. These essential genes encode functions involved in the so-called bifid-shunt, most enzymes related to nucleotide biosynthesis and a range of housekeeping functions. Comparison to the Bifidobacterium and B. breve core genomes highlights a high degree of conservation of essential genes at the species and genus level, while comparison to essential gene datasets from other gut bacteria identified essential genes that appear specific to bifidobacteria. This work establishes a useful molecular tool for scientific discovery of bifidobacteria and identifies targets for further studies aimed at characterizing essential functions not previously examined in bifidobacteria.

TnseqDiff: identification of conditionally essential genes in transposon sequencing studies

BACKGROUND

Tn-Seq is a high throughput technique for analysis of transposon mutant libraries to determine conditional essentiality of a gene under an experimental condition. A special feature of the Tn-seq data is that multiple mutants in a gene provides independent evidence to prioritize that gene as being essential. The existing methods do not account for this feature or rely on a high-density transposon library. Moreover, these methods are unable to accommodate complex designs.

RESULTS

The method proposed here is specifically designed for the analysis of Tn-Seq data. It utilizes two steps to estimate the conditional essentiality for each gene in the genome. First, it collects evidence of conditional essentiality for each insertion by comparing read counts of that insertion between conditions. Second, it combines insertion-level evidence for the corresponding gene. It deals with data from both low- and high-density transposon libraries and accommodates complex designs. Moreover, it is very fast to implement. The performance of the proposed method was tested on simulated data and experimental Tn-Seq data from Serratia marcescens transposon mutant library used to identify genes that contribute to fitness in a murine model of infection.

CONCLUSION

We describe a new, efficient method for identifying conditionally essential genes in Tn-Seq experiments with high detection sensitivity and specificity. It is implemented as TnseqDiff function in R package Tnseq and can be installed from the Comprehensive R Archive Network, CRAN.

Combining Genome-Scale Experimental and Computational Methods To Identify Essential Genes in Rhodobacter sphaeroides

Knowledge about the role of genes under a particular growth condition is required for a holistic understanding of a bacterial cell and has implications for health, agriculture, and biotechnology. We developed the Tn-seq analysis software (TSAS) package to provide a flexible and statistically rigorous workflow for the high-throughput analysis of insertion mutant libraries, advanced the knowledge of gene essentiality in R. sphaeroides, and illustrated how Tn-seq data can be used to more accurately identify genes that play important roles in metabolism and other processes that are essential for cellular survival.

Rhodobacter sphaeroides is one of the best-studied alphaproteobacteria from biochemical, genetic, and genomic perspectives. To gain a better systems-level understanding of this organism, we generated a large transposon mutant library and used transposon sequencing (Tn-seq) to identify genes that are essential under several growth conditions. Using newly developed Tn-seq analysis software (TSAS), we identified 493 genes as essential for aerobic growth on a rich medium. We then used the mutant library to identify conditionally essential genes under two laboratory growth conditions, identifying 85 additional genes required for aerobic growth in a minimal medium and 31 additional genes required for photosynthetic growth. In all instances, our analyses confirmed essentiality for many known genes and identified genes not previously considered to be essential. We used the resulting Tn-seq data to refine and improve a genome-scale metabolic network model (GEM) for R. sphaeroides. Together, we demonstrate how genetic, genomic, and computational approaches can be combined to obtain a systems-level understanding of the genetic framework underlying metabolic diversity in bacterial species.

IMPORTANCE Knowledge about the role of genes under a particular growth condition is required for a holistic understanding of a bacterial cell and has implications for health, agriculture, and biotechnology. We developed the Tn-seq analysis software (TSAS) package to provide a flexible and statistically rigorous workflow for the high-throughput analysis of insertion mutant libraries, advanced the knowledge of gene essentiality in R. sphaeroides, and illustrated how Tn-seq data can be used to more accurately identify genes that play important roles in metabolism and other processes that are essential for cellular survival.

Author Video: An author video summary of this article is available.

Transposon insertion sequencing reveals T4SS as the major genetic trait for conjugation transfer of multi-drug resistance pEIB202 from Edwardsiella

Background

Conjugation is a major type of horizontal transmission of genes that involves transfer of a plasmid into a recipient using specific conjugation machinery, which results in an extended spectrum of bacterial antibiotics resistance. However, there is inadequate knowledge about the regulator and mechanisms that control the conjugation processes, especially in an aquaculture environment where a cocktail of antibiotics may be present. Here, we investigated these with pEIB202, a typical multi-drug resistant IncP plasmid encoding tetracycline, streptomycin, sulfonamide and chloramphenicol resistance in fish pathogen Edwardsiella piscicida strain EIB202.

Results

We used transposon insertion sequencing (TIS) to identify genes that are responsible for conjugation transfer of pEIB202. All ten of the plasmid-borne type IV secretion system (T4SS) genes and a putative lipoprotein p007 were identified to play an important role in pEIB202 horizontal transfer. Antibiotics appear to modulate conjugation frequencies by repressing T4SS gene expression. In addition, we identified topA gene, which encodes topoisomerase I, as an inhibitor of pEIB202 transfer. Furthermore, the RNA-seq analysis of the response regulator EsrB encoded on the chromosome also revealed its essential role in facilitating the conjugation by upregulating the T4SS genes.

Conclusions

Collectively, our screens unraveled the genetic basis of the conjugation transfer of pEIB202 and the influence of horizontally acquired EsrB on this process. Our results will improve the understanding of the mechanism of plasmid conjugation processes that facilitate dissemination of antibiotic resistance especially in aquaculture industries.

Phenotypic difference between Δ(srl-recA)306 and ΔrecA::Km elucidated by next-generation sequencing combined with a long-PCR system

Many significant gene mutations in E. coli have contributed to the development of genetics. Among these, a commonly used recA mutation, Δ(srl-recA)306 has been sequenced by a next-generation sequencer combined with a long PCR. An original report described that Δ(srl-recA)306 cells were deleted from srlR to recA genes in their genome. The next-generation sequencer enables more accurate details to be determined. We ask whether both surrounding genes from hypF to norV for srlR and alaS for recA is there first. The long PCR was carried out with primers, norR and alaS, and amplified DNA fragments differed in length from wild to Δ(srl-recA)306 cells, suggesting that an entire Δ(srl-recA)306 mutation was included. Sequences of those DNA fragments indicated that 9147 bp, from srlR to recA including 10 genes, were replaced by a Tn10 DNA sequence. Junction points at both srlR-Tn10 and Tn10-recA were determined precisely. The results indicate that the first 97% of recA gene sequences were lost with a downstream recX gene remaining intact. The phenotypic difference between Δ(srl-recA)306 and ΔrecA::Km is discussed.

Genetic interaction mapping with microfluidic-based single cell sequencing

Genetic interaction mapping is useful for understanding the molecular basis of cellular decision making, but elucidating interactions genome-wide is challenging due to the massive number of gene combinations that must be tested. Here, we demonstrate a simple approach to thoroughly map genetic interactions in bacteria using microfluidic-based single cell sequencing. Using single cell PCR in droplets, we link distinct genetic information into single DNA sequences that can be decoded by next generation sequencing. Our approach is scalable and theoretically enables the pooling of entire interaction libraries to interrogate multiple pairwise genetic interactions in a single culture. The speed, ease, and low-cost of our approach makes genetic interaction mapping viable for routine characterization, allowing the interaction network to be used as a universal read out for a variety of biology experiments, and for the elucidation of interaction networks in non-model organisms.

A cocktail of three virulent bacteriophages prevents Vibrio cholerae infection in animal models

Effective prevention strategies will be essential in reducing disease burden due to bacterial infections. Here we harness the specificity and rapid-acting properties of bacteriophages as a potential prophylaxis therapy for cholera, a severely dehydrating disease caused by Vibrio cholerae. To this end, we test a cocktail of three virulent phages in two animal models of cholera pathogenesis (infant mouse and rabbit models). Oral administration of the phages up to 24 h before V. cholerae challenge reduces colonization of the intestinal tract and prevents cholera-like diarrhea. None of the surviving V. cholerae colonies are resistant to all three phages. Genome sequencing and variant analysis of the surviving colonies indicate that resistance to the phages is largely conferred by mutations in genes required for the production of the phage receptors. For acute infections, such as cholera, phage prophylaxis could provide a strategy to limit the impact of bacterial disease on human health.

Comprehensive Essentiality Analysis of the Mycobacterium tuberculosis Genome via Saturating Transposon Mutagenesis

For decades, identifying the regions of a bacterial chromosome that are necessary for viability has relied on mapping integration sites in libraries of random transposon mutants to find loci that are unable to sustain insertion. To date, these studies have analyzed subsaturated libraries, necessitating the application of statistical methods to estimate the likelihood that a gap in transposon coverage is the result of biological selection and not the stochasticity of insertion. As a result, the essentiality of many genomic features, particularly small ones, could not be reliably assessed. We sought to overcome this limitation by creating a completely saturated transposon library in Mycobacterium tuberculosis In assessing the composition of this highly saturated library by deep sequencing, we discovered that a previously unknown sequence bias of the Himar1 element rendered approximately 9% of potential TA dinucleotide insertion sites less permissible for insertion. We used a hidden Markov model of essentiality that accounted for this unanticipated bias, allowing us to confidently evaluate the essentiality of features that contained as few as 2 TA sites, including open reading frames (ORF), experimentally identified noncoding RNAs, methylation sites, and promoters. In addition, several essential regions that did not correspond to known features were identified, suggesting uncharacterized functions that are necessary for growth. This work provides an authoritative catalog of essential regions of the M. tuberculosis genome and a statistical framework for applying saturating mutagenesis to other bacteria.

IMPORTANCE

Sequencing of transposon-insertion mutant libraries has become a widely used tool for probing the functions of genes under various conditions. The Himar1 transposon is generally believed to insert with equal probabilities at all TA dinucleotides, and therefore its absence in a mutant library is taken to indicate biological selection against the corresponding mutant. Through sequencing of a saturated Himar1 library, we found evidence that TA dinucleotides are not equally permissive for insertion. The insertion bias was observed in multiple prokaryotes and influences the statistical interpretation of transposon insertion (TnSeq) data and characterization of essential genomic regions. Using these insights, we analyzed a fully saturated TnSeq library for M. tuberculosis, enabling us to generate a comprehensive catalog of in vitro essentiality, including ORFs smaller than those found in any previous study, small (noncoding) RNAs (sRNAs), promoters, and other genomic features.

Genome-Wide Transposon Mutagenesis in Mycobacterium tuberculosis and Mycobacterium smegmatis

TnSeq, or transposon (Tn) insertion sequencing, is a powerful method for identifying the essential-as well as conditionally essential-regions in a genome, both coding and noncoding. The advent of accessible massively parallel DNA sequencing technologies in particular has resulted in the increased use of TnSeq-based approaches to elucidate various aspects of bacterial physiology and metabolism. Moreover, the availability of detailed protocols has enabled even nonspecialist laboratories to adapt and develop TnSeq approaches to address specific research questions. In this chapter, we describe a recently modified experimental protocol used in our laboratory for TnSeq in the major human pathogen, Mycobacterium tuberculosis, as well as the related non-pathogenic mycobacterium, M. smegmatis. The method, which was developed in close consultation with pioneers in the field of mycobacterial genetics, includes the steps involved in preparing a phage stock, generating a mutant library, selection of the library under a specific experimental condition, isolation of genomic DNA from the pooled population of mutants, amplification of the sites of Tn insertion and, finally, determining the essential genomic regions by next-generation sequencing.

PIMMS (Pragmatic Insertional Mutation Mapping System) Laboratory Methodology a Readily Accessible Tool for Identification of Essential Genes in Streptococcus

The Pragmatic Insertional Mutation Mapping (PIMMS) laboratory protocol was developed alongside various bioinformatics packages (Blanchard et al., 2015) to enable detection of essential and conditionally essential genes in Streptococcus and related bacteria. This extended the methodology commonly used to locate insertional mutations in individual mutants to the analysis of mutations in populations of bacteria. In Streptococcus uberis, a pyogenic Streptococcus associated with intramammary infection and mastitis in ruminants, the mutagen pGhost9:ISS1 was shown to integrate across the entire genome. Analysis of >80,000 mutations revealed 196 coding sequences, which were not be mutated and a further 67 where mutation only occurred beyond the 90th percentile of the coding sequence. These sequences showed good concordance with sequences within the database of essential genes and typically matched sequences known to be associated with basic cellular functions. Due to the broad utility of this mutagen and the simplicity of the methodology it is anticipated that PIMMS will be of value to a wide range of laboratories in functional genomic analysis of a wide range of Gram positive bacteria (Streptococcus, Enterococcus, and Lactococcus) of medical, veterinary, and industrial significance.

No-Go Zones for Mariner Transposition

The property of transposons to randomly insert into target DNA has long been exploited for generalized mutagenesis and forward genetic screens. Newer applications that monitor the relative abundance of each transposon insertion in large libraries of mutants can be used to evaluate the roles in cellular fitness of all genes of an organism, provided that transposition is in fact random across all genes. In a recent article, Kimura and colleagues identified an important exception to the latter assumption [S. Kimura, T. P. Hubbard, B. M. Davis, M. K. Waldor, mBio 7(4):e01351-16, 2016, doi:10.1128/mBio.01351-16]. They provide evidence that the Mariner transposon exhibits locus-specific site preferences in the presence of the histone-like nucleoid structuring protein H-NS. This effect was shown to bias results for important virulence loci in Vibrio cholerae and to result in misidentification of genes involved in growth in vitro. Fortunately, the bulk of this bacterium’s genome was unaffected by this bias, and recognizing the H-NS effect allows filtering to improve the accuracy of the results.

Combining Shigella Tn-seq data with gold-standard E. coli gene deletion data suggests rare transitions between essential and non-essential gene functionality

BACKGROUND

Gene essentiality - whether or not a gene is necessary for cell growth - is a fundamental component of gene function. It is not well established how quickly gene essentiality can change, as few studies have compared empirical measures of essentiality between closely related organisms.

RESULTS

Here we present the results of a Tn-seq experiment designed to detect essential protein coding genes in the bacterial pathogen Shigella flexneri 2a 2457T on a genome-wide scale. Superficial analysis of this data suggested that 481 protein-coding genes in this Shigella strain are critical for robust cellular growth on rich media. Comparison of this set of genes with a gold-standard data set of essential genes in the closely related Escherichia coli K12 BW25113 revealed that an excessive number of genes appeared essential in Shigella but non-essential in E. coli. Importantly, and in converse to this comparison, we found no genes that were essential in E. coli and non-essential in Shigella, implying that many genes were artefactually inferred as essential in Shigella. Controlling for such artefacts resulted in a much smaller set of discrepant genes. Among these, we identified three sets of functionally related genes, two of which have previously been implicated as critical for Shigella growth, but which are dispensable for E. coli growth.

CONCLUSIONS

The data presented here highlight the small number of protein coding genes for which we have strong evidence that their essentiality status differs between the closely related bacterial taxa E. coli and Shigella. A set of genes involved in acetate utilization provides a canonical example. These results leave open the possibility of developing strain-specific antibiotic treatments targeting such differentially essential genes, but suggest that such opportunities may be rare in closely related bacteria.

The Nucleoid Binding Protein H-NS Biases Genome-Wide Transposon Insertion Landscapes

Transposon insertion sequencing (TIS; also known as TnSeq) is a potent approach commonly used to comprehensively define the genetic loci that contribute to bacterial fitness in diverse environments. A key presumption underlying analyses of TIS datasets is that loci with a low frequency of transposon insertions contribute to fitness. However, it is not known whether factors such as nucleoid binding proteins can alter the frequency of transposon insertion and thus whether TIS output may systematically reflect factors that are independent of the role of the loci in fitness. Here, we investigated whether the histone-like nucleoid structuring (H-NS) protein, which preferentially associates with AT-rich sequences, modulates the frequency of Mariner transposon insertion in the Vibrio cholerae genome, using comparative analysis of TIS results from wild-type (wt) and Δhns V. cholerae strains. These analyses were overlaid on gene classification based on GC content as well as on extant genome-wide identification of H-NS binding loci. Our analyses revealed a significant dearth of insertions within AT-rich loci in wt V. cholerae that was not apparent in the Δhns insertion library. Additionally, we observed a striking correlation between genetic loci that are overrepresented in the Δhns insertion library relative to their insertion frequency in wt V. cholerae and loci previously found to physically interact with H-NS. Collectively, our findings reveal that factors other than genetic fitness can systematically modulate the frequency of transposon insertions in TIS studies and add a cautionary note to interpretation of TIS data, particularly for AT-rich sequences.

IMPORTANCE

Transposon insertion sequencing (TIS) is often used to assess the relative frequency with which genetic loci can be disrupted, which is taken as an indicator of their importance for bacterial fitness. Here, we report that biological factors other than the relative levels of fitness of insertion mutants can influence TIS output. We found that the presence of the DNA binding protein H-NS, which preferentially recognizes AT-rich sequences, is linked to significant underrepresentation of mutations within AT-rich loci in transposon insertion libraries. Furthermore, there is a marked correspondence between loci bound by H-NS and loci with an increased frequency of disruption in a Δhns insertion library relative to a wt library. Our data suggest that factors other than genetic fitness (e.g., DNA binding proteins such as H-NS) can systematically modulate the frequency of transposon insertions in TIS studies and add a note of caution for interpretation of TIS data.

Multidrug Intrinsic Resistance Factors in Staphylococcus aureus Identified by Profiling Fitness within High-Diversity Transposon Libraries

Staphylococcus aureus is a leading cause of life-threatening infections worldwide. The MIC of an antibiotic against S. aureus, as well as other microbes, is determined by the affinity of the antibiotic for its target in addition to a complex interplay of many other cellular factors. Identifying nontarget factors impacting resistance to multiple antibiotics could inform the design of new compounds and lead to more-effective antimicrobial strategies. We examined large collections of transposon insertion mutants in S. aureus using transposon sequencing (Tn-Seq) to detect transposon mutants with reduced fitness in the presence of six clinically important antibiotics—ciprofloxacin, daptomycin, gentamicin, linezolid, oxacillin, and vancomycin. This approach allowed us to assess the relative fitness of many mutants simultaneously within these libraries. We identified pathways/genes previously known to be involved in resistance to individual antibiotics, including graRS and vraFG (graRS/vraFG), mprF, and fmtA, validating the approach, and found several to be important across multiple classes of antibiotics. We also identified two new, previously uncharacterized genes, SAOUHSC_01025 and SAOUHSC_01050, encoding polytopic membrane proteins, as important in limiting the effectiveness of multiple antibiotics. Machine learning identified similarities in the fitness profiles of graXRS/vraFG, SAOUHSC_01025, and SAOUHSC_01050 mutants upon antibiotic treatment, connecting these genes of unknown function to modulation of crucial cell envelope properties. Therapeutic strategies that combine a known antibiotic with a compound that targets these or other intrinsic resistance factors may be of value for enhancing the activity of existing antibiotics for treating otherwise-resistant S. aureus strains.

Bacterial resistance to every major class of antibiotics has emerged, and we are entering a “post-antibiotic era” where relatively minor infections can lead to serious complications or even death. The utility of an antibiotic for a specific pathogen is limited by both intrinsic and acquired factors. Identifying the repertoire of intrinsic resistance factors of an antibiotic for Staphylococcus aureus, a leading cause of community- and hospital-acquired infections, would inform the design of new drugs as well as the identification of compounds that enhance the activity of existing drugs. To identify factors that limit the activity of antibiotics against S. aureus, we used Tn-Seq to simultaneously assess fitness of transposon mutants in every nonessential gene in the presence of six clinically important antibiotics. This work provides an efficient approach for identifying promising targets for drugs that can enhance susceptibility or restore sensitivity to existing antibiotics.