The RNA-binding protein hnRNPLL induces a T cell alternative splicing program delineated by differential intron retention in polyadenylated RNA

Background

Retention of a subset of introns in spliced polyadenylated mRNA is emerging as a frequent, unexplained finding from RNA deep sequencing in mammalian cells.

Results

Here we analyze intron retention in T lymphocytes by deep sequencing polyadenylated RNA. We show a developmentally regulated RNA-binding protein, hnRNPLL, induces retention of specific introns by sequencing RNA from T cells with an inactivating Hnrpll mutation and from B lymphocytes that physiologically downregulate Hnrpll during their differentiation. In Ptprc mRNA encoding the tyrosine phosphatase CD45, hnRNPLL induces selective retention of introns flanking exons 4 to 6; these correspond to the cassette exons containing hnRNPLL binding sites that are skipped in cells with normal, but not mutant or low, hnRNPLL. We identify similar patterns of hnRNPLL-induced differential intron retention flanking alternative exons in 14 other genes, representing novel elements of the hnRNPLL-induced splicing program in T cells. Retroviral expression of a normally spliced cDNA for one of these targets, Senp2, partially corrects the survival defect of Hnrpll-mutant T cells. We find that integrating a number of computational methods to detect genes with differentially retained introns provides a strategy to enrich for alternatively spliced exons in mammalian RNA-seq data, when complemented by RNA-seq analysis of purified cells with experimentally perturbed RNA-binding proteins.

Conclusions

Our findings demonstrate that intron retention in mRNA is induced by specific RNA-binding proteins and suggest a biological significance for this process in marking exons that are poised for alternative splicing.     

A class of human exons with predicted distant branch points revealed by analysis of AG dinucleotide exclusion zones

Background

The three consensus elements at the 3' end of human introns - the branch point sequence, the polypyrimidine tract, and the 3' splice site AG dinucleotide - are usually closely spaced within the final 40 nucleotides of the intron. However, the branch point sequence and polypyrimidine tract of a few known alternatively spliced exons lie up to 400 nucleotides upstream of the 3' splice site. The extended regions between the distant branch points (dBPs) and their 3' splice site are marked by the absence of other AG dinucleotides. In many cases alternative splicing regulatory elements are located within this region.

Results

We have applied a simple algorithm, based on AG dinucleotide exclusion zones (AGEZ), to a large data set of verified human exons. We found a substantial number of exons with large AGEZs, which represent candidate dBP exons. We verified the importance of the predicted dBPs for splicing of some of these exons. This group of exons exhibits a higher than average prevalence of observed alternative splicing, and many of the exons are in genes with some human disease association.

Conclusion

The group of identified probable dBP exons are interesting first because they are likely to be alternatively spliced. Second, they are expected to be vulnerable to mutations within the entire extended AGEZ. Disruption of splicing of such exons, for example by mutations that lead to insertion of a new AG dinucleotide between the dBP and 3' splice site, could be readily understood even though the causative mutation might be remote from the conventional locations of splice site sequences.     

Transcriptome-wide modulation of splicing by the exon junction complex

Background

The exon junction complex (EJC) is a dynamic multi-protein complex deposited onto nuclear spliced mRNAs upstream of exon-exon junctions. The four core proteins, eIF4A3, Magoh, Y14 and MLN51, are stably bound to mRNAs during their lifecycle, serving as a binding platform for other nuclear and cytoplasmic proteins. Recent evidence has shown that the EJC is involved in the splicing regulation of some specific events in both Drosophila and mammalian cells.

Results

Here, we show that knockdown of EJC core proteins causes widespread alternative splicing changes in mammalian cells. These splicing changes are specific to EJC core proteins, as knockdown of eIF4A3, Y14 and MLN51 shows similar splicing changes, and are different from knockdown of other splicing factors. The splicing changes can be rescued by a siRNA-resistant form of eIF4A3, indicating an involvement of EJC core proteins in regulating alternative splicing. Finally, we find that the splicing changes are linked with RNA polymerase II elongation rates.

Conclusion

Taken together, this study reveals that the coupling between EJC proteins and splicing is broader than previously suspected, and that a possible link exists between mRNP assembly and splice site recognition.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0551-7) contains supplementary material, which is available to authorized users.     

HRP-2, the Caenorhabditis elegans Homolog of Mammalian Heterogeneous Nuclear Ribonucleoproteins Q and R, Is an Alternative Splicing Factor That Binds to UCUAUC Splicing Regulatory Elements

Alternative splicing is regulated by cis sequences in the pre-mRNA that serve as binding sites for trans-acting alternative splicing factors. In a previous study, we used bioinformatics and molecular biology to identify and confirm that the intronic hexamer sequence UCUAUC is a nematode alternative splicing regulatory element. In this study, we used RNA affinity chromatography to identify trans factors that bind to this sequence. HRP-2, the Caenorhabditis elegans homolog of human heterogeneous nuclear ribonucleoproteins Q and R, binds to UCUAUC in the context of unc-52 intronic regulatory sequences as well as to RNAs containing tandem repeats of this sequence. The three Us in the hexamer are the most important determinants of this binding specificity. We demonstrate, using RNA interference, that HRP-2 regulates the alternative splicing of two genes, unc-52 and lin-10, both of which have cassette exons flanked by an intronic UCUAUC motif. We propose that HRP-2 is a protein responsible for regulating alternative splicing through binding interactions with the UCUAUC sequence.Alternative pre-mRNA splicing is a mechanism for generating multiple mRNA isoforms from a single gene. This process can allow a gene to encode for more than one protein isoform. For some genes, it is a mechanism for regulating message stability through production of alternative mRNA isoforms that are substrates for the nonsense-mediated mRNA decay pathway (1). The majority of human genes undergo alternative splicing (2), and the process can be regulated in tissue-specific and developmental stage-specific manners. Current models propose that cis elements on the pre-mRNA, in exons and introns, serve as recognition sites for trans-acting protein factors that bind to the pre-mRNA and regulate assembly of the splicing machinery, thus regulating splice site choice (3).In recent years, a number of groups have employed bioinformatics techniques to identify cis splicing regulatory elements (4). These techniques include using multiple interspecies sequence alignments to identify conserved intronic regions, identification of short sequences in exons that are bounded by weak consensus splice sites, and identification of common intronic sequences flanking similarly regulated alternative exons (5–9). These efforts have added many new sequences to the list of known and potential splicing regulators. The identification of the protein factor partners for these sequences will be important for understanding their function in alternative splicing regulation.Experimental approaches have identified alternative splicing factors that interact with specific cis elements (10), but the number of trans factors discovered still lags behind the number of newly identified cis element partners. Some examples of well-characterized cis element/trans-acting factor interactions include the NOVA K homology domain splicing factor binding to the sequence UCAY (11), the FOX splicing factors binding to the sequence UGCAUG (12–14), and hnRNP³ F/H proteins binding to the sequence GGGG (15, 16). By using cross-linking immunoprecipitation followed by large scale sequencing, entire catalogs of RNAs that the splicing factors NOVA, SF2/ASF, and FOX2 bind to in vivo have been determined (17–19). These approaches have led to models for how the proteins binding to their cis regulatory elements may alter splicing. These models include a role for the relative position of a cis element to an alternative cassette exon in determining alternative exon inclusion or skipping (18, 19).In a previous bioinformatics analysis of evolutionarily conserved intronic sequences flanking alternatively spliced exons, we identified the hexamer sequence UCUAUC as a novel splicing regulatory element (8). UCUAUC is found flanking both sides of alternative exon 16 of the unc-52 gene of Caenorhabditis elegans. Genetic analysis of a class of viable unc-52 mutants led to the discovery that exons 16–18 are alternative cassette exons and that every combination of skipping and inclusion of these three exons occurs (20). This splicing is regulated by the alternative splicing factor MEC-8 (21). Fig. 1A shows a schematic diagram of the alternatively spliced region of unc-52, with the MEC-8-enhanced alternative splicing events indicated. Using an unc-52 splicing reporter trans gene containing alternative exons 15–19, we previously reported that alternative splicing is regulated by the intronic motif UCUAUC in the intron downstream of exon 16 (8). In addition we showed that this element works cooperatively with a UGCAUG hexamer (the consensus FOX-1-binding site) in the upstream intron to regulate alternative splicing (8).Open in a separate windowFIGURE 1.RNA affinity chromatography identifies HRP-2 as binding to UCUAUC elements. A, schematic representation of the alternatively spliced region of unc-52 (adapted from Ref. 21). The alternative splicing events promoted by MEC-8 are indicated by bold lines. The lines next to introns 15 and 16 are the sites of the UCUAUC elements in those introns whose sequences were used in the RNA affinity chromatography. B, table showing sequences of RNAs immobilized to beads in the RNA affinity chromatography experiment. C, Coomassie-stained SDS-PAGE analysis of RNA affinity chromatography. C. elegans embryo extract was incubated with the different immobilized RNA substrates listed on top of the gel. Proteins identified by mass spectrometry are listed to the right of the gel, with arrows pointing to coincident protein bands. D, the left panel shows the silver stain result for the RNA affinity chromatography experiment. Each lane represents a different immobilized substrate, as indicated above. The band corresponding to HRP-2 is indicated by an arrow. The right panel is an immunoblot of the same gel using anti-HRP-2 polyclonal antibody. E, anti-HRP-2 immunoblot of an RNA affinity chromatography experiment for the indicated substrates.In this study, we report the results of a biochemical identification of a protein factor from C. elegans that binds to the UCUAUC intronic splicing regulatory element. We transcribed different short RNA sequences containing the UCUAUC element in its native intronic context, or as part of a repeating unit, and immobilized these onto agarose beads. After passing embryo extracts across these beads, we found that the protein HRP-2, the C. elegans homolog of the mammalian hnRNP Q/R proteins, binds to this sequence with high affinity. By using RNAi to reduce the level of HRP-2 in worms, we observed changes in alternative splicing of unc-52 and lin-10, two genes that contain UCUAUC elements in introns flanking alternative exons. We propose that HRP-2 is an alternative splicing factor that works through the UCUAUC intronic elements to regulate alternative splicing.     

Functional coordination of alternative splicing in the mammalian central nervous system

Background

Alternative splicing (AS) functions to expand proteomic complexity and plays numerous important roles in gene regulation. However, the extent to which AS coordinates functions in a cell and tissue type specific manner is not known. Moreover, the sequence code that underlies cell and tissue type specific regulation of AS is poorly understood.

Results

Using quantitative AS microarray profiling, we have identified a large number of widely expressed mouse genes that contain single or coordinated pairs of alternative exons that are spliced in a tissue regulated fashion. The majority of these AS events display differential regulation in central nervous system (CNS) tissues. Approximately half of the corresponding genes have neural specific functions and operate in common processes and interconnected pathways. Differential regulation of AS in the CNS tissues correlates strongly with a set of mostly new motifs that are predominantly located in the intron and constitutive exon sequences neighboring CNS-regulated alternative exons. Different subsets of these motifs are correlated with either increased inclusion or increased exclusion of alternative exons in CNS tissues, relative to the other profiled tissues.

Conclusion

Our findings provide new evidence that specific cellular processes in the mammalian CNS are coordinated at the level of AS, and that a complex splicing code underlies CNS specific AS regulation. This code appears to comprise many new motifs, some of which are located in the constitutive exons neighboring regulated alternative exons. These data provide a basis for understanding the molecular mechanisms by which the tissue specific functions of widely expressed genes are coordinated at the level of AS.