Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements.

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

Identification of downstream effectors of retinoic acid specifying the zebrafish pancreas by integrative genomics.

Retinoic acid (RA) is a key signal for the specification of the pancreas. Still, the gene regulatory cascade triggered by RA in the endoderm remains poorly characterized. In this study, we investigated this regulatory network in zebrafish by combining RNA-seq, RAR ChIP-seq and ATAC-seq assays. By analysing the effect of RA and of the RA receptor (RAR) inverse-agonist BMS493 on the transcriptome and on the chromatin accessibility of endodermal cells, we identified a large set of genes and regulatory regions regulated by RA signalling. RAR ChIP-seq further defined the direct RAR target genes in zebrafish, including hox genes as well as several pancreatic regulators like mnx1, insm1b, hnf1ba and gata6. Comparison of zebrafish and murine RAR ChIP-seq data highlighted the conserved direct target genes and revealed that some RAR sites are under strong evolutionary constraints. Among them, a novel highly conserved RAR-induced enhancer was identified downstream of the HoxB locus and driving expression in the nervous system and in the gut in a RA-dependent manner. Finally, ATAC-seq data unveiled the role of the RAR-direct targets Hnf1ba and Gata6 in opening chromatin at many regulatory loci upon RA treatment.

Identification of novel potential interaction partners of UDP-galactose (SLC35A2), UDP-N-acetylglucosamine (SLC35A3) and an orphan (SLC35A4) nucleotide sugar transporters.

Nucleotide sugar transporters (NSTs) are ER and Golgi-resident members of the solute carrier 35 (SLC35) family which supply substrates for glycosylation by exchanging lumenal nucleotide monophosphates for cytosolic nucleotide sugars. Defective NSTs have been associated with congenital disorders of glycosylation (CDG), however, molecular basis of many types of CDG remains poorly characterized. To better understand the biology of NSTs, we identified potential interaction partners of UDP-galactose transporter (SLC35A2), UDP-N-acetylglucosamine transporter (SLC35A3) and an orphan nucleotide sugar transporter SLC35A4 of to date unassigned specificity. For this purpose, each of the SLC35A2-A4 proteins was used as a bait in four independent pull-down experiments and the identity of the immunoprecipitated material was discovered using MS techniques. From the candidate list obtained, we selected a few for which the interaction was confirmed in vitro using the NanoBiT system, a split luciferase-based luminescent technique. NSTs have been shown to interact with two ATPases (ATP2A2, ATP2C1), Golgi pH regulator B (GPR89B) and calcium channel (TMCO1), which may reflect the regulation of glycosylation by ion homeostasis, and with basigin (BSG). Our findings provide a starting point for the NST interaction network discovery in order to better understand how glycosylation is regulated and linked to other cellular processes. SIGNIFICANCE: Despite the facts that nucleotide sugar transporters are a key component of the protein glycosylation machinery, and deficiencies in their activity underlie serious metabolic diseases, biology, function and regulation of these essential proteins remain enigmatic. In this study we have advanced the field by identifying sets of new potential interaction partners for UDP-galactose transporter (SLC35A2), UDP-N-acetylglucosamine transporter (SLC35A3) and an orphan transporter SLC35A4 of yet undefined role. Several of these new interactions were additionally confirmed in vitro using the NanoBiT system, a split luciferase complementation assay. This work is also significant in that it addresses the overall challenge of discovering membrane protein interaction partners by a detailed comparison of 4 different co-immunoprecipitation strategies and by custom sample preparation and data processing workflows.

Germ cell differentiation requires Tdrd7-dependent chromatin and transcriptome reprogramming marked by germ plasm relocalization.

In many animal models, primordial germ cell (PGC) development depends on maternally deposited germ plasm, which prevents somatic cell fate. Here, we show that PGCs respond to regulatory information from the germ plasm in two distinct phases using two distinct mechanisms in zebrafish. We demonstrate that PGCs commence zygotic genome activation together with the somatic blastocysts with no demonstrable differences in transcriptional and chromatin opening. Unexpectedly, both PGC and somatic blastocysts activate germ-cell-specific genes, which are only stabilized in PGCs by cytoplasmic germ plasm determinants. Disaggregated perinuclear relocalization of germ plasm during PGC migration is regulated by the germ plasm determinant Tdrd7 and is coupled to dramatic divergence between PGC and somatic transcriptomes. This transcriptional divergence relies on PGC-specific cis-regulatory elements characterized by promoter-proximal distribution. We show that Tdrd7-dependent reconfiguration of chromatin accessibility is required for elaboration of PGC fate but not for PGC migration.

Dual-initiation promoters with intertwined canonical and TCT/TOP transcription start sites diversify transcript processing.

Variations in transcription start site (TSS) selection reflect diversity of preinitiation complexes and can impact on post-transcriptional RNA fates. Most metazoan polymerase II-transcribed genes carry canonical initiation with pyrimidine/purine (YR) dinucleotide, while translation machinery-associated genes carry polypyrimidine initiator (5'-TOP or TCT). By addressing the developmental regulation of TSS selection in zebrafish we uncovered a class of dual-initiation promoters in thousands of genes, including snoRNA host genes. 5'-TOP/TCT initiation is intertwined with canonical initiation and used divergently in hundreds of dual-initiation promoters during maternal to zygotic transition. Dual-initiation in snoRNA host genes selectively generates host and snoRNA with often different spatio-temporal expression. Dual-initiation promoters are pervasive in human and fruit fly, reflecting evolutionary conservation. We propose that dual-initiation on shared promoters represents a composite promoter architecture, which can function both coordinately and divergently to diversify RNAs.

Amphioxus functional genomics and the origins of vertebrate gene regulation.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that-in vertebrates-over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations.

Embryonic stem cell-specific microRNAs contribute to pluripotency by inhibiting regulators of multiple differentiation pathways.

The findings that microRNAs (miRNAs) are essential for early development in many species and that embryonic miRNAs can reprogram somatic cells into induced pluripotent stem cells suggest that these miRNAs act directly on transcriptional and chromatin regulators of pluripotency. To elucidate the transcription regulatory networks immediately downstream of embryonic miRNAs, we extended the motif activity response analysis approach that infers the regulatory impact of both transcription factors (TFs) and miRNAs from genome-wide expression states. Applying this approach to multiple experimental data sets generated from mouse embryonic stem cells (ESCs) that did or did not express miRNAs of the ESC-specific miR-290-295 cluster, we identified multiple TFs that are direct miRNA targets, some of which are known to be active during cell differentiation. Our results provide new insights into the transcription regulatory network downstream of ESC-specific miRNAs, indicating that these miRNAs act on cell cycle and chromatin regulators at several levels and downregulate TFs that are involved in the innate immune response.

The enhancer and promoter landscape of human regulatory and conventional T-cell subpopulations.

CD4(+)CD25(+)FOXP3(+) human regulatory T cells (Tregs) are essential for self-tolerance and immune homeostasis. Here, we describe the promoterome of CD4(+)CD25(high)CD45RA(+) naïve and CD4(+)CD25(high)CD45RA(-) memory Tregs and their CD25(-) conventional T-cell (Tconv) counterparts both before and after in vitro expansion by cap analysis of gene expression (CAGE) adapted to single-molecule sequencing (HeliScopeCAGE). We performed comprehensive comparative digital gene expression analyses and revealed novel transcription start sites, of which several were validated as alternative promoters of known genes. For all in vitro expanded subsets, we additionally generated global maps of poised and active enhancer elements marked by histone H3 lysine 4 monomethylation and histone H3 lysine 27 acetylation, describe their cell type-specific motif signatures, and evaluate the role of candidate transcription factors STAT5, FOXP3, RUNX1, and ETS1 in both Treg- and Tconv-specific enhancer architectures. Network analyses of gene expression data revealed additional candidate transcription factors contributing to cell type specificity and a transcription factor network in Tregs that is dominated by FOXP3 interaction partners and targets. In summary, we provide a comprehensive and easily accessible resource of gene expression and gene regulation in human Treg and Tconv subpopulations.

Transcription and enhancer profiling in human monocyte subsets.

Human blood monocytes comprise at least 3 subpopulations that differ in phenotype and function. Here, we present the first in-depth regulome analysis of human classical (CD14(++)CD16(-)), intermediate (CD14(+)CD16(+)), and nonclassical (CD14(dim)CD16(+)) monocytes. Cap analysis of gene expression adapted to Helicos single-molecule sequencing was used to map transcription start sites throughout the genome in all 3 subsets. In addition, global maps of H3K4me1 and H3K27ac deposition were generated for classical and nonclassical monocytes defining enhanceosomes of the 2 major subsets. We identified differential regulatory elements (including promoters and putative enhancers) that were associated with subset-specific motif signatures corresponding to different transcription factor activities and exemplarily validated novel downstream enhancer elements at the CD14 locus. In addition to known subset-specific features, pathway analysis revealed marked differences in metabolic gene signatures. Whereas classical monocytes expressed higher levels of genes involved in carbohydrate metabolism, priming them for anaerobic energy production, nonclassical monocytes expressed higher levels of oxidative pathway components and showed a higher mitochondrial routine activity. Our findings describe promoter/enhancer landscapes and provide novel insights into the specific biology of human monocyte subsets.

A promoter-level mammalian expression atlas.

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

ISMARA: automated modeling of genomic signals as a democracy of regulatory motifs.

Accurate reconstruction of the regulatory networks that control gene expression is one of the key current challenges in molecular biology. Although gene expression and chromatin state dynamics are ultimately encoded by constellations of binding sites recognized by regulators such as transcriptions factors (TFs) and microRNAs (miRNAs), our understanding of this regulatory code and its context-dependent read-out remains very limited. Given that there are thousands of potential regulators in mammals, it is not practical to use direct experimentation to identify which of these play a key role for a particular system of interest. We developed a methodology that models gene expression or chromatin modifications in terms of genome-wide predictions of regulatory sites and completely automated it into a web-based tool called ISMARA (Integrated System for Motif Activity Response Analysis). Given only gene expression or chromatin state data across a set of samples as input, ISMARA identifies the key TFs and miRNAs driving expression/chromatin changes and makes detailed predictions regarding their regulatory roles. These include predicted activities of the regulators across the samples, their genome-wide targets, enriched gene categories among the targets, and direct interactions between the regulators. Applying ISMARA to data sets from well-studied systems, we show that it consistently identifies known key regulators ab initio. We also present a number of novel predictions including regulatory interactions in innate immunity, a master regulator of mucociliary differentiation, TFs consistently disregulated in cancer, and TFs that mediate specific chromatin modifications.

An epigenetic profile of early T-cell development from multipotent progenitors to committed T-cell descendants.

Cellular differentiation of the T-cell branch of the immune system begins with the HSC, which undergoes a series of stages characterized by progressive restriction in multipotency and acquisition of specific lineage identity At the molecular level, the restriction of cell potential, commitment, and differentiation to a specific lineage is achieved through the coordinated control of gene expression and epigenetic mechanisms. Here, we analyzed and compared the gene expression profiles and the genome-wide histone modification marks H3K4me3 (H3 lysine 4 trimethylation) and H3K27me3 (H3 lysine 27 trimethylation) in (i) in vitro propagated HSCs, (ii) in vitro generated and propagated pro-T cells derived from these stem cells, and (iii) double-positive thymocytes derived from these pro-T cells after injection into Rag-deficient mice. The combined analyses of the different datasets in this unique experimental system highlighted the importance of both transcriptional and epigenetic repression in shaping the early phases of T-cell development.

Sox4 is a master regulator of epithelial-mesenchymal transition by controlling Ezh2 expression and epigenetic reprogramming.

Gene expression profiling has uncovered the transcription factor Sox4 with upregulated activity during TGF-ß-induced epithelial-mesenchymal transition (EMT) in normal and cancerous breast epithelial cells. Sox4 is indispensable for EMT and cell survival in vitro and for primary tumor growth and metastasis in vivo. Among several EMT-relevant genes, Sox4 directly regulates the expression of Ezh2, encoding the Polycomb group histone methyltransferase that trimethylates histone 3 lysine 27 (H3K27me3) for gene repression. Ablation of Ezh2 expression prevents EMT, whereas forced expression of Ezh2 restores EMT in Sox4-deficient cells. Ezh2-mediated H3K27me3 marks associate with key EMT genes, representing an epigenetic EMT signature that predicts patient survival. Our results identify Sox4 as a master regulator of EMT by governing the expression of the epigenetic modifier Ezh2.

Parity induces differentiation and reduces Wnt/Notch signaling ratio and proliferation potential of basal stem/progenitor cells isolated from mouse mammary epithelium.

INTRODUCTION: Early pregnancy has a strong protective effect against breast cancer in humans and rodents, but the underlying mechanism is unknown. Because breast cancers are thought to arise from specific cell subpopulations of mammary epithelia, we studied the effect of parity on the transcriptome and the differentiation/proliferation potential of specific luminal and basal mammary cells in mice. METHODS: Mammary epithelial cell subpopulations (luminal Sca1-, luminal Sca1+, basal stem/progenitor, and basal myoepithelial cells) were isolated by flow cytometry from parous and age-matched virgin mice and examined by using a combination of unbiased genomics, bioinformatics, in vitro colony formation, and in vivo limiting dilution transplantation assays. Specific findings were further investigated with immunohistochemistry in entire glands of parous and age-matched virgin mice. RESULTS: Transcriptome analysis revealed an upregulation of differentiation genes and a marked decrease in the Wnt/Notch signaling ratio in basal stem/progenitor cells of parous mice. Separate bioinformatics analyses showed reduced activity for the canonical Wnt transcription factor LEF1/TCF7 and increased activity for the Wnt repressor TCF3. This finding was specific for basal stem/progenitor cells and was associated with downregulation of potentially carcinogenic pathways and a reduction in the proliferation potential of this cell subpopulation in vitro and in vivo. As a possible mechanism for decreased Wnt signaling in basal stem/progenitor cells, we found a more than threefold reduction in the expression of the secreted Wnt ligand Wnt4 in total mammary cells from parous mice, which corresponded to a similar decrease in the proportion of Wnt4-secreting and estrogen/progesterone receptor-positive cells. Because recombinant Wnt4 rescued the proliferation defect of basal stem/progenitor cells in vitro, reduced Wnt4 secretion appears to be causally related to parity-induced alterations of basal stem/progenitor cell properties in mice. CONCLUSIONS: By revealing that parity induces differentiation and downregulates the Wnt/Notch signaling ratio and the in vitro and in vivo proliferation potential of basal stem/progenitor cells in mice, our study sheds light on the long-term consequences of an early pregnancy. Furthermore, it opens the door to future studies assessing whether inhibitors of the Wnt pathway may be used to mimic the parity-induced protective effect against breast cancer.

Modeling of epigenome dynamics identifies transcription factors that mediate Polycomb targeting.

Although changes in chromatin are integral to transcriptional reprogramming during cellular differentiation, it is currently unclear how chromatin modifications are targeted to specific loci. To systematically identify transcription factors (TFs) that can direct chromatin changes during cell fate decisions, we model the relationship between genome-wide dynamics of chromatin marks and the local occurrence of computationally predicted TF binding sites. By applying this computational approach to a time course of Polycomb-mediated H3K27me3 marks during neuronal differentiation of murine stem cells, we identify several motifs that likely regulate the dynamics of this chromatin mark. Among these, the sites bound by REST and by the SNAIL family of TFs are predicted to transiently recruit H3K27me3 in neuronal progenitors. We validate these predictions experimentally and show that absence of REST indeed causes loss of H3K27me3 at target promoters in trans, specifically at the neuronal progenitor state. Moreover, using targeted transgenic insertion, we show that promoter fragments containing REST or SNAIL binding sites are sufficient to recruit H3K27me3 in cis, while deletion of these sites results in loss of H3K27me3. These findings illustrate that the occurrence of TF binding sites can determine chromatin dynamics. Local determination of Polycomb activity by REST and SNAIL motifs exemplifies such TF based regulation of chromatin. Furthermore, our results show that key TFs can be identified ab initio through computational modeling of epigenome data sets using a modeling approach that we make readily accessible.

SwissRegulon, a database of genome-wide annotations of regulatory sites: recent updates.

Identification of genomic regulatory elements is essential for understanding the dynamics of cellular processes. This task has been substantially facilitated by the availability of genome sequences for many species and high-throughput data of transcripts and transcription factor (TF) binding. However, rigorous computational methods are necessary to derive accurate genome-wide annotations of regulatory sites from such data. SwissRegulon (http://swissregulon.unibas.ch) is a database containing genome-wide annotations of regulatory motifs, promoters and TF binding sites (TFBSs) in promoter regions across model organisms. Its binding site predictions were obtained with rigorous Bayesian probabilistic methods that operate on orthologous regions from related genomes, and use explicit evolutionary models to assess the evidence of purifying selection on each site. New in the current version of SwissRegulon is a curated collection of 190 mammalian regulatory motifs associated with ∼340 TFs, and TFBS annotations across a curated set of ∼35 000 promoters in both human and mouse. Predictions of TFBSs for Saccharomyces cerevisiae have also been significantly extended and now cover 158 of yeast's ∼180 TFs. All data are accessible through both an easily navigable genome browser with search functions, and as flat files that can be downloaded for further analysis.

The corepressor NCoR1 antagonizes PGC-1α and estrogen-related receptor α in the regulation of skeletal muscle function and oxidative metabolism.

Skeletal muscle exhibits a high plasticity and accordingly can quickly adapt to different physiological and pathological stimuli by changing its phenotype largely through diverse epigenetic mechanisms. The nuclear receptor corepressor 1 (NCoR1) has the ability to mediate gene repression; however, its role in regulating biological programs in skeletal muscle is still poorly understood. We therefore studied the mechanistic and functional aspects of NCoR1 function in this tissue. NCoR1 muscle-specific knockout mice exhibited a 7.2% higher peak oxygen consumption (VO(2peak)), a 11% reduction in maximal isometric force, and increased ex vivo fatigue resistance during maximal stimulation. Interestingly, global gene expression analysis revealed a high overlap between the effects of NCoR1 deletion and peroxisome proliferator-activated receptor gamma (PPARγ) coactivator 1α (PGC-1α) overexpression on oxidative metabolism in muscle. Importantly, PPARß/δ and estrogen-related receptor α (ERRα) were identified as common targets of NCoR1 and PGC-1α with opposing effects on the transcriptional activity of these nuclear receptors. In fact, the repressive effect of NCoR1 on oxidative phosphorylation gene expression specifically antagonizes PGC-1α-mediated coactivation of ERRα. We therefore delineated the molecular mechanism by which a transcriptional network controlled by corepressor and coactivator proteins determines the metabolic properties of skeletal muscle, thus representing a potential therapeutic target for metabolic diseases.

Adipose tissue microRNAs as regulators of CCL2 production in human obesity.

In obesity, white adipose tissue (WAT) inflammation is linked to insulin resistance. Increased adipocyte chemokine (C-C motif) ligand 2 (CCL2) secretion may initiate adipose inflammation by attracting the migration of inflammatory cells into the tissue. Using an unbiased approach, we identified adipose microRNAs (miRNAs) that are dysregulated in human obesity and assessed their possible role in controlling CCL2 production. In subcutaneous WAT obtained from 56 subjects, 11 miRNAs were present in all subjects and downregulated in obesity. Of these, 10 affected adipocyte CCL2 secretion in vitro and for 2 miRNAs (miR-126 and miR-193b), regulatory circuits were defined. While miR-126 bound directly to the 3'-untranslated region of CCL2 mRNA, miR-193b regulated CCL2 production indirectly through a network of transcription factors, many of which have been identified in other inflammatory conditions. In addition, overexpression of miR-193b and miR-126 in a human monocyte/macrophage cell line attenuated CCL2 production. The levels of the two miRNAs in subcutaneous WAT were significantly associated with CCL2 secretion (miR-193b) and expression of integrin, α-X, an inflammatory macrophage marker (miR-193b and miR-126). Taken together, our data suggest that miRNAs may be important regulators of adipose inflammation through their effects on CCL2 release from human adipocytes and macrophages.

Tyrosine phosphatase SHP2 promotes breast cancer progression and maintains tumor-initiating cells via activation of key transcription factors and a positive feedback signaling loop.

New cancer therapies are likely to arise from an in-depth understanding of the signaling networks influencing tumor initiation, progression and metastasis. We show a fundamental role for Src-homology 2 domain-containing phosphatase 2 (SHP2) in these processes in human epidermal growth factor receptor 2 (HER2)-positive and triple-negative breast cancers. Knockdown of SHP2 eradicated breast tumor-initiating cells in xenograft models, and SHP2 depletion also prevented invasion in three-dimensional cultures and in a transductal invasion assay in vivo. Notably, SHP2 knockdown in established breast tumors blocked their growth and reduced metastasis. Mechanistically, SHP2 activated stemness-associated transcription factors, including v-myc myelocytomatosis viral oncogene homolog (c-Myc) and zinc finger E-box binding homeobox 1 (ZEB1), which resulted in the repression of let-7 microRNA and the expression of a set of 'SHP2 signature' genes. We found these genes to be simultaneously activated in a large subset of human primary breast tumors that are associated with invasive behavior and poor prognosis. These results provide new insights into the signaling cascades influencing tumor-initiating cells as well as a rationale for targeting SHP2 in breast cancer.

The snoRNA MBII-52 (SNORD 115) is processed into smaller RNAs and regulates alternative splicing.

The loss of HBII-52 and related C/D box small nucleolar RNA (snoRNA) expression units have been implicated as a cause for the Prader-Willi syndrome (PWS). We recently found that the C/D box snoRNA HBII-52 changes the alternative splicing of the serotonin receptor 2C pre-mRNA, which is different from the traditional C/D box snoRNA function in non-mRNA methylation. Using bioinformatic predictions and experimental verification, we identified five pre-mRNAs (DPM2, TAF1, RALGPS1, PBRM1 and CRHR1) containing alternative exons that are regulated by MBII-52, the mouse homolog of HBII-52. Analysis of a single member of the MBII-52 cluster of snoRNAs by RNase protection and northern blot analysis shows that the MBII-52 expressing unit generates shorter RNAs that originate from the full-length MBII-52 snoRNA through additional processing steps. These novel RNAs associate with hnRNPs and not with proteins associated with canonical C/D box snoRNAs. Our data indicate that not a traditional C/D box snoRNA MBII-52, but a processed version lacking the snoRNA stem is the predominant MBII-52 RNA missing in PWS. This processed snoRNA functions in alternative splice-site selection. Its substitution could be a therapeutic principle for PWS.