Reducing the structure bias of RNA-Seq reveals a large number of non-annotated non-coding RNA.

The study of RNA expression is the fastest growing area of genomic research. However, despite the dramatic increase in the number of sequenced transcriptomes, we still do not have accurate estimates of the number and expression levels of non-coding RNA genes. Non-coding transcripts are often overlooked due to incomplete genome annotation. In this study, we use annotation-independent detection of RNA reads generated using a reverse transcriptase with low structure bias to identify non-coding RNA. Transcripts between 20 and 500 nucleotides were filtered and crosschecked with non-coding RNA annotations revealing 111 non-annotated non-coding RNAs expressed in different cell lines and tissues. Inspecting the sequence and structural features of these transcripts indicated that 60% of these transcripts correspond to new snoRNA and tRNA-like genes. The identified genes exhibited features of their respective families in terms of structure, expression, conservation and response to depletion of interacting proteins. Together, our data reveal a new group of RNA that are difficult to detect using standard gene prediction and RNA sequencing techniques, suggesting that reliance on actual gene annotation and sequencing techniques distorts the perceived architecture of the human transcriptome.

snoDB: an interactive database of human snoRNA sequences, abundance and interactions.

Small nucleolar RNAs (snoRNAs) are an abundant type of non-coding RNA with conserved functions in all known eukaryotes. Classified into two main families, the box C/D and H/ACA snoRNAs, they enact their most well characterized role of guiding site specific modifications in ribosomal RNA, through the formation of specific ribonucleoprotein complexes, with fundamental implications in ribosome biogenesis. However, it is becoming increasingly clear that the landscape of snoRNA cellular functionality is much broader than it once seemed with novel members, non-uniform expression patterns, new and diverse targets as well as several emerging non-canonical functions ranging from the modulation of alternative splicing to the regulation of chromatin architecture. In order to facilitate the further characterization of human snoRNAs in a holistic manner, we introduce an online interactive database tool: snoDB. Its purpose is to consolidate information on human snoRNAs from different sources such as sequence databases, target information, both canonical and non-canonical from the literature and from high-throughput RNA-RNA interaction datasets, as well as high-throughput sequencing data that can be visualized interactively.

Introns regulate the production of ribosomal proteins by modulating splicing of duplicated ribosomal protein genes.

Most budding yeast introns exist in the many duplicated ribosomal protein genes (RPGs) and it has been posited that they remain there to modulate the expression of RPGs and cell growth in response to stress. However, the mechanism by which introns regulate the expression of RPGs and their impact on the synthesis of ribosomal proteins remain unclear. In this study, we show that introns determine the ratio of ribosomal protein isoforms through asymmetric paralog-specific regulation of splicing. Exchanging the introns and 3' untranslated regions of the duplicated RPS9 genes altered the splicing efficiency and changed the ratio of the ribosomal protein isoforms. Mutational analysis of the RPS9 genes indicated that splicing is regulated by variations in the intron structure and the 3' untranslated region. Together these data suggest that preferential splicing of duplicated RPGs provides a means for adjusting the ratio of different ribosomal protein isoforms, while maintaining the overall expression level of each ribosomal protein.

Transcriptome wide annotation of eukaryotic RNase III reactivity and degradation signals.

Detection and validation of the RNA degradation signals controlling transcriptome stability are essential steps for understanding how cells regulate gene expression. Here we present complete genomic and biochemical annotations of the signals required for RNA degradation by the dsRNA specific ribonuclease III (Rnt1p) and examine its impact on transcriptome expression. Rnt1p cleavage signals are randomly distributed in the yeast genome, and encompass a wide variety of sequences, indicating that transcriptome stability is not determined by the recurrence of a fixed cleavage motif. Instead, RNA reactivity is defined by the sequence and structural context in which the cleavage sites are located. Reactive signals are often associated with transiently expressed genes, and their impact on RNA expression is linked to growth conditions. Together, the data suggest that Rnt1p reactivity is triggered by malleable RNA degradation signals that permit dynamic response to changes in growth conditions.

Cancer-Associated Perturbations in Alternative Pre-messenger RNA Splicing.

For most of our 25,000 genes, the removal of introns by pre-messenger RNA (pre-mRNA) splicing represents an essential step toward the production of functional messenger RNAs (mRNAs). Alternative splicing of a single pre-mRNA results in the production of different mRNAs. Although complex organisms use alternative splicing to expand protein function and phenotypic diversity, patterns of alternative splicing are often altered in cancer cells. Alternative splicing contributes to tumorigenesis by producing splice isoforms that can stimulate cell proliferation and cell migration or induce resistance to apoptosis and anticancer agents. Cancer-specific changes in splicing profiles can occur through mutations that are affecting splice sites and splicing control elements, and also by alterations in the expression of proteins that control splicing decisions. Recent progress in global approaches that interrogate splicing diversity should help to obtain specific splicing signatures for cancer types. The development of innovative approaches for annotating and reprogramming splicing events will more fully establish the essential contribution of alternative splicing to the biology of cancer and will hopefully provide novel targets and anticancer strategies. Metazoan genes are usually made up of several exons interrupted by introns. The introns are removed from the pre-mRNA by RNA splicing. In conjunction with other maturation steps, such as capping and polyadenylation, the spliced mRNA is then transported to the cytoplasm to be translated into a functional protein. The basic mechanism of splicing requires accurate recognition of each extremity of each intron by the spliceosome. Introns are identified by the binding of U1 snRNP to the 5' splice site and the U2AF65/U2AF35 complex to the 3' splice site. Following these interactions, other proteins and snRNPs are recruited to generate the complete spliceosomal complex needed to excise the intron. While many introns are constitutively removed by the spliceosome, other splice junctions are not used systematically, generating the phenomenon of alternative splicing. Alternative splicing is therefore the process by which a single species of pre-mRNA can be matured to produce different mRNA molecules (Fig. 1). Depending on the number and types of alternative splicing events, a pre-mRNA can generate from two to several thousands different mRNAs leading to the production of a corresponding number of proteins. It is now believed that the expression of at least 70 % of human genes is subjected to alternative splicing, implying an enormous contribution to proteomic diversity, and by extension, to the development and the evolution of complex animals. Defects in splicing have been associated with human diseases (Caceres and Kornblihtt, Trends Genet 18(4):186-93, 2002, Cartegni et al., Nat Rev Genet 3(4):285-98, 2002, Pagani and Baralle, Nat Rev Genet 5(5):389-96, 2004), including cancer (Brinkman, Clin Biochem 37(7):584-94, 2004, Venables, Bioessays 28(4):378-86, 2006, Srebrow and Kornblihtt, J Cell Sci 119(Pt 13):2635-2641, 2006, Revil et al., Bull Cancer 93(9):909-919, 2006, Venables, Transworld Res Network, 2006, Pajares et al., Lancet Oncol 8(4):349-57, 2007, Skotheim and Nees, Int J Biochem Cell Biol 39:1432-1449, 2007). Numerous studies have now confirmed the existence of specific differences in the alternative splicing profiles between normal and cancer tissues. Although there are a few cases where specific mutations are the primary cause for these changes, global alterations in alternative splicing in cancer cells may be primarily derived from changes in the expression of RNA-binding proteins that control splice site selection. Overall, these cancer-specific differences in alternative splicing offer an immense potential to improve the diagnosis and the prognosis of cancer. This review will focus on the functional impact of cancer-associated alternative splicing variants, the molecular determinants that alter the splicing decisions in cancer cells, and future therapeutic strategies.

Proteins associated with the exon junction complex also control the alternative splicing of apoptotic regulators.

Several apoptotic regulators, including Bcl-x, are alternatively spliced to produce isoforms with opposite functions. We have used an RNA interference strategy to map the regulatory landscape controlling the expression of the Bcl-x splice variants in human cells. Depleting proteins known as core (Y14 and eIF4A3) or auxiliary (RNPS1, Acinus, and SAP18) components of the exon junction complex (EJC) improved the production of the proapoptotic Bcl-x(S) splice variant. This effect was not seen when we depleted EJC proteins that typically participate in mRNA export (UAP56, Aly/Ref, and TAP) or that associate with the EJC to enforce nonsense-mediated RNA decay (MNL51, Upf1, Upf2, and Upf3b). Core and auxiliary EJC components modulated Bcl-x splicing through different cis-acting elements, further suggesting that this activity is distinct from the established EJC function. In support of a direct role in splicing control, recombinant eIF4A3, Y14, and Magoh proteins associated preferentially with the endogenous Bcl-x pre-mRNA, interacted with a model Bcl-x pre-mRNA in early splicing complexes, and specifically shifted Bcl-x alternative splicing in nuclear extracts. Finally, the depletion of Y14, eIF4A3, RNPS1, SAP18, and Acinus also encouraged the production of other proapoptotic splice variants, suggesting that EJC-associated components are important regulators of apoptosis acting at the alternative splicing level.

Detection and characterization of transcription termination.

In most eukaryotes, the generation of the 3' end and transcription termination are initiated by cleavage of the pre-mRNA upstream of the polyadenylation site. This cleavage initiates 5'-3' degradation of the 3' end cleavage product by the exoribonuclease Rat1p leading to the dissociation of the RNA polymerase II (RNAPII) complex. The Rat1p-dependent transcription termination was also shown to be initiated by a polyadenylation-independent cleavage performed by the double-stranded RNA-specific ribonuclease (RNase) III (Rnt1p) suggesting that the majority of transcription termination events are RNase dependent. Therefore, it became essential for future studies on transcription termination to carefully consider both the nature of the RNase-dependent RNA transcripts and the association pattern of the RNAPII with the transcriptional unit. Here, we present methods allowing the evaluation of the impact of yeast RNases on the 3' end formation and their contribution to transcription termination. Northern blot analysis of transcripts generated downstream of known genes in the absence of RNases identifies potential transcription termination sites while chromatin immunoprecipitation of RNAPII differentiates between termination- and transcription-independent processing events.

Introns within ribosomal protein genes regulate the production and function of yeast ribosomes.

In budding yeast, the most abundantly spliced pre-mRNAs encode ribosomal proteins (RPs). To investigate the contribution of splicing to ribosome production and function, we systematically eliminated introns from all RP genes to evaluate their impact on RNA expression, pre-rRNA processing, cell growth, and response to stress. The majority of introns were required for optimal cell fitness or growth under stress. Most introns are found in duplicated RP genes, and surprisingly, in the majority of cases, deleting the intron from one gene copy affected the expression of the other in a nonreciprocal manner. Consistently, 70% of all duplicated genes were asymmetrically expressed, and both introns and gene deletions displayed copy-specific phenotypic effects. Together, our results indicate that splicing in yeast RP genes mediates intergene regulation and implicate the expression ratio of duplicated RP genes in modulating ribosome function.

Alternative splicing of SYK regulates mitosis and cell survival.

Most human genes produce multiple mRNA isoforms through alternative splicing. However, the biological relevance of most splice variants remains unclear. In this study, we evaluated the functional impact of alternative splicing in cancer cells. We modulated the splicing pattern of 41 cancer-associated splicing events and scored the effects on cell growth, viability and apoptosis, identifying three isoforms essential for cell survival. Specifically, changing the splicing pattern of the spleen tyrosine kinase gene (SYK) impaired cell-cycle progression and anchorage-independent growth. Notably, exposure of cancer cells to epithelial growth factor modulated the SYK splicing pattern to promote the pro-survival isoform that is associated with cancer tissues in vivo. The data suggest that splicing of selected genes is specifically modified during tumor development to allow the expression of isoforms that promote cancer cell survival.

High-throughput quantification of splicing isoforms.

Most human messenger RNAs (mRNAs) are alternatively spliced and many exhibit disease-specific splicing patterns. However, the contribution of most splicing events to the development and maintenance of human diseases remains unclear. As the contribution of alternative splicing events to diagnosis and prognosis is becoming increasingly recognized, it becomes important to develop precise methods to quantify the abundance of these isoforms in clinical samples. Here we present a pipeline for real-time PCR annotation of splicing events (RASE) that allows accurate identification of a large number of splicing isoforms in human tissues. The RASE automatically designed specific primer pairs for 81% of all alternative splicing events in the NCBI build 36 database. Experimentally, the majority of the RASE designed primers resulted in isoform-specific amplification suitable for quantification in human cell lines or in formalin-fixed, paraffin-embedded (FFPE) RNA extract. Using this pipeline it is now possible to rapidly identify splicing isoform signatures in different types of human tissues or to validate complete sets of data generated by microarray expression profiling and deep sequencing techniques.

Yeast RNase III triggers polyadenylation-independent transcription termination.

Transcription termination of messenger RNA (mRNA) is normally achieved by polyadenylation followed by Rat1p-dependent 5'-3' exoribonuleolytic degradation of the downstream transcript. Here we show that the yeast ortholog of the dsRNA-specific ribonuclease III (Rnt1p) may trigger Rat1p-dependent termination of RNA transcripts that fail to terminate near polyadenylation signals. Rnt1p cleavage sites were found downstream of several genes, and the deletion of RNT1 resulted in transcription readthrough. Inactivation of Rat1p impaired Rnt1p-dependent termination and resulted in the accumulation of 3' end cleavage products. These results support a model for transcription termination in which cotranscriptional cleavage by Rnt1p provides access for exoribonucleases in the absence of polyadenylation signals.

Control of alternative splicing through siRNA-mediated transcriptional gene silencing.

When targeting promoter regions, small interfering RNAs (siRNAs) trigger a previously proposed pathway known as transcriptional gene silencing by promoting heterochromatin formation. Here we show that siRNAs targeting intronic or exonic sequences close to an alternative exon regulate the splicing of that exon. The effect occurred in hepatoma and HeLa cells with siRNA antisense strands designed to enter the silencing pathway, suggesting hybridization with nascent pre-mRNA. Unexpectedly, in HeLa cells the sense strands were also effective, suggesting that an endogenous antisense transcript, detectable in HeLa but not in hepatoma cells, acts as a target. The effect depends on Argonaute-1 and is counterbalanced by factors favoring chromatin opening or transcriptional elongation. The increase in heterochromatin marks (dimethylation at Lys9 and trimethylation at Lys27 of histone H3) at the target site, the need for the heterochromatin-associated protein HP1alpha and the reduction in RNA polymerase II processivity suggest a mechanism involving the kinetic coupling of transcription and alternative splicing.

Cancer-associated regulation of alternative splicing.

Alternative splicing of pre-mRNA increases the diversity of protein functions. Here we show that about half of all active alternative splicing events in ovarian and breast tissues are changed in tumors, and many seem to be regulated by a single factor; sequence analysis revealed binding sites for the RNA binding protein FOX2 downstream of one-third of the exons skipped in cancer. High-resolution analysis of FOX2 binding sites defined the precise positions relative to alternative exons at which the protein may function as either a silencer or an enhancer. Most of the identified targets were shifted in the same direction by FOX2 depletion in cell lines as they were in breast and ovarian cancer tissues. Notably, we found expression of FOX2 itself is downregulated in ovarian cancer and its splicing is altered in breast cancer samples. These results suggest that the decreased expression of FOX2 in cancer tissues modulates splicing and controls proliferation.

Identification of alternative splicing markers for breast cancer.

Breast cancer is the most common cause of cancer death among women under age 50 years, so it is imperative to identify molecular markers to improve diagnosis and prognosis of this disease. Here, we present a new approach for the identification of breast cancer markers that does not measure gene expression but instead uses the ratio of alternatively spliced mRNAs as its indicator. Using a high-throughput reverse transcription-PCR-based system for splicing annotation, we monitored the alternative splicing profiles of 600 cancer-associated genes in a panel of 21 normal and 26 cancerous breast tissues. We validated 41 alternative splicing events that significantly differed in breast tumors relative to normal breast tissues. Most cancer-specific changes in splicing that disrupt known protein domains support an increase in cell proliferation or survival consistent with a functional role for alternative splicing in cancer. In a blind screen, a classifier based on the 12 best cancer-associated splicing events correctly identified cancer tissues with 96% accuracy. Moreover, a subset of these alternative splicing events could order tissues according to histopathologic grade, and 5 markers were validated in a further blind set of 19 grade 1 and 19 grade 3 tumor samples. These results provide a simple alternative for the classification of normal and cancerous breast tumor tissues and underscore the putative role of alternative splicing in the biology of cancer.

Multiple and specific mRNA processing targets for the major human hnRNP proteins.

Alternative splicing is a key mechanism regulating gene expression, and it is often used to produce antagonistic activities particularly in apoptotic genes. Heterogeneous nuclear ribonucleoparticle (hnRNP) proteins form a family of RNA-binding proteins that coat nascent pre-mRNAs. Many but not all major hnRNP proteins have been shown to participate in splicing control. The range and specificity of hnRNP protein action remain poorly documented, even for those affecting splice site selection. We used RNA interference and a reverse transcription-PCR screening platform to examine the implications of 14 of the major hnRNP proteins in the splicing of 56 alternative splicing events in apoptotic genes. Out of this total of 784 alternative splicing reactions tested in three human cell lines, 31 responded similarly to a knockdown in at least two different cell lines. On the other hand, the impact of other hnRNP knockdowns was cell line specific. The broadest effects were obtained with hnRNP K and C, two proteins whose role in alternative splicing had not previously been firmly established. Different hnRNP proteins affected distinct sets of targets with little overlap even between closely related hnRNP proteins. Overall, our study highlights the potential contribution of all of these major hnRNP proteins in alternative splicing control and shows that the targets for individual hnRNP proteins can vary in different cellular contexts.

Anticancer drugs affect the alternative splicing of Bcl-x and other human apoptotic genes.

Inducing an apoptotic response is the goal of most current chemotherapeutic interventions against cancer. However, little is known about the effect of chemotherapeutic agents on the alternative splicing of apoptotic genes. Here, we have tested 20 of the mainstream anticancer drugs for their ability to influence the production of Bcl-x splice isoforms. We find that many drugs shift splicing toward the proapoptotic Bcl-x(S) splice variant in 293 cells. The drugs modulate splicing decisions most likely through signaling events because the splicing switch is not compromised by inhibiting de novo protein synthesis or the activity of caspases. Several drugs also shift Bcl-x splicing in cancer cell lines (MCF-7, HeLa, PC-3, PA-1, and SKOV-3), but the set of active drugs varies between cell lines. We also examined the effect of anticancer agents on the alternative splicing of 95 other human apoptotic genes in different cell lines. Almost every drug can alter a subset of alternative splicing events in each cell line. Although drugs of the same class often influence the alternative splicing of the same units in individual cell lines, these units differ considerably between cell lines, indicating cell line-specific differences in the pathways that control splicing.

A G-tract element in apoptotic agents-induced alternative splicing.

Alternative splicing of a single pre-mRNA transcript can produce protein isoforms that promote either cell growth or death. Here we show that Ro-31-8220 (Ro), an apoptotic agent that inhibits protein kinase C and activates the c-Jun N terminal kinase, decreased the proportion of the cell growth-promoting Bcl-xL splice variant. Targeted mutagenesis analyses narrowed down a critical sequence to a 16-nt G-tract element (Gt16). Transferring this element to a heterologous gene conferred Ro response on an otherwise constitutive exon. The Ro effect was reduced by okadaic acid, an inhibitor of protein phosphatases PP1 and PP2A, in a concentration-dependent manner. Search in the human genome followed by RT-PCR identified a group of genes that contain similar exonic G-tract elements and are responsive to Ro. Moreover, the Gt16 element also mediates the regulation of alternative splicing by other cell apoptosis-inducers particularly retinoic acid. Therefore, the G-tract element likely plays a role in the apoptotic agents-induced alternative splicing of a group of genes. The functions of these genes imply that this regulation will have impact on cell growth/death.

Modern origin of numerous alternatively spliced human introns from tandem arrays.

Despite the widespread occurrence of spliceosomal introns in the genomes of higher eukaryotes, their origin remains controversial. One model proposes that the duplication of small genomic portions could have provided the boundaries for new introns. If this mechanism has occurred recently, the 5' and 3' boundaries of each resulting intron should display distinctive sequence similarity. Here, we report that the human genome contains an excess of introns with perfect matching sequences at boundaries. One-third of these introns interrupt the protein-coding sequences of known genes. Introns with the best-matching boundaries are invariably found in tandem arrays of direct repeats. Sequence analysis of the arrays indicates that many intron-breeding repeats have disseminated in several genes at different times during human evolution. A comparison with orthologous regions in mouse and chimpanzee suggests a young age for the human introns with the most-similar boundaries. Finally, we show that these human introns are alternatively spliced with exceptionally high frequency. Our study indicates that genomic duplication has been an important mode of intron gain in mammals. The alternative splicing of transcripts containing these intron-breeding repeats may provide the plasticity required for the rapid evolution of new human proteins.

RNase III-dependent regulation of yeast telomerase.

In bakers' yeast, in vivo telomerase activity requires a ribonucleoprotein (RNP) complex with at least four associated proteins (Est2p, Est1p, Est3p, and Cdc13p) and one RNA species (Tlc1). The function of telomerase in maintaining chromosome ends, called telomeres, is tightly regulated and linked to the cell cycle. However, the mechanisms that regulate the expression of individual components of telomerase are poorly understood. Here we report that yeast RNase III (Rnt1p), a double-stranded RNA-specific endoribonuclease, regulates the expression of telomerase subunits and is required for maintaining normal telomere length. Deletion or inactivation of RNT1 induced the expression of Est1, Est2, Est3, and Tlc1 RNAs and increased telomerase activity, leading to elongation of telomeric repeat tracts. In silico analysis of the different RNAs coding for the telomerase subunits revealed a canonical Rnt1p cleavage site near the 3' end of Est1 mRNA. This predicted structure was cleaved by Rnt1p and its disruption abolished cleavage in vitro. Mutation of the Rnt1p cleavage signal in vivo impaired the cell cycle-dependent degradation of Est1 mRNA without affecting its steady-state level. These results reveal a new mechanism that influences telomeres length by controlling the expression of the telomerase subunits.