小黑麥SWEET家族基因鑒定及其在不同逆境下表達模式分析

李 根, 牛奎舉

(甘肅農業大學草業學院/草業生態系統教育部重點實驗室/甘肅省草業工程實驗室/中-美草地畜牧可持續發展研究中心,甘肅 蘭州 730070)

在植物體可以利用的養分中,糖是生物體主要的碳源和能量來源[1],同時也是一種重要的信號分子。在高等植物的葉肉細胞中,葉綠體可以通過光合作用將二氧化碳轉化為有機碳固定,使光合作用產生糖[2-4]。然而,糖不能通過植物細胞膜系統獨立運輸,需要適當的糖轉運體的協助[2],目前主要發現有三類蛋白,即單糖轉運體(Monosaccharide transporters,MSTs)[5],蔗糖轉運體(Sucrose transporters,SUTs)[6-7]和糖外排轉運蛋白(Sugars will eventually be exported transporters,SWEETs)[8]。在過去的20年里,糖轉運蛋白的前兩個基因家族MSTs和SUTs在高等植物中得到了廣泛的研究[9],而SWEETs基因家族是近些年來被確定為糖外排體,其成員在已糖或蔗糖跨細胞膜轉運過程中發揮了重要作用[10],但是人們對其功能了解尚不完全。

SWEETs廣泛存在于真核與原核生物,是一類不依賴于環境pH值,僅利用細胞內外的濃度差來實現糖類的雙向可逆轉運功能的糖轉運蛋白家族[10]。真核SWEETs蛋白通常由七個α-螺旋跨膜結構域(7α-helical transmembrane domain,7-TMs)組成,這些結構域由兩個3-TM結構域(包含兩個保守的MtN3/saliva)串聯重復組成,由一個保守程度較低的跨膜螺旋連接,并使其分開,該跨膜螺旋結構通常被稱為3-1-3TM SWEET結構[11]。研究發現,SWEETs糖轉運蛋白家族成員能夠參與調控植物對非生物脅迫的響應過程。擬南芥(Arabidopsisthaliana)AtSWEET15在葉片衰老和滲透脅迫(包括高鹽、寒冷和干旱)會受到脫落酸(Abscisic acid,ABA)信號誘導,過表達AtSWEET15會使植株葉片衰老加速,對高鹽脅迫敏感,而缺乏AtSWEET15的突變株對高鹽敏感性降低[12]。在寒冷脅迫下,在擬南芥中過表達茶樹[Camelliasinensis(L.)O. Ktze.]CsSWEET1a和CsSWEET17,轉基因株系與野生型相比,相對電導率顯著降低,蔗糖轉化酶AtCWINVs和AtVACINVs基因的表達被抑制,從而有效提高擬南芥抗凍性[13]。擬南芥AtSWEET11和AtSWEET12與寒冷及干旱脅迫相應有著密切聯系,AtSWEET11,AtSWEET12雙突變體比野生型和單突變體表現出更強的抗凍性[14]。

小黑麥(×Triticosecale)是一種由人類創制的自花授粉的谷類作物,是由小麥(Triticumspp.)和黑麥(Secalecereale)雜交而成的新物種,在世界各地廣泛種植。它結合了小麥優越的農業形態、最終利用品質特征和黑麥的強適應性、高活力、對非生物和生物脅迫的優質抗性,作為優良的一年生糧飼兼用作物在世界范圍內廣泛種植[15-17],為緩解牧區草畜矛盾,保護生態環境,促進草地畜牧業可持續發展發揮了重要作用。但其在生長發育過程中也同樣易受到環境脅迫,嚴重時會影響牧草的飼用品質及產量。因此,培育高產抗逆的優質小黑麥品種十分必要。隨著現代分子生物學和生物信息學技術的發展,草類植物響應環境脅迫的機理研究已從最初的形態、生理水平逐漸發展到轉錄、翻譯及大分子物質網絡互作調控等分子機制層面[18]。截至目前,很多研究表明SWEETs糖轉運蛋白與生物及非生物脅迫響應聯系密切,但小黑麥SWEETs家族參與響應非生物脅迫的研究較少。因此,本研究對小黑麥SWEETs基因家族進行了鑒定,并系統分析了SWEET基因的結構特點、保守基序、蛋白質結構等基本特征,在此基礎上進一步分析了TwSWEETs在近源及模式植物間的進化關系,同時對其在小黑麥葉片和根部在不同逆境脅迫下的基因表達情況進行了深入分析。本研究為進一步了解SWEETs基因在調節植物發育過程中糖的積累和分配及它們在非生物脅迫中發揮的作用奠定理論基礎,進而有助于優良基因資源的挖掘和新品種培育。

1 材料與方法

1.1 試驗材料

選用‘石大一號’品種小黑麥,由本實驗室保存并提供。

1.2 試驗方法

1.2.1小黑麥TwSWEETs基因家族的鑒定 擬南芥、水稻(OryzasativaL.)、小麥(TriticumaestivumL.)和玉米(ZeamaysL.)的SWEETs基因蛋白序列可直接從NCBI數據庫中下載(https://www.ncbi.nlm.nih.gov/),并根據擬南芥和水稻的SWEET蛋白序列作為鑒定小黑麥TwSWEETs的請求序列,以本課題組已有的多個小黑麥轉錄組數據作為本地數據庫。使用BioEdit Sequence Alignment Editor軟件進行搜索比對,Program選擇‘tblastn’,Expectation Value:E-20,搜索比對小黑麥的SWEET基因家族蛋白序列。

1.2.2小黑麥TwSWEETs蛋白序列理化性質分析 使用在線ExPAsy-ProtParam工具(https://web.expasy.org/protparam/)分析小黑麥TwSWEETs蛋白序列的理化性質,包括氨基酸數量、分子量、理論等電點、不穩定性指數、脂肪族指數和總平均親水性[19];使用TMHMM在線工具預測跨膜螺旋結構(TMHMM 2.0-DTU Health Tech-Bioinformatic Services)[20];使用ProtComp v. 9.0在線工具[Plant-mPLoc server(sjtu.edu.cn)]初步預測亞細胞定位[21];使用SOPMA在線工具對蛋白質的二級結構進行預測[NPS@:SOPMA secondary structure prediction (ibcp.fr)];使用SWEISS MODEL在線工具對蛋白質的三級結構進行預測(https://swissmodel.expasy.org/)[22-26]。

1.2.3小黑麥TwSWEETs蛋白序列系統進化樹及保守基序分析 利用在線工具MEME(http://meme-suite.org/)對小黑麥的SWEET糖轉運蛋白的氨基酸序列進行保守基序分析[27]。利用MEGA 11.0軟件中的最大似然法構建16條小黑麥、59條小麥、17條擬南芥、21條水稻和29條玉米的基序系統進化樹。

1.2.4小黑麥幼苗生長條件及非生物脅迫處理 采用土培法育苗,將小黑麥‘石大一號’種子用無菌水過夜浸泡,沖洗數次后進行消毒處理,然后種植于育苗缽中,放置在甘肅農業大學草業學院生長室中。晝夜溫度為25℃/20℃,每天光周期循環14 h,相對濕度為50%。待生長14天后,選擇長勢均一的幼苗進行非生物脅迫處理。試驗設置3組處理,每組處理4個重復。對照處理(CK)為正常純凈水澆灌;干旱處理為20% PEG6000澆灌種植土壤,處理后放置于生長室中;低溫處理設置溫度為4℃,放置于低溫光照培養箱。分別于處理3 h,6 h,12 h和24 h采集葉片和根系并放入液氮中快速冷凍,保存于-80℃冰箱。

1.2.5RNA提取和實時定量PCR(qRT-PCR)分析 RNA提取使用天根生物科技有限公司(中國)的總RNA提取試劑盒(離心柱型),貨號#DP419。然后使用寶生物(Takara)大連有限公司的PrimerScriptTMRT reagent Kit with gDNA Eraser(Perfect Real Time)反轉錄試劑盒合成cDNA,貨號R047A。具體試驗步驟按照制造商的說明進行。

使用在線工具Primer BLAST(https://www.ncbi.nlm.nih.gov/tools/primer-blast/)設計了小黑麥TwSWEETs基因的特異引物,由西安擎科生物科技有限公司合成。qRT-PCR使用天根生化科技(北京)有限公司的SuperReal RreMix Plus(SYBR Green)熒光定量試劑盒,貨號FP206-02。20 μL反應總體積包括5 μL稀釋模板(20 μL cDNA用180 μL ddH2O稀釋20倍)、3 μL ddH2O,1 μL上游和1 μL下游引物(10 μmol·L-1)和10 μL 2×SuperReal PreMix Plus(SYBR Green)。反應過程在LightCycler?96儀器(中國上海Roche)上進行。反應程序為:95℃預變性15 min,PCR反應為95℃變性10 s和58℃退火30 s,40次循環,熔解曲線分析保留默認,為95℃持續1 s,65℃持續15 s,95℃持續1 s。每個處理設置三次生物學重復,每個生物學重復設置4次技術重復,去掉極端數值。使用2-ΔΔCt方法計算相對基因表達值[28]。選用Actin(ACT)作內參基因,除TwSWEET6b-1和TwSWEET6b-2沒有合適的引物以外,其余14條小黑麥TwSWEETs基因特異性引物如表1所示。

表1 qRT-PCR引物序列表Table 1 Primers used in this study for qPCR

1.2.6數據分析 采用SPSS 24.0 Windows版(SPSS Inc.,USA)進行數據分析,用平均值和標準誤表示測定結果,對照組與不同時間的處理組采用t檢驗,采用GraphPad Prism 8.0.2軟件繪制柱狀圖。

2 結果與分析

2.1 小黑麥TwSWEET基因家族的鑒定及理化性質分析

通過對比分析在小黑麥本地轉錄組數據庫中獲得16條序列,并將其命名為TwSWEET1a,TwSWEET1b,TwSWEET2a,TwSWEET2b,TwSWEET3a,TwSWEET4,TwSWEET6a,TwSWEET6b-1,TwSWEET6b-2,TwSWEET11,TwSWEET12,TwSWEET13,TwSWEET14,TwSWEET15,TwSWEET16和TwSWEET17。

如表2所示,16條小黑麥TwSWEETs蛋白序列編碼的氨基酸數目在231(TwSWEET2b)到311(TwSWEET12)之間,平均約270個氨基酸,分子質量在25.428 KD到33.913 KD之間。除TwSWEET1a/12/17的理論等電點(Theoretical isoelectric point,pI)小于7外,其余蛋白的pI都大于7,屬于堿性蛋白質。根據蛋白質不穩定系數是否大于40可將其分為穩定蛋白(小于40)和不穩定蛋白(大于40)。分析發現,16條小黑麥TwSWEETs蛋白中有10條屬于穩定蛋白,6種屬于不穩定蛋白。另外,16條小黑麥TwSWEETs蛋白脂肪指數在108.86到123.40之間,親水性(Grand average of hydropathicity,GRAVY)均大于零,說明它們都是疏水性蛋白。對TwSWEETs蛋白的跨膜螺旋結構(Transmembrane helices,TMH)預測發現,小黑麥TwSWEETs家族包含6～7個TMHs,約70%的成員擁有7個TMHs;此外,對這些蛋白的亞細胞定位初步預測結果發現,它們多數都定位于質膜(Plasma membrane,PM)上,少數SWEETs基因存在多種定位結果。

表2 小黑麥TwSWEETs蛋白序列理化性質分析Table 2 Physicochemical properties analysis of TwSWEETs protein sequence in triticale

2.2 小黑麥TwSWEETs編碼的蛋白質二級和三級結構預測

如表3和圖1所示,16條TwSWEETs編碼的蛋白質二級結構均由α-螺旋結構、無規則卷曲結構、β-轉角結構和延伸鏈組成,但各結構所占比例有所差異。在所有TwSWEETs蛋白質中,絕大多數蛋白質α-螺旋結構相比其他結構所占比例最高,除TwSWEET3a/4/6a/6b-1/11/12/13外,其余TwSWEETs的α-螺旋結構所占比例均高于40%。無規則卷曲結構所占比例在27.71%～40.84%之間,平均為34.45%,β-轉角結構在四種結構中所占比例最低,其中,TwSWEET11的β-轉角結構占比為2.64%,在所有小黑麥TwSWEETs家族成員中占比最低,TwSWEET16的β-轉角結構占比最高,為5.02%。延伸鏈在16條TwSWEETs家族成員中的占比有所浮動,在12.37%～25.91%之間。綜合分析發現,α-螺旋結構和無規則卷曲結構是小黑麥TwSWEETs家族成員基因結構的主要構成元件。

圖1 小黑麥TwSWEETs編碼的蛋白質二級結構Fig.1 Protein secondary structure encoded by TwSWEETs in triticale注:藍色為α-螺旋;綠色為β-折疊;黃色為無規則卷曲;紅色為延伸鏈Note:Blue,Alpha helix;green,Beta turn;yellow,Random coil;red,Extended strand

表3 小黑麥TwSWEETs編碼的蛋白質二級結構預測Table 3 Prediction of protein secondary structure encoded by triticale TwSWEETs

小黑麥TwSWEETs家族成員蛋白質三級結構預測結果表明(圖2),TwSWEET1a和TwSWEET1b的三級結構基本相同,TwSWEET2a和TwSWEET2b的三級結構預測結果也基本一致。TwSWEET6a/6b-1/6b-2的三級結構有較大的相似程度。TwSWEET12/13/14/15/16/17的蛋白質三級結構差異較大,但α-螺旋結構仍然是主要的組成結構原件。

2.3 小黑麥TwSWEETs蛋白序列的保守基序(motif)分析

基于保守基序分析發現,16小黑麥TwSWEETs蛋白由10個motif組成(圖3),16條TwSWEETs蛋白質均含有motif 1/2/3/4/5,保守性較強。motif 6在TwSWEET6a/6b-1/6b-2中出現。motif 7在TwSWEET1a/1b/4/6a/6b-1/6b-2/13/14中出現,但位置有所不同。motif 8在TwSWEET4/6a/6b-1/6b-2中出現。

如圖4所示,基序的長度范圍在8～50個氨基酸之間,其中motif1和motif2的基序長度均為50個氨基酸,且均含多個保守的氨基酸殘基位點,如motif 1中的甘氨酸(G,1,3)、脯氨酸(P,15,35)、亮氨酸(L,40)、酪氨酸(Y,36,49)和蘇氨酸(T,18);motif 2中的脯氨酸(P,3,20,49)、纈氨酸(V,10,16)、酪氨酸(Y,37)、谷氨酸(E,17)、亮氨酸(L,23)、絲氨酸(S,24)和天冬氨酸(D,43)。motif 3和motif 6的基序寬度均為29個氨基酸,motif4,5,7,8,9,10的基序寬度分別為22,21,8,10,18,17個氨基酸。其中,motif6,motif 8和motif 10含有較多個保守的氨基酸殘基位點,如motif6中的賴氨酸(K,1,2,3)、纈氨酸(V,6,23,25,27)和谷氨酸(E,7,28)等;motif8中的甲硫氨酸(M,1)、纈氨酸(V,2)和絲氨酸(S,3)等;motif10中的纈氨酸(V,6,8,9,12)、甘氨酸(G,11,15)和亮氨酸(L,13,14)等。此外,motif3/4/5/7/9所含有的保守氨基酸殘基位點較少,其中motif3/4/5各含有1個保守位點,motif7和motif9各含有2個保守位點。motif1/2/6/8/10包含多個保守的氨基酸殘基,由此推斷它們可能在小黑麥的糖轉運過程中發揮重要作用。

圖4 MEME程序鑒定的TwSWEETs蛋白的十個基序的序列標志Fig.4 Sequence markers of ten motifs of TwSWEETs protein identified by MEME program注:在每個位點都以字母(氨基酸殘基)排列的方式展示。氨基酸堆疊的總高度以位表示基序中該位點的“信息含量”。堆棧中每個字母的高度代表該位置的氨基酸的保守程度。X軸為基序的長度;Y軸為該位點的字母Note:Each site is displayed in alphabetical (amino acid residue) order. The total height of the amino acid stacking represents the information content of the site in the motif. The height of each letter in the stack represents the degree of conservation of the amino acids at that position. The X axis is the length of the motif, the Y axis is the site of the letter

2.4 小黑麥、小麥、擬南芥、水稻和玉米SWEETs糖轉運蛋白序列的系統發育分析

利用最大似然法(ML)構建了小黑麥、擬南芥、水稻和玉米的SWEETs蛋白的系統發育樹(圖5)。結果表明,小黑麥TwSWEETs蛋白家族成員可分為四個主要分支(Clade I,Clade II,Clade III和Clade IV),其中,Clade I包括TwSWEET1a/1b/2a/2b/3a,Clade II包括TwSWEET4/6a/6b-1/6b-2,Clade III包括Tw-SWEET11/12/13/14/15,Clade IV包括TwSWEE-T16和TwSWEET17。除此之外,16條小麥的TwSWEETs與單子葉植物小麥TaSWEETs、水稻OsSWEETs和玉米ZmSWEETs具有較高的同源性,而雙子葉植物擬南芥的多條AtSWEETs蛋白獨立于TwSWEETs,TaSWEETs,OsSWEETs和ZmSWEETs分布,與單子葉植物中的SWEETs基因家族的親緣關系相對較遠。

圖5 小黑麥、小麥、擬南芥、水稻和玉米SWEETs糖轉運蛋白序列的系統發育分析Fig.5 Phylogenetic analysis of SWEETs sugar transporter sequences in triticale,wheat,Arabidopsis,rice and maize注:SWEETs蛋白被分為四個分支Note:SWEETs protein is divided into four branches

2.5 TwSWEETs在非生物脅迫下的表達模式分析

2.5.1干旱脅迫 如圖6所示,干旱處理后,葉片和根中多個TwSWEETs基因的表達水平出現顯著上調或下調。脅迫6 h內,TwSWEET1b/2a/6a基因的表達在葉片中顯著上調,TwSWEET3a基因表達未受影響,隨著脅迫時間的延長,TwSWEET1b和SWEET6a基因的相對表達量仍顯著上調,于脅迫24 h時的相對表達量分別高達4.98和74.65,極顯著高于對照水平(P<0.01),TwSWEET12在葉片中的表達量呈現出干旱脅迫延長后表達量顯著增強的趨勢,而TwSWEET1a/2b/11/13/14/15/16/17在整個脅迫過程中葉片相對表達量均顯著低于對照(除TwSWEET14,12 h;TwSWEET17,3 h)(P<0.05)。根系表達模式分析發現,TwSWEET1b和TwSWEET12基因在脅迫12 h內顯著下調,24 h時極顯著上調表達(P<0.01)。根系TwSWEET4/11基因表達量在脅迫3 h時極顯著上調(P<0.01),后隨脅迫時間延長出現下調或恢復到與對照組相近水平。TwSWEET6a和TwSWEET16基因的相對表達量在脅迫24 h內均顯著高于對照水平,于脅迫 24 h時的相對表達量分別為70.44和3.50。

圖6 干旱脅迫下小黑麥TwSWEETs基因的表達模式分析Fig.6 Expression pattern analysis of TwSWEETs genes in triticale under drought stress注:*和**表示脅迫前后差異顯著(P<0.05,P<0.01)。下同Note:The difference before and after stress was significant (*P<0.05,**P<0.01). The same as below

2.5.2低溫脅迫 如圖7所示,低溫脅迫3 h時,葉片TwSWEET2a/2b/3a/14/16的相對表達量和對照CK相比顯著升高(P<0.05),其中TwSWEET2a/2b/14的相對表達量分別為4.84,10.32,2.32,極顯著高于對照水平(P<0.01),且TwSWEET16在低溫脅迫24 h內均表現為積極的上調響應,葉片TwSWEET2b/16在脅迫24 h后的相對表達量達到峰值,分別約為對照組CK的35倍和15倍(P<0.01);而葉片TwSWEET1a/11/13/17(除TwSWEET1a,3 h;TwSWEET13,3 h;TwSWEET1a,3 h;TwSWEET17,3/h、6 h;)的表達在低溫脅迫下受明顯抑制。低溫脅迫下各基因在根系中的表達分析發現,TwSWEET1a/2b/6a/11/14/16基因在脅迫期間的多個時間段內絕大多數都表現為上調表達,其中,根系TwSWEET6a和TwSWEET16在脅迫24 h時的相對表達量達到峰值,分別約為對照組CK的25倍和17倍,極顯著高于對照水平(P<0.01)。而根系TwSWEET2a/4/12/13/15基因的表達均受到低溫脅迫的顯著抑制(除TwSWEET15,24 h外)。

圖7 低溫脅迫下小黑麥TwSWEETs基因的表達模式分析Fig.7 Expression pattern analysis of TwSWEETs genes in triticale under low temperature stress

3 討論

糖不僅是植物生長和細胞代謝的能量和碳源,也是重要的信號分子[29-30]。SWEETs糖轉運蛋白家族不僅可以控制細胞內外及細胞器間糖的供應和分布,影響植物的“庫-源”運輸,而且在協調植物生長發育和逆境脅迫中也發揮重要作用[31-32]。本研究共鑒定到16條小黑麥TwSWEETs家族蛋白,較小麥、水稻、玉米等禾本科植物的成員數目少一些。一方面可能是因為小黑麥尚無全基因組,試驗所采用的數據庫僅為轉錄組數據庫,可能查找不全;另一方面可能是由于小黑麥自身在進化和親本雜交的過程中基因組發生了丟失現象[33]。對TMHs預測發現,11個TwSWEETs蛋白有7個TMHs結構,約占TwSWEETs基因總數的70%。5個TwSWEETs蛋白的TMHs少于7個,其中,TwSWEET1a/1b/6a/6b-2/17具有6個TMHs。研究表明,基因在不斷的進化過程中會伴隨外顯子化、丟失和插入等情況的發生,這可能是造成SWEET基因的TMHs少于7個的原因,Zhang等人[34]對草地早熟禾(Poapratensis)的SWEET家族分析發現,PpSWEETs家族成員存在4～7個TMHs[34-35],與本文的研究結果一致。

本研究基于SWEETs蛋白家族的系統發育分析,可將其根據蛋白質的同源性分為4個分支。前人研究結果表明,SWEETs家族不同分支運輸不同的糖。Clade I(SWEET1/2/3)和Clade II(SWEET 4/5/6/7/8)主要用于運輸葡萄糖和果糖,也可以運輸蔗糖;Clade III(SWEET9/10/11/12/13/14/15)主要用于蔗糖的高效轉運[36],最近的一項研究結果表明,如SWEET14和SWEET15能夠參與激素信號調控,可調節對赤霉素的反應[37];Clade IV(SWEET16/17)是一種液泡膜轉運體,研究發現其主要用于運輸葉片和根液泡膜中果糖[38]。本研究結果發現,小黑麥TwSWEETs家族也被分為4個分支,與前人研究結果一致,通常情況下,相似的結構意味著相似的功能,因此可推斷出小黑麥中不同蛋白成員轉運的糖類及功能不同,但每條蛋白具體的功能還需進一步深入研究。此外,研究發現,小黑麥TwSWEETs與禾本科其他作物同源性較高,且與雙子葉植物擬南芥相比,禾本科作物均沒有SWEET8和SWEET9,說明SWEET基因家族在單、雙子葉植物進化上存在差異。

糖是生物體維持正常生長發育、繁殖、代謝活動不可缺少的基礎組分。非生物脅迫會干擾植物體的代謝和光合活動,導致糖穩態被破壞。多項研究結果表明,植物在脅迫環境下,體內會積累大量可溶性糖,以此來降低體內滲透勢,維持正常生長,提高了植物抗逆能力[39]。SWEETs糖轉運蛋白能夠參與細胞器間、細胞間和組織器官間的糖轉運過程,同時積極參與植物對非生物脅迫的響應[40]。

姚利娜對茶樹中發現,AtSWEET4過表達植株中也積累了大量的葡萄糖和果糖,提高了植物的抗凍性[13]。在本研究中,干旱脅迫3 h時,葉片和根系TwSWEET4均出現顯著上調表達,低溫脅迫下,TwSWEET4表達均受到抑制,表明TwSWEET4參與了干旱和低溫脅迫的應答,但響應機制有所差異。Chen等對擬南芥進行滲透脅迫發現,其蔗糖轉運蛋白AtSWEET11/12會在干旱和ABA處理下被迅速磷酸化,該磷酸化會增強SWEETs蛋白的寡聚化和蔗糖轉運活性,從而使植物的根中含有更高的蔗糖含量,促進了根部生長,從而增強了植物在干旱脅下干旱抗性和增加了根冠比[31]。本研究結果發現,低溫脅迫下葉片TwSWEET11相對表達量相比對照顯著下調,而根系中的顯著上調,表明根系可能積累了大量蔗糖從而增加了小黑麥的抗寒性。

徐磊等人[40]研究結果表明,小麥TaSWEET6基因在低溫、干旱、鹽脅迫處理下均上調表達,本研究結果發現,干旱及低溫脅迫處理12 h后,TwSWEET6a在葉片及根部的表達水平顯著升高,極顯著高于對照水平,與小麥中的研究結果基本一致。此外,其他課題組有研究發現,干旱脅迫能夠上調擬南芥葉片和根中AtSWEET11/12的表達水平,促進了從葉到根的碳循環[38],本研究中小黑麥葉片TwSWEET12基因在干旱和低溫脅迫過程中表現為先降低再升高(低溫脅迫時在12 h后相對表達量會再次降低),但是TwSWEET11并沒有與TwSWEET12保持顯著的一致性,此外葉片中TwSWEET11在干旱和低溫脅長時間脅迫后期表達水平持續下調,這表明TwSWEET11/12基因在作物小黑麥中的調控模式有別于模式植物擬南芥。

本研究發現TwSWEET2a/6a/12/16在干旱和低溫處理后,葉片及根部表達量會發生急劇升高或下降,可將它們作為研究小黑麥對非生物脅迫脅迫響應的候選基因,如TwSWEET2a/6a/12/16可作為干旱脅迫候選基因,TwSWEET6a/16可作為寒冷脅迫候選基因。本研究為進一步揭示作物逆境生理代謝、信號轉導及優良品種選育提供了堅實的基礎。

4 結論

本研究從小黑麥轉錄組數據中共鑒定得到16條TwSWEETs糖轉運蛋白。所有TwSWEETs糖轉運蛋白均是疏水性蛋白,包含6～7個TMHs,亞細胞定位預測結果顯示多數TwSWEETs蛋白定位在質膜上,少數TwSWEETs存在多種定位結果。蛋白質二級和三級結構預測結果顯示α-螺旋結構是主要的組成原件。系統發育樹分析將其分為四個亞族,家族成員與禾本科多個作物的SWEET蛋白家族同源性高,保守性較強。TwSWEETs基因積極參與了小黑麥對非生物脅迫的響應過程,其中,TwSWEET2a/6a/12/16在干旱或低溫脅迫下呈現出顯著差異表達,可將它們作為研究小黑麥對環境脅迫反應的重要候選基因。