Surviving the Potato Stems: Differences in Genes Required for Fitness by Dickeya dadantii and Dickeya dianthicola

Bacteria belonging to the genus Dickeya cause blackleg and soft rot symptoms on many plant hosts, including potato. Although there is considerable knowledge about the genetic determinants that allow Dickeya to colonize host plants, as well as the genes that contribute to virulence, much is still unknown. To identify the genes important for fitness in potato stems, we constructed and evaluated randomly barcoded transposon mutant (RB-TnSeq) libraries of Dickeya dadantii and Dickeya dianthicola. We identified 169 and 157 genes important for growth in D. dadantii and D. dianthicola in stems, respectively. This included genes related to metabolic pathways, chemotaxis and motility, transcriptional regulation, transport across membranes, membrane biogenesis, detoxification mechanisms, and virulence-related genes, including a potential virulence cluster srfABC, c-di-GMP modulating genes, and pectin degradation genes. When we compared the results of the stem assay with other datasets, we identified genes important for growth in stems versus tubers and in vitro conditions. Additionally, our data showed differences in fitness determinants for D. dadantii and D. dianthicola. These data provide important insights into the mechanisms used by Dickeya when interacting with and colonizing plants and thus might provide targets for management.

Massively parallel mutant selection identifies genetic determinants of Pseudomonas aeruginosa colonization of Drosophila melanogaster

Pseudomonas aeruginosa is recognized for its ability to colonize diverse habitats and cause disease in a variety of hosts, including plants, invertebrates, and mammals. Understanding how this bacterium is able to occupy wide-ranging niches is important for deciphering its ecology. We used transposon sequencing [Tn-Seq, also known as insertion sequencing (INSeq)] to identify genes in P. aeruginosa that contribute to fitness during the colonization of Drosophila melanogaster. Our results reveal a suite of critical factors, including those that contribute to polysaccharide production, DNA repair, metabolism, and respiration. Comparison of candidate genes with fitness determinants discovered in previous studies on P. aeruginosa identified several genes required for colonization and virulence determinants that are conserved across hosts and tissues. This analysis provides evidence for both the conservation of function of several genes across systems, as well as host-specific functions. These findings, which represent the first use of transposon sequencing of a gut pathogen in Drosophila, demonstrate the power of Tn-Seq in the fly model system and advance the existing knowledge of intestinal pathogenesis by D. melanogaster, revealing bacterial colonization determinants that contribute to a comprehensive portrait of P. aeruginosa lifestyles across habitats.IMPORTANCEDrosophila melanogaster is a powerful model for understanding host-pathogen interactions. Research with this system has yielded notable insights into mechanisms of host immunity and defense, many of which emerged from the analysis of bacterial mutants defective for well-characterized virulence factors. These foundational studies-and advances in high-throughput sequencing of transposon mutants-support unbiased screens of bacterial mutants in the fly. To investigate mechanisms of host-pathogen interplay and exploit the tractability of this model host, we used a high-throughput, genome-wide mutant analysis to find genes that enable the pathogen P. aeruginosa to colonize the fly. Our analysis reveals critical mediators of P. aeruginosa establishment in its host, some of which are required across fly and mouse systems. These findings demonstrate the utility of massively parallel mutant analysis and provide a platform for aligning the fly toolkit with comprehensive bacterial genomics.

Surface colonization by Flavobacterium johnsoniae promotes its survival in a model microbial community

Flavobacterium johnsoniae is a ubiquitous soil and rhizosphere bacterium, but despite its abundance, the factors contributing to its success in communities are poorly understood. Using a model microbial community, The Hitchhikers of the Rhizosphere (THOR), we determined the effects of colonization on the fitness of F. johnsoniae in the community. Insertion sequencing, a massively parallel transposon mutant screen, on sterile sand identified 25 genes likely to be important for surface colonization. We constructed in-frame deletions of candidate genes predicted to be involved in cell membrane biogenesis, motility, signal transduction, and transport of amino acids and lipids. All mutants poorly colonized sand, glass, and polystyrene and produced less biofilm than the wild type, indicating the importance of the targeted genes in surface colonization. Eight of the nine colonization-defective mutants were also unable to form motile biofilms or zorbs, thereby suggesting that the affected genes play a role in group movement and linking stationary and motile biofilm formation genetically. Furthermore, we showed that the deletion of colonization genes in F. johnsoniae affected its behavior and survival in THOR on surfaces, suggesting that the same traits are required for success in a multispecies microbial community. Our results provide insight into the mechanisms of surface colonization by F. johnsoniae and form the basis for further understanding its ecology in the rhizosphere.

IMPORTANCE

Microbial communities direct key environmental processes through multispecies interactions. Understanding these interactions is vital for manipulating microbiomes to promote health in human, environmental, and agricultural systems. However, microbiome complexity can hinder our understanding of the underlying mechanisms in microbial community interactions. As a first step toward unraveling these interactions, we explored the role of surface colonization in microbial community interactions using The Hitchhikers Of the Rhizosphere (THOR), a genetically tractable model community of three bacterial species, Flavobacterium johnsoniae, Pseudomonas koreensis, and Bacillus cereus. We identified F. johnsoniae genes important for surface colonization in solitary conditions and in the THOR community. Understanding the mechanisms that promote the success of bacteria in microbial communities brings us closer to targeted manipulations to achieve outcomes that benefit agriculture, the environment, and human health.

Restriction of Arginine Induces Antibiotic Tolerance in Staphylococcus aureus

Staphylococcus aureus is responsible for a substantial number of invasive infections globally each year. These infections are problematic because they are frequently recalcitrant to antibiotic treatment, particularly when they are caused by Methicillin-Resistant Staphylococcus aureus (MRSA). Antibiotic tolerance, the ability for bacteria to persist despite normally lethal doses of antibiotics, is responsible for most antibiotic treatment failure in MRSA infections. To understand how antibiotic tolerance is induced, S. aureus biofilms exposed to multiple anti-MRSA antibiotics (vancomycin, ceftaroline, delafloxacin, and linezolid) were examined using both quantitative proteomics and transposon sequencing. These screens indicated that arginine metabolism is involved in antibiotic tolerance within a biofilm and led to the hypothesis that depletion of arginine within S. aureus communities can induce antibiotic tolerance. Consistent with this hypothesis, inactivation of argH, the final gene in the arginine synthesis pathway, induces antibiotic tolerance under conditions in which the parental strain is susceptible to antibiotics. Arginine restriction was found to induce antibiotic tolerance via inhibition of protein synthesis. Finally, although S. aureus fitness in a mouse skin infection model is decreased in an argH mutant, its ability to survive in vivo during antibiotic treatment with vancomycin is enhanced, highlighting the relationship between arginine metabolism and antibiotic tolerance during S. aureus infection. Uncovering this link between arginine metabolism and antibiotic tolerance has the potential to open new therapeutic avenues targeting previously recalcitrant S. aureus infections.

Comprehensive essentiality analysis of the Mycobacterium kansasii genome by saturation transposon mutagenesis and deep sequencing

Mycobacterium kansasii (Mk) is an opportunistic pathogen that is frequently isolated from urban water systems, posing a health risk to susceptible individuals. Despite its ability to cause tuberculosis-like pulmonary disease, very few studies have probed the genetics of this opportunistic pathogen. Here, we report a comprehensive essentiality analysis of the Mk genome. Deep sequencing of a high-density library of Mk Himar1 transposon mutants revealed that 86.8% of the chromosomal thymine-adenine (TA) dinucleotide target sites were permissive to insertion, leaving 13.2% TA sites unoccupied. Our analysis identified 394 of the 5,350 annotated open reading frames (ORFs) as essential. The majority of these essential ORFs (84.8%) share essential mutual orthologs with Mycobacterium tuberculosis (Mtb). A comparative genomics analysis identified 139 Mk essential ORFs that share essential orthologs in four other species of mycobacteria. Thirteen Mk essential ORFs share orthologs in all four species that were identified as being not essential, while only two Mk essential ORFs are absent in all species compared. We used the essentiality data and a comparative genomics analysis reported here to highlight differences in essentiality between candidate Mtb drug targets and the corresponding Mk orthologs. Our findings suggest that the Mk genome encodes redundant or additional pathways that may confound validation of potential Mtb drugs and drug target candidates against the opportunistic pathogen. Additionally, we identified 57 intergenic regions containing four or more consecutive unoccupied TA sites. A disproportionally large number of these regions were located upstream of pe/ppe genes. Finally, we present an essentiality and orthology analysis of the Mk pRAW-like plasmid, pMK1248. IMPORTANCE Mk is one of the most common nontuberculous mycobacterial pathogens associated with tuberculosis-like pulmonary disease. Drug resistance emergence is a threat to the control of Mk infections, which already requires long-term, multidrug courses. A comprehensive understanding of Mk biology is critical to facilitate the development of new and more efficacious therapeutics against Mk. We combined transposon-based mutagenesis with analysis of insertion site identification data to uncover genes and other genomic regions required for Mk growth. We also compared the gene essentiality data set of Mk to those available for several other mycobacteria. This analysis highlighted key similarities and differences in the biology of Mk compared to these other species. Altogether, the genome-wide essentiality information generated and the results of the cross-species comparative genomics analysis represent valuable resources to assist the process of identifying and prioritizing potential Mk drug target candidates and to guide future studies on Mk biology.

Genome-wide screen identifies host loci that modulate Mycobacterium tuberculosis fitness in immunodivergent mice

Genetic differences among mammalian hosts and among strains of Mycobacterium tuberculosis (Mtb) are well-established determinants of tuberculosis (TB) patient outcomes. The advent of recombinant inbred mouse panels and next-generation transposon mutagenesis and sequencing approaches has enabled dissection of complex host-pathogen interactions. To identify host and pathogen genetic determinants of Mtb pathogenesis, we infected members of the highly diverse BXD family of strains with a comprehensive library of Mtb transposon mutants (TnSeq). Members of the BXD family segregate for Mtb-resistant C57BL/6J (B6 or B) and Mtb-susceptible DBA/2J (D2 or D) haplotypes. The survival of each bacterial mutant was quantified within each BXD host, and we identified those bacterial genes that were differentially required for Mtb fitness across BXD genotypes. Mutants that varied in survival among the host family of strains were leveraged as reporters of "endophenotypes," each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted quantitative trait loci (QTL) mapping of these bacterial fitness endophenotypes and identified 140 host-pathogen QTL (hpQTL). We located a QTL hotspot on chromosome 6 (75.97-88.58 Mb) associated with the genetic requirement of multiple Mtb genes: Rv0127 (mak), Rv0359 (rip2), Rv0955 (perM), and Rv3849 (espR). Together, this screen reinforces the utility of bacterial mutant libraries as precise reporters of the host immunological microenvironment during infection and highlights specific host-pathogen genetic interactions for further investigation. To enable downstream follow-up for both bacterial and mammalian genetic research communities, all bacterial fitness profiles have been deposited into GeneNetwork.org and added into the comprehensive collection of TnSeq libraries in MtbTnDB.

Identification of Genes Required for Swarming Motility in Bacillus subtilis Using Transposon Mutagenesis and High-Throughput Sequencing (TnSeq)

Bacillus subtilis exhibits swarming motility, a flagellar-mediated form of surface motility. Here, we use transposon mutagenesis and sequencing (TnSeq) to perform a high-throughput screen for candidate genes required for swarming. The TnSeq approach identified all of the known genes required for flagellar biosynthesis and nearly all of the previously reported regulators that promote swarming. Moreover, we identified an additional 36 genes that improve swarming and validated them individually. Among these, two mutants with severe defects were recovered, including fliT, required for flagellar biosynthesis, and a gene of unknown function, yolB, whose defect could not be attributed to a lack of flagella. In addition to discovering additional genes required for B. subtilis swarming, our work validates TnSeq as a powerful approach for comprehensively identifying genes important for nonessential processes such as colony expansion on plates. IMPORTANCE In TnSeq, transposons are randomly inserted throughout the chromosome at a population level, but insertions that disrupt genes of essential function cause strains that carry them to fall out of the population and appear underrepresented at the sequence level. Here, we apply TnSeq to the nonessential phenotype of motility in B. subtilis and spatially select for cells proficient in swarming. We find that insertions in nearly all genes previously identified as required for swarming are underrepresented in TnSeq analysis, and we identify 36 additional genes that enhance swarming. We demonstrate that TnSeq is a powerful tool for the genetic analysis of motility and likely other nonlethal screens for which enrichment is strong.

Analysis of Gene Essentiality from TnSeq Data Using Transit

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

Probing Differences in Gene Essentiality Between the Human and Animal Adapted Lineages of the Mycobacterium tuberculosis Complex Using TnSeq

Members of the Mycobacterium tuberculosis complex (MTBC) show distinct host adaptations, preferences and phenotypes despite being >99% identical at the nucleic acid level. Previous studies have explored gene expression changes between the members, however few studies have probed differences in gene essentiality. To better understand the functional impacts of the nucleic acid differences between Mycobacterium bovis and Mycobacterium tuberculosis, we used the Mycomar T7 phagemid delivery system to generate whole genome transposon libraries in laboratory strains of both species and compared the essentiality status of genes during growth under identical in vitro conditions. Libraries contained insertions in 54% of possible TA sites in M. bovis and 40% of those present in M. tuberculosis, achieving similar saturation levels to those previously reported for the MTBC. The distributions of essentiality across the functional categories were similar in both species. 527 genes were found to be essential in M. bovis whereas 477 genes were essential in M. tuberculosis and 370 essential genes were common in both species. CRISPRi was successfully utilised in both species to determine the impacts of silencing genes including wag31, a gene involved in peptidoglycan synthesis and Rv2182c/Mb2204c, a gene involved in glycerophospholipid metabolism. We observed species specific differences in the response to gene silencing, with the inhibition of expression of Mb2204c in M. bovis showing significantly less growth impact than silencing its orthologue (Rv2182c) in M. tuberculosis. Given that glycerophospholipid metabolism is a validated pathway for antimicrobials, our observations suggest that target vulnerability in the animal adapted lineages cannot be assumed to be the same as the human counterpart. This is of relevance for zoonotic tuberculosis as it implies that the development of antimicrobials targeting the human adapted lineage might not necessarily be effective against the animal adapted lineage. The generation of a transposon library and the first reported utilisation of CRISPRi in M. bovis will enable the use of these tools to further probe the genetic basis of survival under disease relevant conditions.

A Genome-Scale Antibiotic Screen in Serratia marcescens Identifies YdgH as a Conserved Modifier of Cephalosporin and Detergent Susceptibility

Serratia marcescens, a member of the order Enterobacterales, is adept at colonizing health care environments and is an important cause of invasive infections. Antibiotic resistance is a daunting problem in S. marcescens because, in addition to plasmid-mediated mechanisms, most isolates have considerable intrinsic resistance to multiple antibiotic classes. To discover endogenous modifiers of antibiotic susceptibility in S. marcescens, a high-density transposon insertion library was subjected to sub-MICs of two cephalosporins, cefoxitin, and cefepime, as well as the fluoroquinolone ciprofloxacin. Comparisons of transposon insertion abundance before and after antibiotic exposure identified hundreds of potential modifiers of susceptibility to these agents. Using single-gene deletions, we validated several candidate modifiers of cefoxitin susceptibility and chose ydgH, a gene of unknown function, for further characterization. In addition to cefoxitin, deletion of ydgH in S. marcescens resulted in decreased susceptibility to multiple third-generation cephalosporins and, in contrast, to increased susceptibility to both cationic and anionic detergents. YdgH is highly conserved throughout the Enterobacterales, and we observed similar phenotypes in Escherichia coli O157:H7 and Enterobacter cloacae mutants. YdgH is predicted to localize to the periplasm, and we speculate that it may be involved there in cell envelope homeostasis. Collectively, our findings provide insight into chromosomal mediators of antibiotic resistance in S. marcescens and will serve as a resource for further investigations of this important pathogen.

A Method of Transposon Insertion Sequencing in Comprehensively Identifying Vibrio vulnificus Genes Required for Growth in Human Serum

One of the most powerful approaches to detect the loci that enable a pathogen to cause disease is the creation of a high-density transposon mutant library by transposon insertion sequencing (TIS) and the screening of the library using an adequate in vivo and/or ex vivo model of the disease. Here we describe the procedure for detection of the putative loci required for a septicemic pathogen to cause sepsis in humans by using TIS plus an ex vivo model of septicaemia: to grow the pathogen in fresh and inactivated human serum. We selected V. vulnificus because it is a highly invasive pathogen capable of spreading from an infection site to the bloodstream, causing sepsis and death in less than 24 h. To survive and proliferate in blood (or host serum), the pathogen requires mechanisms to overcome the innate immune defenses and metabolic limitations of this host niche. Initially, genes under-represented for insertions can be used to estimate the V. vulnificus essential gene set. Analysis of the relative abundance of insertion mutants in the library after exposure to serum would detect which genes are essential for the pathogen to overcome the diverse limitations imposed by serum.

Modeling Site-Specific Nucleotide Biases Affecting Himar1 Transposon Insertion Frequencies in TnSeq Data Sets

TnSeq is a widely used methodology for determining gene essentiality, conditional fitness, and genetic interactions in bacteria. The Himar1 transposon is restricted to insertions at TA dinucleotides, but otherwise, few site-specific biases have been identified. As a result, most analytical approaches assume that insertions are expected to be randomly distributed among TA sites in nonessential regions. However, through analysis of Himar1 transposon libraries in Mycobacterium tuberculosis, we demonstrate that there are site-specific biases that affect the frequency of insertion of the Himar1 transposon at different TA sites. We use machine learning and statistical models to characterize patterns in the nucleotides surrounding TA sites that correlate with high or low insertion counts. We then develop a quantitative model based on these patterns that can be used to predict the expected counts at each TA site based on nucleotide context, which can explain up to half of the variance in insertion counts. We show that these insertion preferences exist in Himar1 TnSeq data sets from other mycobacterial and nonmycobacterial species. We present an improved method for identification of essential genes, called TTN-Fitness, that can better distinguish true biological fitness effects by comparing observed counts to expected counts based on our site-specific model of insertion preferences. Compared to previous essentiality methods, TTN-Fitness can make finer distinctions among genes whose disruption causes a fitness defect (or advantage), separating them out from the large pool of nonessentials, and is able to classify many smaller genes (with few TA sites) that were previously characterized as uncertain. IMPORTANCE When using the Himar1 transposon to create transposon insertion mutant libraries, it is known that the transposon is restricted to insertions at TA dinucleotide sites throughout the genome, and the absence of insertions is used to infer which genes are essential (or conditionally essential) in a bacterial organism. It is widely assumed that insertions in nonessential regions are otherwise random, and this assumption is used as the basis of several methods for statistical analysis of TnSeq data. In this paper, we show that the nucleotide sequence surrounding TA sites influences the magnitude of insertions, and these Himar1 insertion preferences (sequence biases) can partially explain why some sites have higher counts than others. We use this predictive model to make improved estimates of the fitness effects of genes, which help make finer distinctions of the phenotype and biological consequences of disruption of nonessential genes.

Is amplification bias consequential in transposon sequencing (TnSeq) assays? A case study with a Staphylococcus aureus TnSeq library subjected to PCR-based and amplification-free enrichment methods

As transposon sequencing (TnSeq) assays have become prolific in the microbiology field, it is of interest to scrutinize their potential drawbacks. TnSeq data consist of millions of nucleotide sequence reads that are generated by PCR amplification of transposon-genomic junctions. Reads mapping to the junctions are enumerated thus providing information on the number of transposon insertion mutations in each individual gene. Here we explore the possibility that PCR amplification of transposon insertions in a TnSeq library skews the results by introducing bias into the detection and/or enumeration of insertions. We compared the detection and frequency of mapped insertions when altering the number of PCR cycles, and when including a nested PCR, in the enrichment step. Additionally, we present nCATRAs - a novel, amplification-free TnSeq method where the insertions are enriched via CRISPR/Cas9-targeted transposon cleavage and subsequent Oxford Nanopore MinION sequencing. nCATRAs achieved 54 and 23% enrichment of the transposons and transposon-genomic junctions, respectively, over background genomic DNA. These PCR-based and PCR-free experiments demonstrate that, overall, PCR amplification does not significantly bias the results of TnSeq insofar as insertions in the majority of genes represented in our library were similarly detected regardless of PCR cycle number and whether or not PCR amplification was employed. However, the detection of a small subset of genes which had been previously described as essential is sensitive to the number of PCR cycles. We conclude that PCR-based enrichment of transposon insertions in a TnSeq assay is reliable, but researchers interested in profiling putative essential genes should carefully weigh the number of amplification cycles employed in their library preparation protocols. In addition, nCATRAs is comparable to traditional PCR-based methods (Kendall's correlation=0.896-0.897) although the latter remain superior owing to their accessibility and high sequencing depth.

Reproducible and accessible analysis of transposon insertion sequencing in Galaxy for qualitative essentiality analyses

BACKGROUND

Significant progress has been made in advancing and standardizing tools for human genomic and biomedical research. Yet, the field of next-generation sequencing (NGS) analysis for microorganisms (including multiple pathogens) remains fragmented, lacks accessible and reusable tools, is hindered by local computational resource limitations, and does not offer widely accepted standards. One such "problem areas" is the analysis of Transposon Insertion Sequencing (TIS) data. TIS allows probing of almost the entire genome of a microorganism by introducing random insertions of transposon-derived constructs. The impact of the insertions on the survival and growth under specific conditions provides precise information about genes affecting specific phenotypic characteristics. A wide array of tools has been developed to analyze TIS data. Among the variety of options available, it is often difficult to identify which one can provide a reliable and reproducible analysis.

RESULTS

Here we sought to understand the challenges and propose reliable practices for the analysis of TIS experiments. Using data from two recent TIS studies, we have developed a series of workflows that include multiple tools for data de-multiplexing, promoter sequence identification, transposon flank alignment, and read count repartition across the genome. Particular attention was paid to quality control procedures, such as determining the optimal tool parameters for the analysis and removal of contamination.

CONCLUSIONS

Our work provides an assessment of the currently available tools for TIS data analysis. It offers ready to use workflows that can be invoked by anyone in the world using our public Galaxy platform ( https://usegalaxy.org ). To lower the entry barriers, we have also developed interactive tutorials explaining details of TIS data analysis procedures at https://bit.ly/gxy-tis .

Quantification of Brucella abortus population structure in a natural host

Cattle are natural hosts of the intracellular pathogen Brucella abortus, which inflicts a significant burden on the health and reproduction of these important livestock. The primary routes of infection in field settings have been described, but it is not known how the bovine host shapes the structure of B. abortus populations during infection. We utilized a library of uniquely barcoded B. abortus strains to temporally and spatially quantify population structure during colonization of cattle through a natural route of infection. Introducing 10⁸ bacteria from this barcoded library to the conjunctival mucosa resulted in expected levels of local lymph node colonization at a 1-wk time point. We leveraged variance in strain abundance in the library to demonstrate that only 1 in 10,000 brucellae introduced at the site of infection reached a parotid lymph node. Thus, cattle restrict the overwhelming majority of B. abortus introduced via the ocular conjunctiva at this dose. Individual strains were spatially restricted within the host tissue, and the total B. abortus census was dominated by a small number of distinct strains in each lymph node. These results define a bottleneck that B. abortus must traverse to colonize local lymph nodes from the conjunctival mucosa. The data further support a model in which a small number of spatially isolated granulomas founded by unique strains are present at 1 wk postinfection. These experiments demonstrate the power of barcoded transposon tools to quantify infection bottlenecks and to define pathogen population structure in host tissues.

Genome-Wide Mutagenesis Identifies Factors Involved in Enterococcus faecalis Vaginal Adherence and Persistence

Enterococcus faecalis is a Gram-positive commensal bacterium native to the gastrointestinal tract and an opportunistic pathogen of increasing clinical concern. E. faecalis also colonizes the female reproductive tract, and reports suggest vaginal colonization increases following antibiotic treatment or in patients with aerobic vaginitis. Currently, little is known about specific factors that promote E. faecalis vaginal colonization and subsequent infection. We modified an established mouse vaginal colonization model to explore E. faecalis vaginal carriage and demonstrate that both vancomycin-resistant and -sensitive strains colonize the murine vaginal tract. Following vaginal colonization, we observed E. faecalis in vaginal, cervical, and uterine tissue. A mutant lacking endocarditis- and biofilm-associated pili (Ebp) exhibited a decreased ability to associate with human vaginal and cervical cells in vitro but did not contribute to colonization in vivo Thus, we screened a low-complexity transposon (Tn) mutant library to identify novel genes important for E. faecalis colonization and persistence in the vaginal tract. This screen revealed 383 mutants that were underrepresented during vaginal colonization at 1, 5, and 8 days postinoculation compared to growth in culture medium. We confirmed that mutants deficient in ethanolamine catabolism or in the type VII secretion system were attenuated in persisting during vaginal colonization. These results reveal the complex nature of vaginal colonization and suggest that multiple factors contribute to E. faecalis persistence in the reproductive tract.

Redefining genetic essentiality in Mycobacterium tuberculosis

Sequencing transposon mutant libraries have been pivotal in annotating essential and non-essential genes in bacteria. This is particularly very helpful in the case of Mycobacterium tuberculosis with a large part of its genome without known function. It is not known whether there are any variations in the essentiality states as a function of optimal growth in the absence of any selection pressure. We here grow a high-density mutant library of M. tuberculosis through serial cultures and monitor the temporal fluctuations in insertion frequencies across all TA dinucleotides in the genome. Genes that cause morphological and physiological heterogeneity or enable metabolic bypass were found to gradually lose insertions, while genes comprising the toxin-antitoxin systems were found to get enriched with insertions during growth in nutrient replete conditions. High levels of fluctuations were observed in genes involved in cell wall and cell processes, intermediary metabolism, and genes involved in virulence, suggesting new modes of adaptation undertaken by the mutants. We also report the essentiality status of several newly annotated genetic features.

An improved conjugation method for Pseudomonas syringae

In order to achieve saturating transposon mutagenesis of the genome of plant pathogenic strains of Pseudomonas syringae we needed to improve plasmid conjugation frequency. Manipulation of the growth stage of donor and recipient cells allowed the required increase in frequency and facilitated conjugation of otherwise recalcitrant strains.

Co-catabolism of arginine and succinate drives symbiotic nitrogen fixation

Biological nitrogen fixation emerging from the symbiosis between bacteria and crop plants holds promise to increase the sustainability of agriculture. One of the biggest hurdles for the engineering of nitrogen-fixing organisms is an incomplete knowledge of metabolic interactions between microbe and plant. In contrast to the previously assumed supply of only succinate, we describe here the CATCH-N cycle as a novel metabolic pathway that co-catabolizes plant-provided arginine and succinate to drive the energy-demanding process of symbiotic nitrogen fixation in endosymbiotic rhizobia. Using systems biology, isotope labeling studies and transposon sequencing in conjunction with biochemical characterization, we uncovered highly redundant network components of the CATCH-N cycle including transaminases that interlink the co-catabolism of arginine and succinate. The CATCH-N cycle uses N2 as an additional sink for reductant and therefore delivers up to 25% higher yields of nitrogen than classical arginine catabolism-two alanines and three ammonium ions are secreted for each input of arginine and succinate. We argue that the CATCH-N cycle has evolved as part of a synergistic interaction to sustain bacterial metabolism in the microoxic and highly acid environment of symbiosomes. Thus, the CATCH-N cycle entangles the metabolism of both partners to promote symbiosis. Our results provide a theoretical framework and metabolic blueprint for the rational design of plants and plant-associated organisms with new properties to improve nitrogen fixation.

The type IV pilin PilA couples surface attachment and cell-cycle initiation in Caulobacter crescentus

Understanding how bacteria colonize surfaces and regulate cell-cycle progression in response to cellular adhesion is of fundamental importance. Here, we use transposon sequencing in conjunction with fluorescence resonance energy transfer (FRET) microscopy to uncover the molecular mechanism for how surface sensing drives cell-cycle initiation in Caulobacter crescentus We identify the type IV pilin protein PilA as the primary signaling input that couples surface contact to cell-cycle initiation via the second messenger cyclic di-GMP (c-di-GMP). Upon retraction of pili filaments, the monomeric pilin reservoir in the inner membrane is sensed by the 17-amino acid transmembrane helix of PilA to activate the PleC-PleD two-component signaling system, increase cellular c-di-GMP levels, and signal the onset of the cell cycle. We termed the PilA signaling sequence CIP for "cell-cycle initiating pilin" peptide. Addition of the chemically synthesized CIP peptide initiates cell-cycle progression and simultaneously inhibits surface attachment. The broad conservation of the type IV pili and their importance in pathogens for host colonization suggests that CIP peptide mimetics offer strategies to inhibit surface sensing, prevent biofilm formation and control persistent infections.

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.

Genome-wide assessment of antimicrobial tolerance in Yersinia pseudotuberculosis under ciprofloxacin stress

Yersinia pseudotuberculosis is a Gram-negative bacterium capable of causing gastrointestinal infection and is closely related to the highly virulent plague bacillus Yersinia pestis. Infections by both species are currently treatable with antibiotics such as ciprofloxacin, a quinolone-class drug of major clinical importance in the treatment of many other infections. Our current understanding of the mechanism of action of ciprofloxacin is that it inhibits DNA replication by targeting DNA gyrase, and that resistance is primarily due to mutation of this target site, along with generic efflux and detoxification strategies. We utilized transposon-directed insertion site sequencing (TraDIS or TnSeq) to identify the non-essential chromosomal genes in Y. pseudotuberculosis that are required to tolerate sub-lethal concentrations of ciprofloxacin in vitro. As well as highlighting recognized antibiotic resistance genes, we provide evidence that multiple genes involved in regulating DNA replication and repair are central in enabling Y. pseudotuberculosis to tolerate the antibiotic, including DksA (yptb0734), a regulator of RNA polymerase, and Hda (yptb2792), an inhibitor of DNA replication initiation. We furthermore demonstrate that even at sub-lethal concentrations, ciprofloxacin causes severe cell-wall stress, requiring lipopolysaccharide lipid A, O-antigen and core biosynthesis genes to resist the sub-lethal effects of the antibiotic. It is evident that coping with the consequence(s) of antibiotic-induced stress requires the contribution of scores of genes that are not exclusively engaged in drug resistance.

Comprehensive Functional Analysis of the Enterococcus faecalis Core Genome Using an Ordered, Sequence-Defined Collection of Insertional Mutations in Strain OG1RF

The robust ability of Enterococcus faecalis to survive outside the host and to spread via oral-fecal transmission and its high degree of intrinsic and acquired antimicrobial resistance all complicate the treatment of hospital-acquired enterococcal infections. The conserved E. faecalis core genome serves as an important genetic scaffold for evolution of this bacterium in the modern health care setting and also provides interesting vaccine and drug targets. We used an innovative pooling/sequencing strategy to map a large collection of arrayed transposon insertions in E. faecalis OG1RF and generated an arrayed library of defined mutants covering approximately 70% of the OG1RF genome. Then, we performed high-throughput transposon sequencing experiments using this library to determine core genomic determinants of bile resistance in OG1RF. This collection is a valuable resource for comprehensive, functional enterococcal genomics using both traditional and high-throughput approaches and enables immediate recovery of mutants of interest.

Enterococcus faecalis is a common commensal bacterium in animal gastrointestinal (GI) tracts and a leading cause of opportunistic infections of humans in the modern health care setting. E. faecalis OG1RF is a plasmid-free strain that contains few mobile elements yet retains the robust survival characteristics, intrinsic antibiotic resistance, and virulence traits characteristic of most E. faecalis genotypes. To facilitate interrogation of the core enterococcal genetic determinants for competitive fitness in the GI tract, biofilm formation, intrinsic antimicrobial resistance, and survival in the environment, we generated an arrayed, sequence-defined set of chromosomal transposon insertions in OG1RF. We used an orthogonal pooling strategy in conjunction with Illumina sequencing to identify a set of mutants with unique, single Himar-based transposon insertions. The mutants contained insertions in 1,926 of 2,651 (72.6%) annotated open reading frames and in the majority of hypothetical protein-encoding genes and intergenic regions greater than 100 bp in length, which could encode small RNAs. As proof of principle of the usefulness of this arrayed transposon library, we created a minimal input pool containing 6,829 mutants chosen for maximal genomic coverage and used an approach that we term SMarT (sequence-defined marinertechnology) transposon sequencing (TnSeq) to identify numerous genetic determinants of bile resistance in E. faecalis OG1RF. These included several genes previously associated with bile acid resistance as well as new loci. Our arrayed library allows functional screening of a large percentage of the genome with a relatively small number of mutants, reducing potential effects of bottlenecking, and enables immediate recovery of mutants following competitions.

IMPORTANCE The robust ability of Enterococcus faecalis to survive outside the host and to spread via oral-fecal transmission and its high degree of intrinsic and acquired antimicrobial resistance all complicate the treatment of hospital-acquired enterococcal infections. The conserved E. faecalis core genome serves as an important genetic scaffold for evolution of this bacterium in the modern health care setting and also provides interesting vaccine and drug targets. We used an innovative pooling/sequencing strategy to map a large collection of arrayed transposon insertions in E. faecalis OG1RF and generated an arrayed library of defined mutants covering approximately 70% of the OG1RF genome. Then, we performed high-throughput transposon sequencing experiments using this library to determine core genomic determinants of bile resistance in OG1RF. This collection is a valuable resource for comprehensive, functional enterococcal genomics using both traditional and high-throughput approaches and enables immediate recovery of mutants of interest.

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.

Transposon Insertion Site Sequencing for Synthetic Lethal Screening

Transposon insertion site sequencing (TIS) permits genome-wide, quantitative fitness assessment of individual genomic loci. In addition to the identification of essential genes in given growth conditions, TIS enables the elucidation of genetic networks such as synthetic lethal or suppressor gene combinations. Therefore, TIS becomes an exceptionally powerful tool for the high-throughput determination of genotype-phenotype relationships in bacteria. Here, we describe a protocol for the generation of high-density transposon insertion libraries and subsequent preparation of DNA samples for Illumina sequencing using the Gram-negative bacterium Vibrio cholerae as an example.

Functional genomics of lipid metabolism in the oleaginous yeast Rhodosporidium toruloides

The basidiomycete yeast Rhodosporidium toruloides (also known as Rhodotorula toruloides) accumulates high concentrations of lipids and carotenoids from diverse carbon sources. It has great potential as a model for the cellular biology of lipid droplets and for sustainable chemical production. We developed a method for high-throughput genetics (RB-TDNAseq), using sequence-barcoded Agrobacterium tumefaciens T-DNA insertions. We identified 1,337 putative essential genes with low T-DNA insertion rates. We functionally profiled genes required for fatty acid catabolism and lipid accumulation, validating results with 35 targeted deletion strains. We identified a high-confidence set of 150 genes affecting lipid accumulation, including genes with predicted function in signaling cascades, gene expression, protein modification and vesicular trafficking, autophagy, amino acid synthesis and tRNA modification, and genes of unknown function. These results greatly advance our understanding of lipid metabolism in this oleaginous species and demonstrate a general approach for barcoded mutagenesis that should enable functional genomics in diverse fungi.

The fungus Rhodosporidium toruloides can grow on substances extracted from plant matter that is inedible to humans such as corn stalks, wood pulp, and grasses. Under some growth conditions, the fungus can accumulate massive stores of hydrocarbon-rich fats and pigments. A community of scientists and engineers has begun genetically modifying R. toruloides to convert these naturally produced fats and pigments into fuels, chemicals and medicines. These could form sustainable replacements for products made from petroleum or harvested from threatened animal and plant species.

Fungi, plants, animals and other eukaryotes store fat in specialized compartments called lipid droplets. The genes that control the metabolism – the production, use and storage – of fat in lipid bodies have been studied in certain eukaryotes, including species of yeast. However, R. toruloides is only distantly related to the most well-studied of these species. This means that we cannot be certain that a gene will play the same role in R. toruloides as in those species.

To assemble the most comprehensive list possible of the genes in R. toruloides that affect the production, use, or storage of fat in lipid bodies, Coradetti, Pinel et al. constructed a population of hundreds of thousands of mutant fungal strains, each with its own unique DNA ‘barcode’. The effects that mutations in over 6,000 genes had on growth and fat accumulation in these fungi were measured simultaneously in several experiments. This general approach is not new, but technical limitations had, until now, restricted its use in fungi to a few species.

Coradetti, Pinel et al. identified hundreds of genes that affected the ability of R. toruloides to metabolise fat. Many of these genes were related to genes with known roles in fat metabolism in other eukaryotes. Other genes are involved in different cell processes, such as the recycling of waste products in the cell. Their identification adds weight to the view that the links between these cellular processes and fat metabolism are deep and widespread amongst eukaryotes. Finally, some of the genes identified by Coradetti, Pinel et al. are not closely related to any well-studied genes. Further study of these genes could help us to understand why R. toruloides can accumulate much larger amounts of fat than most other fungi.

The methods developed by Coradetti, Pinel et al. should be possible to implement in many species of fungi. As a result these techniques may eventually contribute to the development of new treatments for human fungal diseases, the protection of important food crops, and a deeper understanding of the roles various fungi play in the broader ecosystem.

Transposon Sequencing Uncovers an Essential Regulatory Function of Phosphoribulokinase for Methylotrophy

Methylotrophy is the ability of organisms to grow at the expense of reduced one-carbon compounds, such as methanol or methane. Here, we used transposon sequencing combining hyper-saturated transposon mutagenesis with high-throughput sequencing to define the essential methylotrophy genome of Methylobacterium extorquens PA1, a model methylotroph. To distinguish genomic regions required for growth only on methanol from general required genes, we contrasted growth on methanol with growth on succinate, a non-methylotrophic reference substrate. About 500,000 insertions were mapped for each condition, resulting in a median insertion distance of five base pairs. We identified 147 genes and 76 genes as specific for growth on methanol and succinate, respectively, and a set of 590 genes as required under both growth conditions. For the integration of metabolic functions, we reconstructed a genome-scale metabolic model and performed in silico essentiality analysis. In total, the approach uncovered 95 genes not previously described as crucial for methylotrophy, including genes involved in respiration, carbon metabolism, transport, and regulation. Strikingly, regardless of the absence of the Calvin cycle in the methylotroph, the screen led to the identification of the gene for phosphoribulokinase as essential during growth on methanol, but not during growth on succinate. Genetic experiments in addition to metabolomics and proteomics revealed that phosphoribulokinase serves a key regulatory function. Our data support a model according to which ribulose-1,5-bisphosphate is an essential metabolite that induces a transcriptional regulator driving one-carbon assimilation.

Genome-Wide Survey of Pseudomonas aeruginosa PA14 Reveals a Role for the Glyoxylate Pathway and Extracellular Proteases in the Utilization of Mucin

Chronic airway infections by the opportunistic pathogen Pseudomonas aeruginosa are a major cause of mortality in cystic fibrosis (CF) patients. Although this bacterium has been extensively studied for its virulence determinants, biofilm growth, and immune evasion mechanisms, comparatively little is known about the nutrient sources that sustain its growth in vivo Respiratory mucins represent a potentially abundant bioavailable nutrient source, although we have recently shown that canonical pathogens inefficiently use these host glycoproteins as a growth substrate. However, given that P. aeruginosa, particularly in its biofilm mode of growth, is thought to grow slowly in vivo, the inefficient use of mucin glycoproteins may be relevant to its persistence within the CF airways. To this end, we used whole-genome fitness analysis, combining transposon mutagenesis with high-throughput sequencing, to identify genetic determinants required for P. aeruginosa growth using intact purified mucins as a sole carbon source. Our analysis reveals a biphasic growth phenotype, during which the glyoxylate pathway and amino acid biosynthetic machinery are required for mucin utilization. Secondary analyses confirmed the simultaneous liberation and consumption of acetate during mucin degradation and revealed a central role for the extracellular proteases LasB and AprA. Together, these studies describe a molecular basis for mucin-based nutrient acquisition by P. aeruginosa and reveal a host-pathogen dynamic that may contribute to its persistence within the CF airways.

Gene Transfer Agent Promotes Evolvability within the Fittest Subpopulation of a Bacterial Pathogen

The Bartonella gene transfer agent (BaGTA) is an archetypical example for domestication of a phage-derived element to permit high-frequency genetic exchange in bacterial populations. Here we used multiplexed transposon sequencing (TnSeq) and single-cell reporters to globally define the core components and transfer dynamics of BaGTA. Our systems-level analysis has identified inner- and outer-circle components of the BaGTA system, including 55 regulatory components, as well as an additional 74 and 107 components mediating donor transfer and recipient uptake functions. We show that the stringent response signal guanosine-tetraphosphate (ppGpp) restricts BaGTA induction to a subset of fast-growing cells, whereas BaGTA particle uptake depends on a functional Tol-Pal trans-envelope complex that mediates outer-membrane invagination upon cell division. Our findings suggest that Bartonella evolved an efficient strategy to promote genetic exchange within the fittest subpopulation while disfavoring exchange of deleterious genetic information, thereby facilitating genome integrity and rapid host adaptation.

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.

Comparison of inherently essential genes of Porphyromonas gingivalis identified in two transposon-sequencing libraries

Porphyromonas gingivalis is a Gram-negative anaerobe and keystone periodontal pathogen. A mariner transposon insertion mutant library has recently been used to define 463 genes as putatively essential for the in vitro growth of P. gingivalis ATCC 33277 in planktonic culture (Library 1). We have independently generated a transposon insertion mutant library (Library 2) for the same P. gingivalis strain and herein compare genes that are putatively essential for in vitro growth in complex media, as defined by both libraries. In all, 281 genes (61%) identified by Library 1 were common to Library 2. Many of these common genes are involved in fundamentally important metabolic pathways, notably pyrimidine cycling as well as lipopolysaccharide, peptidoglycan, pantothenate and coenzyme A biosynthesis, and nicotinate and nicotinamide metabolism. Also in common are genes encoding heat-shock protein homologues, sigma factors, enzymes with proteolytic activity, and the majority of sec-related protein export genes. In addition to facilitating a better understanding of critical physiological processes, transposon-sequencing technology has the potential to identify novel strategies for the control of P. gingivalis infections. Those genes defined as essential by two independently generated TnSeq mutant libraries are likely to represent particularly attractive therapeutic targets.

Transposon insertion mapping with PIMMS - Pragmatic Insertional Mutation Mapping System

The PIMMS (Pragmatic Insertional Mutation Mapping System) pipeline has been developed for simple conditionally essential genome discovery experiments in bacteria. Capable of using raw sequence data files alongside a FASTA sequence of the reference genome and GFF file, PIMMS will generate a tabulated output of each coding sequence with corresponding mapped insertions accompanied with normalized results enabling streamlined analysis. This allows for a quick assay of the genome to identify conditionally essential genes on a standard desktop computer prioritizing results for further investigation.

Availability: The PIMMS script, manual and accompanying test data is freely available at https://github.com/ADAC-UoN/PIMMS