Multipartite Genome of Lyme Disease Borrelia: Structure, Variation and Prophages.

All members of the Borrelia genus that have been examined harbour a linear chromosome that is about 900 kbp in length, as well as a plethora of both linear and circular plasmids in the 5-220 kbp size range. Genome sequences for 27 Lyme disease Borrelia isolates have been determined since the elucidation of the B. burgdorferi B31 genome sequence in 1997. The chromosomes, which carry the vast majority of the housekeeping genes, appear to be very constant in gene content and organization across all Lyme disease Borrelia species. The content of the plasmids, which carry most of the genes that encode the differentially expressed surface proteins that interact with the spirochete's arthropod and vertebrate hosts, is much more variable. Lyme disease Borrelia isolates carry between 7-21 different plasmids, ranging in size from 5-84 kbp. All strains analyzed to date harbor three plasmids, cp26, lp54 and lp17. The plasmids are unusual, as compared to most bacterial plasmids, in that they contain many paralogous sequences, a large number of pseudogenes, and, in some cases, essential genes. In addition, a number of the plasmids have features indicating that they are prophages. Numerous methods have been developed for Lyme disease Borrelia strain typing. These have proven valuable for clinical and epidemiological studies, as well as phylogenomic and population genetic analyses. Increasingly, these approaches have been displaced by whole genome sequencing techniques. Some correlations between genome content and pathogenicity have been deduced, and comparative whole genome analyses promise future progress in this arena.

BpWrapper: BioPerl-based sequence and tree utilities for rapid prototyping of bioinformatics pipelines.

BACKGROUND: Automated bioinformatics workflows are more robust, easier to maintain, and results more reproducible when built with command-line utilities than with custom-coded scripts. Command-line utilities further benefit by relieving bioinformatics developers to learn the use of, or to interact directly with, biological software libraries. There is however a lack of command-line utilities that leverage popular Open Source biological software toolkits such as BioPerl ( http://bioperl.org ) to make many of the well-designed, robust, and routinely used biological classes available for a wider base of end users. RESULTS: Designed as standard utilities for UNIX-family operating systems, BpWrapper makes functionality of some of the most popular BioPerl modules readily accessible on the command line to novice as well as to experienced bioinformatics practitioners. The initial release of BpWrapper includes four utilities with concise command-line user interfaces, bioseq, bioaln, biotree, and biopop, specialized for manipulation of molecular sequences, sequence alignments, phylogenetic trees, and DNA polymorphisms, respectively. Over a hundred methods are currently available as command-line options and new methods are easily incorporated. Performance of BpWrapper utilities lags that of precompiled utilities while equivalent to that of other utilities based on BioPerl. BpWrapper has been tested on BioPerl Release 1.6, Perl versions 5.10.1 to 5.25.10, and operating systems including Apple macOS, Microsoft Windows, and GNU/Linux. Release code is available from the Comprehensive Perl Archive Network (CPAN) at https://metacpan.org/pod/Bio::BPWrapper . Source code is available on GitHub at https://github.com/bioperl/p5-bpwrapper . CONCLUSIONS: BpWrapper improves on existing sequence utilities by following the design principles of Unix text utilities such including a concise user interface, extensive command-line options, and standard input/output for serialized operations. Further, dozens of novel methods for manipulation of sequences, alignments, and phylogenetic trees, unavailable in existing utilities (e.g., EMBOSS, Newick Utilities, and FAST), are provided. Bioinformaticians should find BpWrapper useful for rapid prototyping of workflows on the command-line without creating custom scripts for comparative genomics and other bioinformatics applications.

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.

Genotyping and Quantifying Lyme Pathogen Strains by Deep Sequencing of the Outer Surface Protein C (ospC) Locus.

A mixed infection of a single tick or host by Lyme disease spirochetes is common and a unique challenge for the diagnosis, treatment, and surveillance of Lyme disease. Here, we describe a novel protocol for differentiating Lyme strains on the basis of deep sequencing of the hypervariable outer surface protein C locus (ospC). Improving upon the traditional DNA-DNA hybridization method, the next-generation sequencing-based protocol is high throughput, quantitative, and able to detect new pathogen strains. We applied the method to more than one hundred infected Ixodes scapularis ticks collected from New York State, USA, in 2015 and 2016. An analysis of strain distributions within individual ticks suggests an overabundance of multiple infections by five or more strains, inhibitory interactions among coinfecting strains, and the presence of a new strain closely related to Borreliella bissettiae A supporting bioinformatics pipeline has been developed. The newly designed pair of universal ospC primers target intergenic sequences conserved among all known Lyme pathogens. The protocol could be used for culture-free identification and quantification of Lyme pathogens in wildlife and potentially in clinical specimens.

Proteome-wide analysis of mutant p53 targets in breast cancer identifies new levels of gain-of-function that influence PARP, PCNA, and MCM4.

The gain-of-function mutant p53 (mtp53) transcriptome has been studied, but, to date, no detailed analysis of the mtp53-associated proteome has been described. We coupled cell fractionation with stable isotope labeling with amino acids in cell culture (SILAC) and inducible knockdown of endogenous mtp53 to determine the mtp53-driven proteome. Our fractionation data highlight the underappreciated biology that missense mtp53 proteins R273H, R280K, and L194F are tightly associated with chromatin. Using SILAC coupled to tandem MS, we identified that R273H mtp53 expression in MDA-MB-468 breast cancer cells up- and down-regulated multiple proteins and metabolic pathways. Here we provide the data set obtained from sequencing 73,154 peptide pairs that then corresponded to 3,010 proteins detected under reciprocal labeling conditions. Importantly, the high impact regulated targets included the previously identified transcriptionally regulated mevalonate pathway proteins but also identified two new levels of mtp53 protein regulation for nontranscriptional targets. Interestingly, mtp53 depletion profoundly influenced poly(ADP ribose) polymerase 1 (PARP1) localization, with increased cytoplasmic and decreased chromatin-associated protein. An enzymatic PARP shift occurred with high mtp53 expression, resulting in increased poly-ADP-ribosylated proteins in the nucleus. Mtp53 increased the level of proliferating cell nuclear antigen (PCNA) and minichromosome maintenance 4 (MCM4) proteins without changing the amount of pcna and mcm4 transcripts. Pathway enrichment analysis ranked the DNA replication pathway above the cholesterol biosynthesis pathway as a R273H mtp53 activated proteomic target. Knowledge of the proteome diversity driven by mtp53 suggests that DNA replication and repair pathways are major targets of mtp53 and highlights consideration of combination chemotherapeutic strategies targeting cholesterol biosynthesis and PARP inhibition.

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.

BorreliaBase: a phylogeny-centered browser of Borrelia genomes.

BACKGROUND: The bacterial genus Borrelia (phylum Spirochaetes) consists of two groups of pathogens represented respectively by B. burgdorferi, the agent of Lyme borreliosis, and B. hermsii, the agent of tick-borne relapsing fever. The number of publicly available Borrelia genomic sequences is growing rapidly with the discovery and sequencing of Borrelia strains worldwide. There is however a lack of dedicated online databases to facilitate comparative analyses of Borrelia genomes. DESCRIPTION: We have developed BorreliaBase, an online database for comparative browsing of Borrelia genomes. The database is currently populated with sequences from 35 genomes of eight Lyme-borreliosis (LB) group Borrelia species and 7 Relapsing-fever (RF) group Borrelia species. Distinct from genome repositories and aggregator databases, BorreliaBase serves manually curated comparative-genomic data including genome-based phylogeny, genome synteny, and sequence alignments of orthologous genes and intergenic spacers. CONCLUSIONS: With a genome phylogeny at its center, BorreliaBase allows online identification of hypervariable lipoprotein genes, potential regulatory elements, and recombination footprints by providing evolution-based expectations of sequence variability at each genomic locus. The phylo-centric design of BorreliaBase (http://borreliabase.org) is a novel model for interactive browsing and comparative analysis of bacterial genomes online.

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.

Whole-genome sequences of Borrelia bissettii, Borrelia valaisiana, and Borrelia spielmanii.

It has been known for decades that human Lyme disease is caused by the three spirochete species Borrelia burgdorferi, Borrelia afzelii, and Borrelia garinii. Recently, Borrelia valaisiana, Borrelia spielmanii, and Borrelia bissettii have been associated with Lyme disease. We report the complete genome sequences of B. valaisiana VS116, B. spielmanii A14S, and B. bissettii DN127.

Accelerated and adaptive evolution of yeast sexual adhesins.

There is a recent emergence of interest in the genes involved in gametic recognition as drivers of reproductive isolation. The recent population genomic sequencing of two species of sexually primitive yeasts (Liti G, Carter DM, Moses AM, Warringer J, Parts L, James SA, Davey RP, Roberts IN, Burt A, Koufopanou V et al. [23 co-authors]. 2009. Population genomics of domestic and wild yeasts. Nature 458:337-341.) has provided data for systematic study of the roles these genes play in the early evolution of sex and speciation. Here, we discovered that among genes encoding cell surface proteins, the sexual adhesin genes have evolved significantly more rapidly than others, both within and between Saccharomyces cerevisiae and its closest relative S. paradoxus. This result was supported by analyses using the PAML pairwise model, a modified McDonald-Kreitman test, and the PAML branch model. Moreover, using a combination of a new statistic of neutrality, an information theory-based measure of evolutionary variability, and functional characterization of amino acid changes, we found that a higher proportion of amino acid changes are fixed in the sexual adhesins than in other proteins and a greater proportion of the fixed amino acid changes either between the two species or the two subgroups of S. paradoxus are functionally dissimilar or radically different. These results suggest that the accelerated evolution of sexual adhesin genes may facilitate speciation, or incipient speciation, and promote sexual selection in general.

A Bioinformatics Module for Use in an Introductory Biology Laboratory.

Since biomedical science has become increasingly data-intensive, acquisition of computational and quantitative skills by science students has become more important. For non-science students, an introduction to biomedical databases and their applications promotes the development of a scientifically literate population. Because typical college introductory biology laboratories do not include experiences of this type, we present a bioinformatics module that can easily be included in a 90-minute session of a biology course for both majors and non-majors. Students completing this computational, inquiry-based module observed the value of computer-assisted analysis. The module gave students an understanding of how to read files in a biological database (GenBank) and how to use a software tool (BLAST) to mine the database.

Maximum antigen diversification in a lyme bacterial population and evolutionary strategies to overcome pathogen diversity.

Natural populations of pathogens and their hosts are engaged in an arms race in which the pathogens diversify to escape host immunity while the hosts evolve novel immunity. This co-evolutionary process poses a fundamental challenge to the development of broadly effective vaccines and diagnostics against a diversifying pathogen. Based on surveys of natural allele frequencies and experimental immunization of mice, we show high antigenic specificities of natural variants of the outer surface protein C (OspC), a dominant antigen of a Lyme Disease-causing bacterium (Borrelia burgdorferi). To overcome the challenge of OspC antigenic diversity to clinical development of preventive measures, we implemented a number of evolution-informed strategies to broaden OspC antigenic reactivity. In particular, the centroid algorithm-a genetic algorithm to generate sequences that minimize amino-acid differences with natural variants-generated synthetic OspC analogs with the greatest promise as diagnostic and vaccine candidates against diverse Lyme pathogen strains co-existing in the Northeast United States. Mechanistically, we propose a model of maximum antigen diversification (MAD) mediated by amino-acid variations distributed across the hypervariable regions on the OspC molecule. Under the MAD hypothesis, evolutionary centroids display broad cross-reactivity by occupying the central void in the antigenic space excavated by diversifying natural variants. In contrast to vaccine designs based on concatenated epitopes, the evolutionary algorithms generate analogs of natural antigens and are automated. The novel centroid algorithm and the evolutionary antigen designs based on consensus and ancestral sequences have broad implications for combating diversifying pathogens driven by pathogen-host co-evolution.

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.

Whole genome sequence of an unusual Borrelia burgdorferi sensu lato isolate.

Human Lyme disease is caused by a number of related Borrelia burgdorferi sensu lato species. We report here the complete genome sequence of Borrelia sp. isolate SV1 from Finland. This isolate is to date the closest known relative of B. burgdorferi sensu stricto, but it is sufficiently genetically distinct from that species that it and its close relatives warrant its candidacy for new-species status. We suggest that this isolate should be named "Borrelia finlandensis."

Whole-genome sequences of two Borrelia afzelii and two Borrelia garinii Lyme disease agent isolates.

Human Lyme disease is commonly caused by several species of spirochetes in the Borrelia genus. In Eurasia these species are largely Borrelia afzelii, B. garinii, B. burgdorferi, and B. bavariensis sp. nov. Whole-genome sequencing is an excellent tool for investigating and understanding the influence of bacterial diversity on the pathogenesis and etiology of Lyme disease. We report here the whole-genome sequences of four isolates from two of the Borrelia species that cause human Lyme disease, B. afzelii isolates ACA-1 and PKo and B. garinii isolates PBr and Far04.

Mutational and phylogenetic analyses of the mycobacterial mbt gene cluster.

The mycobactin siderophore system is present in many Mycobacterium species, including M. tuberculosis and other clinically relevant mycobacteria. This siderophore system is believed to be utilized by both pathogenic and nonpathogenic mycobacteria for iron acquisition in both in vivo and ex vivo iron-limiting environments, respectively. Several M. tuberculosis genes located in a so-called mbt gene cluster have been predicted to be required for the biosynthesis of the core scaffold of mycobactin based on sequence analysis. A systematic and controlled mutational analysis probing the hypothesized essential nature of each of these genes for mycobactin production has been lacking. The degree of conservation of mbt gene cluster orthologs remains to be investigated as well. In this study, we sought to conclusively establish whether each of nine mbt genes was required for mycobactin production and to examine the conservation of gene clusters orthologous to the M. tuberculosis mbt gene cluster in other bacteria. We report a systematic mutational analysis of the mbt gene cluster ortholog found in Mycobacterium smegmatis. This mutational analysis demonstrates that eight of the nine mbt genes investigated are essential for mycobactin production. Our genome mining and phylogenetic analyses reveal the presence of orthologous mbt gene clusters in several bacterial species. These gene clusters display significant organizational differences originating from an intricate evolutionary path that might have included horizontal gene transfers. Altogether, the findings reported herein advance our understanding of the genetic requirements for the biosynthesis of an important mycobacterial secondary metabolite with relevance to virulence.

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.

Conserved processes and lineage-specific proteins in fungal cell wall evolution.

The cell wall is a defining organelle that differentiates fungi from its sister clades in the opisthokont superkingdom. With a sensitive technique to align low-complexity protein sequences, we have identified 187 cell wall-related proteins in Saccharomyces cerevisiae and determined the presence or absence of homologs in 17 other fungal genomes. There were both conserved and lineage-specific cell wall proteins, and the degree of conservation was strongly correlated with protein function. Some functional classes were poorly conserved and lineage specific: adhesins, structural wall glycoprotein components, and unannotated open reading frames. These proteins are primarily those that are constituents of the walls themselves. On the other hand, glycosyl hydrolases and transferases, proteases, lipases, proteins in the glycosyl phosphatidyl-inositol-protein synthesis pathway, and chaperones were strongly conserved. Many of these proteins are also conserved in other eukaryotes and are associated with wall synthesis in plants. This gene conservation, along with known similarities in wall architecture, implies that the basic architecture of fungal walls is ancestral to the divergence of the ascomycetes and basidiomycetes. The contrasting lineage specificity of wall resident proteins implies diversification. Therefore, fungal cell walls consist of rapidly diversifying proteins that are assembled by the products of an ancestral and conserved set of genes.

Disruption of the p53-Mdm2 complex by Nutlin-3 reveals different cancer cell phenotypes.

INTRODUCTION: Mdm2 inhibits p53 transactivation by forming a p53-Mdm2 complex on chromatin. Upon DNA damage-induced complex disruption, such latent p53 can be activated, but in cells overexpressing Mdm2 because of a homozygous single nucleotide polymorphism at position 309 (T --> G) of mdm2, the complex is highly stable and cannot be disrupted by DNA damage, rendering p53 inactive. METHODS: To determine whether the p53 response phenotype is influenced differentially in cells with variable mdm2 genotypes, we compared responses to DNA damage and targeted p53-Mdm2 complex disruption by Nutlin-3 in the following wild-type p53 human cancer cell lines: A875 and CCF-STTG-1 (G/G for mdm2 SNP309), SJSA-1 (mdm2 genomic amplification and T/T for mdm2 SNP309), MCF-7 (estrogen-induced Mdm2 overexpression and T/G for mdm2 SNP309), ML-1 and H460 (T/T for mdm2 SNP309), and K562 (p53-null and T/G for mdm2 SNP309). We also examined mdm2 gene-splicing patterns in these lines by cloning and sequencing analyses. RESULTS: While Mdm2-overexpressing G/G cells were resistant to p53 activation by DNA damage, they were sensitive to Nutlin-3. Strikingly, the p53 G1 checkpoint in G/G cells was activated by Nutlin-3 but not by etoposide, whereas in other Mdm2-overexpressing cells, both drugs activated p53 and subsequent G1 arrest or apoptosis. cDNA clones lacking exons 5-9 were generated at a high frequency in cells overexpressing Mdm2. CONCLUSION: Nutlin-3 and DNA damage distinguish a differential phenotype in human cancer cells with G/G mdm2 SNP309 from other Mdm2 overexpressers. Categorization of the Mdm2 isoforms produced and their influence on p53 activity will help in characterization and treatment development for different cancers.

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.