Unprecedented genetic variability of PFam54 paralogs among Eurasian Lyme borreliosis-causing spirochetes.

Lyme borreliosis (LB) is the most common vector-borne disease in the Northern Hemisphere caused by spirochetes belonging to the Borrelia burgdorferi sensu lato (Bbsl) complex. Borrelia spirochetes circulate in obligatory transmission cycles between tick vectors and different vertebrate hosts. To successfully complete this complex transmission cycle, Bbsl encodes for an arsenal of proteins including the PFam54 protein family with known, or proposed, influences to reservoir host and/or vector adaptation. Even so, only fragmentary information is available regarding the naturally occurring level of variation in the PFam54 gene array especially in relation to Eurasian-distributed species. Utilizing whole genome data from isolates (n = 141) originated from three major LB-causing Borrelia species across Eurasia (B. afzelii, B. bavariensis, and B. garinii), we aimed to characterize the diversity of the PFam54 gene array in these isolates to facilitate understanding the evolution of PFam54 paralogs on an intra- and interspecies level. We found an extraordinarily high level of variation in the PFam54 gene array with 39 PFam54 paralogs belonging to 23 orthologous groups including five novel paralogs. Even so, the gene array appears to have remained fairly stable over the evolutionary history of the studied Borrelia species. Interestingly, genes outside Clade IV, which contains genes encoding for proteins associated with Borrelia pathogenesis, more frequently displayed signatures of diversifying selection between clades that differ in hypothesized vector or host species. This could suggest that non-Clade IV paralogs play a more important role in host and/or vector adaptation than previously expected, which would require future lab-based studies to validate.

A high fidelity approach to assembling the complex Borrelia genome.

BACKGROUND: Bacteria of the Borrelia burgdorferi sensu lato (s.l.) complex can cause Lyme borreliosis. Different B. burgdorferi s.l. genospecies vary in their host and vector associations and human pathogenicity but the genetic basis for these adaptations is unresolved and requires completed and reliable genomes for comparative analyses. The de novo assembly of a complete Borrelia genome is challenging due to the high levels of complexity, represented by a high number of circular and linear plasmids that are dynamic, showing mosaic structure and sequence homology. Previous work demonstrated that even advanced approaches, such as a combination of short-read and long-read data, might lead to incomplete plasmid reconstruction. Here, using recently developed high-fidelity (HiFi) PacBio sequencing, we explored strategies to obtain gap-free, complete and high quality Borrelia genome assemblies. Optimizing genome assembly, quality control and refinement steps, we critically appraised existing techniques to create a workflow that lead to improved genome reconstruction. RESULTS: Despite the latest available technologies, stand-alone sequencing and assembly methods are insufficient for the generation of complete and high quality Borrelia genome assemblies. We developed a workflow pipeline for the de novo genome assembly for Borrelia using several types of sequence data and incorporating multiple assemblers to recover the complete genome including both circular and linear plasmid sequences. CONCLUSION: Our study demonstrates that, with HiFi data and an ensemble reconstruction pipeline with refinement steps, chromosomal and plasmid sequences can be fully resolved, even for complex genomes such as Borrelia. The presented pipeline may be of interest for the assembly of further complex microbial genomes.

German Ixodes inopinatus samples may not actually represent this tick species.

Ticks are important vectors of human and animal pathogens, but many questions remain unanswered regarding their taxonomy. Molecular sequencing methods have allowed research to start understanding the evolutionary history of even closely related tick species. Ixodes inopinatus is considered a sister species and highly similar to Ixodes ricinus, an important vector of many tick-borne pathogens in Europe, but identification between these species remains ambiguous with disagreement on the geographic extent of I. inopinatus. In 2018-2019, 1583 ticks were collected from breeding great tits (Parus major) in southern Germany, of which 45 were later morphologically identified as I. inopinatus. We aimed to confirm morphological identification using molecular tools. Utilizing two genetic markers (16S rRNA, TROSPA) and whole genome sequencing of specific ticks (n = 8), we were able to determine that German samples, morphologically identified as I. inopinatus, genetically represent I. ricinus regardless of previous morphological identification, and most likely are not I. ricinus/I. inopinatus hybrids. Further, our results showed that the entire mitochondrial genome, let alone singular mitochondrial genes (i.e., 16S), is unable to distinguish between I. ricinus and I. inopinatus. Our results suggest that I. inopinatus is geographically isolated as a species (northern Africa and potentially southern Spain and Portugal) and brings into question whether I. inopinatus exists in central Europe. Our results highlight the probable existence of I. inopinatus and the power of utilizing genomic data in answering questions regarding tick taxonomy.

Candidatus Lariskella arthopodarum endosymbiont is the main factor differentiating the microbiome communities of female and male Borrelia-positive Ixodes persulcatus ticks.

Ixodes persulcatus, a hard-bodied tick species primarily found in Asia and Eastern Europe, is a vector of pathogens to human and livestock hosts. Little research has been done on the microbiome of this species, especially using individual non-pooled samples and comparing different geographical locations. Here, we use 16S rRNA amplicon sequencing to determine the individual microbial composition of 85 Borrelia-positive I. persulcatus from the Japanese islands of Hokkaido and Honshu. The resulting data (164 unique OTUs) were further analyzed to compare the makeup and diversity of the microbiome by sex and location, as well as to determine the presence of human pathogens. We found that, while location had little influence, the diversity of I. persulcatus microbiome was predominantly dependent on sex. Males were seen to have higher microbiome diversity than females, likely due to the high presence of endosymbiotic Candidatus Lariskella arthropodarum within the female microbial communities. Furthermore, high read counts for five genera containing potentially human pathogenic species were detected among both male and female microbiomes: Ehrlichia, Borrelia, Rickettsia, Candidatus Neoehrlichia and Burkholderia and co-infections between different pathogens were frequent. We conclude that the microbiome of I. persulcatus depends mainly on sex and not geographical location and that the major difference between sexes is due to the high abundance of Ca. L. arthropodarum in females. We also stress the importance of this tick species as a vector of potential human pathogens frequently found in co-infections.

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

Utilizing Two Borrelia bavariensis Isolates Naturally Lacking the PFam54 Gene Array To Elucidate the Roles of PFam54-Encoded Proteins.

Lyme borreliosis is the most common vector-borne disease in the Northern Hemisphere, caused by spirochetes belonging to the Borrelia burgdorferi sensu lato species complex, which are transmitted by ixodid ticks. B. burgdorferi sensu lato species produce a family of proteins on the linear plasmid 54 (PFam54), some of which confer the functions of cell adhesion and inactivation of complement, the first line of host defense. However, the impact of PFam54 in promoting B. burgdorferi sensu lato pathogenesis remains unclear because of the hurdles to simultaneously knock out all PFam54 proteins in a spirochete. Here, we describe two Borrelia bavariensis strains, PBN and PNi, isolated from patients naturally lacking PFam54 but maintaining the rest of the genome with greater than 95% identity to the reference B. bavariensis strain, PBi. We found that PBN and PNi less efficiently survive in human serum than PBi. Such defects were restored by introducing two B. bavariensis PFam54 recombinant proteins, BGA66 and BGA71, confirming the role of these proteins in providing complement evasion of B. bavariensis. Further, we found that all three strains remain detectable in various murine tissues 21 days post-subcutaneous infection, supporting the nonessential role of B. bavariensis PFam54 in promoting spirochete persistence. This study identified and utilized isolates deficient in PFam54 to associate the defects with the absence of these proteins, building the foundation to further study the role of each PFam54 protein in contributing to B. burgdorferi sensu lato pathogenesis. IMPORTANCE To establish infections, Lyme borreliae utilize various means to overcome the host's immune system. Proteins encoded by the PFam54 gene array play a role in spirochete survival in vitro and in vivo. Moreover, this gene array has been described in all currently available Lyme borreliae genomes. By investigating the first two Borrelia bavariensis isolates naturally lacking the entire PFam54 gene array, we showed that both patient isolates display an increased susceptibility to human serum, which can be rescued in the presence of two PFam54 recombinant proteins. However, both isolates remain infectious to mice after intradermal inoculation, suggesting the nonessential role of PFam54 during the long-term, but may differ slightly in the colonization of specific tissues. Furthermore, these isolates show high genomic similarity to type strain PBi (>95%) and could be used in future studies investigating the role of each PFam54 protein in Lyme borreliosis pathogenesis.

The Population Structure of Borrelia lusitaniae Is Reflected by a Population Division of Its Ixodes Vector.

Populations of vector-borne pathogens are shaped by the distribution and movement of vector and reservoir hosts. To study what impact host and vector association have on tick-borne pathogens, we investigated the population structure of Borrelia lusitaniae using multilocus sequence typing (MLST). Novel sequences were acquired from questing ticks collected in multiple North African and European locations and were supplemented by publicly available sequences at the Borrelia Pubmlst database (accessed on 11 February 2020). Population structure of B. lusitaniae was inferred using clustering and network analyses. Maximum likelihood phylogenies for two molecular tick markers (the mitochondrial 16S rRNA locus and a nuclear locus, Tick-receptor of outer surface protein A, trospA) were used to confirm the morphological species identification of collected ticks. Our results confirmed that B. lusitaniae does indeed form two distinguishable populations: one containing mostly European samples and the other mostly Portuguese and North African samples. Of interest, Portuguese samples clustered largely based on being from north (European) or south (North African) of the river Targus. As two different Ixodes species (i.e., I. ricinus and I. inopinatus) may vector Borrelia in these regions, reference samples were included for I. inopinatus but did not form monophyletic clades in either tree, suggesting some misidentification. Even so, the trospA phylogeny showed a monophyletic clade containing tick samples from Northern Africa and Portugal south of the river Tagus suggesting a population division in Ixodes on this locus. The pattern mirrored the clustering of B. lusitaniae samples, suggesting a potential co-evolution between tick and Borrelia populations that deserve further investigation.

Spatial variability in prevalence and genospecies distributions of Borrelia burgdorferi sensu lato from ixodid ticks collected in southern Germany.

Lyme borreliosis (LB) is the most common arthropod-borne disease in Europe and North America and is caused by members of the Borrelia burgdorferi sensu lato (Bbsl) species complex. These bacteria are transmitted by ixodid tick vectors and therefore human LB risk is influenced by the prevalence and distribution of Bbsl genospecies within tick vectors throughout the wild. These distributions can easily change over spatiotemporal scales and, to understand LB risk fully, up to date information on prevalence and distribution of Bbsl is required. The last survey of Bbsl in southern Germany, including parts of the Munich metropolitan area, was completed in 2006 and new data is needed. Ixodid ticks were collected in seven plots located in and around Munich, Germany, from March to July 2019 and were screened for Bbsl. Borrelia burgdorferi s. l. positive ticks (52 adults, 158 nymphs) were found in all plots and adults (0-61.5 % Bbsl positive/plot) and nymphs (17.4-59.5 % Bbsl positive/plot) did not differ significantly in their overall Bbsl prevalence. The number of Bbsl positive nymphs did vary significantly between plots but the number of positive adults did not. In total, six Bbsl genospecies were located with B. afzelii and B. garinii dominating. Additionally, the relapsing-fever species B. miyamotoi was found in two sampling plots. Our results highlight the variability in Bbsl prevalence and genospecies distribution over short geographic distances and aid in understanding LB risk in and around the Munich metropolitan area.

Ticks (Acari: Ixodidae) on birds migrating to the island of Ponza, Italy, and the tick-borne pathogens they carry.

Seasonal migration of birds between breeding and wintering areas can facilitate the spread of tick species and tick-borne diseases. In this study, 151 birds representing 10 different bird species were captured on Ponza Island, an important migratory stopover off the western coast of Italy and screened for tick infestation. Ticks were collected and identified morphologically. Morphological identification was supported through sequencing a fragment of the 16S mitochondrial gene. In total, 16 captured birds carried ticks from four tick species: Hyalomma rufipes (n = 14), Amblyomma variegatum (n = 1), Amblyomma sp. (n = 1), and Ixodes ventalloi (n = 2). All specimens were either larvae (n = 2) or nymphs (n = 16). All ticks were investigated for tick-borne pathogens using published molecular methods. Rickettsia aeschlimannii was detected in six of the 14 collected H. rufipes ticks. Additionally, the singular A. variegatum nymph tested positive for R. africae. In all 14 H. rufipes specimens (2 larvae and 12 nymphs), Francisella-like endosymbionts were detected. Four H. rufipes ticks tested positive for Borrelia burgdorferi sensu lato in a screening PCR but did not produce sufficient amplicon amounts for species identification. All ticks tested negative for tick-borne encephalitis virus, Crimean-Congo hemorrhagic fever virus, Coxiella burnetii, Coxiella-like organisms, Babesia spp., and Theileria spp. This study confirms the role of migratory birds in the spread and establishment of both exotic tick species and tick-borne pathogens outside their endemic range.

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

BACKGROUND: Borrelia 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. RESULTS: We 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). CONCLUSIONS: We 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.

Longitudinal study of prevalence and spatio-temporal distribution of Borrelia burgdorferi sensu lato in ticks from three defined habitats in Latvia, 1999-2010.

Members of the Borrelia burgdorferi sensu lato (s.l.) species complex are known to cause human Lyme borreliosis. Because of longevity of some reservoir hosts and the Ixodes tick vectors' life cycle, long-term studies are required to better understand species and population dynamics of these bacteria in their natural habitats. Ticks were collected between 1999 and 2010 in three ecologically different habitats in Latvia. We used multilocus sequence typing utilizing eight chromosomally located housekeeping genes to obtain information about species and population fluctuations and/or stability of B. burgdorferi s.l. in these habitats. The average prevalence over all years was 18.9%. From initial high-infection prevalences of 25.5%, 33.1% and 31.8%, from 2002 onwards the infection rates steadily decreased to 7.3%. Borrelia afzelii and Borrelia garinii were the most commonly found genospecies but striking local differences were obvious. In one habitat, a significant shift from rodent-associated to bird-associated Borrelia species was noted whilst in the other habitats, Borrelia species composition was relatively stable over time. Sequence types (STs) showed a random spatial and temporal distribution. These results demonstrated that there are temporal regional changes and extrapolations from one habitat to the next are not possible.