IHH enhancer variant within neighboring NHEJ1 intron causes microphthalmia anophthalmia and coloboma.

Genomic sequences residing within introns of few genes have been shown to act as enhancers affecting expression of neighboring genes. We studied an autosomal recessive phenotypic continuum of microphthalmia, anophthalmia and ocular coloboma, with no apparent coding-region disease-causing mutation. Homozygosity mapping of several affected Jewish Iranian families, combined with whole genome sequence analysis, identified a 0.5 Mb disease-associated chromosome 2q35 locus (maximal LOD score 6.8) harboring an intronic founder variant in NHEJ1, not predicted to affect NHEJ1. The human NHEJ1 intronic variant lies within a known specifically limb-development enhancer of a neighboring gene, Indian hedgehog (Ihh), known to be involved in eye development in mice and chickens. Through mouse and chicken molecular development studies, we demonstrated that this variant is within an Ihh enhancer that drives gene expression in the developing eye and that the identified variant affects this eye-specific enhancer activity. We thus delineate an Ihh enhancer active in mammalian eye development whose variant causes human microphthalmia, anophthalmia and ocular coloboma. The findings highlight disease causation by an intronic variant affecting the expression of a neighboring gene, delineating molecular pathways of eye development.

Deciphering transcription factors and their corresponding regulatory elements during inhibitory interneuron differentiation using deep neural networks.

During neurogenesis, the generation and differentiation of neuronal progenitors into inhibitory gamma-aminobutyric acid-containing interneurons is dependent on the combinatorial activity of transcription factors (TFs) and their corresponding regulatory elements (REs). However, the roles of neuronal TFs and their target REs in inhibitory interneuron progenitors are not fully elucidated. Here, we developed a deep-learning-based framework to identify enriched TF motifs in gene REs (eMotif-RE), such as poised/repressed enhancers and putative silencers. Using epigenetic datasets (e.g., ATAC-seq and H3K27ac/me3 ChIP-seq) from cultured interneuron-like progenitors, we distinguished between active enhancer sequences (open chromatin with H3K27ac) and non-active enhancer sequences (open chromatin without H3K27ac). Using our eMotif-RE framework, we discovered enriched motifs of TFs such as ASCL1, SOX4, and SOX11 in the active enhancer set suggesting a cooperativity function for ASCL1 and SOX4/11 in active enhancers of neuronal progenitors. In addition, we found enriched ZEB1 and CTCF motifs in the non-active set. Using an in vivo enhancer assay, we showed that most of the tested putative REs from the non-active enhancer set have no enhancer activity. Two of the eight REs (25%) showed function as poised enhancers in the neuronal system. Moreover, mutated REs for ZEB1 and CTCF motifs increased their in vivo activity as enhancers indicating a repressive effect of ZEB1 and CTCF on these REs that likely function as repressed enhancers or silencers. Overall, our work integrates a novel framework based on deep learning together with a functional assay that elucidated novel functions of TFs and their corresponding REs. Our approach can be applied to better understand gene regulation not only in inhibitory interneuron differentiation but in other tissue and cell types.

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.

MusMorph, a database of standardized mouse morphology data for morphometric meta-analyses.

Complex morphological traits are the product of many genes with transient or lasting developmental effects that interact in anatomical context. Mouse models are a key resource for disentangling such effects, because they offer myriad tools for manipulating the genome in a controlled environment. Unfortunately, phenotypic data are often obtained using laboratory-specific protocols, resulting in self-contained datasets that are difficult to relate to one another for larger scale analyses. To enable meta-analyses of morphological variation, particularly in the craniofacial complex and brain, we created MusMorph, a database of standardized mouse morphology data spanning numerous genotypes and developmental stages, including E10.5, E11.5, E14.5, E15.5, E18.5, and adulthood. To standardize data collection, we implemented an atlas-based phenotyping pipeline that combines techniques from image registration, deep learning, and morphometrics. Alongside stage-specific atlases, we provide aligned micro-computed tomography images, dense anatomical landmarks, and segmentations (if available) for each specimen (N = 10,056). Our workflow is open-source to encourage transparency and reproducible data collection. The MusMorph data and scripts are available on FaceBase ( www.facebase.org , https://doi.org/10.25550/3-HXMC ) and GitHub ( https://github.com/jaydevine/MusMorph ).

Control of meiotic chromosomal bouquet and germ cell morphogenesis by the zygotene cilium.

A hallmark of meiosis is chromosomal pairing, which requires telomere tethering and rotation on the nuclear envelope through microtubules, driving chromosome homology searches. Telomere pulling toward the centrosome forms the "zygotene chromosomal bouquet." Here, we identified the "zygotene cilium" in oocytes. This cilium provides a cable system for the bouquet machinery and extends throughout the germline cyst. Using zebrafish mutants and live manipulations, we demonstrate that the cilium anchors the centrosome to counterbalance telomere pulling. The cilium is essential for bouquet and synaptonemal complex formation, oogenesis, ovarian development, and fertility. Thus, a cilium represents a conserved player in zebrafish and mouse meiosis, which sheds light on reproductive aspects in ciliopathies and suggests that cilia can control chromosomal dynamics.

Cytokinetic abscission is part of the midblastula transition in early zebrafish embryogenesis.

Animal cytokinesis ends with the formation of a thin intercellular membrane bridge that connects the two newly formed sibling cells, which is ultimately resolved by abscission. While mitosis is completed within 15 min, the intercellular bridge can persist for hours, maintaining a physical connection between sibling cells and allowing exchange of cytosolic components. Although cell-cell communication is fundamental for development, the role of intercellular bridges during embryogenesis has not been fully elucidated. In this work, we characterized the spatiotemporal characteristics of the intercellular bridge during early zebrafish development. We found that abscission is delayed during the rapid division cycles that occur in the early embryo, giving rise to the formation of interconnected cell clusters. Abscission was accelerated when the embryo entered the midblastula transition (MBT) phase. Components of the ESCRT machinery, which drives abscission, were enriched at intercellular bridges post-MBT and, interfering with ESCRT function, extended abscission beyond MBT. Hallmark features of MBT, including transcription onset and cell shape modulations, were more similar in interconnected sibling cells compared to other neighboring cells. Collectively, our findings suggest that delayed abscission in the early embryo allows clusters of cells to coordinate their behavior during embryonic development.

PICOT binding to chromatin-associated EED negatively regulates cyclin D2 expression by increasing H3K27me3 at the CCND2 gene promoter.

Protein kinase C (PKC)-interacting cousin of thioredoxin (PICOT; also termed glutaredoxin 3 (Grx3; Glrx3)) is a ubiquitous protein that can interact with the embryonic ectoderm development (EED) protein via each of its two C-terminal PICOT/Grx homology domains. Since EED is a Polycomb-Group protein and a core component of the polycomb repressive complex 2 (PRC2), we tested the involvement of PICOT in the regulation of PRC2-mediated H3 lysine 27 trimethylation (H3K27me3), transcription and translation of selected PRC2 target genes. A fraction of the cellular PICOT protein was found in the nuclei of leukemia cell lines, where it was associated with the chromatin. In addition, PICOT coimmunoprecipitated with chromatin-residing EED derived from Jurkat and COS-7 cell nuclei. PICOT knockdown led to a reduced H3K27me3 mark and a decrease in EED and EZH2 at the CCND2 gene promoter. In agreement, PICOT-deficient T cells exhibited a significant increase in CCND2 mRNA and protein expression. Since elevated expression levels of PICOT were reported in several different tumors and correlated in the current studies with decreased transcription and translation of the CCND2 gene, we tested whether this opposite correlation exists in human cancers. Data from the Cancer Genome Atlas (TCGA) database indicated statistically significant negative correlation between PICOT and CCND2 in eight different human tumors where the highest correlation was in lung (p = 8.67E-10) and pancreatic (p = 1.06E-5) adenocarcinoma. Furthermore, high expression of PICOT and low expression of CCND2 correlated with poor patient survival in five different types of human tumors. The results suggest that PICOT binding to chromatin-associated EED modulates the H3K27me3 level at the CCND2 gene promoter which may be one of the potential mechanisms for regulation of cyclin D2 expression in tumors. These findings also indicate that a low PICOT/CCND2 expression ratio might serve as a good predictor of patient survival in selected human cancers.

Mutations in the microtubule-associated protein MAP11 (C7orf43) cause microcephaly in humans and zebrafish.

Microtubule associated protein 11 (MAP11, previously termed C7orf43) encodes a highly conserved protein whose function is unknown. Through genome-wide linkage analysis combined with whole exome sequencing, we demonstrate that human autosomal recessive primary microcephaly is caused by a truncating mutation in MAP11. Moreover, homozygous MAP11-orthologue CRISPR/Cas9 knock-out zebrafish presented with microcephaly and decreased neuronal proliferation, recapitulating the human phenotype. We demonstrate that MAP11 is ubiquitously transcribed with high levels in brain and cerebellum. Immunofluorescence and co-immunoprecipitation studies in SH-SY5Y cells showed that MAP11 associates with mitotic spindles, co-localizing and physically associating with α-tubulin during mitosis. MAP11 expression precedes α-tubulin in gap formation of cell abscission at the midbody and is co-localized with PLK1, a key regulator of cytokinesis, at the edges of microtubule extensions of daughter cells post cytokinesis abscission, implicating a role in mitotic spindle dynamics and in regulation of cell abscission during cytokinesis. Finally, lentiviral-mediated silencing of MAP11 diminished SH-SY5Y cell viability, reducing proliferation rather than affecting apoptosis. Thus, MAP11 encodes a microtubule-associated protein that plays a role in spindle dynamics and cell division, in which mutations cause microcephaly in humans and zebrafish.

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.

Functional characterization of the ZEB2 regulatory landscape.

Zinc finger E-box-binding homeobox 2 (ZEB2) is a key developmental regulator of the central nervous system (CNS). Although the transcriptional regulation of ZEB2 is essential for CNS development, the elements that regulate ZEB2 expression have yet to be identified. Here, we identified a proximal regulatory region of ZEB2 and characterized transcriptional enhancers during neuronal development. Using chromatin immunoprecipitation sequencing for active (H3K27ac) and repressed (H3K27me3) chromatin regions in human neuronal progenitors, combined with an in vivo zebrafish enhancer assay, we functionally characterized 18 candidate enhancers in the ZEB2 locus. Eight enhancers drove expression patterns that were specific to distinct mid/hindbrain regions (ZEB2#e3 and 5), trigeminal-like ganglia (ZEB2#e6 and 7), notochord (ZEB2#e2, 4 and 12) and whole brain (ZEB2#e14). We further dissected the minimal sequences that drive enhancer-specific activity in the mid/hindbrain and notochord. Using a reporter assay in human cells, we showed an increased activity of the minimal notochord enhancer ZEB2#e2 in response to AP-1 and DLX1/2 expressions, while repressed activity of this enhancer was seen in response to ZEB2 and TFAP2 expressions. We showed that Dlx1 but not Zeb2 and Tfap2 occupies Zeb2#e2 enhancer sequence in the mouse notochord at embryonic day 11.5. Using CRISPR/Cas9 genome editing, we deleted the ZEB2#e2 region, leading to reduction of ZEB2 expression in human cells. We thus characterized distal transcriptional enhancers and trans-acting elements that govern regulation of ZEB2 expression during neuronal development. These findings pave the path toward future analysis of the role of ZEB2 regulatory elements in neurodevelopmental disorders, such as Mowat-Wilson syndrome.

Unraveling the transcriptional regulation of TWIST1 in limb development.

The transcription factor TWIST1 plays a vital role in mesoderm development, particularly in limb and craniofacial formation. Accordingly, haploinsufficiency of TWIST1 can cause limb and craniofacial malformations as part of Saethre-Chotzen syndrome. However, the molecular basis of TWIST1 transcriptional regulation during development has yet to be elucidated. Here, we characterized active enhancers in the TWIST1-HDAC9 locus that drive transcription in the developing limb and branchial arches. Using available p300 and H3K27ac ChIP-seq data, we identified 12 enhancer candidates, located both within and outside the coding sequences of the neighboring gene, Histone deacetyase 9 (HDAC9). Using zebrafish and mouse enhancer assays, we showed that eight of these candidates have limb/fin and branchial arch enhancer activity that resemble Twist1 expression. Using 4C-seq, we showed that the Twist1 promoter region interacts with three enhancers (eTw-5, 6, 7) in the limb bud and branchial arch of mouse embryos at day 11.5. Furthermore, we found that two transcription factors, LMX1B and TFAP2, bind these enhancers and modulate their enhancer activity. Finally, using CRISPR/Cas9 genome editing, we showed that homozygous deletion of eTw5-7 enhancers reduced Twist1 expression in the limb bud and caused pre-axial polydactyly, a phenotype observed in Twist1+/- mice. Taken together, our findings reveal that each enhancer has a discrete activity pattern, and together comprise a spatiotemporal regulatory network of Twist1 transcription in the developing limbs/fins and branchial arches. Our study suggests that mutations in TWIST1 enhancers could lead to reduced TWIST1 expression, resulting in phenotypic outcome as seen with TWIST1 coding mutations.

VPS4 is a dynamic component of the centrosome that regulates centrosome localization of γ-tubulin, centriolar satellite stability and ciliogenesis.

The hexameric AAA ATPase VPS4 facilitates ESCRT III filament disassembly on diverse intracellular membranes. ESCRT III components and VPS4 have been localized to the ciliary transition zone and spindle poles and reported to affect centrosome duplication and spindle pole stability. How the canonical ESCRT pathway could mediate these events is unclear. We studied the association of VPS4 with centrosomes and found that GFP-VPS4 was a dynamic component of both mother and daughter centrioles. A mutant, VPS4EQ, which can't hydrolyze ATP, was less dynamic and accumulated at centrosomes. Centrosome localization of the VPS4EQ mutant, caused reduced γ-tubulin levels at centrosomes and consequently decreased microtubule growth and altered centrosome positioning. In addition, preventing VPS4 ATP hydrolysis nearly eliminated centriolar satellites and paused ciliogensis after formation of the ciliary vesicle. Zebrafish embryos injected with GFP-VPS4EQ mRNA were less viable, exhibited developmental defects and had fewer cilia in Kupffer's vesicle. Surprisingly, ESCRT III proteins seldom localized to centrosomes and their depletion did not lead to these phenotypes. Our data support an ESCRT III-independent function for VPS4 at the centrosome and reveal that this evolutionary conserved AAA ATPase influences diverse centrosome functions and, as a result, global cellular architecture and development.

ALFY-Controlled DVL3 Autophagy Regulates Wnt Signaling, Determining Human Brain Size.

Primary microcephaly is a congenital neurodevelopmental disorder of reduced head circumference and brain volume, with fewer neurons in the cortex of the developing brain due to premature transition between symmetrical and asymmetrical cellular division of the neuronal stem cell layer during neurogenesis. We now show through linkage analysis and whole exome sequencing, that a dominant mutation in ALFY, encoding an autophagy scaffold protein, causes human primary microcephaly. We demonstrate the dominant effect of the mutation in drosophila: transgenic flies harboring the human mutant allele display small brain volume, recapitulating the disease phenotype. Moreover, eye-specific expression of human mutant ALFY causes rough eye phenotype. In molecular terms, we demonstrate that normally ALFY attenuates the canonical Wnt signaling pathway via autophagy-dependent removal specifically of aggregates of DVL3 and not of Dvl1 or Dvl2. Thus, autophagic attenuation of Wnt signaling through removal of Dvl3 aggregates by ALFY acts in determining human brain size.

UNC80 mutation causes a syndrome of hypotonia, severe intellectual disability, dyskinesia and dysmorphism, similar to that caused by mutations in its interacting cation channel NALCN.

BACKGROUND: A syndrome of profound hypotonia, intellectual disability, intrauterine growth retardation with subsequent failure to thrive, dyskinesia and epilepsy was diagnosed in Bedouin Israeli families. Mild dysmorphism was evident: plagiocephaly, broad forehead with prominent nose, smooth philtrum and congenital esotropia. We set out to decipher the molecular basis of this syndrome. METHODS: Genome-wide linkage analysis and fine mapping were done. Whole exome sequencing data were filtered for candidate variants within locus. Validation and segregation of the mutation was assayed via Sanger sequencing. UNC80 expression pattern was analysed through reverse transcription PCR. RESULTS: Homozygosity mapping followed by fine mapping identified a 7.5 Mb disease-associated locus (logarithm of odds score 3.5) on chromosome 2. Whole exome and Sanger sequencing identified a single homozygous nonsense mutation within this locus, segregating within the families as expected for recessive heredity and not found in a homozygous state in 150 Bedouin controls: c.151C>T, p.(R51*) in UNC80. CONCLUSIONS: The syndrome described is caused by a mutation in UNC80, truncating most of the 3258 amino acids highly conserved encoded protein, that has no known motifs. UNC80 bridges between UNC79 and the cation channel NALCN, enabling NALCN's role in basal Na(+) leak conductance in neurons, essential for neuronal function. The phenotype caused by the UNC80 mutation resembles that previously described for homozygous NALCN mutations.

Uncovering the Role of Hypermethylation by CTG Expansion in Myotonic Dystrophy Type 1 Using Mutant Human Embryonic Stem Cells.

CTG repeat expansion in DMPK, the cause of myotonic dystrophy type 1 (DM1), frequently results in hypermethylation and reduced SIX5 expression. The contribution of hypermethylation to disease pathogenesis and the precise mechanism by which SIX5 expression is reduced are unknown. Using 14 different DM1-affected human embryonic stem cell (hESC) lines, we characterized a differentially methylated region (DMR) near the CTGs. This DMR undergoes hypermethylation as a function of expansion size in a way that is specific to undifferentiated cells and is associated with reduced SIX5 expression. Using functional assays, we provide evidence for regulatory activity of the DMR, which is lost by hypermethylation and may contribute to DM1 pathogenesis by causing SIX5 haplo-insufficiency. This study highlights the power of hESCs in disease modeling and describes a DMR that functions both as an exon coding sequence and as a regulatory element whose activity is epigenetically hampered by a heritable mutation.

Dual Function of DNA Sequences: Protein-Coding Sequences Function as Transcriptional Enhancers.

Most of our genome comprises noncoding sequences that include diverse transcriptional regulatory elements, such as enhancers, while only ~1.5% of the genome codes for proteins. Nevertheless, DNA sequences that code for protein (exons) can also function as enhancers (eExons) that regulate transcription. Mutations in eExons can lead to multiple phenotypes due to their dual function. The prevalence of protein-coding sequences that possess transcriptional regulatory function (such as eExons) and the consequences of their mutations are not well described. Using advanced sequencing technologies, protein-coding sequences were analyzed for their potential regulatory function in mammalian cells and found to be overrepresented in the genome (>6%). Dissection of the enhancer activity of eExons at single nucleotide resolution in liver cells has demonstrated that: (1) most nucleotide changes with high impact effect are deleterious; (2) deleterious enhancer mutations are correlated with the location of transcription factor-binding sites; (3) synonymous and non-synonymous mutations have similar effects on enhancer activity; and (4) the transcription factor repertoire that controls the activity of enhancers differs across cell types, indicating differences in deleterious mutation profiles. Thus, eExon mutations can disrupt both protein structure and enhancer activity with differential effect across cell types, suggesting that a mutation in a gene could cause a phenotype that has nothing to do with its protein-coding function but is due to its additional hidden regulatory function.

Systematic dissection of coding exons at single nucleotide resolution supports an additional role in cell-specific transcriptional regulation.

In addition to their protein coding function, exons can also serve as transcriptional enhancers. Mutations in these exonic-enhancers (eExons) could alter both protein function and transcription. However, the functional consequence of eExon mutations is not well known. Here, using massively parallel reporter assays, we dissect the enhancer activity of three liver eExons (SORL1 exon 17, TRAF3IP2 exon 2, PPARG exon 6) at single nucleotide resolution in the mouse liver. We find that both synonymous and non-synonymous mutations have similar effects on enhancer activity and many of the deleterious mutation clusters overlap known liver-associated transcription factor binding sites. Carrying a similar massively parallel reporter assay in HeLa cells with these three eExons found differences in their mutation profiles compared to the liver, suggesting that enhancers could have distinct operating profiles in different tissues. Our results demonstrate that eExon mutations could lead to multiple phenotypes by disrupting both the protein sequence and enhancer activity and that enhancers can have distinct mutation profiles in different cell types.