How to link the assembly accession with the chromosome accession for prokaryotic representative genomes
1
0
Entering edit mode
2.9 years ago
wangdp123 ▴ 250

Dear colleague,

I am working on the analysis of prokaryotic genomes from NCBI genome database.

  1. Downloaded a file called prok_representative_genomes.txt from the following file ftp://ftp.ncbi.nlm.nih.gov/genomes/GENOME_REPORTS/prok_representative_genomes.txt

After opening the file, we could see one column named "Chromosome RefSeq". (e.g., NZ_AQXM00000000)

  1. Download all protein sequences for all bacterial genomes from https://www.ncbi.nlm.nih.gov/assembly.

Each file has a name like "GCF_000834735.1_ASM83473v1_protein.faa.gz".

It is odd that the two datasets use different accession number system. In this case, how to identify if the genome of the proteins is annotated in the file "prok_representative_genomes.txt"?

My aim is to retrieve all the protein sequences for the genomes listed in the file "prok_representative_genomes.txt".

Thanks a lot,

Kind regards

Tom

NCBI prokaryotic genomes • 678 views
ADD COMMENT
0
Entering edit mode

Thanks,

Does anybody know why there are some 0 values in the "Genomes" column of the file prok_representative_genomes.txt? For example, "Acetobacter aceti".

ADD REPLY
0
Entering edit mode

The README says it's manually updated. I tend to believe you'll see a more complete picture when using the prokaryotes.txt (computationally updated...)

ADD REPLY
1
Entering edit mode
2.9 years ago
vkkodali ★ 2.7k

In the same FTP path, there is another file prokaryotes.txt that can be helpful. Join the two files on the organism name (column 3 of prok_representative_genomes.txt and column 1 of prokaryotes.txt). Then, extract the path to the assembly folder on the FTP site (column 21 of prokaryotes.txt file) and download the *.protein.faa.gz files.

Accessions with the NZ_ prefix are RefSeq Chromosome accessions whereas the ones with GCF_ prefix are assembly accessions. An assembly includes both chromosome and plasmid(s).

ADD COMMENT
0
Entering edit mode

As vkkodali wrote, use the prokayotes.txt. However, there's no need to join, one can filter on the Reference column. The NCBI discriminates between manually selected reference (exist for common lab organisms) and computationally selected representative genomes (exists for each species, almost).

ADD REPLY

Login before adding your answer.

Traffic: 1958 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