Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428.

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Mono and biallelic variants in HCN2 cause severe neurodevelopmental disorders.

Hyperpolarization activated Cyclic Nucleotide (HCN) gated channels are crucial for various neurophysiological functions, including learning and sensory functions, and their dysfunction are responsible for brain disorders, such as epilepsy. To date, HCN2 variants have only been associated with mild epilepsy and recently, one monoallelic missense variant has been linked to developmental and epileptic encephalopathy. Here, we expand the phenotypic spectrum of HCN2- related disorders by describing twenty-one additional individuals from fifteen unrelated families carrying HCN2 variants. Seventeen individuals had developmental delay/intellectual disability (DD/ID), two had borderline DD/ID, and one had borderline DD. Ten individuals had epilepsy with DD/ID, with median age of onset of 10 months, and one had epilepsy with normal development. Molecular diagnosis identified thirteen different pathogenic HCN2 variants, including eleven missense variants affecting highly conserved amino acids, one frameshift variant, and one in-frame deletion. Seven variants were monoallelic of which five occurred de novo, one was not maternally inherited, one was inherited from a father with mild learning disabilities, and one was of unknown inheritance. The remaining six variants were biallelic, with four homozygous and two compound heterozygous variants. Functional studies using two-electrode voltage-clamp recordings in Xenopus laevis oocytes were performed on three monoallelic variants, p.(Arg324His), p.(Ala363Val), and p.(Met374Leu), and three biallelic variants, p.(Leu377His), p.(Pro493Leu) and p.(Gly587Asp). The p.(Arg324His) variant induced a strong increase of HCN2 conductance, while p.(Ala363Val) and p.(Met374Leu) displayed dominant negative effects, leading to a partial loss of HCN2 channel function. By confocal imaging, we found that the p.(Leu377His), p.(Pro493Leu) and p.(Gly587Asp) pathogenic variants impaired membrane trafficking, resulting in a complete loss of HCN2 elicited currents in Xenopus oocytes. Structural 3D-analysis in depolarized and hyperpolarized states of HCN2 channels, revealed that the pathogenic variants p.(His205Gln), p.(Ser409Leu), p.(Arg324Cys), p.(Asn369Ser) and p.(Gly460Asp) modify molecular interactions altering HCN2 function. Taken together, our data broadens the clinical spectrum associated with HCN2 variants, and disclose that HCN2 is involved in developmental encephalopathy with or without epilepsy.

The Healthcare and Societal Costs of Familial Intellectual Disability.

Most of the studies on the cost of intellectual disability are limited to a healthcare perspective or cohorts composed of individuals where the etiology of the condition is a mixture of genetic and non-genetic factors. When used in policy development, these can impact the decisions made on the optimal allocation of resources. In our study, we have developed a static microsimulation model to estimate the healthcare, societal, and lifetime cost of individuals with familial intellectual disability, an inheritable form of the condition, to families and government. The results from our modeling show that the societal costs outweighed the health costs (approximately 89.2% and 10.8%, respectively). The lifetime cost of familial intellectual disability is approximately AUD 7 million per person and AUD 10.8 million per household. The lifetime costs to families are second to those of the Australian Commonwealth government (AUD 4.2 million and AUD 9.3 million per household, respectively). These findings suggest that familial intellectual disability is a very expensive condition, representing a significant cost to families and government. Understanding the drivers of familial intellectual disability, especially societal, can assist us in the development of policies aimed at improving health outcomes and greater access to social care for affected individuals and their families.

Biallellic variants in CACNA1S cause fetal akinesia sequence, progressive hydrops and stillbirth.

Fetal arthrogryposis is a well-recognised ultrasonographic phenotype, caused by both genetic, maternal and extrinsic factors. When present with fetal growth restriction, pulmonary hypoplasia and multiple joint contractures, it is often referred to as fetal akinesia deformation sequence (FADS). Historically, elucidating genetic causes of arthryogryposis/FADS has been challenging; there are now more than 150 genes known to cause arthrogryposis through myopathic, neuromuscular and metabolic pathways affecting fetal movement. FADS is associated with over 400 medical conditions making prenatal diagnosis challenging. Here we present a case of FADS diagnosed at 19 weeks gestation with progression to severe fetal hydrops and stillbirth at 26-weeks gestation. Initial investigations including combined first trimester screening, TORCH (infection) screen and chromosomal microarray were normal. Trio whole exome sequencing (WES) detected compound heterozygous likely pathogenic CACNA1S gene variants associated with autosomal dominant (AD) and autosomal recessive (AR) congenital myopathy and FADS. To our knowledge, this is the first prenatal diagnosis of this condition.

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

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

Integration of EpiSign, facial phenotyping, and likelihood ratio interpretation of clinical abnormalities in the re-classification of an ARID1B missense variant.

Heterozygous ARID1B variants result in Coffin-Siris syndrome. Features may include hypoplastic nails, slow growth, characteristic facial features, hypotonia, hypertrichosis, and sparse scalp hair. Most reported cases are due to ARID1B loss of function variants. We report a boy with developmental delay, feeding difficulties, aspiration, recurrent respiratory infections, slow growth, and hypotonia without a clinical diagnosis, where a previously unreported ARID1B missense variant was classified as a variant of uncertain significance. The pathogenicity of this variant was refined through combined methodologies including genome-wide methylation signature analysis (EpiSign), Machine Learning (ML) facial phenotyping, and LIRICAL. Trio exome sequencing and EpiSign were performed. ML facial phenotyping compared facial images using FaceMatch and GestaltMatcher to syndrome-specific libraries to prioritize the trio exome bioinformatic pipeline gene list output. Phenotype-driven variant prioritization was performed with LIRICAL. A de novo heterozygous missense variant, ARID1B p.(Tyr1268His), was reported as a variant of uncertain significance. The ACMG classification was refined to likely pathogenic by a supportive methylation signature, ML facial phenotyping, and prioritization through LIRICAL. The ARID1B genotype-phenotype has been expanded through an extended analysis of missense variation through genome-wide methylation signatures, ML facial phenotyping, and likelihood-ratio gene prioritization.

Identification of Regulatory Molecular "Hot Spots" for LH/PLOD Collagen Glycosyltransferase Activity.

Hydroxylysine glycosylations are post-translational modifications (PTMs) essential for the maturation and homeostasis of fibrillar and non-fibrillar collagen molecules. The multifunctional collagen lysyl hydroxylase 3 (LH3/PLOD3) and the collagen galactosyltransferase GLT25D1 are the human enzymes that have been identified as being responsible for the glycosylation of collagen lysines, although a precise description of the contribution of each enzyme to these essential PTMs has not yet been provided in the literature. LH3/PLOD3 is thought to be capable of performing two chemically distinct collagen glycosyltransferase reactions using the same catalytic site: an inverting beta-1,O-galactosylation of hydroxylysines (Gal-T) and a retaining alpha-1,2-glucosylation of galactosyl hydroxylysines (Glc-T). In this work, we have combined indirect luminescence-based assays with direct mass spectrometry-based assays and molecular structure studies to demonstrate that LH3/PLOD3 only has Glc-T activity and that GLT25D1 only has Gal-T activity. Structure-guided mutagenesis confirmed that the Glc-T activity is defined by key residues in the first-shell environment of the glycosyltransferase catalytic site as well as by long-range contributions from residues within the same glycosyltransferase (GT) domain. By solving the molecular structures and characterizing the interactions and solving the molecular structures of human LH3/PLOD3 in complex with different UDP-sugar analogs, we show how these studies could provide insights for LH3/PLOD3 glycosyltransferase inhibitor development. Collectively, our data provide new tools for the direct investigation of collagen hydroxylysine PTMs and a comprehensive overview of the complex network of shapes, charges, and interactions that enable LH3/PLOD3 glycosyltransferase activities, expanding the molecular framework and facilitating an improved understanding and manipulation of glycosyltransferase functions in biomedical applications.

Stretch-activated ion channel TMEM63B associates with developmental and epileptic encephalopathies and progressive neurodegeneration.

By converting physical forces into electrical signals or triggering intracellular cascades, stretch-activated ion channels allow the cell to respond to osmotic and mechanical stress. Knowledge of the pathophysiological mechanisms underlying associations of stretch-activated ion channels with human disease is limited. Here, we describe 17 unrelated individuals with severe early-onset developmental and epileptic encephalopathy (DEE), intellectual disability, and severe motor and cortical visual impairment associated with progressive neurodegenerative brain changes carrying ten distinct heterozygous variants of TMEM63B, encoding for a highly conserved stretch-activated ion channel. The variants occurred de novo in 16/17 individuals for whom parental DNA was available and either missense, including the recurrent p.Val44Met in 7/17 individuals, or in-frame, all affecting conserved residues located in transmembrane regions of the protein. In 12 individuals, hematological abnormalities co-occurred, such as macrocytosis and hemolysis, requiring blood transfusions in some. We modeled six variants (p.Val44Met, p.Arg433His, p.Thr481Asn, p.Gly580Ser, p.Arg660Thr, and p.Phe697Leu), each affecting a distinct transmembrane domain of the channel, in transfected Neuro2a cells and demonstrated inward leak cation currents across the mutated channel even in isotonic conditions, while the response to hypo-osmotic challenge was impaired, as were the Ca2+ transients generated under hypo-osmotic stimulation. Ectopic expression of the p.Val44Met and p.Gly580Cys variants in Drosophila resulted in early death. TMEM63B-associated DEE represents a recognizable clinicopathological entity in which altered cation conductivity results in a severe neurological phenotype with progressive brain damage and early-onset epilepsy associated with hematological abnormalities in most individuals.

De novo MCM6 variants in neurodevelopmental disorders: a recognizable phenotype related to zinc binding residues.

The minichromosome maintenance (MCM) complex acts as a DNA helicase during DNA replication, and thereby regulates cell cycle progression and proliferation. In addition, MCM-complex components localize to centrosomes and play an independent role in ciliogenesis. Pathogenic variants in genes coding for MCM components and other DNA replication factors have been linked to growth and developmental disorders as Meier-Gorlin syndrome and Seckel syndrome. Trio exome/genome sequencing identified the same de novo MCM6 missense variant p.(Cys158Tyr) in two unrelated individuals that presented with overlapping phenotypes consisting of intra-uterine growth retardation, short stature, congenital microcephaly, endocrine features, developmental delay and urogenital anomalies. The identified variant affects a zinc binding cysteine in the MCM6 zinc finger signature. This domain, and specifically cysteine residues, are essential for MCM-complex dimerization and the induction of helicase activity, suggesting a deleterious effect of this variant on DNA replication. Fibroblasts derived from the two affected individuals showed defects both in ciliogenesis and cell proliferation. We additionally traced three unrelated individuals with de novo MCM6 variants in the oligonucleotide binding (OB)-fold domain, presenting with variable (neuro)developmental features including autism spectrum disorder, developmental delay, and epilepsy. Taken together, our findings implicate de novo MCM6 variants in neurodevelopmental disorders. The clinical features and functional defects related to the zinc binding residue resemble those observed in syndromes related to other MCM components and DNA replication factors, while de novo OB-fold domain missense variants may be associated with more variable neurodevelopmental phenotypes. These data encourage consideration of MCM6 variants in the diagnostic arsenal of NDD.

Dominant-negative variants in CBX1 cause a neurodevelopmental disorder.

PURPOSE: This study aimed to establish variants in CBX1, encoding heterochromatin protein 1ß (HP1ß), as a cause of a novel syndromic neurodevelopmental disorder. METHODS: Patients with CBX1 variants were identified, and clinician researchers were connected using GeneMatcher and physician referrals. Clinical histories were collected from each patient. To investigate the pathogenicity of identified variants, we performed in vitro cellular assays and neurobehavioral and cytological analyses of neuronal cells obtained from newly generated Cbx1 mutant mouse lines. RESULTS: In 3 unrelated individuals with developmental delay, hypotonia, and autistic features, we identified heterozygous de novo variants in CBX1. The identified variants were in the chromodomain, the functional domain of HP1ß, which mediates interactions with chromatin. Cbx1 chromodomain mutant mice displayed increased latency-to-peak response, suggesting the possibility of synaptic delay or myelination deficits. Cytological and chromatin immunoprecipitation experiments confirmed the reduction of mutant HP1ß binding to heterochromatin, whereas HP1ß interactome analysis demonstrated that the majority of HP1ß-interacting proteins remained unchanged between the wild-type and mutant HP1ß. CONCLUSION: These collective findings confirm the role of CBX1 in developmental disabilities through the disruption of HP1ß chromatin binding during neurocognitive development. Because HP1ß forms homodimers and heterodimers, mutant HP1ß likely sequesters wild-type HP1ß and other HP1 proteins, exerting dominant-negative effects.

Biallelic variants in TUBGCP6 result in microcephaly and chorioretinopathy 1: Report of four cases and a literature review.

Autosomal recessive microcephaly and chorioretinopathy-1 (MCCRP1) is a rare Mendelian disorder resulting from biallelic loss of function variants in Tubulin-Gamma Complex Associated Protein 6 (TUBGCP6, MIM#610053). Clinical features of this disorder include microcephaly, cognitive impairment, dysmorphic features, and variable ophthalmological anomalies including chorioretinopathy. Microcephaly can be recognized prenatally and visual impairment becomes evident during the first year of life. The clinical presentation resembles the findings in some acquired conditions such as congenital toxoplasmosis and cytomegalovirus infections; thus, it is important to recognize and diagnose this syndrome in view of its impact on patient health management and familial reproductive plans. To date, only seven molecularly confirmed patients from five unrelated families have been reported. We report an additional four unrelated patients with TUBGCP6 variants including one prenatal diagnosis and review the clinical phenotypes and genotypes of all the known cases. This report expands the molecular and phenotypic spectrum of TUBGCP6 and includes additional prenatal findings associated with MCCRP1.

Australian Genomics: Outcomes of a 5-year national program to accelerate the integration of genomics in healthcare.

Australian Genomics is a national collaborative partnership of more than 100 organizations piloting a whole-of-system approach to integrating genomics into healthcare, based on federation principles. In the first five years of operation, Australian Genomics has evaluated the outcomes of genomic testing in more than 5,200 individuals across 19 rare disease and cancer flagship studies. Comprehensive analyses of the health economic, policy, ethical, legal, implementation and workforce implications of incorporating genomics in the Australian context have informed evidence-based change in policy and practice, resulting in national government funding and equity of access for a range of genomic tests. Simultaneously, Australian Genomics has built national skills, infrastructure, policy, and data resources to enable effective data sharing to drive discovery research and support improvements in clinical genomic delivery.

Bi-allelic TTI1 variants cause an autosomal-recessive neurodevelopmental disorder with microcephaly.

Telomere maintenance 2 (TELO2), Tel2 interacting protein 2 (TTI2), and Tel2 interacting protein 1 (TTI1) are the three components of the conserved Triple T (TTT) complex that modulates activity of phosphatidylinositol 3-kinase-related protein kinases (PIKKs), including mTOR, ATM, and ATR, by regulating the assembly of mTOR complex 1 (mTORC1). The TTT complex is essential for the expression, maturation, and stability of ATM and ATR in response to DNA damage. TELO2- and TTI2-related bi-allelic autosomal-recessive (AR) encephalopathies have been described in individuals with moderate to severe intellectual disability (ID), short stature, postnatal microcephaly, and a movement disorder (in the case of variants within TELO2). We present clinical, genomic, and functional data from 11 individuals in 9 unrelated families with bi-allelic variants in TTI1. All present with ID, and most with microcephaly, short stature, and a movement disorder. Functional studies performed in HEK293T cell lines and fibroblasts and lymphoblastoid cells derived from 4 unrelated individuals showed impairment of the TTT complex and of mTOR pathway activity which is improved by treatment with Rapamycin. Our data delineate a TTI1-related neurodevelopmental disorder and expand the group of disorders related to the TTT complex.

A mutational hotspot in AMOTL1 defines a new syndrome of orofacial clefting, cardiac anomalies, and tall stature.

AMOTL1 encodes angiomotin-like protein 1, an actin-binding protein that regulates cell polarity, adhesion, and migration. The role of AMOTL1 in human disease is equivocal. We report a large cohort of individuals harboring heterozygous AMOTL1 variants and define a core phenotype of orofacial clefting, congenital heart disease, tall stature, auricular anomalies, and gastrointestinal manifestations in individuals with variants in AMOTL1 affecting amino acids 157-161, a functionally undefined but highly conserved region. Three individuals with AMOTL1 variants outside this region are also described who had variable presentations with orofacial clefting and multi-organ disease. Our case cohort suggests that heterozygous missense variants in AMOTL1, most commonly affecting amino acid residues 157-161, define a new orofacial clefting syndrome, and indicates an important functional role for this undefined region.

The role of exome sequencing in childhood interstitial or diffuse lung disease.

BACKGROUND: Children's interstitial and diffuse lung disease (chILD) is a complex heterogeneous group of lung disorders. Gene panel approaches have a reported diagnostic yield of ~ 12%. No data currently exist using trio exome sequencing as the standard diagnostic modality. We assessed the diagnostic utility of using trio exome sequencing in chILD. We prospectively enrolled children meeting specified clinical criteria between 2016 and 2020 from 16 Australian hospitals. Exome sequencing was performed with analysis of an initial gene panel followed by trio exome analysis. A subset of critically ill infants underwent ultra-rapid trio exome sequencing as first-line test. RESULTS: 36 patients [median (range) age 0.34 years (0.02-11.46); 11F] were recruited from multiple States and Territories. Five patients had clinically significant likely pathogenic/pathogenic variants (RARB, RPL15, CTCF, RFXANK, TBX4) and one patient had a variant of uncertain significance (VIP) suspected to contribute to their clinical phenotype, with VIP being a novel gene candidate. CONCLUSIONS: Trio exomes (6/36; 16.7%) had a better diagnostic rate than gene panel (1/36; 2.8%), due to the ability to consider a broader range of underlying conditions. However, the aetiology of chILD in most cases remained undetermined, likely reflecting the interplay between low penetrant genetic and environmental factors.

Heterozygous NOTCH1 Variants Cause CNS Immune Activation and Microangiopathy.

NOTCH1 belongs to the NOTCH family of proteins that regulate cell fate and inflammatory responses. Somatic and germline NOTCH1 variants have been implicated in cancer, Adams-Oliver syndrome, and cardiovascular defects. We describe 7 unrelated patients grouped by the presence of leukoencephalopathy with calcifications and heterozygous de novo gain-of-function variants in NOTCH1. Immunologic profiling showed upregulated CSF IP-10, a cytokine secreted downstream of NOTCH1 signaling. Autopsy revealed extensive leukoencephalopathy and microangiopathy with vascular calcifications. This evidence implicates that heterozygous gain-of-function variants in NOTCH1 lead to a chronic central nervous system (CNS) inflammatory response resulting in a calcifying microangiopathy with leukoencephalopathy. ANN NEUROL 2022;92:895-901.

Whole exome and genome sequencing in mendelian disorders: a diagnostic and health economic analysis.

Whole genome sequencing (WGS) improves Mendelian disorder diagnosis over whole exome sequencing (WES); however, additional diagnostic yields and costs remain undefined. We investigated differences between diagnostic and cost outcomes of WGS and WES in a cohort with suspected Mendelian disorders. WGS was performed in 38 WES-negative families derived from a 64 family Mendelian cohort that previously underwent WES. For new WGS diagnoses, contemporary WES reanalysis determined whether variants were diagnosable by original WES or unique to WGS. Diagnostic rates were estimated for WES and WGS to simulate outcomes if both had been applied to the 64 families. Diagnostic costs were calculated for various genomic testing scenarios. WGS diagnosed 34% (13/38) of WES-negative families. However, contemporary WES reanalysis on average 2 years later would have diagnosed 18% (7/38 families) resulting in a WGS-specific diagnostic yield of 19% (6/31 remaining families). In WES-negative families, the incremental cost per additional diagnosis using WGS following WES reanalysis was AU$36,710 (£19,407;US$23,727) and WGS alone was AU$41,916 (£22,159;US$27,093) compared to WES-reanalysis. When we simulated the use of WGS alone as an initial genomic test, the incremental cost for each additional diagnosis was AU$29,708 (£15,705;US$19,201) whereas contemporary WES followed by WGS was AU$36,710 (£19,407;US$23,727) compared to contemporary WES. Our findings confirm that WGS is the optimal genomic test choice for maximal diagnosis in Mendelian disorders. However, accepting a small reduction in diagnostic yield, WES with subsequent reanalysis confers the lowest costs. Whether WES or WGS is utilised will depend on clinical scenario and local resourcing and availability.

De Novo ZMYND8 variants result in an autosomal dominant neurodevelopmental disorder with cardiac malformations.

PURPOSE: ZMYND8 encodes a multidomain protein that serves as a central interactive hub for coordinating critical roles in transcription regulation, chromatin remodeling, regulation of super-enhancers, DNA damage response and tumor suppression. We delineate a novel neurocognitive disorder caused by variants in the ZMYND8 gene. METHODS: An international collaboration, exome sequencing, molecular modeling, yeast two-hybrid assays, analysis of available transcriptomic data and a knockdown Drosophila model were used to characterize the ZMYND8 variants. RESULTS: ZMYND8 variants were identified in 11 unrelated individuals; 10 occurred de novo and one suspected de novo; 2 were truncating, 9 were missense, of which one was recurrent. The disorder is characterized by intellectual disability with variable cardiovascular, ophthalmologic and minor skeletal anomalies. Missense variants in the PWWP domain of ZMYND8 abolish the interaction with Drebrin and missense variants in the MYND domain disrupt the interaction with GATAD2A. ZMYND8 is broadly expressed across cell types in all brain regions and shows highest expression in the early stages of brain development. Neuronal knockdown of the DrosophilaZMYND8 ortholog results in decreased habituation learning, consistent with a role in cognitive function. CONCLUSION: We present genomic and functional evidence for disruption of ZMYND8 as a novel etiology of syndromic intellectual disability.