桃溪蜜柚是江西省地方特色柚类——井冈蜜柚的代表性品种之一,成熟期为每年9月中下旬,具有早结丰产、果肉脆嫩、化渣,汁多,风味甜爽适口等突出优点[1],近年来栽培面积和产量不断上升[2]。尽管栽培体量大,但由于生产上对其采收成熟度及其贮藏特性认识不够,采后处理措施不当等原因,造成了严重的经济损失。因此,利用科学的方法合理地评价不同成熟度的桃溪蜜柚果实品质及耐贮性差异,可为适时采收、贮藏、出库、销售等决策提供理论参考。
前期研究对不同采收成熟度的3个井冈蜜柚代表品种——桃溪蜜柚、金沙柚和安福金兰柚果实品质相关指标做了检测,并比较分析了不同成熟度果实在贮藏期间的品质变化[3-5]。发现不同采收成熟度的果实其腐烂率、失重率、色泽、风味、抗氧化性能等均表现出明显的差异,总体应根据鲜销、短期贮藏或长期贮藏等目的选择合适的采收成熟度[5]。但由于评价贮藏品质的指标多样,如外观品质、内在风味品质及抗性生理等相关理化指标体系,且每个指标代表的生物学含义及其变化规律纷繁芜杂,独立分析某一或某些相关指标往往缺乏系统性,难以全面直观地反映整体差异情况[3]。主成分分析法(principal components analysis,PCA)是一种通过数据压缩对多维数据进行降维的多元统计分析方法,它的优势在于操作简单,没有参数限制,应用范围较广,常被用于人脸识别、图像压缩、特征提取等领域[6],近年来,这一方法在果蔬品质分析中的应用也越来越多[7-11]。
选择适宜成熟度进行采收、贮藏和出库是保证桃溪蜜柚果实品质和实现产值提升的首要步骤。前期李宏祥等[4]采集了盛花后 178 d(成熟度Ⅰ)、盛花后 192 d(成熟度Ⅱ)和盛花后 208 d(成熟度Ⅲ)的桃溪蜜柚果实,并于常温贮藏后0、30、60、90、120、150和180 d,定期检测了失重率、腐烂率、L*、a*、b*、C*、Ho、色泽指数(color contribution index,CCI)、总糖、可溶性固形物(total soluble solids ,TSS)、可滴定酸(titratable acidity,TA)固酸比(TSS/TA)、抗坏血酸、可食率、出汁率、MDA含量、多酚氧化酶(polyphenol oxidase,PPO)活性、总酚含量、总黄酮含量等19个贮藏特性相关指标。但由于实验因素水平和检测指标较多的原因(采收成熟度3个水平、贮藏周期7个水平、贮藏特性指标19个),导致数据量大且关系复杂,尽管逐一对各指标变化趋势进行了具体分析,但未能明确阐释采收成熟度对贮藏特性的总体影响。因此,为进一步了解盛花后 178、192和208 d采收的桃溪蜜柚贮藏特性的差异,本试验通过PCA将19个指标浓缩为5个主成分,在此基础上分析了品质指标变化规律,并对3个采收成熟度的果实品质进行了综合评分。以分值为依据,对不同成熟度果实的适宜贮藏周期和出库时间提出了意见。
实验原始数据来源于李宏祥等[4]已发表的文章,详见表1。
表1 不同采收成熟度桃溪蜜柚果实贮藏性和品质指标动态变化
Table 1 Dynamic changes of storability and quality indexes for Taoxi pomelo at different harvesting maturities
采收期贮藏期/d失重率/%腐烂率/%L∗a∗b∗C∗HoCCI总糖/(mg·g-1)TSS/%Ⅰ00.00±0.000.00±0.0062.71±0.86c-17.16±0.66c45.89±0.89b49.00±0.64b110.55±1.09a-5.97±0.42c73.12±2.21a10.93±0.06aⅠ301.76±0.21b0.00±0.0066.33±0.99c-13.97±0.21c49.61±0.38a51.54±0.32a105.74±0.32a-4.24±0.07c67.49±2.20b10.30±0.10aⅠ603.47±0.40a0.00±0.0073.76±0.57b-4.12±0.17c53.15±2.14b53.30±2.15b94.43±0.04a-1.05±0.02c69.03±1.09b10.73±0.06bⅠ904.54±0.52a0.00±0.0074.35±0.92a-1.62±0.21b60.16±1.52b60.19±1.51b91.56±0.25b-0.36±0.05b65.49±0.93c11.93±0.06bⅠ1205.28±0.58a1.00±1.73b73.38±0.39b1.20±0.21b67.06±2.93a67.07±2.93a88.97±0.22b0.25±0.05b67.57±0.23b11.33±0.06cⅠ1505.79±0.62ab1.00±1.73b72.90±0.89a3.15±0.43b67.06±1.85a67.13±1.83a87.30±0.43b0.65±0.11b68.03±0.83b12.50±0.00aⅠ1806.79±0.73a1.00±1.73c68.49±1.57b4.16±0.27a61.04±1.48a61.18±1.49a86.12±0.20c1.00±0.05a54.84±1.06a10.17±0.06bⅡ00.00±0.000.00±0.0070.13±0.20b-10.59±0.20b50.39±1.05a51.49±0.98a101.86±0.47b-3.00±0.13b65.64±2.87b11.17±0.06aⅡ302.23±0.04a2.02±1.75a72.20±0.45b-6.44±0.37b48.35±2.46a48.78±2.42a97.59±0.67b-1.85±0.17b66.80±2.60b11.97±0.06bⅡ603.54±0.08a3.02±0.02a72.37±0.69c0.49±0.17a58.75±2.27a58.75±2.28a89.53±0.15c0.11±0.04a65.87±1.75c12.23±0.06bⅡ904.59±0.09a3.02±0.02a72.41±1.19b3.39±0.18a64.97±1.82a65.06±1.81a87.02±0.23c0.72±0.06a69.58±0.93b12.50±0.10bⅡ1205.37±0.09a6.03±3.08a72.50±0.99b4.04±0.35a66.19±2.15a66.31±2.16a86.52±0.21c0.84±0.04a65.56±1.31b12.43±0.06bⅡ1505.95±0.09a6.03±3.08a71.64±0.72a3.87±0.15a64.33±1.04a64.45±1.03a86.57±0.18c0.84±0.04a75.98±2.49a11.93±0.06aⅡ1807.16±0.06a6.03±3.08a71.38±0.91a3.20±0.23a62.56±0.66a62.64±0.66a87.10±0.20b0.72±0.05b72.51±2.50a11.23±0.06aⅢ00.00±0.000.00±0.0076.02±0.77b-5.36±0.12a44.88±1.08b45.20±1.08c96.83±0.06c-1.57±0.01a65.33±1.86b10.33±0.06bⅢ301.71±0.10b0.00±0.0075.18±1.18b-3.35±0.12a42.88±2.20b43.01±2.18b94.51±0.38c-1.04±0.07a77.45±0.93a12.10±0.00bⅢ602.89±0.13b0.00±0.0076.47±0.13a-2.90±0.11b50.26±1.66b50.34±1.67b93.52±0.05b-0.75±0.01b74.59±1.09a11.60±0.00aⅢ903.68±0.16b0.00±0.0075.94±0.73a-2.72±0.04c51.86±1.93c51.93±1.93c93.03±0.08a-0.69±0.01c85.16±1.45a12.47±0.06aⅢ1204.17±0.23b1.00±1.73b75.25±0.36a-1.69±0.07c57.30±1.76b57.32±1.76b91.71±0.07a-0.39±0.02c78.91±1.52a12.57±0.06aⅢ1505.01±0.28b1.01±1.75b72.57±0.50a-0.43±0.09c56.40±3.08b56.40±3.09b90.47±0.05a-0.10±0.02c74.28±1.90a11.93±0.06aⅢ1806.71±0.32a3.98±3.45b71.49±1.70a1.92±0.77b60.64±2.10a60.67±2.09a88.23±0.76a0.45±0.19c64.23±2.43b10.67±0.06a采收期贮藏期/dTA/%固酸比/%抗坏血酸/[mg·(100 g)-1FW]可食率/%出汁率/%MDA/(mmol·g-1 W)PPO/(U·g-1 W)总酚/(U·g-1 )总黄酮/(U·g-1 )Ⅰ00.50±0.01b21.96±0.62a58.51±0.60a34.87±0.71a56.72±3.11a0.01±0.00ab1.23±0.07b171.23±1.73a163.40±0.00aⅠ300.55±0.03b18.69±1.07a57.96±0.60a36.46±2.62b54.50±4.77a0.01±0.00b1.44±0.04a164.97±0.92b152.86±10.55aⅠ600.71±0.03b15.26±0.58b52.28±0.60b37.43±1.26a54.86±6.27a0.02±0.00b1.47±0.04a172.12±0.93b170.43±12.17aⅠ900.75±0.03a15.98±0.56b52.41±0.60c41.01±4.29a39.93±2.84b0.02±0.00c1.19±0.08b170.15±1.97c247.73±65.83aⅠ1200.70±0.02a16.25±0.29c54.15±0.00a35.46±6.73a63.35±2.80a0.02±0.00b1.27±0.04a168.02±0.62c159.89±6.09aⅠ1500.67±0.03a18.69±0.73b51.02±0.00a37.56±2.36a50.21±4.39b0.03±0.00ab1.33±0.05a207.48±2.09a254.76±32.21aⅠ1800.46±0.04a22.26±1.77b43.38±0.60b36.17±1.60b37.84±4.27b0.03±0.00c1.19±0.10a192.85±3.82a216.11±0.00aⅡ00.52±0.02b21.34±0.51a51.37±0.60b38.86±1.68a56.67±11.76a0.02±0.00a1.50±0.12a164.50±0.42b149.34±12.18aⅡ300.72±0.01a16.56±0.13b53.80±0.61c39.01±0.89ab57.67±4.10a0.01±0.00a1.19±0.04b168.76±0.91a138.80±6.09abⅡ600.84±0.05a14.51±0.92b59.00±0.60a39.63±6.65a42.36±9.50a0.02±0.00a0.89±0.09b178.93±0.51a163.40±10.54aⅡ900.79±0.03a15.86±0.42b57.62±0.60a36.52±2.19a43.82±1.79b0.03±0.00a1.11±0.05b174.08±0.64b138.80±37.02bⅡ1200.66±0.01b18.79±0.13b54.15±1.05a38.71±1.10a34.63±7.40b0.02±0.00a1.04±0.14b177.21±1.12b159.89±6.09aⅡ1500.59±0.01b20.28±0.26a49.63±0.60b36.85±7.24a48.92±5.30b0.02±0.00b1.31±0.06a188.59±0.18b173.94±27.89bⅡ1800.36±0.01b30.97±1.05a45.12±0.61a39.47±1.55a38.97±10.40b0.03±0.00a1.19±0.05a193.87±3.16a216.11±21.09aⅢ00.57±0.01a18.11±0.25b46.86±0.00c36.94±3.15a55.83±4.91a0.01±0.00b1.34±0.04ab149.68±1.15c96.64±6.09bⅢ300.70±0.02a17.26±0.36b56.23±0.00b41.47±0.90a60.73±2.57a0.01±0.00ab1.23±0.10b167.67±0.71a128.26±6.09bⅢ600.62±0.01c18.62±0.36a59.70±0.60a40.33±3.05a44.00±1.31a0.02±0.00b1.43±0.08a171.19±0.62b135.28±6.09bⅢ900.69±0.01b18.05±0.28a54.50±0.60b36.90±3.99a53.86±2.74a0.02±0.00b1.34±0.07a180.85±2.67a156.37±16.10bⅢ1200.61±0.01c20.68±0.51a55.54±0.60a37.21±3.84a44.28±3.95b0.01±0.00c1.38±0.07a183.31±1.49a173.94±18.26aⅢ1500.58±0.02b20.55±0.53a50.98±0.59a41.64±1.65a62.84±4.06a0.03±0.00a1.38±0.08a186.32±1.53b177.45±12.17bⅢ1800.36±0.00b29.36±0.23a45.59±1.04a39.51±1.03a57.19±7.34a0.03±0.00b1.30±0.06a163.56±2.08b166.91±12.17b
注:同一列不同小写字母表示具有显著性差异(P<0.05)
所有原始数据利用Microsoft Excel 2010进行预处理;利用SPSS 26.0 对数据进行单因素方差分析及主成分分析;利用Graphpad Prism 8.0对数据统计分析和作图。
盛花后 178 d(成熟度Ⅰ)、192 d(成熟度Ⅱ)和208 d(成熟度Ⅲ)的桃溪蜜柚果实,常温贮藏0、30、60、90、120、150和180 d内的失重率、腐烂率、L*、a*、b*、C*、Ho、CCI、总糖、TSS、TA、固酸比、抗坏血酸、可食率、出汁率、MDA含量、PPO活性、总酚含量、总黄酮含量等19个贮藏特性相关指标原始数据如表1所示。
为消除所检测的19个品质指标原始数据量纲和数量级的差异,首先利用SPSS 26.0将21组(3个采收成熟度×7个贮藏时期)原始数据进行均一化预处理。均一化公式:X=(原始值-平均值)/标准偏差,最终使各个指标的平均数为0,标准偏差为1。各指标数据均一化结果如表2所示。
表2 各品质指标均一化后的值
Table 2 The normalized value of each index
成熟度贮藏期/d失重率腐烂率L∗a∗b∗C∗HoCCI总糖TSSTA固酸比抗坏血酸可食率出汁率MDA含量PPO活性总酚含量总黄酮含量Ⅰ0-1.72-0.77-2.89-2.52-1.35-1.052.67-2.910.5-0.8-0.910.581.19-1.690.72-1.1-0.28-0.37-0.14Ⅰ30-0.93-0.77-1.8-1.98-0.87-0.71.95-1.95-0.37-1.58-0.52-0.21.08-0.880.47-1.221.13-0.85-0.41Ⅰ60-0.17-0.770.45-0.33-0.42-0.470.24-0.18-0.13-1.050.7-1.02-0.12-0.380.51-0.521.31-0.30.05Ⅰ900.31-0.770.630.090.490.46-0.190.21-0.680.451.01-0.85-0.091.44-1.2-0.15-0.57-0.462.08Ⅰ1200.64-0.310.340.561.381.39-0.580.55-0.36-0.30.62-0.780.27-1.391.48-0.290-0.62-0.23Ⅰ1500.87-0.310.190.891.381.4-0.830.77-0.281.150.42-0.2-0.39-0.32-0.030.750.372.442.26Ⅰ1801.32-0.31-1.141.060.60.6-1.010.96-2.32-1.75-1.250.66-2-1.03-1.441.16-0.571.311.25Ⅱ0-1.72-0.77-0.64-1.42-0.77-0.711.36-1.26-0.65-0.51-0.730.43-0.310.340.72-0.781.51-0.89-0.51Ⅱ30-0.720.16-0.02-0.72-1.03-1.080.72-0.62-0.470.490.83-0.710.20.420.83-0.9-0.56-0.56-0.78Ⅱ60-0.130.620.030.440.310.27-0.490.47-0.620.821.77-1.21.30.73-0.930.29-2.610.23-0.14Ⅱ900.340.620.050.931.111.12-0.870.81-0.051.151.35-0.881.01-0.85-0.761.05-1.11-0.15-0.78Ⅱ1200.692.020.071.041.271.29-0.950.87-0.661.070.36-0.180.270.27-1.810.49-1.560.1-0.23Ⅱ1500.942.02-0.191.011.031.04-0.940.870.940.45-0.210.18-0.68-0.68-0.170.420.230.980.14Ⅱ1801.482.02-0.270.890.80.8-0.860.810.4-0.42-1.982.74-1.630.65-1.311.87-0.541.391.25Ⅲ0-1.72-0.771.14-0.54-1.48-1.560.61-0.46-0.7-1.54-0.37-0.34-1.26-0.640.62-1.650.42-2.04-1.89Ⅲ30-0.95-0.770.88-0.2-1.74-1.860.26-0.171.160.660.65-0.540.711.671.18-1.07-0.33-0.64-1.06Ⅲ60-0.43-0.771.27-0.13-0.79-0.870.08-0.010.730.030.02-0.221.451.09-0.74-0.471.08-0.37-0.87Ⅲ90-0.07-0.771.11-0.1-0.58-0.650.030.022.351.110.57-0.350.35-0.650.390.60.450.38-0.32Ⅲ1200.15-0.310.910.080.120.08-0.170.191.391.24-0.050.280.57-0.5-0.71-0.930.720.570.14Ⅲ1500.52-0.310.090.290-0.05-0.350.350.680.45-0.290.25-0.391.751.420.90.740.80.23Ⅲ1801.281.07-0.230.680.550.53-0.690.66-0.87-1.13-22.36-1.530.670.771.560.17-0.96-0.04
2.2.1 主成分提取
为判断均一化后的21组原始数据是否适合PCA聚类分析,首先对均一化后的数据进行KMO和Bartlett检验。结果见表3,尽管Bartlett球形度检验对应的P值(8.6 E-55)<0.05,说明此组数据适合进行主成分分析,但KMO值仅为 0.322,且降维旋转在25次迭代之后无法收敛。因此,为提高数据模型的可靠性,剔除“可食率”和“出汁率”2个统计学上差异较小的指标,将剩余17个指标进行下一步主成分分析。
17个维度的品质指标数据KMO和巴利特检验结果如表3所示,KMO值为0.487,显著性为0.000,说明这17个指标可以较好地进行主成分分析。
表3 KMO和巴利特检验
Table 3 KMO and Bartlett test
19个指标的KMO和巴利特检验17个指标的KMO和巴利特检验KMO 取样适切性量数0.322KMO 取样适切性量数0.487Bartlett 的球形度检验近似卡方636.40Bartlett 的球形度检验近似卡方590.26df171df136Sig.8.6E-55Sig.0.000
利用PCA降维,将贮藏期品质变化的总和划分为若干个主成分,每个主成分的特征值均反映了原始变量信息的数量,表4展示了方差贡献率较大的几个主成分。其中主成分1(PC1)和主成分2(PC2)中的方差贡献率最大,其特征值分别为7.73和3.42,分别可以解释总方差45.46%和20.12%的信息量,累积解释总方差65.58%的信息量。
表4 主成分的特征值和方差贡献率
Table 4 Eigenvalue of the principal components and
their contribution and cumulative contribution
主成分特征值方差贡献率/%累计贡献率/%17.72845.45945.45923.42120.12165.5831.6119.47975.05941.4848.72883.78751.1796.93390.71960.5583.2893.999
2.2.2 基于主成分得分分析不同贮藏周期及果实品质指标变化
根据各个品质指标在PC1和PC2上的得分绘制二维散点图(图1),分布在PC1 [0,1] 区间的指标有:失重率、腐烂率、L*、a*、b*、C*、CCI、TSS、固酸比、MDA、总酚、总黄酮,代表了较高的失重率、腐烂率、膜脂过氧化、总酚、总黄酮和较好的色泽;分布在PC1 [-1,0] 区间的指标有Ho、总糖、TA、抗坏血酸、PPO活性,代表了较低的PPO活性和抗坏血酸含量。分布在PC2 [0,1] 区间的指标有L*、a*、b*、C*、CCI、总糖、TSS、TA、抗坏血酸,代表了较高的糖酸含量和较好的色泽;分布在PC2 [-1,0] 区间的指标有失重率、腐烂率、Ho、固酸比、MDA、PPO活性、总酚、总黄酮,代表了较低的失重率、腐烂率、固酸比、膜脂过氧化、总酚、总黄酮。
将3个采收成熟度果实在不同贮藏时期的得分情况映射到图1中。可以看出,3个采收成熟度的0 d均分布于第3象限,180 d均分布于第4象限,两者分别游离于其他几个时期,且附近无品质参数分布,说明不同成熟度果实在采收初期(0 d)和贮藏末期(180 d)品质总体差异较小。成熟度Ⅰ的 30 d与0 d聚在一起,成熟度Ⅱ和Ⅲ的30 d则均分布于第2象限,说明成熟度Ⅰ果实贮藏30 d与贮藏0 d品质相近,贮藏30 d时成熟度Ⅱ和Ⅲ果实品质间差异较小。成熟度Ⅰ、Ⅱ和Ⅲ果实在贮藏60~150 d 各自聚为一簇(如图1绿色、红色和黄色虚线范围所示),且均与30 d距离较远(成熟度Ⅲ除外);值得注意的是,在60~150 d内,成熟度Ⅱ和Ⅲ两簇截然分离,两者品质差异主要表现在PC1上;而成熟度Ⅰ则与Ⅱ和Ⅲ各有交叉,与其他2个成熟度果实的差异主要体现在PC2上,结合品质指标分布可知,成熟度Ⅰ果实拥有较低的可滴定酸和固酸比。总体而言,采收成熟度Ⅰ表现出较低的PC1和PC2;采收成熟度Ⅱ表现出较低的PC1和较高的PC2;采收成熟度Ⅲ表现出较高的PC1和较低的PC2。结合每个主成分的代表性品质指标及其生物学意义可知,成熟度Ⅱ果实拥有最好的糖酸风味和色泽,但腐烂率和失重率也相对较高。
图1 不同采收成熟度桃溪蜜柚品质参数
和贮藏期PCA聚类图
Fig.1 PCA clusters of quality parameters and storage periods
of Taoxi Pomelo with different harvesting maturities
2.2.3 基于PCA比较3个采收成熟度果实贮藏品质变化差异
为分析不同采收成熟度的桃溪蜜柚果实在贮藏周期内品质变化趋势,分别将采收成熟度Ⅰ、Ⅱ和Ⅲ的品质参数进行PCA聚类,结果如图2所示。每个采收成熟度果实在不同贮藏期品质性状的总差异主要集中在前2个维度主成分1(PC1)和主成分2(PC2)上,其中采收成熟度Ⅰ的品质参数降维后PC1(45.12%)和PC2(20.91%)代表了总差异的66.03%;采收成熟度Ⅱ的品质参数降维后PC1(49.87%)和PC2(29.23%)代表了总差异的79.10%;采收成熟度Ⅲ的品质参数降维后PC1(52.74%)和PC2(23.36%)代表了总差异的76.10%。
3个采收成熟度果实的贮藏品质差异在PC1上(特征值差异幅度>0.7)主要体现为参数:L*、总糖、TA、固酸比、抗坏血酸的差异;在PC2上(特征值差异幅度>1.2)主要体现为参数:L*、TA、TSS、固酸比、抗坏血酸、总酚的差异;总黄酮、PPO差异较小;而腐烂率、失重率、a*、b*、C*、Ho、CCI 和MDA含量变化规律相对一致。这些结果表明采收成熟度不同,果实光泽度、糖酸风味和抗氧化能力有较大差异,但着色变化差异相对较小。
对3个采收成熟度果实贮藏期综合表现进行聚类发现:成熟度Ⅰ和Ⅱ果实品质变化模式相近,7个贮藏时间点均可聚成3簇(图2-a和2-b):0~30 d、60~90 d 和 120~150 d,表明成熟度Ⅰ和Ⅱ果实贮藏品质变化模式均可分为4个阶段:贮藏初期(0~30 d)、贮藏中期(60~90 d)、贮藏晚期(120~150 d)、贮藏终期(180 d)。成熟度Ⅲ 的7个贮藏时间点可聚成两大簇(图2-c):30~90 d、120~150 d,而0 d和180 d分别游离于两簇之外,表明成熟度Ⅲ果实贮藏品质变化模式可分为0、30~60、90~150和180 d 4个阶段。对成熟度Ⅰ而言,抗坏血酸、总糖:分布在60~90 d附近;TA、TSS、L*、总酚4个指标分布在120~150 d附近。对成熟度Ⅱ而言:指标Ho和PPO分布在0~30 d附近;指标TA、抗坏血酸、TSS、L*分布在60~90 d附近;指标总糖、总酚分布在120~150 d附近。对成熟度Ⅲ而言:指标L*、TA、抗坏血酸、总糖、TSS分布在30~90 d附近;总酚分布在120~150 d附近。固酸比指标游离于各个贮藏期外。以上结果说明:TSS、TA、L*变化由成熟度Ⅰ的120~150 d提前至成熟度Ⅱ和Ⅲ的60~90 d。总糖变化在成熟度Ⅱ中推迟至120~150 d。3个采收成熟度果实的抗坏血酸变化主要发生在60~90 d内,总酚变化主要发生在120~150 d。
a-成熟度Ⅰ;a-成熟度Ⅱ;c-成熟度Ⅲ
图2 不同采收成熟度桃溪蜜柚品质参数和贮藏周期PCA聚类图
Fig.2 PCA cluster diagram of the quality parameters and storage periods of Taoxi pomelo with different maturities
为直观地展示3个采收成熟度的桃溪蜜柚贮藏性能的优劣,根据PCA得到的数据计算3个采收成熟度果实品质的综合得分。首先提取PCA特征值>1的前5个主成分作为数据分析的有效成分,其累计方差贡献率达到90.72%,基本反映了所有变量的初始信息。然后根据17个品质指标的荷载系数及特征值计算得到17个因子的特征向量(表5)。以特征向量为系数构建前5个主成分的线性方程,如下所示:
Y1=0.342X1+0.249X2+0.070X3+0.335X4+…+0.313X14-0.146X15+0.247X16+0.217X17
Y2=-0.058X1-0.059X2+0.333X3+0.111X4+…-0.112X14-0.160X15-0.005X16-0.116X17
Y3=0.046X1-0.082X2+0.524X3+0.143X4+…+0.018X14+0.444X15-0.044X16-0.183X17
Y4=0.046X1-0.051X2-0.220X3-0.174X4+…+0.090X14+0.200X15+0.473X16+0.333X17
Y5=-0.052X1+0.565X2-0.123X3-0.001X4+…+0.117X14-0.381X15-0.174X16-0.460X17
将前5个主成分的方差贡献率(表4)作为系数代入,得到综合得分评价方程为:Y=0.454 59Y1+0.201 21Y2+0.094 79Y3+0.087 28Y4+0.069 33 Y5。然后将不同采收成熟度果实贮藏期的品质指标测定值标准化后代入上述方程,计算不同采收成熟度桃溪蜜柚在每个贮藏周期内的品质得分Y,Y值越大品质越好。3个成熟度果实品质得分变化趋势如图3所示。
图3 不同成熟度桃溪蜜柚在常温贮藏期间品质综合得分
Fig.3 Comprehensive quality scores of Taoxi pomelo with
different maturities during ambient storage
3个采收成熟度果实贮藏品质综合得分整体呈上升趋势,说明桃溪蜜柚采后存在后熟过程。除成熟度Ⅲ 150 d外,采收成熟度Ⅱ得分在60~180 d内均显著高于其他2个时期,表明成熟度Ⅱ果实的贮藏表现总体最优。0~30 d内,三者综合得分均为负值;贮藏60 d时,成熟度Ⅱ的综合得分达到正值(0.77);90 d时,成熟度Ⅰ和Ⅲ的综合得分达到正值,此时Ⅰ、Ⅱ和Ⅲ得分分别为1.28、0.34和0.44,随后略有增加,说明贮藏90 d左右品质已达到最佳。根据3个采收成熟度的综合得分可知:桃溪蜜柚采后存在后熟现象,盛花后192 d(成熟度Ⅱ)为桃溪蜜柚的最佳采收成熟度,该时期采收的果实在贮藏90 d时品质达最佳,但考虑到长时间贮藏腐烂率不断增加,因此宜在贮藏90 d后及时出库销售。
表5 主成分的荷载矩阵和特征向量
Table 5 Principal loading matrix and component eigenvectors
指标名称主成分1主成分2主成分3主成分4主成分5荷载系数特征向量荷载系数特征向量荷载系数特征向量荷载系数特征向量荷载系数特征向量失重率0.951 0.342 -0.108 -0.058 0.058 0.046 0.056 0.046 -0.056 -0.052腐烂率0.693 0.249 -0.110 -0.059 -0.104 -0.082 -0.062 -0.051 0.614 0.565L∗0.195 0.070 0.616 0.333 0.665 0.524 -0.268 -0.220 -0.134 -0.123a∗0.930 0.335 0.206 0.111 0.181 0.143 -0.212 -0.174 -0.001 -0.001b∗0.902 0.324 0.047 0.025 -0.236 -0.186 0.022 0.018 -0.094 -0.087C∗0.874 0.314 0.007 0.004 -0.296 -0.233 0.058 0.048 -0.084 -0.077Ho-0.924 -0.332 -0.223 -0.121 -0.201 -0.158 0.210 0.172 0.030 0.028CCI0.897 0.323 0.245 0.132 0.252 0.199 -0.245 -0.201 -0.051 -0.047总糖-0.170 -0.061 0.398 0.215 0.514 0.405 0.632 0.519 0.300 0.276TSS0.328 0.118 0.777 0.420 0.032 0.025 0.387 0.318 0.184 0.169TA-0.063 -0.023 0.931 0.503 -0.255 -0.201 -0.097 -0.080 -0.160 -0.147固酸比0.301 0.108 -0.800 -0.433 0.265 0.209 0.193 0.158 0.329 0.303抗坏血酸-0.465 -0.167 0.645 0.349 -0.316 -0.249 0.298 0.245 0.175 0.161MDA0.871 0.313 -0.207 -0.112 0.023 0.018 0.110 0.090 0.127 0.117PPO-0.407 -0.146 -0.296 -0.160 0.564 0.444 0.244 0.200 -0.414 -0.381总酚0.688 0.247 -0.009 -0.005 -0.056 -0.044 0.576 0.473 -0.189 -0.174总黄酮0.602 0.217 -0.214 -0.116 -0.232 -0.183 0.406 0.333 -0.500 -0.460
随着检测手段的进步,各种果蔬品质指标的测定越来越高效,获得的数据量越来越大,相应地在多变量数据分析中,主成分分析法的应用不可或缺[6,12-15],利用主成分分析可更加简单有效地呈现原多维数据的差异性[8]。
前期李宏祥等采集了盛花后178、192和208 d(成熟度Ⅰ、Ⅱ和Ⅲ)的桃溪蜜柚果实,定期检测了贮藏180 d过程中的19个品质指标,并对这19个品质指标在3种成熟度果实中的差异进行了逐一比较分析,最终认为花后 208 d(采收期Ⅲ)是桃溪蜜柚的最佳采收成熟度,并建议贮存120 d时上市,最晚不要超过150 d[4]。但由于常规ANOVA技术仅能分析不同成熟度果实单个品质指标间的差异,并不能综合表现各个品质指标间的内在关联性和总体差异,这可能导致对桃溪蜜柚的最佳采收成熟度和最佳贮藏期的判断有一定主观偏差。
本文利用主成分分析法,对3个采收成熟度、7个贮藏周期的19个品质参数相关数据重新进行了分析。将19个品质指标相关信息浓缩成5个主成分,代表了总差异的90.72%(表4)。通过前2个主成分(代表总差异的65.58%)直观地展示了3个采收成熟度的桃溪蜜柚贮藏品质的总体差异和其代表性品质指标(图1)。基于品质参数和贮藏周期主成分得分,两两比较不同采收成熟度果实的贮藏品质,展示了不同成熟度果实品质变化的时间模式差异,及其代表性品质指标(图2)。参照金沙柚[5]和金兰柚[3]的PCA综合得分法,利用前5个主成分对不同采收成熟度的桃溪蜜柚贮藏品质进行综合打分,可以明显看出果实品质得分前期呈上升趋势,这与金沙柚和金兰柚结果类似,说明桃溪蜜柚也具有后熟现象在贮藏早期品质逐步改善。但与其他2个井冈蜜柚后期品质得分下降不同的是,桃溪蜜柚的综合得分在90 d后保持稳定,造成这一现象的原因有待进一步探讨。但与李宏祥等[4]结论不同的是:利用PCA综合得分研究发现,成熟度Ⅱ得分在60~180 d内(>1)显著高于成熟度Ⅰ和Ⅲ(<0.7)(图3),据此判断成熟度Ⅱ(盛花后192 d)采收的桃溪蜜柚果实整体品质优于其他2个时期,且在贮藏90 d后达到最佳品质。以PCA综合得分为依据,选择桃溪蜜柚的最适采收期和最佳贮藏周期摒弃了主观偏差,使结论更为客观。
此外,研究表明不同采收成熟度的果实外观、糖酸风味、抗性品质等与耐贮藏性相关的指标有较大差异[9,16-17],适宜的采收期有利于维持果实品质和耐贮藏性[18-21]。本课题组对不同采收成熟度的井冈蜜柚果实贮藏期内品质差异做了相关研究报道[3-5],也发现采收成熟度不同,果实着色、糖酸风味、抗氧化性、腐烂率等贮藏性能差异较大。根据PCA聚类结果可知,成熟度Ⅱ和Ⅲ品质差异明显,在PC1上截然分离;成熟度Ⅰ则介于Ⅱ和Ⅲ之间,分别与Ⅱ和Ⅲ交叉于PC1的不同区域(图1)。与之对应的品质指标聚类结果则进一步展示了不同成熟度桃溪蜜柚果实品质的差异,总体来讲成熟度Ⅰ果实可滴定酸和固酸比较低;成熟度Ⅱ果实糖酸、风味和色泽较好,但腐烂率和失重率也高;与成熟度Ⅱ相反成熟度Ⅲ果实糖酸、风味和色泽较差,但腐烂率和失重率较低(图1)。不同成熟度果实的品质指标在不同贮藏期内的变化模式差异主要体现为:TSS、TA和L*变化由成熟度Ⅰ的120~150 d提前至成熟度Ⅱ和Ⅲ的60~90 d。总糖在成熟度Ⅱ中推迟至120~150 d,这进一步证明了上述观点。
基于主成分分析,桃溪蜜柚存在后熟现象,且采收成熟度不同果实光泽度、糖酸风味和抗氧化能力差异显著。盛花后192 d采收的(成熟度Ⅱ)果实贮藏表现最佳,为桃溪蜜柚的最适采收成熟度,该时期采收的果实在贮藏90 d时品质达最佳,但考虑到更长时间贮藏腐烂增加,因此宜在贮藏90 d后及时出库销售。
[1] 江西省质量技术监督局.DB36/T 810-2014,井冈蜜柚商品果[S].南昌:吉安市果业局,2015.
Jiangxi Provincial Bureau of Quality and Technical Supervision.DB36T 810—2014,Jinggang honey pomelo commodity fruit[S].Nanchang:Ji’an Fruit Industry Bureau,2015.
[2] 王芳飞.吉安市井冈蜜柚产业发展现状及建议[J].现代农业科技,2020,8:94-95.
WANG F F.The status and suggestions on the development of Jinggang honey pomelo industry[J].Modern Agricultural Science and Technology,2020,8:94-95.
[3] 林雄,李宏祥,张亚男,等.主成分分析法综合评价采收成熟度对金兰柚贮藏特性的影响[J].食品与发酵工业,2020,46(9):217-224.
LIN X,LI H X,Z Y N,et al.Effects of harvest maturity on storage properties of Jinlan pomelo based on principal component analysis[J].Food and Fermentation Industries,2020,46(9):217-224.
[4] 李宏祥,马巧利,林雄,等.采收成熟度对桃溪蜜柚贮藏品质及抗氧化性的影响[J].食品与发酵工业,2019,45(13):191-198.
LI H X,MA Q L,LIN X,et al.Effects of harvest maturity on storage quality and antioxidative capacity of Taoxi pomelo[J].Food and Fermentation Industries,2019,45(13):191-198.
[5] 李宏祥,马巧利,林雄,等.PCA综合分析采收成熟度对金沙柚贮藏品质的影响[J].食品工业科技,2019,40(18):255-262;272.
LI H X,MA Q L,LIN X,et al.Comprehensively analyzing the effect of harvest maturity on storage quality of Jinsha pomelo based on PCA[J].Science and Technology of Food Industry,2019,40(18):255-262;272.
[6] JOLLIFFE I T,CADIMA J.Principal component analysis:A review and recent developments[J].Philosophical Transactions of the Royal Society Mathematical Physical & Engineering Sciences,2016,374(2 065):20 150 202.
[7] 李玉,张铭容,董红敏,等.1-MCP对佛手瓜低温贮藏品质影响的主成分分析和综合评价[J].食品与发酵工业,2015,41(11):204-209.
LI Y,ZHANG M R,DONG H M,et al.Principal component analysis and comprehensive assessment on influence of 1-MCP on chayote fruit quality during cold storage[J].Food and Fermentation Industries,2015,41(11):204-209.
[8] ALLEGRA A,GALLOTTA A,CARIMI F,et al.Metabolic profiling and post-harvest behavior of “dottato” fig (Ficus carica L.) fruit covered with an edible coating from O.Ficus-indica[J].Frontiers in Plant Science,2018,9:1321.DOI:10.3389/fpls.2018.01321..
[9] KIM T J,HYEON H,PARK N I,et al.A high-throughput platform for interpretation of metabolite profile data from pepper (Capsicum) fruits of 13 phenotypes associated with different fruit maturity states[J].Food Chemistry,2020,331:127 286.
[10] 高歌,庞雪莉,刘海华,等.基于GC-MS-O香气成分分析和多元统计分析的柚子品种鉴别[J].中国食品学报,2020,20(5):283-292.
GAO G,PANG X L,LIU H H,et al.Volatiles identification of pomelo based on GC-MS-O and multivariate statistical analysis[J].Journal of Chinese Institute of Food Science and Technology,2020,20(5):283-292.
[11] 李俭,周小娟,钟灶发,等.南康不同柚子品种果实品质分析[J].安徽农业科学,2019,47(14):216-219.
LI J,ZHOU X J,ZHONG Z F,et al.Fruit quality analysis of different grapefruit varieties in Nankang[J].Journal of Anhui Agricultural Sciences,2019.47(14):216-219.
[12] TANG N,AN J,DENG W,et al.Metabolic and transcriptional regulatory mechanism associated with postharvest fruit ripening and senescence in cherry tomatoes[J].Postharvest Biology and Technology,2020,168:111 274.
[13] DONADEL J Z,THEWES F R,ANESE R D O,et al.Key volatile compounds of ‘Fuji Kiku’ apples as affected by the storage conditions and shelf life:Correlation between volatile emission by intact fruit and juice extracted from the fruit[J].Food Research International,2019,125:108 625.
[14] TANG N,CHEN N,HU N,et al.Comparative metabolomics and transcriptomic profiling reveal the mechanism of fruit quality deterioration and the resistance of citrus fruit against Penicillium digitatum[J].Postharvest Biology and Technology,2018,145:61-73.
[15] FABRONI S,AMENTA M,TIMPANARO N,et al.Change in taste-altering non-volatile components of blood and common orange fruit during cold storage[J].Food Research International,2020,131:108 916.
[16] RIBEIRO B S,DE FREITAS S T.Maturity stage at harvest and storage temperature to maintain postharvest quality of acerola fruit[J].Scientia Horticulturae,2020,260:108 901.
[17] REICHEL M,CARLE R,SRUAMSIRI P,et al.Influence of harvest maturity on quality and shelf-life of litchi fruit (Litchi chinensis Sonn.)[J].Postharvest Biology and Technology,2010,57(3):162-175.
[18] ECHEVERRA G,CANTN C,M,ORTIZ A,et al.The impact of maturity,storage temperature and storage duration on sensory quality and consumer satisfaction of ‘Big Top®’ nectarines[J].Scientia Horticulturae,2015,190:179-186.
[19] SUN Y D,SINGH Z,TOKALA V Y,et al.Harvest maturity stage and cold storage period influence lemon fruit quality[J].Scientia Horticulturae,2019,249:322-328.
[20] CHEN S W,HSU M,FANG H,et al.Effect of harvest season,maturity and storage temperature on storability of carambola ‘Honglong’ fruit[J].Scientia Horticulturae,2017,220:42-51.
[21] HORVITZ S,CHANAGUANO D,AROZARENA I,et al.Andean blackberries (Rubus glaucus Benth) quality as affected by harvest maturity and storage conditions[J].Scientia Horticulturae,2017,226:293-301.