Multilocus sequence analysis of Borrelia burgdorferi s.l. in Russia

Genome diversity of Borrelia garinii in marine transmission cycles does not match host associations but reflects the strains evolutionary history

Borrelia burgdorferi sensu lato is a species complex of spirochetal bacteria that occupy different ecological niches which is reflected in their reservoir host- and vector-associations. Borrelia genomes possess numerous linear and circular plasmids. Proteins encoded by plasmid genes play a major role in host- and vector-interaction and are important for Borrelia niche adaptation. However, the plasmid composition and therewith the gene repertoire may vary even in strains of a single species. Borrelia garinii, one of the six human pathogenic species, is common in Europe (vector Ixodes ricinus), Asia (vector Ixodes persulcatus) and in marine birds (vector Ixodes uriae). For the latter, only a single culture isolate (Far04) and its genome were previously available. The genome was rather small containing only one circular and six linear plasmids with a notable absence of cp32 plasmids. To further investigate B. garinii from marine transmission cycles and to explore i) whether the small number of plasmids found in isolate Far04 is a common feature in B. garinii from marine birds and presents an adaptation to this particular niche and ii) whether there may be a correlation between genome type and host species, we initiated in vitro cultures from live I. uriae collected in 2017 and 2018 from marine avian hosts and their nests. Hosts included common guillemots, Atlantic Puffin, razorbill, and kittiwake. We obtained 17 novel isolates of which 10 were sequenced using Illumina technology, one also with Pacific Bioscience technology. The 10 genomes segregated into five different genome types defined by plasmid types (based on PFam32 loci). We show that the genomes of seabird associated B. garinii contain fewer plasmids (6-9) than B. garinii from terrestrial avian species (generally ≥10), potentially suggesting niche adaptation. However, genome type did not match an association with the diverse avian seabird hosts investigated but matched the clonal complex they originated from, perhaps reflecting the isolates evolutionary history. Questions that should be addressed in future studies are (i) how is plasmid diversity related to host- and/or vector adaptation; (ii) do the different seabird species differ in reservoir host competence, and (iii) can the genome types found in seabirds use terrestrial birds as reservoir hosts.

Allelic Variants of P66 Gene in Borrelia bavariensis Isolates from Patients with Ixodid Tick-Borne Borreliosis

Protein P66 is one of the crucial virulence factors of Borrelia, inducing the production of specific antibodies in patients with ixodid tick-borne borreliosis (ITBB). Various species of Borrelia are characterized by genetic variability of the surface-exposed loop of P66. However, little is known about this variability in Borrelia bavariensis. Here we describe the variability of the nucleotide sequences of P66 gene locus in isolates of B. bavariensis. Analysis of nucleotide sequences of P66 in 27 isolates of B. bavariensis from ITBB patients revealed three allelic variants of this gene. The alignment score of amino acid sequences in the isolates showed amino acid replacements in various positions confirming the presence of three allelic variants. Two of them are characteristic only for some isolates of B. bavariensis of the Eurasian gene pool from various parts of the geographic ranges of B. bavariensis from various samples. At least three allelic variants of P66 B. bavariensis have been identified, which have different amino acid expression, occur with different frequency in ITBB patients and, presumably, can have different effects on the course of the infection.

Genotyping of tick-borne infections and determination of tick species attacking the humans in Novosibirsk and its suburbs

The south of Western Siberia traditionally belongs to the regions with an unfavorable epidemiological situation for infections transmitted by ticks. Approximately 20 thousand people annually seek medical care after a tick bite in Novosibirsk and its suburbs. The purpose of this work was to determine the species composition of ticks attacking humans in the territory of Novosibirsk and its suburbs by determining the nucleotide sequence of the fragment of the tick mitochondrial gene cytochrome oxidase COI during the spring-summer season of 2018, determination of infection of individual ticks with viral and bacterial tick-borne infections, followed by genotyping of identified pathogens by PCR methods, by determination of nucleotide sequences and phylogenetic analysis. It is established that ticks of five species attack the humans: Ixodes persulcatus, Ixodes pavlovskyi, Dermacentor reticulatus, Dermacentor marginatus, Dermacentor nuttali in the Novosibirsk metropolis. The majority of tick attacks on humans are associated with I. pavlovskyi (43.6%) and D. reticulatus (41.2%) ticks. The genetic material of the tick-borne encephalitis virus was detected in 3.6%, borrelia - 13.8% and rickettsia - 23.1% individual ticks collected from humans. Genotyping of the tick-borne encephalitis virus (TBEV) in the collected ticks showed the presence of the Siberian genotype of TBEV in 83% cases, and the Far Eastern genotype of TBEV in 17% of ticks. The Kemerovo virus could not be detected. Genetic material from Borrelia spp. and Rickettsia spp. has been detected in all five tick species attacking the local human population. The bulk of the identified Borrelia spp. isolates were genotyped as Borrelia garinii (86%) and 13% of the isolates as Borrelia afzelii species. The genetic material of Borrelia miyamotoi was found in one individual tick. Ticks of the genus Dermacentor have been found to infect by Rickettsia raoultii, and ticks of the genus Ixodes are overwhelmingly infected with Rickettsia tarasevichiae. Moreover, in single I. persulcatus tick was infected by Rickettsia helvetica. Thus, the human population of Novosibirsk and its suburbs is attacking by five species of ixodid ticks with the predominant dominance of I. pavlovskyi and D. reticulatus ticks. The genetic material of the following tick-borne infections was found: Siberian and Far Eastern genotype of TBEV, B. garinii, B. afzelii, B. miyamotoi, R. raoulti, R. tarasevichiae and R. helvetica in the individual ticks collected from humans.

Ixodes persulcatus/pavlovskyi natural hybrids in Siberia: Occurrence in sympatric areas and infection by a wide range of tick-transmitted agents

Ixodes persulcatus and Ixodes pavlovskyi ticks, two closely related species of the I. ricinus - I. persulcatus group, are widely distributed in the southern part of Western Siberia. Recently, the existence of natural hybrids of I. persulcatus and I. pavlovskyi ticks has been demonstrated. The aim of this study was to evaluate the abundance of I. persulcatus/pavlovskyi hybrids in several locations with different ratios of parental tick species and to investigate the prevalence and genetic variability of a wide range of infectious agents in these hybrids compared to the parental tick species. Natural hybrids of I. persulcatus and I. pavlovskyi ticks were identified in all examined locations in Altai and Novosibirsk, Western Siberia, Russia. The abundance of hybrids varied from 7% to 40% in different locations and was maximal in a location with similar proportions of I. persulcatus and I. pavlovskyi ticks. For the first time, it was shown that hybrids can be infected with the same agents as their parental tick species: tick-borne encephalitis and Kemerovo viruses, Borrelia afzelii, Borrelia bavariensis, Borrelia garinii, Borrelia miyamotoi, Rickettsia helvetica, Rickettsia raoultii, Rickettsia sibirica, "Candidatus Rickettsia tarasevichiae", Anaplasma phagocytophilum, Ehrlichia muris, "Candidatus Neoehrlichia mikurensis", and Babesia microti. The prevalence of most bacterial agents in hybrids was intermediate compared to their parental tick species. Most genetic variants of the identified agents have been previously found in the parental tick species. Wide distribution of I. persulcatus/pavlovskyi natural hybrids implies that I. persulcatus, I. pavlovskyi and their hybrids coexist in all I. persulcatus - I. pavlovskyi sympatric areas.

Multilocus sequence analysis of Borrelia burgdorferi sensu lato isolates from Western Siberia, Russia and Northern Mongolia

Lyme borreliosis (LB) is the most frequently recorded tick-transmitted disease in Eurasia. Tomsk Province, Western Siberia in Russia and Selenge Aimag in Northern Mongolia are leading regions in the LB incidence rate in these countries. Spirochaetes of the Borrelia burgdorferi sensu lato (s.l.) complex isolated from Ixodes ticks from Tomsk Province (n = 56) and Ixodes persulcatus ticks from Selenge Aimag (n = 5) were genetically characterized using Multi Locus Sequence Typing (MLST), analysis of the 5S23S rRNA intergenic spacer (IGS) amplicons, and p83/100 gene sequencing. According to MLST, B. afzelii (n = 26), B. bavariensis (n = 23), B. garinii (n = 11), and B. valaisiana (n = 1) isolates were detected in Tomsk Province, while B. afzelii and B. bavariensis isolates were identified in Selenge Aimag. Of the 32 revealed sequence types (ST), 21 STs were new and 14 of the new STs belonged to B. afzelii. Several STs of B. afzelii, B. garinii and B. valaisiana identified in this study clustered with European STs found in I. ricinus ticks. Analysis of the 5S23S IGS demonstrated that the studied Borrelia strains showed RFLP pattern characteristic for the following 5S23S IGS types: VS461 (B. afselii), NT29 (B. bavariensis), 20047 (B. bavariensis and B. garinii), VS116 (B. valaisiana), and three new groups (B. afzelii and B. bavariensis). Notably, this is the first report of Asian B. bavariensis possessing a 5S23S IGS RFLP pattern identical to 20047, and analysis of the 5S23S IGS did not provide correct determination of Borrelia species occurring in Asia. Genotyping of Borrelia strains using the clpA, pepX, and p83/100 genes demonstrated the same result as genotyping based on MLST; and further investigations are required to confirm that these three genetic loci could be used for determination of bacterial species from the B. burgdorferi s.l. complex because data based on single loci may be misleading.

Evidence for Borrelia bavariensis Infections of Ixodes uriae within Seabird Colonies of the North Atlantic Ocean

The first report of members of the spirochete genus Borrelia in the seabird tick, Ixodes uriae, and seabird colonies occurred during the early 1990s. Since then, Borrelia spp. have been detected in these ticks and seabird colonies around the world. To date, the primary species detected has been Borrelia garinii, with rare occurrences of Borrelia burgdorferi sensu stricto and Borrelia lusitaniae. During our research on Borrelia and I. uriae in seabird colonies of Newfoundland and Labrador, Canada, we have identified Borrelia bavariensis in I. uriae To our knowledge, B. bavariensis has previously been found only in the Eurasian tick species Ixodes persulcatus and Ixodes ricinus, and it was believed to be a rodent-specific Borrelia ecotype. We found B. bavariensis within I. uriae from three seabird colonies over three calendar years. We also reanalyzed B. garinii sequences collected from I. uriae from Eurasian seabird colonies and determined that sequences from two Russian seabird colonies likely also represent B. bavariensis The Canadian B. bavariensis sequences from I. uriae analyzed in this study cluster with previously described sequences from Asia. Overall, this is an important discovery that illustrates and expands the range of hosts and vectors for B. bavariensis, and it raises questions regarding the possible mechanisms of pathogen dispersal from Asia to North America.IMPORTANCE To our knowledge, this is the first documentation of B. bavariensis outside Eurasia. Additionally, the bacterium was found in a marine ecosystem involving the seabird tick I. uriae and its associated seabird hosts. This indicates that the epizootiology of B. bavariensis transmission is much different from what had been described, with this species previously believed to be a rodent-specific ecotype, and it indicates that this pathogen is established, or establishing, much more widely.

Recurrent evolution of host and vector association in bacteria of the Borrelia burgdorferi sensu lato species complex

Background

The Borrelia burgdorferi sensu lato (s.l.) species complex consists of tick-transmitted bacteria and currently comprises approximately 20 named and proposed genospecies some of which are known to cause Lyme Borreliosis. Species have been defined via genetic distances and ecological niches they occupy. Understanding the evolutionary relationship of species of the complex is fundamental to explaining patterns of speciation. This in turn forms a crucial basis to frame testable hypotheses concerning the underlying processes including host and vector adaptations.

Results

Illumina Technology was used to obtain genome-wide sequence data for 93 strains of 14 named genospecies of the B. burgdorferi species complex and genomic data already published for 18 additional strain (including one new species) was added. Phylogenetic reconstruction based on 114 orthologous single copy genes shows that the genospecies represent clearly distinguishable taxa with recent and still ongoing speciation events apparent in Europe and Asia. The position of Borrelia species in the phylogeny is consistent with host associations constituting a major driver for speciation. Interestingly, the data also demonstrate that vector associations are an additional driver for diversification in this tick-borne species complex. This is particularly obvious in B. bavariensis, a rodent adapted species that has diverged from the bird-associated B. garinii most likely in Asia. It now consists of two populations one of which most probably invaded Europe following adaptation to a new vector (Ixodes ricinus) and currently expands its distribution range.

Conclusions

The results imply that genotypes/species with novel properties regarding host or vector associations have evolved recurrently during the history of the species complex and may emerge at any time. We suggest that the finding of vector associations as a driver for diversification may be a general pattern for tick-borne pathogens. The core genome analysis presented here provides an important source for investigations of the underlying mechanisms of speciation in tick-borne pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3016-4) contains supplementary material, which is available to authorized users.

NGS population genetics analyses reveal divergent evolution of a Lyme Borreliosis agent in Europe and Asia

Borrelia bavariensis is a recently described agent of Lyme disease within the B. burgdorferi sensu lato species complex and exhibits a strong capacity for human pathogenicity. B. bavariensis strains are widely distributed in Eurasia spanning the distribution range of the tick vectors Ixodes persulcatus and I. ricinus. It has been suggested that B. bavariensis forms two populations, one of which arose through vector adaptation and geographic expansion. We have performed phylogenetic and population genetic analyses with next-generation sequencing data of 26 strains of B. bavariensis targeting the main linear chromosome and two plasmids (lp54, cp26). A very low number of single nucleotide polymorphisms (SNPs) was found in the European population and a deep branching pattern between European and Asian B. bavariensis was observed in all phylogenies. The results confirm the population structure of B. bavariensis and strongly support the hypothesis of clonal expansion of the European population of B. bavariensis. In addition, signals of positive selection identified in the populations further support the hypothesis that the European population of B. bavariensis likely underwent vector adaptation in its recent evolutionary history. Identified genes represent promising candidates for experimental vector adaptation studies. Thus, this species forms a very good model to study vector adaptation, which is known to play an important role in the geographic distribution of B. burgdorferi. Analysis of well known virulence determinants that are attributed to severity of clinical manifestation in B. burgdorferi s.s. revealed no variation within the European population of B. bavariensis, underlining the importance of including various Borrelia species into investigations that aim to understand the pathogenesis of Lyme disease agents.

Borrelia spirochetes in Russia: Genospecies differentiation by real-time PCR

Spirochetes of the Borrelia burgdorferi sensu lato complex are the causative agent of Lyme borreliosis which is widespread in Russia. Nowadays, three clinically important B. burgdorferi s.l. genospecies, B. afzelii, B. garinii, B. bavariensis sp. nov., can be found in Russia, as well as B. miyamotoi, which belongs to the tick-borne relapsing fever group of spirochetes. Several techniques have been developed to differentiate Borrelia genospecies. However, most of them do not allow detection of all of these genospecies simultaneously. Also, no method based on the RT-PCR TaqMan approach has been proposed to differentiate the genetically closely related species B. bavariensis and B. garinii. In the present paper, we investigated two species of ticks, I. persulcatus and I. pavlovskyi (1343 and 92 adults, respectively). Two sets of primers and probes for RT-PCR, with uvrA, glpQ and nifS genes as targets, were designed to detect four Borrelia genospecies in positive samples. The average prevalence of Borrelia sp. was about 40%, with B. afzelii as the most prevalent genospecies. Mixed infections of B. bavariensis and B. garinii were found to be extremely rare. While B. bavariensis was predominant in I. persulcatus, I. pavlovskyi ticks were infected exclusively by B. garinii. The proposed technique proved to be efficient in selection of individual Borrelia species for further genetic analysis, in particular, for multilocus sequence typing. Also, it could be applied for the differentiation of Borrelia genospecies in clinical material.

Molecular Typing of Borrelia burgdorferi

Borrelia burgdorferi sensu lato is a group of spirochetes belonging to the genus Borrelia in the family of Spirochaetaceae. The spirochete is transmitted between reservoirs and hosts by ticks of the family Ixodidae. Infection with B. burgdorferi in humans causes Lyme disease or Lyme borreliosis. Currently, 20 Lyme disease-associated Borrelia species and more than 20 relapsing fever-associated Borrelia species have been described. Identification and differentiation of different Borrelia species and strains is largely dependent on analyses of their genetic characteristics. A variety of molecular techniques have been described for Borrelia isolate speciation, molecular epidemiology, and pathogenicity studies. In this unit, we focus on three basic protocols, PCR-RFLP-based typing of the rrs-rrlA and rrfA-rrlB ribosomal spacer, ospC typing, and MLST. These protocols can be employed alone or in combination for characterization of B. burgdorferi isolates or directly on uncultivated organisms in ticks, mammalian host reservoirs, and human clinical specimens.