摘要：大豆是世界四大糧食作物之一，與小麥、水稻和玉米并列。大豆胞囊線蟲（SCN，Heterodera glycines Ichinohe）是全球大豆最有害的病原體之一。利用大豆胞囊線蟲3號生理小種（SCN 3）、感病品種綏農14、抗病品種東農L10為材料，進行轉錄組分析，篩選出差異表達基因GmWRKY-3。為探究大豆GmWRKY-3基因的生物學功能及潛在用途，克隆GmWRKY-3基因，對GmWRKY-3基因的理化性質、編碼的蛋白質的二級結構和三級結構以及系統(tǒng)進化等生物信息學分析；構建pCAMBIA1302-GmWRKY-3利用農桿菌侵染煙草，并進行亞細胞定位分析；構建pCAMBIA3300-GmWRKY-3植物重組質粒，轉化農桿菌K599后瞬時侵染大豆毛狀根，初步探究GmWRKY-3基因的功能。結果表明，GmWRKY-3基因編碼的蛋白分子量約為 24.93 ku，等電點為8.46，分子式為C1 086H1 722N312O343S9，不穩(wěn)定系數(shù)為67.73，屬于不穩(wěn)定蛋白質。二級結構以無規(guī)則卷曲為主，三級結構由1個α-螺旋和5個β-折疊組成；該基因含有多個光響應元件和多種激素順式作用元件等。亞細胞定位結果表明，GmWRKY-3基因主要定位于細胞核。農桿菌瞬時侵染大豆毛狀根結果表明，與pCAMBIA3300空載體轉基因植株相比，pCAMBIA3300-GmWRKY-3重組質粒轉基因受體植株東農50中根部線蟲數(shù)量相對較少。上述結果表明，GmWRKY-3基因對SCN有一定的抵抗能力。

關鍵詞：大豆；大豆胞囊線蟲；GmWRKY-3；亞細胞定位

中圖分類號：S435.651" 文獻標志碼：A

文章編號：1002-1302（2024）23-0030-08

劉" 芳，曲" 碩，孫浩文，等. 大豆胞囊線蟲相關基因GmWRKY-3克隆及功能探究［J］. 江蘇農業(yè)科學，2024，52（23）：30-38.

doi：10.15889/j.issn.1002-1302.2024.23.004

收稿日期：2023-10-18

基金項目：黑龍江省重點基金（編號：ZD2022C002）；黑龍江省重點研發(fā)項目（編號：JD22A015）；國家自然科學基金（編號：U22A20473）；國家現(xiàn)代農業(yè)崗位體系項目（編號：CARS-04-PS07）；東北農業(yè)大學科研項目（編號：NEAU2023QNLJ-003）。

作者簡介：劉" 芳（2000—），女，黑龍江齊齊哈爾人，碩士研究生，主要從事大豆作物遺傳育種，E-mail：18846913802@139.com；共同第一作者：曲" 碩（1994—），男，吉林長春人，博士研究生，主要從事大豆作物遺傳育種，E-mail：531140417@qq.com。

通信作者：韓英鵬，博士，教授，博士生導師，主要從事大豆基因組學遺傳育種研究，E-mail：hyp234286@aliyun.com；趙" 雪，博士，研究員，博士生導師，主要從事大都基因組學和分子育種研究，E-mail：xuezhao@neau.edu.cn。

大豆［Glycine max （Linn.） Merr.］原產于我國，為豆科大豆屬一年生草本植物，是我國重要的糧食和油料作物［1］。生物脅迫是造成大豆減產的主要因素之一，其中大豆胞囊線蟲（SCN）病是大豆主要的病蟲害之一［2］。大豆胞囊線蟲在世界各地的大豆主產區(qū)均有廣泛分布，在我國首次發(fā)現(xiàn)于東北地區(qū)［3］。SCN平均每年造成的大豆產量損失為 10%～30%［4-6］。目前主要的防治措施為輪作和施用農藥，但農藥嚴重影響生態(tài)環(huán)境，因此提高大豆對SCN的抵抗能力迫在眉睫。

大豆胞囊線蟲的生命周期為3～4周，分為3個主要階段，即卵期、幼蟲期和成蟲期，其中幼蟲期又可以分為4個階段（J1～J4）［7］。卵在胞囊內形成，卵蛻皮進入J1階段，J1階段蛻皮進入蠕蟲狀的J2階段，J2階段是線蟲唯一具有侵害大豆根系能力的階段［8］。SCN通過表皮細胞侵染大豆根部，使其不能正常生長，進而造成大豆減產。合胞體是SCN和大豆宿主之間的生物營養(yǎng)界面，它的功能必須保持約3周的時間，線蟲才能通過J3、J4階段的生長發(fā)育，形成成年雌性或雄性。雌性SCN的體內攜帶幾百個卵，在根表皮生長，同時在根維管圓柱體上保持不斷擴大的合胞體攝食位點。雄性SCN放棄合胞體，離開根部，遷移給雌性受精。受精的雌性SCN通常會釋放一些卵，這些卵會釋放后不久孵化，從而在同一生長季節(jié)實現(xiàn)線蟲種群增加和植物損傷。但仍有幾百個卵通常保留在雌性體內，隨著雌性線蟲死亡，形成硬化的胞囊。目前大豆對SCN的主要抗性位點為Rhg1和Rhg4，伴隨著不同生理小種的突變和環(huán)境的變化，發(fā)現(xiàn)更多的抗性位點具有重要意義。

WRKY轉錄因子是植物中最大的轉錄調節(jié)因子家族之一，是調節(jié)許多植物過程信號網(wǎng)的組成部分［9］。WRKY在植物中的表達量非常豐富，與植物的生長發(fā)育密切相關，其中最主要的作用是調節(jié)植物的免疫與病原體發(fā)育。目前，距第1次發(fā)現(xiàn)該轉錄因子已超過20年。近些年關于WRKY的研究越來越廣泛，如在耐旱、耐鹽、抗病等方面都有所研究，且在煙草、馬鈴薯、擬南芥、水稻等植物中均已被鑒定和分析［10-13］。WRKY第1次是在甘薯中被鑒定出來的［14］。WRKY功能主要是由其結構域決定的，長度約為60個氨基酸殘基，WRKYGQK位于其N末端，其C端還具有鋅指結構。

WRKY基因在很多植物中都發(fā)揮著重要功能，但相關研究主要集中在擬南芥和一些其他作物中，而在大豆胞囊線蟲上的研究較為少見。楊雪梅等發(fā)現(xiàn)，擬南芥WRKY51轉錄因子可以抑制RPW8.1啟動子的活性，通過敲除該基因可以調控RPW8.1啟動子的活性，從而影響植物對白粉病的抗性，說明WRKY51介導基本防御反應［15］。Meher等在對稻瘟病抗性基因進行研究時發(fā)現(xiàn)，OsWRKY蛋白可與2個稻瘟病抗病基因Pi2和Pi54的啟動子區(qū)結合，對稻瘟病的發(fā)展具有重要意義［16］。劉麗麗等采用VIGS策略，沉默基因SlWRKY80，發(fā)現(xiàn)植株根部對線蟲的敏感性增強，根結數(shù)量增多，植株病情指數(shù)由感病態(tài)變?yōu)楦吒袘B(tài)，說明SlWRKY80可能作為正調控因子參與番茄對根結線蟲的抵抗［17］。Wang等在對柑橘的研究中發(fā)現(xiàn)，CsWRKY33定位于細胞核，瞬時過表達CsWRKY33可顯著增強寄主對綠霉病的抗性［18］。通過過往的研究發(fā)現(xiàn)，WRKY轉錄因子對植物抵抗逆境脅迫具有重要意義。

目前，關于WRKY基因家族在SCN脅迫下發(fā)揮的功能和機制還沒有相關研究報道。因此對于抗SCN的研究提供了重要的基因資源。本研究利用大豆胞囊線蟲3號生理小種（SCN3）、感病品種綏農14、抗病品種東農L10為材料，進行轉錄組分析，篩選出差異表達基因GmWRKY-3，通過克隆GmWRKY-3基因，對GmWRKY-3基因編碼的蛋白質進行理化性質和蛋白質二級、三級結構分析以及系統(tǒng)進化樹生物信息學分析。構建pCAMBIA1302-GmWRKY-3重組質粒進行亞細胞定位分析，并構建pCAMBIA3300-GmWRKY-3重組質粒，轉化農桿菌K599，對大豆毛狀根進行瞬時侵染，初步探究GmWRKY-3基因的功能。

1" 材料與方法

1.1" 試驗材料

本研究于2023年在東北農業(yè)大學大豆遺傳育種實驗室進行。大豆感病品種綏農14、東農50和抗病品種東農L10均來自東北農業(yè)大學向陽基地（45.76°N，126.93°E）。大豆胞囊線蟲3號生理小種取自黑龍江省大慶市安達試驗基地。具體鑒定方法參照文獻［19］。

1.2" 試驗方法

1.2.1" GmWRKY-3基因克隆

利用TRIzol法提取抗病大豆品種東農L10 的mRNA，利用TOYOBO ReverTra Ace qPCR RT Master Mix with gDNA Remover試劑盒反轉錄合成 cDNA，置于-20 ℃冰箱保存?zhèn)溆谩?/p>

1.2.2" GmWRKY-3基因擴增

在Phytozome數(shù)據(jù)庫（https：//phytozome-next.jgi.doe.gov/）中查找GmWRKY-3基因全長序列，根據(jù)CDS序列設計引物（表1）。反應體系為20 μL，其中上下游引物（10 μmol/L）各1 μL、模板2 μL、KOD ONE 10 μL、ddH2O 6 μL。反應程序：98 ℃預變性3 min；98 ℃變性10 s，58 ℃退火5 s，68 ℃延伸5 s，35個循環(huán)；68 ℃終延伸10 min。擴增程序結束后，取4 μL PCR產物用2%瓊脂糖凝膠進行電泳檢測。

1.2.3" GmWRKY-3基因亞細胞定位分析

根據(jù)GmWRKY-3基因的CDS序列，選取pCAMBIA1302載體的NcoⅠ酶切位點側翼各15 bp同源序列進行引物修飾，利用Primer 5設計基因亞細胞定位引物（表1）。以抗病品種東農L10 的cDNA為模板，進行PCR擴增，反應體系同“1.2.2”節(jié)。參照南京諾唯贊生物科技股份有限公司的ClonExpress Ⅱ One Step Cloning Kit試劑盒的操作步驟，應用同源重組技術，將目的片段與酶切后的pCAMBIA1302載體進行in-fusion連接，取GmWRKY-3-sub PCR擴增產物1 μL、pCAMBIA1302載體 3 μL、5×Mix Buffer 4 μL、酶2 μL、ddH2O 10 μL 混勻，30 ℃連接30 min。將連接好的質粒轉化至大腸桿菌DH5α中，經菌落PCR擴增及測序，確定已成功將GmWRKY-3連接入pCAMBIA1302載體。參照說明書，將pCAMBIA1302-GmWRKY-3轉化至農桿菌GV3101中，涂布于含有Kana抗性的LB固體培養(yǎng)基上，挑取的單菌落轉移至含有液體LB的2 mL EP管中，之后轉移至50 mL液體LB培養(yǎng)基中進行擴大培養(yǎng)；室溫收集菌體；使用MES重懸至D600 nm=1.0，室溫黑暗靜置2～3 h；注射至 4 周齡煙草葉片背部，25 ℃ 黑暗培養(yǎng)12 h后，正常培養(yǎng) 2 d。采用注射器在侵染葉片上選取4～6個葉圓片進行制片；隨后在倒置熒光顯微鏡下檢測綠色熒光信號。

1.3" 基因生物信息學分析

1.3.1" 理化性質及編碼蛋白質分析

利用Expasy（https：//web.expasy.org/protparam/）分析GmWRKY-3轉錄因子的理化性質；使用SOPMA（https：//npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=npsa_sopma.html）分析GmWRKY-3轉錄因子編碼蛋白的二級結構，用SWISS-MODEL（https：//www.swissmodel.expasy.org）分析蛋白的三級結構。研究差異表達基因編碼蛋白的二級、三級空間結構。

1.3.2" 啟動子區(qū)域順式作用元件分析

利用PlantCARE分析GmWRKY-3轉錄因子啟動子區(qū)域順式作用元件。

1.3.3" 系統(tǒng)發(fā)育樹分析

利用MEGA 10.0軟件對基因編碼蛋白進行比對，采用鄰接法（NJ法）構建系統(tǒng)發(fā)育樹，分析基因在不同物種中的分布情況和親緣遠近關系。

1.4" pCAMBIA3300-GmWRKY-3重組質粒連接及測序

根據(jù)GmWRKY-3基因的CDS序列、pCAMBIA3300載體的XbaⅠ酶切位點側翼各15 bp同源序列進行引物修飾，利用Primer 5對GmWRKY-3基因進行引物設計，結果見表1。PCR擴增、目的片段與載體的連接、重組質粒的轉化與鑒定過程參照“1.2.2”“1.2.3”節(jié)。

1.5" pCAMBIA3300-GmWRKY-3農桿菌瞬時侵染

提取測序正確的大腸桿菌pCAMBIA3300-GmWRKY-3的質粒，并轉化至農桿菌K599中。初步鑒定過程參照“1.2.3”節(jié)。選取籽粒大小均勻、完整、無病斑的感病品種綏農14與東農50，播種至蛭石中。蘸取上述發(fā)根農桿菌K599菌液涂布于LB固體培養(yǎng)基上，置于28" ℃培養(yǎng)箱中過夜培養(yǎng)，將菌落輕輕刮下，使用共培養(yǎng)液體培養(yǎng)基重懸菌體，直至D600 nm=0.8，得到菌液侵染液。選取3～4 d苗齡的幼苗，使用1 mL注射器，在幼苗子葉節(jié)點或近端下胚軸處注射菌液，重復2～3次，進行菌液接種，并在傷口下端約1 cm處剪斷，放回蛭石中培養(yǎng)2～3周，取適量根部組織利用Bar試紙條進行轉基因根系鑒定。將在Bar試紙條水平上檢測為轉基因陽性的植株，保留陽性根部組織，移栽至接種大豆胞囊線蟲卵懸液的病土中，繼續(xù)培養(yǎng)2～3周。

1.6" 陽性根表型鑒定

以感病品種東農50為受體試驗材料，將含有pCAMBIA3300-GmWRKY-3重組載體與pCAMBIA3300空載體的K599農桿菌菌液接種于根系生長點，接種15 d后，進行側根取樣，采用Bar試紙條進行鑒定，根據(jù)Bar試紙條驗證結果，選取轉基因陽性根接種HG type 0（3號生理小種）。接種SCN 15 d后，采用酸性品紅染色法分別對2個試驗組合（接種pCAMBIA3300-GmWRKY-3重組載體和pCAMBIA3300空載體農桿菌菌液）發(fā)生突變和未發(fā)生突變的大豆根部進行染色，在體式顯微鏡下觀察線蟲侵染情況，每株取3個側根作為生物學重復。

2" 結果與分析

2.1" GmWRKY-3基因克隆

以抗病品種東農L10大豆根部的cDNA為模板進行PCR特異性擴增后，將擴增產物利用2%瓊脂糖凝膠進行電泳檢測，克隆載體測序結果（圖1）表明，目的片段大小為663 bp，該序列長度與GmWRKY-3 CDS序列相符。

2.2" 基因序列生物信息學分析

2.2.1" GmWRKY-3基因理化性質及編碼蛋白質分析

使用Phytozome在線網(wǎng)站查找GmWRKY-3的基因信息，該基因的CDS序列長度為663 bp，編碼220個氨基酸。利用ExPASy在線工具對GmWRKY-3編碼蛋白的基本理化特性進行分析，結果顯示，該

蛋白分子量約為24.93 ku，等電點為8.46，分子式為C1 086H1 722N312O343S9，脂肪指數(shù)為51.37，不穩(wěn)定系數(shù)為67.73，表明該蛋白為不穩(wěn)定蛋白質。該蛋白質中絲氨酸個數(shù)最多，占總氨基酸數(shù)量的10.0%，其次為谷氨酸、蘇氨酸，占比均在8.0%以上，最少的是蛋氨酸，占比為1.4%。通過SOPMA網(wǎng)站對GmWRKY-3編碼的氨基酸序列進行二級結構預測，結果（圖2）顯示，不同結構類型的占比分別為無規(guī)則卷曲48.64%、α-螺旋32.27%、延伸鏈15.91%、β-轉角3.18%。使用SWISS-MODEL對三級結構進行同源建模，預測該蛋白的結構域主要由1個α-螺旋和5個β-折疊組成。

2.2.2" 啟動子區(qū)域順式作用元件分析

使用PlantCARE啟動子預測網(wǎng)站分析GmWRKY-3基因的順式作用元件，結果（表2）表明，上游區(qū)域含有MYB轉錄因子結合位點，G-box、TCT-motif光響應元件，MBS參與干旱誘導的MYB結合位點，ABRE、Myb等脅迫相關激素的響應元件。MYB轉錄起始關鍵調控位點是最多作用元件。根據(jù)以上結果推測，GmWRKY-3基因的表達可能受到光、脅迫、激素等多種信號的調控，從而影響植物抵抗生物脅迫的能力。

2.2.3" 系統(tǒng)進化樹分析

為了解GmWRKY-3的進化關系，從NCBI（https：//www.ncbi.nlm.nih.gov/）數(shù)據(jù)庫中篩選出與GmWRKY-3基因編碼蛋白同源性較高的多種植物的WRKY蛋白序列，包括向日葵（sunflower）、傘形葡萄（umbrella grape）、角荬菜（endive）、菜豆（pea bean）、三葉草（whitetip clover）、日本蓮花（Japanese lotus）、白蕨（white fern）、龍眼雙果（longan double fruit）、本氏煙草（ben tobacco）等，同時利用MEGA 10.0軟件構建與大豆同源度較高的物種間的系統(tǒng)進化樹， 結果（圖3）表明，大豆的GmWRKY-3基因與同物種大豆WRKY40轉錄因子（XP 025983959.1）親緣關系最近，不同物種之間與野生大豆（RZC15262.1）、赤豆（XP 017429633.1 ）、洋紫荊（KAI4337249.1）、木豆（XP 020206663.1、KYP34952.1）等聚為同一分支，親緣關系較近，同源性較高。

2.3" pCAMBIA1302-GmWRKY-3載體構建及亞細胞定位

將重組質粒pCAMBIA1302-GmWRKY-3連接轉化至DH5α后，在具有Kana抗性的LB培養(yǎng)基上進行篩選培養(yǎng)，使用引物GmWRKY-3-1302-F和 GmWRKY-3-1302-R（表1）對陽性單菌落進行特異性擴增，對獲得含有目的條帶的陽性菌液（圖4）進行測序，并保存至-80 ℃冰箱中。

將含有重組質粒pCAMBIA1302-GmWRKY-3的農桿菌K599菌液注射入本氏煙草中。結果顯示，在瞬時表達pCAMBIA1302-GmWRKY-3融合蛋白的煙草葉片細胞中，細胞核和細胞膜上均能檢測到熒光信號，但主要定位于細胞核；而單獨表達

pCAMBIA1302的對照細胞，熒光信號分布于整個細胞（圖5）?？梢奊mWRKY-3蛋白屬于核定位蛋白。

2.4" pCAMBIA3300-GmWRKY-3載體構建

將重組質粒pCAMBIA3300-GmWRKY-3連接轉化至DH5α后，在具有Kana抗性的培養(yǎng)基上進行篩選培養(yǎng)，使用引物GmWRKY-3-3300-F和 GmWRKY-3-3300-R（表1）對陽性單菌落進行特異性擴增，對獲得的目的條帶正確的菌液（圖6）進行測序，并保存于-80 ℃冰箱內。

2.5" GmWRKY-3基因抗病性驗證

分別選取轉入pCAMBIA3300-GmWRKY-3基因與pCAMBIA3300空載體的植株，利用Bar試紙條在蛋白水平進行檢測，結果（圖7）表明，pCAMBIA3300-GmWRKY-3基因與pCAMBIA3300空載體均成功轉入植株。由圖8可知，與對照相比，pCAMBIA3300-GmWRKY-3轉基因植株在莖、葉水平上無明顯差別，但對照組毛狀根明顯少于pCAMBIA3300-GmWRKY-3轉基因植株。接種SCN 15 d后，對大豆根部進行酸性品紅染色，并在20×體式顯微鏡下觀察線蟲侵染情況。圖9-a～圖9-c 為pCAMBIA3300-GmWRKY-3轉基因陽性植株根中大豆胞囊線蟲侵染情況，平均每個側根有 1～5頭大豆胞囊線蟲，根部線蟲數(shù)量平均值為 4.30個/cm；圖9-d為pCAMBIA3300空載體對照植株根部線蟲侵染情況，根部線蟲數(shù)量平均值為6.30個/cm，可見轉基因陽性植株根部線蟲數(shù)量明顯少于對照植株，初步證實，GmWRKY-3基因對大豆胞囊線蟲病有一定的抵抗能力。

3" 討論

大豆最早在我國被發(fā)現(xiàn)，隨后傳入其他國家，

我國大豆胞囊線蟲主要分布于東北地區(qū)和黃淮海地區(qū)。大豆胞囊線蟲于1954年在美國北卡羅來納州第1次被發(fā)現(xiàn)，是一種土傳性病害，在土壤中沒有寄主存在時，在某些儲存條件下可以存活9年［20］。由于鑒別的寄主不同，大豆胞囊線蟲被分為16個生理小種， 我國主要有11個生理小種，其中關于3號生理小種的研究最為廣泛，4號生理小種的侵染能力最強［21］。

目前，選育并種植抗大豆胞囊線蟲病大豆品種是防治SCN最經濟有效的措施，而篩選抗源是抗性品種選育的基礎。Concibido等使用遺傳標記，分析了3種分離大豆F2代種群，并將結果與之前涉及PI 209332的研究結果進行了比較，總結了有關SCN相互作用的廣泛遺傳數(shù)據(jù)，為SCN研究人員和大豆育種者提供了資源［22］。有研究表明，在不同的大豆產區(qū)，SCN也可以在PI 88788、PI 437654、Peking上繁殖，嚴重影響產量及其他農藝性狀。因此迫切需要從不同抗性來源培育穩(wěn)定的抗性品種。近年來，關于大豆胞囊線蟲抗性基因的篩選取得了重大進展，Rhg1（18號染色體）和Rhg4（8號染色體）為主要的抗性位點。來自PI 88788的Rhg1位點是商業(yè)品種中可用的抗性的主要來源，可以最大限度地減少侵染SCN的產量損失。Rhg1包含3個主要的基因，其中1個為t-SNARE，是大豆Rhg1基因抵抗胞囊線蟲的關鍵決斷因素［23］。為抵御細菌、真菌和線蟲等病原體的入侵，許多植物利用囊泡運輸來應對免疫反應，t-SNARE編碼的蛋白是囊泡運輸?shù)闹饕鞍字唬?4］。來自Peking的Rhg4位點主要抵抗3號生理小種和14號生理小種。目前大多數(shù)抗性品種的抗性均來源于PI88788、Peking和PI437654基因型背景，主要抵抗3號生理小種［22，25］，因此，挖掘新的抗性基因至關重要。

同時研究發(fā)現(xiàn)，線蟲的取食位點還受信號激素的影響，Kammerhofer等于2015年發(fā)現(xiàn)，擬南芥中的茉莉酸（JA）對甜菜胞囊線蟲（Heterodera schachtii）取食起著防御作用，乙烯則對線蟲取食起著積極作用，水楊酸（SA）抑制合胞體的形成［26］；Nahar等于2011年在對水稻的研究中發(fā)現(xiàn)，外源施加乙烯利和茉莉酸甲酯（MeJA）后，線蟲在抵消根系防御途徑方面效果較差，從而增強植物對線蟲的抵抗能力［27］。胡巖峰等于2018年研究發(fā)現(xiàn)了茉莉酸（JA）和水楊酸（SA）對SCN的防御效應和作用機制［28］。SCN是植物病原體中最具挑戰(zhàn)性的病原體之一，主要原因在于胞囊在土壤中存活時間長達數(shù)十年，它的卵在胞囊內孵化，可免受破壞性環(huán)境的影響，包括免受迄今為止已知的最具侵略性（不可接受的侵略性）殺蟲劑的侵害［29］。一旦SCN在農田中建立，就很難根除。遠離大豆輪作5～7年是遠遠不夠的，在下一個種植季節(jié)仍然會有未孵化的SCN卵孵化，并且在有利于發(fā)病的環(huán)境范圍內廣泛傳播。SCN侵擾的地理范圍延伸到美國中部和加拿大的主要大豆產區(qū)以及巴西和阿根廷，現(xiàn)在在我國也很常見［30］。

WRKY蛋白可以根據(jù)WRKY結構域的數(shù)量和鋅指基序的類型分為Ⅰ～Ⅲ組。WRKY也可以根據(jù)其基因內含子類型分為2組，即R型或V型。WRKY轉錄因子可以通過結合順式作用元件激活下游基因，涉及生理變化以及對生物和非生物脅迫的反應。WRKY轉錄因子是植物先天免疫系統(tǒng)許多方面的核心組成部分，包括MTI或PTI、ETI、基礎防御和全身性抵抗。WRKY轉錄因子參與種子發(fā)育、種子休眠和發(fā)芽、開花、衰老、代謝途徑以及毛狀體形態(tài)發(fā)生和植物生長。它們主要通過與特定基因啟動子中存在的保守W-box元件結合來發(fā)揮調節(jié)作用。此外，WRKY轉錄因子，尤其是Ⅲ組，可參與對食草動物、病原體和線蟲的反應。現(xiàn)在有許多例子證明，WRKY基因的過表達或敲除對植物防御能力有影響，Grunewald等已初步探索了AtWRKY23在擬南芥中對線蟲的抵抗能力，為后續(xù)從WRKY轉錄因子中挖掘大豆胞囊線蟲抗性基因提供了試驗基礎［31］。

4" 結論

本研究成功克隆了大豆GmWRKY-3基因，該基因的CDS序列長度為663 bp，共編碼220個氨基酸。生物信息學分析結果表明，GmWRKY-3基因編碼的蛋白質屬于不穩(wěn)定蛋白質。二級結構以無規(guī)則卷曲為主，三級結構由1個α-螺旋和5個 β-折疊組成。由啟動子順式作用元件分析結果可知，該基因含有多個光響應元件和多種激素順式作用元件等。為了進一步研究其功能，本研究構建pCAMBIA1302-GmWRKY-3重組質粒，并進行亞細胞定位，在倒置熒光顯微鏡下觀察發(fā)現(xiàn)，GmWRKY-3基因主要在細胞核中行使功能。本研究構建了植物過表達載體pCAMBIA3300-GmWRKY-3，利用發(fā)根農桿菌侵染大豆毛狀根，結果表明，pCAMBIA3300-GmWRKY-3轉基因陽性植株根中線蟲數(shù)量平均值為4.30個/cm；pCAMBIA3300空載體對照植株根部線蟲數(shù)量平均值為6.30個/cm，說明轉基因陽性植株根部雌蟲平均數(shù)明顯少于對照植株，初步證實GmWRKY-3對大豆胞囊線蟲有一定的抵抗能力。結果可為進一步探究在SCN脅迫下，GmWRKY2-3增強大豆抵抗能力的相關研究提供重要的理論依據(jù)。

參考文獻：

［1］曲" 碩，劉" 芳，宋庚晨，等. 大豆抗胞囊線蟲鋅指蛋白基因GmC_2H_2-2like克隆與初步功能分析［J］. 華北農學報，2024，39（3）：195-207.

［2］Kang H，Ko H，Kim D，et al. Occurrence of the white soybean cyst nematode，Heterodera sojae，and the soybean cyst nematode，H. glycines，in Korea［J］. Plant Disease，2021，105（1）：31-33.

［3］孟凡立，于瑾瑤，李春杰，等. 東北地區(qū)大豆胞囊線蟲病發(fā)生和防控技術研究進展［J］. 東北農業(yè)大學學報，2022，53（1）：87-94.

［4］Bradley C A，Allen T W，Sisson A J，et al. Soybean yield loss estimates due to diseases in the United States and Ontario，Canada，from 2015 to 2019［J］. Plant Health Progress，2021，22（4）：483-495.

［5］Peng D L，Jiang R，Peng H，et al. Soybean cyst nematodes：a destructive threat to soybean production in China［J］. Phytopathology Research，2021，3（1）：1-16.

［6］Arjoune Y，Sugunaraj N，Peri S，et al. Soybean cyst nematode detection and management：a review［J］. Plant Methods，2022，18（1）：110.

［7］Niblack T L，Lambert K N，Tylka G L. A model plant pathogen from the Kingdom Animalia：Heterodera glycines，the soybean cyst nematode［J］. Annual Review of Phytopathology，2006，44：283-303.

［8］Niblack T L. Soybean cyst nematode management reconsidered［J］. Plant Disease，2005，89（10）：1020-1026.

［9］Rushton P J，Somssich I E，Ringler P，et al. WRKY transcription factors［J］. Trends in Plant Science，2010，15（5）：247-258.

［10］賈振宇. 煙草轉錄因子NtWRKY65在干旱和低氮脅迫中的功能研究［D］. 鄭州：河南農業(yè)大學，2022.

［11］Zhang C，Wang D D，Yang C H，et al. Genome-wide identification of the potato WRKY transcription factor family［J］. PLoS One，2017，12（7）：e0181573.

［12］王藝橋，趙新杰，牛芳芳，等. 擬南芥中WRKY31轉錄因子的轉錄活性與互作蛋白分析［J］. 生物技術通報，2018，34（5）：101-109.

［13］郭亦非，李" 培，鄒怡婷，等. 水稻OsWRKY78轉錄因子響應鹽脅迫的表達與功能研究［J］. 揚州大學學報（農業(yè)與生命科學版），2019，40（2）：18-24.

［14］郝青婷，高" 偉，閆虎斌，等. 綠豆WRKY基因家族的全基因組鑒定及生物信息學分析［J］. 西北農林科技大學學報（自然科學版），2023，51（5）：59-71，81.

［15］楊雪梅，趙志學，曹小龍，等. 擬南芥WRKY51參與調控RPW8.1介導的基礎防御反應［C］//彭友良，宋寶安. 植物病理科技創(chuàng)新與綠色防控：中國植物病理學會2021年學術年會論文集. 貴陽：中國植物病理學會，2021：463.

［16］Meher J，Sarkar A，Sarma B K. Binding of stress-responsive OsWRKY proteins through WRKYGQK heptapeptide residue with the promoter region of two rice blast disease resistance genes Pi2 and Pi54 is important for development of blast resistance［J］. 3 Biotech，2023，13（9）：294.

［17］劉麗麗，羅明銀，陳銀霞，等. SlWRKYs在番茄抗根結線蟲病過程中初步功能分析與驗證［J/OL］. 分子植物育種，2023：1-11（2023-03-07）［2024-09-01］. http：//kns.cnki.net/kcms/detail/46.1068.S.20230307.1349.002.html.

［18］Wang W J，Li T，Chen Q，et al. Transcriptional regulatory mechanism of a variant transcription factor CsWRKY23 in Citrus fruit resistance to Penicillium digitatum［J］. Food Chemistry，2023，413：135573.

［19］Riggs R D，Schmitt D P.Complete characterization of the race scheme for Heterodera glycines［J］. Journal of Nematology，1988，20（3）：392-395.

［20］Tylka G L，Marett C C.Distribution of the soybean cyst nematode，Heterodera glycines，in the United States and Canada：1954 to 2014［J］. Plant Health Progress，2014，15（2）：85-87.

［21］練" 云，王金社，李海朝，等. 黃淮大豆主產區(qū)大豆胞囊線蟲生理小種分布調查［J］. 作物學報，2016，42（10）：1479-1486.

［22］Concibido V C，Lange D A，Denny R L，et al. Genome mapping of soybean cyst nematode resistance genes in ‘Peking’，PI 90763，and PI 88788 using DNA markers［J］. Crop Science，1997，37（1）：258-264.

［23］Dong J，Zielinski R E，Hudson M E.T-SNAREs bind the Rhg1 α-SNAP and mediate soybean cyst nematode resistance［J］. Plant Journal，2020，104（2）：318-331.

［24］Jones J D，Dangl J L.The plant immune system［J］. Nature，2006，444（7117）：323-329.

［25］Concibido V C，Diers B W，Arelli P R.A decade of QTL mapping for cyst nematode resistance in soybean［J］. Crop Science，2004，44（4）：1121-1131.

［26］Kammerhofer N，Radakovic Z，Regis J M A，et al. Role of stress-related hormones in plant defence during early infection of the cyst nematode Heterodera schachtii in Arabidopsis［J］. New Phytologist，2015，207（3）：778-789.

［27］Nahar K，Kyndt T，de Vleesschauwer D，et al. The jasmonate pathway is a key player in systemically induced defense against root knot nematodes in rice［J］. Plant Physiology，2011，157（1）：305-316.

［28］胡巖峰，華" 萃，陳" 雙. 植物防御激素茉莉酸和水楊酸對大豆胞囊線蟲的誘導抗性研究［R］. 長春：中國科學院東北地理與農業(yè)生態(tài)研究所農業(yè)技術中心，2018-07-26.

［29］Winter Shawn M J，Rajcan I，Shelp B J. Soybean cyst nematode：challenges and opportunities［J］. Canadian Journal of Plant Science，2006，86：25-32.

［30］Tylka G L，Marett C C. Known distribution of the soybean cyst nematode，Heterodera glycines，in the United States and Canada，1954 to 2017［J］. Plant Health Progress，2017，18：167-168.

［31］Grunewald W，Karimi M，Wieczorek K，et al. A role for AtWRKY23 in feeding site establishment of plant-parasitic nematodes［J］. Plant Physiology，2008，148（1）：358-368.