Search

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.

Background: Finland repeatedly reports some of the highest incidences of tularaemia worldwide. To determine genetic diversity of the aetiologic agent of tularaemia, Francisella tularensis, a total of 76 samples from humans (n = 15) and animals (n = 61) were analysed., Methods: We used CanSNPs and canINDEL hydrolysis or TaqMan MGB probes for the analyses, either directly from the clinical tissue samples (n = 21) or from bacterial isolates (n = 55)., Results: The genotypes of the strains were assigned to three previously described basal subspecies holarctica clades. The majority of strains (n = 67) were assigned to B.12, a clade reported to dominate in Scandinavia and Eastern Europe. A single strain was assigned to clade B.4, previously reported from North America, Europe and China. The remaining strains (n = 8) were members of clade B.6. Importantly, new diversity was discovered in clade B.6. We describe two newly designed TaqMan MGB probe assays for this new B.6 subclade B.70, and its previously identified sister clade B.11, a clade dominantly found in Western Europe., Conclusions: The high genetic diversity of F. tularensis subspecies holarctica present in Finland is consistent with previous findings in Sweden. The results suggest a northern and southern division of the B.6 subclade B.10, where B.11 predominates in Western and Central Europe and B.70 is found in Fennoscandia. Further research is required to define whether the vast diversity of genotypes found is related to different habitats or reservoir species, their different postglacial immigration routes to Fennoscandia, or dynamics of the reservoir species.