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.     

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.     

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.     

Structure and biosynthesis of free lipid A molecules that replace lipopolysaccharide in Francisella tularensis subsp. novicida

Francisella tularensis subsp. novicida U112 phospholipids, extracted without hydrolysis, consist mainly of phosphatidylethanolamine, phosphatidylglycerol, phosphatidylcholine, and two lipid A species, designated A1 and A2. These lipid A species, present in a ratio of 7:1, comprise 15% of the total phospholipids, as judged by 32Pi labeling. Although lipopolysaccharide is detectable in F. tularensis subsp. novicida U112, less than 5% of the total lipid A is covalently linked to it. A1 and A2 were analyzed by electrospray ionization and matrix-assisted laser desorption ionization mass spectrometry, gas chromatography/mass spectrometry, and NMR spectroscopy. Both compounds are disaccharides of glucosamine, acylated with primary 3-hydroxystearoyl chains at positions 2, 3, and 2' and a secondary palmitoyl residue at position 2'. Minor isobaric species and some lipid A molecules containing a 3-hydroxypalmitoyl chain in place of 3-hydroxystearate are also present. The 4'- and 3'-positions of A1 and A2 are not derivatized, and 3-deoxy-d-manno-octulosonic acid (Kdo) is not detectable. The 1-phosphate groups of both A1 and A2 are modified with an alpha-linked galactosamine residue, as shown by NMR spectroscopy and gas chromatography/mass spectrometry. An alpha-linked glucose moiety is attached to the 6'-position of A2. The lipid A released by mild acid hydrolysis of F. tularensis subsp. novicida lipopolysaccharide consists solely of component A1. F. tularensis subsp. novicida mutants lacking the arnT gene do not contain a galactosamine residue on their lipid A. Formation of free lipid A in F. tularensis subsp. novicida might be initiated by an unusual Kdo hydrolase present in the membranes of this organism.     

Comparative proteomic profiling of culture filtrate proteins of less and highly virulent Francisella tularensis strains

The facultative intracellular bacterium Francisella tularensis is the causal agent of the serious infectious disease tularemia. Despite the dynamic progress, which has been made in last few years, important questions regarding Francisella pathogenicity still remain to be answered. Generally, secreted proteins play an important role in pathogenicity of intracellular microbes. In this study, we investigated the protein composition of the culture filtrate proteins of highly virulent F. tularensis subsp. tularensis, strain SCHU S4 and attenuated F. tularensis subsp. holarctica, live vaccine strain using a comparative proteomic analysis. The majority of proteins identified in this study have been implicated in virulence mechanisms of other pathogens, and several have been categorized as having moonlighting properties; those that have more than one unrelated function. This profiling study of secreted proteins resulted in the unique detection of acid phosphatase (precursor) A (AcpA), β-lactamase, and hypothetical protein FTT0484 in the highly virulent strain SCHU S4 secretome. The release of AcpA may be of importance for F. tularensis subsp. tularensis virulence due to the recently described AcpA role in the F. tularensis escape from phagosomes.     

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.     

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).     

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.     

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.     

Proteomic analysis of antibody response in a case of laboratory-acquired infection with<Emphasis Type="Italic">Francisella tularensis</Emphasis> subsp.<Emphasis Type="Italic">tularensis</Emphasis>

Immunoproteomic analysis was applied to study the immunoreactivity of serum samples collected at different time points from a laboratory assistant accidentally infected with highly virulent strain of Francisella tularensis subsp. tularensis. Immunoblotting showed that the spectrum of F. tularensis antigens recognized specifically by immune sera remained with the exception for 1 antigen stable for up to 16 years after infection. Using immunoproteomics approach 10 immunoreactive antigens were successfully identified. Several new immunogenic F. tularensis proteins were described for the first time.     

Phase variations of Francisella tularensis lipopolysaccharide in human infection and immunization

The comparative study of the specificity of antibodies in human sera after tularemia infection and immunization with live tularemia infection was carried out with the use of passive hemagglutination and immunoblotting techniques. The sera of tularemia patients contained two different types of immunoglobulins: strictly specific to the antigenic epitopes of F. tularensis Iipopolysaccharide (LPS) and strictly specific to F. tularensis subsp. novicida LPS. Such phenomenon may be due to phase variations of the antigenic structure of F. tularensis LPS in the body of a slightly susceptible host. The immune sera of vaccinated were found to contain antibodies, strictly specific only to F. tularensis LPS. At the same time in one vaccinee by the presence of pronounced postvaccinal reactions was found sharply defined interaction between serum imunoglobulins and F. tularensis subsp. novicida LPS. As the result, the data on the possibility of the antigenic modification of F. tularensis in tularemia infection in humans were obtained. At the same time antigenic epitopes, characteristic of faintly pathogenic and closely related F. tularensis novicida LPS, appeared in the structure of F. tularensis LPS.     

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.     

Sensitivity spectrum of Francisella tularensis to antibiotics and synthetic antibacterial drugs

Sensitivity of 6 F. tularensis strains to 57 antibiotics and synthetic antibacterial drugs was studied. It was shown that the strains were highly sensitive to aminoglycosides, tetracyclines, anzamycins, quinolones, chloramphenicol, nitrofurantoin, nitroxoline, novobiocin and fusidin and resistant to penicillins, cephalosporins, polypeptides, vancomycin and sulfanylamides. The interrace differences in F. tularensis could be detected only by sensitivity to erythromycin, oleandomycin and spiramycin. There was observed no cross resistance to streptomycin and other aminoglycosides in F. tularensis. Assay of F. tularensis sensitivity to antibacterial drugs of various groups with the rapid photometric procedure and the agar diffusion method revealed that the results were highly comparable.     

A Francisella tularensis pathogenicity island required for intramacrophage growth

Francisella tularensis is a gram-negative, facultative intracellular pathogen that causes the highly infectious zoonotic disease tularemia. We have discovered a ca. 30-kb pathogenicity island of F. tularensis (FPI) that includes four large open reading frames (ORFs) of 2.5 to 3.9 kb and 13 ORFs of 1.5 kb or smaller. Previously, two small genes located near the center of the FPI were shown to be needed for intramacrophage growth. In this work we show that two of the large ORFs, located toward the ends of the FPI, are needed for virulence. Although most genes in the FPI encode proteins with amino acid sequences that are highly conserved between high- and low-virulence strains, one of the FPI genes is present in highly virulent type A F. tularensis, absent in moderately virulent type B F. tularensis, and altered in F. tularensis subsp. novicida, which is highly virulent for mice but avirulent for humans. The G+C content of a 17.7-kb stretch of the FPI is 26.6%, which is 6.6% below the average G+C content of the F. tularensis genome. This extremely low G+C content suggests that the DNA was imported from a microbe with a very low G+C-containing chromosome.     

Taurine: new implications for an old amino acid

We describe here a technique for allelic exchange in Francisella tularensis subsp. novicida utilizing polymerase chain reaction (PCR) products. Linear PCR fragments containing gene deletions with an erythromycin resistance cassette insertion were transformed into F. tularensis. The subsequent ErmR progeny were found to have undergone allelic exchange at the correct location in the genome; the minimum flanking homology necessary was 500 bp. This technique was used to create mglA, iglC, bla, and tul4 mutants in F. tularensis subsp. novicida strains. The mglA and iglC mutants were defective for intramacrophage growth, and the tul4 mutant lacked detectable Tul4 by Western immunoblot, as expected. Interestingly, the bla mutant maintained resistance to ampicillin, indicating the presence of multiple ampicillin resistance genes in F. tularensis.     

Genome characterisation of the genus Francisella reveals insight into similar evolutionary paths in pathogens of mammals and fish

ABSTRACT: BACKGROUND: Prior to this study, relatively few strains of Francisella had been genome-sequenced. Previously published Francisella genome sequences were largely restricted to the zoonotic agent F. tularensis. Only limited data were available for other members of the Francisella genus, including F. philomiragia, an opportunistic pathogen of humans, F. noatunensis, a serious pathogen of farmed fish, and other less well described endosymbiotic species. RESULTS: We determined the phylogenetic relationships of all known Francisella species, including some for which the phylogenetic positions were previously uncertain. The genus Francisella could be divided into two main genetic clades: one included F. tularensis, F. novicida, F. hispaniensis and Wolbachia persica, and another included F. philomiragia and F. noatunensis. Some Francisella species were found to have significant recombination frequencies. However, the fish pathogen F. noatunensis subsp. noatunensis was an exception due to it exhibiting a highly clonal population structure similar to the human pathogen F. tularensis. CONCLUSIONS: The genus Francisella can be divided into two main genetic clades occupying both terrestrial and marine habitats. However, our analyses suggest that the ancestral Francisella species originated in a marine habitat. The observed genome to genome variation in gene content and IS elements of different species supports the view that similar evolutionary paths of host adaptation developed independently in F. tularensis (infecting mammals) and F. noatunensis subsp. noatunensis (infecting fish).     

Direct isolation of Francisella spp. from environmental samples

Aims:  To develop a selective medium for isolation of F. tularensis, F. novicida and F. philomiragia from environmental samples.
Methods and Results:  A selective media, cysteine heart agar with 9% chocolatized sheep blood, containing polymyxin B, amphotericin B, cyclohexamide, cefepime and vancomycin (CHAB-PACCV) was developed and evaluated for growth of Francisella spp. No differences were observed in recovered colony forming units (CFUs) for F. tularensis , F. novicida and F. philomiragia on CHAB-PACCV vs nonselective CHAB. Growth of non- Francisella species was inhibited on CHAB-PACCV. When environmental samples were cultured on CHAB and CHAB-PACCV, only CHAB-PACCV allowed isolation of Francisella spp. Three new Francisella strains were isolated directly from seawater and seaweed samples by culture on CHAB-PACCV.
Conclusions:  CHAB-PACCV can be used for direct isolation of Francisella spp from environmental samples.
Significance and Impact of the Study:  Francisella spp. show a close association with environmental sources. Future utilization of CHAB-PACCV for isolation of Francisella spp. directly from environmental samples should prove valuable for investigating outbreaks and human infections attributed to environmental exposure.     

Identification of Francisella tularensis in natural water samples by PCR

Abstract Francisella tularensis , the causative agent of the epizootic disease tularemia in mammals, can be isolated from mud and water. To study the spread and persistence of Francisella tularensis in water, different strategies for pre-treatment of natural water samples prior to identification of the bacterium by polymerase chain reaction (PCR) were evaluated. A method for handling of samples taken from natural waters was developed. Applied on natural water samples amended with F. tularensis , the method rendered identification by PCR reproducible and it resulted in an amplified Francisella -specific product in all samples from natural waters tested. In addition, by employing primers targeting conserved regions of the 16S rDNA the presence of bacteria was demonstrated in all samples investigated. The results presented will, in combination with other techniques that allow identification, improve studies on the epizootiology and epidemiology of the genus Francisella .     

Independent Measurement of Ampicillin and Cloxacillin in Mixtures

A simplified method of assaying mixtures of ampicillin and cloxacillin is presented. The method consists of cylinder-plate assays of appropriately diluted samples in three assay systems in which: (i) cloxacillin is measured, utilizing agar seeded with a penicillinase-producing Staphylococcus aureus; (ii) ampicillin is measured after the greatest dilution practicable, by use of agar seeded with Sarcina lutea; and (iii) the proportion of activity found in the second assay that remains after treatment with staphylococcal penicillinase is determined on agar seeded with S. lutea. The error was 4% or less when each system was run in quadruplicate, except when there was 10 or more times as much cloxacillin as ampicillin, in which case the greatest error was 12%. Mixtures of methicillin with benzylpenicillin or with ampicillin can also be quantitated by this method.     

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.