Rare Variants Analyses Suggest Novel Cleft Genes in the African Population.

Non-syndromic orofacial clefts (NSOFCs) are common birth defects with a complex etiology. While over 60 common risk loci have been identified, they explain only a small proportion of the heritability for NSOFC. Rare variants have been implicated in the missing heritability. Thus, our study aimed to identify genes enriched with nonsynonymous rare coding variants associated with NSOFCs. Our sample included 814 non-syndromic cleft lip with or without palate (NSCL/P), 205 non-syndromic cleft palate only (NSCPO), and 2150 unrelated control children from Nigeria, Ghana, and Ethiopia. We conducted a gene-based analysis separately for each phenotype using three rare-variants collapsing models: (1) protein-altering (PA), (2) missense variants only (MO); and (3) loss of function variants only (LOFO). Subsequently, we utilized relevant transcriptomics data to evaluate associated gene expression and examined their mutation constraint using the gnomeAD database. In total, 13 genes showed suggestive associations (p = E-04). Among them, eight genes (ABCB1, ALKBH8, CENPF, CSAD, EXPH5, PDZD8, SLC16A9, and TTC28) were consistently expressed in relevant mouse and human craniofacial tissues during the formation of the face, and three genes (ABCB1, TTC28, and PDZD8) showed statistically significant mutation constraint. These findings underscore the role of rare variants in identifying candidate genes for NSOFCs.

CRISPR Activation Reverses Haploinsufficiency and Functional Deficits Caused by TTN Truncation Variants.

BACKGROUND: TTN truncation variants (TTNtvs) are the most common genetic lesion identified in individuals with dilated cardiomyopathy, a disease with high morbidity and mortality rates. TTNtvs reduce normal TTN (titin) protein levels, produce truncated proteins, and impair sarcomere content and function. Therapeutics targeting TTNtvs have been elusive because of the immense size of TTN, the rarity of specific TTNtvs, and incomplete knowledge of TTNtv pathogenicity. METHODS: We adapted CRISPR activation using dCas9-VPR to functionally interrogate TTNtv pathogenicity and develop a therapeutic in human cardiomyocytes and 3-dimensional cardiac microtissues engineered from induced pluripotent stem cell models harboring a dilated cardiomyopathy-associated TTNtv. We performed guide RNA screening with custom TTN reporter assays, agarose gel electrophoresis to quantify TTN protein levels and isoforms, and RNA sequencing to identify molecular consequences of TTN activation. Cardiomyocyte epigenetic assays were also used to nominate DNA regulatory elements to enable cardiomyocyte-specific TTN activation. RESULTS: CRISPR activation of TTN using single guide RNAs targeting either the TTN promoter or regulatory elements in spatial proximity to the TTN promoter through 3-dimensional chromatin interactions rescued TTN protein deficits disturbed by TTNtvs. Increasing TTN protein levels normalized sarcomere content and contractile function despite increasing truncated TTN protein. In addition to TTN transcripts, CRISPR activation also increased levels of myofibril assembly-related and sarcomere-related transcripts. CONCLUSIONS: TTN CRISPR activation rescued TTNtv-related functional deficits despite increasing truncated TTN levels, which provides evidence to support haploinsufficiency as a relevant genetic mechanism underlying heterozygous TTNtvs. CRISPR activation could be developed as a therapeutic to treat a large proportion of TTNtvs.

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.

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.

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.

Identifying key underlying regulatory networks and predicting targets of orphan C/D box SNORD116 snoRNAs in Prader-Willi syndrome.

Prader-Willi syndrome (PWS) is a rare neurodevelopmental disorder characterized principally by initial symptoms of neonatal hypotonia and failure-to-thrive in infancy, followed by hyperphagia and obesity. It is well established that PWS is caused by loss of paternal expression of the imprinted region on chromosome 15q11-q13. While most PWS cases exhibit megabase-scale deletions of the paternal chromosome 15q11-q13 allele, several PWS patients have been identified harboring a much smaller deletion encompassing primarily SNORD116. This finding suggests SNORD116 is a direct driver of PWS phenotypes. The SNORD116 gene cluster is composed of 30 copies of individual SNORD116 C/D box small nucleolar RNAs (snoRNAs). Many C/D box snoRNAs have been shown to guide chemical modifications of other RNA molecules, often ribosomal RNA (rRNA). However, SNORD116 snoRNAs are termed 'orphans' because no verified targets have been identified and their sequences show no significant complementarity to rRNA. It is crucial to identify the targets and functions of SNORD116 snoRNAs because all reported PWS cases lack their expression. To address this, we engineered two different deletions modelling PWS in two distinct human embryonic stem cell (hESC) lines to control for effects of genetic background. Utilizing an inducible expression system enabled quick, reproducible differentiation of these lines into neurons. Systematic comparisons of neuronal gene expression across deletion types and genetic backgrounds revealed a novel list of 42 consistently dysregulated genes. Employing the recently described computational tool snoGloBe, we discovered these dysregulated genes are significantly enriched for predicted SNORD116 targeting versus multiple control analyses. Importantly, our results showed it is critical to use multiple isogenic cell line pairs, as this eliminated many spuriously differentially expressed genes. Our results indicate a novel gene regulatory network controlled by SNORD116 is likely perturbed in PWS patients.

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.

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.

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.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes.

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

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.

Zebrafish models of alx-linked frontonasal dysplasia reveal a role for Alx1 and Alx3 in the anterior segment and vasculature of the developing eye.

The cellular and genetic mechanisms that coordinate formation of facial sensory structures with surrounding skeletal and soft tissue elements remain poorly understood. Alx1, a homeobox transcription factor, is a key regulator of midfacial morphogenesis. ALX1 mutations in humans are linked to severe congenital anomalies of the facial skeleton (frontonasal dysplasia, FND) with malformation or absence of eyes and orbital contents (micro- and anophthalmia). Zebrafish with loss-of-function alx1 mutations develop with craniofacial and ocular defects of variable penetrance, likely due to compensatory upregulation in expression of a paralogous gene, alx3. Here we show that zebrafish alx1;alx3 mutants develop with highly penetrant cranial and ocular defects that resemble human ALX1-linked FND. alx1 and alx3 are expressed in anterior cranial neural crest (aCNC), which gives rise to the anterior neurocranium (ANC), anterior segment structures of the eye and vascular pericytes. Consistent with a functional requirement for alx genes in aCNC, alx1; alx3 mutants develop with nearly absent ANC and grossly aberrant hyaloid vasculature and ocular anterior segment, but normal retina. In vivo lineage labeling identified a requirement for alx1 and alx3 during aCNC migration, and transcriptomic analysis suggested oxidative stress response as a key target mechanism of this function. Oxidative stress is a hallmark of fetal alcohol toxicity, and we found increased penetrance of facial and ocular malformations in alx1 mutants exposed to ethanol, consistent with a protective role for alx1 against ethanol toxicity. Collectively, these data demonstrate a conserved role for zebrafish alx genes in controlling ocular and facial development, and a novel role in protecting these key midfacial structures from ethanol toxicity during embryogenesis. These data also reveal novel roles for alx genes in ocular anterior segment formation and vascular development and suggest that retinal deficits in alx mutants may be secondary to aberrant ocular vascularization and anterior segment defects. This study establishes robust zebrafish models for interrogating conserved genetic mechanisms that coordinate facial and ocular development, and for exploring gene--environment interactions relevant to fetal alcohol syndrome.

Integration of multimodal data in the developing tooth reveals candidate regulatory loci driving human odontogenic phenotypes.

Human odontogenic aberrations such as abnormal tooth number and delayed tooth eruption can occur as a symptom of rare syndromes or, more commonly, as nonsyndromic phenotypes. These phenotypes can require extensive and expensive dental treatment, posing a significant burden. While many dental phenotypes are heritable, most nonsyndromic cases have not been linked to causal genes. We demonstrate the novel finding that common sequence variants associated with human odontogenic phenotypes are enriched in developmental craniofacial enhancers conserved between human and mouse. However, the bulk nature of these samples obscures if this finding is due to the tooth itself or the surrounding tissues. We therefore sought to identify enhancers specifically active in the tooth anlagen and quantify their contribution to the observed genetic enrichments. We systematically identified 22,001 conserved enhancers active in E13.5 mouse incisors using ChIP-seq and machine learning pipelines and demonstrated biologically relevant enrichments in putative target genes, transcription factor binding motifs, and in vivo activity. Multi-tissue comparisons of human and mouse enhancers revealed that these putative tooth enhancers had the strongest enrichment of odontogenic phenotype-associated variants, suggesting a role for dysregulation of tooth developmental enhancers in human dental phenotypes. The large number of these regions genome-wide necessitated prioritization of enhancer loci for future investigations. As enhancers modulate gene expression, we prioritized regions based on enhancers' putative target genes. We predicted these target genes and prioritized loci by integrating chromatin state, bulk gene expression and coexpression, GWAS variants, and cell type resolved gene expression to generate a prioritized list of putative odontogenic phenotype-driving loci active in the developing tooth. These genomic regions are of particular interest for downstream experiments determining the role of specific dental enhancer:gene pairs in odontogenesis.

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.

Massively parallel discovery of human-specific substitutions that alter enhancer activity.

Genetic changes that altered the function of gene regulatory elements have been implicated in the evolution of human traits such as the expansion of the cerebral cortex. However, identifying the particular changes that modified regulatory activity during human evolution remain challenging. Here we used massively parallel enhancer assays in neural stem cells to quantify the functional impact of >32,000 human-specific substitutions in >4,300 human accelerated regions (HARs) and human gain enhancers (HGEs), which include enhancers with novel activities in humans. We found that >30% of active HARs and HGEs exhibited differential activity between human and chimpanzee. We isolated the effects of human-specific substitutions from background genetic variation to identify the effects of genetic changes most relevant to human evolution. We found that substitutions interacted in both additive and nonadditive ways to modify enhancer function. Substitutions within HARs, which are highly constrained compared to HGEs, showed smaller effects on enhancer activity, suggesting that the impact of human-specific substitutions is buffered in enhancers with constrained ancestral functions. Our findings yield insight into how human-specific genetic changes altered enhancer function and provide a rich set of candidates for studies of regulatory evolution in humans.

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate.

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

An Irf6-Esrp1/2 regulatory axis controls midface morphogenesis in vertebrates.

Irf6 and Esrp1 are important for palate development across vertebrates. In zebrafish, we found that irf6 regulates the expression of esrp1 We detailed overlapping Irf6 and Esrp1/2 expression in mouse orofacial epithelium. In zebrafish, irf6 and esrp1/2 share expression in periderm, frontonasal ectoderm and oral epithelium. Genetic disruption of irf6 and esrp1/2 in zebrafish resulted in cleft of the anterior neurocranium. The esrp1/2 mutant also developed cleft of the mouth opening. Lineage tracing of cranial neural crest cells revealed that the cleft resulted not from migration defect, but from impaired chondrogenesis. Analysis of aberrant cells within the cleft revealed expression of sox10, col1a1 and irf6, and these cells were adjacent to krt4+ and krt5+ cells. Breeding of mouse Irf6; Esrp1; Esrp2 compound mutants suggested genetic interaction, as the triple homozygote and the Irf6; Esrp1 double homozygote were not observed. Further, Irf6 heterozygosity reduced Esrp1/2 cleft severity. These studies highlight the complementary analysis of Irf6 and Esrp1/2 in mouse and zebrafish, and identify a unique aberrant cell population in zebrafish expressing sox10, col1a1 and irf6 Future work characterizing this cell population will yield additional insight into cleft pathogenesis.

ALX1-related frontonasal dysplasia results from defective neural crest cell development and migration.

A pedigree of subjects presented with frontonasal dysplasia (FND). Genome sequencing and analysis identified a p.L165F missense variant in the homeodomain of the transcription factor ALX1 which was imputed to be pathogenic. Induced pluripotent stem cells (iPSC) were derived from the subjects and differentiated to neural crest cells (NCC). NCC derived from ALX1L165F/L165F iPSC were more sensitive to apoptosis, showed an elevated expression of several neural crest progenitor state markers, and exhibited impaired migration compared to wild-type controls. NCC migration was evaluated in vivo using lineage tracing in a zebrafish model, which revealed defective migration of the anterior NCC stream that contributes to the median portion of the anterior neurocranium, phenocopying the clinical presentation. Analysis of human NCC culture media revealed a change in the level of bone morphogenic proteins (BMP), with a low level of BMP2 and a high level of BMP9. Soluble BMP2 and BMP9 antagonist treatments were able to rescue the defective migration phenotype. Taken together, these results demonstrate a mechanistic requirement of ALX1 in NCC development and migration.

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.