For many infections transmitting to humans from reservoirs in nature, disease dispersal patterns over space and time are largely unknown. Here, a reversed genomics approach helped us understand disease dispersal and yielded insight into evolution and biological properties of Francisella tularensis, the bacterium causing tularemia. We whole-genome sequenced 67 strains and characterized by single-nucleotide polymorphism assays 138 strains, collected from individuals infected 1947-2012 across Western Europe. We used the data for phylogenetic, population genetic and geographical network analyses. All strains (n=205) belonged to a monophyletic population of recent ancestry not found outside Western Europe. Most strains (n=195) throughout the study area were assigned to a star-like phylogenetic pattern indicating that colonization of Western Europe occurred via clonal expansion. In the East of the study area, strains were more diverse, consistent with a founder population spreading from east to west. The relationship of genetic and geographic distance within the F. tularensis population was complex and indicated multiple long-distance dispersal events. Mutation rate estimates based on year of isolation indicated null rates; in outbreak hotspots only, there was a rate of 0.4 mutations/genome/year. Patterns of nucleotide substitution showed marked AT mutational bias suggestive of genetic drift. These results demonstrate that tularemia has moved from east to west in Europe and that F. tularensis has a biology characterized by long-range geographical dispersal events and mostly slow, but variable, replication rates. The results indicate that mutation-driven evolution, a resting survival phase, genetic drift and long-distance geographical dispersal events have interacted to generate genetic diversity within this species.
This is an open access article under the terms of the Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution and reproduction in any medium, provided the original author and source are credited.
AfsetJ. E.,
LarssenK. W.,
BerghK.,
LärkerydA.,
SjödinA.,
JohanssonA.,
ForsmanM.2015; Phylogeographical pattern of Francisella tularensis in a nationwide outbreak of tularaemia in Norway, 2011. Euro Surveill 20:9–14 [View Article][PubMed]
BirdsellD. N.,
PearsonT.,
PriceE. P.,
HornstraH. M.,
NeraR. D.,
StoneN.,
GruendikeJ.,
KaufmanE. L.,
PettusA. H. et al.2012; Melt analysis of mismatch amplification mutation assays (Melt-MAMA): a functional study of a cost-effective SNP genotyping assay in bacterial models. PLoS One 7:e32866 [View Article][PubMed]
BurrowsS. M.,
ElbertW.,
LawrenceM. G.,
PöschlU.2009; Bacteria in the global atmosphere – Part 1: Review and synthesis of literature data for different ecosystems. Atmos Chem Phys Discuss 9:10777–10827 [View Article]
CuiY.,
YuC.,
YanY.,
LiD.,
LiY.,
JombartT.,
WeinertL. A.,
WangZ.,
GuoZ. et al.2013; Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis
. Proc Natl Acad Sci U S A 110:577–582 [View Article][PubMed]
DempseyM. P.,
DobsonM.,
ZhangC.,
ZhangM.,
LionC.,
Gutiérrez-MartínC. B.,
IwenP. C.,
FeyP. D.,
OlsonM. E. et al.2007; Genomic deletion marking an emerging subclone of Francisella tularensis subsp. holarctica in France and the Iberian Peninsula. Appl Environ Microbiol 73:7465–7470 [View Article][PubMed]
DennisD. T.,
InglesbyT.V,
HendersonD. A.,
BartlettJ. G.,
AscherM. S.,
EitzenE.,
FineA. D.,
FriedlanderA. M.,
HauerJ. et al.2001; Tularemia as a biological weapon: medical and public health management. JAMA 285:2763–2773[PubMed]
DrummondA. J.,
SuchardM. A.,
XieD.,
RambautA.2012; Bayesian phylogenetics with BEAUti and the beast 1.7. Mol Biol Evol 29:1969–1973 [View Article][PubMed]
FerrettiM.,
PaciG.,
PorriniS.,
GalardiL.,
BagliaccaM.2010; Habitat use and home range traits of resident and relocated hares (Lepus europaeus, Pallas). Ital J Anim Sci 9:e54 [View Article]
GutiérrezM. P.,
BratosM. A.,
GarroteJ.I,
DueñasA.,
AlmarazA.,
AlamoR.,
Rodríguez MarcosH.,
Rodríguez RecioM. J.,
MuñozM. F. et al.2003; Serologic evidence of human infection by Francisella tularensis in the population of Castilla y León (Spain) prior to 1997. FEMS Immunol Med Microbiol 35:165–169 [View Article][PubMed]
HedgeJ.,
WilsonD. J.2014; Bacterial phylogenetic reconstruction from whole genomes is robust to recombination but demographic inference is not. MBio 5:e02158 [View Article][PubMed]
JohanssonA.,
LärkerydA.,
WiderströmM.,
MörtbergS.,
MyrtännäsK.,
ÖhrmanC.,
BirdsellD.,
KeimP.,
WagnerD. M. et al.2014; An outbreak of respiratory tularemia caused by diverse clones of Francisella tularensis
. Clin Infect Dis 59:1546–1553 [View Article][PubMed]
JusatzH. J.1952a; [Second report on the propagation of tularemia into middle and western Europe in present time; geomedical investigations on the development in the last decade and epidemiological prognosis]. Z Hyg Infektionskr 134:350–374[PubMed]
JusatzH. J.1961; The geographical distribution of tularemia throughout the world, 1911–1959. In Welt-Suchen Atlas pp. 7–12 Edited by
RodenwaldtE.
Hamburg: Falk-Verlag;
Low-DécarieE.,
KolberM.,
HommeP.,
LofanoA.,
DumbrellA.,
GonzalezA.,
BellG.2015; Community rescue in experimental metacommunities. Proc Natl Acad Sci U S A 112:14307–14312 [View Article][PubMed]
MoranN. A.,
McCutcheonJ. P.,
NakabachiA.2008; Genomics and evolution of heritable bacterial symbionts. Annu Rev Genet 42:165–190 [View Article][PubMed]
NemergutD. R.,
SchmidtS. K.,
FukamiT.,
O'NeillS. P.,
BilinskiT. M.,
StanishL. F.,
KnelmanJ. E.,
DarcyJ. L.,
LynchR. C. et al.2013; Patterns and processes of microbial community assembly. Microbiol Mol Biol Rev 77:342–356 [View Article][PubMed]
NguyenT. M.,
IlefD.,
JarraudS.,
RouilL.,
CampeseC.,
CheD.,
HaeghebaertS.,
GaniayreF.,
MarcelF. et al.2006; A community-wide outbreak of legionnaires disease linked to industrial cooling towers-how far can contaminated aerosols spread?. J Infect Dis 193:102–111 [View Article][PubMed]
PearsonT.,
BuschJ. D.,
RavelJ.,
ReadT. D.,
RhotonS. D.,
U'RenJ. M.,
SimonsonT. S.,
KachurS. M.,
LeademR. R. et al.2004; Phylogenetic discovery bias in Bacillus anthracis using single-nucleotide polymorphisms from whole-genome sequencing. Proc Natl Acad Sci U S A 101:13536–13541 [View Article][PubMed]
PiloP.,
JohanssonA.,
FreyJ.2009; Identification of Francisella tularensis cluster in central and western Europe. Emerg Infect Dis 15:2049–2051 [View Article][PubMed]
RametteA.,
TiedjeJ. M.2007; Multiscale responses of microbial life to spatial distance and environmental heterogeneity in a patchy ecosystem. Proc Natl Acad Sci U S A 104:2761–2766 [View Article][PubMed]
RinaldiA.,
CervioC.,
FrittoliM.,
MandelliG.1964; Descrizione di un focolaio di tularemia in Italia (nota preliminare). [Description of an outbreak of tularaemia (preliminary note)]. Sel. Vet 5:353–363
RochaE. P.,
SmithJ. M.,
HurstL. D.,
HoldenM. T.,
CooperJ. E.,
SmithN. H.,
FeilE. J.2006; Comparisons of dN/dS are time dependent for closely related bacterial genomes. J Theor Biol 239:226–235 [View Article][PubMed]
SimpsonJ. T.,
WongK.,
JackmanS. D.,
ScheinJ. E.,
JonesS. J.,
BirolI.2009; ABySS: a parallel assembler for short read sequence data. Genome Res 19:1117–1123 [View Article][PubMed]
SmithA. T.,
JohnstonC. H.2008; Lepus europaeus 2008. Available online: www.iucnredlist.org [Accessed date on 20 November 2016]
SmithD. J.,
TimonenH. J.,
JaffeD. A.,
GriffinD. W.,
BirmeleM. N.,
PerryK. D.,
WardP. D.,
RobertsM. S.2013; Intercontinental dispersal of bacteria and archaea by transpacific winds. Appl Environ Microbiol 79:1134–1139 [View Article][PubMed]
VoglerA. J.,
BirdsellD.,
PriceL. B.,
BowersJ. R.,
Beckstrom-SternbergS. M.,
AuerbachR. K.,
Beckstrom-SternbergJ. S.,
JohanssonA.,
ClareA. et al.2009; Phylogeography of Francisella tularensis: global expansion of a highly fit clone. J Bacteriol 191:2474–2484 [View Article][PubMed]
WangX.,
Van NostrandJ. D.,
DengY.,
LuX.,
WangC.,
ZhouJ.,
HanX.2015; Scale-dependent effects of climate and geographic distance on bacterial diversity patterns across northern China's grasslands. FEMS Microbiol Ecol 91:fiv133 [View Article]
Antwerpen, M. H., Schacht, E., Kaysser, P. & Splettstoesser, W. D. Complete Genome Sequence of a Francisella tularensis subsp. holarctica Strain from Germany Causing Lethal Infection in Common Marmosets. Genome Announc, 1. GenBank Accession #: PRJNA175244 2013
Atkins, L. M., Holder, M. E., Ajami, N. J., Metcalf, G. A., Weissenberger, G. M., Wang, M., Vee, V., Han, Y., Muzny, D. M., Gibbs, R. A. & other authors. High-Quality Draft Genome Sequence of Francisella tularensis subsp. holarctica Strain OR96-0246. Genome Announc. 3. GenBank Accession #: PRJNA30669 2015
Barabote, R. D., Xie, G., Brettin, T. S., Hinrichs, S. H., Fey, P. D., Jay, J. J., Engle, J. L., Godbole, S. D., Noronha, J. M., Scheuermann, R. H. & other authors. Complete genome sequence of Francisella tularensis subspecies holarctica FTNF002-00. PLoS One, 4, e7041. GenBank Accession #: PRJNA20197 2009
Baylor College of Medicine (BCM), USA. GenBank Accession #s: PRJNA30633 & PRJNA30635 2008
Coolen, J. P., Sjödin, A., Maraha, B., Hajer, G. F., Forsman, M., Verspui, E., Frenay, H. M., Notermans, D. W., de Vries, M. C., Reubsaet, F. A. & other authors. Draft genome sequence of Francisella tularensis subsp. holarctica BD11-00177. Stand Genomic Sci, 8, 539-47. GenBank Accession #: PRJNA177784 2013
Johansson, A., Lärkeryd, A., Widerström, M., Mörtberg, S., Myrtennäs, K., Öhrman, C., Birdsell, D., Keim, P., Wagner, D. M., Forsman, M. & other authors. An outbreak of respiratory tularemia caused by diverse clones of Francisella tularensis. Clin Infect Dis, 59, 1546-53. EMBL Nucleotide Sequence Database (ENA) Accession #s: ERS353713 & ERS353729 2014
Karlsson, E., Svensson, K., Lindgren, P., Byström, M., Sjödin, A., Forsman, M. & Johansson, A. The phylogeographic pattern of Francisella tularensis in Sweden indicates a Scandinavian origin of Eurosiberian tularaemia. Environmental Microbiology, 15, 634-645. GenBank Accession #: PRJNA89145 2013
La Scola, B., Elkarkouri, K., Li, W., Wahab, T., Fournous, G., Rolain, J. M., Biswas, S., Drancourt, M., Robert, C., Audic, S. & other authors. Rapid comparative genomic analysis for clinical microbiology: the Francisella tularensis paradigm. Genome Res, 18, 742-50. GenBank Accession #: PRJNA19645 2008
Petrosino, J. F., Xiang, Q., Karpathy, S. E., Jiang, H., Yerrapragada, S., Liu, Y., Gioia, J., Hemphill, L., Gonzalez, A., Raghavan, T. M. & other authors. Chromosome rearrangement and diversification of Francisella tularensis revealed by the type B (OSU18) genome sequence. J Bacteriol, 188, 6977-85. GenBank Accession #: PRJNA17265 2006
Sjödin, A., Svensson, K., Öhrman, C., Ahlinder, J., Lindgren, P., Duodu, S., Johansson, A., Colquhoun, D. J., Larsson, P. & Forsman, M. Genome characterisation of the genus Francisella reveals insight into similar evolutionary paths in pathogens of mammals and fish. BMC Genomics, 13, 268. GenBank Accession #: PRJNA73369 2012
Svensson, K., Sjödin, A., Byström, M., Granberg, M., Brittnacher, M. J., Rohmer, L., Jacobs, M. A., Sims-Day, E. H., Levy, R., Zhou, Y. & other authors. Genome sequence of Francisella tularensis subspecies holarctica strain FSC200, isolated from a child with tularemia. J Bacteriol, 194, 6965-6. GenBank Accession #: PRJNA16087 2012
Swedish Defence Research Agency (FOI), Sweden. GenBank Accession #s: PRJNA285145 & PRJNA355394 2016
Wittwer, M. Spiez laboratory, Switzerland. GenBank Accession #: PRJNA286987 2016