Molecular evolutionary consequences of niche restriction in Francisella tularensis, a facultative intracellular pathogen

Genetic Homogeneity of Francisella tularensis subsp. mediasiatica Strains in Kazakhstan

Tularemia is an acute febrile disease caused by the Gram-negative bacillus Francisella tularensis. Based on genetic and phenotypic characteristics, three subspecies are distinguished: tularensis, holarctica, and mediasiatica. F. tularensis subsp. mediasiatica remains the least studied subspecies. Over the past decade, new foci of distribution of F. tularensis subsp. mediasiatica have been discovered in Russia (Siberia), expanding the possible distribution area by thousands of kilometers. This article provides whole genome single nucleotide polymorphism (wgSNP) and polymorphic tandem repeats (MLVA) analyses of 28 mediasiatica strains isolated between 1965 and 2004 in Kazakhstan. Despite high genetic homogeneity, MLVA with eleven loci (MLVA11) demonstrates a high discriminatory ability (diversity index, 0.9497). The topological structure of the trees based on wgSNP and MLVA is not comparable; however, clustering remains congruent for most outbreaks, with the exception of two strains from one outbreak that are identical in terms of wgSNP but differ at three tandem repeat loci. Based on wgSNP, the strains are assigned to one of the three currently known mediasiatica sublineages, lineage M.I, together with other historical strains maintained in collections in Russia and Sweden. wgSNP shows limited previously unknown genetic diversity, with the M.I lineage size being only 118 SNPs. The wgSNP genotype is not strongly correlated with year and place of isolation.

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.

Identification of Francisella tularensis Subspecies in a Clinical Setting Using MALDI-TOF MS: An In-House Francisella Library and Biomarkers

Francisella tularensis is a zoonotic bacterium that is endemic in large parts of the world. It is absent in the standard library of the most applied matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) systems: the Vitek MS and the Bruker Biotyper system. The additional Bruker MALDI Biotyper Security library contains F. tularensis without subspecies differentiation. The virulence of F. tularensis differs between the subspecies. The F. tularensis subspecies (ssp.) tularensis is highly pathogenic, whereas the subspecies holarctica displays lower virulence and subspecies novicida and F. tularensis ssp. mediasiatica are hardly virulent. To differentiate the Francisellaceae and the F. tularensis-subspecies, an in-house Francisella library was built with the Bruker Biotyper system and validated together with the existing Bruker databases. In addition, specific biomarkers were defined based on the main spectra of the Francisella strains supplemented with in silico genome data. Our in-house Francisella library accurately differentiates the F. tularensis subspecies and the other Francisellaceae. The biomarkers correctly differentiate the various species within the genus Francisella and the F. tularensis subspecies. These MALDI-TOF MS strategies can successfully be applied in a clinical laboratory setting as a fast and specific method to identify F. tularensis to subspecies level.

Development and evaluation of a multiplex droplet digital polymerase chain reaction method for simultaneous detection of five biothreat pathogens

Biothreat agents pose a huge threat to human and public health, necessitating the development of rapid and highly sensitive detection approaches. This study establishes a multiplex droplet digital polymerase chain reaction (ddPCR) method for simultaneously detecting five high-risk bacterial biothreats: Yersinia pestis, Bacillus anthracis, Brucella spp., Burkholderia pseudomallei, and Francisella tularensis. Unlike conventional multiplex real-time PCR (qPCR) methods, the multiplex ddPCR assay was developed using two types of probe fluorophores, allowing the assay to perform with a common two-color ddPCR system. After optimization, the assay performance was evaluated, showing a lower limit of detection (LOD) (0.1–1.0 pg/μL) and good selectivity for the five bacteria targets. The multiplex assay’s ability to simultaneously detect two or more kinds of targets in a sample was also demonstrated. The assay showed strong sample tolerance when testing simulated soil samples; the LOD for bacteria in soil was 2 × 10²–2 × 10³ colony-forming unit (CFU)/100 mg soil (around 5–50 CFU/reaction), which was 10-fold lower than that of the single-target qPCR method. When testing simulated soil samples at bacterial concentrations of 2 × 10³–2 × 10⁴ CFU/100 mg soil, the assay presented a higher sensitivity (100%, 35/35) than that of the qPCR method (65.71%, 23/35) and a good specificity (100%, 15/15). These results suggest that the developed 5-plex ddPCR method is more sensitive than conventional qPCR methods and is potentially suitable for rapidly detecting or screening the five selected bacterial biothreats in suspicious samples.

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.

Multilocus sequence analysis reveals different lineages of Pseudomonas anguilliseptica associated with disease in farmed lumpfish (Cyclopterus lumpus L.)

The bacterium Pseudomonas anguilliseptica has in recent years emerged as a serious threat to production of lumpfish in Norway. Little is known about the population structure of this bacterium despite its association with disease in a wide range of different fish species throughout the world. The phylogenetic relationships between 53 isolates, primarily derived from diseased lumpfish, but including a number of reference strains from diverse geographical origins and fish species, were reconstructed by Multi-Locus Sequence Analysis (MLSA) using nine housekeeping genes (rpoB, atpD, gyrB, rpoD, ileS, aroE, carA, glnS and recA). MLSA revealed a high degree of relatedness between the studied isolates, altough the seven genotypes identified formed three main phylogenetic lineages. While four genotypes were identified amongst Norwegian lumpfish isolates, a single genotype dominated, irrespective of geographic origin. This suggests the existence of a dominant genotype associated with disease in production of lumpfish in Norwegian aquaculture. Elucidation of the population structure of the bacterium has provided valuable information for potential future vaccine development.

Development, Phenotypic Characterization and Genomic Analysis of a Francisella tularensis Panel for Tularemia Vaccine Testing

Francisella tularensis is one of several biothreat agents for which a licensed vaccine is needed to protect against this pathogen. To aid in the development of a vaccine protective against pneumonic tularemia, we generated and characterized a panel of F. tularensis isolates that can be used as challenge strains to assess vaccine efficacy. Our panel consists of both historical and contemporary isolates derived from clinical and environmental sources, including human, tick, and rabbit isolates. Whole genome sequencing was performed to assess the genetic diversity in comparison to the reference genome F. tularensis Schu S4. Average nucleotide identity analysis showed >99% genomic similarity across the strains in our panel, and pan-genome analysis revealed a core genome of 1,707 genes, and an accessory genome of 233 genes. Three of the strains in our panel, FRAN254 (tick-derived), FRAN255 (a type B strain), and FRAN256 (a human isolate) exhibited variation from the other strains. Moreover, we identified several unique mutations within the Francisella Pathogenicity Island across multiple strains in our panel, revealing unexpected diversity in this region. Notably, FRAN031 (Scherm) completely lacked the second pathogenicity island but retained virulence in mice. In contrast, FRAN037 (Coll) was attenuated in a murine pneumonic tularemia model and had mutations in pdpB and iglA which likely led to attenuation. All of the strains, except FRAN037, retained full virulence, indicating their effectiveness as challenge strains for future vaccine testing. Overall, we provide a well-characterized panel of virulent F. tularensis strains that can be utilized in ongoing efforts to develop an effective vaccine against pneumonic tularemia to ensure protection is achieved across a range F. tularensis strains.

Long-Term Survival of Virulent Tularemia Pathogens outside a Host in Conditions That Mimic Natural Aquatic Environments

Francisella tularensis, the causative agent of the zoonotic disease tularemia, can cause seasonal outbreaks of acute febrile illness in humans with disease peaks in late summer to autumn. Interestingly, its mechanisms for environmental persistence between outbreaks are poorly understood. One hypothesis is that F. tularensis forms biofilms in aquatic environments. We utilized two fully virulent wild-type strains: FSC200 (Francisella tularensis subsp. holarctica) and Schu S4 (Francisella tularensis subsp. tularensis) and three control strains, the attenuated live vaccine strain (LVS; F. tularensis subsp. holarctica), a Schu S4 ΔwbtI mutant that is documented to form biofilms, and the low-virulence strain U112 of the closely related species Francisella novicida Strains were incubated in saline solution (0.9% NaCl) microcosms for 24 weeks at both 4°C and 20°C, whereupon viability and biofilm formation were measured. These temperatures were selected to approximate winter and summer temperatures of fresh water in Scandinavia, respectively. U112 and Schu S4 ΔwbtI formed biofilms, but F. tularensis strains FSC200 and Schu S4 and the LVS did not. All strains exhibited prolonged viability at 4°C compared to 20°C. U112 and FSC200 displayed remarkable long-term persistence at 4°C, with only 1- and 2-fold log reductions, respectively, of viable cells after 24 weeks. Schu S4 exhibited lower survival, yielding no viable cells by week 20. At 24 weeks, cells from FSC200, but not from Schu S4, were still fully virulent in mice. Taken together, these results demonstrate biofilm-independent, long-term survival of pathogenic F. tularensis subsp. holarctica in conditions that mimic overwinter survival in aquatic environments.IMPORTANCE Tularemia, a disease caused by the environmental bacterium Francisella tularensis, is characterized by acute febrile illness. F. tularensis is highly infectious: as few as 10 organisms can cause human disease. Tularemia is not known to be spread from person to person. Rather, all human infections are independently acquired from the environment via the bite of blood-feeding arthropods, ingestion of infected food or water, or inhalation of aerosolized bacteria. Despite the environmental origins of human disease events, the ecological factors governing the long-term persistence of F. tularensis in nature between seasonal human outbreaks are poorly understood. The significance of our research is in identifying conditions that promote long-term survival of fully virulent F. tularensis outside a mammalian host or insect vector. These conditions are similar to those found in natural aquatic environments in winter and provide important new insights on how F. tularensis may persist long-term in the environment.

Characterization of Tularemia Cases in Slovenia with Multiple-Locus Variable-Number Tandem Repeat Analysis

Francisella tularensis is the etiologic agent of tularemia, a bacterial zoonotic disease. The genome of F. tularensis shows a recent evolutionary change, especially in reservoirs. Variable number of tandem repeats (VNTR) is described as a high-speed molecular clock and can thus be used as a high-resolution typing system. The main objective of our study was to investigate the molecular diversity of F. tularensis strains and reveal possible sources of infection. Using real-time PCR targeting the ISFtu2 region, we successfully amplified targeted DNA in 13/31 Slovenian patients with a clinical diagnosis of tularemia, and with PCR targeting the fopA gene, we obtained 11/13 PCR products. Sequencing revealed that all samples were identified as F. tularensis subsp. holarctica. We successfully obtained one F. tularensis isolate from a lymph node aspirate by culture on chocolate agar. Our isolate was clustered into major clade B12 (subclade B43). We optimized VNTR typing to be used directly on clinical samples. Multiple-locus VNTR analysis (MLVA) revealed five unique MLVA types; 45.5% samples had the same MLVA type, another 27.3% shared a different MLVA type, and each of the remaining had a unique MLVA type. Most samples differed at only two VNTR markers (Ft-M03 and Ft-M06). Additionally, we investigated samples from small mammals (n = 532) and Ixodes ricinus ticks (n = 232) captured in the same geographical area in which patients with tularemia were found. No F. tularensis DNA was detected in samples of small mammals or I. ricinus ticks. The diversity of MLVA types in Slovenia was high, despite the small region, but most of the samples from the same region shared the same MLVA type. Our results suggest that MLVA is a useful tool for quick molecular characterization of F. tularensis directly from patient samples, especially when investigating geographically localized outbreaks.

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.

Evaluation of the Biotoxis qPCR Detection Kit for Francisella tularensis Detection in Clinical and Environmental Samples

Rapid and reliable detection and identification of Francisella tularensis (a tier 1 select agent) are of primary interest for both medical and biological threat surveillance purposes. The Biotoxis qPCR detection kit is a real-time quantitative PCR (qPCR) assay designed for the detection of Bacillus anthracis, Yersinia pestis, and F. tularensis in environmental or biological samples. Here, we evaluated its performance for detecting F. tularensis in comparison to previously validated qPCR assays. The Biotoxis qPCR was positive for 87/87 F. tularensis subsp. holarctica (type B) strains but also for F. tularensis subsp. novicida It was negative for Francisella philomiragia and 24/24 strains belonging to other bacterial species. For 31 tularemia clinical specimens, the Biotoxis qPCR displayed a sensitivity between 90.32% and 96.55%, compared to qPCR tests targeting ISFtu2 or a type B-specific DNA sequence, respectively. All 30 nontularemia clinical specimens were Biotoxis qPCR negative. For water samples, the Biotoxis qPCR limit of detection was 1,000 CFU/liter of F. tularensis For 57 environmental water samples collected in France, the Biotoxis qPCR was positive for 6/15 samples positive for ISFtu2 qPCR and 4/4 positive for type B qPCR. In conclusion, the Biotoxis qPCR detection kit demonstrated good performances for F. tularensis detection in various biological and environmental samples, although cross-amplification of F. tularensis subsp. novicida must be considered. This plate format assay could be useful to test a large number of clinical or environmental specimens, especially in the context of natural or intentional tularemia outbreaks.

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.

Using affinity propagation clustering for identifying bacterial clades and subclades with whole-genome sequences of Francisella tularensis

By combining a reference-independent SNP analysis and average nucleotide identity (ANI) with affinity propagation clustering (APC), we developed a significantly improved methodology allowing resolving phylogenetic relationships, based on objective criteria. These bioinformatics tools can be used as a general ruler to determine phylogenetic relationships and clustering of bacteria, exemplary done with Francisella (F.) tularensis. Molecular epidemiology of F. tularensis is currently assessed mostly based on laboratory methods and molecular analysis. The high evolutionary stability and the clonal nature makes Francisella ideal for subtyping with single nucleotide polymorphisms (SNPs). Sequencing and real-time PCR can be used to validate the SNP analysis. We investigate whole-genome sequences of 155 F. tularensis subsp. holarctica isolates. Phylogenetic testing was based on SNPs and average nucleotide identity (ANI) as reference independent, alignment-free methods taking small-scale and large-scale differences within the genomes into account. Especially the whole genome SNP analysis with kSNP3.0 allowed deciphering quite subtle signals of systematic differences in molecular variation. Affinity propagation clustering (APC) resulted in three clusters showing the known clades B.4, B.6, and B.12. These data correlated with the results of real-time PCR assays targeting canSNPs loci. Additionally, we detected two subtle sub-clusters. SplitsTree was used with standard-setting using the aligned SNPs from Parsnps. Together APC, HierBAPS, and SplitsTree enabled us to generate hypotheses about epidemiologic relationships between bacterial clusters and describing the distribution of isolates. Our data indicate that the choice of the typing technique can increase our understanding of the pathogenesis and transmission of diseases with the eventual for prevention. This is opening perspectives to be applied to other bacterial species. The data provide evidence that Germany might be the collision zone where the clade B.12, also known as the East European clade, overlaps with the clade B.6, also known as the Iberian clade. Described methods allow generating a new, more detailed perspective for F. tularensis subsp. holarctica phylogeny. These results may encourage to determine phylogenetic relationships and clustering of other bacteria the same way.

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.

Reclassification of Francisella noatunensis subsp. orientalis Ottem et al. 2009 as Francisella orientalis sp. nov., Francisella noatunensis subsp. chilensis subsp. nov. and emended description of Francisella noatunensis

Francisella noatunensis is a fastidious facultative intracellular bacterial pathogen that causes 'piscine francisellosis', a serious disease affecting both marine and fresh water farmed and wild fish worldwide. Currently two F. noatunensis subspecies are recognized, i.e. F. noatunensis subsp. noatunensis and F. noatunensis subsp. orientalis. In the present study, the taxonomy of F. noatunensis was revisited using a polyphasic approach, including whole genome derived parameters such as digital DNA-DNA hybridization, whole genome average nucleotide identity (wg-ANIm), whole genome phylogenetic analysis, whole genome G+C content, metabolic fingerprinting and chemotaxonomic analyses. The results indicated that isolates belonging to F. noatunensis subsp. orientalis represent a phenotypically and genetically homogenous taxon, clearly distinguishable from F. noatunensis subsp. noatunensis that fulfils requirements for separate species status. We propose, therefore, elevation of F. noatunensis subsp. orientalis to the species rank as Francisella orientalis sp. nov. with the type strain remaining as Ehime-1^T (DSM 21254^T=LMG 24544^T). Furthermore, we identified sufficient phenotypic and genetic differences between F. noatunensis subsp. noatunensis recovered from diseased farmed Atlantic salmon in Chile and those isolated from wild and farmed Atlantic cod in Northern Europe to warrant proposal of the Chilean as a novel F. noatunensis subspecies, i.e. Francisella noatunensis subsp. chilensis subsp. nov. with strain PQ1106^T (CECT 9798^T=NCTC14375^T) as the type strain. Finally, we emend the description of F. noatunensis by including further metabolic information and the description of atypical strains.

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.

A predicted Francisella tularensis DXD-motif glycosyltransferase blocks immune activation

Pathogens enhance their survival during infections by manipulating host defenses. Francisella tularensis evades innate immune responses, which we have found to be dependent on an understudied gene ybeX (FTL_0883/FTT_0615c). To understand the function of YbeX, we sought protein interactors in F. tularensis subsp. holarctica live vaccine strain (LVS). An unstudied Francisella protein co-immunoprecipitated with recombinant YbeX, which is a predicted glycosyltransferase with a DXD-motif. There are up to four genomic copies of this gene with identical sequence in strains of F. tularensis pathogenic to humans, despite ongoing genome decay. Disruption mutations were generated by intron insertion into all three copies of this glycosyltransferase domain containing gene in LVS, gdcA1-3. The resulting strains stimulated more cytokines from macrophages in vitro than wild-type LVS and were attenuated in two in vivo infection models. GdcA was released from LVS during culture and was sufficient to block NF-κB activation when expressed in eukaryotic cells. When co-expressed in zebrafish, GdcA and YbeX were synergistically lethal to embryo development. Glycosyltransferases with DXD-motifs are found in a variety of pathogens including NleB, an Escherichia coli type-III secretion system effector that inhibits NF-κB by antagonizing death receptor signaling. To our knowledge, GdcA is the first DXD-motif glycosyltransferase that inhibits NF-κB in immune cells. Together, these findings suggest DXD-motif glycosyltransferases may be a conserved virulence mechanism used by pathogenic bacteria to remodel host defenses.

Impact of Homologous Recombination on the Evolution of Prokaryotic Core Genomes

Microbial populations exchange genetic material through a process called homologous recombination. Although this process has been studied in particular organisms, we lack an understanding of its differential impact over the genome and across microbes with different life-styles. We used a common analytical framework to assess this process in a representative set of microorganisms. Our results uncovered important trends. First, microbes with different lifestyles are differentially impacted, with endosymbionts and obligate pathogens being those less prone to undergo this process. Second, certain genetic elements such as restriction-modification systems seem to be associated with higher rates of recombination. Most importantly, recombined genomes show the footprints of natural selection in which recombined regions preferentially contain genes that can be related to specific ecological adaptations. Taken together, our results clarify the relative contributions of factors modulating homologous recombination and show evidence for a clear a role of this process in shaping microbial genomes and driving ecological adaptations.

Homologous recombination (HR) enables the exchange of genetic material between and within species. Recent studies suggest that this process plays a major role in the microevolution of microbial genomes, contributing to core genome homogenization and to the maintenance of cohesive population structures. However, we still have a very poor understanding of the possible adaptive roles of intraspecific HR and of the factors that determine its differential impact across clades and lifestyles. Here we used a unified methodological framework to assess HR in 338 complete genomes from 54 phylogenetically diverse and representative prokaryotic species, encompassing different lifestyles and a broad phylogenetic distribution. Our results indicate that lifestyle and presence of restriction-modification (RM) machineries are among the main factors shaping HR patterns, with symbionts and intracellular pathogens having the lowest HR levels. Similarly, the size of exchanged genomic fragments correlated with the presence of RM and competence machineries. Finally, genes exchanged by HR showed functional enrichments which could be related to adaptations to different environments and ecological strategies. Taken together, our results clarify the factors underlying HR impact and suggest important adaptive roles of genes exchanged through this mechanism. Our results also revealed that the extent of genetic exchange correlated with lifestyle and some genomic features. Moreover, the genes in exchanged regions were enriched for functions that reflected specific adaptations, supporting identification of HR as one of the main evolutionary mechanisms shaping prokaryotic core genomes.

ISEScan: automated identification of insertion sequence elements in prokaryotic genomes

Motivation

The insertion sequence (IS) elements are the smallest but most abundant autonomous transposable elements in prokaryotic genomes, which play a key role in prokaryotic genome organization and evolution. With the fast growing genomic data, it is becoming increasingly critical for biology researchers to be able to accurately and automatically annotate ISs in prokaryotic genome sequences. The available automatic IS annotation systems are either providing only incomplete IS annotation or relying on the availability of existing genome annotations. Here, we present a new IS elements annotation pipeline to address these issues.

Results

ISEScan is a highly sensitive software pipeline based on profile hidden Markov models constructed from manually curated IS elements. ISEScan performs better than existing IS annotation systems when tested on prokaryotic genomes with curated annotations of IS elements. Applying it to 2784 prokaryotic genomes, we report the global distribution of IS families across taxonomic clades in Archaea and Bacteria.

Availability and implementation

ISEScan is implemented in Python and released as an open source software at https://github.com/xiezhq/ISEScan.

Contact

hatang@indiana.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Challenges of Francisella classification exemplified by an atypical clinical isolate

The accumulation of sequenced Francisella strains has made it increasingly apparent that the 16S rRNA gene alone is not enough to stratify the Francisella genus into precise and clinically useful classifications. Continued whole-genome sequencing of isolates will provide a larger base of knowledge for targeted approaches with broad applicability. Additionally, examination of genomic information on a case-by-case basis will help resolve outstanding questions regarding strain stratification. We report the complete genome sequence of a clinical isolate, designated here as F. novicida-like strain TCH2015, acquired from the lymph node of a 6-year-old male. Two features were atypical for F. novicida: exhibition of functional oxidase activity and additional gene content, including proposed virulence determinants. These differences, which could potentially impact virulence and clinical diagnosis, emphasize the need for more comprehensive methods to profile Francisella isolates. This study highlights the value of whole-genome sequencing, which will lead to a more robust database of environmental and clinical genomes and inform strategies to improve detection and classification of Francisella strains.

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.

Signatures of Selection at Drug Resistance Loci in Mycobacterium tuberculosis

Tuberculosis (TB) is the leading cause of death by an infectious disease, and global TB control efforts are increasingly threatened by drug resistance in Mycobacterium tuberculosis. Unlike most bacteria, where lateral gene transfer is an important mechanism of resistance acquisition, resistant M. tuberculosis arises solely by de novo chromosomal mutation. Using whole-genome sequencing data from two natural populations of M. tuberculosis, we characterized the population genetics of known drug resistance loci using measures of diversity, population differentiation, and convergent evolution. We found resistant subpopulations to be less diverse than susceptible subpopulations, consistent with ongoing transmission of resistant M. tuberculosis. A subset of resistance genes ("sloppy targets") were characterized by high diversity and multiple rare variants; we posit that a large genetic target for resistance and relaxation of purifying selection contribute to high diversity at these loci. For "tight targets" of selection, the path to resistance appeared narrower, evidenced by single favored mutations that arose numerous times in the phylogeny and segregated at markedly different frequencies in resistant and susceptible subpopulations. These results suggest that diverse genetic architectures underlie drug resistance in M. tuberculosis and that combined approaches are needed to identify causal mutations. Extrapolating from patterns observed for well-characterized genes, we identified novel candidate variants involved in resistance. The approach outlined here can be extended to identify resistance variants for new drugs, to investigate the genetic architecture of resistance, and when phenotypic data are available, to find candidate genetic loci underlying other positively selected traits in clonal bacteria. IMPORTANCEMycobacterium tuberculosis, the causative agent of tuberculosis (TB), is a significant burden on global health. Antibiotic treatment imposes strong selective pressure on M. tuberculosis populations. Identifying the mutations that cause drug resistance in M. tuberculosis is important for guiding TB treatment and halting the spread of drug resistance. Whole-genome sequencing (WGS) of M. tuberculosis isolates can be used to identify novel mutations mediating drug resistance and to predict resistance patterns faster than traditional methods of drug susceptibility testing. We have used WGS from natural populations of drug-resistant M. tuberculosis to characterize effects of selection for advantageous mutations on patterns of diversity at genes involved in drug resistance. The methods developed here can be used to identify novel advantageous mutations, including new resistance loci, in M. tuberculosis and other clonal pathogens.

Multilocus sequence analysis of homologous recombination and diversity in Arthrobacter sensu lato named species and glacier-inhabiting strains

Members of the bacterial genus Arthrobacter sensu lato are Gram-positive actinomycetes distributed worldwide and found in numerous environments including soil, water, glacier ice, and sewage. Homologous recombination is an important driving force in bacterial evolution, but its impact on Arthrobacter sensu lato evolution is poorly understood. We evaluated homologous recombination among 41 Arthrobacter sensu lato named species, using multilocus sequence analysis (MLSA). A high level of recombination was found, associated with strong diversification and a reticulate evolutionary pattern of Arthrobacter sensu lato. We also collected a total of 31 cold-adapted Arthrobacter sensu lato strains from two cold glaciers located in northwest China and two temperate glaciers in southwest China, and evaluated their diversity and population structure by MLSA. The glacier strains displayed high diversity, but rates of recombination among the four glacier groups were quite low, indicating that barriers to homologous recombination formed in the past among the populations on different glaciers. Our findings indicate that historical glaciation events shaped the contemporary distributions, taxonomic relationships, and phylogeographic patterns of Arthrobacter sensu lato species on glaciers.

Adaptation in a Fibronectin Binding Autolysin of Staphylococcus saprophyticus

Human-pathogenic bacteria are found in a variety of niches, including free-living, zoonotic, and microbiome environments. Identifying bacterial adaptations that enable invasive disease is an important means of gaining insight into the molecular basis of pathogenesis and understanding pathogen emergence. Staphylococcus saprophyticus, a leading cause of urinary tract infections, can be found in the environment, food, animals, and the human microbiome. We identified a selective sweep in the gene encoding the Aas adhesin, a key virulence factor that binds host fibronectin. We hypothesize that the mutation under selection (aas_2206A>C) facilitates colonization of the urinary tract, an environment where bacteria are subject to strong shearing forces. The mutation appears to have enabled emergence and expansion of a human-pathogenic lineage of S. saprophyticus. These results demonstrate the power of evolutionary genomic approaches in discovering the genetic basis of virulence and emphasize the pleiotropy and adaptability of bacteria occupying diverse niches. IMPORTANCEStaphylococcus saprophyticus is an important cause of urinary tract infections (UTI) in women; such UTI are common, can be severe, and are associated with significant impacts to public health. In addition to being a cause of human UTI, S. saprophyticus can be found in the environment, in food, and associated with animals. After discovering that UTI strains of S. saprophyticus are for the most part closely related to each other, we sought to determine whether these strains are specially adapted to cause disease in humans. We found evidence suggesting that a mutation in the gene aas is advantageous in the context of human infection. We hypothesize that the mutation allows S. saprophyticus to survive better in the human urinary tract. These results show how bacteria found in the environment can evolve to cause disease.

Phenetic Comparison of Prokaryotic Genomes Using k-mers

Bacterial genomics studies are getting more extensive and complex, requiring new ways to envision analyses. Using the Ray Surveyor software, we demonstrate that comparison of genomes based on their k-mer content allows reconstruction of phenetic trees without the need of prior data curation, such as core genome alignment of a species. We validated the methodology using simulated genomes and previously published phylogenomic studies of Streptococcus pneumoniae and Pseudomonas aeruginosa. We also investigated the relationship of specific genetic determinants with bacterial population structures. By comparing clusters from the complete genomic content of a genome population with clusters from specific functional categories of genes, we can determine how the population structures are correlated. Indeed, the strain clustering based on a subset of k-mers allows determination of its similarity with the whole genome clusters. We also applied this methodology on 42 species of bacteria to determine the correlational significance of five important bacterial genomic characteristics. For example, intrinsic resistance is more important in P. aeruginosa than in S. pneumoniae, and the former has increased correlation of its population structure with antibiotic resistance genes. The global view of the pangenome of bacteria also demonstrated the taxa-dependent interaction of population structure with antibiotic resistance, bacteriophage, plasmid, and mobile element k-mer data sets.

Expanding Francisella models: Pairing up the soil amoeba Dictyostelium with aquatic Francisella

The bacterial genus Francisella comprises highly pathogenic species that infect mammals, arthropods, fish and protists. Understanding virulence and host defense mechanisms of Francisella infection relies on multiple animal and cellular model systems. In this review, we want to summarize the most commonly used Francisella host model platforms and highlight novel, alternative model systems using aquatic Francisella species. Established mouse and macrophage models contributed extensively to our understanding of Francisella infection. However, murine and human cells display significant differences in their response to Francisella infection. The zebrafish and the amoeba Dictyostelium are well-established model systems for host-pathogen interactions and open up opportunities to investigate bacterial virulence and host defense. Comparisons between model systems using human and fish pathogenic Francisella species revealed shared virulence strategies and pathology between them. Hence, zebrafish and Dictyostelium might complement current model systems to find new vaccine candidates and contribute to our understanding of Francisella infection.

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.

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.

A molecular portrait of maternal sepsis from Byzantine Troy

Pregnancy complications are poorly represented in the archeological record, despite their importance in contemporary and ancient societies. While excavating a Byzantine cemetery in Troy, we discovered calcified abscesses among a woman’s remains. Scanning electron microscopy of the tissue revealed ‘ghost cells’, resulting from dystrophic calcification, which preserved ancient maternal, fetal and bacterial DNA of a severe infection, likely chorioamnionitis. Gardnerella vaginalis and Staphylococcus saprophyticus dominated the abscesses. Phylogenomic analyses of ancient, historical, and contemporary data showed that G. vaginalis Troy fell within contemporary genetic diversity, whereas S. saprophyticus Troy belongs to a lineage that does not appear to be commonly associated with human disease today. We speculate that the ecology of S. saprophyticus infection may have differed in the ancient world as a result of close contacts between humans and domesticated animals. These results highlight the complex and dynamic interactions with our microbial milieu that underlie severe maternal infections.

DOI:http://dx.doi.org/10.7554/eLife.20983.001

Why and how have some bacteria evolved to cause illness in humans? One way to study bacterial evolution is to search for ancient samples of bacteria and use DNA sequencing technology to investigate how modern bacteria have changed from their ancestors. Understanding the evolution process may help researchers to understand how some bacteria become resistant to the antibiotics designed to kill them.

Complications that occur during pregnancy, including bacterial infections, have long been a major cause of death for women. Now, Devault, Mortimer et al. have been able to sequence the DNA of bacteria found in tissue collected from a woman buried 800 years ago in a cemetery in Troy. Some of the woman’s tissues had been well preserved because they had calcified (probably as the result of infection), which preserved their structure in a mineralized layer. Two mineralized “nodules” in the body appear to be the remains of abscesses. Some of the human DNA in the nodules came from a male, suggesting that the woman was pregnant with a boy and that the abscesses formed in placental tissue.

Sequencing the DNA of the bacteria in the abscess allowed Devault, Mortimer et al. to diagnose the woman’s infection, which was caused by two types of bacteria. One species, called Gardnerella vaginalis, is found in modern pregnancy-related infections. The DNA of the ancient samples was similar to that of modern bacteria. The other bacteria species was an ancient form of Staphylococcus saprophyticus, a type of bacteria that causes urinary tract infections. However, the DNA of the ancient S. saprophyticus bacteria is quite different to that of the bacteria found in modern humans. Instead, their DNA sequence appears more similar to forms of the bacteria that infect currently livestock. As humans lived closely with their livestock at the time the woman lived, her infection may be due to a type of bacteria that passed easily between humans and animals.

Overall, the results suggest that the disease-causing properties of bacteria can arise from a wide range of sources. In addition, Devault, Mortimer et al. have demonstrated that certain types of tissue found in archeological remains are a potential gold mine of information about the evolution of bacteria and other microbes found in the human body.

DOI:http://dx.doi.org/10.7554/eLife.20983.002

Genome-scale metabolic reconstruction for the insidious bacterium in aquaculture Piscirickettsia salmonis

Piscirickettsia salmonis is a fish bacterium that causes the disease piscirickettsiosis in salmonids. This pathology is partially controlled by vaccines. The lack of knowledge has hindered its culture on laboratory and industrial scale. The study describes the metabolic phenotype of P. salmonis in culture. This study presents the first genome-scale model (iPF215) of the LF-89 strain of P. salmonis, describing the central metabolic pathway, biosynthesis and molecule degradation and transport mechanisms. The model was adjusted with experiment data, allowing the identification of the capacities that were not predicted by the automatic annotation of the genome sequences. The iPF215 model is comprised of 417 metabolites, 445 reactions and 215 genes, was used to reproduce the growth of P. salmonis (μ_max 0.052±0.005h^-1). The metabolic reconstruction of the P. salmonis LF-89 strain obtained in this research provides a baseline that describes the metabolic capacities of the bacterium and is the basis for developing improvements to its cultivation for vaccine formulation.

Sensitive Detection of Francisella tularensis Directly from Whole Blood by Use of the GeneXpert System

Francisella tularensis is a potential bioterrorism agent that is highly infectious at very low doses. Diagnosis of tularemia by blood culture and nucleic acid-based diagnostic tests is insufficiently sensitive. Here, we demonstrate a highly sensitive F. tularensis assay that incorporates sample processing and detection into a single cartridge suitable for point-of-care detection. The assay limit of detection (LOD) and dynamic range were determined in a filter-based cartridge run on the GeneXpert system. F. tularensis DNA in buffer or CFU of F. tularensis was spiked into human or macaque blood. To simulate detection in human disease, the assay was tested on blood drawn from macaques infected with F. tularensis Schu S4 at daily intervals. Assay detection was compared to that with a conventional quantitative PCR (qPCR) assay and blood culture. The assay LOD was 0.1 genome equivalents (GE) per reaction and 10 CFU/ml F. tularensis in both human and macaque blood. In infected macaques, the assay detected F. tularensis on days 1 to 4 postinfection in 21%, 17%, 60%, and 83% of macaques, respectively, compared to conventional qPCR positivity rates of 0%, 0%, 30%, and 100% and CFU detection of blood culture at 0%, 0%, 0%, and 10% positive, respectively. Assay specificity was 100%. The new cartridge-based assay can rapidly detect F. tularensis in bloodstream infections directly in whole blood at the early stages of infection with a sensitivity that is superior to that of other methods. The simplicity of the automated testing procedures may make this test suitable for rapid point-of-care detection.

Long-range dispersal moved Francisella tularensis into Western Europe from the East

For many infections transmitting to humans from reservoirs in nature, disease dispersal patterns over space and time are largely unknown. Here, a reversed genomics approach helped us understand disease dispersal and yielded insight into evolution and biological properties of Francisella tularensis, the bacterium causing tularemia. We whole-genome sequenced 67 strains and characterized by single-nucleotide polymorphism assays 138 strains, collected from individuals infected 1947-2012 across Western Europe. We used the data for phylogenetic, population genetic and geographical network analyses. All strains (n=205) belonged to a monophyletic population of recent ancestry not found outside Western Europe. Most strains (n=195) throughout the study area were assigned to a star-like phylogenetic pattern indicating that colonization of Western Europe occurred via clonal expansion. In the East of the study area, strains were more diverse, consistent with a founder population spreading from east to west. The relationship of genetic and geographic distance within the F. tularensis population was complex and indicated multiple long-distance dispersal events. Mutation rate estimates based on year of isolation indicated null rates; in outbreak hotspots only, there was a rate of 0.4 mutations/genome/year. Patterns of nucleotide substitution showed marked AT mutational bias suggestive of genetic drift. These results demonstrate that tularemia has moved from east to west in Europe and that F. tularensis has a biology characterized by long-range geographical dispersal events and mostly slow, but variable, replication rates. The results indicate that mutation-driven evolution, a resting survival phase, genetic drift and long-distance geographical dispersal events have interacted to generate genetic diversity within this species.

Murine survival of infection with Francisella novicida and protection against secondary challenge is critically dependent on B lymphocytes

Respiratory infection of mice with Francisella novicida has recently been used as a model for the highly virulent human pathogen Francisella tularensis. Similar to F. tularensis, even small doses of F. novicida administered by respiratory routes are lethal for inbred laboratory mice. This feature obviously limits study of infection-induced immunity. Parenteral sublethal infections of mice with F. novicida are feasible, but the resulting immune responses are incompletely characterized. Here we use parenteral intradermal (i.d.) and intraperitoneal (i.p.) F. novicida infections of C57BL/6J mice to determine the role of B cells in controlling primary and secondary F. novicida infections. Despite developing comparable levels of F. novicida-primed T cells, B cell knockout mice were much more susceptible to both primary i.d. infection and secondary i.p. challenge than wild type (normal) C57BL/6J mice. Transfer of F. novicida-immune sera to either wild type C57BL/6J mice or to B cell knockout mice did not appreciably impact survival of subsequent lethal F. novicida challenge. However, F. novicida-immune mice that were depleted of T cells after priming but just before challenge survived and cleared secondary i.p. F. novicida challenge. Collectively these results indicate that B cells, if not serum antibodies, play a major role in controlling F. novicida infections in mice.

Inflammasome-Independent NLRP3 Restriction of a Protective Early Neutrophil Response to Pulmonary Tularemia

Francisella tularensis (Ft) causes a frequently fatal, acute necrotic pneumonia in humans and animals. Following lethal Ft infection in mice, infiltration of the lungs by predominantly immature myeloid cells and subsequent myeloid cell death drive pathogenesis and host mortality. However, following sub-lethal Ft challenge, more mature myeloid cells are elicited and are protective. In addition, inflammasome-dependent IL-1β and IL-18 are important for protection. As Nlrp3 appears dispensable for resistance to infection with Francisella novicida, we considered its role during infection with the virulent Type A strain SchuS4 and the attenuated Type B live vaccine strain LVS. Here we show that both in vitro macrophage and in vivo IL-1β and IL-18 responses to Ft LVS and SchuS4 involve both the Aim2 and Nlrp3 inflammasomes. However, following lethal infection with Francisella, IL-1r-, Caspase-1/11-, Asc- and Aim2-deficient mice exhibited increased susceptibility as expected, while Nlrp3-deficient mice were more resistant. Despite reduced levels of IL-1β and IL-18, in the absence of Nlrp3, Ft infected mice have dramatically reduced lung pathology, diminished recruitment and death of immature myeloid cells, and reduced bacterial burden in comparison to wildtype and inflammasome-deficient mice. Further, increased numbers of mature neutrophil appear in the lung early during lethal Ft infection in Nlrp3-deficient mice. Finally, Ft infection induces myeloid and lung stromal cell death that in part requires Nlrp3, is necrotic/necroptotic in nature, and drives host mortality. Thus, Nlrp3 mediates an inflammasome-independent process that restricts the appearance of protective mature neutrophils and promotes lethal necrotic lung pathology.

Introduction and persistence of tularemia in Bulgaria

Introduction

Outbreaks of the zoonotic disease tularemia occurred in north-east Bulgaria in the 1960s. Then came 30 years of epidemiological silence until new outbreaks occurred in west Bulgaria in the 1990s. To investigate how bacterial strains of Francisella tularensis causing tularemia in wildlife and humans in the 1960s and the 1990s were related, we explored their genetic diversity.

Material and methods

Ten F. tularensis genomes from the 1960s (n=3) and the 1990s (n=7) were sequenced, assigned to canonical single-nucleotide polymorphism (canSNP) clades, and compared to reference genomes. We developed four new canSNP polymerase chain reaction (PCR) assays based on the genome sequence information.

Results and discussion

The genetic analysis showed that the outbreaks in the 1960s as well as in the 1990s involved multiple clones and new genetic diversity. The smallest genetic difference found between any of the Bulgarian strains was five SNPs between the strains L2 and 81 isolated 43 years apart, indicating that F. tularensis may persist locally over long time periods without causing outbreaks. The existence of genetically highly similar strain-pairs isolated the same year in the same area from different hosts supports a hypothesis of local expansion of clones during outbreaks. Close relationship (two SNPs) was found between one strain isolated 1961 in northeast Bulgaria and one strain isolated 5 years before in USSR. Historical data coinciding with the actual time point describe the introduction of water rats from USSR into the Bulgarian outbreak area, which may explain the close genetic relationship and the origin of the outbreak.

Conclusion

Genome analysis of strains from two outbreaks in the 1960s and the 1990s provided valuable information on the genetic diversity and persistence of F. tularensis in Bulgaria.

Francisella tularensis IglG Belongs to a Novel Family of PAAR-Like T6SS Proteins and Harbors a Unique N-terminal Extension Required for Virulence

The virulence of Francisella tularensis, the etiological agent of tularemia, relies on an atypical type VI secretion system (T6SS) encoded by a genomic island termed the Francisella Pathogenicity Island (FPI). While the importance of the FPI in F. tularensis virulence is clearly established, the precise role of most of the FPI-encoded proteins remains to be deciphered. In this study, using highly virulent F. tularensis strains and the closely related species F. novicida, IglG was characterized as a protein featuring a unique α-helical N-terminal extension and a domain of unknown function (DUF4280), present in more than 250 bacterial species. Three dimensional modeling of IglG and of the DUF4280 consensus protein sequence indicates that these proteins adopt a PAAR-like fold, suggesting they could cap the T6SS in a similar way as the recently described PAAR proteins. The newly identified PAAR-like motif is characterized by four conserved cysteine residues, also present in IglG, which may bind a metal atom. We demonstrate that IglG binds metal ions and that each individual cysteine is required for T6SS-dependent secretion of IglG and of the Hcp homologue, IglC and for the F. novicida intracellular life cycle. In contrast, the Francisella-specific N-terminal α-helical extension is not required for IglG secretion, but is critical for F. novicida virulence and for the interaction of IglG with another FPI-encoded protein, IglF. Altogether, our data suggest that IglG is a PAAR-like protein acting as a bi-modal protein that may connect the tip of the Francisella T6SS with a putative T6SS effector, IglF.

Francisella tularensis is a highly pathogenic bacterium causing tularemia. Its ability to cause disease is linked to its ability to replicate in the macrophage cytosol. The intracellular life cycle of Francisella is controlled by a type VI secretion system (T6SS), which is thought to inject effectors into the host cell to allow bacterial escape into the host cytosol. The molecular mechanisms behind this process are still largely unclear. In this work, we identify IglG as a protein with two important domains, one conserved in proteins from more than 250 bacterial species (DUF4280, renamed here as PAAR-like domain) and one specific for the Francisella genus. Using protein sequence analysis and three-dimensional structure predictions, comparative modeling and biochemistry approaches, our data demonstrate that IglG is a metal-binding protein that based on its PAAR-like domain might cap the VgrG spike of the T6SS and act as a membrane-puncturing protein. Furthermore, we identified that the Francisella-specific domain is directly involved in forming a protein complex with another virulence protein, IglF. This work, in addition to enhancing the molecular understanding of the Francisella T6SS, defines the features of the conserved DUF4280, a novel PAAR-like domain involved in type VI secretion (T6S) of many bacterial species.

Clonality of erythromycin resistance in Francisella tularensis

OBJECTIVES

We analysed diverse strains of Francisella tularensis subsp. holarctica to assess if its division into biovars I and II is associated with specific mutations previously linked to erythromycin resistance and to determine the distribution of this resistance trait across this subspecies.

METHODS

Three-hundred and fourteen F. tularensis subsp. holarctica strains were tested for erythromycin susceptibility and whole-genome sequences for these strains were examined for SNPs in genes previously associated with erythromycin resistance. Each strain was assigned to a global phylogenetic framework using genome-wide canonical SNPs. The contribution of a specific SNP to erythromycin resistance was examined using allelic exchange. The geographical distribution of erythromycin-resistant F. tularensis strains was further investigated by literature search.

RESULTS

There was a perfect correlation between biovar II strains (erythromycin resistance) and the phylogenetic group B.12. Only B.12 strains had an A → C SNP at position 2059 in the three copies of the rrl gene. Introducing 2059C into an rrl gene of an erythromycin-susceptible F. tularensis strain resulted in resistance. An additional 1144 erythromycin-resistant strains were identified from the scientific literature, all of them from Eurasia.

CONCLUSIONS

Erythromycin resistance in F. tularensis is caused by an A2059C rrl gene mutation, which exhibits a strictly clonal inheritance pattern found only in phylogenetic group B.12. This group is an extremely successful clone, representing the most common type of F. tularensis throughout Eurasia.

Tularemia vaccines

Francisella tularensis is the causative agent of the potentially lethal disease tularemia. Due to a low infectious dose and ease of airborne transmission, Francisella is classified as a category A biological agent. Despite the possible risk to public health, there is no safe and fully licensed vaccine. A potential vaccine candidate, an attenuated live vaccine strain, does not fulfil the criteria for general use. In this review, we will summarize existing and new candidates for live attenuated and subunit vaccines.

Towards Development of Improved Serodiagnostics for Tularemia by Use of Francisella tularensis Proteome Microarrays

Tularemia in humans is caused mainly by two subspecies of the Gram-negative facultative anaerobe Francisella tularensis: F. tularensis subsp. tularensis (type A) and F. tularensis subsp. holarctica (type B). The current serological test for tularemia is based on agglutination of whole organisms, and the reactive antigens are not well understood. Previously, we profiled the antibody responses in type A and B tularemia cases in the United States using a proteome microarray of 1,741 different proteins derived from the type A strain Schu S4. Fifteen dominant antigens able to detect antibodies to both types of infection were identified, although these were not validated in a different immunoassay format. Since type A and B subspecies are closely related, we hypothesized that Schu S4 antigens would also have utility for diagnosing type B tularemia caused by strains from other geographic locations. To test this, we probed the Schu S4 array with sera from 241 type B tularemia cases in Spain. Despite there being no type A strains in Spain, we confirmed the responses against some of the same potential serodiagnostic antigens reported previously, as well as determined the responses against additional potential serodiagnostic antigens. Five potential serodiagnostic antigens were evaluated on immunostrips, and two of these (FTT1696/GroEL and FTT0975/conserved hypothetical protein) discriminated between the Spanish tularemia cases and healthy controls. We conclude that antigens from the type A strain Schu S4 are suitable for detection of antibodies from patients with type B F. tularensis infections and that these can be used for the diagnosis of tularemia in a deployable format, such as the immunostrip.

Francisella Inflammasomes: Integrated Responses to a Cytosolic Stealth Bacterium

Francisella tularensis is a facultative intracellular bacterium causing tularemia, a zoonotic disease. Francisella replicates in the macrophage cytosol and eventually triggers cytosolic immune responses. In murine macrophages, Francisella novicida and Francisella tularensis live vaccine strain lyse in the host cytosol and activate the cytosolic DNA receptor Aim2. Here, we review the mechanisms leading or contributing to Aim2 inflammasome activation, including the role of TLRs and of IFN signaling and the implication of the guanylate-binding proteins 2 and 5 in triggering cytosolic bacteriolysis. Furthermore, we present how this cytosolic Gram-negative bacterium escapes recognition by caspase-11 but can trigger a non-canonical caspase-8 inflammasome. In addition, we highlight the differences in inflammasome activation in murine and human cells with pyrin, NLRP3, and AIM2 involved in sensing Francisella in human phagocytes. From a bacterial prospective, we describe the hiding strategy of Francisella to escape recognition by innate sensors and to resist to bacteriolysis in the host cytosol. Finally, we discuss the inability of the inflammasome sensors to detect F. tularensis subspecies tularensis strains, making them highly pathogenic stealth microbes.

Scaffolding of a bacterial genome using MinION nanopore sequencing

Second generation sequencing has revolutionized genomic studies. However, most genomes contain repeated DNA elements that are longer than the read lengths achievable with typical sequencers, so the genomic order of several generated contigs cannot be easily resolved. A new generation of sequencers offering substantially longer reads is emerging, notably the Pacific Biosciences (PacBio) RS II system and the MinION system, released in early 2014 by Oxford Nanopore Technologies through an early access program. The latter has highly advantageous portability and sequences samples by measuring changes in ionic current when single-stranded DNA molecules are translocated through nanopores. We show that the MinION system produces long reads with high mapability that can be used for scaffolding bacterial genomes, despite currently producing substantially higher error rates than PacBio reads. With further development we anticipate that MinION will be useful not only for assembling genomes, but also for rapid detection of organisms, potentially in the field.

Genomic analyses of Francisella tularensis strains confirm disease transmission from drinking water sources, Turkey, 2008, 2009 and 2012

Waterborne epidemics of tularaemia caused by Francisella tularensis are increasingly reported in Turkey. We have used whole genome sequencing to investigate if F. tularensis isolated from patients could be traced back to drinking water sources. Tonsil swabs from 33 patients diagnosed with oropharyngeal tularaemia in three outbreaks and 140 water specimens were analysed. F. tularensis subsp. holarctica was confirmed by microagglutination and PCR in 12 patients and five water specimens. Genomic analysis of three pairs of patient and water isolates from outbreaks in Sivas, Çorum, and Kocaeli showed the isolates to belong to two new clusters of the F. tularensis B.12 genetic clade. The clusters were defined by 19 and 15 single nucleotide polymorphisms (SNPs) in a multiple alignment based on 507 F. tularensis genomes. One synonymous SNP was chosen as a new canonical SNP (canSNP) for each cluster for future use in diagnostic assays. No SNP was identified between the genomes from the patient–water pair of isolates from Kocaeli, one SNP between the pair of isolates from Sivas, whereas the pair from Çorum differed at seven SNPs. These results illustrate the power of whole genome sequencing for tracing F. tularensis patient isolates back to their environmental source.

Francisella tularensis Subtype A.II Genomic Plasticity in Comparison with Subtype A.I

Although Francisella tularensis is considered a monomorphic intracellular pathogen, molecular genotyping and virulence studies have demonstrated important differences within the tularensis subspecies (type A). To evaluate genetic variation within type A strains, sequencing and assembly of a new subtype A.II genome was achieved for comparison to other completed F. tularensis type A genomes. In contrast with the F. tularensis A.I strains (SCHU S4, FSC198, NE061598, and TI0902), substantial genomic variation was observed between the newly sequenced F. tularensis A.II strain (WY-00W4114) and the only other publically available A.II strain (WY96-3418). Genome differences between WY-00W4114 and WY96-3418 included three major chromosomal translocations, 1580 indels, and 286 nucleotide substitutions of which 159 were observed in predicted open reading frames and 127 were located in intergenic regions. The majority of WY-00W4114 nucleotide deletions occurred in intergenic regions, whereas most of the insertions and substitutions occurred in predicted genes. Of the nucleotide substitutions, 48 (30%) were synonymous and 111 (70%) were nonsynonymous. WY-00W4114 and WY96-3418 nucleotide polymorphisms were predominantly G/C to A/T allelic mutations, with WY-00W4114 having more A+T enrichment. In addition, the A.II genomes contained a considerably higher number of intact genes and longer repetitive sequences, including transposon remnants than the A.I genomes. Together these findings support the premise that F. tularensis A.II may have a fitness advantage compared to the A.I subtype due to the higher abundance of functional genes and repeated chromosomal sequences. A better understanding of the selective forces driving F. tularensis genetic diversity and plasticity is needed.

Rapid high resolution genotyping of Francisella tularensis by whole genome sequence comparison of annotated genes ("MLST+")

The zoonotic disease tularemia is caused by the bacterium Francisella tularensis. This pathogen is considered as a category A select agent with potential to be misused in bioterrorism. Molecular typing based on DNA-sequence like canSNP-typing or MLVA has become the accepted standard for this organism. Due to the organism's highly clonal nature, the current typing methods have reached their limit of discrimination for classifying closely related subpopulations within the subspecies F. tularensis ssp. holarctica. We introduce a new gene-by-gene approach, MLST+, based on whole genome data of 15 sequenced F. tularensis ssp. holarctica strains and apply this approach to investigate an epidemic of lethal tularemia among non-human primates in two animal facilities in Germany. Due to the high resolution of MLST+ we are able to demonstrate that three independent clones of this highly infectious pathogen were responsible for these spatially and temporally restricted outbreaks.

Importance of host cell arginine uptake in Francisella phagosomal escape and ribosomal protein amounts

Upon entry into mammalian host cells, the pathogenic bacterium Francisella must import host cell arginine to multiply actively in the host cytoplasm. We identified and functionally characterized an arginine transporter (hereafter designated ArgP) whose inactivation considerably delayed bacterial phagosomal escape and intracellular multiplication. Intramacrophagic growth of the ΔargP mutant was fully restored upon supplementation of the growth medium with excess arginine, in both F. tularensis subsp. novicida and F. tularensis subsp. holarctica LVS, demonstrating the importance of arginine acquisition in these two subspecies. High-resolution mass spectrometry revealed that arginine limitation reduced the amount of most of the ribosomal proteins in the ΔargP mutant. In response to stresses such as nutritional limitation, repression of ribosomal protein synthesis has been observed in all kingdoms of life. Arginine availability may thus contribute to the sensing of the intracellular stage of the pathogen and to trigger phagosomal egress. All MS data have been deposited in the ProteomeXchange database with identifier PXD001584 (http://proteomecentral.proteomexchange.org/dataset/PXD001584).