Unconstrained mining of transcript data reveals increased alternative splicing complexity in the human transcriptome

Alternative splicing of clock transcript mediates the response of circadian clocks to temperature changes

Circadian clocks respond to temperature changes over the calendar year, allowing organisms to adjust their daily biological rhythms to optimize health and fitness. In Drosophila, seasonal adaptations and temperature compensation are regulated by temperature-sensitive alternative splicing (AS) of period (per) and timeless (tim) genes that encode key transcriptional repressors of clock gene expression. Although clock (clk) gene encodes the critical activator of clock gene expression, AS of its transcripts and its potential role in temperature regulation of clock function have not been explored. We therefore sought to investigate whether clk exhibits AS in response to temperature and the functional changes of the differentially spliced transcripts. We observed that clk transcripts indeed undergo temperature-sensitive AS. Specifically, cold temperature leads to the production of an alternative clk transcript, hereinafter termed clk-cold, which encodes a CLK isoform with an in-frame deletion of four amino acids proximal to the DNA binding domain. Notably, serine 13 (S13), which we found to be a CK1α-dependent phosphorylation site, is among the four amino acids deleted in CLK-cold protein. Using a combination of transgenic fly, tissue culture, and in vitro experiments, we demonstrated that upon phosphorylation at CLK(S13), CLK-DNA interaction is reduced, thus decreasing CLK occupancy at clock gene promoters. This is in agreement with our findings that CLK occupancy at clock genes and transcriptional output are elevated at cold temperature, which can be explained by the higher amounts of CLK-cold isoforms that lack S13 residue. This study provides new insights into the complex collaboration between AS and phospho-regulation in shaping temperature responses of the circadian clock.

Spliceostatin C, a component of a microbial bioherbicide, is a potent phytotoxin that inhibits the spliceosome

Spliceostatin C (SPC) is a component of a bioherbicide isolated from the soil bacterium Burkholderia rinojensis. The chemical structure of SPC closely resembles spliceostatin A (SPA) which was characterized as an anticancer agent and splicing inhibitor. SPC inhibited the growth of Arabidopsis thaliana seedlings with an IC50 value of 2.2 µM. The seedlings exposed to SPC displayed a significant response with decreased root length and number and inhibition of gravitropism. Reverse transcriptase semi-quantitative PCR (RT-sqPCR) analyses of 19 selected genes demonstrated the active impact of SPC on the quality and quantity of transcripts that underwent intron rearrangements as well as up or down expression upon exposure to SPC. Qualitative and quantitative proteomic profiles identified 66 proteins that were significantly affected by SPC treatment. Further proteomics data analysis revealed that spliceostatin C induces hormone-related responses in Arabidopsis seedlings. In silico binding studies showed that SPC binds to a pocket between the SF3B3 and PF5A of the spliceosome.

Intron retention by a novel intronic mutation in DKC1 gene caused recurrent still birth and early death in a Chinese family

Background

DKC1, the dyskerin encoding gene, functions in telomerase activity and telomere maintenance. DKC1 mutations cause a multisystem disease, dyskeratosis congenita (DC), which is associated with immunodeficiency and bone marrow failure.

Methods

In this research, we reported a novel intronic mutation of DKC1 causing dyskerin functional loss in a Chinese family. Whole exome sequence (WES) of the proband and validation by sanger sequencing help us identify a pathogenic DKC1 mutation. Minigene splicing assays were performed to evaluate functional change of DKC1.

Results

A pathogenic DKC1 intronic mutation(c.84 + 7A > G) was identified in the proband, which was inherited from heterozygous mother and not reported before. We detected the novel transcript with a 7 bp intron retention through minigene splicing assay. The newly spliced transcript is so short that would be degraded by nonsense‐mediated mRNA decay in vitro and we infer that the novel DKC1 mutation would influences normal physiological function of dyskerin.

Conclusions

Our study identified a novel intronic mutation, which expands the spectrum of pathogenic DKC1 gene mutations and can be used in molecular diagnosis. The mutant allele was transmitted to the next generation with high frequency in the family and causes still birth or early death.

RNA-Seq is not required to determine stable reference genes for qPCR normalization

Assessment of differential gene expression by qPCR is heavily influenced by the choice of reference genes. Although numerous statistical approaches have been proposed to determine the best reference genes, they can give rise to conflicting results depending on experimental conditions. Hence, recent studies propose the use of RNA-Seq to identify stable genes followed by the application of different statistical approaches to determine the best set of reference genes for qPCR data normalization. In this study, however, we demonstrate that the statistical approach to determine the best reference genes from commonly used conventional candidates is more important than the preselection of ‘stable’ candidates from RNA-Seq data. Using a qPCR data normalization workflow that we have previously established; we show that qPCR data normalization using conventional reference genes render the same results as stable reference genes selected from RNA-Seq data. We validated these observations in two distinct cross-sectional experimental conditions involving human iPSC derived microglial cells and mouse sciatic nerves. These results taken together show that given a robust statistical approach for reference gene selection, stable genes selected from RNA-Seq data do not offer any significant advantage over commonly used reference genes for normalizing qPCR assays.

RTqPCR is a powerful technique that is widely used to quantify gene expression in research and diagnostics of different diseases. The technique involves making multiple copies (amplification) of a specific target DNA. The amplified target DNA binds to a molecule that emits fluorescence upon binding. The extent of fluorescence correlates to the amount of DNA present. To precisely quantify this fluorescence (and thus the quantities of target DNA), internal control genes also called as reference genes need to be determined. Such genes, in principle, do not have varied expression across samples and would exhibit the same fluorescence in all samples. They can thus be used to normalize the expression of the Target DNA. Unfortunately, choosing the right reference gene is very tricky and poor choice of reference genes results in unreliable data both in research and in diagnostics. In this study, we validate a statistical approach to find stably expressed reference genes for any experimental setting using a given set of candidates. We compare our approach to RNA sequencing which quantifies the expression of thousands of genes at the same time. We highlight the advantages of our approach which is cost effective and saves a lot of time when compared to sequencing.

Alternative splicing level related to intron size and organism complexity

Background

Alternative splicing is the process of selecting different combinations of splice sites to produce variably spliced mRNAs. However, the relationships between alternative splicing prevalence and level (ASP/L) and variations of intron size and organism complexity (OC) remain vague. Here, we developed a robust protocol to analyze the relationships between ASP/L and variations of intron size and OC. Approximately 8 Tb raw RNA-Seq data from 37 eumetazoan species were divided into three sets of species based on variations in intron size and OC.

Results

We found a strong positive correlation between ASP/L and OC, but no correlation between ASP/L and intron size across species. Surprisingly, ASP/L displayed a positive correlation with mean intron size of genes within individual genomes. Moreover, our results revealed that four ASP/L-related pathways contributed to the differences in ASP/L that were associated with OC. In particular, the spliceosome pathway displayed distinct genomic features, such as the highest gene expression level, conservation level, and fraction of disordered regions. Interestingly, lower or no obvious correlations were observed among these genomic features.

Conclusions

The positive correlation between ASP/L and OC ubiquitously exists in eukaryotes, and this correlation is not affected by the mean intron size of these species. ASP/L-related splicing factors may play an important role in the evolution of OC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08172-2.

AI applications in functional genomics

We review the current applications of artificial intelligence (AI) in functional genomics. The recent explosion of AI follows the remarkable achievements made possible by "deep learning", along with a burst of "big data" that can meet its hunger. Biology is about to overthrow astronomy as the paradigmatic representative of big data producer. This has been made possible by huge advancements in the field of high throughput technologies, applied to determine how the individual components of a biological system work together to accomplish different processes. The disciplines contributing to this bulk of data are collectively known as functional genomics. They consist in studies of: i) the information contained in the DNA (genomics); ii) the modifications that DNA can reversibly undergo (epigenomics); iii) the RNA transcripts originated by a genome (transcriptomics); iv) the ensemble of chemical modifications decorating different types of RNA transcripts (epitranscriptomics); v) the products of protein-coding transcripts (proteomics); and vi) the small molecules produced from cell metabolism (metabolomics) present in an organism or system at a given time, in physiological or pathological conditions. After reviewing main applications of AI in functional genomics, we discuss important accompanying issues, including ethical, legal and economic issues and the importance of explainability.

FLAME: long-read bioinformatics tool for comprehensive spliceome characterization

Comprehensive characterization of differentially spliced RNA transcripts with nanopore sequencing is limited by bioinformatics tools that are reliant on existing annotations. We have developed FLAME, a bioinformatics pipeline for alternative splicing analysis of gene-specific or transcriptome-wide long-read sequencing data. FLAME is a Python-based tool aimed at providing comprehensible quantification of full-length splice variants, reliable de novo recognition of splice sites and exons, and representation of consecutive exon connectivity in the form of a weighted adjacency matrix. Notably, this workflow circumvents issues related to inadequate reference annotations and allows for incorporation of short-read sequencing data to improve the confidence of nanopore sequencing reads. In this study, the Epstein-Barr virus long noncoding RNA RPMS1 was used to demonstrate the utility of the pipeline. RPMS1 is ubiquitously expressed in Epstein-Barr virus associated cancer and known to undergo ample differential splicing. To fully resolve the RPMS1 spliceome, we combined gene-specific nanopore sequencing reads from a primary gastric adenocarcinoma and a nasopharyngeal carcinoma cell line with matched publicly available short-read sequencing data sets. All previously reported splice variants, including putative ORFs, were detected using FLAME. In addition, 32 novel exons, including two intron retentions and a cassette exon, were discovered within the RPMS1 gene.

Retracing Schwann Cell Developmental Transitions in Embryonic Dissociated DRG/Schwann Cell Cocultures in Mice

Embryonic Dissociated Dorsal Root Ganglia (DRG) cultures are often used to investigate the role of novel molecular pathways or drugs in Schwann cell development and myelination. These cultures largely recapitulate the order of cellular and molecular events that occur in Schwann cells of embryonic nerves. However, the timing of Schwann cell developmental transitions, notably the transition from Schwann Cell Precursors (SCP) to immature Schwann cells (iSC) and then to myelinating Schwann cells, has not been estimated so far in this culture system. In this study, we determined the expression profiles of Schwann cell developmental genes during the first week of culture and then compared our data to the expression profiles of these genes in developing spinal nerves. This helped in identifying that SCP transition into iSC between the 5th and 7th day in vitro. Furthermore, we also investigated the transition of immature cells into pro-myelinating and myelinating Schwann cells upon the induction of myelination in vitro. Our results suggest that Schwann cell differentiation beyond the immature stage can be observed as early as 4 days post the induction of myelination in cocultures. Finally, we compared the myelinating potential of coculture-derived Schwann cell monocultures to cultures established from neonatal sciatic nerves and found that both these culture systems exhibit similar myelinating phenotypes. In effect, our results allow for a better understanding and interpretation of coculture experiments especially in studies that aim to elucidate the role of a novel actor in Schwann cell development and myelination.

Isoform-resolved correlation analysis between mRNA abundance regulation and protein level degradation

Profiling of biological relationships between different molecular layers dissects regulatory mechanisms that ultimately determine cellular function. To thoroughly assess the role of protein post‐translational turnover, we devised a strategy combining pulse stable isotope‐labeled amino acids in cells (pSILAC), data‐independent acquisition mass spectrometry (DIA‐MS), and a novel data analysis framework that resolves protein degradation rate on the level of mRNA alternative splicing isoforms and isoform groups. We demonstrated our approach by the genome‐wide correlation analysis between mRNA amounts and protein degradation across different strains of HeLa cells that harbor a high grade of gene dosage variation. The dataset revealed that specific biological processes, cellular organelles, spatial compartments of organelles, and individual protein isoforms of the same genes could have distinctive degradation rate. The protein degradation diversity thus dissects the corresponding buffering or concerting protein turnover control across cancer cell lines. The data further indicate that specific mRNA splicing events such as intron retention significantly impact the protein abundance levels. Our findings support the tight association between transcriptome variability and proteostasis and provide a methodological foundation for studying functional protein degradation.

Systematic Analyses of the Differentially Expressed Alternative Splicing Events in Gastric Cancer and Its Clinical Significance

Accumulation of evidence has indicated a close relationship between alternative splicing (AS) and gastric cancer (GC), whereas systematic analyses of the differentially expressed AS events (DEAS) between GC and normal tissues are lacking. RNA-Seq data and the corresponding clinical information were downloaded from TCGA GC cohort. The percent spliced-in (PSI) value calculated in the GC tissues and normal tissues was employed to quantify the DEAS. Further, survival-associated DEAS and DEAS signatures were identified by univariate and multivariate cox regression analyses. To evaluate the association between DEAS and patients' clinical features, Kaplan-Meier analysis, receiver operator characteristic (ROC) curve, Cox proportional regression and nomograms incorporating the DEAS signatures were performed. DEAS and their splicing networks were finally analyzed by bioinformatics methods. In addition, we use the method of random grouping to divide the samples into the training group and the test group. The final results of the two groups are consistent. After strict filtering, a total of 44,935 AS events were identified, among which 11,141 DEAS were preliminarily screened from 5032 genes. A total of 454 DEAS was associated with OS, and 872 DEAS were associated with DFS. The final prognostic signatures were constructed from the survival-associated DEAS with an area under the receiver operating characteristic (ROC) curve (AUC) greater than 0.6. Only ES in ABI1 was simultaneously associated with OS and DFS. Finally, we identified the splicing correlation network between the prognostic splicing factors (SF) and DEAS in GC. Our study provided a systematic portrait of survival-associated DEAS in GC and uncovered splicing networks that are valuable in deciphering the underlying mechanisms of AS in GC.

Tissue and cancer-specific expression of DIEXF is epigenetically mediated by an Alu repeat

Alu repeats constitute a major fraction of human genome and for a small subset of them a role in gene regulation has been described. The number of studies focused on the functional characterization of particular Alu elements is very limited. Most Alu elements are DNA methylated and then assumed to lie in repressed chromatin domains. We hypothesize that Alu elements with low or variable DNA methylation are candidates for a functional role. In a genome-wide study in normal and cancer tissues, we pinpointed an Alu repeat (AluSq2) with differential methylation located upstream of the promoter region of the DIEXF gene. DIEXF encodes a highly conserved factor essential for the development of zebrafish digestive tract. To characterize the contribution of the Alu element to the regulation of DIEXF we analysed the epigenetic landscapes of the gene promoter and flanking regions in different cell types and cancers. Alternate epigenetic profiles (DNA methylation and histone modifications) of the AluSq2 element were associated with DIEXF transcript diversity as well as protein levels, while the epigenetic profile of the CpG island associated with the DIEXF promoter remained unchanged. These results suggest that AluSq2 might directly contribute to the regulation of DIEXF transcription and protein expression. Moreover, AluSq2 was DNA hypomethylated in different cancer types, pointing out its putative contribution to DIEXF alteration in cancer and its potential as tumoural biomarker.

Splicing the Clock to Maintain and Entrain Circadian Rhythms

Circadian clocks drive daily rhythms of physiology and behavior in multiple organisms and synchronize these rhythms to environmental cycles of light and temperature. The basic mechanism of the clock consists of a transcription-translation feedback loop, in which key clock proteins negatively regulate their own transcription. Although much of the focus with respect to clock mechanisms has been on the regulation of transcription and on the stability and activity of clock proteins, it is clear that other regulatory processes also have to be involved to explain aspects of clock function. Here, we review the role of alternative splicing in circadian clocks. Starting with a discussion of the Drosophila clock and then extending to other major circadian model systems, we describe how the control of alternative splicing enables organisms to maintain their circadian clocks as well as to respond to environmental inputs, in particular to temperature changes.

A novel L1CAM isoform with angiogenic activity generated by NOVA2-mediated alternative splicing

The biological players involved in angiogenesis are only partially defined. Here, we report that endothelial cells (ECs) express a novel isoform of the cell-surface adhesion molecule L1CAM, termed L1-ΔTM. The splicing factor NOVA2, which binds directly to L1CAM pre-mRNA, is necessary and sufficient for the skipping of L1CAM transmembrane domain in ECs, leading to the release of soluble L1-ΔTM. The latter exerts high angiogenic function through both autocrine and paracrine activities. Mechanistically, L1-ΔTM-induced angiogenesis requires fibroblast growth factor receptor-1 signaling, implying a crosstalk between the two molecules. NOVA2 and L1-ΔTM are overexpressed in the vasculature of ovarian cancer, where L1-ΔTM levels correlate with tumor vascularization, supporting the involvement of NOVA2-mediated L1-ΔTM production in tumor angiogenesis. Finally, high NOVA2 expression is associated with poor outcome in ovarian cancer patients. Our results point to L1-ΔTM as a novel, EC-derived angiogenic factor which may represent a target for innovative antiangiogenic therapies.

Growing tumors stimulate the formation of new blood vessels to supply the oxygen and nutrients the cancerous cells need to stay alive. Stopping tumors from forming the blood vessels could therefore help us to treat cancer. To do so, we need to understand how different proteins control when and how blood vessels develop.

Cells make proteins by first ‘transcribing’ genes to form RNA molecules. In many cases, the RNA then goes through a process called alternative splicing. Proteins known as splicing factors cut out different segments of the RNA molecule and stick together the remaining segments to form templates for protein production. This enables a single gene to produce many different variants of a protein.

Angiolini, Belloni, Giordano et al. have now studied mouse and human versions of the cells that line the blood vessels grown by tumors. This revealed that a splicing factor called NOVA2 targets a protein called L1CAM, which is normally responsible for gluing adjacent cells together. Angiolini et al. found that NOVA2 splices L1CAM into a form not seen before. Instead of remaining anchored to cell surfaces, the newly identified form of L1CAM is released into the blood circulation, where it stimulates new blood vessels to grow.

Samples taken from the blood vessels of human ovarian tumors showed high levels of both NOVA2 and the modified form of L1CAM, while blood vessels in healthy tissue contain no, or very low levels of both proteins. Therefore, if the new form of L1CAM can be detected in the blood, it could be used to help cancer diagnosis, and to indicate which patients would benefit from treatments that restrict the growth of blood vessels in tumors. Further work is now needed to explore these possibilities.

The cyclin-dependent kinase G group defines a thermo-sensitive alternative splicing circuit modulating the expression of Arabidopsis ATU2AF65A

The ability to adapt growth and development to temperature variations is crucial to generate plant varieties resilient to predicted temperature changes. However, the mechanisms underlying plant response to progressive increases in temperature have just started to be elucidated. Here, we report that the cyclin-dependent kinase G1 (CDKG1) is a central element in a thermo-sensitive mRNA splicing cascade that transduces changes in ambient temperature into differential expression of the fundamental spliceosome component, ATU2AF65A. CDKG1 is alternatively spliced in a temperature-dependent manner. We found that this process is partly dependent on both the cyclin-dependent kinase G2 (CDKG2) and the interacting co-factor CYCLIN L1 (CYCL1), resulting in two distinct messenger RNAs. The relative abundance of both CDKG1 transcripts correlates with ambient temperature and possibly with different expression levels of the associated protein isoforms. Both CDKG1 alternative transcripts are necessary to fully complement the expression of ATU2AF65A across the temperature range. Our data support a previously unidentified temperature-dependent mechanism based on the alternative splicing (AS) of CDKG1 and regulated by CDKG2 and CYCL1. We propose that changes in ambient temperature affect the relative abundance of CDKG1 transcripts, and this in turn translates into differential CDKG1 protein expression coordinating the AS of ATU2AF65A.

Alternative splicing affects the targeting sequence of peroxisome proteins in Arabidopsis

A systematic analysis of the Arabidopsis genome in combination with localization experiments indicates that alternative splicing affects the peroxisomal targeting sequence of at least 71 genes in Arabidopsis. Peroxisomes are ubiquitous eukaryotic cellular organelles that play a key role in diverse metabolic functions. All peroxisome proteins are encoded by nuclear genes and target to peroxisomes mainly through two types of targeting signals: peroxisomal targeting signal type 1 (PTS1) and PTS2. Alternative splicing (AS) is a process occurring in all eukaryotes by which a single pre-mRNA can generate multiple mRNA variants, often encoding proteins with functional differences. However, the effects of AS on the PTS1 or PTS2 and the targeting of the protein were rarely studied, especially in plants. Here, we systematically analyzed the genome of Arabidopsis, and found that the C-terminal targeting sequence PTS1 of 66 genes and the N-terminal targeting sequence PTS2 of 5 genes are affected by AS. Experimental determination of the targeting of selected protein isoforms further demonstrated that AS at both the 5' and 3' region of a gene can affect the inclusion of PTS2 and PTS1, respectively. This work underscores the importance of AS on the global regulation of peroxisome protein targeting.

Alternative splicing shapes transcriptome but not proteome diversity in Physcomitrella patens

Alternative splicing (AS) can significantly impact the transcriptome and proteome of a eukaryotic cell. Here, using transcriptome and proteome profiling data, we analyzed AS in two life forms of the model moss Physcomitrella patens, namely protonemata and gametophores, as well as in protoplasts. We identified 12 043 genes subject to alternative splicing and analyzed the extent to which AS contributes to proteome diversity. We could distinguish a few examples that unambiguously indicated the presence of two or more splice isoforms from the same locus at the proteomic level. Our results indicate that alternative isoforms have a small effect on proteome diversity. We also revealed that mRNAs and pre-mRNAs have thousands of complementary binding sites for long non-coding RNAs (lncRNAs) that may lead to potential interactions in transcriptome. This finding points to an additional level of gene expression and AS regulation by non-coding transcripts in Physcomitrella patens. Among the differentially expressed and spliced genes we found serine/arginine-rich (SR) genes, which are known to regulate AS in cells. We found that treatment with abscisic (ABA) and methyl jasmonic acids (MeJA) led to an isoform-specific response and suggested that ABA in gametophores and MeJA in protoplasts regulate AS and the transcription of SR genes.

Alternative Splicing May Not Be the Key to Proteome Complexity

Alternative splicing is commonly believed to be a major source of cellular protein diversity. However, although many thousands of alternatively spliced transcripts are routinely detected in RNA-seq studies, reliable large-scale mass spectrometry-based proteomics analyses identify only a small fraction of annotated alternative isoforms. The clearest finding from proteomics experiments is that most human genes have a single main protein isoform, while those alternative isoforms that are identified tend to be the most biologically plausible: those with the most cross-species conservation and those that do not compromise functional domains. Indeed, most alternative exons do not seem to be under selective pressure, suggesting that a large majority of predicted alternative transcripts may not even be translated into proteins.

Low Dose Iron Treatments Induce a DNA Damage Response in Human Endothelial Cells within Minutes

BACKGROUND

Spontaneous reports from patients able to report vascular sequelae in real time, and recognition that serum non transferrin bound iron may reach or exceed 10μmol/L in the blood stream after iron tablets or infusions, led us to hypothesize that conventional iron treatments may provoke acute vascular injury. This prompted us to examine whether a phenotype could be observed in normal human endothelial cells treated with low dose iron.

METHODOLOGY

Confluent primary human endothelial cells (EC) were treated with filter-sterilized iron (II) citrate or fresh media for RNA sequencing and validation studies. RNA transcript profiles were evaluated using directional RNA sequencing with no pre-specification of target sequences. Alignments were counted for exons and junctions of the gene strand only, blinded to treatment types.

PRINCIPAL FINDINGS

Rapid changes in RNA transcript profiles were observed in endothelial cells treated with 10μmol/L iron (II) citrate, compared to media-treated cells. Clustering for Gene Ontology (GO) performed on all differentially expressed genes revealed significant differences in biological process terms between iron and media-treated EC, whereas 10 sets of an equivalent number of randomly selected genes from the respective EC gene datasets showed no significant differences in any GO terms. After 1 hour, differentially expressed genes clustered to vesicle mediated transport, protein catabolism, and cell cycle (Benjamini p = 0.0016, 0.0024 and 0.0032 respectively), and by 6 hours, to cellular response to DNA damage stimulus most significantly through DNA repair genes FANCG, BLM, and H2AFX. Comet assays demonstrated that 10μM iron treatment elicited DNA damage within 1 hour. This was accompanied by a brisk DNA damage response pulse, as ascertained by the development of DNA damage response (DDR) foci, and p53 stabilization.

SIGNIFICANCE

These data suggest that low dose iron treatments are sufficient to modify the vascular endothelium, and induce a DNA damage response.

Evolutionary Character of Alternative Splicing in Plants

Alternative splicing (AS) is one of the most important ways to enhance the functional diversity of genes. Huge amounts of data have been produced by microarray, expressed sequence tag, and RNA-seq, and plenty of methods have been developed specifically for this task. The most frequently asked questions in previous research were as follows. What is the content rate of AS genes among the whole gene set? How many AS types are presented in the genome, and which type is dominant? How about the conservation ability of AS among different species? Which kinds of isoforms from some genes have the environmental response to help individual adaptation? Based on this background, we collected analysis results from 17 species to try to map out the landscape of AS studies in plants. We have noted the shortages of previous results, and we appeal to all scientists working in the AS field to make a standard protocol so that analyses between different projects are comparable.

JuncDB: an exon-exon junction database

Intron positions upon the mRNA transcript are sometimes remarkably conserved even across distantly related eukaryotic species. This has made the comparison of intron-exon architectures across orthologous transcripts a very useful tool for studying various evolutionary processes. Moreover, the wide range of functions associated with introns may confer biological meaning to evolutionary changes in gene architectures. Yet, there is currently no database that offers such comparative information. Here, we present JuncDB (http://juncdb.carmelab.huji.ac.il/), an exon-exon junction database dedicated to the comparison of architectures between orthologous transcripts. It covers nearly 40,000 sets of orthologous transcripts spanning 88 eukaryotic species. JuncDB offers a user-friendly interface, access to detailed information, instructive graphical displays of the comparative data and easy ways to download data to a local computer. In addition, JuncDB allows the analysis to be carried out either on specific genes, or at a genome-wide level for any selected group of species.

RNA Dynamics in the Control of Circadian Rhythm

The circadian oscillator is based on transcription-translation feedback loops that generate 24 h oscillations in gene expression. Although circadian regulation of mRNA expression at the transcriptional level is one of the most important steps for the generation of circadian rhythms within the cell, multiple lines of evidence point to a disconnect between transcript oscillation and protein oscillation. This can be explained by regulatory RNA-binding proteins acting on the nascent transcripts to modulate their processing, export, translation and degradation rates. In this chapter we will review what is known about the different steps involved in circadian gene expression from transcription initiation to mRNA stability and translation efficiency. The role of ribonucleoprotein particles in the generation of rhythmic gene expression is only starting to be elucidated, but it is likely that they cooperate with the basal transcriptional machinery to help to maintain the precision of the clock under diverse cellular and environmental conditions.

Comparative analyses between retained introns and constitutively spliced introns in Arabidopsis thaliana using random forest and support vector machine

One of the important modes of pre-mRNA post-transcriptional modification is alternative splicing. Alternative splicing allows creation of many distinct mature mRNA transcripts from a single gene by utilizing different splice sites. In plants like Arabidopsis thaliana, the most common type of alternative splicing is intron retention. Many studies in the past focus on positional distribution of retained introns (RIs) among different genic regions and their expression regulations, while little systematic classification of RIs from constitutively spliced introns (CSIs) has been conducted using machine learning approaches. We used random forest and support vector machine (SVM) with radial basis kernel function (RBF) to differentiate these two types of introns in Arabidopsis. By comparing coordinates of introns of all annotated mRNAs from TAIR10, we obtained our high-quality experimental data. To distinguish RIs from CSIs, We investigated the unique characteristics of RIs in comparison with CSIs and finally extracted 37 quantitative features: local and global nucleotide sequence features of introns, frequent motifs, the signal strength of splice sites, and the similarity between sequences of introns and their flanking regions. We demonstrated that our proposed feature extraction approach was more accurate in effectively classifying RIs from CSIs in comparison with other four approaches. The optimal penalty parameter C and the RBF kernel parameter in SVM were set based on particle swarm optimization algorithm (PSOSVM). Our classification performance showed F-Measure of 80.8% (random forest) and 77.4% (PSOSVM). Not only the basic sequence features and positional distribution characteristics of RIs were obtained, but also putative regulatory motifs in intron splicing were predicted based on our feature extraction approach. Clearly, our study will facilitate a better understanding of underlying mechanisms involved in intron retention.

Evolution of gene structural complexity: an alternative-splicing-based model accounts for intron-containing retrogenes

The structure of eukaryotic genes evolves extensively by intron loss or gain. Previous studies have revealed two models for gene structure evolution through the loss of introns: RNA-based gene conversion, dubbed the Fink model and retroposition model. However, retrogenes that experienced both intron loss and intron-retaining events have been ignored; evolutionary processes responsible for the variation in complex exon-intron structure were unknown. We detected hundreds of retroduplication-derived genes in human (Homo sapiens), fly (Drosophila melanogaster), rice (Oryza sativa), and Arabidopsis (Arabidopsis thaliana) and categorized them either as duplicated genes that have all introns lost or as duplicated genes that have at least lost one and retained one intron compared with the parental copy (intron-retaining [IR] type). Our new model attributes intron retention alternative splicing to the generation of these IR-type gene pairs. We presented 25 parental genes that have an intron retention isoform and have retained introns in the same locations in the IR-type duplicate genes, which directly support our hypothesis. Our alternative-splicing-based model in conjunction with the retroposition and Fink models can explain the IR-type gene observed. We discovered a greater percentage of IR-type genes in plants than in animals, which may be due to the abundance of intron retention cases in plants. Given the prevalence of intron retention in plants, this new model gives a support that plant genomes have very complex gene structures.

Correcting for differential transcript coverage reveals a strong relationship between alternative splicing and organism complexity

What at the genomic level underlies organism complexity? Although several genomic features have been associated with organism complexity, in the case of alternative splicing, which has long been proposed to explain the variation in complexity, no such link has been established. Here, we analyzed over 39 million expressed sequence tags available for 47 eukaryotic species with fully sequenced genomes to obtain a comparable index of alternative splicing estimates, which corrects for the distorting effect of a variable number of transcripts per species—an important obstacle for comparative studies of alternative splicing. We find that alternative splicing has steadily increased over the last 1,400 My of eukaryotic evolution and is strongly associated with organism complexity, assayed as the number of cell types. Importantly, this association is not explained as a by-product of covariance between alternative splicing with other variables previously linked to complexity including gene content, protein length, proteome disorder, and protein interactivity. In addition, we found no evidence to suggest that the relationship of alternative splicing to cell type number is explained by drift due to reduced N_e in more complex species. Taken together, our results firmly establish alternative splicing as a significant predictor of organism complexity and are, in principle, consistent with an important role of transcript diversification through alternative splicing as a means of determining a genome’s functional information capacity.

Intron retention: a human DKC1 gene common splicing event

Identification of alternatively spliced transcripts produced by a gene is a crucial step in deciphering the bulk of its biological roles and the overall processes that regulate its activity. By using a combination of bioinformatic and molecular approaches we identified, cloned, and characterized 3 novel alternative splice isoforms derived from human dyskeratosis congenita 1 (hDKC1), an essential human gene causative of the X-linked dyskeratosis congenita disease and involved in multiple functions related to cell growth, proliferation, and telomere maintenance. Expression of the new isoforms, all characterized by intron retention, was confirmed by RT-PCR in a panel of diverse cell lines and normal human tissues, and despite the presence of premature termination codons, was not down-regulated by the mechanism of nonsense-mediated decay. Accumulation of these transcripts fluctuated distinctly in the diverse tissues and during in vitro differentiation of Caco2 cells, suggesting that their ratio may contribute to the gene functional diversity across different cell types. Intriguingly, the structure of one isoform leads to exonize an intronically encoded small nucleolar RNA (snoRNA), highlighting an additional layer of complexity that can contribute to overall gene regulation.

ExoLocator--an online view into genetic makeup of vertebrate proteins

ExoLocator (http://exolocator.eopsf.org) collects in a single place information needed for comparative analysis of protein-coding exons from vertebrate species. The main source of data—the genomic sequences, and the existing exon and homology annotation—is the ENSEMBL database of completed vertebrate genomes. To these, ExoLocator adds the search for ostensibly missing exons in orthologous protein pairs across species, using an extensive computational pipeline to narrow down the search region for the candidate exons and find a suitable template in the other species, as well as state-of-the-art implementations of pairwise alignment algorithms. The resulting complements of exons are organized in a way currently unique to ExoLocator: multiple sequence alignments, both on the nucleotide and on the peptide levels, clearly indicating the exon boundaries. The alignments can be inspected in the web-embedded viewer, downloaded or used on the spot to produce an estimate of conservation within orthologous sets, or functional divergence across paralogues.

Regulation of alternative splicing by the circadian clock and food related cues

Background

The circadian clock orchestrates daily rhythms in metabolism, physiology and behaviour that allow organisms to anticipate regular changes in their environment, increasing their adaptation. Such circadian phenotypes are underpinned by daily rhythms in gene expression. Little is known, however, about the contribution of post-transcriptional processes, particularly alternative splicing.

Results

Using Affymetrix mouse exon-arrays, we identified exons with circadian alternative splicing in the liver. Validated circadian exons were regulated in a tissue-dependent manner and were present in genes with circadian transcript abundance. Furthermore, an analysis of circadian mutant Vipr2^-/- mice revealed the existence of distinct physiological pathways controlling circadian alternative splicing and RNA binding protein expression, with contrasting dependence on Vipr2-mediated physiological signals. This view was corroborated by the analysis of the effect of fasting on circadian alternative splicing. Feeding is an important circadian stimulus, and we found that fasting both modulates hepatic circadian alternative splicing in an exon-dependent manner and changes the temporal relationship with transcript-level expression.

Conclusions

The circadian clock regulates alternative splicing in a manner that is both tissue-dependent and concurrent with circadian transcript abundance. This adds a novel temporal dimension to the regulation of mammalian alternative splicing. Moreover, our results demonstrate that circadian alternative splicing is regulated by the interaction between distinct physiological cues, and illustrates the capability of single genes to integrate circadian signals at different levels of regulation.

The dynamics of the alternatively spliced NOL7 gene products and role in nucleolar architecture

Three alternatively spliced forms of the human NOL7 gene coding for relatively small proteins were identified. The two shorter forms were generated by intron retention events, and each isoform was differently localized within the cell. The NOL7-SP1 long form (29 kD) localized to the nucleolus, SP2 was nucleoplasmic, while SP3 was distributed throughout the whole cell. NOL7-SP1 was confined to the nucleolar granular component, and during cell division disassociated from the nucleolus. Knockdown of NOL7-SP1 levels abrogated nucleolar architecture, in particular the internal regions, and reduced cell proliferation. Analysis of the nucleolar dynamics of the SP1 protein during interphase showed nucleolar high binding affinity. Dissection of protein domains showed that nucleolar targeting was mediated by a unique C-terminal nucleolar localization sequence (NoLS). However, this sequence was not sufficient for conferring high binding affinity, which required additional regions of the protein. Our analysis shows that NOL7 is important for maintaining internal nucleolar structure and cell growth rates, and that while specific protein localization can be obtained by specific short localization motifs, nucleolar residency through binding must be mediated by a synergistic combination of protein modules.

The Tpr protein regulates export of mRNAs with retained introns that traffic through the Nxf1 pathway

Post-transcriptional regulation of mRNA includes restriction mechanisms to prevent export and expression of mRNAs that are incompletely spliced. Here we present evidence that the mammalian protein Tpr is involved in this restriction. To study the role of Tpr in export of mRNA with retained introns, we used reporters in which the mRNA was exported either via the Nxf1/Nxt1 pathway using a CTE or via the Crm1 pathway using Rev/RRE. Our data show that even modest knockdown of Tpr using RNAi leads to a significant increase in export and translation from the mRNA containing the CTE. In contrast, Tpr perturbation has no effect on export of mRNA containing the RRE, either in the absence or presence of Rev. Also, no effects were observed on export of a completely spliced mRNA. Taken together, our results indicate that Tpr plays an important role in quality control of mRNA trafficked on the Nxf1 pathway.

Dynamic usage of alternative splicing exons during mouse retina development

Alternative processing of pre-mRNA plays an important role in protein diversity and biological function. Previous studies on alternative splicing (AS) often focused on the spatial patterns of protein isoforms across different tissues. Here we studied dynamic usage of AS across time, during murine retina development. Over 7000 exons showed dynamical changes in splicing, with differential splicing events occurring more frequently in early development. The overall splicing patterns for exclusive and inclusive exons show symmetric trends and genes with symmetric splicing patterns that tend to have similar biological functions. Furthermore, we observed that within the retina, retina-enriched genes that are preferentially expressed at the adult stage tend to have more dynamically spliced exons compared to other genes, suggesting that genes maintaining retina homeostasis also play an important role in development via a series of AS events. Interestingly, the transcriptomes of retina-enriched genes largely reflect the retinal developmental process. Finally, we identified a number of candidate cis-regulatory elements for retinal AS by analyzing the relative occurrence of sequence motifs in exons or flanking introns. The occurrence of predicted regulatory elements showed strong correlation with the expression level of known RNA binding proteins, suggesting the high quality of the identified cis-regulatory elements.

A dual platform approach to transcript discovery for the planarian Schmidtea mediterranea to establish RNAseq for stem cell and regeneration biology

The use of planarians as a model system is expanding and the mechanisms that control planarian regeneration are being elucidated. The planarian Schmidtea mediterranea in particular has become a species of choice. Currently the planarian research community has access to this whole genome sequencing project and over 70,000 expressed sequence tags. However, the establishment of massively parallel sequencing technologies has provided the opportunity to define genetic content, and in particular transcriptomes, in unprecedented detail. Here we apply this approach to the planarian model system. We have sequenced, mapped and assembled 581,365 long and 507,719,814 short reads from RNA of intact and mixed stages of the first 7 days of planarian regeneration. We used an iterative mapping approach to identify and define de novo splice sites with short reads and increase confidence in our transcript predictions. We more than double the number of transcripts currently defined by publicly available ESTs, resulting in a collection of 25,053 transcripts described by combining platforms. We also demonstrate the utility of this collection for an RNAseq approach to identify potential transcripts that are enriched in neoblast stem cells and their progeny by comparing transcriptome wide expression levels between irradiated and intact planarians. Our experiments have defined an extensive planarian transcriptome that can be used as a template for RNAseq and can also help to annotate the S. mediterranea genome. We anticipate that suites of other 'omic approaches will also be facilitated by building on this comprehensive data set including RNAseq across many planarian regenerative stages, scenarios, tissues and phenotypes generated by RNAi.

Steroid receptor RNA activator protein binds to and counteracts SRA RNA-mediated activation of MyoD and muscle differentiation

The steroid receptor RNA activator (SRA) has the unusual property to function as both a non-coding RNA (ncRNA) and a protein SRAP. SRA ncRNA is known to increase the activity of a range of nuclear receptors as well as the master regulator of muscle differentiation MyoD. The contribution of SRA to either a ncRNA or a protein is influenced by alternative splicing of the first intron, the retention of which disrupts the SRAP open reading frame. We reported here that the ratio between non-coding and coding SRA isoforms increased during myogenic differentiation of human satellite cells but not myotonic dystrophy patient satellite cells, in which differentiation capacity is affected. Using constructs that exclusively produce SRA ncRNA or SRAP, we demonstrated that whereas SRA ncRNA was indeed an enhancer of myogenic differentiation and myogenic conversion of non-muscle cells through the co-activation of MyoD activity, SRAP prevented this SRA RNA-dependant co-activation. Interestingly, the SRAP inhibitory effect is mediated through the interaction of SRAP with its RNA counterpart via its RRM-like domain interacting with the functional sub-structure of SRA RNA, STR7. This study thus provides a new model for SRA-mediated regulation of MyoD transcriptional activity in the promotion of normal muscle differentiation, which takes into account the nature of SRA molecules present.