Phenotypic and genotypic discrimination of Francisella tularensis ssp. holarctica clades

Francisella tularensis is the causative agent of tularemia, a zoonotic disease with a wide host range. F. tularensis ssp. holarctica (Fth) is of clinical relevance for European countries, including Germany. Whole genome sequencing methods, including canonical Single Nucleotide Polymorphism (canSNP) typing and whole genome SNP typing, have revealed that European Fth strains belong to a few monophyletic populations. The majority of German Fth isolates belong to two basal phylogenetic clades B.6 (biovar I) and B.12 (biovar II). Strains of B.6 and B.12 seem to differ in their pathogenicity, and it has been shown that strains of biovar II are resistant against erythromycin. In this study, we present data corroborating our previous data demonstrating that basal clade B.12 can be divided into clades B.71 and B.72. By applying phylogenetic whole genome analysis as well as proteome analysis, we could verify that strains of these two clades are distinct from one another. This was confirmed by measuring the intensity of backscatter light on bacteria grown in liquid media. Strains belonging to clades B.6, B.71 or B.72 showed clade-specific backscatter growth curves. Furthermore, we present the whole genome sequence of strain A-1341, as a reference genome of clade B.71, and whole proteomes comparison of Fth strains belonging to clades B.6, B.71 and B.72. Further research is necessary to investigate phenotypes and putative differences in pathogenicity of the investigated different clades of Fth to better understand the relationship between observed phenotypes, pathogenicity and distribution of Fth strains.

Genomic characterization of Francisella tularensis and other diverse Francisella species from complex samples

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.

Proteomic Signatures of Antimicrobial Resistance in Yersinia pestis and Francisella tularensis

Antimicrobial resistance (AMR) is a well-recognized, widespread, and growing issue of concern. With increasing incidence of AMR, the ability to respond quickly to infection with or exposure to an AMR pathogen is critical. Approaches that could accurately and more quickly identify whether a pathogen is AMR also are needed to more rapidly respond to existing and emerging biological threats. We examined proteins associated with paired AMR and antimicrobial susceptible (AMS) strains of Yersinia pestis and Francisella tularensis, causative agents of the diseases plague and tularemia, respectively, to identify whether potential existed to use proteins as signatures of AMR. We found that protein expression was significantly impacted by AMR status. Antimicrobial resistance-conferring proteins were expressed even in the absence of antibiotics in growth media, and the abundance of 10–20% of cellular proteins beyond those that directly confer AMR also were significantly changed in both Y. pestis and F. tularensis. Most strikingly, the abundance of proteins involved in specific metabolic pathways and biological functions was altered in all AMR strains examined, independent of species, resistance mechanism, and affected cellular antimicrobial target. We have identified features that distinguish between AMR and AMS strains, including a subset of features shared across species with different resistance mechanisms, which suggest shared biological signatures of resistance. These features could form the basis of novel approaches to identify AMR phenotypes in unknown strains.

FlexTaxD: flexible modification of taxonomy databases for improved sequence classification

SUMMARY

The Flexible Taxonomy Database (FlexTaxD) framework provides a method for modification and merging official and custom taxonomic databases to create improved databases. Using such databases will increase accuracy and precision of existing methods to classify sequence reads.

AVAILABILITY AND IMPLEMENTATION

Source code is freely available at https://github.com/FOI-Bioinformatics/flextaxd and installable through Bioconda.

Reorganized Genomic Taxonomy of Francisellaceae Enables Design of Robust Environmental PCR Assays for Detection of Francisella tularensis

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.

Biological amplification of low frequency mutations unravels laboratory culture history of the bio-threat agent Francisella tularensis

Challenges of investigating a suspected bio attack include establishing if microorganisms have been cultured to produce attack material and to identify their source. Addressing both issues, we have investigated genetic variations that emerge during laboratory culturing of the bacterial pathogen Francisella tularensis. Key aims were to identify genetic variations that are characteristic of laboratory culturing and explore the possibility of using biological amplification to identify genetic variation present at exceedingly low frequencies in a source sample. We used parallel serial passage experiments and high-throughput sequencing of F. tularensis to explore the genetic variation. We found that during early laboratory culture passages of F. tularensis, gene duplications emerged in the pathogen genome followed by single-nucleotide polymorphisms in genes for bacterial capsule synthesis. Based on a biological enrichment scheme and the use of high-throughput sequencing, we identified genetic variation that likely pre-existed in a source sample. The results support that capsule synthesis gene mutations are common during laboratory culture, and that a biological amplification strategy is useful for linking a F. tularensis sample to a specific laboratory variant among many highly similar variants.

A likelihood ratio-based approach for improved source attribution in microbiological forensic investigations

A common objective in microbial forensic investigations is to identify the origin of a recovered pathogenic bacterium by DNA sequencing. However, there is currently no consensus about how degrees of belief in such origin hypotheses should be quantified, interpreted, and communicated to wider audiences. To fill this gap, we have developed a concept based on calculating probabilistic evidential values for microbial forensic hypotheses. The likelihood-ratio method underpinning this concept is widely used in other forensic fields, such as human DNA matching, where results are readily interpretable and have been successfully communicated in juridical hearings. The concept was applied to two case scenarios of interest in microbial forensics: (1) identifying source cultures among series of very similar cultures generated by parallel serial passage of the Tier 1 pathogen Francisella tularensis, and (2) finding the production facilities of strains isolated in a real disease outbreak caused by the human pathogen Listeria monocytogenes. Evidence values for the studied hypotheses were computed based on signatures derived from whole genome sequencing data, including deep-sequenced low-frequency variants and structural variants such as duplications and deletions acquired during serial passages. In the F. tularensis case study, we were able to correctly assign fictive evidence samples to the correct culture batches of origin on the basis of structural variant data. By setting up relevant hypotheses and using data on cultivated batch sources to define the reference populations under each hypothesis, evidential values could be calculated. The results show that extremely similar strains can be separated on the basis of amplified mutational patterns identified by high-throughput sequencing. In the L. monocytogenes scenario, analyses of whole genome sequence data conclusively assigned the clinical samples to specific sources of origin, and conclusions were formulated to facilitate communication of the findings. Taken together, these findings demonstrate the potential of using bacterial whole genome sequencing data, including data on both low frequency SNP signatures and structural variants, to calculate evidence values that facilitate interpretation and communication of the results. The concept could be applied in diverse scenarios, including both epidemiological and forensic source tracking of bacterial infectious disease outbreaks.

Antibiotic susceptibility in vitro of Francisella tularensis subsp. holarctica isolates from Germany

Background

Tularaemia is a zoonotic disease caused by the bacterium Francisella tularensis. In Germany, the disease is still rare (e.g. 34 human cases reported in 2015). There is a lack of data about the susceptibility of F. tularensis strains to antibiotics, because many cases are diagnosed using serological assays only.

Objectives

The antibiotic susceptibility in vitro of F. tularensis subsp. holarctica strains isolated in Germany was assessed to determine whether the currently recommended empirical therapy is still adequate.

Methods

A total of 128 F. tularensis strains were investigated that were collected between 2005 and 2014 in Germany from wild animals, ticks and humans. All isolates were genotyped using real-time PCR assays targeting canonical SNPs, and antibiotic susceptibility was tested using MIC test strips on agar plates. MIC values were interpreted using CLSI breakpoints.

Results

The strains were susceptible to antibiotics commonly recommended for tularaemia therapy, i.e. aminoglycosides (MIC90 values: gentamicin 1 mg/L; streptomycin 4.0 mg/L), tetracyclines (MIC90 values: tetracycline 0.5 mg/L; doxycycline 1.5 mg/L) and quinolones (MIC90 value: ciprofloxacin 0.064 mg/L). Chloramphenicol (MIC90 value: 3.0 mg/L) may be of value in treatment of tularaemia meningitis. Ninety-four isolates were susceptible to erythromycin, which defines biovar I (genotypes B.4 and B.6); 34 were resistant (biovar II; genotype B.12).

Conclusions

The F. tularensis isolates investigated in this study showed the typical antibiotic susceptibility pattern that was previously observed in other countries. Therefore, recommendations for empirical antibiotic therapy of tularaemia can remain unchanged. However, antibiotic susceptibility testing of clinical isolates should be performed whenever possible.

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.

Introduction and persistence of tularemia in Bulgaria

Introduction

Outbreaks of the zoonotic disease tularemia occurred in north-east Bulgaria in the 1960s. Then came 30 years of epidemiological silence until new outbreaks occurred in west Bulgaria in the 1990s. To investigate how bacterial strains of Francisella tularensis causing tularemia in wildlife and humans in the 1960s and the 1990s were related, we explored their genetic diversity.

Material and methods

Ten F. tularensis genomes from the 1960s (n=3) and the 1990s (n=7) were sequenced, assigned to canonical single-nucleotide polymorphism (canSNP) clades, and compared to reference genomes. We developed four new canSNP polymerase chain reaction (PCR) assays based on the genome sequence information.

Results and discussion

The genetic analysis showed that the outbreaks in the 1960s as well as in the 1990s involved multiple clones and new genetic diversity. The smallest genetic difference found between any of the Bulgarian strains was five SNPs between the strains L2 and 81 isolated 43 years apart, indicating that F. tularensis may persist locally over long time periods without causing outbreaks. The existence of genetically highly similar strain-pairs isolated the same year in the same area from different hosts supports a hypothesis of local expansion of clones during outbreaks. Close relationship (two SNPs) was found between one strain isolated 1961 in northeast Bulgaria and one strain isolated 5 years before in USSR. Historical data coinciding with the actual time point describe the introduction of water rats from USSR into the Bulgarian outbreak area, which may explain the close genetic relationship and the origin of the outbreak.

Conclusion

Genome analysis of strains from two outbreaks in the 1960s and the 1990s provided valuable information on the genetic diversity and persistence of F. tularensis in Bulgaria.

Hare-to-human transmission of Francisella tularensis subsp. holarctica, Germany

In November 2012, a group of 7 persons who participated in a hare hunt in North Rhine-Westphalia, Germany, acquired tularemia. Two F. tularensis subsp. holarctica isolates were cultivated from human and hare biopsy material. Both isolates belonged to the FTN002–00 genetic subclade (derived for single nucleotide polymorphisms B.10 and B.18), thus indicating likely hare-to-human transmission.

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.

Phylogeography of Francisella tularensis subspecies holarctica in Finland, 1993-2011

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.

Genome sequence of Coxiella burnetii strain Namibia

We present the whole genome sequence and annotation of the Coxiella burnetii strain Namibia. This strain was isolated from an aborting goat in 1991 in Windhoek, Namibia. The plasmid type QpRS was confirmed in our work. Further genomic typing placed the strain into a unique genomic group. The genome sequence is 2,101,438 bp long and contains 1,979 protein-coding and 51 RNA genes, including one rRNA operon. To overcome the poor yield from cell culture systems, an additional DNA enrichment with whole genome amplification (WGA) methods was applied. We describe a bioinformatics pipeline for improved genome assembly including several filters with a special focus on WGA characteristics.

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.

CanSNPer: a hierarchical genotype classifier of clonal pathogens

SUMMARY

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.

AVAILABILITY AND IMPLEMENTATION

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

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.

An improved multiple-locus variable-number of tandem repeat analysis (MLVA) for the fish pathogen Francisella noatunensis using capillary electrophoresis

BACKGROUND

Francisellosis, caused by the bacterium Francisella noatunensis subsp. noatunensis, remains a serious threat to Atlantic cod (Gadhus morhua) farming in Norway and potentially in other countries. As outbreak strains appear clonal in population structure, access to highly discriminatory typing tools is critical for understanding the epidemiology of francisellosis infections in aquaculture. In this study, a simplified multiple-locus variable-number of tandem repeat analysis (MLVA) targeting five highly polymorphic variable number of tandem repeat (VNTR) loci in a single multiplex PCR was developed to rapidly discriminate between outbreak strains.

RESULTS

The assay resulted in identification of at least 13 different allelic profiles or subpopulations among 91 F. noatunensis isolates from farmed cod in Norway. The VNTR loci appear relatively stable, with isolates originating from individual outbreaks showing identical MLVA profiles following repeated passage. MLVA displayed greater discriminatory power than pulse-field gel electrophoresis (PFGE). Both MLVA and PFGE show good epidemiological concordance by their abilities to separate outbreak strains from epidemiologically unrelated isolates.

CONCLUSIONS

The MLVA method presented here is robust, easy to perform and provides a good alternative to other typing systems for F. noatunensis subsp. noatunensis and epidemiological study of francisellosis in cod.

German Francisella tularensis isolates from European brown hares (Lepus europaeus) reveal genetic and phenotypic diversity

BACKGROUND

Tularemia 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.

RESULTS

Cultivation 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.

CONCLUSIONS

F. 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.