IRFinder: assessing the impact of intron retention on mammalian gene expression
- Robert Middleton†1,
- Dadi Gao†1, 2, 3, 4,
- Aubin Thomas5,
- Babita Singh6,
- Amy Au4, 7,
- Justin J-L Wong4, 7, 8,
- Alexandra Bomane9,
- Bertrand Cosson9,
- Eduardo Eyras6, 10,
- John E. J. Rasko4, 7, 11 and
- William Ritchie1, 5, 12Email authorView ORCID ID profile
© The Author(s). 2017
Received: 8 December 2016
Accepted: 27 February 2017
Published: 15 March 2017
Intron retention (IR) occurs when an intron is transcribed into pre-mRNA and remains in the final mRNA. We have developed a program and database called IRFinder to accurately detect IR from mRNA sequencing data. Analysis of 2573 samples showed that IR occurs in all tissues analyzed, affects over 80% of all coding genes and is associated with cell differentiation and the cell cycle. Frequently retained introns are enriched for specific RNA binding protein sites and are often retained in clusters in the same gene. IR is associated with lower protein levels and intron-retaining transcripts that escape nonsense-mediated decay are not actively translated.
Alternative splicing (AS) affects up to 95% of multi-exonic genes in humans . The three main types of AS are exon skipping, alternative 5′ or 3′ usage and intron retention (IR). IR occurs when an intron is transcribed into pre-mRNA and remains in the final mRNA. It constitutes a class of AS that is often neglected because these events are difficult to measure reliably. IR can introduce functional elements within mRNAs  or alternatively may lead to the introduction of premature termination codons, resulting in degradation of the mRNA by a surveillance mechanism called nonsense-mediated decay (NMD) . This process of IR followed by NMD can downregulate up to 35% of alternatively spliced transcripts in mammals . The NMD pathway is absolutely essential for post-implantation embryonic development as shown by Upf1 nullizygosity . However, an obvious consequence of NMD is that mRNAs containing introns with premature termination codons are degraded and therefore difficult to quantify. Consequently, the role of IR in specific eukaryotic biological pathways has been poorly defined prior to the availability of ultra-deep sequencing technologies.
We recently discovered that IR combined with the NMD pathway is not a by-product of faulty splicing but rather a major driver of the cellular differentiation of granulocytes . In this pioneering study we developed an approach to correctly identify introns that are differentially retained during granulocytic differentiation from sequencing data. Our approach was subsequently used to uncover the role of IR in stem cell reprogramming where it regulated demethylation genes at specific stages of reprogramming . Since then, other publications have highlighted the importance of IR in gene regulation , differentiation , and cancer . Despite increasing evidence that IR can regulate hundreds of genes in numerous systems, current studies still fail to identify IR events in their transcriptomic data.
Here we have developed a significantly enhanced program in terms of sensitivity and speed for detecting retained introns and filtering samples that are inappropriate for IR analysis in terms of library preparation and quality. This novel program was validated using quantitative RT-PCR (RT-qPCR) against retained introns and a NMD knockdown experiment. IRFinder correctly identified IR events and measured the ratio of retained introns to correctly spliced introns with great accuracy. Using IRFinder, we analyzed 3435 human samples, of which 2573 were suitable for analysis. We found that IR occurs in hundreds of genes in all tissues analyzed, affects over 80% of all coding genes, and is associated with cell differentiation and cell cycle. Retained introns in the same genes are frequently adjacent to each other and intron-retaining genes cluster closely on the genome, suggesting a global mechanism that regulates multiple introns simultaneously. By analyzing mass spectrometry and ribosome sequencing data we discovered that intron-retaining genes have lower protein output and that IR transcripts that escape NMD are not actively translated. Finally, by comparing introns that are frequently retained amongst the 2573 samples with introns that are rarely retained, we discovered a distinct primary sequence signature amongst frequently retained introns that is enriched in RNA binding protein sites. These proteins modulate the level of IR and thus the level of repression of frequently retained introns. Our program to analyze IR from sequencing data is available at GitHub (https://github.com/williamritchie/IRFinder) and a database of IR in over 2000 human samples is freely available at IRBase (http://mimirna.centenary.org.au/irfinder/database/).
Fair measurement of intronic expression is challenged by numerous factors. Within introns, highly expressed features such as snoRNAs, microRNAs, or unannotated exons may erroneously inflate count-based measures of intronic expression. Conversely, low complexity regions, common in introns, prevent unique mapping of reads. Because retained introns are generally expressed at a fraction of their flanking exons, uncorrected biases can massively disrupt IR estimation.
IRFinder implements an end-to-end analysis of retained introns from mRNA sequencing data in multiple species. It includes alignment via the STAR algorithm, quality controls on the sample analyzed, IR detection, and quantification and statistics for comparing IR levels between multiple samples. We provide standalone scripts for each of these steps so they can be used independently and provide command line tools to chain them together for complete analysis. On pre-aligned sequencing data, our program can run on a desktop computer with less than 2 GB of memory and takes approximately 10 minutes to detect IR events. Because we use STAR to align reads, our end-to-end analysis with raw reads requires at least 48 GB of memory and depends mainly on STAR runtime.
Our tool facilitates the analysis of large amounts of online data by automatically testing samples for their suitability for IR detection. Unsuitable samples either have high levels of DNA contamination or have been mislabeled as mRNA sequencing when in fact they are other types of experiments such as genome sequencing or ChIP-seq. Although this information should be available from online repositories, we found that only 68% (2573/3774) were suitable for analysis in this study. This was mainly due to the incorrect use of the term “mRNA-Seq”, which was frequently used in whole RNA experiments and CHIP-Seq experiments. Our approach also uses a series of programmatically fast steps written in C++ that automatically detect and trim adapters from sequencing reads. Of the 3435 initial samples we analyzed, 3096 (91%) still had adapters in the sequencing reads despite having been already processed by an adaptor-trimming algorithm.
IRFinder was capable of estimating IR events with low coverage or low mappability as confirmed by RT-qPCR (Additional file 1: Figure S1a–e). When compared with the currently available tools MISO and DEXseq, our IRFinder had higher accuracy and precision (Additional file 2: Text, Figure S2, Tables S2 and S3).
The IRFinder algorithm and instruction manual are available at GitHub (https://github.com/williamritchie/IRFinder).
IRFinder detects IR in over 80% of coding genes and in all samples tested
To determine how consistent IR levels were within the same tissue type we calculated the average IR ratio of retained introns in eight distinct tissue types for which we had over three replicates. We found that there was a significant overlap in genes for which we detected retained introns between each replicate of the same tissue (p < 2.10E-6, hypergeometric test) and found that the average IR ratio was homogenous within each tissue (Fig. 1b).
Having observed that numerous genes had multiple retained introns and that these were often adjacent to each other (separated by only one exon), we calculated the distance between genes with retained introns (Fig. 1c) and how often retained introns within the same gene were adjacent to each other (Fig. 1d). We found that adjacent introns within the same gene were more frequently retained than expected and that genes with retained introns were significantly closer to each other than the same number of randomly selected expressed genes (Mann–Whitney test p = 3.10-3; Additional file 3). Although in this study we were able to confirm that retained introns were generally surrounded by weaker splice sites  (Additional file 4: Text and Figure S3), the above results on IR clusters indicate that IR is not solely regulated at each splicing junction but is also regulated by a more global means of regulation that encompasses multiple introns and even multiple genes. Multiple retained introns within the same gene could be regulated by transcription rate through the gene in agreement with the recent finding that individual retained introns are associated with an accumulation of the elongating form of RNA polymerase II .
We created an online database of IR calculated in these samples, which is available at IRBase (http://mimirna.centenary.org.au/irfinder/database).
IR is associated with reduced protein output
To determine whether retained introns were translated, we analyzed ribosome profiling data based on deep sequencing of ribosome-protected mRNA fragments using cycloheximide to stabilize all translating ribosomes (CHX Ribo-seq). These data were coupled with mRNA sequencing and quantitative translation initiation sequencing (QTI-seq), which uses lactimidomycin to preserve initiating ribosomes and puromycin to deplete elongating ribosomes . These data allowed us to simultaneously assess translation from all translating ribosomes (CHX Ribo-seq), or specifically from initiating ribosomes (QTI-seq), and IR in HEK293 cells (Fig. 2b; Additional file 5). We found 429 IR events (IR ratio >10%) in this dataset. None of these displayed a signal from the CHX Ribo-seq or QTI-seq above background levels. This suggests that even though intron-retaining transcripts are polyadenylated, the translation machinery does generally not translate the IR elements either because the ribosome dissociates before it reaches the IR sequence or because intron-retaining transcripts are made unavailable.
IR is associated with cell differentiation and the cell cycle
hiPSCs can be differentiated into highly homogenous neurons within 4 days by overexpressing the transcription factors Neurogenin-1 and Neurogenin-2 . This allowed the authors to define networks of expression at different stages of neural differentiation. We reanalyzed their data using IRFinder and discovered 120 alternatively retained introns between days 1 and 4 of neurogenin expression (Fig. 3a; Additional file 6: Table S4). These IR events clearly marked the different stages of neuron differentiation. Similar to other studies performed in various tissues [6, 17], we found that IR in the late stage of differentiation (day 4) was enriched in splicing factor genes (p = 2.5 × 10−5). We also found numerous genes involved in neurogenesis with high levels of IR amongst which were the neural transcription factors SOX15 and POU4F1. SOX15 is known to be highly expressed in undifferentiated cells and repressed upon neuron differentiation . In agreement with this, we found that mRNA levels of SOX15 dropped dramatically from 9.4 to 1.3 reads per kilobase per million mappable reads (RPKM) during neurogenin-induced differentiation (Fig. 3b). Interestingly, IR levels increased dramatically from 0 to 36% at these stages, indicating that a large proportion of transcripts retained introns. These retained introns had in-frame stop codons and were thus susceptible to degradation via the NMD pathway. Another gene, POU4F1, is activated following the neurogenic phase to regulate gene expression during neuron sensory differentiation . In agreement, we found that when POU4F1 expression is first induced the intron is mostly retained, but as gene expression increases IR decreases. Here also, the retained intron had an in-frame stop codon and IR coupled with NMD could account for the dramatic change in expression of POU4F1.
We then analyzed another set of data where the authors used genetically matched hESCs and hiPSCs derived from fibroblasts using the Sendai virus reprogramming method to prove that these cells are equivalent at a gene transcription level . To assess whether IR was also equivalent between these cells, we applied IRFinder to the sequencing data from this study (Fig. 3c; Additional file 6: Table S5). We found distinct patterns of IR for hESCs, hiPSCs, and fibroblasts. In agreement with this previous study, we found that fibroblasts derived from hESCs and hiPSCs had similar IR patterns. Moreover, we were able to find Gene Ontology (GO) categories significantly associated with the differential IR (945 introns; Fig. 3d). The most significant of these were epidermal growth factor (EGF), Laminin G, SH3 homology, and plasma membrane. EGF is essential for fibroblast growth and differentiation and has been described as such for over four decades . SH3 homology domains interact with long distance proteins that regulate the cytoskeleton . Interestingly, we found 47 genes with differential IR between iPSC and hESC fibroblasts that had not been detected in the initial study (Additional file 6: Table S6).
Finally, we investigated whether IR changed dynamically during the cell cycle in CD4+ T cells. IR has recently been shown to be an integral regulator of T-cell activation  and a recent study quantified the effects of the cell cycle on a population of naive CD4+ T-helper cells that were induced to differentiate toward a TH2 subtype . They used an initial set of 892 cell cycle genes to determine cell cycle heterogeneity. We reanalyzed the sequencing data using IRFinder and discovered 969 differentially retained introns with distinct patterns of retention for each stage of the cell cycle (Fig. 3e). These introns were retained from genes enriched for phosphoproteins and the cell cycle (p = 2E-17 and p = 8E-6 after Benjamini–Hochberg correction, respectively (Fig. 3f). Surprisingly, only 97 out of the 969 introns that we found to be alternatively retained belonged to genes used to determine cell cycle signatures in the initial study. This indicates that 874 retained introns could be used as new markers of cell cycle stages in sequencing data. When we measured exonic expression of genes with retained introns, the differential expression patterns were lost (Additional file 6: Figure S4), indicating that differential IR did not mirror global mRNA expression and thus could be used as a complementary analysis to discover cell cycle markers. Interestingly, genes from the mRNA splicing GO category were also found to be enriched, leading to the hypothesis of differential regulation of splicing factors and sequence-specific regulation of IR events.
Frequently retained introns are enriched for a subset of RNA binding sites
In this study we developed bioinformatics tools to study the impact of intron retention on gene regulation. Because of their length and low complexity, intron expression is difficult to measure and is affected by poorly annotated experiments. Our validated approach allowed us to measure IR in over 2573 samples to gain unique insight into the widespread nature of IR and how it affects gene regulation. We found that IR was prevalent in all samples we analyzed and that it affected over 80% of all protein coding genes making, IR a major regulator of gene expression. IR events cluster together within the same transcript and genes with IR events are closer than expected, indicating a global mechanism that regulates multiple IR events. This is in agreement with previous studies linking IR to transcription speed .
Given the impact of IR coupled with NMD on protein output , we compared protein and mRNA levels in nine tissue types and found that IR genes were significantly associated with lower protein levels. We also found no convincing evidence of translation from retained introns that escaped NMD and that were thus measurable in our samples. These transcripts may have escaped NMD because of the inefficient recruitment of UPF proteins (1, 2, or 3) to the terminating ribosome or inefficient degradation after the release of the ribosome or because they were not efficiently exported from the nucleus [17, 28, 29]. The commonly accepted role of NMD as a surveillance mechanism is supported by evidence that it prevents deleterious proteins from being created from mis-spliced transcripts in disease . However, in this study and our previous study in granulopoiesis we found that none of the IR events detected created any protein products. In our model, NMD coupled with IR is a regulator of gene expression. This becomes even more apparent given that the SR proteins SRSF1 and SRSF7, which we discovered here to regulate IR levels, are also involved in NMD . Not only do they enhance NMD activity but they are also associated with the exon junction complex core factors, essential for recognition by the NMD pathway . In this context, the NMD machinery and SR proteins would become essential elements for a widespread mechanism of gene regulation via IR by modulating IR levels but also the efficiency of NMD that degrades IR transcripts.
We found multiple examples where IR of functionally related genes was specifically timed during differentiation and the cell cycle. Specific subsets of genes were subject to IR at different time points. Analyses of deep sequencing data in time series give a snapshot of the transcriptome. The analysis of IR may, however, be more informative than just the observation of individual transcript levels at a given time. In the data we analyzed and in our previous studies , we found that IR transcripts were often actively transcribed but with lower than expected protein output due to NMD. Thus, in many cases IR events often highlight genes that were required in a previous stage of differentiation or the cell cycle but are being actively downregulated. Accordingly, they give insight into the previous and future states of transcripts’ fates.
Novel IRFinder algorithm
The IRFinder algorithm is an improved and extended version of our first algorithm  and includes now tools for preparing the genome, cleaning data, and testing the suitability of a given sample for IR analysis. Many of the techniques used here can be used for analyzing other forms of splicing; however, our automated pipeline has been tuned for IR detection. Many of the cleaning and auto-detection tools described below may seem redundant given the extensive annotation of publicly available sequencing data. As we found, however, numerous samples that we collected weren’t cleaned extensively and in some cases may have been misannotated.
A STAR  genome index is built with a user supplied genome fasta file and annotation gtf file. An automated process allows ready use of other Ensembl genomes. All potential introns are extracted from the gtf file, being the region between two exons in any transcript. Regions covered by a gtf feature within each intron are then excluded as they are likely to confound accurate measurement of the true intron level. Excluded features are all annotations in the gtf file except those marked “retained_intron”. For directional sequencing, only features on the same strand as the intron are excluded. For non-directional sequencing, exclusions are omnidirectional.
Regions of poor unique mappability are determined by mapping synthetic reads to the genome. Synthetic reads are 70 bp single end, stepped at 10 bp across the entire reference genome. Every second read is reverse complemented. A single base error is generated in the center of each read. This error is generated in a deterministic manner, allowing reproducible results when the same input files and the same software versions are used. Reads uniquely mapping to the correct location are tallied. Any 70-bp stretch without at least five unique reads is considered poorly mappable. Poorly mappable regions are excluded from the measurable intron area regardless of strand/direction.
Data preparation and quality controls
Illumina sequencing, when performed on a standard paired-end library, generates two reads commencing at opposite ends of the insert and proceeding towards each other. If either read is longer than the insert, the read will continue into the sequencing adaptor. The adaptor is not reached until the entire insert is sequenced. As such, for pairs with adaptor contamination, each read commences with the insert as a reverse complementary intersection and completes with adaptor sequence at the 3′ end. We made use of this feature to automatically detect adaptors and trim them.
Our automatic detection algorithm takes the first 250,000 read pairs and performs a gapless alignment of the two reads against each other. Alignments considered have a reverse-complementary part commencing at the 5′ end with an overhanging 3′ end. From pairs with a best alignment of at least 90%, the non-overlapping components are stored as potential adaptors. This list of potential adaptors is analyzed independently for both forward and reverse reads.
Automatic detection of adaptors is implemented in a standalone PERL program and is available in the main IRFinder package. It can be readily used without dependencies both as part of this package or standalone.
Trimming of the adaptor, once identified, is conducted as a streaming process along with IRFinder’s mapping, count, and sort functions. The streaming process ensures unnecessary temporary files are not produced, saving both disk space and the disk I/O performance impact. Further, the streaming design substantially reduces real time to result.
Trimming uses the known adaptor sequence expected on each pair along with the complementary overlapping portion of the reads. Use of both the known adaptor sequences and the overlapping complementary section allows precise trimming, even when only one base of an adaptor is present. Trimming algorithms not using the reverse complementary segment have insufficient information for accuracy and thus must either over- or under-trim fragments containing only a short length of adaptor.
Where sequencing data are single-end only, STAR’s built-in trimming function is utilized.
Adaptor trimming is implemented as a C++ program. It can be used as an integrated part of this pipeline or standalone.
Auto-detection of directionality
Antisense RNA can confound the calculation of IR levels, so we developed a method to automatically detect if a library was prepared using a directional protocol such as ScriptSeq or using dUTPs. These protocols enable bioinformatics analyses to correctly attribute reads to a transcript or an overlapping antisense transcript. To detect directionality, IRFinder measures the coverage across splice junctions. For each splice junction crossed by more than eight reads, if one direction is more than fourfold the other, it is counted as evidence of directionality; if not, it is evidence against directionality. If the directional score is at least 90%, the directional analysis is output for this sample. In practice the directional score is well over 99% for directional data.
Quality control of the sample
IRFinder automatically detects samples that are not suitable for IR analysis. These samples either have high levels of DNA contamination or have been mislabeled as mRNA sequencing when in fact they are other types of experiments such as genome re-sequencing or ChIP-seq. Both DNA contamination and incorrectly labeled samples can be detected by calculating the ratio of the number of reads that map to intergenic regions to the number that maps to coding regions. If this ratio is higher than 10%, IRFinder emits a warning that this sample may not be suitable for IR detection.
To ensure reads that map to introns do not come from unprocessed mRNA, the sequencing libraries must be enriched for polyadenylated RNA. To detect RNA sequencing experiments for which the library was not enriched for mature, polyadenylated mRNAs, IRFinder counts the number of reads that map to a list of non-polyadenylated genes such as small nucleolar RNAs or specific histone genes (Additional file 9: Table S8). If the sample has not been poly(A) enriched, this read count will be much higher than expected as shown by a comparison of poly(A)- and non-poly(A)-enriched libraries (Additional file 9: Figure S6, Table S9). IRFinder will warn the user if this read count is higher than 0.01%. This threshold clearly distinguishes poly(A)- from non-poly(A)-enriched samples.
Measuring intron retention
Reads are mapped to the reference genome by STAR using default scoring parameters but excluding multi-mapping reads from the output. The default scoring parameters suit the detection of IR as they favor neither mapping from exon to exon across splice junctions nor mapping from exons into introns.
Where intron coverage is less than 1, the proportion of bases covered is used to calculate intronic abundance. The proportion of bases cannot be used as a surrogate for measuring IR ratios but is rather a means of avoiding null values when comparing an intron across multiple samples. We recommend filtering out IR candidates with coverage less than three reads across the entire measurable intron.
Detecting changes in IR between samples
Where a pair of samples or groups is given as input to IRFinder, the final output is a listing of introns most significantly changed. Each intron can be evaluated using a Bayesian statistic adapted for digital counts  or DESeq2  for which we have coded direct plug-ins. Significance of the change in IR ratio will be displayed in the results. Other statistics can easily be implemented on IRFinder output given the extensive information provided.
To validate the results of our algorithm we used RT-qPCR to measure intron retention in two cases where IR detection can be problematic. The first validation was performed on low IR levels, the second on a long intron with regions of low mappability. Total RNA was extracted from cells using Trizol (Ambion), treated with DNaseI (Life Technologies), and converted to cDNA using Superscript III (Life Technologies). Quantification of intronic expression shown in Additional file 1 was performed with normalization against the expression of adjacent exons using RT-qPCR on the CFX96 Real-Time System (Biorad). Average deltaCT values are graphed for each cell line from three experiments. Sequences and locations of forward and reverse primers used to amplify the selected retained introns are listed in Additional file 1.
Comparison with existing tools
To compare the efficiency of IRFinder with other existing tools that can detect IR events, we measured the average change in expression of retained introns predicted by IRFinder, MISO , and DEXseq  in a Upf2 knockout model .
Total RNA-seq of liver tissue from wild-type (Upf2 fl/fl ) and Upf2 knockout (Upf2 fl/fl ;Mx1Cre) mouse (GSE26561) was downloaded from the Gene Expression Omnibus (GEO). The single-end reads from both samples were mapped to the transcriptome annotation version 73 of Ensembl mouse genome GRCm38 using STAR (version and options described in IRFinder methods).
IRFinder with default parameters was applied to determine differentially retained introns between the wild-type and UPF2 knockout samples.
DEXSeq 1.14.0 with default parameters was applied to determine differential intron usage between the wild-type and UPF2 knockout samples. Each intron with q value (false discovery rate-adjusted p value) less than 0.05 in the DEXSeq report was recognized as differentially retained.
MISO 0.5.3 with default parameters was applied to determine differentially retained introns between the wild-type and UPF2 knockout samples. Each differentially retained intron in the MISO report met the following criteria: 1) the Bayesian factor was above 19 (the likelihood that this intron is retained is 19 times higher than that of not being retained); 2) at least one read covered this intronic region; 3) at least ten reads covered the two flanking exons of this intron.
We thank Annie Varrault, Tristan Bouschet, and Laurent Journot for their input on the manuscript figures and text.
This work was supported by the Agence Nationale de la Recherche (ANR 143683); the National Health and Medical Research Council (grant #1061906, #1080530, #1128175, #1126306). BS and EE were supported by grants BIO2014-52566-R from the MINECO (Spanish Government) and FEDER funds, by AGAUR (2014-SGR1121), and by the Sandra Ibarra Foundation for Cancer (FSI2013).
Availability of data and requirements
IRFinder is freely available under the MIT License with DOI doi.org/10.5281/zenodo.265180 at https://github.com/williamritchie/IRFinder. IRFinder requires a UNIX based system with standard LINUX utilities.
WR and RM designed the IRFinder software; RM coded it. AT tested, found limitations, suggested modifications, and maintained the software. DG created the database. WR designed the study. JR, JW, AA, EE, BS, BC, AB, DG, and WR performed experiments, generated figures, and helped to write the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
This study required no ethics approval.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet. 2008;40:1413–5.View ArticlePubMedGoogle Scholar
- Buckley PT, Lee MT, Sul J-Y, Miyashiro KY, Bell TJ, Fisher SA, Kim J, Eberwine J. Cytoplasmic intron sequence-retaining transcripts can be dendritically targeted via ID element retrotransposons. Neuron. 2011;69:877–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Belgrader P, Cheng J, Zhou X, Stephenson LS, Maquat LE. Mammalian nonsense codons can be cis effectors of nuclear mRNA half-life. Mol Cell Biol. 1994;14:8219–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Baek D, Green P. Sequence conservation, relative isoform frequencies, and nonsense-mediated decay in evolutionarily conserved alternative splicing. Proc Natl Acad Sci U S A. 2005;102:12813–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Hwang J, Maquat LE. Nonsense-mediated mRNA decay (NMD) in animal embryogenesis: to die or not to die, that is the question. Curr Opin Genet Dev. 2011;21:422–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Wong JJ, Ritchie W, Ebner OA, Selbach M, Wong JW, Huang Y, Gao D, Pinello N, Gonzalez M, Baidya K, et al. Orchestrated intron retention regulates normal granulocyte differentiation. Cell. 2013;154:583–95.View ArticlePubMedGoogle Scholar
- Hussein SM, Puri MC, Tonge PD, Benevento M, Corso AJ, Clancy JL, Mosbergen R, Li M, Lee DS, Cloonan N, et al. Genome-wide characterization of the routes to pluripotency. Nature. 2014;516:198–206.View ArticlePubMedGoogle Scholar
- Braunschweig U, Barbosa-Morais NL, Pan Q, Nachman EN, Alipanahi B, Gonatopoulos-Pournatzis T, Frey B, Irimia M, Blencowe BJ. Widespread intron retention in mammals functionally tunes transcriptomes. Genome Res. 2014;24:1774–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Llorian M, Gooding C, Bellora N, Hallegger M, Buckroyd A, Wang X, Rajgor D, Kayikci M, Feltham J, Ule J, et al. The alternative splicing program of differentiated smooth muscle cells involves concerted non-productive splicing of post-transcriptional regulators. Nucleic Acids Res. 2016;44(18):8933–8950.Google Scholar
- Jung H, Lee D, Lee J, Park D, Kim YJ, Park WY, Hong D, Park PJ, Lee E. Intron retention is a widespread mechanism of tumor-suppressor inactivation. Nat Genet. 2015;47:1242–8.View ArticlePubMedGoogle Scholar
- Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson A, Kampf C, Sjostedt E, Asplund A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347:1260419.View ArticlePubMedGoogle Scholar
- Gao X, Wan J, Liu B, Ma M, Shen B, Qian SB. Quantitative profiling of initiating ribosomes in vivo. Nat Methods. 2015;12:147–53.View ArticlePubMedGoogle Scholar
- Edwards CR, Ritchie W, Wong JJ, Schmitz U, Middleton R, An X, Mohandas N, Rasko JE, Blobel GA. A dynamic intron retention program in the mammalian megakaryocyte and erythrocyte lineages. Blood. 2016;127:e24–e34.Google Scholar
- Choi J, Lee S, Mallard W, Clement K, Tagliazucchi GM, Lim H, Choi IY, Ferrari F, Tsankov AM, Pop R, et al. A comparison of genetically matched cell lines reveals the equivalence of human iPSCs and ESCs. Nat Biotechnol. 2015;33:1173–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Buettner F, Natarajan KN, Casale FP, Proserpio V, Scialdone A, Theis FJ, Teichmann SA, Marioni JC, Stegle O. Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells. Nat Biotechnol. 2015;33:155–60.View ArticlePubMedGoogle Scholar
- Busskamp V, Lewis NE, Guye P, Ng AH, Shipman SL, Byrne SM, Sanjana NE, Murn J, Li Y, Li S, et al. Rapid neurogenesis through transcriptional activation in human stem cells. Mol Syst Biol. 2014;10:760.View ArticlePubMedPubMed CentralGoogle Scholar
- Pimentel H, Parra M, Gee SL, Mohandas N, Pachter L, Conboy JG. A dynamic intron retention program enriched in RNA processing genes regulates gene expression during terminal erythropoiesis. Nucleic Acids Res. 2016;44:838–51.View ArticlePubMedGoogle Scholar
- Maruyama M, Ichisaka T, Nakagawa M, Yamanaka S. Differential roles for Sox15 and Sox2 in transcriptional control in mouse embryonic stem cells. J Biol Chem. 2005;280:24371–9.View ArticlePubMedGoogle Scholar
- Dykes IM, Tempest L, Lee SI, Turner EE. Brn3a and Islet1 act epistatically to regulate the gene expression program of sensory differentiation. J Neurosci. 2011;31:9789–99.View ArticlePubMedPubMed CentralGoogle Scholar
- Lembach KJ. Induction of human fibroblast proliferation by epidermal growth factor (EGF): enhancement by an EGF-binding arginine esterase and by ascorbate. Proc Natl Acad Sci U S A. 1976;73:183–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Tian L, Chen L, McClafferty H, Sailer CA, Ruth P, Knaus HG, Shipston MJ. A noncanonical SH3 domain binding motif links BK channels to the actin cytoskeleton via the SH3 adapter cortactin. FASEB J. 2006;20:2588–90.View ArticlePubMedGoogle Scholar
- Ni T, Yang W, Han M, Zhang Y, Shen T, Nie H, Zhou Z, Dai Y, Yang Y, Liu P, et al. Global intron retention mediated gene regulation during CD4+ T cell activation. Nucleic Acids Res. 2016;44:6817–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Del Gatto-Konczak F, Bourgeois CF, Le Guiner C, Kister L, Gesnel MC, Stevenin J, Breathnach R. The RNA-binding protein TIA-1 is a novel mammalian splicing regulator acting through intron sequences adjacent to a 5′ splice site. Mol Cell Biol. 2000;20:6287–99.View ArticlePubMedPubMed CentralGoogle Scholar
- Linares AJ, Lin CH, Damianov A, Adams KL, Novitch BG, Black DL. The splicing regulator PTBP1 controls the activity of the transcription factor Pbx1 during neuronal differentiation. Elife. 2015;4:e09268.View ArticlePubMedPubMed CentralGoogle Scholar
- Rogelj B, Easton LE, Bogu GK, Stanton LW, Rot G, Curk T, Zupan B, Sugimoto Y, Modic M, Haberman N, et al. Widespread binding of FUS along nascent RNA regulates alternative splicing in the brain. Sci Rep. 2012;2:603.View ArticlePubMedPubMed CentralGoogle Scholar
- Twyffels L, Gueydan C, Kruys V. Shuttling SR proteins: more than splicing factors. FEBS J. 2011;278:3246–55.View ArticlePubMedGoogle Scholar
- Lareau LF, Inada M, Green RE, Wengrod JC, Brenner SE. Unproductive splicing of SR genes associated with highly conserved and ultraconserved DNA elements. Nature. 2007;446:926–9.View ArticlePubMedGoogle Scholar
- Yap K, Lim ZQ, Khandelia P, Friedman B, Makeyev EV. Coordinated regulation of neuronal mRNA steady-state levels through developmentally controlled intron retention. Genes Dev. 2012;26:1209–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Boutz PL, Bhutkar A, Sharp PA. Detained introns are a novel, widespread class of post-transcriptionally spliced introns. Genes Dev. 2015;29:63–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Feng Q, Snider L, Jagannathan S, Tawil R, van der Maarel SM, Tapscott SJ and Bradley RK. A feedback loop between nonsense-mediated decay and the retrogene DUX4 in facioscapulohumeral muscular dystrophy. eLife. 2015;4:e04996.Google Scholar
- Singh G, Kucukural A, Cenik C, Leszyk JD, Shaffer SA, Weng Z, Moore MJ. The cellular EJC interactome reveals higher-order mRNP structure and an EJC-SR protein nexus. Cell. 2012;151:750–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.View ArticlePubMedGoogle Scholar
- Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7:986–95.PubMedGoogle Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.View ArticlePubMedPubMed CentralGoogle Scholar
- Katz Y, Wang ET, Airoldi EM, Burge CB. Analysis and design of RNA sequencing experiments for identifying isoform regulation. Nat Methods. 2010;7:1009–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Anders S, Reyes A, Huber W. Detecting differential usage of exons from RNA-seq data. Genome Res. 2012;22:2008–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Weischenfeldt J, Waage J, Tian G, Zhao J, Damgaard I, Jakobsen JS, Kristiansen K, Krogh A, Wang J, Porse BT. Mammalian tissues defective in nonsense-mediated mRNA decay display highly aberrant splicing patterns. Genome Biol. 2012;13:R35.View ArticlePubMedPubMed CentralGoogle Scholar