Rapid identification of methylase specificity (RIMS-seq) jointly identifies methylated motifs and generates shotgun sequencing of bacterial genomes

Simultaneous assessment of human genome and methylome data in a single experiment using limited deamination of methylated cytosine

Multiomics require concerted recording of independent information, ideally from a single experiment. In this study, we introduce RIMS-seq2, a high-throughput technique to simultaneously sequence genomes and overlay methylation information while requiring only a small modification of the experimental protocol for high-throughput DNA sequencing to include a controlled deamination step. Importantly, the rate of deamination of 5-methylcytosine is negligible and thus does not interfere with standard DNA sequencing and data processing. Thus, RIMS-seq2 libraries from whole- or targeted-genome sequencing show the same germline variation calling accuracy and sensitivity compared with standard DNA-seq. Additionally, regional methylation levels provide an accurate map of the human methylome.

Proxi-RIMS-seq2 applied to native microbiomes uncovers hundreds of known and novel m5C methyltransferase specificities

Methylation patterns in bacteria can be used to study Restriction-Modification (RM) or other defense systems with novel properties. While m4C and m6A methylation is well characterized mainly through PacBio sequencing, the landscape of m5C methylation is under-characterized. To bridge this gap, we performed RIMS-seq2 on microbiomes composed of resolved assemblies of distinct genomes through proximity ligation. This high-throughput approach enables the identification of m5C methylated motifs and links them to cognate methyltransferases directly on native microbiomes without the need to isolate bacterial strains. Methylation patterns can also be identified on viral DNA and compared to host DNA, strengthening evidence for virus-host interaction. Applied to three different microbiomes, the method unveils over 1900 motifs that were deposited in REBASE. The motifs include a novel 8-base recognition site (CATm5CGATG) that was experimentally validated by characterizing its cognate methyltransferase. Our findings suggest that microbiomes harbor arrays of untapped m5C methyltransferase specificities, providing insights to bacterial biology and biotechnological applications.

Single-base resolution quantitative genome methylation analysis in the model bacterium Helicobacter pylori by enzymatic methyl sequencing (EM-Seq) reveals influence of strain, growth phase, and methyl homeostasis

BACKGROUND

Bacterial epigenetics is a rapidly expanding research field. DNA methylation by diverse bacterial methyltransferases (MTases) contributes to genomic integrity and replication, and many recent studies extended MTase function also to global transcript regulation and phenotypic variation. Helicobacter pylori is currently one of those bacterial species which possess the highest number and the most variably expressed set of DNA MTases. Next-generation sequencing technologies can directly detect DNA base methylation. However, they still have limitations in their quantitative and qualitative performance, in particular for cytosine methylation.

RESULTS

As a complementing approach, we used enzymatic methyl sequencing (EM-Seq), a technology recently established that has not yet been fully evaluated for bacteria. Thereby, we assessed quantitatively, at single-base resolution, whole genome cytosine methylation for all methylated cytosine motifs in two different H. pylori strains and isogenic MTase mutants. EM-Seq reliably detected both m5C and m4C methylation. We demonstrated that three different active cytosine MTases in H. pylori provide considerably different levels of average genome-wide single-base methylation, in contrast to isogenic mutants which completely lost specific motif methylation. We found that strain identity and changed environmental conditions, such as growth phase and interference with methyl donor homeostasis, significantly influenced quantitative global and local genome-wide methylation in H. pylori at specific motifs. We also identified significantly hyper- or hypo-methylated cytosines, partially linked to overlapping MTase target motifs. Notably, we revealed differentially methylated cytosines in genome-wide coding regions under conditions of methionine depletion, which can be linked to transcript regulation.

CONCLUSIONS

This study offers new knowledge on H. pylori global and local genome-wide methylation and establishes EM-Seq for quantitative single-site resolution analyses of bacterial cytosine methylation.

A Survey of Archaeal Restriction-Modification Systems

When compared with bacteria, relatively little is known about the restriction-modification (RM) systems of archaea, particularly those in taxa outside of the haloarchaea. To improve our understanding of archaeal RM systems, we surveyed REBASE, the restriction enzyme database, to catalog what is known about the genes and activities present in the 519 completely sequenced archaeal genomes currently deposited there. For 49 (9.4%) of these genomes, we also have methylome data from Single-Molecule Real-Time (SMRT) sequencing that reveal the target recognition sites of the active m6A and m4C DNA methyltransferases (MTases). The gene-finding pipeline employed by REBASE is trained primarily on bacterial examples and so will look for similar genes in archaea. Nonetheless, the organizational structure and protein sequence of RM systems from archaea are highly similar to those of bacteria, with both groups acquiring systems from a shared genetic pool through horizontal gene transfer. As in bacteria, we observe numerous examples of "persistent" DNA MTases conserved within archaeal taxa at different levels. We experimentally validated two homologous members of one of the largest "persistent" MTase groups, revealing that methylation of C(m5C)WGG sites may play a key epigenetic role in Crenarchaea. Throughout the archaea, genes encoding m6A, m4C, and m5C DNA MTases, respectively, occur in approximately the ratio 4:2:1.

Shaping eukaryotic epigenetic systems by horizontal gene transfer

DNA methylation constitutes one of the pillars of epigenetics, relying on covalent bonds for addition and/or removal of chemically distinct marks within the major groove of the double helix. DNA methyltransferases, enzymes which introduce methyl marks, initially evolved in prokaryotes as components of restriction-modification systems protecting host genomes from bacteriophages and other invading foreign DNA. In early eukaryotic evolution, DNA methyltransferases were horizontally transferred from bacteria into eukaryotes several times and independently co-opted into epigenetic regulatory systems, primarily via establishing connections with the chromatin environment. While C5-methylcytosine is the cornerstone of plant and animal epigenetics and has been investigated in much detail, the epigenetic role of other methylated bases is less clear. The recent addition of N4-methylcytosine of bacterial origin as a metazoan DNA modification highlights the prerequisites for foreign gene co-option into the host regulatory networks, and challenges the existing paradigms concerning the origin and evolution of eukaryotic regulatory systems.

Complete genome assembly and methylome dissection of Methanococcus aeolicus PL15/H p

Although restriction-modification systems are found in both Eubacterial and Archaeal kingdoms, comparatively less is known about patterns of DNA methylation and genome defense systems in archaea. Here we report the complete closed genome sequence and methylome analysis of Methanococcus aeolicus PL15/H p , a strain of the CO2-reducing methanogenic archaeon and a commercial source for MaeI, MaeII, and MaeIII restriction endonucleases. The M. aeolicus PL15/H p genome consists of a 1.68 megabase circular chromosome predicted to contain 1,615 protein coding genes and 38 tRNAs. A combination of methylome sequencing, homology-based genome annotation, and recombinant gene expression identified five restriction-modification systems encoded by this organism, including the methyltransferase and site-specific endonuclease of MaeIII. The MaeIII restriction endonuclease was recombinantly expressed, purified and shown to have site-specific DNA cleavage activity in vitro.

Loss of p53 Concurrent with RAS and TERT Activation Induces Glioma Formation

There is an ongoing debate regarding whether gliomas originate due to functional or genetic changes in neural stem cells (NSCs). Genetic engineering has made it possible to use NSCs to establish glioma models with the pathological features of human tumors. Here, we found that RAS, TERT, and p53 mutations or abnormal expression were associated with the occurrence of glioma in the mouse tumor transplantation model. Moreover, EZH2 palmitoylation mediated by ZDHHC5 played a significant role in this malignant transformation. EZH2 palmitoylation activates H3K27me3, which in turn decreases miR-1275, increases glial fibrillary acidic protein (GFAP) expression, and weakens the binding of DNA methyltransferase 3A (DNMT3A) to the OCT4 promoter region. Thus, these findings are significant because RAS, TERT, and p53 oncogenes in human neural stem cells are conducive to a fully malignant and rapid transformation, suggesting that gene changes and specific combinations of susceptible cell types are important factors in determining the occurrence of gliomas.

Tuning of Gene Expression in Clostridium phytofermentans Using Synthetic Promoters and CRISPRi

Control of gene expression is fundamental to cell engineering. Here we demonstrate a set of approaches to tune gene expression in Clostridia using the model Clostridium phytofermentans. Initially, we develop a simple benchtop electroporation method that we use to identify a set of replicating plasmids and resistance markers that can be cotransformed into C. phytofermentans. We define a series of promoters spanning a >100-fold expression range by testing a promoter library driving the expression of a luminescent reporter. By insertion of tet operator sites upstream of the reporter, its expression can be quantitatively altered using the Tet repressor and anhydrotetracycline (aTc). We integrate these methods into an aTc-regulated dCas12a system with which we show in vivo CRISPRi-mediated repression of reporter and fermentation genes in C. phytofermentans. Together, these approaches advance genetic transformation and experimental control of gene expression in Clostridia.

Methyl-SNP-seq reveals dual readouts of methylome and variome at molecule resolution while enabling target enrichment

Covalent modifications of genomic DNA are crucial for most organisms to survive. Amplicon-based high-throughput sequencing technologies erase all DNA modifications to retain only sequence information for the four canonical nucleobases, necessitating specialized technologies for ascertaining epigenetic information. To also capture base modification information, we developed Methyl-SNP-seq, a technology that takes advantage of the complementarity of the double helix to extract the methylation and original sequence information from a single DNA molecule. More specifically, Methyl-SNP-seq uses bisulfite conversion of one of the strands to identify cytosine methylation while retaining the original four-bases sequence information on the other strand. As both strands are locked together to link the dual readouts on a single paired-end read, Methyl-SNP-seq allows detecting the methylation status of any DNA even without a reference genome. Because one of the strands retains the original four nucleotide composition, Methyl-SNP-seq can also be used in conjunction with standard sequence-specific probes for targeted enrichment and amplification. We show the usefulness of this technology in a broad spectrum of applications ranging from allele-specific methylation analysis in humans to identification of methyltransferase specificity in complex bacterial communities.