3D RNA-seq: a powerful and flexible tool for rapid and accurate differential expression and alternative splicing analysis of RNA-seq data for biologists.

RNA-sequencing (RNA-seq) analysis of gene expression and alternative splicing should be routine and robust but is often a bottleneck for biologists because of different and complex analysis programs and reliance on specialized bioinformatics skills. We have developed the '3D RNA-seq' App, an R shiny App and web-based pipeline for the comprehensive analysis of RNA-seq data from any organism. It represents an easy-to-use, flexible and powerful tool for analysis of both gene and transcript-level gene expression to identify differential gene/transcript expression, differential alternative splicing and differential transcript usage (3D) as well as isoform switching from RNA-seq data. 3D RNA-seq integrates state-of-the-art differential expression analysis tools and adopts best practice for RNA-seq analysis. The program is designed to be run by biologists with minimal bioinformatics experience (or by bioinformaticians) allowing lab scientists to analyse their RNA-seq data. It achieves this by operating through a user-friendly graphical interface which automates the data flow through the programs in the pipeline. The comprehensive analysis performed by 3D RNA-seq is extremely rapid and accurate, can handle complex experimental designs, allows user setting of statistical parameters, visualizes the results through graphics and tables, and generates publication quality figures such as heat-maps, expression profiles and GO enrichment plots. The utility of 3D RNA-seq is illustrated by analysis of data from a time-series of cold-treated Arabidopsis plants and from dexamethasone-treated male and female mouse cortex and hypothalamus data identifying dexamethasone-induced sex- and brain region-specific differential gene expression and alternative splicing.

Evolutionary relationships among barley and Arabidopsis core circadian clock and clock-associated genes.

The circadian clock regulates a multitude of plant developmental and metabolic processes. In crop species, it contributes significantly to plant performance and productivity and to the adaptation and geographical range over which crops can be grown. To understand the clock in barley and how it relates to the components in the Arabidopsis thaliana clock, we have performed a systematic analysis of core circadian clock and clock-associated genes in barley, Arabidopsis and another eight species including tomato, potato, a range of monocotyledonous species and the moss, Physcomitrella patens. We have identified orthologues and paralogues of Arabidopsis genes which are conserved in all species, monocot/dicot differences, species-specific differences and variation in gene copy number (e.g. gene duplications among the various species). We propose that the common ancestor of barley and Arabidopsis had two-thirds of the key clock components identified in Arabidopsis prior to the separation of the monocot/dicot groups. After this separation, multiple independent gene duplication events took place in both monocot and dicot ancestors.

An hnRNP-like RNA-binding protein affects alternative splicing by in vivo interaction with transcripts in Arabidopsis thaliana.

Alternative splicing (AS) of pre-mRNAs is an important regulatory mechanism shaping the transcriptome. In plants, only few RNA-binding proteins are known to affect AS. Here, we show that the glycine-rich RNA-binding protein AtGRP7 influences AS in Arabidopsis thaliana. Using a high-resolution RT-PCR-based AS panel, we found significant changes in the ratios of AS isoforms for 59 of 288 analyzed AS events upon ectopic AtGRP7 expression. In particular, AtGRP7 affected the choice of alternative 5' splice sites preferentially. About half of the events are also influenced by the paralog AtGRP8, indicating that AtGRP7 and AtGRP8 share a network of downstream targets. For 10 events, the AS patterns were altered in opposite directions in plants with elevated AtGRP7 level or lacking AtGRP7. Importantly, RNA immunoprecipitation from plant extracts showed that several transcripts are bound by AtGRP7 in vivo and indeed represent direct targets. Furthermore, the effect of AtGRP7 on these AS events was abrogated by mutation of a single arginine that is required for its RNA-binding activity. This indicates that AtGRP7 impacts AS of these transcripts via direct interaction. As several of the AS events are also controlled by other splicing regulators, our data begin to provide insights into an AS network in Arabidopsis.

Dual functionality of a plant U-rich intronic sequence element.

In potato invertase genes, the constitutively included, 9-nucleotide (nt)-long mini-exon requires a strong branchpoint and U-rich polypyrimidine tract for inclusion. The strength of these splicing signals was demonstrated by greatly enhanced splicing of a poorly spliced intron and by their ability to support splicing of an artificial mini-exon, following their introduction. Plant introns also require a second splicing signal, UA-rich intronic elements, for efficient intron splicing. Mutation of the branchpoint caused loss of mini-exon inclusion without loss of splicing enhancement, showing that the same U-rich sequence can function as either a polypyrimidine tract or a UA-rich intronic element. The distinction between the splicing signals depended on intron context (the presence or absence of an upstream, adjacent and functional branchpoint), and on the sequence context of the U-rich elements. Polypyrimidine tracts tolerated C residues while UA-rich intronic elements tolerated As. Thus, in plant introns, U-rich splicing elements can have dual roles as either a general plant U-rich splicing signal or a polypyrimidine tract. Finally, overexpression of two different U-rich binding proteins enhanced intron recognition significantly. These results highlight the importance of co-operation between splicing signals, the importance of other nucleotides within U-rich elements for optimal binding of competing splicing factors and effects on splicing efficiency of U-rich binding proteins.

Mutational analysis of a plant branchpoint and polypyrimidine tract required for constitutive splicing of a mini-exon.

The branchpoint sequence and associated polypyrimidine tract are firmly established splicing signals in vertebrates. In plants, however, these signals have not been characterized in detail. The potato invertase mini-exon 2 (9 nt) requires a branchpoint sequence positioned around 50 nt upstream of the 5' splice site of the neighboring intron and a U11 element found adjacent to the branchpoint in the upstream intron (Simpson et al., RNA, 2000, 6:422-433). Utilizing the sensitivity of this plant splicing system, these elements have been characterized by systematic mutation and analysis of the effect on inclusion of the mini-exon. Mutation of the branchpoint sequence in all possible positions demonstrated that branchpoints matching the consensus, CURAY, were most efficient at supporting splicing. Branchpoint sequences that differed from this consensus were still able to permit mini-exon inclusion but at greatly reduced levels. Mutation of the downstream U11 element suggested that it functioned as a polypyrimidine tract rather than a UA-rich element, common to plant introns. The minimum sequence requirement of the polypyrimidine tract for efficient splicing was two closely positioned groups of uridines 3-4 nt long (<6 nt apart) that, within the context of the mini-exon system, required being close (<14 nt) to the branchpoint sequence. The functional characterization of the branchpoint sequence and polypyrimidine tract defines these sequences in plants for the first time, and firmly establishes polypyrimidine tracts as important signals in splicing of at least some plant introns.