甘肅、青海地區(qū)小麥條銹菌監(jiān)測及群體遺傳多樣性分析

黃苗苗，陳萬權(quán)，曹世勤，孫振宇，賈秋珍，高利，劉博，劉太國

（1甘肅農(nóng)業(yè)大學(xué)植物保護學(xué)院，蘭州730070；2中國農(nóng)業(yè)科學(xué)院植物保護研究所植物病蟲害生物學(xué)國家重點實驗室，北京 100193；3農(nóng)業(yè)農(nóng)村部農(nóng)產(chǎn)品質(zhì)量安全生物性危害因子（植物源）控制重點實驗室，北京 100193；4農(nóng)業(yè)農(nóng)村部國家植物保護甘谷觀測實驗站，甘肅甘谷741200；5甘肅省農(nóng)業(yè)科學(xué)院植物保護研究所，蘭州730070）

0 引言

【研究意義】小麥條銹病是一種氣傳病害，是危害我國小麥生產(chǎn)的主要流行病害之一，在世界各主要麥區(qū)也時有發(fā)生，其病原菌為條形柄銹菌小麥?；停≒uccinia striiformisf. sp.tritici，Pst）。小麥條銹菌主要在緯度較高或高海拔地區(qū)越夏，具有喜涼怕熱、遠距離氣流傳播、分布范圍廣、流行頻率高、暴發(fā)性強、流行速度快、危害程度重等特點[1]。該病害流行時，可造成10%—70%的產(chǎn)量損失，嚴(yán)重時可至絕收[2-3]，長期威脅著我國小麥主產(chǎn)區(qū)，包括西北、西南、黃淮海和長江中下游地區(qū)。新中國成立以來，小麥條銹病在我國年均發(fā)生400萬公頃左右，特別是1950、1964、1990、2002和2017年的5次大流行，小麥產(chǎn)量損失達1 380萬噸[4]。我國小麥條銹病的病菌來源分為秋季菌源基地和春季菌源基地，秋季菌源基地為甘肅隴南隴東、川西北、青海東部等，向條銹病菌冬繁區(qū)秋苗提供菌源[5]。前人在甘肅小麥條銹菌流行關(guān)系研究中，對與青海菌源的關(guān)系少有涉及，沒有將甘肅和青海的菌源作為研究主體，存在樣本量較小、條銹菌的標(biāo)樣未采集于病害流行初期等缺陷。因此，明確甘肅、青海兩省春季流行期間的菌源傳播以及生殖方式對條銹菌大區(qū)流行預(yù)測及確定越夏初始菌源來源等具有重要意義?！厩叭搜芯窟M展】我國西北、華北、西南主要流行區(qū)的條銹菌群體遺傳多樣性和遺傳結(jié)構(gòu)已進行了廣泛而深入的研究，確定了我國小麥條銹菌具有較高的遺傳多樣性水平，尤其是甘肅隴南地區(qū)，不同流行區(qū)群體間存在著一定的基因交流[6-11]。西北地區(qū)的甘肅和青海是我國小麥條銹菌最重要的越夏區(qū)，條銹菌新小種的首次發(fā)現(xiàn)幾乎均在這一區(qū)域。該區(qū)域秋季條銹菌菌源可向東部廣大冬麥區(qū)傳播，對周邊地區(qū)有重要影響[1]。關(guān)于青海、甘肅、新疆、西藏4?。ㄗ灾螀^(qū)）的小麥條銹菌毒性及遺傳多樣性的研究表明，青海和甘肅的條銹菌群體相似性較高，存在一定的菌源交換關(guān)系[12]。青海與甘肅條銹菌群體的菌源交流主要存在于隴南地區(qū)，基因流較強；與臨夏群體的菌源交流則受到一定程度的限制，基因流相對較弱[10]。因此，學(xué)者們多認為青海春季、夏季菌源全部由甘肅隴南傳播到青海[13]?！颈狙芯壳腥朦c】雖然前人對甘肅、青海地區(qū)小麥條銹病做了相關(guān)研究，但關(guān)于小麥條銹菌在甘肅、青海地區(qū)周年循環(huán)的研究，包括越夏、越冬以及春夏季菌源傳播，特別是春季流行期間甘肅、青海菌源傳播以及兩地可能存在的生殖模式等問題，尚缺乏一些直接的證據(jù)?！緮M解決的關(guān)鍵問題】以采自自然發(fā)病條件下甘肅、青海地區(qū)小麥變異觀察圃材料上的小麥條銹菌群體為對象，利用微衛(wèi)星（simple sequence repeat，SSR）分子標(biāo)記的方法研究甘肅和青海兩省6個市（縣）小麥條銹菌春季流行傳播路線及群體遺傳多樣性和群體的生殖模式，闡明甘肅和青海兩省小麥條銹菌的群體遺傳關(guān)系，明確兩省小麥條銹菌之間的傳播流行關(guān)系及菌源交換規(guī)律，進而為兩省小麥條銹病的預(yù)測預(yù)報、確定越夏初始菌源來源和有效治理提供理論依據(jù)。

1 材料與方法

1.1 試驗點及標(biāo)樣采集

選擇歷年種植小麥、條銹病常發(fā)生的地區(qū)作為調(diào)查和研究區(qū)域。甘肅省共4個試驗點：隴南市文縣、隴東平?jīng)鍪嗅轻紖^(qū)、中部麥區(qū)定西市臨洮縣、臨夏州臨夏縣；青海省2個試驗點：西寧市城北區(qū)、海東市互助縣。

供試6套82份小麥品種組成的抗病性變異觀察圃材料由中國農(nóng)業(yè)科學(xué)院植物保護研究所麥類真菌病害研究組提供，根據(jù)各地小麥播種適期，按照編號順序依次種植在試驗田中。每個小麥材料種植1行，行長2 m，行距0.33 m，每隔一段距離（一般18行）設(shè)置當(dāng)?shù)乜?、感病品種各1行作為對照。在試驗區(qū)周圍播種2—3行銘賢169（高度感染條銹病品種）作為誘發(fā)品種。種植管理按照當(dāng)?shù)卮筇锍Ｒ?guī)栽培管理方法進行，同時注意以下幾點：（1）保持水肥適中，防止土壤干旱和植株倒伏；（2）小麥生育期內(nèi)禁用任何殺菌劑；（3）不可進行人工接種，并盡可能遠離接種試驗地；（4）觀察記載自然發(fā)病結(jié)果（始發(fā)期、盛發(fā)期）及各品種的小麥生長狀況。

標(biāo)樣采集：從6個試驗點82份變異觀察圃上共收集到551份始發(fā)期小麥條銹病標(biāo)樣（表1）。每個試驗點調(diào)查病害的普遍率、嚴(yán)重度和反應(yīng)型[14]并采集條形柄銹菌夏孢子堆之間分界明顯，無其他病原菌污染的葉片作為發(fā)病標(biāo)樣，用單張吸水紙包裹葉片，并單個放入干燥的信封中，對每個信封進行編號，并標(biāo)明地名（市、縣、鄉(xiāng)鎮(zhèn)）、小麥品種及采集時間，帶回實驗室后4℃干燥器中保存?zhèn)溆谩?/p>

表1 甘肅、青海兩省6個觀察圃內(nèi)小麥條銹病自然發(fā)病調(diào)查結(jié)果及采集的551份標(biāo)樣信息Table 1 The collection information of the 551 Pst samples from six regions in Gansu and Qinghai provinces and the incidence results under the natural infection by Pst

續(xù)表1 Continued table 1

續(xù)表1 Continued table 1

1.2 小麥條銹菌基因組 DNA的提取及微衛(wèi)星位點的PCR擴增

對采集到的小麥條銹菌標(biāo)樣，小心將單一病斑剪下裝入2 mL離心管中，按采集地點和品種進行編號，試驗過程要嚴(yán)防交叉污染?；蚪MDNA的提取參照天根植物 DNA提取試劑盒的方法進行，使用Nano Drop（ND-1000，USA）測量DNA的濃度，并將 DNA濃度稀釋至 50 ng·μL-1后置于-20℃冰箱保存?zhèn)溆谩?/p>

用于SSR分子標(biāo)記的15對引物[2,15]（表2）均由上海生工合成。PCR反應(yīng)體系及程序參照文獻[2,15]，對擴增產(chǎn)物進行 2%的瓊脂糖凝膠電泳后篩選符合要求的反應(yīng)產(chǎn)物，使用3500XL Genetic Analyzer（Applied Biosystems，USA）對擴增子進行檢測。

表2 本試驗使用的15對SSR小麥條銹菌引物信息Table 2 Information of the 15 SSR primers for Pst used in this study

3500XL Genetic Analyzer上機前的準(zhǔn)備：按照1∶50的比例對PCR產(chǎn)物進行稀釋，吸取1 μL稀釋液加入9 μL的Hi-Di+Liz混合液（Hi-Di∶Liz = 1 000∶15），混勻后全部加入3500XL Genetic Analyzer配套的96孔板中，貼膠貼，在振蕩儀上充分混勻后，2 000 r/min離心1 min收集待測反應(yīng)液于管底，95℃變性5 min，冰浴5 min，隨后上機。

1.3 數(shù)據(jù)處理

擴增子長度由3500XL Genetic Analyzer（Applied Biosystems，USA）檢測分析后，用GeneMaker（V2.7.0）（Soft Genetics）讀取數(shù)據(jù)。然后將數(shù)據(jù)按照二倍體的格式錄入Excel表，利用GenAlEx 6.5[16]將數(shù)據(jù)轉(zhuǎn)化為可供R軟件POPPR v2.5.0分析的CSV數(shù)據(jù)格式，超過一半位點有缺失數(shù)據(jù)的樣品去除。

利用GenAlEx計算以下遺傳參數(shù)：（1）分子方差分析（analysis of molecular variation，AMOVA）；（2）群體遺傳分化系數(shù)（FST）；（3）基因流（gene flow value，Nm）。

利用 POPPR v2.5.0[17]計算以下參數(shù)：（1）觀察到的等位基因數(shù)量（Na）；（2）均勻度（Evenness）；（3）位點缺失比率（locus missing ratio）；（4）每個采樣試驗點的樣本個數(shù)（N）；（5）特定個體在特定位點上攜帶的等位基因的組合，攜帶相同等位基因組合的個體稱為具有相同的多位點基因型（multi locus genotype，MLG）；（6）基因型多樣性（G）；（7）標(biāo)準(zhǔn)化關(guān)聯(lián)指數(shù)（rbarD）[18]；（8）基因型累積曲線（genotype accumulating curve）：反映位點是否足以進行群體分析，基因型累積曲線達到 100%，代表所選擇的位點數(shù)足夠進行群體分析[19]；（9）采用Bruvo’s distance繪制不同個體間的最小時空網(wǎng)絡(luò)圖（minimum spanning networks，MSN）[20]；（10）主成分判別分析（discriminant analysis of principal components，DAPC）[21]。

(2)根據(jù)基巖面等值線圖及現(xiàn)場地面調(diào)查可知場地溶洞、溶槽、石芽、漏斗等巖溶形態(tài)造成穩(wěn)定持力層判斷較困難；易出現(xiàn)部分基礎(chǔ)位于隱伏陡崖上等不利情況。

2 結(jié)果

2.1 甘肅、青海地區(qū)6個試驗點小麥條銹病監(jiān)測

田間調(diào)查結(jié)果顯示，82份全國變異觀察圃材料在甘肅地區(qū)發(fā)病比青海地區(qū)嚴(yán)重，其中在臨夏地區(qū)發(fā)病材料最多（54份），其次是平?jīng)觯?4份）、臨洮（40份）和文縣（13份）；有11份材料（Triticum speltaalbum、Joss Cambier、Armada、川育 23、Avocet S*/Yr24、Avocet S*/Yr26、Jupateco R Yr18、TcLr38、ISr6 Ra、ISr11 Ra、Trident Sr38）在6個地區(qū)均未發(fā)病，占供試材料的13.41%，其余71份材料至少在1個試驗點發(fā)??；青?；ブ?14份材料發(fā)病，西寧僅 4份材料發(fā)??；青海地區(qū)發(fā)病的18份材料僅1份材料（銘賢169）在甘肅4個地區(qū)都發(fā)病，其他17份材料均未在甘肅文縣發(fā)病，但在平?jīng)龊团R夏地區(qū)部分材料發(fā)?。ū?1），顯示文縣和青海群體小麥條銹菌毒性有一定差異。

2.2 小麥條銹菌遺傳多樣性水平

15對引物組合共擴增出81個位點，每對引物組合產(chǎn)生的多態(tài)性位點為2—12個。均勻度值在0.2719—0.8648，表明等位基因頻率在樣本之間分布不均，其中引物RJ21的均勻度接近0.5，即該引物的等位基因在樣本中均等分布。各個位點的缺失比率在 0—9.62%（表3）。

表3 15對引物的遺傳多樣性信息Table 3 Information of genetic diversity about 15 SSR primers

551份樣本克隆矯正后，共鑒定出505個MLG，甘肅和青海群體總的基因型多樣性（0.917）高，群體的基因型多樣性在 0.841—0.974，其中平?jīng)鋈后w的基因型多樣性最高，互助群體次之，甘肅臨洮群體的最低（表4），可以看出，6個群體的基因型多樣性均比較高。

2.3 甘肅與青海小麥條銹菌群體菌源關(guān)系

在505個MLG中（表4），僅有32個MLG被克隆并進行了2—6次重新采樣，各群體內(nèi)部均檢測到重新采樣的MLG。甘肅、青海群體分別有29和3個MLG被重新采樣，其中MLG-184在臨洮的群體中重新采樣6次，MLG-141在臨夏的群體中重新采樣5次，沒有重新采樣的MLG有473個，表明甘肅和青海小麥條銹菌群體中存在明顯的遺傳差異，MLG多樣性與小麥條銹菌流行期的主要無性繁殖模式形成了明顯的對比，即無性繁殖在群體中所占的比例不高。

青海西寧群體與甘肅平?jīng)鋈后w之間的基因流最大，臨夏次之，表明西寧群體與平?jīng)鋈后w菌源交流最頻繁，與臨夏群體菌源交流較頻繁；其與文縣和臨洮群體的基因流小且群體遺傳分化系數(shù)大（0.234和0.246），說明西寧群體與文縣、臨洮群體之間有較大的遺傳分化，菌源之間交流不頻繁。青?；ブ后w與臨夏群體之間基因流最大，菌源交流最頻繁；與文縣群體之間基因流最小且群體遺傳分化系數(shù)最大，說明互助群體與文縣群體之間遺傳分化大，菌源交流不頻繁?？赏茢?，青海地區(qū)小麥條銹菌的菌源大多數(shù)來源于甘肅平?jīng)龊团R夏地區(qū)的傳播。甘肅的4個群體之間，平?jīng)鋈后w與臨夏群體之間的基因流最大，菌源交流最頻繁（表5）。

表4 重采基因型及克隆校正群體的基因型多樣性統(tǒng)計Table 4 Resampled MLGs and genotypic diversity statistics of the clone-corrected data within all populations examined in this study

表5 甘肅和青海兩省6個小麥條銹菌群體成對的群體遺傳分化系數(shù)（對角線下部）和基因流（對角線上部）Table 5 Pairwise population FST values (lower diagonal) and Nm (upper diagonal) among six populations of Pst collected from Gansu and Qinghai provinces

最小時空網(wǎng)絡(luò)圖（MSN）結(jié)果表明，兩省6個群體相互之間未發(fā)現(xiàn)共享MLG。臨洮群體的MLG與其他5個群體之間的遺傳關(guān)系較小，從而形成了一個獨特的MLG基因簇。甘肅省4個小麥條銹菌群體之間存在不同程度的遺傳變異；西寧群體的 MLG緊密聚集在一起，與甘肅臨夏和平?jīng)鋈后w的部分個體菌源遺傳距離相近，支持菌源從甘肅的臨夏或平?jīng)鰝鞑サ轿鲗?，隨后在西寧群體產(chǎn)生了克隆分支。這說明盡管種群內(nèi)部存在多樣性，但個體親緣關(guān)系仍然密切。文縣群體同樣緊密聚集在一起，并與青海西寧群體之間的遺傳距離較遠，來自甘肅臨夏、平?jīng)龊颓嗪；ブ脑S多MLG非常接近，表明這些群體之間遺傳關(guān)系更密切（圖1）。

同樣，非參數(shù)主成分分析（DAPC）驗證了MSN的結(jié)果，青?；ブ臀鲗幍娜后w與來自于甘肅平?jīng)龊团R夏的群體之間菌源關(guān)系最密切，差異最??；與臨洮群體遺傳距離相對較遠且臨洮群體為相對獨立；文縣群體則是一個完全獨立的群體，與其他5個群體之間的差異最大（圖 2）。DAPC的結(jié)果以及群體間較大的 FST（0.131）（表 6）再次說明小麥條銹菌各群體之間存在遺傳變異。

圖1 根據(jù)Bruvo’s距離進行的所有小麥條銹菌分離株的最小網(wǎng)絡(luò)圖Fig. 1 Minimum spanning network (MSN) of all isolates of Pst based on Bruvo’s distance

圖2 不同小麥條銹菌群體主成分判別分析圖Fig. 2 The discriminant analysis of principal components (DAPC) of different Pst populations

2.4 甘肅、青海小麥條銹菌群體繁殖方式分析

結(jié)果顯示，等位基因與個體以及群體之間的基因座相關(guān)聯(lián)這一假說得到了重要支持，整體上甘肅和青海小麥條銹菌群體表現(xiàn)出無性繁殖特征（表 4）。而來自于甘肅文縣、臨夏的群體以及青海西寧的群體具有與其他地區(qū)群體不同的特征，這3個群體存在不顯著的rbarD值（表4、圖3）表示連鎖平衡，可能是有性生殖群體，尤其文縣群體（rbarD=0.0139，P=0.186）顯示出明顯的有性重組特征，其他3個群體均為無性繁殖模式。

圖3 甘肅、青海3個小麥條銹菌有性群體連鎖不平衡分析Fig. 3 Linkage disequilibrium analysis of three Pst populations from Gansu and Qinghai

3 討論

甘肅、青海地處西北內(nèi)陸，地域遼闊，地形地貌復(fù)雜多樣，區(qū)域間氣候差異大，形成了小麥種植呈垂直分布，冬、春小麥成片交錯種植的農(nóng)田生態(tài)景觀，是小麥條銹病周年流行的適宜場所，使其成為我國小麥條銹菌重要的越夏區(qū)。小麥條銹菌除了在甘肅、青海區(qū)域內(nèi)往返傳播循環(huán)外，還通過孢子的遠程傳播危及到我國東部廣大麥區(qū)[1]。條銹菌的遠距離傳播一直是植物病害流行學(xué)中的熱點研究領(lǐng)域之一，為大區(qū)病害流行規(guī)律的研究提供了極具價值的材料，也極大地豐富了流行學(xué)理論，可從根本上提供地域性病害防治與品種抗病性布局的科學(xué)依據(jù)。因此通過分子標(biāo)記方法明確甘肅、青海小麥條銹菌傳播及遺傳多樣性水平，對該地區(qū)小麥條銹病的準(zhǔn)確預(yù)測和指導(dǎo)防治具有重要作用。

表6 甘肅、青海兩省6個小麥條銹菌群體分子方差分析Table 6 Analysis of molecular variance (AMOVA) for six populations of Pst sampled from Gansu and Qinghai provinces

本研究中，春季流行時小麥條銹菌在甘肅和青海之間存在廣泛的基因交流。根據(jù)田間實地調(diào)查條銹病初始發(fā)病時間，甘肅4個地區(qū)比青海2個地區(qū)發(fā)病早，且發(fā)病程度嚴(yán)重。最小時空網(wǎng)絡(luò)圖和主成分判別分析結(jié)果顯示甘肅省4個群體之間，文縣和臨洮群體是2個相對獨立的群體，臨夏和平?jīng)鋈后w菌源關(guān)系最近，交流最頻繁，且平?jīng)龅貐^(qū)發(fā)病時間早于臨夏，兩群體之間存在很大的基因流，可推斷菌源從平?jīng)鰝鞑サ脚R夏，然后進行繁殖存活。文縣群體和其他群體不同的原因主要是甘肅省春季發(fā)病初始菌源不同，甘肅省春季菌源結(jié)構(gòu)復(fù)雜，有來自于甘肅省本地的越冬菌源，也有其他地方的菌源。2018年4月中旬，筆者在甘肅省春季病害實地調(diào)查中發(fā)現(xiàn)與隴南文縣相鄰的陜西漢中寧強、略陽的小麥生長期為灌漿期，其條銹病已進入暴發(fā)期；而甘肅天水、定西地區(qū)的小麥僅處于苗期，其條銹病呈始發(fā)或未發(fā)狀態(tài)，平均病田率、病點率和病葉率均輕于陜西寧強和略陽。陜西黃陵以南是條銹菌在西北的主要越冬場所，也是條銹病的春季流行區(qū)[1]。另外，前人研究也表明隴南地區(qū)菌源可能來自于四川、云南的越冬菌源或陜西的春季流行菌源，而不是當(dāng)?shù)囟熬碵22]。同時，隴中、隴東地區(qū)小麥種植以甘肅育成小麥品種為主，而隴南以四川及20世紀(jì)70年代引進的抗病品種為主[23-24]，差異較大的寄主造成甘肅條銹菌群體出現(xiàn)一定差異。隴南地區(qū)在小麥條銹病綜合治理中引進了較多抗病品種，使得小麥群體的抗病基因趨于多樣，并對抗病品種進行合理布局和輪換，所以多元化的寄主抗病基因和較低的選擇壓力也是隴南地區(qū)條銹菌群體結(jié)構(gòu)不同的主要成因。

本研究中群體差異主要來自各群體內(nèi)部個體之間，這與之前研究各流行區(qū)群體遺傳分化結(jié)果一致[6-8,10]?；蛄鞣磻?yīng)群體的遺傳分化水平，基因流大表示群體間的遺傳分化小。本研究中，群體間成對的基因流均＞1（表5），表明小麥條銹菌在各群體間有一定的基因交流，與LU等[10]研究不同，未發(fā)現(xiàn)6個群體之間共享任何基因型，僅在各個群體內(nèi)部存在少數(shù)基因型重新采樣，青海群體的基因交流不如甘肅群體頻繁，較高的基因型多樣性與小麥條銹菌流行周期的主要克隆性形成了明顯的對比。

青海互助和西寧群體與甘肅隴東平?jīng)龊碗]中臨夏群體菌源之間的基因流要大于其與隴南文縣群體之間的基因流，本研究表明青海的群體與甘肅臨夏、平?jīng)鋈后w之間的遺傳距離比隴南文縣的近，文縣群體是一個相對獨立的群體。青海地區(qū)發(fā)病的變異觀察圃材料，除了銘賢 169，其他材料均未在隴南文縣地區(qū)發(fā)病，也是群體有差異的證據(jù)之一。而LU等研究認為青海春季菌源主要自隴南地區(qū)傳播而來，兩地群體間基因流較強，菌源交流廣泛；與臨夏群體的菌源交流則受到一定程度的限制，基因流相對較弱[10]，這與本研究結(jié)果不同。造成本研究結(jié)果的原因，一是與各地區(qū)的地形地貌特征和區(qū)域間距離相關(guān)，文縣地區(qū)高山林立，山川縱橫，生態(tài)環(huán)境較為封閉，病菌與外地交流需借助較強的上升氣流和適當(dāng)?shù)娘L(fēng)向才可實現(xiàn)，且氣候溫和，冬季比較溫暖，夏季涼爽，有利于病原菌完成周年循環(huán)；二是環(huán)境的改變，隨著近年來溫度的升高，菌株發(fā)生變異，加之種植制度的變化等原因；三是有性繁殖群體存在的影響?？偟膩碚f，導(dǎo)致隴南文縣地區(qū)小麥條銹菌群體遺傳多樣性與其他地區(qū)不同的主要因素是地理氣候特征，其次，有性生殖的存在也有一定的影響[10,25-26]。

小麥條銹病為氣傳病害，條銹菌夏孢子在風(fēng)力作用下可大范圍傳播。其中青?；ブ臀鲗幦后w之間基因交流不頻繁，原因可能是這些地區(qū)有一部分菌源來源于青海已發(fā)病地區(qū)的菌源傳播。在早春病情調(diào)查中發(fā)現(xiàn)，各地區(qū)小麥品種不盡一致，而青海尖扎、貴德地區(qū)的小麥條銹病發(fā)病時間略早于甘肅臨夏等相鄰地區(qū)，青海尖扎地區(qū)5月21日初始發(fā)病，發(fā)病嚴(yán)重度和普遍率較低（3/10/10），甘肅臨夏地區(qū)6月 16日發(fā)病程度（4/40/60） 比尖扎嚴(yán)重，青海貴德6月14日發(fā)病（4/20/40），青海貴德和尖扎兩地群體菌源來源也不相同。有研究表明，青海一些地區(qū)小麥條銹菌可以少量越冬[27]。所以青海春季流行的初始菌源不一定都來自于甘肅省，兩省菌源有一定的交流。由此推斷，青海小麥條銹菌越冬菌源在春季病害流行具有一定的作用。

本研究中，兩省小麥條銹菌群體遺傳多樣性均較高，這與樣本來源于地理跨度較大的不同縣（市）有一定關(guān)系，是病菌受復(fù)雜多樣生態(tài)環(huán)境自然選擇的結(jié)果[28]。同時，樣本來自不同的變異觀察圃品種，小麥品種擁有多樣化的抗性背景，寄主對病原菌群體結(jié)構(gòu)產(chǎn)生了重新配置的作用[29-30]。另外，甘肅和青海地區(qū)是地形相似的自然區(qū)域，或為高寒冬春麥區(qū)，或為山川交錯、冬春麥垂直分布的麥區(qū)，為條銹菌的生長和越夏提供了適宜的氣候條件和生態(tài)環(huán)境，條銹菌在該地有較長的進化史。因此，這些地區(qū)的條銹菌群體有較高的遺傳多樣性和重組事件。6個群體間沒有共享基因型存在，也說明條銹菌可通過遺傳重組、突變、寄主選擇及對當(dāng)?shù)厣鷳B(tài)環(huán)境的適應(yīng)，在各地區(qū)進行了相對獨立的進化。

根據(jù)哈迪-溫伯格定律，無性繁殖的群體由于基因不是完全獨立的遺傳會出現(xiàn)顯著的連鎖不平衡現(xiàn)象，有性生殖的群體因配子的隨機交配使得基因座間在遺傳過程中沒有相關(guān)性，處于連鎖平衡狀態(tài)。甘肅、青海省6個地區(qū)的群體中，甘肅文縣、臨夏的群體以及青海西寧的群體存在連鎖平衡現(xiàn)象（圖 3），說明這些地區(qū)存在有性重組的可能。ZHAO等[31]在甘肅省從銹菌侵染的小檗上分離純化的2個銹子器可以產(chǎn)生條銹菌夏孢子，為甘肅小麥條銹菌群體存在有性生殖提供了最直接的證據(jù)，前人也報道了甘肅地區(qū)存在有性重組現(xiàn)象[9,25-26]。但是關(guān)于青海地區(qū)存在有性繁殖現(xiàn)象少有報道。青海地區(qū)存在大量的轉(zhuǎn)主寄主小檗[32]，且種類豐富。在5—7月份，青海大多數(shù)小檗可在自然條件下受到銹菌侵染，能產(chǎn)生銹子器，并大量分布在冬麥和春麥田邊。青海東部地區(qū)廣泛分布的鮮黃小檗、直穗小檗、置疑小檗、甘肅小檗、近似小檗、細葉小檗、匙葉小檗、松潘小檗、堆花小檗和秦嶺小檗10 種小檗可以作為小麥條銹菌的轉(zhuǎn)主寄主[31,33-34]，這是青海群體可以發(fā)生有性生殖的必要條件之一。隨著全球氣候的變暖，耕作制度的變革，青海省東部地區(qū)冬小麥種植面積不斷擴大，使得小麥條銹病的發(fā)生流行規(guī)律發(fā)生顯著變化。甘肅和青海小檗種類繁多，為小麥條銹菌有性生殖發(fā)生提供了基礎(chǔ)條件。本研究發(fā)現(xiàn)青海西寧地區(qū)可能存在有性生殖，群體具有很高的遺傳多樣性，趙杰等研究表明青海地區(qū)小檗分布較為廣泛，春季小檗銹病發(fā)生非常嚴(yán)重[35]，但目前尚未確定青海其他地區(qū)是否存在小麥條銹菌的有性生殖，需要結(jié)合更大的樣本量和更多的采樣點進一步深入研究證實。

4 結(jié)論

小麥條銹病春季流行期，甘肅和青海群體均具有較高的基因型多樣性，甘肅地區(qū)與青海東部地區(qū)的傳播路線以甘肅平?jīng)?、臨夏到青海的傳播為主，甘肅文縣到青海的傳播為輔。甘肅文縣、臨夏和青海西寧 3個群體存在有性生殖現(xiàn)象，對該地區(qū)條銹菌豐富的遺傳多樣性的形成具有一定作用。