賈玉庫 高宏歡 馮健超 郝紫瑞 王晨陽 謝迎新 郭天財 馬冬云

研究簡報

小麥G2-like轉(zhuǎn)錄因子家族基因鑒定與表達模式分析

賈玉庫**高宏歡**馮健超 郝紫瑞 王晨陽 謝迎新 郭天財 馬冬云*

河南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院 / 國家小麥工程技術(shù)研究中心/ 河南省小麥技術(shù)創(chuàng)新中心, 河南鄭州 450046

()轉(zhuǎn)錄因子, 是屬于MYB類轉(zhuǎn)錄因子中GARP超家族的成員, 在調(diào)節(jié)葉綠體發(fā)育中起重要作用。本研究利用生物信息學(xué)方法對小麥基因進行了全基因組鑒定, 并對其理化性質(zhì)、亞細胞定位、啟動子順式作用元件及對非生物脅迫和激素的響應(yīng)模式進行了分析。從小麥中共鑒定出87個基因, 不均勻的分布在小麥21條染色體上, 系統(tǒng)發(fā)育分析將這些基因分為14個亞族。蛋白質(zhì)二級結(jié)構(gòu)預(yù)測表明, 小麥基因的氨基酸序列均以α螺旋和隨機卷曲為主要結(jié)構(gòu)。啟動子順式作用元件分析表明其上游2 kb區(qū)域含有7種(P-box、SpI、LTR、ABRE、MBS、TGA-element和AE-box)與逆境脅迫誘導(dǎo)相關(guān)順式作用元件。其中含有順式調(diào)控元件結(jié)合位點最多, 一共有18個結(jié)合位點。qRT-PCR驗證發(fā)現(xiàn),、、和在PEG和鹽脅迫下以及GA、IAA和ABA激素誘導(dǎo)下表達量顯著上調(diào), 這些基因可能介導(dǎo)了小麥對多種非生物逆境的響應(yīng)。

小麥; G2-like轉(zhuǎn)錄因子; 生物信息學(xué); 逆境脅迫; 表達分析

小麥是我國最主要的糧食作物之一, 其產(chǎn)量高低對國家糧食安全至關(guān)重要。近年來, 隨著氣候變化和農(nóng)業(yè)生態(tài)環(huán)境污染, 各種環(huán)境脅迫包括低溫[1-2]、高溫[3]、干旱[4-5]和鹽脅迫[6-7]等嚴重影響小麥生長。因此, 提升小麥在逆境條件下的防御水平, 對小麥增產(chǎn)增收意義重大。

植物中轉(zhuǎn)錄因子在逆境脅迫下可以激活或者抑制某些基因的轉(zhuǎn)錄, 影響其蛋白在植物體內(nèi)的表達量及功能, 從而在植物逆境應(yīng)答和生長發(fā)育中起重要作用[8-9]。(或)轉(zhuǎn)錄因子是屬于MYB類轉(zhuǎn)錄因子中GARP超家族的成員, 廣泛存在于陸地植物中。首先在玉米中被發(fā)現(xiàn)[10], 隨后在擬南芥、水稻、辣椒和番茄中均有存在, 在葉綠體發(fā)育、果實品質(zhì)、植物衰老及響應(yīng)逆境脅迫等方面均起作用[11-13]。e家族成員的特征是在過渡和早期成熟階段調(diào)節(jié)葉綠體的形成[14-15], 對葉綠體的發(fā)育是不可或缺的基因。玉米和這一對同源基因功能基本相同, 分別在葉肉細胞和維管束鞘細胞中表達[16]。過表達可以增強調(diào)控葉綠體的發(fā)育相關(guān)基因的表達, 進而影響成熟果實中的糖類和類胡蘿卜素的含量[13]。

基因沉默后的番茄植株表現(xiàn)出對干旱和低溫的抗性下降[17]。擬南芥在低溫下基因通過抑制花青素合成基因和抑制花青素的積累并增強黃酮和香豆素苷的積累, 從而對低溫產(chǎn)生抗性[18]。對玉米植株在低溫(4℃)和干旱脅迫處理下的轉(zhuǎn)錄因子家族的20個基因表達量分析, 發(fā)現(xiàn)部分基因的表達與植株抗逆性有關(guān)[19]。同時轉(zhuǎn)錄因子影響控制氣孔開閉和鉀離子運動相關(guān)基因的表達, 從而對氣孔起到調(diào)節(jié)作用[20]。植物激素如生長素和油菜素類固醇與果實葉綠體的發(fā)育有關(guān), 且與的表達上調(diào)有關(guān)[21]。目前關(guān)于轉(zhuǎn)錄因子在小麥中的研究尚未有詳細報道。因此, 通過對小麥基因組中轉(zhuǎn)錄因子進行了系統(tǒng)挖掘, 結(jié)合生物信息學(xué)分析和對部分小麥轉(zhuǎn)錄因子在干旱和鹽脅迫下的表達模式分析, 篩選抗逆相關(guān)的家族成員, 為小麥抗逆基因的挖掘和在育種中的應(yīng)用奠定基礎(chǔ)。

1 材料與方法

1.1 小麥G2-like基因的全基因組鑒定、編碼蛋白的理化性質(zhì)、二級結(jié)構(gòu)分析和亞細胞定位分析

從Ensembl Plants數(shù)據(jù)庫(http://plants.ensembl.org/ index.html)獲取小麥蛋白序列, 以玉米G2-like蛋白序列作為query序列進行BLASTP比對(-value<1E–5)獲得同源序列。在Pfam數(shù)據(jù)庫(http://pfam.xfam.org/)下載G2-like保守結(jié)構(gòu)域PF00249.29, 利用HMM3.0進行保守結(jié)構(gòu)域比對; 同時利用Pfam搜索(http://pfam.xfam.org/)、SMART (http://samrt.embl-heidelberg.de/)和HMMER (https://www. ebi.ac.uk/Tools/hmmer/search/hmscan)驗證候選蛋白保守結(jié)構(gòu)域。利用ExPASy-ProtParam (https://web.expasy.org/ protparam/)計算G2-like蛋白的氨基酸長度和分子量等基本理化性質(zhì)。分別利用軟件NPS-SOPMA (https://npsa- prabi.ibcp.fr/cgi-bin/secpred_sopma.pl)和ProtComp 9.0 (https://linux1.Softbe-rry.com/berry.phtml)預(yù)測G2-like的二級結(jié)構(gòu)和亞細胞定位。

1.2 小麥G2-like基因系統(tǒng)發(fā)育進化、基因結(jié)構(gòu)和基序分析

利用Clustal_W工具對G2-like蛋白序列進行多重比對, 結(jié)果放入MEGA 7.0[22]軟件, 采用鄰接法(Neighbor- Jointing, NJ)構(gòu)建系統(tǒng)發(fā)育進化樹(bootstrap, 1000; 基序, 10)。利用TBtools軟件[23]作圖。

1.3 小麥G2-like基因系統(tǒng)發(fā)育進化、染色體定位和共線性分析

根據(jù)IWGSC (http://wheat-urgi.versailles.inra.fr/)小麥基因組數(shù)據(jù)BLASTN構(gòu)建本地數(shù)據(jù)庫, 通過比對獲得基因在染色體上的物理位置。利用MG2C (http://mg2c.iask.in/mg2c_v2.0/)構(gòu)建基因在染色體上的物理圖譜。利用MCScanX[24]進行共線性分析(-value<1E–5), Circos軟件可視化分析[25]。利用KaKs_Calculator軟件計算非同義替代率(Ka)和同義替代率(Ks)[26], 預(yù)測基因CDS區(qū)的適應(yīng)性進化。

1.4 小麥G2-like基因的順式作用元件分析

利用TBtools軟件截取基因CDS上游2.0 kb的DNA序列, 提交至PlantCARE數(shù)據(jù)庫(http://bioinformatics. psb.ugent.be/webtools/plantcare/html)進行啟動子順式作用元件的預(yù)測。

1.5 小麥G2-like基因在不同組織部位、逆境脅迫下的表達分析

利用expVIP數(shù)據(jù)庫(http://www.wheat-expression.com/)和Wheat Exp數(shù)據(jù)庫(https://wheat.pw.usda.gov/WheatExp/)中小麥基因組和在干旱、熱脅迫下的轉(zhuǎn)錄組數(shù)據(jù), 分析基因在不同組織及逆境脅迫下表達模式。

1.6 小麥G2-like轉(zhuǎn)錄因子家族部分基因qRT-PCR驗證

選用河南省主推小麥品種百農(nóng)207為試驗材料, 挑選籽粒飽滿大小一致的種子經(jīng)75%乙醇消毒處理后, 用蒸餾水沖洗5次, 然后置于光照培養(yǎng)箱中。在23℃/18℃ (白天/黑夜)、光/暗周期為16 h/8 h培養(yǎng), 然后挑選發(fā)芽一致的種子采用Hoagland營養(yǎng)液培養(yǎng), 每2 d更換一次培養(yǎng)液。待其長到二葉一心進行干旱(20% PEG-6000)、鹽脅迫(200 mmol L–1NaCl)、以及生長素(100 μmol L–1IAA)、赤霉素(100 μmol L–1GA)和脫落酸(100 μmol L–1ABA)處理, 正常營養(yǎng)液培養(yǎng)作為對照, 所有處理重復(fù)3次。分別在處理0、1、6、12和24 h時取不同處理幼葉, 立即置于液氮中保存。使用TRIZOL (北京全式金生物技術(shù)股份有限公司)提取RNA, 進行反轉(zhuǎn)錄cDNA第一鏈合成, 用于Real-time PCR分析(表1)。在QuantStudio 5 Real-time PCR System (Themo Fisher, BIO INC, CA, 美國)上進行qRT-PCR分析, 反應(yīng)體系與程序參考TB Green Premix ExII (TaKaRa, 大連)試劑盒說明書, 采用2–ΔΔCt法計算相對表達量[27]。

表1 部分G2-like基因熒光定量所用引物

2 結(jié)果與分析

2.1 小麥G2-like基因全基因組鑒定、編碼蛋白的基本理化性質(zhì)、二級結(jié)構(gòu)和亞細胞定位分析

根據(jù)59個玉米基因的蛋白序列, 利用BLASTP和HMM比對搜索, 在小麥基因組數(shù)據(jù)庫中得到110個候選基因。通過Pfam、SMART、HMMER驗證是否有G2-like MYB DNA保守結(jié)構(gòu)域, 共鑒定出87個基因, 據(jù)其在染色體上的位置進行命名(表2)。

蛋白序列一級結(jié)構(gòu)理化性質(zhì)分析顯示, 87個G2-like蛋白序列長度為224～534氨基酸, 分子量為24.03～58.47 kD, 等電點范圍為4.84～9.65, 其中43個蛋白理論等電點在酸性范圍內(nèi), 其余均在堿性范圍內(nèi), 表明G2-like蛋白質(zhì)酸堿分布均勻; 疏水性指數(shù)范圍為–0.945～ –0.319; 不穩(wěn)定系數(shù)范圍為28.00～73.87, 其中有8個G2-like蛋白不穩(wěn)定系數(shù)小于40, 相對來說為穩(wěn)定蛋白質(zhì)。蛋白二級結(jié)構(gòu)均含有α螺旋、延伸鏈、β轉(zhuǎn)角和無規(guī)則卷曲, 但以α螺旋和無規(guī)則卷曲為主要結(jié)構(gòu)。亞細胞定位預(yù)測表明, 大部分位于核膜, 3個位于線粒體膜, 9個位于質(zhì)膜, 1個位于細胞質(zhì)膜。

2.2 G2-like基因家族系統(tǒng)發(fā)育進化、基因結(jié)構(gòu)和基序分析

基于進化關(guān)系、基序和結(jié)構(gòu)分析,基因家族可分為14 (I～XIV)組(圖1-A)。每組基因分布不均勻, 第9組數(shù)量最多(13個), 第3組最少只有2個。一般來說, 聚集在同一分組的基因表現(xiàn)出相似的基序結(jié)構(gòu)和基因結(jié)構(gòu)。利用MEME motif搜索工具, 鑒定出G2-like蛋白10個保守基序(圖1-B)。幾乎每個基因都有基序1和基序2結(jié)構(gòu)域(沒有)。此外, 每個分組都有獨特的基序, 說明轉(zhuǎn)錄因子在進化過程中存在內(nèi)部分化, 這可能是基因功能多樣性的來源。依據(jù)基因組DNA和CDS序列, 獲得了基因的外顯子和內(nèi)含子結(jié)構(gòu)(圖1-C)。所有的e基因都具有完整的基因結(jié)構(gòu), 但是有14個基因沒有UTR結(jié)構(gòu)。當缺失UTR時, 意味著基因的調(diào)控功能(mRNA的穩(wěn)定性、折疊、核內(nèi)運輸相互作用, RNA加工、剪切和翻譯等)會受到影響[28-29]。

2.3 小麥G2-like基因的染色體分布及共線性分析

基因在小麥A、B、D三種染色體上是不均勻分布, 且不同染色體上基因位置和數(shù)目明顯不同(圖2)。其中2A、2B和2D染色體上分布的最少(7個), 在7A、7B和7D染色體上分布的最多(22個)。有6個基因形成3對串聯(lián)重復(fù)基因?qū)?和、和和

通過BLASTP和MCscanX對進行重復(fù)模式分析, 表明小麥基因組內(nèi)存在79,067個共線性區(qū)域(包括Un染色體), 共存在107,545個基因(圖3)。有74個小麥基因定位在共線性區(qū)域內(nèi), 說明這些基因由片段復(fù)制模式產(chǎn)生。同時, 3個基因同時涉及到片段復(fù)制和串聯(lián)重復(fù)事件, 表明這些基因由兩種重復(fù)模式互作擴增而來。家族的同源基因?qū)Φ腒a/Ks值均小于1, 表明這些同源基因?qū)艿截撨x擇作用, 同時其進化過程中受到功能限制作用。

2.4 小麥G2-like基因啟動子的順式作用元件分析

對基因啟動子序列分析表明, 其2 kb上游區(qū)域含有響應(yīng)抗逆脅迫的7種潛在順式作用調(diào)控元件(圖4), 其中與激素、干旱、低溫、鹽脅迫誘導(dǎo)有關(guān)的MBS、P-box、Sp1、TGA-element、LTR、ABRE和AE-box順式作用元件。本研究中62個基因至少含有5個抗逆脅迫順式作用元件結(jié)合位點,含有順式調(diào)控元件結(jié)合位點最多, 一共有18個結(jié)合位點。表明小麥基因家族可能參與多種非生物逆境的響應(yīng)。

圖1 小麥G2-like基因的系統(tǒng)發(fā)育關(guān)系、保守蛋白基序結(jié)構(gòu)和基因結(jié)構(gòu)

A: 利用MEGA7基于小麥G2-like蛋白全長序列構(gòu)建系統(tǒng)發(fā)育樹。B: 小麥G2-like蛋白的基序組成。不同顏色表示motif 1～motif 10; 蛋白質(zhì)的長度可以用底部的刻度來估計。C: 小麥基因的外顯子-內(nèi)含子結(jié)構(gòu)。綠色方框表示未翻譯的5-和3-區(qū)域;黃色框表示外顯子; 黑線表示內(nèi)含子。

A: phylogenetic tree was constructed based on the full-length sequences of wheat G2-like proteins using MEGA 7. B: the motif composition of the wheat G2-like proteins. Different colors indicate motif 1–motif 10. The length of the protein can be estimated using the scale at the bottom. C: the exon-intron structures of wheatgenes. Green boxes indicate untranslated 5-regions and 3-regions; yellow boxes indicate exons; and black lines indicate introns.

圖2 小麥G2-like基因染色體定位

標度代表兆堿基(Mb)。The scale represents Mb.

圖3 小麥G2-like基因的染色體位置及共線性關(guān)系

所有的共線區(qū)域和基因由灰線連接,基因的片段重復(fù)基因?qū)τ眉t線標出。

All syntenic blocks and genes are linked by the grey lines, and segmental duplication pairs ofgenes are highlighted by red lines.

圖4 小麥G2-like基因啟動子區(qū)域與抗逆相關(guān)的7種順式作用元件

順式作用元件用不同的彩色方框表示。

Different colored boxes indicate different cis-acting elements in the scale at the bottom.

2.5 小麥G2-like基因組織表達及逆境脅迫下表達分析

基因的表達模式分析表明(圖5),基因家族可分為5類(A～E), 其中大部分基因在葉中均有表達; 模式A和E中的基因幾乎在根、葉、穗和籽粒的所有時期都有較高水平表達; 模式B中的基因主要在根中表達; 模式D主要是在葉和穗中表達; 在模式C中有12個的基因在不同組織和不同時期的表達量都很低或沒有表達。

對候選小麥基因在干旱、高溫及干旱和高溫雙重脅迫下的表達情況分析(圖6), 發(fā)現(xiàn)大部分基因在脅迫下均有較高的表達量, 表達模式基本相同。、、、、、和在3類脅迫下表達量均較高; 特別是, 其表達量在3種脅迫下最高, 表明這個基因在小麥遭受熱、干旱脅迫時可能發(fā)揮重要作用。、、和在干旱脅迫下的表達量顯著高于熱脅迫, 表明這4個基因?qū)Ω珊得{迫可能更加敏感。但是、、、和在不同的脅迫下的表達量均較低或者沒有, 可能是這5個基因?qū)δ婢趁{迫不敏感。

圖5 不同時期籽粒、葉片、根和穗中87個G2-like基因表達譜熱圖

S: 幼苗期; V: 期; R: 生殖期。紅色或綠色代表每個樣品中每個轉(zhuǎn)錄本的較高或較低的相對豐度。

S: at seedling stage; V: at vegetative stage; R: at reproductive stage. The red or green colors stand for the higher or lower relative abundance of each transcript in each sample.

圖6 小麥G2-like基因在熱、干旱以及熱和干旱共同脅迫下的表達模式

D1、D6:干旱處理1、6 h; H1、H6: 高溫處理1、6 h; D+H1、D+H6: 在高溫和干旱脅迫下處理時間為1、6 h。紅色和綠色分別代表每個樣品中每個轉(zhuǎn)錄本的較高和較低的相對豐度。

D1, D6: drought treatment of 1 h and 6 h; H1, H6: heat treatment of 1 h and 6 h; D+H1, D+H6: heat and drought treatments of 1h and 6 h. The red and green colors represent the higher or lower relative abundance, respectively.

2.6 小麥G2-like轉(zhuǎn)錄因子家族部分基因qRT-PCR驗證分析

為探究基因?qū)Ψ巧锩{迫及激素的響應(yīng), 根據(jù)基因在不同組織部位表達量, 以及啟動子區(qū)域順式調(diào)控元件結(jié)合位點數(shù)量等綜合考慮, 挑選了14個小麥基因進行qRT-PCR表達分析(圖7)。結(jié)果表明, 在鹽脅迫下只有表達量顯著降低, 在脅迫24 h后較CK下降了3倍; 而、、、、、、、和表達量顯著升高, 其中、升高最為明顯, 在脅迫24 h后其表達量相對于CK分別提高了8倍、15倍和11倍。

在干旱脅迫下,、、、和、、和表達量顯著升高, 其中和上升最為顯著, 在脅迫24 h后其表達量相對于CK分別升高了10倍、13倍和12倍; 而在脅迫下上調(diào)最不明顯, 在脅迫24 h后相對于CK上升大致為1.3倍。同時發(fā)現(xiàn)在鹽脅迫下顯著上調(diào), 但在干旱脅迫下其表達量并沒有顯著變化(脅迫12 h除外);在鹽脅迫下顯著下調(diào), 而在干旱脅迫6 h之后其表達量較對照顯著上調(diào), 表明這2個基因在響應(yīng)干旱和鹽脅迫中可能存在差異。

在3種激素誘導(dǎo)下,、、、、、和表達量均顯著升高(圖8), 其中和升高最為明顯, 其表達量相對CK上調(diào)10倍以上。、和表達量隨GA和IAA誘導(dǎo)時間的延長呈現(xiàn)出先升高后降低的趨勢, 且均在誘導(dǎo)1 h之后表達量顯著上調(diào)。但其對ABA誘導(dǎo)的響應(yīng)相對較弱,和分別在ABA誘導(dǎo)6 h和24 h、6 h時表達量較對照顯著上調(diào); 而在ABA處理12 h之后表達量顯著上調(diào)。在IAA處理下表達量顯著升高, 在誘導(dǎo)24 h后較CK上升3倍左右; 而在3種激素處理下表達量均表現(xiàn)出下降的趨勢, 其中在IAA處理下降最明顯, 在誘導(dǎo)24 h后較CK下調(diào)了0.3倍。

3 討論

小麥在生產(chǎn)過程經(jīng)常遭受各種逆境脅迫的影響, 挖掘小麥抗逆基因資源對于小麥抗逆育種有重要作用。植物在面臨逆境脅迫時, 其體內(nèi)轉(zhuǎn)錄因子迅速啟動相應(yīng)基因的表達和逆境信號的傳遞, 啟動抗逆應(yīng)答反應(yīng), 在植物逆境應(yīng)答和生長發(fā)育過程中發(fā)揮重要作用[9]。轉(zhuǎn)錄因子, 是屬于MYB類轉(zhuǎn)錄因子中GARP超家族的成員, 廣泛存在于陸地植物中。本研究在小麥中共鑒定出87個基因, 不均勻的分布在21條染色體上; 通過共線性分析發(fā)現(xiàn)小麥基因家族主要由片段重復(fù)產(chǎn)生, 同時還有3個串聯(lián)重復(fù)基因?qū)? 植物基因家族擴張的主要途徑是串聯(lián)重復(fù)和片段重復(fù)[30], 片段復(fù)制在進化較慢的基因家族中是有規(guī)律出現(xiàn)的。計算Ka/Ks值后發(fā)現(xiàn)2種重復(fù)基因?qū)Ρ戎稻∮?, 表明基因在進化過程中可能是相對保守的, 進化較慢。

圖7 定量實時PCR分析14個G2-like基因?qū)aCl、PEG處理的響應(yīng)

*表示處理間存在顯著差異(< 0.05)。

* indicates significant differences between the stress conditions and the control condition (< 0.05).

(圖8)

圖8 定量實時PCR分析14個G2-like基因?qū)BA、GA和IAA處理的響應(yīng)

*表示處理間存在顯著差異(< 0.05)。

* indicates significant differences between the stress conditions and the control condition at< 0.05.

位于啟動子區(qū)域的順式作用元件通過對不同環(huán)境信號的響應(yīng)來調(diào)控基因轉(zhuǎn)錄, 從而影響植物的生長發(fā)育[31]。在玉米中發(fā)現(xiàn)與激素、干旱和低溫誘導(dǎo)響應(yīng)有關(guān)的MBS、P-box、Sp1、TGA-element和AE-box順式作用元件[32]; 在水稻中發(fā)現(xiàn)的ABRE順式調(diào)控元件, 參與干旱誘導(dǎo)和脫落酸誘導(dǎo)[33-34]; LTR順式調(diào)控元件, 在擬南芥鹽脅迫中發(fā)揮重要作用[35]。本研究發(fā)現(xiàn), 小麥基因啟動子區(qū)域含有多種順式作用元件, 包括MBS、P-box、ABRE等響應(yīng)逆境脅迫元件, 表明小麥基因在參與逆境脅迫中可能發(fā)揮重要作用。已有研究結(jié)果表明基因在植物葉綠體的發(fā)育具有重要作用[36]。對不同組織部位基因表達模式分析表明, 不同基因在小麥葉片中均有表達, 且在干旱、高溫以及干旱和高溫共同脅迫下,基因都有較高的表達量, 表明基因可能在葉綠體發(fā)育和逆境脅迫響應(yīng)中均發(fā)揮著重要作用, 但其具體的調(diào)控機制仍需進一步研究。

玉米中有10個基因在干旱和低溫脅迫下表達量均上調(diào), 其中基因在低溫和干旱脅迫下表達量顯著上調(diào)[19]。劉俊芳等[17]發(fā)現(xiàn)鹽處理后有6個基因表達上調(diào), 干旱處理后有3個基因表達上調(diào), 其中基因在干旱和低溫脅迫中應(yīng)答明顯, 沉默后番茄植株對干旱和低溫的抗性降低。本研究中通過實時熒光定量分析發(fā)現(xiàn), 14個基因中, 幾乎所有的基因在鹽脅迫和干旱脅迫下顯著上調(diào), 表明多數(shù)基因在鹽和干旱脅迫下表現(xiàn)為正向調(diào)控; 同時發(fā)現(xiàn)基因在干旱脅迫下顯著下調(diào); 表明盡管轉(zhuǎn)錄因子的蛋白氨基酸序列相似, 但轉(zhuǎn)錄因子對干旱脅迫的調(diào)控機制不盡相同。其中、和在這兩種脅迫下表達量上調(diào)幅度最高, 可能在小麥抗鹽和抗旱中發(fā)揮重要作用, 有關(guān)其功能驗證及作用機理有待進一步研究。

ABA信號轉(zhuǎn)導(dǎo)途徑是植物響應(yīng)鹽和干旱脅迫的關(guān)鍵途徑之一;雙突變體在擬南芥種子萌發(fā)和幼苗發(fā)育期間表現(xiàn)出對ABA高度敏感, 在幼苗發(fā)育期間也有抗?jié)B透脅迫的作用;基因調(diào)控包括WRKY40在內(nèi)的幾種ABA響應(yīng)基因, 在ABA反應(yīng)中起負調(diào)控作用[37]。本研究中多數(shù)基因在ABA誘導(dǎo)下顯著上調(diào), 與干旱脅迫下表達模式相似, 表明這些轉(zhuǎn)錄因子參與干旱信號轉(zhuǎn)導(dǎo)途徑可能需借助ABA信號途徑; 而表達量對ABA響應(yīng)程度較小, 可能其不依賴ABA信號途徑對抗旱能力進行調(diào)控。如在水稻的研究中表明參與介導(dǎo)植物對ABA和干旱的反應(yīng)[38],而過表達植株對外源ABA沒有響應(yīng)[39]。在沉默了抑制生長素響應(yīng)的基因株系中,的表達升高[40]。本研究14個基因中, 幾乎所有的基因均在IAA和GA誘導(dǎo)下顯著上調(diào), 只有和分別在IAA和GA誘導(dǎo)下顯著下調(diào)表達, 表明不同轉(zhuǎn)錄因子對激素的響應(yīng)模式存在一定差異, 可能與其參與植株生長發(fā)育和逆境脅迫響應(yīng)的調(diào)控模式有關(guān)。多數(shù)基因?qū)γ{迫和激素誘導(dǎo)均明顯響應(yīng), 其中、、和在不同脅迫和激素處理下應(yīng)答最為明顯, 表明這4個基因可能在小麥抗逆應(yīng)答過程中發(fā)揮著重要作用。同時必須指出的是,本研究在逆境脅迫和激素誘導(dǎo)基因表達分析中, 僅選用了1個品種, 且脅迫類型和激素種類較少; 盡管百農(nóng)207為適應(yīng)性好、綜合抗逆能力強的品種, 但選擇抗逆能力不同小麥品種、同時增加逆境脅迫和激素種類可能更有利于抗逆相關(guān)基因鑒定, 關(guān)于基因?qū)Σ煌婢趁{迫和激素誘導(dǎo)的響應(yīng)需進一步驗證。根據(jù)基因表達模式并結(jié)合生信分析,可能在小麥抗逆中發(fā)揮重要作用。通過載體構(gòu)建和遺傳轉(zhuǎn)化, 已經(jīng)成功獲得了10株過表達的轉(zhuǎn)基因株系, 有關(guān)過表達株系在逆境下的表型及相應(yīng)生理機制有待進一步研究。

[1] Janmohannadi M, Zolla L, Rinalducci S M. Low temperature tolerance in plants: changes at the protein level., 2015, 117: 76–89.

[2] Rinalducci S, Egidi M G, Karimzadeh G, Jazii F R, Zolla L. Proteomic analysis of a spring wheat cultivar in response to prolonged cold stress., 2011, 32: 1807–1818.

[3] Gibson L R, Paulsen G. Yield components of wheat grown under high temperature stress during reproductive growth., 1999, 39: 1841–1846.

[4] Farooq M, Hussain M, Siddique K H M. Drought stress in wheat during flowering and grain-filling periods., 2014, 33: 331–349.

[5] Mohammadi R. Efficiency of yield-based drought tolerance indices to identify tolerant genotypes in durum wheat., 2016, 211: 71–89.

[6] Ding H, Ma D, Huang X, Hou J, Wang C, Xie Y, Wang Y, Qin H, Guo T. Exogenous hydrogen sulfide alleviates salt stress by improving antioxidant defenses and the salt overly sensitive pathway in wheat seedlings., 2019, 41: 123.

[7] Saqib M, Akhtar J, Qureshi R H. Pot study on wheat growth in saline and waterlogged compacted soil: I. Grain yield and yield components., 2004, 77: 169–177.

[8] Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AP2/ERF family transcription factors in plant abiotic stress responses., 2012, 1819: 86–96.

[9] Chrispeels H E, Oettinger H, Janvier N, Tague B W., encodingC2H2zinc-finger protein 1, is expressed downstream of photomorphogenic activation., 2000, 42: 279–290.

[10] Jenkins M T. A second gene producing golden plant color in maize., 1926, 60: 484–488.

[11] Brand A, Borovsky Y, Hill T, Rahman K A A, Bellalou A, Van Deynze A, Paran I.regulates natural variation of chlorophyll content and fruit color in pepper fruit., 2014, 127: 2139–2148.

[12] Fitter D W, Martin D J, Copley M J, Scotland R W, Langdale J A.gene pairs regulate chloroplast development in diverse plant species., 2002, 31: 713–727.

[13] Powell A L T, Nguyen C V, Hill T, Cheng K L, Figueroa-Balderas R, Aktas H, Ashrafi H, Pons C, Fernández-Mu?oz R, Vicente A, Lopez-Baltazar J, Barry C S, Liu Y S, Chetelat R, Granell A, Van Deynze A, Giovannoni J J, Bennett A B.encodes atranscription factor regulating tomato fruit chloroplast development., 2012, 336: 1711–1715.

[14] Jarvis P, López-Juez E. Biogenesis and homeostasis of chloroplasts and other plastids., 2013, 14: 787–802.

[15] Li P, Ponnala L, Gandotra N, Wang L, Si Y, Tausta S L, Kebrom T H, Provart N, Patel R, Myers C R, Reidel E J, Turgeon R, Liu P, Sun Q, Nelson T, Brutnell T P. The developmental dynamics of the maize leaf transcriptome., 2010, 42: 1060–1067.

[16] Chang Y M, Liu W Y, Shih A C-C, Shen M N, Lu C H, Lu M-Y J, Yang H W, Wang T Y, Chen S C-C, Chen S M, Li W H, Ku M S B. Characterizing regulatory and functional differentiation between maize mesophyll and bundle sheath cells by transcriptomic analysis., 2012, 160: 165–177.

[17] 劉俊芳. 番茄G2-like轉(zhuǎn)錄因子家族生物信息學(xué)分析及抗逆相關(guān)基因鑒定. 東北農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 黑龍江哈爾濱, 2018.

Liu J F. Bioinformatics Analysis of Tomato G2-like Transcription Factor Family and Identification of Resistance-related Genes. MS Thesis of Northeast Agricultural University, Harbin, Heilongjiang, China, 2018 (in Chinese with English abstract).

[18] Petridis A, D?ll S, Nichelmann L, Bilger W, Mock H P.G2-LIKE FLAVONOID REGULATOR and BRASSINOSTEROID ENHANCED EXPRESSION1 are low-temperature regulators of flavonoid accumulation., 2016, 211: 912–925.

[19] Liu F, Xu Y, Han G, Zhou L, Ali A, Zhu S, Li X. Molecular evolution and genetic variation oftranscription factor genes in maize., 2016, 11: e0161763.

[20] Nagatoshi Y, Mitsuda N, Hayashi M, Inoue S, Okuma E, Kubo A, Murata Y, Seo M, Saji H, Kinoshita T, Oheme-Takagi M. GOLDEN 2-LIKE transcription factors for chloroplast development affect ozone tolerance through the regulation of stomatal movement., 2016, 113: 4218–4223.

[21] Chen M, Ji M, Wen B, Liu L, Li S, Chen X, Gao D, Li L. GOLDEN 2-LIKE transcription factors of plants., 2016, 7: 1509.

[22] Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets., 2016, 33: 1870–1874.

[23] Chen C, Xia R, Chen H, He Y. TBtools, a toolkit for biologists integrating various HTS-data handling tools with a user-friendly interface., 2018, 13: 1194–1202.

[24] Wang Y, Tang H, DeBarry J D, Tan X, Li J, Wang X, Lee T, Jin H, Marler B, Guo H, Kissinger J C, Paterson1 A H. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity., 2012, 40: e49.

[25] Krzywinski M, Schein J, Birol ?, Connors J, Gascoyne R, Horsman D, Jones S J, Marra M A. Circos: an information aesthetic for comparative genomics., 2009, 19: 1639–1645.

[26] Zhang Z, Li J, Zhao X Q, Wang J, Wong G K S, Yu J. KaKs_Calculator: calculating Ka and Ks through model selection and model averaging., 2006, 4: 259–263.

[27] Kang C H, Jung W Y, Kang Y H, Kim J Y, Kim D G, Jeong J C, Baek D W, Jin J B, Lee J Y, Kim M O, Chung W S, Mengiste T, Koiwa H, Kwak S S, Bahk J D, Lee S Y, Nam J S, Yun D J, Cho M J. AtBAG6, a novel calmodulin-binding protein, induces programmed cell death in yeast and plants., 2006, 13: 84–95.

[28] Koundrioukoff S, Polo S, Almouzni G. Interplay between chromatin and cell cycle checkpoints in the context of ATR/ATM- dependent checkpoints., 2004, 3: 969–978.

[29] Meshorer E. Chromatin in embryonic stem cell neuronal differentiation., 2007, 22: 311–319.

[30] Cannon S B, Mitra A, Baumgarten A, Young N D, May G. The roles of segmental and tandem gene duplication in the evolution of large gene families in., 2004, 4: 10.

[31] Kong F, Wang J, Cheng L, Liu S, Wu J, Peng Z, Lu G. Genome-wide analysis of the mitogen-activated protein kinase gene family in., 2012, 499: 108–120.

[32] Tao Y, Wang F, Jia D, Li J, Zhang Y, Jia C, Wang D, Pan H. Cloning and functional analysis of the promoter of a stress- inducible gene () in maize., 2014, 33: 200–208.

[33] Lee S C, Kim S H, Kim S R. Drought induciblepromoter is activated by OsDREB1A and OsDREB1D., 2013, 56: 115–121.

[34] Nakashima K, Jan A, Todaka D, Maruyama K, Goto S, Shinozaki K, Yamaguchi?Shinozaki K. Comparative functional analysis of six drought-responsive promoters in transgenic rice., 2014, 239: 47–60.

[35] 張新宇, 趙蘭杰, 李艷軍. 鹽脅迫對擬南芥基因的誘導(dǎo)表達及其啟動子分析. 西北植物學(xué)報, 2014, 34(1): 54–59.

Zhang X Y, Zhao L J, Li Y J. Salt Stress induced expression and promoter analysis ofgene in, 2014, 34(1): 54–59 (in Chinese with English abstract).

[36] Yasumura Y, Moylan E C, Langdale J A. A conserved transcription factor mediates nuclear control of organelle biogenesis in anciently diverged land plants., 2005, 17: 1894–1907.

[37] Ahmad R, Liu Y, Wang T J, Meng Q, Yin H, Wang X, Wu Y, Nan N, Liu B, Zheng Y X. GOLDEN2-LIKE transcription factors regulateexpression in response to abscisic acid., 2019, 179: 1844–1860.

[38] Yin X, Cui Y, Wang M, Xia X. Overexpression of a novel MYB-related transcription factor,, confers improved drought tolerance and decreased ABA sensitivity in rice., 2017, 490: 1355–1361.

[39] Guo H, Wu T, Li S, He Q, Yang Z, Zhang W, Gan Y, Sun P, Xiang G, Zhang H, Deng H. The methylation patterns and transcriptional responses to chilling stress at the seedling stage in rice., 2019, 20: 5089–5106.

[40] Sagar M, Chervin C, Mila I, Hao Y, Roustan J, Benichou M, Gibon Y, Biais B, Maury P, Latche A, Pech J C, Bouzayen M, Zouine M. SlARF4, an auxin response factor involved in the control of sugar metabolism during tomato fruit development., 2013, 161: 1362–1374.

Genome-wide identification and expression analysis of G2-like transcription factors family genes in wheat

JIA Yu-Ku**, GAO Hong-Huan**, FENG Jian-Chao, HAO Zi-Rui, WANG Chen-Yang, XIE Ying-Xin, GUO Tian-Cai, and MA Dong-Yun*

College of Agronomy, Henan Agricultural University / National Engineering Research Center for Wheat / Technology Innovation Center of Henan Wheat, Zhengzhou 450046, Henan, China

() transcription factor, a member of the GARP superfamily of MYB transcription factors, plays an important role in regulating chloroplast development. In this study, genome-wide identification ofgenes in wheat was carried out by bioinformatics methods, and their physicochemical properties, subcellular localization, cis-acting elements of promoters, and response patterns to abiotic stresses and hormones were analyzed. A total of 87genes were identified from wheat, which distributed in evenly on 21 chromosomes in wheat. Phylogenetic analysis showed that these genes were divided into 14 subfamilies, and fragment replication was the main reason for the expansion of this gene family. The prediction of protein secondary structure revealed that α helix and random curl were the main amino acid sequences ofgene in wheat. Promoter-acting elements showed that there were seven-acting elements (P-box, SpI, LTR, ABRE, MBS, TGA-Element, and AE-box) in 2-kb region upstream of the promoter. Among them,contained the most-regulatory element binding sites with a total of 18 binding sites. The qRT-PCR revealed. that the relative expression levels of,,, andwere significantly up-regulated under PEG and salt stresses, and induced by GA, IAA, and ABA hormones. These genes may mediate the response of wheat plant to various abiotic stresses.

wheat;transcription factors; bioinformatics; adversitystress; the relative expression level

10.3724/SP.J.1006.2023.21036

本研究由國家科技支撐計劃項目(2015BAD26B00)和河南省科技攻關(guān)項目(212102110281)資助。

This study was supported by the National Science and Technology Support Program of China (2015BAD26B00) and the Science and Technology Project of Henan Province (212102110281).

馬冬云, E-mail: xmzxmdy@126.com

**同等貢獻(Contributed equally to this work)

賈玉庫, E-mail: 13633969649@163.com; 高宏歡, E-mail: 2972620298@qq.com

2022-05-13;

2022-07-22;

2022-08-18.

URL: https://kns.cnki.net/kcms/detail/11.1809.S.20220817.1833.011.html

This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).