《Frontiers in Plant Science》：Effect of phosphorus deficiency on biomass and root system architecture in diverse Medicago accessions
编辑推荐：
为应对磷(P)匮乏对豆科作物生长和固氮的制约，研究人员对200份紫花苜蓿复合体种质开展了表型和遗传学研究。他们评估了高低P条件下种质的生长差异，通过k-mer GWAS鉴定了与株高、生物量相关的候选基因，并在20份代表性种质中深入分析了根系构型(RSA)可塑性及磷利用效率(PUE)。结果表明，存在显著的种内变异，低P条件下部分种质可维持稳定的生物量、较高的PUE及适应性RSA特征，为选育磷高效苜蓿品种提供了宝贵的遗传资源和理论基础。
磷，如同植物世界的“能量货币”，对生命活动至关重要，尤其对于依赖生物固氮的豆科植物。紫花苜蓿作为一种高产、高蛋白的多年生豆科牧草，其生长和固氮过程对磷的供应极为敏感。然而，全球许多土壤中磷的有效性很低，成为限制植物生产力的关键因素。更严峻的是，磷肥的利用率普遍偏低，大量未被吸收的磷不仅造成资源浪费，还可能随径流流失，引发水体富营养化。面对磷矿资源有限和农业可持续发展的双重挑战，培育能够高效获取和利用磷的作物品种显得尤为迫切。植物如何应对磷匮乏？不同的紫花苜蓿种质资源在磷胁迫下会表现出怎样的生长差异和适应策略？其背后的遗传基础又是什么？为了回答这些问题，一项题为《Effect of phosphorus deficiency on biomass and root system architecture in diverse Medicago accessions》的研究在《Frontiers in Plant Science》上发布，系统探究了紫花苜蓿复合体种质在磷胁迫下的表型与遗传变异。
研究者们开展了一项整合了温室表型鉴定、基因组关联分析和根系精细剖析的系统性工作。他们首先从一个包含1,234份种质的大型基因库中，基于基因分型测序(GBS)数据筛选出200份遗传背景多样的紫花苜蓿复合体种质，涵盖了栽培种、地方品种和野生材料。在温室实验中，这些种质被置于高磷(HP)和低磷(LP)两种条件下培养，以评估其地上部和地下部生物量、根冠比等性状的差异。为了挖掘控制这些性状的遗传位点，研究采用了不依赖于参考基因组的k-mer全基因组关联分析(k-mer GWAS)方法，直接对GBS原始测序reads中的短序列片段(k-mer)进行关联分析，以捕捉更广泛的遗传变异，包括结构变异。随后，他们根据根干物质对磷供应的响应差异(ΔRDM)，从200份种质中挑选出20份代表性材料（10份高ΔRDM，10份低ΔRDM），进一步测定了其地上部磷含量和磷利用效率(PUE, 即单位磷所生产的地上部干物质量)。最后，利用根窗系统对这些代表性种质的根系构型进行了精细成像和量化分析，获取了总根长、比根长、根系空间分布等十余个关键性状参数。
3.1 盆栽实验中植物在磷缺乏下的表现
研究发现了显著的种内变异。在200份种质中，低磷处理平均导致地上部干物质和根干物质分别下降24%和23%。然而，部分种质表现出了对低磷胁迫的较强耐受性。例如，种质LE 2701在高低磷条件下的地上部和根系干物质均无显著差异。通过计算胁迫耐受指数(STI)，研究者鉴定出了如LE 2673, LE 816, LE 808等在高磷和低磷条件下均能保持较高生物量生产的优异种质。这些结果通过图1的箱线图得到了直观展示，该图清晰对比了高低磷条件下地上部干物质、根干物质和根冠比的分布差异。
3.2 基于k-mer GWAS信号鉴定候选基因
k-mer GWAS分析成功鉴定出与多个性状相关联的基因组区域。例如，在高磷条件下的最大株高关联区域中，发现了编码26S蛋白酶体AAA-ATP酶亚基、光敏色素相互作用bHLH转录因子和微管正端追踪蛋白的候选基因，这些基因与分生组织活性和细胞扩展相关。而在低磷条件下关联到的基因，则更多地与激素信号（如生长素生物合成的YUCCA酶、细胞分裂素激活的LOG酶、茉莉酸代谢的JMT酶）、钙信号（类钙调磷酸酶B钙传感器）以及细胞壁修饰（GT14家族糖基转移酶）等胁迫响应和代谢调整通路有关。这些候选基因列表详见表1，为理解紫花苜蓿响应磷胁迫的分子机制提供了重要线索。
3.3 地上部磷含量与磷利用效率
对20份代表性种质的分析显示，低磷处理显著降低了所有种质的地上部磷含量，但磷利用效率(PUE)却平均增加了38%。这表明植物通过提高内部磷的利用效率来部分补偿吸收的不足。值得注意的是，PUE的变化与生物量的变化并无显著相关性，说明磷吸收（或获取）和内部利用是两个相对独立的过程。高ΔRDM和低ΔRDM种质组在PUE上并无组间差异，但在磷含量对处理的响应上存在交互作用。这些关系通过图3的柱状图得以呈现。
3.4 根窗实验中植物在磷缺乏下的表现
根窗实验证实了低磷对生物量的抑制作用，同时也揭示了根系构型(RSA)的可塑性变化。尽管总根长、根体积、根表面积等反映根系大小的性状在低磷下普遍降低，但反映根系觅食效率的性状，如比根长和比凸包面积，却显著增加。这意味着植物在资源受限时，倾向于投资构建更细、空间探索效率更高的根系。此外，根系在水平方向的扩展（根宽、凸包宽）在低磷下有所增加，而扎根深度未变，这种“浅层化、宽展化”的调整有利于获取富集于土壤表层的难移动磷。不同种质的RSA响应模式存在差异，通过热图聚类分析可将其分为不同的响应类型（图5），展现了遗传背景对根系可塑性的影响。
4 讨论
该研究系统地揭示了紫花苜蓿种质资源应对磷胁迫的多样性与复杂性。在表型层面，研究证实了种内巨大的变异潜力，并鉴定出了一批在低磷下能保持相对稳定生物量、较高PUE和适应性RSA特征的优异种质，例如LE 2673, LE 816, LE 808和LE 2701。这些种质来源于栽培种、地方品种和野生材料，表明有益的磷高效性状存在于广泛的遗传背景中，为育种提供了丰富的素材。在遗传层面，k-mer GWAS首次在该物种中应用于磷响应性状的研究，鉴定出的候选基因主要富集于生长调控、激素信号、胁迫响应和细胞壁代谢等通路，这为解析紫花苜蓿磷效率的分子调控网络奠定了基础。在生理机制层面，研究发现磷吸收（反映在地上部磷含量）和内部磷利用效率(PUE)是相对独立的，植物在低磷下通过提高PUE进行补偿，而根系则通过调整构型（增加比根长、水平扩展）来提升土壤探索效率。这些适应策略共同构成了紫花苜蓿应对磷胁迫的多重防线。
该研究的意义在于，它不仅全面评估了一个大型、多样化的紫花苜蓿种质库的磷效率，将表型变异与遗传关联相结合，还深入剖析了根系构型可塑性这一关键适应性状。研究所采用的k-mer GWAS方法特别适用于紫花苜蓿这种异交、同源四倍体、基因组复杂的物种，能够有效捕捉传统SNP分析可能遗漏的变异。最终，这项研究为后续的功能基因验证、分子标记开发以及旨在培育适合低投入农业系统的磷高效紫花苜蓿品种的育种计划，提供了坚实的表型数据、遗传信息和理论框架。当然，研究也指出了在温室早期生长阶段获得的结果需要在田间和更长生育期进行验证，并且对候选基因的功能解析将是下一步研究的重点。
