Co-evolution of the outer surface protein C gene (ospC) and intraspecific lineages of Borrelia burgdorferi sensu stricto in the northeastern United States

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.

Population dynamics of the Lyme disease bacterium, Borrelia burgdorferi, during rapid range expansion in New York State

Recent changes in climate and human land-use have resulted in alterations of the geographic range of many species, including human pathogens. Geographic range expansion and population growth of human pathogens increase human disease risk. Relatively little empirical work has investigated the impact of range changes on within-population variability, a contributor to both colonization success and adaptive potential, during the precise time in which populations are colonized. This is likely due to the difficulties of collecting appropriate natural samples during the dynamic phase of migration and colonization. We systematically collected blacklegged ticks (Ixodes scapularis) across New York State (NY), USA, between 2006 and 2019, a time period coinciding with a rapid range expansion of ticks and their associated pathogens including Borrelia burgdorferi, the etiological agent of Lyme disease. These samples provide a unique opportunity to investigate the population dynamics of human pathogens as they expand into novel territory. We observed that founder effects were short-lived, as gene flow from long-established populations brought almost all B. burgdorferi lineages to newly colonized populations within just a few years of colonization. By 7 years post-colonization, B. burgdorferi lineage frequency distributions were indistinguishable from long-established sites, indicating that local B. burgdorferi populations experience similar selective pressures despite geographic separation. The B. burgdorferi lineage dynamics elucidate the processes underlying the range expansion and demonstrate that migration into, and selection within, newly colonized sites operate on different time scales.

Whole genome sequence and comparative analysis of Borrelia burgdorferi MM1

Lyme disease is caused by spirochaetes of the Borrelia burgdorferi sensu lato genospecies. Complete genome assemblies are available for fewer than ten strains of Borrelia burgdorferi sensu stricto, the primary cause of Lyme disease in North America. MM1 is a sensu stricto strain originally isolated in the midwestern United States. Aside from a small number of genes, the complete genome sequence of this strain has not been reported. Here we present the complete genome sequence of MM1 in relation to other sensu stricto strains and in terms of its Multi Locus Sequence Typing. Our results indicate that MM1 is a new sequence type which contains a conserved main chromosome and 15 plasmids. Our results include the first contiguous 28.5 kb assembly of lp28-8, a linear plasmid carrying the vls antigenic variation system, from a Borrelia burgdorferi sensu stricto strain.

Reservoir and vector evolutionary pressures shaped the adaptation of Borrelia

The life cycle of spirochetes of the genus Borrelia includes complex networks of vertebrates and ticks. The tripartite association of Borrelia-vertebrate-tick has proved ecologically successful for these bacteria, which have become some of the most prominent tick-borne pathogens in the northern hemisphere. To keep evolutionary pace with its double-host life history, Borrelia must adapt to the evolutionary pressures exerted by both sets of hosts. In this review, we attempt to reconcile functional, phylogenetic, and ecological perspectives to propose a coherent scenario of Borrelia evolution. Available empirical information supports that the association of Borrelia with ticks is very old. The major split between the tick families Argasidae-Ixodidae (dated some 230-290 Mya) resulted in most relapsing fever (Rf) species being restricted to Argasidae and few associated with Ixodidae. A further key event produced the diversification of the Lyme borreliosis (Lb) species: the radiation of ticks of the genus Ixodes from the primitive stock of Ixodidae (around 217 Mya). The ecological interactions of Borrelia demonstrate that Argasidae-transmitted Rf species remain restricted to small niches of one tick species and few vertebrates. The evolutionary pressures on this group are consequently low, and speciation processes seem to be driven by geographical isolation. In contrast to Rf, Lb species circulate in nested networks of dozens of tick species and hundreds of vertebrate species. This greater variety confers a remarkably variable pool of evolutionary pressures, resulting in large speciation of the Lb group, where different species adapt to circulate through different groups of vertebrates. Available data, based on ospA and multilocus sequence typing (including eight concatenated in-house genes) phylogenetic trees, suggest that ticks could constitute a secondary bottleneck that contributes to Lb specialization. Both sets of adaptive pressures contribute to the resilience of highly adaptable meta-populations of bacteria.

Plasmid diversity and phylogenetic consistency in the Lyme disease agent Borrelia burgdorferi

Background

Bacteria from the genus Borrelia are known to harbor numerous linear and circular plasmids. We report here a comparative analysis of the nucleotide sequences of 236 plasmids present in fourteen independent isolates of the Lyme disease agent B. burgdorferi.

Results

We have sequenced the genomes of 14 B. burgdorferi sensu stricto isolates that carry a total of 236 plasmids. These individual isolates carry between seven and 23 plasmids. Their chromosomes, the cp26 and cp32 circular plasmids, as well as the lp54 linear plasmid, are quite evolutionarily stable; however, the remaining plasmids have undergone numerous non-homologous and often duplicative recombination events. We identify 32 different putative plasmid compatibility types among the 236 plasmids, of which 15 are (usually) circular and 17 are linear. Because of past rearrangements, any given gene, even though it might be universally present in these isolates, is often found on different linear plasmid compatibility types in different isolates. For example, the arp gene and the vls cassette region are present on plasmids of four and five different compatibility types, respectively, in different isolates. A majority of the plasmid types have more than one organizationally different subtype, and the number of such variants ranges from one to eight among the 18 linear plasmid types. In spite of this substantial organizational diversity, the plasmids are not so variable that every isolate has a novel version of every plasmid (i.e., there appears to be a limited number of extant plasmid subtypes).

Conclusions

Although there have been many past recombination events, both homologous and nonhomologous, among the plasmids, particular organizational variants of these plasmids correlate with particular chromosomal genotypes, suggesting that there has not been rapid horizontal transfer of whole linear plasmids among B. burgdorferi lineages. We argue that plasmid rearrangements are essentially non-revertable and are present at a frequency of only about 0.65% that of single nucleotide changes, making rearrangement-derived novel junctions (mosaic boundaries) ideal phylogenetic markers in the study of B. burgdorferi population structure and plasmid evolution and exchange.

Electronic supplementary material

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

Host associations and genomic diversity of Borrelia hermsii in an endemic focus of tick-borne relapsing fever in western North America

Background

An unrecognized focus of tick-borne relapsing fever caused by Borrelia hermsii was identified in 2002 when five people became infected on Wild Horse Island in Flathead Lake, Montana. The terrestrial small mammal community on the island is composed primarily of pine squirrels (Tamiasciurus hudsonicus) and deer mice (Peromyscus maniculatus), neither of which was known as a natural host for the spirochete. Thus a 3-year study was performed to identify small mammals as hosts for B. hermsii.

Methods

Small mammals were captured alive on two island and three mainland sites, blood samples were collected and examined for spirochetes, and serological tests performed to detect anti-B. hermsii antibodies. Ornithodoros hermsi ticks were collected and fed on laboratory mice to assess infection. Genomic DNA samples from spirochetes isolated from infected mammals and ticks were analyzed by multilocus sequence typing.

Results

Eighteen pine squirrels and one deer mouse had detectable spirochetemias when captured, from which 12 isolates of B. hermsii were established. Most pine squirrels were seropositive, and the five species of sciurids combined had a significantly higher prevalence of seropositive animals than did the other six small mammal species captured. The greater diversity of small mammals on the mainland in contrast to the islands demonstrated that other species in addition to pine squirrels were also involved in the maintenance of B. hermsii at Flathead Lake. Ornithodoros hermsi ticks produced an additional 12 isolates of B. hermsii and multilocus sequence typing identified both genomic groups of B. hermsii described previously, and identified a new genomic subdivision. Experimental infections of deer mice with two strains of B. hermsii demonstrated that these animals were susceptible to infection with spirochetes belonging to Genomic Group II but not Genomic Group I.

Conclusions

Pine squirrels are the primary hosts for the maintenance of B. hermsii on the islands in Flathead Lake, however serological evidence showed that numerous additional species are also involved on the mainland. Future studies testing the susceptibility of several small mammal species to infection with different genetic types of B. hermsii will help define their role as hosts in this and other endemic foci.

Vaccination with the variable tick protein of the relapsing fever spirochete Borrelia hermsii protects mice from infection by tick-bite

Background

Tick-borne relapsing fevers of humans are caused by spirochetes that must adapt to both warm-blooded vertebrates and cold-blooded ticks. In western North America, most human cases of relapsing fever are caused by Borrelia hermsii, which cycles in nature between its tick vector Ornithodoros hermsi and small mammals such as tree squirrels and chipmunks. These spirochetes alter their outer surface by switching off one of the bloodstream-associated variable major proteins (Vmps) they produce in mammals, and replacing it with the variable tick protein (Vtp) following their acquisition by ticks. Based on this reversion to Vtp in ticks, we produced experimental vaccines comprised on this protein and tested them in mice challenged by infected ticks.

Methods

The vtp gene from two isolates of B. hermsii that encoded antigenically distinct types of proteins were cloned, expressed, and the recombinant Vtp proteins were purified and used to vaccinate mice. Ornithodoros hermsi ticks that were infected with one of the two strains of B. hermsii from which the vtp gene originated were used to challenge mice that received one of the two Vtp vaccines or only adjuvant. Mice were then followed for infection and seroconversion.

Results

The Vtp vaccines produced protective immune responses in mice challenged with O. hermsi ticks infected with B. hermsii. However, polymorphism in Vtp resulted in mice being protected only from the spirochete strain that produced the same Vtp used in the vaccine; mice challenged with spirochetes producing the antigenically different Vtp than the vaccine succumbed to infection.

Conclusions

We demonstrate that by having knowledge of the phenotypic changes made by B. hermsii as the spirochetes are acquired by ticks from infected mammals, an effective vaccine was developed that protected mice when challenged with infected ticks. However, the Vtp vaccines only protected mice from infection when challenged with that strain producing the identical Vtp. A vaccine containing multiple Vtp types may have promise as an oral vaccine for wild mammals if applied to geographic settings such as small islands where the mammal diversity is low and the Vtp types in the B. hermsii population are defined.

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.

Increased diversity of zoonotic pathogens and Borrelia burgdorferi strains in established versus incipient Ixodes scapularis populations across the Midwestern United States

The center of origin theory predicts that genetic diversity will be greatest near a specie's geographic origin because of the length of time for evolution. By corollary, diversity will decrease with distance from the origin; furthermore, invasion and colonization are frequently associated with founder effects that reduce genetic variation in incipient populations. The blacklegged tick, Ixodes scapularis, harbors a suite of zoonotic pathogens, and the geographic range of the tick is expanding in the upper Midwestern United States. Therefore, we posited that diversity of I. scapularis-borne pathogens across its Midwestern range should correlate with the rate of the range expansion of this tick as well as subsequent disease emergence. Analysis of 1565 adult I. scapularis ticks from 13 sites across five Midwestern states revealed that tick infection prevalence with multiple microbial agents (Borrelia burgdorferi, Borrelia miyamotoi, Babesia odocoilei, Babesia microti, and Anaplasma phagocytophilum), coinfections, and molecular genetic diversity of B. burgdorferi all were positively correlated with the duration of establishment of tick populations, and therefore generally support the center of origin - pathogen diversity hypothesis. The observed differences across the gradient of establishment, however, were not strong and were nuanced by the high frequency of coinfections in tick populations at both established and recently-invaded tick populations. These results suggest that the invasion of ticks and their associated pathogens likely involve multiple means of pathogen introduction, rather than the conventionally presented scenario whereby infected, invading ticks are solely responsible for introducing pathogens to naïve host populations.

Inter- and intra-specific pan-genomes of Borrelia burgdorferi sensu lato: genome stability and adaptive radiation

BACKGROUND

Lyme disease is caused by spirochete bacteria from the Borrelia burgdorferi sensu lato (B. burgdorferi s.l.) species complex. To reconstruct the evolution of B. burgdorferi s.l. and identify the genomic basis of its human virulence, we compared the genomes of 23 B. burgdorferi s.l. isolates from Europe and the United States, including B. burgdorferi sensu stricto (B. burgdorferi s.s., 14 isolates), B. afzelii (2), B. garinii (2), B. "bavariensis" (1), B. spielmanii (1), B. valaisiana (1), B. bissettii (1), and B. "finlandensis" (1).

RESULTS

Robust B. burgdorferi s.s. and B. burgdorferi s.l. phylogenies were obtained using genome-wide single-nucleotide polymorphisms, despite recombination. Phylogeny-based pan-genome analysis showed that the rate of gene acquisition was higher between species than within species, suggesting adaptive speciation. Strong positive natural selection drives the sequence evolution of lipoproteins, including chromosomally-encoded genes 0102 and 0404, cp26-encoded ospC and b08, and lp54-encoded dbpA, a07, a22, a33, a53, a65. Computer simulations predicted rapid adaptive radiation of genomic groups as population size increases.

CONCLUSIONS

Intra- and inter-specific pan-genome sizes of B. burgdorferi s.l. expand linearly with phylogenetic diversity. Yet gene-acquisition rates in B. burgdorferi s.l. are among the lowest in bacterial pathogens, resulting in high genome stability and few lineage-specific genes. Genome adaptation of B. burgdorferi s.l. is driven predominantly by copy-number and sequence variations of lipoprotein genes. New genomic groups are likely to emerge if the current trend of B. burgdorferi s.l. population expansion continues.

Multilocus sequence typing of Borrelia burgdorferi suggests existence of lineages with differential pathogenic properties in humans

The clinical manifestations of Lyme disease, caused by Borrelia burgdorferi, vary considerably in different patients, possibly due to infection by strains with varying pathogenicity. Both rRNA intergenic spacer and ospC typing methods have proven to be useful tools for categorizing B. burgdorferi strains that vary in their tendency to disseminate in humans. Neither method, however, is suitable for inferring intraspecific relationships among strains that are important for understanding the evolution of pathogenicity and the geographic spread of disease. In this study, multilocus sequence typing (MLST) was employed to investigate the population structure of B. burgdorferi recovered from human Lyme disease patients. A total of 146 clinical isolates from patients in New York and Wisconsin were divided into 53 sequence types (STs). A goeBURST analysis, that also included previously published STs from the northeastern and upper Midwestern US and adjoining areas of Canada, identified 11 major and 3 minor clonal complexes, as well as 14 singletons. The data revealed that patients from New York and Wisconsin were infected with two distinct, but genetically and phylogenetically closely related, populations of B. burgdorferi. Importantly, the data suggest the existence of B. burgdorferi lineages with differential capabilities for dissemination in humans. Interestingly, the data also indicate that MLST is better able to predict the outcome of localized or disseminated infection than is ospC typing.

Identification of Borrelia burgdorferi ospC genotypes in host tissue and feeding ticks by terminal restriction fragment length polymorphisms

We developed a high-throughput method based on terminal restriction fragment length polymorphisms (T-RFLP) to identify ospC genotypes from field-collected samples of Borrelia burgdorferi. We first validated the method by analyzing B. burgdorferi ospC previously identified by sequencing. We then analyzed and compared ospC genotypes detected from ear biopsy tissue from natural populations of the white-footed mouse, a major B. burgdorferi reservoir host species in the eastern United States, and larval ticks feeding on those individual mice. The T-RFLP method enabled us to distinguish all 17 ospC genotypes tested, as well as mixed samples containing more than one genotype. Analysis costs compare favorably to those of alternative ospC identification methods. The T-RFLP method will facilitate large-scale field studies to advance our understanding of genotype-specific transmission patterns.

Genetics of Borrelia burgdorferi

The spirochetes in the Borrelia burgdorferi sensu lato genospecies group cycle in nature between tick vectors and vertebrate hosts. The current assemblage of B. burgdorferi sensu lato, of which three species cause Lyme disease in humans, originated from a rapid species radiation that occurred near the origin of the clade. All of these species share a unique genome structure that is highly segmented and predominantly composed of linear replicons. One of the circular plasmids is a prophage that exists as several isoforms in each cell and can be transduced to other cells, likely contributing to an otherwise relatively anemic level of horizontal gene transfer, which nevertheless appears to be adequate to permit strong natural selection and adaptation in populations of B. burgdorferi. Although the molecular genetic toolbox is meager, several antibiotic-resistant mutants have been isolated, and the resistance alleles, as well as some exogenous genes, have been fashioned into markers to dissect gene function. Genetic studies have probed the role of the outer membrane lipoprotein OspC, which is maintained in nature by multiple niche polymorphisms and negative frequency-dependent selection. One of the most intriguing genetic systems in B. burgdorferi is vls recombination, which generates antigenic variation during infection of mammalian hosts.

Genome stability of Lyme disease spirochetes: comparative genomics of Borrelia burgdorferi plasmids

Lyme disease is the most common tick-borne human illness in North America. In order to understand the molecular pathogenesis, natural diversity, population structure and epizootic spread of the North American Lyme agent, Borrelia burgdorferi sensu stricto, a much better understanding of the natural diversity of its genome will be required. Towards this end we present a comparative analysis of the nucleotide sequences of the numerous plasmids of B. burgdorferi isolates B31, N40, JD1 and 297. These strains were chosen because they include the three most commonly studied laboratory strains, and because they represent different major genetic lineages and so are informative regarding the genetic diversity and evolution of this organism. A unique feature of Borrelia genomes is that they carry a large number of linear and circular plasmids, and this work shows that strains N40, JD1, 297 and B31 carry related but non-identical sets of 16, 20, 19 and 21 plasmids, respectively, that comprise 33–40% of their genomes. We deduce that there are at least 28 plasmid compatibility types among the four strains. The B. burgdorferi ∼900 Kbp linear chromosomes are evolutionarily exceptionally stable, except for a short ≤20 Kbp plasmid-like section at the right end. A few of the plasmids, including the linear lp54 and circular cp26, are also very stable. We show here that the other plasmids, especially the linear ones, are considerably more variable. Nearly all of the linear plasmids have undergone one or more substantial inter-plasmid rearrangements since their last common ancestor. In spite of these rearrangements and differences in plasmid contents, the overall gene complement of the different isolates has remained relatively constant.

Detection of established virulence genes and plasmids to differentiate Borrelia burgdorferi strains

Borrelia burgdorferi sensu stricto is the major causative agent of Lyme disease in the United States, while B. garinii and B. afzelii are more prevalent in Europe. The highly complex genome of B. burgdorferi is comprised of a linear chromosome and a large number of variably sized linear and circular plasmids. Many plasmids of this spirochete are unstable during its culture in vitro. Given that many of the B. burgdorferi virulence factors identified to date are plasmid encoded, spirochetal plasmid content determination is essential for genetic analysis of Lyme pathogenesis. Although PCR-based assays facilitate plasmid profiling of sequenced B. burgdorferi strains, a rapid genetic content determination strategy for nonsequenced strains has not yet been described. In this study, we combined pulsed-field gel electrophoresis (PFGE) and Southern hybridization for detection of genes encoding known virulence factors, ribosomal RNA gene spacer restriction fragment length polymorphism types (RSTs), ospC group determination, and sequencing of the variable dbpA and ospC genes. We show that two strains isolated from the same tick and both originally named N40 are in fact very distinct. Furthermore, we failed to detect bbk32, which encodes a fibronectin-binding adhesin, in one "N40" strain. Thus, two distinct strains that show different plasmid profiles, as determined by PFGE and PCR, were isolated from the same tick and vary in their ospC and dbpA sequences. However, both belong to group RST3B.

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.

BB0844, an RpoS-regulated protein, is dispensable for Borrelia burgdorferi infectivity and maintenance in the mouse-tick infectious cycle

The genome of Borrelia burgdorferi, the causative agent of Lyme disease, is comprised of a large linear chromosome and numerous smaller linear and circular plasmids. B. burgdorferi exhibits substantial genomic variation, and previous studies revealed genotype-specific variation at the right chromosomal telomere. A correlation has also been established between genotype and invasiveness. The correlation between chromosome length and genotype and between genotype and invasiveness suggested that a gene(s) at the right chromosome telomere may be required for virulence. Of particular interest was bb0844, an RpoS-regulated gene at the right telomere, the expression of which is induced when the spirochete undergoes adaptation to the mammalian host. The structure of the right chromosomal telomere was examined in 53 B. burgdorferi clinical isolates of various genotypes. Four distinct patterns were observed for bb0844: (i) chromosomal localization, (ii) plasmid localization, (iii) presence on both chromosome and plasmid, and (iv) complete absence. These patterns correlated with the B. burgdorferi genotype. On the basis of available sequence data, we propose a mechanism for the genomic rearrangements that accounts for the variability in bb0844 genomic localization. To further explore the role of BB0844 in the spirochete life cycle, a bb0844 deletion mutant was constructed by allelic exchange, and the viability of wild-type and bb0844 deletion mutants was examined in an experimental mouse-tick infection model. The bb0844 mutant was fully infectious in C3H/HeJ mice by either needle inoculation or tick transmission with B. burgdorferi-infected Ixodes scapularis larvae. Naïve larval ticks acquired both wild-type and mutant spirochetes with equal efficiency from B. burgdorferi-infected mice. The results demonstrate that BB0844 is not required for spirochete viability, pathogenicity, or maintenance in the tick vector or the mammalian host. At present, a defined role for BB0844 in B. burgdorferi cannot be ascertained.

Whole-genome sequences of thirteen isolates of Borrelia burgdorferi

Borrelia burgdorferi is a causative agent of Lyme disease in North America and Eurasia. The first complete genome sequence of B. burgdorferi strain 31, available for more than a decade, has assisted research on the pathogenesis of Lyme disease. Because a single genome sequence is not sufficient to understand the relationship between genotypic and geographic variation and disease phenotype, we determined the whole-genome sequences of 13 additional B. burgdorferi isolates that span the range of natural variation. These sequences should allow improved understanding of pathogenesis and provide a foundation for novel detection, diagnosis, and prevention strategies.

The genetic structure of Borrelia afzelii varies with geographic but not ecological sampling scale

The genetic structure of a pathogen is an important determinant of its potential rate of adaptation and can thereby influence the dynamics of host-parasite interactions. We investigated how the genetic structure of Borrelia afzelii varies with geographic and ecological sampling scale. Genetic structure was measured as the degree of linkage disequilibrium (LD) across three loci. To test for the effects of geographic and ecological scale, we calculated LD across or within populations 4-82 km apart and across or within different mammal host species. There was highly significant LD across populations and host species. However, there was also evidence for genome-wide recombination, and the LD largely resulted from epidemic spread of certain haplotypes, rather than lack of recombination. Interestingly, the degree of LD was higher in each population than in the sample as a whole, i.e. LD increased with decreasing geographic scale. In contrast, there was no effect of ecological sampling scale on LD. Strong LD may impede the rate of adaptive evolution. Our results suggest this effect might be particularly strong at a small geographic scale.

Geographic differences in genetic locus linkages for Borrelia burgdorferi

Borrelia burdorferi genotype in the northeastern United States is associated with Lyme borreliosis severity. Analysis of DNA sequences of the outer surface protein C gene and rrs-rrlA intergenic spacer from extracts of Ixodes spp. ticks in 3 US regions showed linkage disequilibrium between the 2 loci within a region but not consistently between regions.

Evolution and distribution of the ospC Gene, a transferable serotype determinant of Borrelia burgdorferi

Borrelia burgdorferi, an emerging bacterial pathogen, is maintained in nature by transmission from one vertebrate host to another by ticks. One of the few antigens against which mammals develop protective immunity is the highly polymorphic OspC protein, encoded by the ospC gene on the cp26 plasmid. Intragenic recombination among ospC genes is known, but the extent to which recombination extended beyond the ospC locus itself is undefined. We accessed and supplemented collections of DNA sequences of ospC and other loci from ticks in three U.S. regions (the Northeast, the Midwest, and northern California); a total of 839 ospC sequences were analyzed. Three overlapping but distinct populations of B. burgdorferi corresponded to the geographic regions. In addition, we sequenced 99 ospC flanking sequences from different lineages and compared the complete cp26 sequences of 11 strains as well as the cp26 bbb02 loci of 56 samples. Besides recombinations with traces limited to the ospC gene itself, there was evidence of lateral gene transfers that involved (i) part of the ospC gene and one of the two flanks or (ii) the entire ospC gene and different lengths of both flanks. Lateral gene transfers resulted in different linkages between the ospC gene and loci of the chromosome or other plasmids. By acquisition of the complete part or a large part of a novel ospC gene, an otherwise adapted strain would assume a new serotypic identity, thereby being comparatively fitter in an area with a high prevalence of immunity to existing OspC types.

The tick-borne zoonosis Lyme borreliosis is increasing in incidence and spreading geospatially in North America. Further understanding of the evolution and genetics of its cause, Borrelia burgdorferi, in its environments fosters progress toward ecologically based control efforts. By means of DNA sequencing of a large sample collection of the pathogen from across the United States, we studied the gene for the bacterium’s highly diverse OspC protein, protective immunity against which develops in animals. We found that the distributions and frequencies of types of OspC genes differed between populations of B. burgdorferi in the Northeast, the Midwest, and California. Over time, OspC genes were transferred between strains through recombinations involving the whole or parts of the gene and one or both flanks. Acquisitions of OspC genes that are novel for the region confer to recipients unique identities to host immune systems and, presumably, selective advantage when immunity to existing types is widespread among hosts.

Evolution of northeastern and midwestern Borrelia burgdorferi, United States

The per capita incidence of human Lyme disease in the northeastern United States is more than twice that in the Midwest. However, the prevalence of Borrelia burgdorferi, the bacterium that causes Lyme disease, in the tick vector is nearly identical in the 2 regions. The disparity in human Lyme disease incidence may result from a disparity in the human invasiveness of the bacteria in the Northeast and Midwest caused by fundamentally different evolutionary histories. B. burgdorferi populations in the Northeast and Midwest are geographically isolated, enabling evolutionary divergence in human invasiveness. However, we found that B. burgdorferi populations in the Northeast and Midwest shared a recent common ancestor, which suggests that substantial evolutionary divergence in human invasiveness has not occurred. We propose that differences in either animal ecology or human behavior are the root cause of the differences in human incidence between the 2 regions.

Identification of residues within ligand-binding domain 1 (LBD1) of the Borrelia burgdorferi OspC protein required for function in the mammalian environment

Borrelia burgdorferi outer surface protein C (ospC) is required for the establishment of infection in mammals. However, its precise function remains controversial. The biologically active form of OspC appears to be a homodimer. Alpha helix 1 and 1' of the apposing monomers form a solvent-accessible pocket at the dimeric interface that presents a putative ligand-binding domain (LBD1). Here we employ site-directed and allelic-exchange mutagenesis to test the hypothesis that LBD1 is a determinant of OspC function in the mammalian environment. Substitution of residues K60, E61 and E63 which line LBD1 resulted in the loss of infectivity or influenced dissemination. Analyses of the corresponding recombinant proteins demonstrated that the loss of function was not due to structural perturbation, impaired dimer formation or the loss of plasminogen binding. This study is the first to assess the involvement of individual residues and domains of OspC in its in vivo function. The data support the hypothesis that OspC interacts with a mammalian derived ligand that is critical for survival during early infection. These results shed new light on the structure-functions relationships of OspC and challenge existing hypotheses regarding OspC function in mammals.

Population structure of the lyme borreliosis spirochete Borrelia burgdorferi in the western black-legged tick (Ixodes pacificus) in Northern California

Factors potentially contributing to the lower incidence of Lyme borreliosis (LB) in the far-western than in the northeastern United States include tick host-seeking behavior resulting in fewer human tick encounters, lower densities of Borrelia burgdorferi-infected vector ticks in peridomestic environments, and genetic variation among B. burgdorferi spirochetes to which humans are exposed. We determined the population structure of B. burgdorferi in over 200 infected nymphs of the primary bridging vector to humans, Ixodes pacificus, collected in Mendocino County, CA. This was accomplished by sequence typing the spirochete lipoprotein ospC and the 16S-23S rRNA intergenic spacer (IGS). Thirteen ospC alleles belonging to 12 genotypes were found in California, and the two most abundant, ospC genotypes H3 and E3, have not been detected in ticks in the Northeast. The most prevalent ospC and IGS biallelic profile in the population, found in about 22% of ticks, was a new B. burgdorferi strain defined by ospC genotype H3. Eight of the most common ospC genotypes in the northeastern United States, including genotypes I and K that are associated with disseminated human infections, were absent in Mendocino County nymphs. ospC H3 was associated with hardwood-dominated habitats where western gray squirrels, the reservoir host, are commonly infected with LB spirochetes. The differences in B. burgdorferi population structure in California ticks compared to the Northeast emphasize the need for a greater understanding of the genetic diversity of spirochetes infecting California LB patients.

A new Borrelia species defined by multilocus sequence analysis of housekeeping genes

Analysis of Lyme borreliosis (LB) spirochetes, using a novel multilocus sequence analysis scheme, revealed that OspA serotype 4 strains (a rodent-associated ecotype) of Borrelia garinii were sufficiently genetically distinct from bird-associated B. garinii strains to deserve species status. We suggest that OspA serotype 4 strains be raised to species status and named Borrelia bavariensis sp. nov. The rooted phylogenetic trees provide novel insights into the evolutionary history of LB spirochetes.

Borrelia burgdorferi genotype predicts the capacity for hematogenous dissemination during early Lyme disease

BACKGROUND

Lyme disease, the most common tickborne disease in the United States, is caused exclusively by Borrelia burgdorferi sensu stricto in North America. The present study evaluated the genotypes of >400 clinical isolates of B. burgdorferi recovered from patients from suburban New York City with early Lyme disease associated with erythema migrans; it is the largest number of borrelial strains from North America ever to be investigated.

METHODS

Genotyping was performed by restriction fragment-length polymorphism polymerase chain reaction analysis of the 16S-23S ribosomal RNA spacer and reverse line blot analysis of the outer surface protein C gene (ospC). For some isolates, DNA sequence analysis was also performed.

RESULTS

The findings showed that the 16S-23S ribosomal spacer and ospC are in strong linkage disequilibrium. Most B. burgdorferi genotypes characterized by either typing method were capable of infecting and disseminating in patients. However, a distinct subset of just 4 of the 16 ospC genotypes identified were responsible for >80% of cases of early disseminated Lyme disease.

CONCLUSIONS

This study identified the B. burgdorferi genotypes that pose the greatest risk of causing hematogenous dissemination in humans. This information should be considered in the future development of diagnostic assays and vaccine preparations.

Profiling the humoral immune response to Borrelia burgdorferi infection with protein microarrays

To determine the cell envelope proteins of Borrelia burgdorferi recognized by immune sera of patients with late Lyme disease, we developed a Borrelia microarray containing proteins encoded by 90 cell envelope genes and their homologs described in the annotated genomic sequence of B. burgdorferi, strain B31. The protein microarray was used to profile the humoral immune response using sera from 13 patients with late Lyme disease and four normal controls. Although there was considerable heterogeneity in the individual sera responses, 25 of the cell envelope proteins were recognized by seven or more samples. Sera from non-infected individuals lacked reactivity against any of the proteins on the array. Among the most antigenic envelope proteins, BLAST search revealed little sequence homology to known microbial proteins from other species. The proteins that were highly seropositive included several members of the Erp gene families, BBA24 (decorin binding protein A (DbpA)) and members of the Borrelia gene family Pfam113 that code for the Mlp lipoprotein gene family. Several novel, uncharacterized B. burgdorferi antigens identified in this study were BBA14, BBG23, BB0108, BB0442 and BBQ03. The accurate diagnosis of Lyme disease depends on correlating objective clinical abnormalities with serological evidence of exposure to B. burgdorferi. A protein array of the envelope proteins of Borrelia burgdorferi may be very useful in specifically identifying patients with Lyme disease. This approach could contribute to a more rapid discovery of antigens not expressed in vitro that may be useful for the development of vaccine and diagnostics.

Variable exposure and immunological response to Lyme disease Borrelia among North Atlantic seabird species

Colonial seabirds often breed in large aggregations. These individuals can be exposed to parasitism by the tick Ixodes uriae, but little is known about the circulation of pathogens carried by this ectoparasite, including Lyme disease Borrelia. Here we investigated the prevalence of antibodies (Ab) against Borrelia burgdorferi sensu lato in seabird species sampled at eight locations across the North Atlantic. Using enzyme-linked immunosorbent assay tests, we found that the prevalence of anti-Borrelia Ab in adult seabirds was 39.6% on average (over 444 individuals), but that it varied among colonies and species. Common guillemots showed higher seroprevalence (77.1%+/-5.9) than black-legged kittiwakes (18.6%+/-6.7) and Atlantic puffins (22.6%+/-6.3). Immunoblot-banding patterns of positive individuals, reflecting the variability of Borrelia antigens against which Ab were produced, also differed among locations and species, and did not tightly match the prevalence of Borrelia phylogroups previously identified in ticks collected from the same host individuals. These results represent the first report of the widespread prevalence of Ab against Borrelia within an assemblage of seabird species and demonstrate that Borrelia is an integrated aspect in the interaction between seabirds and ticks. More detailed studies on the dynamics of Borrelia within and among seabird species at different spatial scales will now be required to better understand the implications of this interaction for seabird ecology and the epidemiology of Lyme disease.

Borrelia burgdorferi sensu stricto is clonal in patients with early Lyme borreliosis

Lyme borreliosis, the most commonly reported vector-borne disease in North America, is caused by the spirochete Borrelia burgdorferi. Given the extensive genetic polymorphism of B. burgdorferi, elucidation of the population genetic structure of the bacterium in clinical samples may be relevant for understanding disease pathogenesis and may have applicability for the development of diagnostic tests and vaccine preparations. In this investigation, the genetic polymorphism of the 16S-23S rRNA (rrs-rrlA) intergenic spacer and ospC was investigated at the sequence level in 127 clinical isolates obtained from patients with early Lyme borreliosis evaluated in suburban New York City. Sixteen distinct rrs-rrlA and 16 distinct ospC alleles were identified, representing virtually all of the genotypes previously found in questing Ixodes scapularis nymphs in this region. In addition, a new ospC group was identified in a single patient. The strong linkage observed between the chromosome-located rrs-rrlA and plasmid-borne ospC genes suggests a clonal structure of B. burgdorferi in these isolates, despite evidence of recombination at ospC.

Wide distribution of a high-virulence Borrelia burgdorferi clone in Europe and North America

The A and B clones of Borrelia burgdorferi sensu stricto, distinguished by outer surface protein C (ospC) gene sequences, are commonly associated with disseminated Lyme disease. To resolve phylogenetic relationships among isolates, we sequenced 68 isolates from Europe and North America at 1 chromosomal locus (16S-23S ribosomal RNA spacer) and 3 plasmid loci (ospC,dbpA, and BBD14). The ospC-A clone appeared to be highly prevalent on both continents, and isolates of this clone were uniform in DNA sequences, which suggests a recent trans-oceanic migration. The genetic homogeneity of ospC-A isolates was confirmed by sequences at 6 additional chromosomal housekeeping loci (gap, alr, glpA, xylB, ackA, and tgt). In contrast, the ospC-B group consists of genotypes distinct to each continent, indicating geographic isolation. We conclude that the ospC-A clone has dispersed rapidly and widely in the recent past. The spread of the ospC-A clone may have contributed, and likely continues to contribute, to the rise of Lyme disease incidence.

MLST of housekeeping genes captures geographic population structure and suggests a European origin of Borrelia burgdorferi

Lyme borreliosis, caused by the tick-borne bacterium Borrelia burgdorferi, has become the most common vector-borne disease in North America over the last three decades. To understand the dynamics of the epizootic spread and to predict the evolutionary trajectories of B. burgdorferi, accurate information on the population structure and the evolutionary relationships of the pathogen is crucial. We, therefore, developed a multilocus sequence typing (MLST) scheme for B. burgdorferi based on eight chromosomal housekeeping genes. We validated the MLST scheme on B. burgdorferi specimens from North America and Europe, comprising both cultured isolates and infected ticks. These data were compared with sequences for the commonly used genetic markers rrs-rrlA intergenic spacer (IGS) and the gene encoding the outer surface protein C (ospC). The study demonstrates that the concatenated sequences of the housekeeping genes of B. burgdorferi provide highly resolved phylogenetic signals and that the housekeeping genes evolve differently compared with the IGS locus and ospC. Using sequence data, the study reveals that North American and European populations of B. burgdorferi correspond to genetically distinct populations. Importantly, the MLST data suggest that B. burgdorferi originated in Europe rather than in North America as proposed previously.

Genomic sequencing of single microbial cells from environmental samples

Recently developed techniques allow genomic DNA sequencing from single microbial cells [Lasken RS: Single-cell genomic sequencing using multiple displacement amplification. Curr Opin Microbiol 2007, 10:510-516]. Here, we focus on research strategies for putting these methods into practice in the laboratory setting. An immediate consequence of single-cell sequencing is that it provides an alternative to culturing organisms as a prerequisite for genomic sequencing. The microgram amounts of DNA required as template are amplified from a single bacterium by a method called multiple displacement amplification (MDA) avoiding the need to grow cells. The ability to sequence DNA from individual cells will likely have an immense impact on microbiology considering the vast numbers of novel organisms, which have been inaccessible unless culture-independent methods could be used. However, special approaches have been necessary to work with amplified DNA. MDA may not recover the entire genome from the single copy present in most bacteria. Also, some sequence rearrangements can occur during the DNA amplification reaction. Over the past two years many research groups have begun to use MDA, and some practical approaches to single-cell sequencing have been developed. We review the consensus that is emerging on optimum methods, reliability of amplified template, and the proper interpretation of 'composite' genomes which result from the necessity of combining data from several single-cell MDA reactions in order to complete the assembly. Preferred laboratory methods are considered on the basis of experience at several large sequencing centers where >70% of genomes are now often recovered from single cells. Methods are reviewed for preparation of bacterial fractions from environmental samples, single-cell isolation, DNA amplification by MDA, and DNA sequencing.

OspC phylogenetic analyses support the feasibility of a broadly protective polyvalent chimeric Lyme disease vaccine

Using available Borrelia outer surface protein C (OspC) sequences, a phylogenetic analysis was undertaken to delineate the number of antigenic domains required for inclusion in a broadly protective, chimeric, OspC-based Lyme disease vaccine. The data indicate that approximately 34 would be required and that an OspC-based vaccinogen is feasible.