low assignment reads % in RNA-seq
2
0
Entering edit mode
5 months ago
parinv ▴ 30

I performed RNA-seq analysis on 8 samples(2 control and 6 samples) . I used featureCounts function (Rsubread package) to count the number of reads. I found that Number of assigned reads were very poor. I am attaching the table here for your reference. I used following function:

fc<- featureCounts(annot$Sample,isPairedEnd=TRUE, GTF.attrType = "transcript_id",GTF.featureType="exon", annot.ext = "GCF_000203855.3_ASM20385v3_genomic.gtf", nthreads = 3, isGTFAnnotationFile = TRUE)

I do not understand why assignment% is low, when alignment % is found higher. Can anyone help me think of possible reasons for why I'm getting such low values for assignment?

Thanks a lot. Parin

rna-seq R • 337 views
ADD COMMENT
1
Entering edit mode

check if multi-mapped reads are counted or not and use correct strandedness. Talk to the core and check if ribodepletion worked for all the samples.

ADD REPLY
0
Entering edit mode

can you post the alignment summary metrics for example like this

ADD REPLY
0
Entering edit mode

I figured out that using GTF.featureType="exon" is giving low assignment %, If I use GTF.featureType="CDS", the assignment % is increasing as shown here. I cannot decide which feature type should I use. I read some answers 1 2 and found that any one can be used. Kindly share some information regarding use of exon or CDS.

Thankyou

ADD REPLY
0
Entering edit mode

This is subjective - pertaining to questions asked in the experiment

Still confused about exons versus CDS

ADD REPLY
2
Entering edit mode
5 months ago

One possible reason is that you may need to indicate the strand information of your library. Besides, this issue has been discussed in multiple threads, like here and here. Those threads and others should contain suggestions that can help you resolve the issue. Good luck.

ADD COMMENT
0
Entering edit mode
5 months ago
sysboolean ▴ 80

Can you post more information about the origin of the samples ? Someone I know once performed RNA-seq on what they thought was human cell culture, but got poor alignment percents like yours. Further investigation (BLASTing a few reads against nr db) showed that the cells were from a non-human mammal. Might help to run Fastqscreen to make sure your samples were not switched during sequencing.

ADD COMMENT

Login before adding your answer.

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