featureCounts unable to assign long reads to one specific isoform
1
0
Entering edit mode
Alex • 0
@ffb4b031
Last seen 8 weeks ago
United States

Hello,

I'm trying to use featureCounts to annotate a long read dataset obtained from nanopore sequencing. I noticed from the results reads that are supposed to be assigned unambiguously to a specific isoform, were actually assigned all different isoforms, clearly ignoring some isoform-specific splice junctions. One example is shown below, where all the reads were assigned to all 3 ACTG1 isoforms.

enter image description here

I was running featureCounts v2.0.6 in a Conda environment with the following command:

featureCounts -t exon -g transcript_id -O -R CORE -L -a ref.gtf -o counts.tab aln.bam

I was wondering what I can do about it, and I can send the dataset if needed.

Thank you!

featureCounts Rsubread • 230 views
ADD COMMENT
1
Entering edit mode
Wei Shi ★ 3.6k
@wei-shi-2183
Last seen 4 days ago
Australia/Melbourne/Olivia Newton-John …

The -O option in your command instructs featureCounts to assign a read to all the transcripts that it overlaps by one or more bases. That is why your reads were assigned to all the isoforms. featureCounts does not consider isoform-specific splice junctions when it performs read assignment. It just performs read assignment based on the overlap between reads and targets.

ADD COMMENT

Login before adding your answer.

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