尖孢镰刀菌寄主专化型电子舌鉴别模型的建立.doc

19

尖孢镰刀菌寄主专化型电子舌性脂肪酸判别鉴别模型的建立 收稿日期： 接受日期：
基金项目：国家高技术研究与发展计划（863）项目（2006AA10A211）；福建省财政专项-福建省农业科学院科技创新团队建设基金(STIF-Y03)；福建省发改委重点项目（闽发改投资[2006]781号）。福建省自然科学基金：非致病性尖孢镰刀菌对西瓜枯萎病的免疫抗病机制研究（2008J0055）
作者简介：蓝江林（1972-），博士，副研究员，研究方向：生物技术与生物防治。电子邮件：lanfz2008@。电话：[phone]。
通讯作者：刘波，博士，研究员。研究方向：生物技术与生物防治。电子邮件：fzliubo@。电话：[phone]。

刘波*，肖荣凤，朱育菁，胡桂萍
（福建省农业科学院农业生物资源研究所，福建 福州 350003）

摘要：应用基于通过微生物细胞脂肪酸成分鉴定的全自动微生物鉴定系统(Sherlock MIS)，检测分析不同寄主尖孢镰刀菌中脂肪酸的生物标记结构组成和的分布特性。，从供试的22个菌株中检测到8种脂肪酸，其中常见脂肪酸有6条，含量60%以上的主要脂肪酸脂肪酸有16:0、18:0 和18:2 CIS 9,12/18:0a 在所有的菌株中都存在，各菌株中3种脂肪酸之和占总脂肪酸含量的60%以上，为主要脂肪酸。对常见脂肪酸进行通过系统聚类分析，可将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+10.40X1+7.96X2+5.40X4+3.42X5，Y3=-173.97+14.85X1+7.25X2+6.58X4+3.53X5，当1≦Yi﹤2时，菌株的寄主为西瓜； 当2≦Yi﹤3时,菌株的寄主为香蕉；当 Yi﹥3时，菌株的寄主为花生，且模型判对的概率为0.90909，说明脂肪酸14:0， 16:0， 18:0， 18:2 CIS 9,12/18:0a是判别尖孢镰刀菌寄主的四个重要指标，同时该模型比较符合实际，有很好的应用价值。

关键词：尖孢镰刀菌；脂肪酸；寄主；判别模型

Abstract:Characteristics of PLFA biomarkers in Fusarium 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 into 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多种植物发病，造成严重减产，甚至绝收。尖孢镰刀菌由致病菌株与无致病菌株组成，种内的致病菌株具有高度的寄主专化性，根据对不同寄主植物的专化性可分为专化型，在同一尖孢镰刀菌不同的的专化型，引起的发病程度完全不同，如尖孢镰刀菌西瓜专化型就不会侵染黄瓜专化型，同样，尖孢镰刀菌古巴专化型只会侵染香蕉而不会侵染花生作物下面。，根据致病力的差异又划分为许多生理小种。尖孢镰刀菌专化型致病型的监测是研究病原菌群体结构的一项重要措施，为病害的防治提供可靠的病原菌信息，并可为抗病育种提供抗性基因。植物病理学家和育种学家认识到病原真菌的变异对寄主植物抗性的有效性和持久性非常重要，病原真菌的基因型与特定的寄主植物基因型相互作用可导致抗病性在很短时间内丧失【BROWN J K，1996】。
然而，鉴定尖孢镰刀菌的专化型，方法复杂、时间很长、实验难度较大。电子舌技术是一种由气相色谱仪和判别模型组成的脂肪酸物质差异分析机构，能识别单一的或混合的脂肪酸溶液。不同的尖孢镰刀菌专化型其脂肪酸种类和组成存在着差异，为电子舌技术作为尖孢镰刀菌专化型识别的工具提供了条件。作者提出利用电子舌技术，“舔”出由气相色谱分析的尖孢镰刀菌脂肪酸组成差异，建立判别模型，鉴别尖孢镰刀菌的专化型。电子舌技术是一种由气相色谱仪和判别模型组成的脂肪酸物质差异分析机构，能识别单一的或混合的脂肪酸溶液。不同的尖孢镰刀菌专化型其脂肪酸种类和组成存在着差异，为电子舌技术作为尖孢镰刀菌专化型识别的工具提供了条件。
尖孢镰刀菌
微生物细胞结构中普遍含有的脂肪酸成分，并与微生物DNA具有高度同源性，各种尖孢镰刀菌微生物专化型具有其特征性的细胞脂肪酸指纹图谱【杜宗敏，2003】，因此，脂肪酸脂类早在60年代脂肪酸就开始作为标志物用于细菌尖孢镰刀菌鉴定了【刘志辉，2005】。同时微生物脂肪酸中含多种有对植物病原菌拮抗作用的抗菌脂肪酸(Choudhury et al., 1994; Benyagoub et al., 1996; Cantrell et al., 2008)。抗菌活性脂肪酸尤其对真菌性病原菌如链格孢属（Alternaria solani）刺盘孢属（Colletotrichum lagenarium） 镰刀菌属（Fusarium 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 et 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尖孢镰刀菌专化型脂肪酸合成一种从Fusarium oxysporum f. sp. lini分离到的脂肪酶对棉籽油中脂肪酸，坚果油中脂肪酸和尖孢镰孢菌分泌的饱和脂肪酸具有水解作用，显示对尖孢镰孢菌专化型分泌的饱和脂肪酸具有优先特异性(Joshi and Dhar, 1987)，形成了对寄主特异性选择(Joshi and Dhar, 1987)。这一原理，为利用电子舌技术分析尖孢镰刀菌专化型提供了基础。目前，
目前脂肪酸生物标记分析较少应用于农业病原真菌的研究，对于尖孢镰刀菌专化型的镰刀菌脂肪酸组成差异研究多集中外界环境如温度，培养基，PHpH值等的影响，而对于由于镰刀菌自身差异，如专化型寄主来源不同引起差异的研究较少。到目前为止，未见关于利用电子舌技术识别尖孢镰刀菌专化型脂肪酸差异的研究的报道。目前利用判别分析法，根据不同寄主镰刀菌的脂肪酸结构特点建立镰刀菌寄主的专化型脂肪酸判别模型的报道还未见。
本研究旨在对本研究室从不同寄主分离获得的尖孢镰刀菌，进行脂肪酸气相色谱典型的模式图谱及菌株脂肪酸峰值变化规律分析，为建立电子舌技术尖孢镰刀菌专化型脂肪酸判别模型提供依据，为仪器化鉴别尖孢镰刀菌专化型菌株的脂肪酸种类组群是否与专化型、生理小种、营养亲和群等生理指标具有相关性奠定基础，并建立尖孢镰刀菌寄主专化型脂肪酸判别模型，为尖孢镰刀菌寄主专化型的判定提供新的思路和方法。
1材料与方法
1.1材料
供试菌株，不同来源分离纯化的尖孢镰刀菌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
永泰
西瓜专化型
致病
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
惠安
？？？？
致病
19. 花生
846
惠安

致病
20. 花生
847
惠安

致病
21. 花生
855
惠安

致病
22. 花生
866
惠安

致病

1.2方法
（1）尖孢镰刀菌脂肪酸的提取和气相色谱分析
脂肪酸提取所用的试剂：包括脂肪酸甲酯混合物标样（由美国MIDI公司提供）、菌落收集培养基（TSBA培养基）、脂肪酸提取所需试剂[皂化试剂、甲基化试剂、萃取试剂和洗涤试剂]（美国Fisher公司提供），除脂肪酸甲酯混合物标样为直接应用外，其它试剂配制由MIDI公司提供；脂肪酸的提取包括菌落收集、细菌脂肪酸皂化释放、脂肪酸的甲基化、脂肪酸甲酯萃取和基本洗涤等步骤，具体操作由MIDI公司提供。在色谱条件下平行分析脂肪酸甲酯混合物标样和待检样本，色谱条件如下：二阶程序升高柱温，170℃起始，5℃/min升至260℃ ，而后40℃/min升温至310℃ ，维持90S；汽化室温度250℃、检测器温度300℃；载气为氢气(2mL/min)、尾吹气为氮气(30mL/min)；柱前压l0.00 psi(1psi= 6.895kPa)；进样量lμl，进样分流比100:l[19]。气相色谱系统采用的是美国Agilent 6890型（Agilent Chromatograph: GC），包括全自动进样装置、石英毛细管柱及氢火焰离子化检测器；分析软件应用美国MIDI公司开发的基于细菌细胞脂肪酸成分鉴定细菌的软件MIS4.5（Microbial Identification System）和 LGS4.5（Library Generation Software）。
（2）尖孢镰刀菌脂肪酸特征图谱分析
选择尖孢镰刀菌脂肪酸指纹图谱匹配率最高菌株的图谱作为模式图谱，说明尖孢镰刀菌脂肪酸种类、特征、种类鉴定参数等。对被测的22株菌株的脂肪酸特征峰的百分比含量，用DPS软件进行方差分析[21]。
（3）尖孢镰刀菌不同寄主脂肪酸生物标记分布
通过MIDI系统检测得到的色谱图和信息表，汇总测试菌株的所有色谱峰(Peak Name)，以Pecent作为含量，分析尖孢镰刀菌不同寄主的脂肪酸分布的种类及含量。
（4）尖孢镰刀菌脂肪酸生物标记聚类分析
以脂肪酸为样本，以菌株为指标，构建聚类分析数据矩阵，数据中心化，以兰氏距离为相关尺度，以类平均法为聚类方法，通过DPS软件对22株尖孢镰刀菌脂肪酸色谱结果进行聚类分析，分析尖孢镰刀菌脂肪酸的同质性。
（5）尖孢镰刀菌脂肪酸生物标记寄主关系
以脂肪酸种类为指标，以菌株的寄主来源为样本，构建分析矩阵，以数据进行规格化转化，马氏距离为相关尺度，以可变平均法进行聚类。分析尖孢镰刀菌脂肪酸与寄主的关系。
（6）尖孢镰刀菌脂肪酸生物标记的主成分分析
以脂肪酸为样本，以菌株为指标，构建矩阵，进行主成分分析。通过主成分分析选择尖孢镰刀菌脂肪酸的特征脂肪酸，以累计特征值大于90%为指标，在因子得分矩阵中，从各脂肪酸组中选择一个因子得分总和最高的脂肪酸，组成尖孢镰刀菌脂肪酸中的特征脂肪酸。
（7）尖孢镰刀菌种下分化脂肪酸型的判别模型的建立
以特征脂肪酸为样本，以，以主成分分析结果得到的选出的特征脂肪酸为指标（Xj），被测的尖孢镰刀菌菌株（Xi）为样本，尖孢镰刀菌寄主级别（Y）为测定值，利用判别分析建立尖孢镰刀菌寄主专化型脂肪酸判别模型，组建判别分析数据矩阵如下：
X11
X12
…
X1j
Y1
X21
X22
…
X2j
Y2
…
…
…
…
…
Xi1
Xi2
…
Xij
Yi

其中j为脂肪酸峰序号，i为寄主类别序号。用Bayes逐步判别方法建模，多组逐步判别分析是按照Bayes准则,用WILKS 统计量从M个因子中，根据各因子判别能力大小的显著性检验结果，精选出若干判别能力大的且配合较好的因子，建立Bayes判别函数。
Bayes判别函数方程：
（i=1, 2, …, p）
其中：

当：，则将个体x划入第i类。判别值的后验概率为：

聚类分析、主成分分析、判别分析采用DPS数据处理系统进行计算[21]。

2结果与分析
2.1尖孢镰刀菌专化型菌株的脂肪酸模式图谱分析
选择尖孢镰刀菌841菌株作为模式菌株。气相色谱典型的模式图见图1。从图中可以看出，除第一个溶剂波峰外，其余每一个波峰代表一种脂肪酸，该菌株具有4个特征性色谱峰，其它波峰为非优势脂肪酸。4个特征性色谱峰的保留时间分别为：（1）C16:0，10.900±0.100min；（2）C18:2 CIS9,12/18:0a，13.900±0.100min；（3）C18:1 CIS 9 (ω9c)，14.400±0.100min；（4）C18:0，14.400±0.100min。在第一个溶剂波峰与C16:0之间还会出现1-2个小峰，在C18:0波峰后部份菌株会出现1-2个小峰。
色谱的有效峰从保留时间1.729起点到保留时间14.46为终点，共7个峰，每个峰代表一类脂肪酸。不同的菌株，峰的数量、峰值大小和脂肪酸种类不同。大部分菌株鉴定用的脂肪酸峰从十四碳脂肪酸（C14:0）开始，到十八碳脂肪酸（C18:0）结束。

图1尖孢镰刀菌841脂肪酸气相色谱典型的模式图
横坐标表示色谱分析时间(min)，纵坐标表示色谱峰高值；第一个色谱峰为溶剂波峰

表2则为该尖孢镰刀菌PFLAs标记的信息表，是将图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系统一般通过各磷脂脂肪酸特征峰的种类和百分比含量进行菌株鉴定，在实验中，我们一般应用Peak Name和Percent来分析不同寄主专化型尖孢镰刀菌的脂肪酸组成特点。

表2 尖孢镰刀菌菌株841脂肪酸气相色谱典型的模式图对应的分析结果
RT
Response
Ar/Ht
RFact
ECL
Peak Name
Percent
Comment1
Comment2
1.729
3.652E+8
0.027
----
7.011
SOLVENT PEAK
----
< min rt

7.648
1638
0.038
1.052
14.000
14:0
0.66
ECL deviates 0.000
Reference -0.002
10.621
3605
0.041
0.887
15.820
16:1 Cis 9 (w 7)
1.23
ECL deviates 0.003

10.927
65136
0.043
0.874
16.000
16:0
21.93
ECL deviates 0.000
Reference -0.003
13.978
126956
0.048
0.777
17.727
18:2 CIS 9,12/18:0a
37.97
ECL deviates 0.007

14.062
107332
0.049
0.774
17.774
Sum In Feature 8
32.00
ECL deviates 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在所有的菌株中都存在，脂肪酸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 尖孢镰刀菌脂肪酸生物标记
寄主
菌株编号
采集地点
脂肪酸含量（%）
13:0
2OH
14:0
15:1
ISO G
16:0
16:1
CIS 9 (W 7)
17:0
18:0
18: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
西瓜
141
永泰
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
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.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.58

2.3基于尖孢镰刀菌专化型的脂肪酸生物标记聚类分析
以脂肪酸为样本，以在菌株中存在情况指标，构建分析矩阵，以兰氏距离为相关尺度，以最长距离法进行聚类，分析结果见表4及图2 。
当λ=20.78时，可将8个脂肪酸分布情况分为3个类群。类群Ⅰ为偶现脂肪酸类型，包括脂肪酸13:0 2OH、15:1 ISO G和17:0，分别仅在1个菌株中出现。类群Ⅱ为不完全分布脂肪酸类型，即该类中的脂肪酸仅在部分被测菌株中出现，它们包括了脂肪酸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聚类过程
脂肪酸

图2 脂肪酸分布情况聚类结果
T
I
J
距离
1
3
1
0.28571
2
6
1
2
3
8
4
7.86566
4
5
2
14.72317
5
7
4
18.04946
6
2
1
21.39759
7
4
1
25.97508

2.4尖孢镰刀菌脂肪酸生物标记寄主关系
以脂肪酸为指标，以寄主来源为样本，构建分析矩阵，数据进行规格化转化，马氏距离为相关尺度，以可变平均法进行聚类，分析结果见表5及图3。当λ=17.22时，可将22株供试菌株分为2个类群。类群Ⅰ包括10株尖孢镰刀菌，6株来源于西瓜，4株来源于香蕉。类群Ⅱ包括5株尖孢镰刀菌，2株来源于香蕉，3株来源于花生。类群Ⅲ包括7株尖孢镰刀菌，其中6株来源于花生，1株来源于香蕉。来源于西瓜和花生的镰刀菌大部分在同一类中，并分布在聚类的上下区域，而来源于香蕉的菌株分布在聚类图的中间区域，并有部分穿插在西瓜和花生中，可见尖孢镰刀菌的寄主也具有过渡寄主，香蕉为尖孢镰刀菌寄主西瓜和花生的过渡性寄主。

表5
T
I
J
距离
G1
λ1=28.70

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
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
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个菌株中分布。类群Ⅲ为完全分布脂肪酸类型，即该类中的脂肪酸在全部被测菌株中出现，它们包括了脂肪酸16:0、18:0和18:2 CIS 9,12/18:0a，在22个菌株中均有分布，且在各菌株中含量之和大于60%，为主要脂肪酸。通过对尖孢镰刀菌脂肪酸生物标记聚类分析可知，3个偶现脂肪酸生物标记，对与寄主的关系不密切，在以下分析尖孢镰刀菌脂肪酸生物标记寄主关系过程，将剔出这3条脂肪酸生物标记。

2.45尖孢镰刀菌专化型脂肪酸生物标记的主成分分析（要用全部脂肪酸数据分析，，，
尖孢镰刀菌脂肪酸生物标记的主成分分析分析结果见表6。前四四个主成分的特征值的累计百分比达93.91%，表明前四个主成分代表了总体信息量的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
48.23
48.23
2
1.03
20.52
68.75
3
0.75
14.97
83.72
4
0.51
10.19
93.91
5
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.03
0.64
0.64
10.9038

2.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株尖孢镰刀菌，1株来源于香蕉，6株来源于花生。在类群I中，还可以分为三个亚群S1、S2、S3；在S1 亚群中，5株来源于西瓜，1株来源于香蕉；在S2亚群中，1株来源于西瓜，3株来源于香蕉；在S3亚群中，2株来源于香蕉，3株来源于花生。从聚类分析的结果可以看出，寄主为西瓜和花生的尖孢镰刀菌是聚类的两个极端，寄主为香蕉的尖孢镰刀菌，分布在聚类图的中间区域，并有部分穿插在西瓜和花生中，可见在尖孢镰刀菌脂肪酸生物标记为指标的进化过程，来自香蕉的尖孢镰刀菌起着在西瓜和花生的尖孢镰刀菌之间的过渡性进化作用。

当λ1=28.70时，可将22株供试菌株分为2个类群。类群Ⅰ包括15株尖孢镰刀菌，6株来源于西瓜，5株来源于香蕉，3株来源于花生。类群Ⅱ包括7株尖孢镰刀菌，1株来源于香蕉，6株来源于花生。在类群I中，还可以分为三个亚群S1、S2、S3；在S1 亚群中，5株来源于西瓜，1株来源于香蕉；在S2亚群中，1株来源于西瓜，3株来源于香蕉；在S3亚群中，2株来源于香蕉，3株来源于花生。从聚类分析的结果可以看出，来致于寄主为西瓜和花生的尖孢镰刀菌是聚类的两个极端，来致于寄主为香蕉的尖孢镰刀菌，分布在聚类图的中间区域，并有部分穿插在西瓜和花生中，可见在尖孢镰刀菌脂肪酸生物标记为指标的进化过程，来致自香蕉的尖孢镰刀菌起着在西瓜和花生的尖孢镰刀菌之间的过渡性进化作用。

表5
T
I
J
距离
S3
S2
S1
G2
G1
λ1=28.70

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
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
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 尖孢镰刀菌寄主专化型判别分析矩阵
寄主
采集地点
菌株
编号
脂肪酸含量（%）
寄主
类别
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.00
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
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