Caroline Öhrman

Genome characterisation of the genus Francisella reveals insight into similar evolutionary paths in pathogens of mammals and fish

BackgroundPrior 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.ResultsWe 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.ConclusionsThe 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).

An Outbreak of Respiratory Tularemia Caused by Diverse Clones of Francisella tularensis

BACKGROUND The bacterium Francisella tularensis is recognized for its virulence, infectivity, genetic homogeneity, and potential as a bioterrorism agent. Outbreaks of respiratory tularemia, caused by inhalation of this bacterium, are poorly understood. Such outbreaks are exceedingly rare, and F. tularensis is seldom recovered from clinical specimens. METHODS A localized outbreak of tularemia in Sweden was investigated. Sixty-seven humans contracted laboratory-verified respiratory tularemia. F. tularensis subspecies holarctica was isolated from the blood or pleural fluid of 10 individuals from July to September 2010. Using whole-genome sequencing and analysis of single-nucleotide polymorphisms (SNPs), outbreak isolates were compared with 110 archived global isolates. RESULTS There were 757 SNPs among the genomes of the 10 outbreak isolates and the 25 most closely related archival isolates (all from Sweden/Finland). Whole genomes of outbreak isolates were >99.9% similar at the nucleotide level and clustered into 3 distinct genetic clades. Unexpectedly, high-sequence similarity grouped some outbreak and archival isolates that originated from patients from different geographic regions and up to 10 years apart. Outbreak and archival genomes frequently differed by only 1-3 of 1 585 229 examined nucleotides. CONCLUSIONS The outbreak was caused by diverse clones of F. tularensis that occurred concomitantly, were widespread, and apparently persisted in the environment. Multiple independent acquisitions of F. tularensis from the environment over a short time period suggest that natural outbreaks of respiratory tularemia are triggered by environmental cues. The findings additionally caution against interpreting genome sequence identity for this pathogen as proof of a direct epidemiological link.

Eight New Genomes and Synthetic Controls Increase the Accessibility of Rapid Melt-MAMA SNP Typing of Coxiella burnetii

The case rate of Q fever in Europe has increased dramatically in recent years, mainly because of an epidemic in the Netherlands in 2009. Consequently, there is a need for more extensive genetic characterization of the disease agent Coxiella burnetii in order to better understand the epidemiology and spread of this disease. Genome reference data are essential for this purpose, but only thirteen genome sequences are currently available. Current methods for typing C. burnetii are criticized for having problems in comparing results across laboratories, require the use of genomic control DNA, and/or rely on markers in highly variable regions. We developed in this work a method for single nucleotide polymorphism (SNP) typing of C. burnetii isolates and tissue samples based on new assays targeting ten phylogenetically stable synonymous canonical SNPs (canSNPs). These canSNPs represent previously known phylogenetic branches and were here identified from sequence comparisons of twenty-one C. burnetii genomes, eight of which were sequenced in this work. Importantly, synthetic control templates were developed, to make the method useful to laboratories lacking genomic control DNA. An analysis of twenty-one C. burnetii genomes confirmed that the species exhibits high sequence identity. Most of its SNPs (7,493/7,559 shared by >1 genome) follow a clonal inheritance pattern and are therefore stable phylogenetic typing markers. The assays were validated using twenty-six genetically diverse C. burnetii isolates and three tissue samples from small ruminants infected during the epidemic in the Netherlands. Each sample was assigned to a clade. Synthetic controls (vector and PCR amplified) gave identical results compared to the corresponding genomic controls and are viable alternatives to genomic DNA. The results from the described method indicate that it could be useful for cheap and rapid disease source tracking at non-specialized laboratories, which requires accurate genotyping, assay accessibility and inter-laboratory comparisons.

Genomic analyses of Francisella tularensis strains confirm disease transmission from drinking water sources, Turkey, 2008, 2009 and 2012.

Waterborne epidemics of tularaemia caused by Francisella tularensis are increasingly reported in Turkey. We have used whole genome sequencing to investigate if F. tularensis isolated from patients could be traced back to drinking water sources. Tonsil swabs from 33 patients diagnosed with oropharyngeal tularaemia in three outbreaks and 140 water specimens were analysed. F. tularensis subsp. holarctica was confirmed by microagglutination and PCR in 12 patients and five water specimens. Genomic analysis of three pairs of patient and water isolates from outbreaks in Sivas, Çorum, and Kocaeli showed the isolates to belong to two new clusters of the F. tularensis B.12 genetic clade. The clusters were defined by 19 and 15 single nucleotide polymorphisms (SNPs) in a multiple alignment based on 507 F. tularensis genomes. One synonymous SNP was chosen as a new canonical SNP (canSNP) for each cluster for future use in diagnostic assays. No SNP was identified between the genomes from the patient–water pair of isolates from Kocaeli, one SNP between the pair of isolates from Sivas, whereas the pair from Çorum differed at seven SNPs. These results illustrate the power of whole genome sequencing for tracing F. tularensis patient isolates back to their environmental source.

Increased knowledge of Francisella genus diversity highlights the benefits of optimised DNA-based assays

BackgroundRecent 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.ResultsIn 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.ConclusionsCurrent 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.

Complete Genome Sequence of Francisella endociliophora Strain FSC1006, Isolated from a Laboratory Culture of the Marine Ciliate Euplotes raikovi

ABSTRACT A strain of Francisella endociliophora was isolated from a laboratory culture of the marine ciliate Euplotes raikovi. Here, we report the complete genome sequence of the bacterial strain FSC1006 (Francisella Strain Collection, Swedish Defence Research Agency, Umeå, Sweden).

Francisella tularensis subsp. tularensis Group A.I, United States

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.

Clonality of erythromycin resistance in Francisella tularensis

OBJECTIVES We analysed diverse strains of Francisella tularensis subsp. holarctica to assess if its division into biovars I and II is associated with specific mutations previously linked to erythromycin resistance and to determine the distribution of this resistance trait across this subspecies. METHODS Three-hundred and fourteen F. tularensis subsp. holarctica strains were tested for erythromycin susceptibility and whole-genome sequences for these strains were examined for SNPs in genes previously associated with erythromycin resistance. Each strain was assigned to a global phylogenetic framework using genome-wide canonical SNPs. The contribution of a specific SNP to erythromycin resistance was examined using allelic exchange. The geographical distribution of erythromycin-resistant F. tularensis strains was further investigated by literature search. RESULTS There was a perfect correlation between biovar II strains (erythromycin resistance) and the phylogenetic group B.12. Only B.12 strains had an A → C SNP at position 2059 in the three copies of the rrl gene. Introducing 2059C into an rrl gene of an erythromycin-susceptible F. tularensis strain resulted in resistance. An additional 1144 erythromycin-resistant strains were identified from the scientific literature, all of them from Eurasia. CONCLUSIONS Erythromycin resistance in F. tularensis is caused by an A2059C rrl gene mutation, which exhibits a strictly clonal inheritance pattern found only in phylogenetic group B.12. This group is an extremely successful clone, representing the most common type of F. tularensis throughout Eurasia.

Long-range dispersal moved Francisella tularensis into Western Europe from the East

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.

Complete Genome Sequence of Francisella guangzhouensis Strain 08HL01032T, Isolated from Air-Conditioning Systems in China.

ABSTRACT We present the complete genome sequence of Francisella guangzhouensis strain 08HL01032T, which consists of one chromosome (1,658,482 bp) and one plasmid (3,045 bp) with G+C contents of 32.0% and 28.7%, respectively.