Conserved sequence elements associated with exon skipping

Cooperation and Competition of RNA Secondary Structure and RNA-Protein Interactions in the Regulation of Alternative Splicing

The regulation of alternative splicing in eukaryotic cells is carried out through the coordinated action of a large number of factors, including RNA-binding proteins and RNA structure. The RNA structure influences alternative splicing by blocking cis-regulatory elements, or bringing them closer or farther apart. In combination with RNA-binding proteins, it generates transcript conformations that help to achieve the necessary splicing outcome. However, the binding of regulatory proteins depends on RNA structure and, vice versa, the formation of RNA structure depends on the interaction with regulators. Therefore, RNA structure and RNA-binding proteins are inseparable components of common regulatory mechanisms. This review highlights examples of alternative splicing regulation by RNA-binding proteins, the regulation through local and long-range RNA structures, as well as how these elements work together, cooperate, and compete.

Relevance and Regulation of Alternative Splicing in Plant Heat Stress Response: Current Understanding and Future Directions

Alternative splicing (AS) is a major mechanism for gene expression in eukaryotes, increasing proteome diversity but also regulating transcriptome abundance. High temperatures have a strong impact on the splicing profile of many genes and therefore AS is considered as an integral part of heat stress response. While many studies have established a detailed description of the diversity of the RNAome under heat stress in different plant species and stress regimes, little is known on the underlying mechanisms that control this temperature-sensitive process. AS is mainly regulated by the activity of splicing regulators. Changes in the abundance of these proteins through transcription and AS, post-translational modifications and interactions with exonic and intronic cis-elements and core elements of the spliceosomes modulate the outcome of pre-mRNA splicing. As a major part of pre-mRNAs are spliced co-transcriptionally, the chromatin environment along with the RNA polymerase II elongation play a major role in the regulation of pre-mRNA splicing under heat stress conditions. Despite its importance, our understanding on the regulation of heat stress sensitive AS in plants is scarce. In this review, we summarize the current status of knowledge on the regulation of AS in plants under heat stress conditions. We discuss possible implications of different pathways based on results from non-plant systems to provide a perspective for researchers who aim to elucidate the molecular basis of AS under high temperatures.

Secondary structures in RNA synthesis, splicing and translation

Even though the functional role of mRNA molecules is primarily decided by the nucleotide sequence, several properties are determined by secondary structure conformations. Examples of secondary structures include long range interactions, hairpins, R-loops and G-quadruplexes and they are formed through interactions of non-adjacent nucleotides. Here, we discuss advances in our understanding of how secondary structures can impact RNA synthesis, splicing, translation and mRNA half-life. During RNA synthesis, secondary structures determine RNA polymerase II (RNAPII) speed, thereby influencing splicing. Splicing is also determined by RNA binding proteins and their binding rates are modulated by secondary structures. For the initiation of translation, secondary structures can control the choice of translation start site. Here, we highlight the mechanisms by which secondary structures modulate these processes, discuss advances in technologies to detect and study them systematically, and consider the roles of RNA secondary structures in disease.

Genome-Wide Analysis of Alternative Splicing (AS) Mechanism Provides Insights into Salinity Adaptation in the Livers of Three Euryhaline Teleosts, including Scophthalmus maximus, Cynoglossus semilaevis and Oncorhynchus mykiss

Simple Summary

Alternative splicing (AS) is a key post-transcriptional regulatory mechanism that acts an important regulator in response to environmental stimuli in organisms. In the present study, 18 RNA-Seq datasets were utilized to investigate the potential roles of AS in response to different salinity environments in the livers of three euryhaline teleosts, including turbot (Scophthalmus maximus), tongue sole (Cynoglossus semilaevis) and steelhead trout (Oncorhynchus mykiss). The results indicated that different salinity environments changed the splicing patterns of numerous RNA splicing regulators, which might affect the splicing decisions of many downstream target genes in response to salinity changes. This study provides preliminary evidence for the important roles of AS events in salinity adaptation in teleosts.

Abstract

Salinity is an important environmental factor that directly affects the survival of aquatic organisms, including fish. However, the underlying molecular mechanism of salinity adaptation at post-transcriptional regulation levels is still poorly understood in fish. In the present study, 18 RNA-Seq datasets were utilized to investigate the potential roles of alternative splicing (AS) in response to different salinity environments in the livers of three euryhaline teleosts, including turbot (Scophthalmus maximus), tongue sole (Cynoglossus semilaevis) and steelhead trout (Oncorhynchus mykiss). A total of 10,826, 10,741 and 10,112 AS events were identified in the livers of the three species. The characteristics of these AS events were systematically investigated. Furthermore, a total of 940, 590 and 553 differentially alternative splicing (DAS) events were determined and characterized in the livers of turbot, tongue sole and steelhead trout, respectively, between low- and high-salinity environments. Functional enrichment analysis indicated that these DAS genes in the livers of three species were commonly enriched in some GO terms and KEGG pathways associated with RNA processing. The most common DAS genes work as RNA-binding proteins and play crucial roles in the regulation of RNA splicing. The study provides new insights into uncovering the molecular mechanisms of salinity adaptation in teleosts.

Reverse complementary matches simultaneously promote both back-splicing and exon-skipping

Background

Circular RNAs (circRNAs) play diverse roles in different biological and physiological environments and are always expressed in a tissue-specific manner. Especially, circRNAs are enriched in the brain tissues of almost all investigated species, including humans, mice, Drosophila, etc. Although circRNAs were found in C. elegans, the neuron-specific circRNA data is not available yet. Exon-skipping is found to be correlated to circRNA formation, but the mechanisms that link them together are not clear.

Results

Here, through large-scale neuron isolation from the first larval (L1) stage of C. elegans followed by RNA sequencing with ribosomal RNA depletion, the neuronal circRNA data in C. elegans were obtained. Hundreds of novel circRNAs were annotated with high accuracy. circRNAs were highly expressed in the neurons of C. elegans and were positively correlated to the levels of their cognate linear mRNAs. Disruption of reverse complementary match (RCM) sequences in circRNA flanking introns effectively abolished circRNA formation. In the zip-2 gene, deletion of either upstream or downstream RCMs almost eliminated the production of both the circular and the skipped transcript. Interestingly, the 13-nt RCM in zip-2 is highly conserved across five nematode ortholog genes, which show conserved exon-skipping patterns. Finally, through in vivo one-by-one mutagenesis of all the splicing sites and branch points required for exon-skipping and back-splicing in the zip-2 gene, I showed that back-splicing still happened without exon-skipping, and vice versa.

Conclusions

Through protocol optimization, total RNA obtained from sorted neurons is increased to hundreds of nanograms. circRNAs highly expressed in the neurons of C. elegans are more likely to be derived from genes also highly expressed in the neurons. RCMs are abundant in circRNA flanking introns, and RCM-deletion is an efficient way to knockout circRNAs. More importantly, these RCMs are not only required for back-splicing but also promote the skipping of exon(s) to be circularized. Finally, RCMs in circRNA flanking introns can directly promote both exon-skipping and back-splicing, providing a new explanation for the correlation between them.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07910-w.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Disease modeling of core pre-mRNA splicing factor haploinsufficiency

The craniofacial disorder mandibulofacial dysostosis Guion-Almeida type is caused by haploinsufficiency of the U5 snRNP gene EFTUD2/SNU114. However, it is unclear how reduced expression of this core pre-mRNA splicing factor leads to craniofacial defects. Here we use a CRISPR-Cas9 nickase strategy to generate a human EFTUD2-knockdown cell line and show that reduced expression of EFTUD2 leads to diminished proliferative ability of these cells, increased sensitivity to endoplasmic reticulum (ER) stress and the mis-expression of several genes involved in the ER stress response. RNA-Seq analysis of the EFTUD2-knockdown cell line revealed transcriptome-wide changes in gene expression, with an enrichment for genes associated with processes involved in craniofacial development. Additionally, our RNA-Seq data identified widespread mis-splicing in EFTUD2-knockdown cells. Analysis of the functional and physical characteristics of mis-spliced pre-mRNAs highlighted conserved properties, including length and splice site strengths, of retained introns and skipped exons in our disease model. We also identified enriched processes associated with the affected genes, including cell death, cell and organ morphology and embryonic development. Together, these data support a model in which EFTUD2 haploinsufficiency leads to the mis-splicing of a distinct subset of pre-mRNAs with a widespread effect on gene expression, including altering the expression of ER stress response genes and genes involved in the development of the craniofacial region. The increased burden of unfolded proteins in the ER resulting from mis-splicing would exceed the capacity of the defective ER stress response, inducing apoptosis in cranial neural crest cells that would result in craniofacial abnormalities during development.

Intronic RNA: Ad'junk' mediator of post-transcriptional gene regulation

RNA splicing, the process through which intervening segments of noncoding RNA (introns) are excised from pre-mRNAs to allow for the formation of a mature mRNA product, has long been appreciated for its capacity to add complexity to eukaryotic proteomes. However, evidence suggests that the utility of this process extends beyond protein output and provides cells with a dynamic tool for gene regulation. In this review, we aim to highlight the role that intronic RNA plays in mediating specific splicing outcomes in pre-mRNA processing, as well as explore an emerging class of stable intronic sequences that have been observed to act in gene expression control. Building from underlying flexibility in both sequence and structure, intronic RNA provides mechanisms for post-transcriptional gene regulation that are amenable to the tissue and condition specific needs of eukaryotic cells. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Circular exonic RNAs: When RNA structure meets topology

Although RNA circularization was first documented in the 1990s, the extent to which it occurs was not known until recent advances in high-throughput sequencing enabled the widespread identification of circular RNAs (circRNAs). Despite this, many aspects of circRNA biogenesis, structure, and function yet remain obscure. This review focuses on circular exonic RNAs, a subclass of circRNAs that are generated through backsplicing. Here, I hypothesize that RNA secondary structure can be the common factor that promotes both exon skipping and spliceosomal RNA circularization, and that backsplicing of double-stranded regions could generate topologically linked circRNA molecules. CircRNAs manifest themselves by the presence of tail-to-head exon junctions, which were previously attributed to post-transcriptional exon permutation and repetition. I revisit these observations and argue that backsplicing does not automatically imply RNA circularization because tail-to-head exon junctions give only local information about transcript architecture and, therefore, they are in principle insufficient to determine globally circular topology. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

C3a and suPAR drive versican V1 expression in tubular cells of focal segmental glomerulosclerosis

Chronic tubulointerstitial injury impacts the prognosis of focal segmental glomerulosclerosis (FSGS). We found that the level of versican V1 was increased in tubular cells of FSGS patients. Tubular cell-derived versican V1 induced proliferation and collagen synthesis by activating the CD44/Smad3 pathway in fibroblasts. Both urine C3a and suPAR were increased and bound to the tubular cells in FSGS patients. C3a promoted the transcription of versican by activating the AKT/β-catenin pathway. C3aR knockout decreased the expression of versican in Adriamycin-treated (ADR-treated) mice. On the other hand, suPAR bound to integrin β6 and activated Rac1, which bound to SRp40 at the 5' end of exon 7 in versican pre-mRNA. This binding inhibited the 3'-end splicing of intron 6 and the base-pair interactions between intron 6 and intron 8, leading to the formation of versican V1. Cotreatment with ADR and suPAR specifically increased the level of versican V1 in tubulointerstitial tissues and caused more obvious interstitial fibrosis in mice than treatment with only ADR. Altogether, our results show that C3a and suPAR drive versican V1 expression in tubular cells by promoting transcription and splicing, respectively, and the increases in tubular cell-derived versican V1 induce interstitial fibrosis by activating fibroblasts in FSGS.

Insights into circular RNA biology

Circular RNAs (circRNAs) are a novel class of non-coding RNA characterized by a covalently closed-loop structure generated through a special type of alternative splicing termed backsplicing. CircRNAs are emerging as a heterogeneous class of molecules involved in modulating gene expression by regulation of transcription, protein and miRNA functions. CircRNA expression is cell type and tissue specific and can be largely independent of the expression level of the linear host gene, indicating that regulation of expression might be an important aspect with regard to control of circRNA function. In this review, a brief introduction to the characteristics that define a circRNA will be given followed by a discussion of putative biogenesis pathways and modulators of circRNA expression as well as of the stage at which circRNA formation takes place. A brief summary of circRNA functions will also be provided and lastly, an outlook with a focus on unanswered questions regarding circRNA biology will be included.

RNA structure in splicing: An evolutionary perspective

Pre-mRNA splicing is a key post-transcriptional regulation process in which introns are excised and exons are ligated together. A novel class of structured intron was recently discovered in fish. Simple expansions of complementary AC and GT dimers at opposite boundaries of an intron were found to form a bridging structure, thereby enforcing correct splice site pairing across the intron. In some fish introns, the RNA structures are strong enough to bypass the need of regulatory protein factors for splicing. Here, we discuss the prevalence and potential functions of highly structured introns. In humans, structured introns usually arise through the co-occurrence of C and G-rich repeats at intron boundaries. We explore the potentially instructive example of the HLA receptor genes. In HLA pre-mRNA, structured introns flank the exons that encode the highly polymorphic β sheet cleft, making the processing of the transcript robust to variants that disrupt splicing factor binding. While selective forces that have shaped HLA receptor are fairly atypical, numerous other highly polymorphic genes that encode receptors contain structured introns. Finally, we discuss how the elevated mutation rate associated with the simple repeats that often compose structured intron can make structured introns themselves rapidly evolving elements.

The role of short RNA loops in recognition of a single-hairpin exon derived from a mammalian-wide interspersed repeat

Splice-site selection is controlled by secondary structure through sequestration or approximation of splicing signals in primary transcripts but the exact role of even the simplest and most prevalent structural motifs in exon recognition remains poorly understood. Here we took advantage of a single-hairpin exon that was activated in a mammalian-wide interspersed repeat (MIR) by a mutation stabilizing a terminal triloop, with splice sites positioned close to each other in a lower stem of the hairpin. We first show that the MIR exon inclusion in mRNA correlated inversely with hairpin stabilities. Employing a systematic manipulation of unpaired regions without altering splice-site configuration, we demonstrate a high correlation between exon inclusion of terminal tri- and tetraloop mutants and matching tri-/tetramers in splicing silencers/enhancers. Loop-specific exon inclusion levels and enhancer/silencer associations were preserved across primate cell lines, in 4 hybrid transcripts and also in the context of a distinct stem, but only if its loop-closing base pairs were shared with the MIR hairpin. Unlike terminal loops, splicing activities of internal loop mutants were predicted by their intramolecular Watson-Crick interactions with the antiparallel strand of the MIR hairpin rather than by frequencies of corresponding trinucleotides in splicing silencers/enhancers. We also show that splicing outcome of oligonucleotides targeting the MIR exon depend on the identity of the triloop adjacent to their antisense target. Finally, we identify proteins regulating MIR exon recognition and reveal a distinct requirement of adjacent exons for C-terminal extensions of Tra2α and Tra2β RNA recognition motifs.

Transposable element-driven transcript diversification and its relevance to genetic disorders

The human genome project and subsequent gene annotation projects have shown that the human genome contains 22,000-25,000 functional genes. Therefore, it is believed that the diversity of protein repertoire is achieved by the alternative splicing (AS) mechanism. Transposable elements (TEs) are mobile in nature and can therefore alter their position in the genome. The insertion of TEs into a new gene region can result in AS of a particular transcript through various mechanisms, including intron retention, and alternative donor or acceptor splice sites. TE-derived AS is thought to have played a part in primate evolution and in hominid radiation. However, TE-derived AS or genetic instability may sometimes result in genetic disorders. For the past two decades, numerous studies have been performed on TEs and their role in genomes. Accumulating evidence shows that the term 'junk DNA', previously used for TEs is a misnomer. Recent research has indicated that TEs may have clinical potential. However, to explore the feasibility of using TEs in clinical practice, additional studies are required. This review summarizes the available literature on TE-derived AS, alternative promoter, and alternative polyadenylation. The review covers the effects of TEs on coding genes and their clinical implications, and provides our perspectives and directions for future research.

Butyrate Induced IGF2 Activation Correlated with Distinct Chromatin Signatures Due to Histone Modification

Histone modification has emerged as a very important mechanism regulating the transcriptional status of the genome. Insulin-like growth factor 2 (IGF2) is a peptide hormone controlling various cellular processes, including proliferation and apoptosis. H19 gene is closely linked to IGF2 gene, and IGF2 and H19 are reciprocally regulated imprinted genes. The epigenetic signature of H19 promoter (hypermethylation) on the paternal allele plays a vital role in allowing the expression of the paternal allele of IGF2.46 Our previous studies demonstrate that butyrate regulates the expression of IGF2 as well as genes encoding IGF Binding proteins. To obtain further understanding of histone modification and its regulatory potentials in controlling IGF2/H19 gene expression, we investigated the histone modification status of some key histones associated with the expression of IGF2/H19 genes in bovine cells using RNA-seq in combination with Chip-seq technology. A high-resolution map of the major chromatin modification at the IGF2/H19 locus induced by butyrate was constructed to illustrate the fundamental association of the chromatin modification landscape that may play a role in the activation of the IGF2 gene. High-definition epigenomic landscape mapping revealed that IGF2 and H19 have distinct chromatin modification patterns at their coding and promoter regions, such as TSSs and TTSs. Moreover, the correlation between the differentially methylated regions (DMRs) of IGF2/H19 locus and histone modification (acetylation and methylation) indicated that epigenetic signatures/markers of DNA methylation, histone methylation and histone acetylation were differentially distributed on the expressed IGF2 and silenced H19 genes. Our evidence also suggests that butyrate-induced regional changes of histone acetylation statusin the upstream regulation domain of H19 may be related to the reduced expression of H19 and strong activation of IGF2. Our results provided insights into the mechanism of butyrate-induced loss of imprinting (LOI) of IGF2 and regulation of gene expression by histone modification.

Eukaryotic TPP riboswitch regulation of alternative splicing involving long-distance base pairing

Thiamin pyrophosphate (TPP) riboswitches are found in organisms from all three domains of life. Examples in bacteria commonly repress gene expression by terminating transcription or by blocking ribosome binding, whereas most eukaryotic TPP riboswitches are predicted to regulate gene expression by modulating RNA splicing. Given the widespread distribution of eukaryotic TPP riboswitches and the diversity of their locations in precursor messenger RNAs (pre-mRNAs), we sought to examine the mechanism of alternative splicing regulation by a fungal TPP riboswitch from Neurospora crassa, which is mostly located in a large intron separating protein-coding exons. Our data reveal that this riboswitch uses a long-distance (∼530-nt separation) base-pairing interaction to regulate alternative splicing. Specifically, a portion of the TPP-binding aptamer can form a base-paired structure with a conserved sequence element (α) located near a 5′ splice site, which greatly increases use of this 5′ splice site and promotes gene expression. Comparative sequence analyses indicate that many fungal species carry a TPP riboswitch with similar intron architecture, and therefore the homologous genes in these fungi are likely to use the same mechanism. Our findings expand the scope of genetic control mechanisms relying on long-range RNA interactions to include riboswitches.

The high kinetic stability of a G-quadruplex limits hnRNP F qRRM3 binding to G-tract RNA

The RNA binding protein heterogeneous nuclear ribonucleoprotein (hnRNP) F is involved in telomeres maintenance and pre-mRNA processing, such as alternative splicing and polyadenylation. It specifically recognizes RNA containing three consecutive guanines (G-tracts) that have the potential to assemble into G-quadruplexes. We have proposed recently that hnRNP F could regulate alternative splicing by remodeling RNA structures, such as G-quadruplexes. However, the exact mechanism of hnRNP F binding to such RNA sequences remains unknown. Here, we have studied the binding of the third RNA binding domain of hnRNP F [quasi-RNA recognition motif 3 (qRRM3)] to G-tract RNA using isothermal titration calorimetry, circular dichroism and nuclear magnetic resonance spectroscopy. Our results show that qRRM3 binds specifically exclusively to single-stranded G-tracts (ssRNA), in contrast to previous reports stating that the G-quadruplex was recognized as well. Furthermore, we demonstrate that the pre-existent ssRNA/G-quadruplex equilibrium slows down the formation of the protein–ssRNA complex. Based on in vitro transcription assays, we show that the rate of the protein–RNA complex formation is faster than that of the G-quadruplex. We propose a model according to which hnRNP F could bind RNA co-transcriptionally and prevents G-quadruplex formation.

Identification, cloning, and functional characterization of the IL-1 receptor antagonist in the chicken reveal important differences between the chicken and mammals

The human IL-1 family contains 11 genes encoded at three separate loci. Nine, including IL-1R antagonist (IL-1RN), are present at a single locus on chromosome 2, whereas IL-18 and IL-33 lie on chromosomes 11 and 9, respectively. There are currently only two known orthologs in the chicken, IL-1β and IL-18, which are encoded on chromosomes 22 and 24, respectively. Two novel chicken IL-1 family sequences were identified from expressed sequence tag libraries, representing secretory and intracellular (icIL-1RN) structural variants of the IL-1RN gene, as seen in mammals. Two further putative splice variants (SVs) of both chicken IL-1RN (chIL-1RN) structural variants were also isolated. Alternative splicing of human icIL-1RN gives three different transcripts; there are no known SVs for human secretory IL-1RN. The chicken icIL-1RN SVs differ from those found in human icIL-1RN in terms of the rearrangements involved. In mammals, IL-1RN inhibits IL-1 activity by physically occupying the IL-1 type I receptor. Both full-length structural variants of chIL-1RN exhibited biological activity similar to their mammalian orthologs in a macrophage cell line bioassay. The four SVs, however, were not biologically active. The chicken IL-1 family is more fragmented in the genome than those of mammals, particularly in that the large multigene locus seen in mammals is absent. This suggests differential evolution of the family since the divergence of birds and mammals from a common ancestor, and makes determination of the full repertoire of chicken IL-1 family members more challenging.

Post-transcriptional exon shuffling events in humans can be evolutionarily conserved and abundant

In silico analyses have established that transcripts from some genes can be processed into RNAs with rearranged exon order relative to genomic structure (post-transcriptional exon shuffling, or PTES). Although known to contribute to transcriptome diversity in some species, to date the structure, distribution, abundance, and functional significance of human PTES transcripts remains largely unknown. Here, using high-throughput transcriptome sequencing, we identify 205 putative human PTES products from 176 genes. We validate 72 out of 112 products analyzed using RT-PCR, and identify additional PTES products structurally related to 61% of validated targets. Sequencing of these additional products reveals GT-AG dinucleotides at >95% of the splice junctions, confirming that they are processed by the spliceosome. We show that most PTES transcripts are expressed in a wide variety of human tissues, that they can be polyadenylated, and that some are conserved in mouse. We also show that they can extend into 5' and 3' UTRs, consistent with formation via trans-splicing of independent pre-mRNA molecules. Finally, we use real-time PCR to compare the abundance of PTES exon junctions relative to canonical exon junctions within the transcripts from seven genes. PTES exon junctions are present at <0.01% to >90% of the levels of canonical junctions, with transcripts from MAN1A2, PHC3, TLE4, and CDK13 exhibiting the highest levels. This is the first systematic experimental analysis of PTES in human, and it suggests both that the phenomenon is much more widespread than previously thought and that some PTES transcripts could be functional.

Mechanistic control of carcinoembryonic antigen-related cell adhesion molecule-1 (CEACAM1) splice isoforms by the heterogeneous nuclear ribonuclear proteins hnRNP L, hnRNP A1, and hnRNP M

Carcinoembryonic antigen-related cell adhesion molecule-1 (CEACAM1) is expressed in a variety of cell types and is implicated in carcinogenesis. Alternative splicing of CEACAM1 pre-mRNA generates two cytoplasmic domain splice variants characterized by the inclusion (L-isoform) or exclusion (S-isoform) of exon 7. Here we show that the alternative splicing of CEACAM1 pre-mRNA is regulated by novel cis elements residing in exon 7. We report the presence of three exon regulatory elements that lead to the inclusion or exclusion of exon 7 CEACAM1 mRNA in ZR75 breast cancer cells. Heterologous splicing reporter assays demonstrated that the maintenance of authentic alternative splicing mechanisms were independent of the CEACAM1 intron sequence context. We show that forced expression of these exon regulatory elements could alter CEACAM1 splicing in HEK-293 cells. Using RNA affinity chromatography, three members of the heterogeneous nuclear ribonucleoprotein family (hnRNP L, hnRNP A1, and hnRNP M) were identified. RNA immunoprecipitation of hnRNP L and hnRNP A1 revealed a binding motif located central and 3' to exon 7, respectively. Depletion of hnRNP A1 or L by RNAi in HEK-293 cells promoted exon 7 inclusion, whereas overexpression led to exclusion of the variable exon. By contrast, overexpression of hnRNP M showed exon 7 inclusion and production of CEACAM1-L mRNA. Finally, stress-induced cytoplasmic accumulation of hnRNP A1 in MDA-MB-468 cells dynamically alters the CEACAM1-S:CEACAM1:L ratio in favor of the l-isoform. Thus, we have elucidated the molecular factors that control the mechanism of splice-site recognition in the alternative splicing regulation of CEACAM1.

RNA secondary structure in mutually exclusive splicing

Mutually exclusive splicing is a regulated means to generate protein diversity, but the underlying mechanisms are poorly understood. Here comparative genome analysis revealed the built-in intronic elements for controlling mutually exclusive splicing of the 14-3-3ξ pre-mRNA. These elements are clade specific but are evolutionarily conserved at the secondary structure level. Combined evidence revealed the triple functions of these inter-intronic RNA pairings in synergistically ensuring the selection of only one of multiple exons, through activation of the proximal variable exon outside the loop by the approximation of cis elements, and simultaneous repression of the exon within the loop, in combination with the physical competition of RNA pairing. Additionally, under this model, we also deciphered a similar structural code in exon clusters 4 and 9 of Dscam (38,016 isoforms) and Mhc (480 isoforms). Our findings suggest a broadly applicable mechanism to ensure mutually exclusive splicing.

A naturally occurring isoform inhibits parathyroid hormone receptor trafficking and signaling

Parathyroid hormone (PTH) regulates calcium homeostasis and bone remodeling through its cognitive receptor (PTHR). We describe here a PTHR isoform harboring an in-frame 42-bp deletion of exon 14 (Δe14-PTHR) that encodes transmembrane domain 7. Δe14-PTHR was detected in human kidney and buccal epithelial cells. We characterized its topology, cellular localization, and signaling, as well as its interactions with PTHR. The C-terminus of the Δe14-PTHR is extracellular, and cell surface expression is strikingly reduced compared with the PTHR. Δe14-PTHR displayed impaired trafficking and accumulated in endoplasmic reticulum. Signaling and activation of cAMP and ERK by Δe14-PTHR was decreased significantly compared with PTHR. Δe14-PTHR acts as a functional dominant-negative by suppressing the action of PTHR. Cells cotransfected with both receptors exhibit markedly reduced PTHR cell membrane expression, colocalization with Δe14-PTHR in endoplasmic reticulum, and diminished cAMP activation and ERK phosphorylation in response to challenge with PTH. Δe14-PTHR forms heterodimers with PTHR, which may account for cytoplasmic retention of PTHR in the presence of Δe14-PTHR. Analysis of the PTHR heteronuclear RNA suggests that base-pair complementarity in introns surrounding exon 14 causes exon skipping and accounts for generation of the Δe14-PTHR isoform. Thus Δe14-PTHR is a poorly functional receptor that acts as a dominant-negative of PTHR trafficking and signaling and may contribute to PTH resistance.

Modulation of alternative splicing by long-range RNA structures in Drosophila

Accurate and efficient recognition of splice sites during pre-mRNA splicing is essential for proper transcriptome expression. Splice site usage can be modulated by secondary structures, but it is unclear if this type of modulation is commonly used or occurs to a significant degree with secondary structures forming over long distances. Using phlyogenetic comparisons of intronic sequences among 12 Drosophila genomes, we elucidated a group of 202 highly conserved pairs of sequences, each at least nine nucleotides long, capable of forming stable stem structures. This set was highly enriched in alternatively spliced introns and introns with weak acceptor sites and long introns, and most occurred over long distances (>150 nucleotides). Experimentally, we analyzed the splicing of several of these introns using mini-genes in Drosophila S2 cells. Wild-type splicing patterns were changed by mutations that opened the stem structure, and restored by compensatory mutations that re-established the base-pairing potential, demonstrating that these secondary structures were indeed implicated in the splice site choice. Mechanistically, the RNA structures masked splice sites, brought together distant splice sites and/or looped out introns. Thus, base-pairing interactions within introns, even those occurring over long distances, are more frequent modulators of alternative splicing than is currently assumed.

Prediction of alternatively skipped exons and splicing enhancers from exon junction arrays

Background

Alternative splicing of exons in a pre-mRNA transcript is an important mechanism which contributes to protein diversity in human. Arrays for detecting alternative splicing are available using several different probe designs, including those based on exon-junctions. In this work, we introduce a new method for predicting alternatively skipped exons from exon-junction arrays. Predictions based on our method are compared against controls and their sequences are analyzed to identify motifs important for regulating alternative splicing.

Results

Our comparison of several alternative methods shows that an exon-skipping score based on neighboring junctions best discriminates between positive and negative controls. Sequence analysis of our predicted exons confirms the presence of known splicing regulatory sequences. In addition, we also derive a set of development-related alternatively spliced genes based on fetal versus adult tissue comparisons and find that our predictions are consistent with their functional annotations. Ab initio motif finding algorithms are applied to identify several motifs that may be relevant for splicing during development.

Conclusion

This work describes a new method for analyzing exon-junction arrays, identifies sequence motifs that are specific for alternative and constitutive splicing and suggests a role for several known splicing factors and their motifs in developmental regulation.

Intronic Alus influence alternative splicing

Examination of the human transcriptome reveals higher levels of RNA editing than in any other organism tested to date. This is indicative of extensive double-stranded RNA (dsRNA) formation within the human transcriptome. Most of the editing sites are located in the primate-specific retrotransposed element called Alu. A large fraction of Alus are found in intronic sequences, implying extensive Alu-Alu dsRNA formation in mRNA precursors. Yet, the effect of these intronic Alus on splicing of the flanking exons is largely unknown. Here, we show that more Alus flank alternatively spliced exons than constitutively spliced ones; this is especially notable for those exons that have changed their mode of splicing from constitutive to alternative during human evolution. This implies that Alu insertions may change the mode of splicing of the flanking exons. Indeed, we demonstrate experimentally that two Alu elements that were inserted into an intron in opposite orientation undergo base-pairing, as evident by RNA editing, and affect the splicing patterns of a downstream exon, shifting it from constitutive to alternative. Our results indicate the importance of intronic Alus in influencing the splicing of flanking exons, further emphasizing the role of Alus in shaping of the human transcriptome.

Pre-mRNA secondary structures influence exon recognition

The secondary structure of a pre-mRNA influences a number of processing steps including alternative splicing. Since most splicing regulatory proteins bind to single-stranded RNA, the sequestration of RNA into double strands could prevent their binding. Here, we analyzed the secondary structure context of experimentally determined splicing enhancer and silencer motifs in their natural pre-mRNA context. We found that these splicing motifs are significantly more single-stranded than controls. These findings were validated by transfection experiments, where the effect of enhancer or silencer motifs on exon skipping was much more pronounced in single-stranded conformation. We also found that the structural context of predicted splicing motifs is under selection, suggesting a general importance of secondary structures on splicing and adding another level of evolutionary constraints on pre-mRNAs. Our results explain the action of mutations that affect splicing and indicate that the structural context of splicing motifs is part of the mRNA splicing code.

Almost all human protein-coding genes contain several exons and introns. Prior to translation, introns have to be removed and exons have to be joined, which happens in a processing step called splicing that generates the mature mRNA. For most genes, certain exons can be either included or excluded from the mature mRNA. It is currently not fully understood which signals are needed to accurately recognize the boundaries of exons in the intron-containing primary transcript. As in transcriptional regulation, enhancer and silencer sequence motifs are crucial for the correct recognition of exons. Splicing regulatory proteins identify these motifs in a sequence-specific manner. In general, these proteins bind to single-stranded RNA. Here, we analyzed local secondary structures of primary transcripts and found that known splicing motifs are preferentially located in a single-stranded context. Experimental tests demonstrated that motifs in single-stranded contexts have a stronger effect on splice site selection than those located in double-stranded regions. These results help to understand the action of human mutations that change the splicing pattern and indicate that local pre-mRNA secondary structures influence exon recognition.

Bioinformatics of alternative splicing and its regulation

The sequencing of the human genome and ensuing wave of data generation have brought new light upon the extent and importance of alternative splicing as an RNA regulatory mechanism. Alternative splicing could potentially explain the complexity of protein repertoire during evolution, and defects in the splicing mechanism are responsible for diseases as complex as cancer. Among the challenges that rise in light of these discoveries are cataloguing splice variation in the human and other eukaryotic genomes, and identifying and characterizing the splicing regulatory elements that control their expression. Bioinformatics efforts tackling these two questions are just at the beginning. This article is a survey of these methods.

An intronic signal for alternative splicing in the human genome

An important level at which the expression of programmed cell death (PCD) genes is regulated is alternative splicing. Our previous work identified an intronic splicing regulatory element in caspase-2 (casp-2) gene. This 100-nucleotide intronic element, In100, consists of an upstream region containing a decoy 3' splice site and a downstream region containing binding sites for splicing repressor PTB. Based on the signal of In100 element in casp-2, we have detected the In100-like sequences as a family of sequence elements associated with alternative splicing in the human genome by using computational and experimental approaches. A survey of human genome reveals the presence of more than four thousand In100-like elements in 2757 genes. These In100-like elements tend to locate more frequent in intronic regions than exonic regions. EST analyses indicate that the presence of In100-like elements correlates with the skipping of their immediate upstream exons, with 526 genes showing exon skipping in such a manner. In addition, In100-like elements are found in several human caspase genes near exons encoding the caspase active domain. RT-PCR experiments show that these caspase genes indeed undergo alternative splicing in a pattern predicted to affect their functional activity. Together, these results suggest that the In100-like elements represent a family of intronic signals for alternative splicing in the human genome.

ADAM15 gene structure and differential alternative exon use in human tissues

Background

ADAM15 is a metalloprotease-disintegrin implicated in ectodomain shedding and cell adhesion. Aberrant ADAM15 expression has been associated with human cancer and other disorders. We have previously shown that the alternative splicing of ADAM15 transcripts is mis-regulated in cancer cells. To gain a better understanding of ADAM15 regulation, its genomic organization and regulatory elements as well as the alternative exon use in human tissues were characterized.

Results

Human ADAM15, flanked by the FLJ32785/DCST1 and ephrin-A4 genes, spans 11.4 kb from the translation initiation codon to the polyadenylation signal, being the shortest multiple-exon ADAM gene. The gene contains 23 exons varying from 63 to 316 bp and 22 introns from 79 to 1283 bp. The gene appeared to have several transcription start sites and their location suggested the promoter location within a CpG island proximal to the translation start. Reporter expression experiments confirmed the location of functional GC-rich, TATAless and CAATless promoter, with the most critical transcription-supporting elements located -266 to -23 bp relative to the translation start. Normal human tissues showed different complex patterns of at least 13 different ADAM15 splice variants arising from the alternative use of the cytosolic-encoding exons 19, 20a/b, and 21a/b. The deduced ADAM15 protein isoforms have different combinations of cytosolic regulatory protein interaction motifs.

Conclusion

Characterization of human ADAM15 gene and identification of elements involved in the regulation of transcription and alternative splicing provide important clues for elucidation of physiological and pathological roles of ADAM15. The present results also show that the alternative exon use is a physiological post-transcriptional mechanism regulating ADAM15 expression in human tissues.

A correlation with exon expression approach to identify cis-regulatory elements for tissue-specific alternative splicing

Correlation of motif occurrences with gene expression intensity is an effective strategy for elucidating transcriptional cis-regulatory logic. Here we demonstrate that this approach can also identify cis-regulatory elements for alternative pre-mRNA splicing. Using data from a human exon microarray, we identified 56 cassette exons that exhibited higher transcript-normalized expression in muscle than in other normal adult tissues. Intron sequences flanking these exons were then analyzed to identify candidate regulatory motifs for muscle-specific alternative splicing. Correlation of motif parameters with gene-normalized exon expression levels was examined using linear regression and linear splines on RNA words and degenerate weight matrices, respectively. Our unbiased analysis uncovered multiple candidate regulatory motifs for muscle-specific splicing, many of which are phylogenetically conserved among vertebrate genomes. The most prominent downstream motifs were binding sites for Fox1- and CELF-related splicing factors, and a branchpoint-like element acuaac; pyrimidine-rich elements resembling PTB-binding sites were most significant in upstream introns. Intriguingly, our systematic study indicates a paucity of novel muscle-specific elements that are dominant in short proximal intronic regions. We propose that Fox and CELF proteins play major roles in enforcing the muscle-specific alternative splicing program, facilitating expression of unique isoforms of cytoskeletal proteins critical to muscle cell function.

RNA structure is a key regulatory element in pathological ATM and CFTR pseudoexon inclusion events

Genomic variations deep in the intronic regions of pre-mRNA molecules are increasingly reported to affect splicing events. However, there is no general explanation why apparently similar variations may have either no effect on splicing or cause significant splicing alterations. In this work we have examined the structural architecture of pseudoexons previously described in ATM and CFTR patients. The ATM case derives from the deletion of a repressor element and is characterized by an aberrant 5′ss selection despite the presence of better alternatives. The CFTR pseudoexon instead derives from the creation of a new 5′ss that is used while a nearby pre-existing donor-like sequence is never selected. Our results indicate that RNA structure is a major splicing regulatory factor in both cases. Furthermore, manipulation of the original RNA structures can lead to pseudoexon inclusion following the exposure of unused 5′ss already present in their wild-type intronic sequences and prevented to be recognized because of their location in RNA stem structures. Our data show that intrinsic structural features of introns must be taken into account to understand the mechanism of pseudoexon activation in genetic diseases. Our observations may help to improve diagnostics prediction programmes and eventual therapeutic targeting.

Complementary intron sequence motifs associated with human exon repetition: a role for intragenic, inter-transcript interactions in gene expression

MOTIVATION

Exon repetition describes the presence of tandemly repeated exons in mRNA in the absence of duplications in the genome. The regulation of this process is not fully understood. We therefore investigated the entire flanking intronic sequences of exons involved in exon repetition for common sequence elements.

RESULTS

A computational analysis of 48 human single exon repetition events identified two common sequence motifs. One of these motifs is pyrimidine-rich and is more common in the upstream intron, whilst the other motif is highly enriched in purines and is more common in the downstream intron. As the two motifs are complementary to each other, they support a model by which exon repetition occurs as a result of trans-splicing between separate pre-mRNA transcripts from the same gene that are brought together during transcription by complementary intronic sequences. The majority of the motif instances overlap with the locations of mobile elements such as Alu elements. We explore the potential importance of complementary intron sequences in a rat gene that undertakes natural exon repetition in a strain specific manner. The possibility that distant complementary sequences can stimulate inter-transcript splicing during transcription suggests an unsuspected new role for potential secondary structures in endogenous genes.

Splicing of a critical exon of human Survival Motor Neuron is regulated by a unique silencer element located in the last intron

Humans have two nearly identical copies of the Survival Motor Neuron (SMN) gene, SMN1 and SMN2. In spinal muscular atrophy (SMA), SMN2 is not able to compensate for the loss of SMN1 due to exclusion of exon 7. Here we describe a novel inhibitory element located immediately downstream of the 5' splice site in intron 7. We call this element intronic splicing silencer N1 (ISS-N1). Deletion of ISS-N1 promoted exon 7 inclusion in mRNAs derived from the SMN2 minigene. Underlining the dominant role of ISS-N1 in exon 7 skipping, abrogation of a number of positive cis elements was tolerated when ISS-N1 was deleted. Confirming the silencer function of ISS-N1, an antisense oligonucleotide against ISS-N1 restored exon 7 inclusion in mRNAs derived from the SMN2 minigene or from endogenous SMN2. Consistently, this oligonucleotide increased the levels of SMN protein in SMA patient-derived cells that carry only the SMN2 gene. Our findings underscore for the first time the profound impact of an evolutionarily nonconserved intronic element on SMN2 exon 7 splicing. Considering that oligonucleotides annealing to intronic sequences do not interfere with exon-junction complex formation or mRNA transport and translation, ISS-N1 provides a very specific and efficient therapeutic target for antisense oligonucleotide-mediated correction of SMN2 splicing in SMA.

Bb2Bb3 regulation of murine Lyme arthritis is distinct from Ncf1 and independent of the phagocyte nicotinamide adenine dinucleotide phosphate oxidase

Several quantitative trait loci regulating murine Lyme arthritis severity have been mapped, including a highly significant linkage found on chromosome 5, termed Bb2Bb3. Within this region, the Ncf1 gene of the phagocyte nicotinamide adenine dinucleotide phosphate (NADPH) oxidase has recently been identified as a major regulator of arthritis severity in rodent models of rheumatoid arthritis, an effect attributed to protective properties of reactive oxygen species. To assess the role of Ncf1 in Lyme arthritis, we introgressed Bb2Bb3 from severely arthritic C3H/He mice onto mildly arthritic C57BL/6 mice. This increased Lyme arthritis severity, whereas the reciprocal transfer conferred protection from disease. A single nucleotide polymorphism was identified in the Ncf1 gene that did not influence the protein sequence or expression of Ncf1. Although polymorphonuclear leukocytes from C57BL/6 mice generated a greater oxidative burst than polymorphonuclear leukocytes from C3H/He mice, studies with the Bb2Bb3 congenic mice demonstrated this difference was not linked to Ncf1 alleles. Furthermore, Lyme arthritis severity was not altered in mice lacking either the Ncf1 or Gp91phox subunits of the NADPH oxidase complex. Together, these results argue that Ncf1 is not a candidate gene for regulation of Lyme arthritis and reveal Lyme arthritis to be independent of NADPH oxidase activity, distinguishing it from other models of rheumatoid arthritis.

The effect of U1 snRNA binding free energy on the selection of 5' splice sites

The importance of U1 snRNA binding free energy in the regulation of alternative splicing has been studied in some genes with site-directed mutagenesis. Here we report a large-scale analysis of its impact on 5' splice site (5'ss) selection in human genome. The results show that free energy exerts different effects on alternative 5'ss choice in different situations and -8.1 kcal/mol is a threshold. When both free energies of two competing 5'ss are larger than -8.1 kcal/mol, the 5'ss with lower free energy is more frequently used. However, in other pairs of 5'ss, lower-free-energy 5'ss does not seem to be favored and even the other 5'ss is used more frequently, which suggests that very low binding free energy would impair splicing. Some observations hold true only for those alternative 5' splicing with short alternative exons (<50nt), which implies a complex mechanism of 5'ss selection involving both U1 snRNA binding free energy and regulatory factors.

Understanding alternative splicing: towards a cellular code

In violation of the 'one gene, one polypeptide' rule, alternative splicing allows individual genes to produce multiple protein isoforms - thereby playing a central part in generating complex proteomes. Alternative splicing also has a largely hidden function in quantitative gene control, by targeting RNAs for nonsense-mediated decay. Traditional gene-by-gene investigations of alternative splicing mechanisms are now being complemented by global approaches. These promise to reveal details of the nature and operation of cellular codes that are constituted by combinations of regulatory elements in pre-mRNA substrates and by cellular complements of splicing regulators, which together determine regulated splicing pathways.

A duplication in the canine beta-galactosidase gene GLB1 causes exon skipping and GM1-gangliosidosis in Alaskan huskies

GM(1)-gangliosidosis is a lysosomal storage disease that is inherited as an autosomal recessive disorder, predominantly caused by structural defects in the beta-galactosidase gene (GLB1). The molecular cause of GM(1)-gangliosidosis in Alaskan huskies was investigated and a novel 19-bp duplication in exon 15 of the GLB1 gene was identified. The duplication comprised positions +1688-+1706 of the GLB1 cDNA. It partially disrupted a potential exon splicing enhancer (ESE), leading to exon skipping in a fraction of the transcripts. Thus, the mutation caused the expression of two different mRNAs from the mutant allele. One transcript contained the complete exon 15 with the 19-bp duplication, while the other transcript lacked exon 15. In the transcript containing exon 15 with the 19-bp duplication a premature termination codon (PTC) appeared, but due to its localization in the last exon of canine GLB1, nonsense-mediated RNA decay (NMD) did not occur. As a consequence of these molecular events two different truncated GLB1 proteins are predicted to be expressed from the mutant GLB1 allele. In heterozygous carrier animals the wild-type allele produces sufficient amounts of the active enzyme to prevent clinical signs of disease. In affected homozygous dogs no functional GLB1 is synthesized and G(M1)-gangliosidosis occurs.

Multiple Products Derived from Two CCL4 Loci: High Incidence of a New Polymorphism in HIV+Patients

Human CCL4/macrophage inflammatory protein (MIP)-1beta and CCL3/MIP-1alpha are two highly related molecules that belong to a cluster of inflammatory CC chemokines located in chromosome 17. CCL4 and CCL3 were formed by duplication of a common ancestral gene, generating the SCYA4 and SCYA3 genes which, in turn, present a variable number of additional non-allelic copies (SCYA4L and SCYA3L1). In this study, we show that both CCL4 loci (SCYA4 and SCYA4L) are expressed and alternatively generate spliced variants lacking the second exon. In addition, we found that the SCYA4L locus is polymorphic and displays a second allelic variant (hereinafter SCYA4L2) with a nucleotide change in the intron 2 acceptor splice site compared with the one described originally (hereinafter SCYA4L1). Therefore, the pattern of SCYA4L2 transcripts is completely different from that of SCYA4L1, since SCYA4L2 uses several new acceptor splice sites and generates nine new mRNAs. Furthermore, we analyzed the contribution of each locus (SCYA4 and SCYA4L1/L2) to total CCL4 expression in human CD8 T cells by RT-amplified fragment length polymorphism and real-time PCR, and we found that L2 homozygous individuals (L2L2) only express half the levels of CCL4 compared with L1L1 individuals. The analysis of transcripts from the SCYA4L locus showed a lower level in L2 homozygous compared with L1 homozygous individuals (12% vs 52% of total CCL4 transcripts). A possible clinical relevance of these CCL4 allelic variants was suggested by the higher frequency of the L2 allele in a group of HIV(+) individuals (n = 175) when compared with controls (n = 220, 28.6% vs 16.6% (p = 0.00016)).

SpliceInfo: an information repository for mRNA alternative splicing in human genome

We have developed an information repository named SpliceInfo to collect the occurrences of the four major alternative-splicing (AS) modes in human genome; these include exon skipping, 5′-alternative splicing, 3′-alternative splicing and intron retention. The dataset is derived by comparing the nucleotide and protein sequences available for a given gene for evidence of AS. Additional features such as the tissue specificity of the mRNA, the protein domain contained by exons, the GC-ratio of exons, the repeats contained within the exons, and the Gene Ontology are annotated computationally for each exonic region that is alternatively spliced. Motivated by a previous investigation of AS-related motifs such as exonic splicing enhancer and exonic splicing silencer, this resource also provides a means of identifying motifs candidates and this should help to identify potential regulatory mechanisms within a particular exonic sequence set and its two flanking intronic sequence sets. This is carried out using motif discovery tools to identify motif candidates related to alternative splicing regulation and together with a secondary structure prediction tool, will help in the identification of the structural properties of such regulatory motifs. The integrated resource is now available on http://SpliceInfo.mbc.NCTU.edu.tw/.

Information for the Coordinates of Exons (ICE): a human splice sites database

We present a comprehensive database, Information for the Coordinates of Exons (ICE), of genomic splice sites (SSs) for 10,803 human genes. ICE contains 91,846 pairs of donor acceptor sites, supported by the alignment of "full-length" human mRNAs (including transcript variants) on human genomic sequences. ICE represents the largest collection of human SSs known to date and provides a significant resource to both molecular biologists and bioinformaticians alike. A user can visualize and extract genomic sequences around SSs of the donor acceptor pairs and can also visualize the primary structure of individual genes. We list in this article the 22 most frequently found canonical and noncanonical splice sites. The top four most represented donor acceptor pairs (GT-AG, GC-AG, AT-AC, and GT-GG) accounted for 99.16% of our data set. In addition, we calculated the SS matrix models for the three most common donor acceptor pairs. The database is focused on providing SSs and surrounding sequence information, associated SS and sequence characteristics, and relation to overall transcript structure. It allows targeted search and presents evidence for the gene structure.

A survey of alternative transcripts of human tissue kallikrein genes

Alternative splicing is prevalent within the human tissue kallikrein gene locus. Aside from being the most important source of protein diversity in eukaryotes, this process plays a significant role in development, physiology and disease. A better understanding of alternative splicing could lead to the use of gene variants as drug targets, therapeutic agents or diagnostic markers. With the rapidly rising number of alternative kallikrein transcripts, classifying new transcripts and piecing together the significance of existing data are becoming increasingly challenging. In this review, we present a systematic analysis of all currently known kallikrein alternative transcripts. By defining a reference form for each of the 15 kallikrein genes (KLK1 to KLK15), we were able to classify alternative splicing patterns. We identified 82 different kallikrein gene transcript forms, including reference forms. Alternative splicing may lead to the synthesis of 56 different protein forms for KLK1-15. In the kallikrein locus, the majority of alternative splicing events occur within the protein-coding region, and to a lesser extent in the 5' untranslated regions (UTRs). The most common alternative splicing event is exon skipping (35%) and the least common events are cryptic exons (3%) and internal exon deletion (3%). Seventy-six percent of kallikrein splice variants that are predicted to encode truncated proteins are the result of frameshifts. Eighty-nine percent of putative proteins encoded by splice variants are predicted to be secreted. Although several reports describe the identification of kallikrein splice variants and their potential clinical utility, this is the first extensive review on this subject. Accumulating evidence suggests that alternative kallikrein forms could be involved in many pathologic conditions or could have practical applications as biomarkers. The organization and analysis of the kallikrein transcripts will facilitate future work in this area and may lead to novel clinical and diagnostic applications.

Evidence for the regulation of alternative splicing via complementary DNA sequence repeats

MOTIVATION

While the mechanism for regulating alternative splicing is poorly understood, secondary structure has been shown to be integral to this process. Due to their propensity for forming complementary hairpin loops and their elevated mutation rates, tandem repeated sequences have the potential to influence splicing regulation.

RESULTS

An analysis of human intronic sequences reveals a strong correlation between alternative splicing and the prevalence of mono- through hexanucleotide tandem repeats that may engage in complementary pairing in introns that flank alternatively spliced exons. While only 44% of the 18 173 genes in the Human Alternative Splicing Database are known to be alternatively spliced, they contain 84% of the 694 237 intronic complementary repeat pairs. Significantly, the normalized frequency and distribution of repeat sequences, independent of their potential for pairing, are indistinguishable between alternatively spliced and non-alternatively spliced genes. Thus, the increased prevalence of repeats with pairing potential in alternatively spliced genes is not merely a consequence of more repeats or repeat composition bias. These results suggest that complementary repeats may play a role in the regulation of alternative splicing.

CONTACT

harold.garner@utsouthwestern.edu.

Over-representation of exonic splicing enhancers in human intronless genes suggests multiple functions in mRNA processing

The human transcriptome is constituted of a great majority of intron-containing and a minority of intron-lacking mRNAs; given the different processing these transcripts undergo, they are expected to carry, intermingled with coding properties, very different editing information. Here we applied a computational approach to compare intronless and intron-containing coding sequences. Hexamer composition comparison allowed the definition of over- and under-represented motifs in intronless genes; surprisingly, experimental testing revealed that intron-lacking coding sequences are enriched rather than depleted in elements with splicing enhancement ability. Similarly, we show evidence that intronless transcripts display a significantly higher frequency of both shuttling and non-shuttling SR protein binding sites compared to intron-containing sequences. These observations suggest that SR proteins (and possibly other splicing factors) play a role in cellular processes distinct from splicing.

Transcription-coupled and splicing-coupled strand asymmetries in eukaryotic genomes

Under no-strand bias conditions, each genomic DNA strand should present equimolarities of A and T and of G and C. Deviations from these rules are attributed to asymmetric properties intrinsic to DNA mutation-repair processes. In bacteria, strand biases are associated with replication or transcription. In eukaryotes, recent studies demonstrate that human genes present transcription-coupled biases that might reflect transcription-coupled repair processes. Here, we study strand asymmetries in intron sequences of evolutionarily distant eukaryotes, and show that two superimposed intron biases can be distinguished. (i) Biases that are maximum at intron extremities and decrease over large distances to zero values in internal regions, possibly reflecting interactions between pre-mRNA and splicing machinery; these extend over approximately 0.5 kb in mammals and Arabidopsis thaliana, and over 1 kb in Caenorhabditis elegans and Drosophila melanogaster. (ii) Biases that are constant along introns, possibly associated with transcription. Strikingly, in C.elegans, these latter biases extend over intergenic regions that separate co-oriented genes. When appropriately examined, all genomes present transcription-coupled excess of T over A in the coding strand. On the opposite, GC skews are either positive (mammals, plants) or negative (invertebrates). These results suggest that transcription-coupled asymmetries result from mutation-repair mechanisms that differ between vertebrates and invertebrates.

hnRNP A1 and the SR proteins ASF/SF2 and SC35 have antagonistic functions in splicing of beta-tropomyosin exon 6B

Mutually exclusive splicing of exons 6A and 6B from the chicken beta-tropomyosin gene involves numerous regulatory sequences. Previously, we identified a G-rich intronic sequence (S3) downstream of exon 6B. This element consists of six G-rich motifs, mutations of which abolish splicing of exon 6B. In this paper, we investigated the cellular factors that bind to this G-rich element. By using RNA affinity chromatography, we identified heterogeneous nuclear ribonucleoprotein (hnRNP) A1, the SR proteins ASF/SF2 and SC35, and hnRNP F/H as specific components that are assembled onto the G-rich element. By using hnRNP A1-depleted HeLa nuclear extract and add-back experiments, we show that hnRNP A1 has a negative effect on splicing of exon 6B. In agreement with in vitro data, artificial recruitment of hnRNP A1, as a fusion with the MS2 coat protein, also represses splicing of exon 6B ex vivo. In contrast, ASF/SF2 and SC35 activate splicing of exon 6B. As observed with other systems, hnRNP A1 counteracts the stimulating effect of the SR proteins. Moreover, cross-linking experiments show that both ASF/SF2 and SC35 are able to displace binding of hnRNP A1 to the G-rich element, suggesting that the binding sites for these proteins are overlapping. These data indicate that the G-rich sequence is a composite element that acts as an enhancer or as a silencer, depending on which proteins bind to them.