Skip directly to search Skip directly to A to Z list Skip directly to page options Skip directly to site content

Volume 16, Number 7—July 2010


Geographic Differences in Genetic Locus Linkages for Borrelia burgdorferi

Bridgit Travinsky, Jonas Bunikis, and Alan G. BarbourComments to Author 
Author affiliations: University of California, Irvine, California, USA

Cite This Article


Highlight and copy the desired format.

EID Travinsky B, Bunikis J, Barbour AG. Geographic Differences in Genetic Locus Linkages for Borrelia burgdorferi. Emerg Infect Dis. 2010;16(7):1147-1150.
AMA Travinsky B, Bunikis J, Barbour AG. Geographic Differences in Genetic Locus Linkages for Borrelia burgdorferi. Emerging Infectious Diseases. 2010;16(7):1147-1150. doi:10.3201/eid1607.091452.
APA Travinsky, B., Bunikis, J., & Barbour, A. G. (2010). Geographic Differences in Genetic Locus Linkages for Borrelia burgdorferi. Emerging Infectious Diseases, 16(7), 1147-1150.


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

Most bacterial pathogens comprise a variety of strains in various proportions. For Borrelia burgdorferi, an agent of Lyme borreliosis, strains differ in their reservoir host preferences (1), propensities to disseminate in humans (2,3), and prevalences in ticks by geographic area (4,5). Strain identification of B. burgdorferi now is predominantly based on DNA sequences of either of 2 genetic loci: 1) the plasmid-borne, highly polymorphic ospC gene, which encodes outer surface protein C (6,7), or 2) the intergenic spacer (IGS) between the rrs and rrlA rDNA, here called IGS1. Other loci for genotyping are the plasmid-borne ospA gene (7) and the rrfA-rrlB rDNA intergenic spacer, here called IGS2 (8). The apparent clonality of B. burgdorferi was justification for inferring strain identity from a single locus (9,10), but the extent of genomewide genetic exchange in this species may have been underestimated (6).

Given reports of an association between disease severity and B. burgdorferi genotype (2,3), prediction of a strain’s virulence potential from its genotype has clinical, diagnostic, and epidemiologic relevance. But is a single locus sufficient for this assessment?

The Study

To investigate this issue, we determined sequences of ospC and IGS1 loci, and in selected cases the ospA and IGS2 loci, in 1,522 DNA extracts from B. burgdorferi–infected Ixodes scapularis nymphs collected from the northeastern, mid-Atlantic, and north-central United States during the summers of 2004, 2005, 2006, and 2007, as described (4,11). We also included results from 214 infected I. pacificus nymphs collected in Mendocino County, California (5); 20 infected I. pacificus adults from Contra Costa County, California (J. Bunikis and A.G. Barbour, unpub. data); and 10 B. burgdorferi genomes (strains B31, ZS7, 156a, 64b, 72a, 118a, WI91-23, 94a, 29805, and CA-11.2a), for which sequences are publicly available ( Multilocus sequence typing (MLST), based on 8 chromosomal housekeeping genes, had been carried out for several strains represented in the extracts (Table) (4,12). The corresponding MLST types of the 10 genome sequences were assigned by reference to a B. burgdorferi MLST database ( (12). For this study, we also determined the MLST type of strain CA8.

The methods for 1) DNA extraction from ticks (11), 2) PCR amplification of ospC, ospA, and IGS1 (7), 3) amplification of IGS2 (8), and 4) amplification of 8 chromosomal loci for MLST (12) have been described. Sequences for both strands were determined from either PCR products or cloned fragments with custom primers (7). We followed the basic nomenclature of Wang et al. (13) until, after exhausting the alphabet, we assigned both a letter and, arbitrarily, the number 3 (e.g., C3) when a new nucleotide sequence differed by >8% from known ospC alleles. We distinguished ospC variants with <1% sequence difference by adding a lowercase letter, e.g., Da and Db. Except for ospC D3 and Oa, novel polymorphisms were confirmed in at least 1 other sample. To simplify IGS1 nomenclature, we numbered types sequentially, beginning with the original 9 types (7); ospA alleles (7) and IGS2 loci were likewise sequentially numbered. The Appendix Table [PDF - 41 KB - 1 page] provides accession numbers for all sequences, as well as original and revised names for IGS1 sequences.


Thumbnail of A) Bayesian and maximum-likelihood phylogenetic inference of outer surface protein C (ospC) gene sequences and B) concatenated multilocus sequence typing (MLST) sequences of Borrelia burgdorferi. Sequences were aligned by codon. Labels at the tips refer to ospC alleles (A) or MLST (ST) and linked ospC alleles (B; Table). Consensus phylograms were the output of the MrBayes version 3.1.2 algorithm ( There were 500,000 generations with the first 1,000 disca

Figure. A) Bayesian and maximum-likelihood phylogenetic inference of outer surface protein C (ospC) gene sequences and B) concatenated multilocus sequence typing (MLST) sequences of Borrelia burgdorferi. Sequences were aligned by codon. Labels...

For 741 Ixodes ticks from northeastern and north-central United States or from northern California, 1 ospC allele was identified and sequenced. In the remaining samples, we found a mixture of strains or evidence of >2 ospC and/or >2 IGS sequences (9). In 678 (91%) of the 741 samples with a single ospC, the allele could be matched with particular IGS1 (Table). We identified 9 unique ospC sequences: Fc, Ob, Ub, A3, B3, C3, D3, E3, and F3, all from the north-central United States. Alleles H3 and I3 of California were recently reported by Girard et al. (5). Of 32 codon-aligned ospC sequences, 6 pairs and 1 trio (Fa, Fb, and Fc) differed in sequence by <1% (Figure, panel A). Nine novel IGS1 sequences, numbered 24–31 and 33, were discovered in samples from which ospC alleles were determined.

When we confined analysis to samples from northeastern states, we confirmed linkage disequilibrium between ospC and IGS1 loci (7,10,14). However, when results from north-central states and California were included, a different picture emerged (Table, Figure, panel B). Most of the ospC alleles showed concordance with the chromosomal loci; monophyletic MLST showed either the same ospC allele or a minor variant of it. However, in several instances, the ospC alleles were linked to different IGS1 sequences, different ospA sequences, and/or different MLST with internal nodes in common. We observed this linkage for ospC alleles A, G, Hb, and N. In the case of ospC Hb, the shared internal node was deep.

We applied the Simpson index of diversity, as implemented by Hunter and Gaston (15), to the data in the Table to compare the discriminatory power (DP) of genotyping on the basis of a combination of ospC and IGS1 sequences with genotyping by 8-locus MLST (12). For double-locus typing, there were 43 types were found for 678 strains; DP value was 0.96. For MLST in this data set, 36 types were shown for 554 strains; DP was 0.95. In the study of Hoen et al. in which selection was made for geographic isolation, 37 types were distributed among 78 strains; DP was 0.97 (4).


Dependence on a single locus for typing may falsely identify different lineages as the same, especially when the samples come from different regions. Other loci may be as informative as ospC or IGS1, but the abundance of extant sequences for these loci justifies their continued use. Uncertainties about the linkage of ospC and IGS1 usually can be resolved by sequencing the ospA allele (Table). IGS2 provided little additional information in this study.

One interpretation of these findings is that lateral gene transfer of all or nearly all of an ospC gene has occurred between different genetic lineages. We previously had not detected recombination at the IGS1 locus on the chromosome (7), but there may be recombination at other chromosomal loci, as well as plasmid loci (6). Besides extending the understanding of the geographic structuring of the B. burgdorferi population, the results indicate that the ospC allele does not fully represent the complexity of B. burgdorferi lineages; thus, inferring phenotypes on the basis of this single locus should be made with caution.

Ms Travinsky is a senior research associate in the Department of Microbiology and Molecular Genetics, University of California, Irvine. Her research interests include the genetic diversity and phylogeography of Borrelia species.


We thank Robert S. Lane for providing strain CA8.

This research was supported by Centers for Disease Control and Prevention Cooperative Agreement CI 00171-01 and National Institutes of Health grant AI065359.


  1. Brisson D, Dykhuizen DE. ospC diversity in Borrelia burgdorferi: different hosts are different niches. Genetics. 2004;168:71322. DOIPubMed
  2. Wormser GP, Brisson D, Liveris D, Hanincova K, Sandigursky S, Nowakowski J, Borrelia burgdorferi genotype predicts the capacity for hematogenous dissemination during early Lyme disease. J Infect Dis. 2008;198:135864. DOIPubMed
  3. Dykhuizen DE, Brisson D, Sandigursky S, Wormser GP, Nowakowski J, Nadelman RB, The propensity of different Borrelia burgdorferi sensu stricto genotypes to cause disseminated infections in humans. Am J Trop Med Hyg. 2008;78:80610.PubMed
  4. Hoen AG, Margos G, Bent SJ, Diuk-Wasser MA, Barbour AG, Kurtenbach K, Phylogeography of Borrelia burgdorferi in the eastern United States reflects multiple independent Lyme disease emergence events. Proc Natl Acad Sci U S A. 2009;106:150138. DOIPubMed
  5. Girard YA, Travinsky B, Schotthoefer A, Federova N, Eisen RJ, Eisen L, Population structure of the Lyme disease spirochete Borrelia burgdorferi in the western black-legged tick (Ixodes pacificus) in northern California. Appl Environ Microbiol. 2009;75:724352. DOIPubMed
  6. Qiu WG, Schutzer SE, Bruno JF, Attie O, Xu Y, Dunn JJ, Genetic exchange and plasmid transfers in Borrelia burgdorferi sensu stricto revealed by three-way genome comparisons and multilocus sequence typing. Proc Natl Acad Sci U S A. 2004;101:141505. DOIPubMed
  7. Bunikis J, Garpmo U, Tsao J, Berglund J, Fish D, Barbour AG. Sequence typing reveals extensive strain diversity of the Lyme borreliosis agents Borrelia burgdorferi in North America and Borrelia afzelii in Europe. Microbiology. 2004;150:174155. DOIPubMed
  8. Derdakova M, Beati L, Pet'ko B, Stanko M, Fish D. Genetic variability within Borrelia burgdorferi sensu lato genospecies established by PCR-single-strand conformation polymorphism analysis of the rrfA-rrlB intergenic spacer in Ixodes ricinus ticks from the Czech Republic. Appl Environ Microbiol. 2003;69:50916. DOIPubMed
  9. Qiu WG, Dykhuizen DE, Acosta MS, Luft BJ. Geographic uniformity of the Lyme disease spirochete (Borrelia burgdorferi) and its shared history with tick vector (Ixodes scapularis) in the northeastern United States. Genetics. 2002;160:83349.PubMed
  10. Hanincova K, Liveris D, Sandigursky S, Wormser GP, Schwartz I. Borrelia burgdorferi sensu stricto is clonal in patients with early Lyme borreliosis. Appl Environ Microbiol. 2008;74:500814. DOIPubMed
  11. Barbour AG, Bunikis J, Travinsky B, Hoen AG, Diuk-Wasser MA, Fish D, Niche partitioning of Borrelia burgdorferi and Borrelia miyamotoi in the same tick vector and mammalian reservoir species. Am J Trop Med Hyg. 2009;81:112031. DOIPubMed
  12. Margos G, Gatewood AG, Aanensen DM, Hanincova K, Terekhova D, Vollmer SA, MLST of housekeeping genes captures geographic population structure and suggests a European origin of Borrelia burgdorferi. Proc Natl Acad Sci U S A. 2008;105:87305. DOIPubMed
  13. Wang IN, Dykhuizen DE, Qiu W, Dunn JJ, Bosler EM, Luft BJ. Genetic diversity of ospC in a local population of Borrelia burgdorferi sensu stricto. Genetics. 1999;151:1530.PubMed
  14. Attie O, Bruno JF, Xu Y, Qiu D, Luft BJ, Qiu WG. Co-evolution of the outer surface protein C gene (ospC) and intraspecific lineages of Borrelia burgdorferi sensu stricto in the northeastern United States. Infect Genet Evol. 2007;7:112. DOIPubMed
  15. Hunter PR, Gaston MA. Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity. J Clin Microbiol. 1988;26:24656.PubMed



Cite This Article

DOI: 10.3201/eid1607.091452

Table of Contents – Volume 16, Number 7—July 2010


Please use the form below to submit correspondence to the authors or contact them at the following address:

Alan G. Barbour, University of California, Irvine, 3012 Hewitt, Irvine, CA 92697-4028, USA

character(s) remaining.

Comment submitted successfully, thank you for your feedback.