巴戟天Morinda officinalis别名巴戟、三蔓草、鸡肠风，是双子叶植物纲Dicotyledoneae、茜草目Rubiales、茜草科Rubiaceae多年生攀援木质藤本植物，是我国“四大南药”之一，以肉质根入药，其主要药效成分为寡糖、多糖蒽醌类、环烯醚萜类等活性物质^[1-2]，具有补肝肾、强筋骨、祛风湿、抗衰老、抗肿瘤和治疗糖尿病、精神疾病以及细菌和病毒感染等作用^[3-7]。巴戟天种植成本较高，且种植过程中存在开垦山林及乱采滥挖现象，这严重破坏了巴戟天的生长环境^[8]。虽然巴戟天是第1批受到《广东省岭南中药材保护条例》保护的岭南药材，但其种植面积在逐年减少。土壤养分是影响巴戟天生长的重要环境因素之一，肉桂Cinnamomum cassia成林速度快，林内土壤动物、微生物种类更为丰富，可有效改善林地内土壤结构，从而提高土壤的肥力，营造出适宜巴戟天生长的种植环境^[9]。

近年来，一些学者从巴戟天的药材化学成分及代谢组学^{[3, 10-12]}、种植条件^[13-14]、种植模式^[15-16]等方面对巴戟天的有效成分含量进行比较研究，发现合理增加肉桂林种植地的光照量，可抑制藤苗生长，从而提高巴戟天肉质根的产量和有效成分积累，有效成分含量均随种植时间的延长呈上升趋势^[16]。林下种植方式是控制种植地光照条件的重要途径之一，而且巴戟天林下种植具有较好的经济效益，解决了坡地连片种植带来的潜在山地环境破坏和水土流失问题，并且提高了巴戟天的产量和药效品质^{[8, 15]}。这一林下生态种植方式为区域中药材的产业化发展提供了新思路。

广东省肇庆市德庆县是巴戟天的主产区，种植巴戟天有着悠久历史，是目前我国种植面积最大、产量最多和品质最好的巴戟天生产基地。本文研究了道地产区土壤养分和生长年限与巴戟天品质的关系，以及林下和非林下种植方式对巴戟天有效成分的影响，以期为巴戟天林下规范化种植提供基础数据和理论支持。

1 材料与方法 1.1 研究区概况

本研究试验地位于广东省肇庆市德庆县，调查采样点分布于N 23°13′30″~N 23°17′37″，E 111°53′52″~E111°55′45″，海拔62~166 m。地处亚热带地区，热量丰富、雨热同期、夏长冬短，年平均温度20~25 ℃，年降雨量600~2 000 mm，土壤以红壤为主，其次是黄壤。德庆县当地农户种植的巴戟天有2个品种，分别是小叶品种‘黑幼籽’与大叶品种‘广宁特’，以小叶品种‘黑幼籽’最为常见。本次研究调查的品种均为‘黑幼籽’，种植密度普遍为株行距20 cm×30 cm，农家肥与大田复合肥为巴戟天施肥最常用肥料。大田复合肥的施肥量较低，通常为224~230 kg·hm⁻²，施肥频率为每年2~3次。

1.2 样品采集

于2018年10月，选取广东省德庆县高良镇以及邻近镇作为研究区域，根据该区域巴戟天种植面积和分布，选取18个具有代表性的非林下巴戟天种植点(表1)，调查分析土壤养分和种植年限对巴戟天寡糖含量的影响。另外，为了调查分析林下和非林下种植巴戟天的品质差异，在其中8个非林下巴戟天种植点的相邻林地上(相关信息见表2)，采集了林下种植点样品，巴戟天品种和种植年限以及土壤、坡度、坡向、海拔等均与对应的非林下种植点一致。在选定的地点，采用随机多点(5~8个点)采样法，在根系外围采集0~20 cm深的土壤，将土壤混匀并用四分法采集分析样品。土壤样品自然风干后，去除植物根系、落叶及石块等杂质，研磨，分别过1.00和0.25 mm筛，保存备用。同时在采样点随机整株挖取5株巴戟天。巴戟天肉质根用自来水清洗3次、去离子水润洗1次后，擦干。去除肉质根的木心部分，烘干，粉碎，过0.25 mm筛，待用。

1.3 样品分析 1.3.1 土壤养分含量分析

土壤样品的养分含量分析按照《土壤农化分析》^[17]进行，具体为：称取10 g过1 mm筛的土壤，加入25 mL无二氧化碳去离子水浸提，利用pH计测定土壤pH；称取1 g过0.25 mm筛的土壤，采用凯氏定氮法测定土壤全氮含量；称取2 g过1 mm筛的土壤，采用碱解扩散法测定土壤碱解氮含量；称取1 g过0.25 mm筛的土壤，采用高氯酸−氢氟酸消煮−钼锑抗比色法测定土壤全磷含量；称取2 g过1 mm筛的土壤，采用氟化铵−盐酸浸提−钼锑抗比色法测定速效磷含量；称取0.25 g过0.25 mm筛的土壤，采用高氯酸−氢氟酸消煮−火焰光度法测定土壤全钾含量；称取5 g过1 mm筛的土壤，采用乙酸铵浸提–火焰光度法测定速效钾含量；称取0.5 g过0.25 mm筛的土壤，采用重铬酸钾氧化外加热法测定土壤有机质含量。

1.3.2 巴戟天肉质根品质成分分析

巴戟天肉质根全氮含量采用H₂SO₄−H₂O₂消煮−奈氏试剂比色法测定；全磷含量采用 H₂SO₄–H₂O₂ 消煮−钼锑抗比色法测定；全钾采用H₂SO₄−H₂O₂消煮−原子吸收法测定。肉质根寡糖含量参照李倩^[18]采用高效液相色谱–蒸发光散射法测定，通过峰面积计算样品中蔗糖、1−蔗果三糖、耐斯糖以及1^F−果呋喃糖基耐斯糖的含量。

1.4 数据分析方法

试验数据采用Excel 2017和SPSS 25.0软件进行整理和统计分析，采用单因素方差分析(One-way ANOVA)，用Duncan’s多重极差检验法比较巴戟天品质指标的差异性；土壤养分与巴戟天品质指标的相关性则采用Pearson法进行分析；林下和非林下品质指标之间的差异性采用配对样本t检验法进行比较。

在采用主成分分析法进行综合评价时，选取反映巴戟天品质综合评估值的n个指标为X₁，X₂，…，X_n，m 个生长年限和种植环境的n项指标构成原始数据矩阵X=[X_ij]_n×m，其中，X_ij为第i个模式的第j项指标数据(i=1, 2, …，n；j=1, 2, …，m)，按公式(1)对原始数据进行标准化，从而消除量纲和数量级的影响，然后确定主成分，从总方差分析表选取累积贡献率≥ 85%的前p个主成分，按公式(2)在建立p个主成分和标准化变量的关系后，按公式(3)确定权重，用第k个主成分的贡献率与选取的p个主成分的总贡献率的比值确定每个主成分的权重，最后按公式(4)构造综合评价函数，根据前p个主成分及其权重构造综合评价函数^[19]。

${X_{ij}}^* = ({X_{ij}} - \overline {{X_j}} )/{S_j},$

(1)

式中，X_ij^*是 X_ij的标准化数据， $\overline {{X_j}} $ 和S_j是第j个指标的平均值和标准差。

${Y_k} = {u_{{k_1}}}{{X}}_1^* + {u_{{k_2}}}{{X}}_2^* + \cdots + {u_{{k_m}}}{{X}}_m^*,$

(2)

式中，Y_k是第k个主成分(k=1, 2, ···，p)，u_k1~u_km分别为第k₁~k_m个主成分的因子荷载，X₁^*~X_m^*均为标准化后的样本矩阵。

${\omega _k} = {\lambda _k}/\sum\limits_{k = 1}^p {{\lambda _k}} ,$

(3)

式中，ω_k为第k个主成分的权重，λ_k为第k个主成分的贡献率。

$F = \sum\limits_{k = 1}^p {{\omega _k}{Y_k}} ,$

(4)

式中，F为不同生长年限和种植环境的巴戟天品质综合评价得分，分值越高表明该条件下的巴戟天品质越好。

2 结果与分析 2.1 道地产区土壤养分含量

广东省德庆县是巴戟天道地产区。由表3可知，所收集到的土壤样品pH在4.08~4.90之间，差异较小，说明巴戟天适合种植在酸性土壤中。土壤有机质质量分数大部分在26~46 g·kg⁻¹之间，11号样品有机质质量分数达到75.15 g·kg⁻¹。土壤全氮、全磷、全钾的含量差异较大，这可能与当地的自然环境、施肥情况密切相关。

2.2 巴戟天品质成分含量

巴戟天药材中全氮质量分数为2.81~12.30 g·kg⁻¹，全磷质量分数为0.12~0.71 g·kg⁻¹，全钾质量分数为6.13~13.20 g·kg⁻¹；巴戟天药材中4种寡糖成分含量差异较大，蔗糖质量分数为34.06~76.16 g·kg⁻¹，1−蔗果三糖质量分数为9.99~59.30 g·kg⁻¹，耐斯糖质量分数为27.93~69.41 g·kg⁻¹，符合2015年版《中华人民共和国药典》巴戟天干燥品耐斯糖含量不得少于2.0%的要求^[20]；1^F−果呋喃糖基耐斯糖质量分数为46.60~97.12 g·kg⁻¹(表4)。

2.3 巴戟天品质与土壤养分含量的相关性分析

从相关性分析结果(表5)可以看出，土壤养分含量与巴戟天营养元素含量密切相关，其中，土壤全磷、速效磷、速效钾含量与巴戟天全氮含量均表现出极显著正相关关系，土壤速效钾含量与巴戟天全磷含量之间以及土壤全钾含量与巴戟天全钾含量之间均具有极显著正相关关系。土壤养分含量与巴戟天寡糖含量密切相关，其中，土壤碱解氮含量与巴戟天蔗糖含量呈极显著正相关关系，但却与耐斯糖和1^F−果呋喃糖基耐斯糖含量呈显著负相关关系。土壤全磷、速效磷、速效钾含量与巴戟天1−蔗果三糖含量均表现出极显著正相关关系，与耐斯糖含量表现出极显著或显著正相关关系。值得注意的是，土壤pH与巴戟天营养元素没有显著相关关系，但与巴戟天寡糖呈显著相关关系，其中，土壤pH与蔗糖含量表现出极显著的负相关关系，与1−蔗果三糖、耐斯糖以及1^F−果呋喃糖基耐斯糖含量均有显著或极显著的正相关关系，说明随着土壤pH的增加蔗糖含量降低，而1−蔗果三糖、耐斯糖与1^F−果呋喃糖基耐斯糖含量却呈增加趋势。另外，土壤有机质含量与巴戟天品质指标之间相关关系不显著，说明土壤有机质含量对巴戟天品质影响不大。

2.4 不同生长年限巴戟天的品质指标含量

本次调查采集的是非林下种植环境的18个样地，巴戟天的生长年限包含1至6年生。由于1年生和6年生的仅有1个样点，故仅分析2~5年生巴戟天肉质根的营养元素和寡糖含量。

2.4.1 不同生长年限巴戟天药材中的营养元素含量

由图1可见，随着巴戟天生长年限的增长，其药用部分肉质根对全氮、全磷、全钾吸收积累量呈现先升高后降低的趋势，整体表现为全钾>全氮>全磷。2年生巴戟天的全氮含量显著低于其他生长年限的，其他处理间的全氮含量没有显著性差异。全磷含量、全钾含量在不同生长年限之间无显著差异，但均在3年生时达到最大值，分别为0.54、11.04 g·kg⁻¹。

2.4.2 不同生长年限巴戟天药材中的寡糖含量

由图2可以看出，巴戟天生长年限不同，寡糖类成分含量各异。随着巴戟天生长年限的增加，蔗糖、1−蔗果三糖含量变化均呈现出先增加后下降的趋势。其中，蔗糖含量为3年生巴戟天最高，质量分数达到64.08 g·kg⁻¹，显著高于其他生长年限的，同时，2年生巴戟天蔗糖含量显著高于4年生和5年生的，而4年生与5年生巴戟天的蔗糖含量之间无显著差异。在1-蔗果三糖含量上，4年生巴戟天显著高于2年生巴戟天(高出76.96%)，4年生巴戟天达到32.68 g·kg⁻¹，与3年生、5年生巴戟天没有显著差异，2年生、3年生与5年生巴戟天之间无显著差异。不同生长年限的巴戟天在耐斯糖和1^F−果呋喃糖基耐斯糖含量上没有显著差异。

2.5 林下与非林下种植巴戟天品质及其主成分分析 2.5.1 林下与非林下种植巴戟天的营养元素和寡糖含量

由于生产上巴戟天收获年龄为4年生，故本研究选择4年生的林下与非林下种植的巴戟天进行品质分析。林下与非林下种植的巴戟天全磷含量与全钾含量差异显著(图3)，林下种植的巴戟天全磷含量比非林下种植巴戟天高44.79%，而全钾含量比非林下巴戟天低13.78%，全氮含量无显著差异。

由图4可以看出，林下与非林下种植的巴戟天除耐斯糖含量没有显著差异外，其他寡糖含量间均有显著差异。其中，林下种植的4年生巴戟天蔗糖和1−蔗果三糖质量分数分别达到53.14和60.22 g·kg⁻¹，分别高出非林下种植巴戟天32.49%和18.68%。而1^F−果呋喃糖基耐斯糖含量却低于非林下巴戟天3.97%。

2.5.2 林下与非林下种植巴戟天品质的主成分分析

主成分分析是对多个变量进行统计分析的一类定量分析方法，通过分析主成分载荷值和得分可以看出不同样品之间的关系，也可以解释样品的特点、分组、相似性及差别^[21-22]。将巴戟天品质指标转化为7个主成分，分析结果如表6。7种指标可以提取出2种主成分，累计贡献率为85.776%。

为了更好地解释各品质指标与成分因子之间的关系，将所提取的主成分因子进行旋转处理，使一个变量在较少的几个因子上有较高的载荷，其载荷值的大小反映各有效成分指标在主成分中的重要程度(表7)。通过主成分分析可知，前2个主成分反映了7个指标的大部分综合信息，用2个主成分评价不同种植方式下巴戟天品质的优劣是可行的。因此，可用营养因子(F₁)、寡糖因子(F₂)代替原来7个指标(X₁~X₇)对巴戟天的品质进行分析，得到巴戟天品质前2个主成分的线性关系式分别为：

F₁=0.41X₁+0.52X₂−0.42X₃+0.36X₄+0.10X₅−0.28X₆−0.42X₇，

F₂= −0.15 X₁+0.17X₂−0.11X₃+0.45X₄+0.59X₅+0.51X₆+0.36X₇。

以每个主成分对应的特征值的方差提取贡献率建立综合评价模型：

F=0.15X₁+0.36X₂−0.28X₃+0.40X₄+0.32X₅+0.08X₆−0.06X₇。

计算林下与非林下种植方式下巴戟天品质的综合评分(表8)，得到不同种植方式下巴戟天品质特征的优劣顺序：7号样地林下种植巴戟天>6号样地林下种植巴戟天>6号样地非林下种植巴戟天>7号样地非林下种植巴戟天。说明，林下种植的巴戟天品质优于同条件下非林下种植的巴戟天。

3 讨论与结论 3.1 讨论

土壤理化性质是反映土壤养分的重要指标，土壤养分是植物生长发育的主要营养来源，同时也是创造优质高产必备的物质条件之一，研究其对巴戟天品质的影响有非常积极的意义。本研究表明，巴戟天品质与土壤养分关系密切，这与刘瑾等^[23-24]的研究结果一致。药理研究显示，巴戟天寡糖具有抗抑郁功效，其主要发挥功能的糖分为四~七糖^[25]，如耐斯糖(简称：四糖)、1−果呋喃糖基耐斯糖(简称：五糖)，土壤碱解氮与巴戟天蔗糖含量呈极显著正相关关系，但与耐斯糖和1^F−果呋喃糖基耐斯糖呈显著负相关关系，故减少土壤碱解氮的含量有利于巴戟天抗抑郁成分的积累。而土壤全磷、速效磷、速效钾对巴戟天1−蔗果三糖、耐斯糖均表现出显著的正相关关系，表明提高土壤磷和钾的有效性，有助于促进巴戟天寡糖累积。土壤的酸碱性对巴戟天品质的影响较大，巴戟天更适合种植于弱酸性土壤中^[26]，与研究发现弱酸性的红壤或棕壤较适宜巴戟天的生长一致^[15]。土壤pH与巴戟天根部营养元素含量无显著相关性，与蔗糖含量呈现出极强的负相关性，而与另外其他3种寡糖之间具有正相关关系，说明当土壤呈弱酸性时，会降低巴戟天根部的蔗糖含量，但会提高1−蔗果三糖、耐斯糖与1^F−果呋喃糖基耐斯糖含量。土壤有机质是土壤肥力的重要组成部分，能够使植株正常生长发育，保证植物在生长过程中的养分需求^[27]，但与巴戟天品质指标没有显著相关性。

不同生长年限的林下种植巴戟天的1−蔗果三糖和耐斯糖含量均大于同年生的非林下巴戟天，说明林下巴戟天从土壤中获取的养分较非林下巴戟天的多，从而使更多的单糖转化为1−蔗果三糖与耐斯糖。随着巴戟天生长年限的增加，巴戟天营养元素及寡糖含量呈现先增加后减少的变化趋势。这与玉竹Polygonatum odoratum、川党参Codonopsis pilosula^[28]等药用植物有效成分积累的趋势一致，这主要是因为巴戟天生长前期(第1~4年)通过植物光合作用不断贮存营养物质，而进入衰老期(第5年)后营养物质的消耗量大于累积量，导致有效成分含量下降。此外，巴戟天药材中发挥抗痴呆、抗抑郁、保护生殖系统等功效的成分主要是低聚糖^[29-31]，而低聚糖需要由单糖在一系列酶的催化下才能逐渐被合成，只有经过种植年限的积累，低聚糖含量才会升高，且巴戟天药材种植4年后其有效药用成分均能满足《中华人民共和国药典》^[20]巴戟天干燥品耐斯糖含量不得少于2.0%的要求，故4年生巴戟天获得的药材质量最佳。

通过对不同种植方式下同年生巴戟天品质特征进行分析，林下种植的4年生巴戟天全磷、蔗糖、1−蔗果三糖含量显著高于非林下种植的，全钾、1^F−果呋喃糖基耐斯糖含量却显著低于非林下种植的，可能是因为植物营养元素与有效成分并不是独立发挥作用，各成分之间具有一定相关性^[32]，从而导致不同种植方式对巴戟天品质特征物质含量有不同的促进作用。树林的郁闭度较大，为林下巴戟天提供了很好的遮荫效果，起到降温保湿的作用^[33]，这是林下种植的巴戟天有效成分含量优于非林下种植巴戟天的原因之一。

3.2 结论

土壤养分对巴戟天肉质根养分和寡糖含量有不同的影响，其中土壤pH与蔗糖含量表现出极显著的负相关性，与 1−蔗果三糖、耐斯糖以及1^F−果呋喃糖基耐斯糖含量呈显著正相关性。土壤全磷、速效磷、速效钾含量与1−蔗果三糖和耐斯糖含量呈显著正相关关系，表明在巴戟天生长过程中可以通过增加土壤速效磷与速效钾含量以提高其品质。不同生长年限的巴戟天肉质根全氮、全磷、全钾含量没有显著差异，而3年生巴戟天蔗糖含量显著高于其他生长年限的巴戟天，4年生巴戟天的1−蔗果三糖含量显著高于2年生巴戟天。林下与非林下种植的巴戟天肉质根全氮、全磷、全钾、蔗糖、1−蔗果三糖、耐斯糖和1^F−果呋喃糖基耐斯糖含量无一致规律，但主成分分析法综合评分结果显示，林下种植巴戟天品质优于非林下种植的巴戟天。