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.

Genetic Traces of the Francisella tularensis Colonization of Spain, 1998-2020.

More than 1000 humans have acquired the febrile disease tularemia in Spain since the first notification of human cases in 1997. We here aimed to study the recent molecular evolution of the causative bacterium Francisella tularensis during disease establishment in Spain. Single-nucleotide polymorphisms (SNPs) and variable-number tandem repeats (VNTRs) were analyzed in whole-genome sequences (WGS) of F. tularensis. Short-read WGS data for 20 F. tularensis strains from humans infected in the periods 2014-2015 and 2018-2020 in Spain were generated. These data were combined with WGS data of 25 Spanish strains from 1998 to 2008 and two reference strains. Capillary electrophoresis data of VNTR genetic regions were generated and compared with the WGS data for the 11 strains from 2014 to 2015. Evolutionary relationships among strains were analyzed by phylogenetic methods. We identified 117 informative SNPs in a 1,577,289-nucleotide WGS alignment of 47 F. tularensis genomes. Forty-five strains from Spain formed a star-like SNP phylogeny with six branches emerging from a basal common node. The most recently evolved genomes formed four additional star-like structures that were derived from four branches of the basal common node. VNTR copy number variation was detected in two out of 10 VNTR regions examined. Genetic clustering of strains by VNTRs agreed with the clustering by SNPs. The SNP data provided higher resolution among strains than the VNTRs data in all but one cases. There was an excellent correlation between VNTR marker sizing by capillary electrophoresis and prediction from WGS data. The genetic data strongly support that tularemia, indeed, emerged recently in Spain. Distinct genetic patterns of local F. tularensis population expansions imply that the pathogen has colonized a previously disease-free geographical area. We also found that genome-wide SNPs provide higher genetic resolution among F. tularensis genomes than the use of VNTRs, and that VNTR copy numbers can be accurately predicted using short-read WGS data.

Molecular Survey of Rickettsia spp., Anaplasma spp., Ehrlichia spp., Bartonella spp., and Borrelia spp. in Fleas and Lice in Ethiopia.

Bacterial arthropod-borne pathogens can often cause fever in Africa, but rural laboratories in these settings are usually too basic to provide a precise picture of their epidemiological impact. Our aim was to determine the prevalence of bacterial pathogens in fleas and lice in a rural area of southeast Ethiopia. Between July and November 2013, we extracted DNA from 91 fleas (Ctenocephalides felis [n = 50; 54.9%], Pulex irritans [n = 37; 40.1%], and C. canis [n = 4; 4.4%] and 30 lice (Pediculus humanus capitis [n = 16; 53.3%] and Pediculus humanus humanus [n = 14; 46.7%]), using two quantitative PCR (qPCR) analyses to look for bacteria from the genera: Anaplasma, Bartonella, Borrelia, Coxiella, Ehrlichia, Francisella, and Rickettsia. Of the 91 fleas analyzed, pathogens were present in 79 (86.8%), including Rickettsia felis (n = 41; 45%), Anaplasma platys (n = 40; 44.0%), Rickettsia monacensis (n = 2; 2.2%), Ehrlichia muris-like agent (n = 1; 1.1%), and Bartonella clarridgeiae (n = 1; 1.1%). P. irritans was the flea species most frequently infected with A. platys (67.7%), followed by C. felis (30.7%) (p < 0.001). Of the 30 lice identified, pathogens were present in 7 (23.3%): Bartonella quintana (n = 4; 16.7%), E. muris (n = 2, 6.7%), and Borrelia recurrentis (n = 1, 3.3%). Thus, in this rural area of Africa, fleas and lice can transmit parasitic pathogens to humans, causing febrile symptoms.

Zoonotic Bacteria in Fleas Parasitizing Common Voles, Northwestern Spain.

We detected Francisella tularensis and Bartonella spp. in fleas parasitizing common voles (Microtus arvalis) from northwestern Spain; mean prevalence was 6.1% for F. tularensis and 51% for Bartonella spp. Contrasted vector-host associations in the prevalence of these bacteria suggest that fleas have distinct roles in the transmission cycle of each pathogen in nature.