郭寶健，李 贏，袁澤宸，呂 超，張新忠，許如根

(江蘇省作物遺傳生理國(guó)家重點(diǎn)實(shí)驗(yàn)室培育點(diǎn)/糧食作物現(xiàn)代產(chǎn)業(yè)技術(shù)協(xié)同創(chuàng)新中心/教育部植物功能基因組學(xué)重點(diǎn)實(shí)驗(yàn)室/揚(yáng)州大學(xué)大麥研究所，江蘇揚(yáng)州 225009)

大麥ARF基因家族的全基因組分析

郭寶健，李 贏，袁澤宸，呂 超，張新忠，許如根

(江蘇省作物遺傳生理國(guó)家重點(diǎn)實(shí)驗(yàn)室培育點(diǎn)/糧食作物現(xiàn)代產(chǎn)業(yè)技術(shù)協(xié)同創(chuàng)新中心/教育部植物功能基因組學(xué)重點(diǎn)實(shí)驗(yàn)室/揚(yáng)州大學(xué)大麥研究所，江蘇揚(yáng)州 225009)

生長(zhǎng)素響應(yīng)因子(auxin response factor,ARF)基因家族是植物特有的轉(zhuǎn)錄因子家族，在調(diào)控植物生長(zhǎng)發(fā)育過(guò)程中起到重要作用。而關(guān)于大麥ARF基因家族全基因組分析的研究尚未見(jiàn)報(bào)道。為進(jìn)一步探討大麥ARF基因的功能，以公布的大麥栽培品種Morex的基因組數(shù)據(jù)為基礎(chǔ)，采用生物信息學(xué)方法鑒定了大麥ARF基因，并對(duì)其結(jié)構(gòu)、染色體分布、蛋白保守結(jié)構(gòu)域、系統(tǒng)進(jìn)化樹(shù)及表達(dá)譜進(jìn)行了分析。結(jié)果表明，共鑒定出了17個(gè)大麥ARF基因，除4號(hào)染色體外，其余6條染色體上均有分布，片段復(fù)制事件是大麥ARF基因家族的主要擴(kuò)張方式；HvARF蛋白均具有保守的B3結(jié)構(gòu)域和Auxin_resp結(jié)構(gòu)域，而Aux/IAA結(jié)構(gòu)域僅存在于12個(gè)HvARF蛋白中，且不同蛋白所含該結(jié)構(gòu)域1～2個(gè)不等；不同植物中ARF基因在功能上具有保守性；不同組織器官中HvARF基因的表達(dá)存在明顯的差異。

大麥；生長(zhǎng)素響應(yīng)因子；進(jìn)化樹(shù)；表達(dá)譜

生長(zhǎng)素響應(yīng)因子(auxin response factor,ARF)基因家族是植物特有的轉(zhuǎn)錄因子家族，調(diào)控生長(zhǎng)素響應(yīng)的基因表達(dá)[1-3]。絕大多數(shù)的ARF蛋白具有3個(gè)保守的結(jié)構(gòu)域，即N-末端B3 DNA結(jié)合域(DNA-binding domain,DBD)，具有DNA結(jié)合活性，能直接與生長(zhǎng)素調(diào)控的下游基因啟動(dòng)子區(qū)的AuxREs結(jié)合；中間區(qū)域(Auxin_resp結(jié)構(gòu)域)，具有激活結(jié)構(gòu)域(AD)或抑制結(jié)構(gòu)域(RD)活性，具有激活或抑制轉(zhuǎn)錄子活性的作用；C末端二聚化結(jié)構(gòu)域(CTD)，即AUX/IAA結(jié)構(gòu)域，與Aux/IAAs蛋白C-末端的結(jié)合域Ⅲ和Ⅳ高度保守，介導(dǎo)Aux/IAAs和ARFs的同源二聚化和異源二聚化的區(qū)域，調(diào)節(jié)ARFs與AuxREs的結(jié)合并調(diào)控下游生長(zhǎng)素響應(yīng)基因的表達(dá)[2,4-5]。

目前，植物中已克隆到大量的ARF基因，部分基因的功能也得到鑒定[3,6-7]。 AtARF5(MONOPTEROS)是植物中第一個(gè)通過(guò)圖位克隆方法鑒定的ARF基因，該基因在胚模式建成、維管組織正常形成過(guò)程中起重要作用[6]。 AtARF7和 AtARF19是報(bào)道最多的ARF基因。擬南芥中， arf7突變體植株主要表現(xiàn)為下胚軸的向光性和對(duì)生長(zhǎng)素運(yùn)輸抑制劑NPA的抗性， arf19突變體植株相對(duì)于野生型的表型沒(méi)有明顯的變化， arf7arf19雙突變除了表現(xiàn)出單突變的表型外，側(cè)根數(shù)目嚴(yán)重減少，即 AtARF7與 AtARF19在功能上有一定的冗余，在擬南芥?zhèn)雀l(fā)生過(guò)程中起重要作用[7-9]。 AtIAA14和 AtLBD16/29與 AtARF7/19共同參與了擬南芥的側(cè)根發(fā)育的調(diào)控網(wǎng)絡(luò)[10-11]。Attia等[12]研究表明，降低 OsARF1基因的表達(dá)量，轉(zhuǎn)基因水稻株系表現(xiàn)出種子活力低、發(fā)育遲緩、葉片短而卷曲和不育，表明 OsARF1對(duì)營(yíng)養(yǎng)器官和生殖器官的發(fā)育都起著重要作用。

許多植物的ARF基因家族已得到系統(tǒng)分析，如擬南芥中有23個(gè)ARF基因、玉米中有31個(gè)ARF基因、水稻中有25個(gè)ARF基因、谷子中有24個(gè)ARF基因[7,13-17]。大麥?zhǔn)鞘澜缟系谒拇蠛坦阮愖魑?，具有廣泛的適應(yīng)性和用途，栽培歷史悠久。2012年，大麥全基因組測(cè)序完成，為大麥分子生物學(xué)研究奠定了良好的基礎(chǔ)[18]。目前，大麥中還沒(méi)有ARF基因家族分析的報(bào)道，也缺少對(duì)該基因家族的生物信息學(xué)分析和表達(dá)分析，因此有必要在大麥基因組測(cè)序完成的基礎(chǔ)上對(duì)ARF基因家族進(jìn)行系統(tǒng)分析。本研究利用生物信息學(xué)方法在全基因組水平上鑒定大麥ARF基因，并對(duì)其基因結(jié)構(gòu)、染色體分布、蛋白保守結(jié)構(gòu)域、分子進(jìn)化及表達(dá)譜進(jìn)行分析，以期為進(jìn)一步探討大麥ARF基因的功能奠定基礎(chǔ)。

1 材料與方法

1.1 大麥ARF基因家族成員的鑒定及其編碼蛋白的理化性質(zhì)分析

利用國(guó)際大麥基因組測(cè)序聯(lián)盟公布的大麥基因組數(shù)據(jù)[18]，以已報(bào)道的擬南芥[7]、水稻[13]和玉米[14-16]的ARF蛋白序列作為引用數(shù)據(jù)，在IPK (http://planttfdb.cbi.pku.edu.cn/)和Gramene((http://gramene.org/)數(shù)據(jù)庫(kù)中進(jìn)行BLASTP比對(duì)(e-value<1e-5)，搜索這些序列的同源蛋白，刪除重復(fù)序列。大麥ARF候選基因的蛋白質(zhì)序列利用在線軟件Pfam (http://www.sanger.ac.uk/Software/Pfam/search.shtml)進(jìn)行保守結(jié)構(gòu)域分析[19]，篩選含有 B3結(jié)構(gòu)域、Auxin_resp結(jié)構(gòu)域和AUX_IAA 結(jié)構(gòu)域的序列，剔除冗余蛋白。借助Expasy (http://www.expasy.org/)在線軟件計(jì)算大麥ARF蛋白的等電點(diǎn)和分子量[20]。

1.2 大麥ARF基因的結(jié)構(gòu)、染色體定位和啟動(dòng)子序列分析

從Gramene數(shù)據(jù)庫(kù)中檢索大麥ARF基因的信息，包括內(nèi)含子、外顯子、染色體上的物理位置信息和啟動(dòng)子序列。采用GSDS (http://gsds.cbi.pku.edu.cn/)分析大麥ARF基因的內(nèi)含子和外顯子結(jié)構(gòu)[21]。利用MapDraw軟件構(gòu)建染色體定位圖[22]。按照大麥ARF基因在染色體上的排列順序進(jìn)行基因命名。在PlantCARE數(shù)據(jù)庫(kù)(http://bioinformatics.psb.ugent.be/webtools/ plantcare/html/)中輸入大麥ARF基因的啟動(dòng)子序列，預(yù)測(cè)生長(zhǎng)素響應(yīng)元件[24]。

1.3 分子進(jìn)化樹(shù)的構(gòu)建

利用Muscle軟件對(duì)大麥、玉米、水稻和擬南芥ARF蛋白序列進(jìn)行多序列比對(duì)，將比對(duì)結(jié)果輸入MEGA 6.0，采用鄰接法(neighbor-joining,NJ) 分別構(gòu)建大麥ARF基因家族進(jìn)化樹(shù)及大麥、玉米、水稻和擬南芥的系統(tǒng)進(jìn)化樹(shù)，其中，校驗(yàn)參數(shù)(bootstrap)設(shè)置為1 000，其余均為默認(rèn)參數(shù)[25]。片段復(fù)制事件參照Maher等[23]的方法，查找旁系同源基因上游和下游各10個(gè)編碼蛋白的基因，如在ARF基因的兩側(cè)還存在旁系同源基因?qū)?，則表明該ARF旁系同源基因?qū)ζ鹪从谄螐?fù)制事件。

1.4 大麥ARF基因的表達(dá)譜分析

利用已公布的大麥RNA-seq數(shù)據(jù)，檢索大麥ARF基因的表達(dá)譜(http://apex.ipk-gatersleben.de/apex/f?p=284:10:6281639160219::NO)[18]。數(shù)據(jù)庫(kù)中提供了大麥栽培品種Morex的8個(gè)組織器官的表達(dá)數(shù)據(jù)，包括種子萌發(fā)后4 d的胚、幼苗的根和嫩芽、幼穗(5 mm和1～1.5 cm)、六葉期的分蘗、開(kāi)花后5 d和15 d的籽粒(去稃殼)，獲取數(shù)據(jù)庫(kù)中8個(gè)組織器官的FPKM(fragments per kilobase of transcript per million mapped reads)值，作為組織表達(dá)譜數(shù)據(jù)，構(gòu)建基因表達(dá)熱圖。圖中用綠色和紅色表示基因表達(dá)的強(qiáng)度，綠色越亮表示基因表達(dá)量越低，紅色越亮表示基因表達(dá)量越高。

2 結(jié)果與分析

2.1 大麥ARF基因的染色體定位及理化性質(zhì)

利用已報(bào)道的31個(gè)玉米ARF基因、25個(gè)水稻ARF基因和23個(gè)擬南芥ARF基因序列在IPK (http://planttfdb.cbi.pku.edu.cn/)和Gramene (http://gramene.org/)數(shù)據(jù)庫(kù)中進(jìn)行BLASTP比對(duì)，得到35個(gè)大麥候選ARF基因。利用Pfam (http://www.sanger.ac.uk/Software/Pfam/search.shtml) 對(duì)獲得的大麥候選ARF基因編碼蛋白質(zhì)序列進(jìn)行保守結(jié)構(gòu)域分析，剔除不含有典型ARF結(jié)構(gòu)域的冗余蛋白后，獲得17個(gè)大麥ARF基因。這17個(gè)大麥ARF基因分布于除4H染色體外的其余6條染色體上，3H和7H上含有的基因數(shù)目最多，均為4個(gè)，其次為2H和6H，均為3個(gè)，1H和5H上含有的ARF基因數(shù)目最少，分別為2個(gè)和1個(gè)(圖1)。按照大麥ARF基因在染色體上的排列順序命名為 HvARF1～ HvARF17(圖1)。序列分析顯示，17個(gè)HvARF基因編碼區(qū)長(zhǎng)度在1 356～3 387 bp之間，編碼長(zhǎng)度為452～1 129個(gè)氨基酸的蛋白質(zhì)，分子量為50.35～125.38 kDa，等電點(diǎn)范圍為5.77～8.92(表1)。

2.2 大麥ARF基因的結(jié)構(gòu)和系統(tǒng)進(jìn)化樹(shù)

利用本研究中大麥ARF全長(zhǎng)蛋白序列構(gòu)建大麥系統(tǒng)進(jìn)化樹(shù)，結(jié)果(圖2A)顯示，在17個(gè)編碼大麥ARF蛋白的基因中檢測(cè)到了5個(gè)旁系同源基因?qū)?，?HvARF1和 HvARF16、 HvARF2和 HvARF7、 HvARF8和 HvARF9、 HvARF11和 HvARF15、 HvARF12和 HvARF17?；蚪Y(jié)構(gòu)分析表明，不同大麥ARF基因所含外顯子差異很大，其中， HvARF3、 HvARF5、 HvARF6、 HvARF9、 HvARF10和 HvARF11均含有14個(gè)外顯子， HvARF4、 HvARF14和 HvARF17均含有13個(gè)外顯子， HvARF8、 HvARF12和 HvARF15均含有12個(gè)外顯子，而 HvARF1、HvARF13和 HvARF16則均含有3個(gè)外顯子， HvARF2和 HvARF7則均含有2個(gè)外顯子(圖2B)。旁系同源基因?qū)Φ幕虿介L(zhǎng)值為100的基因?qū)? HvARF2和 HvARF7、 HvARF8和 HvARF9、 HvARF11和 HvARF15)含有相似數(shù)目的外顯子(圖2B)。片段復(fù)制和串聯(lián)重復(fù)是家族基因擴(kuò)張的主要模式，分析發(fā)現(xiàn)， HvARF1和 HvARF6、 HvARF4和 HvARF5基因?qū)Υ嬖谄螐?fù)制事件；在Gramene數(shù)據(jù)庫(kù)中通過(guò)Gbrowse查看成簇分布的ARF基因家族成員與相鄰蛋白質(zhì)編碼基因的位置關(guān)系，未發(fā)現(xiàn)大麥ARF基因家族成員位于同一個(gè)或相鄰的基因間區(qū)域，即大麥中ARF基因家族的擴(kuò)張不存在串聯(lián)重復(fù)方式，表明片段復(fù)制是大麥ARF基因家族擴(kuò)張的主要方式。大麥(17個(gè))、玉米(31個(gè))、水稻(25個(gè))及擬南芥(23個(gè))ARF基因的系統(tǒng)進(jìn)化樹(shù)分析表明，來(lái)自4個(gè)植物的96個(gè)ARF基因分為6個(gè)類別(ClassⅠ～Ⅵ)[26]，而本研究中得到的17個(gè)大麥ARF基因歸于ClassⅠ～ClassⅤ類別(圖3)。

圖1 大麥ARF基因的染色體定位

表1 大麥ARF家族基因的基本信息

Table 1 The basic information of ARF gene family in barley

基因名稱Genename基因號(hào)GeneID染色體Chromosome物理位置Physicalposition編碼區(qū)長(zhǎng)度Codingsequencelength/bp氨基酸長(zhǎng)度Aminoacidlength/aa分子量MW/Da等電點(diǎn)pIHvARF1MLOC_64795.11H127102527～127107560206468874246.737.59HvARF2MLOC_38232.11H419398338～419403571180060065937.558.77HvARF3MLOC_64596.12H508740133～508745733197465873058.305.81HvARF4MLOC_55345.62H589588539～5895952242871957105593.965.79HvARF5MLOC_51932.12H615547380～615555915250883692948.875.85HvARF6MLOC_65945.13H141906589～141910998228676284195.966.24HvARF7MLOC_66439.13H435834578～435839922164154759634.608.84HvARF8MLOC_68907.13H523997002～524001458224174783941.846.27HvARF9MLOC_18401.13H529062935～529067687247882692329.446.03HvARF10MLOC_63938.15H22704804～22710863267389198194.835.77HvARF11MLOC_73144.46H69851096～6985889633871129125382.045.98HvARF12MLOC_66152.26H124383735～124387564135645250346.196.90HvARF13MLOC_77438.26H247346568～247349350210970376028.868.92HvARF14MLOC_14584.17H202204080～20221199932491083120278.156.00HvARF15MLOC_63193.27H533232921～53324135633031101123570.716.20HvARF16MLOC_69988.17H545375651～545378650209169775629.948.13HvARF17MLOC_58330.27H561567257～5615740832814938103351.366.14

圖2 大麥ARF基因家族的系統(tǒng)進(jìn)化樹(shù)(A)和基因結(jié)構(gòu)(B)

圖3 大麥、玉米、水稻和擬南芥ARF基因的進(jìn)化樹(shù)分析

2.3 大麥ARF基因的啟動(dòng)子元件預(yù)測(cè)結(jié)果

利用PlantCARE軟件對(duì)HvARF基因起始密碼子上游2 kb的DNA序列分析表明，6個(gè)HvARF基因含有生長(zhǎng)素響應(yīng)元件(TGA-element，AACGAC)，其中， HvARF16含有兩個(gè)生長(zhǎng)素響應(yīng)元件，分別位于啟動(dòng)子正義鏈(1 130 bp)和反義鏈(-1 341 bp)； HvARF1、 HvARF7和 HvARF13基因的生長(zhǎng)素響應(yīng)元件均位于正義鏈(1 523 bp、1 340 bp、1 249 bp)； HvARF2和 HvARF15基因的生長(zhǎng)素響應(yīng)元件均位于反義鏈(-1 427 bp、-1 887 bp)。

2.4 大麥ARF蛋白的保守結(jié)構(gòu)域

利用Pfam軟件分析了大麥ARF蛋白的結(jié)構(gòu)域，結(jié)果(圖4)表明，所有大麥ARF蛋白都含有N-端DNA結(jié)合域(B3結(jié)構(gòu)域)，并具有激活或抑制活性的中間區(qū)域(Auxin_resp結(jié)構(gòu)域)，12個(gè)ARF蛋白含有C末端二聚化結(jié)構(gòu)域(AUX/IAA結(jié)構(gòu)域)，其中HvARF8和HvARF9分別含有2個(gè)AUX/IAA結(jié)構(gòu)域。

圖4 大麥ARF蛋白的保守結(jié)構(gòu)域

2.5 大麥ARF基因的表達(dá)譜

栽培品種Morex 8個(gè)組織器官的RNA-seq數(shù)據(jù)分析結(jié)果顯示，17個(gè)HvARF基因中有12個(gè)在8個(gè)組織器官中檢測(cè)到了轉(zhuǎn)錄本。如圖5所示， HvARF3、 HvARF5、 HvARF10、 HvARF12和 HvARF17在幼穗(5 mm和1～1.5 cm)和六葉期的分蘗中具有較高的表達(dá)量，此外， HvARF3在幼苗期的根和嫩芽中也表現(xiàn)出了較高的表達(dá)量； HvARF1基因的轉(zhuǎn)錄本主要集中在六葉期的分蘗和開(kāi)花后5 d的籽粒中，且隨著成熟度的提高，在籽粒中的表達(dá)量下降； HvARF11、 HvARF15和 HvARF16基因主要在六葉期的分蘗中表達(dá)； HvARF4基因的轉(zhuǎn)錄本主要集中在幼穗中； HvARF6和 HvARF13基因?yàn)榻M成型表達(dá)基因，并且在 8個(gè)組織器官中均表現(xiàn)出低表達(dá)模式。

1：萌發(fā)后4天的胚；2：幼苗期的根；3：幼苗期的嫩芽；4：幼穗(5 mm)；5：幼穗(1～1.5 cm)；6：六葉期的分蘗；7：去稃殼籽粒(開(kāi)花后5 d)；8：去稃殼籽粒(開(kāi)花后15 d)。

1:Embryo four-day after germination;2:Roots at seedling stage;3:Shoots at seedling stage;4:Young inflorescences (5 mm);5:Young inflorescences (1～1.5 cm);6:Tillers at six-leaf stage;7:Developing grain with hulls removed (5 DPA);8:Developing grain with hulls removed (15 DPA).

圖5 大麥ARF基因在不同組織器官中的表達(dá)譜

Fig.5 Expression profile ofARFgenes in eight tissues of barley

3 討 論

本研究在大麥全基因組水平上共檢測(cè)到17個(gè)編碼ARF蛋白的基因，少于已報(bào)道的擬南芥(23個(gè))、水稻(25個(gè))、玉米(31個(gè))、谷子(24個(gè))等植物ARF基因家族成員數(shù)目[7,13-17]。片段復(fù)制是大麥ARF基因擴(kuò)張的唯一方式，可能限制了大麥ARF基因家族的擴(kuò)張。

ARF蛋白的系統(tǒng)進(jìn)化樹(shù)分析表明，4個(gè)HvARF蛋白歸于ClassⅠ，其中 HvARF5、 HvARF10和HvARF17蛋白含有B3結(jié)構(gòu)域、轉(zhuǎn)錄激活活性的Auxin_resp結(jié)構(gòu)域和AUX/IAA結(jié)構(gòu)域三個(gè)保守的結(jié)構(gòu)域，同時(shí)編碼這三個(gè)蛋白的基因具有相似的表達(dá)模式，擬南芥中ClassⅠ的ARF蛋白同時(shí)含有三個(gè)保守的結(jié)構(gòu)域，其中AtARF7和AtARF19作為轉(zhuǎn)錄激活子介導(dǎo)了擬南芥?zhèn)雀纬蛇^(guò)程[7,26]，暗示 HvARF5、 HvARF10和 HvARF17基因可能存在相似的功能。本研究中5個(gè)HvARF蛋白歸于ClassⅡ，該類基因數(shù)目最多，除含有B3結(jié)構(gòu)域和AUX/IAA結(jié)構(gòu)域外，還含有轉(zhuǎn)錄抑制活性的Auxin_resp結(jié)構(gòu)域，僅有 HvARF3、 HvARF4和 HvARF9在八個(gè)組織器官中檢測(cè)到轉(zhuǎn)錄本，且表達(dá)模式存在明顯的差異，可能存在功能分化。同類的 AtARF3和 AtARF4編碼了兩個(gè)功能冗余的生長(zhǎng)素響應(yīng)因子， arf3arf4雙突變體所有側(cè)生器官背腹結(jié)構(gòu)受到影響[27]。 HvARF2和 HvARF7基因與 AtARF3和 AtARF4基因同屬于Class Ⅲ，且 HvARF2與 AtARF3均不含有AUX/IAA結(jié)構(gòu)域[4]，表明這類基因在植物生長(zhǎng)發(fā)育過(guò)程中具有不同的功能。ClassⅣ和ClassⅤ的B3結(jié)構(gòu)中含有額外的32～36個(gè)氨基酸，其余的保守結(jié)構(gòu)域分別與ClassⅡ和ClassⅢ相似， AtARF10、 AtARF16和 AtARF17基因的單突變體沒(méi)有明顯的表型，而 arf10arf16雙突變體根冠和根系向地性生長(zhǎng)發(fā)生異常[7,28]，大麥中共有5個(gè)HvARF基因?qū)儆诖祟悇e，但是具體的功能還有待于進(jìn)一步的探討。

本研究利用生物信息學(xué)方法在全基因組水平上分析了大麥ARF基因家族，并對(duì)其基因結(jié)構(gòu)、染色體分布、蛋白保守結(jié)構(gòu)域、分子進(jìn)化及表達(dá)譜進(jìn)行了分析，為進(jìn)一步探討大麥ARF基因的功能奠定一定的基礎(chǔ)，但是關(guān)于大麥ARF基因如何調(diào)控生長(zhǎng)素響應(yīng)的基因表達(dá)參與大麥生長(zhǎng)發(fā)育還需要進(jìn)一步的探討。

[1]GUILFOYLE T J,HAGEN G.Auxin response factors [J].JournalofPlantGrowthRegulation,2001,10:281.

[2]TIWARI S B,HAGEN G,GUILFOYLE T J.The roles of auxin response factor domains in auxin-responsive transcription [J].PlantCell,2003,15(2):533.

[3]ELLIS C M,NAGPAL P,YOUNG J C,etal.AUXINRESPONSEFACTOR1 andAUXINRESPONSEFACTOR2 regulate senescence and floral organ abscission inArabidopsisthaliana[J].Development,2005,132(20):4563.

[4]GUILFOYLE T J,HAGEN G.Auxin response factors [J].CurrentOpinioninPlantBiology,2007,10(5):453.

[5]ABEL S,OELLER P W,THEOLOGIS A.Early auxin induced genes encode short-lived nuclear proteins [J].ProceedingsoftheNationalAcademyofSciences,1994,91:326.

[6]HARDTKE C S,BERLETH T.TheArabidopsisgeneMONOPTEROSencodes a transcription factor mediating embryo axis formation and vascular development [J].EMBOJournal,1998,17(5):1405.

[7]OKUSHIMA Y,OVERVOORDE P J,ARIMA K,etal.Functional genomic analysis of the AUXIN RESPONSE FACTOR gene family members inArabidopsisthaliana:unique and overlapping functions of ARF7 and ARF19 [J].PlantCell,2005,17(2):444.

[8]LISCUM E,BRIGGS W R.Mutations in theNPH1 locus ofArabidopsisdisrupt the perception of phototropic stimuli [J].PlantCell,1995,7(4):473.

[9]STOWE-EVANS E L,HARPER R M,MOTCHOULSKI A V,etal.NPH4,a conditional modulator of auxin-dependent differential growth responses inArabidopsis[J].PlantPhysiology,1998,118(4):1265.

[10]FUKAKI H,NAKAO Y,OKUSHIMA Y,etal.Tissue-specific expression of stabilized SOLITARY-ROOT/IAA14 alters lateral root development inArabidopsis[J].PlantJournal,2005,44:382.

[11]OKUSHIMA Y,FUKAKI H,ONODA M,etal. ARF7 and ARF19 regulate lateral root formation via direct activation ofLBD/ASLgenes inArabidopsis[J].PlantCell,2007,19(1):118.

[12]ATTIA K A,ABDELKHALIK A F,AMMAR M H,etal.Antisense phenotypes reveal a functional expression of OsARF1,an auxin response factor,in transgenic rice [J].CurrentIssuesinMolecularBiology,2009,11:29.

[13]WANG D,PEI K,FU Y,etal.Genome-wide analysis of the auxin response factors (ARF) gene family in rice (Oryzasativa) [J].Gene,2007,394(1):13.

[14]XING H Y,PUDAKE R N,GUO G G,etal.Genome-wide identification and expression profiling of auxin response factor(ARF) gene family in maize [J].BMCGenomics,2011,12:178.[15]LIU Y,JIANG H Y,CHEN W,etal.Genome-wide analysis of the auxin response factor(ARF) gene family in maize(Zeamays) [J].PlantGrowthRegulation,2011,63(3):225.[16]WANG Y J,DENG D X,SHI Y T,etal.Diversification,phylogeny and evolution of auxin response factor (ARF) family:insights gained from analyzing maizeARFgenes [J].MolecularBiologyReports,2012,39(3):2401.

[17] 趙 艷,甕巧云,馬海蓮,等.谷子ARF基因家族的鑒定與生物信息學(xué)分析[J].植物遺傳資源學(xué)報(bào),2016,17(3):547.

ZHAO Y,WENG Q Y,MA H L,etal.Genome-wide identification and bioinformatics analysis of ARF gene family in foxtail milletSetariaitalic[J].JournalofPlantGeneticResources,2016,17(3):547.

[18]The International Barley Genome Sequencing Consortium.A physical,genetic and functional sequence assembly of the barley genome [J].Nature,2012,491:711.

[19]PUNTA M,COGGILL P C,EBERHARDT R Y,etal.The Pfam protein families database [J].NucleicAcidsResearch,2012,40:D290.

[20]ARTIMO P,JONNALAGEDDA M,ARNOLD K,etal.Expasy:sib bioinformatics resource portal [J].NucleicAcidsResearch,2012,40(w1):W597.

[21]HU B,JIN J,GUO A Y,etal.GSDS 2.0:an upgraded gene feature visualization server [J].Bioinformatics,2014,31(8):1296.

[22]LIU R H,MENG J L.MapDraw:a microsoft excel macro for drawing genetic linkage maps based on given genetic linkage data [J].Hereditas,2003,25(3):317.

[23]MAHER C,STEIN L,WARE D.Evolution ofArabidopsismicroRNA families through duplication events [J].GenomeResearch,2006,16(4):510.

[24]LESCOT M,DéHAIS P,THIJS G,etal.PlantCARE,a database of plantcis-acting regulatory elements and a portal to tools forinsilicoanalysis of promoter sequences [J].NucleicAcidsResearch,2002,30:325.

[25]TAMURA K,STECHER G,PETERSON D,etal.MEGA6:molecular evolutionary genetics analysis version 6.0 [J].MolecularBiologyandEvolution,2013,30(4):2725.

[26]FINET C,BERNE-DEDIEU A,SCUTT C P,etal.Evolution of the ARF gene family in land plants:old domains,new tricks [J].MolecularBiologyandEvolution,2013,30(1):45.

[27]PEKKER I,ALVAREZ J P,ESHED Y.Auxin response factors mediateArabidopsisorgan asymmetry via modulation of KANADI activity [J].PlantCell,2005,17(11):2899.

[28]WANG J W,WANG L J,MAO Y B,etal.Control of root cap formation by microRNA-targeted auxin response factors inArabidopsis[J].PlantCell,2005,17(8):2204.

Genome-Wide Analysis of Auxin Response Factor (ARF) Family in Barley

GUO Baojian,LI Ying,YUAN Zechen,Lü Chao,ZHANG Xinzhong,XU Rugen

(Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops/Key Laboratory of Plant Functional Genomics of the Ministry of Education/Barley Research Institution of Yangzhou University,Yangzhou,Jiangsu 225009,China)

Auxin response factor (ARF) gene family is the specific transcription factor in plant,and plays an important role in regulation of plant organ development.However,genome-wide analysis of the ARF gene family has not been reported in barley.In this study,HvARFgenes were identified based on the public genomic database of barley cultivar Morex by using bioinformatics methods.A comprehensive overview of theHvARFgenes was undertaken,including gene structure,chromosome position,conserved domain of proteins,phylogenetic tree and gene expression profiling.The results showed thatHvARFgenes located on all chromosomes except for 4H,and segmental duplication contributed to the expansion of the ARF gene family in barley;all HvARF proteins contain a conserved B3 domain and Auxin_resp domain,and twelve HvARF proteins contain one or two Aux/IAA domains;phylogenetic classification showed that the function ofARFgenes were conserved in different species;and gene expression analysis showed that the expression levels ofHvARFgenes were obviously different at different development stages in barley.

Barley;Auxin response factor;Phylogenetic tree;Expression profiling

時(shí)間：2016-11-04

2016-05-16

2016-10-13

國(guó)家自然科學(xué)基金項(xiàng)目(31401370；31571648)；國(guó)家大麥青稞產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)(CARS-05)；江蘇省農(nóng)業(yè)科技自主創(chuàng)新基金項(xiàng)目(CX(14)2123)；江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程資助項(xiàng)目

E-mail：bjguo@yzu.edu.cn

許如根(E-mail:rgxu@yzu.edu.cn)

S512.3；S330

A

1009-1041(2016)11-1426-07

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/61.1359.S.20161104.0924.004.html