One-pot trimodal mapping of unmethylated, hydroxymethylated, and open chromatin sites unveils distinctive 5hmC roles at dynamic chromatin loci.

We present a method, named Mx-TOP, for profiling of three epigenetic regulatory layers-chromatin accessibility, general DNA modification, and DNA hydroxymethylation-from a single library. The approach is based on chemo-enzymatic covalent tagging of unmodified CG sites and hydroxymethylated cytosine (5hmC) along with GC sites in chromatin, which are then mapped using tag-selective base-resolution TOP-seq sequencing. Our in-depth validation of the approach revealed its sensitivity and informativity in evaluating chromatin accessibility and DNA modification interactions that drive transcriptional regulation. We employed the technology in a study of chromatin and DNA demethylation dynamics during in vitro neuronal differentiation. The study highlighted the involvement of gene body 5hmC in modulating an extensive decoupling between promoter accessibility and transcription. The importance of 5hmC in chromatin remodeling was further demonstrated by the observed resistance of the developmentally acquired open loci to the global 5hmC erasure in neuronal progenitors.

Scikick: A sidekick for workflow clarity and reproducibility during extensive data analysis.

Reproducibility is crucial for scientific progress, yet a clear research data analysis workflow is challenging to implement and maintain. As a result, a record of computational steps performed on the data to arrive at the key research findings is often missing. We developed Scikick, a tool that eases the configuration, execution, and presentation of scientific computational analyses. Scikick allows for workflow configurations with notebooks as the units of execution, defines a standard structure for the project, automatically tracks the defined interdependencies between the data analysis steps, and implements methods to compile all research results into a cohesive final report. Utilities provided by Scikick help turn the complicated management of transparent data analysis workflows into a standardized and feasible practice. Scikick version 0.2.1 code and documentation is available as supplementary material. The Scikick software is available on GitHub (https://github.com/matthewcarlucci/scikick) and is distributed with PyPi (https://pypi.org/project/scikick/) under a GPL-3 license.

Distribution and regulatory roles of oxidized 5-methylcytosines in DNA and RNA of the basidiomycete fungi Laccaria bicolor and Coprinopsis cinerea.

The formation of three oxidative DNA 5-methylcytosine (5mC) modifications (oxi-mCs)-5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC)-by the TET/JBP family of dioxygenases prompted intensive studies of their functional roles in mammalian cells. However, the functional interplay of these less abundant modified nucleotides in other eukaryotic lineages remains poorly understood. We carried out a systematic study of the content and distribution of oxi-mCs in the DNA and RNA of the basidiomycetes Laccaria bicolor and Coprinopsis cinerea, which are established models to study DNA methylation and developmental and symbiotic processes. Quantitative liquid chromatography-tandem mass spectrometry revealed persistent but uneven occurrences of 5hmC, 5fC and 5caC in the DNA and RNA of the two organisms, which could be upregulated by vitamin C. 5caC in RNA (5carC) was predominantly found in non-ribosomal RNA, which potentially includes non-coding, messenger and small RNA species. Genome-wide mapping of 5hmC and 5fC using the single CG analysis techniques hmTOP-seq and foTOP-seq pointed at involvement of oxi-mCs in the regulation of gene expression and silencing of transposable elements. The implicated diverse roles of 5mC and oxi-mCs in the two fungi highlight the epigenetic importance of the latter modifications, which are often neglected in standard whole-genome bisulfite analyses.

Selective chemical tracking of Dnmt1 catalytic activity in live cells.

Enzymatic methylation of cytosine to 5-methylcytosine in DNA is a fundamental epigenetic mechanism involved in mammalian development and disease. DNA methylation is brought about by collective action of three AdoMet-dependent DNA methyltransferases, whose catalytic interactions and temporal interplay are poorly understood. We used structure-guided engineering of the Dnmt1 methyltransferase to enable catalytic transfer of azide tags onto DNA from a synthetic cofactor analog, Ado-6-azide, in vitro. We then CRISPR-edited the Dnmt1 locus in mouse embryonic stem cells to install the engineered codon, which, following pulse internalization of the Ado-6-azide cofactor by electroporation, permitted selective azide tagging of Dnmt1-specific genomic targets in cellulo. The deposited covalent tags were exploited as "click" handles for reading adjoining sequences and precise genomic mapping of the methylation sites. The proposed approach, Dnmt-TOP-seq, enables high-resolution temporal tracking of the Dnmt1 catalysis in mammalian cells, paving the way to selective studies of other methylation pathways in eukaryotic systems.

Multiomics Analysis of Neuroblastoma Cells Reveals a Diversity of Malignant Transformations.

Neuroblastoma (NB) is a pediatric cancer of the developing sympathetic nervous system that exhibits significant variation in the stage of differentiation and cell composition of tumors. Global loss of DNA methylation and genomic 5-hydroxymethylcytosine (5hmC) is a hallmark of human cancers. Here, we used our recently developed single-base resolution approaches, hmTOP-seq and uTOP-seq, for construction of 5hmC maps and identification of large partially methylated domains (PMDs) in different NB cell subpopulations. The 5hmC profiles revealed distinct signatures characteristic to different cell lineages and stages of malignant transformation of NB cells in a conventional and oxygen-depleted environment, which often occurs in tumors. The analysis of the cell-type-specific PMD distribution highlighted differences in global genome organization among NB cells that were ascribed to the same lineage identity by transcriptomic networks. Collectively, we demonstrated a high informativeness of the integrative epigenomic and transcriptomic research and large-scale genome structure in investigating the mechanisms that regulate cell identities and developmental stages of NB cells. Such multiomics analysis, as compared with mutational studies, open new ways for identification of novel disease-associated features which bring prognostic and therapeutic value in treating this aggressive pediatric disease.

Identification of fetal unmodified and 5-hydroxymethylated CG sites in maternal cell-free DNA for non-invasive prenatal testing.

BACKGROUND: Massively parallel sequencing of maternal cell-free DNA (cfDNA) is widely used to test fetal genetic abnormalities in non-invasive prenatal testing (NIPT). However, sequencing-based approaches are still of high cost. Building upon previous knowledge that placenta, the main source of fetal circulating DNA, is hypomethylated in comparison to maternal tissue counterparts of cfDNA, we propose that targeting either unmodified or 5-hydroxymethylated CG sites specifically enriches fetal genetic material and reduces numbers of required analytical sequencing reads thereby decreasing cost of a test. METHODS: We employed uTOPseq and hmTOP-seq approaches which combine covalent derivatization of unmodified or hydroxymethylated CG sites, respectively, with next generation sequencing, or quantitative real-time PCR. RESULTS: We detected increased 5-hydroxymethylcytosine (5hmC) levels in fetal chorionic villi (CV) tissue samples as compared with peripheral blood. Using our previously developed uTOP-seq and hmTOP-seq approaches we obtained whole-genome uCG and 5hmCG maps of 10 CV tissue and 38 cfDNA samples in total. Our results indicated that, in contrast to conventional whole genome sequencing, such epigenomic analysis highly specifically enriches fetal DNA fragments from maternal cfDNA. While both our approaches yielded 100% accuracy in detecting Down syndrome in fetuses, hmTOP-seq maintained such accuracy at ultra-low sequencing depths using only one million reads. We identified 2164 and 1589 placenta-specific differentially modified and 5-hydroxymethylated regions, respectively, in chromosome 21, as well as 3490 and 2002 Down syndrome-specific differentially modified and 5-hydroxymethylated regions, respectively, that can be used as biomarkers for identification of Down syndrome or other epigenetic diseases of a fetus. CONCLUSIONS: uTOP-seq and hmTOP-seq approaches provide a cost-efficient and sensitive epigenetic analysis of fetal abnormalities in maternal cfDNA. The results demonstrated that T21 fetuses contain a perturbed epigenome and also indicated that fetal cfDNA might originate from fetal tissues other than placental chorionic villi. Robust covalent derivatization followed by targeted analysis of fetal DNA by sequencing or qPCR presents an attractive strategy that could help achieve superior sensitivity and specificity in prenatal diagnostics.

A Bisulfite-free Approach for Base-Resolution Analysis of Genomic 5-Carboxylcytosine.

Due to an extreme rarity of 5-carboxylcytosine (5caC) in the mammalian genome, investigation of its role brings a considerable challenge. Methods based on bisulfite sequencing have been proposed for genome-wide 5caC analysis. However, bisulfite-based sequencing of scarcely abundant 5caC demands significant experimental and computational resources, increasing sequencing cost. Here, we present a bisulfite-free approach, caCLEAR, for high-resolution mapping of 5caCGs. The method uses an atypical activity of the methyltransferase eM.SssI to remove a carboxyl group from 5caC, generating unmodified CGs, which are localized by uTOP-seq sequencing. Validation of caCLEAR on model DNA systems and mouse ESCs supports the suitability of caCLEAR for analysis of 5caCGs. The 5caCG profiles of naive and primed pluripotent ESCs reflect their distinct demethylation dynamics and demonstrate an association of 5caC with gene expression. Generally, we demonstrate that caCLEAR is a robust economical approach that could help provide deeper insights into biological roles of 5caC.

Precise genomic mapping of 5-hydroxymethylcytosine via covalent tether-directed sequencing.

5-hydroxymethylcytosine (5hmC) is the most prevalent intermediate on the oxidative DNA demethylation pathway and is implicated in regulation of embryogenesis, neurological processes, and cancerogenesis. Profiling of this relatively scarce genomic modification in clinical samples requires cost-effective high-resolution techniques that avoid harsh chemical treatment. Here, we present a bisulfite-free approach for 5hmC profiling at single-nucleotide resolution, named hmTOP-seq (5hmC-specific tethered oligonucleotide-primed sequencing), which is based on direct sequence readout primed at covalently labeled 5hmC sites from an in situ tethered DNA oligonucleotide. Examination of distinct conjugation chemistries suggested a structural model for the tether-directed nonhomologous polymerase priming enabling theoretical evaluation of suitable tethers at the design stage. The hmTOP-seq procedure was optimized and validated on a small model genome and mouse embryonic stem cells, which allowed construction of single-nucleotide 5hmC maps reflecting subtle differences in strand-specific CG hydroxymethylation. Collectively, hmTOP-seq provides a new valuable tool for cost-effective and precise identification of 5hmC in characterizing its biological role and epigenetic changes associated with human disease.

DiscoRhythm: an easy-to-use web application and R package for discovering rhythmicity.

MOTIVATION: Biological rhythmicity is fundamental to almost all organisms on Earth and plays a key role in health and disease. Identification of oscillating signals could lead to novel biological insights, yet its investigation is impeded by the extensive computational and statistical knowledge required to perform such analysis. RESULTS: To address this issue, we present DiscoRhythm (Discovering Rhythmicity), a user-friendly application for characterizing rhythmicity in temporal biological data. DiscoRhythm is available as a web application or an R/Bioconductor package for estimating phase, amplitude, and statistical significance using four popular approaches to rhythm detection (Cosinor, JTK Cycle, ARSER, and Lomb-Scargle). We optimized these algorithms for speed, improving their execution times up to 30-fold to enable rapid analysis of -omic-scale datasets in real-time. Informative visualizations, interactive modules for quality control, dimensionality reduction, periodicity profiling, and incorporation of experimental replicates make DiscoRhythm a thorough toolkit for analyzing rhythmicity. AVAILABILITY AND IMPLEMENTATION: The DiscoRhythm R package is available on Bioconductor (https://bioconductor.org/packages/DiscoRhythm), with source code available on GitHub (https://github.com/matthewcarlucci/DiscoRhythm) under a GPL-3 license. The web application is securely deployed over HTTPS (https://disco.camh.ca) and is freely available for use worldwide. Local instances of the DiscoRhythm web application can be created using the R package or by deploying the publicly available Docker container (https://hub.docker.com/r/mcarlucci/discorhythm). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Clinical significance of miRNA host gene promoter methylation in prostate cancer.

Only a part of prostate cancer (PCa) patients has aggressive malignancy requiring adjuvant treatment after radical prostatectomy (RP). Biomarkers capable to predict biochemical PCa recurrence (BCR) after RP would significantly improve preoperative risk stratification and treatment decisions. MicroRNA (miRNA) deregulation has recently emerged as an important phenomenon in tumor development and progression, however, the mechanisms remain largely unstudied. In the present study, based on microarray profiling of DNA methylation in 9 pairs of PCa and noncancerous prostate tissues (NPT), host genes of miR-155-5p, miR-152-3p, miR-137, miR-31-5p, and miR-642a, -b were analyzed for promoter methylation in 129 PCa, 35 NPT, and 17 benign prostatic hyperplasia samples (BPH) and compared to the expression of mature miRNAs and their selected targets (DNMT1, KDM1A, and KDM5B). The Cancer Genome Atlas dataset was utilized for validation. Methylation of mir-155, mir-152, and mir-137 host genes was PCa-specific, and downregulation of miR-155-5p significantly correlated with promoter methylation. Higher KDM5B expression was observed in samples with methylated mir-155 or mir-137 promoters, whereas upregulation of KDM1A and DNMT1 was associated with mir-155 and mir-152 methylation status, respectively. Promoter methylation of mir-155, mir-152, and mir-31 was predictive of BCR-free survival in various Cox models and increased the prognostic value of clinicopathologic factors. In conclusion, methylated mir-155, mir-152, mir-137, and mir-31 host genes are promising diagnostic and/or prognostic biomarkers of PCa. Methylation status of particular miRNA host genes as independent variables or in combinations might assist physicians in identifying poor prognosis PCa patients preoperatively.

Tethered Oligonucleotide-Primed Sequencing, TOP-Seq: A High-Resolution Economical Approach for DNA Epigenome Profiling.

Modification of CG dinucleotides in DNA is part of epigenetic regulation of gene function in vertebrates and is associated with complex human disease. Bisulfite sequencing permits high-resolution analysis of cytosine modification in mammalian genomes; however, its utility is often limited due to substantial cost. Here, we describe an alternative epigenome profiling approach, named TOP-seq, which is based on covalent tagging of individual unmodified CG sites followed by non-homologous priming of the DNA polymerase action at these sites to directly produce adjoining regions for their sequencing and precise genomic mapping. Pilot TOP-seq analyses of bacterial and human genomes showed a better agreement of TOP-seq with published bisulfite sequencing maps as compared to widely used MBD-seq and MRE-seq and permitted identification of long-range and gene-level differential methylation among human tissues and neuroblastoma cell types. Altogether, we propose an affordable single CG-resolution technique well suited for large-scale epigenome studies.