The TraDIS toolkit: sequencing and analysis for dense transposon mutant libraries

Clinical and laboratory-induced colistin-resistance mechanisms in Acinetobacter baumannii

The increasing incidence and emergence of multi-drug resistant (MDR) Acinetobacter baumannii has become a major global health concern. Colistin is a historic antimicrobial that has become commonly used as a treatment for MDR A. baumannii infections. The increase in colistin usage has been mirrored by an increase in colistin resistance. We aimed to identify the mechanisms associated with colistin resistance in A. baumannii using multiple high-throughput-sequencing technologies, including transposon-directed insertion site sequencing (TraDIS), RNA sequencing (RNAseq) and whole-genome sequencing (WGS) to investigate the genotypic changes of colistin resistance in A. baumannii. Using TraDIS, we found that genes involved in drug efflux (adeIJK), and phospholipid (mlaC, mlaF and mlaD) and lipooligosaccharide synthesis (lpxC and lpsO) were required for survival in sub-inhibitory concentrations of colistin. Transcriptomic (RNAseq) analysis revealed that expression of genes encoding efflux proteins (adeI, adeC, emrB, mexB and macAB) was enhanced in in vitro generated colistin-resistant strains. WGS of these organisms identified disruptions in genes involved in lipid A (lpxC) and phospholipid synthesis (mlaA), and in the baeS/R two-component system (TCS). We additionally found that mutations in the pmrB TCS genes were the primary colistin-resistance-associated mechanisms in three Vietnamese clinical colistin-resistant A. baumannii strains. Our results outline the entire range of mechanisms employed in A. baumannii for resistance against colistin, including drug extrusion and the loss of lipid A moieties by gene disruption or modification.

Gene fitness landscape of group A streptococcus during necrotizing myositis

Necrotizing fasciitis and myositis are devastating infections characterized by high mortality. Group A streptococcus (GAS) is a common cause of these infections, but the molecular pathogenesis is poorly understood. We report a genome-wide analysis using serotype M1 and M28 strains that identified GAS genes contributing to necrotizing myositis in nonhuman primates (NHP), a clinically relevant model. Using transposon-directed insertion-site sequencing (TraDIS), we identified 126 and 116 GAS genes required for infection by serotype M1 and M28 organisms, respectively. For both M1 and M28 strains, more than 25% of the GAS genes required for necrotizing myositis encode known or putative transporters. Thirteen GAS transporters contributed to both M1 and M28 strain fitness in NHP myositis, including putative importers for amino acids, carbohydrates, and vitamins and exporters for toxins, quorum-sensing peptides, and uncharacterized molecules. Targeted deletion of genes encoding 5 transporters confirmed that each isogenic mutant strain was significantly (P < 0.05) impaired in causing necrotizing myositis in NHPs. Quantitative reverse-transcriptase PCR (qRT-PCR) analysis showed that these 5 genes are expressed in infected NHP and human skeletal muscle. Certain substrate-binding lipoproteins of these transporters, such as Spy0271 and Spy1728, were previously documented to be surface exposed, suggesting that our findings have translational research implications.

Retrospective application of transposon-directed insertion-site sequencing to investigate niche-specific virulence of Salmonella Typhimurium in cattle

Background

Salmonella enterica subspecies enterica is an animal and zoonotic pathogen of global importance. Cattle are a significant reservoir of human non-typhoidal salmonellosis and can suffer enteric and systemic disease owing to the ability of Salmonella to survive within the bovine lymphatic system and intestines. Contamination of food can occur due to the incorporation of contaminated peripheral lymph nodes or by direct contamination of carcasses with gut contents. It is essential to understand the mechanisms used by Salmonella to enter and persist within the bovine lymphatic system and how they differ from those required for intestinal colonization to minimize zoonotic infections.

Results

Transposon-directed insertion site sequencing (TraDIS) was applied to pools of mutants recovered from mesenteric lymph nodes (MLNs) draining the distal ileum of calves after oral inoculation with a library of 8550 random S. Typhimurium mini-Tn5Km2 mutants in pools of 475 mutants per calf. A total of 8315 mutants representing 2852 different genes were detected in MLNs and their in vivo fitness was calculated. Using the same improved algorithm for analysis of transposon-flanking sequences, the identity and phenotype of mutants recovered from the distal ileal mucosa of the same calves was also defined, enabling comparison with previously published data and of mutant phenotypes across the tissues. Phenotypes observed for the majority of mutants were highly significantly correlated in the two tissues. However, 32 genes were identified in which transposon insertions consistently resulted in differential fitness in the ileal wall and MLNs, suggesting niche-specific roles for these genes in pathogenesis. Defined null mutations affecting ptsN and spvC were confirmed to result in tissue-specific phenotypes in calves, thus validating the TraDIS dataset.

Conclusions

This validation of the role of thousands of Salmonella genes and identification of genes with niche-specific roles in a key target species will inform the design of control strategies for bovine salmonellosis and zoonotic infections, for which efficacious and cross-protective vaccines are currently lacking.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5319-0) contains supplementary material, which is available to authorized users.

Adding function to the genome of African Salmonella Typhimurium ST313 strain D23580

Salmonella Typhimurium sequence type (ST) 313 causes invasive nontyphoidal Salmonella (iNTS) disease in sub-Saharan Africa, targeting susceptible HIV+, malarial, or malnourished individuals. An in-depth genomic comparison between the ST313 isolate D23580 and the well-characterized ST19 isolate 4/74 that causes gastroenteritis across the globe revealed extensive synteny. To understand how the 856 nucleotide variations generated phenotypic differences, we devised a large-scale experimental approach that involved the global gene expression analysis of strains D23580 and 4/74 grown in 16 infection-relevant growth conditions. Comparison of transcriptional patterns identified virulence and metabolic genes that were differentially expressed between D23580 versus 4/74, many of which were validated by proteomics. We also uncovered the S. Typhimurium D23580 and 4/74 genes that showed expression differences during infection of murine macrophages. Our comparative transcriptomic data are presented in a new enhanced version of the Salmonella expression compendium, SalComD23580: http://bioinf.gen.tcd.ie/cgi-bin/salcom_v2.pl. We discovered that the ablation of melibiose utilization was caused by three independent SNP mutations in D23580 that are shared across ST313 lineage 2, suggesting that the ability to catabolize this carbon source has been negatively selected during ST313 evolution. The data revealed a novel, to our knowledge, plasmid maintenance system involving a plasmid-encoded CysS cysteinyl-tRNA synthetase, highlighting the power of large-scale comparative multicondition analyses to pinpoint key phenotypic differences between bacterial pathovariants.

The Capsule Regulatory Network of Klebsiella pneumoniae Defined by density-TraDISort

Klebsiella pneumoniae infections affect infants and the immunocompromised, and the recent emergence of hypervirulent and multidrug-resistant K. pneumoniae lineages is a critical health care concern. Hypervirulence in K. pneumoniae is mediated by several factors, including the overproduction of extracellular capsule. However, the full details of how K. pneumoniae capsule biosynthesis is achieved or regulated are not known. We have developed a robust and sensitive procedure to identify genes influencing capsule production, density-TraDISort, which combines density gradient centrifugation with transposon insertion sequencing. We have used this method to explore capsule regulation in two clinically relevant Klebsiella strains, K. pneumoniae NTUH-K2044 (capsule type K1) and K. pneumoniae ATCC 43816 (capsule type K2). We identified multiple genes required for full capsule production in K. pneumoniae, as well as putative suppressors of capsule in NTUH-K2044, and have validated the results of our screen with targeted knockout mutants. Further investigation of several of the K. pneumoniae capsule regulators identified-ArgR, MprA/KvrB, SlyA/KvrA, and the Sap ABC transporter-revealed effects on capsule amount and architecture, serum resistance, and virulence. We show that capsule production in K. pneumoniae is at the center of a complex regulatory network involving multiple global regulators and environmental cues and that the majority of capsule regulatory genes are located in the core genome. Overall, our findings expand our understanding of how capsule is regulated in this medically important pathogen and provide a technology that can be easily implemented to study capsule regulation in other bacterial species.IMPORTANCE Capsule production is essential for K. pneumoniae to cause infections, but its regulation and mechanism of synthesis are not fully understood in this organism. We have developed and applied a new method for genome-wide identification of capsule regulators. Using this method, many genes that positively or negatively affect capsule production in K. pneumoniae were identified, and we use these data to propose an integrated model for capsule regulation in this species. Several of the genes and biological processes identified have not previously been linked to capsule synthesis. We also show that the methods presented here can be applied to other species of capsulated bacteria, providing the opportunity to explore and compare capsule regulatory networks in other bacterial strains and species.

Re-programming of Pseudomonas syringae pv. actinidiae gene expression during early stages of infection of kiwifruit

BACKGROUND

Pseudomonas syringae is a widespread bacterial species complex that includes a number of significant plant pathogens. Amongst these, P. syringae pv. actinidiae (Psa) initiated a worldwide pandemic in 2008 on cultivars of Actinidia chinensis var. chinensis. To gain information about the expression of genes involved in pathogenicity we have carried out transcriptome analysis of Psa during the early stages of kiwifruit infection.

RESULTS

Gene expression in Psa was investigated during the first five days after infection of kiwifruit plantlets, using RNA-seq. Principal component and heatmap analyses showed distinct phases of gene expression during the time course of infection. The first phase was an immediate transient peak of induction around three hours post inoculation (HPI) that included genes that code for a Type VI Secretion System and nutrient acquisition (particularly phosphate). This was followed by a significant commitment, between 3 and 24 HPI, to the induction of genes encoding the Type III Secretion System (T3SS) and Type III Secreted Effectors (T3SE). Expression of these genes collectively accounted for 6.3% of the bacterial transcriptome at this stage. There was considerable variation in the expression levels of individual T3SEs but all followed the same temporal expression pattern, with the exception of hopAS1, which peaked later in expression at 48 HPI. As infection progressed over the time course of five days, there was an increase in the expression of genes with roles in sugar, amino acid and sulfur transport and the production of alginate and colanic acid. These are both polymers that are major constituents of extracellular polysaccharide substances (EPS) and are involved in biofilm production. Reverse transcription-quantitative PCR (RT-qPCR) on an independent infection time course experiment showed that the expression profile of selected bacterial genes at each infection phase correlated well with the RNA-seq data.

CONCLUSIONS

The results from this study indicate that there is a complex remodeling of the transcriptome during the early stages of infection, with at least three distinct phases of coordinated gene expression. These include genes induced during the immediate contact with the host, those involved in the initiation of infection, and finally those responsible for nutrient acquisition.

Genome-wide CRISPR-dCas9 screens in E. coli identify essential genes and phage host factors

High-throughput genetic screens are powerful methods to identify genes linked to a given phenotype. The catalytic null mutant of the Cas9 RNA-guided nuclease (dCas9) can be conveniently used to silence genes of interest in a method also known as CRISPRi. Here, we report a genome-wide CRISPR-dCas9 screen using a starting pool of ~ 92,000 sgRNAs which target random positions in the chromosome of E. coli. To benchmark our method, we first investigate its utility to predict gene essentiality in the genome of E. coli during growth in rich medium. We could identify 79% of the genes previously reported as essential and demonstrate the non-essentiality of some genes annotated as essential. In addition, we took advantage of the intermediate repression levels obtained when targeting the template strand of genes to show that cells are very sensitive to the expression level of a limited set of essential genes. Our data can be visualized on CRISPRbrowser, a custom web interface available at crispr.pasteur.fr. We then apply the screen to discover E. coli genes required by phages λ, T4 and 186 to kill their host, highlighting the involvement of diverse host pathways in the infection process of the three tested phages. We also identify colanic acid capsule synthesis as a shared resistance mechanism to all three phages. Finally, using a plasmid packaging system and a transduction assay, we identify genes required for the formation of functional λ capsids, thus covering the entire phage cycle. This study demonstrates the usefulness and convenience of pooled genome-wide CRISPR-dCas9 screens in bacteria and paves the way for their broader use as a powerful tool in bacterial genomics.

Over the past few years, CRISPR-Cas technologies have emerged as powerful tools to edit genomes and modulate gene expression. They have been applied to perform high-throughput genetic screens with the purpose to understand the function of genes in a systematic manner, but the application of these screens to bacteria have so far remained limited. Here, we present the use of a library of ~92,000 guide RNAs directing the dCas9 protein to silence one by one all the genes in the chromosome of E. coli. To benchmark our method, we first investigate the performance of the technique to identify essential genes, highlighting several non-essential genes also found to be essential by other methods. We then apply our method to detect bacterial genes required by three different bacteriophages to kill E. coli and for the production of functional progeny by phage λ. Our screens highlight previously known and new genetic interactions between phages and their host’s pathways and emphasize the importance of bacterial capsule in the resistance to multiple phages. Altogether, our results demonstrate the usefulness of genome-wide CRISPR-dCas9 screens in bacteria to uncover genes involved in various phenotypes.

A global genomic approach uncovers novel components for twitching motility-mediated biofilm expansion in Pseudomonas aeruginosa

Pseudomonas aeruginosa is an extremely successful pathogen able to cause both acute and chronic infections in a range of hosts, utilizing a diverse arsenal of cell-associated and secreted virulence factors. A major cell-associated virulence factor, the Type IV pilus (T4P), is required for epithelial cell adherence and mediates a form of surface translocation termed twitching motility, which is necessary to establish a mature biofilm and actively expand these biofilms. P. aeruginosa twitching motility-mediated biofilm expansion is a coordinated, multicellular behaviour, allowing cells to rapidly colonize surfaces, including implanted medical devices. Although at least 44 proteins are known to be involved in the biogenesis, assembly and regulation of the T4P, with additional regulatory components and pathways implicated, it is unclear how these components and pathways interact to control these processes. In the current study, we used a global genomics-based random-mutagenesis technique, transposon directed insertion-site sequencing (TraDIS), coupled with a physical segregation approach, to identify all genes implicated in twitching motility-mediated biofilm expansion in P. aeruginosa. Our approach allowed identification of both known and novel genes, providing new insight into the complex molecular network that regulates this process in P. aeruginosa. Additionally, our data suggest that the flagellum-associated gene products have a differential effect on twitching motility, based on whether components are intra- or extracellular. Overall the success of our TraDIS approach supports the use of this global genomic technique for investigating virulence genes in bacterial pathogens.

Rapid, Parallel Identification of Catabolism Pathways of Lignin-Derived Aromatic Compounds in Novosphingobium aromaticivorans

Transposon mutagenesis is a powerful technique in microbial genetics for the identification of genes in uncharacterized pathways. Recently, the throughput of transposon mutagenesis techniques has been dramatically increased through the combination of DNA barcoding and high-throughput sequencing. Here, we show that when applied to catabolic pathways, barcoded transposon libraries can be used to distinguish redundant pathways, decompose complex pathways into substituent modules, discriminate between enzyme homologs, and rapidly identify previously hypothetical enzymes in an unbiased genome-scale search. We used this technique to identify two genes, desC and desD, which are involved in the degradation of the lignin-derived aromatic compound sinapic acid in the nonmodel bacterium Novosphingobium aromaticivorans We show that DesC is a methyl esterase acting on an intermediate formed during sinapic acid catabolism, providing the last enzyme in a proposed catabolic pathway. This approach will be particularly useful in the identification of complete pathways suitable for heterologous expression in metabolic engineering.IMPORTANCE The identification of the genes involved in specific biochemical transformations is a key step in predicting microbial function from nucleic acid sequences and in engineering microbes to endow them with new functions. We have shown that new techniques for transposon mutagenesis can dramatically simplify this process and enable the rapid identification of genes in uncharacterized pathways. These techniques provide the necessary scale to fully elucidate complex biological networks such as those used to degrade mixtures of lignin-derived aromatic compounds.

Transposon Insertion Sequencing Elucidates Novel Gene Involvement in Susceptibility and Resistance to Phages T4 and T7 in Escherichia coli O157

Experiments using bacteriophage (phage) to infect bacterial strains have helped define some basic genetic concepts in microbiology, but our understanding of the complexity of bacterium-phage interactions is still limited. As the global threat of antibiotic resistance continues to increase, phage therapy has reemerged as an attractive alternative or supplement to treating antibiotic-resistant bacterial infections. Further, the long-used method of phage typing to classify bacterial strains is being replaced by molecular genetic techniques. Thus, there is a growing need for a complete understanding of the precise molecular mechanisms underpinning phage-bacterium interactions to optimize phage therapy for the clinic as well as for retrospectively interpreting phage typing data on the molecular level. In this study, a genomics-based fitness assay (TraDIS) was used to identify all host genes involved in phage susceptibility and resistance for a T4 phage infecting Shiga-toxigenic Escherichia coli O157. The TraDIS results identified both established and previously unidentified genes involved in phage infection, and a subset were confirmed by site-directed mutagenesis and phenotypic testing of 14 T4 and 2 T7 phages. For the first time, the entire sap operon was implicated in phage susceptibility and, conversely, the stringent starvation protein A gene (sspA) was shown to provide phage resistance. Identifying genes involved in phage infection and replication should facilitate the selection of bespoke phage combinations to target specific bacterial pathogens.IMPORTANCE Antibiotic resistance has diminished treatment options for many common bacterial infections. Phage therapy is an alternative option that was once popularly used across Europe to kill bacteria within humans. Phage therapy acts by using highly specific viruses (called phages) that infect and lyse certain bacterial species to treat the infection. Whole-genome sequencing has allowed modernization of the investigations into phage-bacterium interactions. Here, using E. coli O157 and T4 bacteriophage as a model, we have exploited a genome-wide fitness assay to investigate all genes involved in defining phage resistance or susceptibility. This knowledge of the genetic determinants of phage resistance and susceptibility can be used to design bespoke phage combinations targeted to specific bacterial infections for successful infection eradication.

Transposon Insertion Sequencing Elucidates Novel Gene Involvement in Susceptibility and Resistance to Phages T4 and T7 in Escherichia coli O157

Experiments using bacteriophage (phage) to infect bacterial strains have helped define some basic genetic concepts in microbiology, but our understanding of the complexity of bacterium-phage interactions is still limited. As the global threat of antibiotic resistance continues to increase, phage therapy has reemerged as an attractive alternative or supplement to treating antibiotic-resistant bacterial infections. Further, the long-used method of phage typing to classify bacterial strains is being replaced by molecular genetic techniques. Thus, there is a growing need for a complete understanding of the precise molecular mechanisms underpinning phage-bacterium interactions to optimize phage therapy for the clinic as well as for retrospectively interpreting phage typing data on the molecular level. In this study, a genomics-based fitness assay (TraDIS) was used to identify all host genes involved in phage susceptibility and resistance for a T4 phage infecting Shiga-toxigenic Escherichia coli O157. The TraDIS results identified both established and previously unidentified genes involved in phage infection, and a subset were confirmed by site-directed mutagenesis and phenotypic testing of 14 T4 and 2 T7 phages. For the first time, the entire sap operon was implicated in phage susceptibility and, conversely, the stringent starvation protein A gene (sspA) was shown to provide phage resistance. Identifying genes involved in phage infection and replication should facilitate the selection of bespoke phage combinations to target specific bacterial pathogens.IMPORTANCE Antibiotic resistance has diminished treatment options for many common bacterial infections. Phage therapy is an alternative option that was once popularly used across Europe to kill bacteria within humans. Phage therapy acts by using highly specific viruses (called phages) that infect and lyse certain bacterial species to treat the infection. Whole-genome sequencing has allowed modernization of the investigations into phage-bacterium interactions. Here, using E. coli O157 and T4 bacteriophage as a model, we have exploited a genome-wide fitness assay to investigate all genes involved in defining phage resistance or susceptibility. This knowledge of the genetic determinants of phage resistance and susceptibility can be used to design bespoke phage combinations targeted to specific bacterial infections for successful infection eradication.

High-throughput analysis of Yersinia pseudotuberculosis gene essentiality in optimised in vitro conditions, and implications for the speciation of Yersinia pestis

Background

Yersinia pseudotuberculosis is a zoonotic pathogen, causing mild gastrointestinal infection in humans. From this comparatively benign pathogenic species emerged the highly virulent plague bacillus, Yersinia pestis, which has experienced significant genetic divergence in a relatively short time span. Much of our knowledge of Yersinia spp. evolution stems from genomic comparison and gene expression studies. Here we apply transposon-directed insertion site sequencing (TraDIS) to describe the essential gene set of Y. pseudotuberculosis IP32953 in optimised in vitro growth conditions, and contrast these with the published essential genes of Y. pestis.

Results

The essential genes of an organism are the core genetic elements required for basic survival processes in a given growth condition, and are therefore attractive targets for antimicrobials. One such gene we identified is yptb3665, which encodes a peptide deformylase, and here we report for the first time, the sensitivity of Y. pseudotuberculosis to actinonin, a deformylase inhibitor. Comparison of the essential genes of Y. pseudotuberculosis with those of Y. pestis revealed the genes whose importance are shared by both species, as well as genes that were differentially required for growth. In particular, we find that the two species uniquely rely upon different iron acquisition and respiratory metabolic pathways under similar in vitro conditions.

Conclusions

The discovery of uniquely essential genes between the closely related Yersinia spp. represent some of the fundamental, species-defining points of divergence that arose during the evolution of Y. pestis from its ancestor. Furthermore, the shared essential genes represent ideal candidates for the development of novel antimicrobials against both species.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1189-5) contains supplementary material, which is available to authorized users.

Use of Proteins Identified through a Functional Genomic Screen To Develop a Protein Subunit Vaccine That Provides Significant Protection against Virulent Streptococcus suis in Pigs

Streptococcus suis is a bacterium that is commonly carried in the respiratory tract and that is also one of the most important invasive pathogens of swine, commonly causing meningitis, arthritis, and septicemia. Due to the existence of many serotypes and a wide range of immune evasion capabilities, efficacious vaccines are not readily available. The selection of S. suis protein candidates for inclusion in a vaccine was accomplished by identifying fitness genes through a functional genomics screen and selecting conserved predicted surface-associated proteins. Five candidate proteins were selected for evaluation in a vaccine trial and administered both intranasally and intramuscularly with one of two different adjuvant formulations. Clinical protection was evaluated by subsequent intranasal challenge with virulent S. suis While subunit vaccination with the S. suis proteins induced IgG antibodies to each individual protein and a cellular immune response to the pool of proteins and provided substantial protection from challenge with virulent S. suis, the immune response elicited and the degree of protection were dependent on the parenteral adjuvant given. Subunit vaccination induced IgG reactive against different S. suis serotypes, indicating a potential for cross protection.

Molecular phenotyping of infection-associated small non-coding RNAs

Infection is a complicated balance, with both pathogen and host struggling to tilt the result in their favour. Bacterial infection biology has relied on forward genetics for many of its advances, defining phenotype in terms of replication in model systems. However, many known virulence factors fail to produce robust phenotypes, particularly in the systems most amenable to genetic manipulation, such as cell-culture models. This has particularly been limiting for the study of the bacterial regulatory small RNAs (sRNAs) in infection. We argue that new sequencing-based technologies can work around this problem by providing a 'molecular phenotype', defined in terms of the specific transcriptional dysregulation in the infection system induced by gene deletion. We illustrate this using the example of our recent study of the PinT sRNA using dual RNA-seq, that is, simultaneous RNA sequencing of host and pathogen during infection. We additionally discuss how other high-throughput technologies, in particular genetic interaction mapping using transposon insertion sequencing, may be used to further dissect molecular phenotypes. We propose a strategy for how high-throughput technologies can be integrated in the study of non-coding regulators as well as bacterial virulence factors, enhancing our ability to rapidly generate hypotheses with regards to their function.This article is part of the themed issue 'The new bacteriology'.

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.

Novel Genes Required for the Fitness of Streptococcus pyogenes in Human Saliva

Streptococcus pyogenes (group A streptococcus [GAS]) causes 600 million cases of pharyngitis each year. Despite this considerable disease burden, the molecular mechanisms used by GAS to infect, cause clinical pharyngitis, and persist in the human oropharynx are poorly understood. Saliva is ubiquitous in the human oropharynx and is the first material GAS encounters in the upper respiratory tract. Thus, a fuller understanding of how GAS survives and proliferates in saliva may provide valuable insights into the molecular mechanisms at work in the human oropharynx. We generated a highly saturated transposon insertion mutant library in serotype M1 strain MGAS2221, a strain genetically representative of a pandemic clone that arose in the 1980s and spread globally. The transposon mutant library was exposed to human saliva to screen for GAS genes required for wild-type fitness in this clinically relevant fluid. Using transposon-directed insertion site sequencing (TraDIS), we identified 92 genes required for GAS fitness in saliva. The more prevalent categories represented were genes involved in carbohydrate transport/metabolism, amino acid transport/metabolism, and inorganic ion transport/metabolism. Using six isogenic mutant strains, we confirmed that each of the mutants was significantly impaired for growth or persistence in human saliva ex vivo. Mutants with an inactivated Spy0644 (sptA) or Spy0646 (sptC) gene had especially severe persistence defects. This study is the first to use of TraDIS to study bacterial fitness in human saliva. The new information we obtained will be valuable for future translational maneuvers designed to prevent or treat human GAS infections. IMPORTANCE The human bacterial pathogen Streptococcus pyogenes (group A streptococcus [GAS]) causes more than 600 million cases of pharyngitis annually worldwide, 15 million of which occur in the United States. The human oropharynx is the primary anatomic site for GAS colonization and infection, and saliva is the first material encountered. Using a genome-wide transposon mutant screen, we identified 92 GAS genes required for wild-type fitness in human saliva. Many of the identified genes are involved in carbohydrate transport/metabolism, amino acid transport/metabolism, and inorganic ion transport/metabolism. The new information is potentially valuable for developing novel GAS therapeutics and vaccine research.

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.

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.

A curated genome-scale metabolic model of Bordetella pertussis metabolism

The Gram-negative bacterium Bordetella pertussis is the causative agent of whooping cough, a serious respiratory infection causing hundreds of thousands of deaths annually worldwide. There are effective vaccines, but their production requires growing large quantities of B. pertussis. Unfortunately, B. pertussis has relatively slow growth in culture, with low biomass yields and variable growth characteristics. B. pertussis also requires a relatively expensive growth medium. We present a new, curated flux balance analysis-based model of B. pertussis metabolism. We enhance the model with an experimentally-determined biomass objective function, and we perform extensive manual curation. We test the model’s predictions with a genome-wide screen for essential genes using a transposon-directed insertional sequencing (TraDIS) approach. We test its predictions of growth for different carbon sources in the medium. The model predicts essentiality with an accuracy of 83% and correctly predicts improvements in growth under increased glutamate:fumarate ratios. We provide the model in SBML format, along with gene essentiality predictions.

Metabolic flux models have been used to understand how organisms adapt their metabolism under different growth conditions, and are finding increasing application in synthetic biology and biotechnology. One barrier to progress in this field is the construction and curation of metabolic flux models for new organisms. Here we present a curated genome-scale metabolic flux model for Bordetella pertussis, the causative agent of whooping cough. Producing vaccines against whooping cough requires growing B. pertussis in large volumes. However, its growth is relatively slow, final yields of biomass are relatively low and growth characteristics can be variable. Understanding B. pertussis metabolism has applications to improving vaccine production, as well as in understanding the basic biology of this organism.

Defining the ABC of gene essentiality in streptococci

BACKGROUND

Utilising next generation sequencing to interrogate saturated bacterial mutant libraries provides unprecedented information for the assignment of genome-wide gene essentiality. Exposure of saturated mutant libraries to specific conditions and subsequent sequencing can be exploited to uncover gene essentiality relevant to the condition. Here we present a barcoded transposon directed insertion-site sequencing (TraDIS) system to define an essential gene list for Streptococcus equi subsp. equi, the causative agent of strangles in horses, for the first time. The gene essentiality data for this group C Streptococcus was compared to that of group A and B streptococci.

RESULTS

Six barcoded variants of pGh9:ISS1 were designed and used to generate mutant libraries containing between 33,000-66,000 unique mutants. TraDIS was performed on DNA extracted from each library and data were analysed separately and as a combined master pool. Gene essentiality determined that 19.5% of the S. equi genome was essential. Gene essentialities were compared to those of group A and group B streptococci, identifying concordances of 90.2% and 89.4%, respectively and an overall concordance of 83.7% between the three species.

CONCLUSIONS

The use of barcoded pGh9:ISS1 to generate mutant libraries provides a highly useful tool for the assignment of gene function in S. equi and other streptococci. The shared essential gene set of group A, B and C streptococci provides further evidence of the close genetic relationships between these important pathogenic bacteria. Therefore, the ABC of gene essentiality reported here provides a solid foundation towards reporting the functional genome of streptococci.

Physical enrichment of transposon mutants from saturation mutant libraries using the TraDISort approach

Transposon-insertion sequencing methods are finding their way into the molecular toolbox of many fields of microbiology. These methods can identify the genomic locations and density of transposon insertions in saturated transposon mutant libraries and can be used to make inferences on gene function. For example, where no insertions or very few insertions are identified within a gene in a mutant library grown under permissive conditions, the gene may be essential. Furthermore, where mutations are enriched or lost in a gene after passaging the library through a selective process, the gene is likely to be involved in phenotypes linked to the process. Typically, a fitness based selection such as a stress condition is used in these experiments and the processed sequencing data are used to identify genes required for fitness under the selection. Our research team recently expanded the utility of the transposon directed insertion sequencing (TraDIS) method by applying a physical separation of a transposon mutant library mediated by fluorescence activated cell sorting, rather than a fitness-based selection. This approach, which we have named "TraDISort" is significant because it allows the study of phenotypes that are not linked to cell survival. The TraDISort approach has a broad range of future applications, in drug development, metabolic engineering and in studies of basic bacterial cell physiology.

Transposase-Mediated Excision, Conjugative Transfer, and Diversity of ICE6013 Elements in Staphylococcus aureus

ICE6013 represents one of two families of integrative conjugative elements (ICEs) identified in the pan-genome of the human and animal pathogen Staphylococcus aureus Here we investigated the excision and conjugation functions of ICE6013 and further characterized the diversity of this element. ICE6013 excision was not significantly affected by growth, temperature, pH, or UV exposure and did not depend on recA The IS30-like DDE transposase (Tpase; encoded by orf1 and orf2) of ICE6013 must be uninterrupted for excision to occur, whereas disrupting three of the other open reading frames (ORFs) on the element significantly affects the level of excision. We demonstrate that ICE6013 conjugatively transfers to different S. aureus backgrounds at frequencies approaching that of the conjugative plasmid pGO1. We found that excision is required for conjugation, that not all S. aureus backgrounds are successful recipients, and that transconjugants acquire the ability to transfer ICE6013 Sequencing of chromosomal integration sites in serially passaged transconjugants revealed a significant integration site preference for a 15-bp AT-rich palindromic consensus sequence, which surrounds the 3-bp target site that is duplicated upon integration. A sequence analysis of ICE6013 from different host strains of S. aureus and from eight other species of staphylococci identified seven divergent subfamilies of ICE6013 that include sequences previously classified as a transposon, a plasmid, and various ICEs. In summary, these results indicate that the IS30-like Tpase functions as the ICE6013 recombinase and that ICE6013 represents a diverse family of mobile genetic elements that mediate conjugation in staphylococci.IMPORTANCE Integrative conjugative elements (ICEs) encode the abilities to integrate into and excise from bacterial chromosomes and plasmids and mediate conjugation between bacteria. As agents of horizontal gene transfer, ICEs may affect bacterial evolution. ICE6013 represents one of two known families of ICEs in the pathogen Staphylococcus aureus, but its core functions of excision and conjugation are not well studied. Here, we show that ICE6013 depends on its IS30-like DDE transposase for excision, which is unique among ICEs, and we demonstrate the conjugative transfer and integration site preference of ICE6013 A sequence analysis revealed that ICE6013 has diverged into seven subfamilies that are dispersed among staphylococci.

Vector Integration Sites Identification for Gene-Trap Screening in Mammalian Haploid Cells

Forward genetic screens using retroviral (or transposon) gene-trap vectors in a haploid genome revolutionized the investigation of molecular networks in mammals. However, the sequencing data generated by Phenotypic interrogation followed by Tag sequencing (PhiT-seq) were not well characterized. The analysis of human and mouse haploid screens allowed us to describe PhiT-seq data and to define quality control steps. Moreover, we identified several blind spots in both haploid genomes where gene-trap vectors can hardly integrate. Integration of transcriptomic data improved the performance of candidate gene identification. Furthermore, we experimented with various statistical tests to account for biological replicates in PhiT-seq and investigated the effect of normalization methods and other parameters on the performance. Finally, we developed: VISITs, a dedicated pipeline for analyzing PhiT-seq data (https://sourceforge.net/projects/visits/).

Transposon insertion libraries for the characterization of mutants from the kiwifruit pathogen Pseudomonas syringae pv. actinidiae

Pseudomonas syringae pv. actinidiae (Psa), the causal agent of kiwifruit canker, is one of the most devastating plant diseases of recent times. We have generated two mini-Tn5-based random insertion libraries of Psa ICMP 18884. The first, a 'phenotype of interest' (POI) library, consists of 10,368 independent mutants gridded into 96-well plates. By replica plating onto selective media, the POI library was successfully screened for auxotrophic and motility mutants. Lipopolysaccharide (LPS) biosynthesis mutants with 'Fuzzy-Spreader'-like morphologies were also identified through a visual screen. The second, a 'mutant of interest' (MOI) library, comprises around 96,000 independent mutants, also stored in 96-well plates, with approximately 200 individuals per well. The MOI library was sequenced on the Illumina MiSeq platform using Transposon-Directed Insertion site Sequencing (TraDIS) to map insertion sites onto the Psa genome. A grid-based PCR method was developed to recover individual mutants, and using this strategy, the MOI library was successfully screened for a putative LPS mutant not identified in the visual screen. The Psa chromosome and plasmid had 24,031 and 1,236 independent insertion events respectively, giving insertion frequencies of 3.65 and 16.6 per kb respectively. These data suggest that the MOI library is near saturation, with the theoretical probability of finding an insert in any one chromosomal gene estimated to be 97.5%. However, only 47% of chromosomal genes had insertions. This surprisingly low rate cannot be solely explained by the lack of insertions in essential genes, which would be expected to be around 5%. Strikingly, many accessory genes, including most of those encoding type III effectors, lacked insertions. In contrast, 94% of genes on the Psa plasmid had insertions, including for example, the type III effector HopAU1. These results suggest that some chromosomal sites are rendered inaccessible to transposon insertion, either by DNA-binding proteins or by the architecture of the nucleoid.

The secondary resistome of multidrug-resistant Klebsiella pneumoniae

Klebsiella pneumoniae causes severe lung and bloodstream infections that are difficult to treat due to multidrug resistance. We hypothesized that antimicrobial resistance can be reversed by targeting chromosomal non-essential genes that are not responsible for acquired resistance but essential for resistant bacteria under therapeutic concentrations of antimicrobials. Conditional essentiality of individual genes to antimicrobial resistance was evaluated in an epidemic multidrug-resistant clone of K. pneumoniae (ST258). We constructed a high-density transposon mutant library of >430,000 unique Tn5 insertions and measured mutant depletion upon exposure to three clinically relevant antimicrobials (colistin, imipenem or ciprofloxacin) by Transposon Directed Insertion-site Sequencing (TraDIS). Using this high-throughput approach, we defined three sets of chromosomal non-essential genes essential for growth during exposure to colistin (n = 35), imipenem (n = 1) or ciprofloxacin (n = 1) in addition to known resistance determinants, collectively termed the "secondary resistome". As proof of principle, we demonstrated that inactivation of a non-essential gene not previously found linked to colistin resistance (dedA) restored colistin susceptibility by reducing the minimum inhibitory concentration from 8 to 0.5 μg/ml, 4-fold below the susceptibility breakpoint (S ≤ 2 μg/ml). This finding suggests that the secondary resistome is a potential target for developing antimicrobial "helper" drugs that restore the efficacy of existing antimicrobials.

A Noise Trimming and Positional Significance of Transposon Insertion System to Identify Essential Genes in Yersinia pestis

Massively parallel sequencing technology coupled with saturation mutagenesis has provided new and global insights into gene functions and roles. At a simplistic level, the frequency of mutations within genes can indicate the degree of essentiality. However, this approach neglects to take account of the positional significance of mutations - the function of a gene is less likely to be disrupted by a mutation close to the distal ends. Therefore, a systematic bioinformatics approach to improve the reliability of essential gene identification is desirable. We report here a parametric model which introduces a novel mutation feature together with a noise trimming approach to predict the biological significance of Tn5 mutations. We show improved performance of essential gene prediction in the bacterium Yersinia pestis, the causative agent of plague. This method would have broad applicability to other organisms and to the identification of genes which are essential for competitiveness or survival under a broad range of stresses.

Fluorescence-Based Flow Sorting in Parallel with Transposon Insertion Site Sequencing Identifies Multidrug Efflux Systems in Acinetobacter baumannii

Multidrug efflux pumps provide clinically significant levels of drug resistance in a number of Gram-negative hospital-acquired pathogens. These pathogens frequently carry dozens of genes encoding putative multidrug efflux pumps. However, it can be difficult to determine how many of these pumps actually mediate antimicrobial efflux, and it can be even more challenging to identify the regulatory proteins that control expression of these pumps. In this study, we developed an innovative high-throughput screening method, combining transposon insertion sequencing and cell sorting methods (TraDISort), to identify the genes encoding major multidrug efflux pumps, regulators, and other factors that may affect the permeation of antimicrobials, using the nosocomial pathogen Acinetobacter baumannii. A dense library of more than 100,000 unique transposon insertion mutants was treated with ethidium bromide, a common substrate of multidrug efflux pumps that is differentially fluorescent inside and outside the bacterial cytoplasm. Populations of cells displaying aberrant accumulations of ethidium were physically enriched using fluorescence-activated cell sorting, and the genomic locations of transposon insertions within these strains were determined using transposon-directed insertion sequencing. The relative abundance of mutants in the input pool compared to the selected mutant pools indicated that the AdeABC, AdeIJK, and AmvA efflux pumps are the major ethidium efflux systems in A. baumannii. Furthermore, the method identified a new transcriptional regulator that controls expression of amvA. In addition to the identification of efflux pumps and their regulators, TraDISort identified genes that are likely to control cell division, cell morphology, or aggregation in A. baumannii.

Transposon-directed insertion sequencing (TraDIS) and related technologies have emerged as powerful methods to identify genes required for bacterial survival or competitive fitness under various selective conditions. We applied fluorescence-activated cell sorting (FACS) to physically enrich for phenotypes of interest within a mutant population prior to TraDIS. To our knowledge, this is the first time that a physical selection method has been applied in parallel with TraDIS rather than a fitness-induced selection. The results demonstrate the feasibility of this combined approach to generate significant results and highlight the major multidrug efflux pumps encoded in an important pathogen. This FACS-based approach, TraDISort, could have a range of future applications, including the characterization of efflux pump inhibitors, the identification of regulatory factors controlling gene or protein expression using fluorescent reporters, and the identification of genes involved in cell replication, morphology, and aggregation.