How to Annotate Long non coding RNA?
2
0
Entering edit mode
5.2 years ago
Hyper_Odin ▴ 320

I am a bit confused about the annotation of Lnc-RNA. I have already downloaded the gencode gtf file and I i am using cufflinks, cuffmerge and cuffdiff. Aftercufflinks i get "transcripts.gtf" and at the end of cuffmerge i get "merged.gtf" Now i am not able to understand the next step. 1) At which file should I apply the cutoff filter of 200bp and then go for CPAT? is it the "transcript.gtf"? 2) If yes, then what will be the command to execute it?

I am reading several papers, but they have not described the procedure clearly! Thanks

RNA-Seq sequencing alignment • 1.1k views
ADD COMMENT
1
Entering edit mode
5.2 years ago

CPAT takes a fasta file in input and gives you an information if a transcript is coding or not . So you should first select the transcript of interest from your gtf then extract the sequence from the genome. Then CPAT

ADD COMMENT
1
Entering edit mode
5.2 years ago
27atcggcta27 ▴ 180

If I remember correctly, cuffmerge gives you a nonredundant list of transcripts between all your samples for you to use as the annotation later on. So I would apply the 200 nc cutoff here (merged.gtf) and then use CPAT. You could also use CNCI or CPC.

ADD COMMENT

Login before adding your answer.

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