李 丹,趙存鵬,劉素恩,王凱輝,張曉慧,趙麗英,郭寶生,耿軍義

(河北省農(nóng)林科學(xué)院 棉花研究所,農(nóng)業(yè)農(nóng)村部黃淮海半干旱區(qū)棉花生物學(xué)與遺傳育種重點(diǎn)實(shí)驗(yàn)室,國家棉花改良中心河北分中心,河北 石家莊 050051)

2-酮戊二酸依賴型氧化酶(2-oxoglutarate-dependent dioxygenase,2OGD)基因超級(jí)家族是參與植物代謝的最大基因家族之一,主要參與天然產(chǎn)物的氧化過程,尤其是萜類化合物[1]。2OGD通常使用亞鐵Fe(Ⅱ)作為輔因子,2OG和O2作為輔助底物,產(chǎn)生副產(chǎn)物CO2和琥珀酸鹽(R+2OG+O2→ R-OH+琥珀酸鹽+CO2)[2]。2OGD的這種催化反應(yīng)與多種生物過程有關(guān),包括脯氨酸羥基化、核酸和組蛋白的氧化去甲基化以及許多代謝物的生物合成[2]。

Kawai等[2]對(duì)6種植物中的479個(gè)2OGD進(jìn)行了系統(tǒng)發(fā)育分類。根據(jù)氨基酸序列的相似性,將其分為3類:DOXA、DOXB和DOXC。DOXA類包括大腸桿菌AlkB的植物同源物,它是參與烷基化核酸和組蛋白氧化去甲基化的2OGD的原型[2-3]。DOXB類在所有植物類群中都是保守的,并參與細(xì)胞壁蛋白質(zhì)合成中的脯氨酸4-羥基化[2,4]。DOXC類參與各種植物化學(xué)物質(zhì)的代謝,包括赤霉素、花青素、類黃酮、生長素、硫代葡萄糖苷等[1-2]。

硫代葡萄糖苷(Glucosinolates,GSL)是主要存在于十字花科植物中的一種次生代謝物質(zhì)。近年研究表明,GSL在十字花科植物應(yīng)對(duì)生物脅迫和非生物脅迫中具有重要作用[5]。GSL可分為3類:脂肪族GSL、吲哚族GSL和芳香族GSL[6]。AOPs是脂肪族GSL的合成基因[7],屬于DOXC類[1],其中AOP2及其同源基因已被證明參與烯基GSL(脂肪族GSL中的一種)生物合成[8-12]。Kliebenstein等[8,10]發(fā)現(xiàn),功能性AOP2的表達(dá)與烯基GSL的積累之間存在完全相關(guān)。

干旱和鹽漬是影響植物生長和農(nóng)作物產(chǎn)量的重要環(huán)境因子[13]。本研究前期測(cè)定了耐旱耐鹽品種冀2658在干旱和鹽脅迫下12,24 h的根蛋白質(zhì)組,發(fā)現(xiàn)GhAOP2-like(2-oxoglutarate-dependent dioxygenase AOP2-like)蛋白在2種處理下均上調(diào)表達(dá),說明GhAOP2-like在陸地棉響應(yīng)干旱和鹽脅迫時(shí)具有重要功能。但是,關(guān)于AOP2的功能研究主要集中在十字花科植物中,錦葵科中尚未見報(bào)道。所以研究陸地棉中GhAOP2-like響應(yīng)干旱和鹽脅迫的機(jī)制具有重要意義,同時(shí)也是對(duì)AOP基因家族功能的完善。為了解析GhAOP2-like響應(yīng)脅迫的機(jī)理,對(duì)GhAOP2-like進(jìn)行了克隆和初步功能研究。

1 材料和方法

1.1 試驗(yàn)材料

試驗(yàn)材料為河北省農(nóng)林科學(xué)院棉花研究所培育的耐旱耐鹽棉花(Gossypiumhirsutum)品種冀2658。

1.2 基因克隆

利用COTTONOMICS(http://cotton.zju.edu.cn/)網(wǎng)站上在線軟件將包含GhAOP2-like序列前后各2 000 bp的一段序列提取出,用軟件Premier 5在GhAOP2-like序列前后設(shè)計(jì)引物,以冀2658的cDNA為模板進(jìn)行基因克隆。引物序列:Forward primer:ATGGGTGTCCAAGCTGAAATT;Reverse primer:TTAATTGGAAATTGGTGGTGAGA。

1.3 生物信息學(xué)分析

GO富集分析用Gene Ontology數(shù)據(jù)庫(http://www.geneontology.org/)。用ExPASy-ProtParam(https://web.expasy.org/protparam/)對(duì)蛋白進(jìn)行理化性質(zhì)分析;用SignalP 5.0 Server(http://www.cbs.dtu.dk/services/SignalP)和TMHMM Server v.2.0(http://www.cbs.dtu.dk/services/TMHMM)分別對(duì)蛋白的信號(hào)肽及跨膜結(jié)構(gòu)進(jìn)行預(yù)測(cè);用SPOMA(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=np%20sa_sopma.html)和Phyre2(http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index)分別對(duì)蛋白進(jìn)行二級(jí)結(jié)構(gòu)和3D模型預(yù)測(cè)。利用PSI(http://bis.zju.edu.cn/psi/)進(jìn)行GhAOP2-like的亞細(xì)胞定位預(yù)測(cè)。利用NCBI上PSI-BLAST Iteration 1程序(https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome)對(duì)序列進(jìn)行Blast比對(duì),用Fast Minimum Evolution方法進(jìn)行聚類分析。

1.4 試驗(yàn)材料的干旱脅迫、鹽脅迫和激素處理

將冀2658種子種植在滅菌的沙子中,待子葉展開,將幼苗輕輕取出,沖洗干凈根部后,將植株移到水培裝置中,培養(yǎng)至二葉一心時(shí)期,取植株的根、莖和葉片進(jìn)行GhAOP2-like組織特異性表達(dá)分析,每個(gè)組織設(shè)置3個(gè)生物學(xué)重復(fù);剩余植株進(jìn)行干旱脅迫、鹽脅迫和激素處理。模擬干旱脅迫用25.16%的PEG-6000溶液,鹽脅迫用150 mmol/L NaCl溶液。為了避免2種脅迫因溶液滲透壓不同引起的試驗(yàn)誤差,本試驗(yàn)所用溶液的濃度在室溫25 ℃下的滲透壓基本相等,PEG-6000溶液的滲透壓計(jì)算見參考文獻(xiàn)[14]。激素處理濃度分別為100 μmol/L ABA、100 μmol/L MeJA、150 μmol/L GA3。取處理后0,6,12,24,48 h的幼苗根部提取RNA進(jìn)行qRT-PCR驗(yàn)證GhAOP2-like表達(dá)量的變化,每個(gè)時(shí)間點(diǎn)取3個(gè)生物學(xué)重復(fù)。

1.5 RNA提取和qRT-PCR驗(yàn)證

用BioFlux Plant Total RNA Extraction Kit(DNA-free)試劑盒(購自BioFlux)按照說明書制備總RNA。第一鏈cDNA合成和qRT-PCR分別采用HiScript Ⅱ Q RT SuperMix for qPCR(+gDNA wiper)試劑盒和ChamQ Universal SYBR qPCR Master Mix(均購自諾唯贊)進(jìn)行,在BIO-RAD CFX96實(shí)時(shí)系統(tǒng)中進(jìn)行qRT-PCR?；虻南鄬?duì)表達(dá)量采用2-ΔΔCt方法分析,試驗(yàn)進(jìn)行3次技術(shù)重復(fù),用軟件SPSS 進(jìn)行顯著性差異分析,以棉花Actin1作為參考基因。定量引物序列為:GhAOP2-like:Forward primer:TGGAAAGGAAACAAAAGAACGC;Reverse primer:CCAAGTCCAAAGCCTTCATACAAC。Actin1:Forward primer:ATCCTCCGTCTTGACCTTG;Reverse primer:TGTCCGTCAGGCAACTCAT。

2 結(jié)果與分析

2.1 GhAOP2-like的篩選

本研究前期利用TMT(Tandem Mass Tags)技術(shù)對(duì)干旱和鹽脅迫12,24 h的根部蛋白質(zhì)組進(jìn)行了分析。結(jié)果發(fā)現(xiàn),GhAOP2-like在2種脅迫下均呈現(xiàn)上調(diào)表達(dá)的趨勢(shì),且差異表達(dá)倍數(shù)較高,表達(dá)水平穩(wěn)定(表1)。GO富集分析顯示,GhAOP2-like參與氧化還原過程。此外,利用Cotton Omics Database(http://cotton.zju.edu.cn/download.html)中TM-1_V2.1-CDS數(shù)據(jù)進(jìn)行Blast比對(duì)獲得TM-1中GhAOP2-like基因序列,TM-1_v2.1_final.annotation中數(shù)據(jù)顯示GhAOP2-like參與萜類合成。

表1 GhAOP2-like在干旱和鹽脅迫下的差異表達(dá)倍數(shù)Tab.1 The FC value of GhAOP2-like under drought and salt stresses

2.2 GhAOP2-like的克隆

利用同源克隆法從冀2658中克隆得到GhAOP2-like的CDS序列(圖1-A),序列全長972 bp,編碼323個(gè)氨基酸。選取TM-1_V2.1-CDS數(shù)據(jù)庫進(jìn)行Blast比對(duì),發(fā)現(xiàn)GhAOP2-like位于D13染色體上,并與G.raimondii-CDS數(shù)據(jù)庫中GhAOP2-like序列完全一致。

A.GhAOP2-like基因PCR擴(kuò)增產(chǎn)物的電泳;B.GhAOP2-like蛋白三維模型。M.DNA分子量標(biāo)記;G.GhAOP2-like 的CDS條帶。A.Electropherogram of PCR amplified product of GhAOP2-like;B.3D model of GhAOP2-like protein.M.DNA Marker;G.CDS sequence of GhAOP2-like.

2.3 GhAOP2-like蛋白的生物信息學(xué)分析

利用ExPASy-ProtParam分析得到GhAOP2-like蛋白的分子式為C1654H2525N429O478S19;分子質(zhì)量為36.676 97 ku;理論等電點(diǎn)為5.20;總親水性平均數(shù)為-0.223;不穩(wěn)定指數(shù)為42.71。GhAOP2-like蛋白包含323個(gè)氨基酸,其中亮氨酸占8.4%,谷氨酸占8.4%,絲氨酸占7.1%,甘氨酸占6.8%。

預(yù)測(cè)顯示,GhAOP2-like不含信號(hào)肽,不含跨膜結(jié)構(gòu)域。二級(jí)結(jié)構(gòu)中,α-螺旋占8.55%,延伸主鏈占32.02%,無規(guī)則卷曲占53.29%,β-折疊占6.14%。3D模型預(yù)測(cè)如圖1-B,模型的置信度達(dá)到100%,氨基酸序列的覆蓋度達(dá)到93%,因此,可以認(rèn)為該蛋白結(jié)構(gòu)的預(yù)測(cè)是可靠的。PSI結(jié)果顯示GhAOP2-like定位在細(xì)胞質(zhì)中。

2.4 GhAOP2-like氨基酸同源序列的聚類分析

利用NCBI上CPSI-BLAST Iteration 1程序?qū)hAOP2-like的氨基酸序列進(jìn)行Blast比對(duì),用Fast Minimum Evolution方法對(duì)比對(duì)到的部分序列進(jìn)行聚類分析。從圖2看出,棉屬與其他物種的AOP序列明顯被分成2組,但是與同樣來自錦葵科的木槿(Hibiscussyriacus)的AOP序列關(guān)系最近。來源于棉屬的AOP1和AOP2序列被明顯分開,但序列之間差異較小。

圖2 GhAOP2-like蛋白與同源序列的聚類分析Fig.2 The cluster analysis of GhAOP2-like protein sequence and its homologous sequences

2.5 GhAOP2-like組織特異性表達(dá)

取植株的根、莖、葉片和發(fā)育中的種子進(jìn)行GhAOP2-like組織特異性表達(dá)分析,如圖3-A所示,GhAOP2-like在根中表達(dá)量最高,推測(cè)GhAOP2-like可能對(duì)根部發(fā)育具有重要作用。為了進(jìn)一步驗(yàn)證GhAOP2-like在不同組織中的表達(dá)量,在Cotton Omics Database搜索了TM-1和海7124中GhAOP2-like的組織表達(dá),發(fā)現(xiàn)GhAOP2-like同樣在根中表達(dá)最高(表2)[15]。

A.GhAOP2-like組織特異性表達(dá)分析;B.干旱處理;C.鹽處理;D.ABA處理;E.GA3處理;F.MeJA處理。不同小寫字母代表不同組織或不同處理時(shí)間的差異顯著性(P<0.05)。A.The tissue-specific expression analysis of GhAOP2-like;B.Drought treatment;C.Salt treatment;D.ABA treatment;E.GA3 treatment;F.MeJA treatment.Different lowercase letters indicate signnificant differences among different tissues or different treatment times(P<0.05).

表2 TM-1和海7124中GhAOP2-like在不同組織中的相對(duì)表達(dá)Tab.2 Relative expression level of GhAOP2-like from TM-1 and Hai 7124 in different tissues

2.6 干旱、鹽脅迫以及激素處理下GhAOP2-like表達(dá)量變化

用qRT-PCR對(duì)干旱、鹽脅迫下冀2658根部GhAOP2-like的表達(dá)量進(jìn)行驗(yàn)證,從圖3-B—F中看出,與對(duì)照相比,GhAOP2-like在干旱和鹽脅迫后表達(dá)量升高,這與蛋白質(zhì)組測(cè)序結(jié)果一致。同樣,在3種激素處理后,GhAOP2-like表達(dá)量也呈現(xiàn)不同程度地上升趨勢(shì),但在MeJA處理下表達(dá)量差異倍數(shù)最高,這可能是植物內(nèi)不同激素間相互調(diào)控導(dǎo)致的。

3 結(jié)論與討論

AOP基因家族參與十字花科植物中脂肪族GSL的合成。不少研究發(fā)現(xiàn),在環(huán)境脅迫下,植物中GSL的積累可參與滲透調(diào)節(jié)、提高水分吸收和運(yùn)輸、增強(qiáng)植物的抗氧化能力[7]。因此,AOP基因?qū)χ参锬婢诚碌纳L發(fā)育具有非常重要的作用。

在擬南芥中,Kliebenstein等[7-8]鑒定出了AOP1、AOP2和AOP33個(gè)基因。其中AOP2能催化產(chǎn)生烯烴基硫苷。另外,Neal等[10]發(fā)現(xiàn),AOP2至少以3個(gè)等位基因的形式存在于不同的擬南芥種質(zhì)資源中,并證明其中一個(gè)等位基因AOP2-2在功能上能夠形成烯基形式硫苷。在甘藍(lán)和芥藍(lán)中,分別克隆了AOP2的同源基因BoGSL-ALK[9-10]和BoaAOP-like[12],它們都能催化產(chǎn)生脂肪族GSL。目前,關(guān)于AOP或AOP-like的功能研究主要集中于十字花科植物中,而且已發(fā)現(xiàn)的含有GSL的16個(gè)科中不包含錦葵科[16]。所以,陸地棉中AOP2-like可能像TM-1_v2.1_final.annotation中注釋的一樣參與其他萜類的合成。

研究發(fā)現(xiàn),植物中GSL代謝應(yīng)答模式根據(jù)脅迫的時(shí)間、強(qiáng)度、植物所處生長時(shí)期而具有多樣性[7]。Salehin等[17]在擬南芥中過表達(dá)AOP2可提高脂肪族GSL含量,并且GSL含量增加可能通過產(chǎn)生ROS促進(jìn)氣孔關(guān)閉從而提高植物耐旱性。鹽脅迫下,脂肪族GSL可作為參與植物滲透調(diào)節(jié)作用的元素影響水平衡[18]。Martínez-Ballesta等[18]研究發(fā)現(xiàn),脂肪族GSL可能與鹽脅迫下根系增殖有關(guān),并且脂肪族GSL含量的變化可以影響擬南芥種子的萌發(fā),改變水通道蛋白的豐度。但在高濃度鹽溶液下,GSL降解酶活性會(huì)提高,GSL的降解產(chǎn)物具有解毒的作用[19]。在芥藍(lán)中,Zheng等[12]克隆了一個(gè)BoaAOP-like基因,能夠調(diào)控脂肪族GSL的合成。BoaAOP-like定位于細(xì)胞質(zhì),與GhAOP2-like亞細(xì)胞定位預(yù)測(cè)結(jié)果一致。BoaAOP-like在根中表達(dá)量最高,并且在BoaAOP-like瞬時(shí)過度表達(dá)的植物中能夠提高AOP2.1的表達(dá)水平。在本研究前期測(cè)定的蛋白質(zhì)組中,鹽脅迫下水通道蛋白豐度確實(shí)發(fā)生改變,其中9個(gè)上調(diào),1個(gè)下調(diào),但這是否與GhAOP2-like的調(diào)控有關(guān)還需進(jìn)一步驗(yàn)證。而在干旱脅迫下,只有1個(gè)水通道蛋白下調(diào),說明GhAOP2-like調(diào)控干旱和鹽脅迫的機(jī)制不同。此外,GhAOP2-like在TM-1和海7124根部的高表達(dá)說明GhAOP2-like對(duì)根部的生長發(fā)育確實(shí)具有重要作用,可能調(diào)控脂肪族GSL類似物或萜類的合成,參與植物中滲透調(diào)節(jié)以維持水分平衡促進(jìn)根系生長,并且在干旱脅迫下參與調(diào)節(jié)ROS的生成。

MYB家族是植物體內(nèi)最大的一類轉(zhuǎn)錄因子,參與多種代謝途徑的調(diào)控[20]。其中MYB28、MYB29、MYB76正調(diào)控脂肪族GSL合成[21]。在甘藍(lán)中,BoAOP2的表達(dá)與BoMYB28、BoMYB29表達(dá)呈正相關(guān)[22]。擬南芥Col-0(擬南芥Col-0天然缺失AtAOP2基因)中過表達(dá)甘藍(lán)BoAOP2導(dǎo)致AtMYB28、AtMYB29的表達(dá)量和硫苷都有所增加。另外將BoAOP2在MYB28/29雙突變體中過表達(dá),與對(duì)照相比,硫苷含量沒有明顯變化,說明BoAOP2過表達(dá)無法在雙突變體中實(shí)現(xiàn)硫苷積累[11]。此外,BoAOP2還可上調(diào)AtMYC2的表達(dá)[11]。研究表明,MYC2除了參與硫苷的合成,也是茉莉酸(Jasmonic acid,JA)信號(hào)通路的主要調(diào)控因子,MYC2參與JA調(diào)控的植物發(fā)育、側(cè)根和不定根形成[23],這可以解釋Martínez-Ballesta等[18]發(fā)現(xiàn)的脂肪族GSL與根系增殖相關(guān)。而且,MYC2也調(diào)節(jié)JA信號(hào)通路與其他植物激素(如脫落酸、水楊酸、赤霉素和生長素)之間的相互作用。由此可見,AOP2參與多個(gè)代謝網(wǎng)絡(luò),可能在植物防御中起關(guān)鍵作用。本研究發(fā)現(xiàn),GhAOP2-like的表達(dá)確實(shí)受到多種激素的調(diào)控,并且在MeJA處理下表達(dá)變化最強(qiáng)烈。所以,和AOP2一樣,GhAOP2-like也可能受到轉(zhuǎn)錄因子MYC2、MYB28、MYB29的調(diào)控,參與多種激素之間的協(xié)同作用,從而提高植物抗旱耐鹽的能力,但還需進(jìn)一步論證。