Question: Ribo-seq vs RNA-seq read count
1
gravatar for ashkan
13 months ago by
ashkan100
ashkan100 wrote:

Hi Guys,

I have RNA-seq and Ribo-seq data. I have align them to the human genome and got the read count per gene. in my output file the number of read count/gene in RNA-seq is lower than Ribo-seq. my question is that should it be higher in RNA-seq or my results are correct? thanks

rna-seq • 803 views
ADD COMMENTlink modified 13 months ago by Devon Ryan73k • written 13 months ago by ashkan100

How many reads per experiment? I mean how deep did you sequence?

I can imagine that RNA-seq would sequence more nt than Ribo-seq, thus with same number of reads you would have more reads per gene with Ribo-seq. But of course it also depends on how many genes were expressed in you RNA-seq experiment...

ADD REPLYlink written 13 months ago by b.nota3.6k

Depends on the sequencing depth I guess. Could you edit your question to add some informations (how did you count the reads, which library prep, sequencing depth, etc...)

ADD REPLYlink modified 13 months ago • written 13 months ago by Nicolas Rosewick5.7k

to prepare the library our sequencing facility used this: TruSeq RNA Library Prep Kit v2 to count the number of reads I used the following command:

python -m HTSeq.scripts.count -f -h bam -r pos -m intersection-nonempty -s no -a 10 accepted_hits.bam gencocd de.v19.annotation.gtf > 285_Ctrl_counts.txt

ADD REPLYlink written 13 months ago by ashkan100
3
gravatar for Devon Ryan
13 months ago by
Devon Ryan73k
Freiburg, Germany
Devon Ryan73k wrote:

You're measuring different things with the two datasets. RNAseq is measuring the amount of RNA present. RiboSeq is measuring the amount of rRNA bound RNA present. The latter is showing either active translation or pausing, the former steady state transcription. For a given gene, there's no a priori reason to expect higher lower relative signal in one dataset versus the other.

Of course, if you're just comparing raw numbers without accounting for the number of alignments or anything else then that's a different issue...

ADD COMMENTlink modified 13 months ago • written 13 months ago by Devon Ryan73k

Hi Davon, to analyse Ribo-seq data I filtered out rRNA, tRNA, adaptors and markers and then align them to the reference genome but for the RNA-seq I aligned to the reference genome without any filtering

ADD REPLYlink written 13 months ago by ashkan100

You might also need to get rid of snRNAs from the RiboSeq dataset. Are you then seeing huge differences between the resulting percentages being included in the counts or between the raw counts? If the latter, how many alignments are going into htseq-count in each?

ADD REPLYlink written 13 months ago by Devon Ryan73k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1380 users visited in the last hour