黄胜佳，叶 霜，刘新亚，席利娟，汪志辉，b，*，付佳玲

(四川农业大学 a. 园艺学院；b. 果蔬研究所，四川 成都 611130)

中国是柑橘果树的重要起源地之一，具有丰富的柑橘种质资源。黄果柑是橘和橙的天然杂交柑[1]，具有特晚熟、极丰产、花果同树等优良性状[2]。柑橘果实富含能参与人体新陈代谢、调节有关生理活动、对人体保健和疾病防治有重要作用的生物活性物质——酚类物质[3-5]。研究表明，酚类物质具有抗炎[6]、抗氧化[7]、抗过敏、抗癌、抗突变、抑菌等重要作用[8]。因此，研究不同肥料对柑橘酚类物质含量的影响，对提高柑橘综合利用价值具有积极意义。

国内外在测定不同柑橘品种、不同部位，以及不同生长时期果实中酚类物质含量等方面已有大量研究，但在环境因子，如土壤营养对柑橘生物活性物质代谢的影响方面的研究还比较少。已有研究表明，适量限制氮、磷肥的施用有利于黄酮类物质的积累[9]。氮素营养的缺乏会导致萜类、酚类等不含氮次生代谢产物的积累，而增加氮素营养会加速含氮次生代谢产物如生物碱、氰苷等的合成[10]。氮肥能显著影响西红柿幼苗中酚的含量[11]，也能影响枸杞中黄酮和类胡萝卜素等生物活性物质的含量。目前，关于肥料对植物酚类物质含量影响的研究大多注重于氮肥用量，而有关有机肥等肥料对其影响的研究报道较少。

灰色关联度分析法是一种针对灰色系统来决定因素主次及其关联程度的评价方法。本研究结合灰色关联度分析法评价不同施肥处理对黄果柑果实酚类物质的影响，以期为增加果实保健价值和提高柑橘的综合利用率提供参考。

1 材料与方法

1.1 试验基地与供试植株

试验于四川省雅安市石棉县小水村黄果柑示范园中进行。选取生长健壮、树势树形基本一致、株行距为3 m × 4 m、中等管理水平的9 a生黄果柑植株27株作为供试植株。

1.2 试验材料

有机肥：发酵后的鸡粪，N+P+K≥8%，有机质含量≥40%。

无机肥：尿素(N≥46.3%)，四川玖源农资化工有限公司生产；过磷酸钙(P2O5≥12.0%)，四川广汉市蜀峰化工有限责任公司生产；俄罗斯钾肥(K2O≥60%)，江西宜春市腾达化工有限责任公司生产。

有机无机复混肥：宝施灵黄果柑专用肥。具体配方为硫酸铵(10%)，尿素(19.6%)，磷酸一铵(15.9%)，硅钙钾肥(2.5%)，硝酸钾(1%)，硫酸钾(19%)，硫酸锌(2%)，钼酸铵(0.02%)，氨基酸(5%)，草炭(2%)，泥炭(2%)，堆肥(4%)，烟渣(2%)，青蒿(3%)，米皮(5%)，蘑菇渣(7%)。

标准品：柚皮素(naringenin)，柚皮苷(naringin)，新橙皮苷(neohesperidin)，咖啡酸(caffeic acid)，芥子酸(sinapic acid)，绿原酸(chlorogenic acid)，阿魏酸(ferulic acid)，没食子酸(gallic acid)，对羟基苯甲酸(p-hydroxybenzoic acid)，均购自Sigma公司，纯度均在98%以上；甲醇、甲酸和乙腈均为色谱纯，购自麦克林公司；水为实验室自制的超纯水。

1.3 试验设计

本试验共选择无机肥(A)、有机肥(B)、有机无机复混肥(C)3类肥料，每类肥料再分别设置低水平(1)、适量(2)、高水平(3)3种用量，共计9个处理，每处理重复3次，以1棵黄果柑作为1个重复。分别于幼果期(5月1日)、果实快速膨大期(8月5日)、转色期(11月10日)各施肥1次。根据文献[12]并结合当地黄果柑园地土壤及施肥情况，设计肥料施用量(表1)，其中，A2、B2、C2均为每种类型肥料中最适黄果柑生长发育的肥料用量。施肥时，围绕树冠滴水线的施肥圈，施肥圈要求深度和宽度达到20～30 cm，以见根不伤根为宜，在施肥圈内放入试验设计用量的肥料，施有机肥时土肥混匀。

1.4 采样

2016-04-15采收，分株单独采摘，每株按树冠外围的东西南北4个方位随机采取中等大小、成熟度相对一致、无病虫害且无霉烂的果实共12个，采回后将果实皮肉分离，置于-80 ℃冷冻保存备用。

1.5 柑橘酚类物质含量测定

1.5.1 样品前处理

从-80 ℃冰箱中取果肉至研钵中，加入液氮研磨成粉末。准确称取0.5 g混合样放入2 mL离心管中，避光条件下加入1.5 mL提取液(甲醇、水、甲酸体积比70%∶28%∶2%)，充分混匀。将离心管放入超声装置中辅助提取，超声30 min后，5 000 r·min-1离心15 min。取出离心管，吸取上清液。用0.45 μm微孔滤膜过滤，贮存于-20 ℃，用于酚酸及类黄酮分析。

表1各处理的施肥安排

Table1Fertilization arrangement of treatments

kg

1.5.2 标准曲线绘制

用色谱纯甲醇分别将每瓶定量的20 mg标准品溶解，定容于10 mL棕色容量瓶中，配制成2 g·L-1的标准品溶液。全程进行避光处理。将上述母液用甲醇分别配制成0.25、0.5、1、2、5、10、25、50、100 mg·L-1的浓度梯度，再配制成混合标准溶液，避光保存于-20 ℃以下环境中。

对不同质量浓度的混合标准溶液进行色谱分析，以标样的质量浓度为横坐标，峰面积为纵坐标，绘制出9种酚类物质的标准曲线，拟算回归方程和相关系数。

1.5.3 色谱条件

使用安捷伦1260液相色谱仪，C18色谱柱(5 μm，4.6 mm×150 mm)；柱温30 ℃；进样量20 μL；流动相A为10%甲酸溶液，流动相B为乙腈；梯度洗脱(洗脱程序如表2所示)；流速1 mL·min-1。

柚皮苷、新橙皮苷、柚皮素在283 nm波长下检测[13]；咖啡酸、芥子酸、绿原酸、阿魏酸在320 nm波长下检测；没食子酸、对羟基苯甲酸在260 nm波长下检测[14]。

表2流动相配比与洗脱时间

Table2Flow matching ratio and elution time

t/minA/%B/%095525851542782250643655955

1.6 数据分析

试验数据用IBM SPSS 19.0软件进行方差分析。综合评价基于Critic和熵权2种客观赋权法进行灰色关联度分析。

2 结果与分析

2.1 酚类物质标准品的色谱分析

图1为混合标准品55 min内分别在320、283、260 nm波长下的色谱图。根据保留时间长短，峰依次为没食子酸、绿原酸、咖啡酸、阿魏酸、芥子酸、柚皮苷、新橙皮苷、柚皮素。由图1可知，8种酚类物质出峰较好，说明本研究所用方法可满足多种酚类物质的同时分离。图2为对羟基苯甲酸在260 nm波长下的色谱图。由于对羟基苯甲酸出峰时间与绿原酸出峰时间相近，故分开测定。

根据色谱图拟合的回归方程及相关系数见表3。可以看出，各标准品的浓度与峰面积呈良好的线性相关关系，相关系数均超过0.999 1。

2.2 不同肥料对黄果柑酚酸类物质含量的影响

由表4可知，不同施肥处理黄果柑酚酸含量组成差异较大。各处理条件下，在黄果柑果肉中均检测到了绿原酸、咖啡酸、阿魏酸、芥子酸、没食子酸和对羟基苯甲酸。绿原酸含量以B2处理最高(4.765 μg·g-1)，显著(P<0.05)高于其他处理，以A3处理最低(2.885 μg·g-1)。咖啡酸含量以B1处理最高(3.529 μg·g-1)，显著(P<0.05)高于其他处理，以C1处理最低(2.692 μg·g-1)。阿魏酸含量以B1处理最高(1.101 μg·g-1)，显著(P<0.05)高于其他处理，以C1处理最低(0.838 μg·g-1)。芥子酸含量以B2处理最高(1.749 μg·g-1)，显著(P<0.05)高于其他处理，以C1处理最低(1.205 μg·g-1)。没食子酸含量以B1处理最高(3.507 μg·g-1)，显著(P<0.05)高于其他处理，以B3处理最低(1.595 μg·g-1)。对羟基苯甲酸含量以B2处理最高(6.261 μg·g-1)，显著(P<0.05)高于其他处理，以C1处理最低(4.262 μg·g-1)。总体来看，酚酸类物质含量最高的主要出现在B1、B2处理。这说明，低水平或适量的有机肥能显著提高果肉中酚酸物质含量。在施用无机肥的3个处理中，A1处理果肉中酚酸含量显著(P<0.05)高于A2、A3处理，且仅次于B1、B2处理。这说明，施用低水平平氮磷的无机肥亦能取得较好效果。

1， 没食子酸； 2， 绿原酸； 3， 咖啡酸； 4， 阿魏酸； 5， 芥子酸； 6， 柚皮苷； 7， 新橙皮苷； 8， 柚皮素。1， Gallic acid; 2， Chlorogenic acid; 3， Caffeic acid; 4， Ferulic acid; 5， Sinapic acid; 6， Naringin; 7， Neohesperidin; 8， Naringenin.图1 混合标样在不同波长下的色谱图Fig.1 Chromatogram of mixed standard samples under different wavelength

1，对羟基苯甲酸 p-hydroxybenzoic acid.图2 对羟基苯甲酸在260 nm波长下的色谱图Fig.2 Chromatogram of p-hydroxybenzoic acid under 260 nm wavelength

表3标准曲线及相关系数

Table3Standard curve and correlation coefficient

样品名称回归方程相关系数SampleRegressionequationCorrelationcoefficient咖啡酸Caffeicacidy=139.21x-86.5330.9996阿魏酸Ferulicacidy=107.09x-27.8130.9999绿原酸Chlorogenicacidy=35.869x-9.5490.9992芥子酸Sinapicacidy=100.16x-11.5231.0000新橙皮苷Neohesperidiny=41.455x-22.920.9993柚皮素Naringeniny=72.394x-53.7710.9991柚皮苷Naringiny=27.071x-0.7560.9999没食子酸Gallicacidy=44.696x+21.6050.9997对羟基苯甲酸p-Hydroxybenzoicacidy=113.19x-16.5810.9996

表4黄果柑果肉中酚酸物质含量

Table4Quantity of phenolic acid in Huangguogan sarcocarp

μg·g-1

表中所有含量均以鲜质量计。同列数据后无相同小写字母的表示处理间差异显著(P<0.05)。下同。

Contents listed in the above table were calculated based on fresh weight. Data within the same column followed by no same letters indicated significant difference atP<0.05. The same as below.

2.3 不同肥料对黄果柑中类黄酮物质含量的影响

由表5可知，各处理果实中均检测到了柚皮苷和新橙皮苷，但均未检测到柚皮素。各处理柚皮苷含量在2.026～4.356 μg·g-1之间，新橙皮苷含量在0.546～0.757 μg·g-1之间。对比果肉中柚皮苷的含量可知，B3处理含量最高，且显著(P<0.05)高于其他处理。同时，以B3、B2和A1处理的新橙皮苷含量较高。

2.4 基于Critic与熵权法的灰色关联度分析

2.4.1 熵权法分析

熵原本是热力学的概念，主要用于形容指标信息的离散程度。一般地，如果某个指标的信息熵(Ej)越小，表明其指标值的变异程度越大，提供的信息量越大，在综合评价中所起的作用越大,其权重也应越大；反之，某指标的信息熵越大，表明其指标值的变异程度越小，提供的信息量越小，在综合评价中所起的作用越小，其权重也应越小。

表5黄果柑果肉中类黄酮含量

Table5Quantity of flavonoid in Huangguogan sarcocarp

μg·g-1

ND, 未检出。

ND, Not detected.

根据熵权法的计算步骤[15]，对不同酚类物质的信息熵和权重熵(ωj)进行测算(因柚皮素在所有处理下都未检测到，故不予考虑)。由表6可知，柚皮苷、没食子酸和绿原酸3个指标的熵值低，差异性系数高，差异梯度较大，因此赋予高权重。柚皮苷权重值最大，为0.296 0；其次为没食子酸(0.292 7)，绿原酸(0.128 4)，阿魏酸权重值最小，为0.034 5。所得指标权重排序为柚皮苷>没食子酸>绿原酸>对羟基苯甲酸>芥子酸>咖啡酸>新橙皮苷>阿魏酸。

2.4.2 Critic法分析

Critic法是由Diakoulaki等[16]提出的另一种客观权重赋权的方法。该方法基于评价指标的对比强度和指标之间的冲突性来综合衡量指标的客观权重。对比强度，是指同一个指标各个评价方案之间取值差距的大小，以标准差的形式来表现；指标之间的冲突性，以指标之间的相关性为基础，若2个评价指标之间具有较强的正相关性，则这2个指标的冲突性较低。用Cj表示第个评价指标所包含的信息量，值越大，说明其所包含的信息量越大，则其权重(ωj)亦越大。

表6各指标的权重(熵权法)

Table6Index weights (Entropy weight method)

样品SampleEjωj绿原酸Chlorogenicacid0.99430.1284咖啡酸Caffeicacid0.99790.0477阿魏酸Ferulicacid0.99850.0345芥子酸Sinapicacid0.99680.0708没食子酸Gallicacid0.98690.2927对羟基苯甲酸p-Hydroxybenzoicacid0.99620.0847柚皮柑Naringin0.98680.2960新橙皮苷Neohesperidin0.99800.0453

由表7可知，没食子酸权重值最大(0.252 2)，其次为柚皮苷(0.206 2)、绿原酸(0.195 8)，阿魏酸权重最小(0.025 1)。所得指标权重排序为没食子酸>柚皮苷>绿原酸>对羟基苯甲酸>咖啡酸>芥子酸>新橙皮苷>阿魏酸。

2.4.3 Critic法与熵权法结合分析

Critic和熵权法计算指标权重都属于客观赋权的方法，其中，熵权法主要表现指标数据的离散性，Critic法则既兼顾到指标数据的相关性，又考虑到数据之间的对比强度。如果将这2种方法相结合，则可以充分考虑到数据的离散性、相关性和对比强度，是一种更加完善的进行客观赋权的方法，计算公式如下：

(1)

式(1)中：ωj是第j个指标权重；qij是指标i与j的相关系数；σj是指标j的标准差；ej是指标j的熵。

基于Critic法与熵权法相结合的权重赋值结果如表8所示，没食子酸权重值最大(0.195 2)，其次为新橙皮苷(0.178 3)、绿原酸(0.140 3)，芥子酸权重值最小(0.075 5)。所得指标权重排序依次为没食子酸>新橙皮苷>绿原酸>阿魏酸>对羟基苯甲酸>柚皮苷>咖啡酸>芥子酸。

表7各指标的权重(Critic法)

Table7Index weights (Critic method)

样品SampleCjωj绿原酸Chlorogenicacid1.36830.1958咖啡酸Caffeicacid0.45210.0647阿魏酸Ferulicacid0.17550.0251芥子酸Sinapicacid0.23980.0343没食子酸Gallicacid1.76250.2522对羟基苯甲酸p-Hydroxybenzoicacid1.34220.1921柚皮柑Naringin1.44090.2062新橙皮苷Neohesperidin0.20710.0296

表8各指标的权重(基于Critic与熵权法的客观赋权法)

Table8Index weights (Objective weighting method based on Critic and Entropy method)

样品Sampleωj绿原酸Chlorogenicacid0.1403咖啡酸Caffeicacid0.0850阿魏酸Ferulicacid0.1172芥子酸Sinapicacid0.0755没食子酸Gallicacid0.1952对羟基苯甲酸p-Hydroxybenzoicacid0.1075柚皮柑Naringin0.1009新橙皮苷Neohesperidin0.1783

2.4.4 灰色关联度分析

将灰色关联分析与客观赋权法相结合，可以克服采用专家赋权法等主观赋权的不确定性。对试验数据进行以下处理：根据实际情况，确定黄果柑主要酚类物质的理想值为参试施肥处理的最高值，组成参考数列X0。为便于分析，采用初值化法对原始数据进行无量纲化处理。求关联系数，进而分别计算求得关联度和关联序。对于不同酚类物质，其相对重要程度不同，采用基于Critic与熵权法的客观赋权法的结果赋予关联系数不同的权重值ωj，加权关联度对各施肥处理进行评价，结果见表9。B1处理的关联度最大，最接近参考数列，即各酚类物质在B1处理下都具有较高含量，其次为B2和A1处理。

表9不同施肥处理的关联度

Table9Relational grade of different fertilizer treatments

处理Treatment关联度RelateddegreeA10.8028A20.6518A30.6187B10.9122B20.8859B30.7447C10.5426C20.5690C30.7505

3 讨论

本研究发现，有机肥、无机肥、有机无机复混肥的施用量对黄果柑酚类物质含量的影响不仅与肥料中氮磷含量密切相关，还与肥料种类相关。单施无机肥时，低水平的氮磷施入量较有利于酚酸和类黄酮物质的积累，这与Copley等[10]、Strissel等[22]的研究结果一致。单施有机肥时，适量或低水平有机肥对增加果实中酚酸含量的效果较好，但当施用高水平有机肥时效果不理想，而适量或高水平有机肥对增加果实中类黄酮物质的效果较好，施用低水平有机肥时效果不明显。有机无机复混肥的用量则以高水平较好，能显著增加果肉中酚酸和类黄酮物质含量。本研究在黄果柑果肉中未检测到柚皮素，与张元梅[14]结果一致。

为准确评价不同施肥处理对黄果柑果实酚类物质的影响，本研究利用基于Critic与熵权法的客观赋权进行灰色关联度分析，结果显示，施用低水平的有机肥对黄果柑果肉中酚类物质含量的提高效果最好，施用适量有机肥和低水平氮磷无机肥的效果次之。因此，在生产实际中，建议适当减少无机肥中氮磷施用量，或适当增施有机肥和有机无机复混肥，以增加黄果柑果肉中酚类物质含量。本研究未涉及不同施肥水平对黄果柑果实发育和果实品质的影响，在今后的研究中，应进一步结合不同施肥水平对酚类物质、果实生长发育及果实品质的影响，全面分析，以确定最适肥料种类和施肥量。

[1] 汪志辉， 刘世福， 严巧巧， 等. 石棉县黄果柑生物学特性调查与差异株系比较[J]. 北方园艺， 2011(14):20-24.

WANG Z H， LIU S F， YAN Q Q， et al. Biological characteristics survey and the differences compare in strains of Huangguogan in Shimian[J].NorthernHorticulture， 2011 (14): 20-24. (in Chinese with English abstract)

[2] 王大华， 李天眷. 黄果柑优良株系：桂晚柑的选育[J]. 中国南方果树， 1991， 20(3): 3-5.

WANG D H， LI T J. Selection of Guiwangan， an excellent line of Huangguogan[J].SouthChinaFruits， 1991， 20(3): 3-5. (in Chinese with English abstract)

[3] 张元梅， 周志钦. 柑橘生物活性物质及其心血管疾病防治作用研究进展[J]. 中药材， 2011， 34(11): 1799-1803.

ZHANG Y M， ZHOU Z Q. Studies on citrus bioactive substances and their prevention and treatment of cardiovascular diseases[J].JournalofChineseMedicinalMaterials， 2011， 34(11): 1799-1803. (in Chinese)

[4] DAI J， MUMPER R J. Plant phenolics: extraction， analysis and their antioxidant and anticancer properties[J].Molecules， 2010， 15(10): 7313-7352.

[5] 徐贵华， 关荣发， 叶兴乾， 等. 不同成熟期蜜桔中酚酸的组成与分布[J]. 食品科学， 2008， 29(2): 137-141.

XU G H， GUAN R F， YE X Q， et al. Composition and distribution of phenolic acids in satsuma mandarin (CitrusunshiuMarc.) during maturity[J].FoodScience， 2008， 29(2): 137-141. (in Chinese with English abstract)

[6] 朱晓艳. 瓯柑果实黄酮类化合物的分离鉴别及相关活性研究[D]. 杭州: 浙江大学， 2012.

ZHU X Y. Separation and identification of flavonoids fromCitrusreticulatecv. Suavissima fruit and their bioactivity study[D]. Hangzhou: Zhejiang University， 2012. (in Chinese with English abstract)

[7] RORELLI A E， GUARDIAT T， JUAREZ A O， et al. Comparative study of flavonoids in experimental models of inflammation[J].PharmacologicalResearch， 2003， 48(6): 601-606.

[8] WILMSEN P K， SPADA D S， SALVADOR M. Antioxidant activity of the flavonoid hesperidin in chemical and biological systems[J].JournalofAgricultural&FoodChemistry， 2005， 53(12): 4757-4761.

[9] 徐文燕， 高徽徽， 何春年. 环境因子对植物黄酮类化合物生物合成的影响[J]. 世界科学技术， 2006， 8(6): 68-72.

XU W Y， GAO H H， HE C N. The influences of environmental factors on flavonoid biosynthesis[J].WorldScienceandTechnology， 2006， 8(6): 68-72. (in Chinese with English abstract)

[10] COPLEY R R， LETUNIC I， BORK P， et al. Genome and protein evolution in eukaryotes[J].CurrentOpinioninChemicalBiology， 2002， 6(1): 39-45.

[11] STOUTM J. Effect of nitrogen availability on expression of constitutive and inducible chemical defenses in tomato，Lycopersiconesculentum[J].JournalofChemicalEcology， 1998， 24(6): 945-963.

[12] 李清南， 曹淑燕， 古咸杰， 等. 两种施肥模式下柑橘树体氮素吸收与运转机制的研究[J]. 浙江农业学报， 2016， 28(1): 51-55.

LI Q N， CAO S Y， GU X J， et al. Study on nitrogen absorption and operation mechanism of citrus tree in two fertilization modes[J].ActaAgriculturaeZhejiangensis， 2016， 28(1): 51-55. (in Chinese with English abstract)

[13] 郑洁， 赵其阳， 张耀海， 等. 超高效液相色谱法同时测定柑橘中主要酚酸和类黄酮物质[J]. 中国农业科学， 2014， 47(23): 4706-4717.

ZHENG J， ZHAO Q Y， ZHANG Y H， et al. Simultaneous determination of main flavonoids and phenolic acids in citrus fruit by ultra performance liquid chromatography[J].ScientiaAgriculturaSinica， 2014， 47(23): 4706-4717. (in Chinese with English abstract)

[14] 张元梅. 中国野生宽皮柑橘酚类物质的检测及其抗氧化活性研究[D]. 重庆:西南大学， 2013.

ZHANG Y M. Phenolic compositions and antioxidant capacities of the Chinese wild mandarin (CitrusreticulateBlanco) fruits[D]. Chongqing: Southwest University， 2013. (in Chinese with English abstract)

[15] 杨阳， 方春晖， 李占超. 基于最大熵的坝体健康评估专家模糊赋权模型[J]. 长江科学学院院报， 2011， 28(3): 20-23.

YANG Y， FANG C H， LI Z C. Experts fuzzy weight model of dam health assessment based on maximum entropy[J].JounalofYangtzeRiverScientificResearchInstitute， 2011， 28(3): 20-23. (in Chinese with English abstract)

[16] DIAKOULAKI D， MAVROTAS G， PAPAYANNAKIS L， et al. Determining objective weights in multiple criteria problems: the critic method[J].Computers&OperationsResearch， 1995， 22(7): 763-770.

[17] STRISSEL T， HALBWIRTH H， HOYER U， et a1. Growth-promoting nitrogen nutrition affects flavonoid biosynthesis in young apple (MalusdomesticaBorkh.) leaves[J].PlantBiology， 2005， 7(6): 677-685.