Megabase Length Hypermutation Accompanies Human Structural Variation at 17p11.2.

DNA rearrangements resulting in human genome structural variants (SVs) are caused by diverse mutational mechanisms. We used long- and short-read sequencing technologies to investigate end products of de novo chromosome 17p11.2 rearrangements and query the molecular mechanisms underlying both recurrent and non-recurrent events. Evidence for an increased rate of clustered single-nucleotide variant (SNV) mutation in cis with non-recurrent rearrangements was found. Indel and SNV formation are associated with both copy-number gains and losses of 17p11.2, occur up to ∼1 Mb away from the breakpoint junctions, and favor C > G transversion substitutions; results suggest that single-stranded DNA is formed during the genesis of the SV and provide compelling support for a microhomology-mediated break-induced replication (MMBIR) mechanism for SV formation. Our data show an additional mutational burden of MMBIR consisting of hypermutation confined to the locus and manifesting as SNVs and indels predominantly within genes.

Single Amino Acid Change Underlies Distinct Roles of H2A.Z Subtypes in Human Syndrome.

The developmental disorder Floating-Harbor syndrome (FHS) is caused by heterozygous truncating mutations in SRCAP, a gene encoding a chromatin remodeler mediating incorporation of histone variant H2A.Z. Here, we demonstrate that FHS-associated mutations result in loss of SRCAP nuclear localization, alter neural crest gene programs in human in vitro models and Xenopus embryos, and cause craniofacial defects. These defects are mediated by one of two H2A.Z subtypes, H2A.Z.2, whose knockdown mimics and whose overexpression rescues the FHS phenotype. Selective rescue by H2A.Z.2 is conferred by one of the three amino acid differences between the H2A.Z subtypes, S38/T38. We further show that H2A.Z.1 and H2A.Z.2 genomic occupancy patterns are qualitatively similar, but quantitatively distinct, and H2A.Z.2 incorporation at AT-rich enhancers and expression of their associated genes are both sensitized to SRCAP truncations. Altogether, our results illuminate the mechanism underlying a human syndrome and uncover selective functions of H2A.Z subtypes during development.

Human CLP1 mutations alter tRNA biogenesis, affecting both peripheral and central nervous system function.

CLP1 is a RNA kinase involved in tRNA splicing. Recently, CLP1 kinase-dead mice were shown to display a neuromuscular disorder with loss of motor neurons and muscle paralysis. Human genome analyses now identified a CLP1 homozygous missense mutation (p.R140H) in five unrelated families, leading to a loss of CLP1 interaction with the tRNA splicing endonuclease (TSEN) complex, largely reduced pre-tRNA cleavage activity, and accumulation of linear tRNA introns. The affected individuals develop severe motor-sensory defects, cortical dysgenesis, and microcephaly. Mice carrying kinase-dead CLP1 also displayed microcephaly and reduced cortical brain volume due to the enhanced cell death of neuronal progenitors that is associated with reduced numbers of cortical neurons. Our data elucidate a neurological syndrome defined by CLP1 mutations that impair tRNA splicing. Reduction of a founder mutation to homozygosity illustrates the importance of rare variations in disease and supports the clan genomics hypothesis.

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.

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

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

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.

Regulation of RAF protein kinases in ERK signalling.

RAF family kinases were among the first oncoproteins to be described more than 30 years ago. They primarily act as signalling relays downstream of RAS, and their close ties to cancer have fuelled a large number of studies. However, we still lack a systems-level understanding of their regulation and mode of action. The recent discovery that the catalytic activity of RAF depends on an allosteric mechanism driven by kinase domain dimerization is providing a vital new piece of information towards a comprehensive model of RAF function. The fact that current RAF inhibitors unexpectedly induce ERK signalling by stimulating RAF dimerization also calls for a deeper structural characterization of this family of kinases.

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.

Defects in diffusion barrier function of ciliary transition zone caused by ciliopathy variations of TMEM218.

Primary cilia are antenna-like structures protruding from the surface of various eukaryotic cells, and have distinct protein compositions in their membranes. This distinct protein composition is maintained by the presence of the transition zone (TZ) at the ciliary base, which acts as a diffusion barrier between the ciliary and plasma membranes. Defects in cilia and the TZ are known to cause a group of disorders collectively called the ciliopathies, which demonstrate a broad spectrum of clinical features, such as perinatally lethal Meckel syndrome (MKS), relatively mild Joubert syndrome (JBTS), and nonsyndromic nephronophthisis (NPHP). Proteins constituting the TZ can be grouped into the MKS and NPHP modules. The MKS module is composed of several transmembrane proteins and three soluble proteins. TMEM218 was recently reported to be mutated in individuals diagnosed as MKS and JBTS. However, little is known about how TMEM218 mutations found in MKS and JBTS affect the functions of cilia. In this study, we found that ciliary membrane proteins were not localized to cilia in TMEM218-knockout cells, indicating impaired barrier function of the TZ. Furthermore, the exogenous expression of JBTS-associated TMEM218 variants but not MKS-associated variants in TMEM218-knockout cells restored the localization of ciliary membrane proteins. In particular, when expressed in TMEM218-knockout cells, the TMEM218(R115H) variant found in JBTS was able to restore the barrier function of cells, whereas the MKS variant TMEM218(R115C) could not. Thus, the severity of symptoms of MKS and JBTS individuals appears to correlate with the degree of their ciliary defects at the cellular level.

CEP41, a ciliopathy-linked centrosomal protein, regulates microtubule assembly and cell proliferation.

Centrosomal proteins play pivotal roles in orchestrating microtubule dynamics, and their dysregulation leads to disorders, including cancer and ciliopathies. Understanding the multifaceted roles of centrosomal proteins is vital to comprehend their involvement in disease development. Here, we report novel cellular functions of CEP41, a centrosomal and ciliary protein implicated in Joubert syndrome. We show that CEP41 is an essential microtubule-associated protein with microtubule-stabilizing activity. Purified CEP41 binds to preformed microtubules, promotes microtubule nucleation and suppresses microtubule disassembly. When overexpressed in cultured cells, CEP41 localizes to microtubules and promotes microtubule bundling. Conversely, shRNA-mediated knockdown of CEP41 disrupts the interphase microtubule network and delays microtubule reassembly, emphasizing its role in microtubule organization. Further, we demonstrate that the association of CEP41 with microtubules relies on its conserved rhodanese homology domain (RHOD) and the N-terminal region. Interestingly, a disease-causing mutation in the RHOD domain impairs CEP41-microtubule interaction. Moreover, depletion of CEP41 inhibits cell proliferation and disrupts cell cycle progression, suggesting its potential involvement in cell cycle regulation. These insights into the cellular functions of CEP41 hold promise for unraveling the impact of its mutations in ciliopathies.

Recapitulating the human segmentation clock with pluripotent stem cells.

Pluripotent stem cells are increasingly used to model different aspects of embryogenesis and organ formation1. Despite recent advances in in vitro induction of major mesodermal lineages and cell types2,3, experimental model systems that can recapitulate more complex features of human mesoderm development and patterning are largely missing. Here we used induced pluripotent stem cells for the stepwise in vitro induction of presomitic mesoderm and its derivatives to model distinct aspects of human somitogenesis. We focused initially on modelling the human segmentation clock, a major biological concept believed to underlie the rhythmic and controlled emergence of somites, which give rise to the segmental pattern of the vertebrate axial skeleton. We observed oscillatory expression of core segmentation clock genes, including HES7 and DKK1, determined the period of the human segmentation clock to be around five hours, and demonstrated the presence of dynamic travelling-wave-like gene expression in in vitro-induced human presomitic mesoderm. Furthermore, we identified and compared oscillatory genes in human and mouse presomitic mesoderm derived from pluripotent stem cells, which revealed species-specific and shared molecular components and pathways associated with the putative mouse and human segmentation clocks. Using CRISPR-Cas9-based genome editing technology, we then targeted genes for which mutations in patients with segmentation defects of the vertebrae, such as spondylocostal dysostosis, have been reported (HES7, LFNG, DLL3 and MESP2). Subsequent analysis of patient-like and patient-derived induced pluripotent stem cells revealed gene-specific alterations in oscillation, synchronization or differentiation properties. Our findings provide insights into the human segmentation clock as well as diseases associated with human axial skeletogenesis.

Characterizing the molecular impact of KMT2D variants on the epigenetic and transcriptional landscapes in Kabuki syndrome.

Kabuki syndrome (KS) is a rare, multisystem disorder with a variable clinical phenotype. The majority of KS is caused by dominant loss-of-function mutations in KMT2D (lysine methyltransferase 2D). KMT2D mediates chromatin accessibility by adding methyl groups to lysine residue 4 of histone 3, which plays a critical role in cell differentiation and homeostasis. The molecular underpinnings of KS remain elusive partly because of a lack of histone modification data from human samples. Consequently, we profiled and characterized alterations in histone modification and gene transcription in peripheral blood mononuclear cells (PBMCs) from 33 patients with KMT2D mutations and 36 unaffected healthy controls. Our analysis identified unique enhancer signatures in H3K4me1 and H3K4me2 in KS compared with controls. Reduced enhancer signals were present for promoter-distal sites of immune-related genes for which co-binding of PBMC-specific transcription factors was predicted; 31% of super-enhancers of normal blood cells overlapped with disrupted enhancers in KS, supporting an association of reduced enhancer activity of immune-related genes with immune deficiency phenotypes. In contrast, increased enhancer signals were observed for promoter-proximal regions of metabolic genes enriched with EGR1 and E2F2 motifs, whose transcriptional levels were significantly increased in KS. Additionally, we identified ~100 de novo enhancers in genes, such as in MYO1F and AGAP2. Together, our results underscore the effect of KMT2D haploinsufficiency on dysregulation of enhancer states and gene transcription and provide a framework for the identification of therapeutic targets and biomarkers in preparation for clinical trial readiness.

COL11A2 as a candidate gene for vertebral malformations and congenital scoliosis.

Human vertebral malformations (VMs) have an estimated incidence of 1/2000 and are associated with significant health problems including congenital scoliosis (CS) and recurrent organ system malformation syndromes such as VACTERL (vertebral anomalies; anal abnormalities; cardiac abnormalities; tracheo-esophageal fistula; renal anomalies; limb anomalies). The genetic cause for the vast majority of VMs are unknown. In a CS/VM patient cohort, three COL11A2 variants (R130W, R1407L and R1413H) were identified in two patients with cervical VM. A third patient with a T9 hemivertebra and the R130W variant was identified from a separate study. These substitutions are predicted to be damaging to protein function, and R130 and R1407 residues are conserved in zebrafish Col11a2. To determine the role for COL11A2 in vertebral development, CRISPR/Cas9 was used to create a nonsense mutation (col11a2L642*) as well as a full gene locus deletion (col11a2del) in zebrafish. Both col11a2L642*/L642* and col11a2del/del mutant zebrafish exhibit vertebral fusions in the caudal spine, which form due to mineralization across intervertebral segments. To determine the functional consequence of VM-associated variants, we assayed their ability to suppress col11a2del VM phenotypes following transgenic expression within the developing spine. While wildtype col11a2 expression suppresses fusions in col11a2del/+ and col11a2del/del backgrounds, patient missense variant-bearing col11a2 failed to rescue the loss-of-function phenotype in these animals. These results highlight an essential role for COL11A2 in vertebral development and support a pathogenic role for two missense variants in CS.

A non-coding insertional mutation of Grhl2 causes gene over-expression and multiple structural anomalies including cleft palate, spina bifida and encephalocele.

Orofacial clefts, including cleft lip and palate (CL/P) and neural tube defects (NTDs) are among the most common congenital anomalies, but knowledge of the genetic basis of these conditions remains incomplete. The extent to which genetic risk factors are shared between CL/P, NTDs and related anomalies is also unclear. While identification of causative genes has largely focused on coding and loss of function mutations, it is hypothesized that regulatory mutations account for a portion of the unidentified heritability. We found that excess expression of Grainyhead-like 2 (Grhl2) causes not only spinal NTDs in Axial defects (Axd) mice but also multiple additional defects affecting the cranial region. These include orofacial clefts comprising midline cleft lip and palate and abnormalities of the craniofacial bones and frontal and/or basal encephalocele, in which brain tissue herniates through the cranium or into the nasal cavity. To investigate the causative mutation in the Grhl2Axd strain, whole genome sequencing identified an approximately 4 kb LTR retrotransposon insertion that disrupts the non-coding regulatory region, lying approximately 300 base pairs upstream of the 5' UTR. This insertion also lies within a predicted long non-coding RNA, oriented on the reverse strand, which like Grhl2 is over-expressed in Axd (Grhl2Axd) homozygous mutant embryos. Initial analysis of the GRHL2 upstream region in individuals with NTDs or cleft palate revealed rare or novel variants in a small number of cases. We hypothesize that mutations affecting the regulation of GRHL2 may contribute to craniofacial anomalies and NTDs in humans.

ANKRD11 pathogenic variants and 16q24.3 microdeletions share an altered DNA methylation signature in patients with KBG syndrome.

Pathogenic variants in ANKRD11 or microdeletions at 16q24.3 are the cause of KBG syndrome (KBGS), a neurodevelopmental syndrome characterized by intellectual disability, dental and skeletal anomalies, and characteristic facies. The ANKRD11 gene encodes the ankyrin repeat-containing protein 11A transcriptional regulator, which is expressed in the brain and implicated in neural development. Syndromic conditions caused by pathogenic variants in epigenetic regulatory genes show unique patterns of DNA methylation (DNAm) in peripheral blood, termed DNAm signatures. Given ANKRD11's role in chromatin modification, we tested whether pathogenic ANKRD11 variants underlying KBGS are associated with a DNAm signature. We profiled whole-blood DNAm in 21 individuals with ANKRD11 variants, 2 individuals with microdeletions at 16q24.3 and 28 typically developing individuals, using Illumina's Infinium EPIC array. We identified 95 differentially methylated CpG sites that distinguished individuals with KBGS and pathogenic variants in ANKRD11 (n = 14) from typically developing controls (n = 28). This DNAm signature was then validated in an independent cohort of seven individuals with KBGS and pathogenic ANKRD11 variants. We generated a machine learning model from the KBGS DNAm signature and classified the DNAm profiles of four individuals with variants of uncertain significance (VUS) in ANKRD11. We identified an intermediate classification score for an inherited missense variant transmitted from a clinically unaffected mother to her affected child. In conclusion, we show that the DNAm profiles of two individuals with 16q24.3 microdeletions were indistinguishable from the DNAm profiles of individuals with pathogenic variants in ANKRD11, and we demonstrate the diagnostic utility of the new KBGS signature by classifying the DNAm profiles of individuals with VUS in ANKRD11.

The Joubert-Meckel-Nephronophthisis Spectrum of Ciliopathies.

The Joubert syndrome (JS), Meckel syndrome (MKS), and nephronophthisis (NPH) ciliopathy spectrum could be the poster child for advances and challenges in Mendelian human genetics over the past half century. Progress in understanding these conditions illustrates many core concepts of human genetics. The JS phenotype alone is caused by pathogenic variants in more than 40 genes; remarkably, all of the associated proteins function in and around the primary cilium. Primary cilia are near-ubiquitous, microtubule-based organelles that play crucial roles in development and homeostasis. Protruding from the cell, these cellular antennae sense diverse signals and mediate Hedgehog and other critical signaling pathways. Ciliary dysfunction causes many human conditions termed ciliopathies, which range from multiple congenital malformations to adult-onset single-organ failure. Research on the genetics of the JS-MKS-NPH spectrum has spurred extensive functional work exploring the broadly important role of primary cilia in health and disease. This functional work promises to illuminate the mechanisms underlying JS-MKS-NPH in humans, identify therapeutic targets across genetic causes, and generate future precision treatments.

Telomeres and telomerase: three decades of progress.

Many recent advances have emerged in the telomere and telomerase fields. This Timeline article highlights the key advances that have expanded our views on the mechanistic underpinnings of telomeres and telomerase and their roles in ageing and disease. Three decades ago, the classic view was that telomeres protected the natural ends of linear chromosomes and that telomerase was a specific telomere-terminal transferase necessary for the replication of chromosome ends in single-celled organisms. While this concept is still correct, many diverse fields associated with telomeres and telomerase have substantially matured. These areas include the discovery of most of the key molecular components of telomerase, implications for limits to cellular replication, identification and characterization of human genetic disorders that result in premature telomere shortening, the concept that inhibiting telomerase might be a successful therapeutic strategy and roles for telomeres in regulating gene expression. We discuss progress in these areas and conclude with challenges and unanswered questions in the field.

The RIPK4-IRF6 signalling axis safeguards epidermal differentiation and barrier function.

The integrity of the mammalian epidermis depends on a balance of proliferation and differentiation in the resident population of stem cells1. The kinase RIPK4 and the transcription factor IRF6 are mutated in severe developmental syndromes in humans, and mice lacking these genes display epidermal hyperproliferation and soft-tissue fusions that result in neonatal lethality2-5. Our understanding of how these genes control epidermal differentiation is incomplete. Here we show that the role of RIPK4 in mouse development requires its kinase activity; that RIPK4 and IRF6 expressed in the epidermis regulate the same biological processes; and that the phosphorylation of IRF6 at Ser413 and Ser424 primes IRF6 for activation. Using RNA sequencing (RNA-seq), histone chromatin immunoprecipitation followed by sequencing (ChIP-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) of skin in wild-type and IRF6-deficient mouse embryos, we define the transcriptional programs that are regulated by IRF6 during epidermal differentiation. IRF6 was enriched at bivalent promoters, and IRF6 deficiency caused defective expression of genes that are involved in the metabolism of lipids and the formation of tight junctions. Accordingly, the lipid composition of the stratum corneum of Irf6-/- skin was abnormal, culminating in a severe defect in the function of the epidermal barrier. Collectively, our results explain how RIPK4 and IRF6 function to ensure the integrity of the epidermis and provide mechanistic insights into why developmental syndromes that are characterized by orofacial, skin and genital abnormalities result when this axis goes awry.

Extending the clinical spectrum of X-linked Tonne-Kalscheuer syndrome (TOKAS): new insights from the fetal perspective.

INTRODUCTION: Tonne-Kalscheuer syndrome (TOKAS) is a recessive X-linked multiple congenital anomaly disorder caused by RLIM variations. Of the 41 patients reported, only 7 antenatal cases were described. METHOD: After the antenatal diagnosis of TOKAS by exome analysis in a family followed for over 35 years because of multiple congenital anomalies in five male fetuses, a call for collaboration was made, resulting in a cohort of 11 previously unpublished cases. RESULTS: We present a TOKAS antenatal cohort, describing 11 new cases in 6 French families. We report a high frequency of diaphragmatic hernia (9 of 11), differences in sex development (10 of 11) and various visceral malformations. We report some recurrent dysmorphic features, but also pontocerebellar hypoplasia, pre-auricular skin tags and olfactory bulb abnormalities previously unreported in the literature. Although no clear genotype-phenotype correlation has yet emerged, we show that a recurrent p.(Arg611Cys) variant accounts for 66% of fetal TOKAS cases. We also report two new likely pathogenic variants in RLIM, outside of the two previously known mutational hotspots. CONCLUSION: Overall, we present the first fetal cohort of TOKAS, describe the clinical features that made it a recognisable syndrome at fetopathological examination, and extend the phenotypical spectrum and the known genotype of this rare disorder.

Missense variants causing Wiedemann-Steiner syndrome preferentially occur in the KMT2A-CXXC domain and are accurately classified using AlphaFold2.

Wiedemann-Steiner syndrome (WDSTS) is a neurodevelopmental disorder caused by de novo variants in KMT2A, which encodes a multi-domain histone methyltransferase. To gain insight into the currently unknown pathogenesis of WDSTS, we examined the spatial distribution of likely WDSTS-causing variants across the 15 different domains of KMT2A. Compared to variants in healthy controls, WDSTS variants exhibit a 61.9-fold overrepresentation within the CXXC domain-which mediates binding to unmethylated CpGs-suggesting a major role for this domain in mediating the phenotype. In contrast, we find no significant overrepresentation within the catalytic SET domain. Corroborating these results, we find that hippocampal neurons from Kmt2a-deficient mice demonstrate disrupted histone methylation (H3K4me1 and H3K4me3) preferentially at CpG-rich regions, but this has no systematic impact on gene expression. Motivated by these results, we combine accurate prediction of the CXXC domain structure by AlphaFold2 with prior biological knowledge to develop a classification scheme for missense variants in the CXXC domain. Our classifier achieved 92.6% positive and 92.9% negative predictive value on a hold-out test set. This classification performance enabled us to subsequently perform an in silico saturation mutagenesis and classify a total of 445 variants according to their functional effects. Our results yield a novel insight into the mechanistic basis of WDSTS and provide an example of how AlphaFold2 can contribute to the in silico characterization of variant effects with very high accuracy, suggesting a paradigm potentially applicable to many other Mendelian disorders.