Is there a way to filter out files by a taxon and a tissue in SRAdb?
2
0
Entering edit mode
6.9 years ago
Pavel Senin ★ 1.9k

Hi guys:

I wonder how to get links for FASTQ files using Bioconductor SRAdb package for specific taxons and tissues? Currently I am trying to get through by using JOINs and fulltext searches with LIKE, but it is so slow to figure it out..., maybe someone already has solved that?

Thanks!

SRAdb • 1.5k views
ADD COMMENT
1
Entering edit mode
6.9 years ago
David W 4.8k

I don't know about SRAdb or what metadata that database contains, but you can use entrez to search SRA with the normal limits (including organism)

library(rentrez)

entrez_search(db="sra", term="Tetrahymena[ORGN]")
# Entrez search result with 113 hits (object contains 20 IDs and no cookie)
 

It looks like tissue isn't a specified term for tissue

 

> entrez_db_searchable("sra")

 

Searchable fields for database 'sra'
  ALL      All terms from all searchable fields 
  UID      Unique number assigned to publication 
  FILT      Limits the records 
  ACCN      Accession number of sequence 
  TITL      Words in definition line 
  PROP      Classification by source qualifiers and molecule type 
  WORD      Free text associated with record 
  ORGN      Scientific and common names of organism, and all higher levels of taxonomy 
  AUTH      Author(s) of publication 
  PDAT      Date sequence added to GenBank 
  MDAT      Date of last update 
  GPRJ      BioProject 
  BSPL      BioSample 

 

But you could search for biosamples that match you query, then find links to SRA?

 

entrez_search(db="biosample", term="Mouse[ORGN] AND Tissue Liver[Attribute]")

Entrez search result with 2934 hits (object contains 20 IDs and no cookie)

 

 

ADD COMMENT
0
Entering edit mode

Hi David: thank you for this excellent example! Right now I am trying to make use of the local SRAdb without much online querying, if not working, that shall be a next step.

ADD REPLY
0
Entering edit mode
6.8 years ago
Lynxoid ▴ 220

Could SHARQ be something like what you're looking for?

ADD COMMENT

Login before adding your answer.

Traffic: 2589 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6