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.

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.

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.

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.

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.

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.

Chromatin connectivity maps reveal dynamic promoter-enhancer long-range associations.

In multicellular organisms, transcription regulation is one of the central mechanisms modelling lineage differentiation and cell-fate determination. Transcription requires dynamic chromatin configurations between promoters and their corresponding distal regulatory elements. It is believed that their communication occurs within large discrete foci of aggregated RNA polymerases termed transcription factories in three-dimensional nuclear space. However, the dynamic nature of chromatin connectivity has not been characterized at the genome-wide level. Here, through a chromatin interaction analysis with paired-end tagging approach using an antibody that primarily recognizes the pre-initiation complexes of RNA polymerase II, we explore the transcriptional interactomes of three mouse cells of progressive lineage commitment, including pluripotent embryonic stem cells, neural stem cells and neurosphere stem/progenitor cells. Our global chromatin connectivity maps reveal approximately 40,000 long-range interactions, suggest precise enhancer-promoter associations and delineate cell-type-specific chromatin structures. Analysis of the complex regulatory repertoire shows that there are extensive colocalizations among promoters and distal-acting enhancers. Most of the enhancers associate with promoters located beyond their nearest active genes, indicating that the linear juxtaposition is not the only guiding principle driving enhancer target selection. Although promoter-enhancer interactions exhibit high cell-type specificity, promoters involved in interactions are found to be generally common and mostly active among different cells. Chromatin connectivity networks reveal that the pivotal genes of reprogramming functions are transcribed within physical proximity to each other in embryonic stem cells, linking chromatin architecture to coordinated gene expression. Our study sets the stage for the full-scale dissection of spatial and temporal genome structures and their roles in orchestrating 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.

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.

Coding exons function as tissue-specific enhancers of nearby genes.

Enhancers are essential gene regulatory elements whose alteration can lead to morphological differences between species, developmental abnormalities, and human disease. Current strategies to identify enhancers focus primarily on noncoding sequences and tend to exclude protein coding sequences. Here, we analyzed 25 available ChIP-seq data sets that identify enhancers in an unbiased manner (H3K4me1, H3K27ac, and EP300) for peaks that overlap exons. We find that, on average, 7% of all ChIP-seq peaks overlap coding exons (after excluding for peaks that overlap with first exons). By using mouse and zebrafish enhancer assays, we demonstrate that several of these exonic enhancer (eExons) candidates can function as enhancers of their neighboring genes and that the exonic sequence is necessary for enhancer activity. Using ChIP, 3C, and DNA FISH, we further show that one of these exonic limb enhancers, Dync1i1 exon 15, has active enhancer marks and physically interacts with Dlx5/6 promoter regions 900 kb away. In addition, its removal by chromosomal abnormalities in humans could cause split hand and foot malformation 1 (SHFM1), a disorder associated with DLX5/6. These results demonstrate that DNA sequences can have a dual function, operating as coding exons in one tissue and enhancers of nearby gene(s) in another tissue, suggesting that phenotypes resulting from coding mutations could be caused not only by protein alteration but also by disrupting the regulation of another gene.

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.

Seborrhea-like dermatitis with psoriasiform elements caused by a mutation in ZNF750, encoding a putative C2H2 zinc finger protein.

We describe an Israeli Jewish Moroccan family presenting with autosomal dominant seborrhea-like dermatosis with psoriasiform elements, including enhanced keratinocyte proliferation, parakeratosis, follicular plugging, Pityrosporum ovale overgrowth and dermal CD4 lymphocyte infiltrate. We mapped the disease gene to a 0.5-cM region overlapping the PSORS2 locus (17q25) and identified a frameshift mutation in ZNF750, which encodes a putative C2H2 zinc finger protein. ZNF750 is normally expressed in keratinocytes but not in fibroblasts and is barely detectable in CD4 lymphocytes.

Functional characterization of tissue-specific enhancers in the DLX5/6 locus.

Disruption of distaless homeobox 5 and 6 (Dlx5/6) in mice results in brain, craniofacial, genital, ear and limb defects. In humans, chromosomal aberrations in the DLX5/6 region, some of which do not encompass DLX5/6, are associated with split hand/foot malformation 1 (SHFM1) as well as intellectual disability, craniofacial anomalies and hearing loss, suggesting that the disruption of DLX5/6 regulatory elements could lead to these abnormalities. Here, we characterized enhancers in the DLX5/6 locus whose tissue-specific expression and genomic location along with previously characterized enhancers correlate with phenotypes observed in individuals with chromosomal abnormalities. By analyzing chromosomal aberrations at 7q21, we refined the minimal SHFM1 critical region and used comparative genomics to select 26 evolutionary conserved non-coding sequences in this critical region for zebrafish enhancer assays. Eight of these sequences were shown to function as brain, olfactory bulb, branchial arch, otic vesicle and fin enhancers, recapitulating dlx5a/6a expression. Using a mouse enhancer assay, several of these zebrafish enhancers showed comparable expression patterns in the branchial arch, otic vesicle, forebrain and/or limb at embryonic day 11.5. Examination of the coordinates of various chromosomal rearrangements in conjunction with the genomic location of these tissue-specific enhancers showed a correlation with the observed clinical abnormalities. Our findings suggest that chromosomal abnormalities that disrupt the function of these tissue-specific enhancers could be the cause of SHFM1 and its associated phenotypes. In addition, they highlight specific enhancers in which mutations could lead to non-syndromic hearing loss, craniofacial defects or limb malformations.

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.

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.

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.

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 ).