張文倩，胡　林，陳　立，曾任森(華南農(nóng)業(yè)大學農(nóng)學院，廣東廣州5064; 中國科學院動物研究所/農(nóng)業(yè)蟲害鼠害綜合治理研究國家重點實驗室，北京000)

基于主成分分析和聚類分析的紅火蟻毒腺生物堿的比較研究

張文倩1，2，胡林2，陳立2，曾任森1
(1華南農(nóng)業(yè)大學農(nóng)學院，廣東廣州510642; 2中國科學院動物研究所/農(nóng)業(yè)蟲害鼠害綜合治理研究國家重點實驗室，北京100101)

摘要:【目的】分析中國和美國2個入侵地紅火蟻Solenopsis invicta毒腺生物堿的主要化學成分及其聚類特征.【方法】以中國廣州、惠州、深圳、珠海和東莞以及美國阿拉巴馬州采樣地共24巢紅火蟻為樣本，提取各個蟻巢毒腺的生物堿流分，分別以紅火蟻順式生物堿流分中的15種生物堿和反式生物堿流分中的7種生物堿相對含量作為變量，建立2個數(shù)據(jù)矩陣，進行主成分分析和聚類分析.【結果和結論】主成分分析結果表明，將紅火蟻順式生物堿流分中的15種生物堿綜合成為2個主成分因子，主要是順式－2－甲基－6－(4'－十三烯)哌啶、2－甲基－6－(6'－十五烯)－Δ1，6－哌啶、2－甲基－6－十五烷基－Δ1，6－哌啶和順式－2－甲基－6－十一烷基哌啶;反式生物堿流分中的7種生物堿綜合成1個主成分因子，主要是反式－2－甲基－6－(4'－十三烯)哌啶和反式－2－甲基－6－(6'－十五烯)哌啶.順式生物堿流分的聚類分析可以將美國種群與中國種群區(qū)分開來，而反式生物堿流分不能.

關鍵詞:紅火蟻;毒腺;生物堿;主成分分析;聚類分析

張文倩，胡林，陳立，等.基于主成分分析和聚類分析的紅火蟻毒腺生物堿的比較研究［J］.華南農(nóng)業(yè)大學學報，2015，36(6):76-81.

優(yōu)先出版時間:2015-10-16

優(yōu)先出版網(wǎng)址: http: / /www.cnki.net/kcms/detail/44.1110.s.20151016.1631.038.html

紅火蟻Solenopsis invicta Buren屬膜翅目Hymenoptera蟻科Formicidae火蟻屬Solenopsis，起源于南美洲的巴拉納河流域［1］.2004年前后于我國臺灣、廣東吳川發(fā)現(xiàn)入侵紅火蟻，短短10年間已遍布廣東省，并相繼蔓延至福建、廣西、湖南、海南、云南、重慶、香港、澳門等地［2-6］.紅火蟻攻擊性強、擴散速度快、為害范圍廣，主要危害農(nóng)作物的根、莖、葉、果實，破壞生態(tài)系統(tǒng)的多樣性，當受到干擾時會蟄刺人、畜，嚴重時會導致人休克甚至死亡［7-9］.

紅火蟻的攻擊性與其毒腺分泌物有重要關系.一般膜翅目昆蟲毒液的主要成分是水溶性蛋白，例如蜜蜂、胡蜂、黃蜂等［10］.紅火蟻毒液中水溶性蛋白極少，主要化學成分是生物堿［10］，儲藏在毒囊中并通過螯針注入被攻擊對象［11-12］.其毒液具有滅菌殺蟲等功能［13］，在與本地螞蟻競爭時具有一定的優(yōu)勢［14］，是紅火蟻能夠成功入侵并定殖的重要武器.國外對紅火蟻工蟻的毒液成分已有較多研究，其主要成分是2－甲基－6－烷基哌啶生物堿混合物，這些哌啶生物堿有2個手性中心，存在立體異構現(xiàn)象，其中，順式生物堿構型為(2R，6S)，反式構型為(2R，6R)［15］.

目前國內對紅火蟻毒腺生物堿已有一些化學成分方面的研究［16-18］，而應用數(shù)學分析法對紅火蟻毒腺生物堿進行總體的比較和綜合研究鮮見報道.本研究應用GC-MS分析了我國的廣州、惠州、深圳、珠海和東莞及美國的阿拉巴馬州6個地區(qū)的紅火蟻毒腺生物堿成分，以鑒定出的順式生物堿流分中的15種生物堿和反式生物堿流分中的7種生物堿的相對含量數(shù)據(jù)作為變量，建立2個矩陣進行主成分分析和聚類分析，找出紅火蟻毒腺生物堿中起主導作用的化學物質，并通過聚類分析指出可以用紅火蟻順式生物堿的相對含量來評價紅火蟻的分布情況，為今后研究該蟻的生態(tài)分布和入侵情況提供新思路.

1　材料與方法

1.1材料及處理

試驗所用數(shù)據(jù)來自官迪等［18］，應用GC-MS分析我國的廣州、惠州、深圳、珠海和東莞以及美國阿拉巴馬州6地24巢紅火蟻的毒腺生物堿成分，鑒定出順式生物堿流分中的15種生物堿和反式生物堿流分中的7種生物堿，各生物堿組分色譜峰積面積的相對百分比就是其相對百分含量.

1.2數(shù)據(jù)處理

以6地各個蟻巢毒腺的順式生物堿流分中的15種生物堿和反式生物堿流分中的7種生物堿相對含量數(shù)據(jù)作為變量，建立數(shù)據(jù)庫，由此構建出順式生物堿流分和反式生物堿流分2個數(shù)據(jù)集，均以蟻巢為數(shù)據(jù)矩陣的行、以生物堿相對含量為列.

利用多元統(tǒng)計分析軟件PC-ORD，對上述2個數(shù)據(jù)集進行主成分分析和聚類分析.其中聚類分析選用Euclidean距離系數(shù)，用WARD法對數(shù)據(jù)矩陣中的各行(即各蟻巢)進行聚類，輸出樹狀圖.

2　結果與分析

2.1紅火蟻的順式生物堿和反式生物堿

通過對24巢紅火蟻浸提物柱層析分析和GCMS鑒定，獲得順式生物堿流分中的15種生物堿成分和反式生物堿流分中的7種生物堿成分(圖1)，各個地區(qū)的毒腺生物堿成分沒有差異.在順式生物堿流分中，順式－2－甲基－6－(4'－十三烯)哌啶、順式－2－甲基－6－(6'－十五烯)哌啶、2－甲基－6 －(6'－十五烯)－Δ1，6哌啶和2－甲基－6－十五烷基－Δ1，6－哌啶的含量較高;在紅火蟻反式生物堿流分中，反式－2－甲基－6－(4'－十三烯)哌啶、反式－2－甲基－6－十三烷基哌啶、反式－2－甲基－6 －(6'－十五烯)哌啶和反式－2－甲基－6－十五烷基哌啶的含量較高.

圖1　紅火蟻工蟻順式和反式生物堿流分總離子色譜圖Fig.1　 Total ion chromatograms of the cis-and trans-alkaloid fractions extracted from red imported fire ant workers from different areas

2.2紅火蟻順式生物堿的主成分分析

對順式生物堿流分中的15種生物堿的主成分分析(表1)表明，前2個主成分累計方差貢獻率達90.245%，剩下的8個主成分因子對紅火蟻順式生物堿的貢獻率很小，可用這2個主成分較好地替代上述15種生物堿.

由表2可見，順式－2－甲基－6－(4'－十三烯)哌啶、2－甲基－6－十五烷基－Δ1，6－哌啶和2－甲基－6－(6'－十五烯)－Δ1，6－哌啶在第1主成分上的載荷分別為－8. 280 7、－6. 604 4、－6. 516 9，它們的載荷值符號同為負值，表明它們之間是負相關關系.主要決定第1主成分大小的為順式－2－甲基－6－(4'－十三烯)哌啶、2－甲基－6－十五烷基－Δ1，6－哌啶和2－甲基－6－(6'－十五烯)－Δ1，6－哌啶.第一主成分反映原始數(shù)據(jù)信息量的74. 392%.順式－2－甲基－6－十一烷基哌啶在第2主成分上的載荷為6. 744 8，說明其主要決定第2主成分大小，其方差貢獻率為15. 853%.

表1　順式生物堿流分的主成分分析Tab.1　 Principal component analysis of the cis-alkaloid fraction

表2　順式生物堿流分各化合物的載荷矩Tab.2 Component matrix of 15 compounds from the cis-alkaloid fraction

2.3紅火蟻反式生物堿的主成分分析

對反式生物堿流分的主成分分析(表3)表明，第1個主成分累計方差貢獻率達92.663%，可以用這一主成分較好地替代上述7種生物堿.

由表4可知，決定第1主成分的主要是反式－2－甲基－6－(4'－十三烯)哌啶和反式－2－甲基－6－(6'－十五烯)哌啶.

2.4紅火蟻生物堿的聚類分析

分別對紅火蟻的順式生物堿流分和反式生物堿流分進行聚類分析，結果如圖2.由圖2A可知，利用紅火蟻順式生物堿流分的2個主成分可以將美國采集的4巢螞蟻聚在一起，并與廣東的5個采集地區(qū)分開來.而紅火蟻反式生物堿流分的聚類結果(圖2B)較雜亂、無規(guī)律.所以，用紅火蟻的順式生物堿流分對紅火蟻的采樣地進行聚類分析更為科學.

表3　反式生物堿流分的主成分分析Tab.3　 Principal component analysis of the trans-alkaloid fraction

表4　反式生物堿流分各化合物的載荷矩Tab.4 Component matrix of 7 compounds from the trans-alkaloid fraction

圖2　紅火蟻順式和反式生物堿流分的聚類分析樹狀圖Fig.2 Dendrogram generated from clustering analysis of 15 compounds from the cis-alkaloid fraction and 7 compounds from the transalkaloid fraction of red imported fire ants

3　討論與結論

MacConnell等［11］利用氣相色譜和質譜分析鑒定了紅火蟻毒液的5種主要生物堿.Chen等［19-20］重新分析了紅火蟻工蟻毒腺的化學成分，主要包括7種順式生物堿和9種反式生物堿，除了已知的7對哌啶生物堿外，還報道了微量的哌啶環(huán)含Δ1，2和Δ1，6雙鍵的哌啶生物堿，推測其為毒液主要哌啶生物堿的生物合成途徑的前體化合物.本研究通過分析采集的24巢紅火蟻，得到相對含量較高的幾個組分，順式生物堿流分中為順式－2－甲基－6－(4'－十三烯)哌啶、順式－2－甲基－6－(6'－十五烯)哌啶、2－甲基－6－(6'－十五烯)－Δ1，6哌啶和2－甲基－6－十五烷基－Δ1，6－哌啶;反式生物堿流分中為反式－2－甲基－6－(4'－十三烯)哌啶、反式－2－甲基－6－十三烷基哌啶、反式－2－甲基－6－(6'－十五烯)哌啶和反式－2－甲基－6－十五烷基哌啶.

主成分分析是把多個指標化為少數(shù)幾個綜合指標的一種統(tǒng)計分析方法，采用一種降維的方法找出幾個綜合因子來代表原來眾多的變量［21］.本研究通過主成分分析將紅火蟻順式生物堿流分中的15種生物堿綜合成為2個主成分，分別稱為主成分1和主成分2，可以用這2個主成分較好的替代紅火蟻順式生物堿流分.在比較分析紅火蟻順式生物堿流分時，可以直接用順式－2－甲基－6－(4'－十三烯)哌啶、2－甲基－6－(6'－十五烯)－Δ1，6－哌啶、2－甲基－6－十五烷基－Δ1，6－哌啶和順式－2－甲基－6－十一烷基哌啶這4個主要化合物來進行分析判斷.而7種反式生物堿流分綜合成了1個主成分，反式－2－甲基－6－(4'－十三烯)哌啶和反式－2－甲基－6－(6'－十五烯)哌啶這2個化合物決定了這一主成分，比較分析紅火蟻反式生物堿流分時即可用這2個化合物進行分析.主成分分析的結果與紅火蟻生物堿相對含量高的化合物基本保持一致.由于紅火蟻的生物堿是2－甲基－6－烷基哌啶類混合物，利用主成分分析把紅火蟻的眾多生物堿綜合成為幾個有代表性的化合物，為紅火蟻的綜合評價提供了很好的理論依據(jù).

對順式生物堿流分進行聚類分析的結果表明，在美國采集的4巢紅火蟻聚成1類，與在中國廣東采集到的紅火蟻區(qū)分開來.中國廣東5個城市采集的20巢紅火蟻聚成3類，但分類沒有明顯規(guī)律.而反式生物堿流分聚類分析的結果不能將美國的紅火蟻與廣東地區(qū)的區(qū)分開來.Yu等［16］定量分析了紅火蟻、黑火蟻和雜合蟻的順式、反式生物堿的種間差異，其中紅火蟻的反式生物堿的含量是順式生物堿的87～248倍，遠遠高于黑火蟻和雜合蟻.中國的紅火蟻是從美國入侵過來的［22］，毒腺分泌的主要成分即反式生物堿不會有差別，而順式生物堿等微量成分有差別是可以理解的.所以，用微量成分即紅火蟻的順式生物堿流分進行聚類分析時，能將美國的紅火蟻和中國的紅火蟻區(qū)分開來，而主要成分即反式生物堿流分聚類分析不能將這2個地區(qū)的紅火蟻區(qū)分開.官迪等［18］采用Tukey-HSD多重比較的方法，分析了中國廣東和美國等6個地區(qū)紅火蟻生物堿相對含量的差異，單個成分之間的比較分析不能代表種間差異.本文則是在官迪等［18］的數(shù)據(jù)基礎上，對6個地區(qū)紅火蟻生物堿的相對含量進行主成分分析和聚類分析，把紅火蟻眾多生物堿綜合成為幾個有代表性的化合物，為紅火蟻的分布和擴散提供了理論依據(jù).此外，用生物堿對紅火蟻種群進行聚類分析，可以從化學進化的角度來解釋紅火蟻的分布.

參考文獻:

［1］TSCHINKEL W R.The fire ants［M］.Cambridge，Massachusetts: Harvard University Press，2006: 723.

［2］ZHANG R，LI Y，LIU N，et al.An overview of the red imported fire ant (Hymenoptera: Formicidae) in mainland China［J］.Fla Entomol，2007，90(4):723-731.

［3］薛大勇，李紅梅，韓紅香，等.紅火蟻在中國的分布區(qū)預測［J］.昆蟲知識，2005，42(1):57-60.

［4］曾玲，陸永躍，何曉芳，等.入侵中國大陸的紅火蟻的鑒定及發(fā)生為害調查［J］.昆蟲知識，2005，42(2):144-148.

［5］WANG L，LU Y Y，XU Y J，et al.The current status of research on Solenopsis invicta Buren (Hymenoptera: Formicidae) in mainland China［J］.Asian Myrmecol，2013 (5):125-138.

［6］陸永躍，梁廣文，曾玲.華南地區(qū)紅火蟻局域和長距離

擴散規(guī)律研究［J］.中國農(nóng)業(yè)科學，2008，41 (4):1053-1063.

［7］VINSON S B.Impact of the invasion of the imported fire ant［J］.Insect Sci，2013，20(4):439-455.

［8］VINSON S B.Invasion of the red imported fire ant (Hymenoptera: Formicidae):Spread，biology，and impact［J］.Am Entomol，1997，43(1):23-39.

［9］陸永躍.中國大陸紅火蟻遠距離傳播速度探討和趨勢預測［J］.廣東農(nóng)業(yè)科學，2014，41(10):70-72.

［10］JONES T H，BLUM M S，F(xiàn)ALES H M.Ant venom alkaloids from Solenopsis and Monorium species: Recent developments［J］.Tetrahedron，1982，38(13):1949-1958.

［11］MACCONNELL J G，BLUM M S，F(xiàn)ALES H M.Alkaloid from fire ant venom: Identification and synthesis［J］.Science，1970，168(3933):840-841.

［12］MACCONNELL J G，BLUM M S，F(xiàn)ALES H M.The chemistry of fire ant venom［J］.Tetrahedron，1971，27 (6):1129-1139.

［13］BLUM M S，WALKER J R，CALLAHAN P S，et al.Chemical，insecticidal and antibiotic properties of fire ant venom［J］.Science，1958，128(3319):306-307.

［14］MORRISON L W.Long-term impacts of an arthropod-community invasion by the imported fire ant，Solenopsis invicta ［J］.Ecology，2002，83(8):2337-2345.

［15］LECLERCQ S，THIRIONET I，BROEDERS F，et al.Absolute configuration of the solenopsins，venom alkaloids of the fire ants［J］.Tetrahedron，1994，50 (28):8465-8478.

［16］YU Y，WEI H，F(xiàn)ADAMIRO H，et al.Quantitative analysis of alkaloidal constituents in imported fire ants by gas chromatography［J］.J Agr Food Chem，2014，62(25):5907-5915.

［17］馬伏寧，曾鑫年，潘達強.多后型紅火蟻工蟻和蟻后毒液生物堿成分的比較［J］.中國農(nóng)學通報，2009，25 (17):57-61.

［18］官迪，廖曉蘭，陳立.中國和美國紅火蟻毒腺生物堿組分的比較分析［J］.昆蟲學報，2013，56(4):365-371.

［19］CHEN L，F(xiàn)ADAMIRO H Y.Re-investigation of venom chemistry of Solenopsis fire ants: I: Identification of novel alkaloids in S.richteri［J］.Toxicon，2009，53(5):469-478.

［20］CHEN L，F(xiàn)ADAMIRO H Y.Re-investigation of venom chemistry of Solenopsis fire ants: II: Identification of novel alkaloids in S.invicta［J］.Toxicon，2009，53(5):479-486.

［21］王芳.主成分分析與因子分析的異同比較及應用［J］.統(tǒng)計教育，2003，5(5):14-17.

［22］ASCUNCE M S，YANG C，OAKEY J，et al.Global invasion history of the fire ant Solenopsis invicta［J］.Science，2011，331(6020):1066-1068.

【責任編輯霍歡】

Comparative studies on alkaloids secreted by vemon gland of Solenopsis invicta based on principal component analysis and cluster analysis

ZHANG Wenqian1，2，HU Lin2，CHEN Li2，ZENG Rensen1
(1 College of Agriculture，South China Agricultural University，Guangzhou 510642，China; 2 State Key Laboratory of Integrated Management of Pest Insects and Rodents/Institute of
Zoology，Chinese Academy of Sciences，Beijing 100101，China)

Abstract:【Objective】To analyze the major components and cluster features of alkaloids，which are secreted by vemon gland of red imported fire ants，Solenopsis invicta，from China and USA.【Method】Twenty four colonies of S.invicta were collected from Guangzhou，Huizhou，Shenzhen，Zhuhai，Dongguan of China and Alabama of USA.Alkaloids of vemon gland were extracted.Two data matrix were established with variable of relative contents of 15 kinds of cis-alkaloids and 7 kinds of trans-alkaloid respectively.The alkaloid fraction was analyzed by principal component analysis and cluster analysis.【Result and conclusion】The results showed that the 15 types of alkaloids from the cis-alkaloid fraction werebook=77,ebook=349consolidated into two principal components，including cis-2-methyl-6-(4'-tridecenyl) piperidine，2-methyl-6-(6'-pentadecenyl) -Δ1，6-piperidine，2-methyl-6-pentadecyl-Δ1，6-piperidine and cis-2-methyl-6-undecylpiperidine seven types of alkaloids from the trans-alkaloid fraction were consolidated into one principal component，including trans-2-methyl-6-(4'-tridecenyl) piperidine and trans-2-methyl-6-(6'-pentadecenyl) piperidine.The results of cluster analyses demonstrate that alkaloids from the cis-alkaloid fraction can distinguish the USA populations from the Chinese populations，but alkaloids from the trans-alkaloid fraction can not.

Key words:Solenopsis invicta; vemon gland; alkaloid; principal component analysis; cluster analysis

基金項目:國家自然科學基金(30970402)

作者簡介:張文倩(1990—)，女，碩士研究生，E-mail: woshizwq@ 126.com;通信作者:曾任森(1965—)，男，教授，博士，E-mail: rszeng@ scau.edu.cn

收稿日期:2014-12-23

文章編號:1001-411X(2015) 06-0076-06

文獻標志碼:A

中圖分類號:S481