Worldwide genetic relationships among Francisella tularensis isolates determined by multiple-locus variable-number tandem repeat analysis

The intracellular bacterium Francisella tularensis is the causative agent of tularemia and poses a serious threat as an agent of bioterrorism. We have developed a highly effective molecular subtyping system from 25 variable-number tandem repeat (VNTR) loci. In our study, multiple-locus VNTR analysis (MLVA) was used to analyze genetic relationships and potential population structure within a global collection of 192 F. tularensis isolates, including representatives from each of the four subspecies. The VNTR loci displayed between 2 and 31 alleles with Nei's diversity values between 0.05 and 0.95. Neighbor-joining cluster analysis of VNTR data revealed 120 genotypes among the 192 F. tularensis isolates, including accurate subspecies identification. F. tularensis subsp. tularensis (type A) isolates showed great diversity at VNTR loci, while F. tularensis subsp. holarctica (type B) isolates showed much lower levels despite a much broader geographical prevalence. The resolution of two distinct clades within F. tularensis subsp. tularensis (designated A.I and A.II) revealed a previously unrecognized genetic division within this highly virulent subspecies. F. tularensis subsp. holarctica appears to have recently spread globally across continents from a single origin, while F. tularensis subsp. tularensis has a long and complex evolutionary history almost exclusively in North America. The sole non-North American type A isolates (Slovakian) were closely related to the SCHU S4 strain. Significant linkage disequilibrium was detected among VNTR loci of F. tularensis consistent with a clonal population structure. Overall, this work greatly augments the study of tularemia ecology and epidemiology, while providing a framework for future forensic analysis of F. tularensis isolates.     

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.     

Comparative characteristics of biological properties of Francisella tularensis strains isolated in the USSR]

With a large collection of the strains of F. tularensis isolated it has been recently shown that cultures belonging to holarctica and mediaasiatica circulate in the endemic foci of the USSR. By their biological and genetic properties the natural strains of F. tularensis were homogeneous and represented type cultures of F. tularensis. Various ecological conditions in the natural environment did not change within the last 20 years the sensitivity of the tularemia microbe to the antibacterial drugs.     

Influence of nutrient status and grazing pressure on the fate of Francisella tularensis in lake water

The natural reservoir of Francisella tularensis , the causative agent of tularaemia, is yet to be identified. We investigated the possibility that Francisella persists in natural aquatic ecosystems between outbreaks. It was hypothesized that nutrient-rich environments, with strong protozoan predation, favour the occurrence of the tularaemia bacterium. To investigate the differences in adaptation to aquatic environments of the species and subspecies of Francisella , we screened 23 strains for their ability to survive grazing by the ciliate Tetrahymena pyriformis . All the Francisella strains tested were consumed at a low rate, although significant differences between subspecies were found. The survival and virulence of gfp -labelled F. tularensis ssp. holarctica were then studied in a microcosm experiment using natural lake water, with varying food web complexities and nutrient availabilities. High nutrient conditions in combination with high abundances of nanoflagellates were found to favour F. tularensis ssp. holarctica . The bacterium was observed both free-living and within the cells of a nanoflagellate. Francisella tularensis entered a viable but nonculturable state during the microcosm experiment. When studied over a longer period of time, F. tularensis ssp. holarctica survived in the lake water, but loss of virulence was not prevented by either high nutrient availability or the presence of predators.     

Targeted inactivation of francisella tularensis genes by group II introns

Studies of the molecular mechanisms of pathogenesis of Francisella tularensis, the causative agent of tularemia, have been hampered by a lack of genetic techniques for rapid targeted gene disruption in the most virulent subspecies. Here we describe efficient targeted gene disruption in F. tularensis utilizing mobile group II introns (targetrons) specifically optimized for F. tularensis. Utilizing a targetron targeted to blaB, which encodes ampicillin resistance, we showed that the system works at high efficiency in three different subspecies: F. tularensis subsp. tularensis, F. tularensis subsp. holarctica, and "F. tularensis subsp. novicida." A targetron was also utilized to inactivate F. tularensis subsp. holarctica iglC, a gene required for virulence. The iglC gene is located within the Francisella pathogenicity island (FPI), which has been duplicated in the most virulent subspecies. Importantly, the iglC targetron targeted both copies simultaneously, resulting in a strain mutated in both iglC genes in a single step. This system will help illuminate the contributions of specific genes, and especially those within the FPI, to the pathogenesis of this poorly studied organism.     

Chromosome rearrangement and diversification of Francisella tularensis revealed by the type B (OSU18) genome sequence

The gamma-proteobacterium Francisella tularensis is one of the most infectious human pathogens, and the highly virulent organism F. tularensis subsp. tularensis (type A) and less virulent organism F. tularensis subsp. holarctica (type B) are most commonly associated with significant disease in humans and animals. Here we report the complete genome sequence and annotation for a low-passage type B strain (OSU18) isolated from a dead beaver found near Red Rock, Okla., in 1978. A comparison of the F. tularensis subsp. holarctica sequence with that of F. tularensis subsp. tularensis strain Schu4 (P. Larsson et al., Nat. Genet. 37:153-159, 2005) highlighted genetic differences that may underlie different pathogenicity phenotypes and the evolutionary relationship between type A and type B strains. Despite extensive DNA sequence identity, the most significant difference between type A and type B isolates is the striking amount of genomic rearrangement that exists between the strains. All but two rearrangements can be attributed to homologous recombination occurring between two prominent insertion elements, ISFtu1 and ISFtu2. Numerous pseudogenes have been found in the genomes and are likely contributors to the difference in virulence between the strains. In contrast, no rearrangements have been observed between the OSU18 genome and the genome of the type B live vaccine strain (LVS), and only 448 polymorphisms have been found within non-transposase-coding sequences whose homologs are intact in OSU18. Nonconservative differences between the two strains likely include the LVS attenuating mutation(s).     

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.     

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.     

Evolution of subspecies of Francisella tularensis

Analysis of unidirectional genomic deletion events and single nucleotide variations suggested that the four subspecies of Francisella tularensis have evolved by vertical descent. The analysis indicated an evolutionary scenario where the highly virulent F. tularensis subsp. tularensis (type A) appeared before the less virulent F. tularensis subsp. holarctica (type B). Compared to their virulent progenitors, attenuated strains of F. tularensis exhibited specific unidirectional gene losses.     

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.     

Population structure of Francisella tularensis

We have sequenced fragments of five metabolic housekeeping genes and two genes encoding outer membrane proteins from 81 isolates of Francisella tularensis, representing all four subspecies. Phylogenetic clustering of gene sequences from F. tularensis subsp. tularensis and F. tularensis subsp. holarctica aligned well with subspecies affiliations. In contrast, F. tularensis subsp. novicida and F. tularensis subsp. mediasiatica were indicated to be phylogenetically incoherent taxa. Incongruent gene trees and mosaic structures of housekeeping genes provided evidence for genetic recombination in F. tularensis.     

Phosphatase and penicillinase activities as stable traits for the differentiation of the racial classification of Francisella tularensis]

In the causative agent of tularemia new markers correlating with different subspecies of this microbe have been detected. Thus, F. tularensis strains belonging to the American and Central Asian subspecies are characterized by phosphatase activity, which makes it possible to use the phosphatase test for their differentiation from the strains of the holarctic variety. F. tularensis subsp. mediasiatica are incapable of producing beta-lactamase which differentiates them from the representatives of the varieties holarctica and tularensis. These newly discovered signs are stable and do not depend on the virulence of the cultures under study and on the conditions of the cultivation of F. tularensis.     

Vntr-genotyping of Francisella tularensis strains isolated in the former USSR territory and some European countries during epizootics in 1988 - 1989

Retrospective VNTR-analysis of 159 Francisella tularensis subsp. holarctica strains isolated in December 1988 - February 1989 in former USSR and some European countries was carried out. Analysis of heterogenic genotypes of strains allow to subdivide them into 30 groups of variants by individual genotypes, while cluster analysis--to subdivide them in 7 clusters with different number of compositions. The predominance of genotype C1 strains isolated on the Rostov and Archangelsk regions and the Crimea was established. F. tularensis strains isolated in winter time 1988 - 1989 in different geographic regions were supposed to be resident cultures typical for their biotope in natural focus of disease.     

Use of acid treatment and a selective medium to enhance the recovery of Francisella tularensis from water

Francisella tularensis has been associated with naturally occurring waterborne outbreaks and is also of interest as a potential biological weapon. Recovery of this pathogen from water using cultural methods is challenging due to the organism's fastidious growth requirements and interference by indigenous bacteria. A 15-min acid treatment procedure prior to culture on a selective agar was evaluated for recovery of F. tularensis from seeded water samples. Mean levels of reduction of virulent strains of F. tularensis subsp. holarctica and subsp. tularensis were less than 20% following acid treatment. The attenuated live vaccine strain (LVS) was less resistant to acid exposure. The acid treatment procedure coupled with plating on cystine heart agar with rabbit blood and antibiotics (CHARBab) allowed the isolation of F. tularensis seeded into five natural water samples.     

Development of multiple-locus variable-number tandem-repeat analysis for rapid genotyping of Ehrlichia ruminantium and its application to infected Amblyomma variegatum collected in heartwater endemic areas in Uganda

The rickettsial bacterium Ehrlichia ruminantium is the causative agent of heartwater, a serious tick-borne disease in ruminants. The genetic diversity of organisms in the field will have implications for cross-protective capacities of any vaccine developed, and for an effective vaccine design strategy proper genotyping and understanding of existing genetic diversity in the field is necessary. We searched for variable-number tandem-repeat (VNTR) loci for use in a multi-locus VNTR analysis (MLVA). Sequencing analysis of 30 potential VNTRs using a panel of 17 reference strains from geographically diverse origins identified 12 VNTRs with allelic profiles differing between strains. Application of MLVA to 38 E. ruminantium-infected Amblyomma variegatum collected from indigenous cattle in 6 different districts of Uganda identified 21 MLVA types. The discriminatory power of MLVA was greater than that of map1 PCR-restriction fragment length polymorphism analysis, with which only 6 genotypes were obtained. The high discriminatory power as well as cost- effective performance of MLVA provide the potential for this technique to be applied in the future with respect to optimizing vaccine trials by identifying local strain diversity, and also raise the possibility of exploring the association between E. ruminantium genotypes and phenotypes such as pathological outcome in the ruminant host.     

Multilocus VNTR-typing of Francisella tularensis strains

In the analysis of F. tularensis genome with the use of the specially developed program "DNA" a great number of loci containing tandem repeats were found. For analysis, 3 of them were selected and designated as FtA, FtB, FtC. The study of DNA of 40 F. tularensis strains in the polymerase chain reaction with specific primers to these loci a great variability in the number of repeats was established, the presence of 17 alleles being found in locus FtA, 5 alleles in locus FtB and 5 alleles in locus FtC. The strains under study formed 24 variants of genotypes, whose occurrence varied from 0.025 to 0.125. Taking into account the variability of the detected loci and a great number of potential loci VNTR in the genome, further development of this method will facilitate the creation of local and general data bases of the strains, thus ensuring more effective genetic typing of F. tularensis.     

Modulation of biogenesis of the Francisella tularensis subsp. novicida-containing phagosome in quiescent human macrophages and its maturation into a phagolysosome upon activation by IFN-gamma

Francisella tularensis is a highly virulent facultative intracellular pathogen that has been categorized as a class A bioterrorism agent, and is classified into four subsp, tularensis, holarctica, mediasiatica and novicida. Although the ability of F. tularensis subsp. novicida to cause tularemia in mice is similar to the virulent subsp. tularensis and holarctica, it is attenuated in humans. It is not known whether attenuation of F. tularensis subsp. novicida in humans is resulting from a different route of trafficking within human macrophages, compared with the tularensis or holarctica subsp. Here we show that in quiescent human monocytes-derived macrophages (hMDMs), the F. tularensis subsp. novicida containing phagosome (FCP) matures into a late endosome-like stage that acquires the late endosomal marker LAMP-2 but does not fuse to lysosomes. This modulation of phagosome biogenesis by F. tularensis is followed by disruption of the phagosome at 4-12 h and subsequent bacterial escape into cytoplasm where the organism replicates. In IFN-gamma-activated hMDMs, intracellular replication of F. tularensis is completely inhibited, and is associated with failure of the organism to escape from the phagosome into the cytoplasm for up to 24 h after infection. In IFN-gamma-activated hMDMs, the FCPs acquire the lysosomal enzymes Cathepsin D, which is excluded in quiescent hMDMs. When the lysosomes of IFN-gamma-activated hMDMs are preload with Texas Red Ovalbumin or BSA-gold, the FCPs acquire both lysosomal tracers. In contrast, both lysosomal tracers are excluded from the FCPs within quiescent hMDMs. We conclude that although F. tularensis subsp. novicida is attenuated in humans, it modulates biogenesis of its phagosome into a late endosome-like compartment followed by bacterial escape into the cytoplasm within quiescent hMDMs, similar to the virulent subsp. tularensis. In IFN-gamma-activated hMDMs, the organism fails to escape into the cytoplasm and its phagosome fuses to lysosomes, similar to inert particles.     

Multiple-locus variable-number of tandem-repeats analysis distinguishes Vibrio parahaemolyticus pandemic O3:K6 strains

A specific serotype of Vibrio parahaemolyticus, O3:K6, has recently been linked to epidemics of gastroenteritis in Southeast Asia, Japan, and North America. These pandemic O3:K6 strains appear to have recently spread across continents from a single origin to reach global coverage, based on profiling of strains by several molecular typing methods. In this study, variable-number tandem repeats (VNTR)-based fingerprinting was applied to clinical and environmental V. parahaemolyticus O3:K6 strains in an attempt to develop a molecular method with increased sensitivity for discriminating strains; the relative discriminatory powers were compared with ribotyping and pulsed-field gel electrophoresis (PFGE). All clinical strains tested were independent human isolates obtained from different outbreaks or from sporadic cases in Tokyo during the period from 1996 to 2003. Multiple-locus VNTR analysis (MLVA) was shown to have high resolution and reproducibility for typing of V. parahaemolyticus clones. MLVA analysis of 28 pandemic V. parahaemolyticus O3:K6 strains isolated from human cases produced 28 distinct VNTR patterns. The VNTR loci displayed between 2 and 15 alleles at each of eight loci with Nei's diversity index ranging from 0.35 and 0.91. These data demonstrated that MLVA is useful for individual strain typing of new O3:K6 strains, which appear to be closely related by other molecular methods.     

Development of a Multiple Loci Variable Number of Tandem Repeats Analysis (MLVA) to Unravel the Intra-Pathovar Structure of Pseudomonas syringae pv. actinidiae Populations Worldwide

The bacterial canker of kiwifruit by Pseudomonas syringae pv. actinidiae is an emblematic example of a catastrophic disease of fruit crops. In 2008 a new, extremely virulent form of the pathogen emerged and rapidly devastated many Actinidia spp. orchards all over the world. In order to understand differences in populations within this pathovar and to elucidate their diffusion and movements on world scale, it is necessary to be able to quickly and on a routine basis compare new isolates with previous records. In this report a worldwide collection of 142 strains was analyzed by MLVA, chosen as investigative technique for its efficacy, reproducibility, simplicity and low cost. A panel of 13 Variable Number of Tandem Repeats (VNTR) loci was identified and used to describe the pathogen population. The MLVA clustering is highly congruent with the population structure as previously established by other molecular approaches including whole genome sequencing and correlates with geographic origin, time of isolation and virulence. For convenience, we divided the VNTR loci in two panels. Panel 1 assay, using six loci, recognizes 23 different haplotypes, clustered into ten complexes with highest congruence with previous classifications. Panel 2, with seven VNTR loci, provides discriminatory power. Using the total set of 13 VNTR loci, 58 haplotypes can be distinguished. The recent hypervirulent type shows very limited diversity and includes, beside the strains from Europe, New Zealand and Chile, a few strains from Shaanxi, China. A broad genetic variability is observed in China, but different types are also retrievable in Japan and Korea. The low virulent strains cluster together and are very different from the other MLVA genotypes. Data were used to generate a public database in MLVAbank. MLVA represents a very promising first-line assay for large-scale routine genotyping, prior to whole genome sequencing of only the most relevant samples.     

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.