Efficient use of a 'dead-end' GA 5' splice site in the human fibroblast growth factor receptor genes

Splicing analysis of STAT3 tandem donor suggests non-canonical binding registers for U1 and U6 snRNAs

Tandem donor splice sites (5'ss) are unique regions with at least two GU dinucleotides serving as splicing cleavage sites. The Δ3 tandem 5'ss are a specific subclass of 5'ss separated by 3 nucleotides which can affect protein function by inserting/deleting a single amino acid. One 5'ss is typically preferred, yet factors governing particular 5'ss choice are not fully understood. A highly conserved exon 21 of the STAT3 gene was chosen as a model to study Δ3 tandem 5'ss splicing mechanisms. Based on multiple lines of experimental evidence, endogenous U1 snRNA most likely binds only to the upstream 5'ss. However, the downstream 5'ss is used preferentially, and the splice site choice is not dependent on the exact U1 snRNA binding position. Downstream 5'ss usage was sensitive to exact nucleotide composition and dependent on the presence of downstream regulatory region. The downstream 5'ss usage could be best explained by two novel interactions with endogenous U6 snRNA. U6 snRNA enables the downstream 5'ss usage in STAT3 exon 21 by two mechanisms: (i) binding in a novel non-canonical register and (ii) establishing extended Watson-Crick base pairing with the downstream regulatory region. This study suggests that U6:5'ss interaction is more flexible than previously thought.

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.

Novel Variant of OCT4, Named OCT4B5, is Highly Expressed in Human Pluripotent Cells

Alternative promoter and alternative splicing are two important mechanisms of gene regulation and protein diversity in different physiological contexts of eukaryotes, especially in stem cells and developmental stages. Pou5f1 gene which codes the stemness marker OCT4, utilizes alternative splicing and promoter mechanisms, which result in generation of multiple spliced variants and subsequently multiple protein isoforms. By far, nine variants of OCT4 (OCT4A, OCT4B, OCT4B1, OCT4B2, OCT4B3, OCT4B4, OCT4C, OCT4C1, and OCT4D) have been introduced. It has been well established that OCT4A plays essential roles in early developmental stages as well as maintenance of stemness in embryonic stem cells (ESCs). However, the roles and functions of other variants and isoforms of OCT4 in biological systems are less appreciated. In this study, we report a new OCT4 variant, designated as OCT4B5. RT-PCR assay on different human cell lines including pluripotent, normal and cancer cells showed that OCT4B5 is expressed at variable level in different cell lines. By semi-quantifying of OCT4B5 expression in pluripotent and differentiated states of NT2 cell lines, we reveal that this variant of OCT4 is highly expressed in undifferentiated state and its expression is down-regulated upon differentiation. Compared to OCT4A which is sharply down-regulated in retinoic acid induced differentiation of NT2 cell line, the expression of OCT4B5 remains at low level in differentiated state. Overall, this study emphasizes the complexity of OCT4 gene expression and regulation in different states of stem cells and physiological contexts. Graphical Abstract.

Profiling of cis- and trans-acting factors supporting noncanonical splice site activation

Recently, by combining transcriptomics with functional splicing reporter assays we were able to identify GT > GC > TT as the three highest ranked dinucleotides of human 5' splice sites (5'ss). Here, we have extended our investigations to the proteomic characterization of nuclear proteins that bind to canonical and noncanonical 5'ss. Surprisingly, we found that U1 snRNP binding to functional 5'ss sequences prevented components of the DNA damage response (DDR) from binding to the RNA, suggesting a close link between spliceosome arrangement and genome stability. We demonstrate that all tested noncanonical 5'ss sequences are bona-fide targets of the U2-type spliceosome and are bound by U1 snRNP, including U1-C, in the presence of splicing enhancers. The quantity of precipitated U1-C protein was similar for all noncanonical 5'ss dinucleotides, so that the highly different 5'ss usage was likely due to a later step after early U1 snRNP binding. In addition, we show that an internal GT at positions +5/+6 can be advantageous for splicing at position +1 of noncanonical splice sites. Likewise, and in agreement with previous observations, splicing inactive U1 snRNP binding sites could serve as splicing enhancers, which may also explain the higher abundance of U1 snRNPs compared to other U snRNPs. Finally, we observe that an arginine-serine (RS)-rich domain recruitment to stem loop I of the U1 snRNA is functionally sufficient to promote exon-definition and upstream 3'ss activation.

Animal, Fungi, and Plant Genome Sequences Harbor Different Non-Canonical Splice Sites

Most protein-encoding genes in eukaryotes contain introns, which are interwoven with exons. Introns need to be removed from initial transcripts in order to generate the final messenger RNA (mRNA), which can be translated into an amino acid sequence. Precise excision of introns by the spliceosome requires conserved dinucleotides, which mark the splice sites. However, there are variations of the highly conserved combination of GT at the 5' end and AG at the 3' end of an intron in the genome. GC-AG and AT-AC are two major non-canonical splice site combinations, which have been known for years. Recently, various minor non-canonical splice site combinations were detected with numerous dinucleotide permutations. Here, we expand systematic investigations of non-canonical splice site combinations in plants across eukaryotes by analyzing fungal and animal genome sequences. Comparisons of splice site combinations between these three kingdoms revealed several differences, such as an apparently increased CT-AC frequency in fungal genome sequences. Canonical GT-AG splice site combinations in antisense transcripts are a likely explanation for this observation, thus indicating annotation errors. In addition, high numbers of GA-AG splice site combinations were observed in Eurytemoraaffinis and Oikopleuradioica. A variant in one U1 small nuclear RNA (snRNA) isoform might allow the recognition of GA as a 5' splice site. In depth investigation of splice site usage based on RNA-Seq read mappings indicates a generally higher flexibility of the 3' splice site compared to the 5' splice site across animals, fungi, and plants.

High-throughput analysis revealed mutations' diverging effects on SMN1 exon 7 splicing

Splicing-affecting mutations can disrupt gene function by altering the transcript assembly. To ascertain splicing dysregulation principles, we modified a minigene assay for the parallel high-throughput evaluation of different mutations by next-generation sequencing. In our model system, all exonic and six intronic positions of the SMN1 gene's exon 7 were mutated to all possible nucleotide variants, which amounted to 180 unique single-nucleotide mutants and 470 double mutants. The mutations resulted in a wide range of splicing aberrations. Exonic splicing-affecting mutations resulted either in substantial exon skipping, supposedly driven by predicted exonic splicing silencer or cryptic donor splice site (5'ss) and de novo 5'ss strengthening and use. On the other hand, a single disruption of exonic splicing enhancer was not sufficient to cause major exon skipping, suggesting these elements can be substituted during exon recognition. While disrupting the acceptor splice site led only to exon skipping, some 5'ss mutations potentiated the use of three different cryptic 5'ss. Generally, single mutations supporting cryptic 5'ss use displayed better pre-mRNA/U1 snRNA duplex stability and increased splicing regulatory element strength across the original 5'ss. Analyzing double mutants supported the predominating splicing regulatory elements' effect, but U1 snRNA binding could contribute to the global balance of splicing isoforms. Based on these findings, we suggest that creating a new splicing enhancer across the mutated 5'ss can be one of the main factors driving cryptic 5'ss use.

Ranking noncanonical 5' splice site usage by genome-wide RNA-seq analysis and splicing reporter assays

Most human pathogenic mutations in 5' splice sites affect the canonical GT in positions +1 and +2, leading to noncanonical dinucleotides. On the other hand, noncanonical dinucleotides are observed under physiological conditions in ∼1% of all human 5'ss. It is therefore a challenging task to understand the pathogenic mutation mechanisms underlying the conditions under which noncanonical 5'ss are used. In this work, we systematically examined noncanonical 5' splice site selection, both experimentally using splicing competition reporters and by analyzing a large RNA-seq data set of 54 fibroblast samples from 27 subjects containing a total of 2.4 billion gapped reads covering 269,375 exon junctions. From both approaches, we consistently derived a noncanonical 5'ss usage ranking GC > TT > AT > GA > GG > CT. In our competition splicing reporter assay, noncanonical splicing was strictly dependent on the presence of upstream or downstream splicing regulatory elements (SREs), and changes in SREs could be compensated by variation of U1 snRNA complementarity in the competing 5'ss. In particular, we could confirm splicing at different positions (i.e., -1, +1, +5) of a splice site for all noncanonical dinucleotides "weaker" than GC. In our comprehensive RNA-seq data set analysis, noncanonical 5'ss were preferentially detected in weakly used exon junctions of highly expressed genes. Among high-confidence splice sites, they were 10-fold overrepresented in clusters with a neighboring, more frequently used 5'ss. Conversely, these more frequently used neighbors contained only the dinucleotides GT, GC, and TT, in accordance with the above ranking.

Noncanonical GA and GG 5' Intron Donor Splice Sites Are Common in the Copepod Eurytemora affinis

The noncanonical 5′ intron donor splice sites GA and GG are exceedingly rare in described eukaryotic genomes; however, they are present in ∼12% of introns in the genome of the copepod Eurytemora affinis. Failure to recognize the high frequency of these donor sites compromised the modeling of genes in this newly sequenced genome, including 10 conserved ionotropic glutamate receptor (GluR) family genes curated herein. These introns appear to have been acquired recently, along with many additional idiosyncratic introns. Their high frequency implies the evolution of modified intron donor splice site recognition in this copepod.

A mechanism for a single nucleotide intron shift

Spliceosomal introns can occupy nearby rather than identical positions in orthologous genes (intron sliding or shifting). Stwintrons are complex intervening sequences, where an ‘internal’ intron interrupts one of the sequences essential for splicing, generating after its excision, a newly formed canonical intron defined as ‘external’. In one experimentally demonstrated configuration, two alternatively excised internal introns, overlapping by one G, disrupt respectively the donor and the acceptor sequence of an external intron, leading to mRNAs encoding identical proteins. In a gene encoding a DHA1 antiporter in Pezizomycotina, we find a variety of predicted intron configurations interrupting the DNA stretch encoding a conserved peptidic sequence. Some sport a stwintron where the internal intron interrupts the donor of the external intron (experimentally confirmed for Aspergillus nidulans). In others, we found and demonstrate (for Trichoderma reesei) alternative, overlapping internal introns. Discordant canonical introns, one nt apart, are present in yet other species, exactly as predicted by the alternative loss of either of the internal introns at the DNA level from an alternatively spliced stwintron. An evolutionary pathway of 1 nt intron shift, involving an alternatively spliced stwintron intermediate is proposed on the basis of the experimental and genomic data presented.

Molecular and functional characterization of the BMPR2 gene in Pulmonary Arterial Hypertension

Pulmonary arterial hypertension is a progressive disease that causes the obstruction of precapillary pulmonary arteries and a sustained increase in pulmonary vascular resistance. The aim was to analyze functionally the variants found in the BMPR2 gene and to establish a genotype-phenotype correlation. mRNA expression studies were performed using pSPL3 vector, studies of subcellular localization were performed using pEGFP-N1 vector and luciferase assays were performed using pGL3-Basic vector. We have identified 30 variants in the BMPR2 gene in 27 of 55 patients. In 16 patients we detected pathogenic mutations. Minigene assays revealed that 6 variants (synonymous, missense) result in splicing defect. By immunofluorescence assay, we observed that 4 mutations affect the protein localization. Finally, 4 mutations located in the 5'UTR region showed a decreased transcriptional activity in luciferase assays. Genotype-phenotype correlation, revealed that patients with pathogenic mutations have a more severe phenotype (sPaP p = 0.042, 6MWT p = 0.041), a lower age at diagnosis (p = 0.040) and seemed to have worse response to phosphodiesterase-5-inhibitors (p = 0.010). Our study confirms that in vitro expression analysis is a suitable approach in order to investigate the phenotypic consequences of the nucleotide variants, especially in cases where the involved genes have a pattern of expression in tissues of difficult access.

Noncanonical registers and base pairs in human 5' splice-site selection

Accurate recognition of splice sites is essential for pre-messenger RNA splicing. Mammalian 5' splice sites are mainly recognized by canonical base-pairing to the 5' end of U1 small nuclear RNA, yet we described multiple noncanonical base-pairing registers by shifting base-pair positions or allowing one-nucleotide bulges. By systematic mutational and suppressor U1 analyses, we prove three registers involving asymmetric loops and show that two-nucleotide bulges but not longer can form in this context. Importantly, we established that a noncanonical uridine-pseudouridine interaction in the 5' splice site/U1 helix contributes to the recognition of certain 5' splice sites. Thermal melting experiments support the formation of noncanonical registers and uridine-pseudouridine interactions. Overall, we experimentally validated or discarded the majority of predicted noncanonical registers, to derive a list of 5' splice sites using such alternative mechanisms that is much different from the original. This study allows not only the mechanistic understanding of the recognition of a wide diversity of mammalian 5' splice sites, but also the future development of better splice-site scoring methods that reliably predict the effects of disease-causing mutations at these sequences.

Therapeutic strategies based on modified U1 snRNAs and chaperones for Sanfilippo C splicing mutations

BACKGROUND

Mutations affecting RNA splicing represent more than 20% of the mutant alleles in Sanfilippo syndrome type C, a rare lysosomal storage disorder that causes severe neurodegeneration. Many of these mutations are localized in the conserved donor or acceptor splice sites, while few are found in the nearby nucleotides.

METHODS

In this study we tested several therapeutic approaches specifically designed for different splicing mutations depending on how the mutations affect mRNA processing. For three mutations that affect the donor site (c.234 + 1G > A, c.633 + 1G > A and c.1542 + 4dupA), different modified U1 snRNAs recognizing the mutated donor sites, have been developed in an attempt to rescue the normal splicing process. For another mutation that affects an acceptor splice site (c.372-2A > G) and gives rise to a protein lacking four amino acids, a competitive inhibitor of the HGSNAT protein, glucosamine, was tested as a pharmacological chaperone to correct the aberrant folding and to restore the normal trafficking of the protein to the lysosome.

RESULTS

Partial correction of c.234 + 1G > A mutation was achieved with a modified U1 snRNA that completely matches the splice donor site suggesting that these molecules may have a therapeutic potential for some splicing mutations. Furthermore, the importance of the splice site sequence context is highlighted as a key factor in the success of this type of therapy. Additionally, glucosamine treatment resulted in an increase in the enzymatic activity, indicating a partial recovery of the correct folding.

CONCLUSIONS

We have assayed two therapeutic strategies for different splicing mutations with promising results for the future applications.

A comprehensive survey of non-canonical splice sites in the human transcriptome

We uncovered the diversity of non-canonical splice sites at the human transcriptome using deep transcriptome profiling. We mapped a total of 3.7 billion human RNA-seq reads and developed a set of stringent filters to avoid false non-canonical splice site detections. We identified 184 splice sites with non-canonical dinucleotides and U2/U12-like consensus sequences. We selected 10 of the herein identified U2/U12-like non-canonical splice site events and successfully validated 9 of them via reverse transcriptase-polymerase chain reaction and Sanger sequencing. Analyses of the 184 U2/U12-like non-canonical splice sites indicate that 51% of them are not annotated in GENCODE. In addition, 28% of them are conserved in mouse and 76% are involved in alternative splicing events, some of them with tissue-specific alternative splicing patterns. Interestingly, our analysis identified some U2/U12-like non-canonical splice sites that are converted into canonical splice sites by RNA A-to-I editing. Moreover, the U2/U12-like non-canonical splice sites have a differential distribution of splicing regulatory sequences, which may contribute to their recognition and regulation. Our analysis provides a high-confidence group of U2/U12-like non-canonical splice sites, which exhibit distinctive features among the total human splice sites.

The early manifestation, tumor-specific occurrence and prognostic significance of TGFBR2 aberrant splicing in oral carcinoma

Alternative splicing is an important mechanism that can disrupt cell cycle control resulting in tumorigenesis. Although many alterations of Transforming Growth Factor Beta (TGFβ) signaling are reported in cancers, the role of splice aberrations in destabilizing this signaling is the least understood mechanism. In this study, we compared TGFBR2 alternative splicing events in potentially malignant oral disorders (PMDs) and oral squamous cell carcinoma (OSCC) samples with those in normal samples. Interestingly, there were five alternatively spliced forms of TGFBR2 with a deficient kinase domain in OSCCs. The TGFBR2 aberrant splicing was tumor-specific, suggesting that selective splicing out of TGFBR2 kinase domain could be a mechanism misused by cancer cells for evading TGFβ signaling-mediated anti-tumor activities. Moreover, these aberrant transcripts were present in PMDs as well, suggesting an early occurrence of these events during oral carcinogenesis and offering the possibility of early diagnosis of malignancy. Furthermore, OSCC patients who harbored these aberrantly spliced transcripts exhibited poor disease free survival (p=0.028) and poor overall survival (p=0.013). Thus, assessing the presence of these TGFBR2 transcripts can serve as a prognostic marker for oral cancer.

Pick one, but be quick: 5' splice sites and the problems of too many choices

Splice site selection is fundamental to pre-mRNA splicing and the expansion of genomic coding potential. 5' Splice sites (5'ss) are the critical elements at the 5' end of introns and are extremely diverse, as thousands of different sequences act as bona fide 5'ss in the human transcriptome. Most 5'ss are recognized by base-pairing with the 5' end of the U1 small nuclear RNA (snRNA). Here we review the history of research on 5'ss selection, highlighting the difficulties of establishing how base-pairing strength determines splicing outcomes. We also discuss recent work demonstrating that U1 snRNA:5'ss helices can accommodate noncanonical registers such as bulged duplexes. In addition, we describe the mechanisms by which other snRNAs, regulatory proteins, splicing enhancers, and the relative positions of alternative 5'ss contribute to selection. Moreover, we discuss mechanisms by which the recognition of numerous candidate 5'ss might lead to selection of a single 5'ss and propose that protein complexes propagate along the exon, thereby changing its physical behavior so as to affect 5'ss selection.

Diverse captive non-human primates with phytanic acid-deficient diets rich in plant products have substantial phytanic acid levels in their red blood cells

BACKGROUND

Humans and rodents with impaired phytanic acid (PA) metabolism can accumulate toxic stores of PA that have deleterious effects on multiple organ systems. Ruminants and certain fish obtain PA from the microbial degradation of dietary chlorophyll and/or through chlorophyll-derived precursors. In contrast, humans cannot derive PA from chlorophyll and instead normally obtain it only from meat, dairy, and fish products.

RESULTS

Captive apes and Old world monkeys had significantly higher red blood cell (RBC) PA levels relative to humans when all subjects were fed PA-deficient diets. Given the adverse health effects resulting from PA over accumulation, we investigated the molecular evolution of thirteen PA metabolism genes in apes, Old world monkeys, and New world monkeys. All non-human primate (NHP) orthologs are predicted to encode full-length proteins with the marmoset Phyh gene containing a rare, but functional, GA splice donor dinucleotide. Acox2, Scp2, and Pecr sequences had amino acid positions with accelerated substitution rates while Amacr had significant variation in evolutionary rates in apes relative to other primates.

CONCLUSIONS

Unlike humans, diverse captive NHPs with PA-deficient diets rich in plant products have substantial RBC PA levels. The favored hypothesis is that NHPs can derive significant amounts of PA from the degradation of ingested chlorophyll through gut fermentation. If correct, this raises the possibility that RBC PA levels could serve as a biomarker for evaluating the digestive health of captive NHPs. Furthermore, the evolutionary rates of the several genes relevant to PA metabolism provide candidate genetic adaptations to NHP diets.

Widespread recognition of 5' splice sites by noncanonical base-pairing to U1 snRNA involving bulged nucleotides

An established paradigm in pre-mRNA splicing is the recognition of the 5' splice site (5'ss) by canonical base-pairing to the 5' end of U1 small nuclear RNA (snRNA). We recently reported that a small subset of 5'ss base-pair to U1 in an alternate register that is shifted by 1 nucleotide. Using genetic suppression experiments in human cells, we now demonstrate that many other 5'ss are recognized via noncanonical base-pairing registers involving bulged nucleotides on either the 5'ss or U1 RNA strand, which we term "bulge registers." By combining experimental evidence with transcriptome-wide free-energy calculations of 5'ss/U1 base-pairing, we estimate that 10,248 5'ss (∼5% of human 5'ss) in 6577 genes use bulge registers. Several of these 5'ss occur in genes with mutations causing genetic diseases and are often associated with alternative splicing. These results call for a redefinition of an essential element for gene expression that incorporates these registers, with important implications for the molecular classification of splicing mutations and for alternative splicing.

FRS2α is essential for the fibroblast growth factor to regulate the mTOR pathway and autophagy in mouse embryonic fibroblasts

Although the fibroblast growth factor (FGF) signaling axis plays important roles in cell survival, proliferation, and differentiation, the molecular mechanism underlying how the FGF elicits these diverse regulatory signals is not well understood. By using the Frs2α null mouse embryonic fibroblast (MEF) in conjunction with inhibitors to multiple signaling pathways, here we report that the FGF signaling axis activates mTOR via the FGF receptor substrate 2α (FRS2α)-mediated PI3K/Akt pathway, and suppresses autophagy activity in MEFs. In addition, the PI3K/Akt pathway regulated mTOR is crucial for the FGF signaling axis to suppress autophagy in MEFs. Since autophagy has been proposed to play important roles in cell survival, proliferation, and differentiation, the findings suggest a novel mechanism for the FGF signaling axis to transmit regulatory signals to downstream effectors.

Characterization of two closely related α-amylase paralogs in the bark beetle, Ips typographus (L.)

Ips typographus (L.), the eight-spined spruce bark beetle, causes severe damage throughout Eurasian spruce forests and suitable nuclear markers are needed in order to study its population structure on a genetic level. Two closely related genes encoding α-amylase in I. typographus were characterized and named AmyA and AmyB. Both α-amylase paralogs consisted of six exons and five introns. AmyA encodes a polypeptide of 483 amino acids, whereas AmyB has two alternative transcripts encoding polypeptides of 483 and 370 amino acids. The expression levels of both genes were high during larval stage and adulthood. The AmyB transcripts were absent in the pupal stage. A modification of the allozyme staining method allowed us to detect two clusters of bands on the electrophoretic gel that may correspond to the two α-amylase genes. There was a correlation between the lack of AmyB expression in pupa and the absence of the fast migrating isozyme cluster at this stage, suggesting that the faster migrating isoforms are products of the AmyB gene, whereas the slowly migrating bands are derived from the AmyA.

Recognition of atypical 5' splice sites by shifted base-pairing to U1 snRNA

Accurate pre-mRNA splicing is critical for gene expression. The 5' splice site (5' ss) — the highly diverse element at the 5' end of introns — is initially recognized via base-pairing to the 5' end of U1 small nuclear RNA (snRNA). However, many natural 5' ss have a very poor match to the consensus sequence, and are predicted to be very weak. Using genetic suppression experiments in human cells, we demonstrate that some atypical 5' ss are actually efficiently recognized by U1, in an alternative base-pairing register that is shifted by one nucleotide. These atypical 5' ss are phylogenetically widespread, and many of them are conserved. Moreover, shifted base-pairing provides an explanation for the effect of a 5' ss mutation associated with pontocerebellar hypoplasia. The unexpected flexibility in 5' ss/U1 base-pairing challenges an established paradigm, and has broad implications for splice-site prediction algorithms and gene-annotation efforts in genome projects.

Comparative analysis of sequence features involved in the recognition of tandem splice sites

Background

The splicing of pre-mRNAs is conspicuously often variable and produces multiple alternatively spliced (AS) isoforms that encode different messages from one gene locus. Computational studies uncovered a class of highly similar isoforms, which were related to tandem 5'-splice sites (5'ss) and 3'-splice sites (3'ss), yet with very sparse anecdotal evidence in experimental studies. To compare the types and levels of alternative tandem splice site exons occurring in different human organ systems and cell types, and to study known sequence features involved in the recognition and distinction of neighboring splice sites, we performed large-scale, stringent alignments of cDNA sequences and ESTs to the human and mouse genomes, followed by experimental validation.

Results

We analyzed alternative 5'ss exons (A5Es) and alternative 3'ss exons (A3Es), derived from transcript sequences that were aligned to assembled genome sequences to infer patterns of AS occurring in several thousands of genes. Comparing the levels of overlapping (tandem) and non-overlapping (competitive) A5Es and A3Es, a clear preference of isoforms was seen for tandem acceptors and donors, with four nucleotides and three to six nucleotides long exon extensions, respectively. A subset of inferred A5E tandem exons was selected and experimentally validated. With the focus on A5Es, we investigated their transcript coverage, sequence conservation and base-paring to U1 snRNA, proximal and distal splice site classification, candidate motifs for cis-regulatory activity, and compared A5Es with A3Es, constitutive and pseudo-exons, in H. sapiens and M. musculus. The results reveal a small but authentic enriched set of tandem splice site preference, with specific distances between proximal and distal 5'ss (3'ss), which showed a marked dichotomy between the levels of in- and out-of-frame splicing for A5Es and A3Es, respectively, identified a number of candidate NMD targets, and allowed a rough estimation of a number of undetected tandem donors based on splice site information.

Conclusion

This comparative study distinguishes tandem 5'ss and 3'ss, with three to six nucleotides long extensions, as having unusually high proportions of AS, experimentally validates tandem donors in a panel of different human tissues, highlights the dichotomy in the types of AS occurring at tandem splice sites, and elucidates that human alternative exons spliced at overlapping 5'ss posses features of typical splice variants that could well be beneficial for the cell.

Widespread and subtle: alternative splicing at short-distance tandem sites

Alternative splicing at donor or acceptor sites located just a few nucleotides apart is widespread in many species. It results in subtle changes in the transcripts and often in the encoded proteins. Several of these tandem splice events contribute to the repertoire of functionally different proteins, whereas many are neutral or deleterious. Remarkably, some of the functional events are differentially spliced in tissues or developmental stages, whereas others exhibit constant splicing ratios, indicating that function is not always associated with differential splicing. Stochastic splice site selection seems to play a major role in these processes. Here, we review recent progress in understanding functional and evolutionary aspects as well as the mechanism of splicing at short-distance tandem sites.

Violating the splicing rules: TG dinucleotides function as alternative 3' splice sites in U2-dependent introns

TG dinucleotides functioning as alternative 3' splice sites were identified and experimentally verified in 36 human genes.

Background

Despite some degeneracy of sequence signals that govern splicing of eukaryotic pre-mRNAs, it is an accepted rule that U2-dependent introns exhibit the 3' terminal dinucleotide AG. Intrigued by anecdotal evidence for functional non-AG 3' splice sites, we carried out a human genome-wide screen.

Results

We identified TG dinucleotides functioning as alternative 3' splice sites in 36 human genes. The TG-derived splice variants were experimentally validated with a success rate of 92%. Interestingly, ratios of alternative splice variants are tissue-specific for several introns. TG splice sites and their flanking intron sequences are substantially conserved between orthologous vertebrate genes, even between human and frog, indicating functional relevance. Remarkably, TG splice sites are exclusively found as alternative 3' splice sites, never as the sole 3' splice site for an intron, and we observed a distance constraint for TG-AG splice site tandems.

Conclusion

Since TGs splice sites are exclusively found as alternative 3' splice sites, the U2 spliceosome apparently accomplishes perfect specificity for 3' AGs at an early splicing step, but may choose 3' TGs during later steps. Given the tiny fraction of TG 3' splice sites compared to the vast amount of non-viable TGs, cis-acting sequence signals must significantly contribute to splice site definition. Thus, we consider TG-AG 3' splice site tandems as promising subjects for studies on the mechanisms of 3' splice site selection.

Fibroblast growth factor receptor 1 (FGFR1) tyrosine phosphorylation regulates binding of FGFR substrate 2alpha (FRS2alpha) but not FRS2 to the receptor

Binding of the fibroblast growth factor (FGF) to the FGF receptor (FGFR) tyrosine kinase leads to receptor tyrosine autophosphorylation as well as phosphorylation of multiple downstream signaling molecules that are recruited to the receptor either by direct binding or through adaptor proteins. The FGFR substrate 2 (FRS2) family consists of two members, FRS2alpha and FRS2beta, and has been shown to recruit multiple signaling molecules, including Grb2 and Shp2, to FGFR1. To better understand how FRS2 interacted with FGFR1, in vivo binding assays with coexpressed FGFR1 and FRS2 recombinant proteins in mammalian cells were carried out. The results showed that the interaction of full-length FRS2alpha, but not FRS2beta, with FGFR1 was enhanced by activation of the receptor kinase. The truncated FRS2alpha mutant that was comprised only of the phosphotyrosine-binding domain (PTB) bound FGFR1 constitutively, suggesting that the C-terminal sequence downstream the PTB domain inhibited the PTB-FGFR1 binding. Inactivation of the FGFR1 kinase and substitutions of tyrosine phosphorylation sites of FGFR1, but not FRS2alpha, reduced binding of FGFR1 with FRS2alpha. The results suggest that although the tyrosine autophosphorylation sites of FGFR1 did not constitute the binding sites for FRS2alpha, phosphorylation of these residues was essential for optimal interaction with FRS2alpha. In addition, it was demonstrated that the Grb2-binding sites of FRS2alpha are essential for mediating signals of FGFR1 to activate the FiRE enhancer of the mouse syndecan 1 gene. The results, for the first time, demonstrate the specific signals mediated by the Grb2-binding sites and further our understanding of FGF signal transmission at the adaptor level.

Aberrant 5' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

Despite a growing number of splicing mutations found in hereditary diseases, utilization of aberrant splice sites and their effects on gene expression remain challenging to predict. We compiled sequences of 346 aberrant 5′splice sites (5′ss) that were activated by mutations in 166 human disease genes. Mutations within the 5′ss consensus accounted for 254 cryptic 5′ss and mutations elsewhere activated 92 de novo 5′ss. Point mutations leading to cryptic 5′ss activation were most common in the first intron nucleotide, followed by the fifth nucleotide. Substitutions at position +5 were exclusively G>A transitions, which was largely attributable to high mutability rates of C/G>T/A. However, the frequency of point mutations at position +5 was significantly higher than that observed in the Human Gene Mutation Database, suggesting that alterations of this position are particularly prone to aberrant splicing, possibly due to a requirement for sequential interactions with U1 and U6 snRNAs. Cryptic 5′ss were best predicted by computational algorithms that accommodate nucleotide dependencies and not by weight-matrix models. Discrimination of intronic 5′ss from their authentic counterparts was less effective than for exonic sites, as the former were intrinsically stronger than the latter. Computational prediction of exonic de novo 5′ss was poor, suggesting that their activation critically depends on exonic splicing enhancers or silencers. The authentic counterparts of aberrant 5′ss were significantly weaker than the average human 5′ss. The development of an online database of aberrant 5′ss will be useful for studying basic mechanisms of splice-site selection, identifying splicing mutations and optimizing splice-site prediction algorithms.

Phylogenetically widespread alternative splicing at unusual GYNGYN donors

BACKGROUND

Splice donor sites have a highly conserved GT or GC dinucleotide and an extended intronic consensus sequence GTRAGT that reflects the sequence complementarity to the U1 snRNA. Here, we focus on unusual donor sites with the motif GYNGYN (Y stands for C or T; N stands for A, C, G, or T).

RESULTS

While only one GY functions as a splice donor for the majority of these splice sites in human, we provide computational and experimental evidence that 110 (1.3%) allow alternative splicing at both GY donors. The resulting splice forms differ in only three nucleotides, which results mostly in the insertion/deletion of one amino acid. However, we also report the insertion of a stop codon in four cases. Investigating what distinguishes alternatively from not alternatively spliced GYNGYN donors, we found differences in the binding to U1 snRNA, a strong correlation between U1 snRNA binding strength and the preferred donor, over-represented sequence motifs in the adjacent introns, and a higher conservation of the exonic and intronic flanks between human and mouse. Extending our genome-wide analysis to seven other eukaryotic species, we found alternatively spliced GYNGYN donors in all species from mouse to Caenorhabditis elegans and even in Arabidopsis thaliana. Experimental verification of a conserved GTAGTT donor of the STAT3 gene in human and mouse reveals a remarkably similar ratio of alternatively spliced transcripts in both species.

CONCLUSION

In contrast to alternative splicing in general, GYNGYN donors in addition to NAGNAG acceptors enable subtle protein variations.

Genomic splice-site analysis reveals frequent alternative splicing close to the dominant splice site

Alternative pre-mRNA splicing may be the most efficient and widespread mechanism to generate multiple protein isoforms from single genes. Here, we describe the genomic analysis of one of the most frequent types of alternative pre-mRNA splicing, alternative 5'- and 3'-splice-site selection. Using an EST-based alternative splicing database recording >47,000 alternative splicing events, we determined the frequency and location of alternative 5'- and 3'-splice sites within the human genome. The most common alternative splice sites used in the human genome are located within 6 nucleotides (nt) of the dominant splice site. We show that the EST database overrepresents alternative splicing events that maintain the reading frame, thus supporting the concept that RNA quality-control steps ensure that mRNAs that encode for potentially harmful protein products are destroyed and do not serve as templates for translation. The most frequent location for alternative 5'-splice sites is 4 nt upstream or downstream from the dominant splice site. Sequence analysis suggests that this preference is a consequence of the U1 snRNP binding sequence at the 5'-splice site, which frequently contains a GU dinucleotide 4 nt downstream from the dominant splice site. Surprisingly, approximately 50% of duplicated 3'-YAG splice junctions are subject to alternative splicing. This high probability of alternative 3'-splice-site activation in close proximity of the dominant 3'-splice site suggests that the second step of the splicing may be prone to violate splicing fidelity.

Genome annotation past, present, and future: How to define an ORF at each locus

Driven by competition, automation, and technology, the genomics community has far exceeded its ambition to sequence the human genome by 2005. By analyzing mammalian genomes, we have shed light on the history of our DNA sequence, determined that alternatively spliced RNAs and retroposed pseudogenes are incredibly abundant, and glimpsed the apparently huge number of non-coding RNAs that play significant roles in gene regulation. Ultimately, genome science is likely to provide comprehensive catalogs of these elements. However, the methods we have been using for most of the last 10 years will not yield even one complete open reading frame (ORF) for every gene--the first plateau on the long climb toward a comprehensive catalog. These strategies--sequencing randomly selected cDNA clones, aligning protein sequences identified in other organisms, sequencing more genomes, and manual curation--will have to be supplemented by large-scale amplification and sequencing of specific predicted mRNAs. The steady improvements in gene prediction that have occurred over the last 10 years have increased the efficacy of this approach and decreased its cost. In this Perspective, I review the state of gene prediction roughly 10 years ago, summarize the progress that has been made since, argue that the primary ORF identification methods we have relied on so far are inadequate, and recommend a path toward completing the Catalog of Protein Coding Genes, Version 1.0.

Extended base pair complementarity between U1 snRNA and the 5' splice site does not inhibit splicing in higher eukaryotes, but rather increases 5' splice site recognition

Spliceosome formation is initiated by the recognition of the 5′ splice site through formation of an RNA duplex between the 5′ splice site and U1 snRNA. We have previously shown that RNA duplex formation between U1 snRNA and the 5′ splice site can protect pre-mRNAs from degradation prior to splicing. This initial RNA duplex must be disrupted to expose the 5′ splice site sequence for base pairing with U6 snRNA and to form the active spliceosome. Here, we investigated whether hyperstabilization of the U1 snRNA/5′ splice site duplex interferes with splicing efficiency in human cell lines or nuclear extracts. Unlike observations in Saccharomyces cerevisiae, we demonstrate that an extended U1 snRNA/5′ splice site interaction does not decrease splicing efficiency, but rather increases 5′ splice site recognition and exon inclusion. However, low complementarity of the 5′ splice site to U1 snRNA significantly increases exon skipping and RNA degradation. Although the splicing mechanisms are conserved between human and S.cerevisiae, these results demonstrate that distinct differences exist in the activation of the spliceosome.

Determinants of the inherent strength of human 5' splice sites

We previously showed that the authentic 5' splice site (5'ss) of the first exon in the human beta-globin gene is intrinsically stronger than a cryptic 5'ss located 16 nucleotides upstream. Here we examined by mutational analysis the contribution of individual 5'ss nucleotides to discrimination between these two 5'ss. Based on the in vitro splicing efficiencies of a panel of 26 wild-type and mutant substrates in two separate 5'ss competition assays, we established a hierarchy of 5'ss and grouped them into three functional subclasses: strong, intermediate, and weak. Competition between two 5'ss from different subclasses always resulted in selection of the 5'ss that belongs to the stronger subclass. Moreover, each subclass has different characteristic features. Strong and intermediate 5'ss can be distinguished by their predicted free energy of base-pairing to the U1 snRNA 5' terminus (DeltaG). Whereas the extent of splicing via the strong 5'ss correlates well with the DeltaG, this is not the case for competition between intermediate 5'ss. Weak 5'ss were used only when the competing authentic 5'ss was inactivated by mutation. These results indicate that extensive complementarity to U1 snRNA exerts a dominant effect for 5'ss selection, but in the case of competing 5'ss with similarly modest complementarity to U1, the role of other 5'ss features is more prominent. This study reveals the importance of additional submotifs present in certain 5'ss sequences, whose characterization will be critical for understanding 5'ss selection in human genes.

hnRNP H binding at the 5' splice site correlates with the pathological effect of two intronic mutations in the NF-1 and TSHbeta genes

We have recently reported a disease-causing substitution (+5G > C) at the donor site of NF-1 exon 3 that produces its skipping. We have now studied in detail the splicing mechanism involved in analyzing RNA-protein complexes at several 5' splice sites. Characteristic protein patterns were observed by pulldown and band-shift/super-shift analysis. Here, we show that hnRNP H binds specifically to the wild-type GGGgu donor sequence of the NF-1 exon 3. Depletion analyses shows that this protein restricts the accessibility of U1 small nuclear ribonucleoprotein (U1snRNA) to the donor site. In this context, the +5G > C mutation abolishes both U1snRNP base pairing and the 5' splice site (5'ss) function. However, exon recognition in the mutant can be rescued by disrupting the binding of hnRNP H, demonstrating that this protein enhances the effects of the +5G > C substitution. Significantly, a similar situation was found for a second disease-causing +5G > A substitution in the 5'ss of TSHbeta exon 2, which harbors a GGgu donor sequence. Thus, the reason why similar nucleotide substitutions can be either neutral or very disruptive of splicing function can be explained by the presence of specific binding signatures depending on local contexts.

Intrinsic differences between authentic and cryptic 5' splice sites

Cryptic splice sites are used only when use of a natural splice site is disrupted by mutation. To determine the features that distinguish authentic from cryptic 5' splice sites (5'ss), we systematically analyzed a set of 76 cryptic 5'ss derived from 46 human genes. These cryptic 5'ss have a similar frequency distribution in exons and introns, and are usually located close to the authentic 5'ss. Statistical analysis of the strengths of the 5'ss using the Shapiro and Senapathy matrix revealed that authentic 5'ss have significantly higher score values than cryptic 5'ss, which in turn have higher values than the mutant ones. beta-Globin provides an interesting exception to this rule, so we chose it for detailed experimental analysis in vitro. We found that the sequences of the beta-globin authentic and cryptic 5'ss, but not their surrounding context, determine the correct 5'ss choice, although their respective scores do not reflect this functional difference. Our analysis provides a statistical basis to explain the competitive advantage of authentic over cryptic 5'ss in most cases, and should facilitate the development of tools to reliably predict the effect of disease-associated 5'ss-disrupting mutations at the mRNA level.