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

A Transposon Screen Identifies Genetic Determinants of Vibrio cholerae Resistance to High-Molecular-Weight Antibiotics

Gram-negative bacteria are notoriously resistant to a variety of high-molecular-weight antibiotics due to the limited permeability of their outer membrane (OM). The basis of OM barrier function and the genetic factors required for its maintenance remain incompletely understood. Here, we employed transposon insertion sequencing to identify genes required for Vibrio cholerae resistance to vancomycin and bacitracin, antibiotics that are thought to be too large to efficiently penetrate the OM. The screen yielded several genes whose protein products are predicted to participate in processes important for OM barrier functions and for biofilm formation. In addition, we identified a novel factor, designated vigA (for vancomycin inhibits growth), that has not previously been characterized or linked to outer membrane function. The vigA open reading frame (ORF) codes for an inner membrane protein, and in its absence, cells became highly sensitive to glycopeptide antibiotics (vancomycin and ramoplanin) and bacitracin but not to other large antibiotics or detergents. In contrast to wild-type (WT) cells, the vigA mutant was stained with fluorescent vancomycin. These observations suggest that VigA specifically prevents the periplasmic accumulation of certain large antibiotics without exerting a general role in the maintenance of OM integrity. We also observed marked interspecies variability in the susceptibilities of Gram-negative pathogens to glycopeptides and bacitracin. Collectively, our findings suggest that the OM barrier is not absolute but rather depends on specific OM-antibiotic interactions.

The essential genome of Streptococcus agalactiae

BACKGROUND

Next-generation sequencing of transposon-genome junctions from a saturated bacterial mutant library (Tn-seq) is a powerful tool that permits genome-wide determination of the contribution of genes to fitness of the organism under a wide range of experimental conditions. We report development, testing, and results from a Tn-seq system for use in Streptococcus agalactiae (group B Streptococcus; GBS), an important cause of neonatal sepsis.

METHODS

Our method uses a Himar1 mini-transposon that inserts at genomic TA dinucleotide sites, delivered to GBS on a temperature-sensitive plasmid that is subsequently cured from the bacterial population. In order to establish the GBS essential genome, we performed Tn-seq on DNA collected from three independent mutant libraries-with at least 135,000 mutants per library-at serial 24 h time points after outgrowth in rich media.

RESULTS

After statistical analysis of transposon insertion density and distribution, we identified 13.5 % of genes as essential and 1.2 % as critical, with high levels of reproducibility. Essential and critical genes are enriched for fundamental cellular housekeeping functions, such as acyl-tRNA biosynthesis, nucleotide metabolism, and glycolysis. We further validated our system by comparing fitness assignments of homologous genes in GBS and a close bacterial relative, Streptococcus pyogenes, which demonstrated 93 % concordance. Finally, we used our fitness assignments to identify signal transduction pathway components predicted to be essential or critical in GBS.

CONCLUSIONS

We believe that our baseline fitness assignments will be a valuable tool for GBS researchers and that our system has the potential to reveal key pathogenesis gene networks and potential therapeutic/preventative targets.

The design and analysis of transposon insertion sequencing experiments

Transposon insertion sequencing (TIS) is a powerful approach that can be extensively applied to the genome-wide definition of loci that are required for bacterial growth under diverse conditions. However, experimental design choices and stochastic biological processes can heavily influence the results of TIS experiments and affect downstream statistical analysis. In this Opinion article, we discuss TIS experimental parameters and how these factors relate to the benefits and limitations of the various statistical frameworks that can be applied to the computational analysis of TIS data.

Genetic analysis of Vibrio parahaemolyticus intestinal colonization

Vibrio parahaemolyticus is the most common cause of seafood-borne gastroenteritis worldwide and a blight on global aquaculture. This organism requires a horizontally acquired type III secretion system (T3SS2) to infect the small intestine, but knowledge of additional factors that underlie V. parahaemolyticus pathogenicity is limited. We used transposon-insertion sequencing to screen for genes that contribute to viability of V. parahaemolyticus in vitro and in the mammalian intestine. Our analysis enumerated and controlled for the host infection bottleneck, enabling robust assessment of genetic contributions to in vivo fitness. We identified genes that contribute to V. parahaemolyticus colonization of the intestine independent of known virulence mechanisms in addition to uncharacterized components of T3SS2. Our study revealed that toxR, an ancestral locus in Vibrio species, is required for V. parahaemolyticus fitness in vivo and for induction of T3SS2 gene expression. The regulatory mechanism by which V. parahaemolyticus ToxR activates expression of T3SS2 resembles Vibrio cholerae ToxR regulation of distinct virulence elements acquired via lateral gene transfer. Thus, disparate horizontally acquired virulence systems have been placed under the control of this ancestral transcription factor across independently evolved human pathogens.

A zero-inflated Poisson model for insertion tolerance analysis of genes based on Tn-seq data

MOTIVATION

Transposon insertion sequencing (Tn-seq) is an emerging technology that combines transposon mutagenesis with next-generation sequencing technologies for the identification of genes related to bacterial survival. The resulting data from Tn-seq experiments consist of sequence reads mapped to millions of potential transposon insertion sites and a large portion of insertion sites have zero mapped reads. Novel statistical method for Tn-seq data analysis is needed to infer functions of genes on bacterial growth.

RESULTS

In this article, we propose a zero-inflated Poisson model for analyzing the Tn-seq data that are high-dimensional and with an excess of zeros. Maximum likelihood estimates of model parameters are obtained using an expectation-maximization (EM) algorithm, and pseudogenes are utilized to construct appropriate statistical tests for the transposon insertion tolerance of normal genes of interest. We propose a multiple testing procedure that categorizes genes into each of the three states, hypo-tolerant, tolerant and hyper-tolerant, while controlling the false discovery rate. We evaluate the proposed method with simulation studies and apply the proposed method to a real Tn-seq data from an experiment that studied the bacterial pathogen, Campylobacter jejuniAvailability and implementation: We provide R code for implementing our proposed method at http://github.com/ffliu/TnSeq A user's guide with example data analysis is also available there.

CONTACT

pliu@iastate.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Normalization of transposon-mutant library sequencing datasets to improve identification of conditionally essential genes

Sequencing of transposon-mutant libraries using next-generation sequencing (TnSeq) has become a popular method for determining which genes and non-coding regions are essential for growth under various conditions in bacteria. For methods that rely on quantitative comparison of counts of reads at transposon insertion sites, proper normalization of TnSeq datasets is vitally important. Real TnSeq datasets are often noisy and exhibit a significant skew that can be dominated by high counts at a small number of sites (often for non-biological reasons). If two datasets that are not appropriately normalized are compared, it might cause the artifactual appearance of Differentially Essential (DE) genes in a statistical test, constituting type I errors (false positives). In this paper, we propose a novel method for normalization of TnSeq datasets that corrects for the skew of read-count distributions by fitting them to a Beta-Geometric distribution. We show that this read-count correction procedure reduces the number of false positives when comparing replicate datasets grown under the same conditions (for which no genuine differences in essentiality are expected). We compare these results to results obtained with other normalization procedures, and show that it results in greater reduction in the number of false positives. In addition we investigate the effects of normalization on the detection of DE genes.

A Cytosine Methyltransferase Modulates the Cell Envelope Stress Response in the Cholera Pathogen [corrected]

DNA methylation is a key epigenetic regulator in all domains of life, yet the effects of most bacterial DNA methyltransferases on cellular processes are largely undefined. Here, we used diverse techniques, including bisulfite sequencing, transcriptomics, and transposon insertion site sequencing to extensively characterize a 5-methylcytosine (5mC) methyltransferase, VchM, in the cholera pathogen, Vibrio cholerae. We have comprehensively defined VchM's DNA targets, its genetic interactions and the gene networks that it regulates. Although VchM is a relatively new component of the V. cholerae genome, it is required for optimal V. cholerae growth in vitro and during infection. Unexpectedly, the usually essential σE cell envelope stress pathway is dispensable in ∆vchM V. cholerae, likely due to its lower activation in this mutant and the capacity for VchM methylation to limit expression of some cell envelope modifying genes. Our work illuminates how an acquired DNA methyltransferase can become integrated within complex cell circuits to control critical housekeeping processes.

Identification of Population Bottlenecks and Colonization Factors during Assembly of Bacterial Communities within the Zebrafish Intestine

The zebrafish, Danio rerio, is a powerful model for studying bacterial colonization of the vertebrate intestine, but the genes required by commensal bacteria to colonize the zebrafish gut have not yet been interrogated on a genome-wide level. Here we apply a high-throughput transposon mutagenesis screen to Aeromonas veronii Hm21 and Vibrio sp. strain ZWU0020 during their colonization of the zebrafish intestine alone and in competition with each other, as well as in different colonization orders. We use these transposon-tagged libraries to track bacterial population sizes in different colonization regimes and to identify gene functions required during these processes. We show that intraspecific, but not interspecific, competition with a previously established bacterial population greatly reduces the ability of these two bacterial species to colonize. Further, using a simple binomial sampling model, we show that under conditions of interspecific competition, genes required for colonization cannot be identified because of the population bottleneck experienced by the second colonizer. When bacteria colonize the intestine alone or at the same time as the other species, we find shared suites of functional requirements for colonization by the two species, including a prominent role for chemotaxis and motility, regardless of the presence of another species.

Zebrafish larvae, which are amenable to large-scale gnotobiotic studies, comprehensive sampling of their intestinal microbiota, and live imaging, are an excellent model for investigations of vertebrate intestinal colonization dynamics. We sought to develop a mutagenesis and tagging system in order to understand bacterial population dynamics and functional requirements during colonization of the larval zebrafish intestine. We explored changes in bacterial colonization dynamics and functional requirements when bacteria colonize a bacterium-free intestine, one previously colonized by their own species, or one colonized previously or simultaneously with a different species. This work provides a framework for rapid identification of colonization factors important under different colonization conditions. Furthermore, we demonstrate that when colonizing bacterial populations are very small, this approach is not accurate because random sampling of the input pool is sufficient to explain the distribution of inserts recovered from bacteria that colonized the intestines.

TRANSIT--A Software Tool for Himar1 TnSeq Analysis

TnSeq has become a popular technique for determining the essentiality of genomic regions in bacterial organisms. Several methods have been developed to analyze the wealth of data that has been obtained through TnSeq experiments. We developed a tool for analyzing Himar1 TnSeq data called TRANSIT. TRANSIT provides a graphical interface to three different statistical methods for analyzing TnSeq data. These methods cover a variety of approaches capable of identifying essential genes in individual datasets as well as comparative analysis between conditions. We demonstrate the utility of this software by analyzing TnSeq datasets of M. tuberculosis grown on glycerol and cholesterol. We show that TRANSIT can be used to discover genes which have been previously implicated for growth on these carbon sources. TRANSIT is written in Python, and thus can be run on Windows, OSX and Linux platforms. The source code is distributed under the GNU GPL v3 license and can be obtained from the following GitHub repository: https://github.com/mad-lab/transit.

Peptidoglycan synthesis in Mycobacterium tuberculosis is organized into networks with varying drug susceptibility

Peptidoglycan (PG), a complex polymer composed of saccharide chains cross-linked by short peptides, is a critical component of the bacterial cell wall. PG synthesis has been extensively studied in model organisms but remains poorly understood in mycobacteria, a genus that includes the important human pathogen Mycobacterium tuberculosis (Mtb). The principle PG synthetic enzymes have similar and, at times, overlapping functions. To determine how these are functionally organized, we carried out whole-genome transposon mutagenesis screens in Mtb strains deleted for ponA1, ponA2, and ldtB, major PG synthetic enzymes. We identified distinct factors required to sustain bacterial growth in the absence of each of these enzymes. We find that even the homologs PonA1 and PonA2 have unique sets of genetic interactions, suggesting there are distinct PG synthesis pathways in Mtb. Either PonA1 or PonA2 is required for growth of Mtb, but both genetically interact with LdtB, which has its own distinct genetic network. We further provide evidence that each interaction network is differentially susceptible to antibiotics. Thus, Mtb uses alternative pathways to produce PG, each with its own biochemical characteristics and vulnerabilities.

A D, D-carboxypeptidase is required for Vibrio cholerae halotolerance

The biological roles of low molecular weight penicillin-binding proteins (LMW PBP) have been difficult to discern in Gram-negative organisms. In Escherichia coli, mutants lacking these proteins often have no phenotype, and cells lacking all seven LMW PBPs remain viable. In contrast, we report here that Vibrio cholerae lacking DacA-1, a PBP5 homologue, displays slow growth, aberrant morphology and altered peptidoglycan (PG) homeostasis in Luria-Bertani (LB) medium, as well as a profound plating defect. DacA-1 alone among V. cholerae's LMW PBPs is critical for bacterial growth; mutants lacking the related protein DacA-2 and/or homologues of PBP4 or PBP7 displayed normal growth and morphology. Remarkably, the growth and morphology of the dacA-1 mutant were unimpaired in LB media containing reduced concentrations of NaCl (100 mM or less), and also within suckling mice, a model host for the study of cholera pathogenesis. Peptidoglycan from the dacA-1 mutant contained elevated pentapeptide levels in standard and low salt media, and comparative analyses suggest that DacA-1 is V. cholerae's principal DD-carboxypeptidase. The basis for the dacA-1 mutant's halosensitivity is unknown; nonetheless, the mutant's survival in biochemically uncharacterized environments (such as the suckling mouse intestine) can be used as a reporter of low Na(+) content.

Infant Rabbit Model for Diarrheal Diseases

Vibrio cholerae is the agent of cholera, a potentially lethal diarrheal disease that remains a significant threat to populations in developing nations. The infant rabbit model of cholera is the only non-surgical small animal model system that closely mimics human cholera. Following orogastric inoculation, V. cholerae colonizes the intestines of infant rabbits, and the animals develop severe cholera-like diarrhea. In this unit, we provide a detailed description of the preparation of the V. cholerae inoculum, the inoculation process and the collection and processing of tissue samples. This infection model is useful for studies of V. cholerae factors and mechanisms that promote its intestinal colonization and enterotoxicity, as well as the host response to infection. The infant rabbit model of cholera enables investigations that will further our understanding of the pathophysiology of cholera and provides a platform for testing new therapeutics.

A Genome-Wide Screen Reveals that the Vibrio cholerae Phosphoenolpyruvate Phosphotransferase System Modulates Virulence Gene Expression

Diverse environmental stimuli and a complex network of regulatory factors are known to modulate expression of Vibrio cholerae's principal virulence factors. However, there is relatively little known about how metabolic factors impinge upon the pathogen's well-characterized cascade of transcription factors that induce expression of cholera toxin and the toxin-coregulated pilus (TCP). Here, we used a transposon insertion site (TIS) sequencing-based strategy to identify new factors required for expression of tcpA, which encodes the major subunit of TCP, the organism's chief intestinal colonization factor. Besides identifying most of the genes known to modulate tcpA expression, the screen yielded ptsI and ptsH, which encode the enzyme I (EI) and Hpr components of the V. cholerae phosphoenolpyruvate phosphotransferase system (PTS). In addition to reduced expression of TcpA, strains lacking EI, Hpr, or the associated EIIA(Glc) protein produced less cholera toxin (CT) and had a diminished capacity to colonize the infant mouse intestine. The PTS modulates virulence gene expression by regulating expression of tcpPH and aphAB, which themselves control expression of toxT, the central activator of virulence gene expression. One mechanism by which PTS promotes virulence gene expression appears to be by modulating the amounts of intracellular cyclic AMP (cAMP). Our findings reveal that the V. cholerae PTS is an additional modulator of the ToxT regulon and demonstrate the potency of loss-of-function TIS sequencing screens for defining regulatory networks.

Characterization of DNA substrate specificities of apurinic/apyrimidinic endonucleases from Mycobacterium tuberculosis

Apurinic/apyrimidinic (AP) endonucleases are key enzymes involved in the repair of abasic sites and DNA strand breaks. Pathogenic bacteria Mycobacterium tuberculosis contains two AP endonucleases: MtbXthA and MtbNfo members of the exonuclease III and endonuclease IV families, which are exemplified by Escherichia coli Xth and Nfo, respectively. It has been shown that both MtbXthA and MtbNfo contain AP endonuclease and 3'→5' exonuclease activities. However, it remains unclear whether these enzymes hold 3'-repair phosphodiesterase and nucleotide incision repair (NIR) activities. Here, we report that both mycobacterial enzymes have 3'-repair phosphodiesterase and 3'-phosphatase, and MtbNfo contains in addition a very weak NIR activity. Interestingly, depending on pH, both enzymes require different concentrations of divalent cations: 0.5mM MnCl2 at pH 7.6 and 10 mM at pH 6.5. MtbXthA requires a low ionic strength and 37 °C, while MtbNfo requires high ionic strength (200 mM KCl) and has a temperature optimum at 60 °C. Point mutation analysis showed that D180 and N182 in MtbXthA and H206 and E129 in MtbNfo are critical for enzymes activities. The steady-state kinetic parameters indicate that MtbXthA removes 3'-blocking sugar-phosphate and 3'-phosphate moieties at DNA strand breaks with an extremely high efficiency (kcat/KM=440 and 1280 μM(-1)∙min(-1), respectively), while MtbNfo exhibits much lower 3'-repair activities (kcat/KM=0.26 and 0.65 μM(-1)∙min(-1), respectively). Surprisingly, both MtbXthA and MtbNfo exhibited very weak AP site cleavage activities, with kinetic parameters 100- and 300-fold lower, respectively, as compared with the results reported previously. Expression of MtbXthA and MtbNfo reduced the sensitivity of AP endonuclease-deficient E. coli xth nfo strain to methylmethanesulfonate and H2O2 to various degrees. Taken together, these data establish the DNA substrate specificity of M. tuberculosis AP endonucleases and suggest their possible role in the repair of oxidative DNA damage generated by endogenous and host- imposed factors.

Essential Role of the ESX-5 Secretion System in Outer Membrane Permeability of Pathogenic Mycobacteria

Mycobacteria possess different type VII secretion (T7S) systems to secrete proteins across their unusual cell envelope. One of these systems, ESX-5, is only present in slow-growing mycobacteria and responsible for the secretion of multiple substrates. However, the role of ESX-5 substrates in growth and/or virulence is largely unknown. In this study, we show that esx-5 is essential for growth of both Mycobacterium marinum and Mycobacterium bovis. Remarkably, this essentiality can be rescued by increasing the permeability of the outer membrane, either by altering its lipid composition or by the introduction of the heterologous porin MspA. Mutagenesis of the first nucleotide-binding domain of the membrane ATPase EccC5 prevented both ESX-5-dependent secretion and bacterial growth, but did not affect ESX-5 complex assembly. This suggests that the rescuing effect is not due to pores formed by the ESX-5 membrane complex, but caused by ESX-5 activity. Subsequent proteomic analysis to identify crucial ESX-5 substrates confirmed that all detectable PE and PPE proteins in the cell surface and cell envelope fractions were routed through ESX-5. Additionally, saturated transposon-directed insertion-site sequencing (TraDIS) was applied to both wild-type M. marinum cells and cells expressing mspA to identify genes that are not essential anymore in the presence of MspA. This analysis confirmed the importance of esx-5, but we could not identify essential ESX-5 substrates, indicating that multiple of these substrates are together responsible for the essentiality. Finally, examination of phenotypes on defined carbon sources revealed that an esx-5 mutant is strongly impaired in the uptake and utilization of hydrophobic carbon sources. Based on these data, we propose a model in which the ESX-5 system is responsible for the transport of cell envelope proteins that are required for nutrient uptake. These proteins might in this way compensate for the lack of MspA-like porins in slow-growing mycobacteria.

A new platform for ultra-high density Staphylococcus aureus transposon libraries

BACKGROUND

Staphylococcus aureus readily develops resistance to antibiotics and achieving effective therapies to overcome resistance requires in-depth understanding of S. aureus biology. High throughput, parallel-sequencing methods for analyzing transposon mutant libraries have the potential to revolutionize studies of S. aureus, but the genetic tools to take advantage of the power of next generation sequencing have not been fully developed.

RESULTS

Here we report a phage-based transposition system to make ultra-high density transposon libraries for genome-wide analysis of mutant fitness in any Φ11-transducible S. aureus strain. The high efficiency of the delivery system has made it possible to multiplex transposon cassettes containing different regulatory elements in order to make libraries in which genes are over- or under-expressed as well as deleted. By incorporating transposon-specific barcodes into the cassettes, we can evaluate how null mutations and changes in gene expression levels affect fitness in a single sequencing data set. Demonstrating the power of the system, we have prepared a library containing more than 690,000 unique insertions. Because one unique feature of the phage-based approach is that temperature-sensitive mutants are retained, we have carried out a genome-wide study of S. aureus genes involved in withstanding temperature stress. We find that many genes previously identified as essential are temperature sensitive and also identify a number of genes that, when disrupted, confer a growth advantage at elevated temperatures.

CONCLUSIONS

The platform described here reliably provides mutant collections of unparalleled genotypic diversity and will enable a wide range of functional genomic studies in S. aureus.

Sequence tag-based analysis of microbial population dynamics

We describe sequence tag-based analysis of microbial populations (STAMP) for characterization of pathogen population dynamics during infection. STAMP analyzes the frequency changes of genetically 'barcoded' organisms to quantify population bottlenecks and infer the founding population size. Analyses of intraintestinal Vibrio cholerae revealed infection-stage and region-specific host barriers to infection and showed unexpected V. cholerae migration counter to intestinal flow. STAMP provides a robust, widely applicable analytical framework for high-confidence characterization of in vivo microbial dissemination.

High-resolution genetic analysis of the requirements for horizontal transmission of the ESBL plasmid from Escherichia coli O104:H4

Horizontal dissemination of the genes encoding extended spectrum beta-lactamases (ESBLs) via conjugative plasmids is facilitating the increasingly widespread resistance of pathogens to beta-lactam antibiotics. However, there is relatively little known about the regulatory factors and mechanisms that govern the spread of these plasmids. Here, we carried out a high-throughput, transposon insertion site sequencing analysis (TnSeq) to identify genes that enable the maintenance and transmission of pESBL, an R64 (IncI1)-related resistance plasmid that was isolated from Escherichia coli O104:H4 linked to a recent large outbreak of gastroenteritis. With a few exceptions, the majority of the genes identified as required for maintenance and transmission of pESBL matched those of their previously defined R64 counterparts. However, our analyses of the high-density transposon insertion library in pESBL also revealed two very short and linked regions that constitute a previously unrecognized regulatory system controlling spread of IncI1 plasmids. In addition, we investigated the function of the pESBL-encoded M.EcoGIX methyltransferase, which is also encoded by many other IncI1 and IncF plasmids. This enzyme proved to protect pESBL from restriction in new hosts, suggesting it aids in expanding the plasmid's host range. Collectively, our work illustrates the power of the TnSeq approach to enable rapid and comprehensive analyses of plasmid genes and sequences that facilitate the dissemination of determinants of antibiotic resistance.