Francisella tularensis is one of the most infectious human pathogens known. In the past, both the former Soviet Union and the US had programs to develop weapons containing the bacterium. We report the complete genome sequence of a highly virulent isolate of F. tularensis (1,892,819 bp). The sequence uncovers previously uncharacterized genes encoding type IV pili, a surface polysaccharide and iron-acquisition systems. Several virulence-associated genes were located in a putative pathogenicity island, which was duplicated in the genome. More than 10% of the putative coding sequences contained insertion-deletion or substitution mutations and seemed to be deteriorating. The genome is rich in IS elements, including IS630 Tc-1 mariner family transposons, which are not expected in a prokaryote. We used a computational method for predicting metabolic pathways and found an unexpectedly high proportion of disrupted pathways, explaining the fastidious nutritional requirements of the bacterium. The loss of biosynthetic pathways indicates that F. tularensis is an obligate host-dependent bacterium in its natural life cycle. Our results have implications for our understanding of how highly virulent human pathogens evolve and will expedite strategies to combat them.

Francisella tularensis is a potent pathogen and a possible bioterrorism agent. Little is known, however, to explain the molecular basis for its virulence and the distinct differences in virulence found between the four recognized subspecies, F. tularensis subsp. tularensis, F. tularensis subsp. mediasiatica, F. tularensis subsp. holarctica, and F. tularensis subsp. novicida. We developed a DNA microarray based on 1,832 clones from a shotgun library used for sequencing of the highly virulent strain F. tularensis subsp. tularensis Schu S4. This allowed a genome-wide analysis of 27 strains representing all four subspecies. Overall, the microarray analysis confirmed a limited genetic variation within the species F. tularensis, and when the strains were compared, at most 3.7% of the probes showed differential hybridization. Cluster analysis of the hybridization data revealed that the causative agents of type A and type B tularemia, i.e., F. tularensis subsp. tularensis and F. tularensis subsp. holarctica, respectively, formed distinct clusters. Despite marked differences in their virulence and geographical origin, a high degree of genomic similarity between strains of F. tularensis subsp. tularensis and F. tularensis subsp. mediasiatica was apparent. Strains from Japan clustered separately, as did strains of F. tularensis subsp. novicida. Eight regions of difference (RD) 0.6 to 11.5 kb in size, altogether comprising 21 open reading frames, were identified that distinguished strains of the moderately virulent subspecies F. tularensis subsp. holarctica and the highly virulent subspecies F. tularensis subsp. tularensis. One of these regions, RD1, allowed for the first time the development of an F. tularensis-specific PCR assay that discriminates each of the four subspecies.

No abstract