Massively parallel disruption of enhancers active in human neural stem cells.

Changes in gene regulation have been linked to the expansion of the human cerebral cortex and to neurodevelopmental disorders, potentially by altering neural progenitor proliferation. However, the effects of genetic variation within regulatory elements on neural progenitors remain obscure. We use sgRNA-Cas9 screens in human neural stem cells (hNSCs) to disrupt 10,674 genes and 26,385 conserved regions in 2,227 enhancers active in the developing human cortex and determine effects on proliferation. Genes with proliferation phenotypes are associated with neurodevelopmental disorders and show biased expression in specific fetal human brain neural progenitor populations. Although enhancer disruptions overall have weaker effects than gene disruptions, we identify enhancer disruptions that severely alter hNSC self-renewal. Disruptions in human accelerated regions, implicated in human brain evolution, also alter proliferation. Integrating proliferation phenotypes with chromatin interactions reveals regulatory relationships between enhancers and their target genes contributing to neurogenesis and potentially to human cortical evolution.

Spatial Multiomics Reveal the Role of Wnt Modulator, Dkk2, in Palatogenesis.

Multiple genetic and environmental etiologies contribute to the pathogenesis of cleft palate, which constitutes the most common among the inherited disorders of the craniofacial complex. Insights into the molecular mechanisms regulating osteogenic differentiation and patterning in the palate during embryogenesis are limited and needed for the development of innovative diagnostics and cures. This study utilized the Pax9-/- mouse model with a consistent phenotype of cleft secondary palate to investigate the role of Pax9 in the process of palatal osteogenesis. While prior research had identified upregulation of Wnt pathway modulators Dkk1 and Dkk2 in Pax9-/- palate mesenchyme, limitations of spatial resolution and technology restricted a more robust analysis. Here, data from single-nucleus transcriptomics and chromatin accessibility assays validated by in situ highly multiplex targeted single-cell spatial profiling technology suggest a distinct relationship between Pax9+ and osteogenic populations. Loss of Pax9 results in spatially restricted osteogenic domains bounded by Dkk2, which normally interfaces with Pax9 in the mesenchyme. These results suggest that Pax9-dependent Wnt signaling modulators influence osteogenic programming during palate formation, potentially contributing to the observed cleft palate phenotype.

Generation of isogenic models of Angelman syndrome and Prader-Willi syndrome in CRISPR/Cas9-engineered human embryonic stem cells.

Angelman Syndrome (AS) and Prader-Willi Syndrome (PWS), two distinct neurodevelopmental disorders, result from loss of expression from imprinted genes in the chromosome 15q11-13 locus most commonly caused by a megabase-scale deletion on either the maternal or paternal allele, respectively. Each occurs at an approximate incidence of 1/15,000 to 1/30,000 live births and has a range of debilitating phenotypes. Patient-derived induced pluripotent stem cells (iPSCs) have been valuable tools to understand human-relevant gene regulation at this locus and have contributed to the development of therapeutic approaches for AS. Nonetheless, gaps remain in our understanding of how these deletions contribute to dysregulation and phenotypes of AS and PWS. Variability across cell lines due to donor differences, reprogramming methods, and genetic background make it challenging to fill these gaps in knowledge without substantially increasing the number of cell lines used in the analyses. Isogenic cell lines that differ only by the genetic mutation causing the disease can ease this burden without requiring such a large number of cell lines. Here, we describe the development of isogenic human embryonic stem cell (hESC) lines modeling the most common genetic subtypes of AS and PWS. These lines allow for a facile interrogation of allele-specific gene regulation at the chromosome 15q11-q13 locus. Additionally, these lines are an important resource to identify and test targeted therapeutic approaches for patients with AS and PWS.

Multimodal spatiotemporal transcriptomic resolution of embryonic palate osteogenesis.

The terminal differentiation of osteoblasts and subsequent formation of bone marks an important phase in palate development that leads to the separation of the oral and nasal cavities. While the morphogenetic events preceding palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms driving the formation of this bony union of the fusing palate. Through bulk, single-nucleus, and spatially resolved RNA-sequencing analyses of the developing secondary palate, we identify a shift in transcriptional programming between embryonic days 14.5 and 15.5 pinpointing the onset of osteogenesis. We define spatially restricted expression patterns of key osteogenic marker genes that are differentially expressed between these developmental timepoints. Finally, we identify genes in the palate highly expressed by palate nasal epithelial cells, also enriched within palatal osteogenic mesenchymal cells. This investigation provides a relevant framework to advance palate-specific diagnostic and therapeutic biomarker discovery.

Rare variant modifier analysis identifies variants in SEC24D associated with orofacial cleft subtypes.

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p = 6.86 [Formula: see text] 10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p = 0.0009), and created binding sites for 23 TFs, including Pax6 (p = 6.12 [Formula: see text] 10-5) and Pax9 (p = 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study indicates that rare genetic variation may contribute to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Rare genetic variants in SEC24D modify orofacial cleft phenotypes.

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p=6.86×10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p=0.0009), and created binding sites for 23 TFs, including Pax6 (p=6.12×10-5) and Pax9 (p= 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study demonstrates that rare genetic variation contributes to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Sarcomere function activates a p53-dependent DNA damage response that promotes polyploidization and limits in vivo cell engraftment.

Human cardiac regeneration is limited by low cardiomyocyte replicative rates and progressive polyploidization by unclear mechanisms. To study this process, we engineer a human cardiomyocyte model to track replication and polyploidization using fluorescently tagged cyclin B1 and cardiac troponin T. Using time-lapse imaging, in vitro cardiomyocyte replication patterns recapitulate the progressive mononuclear polyploidization and replicative arrest observed in vivo. Single-cell transcriptomics and chromatin state analyses reveal that polyploidization is preceded by sarcomere assembly, enhanced oxidative metabolism, a DNA damage response, and p53 activation. CRISPR knockout screening reveals p53 as a driver of cell-cycle arrest and polyploidization. Inhibiting sarcomere function, or scavenging ROS, inhibits cell-cycle arrest and polyploidization. Finally, we show that cardiomyocyte engraftment in infarcted rat hearts is enhanced 4-fold by the increased proliferation of troponin-knockout cardiomyocytes. Thus, the sarcomere inhibits cell division through a DNA damage response that can be targeted to improve cardiomyocyte replacement strategies.

Specific ZNF274 binding interference at SNORD116 activates the maternal transcripts in Prader-Willi syndrome neurons.

Prader-Willi syndrome (PWS) is characterized by neonatal hypotonia, developmental delay and hyperphagia/obesity. This disorder is caused by the absence of paternally expressed gene products from chromosome 15q11-q13. We previously demonstrated that knocking out ZNF274, a Kruppel-associated box-A-domain zinc finger protein capable of recruiting epigenetic machinery to deposit the H3K9me3 repressive histone modification, can activate expression from the normally silent maternal allele of SNORD116 in neurons derived from PWS induced pluripotent stem cells (iPSCs). However, ZNF274 has many other targets in the genome in addition to SNORD116. Depleting ZNF274 will surely affect the expression of other important genes and disrupt other pathways. Here, we used CRISPR/Cas9 to delete ZNF274 binding sites at the SNORD116 locus to determine whether activation of the maternal copy of SNORD116 could be achieved without altering ZNF274 protein levels. We obtained similar activation of gene expression from the normally silenced maternal allele in neurons derived from PWS iPSCs, compared with ZNF274 knockout, demonstrating that ZNF274 is directly involved in the repression of SNORD116. These results suggest that interfering with ZNF274 binding at the maternal SNORD116 locus is a potential therapeutic strategy for PWS.

Evidence against tetrapod-wide digit identities and for a limited frame shift in bird wings.

In crown group tetrapods, individual digits are homologized in relation to a pentadactyl ground plan. However, testing hypotheses of digit homology is challenging because it is unclear whether digits represent distinct and conserved gene regulatory states. Here we show dramatic evolutionary dynamism in the gene expression profiles of digits, challenging the notion that five digits have conserved developmental identities across amniotes. Transcriptomics shows diversity in the patterns of gene expression differentiation of digits, although the anterior-most digit of the pentadactyl limb has a unique, conserved expression profile. Further, we identify a core set of transcription factors that are differentially expressed among the digits of amniote limbs; their spatial expression domains, however, vary between species. In light of these results, we reevaluate the frame shift hypothesis of avian wing evolution and conclude only the identity of the anterior-most digit has shifted position, suggesting a 1,3,4 digit identity in the bird wing.

Multi-focal control of mitochondrial gene expression by oncogenic MYC provides potential therapeutic targets in cancer.

Despite ubiquitous activation in human cancer, essential downstream effector pathways of the MYC transcription factor have been difficult to define and target. Using a structure/function-based approach, we identified the mitochondrial RNA polymerase (POLRMT) locus as a critical downstream target of MYC. The multifunctional POLRMT enzyme controls mitochondrial gene expression, a process required both for mitochondrial function and mitochondrial biogenesis. We further demonstrate that inhibition of this newly defined MYC effector pathway causes robust and selective tumor cell apoptosis, via an acute, checkpoint-like mechanism linked to aberrant electron transport chain complex assembly and mitochondrial reactive oxygen species (ROS) production. Fortuitously, MYC-dependent tumor cell death can be induced by inhibiting the mitochondrial gene expression pathway using a variety of strategies, including treatment with FDA-approved antibiotics. In vivo studies using a mouse model of Burkitt's Lymphoma provide pre-clinical evidence that these antibiotics can successfully block progression of MYC-dependent tumors.

Chromatin immunoprecipitation with fixed animal tissues and preparation for high-throughput sequencing.

Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) is a powerful method used to identify genome-wide binding patterns of transcription factors and distribution of various histone modifications associated with different chromatin states. In most published studies, ChIP-Seq has been performed on cultured cells grown under controlled conditions, allowing generation of large amounts of material in a homogeneous biological state. Although such studies have provided great insight into the dynamic landscapes of animal genomes, they do not allow the examination of transcription factor binding and chromatin states in adult tissues, developing embryonic structures, or tumors. Such knowledge is critical to understanding the information required to create and maintain a complex biological tissue and to identify noncoding regions of the genome directly involved in tissues affected by complex diseases such as autism. Studying these tissue types with ChIP-Seq can be challenging due to the limited availability of tissues and the lack of complex biological states able to be achieved in culture. These inherent differences require alterations of standard cross-linking and chromatin extraction typically used in cell culture. Here we describe a general approach for using small amounts of animal tissue to perform ChIP-Seq directed at histone modifications and transcription factors. Tissue is homogenized before treatment with formaldehyde to ensure proper cross-linking, and a two-step nuclear isolation is performed to increase extraction of soluble chromatin. Small amounts of soluble chromatin are then used for immunoprecipitation (IP) and prepared for multiplexed high-throughput sequencing.