Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events.     

Intronic alternative splicing regulators identified by comparative genomics in nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.     

The tetracycline-responsive promoter contains functional interferon-inducible response elements

Tetracycline (tet)-responsive expression vectors allow controlled inducible expression of proteins in mammalian cells. This system is widely used for experimental research both in vivo and in vitro. In our attempts to use this system to study the antiviral effect of IFNα on hepatitis B virus, we discovered an unexpected feature of the tet-responsive promoter (tet promoter) of the currently available expression vectors. IFNα was found to stimulate tet promoter activity after transient transfection in a dose- and cell type-dependent manner. By sequence inspection, an IFNα-stimulated response element (ISRE)-like sequence was identified in the linker regions located between the heptameric tet operator sequences. Gel shift assays revealed binding of IFN-stimulated gene factors to these sequences, indicating that they mediate the IFNα-mediated promoter stimulation. These data demonstrate an unexpected feature of the tet-responsive expression system which needs to be taken into acount when using this system for analysis of cytokine functions in vitro and in vivo. The data also imply that the tet promoter-based expression system can be rendered non-responsive to IFNα by mutagenesis of the ISREs and this may be essential when considering gene therapy in vivo.     

Quaking and PTB control overlapping splicing regulatory networks during muscle cell differentiation

Alternative splicing contributes to muscle development, but a complete set of muscle-splicing factors and their combinatorial interactions are unknown. Previous work identified ACUAA (“STAR” motif) as an enriched intron sequence near muscle-specific alternative exons such as Capzb exon 9. Mass spectrometry of myoblast proteins selected by the Capzb exon 9 intron via RNA affinity chromatography identifies Quaking (QK), a protein known to regulate mRNA function through ACUAA motifs in 3′ UTRs. We find that QK promotes inclusion of Capzb exon 9 in opposition to repression by polypyrimidine tract-binding protein (PTB). QK depletion alters inclusion of 406 cassette exons whose adjacent intron sequences are also enriched in ACUAA motifs. During differentiation of myoblasts to myotubes, QK levels increase two- to threefold, suggesting a mechanism for QK-responsive exon regulation. Combined analysis of the PTB- and QK-splicing regulatory networks during myogenesis suggests that 39% of regulated exons are under the control of one or both of these splicing factors. This work provides the first evidence that QK is a global regulator of splicing during muscle development in vertebrates and shows how overlapping splicing regulatory networks contribute to gene expression programs during differentiation.     

Cloning and promoter analysis of the chicken interferon regulatory factor-3 gene

Interferon regulatory factors (IRFs) are a family of DNA-binding proteins involved in mediating the cellular response to interferons (IFNs) and viral infection. Although extensively studied in mammals, IRFs of other vertebrates have been less well characterized. Previously, we cloned chicken interferon regulatory factor-3 (chIRF-3) mRNA, which is rapidly and transiently induced by double-stranded (ds)RNA. The chIRF-3 mRNA encodes a protein distinct from any known mammalian IRF. Here, we show that chIRF-3 is activated additively by type I and type II IFNs. To delineate the sequence elements required to regulate chIRF-3 expression, we cloned chlRF-3 and 0.48 kb of 5' flanking sequence. Computer analysis of the proximal promoter revealed three putative binding sites for nuclear factor (NF)-kappaB, two overlapping interferon-stimulated response elements (ISREs), and an interferon gamma activating sequence (GAS). The presence of both GAS and ISRE consensus sequences in the chIRF-3 promoter is unique among IRF family members. Both type I and II IFNs, as well as dsRNA and IRF-1, trans-activate the promoter in short-term transfection experiments. Mutational analysis of the promoter demonstrated that the putative NF-kappaB binding sites are needed for stimulation by dsRNA but not by either type I or type II IFN and that both the overlapping ISREs and GAS are required for full induction by type I or type II IFN.     

Conserved sequences in the final intron of MDM2 are essential for the regulation of alternative splicing of MDM2 in response to stress

Alternative splicing plays a fundamental role in generating proteome diversity and is critical in regulation of eukaryotic gene expression. It is estimated that 50% of disease-causing mutations alter splicing efficiency and/or patterns of splicing. An alternatively spliced form of murine double-minute 2, MDM2-ALT1, is associated with pediatric rhabdomyosarcoma (RMS) at high frequency in primary human tumors and RMS cell lines. We have identified that this isoform can be induced in response to specific types of stress (UV and cisplatin). However, the mechanism of alternative splicing of MDM2 in human cancer is unknown. Using UV and cisplatin to model alternative splicing of the MDM2 gene, we have developed a damage-inducible in vitro splicing system. This system employs an MDM2 minigene that mimics the damage-induced alternative splicing observed in vivo. Using this in vitro splicing system, we have shown that conserved intronic sequences in intron 11 of MDM2 are required for normal splicing. Furthermore, we showed that these intronic elements are also required for the regulated damage-induced alternative splicing of MDM2. The use of this novel damage-inducible system will allow for the systematic identification of regulatory elements and factors involved in the splicing regulation of the MDM2 gene in response to stress. This study has implications for identification of novel intervention points for development of future therapeutics for rhabdomyosarcoma.     

An Alu-derived intronic splicing enhancer facilitates intronic processing and modulates aberrant splicing in ATM

We have previously reported a natural GTAA deletion within an intronic splicing processing element (ISPE) of the ataxia telangiectasia mutated (ATM) gene that disrupts a non-canonical U1 snRNP interaction and activates the excision of the upstream portion of the intron. The resulting pre-mRNA splicing intermediate is then processed to a cryptic exon, whose aberrant inclusion in the final mRNA is responsible for ataxia telangiectasia. We show here that the last 40 bases of a downstream intronic antisense Alu repeat are required for the activation of the cryptic exon by the ISPE deletion. Evaluation of the pre-mRNA splicing intermediate by a hybrid minigene assay indicates that the identified intronic splicing enhancer represents a novel class of enhancers that facilitates processing of splicing intermediates possibly by recruiting U1 snRNP to defective donor sites. In the absence of this element, the splicing intermediate accumulates and is not further processed to generate the cryptic exon. Our results indicate that Alu-derived sequences can provide intronic splicing regulatory elements that facilitate pre-mRNA processing and potentially affect the severity of disease-causing splicing mutations.     

Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans

RNA sequence elements involved in the regulation of pre-mRNA splicing have previously been identified in vertebrate genomes by computational methods. Here, we apply such approaches to predict splicing regulatory elements in Drosophila melanogaster and compare them with elements previously found in the human, mouse, and pufferfish genomes. We identified 99 putative exonic splicing enhancers (ESEs) and 231 putative intronic splicing enhancers (ISEs) enriched near weak 5' and 3' splice sites of constitutively spliced introns, distinguishing between those found near short and long introns. We found that a significant proportion (58%) of fly enhancer sequences were previously reported in at least one of the vertebrates. Furthermore, 20% of putative fly ESEs were previously identified as ESEs in human, mouse, and pufferfish; while only two fly ISEs, CTCTCT and TTATAA, were identified as ISEs in all three vertebrate species. Several putative enhancer sequences are similar to characterized binding-site motifs for Drosophila and mammalian splicing regulators. To provide additional evidence for the function of putative ISEs, we separately identified 298 intronic hexamers significantly enriched within sequences phylogenetically conserved among 15 insect species. We found that 73 putative ISEs were among those enriched in conserved regions of the D. melanogaster genome. The functions of nine enhancer sequences were verified in a heterologous splicing reporter, demonstrating that these sequences are sufficient to enhance splicing in vivo. Taken together, these data identify a set of predicted positive-acting splicing regulatory motifs in the Drosophila genome and reveal regulatory sequences that are present in distant metazoan genomes.     

Functional properties and evolutionary splicing constraints on a composite exonic regulatory element of splicing in CFTR exon 12

In general, splicing regulatory elements are defined as Enhancers or Silencers depending on their positive or negative effect upon exon inclusion. Often, these sequences are usually present separate from each other in exonic/intronic sequences. The Composite Exonic Splicing Regulatory Elements (CERES) represent an extreme physical overlap of enhancer/silencer activity. As a result, when CERES elements are mutated the consequences on the splicing process are difficult to predict. Here, we show that the functional activity of the CERES2 sequence in CFTR exon 12 is regulated by the binding, in very close proximity to each other, of several SR and hnRNP proteins. Moreover, our results show that practically the entire exon 12 sequence context participate in its definition. The consequences of this situation can be observed at the evolutionary level by comparing changes in conservation of different splicing elements in different species. In conclusion, our study highlights how it is increasingly difficult to define many exonic sequences by simply breaking them down in isolated enhancer/silencer or even neutral elements. The real picture is close to one of continuous competition between positive and negative factors where affinity for the target sequences and other dynamic factors decide the inclusion or exclusion of the exon.     

Mutations in U5 snRNA loop 1 influence the splicing of different genes in vivo

The U5 snRNA loop 1 is characterized by the conserved sequence G₁C₂C₃U₄U₅U₆Y₇A₈Y₉ and is essential for the alignment of exons during the second step of pre-mRNA splicing in Saccharo myces cerevisiae. Despite this sequence conservation the size, rather than sequence, of loop 1 is critical for exon alignment in vitro. To determine the in vivo requirements for U5 loop 1 a library of loop 1 sequences was transformed into a yeast strain where the endogenous U5 gene was deleted. Comparison of viable mutations in loop 1 revealed that position 6 was invariant and positions 5 and 7 displayed some sequence conservation. These data indicate positions 5, 6 and 7 in loop 1 are important for U5 function in vivo. A screen for mutations that suppress the temperature-sensitive phenotype of three loop 1 mutants produced eight intragenic suppressors all containing alterations in loop 1. Further analysis of these temperature-sensitive mutants revealed that each displayed distinct cell cycle arrest phenotypes and pre-mRNA splicing inhibition patterns. The cell cycle arrest is likely attributed to inefficient splicing of α-tubulin pre-mRNA in one mutant and actin pre-mRNA in another. These results suggest that various mutations in loop 1 may affect the splicing of different pre-mRNAs in vivo.     

Inference of splicing regulatory activities by sequence neighborhood analysis

Sequence-specific recognition of nucleic-acid motifs is critical to many cellular processes. We have developed a new and general method called Neighborhood Inference (NI) that predicts sequences with activity in regulating a biochemical process based on the local density of known sites in sequence space. Applied to the problem of RNA splicing regulation, NI was used to predict hundreds of new exonic splicing enhancer (ESE) and silencer (ESS) hexanucleotides from known human ESEs and ESSs. These predictions were supported by cross-validation analysis, by analysis of published splicing regulatory activity data, by sequence-conservation analysis, and by measurement of the splicing regulatory activity of 24 novel predicted ESEs, ESSs, and neutral sequences using an in vivo splicing reporter assay. These results demonstrate the ability of NI to accurately predict splicing regulatory activity and show that the scope of exonic splicing regulatory elements is substantially larger than previously anticipated. Analysis of orthologous exons in four mammals showed that the NI score of ESEs, a measure of function, is much more highly conserved above background than ESE primary sequence. This observation indicates a high degree of selection for ESE activity in mammalian exons, with surprisingly frequent interchangeability between ESE sequences.     

Co-regulation of alternative splicing by diverse splicing factors in Caenorhabditis elegans

Regulation of alternative splicing is controlled by pre-mRNA sequences (cis-elements) and trans-acting protein factors that bind them. The combinatorial interactions of multiple protein factors with the cis-elements surrounding a given alternative splicing event lead to an integrated splicing decision. The mechanism of multifactorial splicing regulation is poorly understood. Using a splicing-sensitive DNA microarray, we assayed 352 Caenorhabditis elegans alternative cassette exons for changes in embryonic splicing patterns between wild-type and 12 different strains carrying mutations in a splicing factor. We identified many alternative splicing events that are regulated by multiple splicing factors. Many splicing factors have the ability to behave as splicing repressors for some alternative cassette exons and as splicing activators for others. Unexpectedly, we found that the ability of a given alternative splicing factor to behave as an enhancer or repressor of a specific splicing event can change during development. Our observations that splicing factors can change their effects on a substrate during development support a model in which combinatorial effects of multiple factors, both constitutive and developmentally regulated ones, contribute to the overall splicing decision.     

Finding signals that regulate alternative splicing in the post-genomic era

Alternative splicing of pre-mRNAs is central to the generation of diversity from the relatively small number of genes in metazoan genomes. Auxiliary cis elements and trans-acting factors are required for the recognition of constitutive and alternatively spliced exons and their inclusion in pre-mRNA. Here, we discuss the regulatory elements that direct alternative splicing and how genome-wide analyses can aid in their identification.