A comprehensive transposon mutant library of Francisella novicida, a bioweapon surrogate

rLVS ΔcapB/Yp F1-V single vector platform vaccine expressing Yersinia pestis F1 and LcrV antigens provides complete protection against lethal respiratory challenge with virulent plague bacilli

Yersinia pestis, the causative agent of plague, is classified as a Tier I Select Agent of bioterrorism and is among a few pathogens of high concern as a potential cause of a future pandemic. Currently, there is no licensed vaccine against plague. Previously, we developed a live attenuated vaccine candidate, rLVS ΔcapB/Yp F1-V, that utilizes a highly attenuated capB mutant of Francisella tularensis Live Vaccine Strain as a vector to express a fusion protein of Y. pestis F1 and LcrV antigens. We showed that homologous prime-boost vaccination with this vaccine provided potent protection in mice against lethal respiratory challenge with virulent Y. pestis. Here, we report on the immunogenicity and efficacy of rLVS ΔcapB/Yp F1-V and additional LVS ΔcapB-vectored vaccine candidates in mice. We demonstrate that three homologous prime-boost immunizations with an optimized dose of rLVS ΔcapB/Yp F1-V provided complete protection against pneumonic plague in a stringent mouse model, outperforming other candidates and matching the survival efficacy of the toxic and unlicensed live attenuated Y. pestis EV76 strain vaccine; moreover, mice immunized with the rLVS ΔcapB/Yp F1-V vaccine had minimal weight loss post-challenge that was significantly less than mice immunized with the EV76 vaccine. Protection induced by rLVS ΔcapB/Yp F1-V correlates with F1 and LcrV-specific serum antibody levels. Our results highlight the potential of rLVS ΔcapB/Yp F1-V to address the unmet need for a plague vaccine.

Outer membrane tube formation by Francisella novicida involves extensive envelope modifications and is linked with type VI secretion and alterations to the host phagosomal membrane

Francisella tularensis is a gram-negative, intracellular pathogen that causes the zoonotic disease tularemia. Due to its ease of dissemination and high lethality, F. tularensis is classified as a tier 1 select agent with potential for misuse as a bioweapon. The mechanisms by which Francisella replicates intracellularly and interacts with the host during infection are not well understood. Francisella produces spherical outer membrane vesicles (OMVs) and novel tubular extensions of its cell surface that are also released extracellularly. These OMV and outer membrane tubes (OMTs) contain Francisella virulence factors and are produced in response to amino acid starvation and during infection of macrophages. To investigate how the OMTs are formed, we used cryogenic electron tomography to examine the model Francisella spp., Francisella novicida, during in vitro culture and within the macrophage phagosome. OMT formation involved progressive alterations of the bacterial envelope, resulting in extensions of both the inner and outer membranes. A dynamic cytoplasmic structure was present at the base of the OMT that extended into the tubes during elongation, together with cytoplasmic material. OMT produced within the macrophage phagosome was associated with changes to the phagosomal membrane, suggesting a role in phagosomal escape. Consistent with this, using confocal microscopy, we observed co-localization of the Francisella type VI secretion system with the OMT, both within bacteria and in released tubular vesicles. These findings reveal the cellular transformations that occur during membrane tubulation by Francisella and provide insights into the function of membrane-derived structures during host-pathogen interactions.

IMPORTANCE

Francisella tularensis is an intracellular bacterial pathogen that causes the zoonotic disease tularemia. Following uptake by host cells, the bacteria rapidly escape the phagosome and replicate intracellularly. In previous studies, we found that Francisella produces tubular extensions of its cell surface in response to specific cues and during macrophage infection. In the present study, we used cryogenic electron tomography to examine tube formation by the model Francisella sp., F. novicida. This analysis revealed that tube formation involves extensive bacterial envelope alterations and a dynamic cytoplasmic organelle. Furthermore, tubes produced by bacteria within infected macrophages were associated with the breakdown of the phagosomal membrane. In addition, we found that the Francisella type VI secretion system, which is essential for phagosomal escape, co-localized with the bacterial tubes. These findings reveal the cellular transformations that occur during membrane tubulation by Francisella and suggest a role for the tubes in phagosomal escape.

Mobilizable shuttle vectors with fluorescent markers functional across different species of bacteria

Chromophore-containing proteins (CCPs), including fluorescent and non-fluorescent (chromoproteins), have been widely used in microbiological research. However, several roadblocks often limit their use in non-model bacterial species, including efficient transformation, suitable plasmid origins of replication, and optimal promoter choice. Here, we have engineered a set of 32 shuttle plasmids designed to overcome these roadblocks in an effort to streamline this process for future research. We have selected eight different CCPs: eforCP, YukonOFP, DasherGFP, tinsel Purple, aeBlue, FuGFP, super-folder GFP, and super-folder Cherry2. To broaden the potential host range, we utilized two distinct backbones with p15a either fused to a Francisella origin (FnOri) or to the broad host origin RSF1010 and included a transfer origin (oriT) to facilitate transformation via conjugation. Moreover, we have created versions of each vector, which confer resistance to either kanamycin or chloramphenicol. Lastly, to enable promoter-swapping, we engineered the constitutive pJ23100 promoter element to be flanked by BsaI sites, thereby enabling promoter exchange by the Golden Gate assembly to evaluate CCP expression with different host promoters. To demonstrate the usability of the pKEK-Chrom plasmid series, we evaluated their expression in Escherichia coli, Shewanella oneidensis, and Vibrio alginolyticus. We further demonstrated the utility of promoter swapping in Francisella novicida and validated the functionality of the RSF1010 origin in Acinetobacter baumannii. In summary, the pKEK-Chrom plasmid series provides a palette of different CCPs that streamline their use in non-model gram-negative bacteria.

IMPORTANCE

Chromophore-containing proteins (CCPs), including both fluorescent proteins and pigment-producing (non-fluorescent) chromoproteins, have become invaluable tools for microbial research. However, their successful implementation in understudied bacterial species lacking established genetic tools often requires substantial time and resources. Our goal was to develop a set of plasmid-based vectors that could streamline CCP expression in gram-negative bacteria. To do so, we developed a set of 32 plasmid vectors, the pKEK-Chrom plasmid series, specifically designed to facilitate CCP expression across different bacteria, including Escherichia coli, Vibrio alginolyticus, Shewanella oneidensis, Francisella novicida, and Acinetobacter baumannii.

TolC facilitates the intracellular survival and immunomodulation of Salmonella Typhi in human host cells

Salmonella enterica serovar Typhi (S. Typhi) causes typhoid fever, a systemic infection that affects millions of people worldwide. S. Typhi can invade and survive within host cells, such as intestinal epithelial cells and macrophages, by modulating their immune responses. However, the immunomodulatory capability of S. Typhi in relation to TolC-facilitated efflux pump function remains unclear. The role of TolC, an outer membrane protein that facilitates efflux pump function, in the invasion and immunomodulation of S. Typhi, was studied in human intestinal epithelial cells and macrophages. The tolC deletion mutant of S. Typhi was compared with the wild-type and its complemented strain in terms of their ability to invade epithelial cells, survive and induce cytotoxicity in macrophages, and elicit proinflammatory cytokine production in macrophages. The tolC mutant, which has a defective outer membrane, was impaired in invading epithelial cells compared to the wild-type strain, but the intracellular presence of the tolC mutant exhibited greater cytotoxicity and induced higher levels of proinflammatory cytokines (IL-1β and IL-8) in macrophages compared to the wild-type strain. These effects were reversed by complementing the tolC mutant with a functional tolC gene. Our results suggest that TolC plays a role in S. Typhi to efficiently invade epithelial cells and suppress host immune responses during infection. TolC may be a potential target for the development of novel therapeutics against typhoid fever.

Tolfenpyrad Derivatives Exhibit Potent Francisella-Specific Antibacterial Activity without Toxicity to Mammalian Cells In Vitro

Tularemia is a deadly disease caused by Francisella tularensis, an emerging intracellular bacterial pathogen that can be disseminated rapidly through aerosols and vector-borne transmission. Recent surveillance data demonstrate an increasing incidence in several countries. Although clinical isolates of Francisella strains are sensitive to currently used antibiotics, engineered or horizontal acquisition of antibiotic resistance is a constant threat to public health. Therefore, the identification of antibiotics that target previously undrugged pathways is required to safeguard human health. An environmental pesticide that is registered for use in multiple countries, tolfenpyrad, shows promising activity to block Francisella growth; however, it is not a suitable antimicrobial candidate for use in vivo due to potential toxicity in humans and other animals. In this study, we applied a structure-activity relationship approach to tolfenpyrad to generate compounds with improved antibacterial activity and reduced toxicity. Through screening of a library of derivatives, we identified analogs with improved therapeutic windows compared with tolfenpyrad. Although structural diversity exists among these analogs, they inhibit the growth of Francisella species but not other Gram-negative or Gram-positive species. These compounds block intramacrophage growth of F. novicida and pathogenesis in an in vivo arthropod model of infection. Although the biochemical activity of these drugs is unknown, they appear to target the same pathway as the parent molecule because F. novicida mutants that are resistant to tolfenpyrad are also resistant to its analogs. Taken together, these findings suggest that these tolfenpyrad-derived compounds comprise a new class of Francisella-targeted antimicrobials and merit further evaluation and development.

TolC and EmrA1 contribute to Francisella novicida multidrug resistance and modulation of host cell death

Francisella spp. are Gram-negative, facultative intracellular pathogens. Francisella tularensis causes the human disease tularemia and is considered a biological threat agent due to its high infectivity and virulence. A central aspect of Francisella virulence is its ability to dampen host immune responses. We previously identified the outer membrane channel (OMC) protein TolC as a critical F. tularensis virulence factor required for suppression of apoptotic and proinflammatory responses during macrophage infection. TolC functions as part of multidrug efflux systems and the type I secretion pathway that exports bacterial effector proteins. In these systems, TolC forms tripartite complexes together with an inner membrane transporter and periplasmic membrane fusion protein (MFP). To advance understanding of TolC function in Francisella, we analyzed OMC and MFP homologs in Francisella novicida, a widely used model species that causes a tularemia-like disease in mice. In agreement with the previous F. tularensis studies, all three OMCs present in F. novicida contributed to multidrug resistance, but only TolC was important for suppressing macrophage cell death. In addition, we identified the EmrA1 MFP as important for resisting antimicrobial compounds and dampening host cell death. In contrast to results obtained with F. tularensis, the cell death triggered during infection with the F. novicida tolC and emrA1 mutants was dominated by pyroptosis rather than apoptosis. These data expand our understanding of TolC function in Francisella and underscore both conserved and differential aspects of F. novicida and F. tularensis.

IMPORTANCE

Francisella tularensis is a Gram-negative intracellular bacterial pathogen and causative agent of tularemia. We previously identified the outer membrane channel protein TolC as contributing to antimicrobial resistance and subversion of host responses by F. tularensis. To advance understanding of TolC function in Francisella and to identify components that might work together with TolC, we took advantage of a transposon mutant library in F. novicida, a model species that causes a tularemia-like disease in mice. Our findings identify TolC and the membrane fusion protein EmrA1 as important for both antimicrobial resistance and suppression of macrophage cell death. This study also revealed differences in cell death pathways triggered by F. novicida versus F. tularensis infection that may relate to differences in virulence.

Identifying essential genes in Schaalia odontolytica using a highly-saturated transposon library

The unique epibiotic-parasitic relationship between Nanosynbacter lyticus type strain TM7x, a member of the newly identified Candidate Phyla Radiation, now referred to as Patescibacteria, and its basibiont, Schaalia odontolytica strain XH001 (formerly Actinomyces odontolyticus), require more powerful genetic tools for deeper understanding of the genetic underpinnings that mediate their obligate relationship. Previous studies have mainly characterized the genomic landscape of XH001 during or post TM7x infection through comparative genomic or transcriptomic analyses followed by phenotypic analysis. Comprehensive genetic dissection of the pair is currently cumbersome due to the lack of robust genetic tools in TM7x. However, basic genetic tools are available for XH001 and this study expands the current genetic toolset by developing high-throughput transposon insertion sequencing (Tn-seq). Tn-seq was employed to screen for essential genes in XH001 under laboratory conditions. A highly saturated Tn-seq library was generated with nearly 660,000 unique insertion mutations, averaging one insertion every 2-3 nucleotides. 203 genes, 10.5% of the XH001 genome, were identified as putatively essential.

Integrating bacterial molecular genetics with chemical biology for renewed antibacterial drug discovery

The application of dyes to understanding the aetiology of infection inspired antimicrobial chemotherapy and the first wave of antibacterial drugs. The second wave of antibacterial drug discovery was driven by rapid discovery of natural products, now making up 69% of current antibacterial drugs. But now with the most prevalent natural products already discovered, ∼107 new soil-dwelling bacterial species must be screened to discover one new class of natural product. Therefore, instead of a third wave of antibacterial drug discovery, there is now a discovery bottleneck. Unlike natural products which are curated by billions of years of microbial antagonism, the vast synthetic chemical space still requires artificial curation through the therapeutics science of antibacterial drugs - a systematic understanding of how small molecules interact with bacterial physiology, effect desired phenotypes, and benefit the host. Bacterial molecular genetics can elucidate pathogen biology relevant to therapeutics development, but it can also be applied directly to understanding mechanisms and liabilities of new chemical agents with new mechanisms of action. Therefore, the next phase of antibacterial drug discovery could be enabled by integrating chemical expertise with systematic dissection of bacterial infection biology. Facing the ambitious endeavour to find new molecules from nature or new-to-nature which cure bacterial infections, the capabilities furnished by modern chemical biology and molecular genetics can be applied to prospecting for chemical modulators of new targets which circumvent prevalent resistance mechanisms.

ECDEP: identifying essential proteins based on evolutionary community discovery and subcellular localization

BACKGROUND

In cellular activities, essential proteins play a vital role and are instrumental in comprehending fundamental biological necessities and identifying pathogenic genes. Current deep learning approaches for predicting essential proteins underutilize the potential of gene expression data and are inadequate for the exploration of dynamic networks with limited evaluation across diverse species.

RESULTS

We introduce ECDEP, an essential protein identification model based on evolutionary community discovery. ECDEP integrates temporal gene expression data with a protein-protein interaction (PPI) network and employs the 3-Sigma rule to eliminate outliers at each time point, constructing a dynamic network. Next, we utilize edge birth and death information to establish an interaction streaming source to feed into the evolutionary community discovery algorithm and then identify overlapping communities during the evolution of the dynamic network. SVM recursive feature elimination (RFE) is applied to extract the most informative communities, which are combined with subcellular localization data for classification predictions. We assess the performance of ECDEP by comparing it against ten centrality methods, four shallow machine learning methods with RFE, and two deep learning methods that incorporate multiple biological data sources on Saccharomyces. Cerevisiae (S. cerevisiae), Homo sapiens (H. sapiens), Mus musculus, and Caenorhabditis elegans. ECDEP achieves an AP value of 0.86 on the H. sapiens dataset and the contribution ratio of community features in classification reaches 0.54 on the S. cerevisiae (Krogan) dataset.

CONCLUSIONS

Our proposed method adeptly integrates network dynamics and yields outstanding results across various datasets. Furthermore, the incorporation of evolutionary community discovery algorithms amplifies the capacity of gene expression data in classification.

Diverse roles of low-molecular weight thiol GSH in Francisella's virulence, location sensing and GSH-stealing from host

Low-molecular weight (LMW) thiols, encompassing peptides and small proteins with active cysteine residue(s), are important to bacteria as they are involved in a wide range of redox reactions. They include the tripeptide glutathione (GSH) and the small redox proteins, thioredoxins and glutaredoxins. We review the low MW thiols and related molecules in Francisella species and what role they may play in growth and virulence. Genes for GSH biosynthesis, metabolism and thioredoxins are present in all strains of Francisella, including the fully human-virulent strains. GSH and cysteine (CSH) are the major LMW thiols in Francisella extracts. We explore the potential role of the LMW thiols to overcome the nutritional challenges of intracellular growth (high GSH conditions) as well as the nutritional challenges of planktonic growth (low GSH conditions), and their contribution to Francisella's sensing its environmental location. Francisella may also use GSH as a source of CSH, for which it is auxotrophic. "Glutathione stealing" from the host may be an important part of Francisella's success strategy as a facultative intracellular pathogen both to detect its location and obtain CSH. An understanding of GSH metabolism in Francisella provides insights into the interaction of this pathogen with its host and may reveal additional targets for therapeutic intervention for tularemia infections.

Tolfenpyrad displays Francisella-targeted antibiotic activity that requires an oxidative stress response regulator for sensitivity

IMPORTANCE

Francisella species are highly pathogenic bacteria that pose a threat to global health security. These bacteria can be made resistant to antibiotics through facile methods, and we lack a safe and protective vaccine. Given their history of development as bioweapons, new treatment options must be developed to bolster public health preparedness. Here, we report that tolfenpyrad, a pesticide that is currently in use worldwide, effectively inhibits the growth of Francisella. This drug has an extensive history of use and a plethora of safety and toxicity data, making it a good candidate for development as an antibiotic. We identified mutations in Francisella novicida that confer resistance to tolfenpyrad and characterized a transcriptional regulator that is required for sensitivity to both tolfenpyrad and reactive oxygen species.

Benchmarking tools for detecting longitudinal differential expression in proteomics data allows establishing a robust reproducibility optimization regression approach

Quantitative proteomics has matured into an established tool and longitudinal proteomics experiments have begun to emerge. However, no effective, simple-to-use differential expression method for longitudinal proteomics data has been released. Typically, such data is noisy, contains missing values, and has only few time points and biological replicates. To address this need, we provide a comprehensive evaluation of several existing differential expression methods for high-throughput longitudinal omics data and introduce a Robust longitudinal Differential Expression (RolDE) approach. The methods are evaluated using over 3000 semi-simulated spike-in proteomics datasets and three large experimental datasets. In the comparisons, RolDE performs overall best; it is most tolerant to missing values, displays good reproducibility and is the top method in ranking the results in a biologically meaningful way. Furthermore, RolDE is suitable for different types of data with typically unknown patterns in longitudinal expression and can be applied by non-experienced users.

LI-Detector: a Method for Curating Ordered Gene-Replacement Libraries

In recent years the availability of genome sequence information has grown logarithmically resulting in the identification of a plethora of uncharacterized genes. To address this gap in functional annotation, many high-throughput screens have been devised to uncover novel gene functions. Gene-replacement libraries are one such tool that can be screened in a high-throughput way to link genotype and phenotype and are key community resources. However, for a phenotype to be attributed to a specific gene, there needs to be confidence in the genotype. Construction of large libraries can be laborious and occasionally errors will arise. Here, we present a rapid and accurate method for the validation of any ordered library where a gene has been replaced or disrupted by a uniform linear insertion (LI). We applied our method (LI-detector) to the well-known Keio library of Escherichia coli gene-deletion mutants. Our method identified 3,718 constructed mutants out of a total of 3,728 confirmed isolates, with a success rate of 99.7% for identifying the correct kanamycin cassette position. This data set provides a benchmark for the purity of the Keio mutants and a screening method for mapping the position of any linear insertion, such as an antibiotic resistance cassette in any ordered library. IMPORTANCE The construction of ordered gene replacement libraries requires significant investment of time and resources to create a valuable community resource. During construction, technical errors may result in a limited number of incorrect mutants being made. Such mutants may confound the output of subsequent experiments. Here, using the remarkable E. coli Keio knockout library, we describe a method to rapidly validate the construction of every mutant.

Mycobacterium abscessus Mutants with a Compromised Functional Link between the Type VII ESX-3 System and an Iron Uptake Mechanism Reliant on an Unusual Mycobactin Siderophore

The opportunistic pathogen Mycobacterium abscessus subsp. abscessus (Mab) has become an emerging public health threat due to the increasing number of Mab-associated chronic pulmonary disease cases. Treatment requires multiple drug courses and is often combined with surgical resection. Cure rates are only ~50% due to treatment failure and comorbidities. Deeper understanding of the biology of Mab is required to illuminate potential avenues for the development of better therapeutics against Mab infections. The ESX-3 type VII protein secretion system of Mab has an important role in host inflammatory and pathological responses during infection. In this work, we demonstrate a functional link between ESX-3 and an iron uptake system based on an unusual mycobactin-type siderophore (designated MBT Ab) and exploit this link to implement a large screen for transposon mutants with an impaired ESX-3. Most mutants we identified carry insertions in genes encoding predicted ESX-3 secretion machinery components or potential ESX-3 substrates. The mutants overproduce MBT Ab, a trait consistent with an iron uptake defect. Our characterization of MBT Ab revealed structural features reminiscent of nocardial mycobactin-like compounds with cytotoxicity. This finding raises the possibility that MBT Ab may play roles in pathogenesis unlinked to iron homeostasis. The mutants generated herein will facilitate research to better understand the role of ESX-3 and its interplay with the siderophore system.

The type IV pili component PilO is a virulence determinant of Francisella novicida

Francisella tularensis is a highly pathogenic intracellular bacterium that causes the disease tularemia. While its ability to replicate within cells has been studied in much detail, the bacterium also encodes a less characterised type 4 pili (T4P) system. T4Ps are dynamic adhesive organelles identified as major virulence determinants in many human pathogens. In F. tularensis, the T4P is required for adherence to the host cell, as well as for protein secretion. Several components, including pilins, a pili peptidase, a secretin pore and two ATPases, are required to assemble a functional T4P, and these are encoded within distinct clusters on the Francisella chromosome. While some of these components have been functionally characterised, the role of PilO, if any, still is unknown. Here, we examined the role of PilO in the pathogenesis of F. novicida. Our results show that the PilO is essential for pilus assembly on the bacterial surface. In addition, PilO is important for adherence of F. novicida to human monocyte-derived macrophages, secretion of effector proteins and intracellular replication. Importantly, the pilO mutant is attenuated for virulence in BALB/c mice regardless of the route of infection. Following intratracheal and intradermal infection, the mutant caused no histopathology changes, and demonstrated impaired phagosomal escape and replication within lung liver as well as spleen. Thus, PilO is an essential virulence determinant of F. novicida.

Testing assembly strategies of Francisella tularensis genomes to infer an evolutionary conservation analysis of genomic structures

Background

We benchmarked sequencing technology and assembly strategies for short-read, long-read, and hybrid assemblers in respect to correctness, contiguity, and completeness of assemblies in genomes of Francisella tularensis. Benchmarking allowed in-depth analyses of genomic structures of the Francisella pathogenicity islands and insertion sequences. Five major high-throughput sequencing technologies were applied, including next-generation “short-read” and third-generation “long-read” sequencing methods.

Results

We focused on short-read assemblers, hybrid assemblers, and analysis of the genomic structure with particular emphasis on insertion sequences and the Francisella pathogenicity island. The A5-miseq pipeline performed best for MiSeq data, Mira for Ion Torrent data, and ABySS for HiSeq data from eight short-read assembly methods. Two approaches were applied to benchmark long-read and hybrid assembly strategies: long-read-first assembly followed by correction with short reads (Canu/Pilon, Flye/Pilon) and short-read-first assembly along with scaffolding based on long reads (Unicyler, SPAdes). Hybrid assembly can resolve large repetitive regions best with a “long-read first” approach.

Conclusions

Genomic structures of the Francisella pathogenicity islands frequently showed misassembly. Insertion sequences (IS) could be used to perform an evolutionary conservation analysis. A phylogenetic structure of insertion sequences and the evolution within the clades elucidated the clade structure of the highly conservative F. tularensis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08115-x.

The Biosynthetic Pathway of Ubiquinone Contributes to Pathogenicity of Francisella novicida

Francisella tularensis is the causative agent of tularemia. Because of its extreme infectivity and high mortality rate, this pathogen was classified as a biothreat agent. Francisella spp. are strict aerobes, and ubiquinone (UQ) has been previously identified in these bacteria. While the UQ biosynthetic pathways were extensively studied in Escherichia coli, allowing the identification of 15 Ubi proteins to date, little is known about Francisella spp. In this study, and using Francisella novicida as a surrogate organism, we first identified ubiquinone 8 (UQ₈) as the major quinone found in the membranes of this bacterium. Next, we characterized the UQ biosynthetic pathway in F. novicida using a combination of bioinformatics, genetics, and biochemical approaches. Our analysis disclosed the presence in Francisella of 10 putative Ubi proteins, and we confirmed 8 of them by heterologous complementation in E. coli. The UQ biosynthetic pathways from F. novicida and E. coli share similar patterns. However, differences were highlighted: the decarboxylase remains unidentified in Francisella spp., and homologs of the Ubi proteins involved in the O₂-independent UQ pathway are not present. This is in agreement with the strictly aerobic niche of this bacterium. Next, via two approaches, i.e., the use of an inhibitor (3-amino-4-hydroxybenzoic acid) and a transposon mutant, both of which strongly impair the synthesis of UQ, we demonstrated that UQ is essential for the growth of F. novicida in respiratory medium and contributes to its pathogenicity in Galleria mellonella used as an alternative animal model. IMPORTANCE Francisella tularensis is the causative bacterium of tularemia and is classified as a biothreat agent. Using multidisciplinary approaches, we investigated the ubiquinone (UQ) biosynthetic pathway that operates in F. novicida used as a surrogate. We show that UQ₈ is the major quinone identified in the membranes of Francisella novicida. We identified a new competitive inhibitor that strongly decreased the biosynthesis of UQ. Our demonstration of the crucial roles of UQ for the respiratory metabolism of F. novicida and for the involvement in its pathogenicity in the Galleria mellonella model should stimulate the search for selective inhibitors of bacterial UQ biosynthesis.

Disruption of the Francisella noatunensis subsp. orientalis pdpA Gene Results in Virulence Attenuation and Protection in Zebrafish

Several Francisella spp., including Francisella noatunensis, are regarded as important emerging pathogens of wild and farmed fish. However, very few studies have investigated the virulence factors that allow these bacterial species to be pathogenic in fish. The Francisella pathogenicity island (FPI) is a well-described, gene-dense region encoding major virulence factors for the genus Francisella. pdpA is a member of the pathogenicity-determining protein genes carried by the FPI that are implicated in the ability of the mammalian pathogen Francisella tularensis to escape and replicate in infected host cells. Using a sacB suicide approach, we generated pdpA knockouts to address the role of PdpA as a virulence factor for F. noatunensis. Because polarity can be an issue in gene-dense regions, we generated two different marker-based mutants in opposing polarity (the F. noatunensis subsp. orientalis ΔpdpA1 and ΔpdpA2 strains). Both mutants were attenuated (P < 0.0001) in zebrafish challenges and displayed impaired intracellular replication (P < 0.05) and cytotoxicity (P < 0.05), all of which could be restored to wild-type (WT) levels by complementation for the ΔpdpA1 mutant. Importantly, differences were found for bacterial burden and induction of acute-phase and proinflammatory genes for the F. noatunensis subsp. orientalis ΔpdpA1 and ΔpdpA2 mutants compared to the WT during acute infection. In addition, neither mutant resulted in significant histopathological changes. Finally, immunization with the F. noatunensis subsp. orientalis ΔpdpA1 mutant led to protection (P < 0.012) against an acute 40% lethal dose (LD40) challenge with WT F. noatunensis in the zebrafish model of infection. Taken together, the results from this study further demonstrate physiological similarities within the genus Francisella relative to their phylogenetic relationships and the utility of zebrafish for addressing virulence factors for the genus.

An iteration model for identifying essential proteins by combining comprehensive PPI network with biological information

Background

Essential proteins have great impacts on cell survival and development, and played important roles in disease analysis and new drug design. However, since it is inefficient and costly to identify essential proteins by using biological experiments, then there is an urgent need for automated and accurate detection methods. In recent years, the recognition of essential proteins in protein interaction networks (PPI) has become a research hotspot, and many computational models for predicting essential proteins have been proposed successively.

Results

In order to achieve higher prediction performance, in this paper, a new prediction model called TGSO is proposed. In TGSO, a protein aggregation degree network is constructed first by adopting the node density measurement method for complex networks. And simultaneously, a protein co-expression interactive network is constructed by combining the gene expression information with the network connectivity, and a protein co-localization interaction network is constructed based on the subcellular localization data. And then, through integrating these three kinds of newly constructed networks, a comprehensive protein–protein interaction network will be obtained. Finally, based on the homology information, scores can be calculated out iteratively for different proteins, which can be utilized to estimate the importance of proteins effectively. Moreover, in order to evaluate the identification performance of TGSO, we have compared TGSO with 13 different latest competitive methods based on three kinds of yeast databases. And experimental results show that TGSO can achieve identification accuracies of 94%, 82% and 72% out of the top 1%, 5% and 10% candidate proteins respectively, which are to some degree superior to these state-of-the-art competitive models.

Conclusions

We constructed a comprehensive interactive network based on multi-source data to reduce the noise and errors in the initial PPI, and combined with iterative methods to improve the accuracy of necessary protein prediction, and means that TGSO may be conducive to the future development of essential protein recognition as well.

Prediction of Essential Genes in Comparison States Using Machine Learning

Identifying essential genes in comparison states (EGS) is vital to understanding cell differentiation, performing drug discovery, and identifying disease causes. Here, we present a machine learning method termed Prediction of Essential Genes in Comparison States (PreEGS). To capture the alteration of the network in comparison states, PreEGS extracts topological and gene expression features of each gene in a five-dimensional vector. PreEGS also recruits a positive sample expansion method to address the problem of unbalanced positive and negative samples, which is often encountered in practical applications. Different classifiers are applied to the simulated datasets, and the PreEGS based on the random forests model (PreEGSRF) was chosen for optimal performance. PreEGSRF was then compared with six other methods, including three machine learning methods, to predict EGS in a specific state. On real datasets with four gene regulatory networks, PreEGSRF predicted five essential genes related to leukemia and five enriched KEGG pathways. Four of the predicted essential genes and all predicted pathways were consistent with previous studies and highly correlated with leukemia. With high prediction accuracy and generalization ability, PreEGSRF is broadly applicable for the discovery of disease-causing genes, driver genes for cell fate decisions, and complex biomarkers of biological systems.

Targeted genetic screening in bacteria with a Cas12k-guided transposase

Microbes employ sophisticated cellular networks encoded by complex genomes to rapidly adapt to changing environments. High-throughput genome engineering methods are valuable tools for functionally profiling genotype-phenotype relationships and understanding the complexity of cellular networks. However, current methods either rely on special homologous recombination systems and are thus applicable in only limited bacterial species or can generate only nonspecific mutations and thus require extensive subsequent screening. Here, we report a site-specific transposon-assisted genome engineering (STAGE) method that allows high-throughput Cas12k-guided mutagenesis in various microorganisms, such as Pseudomonas aeruginosa and Klebsiella pneumoniae. Exploiting the powerful STAGE technique, we construct a site-specific transposon mutant library that focuses on all possible transcription factors (TFs) in P. aeruginosa, enabling the comprehensive identification of essential genes and antibiotic-resistance-related factors. Given its broad host range activity and easy programmability, this method can be widely adapted to diverse microbial species for rapid genome engineering and strain evolution.

The pentose phosphate pathway constitutes a major metabolic hub in pathogenic Francisella

Metabolic pathways are now considered as intrinsic virulence attributes of pathogenic bacteria and thus represent potential targets for antibacterial strategies. Here we focused on the role of the pentose phosphate pathway (PPP) and its connections with other metabolic pathways in the pathophysiology of Francisella novicida. The involvement of the PPP in the intracellular life cycle of Francisella was first demonstrated by studying PPP inactivating mutants. Indeed, we observed that inactivation of the tktA, rpiA or rpe genes severely impaired intramacrophage multiplication during the first 24 hours. However, time-lapse video microscopy demonstrated that rpiA and rpe mutants were able to resume late intracellular multiplication. To better understand the links between PPP and other metabolic networks in the bacterium, we also performed an extensive proteo-metabolomic analysis of these mutants. We show that the PPP constitutes a major bacterial metabolic hub with multiple connections to glycolysis, the tricarboxylic acid cycle and other pathways, such as fatty acid degradation and sulfur metabolism. Altogether our study highlights how PPP plays a key role in the pathogenesis and growth of Francisella in its intracellular niche.

Francisella FlmX broadly affects lipopolysaccharide modification and virulence

The outer membrane protects Gram-negative bacteria from the host environment. Lipopolysaccharide (LPS), a major outer membrane constituent, has distinct components (lipid A, core, O-antigen) generated by specialized pathways. In this study, we describe the surprising convergence of these pathways through FlmX, an uncharacterized protein in the intracellular pathogen Francisella. FlmX is in the flippase family, which includes proteins that traffic lipid-linked envelope components across membranes. flmX deficiency causes defects in lipid A modification, core remodeling, and O-antigen addition. We find that an F. tularensis mutant lacking flmX is >1,000,000-fold attenuated. Furthermore, FlmX is required to resist the innate antimicrobial LL-37 and the antibiotic polymyxin. Given FlmX's central role in LPS modification and its conservation in intracellular pathogens Brucella, Coxiella, and Legionella, FlmX may represent a novel drug target whose inhibition could cripple bacterial virulence and sensitize bacteria to innate antimicrobials and antibiotics.

Genetic Determinants of Antibiotic Resistance in Francisella

Tularemia, caused by Francisella tularensis, is endemic to the northern hemisphere. This zoonotic organism has historically been developed into a biological weapon. For this Tier 1, Category A select agent, it is important to expand our understanding of its mechanisms of antibiotic resistance (AMR). Francisella is unlike many Gram-negative organisms in that it does not have significant plasmid mobility, and does not express AMR mechanisms on plasmids; thus plasmid-mediated resistance does not occur naturally. It is possible to artificially introduce plasmids with AMR markers for cloning and gene expression purposes. In this review, we survey both the experimental research on AMR in Francisella and bioinformatic databases which contain genomic and proteomic data. We explore both the genetic determinants of intrinsic AMR and naturally acquired or engineered antimicrobial resistance as well as phenotypic resistance in Francisella. Herein we survey resistance to beta-lactams, monobactams, carbapenems, aminoglycosides, tetracycline, polymyxins, macrolides, rifampin, fosmidomycin, and fluoroquinolones. We also highlight research about the phenotypic AMR difference between planktonic and biofilm Francisella. We discuss newly developed methods of testing antibiotics against Francisella which involve the intracellular nature of Francisella infection and may better reflect the eventual clinical outcomes for new antibiotic compounds. Understanding the genetically encoded determinants of AMR in Francisella is key to optimizing the treatment of patients and potentially developing new antimicrobials for this dangerous intracellular pathogen.

Defined Mutant Library Sequencing (DML-Seq) for Identification of Conditional Essential Genes

Transposon insertion sequencing (TIS) is an emerging technique which utilizes a massive transposon mutant library to screen specific phenotype and determine the conditional essential genetic requirements for bacterial fitness under distinct conditions combined with high-throughput parallel sequencing technology. Compared with a massive mutant library in traditional TIS, the defined mutant library sequencing (DML-Seq) has advantages as: 1) efficient mutagenesis; 2) low bottleneck effects; 3) avoid hotpots caused by screening; 4) can be directly used in the following experiments. Here, we described an optimized procedure of DML-Seq for fitness screen to supply classical TIS using the marine pathogenic bacterium Edwardsiella piscicida as an example.

Utilization of a tetracycline-inducible system for high-level expression of recombinant proteins in Francisella tularensis LVS

Francisella tularensis is a Gram-negative intracellular pathogen causing tularemia. A number of its potential virulence factors have been identified, but their biology and functions are not precisely known. Understanding the biological and immunological functions of these proteins requires adequate genetic tools for homologous and heterologous expression of cloned genes, maintaining both original structure and post-translational modifications. Here, we report the construction of a new multipurpose shuttle plasmid - pEVbr - which can be used for high-level expression in F. tularensis. The pEVbr plasmid has been constructed by modifying the TetR-regulated expression vector pEDL17 (LoVullo, 2012) that includes (i) a strong F. tularensis bfr promoter, and (ii) two tet operator sequences cloned into the promoter. The cloned green fluorescent protein (GFP), used as a reporter, demonstrated almost undetectable basal expression level under uninduced conditions and a highly dynamic dose-dependent response to the inducer. The utility of the system was further confirmed by cloning the gapA and FTT_1676 genes into the pEVbr vector and quantifying proteins expression in F. tularensis LVS, as well as by studying post-translational modification of the cloned genes. This study demonstrates that high levels of recombinant native-like Francisella proteins can be produced in Francisella cells. Hence, this system may be beneficial for the analysis of protein function and the development of new treatments and vaccines.

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.

The Sensor Kinase QseC Regulates the Unlinked PmrA Response Regulator and Downstream Gene Expression in Francisella

The facultative intracellular bacterial pathogen Francisella tularensis is the causative agent of tularemia in humans and animals. Gram-negative bacteria utilize two-component regulatory systems (TCS) to sense and respond to their changing environment. No classical, tandemly arranged sensor kinase and response regulator TCS genes exist in the human virulent Francisella tularensis subsp. tularensis, but orphaned members are present. PmrA is an orphan response regulator responsible for intramacrophage growth and virulence; however, the regulation of PmrA activity is not understood. We and others have shown that PmrA represses the expression of priM, described to encode an antivirulence determinant. By screening a mutant library for increased priM promoter activity, we identified the sensor kinase homolog QseC as an upstream regulator of priM expression, and this regulation is in part dependent upon the aspartate phosphorylation site of PmrA (D51). Several examined environmental signals, including epinephrine, which is reported to activate QseC in other bacteria, do not affect priM expression in a manner dependent on PmrA. Intramacrophage survival assays also question the finding that PriM is an antivirulence factor. Thus, these data suggest that the PmrA-regulated gene priM is modulated by the QseC-PmrA (QseB) TCS in Francisella IMPORTANCE The disease tularemia is caused by the highly infectious Gram-negative pathogen Francisella tularensis This bacterium encodes few regulatory factors (e.g., two-component systems [TCS]). PmrA, required for intramacrophage survival and virulence in the mouse model, is encoded by an orphan TCS response regulator gene. It is unclear how PmrA is responsive to environmental signals to regulate loci, including the PmrA-repressed gene priM We identify an orphan sensor kinase (QseC) that is required for priM repression and further explore both environmental signals that might regulate the QseC-PmrA TCS and the function of PriM.

Expression and protease characterization of a conserved protein YgjD in Vibrio harveyi

The glycopeptidase GCP and its homologue proteins are conserved and essential for survival of bacteria. The ygjD gene (Glycopeptidase homologue) was cloned from Vibrio harveyi strain SF-1. The gene consisted of 1,017 bp, which encodes a 338 amino acid polypeptide. The nucleotide sequence similarity of the ygjD gene with that of V. harveyi FDAARGOS 107 was 95%. The ygjD gene also showed similarities of 68%, 67% and 50% with those of Salmonella enterica, Escherichia coli and Bacillus cereus. The ygjD gene was expressed in E. coli BL21 (DE3) and the recombinant YgjD was purified by Ni2+ affinity chromatography column. The purified YgjD showed a specific 37 kDa band on sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and exhibited protease activities of 59,000 units/mg, 53,700 units/mg and 8,100 units/mg, respectively, on N-Acetyl-L-tyrosine ethyl ester monohydrate (ATEE), N-Benzoyl-L-tyrosine ethyl ester (BTEE) and N-Benzoyl-DL-arginine-4-nitroanilide hydrochloride (BAPNA) substrates. When the conserved amino acids of His111, Glu113 and His115 in the YgjD were replaced with alanine, respectively, the protease activities of the mutants were partly decreased. The two conserved His111 and His115 of YgjD were mutated and the protein lost the protease activity, which implied that the two amino acid played very important roles in maintaining its protease activity. The addition of the purified YgjD to the culture medium of V. harveyi strain SF-1 can effectively promote the bacteria growth. These results indicated that the protease activities may be involved in the survival of bacteria.

Network-based methods for predicting essential genes or proteins: a survey

Genes that are thought to be critical for the survival of organisms or cells are called essential genes. The prediction of essential genes and their products (essential proteins) is of great value in exploring the mechanism of complex diseases, the study of the minimal required genome for living cells and the development of new drug targets. As laboratory methods are often complicated, costly and time-consuming, a great many of computational methods have been proposed to identify essential genes/proteins from the perspective of the network level with the in-depth understanding of network biology and the rapid development of biotechnologies. Through analyzing the topological characteristics of essential genes/proteins in protein-protein interaction networks (PINs), integrating biological information and considering the dynamic features of PINs, network-based methods have been proved to be effective in the identification of essential genes/proteins. In this paper, we survey the advanced methods for network-based prediction of essential genes/proteins and present the challenges and directions for future research.

Francisella novicida Two-Component System Response Regulator BfpR Modulates iglC Gene Expression, Antimicrobial Peptide Resistance, and Biofilm Production

Response regulators are a critical part of the two-component system of gene expression regulation in bacteria, transferring a signal from a sensor kinase into DNA binding activity resulting in alteration of gene expression. In this study, we investigated a previously uncharacterized response regulator in Francisella novicida, FTN_1452 that we have named BfpR (Biofilm-regulating Francisella protein Regulator, FTN_1452). In contrast to another Francisella response regulator, QseB/PmrA, BfpR appears to be a negative regulator of biofilm production, and also a positive regulator of antimicrobial peptide resistance in this bacterium. The protein was crystallized and X-ray crystallography studies produced a 1.8 Å structure of the BfpR N-terminal receiver domain revealing interesting insight into its potential interaction with the sensor kinase. Structural analysis of BfpR places it in the OmpR/PhoP family of bacterial response regulators along with WalR and ResD. Proteomic and transcriptomic analyses suggest that BfpR overexpression affects expression of the critical Francisella virulence factor iglC, as well as other proteins in the bacterium. We demonstrate that mutation of bfpR is associated with an antimicrobial peptide resistance phenotype, a phenotype also associated with other response regulators, for the human cathelicidin peptide LL-37 and a sheep antimicrobial peptide SMAP-29. F. novicida with mutated bfpR replicated better than WT in intracellular infection assays in human-derived macrophages suggesting that the down-regulation of iglC expression in bfpR mutant may enable this intracellular replication to occur. Response regulators have been shown to play important roles in the regulation of bacterial biofilm production. We demonstrate that F. novicida biofilm formation was highly increased in the bfpR mutant, corresponding to altered glycogen synthesis. Waxworm infection experiments suggest a role of BfpR as a negative modulator of iglC expression with de-repression by Mg²⁺. In this study, we find that the response regulator BfpR may be a negative regulator of biofilm formation, and a positive regulator of antimicrobial peptide resistance in F. novicida.

Thioredoxin Modulates Cell Surface Hydrophobicity in Acinetobacter baumannii

Acinetobacter baumannii, a Gram-negative coccobacillus, has become a prevalent nosocomial health threat affecting the majority of hospitals both in the U.S. and around the globe. Microbial cell surface hydrophobicity (CSH) has previously been correlated with virulence, uptake by immune cells, and attachment to epithelial cells. A mutant strain of A. baumannii (ΔtrxA) lacking the redox protein thioredoxin A was found to be more hydrophobic than its wild type (WT) and complemented counterparts, as measured by both Microbial Adhesion to Hydrocarbon (MATH) and salt aggregation. The hydrophobicity of the mutant could be abrogated through treatment with sodium cyanoborohydride (SCBH). This modulation correlated with reduction of disulfide bonds, as SCBH was able to reduce 5,5′-dithio-bis-[2-nitrobenzoic acid] and treatment with the known disulfide reducer, β-mercaptoethanol, also decreased ΔtrxA CSH. Additionally, the ΔtrxA mutant was more readily taken up than WT by J774 macrophages and this differential uptake could be abrogated though SCBH treatment. When partitioned into aqueous and hydrophobic phases, ΔtrxA recovered from the hydrophobic partition was phagocytosed more readily than from the aqueous phase further supporting the contribution of CSH to A. baumannii uptake by phagocytes. A second Gram-negative bacterium, Francisella novicida, also showed the association of TrxA deficiency (FnΔtrxA) with increased hydrophobicity and uptake by J774 cells. We previously have demonstrated that modification of the type IV pilus system (T4P) was associated with the A. baumannii ΔtrxA phenotype, and the Francisella FnΔtrxA mutant also was found to have a marked T4P deficiency. Interestingly, a F. novicida mutant lacking pilT also showed increased hydrophobicity over FnWT. Collective evidence presented in this study suggests that Gram-negative bacterial thioredoxin mediates CSH through multiple mechanisms including disulfide-bond reduction and T4P modulation.

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.

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.

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.

Identifying Bacterial Essential Genes Based on a Feature-Integrated Method

Essential genes are those genes of an organism that are considered to be crucial for its survival. Identification of essential genes is therefore of great significance to advance our understanding of the principles of cellular life. We have developed a novel computational method, which can effectively predict bacterial essential genes by extracting and integrating homologous features, protein domain feature, gene intrinsic features, and network topological features. By performing the principal component regression (PCR) analysis for Escherichia coli MG1655, we established a classification model with the average area under curve (AUC) value of 0.992 in ten times 5-fold cross-validation tests. Furthermore, when employing this new model to a distantly related organism-Streptococcus pneumoniae TIGR4, we still got a reliable AUC value of 0.788. These results indicate that our feature-integrated approach could have practical applications in accurately investigating essential genes from broad bacterial species, and also provide helpful guidelines for the minimal cell.

The Lipid A 1-Phosphatase, LpxE, Functionally Connects Multiple Layers of Bacterial Envelope Biogenesis

Dephosphorylation of the lipid A 1-phosphate by LpxE in Gram-negative bacteria plays important roles in antibiotic resistance, bacterial virulence, and modulation of the host immune system. Our results demonstrate that in addition to removing the 1-phosphate from lipid A, LpxEs also dephosphorylate undecaprenyl pyrophosphate, an important metabolite for the synthesis of the essential envelope components, peptidoglycan and O-antigen. Therefore, LpxEs participate in multiple layers of biogenesis of the Gram-negative bacterial envelope and increase antibiotic resistance. This discovery marks an important step toward understanding the regulation and biogenesis of the Gram-negative bacterial envelope.

Although distinct lipid phosphatases are thought to be required for processing lipid A (component of the outer leaflet of the outer membrane), glycerophospholipid (component of the inner membrane and the inner leaflet of the outer membrane), and undecaprenyl pyrophosphate (C₅₅-PP; precursors of peptidoglycan and O antigens of lipopolysaccharide) in Gram-negative bacteria, we report that the lipid A 1-phosphatases, LpxEs, functionally connect multiple layers of cell envelope biogenesis in Gram-negative bacteria. We found that Aquifex aeolicus LpxE structurally resembles YodM in Bacillus subtilis, a phosphatase for phosphatidylglycerol phosphate (PGP) with a weak in vitro activity on C₅₅-PP, and rescues Escherichia coli deficient in PGP and C₅₅-PP phosphatase activities; deletion of lpxE in Francisella novicida reduces the MIC value of bacitracin, indicating a significant contribution of LpxE to the native bacterial C₅₅-PP phosphatase activity. Suppression of plasmid-borne lpxE in F. novicida deficient in chromosomally encoded C₅₅-PP phosphatase activities results in cell enlargement, loss of O-antigen repeats of lipopolysaccharide, and ultimately cell death. These discoveries implicate LpxE as the first example of a multifunctional regulatory enzyme that orchestrates lipid A modification, O-antigen production, and peptidoglycan biogenesis to remodel multiple layers of the Gram-negative bacterial envelope.

A Francisella novicida Mutant, Lacking the Soluble Lytic Transglycosylase Slt, Exhibits Defects in Both Growth and Virulence

Francisella tularensis is the causative agent of tularemia and has gained recent interest as it poses a significant biothreat risk. F. novicida is commonly used as a laboratory surrogate for tularemia research due to genetic similarity and susceptibility of mice to infection. Currently, there is no FDA-approved tularemia vaccine, and identifying therapeutic targets remains a critical gap in strategies for combating this pathogen. Here, we investigate the soluble lytic transglycosylase or Slt in F. novicida, which belongs to a class of peptidoglycan-modifying enzymes known to be involved in cell division. We assess the role of Slt in biology and virulence of the organism as well as the vaccine potential of the slt mutant. We show that the F. novicida slt mutant has a significant growth defect in acidic pH conditions. Further microscopic analysis revealed significantly altered cell morphology compared to wild-type, including larger cell size, extensive membrane protrusions, and cell clumping and fusion, which was partially restored by growth in neutral pH or genetic complementation. Viability of the mutant was also significantly decreased during growth in acidic medium, but not at neutral pH. Furthermore, the slt mutant exhibited significant attenuation in a murine model of intranasal infection and virulence could be restored by genetic complementation. Moreover, we could protect mice using the slt mutant as a live vaccine strain against challenge with the parent strain; however, we were not able to protect against challenge with the fully virulent F. tularensis Schu S4 strain. These studies demonstrate a critical role for the Slt enzyme in maintaining proper cell division and morphology in acidic conditions, as well as replication and virulence in vivo. Our results suggest that although the current vaccination strategy with F. novicida slt mutant would not protect against Schu S4 challenges, the Slt enzyme could be an ideal target for future therapeutic development.

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.

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.

A model to predict the function of hypothetical proteins through a nine-point classification scoring schema

BACKGROUND

Hypothetical proteins [HP] are those that are predicted to be expressed in an organism, but no evidence of their existence is known. In the recent past, annotation and curation efforts have helped overcome the challenge in understanding their diverse functions. Techniques to decipher sequence-structure-function relationship, especially in terms of functional modelling of the HPs have been developed by researchers, but using the features as classifiers for HPs has not been attempted. With the rise in number of annotation strategies, next-generation sequencing methods have provided further understanding the functions of HPs.

RESULTS

In our previous work, we developed a six-point classification scoring schema with annotation pertaining to protein family scores, orthology, protein interaction/association studies, bidirectional best BLAST hits, sorting signals, known databases and visualizers which were used to validate protein interactions. In this study, we introduced three more classifiers to our annotation system, viz. pseudogenes linked to HPs, homology modelling and non-coding RNAs associated to HPs. We discuss the challenges and performance of these classifiers using machine learning heuristics with an improved accuracy from Perceptron (81.08 to 97.67), Naive Bayes (54.05 to 96.67), Decision tree J48 (67.57 to 97.00), and SMO_npolyk (59.46 to 96.67).

CONCLUSION

With the introduction of three new classification features, the performance of the nine-point classification scoring schema has an improved accuracy to functionally annotate the HPs.

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.

ClpB mutants of Francisella tularensis subspecies holarctica and tularensis are defective for type VI secretion and intracellular replication

Francisella tularensis, a highly infectious, intracellular bacterium possesses an atypical type VI secretion system (T6SS), which is essential for the virulence of the bacterium. Recent data suggest that the HSP100 family member, ClpB, is involved in T6SS disassembly in the subspecies Francisella novicida. Here, we investigated the role of ClpB for the function of the T6SS and for phenotypic characteristics of the human pathogenic subspecies holarctica and tularensis. The ∆clpB mutants of the human live vaccine strain, LVS, belonging to subspecies holarctica, and the highly virulent SCHU S4 strain, belonging to subspecies tularensis, both showed extreme susceptibility to heat shock and low pH, severely impaired type VI secretion (T6S), and significant, but impaired intracellular replication compared to the wild-type strains. Moreover, they showed essentially intact phagosomal escape. Infection of mice demonstrated that both ΔclpB mutants were highly attenuated, but the SCHU S4 mutant showed more effective replication than the LVS strain. Collectively, our data demonstrate that ClpB performs multiple functions in the F. tularensis subspecies holarctica and tularensis and its function is important for T6S, intracellular replication, and virulence.

Single vector platform vaccine protects against lethal respiratory challenge with Tier 1 select agents of anthrax, plague, and tularemia

Bacillus anthracis, Yersinia pestis, and Francisella tularensis are the causative agents of Tier 1 Select Agents anthrax, plague, and tularemia, respectively. Currently, there are no licensed vaccines against plague and tularemia and the licensed anthrax vaccine is suboptimal. Here we report F. tularensis LVS ΔcapB (Live Vaccine Strain with a deletion in capB)- and attenuated multi-deletional Listeria monocytogenes (Lm)-vectored vaccines against all three aforementioned pathogens. We show that LVS ΔcapB- and Lm-vectored vaccines express recombinant B. anthracis, Y. pestis, and F. tularensis immunoprotective antigens in broth and in macrophage-like cells and are non-toxic in mice. Homologous priming-boosting with the LVS ΔcapB-vectored vaccines induces potent antigen-specific humoral and T-cell-mediated immune responses and potent protective immunity against lethal respiratory challenge with all three pathogens. Protection against anthrax was far superior to that obtained with the licensed AVA vaccine and protection against tularemia was comparable to or greater than that obtained with the toxic and unlicensed LVS vaccine. Heterologous priming-boosting with LVS ΔcapB- and Lm-vectored B. anthracis and Y. pestis vaccines also induced potent protective immunity against lethal respiratory challenge with B. anthracis and Y. pestis. The single vaccine platform, especially the LVS ΔcapB-vectored vaccine platform, can be extended readily to other pathogens.

Transcription Elongation Factor GreA Plays a Key Role in Cellular Invasion and Virulence of Francisella tularensis subsp. novicida

Francisella tularensis is a facultative intracellular Gram-negative bacterium that causes the zoonotic disease tularemia. We identified the transcription elongation factor GreA as a virulence factor in our previous study, but its role was not defined. Here, we investigate the effects of the inactivation of the greA gene, generating a greA mutant of F. tularensis subsp. novicida. Inactivation of greA impaired the bacterial invasion into and growth within host cells, and subsequently virulence in mouse infection model. A transcriptomic analysis (RNA-Seq) showed that the loss of GreA caused the differential expression of 196 bacterial genes, 77 of which were identified as virulence factors in previous studies. To confirm that GreA regulates the expression of virulence factors involved in cell invasion by Francisella, FTN_1186 (pepO) and FTN_1551 (ampD) gene mutants were generated. The ampD deletion mutant showed reduced invasiveness into host cells. These results strongly suggest that GreA plays an important role in the pathogenesis of Francisella by affecting the expression of virulence genes and provide new insights into the complex regulation of Francisella infection.

A high-throughput screening campaign to identify inhibitors of DXP reductoisomerase (IspC) and MEP cytidylyltransferase (IspD)

The rise of antibacterial resistance among human pathogens represents a problem that could change the landscape of healthcare unless new antibiotics are developed. The methyl erythritol phosphate (MEP) pathway represents an attractive series of targets for novel antibiotic design, considering each enzyme of the pathway is both essential and has no human homologs. Here we describe a pilot scale high-throughput screening (HTS) campaign against the first and second committed steps in the pathway, catalyzed by DXP reductoisomerase (IspC) and MEP cytidylyltransferase (IspD), using compounds present in the commercially available LOPAC1280 library as well as in an in-house natural product extract library. Hit compounds were characterized to deduce their mechanism of inhibition; most function through aggregation. The HTS workflow outlined here is useful for quickly screening a chemical library, while effectively identifying false positive compounds associated with assay constraints and aggregation.

Zinc Acquisition Mechanisms Differ between Environmental and Virulent Francisella Species

Zinc is an essential nutrient for bacterial growth. Because host cells can restrict pathogen access to zinc as an antimicrobial defense mechanism, intracellular pathogens such as Francisella must sense their environment and acquire zinc in response. In many bacteria, the conserved transcription factor Zur is a key regulator of zinc acquisition. To identify mechanisms of zinc uptake in Francisella novicida U112, transcriptome sequencing of wild-type and putative zur mutant bacteria was performed. Only three genes were confirmed as directly regulated by Zur and zinc limitation by quantitative reverse transcription-PCR. One of these genes, FTN_0879, is predicted to encode a protein with similarity to the zupT family of zinc transporters, which are not typically regulated by Zur. While a putative znuACB operon encoding a high-affinity zinc transporter was identified in U112, expression of this operon was not controlled by Zur or zinc concentration. Disruption of zupT but not znuA in U112 impaired growth under zinc limitation, suggesting that ZupT is the primary mechanism of zinc acquisition under these conditions. In the virulent Francisella tularensis subsp. tularensis Schu S4 strain, zupT is a pseudogene, and attempts to delete znuA were unsuccessful, suggesting that it is essential in this strain. A reverse TetR repression system was used to knock down the expression of znuA in Schu S4, revealing that znuA is required for growth under zinc limitation and contributes to intracellular growth within macrophages. Overall, this work identifies genes necessary for adaptation to zinc limitation and highlights nutritional differences between environmental and virulent Francisella strains.

IMPORTANCEFrancisella tularensis is a tier 1 select agent with a high potential for lethality and no approved vaccine. A better understanding of Francisella virulence factors is required for the development of therapeutics. While acquisition of zinc has been shown to be required for the virulence of numerous intracellular pathogens, zinc uptake has not been characterized in Francisella. This work characterizes the Zur regulon in F. novicida and identifies two transporters that contribute to bacterial growth under zinc limitation. In addition, these data identify differences in mechanisms of zinc uptake and tolerance to zinc limitation between F. tularensis and F. novicida, highlighting the role of znuA in the growth of Schu S4 under zinc limitation.

Amino acid deprivation and central carbon metabolism regulate the production of outer membrane vesicles and tubes by Francisella

Francisella tularensis is a highly virulent Gram-negative bacterial pathogen that causes the zoonotic disease tularemia. F. novicida, a model tularemia strain, produces spherical outer membrane vesicles (OMV), as well as novel tubular vesicles and extensions of the cell surface. These OMV and tubes (OMV/T) are produced in a regulated manner and contain known virulence factors. Mechanisms by which bacterial vesicles are produced and regulated are not well understood. We performed a genetic screen in F. novicida to decipher the molecular basis for regulated OMV/T formation, and identified both hypo- and hyper-vesiculating mutants. Mutations in fumA and tktA, involved in central carbon metabolism, and in FTN_0908 and FTN_1037, of unknown function, resulted in severe defects in OMV/T production. Cysteine deprivation was identified as the signal that triggers OMV/T formation in F. novicida during growth in rich medium. We also found that fully virulent F. tularensis produces OMV/T in a similarly regulated manner. Further analysis revealed that OMV/T production is responsive to deprivation of essential amino acids in addition to cysteine, and that the hypo-vesiculating mutants are defective in responding to this signal. Thus, amino acid starvation, such as encountered by Francisella during host cell invasion, regulates the production of membrane-derived structures.

Bacterial lipoproteins and other factors released by Francisella tularensis modulate human neutrophil lifespan: Effects of a TLR1 SNP on apoptosis inhibition

Francisella tularensis infects several cell types including neutrophils, and aberrant neutrophil accumulation contributes to tissue destruction during tularaemia. We demonstrated previously that F. tularensis strains Schu S4 and live vaccine strain markedly delay human neutrophil apoptosis and thereby prolong cell lifespan, but the bacterial factors that mediate this aspect of virulence are undefined. Herein, we demonstrate that bacterial conditioned medium (CM) can delay apoptosis in the absence of direct infection. Biochemical analyses show that CM contained F. tularensis surface factors as well as outer membrane components. Our previous studies excluded roles for lipopolysaccharide and capsule in apoptosis inhibition, and current studies of [¹⁴C] acetate‐labelled bacteria argue against a role for other bacterial lipids in this process. At the same time, studies of isogenic mutants indicate that TolC and virulence factors whose expression requires FevR or MglA were also dispensable, demonstrating that apoptosis inhibition does not require Type I or Type VI secretion. Instead, we identified bacterial lipoproteins (BLPs) as active factors in CM. Additional studies of isolated BLPs demonstrated dose‐dependent neutrophil apoptosis inhibition via a TLR2‐dependent mechanism that is significantly influenced by a common polymorphism, rs5743618, in human TLR1. These data provide fundamental new insight into pathogen manipulation of neutrophil lifespan and BLP function.

AR-13, a Celecoxib Derivative, Directly Kills Francisella In Vitro and Aids Clearance and Mouse Survival In Vivo

Francisella tularensis (F. tularensis) is the causative agent of tularemia and is classified as a Tier 1 select agent. No licensed vaccine is currently available in the United States and treatment of tularemia is confined to few antibiotics. In this study, we demonstrate that AR-13, a derivative of the cyclooxygenase-2 inhibitor celecoxib, exhibits direct in vitro bactericidal killing activity against Francisella including a type A strain of F. tularensis (SchuS4) and the live vaccine strain (LVS), as well as toward the intracellular proliferation of LVS in macrophages, without causing significant host cell toxicity. Identification of an AR-13-resistant isolate indicates that this compound has an intracellular target(s) and that efflux pumps can mediate AR-13 resistance. In the mouse model of tularemia, AR-13 treatment protected 50% of the mice from lethal LVS infection and prolonged survival time from a lethal dose of F. tularensis SchuS4. Combination of AR-13 with a sub-optimal dose of gentamicin protected 60% of F. tularensis SchuS4-infected mice from death. Taken together, these data support the translational potential of AR-13 as a lead compound for the further development of new anti-Francisella agents.