木豆OFP蛋白基因家族鑒定與表達分析

摘要：了解OFP卵型家族蛋白（OVATE family proteins）在木豆中的特征，為木豆OFP基因的應(yīng)用提供參考。利用生物信息學(xué)方法對木豆OFP蛋白基因進行全基因組鑒定及表達模式分析。結(jié)果表明，木豆含有21個OFP蛋白基因，家族成員間理化特性差異不大，均屬于不穩(wěn)定親水性蛋白，且均含有保守基序；進化樹分析中木豆和擬南芥二者的OFP蛋白相互交錯聚類在一起，親緣關(guān)系較近的OFP蛋白可能具有相似的功能；啟動子順式響應(yīng)元件分析表明木豆OFP基因可能在逆境脅迫響應(yīng)、激素調(diào)節(jié)和生長發(fā)育方面發(fā)揮重要作用；蛋白二級和三級結(jié)構(gòu)均以α-螺旋和無規(guī)則卷曲結(jié)構(gòu)為主要組成成分；組織特異性表達表明OFP基因在頂芽中的表達量較大，在根和葉中的表達量較小，具有組織表達差異；啟動子順式響應(yīng)元件結(jié)合熱脅迫表達分析表明，高溫脅迫下CcOFP1表達量上調(diào)，CcOFP9、CcOFP11表達量下調(diào)。研究結(jié)果為探索OFP蛋白基因在木豆生長發(fā)育和非生物脅迫響應(yīng)方面的生物學(xué)功能提供依據(jù)。

關(guān)鍵詞：木豆；OFP蛋白基因家族；表達分析

中圖分類號：S188；S529.01" 文獻標志碼：A

文章編號：1002-1302（2025）04-0209-09

收稿日期：2024-04-26

基金項目：福建省科技計劃-省屬公益類科研院所基本科研專項（編號：2024R10300010）；福建省農(nóng)業(yè)科學(xué)院對外合作項目（編號：DWHZ-2024-18）；福建省農(nóng)業(yè)科學(xué)院科技創(chuàng)新團隊（編號：CXTD2021011-2）；福建省人民政府、中國農(nóng)業(yè)科學(xué)院“5511”協(xié)同創(chuàng)新工程（編號：XTCXGC2021019）。

作者簡介：李程勛（1991—），男，福建南安人，碩士，助理研究員，主要從事作物遺傳育種及農(nóng)產(chǎn)品天然產(chǎn)物提取研究。E-mail：1219513539@qq.com。

通信作者：李愛萍，研究員，主要從事作物遺傳育種及農(nóng)產(chǎn)品天然產(chǎn)物提取研究。E-mail：apl909@163.com。

木豆（Cajanus cajan），豆科木豆屬，是一種多年生木本作物，別稱鴿豆、樹豆、柳豆等。木豆是世界第六大食用豆類，也是唯一的木本食用豆，木豆含有較高的蛋白質(zhì)，被稱為“窮人的肉食”，是半干旱地區(qū)重要的食用豆類作物［1-2］。木豆還含有豐富的淀粉、氨基酸、纖維素、脂肪、維生素等［3-5］，具有較高的營養(yǎng)價值。木豆葉具有清熱解毒、補腎健骨、止痛止血、降低血液中膽固醇含量等功效［6］，能夠治療骨頭壞死、心血管疾病、肺病、口腔潰瘍等［7-8］。木豆具有良好的耐旱、耐瘠等抗逆性［9］，能適應(yīng)各種類型的土壤，包括沙土、鹽漬土、壤土等［10］，研究木豆生長發(fā)育中抗逆響應(yīng)因子，可為木豆優(yōu)良品種的選育提供參考依據(jù)。熱脅迫會影響植物生物量的累積，造成產(chǎn)量下降，這種現(xiàn)象在熱帶、亞熱帶地區(qū)尤為明顯。當前溫室效應(yīng)越來越嚴重，氣溫升高速度變快，前人研究表明，氣溫每升高3～4 ℃作物產(chǎn)量會下降15%～35%［11］。栽培木豆最佳溫度為25～35 ℃，在木豆生長的季節(jié)中，最高溫度每上升 1 ℃，會造成木豆減產(chǎn)20.8%［12］。因此研究木豆高溫脅迫響應(yīng)因子，為木豆耐熱品種選育打下基礎(chǔ)。

OFP（OVATE family proteins）卵形家族蛋白是植物所特有的一類轉(zhuǎn)錄因子，在番茄中首次被報道［13］，其C端有1段約70個氨基酸的保守結(jié)構(gòu)域，該結(jié)構(gòu)域被命名為OVATE結(jié)構(gòu)域。

目前OFP蛋白的研究在擬南芥（Arabidopsis thaliana）、水稻（Oryza sativa）、小麥（Triticum aestivum）、玉米（Zea mays）和番茄（Solanum lycopersicum）等中均有報道［14-17］，OFP蛋白具有調(diào)控植株生長發(fā)育的多種功能，比如能調(diào)控果實形態(tài)［13］、次生細胞壁形成［18］、胚囊和維管束發(fā)育等［19］。擬南芥中AtOPF1通過抑制赤霉素生物合成途徑中關(guān)鍵酶基因AtGA20ox1的轉(zhuǎn)錄活性，抑制細胞伸長，造成植株矮?。?0］，同時AtOPF1與BLH3（bel1-like homeodomain 3）互作調(diào)控擬南芥從營養(yǎng)生長階段轉(zhuǎn)變到生殖生長階段，與ATH1（homeobox gene 1）互作調(diào)控擬南芥莖的生長、開花時間和花基邊界形成［21］，與KNAT5（knotted1-like homeobox gene 5）互作調(diào)控擬南芥?zhèn)雀鶖?shù)目和根長［22］；AtOPF5蛋白和BLH1-KNAT3蛋白復(fù)合物互作負調(diào)控擬南芥早期胚囊的發(fā)育［23］；AtOFP4與KNAT7互作調(diào)控次生壁的形成［24］；AtOFP15和AtOFP16可能參與調(diào)控擬南芥莢果形態(tài)［25-26］。

OFP蛋白在植物非生物逆境脅迫中也具有很重要的作用，擬南芥中AtOFP8對干旱和高溫敏感，過量表達AtOFP8增強了擬南芥對干旱的耐受性；水稻OsOFP6正向調(diào)控水稻耐冷性［14］；油桐（Vernicia fordii）VfOFP7在低溫下被明顯誘導(dǎo)［27］；在小麥冷馴化研究中得出需要多個TaOFP基因來維持細胞穩(wěn)態(tài)，表明TaOFP在小麥抗寒性的獲得中起到重要作用［28］；多個蘋果（Malus pumila Mill.）OFP家族基因在低溫中表達上調(diào)，OFP家族基因在高低溫、干旱和鹽脅迫中均有一定的響應(yīng)［29］。

OFP基因家族在植物生長發(fā)育和非生物逆境脅迫中起著重要的調(diào)控作用，尤其在高低溫、干旱等非生物脅迫中發(fā)揮著重要的作用，目前有關(guān)木豆OFP基因家族的相關(guān)研究未見報道。本研究通過生物信息學(xué)方法鑒定木豆全基因組中的OFP成員（稱為CcOFP），進一步分析其蛋白理化性質(zhì)、基因結(jié)構(gòu)、保守基序、家族進化關(guān)系、蛋白質(zhì)二級結(jié)構(gòu)、蛋白質(zhì)三級結(jié)構(gòu)和蛋白互作網(wǎng)絡(luò)，并利用轉(zhuǎn)錄組數(shù)據(jù)分析基因組織表達特性和高溫脅迫下的表達模式，為后續(xù)進一步解析CcOFP基因生長發(fā)育調(diào)控機制和耐逆功能提供理論參考依據(jù)。

1 材料與方法

1.1 試驗時間和地點

本試驗于2023年5月在福建省農(nóng)業(yè)科學(xué)院科研綜合實驗樓作物研究所實驗室內(nèi)完成。

1.2 CcOFP基因家族成員鑒定

從NCBI數(shù)據(jù)庫下載木豆基因組和蛋白質(zhì)組數(shù)據(jù)。利用數(shù)據(jù)庫下載OFP的隱馬爾可夫模型，利用hmmsearch得到候選基因家族成員，利用已發(fā)表的擬南芥OFP蛋白基因家族成員的氨基酸序列為探針，進行BLASTp比對獲得同源蛋白序列，結(jié)合hmmsearch和BLASTp結(jié)果，利用NCBI-CDD數(shù)據(jù)庫進一步確認鑒定的OFP蛋白序列中的OFP結(jié)構(gòu)域。利用ExPASy數(shù)據(jù)庫分析蛋白質(zhì)的理化特性。

1.3 系統(tǒng)進化樹比對

提取擬南芥和木豆OFP蛋白氨基酸序列，利用MEGA 7.0軟件構(gòu)建木豆和擬南芥OFP蛋白系統(tǒng)進化樹，采用MUSCLE方法進行多序列比對，采用Neighbor-Joining方法構(gòu)建進化樹，Bootstrap值設(shè)定為1 000，其他參數(shù)采用默認值，利用iToL軟件對進化樹進行美化。

1.4 CcOFP基因結(jié)構(gòu)與氨基酸保守基序分析

利用MEME數(shù)據(jù)庫分析木豆OFP蛋白的氨基酸保守基序，motif值設(shè)為11，其他參數(shù)采用默認值，利用TBtools軟件繪制基因結(jié)構(gòu)和保守基序圖。

1.5 CcOFP基因啟動子順式作用元件分析

利用Plant CARE對木豆OFP基因轉(zhuǎn)錄起始位點上游2 000 bp序列的順勢作用元件進行預(yù)測，利用TBtools軟件繪制啟動子順式作用元件結(jié)構(gòu)圖。

1.6 木豆OFP蛋白質(zhì)二級結(jié)構(gòu)和三級結(jié)構(gòu)分析

利用SOPMA和SWISS-MODEL數(shù)據(jù)庫分別對木豆OFP蛋白二級結(jié)構(gòu)和三級結(jié)構(gòu)進行預(yù)測。

1.7 木豆OFP蛋白互作網(wǎng)絡(luò)分析及可視化

利用String數(shù)據(jù)庫以擬南芥OFP蛋白模型構(gòu)建木豆OFP蛋白互作網(wǎng)絡(luò)［30］，并利用Cytoscape軟件繪圖。

1.8 CcOFP基因組織表達差異和高溫脅迫響應(yīng)分析

通過NCBI數(shù)據(jù)庫下載木豆根、莖、頂芽、葉及42 ℃脅迫下的RNA-Seq數(shù)據(jù)，登錄號見表1和表2，利用TBtools中的Trimmomatic對原始數(shù)據(jù)進行質(zhì)量控制，通過Kallisto分析得到OFP基因的表達量矩陣，利用TPM值進行熱圖繪制［31］。

2 結(jié)果與分析

2.1 木豆OFP蛋白基因家族成員鑒定及其理化特性分析

由表3可知，共鑒定到木豆OFP蛋白基因21個，依次命名為CcOFP1～CcOFP21，CcOFP蛋白序列長度為155～403個氨基酸（aa），根據(jù)數(shù)據(jù)庫中基因位置信息得出，2、3、7、11號染色體上各分布有2個CcOFP基因；1、4、8號染色體上各有1個CcOFP基因，還有10個CcOFP基因未定位。

蛋白質(zhì)的理化性質(zhì)與蛋白質(zhì)功能息息相關(guān)，利用ExPASy數(shù)據(jù)庫的ProParam工具對CcOFP蛋白質(zhì)的理化性質(zhì)進行分析，結(jié)果見表2。CcOFP蛋白等電點在4.55～9.87之間，其中偏堿性（等電點＞7）的CcOFP蛋白有11個，偏酸性（等電點＜7）的有10個；不穩(wěn)定系數(shù)為41.01～75.24，均屬于不穩(wěn)定蛋白；脂溶指數(shù)為52.59～83.45；總平均親水性均＜0，均為親水性蛋白。

2.2 OFP蛋白基因家族的進化分析

將木豆OFP基因家族成員（21個）和擬南芥OFP基因家族成員（18個）的蛋白質(zhì)序列進行系統(tǒng)發(fā)育樹構(gòu)建，結(jié)果見圖1。結(jié)果顯示，39個OFP蛋白成員被分為3個類群，分別為GroupⅠ、GroupⅡ和GroupⅢ。GroupⅠ由2個擬南芥OFP蛋白和1個木豆OFP蛋白組成，GroupⅡ由7個擬南芥OFP蛋白和9個木豆OFP蛋白組成，GroupⅢ由9個擬南芥OFP蛋白和11個木豆OFP蛋白組成。在進化樹中CcOFP16蛋白和AtOFP14蛋白的親緣關(guān)系較近，CcOFP7蛋白和AtOFP13蛋白的親緣關(guān)系較近，CcOFP12蛋白和AtOFP8蛋白的親緣關(guān)系較近，CcOFP17蛋白和AtOFP5蛋白的親緣關(guān)系較近，推測他們可能具有相似的基因功能，在OFP進化樹的分支中木豆和擬南芥二者相互交錯聚類在一起，說明OFP的進化速度較快，且在進化過程中產(chǎn)生了較大的分化。

2.3 CcOFP的基因結(jié)構(gòu)域Motif分析

基于CcOFP基因組和編碼區(qū)序列繪制基因結(jié)構(gòu)圖，結(jié)果見圖2；利用MEME軟件對木豆OFP蛋白保守基序進行預(yù)測，繪制Motif結(jié)構(gòu)圖，結(jié)果見圖3。在CcOFP基因家族中，CcOFP13、CcOFP18各有1個內(nèi)含子，CcOFP17沒有3′和5′非編碼區(qū)，CcOFP1、CcOFP8沒有3′非編碼區(qū)。

保守基序的分析有利于探索蛋白質(zhì)的結(jié)構(gòu)組成，分析結(jié)果表明，木豆OFP蛋白家族具有3～9個保守基序（motif），全部成員均含有motif1、motif2、motif3保守基序，且先后順序都一致，說明OFP蛋白

家族在木豆中比較保守。

2.4 CcOFP基因啟動子順式響應(yīng)元件分析

通過對CcOFP基因啟動子的順式響應(yīng)元件進行分析，得到圖4。CcOFP基因家族啟動子序列含有17類順式響應(yīng)元件，包括4類脅迫響應(yīng)元件：低溫響應(yīng)、厭氧脅迫響應(yīng)、防御和壓力響應(yīng)、干旱響應(yīng)；5類激素響應(yīng)元件：赤霉素響應(yīng)、茉莉酸甲酯響應(yīng)、脫落酸響應(yīng)、水楊酸響應(yīng)、生長素響應(yīng)；8類生長相關(guān)順式作用元件：光響應(yīng)、胚乳表達調(diào)控、晝夜節(jié)律調(diào)控、分生組織表達、種子特異性調(diào)控、柵欄葉肉細胞分化調(diào)控、類黃酮生物合成基因調(diào)控、根特異性元件。其中豐度最大的為光響應(yīng)元件，共有214個響應(yīng)元件，厭氧脅迫響應(yīng)元件和茉莉酸甲酯響應(yīng)元件豐度也較大，分別有31個和30個。表明CcOFP基因在逆境脅迫響應(yīng)、激素調(diào)節(jié)和生長發(fā)育方面發(fā)揮重要作用。

2.5 木豆OFP蛋白的二級結(jié)構(gòu)和三級結(jié)構(gòu)分析

對木豆OFP蛋白的二級結(jié)構(gòu)進行預(yù)測分析，結(jié)果見表4。21個木豆OFP蛋白質(zhì)的二級結(jié)構(gòu)均含有α-螺旋、β-轉(zhuǎn)角、延伸鏈和無規(guī)則卷曲結(jié)構(gòu)，其中α-螺旋結(jié)構(gòu)數(shù)量比例為20.37%～56.77%，β-轉(zhuǎn)角結(jié)構(gòu)數(shù)量比例為0.80%～9.03%， 延伸鏈

結(jié)構(gòu)數(shù)量比例為6.72%～20.44%，無規(guī)則卷曲結(jié)構(gòu)數(shù)量比例為26.45%～69.00%。可見木豆OFP蛋白的二級結(jié)構(gòu)是以α-螺旋和無規(guī)則卷曲結(jié)構(gòu)為主要組成成分。

對木豆OFP蛋白的三級結(jié)構(gòu)進行預(yù)測分析，結(jié)果見圖5。由圖5可見，木豆OFP蛋白的空間結(jié)構(gòu)主要由α-螺旋和無規(guī)則卷曲結(jié)構(gòu)組成，β-轉(zhuǎn)角和延伸鏈結(jié)構(gòu)較少。

2.6 木豆OFP蛋白互作網(wǎng)絡(luò)分析

利用string數(shù)據(jù)庫以擬南芥OFP蛋白模型構(gòu)建木豆OFP蛋白互作網(wǎng)絡(luò)（圖6）。結(jié)果顯示CcOFP6、CcOFP8、CcOFP9、CcOFP10、CcOFP14、CcOFP16、CcOFP17、CcOFP18、CcOFP19、CcOFP20、CcOFP21共11個蛋白參與互作，存在23條互作關(guān)系。

2.7 OFP基因的組織表達和高溫脅迫表達分析

通過對木豆根（S1）、莖（S2）、頂芽（S3）、葉（S4）的RNA-Seq數(shù)據(jù)進行分析，得出OFP基因在不同組織部位的表達水平（圖7）。由圖7可以看出，OFP基因在頂芽中的表達量較大，在根和葉中的表達量較小，其中CcOFP7、CcOFP12和CcOFP14在頂芽中具有較高的表達量，CcOFP15在莖中也有較高的表達量，CcOFP8、CcOFP13、CcOFP21在各個組織中表達量均較低或未表達。

通過木豆42 ℃高溫脅迫下的RNA-Seq數(shù)據(jù)分析，得出OFP基因在高溫脅迫下的表達水平（圖8）。由圖8可見，在高溫脅迫下大部分OFP均有響應(yīng)，其中CcOFP2、CcOFP3、CcOFP6、CcOFP12、CcOFP13、CcOFP18、CcOFP19、CcOFP20、CcOFP21的表達量較小，CcOFP1、CcOFP8具有較為明顯的上調(diào)表達，CcOFP7、CcOFP9、CcOFP10、CcOFP11、CcOFP14具有較為明顯的下調(diào)表達。

3 討論

本研究利用生物信息學(xué)方法鑒定出木豆OFP蛋白家族成員21個，OFP家族的成員數(shù)量在不同物種間存在差異，呈現(xiàn)出單子葉植物的OFP家族成員多于雙子葉植物的趨勢，比如單子葉植物玉米有44個［32］，谷子（Setaria italica）有34個［33］，水稻有31個［34］，雙子葉植物擬南芥有18個［16］，葡萄（Vitis vinifera）有17個［35］， 油桐有12個［27］，推測在植物

進化過程中OFP家族發(fā)生了收縮。進一步通過蛋白質(zhì)理化分析和保守基序分析發(fā)現(xiàn)木豆OFP蛋白均為不穩(wěn)定的親水性蛋白，均含有motif1、motif2、motif3保守基序，蛋白二級和三級結(jié)構(gòu)均以α-螺旋和無規(guī)則卷曲結(jié)構(gòu)為主要組成成分，表明木豆OFP蛋白家族成員間存在保守性，這與程欣在番茄OFP基因家族中的研究結(jié)果［36］一致。

植物逆境響應(yīng)方面，AtOFP8對干旱和高溫敏感，過量表達AtOFP8增強了擬南芥對干旱的耐受性［37］，CcOFP12與AtOFP8結(jié)構(gòu)相近，同時在CcOFP5、CcOFP8、CcOFP11、 CcOFP12、 CcOFP17基

因啟動子上發(fā)現(xiàn)干旱響應(yīng)元件，因此推測CcOFP12基因可能在木豆抗旱過程中起到重要作用。在CcOFP1、CcOFP9、CcOFP11、CcOFP15、CcOFP16和CcOFP18基因啟動子上發(fā)現(xiàn)低溫響應(yīng)元件，推測這些CcOFP基因可能在木豆高低溫脅迫過程中起到重要作用，結(jié)合高溫脅迫表達分析，得出CcOFP1基因啟動子響應(yīng)高溫脅迫信號，基因表達量上調(diào)；CcOFP9和CcOFP11基因啟動子響應(yīng)高溫脅迫信號，基因表達量下調(diào)，其他CcOFP高溫響應(yīng)基因可能通過蛋白互作來響應(yīng)高溫脅迫信號。

生長發(fā)育響應(yīng)方面，本研究OFP進化樹的分支中木豆和擬南芥二者相互交錯聚類在一起，兩者有多個基因具有較近的親緣關(guān)系，因此可能具有相似的功能，比如擬南芥AtOPF5蛋白負調(diào)控擬南芥早期胚囊的發(fā)育［23］，木豆CcOFP17蛋白與AtOPF5蛋白具有較近的親緣關(guān)系，推測CcOFP17蛋白對木豆早期胚囊的發(fā)育起到調(diào)控作用。同時CcOFP基因啟動子也具有光響應(yīng)、胚乳表達調(diào)控、晝夜節(jié)律調(diào)控等多種生長相關(guān)順式作用元件，可以進一步對其進行相應(yīng)的功能特性研究。在組織特異性表達方面，CcOFP7、CcOFP12和CcOFP14在頂芽中具有較高的表達量，CcOFP15在莖中有較高的表達量，表明木豆OFP基因可能主要參與芽的生長調(diào)控。

4 結(jié)論

綜上，本研究利用生物信息學(xué)方法對木豆OFP蛋白基因家族進行了較為全面的分析，鑒定出21個木豆OFP蛋白家族成員，具有較高的保守性，均屬于不穩(wěn)定親水性蛋白，空間結(jié)構(gòu)主要由α-螺旋和無規(guī)則卷曲結(jié)構(gòu)組成。啟動子順式響應(yīng)元件分析結(jié)合高溫脅迫表達分析篩選出高溫響應(yīng)候選基因CcOFP1、CcOFP9和CcOFP11，組織表達分析篩選出頂芽發(fā)育響應(yīng)候選基因CcOFP7、CcOFP12和CcOFP14，莖發(fā)育候選基因CcOFP15。

參考文獻：

［1］Ghosh G，Ganguly S，Purohit A，et al. Transgenic pigeonpea events expressing Cry1Ac and Cry2Aa exhibit resistance to Helicoverpa armigera［J］. Plant Cell Reports，2017，36（7）：1037-1051.

［2］Saxena R K，Obala J，Sinjushin A，et al. Characterization and mapping of Dt1 locus which co-segregates with CcTFL1 for growth habit in pigeon pea［J］. Theoretical and Applied Genetics，2017，130（9）：1773-1784.

［3］Garg N，Singh S. Arbuscular mycorrhiza Rhizophagus irregularis and silicon modulate growth，proline biosynthesis and yield in Cajanus cajan L. millsp. （pigeon pea） genotypes under cadmium and zinc stress［J］. Journal of Plant Growth Regulation，2018，37（1）：46-63.

［4］Ayenan M A T，Ofori K，Ahoton L E，et al. Pigeonpea ［Cajanus cajan （L.） Millsp.］production system，farmers’ preferred traits and implications for variety development and introduction in Benin［J］. Agriculture amp; Food Security，2017，6（1）：48.

［5］Jacob C，Carrasco B，Schwember A R. Advances in breeding and biotechnology of legume crops［J］. Plant Cell，Tissue and Organ Culture，2016，127（3）：561-584.

［6］Dong B Y，Niu L L，Meng D，et al. Genome-wide analysis of MATE transporters and response to metal stress in Cajanus cajan［J］. Journal of Plant Interactions，2019，14（1）：265-275.

［7］Tang R，Tian R H，Cai J Z，et al. Acute and sub-chronic toxicity of Cajanus cajan leaf extracts［J］. Pharmaceutical Biology，2017，55（1）：1740-1746.

［8］Dinore J M，F(xiàn)arooqui M. GC-MS and LC-MS：an integrated approach towards the phytochemical evaluation of methanolic extract of pigeon pea［Cajanus cajan （L.） Millsp.］leaves［J］. Natural Product Research，2022，36（8）：2177-2181.

［9］Nagaraj K，Vanishree S，Muthukumar T. Genotypic variation in response and dependency of Cajanus cajan （L.） Millsp.，on arbuscular mycorrhizal fungi in a tropical Alfisol［J］. Plant Biosystems-an International Journal Dealing with All Aspects of Plant Biology，2021，155（4）：878-890.

［10］李程勛，李愛萍，徐曉俞，等. 淺談木豆的抗逆機制及在福建的應(yīng)用前景［J］. 作物雜志，2018（4）：28-31.

［11］Ramakrishna G，Kaur P，Singh A，et al. Comparative transcriptome analyses revealed different heat stress responses in pigeon pea （Cajanus cajan） and its crop wild relatives［J］. Plant Cell Reports，2021，40（5）：881-898.

［12］Mishra S，Singh R，Kumar R，et al. Impact of climate change on pigeon pea［J］. Economic Affairs，2017，62（3）：455-457.

［13］Liu J P，Van Eck J，Cong B，et al. A new class of regulatory genes underlying the cause of pear-shaped tomato fruit［J］. Proceedings of the National Academy of Sciences of the United States of America，2002，99（20）：13302-13306.

［14］Ma Y M，Yang C，He Y，et al. Rice OVATE family protein 6 regulates plant development and confers resistance to drought and cold stresses［J］. Journal of Experimental Botany，2017，68（17）：4885-4898.

［15］Tang Y，Zhang W，Yin Y L，et al. Expression of ovate family protein 8 affects epicuticular waxes accumulation in Arabidopsis thaliana［J］. Botanical Studies，2018，59（1）：12.

［16］Wang S C，Chang Y，Guo J J，et al. Arabidopsis OVATE family proteins，a novel transcriptional repressor family，control multiple aspects of plant growth and development［J］. PLoS One，2011，6（8）：e23896.

［17］Liu D，Sun W，Yuan Y W，et al. Phylogenetic analyses provide the first insights into the evolution of OVATE family proteins in land plants［J］. Annals of Botany，2014，113（7）：1219-1233.

［18］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 amp; Behavior，2015，10（7）：e1033126.

［19］Schmitz A J，Begcy K，Sarath G，et al. Rice OVATE family protein 2 （OFP2） alters hormonal homeostasis and vasculature development［J］. Plant Science，2015，241：177-188.

［20］Wang S C，Chang Y，Guo J J，et al. Arabidopsis OVATE family protein 1 is a transcriptional repressor that suppresses cell elongation［J］. Plant Journal，2007，50（5）：858-872.

［21］Zhang L G，Zhang X F，Ju H X，et al. OVATE family protein1 interaction with BLH3 regulates transition timing from vegetative to reproductive phase in Arabidopsis［J］. Biochemical and Biophysical Research Communications，2016，470（3）：492-497.

［22］尹艷莉，唐 堯，馮 鵬，等. 擬南芥卵形家族蛋白AtOFP1與KNAT5互作及功能初步分析［J］. 農(nóng)業(yè)生物技術(shù)學(xué)報，2018，26（11）：1865-1871.

［23］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.

［24］Li E Y，Wang S C，Liu Y Y，et al. OVATE family protein4 （OFP4） interaction with KNAT7 regulates secondary cell wall formation in Arabidopsis thaliana［J］. Plant Journal，2011，67（2）：328-341.

［25］蔡 玲. 卵形蛋白家族轉(zhuǎn)錄因子AtOFP15對擬南芥莢果形態(tài)的調(diào)控［D］. 長春：東北師范大學(xué)，2017：32-36.

［26］楊 麗. Ⅲ型卵形蛋白家族轉(zhuǎn)錄因子AtOFP16調(diào)控擬南芥莢果形態(tài)的機制［D］. 長春：東北師范大學(xué)，2017：31-37.

［27］王靜雅，劉文娟，呂祎馨，等. 油桐OFP基因家族的全基因組鑒定及低溫響應(yīng)分析［J］. 植物生理學(xué)報，2020，56（5）：949-960.

［28］徐 渴，王偉偉，武寧靜，等. 小麥OFP基因家族鑒定及其低溫脅迫的表達分析［J］. 農(nóng)業(yè)生物技術(shù)學(xué)報，2021，29（9）：1665-1677.

［29］許瑞瑞，李 睿，王小非，等. 蘋果OFP基因家族的全基因組鑒定與非生物逆境表達分析［J］. 中國農(nóng)業(yè)科學(xué)，2018，51（10）：1948-1959.

［30］楚宗麗，李亮杰，姬 虹，等. 小麥Hsp70基因家族鑒定及蛋白互作網(wǎng)絡(luò)分析［J］. 江蘇農(nóng)業(yè)科學(xué)，2022，50（10）：37-44.

［31］張中起，王 嬌，靳 煒，等. 陸地棉CRK基因家族的鑒定及其表達分析［J］. 中國農(nóng)業(yè)科學(xué)，2018，51（13）：2442-2461.

［32］張彥琴，常建忠，董春林，等. 玉米OFP基因家族全基因組鑒定及表達分析［J］. 分子植物育種，2020，18（20）：6595-6604.

［33］郭 昱，李 娜，李 珊，等. 谷子OFP基因家族的鑒定及表達分析［J］. 山西農(nóng)業(yè)科學(xué)，2020，48（4）：467-476.

［34］Yu H，Jiang W Z，Liu Q，et al. Expression pattern and subcellular localization of the ovate protein family in rice［J］. PLoS One，2015，10（3）：e0118966.

［35］袁 月，張亞光，高世敏，等. 葡萄OVATE基因家族生物信息學(xué)及表達［J］. 中國農(nóng)業(yè)科學(xué)，2016，49（19）：3786-3797.

［36］程 欣. 番茄OFP基因家族生物信息學(xué)分析［J］. 江蘇農(nóng)業(yè)科學(xué)，2021，49（18）：39-48.

［37］劉守銀，唐 堯，韓佳良，等. 卵型家族蛋白轉(zhuǎn)錄因子AtOFP8的抗旱性功能分析［J］. 華北農(nóng)學(xué)報，2020，35（5）：98-106.