Genomic deletion marking an emerging subclone of Francisella tularensis subsp. holarctica in France and the Iberian Peninsula

Francisella tularensis subsp. holarctica is widely disseminated in North America and the boreal and temperate regions of the Eurasian continent. Comparative genomic analyses identified a 1.59-kb genomic deletion specific to F. tularensis subsp. holarctica isolates from Spain and France. Phylogenetic analysis of strains carrying this deletion by multiple-locus variable-number tandem repeat analysis showed that the strains comprise a highly related set of genotypes, implying that these strains were recently introduced or recently emerged by clonal expansion in France and the Iberian Peninsula.     

Phylogeography of Francisella tularensis ssp. holarctica in France

Aim: To investigate the phylogeography of French Francisella tularensis ssp. holarctica isolates. Methods and Results: Canonical SNPs and MLVA were used to genotype 103 French F. tularensis ssp. holarctica isolates. We confirmed the presence of one subclade, the central and western European group (B.Br.FTNF002‐00), and identified four major MLVA genotypes with no obvious geographical differentiation. Conclusions: The lack of geographical resolution among MLVA genotypes suggests rapid dispersal, convergent evolution or a combination of the two. Significance and Impact of the Study: This study expands knowledge of the phylogeography of one of the two dominant European F. tularensis ssp. holarctica subclades and illustrates the need for additional SNP discovery within this subclade.     

Paired-end sequence mapping detects extensive genomic rearrangement and translocation during divergence of Francisella tularensis subsp. tularensis and Francisella tularensis subsp. holarctica populations

Comparative genome hybridization of the Francisella tularensis subsp. tularensis and F. tularensis subsp. holarctica populations have shown that genome content is highly conserved, with relatively few genes in the F. tularensis subsp. tularensis genome being absent in other F. tularensis subspecies. To determine if organization of the genome differs between global populations of F. tularensis subsp. tularensis and F. tularensis subsp. holarctica, we have used paired-end sequence mapping (PESM) to identify regions of the genome where synteny is broken. The PESM approach compares the physical distances between paired-end sequencing reads of a library of a wild-type reference F. tularensis subsp. holarctica strain to the predicted lengths between the reads based on map coordinates of two different F. tularensis genome sequences. A total of 17 different continuous regions were identified in the F. tularensis subsp. holarctica genome (CR(holar)(c)(tica)) which are noncontiguous in the F. tularensis subsp. tularensis genome. Six of the 17 different CR(holarctica) are positioned as adjacent pairs in the F. tularensis subsp. tularensis genome sequence but are translocated in F. tularensis subsp. holarctica, implying that their arrangements are ancestral in F. tularensis subsp. tularensis and derived in F. tularensis subsp. holarctica. PCR analysis of the CR(holarctica) in 88 additional F. tularensis subsp. tularensis and F. tularensis subsp. holarctica isolates showed that the arrangements of the CR(holarctica) are highly conserved, particularly in F. tularensis subsp. holarctica, consistent with the hypothesis that global populations of F. tularensis subsp. holarctica have recently experienced a periodic selection event or they have emerged from a recent clonal expansion. Two unique F. tularensis subsp. tularensis-like strains were also observed which likely are derived from evolutionary intermediates and may represent a new taxonomic unit.     

Re-assembly of the genome of Francisella tularensis subsp. holarctica OSU18

Francisella tularensis is a highly infectious human intracellular pathogen that is the causative agent of tularemia. It occurs in several major subtypes, including the live vaccine strain holarctica (type B). F. tularensis is classified as category A biodefense agent in part because a relatively small number of organisms can cause severe illness. Three complete genomes of subspecies holarctica have been sequenced and deposited in public archives, of which OSU18 was the first and the only strain for which a scientific publication has appeared. We re-assembled the OSU18 strain using both de novo and comparative assembly techniques, and found that the published sequence has two large inversion mis-assemblies. We generated a corrected assembly of the entire genome along with detailed information on the placement of individual reads within the assembly. This assembly will provide a more accurate basis for future comparative studies of this pathogen.     

Borrelia garinii and Francisella tularensis subsp. holarctica detected in migratory shorebirds in Portugal

Migratory shorebirds use, among many, the East Atlantic Flyway that links breeding areas as north as Tundra habitats to aquatic wintering grounds in West Africa. As a consequence, they are potentially important in the spread of global zoonotic diseases transmitted by ticks, such as Lyme borreliosis and tularemia—two diseases previously detected in Portugal. In this study, we looked at the infection status of seven populations of shorebirds during their migration, breeding, or wintering in the Portuguese wetlands to access if they carry these pathogens and to discuss their potential risk in the Portuguese wetlands. A total of 212 migratory shorebirds captured in the Tagus and Sado estuaries; key staging and wintering sites in this flyway and important breeding areas for some species were analyzed for the presence of Borrelia burgdorferi sensu lato and Francisella tularensis. In the present study, B. garinii was identified in seven (3%) specimens (five black-tailed godwits Limosa limosa, one common redshank Tringa totanus, and one little stint Calidris minuta), whereas F. tularensis subsp. holarctica was identified in one (0.4%) little stint. To our knowledge, this is the first evidence that shorebirds that migrate through or winter in Portugal transport these pathogens, potentially contributing for their introduction along the flyway, including the Mediterranean region.     

Comparative proteome analysis of fractions enriched for membrane-associated proteins from Francisella tularensis subsp. tularensis and F. tularensis subsp. holarctica strains

The facultative intracellular pathogen Francisella tularensis is the causative agent of the serious infectious disease tularemia. Despite intensive research, the virulence factors and pathogenetic mechanisms remain largely unknown. To identify novel putative virulence factors, we carried out a comparative proteome analysis of fractions enriched for membrane-associated proteins isolated from the highly virulent subspecies tularensis strain SCHU S4 and three representatives of subspecies holarctica of different virulence including the live vaccine strain. We identified six proteins uniquely expressed and four proteins expressed at significantly higher levels by SCHU S4 compared to the ssp. holarctica strains. Four other protein spots represented mass and charge variants and seven spots were charge variants of proteins occurring in the ssp. holarctica strains. The genes encoding proteins of particular interest were examined by sequencing in order to confirm and explain the findings of the proteome analysis. Our studies suggest that the subspecies tularensis-specific proteins represent novel potential virulence factors.     

Glycosylation of DsbA in Francisella tularensis subsp. tularensis

In Francisella tularensis subsp. tularensis, DsbA has been shown to be an essential virulence factor and has been observed to migrate to multiple protein spots on two-dimensional electrophoresis gels. In this work, we show that the protein is modified with a 1,156-Da glycan moiety in O-linkage. The results of mass spectrometry studies suggest that the glycan is a hexasaccharide, comprised of N-acetylhexosamines, hexoses, and an unknown monosaccharide. Disruption of two genes within the FTT0789-FTT0800 putative polysaccharide locus, including a galE homologue (FTT0791) and a putative glycosyltransferase (FTT0798), resulted in loss of glycan modification of DsbA. The F. tularensis subsp. tularensis ΔFTT0798 and ΔFTT0791::Cm mutants remained virulent in the murine model of subcutaneous tularemia. This indicates that glycosylation of DsbA does not play a major role in virulence under these conditions. This is the first report of the detailed characterization of the DsbA glycan and putative role of the FTT0789-FTT0800 gene cluster in glycan biosynthesis.     

Deletion of IglH in virulent Francisella tularensis subsp. holarctica FSC200 strain results in attenuation and provides protection against the challenge with the parental strain

Francisella tularensis, the causative agent of tularemia, is a highly infectious intracellular pathogen with no licensed vaccine available today. The recent search for genome sequences involved in F. tularensis virulence mechanisms led to the identification of the 30-kb region defined as a Francisella pathogenicity island (FPI). In our previous iTRAQ study we described the concerted upregulation of some FPI proteins in different F. tularensis strains cultivated under stress conditions. Among them we identified the IglH protein whose role in Francisella virulence has not been characterized yet. In this work we deleted the iglH gene in a European clinical isolate of F. tularensis subsp. holarctica FSC200. We showed that the iglH gene is necessary for intracellular growth and escape of F. tularensis from phagosomes. We also showed that the iglH mutant is avirulent in a mouse model of infection and persists in the organs for about three weeks after infection. Importantly, mice vaccinated by infection with the iglH mutant were protected against subcutaneous challenge with the fully virulent parental FSC200 strain. This is the first report of a defined subsp. holarctica FPI deletion strain that provides protective immunity against subsequent subcutaneous challenge with a virulent isolate of F. tularensis subsp. holarctica.     

Genomic Markers for Differentiation of Francisella tularensis subsp. tularensis A.I and A.II Strains

Tularemia is caused by two subspecies of Francisella tularensis, F. tularensis subsp. tularensis (type A) and F. tularensis subsp. holarctica (type B). F. tularensis subsp. tularensis is further subdivided into two genetically distinct populations (A.I and A.II) that differ with respect to geographical location, anatomical source of recovered isolates, and disease outcome. Using two human clinical isolates, suppression subtractive hybridization was performed to identify 13 genomic regions of difference between A.I and A.II strains. Two PCR assays, one to identify A.I and A.II as well as to discriminate between F. tularensis subsp. holarctica and F. novicida and another specific for A.I, were developed. This is the first report to identify and characterize conserved genomic differences between A.I and A.II.     

Characterization of protein glycosylation in Francisella tularensis subsp. holarctica: identification of a novel glycosylated lipoprotein required for virulence

FTH_0069 is a previously uncharacterized strongly immunoreactive protein that has been proposed to be a novel virulence factor in Francisella tularensis. Here, the glycan structure modifying two C-terminal peptides of FTH_0069 was identified utilizing high resolution, high mass accuracy mass spectrometry, combined with in-source CID tandem MS experiments. The glycan observed at m/z 1156 was determined to be a hexasaccharide, consisting of two hexoses, three N-acetylhexosamines, and an unknown monosaccharide containing a phosphate group. The monosaccharide sequence of the glycan is tentatively proposed as X-P-HexNAc-HexNAc-Hex-Hex-HexNAc, where X denotes the unknown monosaccharide. The glycan is identical to that of DsbA glycoprotein, as well as to one of the multiple glycan structures modifying the type IV pilin PilA, suggesting a common biosynthetic pathway for the protein modification. Here, we demonstrate that the glycosylation of FTH_0069, DsbA, and PilA was affected in an isogenic mutant with a disrupted wbtDEF gene cluster encoding O-antigen synthesis and in a mutant with a deleted pglA gene encoding pilin oligosaccharyltransferase PglA. Based on our findings, we propose that PglA is involved in both pilin and general F. tularensis protein glycosylation, and we further suggest an inter-relationship between the O-antigen and the glycan synthesis in the early steps in their biosynthetic pathways.     

Genetic diversity of Francisella tularensis subspecies holarctica strains isolated in Japan

The recently developed MLVA has high discriminatory power for the typing of individual strains or isolates of Francisella tularensis. In the present study, MLVA was applied to 33 Japanese F. tularensis subspecies holarctica strains to examine the genetic diversity of F. tularensis isolated. Among the seven VNTR loci analyzed, Ft-M2, Ft-M10, and Ft-M20 loci showed high genetic polymorphism in Japanese strains, whereas Ft-M3 was most variable in non-Japanese strains. These results provide novel extended information about the genomic diversity among the strains of F. tularensis ssp. holarctica distributed in Japan and enable determination of whether a given isolate is indigenous to Japan by examining these loci using MLVA.     

iTRAQ quantitative analysis of Francisella tularensis ssp. holarctica live vaccine strain and Francisella tularensis ssp. tularensis SCHU S4 response to different temperatures and stationary phases of growth

Proteomics has been shown to significantly contribute to the investigation of the pathogenicity of the extremely infectious bacteria Francisella tularensis. In this study, the authors employed iTRAQ quantitative proteomic analysis in order to monitor alterations in proteomes of F. tularensis ssp. holarctica live vaccine strain and F. tularensis ssp. tularensis SCHU S4 associated with the cultivation at different temperatures or in the stationary phase. Correlated production of the identified proteins studied by the exploratory statistical analysis revealed novel candidates for virulence factors that were regulated in a similar manner to the genes encoded in the Francisella Pathogenicity Island. Moreover, the assessment of the adaptation of live vaccine strain and SCHU S4 strain to the examined stimuli uncovered differences in their physiological responses to the stationary phase of growth.     

Construction of targeted insertion mutations in Francisella tularensis subsp. novicida

Francisella tularensis is one of the most deadly bacterial agents, yet most of the genetic determinants of pathogenesis are still unknown. We have developed an efficient targeted mutagenesis strategy in the model organism F. tularensis subsp. novicida by utilizing universal priming of optimized antibiotic resistance cassettes and splicing by overlap extension (SOE). This process enables fast and efficient construction of targeted insertion mutations in F. tularensis subsp. novicida that have characteristics of nonpolar mutations; optimized targeted mutagenesis strategies will promote the study of this mysterious bacterium and facilitate vaccine development against tularemia. Moreover the general strategy of gene disruption by PCR-based antibiotic resistance cassette insertion is broadly applicable to many bacterial species.     

Biology of Francisella tularensis subspecies holarctica live vaccine strain in the tick vector Dermacentor variabilis

Background

The γ-proteobacterium Francisella tularensis is the etiologic agent of seasonal tick-transmitted tularemia epizootics in rodents and rabbits and of incidental infections in humans. The biology of F. tularensis in its tick vectors has not been fully described, particularly with respect to its quanta and duration of colonization, tissue dissemination, and transovarial transmission. A systematic study of the colonization of Dermacentor variabilis by the F. tularensis subsp. holarctica live vaccine strain (LVS) was undertaken to better understand whether D. variabilis may serve as an inter-epizootic reservoir for F. tularensis.

Methodology/Principal Findings

Colony-reared larva, nymph, and adult D. variabilis were artificially fed LVS via glass capillary tubes fitted over the tick mouthparts, and the level of colonization determined by microbial culture. Larvae and nymphs were initially colonized with 8.8±0.8×10¹ and 1.1±0.03×10³ CFU/tick, respectively. Post-molting, a significant increase in colonization of both molted nymphs and adults occurred, and LVS persisted in 42% of molted adult ticks at 126 days post-capillary tube feeding. In adult ticks, LVS initially colonized the gut, disseminated to hemolymph and salivary glands by 21 days, and persisted up to 165 days. LVS was detected in the salivary secretions of adult ticks after four days post intra-hemocoelic inoculation, and LVS recovered from salivary gland was infectious to mice with an infectious dose 50% of 3 CFU. LVS in gravid female ticks colonized via the intra-hemocoelic route disseminated to the ovaries and then to the oocytes, but the pathogen was not recovered from the subsequently-hatched larvae.

Conclusions/Significance

This study demonstrates that D. variabilis can be efficiently colonized with F. tularensis using artificial methods. The persistence of F. tularensis in D. variabilis suggests that this tick species may be involved in the maintenance of enzootic foci of tularemia in the central United States.     

Large Direct Repeats Flank Genomic Rearrangements between a New Clinical Isolate of Francisella tularensis subsp. tularensis A1 and Schu S4

Francisella tularensis subspecies tularensis consists of two separate populations A1 and A2. This report describes the complete genome sequence of NE061598, an F. tularensis subspecies tularensis A1 isolated in 1998 from a human with clinical disease in Nebraska, United States of America. The genome sequence was compared to Schu S4, an F. tularensis subspecies tularensis A1a strain originally isolated in Ohio in 1941. It was determined that there were 25 nucleotide polymorphisms (22 SNPs and 3 indels) between Schu S4 and NE061598; two of these polymorphisms were in potential virulence loci. Pulsed-field gel electrophoresis analysis demonstrated that NE061598 was an A1a genotype. Other differences included repeat sequences (n = 11 separate loci), four of which were contained in coding sequences, and an inversion and rearrangement probably mediated by insertion sequences and the previously identified direct repeats I, II, and III. Five new variable-number tandem repeats were identified; three of these five were unique in NE061598 compared to Schu S4. Importantly, there was no gene loss or gain identified between NE061598 and Schu S4. Interpretation of these data suggests there is significant sequence conservation and chromosomal synteny within the A1 population. Further studies are needed to determine the biological properties driving the selective pressure that maintains the chromosomal structure of this monomorphic pathogen.     

Comparative proteome analysis of cellular proteins extracted from highly virulent Francisella tularensis ssp. tularensis and less virulent F. tularensis ssp. holarctica and F. tularensis ssp. mediaasiatica

Francisella tularensis is the causative agent of the zoonotic disease tularemia. Four subspecies of this pathogen, namely ssp. tularensis, mediaasiatica, holarctica, and novicida are spread throughout the northern hemisphere. Although there are marked variations in their virulence to mammals, the subspecies are difficult to identify as they are closely genetically related. We carried out the comparative proteome analysis of cellular extracts from isolates representing the highly virulent subspecies tularensis, and the less virulent subspecies mediaasiatica and holarctica in order to identify new diagnostic markers and putative factors of virulence. We identified 27 protein spots that were either specifically present or at significantly higher abundance in ssp. tularensis strains, 22 proteins in ssp. mediaasiatica strains, and 26 proteins in ssp. holarctica strains. Subspecies tularensis-specific proteins might represent putative virulence factors. Of 27 identified tularensis-specific spots 17 represented charge and mass variants of proteins occurring in other subspecies, 7 spots were found to be present at higher abundance, and 3 spots were specifically present in tularensis strains. Amongst them, PilP protein, as a component necessary for the biogenesis of the type IV pilus, virulence and adhesion factor for many human pathogen, was identified. Furthermore, the identification of additional 27 proteins common for ssp. tularensis and mediaasiatica, and 19 proteins shared by ssp. mediaasiatica and holarctica documented apparent closer genetic similarity between ssp. tularensis and mediaasiatica.     

Deletion of ripA alleviates suppression of the inflammasome and MAPK by Francisella tularensis

Francisella tularensis is a facultative intracellular pathogen and potential biothreat agent. Evasion of the immune response contributes to the extraordinary virulence of this organism although the mechanism is unclear. Whereas wild-type strains induced low levels of cytokines, an F. tularensis ripA deletion mutant (LVSΔripA) provoked significant release of IL-1β, IL-18, and TNF-α by resting macrophages. IL-1β and IL-18 secretion was dependent on inflammasome components pyrin-caspase recruitment domain/apoptotic speck-containing protein with a caspase recruitment domain and caspase-1, and the TLR/IL-1R signaling molecule MyD88 was required for inflammatory cytokine synthesis. Complementation of LVSΔripA with a plasmid encoding ripA restored immune evasion. Similar findings were observed in a human monocytic line. The presence of ripA nearly eliminated activation of MAPKs including ERK1/2, JNK, and p38, and pharmacologic inhibitors of these three MAPKs reduced cytokine induction by LVSΔripA. Animals infected with LVSΔripA mounted a stronger IL-1β and TNF-α response than that of mice infected with wild-type live vaccine strain. This analysis revealed novel immune evasive mechanisms of F. tularensis.     

Viability and virulence of Francisella tularensis subsp. Holarctica in water ecosystems (experimental study)

Under conditions of artificial water biocenosis a virulent strain of F. tularensis could be detected in fresh water shrimps and mollusks for about a month, in Conepoda for up to 20 days and in Chydorus sphaericus for up to 7 days from the moment of the aquaria water contamination. In silt F. tularensis could be detected for a longer period (up to 2 months). Daphnia, Oligochaeta and C. sphaericus appeared to be unfavorable environment for this microorganism. The virulence level of F. tularensis microbial cells decreased in paralell with prolongation of their stay in water biocenosis. The presence of water biota favours F. tularensis preservation in water reservoirs for a longer time.