Umeå University's logo

Advances in typing methodologies have recently reformed the field of molecular epidemiology of pathogens. The falling cost of sequencing technologies is creating a deluge of whole genome sequencing data that burdens bioinformatics resources and tool development. In particular, single nucleotide polymorphisms in core genomes of pathogens are recognized as the most important markers for inferring genetic relationships because they are evolutionarily stable and amenable to high-throughput detection methods. Sequence data will provide an excellent opportunity to extend our understanding of infectious disease when the challenge of extracting knowledge from available sequence resources is met. Here, we present an efficient and user-friendly genotype classification pipeline, CanSNPer, based on an easily expandable database of predefined canonical single nucleotide polymorphisms.

We used whole-genome analysis and subsequent characterization of geographically diverse strains using new genetic signatures to identify distinct subgroups within Francisella tularensis subsp. tularensis group A.I: A.I.3, A.I.8, and A.I.12. These subgroups exhibit complex phylogeographic patterns within North America. The widest distribution was observed for A.I.12, which suggests an adaptive advantage.

Previous studies of the causative agent of tularaemia, Francisella tularensis have identified phylogeographic patterns suggestive of environmental maintenance reservoirs. To investigate the phylogeography of tularaemia in Sweden, we selected 163 clinical isolates obtained during 1995-2009 in 10 counties and sequenced one isolate's genome to identify new genetic markers. An improved typing scheme based on two indels and nine SNPs was developed using hydrolysis or TaqMan MGB probe assays. The results showed that much of the known global genetic diversity of F. tularensis subsp. holarctica is present in Sweden. Thirteen of the 163 isolates belonged to a new genetic group that is basal to all other known members of the major genetic clade B.I, which is spread across the Eurosiberian region. One hundred and twenty-five of the 163 Swedish isolates belonged to B.I, but individual clades' frequencies differed from county to county (P < 0.001). Subsequent analyses revealed a correlation between genotype variation over time and recurrent outbreaks at specific places, supporting the 'maintenance reservoir' environmental maintenance hypothesis. Most importantly, the findings reveal the presence of diverse source populations of F. tularensis subsp. holarctica in Sweden and suggest a historical spread of the disease from Scandinavia to other parts of Eurosiberia.

Background: Prior to this study, relatively few strains of Francisella had been genome-sequenced. Previously published Francisella genome sequences were largely restricted to the zoonotic agent F. tularensis. Only limited data were available for other members of the Francisella genus, including F. philomiragia, an opportunistic pathogen of humans, F. noatunensis, a serious pathogen of farmed fish, and other less well described endosymbiotic species. Results: We determined the phylogenetic relationships of all known Francisella species, including some for which the phylogenetic positions were previously uncertain. The genus Francisella could be divided into two main genetic clades: one included F. tularensis, F. novicida, F. hispaniensis and Wolbachia persica, and another included F. philomiragia and F. noatunensis. Some Francisella species were found to have significant recombination frequencies. However, the fish pathogen F. noatunensis subsp. noatunensis was an exception due to it exhibiting a highly clonal population structure similar to the human pathogen F. tularensis. Conclusions: The genus Francisella can be divided into two main genetic clades occupying both terrestrial and marine habitats. However, our analyses suggest that the ancestral Francisella species originated in a marine habitat. The observed genome to genome variation in gene content and IS elements of different species supports the view that similar evolutionary paths of host adaptation developed independently in F. tularensis (infecting mammals) and F. noatunensis subsp. noatunensis (infecting fish).

Here we report the complete, accurate 1.89-Mb genome sequence of Francisella tularensis subsp. holarctica strain FSC200, isolated in 1998 in the Swedish municipality Ljusdal, which is in an area where tularemia is highly endemic. This genome is important because strain FSC200 has been extensively used for functional and genetic studies of Francisella and is well-characterized.

Background: Recent advances in sequencing technologies offer promising tools for generating large numbers of genomes, larger typing databases and improved mapping of environmental bacterial diversity. However, DNA-based methods for the detection of Francisella were developed with limited knowledge about genetic diversity. This, together with the high sequence identity between several Francisella species, means there is a high risk of false identification and detection of the highly virulent pathogen Francisella tularensis. Moreover, phylogenetic reconstructions using single or limited numbers of marker sequences often result in incorrect tree topologies and inferred evolutionary distances. The recent growth in publicly accessible whole-genome sequences now allows evaluation of published genetic markers to determine optimal combinations of markers that minimise both time and laboratory costs. Results: In the present study, we evaluated 38 previously published DNA markers and the corresponding PCR primers against 42 genomes representing the currently known diversity of the genus Francisella. The results highlight that PCR assays for Francisella tularensis are often complicated by low specificity, resulting in a high probability of false positives. A method to select a set of one to seven markers for obtaining optimal phylogenetic resolution or diagnostic accuracy is presented. Conclusions: Current multiple-locus sequence-typing systems and detection assays of Francisella, could be improved by redesigning some of the primers and reselecting typing markers. The use of only a few optimally selected sequence-typing markers allows construction of phylogenetic topologies with almost the same accuracy as topologies based on whole-genome sequences.

The F. tularensis type A strain FSC198 from Slovakia and a second strain FSC043, which has attenuated virulence, are both considered to be derivatives of the North American F. tularensis type A strain SCHU S4. These strains have been propagated under different conditions: the FSC198 has undergone natural propagation in the environment, while the strain FSC043 has been cultivated on artificial media in laboratories. Here, we have compared the genome sequences of FSC198, FSC043, and SCHU S4 to explore the possibility that the contrasting propagation conditions may have resulted in different mutational patterns. We found four insertion/deletion events (INDELs) in the strain FSC043, as compared to the SCHU S4, while no single nucleotide polymorphisms (SNPs) or variable number of tandem repeats (VNTRs) were identified. This result contrasts with previously reported findings for the strain FSC198, where eight SNPs and three VNTR differences, but no INDELs exist as compared to the SCHU S4 strain. The mutations detected in the laboratory and naturally propagated type A strains, respectively, demonstrate distinct patterns supporting that analysis of mutational spectra might be a useful tool to reveal differences in past growth conditions. Such information may be useful to identify leads in a microbial forensic investigation.

A robust, rapid and flexible real-time PCR assay for hierarchical genetic typing of clinical and environmental isolates of Francisella is presented. Typing markers were found by multiple genome and gene comparisons, from which 23 canonical single nucleotide polymorphisms (canSNPs) and 11 canonical insertion-deletion mutations (canINDELs) were selected to provide phylogenetic guidelines for classification from genus to isolate level. The specificity of the developed assay, which uses 68 wells of a 96-well real-time PCR format with a detection limit of 100 pg DNA, was assessed using 62 Francisella isolates of diverse genetic and geographical origins. It was then successfully used for typing 14 F. tularensis subsp. holarctica isolates obtained from tularemia patients in Sweden in 2008 and five more genetically diverse Francisella isolates of global origins. When applied to human ulcer specimens for direct pathogen detection the results were incomplete due to scarcity of DNA, but sufficient markers were identified to detect fine-resolution differences among F. tularensis subsp. holarctica isolates causing infection in the patients. In contrast to other real-time PCR assays for Francisella, which are typically designed for specific detection of a species, subspecies, or strain, this type of assay can be easily tailored to provide appropriate phylogenetic and/or geographical resolution to meet the objectives of the analysis.

This thesis is about analyzing genetic differences among isolates of Francisella tularensis – the tularemia-causing bacterium. To elucidate how these bacterial isolates are related, and their geographical and genetic origins, I have developed typing assays for Francisella and used them to study the epidemiology of tularemia.

Tularemia is an infectious disease of humans and other mammals found throughout the Northern Hemisphere. The severity of the disease depends on the type of F. tularensis causing the infection. In Sweden, as in other countries of Europe and Eurasia, tularemia is caused by F. tularensis subsp. holarctica, while other varieties of the bacterium occur in Middle Asia and North America. It is important to identify a tularemia infection promptly in order to initiate the correct antibiotic treatment. A rapid identification of the causative F. tularensis variety gives additional clinical information. In recent years, several genomes of various Francisella strains have been sequenced, and in this thesis, I have utilized these genomes to identify genetic markers.

In studies reported in the first paper (I) appended to the thesis, we identified and analyzed insertion/deletion mutations (INDELs) inferred to have resulted from a sequence repeat-mediated excision mechanism. We found eight new Regions of Difference (RDs) among Francisella strains. Using RDs together with single nucleotide polymorphisms (SNPs), we were able to predict an evolutionary scenario for F. tularensis in which Francisella novicida was the oldest variety while F. tularensis subsp. holarctica was the youngest. We also found that all virulence-attenuated isolates analyzed had deletions at two specific genetic regions - denoted RD18 and RD19 – suggesting that repeat-mediated excision is a mechanism of attenuation in F. tularensis.

In subsequent studies (presented in paper II), we developed a combined analysis of INDELs lacking flanking repeats and variable number of tandem repeats (VNTRs). Both markers could be assayed using the same analytical equipment. The inclusion of INDELs provided increased phylogenetic robustness compared with the use of VNTRs alone, while still maintaining a high level of genetic resolution.

In analyses described in the next paper (III), we selected INDELs from paper (II) and discovered novel SNPs by DNA comparisons of multiple Francisella strains. Thirty-four phylogenetically informative genetic markers were included in a hierarchical real-time PCR array for rapid and robust characterization of Francisella. We successfully used the assay to genotype 14 F. tularensis isolates from tularemia patients and DNA in six clinical ulcer specimens.

Finally, in paper (IV) we demonstrated a strategy to enhance epidemiological investigations of tularemia by combining GIS-mapping of disease-transmission place collected from patient interviews, with high-resolution genotyping of F. tularensis subsp. holarctica isolates recovered from tularemia patients. We found the geographic distributions of specific F. tularensis subsp. holarctica sub-populations to be highly localized during outbreaks (infections by some genotypes being restricted to areas as small as 2 km²), indicative of a landscape epidemiology of tularemia with distinct point sources of infection.

In conclusion, the results acquired during the studies underlying this thesis contribute to our understanding of the genetic genealogy of tularemia at both global and local outbreak scales.

Summer outbreaks of tularemia that occurred from 1995 through 2005 in 2 locations in Sweden affected 441 persons. We performed an epidemiologic investigation of these outbreaks using a novel strategy, involving high-resolution genotyping of Francisella tularensis isolates obtained from 136 patients (using 18 genetic markers developed from 6 F. tularensis genome sequences) and interviews with the patients. Strong spatial associations were found between F. tularensis subpopulations and the places of disease transmission; infection by some subpopulations occurred within areas as small as 2 km(2), indicating unidentified environmental point sources of tularemia. In both locations, disease clusters were associated with recreational areas beside water, and genetic subpopulations were present throughout the tularemia season and persisted over years. High-resolution genotyping in combination with patients' statements about geographic places of disease transmission provided valuable indications of likely sources of infection and the causal genotypes during these tularemia outbreaks.