本发明涉及水产学领域，更具体地说，涉及一种利用肌肉必需脂肪酸组成区分多个地理群体茎柔鱼的方法。
背景技术：
脂肪酸是一种特殊的化合物，主要以三羧酸甘油酯和磷脂的形式存在，其在生物的摄食活动过程中相对稳定不易变化，可用于分析食物来源，是天然生物标志物之一(许强,杨红生.脂肪酸标志物在海洋生态系统营养关系研究中的应用[j].海洋学报,2011,33(1):1-6)。作为生物标志物，脂肪酸有几大优越性：第一，脂肪酸可以通过形成贮存脂，储存在生物体内，从而可反映一段时期内生物的摄食情况，相对于传统的胃含物分析法，减小了判断生物食性的偶然性；第二，在新陈代谢过程中相对稳定，在摄食和同化过程中结构基本保持不变；第三，生物体内三羧酸甘油酯中的脂肪酸主要来自于所摄入的食物(hendersonrj,sargentjr.lipidbiosynthesisinrainbowtrout,salmogairdnerii,feddietsofdifferinglipidcontent[j].comparativebiochemistry&physiologypartccomparativepharmacology,1981,69(1):31-37)。
海洋生物的脂肪酸分成饱和脂肪酸和不饱和脂肪酸两类。其中n-3和n-6系列高度不饱和脂肪酸只能从饵料中获取，不能自身合成，因此这些脂肪酸被称为必需脂肪酸(olseny.lipidsandessentialfattyacidsinaquaticfoodwebs:whatcanfreshwaterecologistslearnfrommariculture？[m]//lipidsinfreshwaterecosystems.springernewyork,1999:161-202)，这些必需脂肪酸常作为海洋食物关系研究中的脂肪酸标志物(gordonbj,tocherdr,macdonaldfm,etal.effectsofdietsrichinlinoleic(18:2n-6)andα-linolenic(18:3n-3)acidsonthegrowth,lipidclassandfattyacidcompositionsandeicosanoidproductioninjuvenileturbot(scophthalmusmaximusl.)[j].fishphysiologyandbiochemistry,1994,13(2):105；smithgi,athertonp,reedsdn,etal.dietaryomega-3fattyacidsupplementationincreasestherateofmuscleproteinsynthesisinolderadults:arandomizedcontrolledtrial[j].americanjournalofclinicalnutrition,2011,93(2):402-412)。因此可利用脂肪酸标志物法了解茎柔鱼不同地理种群食性和营养地位。
目前利用脂肪酸标志物对茎柔鱼不同地理种群进行摄食生态学研究的文献还未见报道。国内研究主要是对其肌肉和腕足的营养成分分析(杨宪时,王丽丽,李学英等,.秘鲁鱿鱼和日本海鱿鱼营养成分分析与评价[j].现代食品科技,2013,29(9):2247-2251)和食品加工研究(王丽丽.鱿鱼足即食休闲风味食品研究初探[d].大连海洋大学,2014)，国外仅saito等(saitoh,sakaim,wakabayashit.characteristicsofthelipidandfattyacidcompositionsofthehumboldtsquid,dosidicusgigas:thetrophicrelationshipbetweenthesquidanditsprey[j].europeanjournaloflipidscienceandtechnology,2014,116(3):360-366)对秘鲁外海茎柔鱼不同组织和胃含物的脂肪酸组成进行了研究。大量研究发现，利用脂肪酸标志物的示踪作用可以较准确地揭示特定类群在单一或不同生境中的食物来源差异以及某一特定复杂生境中同种类群生物的食物来源情况。
技术实现要素：
本发明的目的在于提供一种利用肌肉必需脂肪酸组成区分多个地理群体茎柔鱼的方法，该方法具有快速、准确性高的优点。
为解决上述问题，本发明采用如下的技术方案：
一种利用肌肉必需脂肪酸组成区分多个地理群体茎柔鱼的方法，其包括如下步骤：
(1)提取茎柔鱼肌肉脂肪酸，得到含有茎柔鱼肌肉脂肪酸的溶液；
(2)利用gc-ms测定肌肉脂肪酸；
(3)分析肌肉脂肪酸数据：
a.定量分析肌肉脂肪酸含量，获得样品中n种脂肪酸的最终浓度(a1……an)，根据公式：an×稀释倍数×原始溶液体积，计算n种脂肪酸的含量，其中原始溶液体积为步骤(1)的含有茎柔鱼肌肉脂肪酸的溶液体积；
b.各脂肪酸的含量除以总脂肪酸含量，计算出各脂肪酸所占含量百分比；
c.筛选出p种n-3和n-6系列不饱和脂肪酸含量百分比(w1……wp)作为该样品的必需脂肪酸分析结果；
(4)计算多个地理群体茎柔鱼的n-3和n-6系列不饱和脂肪酸判别函数系数和常量：
a.以地理群体作为分组变量，n-3和n-6系列不饱和脂肪酸含量百分比作为自变量，对已知来自多个(1……m)地理群体茎柔鱼的p种n-3和n-6系列不饱和脂肪酸含量百分比进行线性判别分析；
b.分析得出多个(1……m)地理群体的p种n-3和n-6系列不饱和脂肪酸的相应判别系数(k1-1……k1-p，k2-1……k2-p，km-1……km-p)和m个地理群体的常量(c1……cm)；
其中n，p，m都为整数，
(5)判别：利用上述步骤(1)-(3)得到待判别的茎柔鱼肌肉中p种n-3和n-6系列不饱和脂肪酸含量百分比(w1……wp)，根据步骤(4)的判别系数和常量计算m个判别函数值(f1……fm)，
其中函数值f1、f2……fm的计算方法为：
f1＝k1-1×w1-1+k1-2×w1-2+……+k1-n×w1-p+c1，
f2＝k2-1×w2-1+k2-2×w2-2+…+k2-p×w2-p+c2，
……
fm＝km-1×wm-1+km-2×wm-2+……+km-p×wm-p+cm，
计算得到m个函数值f1……fm，以最大函数值对应的地理群体表示其所属地理群体(也就是说如果f1值最大，则判断该茎柔鱼是第1个地理群体；如果f2值最大，则判断该茎柔鱼是第1个地理群体，……，如果fm值最大，则判断该茎柔鱼是第m个地理群体。)。
在本发明的一个优选实施方案中，步骤(1)提取茎柔鱼肌肉脂肪酸包括以下步骤：
i.取样：漏斗锁软骨处剪取肌肉，用超纯水清洗样品表面，再去除皮肤，放入离心管中冷冻保存；
ii.冷冻干燥：将冷冻后的样品放入冷冻干燥机中，-50℃至-60℃冷冻干燥24-48小时；
iii.研磨：干燥后的样品研磨成粉末；
iv.提取脂肪酸：利用三氯甲烷-甲醇(体积比1.5-2.5:1)浸泡上述粉末20-40小时，离心，上清液用0.6-1.5％的氯化钠水溶液洗涤，静置1.5-2.5小时后分液，取下层溶液并除去其中的溶剂，得到的脂肪酸在0.3-0.6mol/l氢氧化钠-甲醇溶液中回流8-12min，再加入三氯化硼-甲醇溶液(10-20％三氯化硼的甲醇溶液)回流25-35min，最后加入正己烷回流1-3min，冷却后加入饱和氯化钠水溶液，混匀后静置20-30min，最后得到含有脂肪酸的正己烷上层。
在本发明的一个更优选实施方案中，步骤ii中，将冷冻后的样品放入冷冻干燥机中，-50℃冷冻干燥。
在本发明的一个更优选实施方案中，步骤iii中，干燥后的样品使用冷冻混合球磨仪mm400(retsch，德国)
在本发明的一个更优选实施方案中，步骤iv中，所用三氯甲烷-甲醇中三氯甲烷和甲醇的体积比为2:1。在本发明的一个更优选实施方案中，氯化钠水溶液的质量浓度为0.9％。在本发明的一个更优选实施方案中，步骤iv中，氢氧化钠-甲醇溶液为0.5mol/l氢氧化钠的甲醇溶液。在本发明的一个更优选实施方案中，步骤iv中，三氯化硼-甲醇溶液为14％三氯化硼的甲醇溶液。
在本发明的一个更优选实施方案中，步骤iv中，上清液中加入0.9％的氯化钠水溶液静置2小时后分液。在本发明的一个更优选实施方案中，步骤iv中，得到的脂肪酸在氢氧化钠-甲醇溶液(0.5mol/l氢氧化钠的甲醇溶液)中回流10min。在本发明的一个更优选实施方案中，步骤iv中，加入三氯化硼-甲醇溶液(14％三氯化硼的甲醇溶液)回流30min。在本发明的一个更优选实施方案中，步骤iv中，加入正己烷回流2min。
在本发明的一个更优选实施方案中，步骤iv中氢氧化钠-甲醇溶液、三氯化硼-甲醇溶液和正己烷的用量，按体积比计，为1:1:1。
在本发明的一个更优选实施方案中，步骤(2)gc-ms测定脂肪酸利用含有脂肪酸的正己烷溶液，其中加入十九烷酸甲酯内标。
在本发明的一个更优选实施方案中，十九烷酸甲酯与正己烷溶液的体积比是1:1。
在本发明的一个优选实施方案中，步骤(2)gc-ms测定脂肪酸利用7890b/5977a气质联用仪(agilent，美国)测定，色谱柱为hp-88(60m×0.25mm×0.20μm，agilent，美国)，测样升温程序为初始125℃，以8℃/min升到145℃，保持26min，以2℃/min升到220℃，保持1min，再以1℃/min升到227℃，保持1min，进样口温度250℃，辅助加热器温度280℃，分流比10：1，载气为氦气。
在本发明的一个优选实施方案中，步骤(3)a采用masshunterworkstation根据步骤(2)gc-ms测定的肌肉脂肪酸定量分析脂肪酸含量。
在本发明的一个优选实施方案中，步骤(4)使用spss软件对p种n-3和n-6系列不饱和脂肪酸含量百分比进行线性判别分析。在本发明的一个更优选实施方案中，“使用spss软件对p种n-3和n-6系列不饱和脂肪酸含量百分比进行线性判别分析”是通过将多个地理群体的n-3和n-6系列不饱和脂肪酸含量数据导入spss软件，选择“分析”工具中的“判别”，以地理群体作为分组变量，p种n-3和n-6系列不饱和脂肪酸含量(w1……wp)作为自变量，“统计量”选项中的“函数系数”选择“fisher”,“分类”选项中的“输出”模块选择“摘要表”，“图”模块选择“合并组”，其他均为默认设置，得到多个地理群体的p种n-3和n-6系列不饱和脂肪酸的对应的判别系数和m个地理群体的常量。
在本发明的一个优选实施方案中，n-3和n-6系列不饱和脂肪酸为c18:2n6c、c18:2n6t、c18:3n3、c18:3n6、c20:3n3、c20:3n6、c20:4n6、c20:5n3、c22:2n6和c22:6n3。
本发明中，不同n-3和n-6系列不饱和脂肪酸的判别系数是不同的，并且不同地理群体的同一种n-3和n-6系列不饱和脂肪酸也是不同的。本发明中“判别系数(k1-1……k1-p，k2-1……k2-p，km-1……km-p)”表示不同地理群体的不同n-3和n-6系列不饱和脂肪酸的判别系数。例如，“k1-1”表示m个地理群体中的第1个地理群体茎柔鱼的p种n-3和n-6系列不饱和脂肪酸中的第1种不饱和脂肪酸的判别系数，“km-p”表示m个地理群体中的第m个地理群体茎柔鱼的p种n-3和n-6系列不饱和脂肪酸中的第p种不饱和脂肪酸的判别系数。
本发明中，“常量”对一个地理群体来说是不变的，m种地理群体具有m个常量。本发明中“常量(c1……cm)”表示m个地理群体的m个常量。例如，“c1”表示m个地理群体中的第1个地理群体的常量，“cm”表示m个地理群体中的第m个地理群体的常量。
本发明中，“判别函数值(f1……fm)”是由多种n-3和n-6系列不饱和脂肪酸含量百分比乘以相应地理群体的判别系数得到的和再加上该地理群体的常量得到的，例如“f1”是由多种(p种)n-3和n-6系列不饱和脂肪酸含量百分比乘以相应第1个地理群体的判别系数得到的和再加上第1个地理群体的的常量得到的。“fm”是由多种(p种)n-3和n-6系列不饱和脂肪酸含量百分比乘以相应第m个地理群体的判别系数得到的和再加上第m个地理群体的的常量得到的。
本发明经过大量不同地理群体的茎柔鱼肌肉必需脂肪酸组成比对，创造了一种简单实用的方法，仅通过测量茎柔鱼肌肉必需脂肪酸组成，就能快速，准确地区分鉴别茎柔鱼的地理群体，准确率可高达96％。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体例，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。下面实施例中未注明具体条件的实验方法，通常按照常规条件，或按照制造厂商所建议的条件。
实施例1
一种利用肌肉必需脂肪酸组成区分茎柔鱼地理群体的方法，包括如下步骤：
随机采集来自3个已知地理区域的66尾茎柔鱼，分别为东太平洋赤道海域19尾(0°28'n～1°00's、115°01'w～119°00'w)、秘鲁外海27尾(9°50's～15°42's、79°45'w～85°03'w)和智利外海20尾(37°06's～38°00's、79°01'w～83°00'w)。
(1)提取茎柔鱼肌肉脂肪酸：
i.取样：由胴体腹面剪开，在漏斗锁软骨处使用手术剪刀剪取2×2cm肌肉，用超纯水清洗样品表面，再使用镊子去除皮肤，放入10ml离心管中冷冻保存；
ii.冷冻干燥：将冷冻后的样品放入冷冻干燥机中，-55℃冷冻干燥30小时；
iii.研磨：干燥后的样品使用冷冻混合球磨仪mm400(retsch，德国)研磨成粉末(如果不立即使用，该粉末可装入离心管中，使用真空压缩袋常温保存)；
iv.提取脂肪酸：
取0.2g粉末于15ml离心管中，
加12ml三氯甲烷-甲醇溶液(体积比2:1)，浸泡30小时，
离心后取上清液10ml，加入4ml的0.9％的氯化钠水溶液，静置2小时，
取下层溶液于圆底烧瓶中水浴加热得到脂肪酸，
在圆底烧瓶中加入4ml0.5mol/l的氢氧化钠-甲醇溶液，混合后连接水浴回流装置，水浴加热10min，
加入4ml三氯化硼甲醇溶液(14％三氯化硼的甲醇溶液)水浴回流30min，
加入4ml正己烷水浴回流萃取2min，
冷却后加入10ml氯化钠饱和水溶液，摇晃，再继续加入氯化钠饱和水溶液至烧瓶颈部，静置分层30min，分离得到含有脂肪酸的正己烷层(上层)4ml；
(2)利用gc-ms测定肌肉脂肪酸：
取1ml上层溶液用针式过滤器过滤置于进样瓶中，再加入1ml十九烷酸甲酯溶液，混合后进行gc-ms测定，利用7890b/5977a气质联用仪(agilent，美国)，色谱柱为hp-88(60m×0.25mm×0.20μm，agilent，美国)，测样升温程序为初始125℃，以8℃/min升到145℃，保持26min，以2℃/min升到220℃，保持1min，再以1℃/min升到227℃，保持1min，进样口温度250℃，辅助加热器温度280℃，分流比10:1，载气为氦气；
(3)分析肌肉脂肪酸数据
a.采用masshunterworkstation定量分析脂肪酸含量，获得样品中n种脂肪酸的最终浓度(a1……an)，根据公式an(mg/l)×稀释倍数×0.004l计算n种脂肪酸的含量；
b.各脂肪酸的含量除以总脂肪酸含量，计算出各脂肪酸所占含量百分比(表1)；
c.筛选出10种n-3和n-6系列不饱和脂肪酸(c18:2n6c、c18:2n6t、c18:3n3、c18:3n6、c20:3n3、c20:3n6、c20:4n6、c20:5n3、c22:2n6和c22:6n3)的含量百分比作为该样品的必需脂肪酸分析结果；
表1茎柔鱼脂肪酸数据，其中n.d.表示未检测发现。
(4)计算多个地理群体的n-3和n-6系列不饱和脂肪酸的判别函数系数和常量：
a.以地理群体作为分组变量，10种n-3和n-6系列不饱和脂肪酸含量(w1……w10)作为自变量，使用spss软件将已知来自3个地理群体的66尾茎柔鱼的10种n-3和n-6系列不饱和脂肪酸含量百分比(w1-1……w1-10，w2-1……w2-10，w3-1……w3-10)进行线性判别分析；
进行线性判别分析的具体操作过程为：将3个地理群体的10种n-3和n-6系列不饱和脂肪酸含量数据导入spss软件，选择“分析”工具中的“判别”，以地理群体作为分组变量，10种n-3和n-6系列不饱和脂肪酸含量(w1……w10)作为自变量，“统计量”选项中的“函数系数”选择“fisher”，“分类”选项中的“输出”模块选择“摘要表”，“图”模块选择“合并组”，其他均为默认设置；
b.分析得出3个地理群体茎柔鱼的10种n-3和n-6系列不饱和脂肪酸的判别系数(k1-1……k1-10，k2-1……k2-10，k3-1……k3-10)和常量(c1，c2，c3)，
其中东太平洋赤道海域的k1-1至k1-10分别为24.52、40.85、-23.51、15.89、12.29、-1.13、30.32、13.59、8.96、10.23)，常量c1为-416.35；秘鲁外海的k2-1至k2-10分别为28.20、43.56、-6.51、14.06、6.16、9.51、19.12、11.17、4.88、9.67，常量c2为-342.23；智利外海的k3-1至k3-10分别为14.59、68.39、-17.22、24.31、-11.64、6.49、25.86、11.52、4.96、10.96，常量c3为-439.81；
3种地理群体茎柔鱼的10种n-3和n-6系列不饱和脂肪酸的判别系数(k1-1至k1-10，k2-1至k2-10，k3-1至k3-10)和常量(c1，c2，c3)参见表2；
(5)为同时证明本发明的鉴别准确率，抓取24尾茎柔鱼(东太平洋赤道海域8尾，秘鲁外海8尾，智利外海8尾)，分别标记其实际地理群体后，再运用本发明的方法进行鉴别区分。利用上述步骤(1)-(3)得到24尾待判别的茎柔鱼肌肉n-3和n-6系列不饱和脂肪酸含量百分比(见表3)，利用每尾茎柔鱼n-3和n-6系列不饱和脂肪酸含量百分比和3个地理群体茎柔鱼的判别系数和常量计算得到3个函数值f1、f2和f3，以最大函数值对应的地理群体表示其所属地理群体。也就是说，如果f1值最大，则判断该茎柔鱼是东太平洋赤道海域群体；如果f2值最大，则判断该茎柔鱼是秘鲁外海群体；如果f3值最大，则判断该茎柔鱼是智利外海群体。
例如编号为1茎柔鱼的三个函数值的计算方法为：
函数值f1＝1.37(w1-1)×24.52(k1-1)+1.76(w1-2)×40.85(k1-2)+0.39(w1-3)×(-23.51)(k1-3)+0.71(w1-4)×15.89(k1-4)+1.73(w1-5)×12.29(k1-5)+0.00(w1-6)×(-1.13)(k1-6)+4.88(w1-7)×30.32(k1-7)+8.11(w1-8)×13.59(k1-8)+2.04(w1-9)×8.96(k1-9)+39.60(w1-10)×10.23(k1-10)+(-416.35)(c1)＝394.08
函数值f2＝1.37(w2-1)×28.20(k2-1)+1.76(w2-2)×43.56(k2-2)+0.39(w2-3)×(-6.51)(k2-3)+0.71(w2-4)×14.06(k2-4)+1.73(w2-5)×6.16(k2-5)+0.00(w2-6)×9.51(k2-6)+4.88(w2-7)×19.12(k2-7)+8.11(w2-8)×11.17(k2-8)+2.04(w2-9)×4.88(k2-9)+39.60(w2-10)×9.67(k2-10)+(-342.23)(c2)＝367.95
函数值f3＝1.37(w3-1)×14.59(k3-1)+1.76(w3-2)×68.39(k3-2)+0.39(w3-3)×(-17.22)(k3-3)+0.71(w3-4)×24.31(k3-4)+1.73(w3-5)×(-11.64)(k3-5)+0.00(w3-6)×6.49(k3-6)+4.88(w3-7)×25.86(k3-7)+8.11(w3-8)×11.52(k3-8)+2.04(w3-9)×4.96(k3-9)+39.60(w3-10)×10.96(k3-10)+(-439.81)(c3)＝354.71
对比编号为1茎柔鱼3个判别函数值可知，函数值f1(394.08)>函数值f2(367.95)>函数值f3(354.71)，即1号茎柔鱼属于东太平洋赤道海域。
表2n-3和n-6系列不饱和脂肪酸、相应地理群体的判别系数和常量
从上述结果看，24只茎柔鱼，有23只茎柔鱼所属地理区域判别正确，判别准确率达96％。可见本发明的方法可以为渔业生产和管理机构提供快速、可靠的鉴别方法。所以，渔业生产者在销售或相关管理部门检测茎柔鱼时，只需要从大批量的茎柔鱼中测量其必需脂肪酸含量，即可判别其所属的地理群体。
已经根据本发明的优选生产形态通过实例对本发明进行了没有限制的描述，但是应当理解，在所附权利要求书定义的范围内，专家可以做变化和/或变型，而不脱离相关的保护范围。