How to quantify mRNA transcripts from RNA-Seq data?
3
0
Entering edit mode
5.3 years ago
Neu ▴ 10

Hi, I am analyzing total transcriptome RNA-Seq data. I want to quantify expression of mRNA transcripts. What should be the correct approach. I have aligned the reads with reference genome (sequence file -Mus_musculus.GRCm38.dna.primary_assembly.fa and GTF file - Mus_musculus.GRCm38.95.chr.gtf) using STAR. Please guide me to do the quantification.

Thank you!

RNA-Seq next-gen alignment • 1.7k views
ADD COMMENT
1
Entering edit mode
5.3 years ago
Buffo ★ 2.4k

There are some biotools to quantify gene expression or transcripts, such as ht-seq, or strintie ( from the new protocolo tuxedo). However I recommend you to read about RNAseq-analysis protocols before to continue, in order to know "what" and "why" are you doing.

ADD COMMENT
1
Entering edit mode
5.3 years ago

Gene expression can be quantified either by pseudo-aligners such as kallisto and salmon, which Kristoffer mentioned, or, following read alignment as you just did, using tools such as featureCounts or HTSeq-Counts.

The main difference between the pseudo-aligners and the counting tools is that the counting tools absolutely need the BAM files and there is some evidence that suggests that the results may be more accurate for complex and well-annotated transcriptomes with the pseudo-aligners (see this paper for more details).

The main difference between featureCounts and HTSeq-Counts is speed (featureCounts is way faster).

For more details on how to proceed after the quantification, see Michael Love's tutorial.

ADD COMMENT
0
Entering edit mode
5.3 years ago

I wrote a long section about the different things to think about when quantifying RNASeq data (including tools to use) here.

ADD COMMENT

Login before adding your answer.

Traffic: 1810 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6