李建平，張 鐵，曾文波*

（1.文山學(xué)院環(huán)境與資源學(xué)院，云南 文山 663099；2.文山學(xué)院文山州生物資源開發(fā)研究中心，云南 文山 663099）

蟲草類真菌Cordyceps s.l.感染寄主后形成蟲菌復(fù)合體，包括大多數(shù)蟲草、菌菌復(fù)合體（如大團囊蟲草、頭狀蟲草和分枝蟲草等）[1]，以及重寄生（多頭霉等可寄生于其他蟲草菌上）[2]，蟲草類真菌形態(tài)獨特，昆蟲等被真菌感染后僵死，真菌利用寄主體內(nèi)營養(yǎng)長成子座或孢梗束突出于寄主表面，完整的子實體包括蟲體部分（菌核）和一段似草的部分子座（或孢梗束），因此人們形象地將之稱為蟲草。

“蟬花”（Isaria cicadae）是我國名貴傳統(tǒng)中藥材，最早關(guān)于“蟬花”的記載是成書于公元五世紀的《雷公炮炙論》[3]，現(xiàn)代研究表明，蟬花為真菌Isaria cicadae Miquel感染蟬若蟲后形成的蟲菌復(fù)合體[4]，分布廣泛，包括中國、韓國、日本、澳大利亞、新西蘭、巴西、南亞、歐洲、北美等國家和地區(qū)[4－5]，具有多種藥理作用，在我國的藥用歷史超過1 500年，可作為冬蟲夏草的代用品[3－4]，蟬花是具有巨大開發(fā)價值的野生資源之一。

目前，關(guān)于蟲草類真菌的形態(tài)特征研究多是關(guān)于蟲草子座、菌核、子囊結(jié)構(gòu)[6-12]，以及菌絲體形態(tài)[4]的描述，關(guān)于蟲草居群形態(tài)學(xué)的報道很少，本文以居群為單位，研究蟬花的孢梗束、菌核形態(tài)在居群間和居群內(nèi)及個體間形態(tài)變異式樣。從形態(tài)角度探討蟬花狀變異式樣與生態(tài)適應(yīng)特征，揭示其分化特點，為蟬花資源的開發(fā)利用提供基礎(chǔ)數(shù)據(jù)。

1 材料和測定方法

1.1 材料

用于觀察的材料為新鮮蟲草，包括蟬花10個居群共100個個體，冬蟲夏草（Ophiocordyceps sinensis）5個居群共50個個體，詳見表1。

表1 10個居群蟬花和5個居群冬蟲夏草來源Tab.1 Locality of 10 populations of Isaria cicadae and 5 populations of Ophiocordyceps sinensis

1.2 形態(tài)性狀測定方法

采用生物統(tǒng)計方法[13]，每個居群隨機選取新鮮樣品10個個體，按蟲草不同部位（菌核部分，孢梗束或子座部分），分別測量其重量、周長、長、分枝數(shù)等性狀，其中，重量為樣品干燥后測定值，每個性狀測量3次，取平均值，所測性狀及其對應(yīng)的編號見表2，其中，由于蟬花孢梗束分枝較多，選取一級分枝最長的3枝測量，對于達不到3枝的個體，缺失分枝長BC和分枝周長BZ兩個性狀數(shù)值記為“0”。

表2 蟬花和冬蟲夏草觀察性狀、編號Tab.2 Number of morphological traits of Isaria cicadae and Ophiocordyceps sinensis

1.3 數(shù)據(jù)分析方法

1.3.1 形態(tài)性狀基本統(tǒng)計量

使用SAS 8.1軟件計算蟬花和冬蟲夏草形態(tài)性狀的平均值、標準差、最大值、最小值和變異系數(shù)。標準差StdDev表示形態(tài)性狀值在平均值上下波動大小，標準差越大則反映形態(tài)性狀值離平均值波動越遠；變異系數(shù)CV（標準差除以平均值）表示表型性狀值離散性程度，變異系數(shù)越大則表型性狀值離散程度越大，反映形態(tài)性狀值在居群中的整齊性程度[14]。

1.3.2 形態(tài)性狀巢式方差分析

形態(tài)性狀巢式等級方差分析[15-18]的模型為：

式中：Yijkl為第i個居群第j個單株第k次重復(fù)第l次的觀測值；μ為總體平均值；Si為第i個居群的效應(yīng)（固定）；T(i)j為第i個居群第j個單株的效應(yīng)（隨機）；R(ij)k為第i個居群第j個單株第k次重復(fù)的效應(yīng)（固定）；ε(ijk)l=ε(ijk)=試驗誤差。

為進一步了解蟬花形態(tài)性狀在居群間和個體間的分配情況，利用SAS 8.1進行巢式方差分析，使用的SAS程序[18]如下：

1.3.3 形態(tài)性狀Q-聚類分析

采用AVE聚類法，對10個居群蟬花樣品形態(tài)性狀的平均值進行居群間聚類分析，使用的SAS程序[18]如下：

表3 蟬花與冬蟲夏草物種水平形態(tài)性狀分析結(jié)果Tab.3 The analysis results of morphologic traits of Isaria cicadae and Ophiocordyceps sinensis at species level

2 結(jié)果與分析

2.1 形態(tài)特征基本統(tǒng)計量

變異系數(shù)可表示表型性狀值離散性的程度，變異系數(shù)越大，則性狀值離散程度越大，反之則小。蟬花和冬蟲夏草物種水平形態(tài)性狀統(tǒng)計結(jié)果見表3。

10個居群的蟬花，其菌核長（JC）平均值在3.18 cm～5.15 cm 之間，居群水平 CV 值為 5.93%～10.56%，物種水平 CV值為 14.91%；菌核周長（ZC）平均值在2.98 cm～4.16 cm之間，居群水平CV值為 6.07%～9.97%，物種水平 CV 值為 12.97%，Pop ICKSG的菌核長和周長最小，Pop ICYTS的菌核長和周長最大，Pop ICYTS的寄主為竹蟬，Pop ICKSG的寄主是鳴鳴蟬，Pop ICWYJ的寄主是云南黑蟬，表明蟬花的菌核長和周長主要是與寄主直接相關(guān)的，菌核長和菌核周長這兩個形態(tài)性狀的CV值較小，說明這兩個性狀穩(wěn)定，變異不大。蟬花菌核重量（JZ）平均值在 0.20 g～0.95 g 之間，CV 值為 17.35%～36.63%；孢梗束重量 （ZZ）平均值在 0.09 g～0.44 g之間，CV值為13.40%～45.79%，表明不同產(chǎn)地的蟬花的菌核和孢梗束重量有較大的差異。蟬花孢梗束一級分枝數(shù)（BR）平均值在1.7枝～3.9枝之間，CV值為38.83%～114.32%，表明蟬花孢梗束多分枝，孢梗束的分枝數(shù)變異大。分枝1長（BC－1）平均值在3.05 cm～7.19 cm 之間，分枝 1 周長 （BZ－1）平均值在 1.26 cm～2.31 cm 之間；分枝 2 長 （BC－2）平均值在 1.91 cm～5.36 cm 之間，分枝 2 周長 （BZ－2）平均值在 0.38 cm～0.90 cm 之間；分枝 3 長 （BC－3）平均值在 0.50 cm～2.19 cm 之間，分枝 3周長 （BZ－3）平均值在0.06 cm～0.38 cm之間，蟬花12個形態(tài)性狀的變異系數(shù)（物種水平），分枝3周長（BZ－3）最大，CV值為163.79%，其次為分枝3長度（BC－3），CV值為162.57%，由于個體間分枝數(shù)目不一致，使得分枝相關(guān)的形態(tài)性狀變異最大。

除了與孢梗束分枝數(shù)相關(guān)性狀以外，蟬花孢梗束重量的變異系數(shù)最大（45.88%）；其次是子實體重量QZ（44.29%）；不同居群內(nèi)每一表型性狀的變異系數(shù)均有一定差異，說明蟬花居群間由于遺傳和環(huán)境異質(zhì)性導(dǎo)致居群表型變異存在差異。所研究10個蟬花居群的12個形態(tài)性狀中，與孢梗束分枝數(shù)相關(guān)的形態(tài)性狀變異系數(shù)值均偏大，說明這些性狀在居群內(nèi)個體間分化較大，多態(tài)性豐富，這可能與其適應(yīng)生長環(huán)境有關(guān)。

冬蟲夏草5個居群8個形態(tài)性狀的變異系數(shù)（物種水平）在16.02%～48.44%之間。這些形態(tài)性狀中菌核重量的變異系數(shù)最大（48.44%），不同產(chǎn)地冬蟲夏草的大小和重量等形態(tài)性狀，主要與冬蟲夏草寄主的種類和被感染年齡有關(guān)，子座一級分枝數(shù)（BR）的變異系數(shù)為0，表明該性狀較其他性狀穩(wěn)定，冬蟲夏草子座多單生。

2.2 巢式方差分析

巢式方差分析可以揭示居群間和居群內(nèi)的變異，而且可以劃分居群間和居群內(nèi)變異組分。蟬花形態(tài)性狀的表型分化系數(shù)（Vst）見表4。

表4 蟬花10個居群12個性狀巢式方差分析Tab.4 Nested variation analysis of 12 morphological traits among 10 populations of Isaria cicadae

從表4可知，蟬花菌核部分12個形態(tài)性狀在居群間的平均方差分量占33.79%，居群內(nèi)個體間的平均方差分量占總變異的66.21%，居群內(nèi)個體是形態(tài)性狀變異的主要來源。其中，菌核部分3個形態(tài)性狀在居群間的平均方差分量占69.21%，居群內(nèi)個體間的平均方差分量占總變異的30.79%，居群間的多態(tài)性大于居群內(nèi)的多態(tài)性；蟬花孢梗束部分7個形態(tài)性狀在居群間的平均方差分量占18.02%，居群內(nèi)個體間的平均方差分量占總變異的81.98%，居群內(nèi)的多態(tài)性遠大于居群間的多態(tài)性，與孢梗束分枝等相關(guān)的形態(tài)性狀的差異則主要來自于個體間，同一居群不同個體的孢梗束形態(tài)性狀差異巨大，孢梗束的生長除受到蟬花本身的遺傳因素影響，在生長過程中環(huán)境因素的影響較大，這使得居群內(nèi)個體間的孢梗束形態(tài)變異較菌核相關(guān)形態(tài)性狀豐富。

冬蟲夏草12個性狀中，由于測量個體的分枝數(shù)均為1枝，有效的形態(tài)性狀數(shù)據(jù)有7個，詳見表5。

表5 冬蟲夏草5個居群7個性狀巢式方差分析Tab.5 Nested variation analysis of 7 morphological traits among 5 populations of Ophiocordyceps sinensis

從表5可以看出，這七個性狀均表現(xiàn)為居群間方差分量大于個體間方差分量，表明居群是這些形態(tài)性狀的主要來源，居群間形態(tài)有較大的分化。冬蟲夏草形態(tài)性狀在居群間與居群內(nèi)的分布式樣與蟬花不一樣，冬蟲夏草形態(tài)性狀在居群間的分化大于蟬花居群間的分化。

2.3 Q-聚類分析

10個蟬花居群17個形態(tài)性狀A(yù)VE聚類法的譜系聚類圖見圖1。

圖1 10個蟬花居群17個形態(tài)性狀A(yù)VE聚類法的譜系聚類圖Fig.1 Q-Cluster of 17 morphological traits assayed in 10 populations of Isaria cicadae with average linkage method

由圖1可知，10個居群的蟬花聚成CladeA和CladeB兩支，在CladeB中包括9個居群的蟬花，CladeA中包括1個居群（Pop ICKSG）的蟬花。在CladeB中，產(chǎn)自蘭坪縣的7個居群（Pop ICLHH、I CLHQ、ICLZD、ICLSH、 ICLTL、 ICLHX、ICLTD）先聚成一支，再依次與Pop ICWYJ和Pop ICYTS聚在一起，聚類圖在一定程度上與地理呈相關(guān)性，表明不同居群的蟬花的形態(tài)性狀具有地域性特征。

3 討論

蟲草菌核的大小和重量等性狀主要是受到寄主大小的影響，蟬花的寄主為同翅目的昆蟲蟬，現(xiàn)已報道關(guān)于蟬花的寄主蟬共有8種[19]，而冬蟲夏草的寄主為蝙蝠蛾，已報道的種類有57種[20]，冬蟲夏草寄主的種類遠遠多于蟬花的寄主，這可能是冬蟲夏草形態(tài)性狀居群間分化大于蟬花的原因之一。另外一個可能的原因是：人們在巨大的經(jīng)濟利益驅(qū)動下，對冬蟲夏草進行破壞性的采集，導(dǎo)致冬蟲夏草居群分布點的范圍不斷縮小，這阻礙了冬蟲夏草居群間的基因交流，使得冬蟲夏草居群間分化進一步加大。由于寄主的種類直接決定蟬花和冬蟲夏草的菌核大小，與菌核大小直接相關(guān)的性狀，如菌核長（JC）、菌核周長（ZC）、菌核重量（JZ）等性狀的變異要?。–V值較?。艿讲煌a(chǎn)地的寄主種類和分化的影響，使得菌核大小直接相關(guān)的性狀在居群間形成分化。與孢梗束相關(guān)的形態(tài)性狀，如孢梗束一級分枝數(shù)（BR）、分枝1長（BC-1）、分枝 2長 （BC-2）、分枝3長（BC-3）、分枝1周長（BZ-1）、分枝3周長（BZ-3）、分枝3周長（BZ-3）則變異較大。同一居群內(nèi)不同個體的孢梗束差異較大，可能是由于不同的個體在遺傳上的差異導(dǎo)致的，在培養(yǎng)蟬棒束孢菌絲體的過程中還發(fā)現(xiàn)，不同的菌株在相同的培養(yǎng)條件下，其產(chǎn)孢子、產(chǎn)菌絲、產(chǎn)孢梗束的能力差異明顯，表明其在遺傳結(jié)構(gòu)上可能具有較大的差異，有必要開展蟬花遺傳結(jié)構(gòu)的研究；另外，蟬花孢梗束的形成還受到其生長環(huán)境的影響，如菌核距地表的距離、地表的枯枝落葉、地表植被等。總之，遺傳上的差異和生境的多樣性是蟬花孢梗束形態(tài)變異豐富的來源。豐富的變異一定程度上也反映了蟬花具有較強的適應(yīng)環(huán)境能力，這與蟬花分布廣、適應(yīng)性強[4]的實際情況是一致的。

參考文獻：

[1]梁宗琦.中國真菌志（第三十二卷蟲草屬）[M].北京：科學(xué)出版社，2007.

[2]王元兵.昆明地區(qū)下垂蟲草生物學(xué)研究[D].昆明：云南大學(xué)，2012.

[3]李冰嵐.蟬花的本草學(xué)考證[J].現(xiàn)代應(yīng)用藥學(xué)，1993，10（2）：21-22.

[4]曾文波，常襯心，李建平，等.蟬棒束孢顯微形態(tài)變異樣[J].微生物學(xué)報，2017，57（3）：350-362.

[5]曾文波，王元兵，虞泓，等.蟬棒束孢居群分布及生境調(diào)查[J].中國食用菌，2015，34（4）：50-56.

[6]劉愛英，梁宗琦.中國蟲草屬一新記錄種[J].真菌學(xué)報，1993，12 （3）：248-250.

[7]梁宗琦，劉愛英，劉作易.冬蟲夏草子囊孢子的發(fā)育[J].真菌學(xué)報，1995，14（2）：148-152.

[8]梁宗琦，劉愛英，黃建忠，等.寬闊水自然保護區(qū)的蟲草及其相關(guān)真菌I[J].真菌學(xué)報，1996，15（4）：264-271.

[9]梁宗琦，劉愛英，黃建忠，等.寬闊水自然保護區(qū)的蟲草及其相關(guān)真菌II[J].菌物系統(tǒng)，1997，16（1）：61-67.

[10]劉作易.蟲草屬及其無性型關(guān)系研究[D].武漢：華中農(nóng)業(yè)大學(xué)，1999.

[11]臧穆.我國熱帶真菌兩新種[J].云南植物研究，2001，23（3）：295-297.

[12]楊鐘林，虞泓，陳自宏，等.滇中蛹蟲草居群生物學(xué)及其生態(tài)學(xué)習(xí)性研究[J].中國食用菌，2011，30（5）：43-47.

[13]杜榮騫.生物統(tǒng)計學(xué)[M].北京：高等教育出版社，1990.

[14]高惠璇.實用統(tǒng)計方法與SAS系統(tǒng)[M].北京：北京大學(xué)出版社，2001.

[15]葛頌，洪德元，錢迎倩.遺傳多樣性及其檢測方法，生物多樣性研究的原理與方法 [M].北京：中國科學(xué)技術(shù)出版社，1994.

[16]虞泓，黃瑞復(fù).云南松居群核型變異及其分化研究[J].植物分類學(xué)報，1998，36（3）：222-231.

[17]虞泓，楊彩云，徐正堯.云南松居群花粉形態(tài)多態(tài)性[J].云南大學(xué)學(xué)報：自然科學(xué)版，1999，21（2）：86-89.

[18]何顯靜，王兵益，虞泓.SAS系統(tǒng)在川百合形態(tài)變異式樣研究中的運用[J].云南大學(xué)學(xué)報：自然科學(xué)版，2003，25（增刊）：75-77.

[19]陳祝安，李增智，陳以平.金蟬花[M].北京：中醫(yī)古籍出版社，2014.

[20]Wang XL,Yao YJ.Host insect species of Ophiocordyceps sinensis:a review[J].ZooKeys,2011(127):43-59.