《The Plant Genome》：Insights into the Huanglongbing defense responses in the Australasian citrus species
编辑推荐：
柑橘产业因黄龙病（Huanglongbing，HLB）这一细菌性病害而大幅减产，已在全球范围内造成严重经济损失。长期田间研究表明，柑橘属物种对HLB的耐受性存在差异，这凸显了在育种中引入病害耐受性状的重要性。本研究对柑橘种质资源的系统发育关系和遗传多样性进行了
柑橘产业因黄龙病（Huanglongbing，HLB）这一细菌性病害而大幅减产，已在全球范围内造成严重经济损失。长期田间研究表明，柑橘属物种对HLB的耐受性存在差异，这凸显了在育种中引入病害耐受性状的重要性。本研究对柑橘种质资源的系统发育关系和遗传多样性进行了分析，以支持其在育种计划中的有效利用。研究人员选取了代表不同类群的32份柑橘种质资源，采用系统发育分析、同源纯合性分析与主成分分析（principal component analysis，PCA）进行研究。PCA将栽培品种分为两个明显不同的类群，解释了73.3%的总遗传变异。这一聚类揭示了澳大拉西亚柑橘内部独特的遗传多样性，凸显了其在种质保存和未来研究中的价值。研究人员将被归类为HLB抗性或HLB耐受性的澳大拉西亚种质资源，与感病品种甜橙（Citrus sinensis L. Osbeck，品种Hamlin 1-4-1）进行比较，利用诊断变异（diagnostic variant，DV）鉴定潜在的防御机制。对这些DV中的转座子（transposable element，TE）注释发现，转座子插入发生在基因内部，提示其可能影响防御反应过程中的基因调控。基于抗性组与耐受组之间共享序列比对的功能注释发现，一个仅在感病品种中出现的氨基酸替换可能损害HLB防御功能；而澳大拉西亚柑橘物种中未发生此类替换。受影响的基因与植物防御的关键组分相关，包括信号通路、植物激素调控和次生代谢物生物合成。总体而言，本研究凸显了柑橘种质资源在对抗HLB方面具有宝贵的遗传多样性。
黄龙病（Huanglongbing，HLB）是严重威胁全球柑橘产业的一种细菌性病害，导致美国佛罗里达州近十年柑橘产量下降约90%。商业化栽培品种（如甜橙）遗传背景狭窄，主要源于珠心多胚性和无性繁殖，限制了通过传统育种改良抗病性的能力。澳大拉西亚柑橘物种（原产于澳大利亚、新几内亚岛及附近太平洋岛屿）因长期地理隔离，可能演化出独特的HLB抗性或耐受性机制，是拓宽栽培柑橘遗传基础的宝贵种质资源。然而，目前对澳大拉西亚柑橘与栽培品种间的进化关系及其HLB防御机制的认识仍存在显著不足。为此，研究人员对32份柑橘种质资源开展了群体基因组学研究，通过系统发育分析、同源纯合性（homozygosity by descent，HBD）分析和主成分分析（principal component analysis，PCA）阐明其遗传关系，并利用诊断变异（diagnostic variant，DV）比较抗性/耐受性澳大拉西亚种质与感病甜橙之间的基因差异，以期揭示潜在的HLB防御机制并为未来育种策略提供依据。
研究人员利用Illumina平台对32份柑橘材料进行全基因组测序，其中13份材料的测序数据来自NCBI SRA数据库，其余材料来自美国国家柑橘无性系种质资源库（USDA）、佛罗里达大学（UF）和佛罗里达州植物工业部（DPI）的叶片样本。以C. australis基因组v1.0（UQ）为参考基因组，经变异检测和质量筛选获得116,943个高质量双等位基因SNP。系统发育树采用最大似然法（RAxML）构建，PCA通过PLINK和LEA软件包完成。HLB抗性/耐受性组与感病品种C. sinensis（Hamlin 1-4-1）的DV比较利用SnpEff进行注释，转座子注释采用EDTA流程，功能富集通过KEGG和STRING完成。
结果方面，3.1 Identification of SNP variants and heterozygosity：变异检测共获得116,943个高质量SNP，分布于9条染色体上，其中5号染色体最多（21,623个），9号染色体最少（7,696个）。多态信息含量分析显示PIC为0.37时位点数量最多（103,981个）。所有材料平均观测杂合度为0.07，其中C. jambhiri最高（0.288），C. polyandra最低（0.001），澳大拉西亚物种整体杂合度较低，表明其遗传背景较为纯合。
3.2 Phylogenetic analysis：系统发育分析将32份材料分为两个大类群。第一大类群包含栽培品种和祖先类型，其中Clade I包括宽皮橘-柚复合群及其杂交后代（甜橙、葡萄柚、酸橙），第二簇包含柠檬、来檬和香橼类。第二大类群为澳大拉西亚野生种，形成单系类群（100% bootstrap支持），其中C. polyandra位于最基部，其余物种分为两个姐妹亚分支：一个亚分支包含C. glauca与C. garrowayae，另一较大亚分支按序包含C. papuana、C. warburgiana以及由C. australis复合群与C. inodora、C. australasica构成的核心类群。
3.3 Homozygosity by descent analysis：HBD分析显示不同材料具有独特的HBD谱型。澳大拉西亚物种具有较长HBD片段，其中C. australasica（PI 312872）在R=2类中比例最高，提示近期发生高强度近交。栽培品种则表现较短的纯合片段和较高的R类比例（R≥1024），反映其历史上经历遗传混合与种间杂交。
3.4 Principal component analysis：PCA第一主成分解释73.3%的遗传变异，将栽培品种与澳大拉西亚物种明显分开。进一步的分组PCA显示，澳大拉西亚组内C. polyandra（CREC）与来自澳大利亚的物种分离，反映其独特的新几内亚岛遗传背景；非澳大拉西亚组则形成六个聚类，包括宽皮橘、甜橙、葡萄柚、柚、来檬/柠檬以及枳壳类。
3.5 Functional annotation of the diagnostic variants in the Australasian citrus species：DV注释鉴定出抗性组11,843个和耐受组11,300个高/中影响变异。TE分析发现抗性组68个、耐受组65个与DV重叠的TE插入位点，包含Class I反转录转座子（LTR、gypsy、Copia）和Class II DNA转座子（Mutator、CACTA、PIF Harbinger、Tc1 Mariner），共定位到47个共享基因。非同义变异筛选获得抗性组1,448个和耐受组1,427个基因，经严重度评分后，抗性组有202个高影响和186个中影响基因，耐受组有198个高影响和181个中影响基因，其中311个基因为两组共享。功能注释表明这些基因涉及PTI/ETI信号、细胞壁加固、激素调控和次生代谢等防御通路，感病甜橙中存在氨基酸替换而澳大拉西亚物种中这些基因序列保持完整。
3.6 Identifying the KEGG pathways and gene clusters linked to HLB defense responses in the Australasian citrus species：KEGG通路富集显示MAPK信号通路、α-亚麻酸代谢、植物激素信号转导和苯丙烷生物合成显著富集。STRING聚类分析鉴定出60个基因簇，主要包括37个基因参与的淀粉和蔗糖代谢与核苷酸糖生物合成，23个基因的RNA介导基因沉默和胞嘧啶C-5甲基化，22个基因的苯丙烷生物合成与类黄酮生物合成。
3.7 Functional annotation of the unique DV-associated genes in the resistant and tolerant groups：抗性特异基因主要涉及早期病原识别（Ca2+信号、LRR受体激酶、离子通道）和程序性细胞死亡（metacaspase 4），同时包含细胞壁重塑和抗菌次生代谢物合成基因。耐受特异基因则侧重激素平衡（赤霉素分解、茉莉酸合成、生长素信号）、细胞壁完整性维持、代谢解毒、物质转运和蛋白稳态（chaperonin-60α），在抵御病害的同时维持植物正常生长。
讨论部分，研究人员指出澳大拉西亚柑橘物种与栽培品种的遗传分离源于长期的地理隔离，低杂合度和长HBD片段反映其遗传纯合与近交历史。TE插入可能改变防御相关基因的应激调控网络，而感病品种中多种防御基因的氨基酸替换可能损害PTI/ETI信号、钙-ROS交互和激素串扰，削弱其抵抗HLB的能力。KEGG分析揭示的淀粉和蔗糖代谢基因在细胞壁维持和养分转运中发挥关键作用，其完整保留于澳大拉西亚物种中有助于抵抗病原入侵。结论部分，本研究通过全基因组分析明确了澳大拉西亚柑橘物种的独特遗传地位，鉴定了一系列与HLB防御相关的功能基因，为利用这些野生种质培育HLB抗性新品种提供了重要的分子标记和理论基础。未来研究将致力于开发具有多样化遗传背景的HLB抗性杂交品种，并利用分子标记辅助早期筛选以加速育种进程。
