Revisiting Francisella tularensis subsp. holarctica, Causative Agent of Tularemia in Germany With Bioinformatics: New Insights in Genome Structure, DNA Methylation and Comparative Phylogenetic Analysis

Strengthening the genomic surveillance of Francisella tularensis by using culture-free whole-genome sequencing from biological samples

Introduction

Francisella tularensis is a highly infectious bacterium that causes the zoonotic disease tularemia. The development of genotyping methods, especially those based on whole-genome sequencing (WGS), has recently increased the knowledge on the epidemiology of this disease. However, due to the difficulties associated with the growth and isolation of this fastidious pathogen in culture, the availability of strains and subsequently WGS data is still limited.

Methods

To surpass these constraints, we aimed to implement a culture-free approach to capture and sequence F. tularensis genomes directly from complex samples. Biological samples obtained from 50 common voles and 13 Iberian hares collected in Spain were confirmed as positive for F. tularensis subsp. holarctica and subjected to a WGS target capture and enrichment protocol, using RNA oligonucleotide baits designed to cover F. tularensis genomic diversity.

Results

We obtained full genome sequences of F. tularensis from 13 animals (20.6%), two of which had mixed infections with distinct genotypes, and achieved a higher success rate when compared with culture-dependent WGS (only successful for two animals). The new genomes belonged to different clades commonly identified in Europe (B.49, B.51 and B.262) and subclades. Despite being phylogenetically closely related to other genomes from Spain, the detected clusters were often found in other countries. A comprehensive phylogenetic analysis, integrating 599 F. tularensis subsp. holarctica genomes, showed that most (sub)clades are found in both humans and animals and that closely related strains are found in different, and often geographically distant, countries.

Discussion

Overall, we show that the implemented culture-free WGS methodology yields timely, complete and high-quality genomic data of F. tularensis, being a highly valuable approach to promote and potentiate the genomic surveillance of F. tularensis and ultimately increase the knowledge on the genomics, ecology and epidemiology of this highly infectious pathogen.

Detection and Genotyping of Francisella tularensis in Animal Hosts and Vectors from Six Different Natural Landscape Areas, Gansu Province, China

Objective

Tularemia, also known as hare fever, is caused by the bacterium Francisella tularensis (F. tularensis) transmitted through diseased wild animals, blood sucking insects, or contaminated water or food, which is distributed worldwide. The purpose of this study was to investigate F. tularensis infection in animal hosts and vectors from six different natural landscape areas in Gansu Province and to identify the genotypes of the detected F. tularensis.

Methods

Rodents were captured by snap traps, and ticks were collected by dragging a cloth over the vegetation or from domestic animals. After species identification, DNA was isolated from the captured animals and detected by nested PCR assays targeting the F. tularensis fopA gene. The positive samples were further amplified to discriminate the species, and another two short-sequence tandem repeat regions (SSTR) were amplified to identify their genotypes. All positive fragments were sequenced and analyzed by ClustalX (5.0) and DNAClub software.

Results

A total of 407 rodents of 12 species were captured, among which six rodent species were positive for F. tularensis, with an overall prevalence of 3.93%. The geographical difference in infection rate was statistically significant. At the SSTR9 locus, there were 7 genotypes among positive rodent samples. A total of 1864 ticks were tested for evidence of tularemia by nested PCR assays, 69 of which were positive, with an average positive rate of 3.70% for F. tularensis in ticks. The positive rates were significantly different among different regions. Seven genotypes were identified at the SSTR9 locus, one of which seemed dominant in positive tick samples. All positive samples had the same genotype at the SSTR16 locus.

Conclusion

There is natural infection of F. tularensis among animal vectors and hosts in Gansu Province, with diverse genotypes.

Testing assembly strategies of Francisella tularensis genomes to infer an evolutionary conservation analysis of genomic structures

Background

We benchmarked sequencing technology and assembly strategies for short-read, long-read, and hybrid assemblers in respect to correctness, contiguity, and completeness of assemblies in genomes of Francisella tularensis. Benchmarking allowed in-depth analyses of genomic structures of the Francisella pathogenicity islands and insertion sequences. Five major high-throughput sequencing technologies were applied, including next-generation “short-read” and third-generation “long-read” sequencing methods.

Results

We focused on short-read assemblers, hybrid assemblers, and analysis of the genomic structure with particular emphasis on insertion sequences and the Francisella pathogenicity island. The A5-miseq pipeline performed best for MiSeq data, Mira for Ion Torrent data, and ABySS for HiSeq data from eight short-read assembly methods. Two approaches were applied to benchmark long-read and hybrid assembly strategies: long-read-first assembly followed by correction with short reads (Canu/Pilon, Flye/Pilon) and short-read-first assembly along with scaffolding based on long reads (Unicyler, SPAdes). Hybrid assembly can resolve large repetitive regions best with a “long-read first” approach.

Conclusions

Genomic structures of the Francisella pathogenicity islands frequently showed misassembly. Insertion sequences (IS) could be used to perform an evolutionary conservation analysis. A phylogenetic structure of insertion sequences and the evolution within the clades elucidated the clade structure of the highly conservative F. tularensis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08115-x.

Draft Genome Sequence of Persistent Klebsiella grimontii AT013-Mero-001, Isolated from Human Feces

Here, we report the draft genome sequence of Klebsiella grimontii AT013-Mero-001, which was isolated from feces from a sepsis patient treated with meropenem. This isolate is an antibiotic-susceptible but persistent Enterobacteriaceae strain.

Antimicrobial Susceptibility and Genomic Analysis of Aliarcobacter cibarius and Aliarcobacter thereius, Two Rarely Detected Aliarcobacter Species

Aliarcobacter cibarius and Aliarcobacter thereius are two rarely detected Aliarcobacter species. In the study, we analyzed the antimicrobial susceptibility and provide detailed insights into the genotype and phylogeny of both species using whole-genome sequencing. Thermophilic Campylobacter species are the most common bacterial foodborne pathogens causing gastroenteritis in humans worldwide. The genus Aliarcobacter is part of the Campylobacteraceae family and includes the species Aliarcobacter butzleri, Aliarcobacter cryaerophilus, Aliarcobacter skirrowii, and the rarely described Aliarcobacter cibarius, Aliarcobacter faecis, Aliarcobacter lanthieri, Aliarcobacter thereius, and Acrobarter trophiarum. Aliarcobacter are emergent enteropathogens and potential zoonotic agents. Here, we generated, analyzed, and characterized whole-genome sequences of Aliarcobacter cibarius and Aliarcobacter thereius. They were isolated from water poultry farms in Germany, cultured and identified by MALDI-TOF MS. With PCR the identity was verified. Antibiotic susceptibility testing was carried out with erythromycin, ciprofloxacin, doxycycline, tetracycline, gentamicin, streptomycin, ampicillin, and cefotaxime using the gradient strip method (E-test). Whole-genome sequences were generated including those of reference strains. Complete genomes for six selected strains are reported. These provide detailed insights into the genotype. With these, we predicted in silico known AMR genes, virulence-associated genes, and plasmid replicons. Phenotypic analysis of resistance showed differences between the presence of resistance genes and the prediction of phenotypic resistance profiles. In Aliarcobacter butzleri, the nucleotide sequence of the gyrA gene (DQ464331) can show a signature mutation resulting in an amino acid change T85>I. Acrobarter cibarius and Acrobarter thereius showed the same gene as assessed by similarity annotation of the mutations 254C>G. Most of the isolates were found to be sensitive to ciprofloxacin. The ciprofloxacin-resistant Aliarcobacter thereius isolate was associated with the amino acid change T85>I. But this was not predicted with antibiotic resistance databases, before. Ultimately, a phylogenetic analysis was done to facilitate in future outbreak analysis.

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.

Using affinity propagation clustering for identifying bacterial clades and subclades with whole-genome sequences of Francisella tularensis

By combining a reference-independent SNP analysis and average nucleotide identity (ANI) with affinity propagation clustering (APC), we developed a significantly improved methodology allowing resolving phylogenetic relationships, based on objective criteria. These bioinformatics tools can be used as a general ruler to determine phylogenetic relationships and clustering of bacteria, exemplary done with Francisella (F.) tularensis. Molecular epidemiology of F. tularensis is currently assessed mostly based on laboratory methods and molecular analysis. The high evolutionary stability and the clonal nature makes Francisella ideal for subtyping with single nucleotide polymorphisms (SNPs). Sequencing and real-time PCR can be used to validate the SNP analysis. We investigate whole-genome sequences of 155 F. tularensis subsp. holarctica isolates. Phylogenetic testing was based on SNPs and average nucleotide identity (ANI) as reference independent, alignment-free methods taking small-scale and large-scale differences within the genomes into account. Especially the whole genome SNP analysis with kSNP3.0 allowed deciphering quite subtle signals of systematic differences in molecular variation. Affinity propagation clustering (APC) resulted in three clusters showing the known clades B.4, B.6, and B.12. These data correlated with the results of real-time PCR assays targeting canSNPs loci. Additionally, we detected two subtle sub-clusters. SplitsTree was used with standard-setting using the aligned SNPs from Parsnps. Together APC, HierBAPS, and SplitsTree enabled us to generate hypotheses about epidemiologic relationships between bacterial clusters and describing the distribution of isolates. Our data indicate that the choice of the typing technique can increase our understanding of the pathogenesis and transmission of diseases with the eventual for prevention. This is opening perspectives to be applied to other bacterial species. The data provide evidence that Germany might be the collision zone where the clade B.12, also known as the East European clade, overlaps with the clade B.6, also known as the Iberian clade. Described methods allow generating a new, more detailed perspective for F. tularensis subsp. holarctica phylogeny. These results may encourage to determine phylogenetic relationships and clustering of other bacteria the same way.

Complete genome and plasmid sequences of a multidrug-resistant Salmonella enterica subsp. enterica serovar Panama isolate from a German cattle farm

We describe a rare isolate of Salmonella enterica subsp. enterica serovar Panama with an extended-spectrum β-lactamase (ESBL) profile from a German cattle-fattening farm. Applying two next-generation sequencing methods we generated sequences of the genome as well as the plasmids; assembled the draft genome sequence of Salmonella enterica subsp. enterica serovar Panama isolate 18PM0209. Antimicrobial resistance genes, virulence-associated genes and plasmids were analyzed using bioinformatics. Occurrence of multidrug-resistant Salmonella serovars at cattle-fattening farms indicate the need of enhanced surveillance to prevent further spread of these organisms.

Rapid and Culture Free Identification of Francisella in Hare Carcasses by High-Resolution Tandem Mass Spectrometry Proteotyping

Zoonotic pathogens that can be transmitted via food to humans have a high potential for large-scale emergencies, comprising severe effects on public health, critical infrastructures, and the economy. In this context, the development of laboratory methods to rapidly detect zoonotic bacteria in the food supply chain, including high-resolution mass spectrometry proteotyping are needed. In this work, an optimized sample preparation method for liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based proteome profiling was established for Francisella isolates, and a cluster analysis, as well as a phylogenetic tree, was generated to shed light on evolutionary relationships. Furthermore, this method was applied to tissues of infected hare carcasses from Germany. Even though the non-informative data outnumbered by a manifold the information of the zoonotic pathogen in the resulting proteome profiles, the standardized evaluation of MS data within an established automated analysis pipeline identified Francisella (F.) tularensis and, thus, could be, in principle, an applicable method to monitor food supply chains.

Genome sequence of a pathogenic Corynebacterium ulcerans strain isolated from a wild boar with necrotizing lymphadenitis

Objectives

Corynebacterium ulcerans can colonize a wide variety of animals and also humans are infected, typically by zoonotic transmission. Symptoms range from skin ulcers or systemic infections to diphtheria-like illness. In contrast, Corynebacterium pseudotuberculosis is widely distributed among herds of sheep, goats and other farm animals, where it causes high economic losses due to caseous lymphadenitis. Here we describe the genome sequence of an atypical C. ulcerans strain isolated from a wild boar with necrotizing lymphadenitis. This strain has similarities to C. pseudotuberculosis.

Data description

Genome sequence data of C. ulcerans isolate W25 were generated, analyzed and taxonomical relationship to other Corynebacterium species as well as growth properties of the isolate were characterized. The genome of C. ulcerans W25 comprises 2,550,924 bp with a G+C content of 54.41% and a total of 2376 genes.

Fishing in the Soup - Pathogen Detection in Food Safety Using Metabarcoding and Metagenomic Sequencing

In food safety the detection of food contaminations with pathogenic microorganisms is a race against time and often outpaced by error-prone epidemiological approaches. For evidence-based outbreak investigations fast and reliable techniques and procedures are required to identify the source of infection. Metagenomics has the potential to become a powerful tool in the field of modern food safety, since it allows the detection, identification and characterization of a broad range of pathogens in a single experiment without pre-cultivation within a couple of days. Nevertheless, sample handling, sequencing and data analysis are challenging and can introduce errors and biases into the analysis. In order to evaluate the potential of metagenomics in food safety, we generated a mock community containing DNA of foodborne bacteria. Herewith, we compare the aptitude of the two prevalent approaches - 16S rDNA amplicon sequencing and whole genome shotgun sequencing - for the detection of foodborne bacteria using different parameters during sample preparation, sequencing and data analysis. 16S rDNA sequencing did not only result in high deviations from the expected sample composition on genus and species level, but more importantly lacked the detection of several pathogenic species. While shotgun sequencing is more suitable for species detection, abundance estimation, genome assembly and species characterization, the performance can vary depending on the library preparation kit, which was confirmed for a naturally Francisella tularensis contaminated game meat sample. The application of the Nextera XT DNA Library Preparation Kit for shotgun sequencing did not only result in lower reference genome recovery and coverage, but also in distortions of the mock community composition. For data analysis, we propose a publicly available workflow for pathogen detection and characterization and demonstrate its benefits on the usability of metagenomic sequencing in food safety by analyzing an authentic metagenomic sample.

Antimicrobial Susceptibility and Genomic Structure of Arcobacter skirrowii Isolates

Campylobacter spp. are considered the most common bacterial cause of foodborne gastroenteritis in the world. The family Campylobacteraceae includes the genus Arcobacter with the three species Arcobacter butzleri, Arcobacter cryaerophilus, and Arcobacter skirrowii as emergent enteropathogens and potential zoonotic agents. Here, we characterized genome sequences of Arcobacter that were isolated from water poultry on farms in Germany. Isolates were cultured, identified by MALDI-TOF MS and identification was verified with PCR assays. Antibiotic susceptibility testing of isolates was carried out with erythromycin, ciprofloxacin, doxycycline, tetracycline, gentamicin, and streptomycin using the gradient strip method (E-test). We also sequenced whole genomes and predicted antibiotic resistance determinants, virulence factors, performed a phylogenetic analysis to determine the genetic relatedness of these isolates and searched for plasmids.