油茶是海南岛的传统植物资源，油茶油又被海南人特称为山柚油，海南特有的气候环境和地理特点孕育了丰富而独具特色的油茶资源，所产的油茶油的口感和品质也有别于内地^[1]，海南种植的油茶树种主要是分布于热带和南亚热带地区的特产油茶树种——越南油茶(Camellia vietnamensis T. C. Huang ex Hu)^[2]，该种同时别称为高州油茶、陆川油茶、华南油茶和大果油茶等。海南省现存结果的油茶老林均是越南油茶实生林，新垦成年林也主要以越南油茶实生树为主。

油茶油的产量和品质，决定了油茶在种植和生产利用上的价值和方向^[3]。海南油茶产地遍布全岛，各产区纷纷打造不同油茶油商品品牌，具有不同的产品质量标准，前人探讨过海南个别产区油茶油的农艺性状、经济性状、油茶油理化性质和脂肪酸组成表现特点，为新品种选育提供了一定的理论依据和方法^[4-5]，然而，不同产区的越南油茶在农艺性状和经济性状上究竟有无差异或有怎样的差异，则未见报道。本文对8个新老产区的越南油茶实生林单株样本开展了农艺性状和经济性状的比较研究，以期了解不同产区的越南油茶栽培特点。

在海南省8个越南油茶产区选定正常开花结果单株作为试验材料，除儋州市王五镇埋牛村3株(MN2、MN5、MN6，字母和数字组成的代号表示单株编号，下同)外，其它各产区分别选5株，这些产区包括：保梅岭生态保护区(JB6、JB11、JB18、JB21、JB26)、霸王岭长臂猿生态保护区(QS5、QS6、QS7、QS12、QS15)、白沙县细水乡江排村(XS1、XS7、XS8、XS9、XS10)、澄迈县大丰镇铭缘农业有限公司(MY4、MY8、MY10、MY11、MY13)、澄迈县中兴镇横滩村(HT3、HT4、HT5、HT7、HT8)、澄迈县福山镇长安村(HC2、HC3、HC5、HC7、HC8)和琼海市阳江镇上科村(HD2、HD3、HD4、HD5、HD8)等。各地环境因子和树龄等主要信息如表 1所示。

在寒露至霜降期间先后从琼西南、西北和东取样，在各试验单株树冠中部外围四方随机选当年生春梢30根，采集其中间的无病虫害成熟叶片1片，共计30片；在各试验单株树冠中部外围四方随机采无病虫害的果实30个；按照不同单株分装叶片和果实样品并带回实验室。

叶片长度和宽度使用直尺测量，叶形指数=叶长/叶宽。果实纵径、横径、果皮厚度均采用分辨率为0.01 mm的电子数显游标卡尺测量，果形指数=纵径/横径；采用精确度为0.01 g的电子秤称每个单株果实样品重量，计算出平均单果重。

以石油醚作提取剂，采用索氏提取法提取油脂。

果实充分后熟开裂后，取出种子，采用电子秤称每个单株全部种子样品重量，计算出鲜果出籽率和种子百粒重。将种子放入烘箱在110 ℃下杀青15 min，然后在70 ℃下烘至恒定质量并称种子干重，取出种仁，采用电子秤称种仁干重并计算出干种子出仁率；将干种仁研磨粉碎后，采用1/100天平称取2~5 g种仁干粉，提取油脂，采用1/100天平称重后计算出干种仁含油率；最后换算出干种子含油率、鲜种子含油率和果实含油率。另外，单位树冠面积产量=单株果实产量/树冠垂直投影面积。

过氧化值的测定参照GB/T 5538—2005《动植物油脂过氧化值的测定》^[6]执行，碘值的测定参照GB/T 5532—2008《动植物油脂碘值的测定》^[7]执行，酸价的测定参照GB/T 5530—2005《动植物油脂酸值和酸度测定》^[8]执行，皂化值的测定参照GB/T 5534—1995《动植物油脂皂化值的测定》^[9]执行。

按照国家标准GB/T 17376—2008《动植物油脂脂肪酸甲酯制备》^[10]的要求制备脂肪酸甲酯。采用文献[11]的气相色谱法检测脂肪酸种类及其含量。

采用SAS软件GLM过程作方差分析，用DUNCAN法对不同产区作多重比较分析；去掉差异不显著的指标，对全部单株采用CLUSTER过程和最小距离法作多元聚类分析。

由表 2可知，不同产区的叶片指数和果形指数都无显著差异，其余农艺性状均差异显著。各产区叶片长度以XS最长和以MY最短，XS和MY差异显著，QS、HD与MY差异显著而与XS差异不显著，HT、MN与XS差异显著而与MY差异不显著，JB、HC则与其余所有6个样品无显著差异而居中。叶片宽度以QS、XS、HD最大，以MY最小，其它4个产区居中且与最大、最小两组无显著差异。平均单果重以HT最大和以MN、HD最小，其它5个产区居中且与最大、最小两组差异不显著。果实纵径以JB、XS最大和以MN、HD果实最小，其它4个产区居中且与最大、最小两组均无显著差异；果实横径以HT最大和以HD最小，其它6个产区居中且与最大、最小两组均无显著差异。果皮以XS、HT最厚和以MN、HD最薄，其它4个产区居中且与最大、最小两组均无显著差异。可见，不同产区不影响叶片和果实形状，但影响叶片和果实大小、果皮厚度。即说明叶片和果实形状可能不易受环境影响，应能反映遗传性差异；环境通过影响叶片和果实大小、果皮厚度而影响产量和品质等经济性状。

由表 3可知，不同产区的丰产单株单位树冠面积果实产量、种子数、百粒重等产量性状差异显著。单位树冠面积产量以HD和MY较高，HD与MY差异不显著，其它6个产区均显著低于HD；XS、HC和MN差异不显著且最低；QS、JB和HT均与MY和XS、HC、MN无显著差异。种子数以MN和QS较高，且二者无显著差异，其它6个产区显著低于MN；QS同时显著高于XS、MY、HT、HD和HC，与JB无显著差异；JB与XS、MY、HT和HD无显著差异，而显著高于HC；HC也显著低于XS、MY，而也与HT、HD无显著差异。百粒重以HC和HT较大，二者差异不显著，其余6个产区则显著低于HC且相互间无显著差异；HT与JB、XS、MY和HD无显著差异，但显著高于QS和MN。不同产区的其它产量性状如出籽率、干种子出仁率、干种子含油率、鲜种子含油率、果实含油率等无显著差异。可见，果实的出籽率、干种子出仁率、干种子含油率、鲜种子含油率、果实含油率遗传稳定性好，环境通过影响单位树冠面积果实产量和种子多少、重量而影响最终的油脂产量。

由表 4可知，不同产区油茶油的酸价无显著差异，都符合《油茶籽油》(GB/T11765—2018)所规定的一级油标准，说明各产区油茶油的新鲜度和耐贮藏能力都很好。不同产区油茶油过氧化值差异显著，其中MN最低，XS次之，其它6个产区过氧化值较高，与前两组差异显著，但是都远低于《油茶籽油》(GB/T11765—2018)所规定的一级油阈值标准，说明不同产区油茶油抗氧化能力都很强，具有很高的保健价值和良好的耐贮性，同时在产区之间基于优质油基础之上还可以进一步区分优劣。不同产区油茶油碘值差异显著；HC、MN和MY显著高于XS；HC也显著高于JB，但MN、MY与JB无显著差异，JB与XS也无显著差异；QS、HT和HD与其它所有5个产区均无显著差异；说明不同产区油茶油不饱和脂肪酸的不饱和程度不一样，即脂肪酸碳链不饱和键价态或数量存在显著差异。不同产区油茶油皂化值差异显著，XS显著较高，说明该产区油茶油的相对分子质量较高，不易被人体消化吸收，其它各产区间则无明显差异。总之，不同产区显著影响油茶油理化性质(品质)，不同产区油茶油在品质上均达《油茶籽油》(GB/T11765—2018)所规定的一级产品标准，其中XS、HC、MN和HD等产区在某些品质因素上表现特殊，具有不同的食用营养或加工优点。

对不同产区茶籽油依据其4项理化性质作聚类分析，聚类树形图如图 1所示。全部供试单株可聚为7类(R²=0.906)，其中XS10、XS9、XS8、HT7、MN4和XS7各位一类，其余各地剩余单株聚为一类。XS的5个单株各自分属不同的类，对同一产区而言，显示其品质性状遗传多样性丰富。HT和MN在各自产区中具有少数品质性状遗传特异的单株。其余产区如JB、QS、MY、HC和HD则品质性状遗传多样性趋于单一，并且聚类无产区特征，即相对于环境变化呈现出遗传稳定。

	图 1 不同产区茶油理化性质聚类分析树形图 Fig. 1 Tree diagram by cluster analysis on physicochemical properties of tea oil in different producing regions
图选项

从表 5中可以看出，在不同产区油茶油中共检测到9种脂肪酸，不同产区油茶油的主要脂肪酸组成种类包括棕榈酸、硬脂酸、油酸、亚油酸、亚麻酸等，在个别产区油茶油中还发现了少量肉豆蔻酸、棕榈油酸、花生酸和花生一烯酸，所有产区的油茶油主要脂肪酸组分含量均符合《油茶籽油》(GB/T11765—2018)所规定的质量标准。其中油酸含量最高，不同产区均高于75%，且不同产区差异不显著，说明其在不同产区表现相对稳定。棕榈酸是油茶油中相对含量最高的饱和脂肪酸，各产区含量差异不明显，说明其具有较强的遗传稳定性。肉豆蔻酸、亚油酸含量在不同产区间差异不显著。花生一烯酸、花生酸、硬脂酸在不同产区均相对含量少，但在不同产区间差异显著。MN产区花生一烯酸含量较高，其它各产区差异不显著；花生酸仅在MY产区发现；硬脂酸含量则以QS、HT最高，以HC、HD最低，其它4个产区含量居中且与最高、最低两组差异不显著。说明这些比重极小的脂肪酸在反映不同产区的实生单株遗传差异性上很重要，这些脂肪酸也极有可能是造成不同产区油茶油品质差异的重要因素，同时有可能为不同产区的实生单株进行地理上的分类提供依据。

不同单株茶籽油脂肪酸组成的聚类分析树形图如图 2所示。38个样品可被分为4类(R²=0.909)，即：第Ⅰ类包括JB6、JB26、HD3、QS5、QS6、QS15、XS1、XS7、MY8、MY11、HT7、HT8、MN2、MN4，第Ⅱ类包括JB18、JB21、HD2、HD4、HD5、HD8、HC2、HC3、HC5、HC7、HC8、XS8、XS9、XS10、MY4、MY10、HT3、MN5，第Ⅲ类包括QS7、QS12、HT5，第Ⅳ类：JB11、MY13、HT4。其中HC的5个单株聚进同一类，说明该产区油茶油脂肪酸组成遗传多样性可能很少。HD、QS、XS的5个和MN的3个单株聚进两类，说明这两个产区油茶油脂肪酸组成遗传多样性也可能较少。MY和JB的5个单株聚进三类，说明这两个产区油茶油脂肪酸组成遗传多样可能较丰富。HT的5个单株分属4个不同的类，说明该产区油茶油脂肪酸组成遗传多样性丰富。聚类打破了产区分布关系，不同产区的单株聚进同一类，或者绝大多数同产区不同单株聚进不同的类，说明海南油茶实生单株油茶油脂肪酸组成的遗传多样性丰富，并且基因型与环境互作决定油茶油脂肪酸组成。

	图 2 不同产区油茶油脂脂肪酸组成的聚类分析树形图 Fig. 2 Tree diagram by cluster analysis on fatty acid composition in different producing regions
图选项

本文研究表明，不同产区果实农艺性状和经济性状之间不具备简单的线性关系，由遗传因素和环境共同影响，这与前人在高州油茶(即越南油茶)上的研究结果一致^[3]，8个产区的出籽率、干种子出仁率和鲜果含油率等均符合高州油茶优树评价标准^[12]，其中，MN的种子数、出籽率、含油率在各产区中

都处于较高水平，说明其基因型可能较其它产区优良，对海南油茶新品种选育具有一定的指导意义；但是，其单位树冠面积产量和百粒重处于较低水平，这可能是因为该产区缺乏精细的人工管理；HD的单位树冠面积产量和出籽率均是所有产区中最高，但其出油率较低，可能是该产区年平均温度和降水量较高导致，这与前人研究结果一致，温度与降水是影响含油率的主要因子^[13]，对油茶栽培管理方式的改良具有一定的指导意义。本文研究结果还表明皂化值、碘值和过氧化值在不同油茶产区间存在较为丰富的遗传多样性，具备品质遗传改良的可能性。聚类分析结果表明海南油茶实生单株油茶品质遗传表现较为稳定，油脂肪酸组成的遗传多样性丰富，至于基因型与环境互作如何决定油茶油产量和品质性状等的表现则有必要进行更加深入的研究。

前人在对高州油茶种籽油脂肪酸组分分析中共鉴定出了5种脂肪酸^[3]，而本文检测出9种脂肪酸，说明海南产区越南油茶茶籽油营养更丰富，呈现出鲜明的地方特色。各产区的脂肪酸组成和含量存在一定的差异，与前人在普通油茶(Camellia oleifera Abel)上的研究结果基本一致^[14]，其产生原因包括栽培管理方式、立地条件和气候因子等差异^[15]，本文所涉及产区存在管理和环境因子的差异，在越南油茶上可能也支持这一点。MY和HC的脂肪酸组成均较为丰富，综合分析这可能主要得益于其较为精细的人工栽培管理，说明精细的人工管理能够在油茶造林前期就实现优质栽培的目的。MN、XS、JB和QS的脂肪酸组成较为丰富则可能得益于优良的基因型和海南岛西部、西南部充足的光热资源，在失管状况下显示出这样的优良品质，则说明在这些产区可能存在油茶油品质优良的种质资源，对越南油茶新品种选育很有利^[16]。至于基因型与环境互作如何决定油茶油脂肪酸组成则有必要进行更加深入的研究。

海南产越南油茶茶籽油营养丰富，保健价值高，在海南省已逐渐成为了深受消费者推崇的高端食用植物油，本研究结果亦为其提供了可靠的数据支持和理论依据。基于海南不同产区油茶油均达到《油茶籽油》(GB/T11765—2018)所规定的一级产品质量标准，且不同于岛外其它产区的越南油茶(高州油茶)茶籽油，本岛各产区油茶油又有其营养和综合品质上特殊性，因此，我们应该先申请地理产品标识和原产地保护，再进一步挖掘、强化各产区商品品牌内涵，同时加强对外宣传，开发海南油茶油文化价值，提升海南油茶油的知名度，拓展市场空间^[17]，促进海南油茶产业发展。

环境和基因型通过影响越南油茶农艺性状如果实和种子大小等影响其产量性状，叶片和果实形状遗传稳定而不易受环境影响；环境通过影响单位树冠面积果实产量和种子多少、重量而影响最终的油脂产量；海南不同产区的油茶油在营养、保健和贮藏品质等上均符合《油茶籽油》(GB/T11765—2018)的质量要求，同时部分产区又有其特色，聚类分析结果则表明油茶油品质性状遗传较为稳定；海南省油茶油脂肪酸成分种类较其它地区丰富，基因型和环境的互作效应决定脂肪酸的组成和含量，聚类分析结果与该结果一致；海南产区越南油茶新品种选育的育种目标主要在丰产性、稳产性和抗性等的遗传改良上。