A large-scale binding and functional map of human RNA-binding proteins.

Many proteins regulate the expression of genes by binding to specific regions encoded in the genome1. Here we introduce a new data set of RNA elements in the human genome that are recognized by RNA-binding proteins (RBPs), generated as part of the Encyclopedia of DNA Elements (ENCODE) project phase III. This class of regulatory elements functions only when transcribed into RNA, as they serve as the binding sites for RBPs that control post-transcriptional processes such as splicing, cleavage and polyadenylation, and the editing, localization, stability and translation of mRNAs. We describe the mapping and characterization of RNA elements recognized by a large collection of human RBPs in K562 and HepG2 cells. Integrative analyses using five assays identify RBP binding sites on RNA and chromatin in vivo, the in vitro binding preferences of RBPs, the function of RBP binding sites and the subcellular localization of RBPs, producing 1,223 replicated data sets for 356 RBPs. We describe the spectrum of RBP binding throughout the transcriptome and the connections between these interactions and various aspects of RNA biology, including RNA stability, splicing regulation and RNA localization. These data expand the catalogue of functional elements encoded in the human genome by the addition of a large set of elements that function at the RNA level by interacting with RBPs.

DPH1 syndrome: two novel variants and structural and functional analyses of seven missense variants identified in syndromic patients.

DPH1 variants have been associated with an ultra-rare and severe neurodevelopmental disorder, mainly characterized by variable developmental delay, short stature, dysmorphic features, and sparse hair. We have identified four new patients (from two different families) carrying novel variants in DPH1, enriching the clinical delineation of the DPH1 syndrome. Using a diphtheria toxin ADP-ribosylation assay, we have analyzed the activity of seven identified variants and demonstrated compromised function for five of them [p.(Leu234Pro); p.(Ala411Argfs*91); p.(Leu164Pro); p.(Leu125Pro); and p.(Tyr112Cys)]. We have built a homology model of the human DPH1-DPH2 heterodimer and have performed molecular dynamics simulations to study the effect of these variants on the catalytic sites as well as on the interactions between subunits of the heterodimer. The results show correlation between loss of activity, reduced size of the opening to the catalytic site, and changes in the size of the catalytic site with clinical severity. This is the first report of functional tests of DPH1 variants associated with the DPH1 syndrome. We demonstrate that the in vitro assay for DPH1 protein activity, together with structural modeling, are useful tools for assessing the effect of the variants on DPH1 function and may be used for predicting patient outcomes and prognoses.

Correction: DPH1 syndrome: two novel variants and structural and functional analyses of seven missense variants identified in syndromic patients.

Following the publication of the article, it was noted that the last column in Table 1, the total % should have read 5/8 (62.5) for the 'Epilepsy' row, and not 5.7 (71.4). This has now been amended in the HTML and PDF of the original article.

CeFra-seq reveals broad asymmetric mRNA and noncoding RNA distribution profiles in Drosophila and human cells.

Cells are highly asymmetrical, a feature that relies on the sorting of molecular constituents, including proteins, lipids, and nucleic acids, to distinct subcellular locales. The localization of RNA molecules is an important layer of gene regulation required to modulate localized cellular activities, although its global prevalence remains unclear. We combine biochemical cell fractionation with RNA-sequencing (CeFra-seq) analysis to assess the prevalence and conservation of RNA asymmetric distribution on a transcriptome-wide scale in Drosophila and human cells. This approach reveals that the majority (∼80%) of cellular RNA species are asymmetrically distributed, whether considering coding or noncoding transcript populations, in patterns that are broadly conserved evolutionarily. Notably, a large number of Drosophila and human long noncoding RNAs and circular RNAs display enriched levels within specific cytoplasmic compartments, suggesting that these RNAs fulfill extra-nuclear functions. Moreover, fraction-specific mRNA populations exhibit distinctive sequence characteristics. Comparative analysis of mRNA fractionation profiles with that of their encoded proteins reveals a general lack of correlation in subcellular distribution, marked by strong cases of asymmetry. However, coincident distribution profiles are observed for mRNA/protein pairs related to a variety of functional protein modules, suggesting complex regulatory inputs of RNA localization to cellular organization.

CeFra-seq: Systematic mapping of RNA subcellular distribution properties through cell fractionation coupled to deep-sequencing.

The subcellular trafficking of RNA molecules is a conserved feature of eukaryotic cells and plays key functions in diverse processes implicating polarised cellular activities. Large-scale imaging and subcellular transcriptomic studies suggest that regulated RNA localization is a highly prevalent process that appears to be disrupted in several neuromuscular disorders. These features underline the importance and usefulness of implementing procedures to assess global transcriptome subcellular distribution properties. Here, we present a method combining biochemical fractionation of cells and high-throughput RNA sequencing (CeFra-seq) that enables rapid and efficient systematic mapping of RNA cytotopic distributions in cells. The described procedure involves biochemical fractionation to derive extracts of nuclear, cytosolic, endomembrane, cytoplasmic insoluble and extracellular material from cell culture lines. The RNA content of each fraction can then be profiled by deep-sequencing, revealing global subcellular signatures. We provide a detailed protocol for the CeFra-seq procedure along with relevant validation steps and data analysis guidelines to graphically represent RNA spatial distribution features. As a complement to imaging approaches, CeFra-seq represents a powerful and scalable tool to investigate global alterations in RNA trafficking.

Developmentally regulated elimination of damaged nuclei involves a Chk2-dependent mechanism of mRNA nuclear retention.

The faithful execution of embryogenesis relies on the ability of organisms to respond to genotoxic stress and to eliminate defective cells that could otherwise compromise viability. In syncytial-stage Drosophila embryos, nuclei with excessive DNA damage undergo programmed elimination through an as-yet poorly understood process of nuclear fallout at the midblastula transition. We show that this involves a Chk2-dependent mechanism of mRNA nuclear retention that is induced by DNA damage and prevents the translation of specific zygotic mRNAs encoding key mitotic, cytoskeletal, and nuclear proteins required to maintain nuclear viability. For histone messages, we show that nuclear retention involves Chk2-mediated inactivation of the Drosophila stem loop binding protein (SLBP), the levels of which are specifically depleted in damaged nuclei following Chk2 phosphorylation, an event that contributes to nuclear fallout. These results reveal a layer of regulation within the DNA damage surveillance systems that safeguard genome integrity in eukaryotes.

The many functions of mRNA localization during normal development and disease: from pillar to post.

The regulated intracellular trafficking and localized translation of mRNA molecules represents an important and prevalent mechanism of gene regulation. This process plays a key role in modulating asymmetric protein distribution linked to a wide variety of biological processes in different organisms and cell types. In this review, we begin by discussing the diverse biological functions, advantages, and mechanisms of mRNA localization that have been characterized to date. We then review recent technological innovations in RNA imaging and functional genomics methods that will undoubtedly provide powerful new strategies for the elucidation of mRNA trafficking pathways. Finally, we discuss several examples linking human disease pathogenesis to defects in transcript localization, which further underlines the critical importance of this gene regulatory mechanism.

Transcriptome-wide regulation of pre-mRNA splicing and mRNA localization by muscleblind proteins.

The muscleblind-like (Mbnl) family of RNA-binding proteins plays important roles in muscle and eye development and in myotonic dystrophy (DM), in which expanded CUG or CCUG repeats functionally deplete Mbnl proteins. We identified transcriptome-wide functional and biophysical targets of Mbnl proteins in brain, heart, muscle, and myoblasts by using RNA-seq and CLIP-seq approaches. This analysis identified several hundred splicing events whose regulation depended on Mbnl function in a pattern indicating functional interchangeability between Mbnl1 and Mbnl2. A nucleotide resolution RNA map associated repression or activation of exon splicing with Mbnl binding near either 3' splice site or near the downstream 5' splice site, respectively. Transcriptomic analysis of subcellular compartments uncovered a global role for Mbnls in regulating localization of mRNAs in both mouse and Drosophila cells, and Mbnl-dependent translation and protein secretion were observed for a subset of mRNAs with Mbnl-dependent localization. These findings hold several new implications for DM pathogenesis.

Characterization of the 3p12.3-pcen region associated with tumor suppression in a novel ovarian cancer cell line model genetically modified by chromosome 3 fragment transfer.

The genetic analysis of nontumorigenic radiation hybrids generated by transfer of chromosome 3 fragments into the tumorigenic OV-90 ovarian cancer cell line identified the 3p12.3-pcen region as a candidate tumor suppressor gene (TSG) locus. In the present study, polymorphic microsatellite repeat analysis of the hybrids further defined the 3p12.3-pcen interval to a 16.1 Mb common region containing 12 known or hypothetical genes: 3ptel-ROBO2-ROBO1-GBE1-CADM2-VGLL3-CHMP2B-POU1F1-HTR1F-CGGBP1-ZNF654-C3orf38-EPHA3-3pcen. Seven of these genes, ROBO1, GBE1, VGLL3, CHMP2B, CGGBP1, ZNF654, and C3orf38, exhibited gene expression in the hybrids, placing them as top TSG candidates for further analysis. The expression of all but one (VGLL3) of these genes was also detected in the parental OV-90 cell line. Mutations were not identified in a comparative sequence analysis of the predicted protein coding regions of these candidates in OV-90 and donor normal chromosome 3 contig. However, the nondeleterious sequence variants identified in the transcribed regions distinguished parent of origin alleles for ROBO1, VGLL3, CHMP2B, and CGGBP1 and cDNA sequencing of the hybrids revealed biallelic expression of these genes. Interestingly, underexpression of VGLL3 and ZNF654 were observed in malignant ovarian tumor samples as compared with primary cultures of normal ovarian surface epithelial cells or benign ovarian tumors, and this occurred regardless of allelic content of 3p12.3-pcen. The results taken together suggest that dysregulation of VGLL3 and/or ZNF654 expression may have affected pathways important in ovarian tumorigenesis which was offset by the transfer of chromosome 3 fragments in OV-90, a cell line hemizygous for 3p.