《The Plant Journal》：Telomere-to-telomere genome assembly of the endangered Jacktree (Sinojackia xylocarpa Hu) reveals its centromeric landscape and evolutionary dilemma as a small population
编辑推荐：
许多植物物种的濒危状态通常归因于适合度下降，进而导致自然条件下种群规模的持续缩减。然而，高度重复基因组区域在濒危植物衰退中的作用仍知之甚少。本研究提出了濒危植物秤锤树（Sinojackia xylocarpa Hu）的无缺口端粒到端粒（T2T）参考基因组，其基
许多植物物种的濒危状态通常归因于适合度下降，进而导致自然条件下种群规模的持续缩减。然而，高度重复基因组区域在濒危植物衰退中的作用仍知之甚少。本研究提出了濒危植物秤锤树（Sinojackia xylocarpa Hu）的无缺口端粒到端粒（T2T）参考基因组，其基因组大小为961.94 Mb，N50为78.05 Mb，预测蛋白编码基因36 446个。利用CENH3染色质免疫共沉淀测序（ChIP-seq），研究人员鉴定了全部12条染色体上的着丝粒区域，并揭示秤锤树着丝粒由多种转座元件（TE）家族组成，其中长末端重复（LTR）反转录转座子和LINE元件为主要成分。值得注意的是，1号染色体右臂表现出异常的CENH3富集和类似着丝粒区域的TE积累，暗示形成了着丝粒样染色质结构域。此外，研究人员在自然种群中检测到倒位多态性和倒位杂合子，主要集中于同一染色体臂。这些特征可能与减数分裂异常相关。利用高质量T2T基因组组装，研究人员鉴定了这些小种群面临的基因组侵蚀模式，其特征为与LTR增殖相关的局部重组抑制以及DNA修复基因中有害突变的积累。这种基因组侵蚀形成了进化困境，阻碍了有害等位基因的清除，最终威胁长期适应性和生存。本研究强调了T2T基因组组装在濒危植物保护中的效用。
研究背景方面，小种群不仅是数量描述，更反映种群易受随机事件影响的脆弱状态，形成负反馈循环。小种群受遗传漂变影响，导致遗传多样性快速丧失，伴随近交衰退和有害突变积累，多种内外过程汇聚可形成"灭绝涡旋"（extinction vortex）。基因组侵蚀（genomic erosion）涵盖遗传多样性丧失、遗传负荷增加和适应不良变异积累等一系列遗传威胁，是小种群可能陷入的负反馈循环中的一环。然而，由于端粒、着丝粒及其邻近区域重复序列密度高，历史上完整组装这些区域存在困难，高度重复序列对小种群的进化效应仍知之甚少。近年来T2T组装技术的发展使得完整解析端粒和着丝粒成为可能。秤锤树（Sinojackia xylocarpa Hu）为安息香科（Styracaceae）秤锤树属落叶乔木，分布于长江下游地区，被IUCN列为易危（VU）物种，其分布有限且碎片化，种群均较小。前期研究将濒危原因归因于繁殖能力低下，包括种子休眠延长和减数分裂异常，但秤锤树衰退的遗传原因及当前小种群面临的基因组风险尚未被探索。为此，研究人员组装了秤锤树T2T基因组以探究其小种群的基因组驱动因素，并基于无缺口基因组开展种群水平分析以深入揭示濒危小种群的基因组侵蚀。研究意义在于，这是安息香科首个公开可用的T2T无缺口参考基因组，揭示了濒危植物高度重复区域（如端粒和着丝粒）的精细特征，为更广泛地理解受威胁植物的整体基因组结构奠定了基础。
技术方法方面（不超过250字）：研究整合PacBio HiFi reads（82.53 Gb，85.8×）、ONT超长reads（44.96 Gb，46.3×）和Hi-C reads（93.40 Gb）完成无缺口基因组组装；利用CENH3 ChIP-seq鉴定着丝粒区域；对8个已知种群32个个体进行全基因组重测序；采用PSMC重建种群历史；利用FastEPRR估算重组率；利用SIFT4G预测有害突变；通过SyRi鉴定结构变异。样本队列来源为8个已知秤锤树自然种群。
研究结果方面：第一，秤锤树二倍体T2T基因组组装。通过整合多平台测序数据，研究获得了总长961 939 874 bp的无缺口T2T参考基因组，包含12条染色体，N50为78 050 245 bp，BUSCO完整性达97.30%。第二，基因组注释。预测到36 446个编码基因，基因BUSCO完整性达99.01%，96.24%的基因至少在一個数据库中获得注释；TE覆盖基因组的65.97%，以Gypsy（22.62%）和Copia（8.20%）LTR反转录转座子为主。第三，染色体高度重复序列检测。高重复密度区域集中于前5条染色体，chr1右臂平均重复密度达94.30%；所有染色体均鉴定到端粒；通过CENH3 ChIP-seq鉴定全部12条染色体的着丝粒区域，长度范围为0.94 Mb（chr5）至3.01 Mb（chr1）。第四，比较基因组学分析。16个物种的比较分析鉴定92 380个直系同源基因家族，秤锤树含3 205个特有基因家族，富集于酸性磷酸酶活性及碳代谢等通路；扩张基因家族富集于UDP-糖基转移酶活性和花粉识别等功能。第五，秤锤树基因组进化。系统发育分析揭示安息香科与山茶科亲缘关系最近，分化时间约89.9 Ma；全基因组复制分析推断秤锤树经历两轮WGD事件；鉴定到两个受显著正选择的基因（SxylChr2G042980和SxylChr11G332510）。第六，着丝粒景观。着丝粒区域73.73%由TE组成，LINEs（31.58%）和LTR反转录转座子（23.58%）最为丰富；LTR反转录转座子表现出更强的CENH3富集信号；chr1左右臂呈现不对称CENH3富集模式，右臂富集显著升高且富集区域98.39%位于TE内，提示异位CENH3定位。第七，秤锤树种群遗传学。对32个个体全基因组重测序鉴定67.34 Mb高置信变异；系统发育分析鉴定4个地理上明显的分支；PCA和ADMIXTURE鉴定2个遗传簇及1个混合遗传簇；多数种群LD快速衰减，但PZ种群LD持续性异常；NJ和HM种群FROH最高且核苷酸多样性最低；种群历史显示约20 ka（末次盛冰期）有效种群大小达峰值后持续下降。第八，秤锤树种群染色体倒位。比较两个基因组发现多数染色体存在倒位，chr1右臂PCA将种群分为三个簇，4个个体右臂Ho>0.5，提示存在倒位杂合子个体。第九，基因组重组景观与有害突变积累。全基因组重组率分析揭示广泛区域重组明显受抑；重组率与基因数量正相关（r=0.57），与Gypsy型LTR元件负相关（r=?0.52）；SIFT4G预测68 655个有害突变，在低重组区域（尤其chr1中远端）局部富集；携带有害突变的基因富集于DNA修复（125个基因）、氨基酸代谢过程和DNA损伤应答通路；POLH（SxylChr1G021390）中鉴定到有害突变且未检测到纯合个体；chr2上游端粒倒位区域4个DNA修复基因（RAD51、RFA1、RecA、RAD7）中鉴定到有害突变。
讨论部分总结：在染色体异常方面，chr1右臂高度重复序列（94.30%）及异常CENH3富集提示异位CENH3定位和着丝粒样结构域形成，可能诱导新生着丝粒、干扰正常着丝粒功能，导致染色体分离错误和基因组不稳定；倒位杂合子在减数分裂配对时形成倒位环，臂内倒位杂合子内单交换产生双着丝粒和无着丝粒染色单体，导致染色质桥和无着丝粒片段形成，损害配子形成并降低个体适合度。染色体倒位源于DNA双链断裂（DSBs），高密度重复序列可作为内源性复制或重组压力源诱导DSBs，而DNA修复基因的有害突变（如RAD51）可能降低同源重组修复效率，使DSB修复偏向易错的非同源末端连接（NHEJ）途径，直接导致倒位产生。在小种群自增强进化困境方面，研究人员鉴定出重组受抑的独特基因组架构，chr1特定区域同线性片段少、基因密度低、重复密度高、种群变异显著降低且有害突变比例升高；重组率与Gypsy型LTR密度显著负相关，LTR增殖可能通过促进异染色质形成和改变染色质架构直接抑制局部重组；非同义/同义替换比高于预期，与有害突变积累一致。LTR驱动的重组抑制、缺陷性DNA修复导致倒位和突变积累、小种群效应降低纯化选择三种因素相互强化，形成自增强进化困境：小种群降低自然选择效力，强烈区域重组抑制和DNA修复基因突变导致有害突变积累，进而损害染色体断裂和DNA损伤修复，降低适合度，形成恶性循环。在保护意义方面，NJ和HM种群近交系数最高、遗传多样性最低，应优先保护；由于种群间可能存在未知倒位，异地个体引入须谨慎，同种群个体引入可作为 demographic rescue 的替代策略；目前可人为打破种子休眠以繁殖同种群个体后 reintroduction 至野外；基因编辑技术可用于修复DNA修复基因中的有害突变，但需仔细的生态和伦理评估。研究结论翻译：许多植物物种的濒危状态通常归因于适合度下降，进而导致自然条件下种群规模的持续缩减。然而，高度重复基因组区域在濒危植物衰退中的作用仍知之甚少。本研究提出了濒危植物秤锤树的无缺口T2T参考基因组，其基因组大小为961.94 Mb，N50为78.05 Mb，预测蛋白编码基因36 446个。利用CENH3 ChIP-seq鉴定了全部12条染色体上的着丝粒区域，并揭示秤锤树着丝粒由多种TE家族组成，其中LTR反转录转座子和LINE元件为主要成分。值得注意的是，1号染色体右臂表现出异常的CENH3富集和类似着丝粒区域的TE积累，暗示形成了着丝粒样染色质结构域。此外，在自然种群中检测到倒位多态性和倒位杂合子，主要集中于同一染色体臂。这些特征可能与减数分裂异常相关。利用高质量T2T基因组组装，鉴定了这些小种群面临的基因组侵蚀模式，其特征为与LTR增殖相关的局部重组抑制以及DNA修复基因中有害突变的积累。这种基因组侵蚀形成了进化困境，阻碍了有害等位基因的清除，最终威胁长期适应性和生存。本研究强调了T2T基因组组装在濒危植物保护中的效用。
