The R package and a quick-start vignette is available at http://bioinfo.au.tsinghua.edu.cn/software/degseq

Single-molecule real-time (SMRT) DNA sequencing allows the systematic detection of chemical modifications such as methylation but has not previously been applied on a genome-wide scale. We used this approach to detect 49,311 putative 6-methyladenine (m6A) residues and 1,407 putative 5-methylcytosine (m5C) residues in the genome of a pathogenic Escherichia coli strain. We obtained strand-specific information for methylation sites and a quantitative assessment of the frequency of methylation at each modified position. We deduced the sequence motifs recognized by the methyltransferase enzymes present in this strain without prior knowledge of their specificity. Furthermore, we found that deletion of a phage-encoded methyltransferase-endonuclease (restriction-modification; RM) system induced global transcriptional changes and led to gene amplification, suggesting that the role of RM systems extends beyond protecting host genomes from foreign DNA.

DNA methylation is an important epigenetic mechanism for phenotypic diversification in all forms of life. We previously described remarkable cell-to-cell heterogeneity in epigenetic pattern within a clonal population of Streptococcus pneumoniae, a leading human pathogen. We here report that the epigenetic diversity is caused by extensive DNA inversions among hsdS A, hsdS B, and hsdS C, three methyltransferase hsdS genes in the Spn556II type-I restriction modification (R-M) locus. Because hsdS A encodes the sequence recognition subunit of this type-I R-M DNA methyltransferase, these site-specific recombinations generate pneumococcal cells with variable HsdSA alleles and thereby diverse genome methylation patterns. Most importantly, the DNA methylation pattern specified by the HsdSA1 allele leads to the formation of opaque colonies, whereas the pneumococci lacking HsdSA1 produce transparent colonies. Furthermore, this HsdSA1-dependent phase variation requires intact DNA methylase activity encoded by hsdM in the Spn556II (renamed colony opacity determinant or cod) locus. Thus, the DNA inversion-driven ON/OFF switch of the hsdS A1 allele in the cod locus and resulting epigenetic switch dictate the phase variation between the opaque and transparent phenotypes. Phase variation has been well documented for its importance in pneumococcal carriage and invasive infection, but its molecular basis remains unclear. Our work has discovered a novel epigenetic cause for this significant pathobiology phenomenon in S. pneumoniae. Lastly, our findings broadly represents a significant advancement in our understanding of bacterial R-M systems and their potential in shaping epigenetic and phenotypic diversity of the prokaryotic organisms because similar site-specific recombination systems widely exist in many archaeal and bacterial species.