ncbi-genome-download批量下载基因组数据

本文主要是介绍ncbi-genome-download批量下载基因组数据，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

1. ncbi-genome-download 的下载和安装

ncbi-genome-download 是一个可以直接从NCBI上批量下载序列的软件，支持下载多种格式

利用 conda 对其直接安装 参考

#创建环境
conda create -n ncbi_genome_download
#激活环境
conda activate ncbi_genome_download
# 安装
conda install -c bioconda ncbi-genome-download

2.常用的参数

• -s：选择数据库（genbank，refseq），默认是refseq数据库
• -F：需要下载基因组的格式，可以多种格式同时下载，用逗号隔开，默认是genbank格式
• -l：序列组装程度，可以多种格式同时下载，用逗号隔开
• -g：需要下载序列的属，后面要指定类群，比如bacteria
• S：下载的具体的菌种名称，用逗号隔开，也可以写入一个文件中，一行一个菌种名称
• -o：输出的文件名称
• -r：失败时重新连接的次数，默认是0次
• --flat-output：将下载的文件输入到一个目录中，不创建新的子文件（即下载的数据在指定的文件夹中，每个 Taxonomy ID 一个压缩文件）

3. 批量下载基因组数据

3.1 根据属名下载

将需要下载的属名放置至一个txt文档（换行），利用参数--genera  pant_download.txt  plant

3.2 根据物种的 ID下载

将需要下载的物种的分类 ID 放置至一个txt文档（换行），利用参数--taxids my_taxids.txt 

再加上参数 --assembly-levels 指定下载的基因组的不同类型（包括contig,scaffold,chromosome,all,compete）

不同物种的 taxonomy id查询地址 taxonomy id query

例如 Oryza sativa 的 taxonomy id为4530

3.3 根据物种拉丁名下载

当你有一系列菌种需要下载时，你可以将这一系列菌种名保存到一个txt文件里，每个菌种名为一行，文件名为genera.txt

希望下载这些菌种基因组中的cds序列，并将下载的每个文件放在MyGenera文件夹中，在MyGenera目录下进入终端，运行：

ncbi-genome-download --genera genera.txt bacteria --flat-output --formats cds-fasta

然后，每个物种均会自动下载好指定的基因组类型序列

这样genera文本中的所有菌种的基因组cds序列就一条代码下载完成了。

注意：genera.txt文本中有10个菌种名，而下载了13个文件，说明有的菌种名下面有来自不同上传者提供的基因组信息（即一对多）

如下图，我要下载GCAlist.txt文本中的基因组序列的fasta文件，并保存在Assembly文件夹中，在Assembly目录下进入终端，运行：

ncbi-genome-download --assembly-accessions GCAlist.txt bacteria --section genbank --for

4.帮助查询

查询地址 帮助文档

• 查看版本

ncbi-genome-download -V

• 查看帮助

ncbi-genome-download -h

• 语法格式

ncbi-genome-download [optional arguments] groups

optional arguments为可选参数，详细介绍见下文

groups为物种选择，可选['all', 'archaea', 'bacteria', 'fungi', 'invertebrate', 'metagenomes', 'plant', 'protozoa', 'vertebrate_mammalian', 'vertebrate_other', 'viral']，可选项即为NCBI的FTP下载目录Index of/genomes/refseq和Index of/genomes/genbank下的内容

• 可选参数

--section

指定下载的数据库，可选['refseq', 'genbank']，默认refseq

--formats

指定下载的文件格式，可选['genbank', 'fasta', 'rm', 'features', 'gff', 'protein-fasta', 'genpept', 'wgs', 'cds-fasta', 'rna-fna', 'rna-fasta', 'assembly-report', 'assembly-stats', 'all']，默认genbank

--assembly-levels

指定下载的基因组组装水平，可选['all', 'complete', 'chromosome', 'scaffold', 'contig']，默认all

--genera

根据菌种名下载，后面可接想要下载的菌种名，如--genera 'Rhizobium alamii'

--taxids

根据NCBI taxonomy ID下载，后面可接想要下载的菌种的taxonomy ID，如--taxids '492774'

（还以Rhizobium alamii举例，通过NCBI Taxonomy Browser可以查询到该菌种的txid为492774）

--assembly-accessions

根据assembly accession下载，后面可接想要下载的菌种的assembly accession，如--assembly-accessions ‘GCF_000799895.1’

⚠️注意：因为默认下载的数据库是refseq，所以选择RefSeq assembly accession下载时无需加--section参数即可正常下载，如果要根据GenBank assembly accession下载，请再加上--section genbank。

--output-folder

指定下载目录，后面可接你想要存放的下载目录，如--output-folder ～/Downloads（下载到当前用户的下载文件夹中）

--flat-output

直接将下载的文件放入指定文件夹中，不创建子文件夹

详情参考 ncbi-genome-download工具

5. 核查下载情况

由于利用ncbi-genome-download下载物种的基因组数据时存在未成功下载（后续需自己手动下载）的情况，所有需要将当前目录下的a.genomic.fna.gz文件进行汇总

5.1 利用grep命令查找

随后将sequence_name.txt文件导入excel表中，与自己需要下载的物种进行vlookup函数匹配，找出未成功下载的基因组序列的物种名

grep ".*genomic.fna.gz” ./present dictionary > sequence_name.txt

5.2 将.txt文件转换成.bat文件

另一种方法，比较简单实用：

这篇关于ncbi-genome-download批量下载基因组数据的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

---