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

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.

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.

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 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.

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.

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.

Comparison of methylation episignatures in KMT2B- and KMT2D-related human disorders

Aim & methods: To investigate peripheral blood methylation episignatures in KMT2B-related dystonia (DYT-KMT2B), the authors undertook genome-wide methylation profiling of ∼2 M CpGs using a next-generation sequencing-based assay and compared the findings with those in controls and patients with KMT2D-related Kabuki syndrome type 1 (KS1). Results: A total of 1812 significantly differentially methylated CpG positions (false discovery rate < 0.05) were detected in DYT-KMT2B samples compared with controls. Multi-dimensional scaling analysis showed that the 10 DYT-KMT2B samples clustered together and separately from 29 controls and 10 with pathogenic variants in KMT2D. The authors found that most differentially methylated CpG positions were specific to one disorder and that all (DYT-KMT2B) and most (Kabuki syndrome type 1) methylation alterations in CpG islands were gain of methylation events. Conclusion: Using sensitive methylation profiling methodology, the authors replicated recent reports of a methylation episignature for DYT-KMT2B. These findings will facilitate the development of episignature-based assays to improve diagnostic accuracy.

Delineation of a Phenotype Caused by a KAT6B Missense Variant Not Resembling Say-Barber-Biesecker-Young-Simpson and Genitopatellar Syndromes

Say-Barber-Biesecker-Young-Simpson syndrome (SBBYSS) and genitopatellar syndrome (GPS) are caused by variants of lysine acetyltransferase 6B (KAT6B). These variants tend to occur in the terminal exons of KAT6B. Here, we report a patient with global developmental delay, intellectual disability, autistic behavior, muscular hypotonia, facial dysmorphism, and seizures caused by a novel missense variant in exon 7 of KAT6B. The patient showed a phenotype differing from those of SBBYSS and GPS. We also report patients with missense variants in the proximal exons of KAT6B showing dysmorphic features and autistic behavior not resembling the characteristics of SBBYSS and GPS. Missense variants in the proximal exons of KAT6B may have a dominant negative effect or cause gain of function, leading to unique phenotypes not resembling those of SBBYSS and GPS.

CHD7 regulates bone-fat balance by suppressing PPAR-γ signaling

Chromodomain helicase DNA-binding protein 7 (CHD7), an ATP-dependent eukaryotic chromatin remodeling enzyme, is essential for the development of organs. The mutation of CHD7 is the main cause of CHARGE syndrome, but its function and mechanism in skeletal system remain unclear. Here, we show conditional knockout of Chd7 in bone marrow mesenchymal stem cells (MSCs) and preosteoblasts leads to a pathological phenotype manifested as low bone mass and severely high marrow adiposity. Mechanistically, we identify enhancement of the peroxisome proliferator-activated receptor (PPAR) signaling in Chd7-deficient MSCs. Loss of Chd7 reduces the restriction of PPAR-γ and then PPAR-γ associates with trimethylated histone H3 at lysine 4 (H3K4me3), which subsequently activates the transcription of downstream adipogenic genes and disrupts the balance between osteogenic and adipogenic differentiation. Our data illustrate the pathological manifestations of Chd7 mutation in MSCs and reveal an epigenetic mechanism in skeletal health and diseases.

CHD7 is chromatin remodeler and mutations of CHD7 are the main cause of CHARGE syndrome. Here the authors show that conditional knockout of Chd7 in bone marrow mesenchymal stem cells (MSCs) and pre-osteoblasts leads to a skeletal system development disorder in mice and upregulated PPAR signaling, disrupting the balance between osteogenic and adipogenic differentiation.

Interplay between chromatin marks in development and disease

DNA methylation (DNAme) and histone post-translational modifications (PTMs) have important roles in transcriptional regulation. Although many reports have characterized the functions of such chromatin marks in isolation, recent genome-wide studies reveal surprisingly complex interactions between them. Here, we focus on the interplay between DNAme and methylation of specific lysine residues on the histone H3 tail. We describe the impact of genetic perturbation of the relevant methyltransferases in the mouse on the landscape of chromatin marks as well as the transcriptome. In addition, we discuss the specific neurodevelopmental growth syndromes and cancers resulting from pathogenic mutations in the human orthologues of these genes. Integrating these observations underscores the fundamental importance of crosstalk between DNA and histone H3 methylation in development and disease.

Pseudotime Analysis Reveals Exponential Trends in DNA Methylation Aging with Mortality Associated Timescales

The epigenetic trajectory of DNA methylation profiles has a nonlinear relationship with time, reflecting rapid changes in DNA methylation early in life that progressively slow with age. In this study, we use pseudotime analysis to determine the functional form of these trajectories. Unlike epigenetic clocks that constrain the functional form of methylation changes with time, pseudotime analysis orders samples along a path, based on similarities in a latent dimension, to provide an unbiased trajectory. We show that pseudotime analysis can be applied to DNA methylation in human blood and brain tissue and find that it is highly correlated with the epigenetic states described by the Epigenetic Pacemaker. Moreover, we show that the pseudotime trajectory can be modeled with respect to time, using a sum of two exponentials, with coefficients that are close to the timescales of human age-associated mortality. Thus, for the first time, we can identify age-associated molecular changes that appear to track the exponential dynamics of mortality risk.

Accelerated epigenetic age and shortened telomere length based on DNA methylation in Nicolaides-Baraitser syndrome

BACKGROUND

Nicolaides-Baraitser syndrome (NCBRS) is a rare disorder characterized by neurodevelopmental delays, seizures, and diverse physical characteristics. The DNA methylation (DNAm) profile in NCBRS is significantly different. DNAm is linked to the biological aging of cells and the health risks associated with biological aging. In this study, we examined changes in biological ages in NCBRS to provide insights into the prognosis and health risks of NCBRS.

METHODS

We used a publicly available dataset to examine biological aging in NCBRS using DNAm-based epigenetic ages, such as PhenoAge and GrimAge, as well as DNAm-based estimator of telomere length (DNAmTL). We investigated 12 cases, clinically diagnosed as NCBRS, and 27 controls.

RESULTS

Compared to controls, NCBRS cases exhibited significantly accelerated PhenoAge and GrimAge as well as significantly shortened DNAmTL.

CONCLUSION

These results suggest an acceleration of biological aging in NCBRS and provide insights into the prognosis and health risks of NCBRS.

The utility of DNA methylation signatures in directing genome sequencing workflow: Kabuki syndrome and CDK13-related disorder

Kabuki syndrome (KS) is a neurodevelopmental disorder characterized by hypotonia, intellectual disability, skeletal anomalies, and postnatal growth restriction. The characteristic facial appearance is not pathognomonic for KS as several other conditions demonstrate overlapping features. For 20‐30% of children with a clinical diagnosis of KS, no causal variant is identified by conventional genetic testing of the two associated genes, KMT2D and KDM6A. Here, we describe two cases of suspected KS that met clinical diagnostic criteria and had a high gestalt match on the artificial intelligence platform Face2Gene. Although initial KS testing was negative, genome‐wide DNA methylation (DNAm) was instrumental in guiding genome sequencing workflow to establish definitive molecular diagnoses. In one case, a positive DNAm signature for KMT2D led to the identification of a cryptic variant in KDM6A by genome sequencing; for the other case, a DNAm signature different from KS led to the detection of another diagnosis in the KS differential, CDK13‐related disorder. This approach illustrates the clinical utility of DNAm signatures in the diagnostic workflow for the genome analyst or clinical geneticist—especially for disorders with overlapping clinical phenotypes.

Novel diagnostic DNA methylation episignatures expand and refine the epigenetic landscapes of Mendelian disorders

Overlapping clinical phenotypes and an expanding breadth and complexity of genomic associations are a growing challenge in the diagnosis and clinical management of Mendelian disorders. The functional consequences and clinical impacts of genomic variation may involve unique, disorder-specific, genomic DNA methylation episignatures. In this study, we describe 19 novel episignature disorders and compare the findings alongside 38 previously established episignatures for a total of 57 episignatures associated with 65 genetic syndromes. We demonstrate increasing resolution and specificity ranging from protein complex, gene, sub-gene, protein domain, and even single nucleotide-level Mendelian episignatures. We show the power of multiclass modeling to develop highly accurate and disease-specific diagnostic classifiers. This study significantly expands the number and spectrum of disorders with detectable DNA methylation episignatures, improves the clinical diagnostic capabilities through the resolution of unsolved cases and the reclassification of variants of unknown clinical significance, and provides further insight into the molecular etiology of Mendelian conditions.

OUP accepted manuscript

BACKGROUND

TASP1 encodes an endopeptidase activating histone methyltransferases of the KMT2 family. Homozygous loss-of-function variants in TASP1 have recently been associated with Suleiman-El-Hattab syndrome. We report six individuals with Suleiman-El-Hattab syndrome and provide functional characterization of this novel histone modification disorder in a multi-omics approach.

METHODS

Chromosomal microarray/exome sequencing in all individuals. Western blotting from fibroblasts in two individuals. RNA sequencing and proteomics from fibroblasts in one individual. Methylome analysis from blood in two individuals. Knock-out of tasp1 orthologue in zebrafish and phenotyping.

RESULTS

All individuals had biallelic TASP1 loss-of-function variants and a phenotype including developmental delay, multiple congenital anomalies (including cardiovascular and posterior fossa malformations), a distinct facial appearance and happy demeanor. Western blot revealed absence of TASP1. RNA sequencing/proteomics showed HOX gene downregulation (HOXA4, HOXA7, HOXA1 and HOXB2) and dysregulation of transcription factor TFIIA. A distinct methylation profile intermediate between control and Kabuki syndrome (KMT2D) profiles could be produced. Zebrafish tasp1 knock-out revealed smaller head size and abnormal cranial cartilage formation in tasp1 crispants.

CONCLUSION

This work further delineates Suleiman-El-Hattab syndrome, a recognizable neurodevelopmental syndrome. Possible downstream mechanisms of TASP1 deficiency include perturbed HOX gene expression and dysregulated TFIIA complex. Methylation pattern suggests that Suleiman-El-Hattab syndrome can be categorized into the group of histone modification disorders including Wiedemann-Steiner and Kabuki syndrome.

Deficiency of TET3 leads to a genome-wide DNA hypermethylation episignature in human whole blood

TET3 encodes an essential dioxygenase involved in epigenetic regulation through DNA demethylation. TET3 deficiency, or Beck-Fahrner syndrome (BEFAHRS; MIM: 618798), is a recently described neurodevelopmental disorder of the DNA demethylation machinery with a nonspecific phenotype resembling other chromatin-modifying disorders, but inconsistent variant types and inheritance patterns pose diagnostic challenges. Given TET3's direct role in regulating 5-methylcytosine and recent identification of syndrome-specific DNA methylation profiles, we analyzed genome-wide DNA methylation in whole blood of TET3-deficient individuals and identified an episignature that distinguishes affected and unaffected individuals and those with mono-allelic and bi-allelic pathogenic variants. Validation and testing of the episignature correctly categorized known TET3 variants and determined pathogenicity of variants of uncertain significance. Clinical utility was demonstrated when the episignature alone identified an affected individual from over 1000 undiagnosed cases and was confirmed upon distinguishing TET3-deficient individuals from those with 46 other disorders. The TET3-deficient signature - and the signature resulting from activating mutations in DNMT1 which normally opposes TET3 - are characterized by hypermethylation, which for BEFAHRS involves CpG sites that may be biologically relevant. This work expands the role of epi-phenotyping in molecular diagnosis and reveals genome-wide DNA methylation profiling as a quantitative, functional readout for characterization of this new biochemical category of disease.

Genome sequencing identifies three molecular diagnoses including a mosaic variant in the COL2A1 gene in an individual with Pol III-related leukodystrophy and Feingold syndrome

Undiagnosed genetic disease imposes a significant burden on families and health-care resources, especially in cases with a complex phenotype. Here we present a child with suspected leukodystrophy in the context of additional features, including hearing loss, clinodactyly, rotated thumbs, tapered fingers, and simplified palmar crease. Trio genome sequencing (GS) identified three molecular diagnoses in this individual: compound heterozygous missense variants associated with polymerase III (Pol III)–related leukodystrophy, a 4-Mb de novo copy-number loss including the MYCN gene associated with Feingold syndrome, and a mosaic single-nucleotide variant associated with COL2A1-related disorders. These variants fully account for the individual's features, but also illustrate the potential for superimposed and unclear contributions of multiple diagnoses to an individual's overall presentation. This report demonstrates the advantage of GS in detection of multiple variant types, including low-level mosaic variants, and emphasizes the need for comprehensive genetic analysis and detailed clinical phenotyping to provide individuals and their families with the maximum benefit for clinical care and genetic counseling.