Identifying key underlying regulatory networks and predicting targets of orphan C/D box <i>SNORD116</i> snoRNAs in Prader–Willi syndrome

Prader-Willi syndrome (PWS) is a rare neurodevelopmental disorder characterized by neonatal hypotonia, followed by hyperphagia and obesity. Most PWS cases exhibit megabase-scale deletions of paternally imprinted 15q11-q13 locus. However, several PWS patients have been identified harboring much smaller deletions encompassing the SNORD116 gene cluster, suggesting these genes are direct drivers of PWS phenotypes. This cluster contains 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 RNA molecules, often ribosomal RNA (rRNA). Conversely, SNORD116 snoRNAs show no significant complementarity to rRNA and their targets are unknown. Since many reported PWS cases lack their expression, it is crucial to identify the targets and functions of SNORD116. To address this we modeled PWS in two distinct human embryonic stem cell (hESC) lines with two different sized deletions, differentiated each into neurons, and compared differential gene expression. This analysis identified a novel set of 42 consistently dysregulated genes. These genes were significantly enriched for predicted SNORD116 targeting and we demonstrated impacts on FGF13 protein levels. Our results demonstrate the need for isogenic background comparisons and indicate a novel gene regulatory network controlled by SNORD116 is likely perturbed in PWS patients.