1、1尖孢镰刀菌寄主专化型电子舌性脂肪酸判别鉴别模型的建立 刘波*，肖荣 凤， 朱育菁，胡桂萍（福建省农业科学院农业生物资源研究所，福建 福州 350003）摘要：应用基于通过微生物细胞脂肪酸成分鉴定的全自动微生物鉴定系统(Sherlock MIS)，检测分析不同寄主尖孢镰刀菌中脂肪酸的生物标记结构组成和的分布特性。，从供试的 22 个菌株中检测到 8 种脂肪酸，其中常见脂肪酸有 6条，含量 60%以上的主要脂肪酸 脂肪酸 有 16:0、18:0 和 18:2 CIS 9,12/18:0a 在所有的菌株中都存在，各菌株中 3 种脂肪酸之和占总脂肪酸含量的 60%以上，为主要脂肪酸 。对常见脂肪酸进
2、行 通过系统聚类分析，可将 26 个菌株分为 3类。各类群的划分结果显示仅与脂肪酸生物标记在各菌株中的存在情况相关，而常见脂肪酸与菌株的来源、和生化型、致病型进化演变等无有对应关系。对常见脂肪酸进行主成分分析得到尖孢镰刀菌脂肪酸中的特征脂肪酸有 5 个，分别为 14:0（X1），16:0(X2) ，16:1 CIS 9 (W 7)(X3)，18:0(X4) 和 18:2 CIS 9,12/18:0a(X5)。对尖孢镰刀菌特征脂肪酸进行 Bayes 逐步判别分析得到尖孢镰刀菌寄主判别模型为 Y1=-183.06+12.24X1+ 8.05X2+5.39X4+3.58X5，Y2=-174.09+1
3、0.40X1+7.96X2+5.40X4+3.42X5，Y3=-173.97+14.85X1+7.25X2+6.58X4+3.53X5，当 1Yi2 时，菌株的寄主为西瓜； 当 2Yi3 时,菌株的寄主为香蕉；当 Yi3 时，菌株的寄主为花生，且模型判对的概率为 0.90909，说明脂肪酸 14:0， 16:0， 18:0， 18:2 CIS 9,12/18:0a 是判别尖孢镰刀菌寄主的四个重要指标，同时该模型比较符合实际，有很好的应用价值。关键词：尖孢镰刀菌；脂肪酸；寄主；判别模型Abstract:Characteristics of PLFA biomarkers in Fusarium
4、oxysporum Schl. were study by means of microbe identification system (MIDI system from USA), The results showed that 8 fat acids were identified in 26 strains,there were fat acids 16:00、18:00 and 18:2 CIS 9,12/18:0a in all strains. Using the cluster analysis method the 26 strains were classified int
5、o 4 groups based on the PLFAs biomarkers,and had no association with origin, biochemical variant and pathotype of Fusarium oxysporum Schl.Key: Fusarium oxysporum Schl.；fat acid ；host ; the discriminant model 尖孢镰刀菌（Fusarium oxysporum Schl.）引起的农作枯萎病是一种世界性分布的土传真菌病害，寄主范围广泛，可引起100多种植物发病，造成严重减产，甚至绝收。尖孢镰刀菌
6、由致病菌株与无致病菌株组成，种内的致病菌株具有高度的寄主专化性，根据对不同寄主植物的专化性可分为专化型，在同一尖孢镰刀菌不同的的专化型，引起的发病程度完全不同，如尖孢镰刀菌西瓜专化型就不会侵染黄瓜专化型，同样，尖孢镰刀菌古巴专化型只会侵染香蕉而不会侵染花生作物下面。，根据致病力的差异又划分为许多生理小种。尖孢镰刀菌专化型致病型的监测是研究病原菌群体结构的一收稿日期： 接受日期：基金项目：国家高技术研究与发展计划（863 ）项目（2006AA10A211 ）；福建省财政专项- 福建省农业科学院科技创新团队建设基金(STIF-Y03)；福建省发改委重点项目（闽发改投资 2006781 号）。福建省
7、自然科学基金：非致病性尖孢镰刀菌对西瓜枯萎病的免疫抗病机制研究（2008J0055）作者简介：蓝江林（1972-），博士，副研究员，研究方向：生物技术与生物防治。电子邮件： 。电话：0591-87848933。通讯作者：刘波，博士，研究员。研究方向：生物技术与生物防治。电子邮件：。电话：0591-87864601。2项重要措施，为病害的防治提供可靠的病原菌信息，并可为抗病育种提供抗性基因。植物病理学家和育种学家认识到病原真菌的变异对寄主植物抗性的有效性和持久性非常重要，病原真菌的基因型与特定的寄主植物基因型相互作用可导致抗病性在很短时间内丧失【BROWN J K，1996】。然而，鉴定尖孢镰刀
8、菌的专化型，方法复杂、时间很长、实验难度较大。电子舌技术是一种由气相色谱仪和判别模型组成的脂肪酸物质差异分析机构，能识别单一的或混合的脂肪酸溶液。不同的尖孢镰刀菌专化型其脂肪酸种类和组成存在着差异，为电子舌技术作为尖孢镰刀菌专化型识别的工具提供了条件。作者提出利用电子舌技术，“舔”出由气相色谱分析的尖孢镰刀菌脂肪酸组成差异，建立判别模型，鉴别尖孢镰刀菌的专化型。电子舌技术是一种由气相色谱仪和判别模型组成的脂肪酸物质差异分析机构，能识别单一的或混合的脂肪酸溶液。不同的尖孢镰刀菌专化型其脂肪酸种类和组成存在着差异，为电子舌技术作为尖孢镰刀菌专化型识别的工具提供了条件。尖孢镰刀菌微生物细胞结构中普遍
9、含有的脂肪酸成分，并与微生物 DNA 具有高度同源性，各种尖孢镰刀菌微生物专化型具有其特征性的细胞脂肪酸指纹图谱【杜宗敏，2003】，因此，脂肪酸脂类早在 60年代脂肪酸就开始作为标志物用于细菌尖孢镰刀菌鉴定了【刘志辉，2005】。同时微生物脂肪酸中含多种有对植物病原菌拮抗作用的抗菌脂肪酸( Choudhury et al., 1994; Benyagoub et al., 1996; Cantrell et al., 2008)。抗菌活性脂肪酸尤其对真菌性病原菌如链格孢属（ Alternaria solani）刺盘孢属（Colletotrichum lagenarium） 镰刀菌属（Fusa
10、rium oxysporum f. sp）,镰刀菌属（ Fusarium oxysporum f. sp lycopersici）有抑制作用，对病原菌菌丝体生长和孢子萌发都有很强的抑制作用( Liu et al., 2008)。尖孢镰刀菌（Fusarium oxysporum）的脂肪酸主要是包括 C8 到 C24(Joshi and Mathur, 1987)。，尖孢镰刀菌（Fusarium oxysporum）的脂肪酸组成 其中最主要的成分为 18:1 9 c, 16: 0 和 18:0，.占总脂肪酸的 47%，其他的脂肪酸含量低，主要为 -羟基，支链，二不饱和和奇数脂肪酸( Ellis e
11、t al., 2002)。镰刀菌（Fusarium oxysporum 和 Fusarium roseum）的菌丝体和孢子的脂肪酸包括了主要脂肪酸 18:1, 18:2, 和 16:0 和微量脂肪酸 18:0, 18:3, 和 14:0(RAMBO AND BEAN, 1969)。尖孢镰刀菌专化型脂肪酸的组成和含量不仅与菌株的生长温度有关( BARRAN AND MILLER, 1976)。尖孢镰刀菌专化型，其脂肪酸的组成结构还取决于菌株培养条件的 C 和 N 组成( JOSHI AND MATHUR, 1987)。 JOSHI AND DHAR, 1987 尖孢镰刀菌专化型脂肪酸合成一种从
12、FUSARIUM OXYSPORUM F. SP. LINI 分离到的脂肪酶对棉籽油中脂肪酸，坚果油中脂肪酸和尖孢镰孢菌分泌的饱和脂肪酸具有水解作用，显示对尖孢镰孢菌专化型分泌的饱和脂肪酸具有优先特异性( JOSHI AND DHAR, 1987)，形成了对寄主特异性选择( JOSHI AND DHAR, 1987)。这一原理，为利用电子舌技术分析尖孢镰刀菌专化型提供了基础。目前，目前脂肪酸生物标记分析较少应用于农业病原真菌的研究，对于尖孢镰刀菌专化型的镰刀菌脂肪酸组成差异研究多集中外界环境如温度，培养基，PHpH 值等的影响，而对于由于镰刀菌自身差异，如专化型寄主来源不同引起差异的研究较少。
13、到目前为止，未见关于利用电子舌技术识别尖孢镰刀菌专化型脂肪酸差异的研究的报道。目前利用判别分析法，根据不同寄主镰刀菌的脂肪酸结构特点建立镰刀菌寄主的专化型脂肪酸判别模型的报道还未见。本研究旨在对本研究室从不同寄主分离获得的尖孢镰刀菌，进行脂肪酸气相色谱典型的模式图谱及菌株脂肪酸峰值变化规律分析，为建立电子舌技术尖孢镰刀菌专化型脂肪酸判别模型提供依据，为仪器化鉴别尖孢镰刀菌专化型菌株的脂肪酸种类组群是否与专化型、生理小种、营养亲和群等生理指标具有相关性奠定基础，并建立尖孢镰刀菌寄主专化型脂肪酸判别模型，为尖孢镰刀菌寄主专化型的判定提供新的思路和方法。1 材料与方法1.1 材料3供试菌株，不同来源
14、分离纯化的尖孢镰刀菌 22 株，见表 1。 SDB 培养基：30 g Sabouraud Dextrose Broth（ SDB）+蒸馏水 1000 mL。（SDB 购于 Fisher 公司）。皂化试剂：氢氧化钠 45 g+甲醇150 mL+水 150 mL。甲基化试剂：6N 盐酸 325 mL+甲醇 275 mL。萃取试剂：正已烷 200 mL+甲基叔丁基乙醚 200 mL。洗涤试剂：氢氧化钠 10.8 g+水 900 mL。表 1 供试尖孢镰刀菌菌株寄主 菌株编号 采集地点 专化型 致病型1. 西瓜 134 永泰 西瓜专化型 致病2. 西瓜 136 永泰 西瓜专化型 致病3. 西瓜 137
15、 永泰 西瓜专化型 致病4. 西瓜 138 永泰 西瓜专化型 致病5. 西瓜 140 永泰 西瓜专化型 致病6. 西瓜 141 永泰 西瓜专化型 致病7. 香蕉 378 漳州 古巴专化型 致病8. 香蕉 384 漳州 古巴专化型 致病9. 香蕉 386 漳州 古巴专化型 致病10. 香蕉 387 漳州 古巴专化型 致病11. 香蕉 389 漳州 古巴专化型 致病12. 香蕉 391 漳州 古巴专化型 致病13. 香蕉 392 漳州 古巴专化型 致病14. 花生 816 泉州 致病15. 花生 833 惠安 致病16. 花生 834 惠安 致病17. 花生 836 惠安 致病18. 花生 841
16、 惠安 ？ 致病19. 花生 846 惠安 致病20. 花生 847 惠安 致病21. 花生 855 惠安 致病22. 花生 866 惠安 致病1.2方法（1）尖孢镰刀菌脂肪酸的提取和气相色谱分析脂肪酸提取所用的试剂：包括脂肪酸甲酯混合物标样（由美国 MIDI 公司提供）、菌落收集培养基（TSBA 培养基）、脂肪酸提取所需试剂皂化试剂、甲基化试剂、萃取试剂和洗涤试剂 （美国Fisher 公司提供），除脂肪酸甲酯混合物标样为直接应用外，其它试剂配制由 MIDI 公司提供；脂肪酸的提取包括菌落收集、细菌脂肪酸皂化释放、脂肪酸的甲基化、脂肪酸甲酯萃取和基本洗涤等步骤，具体操作由 MIDI 公司提供。
17、在色谱条件下平行分析脂肪酸甲酯混合物标样和待检样本，色谱条件如下：二阶程序升高柱温，170起始，5/min 升至 260 ，而后 40/min 升温至 310 ，维持 90S；汽化室温度 250、检测器温度 300；载气为氢气(2mL/min)、尾吹气为氮气(30mL/min)；柱前压 l0.00 psi(1psi= 6.895kPa)；进样量 ll，进样分流比 100:l19。气相色谱系统采用的是美国 Agilent 6890 型（Agilent Chromatograph: GC），包括全自动进样装置、石英毛细管柱及氢火焰离子化检测器；分析软件应用美国 MIDI 公司开发的基于细菌细胞脂肪
18、酸成分鉴定细菌的软件 MIS4.5（Microbial Identification System）和 LGS4.5（Library Generation Software）。（2）尖孢镰刀菌脂肪酸特征图谱分析4选择尖孢镰刀菌脂肪酸指纹图谱匹配率最高菌株的图谱作为模式图谱，说明尖孢镰刀菌脂肪酸种类、特征、种类鉴定参数等。对被测的22株菌株的脂肪酸特征峰的百分比含量，用DPS软件进行方差分析 21。（3）尖孢镰刀菌不同寄主脂肪酸生物标记分布通过MIDI系统检测得到的色谱图和信息表，汇总测试菌株的所有色谱峰(Peak Name)，以Pecent 作为含量，分析尖孢镰刀菌不同寄主的脂肪酸分布的种类及
19、含量。（4）尖孢镰刀菌脂肪酸生物标记聚类分析以脂肪酸为样本，以菌株为指标，构建聚类分析数据矩阵，数据中心化，以兰氏距离为相关尺度，以类平均法为聚类方法，通过DPS软件对22株尖孢镰刀菌脂肪酸色谱结果进行聚类分析，分析尖孢镰刀菌脂肪酸的同质性。（5）尖孢镰刀菌脂肪酸生物标记寄主关系以脂肪酸种类为指标，以菌株的寄主来源为样本，构建分析矩阵，以数据进行规格化转化，马氏距离为相关尺度，以可变平均法进行聚类。分析尖孢镰刀菌脂肪酸与寄主的关系。（6）尖孢镰刀菌脂肪酸生物标记的主成分分析以脂肪酸为样本，以菌株为指标，构建矩阵，进行主成分分析。通过主成分分析选择尖孢镰刀菌脂肪酸的特征脂肪酸，以累计特征值大于9
20、0%为指标，在因子得分矩阵中，从各脂肪酸组中选择一个因子得分总和最高的脂肪酸，组成尖孢镰刀菌脂肪酸中的特征脂肪酸。（7）尖孢镰刀菌种下分化脂肪酸型的判别模型的建立以特征脂肪酸为样本，以，以主成分分析结果得到的选出的特征脂肪酸为指标（Xj），被测的尖孢镰刀菌菌株（Xi）为样本，尖孢镰刀菌寄主级别（Y）为测定值，利用判别分析建立尖孢镰刀菌寄主专化型脂肪酸判别模型，组建判别分析数据矩阵如下：X11 X12 X1j Y1 X21 X22 X2j Y2 Xi1 Xi2 XIJ Yi 其中j为脂肪酸峰序号，i 为寄主类别序号。用Bayes逐步判别方法建模，多组逐步判别分析是按照Bayes 准则,用 WIL
21、KS 统计量从 M个因子中，根据各因子判别能力大小的显著性检验结果，精选出若干判别能力大的且配合较好的因子，建立Bayes判别函数。Bayes判别函数方程：（i =1, 2, , p）pjjiii xcxy10)(其中： pkiji10klkjiij wxmnc)1()(当： ，则将个体x划入第i 类。判别值的后验概率为：)()(ma1yxyiipi5phxyiiie1)()(聚类分析、主成分分析、判别分析采用DPS数据处理系统进行计算 21。2 结果与分析2.1 尖孢镰刀菌专化型菌株的 脂肪酸模式图谱分析选择尖孢镰刀菌 841 菌株作为模式菌株。气相色谱典型的模式图见图 1。从图中可以看出，
22、除第一个溶剂波峰外，其余每一个波峰代表一种脂肪酸，该菌株具有 4 个特征性色谱峰，其它波峰为非优势脂肪酸。4 个特征性色谱峰的保留时间分别为：（1）C16:0，10.9000.100min ；（2）C18:2 CIS9,12/18:0a，13.9000.100min ；（3）C18:1 CIS 9 (9c)，14.4000.100min；（4）C18:0， 14.4000.100min。在第一个溶剂波峰与 C16:0 之间还会出现 1-2 个小峰，在 C18:0 波峰后部份菌株会出现 1-2 个小峰。色谱的有效峰从保留时间1.729起点到保留时间14.46 为终点，共 7个峰，每个峰代表一类脂
23、肪酸。不同的菌株，峰的数量、峰值大小和脂肪酸种类不同。大部分菌株鉴定用的脂肪酸峰从十四碳脂肪酸（C14:0 ）开始，到十八碳脂肪酸（C18:0）结束。min2.5 5 7.5 10 12.5 15 17.5pA17.5202.52527.53032.53537.5 FID1 A, (E07130.413B0174679.D) 1.7297.648 10.62110.92713.97814.06214.64图 1 尖孢镰刀菌 841 脂肪酸气相色谱典型的模式图横坐标表示色谱分析时间(min)，纵坐标表示色谱峰高值；第一个色谱峰为溶剂波峰表 2 则为该尖孢镰刀菌 PFLAs 标记的信息表，是将图
24、1 中的气相色谱图信息进行量化，转化成数据形式。其中，retention time (RT)表示各 PLFAs 特征峰的保留时间；Response 表示各 PLFAs 特征峰的峰面积，体现了各 PLFA 的含量； area/ height ratio（Ar/Ht）表示峰面积/ 峰高度比值；RFact 表示相关因子系数；the equivalent chain length (ECL)表示各磷脂脂肪酸的相对链长度；Peak Name 指各磷脂脂肪酸特征峰名称；Percent 是指磷脂脂肪酸各特征峰占总脂肪酸的百分数，Comment1、Comment2 分别表示误差注释和鉴定参考峰注释。MIDI
25、系统一般通过各磷脂脂肪酸特征峰的种类和百分比含量进行菌株鉴定，在实验中，我们一般应用 Peak Name 和 Percent 来分析不同寄主专化型尖孢镰刀菌的脂肪酸组成特点。表2 尖孢镰刀菌菌株 841脂肪酸气相色谱典型的模式图对应的分析结果RT Response Ar/Ht RFact ECL Peak Name Percent Comment1 Comment21.729 3.652E+8 0.027 - 7.011 SOLVENT PEAK - min rt7.648 1638 0.038 1.052 14.000 14:0 0.66 ECL deviates 0.000 Referen
26、ce -0.002610.621 3605 0.041 0.887 15.820 16:1 Cis 9 (w 7) 1.23 ECL deviates 0.00310.927 65136 0.043 0.874 16.000 16:0 21.93 ECL deviates 0.000 Reference -0.00313.978 126956 0.048 0.777 17.727 18:2 CIS 9,12/18:0a 37.97 ECL deviates 0.00714.062 107332 0.049 0.774 17.774 Sum In Feature 8 32.00 ECL devi
27、ates 0.001 18:1 CIS 9 (w 9) 14.464 21115 0.046 0.763 18.000 18:0 6.21 ECL deviates 0.000 Reference -0.003- 107332 - - - Summed Feature 8 32.00 18:1 CIS 9 (w 9) 18:1 (w 8) 2.2 尖孢镰刀菌专化型不同寄主 脂肪酸生物标记测定分布根据上述脂肪酸测定数据的分析结果，将供试的 22 株尖孢镰刀菌脂肪酸分布情况见表 3。共检测到 8 种脂肪酸。其中脂肪酸 16:0、18:0 和 18:2 CIS 9,12/18:0a 在所有的菌株中都
28、存在，脂肪酸 16:0的含量在 16.7%-24.17%之间，脂肪酸 18:0 的含量在 3.98%-10.92%之间，脂肪酸 18:2 CIS 9,12/18:0a 的含量在 32.55%-58.71%之间，各菌株三种脂肪酸之和占总脂肪酸含量的 60%以上，为主要脂肪酸。脂肪酸 16:1 CIS 9 (W 7)在 18 个菌株中存在，含量在 0-1.59%；脂肪酸 14:0 在 12 个菌株中存在，含量在0.53%-1.72%；脂肪酸 13:0 2OH、15:1 ISO G 和 17:0 分别仅在 1 个菌株中检测到，含量均低于 1%。表3 尖孢镰刀菌脂肪酸生物标记脂肪酸含量（%）寄主 菌株
29、编号 采集地点 13:0 2OH 14:015:1 ISO G 16:016:1 CIS 9 (W 7) 17:0 18:018:2 CIS 9,12/18:0a 西瓜 134 永泰 0 0.67 0 22.19 0.98 0 6.28 39.33西瓜 136 永泰 0 0.58 0 22.88 1.00 0 6.66 41.57西瓜 137 永泰 0 0.53 0 23.95 0.82 0 8.39 39.72西瓜 138 永泰 0 0 0 16.99 0.88 0 5.61 48.95西瓜 140 永泰 0.50 1.65 0.90 21.17 1.16 0 3.98 39.63西瓜 14
30、1 永泰 0 0 0 21.20 0.00 0 5.40 49.19香蕉 378 漳州 0 0.64 0 23.52 1.32 0 8.15 32.86香蕉 384 漳州 0 0 0 23.57 0.88 0 8.60 32.55香蕉 386 漳州 0 0 0 20.32 1.09 0.56 6.63 45.75香蕉 387 漳州 0 0 0 23.34 1.17 0 6.78 43.59香蕉 389 漳州 0 0 0 24.17 0.00 0 8.45 36.70香蕉 391 漳州 0 0 0 20.70 1.02 0 5.68 39.10香蕉 392 漳州 0 0 0 18.95 0.75
31、 0 8.62 36.68花生 816 泉州 0 0.55 0 20.81 0.86 0 7.39 35.71花生 833 惠安 0 0.79 0 20.53 1.31 0 10.92 35.04花生 834 惠安 0 0.00 0 22.00 1.59 0 9.80 36.94花生 836 惠安 0 1.56 0 21.93 1.27 0 6.68 41.09花生 841 惠安 0 0.66 0 21.93 1.23 0 6.21 37.97花生 846 惠安 0 0.00 0 21.89 1.53 0 9.65 37.81花生 847 惠安 0 1.11 0 16.70 0.00 0 4.
32、83 58.71花生 855 惠安 0 1.72 0 19.18 0.00 0 5.09 38.59花生 866 惠安 0 0.55 0 17.58 1.04 0 6.30 36.582.3 基于尖孢镰刀菌专化型的 脂肪酸生物标记聚类分析以脂肪酸为样本，以在菌株中存在情况指标，构建分析矩阵，以兰氏距离为相关尺度，以最长距离法进行聚类，分析结果见表 4 及图 2 。当 =20.78 时，可将 8 个脂肪酸分布情况分为 3 个类群。类群为偶现脂肪酸类型，包括脂肪酸 13:0 2OH、15:1 ISO G 和 17:0，分别仅在 1 个菌株中出现。类群为不完全分布脂肪酸类型，即该7类中的脂肪酸仅在部
33、分被测菌株中出现，它们包括了脂肪酸 14:0 和 16:1 CIS 9 (W 7)，分别在 12 个菌株和 18 个菌株中分布。类群为完全分布脂肪酸类型，即该类中的脂肪酸在全部被测菌株中出现，它们包括了脂肪酸 16:0、18:0 和 18:2 CIS 9,12/18:0a，在 22 个菌株中均有分布，且在各菌株中含量之和大于 60%，为主要脂肪酸。通过对尖孢镰刀菌脂肪酸生物标记聚类分析可知，3 个偶现脂肪酸生物标记，对与寄主的关系不密切，在以下分析尖孢镰刀菌脂肪酸生物标记寄主关系过程，将剔出这 3 条脂肪酸生物标记。表 4 聚类过程T I J 距离1 3 1 0.285712 6 1 23 8
34、 4 7.865664 5 2 14.723175 7 4 18.049466 2 1 21.397597 4 1 25.97508 图 2 脂肪酸分布情况聚类结果2.4 尖孢镰刀菌脂肪酸生物标记 寄主关系以脂肪酸为指标，以寄主来源为样本，构建分析矩阵，数据进行规格化转化，马氏距离为相关尺度，以可变平均法进行聚类，分析结果见表 5 及图 3。当 =17.22 时，可将 22 株供试菌株分为 2个类群。类群包括 10 株尖孢镰刀菌，6 株来源于西瓜，4 株来源于香蕉。类群包括 5 株尖孢镰刀菌，2 株来源于香蕉，3 株来源于花生。类群包括 7 株尖孢镰刀菌，其中 6 株来源于花生，1 株来源于香
35、蕉。来源于西瓜和花生的镰刀菌大部分在同一类中，并分布在聚类的上下区域，而来源于香蕉的菌株分布在聚类图的中间区域，并有部分穿插在西瓜和花生中，可见尖孢镰刀菌的寄主也具有过渡寄主，香蕉为尖孢镰刀菌寄主西瓜和花生的过渡性寄主。表 5T I J 距离 1 7 3 0.61 2 2 1 0.83 3 9 6 0.89 4 22 21 0.96 5 17 11 1.29 6 5 4 1.82 7 10 6 2.18 8 20 16 2.25 9 21 13 2.78 脂肪酸 1=28.70G1810 4 1 3.06 11 14 12 3.36 12 8 6 4.19 13 3 1 4.88 14 19
36、16 5.76 15 6 1 6.13 16 18 11 7.27 17 12 11 9.11 18 15 13 11.78 19 11 1 13.41 20 16 13 15.79 21 13 1 28.70 图 3 不同寄主尖孢镰刀菌脂肪酸聚类图当 =20.78 时，可将 8 个脂肪酸分布情况分为 3 个类群。类群为偶现脂肪酸类型，包括脂肪酸 13:0 2OH、15:1 ISO G 和 17:0，分别仅在 1 个菌株中出现。类群为不完全分布脂肪酸类型，即该类中的脂肪酸仅在部分被测菌株中出现，它们包括了脂肪酸 14:0 和 16:1 CIS 9 (W 7)，分别在 12 个菌株和 18 个菌
37、株中分布。类群为完全分布脂肪酸类型，即该类中的脂肪酸在全部被测菌株中出现，它们包括了脂肪酸 16:0、18:0 和 18:2 CIS 9,12/18:0a，在 22 个菌株中均有分布，且在各菌株中含量之和大于 60%，为主要脂肪酸。通过对尖孢镰刀菌脂肪酸生物标记聚类分析可知，3 个偶现脂肪酸生物标记，对与寄主的关系不密切，在以下分析尖孢镰刀菌脂肪酸生物标记寄主关系过程，将剔出这3 条脂肪酸生物标记。2.45 尖孢镰刀菌专化型脂肪酸生物标记的主成分分析（要用全部脂肪酸数据分析，尖孢镰刀菌脂肪酸生物标记的主成分分析分析结果见表 6。前四四个主成分的特征值的累计百分比达 93.91%，表明前四个主成
38、分代表了总体信息量的 90%以上。规格化特征向量见表 7，有 5 个因子，分别为因子 1（14:0），因子 2（16:0 ），因子 3（16:1 CIS 9 (W 7)），因子 4（18:0），因子5（18:2 CIS 9,12/18:0a）；它们的特征向量均值分别为 0.12，0.13、0，0.03，0.16 ，0.30。主成分分析结果得到尖孢镰刀菌脂肪酸中的特征脂肪酸有 55 个，分别为 14:0，16:0 ，16:1 CIS 9 (W 7)，18:0 和18:2 CIS 9,12/18:0a。表6 尖孢镰刀菌脂肪酸生物标记主成分分析特征值No 特征值 百分率% 累计百分率%1 2.41
39、48.23 48.232 1.03 20.52 68.753 0.75 14.97 83.724 0.51 10.19 93.915 0.30 6.09 100.00表7 规格化特征向量:因子 x(1) x(2) x(3) x(4) x(5) 总和均值因子 1 -0.27 0.44 0.42 0.54 -0.52 0.6012 因子 2 0.86 0.13 0.25 -0.22 -0.36 0.6513 因子 3 -0.04 -0.69 0.72 0.08 0.10 0.1503 因子 4 -0.17 0.53 0.50 -0.50 0.43 0.8016 因子 5 0.39 0.18 0.0
40、3 0.64 0.64 10.9038 92.5 基于尖孢镰刀菌脂肪酸生物 标记的尖孢镰刀菌寄主专化型聚类分析关系根据主成分分析的结果，选择主成分分析结果得到尖孢镰刀菌脂肪酸中的 5 个特征脂肪酸，即14:0、16:0、16:1 CIS 9 (W 7)、18:0 和 18:2 CIS 9,12/18:0a，构建矩阵，以脂肪酸为指标，以寄主来源为样本，数据进行规格化转化，马氏距离为相关尺度，以可变平均法进行聚类，分析结果见表 5 及图3。当 1=28.70 时，可将 22 株供试菌株分为 2 个类群。类群包括 15 株尖孢镰刀菌，6 株来源于西瓜，5 株来源于香蕉，3 株来源于花生。类群包括 7
41、 株尖孢镰刀菌，1 株来源于香蕉，6 株来源于花生。在类群 I 中，还可以分为三个亚群 S1、S2 、S3；在 S1 亚群中，5 株来源于西瓜，1 株来源于香蕉；在 S2 亚群中，1 株来源于西瓜，3 株来源于香蕉；在 S3 亚群中，2 株来源于香蕉，3 株来源于花生。从聚类分析的结果可以看出，寄主为西瓜和花生的尖孢镰刀菌是聚类的两个极端，寄主为香蕉的尖孢镰刀菌，分布在聚类图的中间区域，并有部分穿插在西瓜和花生中，可见在尖孢镰刀菌脂肪酸生物标记为指标的进化过程，来自香蕉的尖孢镰刀菌起着在西瓜和花生的尖孢镰刀菌之间的过渡性进化作用。当 1=28.70 时，可将 22 株供试菌株分为 2 个类群。
42、类群包括 15 株尖孢镰刀菌，6 株来源于西瓜，5 株来源于香蕉，3 株来源于花生。类群包括 7 株尖孢镰刀菌，1 株来源于香蕉，6 株来源于花生。在类群 I 中，还可以分为三个亚群 S1、S2 、S3；在 S1 亚群中，5 株来源于西瓜，1 株来源于香蕉；在 S2 亚群中，1 株来源于西瓜，3 株来源于香蕉；在 S3 亚群中，2 株来源于香蕉，3 株来源于花生。从聚类分析的结果可以看出，来致于寄主为西瓜和花生的尖孢镰刀菌是聚类的两个极端，来致于寄主为香蕉的尖孢镰刀菌，分布在聚类图的中间区域，并有部分穿插在西瓜和花生中，可见在尖孢镰刀菌脂肪酸生物标记为指标的进化过程，来致自香蕉的尖孢镰刀菌起着
43、在西瓜和花生的尖孢镰刀菌之间的过渡性进化作用。表 5T I J 距离 1 7 3 0.61 2 2 1 0.83 3 9 6 0.89 4 22 21 0.96 5 17 11 1.29 6 5 4 1.82 1=28.70G1G2S1S3S2107 10 6 2.18 8 20 16 2.25 9 21 13 2.78 10 4 1 3.06 11 14 12 3.36 12 8 6 4.19 13 3 1 4.88 14 19 16 5.76 15 6 1 6.13 16 18 11 7.27 17 12 11 9.11 18 15 13 11.78 19 11 1 13.41 20 16
44、 13 15.79 21 13 1 28.70 图 3 不同寄主尖孢镰刀菌脂肪酸聚类图2.6 不同寄主尖孢镰刀菌专化型 电子舌鉴别脂肪酸生物标记判别模型的建立2.6.1 不同寄主尖孢镰刀菌脂肪酸生物 标记判别模型的原始数据及已知分类以主成分分析结果选出的特征脂肪酸包括 14:0，16:0 ，16:1 CIS 9 (W 7)，18:0，18:2 CIS 9,12/18:0a为指标（Xj ），尖孢镰刀菌的寄主西瓜，香蕉和花生为样本，寄主为西瓜，香蕉和花生，寄主类别（Y ）为测定值，其中定义西瓜为寄主类别 1，香蕉为寄主类别 2，花生为寄主类别 3，组建判别分析数据矩阵见表 8。表8 尖孢镰刀菌寄主
45、专化型判别分析矩阵脂肪酸含量（% ）寄主 采集地点 菌株编号 14:0 16:0 16:1 CIS 9 (W 7) 18:0 18:2 CIS 9,12/18:0a寄主类别西瓜 永泰 134 0.67 22.19 0.98 6.28 39.33 1西瓜 永泰 136 0.58 22.88 1.00 6.66 41.57 1西瓜 永泰 137 0.53 23.95 0.82 8.39 39.72 1西瓜 永泰 138 0.00 16.99 0.88 5.61 48.95 1西瓜 永泰 140 1.65 21.17 1.16 3.98 39.63 1西瓜 永泰 141 0.00 21.20 0.0
46、0 5.40 49.19 1香蕉 漳州 378 0.64 23.52 1.32 8.15 32.86 2香蕉 漳州 384 0.00 23.57 0.88 8.60 32.55 2香蕉 漳州 386 0.00 20.32 1.09 6.63 45.75 2香蕉 漳州 387 0.00 23.34 1.17 6.78 43.59 2香蕉 漳州 389 0.00 24.17 0.00 8.45 36.70 2香蕉 漳州 391 0.00 20.70 1.02 5.68 39.10 2香蕉 漳州 392 0.00 18.95 0.75 8.62 36.68 2花生 泉州 816 0.55 20.81
47、 0.86 7.39 35.71 3花生 惠安 833 0.79 20.53 1.31 10.92 35.04 3花生 惠安 834 0.00 22.00 1.59 9.80 36.94 3花生 惠安 836 1.56 21.93 1.27 6.68 41.09 3花生 惠安 841 0.66 21.93 1.23 6.21 37.97 3花生 惠安 846 0.00 21.89 1.53 9.65 37.81 3花生 惠安 847 1.11 16.70 0.00 4.83 58.71 311寄主 采集地点 菌株编号 脂肪酸含量（ %） 寄主类别14:0 16:0 16:1 CIS 9 (W
48、7) 18:0 18:2 CIS 9,12/18:0a花生 惠安 855 1.72 19.18 0.00 5.09 38.59 3花生 惠安 866 0.55 17.58 1.04 6.30 36.58 3按照 Bayes 准则，采用 Bayes 逐步判别分析，当 Fa 临界值 =0.40 时，尖孢镰刀菌寄主专化型判别分析特征脂肪酸 14:0，16:0，16:1 CIS 9 (W 7)，18:0，18:2 CIS 9,12/18:0a 在三个寄主专化型类别中的平均值各个因子分类均值见表 9.。第一类即菌株寄主为西瓜时， 五个因子（X1 ，X2，X3，X4，X5）即脂肪酸 14:0， 16:0，
49、16:1 CIS 9 (W 7)，18:0，18:2 CIS 9,12/18:0a 的均值分别为0.57，21.40 ，0.81，6.05，43.07；第二类，即菌株寄主为香蕉时，脂肪酸 14:0，16:0，16:1 CIS 9 (W 7)，18:0，18:2 CIS 9,12/18:0a 的均值为 0.09 ，22.08，0.89，7.56 ，38.18；第三类，即菌株寄主为花生时，脂肪酸 14:0，16:0，16:1 CIS 9 (W 7)，18:0，18:2 CIS 9,12/18:0a 的均值 0.77，20.28，0.98 ，7.43，39.83。各个因子类内离差（W）矩阵和总离差（T）矩阵见表 10 和表 11，。表 9 尖孢镰刀菌寄主专化型判别分析各个因子分类均值:寄主 14:0 16:0 16:1 CIS 9 (W 7) 18:0 18:2 CIS 9,12/18:0a西瓜 0.57 21.40 0.81 6.05 43.07 香蕉 0.09 22.08 0.89 7.
