盧苗 李佩 榮鈺瑩 張夢(mèng)涵 賈鵬 欒好安 齊國(guó)輝 張雪梅 董慶龍

摘要：【目的】KNOTTED1 like homeobox（KNOX）蛋白在植物生長(zhǎng)發(fā)育等多個(gè)生物學(xué)過程中發(fā)揮著重要作用。以紫弘富士為材料，分離了多個(gè)蘋果（Malus domestica）KONX基因，研究其結(jié)構(gòu)域、進(jìn)化分析、組織表達(dá)、非生物脅迫響應(yīng)及其與MdOFP相互作用情況?！痉椒ā渴褂肦T-PCR技術(shù)克隆獲得7 個(gè)MdKNOX基因并進(jìn)行生物信息學(xué)分析。利用Array 技術(shù)檢測(cè)MdKNOX基因在蘋果不同組織中的表達(dá)模式。使用RT-qPCR技術(shù)檢測(cè)MdKNOX基因在鹽脅迫和滲透脅迫下的表達(dá)模式。通過Y2H實(shí)驗(yàn)檢測(cè)了MdKNOX蛋白與MdOFP蛋白的互作情況。【結(jié)果】測(cè)序結(jié)果表明，獲得了7 個(gè)KNOX 轉(zhuǎn)錄因子cDNA：MdKNOX1、MdKNOX2、MdKNOX5、MdKNOX10、MdKNOX13、MdKNOX16 和Md-KNOX22（GenBank 登錄號(hào)：MG021644～MG021650）。結(jié)構(gòu)域分析表明，獲得的7 個(gè)MdKNOX蛋白序列均含有MEINOX、HD和ELK 結(jié)構(gòu)域。進(jìn)化分析表明，MdKNOX1、MdKNOX2 和MdKNOX5 屬于KNOX Ⅱ亞組；MdKNOX10、MdKNOX13、MdKNOX16 和MdKNOX22 屬于KNOXⅠ亞組。啟動(dòng)子順式作用元件分析結(jié)果表明，7 個(gè)MdKNOX基因啟動(dòng)子上包含多個(gè)順式作用元件。Array 分析結(jié)果顯示，MdKNOX基因具有不同的組織表達(dá)模式。實(shí)時(shí)熒光定量PCR分析結(jié)果表明，鹽脅迫處理下，MdKNOX13 的相對(duì)表達(dá)水平上調(diào)，而MdKNOX1、MdKNOX2 和MdKNOX5 的相對(duì)表達(dá)水平下調(diào)；滲透脅迫處理下，MdKNOX2 的轉(zhuǎn)錄水平下調(diào)。酵母雙雜交結(jié)果顯示，MdKNOX1/22 蛋白能與MdOFP6 蛋白相互作用，MdKNOX5 蛋白能與多個(gè)MdOFP 蛋白相互作用，且MdKNOX5 蛋白與MdOFP 蛋白相互作用，HD區(qū)域是互作必需的?！窘Y(jié)論】這些結(jié)果為蘋果KNOX轉(zhuǎn)錄因子在生長(zhǎng)、發(fā)育和逆境下生物學(xué)功能的解析、調(diào)控網(wǎng)絡(luò)的構(gòu)建提供了強(qiáng)有力的理論基礎(chǔ)和參考。

關(guān)鍵詞：蘋果；KNOX轉(zhuǎn)錄因子；基因克??；表達(dá)分析；蛋白互作

中圖分類號(hào)：S661.1 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1009-9980（2023）05-0841-11

自從在玉米中發(fā)現(xiàn)knotted1（kn1）基因以來，已在多個(gè)植物中發(fā)現(xiàn)KNOX（KNOTTED1 like homeobox）基因[1]。KNOX 轉(zhuǎn)錄因子屬于TALE（threeamino-acid loop extension）基因家族中的亞家族，在植物中包含多個(gè)家族成員，其顯著的特征是含有同源異型盒結(jié)構(gòu)域（homeodomain，HD），此外還含有KNOX1、KNOX2 和ELK 結(jié)構(gòu)域[2]。其中HD 結(jié)構(gòu)域參與蛋白互作的形成和與DNA 的結(jié)合[3]。KNOX1 和KNOX2 結(jié)構(gòu)域被稱為MEINOX 結(jié)構(gòu)域，對(duì)轉(zhuǎn)基因植株表型的變化起到重要作用[3]。

ELK結(jié)構(gòu)域編碼核定位信號(hào)，在與其他蛋白互作和轉(zhuǎn)錄抑制過程中也起到一定作用[3]。多項(xiàng)研究表明，KNOX 蛋白能與BLH 蛋白相互作用形成不同組合的異質(zhì)二聚體，某些異質(zhì)二聚體還可與MADS-box 或者OVATE family proteins（OFP）轉(zhuǎn)錄因子相互作用形成功能復(fù)合體，調(diào)控植物的生長(zhǎng)和發(fā)育過程[4]?；贙NOX 基因的序列相似性、結(jié)構(gòu)特征、系統(tǒng)進(jìn)化關(guān)系以及表達(dá)模式，可將KNOX蛋白分為2 個(gè)亞家族：Class Ⅰ（KNOX Ⅰ）和Class Ⅱ（KNOX Ⅱ）[2，5-8]。進(jìn)一步的研究表明，在蒺藜苜蓿（Medicago truncatula L.）、番茄（Solanum lycopersicumL.）和擬南芥（Arabidopsis thaliana L. Heynh.）中鑒定到一個(gè)新的KNOX蛋白亞家族KNOXM，與KNOX Ⅰ和KNOX Ⅱ亞家族不同的是其本身缺少ELK 和HD結(jié)構(gòu)域[9-11]。KNOX Ⅰ能夠進(jìn)一步分為3 個(gè)亞組：STM-like、KNAT2/6- like 和KNAT1/BPlike，而KNOX Ⅱ 分為2 個(gè)亞組：KNAT7- like 和KANT3/4/5-like[1]。通過擬南芥突變體研究發(fā)現(xiàn)，KNOX基因不僅調(diào)控植物細(xì)胞增殖和組織分化，還參與調(diào)節(jié)細(xì)胞分裂素和赤霉素信號(hào)途徑[1，7，12]。擬南芥KNOX Ⅰ亞組包含4 個(gè)基因：SHOOTMERISTMELESS（STM）以及KNAT1/2/6。其中，STM 基因在莖頂端分生組織（shoot apical meristem，SAM）早期胚胎發(fā)生期間表達(dá)，對(duì)SAM 形成的起始起到重要作用[12]。KNAT1 和KNAT6 在SAM 功能的形成和花序發(fā)育過程中扮演著重要角色[13-14]。KNAT2對(duì)花模型的調(diào)控起到重要作用[15]。相比KNOX Ⅰ的廣泛研究，KNOX Ⅱ已知突變體缺乏表型，研究相對(duì)較少，在擬南芥中發(fā)現(xiàn)KNAT7 基因在調(diào)控次級(jí)細(xì)胞壁合成轉(zhuǎn)錄網(wǎng)絡(luò)中起到一定作用[16-17]。最近研究表明，KNAT3/4/5 這3 個(gè)基因與KNOX Ⅰ成員起到相反的作用，功能冗余地促進(jìn)結(jié)瘤器官的分化[8]。此外，MtKNAT3/4/5-like 基因參與共生根瘤的發(fā)育，并可調(diào)節(jié)豆科植物根瘤細(xì)胞分裂素生物合成[18]。在蘋果（Malus domestica Borkh.）中發(fā)現(xiàn)Md-KNOX15 可通過調(diào)節(jié)赤霉素水平調(diào)控蘋果株高和開花[19]。MdKNOX19 能夠靶向ABI5 調(diào)節(jié)ABA 敏感性和種子萌發(fā)[20]。

Jia 等[21]通過對(duì)蘋果基因組數(shù)據(jù)庫(kù)序列比對(duì)，發(fā)現(xiàn)蘋果基因組中含有22 個(gè)MdKNOX基因。去除已克隆的MdKNOX15 和MdKNOX19 基因，筆者獲得了7 個(gè)蘋果MdKNOX基因，進(jìn)一步對(duì)它們進(jìn)行了生物信息學(xué)、組織器官表達(dá)、鹽和滲透脅迫表達(dá)以及蛋白相互作用分析，為蘋果KNOX轉(zhuǎn)錄因子在生長(zhǎng)、發(fā)育和逆境下生物學(xué)功能的解析、調(diào)控網(wǎng)絡(luò)的構(gòu)建提供了強(qiáng)有力的理論基礎(chǔ)和參考。

1 材料和方法

1.1 植物材料與處理方法

鹽和滲透脅迫表達(dá)分析供試材料為嘎拉組培苗。嘎拉組培苗培養(yǎng)條件和處理方法參照李慧峰等[22]的文獻(xiàn)描述進(jìn)行操作。將嘎拉組培苗放在液氮中速凍以備提取RNA。

1.2 MdKNOX基因克隆

MdKNOX 基因克隆所采用的模板為完全展開的紫弘富士蘋果葉片，采用改良熱硼酸法進(jìn)行RNA的提取。嘎拉組培苗RNA的提取采用QIAGEN 公司的RNA提取試劑盒RNeasy Plant Mini Kit （貨號(hào)：74903）。cDNA 的合成采用TaRaKa 公司的反轉(zhuǎn)錄試劑盒PrimeScriptTM 1st Strand cDNA Synthesis Kit（貨號(hào)：6110A）。依據(jù)蘋果MdKNOX核苷酸序列使用Primer 3 在線軟件（https：//primer3.ut.ee/）設(shè)計(jì)Md-KNOX 基因特異性引物（表1）進(jìn)行RT-PCR 擴(kuò)增。PCR 反應(yīng)條件參照Dong 等[23]文獻(xiàn)描述進(jìn)行操作。對(duì)PCR反應(yīng)液進(jìn)行膠回收后，將目的基因克隆片段導(dǎo)入到pMD18-T 克隆載體上轉(zhuǎn)化Escherichia coli DH5α 感受態(tài)細(xì)胞，涂板后篩選陽(yáng)性克隆，把搖菌送到公司進(jìn)行測(cè)序。

1.3 蘋果KNOX轉(zhuǎn)錄因子蛋白序列、進(jìn)化關(guān)系和順式作用元件分析

使用蘋果、擬南芥和水稻全長(zhǎng)KNOX氨基酸序列，采用MEGA 6 軟件（http：//www.megasoftware.net）進(jìn)行進(jìn)化樹分析。軟件具體參數(shù)：NJ方法、pairwisedeletion、Poission correction 和bootstrap（1000repeat）[24]。利用序列分析軟件DNAMAN 6.0.3.99 分析蘋果KNOX全長(zhǎng)蛋白序列。MdKNOX基因啟動(dòng)子上含有的順式作用元件利用PlantCARE 網(wǎng)站（http：//bioinformatics.psb.ugent.be/webtools/plantcare/html/）進(jìn)行分析。

1.4 MdKNOX的表達(dá)分析

根據(jù)NCBI 網(wǎng)站中GEO數(shù)據(jù)庫(kù)GSE42873 生成蘋果MdKNOX 組織器官表達(dá)數(shù)據(jù)。該轉(zhuǎn)錄譜（Array）包含10 個(gè)不同蘋果基因型的16 個(gè)器官組織數(shù)據(jù)。通過TIGR MeV v4.8.1 軟件進(jìn)行表達(dá)熱圖的繪制。MdKNOX 基因在鹽和滲透脅迫下的表達(dá)分析采用實(shí)時(shí)熒光定量PCR（RT-qPCR）進(jìn)行檢測(cè)，使用MdMDH基因作為蘋果內(nèi)參基因，相對(duì)表達(dá)水平結(jié)果采用2–ΔΔCT法進(jìn)行計(jì)算[25]。RT-qPCR 反應(yīng)體系和程序參照李慧峰等[22]的文獻(xiàn)描述進(jìn)行操作。Md-KNOX基因的熒光定量引物見表1。

1.5 酵母雙雜交（Y2H）實(shí)驗(yàn)

Y2H方法參照Clontech 公司公布的方法進(jìn)行操作。將MdKNOX1、MdKNOX5 和MdKNOX22全長(zhǎng)編碼框、MdKNOX5 刪除MEINOX 結(jié)構(gòu)域片段和Md-KNOX5 刪除HD 結(jié)構(gòu)域片段插入到pGAD424（GAL4 activation domain，AD）捕獲載體中；將MdOFP1/4/5/6/11/13/14/16/19/20 全長(zhǎng)編碼框克隆到pGBT9（GAL4 DNA-binding domain，BD）誘餌載體中[26]。將不同融合載體以AD和BT配對(duì)形式分別轉(zhuǎn)化酵母菌株Y2HGold，然后涂布在二缺培養(yǎng)基SD/-Trp/-Leu和四缺培養(yǎng)基SD/-Trp/-Leu/-His/-Ade+X-α-Gal 中進(jìn)行培養(yǎng)，通過觀察菌體是否變藍(lán)判斷蛋白是否存在相互作用。采用pGBT9-MdOFP 和pGBT9-MdWKRY52 融合載體與pGAD424 空載體質(zhì)粒共轉(zhuǎn)化的酵母菌作為陰性對(duì)照；采用pGBT9-Md-WRKY52 和pGAD424-MdVQ10 融合載體共轉(zhuǎn)化的酵母菌作為陽(yáng)性對(duì)照[23]。

1.6 數(shù)據(jù)分析

MdKNOX 基因在鹽和滲透脅迫處理下的表達(dá)水平采用IBM SPSS Statistics v.20 軟件中One-wayANOVA 方法及Duncan 檢驗(yàn)（p＜0.05）進(jìn)行差異顯著性分析。

2 結(jié)果與分析

2.1 蘋果MdKNOX基因的克隆

使用紫弘富士葉片的cDNA 作為PCR 反應(yīng)模板，進(jìn)行RT-PCR擴(kuò)增。測(cè)序結(jié)果顯示，7 個(gè)蘋果Md-KNOX基因克隆成功。根據(jù)前人對(duì)MdKNOX基因家族鑒定的研究結(jié)果，對(duì)7 個(gè)蘋果MdKNOX基因進(jìn)行命名，各個(gè)基因的GenBank 登錄號(hào)、基因組ID 號(hào)、染色體定位、開放閱讀框、分子質(zhì)量以及等電點(diǎn)信息詳見表2。DNAMAN 軟件分析結(jié)果顯示，7 個(gè)Md-KNOX轉(zhuǎn)錄因子含有典型的KNOX1、KNOX2、ELK和HD結(jié)構(gòu)域（圖1）。

2.2 蘋果KNOX轉(zhuǎn)錄因子的進(jìn)化分析

利用進(jìn)化樹分析軟件MEGA6 對(duì)9 個(gè)蘋果Md-KNOX（7 個(gè)本文克隆的MdKNOX、MdKNOX15 和MdKNOX19）、9 個(gè)擬南芥KNOX和11 個(gè)水稻（Oryzasativa L.）KNOX 全長(zhǎng)蛋白序列進(jìn)行進(jìn)化分析。圖2 結(jié)果顯示，KNOX成員被清楚地分為3 個(gè)亞家族（KNOX Ⅰ、KNOX Ⅱ和KNOX M），其中KNOX Ⅰ中可進(jìn)一步分為3 個(gè)亞組（STM-like、KNAT2/6-like和KNAT1/BP- like），KNOX Ⅱ 分為2 個(gè)亞組（KNAT7-like 和KANT3/4/5-like）。STM-like 亞組中包含MdKNOX10 和MdKNOX22；KNAT2/6- like 亞組中包含MdKNOX13、MdKNOX15 和Md-KNOX16；KNAT7- like 亞組中包含MdKNOX2 和MdKNOX5；KNAT3/4/5-like 亞組中包含MdKNOX1和MdKNOX19。

2.3 蘋果MdKNOX啟動(dòng)子序列分析

下載已克隆MdKNOX 基因翻譯起始位點(diǎn)上游1500 bp 序列，獲得了7 個(gè)MdKNOX基因的啟動(dòng)子序列。通過PlantCARE數(shù)據(jù)庫(kù)分析啟動(dòng)子序列上的順式作用元件，除了光響應(yīng)元件，還有多個(gè)逆境和激素響應(yīng)元件（圖3）。在7 個(gè)MdKNOX基因的啟動(dòng)子上總共存在13 種不同類型的順式作用元件，分別是對(duì)響應(yīng)逆境的低氧、低溫、熱、干旱、機(jī)械創(chuàng)傷、病原菌和激素響應(yīng)的茉莉酸甲酯、水楊酸、生長(zhǎng)素、赤霉素、乙烯和ABA 等順式作用元件（圖3）。這些結(jié)果表明，MdKNOX啟動(dòng)子上的順式作用元件可能在蘋果生長(zhǎng)和發(fā)育以及逆境脅迫響應(yīng)中起到重要的作用。

2.4 蘋果MdKNOX基因的表達(dá)分析

利用NCBI 網(wǎng)站的GEO 數(shù)據(jù)庫(kù)下載蘋果16 個(gè)不同組織的轉(zhuǎn)錄譜（GSE42873），檢測(cè)MdKNOX 在不同組織中的表達(dá)情況（圖4）。結(jié)果顯示，Md-KNOX基因在被檢測(cè)的組織中均有表達(dá)。其中，Md-KNOX1、MdKNOX2 和MdKNOX5 在被檢測(cè)的組織中有相對(duì)較高的表達(dá)水平；在葉_M49（Leaf M49）、花_M74（Flower M74）、果實(shí)_M20 花后100 d [FruitM20（100DAM）] 和果實(shí)_M20 收獲期[Fruit M20（harvest）] 組織中，MdKNOX10 、MdKNOX13、Md-KNOX16 和MdKNOX22 基因有相對(duì)較高的表達(dá)水平（圖4）。

進(jìn)一步利用熒光實(shí)時(shí)定量PCR分析鹽和甘露醇處理下嘎拉組培苗中MdKNOX基因的表達(dá)情況，結(jié)果顯示，在正常條件下，MdKNOX的相對(duì)表達(dá)水平無明顯變化（圖5-A）；在NaCl處理?xiàng)l件下，MdKNOX13與對(duì)照相比，相對(duì)表達(dá)水平升高，在處理48 h后相對(duì)表達(dá)水平達(dá)到對(duì)照的1.96 倍，MdKNOX1、MdKNOX2和MdKNOX5 的相對(duì)表達(dá)量與對(duì)照相比下降，Md-KNOX1 在處理24 h 后相對(duì)表達(dá)水平為對(duì)照的0.28倍；MdKNOX2 和MdKNOX5 在處理48 h 后，相對(duì)表達(dá)水平分別是對(duì)照的0.32和0.29倍（圖5-B）；在甘露醇處理?xiàng)l件下，MdKNOX2 與對(duì)照相比，表達(dá)水平降低，在48 h時(shí)為對(duì)照的0.47倍（圖5-C）。

2.5 MdKNOX蛋白與MdOFP蛋白相互作用

筆者前期研究發(fā)現(xiàn)蘋果基因組中存在26 個(gè)MdOFP 基因家族成員，并且成功克隆了10 個(gè)MdOFP 基因（MdOFP1/4/5/6/11/13/14/16/19/20）[26]。本研究利用酵母雙雜交試驗(yàn)，檢測(cè)了MdKNOX 蛋白與這10 個(gè)MdOFP 蛋白的相互作用情況。3 個(gè)MdKNOX 全長(zhǎng)cDNA 序列也融合到AD 捕獲載體上，同時(shí)10 個(gè)MdOFP全長(zhǎng)cDNA序列融合到BD誘餌載體中[26]。將AD-MdKNOX 和BD-MdOFP 融合載體共轉(zhuǎn)化到酵母細(xì)胞，通過觀察β-galactosidase 活性來測(cè)試LacZ 報(bào)告基因的表達(dá)。如圖6-A所示，作為陰性對(duì)照，10 個(gè)BD-MdOFP 融合誘餌載體與空AD捕獲載體在四缺培養(yǎng)基上不變藍(lán)，說明這10 個(gè)MdOFP 蛋白沒有轉(zhuǎn)錄自激活活性。MdKNOX1 和MdKNOX22 蛋白能與MdOFP6 蛋白相互作用，MdKNOX5 蛋白能與MdOFP1、MdOFP4、MdOFP14和MdOFP16 蛋白相互作用（圖6-A）。為了檢測(cè)MdKNOX 蛋白哪一段區(qū)域?qū)τ贛dOFP 蛋白相互作用是必需的，構(gòu)建了MdKNOX10 刪除MEINOX結(jié)構(gòu)域捕獲載體（AD-MdKNOX10△MEINOX）和MdKNOX10 刪除HD 結(jié)構(gòu)域捕獲載體（AD-Md-KNOX10△HD）（圖6-B）。如圖6-C 所示，作為陽(yáng)性對(duì)照，MdVQ10 蛋白和MdWRKY52 蛋白相互作用變藍(lán)[23]。MdKNOX5 △ MEINOX 與MdOFP1、MdOFP4、MdOFP14 和MdOFP16 蛋白相互作用，而MdKNOX5△HD與這些蛋白失去相互作用能力，說明MdKNOX5 與MdOFP1、MdOFP4、MdOFP14 和MdOFP16 相互作用，HD區(qū)域是互作必需的。

3 討論

根據(jù)植物KNOX 蛋白的高度保守結(jié)構(gòu)域和蘋果基因組數(shù)據(jù)庫(kù)公布的數(shù)據(jù)，Jia 等[22]對(duì)蘋果KNOX基因進(jìn)行了基因組范圍內(nèi)的家族鑒定，篩選出了22個(gè)MdKNOX 基因，隨后對(duì)MdKNOX15 和Md-KNOX19轉(zhuǎn)錄因子進(jìn)行了功能鑒定[19-20]。本研究根據(jù)前人研究結(jié)果，克隆了7個(gè)MdKNOX基因，對(duì)結(jié)構(gòu)域、進(jìn)化分類和啟動(dòng)子順式作用元件進(jìn)行了分析，檢測(cè)了組織器官表達(dá)、非生物脅迫應(yīng)答以及與MdOFP蛋白互作模式，表明這些MdKNOX基因在蘋果調(diào)控生長(zhǎng)、發(fā)育以及應(yīng)對(duì)非生物脅迫過程中起到不同的作用。

研究表明，KNOX蛋白能夠與BLH形成異質(zhì)二聚體或者與OFP轉(zhuǎn)錄因子相互作用或者它們3 者形成功能復(fù)合體參與調(diào)控植物生長(zhǎng)和發(fā)育[12，16，27]。例如，KANT3 與BLH1 形成異質(zhì)二聚體，AtOFP5 介導(dǎo)它們活性的抑制參與到胚囊的正常發(fā)育以及細(xì)胞命運(yùn)的決定[28]；KNAT7 與OFP4 相互作用，增強(qiáng)KNAT7的轉(zhuǎn)錄抑制活性調(diào)控次生細(xì)胞壁的形成[29]。此外，KNAT7 蛋白還可與BLH6 蛋白相互作用，通過抑制homeodomain-leucine zipper transcription factor REVOLUTA/INTERFASCICULAR FIBERLESS1 （REV/IFL1）的表達(dá)，進(jìn)而調(diào)控次生細(xì)胞壁的形成[16]。進(jìn)一步研究表明，KNAT7-BLH6 異質(zhì)二聚體還可與AtOFP1 和AtOFP4 相互作用形成功能復(fù)合體，參與到次生細(xì)胞壁的形成[30]；水稻OsKNAT7 可與Os-OFP2、BLH6-like 和BLH-like2 相互作用調(diào)控維管發(fā)育[31]。水稻KNOX 蛋白OSH15 與BEL-like home-odomain 蛋白SH5 相互作用形成二聚體，通過抑制木質(zhì)素生物合成基因改善落粒性[32]；玉米（Zea maysL.）KNOTTED1 與BLH12 和BLH14 相互作用，在莖中的脈管和節(jié)間結(jié)構(gòu)中起到重要作用[27]。此外，KNAT7 可與MYB75 相互作用在莖和種皮中對(duì)細(xì)胞壁的形成起到重要的調(diào)控作用[33-34]。在本文中，通過酵母雙雜交實(shí)驗(yàn)，發(fā)現(xiàn)MdKNOX1 和MdKNOX22蛋白能與MdOFP6 蛋白相互作用，MdKNOX5 蛋白能與多個(gè)MdOFP蛋白相互作用。MdKNOX蛋白與MdOFP 之間的這種不同的相互作用可能影響到蘋果可能的BLH-KNOX異質(zhì)二聚體的活性，從而改變它們所調(diào)控靶基因的表達(dá)水平，進(jìn)而調(diào)控蘋果的生長(zhǎng)和發(fā)育，但還需進(jìn)一步的試驗(yàn)來分析這種潛在的相關(guān)分子機(jī)制。

盡管KNOX蛋白能夠調(diào)控植物生長(zhǎng)和發(fā)育等多個(gè)方面，但其對(duì)非生物和生物脅迫的響應(yīng)機(jī)制還研究甚少。最近研究表明，在干旱、鹽和冷處理?xiàng)l件下，鷹嘴豆（Cicer arietinum Linn.）根和莖中部分KNOX 基因響應(yīng)脅迫處理應(yīng)答[35]。在干旱處理下，大豆少數(shù)KNOX基因受脅迫響應(yīng)；在病原菌侵染下，大豆KNOX 基因Glyma17g14180、Glyma04g06810、Glyma09g01000 和Glyma14gg37550 受到不同的響應(yīng)[35]。在本文中，MdKNOX13 受鹽脅迫誘導(dǎo)；Md-KNOX1、MdKNOX2 和MdKNOX5 受鹽脅迫下調(diào)；MdKNOX2 受甘露醇脅迫下調(diào)，并且這些MdKNOX基因啟動(dòng)子上含有多個(gè)逆境脅迫順式作用元件，表明MdKNOX蛋白可能在應(yīng)對(duì)蘋果非生物脅迫中起到一定的作用。相信隨著關(guān)于KNOX基因研究的深入，其在非生物脅迫中作用和機(jī)制也會(huì)日益清晰。

4 結(jié)論

本研究克隆獲得了蘋果中7 個(gè)MdKNOX基因，發(fā)現(xiàn)均含有保守的KNOX1、KNOX2、ELK和HD結(jié)構(gòu)域。通過進(jìn)化分析表明7 個(gè)MdKNOX轉(zhuǎn)錄因子分別屬于KNOX Ⅰ亞組和KNOX Ⅱ亞組。Array 結(jié)果發(fā)現(xiàn)，MdKNOX基因在16 個(gè)組織中均有不同的表達(dá)水平。RT-qPCR 結(jié)果表明，MdKNOX13 的相對(duì)表達(dá)水平受到鹽脅迫誘導(dǎo)，而MdKNOX1、MdKNOX2和MdKNOX5 的相對(duì)表達(dá)水平下調(diào)；MdKNOX2 的轉(zhuǎn)錄水平在滲透脅迫處理下下調(diào)表達(dá)。酵母雙雜交分析表明，MdKNOX蛋白能夠與多個(gè)MdOFP 蛋白相互作用，且HD區(qū)域是互作必需的。這些結(jié)果為蘋果KNOX轉(zhuǎn)錄因子在生長(zhǎng)、發(fā)育和逆境下生物學(xué)功能的解析、調(diào)控網(wǎng)絡(luò)的構(gòu)建提供了強(qiáng)有力的理論基礎(chǔ)和參考。

參考文獻(xiàn)References：

[1] HAMANT O，PAUTOT V. Plant development：A tale story[J].

Comptes Rendus Biologies，2010，333（4）：371-381.

[2] MUKHERJEE K，BROCCHIERI L，BURGLIN T R. A comprehensive

classification and evolutionary analysis of plant homeobox

genes[J]. Molecular Biology and Evolution，2009，26（12）：

2775-2794.

[3] SAKAMOTO T，NISHIMURA A，TAMAOKI M，KUBA M，

TANAKA H，IWAHORI S，MATSUOKA M. The conserved

KNOX domain mediates specificity of tobacco KNOTTED1-

type homeodomain proteins[J]. The Plant Cell，1999，11（8）：

1419-1431.

[4] ARNAUD N，PAUTOT V. Ring the BELL and tie the KNOX：

Roles for TALEs in gynoecium development[J]. Frontiers in

Plant Science，2014，5：93.

[5] B?RGLIN T R. Analysis of TALE superclass homeobox genes

（MEIS，PBC，KNOX，Iroquois，TGIF） reveals a novel domain

conserved between plants and animals[J]. Nucleic Acids Research，

1997，25（21）：4173-4180.

[6] BELLAOUI M，PIDKOWICH M S，SAMACH A，KUSHALAPPA

K，KOHALMI S E，MODRUSAN Z，CROSBY W L，

HAUGHN G W. The Arabidopsis BELL1 and KNOX TALE homeodomain

proteins interact through a domain conserved between

plants and animals[J]. The Plant Cell，2001，13（11）：2455-

2470.

[7] HAKE S，SMITH H M S，HOLTAN H，MAGNANI E，MELE G，

RAMIREZ J. The role of KNOX genes in plant development[J].

Annual Review of Cell and Developmental Biology，2004，20：

125-151.

[8] FURUMIZU C，ALVAREZ J P，SAKAKIBARA K，BOWMAN

J L. Antagonistic roles for KNOX1 and KNOX2 genes in patterning

the land plant body plan following an ancient gene duplication[

J]. PLoS Genetics，2015，11（2）：e1004980.

[9] KIMURA S，KOENIG D，KANG J L，YOONG F Y，SINHA N.

Natural variation in leaf morphology results from mutation of a

novel KNOX gene[J]. Current Biology，2008，18（9）：672-677.

[10] MAGNANI E，HAKE S. KNOX lost the OX：The Arabidopsis

KNATM gene defines a novel class of KNOX transcriptional regulators

missing the homeodomain[J]. The Plant Cell，2008，20

（4）：875-887.

[11] PENG J L，YU J B，WANG H L，GUO Y Q，LI G M，BAI G H，

CHEN R J. Regulation of compound leaf development in Medicago

truncatula by fused compound leaf1，a class M KNOX gene[J].

The Plant Cell，2011，23（11）：3929-3943.

[12] HAY A，TSIANTIS M. KNOX genes：Versatile regulators of

plant development and diversity[J]. Development，2010，137

（19）：3153-3165.

[13] DOUGLAS S J，CHUCK G，DENGLER R E，PELECANDA L，

RIGGS C D. KNAT1 and ERECTA regulate inflorescence architecture

in Arabidopsis[J]. The Plant Cell，2002，14（3）：547-558.

[14] RAGNI L，BELLES-BOIX E，G?NL M，PAUTOT V. Interaction

of KNAT6 and KNAT2 with BREVIPEDICELLUS and PENNYWISE

in Arabidopsis inflorescences[J]. The Plant Cell，2008，

20（4）：888-900.

[15] LI Y，PI L M，HUANG H，XU L. ATH1 and KNAT2 proteins

act together in regulation of plant inflorescence architecture[J].

Journal of Experimental Botany，2012，63（3）：1423-1433.

[16] LIU Y Y，YOU S J，TAYLOR-TEEPLES M，LI W L，SCHUETZ

M，BRADY S M，DOUGLAS C J. BEL1- LIKE HOMEODOMAIN6

and KNOTTED ARABIDOPSIS THALIANA7 interact

and regulate secondary cell wall formation via repression of

REVOLUTA[J]. The Plant Cell，2015，26（12）：4843-4861.

[17] HE J B，ZHAO X H，DU P Z，ZENG W，BEAHAN C T，

WANG Y Q，LI H L，BACIC A，WU A M. KNAT7 positively

regulates xylan biosynthesis by directly activating IRX9 expression

in Arabidopsis[J]. Journal of Integrative Plant Biology，

2018，60（6）：514-528.

[18] DI GIACOMO E，LAFFONT C，SCIARRA F，IANNELLI MA，

FRUGIER F，F(xiàn)RUGIS G. KNAT3/4/5-like class 2 KNOX transcription

factors are involved in Medicago truncatula symbiotic

nodule organ development[J]. New Phytologist，2017，213（2）：

822-837.

[19] JIA P，XING L B，ZHANG C G，CHEN H，LI Y M，ZHANG D，

MAJ J，ZHAO C P，HAN MY，REN X L，AN N. MdKNOX15，

a class I knotted-like transcription factor of apple，controls flowering

and plant height by regulating GA levels through promoting

the MdGA2ox7 transcription[J]. Environmental and Experimental

Botany，2021，185：104411.

[20] JIA P，XING L B，ZHANG C G，ZHANG D，MA J J，ZHAO C

P，HAN M Y，REN X L，AN N. MdKNOX19，a class Ⅱ knotted-

like transcription factor of apple，plays roles in ABA signalling/

sensitivity by targeting ABI5 during organ development[J].

Plant Science，2021，302：110701.

[21] JIA P，ZHANG C G，XING L B，LI Y M，SHAH K，ZUO X Y，

ZHANG D，AN N，HAN MY，REN X L. Genome-wide identification

of the MdKNOX gene family and characterization of its

transcriptional regulation in Malus domestica[J]. Frontiers in

Plant Science，2020，11：128.

[22] 李慧峰，張文芹，董慶龍，王小非，冉昆. 蘋果生長(zhǎng)素響應(yīng)因子

（MdARF）基因克隆與表達(dá)分析[J]. 果樹學(xué)報(bào)，2018，35（10）：

1170-1181.

LI Huifeng，ZHANG Wenqin，DONG Qinglong，WANG Xiaofei，

RAN Kun. Cloning，sequencing and expression analysis of

auxin response factors （MdARF） in apple[J]. Journal of Fruit

Science，2018，35（10）：1170-1181.

[23] DONG Q L，ZHAO S，DUAN D Y，TIAN Y，WANG Y P，

MAO K，ZHOU Z S，MA F W. Structural and functional analyses

of genes encoding VQ proteins in apple[J]. Plant Science，

2018，272：208-219.

[24] TAMURA K，STECHER G，PETERSON D，F(xiàn)ILIPSKI A，KUMAR

S. MEGA6：Molecular evolutionary genetics analysis version

6.0[J]. Molecular Biology and Evolution，2013，30（12）：

2725-2729.

[25] LIVAK K J，SCHMITTGEN T D. Analysis of relative gene expression

data using real- time quantitative PCR and the 2- ΔΔCT

method[J]. Methods，2001，25（4）：402-408.

[26] LI H F，DONG Q L，ZHAO Q，RAN K. Genome-wide identification，

expression profiling，and protein-protein interaction properties

of ovate family proteins in apple[J]. Tree Genetics & Genomes，

2019，15（3）：45.

[27] TSUDA K，ABRAHAM-JUAREZ M J，MAENO A，DONG Z

B，AROMDEE D，MEELEY R，SHIROISHI T，NONOMURA

K I，HAKE S. KNOTTED1 cofactors，BLH12 and BLH14，regulate

internode patterning and vein anastomosis in maize[J]. The

Plant Cell，2017，29（5）：1105-1118.

[28] PAGNUSSAT G C，YU H J，SUNDARESAN V. Cell-fate switch

of synergid to egg cell in Arabidopsis eostre mutant embryo sacs

arises from misexpression of the BEL1-like homeodomain gene

BLH1[J]. The Plant Cell，2007，19（11）：3578-3592.

[29] LI E Y，WANG S C，LIU Y Y，CHEN J G，DOUGLAS C J.

OVATE FAMILY PROTEIN4 （OFP4） interaction with KNAT7

regulates secondary cell wall formation in Arabidopsis thaliana[

J]. The Plant Journal，2011，67（2）：328-341.

[30] LIU Y Y，DOUGLAS C J. A role for OVATE FAMILY PROTEIN1

（OFP1） and OFP4 in a BLH6-KNAT7 multi-protein complex

regulating secondary cell wall formation in Arabidopsis

thaliana[J]. Plant Signaling & Behavior， 2015， 10（7）：

e1033126.

[31] SCHMITZ A J，BEGCY K，SARATH G，WALIA H. Rice Ovate

Family Protein 2 （OFP2） alters hormonal homeostasis and vasculature

development[J]. Plant Science，2015，241：177-188.

[32] YOON J，CHO L H，ANTT H W，KOH H J，AN G. KNOX protein

OSH15 induces grain shattering by repressing lignin biosynthesis

genes[J]. Plant Physiology，2017，174（1）：312-325.

[33] BHARGAVA A，MANSFIELD S D，HALL H C，DOUGLAS C

J，ELLIS B E. MYB75 functions in regulation of secondary cell

wall formation in the Arabidopsis inflorescence stem[J]. Plant

Physiology，2010，154（3）：1428-1438.

[34] BHARGAVA A，AHAD A，WANG S C，MANSFIELD S D，

HAUGHN G W，DOUGLAS C J，ELLIS B E. The interacting

MYB75 and KNAT7 transcription factors modulate secondary

cell wall deposition both in stems and seed coat in Arabidopsis[J].

Planta，2013，237（5）：1199-1211.

[35] BHATTACHARJEE A，GHANGAL R，GARG R，JAIN M. Genome-

wide analysis of homeobox gene family in legumes：Identification，

gene duplication and expression profiling[J]. PLoS

One，2015，10（3）：e0119198.