GTDB基因组分类数据库及GTDB-Tk工具

Introduction

GTDB（Genome Taxonomy
Database）是一个专注于基因组分类学的数据库，旨在提供高质量的细菌、古细菌、叶绿体和线粒体的分类信息。通过利用全基因组信息，GTDB提供了多级别的分类，包括域、界、纲、目、科、属等级别。其体系发育分类体系使得研究职员能够深入相识生物的演化和功能。GTDB重要以细菌中普遍存在的120个单拷贝蛋白质（bac120）和古菌中53个（arc53，从R07-RS207开始）标志蛋白为基础，在对多分组种别消歧后，根据相对演化散度标准化和分级，得到基因组分类结果。
GTDB定期更新，与其他生物信息数据库整合，为基因组学研究和微生物生态学等范畴的研究提供了全面且可靠的分类信息。GTDB数据库已经于2018/2020连续发表两篇Nature
Biotechnology，2021年又发表在Nature Microbiology和Nucleic Acids
Research上，是现在最权势巨子的微生物分类数据库之一。GTDB配套软件在2019/2022发表了两篇Bioinformatcs，时提供了物种分类、进化树构建的全套工具。
官网：https://gtdb.ecogenomic.org/

GTDB数据库

可以看到上面的官网主页，现在最新版的GTDB数据库包含了584,382个细菌基因组和12,477个古菌基因组（2024.04.24
日，Release 09-RS220），在属、种分辨率水平上描述了190多个门。
Browsers

GTDB支持以Alphabetical
table和Tree两种形式来探索或查询整个GTDB数据库，包括分类信息和基因组信息。
此外，可以在Taxon History输入 GTDB 或 NCBI
分类单元，以可视化分配给该分类单元的基因组在版本之间怎样变化，好比p__Proteobacteria从Release
214开始改名为p__Pseudomonadota，在不停壮大的同时也有一部分genome被分配给了其他门。
Tools

GTDB提供了一些配套的工具：

• GTDB-Tk：下一章节详细先容
  
• Third-party tools：AnnoTree，IDTAXA，SingleM，sourmash
  
• Third-party GTDB
  databases：HUMANn2，Kraken2，Bracken2，Centrifuge，Kraken，DADA2等
  
• Third-party
  scripts：tax_from_gtdb.py，genome_updater，gtdb_to_taxdump
  
• FastANI calculator：可用于原核生物 NCBI
  基因组的ANI比力，一次最多可以举行 1,000 次成对比力，比方下图。
  

Downloads

该页面提供了GTDB数据库和配套软件的下载链接。
Statistics

统计页面提供了各个版本的GTDB数据库的各类基础信息：

• Taxon overview：GTDB R220 包含 596,859 个基因组，组织成 113,104
  个物种簇。
  
• Species overview：GTDB R220 由 584,382 个细菌和 12,477
  个古细菌基因组构成，分为 107,235 个细菌和 5,869 个古细菌物种簇。
  
• Genome categories：GTDB
  分类群由分离基因组（isolate）、宏基因组组装基因组 (MAG)
  和单扩增基因组 (SAG) 构成。
  
• GTDB species representatives：每个 GTDB 物种簇都由一个基因组表示。
  
• Quality of GTDB representative genomes：被选为 GTDB
  物种代表的基因组的质量。利用 CheckM 估计基因组完备性和污染，并根据
  MIMAG
  基因组标准举行着色。一般来说，代表性基因组仅限于满足completeness - 5*contamination >50，存在一些例外情况，如下面所示。
  
  
  
  
• Taxa with the largest number of species：给出了每个分类等级包含最多
  GTDB
  物种簇的分类单元，门水平是常见的Pseudomonadota，Bacillota，Bacteroidota，Actinomycetota等。
  
• Taxa with the largest number of sequenced genomes，与6类似。
  
• Relative evolutionary
  divergence（RED）：提供了从门到属的每个分类等级的类群的相对进化分歧（RED）。
  RED 值提供了相对时间的操纵近似值，此中现有分类群存在于当前
  (RED=1)，末了一个共同先人出现在已往的固定时间
  (RED=0)，内部节点根据以下值在这些值之间线性插值谱系特定的进化速率。
  
• Comparison of GTDB and NCBI taxa：GTDB
  物种代表性基因组和所有具有分配的 NCBI 分类法的 GTDB 基因组的 GTDB 和
  NCBI
  分类学分配的比力。对于每个分类等级，如果一个分类单元在两个分类标准中名称相同，则该分类单元被分类为unchange；如果
  GTDB 分类标准提供了 NCBI
  分类标准中缺少的名称信息，则该分类单元被分类为passive change；如果两个分类标准之间的名称差别，则该分类单元被分类为actively changed。可以看到，在genome水平上GTDB和NCBI还是有很多名称不对应的，可以思量利用gtdb_to_taxdump举行转化。
  
  
  
  
• Genomic statistics：GTDB 物种代表性基因组和 GTDB
  中所有基因组的关键基因组统计数据。
  
  
  
• Nomenclatural types per
  rank：显示了每个分类等级的占位符与拉丁化名称的细分。
  

Forum

GTDB还专门搭建了一个用户社区和讨论区：

• Announcement：用于公布最新 GTDB
  新闻、即将举行的会议和以分类学为重点的手稿的开放论坛。
  
• Taxonomy
  (Bacteria/Archaea)：这个分类专注于讨论与GTDB分类法相关的内容。用户在这里讨论分类法的具体细节、与其他数据库（如NCBI）的比力、分类标准的更新等。
  
• GTDB-Tk：专注于GTDB-Tk工具包的利用和讨论。用户可以在这里提出关于GTDB-Tk的利用问题、分享利用履历、讨论工具包的功能和性能等。
  
• GTDB：与基因组分类数据库 (GTDB) 相关的一般问题。
  
• Website feedback：讨论该网站、其组织、其运作方式以及怎样改进它。
  

Help

资助页面也提供了一些非常有效的信息。

• About：GTDB简要先容，管理团队以及发表论文
  
• FAQ：列出了一些常见问题，比方：
  
  
  
  • 为什么门名称的后缀更改为-ota？Whitman et al. (2018)
    建议像其他等级一样规范门等级的后缀。
    
  • 为什么一些family和更高级别的名字以字母后缀结尾？属等级以上的分类单元名称附加字母后缀表示属于以下种别的类群：
    
    
    • GTDB 参考树中不是单系类群，但存在其他证据表明它们是单系类群；
      
    • 在版本之间位置不稳定的组。
      
    
    
  
  
• Methods：对GTDB数据库构建过程中的方法举行具体描述，比方：
  
  
  
  • Gene
    identification，利用Prodigal判断基因，HMMER判断标志基因（来自Pfam和TIGRFAM）并序列比对。
    
  • Multiple sequence alignment，细菌和古菌多重序列比对 (MSA) 分别由
    120 个 (bac120) 或 53 个 (arc53) 体系发育信息标志串联而成。
    
  • Tree inference，细菌参考树是在 WAG 模子下利用
    FastTree。古菌参考树是根据 PMSF 模子下的 IQ-Tree，并利用
    FastTree推断初始树以规复单系谱系。
    
  • Identifying 16S rRNA sequences，利用 nhmmer以及来自 RFAM
    数据库的 16S rRNA 模子（RF00177）。
    
  • Average nucleotide identity，skani（从版本 220
    开始），之前是FastANI。
    
  • RED
    normalization：相对进化分歧用于资助刻画物种等级之上的分类群，RED的概念是相划一级的类群在已往应该同时出现。RED
    提供了相对时间的操纵近似值，此中现有分类群存在于当前 (RED =
    1)，末了一个共同先人出现在已往的固定时间 (RED =
    0)，内部节点根据谱系在这些值之间线性插值-特定的进化速率。
    
  
  
• API：功能较多，但文档不是很详细，未来也许会用上
  

GTDB-Tk

GTDB-Tk 是一个软件工具包，用于根据基因组数据库分类法 (GTDB)
对细菌和古细菌基因组举行客观分类。
它旨在与最新希望相结合，答应直接从情况样本中获得数百或数千个宏基因组组装基因组
(MAG)，还可以应用于分离和单细胞基因组。
所在：https://github.com/Ecogenomics/GTDBTk
官方文档：https://ecogenomics.github.io/GTDBTk/
Install

起主要安装软件和下载对应的Reference
data，最新的R220数据库必要利用v2.4.0的GTDB-Tk（2024.04.24）。如果要和之前项目划一，可以利用R214版本数据库对应v2.3.2的GTDB-Tk。
1.conda

conda安装比力简单，但最好新建一个情况，可能会比力慢（mamba代替conda会快一点）：

1. # using conda
2. conda create -n gtdbtk-2.3.2 -c conda-forge -c bioconda gtdbtk=2.3.2
3. # 下载 GTDB-Tk 参考数据并为其命名
4. download-db.sh
5. # 激活GTDB-Tk conda环境
6. conda 激活 gtdbtk-2.3.2
7. # 将环境变量设置为包含GTDB-Tk参考数据的目录
8. conda env config vars set GTDBTK_DATA_PATH="/path/to/release/package/";

复制代码

2.pip

pip可以快速安装gtdbtk包，但是它的第三方依赖项必要手动安装一下
（不是很复杂，我用的pip安装，因为不想要太多的conda情况，太冗余了）：

1. python -m pip install gtdbtk==2.3.2

复制代码

第三方依赖项：

• Prodigal（原核基因识别和翻译起始位点判断）
  
• HMMER（加速族HMM 搜索）
  
• pplacer（线性时间最大似然和贝叶斯体系发育将序列放置到固定参考树上）
  
• skani (v2.4.0开始依赖，通过 skani
  的希奇链接举行快速、稳健的宏基因组序列比力，比fastANI还快，下次有机会讲一下这个)
  
• 
  
• FastTree（用于大对齐的近似最大似然树）
  
• Mash（快速估计基因组和宏基因组距离）
  

这些都是很常用的软件，没有的话用conda安装一下。
GTDB-Tk 必要将 GTDBTK_DATA_PATH 情况变量设置为包含未存档的 GTDB-Tk
参考数据的目录。

1. echo 'export GTDBTK_DATA_PATH=/path/to/release/package/' >> ~/.bashrc
2. source ~/.bashrc

复制代码

3.docker

如果用过docker的话，这个也非常简单：

1. #下载解压的 GTDB-Tk 参考数据到您选择的目录，例如：/host/release_data/
2. mv /path/to/release/package/ /host/release_data/
3. #为 GTDB-Tk 创建本地 I/O 目录
4. mkdir /host/gtdbtk_io/
5. # 运行容器
6. docker run -v /host/gtdbtk_io:/data -v /host/release_data:/refdata ecogenomic/gtdbtk classify_wf --genome_dir /data/genomes --out_dir /data/output

复制代码

Reference data

GTDB-Tk 必要约 110G
的外部数据，必要下载息争档。直接wget很慢的话，也可以先下载到本地，再通过FTP软件或scp上传到服务器：
数据所在：https://data.gtdb.ecogenomic.org/releases/latest/auxillary_files/gtdbtk_package/
对于完备包(full package)：

1. wget https://data.ace.uq.edu.au/public/gtdb/data/releases/latest/auxillary_files/gtdbtk_package/full_package/gtdbtk_data.tar.gz
2. tar xvzf gtdbtk_data.tar.gz

复制代码

对于拆分包(split package)：

• 确保 GTDB-Tk 发布数据的所有部分都位于同一目录中。
  
• 打开终端或下令提示符。
  
• 导航到包含 GTDB-Tk 发布数据部分的目录。
  
• 利用以下下令将所有部分毗连到一个存档中：
  cat gtdbtk_r220_data.tar.gz.part_* > gtdbtk_r220_data.tar.gz
  
• 下令执行完成后，将在同一目录中拥有一个名为“gtdbtk_r220_data.tar.gz”的存档文件。
  

Check

检查数据库和软件，结果显示OK和Done表示正常（检查fastani这一步要好久，文件非常多。以是如果三方软件都OK的话也可以直接运行下一个示例测试）:

1. gtdbtk check_install

复制代码

1. [2024-05-20 14:10:37] INFO: GTDB-Tk v2.3.2[2024-05-20 14:10:37] INFO: gtdbtk check_install
2. [2024-05-20 14:10:37] INFO: Using GTDB-Tk reference data version r214: /data/home/jianglab/share/GTDB_r214/release214/[2024-05-20 14:10:37] INFO: Running install verification[2024-05-20 14:10:37] INFO: Checking that all third-party software are on the system path:[2024-05-20 14:10:37] INFO:          |-- FastTree         OK[2024-05-20 14:10:37] INFO:          |-- FastTreeMP       OK[2024-05-20 14:10:37] INFO:          |-- fastANI          OK[2024-05-20 14:10:37] INFO:          |-- guppy            OK[2024-05-20 14:10:37] INFO:          |-- hmmalign         OK[2024-05-20 14:10:37] INFO:          |-- hmmsearch        OK[2024-05-20 14:10:37] INFO:          |-- mash             OK[2024-05-20 14:10:37] INFO:          |-- pplacer          OK[2024-05-20 14:10:37] INFO:          |-- prodigal         OK[2024-05-20 14:10:37] INFO: Checking integrity of reference package: /data/home/jianglab/share/GTDB_r214/release214/[2024-05-20 14:10:38] INFO:          |-- pplacer          OK[2024-05-20 14:10:38] INFO:          |-- masks            OK[2024-05-20 14:10:39] INFO:          |-- markers          OK[2024-05-20 14:10:39] INFO:          |-- radii            OK[2024-05-20 14:10:41] INFO:          |-- msa              OK[2024-05-20 14:10:41] INFO:          |-- metadata         OK[2024-05-20 14:10:41] INFO:          |-- taxonomy         OK[2024-05-20 15:13:45] INFO:          |-- fastani          OK[2024-05-20 15:13:45] INFO:          |-- mrca_red         OK[2024-05-20 15:13:45] INFO: Done.

复制代码

注意GTDB-Tk的硬件要求，关键是Memory一定要给足够，否则会程序会中途终止：
DomainMemoryStorageTimeArchaea~60 GB~106 GB~90 minutes / 1,000 genomes @ 64 CPUsBacteria~90GB (545 GB when using –full_tree)~106 GB~90 minutes / 1,000 genomes @ 64 CPUs 测试流程，无报错确定程序的软件和数据库正常（此示例是3个古菌genome，建议分配60GB内存来运行，别的gtdbtk_test必须是不存在的或空文件夹）:

1. gtdbtk test --out_dir gtdbtk_test

复制代码

1. [2024-05-20 14:34:51] INFO: GTDB-Tk v2.3.2[2024-05-20 14:34:51] INFO: gtdbtk test --out_dir gtdbtk_test
2. [2024-05-20 14:34:51] INFO: Using GTDB-Tk reference data version r214: /data/home/jianglab/share/GTDB_r214/release214/[2024-05-20 14:34:51] INFO: Command: gtdbtk classify_wf --skip_ani_screen --genome_dir gtdbtk_test/genomes --out_dir gtdbtk_test/output --cpus 1 -f    <TEST OUTPUT>     [2024-05-20 14:34:52] INFO: gtdbtk classify_wf --skip_ani_screen --genome_dir gtdbtk_test/genomes --out_dir gtdbtk_test/        <TEST OUTPUT>     [2024-05-20 14:34:52] INFO: Using GTDB-Tk reference data version r214: /data/home/jianglab/share/GTDB_r214/release214/          <TEST OUTPUT>     [2024-05-20 14:42:06] INFO: Note that Tk classification mode is insufficient for publication of new taxonomic designatio        <TEST OUTPUT>     [2024-05-20 14:42:06] INFO: Done.     [2024-05-20 14:42:09] INFO: Test has successfully finished.

复制代码

我分配了4核共60GB内存运行gtdbtk test，约439s完成，消耗了~45GB内存。
利用方法

先看看help文档：

1. $ gtdbtk -h
2. ## GTDB-Tk v2.4.0 帮助文档
4.   Workflows:
5.     classify_wf -> 通过在 GTDB 参考树中的位置对基因组进行分类
6.                      (ani_screening -> identify -> align -> classify)
7.     de_novo_wf  -> 推断新的树并使用 GTDB 分类进行修饰
8.                      (identify -> align -> infer -> root -> decorate)
10.   Methods:
11.     identify -> 识别基因组中的标记基因
12.     align    -> 创建多序列比对
13.     classify -> 确定基因组的分类
14.     infer    -> 从多序列比对中推断树
15.     root     -> 使用外群对树进行根定
16.     decorate -> 用 GTDB 分类修饰树
18.   Tools:
19.     infer_ranks        -> 利用RED建立内部节点的分类等级
20.     ani_rep            -> 计算ANI到GTDB的代表性基因组
21.     trim_msa           -> 基于掩码修剪未修剪的MSA文件
22.     export_msa         -> 导出未修剪的古生菌或细菌MSA文件
23.     remove_labels      -> 从Newick树中删除标签(引导值，节点标签)
24.     convert_to_itol    -> 将GTDB-Tk Newick树转换为iTOL树
25.     convert_to_species -> 将GTDB基因组id转换为GTDB物种名称
28.   Testing:
29.     test          -> 用3个古细菌基因组验证classify_wf管道
30.     check_install -> 验证第三方程序和GTDB参考包

复制代码

对于每一个模块可以单独检察help，比方最重要的classify_wf工作流：

1. $ gtdbtk classify_wf -h
2. usage: gtdbtk classify_wf (--genome_dir GENOME_DIR | --batchfile BATCHFILE) --out_dir OUT_DIR
3.                           (--skip_ani_screen | --mash_db MASH_DB) [--no_mash] [--mash_k MASH_K] [--mash_s MASH_S]
4.                           [--mash_v MASH_V] [--mash_max_distance MASH_MAX_DISTANCE] [-f] [-x EXTENSION]
5.                           [--min_perc_aa MIN_PERC_AA] [--prefix PREFIX] [--genes] [--cpus CPUS] [--pplacer_cpus PPLACER_CPUS]
6.                           [--force] [--scratch_dir SCRATCH_DIR] [--write_single_copy_genes] [--keep_intermediates]
7.                           [--min_af MIN_AF] [--tmpdir TMPDIR] [--debug] [-h]
9. mutually exclusive required arguments:
10.   --genome_dir GENOME_DIR 包含FASTA格式基因组文件的目录
11.   --batchfile BATCHFILE 描述基因组的文件路径 - 2或3列的制表符分隔文件（FASTA文件、基因组ID、翻译表[可选]）
13. required named arguments:
14.   --out_dir OUT_DIR     输出文件的目录
16. mutually exclusive required arguments:
17.   --skip_ani_screen     跳过ani筛选步骤，使用mash和skani对基因组进行分类。（默认：False）
18.   --mash_db MASH_DB     保存/读取Mash参考草图数据库（.msh）的路径
20. optional Mash arguments:
21.   --no_mash             跳过使用Mash进行的基因组预过滤（默认：False）
22.   --mash_k MASH_K       k-mer大小[1-32]（默认：16）
23.   --mash_s MASH_S       最大非冗余哈希数（默认：5000）
24.   --mash_v MASH_V       保留的最大p值[0-1]（默认：1.0）
25.   --mash_max_distance MASH_MAX_DISTANCE 选择潜在GTDB基因组作为用户基因组代表的最大Mash距离。（默认：0.15）
27. optional arguments:
28.   -f，--full_tree       在分类步骤中使用未分割的细菌树；这是原始的GTDB-Tk方法（版本<2），
29.                         需要超过320 GB的RAM来加载参考树（默认：False）
30.   -x，--extension EXTENSION 要处理的文件扩展名，gz = gzipped（默认：fna）
31.   --min_perc_aa MIN_PERC_AA
32.                         排除在MSA中AA百分比不足此值的基因组（包含界限）（默认：10）
33.   --prefix PREFIX       所有输出文件的前缀（默认：gtdbtk）
34.   --genes               表示输入文件包含预测的蛋白质（作为氨基酸），跳过基因调用。
35.                         警告：此标志将跳过ANI比较步骤（ani筛选和分类）。（默认：False）
36.   --cpus CPUS           要使用的CPU数量（默认：1）
37.   --pplacer_cpus PPLACER_CPUS pplacer放置过程中使用的CPU数量
38.   --force               如果单个基因组发生错误，则继续处理（默认：False）
39.   --scratch_dir SCRATCH_DIR
40.                         通过写入磁盘来减少pplacer内存使用（较慢）
41.   --write_single_copy_genes
42.                         输出未对齐的单拷贝标记基因（默认：False）
43.   --keep_intermediates  在最终目录中保留中间文件（默认：False）
44.   --min_af MIN_AF       将基因组分配到物种簇的最小对齐比例（默认：0.5）
45.   --tmpdir TMPDIR       指定临时文件的备用目录（默认：/tmp）
46.   --debug               创建用于调试目的的中间文件（默认：False）

复制代码

classify_wf 注释物种

classify_wf全称为 Classify workflow，即分类工作流程。
分类工作流程由四个步骤构成：ani_screen（v2.4.0开始加入）, identify, align, classify.

• ani_screen 步骤将用户基因组与由所有 GTDB 代表性基因组构成的 Mash 数据库举行比力，然后利用 skani 验证最佳的 mash 掷中。利用 skani 分类的用户基因组不会运行管道的其余部分（识别、对齐、分类），并在摘要文件中报告。
  
• identify步骤利用 Prodigal 调用基因，并利用 HMM 模子和 HMMER 包来识别用于体系发育推断的 120 个细菌和 53 个古细菌标志基因 (Parks et al., 2018)。通过将标志基因与其各自的 HMM 模子举行比对来获得多重序列比对 (MSA)。
  
• align步骤毗连比对的标志基因，并将毗连的 MSA 过滤为约莫 5,000 个氨基酸。
  
• 末了，classify步骤利用 pplacer 查找 GTDB-Tk 参考树中每个基因组的最大似然位置。 GTDB-Tk 根据每个基因组在参考树中的位置、相对进化分歧和/或与参考基因组的平均核苷酸同一性 (ANI) 对每个基因组举行分类。
  

1. # 输入基因组目录
2. gtdbtk classify_wf --genome_dir genomes/ --out_dir classify_wf_out --cpus 16

复制代码

可以根据下图估计自己的任务运行时间，评估资源分配：

output

GTDB-Tk 在classify_wf管道的各个阶段输出许多文件，具体格式和示例可以在官网检察：
https://ecogenomics.github.io/GTDBTk/files/index.html
结果目录(classify_wf_out/output)中重要结果文件 (ar53为古菌相关文件，bac120为细菌相关文件)：

• classify/gtdbtk.ar53.classify.tree 树文件，Newick格式，可以用于体系发育树的绘制。
  
• classify/gtdbtk.ar53.summary.tsv 分类结果文件，重要结果，可以用于体系发育树的注释。
  
• identify/gtdbtk.ar53.markers_summary.tsv 利用标志基因情况
  
• align/gtdbtk.ar53.msa.fasta.gz 多序列对齐文件(输入+参考)
  
• align/gtdbtk.ar53.user_msa.fasta.gz 多序列对齐文件(输入)
  

绘制体系发育树可以参考我之前写的R绘制精致的进化树（基础）和R绘制精致的进化树（进阶），大概利用iPhylo网站来交互式绘图：iPhylo 生成并绘制精致的分类树。
分步示例

获取数据

此示例将利用以下两个基因组：基因组A：GCF_003947435.1 [GTDB / NCBI]，基因组B：GCA_002011125.1 [GTDB / NCBI]

1. # 获取基因组。
2. mkdir -p gtdbtk_test_steps/genomes
3. wget -q https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/003/947/435/GCF_003947435.1_ASM394743v1/GCF_003947435.1_ASM394743v1_genomic.fna.gz -O gtdbtk_test_steps/genomes/genome_a.fna.gz
4. wget -q https://ftp.ncbi.nlm.nih.gov/genomes/all/GCA/002/011/125/GCA_002011125.1_ASM201112v1/GCA_002011125.1_ASM201112v1_genomic.fna.gz -O gtdbtk_test_steps/genomes/genome_b.fna.gz

复制代码

基因判断 (identify)

可以通过classify_wf单步完成，但是偶然候分步流程更恰当处理大数据（可以自行分配每一步的资源）：

1. ls -l gtdbtk_test_steps/genomes
2. # –extension，gz是基因组文件的后缀名
3. gtdbtk identify --genome_dir gtdbtk_test_steps/genomes --out_dir gtdbtk_test_steps/identify --extension gz --cpus 2

复制代码

1. [2024-05-20 15:25:20] INFO: GTDB-Tk v2.3.2
2. [2024-05-20 15:25:20] INFO: gtdbtk identify --genome_dir gtdbtk_test_steps/genomes --out_dir gtdbtk_test_steps/identify --extension gz --cpus 2
3. [2024-05-20 15:25:20] INFO: Using GTDB-Tk reference data version r214: /data/home/jianglab/share/GTDB_r214/release214/
4. [2024-05-20 15:25:20] INFO: Identifying markers in 2 genomes with 2 threads.
5. [2024-05-20 15:25:20] TASK: Running Prodigal V2.6.3 to identify genes.
6. [2024-05-20 15:25:30] INFO: Completed 2 genomes in 9.31 seconds (4.65 seconds/genome).
7. [2024-05-20 15:25:30] TASK: Identifying TIGRFAM protein families.
8. [2024-05-20 15:25:38] INFO: Completed 2 genomes in 7.16 seconds (3.58 seconds/genome).
9. [2024-05-20 15:25:38] TASK: Identifying Pfam protein families.
10. [2024-05-20 15:25:38] INFO: Completed 2 genomes in 0.59 seconds (3.37 genomes/second).
11. [2024-05-20 15:25:38] INFO: Annotations done using HMMER 3.3.2 (Nov 2020).
12. [2024-05-20 15:25:38] TASK: Summarising identified marker genes.
13. [2024-05-20 15:25:38] INFO: Completed 2 genomes in 0.03 seconds (75.34 genomes/second).
14. [2024-05-20 15:25:38] INFO: Done.

复制代码

耗时40s，内存消耗131.76MB。
产生idetity/文件夹，关注摘要文件gtdbtk_test_steps/identify/identify/gtdbtk.ar53.markers_summary.tsv，此中详细先容了从古细菌 53 或细菌 120 标志组中识别出的标志。
还可以在每个基因组对应的中间文件目录下找到被调用的基因和标志信息:

1. ls gtdbtk_test_steps/identify/identify/intermediate_results/marker_genes/genome_a.fna/

复制代码

基因比对 (align)

比对步骤将对齐所有已识别的标志，确定最可能的域并输出串联的MSA:

1. gtdbtk align --identify_dir gtdbtk_test_steps/identify --out_dir gtdbtk_test_steps/align --cpus 2

复制代码

1. [2024-05-20 15:32:01] INFO: GTDB-Tk v2.3.2[2024-05-20 15:32:01] INFO: gtdbtk align --identify_dir gtdbtk_test_steps/identify --out_dir gtdbtk_test_steps/align --cpus 2
2. [2024-05-20 15:32:01] INFO: Using GTDB-Tk reference data version r214: /data/home/jianglab/share/GTDB_r214/release214/[2024-05-20 15:32:03] INFO: Aligning markers in 2 genomes with 2 CPUs.[2024-05-20 15:32:04] INFO: Processing 2 genomes identified as archaeal.[2024-05-20 15:32:06] INFO: Read concatenated alignment for 4,416 GTDB genomes.[2024-05-20 15:32:06] TASK: Generating concatenated alignment for each marker.[2024-05-20 15:32:06] INFO: Completed 2 genomes in 0.02 seconds (130.19 genomes/second).[2024-05-20 15:32:06] TASK: Aligning 52 identified markers using hmmalign 3.3.2 (Nov 2020).[2024-05-20 15:32:06] INFO: Completed 52 markers in 0.43 seconds (121.60 markers/second).[2024-05-20 15:32:06] TASK: Masking columns of archaeal multiple sequence alignment using canonical mask.[2024-05-20 15:32:10] INFO: Completed 4,418 sequences in 3.59 seconds (1,231.58 sequences/second).[2024-05-20 15:32:10] INFO: Masked archaeal alignment from 13,540 to 10,135 AAs.[2024-05-20 15:32:10] INFO: 0 archaeal user genomes have amino acids in <10.0% of columns in filtered MSA.[2024-05-20 15:32:10] INFO: Creating concatenated alignment for 4,418 archaeal GTDB and user genomes.[2024-05-20 15:32:12] INFO: Creating concatenated alignment for 2 archaeal user genomes.[2024-05-20 15:32:12] INFO: Done.

复制代码

耗时10s，内存忽略不计（因为这里只有两个genome，输入的多了还是要思量这步的内存的）
产生align/文件夹，如果基因组识别出的标志数量较少，它将被排除在这一步的分析之外。如果是这种情况，将会出现Warning。
根据基因组Domain的差别，将出现 ar53 或 bac120 前缀文件。

• gtdbtk.ar53.msa.fasta.gz：此中包含用户基因组和 GTDB基因组的 MSA。
  
• gtdbtk.ar53.user_msa.fasta.gz：仅包含用户基因组
  

物种分类 (classify)

分类步骤会将基因组放入参考树中，然后确定其最可能的分类。

1. gtdbtk classify --genome_dir gtdbtk_test_steps/genomes --align_dir gtdbtk_test_steps/align --out_dir gtdbtk_test_steps/classify -x gz --cpus 2 --skip_ani_screen

复制代码

1. [2024-05-20 15:42:48] INFO: GTDB-Tk v2.3.2[2024-05-20 15:42:48] INFO: gtdbtk classify --genome_dir gtdbtk_test_steps/genomes --align_dir gtdbtk_test_steps/align --out_dir gtdbtk_test_steps/classify -x gz --cpus 2 --skip_ani_screen
2. [2024-05-20 15:42:48] INFO: Using GTDB-Tk reference data version r214: /data/home/jianglab/share/GTDB_r214/release214/[2024-05-20 15:42:49] TASK: Placing 2 archaeal genomes into reference tree with pplacer using 2 CPUs (be patient).[2024-05-20 15:42:49] INFO: pplacer version: v1.1.alpha19-0-g807f6f3[2024-05-20 15:48:46] INFO: Calculating RED values based on reference tree.[2024-05-20 15:48:48] TASK: Traversing tree to determine classification method.[2024-05-20 15:48:48] INFO: Completed 2 genomes in 0.00 seconds (8,621.39 genomes/second).[2024-05-20 15:48:48] TASK: Calculating average nucleotide identity using FastANI (v1.33).[2024-05-20 15:48:52] INFO: Completed 6 comparisons in 2.69 seconds (2.23 comparisons/second).[2024-05-20 15:48:52] INFO: 2 genome(s) have been classified using FastANI and pplacer.[2024-05-20 15:48:52] INFO: Note that Tk classification mode is insufficient for publication of new taxonomic designations. New designations should be based on one or more de novo trees, an example of which can be produced by Tk in de novo mode.[2024-05-20 15:48:52] INFO: Done.

复制代码

耗时368s，内存消耗45GB（因为这里只有两个古菌genome，如果还有细菌的话应该分配100GB左右）
产生align/文件夹，摘要文件gtdbtk.ar53.summary.tsv和包含这些基因组的参考树gtdbtk.ar53.classify.tree。
注意事项

• 为什么 GTDB-Tk 和 NCBI 或 Silva 分类名称之间的命名体系存在差异？
  

GTDB-Tk 利用 GTDB 分类法 (https://gtdb.ecogenomic.org/)。此分类法与NCBI 和 Silva 类似，但不完全相同。
在许多情况下，GTDB分类法更严格地遵循等级后缀的命名规则，好比GTDB有 Nitrospirota 而没有Nitrospirae。

• 可以将细菌树和古细菌树合并成一棵树吗？
  

细菌树和古细菌树是从差别的标志基因推断出来的。现在，这些树的精确生根仍然是一个开放的研究范畴。
GTDB-Tk 不提供合并树的工具，但可以通过操纵 Newick文件来人为地组合它们。

• GTDB-Tk 怎样利用平均核苷酸同一性（ANI）验证物种分配？
  

GTDB-Tk 利用 skani（在 v2.3.2 之前一直利用 fastANI）来估计基因组之间的ANI。
仅当基因组之间的 ANI 在物种 ANI 限定半径内（通常为 95%）且比对分数(AF)>=0.5时，查询基因组才被分类为与参考基因组属于同一物种。
在某些情况下，查询基因组的体系发育位置可能不支持物种分配。

• 互斥选项 --mash_db 和 --skip_ani_screen 之间有什么区别？
  

从 GTDB-Tk v2.2+ 开始，classify_wf和分类函数必要额外的参数才能运行：--mash_db 或 --skip_ani_screen。
利用新版本的 GTDB-Tk v2.4.0，分类管道（classify_wf 和 classify）的第一阶段是将所有效户基因组与所有参考基因组举行比力，并在可能的情况下根据ANI 匹配对它们举行注释。
利用 --mash_db 选项将向 GTDB-Tk 指示 ANI筛选所需的草图 Mash 数据库的路径。
如果没有可用的数据库（即这是第一次运行分类），则 --mash_db选项将绘制一个可用于后续调用的新 Mash 数据库。 --skip_ani_screen 选项将跳过预筛选步骤并对所有基因组举行分类，类似于之前版本的 GTDB-Tk。

关注公众号 ‘bio llbug’,获取最新推送。

免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！更多信息从访问主页：qidao123.com:ToB企服之家，中国第一个企服评测及商务社交产业平台。