This is a test version of Biostars. For the public version, visit https://www.biostars.org.
A script for extracting information related to a list of gene names from a file

Dear Biostars,

I have a text file containing several rows and columns like this:

"Gene Name"      "Gene Id"       "description"         "GO"
A                1               phosphatase           GO:001256
B                2               synthesize            GO:013154
C                3               methylase             GO:000054
D                4               kinase                GO:001254
E                5               oxigenase             GO:001354
F                6               synthesize            GO:001254

In addition, I have another text file just containing one column and several rows like this:

Gene Name
A
D
C
B

I need to extract the rows of file 1 that contain gene names listed in file 2.

Does anybody have any idea how to do that?

PS: I know how to do that by excel but it does not work with huge rows of information.

Thank you

rna-seq data-mining gene
• 8,287 views • link
updated by Ram Ram • written by ahmad.iut ahmad.iut

2 answers

using linux:

join -1 1 -2 1 <(sort -k1,1 file1.txt) <(sort -k1,1 file2.txt)  > joined.txt

using knime.org:

load both files (Read File) in two tables and join using a "Join" node https://www.knime.org/files/nodedetails/_manipulation_column_column_split_combine_Joiner.html

• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Dear Lindenbaum,

The command worked perfectly. Thank you very much

• 0 views • link
written by ahmad.iut ahmad.iut

thanks a lot, this is very useful also for my problem. Just a question, is it possible to include the header as well? adding --header is not working.

• 0 views • link
written by User000 User000

adding --header is not working.

should work. check the input files order, check the header is the very first line of both files.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Actually you are right, it works, but the header is in the middle not on the top... is there a way to keep it on the top? thanks

• 0 views • link
written by User000 User000

use sed to change the header into someting that should be at the top after sort. Somtehing like 's/^chromosome/00000chromosome/'

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

You can do it with R, with the subset function works pretty intuitively.

• 0 views • link
written by Benn Benn

In R:

file1<-read.table("file1.txt", sep="\t", header=T)
file2<-read.table("file2.txt", sep="\t", header=T)
Selection<-file1[file1$"Gene name" %in% file2$"Gene Name",]

You don't even have to use subset function

• 0 views • link
updated by Ram Ram • written by Benn Benn

Dear Nota,

Thank you for your answer. These command in R just gave me the headers:

Gene.Name   Gene.Id     description GO
• 0 views • link
updated by Ram Ram • written by ahmad.iut ahmad.iut

OK, R substitutes the spaces in the header to dots.

So you can use:

Selection<-file1[file1$Gene.Name %in% file2$Gene.Name,]
• 0 views • link
updated by Ram Ram • written by Benn Benn

Thank you so much Nota, It worked well. the problem was the spaces in headers (like Gene Name).

• 0 views • link
written by ahmad.iut ahmad.iut

Log in to answer this question.

More posts like this