Search

Abstract Introduction Drug courts are criminal justice programs to divert people with substance use disorders from incarceration into treatment. Drug courts have become increasingly popular in the US and other countries. However, their effectiveness in reducing important public health outcomes such as recidivism and substance-related health harms remains ambiguous and contested. We used nationwide register data from Sweden to evaluate the effectiveness of contract treatment sanction, the Swedish version of drug court, in reducing substance misuse, adverse somatic and mental health outcomes, and recidivism. Methods In this prospective cohort study, two quasi-experimental designs were used: difference-in-differences and the within-individual design. In the latter, we compared the risk of outcomes during time on contract treatment to, 1) parole after imprisonment and, 2) probation. Results The cohort included 11,893 individuals (13% women) who underwent contract treatment. Contract treatment was associated with a reduction of 7 percentage points (95% CI: -.088, -.055) in substance misuse, 5 percentage points (-.064, -.034) in adverse mental health events, 9 percentage points (-.113, -.076) in adverse somatic health events, and 3 fewer charges (-3.16, -2.85) for crime in difference-in-differences analyses. Within-individual associations suggested that the same individual had longer times-to-event for all outcomes during contract treatment than on parole or on probation. Conclusions Contract treatment is an effective intervention from both public health and criminal justice perspective. Our findings suggest that it is a superior alternative to incarceration in its target group. Further, we find that an implementation approach that is less punitive and more inclusive than what is typical in the US can be successful.

Abstract Monochloramine is used to regulate microbial regrowth in drinking water distribution systems (DWDS) but produces carcinogenic disinfection byproducts and constitutes a source of energy for nitrifying bacteria. This study followed biofilm-dispersed microbial communities of a full-scale DWDS distributing ultrafiltered water over three years, before and after removal of monochloramine. Communities were described using flow cytometry and amplicon sequencing, including full-length 16S rRNA gene sequencing. Removal of monochloramine increased total cell counts by up to 440%. Increased abundance of heterotrophic bacteria was followed by emergence of the predatory bacteria Bdellovibrio, and a community potentially metabolizing small organic compounds replaced the nitrifying core community. No increased abundance of Mycobacterium or Legionella was observed. Co-occurrence analysis identified a network of Nitrosomonas, Nitrospira, Sphingomonas and Hyphomicrobium, suggesting that monochloramine supported this biofilm community. While some species expanded into the changed niche, no immediate biological risk to consumers was indicated within the DWDS.

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

Abstract Background Surface raw water used as a source for drinking water production is a critical resource, sensitive to contamination. We conducted a study on Swedish raw water sources, aiming to identify mutually co‐occurring metacommunities of bacteria, and environmental factors driving such patterns. Methods The water sources were different regarding nutrient composition, water quality, and climate characteristics, and displayed various degrees of anthropogenic impact. Water inlet samples were collected at six drinking water treatment plants over 3 years, totaling 230 samples. The bacterial communities of DNA sequenced samples (n = 175), obtained by 16S metabarcoding, were analyzed using a joint model for taxa abundance. Results Two major groups of well‐defined metacommunities of microorganisms were identified, in addition to a third, less distinct, and taxonomically more diverse group. These three metacommunities showed various associations to the measured environmental data. Predictions for the well‐defined metacommunities revealed differing sets of favored metabolic pathways and life strategies. In one community, taxa with methanogenic metabolism were common, while a second community was dominated by taxa with carbohydrate and lipid‐focused metabolism. Conclusion The identification of ubiquitous persistent co‐occurring bacterial metacommunities in freshwater habitats could potentially facilitate microbial source tracking analysis of contamination issues in freshwater sources.

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

ObjectivesTo examine differences in recidivism rates between different prisons using two designs-between-individual and within-individual-to account for confounding factors.MethodsWe examined recidivism rates among 37,891 individuals released from 44 Swedish prisons in three security levels, and who were followed from 2006 to 2013. We used longitudinal data from nationwide registers, including all convictions from district courts. First, we applied a between-individual design (Cox proportional hazards regression), comparing reconviction rates between individuals released from prisons within the same security level, while adjusting for a range of individual-level covariates. Second, we applied a within-individual design (stratified Cox proportional hazards regression), comparing rates of reconviction within the same individuals, i.e., we compared rates after release from one prison to the rates in the same individual after release from another prison, thus adjusting for all time-invariant confounders within each individual (e.g. genetics and early environment). We also adjusted for a range of time-varying individual-level covariates.ResultsResults showed differences in the hazard of recidivism between different prisons in between-individual analyses, with hazards ranging from 1.22 (1.05-1.43) to 4.99 (2.44-10.21). Results from within-individual analyses, which further adjusted for all time-invariant confounders, showed minimal differences between prisons, with hazards ranging from 0.95 (0.87-1.05) to 1.05 (0.95-1.16). Only small differences were found when violent and non-violent crimes were analyzed separately.ConclusionsThe study highlights the importance of research designs that more fully adjust for individual-level confounding factors to avoid over-interpretation of the variability in comparisons across prisons.

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

Difficulty in controlling SARS-CoV-2 transmission made the ability to inactivate viruses in aerosols and fomites to be an important and attractive risk reduction measure. Evidence that light frequencies have the ability to inhibit microorganisms has already been reported by many studies which, however, focused on ultraviolet (UV) wavelengths, which are known to induce potential injury in humans. In the present study, the effect on suspensions of SARS-CoV-2 of a Light Emitting Diode (LED) device capable of radiating frequencies in the non-hazardous visible light spectrum (VIS) was investigated. In order to evaluate the efficiency of viral inactivation, plaque assay and western blot of viral proteins were performed. The observed results showed a significant reduction in infectious particles that had been exposed to the LED irradiation of visible light. Furthermore, the analysis of the intracellular expression of viral proteins confirmed the inactivating effect of this irradiation technology. This in vitro study revealed for the first time the inactivation of SARS-CoV-2 through LED irradiation with multiple wavelengths of the visible spectrum. However additional and more in-depth studies can aim to demonstrate the data obtained during these experiments in different matrices, in mutable environmental conditions and on other respiratory viruses such as the influenza virus. The type of LED technology can decisively contribute on reducing virus transmission through the continuous sanitation of common environments without risks for humans and animals.

The migratory behavior of wild birds contributes to the geographical spread of ticks and their microorganisms. In this study, we aimed to investigate the dispersal and co-occurrence of Francisella and spotted fever group Rickettsia (SFGR) in ticks infesting birds migrating northward in the African-Western Palaearctic region (AWPR). Birds were trapped with mist nests across the Mediterranean basin during the 2014 and 2015 spring migration. In total, 575 ticks were collected from 244 birds. We screened the ticks for the species Francisella tularensis, the genus Francisella, and SFGR by microfluidic real-time PCR. Confirmatory analyses and metagenomic sequencing were performed on tick samples that putatively tested positive for F. tularensis during initial screenings. Hyalomma rufipes was the most common tick species and had a high prevalence of Francisella, including co-occurrence of Francisella and SFGR. Metagenomic analysis of total DNA extracted from two H. rufipes confirmed the presence of Francisella, Rickettsia, and Midichloria. Average nucleotide identity and phylogenetic inference indicated the highest identity of the metagenome-assembled genomes to a Francisella-like endosymbiont (FLE), Rickettsia aeschlimannii, and Midichloria mitochondrii. The results of this study suggest that (i) FLE- and SFGR-containing ticks are dispersed by northbound migratory birds in the AWPR, (ii) H. rufipes likely is not involved in transmission of F. tularensis in the AWPR, and (iii) a dual endosymbiosis of FLEs and Midichloria may support some of the nutritional requirements of H. rufipes.

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.

Sequence-based typing of Francisella tularensis has led to insights in the evolutionary developments of tularemia. In Europe, two major basal clades of F. tularensis subsp. holarctica exist, with a distinct geographical distribution. Basal clade B.6 is primarily found in Western Europe, while basal clade B.12 occurs predominantly in the central and eastern parts of Europe. There are indications that tularemia is geographically expanding and that strains from the two clades might differ in pathogenicity, with basal clade B.6 strains being potentially more virulent than basal clade B.12. This study provides information on genotypes detected in the Netherlands during 2011–2017. Data are presented for seven autochthonous human cases and for 29 European brown hares (Lepus europaeus) with laboratory confirmed tularemia. Associated disease patterns are described for 25 European brown hares which underwent post-mortem examination. The basal clades B.6 and B.12 are present both in humans and in European brown hares in the Netherlands, with a patchy geographical distribution. For both genotypes the main pathological findings in hares associated with tularemia were severe (sub)acute necrotizing hepatitis and splenitis as well as necrotizing lesions and hemorrhages in several other organs. Pneumonia was significantly more common in the B.6 than in the B.12 cases. In conclusion, the two major basal clades present in different parts in Europe are both present in the Netherlands. In hares found dead, both genotypes were associated with severe acute disease affecting multiple organs. Hepatitis and splenitis were common pathological findings in hares infected with either genotype, but pneumonia occurred significantly more frequently in hares infected with the B.6 genotype compared to hares infected with the B.12 genotype.

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.

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—obtained mainly from hares (n = 233), other vertebrate animals, and ticks, but also from humans (n = 3)—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.

Microbial source tracking (MST) analysis is essential to identifying and mitigating the fecal pollution of water resources. The signature-based MST method uses a library of sequences to identify contaminants based on operational taxonomic units (OTUs) that are unique to a certain source. However, no clear guidelines for how to incorporate OTU overlap or natural variation in the raw water bacterial community into MST analyses exist. We investigated how the inclusion of bacterial overlap between sources in the library affects source prediction accuracy. To achieve this, large-scale sampling – including feces from seven species, raw sewage, and raw water samples from water treatment plants – was followed by 16S rRNA amplicon sequencing. The MST library was defined using three settings: (i) no raw water communities represented; (ii) raw water communities selected through clustering analysis; and (iii) local water communities collected across consecutive years. The results suggest that incorporating either the local background or representative bacterial composition improves MST analyses, as the results were positively correlated to measured levels of fecal indicator bacteria and the accuracy at which OTUs were assigned to the correct contamination source increased fourfold. Using the proportion of OTUs with high source origin probability, underpinning a contaminating signal, is a solid foundation in a framework for further deciphering and comparing contaminating signals derived in signature-based MST approaches. In conclusion, incorporating background bacterial composition of water in MST can improve mitigation efforts for minimizing the spread of pathogenic and antibiotic resistant bacteria into essential freshwater resources.

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.

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.

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.

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 identified mosquito species. Further analysis confirmed the presence of F. tularensis holarctica–specific 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.

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.

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.

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

We report here the first case of indigenous tularemia detected in The Netherlands, a nonendemic country, since 1953. Whole genome DNA sequence analysis assigned the isolate BD11-00177 to the genomic group B.FTNF002-00, which previously has been exclusively reported from Spain, France, Italy, Switzerland, and Germany. The patient had not been abroad for years, which implies that this is an indigenous infection. The current case might predict an upcoming distribution of Francisella tularensis holarctica genomic group B.FTNF002-00 in Europe.

The F. tularensis type A strain FSC198 from Slovakia and a second strain FSC043, which has attenuated virulence, are both considered to be derivatives of the North American F. tularensis type A strain SCHU S4. These strains have been propagated under different conditions: the FSC198 has undergone natural propagation in the environment, while the strain FSC043 has been cultivated on artificial media in laboratories. Here, we have compared the genome sequences of FSC198, FSC043, and SCHU S4 to explore the possibility that the contrasting propagation conditions may have resulted in different mutational patterns. We found four insertion/deletion events (INDELs) in the strain FSC043, as compared to the SCHU S4, while no single nucleotide polymorphisms (SNPs) or variable number of tandem repeats (VNTRs) were identified. This result contrasts with previously reported findings for the strain FSC198, where eight SNPs and three VNTR differences, but no INDELs exist as compared to the SCHU S4 strain. The mutations detected in the laboratory and naturally propagated type A strains, respectively, demonstrate distinct patterns supporting that analysis of mutational spectra might be a useful tool to reveal differences in past growth conditions. Such information may be useful to identify leads in a microbial forensic investigation.

A robust, rapid and flexible real-time PCR assay for hierarchical genetic typing of clinical and environmental isolates of Francisella is presented. Typing markers were found by multiple genome and gene comparisons, from which 23 canonical single nucleotide polymorphisms (canSNPs) and 11 canonical insertion-deletion mutations (canINDELs) were selected to provide phylogenetic guidelines for classification from genus to isolate level. The specificity of the developed assay, which uses 68 wells of a 96-well real-time PCR format with a detection limit of 100 pg DNA, was assessed using 62 Francisella isolates of diverse genetic and geographical origins. It was then successfully used for typing 14 F. tularensis subsp. holarctica isolates obtained from tularemia patients in Sweden in 2008 and five more genetically diverse Francisella isolates of global origins. When applied to human ulcer specimens for direct pathogen detection the results were incomplete due to scarcity of DNA, but sufficient markers were identified to detect fine-resolution differences among F. tularensis subsp. holarctica isolates causing infection in the patients. In contrast to other real-time PCR assays for Francisella, which are typically designed for specific detection of a species, subspecies, or strain, this type of assay can be easily tailored to provide appropriate phylogenetic and/or geographical resolution to meet the objectives of the analysis.

Francisella tularensis is a potent mammalian pathogen well adapted to intracellular habitats, whereas F. novicida and F. philomiragia are less virulent in mammals and appear to have less specialized lifecycles. We explored adaptations within the genus that may be linked to increased host association, as follows. First, we determined the genome sequence of F. tularensis subsp. mediasiatica, the only subspecies that had not been previously sequenced. This genome, and those of 12 other F. tularensis isolates, were then compared to the genomes of F. novicida (three isolates) and F. philomiragia (one isolate). Signs of homologous recombination were found in approximately 19.2% of F. novicida and F. philomiragia genes, but none among F. tularensis genomes. In addition, random insertions of insertion sequence elements appear to have provided raw materials for secondary adaptive mutations in F. tularensis, e.g. for duplication of the Francisella Pathogenicity Island and multiplication of a putative glycosyl transferase gene. Further, the five major genetic branches of F. tularensis seem to have converged along independent routes towards a common gene set via independent losses of gene functions. Our observations suggest that despite an average nucleotide identity of >97%, F. tularensis and F. novicida have evolved as two distinct population lineages, the former characterized by clonal structure with weak purifying selection, the latter by more frequent recombination and strong purifying selection. F. tularensis and F. novicida could be considered the same bacterial species, given their high similarity, but based on the evolutionary analyses described in this work we propose retaining separate species names.

Tularemia is a geographically widespread, severely debilitating, and occasionally lethal disease in humans. It is caused by infection by a gram-negative bacterium, Francisella tularensis. In order to better understand its potency as an etiological agent as well as its potential as a biological weapon, we have completed draft assemblies and report the first complete genomic characterization of five strains belonging to the following different Francisella subspecies (subsp.): the F. tularensis subsp. tularensis FSC033, F. tularensis subsp. holarctica FSC257 and FSC022, and F. tularensis subsp. novicida GA99-3548 and GA99-3549 strains. Here, we report the sequencing of these strains and comparative genomic analysis with recently available public Francisella sequences, including the rare F. tularensis subsp. mediasiatica FSC147 strain isolate from the Central Asian Region. We report evidence for the occurrence of large-scale rearrangement events in strains of the holarctica subspecies, supporting previous proposals that further phylogenetic subdivisions of the Type B clade are likely. We also find a significant enrichment of disrupted or absent ORFs proximal to predicted breakpoints in the FSC022 strain, including a genetic component of the Type I restriction-modification defense system. Many of the pseudogenes identified are also disrupted in the closely related rarely human pathogenic F. tularensis subsp. mediasiatica FSC147 strain, including modulator of drug activity B (mdaB) (FTT0961), which encodes a known NADPH quinone reductase involved in oxidative stress resistance. We have also identified genes exhibiting sequence similarity to effectors of the Type III (T3SS) and components of the Type IV secretion systems (T4SS). One of the genes, msrA2 (FTT1797c), is disrupted in F. tularensis subsp. mediasiatica and has recently been shown to mediate bacterial pathogen survival in host organisms. Our findings suggest that in addition to the duplication of the Francisella Pathogenicity Island, and acquisition of individual loci, adaptation by gene loss in the more recently emerged tularensis, holarctica, and mediasiatica subspecies occurred and was distinct from evolutionary events that differentiated these subspecies, and the novicida subspecies, from a common ancestor. Our findings are applicable to future studies focused on variations in Francisella subspecies pathogenesis, and of broader interest to studies of genomic pathoadaptation in bacteria.

Francisella tularensis causes tularaemia, a life-threatening zoonosis, and has potential as a biowarfare agent. F. tularensis subsp. tularensis, which causes the most severe form of tularaemia, is usually confined to North America. However, a handful of isolates from this subspecies was obtained in the 1980s from ticks and mites from Slovakia and Austria. Our aim was to uncover the origins of these enigmatic European isolates.We determined the complete genome sequence of FSC198, a European isolate of F. tularensis subsp. tularensis, by whole-genome shotgun sequencing and compared it to that of the North American laboratory strain Schu S4. Apparent differences between the two genomes were resolved by re-sequencing discrepant loci in both strains. We found that the genome of FSC198 is almost identical to that of Schu S4, with only eight SNPs and three VNTR differences between the two sequences. Sequencing of these loci in two other European isolates of F. tularensis subsp. tularensis confirmed that all three European isolates are also closely related to, but distinct from Schu S4.The data presented here suggest that the Schu S4 laboratory strain is the most likely source of the European isolates of F. tularensis subsp. tularensis and indicate that anthropogenic activities, such as movement of strains or animal vectors, account for the presence of these isolates in Europe. Given the highly pathogenic nature of this subspecies, the possibility that it has become established wild in the heartland of Europe carries significant public health implications.