基于非靶向代谢组学与气相色谱–离子迁移谱（GC–IMS）解析武夷与白沙产区奇兰鲜叶品质相关差异

《Foods》：Untargeted Metabolomics and GC–IMS Reveal Quality-Related Differences in Qilan Fresh Tea Leaves Grown in Wuyi and Baisha

　　

在这项研究中，研究人员系统比较了同一 Camellia sinensis (L.) O. Kuntze cv. Qilan 品种鲜叶，其分别种植于海拔相似但纬度差异明显的福建武夷山（2L-FQ栽培点）和海南白沙（2L-HQ栽培点），并采用环境表征、酶活性测定、气相色谱–离子迁移谱（gas chromatography–ion mobility spectrometry，GC–IMS）以及非靶向代谢组学（untargeted metabolomics）进行分析。两组鲜叶在相同标准化加工条件下制成的红茶，进一步通过电子鼻（electronic nose）和电子舌（electronic tongue）分析进行评价。相应地，2L-HQ叶片比2L-FQ叶片更长、更宽，水分含量和水活度更高。2L-HQ较低的L*、b*和色度（chroma）值以及ΔE*_ab为6.54，表明其叶色更深，并伴随更高的叶绿素和类胡萝卜素浓度。与2L-FQ相比，2L-HQ茶多酚含量高11.8%，可溶性多糖高26.4%，咖啡因高11.4%，游离氨基酸高10.6%。非没食子酰化儿茶素（non-gallated catechins）中的儿茶素（catechin）、表儿茶素（epicatechin）、表没食子儿茶素（epigallocatechin）和没食子儿茶素（gallocatechin）在2L-HQ中显著富集。2L-HQ的DPPH（2,2-diphenyl-1-picrylhydrazyl）、ABTS（2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid)）和FRAP（ferric-reducing antioxidant power）抗氧化能力分别高6.1%、4.9%和2.2%。脂氧合酶（lipoxygenase，LOX）、多酚氧化酶（polyphenol oxidase，PPO）和过氧化物酶（peroxidase，POD）活性分别增加137.0%、98.0%和93.8%，而β-葡萄糖苷酶（β-glucosidase）活性保持不变。GC–IMS解析出61种挥发性化合物，其中46种在两个产地间存在显著差异。2L-FQ富含若干芳香醛、β-苯乙醇（β-phenyl ethanol）和果味酯类，而2L-HQ含有较高水平的选定醇类、不饱和醛、呋喃类化合物和脂质衍生酯类。非靶向代谢组学保留4359个代谢特征，其中4356个为两组共有。多重检验校正后，采用VIP > 1、FDR校正p < 0.05和|log₂FC| > 1鉴定出793个差异代谢特征，其中相对2L-FQ，2L-HQ中420个特征丰度更高、373个特征丰度更低。这些差异特征主要与碳水化合物、氨基酸、类黄酮相关、脂质和辅因子代谢相关。总体而言，区域栽培与叶片形态、水合、色素沉着、滋味相关组成、抗氧化能力、氧化酶潜力、挥发性指纹和代谢分配的协调差异相关，为奇兰茶叶区域适配加工提供了数据基础。

研究背景方面，茶是全世界广泛消费的非酒精饮料，鲜叶是茶叶加工的主要原料，其理化性质和代谢谱构成最终茶叶品质的化学基础。儿茶素、咖啡因、氨基酸、有机酸、糖和类黄酮等非挥发性化合物参与苦味、涩味、鲜味、甜味和回甘形成；来自脂质代谢、类胡萝卜素降解、糖苷结合前体和美拉德反应相关途径的挥发性化合物塑造茶叶香气。因此，解析鲜叶理化性质和代谢物对评价加工适宜性和风味形成能力十分重要。除品种和加工条件外，生长环境进一步塑造茶叶理化与代谢特征。纬度、海拔等地理因素可影响温度、相对湿度、水分有效性和光照，并与土壤性质共同影响茶次生代谢及多酚、氨基酸、糖和香气相关前体水平。Tian et al. 用代谢组学方法表明，栽培海拔影响茶叶代谢谱，尤其是类黄酮和氨基酸代谢相关代谢物。因此，比较同一品种在对比区域环境下的表现，可减少遗传变异并评估环境相关鲜叶品质变化。奇兰（Camellia sinensis (L.) O. Kuntze cv. Qilan）是乌龙茶品种，具有形成花香和果香的潜力。以往研究多集中于加工茶或加工诱导的香气变化，对跨区域栽培环境如何重塑奇兰鲜叶理化属性、挥发性指纹和非挥发性代谢谱的系统信息仍有限。GC–IMS适合快速挥发性指纹分析和样品判别，非靶向代谢组学可全面分析整体代谢变化；与酶活性测定和靶向品质分析结合，可将挥发性指纹与上游代谢变化及鲜叶品质性状联系起来。因此，研究人员比较福建武夷与海南白沙的气温、相对湿度、降水、太阳辐射和风速，以及奇兰鲜叶品质性状、酶活性、挥发性指纹和代谢谱，旨在识别与滋味相关化合物和香气前体相关的关键代谢变化，为奇兰茶叶区域特异加工和利用提供科学依据。

关键技术方法方面，样本来源于海南农垦白沙茶业有限公司（Hainan State Farms Baisha Tea Co., Ltd.，白沙黎族自治县）和Mount Wuyi Yiqingerming Tea Industry Co., Ltd.（武夷山市）。研究人员采集同一奇兰品种、13年树龄、一芽二叶鲜叶，采用五点采样，每点五株、每株四梢，混合为生物重复；两个产地分别记为2L-FQ和2L-HQ。环境数据来自NASA POWER。关键技术包括环境表征；叶长、叶宽、叶厚、色差、水分含量、水活度（water activity，a_w）、叶绿素和类胡萝卜素测定；总多酚、可溶性多糖、游离氨基酸及HPLC（high-performance liquid chromatography，高效液相色谱）分析没食子酸、咖啡因和儿茶素；DPPH、ABTS和FRAP抗氧化测定；LOX、PPO、POD和β-葡萄糖苷酶活性测定；HS–GC–IMS挥发性指纹；UPLC–MS/MS（ultra-performance liquid chromatography–tandem mass spectrometry，超高效液相色谱–串联质谱）非靶向代谢组学；标准化加工红茶后电子鼻和电子舌分析。统计采用t检验、PCA（principal component analysis，主成分分析）、OPLS-DA（orthogonal partial least squares-discriminant analysis，正交偏最小二乘判别分析）、200次置换检验、层次聚类、火山图和KEGG（Kyoto Encyclopedia of Genes and Genomes，京都基因与基因组百科全书）富集。

研究结果如下。3.1 环境背景与鲜叶表型性状：通过环境数据、表型和色素测定，两个样点海拔相似，主要差异为水热和辐射条件；2L-HQ栽培点更暖、更湿，平均气温、相对湿度和总降水更高，太阳辐射低于2L-FQ，表层和根区土壤湿度略高。2L-HQ叶片更长、更宽，叶面积更大，水分含量和水活度更高；其L*、b*和c*值较低，ΔE*_ab为6.54，叶色更深，叶绿素a、叶绿素b、总叶绿素和类胡萝卜素含量更高。结论是区域环境相关鲜叶更大、更水合且色素积累更高。3.2 主要品质相关成分：通过靶向成分和抗氧化测定，2L-HQ茶多酚为219.67 ± 2.46 mg g^?1 DW，较2L-FQ的196.42 ± 1.75 mg g^?1 DW高11.8%；可溶性多糖、咖啡因和游离氨基酸更高；C、EC、EGC和GC显著富集，CG较低，ECG、GCG和EGCG无显著差异；DPPH、ABTS和FRAP抗氧化能力较高。结论是2L-HQ具有更高酚类浓度和不同儿茶素谱，不能直接等同于涩味增强，其滋味基础更复杂。3.3 酶活性谱：通过酶活性测定，LOX、PPO和POD活性分别增加137.0%、98.0%和93.8%，β-葡萄糖苷酶活性无显著变化。结论是区域差异集中在脂质氧化和酚氧化，而非糖苷结合香气前体水解。3.4 GC–IMS挥发性指纹：通过HS–GC–IMS，共暂定鉴定61种挥发性化合物，其中46种显著不同，35种在2L-FQ中更高，11种在2L-HQ中更高。2L-FQ富含2-己烯醛、苯乙醛、壬醛、(Z)-6-壬烯醛、香茅醛、β-苯乙醇、辛酸甲酯、丁酸丁酯、2-甲基丙基丁酸酯、茶螺烷、吲哚、2-辛酮、2-戊基呋喃、香芹酚和3-甲基环戊烷-1,2-二酮等；2L-HQ富含2,3-丁二醇、1-戊烯-3-醇、二氢芳樟醇、(E,E)-2,4-庚二烯醛、(E,E)-2,4-壬二烯醛、(Z)-3-己烯基异戊酸酯、乙酸冰片酯、4-甲基苯乙酮和4-羟基-5-乙基-2-甲基-3(2H)-呋喃酮等。PCA中PC1、PC2和PC3分别解释74.9%、13.9%和6.8%，累计95.6%。结论是两地共享挥发本底，区域主要改变相对丰度。3.5 非靶向代谢组学揭示整体代谢差异：通过UPLC–MS/MS，保留4359个代谢特征，其中4356个共有，仅2个为2L-HQ独有、1个为2L-FQ独有。注释代谢物中，有机酸及衍生物占26.21%，类黄酮占13.13%，苯环型占12.55%，脂质及类脂分子占7.32%，氨基酸及衍生物占3.72%，核苷酸及衍生物占3.72%，苯丙素和聚酮占3.60%。PCA中PC1和PC2分别解释53.70%和13.64%，累计67.34%；OPLS-DA预测成分解释55.6%，第一正交成分解释11.9%，原模型R²X = 0.674、Q² = 0.965，200次置换中0/200超过原模型，经验p < 0.005。差异特征793个，其中420个在2L-HQ中更高、373个更低。VIP排名中暂定标注包括L-苯丙氨酸、吲哚-3-甲醛、3-吲哚丙烯酸、3,7-二甲基-1,6-辛二烯-3-醇和methyl 7-epi-12-hydroxyjasmonate glucoside等。结论是区域差异主要是共享代谢特征的定量重编程。3.6 差异代谢通路富集：通过KEGG富集，差异代谢特征主要关联碳水化合物代谢、氨基酸代谢、类黄酮相关生物合成、氧化还原调节和辅因子代谢。富集因子最高为淀粉和蔗糖代谢，其次为硫辛酸代谢、膦酸和次膦酸代谢、半乳糖代谢和丙酮酸代谢；还涉及酪氨酸代谢、缬氨酸亮氨酸异亮氨酸生物合成、维生素B6代谢、辅因子生物合成、异黄酮和花青素生物合成、抗坏血酸和醛酸代谢等。结论是区域分化涉及碳代谢、氨基酸代谢、类黄酮相关生物合成和氧化还原/辅因子代谢的协调变化，但富集分析不指示通路激活或抑制。3.7 电子鼻与电子舌：通过标准化加工红茶电子鼻和电子舌分析，电子鼻雷达图整体响应模式相似，但个别传感器响应强度不同，尤其W1W在2L-HQ中更高，电子鼻PCA可区分成品茶；电子舌雷达图中2L-HQ略高，电子舌PCA中PC1和PC2分别解释57.3%和28.3%，累计85.6%，2L-FQ位于PC1负侧，2L-HQ位于正侧。结论是加工后仪器响应差异仍可检测。

讨论与结论部分，原文结论可翻译为：本研究在海拔、品种、树龄和采摘标准可比的条件下，系统表征了福建武夷山市和海南白沙黎族自治县栽培奇兰鲜叶的区域相关差异。与2L-FQ相比，2L-HQ在更暖、更湿和较低辐射环境下发育，叶片更大、更水合、颜色更深、色素含量更高。白沙鲜叶还含有更高水平茶多酚、可溶性多糖、咖啡因、游离氨基酸和若干非没食子酰化儿茶素，并具有更高体外抗氧化能力。LOX、PPO和POD活性显著增加，表明更强的脂质和酚氧化潜力，而β-葡萄糖苷酶活性不变。GC–IMS检测到61种挥发性化合物，其中46种在两地间显著不同。2L-FQ以若干芳香醛、β-苯乙醇和果味酯类较高为特征，而2L-HQ积累较高水平选定醇类、不饱和醛、呋喃类化合物和脂质衍生酯类。非靶向代谢组学保留4359个代谢特征，4356个共有；793个差异特征产生清晰代谢分离，主要关联碳水化合物、氨基酸、类黄酮相关、脂质、氧化还原和辅因子代谢。这些发现表明，奇兰鲜叶区域分化主要来自一个大体共享化学库的定量重分配，而非形成大量产地特异化合物。电子鼻和电子舌结果进一步提示成品具有可区分的香气和滋味响应指纹。总体而言，对比环境与叶片生理、滋味相关成分、氧化酶能力、挥发组成和代谢分配的协调变化相关，为区域适配加工策略提供基础。然而，研究仅涉及两个栽培点和一个采样时期，多季节和多地点研究不能分离单个环境驱动因子、验证通路激活或建立优越感官品质。需要在多个地点和季节重复、靶向化学确认、控制加工比较和标准化感官评价，以评估这些发现的普遍性和实际意义。结合关键代谢物靶向验证和成品茶标准化感官评价，才能确定这些差异的环境驱动因素和实际品质后果。