Non-coding region variants upstream of MEF2C cause severe developmental disorder through three distinct loss-of-function mechanisms.

Clinical genetic testing of protein-coding regions identifies a likely causative variant in only around half of developmental disorder (DD) cases. The contribution of regulatory variation in non-coding regions to rare disease, including DD, remains very poorly understood. We screened 9,858 probands from the Deciphering Developmental Disorders (DDD) study for de novo mutations in the 5' untranslated regions (5' UTRs) of genes within which variants have previously been shown to cause DD through a dominant haploinsufficient mechanism. We identified four single-nucleotide variants and two copy-number variants upstream of MEF2C in a total of ten individual probands. We developed multiple bespoke and orthogonal experimental approaches to demonstrate that these variants cause DD through three distinct loss-of-function mechanisms, disrupting transcription, translation, and/or protein function. These non-coding region variants represent 23% of likely diagnoses identified in MEF2C in the DDD cohort, but these would all be missed in standard clinical genetics approaches. Nonetheless, these variants are readily detectable in exome sequence data, with 30.7% of 5' UTR bases across all genes well covered in the DDD dataset. Our analyses show that non-coding variants upstream of genes within which coding variants are known to cause DD are an important cause of severe disease and demonstrate that analyzing 5' UTRs can increase diagnostic yield. We also show how non-coding variants can help inform both the disease-causing mechanism underlying protein-coding variants and dosage tolerance of the gene.

Reassessing the association: Evaluation of a polyalanine deletion variant of RUNX2 in non-syndromic sagittal and metopic craniosynostosis.

The RUNT-related transcription factor RUNX2 plays a critical role in osteoblast differentiation, and alterations to gene dosage cause distinct craniofacial anomalies. Uniquely amongst the RUNT-related family, vertebrate RUNX2 encodes a polyglutamine/polyalanine repeat (Gln23-Glu-Ala17 in humans), with the length of the polyalanine component completely conserved in great apes. Surprisingly, a frequent 6-amino acid deletion polymorphism, p.(Ala84_Ala89)del, occurs in humans (termed 11A allele), and a previous association study (Cuellar et al. Bone 137:115395;2020) reported that the 11A variant was significantly more frequent in non-syndromic sagittal craniosynostosis (nsSag; allele frequency [AF] = 0.156; 95% confidence interval [CI] 0.126-0.189) compared to non-syndromic metopic craniosynostosis (nsMet; AF = 0.068; 95% CI 0.045-0.098). However, the gnomAD v.2.1.1 control population used by Cuellar et al. did not display Hardy-Weinberg equilibrium, hampering interpretation. To re-examine this association, we genotyped the RUNX2 11A polymorphism in 225 individuals with sporadic nsSag as parent-child trios and 164 singletons with sporadic nsMet, restricting our analysis to individuals of European ancestry. We compared observed allele frequencies to the non-transmitted alleles in the parent-child trios, and to the genome sequencing data from gnomAD v.4, which display Hardy-Weinberg equilibrium. Observed AFs (and 95% CI) were 0.076 (0.053-0.104) in nsSag and 0.082 (0.055-0.118) in nsMet, compared with 0.062 (0.042-0.089) in non-transmitted parental alleles and 0.065 (0.063-0.067) in gnomAD v.4.0.0 non-Finnish European control genomes. In summary, we observed a non-significant excess, compared to gnomAD data, of 11A alleles in both nsSag (relative risk 1.18, 95% CI 0.83-1.67) and nsMet (relative risk 1.29, 95% CI 0.87-1.92), but we did not replicate the much higher excess of RUNX2 11A alleles in nsSag previously reported (p = 0.0001).

De Novo Variants in CNOT1, a Central Component of the CCR4-NOT Complex Involved in Gene Expression and RNA and Protein Stability, Cause Neurodevelopmental Delay.

CNOT1 is a member of the CCR4-NOT complex, which is a master regulator, orchestrating gene expression, RNA deadenylation, and protein ubiquitination. We report on 39 individuals with heterozygous de novo CNOT1 variants, including missense, splice site, and nonsense variants, who present with a clinical spectrum of intellectual disability, motor delay, speech delay, seizures, hypotonia, and behavioral problems. To link CNOT1 dysfunction to the neurodevelopmental phenotype observed, we generated variant-specific Drosophila models, which showed learning and memory defects upon CNOT1 knockdown. Introduction of human wild-type CNOT1 was able to rescue this phenotype, whereas mutants could not or only partially, supporting our hypothesis that CNOT1 impairment results in neurodevelopmental delay. Furthermore, the genetic interaction with autism-spectrum genes, such as ASH1L, DYRK1A, MED13, and SHANK3, was impaired in our Drosophila models. Molecular characterization of CNOT1 variants revealed normal CNOT1 expression levels, with both mutant and wild-type alleles expressed at similar levels. Analysis of protein-protein interactions with other members indicated that the CCR4-NOT complex remained intact. An integrated omics approach of patient-derived genomics and transcriptomics data suggested only minimal effects on endonucleolytic nonsense-mediated mRNA decay components, suggesting that de novo CNOT1 variants are likely haploinsufficient hypomorph or neomorph, rather than dominant negative. In summary, we provide strong evidence that de novo CNOT1 variants cause neurodevelopmental delay with a wide range of additional co-morbidities. Whereas the underlying pathophysiological mechanism warrants further analysis, our data demonstrate an essential and central role of the CCR4-NOT complex in human brain development.

Pathogenic variants in the paired-related homeobox 1 gene (PRRX1) cause craniosynostosis with incomplete penetrance.

PURPOSE: Studies have previously implicated PRRX1 in craniofacial development, including demonstration of murine Prrx1 expression in the preosteogenic cells of the cranial sutures. We investigated the role of heterozygous missense and loss-of-function (LoF) variants in PRRX1 associated with craniosynostosis. METHODS: Trio-based genome, exome, or targeted sequencing were used to screen PRRX1 in patients with craniosynostosis; immunofluorescence analyses were used to assess nuclear localization of wild-type and mutant proteins. RESULTS: Genome sequencing identified 2 of 9 sporadically affected individuals with syndromic/multisuture craniosynostosis, who were heterozygous for rare/undescribed variants in PRRX1. Exome or targeted sequencing of PRRX1 revealed a further 9 of 1449 patients with craniosynostosis harboring deletions or rare heterozygous variants within the homeodomain. By collaboration, 7 additional individuals (4 families) were identified with putatively pathogenic PRRX1 variants. Immunofluorescence analyses showed that missense variants within the PRRX1 homeodomain cause abnormal nuclear localization. Of patients with variants considered likely pathogenic, bicoronal or other multisuture synostosis was present in 11 of 17 cases (65%). Pathogenic variants were inherited from unaffected relatives in many instances, yielding a 12.5% penetrance estimate for craniosynostosis. CONCLUSION: This work supports a key role for PRRX1 in cranial suture development and shows that haploinsufficiency of PRRX1 is a relatively frequent cause of craniosynostosis.

Unexpected role of SIX1 variants in craniosynostosis: expanding the phenotype of SIX1-related disorders.

BACKGROUND: Pathogenic heterozygous SIX1 variants (predominantly missense) occur in branchio-otic syndrome (BOS), but an association with craniosynostosis has not been reported. METHODS: We investigated probands with craniosynostosis of unknown cause using whole exome/genome (n=628) or RNA (n=386) sequencing, and performed targeted resequencing of SIX1 in 615 additional patients. Expression of SIX1 protein in embryonic cranial sutures was examined in the Six1nLacZ/+ reporter mouse. RESULTS: From 1629 unrelated cases with craniosynostosis we identified seven different SIX1 variants (three missense, including two de novo mutations, and four nonsense, one of which was also present in an affected twin). Compared with population data, enrichment of SIX1 loss-of-function variants was highly significant (p=0.00003). All individuals with craniosynostosis had sagittal suture fusion; additionally four had bilambdoid synostosis. Associated BOS features were often attenuated; some carrier relatives appeared non-penetrant. SIX1 is expressed in a layer basal to the calvaria, likely corresponding to the dura mater, and in the mid-sagittal mesenchyme. CONCLUSION: Craniosynostosis is associated with heterozygous SIX1 variants, with possible enrichment of loss-of-function variants compared with classical BOS. We recommend screening of SIX1 in craniosynostosis, particularly when sagittal±lambdoid synostosis and/or any BOS phenotypes are present. These findings highlight the role of SIX1 in cranial suture homeostasis.

Missense Mutations in NKAP Cause a Disorder of Transcriptional Regulation Characterized by Marfanoid Habitus and Cognitive Impairment.

NKAP is a ubiquitously expressed nucleoplasmic protein that is currently known as a transcriptional regulatory molecule via its interaction with HDAC3 and spliceosomal proteins. Here, we report a disorder of transcriptional regulation due to missense mutations in the X chromosome gene, NKAP. These mutations are clustered in the C-terminal region of NKAP where NKAP interacts with HDAC3 and post-catalytic spliceosomal complex proteins. Consistent with a role for the C-terminal region of NKAP in embryogenesis, nkap mutant zebrafish with a C-terminally truncated NKAP demonstrate severe developmental defects. The clinical features of affected individuals are highly conserved and include developmental delay, hypotonia, joint contractures, behavioral abnormalities, Marfanoid habitus, and scoliosis. In affected cases, transcriptome analysis revealed the presence of a unique transcriptome signature, which is characterized by the downregulation of long genes with higher exon numbers. These observations indicate the critical role of NKAP in transcriptional regulation and demonstrate that perturbations of the C-terminal region lead to developmental defects in both humans and zebrafish.

Bi-allelic Loss-of-Function CACNA1B Mutations in Progressive Epilepsy-Dyskinesia.

The occurrence of non-epileptic hyperkinetic movements in the context of developmental epileptic encephalopathies is an increasingly recognized phenomenon. Identification of causative mutations provides an important insight into common pathogenic mechanisms that cause both seizures and abnormal motor control. We report bi-allelic loss-of-function CACNA1B variants in six children from three unrelated families whose affected members present with a complex and progressive neurological syndrome. All affected individuals presented with epileptic encephalopathy, severe neurodevelopmental delay (often with regression), and a hyperkinetic movement disorder. Additional neurological features included postnatal microcephaly and hypotonia. Five children died in childhood or adolescence (mean age of death: 9 years), mainly as a result of secondary respiratory complications. CACNA1B encodes the pore-forming subunit of the pre-synaptic neuronal voltage-gated calcium channel Cav2.2/N-type, crucial for SNARE-mediated neurotransmission, particularly in the early postnatal period. Bi-allelic loss-of-function variants in CACNA1B are predicted to cause disruption of Ca2+ influx, leading to impaired synaptic neurotransmission. The resultant effect on neuronal function is likely to be important in the development of involuntary movements and epilepsy. Overall, our findings provide further evidence for the key role of Cav2.2 in normal human neurodevelopment.

Identifying phenotypic expansions for congenital diaphragmatic hernia plus (CDH+) using DECIPHER data.

Congenital diaphragmatic hernia (CDH) can occur in isolation or in conjunction with other birth defects (CDH+). A molecular etiology can only be identified in a subset of CDH cases. This is due, in part, to an incomplete understanding of the genes that contribute to diaphragm development. Here, we used clinical and molecular data from 36 individuals with CDH+ who are cataloged in the DECIPHER database to identify genes that may play a role in diaphragm development and to discover new phenotypic expansions. Among this group, we identified individuals who carried putatively deleterious sequence or copy number variants affecting CREBBP, SMARCA4, UBA2, and USP9X. The role of these genes in diaphragm development was supported by their expression in the developing mouse diaphragm, their similarity to known CDH genes using data from a previously published and validated machine learning algorithm, and/or the presence of CDH in other individuals with their associated genetic disorders. Our results demonstrate how data from DECIPHER, and other public databases, can be used to identify new phenotypic expansions and suggest that CREBBP, SMARCA4, UBA2, and USP9X play a role in diaphragm development.

Dissection of contiguous gene effects for deletions around ERF on chromosome 19.

Heterozygous intragenic loss-of-function mutations of ERF, encoding an ETS transcription factor, were previously reported to cause a novel craniosynostosis syndrome, suggesting that ERF is haploinsufficient. We describe six families harboring heterozygous deletions including, or near to, ERF, of which four were characterized by whole-genome sequencing and two by chromosomal microarray. Based on the severity of associated intellectual disability (ID), we identify three categories of ERF-associated deletions. The smallest (32 kb) and only inherited deletion included two additional centromeric genes and was not associated with ID. Three larger deletions (264-314 kb) that included at least five further centromeric genes were associated with moderate ID, suggesting that deletion of one or more of these five genes causes ID. The individual with the most severe ID had a more telomerically extending deletion, including CIC, a known ID gene. Children found to harbor ERF deletions should be referred for craniofacial assessment, to exclude occult raised intracranial pressure.

Histone Lysine Methylases and Demethylases in the Landscape of Human Developmental Disorders.

Histone lysine methyltransferases (KMTs) and demethylases (KDMs) underpin gene regulation. Here we demonstrate that variants causing haploinsufficiency of KMTs and KDMs are frequently encountered in individuals with developmental disorders. Using a combination of human variation databases and existing animal models, we determine 22 KMTs and KDMs as additional candidates for dominantly inherited developmental disorders. We show that KMTs and KDMs that are associated with, or are candidates for, dominant developmental disorders tend to have a higher level of transcription, longer canonical transcripts, more interactors, and a higher number and more types of post-translational modifications than other KMT and KDMs. We provide evidence to firmly associate KMT2C, ASH1L, and KMT5B haploinsufficiency with dominant developmental disorders. Whereas KMT2C or ASH1L haploinsufficiency results in a predominantly neurodevelopmental phenotype with occasional physical anomalies, KMT5B mutations cause an overgrowth syndrome with intellectual disability. We further expand the phenotypic spectrum of KMT2B-related disorders and show that some individuals can have severe developmental delay without dystonia at least until mid-childhood. Additionally, we describe a recessive histone lysine-methylation defect caused by homozygous or compound heterozygous KDM5B variants and resulting in a recognizable syndrome with developmental delay, facial dysmorphism, and camptodactyly. Collectively, these results emphasize the significance of histone lysine methylation in normal human development and the importance of this process in human developmental disorders. Our results demonstrate that systematic clinically oriented pathway-based analysis of genomic data can accelerate the discovery of rare genetic disorders.

Mutations in the BAF-Complex Subunit DPF2 Are Associated with Coffin-Siris Syndrome.

Variants affecting the function of different subunits of the BAF chromatin-remodelling complex lead to various neurodevelopmental syndromes, including Coffin-Siris syndrome. Furthermore, variants in proteins containing PHD fingers, motifs recognizing specific histone tail modifications, have been associated with several neurological and developmental-delay disorders. Here, we report eight heterozygous de novo variants (one frameshift, two splice site, and five missense) in the gene encoding the BAF complex subunit double plant homeodomain finger 2 (DPF2). Affected individuals share common clinical features described in individuals with Coffin-Siris syndrome, including coarse facial features, global developmental delay, intellectual disability, speech impairment, and hypoplasia of fingernails and toenails. All variants occur within the highly conserved PHD1 and PHD2 motifs. Moreover, missense variants are situated close to zinc binding sites and are predicted to disrupt these sites. Pull-down assays of recombinant proteins and histone peptides revealed that a subset of the identified missense variants abolish or impaire DPF2 binding to unmodified and modified H3 histone tails. These results suggest an impairment of PHD finger structural integrity and cohesion and most likely an aberrant recognition of histone modifications. Furthermore, the overexpression of these variants in HEK293 and COS7 cell lines was associated with the formation of nuclear aggregates and the recruitment of both wild-type DPF2 and BRG1 to these aggregates. Expression analysis of truncating variants found in the affected individuals indicated that the aberrant transcripts escape nonsense-mediated decay. Altogether, we provide compelling evidence that de novo variants in DPF2 cause Coffin-Siris syndrome and propose a dominant-negative mechanism of pathogenicity.

Genotype-Phenotype Correlation in NF1: Evidence for a More Severe Phenotype Associated with Missense Mutations Affecting NF1 Codons 844-848.

Neurofibromatosis type 1 (NF1), a common genetic disorder with a birth incidence of 1:2,000-3,000, is characterized by a highly variable clinical presentation. To date, only two clinically relevant intragenic genotype-phenotype correlations have been reported for NF1 missense mutations affecting p.Arg1809 and a single amino acid deletion p.Met922del. Both variants predispose to a distinct mild NF1 phenotype with neither externally visible cutaneous/plexiform neurofibromas nor other tumors. Here, we report 162 individuals (129 unrelated probands and 33 affected relatives) heterozygous for a constitutional missense mutation affecting one of five neighboring NF1 codons-Leu844, Cys845, Ala846, Leu847, and Gly848-located in the cysteine-serine-rich domain (CSRD). Collectively, these recurrent missense mutations affect ∼0.8% of unrelated NF1 mutation-positive probands in the University of Alabama at Birmingham (UAB) cohort. Major superficial plexiform neurofibromas and symptomatic spinal neurofibromas were more prevalent in these individuals compared with classic NF1-affected cohorts (both p < 0.0001). Nearly half of the individuals had symptomatic or asymptomatic optic pathway gliomas and/or skeletal abnormalities. Additionally, variants in this region seem to confer a high predisposition to develop malignancies compared with the general NF1-affected population (p = 0.0061). Our results demonstrate that these NF1 missense mutations, although located outside the GAP-related domain, may be an important risk factor for a severe presentation. A genotype-phenotype correlation at the NF1 region 844-848 exists and will be valuable in the management and genetic counseling of a significant number of individuals.

Evaluating the performance of a clinical genome sequencing program for diagnosis of rare genetic disease, seen through the lens of craniosynostosis.

PURPOSE: Genome sequencing (GS) for diagnosis of rare genetic disease is being introduced into the clinic, but the complexity of the data poses challenges for developing pipelines with high diagnostic sensitivity. We evaluated the performance of the Genomics England 100,000 Genomes Project (100kGP) panel-based pipelines, using craniosynostosis as a test disease. METHODS: GS data from 114 probands with craniosynostosis and their relatives (314 samples), negative on routine genetic testing, were scrutinized by a specialized research team, and diagnoses compared with those made by 100kGP. RESULTS: Sixteen likely pathogenic/pathogenic variants were identified by 100kGP. Eighteen additional likely pathogenic/pathogenic variants were identified by the research team, indicating that for craniosynostosis, 100kGP panels had a diagnostic sensitivity of only 47%. Measures that could have augmented diagnoses were improved calling of existing panel genes (+18% sensitivity), review of updated panels (+12%), comprehensive analysis of de novo small variants (+29%), and copy-number/structural variants (+9%). Recent NHS England recommendations that partially incorporate these measures should achieve 85% overall sensitivity (+38%). CONCLUSION: GS identified likely pathogenic/pathogenic variants in 29.8% of previously undiagnosed patients with craniosynostosis. This demonstrates the value of research analysis and the importance of continually improving algorithms to maximize the potential of clinical GS.

DLG4-related synaptopathy: a new rare brain disorder.

PURPOSE: Postsynaptic density protein-95 (PSD-95), encoded by DLG4, regulates excitatory synaptic function in the brain. Here we present the clinical and genetic features of 53 patients (42 previously unpublished) with DLG4 variants. METHODS: The clinical and genetic information were collected through GeneMatcher collaboration. All the individuals were investigated by local clinicians and the gene variants were identified by clinical exome/genome sequencing. RESULTS: The clinical picture was predominated by early onset global developmental delay, intellectual disability, autism spectrum disorder, and attention deficit-hyperactivity disorder, all of which point to a brain disorder. Marfanoid habitus, which was previously suggested to be a characteristic feature of DLG4-related phenotypes, was found in only nine individuals and despite some overlapping features, a distinct facial dysmorphism could not be established. Of the 45 different DLG4 variants, 39 were predicted to lead to loss of protein function and the majority occurred de novo (four with unknown origin). The six missense variants identified were suggested to lead to structural or functional changes by protein modeling studies. CONCLUSION: The present study shows that clinical manifestations associated with DLG4 overlap with those found in other neurodevelopmental disorders of synaptic dysfunction; thus, we designate this group of disorders as DLG4-related synaptopathy.

Expanding the phenotype of biallelic loss-of-function variants in the NSUN2 gene: Description of four individuals with juvenile cataract, chronic nephritis, or brain anomaly as novel complications.

The NSUN2 gene encodes a tRNA cytosine methyltransferase that functions in the maturation of leucyl tRNA (Leu) (CAA) precursors, which is crucial for the anticodon-codon pairing and correct translation of mRNA. Biallelic loss of function variants in NSUN2 are known to cause moderate to severe intellectual disability. Microcephaly, postnatal growth retardation, and dysmorphic facial features are common complications in this genetic disorder, and delayed puberty is occasionally observed. Here, we report four individuals, two sets of siblings, with biallelic loss-of-function variants in the NSUN2 gene. The first set of siblings have compound heterozygous frameshift variants: c.546_547insCT, p.Met183Leufs*13; c.1583del, p.Pro528Hisfs*19, and the other siblings carry a homozygous frameshift variant: c.1269dup, p.Val424Cysfs*14. In addition to previously reported clinical features, the first set of siblings showed novel complications of juvenile cataract and chronic nephritis. The other siblings showed hypomyelination and simplified gyral pattern in neuroimaging. NSUN2-related intellectual disability is a very rare condition, and less than 20 cases have been reported previously. Juvenile cataract, chronic nephritis, and brain anomaly shown in the present patients have not been previously described. Our report suggests clinical diversity of NSUN2-related intellectual disability.

A genome-wide association study implicates the BMP7 locus as a risk factor for nonsyndromic metopic craniosynostosis.

Our previous genome-wide association study (GWAS) for sagittal nonsyndromic craniosynostosis (sNCS) provided important insights into the genetics of midline CS. In this study, we performed a GWAS for a second midline NCS, metopic NCS (mNCS), using 215 non-Hispanic white case-parent triads. We identified six variants with genome-wide significance (P ≤ 5 × 10-8): rs781716 (P = 4.71 × 10-9; odds ratio [OR] = 2.44) intronic to SPRY3; rs6127972 (P = 4.41 × 10-8; OR = 2.17) intronic to BMP7; rs62590971 (P = 6.22 × 10-9; OR = 0.34), located ~ 155 kb upstream from TGIF2LX; and rs2522623, rs2573826, and rs2754857, all intronic to PCDH11X (P = 1.76 × 10-8, OR = 0.45; P = 3.31 × 10-8, OR = 0.45; P = 1.09 × 10-8, OR = 0.44, respectively). We performed a replication study of these variants using an independent non-Hispanic white sample of 194 unrelated mNCS cases and 333 unaffected controls; only the association for rs6127972 (P = 0.004, OR = 1.45; meta-analysis P = 1.27 × 10-8, OR = 1.74) was replicated. Our meta-analysis examining single nucleotide polymorphisms common to both our mNCS and sNCS studies showed the strongest association for rs6127972 (P = 1.16 × 10-6). Our imputation analysis identified a linkage disequilibrium block encompassing rs6127972, which contained an enhancer overlapping a CTCF transcription factor binding site (chr20:55,798,821-55,798,917) that was significantly hypomethylated in mesenchymal stem cells derived from fused metopic compared to open sutures from the same probands. This study provides additional insights into genetic factors in midline CS.

De Novo Missense Mutations in DHX30 Impair Global Translation and Cause a Neurodevelopmental Disorder.

DHX30 is a member of the family of DExH-box helicases, which use ATP hydrolysis to unwind RNA secondary structures. Here we identified six different de novo missense mutations in DHX30 in twelve unrelated individuals affected by global developmental delay (GDD), intellectual disability (ID), severe speech impairment and gait abnormalities. While four mutations are recurrent, two are unique with one affecting the codon of one recurrent mutation. All amino acid changes are located within highly conserved helicase motifs and were found to either impair ATPase activity or RNA recognition in different in vitro assays. Moreover, protein variants exhibit an increased propensity to trigger stress granule (SG) formation resulting in global translation inhibition. Thus, our findings highlight the prominent role of translation control in development and function of the central nervous system and also provide molecular insight into how DHX30 dysfunction might cause a neurodevelopmental disorder.

YY1 Haploinsufficiency Causes an Intellectual Disability Syndrome Featuring Transcriptional and Chromatin Dysfunction.

Yin and yang 1 (YY1) is a well-known zinc-finger transcription factor with crucial roles in normal development and malignancy. YY1 acts both as a repressor and as an activator of gene expression. We have identified 23 individuals with de novo mutations or deletions of YY1 and phenotypic features that define a syndrome of cognitive impairment, behavioral alterations, intrauterine growth restriction, feeding problems, and various congenital malformations. Our combined clinical and molecular data define "YY1 syndrome" as a haploinsufficiency syndrome. Through immunoprecipitation of YY1-bound chromatin from affected individuals' cells with antibodies recognizing both ends of the protein, we show that YY1 deletions and missense mutations lead to a global loss of YY1 binding with a preferential retention at high-occupancy sites. Finally, we uncover a widespread loss of H3K27 acetylation in particular on the YY1-bound enhancers, underscoring a crucial role for YY1 in enhancer regulation. Collectively, these results define a clinical syndrome caused by haploinsufficiency of YY1 through dysregulation of key transcriptional regulators.

Correction: SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

PURPOSE: Enrichment of heterozygous missense and truncating SMAD6 variants was previously reported in nonsyndromic sagittal and metopic synostosis, and interaction of SMAD6 variants with a common polymorphism nearBMP2 (rs1884302) was proposed to contribute to inconsistent penetrance. We determined the occurrence of SMAD6 variants in all types of craniosynostosis, evaluated the impact of different missense variants on SMAD6 function, and tested independently whether rs1884302 genotype significantly modifies the phenotype. METHODS: We performed resequencing of SMAD6 in 795 unsolved patients with any type of craniosynostosis and genotyped rs1884302 in SMAD6-positive individuals and relatives. We examined the inhibitory activity and stability of SMAD6 missense variants. RESULTS: We found 18 (2.3%) different rare damaging SMAD6 variants, with the highest prevalence in metopic synostosis (5.8%) and an 18.3-fold enrichment of loss-of-function variants comparedwith gnomAD data (P < 10-7). Combined with eight additional variants, ≥20/26 were transmitted from an unaffected parent but rs1884302 genotype did not predict phenotype. CONCLUSION: Pathogenic SMAD6 variants substantially increase the risk of both nonsyndromic and syndromic presentations of craniosynostosis, especially metopic synostosis. Functional analysis is important to evaluate missense variants. Genotyping of rs1884302 is not clinically useful. Mechanisms to explain the remarkable diversity of phenotypes associated with SMAD6 variants remain obscure.