Widespread intron retention diversifies most cancer transcriptomes

Cajal body function in genome organization and transcriptome diversity

Nuclear bodies contribute to non-random organization of the human genome and nuclear function. Using a major prototypical nuclear body, the Cajal body, as an example, we suggest that these structures assemble at specific gene loci located across the genome as a result of high transcriptional activity. Subsequently, target genes are physically clustered in close proximity in Cajal body-containing cells. However, Cajal bodies are observed in only a limited number of human cell types, including neuronal and cancer cells. Ultimately, Cajal body depletion perturbs splicing kinetics by reducing target small nuclear RNA (snRNA) transcription and limiting the levels of spliceosomal snRNPs, including their modification and turnover following each round of RNA splicing. As such, Cajal bodies are capable of shaping the chromatin interaction landscape and the transcriptome by influencing spliceosome kinetics. Future studies should concentrate on characterizing the direct influence of Cajal bodies upon snRNA gene transcriptional dynamics. Also see the video abstract here.

An NXF1 mRNA with a retained intron is expressed in hippocampal and neocortical neurons and is translated into a protein that functions as an Nxf1 cofactor

A small Nxf1 protein, expressed from an NXF1 mRNA with a retained intron is highly expressed in rodent hippocampal and neocortical neurons, colocalizes with Staufen2 proteins in neuronal RNA granules, is present in polysomes, and replaces Nxt1 as an Nxf1 cofactor in export and expression of mRNA with retained introns.

The Nxf1 protein is a major nuclear export receptor for the transport of mRNA, and it also is essential for export of retroviral mRNAs with retained introns. In the latter case, it binds to RNA elements known as constitutive transport elements (CTEs) and functions in conjunction with a cofactor known as Nxt1. The NXF1 gene also regulates expression of its own intron-containing RNA through the use of a functional CTE within intron 10. mRNA containing this intron is exported to the cytoplasm, where it can be translated into the 356–amino acid short Nxf1(sNxf1) protein, despite the fact that it is a prime candidate for nonsense-mediated decay (NMD). Here we demonstrate that sNxf1 is highly expressed in nuclei and dendrites of hippocampal and neocortical neurons in rodent brain. Additionally, we show that sNxf1 localizes in RNA granules in neurites of differentiated N2a mouse neuroblastoma cells, where it shows partial colocalization with Staufen2 isoform SS, a protein known to play a role in dendritic mRNA trafficking. We also show that sNxf1 forms heterodimers in conjunction with the full-length Nxf1 and that sNxf1 can replace Nxt1 to enhance the expression of CTE-containing mRNA and promote its association with polyribosomes.

The bromodomain protein BRD4 regulates splicing during heat shock

The cellular response to heat stress is an ancient and evolutionarily highly conserved defence mechanism characterised by the transcriptional up-regulation of cyto-protective genes and a partial inhibition of splicing. These features closely resemble the proteotoxic stress response during tumor development. The bromodomain protein BRD4 has been identified as an integral member of the oxidative stress as well as of the inflammatory response, mainly due to its role in the transcriptional regulation process. In addition, there are also several lines of evidence implicating BRD4 in the splicing process. Using RNA-sequencing we found a significant increase in splicing inhibition, in particular intron retentions (IR), following heat treatment in BRD4-depleted cells. This leads to a decrease of mRNA abundancy of the affected transcripts, most likely due to premature termination codons. Subsequent experiments revealed that BRD4 interacts with the heat shock factor 1 (HSF1) such that under heat stress BRD4 is recruited to nuclear stress bodies and non-coding SatIII RNA transcripts are up-regulated. These findings implicate BRD4 as an important regulator of splicing during heat stress. Our data which links BRD4 to the stress induced splicing process may provide novel mechanisms of BRD4 inhibitors in regard to anti-cancer therapies.

ExSurv: A Web Resource for Prognostic Analyses of Exons Across Human Cancers Using Clinical Transcriptomes

Survival analysis in biomedical sciences is generally performed by correlating the levels of cellular components with patients' clinical features as a common practice in prognostic biomarker discovery. While the common and primary focus of such analysis in cancer genomics so far has been to identify the potential prognostic genes, alternative splicing - a posttranscriptional regulatory mechanism that affects the functional form of a protein due to inclusion or exclusion of individual exons giving rise to alternative protein products, has increasingly gained attention due to the prevalence of splicing aberrations in cancer transcriptomes. Hence, uncovering the potential prognostic exons can not only help in rationally designing exon-specific therapeutics but also increase specificity toward more personalized treatment options. To address this gap and to provide a platform for rational identification of prognostic exons from cancer transcriptomes, we developed ExSurv (https://exsurv.soic.iupui.edu), a web-based platform for predicting the survival contribution of all annotated exons in the human genome using RNA sequencing-based expression profiles for cancer samples from four cancer types available from The Cancer Genome Atlas. ExSurv enables users to search for a gene of interest and shows survival probabilities for all the exons associated with a gene and found to be significant at the chosen threshold. ExSurv also includes raw expression values across the cancer cohort as well as the survival plots for prognostic exons. Our analysis of the resulting prognostic exons across four cancer types revealed that most of the survival-associated exons are unique to a cancer type with few processes such as cell adhesion, carboxylic, fatty acid metabolism, and regulation of T-cell signaling common across cancer types, possibly suggesting significant differences in the posttranscriptional regulatory pathways contributing to prognosis.

Genome stability versus transcript diversity

Our genome is protected from the introduction of mutations by high fidelity replication and an extensive network of DNA damage response and repair mechanisms. However, the expression of our genome, via RNA and protein synthesis, allows for more diversity in translating genetic information. In addition, the splicing process has become less stringent over evolutionary time allowing for a substantial increase in the diversity of transcripts generated. The result is a diverse transcriptome and proteome that harbor selective advantages over a more tightly regulated system. Here, we describe mechanisms in place that both safeguard the genome and promote translational diversity, with emphasis on post-transcriptional RNA processing.

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.

Global intron retention mediated gene regulation during CD4+ T cell activation

T cell activation is a well-established model for studying cellular responses to exogenous stimulation. Using strand-specific RNA-seq, we observed that intron retention is prevalent in polyadenylated transcripts in resting CD4⁺ T cells and is significantly reduced upon T cell activation. Several lines of evidence suggest that intron-retained transcripts are less stable than fully spliced transcripts. Strikingly, the decrease in intron retention (IR) levels correlate with the increase in steady-state mRNA levels. Further, the majority of the genes upregulated in activated T cells are accompanied by a significant reduction in IR. Of these 1583 genes, 185 genes are predominantly regulated at the IR level, and highly enriched in the proteasome pathway, which is essential for proper T cell proliferation and cytokine release. These observations were corroborated in both human and mouse CD4⁺ T cells. Our study revealed a novel post-transcriptional regulatory mechanism that may potentially contribute to coordinated and/or quick cellular responses to extracellular stimuli such as an acute infection.

The alternative splicing program of differentiated smooth muscle cells involves concerted non-productive splicing of post-transcriptional regulators

Alternative splicing (AS) is a key component of gene expression programs that drive cellular differentiation. Smooth muscle cells (SMCs) are important in the function of a number of physiological systems; however, investigation of SMC AS has been restricted to a handful of events. We profiled transcriptome changes in mouse de-differentiating SMCs and observed changes in hundreds of AS events. Exons included in differentiated cells were characterized by particularly weak splice sites and by upstream binding sites for Polypyrimidine Tract Binding protein (PTBP1). Consistent with this, knockdown experiments showed that that PTBP1 represses many smooth muscle specific exons. We also observed coordinated splicing changes predicted to downregulate the expression of core components of U1 and U2 snRNPs, splicing regulators and other post-transcriptional factors in differentiated cells. The levels of cognate proteins were lower or similar in differentiated compared to undifferentiated cells. However, levels of snRNAs did not follow the expression of splicing proteins, and in the case of U1 snRNP we saw reciprocal changes in the levels of U1 snRNA and U1 snRNP proteins. Our results suggest that the AS program in differentiated SMCs is orchestrated by the combined influence of auxiliary RNA binding proteins, such as PTBP1, along with altered activity and stoichiometry of the core splicing machinery.

Functional genomics analyses of RNA-binding proteins reveal the splicing regulator SNRPB as an oncogenic candidate in glioblastoma

BACKGROUND

Glioblastoma (GBM) is the most common and aggressive type of brain tumor. Currently, GBM has an extremely poor outcome and there is no effective treatment. In this context, genomic and transcriptomic analyses have become important tools to identify new avenues for therapies. RNA-binding proteins (RBPs) are master regulators of co- and post-transcriptional events; however, their role in GBM remains poorly understood. To further our knowledge of novel regulatory pathways that could contribute to gliomagenesis, we have conducted a systematic study of RBPs in GBM.

RESULTS

By measuring expression levels of 1542 human RBPs in GBM samples and glioma stem cell samples, we identified 58 consistently upregulated RBPs. Survival analysis revealed that increased expression of 21 RBPs was also associated with a poor prognosis. To assess the functional impact of those RBPs, we modulated their expression in GBM cell lines and performed viability, proliferation, and apoptosis assays. Combined results revealed a prominent oncogenic candidate, SNRPB, which encodes core spliceosome machinery components. To reveal the impact of SNRPB on splicing and gene expression, we performed its knockdown in a GBM cell line followed by RNA sequencing. We found that the affected genes were involved in RNA processing, DNA repair, and chromatin remodeling. Additionally, genes and pathways already associated with gliomagenesis, as well as a set of general cancer genes, also presented with splicing and expression alterations.

CONCLUSIONS

Our study provides new insights into how RBPs, and specifically SNRPB, regulate gene expression and directly impact GBM development.

SURVIV for survival analysis of mRNA isoform variation

The rapid accumulation of clinical RNA-seq data sets has provided the opportunity to associate mRNA isoform variations to clinical outcomes. Here we report a statistical method SURVIV (Survival analysis of mRNA Isoform Variation), designed for identifying mRNA isoform variation associated with patient survival time. A unique feature and major strength of SURVIV is that it models the measurement uncertainty of mRNA isoform ratio in RNA-seq data. Simulation studies suggest that SURVIV outperforms the conventional Cox regression survival analysis, especially for data sets with modest sequencing depth. We applied SURVIV to TCGA RNA-seq data of invasive ductal carcinoma as well as five additional cancer types. Alternative splicing-based survival predictors consistently outperform gene expression-based survival predictors, and the integration of clinical, gene expression and alternative splicing profiles leads to the best survival prediction. We anticipate that SURVIV will have broad utilities for analysing diverse types of mRNA isoform variation in large-scale clinical RNA-seq projects.

Lessons from non-canonical splicing

Recent improvements in experimental and computational techniques that are used to study the transcriptome have enabled an unprecedented view of RNA processing, revealing many previously unknown non-canonical splicing events. This includes cryptic events located far from the currently annotated exons and unconventional splicing mechanisms that have important roles in regulating gene expression. These non-canonical splicing events are a major source of newly emerging transcripts during evolution, especially when they involve sequences derived from transposable elements. They are therefore under precise regulation and quality control, which minimizes their potential to disrupt gene expression. We explain how non-canonical splicing can lead to aberrant transcripts that cause many diseases, and also how it can be exploited for new therapeutic strategies.

Experimental approaches to studying the nature and impact of splicing variation in zebrafish

From a fixed number of genes carried in all cells, organisms create considerable diversity in cellular phenotype through differential regulation of gene expression. One prevalent source of transcriptome diversity is alternative pre-mRNA splicing, which is manifested in many different forms. Zebrafish models of splicing dysfunction due to mutated spliceosome components provide opportunity to link biochemical analyses of spliceosome structure and function with whole organism phenotypic outcomes. Drawing from experience with two zebrafish mutants: cephalophŏnus (a prpf8 mutant, isolated for defects in granulopoiesis) and caliban (a rnpc3 mutant, isolated for defects in digestive organ development), we describe the use of glycerol gradient sedimentation and native gel electrophoresis to resolve components of aberrant splicing complexes. We also describe how RNAseq can be employed to examine relatively rare alternative splicing events including intron retention. Such experimental approaches in zebrafish can promote understanding of how splicing variation and dysfunction contribute to phenotypic diversity and disease pathogenesis.

A dynamic intron retention program in the mammalian megakaryocyte and erythrocyte lineages

Intron retention (IR) is a form of alternative splicing that can impact mRNA levels through nonsense-mediated decay or by nuclear mRNA detention. A complex, dynamic IR pattern has been described in maturing mammalian granulocytes, but it is unknown whether IR occurs broadly in other hematopoietic lineages. We globally assessed IR in primary maturing mammalian erythroid and megakaryocyte (MK) lineages as well as their common progenitor cells (MEPs). Both lineages exhibit an extensive differential IR program involving hundreds of introns and genes with an overwhelming loss of IR in erythroid cells and MKs compared to MEPs. Moreover, complex IR patterns were seen throughout murine erythroid maturation. Similarly complex patterns were observed in human erythroid differentiation, but not involving the murine orthologous introns or genes. Despite the common origin of erythroid cells and MKs, and overlapping gene expression patterns, the MK IR program is entirely distinct from that of the erythroid lineage with regards to introns, genes, and affected gene ontologies. Importantly, our results suggest that IR serves to broadly regulate mRNA levels. These findings highlight the importance of this understudied form of alternative splicing in gene regulation and provide a useful resource for studies on gene expression in the MK and erythroid lineages.

Using high-throughput sequencing transcriptome data for INDEL detection: challenges for cancer drug discovery

INTRODUCTION

A cancer cell is a mosaic of genomic and epigenomic alterations. Distinct cancer molecular signatures can be observed depending on tumor type or patient genetic background. One type of genomic alteration is the insertion and/or deletion (INDEL) of nucleotides in the DNA sequence, which may vary in length, and may change the encoded protein or modify protein domains. INDELs are associated to a large number of diseases and their detection is done based on low-throughput techniques. However, high-throughput sequencing has also started to be used for detection of novel disease-causing INDELs. This search may identify novel drug targets.

AREAS COVERED

This review presents examples of using high-throughput sequencing (DNA-Seq and RNA-Seq) to investigate the incidence of INDELs in coding regions of human genes. Some of these examples successfully utilized RNA-Seq to identify INDELs associated to diseases. In addition, other studies have described small INDELs related to chemo-resistance or poor outcome of patients, while structural variants were associated with a better clinical outcome.

EXPERT OPINION

On average, there is twice as much RNA-Seq data available at the most used repositories for such data compared to DNA-Seq. Therefore, using RNA-Seq data is a promising strategy for studying cancer samples with unknown mechanisms of drug resistance, aiming at the discovery of proteins with potential as novel drug targets.

Methylome-wide Sequencing Detects DNA Hypermethylation Distinguishing Indolent from Aggressive Prostate Cancer

A critical need in understanding the biology of prostate cancer is characterizing the molecular differences between indolent and aggressive cases. Because DNA methylation can capture the regulatory state of tumors, we analyzed differential methylation patterns genome-wide among benign prostatic tissue and low-grade and high-grade prostate cancer and found extensive, focal hypermethylation regions unique to high-grade disease. These hypermethylation regions occurred not only in the promoters of genes but also in gene bodies and at intergenic regions that are enriched for DNA-protein binding sites. Integration with existing RNA-sequencing (RNA-seq) and survival data revealed regions where DNA methylation correlates with reduced gene expression associated with poor outcome. Regions specific to aggressive disease are proximal to genes with distinct functions from regions shared by indolent and aggressive disease. Our compendium of methylation changes reveals crucial molecular distinctions between indolent and aggressive prostate cancer.

A dynamic intron retention program enriched in RNA processing genes regulates gene expression during terminal erythropoiesis

Differentiating erythroblasts execute a dynamic alternative splicing program shown here to include extensive and diverse intron retention (IR) events. Cluster analysis revealed hundreds of developmentally-dynamic introns that exhibit increased IR in mature erythroblasts, and are enriched in functions related to RNA processing such as SF3B1 spliceosomal factor. Distinct, developmentally-stable IR clusters are enriched in metal-ion binding functions and include mitoferrin genes SLC25A37 and SLC25A28 that are critical for iron homeostasis. Some IR transcripts are abundant, e.g. comprising ∼50% of highly-expressed SLC25A37 and SF3B1 transcripts in late erythroblasts, and thereby limiting functional mRNA levels. IR transcripts tested were predominantly nuclear-localized. Splice site strength correlated with IR among stable but not dynamic intron clusters, indicating distinct regulation of dynamically-increased IR in late erythroblasts. Retained introns were preferentially associated with alternative exons with premature termination codons (PTCs). High IR was observed in disease-causing genes including SF3B1 and the RNA binding protein FUS. Comparative studies demonstrated that the intron retention program in erythroblasts shares features with other tissues but ultimately is unique to erythropoiesis. We conclude that IR is a multi-dimensional set of processes that post-transcriptionally regulate diverse gene groups during normal erythropoiesis, misregulation of which could be responsible for human disease.

The spliceosome, a potential Achilles heel of MYC-driven tumors

Alterations in RNA splicing are frequent in human tumors. Two recent studies of lymphoma and breast cancer have identified components of the spliceosome - the core splicing machinery - that are essential for malignant transformation driven by the transcription factor MYC. These findings provide a direct link between MYC and RNA splicing deregulation, and raise the exciting possibility of targeting spliceosome components in MYC-driven tumors.

Aberrant RNA splicing in cancer; expression changes and driver mutations of splicing factor genes

Alternative splicing is a widespread process contributing to structural transcript variation and proteome diversity. In cancer, the splicing process is commonly disrupted, resulting in both functional and non-functional end-products. Cancer-specific splicing events are known to contribute to disease progression; however, the dysregulated splicing patterns found on a genome-wide scale have until recently been less well-studied. In this review, we provide an overview of aberrant RNA splicing and its regulation in cancer. We then focus on the executors of the splicing process. Based on a comprehensive catalog of splicing factor encoding genes and analyses of available gene expression and somatic mutation data, we identify cancer-associated patterns of dysregulation. Splicing factor genes are shown to be significantly differentially expressed between cancer and corresponding normal samples, and to have reduced inter-individual expression variation in cancer. Furthermore, we identify enrichment of predicted cancer-critical genes among the splicing factors. In addition to previously described oncogenic splicing factor genes, we propose 24 novel cancer-critical splicing factors predicted from somatic mutations.