Genomic DNA Methylation Signatures Enable Concurrent Diagnosis and Clinical Genetic Variant Classification in Neurodevelopmental Syndromes

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

Identification of a DNA methylation episignature for recurrent constellations of embryonic malformations

The term "recurrent constellations of embryonic malformations" (RCEM) is used to describe a number of multiple malformation associations that affect three or more body structures. The causes of these disorders are currently unknown, and no diagnostic marker has been identified. Consequently, providing a definitive diagnosis in suspected individuals is challenging. In this study, genome-wide DNA methylation analysis was conducted on DNA samples obtained from the peripheral blood of 53 individuals with RCEM characterized by clinical features recognized as VACTERL and/or oculoauriculovertebral spectrum association. We identified a common DNA methylation episignature in 40 out of the 53 individuals. Subsequently, a sensitive and specific binary classifier was developed based on the DNA methylation episignature. This classifier can facilitate the use of RCEM episignature as a diagnostic biomarker in a clinical setting. The study also investigated the functional correlation of RCEM DNA methylation relative to other genetic disorders with known episignatures, highlighting the common genomic regulatory pathways involved in the pathophysiology of RCEM.

Loss-of-function of the Zinc Finger Homeobox 4 (ZFHX4) gene underlies a neurodevelopmental disorder

8q21.11 microdeletions encompassing the gene encoding transcription factor ZFHX4, have previously been associated by us with a syndromic form of intellectual disability, hypotonia, decreased balance and hearing loss. Here, we report on 57 individuals, 52 probands and 5 affected family members, with protein truncating variants (n=36), (micro)deletions (n=20) or an inversion (n=1) affecting ZFHX4 with variable developmental delay and intellectual disability, distinctive facial characteristics, morphological abnormalities of the central nervous system, behavioral alterations, short stature, hypotonia, and occasionally cleft palate and anterior segment dysgenesis. The phenotypes associated with 8q21.11 microdeletions and ZFHX4 intragenic loss-of-function variants largely overlap, identifying ZFHX4 as the main driver for the microdeletion syndrome, although leukocyte-derived DNA shows a mild common methylation profile for (micro)deletions only. We identify ZFHX4 as a transcription factor that is increasingly expressed during human brain development and neuronal differentiation. Furthermore, ZFHX4 interacting factors identified via IP-MS in neural progenitor cells, suggest an important role for ZFHX4 in cellular and developmental pathways, especially during histone modifications, cytosolic transport and development. Additionally, using CUT&RUN, we observed that ZFHX4 binds with the promoter regions of genes with crucial roles in embryonic, neuron and axon development. Since loss-of-function variants in ZFHX4 are found with consistent dysmorphic facial features, we investigated whether the disruption of zfhx4 causes craniofacial abnormalities in zebrafish. First-generation (F0) zfhx4 crispant zebrafish, (mosaic) mutant for zfhx4 loss-of-function variants, have significantly shorter Meckel's cartilages and smaller ethmoid plates compared to control zebrafish. Furthermore, behavioral assays show a decreased movement frequency in the zfhx4 crispant zebrafish in comparison with control zebrafish larvae. Although further research is needed, our in vivo work suggests a role for zfhx4 in facial skeleton patterning, palatal development and behavior.

Diagnostic utility of DNA methylation analysis in genetically unsolved pediatric epilepsies and CHD2 episignature refinement

Sequence-based genetic testing identifies causative variants in ~ 50% of individuals with developmental and epileptic encephalopathies (DEEs). Aberrant changes in DNA methylation are implicated in various neurodevelopmental disorders but remain unstudied in DEEs. We interrogate the diagnostic utility of genome-wide DNA methylation array analysis on peripheral blood samples from 582 individuals with genetically unsolved DEEs. We identify rare differentially methylated regions (DMRs) and explanatory episignatures to uncover causative and candidate genetic etiologies in 12 individuals. Using long-read sequencing, we identify DNA variants underlying rare DMRs, including one balanced translocation, three CG-rich repeat expansions, and four copy number variants. We also identify pathogenic variants associated with episignatures. Finally, we refine the CHD2 episignature using an 850 K methylation array and bisulfite sequencing to investigate potential insights into CHD2 pathophysiology. Our study demonstrates the diagnostic yield of genome-wide DNA methylation analysis to identify causal and candidate variants as 2% (12/582) for unsolved DEE cases.

Next-generation phenotyping integrated in a national framework for patients with ultrarare disorders improves genetic diagnostics and yields new molecular findings

Individuals with ultrarare disorders pose a structural challenge for healthcare systems since expert clinical knowledge is required to establish diagnoses. In TRANSLATE NAMSE, a 3-year prospective study, we evaluated a novel diagnostic concept based on multidisciplinary expertise in Germany. Here we present the systematic investigation of the phenotypic and molecular genetic data of 1,577 patients who had undergone exome sequencing and were partially analyzed with next-generation phenotyping approaches. Molecular genetic diagnoses were established in 32% of the patients totaling 370 distinct molecular genetic causes, most with prevalence below 1:50,000. During the diagnostic process, 34 novel and 23 candidate genotype-phenotype associations were identified, mainly in individuals with neurodevelopmental disorders. Sequencing data of the subcohort that consented to computer-assisted analysis of their facial images with GestaltMatcher could be prioritized more efficiently compared with approaches based solely on clinical features and molecular scores. Our study demonstrates the synergy of using next-generation sequencing and phenotyping for diagnosing ultrarare diseases in routine healthcare and discovering novel etiologies by multidisciplinary teams.

The missing link: ARID1B non-truncating variants causing Coffin-Siris syndrome due to protein aggregation

ARID1B is the most frequently mutated gene in Coffin-Siris syndrome (CSS). To date, the vast majority of causative variants reported in ARID1B are truncating, leading to nonsense-mediated mRNA decay. In the absence of experimental data, only few ARID1B amino acid substitutions have been classified as pathogenic, mainly based on clinical data and their de novo occurrence, while most others are currently interpreted as variants of unknown significance. The present study substantiates the pathogenesis of ARID1B non-truncating/NMD-escaping variants located in the SMARCA4-interacting EHD2 and DNA-binding ARID domains. Overexpression assays in cell lines revealed that the majority of EHD2 variants lead to protein misfolding and formation of cytoplasmic aggresomes surrounded by vimentin cage-like structures and co-localizing with the microtubule organisation center. ARID domain variants exhibited not only aggresomes, but also nuclear aggregates, demonstrating robust pathological effects. Protein levels were not compromised, as shown by quantitative western blot analysis. In silico structural analysis predicted the exposure of amylogenic segments in both domains due to the nearby variants, likely causing this aggregation. Genome-wide transcriptome and methylation analysis in affected individuals revealed expression and methylome patterns consistent with those of the pathogenic haploinsufficiency ARID1B alterations in CSS cases. These results further support pathogenicity and indicate two approaches for disambiguation of such variants in everyday practice. The few affected individuals harbouring EHD2 non-truncating variants described to date exhibit mild CSS clinical traits. In summary, this study paves the way for the re-evaluation of previously unclear ARID1B non-truncating variants and opens a new era in CSS genetic diagnosis.

Discovery of DNA methylation signature in the peripheral blood of individuals with history of antenatal exposure to valproic acid

PURPOSE

Valproic acid or valproate is an effective antiepileptic drug; however, embryonic exposure to valproate can result in a teratogenic disorder referred to as fetal valproate syndrome (OMIM #609442). Currently there are no diagnostic biomarkers for the condition. This study aims to define an episignature biomarker for teratogenic antenatal exposure to valproate.

METHODS

DNA extracted from peripheral blood of individuals with teratogenic antenatal exposure to valproate was processed using DNA methylation microarrays. Subsequently, methylation profiling and construction of support vector machine classifiers were performed in R.

RESULTS

Genomic DNA methylation analysis was applied, and a distinct DNA methylation profile was identified in the majority of affected individuals. This profile was used to develop a diagnostic episignature classifier. The valproate exposure episignature exhibited high sensitivity and specificity relative to a large reference data set of unaffected controls and individuals with a wide spectrum of syndromic disorders with episignatures. Gene set enrichment analysis demonstrated an enrichment for terms associated with cell adhesion, including significant overrepresentation of the cadherin superfamily.

CONCLUSION

This study provides evidence of a robust peripheral blood-based diagnostic epigenetic biomarker for a prenatal teratogenic disorder.

DNA methylation episignature and comparative epigenomic profiling for Pitt-Hopkins syndrome caused by TCF4 variants

Pitt-Hopkins syndrome (PTHS) is a neurodevelopmental disorder caused by pathogenic variants in TCF4, leading to intellectual disability, specific morphological features, and autonomic nervous system dysfunction. Epigenetic dysregulation has been implicated in PTHS, prompting the investigation of a DNA methylation (DNAm) "episignature" specific to PTHS for diagnostic purposes and variant reclassification and functional insights into the molecular pathophysiology of this disorder. A cohort of 67 individuals with genetically confirmed PTHS and three individuals with intellectual disability and a variant of uncertain significance (VUS) in TCF4 were studied. The DNAm episignature was developed with an Infinium Methylation EPIC BeadChip array analysis using peripheral blood cells. Support vector machine (SVM) modeling and clustering methods were employed to generate a DNAm classifier for PTHS. Validation was extended to an additional cohort of 11 individuals with PTHS. The episignature was assessed in relation to other neurodevelopmental disorders and its specificity was examined. A specific DNAm episignature for PTHS was established. The classifier exhibited high sensitivity for TCF4 haploinsufficiency and missense variants in the basic-helix-loop-helix domain. Notably, seven individuals with TCF4 variants exhibited negative episignatures, suggesting complexities related to mosaicism, genetic factors, and environmental influences. The episignature displayed degrees of overlap with other related disorders and biological pathways. This study defines a DNAm episignature for TCF4-related PTHS, enabling improved diagnostic accuracy and VUS reclassification. The finding that some cases scored negatively underscores the potential for multiple or nested episignatures and emphasizes the need for continued investigation to enhance specificity and coverage across PTHS-related variants.

MSL2 variants lead to a neurodevelopmental syndrome with lack of coordination, epilepsy, specific dysmorphisms, and a distinct episignature

Epigenetic dysregulation has emerged as an important etiological mechanism of neurodevelopmental disorders (NDDs). Pathogenic variation in epigenetic regulators can impair deposition of histone post-translational modifications leading to aberrant spatiotemporal gene expression during neurodevelopment. The male-specific lethal (MSL) complex is a prominent multi-subunit epigenetic regulator of gene expression and is responsible for histone 4 lysine 16 acetylation (H4K16ac). Using exome sequencing, here we identify a cohort of 25 individuals with heterozygous de novo variants in MSL complex member MSL2. MSL2 variants were associated with NDD phenotypes including global developmental delay, intellectual disability, hypotonia, and motor issues such as coordination problems, feeding difficulties, and gait disturbance. Dysmorphisms and behavioral and/or psychiatric conditions, including autism spectrum disorder, and to a lesser extent, seizures, connective tissue disease signs, sleep disturbance, vision problems, and other organ anomalies, were observed in affected individuals. As a molecular biomarker, a sensitive and specific DNA methylation episignature has been established. Induced pluripotent stem cells (iPSCs) derived from three members of our cohort exhibited reduced MSL2 levels. Remarkably, while NDD-associated variants in two other members of the MSL complex (MOF and MSL3) result in reduced H4K16ac, global H4K16ac levels are unchanged in iPSCs with MSL2 variants. Regardless, MSL2 variants altered the expression of MSL2 targets in iPSCs and upon their differentiation to early germ layers. Our study defines an MSL2-related disorder as an NDD with distinguishable clinical features, a specific blood DNA episignature, and a distinct, MSL2-specific molecular etiology compared to other MSL complex-related disorders.

Multi-omics based artificial intelligence for cancer research

With significant advancements of next generation sequencing technologies, large amounts of multi-omics data, including genomics, epigenomics, transcriptomics, proteomics, and metabolomics, have been accumulated, offering an unprecedented opportunity to explore the heterogeneity and complexity of cancer across various molecular levels and scales. One of the promising aspects of multi-omics lies in its capacity to offer a holistic view of the biological networks and pathways underpinning cancer, facilitating a deeper understanding of its development, progression, and response to treatment. However, the exponential growth of data generated by multi-omics studies present significant analytical challenges. Processing, analyzing, integrating, and interpreting these multi-omics datasets to extract meaningful insights is an ambitious task that stands at the forefront of current cancer research. The application of artificial intelligence (AI) has emerged as a powerful solution to these challenges, demonstrating exceptional capabilities in deciphering complex patterns and extracting valuable information from large-scale, intricate omics datasets. This review delves into the synergy of AI and multi-omics, highlighting its revolutionary impact on oncology. We dissect how this confluence is reshaping the landscape of cancer research and clinical practice, particularly in the realms of early detection, diagnosis, prognosis, treatment and pathology. Additionally, we elaborate the latest AI methods for multi-omics integration to provide a comprehensive insight of the complex biological mechanisms and inherent heterogeneity of cancer. Finally, we discuss the current challenges of data harmonization, algorithm interpretability, and ethical considerations. Addressing these challenges necessitates a multidisciplinary collaboration, paving the promising way for more precise, personalized, and effective treatments for cancer patients.

The detection of a strong episignature for Chung-Jansen syndrome, partially overlapping with Börjeson-Forssman-Lehmann and White-Kernohan syndromes

Chung-Jansen syndrome is a neurodevelopmental disorder characterized by intellectual disability, behavioral problems, obesity and dysmorphic features. It is caused by pathogenic variants in the PHIP gene that encodes for the Pleckstrin homology domain-interacting protein, which is part of an epigenetic modifier protein complex. Therefore, we hypothesized that PHIP haploinsufficiency may impact genome-wide DNA methylation (DNAm). We assessed the DNAm profiles of affected individuals with pathogenic and likely pathogenic PHIP variants with Infinium Methylation EPIC arrays and report a specific and sensitive DNAm episignature biomarker for Chung-Jansen syndrome. In addition, we observed similarities between the methylation profile of Chung-Jansen syndrome and that of functionally related and clinically partially overlapping genetic disorders, White-Kernohan syndrome (caused by variants in DDB1 gene) and Börjeson-Forssman-Lehmann syndrome (caused by variants in PHF6 gene). Based on these observations we also proceeded to develop a common episignature biomarker for these disorders. These newly defined episignatures can be used as part of a multiclass episignature classifier for screening of affected individuals with rare disorders and interpretation of genetic variants of unknown clinical significance, and provide further insights into the common molecular pathophysiology of the clinically-related Chung-Jansen, Börjeson-Forssman-Lehmann and White-Kernohan syndromes.

Diagnostic utility and reporting recommendations for clinical DNA methylation episignature testing in genetically undiagnosed rare diseases

PURPOSE

This study aims to assess the diagnostic utility and provide reporting recommendations for clinical DNA methylation episignature testing based on the cohort of patients tested through the EpiSign Clinical Testing Network.

METHODS

The EpiSign assay utilized unsupervised clustering techniques and a support vector machine-based classification algorithm to compare each patient's genome-wide DNA methylation profile with the EpiSign Knowledge Database, yielding the result that was reported. An international working group, representing distinct EpiSign Clinical Testing Network health jurisdictions, collaborated to establish recommendations for interpretation and reporting of episignature testing.

RESULTS

Among 2399 cases analyzed, 1667 cases underwent a comprehensive screen of validated episignatures, imprinting, and promoter regions, resulting in 18.7% (312/1667) positive reports. The remaining 732 referrals underwent targeted episignature analysis for assessment of sequence or copy-number variants (CNVs) of uncertain significance or for assessment of clinical diagnoses without confirmed molecular findings, and 32.4% (237/732) were positive. Cases with detailed clinical information were highlighted to describe various utility scenarios for episignature testing.

CONCLUSION

Clinical DNA methylation testing including episignatures, imprinting, and promoter analysis provided by an integrated network of clinical laboratories enables test standardization and demonstrates significant diagnostic yield and clinical utility beyond DNA sequence analysis in rare diseases.

DNA methylation signatures for chromatinopathies: current challenges and future applications

Pathogenic variants in genes that encode epigenetic regulators are the cause for more than 100 rare neurodevelopmental syndromes also termed "chromatinopathies". DNA methylation signatures, syndrome-specific patterns of DNA methylation alterations, serve as both a research avenue for elucidating disease pathophysiology and a clinical diagnostic tool. The latter is well established, especially for the classification of variants of uncertain significance (VUS). In this perspective, we describe the seminal DNA methylation signature research in chromatinopathies; the complex relationships between genotype, phenotype and DNA methylation, and the future applications of DNA methylation signatures.

Computational deconvolution of DNA methylation data from mixed DNA samples

In this review, we provide a comprehensive overview of the different computational tools that have been published for the deconvolution of bulk DNA methylation (DNAm) data. Here, deconvolution refers to the estimation of cell-type proportions that constitute a mixed sample. The paper reviews and compares 25 deconvolution methods (supervised, unsupervised or hybrid) developed between 2012 and 2023 and compares the strengths and limitations of each approach. Moreover, in this study, we describe the impact of the platform used for the generation of methylation data (including microarrays and sequencing), the applied data pre-processing steps and the used reference dataset on the deconvolution performance. Next to reference-based methods, we also examine methods that require only partial reference datasets or require no reference set at all. In this review, we provide guidelines for the use of specific methods dependent on the DNA methylation data type and data availability.

Identification of DNA methylation episignature for the intellectual developmental disorder, autosomal dominant 21 syndrome, caused by variants in the CTCF gene

PURPOSE

The main objective of this study was to assess clinical features and genome-wide DNA methylation profiles in individuals affected by intellectual developmental disorder, autosomal dominant 21 (IDD21) syndrome, caused by variants in the CCCTC-binding factor (CTCF) gene.

METHODS

DNA samples were extracted from peripheral blood of 16 individuals with clinical features and genetic findings consistent with IDD21. DNA methylation analysis was performed using the Illumina Infinium Methylation EPIC Bead Chip microarrays. The methylation levels were fitted in a multivariate linear regression model to identify the differentially methylated probes. A binary support vector machine classification model was constructed to differentiate IDD21 samples from controls.

RESULTS

We identified a highly specific, reproducible, and sensitive episignature associated with CTCF variants. Six variants of uncertain significance were tested, of which 2 mapped to the IDD21 episignature and clustered alongside IDD21 cases in both heatmap and multidimensional scaling plots. Comparison of the genomic DNA methylation profile of IDD21 with that of 56 other neurodevelopmental disorders provided insights into the underlying molecular pathophysiology of this disorder.

CONCLUSION

The robust and specific CTCF/IDD21 episignature expands the growing list of neurodevelopmental disorders with distinct DNA methylation profiles, which can be applied as supporting evidence in variant classification.

DNA methylation episignature, extension of the clinical features, and comparative epigenomic profiling of Hao-Fountain syndrome caused by variants in USP7

PURPOSE

Hao-Fountain syndrome (HAFOUS) is a neurodevelopmental disorder caused by pathogenic variants in USP7. HAFOUS is characterized by developmental delay, intellectual disability, speech delay, behavioral abnormalities, autism spectrum disorder, seizures, hypogonadism, and mild dysmorphic features. We investigated the phenotype of 18 participants with HAFOUS and performed DNA methylation (DNAm) analysis, aiming to generate a diagnostic biomarker. Furthermore, we performed comparative analysis with known episignatures to gain more insight into the molecular pathophysiology of HAFOUS.

METHODS

We assessed genomic DNAm profiles of 18 individuals with pathogenic variants and variants of uncertain significance (VUS) in USP7 to map and validate a specific episignature. The comparison between the USP7 cohort and 56 rare genetic disorders with earlier reported DNAm episignatures was performed with statistical and functional correlation.

RESULTS

We mapped a sensitive and specific DNAm episignature for pathogenic variants in USP7 and utilized this to reclassify the VUS. Comparative epigenomic analysis showed evidence of HAFOUS similarity to a number of other rare genetic episignature disorders.

CONCLUSION

We discovered a sensitive and specific DNAm episignature as a robust diagnostic biomarker for HAFOUS that enables VUS reclassification in USP7. We also expand the phenotypic spectrum of 9 new and 5 previously reported individuals with HAFOUS.

Applications of Machine Learning (ML) and Mathematical Modeling (MM) in Healthcare with Special Focus on Cancer Prognosis and Anticancer Therapy: Current Status and Challenges

The use of data-driven high-throughput analytical techniques, which has given rise to computational oncology, is undisputed. The widespread use of machine learning (ML) and mathematical modeling (MM)-based techniques is widely acknowledged. These two approaches have fueled the advancement in cancer research and eventually led to the uptake of telemedicine in cancer care. For diagnostic, prognostic, and treatment purposes concerning different types of cancer research, vast databases of varied information with manifold dimensions are required, and indeed, all this information can only be managed by an automated system developed utilizing ML and MM. In addition, MM is being used to probe the relationship between the pharmacokinetics and pharmacodynamics (PK/PD interactions) of anti-cancer substances to improve cancer treatment, and also to refine the quality of existing treatment models by being incorporated at all steps of research and development related to cancer and in routine patient care. This review will serve as a consolidation of the advancement and benefits of ML and MM techniques with a special focus on the area of cancer prognosis and anticancer therapy, leading to the identification of challenges (data quantity, ethical consideration, and data privacy) which are yet to be fully addressed in current studies.

A mouse model of Weaver syndrome displays overgrowth and excess osteogenesis reversible with KDM6A/6B inhibition

Weaver syndrome is a Mendelian disorder of the epigenetic machinery (MDEM) caused by germline pathogenic variants in EZH2, which encodes the predominant H3K27 methyltransferase and key enzymatic component of Polycomb repressive complex 2 (PRC2). Weaver syndrome is characterized by striking overgrowth and advanced bone age, intellectual disability, and distinctive facies. We generated a mouse model for the most common Weaver syndrome missense variant, EZH2 p.R684C. Ezh2R684C/R684C mouse embryonic fibroblasts (MEFs) showed global depletion of H3K27me3. Ezh2R684C/+ mice had abnormal bone parameters, indicative of skeletal overgrowth, and Ezh2R684C/+ osteoblasts showed increased osteogenic activity. RNA-Seq comparing osteoblasts differentiated from Ezh2R684C/+, and Ezh2+/+ BM-mesenchymal stem cells (BM-MSCs) indicated collective dysregulation of the BMP pathway and osteoblast differentiation. Inhibition of the opposing H3K27 demethylases KDM6A and KDM6B substantially reversed the excessive osteogenesis in Ezh2R684C/+ cells both at the transcriptional and phenotypic levels. This supports both the ideas that writers and erasers of histone marks exist in a fine balance to maintain epigenome state and that epigenetic modulating agents have therapeutic potential for the treatment of MDEMs.

Clinico-biological refinement of BCL11B-related disorder and identification of an episignature: A series of 20 unreported individuals

PURPOSE

BCL11B-related disorder (BCL11B-RD) arises from rare genetic variants within the BCL11B gene, resulting in a distinctive clinical spectrum encompassing syndromic neurodevelopmental disorder, with or without intellectual disability, associated with facial features and impaired immune function. This study presents an in-depth clinico-biological analysis of 20 newly reported individuals with BCL11B-RD, coupled with a characterization of genome-wide DNA methylation patterns of this genetic condition.

METHODS

Through an international collaboration, clinical and molecular data from 20 individuals were systematically gathered, and a comparative analysis was conducted between this series and existing literature. We further scrutinized peripheral blood DNA methylation profile of individuals with BCL11B-RD, contrasting them with healthy controls and other neurodevelopmental disorders marked by established episignature.

RESULTS

Our findings unveil rarely documented clinical manifestations, notably including Rubinstein-Taybi-like facial features, craniosynostosis, and autoimmune disorders, all manifesting within the realm of BCL11B-RD. We refine the intricacies of T cell compartment alterations of BCL11B-RD, revealing decreased levels naive CD4+ T cells and recent thymic emigrants while concurrently observing an elevated proportion of effector-memory expressing CD45RA CD8+ T cells (TEMRA). Finally, a distinct DNA methylation episignature exclusive to BCL11B-RD is unveiled.

CONCLUSION

This study serves to enrich our comprehension of the clinico-biological landscape of BCL11B-RD, potentially furnishing a more precise framework for diagnosis and follow-up of individuals carrying pathogenic BCL11B variant. Moreover, the identification of a unique DNA methylation episignature offers a valuable diagnosis tool for BCL11B-RD, thereby facilitating routine clinical practice by empowering physicians to reevaluate variants of uncertain significance within the BCL11B gene.

Solving the unsolved genetic epilepsies: Current and future perspectives

Many patients with epilepsy undergo exome or genome sequencing as part of a diagnostic workup; however, many remain genetically unsolved. There are various factors that account for negative results in exome/genome sequencing for patients with epilepsy: (1) the underlying cause is not genetic; (2) there is a complex polygenic explanation; (3) the illness is monogenic but the causative gene remains to be linked to a human disorder; (4) family segregation with reduced penetrance; (5) somatic mosaicism or the complexity of, for example, a structural rearrangement; or (6) limited knowledge or diagnostic tools that hinder the proper classification of a variant, resulting in its designation as a variant of unknown significance. The objective of this review is to outline some of the diagnostic options that lie beyond the exome/genome, and that might become clinically relevant within the foreseeable future. These options include: (1) re-analysis of older exome/genome data as knowledge increases or symptoms change; (2) looking for somatic mosaicism or long-read sequencing to detect low-complexity repeat variants or specific structural variants missed by traditional exome/genome sequencing; (3) exploration of the non-coding genome including disruption of topologically associated domains, long range non-coding RNA, or other regulatory elements; and finally (4) transcriptomics, DNA methylation signatures, and metabolomics as complementary diagnostic methods that may be used in the assessment of variants of unknown significance. Some of these tools are currently not integrated into standard diagnostic workup. However, it is reasonable to expect that they will become increasingly available and improve current diagnostic capabilities, thereby enabling precision diagnosis in patients who are currently undiagnosed.

Epigenotype-genotype-phenotype correlations in SETD1A and SETD2 chromatin disorders

Germline pathogenic variants in two genes encoding the lysine-specific histone methyltransferase genes SETD1A and SETD2 are associated with neurodevelopmental disorders (NDDs) characterized by developmental delay and congenital anomalies. The SETD1A and SETD2 gene products play a critical role in chromatin-mediated regulation of gene expression. Specific methylation episignatures have been detected for a range of chromatin gene-related NDDs and have impacted clinical practice by improving the interpretation of variant pathogenicity. To investigate if SETD1A and/or SETD2-related NDDs are associated with a detectable episignature, we undertook targeted genome-wide methylation profiling of > 2 M CpGs using a next-generation sequencing-based assay. A comparison of methylation profiles in patients with SETD1A variants (n = 6) did not reveal evidence of a strong methylation episignature. A review of the clinical and genetic features of the SETD2 patient group revealed that, as reported previously, there were phenotypic differences between patients with truncating mutations (n = 4, Luscan-Lumish syndrome; MIM:616831) and those with missense codon 1740 variants [p.Arg1740Trp (n = 4) and p.Arg1740Gln (n = 2)]. Both SETD2 subgroups demonstrated a methylation episignature, which was characterized by hypomethylation and hypermethylation events, respectively. Within the codon 1740 subgroup, both the methylation changes and clinical phenotype were more severe in those with p.Arg1740Trp variants. We also noted that two of 10 cases with a SETD2-NDD had developed a neoplasm. These findings reveal novel epigenotype-genotype-phenotype correlations in SETD2-NDDs and predict a gain-of-function mechanism for SETD2 codon 1740 pathogenic variants.

Peripheral blood DNA methylation and neuroanatomical responses to HDACi treatment that rescues neurological deficits in a Kabuki syndrome mouse model

BACKGROUND

Recent findings from studies of mouse models of Mendelian disorders of epigenetic machinery strongly support the potential for postnatal therapies to improve neurobehavioral and cognitive deficits. As several of these therapies move into human clinical trials, the search for biomarkers of treatment efficacy is a priority. A potential postnatal treatment of Kabuki syndrome type 1 (KS1), caused by pathogenic variants in KMT2D encoding a histone-lysine methyltransferase, has emerged using a mouse model of KS1 (Kmt2d+/βGeo). In this mouse model, hippocampal memory deficits are ameliorated following treatment with the histone deacetylase inhibitor (HDACi), AR-42. Here, we investigate the effect of both Kmt2d+/βGeo genotype and AR-42 treatment on neuroanatomy and on DNA methylation (DNAm) in peripheral blood. While peripheral blood may not be considered a "primary tissue" with respect to understanding the pathophysiology of neurodevelopmental disorders, it has the potential to serve as an accessible biomarker of disease- and treatment-related changes in the brain.

METHODS

Half of the KS1 and wildtype mice were treated with 14 days of AR-42. Following treatment, fixed brain samples were imaged using MRI to calculate regional volumes. Blood was assayed for genome-wide DNAm at over 285,000 CpG sites using the Illumina Infinium Mouse Methylation array. DNAm patterns and brain volumes were analyzed in the four groups of animals: wildtype untreated, wildtype AR-42 treated, KS1 untreated and KS1 AR-42 treated.

RESULTS

We defined a DNAm signature in the blood of KS1 mice, that overlapped with the human KS1 DNAm signature. We also found a striking 10% decrease in total brain volume in untreated KS1 mice compared to untreated wildtype, which correlated with DNAm levels in a subset KS1 signature sites, suggesting that disease severity may be reflected in blood DNAm. Treatment with AR-42 ameliorated DNAm aberrations in KS1 mice at a small number of signature sites.

CONCLUSIONS

As this treatment impacts both neurological deficits and blood DNAm in mice, future KS clinical trials in humans could be used to assess blood DNAm as an early biomarker of therapeutic efficacy.

Diagnostic Utility of Genome-wide DNA Methylation Analysis in Genetically Unsolved Developmental and Epileptic Encephalopathies and Refinement of a CHD2 Episignature

Sequence-based genetic testing currently identifies causative genetic variants in ∼50% of individuals with developmental and epileptic encephalopathies (DEEs). Aberrant changes in DNA methylation are implicated in various neurodevelopmental disorders but remain unstudied in DEEs. Rare epigenetic variations ("epivariants") can drive disease by modulating gene expression at single loci, whereas genome-wide DNA methylation changes can result in distinct "episignature" biomarkers for monogenic disorders in a growing number of rare diseases. Here, we interrogate the diagnostic utility of genome-wide DNA methylation array analysis on peripheral blood samples from 516 individuals with genetically unsolved DEEs who had previously undergone extensive genetic testing. We identified rare differentially methylated regions (DMRs) and explanatory episignatures to discover causative and candidate genetic etiologies in 10 individuals. We then used long-read sequencing to identify DNA variants underlying rare DMRs, including one balanced translocation, three CG-rich repeat expansions, and two copy number variants. We also identify pathogenic sequence variants associated with episignatures; some had been missed by previous exome sequencing. Although most DEE genes lack known episignatures, the increase in diagnostic yield for DNA methylation analysis in DEEs is comparable to the added yield of genome sequencing. Finally, we refine an episignature for CHD2 using an 850K methylation array which was further refined at higher CpG resolution using bisulfite sequencing to investigate potential insights into CHD2 pathophysiology. Our study demonstrates the diagnostic yield of genome-wide DNA methylation analysis to identify causal and candidate genetic causes as ∼2% (10/516) for unsolved DEE cases.

Functional Insight into and Refinement of the Genomic Boundaries of the JARID2-Neurodevelopmental Disorder Episignature

JARID2 (Jumonji, AT-rich interactive domain 2) haploinsufficiency is associated with a clinically distinct neurodevelopmental syndrome. It is characterized by intellectual disability, developmental delay, autistic features, behavior abnormalities, cognitive impairment, hypotonia, and dysmorphic features. JARID2 acts as a transcriptional repressor protein that is involved in the regulation of histone methyltransferase complexes. JARID2 plays a role in the epigenetic machinery, and the associated syndrome has an identified DNA methylation episignature derived from sequence variants and intragenic deletions involving JARID2. For this study, our aim was to determine whether patients with larger deletions spanning beyond JARID2 present a similar DNA methylation episignature and to define the critical region involved in aberrant DNA methylation in 6p22-p24 microdeletions. We examined the DNA methylation profiles of peripheral blood from 56 control subjects, 13 patients with (likely) pathogenic JARID2 variants or patients carrying copy number variants, and three patients with JARID2 VUS variants. The analysis showed a distinct and strong differentiation between patients with (likely) pathogenic variants, both sequence and copy number, and controls. Using the identified episignature, we developed a binary model to classify patients with the JARID2-neurodevelopmental syndrome. DNA methylation analysis indicated that JARID2 is the driver gene for aberrant DNA methylation observed in 6p22-p24 microdeletions. In addition, we performed analysis of functional correlation of the JARID2 genome-wide methylation profile with the DNA methylation profiles of 56 additional neurodevelopmental disorders. To conclude, we refined the critical region for the presence of the JARID2 episignature in 6p22-p24 microdeletions and provide insight into the functional changes in the epigenome observed when regulation by JARID2 is lost.

Germline pathogenic variants in HNRNPU are associated with alterations in blood methylome

HNRNPU encodes a multifunctional RNA-binding protein that plays critical roles in regulating pre-mRNA splicing, mRNA stability, and translation. Aberrant expression and dysregulation of HNRNPU have been implicated in various human diseases, including cancers and neurological disorders. We applied a next generation sequencing based assay (EPIC-NGS) to investigate genome-wide methylation profiling for >2 M CpGs for 7 individuals with a neurodevelopmental disorder associated with HNRNPU germline pathogenic loss-of-function variants. Compared to healthy individuals, 227 HNRNPU-associated differentially methylated positions were detected. Both hyper- and hypomethylation alterations were identified but the former predominated. The identification of a methylation episignature for HNRNPU-associated neurodevelopmental disorder (NDD) implicates HNPRNPU-related chromatin alterations in the aetiopathogenesis of this disorder and suggests that episignature profiling should have clinical utility as a predictor for the pathogenicity of HNRNPU variants of uncertain significance. The detection of a methylation episignaure for HNRNPU-associated NDD is consistent with a recent report of a methylation episignature for HNRNPK-associated NDD.

Redefining cerebral palsies as a diverse group of neurodevelopmental disorders with genetic aetiology

Cerebral palsy is a clinical descriptor covering a diverse group of permanent, non-degenerative disorders of motor function. Around one-third of cases have now been shown to have an underlying genetic aetiology, with the genetic landscape overlapping with those of neurodevelopmental disorders including intellectual disability, epilepsy, speech and language disorders and autism. Here we review the current state of genomic testing in cerebral palsy, highlighting the benefits for personalized medicine and the imperative to consider aetiology during clinical diagnosis. With earlier clinical diagnosis now possible, we emphasize the opportunity for comprehensive and early genomic testing as a crucial component of the routine diagnostic work-up in people with cerebral palsy.

Episignature analysis of moderate effects and mosaics

DNA methylation classifiers ("episignatures") help to determine the pathogenicity of variants of uncertain significance (VUS). However, their sensitivity is limited due to their training on unambiguous cases with strong-effect variants so that the classification of variants with reduced effect size or in mosaic state may fail. Moreover, episignature evaluation of mosaics as a function of their degree of mosaicism has not been developed so far. We improved episignatures with respect to three categories. Applying (i) minimum-redundancy-maximum-relevance feature selection we reduced their length by up to one order of magnitude without loss of accuracy. Performing (ii) repeated re-training of a support vector machine classifier by step-wise inclusion of cases in the training set that reached probability scores larger than 0.5, we increased the sensitivity of the episignature-classifiers by 30%. In the newly diagnosed patients we confirmed the association between DNA methylation aberration and age at onset of KMT2B-deficient dystonia. Moreover, we found evidence for allelic series, including KMT2B-variants with moderate effects and comparatively mild phenotypes such as late-onset focal dystonia. Retrained classifiers also can detect mosaics that previously remained below the 0.5-threshold, as we showed for KMT2D-associated Kabuki syndrome. Conversely, episignature-classifiers are able to revoke erroneous exome calls of mosaicism, as we demonstrated by (iii) comparing presumed mosaic cases with a distribution of artificial in silico-mosaics that represented all the possible variation in degree of mosaicism, variant read sampling and methylation analysis.

Beyond the exome: What's next in diagnostic testing for Mendelian conditions

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order, and emerging technologies, such as optical genome mapping and long-read DNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to research consortia focused on elucidating the underlying cause of rare unsolved genetic disorders.

Identification of a DNA methylation signature for Renpenning syndrome (RENS1), a spliceopathy

The challenges and ambiguities in providing an accurate diagnosis for patients with neurodevelopmental disorders have led researchers to apply epigenetics as a technique to validate the diagnosis provided based on the clinical examination and genetic testing results. Genome-wide DNA methylation analysis has recently been adapted for clinical testing of patients with genetic neurodevelopmental disorders. In this paper, preliminary data demonstrating a DNA methylation signature for Renpenning syndrome (RENS1 - OMIM 309500), which is an X-linked recessive neurodevelopmental disorder caused by variants in polyglutamine-binding protein 1 (PQBP1) is reported. The identified episignature was then utilized to construct a highly sensitive and specific binary classification model. Besides providing evidence for the existence of a DNA methylation episignature for Renpenning syndrome, this study increases the knowledge of the molecular mechanisms related to the disease. Moreover, the availability of more subjects in future may facilitate the establishment of an episignature that can be utilized for diagnosis in a clinical setting and for reclassification of variants of unknown clinical significance.

Novel mouse model of Weaver syndrome displays overgrowth and excess osteogenesis reversible with KDM6A/6B inhibition

Weaver syndrome is a Mendelian disorder of the epigenetic machinery (MDEM) caused by germline pathogenic variants in EZH2, which encodes the predominant H3K27 methyltransferase and key enzymatic component of Polycomb repressive complex 2 (PRC2). Weaver syndrome is characterized by striking overgrowth and advanced bone age, intellectual disability, and distinctive facies. We generated a mouse model for the most common Weaver syndrome missense variant, EZH2 p.R684C. Ezh2R684C/R684C mouse embryonic fibroblasts (MEFs) showed global depletion of H3K27me3. Ezh2R684C/+ mice had abnormal bone parameters indicative of skeletal overgrowth, and Ezh2R684C/+ osteoblasts showed increased osteogenic activity. RNA-seq comparing osteoblasts differentiated from Ezh2R684C/+ and Ezh2+/+ bone marrow mesenchymal stem cells (BM-MSCs) indicated collective dysregulation of the BMP pathway and osteoblast differentiation. Inhibition of the opposing H3K27 demethylases Kdm6a/6b substantially reversed the excessive osteogenesis in Ezh2R684C/+ cells both at the transcriptional and phenotypic levels. This supports both the ideas that writers and erasers of histone marks exist in a fine balance to maintain epigenome state, and that epigenetic modulating agents have therapeutic potential for the treatment of MDEMs.

Success and Pitfalls of Genetic Testing in Undiagnosed Diseases: Whole Exome Sequencing and Beyond

Novel approaches to uncover the molecular etiology of neurodevelopmental disorders (NDD) are highly needed. Even using a powerful tool such as whole exome sequencing (WES), the diagnostic process may still prove long and arduous due to the high clinical and genetic heterogeneity of these conditions. The main strategies to improve the diagnostic rate are based on family segregation, re-evaluation of the clinical features by reverse-phenotyping, re-analysis of unsolved NGS-based cases and epigenetic functional studies. In this article, we described three selected cases from a cohort of patients with NDD in which trio WES was applied, in order to underline the typical challenges encountered during the diagnostic process: (1) an ultra-rare condition caused by a missense variant in MEIS2, identified through the updated Solve-RD re-analysis; (2) a patient with Noonan-like features in which the NGS analysis revealed a novel variant in NIPBL causing Cornelia de Lange syndrome; and (3) a case with de novo variants in genes involved in the chromatin-remodeling complex, for which the study of the epigenetic signature excluded a pathogenic role. In this perspective, we aimed to (i) provide an example of the relevance of the genetic re-analysis of all unsolved cases through network projects on rare diseases; (ii) point out the role and the uncertainties of the reverse phenotyping in the interpretation of the genetic results; and (iii) describe the use of methylation signatures in neurodevelopmental syndromes for the validation of the variants of uncertain significance.

Epigenomic signatures reveal mechanistic clues and predictive markers for autism spectrum disorder

Autism spectrum disorder (ASD) comprises a heterogeneous group of neurodevelopmental outcomes in children with a commonality in deficits in social communication and language combined with repetitive behaviors and interests. The etiology of ASD is heterogeneous, as several hundred genes have been implicated as well as multiple in utero environmental exposures. Over the past two decades, epigenetic investigations, including DNA methylation, have emerged as a novel way to capture the complex interface of multivariate ASD etiologies. More recently, epigenome-wide association studies using human brain and surrogate accessible tissues have revealed some convergent genes that are epigenetically altered in ASD, many of which overlap with known genetic risk factors. Unlike transcriptomes, epigenomic signatures defined by DNA methylation from surrogate tissues such as placenta and cord blood can reflect past differences in fetal brain gene transcription, transcription factor binding, and chromatin. For example, the discovery of NHIP (neuronal hypoxia inducible, placenta associated) through an epigenome-wide association in placenta, identified a common genetic risk for ASD that was modified by prenatal vitamin use. While epigenomic signatures are distinct between different genetic syndromic causes of ASD, bivalent chromatin and some convergent gene pathways are consistently epigenetically altered in both syndromic and idiopathic ASD, as well as some environmental exposures. Together, these epigenomic signatures hold promising clues towards improved early prediction and prevention of ASD as well genes and gene pathways to target for pharmacological interventions. Future advancements in single cell and multi-omic technologies, machine learning, as well as non-invasive screening of epigenomic signatures during pregnancy or newborn periods are expected to continue to impact the translatability of the recent discoveries in epigenomics to precision public health.

In individuals with Williams syndrome, dysregulation of methylation in non-coding regions of neuronal and oligodendrocyte DNA is associated with pathology and cortical development

Williams syndrome (WS) is a neurodevelopmental disorder caused by a heterozygous micro-deletion in the WS critical region (WSCR) and is characterized by hyper-sociability and neurocognitive abnormalities. Nonetheless, whether and to what extent WSCR deletion leads to epigenetic modifications in the brain and induces pathological outcomes remains largely unknown. By examining DNA methylation in frontal cortex, we revealed genome-wide disruption in the methylome of individuals with WS, as compared to typically developed (TD) controls. Surprisingly, differentially methylated sites were predominantly annotated as introns and intergenic loci and were found to be highly enriched around binding sites for transcription factors that regulate neuronal development, plasticity and cognition. Moreover, by utilizing enhancer-promoter interactome data, we confirmed that most of these loci function as active enhancers in the human brain or as target genes of transcriptional networks associated with myelination, oligodendrocyte (OL) differentiation, cognition and social behavior. Cell type-specific methylation analysis revealed aberrant patterns in the methylation of active enhancers in neurons and OLs, and important neuron-glia interactions that might be impaired in individuals with WS. Finally, comparison of methylation profiles from blood samples of individuals with WS and healthy controls, along with other data collected in this study, identified putative targets of endophenotypes associated with WS, which can be used to define brain-risk loci for WS outside the WSCR locus, as well as for other associated pathologies. In conclusion, our study illuminates the brain methylome landscape of individuals with WS and sheds light on how these aberrations might be involved in social behavior and physiological abnormalities. By extension, these results may lead to better diagnostics and more refined therapeutic targets for WS.

DNA methylation signature classification of rare disorders using publicly available methylation data

Disease-specific DNA methylation patterns (DNAm signatures) have been established for an increasing number of genetic disorders and represent a valuable tool for classification of genetic variants of uncertain significance (VUS). Sample size and batch effects are critical issues for establishing DNAm signatures, but their impact on the sensitivity and specificity of an already established DNAm signature has not previously been tested. Here, we assessed whether publicly available DNAm data can be employed to generate a binary machine learning classifier for VUS classification, and used variants in KMT2D, the gene associated with Kabuki syndrome, together with an existing DNAm signature as proof-of-concept. Using publicly available methylation data for training, a classifier for KMT2D variants was generated, and individuals with molecularly confirmed Kabuki syndrome and unaffected individuals could be correctly classified. The present study documents the clinical utility of a robust DNAm signature even for few affected individuals, and most importantly, underlines the importance of data sharing for improved diagnosis of rare genetic disorders.

Beyond the exome: what's next in diagnostic testing for Mendelian conditions

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order and emerging technologies, such as optical genome mapping and long-read DNA or RNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to a consortium such as GREGoR, which is focused on elucidating the underlying cause of rare unsolved genetic disorders.

Delineation of a KDM2B-related neurodevelopmental disorder and its associated DNA methylation signature

PURPOSE

Pathogenic variants in genes involved in the epigenetic machinery are an emerging cause of neurodevelopment disorders (NDDs). Lysine-demethylase 2B (KDM2B) encodes an epigenetic regulator and mouse models suggest an important role during development. We set out to determine whether KDM2B variants are associated with NDD.

METHODS

Through international collaborations, we collected data on individuals with heterozygous KDM2B variants. We applied methylation arrays on peripheral blood DNA samples to determine a KDM2B associated epigenetic signature.

RESULTS

We recruited a total of 27 individuals with heterozygous variants in KDM2B. We present evidence, including a shared epigenetic signature, to support a pathogenic classification of 15 KDM2B variants and identify the CxxC domain as a mutational hotspot. Both loss-of-function and CxxC-domain missense variants present with a specific subepisignature. Moreover, the KDM2B episignature was identified in the context of a dual molecular diagnosis in multiple individuals. Our efforts resulted in a cohort of 21 individuals with heterozygous (likely) pathogenic variants. Individuals in this cohort present with developmental delay and/or intellectual disability; autism; attention deficit disorder/attention deficit hyperactivity disorder; congenital organ anomalies mainly of the heart, eyes, and urogenital system; and subtle facial dysmorphism.

CONCLUSION

Pathogenic heterozygous variants in KDM2B are associated with NDD and a specific epigenetic signature detectable in peripheral blood.

Molecular characterization of an embryonal rhabdomyosarcoma occurring in a patient with Kabuki syndrome: report and literature review in the light of tumor predisposition syndromes

Kabuki syndrome is a well-recognized syndrome characterized by facial dysmorphism and developmental delay/intellectual disability and in the majority of patients a germline variant in KMT2D is found. As somatic KMT2D variants can be found in 5-10% of tumors a tumor predisposition in Kabuki syndrome is discussed. So far less than 20 patients with Kabuki syndrome and a concomitant malignancy have been published. Here we report on a female patient with Kabuki syndrome and a c.2558_2559delCT germline variant in KMT2D who developed an embryonal rhabdomyosarcoma (ERMS) at 10 years. On tumor tissue we performed DNA-methylation profiling and exome sequencing (ES). Copy number analyses revealed aneuploidies typical for ERMS including (partial) gains of chromosomes 2, 3, 7, 8, 12, 15, and 20 and 3 focal deletions of chromosome 11p. DNA methylation profiling mapped the case to ERMS by a DNA methylation-based sarcoma classifier. Sequencing suggested gain of the wild-type KMT2D allele in the trisomy 12. Including our patient literature review identified 18 patients with Kabuki syndrome and a malignancy. Overall, the landscape of malignancies in patients with Kabuki syndrome was reminiscent of that of the pediatric population in general. Histopathological and molecular data were only infrequently reported and no report included next generation sequencing and/or DNA-methylation profiling. Although we found no strong arguments pointing towards KS as a tumor predisposition syndrome, based on the small numbers any relation cannot be fully excluded. Further planned studies including profiling of additional tumors and long term follow-up of KS-patients into adulthood could provide further insights.

The discovery of the DNA methylation episignature for Duchenne muscular dystrophy

Duchenne Muscular Dystrophy (DMD) is an X-linked recessive neuromuscular disorder characterized by progressive muscle weakness due to loss of function mutations in the dystrophin gene. Variation in clinical presentation, the rate of disease progression, and treatment responsiveness have been observed amongst DMD patients, suggesting that factors beyond the loss of dystrophin may contribute to DMD pathophysiology. Epigenetic mechanisms are becoming recognized as important factors implicated in the etiology and progression of various diseases. A growing number of genetic syndromes have been associated with unique genomic DNA methylation patterns (called "episignatures") that can be used for diagnostic testing and as disease biomarkers. To further investigate DMD pathophysiology, we assessed the genome-wide DNA methylation profiles of peripheral blood from 36 patients with DMD using the combination of Illumina Infinium Methylation EPIC bead chip array and EpiSign technology. We identified a unique episignature for DMD that whose specificity was confirmed in relation other neurodevelopmental disorders with known episignatures. By modeling the DMD episignature, we developed a new DMD episignature biomarker and provided novel insights into the molecular pathogenesis of this disorder, which have the potential to advance more effective, personalized approaches to DMD care.

DNA methylation episignature for Witteveen-Kolk syndrome due to SIN3A haploinsufficiency

PURPOSE

Witteveen-Kolk syndrome (WITKOS) is a rare, autosomal dominant neurodevelopmental disorder caused by heterozygous loss-of-function alterations in the SIN3A gene. WITKOS has variable expressivity that commonly overlaps with other neurodevelopmental disorders. In this study, we characterized a distinct DNA methylation epigenetic signature (episignature) distinguishing WITKOS from unaffected individuals as well as individuals with other neurodevelopmental disorders with episignatures and described 9 previously unpublished individuals with SIN3A haploinsufficiency.

METHODS

We studied the phenotypic characteristics and the genome-wide DNA methylation in the peripheral blood samples of 20 individuals with heterozygous alterations in SIN3A. A total of 14 samples were used for the identification of the episignature and building of a predictive diagnostic biomarker, whereas the diagnostic model was used to investigate the methylation pattern of the remaining 6 samples.

RESULTS

A predominantly hypomethylated DNA methylation profile specific to WITKOS was identified, and the classifier model was able to diagnose a previously unresolved test case. The episignature was sensitive enough to detect individuals with varying degrees of phenotypic severity carrying SIN3A haploinsufficient variants.

CONCLUSION

We identified a novel, robust episignature in WITKOS due to SIN3A haploinsufficiency. This episignature has the potential to aid identification and diagnosis of individuals with WITKOS.

Genetically unresolved case of Rauch-Steindl syndrome diagnosed by its wolf-hirschhorn associated DNA methylation episignature

Wolf-Hirschhorn syndrome (WHS) is caused by deletion of a critical region of the short arm of chromosome 4. Clinical features of WHS include distinct dysmorphic facial features, growth restriction, developmental delay, intellectual disability, epilepsy, and other malformations. The NSD2 gene localizes within this critical region along with several other genes. Pathogenic variants in NSD2 cause Rauch-Steindl (RAUST) syndrome. Clinical features of RAUST syndrome partially overlap with WHS, however epilepsy and the recognizable facial gestalt are not observed. Here, we report a case of a young boy who presented with developmental delay, dysmorphic features and short stature. After negative chromosomal microarray and whole exome sequencing, genomic DNA methylation episignature analysis was performed. Episignatures are sensitive and specific genome-wide DNA methylation patterns associated with a growing number of rare disorders. The patient was positive for the WHS episignature. Reanalysis of the patient's exome data identified a previously undetected frameshift variant in NSD2, leading to a diagnosis of RAUST. This report demonstrates the clinical utility of DNA methylation episignature analysis for unresolved patients, and provides insight into the overlapping pathology between WHS and RAUST as demonstrated by the similarities in their genomic DNA methylation profiles.

DNA Methylation Analysis Reveals Distinct Patterns in Satellite Cell-Derived Myogenic Progenitor Cells of Subjects with Spastic Cerebral Palsy

Spastic type cerebral palsy (CP) is a complex neuromuscular disorder that involves altered skeletal muscle microanatomy and growth, but little is known about the mechanisms contributing to muscle pathophysiology and dysfunction. Traditional genomic approaches have provided limited insight regarding disease onset and severity, but recent epigenomic studies indicate that DNA methylation patterns can be altered in CP. Here, we examined whether a diagnosis of spastic CP is associated with intrinsic DNA methylation differences in myoblasts and myotubes derived from muscle resident stem cell populations (satellite cells; SCs). Twelve subjects were enrolled (6 CP; 6 control) with informed consent/assent. Skeletal muscle biopsies were obtained during orthopedic surgeries, and SCs were isolated and cultured to establish patient-specific myoblast cell lines capable of proliferation and differentiation in culture. DNA methylation analyses indicated significant differences at 525 individual CpG sites in proliferating SC-derived myoblasts (MB) and 1774 CpG sites in differentiating SC-derived myotubes (MT). Of these, 79 CpG sites were common in both culture types. The distribution of differentially methylated 1 Mbp chromosomal segments indicated distinct regional hypo- and hyper-methylation patterns, and significant enrichment of differentially methylated sites on chromosomes 12, 13, 14, 15, 18, and 20. Average methylation load across 2000 bp regions flanking transcriptional start sites was significantly different in 3 genes in MBs, and 10 genes in MTs. SC derived MBs isolated from study participants with spastic CP exhibited fundamental differences in DNA methylation compared to controls at multiple levels of organization that may reveal new targets for studies of mechanisms contributing to muscle dysregulation in spastic CP.

Episignature Mapping of TRIP12 Provides Functional Insight into Clark-Baraitser Syndrome

Clark-Baraitser syndrome is a rare autosomal dominant intellectual disability syndrome caused by pathogenic variants in the TRIP12 (Thyroid Hormone Receptor Interactor 12) gene. TRIP12 encodes an E3 ligase in the ubiquitin pathway. The ubiquitin pathway includes activating E1, conjugating E2 and ligating E3 enzymes which regulate the breakdown and sorting of proteins. This enzymatic pathway is crucial for physiological processes. A significant proportion of TRIP12 variants are currently classified as variants of unknown significance (VUS). Episignatures have been shown to represent a powerful diagnostic tool to resolve inconclusive genetic findings for Mendelian disorders and to re-classify VUSs. Here, we show the results of DNA methylation episignature analysis in 32 individuals with pathogenic, likely pathogenic and VUS variants in TRIP12. We identified a specific and sensitive DNA methylation (DNAm) episignature associated with pathogenic TRIP12 variants, establishing its utility as a clinical biomarker for Clark-Baraitser syndrome. In addition, we performed analysis of differentially methylated regions as well as functional correlation of the TRIP12 genome-wide methylation profile with the profiles of 56 additional neurodevelopmental disorders.

Disease classification for whole-blood DNA methylation: Meta-analysis, missing values imputation, and XAI

BACKGROUND

DNA methylation has a significant effect on gene expression and can be associated with various diseases. Meta-analysis of available DNA methylation datasets requires development of a specific workflow for joint data processing.

RESULTS

We propose a comprehensive approach of combined DNA methylation datasets to classify controls and patients. The solution includes data harmonization, construction of machine learning classification models, dimensionality reduction of models, imputation of missing values, and explanation of model predictions by explainable artificial intelligence (XAI) algorithms. We show that harmonization can improve classification accuracy by up to 20% when preprocessing methods of the training and test datasets are different. The best accuracy results were obtained with tree ensembles, reaching above 95% for Parkinson's disease. Dimensionality reduction can substantially decrease the number of features, without detriment to the classification accuracy. The best imputation methods achieve almost the same classification accuracy for data with missing values as for the original data. XAI approaches have allowed us to explain model predictions from both populational and individual perspectives.

CONCLUSIONS

We propose a methodologically valid and comprehensive approach to the classification of healthy individuals and patients with various diseases based on whole-blood DNA methylation data using Parkinson's disease and schizophrenia as examples. The proposed algorithm works better for the former pathology, characterized by a complex set of symptoms. It allows to solve data harmonization problems for meta-analysis of many different datasets, impute missing values, and build classification models of small dimensionality.

Identification of unique DNA methylation sites in Kabuki syndrome using whole genome bisulfite sequencing and targeted hybridization capture followed by enzymatic methylation sequencing

BACKGROUND

Kabuki syndrome (KS) is a congenital malformation syndrome caused by mutations in the KMT2D and KDM6A genes that encode histone modification enzymes. Although KS is considered a single gene disorder, its symptoms vary widely. Recently, disease-specific DNA methylation patterns, or episignatures, have been recognized and used as a diagnostic tool for KS. Because of various crosstalk mechanisms between histone modifications and DNA methylation, DNA methylation analysis may have high potential for investigations into the pathogenesis of KS.

RESULTS

In this study, we investigated altered CpG-methylation sites that were specific to KS to find important genes associated with the various phenotypes or pathogenesis of KS. Whole genome bisulfite sequencing (WGBS) was performed to select target CpG islands, and enzymatic conversion technology was applied after hybridization capture to confirm KS-specific episignatures of 130 selected differently methylated target regions (DMTRs) in DNA samples from the 65 participants, 31 patients with KS and 34 unaffected individuals, in this study. We identified 26 candidate genes in 22 DMTRs that may be associated with KS. Our results indicate that disease-specific methylation sites can be identified from a small number of WGBS samples, and hybridization capture followed by enzymatic methylation sequencing can simultaneously test the sites.

CONCLUSIONS

Although DNA methylation can be tissue-specific, our results suggest that methylation profiling of DNA extracted from peripheral blood may be a powerful approach to study the pathogenesis of diseases.

TMEM63C mutations cause mitochondrial morphology defects and underlie hereditary spastic paraplegia

The hereditary spastic paraplegias (HSP) are among the most genetically diverse of all Mendelian disorders. They comprise a large group of neurodegenerative diseases that may be divided into 'pure HSP' in forms of the disease primarily entailing progressive lower-limb weakness and spasticity, and 'complex HSP' when these features are accompanied by other neurological (or non-neurological) clinical signs. Here, we identified biallelic variants in the transmembrane protein 63C (TMEM63C) gene, encoding a predicted osmosensitive calcium-permeable cation channel, in individuals with hereditary spastic paraplegias associated with mild intellectual disability in some, but not all cases. Biochemical and microscopy analyses revealed that TMEM63C is an endoplasmic reticulum-localized protein, which is particularly enriched at mitochondria-endoplasmic reticulum contact sites. Functional in cellula studies indicate a role for TMEM63C in regulating both endoplasmic reticulum and mitochondrial morphologies. Together, these findings identify autosomal recessive TMEM63C variants as a cause of pure and complex HSP and add to the growing evidence of a fundamental pathomolecular role of perturbed mitochondrial-endoplasmic reticulum dynamics in motor neurone degenerative diseases.

Functional correlation of genome-wide DNA methylation profiles in genetic neurodevelopmental disorders

An expanding range of genetic syndromes are characterized by genome-wide disruptions in DNA methylation profiles referred to as episignatures. Episignatures are distinct, highly sensitive and specific biomarkers that have recently been applied in clinical diagnosis of genetic syndromes. Episignatures are contained within the broader disorder-specific genome-wide DNA methylation changes which can share significant overlap amongst different conditions. In this study we performed functional genomic assessment and comparison of disorder-specific and overlapping genome-wide DNA methylation changes related to 65 genetic syndromes with previously described episignatures. We demonstrate evidence of disorder-specific and recurring genome-wide differentially methylated probes (DMPs) and regions (DMRs). The overall distribution of DMPs and DMRs across the majority of the neurodevelopmental genetic syndromes analyzed showed substantial enrichment in gene promoters and CpG islands, and under-representation of the more variable intergenic regions. Analysis showed significant enrichment of the DMPs and DMRs in gene pathways and processes related to neurodevelopment, including neurogenesis, synaptic signaling and synaptic transmission. This study expands beyond the diagnostic utility of DNA methylation episignatures by demonstrating correlation between the function of the mutated genes and the consequent genomic DNA methylation profiles as a key functional element in the molecular etiology of genetic neurodevelopmental disorders. This article is protected by copyright. All rights reserved.

DNA Methylation Signature for JARID2-Neurodevelopmental Syndrome

JARID2 (Jumonji, AT Rich Interactive Domain 2) pathogenic variants cause a neurodevelopmental syndrome, that is characterized by developmental delay, cognitive impairment, hypotonia, autistic features, behavior abnormalities and dysmorphic facial features. JARID2 encodes a transcriptional repressor protein that regulates the activity of various histone methyltransferase complexes. However, the molecular etiology is not fully understood, and JARID2-neurodevelopmental syndrome may vary in its typical clinical phenotype. In addition, the detection of variants of uncertain significance (VUSs) often results in a delay of final diagnosis which could hamper the appropriate care. In this study we aim to detect a specific and sensitive DNA methylation signature for JARID2-neurodevelopmental syndrome. Peripheral blood DNA methylation profiles from 56 control subjects, 8 patients with (likely) pathogenic JARID2 variants and 3 patients with JARID2 VUSs were analyzed. DNA methylation analysis indicated a clear and robust separation between patients with (likely) pathogenic variants and controls. A binary model capable of classifying patients with the JARID2-neurodevelopmental syndrome was constructed on the basis of the identified episignature. Patients carrying VUSs clustered with the control group. We identified a distinct DNA methylation signature associated with JARID2-neurodevelopmental syndrome, establishing its utility as a biomarker for this syndrome and expanding the EpiSign diagnostic test.

DNA Methylation Episignatures in Neurodevelopmental Disorders Associated with Large Structural Copy Number Variants: Clinical Implications

Large structural chromosomal deletions and duplications, referred to as copy number variants (CNVs), play a role in the pathogenesis of neurodevelopmental disorders (NDDs) through effects on gene dosage. This review focuses on our current understanding of genomic disorders that arise from large structural chromosome rearrangements in patients with NDDs, as well as difficulties in overlap of clinical presentation and molecular diagnosis. We discuss the implications of epigenetics, specifically DNA methylation (DNAm), in NDDs and genomic disorders, and consider the implications and clinical impact of copy number and genomic DNAm testing in patients with suspected genetic NDDs. We summarize evidence of global methylation episignatures in CNV-associated disorders that can be used in the diagnostic pathway and may provide insights into the molecular pathogenesis of genomic disorders. Finally, we discuss the potential for combining CNV and DNAm assessment into a single diagnostic assay.

NSD1 mutations deregulate transcription and DNA methylation of bivalent developmental genes in Sotos syndrome

Sotos syndrome (SS), the most common overgrowth with intellectual disability (OGID) disorder, is caused by inactivating germline mutations of NSD1, which encodes a histone H3 lysine 36 methyltransferase. To understand how NSD1 inactivation deregulates transcription and DNA methylation (DNAm), and to explore how these abnormalities affect human development, we profiled transcription and DNAm in SS patients and healthy control individuals. We identified a transcriptional signature that distinguishes individuals with SS from controls and was also deregulated in NSD1-mutated cancers. Most abnormally expressed genes displayed reduced expression in SS; these downregulated genes consisted mostly of bivalent genes and were enriched for regulators of development and neural synapse function. DNA hypomethylation was strongly enriched within promoters of transcriptionally deregulated genes: overexpressed genes displayed hypomethylation at their transcription start sites while underexpressed genes featured hypomethylation at polycomb binding sites within their promoter CpG island shores. SS patients featured accelerated molecular aging at the levels of both transcription and DNAm. Overall, these findings indicate that NSD1-deposited H3K36 methylation regulates transcription by directing promoter DNA methylation, partially by repressing polycomb repressive complex 2 (PRC2) activity. These findings could explain the phenotypic similarity of SS to OGID disorders that are caused by mutations in PRC2 complex-encoding genes.

Altered DNA methylation at age-associated CpG sites in children with growth disorders: impact on age estimation?

Age estimation based on DNA methylation (DNAm) can be applied to children, adolescents and adults, but many CG dinucleotides (CpGs) exhibit different kinetics of age-associated DNAm across these age ranges. Furthermore, it is still unclear how growth disorders impact epigenetic age predictions, and this may be particularly relevant for a forensic application. In this study, we analyzed buccal mucosa samples from 95 healthy children and 104 children with different growth disorders. DNAm was analysed by pyrosequencing for 22 CpGs in the genes PDE4C, ELOVL2, RPA2, EDARADD and DDO. The relationship between DNAm and age in healthy children was tested by Spearman's rank correlation. Differences in DNAm between the groups "healthy children" and the (sub-)groups of children with growth disorders were tested by ANCOVA. Models for age estimation were trained (1) based on the data from 11 CpGs with a close correlation between DNAm and age (R ≥ 0.75) and (2) on five CpGs that also did not present significant differences in DNAm between healthy and diseased children. Statistical analysis revealed significant differences between the healthy group and the group with growth disorders (11 CpGs), the subgroup with a short stature (12 CpGs) and the non-short stature subgroup (three CpGs). The results are in line with the assumption of an epigenetic regulation of height-influencing genes. Age predictors trained on 11 CpGs with high correlations between DNAm and age revealed higher mean absolute errors (MAEs) in the group of growth disorders (mean MAE 2.21 years versus MAE 1.79 in the healthy group) as well as in the short stature (sub-)groups; furthermore, there was a clear tendency for overestimation of ages in all growth disorder groups (mean age deviations: total growth disorder group 1.85 years, short stature group 1.99 years). Age estimates on samples from children with growth disorders were more precise when using a model containing only the five CpGs that did not present significant differences in DNAm between healthy and diseased children (mean age deviations: total growth disorder group 1.45 years, short stature group 1.66 years). The results suggest that CpGs in genes involved in processes relevant for growth and development should be avoided in age prediction models for children since they may be sensitive for alterations in the DNAm pattern in cases of growth disorders.