《Frontiers in Plant Science》：Genome-wide subgenome-resolved analysis validates chromosome 4 differentiation and prioritizes introgressed Coffea arabica accessions
编辑推荐：
小粒咖啡（Coffea arabica）中源自帝汶杂种（Timor hybrid）的染色体4渐渗已被确立，但在明确的亚基因组过滤下，种质优先排序的稳健性以及栽培–渐渗分化在canephora来源亚基因组（sgC）与eugenioides来源亚基因组（sgE）间
小粒咖啡（Coffea arabica）中源自帝汶杂种（Timor hybrid）的染色体4渐渗已被确立，但在明确的亚基因组过滤下，种质优先排序的稳健性以及栽培–渐渗分化在canephora来源亚基因组（sgC）与eugenioides来源亚基因组（sgE）间的相对强度仍不清楚。研究人员利用严格contig水平亚基因组过滤、仅限Arabica群体结构分析、SNP面板敏感性检验、全基因组分化扫描、置换检验及直接序列比对，重新分析了44份咖啡种质的公共基因组资源。群体结构与种质排序在不同标记密度和随机种子下均稳定，且始终保留相同的6份渐渗参照。染色体4在两个亚基因组中均排名第一，sgC信号强而sgE信号明显更弱；独立基线面板置换检验支持两个染色体4关联信号。直接比对支持其与预期染色体4假分子对应，同时显示来源覆盖不完全且精确边界未解决。比对支持区块含88个sgC与62个sgE临时防御、信号与调控相关注释。这些结果为优先筛选渐渗种质及表型关联咖啡研究中候选染色体4区域提供了全基因组定量验证框架，但不暗示等价渐渗、精确liftover或因果抗性基因。
研究背景与问题提出
小粒咖啡（Coffea arabica）为异源四倍体，由近缘于C. canephora与C. eugenioides的祖先杂交形成，现代种质携带canephora来源亚基因组（sgC）与eugenioides来源亚基因组（sgE）。其近期起源、少数栽培谱系扩散及后续瓶颈导致遗传基础异常狭窄，制约了气候变化与病害压力下的改良。源自C. arabica–C. canephora杂交的Híbrido de Timor及相关材料向Catimor、Sarchimor等育种背景贡献了抗病多样性，但可能伴随品质或农艺权衡。Saloj?rvi等人（2024）确立了Timor衍生的染色体4渐渗区域主要位于sgC并与抗性基因空间关联。然而在异源多倍体中，同源交换、结构变异与分析错配会扭曲亚基因组推断；原研究未量化四个互补问题：严格contig水平过滤后染色体4是否仍为最强信号；结果在不同SNP面板与种子下是否稳定；哪些渐渗种质可重复优先；临时注释是否获直接序列比对而非比例坐标缩放支持。本研究即在明确亚基因组过滤下提供定量验证层。
主要技术方法概述
研究人员复用Saloj?rvi等人（2024）公共44份种质（39份C. arabica分栽培、渐渗、野生，加3份C. canephora与2份C. eugenioides参照）的sgC/sgE双VCF与ET-39组装。方法包括：contig水平严格过滤仅留label_sg_C/sg_E的11条假分子；双等位SNP位点检出率≥0.80、MAF≥0.05；固定12,000 SNP水库抽样（seed 20250416）作等密度基线，并测6k–48k多面板五种子稳定性；仅Arabica的39份PCA定群内结构；全合格变异做全基因组|ΔAF|扫描与1 Mb窗口统计；严格基线内栽培/渐渗标签500次置换算经验P；前5主成分算种质亲和分数排序；将基线染色体4来源区间分200 kb块用BLAST+ megablast比对ET-39两亚基因组，保留≥1 kb、≥85%一致HSP；提取比对支持区块重叠GFF3基因作临时候选注释。
研究结果
3.1 严格亚基因组过滤后仅Arabica群体结构
排除5个二倍体参照后，sgC与sgE的Arabica-only PCA均分离主要渐渗参照与栽培/野生组；全面板与仅Arabica的Arabica距离中位数相关0.860，但栽培–渐渗质心分离在sgC增约31%、sgE增约83%，故主推断用Arabica-only。
3.2 全合格扫描支持强sgC与较弱sgE染色体4信号
染色体4在两亚基因组全基因组扫描均排第一。sgC中染色体4均|ΔAF|=0.3878，46,552个标记|ΔAF|≥0.90，最大均窗14–15 Mb（0.5450）；sgE均|ΔAF|=0.0843，642个极端标记，最大窗11–12 Mb（0.2137）。12,000-SNP基线500置换中，sgC观测0.5738（14–15 Mb）>null 99百分位0.3482（P=0.001996）；sgE观测0.2232（17–18 Mb）>0.1656（P=0.003992）。结论：两亚基因组均现染色体4关联分化，但幅度强不对称，不支持等价渐渗。
3.3 SNP面板敏感性与种质优先可重复
PCA结构跨面板/种子稳定；12,000 SNP为两亚基因组五种子均取染色体4第一的最小密度，6,000不稳尤其sgE，24k/48k不改结论；全合格与12k基线Arabica距离相关sgC 0.9955、sgE 0.9996。6份渐渗参照在所有条件居前六：Híbrido de Timor、Arabusta Pizze、Iapar 59、Oro Azteca、IPR99、Costa Rica 95；亲和分中位Spearman 0.9729。Híbrido de Timor在sgC 22条件均第一（70.66），Arabusta Pizze在sgE第一（44.43）sgC第二；后四份位次稳定。分数只反映相对质心位置非田间抗性。
3.4 直接比对支持染色体级对应与临时注释资源
来源块最高分均落到预期ET-39染色体4；sgC覆盖70.29%加权一致99.643%，sgE覆盖70.39%一致99.653%；源-靶顺序相关sgC 0.9993、sgE 0.9999；最大非4染色体覆盖<4%。不支持单bp liftover。比对支持区块含sgC 900、sgE 692基因模型，其中88、62满足临时候选（防御受体、激酶信号、RPP13-like/RPP8-like、RLK、WNK3、JA2L、NAC、LRR等），不作因果断言。
讨论部分总结
研究不重发现Timor染色体4渐渗，而添定量验证：严格过滤、全合格扫描、独立置换、多面板敏感、Arabica-only定序、稳定排序、直接比对注释。明确过滤必要因公共VCF混标；sgC强信号获全合格与置换双支持，sgE虽第一但效应远小，解为“同对比组较弱染色体4关联群体分化”而非直接C. canephora渐渗或物理连锁。6份种质跨条件稳居前六且背景合理（自然Timor供体、Arabusta四倍杂种、Timor衍生育成种），但分数非表型预测。直接比对以~70%覆盖与近单调顺序优于比例缩放，但30%无覆盖与少量反向块使精确边界未决。局限含继承公共VCF上游调用、操作区间随密度变、二倍体参照影响几何、比对非完整liftover、注释属假设生成。
结论部分翻译
严格亚基因组过滤、全合格全基因组扫描与严格12,000-SNP基线独立置换检验确认染色体4为sgC栽培–渐渗分化主导信号，sgE为较弱但显著信号。结果跨SNP面板与种子稳定，6份主要渐渗参照在所有敏感条件下保留。直接序列比对支持染色体级对应与比对支持临时注释资源，但不支持精确区间边界或因果抗性基因。该分析将强sgC与较弱sgE染色体4信号联系至6份一致优先渐渗参照：Híbrido de Timor与Arabusta Pizze分别提供最强sgC、sgE亲和锚点，Iapar 59、Oro Azteca、IPR99、Costa Rica 95提供稳定衍生育种参照以供未来染色体4单倍型与表型比较。这些排序反映基因组亲和性而非实测田间抗性或农艺表现。论文发表于《Frontiers in Plant Science》。
