This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Filter out data for coding genes

Hello all

I have obtained RNA-seq data from TCGA. The data contains about 56,000 genes (gene symbols and ensemble ids) that contains pseudo, uncharacterized, non-coding and coding genes. I am looking for the efficient way to filter out the data for only coding genes. Can anyone suggest?

Regards Shivangi Agarwal

coding rna-seq
• 2,603 views • link
written by shivangi.agarwal800 shivangi.agarwal800

In what format do you have the data?

• 1 views • link
written by ATpoint ATpoint

I have patient'swise data which contains expression value (FPKM) for each of the gene in each of the patient. (FPKM.txt files)

• 1 views • link
written by shivangi.agarwal800 shivangi.agarwal800

Then I would get a GTF file for human of the correct assembly (hg19, hg38, don't know what TCGA used, e.g. from GENCODE), and then filter for genes that are coding (CDS in the GTF). From this, extract the gene names and use these to subset your data. This all can be done with Unix tools like awk.

• 1 views • link
written by ATpoint ATpoint

0 answers

No answers yet.

Log in to answer this question.

More posts like this