This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Exclude non-coding genes

Hello, I have HT-Seq read count with around 60000 genes (with ENSEMBLE ID). I only want to keep protein coding genes. How can I remove non-coding genes?

Thanks

rna-seq gene
• 2,801 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by Rob • Rob

1 answer

go to http://www.ensembl.org/biomart/martview

database: gene / human gene

attribute : gene stable id

gene: -> gene_type "protein coding"

export as a text file

filter your list with grep -w -f coding.txt you_ids.txt

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thanks a lot Pierre It was awesome and so easy. I just do not understand your last point: "filter your list with grep -w -f coding.txt you_ids.txt"

how should I do this? sorry I am new to this field and not familiar with things and need more explanation.

• 0 views • link
written by Rob • Rob

You can do this through Terminal if you are on a Mac/Linux as grep is used for Unix operating systems.

  1. First cd to the folder/directory your files (the protein coding genes and your list containing the 6000 genes) are in:

    cd /Users/your/directory/

  2. Run the grep command to filter your your_ids.txt (6000 genes file) using coding.txt (list from Ensembl).

    grep -w -f coding.txt your_ids.txt

  3. This will output to terminal, you can save the output to a file with > output.txt

    grep -w -f coding.txt your_ids.txt > output.txt

• 0 views • link
written by Tota • Tota

Hello I am using windows system. Is there any method for wimdows?

• 0 views • link
written by Rob • Rob

Hello, I'm not familiar with Windows. There is the option to install WSL (Windows Subsystem for Linux) which lets you run the above commands on your Windows machine. Perhaps someone else can suggest a better alternative.

• 0 views • link
written by Tota • Tota

Log in to answer this question.

More posts like this