Loss-of-function of the Zinc Finger Homeobox 4 (ZFHX4) gene underlies a neurodevelopmental disorder.

8q21.11 microdeletions encompassing the gene encoding transcription factor ZFHX4, have previously been associated by us with a syndromic form of intellectual disability, hypotonia, decreased balance and hearing loss. Here, we report on 57 individuals, 52 probands and 5 affected family members, with protein truncating variants (n=36), (micro)deletions (n=20) or an inversion (n=1) affecting ZFHX4 with variable developmental delay and intellectual disability, distinctive facial characteristics, morphological abnormalities of the central nervous system, behavioral alterations, short stature, hypotonia, and occasionally cleft palate and anterior segment dysgenesis. The phenotypes associated with 8q21.11 microdeletions and ZFHX4 intragenic loss-of-function variants largely overlap, identifying ZFHX4 as the main driver for the microdeletion syndrome, although leukocyte-derived DNA shows a mild common methylation profile for (micro)deletions only. We identify ZFHX4 as a transcription factor that is increasingly expressed during human brain development and neuronal differentiation. Furthermore, ZFHX4 interacting factors identified via IP-MS in neural progenitor cells, suggest an important role for ZFHX4 in cellular and developmental pathways, especially during histone modifications, cytosolic transport and development. Additionally, using CUT&RUN, we observed that ZFHX4 binds with the promoter regions of genes with crucial roles in embryonic, neuron and axon development. Since loss-of-function variants in ZFHX4 are found with consistent dysmorphic facial features, we investigated whether the disruption of zfhx4 causes craniofacial abnormalities in zebrafish. First-generation (F0) zfhx4 crispant zebrafish, (mosaic) mutant for zfhx4 loss-of-function variants, have significantly shorter Meckel's cartilages and smaller ethmoid plates compared to control zebrafish. Furthermore, behavioral assays show a decreased movement frequency in the zfhx4 crispant zebrafish in comparison with control zebrafish larvae. Although further research is needed, our in vivo work suggests a role for zfhx4 in facial skeleton patterning, palatal development and behavior.

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.

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.

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.

A novel neurodevelopmental syndrome caused by loss-of-function of the Zinc Finger Homeobox 3 (ZFHX3) gene.

Neurodevelopmental disorders (NDDs) result from impaired development and functioning of the brain. Here, we identify loss-of-function variation in ZFHX3 as a novel cause for syndromic intellectual disability (ID). ZFHX3, previously known as ATBF1, is a zinc-finger homeodomain transcription factor involved in multiple biological processes including cell differentiation and tumorigenesis. Through international collaboration, we collected clinical and morphometric data (Face2Gene) of 41 individuals with protein truncating variants (PTVs) or (partial) deletions of ZFHX3 . We used data mining, RNA and protein analysis to identify the subcellular localization and spatiotemporal expression of ZFHX3 in multiple in vitro models. We identified the DNA targets of ZFHX3 using ChIP seq. Immunoprecipitation followed by mass spectrometry indicated potential binding partners of endogenous ZFHX3 in neural stem cells that were subsequently confirmed by reversed co-immunoprecipitation and western blot. We evaluated a DNA methylation profile associated with ZFHX3 haploinsufficiency using DNA methylation analysis on whole blood extracted DNA of six individuals with ZFHX3 PTVs and four with a (partial) deletion of ZFHX3 . A reversed genetic approach characterized the ZFHX3 orthologue in Drosophila melanogaster . Loss-of-function variation of ZFHX3 consistently associates with (mild) ID and/or behavioural problems, postnatal growth retardation, feeding difficulties, and recognizable facial characteristics, including the rare occurrence of cleft palate. Nuclear abundance of ZFHX3 increases during human brain development and neuronal differentiation in neural stem cells and SH-SY5Y cells, ZFHX3 interacts with the chromatin remodelling BRG1/Brm-associated factor complex and the cleavage and polyadenylation complex. In line with a role for chromatin remodelling, ZFHX3 haploinsufficiency associates with a specific DNA methylation profile in leukocyte-derived DNA. The target genes of ZFHX3 are implicated in neuron and axon development. In Drosophila melanogaster , z fh2, considered to be the ZFHX3 orthologue, is expressed in the third instar larval brain. Ubiquitous and neuron-specific knockdown of zfh2 results in adult lethality underscoring a key role for zfh2 in development and neurodevelopment. Interestingly, ectopic expression of zfh2 as well as ZFHX3 in the developing wing disc results in a thoracic cleft phenotype. Collectively, our data shows that loss-of-function variants in ZFHX3 are a cause of syndromic ID, that associates with a specific DNA methylation profile. Furthermore, we show that ZFHX3 participates in chromatin remodelling and mRNA processing.

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.

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.

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.

Expanding the Phenotype of B3GALNT2-Related Disorders.

Dystroglycanopathies are a group of congenital muscular dystrophies (CMDs) that include a broad phenotypic spectrum ranging from late-onset limb-girdle muscular dystrophy to severe muscle-eye-brain disease, Walker-Warburg syndrome, and Fukuyama congenital muscular dystrophy. In addition to clinical heterogeneity, CMDs are characterized by genetic heterogeneity. To date, 18 genes have been associated with CMDs. One of them is B3GALNT2, which encodes the ß-1,3-N-acetylgalactosaminyltransferase 2 that glycosylates α-dystroglycan. In this study, using exome sequencing, we identify a homozygous frameshift variant in B3GALNT2 due to a mixed uniparental disomy of chromosome 1 in a 7-year-old girl with global developmental delay, severely delayed active language development, and autism spectrum disorder but without any symptoms of muscular dystrophy. In addition to this case, we also provide an overview of all previously reported cases, further expanding the phenotypic spectrum.

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.