劉志敏，劉文獻(xiàn)，賈喜濤，張正社，王彥榮

(草地農(nóng)業(yè)生態(tài)系統(tǒng)國家重點(diǎn)實(shí)驗(yàn)室，蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院，甘肅 蘭州 730020)

蒺藜苜蓿LEA基因家族全基因組分析

劉志敏，劉文獻(xiàn)，賈喜濤，張正社，王彥榮

(草地農(nóng)業(yè)生態(tài)系統(tǒng)國家重點(diǎn)實(shí)驗(yàn)室，蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院，甘肅 蘭州 730020)

胚胎發(fā)育晚期豐富蛋白(LEA)是一類在種子胚胎發(fā)育后期特異表達(dá)并受發(fā)育階段及脫水信號(hào)調(diào)節(jié)的脫水保護(hù)蛋白，在響應(yīng)植物干旱、低溫、高鹽等逆境脅迫中具有重要功能。本研究利用生物信息學(xué)手段，在全基因組水平對蒺藜苜蓿(Medicagotruncatula)LEA基因家族進(jìn)行分析，并對其進(jìn)化、基因結(jié)構(gòu)、進(jìn)化壓力、染色體定位及基因表達(dá)模式進(jìn)行了系統(tǒng)分析。本研究共篩選出23個(gè)蒺藜苜蓿LEA家族基因，可分為8個(gè)亞家族?；蚨ㄎ唤Y(jié)果表明，23個(gè)蒺藜苜蓿LEA基因分布于除6號(hào)染色體外的其他7條染色體上，但分布不均勻；家族成員外顯子數(shù)目都不超過兩個(gè)，結(jié)構(gòu)簡單。不同組織表達(dá)譜分析結(jié)果顯示，LEA家族基因具有不同的組織表達(dá)模式，并受干旱逆境脅迫調(diào)控。本研究結(jié)果可為蒺藜苜蓿LEA基因家族的功能分析奠定基礎(chǔ)。

蒺藜苜蓿；LEA基因家族；系統(tǒng)進(jìn)化；表達(dá)分析

胚胎發(fā)育晚期豐富蛋白(Late Embryogenesis Abundant Proteins，LEA)最早發(fā)現(xiàn)于棉花(Gossypiumhirsutum)子葉中，因其在植物胚胎發(fā)育后期大量表達(dá)而得名[1-2]。LEA基因家族分布廣泛，不僅存在于被子植物、裸子植物和藻類中[3-4]，在線蟲(Aphelenchusavenae)[5]、輪足蟲(Rotifer)[6]和藍(lán)細(xì)菌(Cyanobacteria)[7]中也普遍存在。根據(jù)蛋白質(zhì)序列同源性及其特殊的基元序列，可將LEA蛋白分為不同類群[8-11]。研究發(fā)現(xiàn)，不同LEA蛋白都具有高親水性和熱穩(wěn)定性，并富含甘氨酸或其他氨基酸類(如丙氨酸、絲氨酸和蘇氨酸)[4]。LEA基因家族系統(tǒng)進(jìn)化分析最先在擬南芥(Arabidopsisthaliana)中被報(bào)道，其共有51個(gè)成員[12]，隨后在水稻(Oryzasativa)[13]、大豆(Glycinemax)[14]和楊樹(Populus)[15]中也發(fā)現(xiàn)分別含有34、35和40個(gè)LEA基因。

植物L(fēng)EA蛋白及其轉(zhuǎn)錄本累積水平主要在種子胚胎發(fā)育后期、種子脫水之前顯著提高；當(dāng)植物處于滲透脅迫、脫水以及低溫脅迫條件下，該類蛋白也會(huì)在植物營養(yǎng)組織中高水平累積，與植物耐脫水性密切相關(guān)，具有保護(hù)植物組織免受水分脅迫傷害的功能[16-19]。前人研究結(jié)果顯示，LEA基因家族受多種逆境脅迫調(diào)控表達(dá)[20]。例如，在干旱和低溫脅迫條件下，青楊(P.cathayana)LEA2、LEA3基因在葉片中的表達(dá)量顯著提高[21]。在干旱條件下，擬南芥營養(yǎng)組織中LEA蛋白累積量升高[22]；當(dāng)擬南芥AtLEA4基因在擬南芥中過表達(dá)后，其抗旱性相比野生型顯著提高，AtLEA4基因表達(dá)量降低則導(dǎo)致擬南芥對干旱脅迫更為敏感[23]。應(yīng)用RNA干擾技術(shù)誘導(dǎo)擬南芥LEA14及其同源脫水蛋白(LEA34和LEA51)表達(dá)量降低后，種子貯存壽命顯著降低，說明LEA蛋白對種子貯存具有十分重要的作用[24]。小麥(Triticumaestivum)TaLEA2和TaLEA3分別在酵母細(xì)胞中過表達(dá)后，相比之對照，轉(zhuǎn)化酵母細(xì)胞的抗干旱和耐低溫脅迫能力都顯著增強(qiáng)[25]。在中華羊草(Leymuschinensis)中過表達(dá)TaLEA3后，轉(zhuǎn)基因植株的抗旱性顯著提高[26]。小麥的PMA1959(LEA1)和PMA80(LEA2)基因在水稻中過表達(dá)可明顯提高轉(zhuǎn)基因植株在干旱脅迫下的葉綠素含量和根系的重量[27]。另外，在草莓(Fragaria)中過表達(dá)小麥LEA(WCOR410)基因后能顯著提高草莓葉片的抗寒性[28]，而將大麥(Hordeumvulgare)LEA(HVA1)分別在水稻和小麥中過表達(dá)后能同時(shí)提高轉(zhuǎn)基因植物的抗旱性和耐鹽性[29-30]。

蒺藜苜蓿(Medicagotruncatula)是紫花苜蓿(M.sativa)的近緣種，具有生長周期短、自花授粉、基因組小且遺傳轉(zhuǎn)化效率高等特點(diǎn)，使其成為研究豆科，特別是苜蓿屬植物，基因組學(xué)和生物學(xué)的模式植物[31]。蒺藜苜蓿全基因組測序已經(jīng)完成，使得利用生物信息學(xué)手段研究該物種基因家族系統(tǒng)演化及其功能成為可能，同時(shí)借鑒蒺藜苜蓿研究及其基因資源對促進(jìn)豆類作物遺傳研究具有重要意義[32]。截至目前，在全基因組水平對蒺藜苜蓿LEA基因家族的研究尚未見報(bào)道。本研究利用生物信息學(xué)手段，從全基因組水平對蒺藜苜蓿LEA基因家族的結(jié)構(gòu)特征、進(jìn)化關(guān)系、染色體位置及表達(dá)模式進(jìn)行系統(tǒng)分析，以期為進(jìn)一步研究LEA蛋白家族功能和豆科牧草作物遺傳改良奠定基礎(chǔ)。

1 材料與方法

1.1 材料

蒺藜苜?；蚪M序列下載于JCVI數(shù)據(jù)庫(http://jcvi.org/medicago/index.php)，擬南芥LEA基因序列和蛋白質(zhì)序列下載于Phytozome數(shù)據(jù)庫(http://phytozome.jgi.doe.gov/pz/portal.html)[33]。蒺藜苜蓿芯片探針數(shù)據(jù)來自Affymetrix公司(http://www.affymetrix.com/analysis/index.affx)[34]，所得探針于諾貝爾基金會(huì)(Noble Foundation)數(shù)據(jù)庫(http://mtgea.noble.org/v3/)[35-36]搜索LEA基因在各器官中的表達(dá)量，進(jìn)行基因表達(dá)分析。

1.2 蒺藜苜蓿LEA基因家族成員鑒定

以蒺藜苜蓿全基因組序列構(gòu)建本地BLAST數(shù)據(jù)庫，利用已鑒定出的51個(gè)擬南芥LEA基因家族成員執(zhí)行本地BLAST(P-value=le-003)搜索同源序列。利用Pfam數(shù)據(jù)庫[37]對獲得的蒺藜苜蓿同源LEA基因進(jìn)行蛋白結(jié)構(gòu)預(yù)測，篩選含有LEA典型結(jié)構(gòu)域的蛋白序列。利用Expasy蛋白質(zhì)組數(shù)據(jù)庫(http://web.expasy.org/protparam/)對蒺藜苜蓿LEA蛋白進(jìn)行蛋白質(zhì)分子量、理論等電點(diǎn)和總平均親水性預(yù)測[38]。

1.3 蒺藜苜蓿LEA基因染色體定位

利用Phytozome植物全基因組數(shù)據(jù)庫(http://phytozome.jgi.doe.gov/pz/portal.html)確定蒺藜苜蓿LEA基因染色體位置，利用Mapchart 2.2工具[39]進(jìn)行染色體定位作圖。

1.4 蒺藜苜蓿LEA基因家族成員系統(tǒng)進(jìn)化樹構(gòu)建及基因結(jié)構(gòu)分析

利用MEGA 6.0[40]對蒺藜苜蓿LEA蛋白進(jìn)行多序列聯(lián)配比對分析，采用鄰接法構(gòu)建系統(tǒng)進(jìn)化樹。系統(tǒng)進(jìn)化樹選用p-distance模型，校驗(yàn)參數(shù)Bootstrap設(shè)為1 000次重復(fù)。利用在線工具GSDS 2.0(http://gsds.cbi.pku.edu.cn/)[41]分析LEA基因家族成員基因結(jié)構(gòu)。

1.5 蒺藜苜蓿LEA基因進(jìn)化壓力分析

利用在線工具Clustal Omega[42]和PAL2NAL[43]計(jì)算核苷酸的同義替換率Ks(synonymous substitution rate or synonymous substitutions per synonymous site)、非同義替換率Ka(nonsynonymous substitution rate or nonsynonymous substitutions per nonsynonymous site)及兩者比率Ka/Ks。

1.6 蒺藜苜蓿LEA基因表達(dá)分析

利用Affymetrix基因探針數(shù)據(jù)庫搜索蒺藜苜蓿LEA基因相匹配的探針，從蒺藜苜?；蛐酒脚_(tái)(http://mtgea.noble.org/v3/)獲得對應(yīng)探針表達(dá)量信息，并使用MeV4[44]軟件對蒺藜苜蓿LEA基因在不同組織、不同時(shí)期、不同處理的表達(dá)芯片數(shù)據(jù)進(jìn)行數(shù)據(jù)聚類分析和作圖。

2 結(jié)果與分析

2.1 蒺藜苜蓿LEA基因家族的鑒定和分類

通過BLAST比對分析，本研究共初步篩選出32個(gè)候選蒺藜苜蓿LEA基因。將所得序列在Pfam網(wǎng)站進(jìn)行基因結(jié)構(gòu)域分析并刪除不含LEA結(jié)構(gòu)域的基因，最終得到23個(gè)蒺藜苜蓿LEA基因。根據(jù)結(jié)構(gòu)域可將23個(gè)LEA基因分為8個(gè)亞家族(表1)，包括4個(gè)LEA1、2個(gè)LEA2、2個(gè)LEA3、5個(gè)LEA4、1個(gè)LEA5、2個(gè)LEA6、4個(gè)SMP和3個(gè)Dehydrin。利用在線工具ExPASy對蒺藜苜蓿LEA蛋白質(zhì)分子量、理論等電點(diǎn)系數(shù)和總平均親水性進(jìn)行分析發(fā)現(xiàn)，23個(gè)LEA蛋白中分子量最大為55 384.9 D(LEA4-5)，最小為6 740.2 D(LEA4-4)；等電點(diǎn)的范圍從4.54(LEA7-4)到9.55(LEA1-4)不等；所有LEA蛋白除LEA2-2外均表現(xiàn)為親水性(表1)。

2.2 蒺藜苜蓿LEA基因家族染色體定位

蒺藜苜蓿LEA基因染色體定位結(jié)果顯示，23個(gè)基因分布于除第6條染色體以外的其他7條染色體上，但分布不均一(圖1)。其中，第2、4和7染色體上分布最多，各有5個(gè)LEA基因；其次為第1和第3染色體，各3個(gè)；第5和第8染色體上最少，各1個(gè)。LEA1-2/LEA1-3/LEA1-4、LEA4-2/LEA4-3、LEA7-2/LEA7-3及LEA8-1/LEA8-2在染色體上連鎖在一起，屬于旁系同源基因，其中LEA7-2/LEA7-3的步長值最高(100)。

2.3 蒺藜苜蓿LEA家族基因系統(tǒng)進(jìn)化及其基因結(jié)構(gòu)分析

為研究蒺藜苜蓿LEA基因系統(tǒng)進(jìn)化關(guān)系，對篩選出的蒺藜苜蓿LEA蛋白進(jìn)行了系統(tǒng)進(jìn)化樹構(gòu)建(圖2)。分析結(jié)果顯示，23個(gè)蒺藜苜蓿LEA基因分布在6個(gè)無根系統(tǒng)進(jìn)化樹，其中LEA1、LEA5和LEA8三個(gè)亞家族分布在同一分支上，并明顯地分為3個(gè)次亞家族；同時(shí)，其他5個(gè)亞家族也可明顯分為2個(gè)次亞家族且含有不同基因數(shù)。

利用GSDS在線工具對LEA基因結(jié)構(gòu)進(jìn)行作圖分析。結(jié)果顯示，蒺藜苜蓿LEA基因結(jié)構(gòu)相對簡單，所含基因內(nèi)含子的數(shù)量均不多于兩個(gè)。其中LEA6-1、LEA6-2和LEA8-3不含有內(nèi)含子，LEA4-4、LEA7-1和LEA7-3均各含有兩個(gè)內(nèi)含子，剩余17個(gè)基因只含有1個(gè)內(nèi)含子(圖3)。

2.4 蒺藜苜蓿LEA基因進(jìn)化選擇壓力分析

同義替換(Synonymous Substitution)和非同義替換(Nonsynonymous Substitution)是生物進(jìn)化過程中核苷酸替換的兩種形式，Ka、Ks及Ka/Ks能反映出基因進(jìn)化過程中受到的選擇壓力。對6對旁系同源基因的分析結(jié)果顯示(表2)，6對旁系同源基因Ka/Ks值在0.041 9到0.614 8之間，均小于1。其中，Ka/Ks最大值為LEA4-2和LEA4-3(0.614 8)；Ka/Ks最小值為LEA3-1和LEA3-2(0.041 9)，表明蒺藜苜蓿LEA基因進(jìn)化時(shí)存在負(fù)向選擇作用[45]。

表1 蒺藜苜蓿LEA基因家族信息

2.5 蒺藜苜蓿LEA基因表達(dá)模式分析

利用MeV軟件對蒺藜苜蓿LEA家族基因表達(dá)模式進(jìn)行分析，綠-黑-紅3個(gè)顏色表示基因的表達(dá)強(qiáng)度。綠色越強(qiáng)代表表達(dá)量越小，紅色越強(qiáng)代表表達(dá)量越大，黑色為中間值。本研究對蒺藜苜蓿LEA基因在不同組織(葉、葉柄、莖、花、果莢、種衣和根)、不同時(shí)期及不同水分脅迫處理的31個(gè)表達(dá)芯片數(shù)據(jù)進(jìn)行了分析。

研究結(jié)果顯示，除LEA4-3僅在成熟種子中表達(dá)外，其他LEA基因在種子和種衣中都有較高的表達(dá)量，并具有隨著種子成熟度的增加表達(dá)量呈增高的趨勢(圖4)。其中，LEA1-2、LEA2-1、LEA2-2、LEA3-1、LEA4-4、LEA7-4、LEA8-2和LEA8-3在葉、葉柄、莖、花、果莢、種衣、根和種子中組成型表達(dá)，但LEA1-2在種衣和根,LEA3-1在果莢以及LEA7-4在根中表達(dá)量相對較低。而其他15個(gè)LEA基因僅在種衣和種子中特異表達(dá)。

圖1 蒺藜苜蓿LEA基因染色體定位Fig.1 The chromosome location of the LEA genes in Medicago truncatula

圖2 蒺藜苜蓿LEA基因系統(tǒng)進(jìn)化樹Fig.2 The unrooted dendrogram of all Medicago truncatula LEA genes

對水分脅迫處理的幼苗和根進(jìn)行基因表達(dá)芯片數(shù)據(jù)分析的結(jié)果顯示，LEA5、LEA7-1、LEA7-2、LEA7-3、LEA1-1、LEA3-1、LEA2-1、LEA2-2、LEA8-2、LEA8-3、LEA4-3和LEA4-4表達(dá)量均無明顯變化，LEA4-2、LEA1-3、LEA4-5、LEA7-4和LEA8-1基因在幼苗和根中的表達(dá)量隨著干旱脅迫的加強(qiáng)表達(dá)量升高，恢復(fù)灌水后表達(dá)量急劇下降，LEA1-1、LEA1-4、LEA1-2、LEA4-1、LEA3-2、LEA6-1和LEA6-2僅在根中的表達(dá)受水分脅迫的調(diào)控，恢復(fù)灌水后呈明顯的下降趨勢；而LEA3-1在幼苗和根中隨著干旱脅迫處理和恢復(fù)灌水后，基因表達(dá)量表現(xiàn)為先下降后升高的趨勢。綜合分析蒺藜苜蓿LEA基因在不同組織和不同脅迫處理后的表達(dá)量表明，LEA5、LEA7-1、LEA7-2和LEA7-3在根和幼苗中表達(dá)量極低，卻不受干旱脅迫的調(diào)控，部分基因，如LEA2-1、LEA4-4、LEA2-2、LEA8-2和LEA8-3表達(dá)量無明顯變化；而部分基因表達(dá)量變化較大，例如在幼苗期干旱處理和根的形態(tài)建成過程中水分脅迫時(shí)，LEA3-2、LEA1-4、LEA6-1、LEA6-2、LEA4-2、LEA1-3、LEA4-5、LEA7-4和LEA8-1的表達(dá)量均增加，但LEA4-1和LEA1-1只在根形態(tài)建成過程中受水分脅迫時(shí)表達(dá)量增加而在幼苗期干旱處理中表達(dá)量減低；此外，LEA4-3在幼苗期干旱處理下特異表達(dá)，表明這些基因可能在植株幼苗期對干旱生境的抗逆響應(yīng)具有重要作用。

圖3 蒺藜苜蓿LEA家族基因結(jié)構(gòu)Fig.3 The gene structure of LEA gene family in Medicago truncatula

注：黃色框、黑色線條、藍(lán)色框分別表示外顯子、內(nèi)含子和非編碼區(qū)。

Note：The yellow boxes, black lines and blue boxes represent the exon, intron and untranslated region, respectively.

表2 蒺藜苜蓿LEA基因核苷酸替換率

3 討論與結(jié)論

隨著植物基因組研究的深入，利用比較基因組學(xué)方法分析和預(yù)測基因功能已成為目前植物基因功能研究的熱點(diǎn)之一。目前，在高等植物中已有的相關(guān)研究包括番茄(Lycopersiconesculentum)的LBD基因家族[46]、蒺藜苜蓿WRKY功能基因[47]以及大豆[14]、擬南芥[12]LEA蛋白家族分析等，對其種類、數(shù)目和進(jìn)化特點(diǎn)的研究為各基因家族功能分析奠定了理論基礎(chǔ)。LEA蛋白是植物體內(nèi)廣泛存在的一系列蛋白質(zhì)，在植物抵御干旱、凍害等非生物脅迫方面發(fā)揮著重要作用，但蒺藜苜蓿LEA基因家族及其功能分析的研究尚未見報(bào)道。本研究通過生物信息學(xué)方法分析LEA家族基因，共鑒定出8個(gè)亞家族23個(gè)蒺藜苜蓿LEA蛋白家族基因。在擬南芥中，已鑒定出9個(gè)亞家族51個(gè)基因[12]，而水稻LEA蛋白家族基因則被分為7個(gè)亞家族34個(gè)基因[13]。在大豆中，已鑒定的LEA基因分為8個(gè)亞家族35個(gè)基因[14]，與蒺藜苜蓿中的基因分類數(shù)目一致。另外，LEA基因結(jié)構(gòu)分析表明，多數(shù)蒺藜苜蓿LEA基因只含有一個(gè)或不含有內(nèi)含子，只有3個(gè)LEA基因含兩個(gè)內(nèi)含子。大豆LEA基因中也只含有較少的內(nèi)含子，其中，只有兩個(gè)基因含有多個(gè)內(nèi)含子，其他都只含有一個(gè)或者無內(nèi)含子[14]，表明LEA基因家族的基因結(jié)構(gòu)相對較簡單。

圖4 蒺藜苜蓿LEA基因表達(dá)模式Fig.4 The expression profiles of Medicago truncatula LEA genes

注：圖中用綠-黑-紅三色代表基因表達(dá)水平。綠色越亮表達(dá)越弱，紅色越亮表達(dá)越強(qiáng)。

Note：The gene expression level is displayed with green-black-red scheme in the picture. The brighter the green, the weaker the expression; the brighter the red, the stronger the expression.

蒺藜苜蓿中LEA家族基因分布差異較大，其中LEA4基因數(shù)最多，含有5個(gè)；而LEA5基因數(shù)最少，僅有1個(gè)?；蚨ㄎ粩?shù)據(jù)顯示，第6條染色體上無LEA基因，最多的第2、4和7條染色體上各含有5個(gè)LEA基因。染色體定位和系統(tǒng)進(jìn)化分析表明，有6對旁系同源基因于染色體上緊密連鎖在一起，可能是基因組串聯(lián)復(fù)制的結(jié)果[48]。另外，通過系統(tǒng)進(jìn)化樹分析得出，LEA基因在進(jìn)化上表現(xiàn)出差異性，暗示其基因功能的多樣性[49-50]。

根據(jù)蒺藜苜蓿LEA基因表達(dá)信息的分析結(jié)果發(fā)現(xiàn)，該LEA家族基因在組織表達(dá)上具有一定的特異性，且與耐旱性密切相關(guān)。15個(gè)LEA基因在種子發(fā)育時(shí)期特異表達(dá)，在種皮中也有大量表達(dá)，說明LEA基因家族對植物的種子發(fā)育具有較大影響；另外5個(gè)基因(LEA2-1、LEA2-2、LEA4-4、LEA8-2和LEA8-3)則在蒺藜苜蓿所有的組織中都表達(dá)，而LEA1-3和LEA1-4在葉和花中的表達(dá)量相較之其他非種子的表達(dá)量較高，從而使得該家族的基因功能表顯示出一定的多樣性和復(fù)雜性。有研究表明，LEA基因受干旱脅迫調(diào)控表達(dá)。例如，水稻OsLEA19a基因在幼苗中不表達(dá)，當(dāng)受到干旱和高鹽處理后，其表達(dá)量快速增加，并在處理6 h后達(dá)到最高值[51]。與此對應(yīng)，在干旱脅迫下，13個(gè)基因(LEA1-2、LEA1-3、LEA1-4、LEA3-1、LEA3-2、LEA4-1、LEA4-2、LEA4-3、LEA4-5、LEA6-1、LEA6-2、LEA7-4和LEA8-1)有表達(dá)上調(diào)和下調(diào)的變化，說明其可能在植物響應(yīng)干旱脅迫中具有重要作用。

運(yùn)用生物信息學(xué)手段，對蒺藜苜蓿LEA基因家族全基因組進(jìn)行分析，共鑒定出23個(gè)家族基因。23個(gè)家族基因可分為8個(gè)亞家族，并分布在除第6條染色體外的其他7條染色體上，且其基因中內(nèi)含子數(shù)目較少，結(jié)構(gòu)簡單?；蚪M織表達(dá)譜分析結(jié)果顯示，LEA家族基因組織表達(dá)模式差異較大，并受干旱逆境脅迫調(diào)控。

[1] Dure L,Greenway S C,Galau G A.Developmental biochemistry of cottonseed embryogenesis and germination-changing messenger ribonucleic-acid populations as shown by in vitro and in vivo protein-synthesis[J].Biochemistry,1981,20(14):4162-4168.

[2] Galau G A,Hughes D W,Dure L.Abscisic-acid induction of cloned cotton late embryogenesis-abundant (LEA) messengerrnas[J].Plant Molecular Biology,1986,7(3):155-170.

[3] Campos F,Zamudio F,Covarrubias A A.Two different late embryogenesis abundant proteins fromArabidopsisthalianacontain specific domains that inhibitEscherichiacoligrowth[J].Biochemical and Biophysical Research Communications,2006,342(2):406-413.

[4] Velten J,Oliver M J.Tr288:A rehydrin with a dehydrin twist[J].Plant Molecular Biology,2001,45(6):713-722.

[5] Solomon A,Salomon R,Paperna I,Glazer I.Desiccation stress of entomopathogenic nematodes induces the accumulation of a novel heat-stable protein[J].Parasitology,2000,121(4):409-416.

[6] Tunnacliffe A,Lapinski J,Mcgee B.A putative LEA protein,but no trehalose,is present in anhydrobiotic bdelloid rotifers[J].Hydrobiologia,2005,546:315-321.

[7] Shinozaki K,Yamaguchi-Shinozaki K.Gene networks involved in drought stress response and tolerance[J].Journal of Experimental Botany,2007,58(2):221-227.

[8] Bray E A.Molecular responses to water-deficit[J].Plant Physiology,1993,103(4):1035-1040.

[9] Wise M J.LEAping to conclusions:A computational reanalysis of late embryogenesis abundant proteins and their possible roles[J].BMC Bioinformatics,2003,4:52.

[10] Dure L,Crouch M,Harada J,Ho T,Mundy J,Quatrano R,Thomas T,Sung Z R.Common amino-acid sequence domains among the lea proteins of higher-plants[J].Plant Molecular Biology,1989,12(5):475-486.

[11] Dure L.A repeating 11-mer amino-acid motif and plant desiccation[J].Plant Journal,1993,3(3):363-369.

[12] Hundertmark M,Hincha D K.LEA (Late Embryogenesis Abundant) proteins and their encoding genes inArabidopsisthaliana[J].BMC Genomics,2008,9(1):118.

[13] Wang X,Zhu H,Jin G,Liu H,Wu W,Zhu J.Genome-scale identification and analysis ofLEAgenes in rice (OryzasativaL.)[J].Plant Science,2007,172(2):414-420.

[14] 李樂,許紅亮,楊興露,李雅軒,胡英考.大豆LEA基因家族全基因組鑒定、分類和表達(dá)[J].中國農(nóng)業(yè)科學(xué),2011,44(19):3945-3954.

[15] Lan T,Gao J,Zeng Q.Genome-wide analysis of the LEA (late embryogenesis abundant) protein gene family inPopulustrichocarpa[J].Tree Genetics & Genomes,2013,9(1):253-264.

[16] Steponkus P L,Uemura M,Joseph R A,Gilmour S J,Thomashow M F.Mode of action of theCOR15agene on the freezing tolerance ofArabidopsisthaliana[J].Proceedings of the National Academy of Sciences of the United States of America,1998,95(24):14570-14575.

[17] Grelet J,Benamar A,Teyssier E,Avelange-Macherel M H,Grunwald D,Macherel D.Identification in pea seed mitochondria of a late-embryogenesis abundant protein able to protect enzymes from drying[J].Plant Physiology,2005,137(1):157-167.

[18] Skriver K,Mundy J.Gene expression in response to abscisic acid and osmotic stress[J].Plant Cell,1990,2(6):503-512.

[19] Close T J.Dehydrins:A commonality in the response of plants to dehydration and low temperature[J].Physiologia Plantarum,1997,100(2):291-296.

[20] 王康,朱慧森,董寬虎.植物L(fēng)EA蛋白及其基因家族成員PM16研究進(jìn)展[J].草業(yè)科學(xué),2008,25(2):97-102.

[21] 魯松.青楊(PopuluscathayanaRehd.)對逆境脅迫的反應(yīng)及l(fā)ea基因的表達(dá)分析[D].成都：中國科學(xué)院研究生院,2009.

[22] Bies-Etheve N,Gaubier-Comella P,Debures A,Lasserre E,Jobet E,Raynal M,Cooke R,Delsenyet M.Inventory,evolution and expression profiling diversity of the LEA (late embryogenesis abundant) protein gene family inArabidopsisthaliana[J].Plant Molecular Biology,2008,67(1-2):107-124.

[23] Olvera-Carrillo Y,Campos F,Reyes J L,Garciarrubio A,Covarrubias A A.Functional analysis of the group 4 late embryogenesis abundant proteins reveals their relevance in the adaptive response during water deficit inArabidopsis[J].Plant Physiology,2010,154(1):373-390.

[24] Hundertmark M,Buitink J,Leprince O,Hincha D K.The reduction of seed-specific dehydrins reduces seed longevity inArabidopsisthaliana[J].Seed Science Research,2011,21(3):165-173.

[25] Yu J N,Zhang J S,Shan L,Chen S Y.Two new group 3LEAgenes of wheat and their functional analysis in yeast[J].Journal of Integrative Plant Biology,2005,47(11):1372-1381.

[26] Wang L,Li X,Chen S,Liu G.Enhanced drought tolerance in transgenicLeymuschinensisplants with constitutively expressed wheatTaLEA3[J].Biotechnology Letters,2009,31(2):313-319.

[27] Cheng Z Q,Targolli J,Huang X Q,Wu R.WheatLEAgenes,PMA80 andPMA1959,enhance dehydration tolerance of transgenic rice (OryzasativaL.)[J].Molecular Breeding,2002,10(1-2):71-82.

[28] Houde M,Dallaire S,N'Dong D,Sarhan F.Overexpression of the acidic dehydrinWCOR410 improves freezing tolerance in transgenic strawberry leaves[J].Plant Biotechnology Journal,2004,2(5):381-387.

[29] Sivamani E,Bahieldin A,Wraith J M,Al-Niemi T,Dyer W E,Ho T,Qu R D.Improved biomass productivity and water use efficiency under water deficit conditions in transgenic wheat constitutively expressing the barleyHVA1 gene[J].Plant Science,2000,155(1):1-9.

[30] Xu D P,Duan X L,Wang B Y,Hong B M,Ho T,Wu R.Expression of a late embryogenesis abundant protein gene,HVA1,from barley confers tolerance to water deficit and salt stress in transgenic rice[J].Plant Physiology,1996,110(1):249-257.

[31] 魏臻武,蓋鈞鎰.豆科模式植物——蒺藜苜蓿[J].草業(yè)學(xué)報(bào),2008,17(1):114-120.

[32] 屠德鵬,魏臻武,武自念,雷艷芳,張棟,邱偉偉.蒺藜苜蓿EST-SSRs分布特征及標(biāo)記的開發(fā)[J].草業(yè)科學(xué),2011,28(5):746-752.

[33] Goodstein D M,Shu S,Howson R,Neupane R,Hayes R D,Fazo J,Mitros T,Dirks W,Hellsten U,Putnam N,Rokhsar D S.Phytozome:A comparative platform for green plant genomics[J].Nucleic Acids Research,2012,40(D1):D1178-D1186.

[34] Liu G,Loraine A E,Shigeta R,Cline M,Cheng J,Valmeekam V,Sun S,Kulp D,Siani-Rose M A.NetAffx:Affymetrix probesets and annotations[J].Nucleic Acids Research,2003,31(1):82-86.

[35] He J,Benedito V A,Wang M,Murray J D,Zhao P X,Tang Y,Udvardi M K.TheMedicagotruncatulagene expression atlas web server[J].BMC Bioinformatics,2009,10:441.

[36] Benedito V A,Torres-Jerez I,Murray J D,Andriankaja A,Allen S,Kakar K,Wandrey M,Verdier J,Zuber H,Ott T,Moreau S,Niebel A,Frickey T,Weiller G,He J,Dai X,Zhao P X,Tang Y,Udvardi M K.A gene expression atlas of the model legumeMedicagotruncatula[J].The Plant Journal,2008,55(3):504-513.

[37] Finn R D,Mistry J,Schuster-Bockler B,Griffiths-Jones S,Hollich V,Lassmann T,Moxon S,Marshall M,Khanna A,Durbin R,Eddy S R,Sonnhammer E L,Bateman A.Pfam:Clans,web tools and services[J].Nucleic Acids Research,2006,34:D247-D251.

[38] Artimo P,Jonnalagedda M,Arnold K,Baratin D,Csardi G,de Castro E,Duvaud S,Flegel V,Fortier A,Gasteiger E,Grosdidier A,Hernandez C,Ioannidis V,Kuznetsov D,Liechti R,Moretti S,Mostaguir K,Redaschi N,Rossier G,Xenarios I,Stockinger H.ExPASy:SIB bioinformatics resource portal[J].Nucleic Acids Research,2012,40:W597-W603.

[39] Voorrips R E.MapChart:Software for the graphical presentation of linkage maps and QTLs[J].The Journal of Heredity,2002,93(1):77-78.

[40] Tamura K,Stecher G,Peterson D,Filipski A,Kumar S.MEGA6:Molecular evolutionary genetics analysis version 6.0[J].Molecular Biology and Evolution,2013,30(12):2725-2729.

[41] 郭安源,朱其慧,陳新,羅靜初.GSDS:基因結(jié)構(gòu)顯示系統(tǒng)[J].遺傳,2007,29(8):1023-1026.

[42] Sievers F,Wilm A,Dineen D,Gibson T J,Karplus K,Li W,Lopez R,Mcwilliam H,Remmert M,Soding J,Thompson J D,Higgins D G.Fast,scalable generation of high-quality protein multiple sequence alignments using Clustal Omega[J].Molecular Systems Biology,2011,7:539.

[43] Suyama M,Torrents D,Bork P.PAL2NAL:Robust conversion of protein sequence alignments into the corresponding codon alignments[J].Nucleic Acids Research,2006,34:W609-W612.

[44] Saeed A I,Bhagabati N K,Braisted J C,Liang W,Sharov V,Howe E A,Li J,Thiagarajan M,White J A,Quackenbush J.TM4 microarray software suite[J].Methods in Enzymology,2006,411:134-193.

[45] Yang Z,Bielawski J P.Statistical methods for detecting molecular adaptation[J].Trends in Ecology & Evolution,2000,15(12):496-503.

[46] 王小非,劉鑫,蘇玲,孫永江,張世忠,郝玉金,由春香.番茄LBD基因家族的全基因組序列鑒定及其進(jìn)化和表達(dá)分析[J].中國農(nóng)業(yè)科學(xué),2013,46(12):2501-2513.

[47] 江騰,林勇祥,劉雪,江海洋,朱蘇文.苜蓿全基因組WRKY轉(zhuǎn)錄因子基因的分析[J].草業(yè)學(xué)報(bào),2011,20(3):211-218.

[48] Nei M.Gene duplication and nucleotide substitution in evolution[J].Nature,1969,221:5175-5177.

[49] 周麗,姬向楠,何非,潘秋紅,段長青.植物胚胎發(fā)育晚期豐富蛋白的結(jié)構(gòu)與功能[J].熱帶生物學(xué)報(bào),2012,3(2):191-196.

[50] 陳雷,李磊,李金花,李永春.植物L(fēng)EA蛋白及其功能[J].中國農(nóng)學(xué)通報(bào),2009(24):143-146.

[51] 胡廷章,吳應(yīng)梅,陳再剛,黃小云.水稻OsLEA19a基因的克隆、表達(dá)及生物信息學(xué)分析[J].華北農(nóng)學(xué)報(bào),2011,26(6):73-78.

(責(zé)任編輯 武艷培)

Genome-wide analysis ofLEAgene family inMedicagotruncatula

LIU Zhi-min, LIU Wen-xian, JIA Xi-tao, ZHANG Zheng-she, WANG Yan-rong

(State Key Laboratory of Grassland Agro-ecosystems, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020, China)

Late embryogenesis abundant (LEA) family proteins are expressed abundantly at the late stage of embryonic development of plant seeds, and play important roles in responding to various environmental stress, such as drought, low temperature and high salinity. At present, the phylogenetic studies ofLEAgene family inMedicagotruncatulahave not been reported yet. In this study, with the application of bioinformatics methods, theLEAfamily genes ofM.truncatulawere identified through the whole genome, and the system evolution, gene structure, evolutionary pressure, chromosomal location and gene expression patterns were further analyzed. A total of 23LEAgenes were systematically identified fromM.truncatulaand classified into 8 subfamilies. Gene location results showed that these 23LEAgenes were distributed unevenly on 7 chromosomes; The exon numbers of all the genes were no more than two, indicating a simple gene structure of this gene family. The expression profiles ofM.truncatulaLEAgenes showed a characteristic of temporal and tissue specific, and regulated by drought stress. This study provides a theoretical foundation for verifying the function ofLEAgenes inM.truncatula.

Medicagotruncatula;LEAgene family; system evolution; expression analysis

LIU Wen-xian E-mail: liuwx@lzu.edu.cn

10.11829j.issn.1001-0629.2014-0418

2014-09-16 接受日期：2014-12-11

國家草品種區(qū)域試驗(yàn)項(xiàng)目——苜蓿品種“三性”測試(農(nóng)財(cái)發(fā)[2013]49號(hào))；長江學(xué)者和創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃(IRT13019)；蘭州大學(xué)中央高?；A(chǔ)研究基金項(xiàng)目(lzujbky-2013-82)

劉志敏(1990-)，男，江蘇贛榆人，在讀碩士生，主要從事草類種質(zhì)資源與育種研究。E-mail:liuzhm2013@lzu.edu.cn

劉文獻(xiàn)(1981-)，男，河南開封人，副教授，碩導(dǎo)，博士，主要從事草類種質(zhì)資源和分子生物學(xué)研究。E-mail:liuwx@lzu.edu.cn

S551+.703;Q943.2

A

1001-0629(2015)03-0382-10

劉志敏，劉文獻(xiàn)，賈喜濤，張正社，王彥榮.蒺藜苜蓿LEA基因家族全基因組分析[J].草業(yè)科學(xué),2015,32(3):382-391.

LIU Zhi-min,LIU Wen-xian,JIA Xi-tao,ZHANG Zheng-she,WANG Yan-rong.Genome-wide analysis ofLEAgene family inMedicagotruncatula[J].Pratacultural Science，2015,32(3)：382-391.