甘肅玉米莖腐病病菌禾谷鐮孢復合種遺傳多樣性

郭 成，王寶寶，鄧 東，張小杰，王春明，周天旺，李敏權，郭致杰

(1.甘肅省農業科學院 植物保護研究所，蘭州 730070；2. 甘肅農業大學 草業學院，蘭州 730070；3. 石家莊市農林科學研究院，石家莊 050041；4. 中國農業科學院 作物科學研究所，北京 100081；5. 甘肅農業大學 植物保護學院，蘭州 730070)

玉米是中國種植面積最大的農作物[1]。在2021年“秸稈畜牧業”以及“糧改飼”的形勢和政策下，糧飼兼用型玉米具較大的發展潛力，能實現糧食和飼料有效補給，是國內糧食安全生產的關鍵環節[2]。一般年份，莖腐病田間發病率在5%～10%，重發年份達到20%～30%，高感品種的發病率可高達40%～100%，發病率每增加1%，單穗籽粒損失率約提高0.485%，千粒質量損失率約提高0.304%，若發病嚴重將引起30%以上的產量損失[3]。受玉米連年種植、秸稈還田以及氣候逐年變暖等因素影響，田間土壤中病原菌基數逐年上升，玉米莖腐病呈逐年加重趨勢[4]。據郭成等[5]連續多年調查發現，在甘肅及周邊地區玉米莖腐病自然發生率偏高，玉米莖腐病對產量的制約已成為西北地區以及其他玉米主產區生產上的重大問題。

玉米莖腐病的致病菌種群結構和數量較為復雜，目前已報道玉米莖腐病的病原菌多達30種，其中禾谷鐮孢(Fusariumgraminearum)和腫囊腐霉(Pythiuminflatum)為國內玉米莖腐病的兩大主要致病菌[5]。禾谷鐮孢復合種種群組成復雜，已探明包含超過10個菌種，其對環境有廣泛的適應性，不僅能夠侵染多種水果、蔬菜、花卉，還引起小麥和大豆、馬鈴薯、玉米等多種作物病害的發生[6-11]。病原菌群體遺傳多樣性越豐富表明病原菌對環境的適應力越強，較高遺傳多樣性的致病菌具有較強的進化和生存優勢[12]。ISSR、VNTR、RAPD、AFLP和SSR等分子標記技術廣泛應用于禾谷鐮孢復合種群體遺傳結構分析[13-18]。禾谷鐮孢復合種侵染作物后，在寄主體內產生多種次生代謝物和毒素，這些毒素貯藏在植物體內潛在危害巨大[19-20]。

研究病原菌群體結構及其遺傳多樣性特征是應對植物病害的發生、流行及其防治的關鍵[21]。解析致病菌遺傳分化特征及其致病力是抗病品種選育及病害防控的研究基礎[22-23]。有效控制玉米莖腐病對保障中國玉米產業可持續發展具有重要意義。目前缺乏對甘肅省玉米莖腐病致病菌遺傳多樣性的系統報道，本研究以甘肅省玉米莖腐病優勢病原菌禾谷鐮孢復合種為對象，研究其在4大生態區的分布特征及遺傳多樣性，該結果將為甘肅省玉米生產提供風險警示，為玉米鐮孢莖腐病防控提供基礎信息。

1 材料與方法

1.1 菌株信息

從甘肅省4大生態區采集莖腐病樣，采用組織分離法，即在樣品病健交界處用剪刀剪取約 (0.2×0.5) cm大小的玉米莖稈及維管束組織，參照文獻[24-25]進行分離、培養和純化，獲得菌株。經致病性測定后，采用形態學觀察結合TEF-1α基因測序并將所得序列在NCBI中心的BLAST模塊進行比對分析[18]，鑒定到具體種，從中選取114株禾谷鐮孢復合種(FGSC)菌種，平均分布在4個生態區。菌株分布地點見表1。

表1 禾谷鐮孢復合種采集信息Table 1 Distribution and number of F.graminearum species complex

1.2 供試引物

遺傳多樣性引物采用任旭等[17]和Qu等[26]使用的8對VNTR(序號g1～8)和10對SSR引物(序號g9～18)，引物由上海生工生物工程有限公司合成(表 2)。

表2 禾谷鐮孢復合種遺傳多樣性引物Table 2 Genetic diversity primer of F.graminearum species complex

1.3 PCR 擴增體系

參照賀小倫等[27]的方法以25 μL反應體系進行PCR擴增，退火溫度按照引物進行相應的 設定。

1.4 凝膠電泳

對PCR擴增后的產物用6%的變性聚丙烯酰胺凝膠，以恒定功率75 W，電泳90 min，將凝膠置于顯影液中進行條帶自顯影，條帶顯影后用蒸餾水沖去殘液，將凝膠平鋪在燈箱上統計條帶，并拍照。

1.5 電泳譜帶的統計及數據分析

統計好電泳DNA圖譜的條帶信息，有條帶，統計為“1”，無條帶，統計為“0”，進而將多態性位點，也以0和1的方式，轉換為二進制數據。利用軟件POPGENE來計算引物的遺傳距離和遺傳相似系數、多態位點數、Nei’s基因多樣性指數(H)、群體間的Shannon’s多樣性指數(I)[28]。利用軟件NTSYSpc 2.11的非加權組平均法進行不同地區的菌群的SAHN聚類，構建聚類樹。

2 結果與分析

2.1 TEF-1α基因序列聚類結果

禾谷鐮孢復合種經形態學觀察結合TEF-1α基因測序，在114株中鑒定出16株禾谷鐮孢和98株布氏鐮孢(F.boothii)，選取部分菌株構建聚類樹(圖1)。

圖中節點處的數字表示每100次聚類分析中能聚到同一個分支的次數

統計發現，在4大地區中均以布氏鐮孢為主。隴南地區僅有布氏鐮孢分布，23株，未分離得到禾谷鐮孢；河西走廊布氏鐮孢有27株，禾谷鐮孢有3株；隴東地區布氏鐮孢有20株，禾谷鐮孢有11株；隴中地區布氏鐮孢有28株，禾谷鐮孢有2株(表3)。

表3 4個生態區禾谷鐮孢復合種分布Table 3 Distribution of F.graminearum complex species in four regions

(續表3 Continued table 3)

2.2 18對引物的擴增結果

對VNTR和SSR引物擴增結果表明，114株禾谷鐮孢復合種菌株多態性位點數為26。擴增片段集中在50～300 bp，其中引物擴增出的條帶數均在1條以上，表明其多態性可觀。

2.3 禾谷鐮孢復合種不同地理種群遺傳變異和遺傳分化

POPGENE分析結果顯示，4個生態區的禾谷鐮孢復合種群體內的遺傳分化系數為 0.924 9，而其種群間的遺傳分化系數是0.075 1，群體內基因流為0.02，種群間的基因流為 6.153 5，說明禾谷鐮孢復合種遺傳變異主要來源于群體內。禾谷鐮孢復合種群體總基因多樣度為0.425 9，各地理群體內基因多樣度為0.393 9，不同地區間基因多樣度為0.032，綜上，來自同一個生態區禾谷鐮孢復合種菌株間攜帶的基因最為接近，遺傳變異主要來源于群體內。

Shannon’s信息指數和Nei’s基因多樣性指數均以隴中為最高，分別為0.608 7和0.422 9，表明該地區菌株的基因豐富度最高；隴東地區的次之，分別為0.601 7和0.416 8；而隴南地區的則分別為0.564 1和0.392 6；河西走廊的最低，分別為0.501 7和0.343 2，說明基因多樣性與地理區域有較大的相關性。4大生態區Shannon’s信息指數均在0.50以上，Nei’s基因多樣性指數均在0.34以上，表明甘肅省禾谷鐮孢復合種菌群的遺傳多樣性十分豐富。如表4所示，4個生態區禾谷鐮孢復合種的有效等位基因數(Ne)在 1.612 6～1.775 7，平均有效等位基因數為 1.714 0；多態位點百分率在88.46%～100%，平均為95.19%；多態位點數在23～26，平均為 24.75。研究結果表明，甘肅隴東地區和隴中地區菌株遺傳多樣性最為豐富。

2.4 禾谷鐮孢復合種種群的遺傳相似度和遺傳距離

如表4所示，隴中地區和隴東地區的遺傳距離最小，為0.033 1，兩者的遺傳相似度最大，為0.967 4，說明兩地禾谷鐮孢復合種菌株基因組成最接近；而隴東和隴南地區的遺傳相似度最小，為0.888 0，其遺傳距離最大，為0.118 8，說明兩地菌株間差異性最大，親緣關系最遠。4個地區間遺傳相似度為0.888 0～0.967 4，Nei’s遺傳距離為0.033 1～0.118 8(表5)，4個生態地區的菌株均具有豐富的遺傳多樣性。

表4 4個不同生態區禾谷鐮孢復合種種群的遺傳變異參數Table 4 Genetic diversity parameters of F.graminearum species complex among four different geographical populations

表5 甘肅不同生態區禾谷鐮孢復合種種群的遺傳相似度和遺傳距離Table 5 Genetic similarity coefficient and genetic distance among different geographical groups of F.graminearum species complex among four different geographical populations in Gansu

2.5 聚類分析

依據遺傳相似度，用NTSYSpc version 2.11軟件對4個地區種群進行聚類分析(圖2)，當遺傳相似系數為0.938 0時，4個生態群體被聚為3個分支。隴南地區與河西走廊的菌株各自聚為一大類群，分別為第Ⅰ類群和第Ⅱ類群；而隴東地區和隴中地區的菌株共同聚為一大類，為第Ⅲ類群。結果表明，隴東地區和隴中地區禾谷鐮孢復合種菌群親緣關系最為接近，4個生態區禾谷鐮孢復合種群體按其親緣關系可分為3大類，其遺傳親緣關系與地理來源有一定的相關性。

圖2 禾谷鐮孢復合種4個地理種群的UPGMA聚類圖Fig.2 The UPGMA cluster map of 4 geographical populations of F.graminearum species complex

2.6 禾谷鐮孢復合種種群變異分析

禾谷鐮孢復合種不同地理種群的基因變異方差分析表明，在3大類群間的遺傳變異為5.5%，4個種群間的變異占總變異的3.79%，種群內的變異占總變異的90.71%，同時P<0.001 0(表6)，表明禾谷鐮孢復合種種群的遺傳變異主要發生于種群內部。

表6 禾谷鐮孢復合種地理種群的種群變異分析Table 6 Analysis of molecular variance of F.graminearum species complex for geographical populations

3 結論與討論

為明確甘肅4大生態區玉米鐮孢莖腐病優勢病原菌禾谷鐮孢復合種遺傳多樣性，應用8個VNTR標記和10個SSR標記對自甘肅省4大生態區選取的 114 株禾谷鐮孢復合種菌株進行分析，按遺傳相似度劃分為3大類群，表明禾谷鐮孢復合種的親緣關系與菌株的地理來源有一定相關性。Qu等[26]研究表明，來自中國不同地區的禾谷鐮孢菌群體存在高度的遺傳多樣性，整個群體可分為兩個類群，兩個類群有不同的地理分布特征，每個類群又可進一步劃分為幾個亞類。與賀小倫等[27]對來自河南和河北冬小麥區收集的假禾谷鐮孢菌株(F.pseudograminearum)進行遺傳多樣性研究，發現相似系數為0.966時，來自6個地理群體的166株菌被劃分成 2 個大的類群；馬紅霞等[28]以相似系數 0.56為標準，將來自 11 個省(自治區)禾谷鐮孢復合種的92 個菌株劃分成 9 個類群；任旭等[17]將46 株禾谷鐮孢復合種菌株以遺傳相似系數0.45為標準劃分為3個類群；何婧[29]在相似系數為0.682時，將115株禾谷鐮孢劃分為2個類群等研究結果相吻合，進一步說明禾谷鐮孢復合種群體具有豐富遺傳多樣性且存在較明顯的生態地理類型。

本研究發現，甘肅4個不同生態區禾谷鐮孢復合種種群遺傳變異主要來自種群內部。許多學者采用多種分子技術對禾谷鐮孢的遺傳多樣性進行分析，得出的結果與本研究一致，如陳捷等[30]應用RAPD技術明確禾谷鐮孢菌群間的遺傳變異性，易受地域、環境和寄生部位選擇作用的影響，出現明顯的遺傳分化特征； 紀莉景[31]利用TEF-1α基因測序結合RAPD技術發現中國部分地區玉米穗腐病和小麥赤霉病致病菌包括禾谷鐮孢、亞洲鐮孢(F.asiaticum)和布氏鐮孢3個類群，且不同地理來源的類群內的菌株均存在豐富的遺傳多樣性；何婧[29]使用ISSR標記對玉米莖腐病致病菌禾谷鐮孢分析發現，禾谷鐮孢即使在同一聚類分支或是同一地區的菌株在遺傳多樣性上也存在差異，表明禾谷鐮孢菌株間有明顯的生理分化； Talas等[32]利用高通量測序技術對德國13塊試驗田的213個禾谷鐮孢菌株進行群體基因組學分析，發現田塊內群體具有較高的遺傳多樣性。

本研究首次闡釋甘肅省玉米莖腐病優勢致病菌禾谷鐮孢復合種具有較為豐富的遺傳多樣性，根據遺傳相似度將4個生態區的菌群劃分為3大類群，且其種群遺傳變異主要來自種群內部，但其致病力、產毒類型和種群地理分布特征三者之間是否存在相關性有待進一步研究。