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

Long-chain unsaturated fatty acids sensor controlling the type III/VI secretion system is essential for Edwardsiella piscicida infection

Edwardsiella piscicida is an acute marine pathogen that causes severe damage to the aquaculture industry worldwide. The pathogenesis of E. piscicida is dependent mainly on the type III secretion system (T3SS) and type VI secretion system (T6SS), both of which are critically regulated by EsrB and EsrC. In this study, we revealed that fatty acids influence T3SS expression. Unsaturated fatty acids (UFAs), but not saturated fatty acids (SFAs), directly interact with EsrC, which abolishes the function of EsrC and results in the turn-off of T3/T6SS. Moreover, during the in vivo colonization of E. piscicida, host fatty acids were observed to be transported into E. piscicida through FadL and to modulate the expression of T3/T6SS. Furthermore, the esrCR38G mutant blocked the interaction between EsrC and UFAs, leading to dramatic growth defects in DMEM and impaired colonization in HeLa cells and zebrafish. In conclusion, this study revealed that the interaction between UFAs and EsrC to turn off T3/T6SS expression is essential for E. piscicida infection.

Fur-mediated regulation of hydrogen sulfide synthesis, stress response, and virulence in Edwardsiella piscicida

The production of endogenous hydrogen sulfide (H2S) is an important phenotype of bacteria. H2S plays an important role in bacterial resistance to ROS and antibiotics, which significantly contributes to bacterial pathogenicity. Edwardsiella piscicida, the Gram-negative pathogen causing fish edwardsiellosis, has been documented to produce hydrogen sulfide. In the study, we revealed that Ferric uptake regulator (Fur) controlled H2S synthesis by activating the expression of phsABC operon. Besides, Fur participated in the bacterial defense against ROS and cationic antimicrobial peptides and modulated T3SS expression. Furthermore, the disruption of fur exhibited a significant in vivo colonization defect. Collectively, our study demonstrated the regulation of Fur in H2S synthesis, stress response, and virulence, providing a new perspective for better understanding the pathogenesis of Edwardsiella.

Complete genome sequence and genome-wide transposon mutagenesis enable the determination of genes required for sodium hypochlorite tolerance and drug resistance in pathogen Aeromonas veronii GD2019

Aeromonas veronii, a significant pathogen in aquatic environments, poses a substantial threat to both human and animal health, particularly in aquaculture. In this study, we isolated A. veronii strain GD2019 from diseased largemouth bass (Micropterus salmoides) during a severe outbreak of aeromonad septicemia in Guangdong Province, China. The complete genome sequence of A. veronii GD2019 revealed that GD2019 contains a single chromosome of 4703,168 bp with an average G+C content of 58.3%. Phylogenetic analyses indicated that GD2019 forms a separate sub-branch in A. veronii and comparative genomic analyses identified the existence of an intact Type III secretion system. Moreover, to investigate the genes that are required for the conditional fitness of A. veronii under various stresses, a high-density transposon insertion library in GD2019 was generated by a Tn5-based transposon and covers 6311 genomic loci including 4155 genes and 2156 intergenic regions. Leveraging this library, 630 genes were classified as essential genes for growth in rich-nutrient LB medium. Furthermore, the genes GE001863/NtrC and GE002550 were found to confer tolerance to sodium hypochlorite in A. veronii. GE002562 and GE002614 were associated with the resistance to carbenicillin. Collectively, our results provide abundant genetic information on A. veronii, shedding light on the pathogenetic mechanisms of Aeromonas.

Lysine acetylation regulates the AT-rich DNA possession ability of H-NS

H-NS, the histone-like nucleoid-structuring protein in bacteria, regulates the stability of the bacterial genome by inhibiting the transcription of horizontally transferred genes, such as the type III and type VI secretion systems (T3/T6SS). While eukaryotic histone posttranslational modifications (PTMs) have been extensively studied, little is known about prokaryotic H-NS PTMs. Here, we report that the acetylation of H-NS attenuates its ability to silence horizontally transferred genes in response to amino acid nutrition and immune metabolites. Moreover, LC-MS/MS profiling showed that the acetyllysine sites of H-NS and K120 are indispensable for its DNA-binding ability. Acetylation of K120 leads to a low binding affinity for DNA and enhances T3/T6SS expression. Furthermore, acetylation of K120 impairs the AT-rich DNA recognition ability of H-NS. In addition, lysine acetylation in H-NS modulates in vivo bacterial virulence. These findings reveal the mechanism underlying H-NS PTMs and propose a novel mechanism by which bacteria counteract the xenogeneic silencing of H-NS.

Type III secretion system effector YfiD inhibits the activation of host poly(ADP-ribose) polymerase-1 to promote bacterial infection

Modulation of cell death is a powerful strategy employed by pathogenic bacteria to evade host immune clearance and occupy profitable replication niches during infection. Intracellular pathogens employ the type III secretion system (T3SS) to deliver effectors, which interfere with regulated cell death pathways to evade immune defenses. Here, we reveal that poly(ADP-ribose) polymerase-1 (PARP1)-dependent cell death restrains Edwardsiella piscicida's proliferation in mouse monocyte macrophages J774A.1, of which PARP1 activation results in the accumulation of poly(ADP-ribose) (PAR) and enhanced inflammatory response. Moreover, E. piscicida, an important intracellular pathogen, leverages a T3SS effector YfiD to impair PARP1's activity and inhibit PAR accumulation. Once translocated into the host nucleus, YfiD binds to the ADP-ribosyl transferase (ART) domain of PARP1 to suppress its PARylation ability as the pharmacological inhibitor of PARP1 behaves. Furthermore, the interaction between YfiD and ART mainly relies on the complete unfolding of the helical domain, which releases the inhibitory effect on ART. In addition, YfiD impairs the inflammatory response and cell death in macrophages and promotes in vivo colonization and virulence of E. piscicida. Collectively, our results establish the functional mechanism of YfiD as a potential PARP1 inhibitor and provide more insights into host defense against bacterial infection.

Analysis of key genes for the survival of Pantoea agglomerans under nutritional stress

The absolute amount of nutrients on plant leaves is usually low, and the growth of epiphytic bacteria is typically limited by nutrient content. Thus, is of great significance to study the survival mechanism of epiphytes under nutritional stress for plant disease control. In this paper, Pantoea agglomerans CHTF15 isolated from walnut leaves was used to detect the key genes for the survival of the bacterium under simulated nutrient stress in artificial medium. Genome sequencing was combined with transposon insertion sequencing (Tn-seq) for the detection technique. A total of 105 essential genes were screened from the whole genome. The genes were mainly related to the nucleotide metabolism, protein metabolism, biological oxidation and the gene repair of bacteria analyzed by gene ontology (GO) enrichment analysis. Volcano map analysis demonstrated that the functions of the 15 genes with the most significant differences were generally related to the synthesis of amino acids or proteins, the nucleotide metabolism and homologous recombination and repair. Competitive index analysis revealed that the deletion of the genes dksA and epmA regulating protein synthesis, the gene ribB involved in the nucleotide metabolism and the gene xerD involved in recombination repair induced a significant reduction in the survival ability of the corresponding mutants in the 0.10 % YEP medium and the walnut leaf surface. The results act as a foundation for further in-depth research on the infection process and the mechanisms of pathogenic bacteria.

Identification of a novel transcriptional regulator, CorR, for copper stress response in Edwardsiella piscicida

Copper plays a vital role in the host-pathogen interface, potentially making components of the bacterial copper response suitable targets for the development of innovative antimicrobial strategies. The anti-copper arsenal of intracellular pathogens has expanded as an adaptation to survive copper toxicity in order to escape intracellular killing by the host immune system. Herein, we employed transposon insertion sequencing to investigate the genetic mechanisms underlying the survival of Edwardsiella piscicida under copper stress. A novel transcriptional regulator, ETAE_2324 (named CorR), was identified to participate in the response to copper ions by controlling the expression of copA, the core component of cytoplasmic copper homeostasis. Furthermore, CorR regulated the expression of virulent determinant eseB, influencing the in vivo colonization of E. piscicida. Collectively, our results contribute to the comprehension of the underlying mechanism of the adaption of intracellular pathogens to copper stress during bacterial infections.IMPORTANCECopper ions play a pivotal role in the interaction between bacteria and the host during infection. The host's innate immune system employs copper ions for their bactericidal properties, thereby making bacterial copper tolerance a crucial determinant of virulence. Edwardsiella piscicida, a significant marine pathogen, has caused substantial losses in the global aquaculture industry. To comprehensively investigate how E. piscicida responds to copper stress, we utilized transposon insertion sequencing to explore genes associated with copper tolerance in culture media containing different concentrations of copper ions. A novel transcriptional regulator, CorR, was identified to respond to copper ions and regulates the expression of crucial components of copper homeostasis CopA, along with the essential virulence factor EseB. These findings offer valuable insights into the underlying mechanisms that govern bacterial copper tolerance and present novel perspectives for the development of vaccines and therapeutic strategies targeting E. piscicida.

Versatile lifestyles of Edwardsiella: Free-living, pathogen, and core bacterium of the aquatic resistome

Edwardsiella species in aquatic environments exist either as individual planktonic cells or in communal biofilms. These organisms encounter multiple stresses, include changes in salinity, pH, temperature, and nutrients. Pathogenic species such as E. piscicida, can multiply within the fish hosts. Additionally, Edwardsiella species (E. tarda), can carry antibiotic resistance genes (ARGs) on chromosomes and/or plasmids, that can be transmitted to the microbiome via horizontal gene transfer. E. tarda serves as a core in the aquatic resistome. Edwardsiela uses molecular switches (RpoS and EsrB) to control gene expression for survival in different environments. We speculate that free-living Edwardsiella can transition to host-living and vice versa, using similar molecular switches. Understanding such transitions can help us understand how other similar aquatic bacteria switch from free-living to become pathogens. This knowledge can be used to devise ways to slow down the spread of ARGs and prevent disease outbreaks in aquaculture and clinical settings.

Bridging the Gap Between Environmental Adversity and Neuropsychiatric Disorders: The Role of Transposable Elements

Long regarded as “junk DNA,” transposable elements (TEs) have recently garnered much attention for their role in promoting genetic diversity and plasticity. While many processes involved in mammalian development require TE activity, deleterious TE insertions are a hallmark of several psychiatric disorders. Moreover, stressful events including exposure to gestational infection and trauma, are major risk factors for developing psychiatric illnesses. Here, we will provide evidence demonstrating the intersection of stressful events, atypical TE expression, and their epigenetic regulation, which may explain how neuropsychiatric phenotypes manifest. In this way, TEs may be the “bridge” between environmental perturbations and psychopathology.

FabR senses long-chain unsaturated fatty acids to control virulence in pathogen Edwardsiella piscicida

Long-chain unsaturated fatty acids (UFAs) can serve as nutrient sources or building blocks for bacterial membranes. However, little is known about how UFAs may be incorporated into the virulence programs of pathogens. A previous investigation identified FabR as a positive regulator of virulence gene expression in Edwardsiella piscicida. Here, chromatin immunoprecipitation-sequencing coupled with RNA-seq analyses revealed that 10 genes were under the direct control of FabR, including fabA, fabB, and cfa, which modulate the composition of UFAs. The binding of FabR to its target DNA was facilitated by oleoyl-CoA and inhibited by stearoyl-CoA. In addition, analyses of enzyme mobility shift assay and DNase I footprinting with wild-type and a null mutant (F131A) of FabR demonstrated crucial roles of FabR in binding to the promoters of fabA, fabB, and cfa. Moreover, FabR also binds to the promoter region of the virulence regulator esrB for its activation, facilitating the expression of the type III secretion system (T3SS) in response to UFAs. Furthermore, FabR coordinated with RpoS to modulate the expression of T3SS. Collectively, our results elucidate the molecular machinery of FabR regulating bacterial fatty acid composition and virulence in enteric pathogens, further expanding our knowledge of its crucial role in host-pathogen interactions.

Long-read sequencing for identification of insertion sites in large transposon mutant libraries

Transposon insertion site sequencing (TIS) is a powerful method for associating genotype to phenotype. However, all TIS methods described to date use short nucleotide sequence reads which cannot uniquely determine the locations of transposon insertions within repeating genomic sequences where the repeat units are longer than the sequence read length. To overcome this limitation, we have developed a TIS method using Oxford Nanopore sequencing technology that generates and uses long nucleotide sequence reads; we have called this method LoRTIS (Long-Read Transposon Insertion-site Sequencing). LoRTIS enabled the unique localisation of transposon insertion sites within long repetitive genetic elements of E. coli, such as the transposase genes of insertion sequences and copies of the ~ 5 kb ribosomal RNA operon. We demonstrate that LoRTIS is reproducible, gives comparable results to short-read TIS methods for essential genes, and better resolution around repeat elements. The Oxford Nanopore sequencing device that we used is cost-effective, small and easily portable. Thus, LoRTIS is an efficient means of uniquely identifying transposon insertion sites within long repetitive genetic elements and can be easily transported to, and used in, laboratories that lack access to expensive DNA sequencing facilities.

Dysregulation of Cytosolic c-di-GMP in Edwardsiella piscicida Promotes Cellular Non-Canonical Ferroptosis

Programmed cell death plays an important role in modulating host immune defense and pathogen infection. Ferroptosis is a type of inflammatory cell death induced by intracellular iron-dependent accumulation of toxic lipid peroxides. Although ferroptosis has been associated with cancer and other sterile diseases, very little is known about the role of ferroptosis in modulating host-pathogen interactions. We show that accumulation of the secondary messenger bis-(3′,5′)-cyclic dimeric GMP (c-di-GMP) in the pathogenic bacterium Edwardsiella piscicida (E. piscicida) triggers a non-canonical ferroptosis pathway in infected HeLa cells. Moreover, we observed that the dysregulation of c-di-GMP in E. piscicida promotes iron accumulation, mitochondrial dysfunction, and production of reactive oxygen species, all of which that can be blocked by iron chelator. Importantly, unlike classical ferroptosis that is executed via excess lipid peroxidation, no lipid peroxidation was detected in the infected cells. Furthermore, lipoxygenases inhibitors and lipophilic antioxidants are not able to suppress morphological changes and cell death induced by E. piscicida mutant producing excess c-di-GMP, and this c-di-GMP dysregulation attenuates bacterial virulence in vivo. Collectively, our results reveal a novel non-canonical ferroptosis pathway mediated by bacterial c-di-GMP and provide evidence for a role of ferroptosis in the regulation of pathogen infection.

Interplay between ferric uptake regulator Fur and horizontally acquired virulence regulator EsrB coordinates virulence gene expression in Edwardsiella piscicida

Edwardsiella piscicida mediates hemorrhagic septicemia and is a leading pathogen of fish. E. piscicida invades and colonizes macrophages using type III and VI secretion systems (T3/T6SS) that are controlled by a two-component system (TCS) EsrA-EsrB. Iron acquisition is essential for E. piscicida pathogenesis and coordination between iron and TCS signaling in modulating bacterial virulence is not well understood. Here, we show that iron uptake systems are co-regulated by ferric uptake regulator (Fur) in E. piscicida. Fur bound to 98 genes that harbored conserved Fur-box to globally control the expression of ∼755 genes, including those encoding iron uptake systems, T3/T6SS, and Icc, cAMP phosphodiesterase that represses biofilm formation. Additionally, Fur, in complex with iron, bound to the esrB promoter to repress expression and ultimately attenuated virulence. Conversely, EsrB activated the expression of T3/T6SS and iron uptake systems to mitigate a shortage of intracellular iron during iron scarcity. Furthermore, EsrB directly bound to and activated the fur promoter, leading to Fur-ferrous ion-dependent esrB repression in the presence of iron. Finally, Fur-EsrB interplay was essential for bacterial fitness during in vivo infection and survival in seawater environments. Collectively, we highlight the mechanisms that underlie the reciprocal regulatory networks of iron homeostasis and virulence systems in E. piscicida.

Phosphate transport system mediates the resistance of Enterococcus faecalis to multidrug

Enterococcus faecalis, a severe nosocomial and community opportunistic pathogen, is difficult to control due to its multidrug resistance. Through heredity and the recombination of intrinsic resistance genes and horizontally acquired resistance genes, E. faecalis can rapidly evolve drug resistance. Nisin, an important antimicrobial peptide, is extensively employed in the healthcare and food industries to inhibit Gram-positive bacteria and may induce the emergence of nisin-resistant bacteria worldwide. However, the mechanism governing nisin resistance in E. faecalis has not been fully elucidated. This study utilizes transposon insertion sequencing (TIS) to comprehensively explore novel genes related to nisin resistance. According to the analysis of TIS results, hundreds of genes appear to be essential for nisin resistance in E. faecalis. The phosphate transport system (OG1RF_10018-10021, named PTS), which is screened by TIS results, enhances the resistance of E. faecalis to nisin, the mechanism of which may be involved in potA and/or OG1RF_10526 (hypothetical gene). Meanwhile, PTS also strongly represses the biosynthesis of ribosomes to increase the sensitivity of E. faecalis to gentamycin. In addition, the overexpression of PTS increases the sensitivity of E. faecalis to daptomycin, the mechanism of which is independent of the LiaFSR system. This study first demonstrated that E. faecalis utilizes PTS to mediate the resistance to multidrug, which may help to elucidate the mechanism governing drug resistance and to establish guidelines for the treatment of infectious diseases caused by E. faecalis.

Hijacking Transposable Elements for Saturation Mutagenesis in Fungi

Transposable elements are present in almost all known genomes, these endogenous transposons have recently been referred to as the mobilome. They are now increasingly used in research in order to make extensive mutant libraries in different organisms. Fungi are an essential part of our lives on earth, they influence the availability of our food and they live inside our own bodies both as commensals and pathogenic organisms. Only few fungal species have been studied extensively, mainly due to the lack of appropriate molecular genetic tools. The use of transposon insertion libraries can however help to rapidly advance our knowledge of (conditional) essential genes, compensatory mutations and drug target identification in fungi. Here we give an overview of some recent developments in the use of different transposons for saturation mutagenesis in different fungi.

Gene Silencing Through CRISPR Interference in Bacteria: Current Advances and Future Prospects

Functional genetic screening is an important method that has been widely used to explore the biological processes and functional annotation of genetic elements. CRISPR/Cas (Clustered regularly interspaced short palindromic repeat sequences/CRISPR-associated protein) is the newest tool in the geneticist's toolbox, allowing researchers to edit a genome with unprecedented ease, accuracy, and high-throughput. Most recently, CRISPR interference (CRISPRi) has been developed as an emerging technology that exploits the catalytically inactive Cas9 (dCas9) and single-guide RNA (sgRNA) to repress sequence-specific genes. In this review, we summarized the characteristics of the CRISPRi system, such as programmable, highly efficient, and specific. Moreover, we demonstrated its applications in functional genetic screening and highlighted its potential to dissect the underlying mechanism of pathogenesis. The recent development of the CRISPRi system will provide a high-throughput, practical, and efficient tool for the discovery of functionally important genes in bacteria.

Genome-Wide Identification of Genes Involved in Acid Stress Resistance of Salmonella Derby

Resistance to and survival under acidic conditions are critical for Salmonella to infect the host. As one of the most prevalent serotypes identified in pigs and humans, how S. Derby overcomes acid stress remains unclear. Here, we de novo sequenced the genome of a representative S. Derby strain 14T from our S. Derby strain stock and identified its acid resistance-associated genes using Tn-seq analysis. A total of 35 genes, including those belonging to two-component systems (TCS) (cpxAR), the CRISPR-Cas system (casCE), and other systems, were identified as essential for 14T to survive under acid stress. The results demonstrated that the growth curve and survival ability of ΔcpxA and ΔcpxR were decreased under acid stress, and the adhesion and invasion abilities to the mouse colon cancer epithelial cells (MC38) of ΔcpxR were also decreased compared with the wild type strain, suggesting that the TCS CpxAR plays an essential role in the acid resistance and virulence of S. Derby. Also, CasC and CasE were found to be responsible for acid resistance in S. Derby. Our results indicate that acid stress induces multiple genes’ expression to mediate the acid resistance of S. Derby and enhance its pathogenesis during an infection.

Identification and Characterization of EvpQ, a Novel T6SS Effector Encoded on a Mobile Genetic Element in Edwardsiella piscicida

In this study, a hypothetical protein (ORF02740) secreted by Edwardsiella piscicida was identified. We renamed the ORF02740 protein as EvpQ, which is encoded by a mobile genetic element (MGE) in E. piscicida genome. The evpQ gene is spaced by 513 genes from type VI secretion system (T6SS) gene cluster. Low GC content, three tRNA, and three transposase genes nearby evpQ define this MGE that evpQ localizes as a genomic island. Sequence analysis reveals that EvpQ shares a conserved domain of C70 family cysteine protease and shares 23.91% identity with T3SS effector AvrRpt2 of phytopathogenic Erwinia amylovora. Instead, EvpQ of E. piscicida is proved to be secreted at a T6SS-dependent manner, and it can be translocated into host cells. EvpQ is thereof a novel T6SS effector. Significantly decreased competitive index of ΔevpQ strain in blue gourami fish (0.53 ± 0.27 in head kidney and 0.44 ± 0.19 in spleen) indicates that EvpQ contributes to the pathogenesis of E. piscicida. At 8-, 18-, and 24-h post-subculture into DMEM, the transcription of evpQ was found to be negatively regulated by Fur and positively regulated by EsrC, and the steady-state protein levels of EvpQ are negatively controlled by RpoS. Our study lays a foundation for further understanding the pathogenic role of T6SS in edwardsiellosis.

Genome-Wide fitness analysis of group B Streptococcus in human amniotic fluid reveals a transcription factor that controls multiple virulence traits

Streptococcus agalactiae (group B Streptococcus; GBS) remains a dominant cause of serious neonatal infections. One aspect of GBS that renders it particularly virulent during the perinatal period is its ability to invade the chorioamniotic membranes and persist in amniotic fluid, which is nutritionally deplete and rich in fetal immunologic factors such as antimicrobial peptides. We used next-generation sequencing of transposon-genome junctions (Tn-seq) to identify five GBS genes that promote survival in the presence of human amniotic fluid. We confirmed our Tn-seq findings using a novel CRISPR inhibition (CRISPRi) gene expression knockdown system. This analysis showed that one gene, which encodes a GntR-class transcription factor that we named MrvR, conferred a significant fitness benefit to GBS in amniotic fluid. We generated an isogenic targeted deletion of the mrvR gene, which had a growth defect in amniotic fluid relative to the wild type parent strain. The mrvR deletion strain also showed a significant biofilm defect in vitro. Subsequent in vivo studies showed that while the mutant was able to cause persistent murine vaginal colonization, pregnant mice colonized with the mrvR deletion strain did not develop preterm labor despite consistent GBS invasion of the uterus and the fetoplacental units. In contrast, pregnant mice colonized with wild type GBS consistently deliver prematurely. In a sepsis model the mrvR deletion strain showed significantly decreased lethality. In order to better understand the mechanism by which this newly identified transcription factor controls GBS virulence, we performed RNA-seq on wild type and mrvR deletion GBS strains, which revealed that the transcription factor affects expression of a wide range of genes across the GBS chromosome. Nucleotide biosynthesis and salvage pathways were highly represented among the set of differentially expressed genes, suggesting that MrvR may be involved in regulating nucleotide availability.

Group B Streptococcus (GBS) is a species of Gram-positive bacteria that often colonizes the healthy adult intestinal and reproductive tracts without causing serious symptoms. During pregnancy, however, GBS can invade the pregnant uterus, where it can cause infection of the placenta, fetal membranes, and fetus—a condition known as chorioamnionitis. Chorioamnionitis is associated with serious adverse pregnancy outcomes, including stillbirth, preterm labor, and severe infection of the newborn. GBS can survive in human amniotic fluid, which is low in bacterial nutrients and contains immune molecules that limit microbial persistence, and this ability likely contributes to GBS chorioamnionitis. This study is focused on a single GBS gene that encodes a genetic regulator we called MrvR, which we show is important for GBS resistance to human amniotic fluid. Using a series of genetic techniques combined with animal models of GBS colonization and infection, we show that MrvR also plays a key role in allowing GBS to invade the bloodstream and trigger the inflammatory responses that lead to preterm labor and stillbirth. The study concludes with a survey of other GBS genes whose activity is regulated by MrvR, which seems to be an important contributor to GBS virulence.

Bacterial infection reinforces host metabolic flux from arginine to spermine for NLRP3 inflammasome evasion

Hosts recognize cytosolic microbial infection via the nucleotide-binding domain-like receptor (NLR) protein family, triggering inflammasome complex assembly to provoke pyroptosis or cytokine-related caspase-1-dependent antimicrobial responses. Pathogens have evolved diverse strategies to antagonize inflammasome activation. Here, Edwardsiella piscicida gene-defined transposon library screening for lactate dehydrogenase (LDH) release in nlrc4^-/- bone marrow-derived macrophages (BMDMs) demonstrates that genes clustered in the bacterial arginine metabolism pathway participate in NLRP3 inflammasome inhibition. Blocking arginine uptake or putrescine export significantly relieves NLRP3 inflammasome inhibition, indicating that this bacterium rewires its arginine metabolism network during infection. Moreover, intracellular E. piscicida recruits the host arginine importer (mCAT-1) and putrescine exporter (Oct-2) to bacterium-containing vacuoles, accompanied by reduced arginine and accumulated cytosolic spermine. Neutralizing E. piscicida-induced cytosolic spermine enhancement by spermine synthetase or extracellular spermine significantly alters NLRP3 inflammasome activation. Importantly, accumulated cytosolic spermine inhibits K⁺ efflux-dependent NLRP3 inflammasome activation. These data highlight the mechanism of bacterial gene-mediated arginine metabolism control for NLRP3 inflammasome evasion.

Characterization of a novel live attenuated Edwardsiella piscicida vaccine based on the overexpressed type III secretion system and systematic deletion of the associated effectors

Edwardsiella piscicida causes edwardsiellosis in a variety of fish species and leads to tremendous economic losses in the global aquaculture industries. Thus, effective and safe prevention and control of this bacterium are urgently needed to combat the related infections. Live attenuated vaccines (LAVs) effectively prevent infectious diseases. However, most of the existing E. piscicida LAVs are based on the deletion of genes encoding the translocon components of the type III secretion system (T3SS), the core virulence system, which is the most prominent protective bacterial antigen with the strongest immunogenicity. In this study, we systematically deleted all of the 9 established T3SS effectors in E. piscicida (aka 9Δ) and the rpoS gene encoding the alternative sigma factor, the esrB repressor (10Δ), then we overexpressed esrB and T3SS in E. piscicida to obtain the recombinant strain 10Δ/esrBOE. The modified strains 10Δ and 10Δ/esrBOE exhibited severe attenuation and in vivo colonization defects. Additionally, vaccination by intraperitoneal injection with 10Δ and 10Δ/esrBOE could significantly upregulate the expression of the antigen recognition related gene (TLR5) and the adaptive immune response-related gene (MHC II) in the spleen/kidney of turbot fish, and it also enhanced the hosts' serum bactericidal capacity. Finally, vaccination with 10Δ/esrBOE led to increased immune protection against the challenge of wild type E. piscicida EIB202 in turbot fish. Collectively, these findings demonstrated that 10Δ/esrBOE was a novel LAV strain and therefore a potential novel strategy for the construction of LAVs against bacterial pathogens.

MviN mediates the regulation of environmental osmotic pressure on esrB to control the virulence in the marine fish pathogen Edwardsiella piscicida

Edwardsiella piscicida is a notorious pathogen infecting diverse kinds of fish and causes substantial economic losses in the global aquaculture industries. The EsrA-EsrB two-component system plays a critical role in the regulation of virulence genes expression, including type III and type VI secretion systems (T3/T6SSs). In this study, the putative regulators of esrB were screened by the transposon insertion sequencing (TIS) technology. As a result, MviN, a lipid II flippase, was identified as a modulator to upregulate esrB and downstream T3/T6SS gene expression in the earlier growth phases while downregulate esrB at the later stages. Complement or overexpression of the mviN restored the esrB as well as T3/T6SS expression in the ΔmviN mutant strain. Moreover, MviN also mediated the regulation of environmental osmotic pressure on the expression of esrB. MviN was also found to significantly influence the in vivo colonization of E. piscicida in turbot. Collectively, this study enhanced our understanding of pathogenesis and virulence regulatory network of E. piscicida.

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

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

Application of Transposon Insertion Sequencing to Agricultural Science

Many plant-associated bacteria have the ability to positively affect plant growth and there is growing interest in utilizing such bacteria in agricultural settings to reduce reliance on pesticides and fertilizers. However, our capacity to utilize microbes in this way is currently limited due to patchy understanding of bacterial-plant interactions at a molecular level. Traditional methods of studying molecular interactions have sought to characterize the function of one gene at a time, but the slow pace of this work means the functions of the vast majority of bacterial genes remain unknown or poorly understood. New approaches to improve and speed up investigations into the functions of bacterial genes in agricultural systems will facilitate efforts to optimize microbial communities and develop microbe-based products. Techniques enabling high-throughput gene functional analysis, such as transposon insertion sequencing analyses, have great potential to be widely applied to determine key aspects of plant-bacterial interactions. Transposon insertion sequencing combines saturation transposon mutagenesis and high-throughput sequencing to simultaneously investigate the function of all the non-essential genes in a bacterial genome. This technique can be used for both in vitro and in vivo studies to identify genes involved in microbe-plant interactions, stress tolerance and pathogen virulence. The information provided by such investigations will rapidly accelerate the rate of bacterial gene functional determination and provide insights into the genes and pathways that underlie biotic interactions, metabolism, and survival of agriculturally relevant bacteria. This knowledge could be used to select the most appropriate plant growth promoting bacteria for a specific set of conditions, formulating crop inoculants, or developing crop protection products. This review provides an overview of transposon insertion sequencing, outlines how this approach has been applied to study plant-associated bacteria, and proposes new applications of these techniques for the benefit of agriculture.

Edwardsiella piscicida: A versatile emerging pathogen of fish

Edwardsiella piscicida is an Enterobacteriaceae that is abundant in water and causes food and waterborne infections in fish, animals, and humans. The bacterium causes Edwardsiellosis in farmed fish and can lead to severe economic losses in aquaculture worldwide. E. piscicida is an intracellular pathogen that can also cause systemic infection. Type III and type VI secretion systems are the bacterium’s most lethal weapons against host defenses. It also possesses multi-antibiotic resistant genes and is selected and enriched in the environment due to the overuse of antibiotics. Therefore, the bacterium has great potential to contribute to the evolution of the resistome. All these properties have made this bacterium a perfect model to study bacteria virulence mechanisms and the spread of antimicrobial genes in the environment. We summarize recent advance in E. piscicida biology and provide insights into future research in virulence mechanisms, vaccine development and novel therapeutics.

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

Background

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

Results

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

Conclusions

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

A Bacterial Pathogen Senses Host Mannose to Coordinate Virulence

Bacterial pathogens are thought to activate expression of virulence genes upon detection of host-associated cues, but identification of the nature of such signals has proved difficult. We generated a genome-scale defined transposon mutant library in Edwardsiella piscicida, an important fish pathogen, to quantify the fitness of insertion mutants for intracellular growth in macrophages and in turbot (Scophthalmus maximus). These screens identified EvrA, a transcription activator that induces expression of esrB, a key virulence regulator. EvrA is directly bound and activated by mannose-6-phosphate (man-6P) derived from actively imported mannose. Mutants lacking EvrA or expressing an EvrA unable to bind man-6P were similarly attenuated in turbot. Exogenously added mannose promoted E. piscicida virulence, and high levels of mannose were detected in fish tissue. Together, these observations reveal that binding of a host-derived sugar to a transcription factor can facilitate pathogen sensing of the host environment and trigger virulence programs.

EsrB negatively regulates expression of the glutamine sythetase GlnA in the fish pathogen Edwardsiella piscicida

Edwardsiella piscicida is a gram-negative bacterial pathogen invading a wide range of fish species. Response regulator EsrB is essential for the activation of type III and type VI secretion systems (T3/T6SS). In this study, proteomes of the wild-type E. piscicida EIB202 and ΔesrB mutant strains were compared to identify the regulon components of EsrB cultured in DMEM allowing T3/T6SS expression. As a result, 19 proteins showed different expression, which were identified to be associated with T3/T6SS, related to amino acid transport and metabolism, and energy production. Particularly, GlnA, a glutamine synthetase essential for ammonia assimilation and glutamine biosynthesis from glutamate, was found to be regulated negatively by EsrB. Moreover, GlnA affected bacterial growth in vitro and bacterial colonization in vivo. Collectively, our results indicated that EsrB plays important roles in regulating the expression of metabolic pathways and virulence genes, including glutamine biosynthesis in E. piscicida during infection.

Edwardsiella induces microtubule-severing in host epithelial cells

Edwardsiella bacteria cause economic losses to a variety of commercially important fish globally. Human infections are rare and result in a gastroenteritis-like illness. Because these bacteria are evolutionarily related to other Enterobacteriaceae and the host cytoskeleton is a common target of enterics, we hypothesized that Edwardsiella may cause similar phenotypes. Here we use HeLa and Caco-2 infection models to show that microtubules are severed during the late infections. This microtubule alteration phenotype was not dependant on the type III or type VI secretion system (T3SS and T6SS) of the bacteria as ΔT3SS and ΔT6SS mutants of E. piscicida EIB202 and E. tarda ATCC15947 that lacks both also caused microtubule disassembly. Immunolocalization experiments showed the host katanin catalytic subunits A1 and A like 1 proteins at regions of microtubule severing, suggesting their involvement in the microtubule disassembly events. To identify bacterial components involved in this phenotype, we screened a 2,758 transposon library of E. piscicida EIB202 and found that 4 single mutations in the atpFHAGDC operon disrupted microtubule disassembly in HeLa cells. We then constructed three atp deletion mutants; they all could not disassemble host microtubules. This work provides the first clear evidence of host cytoskeletal alterations during Edwardsiella infections.

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

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

Identification and characterization of FtsH mediating in vivo colonization and stress adaptation in the fish pathogen Edwardsiella piscicida

Edwardsiella piscicida is an important pathogenic enteric bacterium of fish. FtsH is a unique membrane-anchored AAA + protease that regulates protein homeostasis in bacteria. In cooperation with modulators HflK and HflC, FtsH is essential in enteric bacteria and controls the response to environmental stresses. Here, we used in vivo pattern analysis of conditional essentiality (PACE) and identified that ftsH and hflK/C were associated with impaired in vivo colonization in Edw. piscicida and attenuated internalization ability of ZF4 cells. The ftsH mutant displayed increased survival during prolonged treatment of starvation and high osmotic stresses in Edw. piscicida. Further analysis showed that the disruption of ftsH resulted in the overproduction of the established substrate LpxC, which is responsible for the synthesis of LPS (lipopolysaccharide), as well as the substrate YfgM, which is involved in high osmolality tolerance during stationary phase. However, the inconsistency in the abilities of the ftsH and hflK/C mutants to achieve YfgM-based osmotic resistance indicated that there might be multiple, while distinctive, pathways controlled by FtsH and the associated modulator proteins HflK/C. This investigation revealed the unique functions of FtsH and its modulator HflK/C in Edw. piscicida.

Pattern analysis of conditional essentiality (PACE)-based heuristic identification of an in vivo colonization determinant as a novel target for the construction of a live attenuated vaccine against Edwardsiella piscicida

Edwardsiella piscicida is the aetiological agent of fish edwardsiellosis, causing huge economic losses in aquaculture industries. The use of a live attenuated vaccine (LAV) will be an effective strategy to control the disease in farmed fish. Thus, methods facilitating exploration of targets used for construction of an LAV will be of great significance. Previously, we devised an algorithm termed pattern analysis of conditional essentiality (PACE) to perform genome-wide analysis of the temporal dynamic behaviour of E. piscicida mutants colonizing turbot. Here, we correlated the conditional essentiality patterns of the PACE-derived colonization determinants with that of the aroC gene encoding chorismate synthase, the established target for LAV construction in E. piscicida, and identified ETAE_0023 as a novel valuable LAV target. ETAE_0023 encodes an uncharacterized DcrB family protein. Deletion of ETAE_0023 dramatically impaired E. piscicida invasion capability in ZF4 cells as well as colonization in fish and resulted in in vivo clearance at ∼30 days post-infection. ΔETAE_0023 showed an ∼2500-fold higher 50% lethal dose (LD₅₀) than that of the wild type strain. Vaccination with ΔETAE_0023 by intraperitoneal (i.p.) injection upregulated expression of immune factors, i.e., IL-1β, IgM, MHC-I and MHC-II, and produced significantly high levels of E. piscicida-specific IgM as well as serum bactericidal capacities in turbot. Moreover, a single i.p. inoculation with ΔETAE_0023 generated significant protection comparable to the established WED LAV strain in turbot against challenge with the wild type strain after 5 weeks of vaccination. Taken together, we demonstrated a PACE-based method for heuristic identification of targets for LAV construction and presented ΔETAE_0023 as a new LAV candidate against edwardsiellosis.

Genome-Wide Identification of Fitness Factors in Seawater for Edwardsiella piscicida

Marine pathogens are transmitted from one host to another through seawater. Therefore, it is important for marine pathogens to maintain survival or growth in seawater. However, little is known about how marine pathogens adapt to living in seawater environments. Here, transposon insertion sequencing was performed to explore the genetic determinants of Edwardsiella piscicida survival in seawater at 16 and 28°C. Seventy-one mutants with mutations mainly in metabolism-, transportation-, and type III secretion system (T3SS)-related genes showed significantly increased or impaired fitness in 16°C water. In 28°C seawater, 63 genes associated with transcription and translation, as well as energy production and conversion, were essential for optimal survival of the bacterium. In particular, 11 T3SS-linked mutants displayed enhanced fitness in 16°C seawater but not in 28°C seawater. In addition, 13 genes associated with oxidative phosphorylation and 4 genes related to ubiquinone synthesis were identified for survival in 28°C seawater but not in 16°C seawater, which suggests that electron transmission and energy-producing aerobic respiration chain factors are indispensable for E. piscicida to maintain survival in higher-temperature seawater. In conclusion, we defined genes and processes related to metabolism and virulence that operate in E. piscicida to facilitate survival in low- and high-temperature seawater, which may underlie the infection outbreak mechanisms of E. piscicida and facilitate the development of improved vaccines against marine pathogens.IMPORTANCE Edwardsiella piscicida is one of the most important marine pathogens and causes serious edwardsiellosis in farmed fish during the summer-autumn seasonal changes, resulting in enormous losses to aquaculture industries worldwide. Survival and transmission of the pathogen in seawater are critical steps that increase the risk of outbreaks. To investigate the mechanism of survival in seawater for this marine pathogen, we used transposon insertion sequencing analysis to explore the fitness determinants in summer and autumn seawater. Approximately 127 genes linked to metabolism and virulence, as well as other processes, were revealed in E. piscicida to contribute to better adaptations to the seasonal alternations of seawater environments; these genes provide important insights into the infection outbreak mechanisms of E. piscicida and potential improved treatments or vaccines against marine pathogens.

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

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

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

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

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

Sensing of cytosolic LPS through caspy2 pyrin domain mediates noncanonical inflammasome activation in zebrafish

The noncanonical inflammasome is critical for cytosolic sensing of Gram-negative pathogens. Here, we show that bacterial infection induces caspy2 activation in zebrafish fibroblasts, which mediates pyroptosis via a caspase-5-like activity. Zebrafish caspy2 binds directly to lipopolysaccharide via the N-terminal pyrin death domain, resulting in caspy2 oligomerization, which is critical for pyroptosis. Furthermore, we show that caspy2 is highly expressed in the zebrafish gut and is activated during infection. Knockdown of caspy2 expression impairs the ability of zebrafish to restrict bacterial invasion in vivo, and protects larvae from lethal sepsis. Collectively, our results identify a crucial event in the evolution of pattern recognition into the death domain superfamily-mediated intracellular lipopolysaccharide-sensing pathway in innate immunity.

In humans, caspase-5 is an LPS sensor that can induce gasdermin D cleavage and pyroptosis. Here, the authors show that zebrafish caspy2 is a functional homolog as it also senses cytosolic LPS to activate the noncanonical inflammasome and to protect against bacterial infection, but it does so via pyrin death domain interactions.