Next Article in Journal
Hallmarks of Basidiomycete Soft- and White-Rot in Wood-Decay -Omics Data of Two Armillaria Species
Previous Article in Journal
HbxB Is a Key Regulator for Stress Response and β-Glucan Biogenesis in Aspergillus nidulans
Previous Article in Special Issue
Genotyping of Francisella tularensis subsp. holarctica from Hares in Germany
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 1
10.3390/microorganisms9010146
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Caroline Öhrman
1,
Jason W. Sahl
2,
Andreas Sjödin
1,
Ingrid Uneklint
1,
Rebecca Ballard
2,
Linda Karlsson
1,
Ryelan F. McDonough
2,
David Sundell
1,
Kathleen Soria
2,
Stina Bäckman
1,
Kitty Chase
3,
Björn Brindefalk
1,
Shanmuga Sozhamannan
4,
Adriana Vallesi
5,
Emil Hägglund
1,
Jose Gustavo Ramirez-Paredes
6,
Johanna Thelaus
1,
Duncan Colquhoun
7,
Kerstin Myrtennäs
1,
Dawn Birdsell
2,
Anders Johansson
8,
David M. Wagner
2 and
Mats Forsman
1,*add Show full author list remove Hide full author list
1
CBRN Defence and Security, Swedish Defence Research Agency, FOI, SE 901 82 Umeå, Sweden
2
Pathogen and Microbiome Institute, Northern Arizona University, Flagstaff, AZ 86011, USA
3
US Army Medical Research Institute, Fort Detrick, MD 21702, USA
4
Logistics Management Institute supporting Defense Biological Product Assurance Office (DBPAO) Joint Project Lead, CBRND Enabling Biotechnologies (JPL CBRND EB), Frederick, MD 21702, USA
5
School of Biosciences and Veterinary Medicine, University of Camerino, 62032 Camerino, Italy
6
Ridgeway Biologicals Limited a Ceva Santé Animale Company, Units 1-3 Old Station Business Park, Compton, Berkshire, England RG20 6NE, UK
7
Fish Health Research Group, Norwegian Veterinary Institute, Oslo, Pb 750 Sentrum, 23 N-0106 Oslo, Norway
8
Department of Clinical Microbiology, Clinical Bacteriology, and Laboratory for Molecular Infection Medicine Sweden (MIMS), Umeå University, SE-901 85 Umeå, Sweden
*
Author to whom correspondence should be addressed.
Microorganisms 2021, 9(1), 146; https://doi.org/10.3390/microorganisms9010146
Submission received: 21 December 2020 / Revised: 1 January 2021 / Accepted: 4 January 2021 / Published: 11 January 2021
(This article belongs to the Special Issue Epidemiology of Tularemia and Francisella tularensis)
Browse Figures
Versions Notes

Abstract

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

1. Introduction

The family Francisellaceae consists of a group of Gram-negative, non-motile, aerobic bacteria with diverse life cycles. In recent years, this family has been expanded vastly and new genera as well as many new species have been included and characterized. The family is currently comprised of three genera: Francisella, Allofrancisella [1], and Pseudofrancisella [2]. The genus Pseudofrancisella contains a single validly published species isolated from estuarine seawater [2]. The genus Allofrancisella contains three validly published species, of which all have been isolated from water samples collected from cooling towers of air conditioning systems [1]. Despite this diversity within the Francisellaceae, and nine different Francisella species (Table A1), only the three F. tularensis subspecies, i.e., subspecies tularensis (type A), subspecies holarctica (type B), and subspecies mediasiatica [3], are of significant clinical relevance to humans. The subspecies holarctica and tularensis cause tularemia in humans [4]. Knowledge relating to the virulence of F. tularensis subspecies mediasiatica in humans, which is primarily found in Central Asia, is limited; however, an intermediate position between the other two subspecies in terms of virulence is assumed [5,6]. F. tularensis subspecies tularensis is highly infectious in humans and a number of other mammalian species, and is most feared for its potential as a biological weapon [7] or bioterrorism agent [8].
F. noatunensis, F. orientalis, F. marina [9] and F. halioticida [10] are the etiological agents of francisellosis [11,12,13] in cold water fish species, warm water fish species, and mollusks, respectively [10,14]. The genus Francisella also includes opportunistic pathogens that can cause disease in immunocompromised humans (i.e., F. novicida, F. hispaniensis, and F. philomiragia) [4,15,16,17,18,19]. In addition, the genus Francisella includes Francisella-like endosymbionts (FLE) of ticks and ciliates (i.e., F. persica, F. endociliophora, and F. adeliensis) [20,21,22]. Related ‘environmental’ taxa, including F. salina, F. salimarina, F. uliginis, F. frigiditurris, Allofrancisella spp., and Pseudofrancisella aestuarii [1,2,23,24,25,26,27,28,29,30], have been isolated or sequenced directly from samples of brackish water, sea water, water from air conditioning systems, cooling towers, and/or other environmental sources, including natural warm springs or ice machines, suggesting that Francisellaceae inhabit a wide diversity of ecological niches. This group includes both opportunistic pathogens and non-pathogenic species that may become aerosolized from air-cooling systems and also be captured on air filters [1,31,32]. Of note, the BioWatch program in the USA, which tests for the presence of pathogenic organisms on air filters, reported 149 false positive results between 2003 and 2014, of which the majority were associated with molecular assays for F. tularensis [32]. Given the ubiquity of diverse Francisella species in the environment, and the fact that many of them can apparently become aerosolized, it seems likely that the presence of these other species on air filters are the source of these false positive BioWatch results.
There are multiple natural sources of transmission of F. tularensis to humans including contaminated water and food; direct contact with infected animals such as rabbits, hares, squirrels, voles and other rodents; and bites from ticks, flies and mosquitoes or inhalation of aerosol [33]. However, the relative frequency of the different modes of transmission of F. tularensis to humans vary according to the geographic areas considered and local ecology. F. tularensis is well known for its environmental persistence and irregularity of outbreaks, sometimes separated by several years, that can occur in closely confined geographical regions (aka foci). The principal natural reservoir is largely unknown, although studies have shown that the organism may persist for more than 16 months in water or mud [34]. Hence, from a surveillance perspective, vectors and complex air, water, and sediment samples are all relevant. Analysis of vectors and complex environmental samples is challenging in relation to distinguishing the human disease-causing F. tularensis from genetically similar species that may also be present in the samples.
Many Francisella species are difficult or impossible to culture directly from the environment, which has limited our understanding of this genus. This is even true for F. tularensis, where cultivation of the agent from the environment has traditionally been performed by inoculation of environmental samples into laboratory animals, because direct isolation by culturing of environmental samples is commonly compromised by overgrowth of other bacterial species on most growth media. In recent years, F. tularensis (type B) was isolated directly from drinking water in Turkey by filtering samples of spring water through cellulose acetate membranes before cultivation on selective media [35]; however, this seems to be a rare exception. The FLEs, especially, are notoriously difficult to cultivate directly from their hosts, although there are some examples of successful isolations [22]. In addition, advances in culture-independent approaches and direct DNA sequencing of complex environmental samples have revealed an extensive diversity of Francisella near-neighbor species [36,37,38]. Thus, many more environmental Francisella species likely exist that remain to be identified and explored.
F. tularensis is a very monomorphic ISspecies with a very small accessory genome [39,40], possibly because its evolutionary history entailed a transition from a free-living bacterium to a host-associated pathogen [41]. The F. tularensis genome consists of a single circular chromosome of about 1.9 million base pairs with a very high AT content (~68%) [42]. Differences between F. tularensis and other Francisellaceae species that are primarily isolated from environmental sources include metabolic competence, which is higher among the environmental species, and signs of ongoing genome erosion in F. tularensis that has resulted in the deletion of metabolic pathway components (i.e., genes involved in amino acid biosynthesis) [41,42]. F. tularensis is also characterized by a marked increase in insertion sequences [39]; large numbers of insertion sequences are often associated with an evolutionary process resulting in the loss of gene functionality. Interestingly, the etiological agent of francisellosis in cold water fish, F. noatunensis, also contains large numbers of insertion sequences that are similar to those reported for F. tularensis, and the high degree of clonality that characterizes F. tularensis [41] is also found in F. orientalis [39]. Despite infecting different species (fish vs. mammals), it appears that these pathogenic lineages share certain important evolutionary features. The chromosome of F. tularensis has a duplicated Francisella pathogenicity island (FPI). In contrast, the genome of opportunistic F. novicida contains a single FPI. F. tularensis harbors a degenerated CRISPR/Cas system [43], whereas the more metabolically versatile F. novicida possesses a well-functioning Type II CRISPR/Cas system [44]. The primary role for CRISPR/Cas is in defense against invading foreign nucleic acids, a function that thus seems to be lost in F. tularensis. This might suggest a life cycle with little competition and/or association with other microorganisms. In support of this notion, there is an absence of recombination in F. tularensis strains, especially compared with F. novicida, which has a significant rate of recombination [39].
To further elucidate the genomic differences between the non-disease-causing relatives, or near-neighbors, within the family Francisellaceae, and the human disease-causing F. tularensis, we have analyzed a dataset comprising 499 whole-genome sequences, including 172 environmental F. tularensis near-neighbors and 327 F. tularensis genome sequences. We identified the very few genomic regions unique to F. tularensis and demonstrated the utility of specific PCR assays based on these unique regions for the robust detection and identification of F. tularensis.

2. Materials and Methods

2.1. Library Preparation

DNA libraries were generated at Northern Arizona University (NAU), the Swedish Defense Research agency (FOI), and SciLifeLabs in Uppsala, Sweden. At NAU, DNA Libraries were prepared using the KAPA Library Preparation Kits with SRPI Solution and Standard PCR Library Amplification/Illumina series (KAPA Biosystems, Code KK8232) using the following modifications: adapters and 8 bp index oligos from IDT® (Integrated DNA Technologies, San Diego, CA, USA) based on Kozarewa and Turner, 2011 [45], were used in place of those supplied in the KAPA preparation kit. The DNA sample (~1 µg) was fragmented using a SonicMan (Matrical) and fragment targets of 600–650 bp were selected using Agencourt AMPure XP beads (Beckman Coulter, Code A63882). At FOI, DNA libraries were prepared from 1 ng DNA following the Nextera® XT protocol (Illumina, San Diego, CA, USA). At SciLifeLabs, 35 DNA libraries were prepared from 1 µg DNA using TruSeq DNA sample prep kit v2 (Illumina, cat# FC-121-2001/2002) following the manufacturer’s instructions.

2.2. Genome Sequencing

In total, 217 isolates were sequenced with Illumina sequencing technologies using MiSeq, HiSeq2000, or HiSeq2500 instruments. MiSeq sequencing was performed using the 500-cycle MiSeq Reagent Kit V2 or the 600-cycle MiSeq Reagent Kit V3 at FOI or NAU following the reagent protocol. DNA libraries prepared with TruSeq kits were run on HiSeq sequencing platform located at SciLife. Sequencing was performed with paired end 200-cycles with V3 Reagents. For each sample, information about the number of cycles and instrument platform is available at SRA (see accession numbers in Table S1).

2.3. Genome Assembly

Samples sequenced at FOI were assembled using abyss-pe 2.2.2 [46] using an k-value of 51. Contigs shorter than 500 bp were removed and the assembly was polished using pilon v1.22 [47]. The assemblies were checked for contamination with kraken2 v2.07 using a kraken-standard database (all available genomes from Archaea, Bacteria, human, UniVec_Core and Viral). Samples sequenced at NAU were adapter-trimmed with bbduk v38.86 (https://sourceforge.net/projects/bbmap/) and assembled with SPAdes v3.13.1 [48].

2.4. Genome Sequences

The complete dataset used for analysis consisted of 499 genome sequences; 217 newly sequenced and 282 downloaded from NCBI (Table S1). Six additional genomes (GCF_003428165.1, GCF_003574425.1, GCF_002803295.2, GCF_003428155.1, GCF_003574485.1, GCA_006227905.1), reclassified as part of the Francisellaceae family by Genome Taxonomy Database (GTDB Release 05-RS95, 17th July 2020) [49], and Piscirickettsia salmonis LF 89 (GCF_000297215.2), which is genetically similar to Francisellaceae [39], were also included in the initial phylogenetic tree. Recommended representative genome assemblies are listed in Table S4.

2.5. Phylogenetic Trees

Phylogenies were made by first aligning each genome to the reference strain SCHU S4 (GCA_000008985.1) using progressive mauve (snapshot_2015_02_13) [50] with default settings, converting them into the FASTA format according to the reference coordinate system, and merging them into one multi FASTA. All SNPs were then extracted using SNP-sites (v2.51) keeping only ACTG sites [51]. Ten starting trees were calculated using RAxML-HPC-PTHREADS-SSE3 [52] with model GTRGAMMA (--no-bfgs) and 100 bootstraps were calculated for the best ML tree using the same model. From the bootstrap replicate tree, bipartitions on the best maximum likelihood tree were drawn using model GTRCAT. The tree was visualized in FigTree 1.4.3 and rooted on Fangia hongkongensis for the complete dataset (Figure 1), Francisella novicida strain U112 (CP000439.1) for the clade 1 subset (Figure 3), and Piscirickettsia salmonis LF 89 GCF_000297215.2 for Francisellaceace (Figure A1).

2.6. Average Nucleotide Identity

Average nucleotide identity (ANI) was calculated for the dataset consisting of 499 genomes using pyani v0.2.10 [53] and the ANIb (blast) method. The correlation table was aggregated per species by mean first on the x-axis and then on the y-axis using the aggregate function in the R package stats. The resulting table was visualized using function heatmap.2 from the R package ggplot2 (v3.1).

2.7. Genome Taxonomy of Francisellaceae

Information downloaded from Genome Taxonomy Database (GTDB) was incorporated into the phylogenetic tree together with ANI values. The taxonomy from GTDB was compared to the original NCBI taxonomy assignment and type strains not yet included in GTDB were used to update names of unassigned species in the GTDB backbone.

2.8. Identification of Clade-Specific Regions

For each genome assembly, coding region sequences (CDSs) were predicted with Prodigal v2.63 [54]; intergenic regions > 50 nucleotides were also extracted from genome assemblies. All regions were concatenated and clustered with cd-hit v4.8.1 [55] at an identification threshold of 90%. Each representative sequence from each CDS was aligned against all queried genomes with BLAST v2.9.0+ [56] and the blast score ratio (BSR) [57] was calculated. The BSR value is calculated by dividing the query alignment bit score by the reference alignment bit score; a BSR value of 0.8 is approximate to 80% identity over 100% of the alignment length. CDSs unique to each group were identified by the compare_BSR.py script; this requires that a unique region has a BSR value of ≥0.9 in all target genomes and a BSR value of <0.4 in all non-target genomes. The visualization of unique regions was performed by iTOL v5.6.1 [58].

2.9. Design of TaqMan qPCR Assays

Two of the unique F. tularensis regions were of a sufficient length (>300 nts) and associated with coding regions (FTS_0772, FTS_1286). Primers were designed to these targets with Primer3 v2.3.6 [59]; the specificity of primer sequences was determined through the in silico PCR tool, VIPR v0.0.1 (https://github.com/TGenNorth/vipr). The primer sequences were evaluated against NCBI BLAST databases (nr/nt and wgs) to check for potential cross-amplification in other isolated and uncultivated bacteria.

2.10. In Silico Evaluation of Francisella Primers

Previously published primer/probe sets for F. tularensis detection were processed across all analyzed genomes (n = 499) in this study using VIPR. Sensitivity and specificity were calculated using F. tularensis as the target species (Table S2). A subset of the highest scoring assays was selected, and targeted genome regions were extracted from the FNN dataset and aligned to evaluate the theoretical performance and uniqueness of these assays.

2.11. Validation of Francisella tularensis Specific Assays

We tested three potential F. tularensis-specific TaqMan assays on a diverse panel of Francisella DNA samples to validate their specificity to F. tularensis. Two of these assays were previously published [60,61] and a third assay was developed in this study (Table 1). The latter assay (Ft-sp.FTS_0772) targets a F. tularensis-specific signature identified in this study through in silico analysis of diverse Francisella genomes (Table S1). This novel assay was designed as a TaqMan real-time PCR assay using Primer Express software (Thermo Fisher). A minor groove binding moiety was attached to the short probe to increase the probe’s melting temperature and stabilize probe-target hybridization. Each TaqMan real-time PCR assay was run with a 10 µL reaction volume and comprised of the following components: 1× TaqMan Universal PCR master mix (Life Technologies, Applied Biosystems, Foster City, CA), primers (Integrated DNA Technologies, San Diego, CA, USA) and probe (Life Technologies, Applied Biosystems, Foster City, CA) (Table 1), and 1 μL of DNA template (~1 ng). These assays were run on Life Technologies Quant Studio Instruments Flex Real-Time PCR System under conditions described (Table 1) using the following thermal cycle conditions: 50 °C for 2 min, 95 °C for 10 min, and 50 cycles of 95 °C for 15 s and 60 or 63 °C for 1 min (Table 1). Each qPCR experiment included two negative controls and two positive control DNA samples per assay. The sensitivity of the PCR reactions (Ft-sp.FTT0376, Ft-sp.3Pan and Ft-sp.FTS_0772) was evaluated using dilution series with concentrations of F. tularensis DNA between 10 ng and 1 fg.

2.12. Rationale for Francisella Panel Selection

We selected 84 Francisellaceae DNA species/strains (Table S1) to test the specificity of three potential F. tularensis-specific TaqMan assays (Table 1), including 40 F. tularensis and 44 non Francisella tularensis that were selected based on the diversity (genetic and geographic) and availability of DNA samples. Representative strains were included for most groups (Figure 1 and Table 5), but a few groups were not represented due to a lack of availability of DNA samples. The Francisellaceae panel included strains characterized with whole genome sequencing (n = 65) as well as isolates for which whole genome sequences were not available (n = 19) (Table S3). The unsequenced isolates were genetically assigned to species using previously published genotyping methods [62]. To demonstrate that all DNAs included in the panel supported PCR, we first screened the panel using an assay that targets the bacterial 16S rRNA gene [63].

3. Results

3.1. Organisation of Genera and Species Based on Whole-Genome Analysis

Based on the GTDB definition [49], Francisellaceae includes a total of nine genera divided into three major branches (Figure A1). The dataset included in this study covers the entire clade in which the genera Francisella, Parafrancisella, Allofrancisella, and Pseuodfrancisella are located (with addition of the selected outgroup genome of Fangia hongkongensis) and is therefore regarded as the “Francisella near neighbor” (FNN) dataset (Figure 1). The remaining four genera within the Francisellaceae are Caedibacter, represented with three genomes, and the three genera with single sequenced representatives: Fastidiosibacter, QLIT01 (a currently unnamed genus), and M0027 (another currently unnamed genus) [49].
The genus Francisella consists of 11 different species according to the GTDB definition. The classification of F. novicida as a separate species or subspecies of F. tularensis has been debated [19,64,65]. The recognized species F. philomiragia and F. noatunensis have been considered synonymous in GTDB, and according to their classification, F. noatunensis should rather be considered as a subspecies of F. philomiragia. However, as illustrated in Figure 1, F. noatunensis represents a monophyletic sister clade population to F. philomiragia. F. orientalis and GA01-2794 are also regarded as a single species, which correlates with ANI values of 95%.

3.2. Phylogenetic Tree

In the phylogenetic tree (Figure 1), the genus Francisella is divided into four separate clades (clades 1–4) consistent with results obtained from GTDB analysis. Clades 1 and 2 have been previously described [39]. In clade 1, the Francisella isolate TX07-6608 [66] is closely related to F. novicida and F. tularensis within the sub-clade 1.1. In sub-clade 1.2, considerable diversity is identified among FLEs, and F. opportunistica is a new addition. In clade 2, F. salimarina is a new addition and comprises sub-clade 2.2, whereas sub-clade 2.1 includes F. philomiragia, F. noatunensis, F. orientalis, and the new candidate species GA01-2794. Clade 3 contains the single genome of F. endociliophora as the only representative. Clade 4 consists of recently discovered taxa, each of which represents a novel species. The candidate genus Parafrancisella is represented by the single species Parafrancisella adeliensis. The genus Allofrancisella consists of three species, A. inopinata, A. frigidaquae, and A. guangzhouensis, whereas the genus Pseudofrancisella consists of two species: P. aestuarii and F. frigiditurris. Based on the extended Francisellaceae dataset (Figure A1), the order between Parafrancisella and Allofrancisella is uncertain.

3.3. Average Nucleotide Identity

ANI values ranged from 73.6% between Parafrancisella adeliensis and Pseudofrancisella aestuarii to 98.1% between F. tularensis and F. novicida (Figure 2). The mean inter-clade ANI values varied among clades 1–4, with clade 1 most similar to clade 2 (78.7–81.8%) and least similar to clade 4 (77.7–78.9). The intra-clade ANI value varied between 83.8% (clade 4) and 98.1% (clade 1), with the lowest value of 77.6% observed between FLE and P. aestuarii SYW9. Mean ANI values between Francisella and the other genera (Allofrancisella, Parafrancisella and Pseudofrancisella) varied between 73.8% and 76.9%.

3.4. Subspecies and Major Genotypes within Francisella tularensis

Francisella tularensis is divided into three subspecies: tularensis, mediasiatica, and holarctica (Figure 3); F. novicida is regarded as a separate Francisella species and not a subspecies of F. tularensis. F. tularensis subspecies tularensis (Type A) is further separated into two distinct genetic groups, A.I and A.II, with three major branches in A.I [67,68] and two major branches in A.II [69,70]. F. tularensis subspecies mediasiatica is divided into clades, M.I and M.II, wherein M.II is represented by three recently available genomes from Russia [71]. A third clade, M.III, has been proposed but no genomes are available for this clade [5]. F. tularensis subspecies holarctica (Type B) is divided into two major branches B.16 and B.2, wherein the latter is further divided into the four branches B.4, B.6, B.12, and B.159 according to canSNP nomenclature [72]. Most of the known isolates of F. tularensis subspecies holarctica belong to the three branches B.4, B.6, and B.12 [73,74,75,76]. Branch B.16 is usually named “biovar japonica” after where it was first discovered The clade with canSNP designation B.159 is represented by a single strain, F0835, to date only described from California, USA [69]. This strain, phylogenetically assigned to clade B.2 (Figure 3), falls between B.16 and the other three major clades in the F. tularensis subspecies tularensis (B.12, B.6 and B.4) phylogeny.

3.5. Francisella Tularensis Specific Genomic Regions

A total of six unique genomic regions were identified within the F. tularensis core genome that were absent in all near neighbor species, based upon a BSR value ≥ 0.9 in all F. tularensis genomes and a BSR value < 0.4 in all near-neighbor genomes (Table 2). Annotation of these regions, performed using Prodigal, could not match the full coding regions with regions annotated in GenBank; the specific coordinates for each unique region are, however, reported across two reference genomes (Table 3).

3.6. Genomic Regions Specific to Individual Species

To identify species-specific genomic regions, a comparative BSR analysis was performed. The results identify genomic regions unique for each clade (Table 2). The pan genome of fish pathogens (F. orientalis and F. noatunensis) includes the lowest number of unique gene variants, whereas F. novicida and F. philomiragia includes the largest number of unique gene variants.

3.7. Genomic Regions Specific to Francisella Tularensis Subspecies and Major Phylogentic Groups

We also attempted to identify genomic regions unique to each of the three F. tularensis subspecies, as well as major phylogenetic groups within those subspecies. In this analysis, we did not implement a hard threshold for gene absence but, rather, searched for consistently more conserved regions between groups. The results demonstrate variable conservation across the three subspecies (Figure 4, Table 4), suggesting differentiation of these subspecies over time. Evolution in F. tularensis is mainly driven by gene loss [77], so the presence of genes unique to a particular subspecies could be associated with the persistence of ancestral regions in one subspecies that were subsequently lost in the other subspecies.

3.8. In Silico Evaluation of PCR Assay Specificity

PCR primers were screened against all genome assemblies using VIPR, an in silico PCR tool. The number of positive calls was determined by the identification of a predicted amplicon size. A limitation of this approach is that VIPR requires exact matches of primers and probe. In practice, a single nucleotide mismatch may still work in the laboratory. Sensitivity and specificity were calculated using the assumption that assays should only target F. tularensis and miss all near neighbors, including the closely related F. novicida. The results demonstrate that only four assays (Ft-sp.FTT0376, Ft-sp.3Pan, Ft-sp.FTT0523, and Ft-sp.FTS_0772) demonstrated 100% sensitivity and specificity towards F. tularensis (Supplementary Table S3), including the assay designed in this study. The PCR primers designed in this study show no cross-amplification when evaluated against sequences from other bacteria.
The recently published Ft-sp.3Pan assay [61] targets a genomic region that is not unique to F. tularensis. Thus, the specificity of this assay to F. tularensis is based upon utilization of primers designed for a portion of this region that differs between F. tularensis and near-neighbor species by the presence of several SNPs. Of note, the targeted region is closely located to a mobile element (ISFtu1), which could complicate assay evaluation via in silico analyses, especially when using incomplete, unclosed genomes. At least four non-F. tularensis genomes in the FNN dataset contain the targeted region. These genomes (two F. philomiragia, one F. novicida, and one F. salimarina), collectively, have one mismatch in the forward primer, three mismatches in the reverse primer, and one mismatch in the probe sequence. The Ft-spFTT0376 assay [60] targets a region in F. sp. TX07-6608 with conserved primer regions that could potentially cause positive results. Many other Francisellaceae species have alignments to this region but there are no perfect primer matches, so it is unclear if any of these species could cause positive results with this assay. The genomic target of assay Ft-sp.FTT0523 [60] is shared between F. tularensis and some strains of F. novicida, with, collectively, two mismatches in the forward primer, three mismatches in the reverse primer, and five mismatches in probe sequences. It should be noted that only F. tularensis subspecies tularensis subtype A1 has no mismatches in the probe sequence of this assay. F. tularensis subspecies holarctica, F. tularensis subspecies mediasiatica, and F. tularensis subspecies subtype A2 contain one mismatch in the probe sequence.

3.9. Validation of Francisella tularensis Specific Assays

The three potential F. tularensis-specific TaqMan assays all robustly amplified F. tularensis DNA but varied in exclusion of near-neighbors (Table 5). Two assays (Ft-sp.FTT0376 and Ft-sp.FTS_0772) were highly specific to F. tularensis. Among 84 diverse strains in the Francisellaceae DNA panel, all 40 F. tularensis DNAs yielded positive results, whereas all 44 near-neighbor strains yielded negative results with these two assays. These results indicate that both of these assays are both highly sensitive and specific for F. tularensis. The Ft-sp.3pan assay was also very sensitive, yielding positive results for the 40 F. tularensis DNAs, but demonstrated less specificity compared to the other two assays as indicated by robust amplification of DNA from two F. philomiragia near-neighbor strains (Table 5). The failure of this assay on these two near-neighbor DNA samples was not due to poor DNA quality, as they yielded robust amplification with 16S qPCR. Taken together, these data support that two assays, Ft-sp.FTS_0772 and Ft-sp.FTT0376, are highly sensitive and specific to F. tularensis (Table 5). The sensitivity of the three evaluated PCR reactions (Ft-sp.FTT0376, Ft-sp.3Pan and Ft-sp.FTS_0772) are between 1 and 10 fg.

4. Discussion

This study represents the most comprehensive description to date of genomic taxonomy and diversity for the family Francisellaceae. We analyzed a whole genome dataset consisting of 327 F. tularensis genomes and 172 non-F. tularensis Francisellaceae genomes. We recognized nine genera within the family Francisellaceae, including four in the Francisella near neighbor (FNN) dataset (Francisella, Allofrancisella, Pseudofrancisella, and Parafrancisella). Based on genetic distances, the genus Francisella is best divided into four major clades (Figure 1). These new data were used to increase and validate the specificity of DNA-based assays to detect F. tularensis. A genomic region unique to F. tularensis was used to design a novel PCR assay (Ft-sp.FTS_0772) that, against a diagnostic validation panel as well as the in silico testing represented by diversity within the FNN dataset, demonstrated superior or equivalent specificity compared to previously published assays.
The genus Francisella was previously split into two main genetic clades [39], corresponding to clades 1 and 2 in this study. Briefly, F. tularensis strains, which can infect mammals, were assigned to clade 1, whereas F. noatunensis and F. orientalis strains, which can infect fish, were assigned to clade 2 [39]. In the expanded FNN genome dataset reported here, four major clades are recognized within genus Francisella, but it should be noted that the clade division of strains infecting mammals as clade 1 and strains infecting fish as clade 2 still remains valid, although more diversity have been added to both clades. A recent global phylogenetic analysis of the genus Francisella based on 63 genomes [40] divided the genus into three major clades (A-C), roughly corresponding to clades 1, 2, and 4 presented here. However, the rooting used in this analysis resulted in F. persica, F. opportunistica, and F. hispaniensis not being assigned to any major clade, and this analysis did not include all of the species and genera included in this current study.
In the four-clade phylogenetic structure for Francisella presented in this study, clade 1 contains F. tularensis together with F. hispaniensis in one sub-clade and the FLEs and F. opportunistic in the other sub-clade. In general, the majority of clade 1 strains are either endosymbionts (FLEs) of ticks, or facultative intracellular pathogens almost exclusively isolated from humans or other warm-blooded animals. The FLEs only replicate intracellularly and can be transmitted transovarially [38,78,79]. Unlike F. tularensis, FLEs do not grow in cell-free media and their transmission to and virulence in humans is unknown [80,81]. However, the single isolate of F. persica seems to be a transitional form in terms of ability to grow on cell-free media, since it can be cultivated on agar plates but grows extremely slowly (14–20 days). On the other hand, not all members of clade 1 originate from mammals or vectors. Exceptions include a new member of clade 1, Francisella isolate TX07-6608, that was directly isolated from seawater [66] and the first isolate of F. novicida, an exceedingly rare opportunistic human pathogen, that was directly isolated also from salt water [15]. Clearly, F. novicida show distinct differences to F. tularensis in clinical, ecological, genomic and virulence properties [82]. This, in combination with the genome clustering of the 37 strains of F. novicida shown in this study, supports maintaining F. novicida and F. tularensis as separate species.
Although none of the three human pathogenic F. tularensis subspecies have been directly associated with seawater, this capability is present in other members of clade 1. A rare case of F. hispaniensis infection was acquired from seawater [28]. Likewise, the first reported Francisella strain isolated from the Southern Hemisphere originated from a patient infected in brackish water in Australia (F. hispaniensis strain 3523) [39,83]. Interestingly, it has been shown that F. tularensis subspecies tularensis shows enhanced survival in brackish water, suggesting that salt-influenced (marine) environments may promote its survival [84]. Together, these findings are in agreement with a previous suggestion that F. tularensis may have evolved in and from the marine environment [42].
Clade 2 contains the fish pathogens F. noatunensis and F. orientalis, as well as F. philomiragia and F. salimarina. Strains of F. philomiragia cause disease primarily in immunocompromised hosts and/or salt-water near-drowning victims [18,85]. All known strains of F. salimarina have been directly isolated from seawater and, as such, their potential clinical relevance is as yet unknown [86]. Although very rare, the majority of human infections caused by non-F. tularensis members of clade 1 (F. novicida, F. hispaniensis) and some members of clade 2 (F. philomiragia) tend to occur in immunocompromised individuals following exposure to salt or brackish water [87].
Clade 3 is represented by a single isolate of F. endociliophora, an endocytobiont of the marine ciliate species Euplotes [88]. F. endociliophora was directly isolated and cultured in cell free media at room temperature from E. raikovi, a temperate water ciliate species [20]. Another endocytobiont of Euplotes, P. adeliensis, the only isolate of the candidate genus Parafrancisella, was isolated in a similar manner from E. petzi [21]. It is well established that the virulent F. tularensis subspecies are deficient in biosynthesis pathways for many amino acids, a trait that is shared by both these endocytobionts (F. endociliophora and P. adeliensis), suggesting a shared host dependency for supply of amino acids. This further strengthens the assumption of a marine origin for F. tularensis. Strains that constitute clade 4 in the present dataset (this study) are all non-clinical isolates, including isolates from seawater, cooling water systems, and abalone [23,26]. Thus, clade 4 contains what appears to be free-living bacteria in the environment.
Within the dataset analyzed (499 genomes), we identified only six short sequence regions apparently unique to F. tularensis. We used a strict definition for a region to be called unique, which may have limited the comprehensive identification of unique signatures. Another explanation why there are so few regions unique to the virulent F. tularensis is that the evolution of F. tularensis is characterized by genome erosion (Figure 4), so that rare unique sequence regions are remnants that have persisted through the evolution of F. tularensis, which suggest that these regions may play an important role to the lifecycle of F. tularensis. Two of the regions are intergenic and three are partial loci of either pseudogenes or hypothetical proteins. The 6th, region 4, constitutes a partial domain of the adenosine deaminase gene (FTT 0939), a gene encoding a key enzyme in purine metabolism, but this gene is predicted to be functional only in F. tularensis subspecies tularensis [89,90]. Most bacterial pathogens causing systemic infections seems to lack a functional adenosine deaminase gene, an absence that is probably compensated by host metabolic pathways [91]. The reason for the presence of an apparently intact adenosine deaminase gene only in F. tularensis subspecies tularensis is unknown.
The main concern with F. tularensis-specific assays is not false-negatives, as F. tularensis is well characterized and highly monomorphic. Rather, the challenge is the unknown diversity that is assumed to exist amongst uncultivable bacterial species and especially environmental members within clade 1 that remain to be explored, as indicated by the recently characterized F. opportunistica and F. hispaniensis-like strains (Figure 1). Region 1 is the target for the new PCR assay developed in this study (Ft-sp.FTS_0772) that is highly specific for F. tularensis. This unique sequence region is located in the terminal remnant of the hsdS gene and originally encoded a type 1 restriction enzyme, congruent with the gene erosion, resulting in the significant loss of restriction capability in virulent F. tularensis [92]. We also screened previously published PCR assay primers that were developed for specific identification of F. tularensis against all genome assemblies in our dataset. Only two assays, Ft-sp3Pan [61] and Ft-sp.FTT0376 [60], demonstrated high sensitivity and specificity for F. tularensis in these in silico analyses. However, a more detailed alignment showed theoretical specificity issues with both of these assays. It should be noted that these alignments per se would not disqualify these assays; high specificity could in practice still be obtained. Therefore, these two assays were included, as well as the assay designed in this study, in the wet-bench validation against the Francisella strain diversity panel. The results reveal that two of the assays, Ft-sp.FTS_0772 (this study) and Ft-sp.FTT0376 [60], are both highly sensitive and specific to F. tularensis. A caveat in this validation was, however, the lack of strain TX07-6608 in the panel. This taxon could, at least theoretically, compromise the specificity of the Ft-sp.FTT0376 assay. Another finding was that the Ft-sp 3pan assay did not exclude all near-neighbors on our diversity panel, as two F. philomiragia strains (F1091 and F1093) amplified robustly. In silico analysis revealed that strains F1091 and F1093 were the only examples of the tested 23 F. philomiragia strains in the wet-bench validation panel that harbor the target genome region of the Ft-sp 3pan assay. The targeted region is closely located to a mobile element (ISFtu1), which may explain the absence of this region in many of the strains. Taken together, we would rather recommend using the PCR assay designed in this study, or Ft-sp.FTT0376, for detection of F. tularensis in environmental and other complex samples.
This study dramatically expands the genome space of Francisellaceae. Previous diagnostics for F. tularensis have suffered from a largely shared but unexplored Francisellaceae pan genome, resulting in false positive results, which is problematic for routine biothreat surveillance. The unique signatures identified in this study may also be shared by unexplored, unculturable, and previously undetected Francisella near neighbor species that have yet to be characterized. Continued culture-independent characterization of these cryptic species will expand our understanding of the diversity of these lineages.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/9/1/146/s1, Table S1: Metadata describing included Francisellaceae isolates, Table S2: In silico validation of PCR assays, Table S3: Metadata and performance of validation panel, Table S4: Recommended genome assemblies for Francisellaceae species and subspecies.

Author Contributions

Conceptualization, D.M.W., J.W.S., M.F. and A.S.; methodology, J.W.S., D.B., C.Ö. and A.S.; software, J.W.S., C.Ö. and A.S.; validation, J.W.S., D.B., R.B., K.S., R.F.M., S.B. and L.K.; formal analysis, J.W.S., D.B., R.B., C.Ö., I.U., E.H., A.S. and J.T.; investigation, D.B., R.B., L.K., S.B. and J.T.; resources, S.S., K.C., J.G.R.-P., A.J. and A.V.; data curation, J.W.S., D.B., K.S., C.Ö., A.S. and I.U.; writing—original draft preparation, J.W.S., D.B., A.S., C.Ö. and M.F.; writing—review and editing, D.M.W., J.W.S., D.B., S.S., A.S., C.Ö., M.F., D.S., B.B., J.T., D.C., K.M. and A.J.; visualization, J.W.S., D.S. and C.Ö.; supervision, D.M.W., A.S. and M.F.; project administration: D.M.W., A.S., M.F.; funding acquisition: D.M.W., S.S., A.J., M.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Swedish Civil Contingencies Agency [grant number TA 014-2010-01] and the US Department of Homeland Security’s Science and Technology Directorate [award number HSHQDC-17-C-B0021] pursuant to the agreement between the Kingdom of Sweden and the US government on Cooperation in Science and Technology for Homeland Security Matters.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data generated in this study are openly available in NCBI project accession PRJNA657836. Public available genome data presented in this study is listed in Supplementary Table S1.

Acknowledgments

Shotgun sequencing was performed by the SNP&SEQ Technology Platform in Uppsala. The facility is part of the National Genomics Infrastructure (NGI) Sweden and Science for Life Laboratory. The SNP&SEQ Platform is supported by the Swedish Research Council and the Knut and Alice Wallenberg Foundation.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table
Table A1. Detailed overview of species within family Francisellaceae following NCBI taxonomy.
Table A1. Detailed overview of species within family Francisellaceae following NCBI taxonomy.
Francisellaceae genusFrancisellaceae speciesHeterotypic SynonymsType StrainsReferences
FrancisellaFrancisella tularensis [93]
Francisella tularensis subspecies tularensis SchuS4[17,42,94,95]
Francisella tularensis subspecies holarcticaFrancisella tularensis var. palaearcticaFSC257; GIEM 503[17,95,96]
FSC200
Francisella tularensis subspecies mediasiatica FSC147[41,95,97,98]
Francisella tularensis subspecies novicidaF. novicida, Pasteurella novicidaATCC 15482; CCUG 33449; CIP 56.12[15,17,19]
U112
Francisella hispaniensis CCUG 58020; DSM 22475; F62; FhSp1; FnSp1; FSC454[19]
Francisella persica (FLE)Wolbachia persicaATCC VR-331; DSM 101678; FSC845[22,99,100,101]
Francisella opportunistica ATCC BAA-2974; DSM 107100; PA05-1188[23,102,103]
Francisella philomiragiaF. philomiragia subspecies philomiragia, Yersinia philomiragiaATCC 25015; CCUG 19700; CCUG 4992; CIP 82.98; DSM 7535[18,82,104]
Francisella noatunensisF. piscicida; F. philomiragia subspecies noatunensis Francisella noatunensis subspecies noatunensis, 2005/50/F292-6; DSM 23596; LMG 23800; NCIMB 14265[11,12,13,105]
Francisella orientalisF. victoria, F. asiatica, Francisella noatunensis subspecies orientalisDSM 21254; Ehime-1; LMG 24544[11,12,13]
CCUG 56120T; DSM 21294; PQ 1104
Francisella salina TX07-7308[23,66]
Francisella marina E103-15, E95-16[9]
Francisella salimarina CGMCC 1.17031; NBRC 113781; SYSU SYW-1[83]
Francisella endociliophoraF. noatunensis subspecies endociliophora [20,86]
Francisella halioticida DSM 23729; LMG 26062; Shimane-1[10,26]
Francisella sp. LA11-2445 LA11-2445[25]
Francisella uliginis TX07-7310[23,66]
Francisella adeliensis FSC1327[21]
Francisella frigiditurris CA97-1460[23]
AllofrancisellaAllofrancisella inopinata DSM 101834; KCTC 42968; SYSU 23; SYSU YG23[1,27]
Allofrancisella frigidaquae DSM 101835; KCTC 42969; SYSU 10HL1970[1,24,27]
Allofrancisella guanzhousensisF. guanzhousensis, F. cantonensis08HL01032; ATCC BAA-2361; CCUG 60119; NCTC 13503[1,31,106]
PseudofrancisellaPseudofrancisella aestuarii CGMCC 1.13718; KCTC 52557; SYSU WZ-2[2]
Microorganisms 09 00146 g0a1 550
Figure A1. Phylogenetic tree of family Francisellaceae.
Figure A1. Phylogenetic tree of family Francisellaceae.
Microorganisms 09 00146 g0a1

References

  1. Qu, P.-H.; Li, Y.; Salam, N.; Chen, S.; Liu, L.; Gu, Q.; Fang, B.-Z.; Xiao, M.; Li, M.; Chen, C.; et al. Allofrancisella inopinata gen. nov., sp. nov. and Allofrancisella frigidaquae sp. nov., isolated from water-cooling systems, and transfer of Francisella guangzhouensis Qu et al. 2013 to the new genus as Allofrancisella guangzhouensis comb. nov. Int. J. Syst. Evol. Microbiol. 2016, 66, 4832–4838. [Google Scholar] [CrossRef] [PubMed]
  2. Zheng, M.-L.; Jiao, J.-Y.; Dong, L.; Han, M.-X.; Li, L.-H.; Xiao, M.; Chen, C.; Qu, P.-H.; Li, W.-J. Pseudofrancisella aestuarii gen. nov., sp. nov., a novel member of the family Francisellaceae isolated from estuarine seawater. Antonie Van Leeuwenhoek 2019, 112, 877–886. [Google Scholar] [CrossRef] [PubMed]
  3. Colquhoun, D.J.; Larsson, P.; Duodu, S.; Forsman, M. The Family Francisellaceae. Prokaryotes 2014, 287–314. Available online: https://link.springer.com/referenceworkentry/10.1007/978-3-642-38922-1_236 (accessed on 1 December 2020).
  4. Penn, R.L. Francisella tularensis (Tularemia). In Mandell, Douglas, and Bennett’s Principles and Practice of Infectious Diseases; Elsevier: Amsterdam, The Netherlands, 2010; pp. 2927–2937. [Google Scholar]
  5. Timofeev, V.; Bakhteeva, I.; Titareva, G.; Kopylov, P.; Christiany, D.; Mokrievich, A.; Dyatlov, I.; Vergnaud, G. Russian isolates enlarge the known geo-graphic diversity of Francisella tularensis subsp. mediasiatica. PLoS ONE 2017, 12, e0183714. Available online: http://www.ncbi.nlm.nih.gov/pubmed/28873421 (accessed on 1 November 2020). [CrossRef] [PubMed] [Green Version]
  6. Timofeev, V.; Titareva, G.; Bahtejeva, I.; Kombarova, T.; Kravchenko, T.; Mokrievich, A.; Dyatlov, I. The Comparative Virulence of Fran-cisella tularensis Subsp. mediasiatica for Vaccinated Laboratory Animals. Microorganisms 2020, 8, 1403. Available online: https://www.mdpi.com/2076-2607/8/9/1403 (accessed on 1 November 2020). [CrossRef] [PubMed]
  7. Dennis, D.T.; Inglesby, T.V.; Henderson, D.A.; Bartlett, J.G.; Ascher, M.S.; Eitzen, E.; Fine, A.D.; Friedlander, A.M.; Hauer, J.; Layton, M.; et al. Tularemia as a biological weapon: Medical and public health management. JAMA Am. Med Assoc. 2001, 285, 2763–2773. Available online: http://jama.ama-assn.org/content/285/21/2763.abstract (accessed on 1 November 2020). [CrossRef] [PubMed]
  8. Maurin, M. Francisella tularensis as a potential agent of bioterrorism? Expert Rev. Anti. Infect. Ther. 2015, 13, 141–144. Available online: https://pubmed.ncbi.nlm.nih.gov/25413334/ (accessed on 1 November 2020).
  9. Soto, E.; Griffin, M.J.; Morales, J.A.; Calvo, E.B.; Sebastião, F.D.A.; Porras, A.L.; Víquez-Rodríguez, X.; Reichley, S.R.; Rosser, T.G.; Ware, C.; et al. Francisella marinasp. nov., Etiologic Agent of Systemic Disease in Cultured Spotted Rose Snapper (Lutjanus guttatus) in Central America. Appl. Environ. Microbiol. 2018, 84, 1–16. [Google Scholar] [CrossRef] [Green Version]
  10. Brevik, Ø.J.; Ottem, K.; Kamaishi, T.; Watanabe, K.; Nylund, A. Francisella halioticida sp. nov., a pathogen of farmed giant abalone (Haliotis gigantea) in Japan. J. Appl. Microbiol. 2011, 111, 1044–1056. [Google Scholar] [CrossRef]
  11. Ottem, K.; Nylund, A.; Karlsbakk, E.; Friis-Møller, A.; Kamaishi, T. Elevation of Francisella philomiragia subsp. noatunensis Mikalsen et al. (2007) to Francisella noatunensis comb. nov.[syn. Francisella piscicida Ottem et al.(2008) syn. nov.] and characterization of Francisella noatunensis subsp. orientalis subsp. nov., two important fish pathogens. J. Appl. Microbiol. 2009, 106, 1231–1243. Available online: http://www3.interscience.wiley.com/cgi-bin/abstract/121660927/ABSTRACT (accessed on 1 November 2020).
  12. Mikalsen, J.; Colquhoun, D.J. Francisella asiatica sp. nov. isolated from farmed tilapia (Oreochromis sp.) and elevation of Francisella philomiragia subsp. noatunensis to species rank as Francisella noatunensis comb. nov., sp. nov. Int. J. Syst. Evol. Microbiol. 2009. Available online: http://ijs.sgmjournals.org/content/ijs.0.002139-0v1/.abstract (accessed on 1 November 2020). [CrossRef] [PubMed]
  13. Ramirez-Paredes, J.G.; Larsson, P.; Thompson, K.D.; Penman, D.J.; Busse, H.-J.; Öhrman, C.; Sjödin, A.; Soto, E.; Richards, R.H.; Adams, A.; et al. Reclassification of Francisella noatunensis subsp. orientalis Ottem et al. 2009 as Francisella orientalis sp. nov., Francisella noatunensis subsp. chilensis subsp. nov. and emended description of Francisella noatunensis. Int. J. Syst. Evol. Microbiol. 2020, 70, 2034–2048. [Google Scholar] [CrossRef] [PubMed]
  14. Colquhoun, D.J.; Duodu, S. Francisella infections in farmed and wild aquatic organisms. Veter-Res. 2011, 42, 47. [Google Scholar] [CrossRef] [Green Version]
  15. Larson, C.L.; Wicht, W.; Jellison, W.L. A new organism resembling P. tularensis isolated from water. Public Health Rep. 1955, 70, 253–258. Available online: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2024510/ (accessed on 1 November 2020). [CrossRef]
  16. Aravena-Román, M.; Merritt, A.J.; Inglis, T.J. First case of Francisella bacteraemia in Western Australia. New Microbes New Infect. 2015, 8, 75–77. [Google Scholar] [CrossRef] [Green Version]
  17. Olsufiev, N.G.; Emelyanova, O.S.; Dunayeva, T.N. Comparative study of strains of B. tularense in the old and new world and their taxonomy. J. Hyg. Epidemiol. Microbiol. Immunol. 1959, 3, 138–149. [Google Scholar]
  18. Hollis, D.G.; Weaver, R.; Steigerwalt, A.G.; Wenger, J.D.; Moss, C.W.; Brenner, D.J. Francisella philomiragia comb. nov. (formerly Yersinia philomiragia) and Francisella tularensis biogroup novicida (formerly Francisella novicida) associated with human disease. J. Clin. Microbiol. 1989, 27, 1601–1608. [Google Scholar] [CrossRef] [Green Version]
  19. Huber, B.; Escudero, R.; Busse, H.-J.; Seibold, E.; Scholz, H.C.; Anda, P.; Kämpfer, P.; Splettstoesser, W.D. Description of Francisella hispaniensis sp. nov., isolated from human blood, reclassification of Francisella novicida (Larson et al. 1955) Olsufiev et al. 1959 as Francisella tularensis subsp. novicida comb. nov. and emended description of the genus Francisella. Int. J. Syst. Evol. Microbiol. 2010, 60, 1887–1896. [Google Scholar] [CrossRef]
  20. Sjödin, A.; Öhrman, C.; Bäckman, S.; Lärkeryd, A.; Granberg, M.; Lundmark, E.; Karlsson, E.; Nilsson, E.; Vallesi, A.; Tellgren-Roth, C.; et al. Complete Genome Sequence of Francisella endociliophora Strain FSC1006, Isolated from a Laboratory Culture of the Marine Ciliate Euplotes raikovi. Genome Announc. 2014, 2, e01227-14. [Google Scholar] [CrossRef] [Green Version]
  21. Vallesi, A.; Sjödin, A.; Petrelli, D.; Luporini, P.; Taddei, A.R.; Thelaus, J.; Öhrman, C.; Nilsson, E.; Di Giuseppe, G.; Gutiérrez, G.; et al. A New Species of the γ-Proteobacterium Francisella, F. adeliensis sp. Nov., Endocytobiont in an Antarctic Marine Ciliate and Potential Evolutionary Forerunner of Pathogenic Species. Microb. Ecol. 2019, 77, 587–596. [Google Scholar] [CrossRef]
  22. Larson, M.A.; Nalbantoglu, U.; Sayood, K.; Zentz, E.B.; Cer, R.Z.; Iwen, P.C.; Francesconi, S.C.; Bishop-Lilly, K.A.; Mokashi, V.P.; Sjöstedt, A.; et al. Reclassification of Wolbachia persica as Francisella persica comb. nov. and emended description of the family Francisellaceae. Int. J. Syst. Evol. Microbiol. 2016, 66, 1200–1205. [Google Scholar] [CrossRef] [PubMed]
  23. Challacombe, J.F.; Petersen, J.M.; Gallegos-Graves, L.V.; Hodge, D.; Pillai, S.; Kuske, C.R. Whole-Genome Relationships among Francisella Bacteria of Diverse Origins Define New Species and Provide Specific Regions for Detection. Appl. Environ. Microbiol. 2016, 83, 1–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Rydzewski, K.; Schulz, T.; Brzuszkiewicz, E.; Holland, G.; Lück, C.; Fleischer, J.; Grunow, R.; Heuner, K. Genome sequence and phenotypic analysis of a first German Francisella sp. isolate (W12-1067) not belonging to the species Francisella tularensis. BMC Microbiol. 2014, 14, 169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Respicio-Kingry, L.; Byrd, L.; Allison, A.; Brett, M.; Scott-Waldron, C.; Galliher, K.; Hannah, P.; Mead, P.; Petersen, J.M. Cutaneous Infection Caused by a Novel Francisella sp. J. Clin. Microbiol. 2013, 51, 3456–3460. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Uneklint, I.; Öhrman, C.; Karlsson, L.; Byström, M.; Hägglund, E.; Forsman, M.; Sjödin, A. Complete Genome Sequence of Francisella hali-oticida Type Strain DSM 23729 (FSC1005). Microbiol. Resour. Announc. 2020, 9, e01418-17. Available online: https://mra.asm.org/content/9/37/e00541-20 (accessed on 1 November 2020). [CrossRef] [PubMed]
  27. Öhrman, C.; Uneklint, I.; Karlsson, L.; Svensson, D.; Forsman, M.; Sjödin, A. Complete Genome Sequences of Allofrancisella inopinata SYSU YG23 and Allofrancisella frigidaquae SYSU 10HL1970, Isolated from Water from Cooling Systems in China. Microbiol. Resour. Announc. 2020, 9. Available online: https://mra.asm.org/content/9/48/e00554-20 (accessed on 1 November 2020).
  28. Zhou, H.; Yang, Q.; Shen, L.; Yao, Y.; Xu, J.; Ye, J.; Wu, X.; Yu, Y.; Li, Z.; Zhou, J.; et al. Seawater-Associated Highly Pathogenic Francisella hispaniensis Infections Causing Multiple Organ Failure. Emerg. Infect. Dis. 2020, 26, 2424–2428. [Google Scholar] [CrossRef]
  29. Whitehouse, C.; Kesterson, K.; Duncan, D.; Eshoo, M.; Wolcott, M. Identification and characterization of Francisella species from natural warm springs in Utah, USA. Lett. Appl. Microbiol. 2012, 54, 313–324. [Google Scholar] [CrossRef]
  30. Brett, M.; Respicio-Kingry, L.B.; Yendell, S.; Ratard, R.; Hand, J.; Balsamo, G.; Scott-Waldron, C.; O’Neal, C.; Kidwell, D.; Yockey, B.; et al. Outbreak of Francisella novicida Bacteremia Among Inmates at a Louisiana Correctional Facility. Clin. Infect. Dis. 2014, 59, 826–833. [Google Scholar] [CrossRef] [Green Version]
  31. Qu, P.-H.; Chen, S.; Scholz, H.C.; Busse, H.-J.; Gu, Q.; Kämpfer, P.; Foster, J.T.; Glaeser, S.P.; Chen, C.; Yang, Z.-C. Francisella guangzhouensis sp. nov., isolated from air-conditioning systems. Int. J. Syst. Evol. Microbiol. 2013, 63, 3628–3635. [Google Scholar] [CrossRef]
  32. Person, T.; Currie, C. Biosurveillance: DHS Should Not Pursue BioWatch Upgrades or Enhancements Until System Capabilities Are Established; United States Government Accountability Office: Washington, DC, USA, 2015. Available online: https://www.gao.gov/assets/680/673309.pdf (accessed on 1 November 2020).
  33. Friend, M. Tularemia. Acta Vet. Scand. 2016. Available online: http://wildpro.twycrosszoo.org/000ADOBES/Lagomorphs/B527Circ1297Tularemia.pdf (accessed on 1 November 2020).
  34. Parker, R.R.; Steinhaus, E.; Kohls, G.M.; Jellison, W.L. Contamination of natural waters and mud with Pasteurella tularensis and tularemia in beavers and muskrats in the northwestern United States. Bull. Natl. Institutes Heal. (U.S.) 1951, 193, 1–161. [Google Scholar]
  35. Karadenizli, A.; Forsman, M.; Şimşek, H.; Taner, M.; Öhrman, C.; Myrtennäs, K.; Lärkeryd, A.; Johansson, A.; Özdemir, L.; Sjödin, A. Genomic analyses of Francisella tularensis strains confirm disease transmission from drinking water sources, Turkey, 2008, 2009 and 2012. Eurosurveillance 2015, 20. Available online: http://www.ncbi.nlm.nih.gov/pubmed/26062561 (accessed on 1 November 2020). [CrossRef] [PubMed]
  36. Kuske, C.R.; Barns, S.M.; Grow, C.C.; Merrill, L.; Dunbar, J. Environmental Survey for Four Pathogenic Bacteria and Closely Related Species Using Phylogenetic and Functional Genes*. J. Forensic. Sci. 2006, 51, 548–558. [Google Scholar] [CrossRef] [PubMed]
  37. Barns, S.M.; Grow, C.C.; Okinaka, R.T.; Keim, P.; Kuske, C.R. Detection of diverse new Francisella-like bacteria in environmental sam-ples. Appl. Environ. Microbiol. 2005, 71, 5494–5500. Available online: http://aem.asm.org/content/71/9/5494.abstract (accessed on 1 November 2020). [CrossRef] [PubMed] [Green Version]
  38. Gerhart, J.G.; Auguste Dutcher, H.; Brenner, A.E.; Moses, A.S.; Grubhoffer, L.; Raghavan, R. Multiple Acquisitions of Patho-gen-Derived Francisella Endosymbionts in Soft Ticks. Genome Biol. Evol. 2018, 10, 607–615. [Google Scholar] [CrossRef] [Green Version]
  39. Sjödin, A.; Svensson, K.; Öhrman, C.; Ahlinder, J.; Lindgren, P.; Duodu, S.; Johansson, A.; Colquhoun, D.J.; Larsson, P.; Forsman, M. Genome characterisation of the genus Francisella reveals insight into similar evolutionary paths in pathogens of mammals and fish. BMC Genom. 2012, 13, 268. [Google Scholar] [CrossRef] [Green Version]
  40. Kumar, R.; Bröms, J.E.; Sjöstedt, A. Exploring the Diversity within the Genus Francisella—An Integrated Pan-Genome and Ge-nome-Mining Approach. Front. Microbiol. 2020, 11. Available online: https://www.frontiersin.org/article/10.3389/fmicb.2020.01928/full (accessed on 1 November 2020). [CrossRef]
  41. Larsson, P.; Elfsmark, D.; Svensson, K.; Wikström, P.; Forsman, M.; Brettin, T.; Keim, P.; Johansson, A. Molecular Evolutionary Consequences of Niche Restriction in Francisella tularensis, a Facultative Intracellular Pathogen. PLoS Pathog. 2009, 5, e1000472. [Google Scholar] [CrossRef] [Green Version]
  42. Larsson, P.; Oyston, P.C.F.; Chain, P.; Chu, M.C.; Duffield, M.; Fuxelius, H.-H.; Garcia, E.; Hälltorp, G.; Johansson, D.; Isherwood, K.E.; et al. The complete genome sequence of Francisella tu-larensis, the causative agent of tularemia. Nat. Genet. 2005, 37, 153–159. Available online: http://www.ncbi.nlm.nih.gov/pubmed/15640799 (accessed on 1 November 2020). [CrossRef] [Green Version]
  43. Sampson, T.R.; Weiss, D.S. Degeneration of a CRISPR/Cas system and its regulatory target during the evolution of a pathogen. RNA Biol. 2013, 10, 1618–1622. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Schunder, E.; Rydzewski, K.; Grunow, R.; Heuner, K. First indication for a functional CRISPR/Cas system in Francisella tularensis. Int. J. Med. Microbiol. 2013, 303, 51–60. [Google Scholar] [CrossRef] [Green Version]
  45. Kozarewa, I.; Turner, D.J. Amplification-Free Library Preparation for Paired-End Illumina Sequencing. Comput. Biol. 2011, 733, 257–266. [Google Scholar]
  46. Jackman, S.D.; Vandervalk, B.P.; Mohamadi, H.; Chu, J.; Yeo, S.; Hammond, S.A.; Jahesh, G.; Khan, H.; Coombe, L.; Warren, R.L.; et al. ABySS 2.0: Resource-efficient assembly of large genomes using a Bloom filter. Genome Res. 2017, 27, 768–777. [Google Scholar] [CrossRef] [Green Version]
  47. Walker, B.J.; Abeel, T.; Shea, T.; Priest, M.; Abouelliel, A.; Sakthikumar, S.; Cuomo, C.A.; Zeng, Q.; Wortman, J.; Young, S.K.; et al. Pilon: An Integrated Tool for Comprehensive Microbial Variant Detection and Genome Assembly Improvement. PLoS ONE 2014, 9, e112963. [Google Scholar] [CrossRef]
  48. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. Available online: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3342519/ (accessed on 1 November 2020). [CrossRef] [Green Version]
  49. Parks, D.H.; Chuvochina, M.; Chaumeil, P.-A.; Rinke, C.; Mussig, A.J.; Hugenholtz, P. A complete domain-to-species taxonomy for Bacteria and Archaea. Nat. Biotechnol. 2020, 38, 1–8. [Google Scholar] [CrossRef] [PubMed]
  50. Darling, A.E.; Mau, B.; Perna, N.T. progressiveMauve: Multiple Genome Alignment with Gene Gain, Loss and Rearrangement. PLoS ONE 2010, 5, e11147. [Google Scholar] [CrossRef] [Green Version]
  51. Page, A.J.; Taylor, B.; Delaney, A.J.; Soares, J.; Seemann, T.; Keane, J.A.; Harris, S.R. SNP-sites: Rapid efficient extraction of SNPs from multi-FASTA alignments. Microb. Genom. 2016, 2, e000056. [Google Scholar] [CrossRef] [Green Version]
  52. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  53. Pritchard, L.; Glover, R.H.; Humphris, S.; Elphinstone, J.G.; Toth, I.K. Genomics and taxonomy in diagnostics for food security: Soft-rotting enterobacterial plant pathogens. Anal. Methods R. Soc. Chem. 2016, 8, 12–24. Available online: www.rsc.org/methods (accessed on 1 November 2020). [CrossRef]
  54. Hyatt, D.; Chen, G.-L.; Locascio, P.F.; Land, M.L.; Larimer, F.W.; Hauser, L.J. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010, 11, 119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 2012, 28, 3150–3152. [Google Scholar] [CrossRef] [PubMed]
  56. Camacho, C.E.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.S.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Rasko, D.; Myers, G.S.A.; Ravel, J. Visualization of comparative genomic analyses by BLAST score ratio. BMC Bioinform. 2005, 6, 2. [Google Scholar] [CrossRef] [Green Version]
  58. Letunic, I.; Bork, P. Interactive tree of life (iTOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016, 44, W242–W245. [Google Scholar] [CrossRef]
  59. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S. Primer3—New capabilities and interfaces. Nucleic Acids Res. 2012, 40, e115. [Google Scholar] [CrossRef] [Green Version]
  60. Mitchell, J.L.; Chatwell, N.; Christensen, D.; Diaper, H.; Minogue, T.D.; Parsons, T.M.; Walker, B.; Weller, S.A. Development of real-time PCR assays for the specific detection of Francisella tularensis ssp. tularensis, holarctica and mediaasiatica. Mol. Cell. Probes 2010, 24, 72–76. [Google Scholar] [CrossRef]
  61. Larson, M.A.; Sayood, K.; Bartling, A.M.; Meyer, J.R.; Starr, C.; Baldwin, J.; Dempsey, M.P. Differentiation of Francisella tularensis Subspecies and Subtypes. J. Clin. Microbiol 2020, 58, 1–16. Available online: https://jcm.asm.org/content/58/4/e01495-19 (accessed on 1 November 2020). [CrossRef] [Green Version]
  62. Birdsell, D.N.; Vogler, A.J.; Buchhagen, J.; Clare, A.; Kaufman, E.; Naumann, A.; Driebe, E.; Wagner, D.M.; Keim, P.S. TaqMan Real-Time PCR Assays for Sin-gle-Nucleotide Polymorphisms Which Identify Francisella tularensis and Its Subspecies and Subpopulations. PLoS ONE 2014, 9, e107964. Available online: http://www.ncbi.nlm.nih.gov/pubmed/25238067 (accessed on 24 September 2014). [CrossRef]
  63. Liu, C.M.; Aziz, M.; Kachur, S.; Hsueh, P.-R.; Huang, Y.-T.; Keim, P.; Price, L.B. BactQuant: An enhanced broad-coverage bacterial quantitative real-time PCR assay. BMC Microbiol. 2012, 12, 56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Johansson, A.; Celli, J.; Conlan, W.; Elkins, K.L.; Forsman, M.; Keim, P.S.; Larsson, P.; Manoil, C.; Nano, F.E.; Petersen, J.M.; et al. Objections to the transfer of Francisella novicida to the subspecies rank of Francisella tularensis. Int. J. Syst. Evol. Microbiol. 2010, 60, 1717–1718. [Google Scholar] [CrossRef] [PubMed]
  65. Busse, H.-J.; Scholz, H.C.; Escudero, R.; Anda, P.; Kämpfer, P.; Splettstoesser, W.D.; Huber, B.; Seibold, E. Objections to the transfer of Francisella novicida to the subspecies rank of Francisella tularensis—Response to Johansson et al. Int. J. Syst. Evol. Microbiol. 2010, 60, 1718–1720. [Google Scholar] [CrossRef]
  66. Petersen, J.M.M.; Carlson, J.; Yockey, B.; Pillai, S.; Kuske, C.; Garbalena, G.; Pottumarthy, S.; Chalcraft, L. Direct isolation of Francisella spp. from environmental samples. Lett. Appl. Microbiol. 2009, 48, 663–667. Available online: http://www.ingentaconnect.com/content/bsc/lappm/2009/00000048/00000006/art00003 (accessed on 1 November 2020). [PubMed]
  67. Birdsell, D.N.; Johansson, A.; Öhrman, C.; Kaufman, E.; Molins, C.; Pearson, T. Francisella tularensis subsp. tularensis Group A.I, United States. Emerg. Infect. Dis. 2014, 20, 861–865. Available online: http://wwwnc.cdc.gov/eid/article/20/5/13-1559_article.htm (accessed on 1 November 2020). [CrossRef] [Green Version]
  68. Svensson, K.; Granberg, M.; Karlsson, L.; Neubauerova, V.; Forsman, M.; Johansson, A. A real-time PCR array for hierarchical identification of Francisella isolates. PLoS ONE 2009, 4, e8360. Available online: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2793073&tool=pmcentrez&rendertype=abstract (accessed on 1 November 2020). [CrossRef]
  69. Vogler, A.J.; Birdsell, D.; Price, L.B.; Bowers, J.R.; Beckstrom-Sternberg, S.M.; Auerbach, R.K.; Beckstrom-Sternberg, J.S.; Johansson, A.; Clare, A.; Buchhagen, J.L.; et al. Phylogeography of Francisella tu-larensis: Global expansion of a highly fit clone. J. Bacteriol. 2009, 191, 2474–2484. Available online: http://jb.asm.org/content/191/8/2474.abstract (accessed on 1 November 2020). [CrossRef] [Green Version]
  70. Birdsell, D.N.; Yaglom, H.; Rodriguez, E.; Engelthaler, D.M.; Maurer, M.; Gaither, M.; Vinocur, J.; Weiss, J.; Terriquez, J.; Komatsu, K.; et al. Phylogenetic Analysis of Francisella tularensis Group A.II Isolates from 5 Patients with Tularemia, Arizona, USA, 2015–2017. Emerg. Infect. Dis. 2019, 25, 944–946. [Google Scholar] [CrossRef]
  71. Mokrievich, A.N.; Kislichkina, A.A.; Kudryavtseva, T.Y.; Mironova, R.I.; Vakhrameeva, G.M.; Shishkova, N.A.; Timofeev, V.S.; Bogun, A.G.; Pavlov, V.M.; Dyatlov, I.A.; et al. Draft Genome Sequences of Three Francisella tularensis subsp. mediasiatica Strains Isolated in the Altai Territory, Russian Federation. Microbiol Resour Announc [Internet]. American Society for Microbiology; 2020 [cited 2020 Dec 2];9. Available online: https://doi.org/10.1371/journal.pone.0183714 (accessed on 1 November 2020).
  72. Lärkeryd, A.; Myrtennäs, K.; Karlsson, E.; Dwibedi, C.K.; Forsman, M.; Larsson, P.; Johansson, A.; Sjödin, A. CanSNPer: A hierarchical genotype classifier of clonal pathogens. Bioinformatics. 2014, 30, 1762–1764. Available online: http://cancerres.aacrjournals.org/cgi/doi/10.1158/0008-5472.CAN-15-3332 (accessed on 1 November 2020).
  73. Chanturia, G.; Birdsell, D.N.; Kekelidze, M.; Zhgenti, E.; Babuadze, G.; Tsertsvadze, N.; Tsanava, S.; Imnadze, P.; Beckstrom-Sternberg, S.M.; Beckstrom-Sternberg, J.; et al. Phylogeography of Francisella tularensis subspecies holarctica from the country of Georgia. BMC Microbiol. 2011, 11, 139. [Google Scholar] [CrossRef] [Green Version]
  74. Karlsson, E.; Svensson, K.; Lindgren, P.; Byström, M.; Sjödin, A.; Forsman, M.; Johansson, A. The phylogeographic pattern ofFrancisella tularensisin Sweden indicates a Scandinavian origin of Eurosiberian tularaemia. Environ. Microbiol. 2012, 15, 634–645. [Google Scholar] [CrossRef] [PubMed]
  75. Johansson, A.; Lärkeryd, A.; Widerström, M.; Mörtberg, S.; Myrtännäs, K.; Ohrman, C.; Birdsell, D.; Keim, P.; Wagner, D.M.; Forsman, M.; et al. A respiratory tularemia outbreak caused by diverse clones of Francisella tularensis. Clin. Infect. Dis. 2014, 59, 1546–1553. Available online: http://cid.oxfordjournals.org/content/early/2014/08/05/cid.ciu621.abstract (accessed on 1 November 2020). [CrossRef] [PubMed] [Green Version]
  76. Gyuranecz, M.; Szeredi, L.; Makrai, L.; Fodor, L.; Meszaros, A.R.; Szépe, B.; Füleki, M.; Erdélyi, K. Tularemia of European Brown Hare (Lepus euro-paeus). Vet. Pathol. Online 2010, 47, 958–963. Available online: http://vet.sagepub.com/content/47/5/958.abstract (accessed on 1 November 2020). [CrossRef] [PubMed]
  77. Svensson, K.; Larsson, P.; Johansson, D.; Byström, M.; Forsman, M.; Johansson, A. Evolution of Subspecies of Francisella tularensis. J. Bacteriol. 2005, 187, 3903–3908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Gerhart, J.G.; Moses, A.S.; Raghavan, R. A Francisella-like endosymbiont in the Gulf Coast tick evolved from a mammalian patho-gen. Sci. Rep. 2016, 6, 33670. Available online: http://www.ncbi.nlm.nih.gov/pubmed/27645766 (accessed on 1 November 2020). [CrossRef] [PubMed]
  79. Liu, J.N.; Yu, Z.J.; Liu, L.M.; Li, N.X.; Wang, R.R.; Zhang, C.M.; Liu, J.Z. Identification, Distribution and Population Dynamics of Francisella-like Endosymbiont in Haemaphysalis doenitzi (Acari: Ixodidae). Sci. Rep. 2016, 6, 35178. Available online: http://www.ncbi.nlm.nih.gov/pubmed/27731377%5Cnhttp://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=PMC5059625 (accessed on 1 November 2020). [CrossRef] [Green Version]
  80. Angelina, W.-F.; Zając, V.; Sawczyn, A.; Cisak, E.; Sroka, J.; Dutkiewicz, J. Occurrence of Francisella spp. in Dermacentor reticulatus and Ixodes ricinus ticks collected in eastern Poland. Ticks Tick-Borne Dis. 2015, 6, 253–257. [Google Scholar] [CrossRef]
  81. Ivanov, I.N.; Mitkova, N.; Reye, A.L.; Hübschen, J.M.; Vatcheva-Dobrevska, R.S.; Dobreva, E.G.; Kantardjiev, T.V.; Muller, C.P. Detection of New Francisella-Like Tick Endosymbionts in Hyalomma spp. and Rhipicephalus spp. (Acari: Ixodidae) from Bulgaria. Appl. Environ. Microbiol. 2011, 77, 5562–5565. [Google Scholar] [CrossRef] [Green Version]
  82. Kingry, L.; Petersen, J.M. Comparative review of Francisella tularensis and Francisella novicida. Front. Cell. Infect. Microbiol. 2014, 4, 35. [Google Scholar] [CrossRef]
  83. Whipp, M.J.; Davis, J.M.; Lum, G.; De Boer, J.; Zhou, Y.; Bearden, S.W.; Petersen, J.M.; Chu, M.C.; Hogg, G. Characterization of a novicida-like subspecies of Francisella tularensis isolated in Australia. J. Med. Microbiol. 2003, 52, 839–842. [Google Scholar] [CrossRef]
  84. Berrada, Z.L.; Telford, S.R., III. Survival of Francisella tularensis Type A in brackish-water. Arch. Microbiol. 2011, 193, 223–226. Available online: http://www.springerlink.com/content/f277542lg8860482 (accessed on 1 November 2020). [CrossRef] [PubMed] [Green Version]
  85. Wenger, J.D.; Hollis, D.G.; Weaver, R.E.; Baker, C.N.; Brown, G.R.; Brenner, D.J.; Broome, C.V. Infection Caused byFrancisella philomiragia(FormerlyYersinia philomiragia). Ann. Intern. Med. 1989, 110, 888–892. [Google Scholar] [CrossRef] [PubMed]
  86. Li, L.-H.; Luo, H.-M.; Feng, J.-H.; Ming, Y.-Z.; Zheng, M.-L.; Deng, G.-Y.; Chen, C.; Li, W.; Xiao, M.; Qu, P. Francisella salimarina sp. nov., isolated from coastal seawater. Int. J. Syst. Evol. Microbiol. 2020, 70, 1–9. [Google Scholar] [CrossRef] [PubMed]
  87. Hennebique, A.; Boisset, S.; Maurin, M. Tularemia as a waterborne disease: A review. Emerg. Microbes Infect. 2019, 8, 1027–1042. [Google Scholar] [CrossRef] [PubMed]
  88. Schrallhammer, M.; Schweikert, M.; Vallesi, A.; Verni, F.; Petroni, G. Detection of a Novel Subspecies of Francisella noatunensis as Endosymbiont of the Ciliate Euplotes raikovi. Microb. Ecol. 2010, 61, 455–464. [Google Scholar] [CrossRef]
  89. Rohmer, L.; Fong, C.; Abmayr, S.; Wasnick, M.; Freeman, T.J.L.; Radey, M.; Guina, T.; Svensson, K.; Hayden, H.S.; Jacobs, M.A.; et al. Comparison of Francisella tularensis genomes reveals evolutionary events associated with the emergence of human pathogenic strains. Genome Biol. 2007, 8, R102–R116. [Google Scholar] [CrossRef] [Green Version]
  90. Champion, M.D.; Zeng, Q.; Nix, E.B.; Nano, F.E.; Keim, P.; Kodira, C.D.; Borowsky, M.; Young, S.; Koehrsen, M.; Engels, R.; et al. Comparative Genomic Characterization of Francisella tularensis Strains Belonging to Low and High Virulence Subspecies. PLoS Pathog. 2009, 5, e1000459. [Google Scholar] [CrossRef]
  91. Skaldin, M.; Tuittila, M.; Zavialov, A.V.; Zavialov, A.V. Secreted Bacterial Adenosine Deaminase Is an Evolutionary Precursor of Adenosine Deaminase Growth Factor. Mol. Biol 2018, 35, 2851–2861. Available online: https://academic.oup.com/mbe/article/35/12/2851/5133304 (accessed on 1 November 2020). [CrossRef]
  92. Gallagher, L.A.; McKevitt, M.; Ramage, E.R.; Manoil, C. Genetic Dissection of the Francisella novicida Restriction Barrier. J. Bacteriol. 2008, 190, 7830–7837. [Google Scholar] [CrossRef] [Green Version]
  93. Dorofe’ev, K.A. Classification of the causative agent of tularemia. Symp. Res. Work Inst. Epidemiol. Microbiol. Chita 1947, 170–180. [Google Scholar]
  94. McCoy, G.G.W.; Chapin, C.W.C. Further observations on a plague-like disease of rodents with a preliminary note on the causative agent, Bacterium tularense. J Infect Dis. 1912, 10, 61–72. Available online: http://www.jstor.org/stable/10.2307/30071893 (accessed on 1 November 2020). [CrossRef]
  95. Olsufiev, N.G.; Meshcheryakova, I.S. Subspecific Taxonomy of Francisella tularensis McCoy and Chapin 1912. Int. J. Syst. Bacteriol. 1983, 33, 872–874. Available online: http://ijs.sgmjournals.org/content/33/4/872.abstract (accessed on 1 November 2020). [CrossRef] [Green Version]
  96. Svensson, K.; Sjödin, A.; Byström, M.; Granberg, M.; Brittnacher, M.J.; Rohmer, L.; Jacobs, M.A.; Sims-Day, E.H.; Levy, R.; Zhou, Y.; et al. Genome sequence of Francisella tularensis subspecies holarctica strain FSC200, isolated from a child with tularemia. J. Bacteriol. 2012, 194, 6965–6966. Available online: http://jb.asm.org/content/194/24/6965.abstract (accessed on 1 November 2020). [CrossRef] [PubMed] [Green Version]
  97. Aikimbaev, M. Taxonomy of Genus Francisella. Rep. Acad. Sci. Kaz. SSR Ser. Biol. 1966, 5, 42–44. [Google Scholar]
  98. Olsufjev, N.G.; Meshcheryakova, I.S. Infraspecific taxonomy of tularemia agent Francisella tularensis McCoy et Chapin. J. Hyg. Epidemiol. Microbiol. Immunol. 1982, 26, 291–299. Available online: http://www.ncbi.nlm.nih.gov/pubmed/7142690 (accessed on 1 November 2020).
  99. Jr, E.C.S.; Weiss, E.; Suitor, E. Isolation of a Rickettsialike Microorganism (Wolbachia persica, n. sp.) from Argas persicus (Oken). J. Infect. Dis. 1961, 108, 95–106. Available online: http://jid.oxfordjournals.org/content/108/1/95.full.pdf (accessed on 1 November 2020).
  100. Forsman, M.; Sandström, G.; Sjöstedt, A. Analysis of 16S ribosomal DNA sequences of Francisella strains and utilization for determination of the phylogeny of the genus and for identification of strains by PCR. Int. J. Syst. Bacteriol. 1994, 44, 38–46. Available online: http://ijs.sgmjournals.org/content/44/1/38.abstract (accessed on 1 November 2020). [CrossRef]
  101. Noda, H.; Munderloh, U.G.; Kurtti, T.J. Endosymbionts of ticks and their relationship to Wolbachia spp. and tick-borne pathogens of humans and animals. Appl. Environ. Microbiol. 1997, 63, 3926–3932. Available online: http://aem.asm.org/cgi/content/abstract/63/10/3926 (accessed on 1 November 2020). [CrossRef] [Green Version]
  102. Kugeler, K.J.; Mead, P.S.; McGowan, K.L.; Burnham, J.M.; Hogarty, M.D.; Ruchelli, E.; Pollard, K.; Husband, B.; Conley, C.; Rivera, T.; et al. Isolation and characterization of a novel Francisella sp. from human cerebrospinal fluid and blood. J. Clin. Microbiol. 2008, 46, 2428–2431. [Google Scholar] [CrossRef] [Green Version]
  103. Dietrich, E.A.; Kingry, L.C.; Kugeler, K.J.; Levy, C.; Yaglom, H.; Young, J.W.; Mead, P.S.; Petersen, J.M. Francisella opportunistica sp. Nov., isolated from human blood and cerebrospinal fluid. Int. J. Syst. Evol. Microbiol. 2020, 70, 1145–1151. [Google Scholar] [CrossRef]
  104. Jensen, W.I.; Owen, C.R.; Jellison, W.L. Yersinia philomiragia sp. n., a new member of the Pasteurella group of bacteria, naturally pathogenic for the muskrat (Ondatra zibethica). J. Bacteriol. 1969, 100, 1237–1241. Available online: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=250302&tool=pmcentrez&rendertype=abstract (accessed on 1 November 2020). [CrossRef] [PubMed] [Green Version]
  105. Ottem, K.F.; Nylund, A.; Karlsbakk, E.; Friis-Møller, A.; Krossøy, B.; Knappskog, D. New species in the genus Francisella (Gammaproteobacteria; Francisellaceae); Francisella piscicida sp. nov. isolated from cod (Gadus morhua). Arch Microbiol. 2007, 188, 547–550. Available online: http://www.springerlink.com/content/q464n38w86k248l4 (accessed on 1 November 2020). [CrossRef] [PubMed]
  106. Svensson, D.; Öhrman, C.; Bäckman, S.; Karlsson, E.; Nilsson, E.; Byström, M.; Lärkeryd, A.; Myrtennäs, K.; Stenberg, P.; Qu, P.H.; et al. Complete Genome Sequence of Francisella guangzhouensis Strain 08HL01032T, Isolated from Air-Conditioning Systems in China. Genome Announc. 2015, 3, e00024-15. Available online: http://genomea.asm.org/lookup/doi/10.1128/genomeA.00024-15 (accessed on 1 November 2020). [CrossRef] [PubMed] [Green Version]
Microorganisms 09 00146 g001 550
Figure 1. Whole genome maximum likelihood phylogeny of the complete Francisella near neighbor (FNN) dataset calculated with RAxML. Genus Francisella is divided into four separate clades indicated with different colors and numbered 1–4. The four most proximate genera to Francisella are also presented (the complete phylogenetic tree of Francisellaceae is presented in Figure A1).
Figure 1. Whole genome maximum likelihood phylogeny of the complete Francisella near neighbor (FNN) dataset calculated with RAxML. Genus Francisella is divided into four separate clades indicated with different colors and numbered 1–4. The four most proximate genera to Francisella are also presented (the complete phylogenetic tree of Francisellaceae is presented in Figure A1).
Microorganisms 09 00146 g001
Microorganisms 09 00146 g002 550
Figure 2. Clustered heatmap visualizing ANI values between Francisella, Allofrancisella, Parafrancisella, and Pseudofrancisella species.
Figure 2. Clustered heatmap visualizing ANI values between Francisella, Allofrancisella, Parafrancisella, and Pseudofrancisella species.
Microorganisms 09 00146 g002
Microorganisms 09 00146 g003 550
Figure 3. Phylogenetic tree of F. tularensis visualizing subspecies and major genetic subgroups.
Figure 3. Phylogenetic tree of F. tularensis visualizing subspecies and major genetic subgroups.
Microorganisms 09 00146 g003
Microorganisms 09 00146 g004 550
Figure 4. Visualization of genomic regions that are differentially present/absent across the three subspecies of F. tularensis and/or major phylogenetic groups within the subspecies. A phylogeny is presented in the middle of the figure, and the multiple circles around the phylogeny correspond to different genomic regions. The color of each cell in the circles corresponds to the BSR value for each genomic component in each examined whole genome sequence.
Figure 4. Visualization of genomic regions that are differentially present/absent across the three subspecies of F. tularensis and/or major phylogenetic groups within the subspecies. A phylogeny is presented in the middle of the figure, and the multiple circles around the phylogeny correspond to different genomic regions. The color of each cell in the circles corresponds to the BSR value for each genomic component in each examined whole genome sequence.
Microorganisms 09 00146 g004
Table
Table 1. Primers and TagMan-MGB probes for real-time assays targeting F. tularensis.
Table 1. Primers and TagMan-MGB probes for real-time assays targeting F. tularensis.
AssayRegion in Genome of
SCHU S4 a		Sequences of Forward and Reverse Primers b and Probes cPrimer (µM) dMGB eProbe (µM) fTm °C g
Ft-sp. FTT0376377718–377824F CCATATCACTGGCTTTGCTAGACTAGT
R TGTTGGCAAAAGCTAAAGAGTCTAAA
FAM-AAATTATAAAACCAAACCCAGACCTTCAAACCACA		0.9no0.260
Ft-sp.3Pan8967–9049F TTTACACCCGTCTCCGTTAGT
R CTCTTAAGGATGCAATTTGGGATT
FAM-AAGAGGCAAAGCTGGAATTACACTCTCTC		0.9no0.260
Ft-sp. FTS_07721187646–1187807F CAAGGTAAAGAAATTAAGAGTAGTAAAGTTGAATTC
R ATTTAATCTAGTATTATCAATTGGGTAAAAAGGTA
FAM-GATGTGGCAACAACTGAAAT		0.9yes0.263
a Position coordinates of assay target within reference SchuS4 whole genome sequence (AJ749949.2); b F: forward primer, R: reverse primer; c 6-carboxyfluorescein (FAM) fluorophore and a 3= black hole quencher-1 (BHQ-1) quencher; d Primer concentration (μM) per reaction; e Minor Groove Binding moiety on probe; f Probe concentration (μM) per reaction; g Real-time PCR annealing temperature.
Table
Table 2. Overview of genome distribution within clades and species, including detailed information about the number of unique regions for each species and pan/core genome size.
Table 2. Overview of genome distribution within clades and species, including detailed information about the number of unique regions for each species and pan/core genome size.
GenusCladeSpecies#Genomes#genomes#Unique RegionsPan-Genome Size (#Genes)Core-Genome Size (#Genes)
Francisellaclade 1F. tularensis374327624011541
F. novicida37229491499
F. sp. TX07-66081250N/AN/A
F. hispaniensis312120151599
FLE and F. persica317022931006
F. opportunistica346417511713
clade 2F. philomiragia113471029921611
F. noatunensis237419441766
F. sp. GA01-27941180N/AN/A
F. orientalis303622872025
F. salimarina719123711786
clade 3F. endociliophora111340N/AN/A
clade 4F. halioticida411190N/AN/A
F. uliginis1736N/AN/A
F. sp. LA11-24451587N/AN/A
F. sp. SYW 911179N/AN/A
Parafrancisella F. adeliensis33187419361925
Allofrancisella A. inopinata71377N/AN/A
A. frigidaquae221114791356
A. guangzhouensis471515701455
Pseudofrancisella P. aestuarii21375N/AN/A
F. frigiditurris1406N/AN/A
Table
Table 3. Detailed information regarding the six genomic regions unique to F. tularensis based upon in silico analyses.
Table 3. Detailed information regarding the six genomic regions unique to F. tularensis based upon in silico analyses.
RegionLVS * CoordinatesSCHU S4 ** CoordinatesMarker LengthAnnotationNotes
11490829–14911611187496–1187828333hsdS, type I restriction modification domain proteinpartial locus
284841–85406378437–377872567FTT_0376c, hypothetical membrane proteinpartial locus
31673808–1673858547774–54782451FTT_0525, conserved hypothetical protein
443283–43397951193–951307114add1, adenosine deaminasepartial locus
543766–43874950716–950824109intergenic, region between bioA and add1
6135869–136156765604–765891288FTT_0742, hypothetical proteinpartial locus
* CP009694.1, ** GCA_000008985.1.
Table
Table 4. Number of genomic regions that are unique to or have been lost across the three subspecies of F. tularensis and/or major phylogenetic groups within the subspecies.
Table 4. Number of genomic regions that are unique to or have been lost across the three subspecies of F. tularensis and/or major phylogenetic groups within the subspecies.
SpeciesSubspeciesGenetic Group#Unique Regions#Regions LostNumber of F. tularensis Subspecies Region is Present in
F. tularensistularensisA 7013
A13902
A23102
holarcticaB (excluding japonica)25000
japonica512
B including japonica25572
mediasiaticaM732
Table
Table 5. Assay performance for isolates included in the validation panel. Number of isolates indicates the number of strains per species/subspecies included in the validation panel. Numbers in the other columns indicate the number of positive results (i.e., purported F. tularensis-positive results) for each of the three examined PCR assays.
Table 5. Assay performance for isolates included in the validation panel. Number of isolates indicates the number of strains per species/subspecies included in the validation panel. Numbers in the other columns indicate the number of positive results (i.e., purported F. tularensis-positive results) for each of the three examined PCR assays.
#IsolatesFt-sp.FTT0376Ft-sp.3Pan Ft-sp.FTS_0772
Francisella tularensis subspecies tularensis25252525
Francisella tularensis subspecies holarctica13131313
Francisella tularensis subspecies mediasiatica2222
Francisella novicida7
Francisella hispaniensis2
Francisella persica1
Francisella FLE1
Francisella opportunistica1
Francisella philomiragia23 2
Francisella halioticida1
Francisella noatunensis1
Francisella orientalis1
Francisella salimarina1
Francisella endociliophora1
Parafrancisella adeliensis1
Allofrancisella inopinata1
Allofrancisella frigidaquae1
Allofrancisella guangzhouensis1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Öhrman, C.; Sahl, J.W.; Sjödin, A.; Uneklint, I.; Ballard, R.; Karlsson, L.; McDonough, R.F.; Sundell, D.; Soria, K.; Bäckman, S.; et al. Reorganized Genomic Taxonomy of Francisellaceae Enables Design of Robust Environmental PCR Assays for Detection of Francisella tularensis. Microorganisms 2021, 9, 146. https://doi.org/10.3390/microorganisms9010146

AMA Style

Öhrman C, Sahl JW, Sjödin A, Uneklint I, Ballard R, Karlsson L, McDonough RF, Sundell D, Soria K, Bäckman S, et al. Reorganized Genomic Taxonomy of Francisellaceae Enables Design of Robust Environmental PCR Assays for Detection of Francisella tularensis. Microorganisms. 2021; 9(1):146. https://doi.org/10.3390/microorganisms9010146

Chicago/Turabian Style

Öhrman, Caroline, Jason W. Sahl, Andreas Sjödin, Ingrid Uneklint, Rebecca Ballard, Linda Karlsson, Ryelan F. McDonough, David Sundell, Kathleen Soria, Stina Bäckman, and et al. 2021. "Reorganized Genomic Taxonomy of Francisellaceae Enables Design of Robust Environmental PCR Assays for Detection of Francisella tularensis" Microorganisms 9, no. 1: 146. https://doi.org/10.3390/microorganisms9010146

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop