A comparison of dense transposon insertion libraries in the Salmonella serovars Typhi and Typhimurium

High-throughput fitness experiments reveal specific vulnerabilities of human-adapted Salmonella during stress and infection

Salmonella enterica is comprised of genetically distinct 'serovars' that together provide an intriguing model for exploring the genetic basis of pathogen evolution. Although the genomes of numerous Salmonella isolates with broad variations in host range and human disease manifestations have been sequenced, the functional links between genetic and phenotypic differences among these serovars remain poorly understood. Here, we conduct high-throughput functional genomics on both generalist (Typhimurium) and human-restricted (Typhi and Paratyphi A) Salmonella at unprecedented scale in the study of this enteric pathogen. Using a comprehensive systems biology approach, we identify gene networks with serovar-specific fitness effects across 25 host-associated stresses encountered at key stages of human infection. By experimentally perturbing these networks, we characterize previously undescribed pseudogenes in human-adapted Salmonella. Overall, this work highlights specific vulnerabilities encoded within human-restricted Salmonella that are linked to the degradation of their genomes, shedding light into the evolution of this enteric pathogen.

Identifying the suite of genes central to swimming in the biocontrol bacterium Pseudomonas protegens Pf-5

Swimming motility is a key bacterial trait, important to success in many niches. Biocontrol bacteria, such as Pseudomonas protegens Pf-5, are increasingly used in agriculture to control crop diseases, where motility is important for colonization of the plant rhizosphere. Swimming motility typically involves a suite of flagella and chemotaxis genes, but the specific gene set employed for both regulation and biogenesis can differ substantially between organisms. Here we used transposon-directed insertion site sequencing (TraDIS), a genome-wide approach, to identify 249 genes involved in P. protegens Pf-5 swimming motility. In addition to the expected flagella and chemotaxis, we also identified a suite of additional genes important for swimming, including genes related to peptidoglycan turnover, O-antigen biosynthesis, cell division, signal transduction, c-di-GMP turnover and phosphate transport, and 27 conserved hypothetical proteins. Gene knockout mutants and TraDIS data suggest that defects in the Pst phosphate transport system lead to enhanced swimming motility. Overall, this study expands our knowledge of pseudomonad motility and highlights the utility of a TraDIS-based approach for analysing the functions of thousands of genes. This work sets a foundation for understanding how swimming motility may be related to the inconsistency in biocontrol bacteria performance in the field.

Defying the odds: Determinants of the antimicrobial response of Salmonella Typhi and their interplay

Salmonella Typhi, the invasive serovar of S. enterica subspecies enterica, causes typhoid fever in healthy human hosts. The emergence of antibiotic-resistant strains has consistently challenged the successful treatment of typhoid fever with conventional antibiotics. Antimicrobial resistance (AMR) in Salmonella is acquired either by mutations in the genomic DNA or by acquiring extrachromosomal DNA via horizontal gene transfer. In addition, Salmonella can form a subpopulation of antibiotic persistent (AP) cells that can survive at high concentrations of antibiotics. These have reduced the effectiveness of the first and second lines of antibiotics used to treat Salmonella infection. The recurrent and chronic carriage of S. Typhi in human hosts further complicates the treatment process, as a remarkable shift in the immune response from pro-inflammatory Th1 to anti-inflammatory Th2 is observed. Recent studies have also highlighted the overlap between AP, persistent infection (PI) and AMR. These incidents have revealed several areas of research. In this review, we have put forward a timeline for the evolution of antibiotic resistance in Salmonella and discussed the different mechanisms of the same availed by the pathogen at the genotypic and phenotypic levels. Further, we have presented a detailed discussion on Salmonella antibiotic persistence (AP), PI, the host and bacterial virulence factors that can influence PI, and how both AP and PI can lead to AMR.

Identification of essential genes in Coxiella burnetii

Coxiella burnetii is an intracellular pathogen responsible for causing Q fever in humans, a disease with varied presentations ranging from a mild flu-like sickness to a debilitating illness that can result in endocarditis. The intracellular lifestyle of C. burnetii is unique, residing in an acidic phagolysosome-like compartment within host cells. An understanding of the core molecular biology of C. burnetii will greatly increase our understanding of C. burnetii growth, survival and pathogenesis. We used transposon-directed insertion site sequencing (TraDIS) to reveal C. burnetii Nine Mile Phase II genes fundamental for growth and in vitro survival. Screening a transposon library containing >10 000 unique transposon mutants revealed 512 predicted essential genes. Essential routes of synthesis were identified for the mevalonate pathway, as well as peptidoglycan and biotin synthesis. Some essential genes identified (e.g. predicted type IV secretion system effector genes) are typically considered to be associated with C. burnetii virulence, a caveat concerning the axenic media used in the study. Investigation into the conservation of the essential genes identified revealed that 78 % are conserved across all C. burnetii strains sequenced to date, which probably play critical functions. This is the first report of a whole genome transposon screen in C. burnetii that has been undertaken for the identification of essential genes.

Cre/lox-Mediated CRISPRi Library Reveals Core Genome of a Type I Methanotroph Methylotuvimicrobium buryatense 5GB1C

Methanotrophs play key roles in global methane cycling and are promising platforms for methane bioconversion. However, major gaps existing in fundamental knowledge undermines understanding of these methane-consuming microorganisms. To associate genes with a phenotype at the genome-wide level, we developed a Cre/lox-mediated method for constructing a large-scale CRISPRi library in a model methanotroph Methylotuvimicrobium buryatense 5GB1C. The efficiency of this Cre mediated integration method was up to a level of 10⁵ CFU/μg DNA. Targeting 4,100 predicted protein-coding genes, our CRISPRi pooled screening uncovered 788 core genes for the growth of strain 5GB1C using methane. The core genes are highly consistent with the gene knockout results, indicating the reliability of the CRISPRi screen. Insights from the core genes include that annotated isozymes generally exist in metabolic pathways and many core genes are hypothetical genes. This work not only provides functional genomic data for both fundamental research and metabolic engineering of methanotrophs, but also offers a method for CRISPRi library construction. IMPORTANCE Due to their key role in methane cycling and their industrial potential, methanotrophs have drawn increasing attention. Genome-wide experimental approaches for gene-phenotype mapping accelerate our understanding and engineering of a bacterium. However, these approaches are still unavailable in methanotrophs. This work has two significant implications. First, the core genes identified here provide functional genetic basics for complete reconstruction of the metabolic network and afford more clues for knowledge gaps. Second, the Cre-mediated knock-in method developed in this work enables large-scale DNA library construction in methanotrophs; the CRISPRi library can be used to screen the genes associated with special culture conditions.

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

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

Dual transcriptome based reconstruction of Salmonella-human integrated metabolic network to screen potential drug targets

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a highly adaptive pathogenic bacteria with a serious public health concern due to its increasing resistance to antibiotics. Therefore, identification of novel drug targets for S. Typhimurium is crucial. Here, we first created a pathogen-host integrated genome-scale metabolic network by combining the metabolic models of human and S. Typhimurium, which we further tailored to the pathogenic state by the integration of dual transcriptome data. The integrated metabolic model enabled simultaneous investigation of metabolic alterations in human cells and S. Typhimurium during infection. Then, we used the tailored pathogen-host integrated genome-scale metabolic network to predict essential genes in the pathogen, which are candidate novel drug targets to inhibit infection. Drug target prioritization procedure was applied to these targets, and pabB was chosen as a putative drug target. It has an essential role in 4-aminobenzoic acid (PABA) synthesis, which is an essential biomolecule for many pathogens. A structure based virtual screening was applied through docking simulations to predict candidate compounds that eliminate S. Typhimurium infection by inhibiting pabB. To our knowledge, this is the first comprehensive study for predicting drug targets and drug like molecules by using pathogen-host integrated genome-scale models, dual RNA-seq data and structure-based virtual screening protocols. This framework will be useful in proposing novel drug targets and drugs for antibiotic-resistant pathogens.

Genome-Wide Investigation of Pasteurella multocida Identifies the Stringent Response as a Negative Regulator of Hyaluronic Acid Capsule Production

The bacterial pathogen P. multocida can cause serious disease in production animals, including fowl cholera in poultry, hemorrhagic septicemia in cattle and buffalo, atrophic rhinitis in pigs, and respiratory diseases in a range of livestock. P. multocida produces a capsule that is essential for systemic disease, but the complete mechanisms underlying synthesis and regulation of capsule production are not fully elucidated. A whole-genome analysis using TraDIS was undertaken to identify genes essential for growth in rich media and to obtain a comprehensive characterization of capsule production.

New Method for Genome-Scale Functional Genomic Study in Bacteria with Superior Performance: CRISPR Interference Screen

High-throughput genetic screens based on CRISPR/Cas9 technology are powerful tools to genome-wide identify gene function and genotype-phenotype association. Here, we describe a detailed protocol for conducting and evaluating pooled CRISPR screens interfering with gene expression in Escherichia coli. We provide step-by-step instructions for guide RNA library design and construction, genome-scale screening and next-generation sequencing data processing. This tool outperforms transposon sequencing (Tn-seq) with similar library sizes and short gene length. The workflow can be used in follow-up studies implemented in other bacteria systems.

Genetic Determinants in Salmonella enterica Serotype Typhimurium Required for Overcoming In Vitro Stressors in the Mimicking Host Environment

Salmonella enterica serotype Typhimurium, a nontyphoidal Salmonella (NTS), results in a range of enteric diseases, representing a major disease burden worldwide. There is still a significant portion of Salmonella genes whose mechanistic basis to overcome host innate defense mechanisms largely remains unknown. Here, we have applied transposon insertion sequencing (Tn-seq) method to unveil the genetic factors required for the growth or survival of S. Typhimurium under various host stressors simulated in vitro. A highly saturating Tn5 library of S. Typhimurium 14028s was subjected to selection during growth in the presence of short-chain fatty acid (100 mM propionate), osmotic stress (3% NaCl), or oxidative stress (1 mM H2O2) or survival in extreme acidic pH (30 min in pH 3) or starvation (12 days in 1× phosphate-buffered saline [PBS]). We have identified a total of 339 conditionally essential genes (CEGs) required to overcome at least one of these conditions mimicking host insults. Interestingly, all eight genes encoding FoF1-ATP synthase subunit proteins were required for fitness in all five stresses. Intriguingly, a total of 88 genes in Salmonella pathogenicity islands (SPI), including SPI-1, SPI-2, SPI-3, SPI-5, SPI-6, and SPI-11, are also required for fitness under the in vitro conditions. Additionally, by comparative analysis of the genes identified in this study and the genes previously shown to be required for in vivo fitness, we identified novel genes (marBCT, envF, barA, hscA, rfaQ, rfbI, and the genes encoding putative proteins STM14_1138, STM14_3334, STM14_4825, and STM_5184) that have compelling potential for the development of vaccines and antibacterial drugs to curb Salmonella infection. IMPORTANCE Salmonella enterica serotype Typhimurium is a major human bacterial pathogen that enters the food chain through meat animals asymptomatically carrying this pathogen. Despite the rich genome sequence data, a significant portion of Salmonella genes remain to be characterized for their potential contributions to virulence. In this study, we used transposon insertion sequencing (Tn-seq) to elucidate the genetic factors required for growth or survival under various host stressors, including short-chain fatty acids, osmotic stress, oxidative stress, extreme acid, and starvation. Among the total of 339 conditionally essential genes (CEGs) that are required under at least one of these five stress conditions were 221 previously known virulence genes required for in vivo fitness during infection in at least one of four animal species, including mice, chickens, pigs, and cattle. This comprehensive map of virulence phenotype-genotype in S. Typhimurium provides a roadmap for further interrogation of the biological functions encoded by the genome of this important human pathogen to survive in hostile host environments.

Essential Genes of Vibrio anguillarum and Other Vibrio spp. Guide the Development of New Drugs and Vaccines

Essential genes in bacterial pathogens are potential drug targets and vaccine candidates because disrupting their function is lethal. The development of new antibiotics, in addition to effective prevention measures such as vaccination, contributes to addressing the global problem of bacterial antibiotic resistance. The aim of this present study was to determine the essential genes of Vibrio anguillarum, a bacterial pathogen of aquatic animals, as a means to identify putative targets for novel drugs and to assist the prioritisation of potential vaccine candidates. Essential genes were characterised by a Tn-seq approach using the TnSC189 mariner transposon to construct a library of 52,662 insertion mutants. In total, 329 essential genes were identified, with 34.7% found within the core genome of this species; each of these genes represents a strong potential drug target. Seven essential gene products were predicted to reside in the cell membrane or be released extracellularly, thus serving as putative vaccine candidates. Comparison to essential gene data from five other studies of Vibrio species revealed 13 proteins to be conserved across the studies, while 25 genes were specific to V. anguillarum and not found to be essential in the other Vibrio spp. This study provides new information on the essential genes of Vibrio species and the methodology may be applied to other pathogens to guide the development of new drugs and vaccines, which will assist efforts to counter antibiotic resistance.

Comparative genome analysis of Salmonella enterica serovar Gallinarum biovars Pullorum and Gallinarum decodes strain specific genes

Salmonella enterica serovar Gallinarum biovar Pullorum (bvP) and biovar Gallinarum (bvG) are the etiological agents of pullorum disease (PD) and fowl typhoid (FT) respectively, which cause huge economic losses to poultry industry especially in developing countries including India. Vaccination and biosecurity measures are currently being employed to control and reduce the S. Gallinarum infections. High endemicity, poor implementation of hygiene and lack of effective vaccines pose challenges in prevention and control of disease in intensively maintained poultry flocks. Comparative genome analysis unravels similarities and dissimilarities thus facilitating identification of genomic features that aids in pathogenesis, niche adaptation and in tracing of evolutionary history. The present investigation was carried out to assess the genotypic differences amongst S.enterica serovar Gallinarum strains including Indian strain S. Gallinarum Sal40 VTCCBAA614. The comparative genome analysis revealed an open pan-genome consisting of 5091 coding sequence (CDS) with 3270 CDS belonging to core-genome, 1254 CDS to dispensable genome and strain specific genes i.e. singletons ranging from 3 to 102 amongst the analyzed strains. Moreover, the investigated strains exhibited diversity in genomic features such as virulence factors, genomic islands, prophage regions, toxin-antitoxin cassettes, and acquired antimicrobial resistance genes. Core genome identified in the study can give important leads in the direction of design of rapid and reliable diagnostics, and vaccine design for effective infection control as well as eradication. Additionally, the identified genetic differences among the S. enterica serovar Gallinarum strains could be used for bacterial typing, structure based inhibitor development by future experimental investigations on the data generated.

SorTn-seq: a high-throughput functional genomics approach to discovering regulators of bacterial gene expression

We recently developed a high-throughput functional genomics approach, named 'SorTn-seq', to identify factors affecting expression of any gene of interest in bacteria. Our approach facilitates high-throughput screening of complex mutant pools, a task previously hindered by a lack of suitable techniques. SorTn-seq combines high-density, Tn5-like transposon mutagenesis with fluorescence-activated cell sorting of a strain harboring a promoter-fluorescent reporter fusion, to isolate mutants with altered gene expression. The transposon mutant pool is sorted into different bins on the basis of fluorescence, and mutants are deep-sequenced to identify transposon insertions. DNA is prepared for sequencing by using commercial kits augmented with custom primers, enhancing ease of use and reproducibility. Putative regulators are identified by comparing the number of insertions per genomic feature in the different sort bins, by using existing bioinformatic pipelines and software packages. SorTn-seq can be completed in 1-2 weeks and requires general microbiology skills and basic flow cytometry experience.

Functional analysis of colonization factor antigen I positive enterotoxigenic Escherichia coli identifies genes implicated in survival in water and host colonization

Enterotoxigenic Escherichia coli (ETEC) expressing the colonization pili CFA/I are common causes of diarrhoeal infections in humans. Here, we use a combination of transposon mutagenesis and transcriptomic analysis to identify genes and pathways that contribute to ETEC persistence in water environments and colonization of a mammalian host. ETEC persisting in water exhibit a distinct RNA expression profile from those growing in richer media. Multiple pathways were identified that contribute to water survival, including lipopolysaccharide biosynthesis and stress response regulons. The analysis also indicated that ETEC growing in vivo in mice encounter a bottleneck driving down the diversity of colonizing ETEC populations.

Elucidating Essential Genes in Plant-Associated Pseudomonas protegens Pf-5 Using Transposon Insertion Sequencing

Gene essentiality studies have been performed on numerous bacterial pathogens, but essential gene sets have been determined for only a few plant-associated bacteria. Pseudomonas protegens Pf-5 is a plant-commensal, biocontrol bacterium that can control disease-causing pathogens on a wide range of crops. Work on Pf-5 has mostly focused on secondary metabolism and biocontrol genes, but genome-wide approaches such as high-throughput transposon mutagenesis have not yet been used for this species. In this study, we generated a dense P. protegens Pf-5 transposon mutant library and used transposon-directed insertion site sequencing (TraDIS) to identify 446 genes essential for growth on rich media. Genes required for fundamental cellular machinery were enriched in the essential gene set, while genes related to nutrient biosynthesis, stress responses, and transport were underrepresented. The majority of Pf-5 essential genes were part of the P. protegens core genome. Comparison of the essential gene set of Pf-5 with those of two plant-associated pseudomonads, P. simiae and P. syringae, and the well-studied opportunistic human pathogen P. aeruginosa PA14 showed that the four species share a large number of essential genes, but each species also had uniquely essential genes. Comparison of the Pf-5 in silico-predicted and in vitro-determined essential gene sets highlighted the essential cellular functions that are over- and underestimated by each method. Expanding essentiality studies into bacteria with a range of lifestyles may improve our understanding of the biological processes important for bacterial survival and growth.IMPORTANCE Essential genes are those crucial for survival or normal growth rates in an organism. Essential gene sets have been identified in numerous bacterial pathogens but only a few plant-associated bacteria. Employing genome-wide approaches, such as transposon insertion sequencing, allows for the concurrent analyses of all genes of a bacterial species and rapid determination of essential gene sets. We have used transposon insertion sequencing to systematically analyze thousands of Pseudomonas protegens Pf-5 genes and gain insights into gene functions and interactions that are not readily available using traditional methods. Comparing Pf-5 essential genes with those of three other pseudomonads highlights how gene essentiality varies between closely related species.

DEG 15, an update of the Database of Essential Genes that includes built-in analysis tools

Essential genes refer to genes that are required by an organism to survive under specific conditions. Studies of the minimal-gene-set for bacteria have elucidated fundamental cellular processes that sustain life. The past five years have seen a significant progress in identifying human essential genes, primarily due to the successful use of CRISPR/Cas9 in various types of human cells. DEG 15, a new release of the Database of Essential Genes (www.essentialgene.org), has provided major advancements, compared to DEG 10. Specifically, the number of eukaryotic essential genes has increased by more than fourfold, and that of prokaryotic ones has more than doubled. Of note, the human essential-gene number has increased by more than tenfold. Moreover, we have developed built-in analysis modules by which users can perform various analyses, such as essential-gene distributions between bacterial leading and lagging strands, sub-cellular localization distribution, enrichment analysis of gene ontology and KEGG pathways, and generation of Venn diagrams to compare and contrast gene sets between experiments. Additionally, the database offers customizable BLAST tools for performing species- and experiment-specific BLAST searches. Therefore, DEG comprehensively harbors updated human-curated essential-gene records among prokaryotes and eukaryotes with built-in tools to enhance essential-gene analysis.

Fundamental differences in physiology of Bordetella pertussis dependent on the two-component system Bvg revealed by gene essentiality studies

The identification of genes essential for a bacterium's growth reveals much about its basic physiology under different conditions. Bordetella pertussis, the causative agent of whooping cough, adopts both virulent and avirulent states through the activity of the two-component system, Bvg. The genes essential for B. pertussis growth in vitro were defined using transposon sequencing, for different Bvg-determined growth states. In addition, comparison of the insertion indices of each gene between Bvg phases identified those genes whose mutation exerted a significantly different fitness cost between phases. As expected, many of the genes identified as essential for growth in other bacteria were also essential for B. pertussis. However, the essentiality of some genes was dependent on Bvg. In particular, a number of key cell wall biosynthesis genes, including the entire mre/mrd locus, were essential for growth of the avirulent (Bvg minus) phase but not the virulent (Bvg plus) phase. In addition, cell wall biosynthesis was identified as a fundamental process that when disrupted produced greater fitness costs for the Bvg minus phase compared to the Bvg plus phase. Bvg minus phase growth was more susceptible than Bvg plus phase growth to the cell wall-disrupting antibiotic ampicillin, demonstrating the increased susceptibility of the Bvg minus phase to disruption of cell wall synthesis. This Bvg-dependent conditional essentiality was not due to Bvg-regulation of expression of cell wall biosynthesis genes; suggesting that this fundamental process differs between the Bvg phases in B. pertussis and is more susceptible to disruption in the Bvg minus phase. The ability of a bacterium to modify its cell wall synthesis is important when considering the action of antibiotics, particularly if developing novel drugs targeting cell wall synthesis.

Environmental flexibility does not explain metabolic robustness

Cells show remarkable resilience against genetic and environmental perturbations. However, its evolutionary origin remains obscure. In order to leverage methods of systems biology for examining cellular robustness, a computationally accessible way of quantification is needed. Here, we present an unbiased metric of structural robustness in genome-scale metabolic models based on concepts prevalent in reliability engineering and fault analysis. The probability of failure (PoF) is defined as the (weighted) portion of all possible combinations of loss-of-function mutations that disable network functionality. It can be exactly determined if all essential reactions, synthetic lethal pairs of reactions, synthetic lethal triplets of reactions etc. are known. In theory, these minimal cut sets (MCSs) can be calculated for any network, but for large models the problem remains computationally intractable. Herein, we demonstrate that even at the genome scale only the lowest-cardinality MCSs are required to efficiently approximate the PoF with reasonable accuracy. Building on an improved theoretical understanding, we analysed the robustness of 489 E. coli, Shigella, Salmonella, and fungal genome-scale metabolic models (GSMMs). In contrast to the popular "congruence theory", which explains the origin of genetic robustness as a byproduct of selection for environmental flexibility, we found no correlation between network robustness and the diversity of growth-supporting environments. On the contrary, our analysis indicates that amino acid synthesis rather than carbon metabolism dominates metabolic robustness.

Functional genomics reveals the toxin-antitoxin repertoire and AbiE activity in Serratia

Bacteriophage defences are divided into innate and adaptive systems. Serratia sp. ATCC 39006 has three CRISPR-Cas adaptive immune systems, but its innate immune repertoire is unknown. Here, we re-sequenced and annotated the Serratia genome and predicted its toxin-antitoxin (TA) systems. TA systems can provide innate phage defence through abortive infection by causing infected cells to 'shut down', limiting phage propagation. To assess TA system function on a genome-wide scale, we utilized transposon insertion and RNA sequencing. Of the 32 TA systems predicted bioinformatically, 4 resembled pseudogenes and 11 were demonstrated to be functional based on transposon mutagenesis. Three functional systems belonged to the poorly characterized but widespread, AbiE, abortive infection/TA family. AbiE is a type IV TA system with a predicted nucleotidyltransferase toxin. To investigate the mode of action of this toxin, we measured the transcriptional response to AbiEii expression. We observed dysregulated levels of tRNAs and propose that the toxin targets tRNAs resulting in bacteriostasis. A recent report on a related toxin shows this occurs through addition of nucleotides to tRNA(s). This study has demonstrated the utility of functional genomics for probing TA function in a high-throughput manner, defined the TA repertoire in Serratia and shown the consequences of AbiE induction.

Bacterial Genetic Engineering by Means of Recombineering for Reverse Genetics

Serving a robust platform for reverse genetics enabling the in vivo study of gene functions primarily in enterobacteriaceae, recombineering -or recombination-mediated genetic engineering-represents a powerful and relative straightforward genetic engineering tool. Catalyzed by components of bacteriophage-encoded homologous recombination systems and only requiring short ∼40–50 base homologies, the targeted and precise introduction of modifications (e.g., deletions, knockouts, insertions and point mutations) into the chromosome and other episomal replicons is empowered. Furthermore, by its ability to make use of both double- and single-stranded linear DNA editing substrates (e.g., PCR products or oligonucleotides, respectively), lengthy subcloning of specific DNA sequences is circumvented. Further, the more recent implementation of CRISPR-associated endonucleases has allowed for more efficient screening of successful recombinants by the selective purging of non-edited cells, as well as the creation of markerless and scarless mutants. In this review we discuss various recombineering strategies to promote different types of gene modifications, how they are best applied, and their possible pitfalls.

Genomic Profiling Reveals Distinct Routes To Complement Resistance in Klebsiella pneumoniae

The serum complement system is a first line of defense against bacterial invaders. Resistance to killing by serum enhances the capacity of Klebsiella pneumoniae to cause infection, but it is an incompletely understood virulence trait. Identifying and characterizing the factors responsible for preventing activation of, and killing by, serum complement could inform new approaches to treatment of K. pneumoniae infections. Here, we used functional genomic profiling to define the genetic basis of complement resistance in four diverse serum-resistant K. pneumoniae strains (NTUH-K2044, B5055, ATCC 43816, and RH201207), and explored their recognition by key complement components.

The serum complement system is a first line of defense against bacterial invaders. Resistance to killing by serum enhances the capacity of Klebsiella pneumoniae to cause infection, but it is an incompletely understood virulence trait. Identifying and characterizing the factors responsible for preventing activation of, and killing by, serum complement could inform new approaches to treatment of K. pneumoniae infections. Here, we used functional genomic profiling to define the genetic basis of complement resistance in four diverse serum-resistant K. pneumoniae strains (NTUH-K2044, B5055, ATCC 43816, and RH201207), and explored their recognition by key complement components. More than 90 genes contributed to resistance in one or more strains, but only three, rfaH, lpp, and arnD, were common to all four strains. Deletion of the antiterminator rfaH, which controls the expression of capsule and O side chains, resulted in dramatic complement resistance reductions in all strains. The murein lipoprotein gene lpp promoted capsule retention through a mechanism dependent on its C-terminal lysine residue; its deletion led to modest reductions in complement resistance. Binding experiments with the complement components C3b and C5b-9 showed that the underlying mechanism of evasion varied in the four strains: B5055 and NTUH-K2044 appeared to bypass recognition by complement entirely, while ATCC 43816 and RH201207 were able to resist killing despite being associated with substantial levels of C5b-9. All rfaH and lpp mutants bound C3b and C5b-9 in large quantities. Our findings show that, even among this small selection of isolates, K. pneumoniae adopts differing mechanisms and utilizes distinct gene sets to avoid complement attack.

A decade of advances in transposon-insertion sequencing

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

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

Selection-Driven Gene Inactivation in Salmonella

Bacterial genes are sometimes found to be inactivated by mutation. This inactivation may be observable simply because selection for function is intermittent or too weak to eliminate inactive alleles quickly. Here, I investigate cases in Salmonella enterica where inactivation is instead positively selected. These are identified by a rate of introduction of premature stop codons to a gene that is higher than expected under selective neutrality, as assessed by comparison to the rate of synonymous changes. I identify 84 genes that meet this criterion at a 10% false discovery rate. Many of these genes are involved in virulence, motility and chemotaxis, biofilm formation, and resistance to antibiotics or other toxic substances. It is hypothesized that most of these genes are subject to an ongoing process in which inactivation is favored under rare conditions, but the inactivated allele is deleterious under most other conditions and is subsequently driven to extinction by purifying selection.

TraDIS-Xpress: a high-resolution whole-genome assay identifies novel mechanisms of triclosan action and resistance

Understanding the genetic basis for a phenotype is a central goal in biological research. Much has been learnt about bacterial genomes by creating large mutant libraries and looking for conditionally important genes. However, current genome-wide methods are largely unable to assay essential genes which are not amenable to disruption. To overcome this limitation, we developed a new version of "TraDIS" (transposon directed insertion-site sequencing) that we term "TraDIS-Xpress" that combines an inducible promoter into the transposon cassette. This allows controlled overexpression and repression of all genes owing to saturation of inserts adjacent to all open reading frames as well as conventional inactivation. We applied TraDIS-Xpress to identify responses to the biocide triclosan across a range of concentrations. Triclosan is endemic in modern life, but there is uncertainty about its mode of action with a concentration-dependent switch from bacteriostatic to bactericidal action unexplained. Our results show a concentration-dependent response to triclosan with different genes important in survival between static and cidal exposures. These genes include those previously reported to have a role in triclosan resistance as well as a new set of genes, including essential genes. Novel genes identified as being sensitive to triclosan exposure include those involved in barrier function, small molecule uptake, and integrity of transcription and translation. We anticipate the approach we show here, by allowing comparisons across multiple experimental conditions of TraDIS data, and including essential genes, will be a starting point for future work examining how different drug conditions impact bacterial survival mechanisms.

The fitness landscape of the African Salmonella Typhimurium ST313 strain D23580 reveals unique properties of the pBT1 plasmid

We have used a transposon insertion sequencing (TIS) approach to establish the fitness landscape of the African Salmonella enterica serovar Typhimurium ST313 strain D23580, to complement our previous comparative genomic and functional transcriptomic studies. We used a genome-wide transposon library with insertions every 10 nucleotides to identify genes required for survival and growth in vitro and during infection of murine macrophages. The analysis revealed genomic regions important for fitness under two in vitro growth conditions. Overall, 724 coding genes were required for optimal growth in LB medium, and 851 coding genes were required for growth in SPI-2-inducing minimal medium. These findings were consistent with the essentiality analyses of other S. Typhimurium ST19 and S. Typhi strains. The global mutagenesis approach also identified 60 sRNAs and 413 intergenic regions required for growth in at least one in vitro growth condition. By infecting murine macrophages with the transposon library, we identified 68 genes that were required for intra-macrophage replication but did not impact fitness in vitro. None of these genes were unique to S. Typhimurium D23580, consistent with a high conservation of gene function between S. Typhimurium ST313 and ST19 and suggesting that novel virulence factors are not involved in the interaction of strain D23580 with murine macrophages. We discovered that transposon insertions rarely occurred in many pBT1 plasmid-encoded genes (36), compared with genes carried by the pSLT-BT virulence plasmid and other bacterial plasmids. The key essential protein encoded by pBT1 is a cysteinyl-tRNA synthetase, and our enzymological analysis revealed that the plasmid-encoded CysRSpBT1 had a lower ability to charge tRNA than the chromosomally-encoded CysRSchr enzyme. The presence of aminoacyl-tRNA synthetases in plasmids from a range of Gram-negative and Gram-positive bacteria suggests that plasmid-encoded essential genes are more common than had been appreciated.

ePath: an online database towards comprehensive essential gene annotation for prokaryotes

Experimental techniques for identification of essential genes (EGs) in prokaryotes are usually expensive, time-consuming and sometimes unrealistic. Emerging in silico methods provide alternative methods for EG prediction, but often possess limitations including heavy computational requirements and lack of biological explanation. Here we propose a new computational algorithm for EG prediction in prokaryotes with an online database (ePath) for quick access to the EG prediction results of over 4,000 prokaryotes ( https://www.pubapps.vcu.edu/epath/ ). In ePath, gene essentiality is linked to biological functions annotated by KEGG Ortholog (KO). Two new scoring systems, namely, E_score and P_score, are proposed for each KO as the EG evaluation criteria. E_score represents appearance and essentiality of a given KO in existing experimental results of gene essentiality, while P_score denotes gene essentiality based on the principle that a gene is essential if it plays a role in genetic information processing, cell envelope maintenance or energy production. The new EG prediction algorithm shows prediction accuracy ranging from 75% to 91% based on validation from five new experimental studies on EG identification. Our overall goal with ePath is to provide a comprehensive and reliable reference for gene essentiality annotation, facilitating the study of those prokaryotes without experimentally derived gene essentiality information.

Functional analysis of Salmonella Typhi adaptation to survival in water

Contaminated water is a major risk factor associated with the transmission of Salmonella enterica serovar Typhi (S. Typhi), the aetiological agent of human typhoid. However, little is known about how this pathogen adapts to living in the aqueous environment. We used transcriptome analysis (RNA‐seq) and transposon mutagenesis (TraDIS) to characterize these adaptive changes and identify multiple genes that contribute to survival. Over half of the genes in the S. Typhi genome altered expression level within the first 24 h following transfer from broth culture to water, although relatively few did so in the first 30 min. Genes linked to central metabolism, stress associated with arrested proton motive force and respiratory chain factors changed expression levels. Additionally, motility and chemotaxis genes increased expression, consistent with a scavenging lifestyle. The viaB‐associated gene tviC encoding a glcNAc epimerase that is required for Vi polysaccharide biosynthesis was, along with several other genes, shown to contribute to survival in water. Thus, we define regulatory adaptation operating in S. Typhi that facilitates survival in water.

Pyrimidine biosynthesis in pathogens - Structures and analysis of dihydroorotases from Yersinia pestis and Vibrio cholerae

The de novo pyrimidine biosynthesis pathway is essential for the proliferation of many pathogens. One of the pathway enzymes, dihydroorotase (DHO), catalyzes the reversible interconversion of N-carbamoyl-l-aspartate to 4,5-dihydroorotate. The substantial difference between bacterial and mammalian DHOs makes it a promising drug target for disrupting bacterial growth and thus an important candidate to evaluate as a response to antimicrobial resistance on a molecular level. Here, we present two novel three-dimensional structures of DHOs from Yersinia pestis (YpDHO), the plague-causing pathogen, and Vibrio cholerae (VcDHO), the causative agent of cholera. The evaluations of these two structures led to an analysis of all available DHO structures and their classification into known DHO types. Comparison of all the DHO active sites containing ligands that are listed in DrugBank was facilitated by a new interactive, structure-comparison and presentation platform. In addition, we examined the genetic context of characterized DHOs, which revealed characteristic patterns for different types of DHOs. We also generated a homology model for DHO from Plasmodium falciparum.

Identification of common highly expressed genes of Salmonella Enteritidis by in silico prediction of gene expression and in vitro transcriptomic analysis

Chickens are the reservoir host of Salmonella Enteritidis. Salmonella Enteritidis colonizes the gastro-intestinal tract of chickens and replicates within macrophages without causing clinically discernable illness. Persistence of S. Enteritidis in the hostile environments of intestinal tract and macrophages allows it to disseminate extra-intestinally to liver, spleen, and reproductive tract. Extra-intestinal dissemination into reproductive tract leads to contamination of internal contents of eggs, which is a major risk factor for human infection. Understanding the genes that contribute to S. Enteritidis persistence in the chicken host is central to elucidate the genetic basis of the unique pathobiology of this public health pathogen. The aim of this study was to identify a succinct set of genes associated with infection-relevant in vitro environments to provide a rational foundation for subsequent biologically-relevant research. We used in silico prediction of gene expression and RNA-seq technology to identify a core set of 73 S. Enteritidis genes that are consistently highly expressed in multiple S. Enteritidis strains cultured at avian physiologic temperature under conditions that represent intestinal and intracellular environments. These common highly expressed (CHX) genes encode proteins involved in bacterial metabolism, protein synthesis, cell-envelope biogenesis, stress response, and a few proteins with uncharacterized functions. Further studies are needed to dissect the contribution of these CHX genes to the pathobiology of S. Enteritidis in the avian host. Several of the CHX genes could serve as promising targets for studies towards the development of immunoprophylactic and novel therapeutic strategies to prevent colonization of chickens and their environment with S. Enteritidis.

Ribosomal maturation factor (RimP) is essential for survival of nontuberculous mycobacteria Mycobacterium fortuitum under in vitro acidic stress conditions

Mycobacterium fortuitum is an important human pathogenic NTM, which resists stress conditions inside macrophages by exploitation of specific genes. TnphoA-based transposon mutagenesis was employed to identify membrane genes responsible for survival of M. fortuitum under such stress conditions. A library of about 450 mutants was constructed after electroporation of vector pRT291 into wild-type M. fortuitum. On the basis of blue color development and alkaline phosphatase assay, 20 mutants were shortlisted to screen for growth and survival under acidic stress at pH 6.5, 5.5, 4.5, and 3.5. Mutant MT727 showed reduced growth and survival under acidic stress. The acid susceptible mutant MT727 was subjected to other in vitro stress conditions prevalent inside macrophages including oxidative, nutrient starvation and nitrosative stress. However, the mutant showed no appreciable difference in growth behavior under oxidative, nutrient starvation and nitrosative stress conditions in comparison to the wild type. Genomic and bioinformatics analysis of MT727 led to identification of putative ribosomal maturation factor RimP of M. fortuitum to be affected by mutagenesis, showing closest homology to M. abscessus RimP. In silico functional interaction of RimP protein using STRING database showed its interaction with proteins of ribosomal assembly and maturation. Results indicate role of rimP gene in survival of M. fortuitum under acidic stress conditions which may be further explored for use as a potential drug target against M. fortuitum and other mycobacterial infections.

Global Analysis of Genes Essential for Francisella tularensis Schu S4 Growth In Vitro and for Fitness during Competitive Infection of Fischer 344 Rats

The highly virulent intracellular pathogen Francisella tularensis is a Gram-negative bacterium that has a wide host range, including humans, and is the causative agent of tularemia. To identify new therapeutic drug targets and vaccine candidates and investigate the genetic basis of Francisella virulence in the Fischer 344 rat, we have constructed an F. tularensis Schu S4 transposon library. This library consists of more than 300,000 unique transposon mutants and represents a transposon insertion for every 6 bp of the genome. A transposon-directed insertion site sequencing (TraDIS) approach was used to identify 453 genes essential for growth in vitro Many of these essential genes were mapped to key metabolic pathways, including glycolysis/gluconeogenesis, peptidoglycan synthesis, fatty acid biosynthesis, and the tricarboxylic acid (TCA) cycle. Additionally, 163 genes were identified as required for fitness during colonization of the Fischer 344 rat spleen. This in vivo selection screen was validated through the generation of marked deletion mutants that were individually assessed within a competitive index study against the wild-type F. tularensis Schu S4 strain.IMPORTANCE The intracellular bacterial pathogen Francisella tularensis causes a disease in humans characterized by the rapid onset of nonspecific symptoms such as swollen lymph glands, fever, and headaches. F. tularensis is one of the most infectious bacteria known and following pulmonary exposure can have a mortality rate exceeding 50% if left untreated. The low infectious dose of this organism and concerns surrounding its potential as a biological weapon have heightened the need for effective and safe therapies. To expand the repertoire of targets for therapeutic development, we initiated a genome-wide analysis. This study has identified genes that are important for F. tularensis under in vitro and in vivo conditions, providing candidates that can be evaluated for vaccine or antibacterial development.

Transposon Insertion Site Sequencing in a Urinary Tract Model

Transposon sequencing (Tn-seq) is a technique that combines quantitative next-generation sequencing and a saturating transposon mutant library for an organism of interest, and ultimately allows for quantitation of the relative abundance of all of the mutants under a given condition, such as during experimental infection. The massively parallel sequencing capabilities of this technique provide a significant advance over more traditional methods of screening transposon mutant pools or individually determining the fitness contribution of genes of interest. Here, we describe a method for generating a genome-saturating transposon mutant library in Proteus mirabilis, determining the appropriate number of mutants for inoculation in an experimental infection model, preparing transposon insertion junctions for Illumina sequencing, and downstream analysis of mapped DNA sequencing reads for estimation of the contribution of each gene in the genome to fitness during infection.

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.

Metabolic models and gene essentiality data reveal essential and conserved metabolism in prokaryotes

Essential metabolic reactions are shaping constituents of metabolic networks, enabling viable and distinct phenotypes across diverse life forms. Here we analyse and compare modelling predictions of essential metabolic functions with experimental data and thereby identify core metabolic pathways in prokaryotes. Simulations of 15 manually curated genome-scale metabolic models were integrated with 36 large-scale gene essentiality datasets encompassing a wide variety of species of bacteria and archaea. Conservation of metabolic genes was estimated by analysing 79 representative genomes from all the branches of the prokaryotic tree of life. We find that essentiality patterns reflect phylogenetic relations both for modelling and experimental data, which correlate highly at the pathway level. Genes that are essential for several species tend to be highly conserved as opposed to non-essential genes which may be conserved or not. The tRNA-charging module is highlighted as ancestral and with high centrality in the networks, followed closely by cofactor metabolism, pointing to an early information processing system supplied by organic cofactors. The results, which point to model improvements and also indicate faults in the experimental data, should be relevant to the study of centrality in metabolic networks and ancient metabolism but also to metabolic engineering with prokaryotes.

If we tried to list every known chemical reaction within an organism–human, plant or even bacteria–we would get quite a long and confusing read. But when this information is represented in so-called genome-scale metabolic networks, we have the means to access computationally each of those reactions and their interconnections. Some parts of the network have alternatives, while others are unique and therefore can be essential for growth. Here, we simulate growth and compare essential reactions and genes for the simplest type of unicellular species–prokaryotes–to understand which parts of their metabolism are universally essential and potentially ancestral. We show that similar patterns of essential reactions echo phylogenetic relationships (this makes sense, as the genome provides the building plan for the enzymes that perform those reactions). Our computational predictions correlate strongly with experimental essentiality data. Finally, we show that a crucial step of protein synthesis (tRNA charging) and the synthesis and transformation of small molecules that enzymes require (cofactors) are the most essential and conserved parts of metabolism in prokaryotes. Our results are a step further in understanding the biology and evolution of prokaryotes but can also be relevant in applied studies including metabolic engineering and antibiotic design.

Iron-dependent essential genes in Salmonella Typhimurium

Background

The molecular mechanisms underlying bacterial cell death due to stresses or bactericidal antibiotics are complex and remain puzzling. Due to the current crisis of antibiotic resistance, development of effective antibiotics is urgently required. Previously, it has been shown that iron is required for effective killing of bacterial cells by numerous bactericidal antibiotics.

Results

We investigated the death or growth inhibition of S. Typhimurium under iron-restricted conditions, following disruption of essential genes, by transposon mutagenesis using transposon sequencing (Tn-seq). Our high-resolution Tn-seq analysis revealed that transposon mutants of S. Typhimurium with insertions in essential genes escaped immediate killing or growth inhibition under iron-restricted conditions for approximately one-third of all previously known essential genes. Based on this result, we classified all essential genes into two categories, iron-dependent essential genes, for which the insertion mutants can grow slowly if iron is restricted, and iron-independent essential genes, for which the mutants become nonviable regardless of iron concentration. The iron-dependency of the iron-dependent essential genes was further validated by the fact that the relative abundance of these essential gene mutants increased further with more severe iron restrictions. Our unexpected observation can be explained well by the common killing mechanisms of bactericidal antibiotics via production of reactive oxygen species (ROS). In this model, iron restriction would inhibit production of ROS, leading to reduced killing activity following blocking of essential gene functions. Interestingly, the targets of most antibiotics currently in use clinically are iron-dependent essential genes.

Conclusions

Our result suggests that targeting iron-independent essential genes may be a better strategy for future antibiotic development, because blocking their essential gene functions would lead to immediate cell death regardless of the iron concentration. This work expands our knowledge on the role of iron to a broad range of essential functions and pathways, providing novel insights for development of more effective antibiotics.

Electronic supplementary material

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

Pooled CRISPR interference screening enables genome-scale functional genomics study in bacteria with superior performance

To fully exploit the microbial genome resources, a high-throughput experimental platform is needed to associate genes with phenotypes at the genome level. We present here a novel method that enables investigation of the cellular consequences of repressing individual transcripts based on the CRISPR interference (CRISPRi) pooled screening in bacteria. We identify rules for guide RNA library design to handle the unique structure of prokaryotic genomes by tiling screening and construct an E. coli genome-scale guide RNA library (~60,000 members) accordingly. We show that CRISPRi outperforms transposon sequencing, the benchmark method in the microbial functional genomics field, when similar library sizes are used or gene length is short. This tool is also effective for mapping phenotypes to non-coding RNAs (ncRNAs), as elucidated by a comprehensive tRNA-fitness map constructed here. Our results establish CRISPRi pooled screening as a powerful tool for mapping complex prokaryotic genetic networks in a precise and high-throughput manner.

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.

TnSeq of Mycobacterium tuberculosis clinical isolates reveals strain-specific antibiotic liabilities

Once considered a phenotypically monomorphic bacterium, there is a growing body of work demonstrating heterogeneity among Mycobacterium tuberculosis (Mtb) strains in clinically relevant characteristics, including virulence and response to antibiotics. However, the genetic and molecular basis for most phenotypic differences among Mtb strains remains unknown. To investigate the basis of strain variation in Mtb, we performed genome-wide transposon mutagenesis coupled with next-generation sequencing (TnSeq) for a panel of Mtb clinical isolates and the reference strain H37Rv to compare genetic requirements for in vitro growth across these strains. We developed an analytic approach to identify quantitative differences in genetic requirements between these genetically diverse strains, which vary in genomic structure and gene content. Using this methodology, we found differences between strains in their requirements for genes involved in fundamental cellular processes, including redox homeostasis and central carbon metabolism. Among the genes with differential requirements were katG, which encodes the activator of the first-line antitubercular agent isoniazid, and glcB, which encodes malate synthase, the target of a novel small-molecule inhibitor. Differences among strains in their requirement for katG and glcB predicted differences in their response to these antimicrobial agents. Importantly, these strain-specific differences in antibiotic response could not be predicted by genetic variants identified through whole genome sequencing or by gene expression analysis. Our results provide novel insight into the basis of variation among Mtb strains and demonstrate that TnSeq is a scalable method to predict clinically important phenotypic differences among Mtb strains.

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

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

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'.

Role of sapA and yfgA in Susceptibility to Antibody-Mediated Complement-Dependent Killing and Virulence of Salmonella enterica Serovar Typhimurium

The ST313 pathovar of Salmonella enterica serovar Typhimurium contributes to a high burden of invasive disease among African infants and HIV-infected adults. It is characterized by genome degradation (loss of coding capacity) and has increased resistance to antibody-dependent complement-mediated killing compared with enterocolitis-causing strains of S. Typhimurium. Vaccination is an attractive disease-prevention strategy, and leading candidates focus on the induction of bactericidal antibodies. Antibody-resistant strains arising through further gene deletion could compromise such a strategy. Exposing a saturating transposon insertion mutant library of S. Typhimurium to immune serum identified a repertoire of S. Typhimurium genes that, when interrupted, result in increased resistance to serum killing. These genes included several involved in bacterial envelope biogenesis, protein translocation, and metabolism. We generated defined mutant derivatives using S. Typhimurium SL1344 as the host. Based on their initial levels of enhanced resistance to killing, yfgA and sapA mutants were selected for further characterization. The S. Typhimurium yfgA mutant lost the characteristic Salmonella rod-shaped appearance, exhibited increased sensitivity to osmotic and detergent stress, lacked very long lipopolysaccharide, was unable to invade enterocytes, and demonstrated decreased ability to infect mice. In contrast, the S. Typhimurium sapA mutants had similar sensitivity to osmotic and detergent stress and lipopolysaccharide profile and an increased ability to infect enterocytes compared with the wild type, but it had no increased ability to cause in vivo infection. These findings indicate that increased resistance to antibody-dependent complement-mediated killing secondary to genetic deletion is not necessarily accompanied by increased virulence and suggest the presence of different mechanisms of antibody resistance.

Functional Profiling of a Plasmodium Genome Reveals an Abundance of Essential Genes

The genomes of malaria parasites contain many genes of unknown function. To assist drug development through the identification of essential genes and pathways, we have measured competitive growth rates in mice of 2,578 barcoded Plasmodium berghei knockout mutants, representing >50% of the genome, and created a phenotype database. At a single stage of its complex life cycle, P. berghei requires two-thirds of genes for optimal growth, the highest proportion reported from any organism and a probable consequence of functional optimization necessitated by genomic reductions during the evolution of parasitism. In contrast, extreme functional redundancy has evolved among expanded gene families operating at the parasite-host interface. The level of genetic redundancy in a single-celled organism may thus reflect the degree of environmental variation it experiences. In the case of Plasmodium parasites, this helps rationalize both the relative successes of drugs and the greater difficulty of making an effective vaccine.

•

Two-thirds of Plasmodium berghei genes contribute to normal blood stage growth

•

The core genome of malaria parasites is highly optimized for rapid host colonization

•

Essential parasite genes and pathways are identified for drug target prioritization

•

Low functional redundancy reflects the constant environment encountered by a parasite

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.

Small RNAs in Bacterial Virulence and Communication

Bacterial pathogens must endure or adapt to different environments and stresses during transmission and infection. Posttranscriptional gene expression control by regulatory RNAs, such as small RNAs and riboswitches, is now considered central to adaptation in many bacteria, including pathogens. The study of RNA-based regulation (riboregulation) in pathogenic species has provided novel insight into how these bacteria regulate virulence gene expression. It has also uncovered diverse mechanisms by which bacterial small RNAs, in general, globally control gene expression. Riboregulators as well as their targets may also prove to be alternative targets or provide new strategies for antimicrobials. In this article, we present an overview of the general mechanisms that bacteria use to regulate with RNA, focusing on examples from pathogens. In addition, we also briefly review how deep sequencing approaches have aided in opening new perspectives in small RNA identification and the study of their functions. Finally, we discuss examples of riboregulators in two model pathogens that control virulence factor expression or survival-associated phenotypes, such as stress tolerance, biofilm formation, or cell-cell communication, to illustrate how riboregulation factors into regulatory networks in bacterial pathogens.

The Siderophore Product Ornibactin Is Required for the Bactericidal Activity of Burkholderia contaminans MS14

Burkholderia contaminans MS14 was isolated from soil in Mississippi. When it is cultivated on nutrient broth-yeast extract agar, the colonies exhibit bactericidal activity against a wide range of plant-pathogenic bacteria. A bacteriostatic compound with siderophore activity was successfully purified and was determined by nuclear magnetic resonance spectroscopy to be ornibactin. Isolation of the bactericidal compound has not yet been achieved; therefore, the exact nature of the bactericidal compound is still unknown. During an attempt to isolate the bactericidal compound, an interesting relationship between the production of ornibactin and the bactericidal activity of MS14 was characterized. Transposon mutagenesis resulted in two strains that lost bactericidal activity, with insertional mutations in a nonribosomal peptide synthetase (NRPS) gene for ornibactin biosynthesis and a luxR family transcriptional regulatory gene. Coculture of these two mutant strains resulted in restoration of the bactericidal activity. Furthermore, the addition of ornibactin to the NRPS mutant restored the bactericidal phenotype. It has been demonstrated that, in MS14, ornibactin has an alternative function, aside from iron sequestration. Comparison of the ornibactin biosynthesis genes in Burkholderia species shows diversity among the regulatory elements, while the gene products for ornibactin synthesis are conserved. This is an interesting observation, given that ornibactin is thought to have the same defined function within Burkholderia species. Ornibactin is produced by most Burkholderia species, and its role in regulating the production of secondary metabolites should be investigated.IMPORTANCE Identification of the antibacterial product from strain MS14 is not the key feature of this study. We present a series of experiments that demonstrate that ornibactin is directly involved in the bactericidal phenotype of MS14. This observation provides evidence for an alternative function for ornibactin, aside from iron sequestration. Ornibactin should be further evaluated for its role in regulating the biosynthesis of secondary metabolites in other Burkholderia species.

Molecular mechanism of mRNA repression in trans by a ProQ-dependent small RNA

Research into post-transcriptional control of mRNAs by small noncoding RNAs (sRNAs) in the model bacteria Escherichia coli and Salmonella enterica has mainly focused on sRNAs that associate with the RNA chaperone Hfq. However, the recent discovery of the protein ProQ as a common binding partner that stabilizes a distinct large class of structured sRNAs suggests that additional RNA regulons exist in these organisms. The cellular functions and molecular mechanisms of these new ProQ-dependent sRNAs are largely unknown. Here, we report in Salmonella Typhimurium the mode-of-action of RaiZ, a ProQ-dependent sRNA that is made from the 3' end of the mRNA encoding ribosome-inactivating protein RaiA. We show that RaiZ is a base-pairing sRNA that represses in trans the mRNA of histone-like protein HU-α. RaiZ forms an RNA duplex with the ribosome-binding site of hupA mRNA, facilitated by ProQ, to prevent 30S ribosome loading and protein synthesis of HU-α. Similarities and differences between ProQ- and Hfq-mediated regulation will be discussed.

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.