Long-range RNA pairings contribute to mutually exclusive splicing

Mutually exclusive splicing is an important means of increasing the protein repertoire, by which the Down''s syndrome cell adhesion molecule (Dscam) gene potentially generates 38,016 different isoforms in Drosophila melanogaster. However, the regulatory mechanisms remain obscure due to the complexity of the Dscam exon cluster. Here, we reveal a molecular model for the regulation of the mutually exclusive splicing of the serpent pre-mRNA based on competition between upstream and downstream RNA pairings. Such dual RNA pairings confer fine tuning of the inclusion of alternative exons. Moreover, we demonstrate that the splicing outcome of alternative exons is mediated in relative pairing strength-correlated mode. Combined comparative genomics analysis and experimental evidence revealed similar bidirectional structural architectures in exon clusters 4 and 9 of the Dscam gene. Our findings provide a novel mechanistic framework for the regulation of mutually exclusive splicing and may offer potentially applicable insights into long-range RNA–RNA interactions in gene regulatory networks.     

A DNA damage signal activates and derepresses exon inclusion in Drosophila TAF1 alternative splicing

Signal-dependent alternative splicing is important for regulating gene expression in eukaryotes, yet our understanding of how signals impact splicing mechanisms is limited. A model to address this issue is alternative splicing of Drosophila TAF1 pre-mRNA in response to camptothecin (CPT)-induced DNA damage signals. CPT treatment of Drosophila S2 cells causes increased inclusion of TAF1 alternative cassette exons 12a and 13a through an ATR signaling pathway. To evaluate the role of TAF1 pre-mRNA sequences in the alternative splicing mechanism, we developed a TAF1 minigene (miniTAF1) and an S2 cell splicing assay that recapitulated key aspects of CPT-induced alternative splicing of endogenous TAF1. Analysis of miniTAF1 indicated that splice site strength underlies independent and distinct mechanisms that control exon 12a and 13a inclusion. Mutation of the exon 13a weak 5' splice site or weak 3' splice site to a consensus sequence was sufficient for constitutive exon 13a inclusion. In contrast, mutation of the exon 12a strong 5' splice site or moderate 3' splice site to a consensus sequence was only sufficient for constitutive exon 12a inclusion in the presence of CPT-induced signals. Analogous studies of the exon 13 3' splice site suggest that exon 12a inclusion involves signal-dependent pairing between constitutive and alternative splice sites. Finally, intronic elements identified by evolutionary conservation were necessary for full repression of exon 12a inclusion or full activation of exon 13a inclusion and may be targets of CPT-induced signals. In summary, this work defines the role of sequence elements in the regulation of TAF1 alternative splicing in response to a DNA damage signal.     

Massive expansions of Dscam splicing diversity via staggered homologous recombination during arthropod evolution

The arthropod Down syndrome cell adhesion molecule (Dscam) gene can generate tens of thousands of protein isoforms via combinatorial splicing of numerous alternative exons encoding immunoglobulin variable domains organized into three clusters referred to as the exon 4, 6, and 9 clusters. Dscam protein diversity is important for nervous system development and immune functions. We have performed extensive phylogenetic analyses of Dscam from 20 arthropods (each containing between 46 and 96 alternative exons) to reconstruct the detailed history of exon duplication and loss events that built this remarkable system over 450 million years of evolution. Whereas the structure of the exon 4 cluster is ancient, the exon 6 and 9 clusters have undergone massive, independent expansions in each insect lineage. An analysis of nearly 2000 duplicated exons enabled detailed reconstruction of the timing, location, and boundaries of these duplication events. These data clearly show that new Dscam exons have arisen continuously throughout arthropod evolution and that this process is still occurring in the exon 6 and 9 clusters. Recently duplicated regions display boundaries corresponding to a single exon and the adjacent intron. The boundaries, homology, location, clustering, and relative frequencies of these duplication events strongly suggest that staggered homologous recombination is the major mechanism by which new Dscam exons evolve. These data provide a remarkably detailed picture of how complex gene structure evolves and reveal the molecular mechanism behind this process.     

Autoregulation of Fox protein expression to produce dominant negative splicing factors

The Fox proteins are a family of regulators that control the alternative splicing of many exons in neurons, muscle, and other tissues. Each of the three mammalian paralogs, Fox-1 (A2BP1), Fox-2 (RBM9), and Fox-3 (HRNBP3), produces proteins with a single RNA-binding domain (RRM) flanked by N- and C-terminal domains that are highly diversified through the use of alternative promoters and alternative splicing patterns. These genes also express protein isoforms lacking the second half of the RRM (FoxΔRRM), due to the skipping of a highly conserved 93-nt exon. Fox binding elements overlap the splice sites of these exons in Fox-1 and Fox-2, and the Fox proteins themselves inhibit exon inclusion. Unlike other cases of splicing autoregulation by RNA-binding proteins, skipping the RRM exon creates an in-frame deletion in the mRNA to produce a stable protein. These FoxΔRRM isoforms expressed from cDNA exhibit highly reduced binding to RNA in vivo. However, we show that they can act as repressors of Fox-dependent splicing, presumably by competing with full-length Fox isoforms for interaction with other splicing factors. Interestingly, the Drosophila Fox homolog contains a nearly identical exon in its RRM domain that also has flanking Fox-binding sites. Thus, rather than autoregulation of splicing controlling the abundance of the regulator, the Fox proteins use a highly conserved mechanism of splicing autoregulation to control production of a dominant negative isoform.     

Multisite and bidirectional exonic splicing enhancer in CD44 alternative exon v3

The human CD44 gene encodes multiple isoforms of a transmembrane protein that differ in their extracellular domains as a result of alternative splicing of its variable exons. Expression of CD44 is tightly regulated according to the type and physiological status of a cell, with expression of high molecular weight isoforms by inclusion of variable exons and low molecular weight isoforms containing few or no variable exons. Human CD44 variable exon 3 (v3) can follow a specific alternative splicing route different from that affecting other variable exons. Here we map and functionally describe the splicing enhancer element within CD44 exon v3 which regulates its inclusion in the final mRNA. The v3 splicing enhancer is a multisite bipartite element consisting of a tandem nonamer, the XX motif, and an heptamer, the Y motif, located centrally in the exon. Each of the three sites of this multisite enhancer partially retains its splicing enhancing capacity independently from each other in CD44 and shows full enhancing function in gene contexts different from CD44. We further demonstrate that these motifs act cooperatively as at least two motifs are needed to maintain exon inclusion. Their action is differential with respect to the splice-site target abutting v3. The first X motif acts on the 3' splice site, the second X motif acts on both splice sites (as a bidirectional exonic splicing enhancer), and the Y motif acts on the 5' splice site. We also show that the multisite v3 splicing enhancer is functional irrespective of flanking intron length and spatial organization within v3.     

A cell-based screen for splicing regulators identifies hnRNP LL as a distinct signal-induced repressor of CD45 variable exon 4

The human CD45 gene encodes a protein–tyrosine phosphatase that exhibits differential isoform expression in resting and activated T cells due to alternative splicing of three variable exons. Previously, we have used biochemical methods to identify two regulatory proteins, hnRNP L and PSF, which contribute to the activation-induced skipping of CD45 via the ESS1 regulatory element in variable exon 4. Here we report the identification of a third CD45 regulatory factor, hnRNP L-like (hnRNP LL), via a cell-based screen for clonal variants that exhibit an activation-like phenotype of CD45 splicing even under resting conditions. Microarray analysis of two splicing-altered clones revealed increased expression of hnRNP LL relative to wild-type cells. We further demonstrate that both the expression of hnRNP LL protein and its binding to ESS1 are up-regulated in wild-type cells upon activation. Forced overexpression of hnRNP LL in wild-type cells results in an increase in exon repression, while knock-down of hnRNP LL eliminates activation-induced exon skipping. Interestingly, analysis of the binding of hnRNP L and hnRNP LL to mutants of ESS1 reveals that these proteins have overlapping, but distinct binding requirements. Together, these data establish that hnRNP LL plays a critical and unique role in the signal-induced regulation of CD45 and demonstrate the utility of cell-based screens for the identification of novel splicing regulatory factors.     

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

Splicing of vertebrate introns involves recognition of three consensus elements at the 3′ end. The branch point (BP) and polypyrimidine tract (PPT) are usually located within 40 nucleotides (nt) of the 3′ splice site (3′ ss), AG, but can be much more distant. A characteristic of the region between distant BPs (dBPs) and the 3′ ss is the absence of intervening AG dinucleotides, leading to its designation as the “AG exclusion zone” (AGEZ). The human HTR4 gene, which encodes serotonin receptor 4 and has been associated with schizophrenia, bipolar disease, and gastrointestinal disorders, has four exons with extensive AGEZs. We have mapped the BPs for HTR4 exons 3, 4, 5, and g generated by in vitro splicing, and validated them by mutagenesis in exon-trapping vectors. All exons used dBPs up to 273 nt upstream of the exon. Strikingly, exons 4 and 5 used combinations of both distant and conventionally located BPs, suggesting that successful splicing of these exons can occur by distinct pathways. Our results emphasize the importance for single nucleotide polymorphism resequencing projects to take account of potential dBPs, as the extended AGEZs are vulnerable to mutations that could affect splicing itself or regulation of alternative splicing.     

Evolution of multipartite mitochondrial genomes in the booklice of the genus Liposcelis (Psocoptera)

Background

The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes.

Results

To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes.

Conclusions

Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-861) contains supplementary material, which is available to authorized users.     

Alternative splicing of Myb-related genes MYR1 and MYR2 may modulate activities through changes in dimerization,localization, or protein folding

Arabidopsis genes MYR1 and MYR2 are regulators of flowering time under low light intensity. These Myb-related genes are expressed as alternative splice variants affected in their coiled-coil and DNA-binding domains. We tested whether alternative splicing could affect dimerization and localization of MYR1 and MYR2, thereby potentially affecting their activity. Using MYR1 as a model for variants within the coiled-coil region, we detected 2 types of homodimers. For MYR2, alternative splicing in the DNA-binding Myb-like domain abolished the ability of MYR2 to dimerize. Alternative splicing in the coiled-coil domain did not affect nuclear localization, as determined by transient expression in tobacco, while alternative splicing in the DNA-binding domain of MYR2 yielded a distinct intranuclear localization pattern that may reflect changes in phosphorylation-dependent protein folding. Thus alternative splicing of these genes may result in changes in dimerization or protein folding resulting in changes in activity and abundance of MYR1 or MYR2 protein.     

RBM10 regulates alternative splicing

RBM10, originally called S1-1, is a nuclear RNA-binding protein with domains characteristic of RNA processing proteins. It has been reported that RBM10 constitutes spliceosome complexes and that RBM5, a close homologue of RBM10, regulates alternative splicing of apoptosis-related genes, Fas and cFLIP. In this study, we examined whether RBM10 has a regulatory function in splicing similar to RBM5, and determined that it indeed regulates alternative splicing of Fas and Bcl-x genes. RBM10 promotes exon skipping of Fas pre-mRNA as well as selection of an internal 5′-splice site in Bcl-x pre-mRNA. We propose a consensus RBM10-binding sequence at 5′-splice sites of target exons and a mechanistic model of RBM10 action in the alternative splicing.     

Convergent origins and rapid evolution of spliced leader trans-splicing in Metazoa: Insights from the Ctenophora and Hydrozoa

Replacement of mRNA 5′ UTR sequences by short sequences trans-spliced from specialized, noncoding, spliced leader (SL) RNAs is an enigmatic phenomenon, occurring in a set of distantly related animal groups including urochordates, nematodes, flatworms, and hydra, as well as in Euglenozoa and dinoflagellates. Whether SL trans-splicing has a common evolutionary origin and biological function among different organisms remains unclear. We have undertaken a systematic identification of SL exons in cDNA sequence data sets from non-bilaterian metazoan species and their closest unicellular relatives. SL exons were identified in ctenophores and in hydrozoan cnidarians, but not in other cnidarians, placozoans, or sponges, or in animal unicellular relatives. Mapping of SL absence/presence obtained from this and previous studies onto current phylogenetic trees favors an evolutionary scenario involving multiple origins for SLs during eumetazoan evolution rather than loss from a common ancestor. In both ctenophore and hydrozoan species, multiple SL sequences were identified, showing high sequence diversity. Detailed analysis of a large data set generated for the hydrozoan Clytia hemisphaerica revealed trans-splicing of given mRNAs by multiple alternative SLs. No evidence was found for a common identity of trans-spliced mRNAs between different hydrozoans. One feature found specifically to characterize SL-spliced mRNAs in hydrozoans, however, was a marked adenosine enrichment immediately 3′ of the SL acceptor splice site. Our findings of high sequence divergence and apparently indiscriminate use of SLs in hydrozoans, along with recent findings in other taxa, indicate that SL genes have evolved rapidly in parallel in diverse animal groups, with constraint on SL exon sequence evolution being apparently rare.     

Mutation analyses in pedigrees and sporadic cases of ethnic Han Chinese Kallmann syndrome patients

Kallmann syndrome, a form of idiopathic hypogonadotropic hypogonadism, is characterized by developmental abnormalities of the reproductive system and abnormal olfaction. Despite association of certain genes with idiopathic hypogonadotropic hypogonadism, the genetic inheritance and expression are complex and incompletely known. In the present study, seven Kallmann syndrome pedigrees in an ethnic Han Chinese population were screened for genetic mutations. The exons and intron–exon boundaries of 19 idiopathic hypogonadotropic hypogonadism (idiopathic hypogonadotropic hypogonadism)-related genes in seven Chinese Kallmann syndrome pedigrees were sequenced. Detected mutations were also tested in 70 sporadic Kallmann syndrome cases and 200 Chinese healthy controls. In pedigrees 1, 2, and 7, the secondary sex characteristics were poorly developed and the patients’ sense of smell was severely or completely lost. We detected a genetic mutation in five of the seven pedigrees: homozygous KAL1 p.R191ter (pedigree 1); homozygous KAL1 p.C13ter (pedigree 2; a novel mutation); heterozygous FGFR1 p.R250W (pedigree 3); and homozygous PROKR2 p.Y113H (pedigrees 4 and 5). No genetic change of the assayed genes was detected in pedigrees 6 and 7. Among the 70 sporadic cases, we detected one homozygous and one heterozygous PROKR2 p.Y113H mutation. This mutation was also detected heterozygously in 2/200 normal controls and its pathogenicity is likely questionable. The genetics and genotype–phenotype relationships in Kallmann syndrome are complicated. Classical monogenic inheritance does not explain the full range of genetic inheritance of Kallmann syndrome patients. Because of stochastic nature of genetic mutations, exome analyses of Kallmann syndrome patients may provide novel insights.     

Neuropeptide Receptors NPR-1 and NPR-2 Regulate Caenorhabditis elegans Avoidance Response to the Plant Stress Hormone Methyl Salicylate

Alternative splicing is prevalent in plants, but little is known about its regulation in the context of developmental and signaling pathways. We describe here a new factor that influences pre-messengerRNA (mRNA) splicing and is essential for embryonic development in Arabidopsis thaliana. This factor was retrieved in a genetic screen that identified mutants impaired in expression of an alternatively spliced GFP reporter gene. In addition to the known spliceosomal component PRP8, the screen recovered Arabidopsis RTF2 (AtRTF2), a previously uncharacterized, evolutionarily conserved protein containing a replication termination factor 2 (Rtf2) domain. A homozygous null mutation in AtRTF2 is embryo lethal, indicating that AtRTF2 is an essential protein. Quantitative RT-PCR demonstrated that impaired expression of GFP in atrtf2 and prp8 mutants is due to inefficient splicing of the GFP pre-mRNA. A genome-wide analysis using RNA sequencing indicated that 13–16% of total introns are retained to a significant degree in atrtf2 mutants. Considering these results and previous suggestions that Rtf2 represents an ubiquitin-related domain, we discuss the possible role of AtRTF2 in ubiquitin-based regulation of pre-mRNA splicing.     

Structure, Dynamics, and RNA Interaction Analysis of the Human SBDS Protein

Shwachman-Bodian-Diamond syndrome is an autosomal recessive genetic syndrome with pleiotropic phenotypes, including pancreatic deficiencies, bone marrow dysfunctions with increased risk of myelodysplasia or leukemia, and skeletal abnormalities. This syndrome has been associated with mutations in the SBDS gene, which encodes a conserved protein showing orthologs in Archaea and eukaryotes. The Shwachman-Bodian-Diamond syndrome pleiotropic phenotypes may be an indication of different cell type requirements for a fully functional SBDS protein. RNA-binding activity has been predicted for archaeal and yeast SBDS orthologs, with the latter also being implicated in ribosome biogenesis. However, full-length SBDS orthologs function in a species-specific manner, indicating that the knowledge obtained from model systems may be of limited use in understanding major unresolved issues regarding SBDS function, namely, the effect of mutations in human SBDS on its biochemical function and the specificity of RNA interaction. We determined the solution structure and backbone dynamics of the human SBDS protein and describe its RNA binding site using NMR spectroscopy. Similarly to the crystal structures of Archaea, the overall structure of human SBDS comprises three well-folded domains. However, significant conformational exchange was observed in NMR dynamics experiments for the flexible linker between the N-terminal domain and the central domain, and these experiments also reflect the relative motions of the domains. RNA titrations monitored by heteronuclear correlation experiments and chemical shift mapping analysis identified a classic RNA binding site at the N-terminal FYSH (fungal, Yhr087wp, Shwachman) domain that concentrates most of the mutations described for the human SBDS.