-2
$\begingroup$

Are there any databases or tools that contain a list of all RefSeq protein IDs for the interest organism/microorganism, e.g. get a specific microorganism name and return all RefSeq Protein IDs that have been identified?

$\endgroup$
0
$\begingroup$

How about using RefSeq itself? Say, you are interested in all E. coli proteins, you can use the query "Escherichia coli"[Organism] AND RefSeq[Filter] here.

If you know the Taxonomy identifier for the specific organism of interest to you, it can also be used. For example, if you are interested in only "Escherichia coli str. K-12 substr. DH10B", you can use the taxonomy identifier 316385 (obtained from here) in the query txid316385[Organism] AND RefSeq[Filter].

To extract a list of accessions and save them to a file on disk, use the "Send To" menu at the top right corner, choose "File" as destination followed by "Accession List" as the format and download the file. Alternatively, you can use E-utilities API or the EntrezDirect tool on the Unix command line to download these data programmatically.

To address the question about E. coli proteins shared between different species... There's no need to parse FTP files to get all the RefSeq proteins for a bacteria species. While it's true that some RefSeq WP proteins have taxids above species level, like Enterobacterales WP_003861789.1 where the identical protein is found in genomes from multiple species of the order Enterobacterales, those proteins are still retrieved by a species-level organism search. NCBI implemented some extra handling for the RefSeq WP proteins to treat the MULTISPECIES WPs like they belong to all of the species with genomes on which they're annotated. So of the 56947 WPs assigned as "Bacteria", 7401 are found by "Escherichia coli"[Organism], and the species-level search will give you the full set of proteins.

You could also use the AutonomousProtein2Genomic.gz file from the RefSeq FTP path which has a mapping of how all WPs map to specific genomic sequences, but it's a really big file and it would take a long time to grep for "Escherichia coli".

$\endgroup$
3
  • $\begingroup$ Thanks warmly for the response. Actually, I was looking for a database in which I can extract all RefSeq protein IDs as a .csv or .xls file. $\endgroup$ Jun 22 at 17:17
  • $\begingroup$ I would say,suggested query is not working as expected for refseq.Some escherichia coli refseq proteins have “lower” taxonomy rank ,for example “”Escherichia- proteins which shared between different Escherichia species.Technically ,I know only one way for now extract ALL Escherichia Coli refseq proteins $\endgroup$
    – user6690
    Jun 27 at 9:04
  • $\begingroup$ retrieve all corresponding refseq assemblies from ftp and parse them $\endgroup$
    – user6690
    Jun 27 at 9:11

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.