Comparison of Francisella tularensis genomes reveals evolutionary events associated with the emergence of human pathogenic strains

Dual proteomics of infected macrophages reveal bacterial and host players involved in the Francisella intracellular life cycle and cell to cell dissemination by merocytophagy

Bacterial pathogens adapt and replicate within host cells, while host cells develop mechanisms to eliminate them. Using a dual proteomic approach, we characterized the intra-macrophage proteome of the facultative intracellular pathogen, Francisella novicida. More than 900 Francisella proteins were identified in infected macrophages after a 10-h infection. Biotin biosynthesis-related proteins were upregulated, emphasizing the role of biotin-associated genes in Francisella replication. Conversely, proteins encoded by the Francisella pathogenicity island (FPI) were downregulated, supporting the importance of the F. tularensis Type VI Secretion System for vacuole escape, not cytosolic replication. In the host cell, over 300 proteins showed differential expression among the 6200 identified during infection. The most upregulated host protein was cis-aconitate decarboxylase IRG1, known for itaconate production with antimicrobial properties in Francisella. Surprisingly, disrupting IRG1 expression did not impact Francisella's intracellular life cycle, suggesting redundancy with other immune proteins or inclusion in larger complexes. Over-representation analysis highlighted cell-cell contact and actin polymerization in macrophage deregulated proteins. Using flow cytometry and live cell imaging, we demonstrated that merocytophagy involves diverse cell-to-cell contacts and actin polymerization-dependent processes. These findings lay the groundwork for further exploration of merocytophagy and its molecular mechanisms in future research.Data are available via ProteomeXchange with identifier PXD035145.

Phenotypic and transcriptional characterization of F. tularensis LVS during transition into a viable but non-culturable state

Francisella tularensis is a gram-negative, intracellular pathogen which can cause serious, potentially fatal, illness in humans. Species of F. tularensis are found across the Northern Hemisphere and can infect a broad range of host species, including humans. Factors affecting the persistence of F. tularensis in the environment and its epidemiology are not well understood, however, the ability of F. tularensis to enter a viable but non-culturable state (VBNC) may be important. A broad range of bacteria, including many pathogens, have been observed to enter the VBNC state in response to stressful environmental conditions, such as nutrient limitation, osmotic or oxidative stress or low temperature. To investigate the transition into the VBNC state for F. tularensis, we analyzed the attenuated live vaccine strain, F. tularensis LVS grown under standard laboratory conditions. We found that F. tularensis LVS rapidly and spontaneously enters a VBNC state in broth culture at 37°C and that this transition coincides with morphological differentiation of the cells. The VBNC bacteria retained an ability to interact with both murine macrophages and human erythrocytes in in vitro assays and were insensitive to treatment with gentamicin. Finally, we present the first transcriptomic analysis of VBNC F. tularensis, which revealed clear differences in gene expression, and we identify sets of differentially regulated genes which are specific to the VBNC state. Identification of these VBNC specific genes will pave the way for future research aimed at dissecting the molecular mechanisms driving entry into the VBNC state.

Mutation of wbtJ, a N-formyltransferase involved in O-antigen synthesis, results in biofilm formation, phase variation and attenuation in Francisella tularensis

Two clinically important subspecies, Francisella tularensis subsp. tularensis (type A) and F. tularensis subsp. holarctica (type B) are responsible for most tularaemia cases, but these isolates typically form a weak biofilm under in vitro conditions. Phase variation of the F. tularensis lipopolysaccharide (LPS) has been reported in these subspecies, but the role of variation is unclear as LPS is crucial for virulence. We previously demonstrated that a subpopulation of LPS variants can constitutively form a robust biofilm in vitro, but it is unclear whether virulence was affected. In this study, we show that biofilm-forming variants of both fully virulent F. tularensis subspecies were highly attenuated in the murine tularaemia model by multiple challenge routes. Genomic sequencing was performed on these strains, which revealed that all biofilm-forming variants contained a lesion within the wbtJ gene, a formyltransferase involved in O-antigen synthesis. A ΔwbtJ deletion mutant recapitulated the biofilm, O-antigen and virulence phenotypes observed in natural variants and could be rescued through complementation with a functional wbtJ gene. Since the spontaneously derived biofilm-forming isolates in this study were a subpopulation of natural variants, reversion events to the wbtJ gene were detected that eliminated the phenotypes associated with biofilm variants and restored virulence. These results demonstrate a role for WbtJ in biofilm formation, LPS variation and virulence of F. tularensis.

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.

The rLVS ΔcapB/iglABC vaccine provides potent protection in Fischer rats against inhalational tularemia caused by various virulent Francisella tularensis strains

Francisella tularensis is one of the several biothreat agents for which a licensed vaccine is needed. To ensure vaccine protection is achieved across a range of virulent F. tularensis strains, we assembled and characterized a panel of F. tularensis isolates to be utilized as challenge strains. A promising tularemia vaccine candidate is rLVS ΔcapB/iglABC (rLVS), in which the vector is the LVS strain with a deletion in the capB gene and which additionally expresses a fusion protein comprising immunodominant epitopes of proteins IglA, IglB, and IglC. Fischer rats were immunized subcutaneously 1-3 times at 3-week intervals with rLVS at various doses. The rats were exposed to a high dose of aerosolized Type A strain Schu S4 (FRAN244), a Type B strain (FRAN255), or a tick derived Type A strain (FRAN254) and monitored for survival. All rLVS vaccination regimens including a single dose of 107 CFU rLVS provided 100% protection against both Type A strains. Against the Type B strain, two doses of 107 CFU rLVS provided 100% protection, and a single dose of 107 CFU provided 87.5% protection. In contrast, all unvaccinated rats succumbed to aerosol challenge with all of the F. tularensis strains. A robust Th1-biased antibody response was induced in all vaccinated rats against all F. tularensis strains. These results demonstrate that rLVS ΔcapB/iglABC provides potent protection against inhalational challenge with either Type A or Type B F. tularensis strains and should be considered for further analysis as a future tularemia vaccine.

Pathogenicity and virulence of Francisella tularensis

Tularaemia is a zoonotic disease caused by the Gram-negative bacterium, Francisella tularensis. Depending on its entry route into the organism, F. tularensis causes different diseases, ranging from life-threatening pneumonia to less severe ulceroglandular tularaemia. Various strains with different geographical distributions exhibit different levels of virulence. F. tularensis is an intracellular bacterium that replicates primarily in the cytosol of the phagocytes. The main virulence attribute of F. tularensis is the type 6 secretion system (T6SS) and its effectors that promote escape from the phagosome. In addition, F. tularensis has evolved a peculiar envelope that allows it to escape detection by the immune system. In this review, we cover tularaemia, different Francisella strains, and their pathogenicity. We particularly emphasize the intracellular life cycle, associated virulence factors, and metabolic adaptations. Finally, we present how F. tularensis largely escapes immune detection to be one of the most infectious and lethal bacterial pathogens.

Comparative analysis of absent in melanoma 2-inflammasome activation in Francisella tularensis and Francisella novicida

Francisella tularensis is a highly virulent Gram-negative bacterium that causes the fatal zoonotic disease tularemia. The mechanisms and signaling pathways leading to the absent in melanoma 2 (Aim2) inflammasome activation have been elegantly elucidated using Francisella novicida as a model. Although not pathogenic for humans, F. novicida can cause tularemia in mice, and the inflammatory response it triggers is the polar opposite to that observed in mice infected with F. tularensis strains. This study aimed to understand the mechanisms of Aim2 inflammasome activation in F. tularensis-infected macrophages. The results reveal that macrophages infected with the F. tularensis live vaccine strain (LVS) induce lower levels of Aim2-dependent IL-1β than those infected with F. novicida. The suppression/weak activation of Aim2 in F. tularensis LVS-infected macrophages is due to the suppression of the cGAS-STING DNA-sensing pathway. Furthermore, the introduction of exogenous F. tularensis LVS DNA into the cytosol of the F. tularensis LVS-infected macrophages, alone or in conjunction with a priming signal, failed to restore IL-1β levels similar to those observed for F. novicida-infected macrophages. These results indicated that, in addition to the bacterial DNA, DNA from some other sources, specifically from the damaged mitochondria, might contribute to the robust Aim2-dependent IL-1β levels observed in F. novicida-infected macrophages. The results indicate that F. tularensis LVS induces mitophagy that may potentially prevent the leakage of mitochondrial DNA and the subsequent activation of the Aim2 inflammasome. Collectively, this study demonstrates that the mechanisms of Aim2 inflammasome activation established for F. novicida are not operative in F. tularensis.

Why vary what's working? Phase variation and biofilm formation in Francisella tularensis

The notoriety of high-consequence human pathogens has increased in recent years and, rightfully, research efforts have focused on understanding host-pathogen interactions. Francisella tularensis has been detected in an impressively broad range of vertebrate hosts as well as numerous arthropod vectors and single-celled organisms. Two clinically important subspecies, F. tularensis subsp. tularensis (Type A) and F. tularensis subsp. holarctica (Type B), are responsible for the majority of tularemia cases in humans. The success of this bacterium in mammalian hosts can be at least partly attributed to a unique LPS molecule that allows the bacterium to avoid detection by the host immune system. Curiously, phase variation of the O-antigen incorporated into LPS has been documented in these subspecies of F. tularensis, and these variants often display some level of attenuation in infection models. While the role of phase variation in F. tularensis biology is unclear, it has been suggested that this phenomenon can aid in environmental survival and persistence. Biofilms have been established as the predominant lifestyle of many bacteria in the environment, though, it was previously thought that Type A and B isolates of F. tularensis typically form poor biofilms. Recent studies question this ideology as it was shown that alteration of the O-antigen allows robust biofilm formation in both Type A and B isolates. This review aims to explore the link between phase variation of the O-antigen, biofilm formation, and environmental persistence with an emphasis on clinically relevant subspecies and how understanding these poorly studied mechanisms could lead to new medical countermeasures to combat tularemia.

Identification of pyrC gene as an immunosuppressive factor in Francisella novicida infection

Francisella tularensis, a bacterial causative agent of the zoonosis tularemia, is highly pathogenic to humans. The pathogenicity of this bacterium is characterized by intracellular growth in immune cells, like macrophages, and host immune suppression. However, the detailed mechanism of immune suppression by F. tularensis is still unclear. To identify the key factors causing Francisella-mediated immunosuppression, large-scale screening using a transposon random mutant library containing 3552 mutant strains of F. tularensis subsp. novicida (F. novicida) was performed. Thirteen mutants that caused stronger tumor necrosis factor (TNF)-α production in infected U937 human macrophage cells than the wild-type F. novicida strain were isolated. Sequencing analysis of transposon insertion sites revealed 10 genes, including six novel genes, as immunosuppressive factors of Francisella. Among these, the relationship of the pyrC gene, which encodes dihydroorotase in the pyrimidine biosynthesis pathway, with Francisella-mediated immunosuppression was investigated. The pyrC deletion mutant strain (ΔpyrC) induced higher TNF-α production in U937 host cells than the wild-type F. novicida strain. The ΔpyrC mutant strain was also found to enhance host interleukin-1β and interferon (IFN)-β production. The heat-inactivated ΔpyrC mutant strain could not induce host TNF-α production. Moreover, the production of IFN-β resulting from ΔpyrC infection in U937 cells was repressed upon treatment with the stimulator of interferon genes (STING)-specific inhibitor, H-151. These results suggest that pyrC is related to the immunosuppressive activity and pathogenicity of Francisella via the STING pathway.

Atomic Structure of IglD Demonstrates Its Role as a Component of the Baseplate Complex of the Francisella Type VI Secretion System

Francisella tularensis, a Tier 1 select agent of bioterrorism, contains a type VI secretion system (T6SS) encoded within the Francisella pathogenicity island (FPI), which is critical for its pathogenesis. Among the 18 proteins encoded by FPI is IglD, which is essential to Francisella’s intracellular growth and virulence, but neither its location within T6SS nor its functional role has been established. Here, we present the cryoEM structure of IglD from Francisella novicida and show that the Francisella IglD forms a homotrimer that is structurally homologous to the T6SS baseplate protein TssK in Escherichia coli. Each IglD monomer consists of an N-terminal β-sandwich domain, a 4-helix bundle domain, and a flexible C-terminal domain. While the overall folds of IglD and TssK are similar, the two structures differ in three aspects: the relative orientation between their β-sandwich and the 4-helix bundle domains; two insertion loops present in TssK’s β-sandwich domain; and, consequently, a lack of subunit-subunit interaction between insertion loops in the IglD trimer. Phylogenetic analysis indicates that IglD is genetically remote from the TssK orthologs in other T6SSs. While the other components of the Francisella baseplate are unknown, we conducted pulldown assays showing IglJ interacts with IglD and IglH, pointing to a model wherein IglD, IglH, and IglJ form the baseplate of the Francisella T6SS. Alanine substitution mutagenesis further established that IglD’s hydrophobic pocket in the N-terminal β-sandwich domain interacts with two loops of IglJ, reminiscent of the TssK-TssG interaction. These results form a framework for understanding the hitherto unexplored Francisella T6SS baseplate.

Arginine Catabolism and Polyamine Biosynthesis Pathway Disparities Within Francisella tularensis Subpopulations

Francisella tularensis is a highly infectious zoonotic pathogen with as few as 10 organisms causing tularemia, a disease that is fatal if untreated. Although F. tularensis subspecies tularensis (type A) and subspecies holarctica (type B) share over 99.5% average nucleotide identity, notable differences exist in genomic organization and pathogenicity. The type A clade has been further divided into subtypes A.I and A.II, with A.I strains being recognized as some of the most virulent bacterial pathogens known. In this study, we report on major disparities that exist between the F. tularensis subpopulations in arginine catabolism and subsequent polyamine biosynthesis. The genes involved in these pathways include the speHEA and aguAB operons, along with metK. In the hypervirulent F. tularensis A.I clade, such as the A.I prototype strain SCHU S4, these genes were found to be intact and highly transcribed. In contrast, both subtype A.II and type B strains have a truncated speA gene, while the type B clade also has a disrupted aguA and truncated aguB. Ablation of the chromosomal speE gene that encodes a spermidine synthase reduced subtype A.I SCHU S4 growth rate, whereas the growth rate of type B LVS was enhanced. These results demonstrate that spermine synthase SpeE promotes faster replication in the F. tularensis A.I clade, whereas type B strains do not rely on this enzyme for in vitro fitness. Our ongoing studies on amino acid and polyamine flux within hypervirulent A.I strains should provide a better understanding of the factors that contribute to F. tularensis pathogenicity.

The intracellular pathogen Francisella tularensis escapes from adaptive immunity by metabolic adaptation

This study shows that metabolic adaptation allows the intracellular bacterial pathogen Francisella tularensis to escape recognition by the host adaptive immunity.

Intracellular pathogens lose many metabolic genes during their evolution from free-living bacteria, but the pathogenic consequences of their altered metabolic programs on host immunity are poorly understood. Here, we show that a pathogenic strain of Francisella tularensis subsp. tularensis (FT) has five amino acid substitutions in RibD, a converting enzyme of the riboflavin synthetic pathway responsible for generating metabolites recognized by mucosal-associated invariant T (MAIT) cells. Metabolites from a free-living strain, F. tularensis subsp. novicida (FN), activated MAIT cells in a T-cell receptor (TCR)–dependent manner, whereas introduction of FT-type ribD to the free-living strain was sufficient to attenuate this activation in both human and mouse MAIT cells. Intranasal infection in mice showed that the ribD^FT-expressing FN strain induced impaired Th1-type MAIT cell expansion and resulted in reduced bacterial clearance and worsened survival compared with the wild-type free-living strain FN. These results demonstrate that F. tularensis can acquire immune evasion capacity by alteration of metabolic programs during evolution.

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.

Cross-Species Proteomic Comparison of Outer Membrane Vesicles and Membranes of Francisella tularensis subsp. tularensis versus subsp. holarctica

Release of outer membrane vesicles (OMV) is an important phenomenon in Gram-negative bacteria playing multiple roles in their lifestyle, including in relation to virulence and host-pathogen interaction. Francisella tularensis, unlike other bacteria, releases unusually shaped, tubular OMV. We present a proteomic comparison of OMV and membrane fractions from two F. tularensis strains: moderately virulent subsp. holarctica strain FSC200 and highly virulent subsp. tularensis strain SchuS4. Proteomic comparison studies routinely evaluate samples from the same proteome, but sometimes we must compare samples from closely related organisms. This raises quantification issues. We propose a novel approach to cross-species proteomic comparison based on an intersection protein database from the individual single-species databases. This is less prone to quantification errors arising from differences in the sequences. Consecutively comparing subproteomes of OMV and membranes of the two strains allows distinguishing differences in relative protein amounts caused by global expression changes from those caused by preferential protein packing to OMV or membranes. Among the proteins most differently packed into OMV between the two strains, we detected proteins involved in biosynthesis and metabolism of bacterial envelope components like O-antigen, lipid A, phospholipids, and fatty acids, as well as some major structural outer membrane proteins. The data are available via ProteomeXchange with identifier PXD022406.

Reorganized Genomic Taxonomy of Francisellaceae Enables Design of Robust Environmental PCR Assays for Detection of Francisella tularensis

In recent years, an increasing diversity of species has been recognized within the family Francisellaceae. Unfortunately, novel isolates are sometimes misnamed in initial publications or multiple sources propose different nomenclature for genetically highly similar isolates. Thus, unstructured and occasionally incorrect information can lead to confusion in the scientific community. Historically, detection of Francisella tularensis in environmental samples has been challenging due to the considerable and unknown genetic diversity within the family, which can result in false positive results. We have assembled a comprehensive collection of genome sequences representing most known Francisellaceae species/strains and restructured them according to a taxonomy that is based on phylogenetic structure. From this structured dataset, we identified a small number of genomic regions unique to F. tularensis that are putatively suitable for specific detection of this pathogen in environmental samples. We designed and validated specific PCR assays based on these genetic regions that can be used for the detection of F. tularensis in environmental samples, such as water and air filters.

Alternative pathways utilize or circumvent putrescine for biosynthesis of putrescine-containing rhizoferrin

The siderophore rhizoferrin (N¹,N⁴-dicitrylputrescine) is produced in fungi and bacteria to scavenge iron. Putrescine-producing bacterium Ralstonia pickettii synthesizes rhizoferrin and encodes a single nonribosomal peptide synthetase-independent siderophore (NIS) synthetase. From biosynthetic logic, we hypothesized that this single enzyme is sufficient for rhizoferrin biosynthesis. We confirmed this by expression of R. pickettii NIS synthetase in Escherichia coli, resulting in rhizoferrin production. This was further confirmed in vitro using the recombinant NIS synthetase, synthesizing rhizoferrin from putrescine and citrate. Heterologous expression of homologous lbtA from Legionella pneumophila, required for rhizoferrin biosynthesis in that species, produced siderophore activity in E. coli. Rhizoferrin is also synthesized by Francisella tularensis and Francisella novicida, but unlike R. pickettii or L. pneumophila, Francisella species lack putrescine biosynthetic pathways because of genomic decay. Francisella encodes a NIS synthetase FslA/FigA and an ornithine decarboxylase homolog FslC/FigC, required for rhizoferrin biosynthesis. Ornithine decarboxylase produces putrescine from ornithine, but we show here in vitro that FigA synthesizes N-citrylornithine, and FigC is an N-citrylornithine decarboxylase that together synthesize rhizoferrin without using putrescine. We co-expressed F. novicida figA and figC in E. coli and produced rhizoferrin. A 2.1 Å X-ray crystal structure of the FigC N-citrylornithine decarboxylase reveals how the larger substrate is accommodated and how active site residues have changed to recognize N-citrylornithine. FigC belongs to a new subfamily of alanine racemase-fold PLP-dependent decarboxylases that are not involved in polyamine biosynthesis. These data reveal a natural product biosynthetic workaround that evolved to bypass a missing precursor and re-establish it in the final structure.

Control of Francisella tularensis Virulence at Gene Level: Network of Transcription Factors

Regulation of gene transcription is the initial step in the complex process that controls gene expression within bacteria. Transcriptional control involves the joint effort of RNA polymerases and numerous other regulatory factors. Whether global or local, positive or negative, regulators play an essential role in the bacterial cell. For instance, some regulators specifically modify the transcription of virulence genes, thereby being indispensable to pathogenic bacteria. Here, we provide a comprehensive overview of important transcription factors and DNA-binding proteins described for the virulent bacterium Francisella tularensis, the causative agent of tularemia. This is an unexplored research area, and the poorly described networks of transcription factors merit additional experimental studies to help elucidate the molecular mechanisms of pathogenesis in this bacterium, and how they contribute to disease.

Differences in Blood-derived Francisella tularensis Type B Strains from Clinical Cases of Tularemia

Francisella tularensis can cause the zoonotic disease tularemia and is partitioned into subspecies due to differences in chromosomal organization and virulence. The subspecies holarctica (type B) is generally considered more clonal than the other subpopulations with moderate virulence compared to the hypervirulent A.I clade. We performed whole genome sequencing (WGS) on six type B strains isolated from the blood of patients with tularemia within a one-year period from the same United States region, to better understand the associated pathogenicity. The WGS data were compared to the prototype strain for this subspecies, specifically FSC200, which was isolated from a patient with tularemia in Europe. These findings revealed 520–528 single nucleotide polymorphisms (SNPs) between the six United States type B strains compared to FSC200, with slightly higher A+T content in the latter strain. In contrast, comparisons between the six type B isolates showed that five of the six type B isolates had only 4–22 SNPs, while one of the strains had 47–53 SNPs. Analysis of SNPs in the core genome for the six United States type B isolates and the FSC200 strain gave similar results, suggesting that some of these mutations may have been nonsynonymous, resulting in altered protein function and pathogenicity.

Exploring the Diversity Within the Genus Francisella - An Integrated Pan-Genome and Genome-Mining Approach

Pan-genome analysis is a powerful method to explore genomic heterogeneity and diversity of bacterial species. Here we present a pan-genome analysis of the genus Francisella, comprising a dataset of 63 genomes and encompassing clinical as well as environmental isolates from distinct geographic locations. To determine the evolutionary relationship within the genus, we performed phylogenetic whole-genome studies utilizing the average nucleotide identity, average amino acid identity, core genes and non-recombinant loci markers. Based on the analyses, the phylogenetic trees obtained identified two distinct clades, A and B and a diverse cluster designated C. The sizes of the pan-, core-, cloud-, and shell-genomes of Francisella were estimated and compared to those of two other facultative intracellular pathogens, Legionella and Piscirickettsia. Francisella had the smallest core-genome, 692 genes, compared to 886 and 1,732 genes for Legionella and Piscirickettsia respectively, while the pan-genome of Legionella was more than twice the size of that of the other two genera. Also, the composition of the Francisella Type VI secretion system (T6SS) was analyzed. Distinct differences in the gene content of the T6SS were identified. In silico approaches performed to identify putative substrates of these systems revealed potential effectors targeting the cell wall, inner membrane, cellular nucleic acids as well as proteins, thus constituting attractive targets for site-directed mutagenesis. The comparative analysis performed here provides a comprehensive basis for the assessment of the phylogenomic relationship of members of the genus Francisella and for the identification of putative T6SS virulence traits.

Metabolic Profiling of Volatile Organic Compounds (VOCs) Emitted by the Pathogens Francisella tularensis and Bacillus anthracis in Liquid Culture

We conducted comprehensive (untargeted) metabolic profiling of volatile organic compounds (VOCs) emitted in culture by bacterial taxa Francisella tularensis (F. tularensis) subspecies novicida and Bacillus anthracis (B. anthracis) Sterne, surrogates for potential bacterial bioterrorism agents, as well as selective measurements of VOCs from their fully virulent counterparts, F. tularensis subspecies tularensis strain SCHU S4 and B. anthracis Ames. F. tularensis and B. anthracis were grown in liquid broth for time periods that covered logarithmic growth, stationary, and decline phases. VOCs emitted over the course of the growth phases were collected from the headspace above the cultures using solid phase microextraction (SPME) and were analyzed using gas chromatography-mass spectrometry (GC-MS). We developed criteria for distinguishing VOCs originating from bacteria versus background VOCs (originating from growth media only controls or sampling devices). Analyses of collected VOCs revealed methyl ketones, alcohols, esters, carboxylic acids, and nitrogen- and sulfur-containing compounds that were present in the bacterial cultures and absent (or present at only low abundance) in control samples indicating that these compounds originated from the bacteria. Distinct VOC profiles where observed for F. tularensis when compared with B. anthracis while the observed profiles of each of the two F. tularensis and B. anthracis strains exhibited some similarities. Furthermore, the relative abundance of VOCs was influenced by bacterial growth phase. These data illustrate the potential for VOC profiles to distinguish pathogens at the genus and species-level and to discriminate bacterial growth phases. The determination of VOC profiles lays the groundwork for non-invasive probes of bacterial metabolism and offers prospects for detection of microbe-specific VOC biomarkers from two potential biowarfare agents.

Proteome Wide Profiling of N-ε-Lysine Acetylation Reveals a Novel Mechanism of Regulation of the Chitinase Activity in Francisella novicida

Francisella tularensis is a Gram-negative bacterium that causes the zoonotic disease tularemia. The historical development of tularemia as a biological weapon has led to it being characterized by the CDC as a category A biothreat agent. Neither posttranslational modification (PTM) of proteins, in particular lysine acetylation, in Francisella nor its subsequent regulation of the protein activity has been well studied. In this work, we analyze N-ε-lysine acetylation of the F. tularensis ssp. novicida proteome by mass spectrometry for the first time. To create a comprehensive acetylation profile, we enriched protein acetylation using two approaches: (1) the addition of glucose or acetate into the culture medium and (2) direct chemical acetylation of N-ε-lysines with acetyl phosphate. We discovered 280 acetylated proteins with 1178 acetylation sites in the F. tularensis ssp. novicida strain U112. Lysine acetylation is an important PTM that regulates multiple cellular processes in bacteria, including metabolism, transcription, translation, stress response, and protein folding. We discovered that Francisella chitinases A and B are acetylated naturally and when chemically induced by acetyl phosphate. Moreover, chemical overacetylation of chitinases results in silencing of the enzymatic activity. Our findings suggest a novel mechanism of posttranslational regulation of the chitinase activity and that acetylation may play a role in Francisella's regulation of the protein activity.

Molecular Signatures Related to the Virulence of Bacillus cereus Sensu Lato, a Leading Cause of Devastating Endophthalmitis

Bacillus endophthalmitis is a devastating eye infection that causes rapid blindness through extracellular tissue-destructive exotoxins. Despite its importance, knowledge of the phylogenetic relationships and population structure of intraocular Bacillus spp. is lacking. In this study, we sequenced the whole genomes of eight Bacillus intraocular pathogens independently isolated from 8/52 patients with posttraumatic Bacillus endophthalmitis infections in the Eye Hospital of Wenzhou Medical University between January 2010 and December 2018. Phylogenetic analysis revealed that the pathogenic intraocular isolates belonged to Bacillus cereus, Bacillus thuringiensis and Bacillus toyonensis To determine the virulence of the ocular isolates, three representative strains were injected into mouse models, and severe endophthalmitis leading to blindness was observed. Through incorporating publicly available genomes for Bacillus spp., we found that the intraocular pathogens could be isolated independently but displayed a similar genetic context. In addition, our data provide genome-wide support for intraocular and gastrointestinal sources of Bacillus spp. belonging to different lineages. Importantly, we identified five molecular signatures of virulence and motility genes associated with intraocular infection, namely, plcA-2, InhA-3, InhA-4, hblA-5, and fliD using pangenome-wide association studies. The characterization of overrepresented genes in the intraocular isolates holds value to predict bacterial evolution and for the design of future intervention strategies in patients with endophthalmitis.IMPORTANCE In this study, we provided a detailed and comprehensive clinicopathological and pathogenic report of Bacillus endophthalmitis over the 8 years of the study period. We first reported the whole-genome sequence of Bacillus spp. causing devastating endophthalmitis and found that Bacillus toyonensis is able to cause endophthalmitis. Finally, we revealed significant endophthalmitis-associated virulence genes involved in hemolysis, immunity inhibition, and pathogenesis. Overall, as more sequencing data sets become available, these data will facilitate comparative research and will reveal the emergence of pathogenic "ocular bacteria."

Francisella tularensis subsp. holarctica Releases Differentially Loaded Outer Membrane Vesicles Under Various Stress Conditions

Francisella tularensis is a Gram-negative, facultative intracellular bacterium, causing a severe disease called tularemia. It secretes unusually shaped nanotubular outer membrane vesicles (OMV) loaded with a number of virulence factors and immunoreactive proteins. In the present study, the vesicles were purified from a clinical isolate of subsp. holarctica strain FSC200. We here provide a comprehensive proteomic characterization of OMV using a novel approach in which a comparison of OMV and membrane fraction is performed in order to find proteins selectively enriched in OMV vs. membrane. Only these proteins were further considered to be really involved in the OMV function and/or their exceptional structure. OMV were also isolated from bacteria cultured under various cultivation conditions simulating the diverse environments of F. tularensis life cycle. These included conditions mimicking the milieu inside the mammalian host during inflammation: oxidative stress, low pH, and high temperature (42°C); and in contrast, low temperature (25°C). We observed several-fold increase in vesiculation rate and significant protein cargo changes for high temperature and low pH. Further proteomic characterization of stress-derived OMV gave us an insight how the bacterium responds to the hostile environment of a mammalian host through the release of differentially loaded OMV. Among the proteins preferentially and selectively packed into OMV during stressful cultivations, the previously described virulence factors connected to the unique intracellular trafficking of Francisella were detected. Considerable changes were also observed in a number of proteins involved in the biosynthesis and metabolism of the bacterial envelope components like O-antigen, lipid A, phospholipids, and fatty acids. Data are available via ProteomeXchange with identifier PXD013074.

Can Insertion Sequences Proliferation Influence Genomic Plasticity? Comparative Analysis of Acinetobacter baumannii Sequence Type 78, a Persistent Clone in Italian Hospitals

Acinetobacter baumannii is a known opportunistic pathogen, dangerous for public health mostly due to its ability to rapidly acquire antibiotic-resistance traits. Its genome was described as characterized by remarkable plasticity, with a high frequency of homologous recombinations and proliferation of Insertion Sequences (IS). The SMAL pulsotype is an A. baumannii strain currently isolated only in Italy, characterized by a low incidence and a high persistence over the years. In this present work, we have conducted a comparative genomic analysis on this clone. The genome of 15 SMAL isolates was obtained and characterized in comparison with 24 other assemblies of evolutionary related isolates. The phylogeny highlighted the presence of a monophyletic clade (named ST78A), which includes the SMAL isolates. ST78A isolates have a low rate of homologous recombination and low gene content variability when compared to two related clades (ST78B and ST49) and to the most common A. baumannii variants worldwide (International Clones I and II). Remarkably, genomes in the ST78A clade present a high number of IS, including classes mostly absent in the other related genomes. Among these IS, one copy of IS66 was found to interrupt the gene comEC/rec2, involved in the acquisition of exogenous DNA. The genomic characterization of SMAL isolates shed light on the surprisingly low genomic plasticity and the high IS proliferation present in this strain. The interruption of the gene comEC/rec2 by an IS in the SMAL genomes brought us to formulate an evolutionary hypothesis according to which the proliferation of IS is slowing the acquisition of exogenous DNA, thus limiting genome plasticity. Such genomic architecture could explain the epidemiological behavior of high persistence and low incidence of the clone and provides an interesting framework to compare ST78 with the highly epidemic International Clones, characterized by high genomic plasticity.

Screen for fitness and virulence factors of Francisella sp. strain W12-1067 using amoebae

Francisella tularensis is the causative agent of the human disease referred to as tularemia. Other Francisella species are known but less is understood about their virulence factors. The role of environmental amoebae in the life-cycle of Francisella is still under discussion. Francisella sp. strain W12-1067 (F-W12) is an environmental Francisella isolate recently identified in Germany which is negative for the Francisella pathogenicity island, but exhibits a putative alternative type VI secretion system. Putative virulence factors have been identified in silico in the genome of F-W12. In this work, we established a "scatter screen", used earlier for pathogenic Legionella, to verify experimentally and identify candidate fitness factors using a transposon mutant bank of F-W12 and Acanthamoeba lenticulata as host organism. In these experiments, we identified 79 scatter clones (amoeba sensitive), which were further analyzed by an infection assay identifying 9 known virulence factors, but also candidate fitness factors of F-W12 not yet described as fitness factors in Francisella. The majority of the identified genes encoded proteins involved in the synthesis or maintenance of the cell envelope (LPS, outer membrane, capsule) or in the metabolism (glycolysis, gluconeogenesis, pentose phosphate pathway). Further ¹³C-flux analysis of the Tn5 glucokinase mutant strain revealed that the identified gene indeed encodes the sole active glucokinase in F-W12. In conclusion, candidate fitness factors of the new Francisella species F-W12 were identified using the scatter screen method which might also be usable for other Francisella species.

Citryl Ornithine Is an Intermediate in a Three-Step Biosynthetic Pathway for Rhizoferrin in Francisella

The Gram-negative bacterium Francisella tularensis secretes the siderophore rhizoferrin to scavenge necessary iron from the environment. Rhizoferrin, also produced by a variety of fungi and bacteria, comprises two citrate molecules linked by amide bonds to a central putrescine (diaminobutane) moiety. Genetic analysis has determined that rhizoferrin production in F. tularensis requires two enzymes: FslA, a siderophore synthetase of the nonribosomal peptide synthetase-independent siderophore synthetase (NIS) family, and FslC, a pyridoxal-phosphate-dependent decarboxylase. To discern the steps in the biosynthetic pathway, we tested F. tularensis strain LVS and its ΔfslA and ΔfslC mutants for the ability to incorporate potential precursors into rhizoferrin. Unlike putrescine supplementation, supplementation with ornithine greatly enhanced siderophore production by LVS. Radioactivity from L-[U-¹⁴C] ornithine, but not from L-[1-¹⁴C] ornithine, was efficiently incorporated into rhizoferrin by LVS. Although neither the ΔfslA nor the ΔfslC mutant produced rhizoferrin, a putative siderophore intermediate labeled by both [U-¹⁴C] ornithine and [1-¹⁴C] ornithine was secreted by the ΔfslC mutant. Rhizoferrin was identified by liquid chromatography and mass spectrometry in LVS culture supernatants, while citryl-ornithine was detected as the siderophore intermediate in the culture supernatant of the ΔfslC mutant. Our findings support a three-step pathway for rhizoferrin production in Francisella; unlike the fungus Rhizopus delemar, where putrescine functions as a primary precursor for rhizoferrin, biosynthesis in Francisella preferentially starts with ornithine as the substrate for FslA-mediated condensation with citrate. Decarboxylation of this citryl ornithine intermediate by FslC is necessary for a second condensation reaction with citrate to produce rhizoferrin.

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.

Two-Component Systems in Francisella Species

Bacteria alter gene expression in response to changes in their environment through various mechanisms that include signal transduction systems. These signal transduction systems use membrane histidine kinase with sensing domains to mediate phosphotransfer to DNA-binding proteins that alter the level of gene expression. Such regulators are called two-component systems (TCSs). TCSs integrate external signals and information from stress pathways, central metabolism and other global regulators, thus playing an important role as part of the overall regulatory network. This review will focus on the knowledge of TCSs in the Gram-negative bacterium, Francisella tularensis, a biothreat agent with a wide range of potential hosts and a significant ability to cause disease. While TCSs have been well-studied in several bacterial pathogens, they have not been well-studied in non-model organisms, such as F. tularensis and its subspecies, whose canonical TCS content surprisingly ranges from few to none. Additionally, of those TCS genes present, many are orphan components, including KdpDE, QseC, QseB/PmrA, and an unnamed two-component system (FTN_1452/FTN_1453). We discuss recent advances in this field related to the role of TCSs in Francisella physiology and pathogenesis and compare the TCS genes present in human virulent versus. environmental species and subspecies of Francisella.

Human contamination in bacterial genomes has created thousands of spurious proteins

Contaminant sequences that appear in published genomes can cause numerous problems for downstream analyses, particularly for evolutionary studies and metagenomics projects. Our large-scale scan of complete and draft bacterial and archaeal genomes in the NCBI RefSeq database reveals that 2250 genomes are contaminated by human sequence. The contaminant sequences derive primarily from high-copy human repeat regions, which themselves are not adequately represented in the current human reference genome, GRCh38. The absence of the sequences from the human assembly offers a likely explanation for their presence in bacterial assemblies. In some cases, the contaminating contigs have been erroneously annotated as containing protein-coding sequences, which over time have propagated to create spurious protein “families” across multiple prokaryotic and eukaryotic genomes. As a result, 3437 spurious protein entries are currently present in the widely used nr and TrEMBL protein databases. We report here an extensive list of contaminant sequences in bacterial genome assemblies and the proteins associated with them. We found that nearly all contaminants occurred in small contigs in draft genomes, which suggests that filtering out small contigs from draft genome assemblies may mitigate the issue of contamination while still keeping nearly all of the genuine genomic sequences.

BDBM 1.0: A Desktop Application for Efficient Retrieval and Processing of High-Quality Sequence Data and Application to the Identification of the Putative Coffea S-Locus

Nowadays, bioinformatics is one of the most important areas in modern biology and the creation of high-quality scientific software supporting this recent research area is one of the core activities of many researchers. In this context, high-quality sequence datasets are needed to perform inferences on the evolution of species, genes, and gene families, or to get evidence for adaptive amino acid evolution, among others. Nevertheless, sequence data are very often spread over several databases, many useful genomes and transcriptomes are non-annotated, the available annotation is not for the desired coding sequence isoform, and/or is unlikely to be accurate. Moreover, although the FASTA text-based format is quite simple and usable by most software applications, there are a number of issues that may be critical depending on the software used to analyse such files. Therefore, researchers without training in informatics often use a fraction of all available data. The above issues can be addressed using already available software applications, but there is no easy-to-use single piece of software that allows performing all these tasks within the same graphical interface, such as the one here presented, named BDBM (Blast DataBase Manager). BDBM can be used to efficiently get gene sequences from annotated and non-annotated genomes and transcriptomes. Moreover, it can be used to look for alternatives to existing annotations and to easily create reliable custom databases. Such databases are essential to prepare high-quality datasets. The analyses that we have performed on the Coffea canephora genome using BDBM aimed at the identification of the S-locus region (that harbours the genes involved in gametophytic self-incompatibility) led to the conclusion that there are two likely regions, one on chromosome 2 (around region 6600000-6650000), and another on chromosome 5 (around 15830000-15930000). Such findings are discussed in the context of the Rubiaceae gametophytic self-incompatibility evolution.

Complete genome sequencing of sixteen Francisella noatunensis subsp. orientalis isolates: A genomic approach for molecular characterization and spread dynamics of this clonal population

Francisella noatunensis subsp. orientalis (FNO) is an important emerging pathogen associated with disease outbreaks in farm-raised Nile tilapia. FNO genetic diversity using PCR-based typing, no intra-species discrimination was achieved among isolates/strains from different countries, thus demonstrating a clonal behaviour pattern. In this study, we aimed to evaluate the population structure of FNO isolates by comparing whole-genome sequencing data. The analysis of recombination showed that Brazilian isolates group formed a clonal population; whereas other lineages are also supported by this analysis for isolates from foreign countries. The whole-genome multilocus sequence typing (wgMLST) analysis showed varying numbers of dissimilar alleles, suggesting that the Brazilian clonal population are in expansion. Each Brazilian isolate could be identified as a single node by high-resolution gene-by-gene approach, presenting slight genetic differences associated to mutational events. The common ancestry node suggests a single entry into the country before 2012, and the rapid dissemination of this infectious agent may be linked to market sales of infected fingerlings.

Sequence comparison of Francisella tularensis LVS, LVS-G and LVS-R

Francisella tularensis is a gram-negative organism found in many regions of the world. F. tularensis can cause a fatal, febrile illness, although these natural tularemia infections are rare in the United States. However, the development of F. tularensis as a potential weapon of bioterrorism during the Cold War spurred the development of a live attenuated vaccine, LVS, from F. tularensis subsp. holarctica in the 1960s. Two colony morphology variants, LVS-G and LVS-R, were generated from parental LVS by plate passage and by acridine orange mutagenesis, respectively. In vaccinated mice, LVS-G and LVS-R exhibit altered immunogenicity and protective capacities. While the exact nature of the mutations in these strains was unknown, previous studies indicated that both had altered lipopolysaccharide structures. To better understand the impact of these mutations on LVS' immunogenicity, we sequenced the genomes of LVS-G and LVS-R as well as our parental laboratory stock of LVS, originally obtained from ATCC, and compared these to the F. tularensis subsp. holarctica LVS genome currently deposited in GenBank. The results indicate that the genomic sequence of ATCC LVS is nearly identical to that of the human LVS vaccine. Furthermore, a limited number of genomic mutations likely account for the phenotypes of LVS-G and LVS-R.

Galleria mellonella Reveals Niche Differences Between Highly Pathogenic and Closely Related Strains of Francisella spp

Francisella tularensis, a highly virulent bacteria that causes the zoonotic disease tularemia, is considered a potential agent of biological warfare and bioterrorism. Although the host range for several species within the Francisella is known, little is known about the natural reservoirs of various Francisella species. The lack of knowledge regarding the environmental fates of these pathogens greatly reduces the possibilities for microbial risk assessments. The greater wax moth (Galleria mellonella) is an insect of the order Lepidoptera that has been used as an alternative model to study microbial infection during recent years. The aim of this study was to evaluate G. mellonella as a model system for studies of human pathogenic and closely related opportunistic and non-pathogenic strains within the Francisella genus. The employed G. mellonella larvae model demonstrated differences in lethality between human pathogenic and human non-pathogenic or opportunistic Francisella species. The F. novicida, F. hispaniensis and F. philomiragia strains were significantly more virulent in the G. mellonella model than the strains of human pathogens F. t. holarctica and F. t. tularensis. Our data show that G. mellonella is a possible in vivo model of insect immunity for studies of both opportunistic and virulent lineages of Francisella spp., that produces inverse results regarding lethality in G. mellonella and incapacitating disease in humans. The results provide insight into the potential host specificity of F. tularensis and closely related members of the same genus, thus increasing our present understanding of Francisella spp. ecology.

Live Attenuated Tularemia Vaccines for Protection Against Respiratory Challenge With Virulent F. tularensis subsp. tularensis

Francisella tularensis is the causative agent of tularemia and a Tier I bioterrorism agent. In the 1900s, several vaccines were developed against tularemia including the killed "Foshay" vaccine, subunit vaccines comprising F. tularensis protein(s) or lipoproteins(s) in an adjuvant formulation, and the F. tularensis Live Vaccine Strain (LVS); none were licensed in the U.S.A. or European Union. The LVS vaccine retains toxicity in humans and animals-especially mice-but has demonstrated efficacy in humans, and thus serves as the current gold standard for vaccine efficacy studies. The U.S.A. 2001 anthrax bioterrorism attack spawned renewed interest in vaccines against potential biowarfare agents including F. tularensis. Since live attenuated-but not killed or subunit-vaccines have shown promising efficacy and since vaccine efficacy against respiratory challenge with less virulent subspecies holarctica or F. novicida, or against non-respiratory challenge with virulent subsp. tularensis (Type A) does not reliably predict vaccine efficacy against respiratory challenge with virulent subsp. tularensis, the route of transmission and species of greatest concern in a bioterrorist attack, in this review, we focus on live attenuated tularemia vaccine candidates tested against respiratory challenge with virulent Type A strains, including homologous vaccines derived from mutants of subsp. holarctica, F. novicida, and subsp. tularensis, and heterologous vaccines developed using viral or bacterial vectors to express F. tularensis immunoprotective antigens. We compare the virulence and efficacy of these vaccine candidates with that of LVS and discuss factors that can significantly impact the development and evaluation of live attenuated tularemia vaccines. Several vaccines meet what we would consider the minimum criteria for vaccines to go forward into clinical development-safety greater than LVS and efficacy at least as great as LVS, and of these, several meet the higher standard of having efficacy ≥LVS in the demanding mouse model of tularemia. These latter include LVS with deletions in purMCD, sodBFt , capB or wzy; LVS ΔcapB that also overexpresses Type VI Secretion System (T6SS) proteins; FSC200 with a deletion in clpB; the single deletional purMCD mutant of F. tularensis SCHU S4, and a heterologous prime-boost vaccine comprising LVS ΔcapB and Listeria monocytogenes expressing T6SS proteins.

The Ability to Acquire Iron Is Inversely Related to Virulence and the Protective Efficacy of Francisella tularensis Live Vaccine Strain

Francisella tularensis is a highly infectious bacterial pathogen that causes the potentially fatal disease tularemia. The Live Vaccine Strain (LVS) of F. tularensis subsp. holarctica, while no longer licensed as a vaccine, is used as a model organism for identifying correlates of immunity and bacterial factors that mediate a productive immune response against F. tularensis. Recently, it was reported that two biovars of LVS differed in their virulence and vaccine efficacy. Genetic analysis showed that they differ in ferrous iron homeostasis; lower Fe2+ levels contributed to increased resistance to hydrogen peroxide in the vaccine efficacious LVS biovar. This also correlated with resistance to the bactericidal activity of interferon γ-stimulated murine bone marrow-derived macrophages. We have extended these findings further by showing that a mutant lacking bacterioferritin stimulates poor protection against Schu S4 challenge in a mouse model of tularemia. Together these results suggest that the efficacious biovar of LVS stimulates productive immunity by a mechanism that is dependent on its ability to limit the toxic effects of oxidative stress by maintaining optimally low levels of intracellular Fe2+.

Revisiting Francisella tularensis subsp. holarctica, Causative Agent of Tularemia in Germany With Bioinformatics: New Insights in Genome Structure, DNA Methylation and Comparative Phylogenetic Analysis

Francisella (F.) tularensis is a highly virulent, Gram-negative bacterial pathogen and the causative agent of the zoonotic disease tularemia. Here, we generated, analyzed and characterized a high quality circular genome sequence of the F. tularensis subsp. holarctica strain 12T0050 that caused fatal tularemia in a hare. Besides the genomic structure, we focused on the analysis of oriC, unique to the Francisella genus and regulating replication in and outside hosts and the first report on genomic DNA methylation of a Francisella strain. The high quality genome was used to establish and evaluate a diagnostic whole genome sequencing pipeline. A genotyping strategy for F. tularensis was developed using various bioinformatics tools for genotyping. Additionally, whole genome sequences of F. tularensis subsp. holarctica isolates isolated in the years 2008–2015 in Germany were generated. A phylogenetic analysis allowed to determine the genetic relatedness of these isolates and confirmed the highly conserved nature of F. tularensis subsp. holarctica.

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.

An integrative strategy to identify the entire protein coding potential of prokaryotic genomes by proteogenomics

Accurate annotation of all protein-coding sequences (CDSs) is an essential prerequisite to fully exploit the rapidly growing repertoire of completely sequenced prokaryotic genomes. However, large discrepancies among the number of CDSs annotated by different resources, missed functional short open reading frames (sORFs), and overprediction of spurious ORFs represent serious limitations. Our strategy toward accurate and complete genome annotation consolidates CDSs from multiple reference annotation resources, ab initio gene prediction algorithms and in silico ORFs (a modified six-frame translation considering alternative start codons) in an integrated proteogenomics database (iPtgxDB) that covers the entire protein-coding potential of a prokaryotic genome. By extending the PeptideClassifier concept of unambiguous peptides for prokaryotes, close to 95% of the identifiable peptides imply one distinct protein, largely simplifying downstream analysis. Searching a comprehensive Bartonella henselae proteomics data set against such an iPtgxDB allowed us to unambiguously identify novel ORFs uniquely predicted by each resource, including lipoproteins, differentially expressed and membrane-localized proteins, novel start sites and wrongly annotated pseudogenes. Most novelties were confirmed by targeted, parallel reaction monitoring mass spectrometry, including unique ORFs and single amino acid variations (SAAVs) identified in a re-sequenced laboratory strain that are not present in its reference genome. We demonstrate the general applicability of our strategy for genomes with varying GC content and distinct taxonomic origin. We release iPtgxDBs for B. henselae, Bradyrhizobium diazoefficiens and Escherichia coli and the software to generate both proteogenomics search databases and integrated annotation files that can be viewed in a genome browser for any prokaryote.

Intratracheal Inoculation of Fischer 344 Rats with Francisella tularensis

Pulmonary infection with the bacterium Francisella tularensis can lead to the serious and potentially fatal disease, tularemia, in humans. Due to the current lack of an approved tularemia vaccine for humans, research is focused on vaccine development utilizing appropriate animal models. The Fischer 344 rat has emerged as a model that reflects human susceptibility to F. tularensis infection, and thus is an attractive model for tularemia vaccine development. Intratracheal inoculation of the Fischer 344 rat with F. tularensis mimics pulmonary exposure in humans. The successful delivery into the rat trachea is critical for pulmonary delivery. A laryngoscope with illumination is used to properly intubate the tracheae of anesthetized rats; the correct placement within the trachea is determined by a simple device to detect breathing. Following intubation, the F. tularensis culture is delivered in a measured dose via syringe. This technique standardizes pulmonary delivery of F. tularensis within the rat trachea to evaluate vaccine efficacy.

Genome Plasticity and Polymorphisms in Critical Genes Correlate with Increased Virulence of Dutch Outbreak-Related Coxiella burnetii Strains

Coxiella burnetii is an obligate intracellular bacterium and the etiological agent of Q fever. During 2007-2010 the largest Q fever outbreak ever reported occurred in The Netherlands. It is anticipated that strains from this outbreak demonstrated an increased zoonotic potential as more than 40,000 individuals were assumed to be infected. The acquisition of novel genetic factors by these C. burnetii outbreak strains, such as virulence-related genes, has frequently been proposed and discussed, but is not proved yet. In the present study, the whole genome sequence of several Dutch strains (CbNL01 and CbNL12 genotypes), a few additionally selected strains from different geographical locations and publicly available genome sequences were used for a comparative bioinformatics approach. The study focuses on the identification of specific genetic differences in the outbreak related CbNL01 strains compared to other C. burnetii strains. In this approach we investigated the phylogenetic relationship and genomic aspects of virulence and host-specificity. Phylogenetic clustering of whole genome sequences showed a genotype-specific clustering that correlated with the clustering observed using Multiple Locus Variable-number Tandem Repeat Analysis (MLVA). Ortholog analysis on predicted genes and single nucleotide polymorphism (SNP) analysis of complete genome sequences demonstrated the presence of genotype-specific gene contents and SNP variations in C. burnetii strains. It also demonstrated that the currently used MLVA genotyping methods are highly discriminatory for the investigated outbreak strains. In the fully reconstructed genome sequence of the Dutch outbreak NL3262 strain of the CbNL01 genotype, a relatively large number of transposon-linked genes were identified as compared to the other published complete genome sequences of C. burnetii. Additionally, large numbers of SNPs in its membrane proteins and predicted virulence-associated genes were identified in all Dutch outbreak strains compared to the NM reference strain and other strains of the CbNL12 genotype. The presence of large numbers of transposable elements and mutated genes, thereof most likely resulted in high level of genome rearrangements and genotype-specific pathogenicity of outbreak strains. Thus, the epidemic potential of Dutch outbreak strains could be linked to increased genome plasticity and mutations in critical genes involved in virulence and the evasion of the host immune system.

Clonal structure in Ichthyobacterium seriolicida, the causative agent of bacterial haemolytic jaundice in yellowtail, Seriola quinqueradiata, inferred from molecular epidemiological analysis

Bacterial haemolytic jaundice caused by Ichthyobacterium seriolicida has been responsible for mortality in farmed yellowtail, Seriola quinqueradiata, in western Japan since the 1980s. In this study, polymorphic analysis of I. seriolicida was performed using three molecular methods: amplified fragment length polymorphism (AFLP) analysis, multilocus sequence typing (MLST) and multiple-locus variable-number tandem repeat analysis (MLVA). Twenty-eight isolates were analysed using AFLP, while 31 isolates were examined by MLST and MLVA. No polymorphisms were identified by AFLP analysis using EcoRI and MseI, or by MLST of internal fragments of eight housekeeping genes. However, MLVA revealed variation in repeat numbers of three elements, allowing separation of the isolates into 16 sequence types. The unweighted pair group method using arithmetic averages cluster analysis of the MLVA data identified four major clusters, and all isolates belonged to clonal complexes. It is likely that I. seriolicida populations share a common ancestor, which may be a recently introduced strain.

Multi-omics Analysis Sheds Light on the Evolution and the Intracellular Lifestyle Strategies of Spotted Fever Group Rickettsia spp

Arthropod-borne Rickettsia species are obligate intracellular bacteria which are pathogenic for humans. Within this genus, Rickettsia slovaca and Rickettsia conorii cause frequent and potentially severe infections, whereas Rickettsia raoultii and Rickettsia massiliae cause rare and milder infections. All four species belong to spotted fever group (SFG) rickettsiae. However, R. slovaca and R. raoultii cause scalp eschar and neck lymphadenopathy (SENLAT) and are mainly associated with Dermacentor ticks, whereas the other two species cause Mediterranean spotted fever (MSF) and are mainly transmitted by Rhipicephalus ticks. To identify the potential genes and protein profiles and to understand the evolutionary processes that could, comprehensively, relate to the differences in virulence and pathogenicity observed between these four species, we compared their genomes and proteomes. The virulent and milder agents displayed divergent phylogenomic evolution in two major clades, whereas either SENLAT or MSF disease suggests a discrete convergent evolution of one virulent and one milder agent, despite their distant genetic relatedness. Moreover, the two virulent species underwent strong reductive genomic evolution and protein structural variations, as well as a probable loss of plasmid(s), compared to the two milder species. However, an abundance of mobilome genes was observed only in the less pathogenic species. After infecting Xenopus laevis cells, the virulent agents displayed less up-regulated than down-regulated proteins, as well as less number of identified core proteins. Furthermore, their similar and distinct protein profiles did not contain some genes (e.g., ompA/B and rickA) known to be related to rickettsial adhesion, motility and/or virulence, but may include other putative virulence-, antivirulence-, and/or disease-related proteins. The identified evolutionary forces herein may have a strong impact on intracellular expressions and strategies in these rickettsiae, and that may contribute to the emergence of distinct virulence and diseases in humans. Thus, the current multi-omics data provide new insights into the evolution and fitness of SFG virulence and pathogenicity, and intracellular pathogenic bacteria.

Disentangling the effects of selection and loss bias on gene dynamics

We combine mathematical modeling of genome evolution with comparative analysis of prokaryotic genomes to estimate the relative contributions of selection and intrinsic loss bias to the evolution of different functional classes of genes and mobile genetic elements (MGE). An exact solution for the dynamics of gene family size was obtained under a linear duplication-transfer-loss model with selection. With the exception of genes involved in information processing, particularly translation, which are maintained by strong selection, the average selection coefficient for most nonparasitic genes is low albeit positive, compatible with observed positive correlation between genome size and effective population size. Free-living microbes evolve under stronger selection for gene retention than parasites. Different classes of MGE show a broad range of fitness effects, from the nearly neutral transposons to prophages, which are actively eliminated by selection. Genes involved in antiparasite defense, on average, incur a fitness cost to the host that is at least as high as the cost of plasmids. This cost is probably due to the adverse effects of autoimmunity and curtailment of horizontal gene transfer caused by the defense systems and selfish behavior of some of these systems, such as toxin-antitoxin and restriction modification modules. Transposons follow a biphasic dynamics, with bursts of gene proliferation followed by decay in the copy number that is quantitatively captured by the model. The horizontal gene transfer to loss ratio, but not duplication to loss ratio, correlates with genome size, potentially explaining increased abundance of neutral and costly elements in larger genomes.

The Legionella pneumophila GIG operon responds to gold and copper in planktonic and biofilm cultures

Legionella pneumophila contaminates man-made water systems and creates numerous exposure risks for Legionnaires’ Disease. Because copper/silver ionization is commonly used to control L. pneumophila, its mechanisms of metal response and detoxification are of significant interest. Here we describe an L. pneumophila operon with significant similarity to the GIG operon of Cupriavidus metallidurans. The Legionella GIG operon is present in a subset of strains and has been acquired as part of the ICE-βox 65-kB integrative conjugative element. We assessed GIG promoter activity following exposure of L. pneumophila to multiple concentrations of HAuCl₄, CuSO_{4 and} AgNO₃. At 37°C, control stationary phase cultures exhibited GIG promoter activity. This activity increased significantly in response to 20 and 50uM HAuCl₄ and CuSO₄ but not in response to AgNO₃. Conversely, at 26°C, cultures exhibited decreased promoter response to copper. GIG promoter activity was also induced by HAuCl₄ or CuSO₄ during early biofilm establishment at both temperatures. When an L. pneumophila GIG promoter construct was transformed into E. coli DH5α, cultures showed baseline expression levels that did not increase following metal addition. Analysis of L. pneumophila transcriptional regulatory mutants suggested that GIG up-regulation in the presence of metal ions may be influenced by the stationary phase sigma factor, RpoS.

The impact of insertion sequences on bacterial genome plasticity and adaptability

Transposable elements (TE), small mobile genetic elements unable to exist independently of the host genome, were initially believed to be exclusively deleterious genomic parasites. However, it is now clear that they play an important role as bacterial mutagenic agents, enabling the host to adapt to new environmental challenges and to colonize new niches. This review focuses on the impact of insertion sequences (IS), arguably the smallest TE, on bacterial genome plasticity and concomitant adaptability of phenotypic traits, including resistance to antibacterial agents, virulence, pathogenicity and catabolism. The direct consequence of IS transposition is the insertion of one DNA sequence into another. This event can result in gene inactivation as well as in modulation of neighbouring gene expression. The latter is usually mediated by de-repression or by the introduction of a complete or partial promoter located within the element. Furthermore, transcription and transposition of IS are affected by host factors and in some cases by environmental signals offering the host an adaptive strategy and promoting genetic variability to withstand the environmental challenges.

Iron and Virulence in Francisella tularensis

Francisella tularensis, the causative agent of tularemia, is a Gram-negative bacterium that infects a variety of cell types including macrophages, and propagates with great efficiency in the cytoplasm. Iron, essential for key enzymatic and redox reactions, is among the nutrients required to support this pathogenic lifestyle and the bacterium relies on specialized mechanisms to acquire iron within the host environment. Two distinct pathways for iron acquisition are encoded by the F. tularensis genome- a siderophore-dependent ferric iron uptake system and a ferrous iron transport system. Genes of the Fur-regulated fslABCDEF operon direct the production and transport of the siderophore rhizoferrin. Siderophore biosynthesis involves enzymes FslA and FslC, while export across the inner membrane is mediated by FslB. Uptake of the rhizoferrin- ferric iron complex is effected by the siderophore receptor FslE in the outer membrane in a TonB-independent process, and FslD is responsible for uptake across the inner membrane. Ferrous iron uptake relies largely on high affinity transport by FupA in the outer membrane, while the Fur-regulated FeoB protein mediates transport across the inner membrane. FslE and FupA are paralogous proteins, sharing sequence similarity and possibly sharing structural features as well. This review summarizes current knowledge of iron acquisition in this organism and the critical role of these uptake systems in bacterial pathogenicity.

A spontaneous mutation in kdsD, a biosynthesis gene for 3 Deoxy-D-manno-Octulosonic Acid, occurred in a ciprofloxacin resistant strain of Francisella tularensis and caused a high level of attenuation in murine models of tularemia

Francisella tularensis, a gram-negative facultative intracellular bacterial pathogen, is the causative agent of tularemia and able to infect many mammalian species, including humans. Because of its ability to cause a lethal infection, low infectious dose, and aerosolizable nature, F. tularensis subspecies tularensis is considered a potential biowarfare agent. Due to its in vitro efficacy, ciprofloxacin is one of the antibiotics recommended for post-exposure prophylaxis of tularemia. In order to identify therapeutics that will be efficacious against infections caused by drug resistant select-agents and to better understand the threat, we sought to characterize an existing ciprofloxacin resistant (CipR) mutant in the Schu S4 strain of F. tularensis by determining its phenotypic characteristics and sequencing the chromosome to identify additional genetic alterations that may have occurred during the selection process. In addition to the previously described genetic alterations, the sequence of the CipR mutant strain revealed several additional mutations. Of particular interest was a frameshift mutation within kdsD which encodes for an enzyme necessary for the production of 3-Deoxy-D-manno-Octulosonic Acid (KDO), an integral component of the lipopolysaccharide (LPS). A kdsD mutant was constructed in the Schu S4 strain. Although it was not resistant to ciprofloxacin, the kdsD mutant shared many phenotypic characteristics with the CipR mutant, including growth defects under different conditions, sensitivity to hydrophobic agents, altered LPS profiles, and attenuation in multiple models of murine tularemia. This study demonstrates that the KdsD enzyme is essential for Francisella virulence and may be an attractive therapeutic target for developing novel medical countermeasures.

Enhanced survival but not amplification of Francisella spp. in the presence of free-living amoebae

Transmission of Francisella tularensis, the etiologic agent of tularemia, has been associated with various water sources. Survival of many waterborne pathogens within free-living amoeba (FLA) is well documented; however, the role of amoebae in the environmental persistence of F. tularensis is unclear. In this study, axenic FLA cultures of Acanthamoeba castellanii, Acanthamoeba polyphaga, and Vermamoeba vermiformis were each inoculated with virulent strains of F. tularensis (Types A and B), the attenuated live vaccine strain, and Francisella novicida. Experimental parameters included low and high multiplicity of infection and incubation temperatures of 25 and 30 °C for 0-10 days. Francisella spp. survival was enhanced by the presence of FLA; however, bacterial growth and protozoa infectivity were not observed. In contrast, co-infections of A. polyphaga and Legionella pneumophila, used as an amoeba pathogen control, resulted in bacterial proliferation, cytopathic effects, and amoebal lysis. Collectively, even though short-term incubation with FLA was beneficial, the long-term effects on Francisella survival are unknown, especially given the expenditure of available amoebal derived nutrients and the fastidious nature of Francisella spp. These factors have clear implications for the role of FLA in Francisella environmental persistence.

Genomic changes associated with the evolutionary transition of an insect gut symbiont into a blood-borne pathogen

The genus Bartonella comprises facultative intracellular bacteria with a unique lifestyle. After transmission by blood-sucking arthropods they colonize the erythrocytes of mammalian hosts causing acute and chronic infectious diseases. Although the pathogen–host interaction is well understood, little is known about the evolutionary origin of the infection strategy manifested by Bartonella species. Here we analyzed six genomes of Bartonella apis, a honey bee gut symbiont that to date represents the closest relative of pathogenic Bartonella species. Comparative genomics revealed that B. apis encodes a large set of vertically inherited genes for amino acid and cofactor biosynthesis and nitrogen metabolism. Most pathogenic bartonellae have lost these ancestral functions, but acquired specific virulence factors and expanded a vertically inherited gene family for harvesting cofactors from the blood. However, the deeply rooted pathogen Bartonella tamiae has retained many of the ancestral genome characteristics reflecting an evolutionary intermediate state toward a host-restricted intraerythrocytic lifestyle. Our findings suggest that the ancestor of the pathogen Bartonella was a gut symbiont of insects and that the adaptation to blood-feeding insects facilitated colonization of the mammalian bloodstream. This study highlights the importance of comparative genomics among pathogens and non-pathogenic relatives to understand disease emergence within an evolutionary-ecological framework.

Whole-Genome Relationships among Francisella Bacteria of Diverse Origins Define New Species and Provide Specific Regions for Detection

Francisella tularensis is a highly virulent zoonotic pathogen that causes tularemia and, because of weaponization efforts in past world wars, is considered a tier 1 biothreat agent. Detection and surveillance of F. tularensis may be confounded by the presence of uncharacterized, closely related organisms. Through DNA-based diagnostics and environmental surveys, novel clinical and environmental Francisella isolates have been obtained in recent years. Here we present 7 new Francisella genomes and a comparison of their characteristics to each other and to 24 publicly available genomes as well as a comparative analysis of 16S rRNA and sdhA genes from over 90 Francisella strains. Delineation of new species in bacteria is challenging, especially when isolates having very close genomic characteristics exhibit different physiological features-for example, when some are virulent pathogens in humans and animals while others are nonpathogenic or are opportunistic pathogens. Species resolution within Francisella varies with analyses of single genes, multiple gene or protein sets, or whole-genome comparisons of nucleic acid and amino acid sequences. Analyses focusing on single genes (16S rRNA, sdhA), multiple gene sets (virulence genes, lipopolysaccharide [LPS] biosynthesis genes, pathogenicity island), and whole-genome comparisons (nucleotide and protein) gave congruent results, but with different levels of discrimination confidence. We designate four new species within the genus; Francisella opportunistica sp. nov. (MA06-7296), Francisella salina sp. nov. (TX07-7308), Francisella uliginis sp. nov. (TX07-7310), and Francisella frigiditurris sp. nov. (CA97-1460). This study provides a robust comparative framework to discern species and virulence features of newly detected Francisella bacteria.

IMPORTANCE

DNA-based detection and sequencing methods have identified thousands of new bacteria in the human body and the environment. In most cases, there are no cultured isolates that correspond to these sequences. While DNA-based approaches are highly sensitive, accurately assigning species is difficult without known near relatives for comparison. This ambiguity poses challenges for clinical cases, disease epidemics, and environmental surveillance, for which response times must be short. Many new Francisella isolates have been identified globally. However, their species designations and potential for causing human disease remain ambiguous. Through detailed genome comparisons, we identified features that differentiate F. tularensis from clinical and environmental Francisella isolates and provide a knowledge base for future comparison of Francisella organisms identified in clinical samples or environmental surveys.