How does NCBI populate data fro gene entries? I want to get all refseq mRNA..

[beeman]

Hi All,

I'm just wondering if anyone can shed light on how to obtain the latest annotations of a given organism from NCBI, and more specifically how to get all of the current transcript variants that are listed as refseq's..

I'm working on honeybees and I've grabbed the latest gff flatfile annotations from ftp://ftp.ncbi.nih.gov/genomes/Apis_mellifera/GFF/ but they don't contain all of the current refseq transcripts..

For example, the gene cort (http://www.ncbi.nlm.nih.gov/gene/726912) has three transcript variants listed as refseq entries; XM_006557348.1, XM_001122629.3 and XM_006557349.1, but in the gff annotations the only transcript ID is XM_001122629.2...

Is there any way to build a current set of annotations from the data NCBI uses to populate transcripts for gene records?


Thanks

[GenoMax]

You can do this in a couple of different ways. One would be to get the invertebrate RefSeq data files (ftp://ftp.ncbi.nlm.nih.gov/refseq/release/invertebrate/) and parse out Apis mellifera entries.

A simpler way would be to do a search: http://www.ncbi.nlm.nih.gov/nuccore/...D+srcdb_refseq. Change the "Display settings" to indicate the format you want (GenBank, Fasta etc) and then "Send to" a File. I currently see ~40,500 entries for the search above.

Keep in mind that the above list will likely include all of these RefSeq types (http://www.ncbi.nlm.nih.gov/books/NB...ort=objectonly).

[Ohad]

I downloaded ref_Amel_4.5_scaffolds.gff3

And I see XM_006557348.1, XM_001122629.3 and XM_006557349.1 inside...

[beeman]

Thanks for your replies..

Definitely operator error...

I need to start appending version dates as I've been looking at the wrong flatfile ....