CHARGE and Kabuki Syndromes: Gene-Specific DNA Methylation Signatures Identify Epigenetic Mechanisms Linking These Clinically Overlapping Conditions

Artificial intelligence-driven genotype-epigenotype-phenotype approaches to resolve challenges in syndrome diagnostics

BACKGROUND

Decisions to split two or more phenotypic manifestations related to genetic variations within the same gene can be challenging, especially during the early stages of syndrome discovery. Genotype-based diagnostics with artificial intelligence (AI)-driven approaches using next-generation phenotyping (NGP) and DNA methylation (DNAm) can be utilized to expedite syndrome delineation within a single gene.

METHODS

We utilized an expanded cohort of 56 patients (22 previously unpublished individuals) with truncating variants in the MN1 gene and attempted different methods to assess plausible strategies to objectively delineate phenotypic differences between the C-Terminal Truncation (CTT) and N-Terminal Truncation (NTT) groups. This involved transcriptomics analysis on available patient fibroblast samples and AI-assisted approaches, including a new statistical method of GestaltMatcher on facial photos and blood DNAm analysis using a support vector machine (SVM) model.

FINDINGS

RNA-seq analysis was unable to show a significant difference in transcript expression despite our previous hypothesis that NTT variants would induce nonsense mediated decay. DNAm analysis on nine blood DNA samples revealed an episignature for the CTT group. In parallel, the new statistical method of GestaltMatcher objectively distinguished the CTT and NTT groups with a low requirement for cohort number. Validation of this approach was performed on syndromes with known DNAm signatures of SRCAP, SMARCA2 and ADNP to demonstrate the effectiveness of this approach.

INTERPRETATION

We demonstrate the potential of using AI-based technologies to leverage genotype, phenotype and epigenetics data in facilitating splitting decisions in diagnosis of syndromes with minimal sample requirement.

FUNDING

The specific funding of this article is provided in the acknowledgements section.

The Multi-State Epigenetic Pacemaker enables the identification of combinations of factors that influence DNA methylation

Epigenetic clocks, DNA methylation-based predictive models of chronological age, are often utilized to study aging associated biology. Despite their widespread use, these methods do not account for other factors that also contribute to the variability of DNA methylation data. For example, many CpG sites show strong sex-specific or cell-type-specific patterns that likely impact the predictions of epigenetic age. To overcome these limitations, we developed a multidimensional extension of the Epigenetic Pacemaker, the Multi-state Epigenetic Pacemaker (MSEPM). We show that the MSEPM is capable of accurately modeling multiple methylation-associated factors simultaneously, while also providing site-specific models that describe the per site relationship between methylation and these factors. We utilized the MSEPM with a large aggregate cohort of blood methylation data to construct models of the effects of age-, sex-, and cell-type heterogeneity on DNA methylation. We found that these models capture a large faction of the variability at thousands of DNA methylation sites. Moreover, this approach allows us to identify sites that are primarily affected by aging and no other factors. An analysis of these sites reveals that those that lose methylation over time are enriched for CTCF transcription factor chip peaks, while those that gain methylation over time are associated with bivalent promoters of genes that are not expressed in blood. These observations suggest mechanisms that underlie age-associated methylation changes and suggest that age-associated increases in methylation may not have strong functional consequences on cell states. In conclusion, the MSEPM is capable of accurately modeling multiple methylation-associated factors, and the models produced can illuminate site-specific combinations of factors that affect methylation dynamics.

The Arg99Gln Substitution in HNRNPC Is Associated with a Distinctive Clinical Phenotype Characterized by Facial Dysmorphism and Ocular and Cochlear Anomalies

Background/Objectives: Heterozygous variants in the heterogeneous nuclear ribonucleoprotein C gene (HNRNPC) have recently been reported to cause intellectual developmental disorder-74 (MRD74), a neurodevelopmental disorder with no recurrent diagnostic handles. Affected individuals show variable, non-specific, and subtle dysmorphic features. The degree of developmental delay (DD)/intellectual disability (ID) is also wide, ranging from mild to severe. The mutational spectrum is relatively broad with exon deletions and splice site and frameshift variants distributed along the entire length of the gene leading to HNRNPC loss of function. Only two missense changes located within the RNA-binding motif (RBM) and adjacent linker region of the more abundant isoform (Arg64Trp and Arg99Gln) have been described. Notably, the Arg99Gln amino acid substitution was reported in a subject presenting with a more complex and unique clinical phenotype characterized by distinctive facial features, DD/ID, cochlear aplasia, and bilateral colobomatous microphthalmia, suggesting the possible occurrence of phenotypic heterogeneity. Results: Here, we report the second individual carrying the Arg99Gln change in HNRNPC and having clinical features with a significant overlap with the peculiar phenotype of the previously described subject, supporting the occurrence of a genotype-phenotype correlation. Conclusions: Due to the concomitant occurrence of ocular and cochlear involvement as recognizable diagnostic handles, we propose that the HNRNPCArg99Gln-related phenotype should be considered as a potential differential diagnosis in subjects with ID and major signs of CHARGE syndrome not fulfilling the minimum criteria for a clinical diagnosis.

Polycomb-associated and Trithorax-associated developmental conditions-phenotypic convergence and heterogeneity

Polycomb group (PcG) and Trithorax group (TrxG) complexes represent two major components of the epigenetic machinery. This study aimed to delineate phenotypic similarities and differences across developmental conditions arising from rare variants in PcG and TrxG genes, using data-driven approaches. 462 patients with a PcG or TrxG-associated condition were identified in the DECIPHER dataset. We analysed Human Phenotype Ontology (HPO) data to identify phenotypes enriched in this group, in comparison to other monogenic conditions within DECIPHER. We then assessed phenotypic relationships between single gene diagnoses within the PcG and TrxG group, by applying semantic similarity analysis and hierarchical clustering. Finally, we analysed patient-level phenotypic heterogeneity in this group, irrespective of specific genetic diagnosis, by applying the same clustering approach. Collectively, PcG/TrxG diagnoses were associated with increased reporting of HPO terms relating to integument, growth, head and neck, limb and digestive abnormalities. Gene group analysis identified three multi-gene clusters differentiated by microcephaly, limb/digit dysmorphologies, growth abnormalities and atypical behavioural phenotypes. Patient-level analysis identified two large clusters differentiated by neurodevelopmental abnormalities and facial dysmorphologies respectively, as well as smaller clusters associated with more specific phenotypes including behavioural characteristics, eye abnormalities, growth abnormalities and skull dysmorphologies. Importantly, patient-level phenotypic clusters did not align with genetic diagnoses. Data-driven approaches can highlight pathway-level and gene-level phenotypic convergences, and individual-level phenotypic heterogeneities. Future studies are needed to understand the multi-level mechanisms contributing to both convergence and variability within this population, and to extend data collection and analyses to later-emerging health characteristics.

Expanding Upon Genomics in Rare Diseases: Epigenomic Insights

DNA methylation is an essential epigenetic modification that plays a crucial role in regulating gene expression and maintaining genomic stability. With the advancement in sequencing technology, methylation studies have provided valuable insights into the diagnosis of rare diseases through the various identification of episignatures, epivariation, epioutliers, and allele-specific methylation. However, current methylation studies are not without limitations. This mini-review explores the current understanding of DNA methylation in rare diseases, highlighting the key mechanisms and diagnostic potential, and emphasizing the need for advanced methodologies and integrative approaches to enhance the understanding of disease progression and design more personable treatment for patients, given the nature of rare diseases.

Higher order interaction analysis quantifies coordination in the epigenome revealing novel biological relationships in Kabuki syndrome

Complex direct and indirect relationships between multiple variables, termed higher order interactions (HOIs), are characteristics of all natural systems. Traditional differential and network analyses fail to account for the omic datasets richness and miss HOIs. We investigated peripheral blood DNA methylation data from Kabuki syndrome type 1 (KS1) and control individuals, identified 2,002 differentially methylated points (DMPs), and inferred 17 differentially methylated regions, which represent only 189 DMPs. We applied hypergraph models to measure HOIs on all the CpGs and revealed differences in the coordination of DMPs with lower entropy and higher coordination of the peripheral epigenome in KS1 implying reduced network complexity. Hypergraphs also capture epigenomic trans-relationships, and identify biologically relevant pathways that escape the standard analyses. These findings construct the basis of a suitable model for the analysis of organization in the epigenome in rare diseases, which can be applied to investigate mechanism in big data.

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 changes in the genome of patients with hypogonadotropic hypogonadism

Although some Mendelian neurodevelopmental disorders have been shown to entail specific DNA methylation changes designated as epi-signatures, it remains unknown whether epi-signatures are consistent features of other genetic disorders. Here, we analyzed DNA methylation profiles of patients with hypogonadotropic hypogonadism (HH), a rare neuroendocrine disorder typically caused by monogenic or oligogenic mutations. First, we performed microarray-based genome-wide methylation analyses of nine patients with HH due to ANOS1, SOX2, or SOX10 variants and 12 control individuals. The results showed that 1118 probes were differentially methylated in one or more patients. The differentially methylated probes were highly variable among patients. No significant methylation changes were observed in genes functionally associated with ANOS1, SOX2, or SOX10. Then, we performed pyrosequencing of six selected CpG sites in the nine patients and 35 additional HH patients. The results of the patients were compared with those of 48 fertile men. There were no common methylation changes among these patients, with the exception of hypermethylation of two CpG sites in the ZNF245 promoter of three patients. Hypermethylation of the promoter has previously been reported as a very rare epigenetic polymorphism in the general population. These results indicate that genomes of HH patients have considerable DNA methylation changes; however, these changes are more likely to be physiological epigenetic variations than disease-specific epi-signatures. Our data suggest a possible association between hypermethylation of the ZNF254 promoter and HH, which needs to be examined in future studies.

Knockout, Knockdown, and the Schrödinger Paradox: Genetic Immunity to Phenotypic Recapitulation in Zebrafish

Previous research has highlighted significant phenotypic discrepancies between knockout and knockdown approaches in zebrafish, raising concerns about the reliability of these methods. However, our study suggests that these differences are not as pronounced as was once believed. By carefully examining the roles of maternal and zygotic gene contributions, we demonstrate that these factors significantly influence phenotypic outcomes, often accounting for the observed discrepancies. Our findings emphasize that morpholinos, despite their potential off-target effects, can be effective tools when used with rigorous controls. We introduce the concept of graded maternal contribution, which explains how the uneven distribution of maternal mRNA and proteins during gametogenesis impacts phenotypic variability. Our research categorizes genes into three types-susceptible, immune, and "Schrödinger" (conditional)-based on their phenotypic expression and interaction with genetic compensation mechanisms. This distinction provides new insights into the paradoxical outcomes observed in genetic studies. Ultimately, our work underscores the importance of considering both maternal and zygotic contributions, alongside rigorous experimental controls, to accurately interpret gene function and the mechanisms underlying disease. This study advocates for the continued use of morpholinos in conjunction with advanced genetic tools like CRISPR/Cas9, stressing the need for a meticulous experimental design to optimize the utility of zebrafish in genetic research and therapeutic development.

Pathogenic variants in KMT2C result in a neurodevelopmental disorder distinct from Kleefstra and Kabuki syndromes

Trithorax-related H3K4 methyltransferases, KMT2C and KMT2D, are critical epigenetic modifiers. Haploinsufficiency of KMT2C was only recently recognized as a cause of neurodevelopmental disorder (NDD), so the clinical and molecular spectrums of the KMT2C-related NDD (now designated as Kleefstra syndrome 2) are largely unknown. We ascertained 98 individuals with rare KMT2C variants, including 75 with protein-truncating variants (PTVs). Notably, ∼15% of KMT2C PTVs were inherited. Although the most highly expressed KMT2C transcript consists of only the last four exons, pathogenic PTVs were found in almost all the exons of this large gene. KMT2C variant interpretation can be challenging due to segmental duplications and clonal hematopoesis-induced artifacts. Using samples from 27 affected individuals, divided into discovery and validation cohorts, we generated a moderate strength disorder-specific KMT2C DNA methylation (DNAm) signature and demonstrate its utility in classifying non-truncating variants. Based on 81 individuals with pathogenic/likely pathogenic variants, we demonstrate that the KMT2C-related NDD is characterized by developmental delay, intellectual disability, behavioral and psychiatric problems, hypotonia, seizures, short stature, and other comorbidities. The facial module of PhenoScore, applied to photographs of 34 affected individuals, reveals that the KMT2C-related facial gestalt is significantly different from the general NDD population. Finally, using PhenoScore and DNAm signatures, we demonstrate that the KMT2C-related NDD is clinically and epigenetically distinct from Kleefstra and Kabuki syndromes. Overall, we define the clinical features, molecular spectrum, and DNAm signature of the KMT2C-related NDD and demonstrate they are distinct from Kleefstra and Kabuki syndromes highlighting the need to rename this condition.

Genomic Landscape of Branchio-Oto-Renal Syndrome through Whole-Genome Sequencing: A Single Rare Disease Center Experience in South Korea

Branchio-oto-renal (BOR) and branchio-otic (BO) syndromes are characterized by anomalies affecting the ears, often accompanied by hearing loss, as well as abnormalities in the branchial arches and renal system. These syndromes exhibit a broad spectrum of phenotypes and a complex genomic landscape, with significant contributions from the EYA1 gene and the SIX gene family, including SIX1 and SIX5. Due to their diverse phenotypic presentations, which can overlap with other genetic syndromes, molecular genetic confirmation is essential. As sequencing technologies advance, whole-genome sequencing (WGS) is increasingly used in rare disease diagnostics. We explored the genomic landscape of 23 unrelated Korean families with typical or atypical BOR/BO syndrome using a stepwise approach: targeted panel sequencing and exome sequencing (Step 1), multiplex ligation-dependent probe amplification (MLPA) with copy number variation screening (Step 2), and WGS (Step 3). Integrating WGS into our diagnostic pipeline detected structure variations, including cryptic inversion and complex genomic rearrangement, eventually enhancing the diagnostic yield to 91%. Our findings expand the genomic architecture of BOR/BO syndrome and highlight the need for WGS to address the genetic diagnosis of clinically heterogeneous rare diseases.

DNA methylation profiling in Kabuki syndrome: reclassification of germline KMT2D VUS and sensitivity in validating postzygotic mosaicism

Autosomal dominant Kabuki syndrome (KS) is a rare multiple congenital anomalies/neurodevelopmental disorder caused by heterozygous inactivating variants or structural rearrangements of the lysine-specific methyltransferase 2D (KMT2D) gene. While it is often recognizable due to a distinctive gestalt, the disorder is clinically variable, and a phenotypic scoring system has been introduced to help clinicians to reach a clinical diagnosis. The phenotype, however, can be less pronounced in some patients, including those carrying postzygotic mutations. The full spectrum of pathogenic variation in KMT2D has not fully been characterized, which may hamper the clinical classification of a portion of these variants. DNA methylation (DNAm) profiling has successfully been used as a tool to classify variants in genes associated with several neurodevelopmental disorders, including KS. In this work, we applied a KS-specific DNAm signature in a cohort of 13 individuals with KMT2D VUS and clinical features suggestive or overlapping with KS. We succeeded in correctly classifying all the tested individuals, confirming diagnosis for three subjects and rejecting the pathogenic role of 10 VUS in the context of KS. In the latter group, exome sequencing allowed to identify the genetic cause underlying the disorder in three subjects. By testing five individuals with postzygotic pathogenic KMT2D variants, we also provide evidence that DNAm profiling has power to recognize pathogenic variants at different levels of mosaicism, identifying 15% as the minimum threshold for which DNAm profiling can be applied as an informative diagnostic tool in KS mosaics.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically shows phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption-the root of the pathogenesis-is similar in the different disease-relevant cell types. This is possible in principle, because all these cell types are subject to effects from the same causative gene, which has the same kind of function (e.g., methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2 and find that the chromatin accessibility changes in neurons are mostly distinct from changes in B or T cells. This is not because the neuronal accessibility changes occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that regulatory elements disrupted in B/T cells do show chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators and suggest that blood-derived episignatures, although useful diagnostically, may not be well suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Growth deficiency in a mouse model of Kabuki syndrome 2 bears mechanistic similarities to Kabuki syndrome 1

Growth deficiency is a characteristic feature of both Kabuki syndrome 1 (KS1) and Kabuki syndrome 2 (KS2), Mendelian disorders of the epigenetic machinery with similar phenotypes but distinct genetic etiologies. We previously described skeletal growth deficiency in a mouse model of KS1 and further established that a Kmt2d-/- chondrocyte model of KS1 exhibits precocious differentiation. Here we characterized growth deficiency in a mouse model of KS2, Kdm6atm1d/+. We show that Kdm6atm1d/+ mice have decreased femur and tibia length compared to controls and exhibit abnormalities in cortical and trabecular bone structure. Kdm6atm1d/+ growth plates are also shorter, due to decreases in hypertrophic chondrocyte size and hypertrophic zone height. Given these disturbances in the growth plate, we generated Kdm6a-/- chondrogenic cell lines. Similar to our prior in vitro model of KS1, we found that Kdm6a-/- cells undergo premature, enhanced differentiation towards chondrocytes compared to Kdm6a+/+ controls. RNA-seq showed that Kdm6a-/- cells have a distinct transcriptomic profile that indicates dysregulation of cartilage development. Finally, we performed RNA-seq simultaneously on Kmt2d-/-, Kdm6a-/-, and control lines at Days 7 and 14 of differentiation. This revealed surprising resemblance in gene expression between Kmt2d-/- and Kdm6a-/- at both time points and indicates that the similarity in phenotype between KS1 and KS2 also exists at the transcriptional level.

The expanding diagnostic toolbox for rare genetic diseases

Genomic technologies, such as targeted, exome and short-read genome sequencing approaches, have revolutionized the care of patients with rare genetic diseases. However, more than half of patients remain without a diagnosis. Emerging approaches from research-based settings such as long-read genome sequencing and optical genome mapping hold promise for improving the identification of disease-causal genetic variants. In addition, new omic technologies that measure the transcriptome, epigenome, proteome or metabolome are showing great potential for variant interpretation. As genetic testing options rapidly expand, the clinical community needs to be mindful of their individual strengths and limitations, as well as remaining challenges, to select the appropriate diagnostic test, correctly interpret results and drive innovation to address insufficiencies. If used effectively - through truly integrative multi-omics approaches and data sharing - the resulting large quantities of data from these established and emerging technologies will greatly improve the interpretative power of genetic and genomic diagnostics for rare diseases.

The omics era: a nexus of untapped potential for Mendelian chromatinopathies

The OMICs cascade describes the hierarchical flow of information through biological systems. The epigenome sits at the apex of the cascade, thereby regulating the RNA and protein expression of the human genome and governs cellular identity and function. Genes that regulate the epigenome, termed epigenes, orchestrate complex biological signaling programs that drive human development. The broad expression patterns of epigenes during human development mean that pathogenic germline mutations in epigenes can lead to clinically significant multi-system malformations, developmental delay, intellectual disabilities, and stem cell dysfunction. In this review, we refer to germline developmental disorders caused by epigene mutation as "chromatinopathies". We curated the largest number of human chromatinopathies to date and our expanded approach more than doubled the number of established chromatinopathies to 179 disorders caused by 148 epigenes. Our study revealed that 20.6% (148/720) of epigenes cause at least one chromatinopathy. In this review, we highlight key examples in which OMICs approaches have been applied to chromatinopathy patient biospecimens to identify underlying disease pathogenesis. The rapidly evolving OMICs technologies that couple molecular biology with high-throughput sequencing or proteomics allow us to dissect out the causal mechanisms driving temporal-, cellular-, and tissue-specific expression. Using the full repertoire of data generated by the OMICs cascade to study chromatinopathies will provide invaluable insight into the developmental impact of these epigenes and point toward future precision targets for these rare disorders.

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.

A new blood DNA methylation signature for Koolen-de Vries syndrome: Classification of missense KANSL1 variants and comparison to fibroblast cells

Pathogenic variants in KANSL1 and 17q21.31 microdeletions are causative of Koolen-de Vries syndrome (KdVS), a neurodevelopmental syndrome with characteristic facial dysmorphia. Our previous work has shown that syndromic conditions caused by pathogenic variants in epigenetic regulatory genes have identifiable patterns of DNA methylation (DNAm) change: DNAm signatures or episignatures. Given the role of KANSL1 in histone acetylation, we tested whether variants underlying KdVS are associated with a DNAm signature. We profiled whole-blood DNAm for 13 individuals with KANSL1 variants, four individuals with 17q21.31 microdeletions, and 21 typically developing individuals, using Illumina's Infinium EPIC array. In this study, we identified a robust DNAm signature of 456 significant CpG sites in 8 individuals with KdVS, a pattern independently validated in an additional 7 individuals with KdVS. We also demonstrate the diagnostic utility of the signature and classify two KANSL1 VUS as well as four variants in individuals with atypical clinical presentation. Lastly, we investigated tissue-specific DNAm changes in fibroblast cells from individuals with KdVS. Collectively, our findings contribute to the understanding of the epigenetic landscape related to KdVS and aid in the diagnosis and classification of variants in this structurally complex genomic region.

A novel CHD7 variant in a chinese family with CHARGE syndrome

OBJECTIVE

CHARGE syndrome is a rare autosomal dominant (AD) multi-system disorder with a broad and variable clinical manifestation and occurs in approximately 1/10,000 newborns in the world. Mutations in the CHD7 gene are the genetic cause of over 90% of patients with typical CHARGE syndrome. The present study reported a novel variant in the CHD7 gene in a Chinese family with an abnormal fetus.

METHODS

Routine prenatal ultrasound screening showed fetal heart abnormality and left foot varus. Chromosomal microarray analysis (CMA) and fetus-parent whole-exome sequencing (trio-WES) were performed to determine the genetic cause of the fetus. The candidate variant was further verified using Sanger sequencing.

RESULTS

CMA analysis revealed normal results. However, WES analysis identified a de novo heterozygous variant of c.2919_2922del (NM_017780.4) on exon 11 of CHD7 gene, resulting in a premature truncation of the CHD7 protein (p.Gly975*). The variant was classified as Pathogenic (PVS1 + PS2_Moderate + PM2_Supporting) based on the ACMG guidelines. Combined with the clinical phenotype of fetal heart abnormalities, it was confirmed CHARGE syndrome.

CONCLUSION

We identified a novel heterozygous variant c.2919_2922del in CHD7 of a Chinese fetus with CHARGE syndrome, enriching the genotype-phenotype spectrum of CHD7. These results suggest that genetic testing could help facilitate prenatal diagnosis of CHARGE syndrome, thus promoting the appropriate genetic counseling.

Episignatures in practice: independent evaluation of published episignatures for the molecular diagnostics of ten neurodevelopmental disorders

Variants of uncertain significance (VUS) are a significant issue for the molecular diagnosis of rare diseases. The publication of episignatures as effective biomarkers of certain Mendelian neurodevelopmental disorders has raised hopes to help classify VUS. However, prediction abilities of most published episignatures have not been independently investigated yet, which is a prerequisite for an informed and rigorous use in a diagnostic setting. We generated DNA methylation data from 101 carriers of (likely) pathogenic variants in ten different genes, 57 VUS carriers, and 25 healthy controls. Combining published episignature information and new validation data with a k-nearest-neighbour classifier within a leave-one-out scheme, we provide unbiased specificity and sensitivity estimates for each of the signatures. Our procedure reached 100% specificity, but the sensitivities unexpectedly spanned a very large spectrum. While ATRX, DNMT3A, KMT2D, and NSD1 signatures displayed a 100% sensitivity, CREBBP-RSTS and one of the CHD8 signatures reached <40% sensitivity on our dataset. Remaining Cornelia de Lange syndrome, KMT2A, KDM5C and CHD7 signatures reached 70-100% sensitivity at best with unstable performances, suffering from heterogeneous methylation profiles among cases and rare discordant samples. Our results call for cautiousness and demonstrate that episignatures do not perform equally well. Some signatures are ready for confident use in a diagnostic setting. Yet, it is imperative to characterise the actual validity perimeter and interpretation of each episignature with the help of larger validation sample sizes and in a broader set of episignatures.

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.

Identifying epigenetic aging moderators using the epigenetic pacemaker

Epigenetic clocks are DNA methylation-based chronological age prediction models that are commonly employed to study age-related biology. The difference between the predicted and observed age is often interpreted as a form of biological age acceleration, and many studies have measured the impact of environmental and disease-associated factors on epigenetic age. Most epigenetic clocks are fit using approaches that minimize the error between the predicted and observed chronological age, and as a result, they may not accurately model the impact of factors that moderate the relationship between the actual and epigenetic age. Here, we compare epigenetic clocks that are constructed using penalized regression methods to an evolutionary framework of epigenetic aging with the epigenetic pacemaker (EPM), which directly models DNA methylation as a function of a time-dependent epigenetic state. In simulations, we show that the value of the epigenetic state is impacted by factors such as age, sex, and cell-type composition. Next, in a dataset aggregated from previous studies, we show that the epigenetic state is also moderated by sex and the cell type. Finally, we demonstrate that the epigenetic state is also moderated by toxins in a study on polybrominated biphenyl exposure. Thus, we find that the pacemaker provides a robust framework for the study of factors that impact epigenetic age acceleration and that the effect of these factors may be obscured in traditional clocks based on linear regression models.

KMT2D Deficiency Causes Sensorineural Hearing Loss in Mice and Humans

Individuals with Kabuki syndrome type 1 (KS1) often have hearing loss recognized in middle childhood. Current clinical dogma suggests that this phenotype is caused by frequent infections due to the immune deficiency in KS1 and/or secondary to structural abnormalities of the ear. To clarify some aspects of hearing loss, we collected information on hearing status from 21 individuals with KS1 and found that individuals have both sensorineural and conductive hearing loss, with the average age of presentation being 7 years. Our data suggest that while ear infections and structural abnormalities contribute to the observed hearing loss, these factors do not explain all loss. Using a KS1 mouse model, we found hearing abnormalities from hearing onset, as indicated by auditory brainstem response measurements. In contrast to mouse and human data for CHARGE syndrome, a disorder possessing overlapping clinical features with KS and a well-known cause of hearing loss and structural inner ear abnormalities, there are no apparent structural abnormalities of the cochlea in KS1 mice. The KS1 mice also display diminished distortion product otoacoustic emission levels, which suggests outer hair cell dysfunction. Combining these findings, our data suggests that KMT2D dysfunction causes sensorineural hearing loss compounded with external factors, such as infection.

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.

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.

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.

DNA methylation profiles in individuals with rare, atypical 7q11.23 CNVs correlate with GTF2I and GTF2IRD1 copy number

Williams-Beuren syndrome (WBS) and 7q11.23 duplication syndrome (Dup7) are rare neurodevelopmental disorders caused by deletion and duplication of a 1.5 Mb region that includes at least five genes with a known role in epigenetic regulation. We have shown that CNV of this chromosome segment causes dose-dependent, genome-wide changes in DNA methylation, but the specific genes driving these changes are unknown. We measured genome-wide whole blood DNA methylation in six participants with atypical CNV of 7q11.23 (three with deletions and three with duplications) using the Illumina HumanMethylation450k array and compared their profiles with those from groups of individuals with classic WBS or classic Dup7 and with typically developing (TD) controls. Across the top 1000 most variable positions we found that only the atypical rearrangements that changed the copy number of GTF2IRD1 and/or GTF2I (coding for the TFII-IRD1 and TFII-I proteins) clustered with their respective syndromic cohorts. This finding was supported by results from hierarchical clustering across a selection of differentially methylated CpGs, in addition to pyrosequencing validation. These findings suggest that CNV of the GTF2I genes at the telomeric end of the 7q11.23 interval is a key contributor to the large changes in DNA methylation that are seen in blood DNA from our WBS and Dup7 cohorts, compared to TD controls. Our findings suggest that members of the TFII-I protein family are involved in epigenetic processes that alter DNA methylation on a genome-wide level.

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.

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.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically exhibits phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption - the root of the pathogenesis - is similar in the different disease-relevant cell types. This is possible in principle, since all these cell-types are subject to effects from the same causative gene, that has the same kind of function (e.g. methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2, and find that the chromatin accessibility abnormalities in neurons are mostly distinct from those in B or T cells. This is not because the neuronal abnormalities occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that the regions disrupted in B/T cells do exhibit chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators, and suggest that blood-derived "episignatures" may not be well-suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

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.

Identification of a KDM6A somatic mutation responsible for Kabuki syndrome by excluding a conflicting KMT2D germline variant through episignature analysis

Kabuki syndrome (KS) is a congenital disorder caused by mutations in either KMT2D on chromosome 12 or KDM6A on chromosome X, encoding a lysine methyltransferase and a lysine demethylase, respectively. A 9-year-4-month-old male patient with a normal karyotype presented with KS and autism spectrum disorder. Genetic testing for KS was conducted by Sanger sequencing and episignature analysis using DNA methylation array data. The patient had a mosaic stop-gain variant in KDM6A and a heterozygous missense variant (rs201078160) in KMT2D. The KDM6A variant is expected to be deleterious. The KMT2D variant pathogenicity has been inconsistently reported in the ClinVar database. Using biobanking resources, we identified two heterozygous individuals possessing the rs201078160 variant. In a subsequent episignature analysis, the KS patient showed the KS episignature, but two control individuals with the rs201078160 variant did not. Our results indicate that the mosaic stop-gained variant in KDM6A, but not the rs201078160 variant in KMT2D, is responsible for the KS phenotype in the patient. This study further demonstrated the utility of DNA methylation information in diagnosing rare genetic diseases and emphasized the importance of a reference dataset containing both genotype and DNA methylation information.

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.

DNA methylation differences in monozygotic twins with Van der Woude syndrome

Introduction

Van der Woude Syndrome (VWS) is an autosomal dominant disorder responsible for 2% of all syndromic orofacial clefts (OFCs) with IRF6 being the primary causal gene (70%). Cases may present with lip pits and either cleft lip, cleft lip with cleft palate, or cleft palate, with marked phenotypic discordance even among individuals carrying the same mutation. This suggests that genetic or epigenetic modifiers may play additional roles in the syndrome's etiology and variability in expression. We report the first DNA methylation profiling of 2 pairs of monozygotic twins with VWS. Our goal is to explore epigenetic contributions to VWS etiology and variable phenotypic expressivity by comparing DNAm profiles in both twin pairs. While the mutations that cause VWS in these twins are known, the additional mechanism behind their phenotypic risk and variability in expression remains unclear.

Methods

We generated whole genome DNAm data for both twin pairs. Differentially methylated positions (DMPs) were selected based on: (1) a coefficient of variation in DNAm levels in unaffected individuals < 20%, and (2) intra-twin pair absolute difference in DNAm levels >5% (delta beta > | 0.05|). We then divided the DMPs in two subgroups for each twin pair for further analysis: (1) higher methylation levels in twin A (Twin A > Twin B); and (2) higher methylation levels in twin B (Twin B >Twin A).

Results and Discussion

Gene ontology analysis revealed a list of enriched genes that showed significant differential DNAm, including clef-associated genes. Among the cleft-associated genes, TP63 was the most significant hit (p=7.82E-12). Both twin pairs presented differential DNAm levels in CpG sites in/near TP63 (Twin 1A > Twin 1B and Twin 2A < Twin 2B). The genes TP63 and IRF6 function in a biological regulatory loop to coordinate epithelial proliferation and differentiation in a process that is critical for palatal fusion. The effects of the causal mutations in IRF6 can be further impacted by epigenetic dysregulation of IRF6 itself, or genes in its pathway. Our data shows evidence that changes in DNAm is a plausible mechanism that can lead to markedly distinct phenotypes, even among individuals carrying the same mutation.

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.

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.

Identification of differentially methylated regions in rare diseases from a single-patient perspective

BACKGROUND

DNA methylation (5-mC) is being widely recognized as an alternative in the detection of sequence variants in the diagnosis of some rare neurodevelopmental and imprinting disorders. Identification of alterations in DNA methylation plays an important role in the diagnosis and understanding of the etiology of those disorders. Canonical pipelines for the detection of differentially methylated regions (DMRs) usually rely on inter-group (e.g., case versus control) comparisons. However, these tools might perform suboptimally in the context of rare diseases and multilocus imprinting disturbances due to small cohort sizes and inter-patient heterogeneity. Therefore, there is a need to provide a simple but statistically robust pipeline for scientists and clinicians to perform differential methylation analyses at the single patient level as well as to evaluate how parameter fine-tuning may affect differentially methylated region detection.

RESULT

We implemented an improved statistical method to detect differentially methylated regions in correlated datasets based on the Z-score and empirical Brown aggregation methods from a single-patient perspective. To accurately assess the predictive power of our method, we generated semi-simulated data using a public control population of 521 samples and investigated how the size of the control population, methylation difference, and region size affect DMR detection. In addition, we validated the detection of methylation events in patients suffering from rare multi-locus imprinting disturbance and evaluated how this method could complement existing tools in the context of clinical diagnosis.

CONCLUSION

In this study, we present a robust statistical method to perform differential methylation analysis at the single patient level and describe its optimal parameters to increase DMRs identification performance. Finally, we show its diagnostic utility when applied to rare disorders.

Craniofacial and cardiac defects in chd7 zebrafish mutants mimic CHARGE syndrome

Congenital heart defects occur in almost 80% of patients with CHARGE syndrome, a sporadically occurring disease causing craniofacial and other abnormalities due to mutations in the CHD7 gene. Animal models have been generated to mimic CHARGE syndrome; however, heart defects are not extensively described in zebrafish disease models of CHARGE using morpholino injections or genetic mutants. Here, we describe the co-occurrence of craniofacial abnormalities and heart defects in zebrafish chd7 mutants. These mutant phenotypes are enhanced in the maternal zygotic mutant background. In the chd7 mutant fish, we found shortened craniofacial cartilages and extra cartilage formation. Furthermore, the length of the ventral aorta is altered in chd7 mutants. Many CHARGE patients have aortic arch anomalies. It should be noted that the aberrant branching of the first branchial arch artery is observed for the first time in chd7 fish mutants. To understand the cellular mechanism of CHARGE syndrome, neural crest cells (NCCs), that contribute to craniofacial and cardiovascular tissues, are examined using sox10:Cre lineage tracing. In contrast to its function in cranial NCCs, we found that the cardiac NCC-derived mural cells along the ventral aorta and aortic arch arteries are not affected in chd7 mutant fish. The chd7 fish mutants we generated recapitulate some of the craniofacial and cardiovascular phenotypes found in CHARGE patients and can be used to further determine the roles of CHD7.

Generation of DNA Methylation Signatures and Classification of Variants in Rare Neurodevelopmental Disorders Using EpigenCentral

There are more than 700 genes that encode proteins that function in epigenetic regulation and chromatin modification. Germline variants in these genes (typically heterozygous) are associated with rare neurodevelopmental disorders (NDDs) characterized by growth abnormalities and intellectual and developmental delay. Advancements in next-generation sequencing have dramatically increased the detection of pathogenic sequence variants in genes encoding epigenetic machinery associated with NDDs and, concurrently, the number of clinically uninterpretable variants classified as variants of uncertain significance (VUS). Recently, DNA methylation (DNAm) signatures, disorder-specific patterns of DNAm change, have emerged as a functional tool that provides insights into disorder pathophysiology and can classify pathogenicity of variants in NDDs. To date, our group and others have identified DNAm signatures for more than 60 Mendelian neurodevelopmental disorders caused by variants in genes encoding epigenetic machinery. There is broad interest in both the research and clinical communities to develop and catalog DNAm signatures in rare NDDs, but there are challenges in optimizing study design considerations and availability of platforms that integrate bioinformatics tools with the appropriate statistical framework required to analyze genome-wide DNAm data. We previously published EpigenCentral, a platform for analysis of DNAm data in rare NDDs. In this article, we utilize the published Weaver syndrome dataset to provide step-by-step protocols for using EpigenCentral for exploratory analysis to identify DNAm signatures and for classification of NDD variants. We also provide important considerations for experimental design and interpretation of DNAm results. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Exploratory analysis to identify disorder-specific DNAm signatures Basic Protocol 2: Classification of variants associated with neurodevelopmental disorders.

Methylome Profiling in Fabry Disease in Clinical Practice: A Proof of Concept

Anderson−Fabry disease (FD) is an X-linked disease caused by a functional deficit of the α-galactosidase A enzyme. FD diagnosis relies on the clinical manifestations and research of GLA gene mutations. However, because of the lack of a clear genotype/phenotype correlation, FD diagnosis can be challenging. Recently, several studies have highlighted the importance of investigating DNA methylation patterns for confirming the correct diagnosis of different rare Mendelian diseases, but to date, no such studies have been reported for FD. Thus, in the present investigation, we analyzed for the first time the genome-wide methylation profile of a well-characterized cohort of patients with Fabry disease. We profiled the methylation status of about 850,000 CpG sites in 5 FD patients, all carrying the same mutation in the GLA gene (exon 6 c.901C>G) and presenting comparable low levels of α-Gal A activity. We found that, although the whole methylome profile did not discriminate the FD group from the unaffected one, several genes were significantly differentially methylated in Fabry patients. Thus, we provide here a proof of concept, to be tested in patients with different mutations and in a larger cohort, that the methylation state of specific genes can potentially identify Fabry patients and possibly predict organ involvement and disease evolution.

An HNRNPK-specific DNA methylation signature makes sense of missense variants and expands the phenotypic spectrum of Au-Kline syndrome

Au-Kline syndrome (AKS) is a neurodevelopmental disorder associated with multiple malformations and a characteristic facial gestalt. The first individuals ascertained carried de novo loss-of-function (LoF) variants in HNRNPK. Here, we report 32 individuals with AKS (26 previously unpublished), including 13 with de novo missense variants. We propose new clinical diagnostic criteria for AKS that differentiate it from the clinically overlapping Kabuki syndrome and describe a significant phenotypic expansion to include individuals with missense variants who present with subtle facial features and few or no malformations. Many gene-specific DNA methylation (DNAm) signatures have been identified for neurodevelopmental syndromes. Because HNRNPK has roles in chromatin and epigenetic regulation, we hypothesized that pathogenic variants in HNRNPK may be associated with a specific DNAm signature. Here, we report a unique DNAm signature for AKS due to LoF HNRNPK variants, distinct from controls and Kabuki syndrome. This DNAm signature is also identified in some individuals with de novo HNRNPK missense variants, confirming their pathogenicity and the phenotypic expansion of AKS to include more subtle phenotypes. Furthermore, we report that some individuals with missense variants have an "intermediate" DNAm signature that parallels their milder clinical presentation, suggesting the presence of an epi-genotype phenotype correlation. In summary, the AKS DNAm signature may help elucidate the underlying pathophysiology of AKS. This DNAm signature also effectively supported clinical syndrome delineation and is a valuable aid for variant interpretation in individuals where a clinical diagnosis of AKS is unclear, particularly for mild presentations.

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.

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

BACKGROUND

Pathogenic mutations in genes that control chromatin function have been implicated in rare genetic syndromes. These chromatin modifiers exhibit extraordinary diversity in the scale of the epigenetic changes they affect, from single basepair modifications by DNMT1 to whole genome structural changes by PRM1/2. Patterns of DNA methylation are related to a diverse set of epigenetic features across this full range of epigenetic scale, making DNA methylation valuable for mapping regions of general epigenetic dysregulation. However, existing methods are unable to accurately identify regions of differential methylation across this full range of epigenetic scale directly from DNA methylation data.

RESULTS

To address this, we developed DMRscaler, a novel method that uses an iterative windowing procedure to capture regions of differential DNA methylation (DMRs) ranging in size from single basepairs to whole chromosomes. We benchmarked DMRscaler against several DMR callers in simulated and natural data comparing XX and XY peripheral blood samples. DMRscaler was the only method that accurately called DMRs ranging in size from 100 bp to 1 Mb (pearson's r = 0.94) and up to 152 Mb on the X-chromosome. We then analyzed methylation data from rare-disease cohorts that harbor chromatin modifier gene mutations in NSD1, EZH2, and KAT6A where DMRscaler identified novel DMRs spanning gene clusters involved in development.

CONCLUSION

Taken together, our results show DMRscaler is uniquely able to capture the size of DMR features across the full range of epigenetic scale and identify novel, co-regulated regions that drive epigenetic dysregulation in human disease.

X-linked myotubular myopathy is associated with epigenetic alterations and is ameliorated by HDAC inhibition

X-linked myotubular myopathy (XLMTM) is a fatal neuromuscular disorder caused by loss of function mutations in MTM1. At present, there are no directed therapies for XLMTM, and incomplete understanding of disease pathomechanisms. To address these knowledge gaps, we performed a drug screen in mtm1 mutant zebrafish and identified four positive hits, including valproic acid, which functions as a potent suppressor of the mtm1 zebrafish phenotype via HDAC inhibition. We translated these findings to a mouse XLMTM model, and showed that valproic acid ameliorates the murine phenotype. These observations led us to interrogate the epigenome in Mtm1 knockout mice; we found increased DNA methylation, which is normalized with valproic acid, and likely mediated through aberrant 1-carbon metabolism. Finally, we made the unexpected observation that XLMTM patients share a distinct DNA methylation signature, suggesting that epigenetic alteration is a conserved disease feature amenable to therapeutic intervention.

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.

Further Introduction of DNA Methylation (DNAm) Arrays in Regular Diagnostics

Methylation tests have been used for decades in regular DNA diagnostics focusing primarily on Imprinting disorders or specific loci annotated to specific disease associated gene promotors. With the introduction of DNA methylation (DNAm) arrays such as the Illumina Infinium HumanMethylation450 Beadchip array or the Illumina Infinium Methylation EPIC Beadchip array (850 k), it has become feasible to study the epigenome in a timely and cost-effective way. This has led to new insights regarding the complexity of well-studied imprinting disorders such as the Beckwith Wiedemann syndrome, but it has also led to the introduction of tests such as EpiSign, implemented as a diagnostic test in which a single array experiment can be compared to databases with known episignatures of multiple genetic disorders, especially neurodevelopmental disorders. The successful use of such DNAm tests is rapidly expanding. More and more disorders are found to be associated with discrete episignatures which enables fast and definite diagnoses, as we have shown. The first examples of environmentally induced clinical disorders characterized by discrete aberrant DNAm are discussed underlining the broad application of DNAm testing in regular diagnostics. Here we discuss exemplary findings in our laboratory covering this broad range of applications and we discuss further use of DNAm tests in the near future.

Discovery of the First-in-Class G9a/GLP Covalent Inhibitors

The highly homologous protein lysine methyltransferases G9a and GLP, which catalyze mono- and dimethylation of histone H3 lysine 9 (H3K9), have been implicated in various human diseases. To investigate functions of G9a and GLP in human diseases, we and others reported several noncovalent reversible small-molecule inhibitors of G9a and GLP. Here, we report the discovery of the first-in-class G9a/GLP covalent irreversible inhibitors, 1 and 8 (MS8511), by targeting a cysteine residue at the substrate binding site. We characterized these covalent inhibitors in enzymatic, mass spectrometry based and cellular assays and using X-ray crystallography. Compared to the noncovalent G9a/GLP inhibitor UNC0642, covalent inhibitor 8 displayed improved potency in enzymatic and cellular assays. Interestingly, compound 8 also displayed potential kinetic preference for covalently modifying G9a over GLP. Collectively, compound 8 could be a useful chemical tool for studying the functional roles of G9a and GLP by covalently modifying and inhibiting these methyltransferases.

DNA methylation signature associated with Bohring-Opitz syndrome: a new tool for functional classification of variants in ASXL genes

The additional sex combs-like (ASXL) gene family-encoded by ASXL1, ASXL2, and ASXL3-is crucial for mammalian development. Pathogenic variants in the ASXL gene family are associated with three phenotypically distinct neurodevelopmental syndromes. Our previous work has shown that syndromic conditions caused by pathogenic variants in epigenetic regulatory genes show consistent patterns of genome-wide DNA methylation (DNAm) alterations, i.e., DNAm signatures in peripheral blood. Given the role of ASXL1 in chromatin modification, we hypothesized that pathogenic ASXL1 variants underlying Bohring-Opitz syndrome (BOS) have a unique DNAm signature. We profiled whole-blood DNAm for 17 ASXL1 variants, and 35 sex- and age-matched typically developing individuals, using Illumina's Infinium EPIC array. We identified 763 differentially methylated CpG sites in individuals with BOS. Differentially methylated sites overlapped 323 unique genes, including HOXA5 and HOXB4, supporting the functional relevance of DNAm signatures. We used a machine-learning classification model based on the BOS DNAm signature to classify variants of uncertain significance in ASXL1, as well as pathogenic ASXL2 and ASXL3 variants. The DNAm profile of one individual with the ASXL2 variant was BOS-like, whereas the DNAm profiles of three individuals with ASXL3 variants were control-like. We also used Horvath's epigenetic clock, which showed acceleration in DNAm age in individuals with pathogenic ASXL1 variants, and the individual with the pathogenic ASXL2 variant, but not in individuals with ASXL3 variants. These studies enhance our understanding of the epigenetic dysregulation underpinning ASXL gene family-associated syndromes.

SOX11 variants cause a neurodevelopmental disorder with infrequent ocular malformations and hypogonadotropic hypogonadism and with distinct DNA methylation profile

PURPOSE

This study aimed to undertake a multidisciplinary characterization of the phenotype associated with SOX11 variants.

METHODS

Individuals with protein altering variants in SOX11 were identified through exome and genome sequencing and international data sharing. Deep clinical phenotyping was undertaken by referring clinicians. Blood DNA methylation was assessed using Infinium MethylationEPIC array. The expression pattern of SOX11 in developing human brain was defined using RNAscope.

RESULTS

We reported 38 new patients with SOX11 variants. Idiopathic hypogonadotropic hypogonadism was confirmed as a feature of SOX11 syndrome. A distinctive pattern of blood DNA methylation was identified in SOX11 syndrome, separating SOX11 syndrome from other BAFopathies.

CONCLUSION

SOX11 syndrome is a distinct clinical entity with characteristic clinical features and episignature differentiating it from BAFopathies.