Whole genome single nucleotide polymorphism based phylogeny of Francisella tularensis and its application to the development of a strain typing assay

BackgroundA low genetic diversity in Francisella tularensis has been documented. Current DNA based genotyping methods for typing F. tularensis offer a limited and varying degree of subspecies, clade and strain level discrimination power. Whole genome sequencing is the most accurate and reliable method to identify, type and determine phylogenetic relationships among strains of a species. However, lower cost typing schemes are necessary in order to enable typing of hundreds or even thousands of isolates.ResultsWe have generated a high-resolution phylogenetic tree from 40 Francisella isolates, including 13 F. tularensis subspecies holarctica (type B) strains, 26 F. tularensis subsp. tularensis (type A) strains and a single F. novicida strain. The tree was generated from global multi-strain single nucleotide polymorphism (SNP) data collected using a set of six Affymetrix GeneChip® resequencing arrays with the non-repetitive portion of LVS (type B) as the reference sequence complemented with unique sequences of SCHU S4 (type A). Global SNP based phylogenetic clustering was able to resolve all non-related strains. The phylogenetic tree was used to guide the selection of informative SNPs specific to major nodes in the tree for development of a genotyping assay for identification of F. tularensis subspecies and clades. We designed and validated an assay that uses these SNPs to accurately genotype 39 additional F. tularensis strains as type A (A1, A2, A1a or A1b) or type B (B1 or B2).ConclusionWhole-genome SNP based clustering was shown to accurately identify SNPs for differentiation of F. tularensis subspecies and clades, emphasizing the potential power and utility of this methodology for selecting SNPs for typing of F. tularensis to the strain level. Additionally, whole genome sequence based SNP information gained from a representative population of strains may be used to perform evolutionary or phylogenetic comparisons of strains, or selection of unique strains for whole-genome sequencing projects.

[1]  A. Benson,et al.  Genome diversity among regional populations of Francisella tularensis subspecies tularensis and Francisella tularensis subspecies holarctica isolated from the US. , 2004, FEMS microbiology letters.

[2]  P. Oyston.  Francisella tularensis: unravelling the secrets of an intracellular pathogen. , 2008, Journal of medical microbiology.

[3]  A. Vogler,et al.  An optimized, multiplexed multi‐locus variable‐number tandem repeat analysis system for genotyping Francisella tularensis , 2009, Letters in applied microbiology.

[4]  P. Keim,et al.  Francisella tularensis in the United States , 2005, Emerging infectious diseases.

[5]  P. Mead,et al.  Molecular Epidemiology of Francisella tularensis in the United States. , 2009, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[6]  Samuel I. Miller,et al.  Potential Source of Francisella tularensis Live Vaccine Strain Attenuation Determined by Genome Comparison , 2006, Infection and Immunity.

[7]  S Rozen,et al.  Primer3 on the WWW for general users and for biologist programmers. , 2000, Methods in molecular biology.

[8]  A. Sjöstedt,et al.  Worldwide Genetic Relationships among Francisella tularensis Isolates Determined by Multiple-Locus Variable-Number Tandem Repeat Analysis , 2004, Journal of bacteriology.

[9]  M. Forsman,et al.  Canonical Insertion-Deletion Markers for Rapid DNA Typing of Francisella tularensis , 2007, Emerging infectious diseases.

[10]  T. Cebula,et al.  Chips and SNPs, bugs and thugs: a molecular sleuthing perspective. , 2005, Journal of food protection.

[11]  J. Petersen,et al.  Tularemia: emergence/re-emergence. , 2005, Veterinary research.

[12]  Raymond K. Auerbach,et al.  Phylogeography of Francisella tularensis: Global Expansion of a Highly Fit Clone , 2009, Journal of bacteriology.

[13]  P. Mead,et al.  Epidemiologic and Molecular Analysis of Human Tularemia, United States, 1964–2004 , 2006, Emerging infectious diseases.

[14]  R. Fleischmann,et al.  A bioinformatic filter for improved base-call accuracy and polymorphism detection using the Affymetrix GeneChip® whole-genome resequencing platform , 2007, Nucleic Acids Research.

[15]  A. Yamada,et al.  Genetic diversity of Francisella tularensis subspecies holarctica strains isolated in Japan , 2008, Microbiology and immunology.

[16]  P. Keim,et al.  Molecular Epidemiology, Evolution, and Ecology of Francisella , 2007, Annals of the New York Academy of Sciences.

[17]  A. Sjöstedt,et al.  Discrimination of Human Pathogenic Subspecies of Francisella tularensis by Using Restriction Fragment Length Polymorphism , 2003, Journal of Clinical Microbiology.

[18]  M. Forsman,et al.  Evolution of Subspecies of Francisella tularensis , 2005, Journal of Bacteriology.

[19]  A. Sjöstedt,et al.  Evaluation of PCR-Based Methods for Discrimination ofFrancisella Species and Subspecies and Development of a Specific PCR That Distinguishes the Two Major Subspecies of Francisella tularensis , 2000, Journal of Clinical Microbiology.

[20]  Richard W Titball,et al.  Vaccination strategies for Francisella tularensis. , 2005, Advanced drug delivery reviews.

[21]  John P. Huelsenbeck,et al.  MrBayes 3: Bayesian phylogenetic inference under mixed models , 2003, Bioinform..

[22]  E. Ferri,et al.  Comparison of Different PCR Approaches for Typing of Francisella tularensis Strains , 2000, Journal of Clinical Microbiology.

[23]  John P. Huelsenbeck,et al.  MRBAYES: Bayesian inference of phylogenetic trees , 2001, Bioinform..

[24]  A. Sjöstedt,et al.  The development of tools for diagnosis of tularemia and typing of Francisella tularensis , 2004, APMIS : acta pathologica, microbiologica, et immunologica Scandinavica.

[25]  Jonathan P. Bollback,et al.  Bayesian Inference of Phylogeny and Its Impact on Evolutionary Biology , 2001, Science.