Search

Francisella tularensis, the bacterium that causes the zoonosis tularemia, and its genetic near neighbor species, can be difficult or impossible to cultivate from complex samples. Thus, there is a lack of genomic information for these species that has, among other things, limited the development of robust detection assays for F. tularensis that are both specific and sensitive. The objective of this study was to develop and validate approaches to capture, enrich, sequence, and analyze Francisella DNA present in DNA extracts generated from complex samples. RNA capture probes were designed based upon the known pan genome of F. tularensis and other diverse species in the family Francisellaceae. Probes that targeted genomic regions also present in non-Francisellaceae species were excluded, and probes specific to particular Francisella species or phylogenetic clades were identified. The capture-enrichment system was then applied to diverse, complex DNA extracts containing low-level Francisella DNA, including human clinical tularemia samples, environmental samples (i.e., animal tissue and air filters), and whole ticks/tick cell lines, which was followed by sequencing of the enriched samples. Analysis of the resulting data facilitated rigorous and unambiguous confirmation of the detection of F. tularensis or other Francisella species in complex samples, identification of mixtures of different Francisella species in the same sample, analysis of gene content (e.g., known virulence and antimicrobial resistance loci), and high-resolution whole genome-based genotyping. The benefits of this capture-enrichment system include: even very low target DNA can be amplified; it is culture-independent, reducing exposure for research and/or clinical personnel and allowing genomic information to be obtained from samples that do not yield isolates; and the resulting comprehensive data not only provide robust means to confirm the presence of a target species in a sample, but also can provide data useful for source attribution, which is important from a genomic epidemiology perspective., Competing Interests: The authors have declared that no competing interests exist.

Tularaemia, caused by the bacterium Francisella tularensis, is endemic in Sweden and is poorly understood. The aim of this study was to evaluate the effectiveness of three different genetic typing systems to link a genetic type to the source and place of tularemia infection in Sweden. Canonical single nucleotide polymorphisms (canSNPs), MLVA including five variable number of tandem repeat loci and PmeI-PFGE were tested on 127 F. tularensis positive specimens collected from Swedish case-patients. All three typing methods identified two major genetic groups with near-perfect agreement. Higher genetic resolution was obtained with canSNP and MLVA compared to PFGE; F. tularensis samples were first assigned into ten phylogroups based on canSNPs followed by 33 unique MLVA types. Phylogroups were geographically analysed to reveal complex phylogeographic patterns in Sweden. The extensive phylogenetic diversity found within individual counties posed a challenge to linking specific genetic types with specific geographic locations. Despite this, a single phylogroup (B.22), defined by a SNP marker specific to a lone Swedish sequenced strain, did link genetic type with a likely geographic place. This result suggests that SNP markers, highly specific to a particular reference genome, may be found most frequently among samples recovered from the same location where the reference genome originated. This insight compels us to consider whole-genome sequencing (WGS) as the appropriate tool for effectively linking specific genetic type to geography. Comparing the WGS of an unknown sample to WGS databases of archived Swedish strains maximizes the likelihood of revealing those rare geographically informative SNPs.

Francisella tularensis, the etiologic agent of tularemia and a Class A Select Agent, is divided into three subspecies and multiple subpopulations that differ in virulence and geographic distribution. Given these differences, there is a need to rapidly and accurately determine if a strain is F. tularensis and, if it is, assign it to subspecies and subpopulation. We designed TaqMan real-time PCR genotyping assays using eleven single nucleotide polymorphisms (SNPs) that were potentially specific to closely related groups within the genus Francisella, including numerous subpopulations within F. tularensis species. We performed extensive validation studies to test the specificity of these SNPs to particular populations by screening the assays across a set of 565 genetically and geographically diverse F. tularensis isolates and an additional 21 genetic near-neighbor (outgroup) isolates. All eleven assays correctly determined the genetic groups of all 565 F. tularensis isolates. One assay differentiates F. tularensis, F. novicida, and F. hispaniensis from the more genetically distant F. philomiragia and Francisella-like endosymbionts. Another assay differentiates F. tularensis isolates from near neighbors. The remaining nine assays classify F. tularensis-confirmed isolates into F. tularensis subspecies and subpopulations. The genotyping accuracy of these nine assays diminished when tested on outgroup isolates (i.e. non F. tularensis), therefore a hierarchical approach of assay usage is recommended wherein the F. tularensis-specific assay is used before the nine downstream assays. Among F. tularensis isolates, all eleven assays were highly sensitive, consistently amplifying very low concentrations of DNA. Altogether, these eleven TaqMan real-time PCR assays represent a highly accurate, rapid, and sensitive means of identifying the species, subspecies, and subpopulation of any F. tularensis isolate if used in a step-wise hierarchical scheme. These assays would be very useful in clinical, epidemiological, and/or forensic investigations involving F. tularensis.