Completed genome and emergence scenario of the multidrug-resistant nosocomial pathogen Staphylococcus epidermidis ST215.

BACKGROUND: A multidrug-resistant lineage of Staphylococcus epidermidis named ST215 is a common cause of prosthetic joint infections and other deep surgical site infections in Northern Europe, but is not present elsewhere. The increasing resistance among S. epidermidis strains is a global concern. We used whole-genome sequencing to characterize ST215 from healthcare settings. RESULTS: We completed the genome of a ST215 isolate from a Swedish hospital using short and long reads, resulting in a circular 2,676,787 bp chromosome and a 2,326 bp plasmid. The new ST215 genome was placed in phylogenetic context using 1,361 finished public S. epidermidis reference genomes. We generated 10 additional short-read ST215 genomes and 11 short-read genomes of ST2, which is another common multidrug-resistant lineage at the same hospital. We studied recombination's role in the evolution of ST2 and ST215, and found multiple recombination events averaging 30-50 kb. By comparing the results of antimicrobial susceptibility testing for 31 antimicrobial drugs with the genome content encoding antimicrobial resistance in the ST215 and ST2 isolates, we found highly similar resistance traits between the isolates, with 22 resistance genes being shared between all the ST215 and ST2 genomes. The ST215 genome contained 29 genes that were historically identified as virulence genes of S. epidermidis ST2. We established that in the nucleotide sequence stretches identified as recombination events, virulence genes were overrepresented in ST215, while antibiotic resistance genes were overrepresented in ST2. CONCLUSIONS: This study features the extensive antibiotic resistance and virulence gene content in ST215 genomes. ST215 and ST2 lineages have similarly evolved, acquiring resistance and virulence through genomic recombination. The results highlight the threat of new multidrug-resistant S. epidermidis lineages emerging in healthcare settings.

FlexTaxD: flexible modification of taxonomy databases for improved sequence classification.

SUMMARY: The Flexible Taxonomy Database framework provides a method for modification and merging official and custom taxonomic databases to create improved databases. Using such databases will increase accuracy and precision of existing methods to classify sequence reads. AVAILABILITY AND IMPLEMENTATION: Source code is freely available at https://github.com/FOI-Bioinformatics/flextaxd and installable through Bioconda.

Reclassification of Francisella noatunensis subsp. orientalis Ottem et al. 2009 as Francisella orientalis sp. nov., Francisella noatunensis subsp. chilensis subsp. nov. and emended description of Francisella noatunensis.

Francisella noatunensis is a fastidious facultative intracellular bacterial pathogen that causes 'piscine francisellosis', a serious disease affecting both marine and fresh water farmed and wild fish worldwide. Currently two F. noatunensis subspecies are recognized, i.e. F. noatunensis subsp. noatunensis and F. noatunensis subsp. orientalis. In the present study, the taxonomy of F. noatunensis was revisited using a polyphasic approach, including whole genome derived parameters such as digital DNA-DNA hybridization, whole genome average nucleotide identity (wg-ANIm), whole genome phylogenetic analysis, whole genome G+C content, metabolic fingerprinting and chemotaxonomic analyses. The results indicated that isolates belonging to F. noatunensis subsp. orientalis represent a phenotypically and genetically homogenous taxon, clearly distinguishable from F. noatunensis subsp. noatunensis that fulfils requirements for separate species status. We propose, therefore, elevation of F. noatunensis subsp. orientalis to the species rank as Francisella orientalis sp. nov. with the type strain remaining as Ehime-1T (DSM 21254T=LMG 24544T). Furthermore, we identified sufficient phenotypic and genetic differences between F. noatunensis subsp. noatunensis recovered from diseased farmed Atlantic salmon in Chile and those isolated from wild and farmed Atlantic cod in Northern Europe to warrant proposal of the Chilean as a novel F. noatunensis subspecies, i.e. Francisella noatunensis subsp. chilensis subsp. nov. with strain PQ1106T (CECT 9798T=NCTC14375T) as the type strain. Finally, we emend the description of F. noatunensis by including further metabolic information and the description of atypical strains.

A New Species of the γ-Proteobacterium Francisella, F. adeliensis Sp. Nov., Endocytobiont in an Antarctic Marine Ciliate and Potential Evolutionary Forerunner of Pathogenic Species.

The study of the draft genome of an Antarctic marine ciliate, Euplotes petzi, revealed foreign sequences of bacterial origin belonging to the γ-proteobacterium Francisella that includes pathogenic and environmental species. TEM and FISH analyses confirmed the presence of a Francisella endocytobiont in E. petzi. This endocytobiont was isolated and found to be a new species, named F. adeliensis sp. nov.. F. adeliensis grows well at wide ranges of temperature, salinity, and carbon dioxide concentrations implying that it may colonize new organisms living in deeply diversified habitats. The F. adeliensis genome includes the igl and pdp gene sets (pdpC and pdpE excepted) of the Francisella pathogenicity island needed for intracellular growth. Consistently with an F. adeliensis ancient symbiotic lifestyle, it also contains a single insertion-sequence element. Instead, it lacks genes for the biosynthesis of essential amino acids such as cysteine, lysine, methionine, and tyrosine. In a genome-based phylogenetic tree, F. adeliensis forms a new early branching clade, basal to the evolution of pathogenic species. The correlations of this clade with the other clades raise doubts about a genuine free-living nature of the environmental Francisella species isolated from natural and man-made environments, and suggest to look at F. adeliensis as a pioneer in the Francisella colonization of eukaryotic organisms.

Clonality of erythromycin resistance in Francisella tularensis.

OBJECTIVES: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

An outbreak of respiratory tularemia caused by diverse clones of Francisella tularensis.

BACKGROUND: 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. METHODS: 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. RESULTS: 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 were >99.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. CONCLUSIONS: 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.

Francisella tularensis subsp. tularensis group A.I, United States.

We used whole-genome analysis and subsequent characterization of geographically diverse strains using new genetic signatures to identify distinct subgroups within Francisella tularensis subsp. tularensis group A.I: A.I.3, A.I.8, and A.I.12. These subgroups exhibit complex phylogeographic patterns within North America. The widest distribution was observed for A.I.12, which suggests an adaptive advantage.

Genome characterisation of the genus Francisella reveals insight into similar evolutionary paths in pathogens of mammals and fish.

BACKGROUND: Prior to this study, relatively few strains of Francisella had been genome-sequenced. Previously published Francisella genome sequences were largely restricted to the zoonotic agent F. tularensis. Only limited data were available for other members of the Francisella genus, including F. philomiragia, an opportunistic pathogen of humans, F. noatunensis, a serious pathogen of farmed fish, and other less well described endosymbiotic species. RESULTS: We determined the phylogenetic relationships of all known Francisella species, including some for which the phylogenetic positions were previously uncertain. The genus Francisella could be divided into two main genetic clades: one included F. tularensis, F. novicida, F. hispaniensis and Wolbachia persica, and another included F. philomiragia and F. noatunensis.Some Francisella species were found to have significant recombination frequencies. However, the fish pathogen F. noatunensis subsp. noatunensis was an exception due to it exhibiting a highly clonal population structure similar to the human pathogen F. tularensis. CONCLUSIONS: The genus Francisella can be divided into two main genetic clades occupying both terrestrial and marine habitats. However, our analyses suggest that the ancestral Francisella species originated in a marine habitat. The observed genome to genome variation in gene content and IS elements of different species supports the view that similar evolutionary paths of host adaptation developed independently in F. tularensis (infecting mammals) and F. noatunensis subsp. noatunensis (infecting fish).

Increased knowledge of Francisella genus diversity highlights the benefits of optimised DNA-based assays.

BACKGROUND: Recent advances in sequencing technologies offer promising tools for generating large numbers of genomes, larger typing databases and improved mapping of environmental bacterial diversity. However, DNA-based methods for the detection of Francisella were developed with limited knowledge about genetic diversity. This, together with the high sequence identity between several Francisella species, means there is a high risk of false identification and detection of the highly virulent pathogen Francisella tularensis. Moreover, phylogenetic reconstructions using single or limited numbers of marker sequences often result in incorrect tree topologies and inferred evolutionary distances. The recent growth in publicly accessible whole-genome sequences now allows evaluation of published genetic markers to determine optimal combinations of markers that minimise both time and laboratory costs. RESULTS: In the present study, we evaluated 38 previously published DNA markers and the corresponding PCR primers against 42 genomes representing the currently known diversity of the genus Francisella. The results highlight that PCR assays for Francisella tularensis are often complicated by low specificity, resulting in a high probability of false positives. A method to select a set of one to seven markers for obtaining optimal phylogenetic resolution or diagnostic accuracy is presented. CONCLUSIONS: Current multiple-locus sequence-typing systems and detection assays of Francisella, could be improved by redesigning some of the primers and reselecting typing markers. The use of only a few optimally selected sequence-typing markers allows construction of phylogenetic topologies with almost the same accuracy as topologies based on whole-genome sequences.

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

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

Co-Occurrence of Francisella, Spotted Fever Group Rickettsia, and Midichloria in Avian-Associated Hyalomma rufipes.

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.

Complete Genome Sequence of Francisella sp. Strain LA11-2445 (FDC406), a Novel Francisella Species Isolated from a Human Skin Lesion.

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.

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

In recent years, an increasing diversity of species has been recognized within the family Francisellaceae. Unfortunately, novel isolates are sometimes misnamed in initial publications or multiple sources propose different nomenclature for genetically highly similar isolates. Thus, unstructured and occasionally incorrect information can lead to confusion in the scientific community. Historically, detection of Francisella tularensis in environmental samples has been challenging due to the considerable and unknown genetic diversity within the family, which can result in false positive results. We have assembled a comprehensive collection of genome sequences representing most known Francisellaceae species/strains and restructured them according to a taxonomy that is based on phylogenetic structure. From this structured dataset, we identified a small number of genomic regions unique to F. tularensis that are putatively suitable for specific detection of this pathogen in environmental samples. We designed and validated specific PCR assays based on these genetic regions that can be used for the detection of F. tularensis in environmental samples, such as water and air filters.

Complete Genome Sequence of Francisella halioticida Type Strain DSM 23729 (FSC1005).

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

Complete Genome Sequences of Allofrancisella inopinata SYSU YG23 and Allofrancisella frigidaquae SYSU 10HL1970, Isolated from Water from Cooling Systems in China.

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

Genotyping of Francisella tularensis subsp. holarctica from Hares in Germany.

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.

Complete Genome Sequence of Francisella tularensis subsp. holarctica Strain A271_1 (FDC408), Isolated from a Eurasian Beaver (Castor fiber).

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.

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

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

Draft Genome Sequence of Francisella noatunensis subsp. orientalis STIR-GUS-F2f7, a Highly Virulent Strain Recovered from Diseased Red Nile Tilapia Farmed in Europe.

A highly virulent strain of Francisella noatunensis subsp. orientalis, STIR-GUS-F2f7, was isolated from moribund red Nile tilapia (Oreochromis niloticus) farmed in Europe. In this communication, the complete genome sequencing of this bacterium is reported.