湖北省煙草赤星病長柄鏈格孢遺傳多樣性分析

楊濤，黎妍妍，鄭露，黃俊斌，黃凱，郭利，周鵬，袁躍斌，彭五星，李錫宏

1 華中農業大學植物科技學院，湖北省作物病害監測和安全控制重點實驗室，武漢430070；

2 湖北省中煙工業有限責任公司恩施卷煙廠，湖北 恩施 445000；

3 湖北省煙草科學研究院，武漢 430030；

4 湖北省煙草公司十堰市公司，湖北 十堰 442000；

5 湖北省煙草公司襄陽市公司，湖北 襄陽 441100；

6 湖北省煙草公司宜昌市公司，湖北 宜昌 443000；

7 湖北省煙草公司恩施州公司，湖北 恩施 445500

煙草赤星病是危害全世界煙草生產最主要的真菌性病害之一，其病原菌屬于半知菌類鏈格孢屬[1]。我國于1916年首次在北京近郊發現該病，80年代以來，全國煙區有日益嚴重的趨勢，尤以黑龍江省、吉林兩省危害最重。1990年以后，該病在我國迅速發展成為煙草主要病害之一，對煙葉產量、產值有顯著影響[2-5]，煙草赤星病發生嚴重時對煙葉產量的損失達到30%，而對產值損失率最高可以達到90%[6]，據統計，2002年煙草赤星病發生面積2萬公頃，產量損失約2萬公斤，產值損失達到1.8億元[7]。目前國內已報道的引起煙草赤星病的鏈格孢種類有A.alternata、A. longipes、A. tenuissim、A. yaliinficiens和A. tabacicola，其 中A. alternata和A. longipes可引起云南[8]和貴州煙草赤星病[9]，A. tenuissima可引起重慶煙草赤星病[10]，河南省發現的煙草赤星病病原菌包含了已知的5個種[11]，而湖北省煙草赤星病菌由A. tenuissima、A. alternata、A. yaliinficiens和A.longipes4個種引起，其中長柄鏈格孢A. longipes是優勢致病菌，不僅在湖北省各個煙區都有分布，且不同地方的菌株致病力存在差異[12]。

ISSR是Zietkiewicz在SSR基礎上發展起來的一種分子標記技術[13]。該技術已廣泛用于生物遺傳連鎖圖譜的構建、基因定位、種質資源鑒定、分類、進化和遺傳多樣性分析等方面的研究[14-15]。王葉[16]采用ISSR技術對引起棗果病害的鏈格孢屬病原真菌進行遺傳多樣性分析，結果表明在地理區域劃分中濮陽地區和山東地區的菌株能夠大致的分在獨立的組，說明ISSR技術可以用做對鏈格孢病原菌的遺傳多樣性研究。胡中會對赤星病ISSR-PCR正交體系進行了優化[17]，并采用ISSR和RAPD兩種分子標記方法對種群的遺傳多樣性與地理來源的關系進行了分析，結果顯示菌株的遺傳多樣性和地理來源并無明顯的相關性[18]。祖艷青[11]采用ISSR技術對河南省煙草赤星病 鑒 定 為A.alternata、A.longipes、A.tenuissima及A.yaliinficiens的4個種也進行了遺傳多樣性和地理來源的相關性分析，與胡中會研究得到的結果相似。

本研究擬通過ISSR分子標記的方法，以湖北省十堰、襄陽、恩施、宜昌等4個產煙區的優勢致病種—長柄鏈格孢菌為研究對象，對其遺傳多樣性、遺傳結構與分化以及種群間遺傳變異進行了多方位的分析，以期為煙草赤星病的防治奠定理論基礎。

1 材料和方法

1.1 菌株來源

采用單孢分離法，并經過形態學和分子鑒定從煙草赤星病病斑上獲得了34株長柄鏈格孢菌株，菌株編號和信息見表1。

表1 供試菌株來源以及寄主信息Tab.1 Source of test isolate and host information

1.2 基因組DNA提取及檢測

將供試菌株轉入鋪有滅菌玻璃紙的PDA培養基平板上培養，7 d后收集菌絲。采用Huang等[19]CTAB法提取菌株基因組DNA，用0.8%(m/V)的瓊脂糖凝膠電泳檢測 DNA 的完整性，用德國IMPLENP300核酸蛋白超微量紫外分光光度計檢測DNA的純度和濃度，A260/A280在1.90以上，并將DNA用TE稀釋至終濃度為 50 ng·μL-1，放入 -20℃冰箱保存。

1.3 引物的篩選及ISSR-PCR反應體系

所用ISSR引物參考Zietkiewicz[13]、王葉[16]的研究。用擎科生物科技有限公司合成的32條引物進行初步篩選，選取8條多態性好、譜帶清晰、穩定可靠的引物(表2)用于對菌株的遺傳關系研究。采用25 μL的體系進行擴增，10 μmol/L的上下游引物共2 μL，ddH2O 9.5 μL，20 ng/μL 的 DNA 模板 1μL，反應用的PCR Mix12.5 μL。擴增程序為94℃預變性5min；接著進行35個循環：94℃變性30 s，復性45 s，72℃延伸1.5 min；循環結束后72℃延伸10 min。

1.4 遺傳多樣性數據統計與分析

根據ISSR-PCR條帶的有無，將統計結果構建0，1矩陣，并保存為Excel文件格式。利用POPGENE（version 3.2）軟件計算以下參數。

反映基因多少和狀況的遺傳參數：多態性位點百分率P（Proportion of polymorphic loci），等位基因觀測數Na（Observed number of alleles）和有效等位基因數 Ne（Effective number of alleles）。

反映基因多樣性信息的遺傳參數：Shannon 信息指數 I（Shannons information index）[20]，表型多樣性指數[21]，Nei’s基因多樣性指數H（Gene diversity within the species）[20]。

反映遺傳分化程度的參數：總的遺傳多樣性Ht（Gene diversity），群體內的遺傳多樣性Hs（Gene diversity within population），群體間遺傳多樣性Dst（Gene diversity among population），Dst = Ht﹣Hs， 遺 傳分化系數Gst（Coefficient of gene differentiation among population within species），Gst = Dst/Ht。居群每代遷移數Nm（number of migrants per generation）。

反映基因差異的參數：不同群體間的遺傳相似性I（Genetic identity）和遺傳距離D（Genetic distance）。對遺傳距離進行聚類分析并通過Tree plot模塊生成聚類圖。

采用NTSYSpc-2.10軟件SHAN程序中的UPGMA（Unweighted pair-group mean average，非加權算術平均聚類），按照遺傳相似性對群體進行聚類分析。

2 結果與分析

2.1 ISSR-PCR引物擴增多態性

用篩選得到的8條重復性好、多態性豐富的引物擴增34個菌株，部分引物的篩選結果如下圖所示（圖1），引物的擴增結果如圖（圖2）；由表2統計可以看出，8條ISSR引物擴增34個菌株共擴增出116個位點，其中多態性位點99個，多態位點百分率（P）為85.34%。引物UBC810的擴增位點最多，為18個；引物UBC811的擴增位點最少，為12個。

圖1 ISSR引物UBC810、UBC826對部分供試長柄鏈格孢菌株的電泳圖Fig.1 Electropherograms of ISSR primers UBC810, UBC826 on some tested strains of A. longipes

圖2 ISSR引物UBC811擴增菌株的電泳圖Fig.2 Electropherograms of ISSR primers UBC811 on some tested strains of A.longipes

表2 ISSR引物的多態性Tab.2 Polymorphism of ISSR primers

2.2 長柄鏈格孢菌菌株的遺傳多樣性分析

2.2.1 長柄鏈格孢菌不同地理來源菌株的遺傳多樣性由軟件POPGENE（version 3.2）計算得出不同種群菌株的平均多態位點數為55，Nei’s 多樣性指數為0.1118 ～0.2061，Shnanon 信息指數為 0.1972～0.3177；其中十堰菌株的Nei’s多樣性指數為0.2061，Shnanon信息指數為0.3177，為4個地理種群中最高；而宜昌菌株的Nei’s多樣性指數為0.1118，Shnanon 信息指數為0.1689，相對其它地理種群較低（表3）。上述結果表明長柄鏈格孢菌不同地理種群之間有豐富的遺傳多樣性。

表3 基于ISSR標記的長柄鏈格孢菌遺傳多樣性Tab.3 Genetic diversity of A. longipes based on ISSR marker

2.2.2 長柄鏈格孢菌不同地理種群的遺傳距離分析

湖北省煙草赤星病長病鏈格孢菌總基因遺傳多樣性為0.1815，種群內遺傳多樣性為 0.1605，種群間遺傳分化系數為0.1305，居群每代遷移數（Nm）為3.2740（表4）。通過對4個不同產煙區種群的基因分化系數（Gst）分析，種群內遺傳多樣性變異占85.34%，表明菌株的遺傳變異主要存在種群內。4個不同地理種群的居群每代遷移數（Nm）為3.2740＞1，表明在不同地理來源的種群間存在廣泛的基因流動，從而阻礙了由遺傳漂變所導致的遺傳分化。

表4 長柄鏈格孢菌株基因多樣性的 Nei’s 分析Tab.4 Nei’s analysis of genetic diversity of Alternaria longipes strains

種群間的遺傳距離結果見表5。4個地理來源種群的遺傳距離平均值為0.03915，遺傳相似性平均值為0.9616。十堰和恩施的種群遺傳距離最大為0.0503，而襄陽與十堰種群的遺傳距離最小為0.0282。

表5 長柄鏈格孢菌種群間Nei’s遺傳相似性（I，上三角）及遺傳距離（D，下三角）Tab.5 Nei’s genetic similarity (I, upper triangle) and genetic distance (D, lower triangle) between populations of A. longipes

根據種群間遺傳距離結果，按照非加權算術平均（UPGMA）法繪制4個不同種群間聚類分析樹（圖3），4個種群可以分為兩類，十堰和襄陽的聚為一類，恩施和宜昌的聚為一類。

圖3 長柄鏈格孢菌株病菌種群ISSR聚類分析圖Fig.3 ISSR clustering analysis of A. longipes strains

2.2.3 煙草長柄鏈格孢菌（A. longipes）ISSR聚類分析結果

由NTSYSpc-2.10軟件構建UPGAM的聚類分析樹（圖4）可以看出，供試的34株長柄鏈格孢菌株具有較高的遺傳相似性，其相似系數為0.52～0.99。當相似性系數設定為0.79時，可將34個供試菌株分為4個類群：第1類群包括襄陽?？档?個菌株XYBK-910-2；第2個類群包括了31個菌株，是優勢類群，菌株來自襄陽、十堰、恩施和宜昌4個煙葉產區；第3個類群只有來自恩施利川市的1個菌株ESLC-1410-1；第4個類群是來自十堰市鄖西縣的1個菌株SY-YX-920-1。34株供試長柄鏈格孢菌株并沒有以不同地理來源完全聚為不同的組。

圖4 34個長柄鏈格孢菌株的ISSR樹狀聚類圖Fig.4 UPGMA dendrogram of 34 A. longipes strains

3 討論

國內外利用ISSR分子標記技術對鏈格孢菌的遺傳多樣性進行了研究[22-24]，其中祖艷青[11]和胡中會[18]采用該方法研究了赤星病菌遺傳多樣性和地理來源的關系，但并沒有對種群的遺傳分化和種群的遺傳變異做進一步分析，李六英[25]研究了鏈格孢菌和長柄鏈格孢菌的遺傳多樣性及遺傳分化，結果發現鏈格孢菌的遺傳多樣性與地理來源有明顯的相關性，但長柄鏈格孢菌的研究結果相反，其初步分析了2個種的遺傳分化特點，但未對遺傳分化結果和種群的遺傳變異做深入研究。本研究從遺傳多樣性、遺傳結構與分化和種群的遺傳變異多方面對湖北省煙草赤星病優勢致病種進行了分析。

本研究對湖北省4個煙區34個長柄鏈格孢菌的遺傳多樣性參數進行統計分析，結果表明菌株遺傳多樣性處于相對較高的水平，但是不同地理種群間遺傳多樣性水平不一樣，差異并不顯著。另外煙草種植環境和品種的選擇往往能夠誘導赤星病致病菌等位基因的變異，導致遺傳結構的分化，其分化的程度利用遺傳分化指數、遺傳距離等進行度量，數值越大，表明群體間的遺傳分化水平越高[26-28]，本研究中，4個地理來源種群的遺傳分化系數為0.1305，遺傳距離平均值為0.03915，表明種群間產生了一定程度的遺傳分化。此外，Wright[29]認為當種群間的居群每代遷移數（Nm）大于1時，表明群體間的基因流動發生多，若小于1，則基因流動少，遺傳分化明顯；4個不同地理種群的Nm為3.2740＞1，說明4個地理種群存在不同程度的基因交流，可能是由于種植地點都在山區且氣候相似，遺傳分化的程度不明顯。

由NTSYSpc-2.10軟件對34個供試菌株構建UPGAM的聚類分析樹可知，當相似性系數設定為0.79時，可將供試菌株分為4個類群：其中第1、3、4類群的菌株分別來至不同煙區，進一步說明各煙區種群間存在一定程度的遺傳分化；但第2種群的31個菌株來至4個煙葉產區，菌株并沒按照地理來源聚為不同的組，表明菌株之間的遺傳多態性與其地理來源無明顯的相關性。

4 結論

本研究采用ISSR分子標記方法對湖北4個煙區的34株長柄鏈格孢菌進行了遺傳多樣性分析，結果表明種群間有豐富的遺傳多樣性，但菌株之間的遺傳多樣性與其地理來源無顯著相關性；4個煙區的菌株產生了一定程度的遺傳分化，但分化的程度不高。該研究通過分析湖北省煙草長柄鏈格孢菌的遺傳多樣性、遺傳結構分化以及遺傳多樣性與地理來源的關系，為后期煙草赤星病的防治奠定了理論基礎。