Search

The highly infectious and zoonotic pathogen Francisella tularensis is the etiologic agent of tularemia, a potentially fatal disease if untreated. Despite the high average nucleotide identity, which is >99.2% for the virulent subspecies and >98% for all four subspecies, including the opportunistic microbe Francisella tularensis subsp. novicida, there are considerable differences in genetic organization. These chromosomal disparities contribute to the substantial differences in virulence observed between the various F. tularensis subspecies and subtypes., The highly infectious and zoonotic pathogen Francisella tularensis is the etiologic agent of tularemia, a potentially fatal disease if untreated. Despite the high average nucleotide identity, which is >99.2% for the virulent subspecies and >98% for all four subspecies, including the opportunistic microbe Francisella tularensis subsp. novicida, there are considerable differences in genetic organization. These chromosomal disparities contribute to the substantial differences in virulence observed between the various F. tularensis subspecies and subtypes. The methods currently available to genotype F. tularensis cannot conclusively identify the associated subpopulation without using time-consuming testing or complex scoring matrices. To address this need, we developed both single and multiplex quantitative real-time PCR (qPCR) assays that can accurately detect and identify the hypervirulent F. tularensis subsp. tularensis subtype A.I, the virulent F. tularensis subsp. tularensis subtype A.II, F. tularensis subsp. holarctica (also referred to as type B), and F. tularensis subsp. mediasiatica, as well as opportunistic F. tularensis subsp. novicida from each other and near neighbors, such as Francisella philomiragia, Francisella persica, and Francisella-like endosymbionts found in ticks. These fluorescence-based singleplex and non-matrix scoring multiplex qPCR assays utilize a hydrolysis probe, providing sensitive and specific F. tularensis subspecies and subtype identification in a rapid manner. Furthermore, sequencing of the amplified F. tularensis targets provides clade confirmation and informative strain-specific details. Application of these qPCR- and sequencing-based detection assays will provide an improved capability for molecular typing and clinical diagnostics, as well as facilitate the accurate identification and differentiation of F. tularensis subpopulations during epidemiological investigations of tularemia source outbreaks.

Tularemia is a potentially fatal disease that is caused by the highly infectious and zoonotic pathogen Francisella tularensis . Despite the monomorphic nature of sequenced F. tularensis genomes, there is a significant degree of plasticity in the organization of genetic elements. The observed variability in these genomes is due primarily to the transposition of direct repeats and insertion sequence (IS) elements. Since current methods used to genotype F. tularensis are time-consuming and require extensive laboratory resources, IS elements were investigated as a means to subtype this organism. The unique spatial location of specific IS elements provided the basis for the development of a differential IS amplification (DISA) assay to detect and distinguish the more virulent F. tularensis subsp. tularensis (subtypes A.I and A.II) and subsp. holarctica (type B) strains from F. tularensis subsp. novicida and other near neighbors, including Francisella philomiragia and Francisella -like endosymbionts found in ticks. Amplicon sizes and sequences derived from DISA showed heterogeneity within members of the subtype A.I and A.II isolates but not the type B strains. These differences were due to a 312-bp fragment derived from the IS element IS Ftu1 . Analysis of wild-type F. tularensis isolates by DISA correlated with pulsed-field gel electrophoresis genotyping utilizing two different restriction endonucleases and provided rapid results with minimal sample processing. The applicability of this molecular typing assay for environmental studies was demonstrated by the accurate identification and differentiation of tick-borne F. tularensis . The described approach to IS targeting and amplification provides new capability for epidemiological investigations and characterizations of tularemia source outbreaks.

Francisella tularensis contains several highly pathogenic subspecies, including Francisella tularensis subsp. holarctica , whose distribution is circumpolar in the northern hemisphere. The phylogeography of these subspecies and their subclades was examined using whole-genome single nucleotide polymorphism (SNP) analysis, high-density microarray SNP genotyping, and real-time-PCR-based canonical SNP (canSNP) assays. Almost 30,000 SNPs were identified among 13 whole genomes for phylogenetic analysis. We selected 1,655 SNPs to genotype 95 isolates on a high-density microarray platform. Finally, 23 clade- and subclade-specific canSNPs were identified and used to genotype 496 isolates to establish global geographic genetic patterns. We confirm previous findings concerning the four subspecies and two Francisella tularensis subsp. tularensis subpopulations and identify additional structure within these groups. We identify 11 subclades within F. tularensis subsp. holarctica , including a new, genetically distinct subclade that appears intermediate between Japanese F. tularensis subsp. holarctica isolates and the common F. tularensis subsp. holarctica isolates associated with the radiation event (the B radiation) wherein this subspecies spread throughout the northern hemisphere. Phylogenetic analyses suggest a North American origin for this B-radiation clade and multiple dispersal events between North America and Eurasia. These findings indicate a complex transmission history for F. tularensis subsp. holarctica .

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.

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.