DMRscaler: a scale-aware method to identify regions of differential DNA methylation spanning basepair to multi-megabase features

A novel approach toward optimal workflow selection for DNA methylation biomarker discovery

DNA methylation is a major epigenetic modification involved in many physiological processes. Normal methylation patterns are disrupted in many diseases and methylation-based biomarkers have shown promise in several contexts. Marker discovery typically involves the analysis of publicly available DNA methylation data from high-throughput assays. Numerous methods for identification of differentially methylated biomarkers have been developed, making the need for best practices guidelines and context-specific analyses workflows exceedingly high. To this end, here we propose TASA, a novel method for simulating methylation array data in various scenarios. We then comprehensively assess different data analysis workflows using real and simulated data and suggest optimal start-to-finish analysis workflows. Our study demonstrates that the choice of analysis pipeline for DNA methylation-based marker discovery is crucial and different across different contexts.

Comprehensive evaluation of the implementation of episignatures for diagnosis of neurodevelopmental disorders (NDDs)

Episignatures are popular tools for the diagnosis of rare neurodevelopmental disorders. They are commonly based on a set of differentially methylated CpGs used in combination with a support vector machine model. DNA methylation (DNAm) data often include missing values due to changes in data generation technology and batch effects. While many normalization methods exist for DNAm data, their impact on episignature performance have never been assessed. In addition, technologies to quantify DNAm evolve quickly and this may lead to poor transposition of existing episignatures generated on deprecated array versions to new ones. Indeed, probe removal between array versions, technologies or during preprocessing leads to missing values. Thus, the effect of missing data on episignature performance must also be carefully evaluated and addressed through imputation or an innovative approach to episignatures design. In this paper, we used data from patients suffering from Kabuki and Sotos syndrome to evaluate the influence of normalization methods, classification models and missing data on the prediction performances of two existing episignatures. We compare how six popular normalization methods for methylarray data affect episignature classification performances in Kabuki and Sotos syndromes and provide best practice suggestions when building new episignatures. In this setting, we show that Illumina, Noob or Funnorm normalization methods achieved higher classification performances on the testing sets compared to Quantile, Raw and Swan normalization methods. We further show that penalized logistic regression and support vector machines perform best in the classification of Kabuki and Sotos syndrome patients. Then, we describe a new paradigm to build episignatures based on the detection of differentially methylated regions (DMRs) and evaluate their performance compared to classical differentially methylated cytosines (DMCs)-based episignatures in the presence of missing data. We show that the performance of classical DMC-based episignatures suffers from the presence of missing data more than the DMR-based approach. We present a comprehensive evaluation of how the normalization of DNA methylation data affects episignature performance, using three popular classification models. We further evaluate how missing data affect those models' predictions. Finally, we propose a novel methodology to develop episignatures based on differentially methylated regions identification and show how this method slightly outperforms classical episignatures in the presence of missing data.

KAT6A mutations in Arboleda-Tham syndrome drive epigenetic regulation of posterior HOXC cluster

Arboleda-Tham Syndrome (ARTHS) is a rare genetic disorder caused by heterozygous, de novo mutations in Lysine(K) acetyltransferase 6A (KAT6A). ARTHS is clinically heterogeneous and characterized by several common features, including intellectual disability, developmental and speech delay, and hypotonia, and affects multiple organ systems. KAT6A is the enzymatic core of a histone-acetylation protein complex; however, the direct histone targets and gene regulatory effects remain unknown. In this study, we use ARTHS patient (n = 8) and control (n = 14) dermal fibroblasts and perform comprehensive profiling of the epigenome and transcriptome caused by KAT6A mutations. We identified differential chromatin accessibility within the promoter or gene body of 23% (14/60) of genes that were differentially expressed between ARTHS and controls. Within fibroblasts, we show a distinct set of genes from the posterior HOXC gene cluster (HOXC10, HOXC11, HOXC-AS3, HOXC-AS2, and HOTAIR) that are overexpressed in ARTHS and are transcription factors critical for early development body segment patterning. The genomic loci harboring HOXC genes are epigenetically regulated with increased chromatin accessibility, high levels of H3K23ac, and increased gene-body DNA methylation compared to controls, all of which are consistent with transcriptomic overexpression. Finally, we used unbiased proteomic mass spectrometry and identified two new histone post-translational modifications (PTMs) that are disrupted in ARTHS: H2A and H3K56 acetylation. Our multi-omics assays have identified novel histone and gene regulatory roles of KAT6A in a large group of ARTHS patients harboring diverse pathogenic mutations. This work provides insight into the role of KAT6A on the epigenomic regulation in somatic cell types.

Vitorino FN, Amano S, Lee J, Echeverria G, Gomez D, Garcia BA, Arboleda VA. KAT6A mutations in Arboleda-Tham syndrome drive epigenetic regulation of posterior HOXC cluster

Arboleda-Tham Syndrome (ARTHS) is a rare genetic disorder caused by heterozygous, de novo truncating mutations in Lysine(K) acetyltransferase 6A (KAT6A). ARTHS is clinically heterogeneous and characterized by several common features including intellectual disability, developmental and speech delay, hypotonia and affects multiple organ systems. KAT6A is highly expressed in early development and plays a key role in cell-type specific differentiation. KAT6A is the enzymatic core of a histone-acetylation protein complex, however the direct histone targets and gene regulatory effects remain unknown. In this study, we use ARTHS patient (n=8) and control (n=14) dermal fibroblasts and perform comprehensive profiling of the epigenome and transcriptome caused by KAT6A mutations. We identified differential chromatin accessibility within the promoter or gene body of 23%(14/60) of genes that were differentially expressed between ARTHS and controls. Within fibroblasts, we show a distinct set of genes from the posterior HOXC gene cluster (HOXC10, HOXC11, HOXC-AS3, HOXC-AS2, HOTAIR) that are overexpressed in ARTHS and are transcription factors critical for early development body segment patterning. The genomic loci harboring HOXC genes are epigenetically regulated with increased chromatin accessibility, high levels of H3K23ac, and increased gene-body DNA methylation compared to controls, all of which are consistent with transcriptomic overexpression. Finally, we used unbiased proteomic mass spectrometry and identified two new histone post-translational modifications (PTMs) that are disrupted in ARTHS: H2A and H3K56 acetylation. Our multi-omics assays have identified novel histone and gene regulatory roles of KAT6A in a large group of ARTHS patients harboring diverse pathogenic mutations. This work provides insight into the role of KAT6A on the epigenomic regulation in somatic cell types.

Assessing the Differential Methylation Analysis Quality for Microarray and NGS Platforms

Differential methylation (DM) is actively recruited in different types of fundamental and translational studies. Currently, microarray- and NGS-based approaches for methylation analysis are the most widely used with multiple statistical models designed to extract differential methylation signatures. The benchmarking of DM models is challenging due to the absence of gold standard data. In this study, we analyze an extensive number of publicly available NGS and microarray datasets with divergent and widely utilized statistical models and apply the recently suggested and validated rank-statistic-based approach Hobotnica to evaluate the quality of their results. Overall, microarray-based methods demonstrate more robust and convergent results, while NGS-based models are highly dissimilar. Tests on the simulated NGS data tend to overestimate the quality of the DM methods and therefore are recommended for use with caution. Evaluation of the top 10 DMC and top 100 DMC in addition to the not-subset signature also shows more stable results for microarray data. Summing up, given the observed heterogeneity in NGS methylation data, the evaluation of newly generated methylation signatures is a crucial step in DM analysis. The Hobotnica metric is coordinated with previously developed quality metrics and provides a robust, sensitive, and informative estimation of methods' performance and DM signatures' quality in the absence of gold standard data solving a long-existing problem in DM analysis.