Four exons of the serotonin receptor 4 gene are associated with multiple distant branch points

Restriction of an intron size en route to endothermy

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5' splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5' splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

Inducible nuclear import by TetR aptamer-controlled 3' splice site selection

Correct cellular localization is essential for the function of many eukaryotic proteins and hence cell physiology. Here, we present a synthetic genetic device that allows the control of nuclear and cytosolic localization based on controlled alternative splicing in human cells. The device is based on the fact that an alternative 3' splice site is located within a TetR aptamer that in turn is positioned between the branch point and the canonical splice site. The novel splice site is only recognized when the TetR repressor is bound. Addition of doxycycline prevents TetR aptamer binding and leads to recognition of the canonical 3' splice site. It is thus possible to produce two independent splice isoforms. Since the terminal loop of the aptamer may be replaced with any sequence of choice, one of the two isoforms may be extended by the respective sequence of choice depending on the presence of doxycycline. In a proof-of-concept study, we fused a nuclear localization sequence to a cytosolic target protein, thus directing the protein into the nucleus. However, the system is not limited to the control of nuclear localization. In principle, any target sequence can be integrated into the aptamer, allowing not only the production of a variety of different isoforms on demand, but also to study the function of mislocalized proteins. Moreover, it also provides a valuable tool for investigating the mechanism of alternative splicing in human cells.

Most human introns are recognized via multiple and tissue-specific branchpoints

Pineda and Bradley demonstrate that almost all human introns contain multiple branchpoints. Approximately three-quarters of constitutive introns exhibit tissue-specific branchpoint usage.

Although branchpoint recognition is an essential component of intron excision during the RNA splicing process, the branchpoint itself is frequently assumed to be a basal, rather than regulatory, sequence feature. However, this assumption has not been systematically tested due to the technical difficulty of identifying branchpoints and quantifying their usage. Here, we analyzed ∼1.31 trillion reads from 17,164 RNA sequencing data sets to demonstrate that almost all human introns contain multiple branchpoints. This complexity holds even for constitutive introns, 95% of which contain multiple branchpoints, with an estimated five to six branchpoints per intron. Introns upstream of the highly regulated ultraconserved poison exons of SR genes contain twice as many branchpoints as the genomic average. Approximately three-quarters of constitutive introns exhibit tissue-specific branchpoint usage. In an extreme example, we observed a complete switch in branchpoint usage in the well-studied first intron of HBB (β-globin) in normal bone marrow versus metastatic prostate cancer samples. Our results indicate that the recognition of most introns is unexpectedly complex and tissue-specific and suggest that alternative splicing catalysis typifies the majority of introns even in the absence of differences in the mature mRNA.

CK2 regulates 5-HT4 receptor signaling and modulates depressive-like behavior

The serotonergic neurotransmitter system has been widely implicated in the pathophysiology of mood-related disorders such as anxiety and major depressive disorder (MDD). The onset of therapeutic efficacy of traditional antidepressants is delayed by several weeks. The 5-HT₄ receptor has emerged as a new therapeutic target since agonists of this receptor induce rapid antidepressant-like responses in rodents. Here we show that the 5-HT₄ receptor is regulated by CK2, at transcriptional and post-transcriptional levels. We present evidence, in two different CK2α knockout mouse lines, that this regulation is region-specific, with the 5-HT₄ receptor upregulated in prefrontal cortex (PFC) but not striatum or hippocampus where CK2α is also ablated. 5-HT₄ receptor signaling is enhanced in vitro, as evidenced by enhanced cAMP production or receptor plasma membrane localization in the presence of CK2 inhibitor or shRNA targeting CK2α. In vivo, 5-HT₄ receptor signaling is also upregulated since ERK activation is elevated and sensitive to the inverse agonist, GR113808 in the PFC of CK2α KO mice. Behaviorally, KO mice as well as mice with AAV-mediated deletion of CK2α in the PFC show a robust 'anti-depressed-like' phenotype and display an enhanced response to antidepressant treatment when tested in paradigms for mood and anxiety. Importantly, it is sufficient to overexpress the 5-HT₄ receptor in the mPFC to generate mice with a similar 'anti-depressed-like' phenotype. Our findings identify the mPFC as the region that mediates the effect of enhanced 5-HT₄ receptor activity and CK2 as modulator of 5-HT4 receptor levels in this brain region that regulates mood-related phenotypes.

Large-scale analysis of branchpoint usage across species and cell lines

The coding sequence of each human pre-mRNA is interrupted, on average, by 11 introns that must be spliced out for proper gene expression. Each intron contains three obligate signals: a 5′ splice site, a branch site, and a 3′ splice site. Splice site usage has been mapped exhaustively across different species, cell types, and cellular states. In contrast, only a small fraction of branch sites have been identified even once. The few reported annotations of branch site are imprecise as reverse transcriptase skips several nucleotides while traversing a 2–5 linkage. Here, we report large-scale mapping of the branchpoints from deep sequencing data in three different species and in the SF3B1 K700E oncogenic mutant background. We have developed a novel method whereby raw lariat reads are refined by U2snRNP/pre-mRNA base-pairing models to return the largest current data set of branchpoint sequences with quality metrics. This analysis discovers novel modes of U2snRNA:pre-mRNA base-pairing conserved in yeast and provides insight into the biogenesis of intron circles. Finally, matching branch site usage with isoform selection across the extensive panel of ENCODE RNA-seq data sets offers insight into the mechanisms by which branchpoint usage drives alternative splicing.

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.

Intronic features that determine the selection of the 3' splice site

Most eukaryotic primary transcripts include segments, or introns, that will be accurately removed during RNA biogenesis. This process, known as pre-messenger RNA splicing, is catalyzed by the spliceosome, accurately selecting a set of intronic marks from others apparently equivalent. This identification is critical, as incorrectly spliced RNAs can be toxic for the organism. One of these marks, the dinucleotide AG, signals the intronic 3' end, or 3' splice site (ss). In this review we will focus on those intronic features that have an impact on 3' ss selection. These include the location and type of neighboring sequences, and their distance to the 3' end. We will see that their interplay is needed to select the right intronic end, and that this can be modulated by additional intronic elements that contribute to alternative splicing, whereby diverse RNAs can be generated from identical precursors. This complexity, still poorly understood, is fundamental for the accuracy of gene expression. In addition, a clear knowledge of 3' ss selection is needed to fully decipher the coding potential of genomes.

Secondary structure is required for 3' splice site recognition in yeast

Higher order RNA structures can mask splicing signals, loop out exons, or constitute riboswitches all of which contributes to the complexity of splicing regulation. We identified a G to A substitution between branch point (BP) and 3′ splice site (3′ss) of Saccharomyces cerevisiae COF1 intron, which dramatically impaired its splicing. RNA structure prediction and in-line probing showed that this mutation disrupted a stem in the BP-3′ss region. Analyses of various COF1 intron modifications revealed that the secondary structure brought about the reduction of BP to 3′ss distance and masked potential 3′ss. We demonstrated the same structural requisite for the splicing of UBC13 intron. Moreover, RNAfold predicted stable structures for almost all distant BP introns in S. cerevisiae and for selected examples in several other Saccharomycotina species. The employment of intramolecular structure to localize 3′ss for the second splicing step suggests the existence of pre-mRNA structure-based mechanism of 3′ss recognition.

Genome-wide association between branch point properties and alternative splicing

The branch point (BP) is one of the three obligatory signals required for pre-mRNA splicing. In mammals, the degeneracy of the motif combined with the lack of a large set of experimentally verified BPs complicates the task of modeling it in silico, and therefore of predicting the location of natural BPs. Consequently, BPs have been disregarded in a considerable fraction of the genome-wide studies on the regulation of splicing in mammals. We present a new computational approach for mammalian BP prediction. Using sequence conservation and positional bias we obtained a set of motifs with good agreement with U2 snRNA binding stability. Using a Support Vector Machine algorithm, we created a model complemented with polypyrimidine tract features, which considerably improves the prediction accuracy over previously published methods. Applying our algorithm to human introns, we show that BP position is highly dependent on the presence of AG dinucleotides in the 3' end of introns, with distance to the 3' splice site and BP strength strongly correlating with alternative splicing. Furthermore, experimental BP mapping for five exons preceded by long AG-dinucleotide exclusion zones revealed that, for a given intron, more than one BP can be chosen throughout the course of splicing. Finally, the comparison between exons of different evolutionary ages and pseudo exons suggests a key role of the BP in the pathway of exon creation in human. Our computational and experimental analyses suggest that BP recognition is more flexible than previously assumed, and it appears highly dependent on the presence of downstream polypyrimidine tracts. The reported association between BP features and the splicing outcome suggests that this, so far disregarded but yet crucial, element buries information that can complement current acceptor site models.