黑龍江省和海南省稻瘟病菌中AVR-Pita家族的分布及變異分析

劉瑞，趙羽涵，顧欣怡，王艷霞，靳學慧，吳偉懷，張亞玲

黑龍江省和海南省稻瘟病菌中家族的分布及變異分析

劉瑞1，趙羽涵1，顧欣怡1，王艷霞1，靳學慧1，吳偉懷2，張亞玲1

1黑龍江八一農(nóng)墾大學農(nóng)學院/黑龍江省植物抗性研究中心，黑龍江大慶 163319；2中國熱帶農(nóng)業(yè)科學院環(huán)境與植物保護研究所，?？?571101

【目的】通過檢測黑龍江省和海南省不同稻瘟病菌（）菌株群體中家族的分布情況和變異特征，了解其變異類型的致病表型，為區(qū)域內(nèi)抗病種質(zhì)的篩選與選育提供參考?！痉椒ā客ㄟ^參考NCBI中公布的家族序列分別對啟動子區(qū)和CDS區(qū)設計特異性引物共6對，對2020年采自黑龍江省和海南省不同地區(qū)的397個稻瘟病菌單孢菌株提取DNA，進行PCR擴增。通過電泳檢測結果，分別選取囊括不同地區(qū)的代表菌株對擴增片段進行測序。測序結果與NCBI中相應的堿基與氨基酸序列進行比較分析，并利用水稻抗性單基因系，對不同變異類型的稻瘟病菌菌株進行無毒功能驗證?！窘Y果】在PCR電泳檢測結果中，黑龍江省稻瘟病菌菌株攜帶所有待測基因，分布廣泛且檢出頻率較高；而海南省稻瘟病菌菌株中僅攜帶和，并以低頻率集中存在。無毒基因組成分析結果顯示，黑龍江稻區(qū)的稻瘟病菌菌株復雜多樣，攜帶的基因型種類較海南菌株豐富。PCR產(chǎn)物測序檢測結果顯示，家族以點突變、插入和缺失為主要變異類型劃分為19種,且不同稻瘟病菌群體來源的菌株，其變異類型具有特異性。檢測出10種變異類型，其中（1—5）為黑龍江省稻瘟病菌菌株獨有的變異類型，（6—10）為海南省稻瘟病菌菌株獨有。對這10種變異類型經(jīng)功能驗證后發(fā)現(xiàn)，無毒功能均已喪失。檢測出8種變異類型，其中（1—4）為黑龍江省稻瘟病菌菌株獨有，（5—8）為海南省稻瘟病菌菌株獨有。經(jīng)功能驗證后發(fā)現(xiàn)，無毒功能均已喪失。在黑龍江省僅檢測出1種變異類型?！窘Y論】黑龍江和海南稻瘟病菌群體中家族均為突變后的等位基因型存在，經(jīng)致病性鑒定后，檢測出的所有突變類型均不能被相對應抗性基因和Pi-ta所識別，表現(xiàn)為感病。因此，抗性基因和Pi-ta在黑龍江省和海南省對稻瘟病的抗病育種與利用過程中可聚合其他抗性基因應用來保證品種抗病性。同時，不同地理來源的稻瘟病菌菌株群體中家族的分布和變異類型具有特異性。

稻瘟病菌；家族；變異；無毒基因；黑龍江??；海南省

0 引言

【研究意義】由稻瘟病菌（無性態(tài)：）侵染引起的稻瘟病是對水稻生產(chǎn)造成巨大威脅的三大病害之一[1-2]。選育和種植抗性品種被證明是控制稻瘟病最經(jīng)濟有效和環(huán)境無害的途徑。水稻與稻瘟病菌之間符合經(jīng)典的“基因對基因”系統(tǒng)[3]，它解釋了病原菌中的無毒基因（）與水稻中特定的抗性基因（）相對應，以及品種中缺失某抗性基因如何使該品種成為感病品種[4]。但由于稻瘟病菌生理小種的復雜多變及無毒基因的頻繁變異，不同稻瘟病菌群體遺傳結構發(fā)生變化，產(chǎn)生新的無毒基因型，最終導致水稻品種抗性“喪失”[5]。我國地貌紛繁，生態(tài)環(huán)境多樣，致使不同稻瘟病菌菌群中無毒基因結構復雜，遺傳變異類型多樣[6]。因此，監(jiān)測不同地域稻瘟病菌群體中無毒基因的分布及變異情況，了解無毒基因變異對病原物無毒功能喪失的機制，對研究不同稻區(qū)稻瘟病菌菌群的遺傳多樣性和區(qū)域內(nèi)指導抗病育種方向，并實現(xiàn)抗病持久化具有重要意義?！厩叭搜芯窟M展】據(jù)不完全統(tǒng)計，目前已鑒定的無毒基因有40多個，已被克隆的無毒基因有12個，分別為、[7-16]。Orbach等[7]通過PFLP標記克隆了無毒基因，定位在稻瘟病菌第3號染色體近端粒處，該基因編碼一個含223個氨基酸的分泌蛋白，含典型鋅金屬蛋白酶結構。由于端粒具有不穩(wěn)定性，常發(fā)生點突變、插入和缺失等變異來逃脫的識別，使水稻抗性喪失。Jia等[17]通過蛋白質(zhì)印跡證明是首個編碼蛋白與相應抗病基因產(chǎn)物直接互作的無毒基因。Khang等[18]研究發(fā)現(xiàn)，是一個基因家族，包括、和，其中為原來的基因。研究還發(fā)現(xiàn)，分別與和具有92%和71%的DNA序列同一性。和具有與相對應的無毒功能，且兩基因均存在于種群和種群中，而僅存在于種群中。一些菌株中的和位于端粒附近，兩側是反向重復DNA元件，兩基因的頻繁缺失和擴增可能是重復DNA元件發(fā)生重組所引起的。在水稻稻瘟病抗性基因克隆方面，迄今為止已報道了40個抗瘟基因并完成了分子克隆[19]，與Pi-ta均定位于水稻12號染色體的著絲粒區(qū)域。目前雖已成功克隆，但與Pi-ta仍然不能在F2代中有效分離，因此猜測兩基因可能緊密連鎖或等位基因[19]。Bryan等[20]發(fā)現(xiàn)，對識別特異性的差異取決于Pi-ta蛋白之間的單個氨基酸差異。研究還表明等位基因也出現(xiàn)在含有Pi-ta的水稻品種中，實際上就是Pi-ta為至少加上一個緊密連鎖的其他抗性基因。甘玉姿等[21]在菲律賓稻瘟病菌中發(fā)現(xiàn)一個新的等位基因，命名為，研究還發(fā)現(xiàn)和不具有無毒功能。余歡等[22]在國內(nèi)草類寄主梨孢菌中發(fā)現(xiàn)，具有豐富的變異性，在基因上游存在大量序列缺失現(xiàn)象，基因下游存在缺失、部分缺失和多位點現(xiàn)象。【本研究切入點】水稻在中國種植范圍廣，南起海南島，北至黑龍江，分別以種植秈稻和粳稻為主。黑龍江省地處東北早熟單季稻稻作區(qū)，屬于寒溫帶與溫帶大陸性季風氣候；海南省地處華南雙季稻稻作區(qū)，屬于熱帶季風海洋性。受種植區(qū)內(nèi)生態(tài)、地勢等因素影響，區(qū)域內(nèi)稻瘟病菌的無毒基因可能會存在差異。筆者實驗室已報道了黑龍江省和海南省基因家族在稻瘟病菌中的分布及變異情況[23]，但目前針對南北不同地域差異的稻瘟病菌群體中家族的分布變異情況的研究鮮有報道。因此，為探究南北不同稻瘟病菌菌群中家族的分布及變異機制，對采自黑龍江省和海南省不同水稻產(chǎn)區(qū)的稻瘟病菌進行分析?！緮M解決的關鍵問題】通過對黑龍江省和海南省稻瘟病菌中家族進行PCR檢測、基因序列測序以及對基因的不同變異類型進行功能驗證，明確這兩個稻區(qū)內(nèi)稻瘟病菌家族的分布差異及變異特性，為區(qū)域內(nèi)抗病種質(zhì)的篩選與選育提供參考。

1 材料與方法

1.1 試驗材料

稻瘟病菌供試群體均為2020年于黑龍江省7市14縣和海南省5市6縣采集的稻瘟病穗頸瘟標樣經(jīng)分離獲得的397個單孢菌株，采集地點及編號如表1。

表1 黑龍江省和海南省稻瘟病菌菌株

a：黑龍江省菌株數(shù)Number of strains in Heilongjiang Province；b：海南省菌株數(shù)Number of strains in Hainan Province

1.2 稻瘟病菌基因組DNA的提取

將供試的稻瘟病菌單孢菌株置于PDA平板培養(yǎng)基上活化處理，取5塊菌碟放入滅菌的150 mL液體酵母培養(yǎng)基內(nèi)，置于立式恒溫振蕩培養(yǎng)箱搖菌一周，用真空泵抽濾菌絲后分裝于離心管中。使用CTAB[24]法提取菌絲DNA。采用核酸蛋白測定儀檢測提取的DNA濃度和純度，以確保其質(zhì)量，并將DNA原液用TE稀釋成50 ng·μL-1備用，置于-20℃保存。

1.3 引物設計

于NCBI（https://www.ncbi.nlm.nih.gov）上查找稻瘟病菌家族的基因序列，利用Primer premier 5設計3對特異性引物，啟動子區(qū)序列引物參考Khang等[18]，編碼序列（CDS）區(qū)序列引物參考孟峰等[25]，編碼序列（CDS）區(qū)序列引物參考Shi等[26]。6對引物均委托生工生物工程（上海）股份有限公司合成（表2）。

表2 用于擴增無毒基因的引物

1.4 PCR擴增及電泳檢測

PCR反應體系（20 μL）：dNTPs 1.6 μL，10×Buffer 2 μL，正、反引物各0.3 μL，模板DNA 1 μL，rTaq酶0.1 μL，ddH2O補足至20 μL。

PCR擴增程序：95℃預變性4 min，94℃變性45 s，50—60℃退火45 s，72℃延伸30—90 s，35個循環(huán)，72℃延伸5 min，4℃保存，待檢測。

電泳檢測：PCR產(chǎn)物在電泳液為0.5×TBE和瓊脂糖配置好的1%瓊脂糖凝膠中電泳20 min，150 V。在電泳成像儀下利用Image Lab系統(tǒng)中的條帶檢測功能讀取并記錄電泳條帶，根據(jù)電泳圖譜，運用二進制方法，將有條帶的記為“1”，沒有帶的記為“0”。將數(shù)據(jù)制成Microsoft Excel文件，并對文件中的數(shù)據(jù)進行分析。根據(jù)檢測結果計算各稻區(qū)稻瘟病菌中無毒基因家族的檢出頻率。

1.5 無毒基因的基因序列分析

從無毒基因引物擴增的PCR產(chǎn)物中挑選不同帶型不同水稻產(chǎn)區(qū)的部分菌株，送至生工生物工程（上海）股份有限公司測序，測序結果采用DNAMAN軟件對堿基及核苷酸的序列進行比對分析，確定變異類型。

1.6 稻瘟病菌菌株致病性測定

水稻品種：含有、Pi-ta抗性基因的水稻單基因系IRBLta-K1和IRBLta2-Pi，為國際水稻研究所選育；感病對照品種：麗江新團黑谷（LTH）。

菌株活化及產(chǎn)孢培養(yǎng)：在PDA平板培養(yǎng)基上將供試的稻瘟病菌單孢菌株活化處理，5 d后接種到米糠培養(yǎng)基上進行產(chǎn)孢培養(yǎng)，待菌絲長滿培養(yǎng)基表面時，用涂布棒和無菌水刮洗掉氣生菌絲，置于恒溫恒濕環(huán)境中，用藍紫光燈照射3 d，進行產(chǎn)孢培養(yǎng)。

水稻種植：挑選籽粒飽滿的水稻種子用1.25‰的多菌靈溶液室溫下浸泡24 h，用清水沖洗干凈后置于濾紙保濕的培養(yǎng)皿內(nèi)，于30℃的恒溫培養(yǎng)箱催芽。將發(fā)芽后的種子播種于固定在塑料水培盒（尺寸：225 mm×155 mm×55 mm）內(nèi)的定植籃（上端內(nèi)徑34 mm，下端內(nèi)徑28 mm，高45 mm）中，每個定植籃內(nèi)播種10粒，自來水培養(yǎng)3 d后，營養(yǎng)液培養(yǎng)，營養(yǎng)液采用經(jīng)典的霍格蘭（Hoagland）配方，每4 d換一次營養(yǎng)液，待水稻3葉1心時進行人工噴霧接種。

接種調(diào)查：待米糠培養(yǎng)基表面產(chǎn)生灰色霉層后，每個培養(yǎng)皿用5 mL無菌水洗下孢子，用雙層紗布過濾，制成孢子懸浮液，在40倍顯微鏡下每視野孢子量大于40個，加入1%的吐溫20，搖勻后均勻地噴灑在秧苗上，將其移入25℃的遮光保濕棚內(nèi)培養(yǎng)24 h后返回室溫自然培養(yǎng)，每組處理重復2次，接種5 d后調(diào)查植株發(fā)病情況，調(diào)查標準參考靳學慧[27]。

2 結果

2.1 AVR-Pita基因家族在黑龍江省和海南省稻瘟病菌中的分布

利用基因家族啟動子區(qū)和CDS區(qū)序列引物對397個稻瘟病菌菌株DNA進行PCR檢測（圖1），電泳結果顯示，基因家族共出現(xiàn)3種帶型不同的條帶：有帶、無帶和低帶，用“+”“-”和“1”表示。在黑龍江省186個菌株中擴增出特異性片段，檢出頻率為86.51%；在海南省90個菌株中均擴增出特異性片段，檢出頻率為49.45%，其中東方感城、澄邁金江檢出頻率為0。在黑龍江省178個菌株中擴增出特異性片段，檢出頻率為82.79%；在海南省89個菌株中均擴增出特異性片段，檢出頻率為48.90%，其中東方感城、澄邁金江檢出頻率為0。僅在黑龍江省186個菌株中擴增出特異性片段，檢出頻率為86.51%；而海南省菌株中均未增出特異性片段，檢出頻率為0。

A：AVR-Pita1啟動子區(qū)的部分擴增結果Partial amplification results of AVR-Pita1 promoter region；B：AVR-Pita1 CDS區(qū)的部分擴增結果Partial amplification of AVR-Pita1 CDS region；C：AVR-Pita2啟動子區(qū)的部分擴增結果Partial amplification results of AVR-Pita2 promoter region；D：AVR-Pita2 CDS區(qū)的部分擴增結果Partial amplification of AVR-Pita2 CDS region；E：AVR-Pita3啟動子區(qū)的部分擴增結果Partial amplification results of AVR-Pita3 promoter region；F：AVR-Pita3 CDS區(qū)的部分擴增結果Partial amplification of AVR-Pita3 CDS region；M：DNA marker DL2000；1：低帶Low band；+：存在Existence；-：缺失Missing

2.2 AVR-Pita基因家族在黑龍江省和海南省稻瘟病菌中的基因型

基因家族廣泛分布于所有被調(diào)查的稻瘟病菌菌株中，單個基因的存在和不同組合的分析如表3所示。397個菌株中共檢測到7種不同基因型，分別為、、+、+、+和++。黑龍江省菌株含有全部組合類型，其中攜帶家族的3個無毒基因的菌株數(shù)量最多，為155株，占72.09%且均以組合形式存在。而海南菌株僅有2種類型，和+，其中不攜帶待測家族基因的菌株數(shù)量最多，為92株，占50.55%；攜帶+基因型的菌株數(shù)量為90株，占49.45%，單個家族基因不單獨存在。由此可見，黑龍江菌株攜帶的基因型種類較海南菌株豐富。

表3 參試菌株攜帶無毒基因型的檢測結果

：不含待測基因的基因型，其中，f表示不存在Genotypes that do not contain the gene to be tested, where, f stands for free

2.3 AVR-Pita基因家族的序列及致病表型分析

2.3.1的序列及致病表型分析 分別于兩省各選取20個囊括所有地區(qū)，并在啟動子區(qū)和CDS區(qū)均擴增出特異性片段菌株的PCR產(chǎn)物送去測序，將測序結果與的標準序列進行比對。結果顯示，用于測序的40個菌株表現(xiàn)出10種序列類型（表4），第1類同余歡等[22]發(fā)現(xiàn)的變異類型一致（），與標準序列比對，在啟動子區(qū)-1 004 bp至CDS區(qū)+246 bp處存在1 250 bp基因序列的缺失；其余9類的差異主要體現(xiàn)在堿基的替換與插入，即-413（A/G）、-376（G/A）、-246（T/C）、-116（C/A）、-52（A/C）、18-20GTT（Insert）、207（C/T）、245（G/A）、247（A/G）、263（A/G）、312（T/G）、407（A/G）、520（A/G）、575（G/A）/（T/A）、580（C/T）和583（C/G），其中575位點處A堿基突變了2種類型，將這9種變異類型分別命名為、、、、、、、和。

的9類堿基序列中，-（2，3）、-（6，7，8）和-（9，10）變異類型在CDS區(qū)變異位點一致，因此從中選取、和的氨基酸序列參與序列比對。將5類堿基序列所翻譯的氨基酸序列與標準氨基酸序列進行比對（圖2），結果顯示，207（C/T）處為同義突變，其余堿基突變均使氨基酸發(fā)生錯義。在6（L/Insert）、82（S/N）、83（N/D）、88（K/R）、104（N/K）、174（I/V）和192（C/Y）處發(fā)生氨基酸錯義翻譯；在6（L/Insert）、82（S/N）、83（N/D）、88（K/R）、104（N/K）、136（E/G）、174（I/V）和192（C/Y）處發(fā)生氨基酸錯義翻譯；在6（L/Insert）、82（S/N）、83（N/D）、88（K/R）、104（N/K）、136（E/G）和174（I/V）處發(fā)生氨基酸錯義翻譯；在類型的基礎上增加一處，即194（R/W）；在6（L/Insert）、192（F/Y）和195（H/D）處發(fā)生氨基酸錯義翻譯。

表4 部分菌株AVR-Pita1堿基序列比對結果

-代表在該點的堿基與標準序列相同- represents that the base at this position is the same as the standard sequence

圖2 AVR-Pita1氨基酸序列比對分析

對無毒基因型的10種突變類型的菌株接種水稻單基因系K1（）進行功能驗證，以麗江新團黑谷（LTH）為對照。結果顯示（圖3），-（1—10）均不能被識別而表現(xiàn)為有毒性（V），無毒功能喪失。

2.3.2的序列及致病表型分析 分別于兩省各選取20個囊括所有地區(qū)，并在啟動子區(qū)和CDS區(qū)均擴增出特異性片段菌株的PCR產(chǎn)物送去測序，將測序結果與的標準序列進行比對。結果顯示，用于測序的40個菌株表現(xiàn)出8種序列類型（圖4），均與標準序列有多處差異，主要體現(xiàn)在單堿基的替換、插入與缺失，其中，-115、506和575位點處單堿基各突變了2種類型，即-115（C/T）/（A/T）、506（G/A）/（T/A）和575（G/A）/（T/A）。第1類、第2類和第4類各與標準序列相比均有89處差異，將這3類分別命名為、和；第3類和第7類各與標準序列相比均有88處差異，將這2類分別命名為和；第5類、第6類和第8類各與標準序列相比分別有86、85和87處差異，將這3類分別命名為、和。

圖3 AVR-Pita1變異菌株的致病型

圖4 部分菌株AVR-Pita2堿基序列比對結果

將的8類堿基序列所翻譯的氨基酸序列與標準氨基酸序列進行比對（圖5），結果顯示，8種變異類型均與標準氨基酸序列有多處差異，均為堿基突變使氨基酸發(fā)生錯義。其中，88、169和192位點處各錯義突變了2種類型，即88（K/E）/（R/E）、169（G/Q）/（V/Q）和192（C/Y）/（F/Y）。（1，2，4，7）均與標準序列比對有33處差異；（3，8）均與標準序列比對有32處差異；和與標準序列比對分別有31和30處差異。所有變異位點處均發(fā)生氨基酸錯義翻譯。

對無毒基因型的8種突變類型的菌株接種水稻單基因系PiNo.4（Pi-ta）進行功能驗證，以麗江新團黑谷（LTH）為對照。結果顯示（圖6），-（1—8）均不能被Pi-ta識別而表現(xiàn)為有毒（V），無毒功能喪失。

圖5 AVR-Pita2氨基酸序列比對分析

圖6 AVR-Pita2變異菌株的致病型

2.3.3序列分析 因海南菌株中未擴增出目的片段，故于黑龍江菌株中選取20個囊括所有地區(qū)，并在啟動子區(qū)和CDS區(qū)均擴增出特異性片段菌株的PCR產(chǎn)物送去測序，將測序結果與的標準序列進行比對。結果顯示，用于測序的20個菌株僅表現(xiàn)出1種序列類型（表5），與標準序列比對有多處差異，主要體現(xiàn)在單堿基的替換，即55（G/A）、103（C/T）、122（T/C）、135（A/G）、148（C/A）、183（C/A）、267（A/C）、307（A/G）、315（T/C）、319—321（CGT/AAG）、333（C/A）、479（G/A）、484（A/G）和508（G/A）。

將堿基序列所翻譯的氨基酸序列與標準氨基酸序列進行比對（圖7），結果顯示，135（A/G）和315（T/C）處為同義突變，其余堿基突變均使氨基酸發(fā)生錯義，即19（V/I）、35（P/S）、41（V/A）、50（Q/K）、61（D/E）、89（L/F）、103（K/E）、107（R/K）、111（D/E）、160（R/Q）、162（K/E）和170（A/T）。

2.3.4基因家族變異類型占比 在測序結果中，基因家族序列類型可劃分為19種變異類型，其在黑龍江省和海南省占比如圖8所示。-（1—5）僅在黑龍江省存在，占比較均勻，分別為30.00%、25.00%、5.00%、25.00%和15.00%；反之-（6—10）僅在海南省存在，其中-占比最大，為50.00%，其余基因型占比分別為15.00%、10.00%、5.00%和20.00%。-（1—4）僅在黑龍江省存在，其中占比最大，為70.00%，其余基因型占比分別為10.00%、5.00%和15.00%；反之-（5—8）僅在海南省存在，其中占比最大，為55.00%，其余基因型占比分別為20.00%、15.00%和10.00%。僅在黑龍江菌株中檢測到，且所有測序菌株均有同一種變異類型。以上結果說明南北不同稻瘟病菌菌株群體中家族具有顯著差異，并均含有其特異性變異類型，即不同地理及品種來源的菌株存在另一地區(qū)不存在的變異類型，并且不同菌群來源的菌株具有豐富的遺傳多樣性。

表5 部分菌株AVR-Pita3堿基序列比對結果

圖7 AVR-Pita3氨基酸序列比對分析

圖8 AVR-Pita基因家族變異類型比例分析

3 討論

3.1 黑龍江省和海南省稻瘟病菌無毒基因型的分布

黑龍江省和海南省分別是中國粳/秈稻的主產(chǎn)稻區(qū)，由于我國生態(tài)環(huán)境的差異，稻瘟病菌極易通過熟悉的抗性條件衍生出不同的致病型，同時，抗病品種的大面積單一化種植也會對稻瘟病菌造成選擇壓力，促使其發(fā)生變異導致品種抗性“喪失”[28]，進而使水稻大幅減產(chǎn)甚至絕收。Valent等[29]研究表明，水稻品種的持久抗性與所處稻區(qū)的稻瘟病菌無毒基因的變化相關。因此，及時檢測不同稻瘟病菌群體中無毒基因的分布和無毒基因型的組成情況，可為不同稻區(qū)合理進行抗病基因應用提供科學依據(jù)。

本試驗中檢測到大部分黑龍江省水稻種植區(qū)內(nèi)，稻瘟病菌菌株中攜帶所有待測基因，且出現(xiàn)頻率較高，表明家族在該地區(qū)廣泛存在；而少數(shù)海南省水稻種植區(qū)內(nèi)，稻瘟病菌菌株中攜帶和，出現(xiàn)頻率較低，且所有海南菌株中均不攜帶，由于本文供試的海南菌株均采自稻瘟病發(fā)生較重且具有代表性的地區(qū)，因此，推測海南菌株中不存在或含有較少的，表明該家族在海南地區(qū)以低頻率集中存在并具有地域差異性。這些結果表明，黑龍江和海南省粳/秈稻不同稻瘟病菌群體中所檢測的家族分布差異較大，與廖靜靜等[30]對云南省粳/秈型稻瘟病菌中主要無毒基因分布的研究結果相吻合，表明粳/秈型稻瘟病菌無毒基因的分布是稻瘟病菌對水稻亞種適應性進化的結果。對每個菌株的無毒基因組成進行分析，結果顯示，不同稻瘟病菌群體中同一無毒基因型的出現(xiàn)頻率存在差異，與房文文等的研究結果一致[31-32]。所有供試菌株可劃分為7種無毒基因型，黑龍江菌株含有全部組合類型；而海南菌株僅含有2種類型，說明黑龍江稻區(qū)的稻瘟病菌菌株中攜帶的基因組合種類較為豐富，而海南稻區(qū)較為單一。稻瘟病菌群體無毒基因分布形成明顯區(qū)域特征的原因可能是由于各地區(qū)均有其獨具的生態(tài)條件，長期以來經(jīng)自然選擇的結果，使不同地區(qū)稻瘟病菌群體產(chǎn)生了獨特的無毒基因結構。

周江鴻等[33]研究表明，我國稻瘟病菌群體對攜帶抗性基因的品種表現(xiàn)為中等至強毒力水平，因此推測家族存在高變異性，易失去無毒功能。與本研究結果一致，黑龍江和海南稻瘟病菌群體中家族均為突變后的等位基因型存在，經(jīng)致病型鑒定后，檢測出的所有突變類型均不能被相對應抗性基因和Pi-ta所識別，表現(xiàn)為感病。高清等對粳型水稻中抗瘟基因分布和抗病性的研究結果表明，抗性基因在黑龍江省水稻品種中總體應用較廣，聚合其他抗性基因的抗病率較高，認為抗性基因在黑龍江省水稻抗病育種中具有較高貢獻[34-35]。孟峰等研究表明，黑龍江省各稻區(qū)稻瘟病菌中變異類型相對穩(wěn)定，認為在黑龍江省水稻抗病育種中，需要聚合其他抗性基因才能提高抗病率[25,36]，與本文研究結果一致，抗性基因和Pi-ta在黑龍江省水稻的抗病育種與利用過程中需要聚合其他抗性基因應用來保證品種抗病性。董麗英等[37]研究表明，云南省秈稻中和Pi-ta的抗性頻率均高于80%，兩抗性基因可作為抗原在該地區(qū)應用；曹雪琦[38]研究表明，福建省秈稻中Pi-ta的抗性頻率低于30%，攜帶該基因的水稻品種在福建省應慎重推廣；周瑚等[39]研究表明，湖南省秈稻中和出現(xiàn)頻率低，且和Pi-ta的抗性頻率均在40%左右，攜帶兩單基因的水稻品種在湖南省應慎重推廣。本研究中，海南省稻瘟病菌菌株中和出現(xiàn)頻率低，因此抗性基因和Pi-ta在海南省對稻瘟病的抗病育種應用中同樣需要慎重考慮，可聚合其他抗性基因應用來保證品種抗病性。

3.2 黑龍江省和海南省稻瘟病菌無毒基因型的進化

對遺傳變異進行評估是了解基因進化和稻瘟病菌與水稻共同進化的主要分子機制。由于家族基因位于不穩(wěn)定的端粒處，并在稻瘟病菌1、3—7號等染色體上均有發(fā)現(xiàn)，導致該基因具有高度的可變性和多態(tài)性，進而導致更高水平的遺傳多樣性。家族的變異可發(fā)生在基因的CDS區(qū)、5′端及端粒區(qū)[40]，任何通過點突變、插入與缺失等細微的變異都可以影響整個分泌蛋白的結構，對無毒功能造成影響，最終達到改變其致病性的目的。因此，了解不同菌株群體中無毒基因的變異機制對掌握其不同遺傳多樣性和對區(qū)域內(nèi)抗性品種的選育是至關重要的。

通過挑選不同群體中家族的PCR產(chǎn)物測序后發(fā)現(xiàn)，家族在不同群體中有著豐富的遺傳多樣性和特異性。Sirisathaworn等[41]報道了基因突變是影響稻瘟病菌群體遺傳多樣性的唯一因素，隨后作者對泰國稻瘟病菌中3個的編碼序列進行分析發(fā)現(xiàn)，3個均呈現(xiàn)較低水平的遺傳變異，即泰國稻瘟病菌具有低水平的遺傳多樣性。Kasetsomboon等[42]報道了編碼區(qū)存在15個單倍型，且具有高度的遺傳多樣性。本研究結果同樣證實了以上觀點，家族具有豐富的多態(tài)性位點和遺傳多樣性，該家族在啟動子區(qū)和CDS區(qū)存在多處點突變、插入和缺失現(xiàn)象，共鑒定出19種變異類型，其中有18種無毒功能喪失，另外1種是的突變類型，由于缺乏水稻抗性基因，因此還需要后續(xù)工作進一步驗證此類型是否會引起無毒功能發(fā)生改變。轉座子在啟動子區(qū)或CDS區(qū)插入可導致新的毒力等位基因。Zhou等[43]在毒性野生型菌株B2中發(fā)現(xiàn)區(qū)域有轉座子Pot3插入，導致菌株恢復毒性，但本研究所有供測菌株均未檢測到轉座子Pot3插入情況，與孟峰等[25]的研究結果一致。Dai等[44]報道了存在片段缺失和移碼突變現(xiàn)象，共鑒定出27種變異類型。Chuma等[45]發(fā)現(xiàn)和在染色體上的位置及序列經(jīng)常改變，但卻能在7號染色體上以穩(wěn)定形式存在，但本研究中，在多個位點發(fā)生單堿基替換，與前人研究不符，可能是由于稻瘟病菌與水稻協(xié)同進化過程中，水稻不斷進化產(chǎn)生新的抗性基因，病原菌的無毒基因為了逃避抗性基因給予的選擇壓力，自身發(fā)生變異來克服小種抗性。

本研究還發(fā)現(xiàn)，不同地域來源的稻瘟病菌菌株群體遺傳結構存在差異，與肖丹鳳[46]對不同稻作區(qū)稻瘟病菌遺傳分析的研究結果一致。Huang等[47]報道了中國不同群體的稻瘟病菌呈現(xiàn)出互相接近的遺傳多樣性，并且在地理位置上沒有差異。但本研究中，黑龍江和海南稻作區(qū)內(nèi)稻瘟病菌菌株群體之間的變異類型具有區(qū)域特征，即兩稻區(qū)均含有其特異性等位基因型，與Zhang等[48]研究結果一致，中國南方和北方水稻產(chǎn)區(qū)支持兩個截然不同的病原體群體，更加充分驗證了不同稻區(qū)的稻瘟病菌群體間具有豐富的遺傳多樣性。這種特異性的存在，再一次證實了不同稻瘟病菌群體間的遺傳特性與其所處的生態(tài)環(huán)境和耕作栽培等具有密切的相關性。

4 結論

家族在南北不同稻瘟病菌群體中均有分布且分布差異較大，并均含有其特異性變異類型。黑龍江菌株含有全部待測基因并均以高頻率廣泛存在；而海南菌株均以低頻率集中并僅存在和，且黑龍江菌株較海南菌株無毒基因型組成種類豐富。同時家族在不同群體中均含有豐富的多態(tài)性和區(qū)域特征，兩群體中家族均為突變后的等位基因型存在，共有19種變異類型，經(jīng)致病性鑒定后，有18種失去無毒功能。檢測出的所有突變類型均不能被相對應抗性基因和Pi-ta所識別，表現(xiàn)為感病?？剐曰蚝?i>Pi-ta在黑龍江省和海南省對稻瘟病的抗病育種與利用過程中可聚合其他抗性基因應用來保證品種抗病性。因此，不同群體中稻瘟病菌復雜多樣，應繼續(xù)監(jiān)測稻瘟病菌的變化情況，持續(xù)探究其在不同群體中的變異規(guī)律，及時掌握區(qū)域內(nèi)田間稻瘟病菌動態(tài)，不僅對研究不同稻瘟病菌群體遺傳結構的多樣性具有重要意義，同時還可為區(qū)域內(nèi)抗病種質(zhì)的篩選與選育提供參考。

[1] HAYASHI K, YASUDA N, FUJITA Y, KOIZUMI S, YOSHIDA H. Identification of the blast resistance genein rice cultivars using functional markers. Theoretical and Applied Genetics, 2010, 121(7): 1357-1367.

[2] 張亞玲, 高清, 趙羽涵, 劉瑞, 付忠舉, 李雪, 孫宇佳, 靳學慧. 黑龍江省水稻種質(zhì)稻瘟病抗性評價及抗瘟基因結構分析. 中國農(nóng)業(yè)科學, 2022, 55(4): 625-640.

ZHANG Y L, GAO Q, ZHAO Y H, LIU R, FU Z J, LI X, SUN Y J, JIN X H. Evaluation of rice blast resistance and genetic structure analysis of rice germplasm in Heilongjiang province. Scientia Agricultura Sinica, 2022, 55(4): 625-640. (in Chinese)

[3] FLOR H H. Current status of the gene-for-gene concept. Annual Review of Phytopathology, 1971, 9: 275-296.

[4] AZIZI P, RAFII M, ABDULLAH S, NEJAT N, MAZIAH M, HANAFI M, LATIF M, SAHEBI M. Toward understanding of rice innate immunity against. Critical Reviews in Biotechnology, 2016, 36(1): 165-174.

[5] 毛洧, 陳學偉, 王靜. 水稻抗稻瘟病機制的研究進展. 中國科學: 生命科學, 2022, 52(10): 1495-1510.

MAO W, CHEN X W, WANG J. Recent progress on rice resistance to blast disease. Scientia Sinica Vitae, 2022, 52(10): 1495-1510. (in Chinese)

[6] 吳偉懷, 王玲, 程貫忠, 朱有勇, 潘慶華. 稻瘟病菌群體的分子遺傳學研究——廣東省與云南省稻瘟病菌群體遺傳及致病型結構的比較分析. 中國農(nóng)業(yè)科學, 2004, 37(5): 675-680.

WU W H, WANG L, CHENG G Z, ZHU Y Y, PAN Q H. Studies on molecular genetics of rice blast fungus population——comparison of genetic and pathotypic structures of two rice blast fungus populations derived from Guangdong and Yunnan provinces of China. Scientia Agricultura Sinica, 2004, 37(5): 675-680. (in Chinese)

[7] ORBACH M J, FARRALL L, SWEIGARD J A, CHUMLEY F G, VALENT B. A telomeric avirulence gene determines efficacy for the rice blast resistance gene. The Plant Cell, 2000, 12(11): 2019-2032.

[8] Kang S, Sweigard J A, Valent B. Thehost specificity gene family in the blast fungus. Molecular Plant-Microbe Interactions, 1995, 8(6): 939-948.

[9] SWEIGAED J A, CARROLL A M, KANG S, FARRALL L, CHUNLEY F G, VALENT B. Identification, cloning, and characterization of, a gene for host species specificity in the rice blast fungus. The Plant Cell, 1995, 7(8): 1221-1233.

[10] SCHNEIDER D R S, SARAIVA A M, AZZONI A R, MIRANDA H R, de Toledo M A S, PELLOSO A C, SOUZA A P. Overexpression and purification of, a mutant of the effector proteinfrom. Protein Expression and Purification, 2010, 74(1): 24-31.

[11] COLLEMARE J, PIANFETTI M, HOULLE A E, MORIN D, CAMBORDE L, GAGEY M J, BARBISAN C, FUDAL I, LEBRUN M H, BOHNERT H U.avirulence genebelongs to an infection-specific gene cluster involved in secondary metabolism. New Phytologist, 2008, 179(1): 196-208.

[12] FARMAN M L, LEONG S A. Chromosome walking to theavirulence gene of: discrepancy between the physical and genetic maps. Genetics, 1998, 150(3): 1049-1058.

[13] LI W, WANG B H, WU J, LU G D, HU Y J, ZHANG X, ZHANG Z G, ZHAO Q, FENG Q, ZHANG H Y, WANG Z Y, WANG G L, HAN B, WANG Z H, ZHOU B. Theavirulence geneencodes a predicted secreted protein that triggers the immunity in rice mediated by the blast resistance gene. Molecular Plant-Microbe Interactions, 2009, 22(4): 411-420.

[14] YOSHIDA K, SAITOH H, FUJISAWA S, KANZAKI H, MATSUMURA H, YOSHIDA K, TOSA Y, CHUMA I, TAKANO Y, WIN J, KAMOUN S, TERAUCHI R. Association genetics reveals three novel avirulence genes from the rice blast fungal pathogen. The Plant Cell, 2009, 21(5): 1573-1591.

[15] WU J, KOU Y J, BAO J D, LI Y, TANG M Z, ZHU X L, PONAYA A, XIAO G, LI J B, LI C Y,. Comparative genomics identifies theavirulence effectorthat triggers-mediated blast resistance in rice. New Phytologist, 2015, 206(4): 1463-1475.

[16] ZHANG S L, WANG L, WU W H, HE L Y, YANG X F, PAN Q H. Function and evolution ofavirulence generesponding to the rice blast resistance gene. Scientific reports, 2015, 5: 11642.

[17] JIA Y, MCADAMS S A, BRYAN G T, HERSHEY H P, VALENT B. Direct interaction of resistance gene and avirulence gene products confers rice blast resistance. The EMBO Journal, 2000, 19(15): 4004-4014.

[18] KHANG C H, PARK S Y, LEE Y H, VALENT B, KANG S. Genome organization and evolution of theavirulence gene family in thespecies complex. Molecular Plant-Microbe Interactions, 2008, 21(5): 658-670.

[19] 劉松青. 水稻稻瘟病抗性基因Pi-ta的候選基因篩選與初步分析[D]. 武漢: 中南民族大學, 2019.

Liu S Q. Screening and preliminary analysis of the candidate genes for rice blast resistance genePi-ta[D]. Wuhan: South-Central University for Nationalities, 2019. (in Chinese)

[20] Bryan G T, Wu K S, Farrall L, Jia Y, Hershey H P, McAdams S A, Faulk K N, Donaldson G K, Tarchini R, Valent B. A single amino acid difference distinguishes resistant and susceptible alleles of the rice blast resistance gene. The Plant Cell, 2000, 12(11): 2033-2045.

[21] 甘玉姿, 肖貴, 鄧啟云, 吳俊, 柏斌, 盧向陽, 周波. 菲律賓稻瘟病菌生理小種中及其同源基因的序列與功能分析. 中國生物防治學報, 2018, 34(3): 488-498.

GAN Y Z, XIAO G, DENG Q Y, WU J, BAI B, LU X Y, ZHOU B. The sequence and function analysis ofgene family of rice blast fungus,in Philippines. Chinese Journal of Biological Control, 2018, 34(3): 488-498. (in Chinese)

[22] 余歡, 姜華, 王艷麗, 孫國昌. 無毒基因在不同寄主梨孢菌中的變異研究. 浙江農(nóng)業(yè)學報, 2015, 27(8): 1414-1421.

YU H, JIANG H, WANG Y L, SUN G C. Variability of avirulence genes inisolates from different hosts. Acta Agriculturae Zhejiangensis, 2015, 27(8): 1414-1421. (in Chinese)

[23] 劉瑞, 趙羽涵, 付忠舉, 顧欣怡, 王艷霞, 靳學慧, 楊瑩, 吳偉懷, 張亞玲. 黑龍江省和海南省基因家族在稻瘟病菌中的分布及變異. 中國農(nóng)業(yè)科學, 2023, 56(2): 264-274.

LIU R, ZHAO Y H, FU Z J, GU X Y, WANG Y X, JIN X H, YANG Y, WU W H, ZHANG Y L. Distribution and variation ofgene family in ricefrom Heilongjiang province and Hainan province. Scientia Agricultura Sinica, 2023, 56(2): 264-274. (in Chinese)

[24] 李煥宇, 付婷婷, 張云, 呂天佑, 李遠, 徐秉良. 5種方法提取真菌基因組DNA作為PCR模板效果的比較. 中國農(nóng)學通報, 2017, 33(16): 28-35.

LI H Y, FU T T, ZHANG Y, Lü T Y, LI Y, XU B L. Effect comparison of five methods to extract fungal genomic DNA as PCR templates. Chinese Agricultural Science Bulletin, 2017, 33(16): 28-35. (in Chinese)

[25] 孟峰, 張亞玲, 靳學慧. 黑龍江省稻瘟病菌無毒基因及其同源基因的檢測與分析. 中國水稻科學, 2020, 34(2): 143-149.

MENG F, ZHANG Y L, JIN X H. Detection and analysis ofavirulent geneand its homologous genes in Heilongjiang province. Chinese Journal of Rice Science, 2020, 34(2): 143-149. (in Chinese)

[26] SHI N N, RUAN H C, LIU X Z, YANG X J, DAI Y L, GAN L, CHEN F R, DU Y X. Virulence structure ofpopulations from Fujian province, China. Canadian Journal of Plant Pathology, 2018, 40(4): 542-550.

[27] 靳學慧. 農(nóng)業(yè)植物病理學. 赤峰: 內(nèi)蒙古科學技術出版社, 1999.

JIN X H. Agricultural Plant Pathology. Chifeng: Inner Mongolia Science and Technology Press, 1999. (in Chinese)

[28] 李思博, 魏松紅, 王海寧, 羅文芳, 張優(yōu), 劉志恒. 2015-2016年遼寧省稻瘟病菌種群動態(tài)分析. 沈陽農(nóng)業(yè)大學學報, 2017, 48(3): 284-289.

LI S B, WEI S H, WANG H N, LUO W F, ZHANG Y, LIU Z H. Dynamics of rice blast fungus population in Liaoning province in 2015-2016. Journal of Shenyang Agricultural University, 2017, 48(3): 284-289. (in Chinese)

[29] Valent B, Khang C H. Recent advances in rice blast effector research. Current Opinion in Plant Biology, 2010, 13(4): 434-441.

[30] 廖靜靜, 謝華, 王殿東, 何霞紅. 云南省元陽縣秈/粳型稻瘟病菌無毒基因、和多樣性分析. 植物保護學報, 2019, 46(5): 1057-1064.

LIAO J J, XIE H, WANG D D, HE X H. Polymorphism analysis of,andin indica/japonica-borneisolates in Yuanyang county, Yunnan. Journal of Plant Protection, 2019, 46(5): 1057-1064. (in Chinese)

[31] 房文文. 稻瘟菌群體無毒基因型的時空動態(tài)研究[D]. 北京: 中國農(nóng)業(yè)大學, 2018.

Fang W W. Temporal and spatial dynamics of avirulence genotypes in the populations of[D]. Beijing: China Agricultural University, 2018. (in Chinese)

[32] 王世維, 鄭文靜, 趙家銘, 魏松紅, 王妍, 趙寶海, 劉志恒. 遼寧省稻瘟病菌無毒基因型鑒定及分析. 中國農(nóng)業(yè)科學, 2014, 47(3): 462-472.

WANG S W, ZHENG W J, ZHAO J M, WEI S H, WANG Y, ZHAO B H, LIU Z H. Identification and analysis ofavirulence genes in Liaoning province. Scientia Agricultura Sinica, 2014, 47(3): 462-472. (in Chinese)

[33] 周江鴻, 王久林, 蔣琬如, 雷財林, 凌忠專. 我國稻瘟病菌毒力基因的組成及其地理分布. 作物學報, 2003, 29(5): 646-651.

ZHOU J H, WANG J L, JIANG W R, LEI C L, LING Z Z. Virulence genes diversity and geographic distribution ofin China. Acta Agronomica Sinica, 2003, 29(5): 646-651. (in Chinese)

[34] 高清, 張亞玲, 周弋力, 于連鵬, 聶強, 靳學慧. 黑龍江省粳稻品種稻瘟病主效抗性基因鑒定與抗性評價. 作物雜志, 2021(4): 59-66.

GAO Q, ZHANG Y L, ZHOU Y L, YU L P, NIE Q, JIN X H. Identification of major resistance genes and resistance evaluation to rice blast in japonica rice varieties in Heilongjiang province. Crops, 2021(4): 59-66. (in Chinese)

[35] 相亞超, 王麗麗, 徐凡, 馬殿榮. 抗稻瘟病基因在黑龍江水稻資源中的分布. 分子植物育種, 2018, 16(23): 7705-7717.

XIANG Y C, WANG L L, XU F, MA D R. Study on the distribution of rice blast resistant genes in rice resources of Heilongjiang province. Molecular Plant Breeding, 2018, 16(23): 7705-7717. (in Chinese)

[36] 唐雪婷. 黑龍江省四個稻瘟病流行生態(tài)區(qū)稻瘟病菌無毒基因動態(tài)分析[D]. 大慶: 黑龍江八一農(nóng)墾大學, 2021.

TANG X T. Dynamic analysis of avirulence genes ofin four epidemic ecological regions of Heilongjiang province[D]. Daqing: Heilongjang Bayi Agricultural University, 2021. (in Chinese)

[37] 董麗英, 王群, 劉樹芳, 鄭鳳萍, 李迅東, 楊勤忠. 云南省稻瘟病菌群體對稻瘟病抗性單基因系的致病性分析. 西南農(nóng)業(yè)學報, 2012, 25(2): 467-473.

DONG L Y, WANG Q, LIU S F, ZHENG F P, LI X D, YANG Q Z. Pathogenicity analysis ofpopulations of Yunnan on monogenic lines for resistance to rice blast. Southwest China Journal of Agricultural Sciences, 2012, 25(2): 467-473. (in Chinese)

[38] 曹雪琦. 福建省稻瘟病菌田間種群無毒基因的變異檢測[D]. 福建: 福建農(nóng)林大學, 2020.

CAO X Q. Detection of avirulence gene mutations in the field population of rice blast flungus in Fujian province[D]. Fujian: Fujian Agriculture and Forestry University, 2020. (in Chinese)

[39] 周瑚, 任佐華, 王恒滬, 張譯允, 鄒秋霞, 劉二明. 湖南桃江病圃稻瘟病菌的無毒基因及水稻抗瘟單基因聯(lián)合抗性分析. 微生物學通報, 2017, 44(10): 2353-2360.

ZHOU H, REN Z H, WANG H H, ZHANG Y Y, ZOU Q X, LIU E M. Analysis of avirulence genes ofand resistance association of monogene against blast from rice blast nursery in Hunan Taojiang. Microbiology China, 2017, 44(10): 2353-2360. (in Chinese)

[40] 姜華, 余歡, 王艷麗, 孫國昌. 稻瘟病菌無毒基因序列變異研究進展. 浙江農(nóng)業(yè)學報, 2015, 27(3): 512-520.

JIANG H, YU H, WANG Y L, SUN G C. Progress on sequence variation of avirulence genes in the rice blast fungus. Acta Agriculturae Zhejiangensis, 2015, 27(3): 512-520. (in Chinese)

[41] Sirisathaworn T, Srirat T, Longya A, Jantasuriyarat C. Evaluation of mating type distribution and genetic diversity of threeavirulence genes,,and, in Thailand rice blast isolates. Agriculture and Natural Resources, 2017, 51(1): 7-14.

[42] Kasetsomboon T, Kate-Ngam S, Sriwongchai T, Zhou B, Jantasuriyarat C. Sequence variation of avirulence gene. Mycological Progress, 2013, 12(4): 617-628.

[43] Zhou E, Jia Y, Singh P, Correll J C, Lee F N. Instability of theavirulence genealters virulence. Fungal Genetics and biology, 2007, 44(10): 1024-1034.

[44] Dai Y, Jia Y, Correll J, Wang X, Wang Y. Diversification and evolution of the avirulence gene. Fungal Genetics and Biology, 2010, 47(12): 973-980.

[45] Chuma I, Isobe C, Hotta Y, Ibaragi K, Futamata N, Kusaba M, Yoshida K, Terauchi R, Fujita Y, Nakayashiki H, VALENT B, TOSA Y. Multiple translocation of theeffector gene among chromosomes of the rice blast fungusand related species. PLoS Pathogens, 2011, 7(7): e1002147.

[46] 肖丹鳳. 三個不同稻作區(qū)稻瘟病菌致病性與品種互作研究[D]. 北京: 中國農(nóng)業(yè)科學院, 2013.

XIAO D F. Studies on the pathogenicity ofin three different rice growing regions and their interactions to rice varieties[D]. Beijing: Chinese Academy of Agricultural Sciences, 2013. (in Chinese)

[47] Huang J, Si W, Deng Q, Li P, Yang S. Rapid evolution of avirulence genes in rice blast fungus. BMC Genetics, 2014, 15: 45.

[48] ZHANG Y L, ZHU Q L, YAO Y X, ZHAO Z H, CORRELL J C, WANG L, PAN Q H. The race structure of the rice blast pathogen across southern and northeastern China. Rice, 2017, 10(1): 46.

Distribution and Variation Analysis offamily infrom Heilongjiang Province and Hainan Province

1College of Agronomy, Heilongjiang Bayi Agricultural University/Heilongjiang Plant Resistance Research Center, Daqing 163319, Heilongjiang;2Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571101

【Objective】By detecting the distribution and variation characteristics of thefamily in differentstrains from Heilongjiang Province and Hainan Province, the pathogenic phenotypes of the variation types were understood, which provided a reference for the screening and breeding of disease-resistant germplasm in the region.【Method】Six pairs of specific primers were designed for the promoter region and CDS region by referring to thefamily sequence published in NCBI. DNA was extracted from 397 single spore strains ofcollected from different regions of Heilongjiang Province and Hainan Province in 2020, and PCR amplification was performed. By electrophoresis test results, representative strains from different regions were selected to sequence the amplified fragments. The sequencing results were compared with the corresponding base and amino acid sequences in NCBI, and the avirulence function of different variation types ofstrains was verified by using rice resistant single gene lines.【Result】In the PCR electrophoresis results, thestrains in Heilongjiang Province carried all the genes to be tested, which were widely distributed and had a high detection frequency. However, onlyandwere carried bystrains in Hainan Province, and concentrated in low frequency. The results of avirulence gene composition analysis showed that the strains ofin Heilongjiang were complex and diverse, and the genotypes were more abundant than those in Hainan. The sequencing results of PCR products showed that thefamily was divided into 19 types with point mutation, insertion and deletion as the main variation types, and the variation types of strains from differentpopulations were specific. Ten variation types were detected in, among which- (1-5) was a unique variation type ofstrains in Heilongjiang Province, and- (6-10) was a unique variation type ofstrains in Hainan Province. After functional verification of these ten variation types, it was found that the avirulence functions were lost. Eight variation types were detected in, among which- (1-4) was a unique variation type ofstrains in Heilongjiang Province, and- (5-8) was a unique variation type ofstrains in Hainan Province. After functional verification, it was found that the avirulence functions were lost. Only one variation type ()was detectedin【Conclusion】Thefamily in Heilongjiang and Hainan populations were mutated alleles. After pathogenicity identification, all mutation types could not be identified by the corresponding resistance genesandPi-ta. Therefore, the resistance genesandPi-tacan be used to polymerize other resistance genes to ensure the disease resistance of varieties in the process of disease resistance breeding and utilization of rice blast in Heilongjiang Province and Hainan Province. At the same time, the distribution and variation types offamily instrains from different geographical sources are specific.

;family; mutation; avirulence gene; Heilongjiang province; Hainan province

10.3864/j.issn.0578-1752.2023.03.006

2022-09-28；

2022-11-01

國家自然科學基金（U20A2025）、黑龍江省農(nóng)墾總局科技攻關項目（HKKYZD190205）、黑龍江八一農(nóng)墾大學科研啟動項目（XDB201802，XDB201605）

劉瑞，E-mail：2271031165@qq.com。通信作者張亞玲，E-mail：byndzyl@163.com

（責任編輯 岳梅）