Initiation of Parental Genome Reprogramming in Fertilized Oocyte by Splicing Kinase SRPK1-Catalyzed Protamine Phosphorylation.

The paternal genome undergoes a massive exchange of histone with protamine for compaction into sperm during spermiogenesis. Upon fertilization, this process is potently reversed, which is essential for parental genome reprogramming and subsequent activation; however, it remains poorly understood how this fundamental process is initiated and regulated. Here, we report that the previously characterized splicing kinase SRPK1 initiates this life-beginning event by catalyzing site-specific phosphorylation of protamine, thereby triggering protamine-to-histone exchange in the fertilized oocyte. Interestingly, protamine undergoes a DNA-dependent phase transition to gel-like condensates and SRPK1-mediated phosphorylation likely helps open up such structures to enhance protamine dismissal by nucleoplasmin (NPM2) and enable the recruitment of HIRA for H3.3 deposition. Remarkably, genome-wide assay for transposase-accessible chromatin sequencing (ATAC-seq) analysis reveals that selective chromatin accessibility in both sperm and MII oocytes is largely erased in early pronuclei in a protamine phosphorylation-dependent manner, suggesting that SRPK1-catalyzed phosphorylation initiates a highly synchronized reorganization program in both parental genomes.

Pervasive Chromatin-RNA Binding Protein Interactions Enable RNA-Based Regulation of Transcription.

Increasing evidence suggests that transcriptional control and chromatin activities at large involve regulatory RNAs, which likely enlist specific RNA-binding proteins (RBPs). Although multiple RBPs have been implicated in transcription control, it has remained unclear how extensively RBPs directly act on chromatin. We embarked on a large-scale RBP ChIP-seq analysis, revealing widespread RBP presence in active chromatin regions in the human genome. Like transcription factors (TFs), RBPs also show strong preference for hotspots in the genome, particularly gene promoters, where their association is frequently linked to transcriptional output. Unsupervised clustering reveals extensive co-association between TFs and RBPs, as exemplified by YY1, a known RNA-dependent TF, and RBM25, an RBP involved in splicing regulation. Remarkably, RBM25 depletion attenuates all YY1-dependent activities, including chromatin binding, DNA looping, and transcription. We propose that various RBPs may enhance network interaction through harnessing regulatory RNAs to control transcription.

SR proteins collaborate with 7SK and promoter-associated nascent RNA to release paused polymerase.

RNAP II is frequently paused near gene promoters in mammals, and its transition to productive elongation requires active recruitment of P-TEFb, a cyclin-dependent kinase for RNAP II and other key transcription elongation factors. A fraction of P-TEFb is sequestered in an inhibitory complex containing the 7SK noncoding RNA, but it has been unclear how P-TEFb is switched from the 7SK complex to RNAP II during transcription activation. We report that SRSF2 (also known as SC35, an SR-splicing factor) is part of the 7SK complex assembled at gene promoters and plays a direct role in transcription pause release. We demonstrate RNA-dependent, coordinated release of SRSF2 and P-TEFb from the 7SK complex and transcription activation via SRSF2 binding to promoter-associated nascent RNA. These findings reveal an unanticipated SR protein function, a role for promoter-proximal nascent RNA in gene activation, and an analogous mechanism to HIV Tat/TAR for activating cellular genes.

Direct conversion of fibroblasts to neurons by reprogramming PTB-regulated microRNA circuits.

The induction of pluripotency or trans-differentiation of one cell type to another can be accomplished with cell-lineage-specific transcription factors. Here, we report that repression of a single RNA binding polypyrimidine-tract-binding (PTB) protein, which occurs during normal brain development via the action of miR-124, is sufficient to induce trans-differentiation of fibroblasts into functional neurons. Besides its traditional role in regulated splicing, we show that PTB has a previously undocumented function in the regulation of microRNA functions, suppressing or enhancing microRNA targeting by competitive binding on target mRNA or altering local RNA secondary structure. A key event during neuronal induction is the relief of PTB-mediated blockage of microRNA action on multiple components of the REST complex, thereby derepressing a large array of neuronal genes, including miR-124 and multiple neuronal-specific transcription factors, in nonneuronal cells. This converts a negative feedback loop to a positive one to elicit cellular reprogramming to the neuronal lineage.

Regulation of the Hippo-YAP pathway by G-protein-coupled receptor signaling.

The Hippo pathway is crucial in organ size control, and its dysregulation contributes to tumorigenesis. However, upstream signals that regulate the mammalian Hippo pathway have remained elusive. Here, we report that the Hippo pathway is regulated by G-protein-coupled receptor (GPCR) signaling. Serum-borne lysophosphatidic acid (LPA) and sphingosine 1-phosphophate (S1P) act through G12/13-coupled receptors to inhibit the Hippo pathway kinases Lats1/2, thereby activating YAP and TAZ transcription coactivators, which are oncoproteins repressed by Lats1/2. YAP and TAZ are involved in LPA-induced gene expression, cell migration, and proliferation. In contrast, stimulation of Gs-coupled receptors by glucagon or epinephrine activates Lats1/2 kinase activity, thereby inhibiting YAP function. Thus, GPCR signaling can either activate or inhibit the Hippo-YAP pathway depending on the coupled G protein. Our study identifies extracellular diffusible signals that modulate the Hippo pathway and also establishes the Hippo-YAP pathway as a critical signaling branch downstream of GPCR.

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.

The Augmented R-Loop Is a Unifying Mechanism for Myelodysplastic Syndromes Induced by High-Risk Splicing Factor Mutations.

Mutations in several general pre-mRNA splicing factors have been linked to myelodysplastic syndromes (MDSs) and solid tumors. These mutations have generally been assumed to cause disease by the resultant splicing defects, but different mutations appear to induce distinct splicing defects, raising the possibility that an alternative common mechanism is involved. Here we report a chain of events triggered by multiple splicing factor mutations, especially high-risk alleles in SRSF2 and U2AF1, including elevated R-loops, replication stress, and activation of the ataxia telangiectasia and Rad3-related protein (ATR)-Chk1 pathway. We further demonstrate that enhanced R-loops, opposite to the expectation from gained RNA binding with mutant SRSF2, result from impaired transcription pause release because the mutant protein loses its ability to extract the RNA polymerase II (Pol II) C-terminal domain (CTD) kinase-the positive transcription elongation factor complex (P-TEFb)-from the 7SK complex. Enhanced R-loops are linked to compromised proliferation of bone-marrow-derived blood progenitors, which can be partially rescued by RNase H overexpression, suggesting a direct contribution of augmented R-loops to the MDS phenotype.

R-ChIP Using Inactive RNase H Reveals Dynamic Coupling of R-loops with Transcriptional Pausing at Gene Promoters.

R-loop, a three-stranded RNA/DNA structure, has been linked to induced genome instability and regulated gene expression. To enable precision analysis of R-loops in vivo, we develop an RNase-H-based approach; this reveals predominant R-loop formation near gene promoters with strong G/C skew and propensity to form G-quadruplex in non-template DNA, corroborating with all biochemically established properties of R-loops. Transcription perturbation experiments further indicate that R-loop induction correlates to transcriptional pausing. Interestingly, we note that most mapped R-loops are each linked to a nearby free RNA end; by using a ribozyme to co-transcriptionally cleave nascent RNA, we demonstrate that such a free RNA end coupled with a G/C-skewed sequence is necessary and sufficient to induce R-loop. These findings provide a topological solution for RNA invasion into duplex DNA and suggest an order for R-loop initiation and elongation in an opposite direction to that previously proposed.

RBFox2 Binds Nascent RNA to Globally Regulate Polycomb Complex 2 Targeting in Mammalian Genomes.

Increasing evidence suggests that diverse RNA binding proteins (RBPs) interact with regulatory RNAs to regulate transcription. RBFox2 is a well-characterized pre-mRNA splicing regulator, but we now encounter an unexpected paradigm where depletion of this RBP induces widespread increase in nascent RNA production in diverse cell types. Chromatin immunoprecipitation sequencing (ChIP-seq) reveals extensive interaction of RBFox2 with chromatin in a nascent RNA-dependent manner. Bayesian network analysis connects RBFox2 to Polycomb complex 2 (PRC2) and H3K27me3, and biochemical experiments demonstrate the ability of RBFox2 to directly interact with PRC2. Strikingly, RBFox2 inactivation eradicates PRC2 targeting on the majority of bivalent gene promoters and leads to transcriptional de-repression. Together, these findings uncover a mechanism underlying the enigmatic association of PRC2 with numerous active genes, highlight the importance of gene body sequences to gauge transcriptional output, and suggest nascent RNAs as critical signals for transcriptional feedback control to maintain homeostatic gene expression in mammalian genomes.

Context-dependent modulation of Pol II CTD phosphatase SSUP-72 regulates alternative polyadenylation in neuronal development.

Alternative polyadenylation (APA) is widespread in neuronal development and activity-mediated neural plasticity. However, the underlying molecular mechanisms are largely unknown. We used systematic genetic studies and genome-wide surveys of the transcriptional landscape to identify a context-dependent regulatory pathway controlling APA in the Caenorhabditis elegans nervous system. Loss of function in ssup-72, a Ser5 phosphatase for the RNA polymerase II (Pol II) C-terminal domain (CTD), dampens transcription termination at a strong intronic polyadenylation site (PAS) in unc-44/ankyrin yet promotes termination at the weak intronic PAS of the MAP kinase dlk-1. A nuclear protein, SYDN-1, which regulates neuronal development, antagonizes the function of SSUP-72 and several nuclear polyadenylation factors. This regulatory pathway allows the production of a neuron-specific isoform of unc-44 and an inhibitory isoform of dlk-1. Dysregulation of the unc-44 and dlk-1 mRNA isoforms in sydn-1 mutants impairs neuronal development. Deleting the intronic PAS of unc-44 results in increased pre-mRNA processing of neuronal ankyrin and suppresses sydn-1 mutants. These results reveal a mechanism by which regulation of CTD phosphorylation controls coding region APA in the nervous system.

An evolutionarily conserved DNA architecture determines target specificity of the TWIST family bHLH transcription factors.

Basic helix-loop-helix (bHLH) transcription factors recognize the canonical E-box (CANNTG) to regulate gene transcription; however, given the prevalence of E-boxes in a genome, it has been puzzling how individual bHLH proteins selectively recognize E-box sequences on their targets. TWIST is a bHLH transcription factor that promotes epithelial-mesenchymal transition (EMT) during development and tumor metastasis. High-resolution mapping of TWIST occupancy in human and Drosophila genomes reveals that TWIST, but not other bHLH proteins, recognizes a unique double E-box motif with two E-boxes spaced preferentially by 5 nucleotides. Using molecular modeling and binding kinetic analyses, we found that the strict spatial configuration in the double E-box motif aligns two TWIST-E47 dimers on the same face of DNA, thus providing a high-affinity site for a highly stable intramolecular tetramer. Biochemical analyses showed that the WR domain of TWIST dimerizes to mediate tetramer formation, which is functionally required for TWIST-induced EMT. These results uncover a novel mechanism for a bHLH transcription factor to recognize a unique spatial configuration of E-boxes to achieve target specificity. The WR-WR domain interaction uncovered here sets an example of target gene specificity of a bHLH protein being controlled allosterically by a domain outside of the bHLH region.

A tumorigenic index for quantitative analysis of liver cancer initiation and progression.

Primary liver cancer develops from multifactorial etiologies, resulting in extensive genomic heterogeneity. To probe the common mechanism of hepatocarcinogenesis, we interrogated temporal gene expression profiles in a group of mouse models with hepatic steatosis, fibrosis, inflammation, and, consequently, tumorigenesis. Instead of anticipated progressive changes, we observed a sudden molecular switch at a critical precancer stage, by developing analytical platform that focuses on transcription factor (TF) clusters. Coarse-grained network modeling demonstrated that an abrupt transcriptomic transition occurred once changes were accumulated to reach a threshold. Based on the experimental and bioinformatic data analyses as well as mathematical modeling, we derived a tumorigenic index (TI) to quantify tumorigenic signal strengths. The TI is powerful in predicting the disease status of patients with metabolic disorders and also the tumor stages and prognosis of liver cancer patients with diverse backgrounds. This work establishes a quantitative tool for triage of liver cancer patients and also for cancer risk assessment of chronic liver disease patients.

Molecular basis for 5-carboxycytosine recognition by RNA polymerase II elongation complex.

DNA methylation at selective cytosine residues (5-methylcytosine (5mC)) and their removal by TET-mediated DNA demethylation are critical for setting up pluripotent states in early embryonic development. TET enzymes successively convert 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC), with 5fC and 5caC subject to removal by thymine DNA glycosylase (TDG) in conjunction with base excision repair. Early reports indicate that 5fC and 5caC could be stably detected on enhancers, promoters and gene bodies, with distinct effects on gene expression, but the mechanisms have remained elusive. Here we determined the X-ray crystal structure of yeast elongating RNA polymerase II (Pol II) in complex with a DNA template containing oxidized 5mCs, revealing specific hydrogen bonds between the 5-carboxyl group of 5caC and the conserved epi-DNA recognition loop in the polymerase. This causes a positional shift for incoming nucleoside 5'-triphosphate (NTP), thus compromising nucleotide addition. To test the implication of this structural insight in vivo, we determined the global effect of increased 5fC/5caC levels on transcription, finding that such DNA modifications indeed retarded Pol II elongation on gene bodies. These results demonstrate the functional impact of oxidized 5mCs on gene expression and suggest a novel role for Pol II as a specific and direct epigenetic sensor during transcription elongation.

Genome-wide analysis reveals SR protein cooperation and competition in regulated splicing.

SR proteins are well-characterized RNA binding proteins that promote exon inclusion by binding to exonic splicing enhancers (ESEs). However, it has been unclear whether regulatory rules deduced on model genes apply generally to activities of SR proteins in the cell. Here, we report global analyses of two prototypical SR proteins, SRSF1 (SF2/ASF) and SRSF2 (SC35), using splicing-sensitive arrays and CLIP-seq on mouse embryo fibroblasts (MEFs). Unexpectedly, we find that these SR proteins promote both inclusion and skipping of exons in vivo, but their binding patterns do not explain such opposite responses. Further analyses reveal that loss of one SR protein is accompanied by coordinated loss or compensatory gain in the interaction of other SR proteins at the affected exons. Therefore, specific effects on regulated splicing by one SR protein actually depend on a complex set of relationships with multiple other SR proteins in mammalian genomes.

Toxic gain of function from mutant FUS protein is crucial to trigger cell autonomous motor neuron loss.

FUS is an RNA-binding protein involved in amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Cytoplasmic FUS-containing aggregates are often associated with concomitant loss of nuclear FUS Whether loss of nuclear FUS function, gain of a cytoplasmic function, or a combination of both lead to neurodegeneration remains elusive. To address this question, we generated knockin mice expressing mislocalized cytoplasmic FUS and complete FUS knockout mice. Both mouse models display similar perinatal lethality with respiratory insufficiency, reduced body weight and length, and largely similar alterations in gene expression and mRNA splicing patterns, indicating that mislocalized FUS results in loss of its normal function. However, FUS knockin mice, but not FUS knockout mice, display reduced motor neuron numbers at birth, associated with enhanced motor neuron apoptosis, which can be rescued by cell-specific CRE-mediated expression of wild-type FUS within motor neurons. Together, our findings indicate that cytoplasmic FUS mislocalization not only leads to nuclear loss of function, but also triggers motor neuron death through a toxic gain of function within motor neurons.

WNT7A and PAX6 define corneal epithelium homeostasis and pathogenesis.

The surface of the cornea consists of a unique type of non-keratinized epithelial cells arranged in an orderly fashion, and this is essential for vision by maintaining transparency for light transmission. Cornea epithelial cells (CECs) undergo continuous renewal from limbal stem or progenitor cells (LSCs), and deficiency in LSCs or corneal epithelium--which turns cornea into a non-transparent, keratinized skin-like epithelium--causes corneal surface disease that leads to blindness in millions of people worldwide. How LSCs are maintained and differentiated into corneal epithelium in healthy individuals and which key molecular events are defective in patients have been largely unknown. Here we report establishment of an in vitro feeder-cell-free LSC expansion and three-dimensional corneal differentiation protocol in which we found that the transcription factors p63 (tumour protein 63) and PAX6 (paired box protein PAX6) act together to specify LSCs, and WNT7A controls corneal epithelium differentiation through PAX6. Loss of WNT7A or PAX6 induces LSCs into skin-like epithelium, a critical defect tightly linked to common human corneal diseases. Notably, transduction of PAX6 in skin epithelial stem cells is sufficient to convert them to LSC-like cells, and upon transplantation onto eyes in a rabbit corneal injury model, these reprogrammed cells are able to replenish CECs and repair damaged corneal surface. These findings suggest a central role of the WNT7A-PAX6 axis in corneal epithelial cell fate determination, and point to a new strategy for treating corneal surface diseases.

Nuclear matrix factor hnRNP U/SAF-A exerts a global control of alternative splicing by regulating U2 snRNP maturation.

The nuclear matrix-associated hnRNP U/SAF-A protein has been implicated in diverse pathways from transcriptional regulation to telomere length control to X inactivation, but the precise mechanism underlying each of these processes has remained elusive. Here, we report hnRNP U as a regulator of SMN2 splicing from a custom RNAi screen. Genome-wide analysis by CLIP-seq reveals that hnRNP U binds virtually to all classes of regulatory noncoding RNAs, including all snRNAs required for splicing of both major and minor classes of introns, leading to the discovery that hnRNP U regulates U2 snRNP maturation and Cajal body morphology in the nucleus. Global analysis of hnRNP U-dependent splicing by RNA-seq coupled with bioinformatic analysis of associated splicing signals suggests a general rule for splice site selection through modulating the core splicing machinery. These findings exemplify hnRNP U/SAF-A as a potent regulator of nuclear ribonucleoprotein particles in diverse gene expression pathways.