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.

A comprehensive mechanism for 5-carboxylcytosine-induced transcriptional pausing revealed by Markov state models.

RNA polymerase II (Pol II) surveils the genome, pausing as it encounters DNA lesions and base modifications and initiating signals for DNA repair among other important regulatory events. Recent work suggests that Pol II pauses at 5-carboxycytosine (5caC), an epigenetic modification of cytosine, because of a specific hydrogen bond between the carboxyl group of 5caC and a specific residue in fork loop 3 of Pol II. This hydrogen bond compromises productive NTP binding and slows down elongation. Apart from this specific interaction, the carboxyl group of 5caC can potentially interact with numerous charged residues in the cleft of Pol II. However, it is not clear how other interactions between Pol II and 5caC contribute to pausing. In this study, we use Markov state models (a type of kinetic network models) built from extensive molecular dynamics simulations to comprehensively study the impact of 5caC on Pol II translocation. We describe two translocation intermediates with specific interactions that prevent the template base from loading into the Pol II active site. In addition to the previously observed state with 5caC constrained by fork loop 3, we discovered a new intermediate state with a hydrogen bond between 5caC and fork loop 2. Surprisingly, we find that 5caC may curb translocation by suppressing kinking of the helix bordering the active site (the bridge helix) because its high flexibility is critical to translocation. Our work provides new insights into how epigenetic modifications of genomic DNA can modulate Pol II translocation, inducing pauses in transcription.

Evaluation and minimization of Cas9-independent off-target DNA editing by cytosine base editors.

Cytosine base editors (CBEs) enable targeted C•G-to-T•A conversions in genomic DNA. Recent studies report that BE3, the original CBE, induces a low frequency of genome-wide Cas9-independent off-target C•G-to-T•A mutation in mouse embryos and in rice. Here we develop multiple rapid, cost-effective methods to screen the propensity of different CBEs to induce Cas9-independent deamination in Escherichia coli and in human cells. We use these assays to identify CBEs with reduced Cas9-independent deamination and validate via whole-genome sequencing that YE1, a narrowed-window CBE variant, displays background levels of Cas9-independent off-target editing. We engineered YE1 variants that retain the substrate-targeting scope of high-activity CBEs while maintaining minimal Cas9-independent off-target editing. The suite of CBEs characterized and engineered in this study collectively offer ~10-100-fold lower average Cas9-independent off-target DNA editing while maintaining robust on-target editing at most positions targetable by canonical CBEs, and thus are especially promising for applications in which off-target editing must be minimized.

Sequential Assessment of Multiple Epigenetic Modifications of Cytosine in Whole Genomic DNA by Surface Plasmon Resonance.

Since the discovery of the active DNA demethylation pathway in mammals, numerous efforts have been made to distinguish epigenetic cytosine variants, including 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). However, the rapid discrimination of multiple cytosine variants in DNA remains challenging because the conventional assays require time-consuming DNA pretreatments, such as enzymatical digestion and chemical conversion. Here we demonstrated the high-throughput discrimination of four cytosine variants in DNA by using a sequential surface-plasmon-resonance (SPR)-based immunochemical assay. The target DNAs were biotinylated in one step with a bifunctional linker 1 and robustly immobilized on a streptavidin-coated sensor surface to hold them in place during an alkali washing designed to remove residual antibodies. By repeating the injection of antibodies and washing, we achieved a sequential assessment of cytosine variants in identical DNA and identified the yield of in vitro 5mC oxidation in genomic DNA by the ten-eleven translocation 1 (TET1) enzyme. These results demonstrated that our sequential SPR-based immunochemical assay was effective for evaluating multiple epigenetic modifications in a whole genome with a single row operation without time-consuming DNA pretreatments.

Assessment and site-specific manipulation of DNA (hydroxy-)methylation during mouse corticogenesis.

Dynamic changes in DNA (hydroxy-)methylation are fundamental for stem cell differentiation. However, the signature of these epigenetic marks in specific cell types during corticogenesis is unknown. Moreover, site-specific manipulation of cytosine modifications is needed to reveal the significance and function of these changes. Here, we report the first assessment of (hydroxy-)methylation in neural stem cells, neurogenic progenitors, and newborn neurons during mammalian corticogenesis. We found that gain in hydroxymethylation and loss in methylation occur sequentially at specific cellular transitions during neurogenic commitment. We also found that these changes predominantly occur within enhancers of neurogenic genes up-regulated during neurogenesis and target of pioneer transcription factors. We further optimized the use of dCas9-Tet1 manipulation of (hydroxy-)methylation, locus-specifically, in vivo, showing the biological relevance of our observations for Dchs1, a regulator of corticogenesis involved in developmental malformations and cognitive impairment. Together, our data reveal the dynamics of cytosine modifications in lineage-related cell types, whereby methylation is reduced and hydroxymethylation gained during the neurogenic lineage concurrently with up-regulation of pioneer transcription factors and activation of enhancers for neurogenic genes.

Uncovering universal rules governing the selectivity of the archetypal DNA glycosylase TDG.

Thymine DNA glycosylase (TDG) is a pivotal enzyme with dual roles in both genome maintenance and epigenetic regulation. TDG is involved in cytosine demethylation at CpG sites in DNA. Here we have used molecular modeling to delineate the lesion search and DNA base interrogation mechanisms of TDG. First, we examined the capacity of TDG to interrogate not only DNA substrates with 5-carboxyl cytosine modifications but also G:T mismatches and nonmismatched (A:T) base pairs using classical and accelerated molecular dynamics. To determine the kinetics, we constructed Markov state models. Base interrogation was found to be highly stochastic and proceeded through insertion of an arginine-containing loop into the DNA minor groove to transiently disrupt Watson-Crick pairing. Next, we employed chain-of-replicas path-sampling methodologies to compute minimum free energy paths for TDG base extrusion. We identified the key intermediates imparting selectivity and determined effective free energy profiles for the lesion search and base extrusion into the TDG active site. Our results show that DNA sculpting, dynamic glycosylase interactions, and stabilizing contacts collectively provide a powerful mechanism for the detection and discrimination of modified bases and epigenetic marks in DNA.

Dynamics of the excised base release in thymine DNA glycosylase during DNA repair process.

Thymine DNA glycosylase (TDG) initiates base excision repair by cleaving the N-glycosidic bond between the sugar and target base. After catalysis, the release of excised base is a requisite step to terminate the catalytic cycle and liberate the TDG for the following enzymatic reactions. However, an atomistic-level understanding of the dynamics of the product release process in TDG remains unknown. Here, by employing molecular dynamics simulations combined with the Markov State Model, we reveal the dynamics of the thymine release after the excision at microseconds timescale and all-atom resolution. We identify several key metastable states of the thymine and its dominant releasing pathway. Notably, after replacing the TDG residue Gly142 with tyrosine, the thymine release is delayed compared to the wild-type (wt) TDG, as supported by our potential of mean force (PMF) calculations. These findings warrant further experimental tests to potentially trap the excised base in the active site of TDG after the catalysis, which had been unsuccessful by previous attempts. Finally, we extended our studies to other TDG products, including the uracil, 5hmU, 5fC and 5caC bases in order to compare the product release for different targeting bases in the TDG-DNA complex.

Detecting DNA cytosine methylation using nanopore sequencing.

In nanopore sequencing devices, electrolytic current signals are sensitive to base modifications, such as 5-methylcytosine (5-mC). Here we quantified the strength of this effect for the Oxford Nanopore Technologies MinION sequencer. By using synthetically methylated DNA, we were able to train a hidden Markov model to distinguish 5-mC from unmethylated cytosine. We applied our method to sequence the methylome of human DNA, without requiring special steps for library preparation.

Detection of differentially methylated regions from bisulfite-seq data by hidden Markov models incorporating genome-wide methylation level distributions.

BACKGROUND: Detection of differential methylation between biological samples is an important task in bisulfite-seq data analysis. Several studies have attempted de novo finding of differentially methylated regions (DMRs) using hidden Markov models (HMMs). However, there is room for improvement in the design of HMMs, especially on emission functions that evaluate the likelihood of differential methylation at each cytosine site. RESULTS: We describe a new HMM for DMR detection from bisulfite-seq data. Our method utilizes emission functions that combine binomial models for aligned read counts, and beta mixtures for incorporating genome-wide methylation level distributions. We also develop unsupervised learning algorithms to adjust parameters of the beta-binomial models depending on differential methylation types (up, down, and not changed). In experiments on both simulated and real datasets, the new HMM improves DMR detection accuracy compared with HMMs in our previous study. Furthermore, our method achieves better accuracy than other methods using Fisher's exact test and methylation level smoothing. CONCLUSIONS: Our method enables accurate DMR detection from bisulfite-seq data. The implementation of our method is named ComMet, and distributed as a part of Bisulfighter package, which is available at http://epigenome.cbrc.jp/bisulfighter.

JBP1-seq: a fast and efficient method for genome-wide profiling of 5hmC.

We developed a novel approach, J-binding protein 1 sequencing (JBP1-seq), that combines the benefits of an improved recombinant JBP1 protein, Nextera-based library construction, and next-generation sequencing (NGS) for genome-wide profiling of 5-hydroxymethylcytosine (5hmC). Compared with the original JBP1, this new recombinant JBP1 was biotinylated in vivo and conjugated to magnetic beads via biotin-streptavidin interactions. These modifications allowed a more efficient and consistent pull-down of ß-glucosyl-5-hydroxymethylcytosine (ß-glu-5hmC), and sequence-ready libraries can be generated within 4.5h from DNA inputs as low as 50ng. 5hmC enrichment of human brain DNA using the new JBP1 resulted in over 25,000 peaks called, which is significantly higher than the 4003 peaks enriched using the old JBP1. Comparison of the technical duplicates and validations with other platforms indicated the results are reproducible and reliable. Thus, JBP1-seq provides a fast, efficient, and cost-effective method for accurate 5hmC genome-wide profiling.

Application of a low cost array-based technique - TAB-Array - for quantifying and mapping both 5mC and 5hmC at single base resolution in human pluripotent stem cells.

5-hydroxymethylcytosine (5hmC), an oxidized derivative of 5-methylcytosine (5mC), has been implicated as an important epigenetic regulator of mammalian development. Current procedures use DNA sequencing methods to discriminate 5hmC from 5mC, limiting their accessibility to the scientific community. Here we report a method that combines TET-assisted bisulfite conversion with Illumina 450K DNA methylation arrays for a low-cost high-throughput approach that distinguishes 5hmC and 5mC signals at base resolution. Implementing this approach, termed "TAB-array", we assessed DNA methylation dynamics in the differentiation of human pluripotent stem cells into cardiovascular progenitors and neural precursor cells. With the ability to discriminate 5mC and 5hmC, we identified a large number of novel dynamically methylated genomic regions that are implicated in the development of these lineages. The increased resolution and accuracy afforded by this approach provides a powerful means to investigate the distinct contributions of 5mC and 5hmC in human development and disease.

Cytosine unstacking and strand slippage at an insertion-deletion mutation sequence in an overhang-containing DNA duplex.

Base unstacking in template strands, when accompanied by strand slippage, can result in deletion mutations during strand extension by nucleic acid polymerases. In a GCCC mutation hot-spot sequence, which was previously identified to have a 50% probability of causing such mutations during DNA replication by a Y-family polymerase, a single-base deletion mutation could result from such unstacking of any one of its three template cytosines. In this study, the intrinsic energetic differences in unstacking among these three cytosines in a solvated DNA duplex overhang model were examined using umbrella sampling molecular dynamics simulations. The free energy profiles obtained show that cytosine unstacking grows progressively more unfavorable as one moves inside the duplex from the 5'-end of the overhang template strand. Spontaneous strand slippage occurs in response to such base unstacking in the direction of both the major and minor grooves for all three cytosines. Unrestrained simulations run from three distinct strand-slipped states and one non-strand-slipped state suggest that a more duplexlike environment can help stabilize strand slippage. The possible underlying reasons and biological implications of these observations are discussed in the context of nucleic acid replication active site dynamics.

Direct and indirect organogenesis of Clivia miniata and assessment of DNA methylation changes in various regenerated plantlets.

UNLABELLED: Clivia miniata is an important indoor ornamental plant and has been reported to have medicinal value. We developed an efficient in vitro micropropagation protocol from young leaves (indirect organogenesis), young petals (indirect organogenesis) and shoot tips (direct organogenesis) of this plant. Using young leaves and shoot tips as explants, the regeneration frequencies were much higher than those in previous investigation and the regeneration was dependent upon less nutrition. We speculated that the leaf-derived callus can generate amino acids necessary for protein synthesis by itself. We employed the methylation-sensitive amplified polymorphism (MSAP) method to assess cytosine methylation variation in various regenerated plantlets and between organs. The MSAP profiles indicated that the frequency of somaclonal variation in the form of cytosine methylation was highest in petal-derived plantlets followed by secondary leaf-derived, primary leaf-derived and shoot tip-derived plantlets, but the methylation variation in petal-derived plantlets was lower than between petals and leaves of a single plant. The results indicated that the methylation variation in regenerated plantlets was related to the types of explants, regeneration pathways and number of regeneration generations. Two possible factors for the highest somaclonal variation rate in petal-derived plantlets are the callus phase and petal-specific set of epigenetic regulators. The property of meristem integrity can account for the lowest variation rate in shoot tip-derived plantlets. Moreover, the secondary plantlets underwent a longer total period of in vitro culture, which can explain why the methylation variation rate in the secondary plantlets is higher than in the primary ones. KEY MESSAGE: Methylation variation in regenerated plantlets of C. miniata was found to be related to the types of explants, regeneration pathways and number of regeneration generations.

Quantitative measurement of genome-wide DNA methylation by a reliable and cost-efficient enzyme-linked immunosorbent assay technique.

DNA methylation, the conversion of cytosine to 5-methylcytosine, is an important epigenetic modification involved in gene regulation. DNA methylation is essential for normal development whereas abnormal methylation has been implicated in pathological conditions including cancer. To evaluate the extent and variation of genome-wide DNA methylation and its changes during cellular differentiation and tumorgenesis as well as the interplay with histone modifications, accurate and reproducible quantification of the genomic DNA methylation level is required. These measurements have so far been achieved only by sophisticated and costly techniques. Here we report the generation of an enzyme-linked immunosorbent assay (methDNA-ELISA) for the accurate quantification of global DNA methylation levels. The linear region of this methDNA-ELISA ranges from 1 to 10%, making it highly suitable for the typical ranges from 2 to 6% in mammalian genomes. This method requires 10 ng of isolated DNA per sample, thus permitting investigation with minimal amounts of DNA previously not applicable for global DNA methylation analysis, e.g., clinical biopsies or cells collected by microdissection.

Measuring topology of low-intensity DNA methylation sites for high-throughput assessment of epigenetic drug-induced effects in cancer cells.

Epigenetic anti-cancer drugs with demethylating effects have shown to alter genome organization in mammalian cell nuclei. The interest in the development of novel epigenetic drugs has increased the demand for cell-based assays to evaluate drug performance in pre-clinical studies. An imaging-based cytometrical approach that can measure demethylation effects as changes in the spatial nuclear distributions of methylated cytosine and global DNA in cancer cells is introduced in this paper. The cells were studied by immunofluorescence with a specific antibody against 5-methylcytosine (MeC), and 4,6-diamidino-2-phenylindole (DAPI) for delineation of methylated sites and global DNA in nuclei. In the preprocessing step the segmentation of nuclei in three-dimensional images (3-D) is followed by an automated assessment of nuclear DAPI/MeC patterns to exclude dissimilar entities. Next, low-intensity MeC (LIM) and low-intensity DNA (LID) sites of similar nuclei are localized and processed to obtain specific nuclear density profiles. These profiles sampled at half of the total nuclear volume yielded two parameters: LIM(0.5) and LID(0.5). The analysis shows that zebularine and 5-azacytidine-the two tested epigenetic drugs introduce changes in the spatial distribution of low-intensity DNA and MeC signals. LIM(0.5) and LID(0.5) were significantly different (p<0.001) in 5-azacytidine treated (n=660) and zebularine treated (n=496) vs. untreated (n=649) DU145 human prostate cancer cells. In the latter case the LIM sites were predominantly found at the nuclear border, whereas treated populations showed different degrees of increase in LIMs towards the interior nuclear space, in which a large portion of heterochromatin is located. The cell-by-cell evaluation of changes in the spatial reorganization of MeC/DAPI signals revealed that zebularine is a more gentle demethylating agent than 5-azacytidine. Measuring changes in the topology of low-intensity sites can potentially be a valuable component in the high-throughput assessment of demethylation and risk of chromatin reorganization in epigenetic-drug screening tasks.

DNA mismatch correction by Very Short Patch repair may have altered the abundance of oligonucleotides in the E. coli genome.

A base mismatch correction process in E. coli K-12 called Very Short Patch (VSP) repair corrects T:G mismatches to C:G when found in certain sequence contexts. Two of the substrate mismatches (5'-CTWGG/3'-GGW'CC; W = A or T) occur in the context of cytosine methylation in DNA and reduce the mutagenic effects of 5-methylcytosine deamination to thymine. However, VSP repair is also known to repair T:G mismatches that are not expected to arise from 5-methylcytosine deamination (example--CTAG/GGT-C). In these cases, if the original base pair were a T:A, VSP repair would cause a T to C transition. We have carried out Markov chain analysis of an E. coli sequence database to determine if repair at the latter class of sites has altered the abundance of the relevant tetranucleotides. The results are consistent with the prediction that VSP repair would tend to deplete the genome of the 'T' containing sequences (example--CTAG), while enriching it for the corresponding 'C' containing sequences (CCAG). Further, they provide an explanation for the known scarcity of CTAG containing restriction enzyme sites among the genomes of enteric bacteria and identify VSP repair as a force in shaping the sequence composition of bacterial genomes.