澳洲堅(jiān)果MiMYB2基因克隆及結(jié)構(gòu)與功能分析

王文林 陳海生* 鄭樹芳 樊松樂 王立豐 譚秋錦覃振師 黃錫云 賀 鵬 湯秀華 許 鵬

（1. 廣西南亞熱帶農(nóng)業(yè)科學(xué)研究所，龍州 532400；2. 海南大學(xué)熱帶作物學(xué)院，海南省熱帶生物資源可持續(xù)利用重點(diǎn)實(shí)驗(yàn)室，?？?70228；3. 農(nóng)業(yè)農(nóng)村部橡膠樹生物學(xué)與遺傳資源利用重點(diǎn)實(shí)驗(yàn)室；省部共建國家重點(diǎn)實(shí)驗(yàn)室培育基地-海南省熱帶作物栽培生理學(xué)重點(diǎn)實(shí)驗(yàn)室；農(nóng)業(yè)農(nóng)村部儋州熱帶作物科學(xué)觀測實(shí)驗(yàn)站；中國熱帶農(nóng)業(yè)科學(xué)院橡膠研究所，海口 571101）

澳洲堅(jiān)果（Macadamia integrifolia）是一種原產(chǎn)于澳大利亞的常綠喬木果樹，其果仁含有豐富的不飽和脂肪酸、蛋白質(zhì)、氨基酸、可溶性糖、膳食纖維和礦物質(zhì)等營養(yǎng)成分，具有藥用和保健等功能［1］。澳洲堅(jiān)果果仁可加工成飲品、食品和果油等［2-3］。澳洲堅(jiān)果副產(chǎn)品，比如果殼可用于制備活性炭［4］和吸附劑［5］，油粕可用于生產(chǎn)動(dòng)物飼料［1］和多肽［6］等，經(jīng)加工后可提高其經(jīng)濟(jì)附加值。我國澳洲堅(jiān)果品種主要是經(jīng)過試種馴化的國外優(yōu)良品種。國外優(yōu)良品種在我國試種期間，經(jīng)常面臨產(chǎn)量及品質(zhì)降低的問題［7］。這是由于我國澳洲堅(jiān)果主要種植地域分布在廣西、云南、貴州和四川，隨著環(huán)境的變化，各種病蟲害如花疫病、枝枯病、炭疽病、蚜蟲和蛾類等［8］，干旱［9］、磷［10］、低溫［11］等不可預(yù)測的非生物脅迫和生物脅迫對澳洲堅(jiān)果的生長發(fā)育、品質(zhì)和產(chǎn)量均有影響，制約了其產(chǎn)業(yè)的發(fā)展。

MYB轉(zhuǎn)錄因子是一個(gè)植物轉(zhuǎn)錄因子家族。其中MYB 家族成員R2R3-MYBs 參與了多種生物過程，包括生長發(fā)育、激素信號轉(zhuǎn)導(dǎo)、代謝、生物脅迫和非生物脅迫等。例如，AtMYB77 和ARF7 相互作用參與擬南芥生長素的信號轉(zhuǎn)導(dǎo)［12］。AtMYB32參與擬南芥花粉正常的發(fā)育［13］。楊樹MYB165 和MYB194屬于R2R3-MYB家族，啟動(dòng)子激活實(shí)驗(yàn)表明MYB194 和MYB165 抑制原花青素生物合成相關(guān)基因啟動(dòng)子激活，過表達(dá)MYB165 和MYB194 轉(zhuǎn)基因楊樹導(dǎo)致葉片花青素積累減少，根系原花青素含量降低［14］。R2R3-MYB 家族轉(zhuǎn)錄因子成員AtMYB102 對 昆 蟲 草 食 動(dòng) 物 菜 青 蟲 有 抗 性［15］。PbrMYB21 在耐旱性方面發(fā)揮著積極的作用［16］。迄今，澳洲堅(jiān)果轉(zhuǎn)錄調(diào)控方面的研究和MYB 家族轉(zhuǎn)錄因子成員的結(jié)構(gòu)與功能信息研究尚屬空白。研究R2R3-MYB 家族成員的結(jié)構(gòu)與功能為闡明澳洲堅(jiān)果抗逆機(jī)制和提高澳洲堅(jiān)果品質(zhì)和產(chǎn)量，促進(jìn)產(chǎn)業(yè)發(fā)展具有重要意義。本研究從“桂熱一號”品種的葉片中克隆MiMYB2 基因，利用生物信息學(xué)對其結(jié)構(gòu)和功能進(jìn)行分析。澳洲堅(jiān)果R2R3-MYB 家族轉(zhuǎn)錄因子成員的結(jié)構(gòu)和功能的深入研究，將為闡明澳洲堅(jiān)果生長發(fā)育和產(chǎn)量形成機(jī)制提供理論參考。

1 材料與方法

1.1 試驗(yàn)材料

試驗(yàn)材料“桂熱一號”和“695”品種由廣西南亞熱帶農(nóng)業(yè)科學(xué)研究所澳洲堅(jiān)果種質(zhì)圃（E106.79.84，N22.14.34）提 供。“桂 熱 一 號”和“695”澳洲堅(jiān)果品種嫁接苗種植于廣西南亞熱帶農(nóng)業(yè)科學(xué)研究所中，2019 年3 月選取2008 年嫁接種植，樹齡為11 年兩品種澳洲堅(jiān)果成年樹健康植株的枝條、花、葉為材料進(jìn)行轉(zhuǎn)錄組分析。

1.2 試驗(yàn)方法

1.2.1 提取總RNA和cDNA的合成

按照總RNA 提取試劑盒提供的說明書，提取澳洲堅(jiān)果“桂熱一號”品種葉片的總RNA，RNA 和反轉(zhuǎn)錄得到cDNA 的濃度與純度由Thermo Fisher NanoDrop 2000 超微量核酸蛋白分析儀檢測（Gene Company Limited，上海）。按照TaKaRa 反轉(zhuǎn)錄試劑盒的說明書合成cDNA，1%瓊脂糖凝膠電泳檢測RNA完整性。

1.2.2 克隆MiMYB2基因全長

根據(jù)澳洲堅(jiān)果轉(zhuǎn)錄組測序序列，利用primer primer 6.0 軟 件 設(shè) 計(jì) 引 物MiMYB2F，5′-CGGTCCTTCTCAGTTTCTGC-3′，MiMYB2 R，5′-GCATCCTAACATCCCTGGAA-3′，引物序列由公司合成。以澳洲堅(jiān)果“桂熱一號”品種正常葉片的cDNA 為模板，利用PCR 擴(kuò)增MiMYB2 的全長序列，PCR 反應(yīng)體系20.0 μL：69 μg·L-1cDNA 模板1.0 μL，2×EX-Taq PCR Master Mix 10.0 μL，上、下游引物各0.5 μL，ddH2O 補(bǔ)足至20.0 μL。擴(kuò)增程序：95℃預(yù)變性3 min；95℃30 s，55℃50 s，72℃，80 s，32 個(gè)循環(huán)；72℃10 min。將PCR 產(chǎn)物加到0.8%瓊脂糖凝膠中進(jìn)行電泳檢測，切割目的條帶后回收純化，目的片段連接pMD-18T 載體，16℃水浴 4 h 后，轉(zhuǎn)化大腸桿菌DH5α 感受態(tài)細(xì)胞中，恢復(fù)培養(yǎng)后均勻涂布在含羧芐抗性的培養(yǎng)基上，經(jīng)PCR檢測后挑取陽性克隆測序。

1.2.3 MiMYB2生物信息學(xué)分析

利用在線分析工具ProtParam（https：//web.expasy.org/protparam/）、TMPred（https：//embnet.vitalit.ch/software/TMPRED_form.html）、SMART：Main page （http：//smart. embl-heidelberg. de/smart/set_mode.cgi？NORMAL=1）、PHYRE、PlantmPLoc2、PSIPRED V4.0（http：//bioinf.cs.ucl.ac.uk/psipred/）、SignalP-5.0 Server（http：//www.cbs.dtu.dk/services/SignalP/）和 MEME（http：//meme-suite. org/tools/meme）等對MiMYB2 進(jìn)行生物信息學(xué)預(yù)測分析。利用Clustal Omega 在線工具軟件進(jìn)行多序列比對。MEGAX 軟件構(gòu)建系統(tǒng)進(jìn)化樹（neighbor-joining，bootstrap 值設(shè)為1 000）。用于系統(tǒng)進(jìn)化分析的擬南芥和其他物種的蛋白序列均下載自NCBI。

2 結(jié)果與分析

2.1 MiMYB2基因克隆

按照試劑盒操作方法提取澳洲堅(jiān)果“桂熱一號”品種葉片的總RNA，從圖1 可見RNA 28S，18S和5S 條帶清晰完整且沒有發(fā)生降解（見圖1A）。以cDNA 為模板，采用PCR 擴(kuò)增目的基因，結(jié)果如圖1B 所示，擴(kuò)增出了與預(yù)期大小相符的特異性條帶。以澳洲堅(jiān)果“桂熱一號”品種正常葉片的cDNA 為模板，利用PCR 技術(shù)克隆得到MiMYB2 cDNA 基因全長，測序驗(yàn)證后，將cDNA 序列命名為MiMYB2。將MiMYB2 基因序列上傳至NCBI 數(shù)據(jù)庫，GenBank 登陸號為MN254976，MiMYB2 基因的cDNA 全長序列1 210 bp，包含1 002 bp 的開放閱讀框（ORF）如圖2 下劃線標(biāo)記序列，編碼333 個(gè)氨基酸。

圖1 澳洲堅(jiān)果桂熱一號RNA檢測及基因克隆結(jié)果A.RNA電泳圖；B.MiMYB2基因PCR擴(kuò)增結(jié)果（M.Marker DL2000；1. 目的基因擴(kuò)增產(chǎn)物）Fig.1 Total RNA of M.ntegrifolia（GUIRE1）and the PCR amplification products of MiMYB2 geneA. Total RNA of M. integrifolia（GUIRE1）；B. The PCR amplification products of MiMYB2 gene. M：Marker DL200；1. The objective gene amplification products

2.2 MiMYB2蛋白質(zhì)的理化性質(zhì)預(yù)測

利用ProtParam 工具在線分析MiMYB2 蛋白質(zhì)的理化性質(zhì)：分子式C1556H2477N447O511S19，總原子數(shù)5 010，分子量為36.231 61 kDa，等電點(diǎn)5.70，帶負(fù)電荷殘基總數(shù)（Asp+Glu）為42，帶正電殘基總數(shù)（Arg+Lys）為37，不穩(wěn)定系數(shù)為67.68（該蛋白屬于不穩(wěn)定蛋白），脂肪族氨基酸指數(shù)是67.09，總平均親水性為-0.507，推測MiMYB2蛋白是一個(gè)不穩(wěn)定的親水蛋白。

2.3 MiMYB2蛋白質(zhì)的磷酸化位點(diǎn)、跨膜結(jié)構(gòu)及亞細(xì)胞定位預(yù)測

圖2 MiMYB2基因測序結(jié)果Fig.2 Sequence result of MiMYB2 gene

圖3 MiMYB2蛋白的磷酸化位點(diǎn)和跨膜結(jié)構(gòu)預(yù)測結(jié)果A.MiMYB2蛋白的磷酸化位點(diǎn)預(yù)測結(jié)果；B.MiMYB2跨膜結(jié)構(gòu)預(yù)測Fig.3 Prediction of phosphorylation sites and transmembrane structures of MiMYB2 proteinA.Prediction of phosphorylation site of MiMYB2 protein；B.Prediction of transmembrane domain of MiMYB2

利用NetPhos 3.1 在線分析工具預(yù)測MiMYB2蛋白磷酸化位點(diǎn)，結(jié)果表明，MiMYB2 氨基酸序列有可能發(fā)生磷酸化的氨基酸位點(diǎn)共有53 個(gè)，絲氨酸磷酸化位點(diǎn)44個(gè)，蘇氨酸磷酸化位點(diǎn)7個(gè)，酪氨酸磷酸位點(diǎn)2 個(gè)（見圖3A）。根據(jù)MiMYB2 蛋白磷酸化位點(diǎn)預(yù)測結(jié)果，推測MiMYB2 蛋白以絲氨酸為主的磷酸化修飾，調(diào)控澳洲堅(jiān)果的生物功能。利用TMPred 在線工具分析預(yù)測MiMYB2 跨膜結(jié)構(gòu)，結(jié)果表明MiMYB2 蛋白無跨膜結(jié)構(gòu)，推測Mi-MYB2 是非跨膜蛋白。利用SignalP-5.0 Server 分析預(yù)測MiMYB2 蛋白存在信號肽的概率為0.09%（見圖3B），推測MiMYB2 是非分泌型蛋白。利用PlantmPLoc2 預(yù)測分析MiMYB2 蛋白亞細(xì)胞定位，預(yù)測結(jié)果顯示其定位于細(xì)胞核中，其進(jìn)一步證明MiMYB2 編碼一個(gè)轉(zhuǎn)錄因子蛋白。

2.4 MiMYB2 蛋白質(zhì)保守結(jié)構(gòu)域、二級結(jié)構(gòu)、三級結(jié)構(gòu)分析

利用SMART：Main page 分析預(yù)測MiMYB2 蛋白的保守結(jié)構(gòu)域，發(fā)現(xiàn)其存在2 個(gè)SANT 結(jié)構(gòu)域和5 個(gè)低復(fù)雜度其位置排序分別在13-63、66-114、128-143、168-179、214-226、235-246 和253-284 氨基酸處，MiMYB2 氨基酸序列具有兩個(gè)SANT 保守結(jié)構(gòu)域（見圖4B），推測其屬于R2R3-MYB 家族。利用SOPMA 分析MiMYB2 蛋白的二級結(jié)構(gòu)，結(jié)果顯示α 螺旋占37.54%，延長鏈占10.81%，β 折疊占6.31%，無規(guī)卷曲占45.35%（見圖4A）。利用PHYRE 分析MiMYB2 的三級結(jié)構(gòu)如圖4C 所示。據(jù)此推測MiMYB2 是一個(gè)轉(zhuǎn)錄因子并且屬于R2R3-MYB家族成員。

圖4 MiMYB2結(jié)構(gòu)分析A.MiMYB2二級結(jié)構(gòu)預(yù)測；B.MiMYB2保守結(jié)構(gòu)域預(yù)測；C.MiMYB2三級結(jié)構(gòu)預(yù)測Fig.4 Structural analysis of MiMYB2A.PredictionforthesecondarystructureofMiMYB2；B.PredictionoftheconserveddomainsofMiMYB2；C.PredictionforthetertiarystructureofMiMYB2

2.5 MiMYB2系統(tǒng)進(jìn)化分析

如圖5 所示，澳洲堅(jiān)果MiMYB2（M.integrifolia，MN254975）與荷花NnMYB3-like（Nelumbo nucifera，XP_010255443.1）、荷花Nnmyb-related protein 308-like（Nelumbo nucifera，XP_010273102.1）、醉 蝶 花 ThMYB6（Tarenaya hassleriana，XP_010523735.1）、麻 風(fēng) 樹JcMYB6（Jatropha curcas，XP_012083427.1）和哥倫比亞錦葵HuMYB41（Herrania umbratica，XP_021280626.1）進(jìn)行同源性分析，相似性分別為64.79%，65.99%，59.88%，58.08%和56.18%，其與荷花NnMYB3-like 同源性最高。Clustal 進(jìn)行多序列比對分析并根據(jù)圖4B分析的結(jié)果標(biāo)注MiMYB2 保守結(jié)構(gòu)域SANT（見圖5A）。利用MEME 分析澳洲堅(jiān)果MiMYB2 與其在NCBI數(shù)據(jù)庫中比對親緣性較高的MYB蛋白序列，預(yù)測了排名前三的motif 且標(biāo)注了其在相應(yīng)序列的位置，如圖5B 所示澳洲堅(jiān)果MiMYB2 a（77-126）、b（1-50）、c（339-367），同源性越高的序列，其motif 的大小及位置越接近。Motif 與植物的生物學(xué)功能息息相關(guān)。將MiMYB2 編碼的蛋白序列與BLAST 比對結(jié)果中相似性較高的已知其他植物的MYB 序列進(jìn)行系統(tǒng)進(jìn)化分析，表明MiMYB2 與荷花 NnMYB3-like （Nelumbo nucifera， XP_010255443.1）、荷花Nnmyb-related protein 308-like（Nelumbo nucifera，XP_010273102.1）親緣關(guān)系較近（見圖6A）。系統(tǒng)進(jìn)化分析澳洲堅(jiān)果MiMYB2與126個(gè)擬南芥R2R3-MYB家族成員的蛋白序列，如圖6B 結(jié)果顯示MiMYB2 與AtMYB17、AtMYB106、AtMYB91、AtMYB16 親緣關(guān)系較近，AtMYB17、At-MYB106、AtMYB16 屬于擬南芥S9 亞族［17］，推測MiMYB2具有與S9亞族成員相似的功能。

2.6 MiMYB2基因組織表達(dá)模式分析

為研究MiMYB2 基因在“桂熱一號”和“695”品種不同組織中的表達(dá)模式，本研究從“桂熱一號”和“695”品種澳洲堅(jiān)果轉(zhuǎn)錄組數(shù)據(jù)中提取了MiMYB2基因在枝條、花、和葉片的表達(dá)數(shù)據(jù)，分析其在不同品種和不同組織的表達(dá)模式。如圖7 所示，“桂熱一號”品種，MiMYB2 主要在枝條和葉片中表達(dá)，花中的表達(dá)量較低，葉片中MiMYB2 基因表達(dá)量是花的8 倍左右；“695”品種中，MiMYB2 主要在花中表達(dá)，枝條和葉片中表達(dá)量較低，在花中的表達(dá)量是葉片的10 倍左右?！肮馃嵋惶枴焙汀?95”品種同一組織MiMYB2 基因的表達(dá)在花和葉片中相差較大，特別是其在花中的表達(dá)量“695”是“桂熱一號”的13倍左右。

3 討論

圖6 澳洲堅(jiān)果MiMYB2系統(tǒng)進(jìn)化分析A. 植物中部分MYB轉(zhuǎn)錄因子的系統(tǒng)進(jìn)化樹；B. 澳洲堅(jiān)果MiMYB2蛋白與擬南芥中126 R2R3-MYB家族成員蛋白序列的系統(tǒng)樹聚類分析Fig.6 Phylogenetic analysis of M.integrifolia MiMYB2A. Phylogenetic tree of MiMYB2 and other plant MYB transcription factors；B. Phylogenetic tree of MiMYB2 from M. integrifolia and 126 R2R3-MYB protein members from Arabidopsis with indicated GenBank accession numbers

澳洲堅(jiān)果果仁具有豐富的脂肪酸，趙大宣等采用氣相色譜、質(zhì)譜聯(lián)用技術(shù)檢測到13 種脂肪酸成分，其中不飽和脂肪酸（油酸、異油酸、棕櫚酸）總量占了69.12%，飽和脂肪酸總量占30.87%，不飽和脂肪酸對冠心病、糖尿病等疾病具有效用［18］。楊為海等分析了不同種質(zhì)澳洲堅(jiān)果果仁礦質(zhì)元素含量，研究發(fā)現(xiàn)果仁中含有多種常量和微量元素，如K、Mg、Mn、Fe 平均含量分別為3.87、1.85、95.91和90.58 mg·kg-1［19］。澳洲堅(jiān)果果油含有生育酚、角鯊烯等多種抗氧化物質(zhì)［1］。澳洲堅(jiān)果青皮提取的多酚，以Trolox 為陽性對照評價(jià)其抗氧化性，研究發(fā)現(xiàn)其可用于制備天然抗氧化劑［20］。寧平等研究表明澳洲堅(jiān)果果殼可用于制備活性炭。涂行浩等利用微波輻照澳洲堅(jiān)果果殼制備的活性碳產(chǎn)品，其吸附性能超過國家一級標(biāo)椎［4，21］。楊雯等以澳洲堅(jiān)果果殼等為原料，利用超聲-7 堿組合技術(shù)制備的吸附劑可用于處理印染廢水［5］。澳洲堅(jiān)果及其副產(chǎn)品具有較高的營養(yǎng)價(jià)值和經(jīng)濟(jì)價(jià)值。我國澳洲堅(jiān)果品種主要是經(jīng)過試種馴化的國外優(yōu)良品種，優(yōu)良品種離開本土后其產(chǎn)量和品質(zhì)均有下降，不可預(yù)測的生物和非生物脅影響著澳洲堅(jiān)果的品質(zhì)和產(chǎn)量。迄今，澳洲堅(jiān)果的研究主要集中在栽培管理和雜交育種、誘變育種等，然而澳洲堅(jiān)果如何調(diào)控高產(chǎn)優(yōu)質(zhì)的機(jī)制及抗逆機(jī)制尚不清楚且研究較少。分子標(biāo)記育種技術(shù)結(jié)合雜交育種是選育澳洲堅(jiān)果優(yōu)良新品種的重要方法。宮麗丹等［22］研究澳洲堅(jiān)果耗水規(guī)律及灌溉制度，結(jié)果表明不同時(shí)期合理灌溉可提高產(chǎn)量和品質(zhì)。譚秋錦等［23］研究發(fā)現(xiàn)土壤養(yǎng)分因子（氮、鉀、pH、有機(jī)碳）對澳洲堅(jiān)果果實(shí)的品質(zhì)有影響，適當(dāng)提高土壤氮、鉀含量合理控制施肥種類可以提高其品質(zhì)。唐瑩瑩等［24］建立的澳洲堅(jiān)果RAPD 標(biāo)記反應(yīng)體系可用于分析澳洲堅(jiān)果遺傳多樣性及鑒定種質(zhì)資源。孔廣紅等［25］利用Co-γ 射線輻照不同種質(zhì)澳洲堅(jiān)果種子，篩選出了4 株性狀優(yōu)良的突變單株。Nock 等對廣泛種植的澳洲堅(jiān)果品種741 的基因組和轉(zhuǎn)錄組進(jìn)行了測序，將澳洲堅(jiān)果測序結(jié)果與其他5 種真雙子葉植物的基因家族比較，發(fā)現(xiàn)了13 689 個(gè)含有澳洲堅(jiān)果基因的簇和1 005 個(gè)澳洲堅(jiān)果特異簇。通過對參與氰苷生物合成的澳洲堅(jiān)果基因的分析，發(fā)現(xiàn)了幾個(gè)高表達(dá)的候選基因，為病原識別、植物防御和單萜合成中基因家族的功能研究打下了基礎(chǔ)［26］。Abubaker等［27］對埃及栽培的澳洲堅(jiān)果進(jìn)行了DNA 指紋圖譜分析，在DNA 指紋圖譜研究中記錄了16個(gè)不同的RAPD片段。

圖7 MiMYB2基因在不同品種和組織中的表達(dá)分析各柱形圖上用不同大寫字母標(biāo)識表示數(shù)據(jù)間差異極顯著（P＜0.01）Fig.7 MiMYB2 expression in different cultivars and tissuesDifferent uppercase letters above each column means significant at P＜0.01 level

R2R3-MYB 家族轉(zhuǎn)錄因子成員與植物的生長發(fā)育、生物脅迫和非生物脅迫密切相關(guān)。研究澳洲堅(jiān)果R2R3-MYB 家族成員的結(jié)構(gòu)與功能，將為闡明澳洲堅(jiān)果產(chǎn)量調(diào)控機(jī)制具有重要的意義。本研究通過PCR技術(shù)從澳洲堅(jiān)果“桂熱一號”的葉片克隆了MiMYB2 基因，測序分析cDNA 序列全長為1 210 bp，推導(dǎo)的氨基酸序列為1 002 bp ORF，編碼333 個(gè)氨基酸。生物信息學(xué)分析發(fā)現(xiàn)，MiMYB2具有兩個(gè)SANT保守結(jié)構(gòu)域，推測屬于屬于澳洲堅(jiān)果R2R3-MYB 家族；MiMYB2 定位于細(xì)胞核是一個(gè)無跨膜結(jié)構(gòu)、無信號肽的不穩(wěn)定親水蛋白。

MiMYB2 蛋白序列提交到NCBI 數(shù)據(jù)庫進(jìn)行BLAST 比對分析，結(jié)果表明其與荷花NnMYB3-like（Nelumbo nucifera，XP_010255443.1）同源性最高（見圖5A）。Deng 等對荷花（Nelumbo nucifera）中116個(gè)R2R3 MYB基因的基因結(jié)構(gòu)、表達(dá)模式等進(jìn)行分析，R2R3-MYB基因可分為4組并利用熒光定量PCR 分析了13 個(gè)與類黃酮生物合成相關(guān)的MYB 基因在荷花不同組織中的表達(dá)模式［28］。Sun等［29］利用酵母雙雜交試驗(yàn)證明，NnMYB5 與Nnb-HLH1、NlbHLH1 和NnTTG1 相互作用，NnTTG1 也與NnbHLH1 和NlbHLH1 相互作用。在擬南芥中過表達(dá)NnMYB5 導(dǎo)致擬南芥未成熟種子和花梗中的花青素積累，TT19表達(dá)上調(diào)，證明NnMYB5是花青素合成的轉(zhuǎn)錄激活因子。Peng 等［30］對麻瘋樹MYB 家族的全基因組鑒定，鑒定了128 個(gè)MYB 基因包括123 R2R3-MYBs，4 R1R2R3-MYBs 和1 4RMYB，選取JcMYB2 基因進(jìn)行功能研究，發(fā)現(xiàn)過表達(dá)JcMYB2 植物比野生型更耐鹽、耐冷脅迫。Mi-MYB2 氨基酸序列與126 個(gè)擬南芥R2R3-MYB 家族成員的系統(tǒng)進(jìn)化分析表明，MiMYB2 與At-MYB17、AtMYB106、AtMYB91、AtMYB16親緣關(guān)系較近（見圖6B），AtMYB17、AtMYB106、AtMYB16屬于擬南芥S9 亞族，將MiMYB2 聚類為S9 亞族。R2R3-MYB 家族轉(zhuǎn)錄因子成員，AtMYB17 與早期花序發(fā)育和種子萌發(fā)相關(guān)［31］。成熟擬南芥表皮毛的轉(zhuǎn)錄譜分析表明，NOECK 編碼MIXTA-Like 轉(zhuǎn)錄調(diào)控因子MYB106 與細(xì)胞生長有關(guān)［32］。白樺中一種新的R2R3-MYB 轉(zhuǎn)錄因子BpMYB106 通過上調(diào)光合作用和氧化磷酸化途徑相關(guān)基因表達(dá)，提高光合作用和生長速率［33］。擬南芥R2R3-MYB 家族成員MYB16 和MYB106 在生殖器官、營養(yǎng)器官和表皮毛中調(diào)控角質(zhì)層形成［34］。擬南芥S9 亞族成員主要與植物的生長發(fā)育相關(guān)。通過系統(tǒng)進(jìn)化分 析 將MiMYB2 與AtMYB17、AtMYB106、At-MYB16 聚類為S9 亞族。分析“桂熱一號”和“695”品種澳洲堅(jiān)果轉(zhuǎn)錄組數(shù)據(jù)中MiMYB2 基因在不同組織的表達(dá)模式。從圖7可知，“桂熱一號”中，Mi-MYB2 在花中的表達(dá)量最低，然而“695”中花中的表達(dá)量最高。花屬于植物的生殖器官，其形成與生長屬于生殖生長，成花的遲早接的影響了果樹的產(chǎn)量。Cao 等［35］通過熒光定量PCR 檢測了20 個(gè)可能參與木質(zhì)素生物合成的PbMYB 基因的表達(dá)模式，推測PbMYB25 和PbMYB52 是梨果實(shí)發(fā)育過程中參與木質(zhì)素合成調(diào)控的候選基因。蘋果MYB轉(zhuǎn)錄因子MdMYB3 參與花青素生物合成和花發(fā)育的轉(zhuǎn)錄調(diào)控，MdMYB39L 與花粉管生長相關(guān)［36-37］。在煙草中過表達(dá)桃子轉(zhuǎn)錄因子MYB10.1 導(dǎo)致了NtMYB305 的抑制，而NtMYB30 是花發(fā)育所必需的［38］。綜上，推測MiMYB2 參與調(diào)控澳洲堅(jiān)果生長發(fā)育及產(chǎn)量形成機(jī)制，有待深入研究。