How to download gene annotation from NCBI?

[jgarbe]

The NCBI Map Viewer has the latest pig genome build and shows the
locations of all the genes. I would like to download this gene
annotation so I can load it into my own GBrowse genome browser.
So I need the NCBI gene annotation for the latest pig genome build in
gff3 format, and the way to do it seems to be to download an asn.1
file from NCBI, convert it to genbank format, and then use the bioperl
script bp_genbank2gff3.pl to convert from genbank to gff3.

I downloaded the gene annotation for the pig genome from the NCBI ftp site at
ftp://ftp.ncbi.nlm.nih.gov/gene/DATA..._scrofa.ags.gz

I downloaded the asn2gb conversion program from
ftp://ftp.ncbi.nlm.nih.gov/asn1-conv...latform/linux/

I run ./linux.asn2gb -i Sus_scrofa.ags -b T
and get the error "Asn io_failure for input file 'Sus_scrofa.ags'"
I've tried all the options for the -a and -t flags without luck.

I'm able to convert the Sus_scrofa.ags file to xml format using the
gene2xml program, but I don't know of any tool that can convert from
XML to gff3.
I downloaded a genbank format file of pig genes from
ftp://ftp.ncbi.nlm.nih.gov/genomes/S...RNA/rna.gbk.gz but the
file doesn't give chromosome coordinates for the genes, so I can't
make a gff3 file out of it.

Any pointers on how to use the asn tools properly, or how to get NCBI
annotation in gff format in general, would be much appreciated.

Thanks

-John

[darked89]

I managed to run:

./gene2xml.linux -i Sus_scrofa.ags -b T -c T

This prints XML output. Strangely Sus_scrofa.ags had to be gzipped and named Sus_scrofa.ags.gz.

XML to gff convertion should be fairly easy, but I do not know a tool yet. You may check:

http://brendelgroup.org/mespar1/gthxml/gthxmlToGFF.py

[jjw14]

John,

I'm running into the same problem that you had. The NCBI Sus scrofa genome FTP site provides .asn, .fa, .gbk, .gbs, and .mfa files for each chromosome (last updated 10-12-2011).

Were you able to convert the .asn data to .gff3 or .gtf format for annotation? I'd be interested to hear the best method you found for generating the annotation file that corresponds to the most recent S. scrofa genome.

Thanks in advance,
jjw

[jgarbe]

I was not able to figure out how to convert any of the NCBI annotation data into a usable form. I sent an email to NCBI but didn't get a useful reply from them. Thankfully another group has generated a good gene build for Sscr10.2. As described here: http://animalgenome.org/pig/newsletter/No.110.html, you can download annotation at this site: http://gbi.agrsci.dk/pig/sscrofa10_2_annotation/
Alternatively, Ensembl is running 10.2 through their pipeline and should have a gene build available in two or three months. If you can wait that long that would be another good alternative to NCBI's annotation.

[jjw14]

Thanks for the quick reply, John.

No doubt, you've saved me a lot of frustration. I appreciate it.

jjw

[maubp]

Any pointers on how to use the asn tools properly, or how to get NCBI annotation in gff format in general, would be much appreciated.

The NCBI are currently revising all their GFF3 output (it hadn't been compliant with the standards), so this should be much easier now/soon.

Try ftp://ftp.ncbi.nlm.nih.gov/genomes/Sus_scrofa/GFF/ for the NCBI RefSeq annotation of pig Sscrofa10.2

[jjw14]

Thanks Peter,

I'll take a look.

jjw14

[noha osman]

Hi ,all
I need buffalo gff or.ggf3 file from NCBI but I donot know how can get it .
Could anyone help me to know the answer
Thanks

[maubp]

Bubalus bubalis? ftp://ftp.ncbi.nlm.nih.gov/genomes/Bubalus_bubalis/GFF/

[noha osman]

Thanks Peter