The impact of inversions across 33,924 families with rare disease from a national genome sequencing project.

Detection of structural variants (SVs) is currently biased toward those that alter copy number. The relative contribution of inversions toward genetic disease is unclear. In this study, we analyzed genome sequencing data for 33,924 families with rare disease from the 100,000 Genomes Project. From a database hosting >500 million SVs, we focused on 351 genes where haploinsufficiency is a confirmed disease mechanism and identified 47 ultra-rare rearrangements that included an inversion (24 bp to 36.4 Mb, 20/47 de novo). Validation utilized a number of orthogonal approaches, including retrospective exome analysis. RNA-seq data supported the respective diagnoses for six participants. Phenotypic blending was apparent in four probands. Diagnostic odysseys were a common theme (>50 years for one individual), and targeted analysis for the specific gene had already been performed for 30% of these individuals but with no findings. We provide formal confirmation of a European founder origin for an intragenic MSH2 inversion. For two individuals with complex SVs involving the MECP2 mutational hotspot, ambiguous SV structures were resolved using long-read sequencing, influencing clinical interpretation. A de novo inversion of HOXD11-13 was uncovered in a family with Kantaputra-type mesomelic dysplasia. Lastly, a complex translocation disrupting APC and involving nine rearranged segments confirmed a clinical diagnosis for three family members and resolved a conundrum for a sibling with a single polyp. Overall, inversions play a small but notable role in rare disease, likely explaining the etiology in around 1/750 families across heterogeneous clinical cohorts.

Variants in ZFX are associated with an X-linked neurodevelopmental disorder with recurrent facial gestalt.

Pathogenic variants in multiple genes on the X chromosome have been implicated in syndromic and non-syndromic intellectual disability disorders. ZFX on Xp22.11 encodes a transcription factor that has been linked to diverse processes including oncogenesis and development, but germline variants have not been characterized in association with disease. Here, we present clinical and molecular characterization of 18 individuals with germline ZFX variants. Exome or genome sequencing revealed 11 variants in 18 subjects (14 males and 4 females) from 16 unrelated families. Four missense variants were identified in 11 subjects, with seven truncation variants in the remaining individuals. Clinical findings included developmental delay/intellectual disability, behavioral abnormalities, hypotonia, and congenital anomalies. Overlapping and recurrent facial features were identified in all subjects, including thickening and medial broadening of eyebrows, variations in the shape of the face, external eye abnormalities, smooth and/or long philtrum, and ear abnormalities. Hyperparathyroidism was found in four families with missense variants, and enrichment of different tumor types was observed. In molecular studies, DNA-binding domain variants elicited differential expression of a small set of target genes relative to wild-type ZFX in cultured cells, suggesting a gain or loss of transcriptional activity. Additionally, a zebrafish model of ZFX loss displayed an altered behavioral phenotype, providing additional evidence for the functional significance of ZFX. Our clinical and experimental data support that variants in ZFX are associated with an X-linked intellectual disability syndrome characterized by a recurrent facial gestalt, neurocognitive and behavioral abnormalities, and an increased risk for congenital anomalies and hyperparathyroidism.

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.

Rare de novo gain-of-function missense variants in DOT1L are associated with developmental delay and congenital anomalies.

Misregulation of histone lysine methylation is associated with several human cancers and with human developmental disorders. DOT1L is an evolutionarily conserved gene encoding a lysine methyltransferase (KMT) that methylates histone 3 lysine-79 (H3K79) and was not previously associated with a Mendelian disease in OMIM. We have identified nine unrelated individuals with seven different de novo heterozygous missense variants in DOT1L through the Undiagnosed Disease Network (UDN), the SickKids Complex Care genomics project, and GeneMatcher. All probands had some degree of global developmental delay/intellectual disability, and most had one or more major congenital anomalies. To assess the pathogenicity of the DOT1L variants, functional studies were performed in Drosophila and human cells. The fruit fly DOT1L ortholog, grappa, is expressed in most cells including neurons in the central nervous system. The identified DOT1L variants behave as gain-of-function alleles in flies and lead to increased H3K79 methylation levels in flies and human cells. Our results show that human DOT1L and fly grappa are required for proper development and that de novo heterozygous variants in DOT1L are associated with a Mendelian disease.

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.

Null and missense mutations of ERI1 cause a recessive phenotypic dichotomy in humans.

ERI1 is a 3'-to-5' exoribonuclease involved in RNA metabolic pathways including 5.8S rRNA processing and turnover of histone mRNAs. Its biological and medical significance remain unclear. Here, we uncover a phenotypic dichotomy associated with bi-allelic ERI1 variants by reporting eight affected individuals from seven unrelated families. A severe spondyloepimetaphyseal dysplasia (SEMD) was identified in five affected individuals with missense variants but not in those with bi-allelic null variants, who showed mild intellectual disability and digital anomalies. The ERI1 missense variants cause a loss of the exoribonuclease activity, leading to defective trimming of the 5.8S rRNA 3' end and a decreased degradation of replication-dependent histone mRNAs. Affected-individual-derived induced pluripotent stem cells (iPSCs) showed impaired in vitro chondrogenesis with downregulation of genes regulating skeletal patterning. Our study establishes an entity previously unreported in OMIM and provides a model showing a more severe effect of missense alleles than null alleles within recessive genotypes, suggesting a key role of ERI1-mediated RNA metabolism in human skeletal patterning and chondrogenesis.

Rare EIF4A2 variants are associated with a neurodevelopmental disorder characterized by intellectual disability, hypotonia, and epilepsy.

Eukaryotic initiation factor-4A2 (EIF4A2) is an ATP-dependent RNA helicase and a member of the DEAD-box protein family that recognizes the 5' cap structure of mRNAs, allows mRNA to bind to the ribosome, and plays an important role in microRNA-regulated gene repression. Here, we report on 15 individuals from 14 families presenting with global developmental delay, intellectual disability, hypotonia, epilepsy, and structural brain anomalies, all of whom have extremely rare de novo mono-allelic or inherited bi-allelic variants in EIF4A2. Neurodegeneration was predominantly reported in individuals with bi-allelic variants. Molecular modeling predicts these variants would perturb structural interactions in key protein domains. To determine the pathogenicity of the EIF4A2 variants in vivo, we examined the mono-allelic variants in Drosophila melanogaster (fruit fly) and identified variant-specific behavioral and developmental defects. The fruit fly homolog of EIF4A2 is eIF4A, a negative regulator of decapentaplegic (dpp) signaling that regulates embryo patterning, eye and wing morphogenesis, and stem cell identity determination. Our loss-of-function (LOF) rescue assay demonstrated a pupal lethality phenotype induced by loss of eIF4A, which was fully rescued with human EIF4A2 wild-type (WT) cDNA expression. In comparison, the EIF4A2 variant cDNAs failed or incompletely rescued the lethality. Overall, our findings reveal that EIF4A2 variants cause a genetic neurodevelopmental syndrome with both LOF and gain of function as underlying mechanisms.

The clinical and molecular spectrum of the KDM6B-related neurodevelopmental disorder.

De novo variants are a leading cause of neurodevelopmental disorders (NDDs), but because every monogenic NDD is different and usually extremely rare, it remains a major challenge to understand the complete phenotype and genotype spectrum of any morbid gene. According to OMIM, heterozygous variants in KDM6B cause "neurodevelopmental disorder with coarse facies and mild distal skeletal abnormalities." Here, by examining the molecular and clinical spectrum of 85 reported individuals with mostly de novo (likely) pathogenic KDM6B variants, we demonstrate that this description is inaccurate and potentially misleading. Cognitive deficits are seen consistently in all individuals, but the overall phenotype is highly variable. Notably, coarse facies and distal skeletal anomalies, as defined by OMIM, are rare in this expanded cohort while other features are unexpectedly common (e.g., hypotonia, psychosis, etc.). Using 3D protein structure analysis and an innovative dual Drosophila gain-of-function assay, we demonstrated a disruptive effect of 11 missense/in-frame indels located in or near the enzymatic JmJC or Zn-containing domain of KDM6B. Consistent with the role of KDM6B in human cognition, we demonstrated a role for the Drosophila KDM6B ortholog in memory and behavior. Taken together, we accurately define the broad clinical spectrum of the KDM6B-related NDD, introduce an innovative functional testing paradigm for the assessment of KDM6B variants, and demonstrate a conserved role for KDM6B in cognition and behavior. Our study demonstrates the critical importance of international collaboration, sharing of clinical data, and rigorous functional analysis of genetic variants to ensure correct disease diagnosis for rare disorders.

Bi-allelic loss-of-function variants in TMEM147 cause moderate to profound intellectual disability with facial dysmorphism and pseudo-Pelger-Huët anomaly.

The transmembrane protein TMEM147 has a dual function: first at the nuclear envelope, where it anchors lamin B receptor (LBR) to the inner membrane, and second at the endoplasmic reticulum (ER), where it facilitates the translation of nascent polypeptides within the ribosome-bound TMCO1 translocon complex. Through international data sharing, we identified 23 individuals from 15 unrelated families with bi-allelic TMEM147 loss-of-function variants, including splice-site, nonsense, frameshift, and missense variants. These affected children displayed congruent clinical features including coarse facies, developmental delay, intellectual disability, and behavioral problems. In silico structural analyses predicted disruptive consequences of the identified amino acid substitutions on translocon complex assembly and/or function, and in vitro analyses documented accelerated protein degradation via the autophagy-lysosomal-mediated pathway. Furthermore, TMEM147-deficient cells showed CKAP4 (CLIMP-63) and RTN4 (NOGO) upregulation with a concomitant reorientation of the ER, which was also witnessed in primary fibroblast cell culture. LBR mislocalization and nuclear segmentation was observed in primary fibroblast cells. Abnormal nuclear segmentation and chromatin compaction were also observed in approximately 20% of neutrophils, indicating the presence of a pseudo-Pelger-Huët anomaly. Finally, co-expression analysis revealed significant correlation with neurodevelopmental genes in the brain, further supporting a role of TMEM147 in neurodevelopment. Our findings provide clinical, genetic, and functional evidence that bi-allelic loss-of-function variants in TMEM147 cause syndromic intellectual disability due to ER-translocon and nuclear organization dysfunction.

Systematic reanalysis of copy number losses of uncertain clinical significance.

BACKGROUND: Reanalysis of exome/genome data improves diagnostic yield. However, the value of reanalysis of clinical array comparative genomic hybridisation (aCGH) data has never been investigated. Case-by-case reanalysis can be challenging in busy diagnostic laboratories. METHODS AND RESULTS: We harmonised historical postnatal clinical aCGH results from ~16 000 patients tested via our diagnostic laboratory over ~7 years with current clinical guidance. This led to identification of 37 009 copy number losses (CNLs) including 33 857 benign, 2173 of uncertain significance and 979 pathogenic. We found benign CNLs to be significantly less likely to encompass haploinsufficient genes compared with the pathogenic or CNLs of uncertain significance in our database. Based on this observation, we developed a reanalysis pipeline using up-to-date disease association data and haploinsufficiency scores and shortlisted 207 CNLs of uncertain significance encompassing at least one autosomal dominant disease-gene associated with haploinsufficiency or loss-of-function mechanism. Clinical scientist reviews led to reclassification of 15 CNLs of uncertain significance as pathogenic or likely pathogenic. This was ~0.7% of the starting cohort of 2173 CNLs of uncertain significance and 7.2% of 207 shortlisted CNLs. The reclassified CNLs included first cases of CNV-mediated disease for some genes where all previously described cases involved only point variants. Interestingly, some CNLs could not be reclassified because the phenotypes of patients with CNLs seemed distinct from the known clinical features resulting from point variants, thus raising questions about accepted underlying disease mechanisms. CONCLUSIONS: Reanalysis of clinical aCGH data increases diagnostic yield.

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.

Bi-allelic loss-of-function variants in BCAS3 cause a syndromic neurodevelopmental disorder.

BCAS3 microtubule-associated cell migration factor (BCAS3) is a large, highly conserved cytoskeletal protein previously proposed to be critical in angiogenesis and implicated in human embryogenesis and tumorigenesis. Here, we established BCAS3 loss-of-function variants as causative for a neurodevelopmental disorder. We report 15 individuals from eight unrelated families with germline bi-allelic loss-of-function variants in BCAS3. All probands share a global developmental delay accompanied by pyramidal tract involvement, microcephaly, short stature, strabismus, dysmorphic facial features, and seizures. The human phenotype is less severe compared with the Bcas3 knockout mouse model and cannot be explained by angiogenic defects alone. Consistent with being loss-of-function alleles, we observed absence of BCAS3 in probands' primary fibroblasts. By comparing the transcriptomic and proteomic data based on probands' fibroblasts with those of the knockout mouse model, we identified similar dysregulated pathways resulting from over-representation analysis, while the dysregulation of some proposed key interactors could not be confirmed. Together with the results from a tissue-specific Drosophila loss-of-function model, we demonstrate a vital role for BCAS3 in neural tissue development.

Mutation-specific pathophysiological mechanisms define different neurodevelopmental disorders associated with SATB1 dysfunction.

Whereas large-scale statistical analyses can robustly identify disease-gene relationships, they do not accurately capture genotype-phenotype correlations or disease mechanisms. We use multiple lines of independent evidence to show that different variant types in a single gene, SATB1, cause clinically overlapping but distinct neurodevelopmental disorders. Clinical evaluation of 42 individuals carrying SATB1 variants identified overt genotype-phenotype relationships, associated with different pathophysiological mechanisms, established by functional assays. Missense variants in the CUT1 and CUT2 DNA-binding domains result in stronger chromatin binding, increased transcriptional repression, and a severe phenotype. In contrast, variants predicted to result in haploinsufficiency are associated with a milder clinical presentation. A similarly mild phenotype is observed for individuals with premature protein truncating variants that escape nonsense-mediated decay, which are transcriptionally active but mislocalized in the cell. Our results suggest that in-depth mutation-specific genotype-phenotype studies are essential to capture full disease complexity and to explain phenotypic variability.

Discovery of DNA methylation signature in the peripheral blood of individuals with history of antenatal exposure to valproic acid.

PURPOSE: Valproic acid or valproate is an effective antiepileptic drug; however, embryonic exposure to valproate can result in a teratogenic disorder referred to as fetal valproate syndrome (OMIM #609442). Currently there are no diagnostic biomarkers for the condition. This study aims to define an episignature biomarker for teratogenic antenatal exposure to valproate. METHODS: DNA extracted from peripheral blood of individuals with teratogenic antenatal exposure to valproate was processed using DNA methylation microarrays. Subsequently, methylation profiling and construction of support vector machine classifiers were performed in R. RESULTS: Genomic DNA methylation analysis was applied, and a distinct DNA methylation profile was identified in the majority of affected individuals. This profile was used to develop a diagnostic episignature classifier. The valproate exposure episignature exhibited high sensitivity and specificity relative to a large reference data set of unaffected controls and individuals with a wide spectrum of syndromic disorders with episignatures. Gene set enrichment analysis demonstrated an enrichment for terms associated with cell adhesion, including significant overrepresentation of the cadherin superfamily. CONCLUSION: This study provides evidence of a robust peripheral blood-based diagnostic epigenetic biomarker for a prenatal teratogenic disorder.

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

PURPOSE: This study aims to assess the diagnostic utility and provide reporting recommendations for clinical DNA methylation episignature testing based on the cohort of patients tested through the EpiSign Clinical Testing Network. METHODS: The EpiSign assay utilized unsupervised clustering techniques and a support vector machine-based classification algorithm to compare each patient's genome-wide DNA methylation profile with the EpiSign Knowledge Database, yielding the result that was reported. An international working group, representing distinct EpiSign Clinical Testing Network health jurisdictions, collaborated to establish recommendations for interpretation and reporting of episignature testing. RESULTS: Among 2399 cases analyzed, 1667 cases underwent a comprehensive screen of validated episignatures, imprinting, and promoter regions, resulting in 18.7% (312/1667) positive reports. The remaining 732 referrals underwent targeted episignature analysis for assessment of sequence or copy-number variants (CNVs) of uncertain significance or for assessment of clinical diagnoses without confirmed molecular findings, and 32.4% (237/732) were positive. Cases with detailed clinical information were highlighted to describe various utility scenarios for episignature testing. CONCLUSION: Clinical DNA methylation testing including episignatures, imprinting, and promoter analysis provided by an integrated network of clinical laboratories enables test standardization and demonstrates significant diagnostic yield and clinical utility beyond DNA sequence analysis in rare diseases.

Biallelic CRELD1 variants cause a multisystem syndrome, including neurodevelopmental phenotypes, cardiac dysrhythmias, and frequent infections.

PURPOSE: We sought to delineate a multisystem disorder caused by recessive cysteine-rich with epidermal growth factor-like domains 1 (CRELD1) gene variants. METHODS: The impact of CRELD1 variants was characterized through an international collaboration utilizing next-generation DNA sequencing, gene knockdown, and protein overexpression in Xenopus tropicalis, and in vitro analysis of patient immune cells. RESULTS: Biallelic variants in CRELD1 were found in 18 participants from 14 families. Affected individuals displayed an array of phenotypes involving developmental delay, early-onset epilepsy, and hypotonia, with about half demonstrating cardiac arrhythmias and some experiencing recurrent infections. Most harbored a frameshift in trans with a missense allele, with 1 recurrent variant, p.(Cys192Tyr), identified in 10 families. X tropicalis tadpoles with creld1 knockdown displayed developmental defects along with increased susceptibility to induced seizures compared with controls. Additionally, human CRELD1 harboring missense variants from affected individuals had reduced protein function, indicated by a diminished ability to induce craniofacial defects when overexpressed in X tropicalis. Finally, baseline analyses of peripheral blood mononuclear cells showed similar proportions of immune cell subtypes in patients compared with healthy donors. CONCLUSION: This patient cohort, combined with experimental data, provide evidence of a multisystem clinical syndrome mediated by recessive variants in CRELD1.

Dominant Mutations in the Autoimmune Regulator AIRE Are Associated with Common Organ-Specific Autoimmune Diseases.

The autoimmune regulator (AIRE) gene is crucial for establishing central immunological tolerance and preventing autoimmunity. Mutations in AIRE cause a rare autosomal-recessive disease, autoimmune polyendocrine syndrome type 1 (APS-1), distinguished by multi-organ autoimmunity. We have identified multiple cases and families with mono-allelic mutations in the first plant homeodomain (PHD1) zinc finger of AIRE that followed dominant inheritance, typically characterized by later onset, milder phenotypes, and reduced penetrance compared to classical APS-1. These missense PHD1 mutations suppressed gene expression driven by wild-type AIRE in a dominant-negative manner, unlike CARD or truncated AIRE mutants that lacked such dominant capacity. Exome array analysis revealed that the PHD1 dominant mutants were found with relatively high frequency (>0.0008) in mixed populations. Our results provide insight into the molecular action of AIRE and demonstrate that disease-causing mutations in the AIRE locus are more common than previously appreciated and cause more variable autoimmune phenotypes.

Homozygous loss-of-function variants in FILIP1 cause autosomal recessive arthrogryposis multiplex congenita with microcephaly.

Arthrogryposis multiplex congenita forms a broad group of clinically and etiologically heterogeneous disorders characterized by congenital joint contractures that involve at least two different parts of the body. Neurological and muscular disorders are commonly underlying arthrogryposis. Here, we report five affected individuals from three independent families sharing an overlapping phenotype with congenital contractures affecting shoulder, elbow, hand, hip, knee and foot as well as scoliosis, reduced palmar and plantar skin folds, microcephaly and facial dysmorphism. Using exome sequencing, we identified homozygous truncating variants in FILIP1 in all patients. FILIP1 is a regulator of filamin homeostasis required for the initiation of cortical cell migration in the developing neocortex and essential for the differentiation process of cross-striated muscle cells during myogenesis. In summary, our data indicate that bi-allelic truncating variants in FILIP1 are causative of a novel autosomal recessive disorder and expand the spectrum of genetic factors causative of arthrogryposis multiplex congenita.

De novo mutations in the BMP signaling pathway in lambdoid craniosynostosis.

Lambdoid craniosynostosis (CS) is a congenital anomaly resulting from premature fusion of the cranial suture between the parietal and occipital bones. Predominantly sporadic, it is the rarest form of CS and its genetic etiology is largely unexplored. Exome sequencing of 25 kindreds, including 18 parent-offspring trios with sporadic lambdoid CS, revealed a marked excess of damaging (predominantly missense) de novo mutations that account for ~ 40% of sporadic cases. These mutations clustered in the BMP signaling cascade (P = 1.6 × 10-7), including mutations in genes encoding BMP receptors (ACVRL1 and ACVR2A), transcription factors (SOX11, FOXO1) and a transcriptional co-repressor (IFRD1), none of which have been implicated in other forms of CS. These missense mutations are at residues critical for substrate or target sequence recognition and many are inferred to cause genetic gain-of-function. Additionally, mutations in transcription factor NFIX were implicated in syndromic craniosynostosis affecting diverse sutures. Single cell RNA sequencing analysis of the mouse lambdoid suture identified enrichment of mutations in osteoblast precursors (P = 1.6 × 10-6), implicating perturbations in the balance between proliferation and differentiation of osteoprogenitor cells in lambdoid CS. The results contribute to the growing knowledge of the genetics of CS, have implications for genetic counseling, and further elucidate the molecular etiology of premature suture fusion.

Delineation of a Human Mendelian Disorder of the DNA Demethylation Machinery: TET3 Deficiency.

Germline pathogenic variants in chromatin-modifying enzymes are a common cause of pediatric developmental disorders. These enzymes catalyze reactions that regulate epigenetic inheritance via histone post-translational modifications and DNA methylation. Cytosine methylation (5-methylcytosine [5mC]) of DNA is the quintessential epigenetic mark, yet no human Mendelian disorder of DNA demethylation has yet been delineated. Here, we describe in detail a Mendelian disorder caused by the disruption of DNA demethylation. TET3 is a methylcytosine dioxygenase that initiates DNA demethylation during early zygote formation, embryogenesis, and neuronal differentiation and is intolerant to haploinsufficiency in mice and humans. We identify and characterize 11 cases of human TET3 deficiency in eight families with the common phenotypic features of intellectual disability and/or global developmental delay; hypotonia; autistic traits; movement disorders; growth abnormalities; and facial dysmorphism. Mono-allelic frameshift and nonsense variants in TET3 occur throughout the coding region. Mono-allelic and bi-allelic missense variants localize to conserved residues; all but one such variant occur within the catalytic domain, and most display hypomorphic function in an assay of catalytic activity. TET3 deficiency and other Mendelian disorders of the epigenetic machinery show substantial phenotypic overlap, including features of intellectual disability and abnormal growth, underscoring shared disease mechanisms.