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.

HNRNPC haploinsufficiency affects alternative splicing of intellectual disability-associated genes and causes a neurodevelopmental disorder.

Heterogeneous nuclear ribonucleoprotein C (HNRNPC) is an essential, ubiquitously abundant protein involved in mRNA processing. Genetic variants in other members of the HNRNP family have been associated with neurodevelopmental disorders. Here, we describe 13 individuals with global developmental delay, intellectual disability, behavioral abnormalities, and subtle facial dysmorphology with heterozygous HNRNPC germline variants. Five of them bear an identical in-frame deletion of nine amino acids in the extreme C terminus. To study the effect of this recurrent variant as well as HNRNPC haploinsufficiency, we used induced pluripotent stem cells (iPSCs) and fibroblasts obtained from affected individuals. While protein localization and oligomerization were unaffected by the recurrent C-terminal deletion variant, total HNRNPC levels were decreased. Previously, reduced HNRNPC levels have been associated with changes in alternative splicing. Therefore, we performed a meta-analysis on published RNA-seq datasets of three different cell lines to identify a ubiquitous HNRNPC-dependent signature of alternative spliced exons. The identified signature was not only confirmed in fibroblasts obtained from an affected individual but also showed a significant enrichment for genes associated with intellectual disability. Hence, we assessed the effect of decreased and increased levels of HNRNPC on neuronal arborization and neuronal migration and found that either condition affects neuronal function. Taken together, our data indicate that HNRNPC haploinsufficiency affects alternative splicing of multiple intellectual disability-associated genes and that the developing brain is sensitive to aberrant levels of HNRNPC. Hence, our data strongly support the inclusion of HNRNPC to the family of HNRNP-related neurodevelopmental disorders.

EMILIN1 deficiency causes arterial tortuosity with osteopenia and connects impaired elastogenesis with defective collagen fibrillogenesis.

EMILIN1 (elastin-microfibril-interface-located-protein-1) is a structural component of the elastic fiber network and localizes to the interface between the fibrillin microfibril scaffold and the elastin core. How EMILIN1 contributes to connective tissue integrity is not fully understood. Here, we report bi-allelic EMILIN1 loss-of-function variants causative for an entity combining cutis laxa, arterial tortuosity, aneurysm formation, and bone fragility, resembling autosomal-recessive cutis laxa type 1B, due to EFEMP2 (FBLN4) deficiency. In both humans and mice, absence of EMILIN1 impairs EFEMP2 extracellular matrix deposition and LOX activity resulting in impaired elastogenesis, reduced collagen crosslinking, and aberrant growth factor signaling. Collagen fiber ultrastructure and histopathology in EMILIN1- or EFEMP2-deficient skin and aorta corroborate these findings and murine Emilin1-/- femora show abnormal trabecular bone formation and strength. Altogether, EMILIN1 connects elastic fiber network with collagen fibril formation, relevant for both bone and vascular tissue homeostasis.

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.

ASXL3-related disorder: Molecular phenotyping and comprehensive review providing insights into disease mechanism.

ASXL3-related disorder, sometimes referred to as Bainbridge-Ropers syndrome, was first identified as a distinct neurodevelopmental disorder by Bainbridge et al. in 2013. Since then, there have been a number of case series and single case reports published worldwide. A comprehensive review of the literature was carried out. Abstracts were screened, relevant literature was analysed, and descriptions of common phenotypic features were quantified. ASXL3 variants were collated and categorised. Common phenotypic features comprised global developmental delay or intellectual disability (97%), feeding problems (76%), hypotonia (88%) and characteristic facial features (93%). The majority of genetic variants were de novo truncating variants in exon 11 or 12 of the ASXL3 gene. Several gaps in our knowledge of this disorder were identified, namely, underlying pathophysiology and disease mechanism, disease contribution of missense variants, relevance of variant location, prevalence and penetrance data. Clinical information is currently limited by patient numbers and lack of longitudinal data, which this review aims to address.

Clinical and molecular characterization of patients with YWHAG-related epilepsy.

OBJECTIVE: YWHAG variant alleles have been associated with a rare disease trait whose clinical synopsis includes an early onset epileptic encephalopathy with predominantly myoclonic seizures, developmental delay/intellectual disability, and facial dysmorphisms. Through description of a large cohort, which doubles the number of reported patients, we further delineate the spectrum of YWHAG-related epilepsy. METHODS: We included in this study 24 patients, 21 new and three previously described, with pathogenic/likely pathogenic variants in YWHAG. We extended the analysis of clinical, electroencephalographic, brain magnetic resonance imaging, and molecular genetic information to 24 previously published patients. RESULTS: The phenotypic spectrum of YWHAG-related disorders ranges from mild developmental delay to developmental and epileptic encephalopathy (DEE). Epilepsy onset is in the first 2 years of life. Seizure freedom can be achieved in half of the patients (13/24, 54%). Intellectual disability (23/24, 96%), behavioral disorders (18/24, 75%), neurological signs (13/24, 54%), and dysmorphisms (6/24, 25%) are common. A genotype-phenotype correlation emerged, as DEE is more represented in patients with missense variants located in the ligand-binding domain than in those with truncating or missense variants in other domains (90% vs. 19%, p < .001). SIGNIFICANCE: This study suggests that pathogenic YWHAG variants cause a wide range of clinical presentations with variable severity, ranging from mild developmental delay to DEE. In this allelic series, a genotype-phenotype correlation begins to emerge, potentially providing prognostic information for clinical management and genetic counseling.

Familial Bainbridge-Ropers syndrome: Report of familial ASXL3 inheritance and a milder phenotype.

De novo truncating and splicing pathogenic variants in the Additional Sex Combs-Like 3 (ASXL3) gene are known to cause neurodevelopmental delay, intellectual disability, behavioral difficulties, hypotonia, feeding problems and characteristic facial features. We previously reported 45 patients with ASXL3-related disorder including three individuals with a familial variant. Here we report the detailed clinical and molecular characteristics of these three families with inherited ASXL3-related disorder. First, a father and son with c.2791_2792del p.Gln931fs pathogenic variant. The second, a mother, daughter and son with c.4534C > T, p.Gln1512Ter pathogenic variant. The third, a mother and her daughter with c.4441dup, p.Leu1481fs maternally inherited pathogenic variant. This report demonstrates intrafamilial phenotypic heterogeneity and confirms heritability of ASXL3-related disorder.

Expanding the phenotype of SPARC-related osteogenesis imperfecta: clinical findings in two patients with pathogenic variants in SPARC and literature review.

BACKGROUND: Secreted protein, acidic, cysteine rich (SPARC)-related osteogenesis imperfecta (OI), also referred to as OI type XVII, was first described in 2015, since then there has been only one further report of this form of OI. SPARC is located on chromosome 5 between bands q31 and q33. The encoded protein is necessary for calcification of the collagen in bone, synthesis of extracellular matrix and the promotion of changes to cell shape. METHODS: We describe a further two patients with previously unreported homozygous SPARC variants with OI: one splice site; one nonsense pathogenic variant. We present detailed information on the clinical and radiological phenotype and correlate this with their genotype. There are only two previous reports by Mendozo-Londono et al and Hayat et al with clinical descriptions of patients with SPARC variants. RESULTS: From the data we have obtained, common clinical features in individuals with OI type XVII caused by SPARC variants include scoliosis (5/5), vertebral compression fractures (5/5), multiple long bone fractures (5/5) and delayed motor development (3/3). Interestingly, 2/4 patients also had abnormal brain MRI, including high subcortical white matter changes, abnormal fluid-attenuated inversion in the para-atrial white matter and a large spinal canal from T10 to L1. Of significance, both patients reported here presented with significant neuromuscular weakness prompting early workup. CONCLUSION: Common phenotypic expressions include delayed motor development with neuromuscular weakness, scoliosis and multiple fractures. The data presented here broaden the phenotypic spectrum establishing similar patterns of neuromuscular presentation with a presumed diagnosis of 'myopathy'.

Report of two children with global developmental delay in association with de novo TLK2 variant and literature review.

We describe clinical details, including novel findings, of two further children with the newly defined TLK2-related disorder. One patient was recruited to the Deciphering Developmental Delay (DDD) Study to identify underlying etiology of global developmental delay. The other was detected on whole-exome sequencing as part of second line investigations following normal microarray. Both patients were found to have de novo heterozygous pathogenic TLK2 variants. A novel c.6del p.(Glu3Lysfs*) loss-of-function frameshift variant was found in Patient 1. A c.1121+1G>A splice-donor variant was detected in Patient 2. TLK2-related neurodevelopmental disorder is a specific syndrome that has been recently described. Global developmental delay, behavioral problems, gastrointestinal disorders, and typical facial dysmorphism are common features. Neuropsychiatric disorders, ophthalmic, musculoskeletal and cranial abnormalities, as well as short stature, have also all been described. The novel findings we describe include sleep disturbance, nondifferentiation of lateral semi-circular canals (where asymmetric semi-circular canals were a feature in the previous cohort), vesico-ureteric reflux, and bilateral periauricular skin tags. Here, we report a novel TLK2 variant and previously undescribed features of TLK2-related disorder, to expand the clinical phenotype and provide further genotype-phenotype correlation.

Expanding the phenotype of TAB2 variants and literature review.

TAB2 is a gene located on chromosome 6q25.1 and plays a key role in development of the heart. Existing literature describes congenital heart disease as a common recognized phenotype of TAB2 gene variants, with evidence of a distinct syndromic phenotype also existing beyond this. Here we describe 14 newly identified individuals with nine novel, pathogenic TAB2 variants. The majority of individuals were identified through the Deciphering Developmental Disorders study through trio whole exome sequencing. Eight individuals had de novo variants, the other six individuals were found to have maternally inherited, or likely maternally inherited, variants. Five individuals from the same family were identified following cardiac disease gene panel in the proband and subsequent targeted familial gene sequencing. The clinical features of this cohort were compared to the existing literature. Common clinical features include distinctive facial features, growth abnormalities, joint hypermobility, hypotonia, and developmental delay. Newly identified features included feeding difficulties, sleep problems, visual problems, genitourinary abnormality, and other anatomical variations. Here we report 14 new individuals, including novel TAB2 variants, in order to expand the emerging syndromic clinical phenotype and provide further genotype-phenotype correlation.

Expanding the phenotype of HNRNPU-related neurodevelopmental disorder with emphasis on seizure phenotype and review of literature.

Pathogenic variants in heterogeneous nuclear ribonucleoprotein U (HNRNPU) results in a novel neurodevelopmental disorder recently delineated. Here, we report on 17 previously unpublished patients carrying HNRNPU pathogenic variants. All patients were found to harbor de novo loss-of-function variants except for one individual where the inheritance could not be determined, as a parent was unavailable for testing. All patients had seizures which started in early childhood, global developmental delay, intellectual disability, and dysmorphic features. In addition, hypotonia, behavioral abnormalities (such as autistic features, aggression, anxiety, and obsessive-compulsive behaviors), and cardiac (septal defects) and/or brain abnormalities (ventriculomegaly and corpus callosum thinning/agenesis) were frequently observed. We have noted four recurrent variants in the literature (c.1089G>A p.(Trp363*), c.706_707del p.(Glu236Thrfs*6), c.847_857del p.(Phe283Serfs*5), and c.1681dels p.(Gln561Serfs*45)).

Further delineation of phenotypic spectrum of SCN2A-related disorder.

SCN2A-related disorders include intellectual disability, autism spectrum disorder, seizures, episodic ataxia, and schizophrenia. In this study, the phenotype-genotype association in SCN2A-related disorders was further delineated by collecting detailed clinical and molecular characteristics. Using previously proposed genotype-phenotype hypotheses based on variant function and position, the potential of phenotype prediction from the variants found was examined. Patients were identified through the Deciphering Developmental Disorders study and gene matching strategies. Phenotypic information and variant interpretation evidence were collated. Seventeen previously unreported patients and five patients who had been previously reported (but with minimal phenotypic and segregation data) were included (10 males, 12 females; median age 10.5 years). All patients had developmental delays and the majority had intellectual disabilities. Seizures were reported in 15 of 22 (68.2%), four of 22 (18.2%) had autism spectrum disorder and no patients were reported with episodic ataxia. The majority of variants were de novo. One family had presumed gonadal mosaicism. The correlation of the use of sodium channel-blocking antiepileptic drugs with phenotype or genotype was variable. These data suggest that variant type and position alone can provide some predictive information about the phenotype in a proportion of cases, but more precise assessment of variant function is needed for meaningful phenotype prediction.

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.

NBAS Variants Are Associated with Quantitative and Qualitative NK and B Cell Deficiency.

PURPOSE: Biallelic pathogenic NBAS variants manifest as a multisystem disorder with heterogeneous clinical phenotypes such as recurrent acute liver failure, growth retardation, and susceptibility to infections. This study explores how NBAS-associated disease affects cells of the innate and adaptive immune system. METHODS: Clinical and laboratory parameters were combined with functional multi-parametric immunophenotyping methods in fifteen NBAS-deficient patients to discover possible alterations in their immune system. RESULTS: Our study revealed reduced absolute numbers of mature CD56dim natural killer (NK) cells. Notably, the residual NK cell population in NBAS-deficient patients exerted a lower potential for activation and degranulation in response to K562 target cells, suggesting an NK cell-intrinsic role for NBAS in the release of cytotoxic granules. NBAS-deficient NK cell activation and degranulation was normalized upon pre-activation by IL-2 in vitro, suggesting that functional impairment was reversible. In addition, we observed a reduced number of naïve B cells in the peripheral blood associated with hypogammaglobulinemia. CONCLUSION: In summary, we demonstrate that pathogenic biallelic variants in NBAS are associated with dysfunctional NK cells as well as impaired adaptive humoral immunity.

Growth, development, and phenotypic spectrum of individuals with deletions of 2q33.1 involving SATB2.

SATB2-Associated syndrome (SAS) is an autosomal dominant, multisystemic, neurodevelopmental disorder due to alterations in SATB2 at 2q33.1. A limited number of individuals with 2q33.1 contiguous deletions encompassing SATB2 (ΔSAS) have been described in the literature. We describe 17 additional individuals with ΔSAS, review the phenotype of 33 previously published individuals with 2q33.1 deletions (n = 50, mean age = 8.5 ± 7.8 years), and provide a comprehensive comparison to individuals with other molecular mechanisms that result in SAS (non-ΔSAS). Individuals in the ΔSAS group were often underweight for age (20/41 = 49%) with a progressive decline in weight (95% CI = -2.3 to -1.1, p < 0.0001) and height (95% CI = -2.3 to -1.0, p < 0.0001) Z-score means from birth to last available measurement. ΔSAS individuals were often noted to have a broad spectrum of facial dysmorphism. A composite image of ΔSAS individuals generated by automated image analysis was distinct as compared to matched controls and non-ΔSAS individuals. We also present additional genotype-phenotype correlations for individuals in the ΔSAS group such as an increased risk for aortic root/ascending aorta dilation and primary pulmonary hypertension for those individuals with contiguous gene deletions that include COL3A1/COL5A2 and BMPR2, respectively. Based on these findings, we provide additional care recommendations for individuals with ΔSAS variants.

ZMYND11 variants are a novel cause of centrotemporal and generalised epilepsies with neurodevelopmental disorder.

ZMYND11 is the critical gene in chromosome 10p15.3 microdeletion syndrome, a syndromic cause of intellectual disability. The phenotype of ZMYND11 variants has recently been extended to autism and seizures. We expand on the epilepsy phenotype of 20 individuals with pathogenic variants in ZMYND11. We obtained clinical descriptions of 16 new and nine published individuals, plus detailed case history of two children. New individuals were identified through GeneMatcher, ClinVar and the European Network for Therapies in Rare Epilepsy (NETRE). Genetic evaluation was performed using gene panels or exome sequencing; variants were classified using American College of Medical Genetics (ACMG) criteria. Individuals with ZMYND11 associated epilepsy fell into three groups: (i) atypical benign partial epilepsy or idiopathic focal epilepsy (n = 8); (ii) generalised epilepsies/infantile epileptic encephalopathy (n = 4); (iii) unclassified (n = 8). Seizure prognosis ranged from spontaneous remission to drug resistant. Neurodevelopmental deficits were invariable. Dysmorphic features were variable. Variants were distributed across the gene and mostly de novo with no precise genotype-phenotype correlation. ZMYND11 is one of a small group of chromatin reader genes associated in the pathogenesis of epilepsy, and specifically ABPE. More detailed epilepsy descriptions of larger cohorts and functional studies might reveal genotype-phenotype correlation. The epileptogenic mechanism may be linked to interaction with histone H3.3.

Rough endoplasmic reticulum expansion: a consistent finding in a patient cohort with vascular Ehlers-Danlos Syndrome and Osteogenesis Imperfecta.

Vascular Ehlers-Danlos Syndrome (vEDS) and Osteogenesis Imperfecta (OI) are two forms of connective tissue disorders. Previously, transmission electron microscopy of skin biopsies was routinely performed on all patients who were clinically suspected to have vEDS. At present, molecular genetics using genomic DNA extracted from a blood sample is the first line investigation for these patients. However, when variants of uncertain clinical significance are identified on genetic testing and individuals do not have the classical features of OI or vEDS, additional phenotypic information obtained from a skin biopsy can be valuable for contributing to the evidence for re-classifying pathogenicity of variants.We present a cohort of six patients with molecularly confirmed vEDS and one patient with a severe form of OI, who each had expanded (or dilated), protein-filled, rough endoplasmic reticulum identified on transmission electron microscopy. The patients were identified through retrospective screening of medical records, and biopsies were taken between 1999-2016. We discuss the potential role for assessing rough endoplasmic reticulum expansion as a useful tool to allow further phenotyping of these individuals.

ZMYND11-related syndromic intellectual disability: 16 patients delineating and expanding the phenotypic spectrum.

Pathogenic variants in ZMYND11, which acts as a transcriptional repressor, have been associated with intellectual disability, behavioral abnormalities, and seizures. Only 11 affected individuals have been reported to date, and the phenotype associated with pathogenic variants in this gene have not been fully defined. Here, we present 16 additional patients with predicted pathogenic heterozygous variants in including four individuals from the same family, to further delineate and expand the genotypic and phenotypic spectrum of ZMYND11-related syndromic intellectual disability. The associated phenotype includes developmental delay, particularly affecting speech, mild-moderate intellectual disability, significant behavioral abnormalities, seizures, and hypotonia. There are subtle shared dysmorphic features, including prominent eyelashes and eyebrows, a depressed nasal bridge with bulbous nasal tip, anteverted nares, thin vermilion of the upper lip, and wide mouth. Novel features include brachydactyly and tooth enamel hypoplasia. Most identified variants are likely to result in premature truncation and/or nonsense-mediated decay. Two ZMYND11 variants located in the final exon-p.(Gln586*) (likely escaping nonsense-mediated decay) and p.(Cys574Arg)-are predicted to disrupt the MYND-type zinc-finger motif and likely interfere with binding to its interaction partners. Hence, the homogeneous phenotype likely results from a common mechanism of loss-of-function.

Delineation of phenotypes and genotypes related to cohesin structural protein RAD21.

RAD21 encodes a key component of the cohesin complex, and variants in RAD21 have been associated with Cornelia de Lange Syndrome (CdLS). Limited information on phenotypes attributable to RAD21 variants and genotype-phenotype relationships is currently published. We gathered a series of 49 individuals from 33 families with RAD21 alterations [24 different intragenic sequence variants (2 recurrent), 7 unique microdeletions], including 24 hitherto unpublished cases. We evaluated consequences of 12 intragenic variants by protein modelling and molecular dynamic studies. Full clinical information was available for 29 individuals. Their phenotype is an attenuated CdLS phenotype compared to that caused by variants in NIPBL or SMC1A for facial morphology, limb anomalies, and especially for cognition and behavior. In the 20 individuals with limited clinical information, additional phenotypes include Mungan syndrome (in patients with biallelic variants) and holoprosencephaly, with or without CdLS characteristics. We describe several additional cases with phenotypes including sclerocornea, in which involvement of the RAD21 variant is uncertain. Variants were frequently familial, and genotype-phenotype analyses demonstrated striking interfamilial and intrafamilial variability. Careful phenotyping is essential in interpreting consequences of RAD21 variants, and protein modeling and dynamics can be helpful in determining pathogenicity. The current study should be helpful when counseling families with a RAD21 variation.

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.