Search

Tularemia is a widely spread zoonotic disease in the northern hemisphere, caused by the bacterium Francisella tularensis. In humans, tularemia is an acute febrile illness with incidence peaks in late summer to early autumn of outbreak years, but there is no early warning system in place that can reduce the impact of disease by providing timely risk information. In this study, we revisit previously unpublished data on F. tularensis in water, sediment, soil, and small mammals from 1984 in northern Sweden. In addition, we used human case data from the national surveillance system for tularemia in the same year. In the environmental and small mammal material, bank vole (Myodes glareolus) samples from urine and bladder were the only samples that tested positive for F. tularensis. The prevalence of F. tularensis among trapped bank voles was 13.5%, although all six bank voles that were retrieved from owl nest boxes in early May tested positive. Forty-two human tularemia cases were reported from August to December in 1984. Based on these results, we encourage investigating the potential role of tularemia-infected bank voles retrieved from owl nest boxes in spring as an early warning for outbreaks of tularemia among humans in summer and autumn of the same year.

Tularemia, a disease caused by the environmental bacterium Francisella tularensis, is characterized by acute febrile illness. F. tularensis is highly infectious: as few as 10 organisms can cause human disease. Tularemia is not known to be spread from person to person., Francisella tularensis, the causative agent of the zoonotic disease tularemia, can cause seasonal outbreaks of acute febrile illness in humans with disease peaks in late summer to autumn. Interestingly, its mechanisms for environmental persistence between outbreaks are poorly understood. One hypothesis is that F. tularensis forms biofilms in aquatic environments. We utilized two fully virulent wild-type strains: FSC200 (Francisella tularensis subsp. holarctica) and Schu S4 (Francisella tularensis subsp. tularensis) and three control strains, the attenuated live vaccine strain (LVS; F. tularensis subsp. holarctica), a Schu S4 ΔwbtI mutant that is documented to form biofilms, and the low-virulence strain U112 of the closely related species Francisella novicida. Strains were incubated in saline solution (0.9% NaCl) microcosms for 24 weeks at both 4°C and 20°C, whereupon viability and biofilm formation were measured. These temperatures were selected to approximate winter and summer temperatures of fresh water in Scandinavia, respectively. U112 and Schu S4 ΔwbtI formed biofilms, but F. tularensis strains FSC200 and Schu S4 and the LVS did not. All strains exhibited prolonged viability at 4°C compared to 20°C. U112 and FSC200 displayed remarkable long-term persistence at 4°C, with only 1- and 2-fold log reductions, respectively, of viable cells after 24 weeks. Schu S4 exhibited lower survival, yielding no viable cells by week 20. At 24 weeks, cells from FSC200, but not from Schu S4, were still fully virulent in mice. Taken together, these results demonstrate biofilm-independent, long-term survival of pathogenic F. tularensis subsp. holarctica in conditions that mimic overwinter survival in aquatic environments. IMPORTANCE Tularemia, a disease caused by the environmental bacterium Francisella tularensis, is characterized by acute febrile illness. F. tularensis is highly infectious: as few as 10 organisms can cause human disease. Tularemia is not known to be spread from person to person. Rather, all human infections are independently acquired from the environment via the bite of blood-feeding arthropods, ingestion of infected food or water, or inhalation of aerosolized bacteria. Despite the environmental origins of human disease events, the ecological factors governing the long-term persistence of F. tularensis in nature between seasonal human outbreaks are poorly understood. The significance of our research is in identifying conditions that promote long-term survival of fully virulent F. tularensis outside a mammalian host or insect vector. These conditions are similar to those found in natural aquatic environments in winter and provide important new insights on how F. tularensis may persist long-term in the environment.

Here, we present the 2,139,666-bp circular chromosome of Francisella sp. strain LA11-2445 (FDC406), a proposed novel species of Francisella that was isolated from a human cutaneous lesion and is related to Francisella species from marine environments.

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.

Francisella tularensis is the causative agent of the zoonotic disease tularemia. In Germany, most human infections are caused by contact with infected hares. The aim of this study was to characterize Francisella tularensis subsp. holarctica strains isolated from hares in Germany and to develop bioinformatics tools to analyze their genetic relatedness. In total, 257 German isolates&mdash, obtained mainly from hares (n = 233), other vertebrate animals, and ticks, but also from humans (n = 3)&mdash, were analyzed within this study. Publically available sequence data from 49 isolates were used to put our isolates into an epidemiological context and to compare isolates from natural foci and humans. Whole-genome sequences were analyzed using core-genome Multi-Locus-Sequence-Typing, canonical Single Nucleotide Polymorphism (SNP) typing and whole-genome SNP typing. An overall conformity of genotype clustering between the typing methods was found, albeit with a lower resolution for canonical single SNP typing. The subclade distribution, both on local and national levels, among strains from humans and hares was similar, suggesting circulation of the same genotypes both in animals and humans. Whilst close to identical isolates of the same subclade were found distributed over large areas, small geographical foci often harbored members of different subclades. In conclusion, although genomic high-resolution typing was shown to be robust, reproducible and allowed the identification of highly closely related strains, genetic profiling alone is not always conclusive for epidemiological linkage of F. tularensis strains.

Near neighbors to the causative agent of tularemia, Francisella tularensis, isolated from diverse sources, have been reported in recent years. In this announcement, we present the complete sequence of the circular chromosome of one of the closest neighboring genera of Francisella, the type strains of Allofrancisella inopinata and Allofrancisella frigidaquae., Near neighbors to the causative agent of tularemia, Francisella tularensis, isolated from diverse sources, have been reported in recent years. In this announcement, we present the complete sequences of circular chromosomes of one of the closest neighboring genera of Francisella (i.e., the type strains of Allofrancisella inopinata and Allofrancisella frigidaquae).

Here, we report the complete genome sequence of Francisella tularensis subsp. holarctica strain A271_1, isolated from a Eurasian beaver (Castor fiber) in 2012 in the Berlin/Brandenburg region, Germany.

Here, we announce the complete genome sequence of the Francisella halioticida type strain DSM 23729 (FSC1005), isolated from a diseased cultured giant abalone in Japan in 2005. The genome is composed of a 2,197,430-bp-long circular chromosome, with a G+C content of 31.2%.

Water in a full-scale drinking water treatment plant was irradiated with ultraviolet (UV) doses of 250, 400, and 600 J/m2, and the effect on bacterial communities investigated using 16s rRNA gene amplicon sequencing, heterotrophic plate counts (HPCs), coliform, and Escherichia coli counts. The bacteria in the irradiated water were also analyzed following storage for 6 days at 7 °C, to approximate the conditions in the distribution system. The log10 reduction of HPCs at 400 J/m2 was 0.43 ± 0.12. Phylogenetic examination, including DESeq2 analysis, showed that Actinobacteria was more resistant to UV irradiation, whereas Bacteroidetes was sensitive to UV. Phylum Proteobacteria contained monophyletic groups that were either sensitive or resistant to UV exposure. The amplicon sequence variants (ASVs) resistant to UV irradiation had a greater average GC content than the ASVs sensitive to UV, at 55% ± 1.7 (n = 19) and 49% ± 2.5 (n = 16), respectively. Families Chitinophagaceae, Pelagibacteraceae, Holophagaceae, Methylophilaceae, and Cytophagaceae decreased linearly in relative abundance, with increasing UV dose (P Chitinophagaceae, Comamonadaceae, and Flavobacteriaceae families decreased in relative abundance, whereas ACK-M1, Mycobacteriaceae, and Nitrosomonadaceae were increasing in relative abundance. This suggests that the impact of UV irradiation cannot only be considered directly after application but that this treatment step likely continues to influence microbial dynamics throughout the distribution system.

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.

Microorganisms are essential constituents of ecosystems. To improve our understanding of how various factors shape microbial diversity and composition in nature it is important to study how microorganisms vary in space and time. Factors shaping microbial communities in ground level air have been surveyed in a limited number of studies, indicating that geographic location, season and local climate influence the microbial communities. However, few have surveyed more than one location, at high latitude or continuously over more than a year. We surveyed the airborne microbial communities over two full consecutive years in Kiruna, in the Arctic boreal zone, and Ljungbyhed, in the Southern nemoral zone of Sweden, by using a unique collection of archived air filters. We mapped both geographic and seasonal differences in bacterial and fungal communities and evaluated environmental factors that may contribute to these differences and found that location, season and weather influence the airborne communities. Location had stronger influence on the bacterial community composition compared to season, while location and season had equal influence on the fungal community composition. However, the airborne bacterial and fungal diversity showed overall the same trend over the seasons, regardless of location, with a peak during the warmer parts of the year, except for the fungal seasonal trend in Ljungbyhed, which fluctuated more within season. Interestingly, the diversity and evenness of the airborne communities were generally lower in Ljungbyhed. In addition, both bacterial and fungal communities varied significantly within and between locations, where orders like Rhizobiales, Rhodospirillales and Agaricales dominated in Kiruna, whereas Bacillales, Clostridiales and Sordariales dominated in Ljungbyhed. These differences are a likely reflection of the landscape surrounding the sampling sites where the landscape in Ljungbyhed is more homogenous and predominantly characterized by artificial and agricultural surroundings. Our results further indicate that local landscape, as well as seasonal variation, shapes microbial communities in air.

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.

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. Keywords: Francisella; genome; sequencing; Illumina; outbreak; phylogeography; subtyping; canSNP; assay; zoonosis (Published: 27 October 2016) Citation: Infection Ecology and Epidemiology 2016, 6: 32838 - http://dx.doi.org/10.3402/iee.v6.32838

SUMMARYIn Germany tularemia is a re-emerging zoonotic disease. Therefore, we investigated wild animals and environmental water samples for the presence and phylogenetic diversity of Francisella tularensis in the poorly studied Berlin/Brandenburg region. The phylogenomic analysis of three isolates from wild animals revealed three new subclades within the phylogenetic tree of F. tularensis [B.71 from a raccoon dog (Nyctereutes procyonoides); B.74 from a red fox (Vulpes vulpes), and B.75 from a Eurasian beaver (Castor fiber albicus)]. The results from histological, PCR, and genomic investigations on the dead beaver showed that the animal suffered from a systemic infection. Indications were found that the bacteria were released from the beaver carcass into the surrounding environment. We demonstrated unexpectedly high and novel phylogenetic diversity of F. tularensis in Germany and the fact that the bacteria persist in the environment for at least one climatic season. These findings support a broader host species diversity than previously known regarding Germany. Our data further support the assumption derived from previous serological studies of an underestimated frequency of occurrence of the pathogen in the environment and in wild animals. F. tularensis was isolated from animal species not previously reported as natural hosts in Germany.

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

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.

Francisella tularensis, a highly virulent bacteria that causes the zoonotic disease tularemia, is considered a potential agent of biological warfare and bioterrorism. Although the host range for several species within the Francisella is known, little is known about the natural reservoirs of various Francisella species. The lack of knowledge regarding the environmental fates of these pathogens greatly reduces the possibilities for microbial risk assessments. The greater wax moth (Galleria mellonella) is an insect of the order Lepidoptera that has been used as an alternative model to study microbial infection during recent years. The aim of this study was to evaluate G. mellonella as a model system for studies of human pathogenic and closely related opportunistic and non-pathogenic strains within the Francisella genus. The employed G. mellonella larvae model demonstrated differences in lethality between human pathogenic and human non-pathogenic or opportunistic Francisella species. The F. novicida, F. hispaniensis and F. philomiragia strains were significantly more virulent in the G. mellonella model than the strains of human pathogens F. t. holarctica and F. t. tularensis. Our data show that G. mellonella is a possible in vivo model of insect immunity for studies of both opportunistic and virulent lineages of Francisella spp., that produces inverse results regarding lethality in G. mellonella and incapacitating disease in humans. The results provide insight into the potential host specificity of F. tularensis and closely related members of the same genus, thus increasing our present understanding of Francisella spp. ecology.

Real-time polymerase chain reaction (qPCR) is the cornerstone of DNA analysis, enabling detection and quantification of minute nucleic acid amounts. However, PCR-based analysis is limited, in part, by the presence of inhibitors in the samples. PCR inhibition has been viewed solely as failure to efficiently generate amplicons, that is, amplification inhibition. Humic substances (HS) are well-known inhibitors of PCR amplification. Here we show that HS from environmental samples, specifically humic acid (HA), are very potent detection inhibitors, that is, quench the fluorescence signal of double-stranded DNA (dsDNA) binding dyes. HA quenched the fluorescence of the commonly used qPCR dyes EvaGreen, ResoLight, SYBR Green I, and SYTO 82, generating lowered amplification plots, although amplicon production was unaffected. For EvaGreen, 500ng of HA quenched nearly all fluorescence, whereas 1000ng of HA completely inhibited amplification when applying Immolase DNA polymerase with bovine serum albumin (BSA). Fluorescence spectroscopy measurements showed that HA quenching was either static or collisional and indicated that HA bound directly to the dye. Fulvic acid did not act as a qPCR detection inhibitor but inhibited amplification similarly to HA. Hydrolysis probe fluorescence was not quenched by HA. Detection inhibition is an overlooked phenomenon that needs to be considered to allow for development of optimal qPCR assays.

Here, we report a high-quality draft genome sequence of Francisella tularensis subsp. holarctica strain 08T0073, isolated from the cadaver of a wild European hare ( Lepus europaeus ) found near Helmstedt, Lower Saxony, Germany, in 2007. In Germany, infected hares are a major source of tularemia in humans.

In Sweden, mosquitoes are considered the major vectors of the bacterium Francisella tularensis subsp. holarctica, which causes tularaemia. The aim of this study was to investigate whether mosquitoes acquire the bacterium as aquatic larvae and transmit the disease as adults. Mosquitoes sampled in a Swedish area where tularaemia is endemic (A-rebro) were positive for the presence of F. tularensis deoxyribonucleic acid throughout the summer. Presence of the clinically relevant F. tularensis subsp. holarctica was confirmed in 11 out of the 14 mosquito species sampled. Experiments performed using laboratory-reared Aedes aegypti confirmed that F. tularensis subsp. holarctica was transstadially maintained from orally infected larvae to adult mosquitoes and that 25 % of the adults exposed as larvae were positive for the presence of F. tularensis-specific sequences for at least 2 weeks. In addition, we found that F. tularensis subsp. holarctica was transmitted to 58 % of the adult mosquitoes feeding on diseased mice. In a small-scale in vivo transmission experiment with F. tularensis subsp. holarctica-positive adult mosquitoes and susceptible mice, none of the animals developed tularaemia. However, we confirmed that there was transmission of the bacterium to blood vials by mosquitoes that had been exposed to the bacterium in the larval stage. Taken together, these results provide evidence that mosquitoes play a role in disease transmission in part of Sweden where tularaemia recurs.

Francisella tularensis is a facultative intracellular bacterium in the class Gammaproteobacteria. This strain is of interest because it is the etiologic agent of tularemia and a highly virulent category A biothreat agent. Here we describe the draft genome sequence and annotation of Francisella tularensis subsp. holarctica BD11-00177, isolated from the first case of indigenous tularemia detected in The Netherlands since 1953. Whole genome DNA sequence analysis assigned this isolate to the genomic group B.FTNF002–00, which previously has been exclusively reported from Spain, France, Italy, Switzerland and Germany. Automatic annotation of the 1,813,372 bp draft genome revealed 2,103 protein-coding and 46 RNA genes.

Previous studies of the causative agent of tularaemia, Francisella tularensis have identified phylogeographic patterns suggestive of environmental maintenance reservoirs. To investigate the phyloge ...

Specific identification and differentiation of the 2 subspecies of the fish pathogen Francisella noatunensis, namely, F. noatunensis subsp. noatunensis and F. noatunensis subsp. orientalis, remains a major diagnostic challenge. Following whole-genome sequencing and analysis of representatives of all major subclades of the genus Francisella, specific genomic regions were identified for each of the subspecies of this fish pathogen. Two specific real-time quantitative PCR assays, directed at hypothetical genes within these regions were developed. Specificity was confirmed by lack of signal and cross-reactivity with the closest relative, F. philomiragia, and other common bacterial fish pathogens. Both assays, used either as monoplex or multiplex, have a limit of detection of 10 genome equivalents. The quantitative sensitivity of the assays was not affected by the presence of kidney tissues or DNA from Atlantic cod Gadus morhua or tilapia Oreochromis sp.

In Sweden, human cases of tularemia caused by Francisella tularensis holarctica are assumed to be transmitted by mosquitoes, but how mosquito vectors acquire and transmit the bacterium is not clear. To determine how transmission of this bacterium occurs, mosquito larvae were collected in an area where tularemia is endemic, brought to the laboratory, and reared to adults in their original pond water. Screening of adult mosquitoes by real-time PCR demonstrated F. tularensis lpnA sequences in 14 of the 48 mosquito pools tested; lpnA sequences were demonstrated in 6 of 9 identifi ed mosquito species. Further analysis confi rmed the presence of F. tularensis holarctica–specifi c 30-bp deletion region sequences (FtM19inDel) in water from breeding containers and in 3 mosquito species (Aedes sticticus, Ae. vexans, and Ae. punctor) known to take blood from humans. Our results suggest that the mosquitoes that transmit F. tularensis holarctica during tularemia outbreaks acquire the bacterium already as larvae.

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

The natural reservoir of Francisella tularensis, the causative agent of tularaemia, is yet to be identified. We investigated the possibility that Francisella persists in natural aquatic ecosystems between outbreaks. It was hypothesized that nutrient-rich environments, with strong protozoan predation, favour the occurrence of the tularaemia bacterium. To investigate the differences in adaptation to aquatic environments of the species and subspecies of Francisella, we screened 23 strains for their ability to survive grazing by the ciliate Tetrahymena pyriformis. All the Francisella strains tested were consumed at a low rate, although significant differences between subspecies were found. The survival and virulence of gfp-labelled F. tularensis ssp. holarctica were then studied in a microcosm experiment using natural lake water, with varying food web complexities and nutrient availabilities. High nutrient conditions in combination with high abundances of nanoflagellates were found to favour F. tularensis ssp. holarctica. The bacterium was observed both free-living and within the cells of a nanoflagellate. Francisella tularensis entered a viable but nonculturable state during the microcosm experiment. When studied over a longer period of time, F. tularensis ssp. holarctica survived in the lake water, but loss of virulence was not prevented by either high nutrient availability or the presence of predators.

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. 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-10-C-00139;] pursuant to the agreement between the Kingdom of Sweden and the US government on Cooperation in Science and Technology for Homeland Security Matters. The authors R.E. and P.A. were supported by the Ministerio de Economía y Competitividad, Spain [grant number CGL2015-66962-C2-2-R]. We thank Johanna Thelaus, Per Stenberg, Andreas Sjödin, and Petter Lindgren for comments that improved the manuscript. We are grateful to the researchers, physicians, and veterinarians who contributed F. tularensis samples for this study.

By combining analysis of indel markers with multiple-locus variable-number tandem repeat analysis, individual strains were identified., To develop effective and accurate typing of strains of Francisella tularensis, a potent human pathogen and a putative bioterrorist agent, we combined analysis of insertion-deletion (indel) markers with multiple-locus variable-number tandem repeat analysis (MLVA). From 5 representative F. tularensis genome sequences, 38 indel markers with canonical properties, i.e., capable of sorting strains into major genetic groups, were selected. To avoid markers with a propensity for homoplasy, we used only those indels with 2 allelic variants and devoid of substantial sequence repeats. MLVA included sequences with much diversity in copy number of tandem repeats. The combined procedure allowed subspecies division, delineation of clades A.I and A.II of subspecies tularensis, differentiation of Japanese strains from other strains of subspecies holarctica, and high-resolution strain typing. The procedure uses limited amounts of killed bacterial preparations and, because only 1 single analytic method is needed, is time- and cost-effective.

Top-down control of lower trophic levels, e.g., bacteria, has been suggested to increase along aquatic productivity gradients. The response by the bacterial community may be to become more predation resistant in highly productive environments. To test this hypothesis, samples were taken from 20 aquatic systems along a productivity gradient (dissolved organic carbon from 7 to 71 mg/L), during late summer. The results showed that the biomass of bacteria, phytoplankton, and ciliates increased along the gradient (r2 = 0.532, 0.426, and 0.758, P < 0.01, respectively). However, the organisms did not increase equally, and the ratio of protozoan to bacterial biomass showed a 100-fold increase along the gradient. Ciliates dominated the protozoan biomass in the more nutrient-rich waters. The edibility of colony-forming bacteria was tested using a ciliate predator, Tetrahymena pyriformis. Bacterial edibility was found to decrease with increases in nutrient richness and ciliate biomass in the aquatic systems (r2 = 0.358, P < 0.01; r2 = 0.242, P < 0.05, respectively). Quantile regression analysis indicated that the selection pressures on edible bacteria were increasing along the productivity gradient. Thus, inedible forms of bacteria were selected for in the transition from oligotrophic to eutrophic conditions. Isolated bacteria were distributed among the alpha-, beta-, and gamma- Proteobacteria and the Actinobacteria and Firmicutes taxa. We conclude that bacterial predation resistance increases in nutrient-rich waters with high protozoan predation.

Tularemia is a plague-like infection caused by Francisella (F.) tularensis classified as a biological warfare agent. F. tularensis subsp. tularensis is the most virulent subspecies demanding rapid diagnosis. Typing systems for this fastidious bacterium to the subspecies level are laborious and time consuming. Therefore, the aim of this study was to develop a real-time PCR for the rapid and specific identification of F. tularensis subsp. tularensis. The specificity of the assay was determined using a comprehensive panel of Francisella strains, clinically relevant bacteria, and DNA preparations of potential hosts. F. tularensis subsp. tularensis was specifically detected but no other organisms. The range of linearity was determined to be 100 fg to 10 ng, the lower limit of detection was 25 fg of DNA (13 genome equivalents). An internal amplification control PCR system targeting lambda phage DNA was included. Neither the internal amplification control nor host DNA influenced the cycle threshold values obtained for F. tularensis subsp. tularensis. In conclusion, we have developed a highly sensitive and specific assay that can be integrated into real-time PCR-based identification procedures for biological agents. This is a major diagnostic improvement, as all other methods for the specific identification of F. tularensis subsp. tularensis are more time consuming.

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 identifi ed 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 fi ndings 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 defi ne 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.

Analysis of unidirectional genomic deletion events and single nucleotide variations suggested that the four subspecies of Francisella tularensis have evolved by vertical descent. The analysis indicated an evolutionary scenario where the highly virulent F. tularensis subsp. tularensis (type A) appeared before the less virulent F. tularensis subsp. holarctica (type B). Compared to their virulent progenitors, attenuated strains of F. tularensis exhibited specific unidirectional gene losses.

Francisella tularensis live vaccine strain infection of mice has been established as an experimental model of tularemia that is suitable for studies of immune mechanisms against the intracellular pathogen. In this study, the model was used to explore immunogenic repertoire of F. tularensis with the aim of identifying new molecules able to activate the host immune system, potential bacterial markers with vaccine, and diagnostic applications. Immunoproteomic approach based on the combination of two-dimensional gel electrophoresis, immunoblotting, and mass spectrometry was applied. Globally, 36 different proteins were identified, which strongly reacted with sera from experimentally infected mice, including several putative virulence markers of intracellular pathogens as nucleoside diphosphate kinase, isocitrate dehydrogenase, RNA-binding protein Hfq, and molecular chaperone ClpB. Of them, 27 proteins are described for the first time as immunorelevant Francisella proteins. When comparing murine immunoproteome of F. tularensis with our previous data from human patients, 25 of the total of 50 identified murine sera immunoreactive spots were recognized by human sera collected from patients suffering from tularemia, as well. Immune sera from two Lps gene congenic strains of mice, C3H/HeN (Lpsn) and C3H/HeJ (Lpsd), represented murine immunoproteome in this study. The spectrum of immunoreactive spots detected by two-dimensional immunoblotting varied throughout the course of infection depending on murine strain. Nevertheless, the antibody patterns of the two strains showed significant homogeneity in being directed against almost identical subset of antigens.

Francisella tularensis is one of the most infectious human pathogens known. In the past, both the former Soviet Union and the US had programs to develop weapons containing the bacterium. We report the complete genome sequence of a highly virulent isolate of F. tularensis (1,892,819 bp). The sequence uncovers previously uncharacterized genes encoding type IV pili, a surface polysaccharide and iron-acquisition systems. Several virulence-associated genes were located in a putative pathogenicity island, which was duplicated in the genome. More than 10% of the putative coding sequences contained insertion-deletion or substitution mutations and seemed to be deteriorating. The genome is rich in IS elements, including IS630 Tc-1 mariner family transposons, which are not expected in a prokaryote. We used a computational method for predicting metabolic pathways and found an unexpectedly high proportion of disrupted pathways, explaining the fastidious nutritional requirements of the bacterium. The loss of biosynthetic pathways indicates that F. tularensis is an obligate host-dependent bacterium in its natural life cycle. Our results have implications for our understanding of how highly virulent human pathogens evolve and will expedite strategies to combat them.

The intracellular bacterium Francisella tularensis is the causative agent of tularemia and poses a serious threat as an agent of bioterrorism. We have developed a highly effective molecular subtyping system from 25 variable-number tandem repeat (VNTR) loci. In our study, multiple-locus VNTR analysis (MLVA) was used to analyze genetic relationships and potential population structure within a global collection of 192 F. tularensis isolates, including representatives from each of the four subspecies. The VNTR loci displayed between 2 and 31 alleles with Nei's diversity values between 0.05 and 0.95. Neighbor-joining cluster analysis of VNTR data revealed 120 genotypes among the 192 F. tularensis isolates, including accurate subspecies identification. F. tularensis subsp. tularensis (type A) isolates showed great diversity at VNTR loci, while F. tularensis subsp. holarctica (type B) isolates showed much lower levels despite a much broader geographical prevalence. The resolution of two distinct clades within F. tularensis subsp. tularensis (designated A.I and A.II) revealed a previously unrecognized genetic division within this highly virulent subspecies. F. tularensis subsp. holarctica appears to have recently spread globally across continents from a single origin, while F. tularensis subsp. tularensis has a long and complex evolutionary history almost exclusively in North America. The sole non-North American type A isolates (Slovakian) were closely related to the SCHU S4 strain. Significant linkage disequilibrium was detected among VNTR loci of F. tularensis consistent with a clonal population structure. Overall, this work greatly augments the study of tularemia ecology and epidemiology, while providing a framework for future forensic analysis of F. tularensis isolates.

Mosquitoes are thought to function as mechanical vectors of Francisella tularensis subspecies holarctica (F. t. holarctica) causing tularemia in humans. We investigated the clinical relevance of transstadially maintained F. t. holarctica in mosquitoes. Aedes egypti larvae exposed to a fully virulent F. t. holarctica strain for 24 hours, were allowed to develop into adults when they were individually homogenized. Approximately 24% of the homogenates tested positive for F. t. DNA in PCR. Mice injected with the mosquito homogenates acquired tularemia within 5 days. This novel finding demonstrates the possibility of transmission of bacteria by adult mosquitoes having acquired the pathogen from their aquatic larval habitats.

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

To the Editor: Tularemia is a zoonosis caused by the bacterium Francisella tularensis; the main forms of disease that occur in humans are ulceroglandular/glandular, oculoglandular, oropharyngeal, and respiratory. In Turkey, tularemia outbreaks were described as early as 1936–1938 (1), but tularemia was not reportable until 2004. Recently, multiple tularemia outbreaks in Turkey have been described, including in regions where the disease has not been previously reported; it is now considered a reemerging zoonotic disease in Turkey (1).

Second generation sequencing has revolutionized genomic studies. However, most genomes contain repeated DNA elements that are longer than the read lengths achievable with typical sequencers, so the genomic order of several generated contigs cannot be easily resolved. A new generation of sequencers offering substantially longer reads is emerging, notably the Pacific Biosciences (PacBio) RS II system and the MinION system, released in early 2014 by Oxford Nanopore Technologies through an early access program. The latter has highly advantageous portability and sequences samples by measuring changes in ionic current when single-stranded DNA molecules are translocated through nanopores. We show that the MinION system produces long reads with high mapability that can be used for scaffolding bacterial genomes, despite currently producing substantially higher error rates than PacBio reads. With further development we anticipate that MinION will be useful not only for assembling genomes, but also for rapid detection of organisms, potentially in the field.

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.

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

High throughput sequencing (HTS) generates large amounts of high quality sequence data for microbial genomics. The value of HTS for microbial forensics is the speed at which evidence can be collected and the power to characterize microbial-related evidence to solve biocrimes and bioterrorist events. As HTS technologies continue to improve, they provide increasingly powerful sets of tools to support the entire field of microbial forensics. Accurate, credible results allow analysis and interpretation, significantly influencing the course and/or focus of an investigation, and can impact the response of the government to an attack having individual, political, economic or military consequences. Interpretation of the results of microbial forensic analyses relies on understanding the performance and limitations of HTS methods, including analytical processes, assays and data interpretation. The utility of HTS must be defined carefully within established operating conditions and tolerances. Validation is essential in the development and implementation of microbial forensics methods used for formulating investigative leads attribution. HTS strategies vary, requiring guiding principles for HTS system validation. Three initial aspects of HTS, irrespective of chemistry, instrumentation or software are: 1) sample preparation, 2) sequencing, and 3) data analysis. Criteria that should be considered for HTS validation for microbial forensics are presented here. Validation should be defined in terms of specific application and the criteria described here comprise a foundation for investigators to establish, validate and implement HTS as a tool in microbial forensics, enhancing public safety and national security. Published version

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.

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. Contact: adrian.larkeryd@foi.se

A new method for monitoring phenotypic profiles of pure cultures and complex microbial communities was evaluated. The approach was to stain microorganisms with a battery of fluorescent dyes prior to flow cytometry analysis (FCM) and to analyse the data using multivariate methods, including principal component analysis and partial least squares. The FCM method was quantitatively evaluated using different mixtures of pure cultures as well as microbial communities. The results showed that the method could quantitatively and reproducibly resolve both populations and communities of microorganisms with 5% abundance in a diverse microbial background. The feasibility of monitoring complex microbial communities over time during the biodegradation of naphthalene using the FCM method was demonstrated. The biodegradation of naphthalene occurred to differing extents in microcosms representing three different types of aromatic-contaminated groundwater and a sample of bio-basin water. The FCM method distinguished each of these four microbial communities. The phenotypic profiles were compared with genotypic profiles generated by random-amplified polymorphic DNA analysis. The genotypic profiles of the microbial communities described only the microbial composition, and not their functional change, whereas the phenotypic profiles seemed to contain information on both the composition and the functional change of the microorganisms. Furthermore, event analysis of the FCM data showed that microbial communities with initially differing compositions could converge towards a similar composition if they had a capacity for high levels of degradation, whereas microbial communities with similar initial compositions could diverge if they differed in biodegrading ability.

PCR and culture were comparatively evaluated for their abilities to demonstrate Francisella tularensis in wound specimens from tularemia patients during an outbreak in Sweden in 1998. For transport of the specimens used for PCR, a buffer solution containing a nuclease inhibitor was used, and for transport of the specimens used for culture, a commercial transport system was selected after experimental comparison of various systems. Of 40 patients with culture- and/or serology-verified ulceroglandular tularemia, PCR detected F. tularensis DNA in 30 (75%) patients, whereas culture detected bacterial growth in 25 (62%) patients. Compared to data from a previous study, the present inclusion of a nuclease inhibitor in the transport medium did not improve the sensitivity of the PCR, whereas the sensitivity of the culture procedure was significantly increased by selection of the system used for transport. Among eight patients with clinically suspected tularemia but with negative serology and culture, specimens from four patients showed detectable DNA. In three of these patients the diagnosis was verified by the demonstration of an F. tularensis -specific T-cell response in vitro. In conclusion, PCR was more sensitive than culture for demonstration of F. tularensis in wound specimens. Besides, we showed that tularemia may proceed without development of serum antibodies, and in these patients, PCR may be of special importance for verification of the diagnosis.

This study demonstrates a method for describing changes in complex microbial communities without bias in ex situ cultivation. Bulk DNA was extracted and purified at regular intervals from two different model systems over a period of three weeks. Random amplified polymorphic DNA was used to generate DNA fingerprints, the polymerase chain reaction products were separated and detected by an automated sequencer and translated into data points, which were compiled into a matrix and assessed by principal component analysis. The reproducibility of the method was evaluated, and its ability to monitor overall changes in complex microbial communities was demonstrated.

Random Amplified Polymorphic DNA analysis (RAPD) is a methodology that has been used as a tool for monitoring microbial communities. To be useful in this application RAPD, and any other methodology, must show properties that allows for the detection of quantitative changes in composition of the microbiota. Therefore, the objective of this study was to establish whether RAPD possesses such properties. The strategy was to use genomic DNA, extracted from a set of tertiary bacterial mixtures defined according to an experimental mixture design, and containing varying proportions of Escherichia coli, Bacillus subtilis, and Pseudomonas CF600. RAPD-PCR was performed on the mixed DNA extracts and the amplified DNA fragments were separated on sequencing gels to produce genomic fingerprints that were digitized and modeled by Partial Least Squares regression (PLS). Significant predictions were obtained using an external test set for validation, with Root Mean Square Error of Predictions (RMSEP) of 0.21, 0.19 and 0.20 for the proportion of E. coli, B. subtilis and Pseudomonas CF600 respectively. Taken together, the results showed that RAPD patterns quantitatively represented the initial mixture proportions. Therefore, the view that RAPD could be useful for whole microbial community monitoring was strengthened.