Changing geographic ranges of ticks and tick-borne pathogens: drivers, mechanisms and consequences for pathogen diversity

The geographic ranges of ticks and tick-borne pathogens are changing due to global and local environmental (including climatic) changes. In this review we explore current knowledge of the drivers for changes in the ranges of ticks and tick-borne pathogen species and strains via effects on their basic reproduction number (R-0), and the mechanisms of dispersal that allow ticks and tick-borne pathogens to invade suitable environments. Using the expanding geographic distribution of the vectors and agent of Lyme disease as an example we then investigate what could be expected of the diversity of tick-borne pathogens during the process of range expansion, and compare this with what is currently being observed. Lastly we explore how historic population and range expansions and contractions could be reflected in the phylogeography of ticks and tick-borne pathogens seen in recent years, and conclude that combined study of currently changing tick and tick-borne pathogen ranges and diversity, with phylogeographic analysis, may help us better predict future patterns of invasion and diversity

Two Boundaries Separate Borrelia burgdorferi Populations in North America

Understanding the spread of infectious diseases is crucial for implementing effective control measures. For this, it is important to obtain information on the contemporary population structure of a disease agent and to infer the evolutionary processes that may have shaped it. Here, we investigate on a continental scale the population structure of Borrelia burgdorferi, the causative agent of Lyme borreliosis (LB), a tick-borne disease, in North America. We test the hypothesis that the observed d population structure is congruent with recent population expansions and that these were preceded by bottlenecks mostly likely caused by the near extirpation in the 1900s of hosts required for sustaining tick populations. Multilocus sequence typing and complementary population analytical tools were used to evaluate B. burgdorferi samples collected in the Northeastern, Upper Midwestern, and Far-Western United States and Canada. The spatial distribution of sequence types (STs) and inferred population boundaries suggest that the current populations are geographically separated. One major population boundary separated western B. burgdorferi populations transmitted by Ixodes pacificus in California from Eastern populations transmitted by I. scapularis; the other divided Midwestern and Northeastern populations. However, populations from all three regions were genetically closely related. Together, our findings suggest that although the contemporary populations of North American B. burgdorferi now com- prise three geographically separated subpopulations with no or limited gene flow among them, they arose from a common ancestral population. A comparative analysis of the B. burgdorferi outer surface protein C (ospC) gene revealed novel linkages and provides additional insights into the genetic characteristics of strains

A Novel Rapid Sample Preparation Method for MALDI-TOF MS Permits Borrelia burgdorferi Sensu Lato Species and Isolate Differentiation

The genus Borrelia comprises vector-borne bacterial pathogens that can severely affect human and animal health. Members of the Borrelia burgdorferi sensu lato species complex can cause Lyme borreliosis, one of the most common vector-borne diseases in the Northern hemisphere. Besides, members of the relapsing fever group of spirochetes can cause tick-borne relapsing fever in humans and various febrile illnesses in animals in tropical, subtropical and temperate regions. Borrelia spp. organisms are fastidious to cultivate and to maintain in vitro, and therefore, difficult to work with in the laboratory. Currently, borrelia identification is mainly performed using PCR and DNA sequencing methods, which can be complicated/frustrating on complex DNA templates and may still be relatively expensive. Alternative techniques such as matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) are not well established for Borrelia spp., although this technique is currently one of the most used techniques for rapid identification of bacteria in microbiological diagnostic laboratories. This is mainly due to unsatisfactory results obtained by use of simple sample preparation techniques and medium-contamination obscuring the mass spectra. In addition, comprehensive libraries for Borrelia spp. MALDI-TOF MS have yet to be established. In this study, we developed a new filter-based chemical extraction technique that allows measurement of high quality Borrelia spp. spectra from less than 100,000 bacteria per spot in MALDI-TOF MS. We used 49 isolates of 13 different species to produce the largest mass-library for Borrelia spp. so far and to validate the protocol. The library was successfully established and identifies >96% of used isolates correctly to species level. Cluster analysis on the sum spectra was applied to all the different isolates, which resulted in tight cluster generation for most species. Comparative analysis of the generated cluster to a phylogeny based on concatenated multi-locus sequence typing genes provided a surprising homology. Our data demonstrate that the technique described here can be used for fast and reliable species and strain typing within the borrelia complex

Population genetics, taxonomy, phylogeny and evolution of Borrelia burgdorferi sensu lato

In order to understand the population structure and dynamics of bacterial microorganisms, typing systems that accurately reflect the phylogenetic and evolutionary relationship of the agents are required. Over the past 15 years multilocus sequence typing schemes have replaced single locus approaches, giving novel insights into phylogenetic and evolutionary relationships of many bacterial species and facilitating taxonomy. Since 2004, several schemes using multiple loci have been developed to better understand the taxonomy, phylogeny and evolution of Lyme borreliosis spirochetes and in this paper we have reviewed and summarized the progress that has been made for this important group of vector-borne zoonotic bacteria

Candidatus Borrelia kalaharica Detected from a Febrile Traveller Returning to Germany from Vacation in Southern Africa

A 26 year-old female patient presented to the Tropical Medicine outpatient unit of the Ludwig Maximilians-University in Munich with febrile illness after returning from Southern Africa, where she contracted a bite by a large mite-like arthropod, most likely a soft-tick. Spirochetes were detected in Giemsa stained blood smears and treatment was started with doxycycline for suspected tick-borne relapsing fever. The patient eventually recovered after developing a slight Jarisch-Herxheimer reaction during therapy. PCR reactions performed from EDTA-blood revealed a 16S rRNA sequence with 99.4% similarity to both, Borrelia duttonii, and B. parkeri. Further sequences obtained from the flagellin gene (flaB) demonstrated genetic distances of 0.066 and 0.097 to B. parkeri and B. duttonii, respectively. Fragments of the uvrA gene revealed genetic distance of 0.086 to B. hermsii in genetic analysis and only distant relations with classic Old World relapsing fever species. This revealed the presence of a novel species of tick-borne relapsing fever spirochetes that we propose to name "Candidatus Borrelia kalaharica", as it was contracted from an arthropod bite in the Kalahari Desert belonging to both, Botswana and Namibia, a region where to our knowledge no relapsing fever has been described so far. Interestingly, the novel species shows more homology to New World relapsing fever Borrelia such as B. parkeri or B. hermsii than to known Old World species such as B. duttonii or B. crocidurae

First insights in the variability of Borrelia recurrentis genomes

Background: Borrelia recurrentis is the causative agent of louse-borne relapsing fever, endemic to the Horn of Africa. New attention was raised in Europe, with the highest number of cases (n = 45) reported among migrants in 2015 in Germany and sporadically from other European countries. So far only one genome was sequenced, hindering the development of specific molecular diagnostic and typing tools. Here we report on modified culture conditions for B. recurrentis and the intraspecies genome variability of six isolates isolated and cultured in different years in order to explore the possibility to identify new targets for typing and examine the molecular epidemiology of the pathogen. Methodology/Principal findings: Two historical isolates from Ethiopia and four isolates from migrants from Somalia (n = 3) and Ethiopia (n = 1) obtained in 2015 were cultured in MPK-medium supplemented with 50% foetal calf serum. Whole DNA was sequenced using Illumina MiSeq technology and analysed using the CLC Genomics Workbench and SPAdes de novo assembler. Compared to the reference B. recurrentis A1 29-38 SNPs were identified in the genome distributed on the chromosome and plasmids. In addition to that, plasmids of differing length, compared to the available reference genome were identified. Conclusions/Significance: The observed low genetic variability of B. recurrentis isolates is possibly due to the adaptation to a very conserved vector-host (louse-human) cycle, or influenced by the fastidious nature of the pathogen and their resistance to in vitro growth. Nevertheless, isolates obtained in 2015 were bearing the same chromosomal SNPs and could be distinguished from the historical isolates by means of whole genome sequencing, but not hitherto used typing methods. This is the first study examining the molecular epidemiology of B. recurrentis and provides the necessary background for the development of better diagnostic tools

Next generation sequencing analysis of nine Corynebacterium ulcerans isolates reveals zoonotic transmission and a novel putative diphtheria toxin-encoding pathogenicity island

Background: Toxigenic Corynebacterium ulcerans can cause a diphtheria-like illness in humans and have been found in domestic animals, which were suspected to serve as reservoirs for a zoonotic transmission. Additionally, toxigenic C. ulcerans were reported to take over the leading role in causing diphtheria in the last years in many industrialized countries. Methods: To gain deeper insights into the tox gene locus and to understand the transmission pathway in detail, we analyzed nine isolates derived from human patients and their domestic animals applying next generation sequencing and comparative genomics. Results: We provide molecular evidence for zoonotic transmission of C. ulcerans in four cases and demonstrate the superior resolution of next generation sequencing compared to multi-locus sequence typing for epidemiologic research. Additionally, we provide evidence that the virulence of C. ulcerans can change rapidly by acquisition of novel virulence genes. This mechanism is exemplified by an isolate which acquired a prophage not present in the corresponding isolate from the domestic animal. This prophage contains a putative novel virulence factor, which shares high identity with the RhuM virulence factor from Salmonella enterica but which is unknown in Corynebacteria so far. Furthermore, we identified a putative pathogenicity island for C. ulcerans bearing a diphtheria toxin gene. Conclusion: The novel putative diphtheria toxin pathogenicity island could provide a new and alternative pathway for Corynebacteria to acquire a functional diphtheria toxin-encoding gene by horizontal gene transfer, distinct from the previously well characterized phage infection model. The novel transmission pathway might explain the unexpectedly high number of toxigenic C. ulcerans

Out of Asia? Expansion of Eurasian Lyme borreliosis causing genospecies display unique evolutionary trajectories

Vector-borne pathogens exist in obligate transmission cycles between vector and reservoir host species. Host and vector shifts can lead to geographic expansion of infectious agents and the emergence of new diseases in susceptible individuals. Three bacterial genospecies (Borrelia afzelii, Borrelia bavariensis, and Borrelia garinii) predominantly utilize two distinct tick species as vectors in Asia (Ixodes persuicatus) and Europe (Ixodes ricinus). Through these vectors, the bacteria can infect various vertebrate groups (e.g., rodents, birds) including humans where they cause Lyme borreliosis, the most common vector-borne disease in the Northern hemisphere. Yet, how and in which order the three Borrelia genospecies colonized each continent remains unclear including the evolutionary consequences of this geographic expansion. Here, by reconstructing the evolutionary history of 142 Eurasian isolates, we found evidence that the ancestors of each of the three genospecies probably have an Asian origin. Even so, each genospecies studied displayed a unique substructuring and evolutionary response to the colonization of Europe. The pattern of allele sharing between continents is consistent with the dispersal rate of the respective vertebrate hosts, supporting the concept that adaptation of Borrelia genospecies to the host is important for pathogen dispersal. Our results highlight that Eurasian Lyme borreliosis agents are all capable of geographic expansion with host association influencing their dispersal;further displaying the importance of host and vector association to the geographic expansion of vector-borne pathogens and potentially conditioning their capacity as emergent pathogens

High conservation combined with high plasticity: genomics and evolution of Borrelia bavariensis

BackgroundBorrelia bavariensis is one of the agents of Lyme Borreliosis (or Lyme disease) in Eurasia. The genome of the Borrelia burgdorferi sensu lato species complex, that includes B. bavariensis, is known to be very complex and fragmented making the assembly of whole genomes with next-generation sequencing data a challenge.ResultsWe present a genome reconstruction for 33 B. bavariensis isolates from Eurasia based on long-read (Pacific Bioscience, for three isolates) and short-read (Illumina) data. We show that the combination of both sequencing techniques allows proper genome reconstruction of all plasmids in most cases but use of a very close reference is necessary when only short-read sequencing data is available. B. bavariensis genomes combine a high degree of genetic conservation with high plasticity: all isolates share the main chromosome and five plasmids, but the repertoire of other plasmids is highly variable. In addition to plasmid losses and gains through horizontal transfer, we also observe several fusions between plasmids. Although European isolates of B. bavariensis have little diversity in genome content, there is some geographic structure to this variation. In contrast, each Asian isolate has a unique plasmid repertoire and we observe no geographically based differences between Japanese and Russian isolates. Comparing the genomes of Asian and European populations of B. bavariensis suggests that some genes which are markedly different between the two populations may be good candidates for adaptation to the tick vector, (Ixodes ricinus in Europe and I. persulcatus in Asia).ConclusionsWe present the characterization of genomes of a large sample of B. bavariensis isolates and show that their plasmid content is highly variable. This study opens the way for genomic studies seeking to understand host and vector adaptation as well as human pathogenicity in Eurasian Lyme Borreliosis agents.Peer reviewe

Evidence for host-genotype associations of Borrelia burgdorferi sensu stricto

Different genotypes of the agent of Lyme disease in North America, Borrelia burgdorferi sensu stricto, show varying degrees of pathogenicity in humans. This variation in pathogenicity correlates with phylogeny and we have hypothesized that the different phylogenetic lineages in North America reflect adaptation to different host species. In this study, evidence for host species associations of B. burgdorferi genotypes was investigated using 41 B. burgdorferi-positive samples from five mammal species and 50 samples from host-seeking ticks collected during the course of field studies in four regions of Canada: Manitoba, northwestern Ontario, Quebec, and the Maritimes. The B. burgdorferi genotypes in the samples were characterized using three established molecular markers (multi-locus sequence typing [MLST], 16S-23S rrs-rrlA intergenic spacer, and outer surface protein C sequence [ospC] major groups). Correspondence analysis and generalized linear mixed effect models revealed significant associations between B. burgdorferi genotypes and host species (in particular chipmunks, and white-footed mice and deer mice), supporting the hypotheses that host adaptation contributes to the phylogenetic structure and possibly the observed variation in pathogenicity in humans