秦玉芝，邢錚，潘妃，熊興耀,2*

(1.湖南農業(yè)大學園藝園林學院，湖南 長沙 410128；2.中國農業(yè)科學院蔬菜花卉研究所，北京 100081)

植物MYB 類轉錄因子以含有MYB 結構域為共同特征[1]。每個MYB 結構域以螺旋–轉角–螺旋(helix–tum–helix，HTH)的形式結合于目標DNA 大溝處，調控目的基因的轉錄[2]。第1個植物myb 基因ZmMYBc1 是Paz–Ares 于1987年在玉米中發(fā)現(xiàn)的，ZmMYBc1 與動物myb 基因的同源性達40%[3]。不同作物中myb 基因的數(shù)量各異，在玉米中表達的myb 基因約有80個以上[4]，擬南芥中MYB 類轉錄因子有129個[5]，馬鈴薯中MYB 類轉錄因子有117個。MYB 類轉錄因子中以R2R3–MYB 的數(shù)量最多，其功能涉及植物次生代謝調控、細胞形態(tài)發(fā)生、激素刺激和環(huán)境脅迫應答、分生組織形成及細胞周期控制等過程[6]。關于馬鈴薯MYB 類基因克隆和功能分析的報道尚少。在前期研究中，筆者發(fā)現(xiàn)了1個受光影響，并且對馬鈴薯花青素合成關鍵酶StCHS 基因表達具有明顯正向作用的MYB 轉錄因子(GenBank: DQ917781.1)[7]，本研究中通過同源克隆，在馬鈴薯原始栽培種Yan 中獲得了該轉錄因子(StR2R3– MYB1)的完整CDS 序列，并運用現(xiàn)代生物信息學技術對馬鈴薯StR2R3–MYB1 基因序列和蛋白質序列特征進行了比對、分析和預測，對基因表達的亞細胞定位、組織表達特異性及光照應答模式進行了分析，現(xiàn)將結果報道如下。

1 材料與方法

1.1 材料

馬鈴薯原始栽培種 Yan(S. tuberosum subsp. andigena var. yanacochense)由湖南農業(yè)大學馬鈴薯工程技術中心種質資源庫提供。擬南芥原生質體由湖南大學生命科學與技術研究院劉選明先生惠贈。表達載體pEGAD、pEGAD–MYC、PA7–YFP、大腸桿菌DH5α、農桿菌菌種Agl–0 均由湖南大學生命科學與技術研究院劉選明先生惠贈；T 載體pGM–T 購自天根生化公司。

1.2 方法

1.2.1 基因的克隆

根據(jù)表1 設計的引物，以由RNA 逆轉錄所得cDNA 為模板進行RT–PCR 。

表1 StR2R3–MYB1 克隆引物 Table 1 Primer sequences for StR2R3–MYB1 cloning

1.2.2 序列生物信息學分析的相關網(wǎng)站和軟件

馬鈴薯基因組數(shù)據(jù)庫網(wǎng)站：http://www. potatogenome.net/index.php/Main_Page。

蛋白結構域分析網(wǎng)站：http://www.Ncbi.nlm.nih.gov/。

轉錄調控區(qū)域元件分析網(wǎng)站：PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/)。

蛋白一級結構預測網(wǎng)站：ProtParam(http://web. expasy.org/protparam/)。

卷曲螺旋結構預測網(wǎng)站：COILSServer(http: //embnet.vital–it.ch/software/COILS_form.html)。

三級結構預測網(wǎng)站：SWISS–MODEL(http:// swissmodel.expasy.org/)。

跨膜區(qū)域分析網(wǎng)站：TMHMM(http://www.cbs. dtu.dk/services/TMHMM/)。

信號肽分析工具網(wǎng)站：SignalIP3.0(http://www. cbs.dtu.dk/services/SignalP/)。

親疏水性分析網(wǎng)站：http://www.expasy.ch/cgi– bin/protscale.pl。

同源性分析采用GeneDoc 軟件。進化樹構建采用MEGA5.0 軟件。

1.2.3 實時熒光定量PCR

定量PCR 反應程序：95℃預變性10min，95℃變性30 s，58℃退火30 s，72℃延伸20 s， 45個循環(huán)。

內參為ActinB。所有檢測樣品均設置3次重復。定量引物見表2。

表2 StR2R3–MYB1 實時熒光定量PCR 引物 Table 2 qPCR primer sequences for StR2R3–MYB1

1.2.4 擬南芥原生質體轉化

取20μg 重組質粒，加入100 μL 擬南芥原生質體和110 μL PEG 溶液，混勻，放置30min。加入440 μL W5 溶液，23℃、750 r/min 離心1min。

1.2.5 光處理

新繼代的組培苗置于黑暗、22℃培養(yǎng)箱中生長12 d，形成白化苗。將白化苗分別置于白光、藍光、紅光、遠紅光培養(yǎng)箱中，于放置后0、0.5、1、2、4、8、12、24 h 取樣，–80℃保存。藍光、紅光、遠紅光和白光光源分別為LED–藍光(波長為 470 nm)，LED–紅光(波長為 660 nm)，LED–遠紅光(波長為740 nm)和白色熒光燈。

2 結果與分析

2.1 馬鈴薯StR2R3–MYB1 基因的克隆與同源性分析

經RT–PCR 擴增，獲得1 條800 bp 的特異性片段。該片段包含1個全長798 bp 的CDS 序列，編碼265個氨基酸長度的蛋白質，命名為StR2R3– MYB1。馬鈴薯StR2R3–MYB1 蛋白序列與番茄、茄子的同源性高達99%，與擬南芥MYB113 的同源性也在70%以上(圖1～3)。

圖1 StR2R3–MYB1 氨基酸序列系統(tǒng)進化分析結果 Fig. 1 Homology analysis for amino acid sequence of StR2R3 – MYB1

圖2 StR2R3–MYB1 cDNA 序列與其氨基酸序列預測結果 Fig. 2 Prediction for cDNA and amino acid sequence of StR2R3–MYB1

圖3 StR2R3–MYB1 氨基酸序列同源性分析的結果 Fig.3 Homology analysis for amino acid sequence of StR2R3–MYB1

2.2 StR2R3–MYB1 蛋白序列結構的分析

StR2R3–MYB1 蛋白全長265 aa，主要氨基酸為天冬氨酸(9.8%)、亮氨酸(9.4%)、甘氨酸(7.5%)和谷氨酸(6.8%)，正/負電荷殘基數(shù)分別為32 和32，相對分子質量為30.23，分子式為C1316H2072N386O399S17，等電點為6.95，不穩(wěn)定系數(shù)為47.52(＞40)，穩(wěn)定性弱，沒有跨膜信號，蛋白氨基酸序列不存在卷曲螺旋。StR2R3–MYB1 蛋白的疏水區(qū)域分布比較均衡，疏水性氨基酸多于親水性氨基酸，但其N 端具有明顯的親水區(qū)域（圖4）。

2.3 StR2R3–MYB1 核酸序列的結構

StR2R3–MYB1 基因位于馬鈴薯10 號染色體的51 749 197～51 750 359 bp，基因全長1 163 bp，由3個外顯子和2個內含子構成。CDS 序列中含有2個DNA–結合結構域，分別位于第21～71 aa 和第74～122 aa，屬于R2R3 類MYB 轉錄因子，將該基因命名為StR2R3– MYB1。對啟動子順式作用元件[8]進行分析的結果表明，StR2R3–MYB1 基因的調控區(qū)域存在大量與光相關的元件，如G–Box、Box4、BoxⅡ、I–Box、MNF1 等，還存在脫落酸誘導相關元件ABRE、生物鐘相關元件Circadian 和茉莉酸甲酯相關元件CGTCA– motif (圖5)。

圖4 StR2R3-MYB1 蛋白疏水性/親水性預測結果 Fig. 4 StR2R3-MYB1 protein hydrophobic/hydrophilic prediction

2.4 StR2R3–MYB1 蛋白的亞細胞定位分析

將構建好的StR2R3–MYB1 與YFP 融合蛋白表達載體轉化擬南芥原生質體，在激光共聚焦掃描顯微鏡下觀察基因的亞細胞定位情況，可見StR2R3–MYB1∷YFP 融合蛋白定位在細胞核內 (圖6)。

2.5 StR2R3–MYB1 基因表達模式的分析

實時熒光定量PCR 分析結果顯示，StR2R3– MYB1 在根、莖、葉的相對表達量依次升高(表3)。

表3 StR2R3–MYB1 基因在不同組織部位的相對表達量 Table 3 Gene expression pattern analysis for StR2R3–MYB1

2.6 光調節(jié)StR2R3–MYB1 基因的表達

黑暗培養(yǎng)12 d 的馬鈴薯白化苗轉移至白光下生長0、0.5、2、4、8、12、24 h 的StR2R3–MYB1基因表達量隨光處理時間的延長呈上升趨勢(圖7)， StR2R3–MYB1 的表達受到了白光的誘導。分別用紅、藍和遠紅光處理白化苗后，3 種單色光都能誘導StR2R3–MYB1 的表達，其中遠紅光、藍光誘導4 h 時，StR2R3–MYB1 的表達量出現(xiàn)第1個峰值，隨后下降，在12 h 后回升。紅光對StR2R3–MYB1 的誘導作用在8 h 后開始顯現(xiàn)(圖7)。

圖7 光調節(jié)StR2R3–MYB1 基因的表達 Fig. 7 Expression of light regulationgene of StR2R3–MYB1

3 結論與討論

在植物響應逆境脅迫的信號傳導系統(tǒng)中，轉錄因子起著重要的作用。序列分析結果表明，馬鈴薯StR2R3–MYB1 與其他MYB家族轉錄因子在DNA 結合域具有高度保守性，屬于典型的R2R3 類型MYB 轉錄因子。核苷酸序列和氨基酸序列分析結果表明，StR2R3–MYB1 的調控區(qū)域存在大量與光相關的元件，如G–Box、Box4、BoxⅡ、I–Box、MNF1等，還存在脫落酸誘導相關元件ABRE、生物鐘相關元件Circadian、茉莉酸甲酯相關元件CGTCA。StR2R3–MYB1 的轉錄很可能受光信號及環(huán)境脅迫的誘導。StR2R3–MYB1 氨基酸序列上沒有跨膜信號，不存在卷曲螺旋，具有明顯的親水區(qū)域，氨基酸殘基與水分子的相互作用較強，表明該蛋白的穩(wěn)定性較低。亞細胞定位分析發(fā)現(xiàn)StR2R3–MYB1 蛋白在細胞核中表達，這正符合轉錄因子的亞細胞定位特征。系統(tǒng)進化分析結果表明，StR2R3–MYB1 與擬南芥AtMTB113、AtMTB90 和茄子SmMYB 及番茄SlCMYB1 同屬于1個亞族。該亞族MYB 轉錄因子的功能多數(shù)與花青素合成相關，如AtMTB113 與花青素積累有關[8–9]，而AtMTB90 與色素合成相關[10]，SmMYB 與花青素合成相關[11]，SlCMYB1 則與冷誘導、低溫脅迫相關[12]。筆者的前期研究[7]中運用多因子及其互作逐步回歸法建立光照度、溫度和處理時間影響基因表達的模式，發(fā)現(xiàn)環(huán)境溫度與StR2R3–MYB1 和StPAL、StDFR 基因的表達量呈負相關，環(huán)境光照度與StR2R3–MYB1 和StCHS、StDFR的表達呈正相關。馬鈴薯StR2R3–MYB1 轉錄調節(jié)因子對StCHS 基因的表達具有明顯正向影響。本研究中StR2R3–MYB1 蛋白結構分析、蛋白功能預測、基因組織表達特性分析及光處理后基因表達特性分析等研究結果，使前期研究中的推斷StR2R3–MYB1 與受光和低溫調節(jié)的馬鈴薯花青素合成相關[7]得到了進一步的印證。

[1] Thompson M A，Ramsay R G．Myb：An old onco protein with new roles[J]．Bioessays，1995：17(4)：341–350．

[2] Ogata K，Hojo H，Aimoto S，et al．Solution structure of a DNA–binding unit of Myb：A helix-turn-helix-related motif with conserved tryptophans forming a hydrophobic core[J]．Proc Natl Acad Sci USA，1992，89(14)：6428– 6432．

[3] Paz-Ares J，Ghosal D，Wienand U，et al．The regulatory c1 locus of zea mays encodes a protein with homology to myb proto-oncogene products and with structural similarities to transcriptional activators[J]．EMBO J，1987，6(12)：3553–3558．

[4] Rabinowicz P D，Braun E L，Wolfe A D，et al．Maize R2R3-Myb genes ： Sequence analysis reveals amplification in the higher plants[J]．Genetics，1999，153(1)：427–444．

[5] Yanhui C，Xiaoyuan Y，Kun H，et al．The MYB transcription factor superfamily of Arabidopsis ：Expression analysis and phylogenetic comparison with the rice MYB family[J]．Plant Mol Biol，2006，60(1)：107–124．

[6] Knight H，Trewavas A J，Knight M R．Calcium signalling in Arabidopsis thaliana responding to drought and salinity[J]．Plant J，1997，12(5)：1067–1078．

[7] Qin Y Z，Tai G，Xie K Y，et al．Ambient light atersgene expression pattern of enzymes and transcription factor involved in phenylpropanoid metabolic pathway in potato under chilling stress[J] ． Agricultural Science & Technology，2014，15(11)：1899–1904．

[8] Schwob E，Bohm T，Mendenhall M D，et al．The B–type cyclin kinase inhibitor p40SIC1 controls the G1 to S transition in S．cerevisiae [J]．Cell,1994，79(2)：233–244．

[9] Luo Q J，Mittal A，Jia F，et al．An autoregulatory feedback loop involving PAP1 and TAS4 in response to sugars in Arabidopsis[J]．Plant Mol Biol，2012，80(1)：117–129．

[10] Velten J，Cakir C，Cazzonelli C I．A spontaneous dominant-negative mutation within a 35S：AtMYB90 transgene inhibits flower pigment production in tobacco[J]．PLoS One，2010，5(3)：e9917．

[11] 邵文婷，劉楊，韓洪強，等．茄子花青素合成相關基因SmMYB 的克隆與表達分析[J]．園藝學報，2013，40(3)：467–478．

[12] 張欣，程治軍，林啟冰，等．番茄冷誘導基因SlCMYB1的克隆及其在水稻中異源表達研究[J]．作物學報，2011(4)：587–594．