甘肅省玉米大斑病菌交配型組成及遺傳多樣性分析

高建昊 戴子淙 張小杰 周天旺 王春明 洪流 郭成

摘 要 為明確甘肅省玉米大斑病菌交配型組成及遺傳結構，利用交配型鑒定特異引物對采自該省不同地區的玉米大斑病菌進行交配型測定，并采用簡單序列重復區間擴增多態性（ISSR）分子標記技術對供試菌株進行遺傳多樣性分析。結果表明：甘肅省玉米大斑病菌主要包括3種交配型，即A交配型、a交配型和Aa交配型，其比例約為8∶8∶1。24條ISSR引物共擴增出136條條帶，平均每個引物擴增出5.67條，其中多態性條帶134條，多態性條帶比例為98.53%。利用PopGen 32軟件分析發現甘肅3個地區玉米大斑病菌間的親緣關系更近，遺傳距離均小于0.08；隴南地區與陜西省的菌株有更高的遺傳相似度，為0.933 8；用NTSYS-pc軟件對隴東地區的菌株聚類分析，發現ISSR類群的劃分與菌株交配型間并沒有直接關系。說明玉米大斑病菌遺傳類群和交配型之間無明顯相關性，但與菌株地理來源有一定的相關性。

關鍵詞 玉米大斑病菌；交配型；ISSR；遺傳多樣性

玉米大斑病是世界玉米產區普遍發生的一種葉部病害，在玉米整個生育期均可發病，影響玉米的生產，一般年份發病可造成20%左右的產量損失，嚴重時可達到50%以上［1］。近年來，隨著感病品種‘先玉335及其近緣種的廣泛種植，加之氣候及耕作制度的變化，導致玉米大斑病在不同地區和不同年份頻頻發生與流行，現已成為危害甘肅東南部玉米生產的主要葉部病害之一［2］。

大斑凸臍蠕孢菌Exserohilum turcicum （Pass.） Leonard & Suggs［有性態為大斑剛毛座腔菌Setosphaeria turcia （Luttr.） Leonard & Suggs］是引起玉米大斑病的病原菌［2］。該病菌為異宗配合真菌，其有性生殖受MAT基因控制，包括僅含 ?StMAT1-1基因的A交配型、僅含 ???StMAT1-2基因的a交配型和同時包含兩種基因的Aa交配型［3］。在自然條件下不斷地進行有性生殖，病原菌基因有性重組，使得病菌變異頻率頻繁，生理小種也隨之不停變化，對病害的防治與抗病品種造成嚴重威脅［4］。

分子標記能夠在DNA水平上反應病原菌個體及種群間的差異［5-6］，目前已有多種分子標記技術被廣泛應用。ISSR具有可獲信息量高，可重復性強，檢測多態性高且操作簡單等優點［2］ 。谷守芹等［7］利用ISSR分子標記技術分析了玉米大斑病菌的遺傳多樣性；郭麗媛［2］在谷守芹等［7］的基礎上優化了ISSR-PCR分析的最佳反應體系和反應程序，并討論了玉米大斑病菌的遺傳多樣性與菌株地理來源和交配型間的關系；代玉立等［4］利用ISSR法研究了福建省玉米大斑病菌的遺傳多樣性。目前尚未見到對甘肅省玉米大斑病菌遺傳多樣性研究的報道。

根據郭成等［8］的調查，玉米大斑病在甘肅省隴東地區發病嚴重，部分地塊病情級別高達7～9級，平均病株率為87.9%，嚴重危害了當地玉米安全生產。本研究對甘肅省玉米大斑病菌的交配型和遺傳多樣性進行研究，以期為進一步深入研究甘肅省玉米大斑病的流行規律及綜合防控技術提供理論依據。

1 材料與方法

1.1 供試材料

供試病樣：2021-2022年采集甘肅省及陜西省具有典型玉米大斑病癥狀的病葉標樣，風干后帶回實驗室進行病原菌分離。

供試培養基：馬鈴薯葡萄糖瓊脂（potato dextrose agar，PDA）培養基。馬鈴薯200 g，葡萄糖20 g，瓊脂15 g，蒸餾水1 000 mL。

試劑及儀器：2×Power Taq PCR Master Mix、DL1000 DNA Marker和DL2000 DNA Marker均購自近岸蛋白質科技有限公司；TBE購自北京索萊寶科技有限公司；B518229-0100真菌基因組DNA快速抽提試劑盒購自上海生工生物工程股份有限公司；PCR引物均由上海生工生物工程股份有限公司合成。

1.2 試驗方法

1.2.1 病原菌的分離培養 采用組織分離法，剪取病健交界處的組織塊（4 mm×6 mm），于75%的酒精中消毒40 s，然后用無菌水清洗3次，置于PDA平板中。室溫培養3～4 d后轉至PDA上純化2～3次，保存菌株。經分離培養后共獲得60株供試菌株，其中50株甘肅菌株進行交配型鑒定和遺傳多樣性分析，10株陜西菌株僅作為地理種群的對照進行遺傳多樣性分析（表1）。[FL）]

1.2.2 DNA提取 刮取純化的病菌菌絲裝入滅過菌的1.5 mL離心管中，干燥后加入小鋼珠，于液氮中冷卻，用組織研磨機將菌絲研磨成粉末，利用B518229-0100真菌基因組DNA快速抽提試劑盒提取各菌株DNA，測定OD260/OD280，確定其濃度與純度，-20 ℃保存，備用。

1.2.3 病原菌鑒定 將純化后的單菌落接種在PDA平板上，室溫培養7 d，觀察菌落特征，并于普通光學顯微鏡下觀察孢子形態。同時，以提取的DNA為模板，利用引物ITS1（5′- TCCGTAGGTGAACCTGCGG-3′）和ITS4（5′-TCCTCCGCTTATTGATATGC-3′）［6，9］對ITS區序列進行擴增。采用25? μL 的PCR反應體系，即? 12.5? μL 2×Power Taq PCR Master Mix，1? μL上游引物，1? μL下游引物，2? μL DNA模板，ddH2O補足至25? μL。PCR反應程序參考馬周杰等［6］病原菌鑒定的PCR反應程序。PCR擴增產物經1.0%瓊脂糖凝膠電泳檢測拍照后，交于上海生工生物工程股份有限公司測序。將測序結果在GenBank中進行比對，下載相似性大于99%的序列，以蕓薹鏈格孢（Alternaria brassicae）、細極鏈格孢（Alternaria tenuissima）、茄鏈格孢（Alternaria solani）、菜豆殼球孢（Macrophomina phaseolina）等為外類群，以大斑剛毛座腔菌（Setosphearia turcica）和凸臍蠕孢菌（Exserohilum turcicum）為內類群，基于ITS序列構建系統發育樹。

1.2.4 交配型分子鑒定

以提取的50株甘肅省玉米大斑病菌株DNA為模板，參照張泗舉［3］提出的玉米大斑病菌交配型測定引物（表2）進行PCR鑒定。反應體系同“1.2.3”節。參照張泗舉［3］設定PCR反應程序，PCR擴增產物經1.2%瓊脂糖凝膠電泳檢測，根據擴增條帶有無及大小確定交配型。只擴增出800 bp左右目的條帶的菌株為A交配型菌株，只擴增出250 bp左右目的條帶的菌株為a交配型菌株，同時可擴增出2個目的條帶的菌株為Aa交配型菌株，無目的條帶擴增結果的菌株為中性菌株。

1.2.5 ISSR-PCR擴增 參考谷守芹等［7］、郭麗媛［2］和代玉立等［4］優化的ISSR-PCR反應體系和篩選的引物，對所有供試的60株大斑病菌株進行ISSR擴增，引物由上海生工生物工程股份有限公司合成（表3）。反應體系同“1.2.3”節。PCR反應程序參考郭麗媛［2］ISSR-PCR擴增所用程序（退火溫度見表3）。取6? μL擴增產物用的1.5%瓊脂糖凝膠在220 V恒定電壓下電泳約20 min，用凝膠成像系統檢測并拍照。

1.2.6 數據分析 統計電泳條帶，在遷移率相同的位置，有條帶的記為1，無條帶的記為0，將數據結果錄入Excel表中構建原始數據矩陣。利用PopGen? 32計算遺傳距離、遺傳相似度、Neis基因多樣性指數和Shannon信息指數等遺傳變異參數，并利用UPGMA構建基于Neis遺傳距離的聚類圖譜。

2 結果與分析

2.1 病原菌鑒定

所有供試菌株的菌落背面為欖綠色或黑褐色；在光學顯微鏡下觀察（圖1），菌絲為灰黑色，分生孢子梗呈直立狀或彎曲狀，單生或簇生，棕褐色，頂端為淺棕色，有2～6 個橫隔膜，分生孢子為黃褐色，呈長梭形或紡錘形，中下部較寬，有明顯的臍點，具2～8個橫隔膜，大小為（36.56～? 59.14）? μm×（8.19～14.51） μm。

所有供試菌株均用ITS1/ITS4引物擴增出約600 bp的條帶，PCR擴增產物序列結果與GenBank中的玉米大斑病菌的無性態大斑凸臍蠕孢菌（Exserohilum turcicum） 菌株序列（NCBI accession NO. MG727619.1）和有性態大斑剛毛座腔菌（Setosphearia turcica）菌株序列（NCBI accession No. KJ922728.1）進行比對，同源性均達99%以上。利用軟件MEGA 11.0構建系統發育樹（圖2），所有供試菌株均能與E. turcicum和S. turcica聚為一支。因此，所有分離株均被鑒定為大斑剛毛座腔菌（S. turcica）。

2.2 交配型組成分析

利用鑒定交配型的2對特異性引物對甘肅省的50株供試菌株DNA進行PCR擴增，在供試的50株甘肅菌株中，A交配型（僅包含 ?StMAT1-1）共有24株，占菌株總數的48%，a交配型（僅包含 ?StMAT1-2）共有23株，占菌株總數的46%，Aa交配型的（同時包含 ?StMAT1-1和 ?StMAT1-2）共有3株，占菌株總數的6%，未檢測到中性菌株。從甘肅不同地區交配型的分布情況來看，28株來自隴南地區的供試菌株中有17株為A交配型，10株為a交配型，1株為Aa交配型，A交配型和a交配型的比例接近2∶1；15株來自隴東地區的供試菌株中有4株為A交配型，9株為a交配型，2株為Aa交配型，A交配型和a交配型的比例接近1∶2；7株隴中地區的菌株中有3株為A交配型，4株為a交配型， A交配型和a交配型的比例接近1∶1。這表明隴南地區玉米大斑病菌的優勢交配型為A交配型，隴東地區玉米大斑病菌的優勢交配型為a交配型，隴中地區玉米大斑病菌A交配型和a交配型同為優勢交配型。

2.3 ISSR-PCR擴增

對供試的60株玉米大斑病菌進行ISSR-PCR擴增，擴增條帶的大小為100～5 000 bp，24條ISSR引物共擴增出136條條帶，平均每條引物擴增出5.67條條帶，其中多態性的條帶數為134條，多態性條帶比例為98.53%，表明玉米大斑病菌群體存在豐富的遺傳多樣性。將統計的電泳譜帶建立0、1矩陣，利用POPGEN軟件分析，結果顯示（表4），Neis基因多樣性指數和Shannon信息指數均以隴南地區為最高，分別為? 0.394 7和0.565? 5 ，表明該地區菌株的基因豐富度最高；陜西省次之，分別為0.373? 8和0.544 9；而隴東地區分別為0.371 4和0.538 0；隴中地區的則分別為0.329 5和0.474 3，說明玉米大斑病菌基因多樣性與地理區域有較大的相關性。甘肅省3個地區的Shannon信息指數均在0.47以上，Neis基因多樣性指數均在0.32以上，表明甘肅省玉米大斑病菌的遺傳多樣性十分豐富。如表4所示，甘肅省3個地區中有效等位基因數、Neis基因多樣性指數、Shannon信息指數、多態性位點數和多態位點百分率均以隴南地區最高，表明隴南地區菌株遺傳多樣性最豐富。

如表5所示，隴南地區與隴東地區的遺傳距離最小，為0.036 4；遺傳相似度最高，為0.964 3。相比隴東地區和陜西省，隴中地區與隴南地區的遺傳距離更小，為0.043 9；遺傳相似度更高，為? 0.957 1。甘肅3個地區中，隴南地區與陜西地區的遺傳距離最小，為0.068 4；遺傳相似度最高，為0.933 8 。對popgen聚類圖譜分析發現（圖3），在距離2.85處甘肅省3個地區和陜西省可聚為2個類群，分別是甘肅類群和陜西類群；在距離? 1.82處甘肅類群又聚為2個亞類群，其中隴東地區和隴南地區聚為一支，隴中地區聚為一支。以上結果表明，甘肅各地區玉米大斑病菌的親緣關系更近，隴南、隴中、隴東3個地區相比，隴南地區與陜西地區的菌株遺傳相似度更高。說明玉米大斑病菌的遺傳親緣關系與地理來源有相關性。

在用NTSYS-pc軟件對隴東地區的菌株進行聚類分析時發現，不同交配型的菌株可以聚類在同一類群中，說明菌株的交配型與遺傳分化之間并沒有明顯的相關性。

3 討論與結論

在甘肅，玉米大斑病主要危害隴東和隴南地區的玉米生產，是不可忽視的重要病害。中國玉米大斑病菌存在不同的交配類型（A交配型、a交配型和Aa交配型），且可以進行有性生殖，這增加了玉米大斑病菌遺傳變異的幾率。異宗配合真菌的有性生殖受交配基因MAT的控制，交配類型A和a分別由 ?StMAT1-1和 StMAT1-2基因編碼［10］。在本研究中，檢測出甘肅省玉米大斑病菌存在A、a和Aa 3種交配類型，其中A交配型和a交配型比例接近1∶1，與Bunkoed等［11］報道泰國玉米大斑病菌A交配型與a交配型比例接近1∶1的研究結果比較一致；但與孫淑琴等［12］的中國玉米大斑病菌中Aa交配型的兩性菌株數量最大，王利智等［13］的云南省玉米大斑病菌A交配型比例明顯大于a交配型，以及郭麗媛［2］的中國北方玉米產區玉米大斑病菌a交配型比例遠大于A交配型之結果有差異。以上研究結果表明，不同地區存在不同的優勢交配型，且交配型分布比例也存在差異，有些地區交配型組成較為簡單，有些地區交配型組成較為復雜，而不同的交配型有性雜交可能是生理小種趨于復雜的原因之一［12］。因此，監測不同地區玉米大斑病菌的交配型分布及比例，對大斑病菌遺傳分化與病害流行的預測具有重要意義。此外，Nelson［14］發現玉米大斑病菌還存在雙性菌株和中性菌株，認為玉米大斑病菌的交配類型不是簡單的雙極性，但本研究并未發現中性菌株，今后還需擴大樣本容量。

根據郭成等［8］的調查，近年來玉米大斑病在甘肅省發病日益嚴重，除氣候適宜等原因外，可能還與病原菌豐富的遺傳多樣性有關。因此，本研究以大斑病發病最嚴重的隴東地區為中心，采集了隴東地區、隴中地區、隴南地區及臨近甘肅省的陜西省部分地區的大斑病樣用于試驗。本研究對供試的60株玉米大斑病菌進行聚類分析，發現甘肅各地區菌株間的親緣關系更近，其中隴東地區與隴南地區的菌株有更高的遺傳相似度；隴中地區與隴東地區的遺傳相似度及隴中地區與隴南地區的遺傳相似度均高于隴中地區與陜西省的遺傳相似度，表明玉米大斑病菌的變異與地域之間有一定的關系，地域間的距離越近，菌株間的親緣關系越近。這與Ferguson等［15］、 Borchardt等［16］和郭麗媛［2］研究所得玉米大斑病菌的遺傳變異與地域間有一定的關系的結論一致。但與谷守芹等［7］、張明會等［17］、王利智等［13］和馬周杰等［6］分析認為玉米大斑病菌類群劃分與菌株地理來源間的相關性并不顯著的結果不一致。以上結果表明，不同地區的病原菌受氣候、生態環境、人工作業、推廣品種的抗病基因組成等諸多因素影響，其地理來源對遺傳變異的影響表現出差異。本研究的結果表明，在甘肅省內，病原菌的地理來源是影響遺傳多樣性的一個重要因素。同時在本研究中，不同的交配型可以劃分在同一類群中，說明玉米大斑病菌的ISSR多態性與交配型之間的相關性不顯著。這與郭麗媛［2］交配型與菌株遺傳分化的關系不明顯的結果一致；但與谷守芹等［7］和馬周杰等［6］研究發現玉米大斑病菌菌株遺傳類群和交配型在分子水平上存在一定關系的結果不一致。這可能與各地菌株遺傳背景或海拔、氣候等環境因素的差異有關。因河西走廊干旱少雨，大斑病零星發生，甘肅中部干旱，大部分地區大斑病發生偏輕，故所有供試菌株中，主要缺乏甘肅中部地區的菌株，今后還需增大樣本容量，同時進行生理小種的鑒定，進一步確定玉米大斑病菌遺傳多樣性與地域、交配型和生理小種之間的關系以及同一地區的菌株變異趨勢，以期對大斑病的綜合防控起到支撐作用。

參考文獻 Reference：

［1］ 左建英，李育才.山西省玉米大斑病的發生特點和綜合防治技術［J］.內蒙古農業科技，2011（3）：123，128.

ZUO J Y，LI Y C.Occurrence characteristics and integrated control techniques of Exserohilum turcicum in Shanxi province［J］. Journal of Northern Agriculture，2011（3）：123，128.

［2］ 郭麗媛.中國北方玉米產區玉米大斑病菌生理小種監測及遺傳多樣性分析［D］.河北保定：河北農業大學，2013.

GUO L Y.Physiological races monitoring and genetic diversity of Setosphaeria turcica in maize production areas of Northern China［D］.Baoding Hebei：Hebei Agricultural?? University，2013.

［3］ 張泗舉.玉米大斑病菌交配型基因片段的克隆與分析［D］.河北保定：河北農業大學，2008.

ZHANG S J.Cloning and analysis of mating type gene fragments in Setosphaeria turcica［D］.Baoding Hebei：Hebei Agricultural University，2008.

［4］ 代玉立，甘 林，阮宏椿，等.福建省玉米大斑病菌ISSR-PCR反應體系的優化和引物的篩選［J］.福建農業學報，2019，34（7）：810-817.

DAI Y L，GAN L，RUAN H CH，et al.ISSR-PCR optimization and primer selection for analyzing Exserohilum turcicum isolates collected in Fujian［J］.Fujian Journal of Agricultural Sciences， 2019，34（7）：810-817.

［5］ 周艷春.無芒雀麥種質資源遺傳多樣性分析及其核心種質的構建［D］.長春：東北師范大學，2020.

ZHOU Y CH.Genetic diversity analysis and construction of core germplasm in? Bromus inermis Leyss［D］.Changchun：Northeast Normal University，2020.

［6］ 馬周杰，何世道，龐欣宇，等.東北地區玉米大斑病菌交配型組成及遺傳多樣性分析［J］.沈陽農業大學學報，2021，? 52（4）：402-408.

MA ZH J，HE SH D，PANG X Y，et al.Analysis of mating type composition and genetic diversity of Setosphaeria turcica in Northeast China［J］.Journal of Shenyang Agricultural University， 2021，52（4）：402-408.

［7］ 谷守芹，范永山，李 坡，等.玉米大斑病ISSR反應體系的優化和遺傳多樣性分析［J］.植物保護學報，2008，35（5）：427-432.

GU SH Q，FAN Y SH，LI P，et al.Optimization of ISSR reaction genetic diversity analysis of Exserohilum turcicum［J］.Journal of Plant Protection，2008，35（5）：427-432.

［8］ 郭 成，徐生軍，金社林.2015年甘肅玉米病蟲害發生情況調查報告［J］.甘肅農業科技，2016（4）：1-3.

GUO CH，XU SH J，JIN SH L.Investigation report on the occurrence of corn diseases and pests in Gansu in 2015［J］.Gansu Agricultural Science and Technology， 2016（4）：1-3.

［9］ 唐 琳.凸臍蠕孢菌玉米、高粱?；偷闹虏》肿硬町愌芯浚跠］.沈陽：沈陽農業大學，2015.

TANG L.The pathogenic molecular differences of Setosphearia turcica f.sp.zeae and? Setosphearia turcica f.sp.sorghi［D］.Shenyang：Shenyang Agricultural University，2015.

［10］ARIE T，CHRISTIANSEN S K，YODER O C，et al.Efficient cloning of ascomycete mating type genes by PCR amplification of the conserved MAT HMG box［J］.Fungal Genetics and Biology，1997，21（1）：118-130.

［11］ BUNKOED W，KASAM S，CHAIJUCKAM P，et al.Sexual reproduction of Setosphaeria turcica in natural corn fields in Thailand［J］.Agriculture and Natural Resources， 2014，48（2）：175-182.

［12］ 孫淑琴，溫雷蕾，董金皋.玉米大斑病菌的生理小種及交配型測定［J］.玉米科學，2005，13（4）：112-113，123.

SUN SH Q，WEN L L，DONG J G.Identifiction of physiological races and mating type of Exserohilum turcicum［J］.Journal of Maize Sciences，2005，13（4）：112-113，123.

［13］ 王利智，吳景芝，康志鈺，等.云南省玉米大斑病菌的RAPD分析［J］.云南大學學報（自然科學版），2011，? 33（2）：244-248.

WANG L ZH，WU J ZH，KANG ZH Y，et al.Phylogenetic analysis on Exserohilum turcicum isolates from Yunnan province by RAPD［J］. Journal of Yunnan University（Natural Sciences Edition），2011，33（2）：244-248.

［14］ NELSON R R.A major gene locus for compatibility inTrichometasphaeria turcica［J］.Phytopathology，1959，? 49（3）：159-160.

［15］ FERGUSON L M，CARSON M L.Spatial diversity of Setosphaeria turcica sampled from the eastern United States［J］.Phytopathology， 2004，94（8）：892-900.

［16］ BORCHARDT D S，WELZ H G，GEIGER H H.Molecular marker analysis of European Setosphaeria turcica populations［J］.European Journal of Plant Pathology， 1998，104（6）：611-617.

［17］ 張明會，徐秀德，姜 鈺，等.玉米大斑病菌種群遺傳多樣性研究［J］.玉米科學，2009，17（1）：143-146.

ZHANG M H，XU X D，JIANG Y，et al.Genetic diversity of Exserohilum turcicum on maize［J］.Journal of Maize Sciences，2009，17（1）：143-146.

Analysis of Mating Type Composition and Genetic Diversity of Northern Maize Leaf Blight in Gansu

Abstract To clarify the current situation of mating type composition and genetic structure of Exserohilum turcicum in Gansu province，the specific primers were used to identify the mating type of Exserohilum turcicum collected from different regions of Gansu province，and the tested strains were analyzed for genetic diversity using simple sequence repeat interval amplified polymorphism （ISSR） molecular marker technology. The results showed that there were three mating types of the pathogen in Gansu? province，namely “A” mating type （only containing? ?StMAT1-1 gene） ，“a” mating type （only containing ?StMAT1-2 gene）? and “Aa” mating type （containing both ?StMAT1-1 and ?StMAT1-2 genes），its ratio was about 8∶8∶1. A total of 136 bands were amplified from 24 ISSR primers，with an average of 5.67 bands amplified per primer，including 134 polymorphic bands and a 98.53% proportion of polymorphic bands. PopGen32 software was used to analyze the genetic relationship among the three areas in Gansu，and the genetic distance was less than 0.08; the strains from Longnan had a higher genetic similarity （0.933? 8） with those from Shaanxi; NTSYS-pc software was used to cluster the strains in Longdong area，and it was found that there was no direct relationship between ISSR classification and strain mating type. The results indicated that there was no significant correlation between the genetic group and mating type of the pathogen，while there was a certain correlation with the geographical origin of the strain.

Key words Maize leaf spot pathogen; Mating types; ISSR; Genetic diversity