蟬花居群形態變異式樣

李建平，張 鐵，曾文波*

（1.文山學院環境與資源學院，云南 文山 663099；2.文山學院文山州生物資源開發研究中心，云南 文山 663099）

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

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

目前，關于蟲草類真菌的形態特征研究多是關于蟲草子座、菌核、子囊結構[6-12]，以及菌絲體形態[4]的描述，關于蟲草居群形態學的報道很少，本文以居群為單位，研究蟬花的孢梗束、菌核形態在居群間和居群內及個體間形態變異式樣。從形態角度探討蟬花狀變異式樣與生態適應特征，揭示其分化特點，為蟬花資源的開發利用提供基礎數據。

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 形態性狀測定方法

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

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

1.3 數據分析方法

1.3.1 形態性狀基本統計量

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

1.3.2 形態性狀巢式方差分析

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

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

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

1.3.3 形態性狀Q-聚類分析

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

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

2 結果與分析

2.1 形態特征基本統計量

變異系數可表示表型性狀值離散性的程度，變異系數越大，則性狀值離散程度越大，反之則小。蟬花和冬蟲夏草物種水平形態性狀統計結果見表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的寄主是云南黑蟬，表明蟬花的菌核長和周長主要是與寄主直接相關的，菌核長和菌核周長這兩個形態性狀的CV值較小，說明這兩個性狀穩定，變異不大。蟬花菌核重量（JZ）平均值在 0.20 g～0.95 g 之間，CV 值為 17.35%～36.63%；孢梗束重量 （ZZ）平均值在 0.09 g～0.44 g之間，CV值為13.40%～45.79%，表明不同產地的蟬花的菌核和孢梗束重量有較大的差異。蟬花孢梗束一級分枝數（BR）平均值在1.7枝～3.9枝之間，CV值為38.83%～114.32%，表明蟬花孢梗束多分枝，孢梗束的分枝數變異大。分枝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個形態性狀的變異系數（物種水平），分枝3周長（BZ－3）最大，CV值為163.79%，其次為分枝3長度（BC－3），CV值為162.57%，由于個體間分枝數目不一致，使得分枝相關的形態性狀變異最大。

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

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

2.2 巢式方差分析

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

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

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

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

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

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

2.3 Q-聚類分析

10個蟬花居群17個形態性狀AVE聚類法的譜系聚類圖見圖1。

圖1 10個蟬花居群17個形態性狀AVE聚類法的譜系聚類圖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中，產自蘭坪縣的7個居群（Pop ICLHH、I CLHQ、ICLZD、ICLSH、 ICLTL、 ICLHX、ICLTD）先聚成一支，再依次與Pop ICWYJ和Pop ICYTS聚在一起，聚類圖在一定程度上與地理呈相關性，表明不同居群的蟬花的形態性狀具有地域性特征。

3 討論

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

參考文獻：

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

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

[3]李冰嵐.蟬花的本草學考證[J].現代應用藥學，1993，10（2）：21-22.

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

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

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

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

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

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

[10]劉作易.蟲草屬及其無性型關系研究[D].武漢：華中農業大學，1999.

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

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

[13]杜榮騫.生物統計學[M].北京：高等教育出版社，1990.

[14]高惠璇.實用統計方法與SAS系統[M].北京：北京大學出版社，2001.

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

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

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

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

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

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