Searching for splicing motifs

Imaging Intron Evolution

Intron evolution may be readily imaged through the combined use of the “dot plot” function of the NCBI BLAST, aligning two sequences at a time, and the Vertebrate “Multiz” alignment and conservation tool of the UCSC Genome Browser. With the NCBI BLAST, an ideal alignment of two highly conserved sequences generates a diagonal straight line in the plot from the lower left corner to the upper right corner. Gaps in this line correspond to non-conserved sections. In addition, the dot plot of the alignment of a sequence with the same sequence after the removal of the Transposable Elements (TEs) can be observed along the diagonal gaps that correspond to the sites of TE insertion. The UCSC Genome Browser can graph, along the entire sequence of a single gene, the level of overall conservation in vertebrates. This level can be compared with the conservation level of the gene in one or more selected vertebrate species. As an example, we show the graphic analysis of the intron conservation in two genes: the mitochondrial solute carrier 21 (SLC25A21) and the growth hormone receptor (GHR), whose coding sequences are conserved through vertebrates, while their introns show dramatic changes in nucleotide composition and even length. In the SLC25A21, a few short but significant nucleotide sequences are conserved in zebrafish, Xenopus and humans, and the rate of conservation steadily increases from chicken/human to mouse/human alignments. In the GHR, a less conserved gene, the earlier indication of intron conservation is a small signal in chicken/human alignment. The UCSC tool may simultaneously display the conservation level of a gene in different vertebrates, with reference to the level of overall conservation in Vertebrates. It is shown that, at least in SLC25A21, the sites of higher conservation are not always coincident in chicken and zebrafish nor are the sites of higher vertebrate conservation.

OUP accepted manuscript

Alternative splicing (AS) and transcription elongation are vital biological processes, and their dysregulation causes multiple diseases, including tumors. However, the coregulatory mechanism of AS and transcription elongation in tumors remains unclear. This study demonstrates a novel AS pattern of tight junction protein 1 (ZO1) regulated by the RNA polymerase II elongation rate in colorectal cancer (CRC). Glioma tumor suppressor candidate region gene 1 (GLTSCR1) decreases the transcription elongation rate of ZO1 to provide a time window for binding of the splicing factor HuR to the specific motif in intron 22 of ZO1 and spliceosome recognition of the weak 3′ and 5′ splice sites in exon 23 to promote exon 23 inclusion. Since exon 23 inclusion in ZO1 suppresses migration and invasion of CRC cells, our findings suggest a novel potential therapeutic target for CRC.

SPLICELECT™: an adaptable cell surface display technology based on alternative splicing allowing the qualitative and quantitative prediction of secreted product at a single-cell level

We describe a mammalian expression construct (SPLICELECT™) that allows the redirection of a proportion of a secreted protein onto the cell surface using alternative splicing: whereas the majority of the RNA is spliced into a transcript encoding a secreted protein, a weak splice donor site yields a secondary transcript encoding, in addition, a C-terminal transmembrane domain. The different sequence elements can be modified in order to modulate the level of cell surface display and of secretion in an independent manner. In this work, we demonstrated that the cell surface display of stable cell lines is correlated with the level of the secreted protein of interest, but also with the level of heterodimerization in the case of a bispecific antibody. It was also shown that this construct may be useful for rapid screening of multiple antibody candidates in binding assays following transient transfection. Thus, the correlation of product quantity and quality of the secreted and of membrane-displayed product in combination with the flexibility of the construct with regards to cell surface display/secretion levels make SPLICELECT™ a valuable tool with many potential applications, not limited to industrial cell line development or antibody engineering.

ESI mutagenesis: a one-step method for introducing mutations into bacterial artificial chromosomes

A simple and efficient recombineering-based method for introducing point mutations into bacterial artificial chromosomes using an artificial intron cassette.

Bacterial artificial chromosome (BAC)–based transgenes have emerged as a powerful tool for controlled and conditional interrogation of protein function in higher eukaryotes. Although homologous recombination-based recombineering methods have streamlined the efficient integration of protein tags onto BAC transgenes, generating precise point mutations has remained less efficient and time-consuming. Here, we present a simplified method for inserting point mutations into BAC transgenes requiring a single recombineering step followed by antibiotic selection. This technique, which we call exogenous/synthetic intronization (ESI) mutagenesis, relies on co-integration of a mutation of interest along with a selectable marker gene, the latter of which is harboured in an artificial intron adjacent to the mutation site. Cell lines generated from ESI-mutated BACs express the transgenes equivalently to the endogenous gene, and all cells efficiently splice out the synthetic intron. Thus, ESI mutagenesis provides a robust and effective single-step method with high precision and high efficiency for mutating BAC transgenes.

Structural disruption of exonic stem-loops immediately upstream of the intron regulates mammalian splicing

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Large-scale comparative evaluation of user-friendly tools for predicting variant-induced alterations of splicing regulatory elements

Discriminating which nucleotide variants cause disease or contribute to phenotypic traits remains a major challenge in human genetics. In theory, any intragenic variant can potentially affect RNA splicing by altering splicing regulatory elements (SREs). However, these alterations are often ignored mainly because pioneer SRE predictors have proved inefficient. Here, we report the first large-scale comparative evaluation of four user-friendly SRE-dedicated algorithms (QUEPASA, HEXplorer, SPANR, and HAL) tested both as standalone tools and in multiple combined ways based on two independent benchmark datasets adding up to >1,300 exonic variants studied at the messenger RNA level and mapping to 89 different disease-causing genes. These methods display good predictive power, based on decision thresholds derived from the receiver operating characteristics curve analyses, with QUEPASA and HAL having the best accuracies either as standalone or in combination. Still, overall there was a tight race between the four predictors, suggesting that all methods may be of use. Additionally, QUEPASA and HEXplorer may be beneficial as well for predicting variant-induced creation of pseudoexons deep within introns. Our study highlights the potential of SRE predictors as filtering tools for identifying disease-causing candidates among the plethora of variants detected by high-throughput DNA sequencing and provides guidance for their use in genomic medicine settings.

Identification of novel nucleotide sequence variations in an extended region of the bovine leptin gene (LEP) across a variety of cattle breeds from New Zealand and Nigeria

Leptin is mainly secreted by white adipose tissue in animals. Leptin acts by stimulating or inhibiting the release of a neurotransmitter, which eventually results in a decrease in food/feed intake and an increase in energy expenditure. In this investigation, the polymerase chain reaction (PCR) coupled with single-strand conformation polymorphism (SSCP) analysis was used to reveal nucleotide sequence variations in bovine leptin gene (LEP) in 338 cattle of a variety of breeds farmed in New Zealand (NZ) and Nigeria. These included NZ Hereford, Angus, Shorthorn, and crossbred Holstein-Friesian  ×  Jersey cattle and the Nigerian Sokoto Gudali, Red Bororo, White Fulani, and crossbred Holstein-Friesian  ×  White Fulani cattle. Sequence analysis of three regions of bovine LEP that encompassed selected coding and non-coding regions, revealed a total of 12 nucleotide sequence variations (six in exons and six in introns). Of these, three are reported here for the first time, whereas nine have been previously described. Some of the variations identified were common in both the NZ and Nigerian cattle breeds, while others were peculiar to particular breeds from a specific region. The sharing of common variants across different breeds irrespective of geography may indicate an evolutionary relationship, just as the differences within a breed might be attributable to either selective pressure for specific traits or random genetic drift. The detection of both new and previously documented variations in bovine LEP suggests that the gene is highly variable.

5' and 3' splicing signals evolution in vertebrates: Analysis in a conserved gene family

The mitochondrial solute carrier genes (SLC25) are highly conserved during vertebrate evolution. In most SLC25 genes of zebrafish, chicken, mouse, and human, the introns are located at exactly superimposable positions. In these topographically corresponding introns we studied the composition of the initial and terminal hexanucleotides (5'ss and 3'ss) which are instrumental in splicing signaling, focusing on the evolutionary conservation/mutation dynamics of these genetically related sequences. At each position, the per cent conservation of zebrafish individual nucleotides in chicken, mouse and human is proportional to their percent frequency in zebrafish; furthermore, nucleotide mutations are biased in favor of the more represented nucleotides, thus compensating for those highly represented zebrafish nucleotides which have not been conserved. As a result of these evolutionary dynamics, the general nucleotide composition at each position has remained relatively conserved throughout vertebrates. At 5'ss, following the canonical GT, A and G are largely prevailing at position +3, A at +4 and G at +5 (GT[A/G]AGx). At 3'ss, T and C are largely prevailing at positions -6, -5 and -3, preceding the canonical intron terminal AG ([C/T] [C/T]x[C/T]AG). However, the actual composition of the tetranucleotides at 5' and 3' often does not conform to the above scheme. At 5'ss the more canonical sequence is completely expressed in 63% of cases and partially (2 or 1 matches) in 37 % of cases. At 3'ss the more canonical sequence is completely expressed in 71 % of cases and partially (2 or 1 matches) in 29 % of cases. The nucleotide conservation loss (nucleotide mutation) is higher in the evolution from fish to the last common ancestor of birds and mammals (58 %), then diminishes in the successive evolution steps up to the mammalian common ancestor (10 %), and becomes still lower at the divergence of rodents and primates (5 %).

Genetic variations and haplotypic diversity in the Myostatin gene of New Zealand cattle breeds

Myostatin (MSTN) is a circulating factor that is secreted by muscle cells, and that acts upon those cells to inhibit the proliferation of muscle fibres during pre-natal muscle growth. The Polymerase Chain Reaction (PCR) coupled with Single Strand Conformational Polymorphism (SSCP) analysis, was used to reveal variation in the bovine MSTN gene (MSTN) in 722 cattle from a variety of breeds farmed in New Zealand (NZ). These included Hereford, Angus, Charolais, Simmental, Red Poll, South Devon, Shorthorn, Murray Grey, cross-bred Holstein-Friesian × Jersey cattle, and other composite breeds of cattle. Sequence analysis of five regions of MSTN that encompassed coding and non-coding regions of the gene, revealed a total of twelve single-nucleotide substitutions (7 in intron 1 and 5 in a region spanning the intron 2 - exon 3 boundary), and a single nucleotide deletion. Of these 12 substitutions, five are reported here for the first time, whereas seven have been previously described. The deletion c.748-78del, was located in the intron 2 - exon 3 boundary region, and has been reported previously. No nucleotide variation was identified in exons 1, 2 and 3. A total of 18 extended haplotypes were resolved spanning two variable regions (intron 1 and the intron 2 - exon 3 boundary), some of which were common across the breeds, while others were peculiar to particular breeds. The genetic variations identified provide insight into the conserved and polymorphic nature of the coding and non-coding sequences of bovine MSTN respectively, and thus provides a baseline for further study into how variation in the gene might affect growth and carcass traits in NZ cattle.

Case-control study on the association between the GATA2 gene and premature myocardial infarction in the Iranian population

In recent decades, due to the high prevalence of coronary artery disease (CAD) and myocardial infarction (MI), numerous studies have attempted to elucidate genetic contributing factors in these complex disorders. A very interesting gene in this regard is GATA-binding protein 2 (GATA2), an important regulator of various gene expressions in vascular endothelial cells. Accordingly, the association of different GATA2 polymorphisms with CAD and MI has already been evaluated. Rs2713604 is a genetic marker whose association with CAD has not been reproduced in previous studies. Considering the importance of replicating the initial association, the present case-control study aimed to examine the association of this intronic variant with premature MI in a sample of the Iranian population. In this study, 193 participants from Jahrom Hospital (Jahrom, Iran) were consecutively recruited during a 1.5-year period, and, following blood sampling, genomic DNA was extracted. We then proceeded to genotype rs2713604 using the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method and statistically analyzed the data. After adjustment for hyperlipidemia, hypertension, and type 2 diabetes mellitus, the results of the multivariate regression analysis showed no significant association between rs2713604 and premature MI. Interestingly, the risk allele (A-allele) of rs2713604 displayed a slightly higher frequency among controls compared to cases.

Novel diagnostic tool for prediction of variant spliceogenicity derived from a set of 395 combined in silico/in vitro studies: an international collaborative effort

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5′ and 3′ consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

Serine/Arginine-rich protein family of splicing regulators: New approaches to study splice isoform functions

Serine/arginine-rich (SR) proteins are conserved RNA-binding proteins that play major roles in RNA metabolism. They function as molecular adaptors, facilitate spliceosome assembly and modulate constitutive and alternative splicing of pre-mRNAs. Pre-mRNAs encoding SR proteins and many other proteins involved in stress responses are extensively alternatively spliced in response to diverse stresses. Hence, it is proposed that stress-induced changes in splice isoforms contribute to the adaptation of plants to stress responses. However, functions of most SR genes and their splice isoforms in stress responses are not known. Lack of easy and robust tools hindered the progress in this area. Emerging technologies such as CRISPR/Cas9 will facilitate studies of SR function by enabling the generation of single and multiple knock-out mutants of SR subfamily members. Moreover, CRISPR/Cas13 allows targeted manipulation of splice isoforms from SR and other genes in a constitutive or tissue-specific manner to evaluate functions of individual splice variants. Identification of the in vivo targets of SR proteins and their splice variants using the recently developed TRIBE (Targets of RNA-binding proteins Identified By Editing) and other methods will help unravel their mode of action and splicing regulatory elements under various conditions. These new approaches are expected to provide significant new insights into the roles of SRs and splice isoforms in plants adaptation to diverse stresses.

Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

Detection of Viral RNA Splicing in Diagnostic Virology

This chapter is the first one to introduce the detection of viral RNA splicing as a new tool for clinical diagnosis of virus infections. These include various infections caused by influenza viruses, human immunodeficiency viruses (HIV), human T-cell leukemia viruses (HTLV), Torque teno viruses (TTV), parvoviruses, adenoviruses, hepatitis B virus, polyomaviruses, herpesviruses, and papillomaviruses. Detection of viral RNA splicing for active viral gene expression in a clinical sample is a nucleic acid-based detection. The interpretation of the detected viral RNA splicing results is straightforward without concern for carry-over DNA contamination, because the spliced RNA is smaller than its corresponding DNA template.

Although many methods can be used, a simple method to detect viral RNA splicing is reverse transcription-polymerase chain reaction (RT-PCR). In principle, the detection of spliced RNA transcripts by RT-PCR depends on amplicon selection and primer design. The most common approach is the amplification over the intron regions by a set of primers in flanking exons. A larger product than the predicted size of smaller, spliced RNA is in general an unspliced RNA or contaminating viral genomic DNA. A spliced mRNA always gives a smaller RT-PCR product than its unspliced RNA due to removal of intron sequences by RNA splicing. The contaminating viral DNA can be determined by a minus RT amplification (PCR). Alternatively, specific amplification of a spliced RNA can be obtained by using an exon-exon junction primer because the sequence at exon-exon junction is not present in the unspliced RNA nor in viral genomic DNA.

Histone and RNA-binding protein interaction creates crosstalk network for regulation of alternative splicing

Alternative splicing is an essential process in eukaryotes, as it increases the complexity of gene expression by generating multiple proteins from a single pre-mRNA. However, information on the regulatory mechanisms for alternative splicing is lacking, because splicing occurs over a short period via the transient interactions of proteins within functional complexes of the spliceosome. Here, we investigated in detail the molecular mechanisms connecting alternative splicing with epigenetic mechanisms. We identified interactions between histone proteins and splicing factors such as Rbfox2, Rbfox3, and splicing factor proline and glutamine rich protein (SFPQ) by in vivo crosslinking and immunoprecipitation. Furthermore, we confirmed that splicing factors were bound to specific modified residues of histone proteins. Additionally, changes in histone methylation due to histone methyltransferase inhibitor treatment notably affected alternative splicing in selected genes. Therefore, we suggested that there may be crosstalk mechanisms connecting histone modifications and RNA-binding proteins that increase the local concentration of RNA-binding proteins in alternative exon loci of nucleosomes by binding specific modified histone proteins, leading to alternative splicing. This crosstalk mechanism may play a major role in epigenetic processes such as histone modification and the regulation of alternative splicing.

Deep intronic variation in splicing regulatory element of the ERCC8 gene associated with severe but long-term survival Cockayne syndrome

Cockayne syndrome is an autosomal recessive multisystem disorder characterized by intellectual disability, microcephaly, severe growth failure, sensory impairment, peripheral neuropathy, and cutaneous sensitivity. This rare disease is linked to disease-causing variations in the ERCC6 (CSB) and ERCC8 (CSA) genes. Various degrees of severity have been described according to age at onset and survival, without any clear genotype-phenotype correlation. All types of nucleotide changes have been observed in CS genes, including splice variations mainly affecting the splice site consensus sequences. We report here the case of two brothers from a consanguineous family presenting a severe but long-term survival phenotype of Cockayne syndrome. We identified in the patients a homozygous deep intronic nucleotide variation causing the insertion of a cryptic exon in the ERCC8 (CSA) transcript, by modifying intronic regulatory elements important for exon definition. The pathogenesis of the nucleotide variant NG_009289.1(NM_000082.3):c.173+1119G>C was validated in vitro with a reporter minigene system. To our knowledge, these are the first Cockayne patients described with this kind of disease-causing variation, though molecular mechanism underlying early onset symptoms and unexpected slow raise of progression of the disease remain to be elucidated.

The SR protein B52/SRp55 regulates splicing of the period thermosensitive intron and mid-day siesta in Drosophila

Similar to many diurnal animals, Drosophila melanogaster exhibits a mid-day siesta that is more robust as temperature increases, an adaptive response that aims to minimize the deleterious effects from exposure to heat. This temperature-dependent plasticity in mid-day sleep levels is partly based on the thermal sensitive splicing of an intron in the 3′ untranslated region (UTR) of the circadian clock gene termed period (per). In this study, we evaluated a possible role for the serine/arginine-rich (SR) splicing factors in the regulation of dmpi8 splicing efficiency and mid-day siesta. Using a Drosophila cell culture assay we show that B52/SRp55 increases dmpi8 splicing efficiency, whereas other SR proteins have little to no effect. The magnitude of the stimulatory effect of B52 on dmpi8 splicing efficiency is modulated by natural variation in single nucleotide polymorphisms (SNPs) in the per 3′ UTR that correlate with B52 binding levels. Down-regulating B52 expression in clock neurons increases mid-day siesta and reduces dmpi8 splicing efficiency. Our results establish a novel role for SR proteins in sleep and suggest that polymorphisms in the per 3′ UTR contribute to natural variation in sleep behavior by modulating the binding efficiencies of SR proteins.

Exploring the relationship between intron retention and chromatin accessibility in plants

Background

Intron retention (IR) is the most prevalent form of alternative splicing in plants. IR, like other forms of alternative splicing, has an important role in increasing gene product diversity and regulating transcript functionality. Splicing is known to occur co-transcriptionally and is influenced by the speed of transcription which in turn, is affected by chromatin structure. It follows that chromatin structure may have an important role in the regulation of splicing, and there is preliminary evidence in metazoans to suggest that this is indeed the case; however, nothing is known about the role of chromatin structure in regulating IR in plants. DNase I-seq is a useful experimental tool for genome-wide interrogation of chromatin accessibility, providing information on regions of chromatin with very high likelihood of cleavage by the enzyme DNase I, known as DNase I Hypersensitive Sites (DHSs). While it is well-established that promoter regions are highly accessible and are over-represented with DHSs, not much is known about DHSs in the bodies of genes, and their relationship to splicing in general, and IR in particular.

Results

In this study we use publicly available DNase I-seq data in arabidopsis and rice to investigate the relationship between IR and chromatin structure. We find that IR events are highly enriched in DHSs in both species. This implies that chromatin is more open in retained introns, which is consistent with a kinetic model of the process whereby higher speeds of transcription in those regions give less time for the spliceosomal machinery to recognize and splice out those introns co-transcriptionally. The more open chromatin in IR can also be the result of regulation mediated by DNA-binding proteins. To test this, we performed an exhaustive search for footprints left by DNA-binding proteins that are associated with IR. We identified several hundred short sequence elements that exhibit footprints in their DNase I-seq coverage, the telltale sign for binding events of a regulatory protein, protecting its binding site from cleavage by DNase I. A highly significant fraction of those sequence elements are conserved between arabidopsis and rice, a strong indication of their functional importance.

Conclusions

In this study we have established an association between IR and chromatin accessibility, and presented a mechanistic hypothesis that explains the observed association from the perspective of the co-transcriptional nature of splicing. Furthermore, we identified conserved sequence elements for DNA-binding proteins that affect splicing.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4393-z) contains supplementary material, which is available to authorized users.

Saturation mutagenesis reveals manifold determinants of exon definition

To illuminate the extent and roles of exonic sequences in the splicing of human RNA transcripts, we conducted saturation mutagenesis of a 51-nt internal exon in a three-exon minigene. All possible single and tandem dinucleotide substitutions were surveyed. Using high-throughput genetics, 5560 minigene molecules were assayed for splicing in human HEK293 cells. Up to 70% of mutations produced substantial (greater than twofold) phenotypes of either increased or decreased splicing. Of all predicted secondary structural elements, only a single 15-nt stem–loop showed a strong correlation with splicing, acting negatively. The in vitro formation of exon-protein complexes between the mutant molecules and proteins associated with spliceosome formation (U2AF35, U2AF65, U1A, and U1-70K) correlated with splicing efficiencies, suggesting exon definition as the step affected by most mutations. The measured relative binding affinities of dozens of human RNA binding protein domains as reported in the CISBP-RNA database were found to correlate either positively or negatively with splicing efficiency, more than could fit on the 51-nt test exon simultaneously. The large number of these functional protein binding correlations point to a dynamic and heterogeneous population of pre-mRNA molecules, each responding to a particular collection of binding proteins.

The combinatorial control of alternative splicing in C. elegans

Normal development requires the right splice variants to be made in the right tissues at the right time. The core splicing machinery is engaged in all splicing events, but which precise splice variant is made requires the choice between alternative splice sites—for this to occur, a set of splicing factors (SFs) must recognize and bind to short RNA motifs in the pre-mRNA. In C. elegans, there is known to be extensive variation in splicing patterns across development, but little is known about the targets of each SF or how multiple SFs combine to regulate splicing. Here we combine RNA-seq with in vitro binding assays to study how 4 different C. elegans SFs, ASD-1, FOX-1, MEC-8, and EXC-7, regulate splicing. The 4 SFs chosen all have well-characterised biology and well-studied loss-of-function genetic alleles, and all contain RRM domains. Intriguingly, while the SFs we examined have varied roles in C. elegans development, they show an unexpectedly high overlap in their targets. We also find that binding sites for these SFs occur on the same pre-mRNAs more frequently than expected suggesting extensive combinatorial control of splicing. We confirm that regulation of splicing by multiple SFs is often combinatorial and show that this is functionally significant. We also find that SFs appear to combine to affect splicing in two modes—they either bind in close proximity within the same intron or they appear to bind to separate regions of the intron in a conserved order. Finally, we find that the genes whose splicing are regulated by multiple SFs are highly enriched for genes involved in the cytoskeleton and in ion channels that are key for neurotransmission. Together, this shows that specific classes of genes have complex combinatorial regulation of splicing and that this combinatorial regulation is critical for normal development to occur.

Alternative splicing (AS) is a highly regulated process that is crucial for normal development. It requires the core splicing machinery, but the specific choice of splice site during AS is controlled by splicing factors (SFs) such as ELAV or RBFOX proteins that bind to specific sequences in pre-mRNAs to regulate usage of different splice sites. AS varies across the C. elegans life cycle and here we study how diverse SFs combine to regulate AS during C. elegans development. We selected 4 RRM-containing SFs that are all well studied and that have well-characterised loss-of-function genetic alleles. We find that these SFs regulate many of the same targets, and that combinatorial interactions between these SFs affect both individual splicing events and organism-level phenotypes including specific effects on the neuromuscular system. We further show that SFs combine to regulate splicing of an individual pre-mRNA in two distinct modes—either by binding in close proximity or by binding in a defined order on the pre-mRNA. Finally, we find that the genes whose splicing are most likely to be regulated by multiple SFs are genes that are required for the proper function of the neuromuscular system. These genes are also most likely to have changing AS patterns across development, suggesting that their splicing regulation is highly complex and developmentally regulated. Taken together, our data show that the precise splice variant expressed at any point in development is often the outcome of regulation by multiple SFs.

Succession of splicing regulatory elements determines cryptic 5΄ss functionality

A critical step in exon definition is the recognition of a proper splice donor (5΄ss) by the 5’ end of U1 snRNA. In the selection of appropriate 5΄ss, cis-acting splicing regulatory elements (SREs) are indispensable. As a model for 5΄ss recognition, we investigated cryptic 5΄ss selection within the human fibrinogen Bβ-chain gene (FGB) exon 7, where we identified several exonic SREs that simultaneously acted on up- and downstream cryptic 5΄ss. In the FGB exon 7 model system, 5΄ss selection iteratively proceeded along an alternating sequence of U1 snRNA binding sites and interleaved SREs which in principle supported different 3’ exon ends. Like in a relay race, SREs either suppressed a potential 5΄ss and passed the splicing baton on or splicing actually occurred. From RNA-Seq data, we systematically selected 19 genes containing exons with silent U1 snRNA binding sites competing with nearby highly used 5΄ss. Extensive SRE analysis by different algorithms found authentic 5΄ss significantly more supported by SREs than silent U1 snRNA binding sites, indicating that our concept may permit generalization to a model for 5΄ss selection and 3’ exon end definition.

Improving genetic diagnosis in Mendelian disease with transcriptome sequencing

Exome and whole-genome sequencing are becoming increasingly routine approaches in Mendelian disease diagnosis. Despite their success, the current diagnostic rate for genomic analyses across a variety of rare diseases is approximately 25 to 50%. We explore the utility of transcriptome sequencing [RNA sequencing (RNA-seq)] as a complementary diagnostic tool in a cohort of 50 patients with genetically undiagnosed rare muscle disorders. We describe an integrated approach to analyze patient muscle RNA-seq, leveraging an analysis framework focused on the detection of transcript-level changes that are unique to the patient compared to more than 180 control skeletal muscle samples. We demonstrate the power of RNA-seq to validate candidate splice-disrupting mutations and to identify splice-altering variants in both exonic and deep intronic regions, yielding an overall diagnosis rate of 35%. We also report the discovery of a highly recurrent de novo intronic mutation in COL6A1 that results in a dominantly acting splice-gain event, disrupting the critical glycine repeat motif of the triple helical domain. We identify this pathogenic variant in a total of 27 genetically unsolved patients in an external collagen VI-like dystrophy cohort, thus explaining approximately 25% of patients clinically suggestive of having collagen VI dystrophy in whom prior genetic analysis is negative. Overall, this study represents a large systematic application of transcriptome sequencing to rare disease diagnosis and highlights its utility for the detection and interpretation of variants missed by current standard diagnostic approaches.

Depletion of Arabidopsis SC35 and SC35-like serine/arginine-rich proteins affects the transcription and splicing of a subset of genes

Serine/arginine-rich (SR) proteins are important splicing factors which play significant roles in spliceosome assembly and splicing regulation. However, little is known regarding their biological functions in plants. Here, we analyzed the phenotypes of mutants upon depleting different subfamilies of Arabidopsis SR proteins. We found that loss of the functions of SC35 and SC35-like (SCL) proteins cause pleiotropic changes in plant morphology and development, including serrated leaves, late flowering, shorter roots and abnormal silique phyllotaxy. Using RNA-seq, we found that SC35 and SCL proteins play roles in the pre-mRNA splicing. Motif analysis revealed that SC35 and SCL proteins preferentially bind to a specific RNA sequence containing the AGAAGA motif. In addition, the transcriptions of a subset of genes are affected by the deletion of SC35 and SCL proteins which interact with NRPB4, a specific subunit of RNA polymerase II. The splicing of FLOWERING LOCUS C (FLC) intron1 and transcription of FLC were significantly regulated by SC35 and SCL proteins to control Arabidopsis flowering. Therefore, our findings provide mechanistic insight into the functions of plant SC35 and SCL proteins in the regulation of splicing and transcription in a direct or indirect manner to maintain the proper expression of genes and development.

SR proteins were identified to be important splicing factors. This work generated mutants of different subfamilies of the classic Arabidopsis SR proteins. Genetic analysis revealed that loss of the function of SC35/SCL proteins influences the plant development. This study revealed SC35/SCL proteins regulate alternative splicing, preferentially bind a specific RNA motif, interact with NRPB4, and affect the transcription of a subset of genes. This study further revealed that SC35/SCL proteins control flowering by regulating the splicing and transcription of FLC. These results shed light on the functions of SR proteins in plants.

Conservation/Mutation in the Splice Sites of Mitochondrial Solute Carrier Genes of Vertebrates

The "canonical" introns begin by the dinucleotide GT and end by the dinucleotide AG. GT, together with a few downstream nucleotides, and AG, with a few of the immediately preceding nucleotides, are thought to be the strongest splicing signals (5'ss and 3'ss, respectively). We examined the composition of the intronic initial and terminal hexanucleotides of the mitochondrial solute carrier genes (SLC25A's) of zebrafish, chicken, mouse, and human. These genes are orthologous and we selected the transcripts in which the arrangement of exons and introns was superimposable in the species considered. Both 5'ss and 3'ss were highly polymorphic, with 104 and 126 different configurations, respectively, in our sample. In the line of evolution from zebrafish to chicken, as well as in that from zebrafish to mammals, the average nucleotide conservation in the four variable nucleotides was about 50 % at 5' and 40 % at 3'. In the divergent evolution of mouse and human, the conservation was about 80 % at 5' and 70 % at 3'. Despite these changes, the splicing signals remain strong enough to operate at the same site. At both 5' and 3', the frequency of a nucleotide at a given position in the zebrafish sequence is positively correlated with its conservation in chicken and mammals, suggesting that selection continued to operate in birds and mammals along similar lines.

Defective control of pre-messenger RNA splicing in human disease

Examples of associations between human disease and defects in pre-messenger RNA splicing/alternative splicing are accumulating. Although many alterations are caused by mutations in splicing signals or regulatory sequence elements, recent studies have noted the disruptive impact of mutated generic spliceosome components and splicing regulatory proteins. This review highlights recent progress in our understanding of how the altered splicing function of RNA-binding proteins contributes to myelodysplastic syndromes, cancer, and neuropathologies.

mCarts: Genome-Wide Prediction of Clustered Sequence Motifs as Binding Sites for RNA-Binding Proteins

RNA-binding proteins (RBPs) are critical components of post-transcriptional gene expression regulation. However, their binding sites have until recently been difficult to determine due to the apparent low specificity of RBPs for their target transcripts and the lack of high-throughput assays for analyzing binding sites genome wide. Here we present a bioinformatics method for predicting RBP binding motif sites on a genome-wide scale that leverages motif conservation, RNA secondary structure, and the tendency of RBP binding sites to cluster together. A probabilistic model is learned from bona fide binding sites determined by CLIP and applied genome wide to generate high specificity binding site predictions.

Regulation of alternative splicing through coupling with transcription and chromatin structure

Alternative precursor messenger RNA (pre-mRNA) splicing plays a pivotal role in the flow of genetic information from DNA to proteins by expanding the coding capacity of genomes. Regulation of alternative splicing is as important as regulation of transcription to determine cell- and tissue-specific features, normal cell functioning, and responses of eukaryotic cells to external cues. Its importance is confirmed by the evolutionary conservation and diversification of alternative splicing and the fact that its deregulation causes hereditary disease and cancer. This review discusses the multiple layers of cotranscriptional regulation of alternative splicing in which chromatin structure, DNA methylation, histone marks, and nucleosome positioning play a fundamental role in providing a dynamic scaffold for interactions between the splicing and transcription machineries. We focus on evidence for how the kinetics of RNA polymerase II (RNAPII) elongation and the recruitment of splicing factors and adaptor proteins to chromatin components act in coordination to regulate alternative splicing.

A single-nucleotide exon found in Arabidopsis

The presence of introns in gene-coding regions is one of the most mysterious evolutionary inventions in eukaryotic organisms. It has been proposed that, although sequences involved in intron recognition and splicing are mainly located in introns, exonic sequences also contribute to intron splicing. The smallest constitutively spliced exon known so far has 6 nucleotides, and the smallest alternatively spliced exon has 3 nucleotides. Here we report that the Anaphase Promoting Complex subunit 11 (APC11) gene in Arabidopsis thaliana carries a constitutive single-nucleotide exon. In vivo transcription and translation assays performed using APC11-Green Fluorescence Protein (GFP) fusion constructs revealed that intron splicing surrounding the single-nucleotide exon is effective in both Arabidopsis and rice. This discovery warrants attention to genome annotations in the future.

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.

HuR antagonizes the effect of an intronic pyrimidine-rich sequence in regulating WT1 +/-KTS isoforms

WT1 + KTS and -KTS isoforms only differ in 3 amino acids in protein sequence but show significant functional difference. The +/-KTS isoforms were generated by alternative usage of 2 adjacent 5' splice sites at RNA level, however, how these 2 isoforms are regulated is still elusive. Here we report the identification of an intronic pyrimidine-rich sequence that is critical for the ratio of +/-KTS isoforms, deletion or partial replacement of the sequence led to full/significant shift to -KTS isoform. To identify trans-factors that can regulate +/-KTS isoforms via the binding to the element, we performed RNP assembly using in vitro transcribed RNA with or without the pyrimidine-rich sequence. Mass spectrometry analysis of purified RNPs showed that the element associated with many splicing factors. Co-transfection of these factors with WT1 reporter revealed that HuR promoted the production of -KTS isoform at the reporter level. RNA immuno-precipitation experiment indicated that HuR interacted with the pyrimidine-rich element in WT1 intron 9. We further presented evidence that transient or stable over-expression of HuR led to enhanced expression of endogenous -KTS isoform. Moreover, knockdown of HuR resulted in decreased expression of endogenous -KTS isoform in 293T, SW620, SNU-387 and AGS cell lines. Together, these data indicate that HuR binds to the pyrimidine-rich sequence and antagonize its effect in regulating WT1 +/-KTS isoforms.

CoSREM: a graph mining algorithm for the discovery of combinatorial splicing regulatory elements

BACKGROUND

Alternative splicing (AS) is a post-transcriptional regulatory mechanism for gene expression regulation. Splicing decisions are affected by the combinatorial behavior of different splicing factors that bind to multiple binding sites in exons and introns. These binding sites are called splicing regulatory elements (SREs). Here we develop CoSREM (Combinatorial SRE Miner), a graph mining algorithm to discover combinatorial SREs in human exons. Our model does not assume a fixed length of SREs and incorporates experimental evidence as well to increase accuracy. CoSREM is able to identify sets of SREs and is not limited to SRE pairs as are current approaches.

RESULTS

We identified 37 SRE sets that include both enhancer and silencer elements. We show that our results intersect with previous results, including some that are experimental. We also show that the SRE set GGGAGG and GAGGAC identified by CoSREM may play a role in exon skipping events in several tumor samples. We applied CoSREM to RNA-Seq data for multiple tissues to identify combinatorial SREs which may be responsible for exon inclusion or exclusion across tissues.

CONCLUSION

The new algorithm can identify different combinations of splicing enhancers and silencers without assuming a predefined size or limiting the algorithm to find only pairs of SREs. Our approach opens new directions to study SREs and the roles that AS may play in diseases and tissue specificity.

Identifying splicing regulatory elements with de Bruijn graphs

Splicing regulatory elements (SREs) are short, degenerate sequences on pre-mRNA molecules that enhance or inhibit the splicing process via the binding of splicing factors, proteins that regulate the functioning of the spliceosome. Existing methods for identifying SREs in a genome are either experimental or computational. Here, we propose a formalism based on de Bruijn graphs that combines genomic structure, word count enrichment analysis, and experimental evidence to identify SREs found in exons. In our approach, SREs are not restricted to a fixed length (i.e., k-mers, for a fixed k). As a result, we identify 2001 putative exonic enhancers and 3080 putative exonic silencers for human genes, with lengths varying from 6 to 15 nucleotides. Many of the predicted SREs overlap with experimentally verified binding sites. Our model provides a novel method to predict variable length putative regulatory elements computationally for further experimental investigation.

Identification and Validation of Evolutionarily Conserved Unusually Short Pre-mRNA Introns in the Human Genome

According to the length distribution of human introns, there is a large population of short introns with a threshold of 65 nucleotides (nt) and a peak at 85 nt. Using human genome and transcriptome databases, we investigated the introns shorter than 66 nt, termed ultra-short introns, the identities of which are scarcely known. Here, we provide for the first time a list of bona fide human ultra-short introns, which have never been characterized elsewhere. By conducting BLAST searches of the databases, we screened 22 introns (37–65 nt) with conserved lengths and sequences among closely related species. We then provide experimental and bioinformatic evidence for the splicing of 15 introns, of which 12 introns were remarkably G-rich and 9 introns contained completely inefficient splice sites and/or branch sites. These unorthodox characteristics of ultra-short introns suggest that there are unknown splicing mechanisms that differ from the well-established mechanism.

Coupling and coordination in gene expression processes with pre-mRNA splicing

RNA processing is a tightly regulated and highly complex pathway which includes transcription, splicing, editing, transportation, translation and degradation. It has been well-documented that splicing of RNA polymerase II medicated nascent transcripts occurs co-transcriptionally and is functionally coupled to other RNA processing. Recently, increasing experimental evidence indicated that pre-mRNA splicing influences RNA degradation and vice versa. In this review, we summarized the recent findings demonstrating the coupling of these two processes. In addition, we highlighted the importance of splicing in the production of intronic miRNA and circular RNAs, and hence the discovery of the novel mechanisms in the regulation of gene expression.

Splicing of designer exons informs a biophysical model for exon definition

Pre-mRNA molecules in humans contain mostly short internal exons flanked by longer introns. To explain the removal of such introns, exon recognition instead of intron recognition has been proposed. We studied this exon definition using designer exons (DEs) made up of three prototype modules of our own design: an exonic splicing enhancer (ESE), an exonic splicing silencer (ESS), and a Reference Sequence (R) predicted to be neither. Each DE was examined as the central exon in a three-exon minigene. DEs made of R modules showed a sharp size dependence, with exons shorter than 14 nt and longer than 174 nt splicing poorly. Changing the strengths of the splice sites improved longer exon splicing but worsened shorter exon splicing, effectively displacing the curve to the right. For the ESE we found, unexpectedly, that its enhancement efficiency was independent of its position within the exon. For the ESS we found a step-wise positional increase in its effects; it was most effective at the 3' end of the exon. To apply these results quantitatively, we developed a biophysical model for exon definition of internal exons undergoing cotranscriptional splicing. This model features commitment to inclusion before the downstream exon is synthesized and competition between skipping and inclusion fates afterward. Collision of both exon ends to form an exon definition complex was incorporated to account for the effect of size; ESE/ESS effects were modeled on the basis of stabilization/destabilization. This model accurately predicted the outcome of independent experiments on more complex DEs that combined ESEs and ESSs.

Genetic variation and RNA binding proteins: tools and techniques to detect functional polymorphisms

At its most fundamental level the goal of genetics is to connect genotype to phenotype. This question is asked at a basic level evaluating the role of genes and pathways in genetic model organism. Increasingly, this question is being asked in the clinic. Genomes of individuals and populations are being sequenced and compared. The challenge often comes at the stage of analysis. The variant positions are analyzed with the hope of understanding human disease. However after a genome or exome has been sequenced, the researcher is often deluged with hundreds of potentially relevant variations. Traditionally, amino-acid changing mutations were considered the tractable class of disease-causing mutations; however, mutations that disrupt noncoding elements are the subject of growing interest. These noncoding changes are a major avenue of disease (e.g., one in three hereditary disease alleles are predicted to affect splicing). Here, we review some current practices of medical genetics, the basic theory behind biochemical binding and functional assays, and then explore technical advances in how variations that alter RNA protein recognition events are detected and studied. These advances are advances in scale-high-throughput implementations of traditional biochemical assays that are feasible to perform in any molecular biology laboratory. This chapter utilizes a case study approach to illustrate some methods for analyzing polymorphisms. The first characterizes a functional intronic SNP that deletes a high affinity PTB site using traditional low-throughput biochemical and functional assays. From here we demonstrate the utility of high-throughput splicing and spliceosome assembly assays for screening large sets of SNPs and disease alleles for allelic differences in gene expression. Finally we perform three pilot drug screens with small molecules (G418, tetracycline, and valproic acid) that illustrate how compounds that rescue specific instances of differential pre-mRNA processing can be discovered.

Slight variations in the SC35 ESE sequence motif among human chromosomes: a computational approach

Gene expression is initiated by the binding of transcription factors to cis-regulatory modules such as enhancer elements binding to the Serine/Arginine proteins. Recently, we noticed an increased ability to identify the location as well as the motifs of enhancers using genome-wide information on spliceosomal factor occupancy, cofactor recruitment and chromatin modifications. In this study, we have undertaken a large-scale genomic analysis in an attempt to uncover if the exonic splicing enhancer motif binding to the SC35 and the SRp40 SR proteins is conserved among several groups of human genes. For the SRp40, the results showed that the ESE consensus is conserved among human genes. Concerning the SC35 SR protein, results showed an ESE motif conserved among human tissues and between different levels of muscular cell differentiation and within the same chromosome. However, this motif displays subtle discrepancies between genes localized in different chromosomes. These results emphasize the presence of different translational isoforms of the SFRS2 gene encoding for the SC35, or different post-translational protein maturations in different chromosomes, confirming that chromatin structure is another layer of gene regulation. These links between chromatin pattern and splicing give further mechanistic support to functional interconnections between splicing, transcription and chromatin structure, and raise the intriguing possibility of the existence of a memory for splicing patterns to be inherited through epigenetic modifications.

SeeSite: Characterizing Relationships between Splice Junctions and Splicing Enhancers

RNA splicing is a cellular process driven by the interaction between numerous regulatory sequences and binding sites, however, such interactions have been primarily explored by laboratory methods since computational tools largely ignore the relationship between different splicing elements. Current computational methods identify either splice sites or other regulatory sequences, such as enhancers and silencers. We present a novel approach for characterizing co-occurring relationships between splice site motifs and splicing enhancers. Our approach relies on an efficient algorithm for approximately solving Consensus Sequence with Outliers , an NP-complete string clustering problem. In particular, we give an algorithm for this problem that outputs near-optimal solutions in polynomial time. To our knowledge, this is the first formulation and computational attempt for detecting co-occurring sequence elements in RNA sequence data. Further, we demonstrate that SeeSite is capable of showing that certain ESEs are preferentially associated with weaker splice sites, and that there exists a co-occurrence relationship with splice site motifs.

NFAT1 and NFAT3 cooperate with HDAC4 during regulation of alternative splicing of PMCA isoforms in PC12 cells

Background

The bulk of human genes undergo alternative splicing (AS) upon response to physiological stimuli. AS is a great source of protein diversity and biological processes and is associated with the development of many diseases. Pheochromocytoma is a neuroendocrine tumor, characterized by an excessive Ca²⁺-dependent secretion of catecholamines. This underlines the importance of balanced control of calcium transport via regulation of gene expression pattern, including different calcium transport systems, such as plasma membrane Ca²⁺-ATPases (PMCAs), abundantly expressed in pheochromocytoma chromaffin cells (PC12 cells). PMCAs are encoded by four genes (Atp2b1, Atp2b2, Atp2b3, Atp2b4), whose transcript products undergo alternative splicing giving almost 30 variants.

Results

In this scientific report, we propose a novel mechanism of regulation of PMCA alternative splicing in PC12 cells through cooperation of the nuclear factor of activated T-cells (NFAT) and histone deacetylases (HDACs). Luciferase assays showed increased activity of NFAT in PC12 cells, which was associated with altered expression of PMCA. RT-PCR experiments suggested that inhibition of the transcriptional activity of NFAT might result in the rearrangement of PMCA splicing variants in PC12 cells. NFAT inhibition led to dominant expression of 2x/c, 3x/a and 4x/a PMCA variants, while in untreated cells the 2w,z/b, 3z,x/b,c,e,f, and 4x/b variants were found as well. Furthermore, chromatin immunoprecipitation experiments showed that NFAT1-HDAC4 or NFAT3-HDAC4 complexes might be involved in regulation of PMCA2x splicing variant generation.

Conclusions

We suggest that the influence of NFAT/HDAC on PMCA isoform composition might be important for altered dopamine secretion by PC12 cells.

Histone modifications are associated with transcript isoform diversity in normal and cancer cells

Mechanisms that generate transcript diversity are of fundamental importance in eukaryotes. Although a large fraction of human protein-coding genes and lincRNAs produce more than one mRNA isoform each, the regulation of this phenomenon is still incompletely understood. Much progress has been made in deciphering the role of sequence-specific features as well as DNA-and RNA-binding proteins in alternative splicing. Recently, however, several experimental studies of individual genes have revealed a direct involvement of epigenetic factors in alternative splicing and transcription initiation. While histone modifications are generally correlated with overall gene expression levels, it remains unclear how histone modification enrichment affects relative isoform abundance. Therefore, we sought to investigate the associations between histone modifications and transcript diversity levels measured by the rates of transcription start-site switching and alternative splicing on a genome-wide scale across protein-coding genes and lincRNAs. We found that the relationship between enrichment levels of epigenetic marks and transcription start-site switching is similar for protein-coding genes and lincRNAs. Furthermore, we found associations between splicing rates and enrichment levels of H2az, H3K4me1, H3K4me2, H3K4me3, H3K9ac, H3K9me3, H3K27ac, H3K27me3, H3K36me3, H3K79me2, and H4K20me, marks traditionally associated with enhancers, transcription initiation, transcriptional repression, and others. These patterns were observed in both normal and cancer cell lines. Additionally, we developed a novel computational method that identified 840 epigenetically regulated candidate genes and predicted transcription start-site switching and alternative exon splicing with up to 92% accuracy based on epigenetic patterning alone. Our results suggest that the epigenetic regulation of transcript isoform diversity may be a relatively common genome-wide phenomenon representing an avenue of deregulation in tumor development.

Traditionally, the regulation of gene expression was thought to be largely based on DNA and RNA sequence motifs. However, this dogma has recently been challenged as other factors, such as epigenetic patterning of the genome, have become better understood. Sparse but convincing experimental evidence suggests that the epigenetic background, in the form of histone modifications, acts as an additional layer of regulation determining how transcripts are processed. Here we developed a computational approach to investigate the genome-wide prevalence and the level of association between the enrichment of epigenetic marks and transcript diversity generated via alternative transcription start sites and splicing. We found that the role of epigenetic patterning in alternative transcription start-site switching is likely to be the same for all genes whereas the role of epigenetic patterns in splicing is likely gene-specific. Furthermore, we show that epigenetic data alone can be used to predict the inclusion pattern of an exon. These findings have significant implications for a better understanding of the regulation of transcript diversity in humans as well as the modifications arising during tumor development.

T cell activation regulates CD6 alternative splicing by transcription dynamics and SRSF1

The T cell-surface glycoprotein CD6 is a modulator of cellular responses and has been implicated in several autoimmune diseases such as multiple sclerosis, rheumatoid arthritis, and psoriasis. During Ag presentation, CD6 is targeted to the immunological synapse in a ligand binding-dependent manner, in which CD6 domain 3 directly contacts CD166, expressed on the APC. T cell activation results in the induction of CD6Δd3, an alternatively spliced isoform that lacks the ligand-binding domain and thus no longer localizes at the immunological synapse. In this study, we investigated the molecular mechanisms regulating the expression of CD6Δd3 upon human primary T cell activation. Using chromatin immunoprecipitation, we observed an increase in RNA polymerase II occupancy along the CD6 gene and augmented CD6 transcription. We showed that activation leads to transcription-related chromatin modifications, revealed by higher CD6 acetylation levels. Modulation of chromatin conformation using a histone deacetylase inhibitor that increases transcription rate causes an increase of exon 5 skipping. We further showed that the splicing factor SRSF1 binds to a regulatory element in CD6 intron 4, activating exon 5 splicing and promoting exon 5 inclusion. Concomitant with T cell activation-induced exon 5 skipping, we observed a downregulation of SRSF1. Using RNA immunoprecipitation, we showed that in activated T cells, SRSF1 recruitment to the CD6 transcript is impaired by increased chromatin acetylation levels. We propose that upon T cell activation, SRSF1 becomes limiting, and its function in CD6 exon 5 splicing is countered by an increase in CD6 transcription, dependent on chromatin acetylation.

The UPF1 RNA surveillance gene is commonly mutated in pancreatic adenosquamous carcinoma

Pancreatic adenosquamous carcinoma (ASC) is an enigmatic and aggressive tumor that has a worse prognosis and higher metastatic potential than its adenocarcinoma counterpart. Here we report that ASC tumors frequently harbor somatically acquired mutations in the UPF1 gene, which encodes the core component of the nonsense-mediated RNA decay (NMD) pathway. These tumor-specific mutations alter UPF1 RNA splicing and perturb NMD, leading to upregulated levels of NMD substrate mRNAs. UPF1 mutations are, to our knowledge, the first known unique molecular signatures of pancreatic ASC.

A view of pre-mRNA splicing from RNase R resistant RNAs

During pre-mRNA splicing, exons in the primary transcript are precisely connected to generate an mRNA. Intron lariat RNAs are formed as by-products of this process. In addition, some exonic circular RNAs (circRNAs) may also result from exon skipping as by-products. Lariat RNAs and circRNAs are both RNase R resistant RNAs. RNase R is a strong 3' to 5' exoribonuclease, which efficiently degrades linear RNAs, such as mRNAs and rRNAs; therefore, the circular parts of lariat RNAs and the circRNAs can be segregated from eukaryotic total RNAs by their RNase R resistance. Thus, RNase R resistant RNAs could provide unexplored splicing information not available from mRNAs. Analyses of these RNAs identified repeating splicing phenomena, such as re-splicing of mature mRNAs and nested splicing. Moreover, circRNA might function as microRNA sponges. There is an enormous variety of endogenous circRNAs, which are generally synthesized in cells and tissues.

Chronic cocaine-regulated epigenomic changes in mouse nucleus accumbens

BACKGROUND

Increasing evidence supports a role for altered gene expression in mediating the lasting effects of cocaine on the brain, and recent work has demonstrated the involvement of chromatin modifications in these alterations. However, all such studies to date have been restricted by their reliance on microarray technologies that have intrinsic limitations.

RESULTS

We use next generation sequencing methods, RNA-seq and ChIP-seq for RNA polymerase II and several histone methylation marks, to obtain a more complete view of cocaine-induced changes in gene expression and associated adaptations in numerous modes of chromatin regulation in the mouse nucleus accumbens, a key brain reward region. We demonstrate an unexpectedly large number of pre-mRNA splicing alterations in response to repeated cocaine treatment. In addition, we identify combinations of chromatin changes, or signatures, that correlate with cocaine-dependent regulation of gene expression, including those involving pre-mRNA alternative splicing. Through bioinformatic prediction and biological validation, we identify one particular splicing factor, A2BP1(Rbfox1/Fox-1), which is enriched at genes that display certain chromatin signatures and contributes to drug-induced behavioral abnormalities. Together, this delineation of the cocaine-induced epigenome in the nucleus accumbens reveals several novel modes of regulation by which cocaine alters the brain.

CONCLUSIONS

We establish combinatorial chromatin and transcriptional profiles in mouse nucleus accumbens after repeated cocaine treatment. These results serve as an important resource for the field and provide a template for the analysis of other systems to reveal new transcriptional and epigenetic mechanisms of neuronal regulation.

Exonic splicing signals impose constraints upon the evolution of enzymatic activity

Exon splicing enhancers (ESEs) overlap with amino acid coding sequences implying a dual evolutionary selective pressure. In this study, we map ESEs in the placental alkaline phosphatase gene (ALPP), absent in the corresponding exon of the ancestral tissue-non-specific alkaline phosphatase gene (ALPL). The ESEs are associated with amino acid differences between the transcripts in an area otherwise conserved. We switched out the ALPP ESEs sequences with the sequence from the related ALPL, introducing the associated amino acid changes. The resulting enzymes, produced by cDNA expression, showed different kinetic characteristics than ALPL and ALPP. In the organism, this enzyme will never be subjected to selection because gene splicing analysis shows exon skipping due to loss of the ESE. Our data prove that ESEs restrict the evolution of enzymatic activity. Thus, suboptimal proteins may exist in scenarios when coding nucleotide changes and consequent amino acid variation cannot be reconciled with the splicing function.

Global dissection of alternative splicing in paleopolyploid soybean

Alternative splicing (AS) is common in higher eukaryotes and plays an important role in gene posttranscriptional regulation. It has been suggested that AS varies dramatically among species, tissues, and duplicated gene families of different sizes. However, the genomic forces that govern AS variation remain poorly understood. Here, through genome-wide identification of AS events in the soybean (Glycine max) genome using high-throughput RNA sequencing of 28 samples from different developmental stages, we found that more than 63% of multiexonic genes underwent AS. More AS events occurred in the younger developmental stages than in the older developmental stages for the same type of tissue, and the four main AS types, exon skipping, intron retention, alternative donor sites, and alternative acceptor sites, exhibited different characteristics. Global computational analysis demonstrated that the variations of AS frequency and AS types were significantly correlated with the changes of gene features and gene transcriptional level. Further investigation suggested that the decrease of AS within the genome-wide duplicated genes were due to the diminution of intron length, exon number, and transcriptional level. Altogether, our study revealed that a large number of genes were alternatively spliced in the soybean genome and that variations in gene structure and transcriptional level may play important roles in regulating AS.