Dataset Information

Identification and visualization of differential isoform expression in RNA-seq time series.

ABSTRACT: Motivation:As sequencing technologies improve their capacity to detect distinct transcripts of the same gene and to address complex experimental designs such as longitudinal studies, there is a need to develop statistical methods for the analysis of isoform expression changes in time series data. Results:Iso-maSigPro is a new functionality of the R package maSigPro for transcriptomics time series data analysis. Iso-maSigPro identifies genes with a differential isoform usage across time. The package also includes new clustering and visualization functions that allow grouping of genes with similar expression patterns at the isoform level, as well as those genes with a shift in major expressed isoform. Availability and implementation:The package is freely available under the LGPL license from the Bioconductor web site. Contact:mj.nueda@ua.es or aconesa@ufl.edu. Supplementary information:Supplementary data are available at Bioinformatics online.

SUBMITTER: Nueda MJ

PROVIDER: S-EPMC5860359 | biostudies-literature | 2018 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Identification and visualization of differential isoform expression in RNA-seq time series.

Nueda María José MJ Martorell-Marugan Jordi J Martí Cristina C Tarazona Sonia S Conesa Ana A

Bioinformatics (Oxford, England) 20180201 3

<h4>Motivation</h4>As sequencing technologies improve their capacity to detect distinct transcripts of the same gene and to address complex experimental designs such as longitudinal studies, there is a need to develop statistical methods for the analysis of isoform expression changes in time series data.<h4>Results</h4>Iso-maSigPro is a new functionality of the R package maSigPro for transcriptomics time series data analysis. Iso-maSigPro identifies genes with a differential isoform usage across ...[more]

PMID: 28968682

Similar Datasets

Project description:BACKGROUND: RNA-seq is now widely used to quantitatively assess gene expression, expression differences and isoform switching, and promises to deliver results for the entire transcriptome. However, whether the transcriptional state of a gene can be captured accurately depends critically on library preparation, read alignment, expression estimation and the tests for differential expression and isoform switching. There are comparisons available for the individual steps but there is not yet a systematic investigation which specific genes are impacted by biases throughout the entire analysis workflow. It is especially unclear whether for a given gene, with current methods and protocols, expression changes and isoform switches can be detected. RESULTS: For the human genes, we report their detectability under various conditions using different approaches. Overall, we find that the input material has the biggest influence and may, depending on the protocol and RNA degradation, exhibit already strong length-dependent over- and underrepresentation of transcripts. The alignment step aligns for 50% of the isoforms up to 99% of the reads correctly; only in the presence of transcript modifications mainly short isoforms will have a low alignment rate. In our dataset, we found that, depending on the aligner and the input material used, the expression estimation of up to 93% of the genes being accurate within a factor of two; with the deviations being due to ambiguous alignments. Detection of differential expression using a negative-binomial count model works reliably for our simulated data but is dependent on the count accuracy. Interestingly, using the fold-change instead of the p-value as a score for differential expression yields the same performance in the situation of three replicates and the true change being two-fold. Isoform switching is harder to detect and for at least 109 genes the isoform differences evade detection independent of the method used. CONCLUSIONS: RNA-seq is a reliable tool but the repetitive nature of the human genome makes the origin of the reads ambiguous and limits the detectability for certain genes. RNA-seq does not equally well represent isoforms independent of their size which may range from ~200nt to ~100'000nt. Researchers are advised to verify that their target genes do not have extreme properties with respect to repeated regions, GC content, and isoform length and complexity.

Project description:Graphical abstract The cGAS-STING pathway, orchestrating complicated transcriptome-wide immune responses, is essential for host antiviral defense but can also drive immunopathology in severe COVID-19. Here, we performed time-course RNA-Seq experiments to dissect the transcriptome expression dynamics at the gene-isoform level after cGAS-STING pathway activation. The in-depth time-course transcriptome after cGAS-STING pathway activation within 12 h enabled quantification of 48,685 gene isoforms. By employing regression models, we obtained 13,232 gene isoforms with expression patterns significantly associated with the process of cGAS-STING pathway activation, which were named activation-associated isoforms. The combination of hierarchical and k-means clustering algorithms revealed four major expression patterns of activation-associated isoforms, including two clusters with increased expression patterns enriched in cell cycle, autophagy, antiviral innate-immune functions, and COVID-19 coronavirus disease pathway, and two clusters showing decreased expression pattern that mainly involved in ncRNA metabolism, translation process, and mRNA processing. Importantly, by merging four clusters of activation-associated isoforms, we identified three types of genes that underwent isoform usage alteration during the cGAS-STING pathway activation. We further found that genes exhibiting protein-coding and non-protein-coding gene isoform usage alteration were strongly enriched for the factors involved in innate immunity and RNA splicing. Notably, overexpression of an enriched splicing factor, EFTUD2, shifted transcriptome towards the cGAS-STING pathway activated status and promoted protein-coding isoform abundance of several key regulators of the cGAS-STING pathway. Taken together, our results revealed the isoform-level gene expression dynamics of the cGAS-STING pathway and uncovered novel roles of splicing factors in regulating cGAS-STING pathway mediated immune responses.

Dataset Information

Identification and visualization of differential isoform expression in RNA-seq time series.

Publications

Identification and visualization of differential isoform expression in RNA-seq time series.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets