大豆 FIP37/WTAP 同源基因家族分析（V4 版本）
四个 paralog：Glyma.04G186400 / Glyma.05G040200 / Glyma.06G179400 / Glyma.17G086600
基因组版本：Wm82.gnm4.ann1（V4）
数据来源：SoyBase + GlycineMine + CoNekT（表达数据仅 gnm2），检索日期 2025-06-23
关键发现
四个基因编码同一蛋白：FKBP12-interacting protein of 37 kDa（FIP37）——拟南芥 AT3G54170 的大豆同源基因。
这是 m6A RNA 甲基化 writer 复合体的核心组分（对应人类 WTAP），功能极其重要且高度保守。
一、基因基本信息（gnm4.ann1）
| 属性 | Glyma.04G186400 | Glyma.05G040200 | Glyma.06G179400 | Glyma.17G086600 |
|---|---|---|---|---|
| 完整 ID | glyma.Wm82.gnm4.ann1.Glyma.04G186400 | glyma.Wm82.gnm4.ann1.Glyma.05G040200 | glyma.Wm82.gnm4.ann1.Glyma.06G179400 | glyma.Wm82.gnm4.ann1.Glyma.17G086600 |
| 染色体 | Gm04 | Gm05 | Gm06 | Gm17 |
| 位置 (gnm4) | 44,477,203–44,483,937 | 3,703,766–3,711,125 | 15,184,806–15,192,853 | 6,689,072–6,695,659 |
| 链 | (+) | (+) | (-) | (-) |
| 基因长度 | 6,735 bp | 7,360 bp | 8,048 bp | 6,588 bp |
| Child Features | 2 | 4 | 2 | 4 |
| Proteins | 1 | 2 | 1 | 2 |
| Data Sets | 2 | 2 | 2 | 2 |
| Pan Gene Sets | 1 | 1 | 1 | 1 |
| Gene Family | 1 | 1 | 1 | 1 |
| Ontology Annotations | 2 | 2 | 2 | 2 |
| Protein Domains | 0 | 0 | 0 | 0 |
| Overlapping Features | 22 | 35 | 29 | 38 |
| Publications | 1 | 1 | 1 | 1 |
| 蛋白 | FKBP12-interacting protein of 37 kDa-like isoform X1 | |||
| 基因家族 | legfed_v1_0.L_2X1ZY8 | |||
| 物种 | Glycine max cv. Williams 82 | |||
| 丰度 | ⚠️ 低丰度（gnm2 CoNekT） | ⚠️ 低丰度 | ⚠️ 低丰度 | ⚠️ 低丰度 |
Gnm2 Vs Gnm4 关键差异
| 属性 | gnm2 | gnm4 | 变化 |
|---|---|---|---|
| 04G186400 位置 | Gm04:45,656,968–45,664,438 | Gm04:44,477,203–44,483,937 | 偏移 ~1.18 Mb |
| 04G186400 长度 | 7,471 bp | 6,735 bp | −736 bp |
| 05G040200 位置 | Gm05:3,582,202–3,589,663 | Gm05:3,703,766–3,711,125 | 偏移 ~121 kb |
| 06G179400 位置 | Gm06:15,233,148–15,241,225 | Gm06:15,184,806–15,192,853 | 偏移 ~48 kb |
| 17G086600 位置 | Gm17:6,689,133–6,695,722 | Gm17:6,689,072–6,695,659 | 几乎一致 |
| 04G186400 Transcripts | 3 → 3 proteins | 2 → 1 protein | 基因模型更保守 |
| 05G040200 Transcripts | 3 → 3 proteins | 4 → 2 proteins | 转录本增加，蛋白减少 |
| 06G179400 Transcripts | 3 → 3 proteins | 2 → 1 protein | 基因模型更保守 |
| 17G086600 Transcripts | 2 → 2 proteins | 4 → 2 proteins | 转录本增加 |
⚠️ 表达数据（CoNekT）仅在 gnm2 版本可用。gnm4 版本无独立的 CoNekT 表达数据，但基因功能注释（GO/InterPro）跨版本一致。
二、拟南芥同源基因功能（AtFIP37 / AT3G54170）
| 功能维度 | 描述 | 文献 |
|---|---|---|
| 分子功能 | m6A RNA 甲基转移酶复合体核心亚基（MTA/MTB/FIP37/VIR/HAKAI） | Shen et al., 2016 Dev Cell |
| 人类同源 | WTAP（Wilms’ Tumor 1-Associating Protein） | PMC419804 |
| 亚细胞定位 | 细胞核，核散斑（nuclear speckles）富集 | Vespa et al., 2004 Plant Physiol |
| 发育功能 | 茎尖分生组织（SAM）维持；fip37 突变体 SAM 过度增殖 | Shen et al., 2016 |
| 分子表型 | fip37-4 突变体中 m6A 水平降至野生型的 ~20% | Shen et al., 2016 |
| 免疫表型 | m6A writer 突变体（fip37/vir/mta）触发自身免疫：PR1 上调 300+ 倍，对 P. syringae 和 B. cinerea 抗性增强 | PLOS Genetics, 2025 |
| 蛋白互作 | 与 FKBP12 直接物理互作（酵母双杂 + pull-down），与 MTA 互作 | Vespa et al., 2004 |
核心生物学意义：FIP37 通过介导干细胞关键基因（如 STM、WUS）mRNA 的 m6A 修饰来调控其稳定性，维持分生组织稳态。m6A writer 活性丧失 → 防御基因去抑制 → 自身免疫激活。
三、表达谱（CoNekT Gnm2 TPM）
⚠️ 仅 gnm2 版本有此数据，gnm4 对应基因的功能注释相同，表达趋势可参考。
3.1 组织表达对比
| 组织 | 04G186400 | 05G040200 | 06G179400 | 17G086600 |
|---|---|---|---|---|
| Apical Meristem | 3.41 | 1.89 | 2.16 | 2.44 |
| Root Tip | 3.84 | 1.10 | 2.43 | 2.61 |
| Flower | 1.97 | 1.04 | 1.15 | 1.30 |
| Root | 1.36 | 0.93 | 1.73 | 1.83 |
| Nodule | 2.14 | 0.54 | 1.07 | 1.42 |
| Leaves | 1.57 | 0.42 | 1.16 | 0.47 |
| Green Pods | 1.98 | 0.22 | 0.41 | 1.50 |
模式：04G186400 表达量最高，05G040200 最低。四个基因均在根尖和茎尖分生组织中表达最高（与拟南芥 FIP37 在 SAM 中的核心功能一致）。
3.2 根毛侵染实验（HA1）
| 条件 | 04G186400 | 05G040200 | 06G179400 | 17G086600 |
|---|---|---|---|---|
| 12HA1_IN_RH | 2.62 ↑ | 0.72 | 1.21 | 0.81 ↓ |
| 12HA1_UN_RH | 1.32 | 0.71 | 1.22 | 1.15 |
| 24HA1_IN_RH | 2.17 ↑ | 0.50 | 1.08 ↓ | 1.80 ↑ |
| 24HA1_UN_RH | 1.29 | 0.49 | 1.50 | 1.10 |
| 48HA1_IN_RH | 3.21 ↓ | 0.97 ↑ | 1.62 | 1.87 ↓ |
| 48HA1_UN_RH | 3.80 | 0.70 | 1.72 | 2.20 |
| 48HA1_Scrip_Root | 1.02 | 1.41 | 1.38 | 0.71 |
↑ = 侵染组高于对照组；↓ = 侵染组低于对照组
四、共表达网络（CoNekT gnm2，Top 5 PCC）
Glyma.04G186400
| 排名 | 共表达基因 | PCC | 功能 |
|---|---|---|---|
| 1 | Glyma.15G073500 | 0.923 | tRNA 甲基转移酶（RNA 修饰→与 FIP37 功能协同） |
| 2 | Glyma.16G022300 | 0.899 | 精氨酸/丝氨酸富集剪接因子 |
| 3 | Glyma.19G191100 | 0.882 | U2AF 剪接因子小亚基（mRNA 加工） |
| 4 | Glyma.17G055300 | 0.867 | 未知蛋白（线粒体/质体定位） |
Glyma.05G040200
| 排名 | 共表达基因 | PCC | 功能 |
|---|---|---|---|
| 1 | Glyma.04G000700 | 0.928 | LIN52 同源（DRM 复合体/细胞周期） |
| 2 | Glyma.06G053600 | 0.921 | 线粒体分布与形态蛋白 |
| 3 | Glyma.19G259600 | 0.919 | DUF4228 未知功能蛋白（植物特有） |
| 4 | Glyma.04G063400 | 0.917 | 转录调控因子（SAP30 结合蛋白） |
| 5 | Glyma.08G195600 | 0.893 | F-box/RNI 超家族蛋白（泛素化） |
Glyma.06G179400
| 排名 | 共表达基因 | PCC | 功能 |
|---|---|---|---|
| 1 | Glyma.U015000 | 0.915 | 假设蛋白 |
| 2 | Glyma.08G346700 | 0.910 | FKBP 型肽基脯氨酰顺反异构酶 ← FIP37 的互作伙伴! |
| 3 | Glyma.07G252700 | 0.908 | SR 剪接因子 33 |
| 4 | Glyma.18G160500 | 0.900 | 双特异性蛋白磷酸酶 |
| 5 | Glyma.06G151500 | 0.887 | ASHH4 组蛋白甲基转移酶（表观遗传） |
Glyma.17G086600
| 排名 | 共表达基因 | PCC | 功能 |
|---|---|---|---|
| 1 | Glyma.13G214200 | 0.920 | BRCT 结构域蛋白（DNA 损伤修复） |
| 2 | Glyma.02G038600 | 0.917 | MATE 外排转运蛋白（药物/毒素转运） |
| 3 | Glyma.12G074700 | 0.912 | NOT 转录复合体亚基 VIP2 |
| 4 | Glyma.08G347800 | 0.909 | Ras GTPase 激活结合蛋白（核转运因子） |
| 5 | Glyma.19G072400 | 0.904 | COP9 信号体复合体亚基 2 |
🔑 最关键的共表达：06G179400 与 Glyma.08G346700（FKBP 型 PPIase）共表达，这正是 FIP37 的已知互作伙伴 FKBP12 的同源基因——在分子层面验证了该基因家族功能的保守性！
五、在生物胁迫中的角色
5.1 核心机制：m6A 修饰与免疫的平衡
1 | 正常状态 m6A writer 突变/抑制 |
5.2 四个 Paralog 的胁迫表达总结（基于 Gnm2 表达数据）
| 基因 | 侵染响应趋势 | gnm4 位置 | 功能推断 |
|---|---|---|---|
| 04G186400 | 早期侵染↑，晚期↓ | Gm04:44.48 Mb | 主要功能拷贝，可能响应侵染初期信号 |
| 05G040200 | 极低表达，48h 侵染略↑ | Gm05:3.70 Mb | 可能是表达量最低的备份拷贝 |
| 06G179400 | 波动较小 | Gm06:15.18 Mb (-) | 稳定表达，与 FKBP 共表达 |
| 17G086600 | 侵染后持续↓ | Gm17:6.69 Mb (-) | 侵染条件下被下调 |
5.3 关键推论
- 这四个基因是 m6A writer 复合体的核心亚基，不是经典的防御基因
- 低丰度特征 + 分生组织偏好 → 主要功能是发育调控（维持干细胞命运）
- m6A writer 活性下降会触发自身免疫 → 侵染条件下表达下调可能是一种防御激活策略
- 拟南芥中 fip37-4 突变体对 P. syringae DC3000 和 B. cinerea 抗性显著增强
- 04G186400 是四个拷贝中表达最高、响应最活跃的 → 优先进行功能验证
- gnm4 相比 gnm2 基因模型更保守（更少的 transcript/protein isoforms），在进行克隆或 CRISPR 设计时应优先参考 gnm4
六、验证建议
- qRT-PCR：在 P. sojae 或 P. syringae 侵染大豆后 0h/12h/24h/48h 检测四个拷贝的表达
- 亚细胞定位：GFP 融合确认核定位（预期：核散斑）
- 酵母双杂：验证与 GmFKBP12 同源蛋白的互作
- VIGS/RNAi 沉默：分别沉默各拷贝，观察发育表型 + 抗病性变化
- m6A-IP-seq：比较侵染前后 m6A 修饰谱变化
- 重点验证 04G186400 和 06G179400
- 如需克隆或设计 gRNA，优先使用 gnm4 坐标
V4 (Gnm4) 完整标识符速查
| 基因 | 完整 gnm4 ID |
|---|---|
| Glyma.04G186400 | glyma.Wm82.gnm4.ann1.Glyma.04G186400 |
| Glyma.05G040200 | glyma.Wm82.gnm4.ann1.Glyma.05G040200 |
| Glyma.06G179400 | glyma.Wm82.gnm4.ann1.Glyma.06G179400 |
| Glyma.17G086600 | glyma.Wm82.gnm4.ann1.Glyma.17G086600 |
数据来源
| 数据库 | URL | 提供数据 |
|---|---|---|
| GlycineMine | https://mines.legumeinfo.org/glycinemine/ | gnm4 基因模型、位置、注释统计 |
| SoyBase | https://www.soybase.org/ | 跨版本基因搜索 |
| CoNekT | https://conekt.legumeinfo.org/ | gnm2 表达谱 + 共表达网络 |
| 拟南芥 FIP37 文献 | Shen et al., 2016 Dev Cell; Vespa et al., 2004 Plant Physiol | 功能与机制 |
| m6A + 免疫 | PLOS Genetics, 2025; PMC10226572, 2023 | 免疫表型 |
整理日期：2025-06-23
工具链：GlycineMine（gnm4 注释）→ CoNekT（gnm2 表达+共表达）→ Exa 文献搜索
