Context-dependent robustness to 5' splice site polymorphisms in human populations

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

Context matters: Regulation of splice donor usage

Elaborate research on splicing, starting in the late seventies, evolved from the discovery that 5' splice sites are recognized by their complementarity to U1 snRNA towards the realization that RNA duplex formation cannot be the sole basis for 5'ss selection. Rather, their recognition is highly influenced by a number of context factors including transcript architecture as well as splicing regulatory elements (SREs) in the splice site neighborhood. In particular, proximal binding of splicing regulatory proteins highly influences splicing outcome. The importance of SRE integrity especially becomes evident in the light of human pathogenic mutations where single nucleotide changes in SREs can severely affect the resulting transcripts. Bioinformatics tools nowadays greatly assist in the computational evaluation of 5'ss, their neighborhood and the impact of pathogenic mutations. Although predictions are already quite robust, computational evaluation of the splicing regulatory landscape still faces challenges to increase future reliability. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

The Arabidopsis sickle Mutant Exhibits Altered Circadian Clock Responses to Cool Temperatures and Temperature-Dependent Alternative Splicing

The circadian clock allows plants to anticipate and respond to daily changes in ambient temperature. Mechanisms establishing the timing of circadian rhythms in Arabidopsis thaliana through temperature entrainment remain unclear. Also incompletely understood is the temperature compensation mechanism that maintains consistent period length within a range of ambient temperatures. A genetic screen for Arabidopsis mutants affecting temperature regulation of the PSEUDO-RESPONSE REGULATOR7 promoter yielded a novel allele of the SICKLE (SIC) gene. This mutant, sic-3, and the existing sic-1 mutant both exhibit low-amplitude or arrhythmic expression of core circadian clock genes under cool ambient temperature cycles, but not under light-dark entrainment. sic mutants also lengthen free running period in a manner consistent with impaired temperature compensation. sic mutant alleles accumulate LATE ELONGATED HYPOCOTYL (LHY) and CIRCADIAN CLOCK ASSOCIATED1 (CCA1) splice variants, among other alternatively spliced transcripts, which is exacerbated by cool temperatures. The cca1-1 lhy-20 double mutant is epistatic to sic-3, indicating the LHY and CCA1 splice variants are needed for sic-3 circadian clock phenotypes. It is not expected that SIC is directly involved in the circadian clock mechanism; instead, SIC likely contributes to pre-mRNA metabolism, and the splice variants that accumulate in sic mutants likely affect the circadian clock response to cool ambient temperature.

Transcriptome sequencing reveals aberrant alternative splicing in Huntington's disease

Huntington's disease (HD) is an autosomal dominant neurodegenerative disorder caused by a CAG expansion in the gene-encoding Huntingtin (HTT). Transcriptome dysregulation is a major feature of HD pathogenesis, as revealed by a large body of work on gene expression profiling of tissues from human HD patients and mouse models. These studies were primarily focused on transcriptional changes affecting steady-state overall gene expression levels using microarray based approaches. A major missing component, however, has been the study of transcriptome changes at the post-transcriptional level, such as alternative splicing. Alternative splicing is a critical mechanism for expanding regulatory and functional diversity from a limited number of genes, and is particularly complex in the mammalian brain. Here we carried out a deep RNA-seq analysis of the BA4 (Brodmann area 4) motor cortex from seven human HD brains and seven controls to systematically discover aberrant alternative splicing events and characterize potential associated splicing factors in HD. We identified 593 differential alternative splicing events between HD and control brains. Using two expanded panels with a total of 108 BA4 tissues from patients and controls, we identified four splicing factors exhibiting significantly altered expression levels in HD patient brains. Moreover, follow-up molecular analyses of one splicing factor PTBP1 revealed its impact on disease-associated splicing patterns in HD. Collectively, our data provide genomic evidence for widespread splicing dysregulation in HD brains, and suggest the role of aberrant alternative splicing in the pathogenesis of HD.

Discover hidden splicing variations by mapping personal transcriptomes to personal genomes

RNA-seq has become a popular technology for studying genetic variation of pre-mRNA alternative splicing. Commonly used RNA-seq aligners rely on the consensus splice site dinucleotide motifs to map reads across splice junctions. Consequently, genomic variants that create novel splice site dinucleotides may produce splice junction RNA-seq reads that cannot be mapped to the reference genome. We developed and evaluated an approach to identify ‘hidden’ splicing variations in personal transcriptomes, by mapping personal RNA-seq data to personal genomes. Computational analysis and experimental validation indicate that this approach identifies personal specific splice junctions at a low false positive rate. Applying this approach to an RNA-seq data set of 75 individuals, we identified 506 personal specific splice junctions, among which 437 were novel splice junctions not documented in current human transcript annotations. 94 splice junctions had splice site SNPs associated with GWAS signals of human traits and diseases. These involve genes whose splicing variations have been implicated in diseases (such as OAS1), as well as novel associations between alternative splicing and diseases (such as ICA1). Collectively, our work demonstrates that the personal genome approach to RNA-seq read alignment enables the discovery of a large but previously unknown catalog of splicing variations in human populations.

Determinants of the Usage of Splice-Associated cis-Motifs Predict the Distribution of Human Pathogenic SNPs

Where in genes do pathogenic mutations tend to occur and does this provide clues as to the possible underlying mechanisms by which single nucleotide polymorphisms (SNPs) cause disease? As splice-disrupting mutations tend to occur predominantly at exon ends, known also to be hot spots of cis-exonic splice control elements, we examine the relationship between the relative density of such exonic cis-motifs and pathogenic SNPs. In particular, we focus on the intragene distribution of exonic splicing enhancers (ESE) and the covariance between them and disease-associated SNPs. In addition to showing that disease-causing genes tend to be genes with a high intron density, consistent with missplicing, five factors established as trends in ESE usage, are considered: relative position in exons, relative position in genes, flanking intron size, splice sites usage, and phase. We find that more than 76% of pathogenic SNPs are within 3–69 bp of exon ends where ESEs generally reside, this being 13% more than expected. Overall from enrichment of pathogenic SNPs at exon ends, we estimate that approximately 20–45% of SNPs affect splicing. Importantly, we find that within genes pathogenic SNPs tend to occur in splicing-relevant regions with low ESE density: they are found to occur preferentially in the terminal half of genes, in exons flanked by short introns and at the ends of phase (0,0) exons with 3′ non-“AGgt” splice site. We suggest the concept of the “fragile” exon, one home to pathogenic SNPs owing to its vulnerability to splice disruption owing to low ESE density.

Impact of G-quadruplex structures and intronic polymorphisms rs17878362 and rs1642785 on basal and ionizing radiation-induced expression of alternative p53 transcripts

G-quadruplex (G4) structures in intron 3 of the p53 pre-mRNA modulate intron 2 splicing, altering the balance between the fully spliced p53 transcript (FSp53, encoding full-length p53) and an incompletely spliced transcript retaining intron 2 (p53I2 encoding the N-terminally truncated Δ40p53 isoform). The nucleotides forming G4s overlap the polymorphism rs17878362 (A1 wild-type allele, A2 16-base pair insertion) which is in linkage disequilibrium with rs1642785 in intron 2 (c.74+38 G>C). Biophysical and biochemical analyses show rs17878362 A2 alleles form similar G4 structures as A1 alleles although their position is shifted with respect to the intron 2 splice acceptor site. In addition basal FSp53 and p53I2 levels showed allele specific differences in both p53-null cells transfected with reporter constructs or lymphoblastoid cell lines. The highest FSp53 and p53I2 levels were associated with combined rs1642785-GG/rs17878362-A1A1 alleles, whereas the presence of rs1642785-C with either rs17878362 allele was associated with lower p53 pre-mRNA, total TP53, FSp53 and p53I2 levels, due to the lower stability of transcripts containing rs1642785-C. Treatment of lymphoblastoid cell with the G4 binding ligands 360A or PhenDC3 or with ionizing radiation increased FSp53 levels only in cells with rs17878362 A1 alleles, suggesting that under this G4 configuration full splicing is favoured. These results demonstrate the complex effects of intronic TP53 polymorphisms on G4 formation and identify a new role for rs1642785 on mRNA splicing and stability, and thus on the differential expression of isoform-specific transcripts of the TP53 gene.

GLiMMPS: robust statistical model for regulatory variation of alternative splicing using RNA-seq data

To characterize the genetic variation of alternative splicing, we develop GLiMMPS, a robust statistical method for detecting splicing quantitative trait loci (sQTLs) from RNA-seq data. GLiMMPS takes into account the individual variation in sequencing coverage and the noise prevalent in RNA-seq data. Analyses of simulated and real RNA-seq datasets demonstrate that GLiMMPS outperforms competing statistical models. Quantitative RT-PCR tests of 26 randomly selected GLiMMPS sQTLs yielded a validation rate of 100%. As population-scale RNA-seq studies become increasingly affordable and popular, GLiMMPS provides a useful tool for elucidating the genetic variation of alternative splicing in humans and model organisms.

Screening of SLC26A4 gene in autoimmune thyroid diseases

The Pendred syndrome (PS) gene, SLC26A4, was involved in the genetic susceptibility of autoimmune thyroid disease (AITD) in Tunisian population. Recently, functional assays have shown a differential expression of SLC26A4 gene between Graves' disease (GD) and Hashimoto's thyroiditis (HT). Here, by the mean of DHPLC and HRM, we explored the 21 exons and their flanking intronic sequences of 128 patients affected with GD (n = 64) or HT (n = 64). The pathogenic effect of identified variations on splice was investigated using the web server HSF. Eighteen allelic variations were identified and ranged on missense, sens and splice variations. Nine identified variations (c.-66C>G, c.898A>C, c.1002-9A>C, c.1061T>C, c.1544 + 9G>T, c.1545-5T>G, c.1790T>C, c.1826T>G, c.2139T>G) were previously reported in hearing impairment studies. Forty-seven per cent (30/64) of GD patients and 37,5% (24/64) of HT patients present at least one variant in the explored sequences. Moreover, the analysis of the variant distribution between HT (9 (5'UTR), 12 exonic and 13 intronic) and GD (18 (5'UTR), 13 exonic and 5 intronic) patients showed a significant difference (χ² = 6.54, 2df, P = 0.03). Interestingly, missense changes (I300L, p.M283I, F354S and p.L597S) affected conserved residues of pendrin. On the other hand, the HSF analyses ascertain that some variants identified in HT disease are predicted to have a pathogenic effect on splice. In conclusion, our analysis of SLC26A4 sequence variations suggested a distinct genetics basis between HT and GD patients, which should be confirmed on a large cohort.

Global profiling of alternative splicing events and gene expression regulated by hnRNPH/F

In this study, we have investigated the global impact of heterogeneous nuclear Ribonuclear Protein (hnRNP) H/F-mediated regulation of splicing events and gene expression in oligodendrocytes. We have performed a genome-wide transcriptomic analysis at the gene and exon levels in Oli-neu cells treated with siRNA that targets hnRNPH/F compared to untreated cells using Affymetrix Exon Array. Gene expression levels and regulated exons were identified with the GenoSplice EASANA algorithm. Bioinformatics analyses were performed to determine the structural properties of G tracts that correlate with the function of hnRNPH/F as enhancers vs. repressors of exon inclusion. Different types of alternatively spliced events are regulated by hnRNPH/F. Intronic G tracts density, length and proximity to the 5′ splice site correlate with the hnRNPH/F enhancer function. Additionally, 6% of genes are differently expressed upon knock down of hnRNPH/F. Genes that regulate the transition of oligodendrocyte progenitor cells to oligodendrocytes are differentially expressed in hnRNPH/F depleted Oli-neu cells, resulting in a decrease of negative regulators and an increase of differentiation-inducing regulators. The changes were confirmed in developing oligodendrocytes in vivo. This is the first genome wide analysis of splicing events and gene expression regulated by hnRNPH/F in oligodendrocytes and the first report that hnRNPH/F regulate genes that are involved in the transition from oligodendrocyte progenitor cells to oligodendrocytes.

Comprehensive analysis of the genome transcriptome and proteome landscapes of three tumor cell lines

We here present a comparative genome, transcriptome and functional network analysis of three human cancer cell lines (A431, U251MG and U2OS), and investigate their relation to protein expression. Gene copy numbers significantly influenced corresponding transcript levels; their effect on protein levels was less pronounced. We focused on genes with altered mRNA and/or protein levels to identify those active in tumor maintenance. We provide comprehensive information for the three genomes and demonstrate the advantage of integrative analysis for identifying tumor-related genes amidst numerous background mutations by relating genomic variation to expression/protein abundance data and use gene networks to reveal implicated pathways.

Genetic variations and alternative splicing: the Glioma associated oncogene 1, GLI1

Alternative splicing is a post-transcriptional regulatory process that is attaining stronger recognition as a modulator of gene expression. Alternative splicing occurs when the primary RNA transcript is differentially processed into more than one mature RNAs. This is the result of a variable definition/inclusion of the exons, the sequences that are excised from the primary RNA to form the mature RNAs. Consequently, RNA expression can generate a collection of differentially spliced RNAs, which may distinctly influence subsequent biological events, such as protein synthesis or other biomolecular interactions. Still the mechanisms that control exon definition and exon inclusion are not fully clarified. This mini-review highlights advances in this field as well as the impact of single nucleotide polymorphisms in affecting splicing decisions. The Glioma-associated oncogene 1, GLI1, is taken as an example in addressing the role of nucleotide substitutions for splicing regulation.

Transcriptome landscape of the human placenta

BACKGROUND

The placenta is a key component in understanding the physiological processes involved in pregnancy. Characterizing genes critical for placental function can serve as a basis for identifying mechanisms underlying both normal and pathologic pregnancies. Detailing the placental tissue transcriptome could provide a valuable resource for genomic studies related to placental disease.

RESULTS

We have conducted a deep RNA sequencing (RNA-Seq) study on three tissue components (amnion, chorion, and decidua) of 5 human placentas from normal term pregnancies. We compared the placental RNA-Seq data to that of 16 other human tissues and observed a wide spectrum of transcriptome differences both between placenta and other human tissues and between distinct compartments of the placenta. Exon-level analysis of the RNA-Seq data revealed a large number of exons with differential splicing activities between placenta and other tissues, and 79% (27 out of 34) of the events selected for RT-PCR test were validated. The master splicing regulator ESRP1 is expressed at a proportionately higher level in amnion compared to all other analyzed human tissues, and there is a significant enrichment of ESRP1-regulated exons with tissue-specific splicing activities in amnion. This suggests an important role of alternative splicing in regulating gene function and activity in specific placental compartments. Importantly, genes with differential expression or splicing in the placenta are significantly enriched for genes implicated in placental abnormalities and preterm birth. In addition, we identified 604-1007 novel transcripts and 494-585 novel exons expressed in each of the three placental compartments.

CONCLUSIONS

Our data demonstrate unique aspects of gene expression and splicing in placental tissues that provide a basis for disease investigation related to disruption of these mechanisms. These data are publicly available providing the community with a rich resource for placental physiology and disease-related studies.

MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data

Ultra-deep RNA sequencing has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We develop MATS (multivariate analysis of transcript splicing), a Bayesian statistical framework for flexible hypothesis testing of differential alternative splicing patterns on RNA-Seq data. MATS uses a multivariate uniform prior to model the between-sample correlation in exon splicing patterns, and a Markov chain Monte Carlo (MCMC) method coupled with a simulation-based adaptive sampling procedure to calculate the P-value and false discovery rate (FDR) of differential alternative splicing. Importantly, the MATS approach is applicable to almost any type of null hypotheses of interest, providing the flexibility to identify differential alternative splicing events that match a given user-defined pattern. We evaluated the performance of MATS using simulated and real RNA-Seq data sets. In the RNA-Seq analysis of alternative splicing events regulated by the epithelial-specific splicing factor ESRP1, we obtained a high RT–PCR validation rate of 86% for differential exon skipping events with a MATS FDR of <10%. Additionally, over the full list of RT–PCR tested exons, the MATS FDR estimates matched well with the experimental validation rate. Our results demonstrate that MATS is an effective and flexible approach for detecting differential alternative splicing from RNA-Seq data.

Genetic variation of pre-mRNA alternative splicing in human populations

The precise splicing outcome of a transcribed gene is controlled by complex interactions between cis regulatory splicing signals and trans-acting regulators. In higher eukaryotes, alternative splicing is a prevalent mechanism for generating transcriptome and proteome diversity. Alternative splicing can modulate gene function, affect organismal phenotype and cause disease. Common genetic variation that affects splicing regulation can lead to differences in alternative splicing between human individuals and consequently impact expression level or protein function. In several well-documented examples, such natural variation of alternative splicing has indeed been shown to influence disease susceptibility and drug response. With new microarray and sequencing-based genomic technologies that can analyze eukaryotic transcriptomes at the exon or nucleotide level, it has become possible to globally compare the alternative splicing profiles across human individuals in any tissue or cell type of interest. Recent large-scale transcriptome studies using high-density splicing-sensitive microarray and deep RNA sequencing (RNA-Seq) have revealed widespread genetic variation of alternative splicing in humans. In the future, an extensive catalog of alternative splicing variation in human populations will help elucidate the molecular underpinnings of complex traits and human diseases, and shed light on the mechanisms of splicing regulation in human cells.