Intrasplicing coordinates alternative first exons with alternative splicing in the protein 4.1R gene

Splicing the Difference: Harnessing the Complexity of the Transcriptome in Hematopoiesis

Alternative splicing has long been recognized as a powerful tool to expand the diversity of the transcriptome and the proteome. The study of hematopoiesis, from hematopoietic stem cell maintenance and differentiation into committed progenitors to maturation into functional blood cells, has led the field of stem cell research and cellular differentiation for decades. The importance of aberrant splicing due to mutations in cis has been exemplified in thalassemias, resulting from aberrant expression of β-globin. The simultaneous development of increasingly sophisticated technologies, in particular the combination of multicolor flow cytometric cell sorting with bulk and single-cell sequencing, has provided sophisticated insights into the complex regulation of the blood system. The recognition that mutations in key splicing factors drive myeloid malignancies, in particular myelodysplastic syndromes, has galvanized research into alternative splicing in hematopoiesis and its diseases. In this review, we will update the audience on the exciting novel technologies, highlight alternative splicing events and their regulators with essential functions in hematopoiesis, and provide a high-level overview how splicing factor mutations contribute to hematologic malignancies.

The Exon Junction Complex Core Represses Cancer-Specific Mature mRNA Re-splicing: A Potential Key Role in Terminating Splicing

Using TSG101 pre-mRNA, we previously discovered cancer-specific re-splicing of mature mRNA that generates aberrant transcripts/proteins. The fact that mRNA is aberrantly re-spliced in various cancer cells implies there must be an important mechanism to prevent deleterious re-splicing on the spliced mRNA in normal cells. We thus postulated that mRNA re-splicing is controlled by specific repressors, and we searched for repressor candidates by siRNA-based screening for mRNA re-splicing activity. We found that knock-down of EIF4A3, which is a core component of the exon junction complex (EJC), significantly promoted mRNA re-splicing. Remarkably, we could recapitulate cancer-specific mRNA re-splicing in normal cells by knock-down of any of the core EJC proteins, EIF4A3, MAGOH, or RBM8A (Y14), implicating the EJC core as the repressor of mRNA re-splicing often observed in cancer cells. We propose that the EJC core is a critical mRNA quality control factor to prevent over-splicing of mature mRNA.

Unannotated splicing regulatory elements in deep intron space

Deep intron space harbors a diverse array of splicing regulatory elements that cooperate with better-known exon-proximal elements to enforce proper tissue-specific and development-specific pre-mRNA processing. Many deep intron elements have been highly conserved through vertebrate evolution, yet remain poorly annotated in the human genome. Recursive splicing exons (RS-exons) and intraexons promote noncanonical, multistep resplicing pathways in long introns, involving transient intermediate structures that are greatly underrepresented in RNA-seq datasets. Decoy splice sites and decoy exons act at a distance to inhibit splicing catalysis at annotated splice sites, with functional consequences such as exon skipping and intron retention. RNA:RNA bridges can juxtapose distant sequences within or across introns to activate deep intron splicing enhancers and silencers, to loop out exons to be skipped, or to select one member of a mutually exclusive set of exons. Similarly, protein bridges mediated by interactions among transcript-bound RNA binding proteins (RBPs) can modulate splicing outcomes. Experimental disruption of deep intron elements serving any of these functions can abrogate normal splicing, strongly suggesting that natural mutations of deep intron elements can do likewise to cause human disease. Understanding noncanonical splicing pathways and discovering deep intron regulatory signals, many of which map hundreds to many thousands of nucleotides from annotated splice junctions, is of great academic interest for basic scientists studying alternative splicing mechanisms. Hopefully, this knowledge coupled with increased analysis of deep intron sequences will also have important medical applications, as better interpretation of deep intron mutations may reveal new disease mechanisms and suggest new therapies. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing.

Antisense targeting of decoy exons can reduce intron retention and increase protein expression in human erythroblasts

The decoy exon model has been proposed to regulate a subset of intron retention (IR) events involving predominantly larger introns (>1 kb). Splicing reporter studies have shown that decoy splice sites are essential for activity, suggesting that decoys act by engaging intron-terminal splice sites and competing with cross-intron interactions required for intron excision. The decoy model predicts that antisense oligonucleotides may be able to block decoy splice sites in endogenous pre-mRNA, thereby reducing IR and increasing productive gene expression. Indeed, we now demonstrate that targeting a decoy 5' splice site in the O-GlcNAc transferase (OGT) gene reduced IR from ∼80% to ∼20% in primary human erythroblasts, accompanied by increases in spliced OGT RNA and OGT protein expression. The remaining OGT IR was refractory to antisense treatment and might be mediated by independent mechanism(s). In contrast, other retained introns were strongly dependent on decoy function, since antisense targeting of decoy 5' splice sites greatly reduced (SNRNP70) or nearly eliminated (SF3B1) IR in two widely expressed splicing factors, and also greatly reduced IR in transcripts encoding the erythroid-specific structural protein, α-spectrin (SPTA1). These results show that modulating decoy exon function can dramatically alter IR and suggest that dynamic regulation of decoy exons could be a mechanism to fine-tune gene expression post-transcriptionally in many cell types.

A spliceosomal twin intron (stwintron) participates in both exon skipping and evolutionary exon loss

Spliceosomal twin introns (stwintrons) are introns where any of the three consensus sequences involved in splicing is interrupted by another intron (internal intron). In Aspergillus nidulans, a donor-disrupted stwintron (intron-1) is extant in the transcript encoding a reticulon-like protein. The orthologous transcript of Aspergillus niger can be alternatively spliced; the exon downstream the stwintron could be skipped by excising a sequence that comprises this stwintron, the neighbouring intron-2, and the exon bounded by these. This process involves the use of alternative 3' splice sites for the internal intron, the resulting alternative intervening sequence being a longer 3'-extended stwintron. In 29 species of Onygenales, a multi-step splicing process occurs in the orthologous transcript, in which a complex intervening sequence including the stwintron and neigbouring intron-2, generates by three splicing reactions a "second order intron" which must then be excised with a fourth splicing event. The gene model in two species can be envisaged as one canonical intron (intron-1) evolved from this complex intervening sequence of nested canonical introns found elsewhere in Onygenales. Postulated splicing intermediates were experimentally verified in one or more species. This work illustrates a role of stwintrons in both alternative splicing and the evolution of intron structure.

RNA splicing during terminal erythropoiesis

PURPOSE OF REVIEW

Erythroid progenitors must accurately and efficiently splice thousands of pre-mRNAs as the cells undergo extensive changes in gene expression and cellular remodeling during terminal erythropoiesis. Alternative splicing choices are governed by interactions between RNA binding proteins and cis-regulatory binding motifs in the RNA. This review will focus on recent studies that define the genome-wide scope of splicing in erythroblasts and discuss what is known about its regulation.

RECENT FINDINGS

RNA-seq analysis of highly purified erythroblast populations has revealed an extensive program of alternative splicing of both exons and introns. During normal erythropoiesis, stage-specific splicing transitions alter the structure and abundance of protein isoforms required for optimized red cell production. Mutation or deficiency of splicing regulators underlies hematopoietic disease in myelopdysplasia syndrome patients via disrupting the splicing program.

SUMMARY

Erythroid progenitors execute an elaborate alternative splicing program that modulates gene expression posttranscriptionally, ultimately regulating the structure and function of the proteome in a differentiation stage-specific manner during terminal erythropoiesis. This program helps drive differentiation and ensure synthesis of the proper protein isoforms required to produce mechanically stable red cells. Mutation or deficiency of key splicing regulatory proteins disrupts the splicing program to cause disease.

Circulating primitive erythroblasts establish a functional, protein 4.1R-dependent cytoskeletal network prior to enucleating

Hematopoietic ontogeny is characterized by distinct primitive and definitive erythroid lineages. Definitive erythroblasts mature and enucleate extravascularly and form a unique membrane skeleton, composed of spectrin, 4.1R-complex, and ankyrinR-complex components, to survive the vicissitudes of the adult circulation. However, little is known about the formation and composition of the membrane skeleton in primitive erythroblasts, which progressively mature while circulating in the embryonic bloodstream. We found that primary primitive erythroblasts express the major membrane skeleton genes present in similarly staged definitive erythroblasts, suggesting that the composition and formation of this membrane network is conserved in maturing primitive and definitive erythroblasts despite their respective intravascular and extravascular locations. Membrane deformability and stability of primitive erythroblasts, assayed by microfluidic studies and fluorescence imaged microdeformation, respectively, significantly increase prior to enucleation. These functional changes coincide with protein 4.1 R isoform switching and protein 4.1R-null primitive erythroblasts fail to establish normal membrane stability and deformability. We conclude that maturing primitive erythroblasts initially navigate the embryonic vasculature prior to establishing a deformable cytoskeleton, which is ultimately formed prior to enucleation. Formation of an erythroid-specific, protein 4.1R-dependent membrane skeleton is an important feature not only of definitive, but also of primitive, erythropoiesis in mammals.

Cutting a Long Intron Short: Recursive Splicing and Its Implications

Over time eukaryotic genomes have evolved to host genes carrying multiple exons separated by increasingly larger intronic, mostly non-protein-coding, sequences. Initially, little attention was paid to these intronic sequences, as they were considered not to contain regulatory information. However, advances in molecular biology, sequencing, and computational tools uncovered that numerous segments within these genomic elements do contribute to the regulation of gene expression. Introns are differentially removed in a cell type-specific manner to produce a range of alternatively-spliced transcripts, and many span tens to hundreds of kilobases. Recent work in human and fruitfly tissues revealed that long introns are extensively processed cotranscriptionally and in a stepwise manner, before their two flanking exons are spliced together. This process, called "recursive splicing," often involves non-canonical splicing elements positioned deep within introns, and different mechanisms for its deployment have been proposed. Still, the very existence and widespread nature of recursive splicing offers a new regulatory layer in the transcript maturation pathway, which may also have implications in human disease.

Convergent evolution of twintron-like configurations: One is never enough

Introns inserted within introns are commonly referred to as twintrons, however the original definition for twintron implied that splicing of the external member of the twintron could only proceed upon splicing of the internal member. This review examines the various types of twintron-like arrangements that have been reported and assigns them to either nested or twintron categories that are subdivided further into subtypes based on differences of their mode of splicing. Twintron-like arrangements evolved independently by fortuitous events among different types of introns but once formed they offer opportunities for the evolution of new regulatory strategies and/or novel genetic elements.

Lessons from non-canonical splicing

Recent improvements in experimental and computational techniques that are used to study the transcriptome have enabled an unprecedented view of RNA processing, revealing many previously unknown non-canonical splicing events. This includes cryptic events located far from the currently annotated exons and unconventional splicing mechanisms that have important roles in regulating gene expression. These non-canonical splicing events are a major source of newly emerging transcripts during evolution, especially when they involve sequences derived from transposable elements. They are therefore under precise regulation and quality control, which minimizes their potential to disrupt gene expression. We explain how non-canonical splicing can lead to aberrant transcripts that cause many diseases, and also how it can be exploited for new therapeutic strategies.

Recursive splicing in long vertebrate genes

It is generally believed that splicing removes introns as single units from precursor messenger RNA transcripts. However, some long Drosophila melanogaster introns contain a cryptic site, known as a recursive splice site (RS-site), that enables a multi-step process of intron removal termed recursive splicing. The extent to which recursive splicing occurs in other species and its mechanistic basis have not been examined. Here we identify highly conserved RS-sites in genes expressed in the mammalian brain that encode proteins functioning in neuronal development. Moreover, the RS-sites are found in some of the longest introns across vertebrates. We find that vertebrate recursive splicing requires initial definition of an 'RS-exon' that follows the RS-site. The RS-exon is then excluded from the dominant mRNA isoform owing to competition with a reconstituted 5' splice site formed at the RS-site after the first splicing step. Conversely, the RS-exon is included when preceded by cryptic promoters or exons that fail to reconstitute an efficient 5' splice site. Most RS-exons contain a premature stop codon such that their inclusion can decrease mRNA stability. Thus, by establishing a binary splicing switch, RS-sites demarcate different mRNA isoforms emerging from long genes by coupling cryptic elements with inclusion of RS-exons.

A 130-kDa protein 4.1B regulates cell adhesion, spreading, and migration of mouse embryo fibroblasts by influencing actin cytoskeleton organization

Protein 4.1B is a member of protein 4.1 family, adaptor proteins at the interface of membranes and the cytoskeleton. It is expressed in most mammalian tissues and is known to be required in formation of nervous and cardiac systems; it is also a tumor suppressor with a role in metastasis. Here, we explore functions of 4.1B using primary mouse embryonic fibroblasts (MEF) derived from wild type and 4.1B knock-out mice. MEF cells express two 4.1B isoforms: 130 and 60-kDa. 130-kDa 4.1B was absent from 4.1B knock-out MEF cells, but 60-kDa 4.1B remained, suggesting incomplete knock-out. Although the 130-kDa isoform was predominantly located at the plasma membrane, the 60-kDa isoform was enriched in nuclei. 130-kDa-deficient 4.1B MEF cells exhibited impaired cell adhesion, spreading, and migration; they also failed to form actin stress fibers. Impaired cell spreading and stress fiber formation were rescued by re-expression of the 130-kDa 4.1B but not the 60-kDa 4.1B. Our findings document novel, isoform-selective roles for 130-kDa 4.1B in adhesion, spreading, and migration of MEF cells by affecting actin organization, giving new insight into 4.1B functions in normal tissues as well as its role in cancer.

Spliceosome twin introns in fungal nuclear transcripts

The spliceosome is an RNA/protein complex, responsible for intron excision from eukaryotic nuclear transcripts. In bacteria, mitochondria and plastids, intron excision does not involve the spliceosome, but occurs through mechanisms dependent on intron RNA secondary and tertiary structure. For group II/III chloroplast introns, "twintrons" (introns within introns) have been described. The excision of the external intron, and thus proper RNA maturation, necessitates prior removal of the internal intron, which interrupts crucial sequences of the former. We have here predicted analogous instances of spliceosomal twintrons ("stwintrons") in filamentous fungi. In two specific cases, where the internal intron interrupts the donor of the external intron after the first or after the second nucleotide, respectively, we show that intermediates with the sequence predicted by the "stwintron" hypothesis, are produced in the splicing process. This implies that two successive rounds of RNA scanning by the spliceosome are necessary to produce the mature mRNA. The phylogenetic distributions of the stwintrons we have identified suggest that they derive from "late" events, subsequent to the appearance of the host intron. They may well not be limited to fungal nuclear transcripts, and their generation and eventual disappearance in the evolutionary process are relevant to hypotheses of intron origin and alternative splicing.

Nested introns in an intron: evidence of multi-step splicing in a large intron of the human dystrophin pre-mRNA

The mechanisms by which huge human introns are spliced out precisely are poorly understood. We analyzed large intron 7 (110199 nucleotides) generated from the human dystrophin (DMD) pre-mRNA by RT-PCR. We identified branching between the authentic 5' splice site and the branch point; however, the sequences far from the branch site were not detectable. This RT-PCR product was resistant to exoribonuclease (RNase R) digestion, suggesting that the detected lariat intron has a closed loop structure but contains gaps in its sequence. Transient and concomitant generation of at least two branched fragments from nested introns within large intron 7 suggests internal nested splicing events before the ultimate splicing at the authentic 5' and 3' splice sites. Nested splicing events, which bring the authentic 5' and 3' splice sites into close proximity, could be one of the splicing mechanisms for the extremely large introns.

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.

Re-splicing of mature mRNA in cancer cells promotes activation of distant weak alternative splice sites

Transcripts of the human tumor susceptibility gene 101 (TSG101) are aberrantly spliced in many cancers. A major aberrant splicing event on the TSG101 pre-mRNA involves joining of distant alternative 5′ and 3′ splice sites within exon 2 and exon 9, respectively, resulting in the extensive elimination of the mRNA. The estimated strengths of the alternative splice sites are much lower than those of authentic splice sites. We observed that the equivalent aberrant mRNA could be generated from an intron-less TSG101 gene expressed ectopically in breast cancer cells. Remarkably, we identified a pathway-specific endogenous lariat RNA consisting solely of exonic sequences, predicted to be generated by a re-splicing between exon 2 and exon 9 on the spliced mRNA. Our results provide evidence for a two-step splicing pathway in which the initial constitutive splicing removes all 14 authentic splice sites, thereby bringing the weak alternative splice sites into close proximity. We also demonstrate that aberrant multiple-exon skipping of the fragile histidine triad (FHIT) pre-mRNA in cancer cells occurs via re-splicing of spliced FHIT mRNA. The re-splicing of mature mRNA can potentially generate mutation-independent diversity in cancer transcriptomes. Conversely, a mechanism may exist in normal cells to prevent potentially deleterious mRNA re-splicing events.

Deep intron elements mediate nested splicing events at consecutive AG dinucleotides to regulate alternative 3' splice site choice in vertebrate 4.1 genes

Distal intraexon (iE) regulatory elements in 4.1R pre-mRNA govern 3' splice site choice at exon 2 (E2) via nested splicing events, ultimately modulating expression of N-terminal isoforms of cytoskeletal 4.1R protein. Here we explored intrasplicing in other normal and disease gene contexts and found conservation of intrasplicing through vertebrate evolution. In the paralogous 4.1B gene, we identified ∼120 kb upstream of E2 an ultradistal intraexon, iE(B), that mediates intrasplicing by promoting two intricately coupled splicing events that ensure selection of a weak distal acceptor at E2 (E2dis) by prior excision of the competing proximal acceptor (E2prox). Mutating iE(B) in minigene splicing reporters abrogated intrasplicing, as did blocking endogenous iE(B) function with antisense morpholinos in live mouse and zebrafish animal models. In a human elliptocytosis patient with a mutant 4.1R gene lacking E2 through E4, we showed that aberrant splicing is consistent with iE(R)-mediated intrasplicing at the first available exons downstream of iE(R), namely, alternative E5 and constitutive E6. Finally, analysis of heterologous acceptor contexts revealed a strong preference for nested 3' splice events at consecutive pairs of AG dinucleotides. Distal regulatory elements may control intrasplicing at a subset of alternative 3' splice sites in vertebrate pre-mRNAs to generate proteins with functional diversity.

Insights into the Function of the Unstructured N-Terminal Domain of Proteins 4.1R and 4.1G in Erythropoiesis

Membrane skeletal protein 4.1R is the prototypical member of a family of four highly paralogous proteins that include 4.1G, 4.1N, and 4.1B. Two isoforms of 4.1R (4.1R(135) and 4.1R(80)), as well as 4.1G, are expressed in erythroblasts during terminal differentiation, but only 4.1R(80) is present in mature erythrocytes. One goal in the field is to better understand the complex regulation of cell type and isoform-specific expression of 4.1 proteins. To start answering these questions, we are studying in depth the important functions of 4.1 proteins in the organization and function of the membrane skeleton in erythrocytes. We have previously reported that the binding profiles of 4.1R(80) and 4.1R(135) to membrane proteins and calmodulin are very different despite the similar structure of the membrane-binding domain of 4.1G and 4.1R(135). We have accumulated evidence for those differences being caused by the N-terminal 209 amino acids headpiece region (HP). Interestingly, the HP region is an unstructured domain. Here we present an overview of the differences and similarities between 4.1 isoforms and paralogs. We also discuss the biological significance of unstructured domains.

Efficient in vivo manipulation of alternative pre-mRNA splicing events using antisense morpholinos in mice

Mammalian pre-mRNA alternative splicing mechanisms are typically studied using artificial minigenes in cultured cells, conditions that may not accurately reflect the physiological context of either the pre-mRNA or the splicing machinery. Here, we describe a strategy to investigate splicing of normal endogenous full-length pre-mRNAs under physiological conditions in live mice. This approach employs antisense vivo-morpholinos (vMOs) to mask cis-regulatory sequences or to disrupt splicing factor expression, allowing functional evaluation of splicing regulation in vivo. We applied this strategy to gain mechanistic insight into alternative splicing events involving exons 2 and 16 (E2 and E16) that control the structure and function of cytoskeletal protein 4.1R. In several mouse tissues, inclusion of E16 was substantially inhibited by interfering with a splicing enhancer mechanism using a target protector morpholino that blocked Fox2-dependent splicing enhancers in intron 16 or a splice-blocking morpholino that disrupted Fox2 expression directly. For E2, alternative 3'-splice site choice is coordinated with upstream promoter use across a long 5'-intron such that E1A splices almost exclusively to the distal acceptor (E2dis). vMOs were used to test the in vivo relevance of a deep intron element previously proposed to determine use of E2dis via a two-step intrasplicing model. Two independent vMOs designed against this intronic regulatory element inhibited intrasplicing, robustly switching E1A splicing to the proximal acceptor (E2prox). This finding strongly supports the in vivo physiological relevance of intrasplicing. vMOs represent a powerful tool for alternative splicing studies in vivo and may facilitate exploration of alternative splicing networks in vivo.

Coupled transcription-splicing regulation of mutually exclusive splicing events at the 5' exons of protein 4.1R gene

The tightly regulated production of distinct erythrocyte protein 4.1R isoforms involves differential splicing of 3 mutually exclusive first exons (1A, 1B, 1C) to the alternative 3' splice sites (ss) of exon 2'/2. Here, we demonstrate that exon 1 and 2'/2 splicing diversity is regulated by a transcription-coupled splicing mechanism. We also implicate distinctive regulatory elements that promote the splicing of exon 1A to the distal 3' ss and exon 1B to the proximal 3' ss in murine erythroleukemia cells. A hybrid minigene driven by cytomegalovirus promoter mimicked 1B-promoter-driven splicing patterns but differed from 1A-promoter-driven splicing patterns, suggesting that promoter identity affects exon 2'/2 splicing. Furthermore, splicing factor SF2/ASF ultraviolet (UV) cross-linked to the exon 2'/2 junction CAGAGAA, a sequence that overlaps the distal U2AF(35)-binding 3' ss. Consequently, depletion of SF2/ASF allowed exon 1B to splice to the distal 3' ss but had no effect on exon 1A splicing. These findings identify for the first time that an SF2/ASF binding site also can serve as a 3' ss in a transcript-dependent manner. Taken together, our results suggest that 4.1R gene expression involves transcriptional regulation coupled with a complex splicing regulatory network.

Systematic analysis of alternative promoters correlated with alternative splicing in human genes

Interactions between various events are essential for complex and delicate transcriptional regulation. To delineate the features and potential roles of alternative promoters (APs) correlated with alternative splicing (AS), we have systematically analyzed 9908 putative alternative promoters (PAPs) from 3797 human genes. Our results showed that approximately 65% of AS events are associated with PAPs. Intriguingly, PAPs per human AS gene only averaged 2.6 for our dataset, which was significantly lower than previously reported. This seems to imply that the human genome contains a small pool of appropriable PAPs for AS genes. Exploration of the characteristics of PAPs such as CpG islands, TATA boxes, GC-content, transcription factor binding sites (TFBSs) and repetitive elements suggested that, respectively, 87% and 90% of PAPs of human AS genes are CpG- and TATA box-poor. The GC-content is significantly higher in the downstream of transcription start sites (TSSs) than upstream (58% vs. 53%), and there is a strong negative correlation between the GC-content and the number of PAPs. These suggested that GC-content around the TSSs plays an important role in the regulation of AS. Moreover, different APs contain distinct densities of repetitive elements and TFBSs, indicating that such sequences have an intrinsic role in the divergent regulation of PAPs and AS. Substantial difference was also found between human AS genes in terms of PAP numbers. A close connection between PAPs and AS may play a critical role in the choice of APs and regulation of AS genes. Furthermore, the distribution of AS genes on different human chromosomes also influences the numbers of PAPs and isoforms of AS genes. Our results may provide important clues for further studies on regulatory network of transcription-related events.

Studying alternative splicing regulatory networks through partial correlation analysis

The identification of links between exons and their regulators or targets and between co-spliced exons in human, mouse and rat provides novel insights into the alternative splicing regulatory network.

Background

Alternative pre-mRNA splicing is an important gene regulation mechanism for expanding proteomic diversity in higher eukaryotes. Each splicing regulator can potentially influence a large group of alternative exons. Meanwhile, each alternative exon is controlled by multiple splicing regulators. The rapid accumulation of high-throughput data provides us with a unique opportunity to study the complicated alternative splicing regulatory network.

Results

We propose the use of partial correlation analysis to identify association links between exons and their upstream regulators or their downstream target genes (exon-gene links) and links between co-spliced exons (exon-exon links). The partial correlation analysis avoids taking the ratio of two noisy random variables, exon expression level and gene expression level, so that it achieves a higher statistical power. We named this analysis procedure pCastNet (partial Correlation analysis of splicing transcriptome Network). Through studies of known alternative exons, conservation patterns, relative positions, functional annotations, and RT-PCR experiments, we concluded that pCastNet can effectively identify exon-gene or exon-exon links. We further found that gene pairs with exon-gene or exon-exon links tend to have similar functions or are present in the same pathways. More interestingly, gene pairs with exon-gene or exon-exon links tend to share cis-elements in promoter regions and microRNA binding elements in 3' untranslated regions, which suggests the coupling of co-alternative-splicing, co-transcription-factor-binding, and co-microRNA-binding.

Conclusions

Alternative splicing regulatory networks reconstructed by pCastNet can help us better understand the coordinate and combinatorial nature of alternative splicing regulation. The proposed tool can be readily applied to other high-throughput data such as transcriptome sequencing data.

Mutations of LRTOMT, a fusion gene with alternative reading frames, cause nonsyndromic deafness in humans

Many proteins necessary for sound transduction have been identified through positional cloning of genes that cause deafness. We report here that mutations of LRTOMT are associated with profound nonsyndromic hearing loss at the DFNB63 locus on human chromosome 11q13.3-q13.4. LRTOMT has two alternative reading frames and encodes two different proteins, LRTOMT1 and LRTOMT2, detected by protein blot analyses. LRTOMT2 is a putative methyltransferase. During evolution, new transcripts can arise through partial or complete coalescence of genes. We provide evidence that in the primate lineage LRTOMT evolved from the fusion of two neighboring ancestral genes, which exist as separate genes (Lrrc51 and Tomt) in rodents.

Alternatively spliced exon 5 of the FERM domain of protein 4.1R encodes a novel binding site for erythrocyte p55 and is critical for membrane targeting in epithelial cells

Direct physical linkage of MAGUKs to the actin cytoskeleton was first established by the interaction of erythrocyte p55 with the FERM domain of protein 4.1R. Subsequently, it was reported that p55 binds to a 51-amino acid peptide, encoded by exon 10, located within the FERM domain of protein 4.1R. In this study, we investigated the nature of the p55-FERM domain binding interface and show that p55 binds to a second 35-amino acid peptide, encoded by an alternatively spliced exon 5, located within the FERM domain of protein 4.1R. Competition and Surface Plasmon Resonance-binding measurements suggest that the peptides encoded by exons 5 and 10 bind to independent sites within the D5 domain of p55. Interestingly, the full length 135 kDa isoform of protein 4.1R containing both exons 5 and 10 was targeted exclusively to the plasma membrane of epithelial cells whereas the same isoform without exon 5 completely lost its membrane localization capacity. Together, these results indicate that p55 binds to two distinct sites within the FERM domain, and the alternatively spliced exon 5 is necessary for the membrane targeting of protein 4.1R in epithelial cells. Since sequences similar to the exon 5-peptide of protein 4.1R and D5 domain of p55 are conserved in many proteins, our findings suggest that a similar mechanism may govern the membrane targeting of other FERM domain containing proteins.