杧果RAV基因家族的全基因組分析

孫宇 劉志鑫 葉子 羅睿雄 劉曉妹 蒲金基 張賀

摘要： 為了揭示RAV家族基因在杧果中的序列特征及其表達特性，采用生物信息學(xué)方法對杧果RAV家族基因進行序列分析，并通過qRT-PCR技術(shù)研究杧果膠孢炭疽菌和細菌性黑斑病菌侵染過程中該家族基因的相對表達量。結(jié)果表明，從杧果基因組中鑒定出6個RAV基因家族成員，其編碼的蛋白質(zhì)均具有AP2和B3超家族保守域結(jié)構(gòu)，命名為MiRAV1～MiRAV6。MiRAV1和MiRAV5為不穩(wěn)定親水堿性蛋白質(zhì)，MiRAV2和MiRAV3為不穩(wěn)定親水酸性蛋白質(zhì)，MiRAV4為穩(wěn)定親水酸性蛋白質(zhì)，MiRAV6為穩(wěn)定親水堿性蛋白質(zhì)。杧果RAV基因編碼的蛋白質(zhì)主要定位于細胞核中，無規(guī)則卷曲和α-螺旋是6個MiRAVs蛋白二級結(jié)構(gòu)的主要元件。系統(tǒng)進化樹結(jié)合基序分析結(jié)果表明，在杧果與擬南芥、煙草、蘋果、菠蘿和毛果楊中，RAV家族基因具有多樣性，在進化上結(jié)構(gòu)具有保守性。qRT-PCR結(jié)果顯示，膠孢炭疽菌侵染過程中，杧果RAV基因的相對表達量均顯著下調(diào);細菌性黑斑病菌侵染過程中，MiRAV1～ MiRAV6的相對表達量在3 h時均顯著下調(diào)，MiRAV1、MiRAV2、MiRAV3和MiRAV6的相對表達量在6 h時均顯著上調(diào)，MiRAV1、MiRAV5和MiRAV6的相對表達量分別在12 h、24 h、48 h、72 h時均顯著上調(diào)。以上發(fā)現(xiàn)將為研究杧果RAV基因家族成員的功能和作用機制奠定基礎(chǔ)。

關(guān)鍵詞： 杧果;RAV基因家族;全基因組鑒定;表達分析

中圖分類號： S667.7?? 文獻標識碼： A?? 文章編號： 1000-4440（2021）04-0957-11

Genome-wide analysis of the RAV gene family in mango

SUN Yu1，2， LIU Zhi-xin2， YE Zi1，2， LUO Rui-xiong3， LIU Xiao-mei1，2， PU Jin-ji1，2， ZHANG He1，2

（1.College of Plant Protection， Hainan University， Haikou 570228， China;2.Environment and Plant Protection Institute， Chinese Academy of Tropical Agricultural Sciences/Key Laboratory of Integrated Pest Management on Tropical Grops， Ministry of Agriculture and Rural Affairs， Haikou 571101， China;3.Tropical Crops Genetic Resources Institute， Chinese Academy of Tropical Agricultural Sciences， Haikou 571101， China）

Abstract： In order to reveal the sequence and expression characteristics of RAV gene in Mangifera indica， the sequence of the RAV gene family was analyzed using bioinformatics methods， and the relative expression level of RAV gene family was studied by qRT-PCR during the infection of Colletotrichum gloeosporioides and Xanthomonas campestris pv. mangiferaeindicae. The results showed that six RAV gene family members were identified from the mango genome， and the encoded proteins had the conserved domains of the AP2 and B3 superfamily， named MiRAV1-MiRAV6. MiRAV1 and MiRAV5 were unstable hydrophilic basic proteins， MiRAV2 and MiRAV3 were unstable hydrophilic acidic proteins， MiRAV4 was a stable hydrophilic acidic protein， and MiRAV6 was a stable hydrophilic basic protein. The proteins encoded by RAV gene of mango were mainly located in the nucleus， and random coil and α-helix were the main elements of secondary structure of six MiRAVs proteins. Phylogenetic tree construction combined with motif analysis showed that， among mango and Arabidopsis thaliana， tobacco， apple， pineapple and Populus tomentosa， the RAV family genes were diverse and their structure was conservative in evolution. The results of qRT-PCR indicated that the relative expression of RAV gene in mango was significantly down regulated during the infection of Colletotrichum gloeosporioides. During the infection of Xanthomonas campestris pv. mangiferaeindicae， the relative expression levels of MiRAV1-MiRAV6 were significantly down regulated at 3 h， and the relative expression levels of MiRAV1， MiRAV2， MiRAV3 and MiRAV6 were significantly up-regulated at 6 h， and the relative expression levels of MiRAV1， MiRAV5 and MiRAV6 were significantly up-regulated at 12 h， 24 h， 48 h and 72 h， respectively. These results will lay a foundation for further study on the function and mechanism of RAV gene family members in mango.

Key words： mango;RAV gene family;genomic identification;expression analysis

轉(zhuǎn)錄因子，又稱反式作用因子，其具有特殊的結(jié)構(gòu)可調(diào)控基因表達，是植物中重要的調(diào)控因子。植物轉(zhuǎn)錄因子在脅迫反應(yīng)中，對調(diào)控基因表達發(fā)揮重要作用[1-3]，當植物受到逆境脅迫時，轉(zhuǎn)錄因子與相應(yīng)的順式作用元件結(jié)合，可以調(diào)控下游相關(guān)基因的表達 [4]。RAV（Related to AB13/VP1）轉(zhuǎn)錄因子廣泛存在于植物中，該家族同時具有2個保守結(jié)構(gòu)域[5-8]：位于N端，可與GCC-box區(qū)域特異性結(jié)合的AP2結(jié)構(gòu)域;位于C端，可與CACCTG序列特異性結(jié)合的B3結(jié)構(gòu)域[9]。VP1/ABI3以及ARFs轉(zhuǎn)錄因子家族都具有該結(jié)構(gòu)域，主要調(diào)控植物對生長素和脫落酸的反應(yīng)[10]。RAV屬于AP2/ERF家族的一個亞家族，該家族還包含ERF、AP2、DREB亞家族[11]。

迄今為止，RAV基因家族成員已在擬南芥（Arabidopsis thaliana）[12-13]、黃瓜（Cucumis sativus）[14]、棉花（Gossypium）[15]、東方山羊豆（Galega orientalis）[16]、煙草（Nicotiana tabacum）[17]、茶樹（Camellia sinensis）[18]等植物中被陸續(xù)克隆和鑒定。例如：在擬南芥中發(fā)現(xiàn)7個RAV基因家族成員AtRAV1～AtRAV7，AtRAV1和AtRAV2在擬南芥中參與應(yīng)對外界脅迫和葉片衰老 [12-13];黃瓜中CsRAV在各個組織中均有表達[14];棉花中GhRAV表達受黃萎病菌的誘導(dǎo)[15];東方山羊豆中GoRAV表達不僅受低溫和聚乙二醇的誘導(dǎo)，還受茉莉酸甲酯、水楊酸和脫落酸的誘導(dǎo)[16];煙草中GmRAV過量表達會導(dǎo)致葉片失綠，同時抑制植株的生長 [17];茶樹中CsRAV表達受NaCl、乙烯和低溫的誘導(dǎo)[18]。RAV基因廣泛的分布表明它們在植物的生長發(fā)育和抗逆反應(yīng)調(diào)控過程中起著非常重要的作用。

杧果（Mangifera indica），被稱為熱帶水果之王，是熱帶地區(qū)重要的果樹，也是農(nóng)民脫貧致富的“金杧果”和“致富樹”。隨著2020年杧果全基因組測序的完成 [19]，杧果基因功能的研究進程明顯加快，然而目前國內(nèi)外對杧果RAV基因家族的研究還未見報道。為此，本研究對杧果RAV基因家族成員的序列特征及其表達特性進行分析，為杧果基因功能的分析以及杧果的遺傳育種等方面提供參考。

1 材料與方法

1.1 杧果生物脅迫處理

杧果嫁接苗由農(nóng)業(yè)農(nóng)村部儋州杧果種質(zhì)資源圃提供，品種為貴妃杧。將株高約60 cm，均勻一致，樹齡1年的杧果嫁接苗移栽到花盆內(nèi)，以室溫25 ℃，相對濕度70%～90%，光周期12 h（光）/12 h（暗）的條件下進行生物脅迫處理。設(shè)置膠孢炭疽菌和細菌性黑斑病菌 2個處理，每個處理3株苗，處理0 h作為對照。以分生孢子含量為1 ml 2×106 個的膠孢炭疽菌分生孢子懸浮液、含量為2×107CFU 的細菌性黑斑病菌懸浮液分別對杧果嫁接苗葉片均勻噴霧，取樣時間分別為0 h、3 h、6 h、12 h、24 h、48 h、72 h，將杧果葉片剪碎，液氮速凍，-80 ℃保存待用。

1.2 杧果全基因組數(shù)據(jù)的來源

杧果全基因組數(shù)據(jù)來源于NCBI（PRJNA487154），擬南芥、煙草、蘋果（Malus domestica）、菠蘿（Ananas comosus）、毛果楊（Populus trichocarpa）5個物種的基因組和RAV基因信息分別來源于NCBI（https：//www.ncbi.nlm.nih.gov/）和PlantTFDB（http：//planttfdb.gao-lab.org/）數(shù)據(jù)庫。

1.3 杧果RAV基因家族成員的鑒定

從擬南芥數(shù)據(jù)庫獲取AtRAV1～AtRAV7氨基酸序列，并將其作為查詢對象在杧果基因組數(shù)據(jù)庫中執(zhí)行Blastp比較（E≤10-10）。使用Pfam工具（http：//pfam.xfam.org/）檢測得到的杧果蛋白質(zhì)氨基酸序列結(jié)構(gòu)域，進一步去除沒有典型AP2和B3結(jié)構(gòu)域的蛋白質(zhì)氨基酸序列，最終得到所有杧果RAV基因家族成員[20]。其他5個物種中的RAV基因家族成員采用類似的方法進行篩選。

1.4 杧果RAV基因家族的生物信息學(xué)分析

通過在線工具ProtParam預(yù)測RAV相對分子質(zhì)量、等電點、外顯子數(shù)、脂溶指數(shù)、氨基酸數(shù)等理化性質(zhì)[21];利用PSORT在線軟件進行亞細胞定位預(yù)測;通過SOPMA在線軟件預(yù)測蛋白質(zhì)的二級結(jié)構(gòu);利用NCBI Conserved Domain Search預(yù)測RAV的保守結(jié)構(gòu)域，再通過SWISS-MODEL在線工具對蛋白質(zhì)保守域的三級結(jié)構(gòu)進行預(yù)測;利用MEME在線軟件對杧果RAV家族的蛋白質(zhì)保守基序進行分析[22]，再通過DNAMAN 6.0軟件進行高同源蛋白的氨基酸序列比對。

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

利用ClustalW程序?qū)x果和擬南芥、煙草、蘋果、菠蘿、毛果楊的RAV家族成員的氨基酸序列進行多序列比對。獲得RAV蛋白氨基酸序列后使用MEGA 7.0軟件，通過鄰近法構(gòu)建進化樹，其中校驗參數(shù)Bootstrap值設(shè)置為1 000次重復(fù)[23]。

1.6 RNA提取及cDNA合成

杧果葉片總RNA提取參照RNAprep Pure多糖多酚植物總RNA提取試劑盒（天根生化科技有限公司產(chǎn)品）說明書中的方法。利用RevertAid First Strand cDNA Synthesis Kit試劑盒反轉(zhuǎn)錄成第1鏈cDNA [24]。

1.7 杧果RAV基因家族表達分析

利用 QuantStudio 6 Flex的實時熒光定量PCR檢測系統(tǒng)檢測膠孢炭疽菌和細菌性黑斑病菌侵染下杧果RAV家族基因的表達量，引物見表1。反應(yīng)體系為18.0 μl，設(shè)置2個復(fù)孔。反應(yīng)程序參照UltraSYBR Mixture試劑盒說明書（北京康為世紀生物科技有限公司產(chǎn)品）進行設(shè)置。以杧果MiActin作為內(nèi)參基因，0 h的表達量為對照，運用2-△△Ct法進行數(shù)據(jù)統(tǒng)計，確定MiRAV1～MiRAV6的相對表達量[25]，并通過HemI 1.0.3.7軟件繪制熱圖，其中設(shè)置膠孢炭疽菌相對表達量熱圖參數(shù)（條形圖的數(shù)量）為11，小數(shù)點位數(shù)值為2，最小值為-20，階躍值為2;設(shè)置細菌性黑斑病菌相對表達量熱圖參數(shù)（條形圖的數(shù)量）為11，小數(shù)點位數(shù)值為2，最小值為-2，階躍值為2。整個表達矩陣每個值進行對數(shù)轉(zhuǎn)化，底數(shù)為2。

2 結(jié)果與分析

2.1 全基因組水平鑒定杧果RAV基因家族成員

在杧果基因組數(shù)據(jù)庫中進行Blastp檢索，發(fā)現(xiàn)杧果RAV基因家族有6個成員，分別命名為MiRAV1～MiRAV6。通過生物學(xué)在線分析工具ProtParam對杧果RAV基因家族6個成員MiRAV1～MiRAV6編碼的蛋白質(zhì)氨基酸序列的理化性質(zhì)進行了分析，結(jié)果（表2）表明，6個MiRAVs蛋白相對分子質(zhì)量為32 523 340～44 122 280，蛋白質(zhì)等電點為5.83～9.34，外顯子數(shù)均為2，脂溶指數(shù)為58.93～74.31，氨基酸數(shù)為281～392個，不穩(wěn)定指數(shù)為34.14～46.61，親水性平均數(shù)為-0.506～-0.863。因此，預(yù)測這6個MiRAVs蛋白中MiRAV1和MiRAV5為不穩(wěn)定親水堿性蛋白質(zhì)，MiRAV2～MiRAV3為不穩(wěn)定親水酸性蛋白質(zhì)，MiRAV4為穩(wěn)定親水酸性蛋白質(zhì)，MiRAV6為穩(wěn)定親水堿性蛋白質(zhì)。

對杧果RAV家族成員進行亞細胞定位預(yù)測，結(jié)果（表3）顯示，它們主要定位于細胞核中，其次定位于葉綠體和細胞質(zhì)中，推斷該家族成員主要存在于進行光合作用的器官中。其中，只有MiRAV1和MiRAV4沒有定位于葉綠體，MiRAV4和MiRAV6沒有定位于細胞質(zhì)，線粒體、過氧化氫酶體，質(zhì)膜定位的MiRAVs成員最少。

通過SOPMA在線軟件對MiRAV1～MiRAV6蛋白進行二級結(jié)構(gòu)預(yù)測。如圖1所示，6個MiRAVs蛋白無規(guī)則卷曲結(jié)構(gòu)占比為44.58%～54.85%，α-螺旋結(jié)構(gòu)占比為21.88%～31.67%，延伸鏈結(jié)構(gòu)占比為16.86%～19.13%，β-轉(zhuǎn)角結(jié)構(gòu)占比為4.85%～7.56%。因此，預(yù)測無規(guī)則卷曲和α-螺旋是6個MiRAVs蛋白二級結(jié)構(gòu)的主要元件，其次是延伸鏈和β-轉(zhuǎn)角。

通過NCBI Conserved Domain Search對杧果RAV家族蛋白質(zhì)氨基酸序列進行保守結(jié)構(gòu)域分析，發(fā)現(xiàn)MiRAV1～MiRAV6蛋白均含有AP2和B3超家族保守域（圖2）。6個MiRAVs蛋白AP2超家族保守域位置介于18～124;B3超家族保守域位置介于137～291。基于同源建模原理，通過SWISS-MODEL軟件對杧果RAV蛋白的2個保守結(jié)構(gòu)域（AP2、B3）三級結(jié)構(gòu)進行預(yù)測，發(fā)現(xiàn)在AP2保守結(jié)構(gòu)域中，含有1個 α-螺旋和3個反向平行的β-轉(zhuǎn)角;B3保守結(jié)構(gòu)域中，含有1個α-螺旋、2個反向平行的β-折疊和延伸鏈。

2.2 杧果RAV蛋白系統(tǒng)進化樹分析

為了更好地了解杧果RAV基因家族的系統(tǒng)發(fā)育關(guān)系，將杧果（6個）、擬南芥（7個）、煙草（5個）、蘋果（11個）、菠蘿（2個）、毛果楊（5個）RAV基因編碼的蛋白質(zhì)氨基酸序列通過Blast工具進行蛋白質(zhì)氨基酸序列分析，利用ClustalW程序和MEGA 7.0軟件構(gòu)建6個杧果RAV蛋白與其他物種（共5種）植物RAV蛋白系統(tǒng)進化樹（圖3），將36個RAV成員分成3個分支（Class Ⅰ、Class Ⅱ和Class Ⅲ），結(jié)果顯示，6個物種的RAV蛋白家族成員相互嵌合在一起，表明植物RAV在進化過程中相對保守。每個分支中均存在2個杧果RAV成員，其中，Class Ⅰ包含杧果2個、擬南芥2個、毛果楊3個、蘋果5個和煙草1個;Class Ⅱ包含杧果2個、蘋果4個、菠蘿2個和擬南芥5個;Class Ⅲ包含杧果2個、毛果楊2個、煙草4個和蘋果2個。在36個RAV成員中，MiRAV1、MiRAV6分別與MDP0000128924、MDP0000939633親緣關(guān)系最近，MiRAV2、MiRAV3分別與AtRAV1、AtRAV7親緣關(guān)系最近，MiRAV4與MDP0000207722、MDP0000485280親緣關(guān)系最近，MiRAV5與Potri.018G109200.1、Potri.018G109200.2親緣關(guān)系最近。因此，推測RAV家族蛋白質(zhì)不僅具有多樣性，對于親緣關(guān)系較近的蛋白質(zhì)，它們還具有相似或相近的生物學(xué)功能。

2.3 杧果RAV家族蛋白質(zhì)motif和保守結(jié)構(gòu)域分析

基于在線軟件MEME對杧果、擬南芥、煙草、蘋果、菠蘿、毛果楊RAV轉(zhuǎn)錄因子家族的蛋白質(zhì)保守基序進行分析（圖4），在RAV轉(zhuǎn)錄因子家族中鑒定出5個保守基序，主要的保守基序有4個（motif1～motif4），N端均存在保守的motif1，C端均存在保守的motif2。其中，MiRAV4～MiRAV5、AtRAV3～AtRAV4、MDP0000526584、MDP0000534780、MDP0000165802、MDP0000207722、MDP0000485280、Potri.018G109200.1、Potri.018G109200.2、Potri.003G212800.1、XP_016492322.1不含有motif5。通過基序分析發(fā)現(xiàn)，RAV基因家族在進化上結(jié)構(gòu)具有保守性。

通過DNAMAN 6.0軟件對杧果6個RAV蛋白和擬南芥7個、煙草5個、蘋果11個、菠蘿2個、毛果楊5個RAV蛋白進行多序列比對。結(jié)果（圖5）發(fā)現(xiàn)，MiRAV1～MiRAV6與5個物種的RAV轉(zhuǎn)錄因子家族存在2個相似的保守域A和保守域B，保守域A與motif1 基本重疊，保守域B與 motif2、motif3 基本重疊。因此，可以推測motif1、motif2和motif3是構(gòu)成AP2、B3結(jié)構(gòu)域的主要結(jié)構(gòu)。

2.4 杧果RAV基因家族的脅迫響應(yīng)

以杧果MiActin作為內(nèi)參基因，0 h處理作為對照，利用qRT-PCR檢測并分析了2種病原菌侵染過程中杧果RAV家族基因的表達程度。qRT-PCR結(jié)果顯示，MiRAV1～MiRAV6和內(nèi)參基因MiActin的熔解曲線為單峰，熔解溫度分別為81.69 ℃、74.15 ℃、71.79 ℃、73.75 ℃、77.79 ℃、80.52 ℃、83.73 ℃，將qRT-PCR產(chǎn)物經(jīng)2%的瓊脂糖凝膠進行電泳，結(jié)果（圖6）顯示，目標條帶單一，清晰，與預(yù)計大小一致，說明所設(shè)計的引物具有特異性;杧果膠孢炭疽菌侵染過程中（圖7A），MiRAV1～MiRAV6的相對表達量均顯著下調(diào)，其中MiRAV1～MiRAV6相對表達量在6 h時達到最低值，MiRAV1～MiRAV5的相對表達量在72 h時達到最高值，MiRAV6的相對表達量在48 h時達到最高值;杧果細菌性黑斑病菌侵染過程中（圖7B），MiRAV1～ MiRAV6的相對表達量在3 h時均顯著下調(diào)，MiRAV1、MiRAV2、MiRAV3和 MiRAV6的相對表達量在6 h時顯著上調(diào)，MiRAV1、MiRAV5和 MiRAV6的相對表達量分別在12 h、24 h、48 h、72 h時均顯著上調(diào)。結(jié)果表明，在不同病原菌侵染過程中，MiRAV1～ MiRAV6的相對表達量與對照差異顯著。

3 討論

基因組測序技術(shù)的提高和分子生物學(xué)研究的深入，為植物基因家族的鑒定、功能基因的挖掘提供了基礎(chǔ)，大量植物全基因組已經(jīng)被成功測序。RAV家族是一類植物特異的轉(zhuǎn)錄因子，廣泛存在于植物細胞內(nèi)。在擬南芥[12-13]、黃瓜[14]、棉花[15]、東方山羊豆[16]、煙草[17]、茶樹[18]等中關(guān)于RAV已開展相應(yīng)研究，主要體現(xiàn)在參與應(yīng)對外界脅迫和葉片衰老等方面功能，但關(guān)于其杧果RAV轉(zhuǎn)錄因子家族的研究至今仍未見報道。

本研究在杧果全基因組水平鑒定到6個RAV基因，理化性質(zhì)分析結(jié)果表明，MiRAV1和MiRAV5為不穩(wěn)定親水堿性蛋白質(zhì)，MiRAV2和MiRAV3為不穩(wěn)定親水酸性蛋白質(zhì)，MiRAV4為穩(wěn)定親水酸性蛋白質(zhì)，MiRAV6為穩(wěn)定親水堿性蛋白質(zhì)。杧果RAV蛋白主要存在于細胞核中，無規(guī)則卷曲和α-螺旋是6個MiRAVs蛋白二級結(jié)構(gòu)的主要元件，均具有AP2和B3超家族保守域結(jié)構(gòu)，AP2保守域三級結(jié)構(gòu)含有1個 α-螺旋和3個反向平行的β-轉(zhuǎn)角，B3保守域三級結(jié)構(gòu)含有1個 α-螺旋、2個反向平行的 β-折疊和延伸鏈，據(jù)此推斷它們在調(diào)控生物性狀上存在協(xié)同作用。杧果與擬南芥、煙草、蘋果、菠蘿、毛果楊RAV蛋白家族構(gòu)建系統(tǒng)進化樹結(jié)合基序分析結(jié)果顯示，6個物種的RAV蛋白家族成員相互嵌合在一起，表明植物RAV基因在進化過程中相對保守;在模式植物擬南芥中參與應(yīng)對外界脅迫和葉片衰老的基因AtRAV1和AtRAV2被聚類到ClassⅡ中，因此，推斷該ClassⅡ分支下的杧果MiRAV2和MiRAV3基因可能與應(yīng)對外界脅迫和葉片衰老有關(guān)[12-13];對于親緣關(guān)系較近的蛋白質(zhì)，推斷它們還具有相似或相近的生物學(xué)功能;杧果RAV蛋白主要的保守基序有4個（motif1～ motif 4），N端均存在保守的motif1，C端均存在保守的motif2，motif1～ motif 3與盧合均等[15]研究棉花時所鑒定到的motif1～motif3基本一致，motif1與韓林賀等[26]研究普通煙草時所鑒定到的motif1基本一致;杧果與其他5個物種RAV蛋白進行多序列比對結(jié)果顯示，存在2個明顯的保守域A和保守域B，保守域A與motif1基本重疊，保守域B與motif2、motif3基本重疊，因此，可以推測motif1是構(gòu)成AP2的主要結(jié)構(gòu)，motif2和motif3是構(gòu)成B3的主要結(jié)構(gòu)。

中國杧果病害危害嚴重，常見的病害有炭疽病、細菌性黑斑病、白粉病、蒂腐病、瘡痂病等[27-28]，其中，杧果膠孢炭疽病與細菌性黑斑病為害杧果最為嚴重[29-30]。為了解杧果RAV基因在膠孢炭疽菌與細菌性黑斑病菌侵染下相對表達量的變化，對杧果葉片進行膠孢炭疽菌與細菌性黑斑病菌生物脅迫處理，并通過qRT-PCR技術(shù)研究膠孢炭疽菌與細菌性黑斑病菌侵染過程中杧果RAV家族基因的相對表達量。結(jié)果顯示，在不同病原菌侵染過程中，MiRAV1～MiRAV6的相對表達量與對照差異顯著，表明MiRAV1～MiRAV6對生物脅迫有響應(yīng)。RAV基因家族的表達分析研究，更多集中在MeJA[16]、NaCl[18]、ABA和低溫[26]等非生物脅迫上，也有前人在研究中提到RAV基因?qū)γ藁S萎病菌[15]、番茄青枯病菌[31]等生物脅迫的響應(yīng)，本研究通過對杧果接種膠孢炭疽菌與細菌性黑斑病菌后的分析，發(fā)現(xiàn)了此家族的部分成員在膠孢炭疽菌與細菌性黑斑病菌侵染植株之后出現(xiàn)上調(diào)表達的現(xiàn)象， 表明RAV基因可能在杧果提高抗炭疽病和細菌性黑斑病的能力方面起著積極作用，如MiRAV1、MiRAV5、MiRAV6在細菌性黑斑病菌侵染后的12～72 h顯著上調(diào)表達，可以作為后續(xù)研究杧果抗細菌性黑斑病機制的候選基因，這也在分子水平上拓寬了杧果響應(yīng)生物脅迫的機制。

參考文獻：

[1] OHAMA N， SATO H， SHINOZAKI K， et al. Transcriptional regulatory network of plant heat stress response[J]. Trends in Plant Science， 2017， 22（1）： 53-65.

[2] KIM Y S， AN C， PARK S， et al. CAMTA-mediated regulation of salicylic acid immunity pathway genes in arabidopsis exposed to low temperature and pathogen infection[J]. The Plant Cell， 2017， 29（10）： 2465-2477.

[3] PRASAD K V S K， ABDRL-HAMEED A A E， XING D， et al. Global gene expression analysis using RNA-seq uncovered a new role for SR1/CAMTA3 transcription factor in salt stress[J]. Scientific Reports， 2016， 6（1）： 443-448.

[4] SCHWECHHEIMER C， BEVAN M. The regulation of transcription factor activity in plants[J]. Trend in Plant Science， 1998， 3（10）： 278-283.

[5] GIRAUDAT J， HAUGE B M， VALON C， et al. Isolation of the Arabidopsis ABI3 gene by positional cloning[J]. The Plant Cell， 1992， 4（10）： 1251-1261.

[6] SUZUKI M， KAO C Y， MCCARTY D R. The conserved B3 domain of VIVIPAROUS1 has a cooperative DNA binding activity[J]. The Plant Cell Online， 1997， 9（5）： 799-807.

[7] RIECHMANN J L， LAIN S， GARCIA J A. Highlights and prospects of potyvirus molecular biology[J]. Journal of General Virology， 1992， 73（1）： 1-16.

[8] KAGAYA Y， OHMIYA K， HATTORI T. RAV1， a novel DNA-binding protein， binds to bipartiterecognition sequence through two distinct DNA-binding domains uniquely found in higher plants[J]. Nucleic Acids Res， 1999， 27（2）： 470.

[9] RIECHMANN J L， HEARD J， MARTIN G， et al. Arabidopsis transcription factors：genome-widecomparative analysis among eukaryotes[J]. Science， 2000， 290（10）： 2105.

[10]FINKELSTEIN R R， WANG M L， LYNCH T J， et al. The arabidopsis abscisic acid response locusABI4 encodes an APETALA2 domain protein[J]. Plant Cell， 1998， 10（6）： 1043.

[11]孫 堯，孫 鑫，王 雷. 毛果楊RAV/PLC基因生物信息學(xué)分析[J].廣東農(nóng)業(yè)科學(xué)， 2019， 46（12）： 36-41，153.

[12]HU Y X， WANG Y X， LIU X F， et al. Arabidopsis RAV1 is down-regulated by brassinosteroid and may act as a negative regulator during plant development[J].Cell Research， 2004， 14（1）： 8-15.

[13]SOHN K H， LEE S C， JUNG H W， et al. Expression and functional roles of the pepper pathogen-induced transcription factor RAV1 in bacterial disease resistance，and drought and salt stress tolerance[J].Plant Molecular Biology， 2006， 61（6）： 897-915.

[14]蔣倫偉，胡麗芳，袁永成，等. 黃瓜RAV基因家族的全基因組分析[J].中國蔬菜， 2012， 1（16）： 15-21.

[15]盧合均，周忠麗，陳浩東，等. 棉花RAV基因家族的全基因組分析[J].棉花學(xué)報， 2014， 26（6）： 471-482.

[16] CHEN X， WANG Z， WANG X， et al. Isolationand characterization of GoRAV， a novel geneencoding a RAV-type protein in Galegae orientalis[J].Genes Genet Syst， 2009， 84（2）： 101.

[17]SRINIVASAN C， LIU Z R， HEIDMANN I， et al. Heterologous expression of the BABY BOOM AP2/ERF transcription factor enhances the regeneration capacity of tobacco （Nicotiana tabacum L.）[J].Planta， 2007， 225： 341-351.

[18]陳林波，李葉云，王 琴，等. 茶樹冷誘導(dǎo)基因RAV的克隆與表達特性分析[J].植物生理學(xué)通訊， 2010， 46（4）： 354.

[19]PENG W， YING F L， JIAN F H， et al. The genome evolution and domestication of tropical fruit mango[J]. Genome Biology， 2020， 21（1）. Doi：10.1186/s13059-020-01959-8.

[20]FINN R D， COGGILL P， EBERHARDT R Y， et al. The pfam protein families database： towards a more sustainable future[J].Nucleic Acids Research， 2016， 44：279-285.

[21]GASTEIGER E， GATTIKER A， HOOGLAND C， et al. ExPASy： the proteomics server for in-depth protein knowledge and analysis[J].Nucleic Acids Research， 2003， 31（13）： 3784-3788.

[22]BAILEY T L， BODEN M， BUSKE F A， et al. MEME SUITE： tools for motif discovery and searching[J].Nucleic Acids Research， 2009， 37： 202-208.

[23]KUMAR S， STECHER G， TAMUR K. MEGA7： molecular evolutionary genetics analysis Version 7.0 for bigger datasets [J].Molecular Biology and Evolution， 2016， 33（7）： 1870-1874.

[24]夏 楊，蘇初連，晁 駿，等. 菠蘿VOZ轉(zhuǎn)錄因子序列特征及其對非生物脅迫的響應(yīng)[J].西北植物學(xué)報， 2018， 38（7）： 1228-1234.

[25]LIVAK K J， SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2 [-Delta Delta C（T）] method [J].Methods-A Companion to Methods in Enzymology， 2001， 25（4）： 402-408.

[26]韓林賀，徐同慶，任昂彥，等.普通煙草RAV家族基因結(jié)構(gòu)與功能分析[J].基因組學(xué)與應(yīng)用生物學(xué)， 2019， 38（4）： 1658-1665.

[27]蒲金基，張 賀，周文忠. 芒果病害綜合防治技術(shù)[J].中國熱帶農(nóng)業(yè)， 2015（3）： 38-42.

[28]楊鄭州，黃柳芳，謝曉娜，等. 葉疫病菌侵染芒果后葉片細胞壁降解酶活性測定[J]. 江蘇農(nóng)業(yè)科學(xué)，2019，47（10）：138-141.

[29]喻群芳，漆艷香，張輝強，等. 4種生防菌對芒果細菌性黑斑病的田間防效[J].中國熱帶農(nóng)業(yè)， 2019（3）： 23-25.

[30]于海英，蘭建強，王曉燕，等. 芒果膠孢炭疽菌致病性的初步研究[C]//中國植物病理學(xué)會. 中國植物病理學(xué)會2012年學(xué)術(shù)年會論文集. 北京：中國農(nóng)業(yè)科學(xué)技術(shù)出版社， 2012： 210-213.

[31]LI C W， SU R C， CHENG C P， et al. Tomato RAV transcription factor is a pivotal modulator involved in the AP2/EREBP-mediated defense pathway [J].Plant Physiology， 2011， 156（1）：213-227.

（責任編輯：陳海霞）

收稿日期：2020-11-25

基金項目：國家重點研發(fā)專項（2019YFD1000504）;中國熱帶農(nóng)業(yè)科學(xué)院基本科研業(yè)務(wù)費專項（1630042017019）

作者簡介：孫 宇（1995-），男，海南三亞人，碩士研究生，研究方向為植物保護。（E-mail）sunyuqp@qq.com

通訊作者：張 賀，（E-mail）atzzhef@163.com;蒲金基，（E-mail）cataspjj@163.com