劉晶，王俊杰，程璐，陳利青，禾璐，候蕊，韓淵懷,*

(1.山西農業(yè)大學 農學院，山西 太谷 030801； 2.山西農業(yè)大學 生物工程研究所，山西 太谷 030801；3.山西省農業(yè)科學院玉米研究所，山西 忻州 034000； 4.山西省晉城市城區(qū)農業(yè)委員會，山西 晉城 048400)

?

谷子Aux/IAA家族基因干旱脅迫響應表達模式分析

劉晶1，王俊杰2，程璐1，陳利青1，禾璐3，候蕊4，韓淵懷1,2*

(1.山西農業(yè)大學 農學院，山西 太谷 030801； 2.山西農業(yè)大學 生物工程研究所，山西 太谷 030801；3.山西省農業(yè)科學院玉米研究所，山西 忻州 034000； 4.山西省晉城市城區(qū)農業(yè)委員會，山西 晉城 048400)

[目的]Aux/IAA是重要的蛋白質轉錄因子，廣泛參與植物生長素介導的應答作用，同時參與植物的脅迫和防御響應。本文通過分析Aux/IAA家族基因在干旱脅迫中的表達情況，探索谷子抗旱與Aux/IAA家族基因的關系。[方法]本文對谷子Aux/IAA家族基因理化性質、蛋白親疏水性、啟動子功能元件等進行生物信息學分析，并對抗旱(GG)和干旱敏感(JF16)谷子品種在澆水和干旱脅迫下基因的表達譜數(shù)據(jù)進行分析。[結果]谷子Aux/IAA家族基因均含有多個與脅迫相關的功能元件，且5個基因編碼的蛋白均為親水性蛋白；干旱脅迫處理后，GG和JF16中5個基因的表達變化趨勢有所不同，Seita.5G126200、Seita.1G028400、Seita.3G109400表達量均下調，Seita.1G356800基因的表達量明顯上調。[結論]谷子Aux/IAA家族基因參與調控谷子干旱脅迫，其調控過程比較復雜，可作為參與谷子抗旱的候選基因。本研究結果為進一步探究谷子中Aux/IAA與抗旱機制提供一定理論依據(jù)。

谷子；Aux/IAA； 干旱脅迫； 生物信息學分析

谷子(Setariaitalica(L.)Beauv.)，在古代被稱做“百谷之長”，是我國種植最古老的農作物之一，是重要的雜糧作物[1]。我國的谷子種植范圍較廣，主要分布在華北、東北地區(qū)。由于種植地區(qū)不同，各地區(qū)的環(huán)境、氣候條件差異較大，對谷子的生長發(fā)育產(chǎn)生較大影響。

目前，影響谷子產(chǎn)量的眾多因素中，干旱(drought)是限制谷子產(chǎn)量提高的主要因素。因此，了解谷子干旱脅迫條件下的分子機制對提高谷子抗旱性有很大作用。

生長素IAA(Indole-3-acetic acid)是重要的植物激素之一，對植物生長發(fā)育具有重要的調節(jié)作用，其在轉錄水平上的調控受很多蛋白的影響。Aux/IAA(auxin /indoleacetic acids protein)和生長素應答因子ARF(auxin response factor)是2個重要的蛋白質轉錄因子，它們廣泛參與植物生長素介導的應答作用， Aux/IAA蛋白是通過與ARF轉錄因子相互作用介導生長素反應的短壽命轉錄調控因子。通常Aux/IAA家族基因均包含4個保守結構域：結構域 I、II、III、IV[2]。

近年來對Aux/IAA基因家族的研究逐漸成為熱點，雖然對植物生長素的研究主要集中在植物發(fā)育過程中，但最近的一些研究表明，生長素也參與植物的脅迫或防御反應。據(jù)報道，植物感病時內源性激素IAA含量顯著增加[3]，同時，一些生長素調節(jié)基因的表達量也會發(fā)生改變[4]；microRNA介導的生長素信號被抑制時，可以使細菌的耐藥性增強[5]；對水稻生長素相關的基因進行全基因組分析發(fā)現(xiàn)，幾種與生長素相關的基因在各種非生物脅迫條件下的表達存在差異，說明生長素與非生物脅迫信號之間有聯(lián)系[6]；擬南芥處于低溫環(huán)境時，其中Aux/IAA和ARF基因家族成員的表達也會發(fā)生改變[7]。

目前國內外對于谷子中Aux/IAA基因家族干旱脅迫表達分析的研究較少，本研究以耐旱品種勾勾母雞咀(GG)和干旱敏感品種晉汾16(JF16)為材料，種植21 d后，用20%的PEG-6000對谷子幼苗采取干旱脅迫處理，將處理后的幼苗提取RNA，進行表達譜測序分析，找到與脅迫相關的Aux/IAA家族基因，對其理化性質、啟動子元件等進行生物信息學分析，同時對干旱脅迫下耐旱品種GG和干旱敏感品種JF16中不同Aux/IAA家族基因成員的表達量進行分析，為進一步研究谷子Aux/IAA家族基因在干旱脅迫過程中的功能奠定基礎，為闡明谷子在分子水平的抗旱機制提供一定的理論依據(jù)。

1 材料與方法

1.1 試驗材料

試驗所用材料耐旱品種勾勾母雞咀(GG)和干旱敏感品種晉汾16(JF16)，由山西農業(yè)大學農業(yè)生物工程研究所和山西省農業(yè)科學院品種資源研究所提供，在山西農業(yè)大學農學院人工氣候室種植培養(yǎng)。

將GG和JF16種于混合土壤(蛭石∶營養(yǎng)土=1∶3)，在14 h光照(28 ℃)、10 h黑暗(23 ℃)條件下培養(yǎng)21 d后，用20%的PEG-6000對谷子幼苗進行30 min干旱脅迫處理，對照組用蒸餾水處理。將處理后的幼苗隨機選取3株，提取RNA、反轉錄，然后進行表達譜測序分析。

1.2 試驗方法

利用ExPASY(https://web.expasy.org/compute-pi/)分析谷子Aux/IAA家族基因的理化性質；

利用GSDS 2.0(http://gsds.cbi.pku.edu.cn/)分析基因結構；

利用Expasy Protscale(http://web.expasy.org/protscale/)分析Aux/IAA家族基因編碼蛋白的親水/疏水性；

利用Phytozome 12(https://phytozome.jgi.doe.gov/pz/portal.html)谷子基因組數(shù)據(jù)庫查找Aux/IAA家族基因的基因序列信息；

利用PlantCARE(http://bioinformatics.psb.ugent.be/webtools/plantcare/html/)對谷子Aux/IAA家族基因進行啟動子元件分析；

利用MAGA7.0軟件對谷子中找到的5個Aux/IAA家族基因與玉米、高梁等禾本科作物的Aux/IAA基因進行比對，構建系統(tǒng)進化樹。

2 結果與分析

2.1 谷子Aux/IAA家族基因的基因結構與理化性質

在谷子中找到了5個Aux/IAA家族基因，對它們進行基因結構和理化性質的分析，主要包括基因組序列長度、CDS序列長度、在染色體上的位置、氨基酸數(shù)量、等電點和分子量。由表1可見，谷子Aux/IAA家族的5個基因分布在3條不同的染色體上，分別位于1，3，5號染色體的不同位點。他們的基因組序列長度有很大區(qū)別，最小為973 bp，最大的有4 822 bp；有4個蛋白長度在231～277個氨基酸之間，最小(Seita.1G028400)為175個氨基酸，相對應的蛋白分子量在18 914.19(Seita.1G356800)～29 316.79(Seita.1G028400)Da之間，CDS序列長度在528～834 bp之間。等電點分析發(fā)現(xiàn)，除Seita.1G356800等電點的為6.08顯酸性外，其余4個均大于7.5，顯堿性。

表1 谷子Aux/IAA 家族基因的基本信息

2.2Aux/IAA家族基因內含子-外顯子結構組成分析

使用GSDS2.0對Aux/IAA家族基因的CDS序列、基因序列進行結構分析。由圖1可見所有基因都包含內含子和外顯子結構，Seita.1G028400、Seita.3G109400都含有5個內含子，Seita.5G268600含有4個內含子，而Seita.1G356800和Seita.5G126200只含有2個內含子；Seita.5G126200不含有上下游序列，Seita.1G028400、Seita.3G109400有較長的上游序列。

圖1 谷子Aux/IAA 家族基因內含子-外顯子結構Fig.1 Intron-exon structures of Aux/IAA family members in foxtail millet

2.3Aux/IAA家族基因編碼蛋白的親水/疏水性分析

用Expasy Protscale分析Aux/IAA家族基因編碼蛋白的親水/疏水性，結果如圖2所示，Aux/IAA家族基因蛋白氨基酸殘基中疏水性最強的為Seita.1G356800(2.500)，最弱的為Seita.3G109400(2.360)；親水性最強的為Seita.5G126200(-3.780)，最弱的為Seita.3G109400(-3.080)。在5個Aux/IAA家族基因編碼的蛋白中，親水性氨基酸數(shù)量明顯多于疏水性氨基酸的數(shù)目，故推斷5個Aux/IAA家族基因編碼的蛋白均為親水蛋白。

2.4Aux/IAA家族基因的啟動子元件及功能分析

使用在線軟件PlantCARE對所取谷子的Aux/IAA家族基因啟動子上游1 500 bp進行啟動子元件分析。由表2可見，谷子Aux/IAA家族基因都含有至少一個或多個與脅迫相關的調控元件，它們共同調控谷子生長過程對脅迫的應對。5個基因均含有與茉莉酸甲酯(MeJA)相關的調控元件，除Seita.5G126200外均含有與脫落酸(ABA)相關的調控元件，除Seita.1G356800外都含有與水楊酸(SA)相關的調控元件，Seita.1G356800和Seita.1G028400含有與赤霉素(GA)相關的調控元件等。

圖2 Aux/IAA家族基因蛋白的親水/疏水性預測Fig.2 Prediction of hydrophobicity and rophilicity of Aux/IAA family members

2.5Aux/IAA家族基因成員干旱脅迫時在GG和JF中的表達量分析

使用20%的PEG-6000對谷子幼苗進行干旱脅迫處理后，對Aux/IAA家族基因在耐旱品種GG和干旱敏感品種JF16中的表達量進行分析。

由圖3可見，GG和JF16中，Seita.5G126200、Seita.1G028400、Seita.3G109400三個基因在干旱脅迫條件下，表達量都下調；Seita.5G268600基因在正常澆水時，2個品種相比表達水平差異不明顯，在干旱脅迫時，耐旱品種GG中表達量上調，在干旱敏感品種JF16中表達量卻出現(xiàn)下調；Seita.1G356800基因在正常澆水情況下，在GG中表達水平很低，在JF16中表達水平較高；在干旱脅迫時，GG和JF16中Seita.1G356800基因的表達量都明顯上調。

表2 谷子Aux/IAA家族基因啟動子功能元件數(shù)目

圖3 谷子Aux/IAA家族基因干旱脅迫及對照情況下在GG和JF中的表達量Fig.3 The expression of Aux/IAA family members in GG and JF under drought stress and control

2.6 谷子Aux/IAA家族基因系統(tǒng)進化樹

使用MAGA7.0軟件對谷子中找到的5個Aux/IAA家族基因與玉米、高梁等禾本科作物的Aux/IAA基因進行比對，構建系統(tǒng)進化樹，可以直觀的看到5個成員與玉其他植物Aux/IAA基因親緣關系的遠近，了解其進化過程。由圖4可見，Aux/IAA家族5個基因中，Seita.5G268600與Zeamays5親緣關系最近，Seita.5G126200與Zeamays15處在同一分支上，親緣關系較近，而Seita.1G356800與Zeamays9親緣關系較近，且Seita.5G268600、Seita.5G126200、Seita.1G356800處在同一大的分支，表明這3個基因親緣關系較近。Seita.1G028400與Seita.3G109400處在一個大的分支，分別與Zeamays2和Zeamays16聚為一類，說明二者親緣關系較其他3個基因近。

圖4 谷子Aux/IAA家族基因進化樹分析Fig.4 Phylogenetic tree of Aux/IAA family members

3 結論與討論

生長素IAA是植物體內一種非常重要的激素，在植物的生長發(fā)育過程中起重要的調節(jié)作用。Aux/IAA基因是生長素信號轉導過程中重要的調節(jié)轉錄因子，同時它在植物抗逆性方面也起到調節(jié)作用。本試驗通過對谷子幼苗進行20%的PEG-6000干旱脅迫處理，找到5個響應干旱脅迫的Aux/IAA家族基因，對其基因結構、理化性質、與脅迫相關的啟動子原件等進行分析，為研究谷子Aux/IAA家族基因在干旱脅迫等發(fā)育進程中的功能奠定基礎。

通過對Aux/IAA家族5個基因編碼的蛋白質進行親水/疏水性分析，發(fā)現(xiàn)5個蛋白均為親水性蛋白。蛋白質的高親水性使植物即使受到干旱脅迫，蛋白質也能夠代替細胞中的水分子，從而很好的抵御干旱脅迫對植物造成的損傷，在提高植物的抗逆性方面有很大作用[8]。表明，Aux/IAA家族基因對谷子抗旱性的提高方面有影響。

通過對谷子Aux/IAA家族基因啟動子上游1 500 bp進行啟動子元件分析，結果顯示，谷子Aux/IAA家族基因包含許多與脅迫相關的元件：ABRE是與脫落酸(ABA)相關的調控元件，而ABA可以提高植物的抗旱能力；CGTCA-motif和TGACG-motif均是與茉莉酸甲酯(MeJA)相關的調控元件，MeJA能夠激發(fā)植物防御基因的表達，從而提高植物的抗逆能力[9]；TCA-element和TATC-box分別是與水楊酸(SA)和赤霉素(GA)相關的調控元件，據(jù)報道，干旱條件下在玉米幼苗葉表皮噴灑SA可以提高其抗旱性、GA可以提高煙草幼苗及種子的抗旱性[10,11]；MBS(TAACTG)是與干旱誘導相關的調控元件；MBS(CGGTCA)、TC-rich repeats等都是與脅迫相關的調控元件，它們共同作用來調節(jié)谷子生長過程對脅迫的響應。

有研究表明，對不同脅迫下水稻中Aux/IAA家族基因的表達譜進行分析，多個基因表達量出現(xiàn)上調或下調[12]。本試驗通過分析Aux/IAA家族的5個基因在耐旱品種GG和干旱敏感品種JF16中的表達量，發(fā)現(xiàn)在干旱脅迫時，Seita.1G356800基因表達量在2個品種中都上調，說明該基因響應干旱脅迫，提高植物的抗旱能力，與Jain等人的研究結果一致[6]；Seita.5G268600基因干旱脅迫下在耐旱品種GG中表達量上調，在干旱敏感品種JF中表達量下調，說明該基因在2個不同品種中表達響應不同，可能是造成2個品種抗旱能力不同的原因之一；Seita.5G126200、Seita.1G028400、Seita.3G109400三個基因表達量均下調，說明Aux/IAA家族基因在生長素、水楊酸和抗生素脅迫信號通路之間的串擾中起到重要作用[13]。

綜上所述，本文對谷子Aux/IAA家族基因的理化性質、蛋白質親疏水性、所含脅迫相關元件以及系統(tǒng)進化樹進行了生物信息學分析，基于不同品種谷子干旱脅迫下的表達譜數(shù)據(jù)，初步探究了谷子Aux/IAA家族基因在干旱脅迫下的表達量及變化規(guī)律，為研究谷子抗旱分子機制提供了一定的理論依據(jù)。

[1]禾璐．綠色小米品種大青谷米色形成機制研究[D]．太谷：山西農業(yè)大學，2015．

[2]Hagen G，Guilfoyle T．Auxin-responsive gene expression: genes,promoters and regulatory factors[J]．Plant Molecular Biology，2002，49(3 /4)：373-385.

[3]O’Donnell P J，Schmelz E A，Moussatche P，etal．Susceptible to intolerance a range of hormonal actions in a susceptible Arabidopsis pathogen response[J]．Plant Journal，2003，33(2)：245-257．

[4]Dowd C，Wilson I W，McFadden H．Gene expression profile changes in cotton root and hypocotyl tissues in response to infection with Fusarium oxysporum f. sp. vasinfectum[J]．Molecular Plant-Microbe Interactions，2004，17(6)：654-667．

[5]Navarro L，Dunoyer P，Jay F，et al．A plant miRNA contributes to antibacterial resistance by repressing auxin signaling[J]．Science，2006，312(5772)：436-439．

[6]Jain M，Khurana J P．Transcript profiling reveals diverse roles of auxin-responsive genes during reproductive development and abiotic stress in rice[J]．Febs Journal，2009，276(11)：3148-3162．

[7]Hannah M A，Heyer A G，Hincha D K．A global survey of gene regulation during cold acclimation in Arabidopsis thaliana[J]．Plos Genetics，2005，1(2)：e26．

[8]Baker J，Steele C，Dure L．Sequence and characterization of 6 Lea proteins and their genes from cotton[J]．Plant Molecular Biology，1988，11(3)：277-291．

[9]Horbowicz M，Kosson R，Wiczkowski W，et al．The effect of methyl jasmonate on accumulation of 2-phenylethylamine and putrescine in seedlings of common buckwheat (Fagopyrum esculentum)[J]．Acta Physiol Plantarum，2011，33(3)：897-903．

[10]坎平，王莎莎，馬文廣，等. 赤霉素引發(fā)同時提高煙草種子及幼苗抗旱性和抗冷性[J]. 種子，2014，33(2)：30-34，38．

[11]曹翠玲，劉林麗，田強兵．水楊酸對玉米幼苗抗旱性的影響[J]．玉米科學，2004，12(增刊)：103-104．

[12]Song Y，Wang L，Xiong L．Comprehensive expression profiling analysis of OsIAAgene family in developmental processes and in response to phytohormone and stress treatments[J]．Planta，2009，229(3)：577-591．

[13]Cakir B，Kilickaya O，Olcay A C．Genome-wide analysis ofAux/IAAgenes in Vitis vinifera: cloning and expression profiling of a grapeAux/IAAgene in response to phytohormone and abiotic stresses[J]．Acta Physiologiae Plantarum，2013，35(2)：365-377．

(編輯：韓志強)

Expression analysis ofAux/IAAfamily members in response to drought stress in foxtail millet

Liu Jing1, Wang Junjie2, Cheng Lu1, Chen Liqing1, He Lu3, Hou Rui4, Han Yuanhuai1,2*

(1.CollegeofAgriculture，ShanxiAgriculturalUniversity，Taigu030801，China; 2.InstituteofAgriculturalBioengineering，ShanxiAgriculturalUniversity，Taigu030801，China; 3.MaizeResearchInstitutetoShanxiAcademyofAgricultureSciences，Xinzhou034000，China; 4.UrbanAgricultureCommitteeofShanxiJincheng，Jincheng048400，China)

[Objective]Aux/IAA is an important transcription factor protein. It is widely involved in plant auxin-mediated response, plant stress and defense response. The relationship between drought resistance andAux/IAAfamily members was studied by analyzing the expression ofAux/IAAfamily members under drought stress in foxtail millet.[Methods]The bioinformatics analysis was carried out on the physicochemical properties, protein hydrophobicity and cis-acting elements of theAux/IAAfamily members in foxtail millet. Then we analyzed the expression ofAux/IAAfamily members under drought stress.[Results]The result showed that there were a number of stress-related functional elements in theAux/IAAfamily members in foxtail millet, and the proteins they code were all hydrophilic. The expression levels of genes were significantly different between drought-tolerant and drought-sensitive varieties under drought stress. The expression levels ofSeita1.5G126200,Seita.1G028400 andSeita.3G109400 were down-regulated, whereas the expression level ofSeita.1G356800 gene was significantly up-regulated.[Conclusion]TheAux/IAAfamily members in foxtail millet are involved in the regulation of millet drought stress, and the process of regulation is a complex process. Some of them could also be used as candidate genes for drought tolerance in foxtail millet. Theseresult provide a theoretical basis for further study of the mechanism of the relationship between theAux/IAAand drought tolerance in foxtail millet.

Foxtail millet,Aux/IAAfamily members, Drought stress, Bioinformatics analysis

2017-04-21

2017-05-22

劉晶(1994-)，女(漢)，山西晉中人，碩士，研究方向：基因組學與分子育種

*通信作者：韓淵懷，教授，博士生導師，Tel：0354-6287239；E-mail：swgctd@163.com

國家自然科學基金(31401396)；山西農業(yè)大學科技創(chuàng)新基金(20142-07)

S515

A

1671-8151(2017)09-0609-07