Core genome phylogenetic analysis of the avian associated Borrelia turdi indicates a close relationship to Borrelia garinii

Natural selection and recombination at host-interacting lipoprotein loci drive genome diversification of Lyme disease and related bacteria

Lyme disease, caused by spirochetes in the Borrelia burgdorferi sensu lato clade within the Borrelia genus, is transmitted by Ixodes ticks and is currently the most prevalent and rapidly expanding tick-borne disease in Europe and North America. We report complete genome sequences of 47 isolates that encompass all established species in this clade while highlighting the diversity of the widespread human pathogenic species B. burgdorferi. A similar set of plasmids has been maintained throughout Borrelia divergence, indicating that they are a key adaptive feature of this genus. Phylogenetic reconstruction of all sequenced Borrelia genomes revealed the original divergence of Eurasian and North American lineages and subsequent dispersals that introduced B. garinii, B. bavariensis, B. lusitaniae, B. valaisiana, and B. afzelii from East Asia to Europe and B. burgdorferi and B. finlandensis from North America to Europe. Molecular phylogenies of the universally present core replicons (chromosome and cp26 and lp54 plasmids) are highly consistent, revealing a strong clonal structure. Nonetheless, numerous inconsistencies between the genome and gene phylogenies indicate species dispersal, genetic exchanges, and rapid sequence evolution at plasmid-borne loci, including key host-interacting lipoprotein genes. While localized recombination occurs uniformly on the main chromosome at a rate comparable to mutation, lipoprotein-encoding loci are recombination hotspots on the plasmids, suggesting adaptive maintenance of recombinant alleles at loci directly interacting with the host. We conclude that within- and between-species recombination facilitates adaptive sequence evolution of host-interacting lipoprotein loci and contributes to human virulence despite a genome-wide clonal structure of its natural populations.

IMPORTANCE

Lyme disease (also called Lyme borreliosis in Europe), a condition caused by spirochete bacteria of the genus Borrelia, transmitted by hard-bodied Ixodes ticks, is currently the most prevalent and rapidly expanding tick-borne disease in the United States and Europe. Borrelia interspecies and intraspecies genome comparisons of Lyme disease-related bacteria are essential to reconstruct their evolutionary origins, track epidemiological spread, identify molecular mechanisms of human pathogenicity, and design molecular and ecological approaches to disease prevention, diagnosis, and treatment. These Lyme disease-associated bacteria harbor complex genomes that encode many genes that do not have homologs in other organisms and are distributed across multiple linear and circular plasmids. The functional significance of most of the plasmid-borne genes and the multipartite genome organization itself remains unknown. Here we sequenced, assembled, and analyzed whole genomes of 47 Borrelia isolates from around the world, including multiple isolates of the human pathogenic species. Our analysis elucidates the evolutionary origins, historical migration, and sources of genomic variability of these clinically important pathogens. We have developed web-based software tools (BorreliaBase.org) to facilitate dissemination and continued comparative analysis of Borrelia genomes to identify determinants of human pathogenicity.

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.

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.

Cellular and immunological mechanisms influence host-adapted phenotypes in a vector-borne microparasite

Predicting pathogen emergence and spillover risk requires understanding the determinants of a pathogens' host range and the traits involved in host competence. While host competence is often considered a fixed species-specific trait, it may be variable if pathogens diversify across hosts. Balancing selection can lead to maintenance of pathogen polymorphisms (multiple-niche-polymorphism; MNP). The causative agent of Lyme disease, Borrelia burgdorferi (Bb), provides a model to study the evolution of host adaptation, as some Bb strains defined by their outer surface protein C (ospC) genotype, are widespread in white-footed mice and others are associated with non-rodent vertebrates (e.g. birds). To identify the mechanisms underlying potential strain × host adaptation, we infected American robins and white-footed mice, with three Bb strains of different ospC genotypes. Bb burdens varied by strain in a host-dependent fashion, and strain persistence in hosts largely corresponded to Bb survival at early infection stages and with transmission to larvae (i.e. fitness). Early survival phenotypes are associated with cell adhesion, complement evasion and/or inflammatory and antibody-mediated removal of Bb, suggesting directional selective pressure for host adaptation and the potential role of MNP in maintaining OspC diversity. Our findings will guide future investigations to inform eco-evolutionary models of host adaptation for microparasites.

Borreliae Part 1: Borrelia Lyme Group and Echidna-Reptile Group

Simple Summary

Borreliae are spirochaetes, which represent a heterogeneous phylum within bacteria. Spirochaetes are indeed distinguished from other bacteria for their spiral shape, which also characterizes Borreliae. This review describes briefly the organization of the phylum Spirocheteales with a digression about its pathogenicity and historical information about bacteria isolation and characterization. Among spirochaetes, Borrelia genus is here divided into three groups, namely the Lyme group (LG), the Echidna-Reptile group (REPG) and the Relapsing Fever group (RFG). Borreliae Part 1 deals with Lyme group and Echidna-Reptile group Borreliae, while the subject of Borreliae Part 2 is Relapsing Fever group and unclassified Borreliae. Lyme group Borreliae is organized here in sections describing ecology, namely tick vectors and animal hosts, epidemiology, microbiology, and Borrelia genome organization and antigen characterization. Furthermore, the main clinical manifestations in Lyme borreliosis are also described. Although included in the Lyme group due to their particular clinical features, Borrelia causing Baggio Yoshinari syndrome and Borrelia mayonii are described in dedicated paragraphs. The Borrelia Echidna-Reptile group has been recently characterized including spirochaetes that apparently are not pathogenic to humans, but infect reptiles and amphibians. The paragraph dedicated to this group of Borreliae describes their vectors, hosts, geographical distribution and their characteristics.

Abstract

Borreliae are divided into three groups, namely the Lyme group (LG), the Echidna-Reptile group (REPG) and the Relapsing Fever group (RFG). Currently, only Borrelia of the Lyme and RF groups (not all) cause infection in humans. Borreliae of the Echidna-Reptile group represent a new monophyletic group of spirochaetes, which infect amphibians and reptiles. In addition to a general description of the phylum Spirochaetales, including a brief historical digression on spirochaetosis, in the present review Borreliae of Lyme and Echidna-Reptile groups are described, discussing the ecology with vectors and hosts as well as microbiological features and molecular characterization. Furthermore, differences between LG and RFG are discussed with respect to the clinical manifestations. In humans, LG Borreliae are organotropic and cause erythema migrans in the early phase of the disease, while RFG Borreliae give high spirochaetemia with fever, without the development of erythema migrans. With respect of LG Borreliae, recently Borrelia mayonii, with intermediate characteristics between LG and RFG, has been identified. As part of the LG, it gives erythema migrans but also high spirochaetemia with fever. Hard ticks are vectors for both LG and REPG groups, but in LG they are mostly Ixodes sp. ticks, while in REPG vectors do not belong to that genus.

The need to unravel the twisted nature of the Borrelia burgdorferi sensu lato complex across Europe

Lyme borreliosis is a vector-borne infection caused by bacteria under the Borrelia burgdorferi sensu lato complex, both in Europe and North America. Differential gene expression at different times throughout its infectious cycle allows the spirochete to survive very diverse environments within different mammalian hosts as well as the tick vector. To date, the vast majority of data about spirochetal proteins and their functions are from genetic studies carried out on North American strains of a single species, i.e. B. burgdorferi sensu stricto. The whole-genome sequences recently obtained for several European species/strains make it feasible to adapt and use genetic techniques to study inherent differences between them. This review highlights the crucial need to undertake independent studies of genospecies within Europe, given their varying genetic content and pathogenic potential, and differences in clinical manifestation.

Evolutionary Genetics of Borrelia

The genus Borrelia consists of evolutionarily and genetically diverse bacterial species that cause a variety of diseases in humans and domestic animals. These vector-borne spirochetes can be classified into two major evolutionary groups, the Lyme borreliosis clade and the relapsing fever clade, both of which have complex transmission cycles during which they interact with multiple host species and arthropod vectors. Molecular, ecological, and evolutionary studies have each provided significant contributions towards our understanding of the natural history, biology and evolutionary genetics of Borrelia species; however, integration of these studies is required to identify the evolutionary causes and consequences of the genetic variation within and among Borrelia species. For example, molecular and genetic studies have identified the adaptations that maximize fitness components throughout the Borrelia lifecycle and enhance transmission efficacy but provide limited insights into the evolutionary pressures that have produced them. Ecological studies can identify interactions between Borrelia species and the vertebrate hosts and arthropod vectors they encounter and the resulting impact on the geographic distribution and abundance of spirochetes but not the genetic or molecular basis underlying these interactions. In this review we discuss recent findings on the evolutionary genetics from both of the evolutionarily distinct clades of Borrelia species. We focus on connecting molecular interactions to the ecological processes that have driven the evolution and diversification of Borrelia species in order to understand the current distribution of genetic and molecular variation within and between Borrelia species.

Evolutionary ecology of Lyme Borrelia

The bacterial genus, Borrelia, is comprised of vector-borne spirochete species that infect and are transmitted from multiple host species. Some Borrelia species cause highly-prevalent diseases in humans and domestic animals. Evolutionary, ecological, and molecular research on many Borrelia species have resulted in tremendous progress toward understanding the biology and natural history of these species. Yet, many outstanding questions, such as how Borrelia populations will be impacted by climate and land-use change, will require an interdisciplinary approach. The evolutionary ecology research framework incorporates theory and data from evolutionary, ecological, and molecular studies while overcoming common assumptions within each field that can hinder integration across these disciplines. Evolutionary ecology offers a framework to evaluate the ecological consequences of evolved traits and to predict how present-day ecological processes may result in further evolutionary change. Studies of microbes with complex transmission cycles, like Borrelia, which interact with multiple vertebrate hosts and arthropod vectors, are poised to leverage the power of the evolutionary ecology framework to identify the molecular interactions involved in ecological processes that result in evolutionary change. Using existing data, we outline how evolutionary ecology theory can delineate how interactions with other species and the physical environment create selective forces or impact migration of Borrelia populations and result in micro-evolutionary changes. We further discuss the ecological and molecular consequences of those micro-evolutionary changes. While many of the currently outstanding questions will necessitate new experimental designs and additional empirical data, many others can be addressed immediately by integrating existing molecular and ecological data within an evolutionary ecology framework.

Detection of Borrelia DNA in Tick Species Collected from Vegetation and Wild Animals in Fukuoka, Japan

We screened for the presence of Borrelia spp. in ticks collected from vegetation by flagging and from wild animals between May 2017 and November 2018 in Fukuoka, located in the northern Kyushu region of Japan. A total of 1,601 ticks were collected and separated based on morphology into nine species, namely Ixodes turdus, I. ovatus, Amblyomma testudinarium, Haemaphysalis flava, H. formosensis, H. kitaokai, H. longicornis, H. hystricis, and H. megaspinosa. The ticks were segregated into 561 pools and nested PCR was used to detect borrelial DNA. Borrelia turdi and Borrelia sp. HM were identified in two of the 561 pools. This is the first report of the presence of the Lyme disease group of Borrelia and of the relapsing fever group of Borrelia in Fukuoka, Japan.