A population genetic study of Borrelia burgdorferi sensu stricto from eastern Long Island, New York, suggested frequency-dependent selection, gene flow and host adaptation

Characterization and genomic analysis of the Lyme disease spirochete bacteriophage ϕBB-1

Lyme disease is a tick-borne infection caused by the spirochete Borrelia (Borreliella) burgdorferi. Borrelia species have highly fragmented genomes composed of a linear chromosome and a constellation of linear and circular plasmids some of which are required throughout the enzootic cycle. Included in this plasmid repertoire by almost all Lyme disease spirochetes are the 32-kb circular plasmid cp32 prophages that are capable of lytic replication to produce infectious virions called ϕBB-1. While the B. burgdorferi genome contains evidence of horizontal transfer, the mechanisms of gene transfer between strains remain unclear. While we know that ϕBB-1 transduces cp32 and shuttle vector DNA during in vitro cultivation, the extent of ϕBB-1 DNA transfer is not clear. Herein, we use proteomics and long-read sequencing to further characterize ϕBB-1 virions. Our studies identified the cp32 pac region and revealed that ϕBB-1 packages linear cp32s via a headful mechanism with preferentially packaging of plasmids containing the cp32 pac region. Additionally, we find ϕBB-1 packages fragments of the linear chromosome and full-length plasmids including lp54, cp26, and others. Furthermore, sequencing of ϕBB-1 packaged DNA allowed us to resolve the covalently closed hairpin telomeres for the linear B. burgdorferi chromosome and most linear plasmids in strain CA-11.2A. Collectively, our results shed light on the biology of the ubiquitous ϕBB-1 phage and further implicates ϕBB-1 in the generalized transduction of diverse genes and the maintenance of genetic diversity in Lyme disease spirochetes.

Host adaptation drives genetic diversity in a vector-borne disease system

The range of hosts a pathogen can infect is a key trait, influencing human disease risk and reservoir host infection dynamics. Borrelia burgdorferi sensu stricto (Bb), an emerging zoonotic pathogen, causes Lyme disease and is widely considered a host generalist, commonly infecting mammals and birds. Yet the extent of intraspecific variation in Bb host breadth, its role in determining host competence, and potential implications for human infection remain unclear. We conducted a long-term study of Bb diversity, defined by the polymorphic ospC locus, across white-footed mice, passerine birds, and tick vectors, leveraging long-read amplicon sequencing. Our results reveal strong variation in host breadth across Bb genotypes, exposing a spectrum of genotype-specific host-adapted phenotypes. We found support for multiple niche polymorphism, maintaining Bb diversity in nature and little evidence of temporal shifts in genotype dominance, as would be expected under negative frequency-dependent selection. Passerine birds support the circulation of several human-invasive strains (HISs) in the local tick population and harbor greater Bb genotypic diversity compared with white-footed mice. Mouse-adapted Bb genotypes exhibited longer persistence in individual mice compared with nonadapted genotypes. Genotype communities infecting individual mice preferentially became dominated by mouse-adapted genotypes over time. We posit that intraspecific variation in Bb host breadth and adaptation helps maintain overall species fitness in response to transmission by a generalist vector.

Competition between strains of Borrelia afzelii in the host tissues and consequences for transmission to ticks

Pathogen species often consist of genetically distinct strains, which can establish mixed infections or coinfections in the host. In coinfections, interactions between pathogen strains can have important consequences for their transmission success. We used the tick-borne bacterium Borrelia afzelii, which is the most common cause of Lyme disease in Europe, as a model multi-strain pathogen to investigate the relationship between coinfection, competition between strains, and strain-specific transmission success. Mus musculus mice were infected with one or two strains of B. afzelii, strain transmission success was measured by feeding ticks on mice, and the distribution of each strain in six different mouse organs and the ticks was measured using qPCR. Coinfection and competition reduced the tissue infection prevalence of both strains and changed their bacterial abundance in some tissues. Coinfection and competition also reduced the transmission success of the B. afzelii strains from the infected hosts to feeding ticks. The ability of the B. afzelii strains to establish infection in the host tissues was strongly correlated with their transmission success to the tick vector. Our study demonstrates that coinfection and competition between pathogen strains inside the host tissues can have major consequences for their transmission success.

Evidence for an effect of landscape connectivity on Borrelia burgdorferi sensu stricto dispersion in a zone of range expansion

In North America, different strains of the Lyme disease-causing bacterium Borrelia burgdorferi sensu stricto cluster into phylogenetic groups that are associated with different levels of pathogenicity and, for some, specific rodent reservoir hosts. Here we explore whether landscape connectivity, by impacting host dispersal, influences B. burgdorferi s.s. spread patterns. This question is central to modelling spatial patterns of the spread of Lyme disease risk in the zone of northward range-expansion of B. burgdorferi s.s. in southeastern Canada where the study was conducted. We used multi-locus sequence typing (MLST) to characterise B. burgdorferi s.s. in positive ticks collected at 13 sites in southern Quebec, Canada during the early stages of B. burgdorferi s.s. invasion. We used mixed effects logistic regression to investigate whether landscape connectivity (probability of connectivity; PC) affected the probability that samples collected at different sites were of the same strain (MLST sequence type: ST). PC was calculated from a habitat map based on high spatial resolution (15 m) Landsat 8 imagery to identify woodland habitat that are preferred by rodent hosts of B. burgdorferi s.s. There was a significant positive association between the likelihood that two samples were of the same ST and PC, when PC values were grouped into three categories of low, medium and high. When analysing data for individual STs, samples at different sites were significantly more likely to be the same when PC was higher for the rodent-associated ST1. These findings support the hypothesis that dispersion trajectories of B. burgdorferi s.s. in general, and some rodent-associated strains in particular, are at least partly determined by landscape connectivity. This may suggest that dispersion of B. burgdorferi s.s. is more common by terrestrial mammal hosts (which would likely disperse according to landscape connectivity) than by birds, the dispersal of which is likely less constrained by landscape. This study suggests that accounting for landscape connectivity may improve model-based predictions of spatial spread patterns of B. burgdorferi s.s. The findings are consistent with possible past dispersal patterns of B. burgdorferi s.s. as determined by phylogeographic studies.

Primordial origin and diversification of plasmids in Lyme disease agent bacteria

BACKGROUND

With approximately one-third of their genomes consisting of linear and circular plasmids, the Lyme disease agent cluster of species has the most complex genomes among known bacteria. We report here a comparative analysis of plasmids in eleven Borreliella (also known as Borrelia burgdorferi sensu lato) species.

RESULTS

We sequenced the complete genomes of two B. afzelii, two B. garinii, and individual B. spielmanii, B. bissettiae, B. valaisiana and B. finlandensis isolates. These individual isolates carry between seven and sixteen plasmids, and together harbor 99 plasmids. We report here a comparative analysis of these plasmids, along with 70 additional Borreliella plasmids available in the public sequence databases. We identify only one new putative plasmid compatibility type (the 30th) among these 169 plasmid sequences, suggesting that all or nearly all such types have now been discovered. We find that the linear plasmids in the non-B. burgdorferi species have undergone the same kinds of apparently random, chaotic rearrangements mediated by non-homologous recombination that we previously discovered in B. burgdorferi. These rearrangements occurred independently in the different species lineages, and they, along with an expanded chromosomal phylogeny reported here, allow the identification of several whole plasmid transfer events among these species. Phylogenetic analyses of the plasmid partition genes show that a majority of the plasmid compatibility types arose early, most likely before separation of the Lyme agent Borreliella and relapsing fever Borrelia clades, and this, with occasional cross species plasmid transfers, has resulted in few if any species-specific or geographic region-specific Borreliella plasmid types.

CONCLUSIONS

The primordial origin and persistent maintenance of the Borreliella plasmid types support their functional indispensability as well as evolutionary roles in facilitating genome diversity. The improved resolution of Borreliella plasmid phylogeny based on conserved partition-gene clusters will lead to better determination of gene orthology which is essential for prediction of biological function, and it will provide a basis for inferring detailed evolutionary mechanisms of Borreliella genomic variability including homologous gene and plasmid exchanges as well as non-homologous rearrangements.

Negative Frequency-Dependent Selection Is Frequently Confounding

Persistent genetic variation within populations presents an evolutionary problem, as natural selection and genetic drift tend to erode genetic diversity. Models of balancing selection were developed to account for the maintenance of genetic variation observed in natural populations. Negative frequency-dependent selection is a powerful type of balancing selection that maintains many natural polymorphisms, but it is also commonly misinterpreted. This review aims to clarify the processes underlying negative frequency-dependent selection, describe classes of polymorphisms that can and cannot result from these processes, and discuss the empirical data needed to accurately identify processes that generate or maintain diversity in nature. Finally, the importance of accurately describing the processes affecting genetic diversity within populations as it relates to research progress is considered.

Fitness estimates from experimental infections predict the long-term strain structure of a vector-borne pathogen in the field

The populations of many pathogen species consist of a collection of common and rare strains but the factors underlying this strain-specific variation in frequency are often unknown. Understanding frequency variation among strains is particularly challenging for vector-borne pathogens where the strain-specific fitness depends on the performance in both the vertebrate host and the arthropod vector. Two sympatric multiple-strain tick-borne pathogens, Borrelia afzelii and B. garinii, that use the same tick vector, Ixodes ricinus, but different vertebrate hosts were studied. 454-sequencing of the polymorphic ospC gene was used to characterize the community of Borrelia strains in a local population of I. ricinus ticks over a period of 11 years. Estimates of the reproduction number (R0), a measure of fitness, were obtained for six strains of B. afzelii from a previous laboratory study. There was substantial variation in prevalence among strains and some strains were consistently common whereas other strains were consistently rare. In B. afzelii, the strain-specific estimates of R0 in laboratory mice explained over 70% of the variation in the prevalences of the strains in our local population of ticks. Our study shows that laboratory estimates of fitness can predict the community structure of multiple-strain pathogens in the field.

Co-feeding transmission facilitates strain coexistence in Borrelia burgdorferi, the Lyme disease agent

Coexistence of multiple tick-borne pathogens or strains is common in natural hosts and can be facilitated by resource partitioning of the host species, within-host localization, or by different transmission pathways. Most vector-borne pathogens are transmitted horizontally via systemic host infection, but transmission may occur in the absence of systemic infection between two vectors feeding in close proximity, enabling pathogens to minimize competition and escape the host immune response. In a laboratory study, we demonstrated that co-feeding transmission can occur for a rapidly-cleared strain of Borrelia burgdorferi, the Lyme disease agent, between two stages of the tick vector Ixodes scapularis while feeding on their dominant host, Peromyscus leucopus. In contrast, infections rapidly became systemic for the persistently infecting strain. In a field study, we assessed opportunities for co-feeding transmission by measuring co-occurrence of two tick stages on ears of small mammals over two years at multiple sites. Finally, in a modeling study, we assessed the importance of co-feeding on R₀, the basic reproductive number. The model indicated that co-feeding increases the fitness of rapidly-cleared strains in regions with synchronous immature tick feeding. Our results are consistent with increased diversity of B. burgdorferi in areas of higher synchrony in immature feeding – such as the midwestern United States. A higher relative proportion of rapidly-cleared strains, which are less human pathogenic, would also explain lower Lyme disease incidence in this region. Finally, if co-feeding transmission also occurs on refractory hosts, it may facilitate the emergence and persistence of new pathogens with a more limited host range.

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.

Public health impact of strain specific immunity to Borrelia burgdorferi

BACKGROUND

Lyme disease, caused by Borrelia burgdorferi, is the most common tick-borne infection in the United States. Although humans can be infected by at least 16 different strains of B. burgdorferi, the overwhelming majority of infections are due to only four strains. It was recently demonstrated that patients who are treated for early Lyme disease develop immunity to the specific strain of B. burgdorferi that caused their infection. The aim of this study is to estimate the reduction in cases of Lyme disease in the United States that may occur as a result of type specific immunity.

METHODS

The analysis was performed based on three analytical models that assessed the effects of type specific immunity. Observational data on the frequency with which different B. burgdorferi strains cause human infection in culture-confirmed patients with an initial episode of erythema migrans diagnosed between 1991 and 2005 in the Northeastern United States were used in the analyses.

RESULTS

Assuming a reinfection rate of 3 % and a total incidence of Lyme disease per year of 300,000, the estimated number of averted cases of Lyme disease per year ranges from 319 to 2378 depending on the duration of type specific immunity and the model used.

CONCLUSION

Given the assumptions of the analyses, this analysis suggests that type specific immunity is likely to have public health significance in the United States.

Evolution and population genomics of the Lyme borreliosis pathogen, Borrelia burgdorferi

Population genomic studies have the potential to address many unresolved questions about microbial pathogens by facilitating the identification of genes underlying ecologically important traits, such as novel virulence factors and adaptations to humans or other host species. Additionally, this framework improves estimations of population demography and evolutionary history to accurately reconstruct recent epidemics and identify the molecular and environmental factors that resulted in the outbreak. The Lyme disease bacterium, Borrelia burgdorferi, exemplifies the power and promise of the application of population genomics to microbial pathogens. We discuss here the future of evolutionary studies in B. burgdorferi, focusing on the primary evolutionary forces of horizontal gene transfer, natural selection, and migration, as investigations transition from analyses of single genes to genomes.

Co-feeding transmission in Lyme disease pathogens

This review examines the phenomenon of co-feeding transmission in tick-borne pathogens. This mode of transmission is critical for the epidemiology of several tick-borne viruses but its importance for Borrelia burgdorferi sensu lato, the causative agents of Lyme borreliosis, is still controversial. The molecular mechanisms and ecological factors that facilitate co-feeding transmission are therefore examined with particular emphasis on Borrelia pathogens. Comparison of climate, tick ecology and experimental infection work suggests that co-feeding transmission is more important in European than North American systems of Lyme borreliosis, which potentially explains why this topic has gained more traction in the former continent than the latter. While new theory shows that co-feeding transmission makes a modest contribution to Borrelia fitness, recent experimental work has revealed new ecological contexts where natural selection might favour co-feeding transmission. In particular, co-feeding transmission might confer a fitness advantage in the Darwinian competition among strains in mixed infections. Future studies should investigate the ecological conditions that favour the evolution of this fascinating mode of transmission in tick-borne pathogens.

Occurrence and transmission efficiencies of Borrelia burgdorferi ospC types in avian and mammalian wildlife

Borrelia burgdorferi s.s., the bacterium that causes Lyme disease in North America, circulates among a suite of vertebrate hosts and their tick vector. The bacterium can be differentiated at the outer surface protein C (ospC) locus into 25 genotypes. Wildlife hosts can be infected with a suite of ospC types but knowledge on the transmission efficiencies of these naturally infected hosts to ticks is still lacking. To evaluate the occupancy and detection of ospC types in wildlife hosts, we adapted a likelihood-based species patch occupancy model to test for the occurrence probabilities (ψ - "occupancy") and transmission efficiencies (ε - "detection") of each ospC type. We detected differences in ospC occurrence and transmission efficiencies from the null models with HIS (human invasive strains) types A and K having the highest occurrence estimates, but both HIS and non-HIS types having high transmission efficiencies. We also examined ospC frequency patterns with respect to strains known to be invasive in humans across the host species and phylogenetic groups. We found that shrews and to a lesser extent, birds, were important host groups supporting relatively greater frequencies of HIS to non-HIS types. This novel method of simultaneously assessing occurrence and transmission of ospC types provides a powerful tool in assessing disease risk at the genotypic level in naturally infected wildlife hosts and offers the opportunity to examine disease risk at the community level.

Lyme disease risk not amplified in a species-poor vertebrate community: similar Borrelia burgdorferi tick infection prevalence and OspC genotype frequencies

The effect of biodiversity declines on human health is currently debated, but empirical assessments are lacking. Lyme disease provides a model system to assess relationships between biodiversity and human disease because the etiologic agent, Borrelia burgdorferi, is transmitted in the United States by the generalist black-legged tick (Ixodes scapularis) among a wide range of mammalian and avian hosts. The 'dilution effect' hypothesis predicts that species-poor host communities dominated by white-footed mice (Peromyscus leucopus) will pose the greatest human risk because P. leucopus infects the largest numbers of ticks, resulting in higher human exposure to infected I. scapularis ticks. P. leucopus-dominated communities are also expected to maintain a higher frequency of those B. burgdorferi outer surface protein C (ospC) genotypes that this host species more efficiently transmits ('multiple niche polymorphism' hypothesis). Because some of these genotypes are human invasive, an additive increase in human disease risk is expected in species-poor settings. We assessed these theoretical predictions by comparing I. scapularis nymphal infection prevalence, density of infected nymphs and B. burgdorferi genotype diversity at sites on Block Island, RI, where P. leucopus dominates the mammalian host community, to species-diverse sites in northeastern Connecticut. We found no support for the dilution effect hypothesis; B. burgdorferi nymphal infection prevalence was similar between island and mainland and the density of B. burgdorferi infected nymphs was higher on the mainland, contrary to what is predicted by the dilution effect hypothesis. Evidence for the multiple niche polymorphism hypothesis was mixed: there was lower ospC genotype diversity at island than mainland sites, but no overrepresentation of genotypes with higher fitness in P. leucopus or that are more invasive in humans. We conclude that other mechanisms explain similar nymphal infection prevalence in both communities and that high ospC genotype diversity can be maintained in both species-poor and species-rich communities.

Comparative population genomics of the Borrelia burgdorferi species complex reveals high degree of genetic isolation among species and underscores benefits and constraints to studying intra-specific epidemiological processes

Lyme borreliosis, one of the most frequently contracted zoonotic diseases in the Northern Hemisphere, is caused by bacteria belonging to different genetic groups within the Borrelia burgdorferi species complex, which are transmitted by ticks among various wildlife reservoirs, such as small mammals and birds. These features make the Borrelia burgdorferi species complex an attractive biological model that can be used to study the diversification and the epidemiology of endemic bacterial pathogens. We investigated the potential of population genomic approaches to study these processes. Sixty-three strains belonging to three species within the Borrelia burgdorferi complex were isolated from questing ticks in Alsace (France), a region where Lyme disease is highly endemic. We first aimed to characterize the degree of genetic isolation among the species sampled. Phylogenetic and coalescent-based analyses revealed clear delineations: there was a ∼50 fold difference between intra-specific and inter-specific recombination rates. We then investigated whether the population genomic data contained information of epidemiological relevance. In phylogenies inferred using most of the genome, conspecific strains did not cluster in clades. These results raise questions about the relevance of different strategies when investigating pathogen epidemiology. For instance, here, both classical analytic approaches and phylodynamic simulations suggested that population sizes and migration rates were higher in B. garinii populations, which are normally associated with birds, than in B. burgdorferi s.s. populations. The phylogenetic analyses of the infection-related ospC gene and its flanking region provided additional support for this finding. Traces of recombination among the B. burgdorferi s.s. lineages and lineages associated with small mammals were found, suggesting that they shared the same hosts. Altogether, these results provide baseline evidence that can be used to formulate hypotheses regarding the host range of B. burgdorferi lineages based on population genomic data.

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.

Multiple-strain infections of Borrelia afzelii: a role for within-host interactions in the maintenance of antigenic diversity?

Genetically diverse infections are common but little is known about what effects coinfecting strains have on each other in natural host-parasite systems. To explore the nature and consequences of interactions in the wild, we studied the tick-transmitted bacterium Borrelia afzelii in one of its main reservoir hosts, the bank vole Myodes glareolus. We measured overall infection intensity with quantitative polymerase chain reaction (PCR) and resolved the composition of multiple infections using strain-specific PCR assays targeting the ospC gene (which encodes an immunodominant surface protein). We found seven different strains, as defined by ospC genotype. There was little evidence for interactions affecting infection intensities, but strains were highly aggregated (i.e., there were more multiple infections than expected from random co-occurrence). Moreover, there was a positive correlation between the difference at the amino acid level between two OspC types and their degree of association. Overall, the observed patterns suggest that facilitation is more important than competition in this system and that more diverse infections have an advantage in establishing and/or maintaining infection. We propose that this advantage is one of the factors that favors antigenic diversity.

Detection of Borrelia burgdorferi sensu stricto ospC alleles associated with human lyme borreliosis worldwide in non-human-biting tick Ixodes affinis and rodent hosts in Southeastern United States

Comparative analysis of ospC genes from 127 Borrelia burgdorferi sensu stricto strains collected in European and North American regions where Lyme disease is endemic and where it is not endemic revealed a close relatedness of geographically distinct populations. ospC alleles A, B, and L were detected on both continents in vectors and hosts, including humans. Six ospC alleles, A, B, L, Q, R, and V, were prevalent in Europe; 4 of them were detected in samples of human origin. Ten ospC alleles, A, B, D, E3, F, G, H, H3, I3, and M, were identified in the far-western United States. Four ospC alleles, B, G, H, and L, were abundant in the southeastern United States. Here we present the first expanded analysis of ospC alleles of B. burgdorferi strains from the southeastern United States with respect to their relatedness to strains from other North American and European localities. We demonstrate that ospC genotypes commonly associated with human Lyme disease in European and North American regions where the disease is endemic were detected in B. burgdorferi strains isolated from the non-human-biting tick Ixodes affinis and rodent hosts in the southeastern United States. We discovered that some ospC alleles previously known only from Europe are widely distributed in the southeastern United States, a finding that confirms the hypothesis of transoceanic migration of Borrelia species.

Experimental infections of the reservoir species Peromyscus leucopus with diverse strains of Borrelia burgdorferi, a Lyme disease agent

The rodent Peromyscus leucopus is a major natural reservoir for the Lyme disease agent Borrelia burgdorferi and a host for its vector Ixodes scapularis. At various locations in northeastern United States 10 to 15 B. burgdorferi strains coexist at different prevalences in tick populations. We asked whether representative strains of high or low prevalence differed in their infections of P. leucopus. After 5 weeks of experimental infection of groups with each of 6 isolates, distributions and burdens of bacteria in tissues were measured by quantitative PCR, and antibodies to B. burgdorferi were evaluated by immunoblotting and protein microarray. All groups of animals were infected in their joints, ears, tails, and hearts, but overall spirochete burdens were lower in animals infected with low-prevalence strains. Animals were similar regardless of the infecting isolate in their levels of antibodies to whole cells, FlaB, BmpA, and DbpB proteins, and the conserved N-terminal region of the serotype-defining OspC proteins. But there were strain-specific antibody responses to full-length OspC and to plasmid-encoded VlsE, BBK07, and BBK12 proteins. Sequencing of additional VlsE genes revealed substantial diversity within some pairs of strains but near-identical sequences within other pairs, which otherwise differed in their ospC alleles. The presence or absence of full-length bbk07 and bbk12 genes accounted for the differences in antibody responses. We propose that for B. burgdorferi, there is selection in reservoir species for (i) sequence diversity, as for OspC and VlsE, and (ii) the presence or absence of polymorphisms, as for BBK07 and BBK12.

Humans are dead-end hosts for Borrelia agents of Lyme disease (LD), and, thus, irrelevant for the pathogens’ maintenance. Many reports of human cases and laboratory mouse infections exist, but less is known about infection and immunity in natural reservoirs, such as the rodent Peromyscus leucopus. We observed that high- and low-prevalence strains of Borrelia burgdorferi were capable of infecting P. leucopus but elicited different patterns of antibody responses. Antibody reactivities to the VlsE protein were as type-specific as previously characterized reactivities to serotype-defining OspC proteins. In addition, the low-prevalence strains lacked full-length genes for two proteins that (i) are encoded by a virulence-associated plasmid in some high-prevalence strains and (ii) LD patients and field-captured rodents commonly have antibodies to. Immune selection against these genes may have led to null phenotype lineages that can infect otherwise immune hosts but at the cost of reduced fitness and lower prevalence.

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.

Pervasive recombination and sympatric genome diversification driven by frequency-dependent selection in Borrelia burgdorferi, the Lyme disease bacterium

How genomic diversity within bacterial populations originates and is maintained in the presence of frequent recombination is a central problem in understanding bacterial evolution. Natural populations of Borrelia burgdorferi, the bacterial agent of Lyme disease, consist of diverse genomic groups co-infecting single individual vertebrate hosts and tick vectors. To understand mechanisms of sympatric genome differentiation in B. burgdorferi, we sequenced and compared 23 genomes representing major genomic groups in North America and Europe. Linkage analysis of >13,500 single-nucleotide polymorphisms revealed pervasive horizontal DNA exchanges. Although three times more frequent than point mutation, recombination is localized and weakly affects genome-wide linkage disequilibrium. We show by computer simulations that, while enhancing population fitness, recombination constrains neutral and adaptive divergence among sympatric genomes through periodic selective sweeps. In contrast, simulations of frequency-dependent selection with recombination produced the observed pattern of a large number of sympatric genomic groups associated with major sequence variations at the selected locus. We conclude that negative frequency-dependent selection targeting a small number of surface-antigen loci (ospC in particular) sufficiently explains the maintenance of sympatric genome diversity in B. burgdorferi without adaptive divergence. We suggest that pervasive recombination makes it less likely for local B. burgdorferi genomic groups to achieve host specialization. B. burgdorferi genomic groups in the northeastern United States are thus best viewed as constituting a single bacterial species, whose generalist nature is a key to its rapid spread and human virulence.

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.

Low-density microarrays for the detection of Borrelia burgdorferi s.s. (the Lyme disease spirochete) in nymphal Ixodes scapularis

Lyme disease is the most common tick-borne disease in Europe and North America. In the hyperendemic Lyme disease regions of the eastern United States, nymphal Ixodes scapularis are the principal ticks transmitting the Lyme disease spirochete, Borrelia burgdorferi sensu stricto (s.s.). Approximately 25% of questing nymphs in endemic regions are infected with spirochetes. High throughput-sensitive and specific methods for testing nymphal I. scapularis for infection with B. burgdorferi are clearly needed. In the current study, we evaluated whether low-density microarrays could be adapted for the rapid and accurate detection and characterization of spirochetes in nymphal I. scapularis. Three different microarray platforms were developed and tested for the detection of spirochetes in ticks. They could both detect and differentiate different Borrelia genospecies, in one case detecting as few as a single copy of Borrelia DNA.

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.

Comprehensive seroprofiling of sixteen B. burgdorferi OspC: implications for Lyme disease diagnostics design

Early diagnosis of Lyme disease (LD) is critical to successful treatment. However, current serodiagnostic tests do not reliably detect antibodies during early infection. OspC induces a potent early immune response and is also one of the most diverse proteins in the Borrelia proteome. Yet, at least 70% of the amino acid sequence is conserved among all 21 known OspC types. We performed a series of comprehensive seroprofiling studies to select the OspC types that have the most cross-reactive immunodominant epitopes. We found that proteins belonging to seven OspC types detect antibodies from all three infected host species regardless of the OspC genotype of the infecting strain. Although no one OspC type identifies all seropositive human samples, combinations of as few as two OspC proteins identified all patients that had anti-OspC antibodies.

Complex population structure of Lyme borreliosis group spirochete Borrelia garinii in subarctic Eurasia

Borrelia garinii, a causative agent of Lyme borreliosis in Europe and Asia, is naturally maintained in marine and terrestrial enzootic cycles, which primarily involve birds, including seabirds and migratory passerines. These bird groups associate with, correspondingly, Ixodes uriae and Ixodes ricinus ticks, of which the latter species may bite and transmit the infection to humans. Studies of the overlap between these two natural cycles of B. garinii have been limited, in part due to the absence of representative collections of this spirochete's samples, as well as of the lack of reliable measure of the genetic heterogeneity of its strains. As a prerequisite for understanding the epidemiological correlates of the complex maintenance of B. garinii, the present study sought to assess the diversity and phylogenetic relationships of this species' strains from its natural hosts and patients with Lyme borreliosis from subarctic Eurasia. We used sequence typing of the partial rrs-rrl intergenic spacer (IGS) of archived and prospective samples of B. garinii from I. uriae ticks collected predominantly on Commander Islands in North Pacific, as well as on the islands in northern Sweden and arctic Norway. We also typed B. garinii samples from patients with Lyme borreliosis and I. ricinus ticks infesting migratory birds in southern Sweden, or found questing in selected sites on the islands in the Baltic Sea and Lithuania. Fifty-two (68%) of 77 B. garinii samples representing wide geographical range and associated with I. ricinus and infection of humans contributed 12 (60%) of total 20 identified IGS variants. In contrast, the remaining 25 (32%) samples recovered from I. uriae ticks from a few islands accounted for as many as 10 (50%) IGS types, suggesting greater local diversity of B. garinii maintained by seabirds and their ticks. Two IGS variants of the spirochete in common for both tick species were found in I. ricinus larvae from migratory birds, an indication that B. garinii strains are exchanged between different ecological niches. Notably, B. garinii variants associated with I. uriae ticks were found in each of the six clusters, representing two phylogenetic lineages of this species identified among the studied samples. Our findings suggest that B. garinii in subarctic Eurasia comprises two partially overlapping populations with different levels of genetic heterogeneity, presumably, due to distinctive selective pressures on the spirochete in its marine and terrestrial enzootic cycles.

Reviewing molecular adaptations of Lyme borreliosis spirochetes in the context of reproductive fitness in natural transmission cycles

Lyme borreliosis (LB) is caused by a group of pathogenic spirochetes – most often Borrelia burgdorferi, B. afzelii, and B. garinii – that are vectored by hard ticks in the Ixodes ricinus-persulcatus complex, which feed on a variety of mammals, birds, and lizards. Although LB is one of the best-studied vector-borne zoonoses, the annual incidence in North America and Europe leads other vector-borne diseases and continues to increase. What factors make the LB system so successful, and how can researchers hope to reduce disease risk – either through vaccinating humans or reducing the risk of contacting infected ticks in nature? Discoveries of molecular interactions involved in the transmission of LB spirochetes have accelerated recently, revealing complex interactions among the spirochete-tick-vertebrate triad. These interactions involve multiple, and often redundant, pathways that reflect the evolution of general and specific mechanisms by which the spirochetes survive and reproduce. Previous reviews have focused on the molecular interactions or population biology of the system. Here molecular interactions among the LB spirochete, its vector, and vertebrate hosts are reviewed in the context of natural maintenance cycles, which represent the ecological and evolutionary contexts that shape these interactions. This holistic system approach may help researchers develop additional testable hypotheses about transmission processes, interpret laboratory results, and guide development of future LB control measures and management.

Fine-scale phylogeographic structure of Borrelia lusitaniae revealed by multilocus sequence typing

Borrelia lusitaniae is an Old World species of the Lyme borreliosis (LB) group of tick-borne spirochetes and prevails mainly in countries around the Mediterranean Basin. Lizards of the family Lacertidae have been identified as reservoir hosts of B. lusitaniae. These reptiles are highly structured geographically, indicating limited migration. In order to examine whether host geographic structure shapes the evolution and epidemiology of B. lusitaniae, we analyzed the phylogeographic population structure of this tick-borne bacterium using a recently developed multilocus sequence typing (MLST) scheme based on chromosomal housekeeping genes. A total of 2,099 questing nymphal and adult Ixodes ricinus ticks were collected in two climatically different regions of Portugal, being approximately 130 km apart. All ticks were screened for spirochetes by direct PCR. Attempts to isolate strains yielded 16 cultures of B. lusitaniae in total. Uncontaminated cultures as well as infected ticks were included in this study. The results using MLST show that the regional B. lusitaniae populations constitute genetically distinct populations. In contrast, no clear phylogeographic signals were detected in sequences of the commonly used molecular markers ospA and ospC. The pronounced population structure of B. lusitaniae over a short geographic distance as captured by MLST of the housekeeping genes suggests that the migration rates of B. lusitaniae are rather low, most likely because the distribution of mediterranean lizard populations is highly parapatric. The study underlines the importance of vertebrate hosts in the geographic spread of tick-borne microparasites.

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.

Genetic diversity of Borrelia burgdorferi sensu stricto in Peromyscus leucopus, the primary reservoir of Lyme disease in a region of endemicity in southern Maryland

In the north central and northeastern United States, Borrelia burgdorferi sensu stricto, the etiologic agent of Lyme disease (LD), is maintained in an enzootic cycle between the vector, Ixodes scapularis, and the primary reservoir host, Peromyscus leucopus. Genetic diversity of the pathogen based on sequencing of two plasmid-located genes, those for outer surface protein A (ospA) and outer surface protein C (ospC), has been examined in both tick and human specimens at local, regional, and worldwide population scales. Additionally, previous studies have only been conducted with tick or human specimens at the local population level in areas with high LD transmission rates. This study examined the genetic diversity of circulating borreliae in the reservoir population from a large region of the western coastal plains of southern Maryland, where moderate numbers of human LD cases are reported. Six ospA mobility classes, including two that were not previously described, and eight ospC groups were found among the P. leucopus samples. Twenty-five percent of all specimens were infected with more than one ospA or ospC variant. The frequency distribution of variants was homogeneous, both locally and spatially. The spirochete diversity found in Maryland was not as high as that observed among northern tick populations, yet similar genotypes were observed in both populations. These results also show that mice are important for maintaining Borrelia variants, even rare variants, and that reservoir populations should therefore be considered when assessing the diversity of B. burgdorferi.

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

Clinical and tick isolates of Borrelia burgdorferi sensu stricto, the bacterial agent of Lyme disease, from the northeastern United States were sequenced at 12 loci located on the main chromosome and 7 plasmids (lp54, cp26, cp9, lp17, lp25, lp28-2, and lp38). The outer surface protein C gene (ospC) showed the highest number (12) of major alleles (defined as alleles differing by 5% or more in nucleotide sequence) while other ORFs had only two to four major alleles. A non-recombining chromosomal marker, the rrs-rrlA ribosomal RNA spacer, was used to infer the intraspecific phylogeny among these B. burgdorferi isolates. We were thus able to analyze the multilocus genotypes in the context of a B. burgdorferi intraspecific phylogeny. Except for ospC, sequence variations at plasmid-borne loci showed broad inconsistency with the intraspecific phylogeny, supporting DNA exchanges mediated by plasmid transfers. The multilocus linkage frequently observed in B. burgdorferi populations is due more likely to a "founder effect" than to a lack of recombination. The exceptional phylogenetic consistency of ospC, in conjunction with its selectively maintained high intraspecific diversity, suggested a dominant role ospC plays in the initiation and maintenance of adaptive differentiation in B. burgdorferi.

Structure-based Design of a Second-generation Lyme Disease Vaccine Based on a C-terminal Fragment of Borrelia burgdorferi OspA

Here, we describe a structure-based approach to reduce the size of an antigen protein for a subunit vaccine. Our method consists of (i) determining the three-dimensional structure of an antigen, (ii) identifying protective epitopes, (iii) generation of an antigen fragment that contains the protective epitope, and (iv) rational design to compensate for destabilization caused by truncation. Using this approach we have successfully developed a second-generation Lyme disease vaccine. Outer surface protein A (OspA) from the Lyme disease spirochete Borrelia burgdorferi elicits protective immunity that blocks transmission of Borrelia from the tick vector to the vaccinated animal, and thus has been a focus of vaccine development. OspA has two globular domains that are connected via a unique single-layer beta-sheet. All anti-OspA monoclonal antibodies that block Borrelia transmission bind to conformational epitopes in the C-terminal domain of OspA, suggesting the possibility of using the C-terminal domain alone as a recombinant protein-based vaccine. The removal of ineffective parts from the OspA antigen may reduce side effects and lead to a safer vaccine. We prepared a C-terminal fragment of OspA by removing approximately 45% of residues from the N terminus. Although the fragment retained the native conformation and affinity to a protective antibody, its vaccine efficacy and conformational stability were significantly reduced with respect to full-length OspA. We successfully stabilized the fragment by replacing amino acid residues involved in buried salt-bridges with residues promoting hydrophobic interactions. The mutations promoted the vaccine efficacy of the redesigned fragment to a level comparable to that of the full-length protein, demonstrating the importance of the antigen stability for OspA's vaccine efficacy. Our strategy should be useful for further refining OspA-based vaccines and developing recombinant vaccines for other diseases.

ospC diversity in Borrelia burgdorferi: different hosts are different niches

The outer surface protein C (ospC) locus of the Lyme disease bacterium, Borrelia burgdorferi, is at least an order of magnitude more variable than other genes in the species. This variation is classified into 22 ospC major groups, 15 of which are found in the northeastern United States. The frequency distributions of ospC within populations suggest that this locus is under balancing selection. In multiple-niche polymorphism, a type of balancing selection, diversity within a population can be maintained when the environment is heterogeneous and no one genotype has the highest fitness in all environments. Genetically different individuals within vertebrate species and different vertebrate species constitute diverse environments for B. burgdorferi. We examined four important host species of B. burgdorferi and found that the strains that infected each species had different sets of ospC major groups. We found no variation among conspecific hosts in the ospC major groups of their infecting strains. These results suggest multiple niches create balancing selection at the ospC locus.

Sequence typing reveals extensive strain diversity of the Lyme borreliosis agents Borrelia burgdorferi in North America and Borrelia afzelii in Europe

The genetic polymorphism of Borrelia burgdorferi and Borrelia afzelii, two species that cause Lyme borreliosis, was estimated by sequence typing of four loci: the rrs-rrlA intergenic spacer (IGS) and the outer-membrane-protein gene p66 on the chromosome, and the outer-membrane-protein genes ospA and ospC on plasmids. The major sources of DNA for PCR amplification and sequencing were samples of the B. burgdorferi tick vector Ixodes scapularis, collected at a field site in an endemic region of the north-eastern United States, and the B. afzelii vector Ixodes ricinus, collected at a similar site in southern Sweden. The sequences were compared with those of reference strains and skin biopsy isolates, as well as database sequences. For B. burgdorferi, 10-13 alleles for each of the 4 loci, and a total of 9 distinct clonal lineages with linkage of all 4 loci, were found. For B. afzelii, 2 loci, ospC and IGS, were examined, and 11 IGS genotypes, 12 ospC alleles, and a total of 9 linkage groups were identified. The genetic variants of B. burgdorferi and B. afzelii among samples from the field sites accounted for the greater part of the genetic diversity previously reported from larger areas of the north-eastern United States and central and northern Europe. Although ospC alleles of both species had higher nucleotide diversity than other loci, the ospC locus showed evidence of intragenic recombination and was unsuitable for phylogenetic inference. In contrast, there was no detectable recombination at the IGS locus of B. burgdorferi. Moreover, beyond the signature nucleotides that specified 10 IGS genotypes, there were additional nucleotide polymorphisms that defined a total of 24 subtypes. Maximum-likelihood and parsimony cladograms of B. burgdorferi aligned IGS sequences revealed the subtype sequences to be terminal branches of clades, and the existence of at least three monophyletic lineages within B. burgdorferi. It is concluded that B. burgdorferi and B. afzelii have greater genetic diversity than had previously been estimated, and that the IGS locus alone is sufficient for strain typing and phylogenetic studies.

Lyme borreliosis agents and the genetics and sex of their vector, Ixodes ricinus

The tick Ixodes ricinus is responsible for the transmission and maintenance of a wide variety of pathogenic organisms in the Northern Hemisphere, among which Lyme disease represents a major threat to humans. Despite numerous studies, the epidemiology of the different bacterial species responsible for this disease remains unclear. Recent evidence for a sex-biased genetic structure of its European vector leads us to analyse the consequences of this pattern on Borrelia transmission. Here we show that male and female ticks are not equivalently infected by Borrelia burgdorferi, that Borrelia afzelii affects tick migration capabilities, especially for the most vagile sex (i.e., male) and that Lyme borreliosis agents are consequently vectorised in a much more complex way than usually thought. Such results change the epidemiological perception of Lyme borreliosis and suggest new co-evolutionary pathways between the ticks and the borrelia.

'Lyme disease': ancient engine of an unrecognized borreliosis pandemic?

Unexpectedly we have found large numbers of chronically ill Borrelia burgdorferi PCR- and seropositive patients in Houston, Texas, a zoonotically 'non-endemic' area. In order to understand this finding prior to sufficient data availability, we chose to examine critically currently accepted but troublesome 'Lyme disease' concepts. Our method was to analyze each foundation 'Lyme disease' premise within the context of available medical and veterinary literature, then to reconstruct the disease model consistent with the preponderance of that data. We find the present conceptualization of the illness seriously truncated, with a high likelihood of two distinct but connected forms of human B. burgdorferi infection. The yet-unrecognized form appears to have a broader clinical presentation, wider geographic distribution, and vastly greater prevalence. We conclude that 'Lyme disease' currently acknowledges only its zoonosis arm and is a limited conceptualization of a far more pervasive and unrecognized infection state that must be considered a global epidemic.

OspA immunization decreases transmission of Borrelia burgdorferi spirochetes from infected Peromyscus leucopus mice to larval Ixodes scapularis ticks

Recombinant outer surface protein A (OspA) vaccination of wild animal reservoirs has potential application for reducing Borrelia burgdorferi transmission in nature and subsequent risk of human infection. As a major reservoir host, the white-footed mouse (Peromyscus leucopus) is a candidate for a vaccination program designed to reduce infection prevalence in vector ticks. In this study we characterized the effect of various levels of immunization with recombinant OspA-glutathione transferase fusion protein on transmission dynamics from infected P. leucopus to larval ticks. Control mice were vaccinated with glutathione transferase alone. All mice were experimentally infected with B. burgdorferi before vaccination. The immune responses of the immunized mice were assessed by enzyme-linked immunosorbent assay for antibodies to OspA. Transmission of B. burgdorferi from infected mice was determined by xenodiagnosis with uninfected larval ticks. Spirochetes in ticks were counted by direct immunofluorescence assay. The concentration of antibody to OspA increased with each OspA vaccination but most markedly after the first and second vaccinations. In comparison with control mice, there was reduced transmission by OspA-vaccinated mice to uninfected ticks. One, two, or three doses of OspA reduced infection prevalence in xenodiagnostic ticks by 48%, 92%, or 99% and the numbers of spirochetes per tick by 84%, 98%, or 99%, respectively. This study suggests that vaccination of P. leucopus with OspA could reduce transmission to the tick vector in nature despite prior infection of the reservoir host.

Geographic uniformity of the Lyme disease spirochete (Borrelia burgdorferi) and its shared history with tick vector (Ixodes scapularis) in the Northeastern United States

Over 80% of reported cases of Lyme disease in the United States occur in coastal regions of northeastern and mid-Atlantic states. The genetic structure of the Lyme disease spirochete (Borrelia burgdorferi) and its main tick vector (Ixodes scapularis) was studied concurrently and comparatively by sampling natural populations of I. scapularis ticks along the East Coast from 1996 to 1998. Borrelia is genetically highly diverse at the outer surface protein ospC. Since Borrelia is highly clonal, the ospC alleles can be used to define clones. A newly designed reverse line blotting (RLB) assay shows that up to 10 Borrelia clones can infect a single tick. The clone frequencies in Borrelia populations are the same across the Northeast. On the other hand, I. scapularis populations show strong regional divergence (among northeastern, mid-Atlantic, and southern states) as well as local differentiation. The high genetic diversity within Borrelia populations and the disparity in the genetic structure between Borrelia and its tick vector are likely consequences of strong balancing selection on local Borrelia clones. Demographically, both Borrelia and I. scapularis populations in the Northeast show the characteristics of a species that has recently expanded from a population bottleneck. Major geological and ecological events, such as the last glacial maximum (18,000 years ago) and the modern-day expansion of tick habitats, are likely causes of the observed "founder effects" for the two organisms in the Northeast. We therefore conclude that the genetic structure of B. burgdorferi has been intimately shaped by the natural history of its main vector, the northern lineage of I. scapularis ticks.

Antigen polymorphism in Borrelia hermsii, a clonal pathogenic bacterium

The relapsing fever spirochete, Borrelia hermsii, escapes immune selection by alternating expression of surface lipoprotein alleles. The switch results from a duplicative transposition of one of several surface lipoprotein-encoding nucleotide sequences into the singular expression site. These nucleotide sequences constitute a large gene family whose diversity originated, in some cases, before the major divergences of Borrelia species. We have examined the B. hermsii vsp subfamily of alleles, which are carried on linear plasmids within each cell and maintained in several diverse copies as an antigenic archive. Each encodes a distinct serotype-specific protein. We sequenced more than 90% of the alleles within a single strain-B. hermsii strain HS1. A preponderance of allelic mosaicism suggests that intragenic recombination, coupled with selection imposed by host immune response, has driven diversification of the archived ensemble of vsp alleles. The recombinational diversification of vsp alleles generates change in the associated serotypes of the magnitude (30-40% amino acid differentiation) necessary for overcoming cross-reactivity of neutralizing antibodies. We conclude that evolution of vsp has occurred by punctuated occurrence of allelic differentiation, rather than by gradual selection of incremental point mutations that do not meet the threshold for antigenic diversity.

Defining and interpreting intraspecific molecular variation

Defining the extent and character of intraspecific genetic variation provides important information about gene function and organismal history. Powerful tests may be applied to sequenced alleles in order to critically examine whether natural selection is responsible for limiting or elevating intraspecific polymorphism in particular genes. Unconventional patterns of sequence variation and unusual allelic frequency distributions can be used to test whether genes encoding parasite antigens are being diversified by immune selection. The strikingly limited genetic variation in the falciparum malaria genome, and in human chromosomes encoding resistance to severe malaria, date the emergence of this disease to within the last few thousand years, illustrating the power of population genetic analysis to elucidate the history of host-parasite interactions. Coupling phylogenetic and geographic information and analyzing the rate of diversification in intraspecific gene trees provides new and rich sources of information on microbial evolution and epidemiology.

Likelihoods and simulation methods for a class of nonneutral population genetics models

Methods for simulating samples and sample statistics, under mutation-selection-drift equilibrium for a class of nonneutral population genetics models, and for evaluating the likelihood surface, in selection and mutation parameters, are developed and applied for observed data. The methods apply to large populations in settings in which selection is weak, in the sense that selection intensities, like mutation rates, are of the order of the inverse of the population size. General diploid selection is allowed, but the approach is currently restricted to models, such as the infinite alleles model and certain K-models, in which the type of a mutant allele does not depend on the type of its progenitor allele. The simulation methods have considerable advantages over available alternatives. No other methods currently seem practicable for approximating likelihood surfaces.

Distinct combinations of Borrelia burgdorferi sensu lato genospecies found in individual questing ticks from Europe

The genetic diversity of Borrelia burgdorferi sensu lato was assessed in individual adult Ixodes ricinus ticks from Europe by direct PCR amplification of spirochetal DNA followed by genospecies-specific hybridization. Analysis of mixed infections in the ticks showed that B. garinii and B. valaisiana segregate from B. afzelii. This and previous findings suggest that host complement interacts with spirochetes in the tick, thereby playing an important role in the ecology of Lyme borreliosis.

The structure of a local population of phytopathogenic Pseudomonas brassicacearum from agricultural soil indicates development under purifying selection pressure

Among the isolates of a bacterial community from a soil sample taken from an agricultural plot in northern Germany, a population consisting of 119 strains was obtained that was identified by 16S rDNA sequencing and genomic fingerprinting as belonging to the recently described species Pseudomonas brassicacearum. Analysis of the population structure by allozyme electrophoresis (11 loci) and random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR; four primers) showed higher resolution with the latter method. Both methods indicated the presence of three lineages, one of which dominated strongly. Stochastic tests derived from the neutral theory of evolution (including Slatkin's exact test, Watterson's homozygosity test and the Tajima test) indicated that the population had developed under strong purifying selection pressure. The presence of strains clearly divergent from the majority of the population can be explained by in situ evolution or by influx of strains as a result of migration or both. Phytopathogenicity of a P. brassicacearum strain determined with tomato plants reached the level obtained with the type strain of the known pathogen Pseudomonas corrugata. The results show that a selective sweep was identified in a local population. Previously, a local selective sweep had not been seen in several populations of different bacterial species from a variety of environmental habitats.

Genetic diversity of Borrelia burgdorferi sensu lato in ticks from mainland Portugal

To date Borrelia lusitaniae is the only genospecies of Borrelia burgdorferi sensu lato isolated from Ixodes ricinus ticks collected in Portugal and Tunisia. This suggests that the genospecies diversity of B. burgdorferi sensu lato decreases toward the southwestern margin of its Old World subtropical range. In order to further explore the genetic diversity of B. burgdorferi sensu lato from this region, 55 I. ricinus and 27 Hyalomma marginatum questing adults, collected during the spring of 1998 from a sylvatic habitat south of Lisbon, Portugal, were analyzed. Infection prevalences of 75% in I. ricinus ticks and 7% in H. marginatum ticks were detected by a nested PCR that targets the rrf (5S)-rrl (23S) spacer of B. burgdorferi sensu lato. Restriction fragment length polymorphism (RFLP) analysis of the I. ricinus-derived amplicons showed that the sequences in the majority of samples were similar to those of B. lusitaniae type strains (76% for strain PotiB1, 5% for strain PotiB3). Two novel RFLP patterns were obtained from 12% of the samples. The remaining 7% of samples gave mixed RFLP patterns. Phylogenetic analysis of rrf-rrl spacer sequences revealed a diverse population of B. lusitaniae in questing adult I. ricinus ticks (the sequences did not cluster with those of any other genospecies). This population consisted of 10 distinct sequence types, suggesting that multiple strains of B. lusitaniae were present in the local I. ricinus population. We hypothesize that B. lusitaniae has a narrow ecological niche that involves host species restricted to the Mediterranean Basin.

Methods for estimating gene frequencies and detecting selection in bacterial populations

Recent breakthroughs in molecular technology, most significantly the polymerase chain reaction (PCR) and in situ hybridization, have allowed the detection of genetic variation in bacterial communities without prior cultivation. These methods often produce data in the form of the presence or absence of alleles or genotypes, however, rather than counts of alleles. Using relative allele frequencies from presence-absence data as estimates of population allele frequencies tends to underestimate the frequencies of common alleles and overestimate those of rare ones, potentially biasing the results of a test of neutrality in favor of balancing selection. In this study, a maximum-likelihood estimator (MLE) of bacterial allele frequencies designed for use with presence-absence data is derived using an explicit stochastic model of the host infection (or bacterial sampling) process. The performance of the MLE is evaluated using computer simulation and a method is presented for evaluating the fit of estimated allele frequencies to the neutral infinite alleles model (IAM). The methods are applied to estimate allele frequencies at two outer surface protein loci (ospA and ospC) of the Lyme disease spirochete, Borrelia burgdorferi, infecting local populations of deer ticks (Ixodes scapularis) and to test the fit to a neutral IAM.

Population dynamics of a naturally occurring heterogeneous mixture of Borrelia burgdorferi clones

Two unique isolates of Borrelia burgdorferi, differing in plasmid content and outer surface protein C expression, were cultured on sequential captures of a single free-living Peromyscus leucopus mouse and were examined for differences in transmissibility. Both isolates were transmissible from inoculated C.B-17 mice to larval Ixodes scapularis ticks and, subsequently, from infected nymphal ticks to C3H/HeJ mice. Plasmid and protein analyses suggested that the original isolates were a mixed population of B. burgdorferi, and cloning by limiting dilution resulted in the identification of two clonal groups. In addition to being heterogeneous in plasmid and genomic macrorestriction analyses, the clones varied with respect to the electrophoretic mobilities and antigenicity of their OspC proteins, as shown by their reactivity to a panel of monoclonal antibodies. Plasmid analysis of sequential isolates from C3H mice experimentally infected with the primary isolate or various mixtures of its subclones showed an apparently random fluctuation in clonal dominance in the majority of mice. Surprisingly, mice infected with each subclone were permissive to superinfection with the heterologous subclone, despite the presence of anti-B. burgdorferi antibodies at the time of the secondary challenge. These results show conclusively that mice captured at Lyme disease enzootic sites may be infected by mixed populations of genetically and antigenically distinct B. burgdorferi clones and that these infections can be acquired by coinfection or by sequential infection. The lack of cross-immunization between clones existing within a naturally occurring population may play a role in the maintenance of the genetic heterogeneity of B. burgdorferi in nature.

Evolution of the linear DNA replicons of the Borrelia spirochetes

Members of the spirochete genus Borrelia carry numerous linear DNA replicons with covalently closed hairpin telomeres. The genome of one member of this genus, B. burgdorferi B31, has now been completely characterized and contains a linear chromosome, twelve linear plasmids and nine circular extra-chromosomal elements. The phylogenetic position of the Borrelia spirochetes strongly suggests that a progenitor with circular replicons acquired the ability to replicate linear DNA molecules.

Four clones of Borrelia burgdorferi sensu stricto cause invasive infection in humans

Lyme disease begins at the site of a tick bite, producing a primary infection with spread of the organism to secondary sites occurring early in the course of infection. A major outer surface protein expressed by the spirochete early in infection is outer surface protein C (OspC). In Borrelia burgdorferi sensu stricto, OspC is highly variable. Based on sequence divergence, alleles of ospC can be divided into 21 major groups. To assess whether strain differences defined by ospC group are linked to invasiveness and pathogenicity, we compared the frequency distributions of major ospC groups from ticks, from the primary erythema migrans skin lesion, and from secondary sites, principally from blood and spinal fluid. The frequency distribution of ospC groups from ticks is significantly different from that from primary sites, which in turn is significantly different from that from secondary sites. The major groups A, B, I, and K had higher frequencies in the primary sites than in ticks and were the only groups found in secondary sites. We define three categories of major ospC groups: one that is common in ticks but very rarely if ever causes human disease, a second that causes only local infection at the tick bite site, and a third that causes systemic disease. The finding that all systemic B. burgdorferi sensu stricto infections are associated with four ospC groups has importance in the diagnosis, treatment, and prevention of Lyme disease.