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

BACKGROUND: Bacteria from the genus Borrelia are known to harbor numerous linear and circular plasmids. We report here a comparative analysis of the nucleotide sequences of 236 plasmids present in fourteen independent isolates of the Lyme disease agent B. burgdorferi. RESULTS: We have sequenced the genomes of 14 B. burgdorferi sensu stricto isolates that carry a total of 236 plasmids. These individual isolates carry between seven and 23 plasmids. Their chromosomes, the cp26 and cp32 circular plasmids, as well as the lp54 linear plasmid, are quite evolutionarily stable; however, the remaining plasmids have undergone numerous non-homologous and often duplicative recombination events. We identify 32 different putative plasmid compatibility types among the 236 plasmids, of which 15 are (usually) circular and 17 are linear. Because of past rearrangements, any given gene, even though it might be universally present in these isolates, is often found on different linear plasmid compatibility types in different isolates. For example, the arp gene and the vls cassette region are present on plasmids of four and five different compatibility types, respectively, in different isolates. A majority of the plasmid types have more than one organizationally different subtype, and the number of such variants ranges from one to eight among the 18 linear plasmid types. In spite of this substantial organizational diversity, the plasmids are not so variable that every isolate has a novel version of every plasmid (i.e., there appears to be a limited number of extant plasmid subtypes). CONCLUSIONS: Although there have been many past recombination events, both homologous and nonhomologous, among the plasmids, particular organizational variants of these plasmids correlate with particular chromosomal genotypes, suggesting that there has not been rapid horizontal transfer of whole linear plasmids among B. burgdorferi lineages. We argue that plasmid rearrangements are essentially non-revertable and are present at a frequency of only about 0.65% that of single nucleotide changes, making rearrangement-derived novel junctions (mosaic boundaries) ideal phylogenetic markers in the study of B. burgdorferi population structure and plasmid evolution and exchange. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12864-017-3553-5) contains supplementary material, which is available to authorized users.