Architectural limits on split genes

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Regulation of alternative splicing by retrograde and light signals converges to control chloroplast proteins

Retrograde signals sent by chloroplasts control transcription in the nucleus. These signals antagonistically converge with light signals to coordinate the expression of genes involved in chloroplast functioning and seedling development. Although significant advances have been made in understanding the molecular interplay between light and retrograde signals at the transcriptional level, little is known about their interconnection at the post-transcriptional level. By using different publicly available datasets, this study addresses the influence of retrograde signaling on alternative splicing and defines the molecular and biological functions of this regulation. These analyses revealed that alternative splicing mimics transcriptional responses triggered by retrograde signals at different levels. First, both molecular processes similarly depend on the chloroplast-localized pentatricopeptide-repeat protein GUN1 to modulate the nuclear transcriptome. Secondly, as described for transcriptional regulation, alternative splicing coupled with the nonsense-mediated decay pathway effectively downregulates expression of chloroplast proteins in response to retrograde signals. Finally, light signals were found to antagonistically control retrograde signaling-regulated splicing isoforms, which consequently generates opposite splicing outcomes that likely contribute to the opposite roles these signals play in controlling chloroplast functioning and seedling development.

Selection of Olduvai Domains during Evolution: A Role for Primate-Specific Splicing Super-Enhancer and RNA Guanine Quadruplex in Bipartite NBPF Exons

Olduvai protein domains (also known as DUF1220 or NBPF) have undergone the greatest human-specific increase in the copy number of any coding region in the genome. Their repeat number was strongly associated with the evolutionary expansion of brain volumes, neuron counts and cognitive abilities, as well as with disorders of the autistic spectrum. Nevertheless, the domain function and cellular mechanisms underlying the positive selection of Olduvai DNA sequences in higher primates remain obscure. Here, I show that the inclusion of Olduvai exon doublets in mature transcripts is facilitated by a potent splicing enhancer that was created through duplication within the first exon. The enhancer is the strongest among the NBPF transcripts and further promotes the already high splicing activity of the unexpanded first exons of the two-exon domains, safeguarding the expanded Olduvai exon doublets in the mature transcriptome. The duplication also creates a predicted RNA guanine quadruplex that may regulate the access to spliceosomal components of the super-enhancer and influence the splicing of adjacent exons. Thus, positive Olduvai selection during primate evolution is likely to result from a combination of multiple targets in gene expression pathways, including RNA splicing.

Gene architecture directs splicing outcome in separate nuclear spatial regions

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

Transcription and splicing dynamics during early Drosophila development

Widespread cotranscriptional splicing has been demonstrated from yeast to human. However, most studies to date addressing the kinetics of splicing relative to transcription used either Saccharomyces cerevisiae or metazoan cultured cell lines. Here, we adapted native elongating transcript sequencing technology (NET-seq) to measure cotranscriptional splicing dynamics during the early developmental stages of Drosophila melanogaster embryos. Our results reveal the position of RNA polymerase II (Pol II) when both canonical and recursive splicing occur. We found heterogeneity in splicing dynamics, with some RNAs spliced immediately after intron transcription, whereas for other transcripts no splicing was observed over the first 100 nt of the downstream exon. Introns that show splicing completion before Pol II has reached the end of the downstream exon are necessarily intron-defined. We studied the splicing dynamics of both nascent pre-mRNAs transcribed in the early embryo, which have few and short introns, as well as pre-mRNAs transcribed later in embryonic development, which contain multiple long introns. As expected, we found a relationship between the proportion of spliced reads and intron size. However, intron definition was observed at all intron sizes. We further observed that genes transcribed in the early embryo tend to be isolated in the genome whereas genes transcribed later are often overlapped by a neighboring convergent gene. In isolated genes, transcription termination occurred soon after the polyadenylation site, while in overlapped genes, Pol II persisted associated with the DNA template after cleavage and polyadenylation of the nascent transcript. Taken together, our data unravel novel dynamic features of Pol II transcription and splicing in the developing Drosophila embryo.

Exon-independent recruitment of SRSF1 is mediated by U1 snRNP stem-loop 3

SRSF1 protein and U1 snRNPs are closely connected splicing factors. They both stimulate exon inclusion, SRSF1 by binding to exonic splicing enhancer sequences (ESEs) and U1 snRNPs by binding to the downstream 5' splice site (SS), and both factors affect 5' SS selection. The binding of U1 snRNPs initiates spliceosome assembly, but SR proteins such as SRSF1 can in some cases substitute for it. The mechanistic basis of this relationship is poorly understood. We show here by single-molecule methods that a single molecule of SRSF1 can be recruited by a U1 snRNP. This reaction is independent of exon sequences and separate from the U1-independent process of binding to an ESE. Structural analysis and cross-linking data show that SRSF1 contacts U1 snRNA stem-loop 3, which is required for splicing. We suggest that the recruitment of SRSF1 to a U1 snRNP at a 5'SS is the basis for exon definition by U1 snRNP and might be one of the principal functions of U1 snRNPs in the core reactions of splicing in mammals.

Splice site proximity influences alternative exon definition

Alternative splicing enables higher eukaryotes to expand mRNA diversity from a finite number of genes through highly combinatorial splice site selection mechanisms that are influenced by the sequence of competing splice sites, cis-regulatory elements binding trans-acting factors, the length of exons and introns harbouring alternative splice sites and RNA secondary structures at putative splice junctions. To test the hypothesis that the intron definition or exon definition modes of splice site recognition direct the selection of alternative splice patterns, we created a database of alternative splice site usage (ALTssDB). When alternative splice sites are embedded within short introns (intron definition), the 5' and 3' splice sites closest to each other across the intron preferentially pair, consistent with previous observations. However, when alternative splice sites are embedded within large flanking introns (exon definition), the 5' and 3' splice sites closest to each other across the exon are preferentially selected. Thus, alternative splicing decisions are influenced by the intron and exon definition modes of splice site recognition. The results demonstrate that the spliceosome pairs splice sites that are closest in proximity within the unit of initial splice site selection.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

Spliceosomal component PRP-40 is a central regulator of microexon splicing

Microexons (≤27 nt) play critical roles in nervous system development and function but create unique challenges for the splicing machinery. The mechanisms of microexon regulation are therefore of great interest. We performed a genetic screen for alternative splicing regulators in the C. elegans nervous system and identify PRP-40, a core component of the U1 snRNP. RNA-seq reveals that PRP-40 is required for inclusion of alternatively spliced, but not constitutively spliced, exons. PRP-40 is particularly required for inclusion of neuronal microexons, and our data indicate that PRP-40 is a central regulator of microexon splicing. Microexons can be relieved from PRP-40 dependence by artificially increasing exon size or reducing flanking intron size, indicating that PRP-40 is specifically required for microexons surrounded by conventionally sized introns. Knockdown of the orthologous PRPF40A in mouse neuroblastoma cells causes widespread dysregulation of microexons but not conventionally sized exons. PRP-40 regulation of neuronal microexons is therefore a widely conserved phenomenon.

Induction of cryptic pre-mRNA splice-switching by antisense oligonucleotides

Antisense oligomers (AOs) are increasingly being used to modulate RNA splicing in live cells, both for research and for the development of therapeutics. While the most common intended effect of these AOs is to induce skipping of whole exons, rare examples are emerging of AOs that induce skipping of only part of an exon, through activation of an internal cryptic splice site. In this report, we examined seven AO-induced cryptic splice sites in six genes. Five of these cryptic splice sites were discovered through our own experiments, and two originated from other published reports. We modelled the predicted effects of AO binding on the secondary structure of each of the RNA targets, and how these alterations would in turn affect the accessibility of the RNA to splice factors. We observed that a common predicted effect of AO binding was disruption of the exon definition signal within the exon's excluded segment.

Identification of prognostic and bone metastasis​-related alternative splicing signatures in mesothelioma

Mesothelioma (MESO) is an infrequent tumor derived from mesothelial cells of pleura, peritoneum, pericardium, and tunica vaginalis testis. Despite advancement in technologies and better understanding of tumor progression mechanism, the prognosis of MESO remains poor. The role of alternative splicing events (ASEs) in the oncogenesis, tumor metastasis and drug resistance has been widely discussed in multiple cancers. But the prognosis and potential therapeutic value of ASEs in MESO were not clearly studied by now. We constructed a prognostic model using RNA sequencing data and matched ASE data of MESO patients obtained from the TCGA and TCGASpliceSeq database. A total of 3,993 ASEs were identified associated with overall survival using Cox regression analysis. Eight of them were finally figured out to institute the model by lasso regression analysis. The risk score of the model can predict the prognosis independently. Among the identified 390 splicing factors (SF), HSPA1A and DDX3Y was significantly associated with 43 OS-SEs. Among these OS-SEs, SNX5-58744-AT (p = 0.048) and SNX5-58745-AT (p = 0.048) were significantly associated with bone metastasis. Co-expression analysis of signal pathways and SNX5-58744-AT, SNX5-58745-AT was also depicted using GSVA. Finally, we proposed that splicing factor (SF) HSPA1A could regulate SNX5-58744-AT (R = -0.414) and SNX5-58745-AT (R = 0.414) through the pathway "Class I MHC mediated antigen processing and presentation" (R = 0.400). In this way, tumorigenesis and bone metastasis of MESO were controlled.

Widespread Alternative Splicing Changes in Metastatic Breast Cancer Cells

Aberrant alternative splicing (AS) is a hallmark of cancer and a potential target for novel anti-cancer therapeutics. Breast cancer-associated AS events are known to be linked to disease progression, metastasis, and survival of breast cancer patients. To identify altered AS programs occurring in metastatic breast cancer, we perform a global analysis of AS events by using RNA-mediated oligonucleotide annealing, selection, and ligation coupled with next-generation sequencing (RASL-seq). We demonstrate that, relative to low-metastatic, high-metastatic breast cancer cells show different AS choices in genes related to cancer progression. Supporting a global reshape of cancer-related splicing profiles in metastatic breast cancer we found an enrichment of RNA-binding motifs recognized by several splicing regulators, which have aberrant expression levels or activity during breast cancer progression, including SRSF1. Among SRSF1-regulated targets we found DCUN1D5, a gene for which skipping of exon 4 in its pre-mRNA introduces a premature termination codon (PTC), thus generating an unstable transcript degraded by nonsense-mediated mRNA decay (NMD). Significantly, distinct breast cancer subtypes show different DCUN1D5 isoform ratios with metastatic breast cancer expressing the highest level of the NMD-insensitive DCUN1D5 mRNA, thus showing high DCUN1D5 expression levels, which are ultimately associated with poor overall and relapse-free survival in breast cancer patients. Collectively, our results reveal global AS features of metastatic breast tumors, which open new possibilities for the treatment of these aggressive tumor types.

The Alter Retina: Alternative Splicing of Retinal Genes in Health and Disease

Alternative splicing of mRNA is an essential mechanism to regulate and increase the diversity of the transcriptome and proteome. Alternative splicing frequently occurs in a tissue- or time-specific manner, contributing to differential gene expression between cell types during development. Neural tissues present extremely complex splicing programs and display the highest number of alternative splicing events. As an extension of the central nervous system, the retina constitutes an excellent system to illustrate the high diversity of neural transcripts. The retina expresses retinal specific splicing factors and produces a large number of alternative transcripts, including exclusive tissue-specific exons, which require an exquisite regulation. In fact, a current challenge in the genetic diagnosis of inherited retinal diseases stems from the lack of information regarding alternative splicing of retinal genes, as a considerable percentage of mutations alter splicing or the relative production of alternative transcripts. Modulation of alternative splicing in the retina is also instrumental in the design of novel therapeutic approaches for retinal dystrophies, since it enables precision medicine for specific mutations.

NineTeen Complex-subunit Salsa is required for efficient splicing of a subset of introns and dorsal-ventral patterning

The NineTeen Complex (NTC), also known as pre-mRNA-processing factor 19 (Prp19) complex, regulates distinct spliceosome conformational changes necessary for splicing. During Drosophila midblastula transition, splicing is particularly sensitive to mutations in NTC-subunit Fandango, which suggests differential requirements of NTC during development. We show that NTC-subunit Salsa, the Drosophila ortholog of human RNA helicase Aquarius, is rate-limiting for splicing of a subset of small first introns during oogenesis, including the first intron of gurken Germline depletion of Salsa and splice site mutations within gurken first intron impair both adult female fertility and oocyte dorsal-ventral patterning, due to an abnormal expression of Gurken. Supporting causality, the fertility and dorsal-ventral patterning defects observed after Salsa depletion could be suppressed by the expression of a gurken construct without its first intron. Altogether, our results suggest that one of the key rate-limiting functions of Salsa during oogenesis is to ensure the correct expression and efficient splicing of the first intron of gurken mRNA. Retention of gurken first intron compromises the function of this gene most likely because it undermines the correct structure and function of the transcript 5'UTR.

The whale shark genome reveals how genomic and physiological properties scale with body size

The endangered whale shark (Rhincodon typus) is the largest fish on Earth and a long-lived member of the ancient Elasmobranchii clade. To characterize the relationship between genome features and biological traits, we sequenced and assembled the genome of the whale shark and compared its genomic and physiological features to those of 83 animals and yeast. We examined the scaling relationships between body size, temperature, metabolic rates, and genomic features and found both general correlations across the animal kingdom and features specific to the whale shark genome. Among animals, increased lifespan is positively correlated to body size and metabolic rate. Several genomic traits also significantly correlated with body size, including intron and gene length. Our large-scale comparative genomic analysis uncovered general features of metazoan genome architecture: Guanine and cytosine (GC) content and codon adaptation index are negatively correlated, and neural connectivity genes are longer than average genes in most genomes. Focusing on the whale shark genome, we identified multiple features that significantly correlate with lifespan. Among these were very long gene length, due to introns being highly enriched in repetitive elements such as CR1-like long interspersed nuclear elements, and considerably longer neural genes of several types, including connectivity, activity, and neurodegeneration genes. The whale shark genome also has the second slowest evolutionary rate observed in vertebrates to date. Our comparative genomics approach uncovered multiple genetic features associated with body size, metabolic rate, and lifespan and showed that the whale shark is a promising model for studies of neural architecture and lifespan.

Exon size and sequence conservation improves identification of splice-altering nucleotides

Pre-mRNA splicing is regulated through multiple trans-acting splicing factors. These regulators interact with the pre-mRNA at intronic and exonic positions. Given that most exons are protein coding, the evolution of exons must be modulated by a combination of selective coding and splicing pressures. It has previously been demonstrated that selective splicing pressures are more easily deconvoluted when phylogenetic comparisons are made for exons of identical size, suggesting that exon size-filtered sequence alignments may improve identification of nucleotides evolved to mediate efficient exon ligation. To test this hypothesis, an exon size database was created, filtering 76 vertebrate sequence alignments based on exon size conservation. In addition to other genomic parameters, such as splice-site strength, gene position, or flanking intron length, this database permits the identification of exons that are size- and/or sequence-conserved. Highly size-conserved exons are always sequence-conserved. However, sequence conservation does not necessitate exon size conservation. Our analysis identified evolutionarily young exons and demonstrated that length conservation is a strong predictor of alternative splicing. A published data set of approximately 5000 exonic SNPs associated with disease was analyzed to test the hypothesis that exon size-filtered sequence comparisons increase detection of splice-altering nucleotides. Improved splice predictions could be achieved when mutations occur at the third codon position, especially when a mutation decreases exon inclusion efficiency. The results demonstrate that coding pressures dominate nucleotide composition at invariable codon positions and that exon size-filtered sequence alignments permit identification of splice-altering nucleotides at wobble positions.

Combinatorial regulation of alternative splicing

The generation of protein coding mRNAs from pre-mRNA is a fundamental biological process that is required for gene expression. Alternative pre-mRNA splicing is responsible for much of the transcriptomic and proteomic diversity observed in higher order eukaryotes. Aberrations that disrupt regular alternative splicing patterns are known to cause human diseases, including various cancers. Alternative splicing is a combinatorial process, meaning many factors affect which two splice sites are ligated together. The features that dictate exon inclusion are comprised of splice site strength, intron-exon architecture, RNA secondary structure, splicing regulatory elements, promoter use and transcription speed by RNA polymerase and the presence of post-transcriptional nucleotide modifications. A comprehensive view of all of the factors that influence alternative splicing decisions is necessary to predict splicing outcomes and to understand the molecular basis of disease. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Splicing of long non-coding RNAs primarily depends on polypyrimidine tract and 5' splice-site sequences due to weak interactions with SR proteins

Many nascent long non-coding RNAs (lncRNAs) undergo the same maturation steps as pre-mRNAs of protein-coding genes (PCGs), but they are often poorly spliced. To identify the underlying mechanisms for this phenomenon, we searched for putative splicing inhibitory sequences using the ncRNA-a2 as a model. Genome-wide analyses of intergenic lncRNAs (lincRNAs) revealed that lincRNA splicing efficiency positively correlates with 5'ss strength while no such correlation was identified for PCGs. In addition, efficiently spliced lincRNAs have higher thymidine content in the polypyrimidine tract (PPT) compared to efficiently spliced PCGs. Using model lincRNAs, we provide experimental evidence that strengthening the 5'ss and increasing the T content in PPT significantly enhances lincRNA splicing. We further showed that lincRNA exons contain less putative binding sites for SR proteins. To map binding of SR proteins to lincRNAs, we performed iCLIP with SRSF2, SRSF5 and SRSF6 and analyzed eCLIP data for SRSF1, SRSF7 and SRSF9. All examined SR proteins bind lincRNA exons to a much lower extent than expression-matched PCGs. We propose that lincRNAs lack the cooperative interaction network that enhances splicing, which renders their splicing outcome more dependent on the optimality of splice sites.

The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture

Production of most eukaryotic mRNAs requires splicing of introns from pre-mRNA. The splicing reaction requires definition of splice sites, which are initially recognized in either intron-spanning (‘intron definition’) or exon-spanning (‘exon definition’) pairs. To understand how exon and intron length and splice site recognition mode impact splicing, we measured splicing rates genome-wide in Drosophila, using metabolic labeling/RNA sequencing and new mathematical models to estimate rates. We found that the modal intron length range of 60–70 nt represents a local maximum of splicing rates, but that much longer exon-defined introns are spliced even faster and more accurately. We observed unexpectedly low variation in splicing rates across introns in the same gene, suggesting the presence of gene-level influences, and we identified multiple gene level variables associated with splicing rate. Together our data suggest that developmental and stress response genes may have preferentially evolved exon definition in order to enhance the rate or accuracy of splicing.

Splicing and transcription touch base: co-transcriptional spliceosome assembly and function

Several macromolecular machines collaborate to produce eukaryotic messenger RNA. RNA polymerase II (Pol II) translocates along genes that are up to millions of base pairs in length and generates a flexible RNA copy of the DNA template. This nascent RNA harbours introns that are removed by the spliceosome, which is a megadalton ribonucleoprotein complex that positions the distant ends of the intron into its catalytic centre. Emerging evidence that the catalytic spliceosome is physically close to Pol II in vivo implies that transcription and splicing occur on similar timescales and that the transcription and splicing machineries may be spatially constrained. In this Review, we discuss aspects of spliceosome assembly, transcription elongation and other co-transcriptional events that allow the temporal coordination of co-transcriptional splicing.

Estimating the prevalence of functional exonic splice regulatory information

In addition to coding information, human exons contain sequences necessary for correct splicing. These elements are known to be under purifying selection and their disruption can cause disease. However, the density of functional exonic splicing information remains profoundly uncertain. Several groups have experimentally investigated how mutations at different exonic positions affect splicing. They have found splice information to be distributed widely in exons, with one estimate putting the proportion of splicing-relevant nucleotides at >90%. These results suggest that splicing could place a major pressure on exon evolution. However, analyses of sequence conservation have concluded that the need to preserve splice regulatory signals only slightly constrains exon evolution, with a resulting decrease in the average human rate of synonymous evolution of only 1–4%. Why do these two lines of research come to such different conclusions? Among other reasons, we suggest that the methods are measuring different things: one assays the density of sites that affect splicing, the other the density of sites whose effects on splicing are visible to selection. In addition, the experimental methods typically consider short exons, thereby enriching for nucleotides close to the splice junction, such sites being enriched for splice-control elements. By contrast, in part owing to correction for nucleotide composition biases and to the assumption that constraint only operates on exon ends, the conservation-based methods can be overly conservative.

Microexons: discovery, regulation, and function

The importance of RNA splicing in numerous cellular processes is well established. However, an underappreciated aspect is the ability of the spliceosome to recognize a set of very small (3-30 nucleotide, 1-10 amino acid) exons named microexons. Despite their small size, microexons and their regulation through alternative splicing have now been shown to play critical roles in protein and system function. Here we review the discovery of microexons over time and the mechanisms by which their splicing is regulated, including recent progress made through deep RNA sequencing. We also discuss the functional role of microexons in biology and disease. WIREs RNA 2017, 8:e1418. doi: 10.1002/wrna.1418 For further resources related to this article, please visit the WIREs website.

Selection preserves Ubiquitin Specific Protease 4 alternative exon skipping in therian mammals

Ubiquitin specific protease 4 (USP4) is a highly networked deubiquitinating enzyme with reported roles in cancer, innate immunity and RNA splicing. In mammals it has two dominant isoforms arising from inclusion or skipping of exon 7 (E₇). We evaluated two plausible mechanisms for the generation of these isoforms: (A) E₇ skipping due to a long upstream intron and (B) E₇ skipping due to inefficient 5′ splice sites (5′SS) and/or branchpoint sites (BPS). We then assessed whether E₇ alternative splicing is maintained by selective pressure or arose from genetic drift. Both transcript variants were generated from a USP4-E₇ minigene construct with short flanking introns, an observation consistent with the second mechanism whereby differential splice signal strengths are the basis of E₇ skipping. Optimization of the downstream 5′SS eliminated E₇ skipping. Experimental validation of the correlation between 5′SS identity and exon skipping in vertebrates pinpointed the +6 site as the key splicing determinant. Therian mammals invariably display a 5′SS configuration favouring alternative splicing and the resulting isoforms have distinct subcellular localizations. We conclude that alternative splicing of mammalian USP4 is under selective maintenance and that long and short USP4 isoforms may target substrates in various cellular compartments.

Vemurafenib-resistant BRAF selects alternative branch points different from its wild-type BRAF in intron 8 for RNA splicing

One mechanism of resistance of the melanoma-associated BRAF kinase to its small molecule inhibitor vemurafenib is by point mutations in its intron 8 resulting in exons 4–8 skipping. In this report, we carried out in vitro BRAF RNA splicing assays and lariat RT-PCR to map the intron 8 branch points in wild-type and BRAF mutants. We identify multiple branch points (BP) in intron 8 of both wild-type (wt) and vemurafenib-resistant BRAF RNA. In wt BRAF, BPs are located at -29A, -28A and -26A, whereas in a vemurafenib-resistant BRAF splicing mutant, BPs map to -22A, -18A and -15A, proximal to the intron 8 3′ splice site. This finding of a distal-to-proximal shift of the branch point sequence in BRAF splicing in response to point-mutations in intron 8 provides insight into the regulation of BRAF alternative splicing upon vemurafenib resistance.

Alternative RNA splicing of KSHV ORF57 produces two different RNA isoforms

In lytically infected B cells Kaposi sarcoma-associated herpesvirus (KSHV) ORF57 gene encodes two RNA isoforms by alternative splicing of its pre-mRNA, which contains a small, constitutive intron in its 5' half and a large, suboptimal intron in its 3's half. The RNA1 isoform encodes full-length ORF57 and is a major isoform derived from splicing of the constitutive small intron, but retaining the suboptimal large intron as the coding region. A small fraction (<5%) of ORF57 RNA undergoes double splicing to produce a smaller non-coding RNA2 due to lack of a translational termination codon. Both RNAs are cleaved and polyadenylated at the same cleavage site CS83636. The insertion of ORF57 RNA1 into a restriction cutting site in certain mammalian expression vectors activates splicing of the subopitmal intron and produces a truncated ORF57 protein.

The Orthology Clause in the Next Generation Sequencing Era: Novel Reference Genes Identified by RNA-seq in Humans Improve Normalization of Neonatal Equine Ovary RT-qPCR Data

BACKGROUND

Vertebrate evolution is accompanied by a substantial conservation of transcriptional programs with more than a third of unique orthologous genes showing constrained levels of expression. Moreover, there are genes and exons exhibiting excellent expression stability according to RNA-seq data across a panel of eighteen tissues including the ovary (Human Body Map 2.0).

RESULTS

We hypothesized that orthologs of these exons would also be highly uniformly expressed across neonatal ovaries of the horse, which would render them appropriate reference genes (RGs) for normalization of reverse transcription quantitative PCR (RT-qPCR) data in this context. The expression stability of eleven novel RGs (C1orf43, CHMP2A, EMC7, GPI, PSMB2, PSMB4, RAB7A, REEP5, SNRPD3, VCP and VPS29) was assessed by RT-qPCR in ovaries of seven neonatal fillies and compared to that of the expressed repetitive element ERE-B, two universal (OAZ1 and RPS29) and four traditional RGs (ACTB, GAPDH, UBB and B2M). Expression stability analyzed with the software tool RefFinder top ranked the normalization factor constituted of the genes SNRPD3 and VCP, a gene pair that is not co-expressed according to COEXPRESdb and GeneMANIA. The traditional RGs GAPDH, B2M, ACTB and UBB were only ranked 3rd and 12th to 14th, respectively.

CONCLUSIONS

The functional diversity of the novel RGs likely facilitates expression studies over a wide range of physiological and pathological contexts related to the neonatal equine ovary. In addition, this study augments the potential for RT-qPCR-based profiling of human samples by introducing seven new human RG assays (C1orf43, CHMP2A, EMC7, GPI, RAB7A, VPS29 and UBB).

Prp40 and early events in splice site definition

The alternative splicing (AS) of precursor messenger RNA (pre-mRNA) is a tightly regulated process through which introns are removed to leave the resulting exons in the mRNA appropriately aligned and ligated. The AS of pre-mRNA is a key mechanism for increasing the complexity of proteins encoded in the genome. In humans, more than 90% of genes undergo AS, underscoring the importance of this process in RNA biogenesis. As such, AS misregulation underlies multiple human diseases. The splicing reaction is catalyzed by the spliceosome, a highly dynamic complex that assembles at or near the intron/exon boundaries and undergoes sequential conformational and compositional changes during splicing. The initial recognition of splice sites defines the exons that are going to be removed, which is a critical step in the highly regulated splicing process. Although the available lines of evidence are increasing, the molecular mechanisms governing AS, including the initial interactions occurring at intron/exon boundaries, and the factors that modulate these critical connections by functioning as a scaffold for active-site RNAs or proteins, remain poorly understood. In this review, we summarize the major hallmarks of the initial steps in the splicing process and the role of auxiliary factors that contribute to the assembly of the spliceosomal complex. We also discuss the role of the essential yeast Prp40 protein and its mammalian homologs in the specificity of this pre-mRNA processing event. In addition, we provide the first exhaustive phylogenetic analysis of the molecular evolution of Prp40 family members. WIREs RNA 2016, 7:17-32. doi: 10.1002/wrna.1312 For further resources related to this article, please visit the WIREs website.

A diverse epigenetic landscape at human exons with implication for expression

DNA methylation is an important epigenetic marker associated with gene expression regulation in eukaryotes. While promoter methylation is relatively well characterized, the role of intragenic DNA methylation remains unclear. Here, we investigated the relationship of DNA methylation at exons and flanking introns with gene expression and histone modifications generated from a human fibroblast cell-line and primary B cells. Consistent with previous work we found that intragenic methylation is positively correlated with gene expression and that exons are more highly methylated than their neighboring intronic environment. Intriguingly, in this study we identified a unique subset of hypomethylated exons that demonstrate significantly lower methylation levels than their surrounding introns. Furthermore, we observed a negative correlation between exon methylation and the density of the majority of histone modifications. Specifically, we demonstrate that hypo-methylated exons at highly expressed genes are associated with open chromatin and have a characteristic histone code comprised of significantly high levels of histone markings. Overall, our comprehensive analysis of the human exome supports the presence of regulatory hypomethylated exons in protein coding genes. In particular our results reveal a previously unrecognized diverse and complex role of the epigenetic landscape within the gene body.

Splicing of designer exons informs a biophysical model for exon definition

Pre-mRNA molecules in humans contain mostly short internal exons flanked by longer introns. To explain the removal of such introns, exon recognition instead of intron recognition has been proposed. We studied this exon definition using designer exons (DEs) made up of three prototype modules of our own design: an exonic splicing enhancer (ESE), an exonic splicing silencer (ESS), and a Reference Sequence (R) predicted to be neither. Each DE was examined as the central exon in a three-exon minigene. DEs made of R modules showed a sharp size dependence, with exons shorter than 14 nt and longer than 174 nt splicing poorly. Changing the strengths of the splice sites improved longer exon splicing but worsened shorter exon splicing, effectively displacing the curve to the right. For the ESE we found, unexpectedly, that its enhancement efficiency was independent of its position within the exon. For the ESS we found a step-wise positional increase in its effects; it was most effective at the 3' end of the exon. To apply these results quantitatively, we developed a biophysical model for exon definition of internal exons undergoing cotranscriptional splicing. This model features commitment to inclusion before the downstream exon is synthesized and competition between skipping and inclusion fates afterward. Collision of both exon ends to form an exon definition complex was incorporated to account for the effect of size; ESE/ESS effects were modeled on the basis of stabilization/destabilization. This model accurately predicted the outcome of independent experiments on more complex DEs that combined ESEs and ESSs.

RBFOX and PTBP1 proteins regulate the alternative splicing of micro-exons in human brain transcripts

Ninety-four percent of mammalian protein-coding exons exceed 51 nucleotides (nt) in length. The paucity of micro-exons (≤ 51 nt) suggests that their recognition and correct processing by the splicing machinery present greater challenges than for longer exons. Yet, because thousands of human genes harbor processed micro-exons, specialized mechanisms may be in place to promote their splicing. Here, we survey deep genomic data sets to define 13,085 micro-exons and to study their splicing mechanisms and molecular functions. More than 60% of annotated human micro-exons exhibit a high level of sequence conservation, an indicator of functionality. While most human micro-exons require splicing-enhancing genomic features to be processed, the splicing of hundreds of micro-exons is enhanced by the adjacent binding of splice factors in the introns of pre-messenger RNAs. Notably, splicing of a significant number of micro-exons was found to be facilitated by the binding of RBFOX proteins, which promote their inclusion in the brain, muscle, and heart. Our analyses suggest that accurate regulation of micro-exon inclusion by RBFOX proteins and PTBP1 plays an important role in the maintenance of tissue-specific protein-protein interactions.

Pre-mRNA splicing in cancer: the relevance in oncogenesis, treatment and drug resistance

INTRODUCTION

Aberrant pre-mRNA splicing in cancer is emerging as an important determinant of oncogenesis, response to treatment and anticancer drug resistance. At the same time, the spliceosome has become a target for a novel class of pre-clinical chemotherapeutics with a potential future application in cancer treatment. Taken together, these findings offer novel opportunities for the enhancement of the efficacy of cancer therapy.

AREAS COVERED

This review presents a comprehensive overview of the molecular mechanisms involved in splicing and current developments regarding splicing aberrations in relation to several aspects of cancer formation and therapy. Identified mutations in the various components of the spliceosome and their implications for cancer prognosis are delineated. Moreover, the contribution of abnormal splicing patterns as well as deregulated splicing factors to chemoresistance is discussed, along with novel splicing-based therapeutic approaches.

EXPERT OPINION

Significant progress has been made in deciphering the role of splicing factors in cancer including carcinogenesis and drug resistance. Splicing-based prognostic tools as well as therapeutic options hold great potential towards improvements in cancer therapy. However, gaining more in-depth molecular insight into the consequences of mutations in various components of the splicing machinery as well as of cellular effects of spliceosome inhibition is a prerequisite to establish the role of splicing in tumor progression and treatment options, respectively.

Genome-wide survey of Alternative Splicing in Sorghum Bicolor

Sorghum bicolor is a member of grass family which is an attractive model plant for genome study due to interesting genome features like low genome size. In this research, we performed comprehensive investigation of Alternative Splicing and ontology aspects of genes those have undergone these events in sorghum bicolor. We used homology based alignments between gene rich transcripts, represented by tentative consensus (TC) transcript sequences, and genomic scaffolds to deduce the structure of genes and identify alternatively spliced transcripts in sorghum. Using homology mapping of assembled expressed sequence tags with genomics data, we identified 2,137 Alternative Splicing events in S. bicolor. Our study showed that complex events and intron retention are the main types of Alternative Splicing events in S. bicolor and highlights the prevalence of splicing site recognition for definition of introns in this plant. Annotations of the alternatively spliced genes revealed that they represent diverse biological process and molecular functions, suggesting a fundamental role for Alternative Splicing in affecting the development and physiology of S. bicolor.

A mutation in a splicing factor that causes retinitis pigmentosa has a transcriptome-wide effect on mRNA splicing

BACKGROUND

Substantial progress has been made in the identification of sequence elements that control mRNA splicing and the genetic variants in these elements that alter mRNA splicing (referred to as splicing quantitative trait loci - sQTLs). Genetic variants that affect mRNA splicing in trans are harder to identify because their effects can be more subtle and diffuse, and the variants are not co-located with their targets. We carried out a transcriptome-wide analysis of the effects of a mutation in a ubiquitous splicing factor that causes retinitis pigmentosa (RP) on mRNA splicing, using exon microarrays.

RESULTS

Exon microarray data was generated from whole blood samples obtained from four individuals with a mutation in the splicing factor PRPF8 and four sibling controls. Although the mutation has no known phenotype in blood, there was evidence of widespread differences in splicing between cases and controls (affecting approximately 20% of exons). Most probesets with significantly different inclusion (defined as the expression intensity of the exon divided by the expression of the corresponding transcript) between cases and controls had higher inclusion in cases and corresponded to exons that were shorter than average, AT rich, located towards the 5' end of the gene and flanked by long introns. Introns flanking affected probesets were particularly depleted for the shortest category of introns, associated with splicing via intron definition.

CONCLUSIONS

Our results show that a mutation in a splicing factor, with a phenotype that is restricted to retinal tissue, acts as a trans-sQTL cluster in whole blood samples. Characteristics of the affected exons suggest that they are spliced co-transcriptionally and via exon definition. However, due to the small sample size available for this study, further studies are required to confirm the widespread impact of this PRPF8 mutation on mRNA splicing outside the retina.

A day in the life of the spliceosome

One of the most amazing findings in molecular biology was the discovery that eukaryotic genes are discontinuous, with coding DNA being interrupted by stretches of non-coding sequence. The subsequent realization that the intervening regions are removed from pre-mRNA transcripts via the activity of a common set of small nuclear RNAs (snRNAs), which assemble together with associated proteins into a complex known as the spliceosome, was equally surprising. How do cells coordinate the assembly of this molecular machine? And how does the spliceosome accurately recognize exons and introns to carry out the splicing reaction? Insights into these questions have been gained by studying the life cycle of spliceosomal snRNAs from their transcription, nuclear export and re-import to their dynamic assembly into the spliceosome. This assembly process can also affect the regulation of alternative splicing and has implications for human disease.

Global splicing pattern reversion during somatic cell reprogramming

Alternative splicing generates multiple transcripts from a single gene, and cell-type-specific splicing profiles are important for the properties and functions of the cells. Recently, somatic cells have been shown to undergo dedifferentiation after the forced expression of transcription factors. However, it remains unclear whether somatic cell splicing is reorganized during reprogramming. Here, by combining deep sequencing with high-throughput absolute qRT-PCR, we show that somatic splicing profiles revert to pluripotent ones during reprogramming. Remarkably, the splicing pattern in pluripotent stem cells resembles that in testes, and the regulatory regions have specific characteristics in length and sequence. Furthermore, our siRNA screen has identified RNA-binding proteins that regulate splicing events in iPSCs. We have then demonstrated that two of the RNA-binding proteins, U2af1 and Srsf3, play a role in somatic cell reprogramming. Our results indicate that the drastic alteration in splicing represents part of the molecular network involved in the reprogramming process.

A novel otoferlin splice-site mutation in siblings with auditory neuropathy spectrum disorder

We characterize a novel otoferlin mutation discovered in a sibling pair diagnosed with auditory neuropathy spectrum disorder and investigate auditory nerve function through their cochlear implants. Genetic sequencing revealed a homozygous mutation at the otoferlin splice donor site of exon 28 (IVS28 + 1G>T) in both siblings. Functional investigation showed that the intronic sequence between exons 28 and 29 was retained in the mutated minigenes that were expressed in 293T cells. Auditory nerve compound action potential recovery functions in the siblings demonstrated different rates of neural recovery, with sibling AN1 showing rapid recovery (1.14 ms) and AN2 showing average recovery (0.78 ms) compared to subjects with sensorineural hearing loss (average: adults 0.71 ms, children 0.85 ms). Differences in neural recovery were consistent with speech perception differences between the siblings. Genotype information may indicate site of lesion in hearing loss; however, additional, as yet, unknown factors may impact clinical outcomes and must be considered.

A simple model to explain evolutionary trends of eukaryotic gene architecture and expression: how competition between splicing and cleavage/polyadenylation factors may affect gene expression and splice-site recognition in eukaryotes

Enormous phylogenetic variation exists in the number and sizes of introns in protein-coding genes. Although some consideration has been given to the underlying role of the population-genetic environment in defining such patterns, the influence of the intracellular environment remains virtually unexplored. Drawing from observations on interactions between co-transcriptional processes involved in splicing and mRNA 3'-end formation, a mechanistic model is proposed for splice-site recognition that challenges the commonly accepted intron- and exon-definition models. Under the suggested model, splicing factors that outcompete 3'-end processing factors for access to intronic binding sites concurrently favor the recruitment of 3'-end processing factors at the pre-mRNA tail. This hypothesis sheds new light on observations such as the intron-mediated enhancement of gene expression and the negative correlation between intron length and levels of gene expression.

Testing for natural selection in human exonic splicing regulators associated with evolutionary rate shifts

Despite evidence that at the interspecific scale, exonic splicing silencers (ESSs) are under negative selection in constitutive exons, little is known about the effects of slightly deleterious polymorphisms on these splicing regulators. Through the application of a modified version of the McDonald-Kreitman test, we compared the normalized proportions of human polymorphisms and human/rhesus substitutions affecting exonic splicing regulators (ESRs) on sequences of constitutive and alternative exons. Our results show a depletion of substitutions and an enrichment of SNPs associated with ESS gain in constitutive exons. Moreover, we show that this evolutionary pattern is also present in a set of ESRs previously involved in the transition from constitutive to skipped exons in the mammalian lineage. The similarity between these two sets of ESRs suggests that the transition from constitutive to skipped exons in mammals is more frequently associated with the inhibition than with the promotion of splicing signals. This is in accordance with the hypothesis of a constitutive origin of exon skipping and corroborates previous findings about the antagonistic role of certain exonic splicing enhancers.

Interpretation, stratification and evidence for sequence variants affecting mRNA splicing in complete human genome sequences

Information theory-based methods have been shown to be sensitive and specific for predicting and quantifying the effects of non-coding mutations in Mendelian diseases. We present the Shannon pipeline software for genome-scale mutation analysis and provide evidence that the software predicts variants affecting mRNA splicing. Individual information contents (in bits) of reference and variant splice sites are compared and significant differences are annotated and prioritized. The software has been implemented for CLC-Bio Genomics platform. Annotation indicates the context of novel mutations as well as common and rare SNPs with splicing effects. Potential natural and cryptic mRNA splicing variants are identified, and null mutations are distinguished from leaky mutations. Mutations and rare SNPs were predicted in genomes of three cancer cell lines (U2OS, U251 and A431), which were supported by expression analyses. After filtering, tractable numbers of potentially deleterious variants are predicted by the software, suitable for further laboratory investigation. In these cell lines, novel functional variants comprised 6–17 inactivating mutations, 1–5 leaky mutations and 6–13 cryptic splicing mutations. Predicted effects were validated by RNA-seq analysis of the three aforementioned cancer cell lines, and expression microarray analysis of SNPs in HapMap cell lines.

Cotranscriptional splicing efficiency differs dramatically between Drosophila and mouse

Spliceosome assembly and/or splicing of a nascent transcript may be crucial for proper isoform expression and gene regulation in higher eukaryotes. We recently showed that cotranscriptional splicing occurs efficiently in Drosophila, but there are not comparable genome-wide nascent splicing data from mammals. To provide this comparison, we analyze a recently generated, high-throughput sequencing data set of mouse liver nascent RNA, originally studied for circadian transcriptional regulation. Cotranscriptional splicing is approximately twofold less efficient in mouse liver than in Drosophila, i.e., nascent intron levels relative to exon levels are ∼0.55 in mouse versus 0.25 in the fly. An additional difference between species is that only mouse cotranscriptional splicing is optimal when 5'-exon length is between 50 and 500 bp, and intron length does not correlate with splicing efficiency, consistent with exon definition. A similar analysis of intron and exon length dependence in the fly is more consistent with intron definition. Contrasted with these differences are many similarities between the two systems: Alternatively annotated introns are less efficiently spliced cotranscriptionally than constitutive introns, and introns of single-intron genes are less efficiently spliced than introns from multi-intron genes. The most striking common feature is intron position: Cotranscriptional splicing is much more efficient when introns are far from the 3' ends of their genes. Additionally, absolute gene length correlates positively with cotranscriptional splicing efficiency independently of intron location and position, in flies as well as in mice. The gene length and distance effects indicate that more "nascent time" gives rise to greater cotranscriptional splicing efficiency in both systems.

Intron definition and a branch site adenosine at nt 385 control RNA splicing of HPV16 E6*I and E7 expression

HPV16 E6 and E7, two viral oncogenes, are expressed from a single bicistronic pre-mRNA. In this report, we provide the evidence that the bicistronic pre-mRNA intron 1 contains three 5' splice sites (5' ss) and three 3' splice sites (3' ss) normally used in HPV16(+) cervical cancer and its derived cell lines. The choice of two novel alternative 5' ss (nt 221 5' ss and nt 191 5' ss) produces two novel isoforms of E6E7 mRNAs (E6*V and E6*VI). The nt 226 5' ss and nt 409 3' ss is preferentially selected over the other splice sites crossing over the intron to excise a minimal length of the intron in RNA splicing. We identified AACAAAC as the preferred branch point sequence (BPS) and an adenosine at nt 385 (underlined) in the BPS as a branch site to dictate the selection of the nt 409 3' ss for E6*I splicing and E7 expression. Introduction of point mutations into the mapped BPS led to reduced U2 binding to the BPS and thereby inhibition of the second step of E6E7 splicing at the nt 409 3' ss. Importantly, the E6E7 bicistronic RNA with a mutant BPS and inefficient splicing makes little or no E7 and the resulted E6 with mutations of (91)QYNK(94) to (91)PSFW(94) displays attenuate activity on p53 degradation. Together, our data provide structural basis of the E6E7 intron 1 for better understanding of how viral E6 and E7 expression is regulated by alternative RNA splicing. This study elucidates for the first time a mapped branch point in HPV16 genome involved in viral oncogene expression.

Chromatin and epigenetic regulation of pre-mRNA processing

New data are revealing a complex landscape of gene regulation shaped by chromatin states that extend into the bodies of transcribed genes and associate with distinct RNA elements such as exons, introns and polyadenylation sites. Exons are characterized by increased levels of nucleosome positioning, DNA methylation and certain histone modifications. As pre-mRNA splicing occurs co-transcriptionally, changes in the transcription elongation rate or epigenetic marks can influence exon splicing. These new discoveries broaden our understanding of the epigenetic code and ascribe a novel role for chromatin in controlling pre-mRNA processing. In this review, we summarize the recently discovered interplay between the modulation of chromatin states and pre-mRNA processing with the particular focus on how these processes communicate with one another to control gene expression.

Differential GC content between exons and introns establishes distinct strategies of splice-site recognition

During evolution segments of homeothermic genomes underwent a GC content increase. Our analyses reveal that two exon-intron architectures have evolved from an ancestral state of low GC content exons flanked by short introns with a lower GC content. One group underwent a GC content elevation that abolished the differential exon-intron GC content, with introns remaining short. The other group retained the overall low GC content as well as the differential exon-intron GC content, and is associated with longer introns. We show that differential exon-intron GC content regulates exon inclusion level in this group, in which disease-associated mutations often lead to exon skipping. This group's exons also display higher nucleosome occupancy compared to flanking introns and exons of the other group, thus "marking" them for spliceosomal recognition. Collectively, our results reveal that differential exon-intron GC content is a previously unidentified determinant of exon selection and argue that the two GC content architectures reflect the two mechanisms by which splicing signals are recognized: exon definition and intron definition.

Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq

An extensive repertoire of modifications is known to underlie the versatile coding, structural and catalytic functions of RNA, but it remains largely uncharted territory. Although biochemical studies indicate that N(6)-methyladenosine (m(6)A) is the most prevalent internal modification in messenger RNA, an in-depth study of its distribution and functions has been impeded by a lack of robust analytical methods. Here we present the human and mouse m(6)A modification landscape in a transcriptome-wide manner, using a novel approach, m(6)A-seq, based on antibody-mediated capture and massively parallel sequencing. We identify over 12,000 m(6)A sites characterized by a typical consensus in the transcripts of more than 7,000 human genes. Sites preferentially appear in two distinct landmarks--around stop codons and within long internal exons--and are highly conserved between human and mouse. Although most sites are well preserved across normal and cancerous tissues and in response to various stimuli, a subset of stimulus-dependent, dynamically modulated sites is identified. Silencing the m(6)A methyltransferase significantly affects gene expression and alternative splicing patterns, resulting in modulation of the p53 (also known as TP53) signalling pathway and apoptosis. Our findings therefore suggest that RNA decoration by m(6)A has a fundamental role in regulation of gene expression.

Immunoglobulin heavy chain gene regulation through polyadenylation and splicing competition

The immunoglobulin heavy chain (IgH) genes, which encode one of the two chains of antibody molecules, were the first cellular genes shown to undergo developmentally regulated alternative RNA processing. These genes produce two different mRNAs from a single primary transcript. One mRNA is cleaved and polyadenylated at an upstream poly(A) signal while the other mRNA removes this poly(A) signal by RNA splicing and is cleaved and polyadenylated at a downstream poly(A) site. A broad range of studies have been performed to understand the mechanism of IgH RNA processing regulation during B lymphocyte development. The model that has emerged is much more complex than envisioned by the earliest view of regulation through poly(A) signal choice. Regulation requires that the IgH gene contain competing splice and cleavage-polyadenylation reactions with balanced efficiencies. Because non-IgH genes with these structural features also can be regulated, IgH gene-specific sequence elements are not required for regulation. Changes in cleavage-polyadenylation and RNA splicing, as well as pol II elongation, all contribute to IgH developmental RNA processing regulation. Multiple factors are likely involved in the regulation during B lymphocyte maturation. Additional biologically relevant factors that contribute to IgH regulation remain to be identified and incorporated into a mechanistic model for regulation. Much of the work to date confirms the complex nature of IgH mRNA regulation and suggests that a thorough understanding of this control will remain a challenge. However, it is also likely that such understanding will help elucidate novel mechanisms of RNA processing regulation.

Efficient internal exon recognition depends on near equal contributions from the 3' and 5' splice sites

Pre-mRNA splicing is carried out by the spliceosome, which identifies exons and removes intervening introns. In vertebrates, most splice sites are initially recognized by the spliceosome across the exon, because most exons are small and surrounded by large introns. This gene architecture predicts that efficient exon recognition depends largely on the strength of the flanking 3' and 5' splice sites. However, it is unknown if the 3' or the 5' splice site dominates the exon recognition process. Here, we test the 3' and 5' splice site contributions towards efficient exon recognition by systematically replacing the splice sites of an internal exon with sequences of different splice site strengths. We show that the presence of an optimal splice site does not guarantee exon inclusion and that the best predictor for exon recognition is the sum of both splice site scores. Using a genome-wide approach, we demonstrate that the combined 3' and 5' splice site strengths of internal exons provide a much more significant separator between constitutive and alternative exons than either the 3' or the 5' splice site strength alone.