Haploinsufficiency of ZFHX3, encoding a key player in neuronal development, causes syndromic intellectual disability.

Neurodevelopmental disorders (NDDs) result from impaired development and functioning of the brain. Here, we identify loss-of-function (LoF) variation in ZFHX3 as a cause for syndromic intellectual disability (ID). ZFHX3 is a zinc-finger homeodomain transcription factor involved in various biological processes, including cell differentiation and tumorigenesis. We describe 42 individuals with protein-truncating variants (PTVs) or (partial) deletions of ZFHX3, exhibiting variable intellectual disability and autism spectrum disorder, recurrent facial features, relative short stature, brachydactyly, and, rarely, cleft palate. ZFHX3 LoF associates with a specific methylation profile in whole blood extracted DNA. Nuclear abundance of ZFHX3 increases during human brain development and neuronal differentiation. ZFHX3 was found to interact with the chromatin remodeling BRG1/Brm-associated factor complex and the cleavage and polyadenylation complex, suggesting a function in chromatin remodeling and mRNA processing. Furthermore, ChIP-seq for ZFHX3 revealed that it predominantly binds promoters of genes involved in nervous system development. We conclude that loss-of-function variants in ZFHX3 are a cause of syndromic ID associating with a specific DNA methylation profile.

Transcriptional and functional consequences of alterations to MEF2C and its topological organization in neuronal models.

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization.

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy.

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

HDAC9 structural variants disrupting TWIST1 transcriptional regulation lead to craniofacial and limb malformations.

Structural variants (SVs) can affect protein-coding sequences as well as gene regulatory elements. However, SVs disrupting protein-coding sequences that also function as cis-regulatory elements remain largely uncharacterized. Here, we show that craniosynostosis patients with SVs containing the histone deacetylase 9 (HDAC9) protein-coding sequence are associated with disruption of TWIST1 regulatory elements that reside within the HDAC9 sequence. Based on SVs within the HDAC9-TWIST1 locus, we defined the 3'-HDAC9 sequence as a critical TWIST1 regulatory region, encompassing craniofacial TWIST1 enhancers and CTCF sites. Deletions of either Twist1 enhancers (eTw5-7Δ/Δ) or CTCF site (CTCF-5Δ/Δ) within the Hdac9 protein-coding sequence led to decreased Twist1 expression and altered anterior/posterior limb expression patterns of SHH pathway genes. This decreased Twist1 expression results in a smaller sized and asymmetric skull and polydactyly that resembles Twist1+/- mouse phenotype. Chromatin conformation analysis revealed that the Twist1 promoter interacts with Hdac9 sequences that encompass Twist1 enhancers and a CTCF site, and that interactions depended on the presence of both regulatory regions. Finally, a large inversion of the entire Hdac9 sequence (Hdac9 INV/+) in mice that does not disrupt Hdac9 expression but repositions Twist1 regulatory elements showed decreased Twist1 expression and led to a craniosynostosis-like phenotype and polydactyly. Thus, our study elucidates essential components of TWIST1 transcriptional machinery that reside within the HDAC9 sequence. It suggests that SVs encompassing protein-coding sequences could lead to a phenotype that is not attributed to its protein function but rather to a disruption of the transcriptional regulation of a nearby gene.

MYT1L-associated neurodevelopmental disorder: description of 40 new cases and literature review of clinical and molecular aspects.

Pathogenic variants of the myelin transcription factor-1 like (MYT1L) gene include heterozygous missense, truncating variants and 2p25.3 microdeletions and cause a syndromic neurodevelopmental disorder (OMIM#616,521). Despite enrichment in de novo mutations in several developmental disorders and autism studies, the data on clinical characteristics and genotype-phenotype correlations are scarce, with only 22 patients with single nucleotide pathogenic variants reported. We aimed to further characterize this disorder at both the clinical and molecular levels by gathering a large series of patients with MYT1L-associated neurodevelopmental disorder. We collected genetic information on 40 unreported patients with likely pathogenic/pathogenic MYT1L variants and performed a comprehensive review of published data (total = 62 patients). We confirm that the main phenotypic features of the MYT1L-related disorder are developmental delay with language delay (95%), intellectual disability (ID, 70%), overweight or obesity (58%), behavioral disorders (98%) and epilepsy (23%). We highlight novel clinical characteristics, such as learning disabilities without ID (30%) and feeding difficulties during infancy (18%). We further describe the varied dysmorphic features (67%) and present the changes in weight over time of 27 patients. We show that patients harboring highly clustered missense variants in the 2-3-ZNF domains are not clinically distinguishable from patients with truncating variants. We provide an updated overview of clinical and genetic data of the MYT1L-associated neurodevelopmental disorder, hence improving diagnosis and clinical management of these patients.

A neuronal enhancer network upstream of MEF2C is compromised in patients with Rett-like characteristics.

Mutations in myocyte enhancer factor 2C (MEF2C), an important transcription factor in neurodevelopment, are associated with a Rett-like syndrome. Structural variants (SVs) upstream of MEF2C, which do not disrupt the gene itself, have also been found in patients with a similar phenotype, suggesting that disruption of MEF2C regulatory elements can also cause a Rett-like phenotype. To characterize those elements that regulate MEF2C during neural development and that are affected by these SVs, we used genomic tools coupled with both in vitro and in vivo functional assays. Through circularized chromosome conformation capture sequencing (4C-seq) and the assay for transposase-accessible chromatin using sequencing (ATAC-seq), we revealed a complex interaction network in which the MEF2C promoter physically contacts several distal enhancers that are deleted or translocated by disease-associated SVs. A total of 16 selected candidate regulatory sequences were tested for enhancer activity in vitro, with 14 found to be functional enhancers. Further analyses of their in vivo activity in zebrafish showed that each of these enhancers has a distinct activity pattern during development, with eight enhancers displaying neuronal activity. In summary, our results disentangle a complex regulatory network governing neuronal MEF2C expression that involves multiple distal enhancers. In addition, the characterized neuronal enhancers pose as novel candidates to screen for mutations in neurodevelopmental disorders, such as Rett-like syndrome.

Interpreting the impact of noncoding structural variation in neurodevelopmental disorders.

The emergence of novel sequencing technologies has greatly improved the identification of structural variation, revealing that a human genome harbors tens of thousands of structural variants (SVs). Since these SVs primarily impact noncoding DNA sequences, the next challenge is one of interpretation, not least to improve our understanding of human disease etiology. However, this task is severely complicated by the intricacy of the gene regulatory landscapes embedded within these noncoding regions, their incomplete annotation, as well as their dependence on the three-dimensional (3D) conformation of the genome. Also in the context of neurodevelopmental disorders (NDDs), reports of putatively causal, noncoding SVs are accumulating and understanding their impact on transcriptional regulation is presenting itself as the next step toward improved genetic diagnosis.

Phenotypic spectrum of the RBM10-mediated intellectual disability and congenital malformation syndrome beyond classic TARP syndrome features.

Pathogenic variants in the RBM10 gene cause a rare X-linked disorder described as TARP (Talipes equinovarus, Atrial septal defect, Robin sequence, and Persistent left vena cava superior) syndrome. We report two novel patients with truncating RBM10 variants in view of the literature, presenting a total of 26 patients from 15 unrelated families. Our results illustrate the highly pleiotropic nature of RBM10 pathogenic variants, beyond the classic TARP syndrome features. Major clinical characteristics include severe developmental delay, failure to thrive, brain malformations, neurological symptoms, respiratory issues, and facial dysmorphism. Minor features are growth retardation, cardiac, gastrointestinal, limb, and skeletal abnormalities. Additional recurrent features include genital and renal abnormalities as well as hearing and visual impairment. Thus, RBM10 loss of function variants typically cause an intellectual disability and congenital malformation syndrome that requires assessment of multiple organ systems at diagnosis and for which provided clinical features might simplify diagnostic assessment. Furthermore, evidence for an RBM10-related genotype-phenotype correlation is emerging, which can be important for prognosis.

Expanding the molecular spectrum and the neurological phenotype related to CAMTA1 variants.

The CAMTA1-associated phenotype was initially defined in patients with intragenic deletions and duplications who showed nonprogressive congenital ataxia, with or without intellectual disability. Here, we describe 10 individuals with CAMTA1 variants: nine previously unreported (likely) pathogenic variants comprising one missense, four frameshift and four nonsense variants, and one missense variant of unknown significance. Six patients were diagnosed following whole exome sequencing and four individuals with exome-based targeted panel analysis. Most of them present with developmental delay, manifesting in speech and motor delay. Other frequent findings are hypotonia, cognitive impairment, cerebellar dysfunction, oculomotor abnormalities, and behavioral problems. Feeding problems occur more frequently than previously observed. In addition, we present a systematic review of 19 previously published individuals with causal variants, including copy number, truncating, and missense variants. We note a tendency of more severe cognitive impairment and recurrent dysmorphic features in individuals with a copy number variant. Pathogenic variants are predominantly observed in and near the N- and C- terminal functional domains. Clinical heterogeneity is observed, but 3'-terminal variants seem to associate with less pronounced cerebellar dysfunction.

Isolated and Syndromic Retinal Dystrophy Caused by Biallelic Mutations in RCBTB1, a Gene Implicated in Ubiquitination.

Inherited retinal dystrophies (iRDs) are a group of genetically and clinically heterogeneous conditions resulting from mutations in over 250 genes. Here, homozygosity mapping and whole-exome sequencing (WES) in a consanguineous family revealed a homozygous missense mutation, c.973C>T (p.His325Tyr), in RCBTB1. In affected individuals, it was found to segregate with retinitis pigmentosa (RP), goiter, primary ovarian insufficiency, and mild intellectual disability. Subsequent analysis of WES data in different cohorts uncovered four additional homozygous missense mutations in five unrelated families in whom iRD segregates with or without syndromic features. Ocular phenotypes ranged from typical RP starting in the second decade to chorioretinal dystrophy with a later age of onset. The five missense mutations affect highly conserved residues either in the sixth repeat of the RCC1 domain or in the BTB1 domain. A founder haplotype was identified for mutation c.919G>A (p.Val307Met), occurring in two families of Mediterranean origin. We showed ubiquitous mRNA expression of RCBTB1 and demonstrated predominant RCBTB1 localization in human inner retina. RCBTB1 was very recently shown to be involved in ubiquitination, more specifically as a CUL3 substrate adaptor. Therefore, the effect on different components of the CUL3 and NFE2L2 (NRF2) pathway was assessed in affected individuals' lymphocytes, revealing decreased mRNA expression of NFE2L2 and several NFE2L2 target genes. In conclusion, our study puts forward mutations in RCBTB1 as a cause of autosomal-recessive non-syndromic and syndromic iRD. Finally, our data support a role for impaired ubiquitination in the pathogenetic mechanism of RCBTB1 mutations.

De novo missense variants in RAC3 cause a novel neurodevelopmental syndrome.

PURPOSE: RAC3 is an underexamined member of the Rho GTPase gene family that is expressed in the developing brain and linked to key cellular functions. De novo missense variants in the homolog RAC1 were recently associated with developmental disorders. In the RAC subfamily, transforming missense changes at certain shared residues have been observed in human cancers and previously characterized in experimental studies. The purpose of this study was to determine whether constitutional dysregulation of RAC3 is associated with human disease. METHODS: We discovered a RAC3 variant in the index case using genome sequencing, and searched for additional variants using international data-sharing initiatives. Functional effects of the variants were assessed using a multifaceted approach generalizable to most clinical laboratory settings. RESULTS: We rapidly identified five individuals with de novo monoallelic missense variants in RAC3, including one recurrent change. Every participant had severe intellectual disability and brain malformations. In silico protein modeling, and prior in vivo and in situ experiments, supported a transforming effect for each of the three different RAC3 variants. All variants were observed in databases of somatic variation in cancer. CONCLUSIONS: Missense variants in RAC3 cause a novel brain disorder, likely through a mechanism of constitutive protein activation.

Mutations affecting the SAND domain of DEAF1 cause intellectual disability with severe speech impairment and behavioral problems.

Recently, we identified in two individuals with intellectual disability (ID) different de novo mutations in DEAF1, which encodes a transcription factor with an important role in embryonic development. To ascertain whether these mutations in DEAF1 are causative for the ID phenotype, we performed targeted resequencing of DEAF1 in an additional cohort of over 2,300 individuals with unexplained ID and identified two additional individuals with de novo mutations in this gene. All four individuals had severe ID with severely affected speech development, and three showed severe behavioral problems. DEAF1 is highly expressed in the CNS, especially during early embryonic development. All four mutations were missense mutations affecting the SAND domain of DEAF1. Altered DEAF1 harboring any of the four amino acid changes showed impaired transcriptional regulation of the DEAF1 promoter. Moreover, behavioral studies in mice with a conditional knockout of Deaf1 in the brain showed memory deficits and increased anxiety-like behavior. Our results demonstrate that mutations in DEAF1 cause ID and behavioral problems, most likely as a result of impaired transcriptional regulation by DEAF1.

New insights into the phenotype of FARS2 deficiency.

Mutations in FARS2 are known to cause dysfunction of mitochondrial translation due to deficient aminoacylation of the mitochondrial phenylalanine tRNA. Here, we report three novel mutations in FARS2 found in two patients in a compound heterozygous state. The missense mutation c.1082C>T (p.Pro361Leu) was detected in both patients. The mutations c.461C>T (p.Ala154Val) and c.521_523delTGG (p.Val174del) were each detected in one patient. We report abnormal in vitro aminoacylation assays as a functional validation of the molecular genetic findings. Based on the phenotypic data of previously reported subjects and the two subjects reported here, we conclude that FARS2 deficiency can be associated with two phenotypes: (i) an epileptic phenotype, and (ii) a spastic paraplegia phenotype.

Exonic deletions in AUTS2 cause a syndromic form of intellectual disability and suggest a critical role for the C terminus.

Genomic rearrangements involving AUTS2 (7q11.22) are associated with autism and intellectual disability (ID), although evidence for causality is limited. By combining the results of diagnostic testing of 49,684 individuals, we identified 24 microdeletions that affect at least one exon of AUTS2, as well as one translocation and one inversion each with a breakpoint within the AUTS2 locus. Comparison of 17 well-characterized individuals enabled identification of a variable syndromic phenotype including ID, autism, short stature, microcephaly, cerebral palsy, and facial dysmorphisms. The dysmorphic features were more pronounced in persons with 3'AUTS2 deletions. This part of the gene is shown to encode a C-terminal isoform (with an alternative transcription start site) expressed in the human brain. Consistent with our genetic data, suppression of auts2 in zebrafish embryos caused microcephaly that could be rescued by either the full-length or the C-terminal isoform of AUTS2. Our observations demonstrate a causal role of AUTS2 in neurocognitive disorders, establish a hitherto unappreciated syndromic phenotype at this locus, and show how transcriptional complexity can underpin human pathology. The zebrafish model provides a valuable tool for investigating the etiology of AUTS2 syndrome and facilitating gene-function analysis in the future.

Implementation of non-invasive prenatal testing by semiconductor sequencing in a genetic laboratory.

OBJECTIVES: To implement non-invasive prenatal testing (NIPT) for fetal aneuploidies with semiconductor sequencing in an academic cytogenomic laboratory and to evaluate the first 15-month experience on clinical samples. METHODS: We validated a NIPT protocol for cell-free fetal DNA sequencing from maternal plasma for the detection of trisomy 13, 18 and 21 on a semiconductor sequencing instrument. Fetal DNA fraction calculation for all samples and several quality parameters were implemented in the workflow. One thousand eighty-one clinical NIPT samples were analysed, following the described protocol. RESULTS: Non-invasive prenatal testing was successfully implemented and validated on 201 normal and 74 aneuploid samples. From 1081 clinical samples, 17 samples showed an abnormal result: 14 trisomy 21 samples, one trisomy 18 and one trisomy 16 were detected. Also a maternal copy number variation on chromosome 13 was observed, which could potentially lead to a false positive trisomy 13 result. One sex discordant result was reported, possibly attributable to a vanishing twin. Moreover, our combined fetal fraction calculation enabled a more reliable risk estimate for trisomy 13, 18 and 21. CONCLUSIONS: Non-invasive prenatal testing for trisomy 21, 18 and 13 has a very high specificity and sensitivity. Because of several biological phenomena, diagnostic invasive confirmation of abnormal results remains required. © 2016 The Authors. Prenatal Diagnosis published by John Wiley & Sons, Ltd.

Two siblings with homozygous pathogenic splice-site variant in mitochondrial asparaginyl-tRNA synthetase (NARS2).

A homozygous missense mutation (c.822G>C) was found in the gene encoding the mitochondrial asparaginyl-tRNA synthetase (NARS2) in two siblings born to consanguineous parents. These siblings presented with different phenotypes: one had mild intellectual disability and epilepsy in childhood, whereas the other had severe myopathy. Biochemical analysis of the oxidative phosphorylation (OXPHOS) complexes in both siblings revealed a combined complex I and IV deficiency in skeletal muscle. In-gel activity staining after blue native-polyacrylamide gel electrophoresis confirmed the decreased activity of complex I and IV, and, in addition, showed the presence of complex V subcomplexes. Considering the consanguineous descent, homozygosity mapping and whole-exome sequencing were combined revealing the presence of one single missense mutation in the shared homozygous region. The c.822G>C variant affects the 3' splice site of exon 7, leading to skipping of the whole exon 7 and a part of exon 8 in the NARS2 mRNA. In EBV-transformed lymphoblasts, a specific decrease in the amount of charged mt-tRNA(Asn) was demonstrated as compared with controls. This confirmed the pathogenic nature of the variant. To conclude, the reported variant in NARS2 results in a combined OXPHOS complex deficiency involving complex I and IV, making NARS2 a new member of disease-associated aaRS2.

Refinement of the critical 2p25.3 deletion region: the role of MYT1L in intellectual disability and obesity.

PURPOSE: Submicroscopic deletions of chromosome band 2p25.3 are associated with intellectual disability and/or central obesity. Although MYT1L is believed to be a critical gene responsible for intellectual disability, so far no unequivocal data have confirmed this hypothesis. METHODS: In this study we evaluated a cohort of 22 patients (15 sporadic patients and two families) with a 2p25.3 aberration to further refine the clinical phenotype and to delineate the role of MYT1L in intellectual disability and obesity. In addition, myt1l spatiotemporal expression in zebrafish embryos was analyzed by quantitative polymerase chain reaction and whole-mount in situ hybridization. RESULTS: Complete MYT1L deletion, intragenic deletion, or duplication was observed in all sporadic patients, in addition to two patients with a de novo point mutation in MYT1L. The familial cases comprise a 6-Mb deletion in a father and his three children and a 5' MYT1L overlapping duplication in a father and his two children. Expression analysis in zebrafish embryos shows specific myt1l expression in the developing brain. CONCLUSION: Our data strongly strengthen the hypothesis that MYT1L is the causal gene for the observed syndromal intellectual disability. Moreover, because 17 patients present with obesity/overweight, haploinsufficiency of MYT1L might predispose to weight problems with childhood onset.Genet Med 17 6, 460-466.

Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci.

BACKGROUND: Vision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space. RESULTS: To understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4. CONCLUSIONS: Through comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development.

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.