This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to extract Locus tag from GeneIDs in NCBI for soybean database.

Hi all, I have been trying to get information on converting NCBI GeneID to Glyma ID for soybean gene annotation. However it seems like such file does not exist. For instance, GeneID:100790502 search in NCBI shows Locus tag for this id as GLYMA_09G197400 which is what I want to extract but for all 56K soybean genes. Is there a way to query all 56K NCBI format GeneIDs and extract the Locus tag value? That way, I can use Glyma IDs and annotate them using phytozome annotation file. Phytozome already has their genes in Glyma format and I am not sure why NCBI does not have this information in their fasta or gff/gtf files.

Your insight is much appreciated.

Best,

rnaseq soybean tag locus id glyma

2 answers

The file gene_info.gz on the Gene FTP site has this information. Since you are interested in soybean only, you can download the All_Plants.gene_info.gz file from here. On a Unix command line, you can extract these as follows:

zcat All_Plants.gene_info.gz | awk 'BEGIN{FS="\t";OFS="\t"}(($1~/^#/||$1==3847) && $4~/GLYMA/)'

This is great information and saved me a lot of headache. Thanks a lot. Appreciate it. Best.

You can use EntrezDirect. First column as Entrez gene ID and second Locus tag :

$ esearch -db gene -query 100790502  | esummary | xtract -pattern DocumentSummary -element Id,OtherAliases
100790502   GLYMA_09G197400

This may get you most of them. Showing only 10 here (remove | head -10 to get them all) :

$ esearch -db gene -query GLYMA | esummary | xtract -pattern DocumentSummary -element Id,OtherAliases | head -10
547923  GLYMA_13G347600, L-1, Lx1
548076  GLYMA_13G288100
547831  GLYMA_08G341500, KTi, Ti-a, Ti-b, Tia, Tic, Tie
100788438   GLYMA_03G181700, GmPAL1.2, PAL1
547900  GLYMA_03G163500, A2B1a, glycinin
547641  GLYMA_18G023500, RLK-RHG1, rhg1-like, rhg1g, rhg1s
547931  GLYMA_06G301500, BMY1, Gm-BamyDam, Gm-BamyKza
100787872   GLYMA_02G309300, GmPAL3.1
100527427   GLYMA_10G199100, GmLb, N-2, Nodulin-2
547869  GLYMA_15G026300, L-3, LOX1.3, Lx3

This looks like what I needed. I will try this. Thanks a lot @genomax. I really appreciate it. Best.

Log in to answer this question.