Analysis of linear plasmid dimers in Borrelia burgdorferi sensu lato isolates: implications concerning the potential mechanism of linear plasmid replication

Natural selection and recombination at host-interacting lipoprotein loci drive genome diversification of Lyme disease and related bacteria

Lyme disease, caused by spirochetes in the Borrelia burgdorferi sensu lato clade within the Borrelia genus, is transmitted by Ixodes ticks and is currently the most prevalent and rapidly expanding tick-borne disease in Europe and North America. We report complete genome sequences of 47 isolates that encompass all established species in this clade while highlighting the diversity of the widespread human pathogenic species B. burgdorferi. A similar set of plasmids has been maintained throughout Borrelia divergence, indicating that they are a key adaptive feature of this genus. Phylogenetic reconstruction of all sequenced Borrelia genomes revealed the original divergence of Eurasian and North American lineages and subsequent dispersals that introduced B. garinii, B. bavariensis, B. lusitaniae, B. valaisiana, and B. afzelii from East Asia to Europe and B. burgdorferi and B. finlandensis from North America to Europe. Molecular phylogenies of the universally present core replicons (chromosome and cp26 and lp54 plasmids) are highly consistent, revealing a strong clonal structure. Nonetheless, numerous inconsistencies between the genome and gene phylogenies indicate species dispersal, genetic exchanges, and rapid sequence evolution at plasmid-borne loci, including key host-interacting lipoprotein genes. While localized recombination occurs uniformly on the main chromosome at a rate comparable to mutation, lipoprotein-encoding loci are recombination hotspots on the plasmids, suggesting adaptive maintenance of recombinant alleles at loci directly interacting with the host. We conclude that within- and between-species recombination facilitates adaptive sequence evolution of host-interacting lipoprotein loci and contributes to human virulence despite a genome-wide clonal structure of its natural populations.

IMPORTANCE

Lyme disease (also called Lyme borreliosis in Europe), a condition caused by spirochete bacteria of the genus Borrelia, transmitted by hard-bodied Ixodes ticks, is currently the most prevalent and rapidly expanding tick-borne disease in the United States and Europe. Borrelia interspecies and intraspecies genome comparisons of Lyme disease-related bacteria are essential to reconstruct their evolutionary origins, track epidemiological spread, identify molecular mechanisms of human pathogenicity, and design molecular and ecological approaches to disease prevention, diagnosis, and treatment. These Lyme disease-associated bacteria harbor complex genomes that encode many genes that do not have homologs in other organisms and are distributed across multiple linear and circular plasmids. The functional significance of most of the plasmid-borne genes and the multipartite genome organization itself remains unknown. Here we sequenced, assembled, and analyzed whole genomes of 47 Borrelia isolates from around the world, including multiple isolates of the human pathogenic species. Our analysis elucidates the evolutionary origins, historical migration, and sources of genomic variability of these clinically important pathogens. We have developed web-based software tools (BorreliaBase.org) to facilitate dissemination and continued comparative analysis of Borrelia genomes to identify determinants of human pathogenicity.

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 preferential 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.

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.

A high fidelity approach to assembling the complex Borrelia genome

BACKGROUND

Bacteria of the Borrelia burgdorferi sensu lato (s.l.) complex can cause Lyme borreliosis. Different B. burgdorferi s.l. genospecies vary in their host and vector associations and human pathogenicity but the genetic basis for these adaptations is unresolved and requires completed and reliable genomes for comparative analyses. The de novo assembly of a complete Borrelia genome is challenging due to the high levels of complexity, represented by a high number of circular and linear plasmids that are dynamic, showing mosaic structure and sequence homology. Previous work demonstrated that even advanced approaches, such as a combination of short-read and long-read data, might lead to incomplete plasmid reconstruction. Here, using recently developed high-fidelity (HiFi) PacBio sequencing, we explored strategies to obtain gap-free, complete and high quality Borrelia genome assemblies. Optimizing genome assembly, quality control and refinement steps, we critically appraised existing techniques to create a workflow that lead to improved genome reconstruction.

RESULTS

Despite the latest available technologies, stand-alone sequencing and assembly methods are insufficient for the generation of complete and high quality Borrelia genome assemblies. We developed a workflow pipeline for the de novo genome assembly for Borrelia using several types of sequence data and incorporating multiple assemblers to recover the complete genome including both circular and linear plasmid sequences.

CONCLUSION

Our study demonstrates that, with HiFi data and an ensemble reconstruction pipeline with refinement steps, chromosomal and plasmid sequences can be fully resolved, even for complex genomes such as Borrelia. The presented pipeline may be of interest for the assembly of further complex microbial genomes.

Considerations on bacterial nucleoids

The classic genome organization of the bacterial chromosome is normally envisaged with all its genetic markers linked, thus forming a closed genetic circle of duplex stranded DNA (dsDNA) and several proteins in what it is called as "the bacterial nucleoid." This structure may be more or less corrugated depending on the physiological state of the bacterium (i.e., resting state or active growth) and is not surrounded by a double membrane as in eukayotic cells. The universality of the closed circle model in bacteria is however slowly changing, as new data emerge in different bacterial groups such as in Planctomycetes and related microorganisms, species of Borrelia, Streptomyces, Agrobacterium, or Phytoplasma. In these and possibly other microorganisms, the existence of complex formations of intracellular membranes or linear chromosomes is typical; all of these situations contributing to weakening the current cellular organization paradigm, i.e., prokaryotic vs eukaryotic cells.

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.

Comparison of MKP and BSK-H media for the cultivation and isolation of Borrelia burgdorferi sensu lato

The isolation of B. burgdorferi sensu lato requires the use of complex cultivation media. The aim of the study was to compare the usefulness of BSK-H (a commercial medium produced by HiMedia, India) and MKP medium. MKP and BSK-H media were prepared in accordance with the relevant protocols. Borrelia strains and skin culture biopsies were simultaneously inoculated into both media, incubated and checked for growth. Borrelial growth characteristics, isolation rates and characteristics of the isolated borreliae were analysed and compared. Initially, numbers of spirochaetes were higher in BSK-H than in MKP; however, in comparison with MKP, the strains subcultured in BSK-H medium were more frequently irregular, thin and non-motile, and rapidly died. In addition, the borrelial isolation rate from erythema migrans skin samples was higher in MKP than in BSK-H medium (108/171, 63.2% versus 70/171, 40.9%; p<0.0001). The far most frequently isolated species was Borrelia afzelii (92.9% and 97.2% strains isolated from BSK-H and MKP, respectively). Comparison of strains cultured from individual patients in both media showed differences in plasmid contents in 9/46 (19.6%) strain pairs, and protein profiles differed in 30/43 (69.8%) strain pairs, most often in the expression of OspC (in 27/28 patients OspC was expressed only in strains growing in MKP). BSK-H medium supports the growth of borrelial strains but MKP is superior with regard to the isolation rate, morphology and motility of strains. BSK-H medium supports fast initial growth of borreliae but this is followed by rapid deformation and death of the spirochaetes.

Borrelia Genotyping in Lyme Disease

Borrelia burgdorferisensu lato is the causative agent of Lyme borreliosis, multisystem disorder characterized by a wide spectrum of clinical manifestations. Different borrelia species can lead to distinct clinical presentations, but some species were associated with defined clinical manifestation likeBorrelia afzeliiwith skin manifestations,Borrelia gariniiwith central nervous system disorders andBorrelia burgdorferisensu stricto with Lyme arthritis.Ixodesticks represent the main vectors ofB. burgdorferisensu lato; wild animals, lizards and birds are the natural reservoir of borrelia. Genotyping of borrelia strains is of great importance for epidemiological, clinical, and evolutionary studies. Numerous methods are available for the genotyping ofB. burgdorferisensu lato based either on whole genome or PCR based typing. Typing methods differ in their approach and target, many of them were implemented more or less successfully for diagnostic purposes.

Highlights on molecular identification of closely related species

The term "complex" emerged in the literature at the beginning of the genomic era associated to taxonomy and grouping organisms that belong to different species but exhibited similar patterns according to their morphological, physiological and/or other phenotypic features. DNA-DNA hybridization values ~70% and high identity on 16S rRNA gene sequences were recommended for species delineation. Electrophoretic methods showed in some cases to be useful for species identification and population structure but the reproducibility was questionable. Later, the implementation of polyphasic approaches involving phenotypic and molecular methods brought new insights into the analysis of population structure and phylogeny of several "species complexes", allowing the identification of new closely related species. Likewise, the introduction of multilocus sequence typing and sequencing analysis of several genes offered an evolutionary perspective to the term "species complex". Several centres worldwide have recently released increasing genetic information on distinct microbial species. A brief review will be presented to highlight the definition of "species complex" for selected microorganisms, mainly the prokaryotic Acinetobacter calcoaceticus -Acinetobacter baumannii, Borrelia burgdorferi sensu lato, Burkholderia cepacia, Mycobacterium tuberculosis and Nocardia asteroides complexes, and the eukaryotic Aspergillus fumigatus, Leishmania donovani and Saccharomyces sensu stricto complexes. The members of these complexes may show distinct epidemiology, pathogenicity and susceptibility, turning critical their correct identification. Dynamics of prokaryotic and eukaryotic genomes can be very distinct and the term "species complex" should be carefully extended.

Requirements for Borrelia burgdorferi plasmid maintenance

Borrelia burgdorferi has multiple linear and circular plasmids that are faithfully replicated and partitioned as the bacterium grows and divides. The low copy number of these replicons implies that active partitioning contributes to plasmid stability. Analyzing the requirements for plasmid replication and partition in B. burgdorferi is complicated by the complexity of the genome and the possibility that products may act in trans. Consequently, we have studied the replication-partition region (bbb10-13) of the B. burgdorferi 26kb circular plasmid (cp26) in Escherichia coli, by fusion with a partition-defective miniF plasmid. Our analysis demonstrated that bbb10, bbb11, and bbb13 are required for stable miniF maintenance, whereas bbb12 is dispensable. To validate these results, we attempted to inactivate two of these genes in B. burgdorferi. bbb12 mutants were obtained at a typical frequency, suggesting that the bbb12 product is dispensable for cp26 maintenance as well. We could not directly measure cp26 stability in the bbb12 mutant, because cp26 carries essential genes, and bacteria that have lost cp26 are inviable. Conversely, we were unable to inactivate bbb10 on cp26 of B. burgdorferi. Our results suggest that bbb12 is dispensable for cp26 maintenance, whereas bbb10, bbb11, and bbb13 play crucial roles in that process.

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.

Detection of a large linear plasmid in Borrelia spielmanii isolate

Borrelia spielmanii belongs to human pathogenic species within the Borrelia burgdorferi sensu lato complex in Europe, which is a causative agent of Lyme disease. So far, the human disease caused by B. spielmanii has been associated with skin manifestations. The aim of the study was to analyze 4 human B. spielmanii isolates by pulsed-field gel electrophoresis and to localize genes of 3 important Borrelia proteins: OspA, DbpA, and VlsE. The analysis revealed variation within linear plasmid profiles among the strains; isolate PSigII contained a large plasmid of 100 kb compared with a 50 kb plasmid present in the 3 other B. spielmanii isolates, all carried the genes ospA and dbpA. Differences in the size of linear plasmids among the Borrelia strains may be a result of host-pathogen interactions, as the PSigII strain was the only strain of the 4 tested strains to be isolated from a patient with a previous history of Lyme disease, whereas 3 other patients were diagnosed with this disease for the first time.

Transposon insertion reveals pRM, a plasmid of Rickettsia monacensis

Until the recent discovery of pRF in Rickettsia felis, the obligate intracellular bacteria of the genus Rickettsia (Rickettsiales: Rickettsiaceae) were thought not to possess plasmids. We describe pRM, a plasmid from Rickettsia monacensis, which was detected by pulsed-field gel electrophoresis and Southern blot analyses of DNA from two independent R. monacensis populations transformed by transposon-mediated insertion of coupled green fluorescent protein and chloramphenicol acetyltransferase marker genes into pRM. Two-dimensional electrophoresis showed that pRM was present in rickettsial cells as circular and linear isomers. The 23,486-nucleotide (31.8% G/C) pRM plasmid was cloned from the transformant populations by chloramphenicol marker rescue of restriction enzyme-digested transformant DNA fragments and PCR using primers derived from sequences of overlapping restriction fragments. The plasmid was sequenced. Based on BLAST searches of the GenBank database, pRM contained 23 predicted genes or pseudogenes and was remarkably similar to the larger pRF plasmid. Two of the 23 genes were unique to pRM and pRF among sequenced rickettsial genomes, and 4 of the genes shared by pRM and pRF were otherwise found only on chromosomes of R. felis or the ancestral group rickettsiae R. bellii and R. canadensis. We obtained pulsed-field gel electrophoresis and Southern blot evidence for a plasmid in R. amblyommii isolate WB-8-2 that contained genes conserved between pRM and pRF. The pRM plasmid may provide a basis for the development of a rickettsial transformation vector.

Topology-dependent transcription in linear and circular plasmids of the segmented genome ofBorrelia burgdorferi

The segmented genome of Borrelia burgdorferi, a causative agent of Lyme disease, contains a mixture of over 20 linear and circular plasmids. Genes encoding five paralogous families of plasmid replication proteins are located on both circular and linear DNA molecules. The effect of DNA topology on the transcription of replication proteins from two B. burgdorferi plasmids, cp9 and lp17, was examined using quantitative reverse transcription polymerase chain reaction. Circular to linear conversion of a cp9-derived plasmid resulted in a 160-fold decrease in transcript levels of bbc01, believed to encode the replication initiator. A 14.9-fold reduction in plasmid copy number was also observed, resulting in a net 10.7-fold lower transcription level per gene copy on a linear versus a circular plasmid. In contrast, expression of the bbd14 replication initiator for the linear plasmid lp17 was 7.2 times higher per gene copy on a linear versus a circular plasmid. Topology-dependent transcription of these genes may help to block topological interconversions during genome evolution, offers a new avenue for global gene regulation and also has important implications for the design of genetic complementation experiments in B. burgdorferi.

Comparative genome hybridization reveals substantial variation among clinical isolates of Borrelia burgdorferi sensu stricto with different pathogenic properties

Clinical and murine studies suggest that there is a differential pathogenicity of different genotypes of Borrelia burgdorferi, the spirochetal agent of Lyme disease. Comparative genome hybridization was used to explore the relationship between different genotypes. The chromosomes of all studied isolates were highly conserved (>93%) with respect to both sequence and gene order. Plasmid sequences were substantially more diverse. Plasmids lp54, cp26, and cp32 were present in all tested isolates, and their sequences and gene order were conserved. The majority of linear plasmids showed variation both in terms of presence among different isolates and in terms of sequence and gene order. The data strongly imply that all B. burgdorferi clinical isolates contain linear plasmids related to each other, but the structure of these replicons may vary substantially from isolate to isolate. These alterations include deletions and presumed rearrangements that are likely to result in unique plasmid elements in many isolates. There is a strong correlation between complete genome hybridization profiles and other typing methods, which, in turn, also correlate to differences in pathogenicity. Because there is substantially less variation in the chromosomal and circular plasmid portions of the genome, the major differences in open reading frame content and genomic diversity among isolates are linear plasmid driven.

Hairpin telomeres and genome plasticity in Borrelia: all mixed up in the end

Spirochetes of the genus Borrelia have a highly unusual genome structure composed of over 20 replicons. Most of these replicons are linear and terminated by covalently closed hairpin ends or telomeres. Moreover, the linear replicons are affected by extensive DNA rearrangements, including telomere exchanges, DNA duplications, and harbour a large number of pseudogenes. The mechanism for the unusual genome plasticity in the linear replicons has remained elusive. The enzymatic machinery (the telomere resolvase ResT) responsible for generating the hairpin ends from replicative intermediates has recently been shown to also perform a reverse reaction that fuses telomeres on unrelated replicons. Infrequent stabilization of such fusion events over evolutionary time provides the first proposed biochemical mechanism for the DNA rearrangements that are so prominent in the linear replicons of B. burgdorferi.

In vitro culture of Borrelia garinii results in loss of flagella and decreased invasiveness

A virulent, low-passage culture of a tick-derived strain of Borrelia garinii was subjected to serial in vitro passages, from which inoculations were made into C3H/HeN mice. A full display of pathogenicity was observed through passage 4, as measured by cultures of ear punch biopsy samples and internal organs and determination of tibiotarsal joint swelling. Decreased dissemination through skin and infection of internal organs were observed beginning at passage 6. These losses correlated with both the selection of clones harboring 21% less flagella than the parent strain, as seen by electron microscopy, and loss of the motility of the higher passages, as demonstrated by a swarm assay. However, during the chronic phase (3 months after infection), spirochetes were cultured from the bladder and kidney of a mouse inoculated with passage 12. The kidney isolate had the same number of flagella and motility as the original low-passage isolate. Although we can't exclude the possibility that other subtle variations may be arising given the uncloned nature of the isolate, we have found a strong association between loss of flagella and decreased invasiveness. Arthritogenicity progressively decreased with passages, so that only 12.5% of chronically infected mice inoculated with passage 29 still presented with joint swelling, concurrent with a decrease in the staining intensity in a Southern blot with a vlsE-based probe. These results suggest a multifactorial model in which the number of flagella drives the invasiveness of this agent, while plasmid-associated factors are responsible for triggering arthritogenicity.

Analysis of mechanisms associated with loss of infectivity of clonal populations of Borrelia burgdorferi B31MI

Numerous studies have provided suggestive evidence that the loss of plasmids correlates with the loss of infectivity of the Lyme disease spirochetes. In this study we have further investigated this correlation. Clonal populations were obtained from the skin of a mouse infected for 3 months with a clonal population of Borrelia burgdorferi B31MI. The complete plasmid compositions of these populations were determined using a combination of PCR and Southern hybridization. The infectivities of clones differing in plasmid composition were tested using the C3H-HeJ murine model for Lyme disease. While several clones were found to be noninfectious, a correlation between the loss of a specific plasmid and loss of infectivity in the clones analyzed in this report was not observed. While it is clear from recent studies that the loss of some specific plasmids results in attenuated virulence, this study demonstrates that additional mechanisms also contribute to the loss of infectivity.

Telomere resolution in the Lyme disease spirochete

The genus Borrelia includes the causative agents of Lyme disease and relapsing fever. An unusual feature of these bacteria is a genome that includes linear DNA molecules with covalently closed hairpin ends referred to as telomeres. We have investigated the mechanism by which the hairpin telomeres are processed during replication. A synthetic 140 bp sequence having the predicted structure of a replicated telomere was shown to function as a viable substrate for telomere resolution in vivo, and was sufficient to convert a circular replicon to a linear form. Our results suggest that the final step in the replication of linear Borrelia replicons is a site-specific DNA breakage and reunion event to regenerate covalently closed hairpin ends. The telomere substrate described here will be valuable both for in vivo manipulation of linear DNA in Borrelia and for in vitro studies to identify and characterize the telomere resolvase.

Distribution of twelve linear extrachromosomal DNAs in natural isolates of Lyme disease spirochetes

We have analyzed a panel of independent North American isolates of the Lyme disease agent spirochete, Borrelia burgdorferi (sensu stricto), for the presence of linear plasmids with sequence similarities to the 12 linear plasmids present in the B. burgdorferi type strain, isolate B31. The frequency of similarities to probes from each of the 12 B31 plasmids varied from 13 to 100% in the strain panel examined, and these similarities usually reside on plasmids similar in size to the cognate B31 plasmid. Sequences similar to 5 of the 12 B31 plasmids were found in all of the isolates examined, and >66% of the panel members hybridized to probes from 4 other plasmids. Sequences similar to most of the B. burgdorferi B31 plasmid-derived DNA probes used were also found on linear plasmids in the related Eurasian Lyme agents Borrelia garinii and Borrelia afzelii; however, some of these plasmids had uniform but substantially different sizes from their B. burgdorferi counterparts.

Molecular and immunological analyses of the Borrelia turicatae Bdr protein family

Here, we describe the molecular and immunological characterization of the bdr gene family of Borrelia turicatae, a relapsing-fever spirochete. Nine bdr alleles belonging to two different subfamilies were sequenced and localized to linear plasmids. Anti-Bdr antiserum was generated and used to analyze Bdr expression in pre- and postinfection isogenic populations. The analyses presented here provide a detailed characterization of the Bdr proteins in a relapsing-fever spirochete species, enhancing our understanding of these proteins at the genus-wide level.

Evolutionary and molecular analyses of the Borrelia bdr super gene family: delineation of distinct sub-families and demonstration of the genus wide conservation of putative functional domains, structural properties and repeat motifs

B. turicatae, a causative agent of relapsing fever, carries a polymorphic gene family that is homologous to the bdr gene family of the Lyme disease spirochetes (previously referred to as the rep+ or ORF-E gene family). Here we demonstrate that bdr related genes are widely distributed among pathogenic Borrelia species and exist as large, polymorphic, plasmid carried, gene families. Twenty distinct bdr alleles were identified in isolates of the relapsing fever spirochete, B. hermsii, and were localized to linear plasmids. Cloning and sequence analyses demonstrate that the putative Bdr functional domains (i.e. the phosphorylation motifs and the transmembrane C-terminal domain) are conserved across the genus while other regions of these proteins exhibit variability. An assessment of the evolutionary relationships among all known Bdr protein sequences obtained from five pathogenic Borrelia species revealed that there are distinct Bdr sub-families. The recognition of distinct phyletic clusters serves as the basis of a revised and simplified nomenclature for the bdr proteins that can be applied genus wide. At the biological level the delineation of multiple bdr sub-families within isogeneic populations raises the possibility that there may be functional partitioning among alleles. In summary, the distribution and conservation of the Bdr proteins suggests that they are important in the biology/pathogenesis of the Borrelia at the genus wide level.

A bacterial genome in flux: the twelve linear and nine circular extrachromosomal DNAs in an infectious isolate of the Lyme disease spirochete Borrelia burgdorferi

We have determined that Borrelia burgdorferi strain B31 MI carries 21 extrachromosomal DNA elements, the largest number known for any bacterium. Among these are 12 linear and nine circular plasmids, whose sequences total 610 694 bp. We report here the nucleotide sequence of three linear and seven circular plasmids (comprising 290 546 bp) in this infectious isolate. This completes the genome sequencing project for this organism; its genome size is 1 521 419 bp (plus about 2000 bp of undetermined telomeric sequences). Analysis of the sequence implies that there has been extensive and sometimes rather recent DNA rearrangement among a number of the linear plasmids. Many of these events appear to have been mediated by recombinational processes that formed duplications. These many regions of similarity are reflected in the fact that most plasmid genes are members of one of the genome's 161 paralogous gene families; 107 of these gene families, which vary in size from two to 41 members, contain at least one plasmid gene. These rearrangements appear to have contributed to a surprisingly large number of apparently non-functional pseudogenes, a very unusual feature for a prokaryotic genome. The presence of these damaged genes suggests that some of the plasmids may be in a period of rapid evolution. The sequence predicts 535 plasmid genes >/=300 bp in length that may be intact and 167 apparently mutationally damaged and/or unexpressed genes (pseudogenes). The large majority, over 90%, of genes on these plasmids have no convincing similarity to genes outside Borrelia, suggesting that they perform specialized functions.

Genetic analysis of Borrelia garinii OspA serotype 4 strains associated with neuroborreliosis: evidence for extensive genetic homogeneity

Infection with Borrelia garinii outer surface protein (Osp) A serotype 4 strains has been correlated with the development of neuroborreliosis in Lyme borreliosis patients in Europe. OspA serotype 4 isolates have been recovered primarily from human cerebrospinal fluid, suggesting a tropism for this environment. Previous studies with monoclonal antibodies directed against OspA and OspC demonstrated that OspA serotype 4 strains are antigenically closely related. In view of the pronounced antigenic and genetic variability that has been noted in the Osps of other Borrelia isolates, we sought to determine if OspA serotype 4 strains represent a recently emerged clonal lineage of B. garinii. Toward this goal, a representative group of OspA serotype 4 strains was analyzed for traits that typically exhibit hypervariability among isolates that cause Lyme borreliosis. The following criteria were assessed: (i) ospC sequences, (ii) plasmid composition, (iii) genomic restriction fragment length polymorphism (RFLP) patterns, and (iv) the RFLP patterns of the upstream homology box (UHB) element which flanks members of the UHB gene family at their 5' end. Collectively, these analyses demonstrate genetic homogeneity, suggesting that OspA serotype 4 strains are a recently emerged clonal lineage with an apparent tropism for the central nervous system.

The bdr gene families of the Lyme disease and relapsing fever spirochetes: potential influence on biology, pathogenesis, and evolution

Species of the genus Borrelia cause human and animal infections, including Lyme disease, relapsing fever, and epizootic bovine abortion. The borrelial genome is unique among bacterial genomes in that it is composed of a linear chromosome and a series of linear and circular plasmids. The plasmids exhibit significant genetic redundancy and carry 175 paralogous gene families, most of unknown function. Homologous alleles on different plasmids could influence the organization and evolution of the Borrelia genome by serving as foci for interplasmid homologous recombination. The plasmid-carried Borrelia direct repeat (bdr) gene family encodes polymorphic, acidic proteins with putative phosphorylation sites and transmembrane domains. These proteins may play regulatory roles in Borrelia. We describe recent progress in the characterization of the Borrelia bdr genes and discuss the possible influence of this gene family on the biology, pathogenesis, and evolution of the Borrelia genome.

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.

The plasmid prophage N15: a linear DNA with covalently closed ends

Coliphage N15 is a temperate bacteriophage whose prophage is a linear plasmid molecule with covalently closed ends (telomeres). The N15 prophage provided the first example of such DNA in prokaryotes and, up to now, it is the only known example of a linear plasmid in Escherichia coli. The linear N15 mature phage DNA has single-stranded cohesive ends. The phage and plasmid prophage DNAs are circularly permuted. The nucleotide structure of the telomere-forming site tel RL in phage DNA corresponds to the structures of the terminal hairpin loops. It suggests a unique mechanism for conversion of the circular phage DNA to the linear plasmid form, which is performed by the prokaryotic telomerase (protelomerase). The results of a comparison of the protelomerase with integrases lead us to suggest that these proteins may have evolved from a common ancestor. The mechanism of plasmid N15 replication is unknown. We propose that the protelomerase participates in linear plasmid replication, acting as a resolvase of replicative intermediates that are tail-to-tail linear dimers. The sequence analysis of the N15 DNA showed that it represents an evolutionary 'link' between plasmids F, P1, P4 and lambdoid bacteriophages.

Molecular and evolutionary analysis of Borrelia burgdorferi 297 circular plasmid-encoded lipoproteins with OspE- and OspF-like leader peptides

We previously described two OspE and three OspF homologs in Borrelia burgdorferi 297 (D. R. Akins, S. F. Porcella, T. G. Popova, D. Shevchenko, S. I. Baker, M. Li, M. V. Norgard, and J. D. Radolf, Mol. Microbiol. 18:507-520, 1995; D. R. Akins, K. W. Bourell, M. J. Caimano, M. V. Norgard, and J. D. Radolf, J. Clin. Investig. 101:2240-2250, 1998). In this study, we characterized four additional lipoproteins with OspE/F-like leader peptides (Elps) and demonstrated that all are encoded on plasmids homologous to cp32 and cp18 from the B31 and N40 strains, respectively. Statistical analysis of sequence similarities using the binary comparison algorithm revealed that the nine lipoproteins from strain 297, as well as the OspE, OspF, and Erp proteins from the N40 and B31 strains, fall into three distinct families. Based upon the observation that these lipoproteins all contain highly conserved leader peptides, we now propose that the ancestors of each of the three families arose from gene fusion events which joined a common N terminus to unrelated proteins. Additionally, further sequence analysis of the strain 297 circular plasmids revealed that rearrangements appear to have played an important role in generating sequence diversity among the members of these three families and that recombinational events in the downstream flanking regions appear to have occurred independently of those within the lipoprotein-encoding genes. The association of hypervariable regions with genes which are differentially expressed and/or subject to immunological pressures suggests that the Lyme disease spirochete has exploited recombinatorial processes to foster its parasitic strategy and enhance its immunoevasiveness.

The diverse and dynamic structure of bacterial genomes

Bacterial genome sizes, which range from 500 to 10,000 kbp, are within the current scope of operation of large-scale nucleotide sequence determination facilities. To date, 8 complete bacterial genomes have been sequenced, and at least 40 more will be completed in the near future. Such projects give wonderfully detailed information concerning the structure of the organism's genes and the overall organization of the sequenced genomes. It will be very important to put this incredible wealth of detail into a larger biological picture: How does this information apply to the genomes of related genera, related species, or even other individuals from the same species? Recent advances in pulsed-field gel electrophoretic technology have facilitated the construction of complete and accurate physical maps of bacterial chromosomes, and the many maps constructed in the past decade have revealed unexpected and substantial differences in genome size and organization even among closely related bacteria. This review focuses on this recently appreciated plasticity in structure of bacterial genomes, and diversity in genome size, replicon geometry, and chromosome number are discussed at inter- and intraspecies levels.

Characterization of circular plasmid dimers in Borrelia burgdorferi

We have inactivated the ospC, oppAIV, and guaB genes on the 26-kb circular plasmid of Borrelia burgdorferi (cp26) by allelic exchange. On several occasions following such transformations, the cp26 of transformants had an aberrant mobility through agarose gels. Characterization of these cp26 molecules showed that the plasmid had dimerized. These dimers were quite stable during either selective or nonselective passage. Subsequent transformations with dimer DNA supported the hypothesis that in B. burgdorferi, transforming cp26 DNA most likely does not displace the resident homologous plasmid but rather must recombine in order to donate sequences that it carries. These serendipitous findings provide a mechanism for obtaining heterozygous complemented control strains when mutant phenotypes are characterized.

Cloning and molecular characterization of a multicopy, linear plasmid-carried, repeat motif-containing gene from Borrelia turicatae, a causative agent of relapsing fever

Borrelia turicatae is one of several spirochete species that can cause relapsing fever. Here, we describe the identification and characterization of a gene from B. turicatae and other relapsing-fever spirochetes that exhibits homology with the rep+ and ORF-E gene families of the Lyme disease spirochetes. This gene, which we have designated repA, encodes a putative protein of 30.2 kDa with an isoelectric point of 4.69. The central region of RepA harbors a series of amino acid repeat motifs which exhibit homology with casein kinase 2 phosphorylation sites. Through Southern hybridization analyses, we demonstrate that repA (or a closely related sequence) is multicopy in the relapsing-fever spirochetes and is carried on variably sized linear plasmids in both Borrelia parkeri and B. turicatae. Transcriptional analyses demonstrate that repA is expressed, albeit at low levels, during in vitro cultivation of B. turicatae. Transcriptional start site analysis revealed that repA is preceded by a consensus ribosomal binding site and an appropriately spaced promoter element. The sequence conservation, unique features, and multicopy status of repA and its homologs suggest that RepA may play an important genus-wide role in the biology of the Borrelia.

Analysis of the organization of multicopy linear- and circular-plasmid-carried open reading frames in Borrelia burgdorferi sensu lato isolates

Plasmid cp8.3 of Borrelia afzelii IP21 carries several open reading frames (ORFs) and a 184-bp inverted repeat (IR) element. It has been speculated that this plasmid may encode factors involved in virulence or infectivity. In this report, we have characterized the distribution, molecular variability, and organization of ORFs 1, 2, and 4 and the IR elements among isolates of the Borrelia burgdorferi sensu lato complex. ORFs 1 and 2 are contained within a segment of cp8.3 that is bordered by the IR elements, while ORF 4 resides just outside of the IR-bordered region. By PCR, ORF 4 was amplified from most isolates while ORFs 1 and 2 were amplified from only some B. afzelii isolates. However, Southern hybridization analyses with ORF 1, 2, and 4 probes detected related sequences even in some isolates that were PCR negative. The ORF restriction fragment length polymorphism patterns varied widely even among isolates of the same species. Two-dimensional contour-clamped homogeneous electric field-pulsed-field gel electrophoresis and Southern hybridization detected ORF 1-, 2-, and 4-related sequences on linear and circular plasmids. In addition, an ORF 4-related sequence was detected on a previously uncharacterized, circular plasmid that is greater than 70 kb in size. The IR elements originally identified on plasmid cp8.3 of B. afzelii IP21 were also analyzed by Southern hybridization. Related sequences were detected in some but not all B. burgdorferi sensu lato isolates. These sequences are carried on plasmids in addition to cp8.3 in some isolates. Single-primer PCR analyses demonstrated that in some isolates these sequences exist with IR orientation. The data presented here demonstrate that the IR elements and the ORF 1-, 2-, and 4-related sequences are multicopy and are variable in organization and in genomic location among isolates of the B. burgdorferi sensu lato complex. These analyses provide additional evidence for the highly variable organization of the plasmid component of the B. burgdorferi sensu lato genome.

An unusual illegitimate recombination occurs in the linear-plasmid-encoded outer-surface protein A gene of Borrelia afzelii

In this study, we describe an unusual illegitimate recombination in the linear-plasmid-encoded outer-surface protein A gene of Borrelia afzelii. A 96 bp DNA segment was deleted from the ospA structural gene of B. afzelii strain R9. The nature of the rearrangement suggested that it arose by a strand slippage mechanism, which was stimulated by a 18-mer palindromic sequence and 5-mer short direct repeats at both termini of the deleted DNA. The deleted sequence could form a complex hairpin structure suggesting that it may have played important roles in pausing of replication and slippaging of the nascent strand across the replication fork. In addition, the mutant strain was isolated from a chronic Lyme disease patient, implying that the variation mechanism may have been used by the borrelial strain to avoid host immune elimination.

Cloning and sequencing of a species-specific nucleotide fragment of Borrelia burgdorferi sensu stricto, which is repeated in several plasmids of the species

Among the etiological agents of Lyme disease, Borrelia burgdorferi sensu stricto strains carry a 16 kb plasmid, which did not hybridize to plasmids of B. garinii and B. afzelii strains. A 1271 bp DNA fragment of the 16 kb plasmid was cloned. It hybridized to several plasmids of this species (16, 27 and 55 kb). Sequencing of the cloned insert revealed a 327 bp ORF coding for a 14 kDa protein of unknown function, which could be expressed in E. coli. This ORF, conserved among B. burgdorferi sensu stricto strains, was carried by the same three plasmids.

Homology throughout the multiple 32-kilobase circular plasmids present in Lyme disease spirochetes

We have characterized seven different 32-kb circular plasmids carried by Borrelia burgdorferi isolate B31. Restriction endonuclease recognition site mapping and partial sequencing of these plasmids indicated that all seven are probably closely related to each other throughout their lengths and have substantial relationships to cp8.3, an 8.3-kb circular plasmid of B. burgdorferi sensu lato isolate Ip21. With the addition of the seven 32-kb plasmids, this bacterial strain is known to carry at least 10 linear and 9 circular plasmids. Variant cultures of B. burgdorferi B31 lacking one or more of the 32-kb circular plasmids are viable and, at least in some cases, infectious. We have examined a number of different natural isolates of Lyme disease borreliae and found that all of the B. burgdorferi sensu stricto isolates and most of the B. burgdorferi sensu lato isolates tested appear to carry multiple 32-kb circular plasmids related to those of B. burgdorferi B31. The ubiquity of these plasmids suggests that they may be important in the natural life cycle of these organisms. They may be highly conjugative plasmids or prophage genomes, which could prove to be useful in genetically manipulating B. burgdorferi.