lncRNA annotation for featureCouns and STAR
1
0
Entering edit mode
inah • 0
@inah-13176
Last seen 3.5 years ago

Hi,

  I have total RNAseq data and one of my main goals is to map lncRNAs (long noncoding RNAs). I am using STAR for alignment and featureCounts for mapping (counting). For the lncRNAs, I think I need to use the annotation file gencode.v28.long_noncoding_RNAs.gtf.gz.  I will pass this annotation file on to the counter featureCounts.

 

I have already created a Genome Index for STAR but using the different annotation file, with --genomeFastaFiles /home/inah/RefGTF/GRCh38/sequence/Homo_sapiens.GRCh38.dna.primary_assembly.fa  --sjdbGTFfile /home/inah/RefGTF/GRCh38/annotation/Homo_sapiens.GRCh38.85.gtf .

 

My question is: Do I need to rerun STAR to create a Genome Index using the lncRNA annotation file gencode.v28.long_noncoding_RNAs.gtf.gz and then rerun STAR for the alignment before I run featureCounts with the lncRNA annotation file?   Or can I just run featureCounts with the lncRNA annotation file using the STAR alignments obtained with the other annotation file?  In other words, do I need to use the same annotation file for the aligner STAR and the counter featureCounts?

 

Thanks, Ina

 

 

lncRNA featureCounts STAR annotation • 1.2k views
ADD COMMENT
0
Entering edit mode

Hi, I am having similar problem. Did you find any solution?

ADD REPLY
0
Entering edit mode
Wei Shi ★ 3.3k
@wei-shi-2183
Last seen 4 days ago
Australia/Melbourne/Olivia Newton-John …

Aligners like STAR, TopHat and Subread use an annotation to assist the mapping of RNA-seq reads. Since your data is total RNA sequencing data, it makes more sense to me to use a transcriptome annotation in your read mapping instead of a lncRNA annotation. You can then count mapped reads to lncRNAs using featureCounts.

Technically featureCounts works with any alignment results output from the aligner no matter which annotation was used in the mapping (or no annotation was used).

ADD COMMENT

Login before adding your answer.

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