How to use RSEM to quantify genes and transcript expression from STAR outpute file (Aligned.toTranscriptome.out.bam)
0
0
Entering edit mode
6 weeks ago
otieno43 • 0

Have already run normal Star alignment and obtained Aligned.toTranscriptome.out.bam file which should be compatible with RSEM tool. However, when I use RSEM to quantify genes and transcript expression I do not any result. I think I do not get what one part rsem-calculate-expression.

If this is the command:

rsem-calculate-expression --bam --no-bam-output -p 10 \
--paired-end --forward-prob 0 \
RNASEQ_data/star_GM12878_rep1/Aligned.toTranscriptome.out.bam # (star output file)
GENOME_data/rsem/GRCh38 \ #(rsem reference genome)
RNASEQ_data/rsem_GM12878_rep1/rsem \ # What is this suppose to be? I think this is what I am missing)
>& RNASEQ_data/rsem_GM12878_rep1/rsem.log

Any help please.

Erick

RSEM STAR • 397 views
ADD COMMENT
0
Entering edit mode

The "what is this supposed to be" part is the output prefix. Are you using >& when you mean to use &>? In any case, use >rsem.out 2>rsem.err and show us the contents of both files.

ADD REPLY
0
Entering edit mode

Hello Ram,

I have tried what you suggested, I still got an error. The error reads, "-tag XM" failed! Plase check if you provide correct parameters/options for the pipeline!"

rsem-calculate-expression --bam --no-bam-output -p 10 \
--paired-end --forward-prob 0 \
RNASEQ_data/star_GM12878_rep1/Aligned.toTranscriptome.out.bam \
GENOME_data/rsem/GRCh38 \
RNASEQ_data/rsem_GM12878_rep1/rsem \
> RNASEQ_data/rsem_GM12878_rep1/rsem.log 2>RNASEQ_data/rsem_GM12878_rep1/rsem.err
ADD REPLY
0
Entering edit mode

Another error is found in rsem.err file reads, "The SAM/BAM file declares more reference sequences (25014) than RSEM knows (20519)!"

This seems strange because I used the same reference genome (fasta) and annotation (gft) file to generate rsem reference genome and star index. I used stand alone star to build index same to rsem ref genome.

ADD REPLY
0
Entering edit mode

Can you share your STAR command as well? See the STAR command line used by RSEM here: https://github.com/deweylab/RSEM/blob/master/rsem-calculate-expression#L457

You will need to dig deeper to understand why your STAR generated BAM was not compatible with RSEM. I also think that your STAR and RSEM are not being run against the same prepared index or at least not with compatible parameters.

ADD REPLY

Login before adding your answer.

Traffic: 1847 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