I downloaded RNASeq bam files from TCGA data portal. I am interested in the expression analysis at isofrom level and wondering how to translate these bam files into counts at each isoform. I have done some analysis at gene level using star or htseq-count but fresh to isoform analysis. Could somebody have me some suggestions. I knew there are softwares such as rMATs for alternative splicing (AS) analysis but what I wanted are the counts for each isoform rather than the AS analysis.
•
written
by
ddzhangzz
2 answers
There is an example (case study) in the edgeR user's guide. Chapter 4.5. They use featureCounts first instead of htseq.
https://bioconductor.org/packages/release/bioc/vignettes/edgeR/inst/doc/edgeRUsersGuide.pdf
Try sailfish They are not the exact counts but would work for many purposes.
Log in to answer this question.
-
Isoform analysis after quantification with Salmon/Star
written by rajdeepboral00 •
-
is there any way to normalize the counts and compare it for the splicing quantification?
written by PK •
-
Tools for Alternative Splicing Events in RNA-Seq analysis
written by joe_genome •
-
Software/packages for analyzing differential isoforms from RNA-seq data?
written by Kai_Qi
-
Workflow/pipeline for detecting isoforms/alternative splicing
written by Aynur •
-
Methods to convert SOLiD reads to nucleotide space
written by gsc529 •
-
same query but different results (TCGAbiolinks vs TCGA Portal)
written by tyasird •
-
Alternative splicing analysis of targeted RNASeq data
written by ddzhangzz
-
Alternative splicing using rMATS
written by Bionic •
-
How the level 3 TCGA raw counts are calculated for each sample?
written by pbio
You would want to carefully check how multi-mappers were handled if you are starting with BAM files and considering isoforms.