This is a test version of Biostars. For the public version, visit https://www.biostars.org.
A Basic question regarding lncRNA identification pipeline.

Hi,

I have been analyzing RNA-Seq data sets of some Breast cancer cell lines to create a high confidence list of expressed lncRNAs. However as, I am new to NGS, I cannot figure out how do I filter out the known Expressed gene/protein coding transcripts from my annotation file after cufflinks assembly? Are there any specific tools to do the filtering? If anyone could help me regarding this, I will really appreciate your help.

Thanks

R

rna-seq lncrna pipeline
• 2,482 views • link
written by rajeev.vikram • rajeev.vikram

2 answers

http://cole-trapnell-lab.github.io/cufflinks/cuffcompare/

• 0 views • link
written by GouthamAtla GouthamAtla

Thanks, but my question is slightly different, , basically, after top-hat assembly with bowtie2 , I Used RABT assembly in cufflinks and then merged all transcripts (elegant= gtf file of annotated transcripts), then did cuffmerge of the replicates. After running cuffcompare with r- given as annotated gencode assembly, I got the transfrags identified with diff signs (=, c x etc.) Now I want to filter out all transfrags of ‘i’, ‘j’, ‘o’, ‘u’ and ‘x’ option, while making an extra file of known lncRNAs (by matching with bodymap annotated lncRNA.gtf). I am curious if I can do all that in command line in one comment, something like:

awk ‘$22 ~ /j,i,o,u,x/ { print }’..
• 0 views • link
written by rajeev.vikram • rajeev.vikram

Log in to answer this question.

More posts like this