Recruitment of the m6A/m6Am demethylase FTO to target RNAs by the telomeric zinc finger protein ZBTB48.

BACKGROUND: N6-methyladenosine (m6A), the most abundant internal modification on eukaryotic mRNA, and N6, 2'-O-dimethyladenosine (m6Am), are epitranscriptomic marks that function in multiple aspects of posttranscriptional regulation. Fat mass and obesity-associated protein (FTO) can remove both m6A and m6Am; however, little is known about how FTO achieves its substrate selectivity. RESULTS: Here, we demonstrate that ZBTB48, a C2H2-zinc finger protein that functions in telomere maintenance, associates with FTO and binds both mRNA and the telomere-associated regulatory RNA TERRA to regulate the functional interactions of FTO with target transcripts. Specifically, depletion of ZBTB48 affects targeting of FTO to sites of m6A/m6Am modification, changes cellular m6A/m6Am levels and, consequently, alters decay rates of target RNAs. ZBTB48 ablation also accelerates growth of HCT-116 colorectal cancer cells and modulates FTO-dependent regulation of Metastasis-associated protein 1 (MTA1) transcripts by controlling the binding to MTA1 mRNA of the m6A reader IGF2BP2. CONCLUSIONS: Our findings thus uncover a previously unknown mechanism of posttranscriptional regulation in which ZBTB48 co-ordinates RNA-binding of the m6A/m6Am demethylase FTO to control expression of its target RNAs.

Matrix stiffness-dependent regulation of immunomodulatory genes in human MSCs is associated with the lncRNA CYTOR.

Cell-matrix interactions in 3D environments significantly differ from those in 2D cultures. As such, mechanisms of mechanotransduction in 2D cultures are not necessarily applicable to cell-encapsulating hydrogels that resemble features of tissue architecture. Accordingly, the characterization of molecular pathways in 3D matrices is expected to uncover insights into how cells respond to their mechanical environment in physiological contexts, and potentially also inform hydrogel-based strategies in cell therapies. In this study, a bone marrow-mimetic hydrogel was employed to systematically investigate the stiffness-responsive transcriptome of mesenchymal stromal cells. High matrix rigidity impeded integrin-collagen adhesion, resulting in changes in cell morphology characterized by a contractile network of actin proximal to the cell membrane. This resulted in a suppression of extracellular matrix-regulatory genes involved in the remodeling of collagen fibrils, as well as the upregulation of secreted immunomodulatory factors. Moreover, an investigation of long noncoding RNAs revealed that the cytoskeleton regulator RNA (CYTOR) contributes to these 3D stiffness-driven changes in gene expression. Knockdown of CYTOR using antisense oligonucleotides enhanced the expression of numerous mechanoresponsive cytokines and chemokines to levels exceeding those achievable by modulating matrix stiffness alone. Taken together, our findings further our understanding of mechanisms of mechanotransduction that are distinct from canonical mechanotransductive pathways observed in 2D cultures.

High-throughput sensitive screening of small molecule modulators of microexon alternative splicing using dual Nano and Firefly luciferase reporters.

Disruption of alternative splicing frequently causes or contributes to human diseases and disorders. Consequently, there is a need for efficient and sensitive reporter assays capable of screening chemical libraries for compounds with efficacy in modulating important splicing events. Here, we describe a screening workflow employing dual Nano and Firefly luciferase alternative splicing reporters that affords efficient, sensitive, and linear detection of small molecule responses. Applying this system to a screen of ~95,000 small molecules identified compounds that stimulate or repress the splicing of neuronal microexons, a class of alternative exons often disrupted in autism and activated in neuroendocrine cancers. One of these compounds rescues the splicing of several analyzed microexons in the cerebral cortex of an autism mouse model haploinsufficient for Srrm4, a major activator of brain microexons. We thus describe a broadly applicable high-throughput screening system for identifying candidate splicing therapeutics, and a resource of small molecule modulators of microexons with potential for further development in correcting aberrant splicing patterns linked to human disorders and disease.

Global detection of human variants and isoforms by deep proteome sequencing.

An average shotgun proteomics experiment detects approximately 10,000 human proteins from a single sample. However, individual proteins are typically identified by peptide sequences representing a small fraction of their total amino acids. Hence, an average shotgun experiment fails to distinguish different protein variants and isoforms. Deeper proteome sequencing is therefore required for the global discovery of protein isoforms. Using six different human cell lines, six proteases, deep fractionation and three tandem mass spectrometry fragmentation methods, we identify a million unique peptides from 17,717 protein groups, with a median sequence coverage of approximately 80%. Direct comparison with RNA expression data provides evidence for the translation of most nonsynonymous variants. We have also hypothesized that undetected variants likely arise from mutation-induced protein instability. We further observe comparable detection rates for exon-exon junction peptides representing constitutive and alternative splicing events. Our dataset represents a resource for proteoform discovery and provides direct evidence that most frame-preserving alternatively spliced isoforms are translated.

Regulation of alternative polyadenylation by the C2H2-zinc-finger protein Sp1.

Alternative polyadenylation (APA) enhances gene regulatory potential by increasing the diversity of mRNA transcripts. 3' UTR shortening through APA correlates with enhanced cellular proliferation and is a widespread phenomenon in tumor cells. Here, we show that the ubiquitously expressed transcription factor Sp1 binds RNA in vivo and is a common repressor of distal poly(A) site usage. RNA sequencing identified 2,344 genes (36% of the total mapped mRNA transcripts) with lengthened 3' UTRs upon Sp1 depletion. Sp1 preferentially binds the 3' UTRs of such lengthened transcripts and inhibits cleavage at distal sites by interacting with the subunits of the core cleavage and polyadenylation (CPA) machinery. The 3' UTR lengths of Sp1 target genes in breast cancer patient RNA-seq data correlate with Sp1 expression levels, implicating Sp1-mediated APA regulation in modulating tumorigenic properties. Taken together, our findings provide insights into the mechanism for dynamic APA regulation by unraveling a previously unknown function of the DNA-binding transcription factor Sp1.

Autism-Misregulated eIF4G Microexons Control Synaptic Translation and Higher Order Cognitive Functions.

Microexons represent the most highly conserved class of alternative splicing, yet their functions are poorly understood. Here, we focus on closely related neuronal microexons overlapping prion-like domains in the translation initiation factors, eIF4G1 and eIF4G3, the splicing of which is activity dependent and frequently disrupted in autism. CRISPR-Cas9 deletion of these microexons selectively upregulates synaptic proteins that control neuronal activity and plasticity and further triggers a gene expression program mirroring that of activated neurons. Mice lacking the Eif4g1 microexon display social behavior, learning, and memory deficits, accompanied by altered hippocampal synaptic plasticity. We provide evidence that the eIF4G microexons function as a translational brake by causing ribosome stalling, through their propensity to promote the coalescence of cytoplasmic granule components associated with translation repression, including the fragile X mental retardation protein FMRP. The results thus reveal an autism-disrupted mechanism by which alternative splicing specializes neuronal translation to control higher order cognitive functioning.

Efficient and Accurate Quantitative Profiling of Alternative Splicing Patterns of Any Complexity on a Laptop.

Alternative splicing (AS) is a widespread process underlying the generation of transcriptomic and proteomic diversity and is frequently misregulated in human disease. Accordingly, an important goal of biomedical research is the development of tools capable of comprehensively, accurately, and efficiently profiling AS. Here, we describe Whippet, an easy-to-use RNA-seq analysis method that rapidly-with hardware requirements compatible with a laptop-models and quantifies AS events of any complexity without loss of accuracy. Using an entropic measure of splicing complexity, Whippet reveals that one-third of human protein coding genes produce transcripts with complex AS events involving co-expression of two or more principal splice isoforms. We observe that high-entropy AS events are more prevalent in tumor relative to matched normal tissues and correlate with increased expression of proto-oncogenic splicing factors. Whippet thus affords the rapid and accurate analysis of AS events of any complexity, and as such will facilitate future biomedical research.

Multilayered Control of Alternative Splicing Regulatory Networks by Transcription Factors.

Networks of coordinated alternative splicing (AS) events play critical roles in development and disease. However, a comprehensive knowledge of the factors that regulate these networks is lacking. We describe a high-throughput system for systematically linking trans-acting factors to endogenous RNA regulatory events. Using this system, we identify hundreds of factors associated with diverse regulatory layers that positively or negatively control AS events linked to cell fate. Remarkably, more than one-third of the regulators are transcription factors. Further analyses of the zinc finger protein Zfp871 and BTB/POZ domain transcription factor Nacc1, which regulate neural and stem cell AS programs, respectively, reveal roles in controlling the expression of specific splicing regulators. Surprisingly, these proteins also appear to regulate target AS programs via binding RNA. Our results thus uncover a large "missing cache" of splicing regulators among annotated transcription factors, some of which dually regulate AS through direct and indirect mechanisms.

Control of embryonic stem cell self-renewal and differentiation via coordinated alternative splicing and translation of YY2.

Translational control of gene expression plays a key role during the early phases of embryonic development. Here we describe a transcriptional regulator of mouse embryonic stem cells (mESCs), Yin-yang 2 (YY2), that is controlled by the translation inhibitors, Eukaryotic initiation factor 4E-binding proteins (4E-BPs). YY2 plays a critical role in regulating mESC functions through control of key pluripotency factors, including Octamer-binding protein 4 (Oct4) and Estrogen-related receptor-ß (Esrrb). Importantly, overexpression of YY2 directs the differentiation of mESCs into cardiovascular lineages. We show that the splicing regulator Polypyrimidine tract-binding protein 1 (PTBP1) promotes the retention of an intron in the 5'-UTR of Yy2 mRNA that confers sensitivity to 4E-BP-mediated translational suppression. Thus, we conclude that YY2 is a major regulator of mESC self-renewal and lineage commitment and document a multilayer regulatory mechanism that controls its expression.

IGF2BP3 Modulates the Interaction of Invasion-Associated Transcripts with RISC.

Insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3) expression correlates with malignancy, but its role(s) in pathogenesis remains enigmatic. We interrogated the IGF2BP3-RNA interaction network in pancreatic ductal adenocarcinoma (PDAC) cells. Using a combination of genome-wide approaches, we have identified 164 direct mRNA targets of IGF2BP3. These transcripts encode proteins enriched for functions such as cell migration, proliferation, and adhesion. Loss of IGF2BP3 reduced PDAC cell invasiveness and remodeled focal adhesion junctions. Individual nucleotide resolution crosslinking immunoprecipitation (iCLIP) revealed significant overlap of IGF2BP3 and microRNA (miRNA) binding sites. IGF2BP3 promotes association of the RNA-induced silencing complex (RISC) with specific transcripts. Our results show that IGF2BP3 influences a malignancy-associated RNA regulon by modulating miRNA-mRNA interactions.

An extensive program of periodic alternative splicing linked to cell cycle progression.

Progression through the mitotic cell cycle requires periodic regulation of gene function at the levels of transcription, translation, protein-protein interactions, post-translational modification and degradation. However, the role of alternative splicing (AS) in the temporal control of cell cycle is not well understood. By sequencing the human transcriptome through two continuous cell cycles, we identify ~1300 genes with cell cycle-dependent AS changes. These genes are significantly enriched in functions linked to cell cycle control, yet they do not significantly overlap genes subject to periodic changes in steady-state transcript levels. Many of the periodically spliced genes are controlled by the SR protein kinase CLK1, whose level undergoes cell cycle-dependent fluctuations via an auto-inhibitory circuit. Disruption of CLK1 causes pleiotropic cell cycle defects and loss of proliferation, whereas CLK1 over-expression is associated with various cancers. These results thus reveal a large program of CLK1-regulated periodic AS intimately associated with cell cycle control.

Compound heterozygous mutations in the noncoding RNU4ATAC cause Roifman Syndrome by disrupting minor intron splicing.

Roifman Syndrome is a rare congenital disorder characterized by growth retardation, cognitive delay, spondyloepiphyseal dysplasia and antibody deficiency. Here we utilize whole-genome sequencing of Roifman Syndrome patients to reveal compound heterozygous rare variants that disrupt highly conserved positions of the RNU4ATAC small nuclear RNA gene, a minor spliceosome component that is essential for minor intron splicing. Targeted sequencing confirms allele segregation in six cases from four unrelated families. RNU4ATAC rare variants have been recently reported to cause microcephalic osteodysplastic primordial dwarfism, type I (MOPD1), whose phenotype is distinct from Roifman Syndrome. Strikingly, all six of the Roifman Syndrome cases have one variant that overlaps MOPD1-implicated structural elements, while the other variant overlaps a highly conserved structural element not previously implicated in disease. RNA-seq analysis confirms extensive and specific defects of minor intron splicing. Available allele frequency data suggest that recessive genetic disorders caused by RNU4ATAC rare variants may be more prevalent than previously reported.

The alternative splicing factor Nova2 regulates vascular development and lumen formation.

Vascular lumen formation is a fundamental step during angiogenesis; yet, the molecular mechanisms underlying this process are poorly understood. Recent studies have shown that neural and vascular systems share common anatomical, functional and molecular similarities. Here we show that the organization of endothelial lumen is controlled at the post-transcriptional level by the alternative splicing (AS) regulator Nova2, which was previously considered to be neural cell-specific. Nova2 is expressed during angiogenesis and its depletion disrupts vascular lumen formation in vivo. Similarly, Nova2 depletion in cultured endothelial cells (ECs) impairs the apical distribution and the downstream signalling of the Par polarity complex, resulting in altered EC polarity, a process required for vascular lumen formation. These defects are linked to AS changes of Nova2 target exons affecting the Par complex and its regulators. Collectively, our results reveal that Nova2 functions as an AS regulator in angiogenesis and is a novel member of the 'angioneurins' family.

An alternative splicing event amplifies evolutionary differences between vertebrates.

Alternative splicing (AS) generates extensive transcriptomic and proteomic complexity. However, the functions of species- and lineage-specific splice variants are largely unknown. Here we show that mammalian-specific skipping of polypyrimidine tract-binding protein 1 (PTBP1) exon 9 alters the splicing regulatory activities of PTBP1 and affects the inclusion levels of numerous exons. During neurogenesis, skipping of exon 9 reduces PTBP1 repressive activity so as to facilitate activation of a brain-specific AS program. Engineered skipping of the orthologous exon in chicken cells induces a large number of mammalian-like AS changes in PTBP1 target exons. These results thus reveal that a single exon-skipping event in an RNA binding regulator directs numerous AS changes between species. Our results further suggest that these changes contributed to evolutionary differences in the formation of vertebrate nervous systems.

A germline mutation in SRRM2, a splicing factor gene, is implicated in papillary thyroid carcinoma predisposition.

Papillary thyroid carcinoma (PTC) displays strong but so far largely uncharacterized heritability. Here we studied genetic predisposition in a family with six affected individuals. We genotyped all available family members and conducted whole exome sequencing of blood DNA from two affected individuals. Haplotype analysis and other genetic criteria narrowed our list of candidates to a germline variant in the serine/arginine repetitive matrix 2 gene (SRRM2). This heterozygous variant, c.1037C > T (Ser346Phe or S346F; rs149019598) cosegregated with PTC in the family. It was not found in 138 other PTC families. It was found in 7/1,170 sporadic PTC cases and in 0/1,404 controls (p = 0.004). The encoded protein SRRM2 (also called SRm300) is part of the RNA splicing machinery. To evaluate the possibility that the S346F missense mutation affects alternative splicing, we compared RNA-Seq data in leukocytes from three mutation carriers and three controls. Significant differences in alternative splicing were identified for 1,642 exons, of which a subset of 7 exons was verified experimentally. The results confirmed a higher ratio of inclusion of exons in mutation carriers. These data suggest that the S346F mutation in SRRM2 predisposes to PTC by affecting alternative splicing of unidentified downstream target genes.

Myc and SAGA rewire an alternative splicing network during early somatic cell reprogramming.

Embryonic stem cells are maintained in a self-renewing and pluripotent state by multiple regulatory pathways. Pluripotent-specific transcriptional networks are sequentially reactivated as somatic cells reprogram to achieve pluripotency. How epigenetic regulators modulate this process and contribute to somatic cell reprogramming is not clear. Here we performed a functional RNAi screen to identify the earliest epigenetic regulators required for reprogramming. We identified components of the SAGA histone acetyltransferase complex, in particular Gcn5, as critical regulators of reprogramming initiation. Furthermore, we showed in mouse pluripotent stem cells that Gcn5 strongly associates with Myc and that, upon initiation of somatic reprogramming, Gcn5 and Myc form a positive feed-forward loop that activates a distinct alternative splicing network and the early acquisition of pluripotency-associated splicing events. These studies expose a Myc-SAGA pathway that drives expression of an essential alternative splicing regulatory network during somatic cell reprogramming.

RNA splicing. The human splicing code reveals new insights into the genetic determinants of disease.

To facilitate precision medicine and whole-genome annotation, we developed a machine-learning technique that scores how strongly genetic variants affect RNA splicing, whose alteration contributes to many diseases. Analysis of more than 650,000 intronic and exonic variants revealed widespread patterns of mutation-driven aberrant splicing. Intronic disease mutations that are more than 30 nucleotides from any splice site alter splicing nine times as often as common variants, and missense exonic disease mutations that have the least impact on protein function are five times as likely as others to alter splicing. We detected tens of thousands of disease-causing mutations, including those involved in cancers and spinal muscular atrophy. Examination of intronic and exonic variants found using whole-genome sequencing of individuals with autism revealed misspliced genes with neurodevelopmental phenotypes. Our approach provides evidence for causal variants and should enable new discoveries in precision medicine.

A global regulatory mechanism for activating an exon network required for neurogenesis.

The vertebrate and neural-specific Ser/Arg (SR)-related protein nSR100/SRRM4 regulates an extensive program of alternative splicing with critical roles in nervous system development. However, the mechanism by which nSR100 controls its target exons is poorly understood. We demonstrate that nSR100-dependent neural exons are associated with a unique configuration of intronic cis-elements that promote rapid switch-like regulation during neurogenesis. A key feature of this configuration is the insertion of specialized intronic enhancers between polypyrimidine tracts and acceptor sites that bind nSR100 to potently activate exon inclusion in neural cells while weakening 3' splice site recognition and contributing to exon skipping in nonneural cells. nSR100 further operates by forming multiple interactions with early spliceosome components bound proximal to 3' splice sites. These multifaceted interactions achieve dominance over neural exon silencing mediated by the splicing regulator PTBP1. The results thus illuminate a widespread mechanism by which a critical neural exon network is activated during neurogenesis.

Next-generation RNA sequencing of archival formalin-fixed paraffin-embedded urothelial bladder cancer.

UNLABELLED: Molecular profiling of individual cancers is key to personalised medicine. While sequencing technologies have required stringent sample collection and handling, recent technical advances offer sequencing from tissues collected in routine practice and tissues already stored in archives. In this paper, we establish methods for whole-transcriptome RNA sequencing (RNA-seq) from formalin-fixed paraffin-embedded tissues. We obtain average RNA-seq reads of >100 million per sample using the Illumina HiSeq2000 platform. We find high concordance with results from matching fresh frozen samples (>0.8 Spearman correlation). For validation, we compared low- and high-grade bladder cancer transcriptomes in 49 tumour samples after transurethral resection of bladder tumour. We found 947 differentially expressed protein-coding genes. While high-grade lesions exhibited distinct intertumour transcriptome heterogeneity, the transcriptome of low-grade tumours was homogeneous. PATIENT SUMMARY: In this report, we show that it is now possible to use universally available bladder cancer samples that have been fixed in formalin to perform high-quality transcriptome analysis. This ability will facilitate the development of transcriptome-wide tests based on gene expression correlated with clinical outcome.

Widespread intron retention in mammals functionally tunes transcriptomes.

Alternative splicing (AS) of precursor RNAs is responsible for greatly expanding the regulatory and functional capacity of eukaryotic genomes. Of the different classes of AS, intron retention (IR) is the least well understood. In plants and unicellular eukaryotes, IR is the most common form of AS, whereas in animals, it is thought to represent the least prevalent form. Using high-coverage poly(A)(+) RNA-seq data, we observe that IR is surprisingly frequent in mammals, affecting transcripts from as many as three-quarters of multiexonic genes. A highly correlated set of cis features comprising an "IR code" reliably discriminates retained from constitutively spliced introns. We show that IR acts widely to reduce the levels of transcripts that are less or not required for the physiology of the cell or tissue type in which they are detected. This "transcriptome tuning" function of IR acts through both nonsense-mediated mRNA decay and nuclear sequestration and turnover of IR transcripts. We further show that IR is linked to a cross-talk mechanism involving localized stalling of RNA polymerase II (Pol II) and reduced availability of spliceosomal components. Collectively, the results implicate a global checkpoint-type mechanism whereby reduced recruitment of splicing components coupled to Pol II pausing underlies widespread IR-mediated suppression of inappropriately expressed transcripts.