Prokaryote restriction modification (RM) systems serve to protect bacteria from potentially detrimental foreign DNA. Recent evidence suggests that DNA methylation by the methyltransferase (MTase) components of RM systems can also have effects on transcriptome profiles. The type strain of the causative agent of Lyme disease, Borrelia burgdorferi B31, possesses two RM systems with N6-methyladenosine (m6A) MTase activity, which are encoded by the bbe02 gene located on linear plasmid lp25 and bbq67 on lp56. The specific recognition and/or methylation sequences had not been identified for either of these B. burgdorferi MTases, and it was not previously known whether these RM systems influence transcript levels. In the current study, single-molecule real-time sequencing was utilized to map genome-wide m6A sites and to identify consensus modified motifs in wild-type B. burgdorferi as well as MTase mutants lacking either the bbe02 gene alone or both bbe02 and bbq67 genes. Four novel conserved m6A motifs were identified and were fully attributable to the presence of specific MTases. Whole-genome transcriptome changes were observed in conjunction with the loss of MTase enzymes, indicating that DNA methylation by the RM systems has effects on gene expression. Genes with altered transcription in MTase mutants include those involved in vertebrate host colonization (e.g., rpoS regulon) and acquisition by/transmission from the tick vector (e.g., rrp1 and pdeB). The results of this study provide a comprehensive view of the DNA methylation pattern in B. burgdorferi, and the accompanying gene expression profiles add to the emerging body of research on RM systems and gene regulation in bacteria.

Document Type


Publication Date


Notes/Citation Information

Published in Journal of Bacteriology, v. 200, issue 24, e00395-18, p. 1-17.

Copyright © 2018 American Society for Microbiology. All Rights Reserved.

The copyright holder has granted the permission for posting the article here.

Digital Object Identifier (DOI)


Funding Information

This work was funded by NIH/NIGMS P20 GM113123-01 and NIH/NIGMS P20GM104360-01 to C. A. Brissette and NIH R03-AI113648 to B. Stevenson.

Related Content

Sequences have been deposited in the NCBI GEO sequence read archive database under accession number GSE115308.

Supplemental material for this article may be found at https://doi.org/10.1128/JB.00395-18.

zjb999094930s1.xlsx (641 kB)

zjb999094930s2.pdf (143 kB)