沙棘种质的表型与分子分析揭示与微量营养素相关的优良基因型及候选基因

《BMC Plant Biology》：Phenotypic and molecular analysis of seabuckthorn accessions reveal promising genotypes and candidate genes associated with micronutrients

　　

沙棘（沙棘属Hippophae L.，胡颓子科Elaeagnaceae）富含维生素C、花青素以及铁（Fe）和锌（Zn）等矿物质，是应对微量营养素营养不良的潜在候选物种。本研究中，研究人员在巴基斯坦北部吉尔吉特地区（海拔2444–3172米）的70份沙棘种质中发现了Fe（39.7–341.5 mg kg^?1）、Zn（0.6–12.7 mg kg^?1）和维生素C（40–595 mg 100 g^?1）的高度变异性。低海拔地区大量分布的橙/黄色果实表现出最高的Fe浓度（250.03 ± 61.5 mg kg^?1）。尽管植酸（phytic acid, PA）水平较高（785–1811 mg kg^?1），但PA:Fe < 0.8、PA:Zn < 0.04的低比值表明其具有较高的生物有效性。研究共鉴定到18个HrAKR基因，其中10个通过聚合酶链式反应（PCR）在供试种质中得到验证。3个HrAKR基因呈组成型表达，而7个基因在不同组织中差异表达，提示其参与营养素代谢。高海拔地区（Sost）的种质兼具较高的微量营养素水平、更高的生物有效性以及更强的AKR表达。这些发现凸显了 hrAKR 关联性状通过基因聚合（gene pyramiding）提高沙棘营养价值的潜力。尚需进一步的功能研究以明确每个HrAKR基因的功能，从而为沙棘育种计划开发分子标记。

沙棘作为一种营养密集型的被忽视和未充分利用物种（neglected and underutilized species, NUS），在应对全球微量营养素营养不良方面具有重要潜力。然而，尽管沙棘含有丰富的维生素C、花青素、Fe和Zn等微量营养素，其遗传机制尤其是控制微量营养素积累的分子基础仍缺乏深入研究。当前生物强化工作多集中于主粮作物，而对沙棘等NUS的遗传机制表征不足；同时，虽然已有研究涉及沙棘的成分分析、加工利用及抗氧化特性，但关于营养素生物有效性的遗传基础研究相对有限。此外，Fe和Zn的吸收受复杂基因网络调控，涉及醛酮还原酶（aldo-keto-reductase, AKR）等基因家族，但该家族在沙棘中的功能尚未探索。基于此，研究人员开展了这项研究，旨在整合表型、生化和分子数据，为选育高营养型沙棘品种提供遗传基础。该研究发表在《BMC Plant Biology》上，对利用沙棘改善微量营养素营养安全具有重要意义。

研究人员的主要技术方法包括：以巴基斯坦北部吉尔吉特地区五个不同海拔地点的70份沙棘种质为样本队列，进行土壤理化性质分析（包括pH、电导率、有机质及Fe、Zn等养分含量）；对果实形态性状（果色、果径、果积、含水率、鲜重和干重）进行表型鉴定；采用湿法消化结合原子吸收光谱法测定果实Fe和Zn含量；选取10份代表性种质，用2,6-二氯靛酚（DCPIP）滴定法测定抗坏血酸含量，用试剂盒法测定植酸含量并计算PA:Fe和PA:Zn摩尔比以评估生物有效性；通过Pfam数据库鉴定沙棘基因组中的AKR结构域基因，进行多物种（拟南芥、莱茵衣藻、栽培草莓）系统发育分析、保守基序和结构域分析、亚细胞定位预测及蛋白质理化性质分析；设计基因特异性引物，在70份种质中进行PCR验证HrAKR基因的存在/缺失；选取10份种质的果实和果叶组织，通过半定量逆转录PCR（semi-quantitative RT-PCR）分析HrAKR基因的表达模式；运用方差分析（ANOVA）、巢式设计方差分析、遗传参数估算（遗传力、遗传进展、变异系数）、主成分分析（PCA）和Pearson相关性分析进行数据统计。

研究结果部分包含以下方面：

土壤养分特征：五个地点土壤均为壤质、碱性（pH 7.6–8.4），高海拔地区土壤Fe、P、Zn有效性较低且pH更高，低海拔地区土壤更肥沃。

果实颜色分类与微量营养素浓度：70份种质可分为7种果色类型，橙黄色（31%）和红?（15%）、橙色（14%）为主。Fe浓度范围40–342 mg kg^?1，Zn为0.6–12.7 mg kg^?1。橙黄色果实Fe浓度最高（250 ± 61.5 mg kg^?1），显著高于红色及其 Shades。Shishkat（低海拔）和Gulmit、Misgar种质Fe含量较高，Passu地区Zn变异最大。果色与Fe浓度呈显著正相关。

遗传参数估算：方差分析表明所有性状在种质间差异极显著（p < 0.001）。Fe、Zn、果长、鲜重的表型变异系数（PCV）和基因型变异系数（GCV）均较高；Fe和Zn的广义遗传力（H²）高达97%，遗传进展（GAM）分别为71%和较高值，表明这些性状受加性基因控制，选择改良潜力大。果色、果径、干重等受环境中等影响。

抗坏血酸浓度与微量营养素生物有效性：10份核心种质的植酸含量785–1811 mg kg^?1，抗坏血酸155–649 mg 100 g^?1。所有种质PA:Fe < 0.8、PA:Zn < 0.04，低于生物有效性阈值，表明Fe和Zn生物有效性良好。高海拔种质植酸较高，但S8等种质兼具高抗坏血酸、低植酸和优越的生物有效性。

HrAKR基因的鉴定与系统发育关系：从沙棘基因组鉴定到19个AKR结构域基因，其中18个为非冗余序列（HrAKR1–18）。与拟南芥（22个AtAKR）、栽培草莓（102个FaAKR）和莱茵衣藻（19个CrAKR）的系统发育分析将所有AKR蛋白分为6个主要聚类（AKR-1至AKR-6），沙棘成员分布于所有聚类。

亚细胞定位与蛋白理化性质：HrAKR蛋白预测定位于细胞质（7个）、叶绿体（4个，部分兼有细胞质定位）、细胞核（2个）、线粒体（2个）、过氧化物酶体等。蛋白长度121–451个氨基酸，分子量13.4–50.56 kDa。

结构域、基序与保守残基分析：所有HrAKR蛋白仅含NADPH依赖性氧化还原酶结构域。基序分析显示除基序10外高度保守，HrAKR02仅含3个基序，可能存在功能分化或假基因化。序列标志分析揭示Gly173、Gly179、Val180、Ser181、Asn182等高度保守残基。

HrAKR在种质中的转录丰度：10个HrAKR转录本在70份种质中得到PCR验证。HrAKR08在所有种质中均有表达，HrAKR10和14除个别种质外普遍存在。HrAKR03和15仅少数种质表达。同一区域种质的果叶组织表达模式相似，但HrAKR01在果实中表达降低。S8种质具有最多数目的表达转录本，与其高微量营养素和高生物有效性一致。

讨论部分，研究人员首先强调NUS在应对贫困和营养不良中的潜力，指出本研究为沙棘的生物化学和基因表达提供了综合证据。关于果色与Fe关联，研究确认橙黄色果实与高Fe浓度的统计相关性，这与β-胡萝卜素增强非血红素Fe生物有效性的机制相符，但需分子水平验证以确立果色作为营养表型标记的可靠性。海拔梯度分析显示高海拔（Misgar）和低海拔（Shishkat）种质Fe含量均较高，而中海拔（尤其Sost）较低；高海拔土壤养分有效性虽低，但种质可能通过更强的基因型调控实现养分吸收，这与某些喜马拉雅物种报道不同，体现了物种特异性响应。

微量营养素生物有效性方面，尽管植酸普遍被认为抑制矿质元素吸收，但所有种质的PA:Fe和PA:Zn比值均远低于临界阈值，结合高抗坏血酸水平，表明沙棘果实具有优越的Fe和Zn生物有效性。抗坏血酸通过将Fe³⁺还原为Fe²⁺促进Fe吸收，与拟南芥中的发现一致。遗传分析表明Fe、Zn等性状受加性基因主控，适于直接选择或分子标记辅助选择；而果色、果径等受非加性基因影响，建议采用轮回选择或杂种优势育种。

分子层面，10个HrAKR基因在种质中的存在/缺失变异可能源于结构变异、假基因预测或引物结合位点多态性。所有HrAKR蛋白的保守NAD(P)H依赖氧化还原酶结构域与蒺藜苜蓿（Medicago truncatula）、陆地棉（Gossypium hirsutum）和白菜（Brassica rapa）同源蛋白一致。亚细胞定位预测显示细胞质、叶绿体和细胞核为主要定位位点：叶绿体定位的HrAKR10、14等可能参与光合氧化还原平衡或代谢物解毒；细胞核定位的HrAKR08、18可能调控核内氧化还原过程或基因表达；兼具叶绿体和线粒体定位的HrAKR12可能具有多重功能；含核输出信号的HrAKR13提示可能参与核质运输；过氧化物酶体定位的HrAKR15可能参与脂质代谢或活性氧解毒。系统发育中沙棘与拟南芥同源性较高，但HrAKR01、02、09、16与草莓同源蛋白聚类更近，暗示潜在功能相似性。表达分析显示S8种质具有最丰富的HrAKR表达谱，与其营养优势对应，但需定量PCR进一步验证。

研究结论部分翻译如下：本研究首次对沙棘的Fe、Zn和维生素C进行了整合的形态学、生化学和分子学表征。橙黄色果色与高Fe浓度之间的强相关性表明其作为营养富集种质视觉标记的可能潜力，但该统计学关联的遗传基础尚需进一步研究。海拔变异也影响微量营养素积累，最高和最低海拔的种质相对于中海拔具有更高的营养素含量。在分子水平上，研究鉴定到18个HrAKR基因，其中HrAKR03、13和15在种质间差异表达，但需通过定量PCR或功能研究加以验证。关键HrAKR的保守NADPH结构域以及叶绿体和细胞质定位表明其可能参与维生素C生物合成、Fe还原、氧化还原稳态和养分吸收，但需功能研究验证。数量遗传分析表明Fe、Zn、果长和果重主要受加性基因作用控制，可通过直接表型选择或分子标记辅助选择进行改良，但需窄义遗传力估计验证。果座和果色可能对轮回选择或子代测验响应更佳。该综合分析揭示类胡萝卜素丰富的果色、低植酸含量和高抗坏血酸水平是提高沙棘矿质元素（尤其是Fe）生物有效性的关键标志。这些发现为培育营养密集型沙棘品种、促进可持续营养安全计划建立了实践框架。