Previous Article
Journal of Bacteriology, June 2005, p. 3903-3908, Vol. 187, No. 11
0021-9193/05/$08.00+0 doi:10.1128/JB.187.11.3903-3908.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Evolution of Subspecies of Francisella tularensis
Kerstin Svensson,1
Pär Larsson,1,2
Daniel Johansson,1
Mona Byström,1
Mats Forsman,1 and
Anders Johansson1,3*
Department of NBC-Analysis, Swedish Defence Research Agency, SE-901 82 Umeå, Sweden,1
Department of Clinical Microbiology, Bacteriology,2
Department of Clinical Microbiology, Infectious Diseases, Umeå University, SE-901 85 Umeå, Sweden3
Received 26 November 2004/
Accepted 23 February 2005

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

TEXT
Francisella tularensis, the causative agent of tularemia, is
considered a potential bioterrorist agent. Important premises
are an extremely low infectious dose and a potential for airborne
transmission (
6). The two most clinically important entities
of tularemia, type A and type B, correspond to the highly virulent
subspecies
F. tularensis subsp.
tularensis and the moderately
virulent
F. tularensis subsp.
holarctica, respectively (
18,
24). Taxonomic work has identified two additional subspecies,
F. tularensis subsp.
mediasiatica, exhibiting a moderate virulence,
and "
F. tularensis subsp.
novicida," with a low virulence in
animals and humans (
24). The latter subspecies has less-fastidious
extracellular growth requirements than the other subspecies
and a distinct lipopolysaccharide O-antigen (
10,
29).
Based on small subunit RNA sequences, F. tularensis is classified as a member of the
-subgroup of proteobacteria (24). The two species F. tularensis and F. philomiragia and in addition a number of more recently identified tick endosymbionts are the only members of the genus Francisella, which diverges deeply among the
-proteobacteria (22, 27).
Assuming that F. tularensis subspecies have evolved from a common ancestor, identification of genetic differences among subspecies might provide insights into species phylogeny and provides a basis for studies of more functional issues. Although exhibiting differences in virulence, geographical distribution, and a few biochemical tests, F. tularensis subspecies are highly similar in gene content (1). Mapping of genetic differences will allow future exploration of their relationships to functional correlates. In a recent DNA-based study using highly mutable variable-number tandem repeat (VNTR) sequences dispersed over the genome, F. tularensis subsp. tularensis (type A) showed more diversity than F. tularensis subsp. holarctica (type B), suggesting the former subspecies to be evolutionarily older (11). Significant linkage disequilibrium was detected among VNTR loci of F. tularensis, consistent with a predominantly clonal population structure.
The present work was based on findings from a whole-genome microarray study of multiple F. tularensis strains (1). The microarray study identified large size regions of difference (RDs) among F. tularensis strains. Notably, a pair of direct repeated sequences flanked seven of eight identified RDs. We hypothesized that the presence of direct repeats in a genomic region would represent a propensity for interstrain variation of that region. Therefore, direct repeated sequences were searched for in the genome sequence of strain SCHU S4, and when identified, the corresponding genomic regions were assayed for variability among various strains of F. tularensis. In parallel, single nucleotide variations (SNVs) were analyzed by sequencing of internal fragments of genes of F. tularensis and related bacteria of the
-subgroup of proteobacteria.
Bacteria and identification of genomic regions with direct repeats.
Criteria for including F. tularensis strains in this study were to ensure maximum spatial and temporal diversity. Forty-five isolates (Table 1) representing the four F. tularensis subspecies were selected from the Francisella Strain Collection (FSC) in Umeå, Sweden. Bacteria were grown on modified Thayer-Martin agar. Cells were suspended into phosphate-buffered saline and heat killed. DNA was prepared using silica and guanidine isothiocyanate buffer as described previously (23).
A Perl script (detailed in the supplemental text) was devised
to search for direct repeats in the, at that time, unfinished
genome sequence of
F. tularensis strain SCHU S4 (FSC237). The
completed genome sequence is assigned accession no.
AJ749949 in GenBank (
13). More than 70 genomic regions flanked by direct
repeats were evaluated by PCR in a subset of five
F. tularensis isolates that represented each of the four
F. tularensis subspecies
and the Japanese
F. tularensis subsp.
holarctica variant. Genomic
regions that exhibited size polymorphism on agarose gels were
further analyzed in 45 strains (Table
1). PCR amplicons were
sequenced to identify the junctions of the deletion/insertion
events. Primer sequences for PCR amplification were posted in
the supplemental material in Table S1.
Use of unidirectional deletions for phylogenetic analysis.
The use of RDs for phylogenetic analysis relies on an assumption of unidirectional deletion events that eventually become fixed in bacterial populations (phylogenetic lineages) (2, 14). An RD was used for depicting evolution only when it was found to represent a deletion and not an insertion. Three criteria were used for excluding RDs from use in phylogenetic analysis. (i) The first criterion was whether the corresponding region in the SCHU S4 genome exhibited three or more closely localized direct repeats. We reasoned that the presence of more than two repeats might indicate an evolutionary history of genome amplification. Therefore, the SCHU S4 genome was scrutinized for repeats present in the vicinity of RDs (±10,000 bp). For a review of direct repeat mediated excision or amplification in bacteria, see references 19 and 20. (ii) The second criterion was whether the RD was flanked by multicopy insertion sequence elements. This criterion was based on an assumption that the genomic region is potentially prone to recombination events within and among taxa in phylogenetic lineages. (iii) The third criterion was whether the RD showed polymorphism within a subspecies.
SNVs and calculation of phylogenetic trees.
A subset of 15 F. tularensis strains was selected for analysis of SNVs (Table 1). Additionally, DNA was isolated from a freeze-dried ampoule of the type strain of the tick endosymbiont Wolbachia persica (ATCC VR-331) as described previously (23). A set of seven genes was selected, previously used in multiple locus sequence typing of various bacteria: aroA, parC, pgm, tpiA, trpE, atpA, and uup (28). Internal fragments of the genes were PCR amplified. Primer sequences for PCR amplification were posted in supplemental material in Table S1. Amplification was successful for seven gene fragments from each F. tularensis strain and four gene fragments from W. persica (aroA, parC, tpiA, and atpA). PCR amplicons were sequenced on both strands using the amplification primers and standard protocol on an ABI 377 platform (PE Applied Biosystems). Sequences of gene orthologs in the genomes of Yersinia pestis (CO92) and Agrobacterium tumefaciens (C58) were retrieved from GenBank. Phylogenetic analyses were performed using concatenated nucleotide sequences, individual gene fragments, and concatenated in-frame-translated amino acid sequences. Maximum likelihood trees were reconstructed from DNA sequences by using PAUP v. 4.0b10 (25) as described previously (7). The extent of pairwise congruence between individual gene trees was estimated by a method comparing the likelihood values calculated for single gene trees to the values for other gene trees and to 200 randomly generated trees of the same size (7). Congruence of gene trees indicates a common evolutionary history of investigated genes. Phylogenetic reconstruction from amino acid sequences was performed using PHYLIP (v. 3.6b). Phylogenetic methods are further detailed in the supplemental text.
Analysis of RDs among F. tularensis strains.
PCR analysis of the genomes of 45 F. tularensis isolates identified 17 RDs in F. tularensis (Table 1). Generally, PCR amplification showed the presence of one of two alternate fragment sizes for each RD and isolate. The pattern of amplification allowed an immediate subspecies identification. With very few exceptions, all isolates of a subspecies showed the same fragment size for each RD. A polymorphism was present at RD21-22 among isolates of a subspecies and at RD18-19 among variants of a strain.
F. tularensis genomes showing alternate PCR fragment sizes were sequenced across each RD. Alignment of sequences for each RD showed the smaller-sized PCR fragments to contain only one direct repeat and to be missing a large genomic segment compared to the SCHU S4 sequence. The boundary of the missing sequence was defined by the two direct repeats (Fig. 1). Nucleotide degeneracy of flanking repeats allowed further analysis. We found each single repeat present in a smaller PCR fragment to be a composite repeat derived from fusion of a left and right flanking repeat (Fig. 1; also Fig. S1 in the supplemental material). The finding is consistent with a deletion mechanism mediated by recombination of flanking repeats followed by excision of the intervening sequence. The deletion will be irreversible if two direct repeats were originally present and a single composite repeat is left after excision. As an alternative, two direct repeats may mediate genome amplification of intervening sequence as demonstrated for RD22 (Table 1). At this locus, sequence analysis showed duplication (12 isolates) or triplication (3 isolates) of the region flanked by repeats (sequence data not shown).
Genes were truncated in isolates with smaller-sized PCR fragments
at 16/17 RDs (Table
2; also Table S2 in the supplemental material).
Only RD1c represents an intergenic sequence. RD11 and RD16-22
were identified and characterized de novo in this study and
found to affect 12 open reading frames in the
F. tularensis SCHU S4 genome (Table
2). Gene content of RD1-7 has been previously
described (
1) and was further characterized in this study (Table
S2 in the supplemental material). Among deleted genes, some
require discussion with regard to virulence properties of various
isolates. The gene
pdpD of RD6 was recently suggested to play
a role in
F. tularensis virulence (
15). In concordance to previous
analysis of a smaller number of isolates, we confirmed the gene
to be deleted in all type B isolates including the live vaccine
strain (
1,
15). Genes in RD18 and RD19 discriminated among isolates
of individual subspecies. Genes in RD19 represent type IV pili
building block proteins that were absent in, e.g., the live
vaccine strain of
F. tularensis. Homologous genes have been
implied as virulence determinants in other bacteria (
5). Finally,
genes in RD18 represent genes of a novel protein family unique
to
F. tularensis (
13). The latter genes were found to be deleted
in an avirulent phenotype of strain SCHU (FSC043) and in the
live vaccine strain of
F. tularensis.
A deletion-based phylogeny of F. tularensis.
Nine RDs fulfilled our strict selection criteria for a unidirectional
genetic event (see above) and were included in the deletion-based
phylogenetic analysis (Table
1 and Fig.
2). An evolutionary
scenario is suggested where successive losses of genetic material
at RDs parallel the taxonomic units of
F. tularensis at the
subspecies level. The analysis proposes an evolution of
F. tularensis where the highly virulent
F. tularensis subsp.
tularensis (type
A) preceded the appearance of the less virulent
F. tularensis subsp.
holarctica (type B). Our findings partially contrast
to the evolution of
Y. pestis that apparently became more virulent
while losing genetic material in the transition from its ancestor,
Yersinia pseudotuberculosis (
3). The evolutionary transitions
of
F. tularensis subspecies (
F. tularensis subsp.
novicida,
subsp.
tularensis, and subsp.
holarctica) are more complex and
apparently involve both acquisition and loss of virulence since
F. tularensis subsp.
novicida exhibits the lowest virulence
and was found to be evolutionarily oldest. Work on other pathogens
has shown that the evolution of virulence, in particular the
adaptation to living within eukaryotic hosts, may involve both
acquisition and loss of genes (
16), a fact that may explain
the observations.
All nine RDs were present in all
F. tularensis subsp.
tularensis isolates. The deletion data show isolates of
F. tularensis subsp.
mediasiatica to be very similar to isolates of
F. tularensis subsp.
tularensis. Only RD1c showed a deletion event unique
for
F. tularensis subsp.
mediasiatica isolates. Therefore, the
branching order of
F. tularensis subsp.
mediasiatica and
F. tularensis subsp.
tularensis remains unclear. Isolates from
Japan of
F. tularensis subsp.
holarctica showed deletion at
six RDs. Two additional deletions at RDs were found in Eurasian
and North American
F. tularensis subsp.
holarctica isolates.
This suggests that Japanese isolates are evolutionary intermediates
between the highly virulent type A (
F. tularensis subsp.
tularensis)
isolates and the less virulent Eurasian/North American type
B isolates (
F. tularensis subsp.
holarctica). The single isolate
of
F. tularensis subsp.
novicida was chosen for rooting, based
on the most extensive number of SNVs in seven gene fragments
(see below). The deletion analysis per se could not place
F. tularensis subsp.
novicida due to failure of PCR amplification
of three RDs (RD5, RD7, and RD16). Notably, none of 17 RDs showed
evidence of deletion events in
F. tularensis subsp.
novicida.
The evolutionary direction might also be inferred from analysis
of the gene content of RDs. The finding that 16/17 deletions
truncate conserved genes (genes with homologs in other bacteria)
implies deletions but not insertions have occurred during evolution.
(Table
2; also Table S2 in the supplemental material). Because
deletions are expected to represent unidirectional genetic events,
RDs of
F. tularensis are suggested to be valuable evolutionary
markers. In fact, if all RDs with a uniform pattern within a
subspecies were used in phylogenetic analysis the depicted evolutionary
scenario would not change (data not shown).
A phylogeny inferred from analysis of SNVs.
Compared to other genomic characters, SNVs exhibit slow mutation rates, making them valuable for phylogenetic analysis (9, 17). We analyzed seven housekeeping gene fragments (3,135 bp) in 15 F. tularensis isolates and four gene fragments (1,830 bp) of the tick endosymbiont W. persica. A total of 78 SNVs was detected among 15 F. tularensis isolates (see Fig. S2 in the supplemental material). The single isolate of F. tularensis subsp. novicida was the most divergent, showing 41 unique SNVs. Overall, the sample of worldwide F. tularensis isolates exhibited a low level of average pairwise sequence diversity (
) with values in the range of 0.0026 to 0.0101 for individual genes. The pattern of SNVs allowed a straightforward subspecies recognition, each F. tularensis subspecies being identified by several fixed polymorphisms. There was an extensive sequence similarity among isolates of a F. tularensis subspecies.
Four-gene phylogenetic trees were computed to describe the phylogenetic position of F. tularensis and its four subspecies among the gamma proteobacteria (Fig. 3). To root the phylogeny, we used A. tumefaciens as outgroup, as the
-subgroup is considered ancient among gamma proteobacteria (8). Y. pestis of the
-subgroup of proteobacteria was included as a second outgroup. In the phylogeny obtained, the tick endosymbiont W. persica is the first and deepest branch, thus representing an early divergence from a common Francisella ancestor. Among the four F. tularensis subspecies, F. tularensis subsp. novicida appeared to represent the first branching. The result is in agreement to findings of a greater amount of genetic divergence among tick endosymbionts or isolates of F. tularensis subsp. novicida compared to isolates of three other F. tularensis subspecies (4, 12, 22). A greater diversity among isolates likely represents a longer evolutionary history. The four-gene trees demonstrated the split between F. tularensis subsp. tularensis/mediasiatica on one hand and F. tularensis subsp. holarctica isolates from Japan/Eurasia/North America on the other hand. Seven-gene trees illustrated the position of Japanese F. tularensis subsp. holarctica isolates close to but distinct from Eurasian/North American F. tularensis subsp. holarctica isolates (Fig. 3). The five F. tularensis subsp. tularensis isolates were grouped into two distinct clades. Overall, the present SNV and deletion-based phylogenies support the genetic relationships of F. tularensis subspecies calculated from VNTR data, including a genetic division of type A isolates into two distinct clades recently designated A.I and A.II, respectively (Fig. 2 and Fig. 3) (11).
The pairwise statistical comparisons of congruence between gene
trees indicated that recombination across
F. tularensis subspecies
occurs rarely or not at all and justifies the use of concatenated
gene sequences for tree reconstruction (data not shown). Based
on the pattern of deletions and SNVs among isolates, we conclude
that the distribution of fixed characters unique to
F. tularensis subspecies reflects a genetic separation of subspecies and also
that the subspecies have descended clonally.
Deletions that are specific to laboratory strains.
Cell populations within a single F. tularensis strain might carry unique unidirectional deletions, i.e., at RD18-19 (Table 1). The finding of unique deletion characters among bacterial colonies derived from a single bacterial seed stock suggests that regions flanked by direct repeats are prone to mutation in F. tularensis genomes. Under natural evolution, it appears that direct repeat-mediated deletions only rarely became fixed, as illustrated in our analysis across F. tularensis subspecies. According to our results, the majority of observed deletions might be regarded as evolutionary footprints in F. tularensis genomes spanning over wide evolutionary time scales. In contrast, it seems that laboratory culture on artificial media facilitates the fixation of deletion variants. We suggest that deletion at RD18-19 might represent genetic events that occurred at the laboratory during passage or storage of F. tularensis strains. This is supported by the finding of a mix of two distinct cell populations that exhibit these RDs in laboratory stocks of strains FSC074 (RD19) and FSC158 (RD18) (Table 1). Based on sequence analysis, direct repeat-mediated excisions at RD18 and RD19 have occurred also in the genome of the live vaccine strain of F. tularensis. The live vaccine strain was derived during repeated passages on artificial media (26). Thus, it seems plausible that the fixation of a deletion variant (the live vaccine strain) was facilitated by a repeated selection of single colonies during laboratory passages. It remains to be determined if genes located in RD18 and/or RD19 have a role in the attenuation of the live vaccine strain, the Russian vaccine strain 015 (FSC338), and the avirulent phenotype of the SCHU strain (FSC043). All virulence-attenuated strains included in this study showed deletion at RD18 and/or RD19.
Nucleotide sequence accession numbers.
All sequences reported in this paper have been deposited in the GenBank database with accession no. AY794406 to AY794434 (RD2 to RD7, RD11, RD16-22) and AY794435 to AY794543 (SNVs in seven gene fragments). Sequences of RD1 are assigned accession no. AF469614 to AF469619 (1).

ACKNOWLEDGMENTS
This work was supported by funding from the Swedish MoD, project
no. A4854, the Swedish Medical Research Council, Formas, Sweden,
and the Medical Faculty, Umeå University, Umeå,
Sweden.
We are indebted to numerous scientific colleagues for kindly providing strains to the Francisella Strain Collection in Umeå, Sweden.

FOOTNOTES
* Corresponding author. Mailing address: NBC-Analysis, Swedish Defence Research Agency, SE-901 82 Umeå, Sweden. Phone: 46 90 10 66 16. Fax: 46 90 133006. E-mail:
anders.johansson{at}infdis.umu.se.

Supplemental material for this article may be found at http://jb.asm.org/. 

REFERENCES
1 - Broekhuijsen, M., P. Larsson, A. Johansson, M. Byström, U. Eriksson, E. Larsson, R. G. Prior, A. Sjöstedt, R. W. Titball, and M. Forsman. 2003. Genome-wide DNA microarray analysis of Francisella tularensis strains demonstrates extensive genetic conservation within the species but identifies regions that are unique to the highly virulent F. tularensis subsp. tularensis. J. Clin. Microbiol. 41:2924-2931.[Abstract/Free Full Text]
2 - Brosch, R., S. V. Gordon, M. Marmiesse, P. Brodin, C. Buchrieser, K. Eiglmeier, T. Garnier, C. Gutierrez, G. Hewinson, K. Kremer, L. M. Parsons, A. S. Pym, S. Samper, D. van Soolingen, and S. T. Cole. 2002. A new evolutionary scenario for the Mycobacterium tuberculosis complex. Proc. Natl. Acad. Sci. USA 99:3684-3689.[Abstract/Free Full Text]
3 - Chain, P. S., E. Carniel, F. W. Larimer, J. Lamerdin, P. O. Stoutland, W. M. Regala, A. M. Georgescu, L. M. Vergez, M. L. Land, V. L. Motin, R. R. Brubaker, J. Fowler, J. Hinnebusch, M. Marceau, C. Medigue, M. Simonet, V. Chenal-Francisque, B. Souza, D. Dacheux, J. M. Elliott, A. Derbise, L. J. Hauser, and E. Garcia. 2004. Insights into the evolution of Yersinia pestis through whole-genome comparison with Yersinia pseudotuberculosis. Proc. Natl. Acad. Sci. USA 101:13826-13831.[Abstract/Free Full Text]
4 - Clarridge, J. E., III, T. J. Raich, A. Sjöstedt, G. Sandström, R. O. Darouiche, R. M. Shawar, P. R. Georghiou, C. Osting, and L. Vo. 1996. Characterization of two unusual clinically significant Francisella strains. J. Clin. Microbiol. 34:1995-2000.[Abstract]
5 - Craig, L., M. E. Pique, and J. A. Tainer. 2004. Type IV pilus structure and bacterial pathogenicity. Nat. Rev. Microbiol. 2:363-378.[CrossRef][Medline]
6 - Dennis, D. T., T. V. Inglesby, D. A. Henderson, J. G. Bartlett, M. S. Ascher, E. Eitzen, A. D. Fine, A. M. Friedlander, J. Hauer, M. Layton, S. R. Lillibridge, J. E. McDade, M. T. Osterholm, T. O'Toole, G. Parker, T. M. Perl, P. K. Russell, and K. Tonat. 2001. Tularemia as a biological weapon: medical and public health management. JAMA 285:2763-2773.[Abstract/Free Full Text]
7 - Feil, E. J., E. C. Holmes, D. E. Bessen, M. S. Chan, N. P. Day, M. C. Enright, R. Goldstein, D. W. Hood, A. Kalia, C. E. Moore, J. Zhou, and B. G. Spratt. 2001. Recombination within natural populations of pathogenic bacteria: short-term empirical estimates and long-term phylogenetic consequences. Proc. Natl. Acad. Sci. USA 98:182-187.[Abstract/Free Full Text]
8 - Gupta, R., and E. Griffiths. 2002. Critical issues in bacterial phylogeny. Theor. Popul. Biol. 61:423.[CrossRef][Medline]
9 - Gutacker, M. M., J. C. Smoot, C. A. Migliaccio, S. M. Ricklefs, S. Hua, D. V. Cousins, E. A. Graviss, E. Shashkina, B. N. Kreiswirth, and J. M. Musser. 2002. Genome-wide analysis of synonymous single nucleotide polymorphisms in Mycobacterium tuberculosis complex organisms: resolution of genetic relationships among closely related microbial strains. Genetics 162:1533-1543.[Abstract/Free Full Text]
10 - Hollis, D. G., R. E. Weaver, A. G. Steigerwalt, J. D. Wenger, C. W. Moss, and D. J. Brenner. 1989. Francisella philomiragia comb. nov. (formerly Yersinia philomiragia) and Francisella tularensis biogroup novicida (formerly Francisella novicida) associated with human disease. J. Clin. Microbiol. 27:1601-1608.[Abstract/Free Full Text]
11 - Johansson, A., J. Farlow, P. Larsson, M. Dukerich, E. Chambers, M. Byström, J. Fox, M. C. Chu, M. Forsman, A. Sjöstedt, and P. Keim. 2004. Worldwide genetic relationships among Francisella tularensis isolates determined by multiple-locus variable-number tandem repeat analysis. J. Bacteriol. 186:5808-5818.[Abstract/Free Full Text]
12 - Johansson, A., A. Ibrahim, I. Göransson, U. Eriksson, D. Gurycova, J. E. Clarridge III, and A. Sjöstedt. 2000. Evaluation of PCR-based methods for discrimination of Francisella species and subspecies and development of a specific PCR that distinguishes the two major subspecies of Francisella tularensis. J. Clin. Microbiol. 38:4180-4185.[Abstract/Free Full Text]
13 - Larsson, P., P. C. Oyston, P. Chain, M. C. Chu, M. Duffield, H. H. Fuxelius, E. Garcia, G. Halltorp, D. Johansson, K. E. Isherwood, P. D. Karp, E. Larsson, Y. Liu, S. Michell, J. Prior, R. Prior, S. Malfatti, A. Sjöstedt, K. Svensson, N. Thompson, L. Vergez, J. K. Wagg, B. W. Wren, L. E. Lindler, S. G. Andersson, M. Forsman, and R. W. Titball. 2005. The complete genome sequence of Francisella tularensis, the causative agent of tularemia. Nat. Genet. 37:153-159.[CrossRef][Medline]
14 - Mostowy, S., D. Cousins, J. Brinkman, A. Aranaz, and M. A. Behr. 2002. Genomic deletions suggest a phylogeny for the Mycobacterium tuberculosis complex. J. Infect. Dis. 186:74-80.[CrossRef][Medline]
15 - Nano, F. E., N. Zhang, S. C. Cowley, K. E. Klose, K. K. Cheung, M. J. Roberts, J. S. Ludu, G. W. Letendre, A. I. Meierovics, G. Stephens, and K. L. Elkins. 2004. A Francisella tularensis pathogenicity island required for intramacrophage growth. J. Bacteriol. 186:6430-6436.[Abstract/Free Full Text]
16 - Ochman, H., and N. A. Moran. 2001. Genes lost and genes found: evolution of bacterial pathogenesis and symbiosis. Science 292:1096-1099.[Abstract/Free Full Text]
17 - Pearson, T., J. D. Busch, J. Ravel, T. D. Read, S. D. Rhoton, J. M. U'Ren, T. S. Simonson, S. M. Kachur, R. R. Leadem, M. L. Cardon, M. N. Van Ert, L. Y. Huynh, C. M. Fraser, and P. Keim. 2004. Phylogenetic discovery bias in Bacillus anthracis using single-nucleotide polymorphisms from whole-genome sequencing. Proc. Natl. Acad. Sci. USA 101:13536-13541.[Abstract/Free Full Text]
18 - Penn, R. L. 2005. Francisella tularensis (Tularemia), p. 2674-2685. In G. L. Mandell, J. E. Bennet, and R. Dolin (ed.), Mandell, Douglas and Bennett's principles and practice of infectious diseases, 6th ed., vol. 2. Churchill Livingstone, Ltd., Edinburgh, Scotland.
19 - Rocha, E. P. 2003. An appraisal of the potential for illegitimate recombination in bacterial genomes and its consequences: from duplications to genome reduction. Genome Res. 13:1123-1132.[Abstract/Free Full Text]
20 - Romero, D., and R. Palacios. 1997. Gene amplification and genomic plasticity in prokaryotes. Annu. Rev. Genet. 31:91-111.[CrossRef][Medline]
21 - Samrakandi, M. M., C. Zhang, M. Zhang, J. Nietfeldt, J. Kim, P. C. Iwen, M. E. Olson, P. D. Fey, G. E. Duhamel, S. H. Hinrichs, J. D. Cirillo, and A. K. Benson. 2004. Genome diversity among regional populations of Francisella tularensis subspecies tularensis and Francisella tularensis subspecies holarctica isolated from the US. FEMS Microbiol. Lett. 237:9-17.[CrossRef][Medline]
22 - Scoles, G. A. 2004. Phylogenetic analysis of the Francisella-like endosymbionts of Dermacentor ticks. J. Med. Entomol. 41:277-286.[Medline]
23 - Sjöstedt, A., U. Eriksson, L. Berglund, and A. Tärnvik. 1997. Detection of Francisella tularensis in ulcers of patients with tularemia by PCR. J. Clin. Microbiol. 35:1045-1048.[Abstract]
24 - Sjöstedt, A. B. 2005. Francisella, p. 200-210. In D. J. Brenner, N. R. Krieg, J. T. Staley, and G. M. Garrity (ed.), The Proteobacteria, part B. Bergey's manual of systematic bacteriology, 2nd ed., vol. 2. Springer-Verlag, New York, N.Y.
25 - Swofford, D. L. 1998. PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. Sinauer Associates Inc., Sunderland, Mass.
26 - Tigertt, W. D. 1962. Soviet viable Pasteurella tularensis vaccines. Bacteriol. Rev. 26:354-373.[Free Full Text]
27 - Titball, R. W., A. Johansson, and M. Forsman. 2003. Will the enigma of Francisella tularensis virulence soon be solved? Trends Microbiol. 11:118-123.[CrossRef][Medline]
28 - Urwin, R., and M. C. Maiden. 2003. Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol. 11:479-487.[CrossRef][Medline]
29 - Vinogradov, E., W. J. Conlan, J. S. Gunn, and M. B. Perry. 2004. Characterization of the lipopolysaccharide O-antigen of Francisella novicida (U112). Carbohydr. Res. 339:649-654.[CrossRef][Medline]
Journal of Bacteriology, June 2005, p. 3903-3908, Vol. 187, No. 11
0021-9193/05/$08.00+0 doi:10.1128/JB.187.11.3903-3908.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
This article has been cited by other articles:
-
Lindgren, H., Honn, M., Golovlev, I., Kadzhaev, K., Conlan, W., Sjostedt, A.
(2009). The 58-Kilodalton Major Virulence Factor of Francisella tularensis Is Required for Efficient Utilization of Iron. Infect. Immun.
77: 4429-4436
[Abstract]
[Full Text]
-
Salomonsson, E., Kuoppa, K., Forslund, A.-L., Zingmark, C., Golovliov, I., Sjostedt, A., Noppa, L., Forsberg, A.
(2009). Reintroduction of Two Deleted Virulence Loci Restores Full Virulence to the Live Vaccine Strain of Francisella tularensis. Infect. Immun.
77: 3424-3431
[Abstract]
[Full Text]
-
Salomonsson, E., Forsberg, A., Roos, N., Holz, C., Maier, B., Koomey, M., Winther-Larsen, H. C.
(2009). Functional analyses of pilin-like proteins from Francisella tularensis: complementation of type IV pilus phenotypes in Neisseria gonorrhoeae. Microbiology
155: 2546-2559
[Abstract]
[Full Text]
-
Vogler, A. J., Birdsell, D., Price, L. B., Bowers, J. R., Beckstrom-Sternberg, S. M., Auerbach, R. K., Beckstrom-Sternberg, J. S., Johansson, A., Clare, A., Buchhagen, J. L., Petersen, J. M., Pearson, T., Vaissaire, J., Dempsey, M. P., Foxall, P., Engelthaler, D. M., Wagner, D. M., Keim, P.
(2009). Phylogeography of Francisella tularensis: Global Expansion of a Highly Fit Clone. J. Bacteriol.
191: 2474-2484
[Abstract]
[Full Text]
-
Gallagher, L. A., McKevitt, M., Ramage, E. R., Manoil, C.
(2008). Genetic Dissection of the Francisella novicida Restriction Barrier. J. Bacteriol.
190: 7830-7837
[Abstract]
[Full Text]
-
La Scola, B., Elkarkouri, K., Li, W., Wahab, T., Fournous, G., Rolain, J.-M., Biswas, S., Drancourt, M., Robert, C., Audic, S., Lofdahl, S., Raoult, D.
(2008). Rapid comparative genomic analysis for clinical microbiology: The Francisella tularensis paradigm. Genome Res
18: 742-750
[Abstract]
[Full Text]
-
Molins-Schneekloth, C. R., Belisle, J. T., Petersen, J. M.
(2008). Genomic Markers for Differentiation of Francisella tularensis subsp. tularensis A.I and A.II Strains. Appl. Environ. Microbiol.
74: 336-341
[Abstract]
[Full Text]
-
Dempsey, M. P., Dobson, M., Zhang, C., Zhang, M., Lion, C., Gutierrez-Martin, C. B., Iwen, P. C., Fey, P. D., Olson, M. E., Niemeyer, D., Francesconi, S., Crawford, R., Stanley, M., Rhodes, J., Wagner, D. M., Vogler, A. J., Birdsell, D., Keim, P., Johansson, A., Hinrichs, S. H., Benson, A. K.
(2007). Genomic Deletion Marking an Emerging Subclone of Francisella tularensis subsp. holarctica in France and the Iberian Peninsula. Appl. Environ. Microbiol.
73: 7465-7470
[Abstract]
[Full Text]
-
Guina, T., Radulovic, D., Bahrami, A. J., Bolton, D. L., Rohmer, L., Jones-Isaac, K. A., Chen, J., Gallagher, L. A., Gallis, B., Ryu, S., Taylor, G. K., Brittnacher, M. J., Manoil, C., Goodlett, D. R.
(2007). MglA Regulates Francisella tularensis subsp. novicida (Francisella novicida) Response to Starvation and Oxidative Stress. J. Bacteriol.
189: 6580-6586
[Abstract]
[Full Text]
-
Thomas, R. M., Titball, R. W., Oyston, P. C. F., Griffin, K., Waters, E., Hitchen, P. G., Michell, S. L., Grice, I. D., Wilson, J. C., Prior, J. L.
(2007). The Immunologically Distinct O Antigens from Francisella tularensis Subspecies tularensis and Francisella novicida Are both Virulence Determinants and Protective Antigens. Infect. Immun.
75: 371-378
[Abstract]
[Full Text]
-
Hajjar, A. M., Harvey, M. D., Shaffer, S. A., Goodlett, D. R., Sjostedt, A., Edebro, H., Forsman, M., Bystrom, M., Pelletier, M., Wilson, C. B., Miller, S. I., Skerrett, S. J., Ernst, R. K.
(2006). Lack of In Vitro and In Vivo Recognition of Francisella tularensis Subspecies Lipopolysaccharide by Toll-Like Receptors. Infect. Immun.
74: 6730-6738
[Abstract]
[Full Text]
-
Rohmer, L., Brittnacher, M., Svensson, K., Buckley, D., Haugen, E., Zhou, Y., Chang, J., Levy, R., Hayden, H., Forsman, M., Olson, M., Johansson, A., Kaul, R., Miller, S. I.
(2006). Potential Source of Francisella tularensis Live Vaccine Strain Attenuation Determined by Genome Comparison. Infect. Immun.
74: 6895-6906
[Abstract]
[Full Text]
-
Ben Nasr, A., Haithcoat, J., Masterson, J. E., Gunn, J. S., Eaves-Pyles, T., Klimpel, G. R.
(2006). Critical role for serum opsonins and complement receptors CR3 (CD11b/CD18) and CR4 (CD11c/CD18) in phagocytosis of Francisella tularensis by human dendritic cells (DC): uptake of Francisella leads to activation of immature DC and intracellular survival of the bacteria. J. Leukoc. Biol.
80: 774-786
[Abstract]
[Full Text]
-
Dempsey, M. P., Nietfeldt, J., Ravel, J., Hinrichs, S., Crawford, R., Benson, A. K.
(2006). Paired-End Sequence Mapping Detects Extensive Genomic Rearrangement and Translocation during Divergence of Francisella tularensis subsp. tularensis and Francisella tularensis subsp. holarctica Populations.. J. Bacteriol.
188: 5904-5914
[Abstract]
[Full Text]
-
Nubel, U., Reissbrodt, R., Weller, A., Grunow, R., Porsch-Ozcurumez, M., Tomaso, H., Hofer, E., Splettstoesser, W., Finke, E.-J., Tschape, H., Witte, W.
(2006). Population Structure of Francisella tularensis.. J. Bacteriol.
188: 5319-5324
[Abstract]
[Full Text]
-
Twine, S., Bystrom, M., Chen, W., Forsman, M., Golovliov, I., Johansson, A., Kelly, J., Lindgren, H., Svensson, K., Zingmark, C., Conlan, W., Sjostedt, A.
(2005). A Mutant of Francisella tularensis Strain SCHU S4 Lacking the Ability To Express a 58-Kilodalton Protein Is Attenuated for Virulence and Is an Effective Live Vaccine. Infect. Immun.
73: 8345-8352
[Abstract]
[Full Text]
-
Bystrom, M., Bocher, S., Magnusson, A., Prag, J., Johansson, A.
(2005). Tularemia in Denmark: Identification of a Francisella tularensis subsp. holarctica Strain by Real-Time PCR and High-Resolution Typing by Multiple-Locus Variable-Number Tandem Repeat Analysis. J. Clin. Microbiol.
43: 5355-5358
[Abstract]
[Full Text]