224 lines
6.8 KiB
Markdown
Executable File
224 lines
6.8 KiB
Markdown
Executable File
# Ka/Ks 分析图表解读指南
|
||
|
||
## 1. Ka/Ks 分布直方图 (kaks_distribution.png)
|
||
|
||
### 图表说明
|
||
|
||
```
|
||
Ka/Ks Distribution
|
||
Mean Ka/Ks = 3.1038
|
||
|
||
↑ Count
|
||
2.0│ ┌───┐
|
||
│ │ │
|
||
1.5│ │ │
|
||
│ │ │
|
||
1.0│ ┌───┐ ┌───┐ ┌───┐ ┌───┐│ │┌───┐ ┌───┐
|
||
│ │ │ │ │ │ │ │ ││ ││ │ │ │
|
||
0.5│ │ │ │ │ │ │ │ ││ ││ │ │ │
|
||
│ │ │ │ │ │ │ │ ││ ││ │ │ │
|
||
0.0└──┴───┴─┴───┴─┴───┴─────┴───┴┴───┴┴───┴─┴───┴──→ Ka/Ks Ratio
|
||
│ 2 │ 3 │ 4
|
||
│ │ ↑ │
|
||
│ │ 绿色虚线 │
|
||
│ │ (平均值) │
|
||
↑
|
||
红色虚线
|
||
(中性选择线)
|
||
Ka/Ks = 1
|
||
```
|
||
|
||
### 关键元素标注
|
||
|
||
| 元素 | 位置 | 含义 |
|
||
|------|------|------|
|
||
| **红色虚线** | Ka/Ks = 1 | 中性选择线。Ka/Ks = 1 表示基因经历中性进化（同义和非同义替换率相同） |
|
||
| **绿色虚线** | Ka/Ks = 3.1038 | 平均 Ka/Ks 值。所有基因对的 Ka/Ks 平均值 |
|
||
| **蓝色柱子** | 分布区间 | 显示有多少基因对的 Ka/Ks 值落在该区间 |
|
||
|
||
### 选择压力判断标准
|
||
|
||
```
|
||
Ka/Ks < 1 → 纯化选择 (Purifying selection)
|
||
有害突变被淘汰，基因功能保守
|
||
|
||
Ka/Ks = 1 → 中性选择 (Neutral selection)
|
||
突变随机固定，无选择压力
|
||
|
||
Ka/Ks > 1 → 正选择 (Positive selection)
|
||
有利突变被保留，适应性进化
|
||
```
|
||
|
||
### 你的结果解读
|
||
|
||
**Ft.SPL vs Ath.SPL 分析结果：**
|
||
|
||
- ✅ **所有基因对的 Ka/Ks > 1**（都在红色虚线右侧）
|
||
- ✅ **平均 Ka/Ks = 3.1038**（远高于 1）
|
||
- ✅ **分布范围：2.0 - 4.1**
|
||
|
||
**结论**：SPL 基因家族在苦荞和拟南芥之间经历了**强烈的正选择**，表明这些基因在物种分化过程中可能参与了适应性进化。
|
||
|
||
---
|
||
|
||
## 2. Ka vs Ks 散点图
|
||
|
||
### 图表说明
|
||
|
||
```
|
||
↑ Ka (非同义替换)
|
||
│
|
||
600│ ●
|
||
│ ●
|
||
400│ ●
|
||
│ ●
|
||
200│ ●
|
||
│ ●
|
||
0├─────┼─────┼─────┼─────┼─────┼→ Ks (同义替换)
|
||
0 500 1000 1500 2000
|
||
|
||
红线: Ka = Ks (斜率为1)
|
||
颜色: 蓝→绿→红 (Ka/Ks 从小到大)
|
||
```
|
||
|
||
### 关键元素标注
|
||
|
||
| 元素 | 含义 |
|
||
|------|------|
|
||
| **X轴 (Ks)** | 同义替换率（不改变氨基酸的突变） |
|
||
| **Y轴 (Ka)** | 非同义替换率（改变氨基酸的突变） |
|
||
| **红点** | Ka/Ks > 1，正选择 |
|
||
| **绿点** | Ka/Ks ≈ 1，中性选择 |
|
||
| **蓝点** | Ka/Ks < 1，纯化选择 |
|
||
| **红线** | Ka = Ks 参考线 |
|
||
|
||
### 解读方法
|
||
|
||
```
|
||
点在红线上方 → Ka > Ks → Ka/Ks > 1 → 正选择
|
||
点在红线下方 → Ka < Ks → Ka/Ks < 1 → 纯化选择
|
||
点在红线上 → Ka = Ks → Ka/Ks = 1 → 中性选择
|
||
```
|
||
|
||
---
|
||
|
||
## 3. Ks 分布图
|
||
|
||
### 图表说明
|
||
|
||
```
|
||
↑ Count
|
||
│
|
||
3.0│ ┌───┐
|
||
│ │ │
|
||
2.0│ ┌───┐│ │┌───┐
|
||
│ │ ││ ││ │
|
||
1.0│┌──┤ ││ ││ │┌──┐
|
||
││ │ ││ ││ ││ │
|
||
0.0└┴──┴───┴┴───┴┴───┴┴──┴──→ Ks
|
||
0 500 1000 1500 2000
|
||
|
||
Ks 越大 → 分歧时间越长
|
||
```
|
||
|
||
### 关键元素标注
|
||
|
||
| 元素 | 含义 |
|
||
|------|------|
|
||
| **X轴 (Ks)** | 同义替换数 |
|
||
| **Y轴 (Count)** | 基因对数量 |
|
||
| **峰值位置** | 最普遍的 Ks 值 |
|
||
|
||
### 生物学意义
|
||
|
||
Ks 是**分子钟**的基础：
|
||
- Ks 随时间线性累积
|
||
- Ks 不受选择压力影响（同义突变不改变蛋白功能）
|
||
- 用于计算物种分歧时间
|
||
|
||
---
|
||
|
||
## 4. 分歧时间分布图
|
||
|
||
### 图表说明
|
||
|
||
```
|
||
↑ Count
|
||
│
|
||
3.0│ ┌───┐
|
||
│ │ │
|
||
2.0│ ┌───┐│ │┌───┐
|
||
│ │ ││ ││ │
|
||
1.0│┌──┤ ││ ││ │┌──┐
|
||
││ │ ││ ││ ││ │
|
||
0.0└┴──┴───┴┴───┴┴───┴┴──┴──→ Divergence Time (Mya)
|
||
0 50 100 150 200
|
||
(百万年前)
|
||
```
|
||
|
||
### 关键元素标注
|
||
|
||
| 元素 | 含义 |
|
||
|------|------|
|
||
| **X轴 (Mya)** | 百万年前 (Million years ago) |
|
||
| **Y轴 (Count)** | 基因对数量 |
|
||
| **峰值** | 最可能的分歧时间 |
|
||
|
||
### 计算公式
|
||
|
||
```
|
||
分歧时间 (Mya) = Ks / (2 × 中性替换率) / 1,000,000
|
||
|
||
其中：
|
||
- Ks = 同义替换数
|
||
- 中性替换率 = 7.08e-9 (默认值，来自拟南芥)
|
||
- 除以 2 是因为两个物种都在积累突变
|
||
```
|
||
|
||
### 你的结果
|
||
|
||
**Ft vs Ath 分歧时间：约 145 Mya**
|
||
|
||
- 对应地质年代：**早白垩纪 (Early Cretaceous)**
|
||
- 这是十字花科和蓼科植物分化的大致时间
|
||
|
||
---
|
||
|
||
## 5. 结果汇总表
|
||
|
||
| 指标 | 数值 | 生物学意义 |
|
||
|------|------|-----------|
|
||
| **平均 Ka/Ks** | 3.1038 | 强烈的正选择 |
|
||
| **Ka/Ks 范围** | 2.07 - 4.13 | 所有基因对都经历正选择 |
|
||
| **平均 Ks** | 1450.2 | 较高的同义替换率 |
|
||
| **平均分歧时间** | ~145 Mya | 早白垩纪分化 |
|
||
| **基因对数量** | 9对 | 有效的同源基因对 |
|
||
|
||
---
|
||
|
||
## 6. 生物学结论
|
||
|
||
### SPL 基因家族的进化特征
|
||
|
||
1. **强烈的正选择信号**
|
||
- 所有 Ka/Ks > 1，平均 3.1
|
||
- 表明 SPL 基因在物种适应中发挥重要作用
|
||
|
||
2. **古老的分歧时间**
|
||
- ~145 Mya 的分歧时间与十字花科-蓼科分化时间一致
|
||
- 说明 SPL 基因家族在核心真双子叶植物分化前已存在
|
||
|
||
3. **功能分化**
|
||
- 正选择可能导致 SPL 基因在苦荞中获得了新功能
|
||
- 可能与苦荞特有的生物学特性（如黄酮类合成）相关
|
||
|
||
---
|
||
|
||
## 7. 参考文献
|
||
|
||
1. Yang, Z., & Nielsen, R. (2000). Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. *Molecular Biology and Evolution*, 17(1), 32-43.
|
||
|
||
2. Nei, M., & Gojobori, T. (1986). Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. *Molecular Biology and Evolution*, 3(5), 418-426.
|
||
|
||
3. 中性替换率参考：Koch, M. A., et al. (2000). Molecular systematics of the Brassicaceae: evidence from coding plastidic matK and nuclear Chs sequences. *American Journal of Botany*, 87(4), 534-544.
|