This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Obtaining GSM numbers from SRA ids

I have a lot of SRS values from the SRA values/ids for relevant entries in the SRA, and I would like to find the corresponding GSM ids that match to these. I can find them via the SRA webpage manually, but would really like a way to use curl or python (not R) to get these values--is that possible?

rna-seq sra
• 3,434 views • link
updated by Biostar Biostar • written by Dirk Dirk

2 answers

Can you give an example to test?

Take a look at this answer: A: Get a complete GSM-to-SRX/SRR table

• 1 views • link
written by GenoMax GenoMax

Awesome, this looks like the best answer. It's a bit unfortunate you have to load a ~4GB file, but its still completely viable!

• 1 views • link
written by Dirk Dirk

If you don't want to download the entire SRA_Accessions.tab file, you can use Entrez Direct for this: 

esearch -db sra -q 'SRS213308' | elink -db sra -target gds -name sra_gds | esummary | xtract -pattern DocumentSummary -subset DocumentSummary -if Accession -starts-with 'GSM' -first Accession
GSE30198
GSM747492
• 1 views • link
written by vkkodali_ncbi vkkodali_ncbi

Curious why this is picking a GSE accession when you are asking for GSM?

• 1 views • link
written by GenoMax GenoMax

The SRS accession is linked to both the GSM sample accession and the GSE series accessions in GEO Datasets. I could not come up with a good xtract command that will filter out the GSE hit. But it can be easily dropped using unix grep.

• 1 views • link
written by vkkodali_ncbi vkkodali_ncbi

Log in to answer this question.

More posts like this