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.

In recent years, an increasing diversity of species has been recognized within the family Francisellaceae. Unfortunately, novel isolates are sometimes misnamed in initial publications or multiple sources propose different nomenclature for genetically highly similar isolates. Thus, unstructured and occasionally incorrect information can lead to confusion in the scientific community. Historically, detection of Francisella tularensis in environmental samples has been challenging due to the considerable and unknown genetic diversity within the family, which can result in false positive results. We have assembled a comprehensive collection of genome sequences representing most known Francisellaceae species/strains and restructured them according to a taxonomy that is based on phylogenetic structure. From this structured dataset, we identified a small number of genomic regions unique to F. tularensis that are putatively suitable for specific detection of this pathogen in environmental samples. We designed and validated specific PCR assays based on these genetic regions that can be used for the detection of F. tularensis in environmental samples, such as water and air filters.

In Germany tularemia is a re-emerging zoonotic disease. Therefore, we investigated wild animals and environmental water samples for the presence and phylogenetic diversity of Francisella tularensis in the poorly studied Berlin/Brandenburg region. The phylogenomic analysis of three isolates from wild animals revealed three new subclades within the phylogenetic tree of F. tularensis [B.71 from a raccoon dog (Nyctereutes procyonoides); B.74 from a red fox (Vulpes vulpes), and B.75 from a Eurasian beaver (Castor fiber albicus)]. The results from histological, PCR, and genomic investigations on the dead beaver showed that the animal suffered from a systemic infection. Indications were found that the bacteria were released from the beaver carcass into the surrounding environment. We demonstrated unexpectedly high and novel phylogenetic diversity of F. tularensis in Germany and the fact that the bacteria persist in the environment for at least one climatic season. These findings support a broader host species diversity than previously known regarding Germany. Our data further support the assumption derived from previous serological studies of an underestimated frequency of occurrence of the pathogen in the environment and in wild animals. F. tularensis was isolated from animal species not previously reported as natural hosts in Germany.

BackgroundTularemia is a zoonotic disease caused by Francisella tularensis that has been found in many different vertebrates. In Germany most human infections are caused by contact with infected European brown hares (Lepus europaeus). The aim of this study was to elucidate the epidemiology of tularemia in hares using phenotypic and genotypic characteristics of F. tularensis.ResultsCultivation of F. tularensis subsp. holarctica bacteria from organ material was successful in 31 of 52 hares that had a positive PCR result targeting the Ft-M19 locus. 17 isolates were sensitive to erythromycin and 14 were resistant. Analysis of VNTR loci (Ft-M3, Ft-M6 and Ft-M24), INDELs (Ftind33, Ftind38, Ftind49, RD23) and SNPs (B.17, B.18, B.19, and B.20) was shown to be useful to investigate the genetic relatedness of Francisella strains in this set of strains. The 14 erythromycin resistant isolates were assigned to clade B.I, and 16 erythromycin sensitive isolates to clade B.IV and one isolate was found to belong to clade B.II. MALDI-TOF mass spectrometry (MS) was useful to discriminate strains to the subspecies level.ConclusionsF. tularensis seems to be a re-emerging pathogen in Germany. The pathogen can easily be identified using PCR assays. Isolates can also be identified within one hour using MALDI-TOF MS in laboratories where specific PCR assays are not established. Further analysis of strains requires genotyping tools. The results from this study indicate a geographical segregation of the phylogenetic clade B.I and B.IV, where B.I strains localize primarily within eastern Germany and B.IV strains within western Germany. This phylogeographical pattern coincides with the distribution of biovar I (erythromycin sensitive) and biovar II (erythromycin resistance) strains. When time and costs are limiting parameters small numbers of isolates can be analysed using PCR assays combined with DNA sequencing with a focus on genetic loci that are most likely discriminatory among strains found in a specific area. In perspective, whole genome data will have to be investigated especially when terrorist attack strains need to be tracked to their genetic and geographical sources.

All documentation and Python-based source code for the CanSNPer are freely available at http://github.com/adrlar/CanSNPer.