不止于双物种对比：手把手教你用TBtools的‘Unlimited Synteny’功能绘制多物种共线性圈图

多物种基因组共线性圈图绘制实战：TBtools高级可视化技巧解析

在基因组比较研究中，共线性分析是揭示物种间进化关系的重要工具。传统双物种比较已无法满足复杂进化研究的需要，而多物种共线性可视化能同时呈现多个基因组的保守区域与重排事件，为理解大规模基因组演化提供更全面的视角。本文将深入解析如何利用TBtools的"Unlimited Synteny Visualization"功能，制作可直接用于学术发表的多物种共线性圈图。

1. 多物种共线性分析的核心准备

1.1 数据获取与标准化处理

获取高质量的基因组注释文件是共线性分析的基础。Ensembl Plants、Phytozome和NCBI是常用的植物基因组数据来源。以拟南芥(Arabidopsis thaliana)、水稻(Oryza sativa)和油菜(Brassica rapa)三物种比较为例，需要准备：

• 基因组序列文件：FASTA格式的基因组组装序列
• 基因注释文件：GFF3格式的基因结构注释
• 共线性文件：MCScanX生成的Collinearity文件

关键预处理步骤：

# 示例：GFF文件标准化处理 grep -v "^#" Athaliana.gff | awk '$3=="gene"' > Athaliana_genes.gff sed -i 's/ID=gene://g' Athaliana_genes.gff

1.2 多物种共线性分析流程

使用TBtools进行多物种共线性分析时，推荐的分步策略：

1. 两两物种间分别运行MCScanX分析
2. 合并各次分析结果生成综合共线性网络
3. 使用"Advanced Circos"模块进行可视化配置

注意：跨物种比较时建议使用蛋白序列而非DNA序列，可提高远缘物种间的共线性检测灵敏度

2. Unlimited Synteny功能深度配置

2.1 输入文件的结构化组织

TBtools的多物种共线性可视化需要特定格式的输入文件。建议创建如下目录结构：

synteny_project/ ├── input/ │ ├── Athaliana.gff │ ├── Osativa.gff │ ├── Brapa.gff │ └── all.collinearity └── config/ └── color_scheme.txt

关键配置文件示例（color_scheme.txt）：

Athaliana:#FF6B6B Osativa:#4ECDC4 Brapa:#45B7D1 synteny_block:#A5FFD6

2.2 可视化参数优化技巧

在"Unlimited Synteny Visualization"界面中，几个关键参数影响最终输出效果：

实际操作中，可通过"Preview"功能实时调整以下元素：

• 染色体分段策略（按长度或基因密度）
• 共线性区块的配色梯度
• 物种标签的位置偏移量

3. 高级可视化效果实现

3.1 复杂进化事件的突出展示

对于全基因组复制(WGD)或大规模染色体重排等事件，可通过分层着色策略增强可视化效果：

1. 使用不同深浅表示古老/近期的共线性区块
2. 对特定进化分支添加标记符号
3. 用虚线连接线表示可能的非共线性同源基因

# 示例：使用Python预处理共线性区块标签 import pandas as pd df = pd.read_csv('collinearity.csv') df['color'] = df['divergence'].apply( lambda x: '#FF0000' if x<0.5 else '#0000FF') df.to_csv('annotated_collinearity.csv', index=False)

3.2 出版级图表的美学优化

学术期刊对图表有严格要求，TBtools输出的SVG/PDF格式可进一步用Inkscape或Adobe Illustrator调整：

• 字体统一：全图使用期刊指定字体家族
• 图例完善：添加比例尺和进化时间标尺
• 标注清晰：用箭头指示关键共线性断裂点
• 分辨率保障：最终导出600dpi以上的TIFF格式

提示：Nature系列期刊推荐使用CMYK色彩模式，线上发表则可保留RGB模式

4. 结果解读与生物学洞见挖掘

4.1 多维度共线性模式分析

通过三物种比较可识别不同层级的保守基因组区域：

• 核心共线性区块：三物种共有的古老保守区域
• 谱系特异性区块：仅两物种共享的较新保守区域
• 物种特有重排：单个物种特有的基因组结构变异

典型分析流程：

1. 统计各染色体对的共线性基因比例
2. 计算共线性区块的Ka/Ks比值分布
3. 关联已知功能基因家族的位置信息

4.2 进化假说验证案例

以十字花科基因组三倍化事件研究为例，多物种共线性图可揭示：

• 拟南芥与油菜共享的古老六倍化痕迹
• 水稻作为外类群的共线性断裂模式
• 特定基因家族（如MADS-box）的复制保留情况

实际操作中，建议结合PhyloNet或Notung等系统发育工具，将共线性结果置于明确的进化框架下解读。