Human spliceosomal snRNA sequence variants generate variant spliceosomes.

Human pre-mRNA splicing is primarily catalyzed by the major spliceosome, comprising five small nuclear ribonucleoprotein complexes, U1, U2, U4, U5, and U6 snRNPs, each of which contains the corresponding U-rich snRNA. These snRNAs are encoded by large gene families exhibiting significant sequence variation, but it remains unknown if most human snRNA genes are untranscribed pseudogenes or produce variant snRNAs with the potential to differentially influence splicing. Since gene duplication and variation are powerful mechanisms of evolutionary adaptation, we sought to address this knowledge gap by systematically profiling human U1, U2, U4, and U5 snRNA variant gene transcripts. We identified 55 transcripts that are detectably expressed in human cells, 38 of which incorporate into snRNPs and spliceosomes in 293T cells. All U1 snRNA variants are more than 1000-fold less abundant in spliceosomes than the canonical U1, whereas at least 1% of spliceosomes contain a variant of U2 or U4. In contrast, eight U5 snRNA sequence variants occupy spliceosomes at levels of 1% to 46%. Furthermore, snRNA variants display distinct expression patterns across five human cell lines and adult and fetal tissues. Different RNA degradation rates contribute to the diverse steady state levels of snRNA variants. Our findings suggest that variant spliceosomes containing noncanonical snRNAs may contribute to different tissue- and cell-type-specific alternative splicing patterns.

Alternative Splicing and Cleavage of GLUT8.

The GLUT (SLC2) family of membrane-associated transporters are described as glucose transporters. However, this family is divided into three classes and, though the regulated transporter activity of class I proteins is becoming better understood, class III protein functions continue to be obscure. We have cataloged the relative expression and splicing of SLC2 mRNA isomers in tumors and normal tissues, with a focus on breast tumors and cell lines. mRNA for the class III protein GLUT8 is the predominant SLC2 species expressed alongside GLUT1 in many tissues, but GLUT8 mRNA exists mostly as an untranslated splice form in tumors. We confirm that GLUT8 is not presented at the cell surface and does not transport glucose directly. However, we reveal a lysosome-dependent reaction that cleaves the GLUT8 protein and releases the carboxy-terminal peptide to a separate vesicle population. Given the localization of GLUT8 at a major metabolic hub (the late endosomal/lysosomal interface) and its regulated cleavage reaction, we evaluated TXNIP-mediated hexosamine homeostasis and speculate that GLUT8 may function as a sensory component of this reaction.

Coordinated alterations in RNA splicing and epigenetic regulation drive leukaemogenesis.

Transcription and pre-mRNA splicing are key steps in the control of gene expression and mutations in genes regulating each of these processes are common in leukaemia1,2. Despite the frequent overlap of mutations affecting epigenetic regulation and splicing in leukaemia, how these processes influence one another to promote leukaemogenesis is not understood and, to our knowledge, there is no functional evidence that mutations in RNA splicing factors initiate leukaemia. Here, through analyses of transcriptomes from 982 patients with acute myeloid leukaemia, we identified frequent overlap of mutations in IDH2 and SRSF2 that together promote leukaemogenesis through coordinated effects on the epigenome and RNA splicing. Whereas mutations in either IDH2 or SRSF2 imparted distinct splicing changes, co-expression of mutant IDH2 altered the splicing effects of mutant SRSF2 and resulted in more profound splicing changes than either mutation alone. Consistent with this, co-expression of mutant IDH2 and SRSF2 resulted in lethal myelodysplasia with proliferative features in vivo and enhanced self-renewal in a manner not observed with either mutation alone. IDH2 and SRSF2 double-mutant cells exhibited aberrant splicing and reduced expression of INTS3, a member of the integrator complex3, concordant with increased stalling of RNA polymerase II (RNAPII). Aberrant INTS3 splicing contributed to leukaemogenesis in concert with mutant IDH2 and was dependent on mutant SRSF2 binding to cis elements in INTS3 mRNA and increased DNA methylation of INTS3. These data identify a pathogenic crosstalk between altered epigenetic state and splicing in a subset of leukaemias, provide functional evidence that mutations in splicing factors drive myeloid malignancy development, and identify spliceosomal changes as a mediator of IDH2-mutant leukaemogenesis.

RNA components of the spliceosome regulate tissue- and cancer-specific alternative splicing.

Alternative splicing of pre-mRNAs plays a pivotal role during the establishment and maintenance of human cell types. Characterizing the trans-acting regulatory proteins that control alternative splicing has therefore been the focus of much research. Recent work has established that even core protein components of the spliceosome, which are required for splicing to proceed, can nonetheless contribute to splicing regulation by modulating splice site choice. We here show that the RNA components of the spliceosome likewise influence alternative splicing decisions. Although these small nuclear RNAs (snRNAs), termed U1, U2, U4, U5, and U6 snRNA, are present in equal stoichiometry within the spliceosome, we found that their relative levels vary by an order of magnitude during development, across tissues, and across cancer samples. Physiologically relevant perturbation of individual snRNAs drove widespread gene-specific differences in alternative splicing but not transcriptome-wide splicing failure. Genes that were particularly sensitive to variations in snRNA abundance in a breast cancer cell line model were likewise preferentially misspliced within a clinically diverse cohort of invasive breast ductal carcinomas. As aberrant mRNA splicing is prevalent in many cancers, we propose that a full understanding of such dysregulated pre-mRNA processing requires study of snRNAs, as well as protein splicing factors. Together, our data show that the RNA components of the spliceosome are not merely basal factors, as has long been assumed. Instead, these noncoding RNAs constitute a previously uncharacterized layer of regulation of alternative splicing, and contribute to the establishment of global splicing programs in both healthy and malignant cells.

Regulation of alternative mRNA splicing: old players and new perspectives.

Nearly all human multiexon genes are subject to alternative splicing in one or more cell types. The splicing machinery therefore has to select between multiple splice sites in a context-dependent manner, relying on sequence features in cis- and trans-acting splicing regulators that either promote or repress splice site recognition and spliceosome assembly. However, the functional coupling between multiple gene regulatory layers signifies that splicing can also be modulated by transcriptional or epigenetic characteristics. Other, less obvious, aspects of alternative splicing have come to light in recent years, often involving core components of the spliceosome previously thought to perform a basal rather than a regulatory role in splicing. Together this paints a highly dynamic picture of splicing regulation, where the final splice site choice is governed by the entire transcriptional environment of a gene and its cellular context.

RNA splicing factors as oncoproteins and tumour suppressors.

The recent genomic characterization of cancers has revealed recurrent somatic point mutations and copy number changes affecting genes encoding RNA splicing factors. Initial studies of these 'spliceosomal mutations' suggest that the proteins bearing these mutations exhibit altered splice site and/or exon recognition preferences relative to their wild-type counterparts, resulting in cancer-specific mis-splicing. Such changes in the splicing machinery may create novel vulnerabilities in cancer cells that can be therapeutically exploited using compounds that can influence the splicing process. Further studies to dissect the biochemical, genomic and biological effects of spliceosomal mutations are crucial for the development of cancer therapies targeted at these mutations.

Modulation of splicing catalysis for therapeutic targeting of leukemia with mutations in genes encoding spliceosomal proteins.

Mutations in genes encoding splicing factors (which we refer to as spliceosomal genes) are commonly found in patients with myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML). These mutations recurrently affect specific amino acid residues, leading to perturbed normal splice site and exon recognition. Spliceosomal gene mutations are always heterozygous and rarely occur together with one another, suggesting that cells may tolerate only a partial deviation from normal splicing activity. To test this hypothesis, we engineered mice to express a mutated allele of serine/arginine-rich splicing factor 2 (Srsf2(P95H))-which commonly occurs in individuals with MDS and AML-in an inducible, hemizygous manner in hematopoietic cells. These mice rapidly succumbed to fatal bone marrow failure, demonstrating that Srsf2-mutated cells depend on the wild-type Srsf2 allele for survival. In the context of leukemia, treatment with the spliceosome inhibitor E7107 (refs. 7,8) resulted in substantial reductions in leukemic burden, specifically in isogenic mouse leukemias and patient-derived xenograft AMLs carrying spliceosomal mutations. Whereas E7107 treatment of mice resulted in widespread intron retention and cassette exon skipping in leukemic cells regardless of Srsf2 genotype, the magnitude of splicing inhibition following E7107 treatment was greater in Srsf2-mutated than in Srsf2-wild-type leukemia, consistent with the differential effect of E7107 on survival. Collectively, these data provide genetic and pharmacologic evidence that leukemias with spliceosomal gene mutations are preferentially susceptible to additional splicing perturbations in vivo as compared to leukemias without such mutations. Modulation of spliceosome function may thus provide a new therapeutic avenue in genetically defined subsets of individuals with MDS or AML.

Expression of androgen receptor splice variants in clinical breast cancers.

The importance of androgen receptor (AR) signaling is increasingly being recognized in breast cancer, which has elicited clinical trials aimed at assessing the efficacy of androgen deprivation therapy (ADT) for metastatic disease. In prostate cancer, resistance to ADT is frequently associated with the emergence of androgen-independent splice variants of the AR (AR variants, AR-Vs) that lack the LBD and are constitutively active. Women with breast cancer may be prone to a similar phenomenon. Herein, we show that in addition to the prototypical transcript, the AR gene produces a diverse range of AR-V transcripts in primary breast tumors. The most frequently and highly expressed variant was AR-V7 (exons 1/2/3/CE3), which was detectable at the mRNA level in > 50% of all breast cancers and at the protein level in a subset of ERα-negative tumors. Functionally, AR-V7 is a constitutively active and ADT-resistant transcription factor that promotes growth and regulates a transcriptional program distinct from AR in ERα-negative breast cancer cells. Importantly, we provide ex vivo evidence that AR-V7 is upregulated by the AR antagonist enzalutamide in primary breast tumors. These findings have implications for treatment response in the ongoing clinical trials of ADT in breast cancer.

Widespread intron retention diversifies most cancer transcriptomes.

BACKGROUND: Somatic mutations affecting components of the RNA splicing machinery occur with high frequencies across many tumor types. These mutations give rise to distinct alterations in normal splice site and exon recognition, such as unusual 3' splice site preferences, that likely contribute to tumorigenesis. METHODS: We analyzed genome-wide patterns of RNA splicing across 805 matched tumor and normal control samples from 16 distinct cancer types to identify signals of abnormal cancer-associated splicing. RESULTS: We found that abnormal RNA splicing, typified by widespread intron retention, is common across cancers even in the absence of mutations directly affecting the RNA splicing machinery. Almost all liquid and solid cancer types exhibited frequent retention of both alternative and constitutive introns relative to control normal tissues. The sole exception was breast cancer, where intron retention typified adjacent normal rather than cancer tissue. Different introns were preferentially retained in specific cancer types, although a small subset of introns enriched for genes encoding RNA splicing and export factors exhibited frequent retention across diverse cancers. The extent of intron retention correlated with the presence of IDH1 and IDH2 mutations in acute myeloid leukemia and across molecular subtypes in breast cancer. Many introns that were preferentially retained in primary cancers were present at high levels in the cytoplasmic mRNA pools of cancer cell lines. CONCLUSIONS: Our data indicate that abnormal RNA splicing is a common characteristic of cancers even in the absence of mutational insults to the splicing machinery, and suggest that intron-containing mRNAs contribute to the transcriptional diversity of many cancers.

Integrative clinical genomics of advanced prostate cancer.

Toward development of a precision medicine framework for metastatic, castration-resistant prostate cancer (mCRPC), we established a multi-institutional clinical sequencing infrastructure to conduct prospective whole-exome and transcriptome sequencing of bone or soft tissue tumor biopsies from a cohort of 150 mCRPC affected individuals. Aberrations of AR, ETS genes, TP53, and PTEN were frequent (40%-60% of cases), with TP53 and AR alterations enriched in mCRPC compared to primary prostate cancer. We identified new genomic alterations in PIK3CA/B, R-spondin, BRAF/RAF1, APC, ß-catenin, and ZBTB16/PLZF. Moreover, aberrations of BRCA2, BRCA1, and ATM were observed at substantially higher frequencies (19.3% overall) compared to those in primary prostate cancers. 89% of affected individuals harbored a clinically actionable aberration, including 62.7% with aberrations in AR, 65% in other cancer-related genes, and 8% with actionable pathogenic germline alterations. This cohort study provides clinically actionable information that could impact treatment decisions for these affected individuals.

Sample processing obscures cancer-specific alterations in leukemic transcriptomes.

Substantial effort is currently devoted to identifying cancer-associated alterations using genomics. Here, we show that standard blood collection procedures rapidly change the transcriptional and posttranscriptional landscapes of hematopoietic cells, resulting in biased activation of specific biological pathways; up-regulation of pseudogenes, antisense RNAs, and unannotated coding isoforms; and RNA surveillance inhibition. Affected genes include common mutational targets and thousands of other genes participating in processes such as chromatin modification, RNA splicing, T- and B-cell activation, and NF-κB signaling. The majority of published leukemic transcriptomes exhibit signals of this incubation-induced dysregulation, explaining up to 40% of differences in gene expression and alternative splicing between leukemias and reference normal transcriptomes. The effects of sample processing are particularly evident in pan-cancer analyses. We provide biomarkers that detect prolonged incubation of individual samples and show that keeping blood on ice markedly reduces changes to the transcriptome. In addition to highlighting the potentially confounding effects of technical artifacts in cancer genomics data, our study emphasizes the need to survey the diversity of normal as well as neoplastic cells when characterizing tumors.

The shaping and functional consequences of the microRNA landscape in breast cancer.

MicroRNAs (miRNAs) show differential expression across breast cancer subtypes, and have both oncogenic and tumour-suppressive roles. Here we report the miRNA expression profiles of 1,302 breast tumours with matching detailed clinical annotation, long-term follow-up and genomic and messenger RNA expression data. This provides a comprehensive overview of the quantity, distribution and variation of the miRNA population and provides information on the extent to which genomic, transcriptional and post-transcriptional events contribute to miRNA expression architecture, suggesting an important role for post-transcriptional regulation. The key clinical parameters and cellular pathways related to the miRNA landscape are characterized, revealing context-dependent interactions, for example with regards to cell adhesion and Wnt signalling. Notably, only prognostic miRNA signatures derived from breast tumours devoid of somatic copy-number aberrations (CNA-devoid) are consistently prognostic across several other subtypes and can be validated in external cohorts. We then use a data-driven approach to seek the effects of miRNAs associated with differential co-expression of mRNAs, and find that miRNAs act as modulators of mRNA-mRNA interactions rather than as on-off molecular switches. We demonstrate such an important modulatory role for miRNAs in the biology of CNA-devoid breast cancers, a common subtype in which the immune response is prominent. These findings represent a new framework for studying the biology of miRNAs in human breast cancer.

Genome-wide characterization of Foxa2 targets reveals upregulation of floor plate genes and repression of ventrolateral genes in midbrain dopaminergic progenitors.

The transcription factors Foxa1 and Foxa2 promote the specification of midbrain dopaminergic (mDA) neurons and the floor plate. Whether their role is direct has remained unclear as they also regulate the expression of Shh, which has similar roles. We characterized the Foxa2 cis-regulatory network by chromatin immunoprecipitation followed by high-throughput sequencing of mDA progenitors. This identified 9160 high-quality Foxa2 binding sites associated with 5409 genes, providing mechanistic insights into Foxa2-mediated positive and negative regulatory events. Foxa2 regulates directly and positively key determinants of mDA neurons, including Lmx1a, Lmx1b, Msx1 and Ferd3l, while negatively inhibiting transcription factors expressed in ventrolateral midbrain such as Helt, Tle4, Otx1, Sox1 and Tal2. Furthermore, Foxa2 negatively regulates extrinsic and intrinsic components of the Shh signaling pathway, possibly by binding to the same enhancer regions of co-regulated genes as Gli1. Foxa2 also regulates the expression of floor plate factors that control axon trajectories around the midline of the embryo, thereby contributing to the axon guidance function of the floor plate. Finally, this study identified multiple Foxa2-regulated enhancers that are active in the floor plate of the midbrain or along the length of the embryo in mouse and chick. This work represents the first comprehensive characterization of Foxa2 targets in mDA progenitors and provides a framework for elaborating gene regulatory networks in a functionally important progenitor population.

NuRD-mediated deacetylation of H3K27 facilitates recruitment of Polycomb Repressive Complex 2 to direct gene repression.

Pluripotent cells possess the ability to differentiate into any cell type. Commitment to differentiate into specific lineages requires strict control of gene expression to coordinate the downregulation of lineage inappropriate genes while enabling the expression of lineage-specific genes. The nucleosome remodelling and deacetylation complex (NuRD) is required for lineage commitment of pluripotent cells; however, the mechanism through which it exerts this effect has not been defined. Here, we show that histone deacetylation by NuRD specifies recruitment for Polycomb Repressive Complex 2 (PRC2) in embryonic stem (ES) cells. NuRD-mediated deacetylation of histone H3K27 enables PRC2 recruitment and subsequent H3K27 trimethylation at NuRD target promoters. We propose a gene-specific mechanism for modulating expression of transcriptionally poised genes whereby NuRD controls the balance between acetylation and methylation of histones, thereby precisely directing the expression of genes critical for embryonic development.

PeakAnalyzer: genome-wide annotation of chromatin binding and modification loci.

BACKGROUND: Functional genomic studies involving high-throughput sequencing and tiling array applications, such as ChIP-seq and ChIP-chip, generate large numbers of experimentally-derived signal peaks across the genome under study. In analyzing these loci to determine their potential regulatory functions, areas of signal enrichment must be considered relative to proximal genes and regulatory elements annotated throughout the target genome Regions of chromatin association by transcriptional regulators should be distinguished as individual binding sites in order to enhance downstream analyses, such as the identification of known and novel consensus motifs. RESULTS: PeakAnalyzer is a set of high-performance utilities for the automated processing of experimentally-derived peak regions and annotation of genomic loci. The programs can accurately subdivide multimodal regions of signal enrichment into distinct subpeaks corresponding to binding sites or chromatin modifications, retrieve genomic sequences encompassing the computed subpeak summits, and identify positional features of interest such as intersection with exon/intron gene components, proximity to up- or downstream transcriptional start sites and cis-regulatory elements. The software can be configured to run either as a pipeline component for high-throughput analyses, or as a cross-platform desktop application with an intuitive user interface. CONCLUSIONS: PeakAnalyzer comprises a number of utilities essential for ChIP-seq and ChIP-chip data analysis. High-performance implementations are provided for Unix pipeline integration along with a GUI version for interactive use. Source code in C++ and Java is provided, as are native binaries for Linux, Mac OS X and Windows systems.

Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression.

RNA abundance and DNA copy number are routinely measured in high-throughput using microarray and next-generation sequencing (NGS) technologies, and the attributes of different platforms have been extensively analyzed. Recently, the application of both microarrays and NGS has expanded to include microRNAs (miRNAs), but the relative performance of these methods has not been rigorously characterized. We analyzed three biological samples across six miRNA microarray platforms and compared their hybridization performance. We examined the utility of these platforms, as well as NGS, for the detection of differentially expressed miRNAs. We then validated the results for 89 miRNAs by real-time RT-PCR and challenged the use of this assay as a "gold standard." Finally, we implemented a novel method to evaluate false-positive and false-negative rates for all methods in the absence of a reference method.

HTqPCR: high-throughput analysis and visualization of quantitative real-time PCR data in R.

MOTIVATION: Quantitative real-time polymerase chain reaction (qPCR) is routinely used for RNA expression profiling, validation of microarray hybridization data and clinical diagnostic assays. Although numerous statistical tools are available in the public domain for the analysis of microarray experiments, this is not the case for qPCR. Proprietary software is typically provided by instrument manufacturers, but these solutions are not amenable to the tandem analysis of multiple assays. This is problematic when an experiment involves more than a simple comparison between a control and treatment sample, or when many qPCR datasets are to be analyzed in a high-throughput facility. RESULTS: We have developed HTqPCR, a package for the R statistical computing environment, to enable the processing and analysis of qPCR data across multiple conditions and replicates. AVAILABILITY: HTqPCR and user documentation can be obtained through Bioconductor or at http://www.ebi.ac.uk/bertone/software. CONTACT: bertone@ebi.ac.uk

Mapping organelle proteins and protein complexes in Drosophila melanogaster.

Many proteins within eukaryotic cells are organized spatially and functionally into membrane bound organelles and complexes. A protein's location thus provides information about its function. Here, we apply LOPIT, a mass-spectrometry based technique that simultaneously maps proteins to specific subcellular compartments, to Drosophila embryos. We determine the subcellular distribution of hundreds of proteins, and protein complexes. Our results reveal the potential of LOPIT to provide average snapshots of cells.