BackgroundLyme 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).ResultsRobust 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.ConclusionsIntra- 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.

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. G ENETIC discontinuity, the basis of biodiversity, is ubiquitous in prokaryotes as well as in eukaryotes. Most bacterial populations display a highly clonal genetic structure, in which the observable number of multilocus genotypes is far fewer than the number expected under the assumption of free recombination (Maynard Smith et al. 1993). Bacterial clonality was originally thought of as a result of a lack or rarity of recombination among asexually reproducing and independently evolving clones (Ochman and Selander 1984). Since then, molecular surveys of natural bacterial populations using protein electrophoresis, multilocus sequencing typing (MLST), and whole-genome PRJNA3, PRJNA28633, PRJNA19839, PRJNA29359, PRJNA28629, PRJNA29357, PRJNA21003, PRJNA19835, PRJNA28627, PRJNA21001, PRJNA29361, PRJNA28621, PRJNA19837, PRJNA20999, PRJNA28631, PRJNA29363, PRJNA17057, PRJNA19841, PRJNA12554, PRJNA28625, PRJNA29573, PRJNA19843, and PRJNA28635. 1 Present address: Odum School of Ecology, University of Georgia, Athens, GA 30602. sequencing revealed that horizontal genetic exchange is in fact often more frequent than point mutations in bacteria, including species known as strongly clonal (Maynard Smith et al. 1993;Feil and Spratt 2001;Didelot and Maiden 2010;Retc...

Analysis of natural host-parasite relationships reveals the evolutionary forces that shape the delicate and unique specificity characteristic of such interactions. The accessory long gland-reservoir complex of the wasp Leptopilina heterotoma (Figitidae) produces venom with virus-like particles. Upon delivery, venom components delay host larval development and completely block host immune responses. The host range of this Drosophila endoparasitoid notably includes the highly-studied model organism, Drosophila melanogaster. Categorization of 827 unigenes, using similarity as an indicator of putative homology, reveals that approximately 25% are novel or classified as hypothetical proteins. Most of the remaining unigenes are related to processes involved in signaling, cell cycle, and cell physiology including detoxification, protein biogenesis, and hormone production. Analysis of L. heterotoma’s predicted venom gland proteins demonstrates conservation among endo- and ectoparasitoids within the Apocrita (e.g., this wasp and the jewel wasp Nasonia vitripennis) and stinging aculeates (e.g., the honey bee and ants). Enzyme and KEGG pathway profiling predicts that kinases, esterases, and hydrolases may contribute to venom activity in this unique wasp. To our knowledge, this investigation marks the first functional genomic study for a natural parasitic wasp of Drosophila. Our findings will help explain how L. heterotoma shuts down its hosts’ immunity and shed light on the molecular basis of a natural arms race between these insects.