In this study we have investigated the roles of putative R-M genes in the type strain of the Lyme disease spirochete,
B. burgdorferi B31. Although steady progress has been made since
B. burgdorferi was first transformed in 1994, genetic manipulation of
B. burgdorferi remains inefficient, particularly when working with infectious clones that retain most of the plasmids (
50). Bioinformatics analyses, coupled with transformation of strain B31 clones with various plasmid backgrounds, previously pinpointed the
bbe02 and
bbq67 plasmid loci, on lp25 and lp56, respectively, as likely components of the R-M system of
B. burgdorferi (
11,
19,
29,
34). Kawabata and colleagues highlighted the roles of these gene products in limiting stable introduction of exogenous DNA by inactivation of
bbe02 in a
B. burgdorferi clone that lacked linear plasmid lp56, which dramatically increased the transformation efficiency of shuttle vector DNA from
E. coli (
27). They and we have presumed that
bbq67 is the only gene on lp56 that influences transformation phenotype, but a formal demonstration of this is needed. A subsequent study demonstrated that CpG methylation of foreign plasmid DNA enhanced transformation of
B. burgdorferi carrying lp56, consistent with an R-M system encoded by
bbq67 (
13). High-throughput genetic transformation techniques, like transposon mutagenesis, have been effective only in B31 strains that lack lp25 and lp56 (
6,
54).
We extended the observation that
bbe02 and
bbq67 limit stable introduction of foreign DNA into
B. burgdorferi, using three closely related shuttle vectors that differ primarily in the selectable markers they carry (Fig. A). Transformation with all three plasmids prepared from
E. coli was decreased by the presence of
bbe02 and
bbq67 (or the lp25 and lp56 plasmids that carry them) in targeted recipient bacteria, although the pBSV2 shuttle vector appeared to be more sensitive than pKFSS1 and pBSV2G to the barrier presented by
bbe02 (Fig. ). The
bbh09 gene on lp28-3 is a homolog of
bbe02 (
19,
34), yet we have not detected any influence of lp28-3 on shuttle vector transformations in unrelated studies (data not shown). This may indicate that these shuttle vectors do not possess any sites recognized/cleaved by BBH09 or that the enzyme is not active or made in these clones. Transformation studies in other bacteria have shown that methylated plasmid DNA from
E. coli can be restricted unless the recipient bacterium contains a similar methyltransferase activity (
12,
37,
38,
39,
41). Hughes and Johnson (
23) detected Dam methylation of DNA by the relapsing fever spirochete
Borrelia hermsii but did not detect Dam methylation by most
B. burgdorferi strains, including B31. We did not observe a major influence on transformation of
B. burgdorferi by modification of shuttle vector DNA by Dam, Dcm, or Hsd methyltransferases in
E. coli (Fig. ). However, transformation of
B. burgdorferi carrying both
bbe02 and
bbq67 with the shuttle vector pBSV2 was successful only when plasmid DNA was obtained from
E. coli lacking
dam and subsequently modified
in vitro with the CpG methyltransferase M.SssI.
The data presented in this study demonstrate for the first time that the bbe02 and bbq67 gene products both modify endogenous DNA and restrict foreign DNA lacking similar modification in a sequence-specific fashion, as predicted for R-M systems. This conclusion was supported by experiments in which increased transformation efficiencies were obtained when B. burgdorferi was both the source and recipient of transforming shuttle vector DNA (Fig. ). The transformation barrier presented by bbe02 and bbq67 in recipient bacteria was overcome if shuttle vector DNA came from B. burgdorferi donors that also carried these loci. The inferred modification of shuttle vector DNA by bbe02 and bbq67 gene products did not limit transformation into recipient B. burgdorferi, consistent with recognition and cleavage of distinct and unmethylated DNA sequences by these putative R-M systems. The highest transformation efficiencies were obtained with B. burgdorferi recipients that lack bbe02 or bbq67 (Fig. and Table ), suggesting incomplete modification and only partial protection of restriction sites on endogenous shuttle vector DNA by bacteria carrying bbe02 and bbq67. Significantly, transformation of B. burgdorferi containing bbe02 and bbq67 was dramatically enhanced with plasmid DNA obtained from B. burgdorferi relative to shuttle vector DNA from E. coli, even with in vitro CpG modification (Fig. ). Together these results suggest that the bbe02 and bbq67 loci of B. burgdorferi encode R-M enzymes that recognize and cleave distinct unmodified sequences.
NCBI and PROSITE analyses identified bbe02 and bbq67 as homologs of adenine methyltransferase. We have provided the first direct evidence supporting this designation using an antiserum that recognizes DNA containing N6-methyladenine. Antibody binding was observed with DNA from B. burgdorferi containing bbe02 or bbq67 when used in a Southwestern blot format (Fig. ). Dot blot assays with the same antiserum, in which genomic DNA from these strains was focused in a single spot, indicated the presence of weak adenine methylation in the strain lacking both bbe02 and bbq67 (lp56) (data not shown). This result suggests that B. burgdorferi may contain additional adenine methyltransferases but with more limited activities than BBE02 and BBQ67.
Further support for adenine methylation by the
bbq67 gene product was suggested by protection of certain RsaI sites on plasmid DNA in strains containing
bbq67 (Fig. ; see Fig. S3 and S4 in the supplemental material). RsaI cleavage is blocked by adenine methylation of a GT
AC restriction site or by overlapping cytosine methylation of
CGTAC or GTA
CG sequences (
43). The latter sequences, however, are not among the protected RsaI sites present on
B. burgdorferi shuttle vectors (Fig. ). We also have not observed any protection of RsaI sites on shuttle vector DNA following
in vitro CpG or GpC methylation (data not shown). Rocha et al. (
49) suggested that GTAC could be a potential restriction site for endonucleases in
Borrelia based on the palindromic avoidance of this sequence in the
B. burgdorferi genome. The cleavage by RsaI of unprotected GTAC sites on
B. burgdorferi DNA indicates that GTAC cannot comprise the entire sequence recognized and modified by the
bbq67 gene product. A shared feature of the sites protected from RsaI cleavage is a TCGA sequence five bases upstream of the GTAC RsaI site, which suggests that it may contribute to the recognition sequence (Fig. ). We found no evidence that any of the other potential restriction sites mentioned by Rocha and colleagues (
49) could be blocked by BBQ67 or BBE02 modification (data not shown). Modification of GTAC sequences has been reported in another spirochete,
Leptospira (
46).
A search of the NCBI database indicated that all B. burgdorferi strains in the database have a homolog of bbe02, whereas only a subset of these strains have a homolog of bbq67. REBASE, the database for genes that encode R-M enzymes, identified homologs of a cytosine methyltransferase (GCGC) in available B. burgdorferi sequences, including strains 64b, 118, and 156a. B. burgdorferi strain N40 also carries this homolog (S. Casjens, personal communication). No homologs of this methyltransferase are present in B. burgdorferi B31. NCBI databases indicate that those strains possessing a bbq67 homolog do not possess a homolog of the GCGC methyltransferase and vice versa. HhaI, which cleaves GCGC but is blocked by methylation of this sequence, was efficient in digesting DNA from strain B31 but not from strain N40 (data not shown), indicating the occurrence of two different R-M systems in these B. burgdorferi strains.
Modification of endogenous DNA permits bacteria to distinguish self from nonself and block introduction of foreign DNA (
5). A question arises regarding the selective advantage of different R-M genes carried by individuals within the same species. Jeltsch has proposed that diversification of R-M systems within a species leads to different biotypes, which then leads to adaptation to different ecological niches (
24). This protection from foreign DNA, albeit from within the same species, would avoid dilution of advantageous traits for a particular niche that arose after diversifying from a common ancestor. The array of mammalian and avian hosts that
B. burgdorferi sensu lato infects could represent distinct biological niches, but currently there are no data linking particular plasmids (and the R-M systems they carry) with the adaptation of a strain to a specific host. The presence of multiple plasmid-borne R-M systems in the segmented genome of
B. burgdorferi, the heterogeneity in plasmid content among
B. burgdorferi strains, and the horizontal transfer of plasmid sequences between spirochetes, all need to be taken into account when considering the role of R-M systems in
B. burgdorferi.
Transformation studies in other bacteria have clearly shown that endogenously modified plasmid DNA can more efficiently transform strains possessing the same R-M system (
2,
3,
15,
37,
39). The data provided by this study, as well as previous work by other investigators, indicate that a primary limitation to transformation of
B. burgdorferi with shuttle vector DNA from
E. coli is the barrier presented by the endogenous plasmid-borne R-M systems of
B. burgdorferi. This transformation barrier can be overcome by inactivation of the R-M loci or loss of the plasmids that carry these loci in recipient
B. burgdorferi, which can complicate or impede subsequent studies. Obtaining plasmid DNA from an
E. coli host that also expresses the R-M genes of the recipient bacterium is an approach that has been used to introduce foreign DNA into other bacteria (
58). Our results clearly show that for such an approach to work in wild-type
B. burgdorferi strain B31 clone, both the
bbe02 and
bbq67 loci need to be taken into account. Identification of the sequences recognized by BBE02 and BBQ67 in modifying and cleaving DNA would enhance genetic manipulation of infectious
B. burgdorferi clones in the B31 strain background. Future work to address the R-M properties of BBE02 and BBQ67 include overexpressing and purifying recombinant forms of these proteins. Coupled with genome sequences, similar studies with the putative R-M genes in other
B. burgdorferi strains are needed to advance genetic studies in potentially different biotypes of this pathogenic spirochete and to explore the biological significance of gene transfer between strains.