Evaluation of DNA Methylation Episignatures for Diagnosis and Phenotype Correlations in 42 Mendelian Neurodevelopmental Disorders

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

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

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

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

Comprehensive EHMT1 variants analysis broadens genotype-phenotype associations and molecular mechanisms in Kleefstra syndrome

The shift to a genotype-first approach in genetic diagnostics has revolutionized our understanding of neurodevelopmental disorders, expanding both their molecular and phenotypic spectra. Kleefstra syndrome (KLEFS1) is caused by EHMT1 haploinsufficiency and exhibits broad clinical manifestations. EHMT1 encodes euchromatic histone methyltransferase-1-a pivotal component of the epigenetic machinery. We have recruited 209 individuals with a rare EHMT1 variant and performed comprehensive molecular in silico and in vitro testing alongside DNA methylation (DNAm) signature analysis for the identified variants. We (re)classified the variants as likely pathogenic/pathogenic (molecularly confirming Kleefstra syndrome) in 191 individuals. We provide an updated and broader clinical and molecular spectrum of Kleefstra syndrome, including individuals with normal intelligence and familial occurrence. Analysis of the EHMT1 variants reveals a broad range of molecular effects and their associated phenotypes, including distinct genotype-phenotype associations. Notably, we showed that disruption of the "reader" function of the ankyrin repeat domain by a protein altering variant (PAV) results in a KLEFS1-specific DNAm signature and milder phenotype, while disruption of only "writer" methyltransferase activity of the SET domain does not result in KLEFS1 DNAm signature or typical KLEFS1 phenotype. Similarly, N-terminal truncating variants result in a mild phenotype without the DNAm signature. We demonstrate how comprehensive variant analysis can provide insights into pathogenesis of the disorder and DNAm signature. In summary, this study presents a comprehensive overview of KLEFS1 and EHMT1, revealing its broader spectrum and deepening our understanding of its molecular mechanisms, thereby informing accurate variant interpretation, counseling, and clinical management.

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

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

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

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

Genome-wide epigenetic signatures facilitated the variant classification of the PURA gene and uncovered the pathomechanism of PURA-related neurodevelopmental disorders

PURPOSE

Rare genetic variants in the PURA gene cause the PURA-related neurodevelopmental disorder (PURA-NDD), characterized by neonatal abnormalities and developmental delay. Using genome-wide DNA methylation analysis on patients with PURA variants, we aim to establish a PURA-NDD-specific methylation profile and provide further insights on the molecular basis of the PURA-NDD.

METHODS

Twenty three individuals (including 12 unpublished) carrying PURA variants were enrolled. We conducted the Illumina Infinium EPIC microarray analysis in 17 PURA-NDD individuals. In vitro experiments were performed to examine how PURA variants affect Pur-a expression.

RESULTS

Additional phenotypes in 12 newly identified patients were described in this study. Genome-wide DNA methylation analysis unveiled distinctive methylation profiles to PURA-NDD, and the established classifier can reclassify PURA variants of uncertain significance. Patients bearing PURA hapoloinsufficient and missense variants have comparable DNA methylation profiles, and cells expressing these PURA variants showed consistent Pur-a downregulation, suggesting a haploinsufficiency mechanism.

CONCLUSION

Patients with PURA-NDD exhibit a specific episignature, which has potential to aid identification and diagnosis of PURA-NDD patients and offer implications for further functional investigations.

Menke-Hennekam syndrome; delineation of domain-specific subtypes with distinct clinical and DNA methylation profiles

CREB-binding protein (CBP, encoded by CREBBP) and its paralog E1A-associated protein (p300, encoded by EP300) are involved in histone acetylation and transcriptional regulation. Variants that produce a null allele or disrupt the catalytic domain of either protein cause Rubinstein-Taybi syndrome (RSTS), while pathogenic missense and in-frame indel variants in parts of exons 30 and 31 cause phenotypes recently described as Menke-Hennekam syndrome (MKHK). To distinguish MKHK subtypes and define their characteristics, molecular and extended clinical data on 82 individuals (54 unpublished) with variants affecting CBP (n = 71) or p300 (n = 11) (NP_004371.2 residues 1,705-1,875 and NP_001420.2 residues 1,668-1,833, respectively) were summarized. Additionally, genome-wide DNA methylation profiles were assessed in DNA extracted from whole peripheral blood from 54 individuals. Most variants clustered closely around the zinc-binding residues of two zinc-finger domains (ZZ and TAZ2) and within the first α helix of the fourth intrinsically disordered linker (ID4) of CBP/p300. Domain-specific methylation profiles were discerned for the ZZ domain in CBP/p300 (found in nine out of 10 tested individuals) and TAZ2 domain in CBP (in 14 out of 20), while a domain-specific diagnostic episignature was refined for the ID4 domain in CBP/p300 (in 21 out of 21). Phenotypes including intellectual disability of varying degree and distinct physical features were defined for each of the regions. These findings demonstrate existence of at least three MKHK subtypes, which are domain specific (MKHK-ZZ, MKHK-TAZ2, and MKHK-ID4) rather than gene specific (CREBBP/EP300). DNA methylation episignatures enable stratification of molecular pathophysiologic entities within a gene or across a family of paralogous genes.

DNA methylation analysis in patients with neurodevelopmental disorders improves variant interpretation and reveals complexity

Analysis of genomic DNA methylation by generating epigenetic signature profiles (episignatures) is increasingly being implemented in genetic diagnosis. Here we report our experience using episignature analysis to resolve both uncomplicated and complex cases of neurodevelopmental disorders (NDDs). We analyzed 97 NDDs divided into (1) a validation cohort of 59 patients with likely pathogenic/pathogenic variants characterized by a known episignature and (2) a test cohort of 38 patients harboring variants of unknown significance or unidentified variants. The expected episignature was obtained in most cases with likely pathogenic/pathogenic variants (53/59 [90%]), a revealing exception being the overlapping profile of two SMARCB1 pathogenic variants with ARID1A/B:c.6200, confirmed by the overlapping clinical features. In the test cohort, five cases showed the expected episignature, including (1) novel pathogenic variants in ARID1B and BRWD3; (2) a deletion in ATRX causing MRXFH1 X-linked mental retardation; and (3) confirmed the clinical diagnosis of Cornelia de Lange (CdL) syndrome in mutation-negative CdL patients. Episignatures analysis of the in BAF complex components revealed novel functional protein interactions and common episignatures affecting homologous residues in highly conserved paralogous proteins (SMARCA2 M856V and SMARCA4 M866V). Finally, we also found sex-dependent episignatures in X-linked disorders. Implementation of episignature profiling is still in its early days, but with increasing utilization comes increasing awareness of the capacity of this methodology to help resolve the complex challenges of genetic diagnoses.

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

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

Genetic susceptibility to neurodevelopmental conditions associates with neonatal DNA methylation patterns in the general population: an individual participant data meta-analysis

Background

Autism spectrum disorder (ASD), attention-deficit/hyperactivity disorder (ADHD), and schizophrenia (SCZ) are highly heritable and linked to disruptions in foetal (neuro)development. While epigenetic processes are considered an important underlying pathway between genetic susceptibility and neurodevelopmental conditions, it is unclear (i) whether genetic susceptibility to these conditions is associated with epigenetic patterns, specifically DNA methylation (DNAm), already at birth; (ii) to what extent DNAm patterns are unique or shared across conditions, and (iii) whether these neonatal DNAm patterns can be leveraged to enhance genetic prediction of (neuro)developmental outcomes.

Methods

We conducted epigenome-wide meta-analyses of genetic susceptibility to ASD, ADHD, and schizophrenia, quantified using polygenic scores (PGSs) on cord blood DNAm, using four population-based cohorts (n pooled=5,802), all North European. Heterogeneity statistics were used to estimate overlap in DNAm patterns between PGSs. Subsequently, DNAm-based measures of PGSs were built in a target sample, and used as predictors to test incremental variance explained over PGS in 130 (neuro)developmental outcomes spanning birth to 14 years.

Outcomes

In probe-level analyses, SCZ-PGS associated with neonatal DNAm at 246 loci (p<9×10-8), predominantly in the major histocompatibility complex. Functional characterization of these DNAm loci confirmed strong genetic effects, significant blood-brain concordance and enrichment for immune-related pathways. 8 loci were identified for ASD-PGS (mapping to FDFT1 and MFHAS1), and none for ADHD-PGS. Regional analyses indicated a large number of differentially methylated regions for all PGSs (SCZ-PGS: 157, ASD-PGS: 130, ADHD-PGS: 166). DNAm signals showed little overlap between PGSs. We found suggestive evidence that incorporating DNAm-based measures of genetic susceptibility at birth increases explained variance for several child cognitive and motor outcomes over and above PGS.

Interpretation

Genetic susceptibility for neurodevelopmental conditions, particularly schizophrenia, is detectable in cord blood DNAm at birth in a population-based sample, with largely distinct DNAm patterns between PGSs. These findings support an early-origins perspective on schizophrenia.

Funding

HorizonEurope; European Research Council.

mLiftOver: harmonizing data across Infinium DNA methylation platforms

MOTIVATION

Infinium DNA methylation BeadChips are widely used for genome-wide DNA methylation profiling at the population scale. Recent updates to probe content and naming conventions in the EPIC version 2 (EPICv2) arrays have complicated integrating new data with previous Infinium array platforms, such as the MethylationEPIC (EPIC) and the HumanMethylation450 (HM450) BeadChip.

RESULTS

We present mLiftOver, a user-friendly tool that harmonizes probe ID, methylation level, and signal intensity data across different Infinium platforms. It manages probe replicates, missing data imputation, and platform-specific bias for accurate data conversion. We validated the tool by applying HM450-based cancer classifiers to EPICv2 cancer data, achieving high accuracy. Additionally, we successfully integrated EPICv2 healthy tissue data with legacy HM450 data for tissue identity analysis and produced consistent copy number profiles in cancer cells.

AVAILABILITY AND IMPLEMENTATION

mLiftOver is implemented R and available in the Bioconductor package SeSAMe (version 1.21.13+): https://bioconductor.org/packages/release/bioc/html/sesame.html. Analysis of EPIC and EPICv2 platform-specific bias and high-confidence mapping is available at https://github.com/zhou-lab/InfiniumAnnotationV1/raw/main/Anno/EPICv2/EPICv2ToEPIC_conversion.tsv.gz. The source code is available at https://github.com/zwdzwd/sesame/blob/devel/R/mLiftOver.R under the MIT license.

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.

Blepharophimosis with intellectual disability and Helsmoortel-Van Der Aa Syndrome share episignature and phenotype

Blepharophimosis with intellectual disability (BIS) is a recently recognized disorder distinct from Nicolaides-Baraister syndrome that presents with distinct facial features of blepharophimosis, developmental delay, and intellectual disability. BIS is caused by pathogenic variants in SMARCA2, that encodes the catalytic subunit of the superfamily II helicase group of the BRG1 and BRM-associated factors (BAF) forming the BAF complex, a chromatin remodeling complex involved in transcriptional regulation. Individuals bearing variants within the bipartite nuclear localization (BNL) signal domain of ADNP present with the neurodevelopmental disorder known as Helsmoortel-Van Der Aa Syndrome (HVDAS). Distinct DNA methylation profiles referred to as episignatures have been reported in HVDAS and BAF complex disorders. Due to molecular interactions between ADNP and BAF complex, and an overlapping craniofacial phenotype with narrowing of the palpebral fissures in a subset of patients with HVDAS and BIS, we hypothesized the possibility of a common phenotype-specific episignature. A distinct episignature was shared by 15 individuals with BIS-causing SMARCA2 pathogenic variants and 12 individuals with class II HVDAS caused by truncating pathogenic ADNP variants. This represents first evidence of a sensitive phenotype-specific episignature biomarker shared across distinct genetic conditions that also exhibit unique gene-specific episignatures.

Refining the 9q34.3 microduplication syndrome reveals mild neurodevelopmental features associated with a distinct global DNA methylation profile

Precise regulation of gene expression is important for correct neurodevelopment. 9q34.3 deletions affecting the EHMT1 gene result in a syndromic neurodevelopmental disorder named Kleefstra syndrome. In contrast, duplications of the 9q34.3 locus encompassing EHMT1 have been suggested to cause developmental disorders, but only limited information has been available. We have identified 15 individuals from 10 unrelated families, with 9q34.3 duplications <1.5 Mb in size, encompassing EHMT1 entirely. Clinical features included mild developmental delay, mild intellectual disability or learning problems, autism spectrum disorder, and behavior problems. The individuals did not consistently display dysmorphic features, congenital anomalies, or growth abnormalities. DNA methylation analysis revealed a weak DNAm profile for the cases with 9q34.3 duplication encompassing EHMT1, which could segregate the majority of the affected cases from controls. This study shows that individuals with 9q34.3 duplications including EHMT1 gene present with mild non-syndromic neurodevelopmental disorders and DNA methylation changes different from Kleefstra syndrome.

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.

ARID2, a milder cause of Coffin-Siris Syndrome? Broadening the phenotype with 17 additional individuals

Coffin-Siris Syndrome (CSS, MIM 135900) is now a well-described genetic condition caused by pathogenic variants in the Bromocriptine activating factor (BAF) complex, including ARID1B, ARID1A, ARID2, SMARCA4, SMARCE1, SMARCB1, SOX11, SMARCC2, DPF2, and more recently, BICRA. Individuals with CSS have a spectrum of various medical challenges, most often evident at birth, including feeding difficulties, hypotonia, organ-system anomalies, and learning and developmental differences. The classic finding of fifth digit hypo- or aplasia is seen variably. ARID2, previously described, is one of the less frequently observed gene changes in CSS. Although individuals with ARID2 have been reported to have classic features of CSS including hypertrichosis, coarse facial features, short stature, and fifth digit anomalies, as with many of the other CSS genes, there appears to be a spectrum of phenotypes. We report here a cohort of 17 individuals with ARID2 variants from the Coffin-Siris/BAF clinical registry and detail their medical challenges as well as developmental progress. Feeding difficulties, hypotonia, and short stature occur often, and hip dysplasia appears to occur more often than with other genes, however more severe medical challenges such as significant brain and cardiac malformations are rarer. Individuals appear to have mild to moderate intellectual impairment and may carry additional diagnoses such as ADHD. Further phenotypic description of this gene will aid clinicians caring for individuals with this rarer form of CSS.

The Phenotype-Based Approach Can Solve Cold Cases: The Paradigm of Mosaic Mutations of the CREBBP Gene

Rubinstein-Taybi syndrome (RTS) is a rare genetic disorder characterized by intellectual disability, facial dysmorphisms, and enlarged thumbs and halluces. Approximately 55% of RTS cases result from pathogenic variants in the CREBBP gene, with an additional 8% linked to the EP300 gene. Given the close relationship between these two genes and their involvement in epigenomic modulation, RTS is grouped into chromatinopathies. The extensive clinical heterogeneity observed in RTS, coupled with the growing number of disorders involving the epigenetic machinery, poses a challenge to a phenotype-based diagnostic approach for these conditions. Here, we describe the first case of a patient clinically diagnosed with RTS with a CREBBP truncating variant in mosaic form. We also review previously described cases of mosaicism in CREBBP and apply clinical diagnostic guidelines to these patients, confirming the good specificity of the consensus. Nonetheless, these reports raise questions about the potential underdiagnosis of milder cases of RTS. The application of a targeted phenotype-based approach, coupled with high-depth NGS, may enhance the diagnostic yield of whole-exome sequencing (WES) in mild and mosaic conditions.

Diagnosis and management in Rubinstein-Taybi syndrome: first international consensus statement

Rubinstein-Taybi syndrome (RTS) is an archetypical genetic syndrome that is characterised by intellectual disability, well-defined facial features, distal limb anomalies and atypical growth, among numerous other signs and symptoms. It is caused by variants in either of two genes (CREBBP, EP300) which encode for the proteins CBP and p300, which both have a function in transcription regulation and histone acetylation. As a group of international experts and national support groups dedicated to the syndrome, we realised that marked heterogeneity currently exists in clinical and molecular diagnostic approaches and care practices in various parts of the world. Here, we outline a series of recommendations that document the consensus of a group of international experts on clinical diagnostic criteria for types of RTS (RTS1: CREBBP; RTS2: EP300), molecular investigations, long-term management of various particular physical and behavioural issues and care planning. The recommendations as presented here will need to be evaluated for improvements to allow for continued optimisation of diagnostics and care.

MSA: scalable DNA methylation screening BeadChip for high-throughput trait association studies

The Infinium DNA Methylation BeadChips have significantly contributed to population-scale epigenetics research by enabling epigenome-wide trait association discoveries. Here, we design, describe, and experimentally verify a new iteration of this technology, the Methylation Screening Array (MSA), to focus on human trait screening and discovery. This array utilizes extensive data from previous Infinium platform-based epigenome-wide association studies (EWAS). It incorporates knowledge from the latest single-cell and cell type-resolution whole genome methylome profiles. The MSA is engineered to achieve scalable screening of epigenetics-trait association in an ultra-high sample throughput. Our design encompassed diverse human trait associations, including those with genetic, cellular, environmental, and demographical variables and human diseases such as genetic, neurodegenerative, cardiovascular, infectious, and immune diseases. We comprehensively evaluated this array's reproducibility, accuracy, and capacity for cell-type deconvolution and supporting 5-hydroxymethylation profiling in diverse human tissues. Our first atlas data using this platform uncovered the complex chromatin and tissue contexts of DNA modification variations and genetic variants linked to human phenotypes.

Narrowing the diagnostic gap: Genomes, episignatures, long-read sequencing, and health economic analyses in an exome-negative intellectual disability cohort

PURPOSE

Genome sequencing (GS)-specific diagnostic rates in prospective tightly ascertained exome sequencing (ES)-negative intellectual disability (ID) cohorts have not been reported extensively.

METHODS

ES, GS, epigenetic signatures, and long-read sequencing diagnoses were assessed in 74 trios with at least moderate ID.

RESULTS

The ES diagnostic yield was 42 of 74 (57%). GS diagnoses were made in 9 of 32 (28%) ES-unresolved families. Repeated ES with a contemporary pipeline on the GS-diagnosed families identified 8 of 9 single-nucleotide variations/copy-number variations undetected in older ES, confirming a GS-unique diagnostic rate of 1 in 32 (3%). Episignatures contributed diagnostic information in 9% with GS corroboration in 1 of 32 (3%) and diagnostic clues in 2 of 32 (6%). A genetic etiology for ID was detected in 51 of 74 (69%) families. Twelve candidate disease genes were identified. Contemporary ES followed by GS cost US$4976 (95% CI: $3704; $6969) per diagnosis and first-line GS at a cost of $7062 (95% CI: $6210; $8475) per diagnosis.

CONCLUSION

Performing GS only in ID trios would be cost equivalent to ES if GS were available at $2435, about a 60% reduction from current prices. This study demonstrates that first-line GS achieves higher diagnostic rate than contemporary ES but at a higher cost.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Autoimmune cytopenia and Kabuki syndrome in paediatrics: Insights in 11 patients

Kabuki syndrome (KS) is now listed in the Human Inborn Errors of Immunity (IEI) Classification. It is a rare disease caused by KMT2D and KDM6A variants, dominated by intellectual disability and characteristic facial features. Recurrently, pathogenic variants are identified in those genes in patients examined for autoimmune cytopenia (AIC), but interpretation remains challenging. This study aims to describe the genetic diagnosis and the clinical management of patients with paediatric-onset AIC and KS. Among 11 patients with AIC and KS, all had chronic immune thrombocytopenic purpura, and seven had Evans syndrome. All had other associated immunopathological manifestations, mainly symptomatic hypogammaglobinaemia. They had a median of 8 (5-10) KS-associated manifestations. Pathogenic variants were detected in KMT2D gene without clustering, during the immunological work-up of AIC in three cases, and the clinical strategy to validate them is emphasized. Eight patients received second-line treatments, mainly rituximab and mycophenolate mofetil. With a median follow-up of 17 (2-31) years, 8/10 alive patients still needed treatment for AIC. First-line paediatricians should be able to recognize and confirm KS in children with ITP or multiple AIC, to provide early appropriate clinical management and specific long-term follow-up. The epigenetic immune dysregulation in KS opens exciting new perspectives.

Low-input and single-cell methods for Infinium DNA methylation BeadChips

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection P-value calculation achieved higher sensitivities for low-input datasets and was validated in over 100 000 public diverse methylome profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

Heterozygous loss-of-function SMC3 variants are associated with variable growth and developmental features

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 14 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism, reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated an overall milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, and some had alternative symptomatologies with rational biological links to SMC3. Analyses of tumor and model system transcriptomic data and epigenetic data in a subset of cases suggest that SMC3 pLoF variants reduce SMC3 expression but do not strongly support clustering with functional genomic signatures of typical CdLS. Our finding of substantial population-scale LoF intolerance in concert with variable growth and developmental features in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multilayered genomic data paired with careful phenotyping.

Haploinsufficiency of ZFHX3, encoding a key player in neuronal development, causes syndromic intellectual disability

Neurodevelopmental disorders (NDDs) result from impaired development and functioning of the brain. Here, we identify loss-of-function (LoF) variation in ZFHX3 as a cause for syndromic intellectual disability (ID). ZFHX3 is a zinc-finger homeodomain transcription factor involved in various biological processes, including cell differentiation and tumorigenesis. We describe 42 individuals with protein-truncating variants (PTVs) or (partial) deletions of ZFHX3, exhibiting variable intellectual disability and autism spectrum disorder, recurrent facial features, relative short stature, brachydactyly, and, rarely, cleft palate. ZFHX3 LoF associates with a specific methylation profile in whole blood extracted DNA. Nuclear abundance of ZFHX3 increases during human brain development and neuronal differentiation. ZFHX3 was found to interact with the chromatin remodeling BRG1/Brm-associated factor complex and the cleavage and polyadenylation complex, suggesting a function in chromatin remodeling and mRNA processing. Furthermore, ChIP-seq for ZFHX3 revealed that it predominantly binds promoters of genes involved in nervous system development. We conclude that loss-of-function variants in ZFHX3 are a cause of syndromic ID associating with a specific DNA methylation profile.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

Data science using the human epigenome for predicting multifactorial diseases and symptoms

Tweetable abstract This article reviews machine learning models that leverages epigenomic data for predicting multifactorial diseases and symptoms as well as how such models can be utilized to explore new research questions.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Mitotic Block and Epigenetic Repression Underlie Neurodevelopmental Defects and Neurobehavioral Deficits in Congenital Heart Disease

Poor neurodevelopment is often observed with congenital heart disease (CHD), especially with mutations in chromatin modifiers. Here analysis of mice with hypoplastic left heart syndrome (HLHS) arising from mutations in Sin3A associated chromatin modifier Sap130 , and adhesion protein Pcdha9, revealed neurodevelopmental and neurobehavioral deficits reminiscent of those in HLHS patients. Microcephaly was associated with impaired cortical neurogenesis, mitotic block, and increased apoptosis. Transcriptional profiling indicated dysregulated neurogenesis by REST, altered CREB signaling regulating memory and synaptic plasticity, and impaired neurovascular coupling modulating cerebral blood flow. Many neurodevelopmental/neurobehavioral disease pathways were recovered, including autism and cognitive impairment. These same pathways emerged from genome-wide DNA methylation and Sap130 chromatin immunoprecipitation sequencing analyses, suggesting epigenetic perturbation. Mice with Pcdha9 mutation or forebrain-specific Sap130 deletion without CHD showed learning/memory deficits and autism-like behavior. These novel findings provide mechanistic insights indicating the adverse neurodevelopment in HLHS may involve cell autonomous/nonautonomous defects and epigenetic dysregulation and suggest new avenues for therapy.

Diagnostic utility of DNA methylation episignature analysis for early diagnosis of KMT2B-related disorders: case report

The lysine methyltransferase 2B (KMT2B) gene product is important for epigenetic modifications associated with active gene transcription in normal development and in maintaining proper neural function. Pathogenic variants in KMT2B have been associated with childhood-onset Dystonia-28 and Intellectual developmental disorder, autosomal dominant 68 (MRD 68) for cases of neurodevelopmental impairment without dystonia (DYT28; OMIM 617284 and MRD68; OMIM 619934, respectively). Since its first description in 2016, approximately one hundred KMT2B genetic variants have been reported with heterogeneous phenotypes, including atypical patterns of dystonia evolution and non-dystonic neurodevelopmental phenotypes. KMT2B-related disorders share many overlapping phenotypic characteristics with other neurodevelopmental disorders and delayed dystonia, that can appear later in childhood, often delaying clinical diagnosis. Furthermore, conventional genetic testing may not always provide actionable information (e.g., gene panel selection based on early clinical presentation or variants of uncertain significance), which prevents patients and families from obtaining early access to treatments and support. Herein, we describe the early diagnosis of KMT2B-related neurodevelopmental disorder by DNA methylation episignature testing in a 4-year-old patient without features of dystonia at diagnosis, which is reported to develop in more than 80% of KMT2B-related disorder cases. The proband, a 4-year-old female of Jewish-Israeli descent, presented with speech delay, microcephaly, poor weight gain, attention-deficit and hyperactivity disorder, dysmorphism, intellectual disabilities and joint hyperlaxity, but presented no signs of dystonia at initial evaluation. Episignature screening in this pre-symptomatic patient enabled accurate genetic diagnosis and timely and actionable intervention earlier in the natural history of Childhood-onset Dystonia-28.

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.

Next-generation sequencing and bioinformatics in rare movement disorders

The ability to sequence entire exomes and genomes has revolutionized molecular testing in rare movement disorders, and genomic sequencing is becoming an integral part of routine diagnostic workflows for these heterogeneous conditions. However, interpretation of the extensive genomic variant information that is being generated presents substantial challenges. In this Perspective, we outline multidimensional strategies for genetic diagnosis in patients with rare movement disorders. We examine bioinformatics tools and computational metrics that have been developed to facilitate accurate prioritization of disease-causing variants. Additionally, we highlight community-driven data-sharing and case-matchmaking platforms, which are designed to foster the discovery of new genotype-phenotype relationships. Finally, we consider how multiomic data integration might optimize diagnostic success by combining genomic, epigenetic, transcriptomic and/or proteomic profiling to enable a more holistic evaluation of variant effects. Together, the approaches that we discuss offer pathways to the improved understanding of the genetic basis of rare movement disorders.

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.

Clinical heterogeneity of polish patients with KAT6B-related disorder

BACKGROUND

Say-Barber-Biesecker-Young-Simpson (SBBYSS) variant of Ohdo syndrome is a rare, autosomal dominant and clinically heterogenous disorder, caused by pathogenic variants in the KAT6B gene located on chromosome 10q22.2. KAT6B encodes a highly conserved histone acetyltransferase belonging to the MYST family. Currently, diseases caused by pathogenic variants in KAT6B (KAT6B-related disorders) comprise two allelic entities: SBBYSS variant of Ohdo syndrome and genitopatellar syndrome (GPS). Increase in the number of cases with overlapping GPS/SBBYSS phenotype which makes it necessary to redefine this group of phenotypes as KAT6B-related disorders or KAT6B spectrum disorders. Individuals with SBBYSS usually present with facial abnormalities, hypotonia, joint laxity, feeding problems, and long thumbs/great toes. This syndrome also typically involves skeletal problems including patellar hypoplasia/agenesis.

METHODS

Here we report six SBBYS syndrome patients with the same dysmorphic features but a different course of the disease. One known and five novel KATB6 pathogenic variants were identified by molecular diagnostics using Next Generation Sequencing (NGS).

RESULTS

We present a detailed phenotypic analysis of six individuals with KAT6B-related disorders, in whom a heterozygous pathogenic variant in KAT6B gene was found. In all of our patients facial dysmorphism as well as developmental and speech delay were present. Additionally, all but one patients presented with hypotonia, ocular abnormalities and long thumbs. Most of our probands showed blepharophimosis and skeletal (mainly knee) defects. Contrary to previously reported severe patellar defects (hypoplasia/agenesis) anomalies presented by our patients were less severe (dysplasia, habitual dislocation, subluxation) referring to KAT6B-related disorders.

CONCLUSION

While most of the anomalies found in our patients comply with SBBYSS criteria, phenotypic differences in our probands support a broader spectrum of the disease phenotype. To establish the range of this spectrum, a detailed analysis of clinical variability among patients with SBBYSS requires further investigation.

Identification of a robust DNA methylation signature for Fanconi anemia

Fanconi anemia (FA) is a clinically variable and genetically heterogeneous cancer-predisposing disorder representing the most common bone marrow failure syndrome. It is caused by inactivating predominantly biallelic mutations involving >20 genes encoding proteins with roles in the FA/BRCA DNA repair pathway. Molecular diagnosis of FA is challenging due to the wide spectrum of the contributing gene mutations and structural rearrangements. The assessment of chromosomal fragility after exposure to DNA cross-linking agents is generally required to definitively confirm diagnosis. We assessed peripheral blood genome-wide DNA methylation (DNAm) profiles in 25 subjects with molecularly confirmed clinical diagnosis of FA (FANCA complementation group) using Illumina's Infinium EPIC array. We identified 82 differentially methylated CpG sites that allow to distinguish subjects with FA from healthy individuals and subjects with other genetic disorders, defining an FA-specific DNAm signature. The episignature was validated using a second cohort of subjects with FA involving different complementation groups, documenting broader genetic sensitivity and demonstrating its specificity using the EpiSign Knowledge Database. The episignature properly classified DNA samples obtained from bone marrow aspirates, demonstrating robustness. Using the selected probes, we trained a machine-learning model able to classify EPIC DNAm profiles in molecularly unsolved cases. Finally, we show that the generated episignature includes CpG sites that do not undergo functional selective pressure, allowing diagnosis of FA in individuals with reverted phenotype due to gene conversion. These findings provide a tool to accelerate diagnostic testing in FA and broaden the clinical utility of DNAm profiling in the diagnostic setting.

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

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

Precision medicine in rare diseases: What is next?

Molecular diagnostics is a cornerstone of modern precision medicine, broadly understood as tailoring an individual's treatment, follow-up, and care based on molecular data. In rare diseases (RDs), molecular diagnoses reveal valuable information about the cause of symptoms, disease progression, familial risk, and in certain cases, unlock access to targeted therapies. Due to decreasing DNA sequencing costs, genome sequencing (GS) is emerging as the primary method for precision diagnostics in RDs. Several ongoing European initiatives for precision medicine have chosen GS as their method of choice. Recent research supports the role for GS as first-line genetic investigation in individuals with suspected RD, due to its improved diagnostic yield compared to other methods. Moreover, GS can detect a broad range of genetic aberrations including those in noncoding regions, producing comprehensive data that can be periodically reanalyzed for years to come when further evidence emerges. Indeed, targeted drug development and repurposing of medicines can be accelerated as more individuals with RDs receive a molecular diagnosis. Multidisciplinary teams in which clinical specialists collaborate with geneticists, genomics education of professionals and the public, and dialogue with patient advocacy groups are essential elements for the integration of precision medicine into clinical practice worldwide. It is also paramount that large research projects share genetic data and leverage novel technologies to fully diagnose individuals with RDs. In conclusion, GS increases diagnostic yields and is a crucial step toward precision medicine for RDs. Its clinical implementation will enable better patient management, unlock targeted therapies, and guide the development of innovative treatments.

Downstream Assays for Variant Resolution: Epigenetics, RNA Sequnncing, and Metabolomics

As the availability of advanced molecular testing like whole exome and genome sequencing expands, it comes with the added complication of interpreting inconclusive results, including determining the relevance of variants of uncertain significance or failing to find a variant in an otherwise suspected specific genetic disorder. This complication necessitates the use of alternative testing methods to gather more information in support of, or against, a particular genetic diagnosis. Therefore, new genome-wide approaches, including DNA epigenetic testing, RNA sequencing, and metabolomics, are increasingly being used to increase the diagnostic yield when used in conjunction with more conventional genetic tests.

Pathogenic variants in SMARCA1 cause an X-linked neurodevelopmental disorder modulated by NURF complex composition

Pathogenic variants in ATP-dependent chromatin remodeling proteins are a recurrent cause of neurodevelopmental disorders (NDDs). The NURF complex consists of BPTF and either the SNF2H (SMARCA5) or SNF2L (SMARCA1) ISWI-chromatin remodeling enzyme. Pathogenic variants in BPTF and SMARCA5 were previously implicated in NDDs. Here, we describe 40 individuals from 30 families with de novo or maternally inherited pathogenic variants in SMARCA1. This novel NDD was associated with mild to severe ID/DD, delayed or regressive speech development, and some recurrent facial dysmorphisms. Individuals carrying SMARCA1 loss-of-function variants exhibited a mild genome-wide DNA methylation profile and a high penetrance of macrocephaly. Genetic dissection of the NURF complex using Smarca1, Smarca5, and Bptfsingle and double mouse knockouts revealed the importance of NURF composition and dosage for proper forebrain development. Finally, we propose that genetic alterations affecting different NURF components result in a NDD with a broad clinical spectrum.

Heterozygous loss-of-function SMC3 variants are associated with variable and incompletely penetrant growth and developmental features

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 13 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated a milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, some instead having intriguing symptomatologies with rational biological links to SMC3 including bone marrow failure, acute myeloid leukemia, and Coats retinal vasculopathy. Analyses of transcriptomic and epigenetic data suggest that SMC3 pLoF variants reduce SMC3 expression but do not result in a blood DNA methylation signature clustering with that of CdLS, and that the global transcriptional signature of SMC3 loss is model-dependent. Our finding of substantial population-scale LoF intolerance in concert with variable penetrance in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multi-layered genomic data paired with careful phenotyping.

Complexity in Genetic Epilepsies: A Comprehensive Review

Epilepsy is a highly prevalent neurological disorder, affecting between 5-8 per 1000 individuals and is associated with a lifetime risk of up to 3%. In addition to high incidence, epilepsy is a highly heterogeneous disorder, with variation including, but not limited to the following: severity, age of onset, type of seizure, developmental delay, drug responsiveness, and other comorbidities. Variable phenotypes are reflected in a range of etiologies including genetic, infectious, metabolic, immune, acquired/structural (resulting from, for example, a severe head injury or stroke), or idiopathic. This review will focus specifically on epilepsies with a genetic cause, genetic testing, and biomarkers in epilepsy.

Synchronized long-read genome, methylome, epigenome, and transcriptome for resolving a Mendelian condition

Resolving the molecular basis of a Mendelian condition (MC) remains challenging owing to the diverse mechanisms by which genetic variants cause disease. To address this, we developed a synchronized long-read genome, methylome, epigenome, and transcriptome sequencing approach, which enables accurate single-nucleotide, insertion-deletion, and structural variant calling and diploid de novo genome assembly, and permits the simultaneous elucidation of haplotype-resolved CpG methylation, chromatin accessibility, and full-length transcript information in a single long-read sequencing run. Application of this approach to an Undiagnosed Diseases Network (UDN) participant with a chromosome X;13 balanced translocation of uncertain significance revealed that this translocation disrupted the functioning of four separate genes (NBEA, PDK3, MAB21L1, and RB1) previously associated with single-gene MCs. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included nonsense-mediated decay, fusion transcript formation, enhancer adoption, transcriptional readthrough silencing, and inappropriate X chromosome inactivation of autosomal genes. Overall, this highlights the utility of synchronized long-read multi-omic profiling for mechanistically resolving complex phenotypes.

Low-input and single-cell methods for Infinium DNA methylation BeadChips

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection p -values calculation achieved higher sensitivities for low-input datasets and was validated in over 100,000 public datasets with diverse methylation profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

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.

Integrative omics approaches to advance rare disease diagnostics

Over the past decade high-throughput DNA sequencing approaches, namely whole exome and whole genome sequencing became a standard procedure in Mendelian disease diagnostics. Implementation of these technologies greatly facilitated diagnostics and shifted the analysis paradigm from variant identification to prioritisation and evaluation. The diagnostic rates vary widely depending on the cohort size, heterogeneity and disease and range from around 30% to 50% leaving the majority of patients undiagnosed. Advances in omics technologies and computational analysis provide an opportunity to increase these unfavourable rates by providing evidence for disease-causing variant validation and prioritisation. This review aims to provide an overview of the current application of several omics technologies including RNA-sequencing, proteomics, metabolomics and DNA-methylation profiling for diagnostics of rare genetic diseases in general and inborn errors of metabolism in particular.

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

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

Danon Disease: Entire LAMP2 Gene Deletion with Unusual Clinical Presentation-Case Report and Review of the Literature

Danon disease is a rare x-linked dominant multisystemic disorder with a clinical triad of severe cardiomyopathy, skeletal myopathy, and intellectual disability. It is caused by defects in the lysosome-associated membrane protein-2 (LAMP2) gene. Numerous different mutations in the LAMP2 protein have been described. Danon disease is typically lethal by the mid-twenties in male patients due to cardiomyopathy and heart failure. Female patients usually present with milder and variable symptoms. This report describes a 42-year-old father and his 3-year-old daughter presenting with mild manifestations of the disease. The father has normal intellectual development and normal physical activity. At the age of 13, he was diagnosed with mild ventricular pre-excitation known as Wolf-Parkinson-White syndrome (WPWs), very mild and mostly asymptomatic cardiomyopathy and left ventricular hypertrophy, and at about the age of 25 presented with visual impairment due to cone-rod dystrophy. His daughter showed normal development and very mild asymptomatic electrocardiographic WPWs abnormalities with left mild ventricular hypertrophy. Genetic testing revealed an Xq24 microdeletion encompassing the entire LAMP2 gene. Relevant literature was reviewed as a reference for the etiology, diagnosis, treatment and case management.

A systematic review of neuroimaging epigenetic research: calling for an increased focus on development

Epigenetic mechanisms, such as DNA methylation (DNAm), have gained increasing attention as potential biomarkers and mechanisms underlying risk for neurodevelopmental, psychiatric and other brain-based disorders. Yet, surprisingly little is known about the extent to which DNAm is linked to individual differences in the brain itself, and how these associations may unfold across development - a time of life when many of these disorders emerge. Here, we systematically review evidence from the nascent field of Neuroimaging Epigenetics, combining structural or functional neuroimaging measures with DNAm, and the extent to which the developmental period (birth to adolescence) is represented in these studies. We identified 111 articles published between 2011-2021, out of which only a minority (21%) included samples under 18 years of age. Most studies were cross-sectional (85%), employed a candidate-gene approach (67%), and examined DNAm-brain associations in the context of health and behavioral outcomes (75%). Nearly half incorporated genetic data, and a fourth investigated environmental influences. Overall, studies support a link between peripheral DNAm and brain imaging measures, but there is little consistency in specific findings and it remains unclear whether DNAm markers present a cause, correlate or consequence of brain alterations. Overall, there is large heterogeneity in sample characteristics, peripheral tissue and brain outcome examined as well as the methods used. Sample sizes were generally low to moderate (median nall = 98, ndevelopmental = 80), and attempts at replication or meta-analysis were rare. Based on the strengths and weaknesses of existing studies, we propose three recommendations on how advance the field of Neuroimaging Epigenetics. We advocate for: (1) a greater focus on developmentally oriented research (i.e. pre-birth to adolescence); (2) the analysis of large, prospective, pediatric cohorts with repeated measures of DNAm and imaging to assess directionality; and (3) collaborative, interdisciplinary science to identify robust signals, triangulate findings and enhance translational potential.

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.