蔓花生AP2基因家族的生物信息學(xué)分析

王玲 劉曉偉 江納 付春

摘要：AP2作為一類植物轉(zhuǎn)錄因子，在花的發(fā)育和營養(yǎng)器官的形成過程中都起著非常重要的作用。旨在對蔓花生（Arachis duranensis）AP2蛋白的理化性質(zhì)及分子進(jìn)化關(guān)系等進(jìn)行詳細(xì)的生物信息學(xué)分析。結(jié)果表明，蔓花生AP2基因家族成員編碼蛋白亞細(xì)胞定位預(yù)測均在細(xì)胞核中，有11個成員的理論等電點小于7，推測該家族蛋白多為酸性蛋白，有8個成員存在2個保守結(jié)構(gòu)域;該家族成員均無信號肽;除了Aradu.SE6Q0與Aradu.42K79成員中存在跨膜現(xiàn)象外，其他成員鮮少出現(xiàn)跨膜現(xiàn)象。親疏水性預(yù)測結(jié)果表明，蔓花生AP2基因家族編碼蛋白的親疏水性系數(shù)均小于0，推測其為親水蛋白;磷酸化位點主要集中在絲氨酸上;其二級結(jié)構(gòu)以無規(guī)則卷曲和α-螺旋為主。研究還得出，蔓花生AP2基因家族編碼蛋白的主要保守基序是Motif1、Motif2和Motif4;Motif3、Motif4都含有YRG結(jié)構(gòu)域，Motif5含有RAYD結(jié)構(gòu)域;構(gòu)建的系統(tǒng)發(fā)育樹也有力地說明了AP2基因家族的進(jìn)化程度。

關(guān)鍵詞：蔓花生;AP2基因;基因家族;生物信息學(xué)

中圖分類號： S541+.901? 文獻(xiàn)標(biāo)志碼： A

文章編號：1002-1302（2020）14-0065-13

蔓花生（Arachis duranensis），別稱遍地黃金，是豆科蝶形花亞科蔓花生屬的多年生宿根草本植物，原產(chǎn)于亞洲熱帶及南美洲地區(qū)，其根系發(fā)達(dá)，抗逆性強，病蟲害少，是動物的良好牧草[1-6]。蔓花生有一定的耐熱耐旱性、抗寒性、耐陰性、耐貧瘠性[2-3，7-13]，對有害氣體有較強的抗性[14-15]。APETALA2（AP2）是一類植物轉(zhuǎn)錄因子，其結(jié)構(gòu)首次被發(fā)現(xiàn)于擬南芥中，并證實在花的發(fā)育過程中和營養(yǎng)器官中均有表達(dá)[16-20]。該類因子由60～70個氨基酸殘基組成，并參與DNA的結(jié)合過程[18，21-22]。AP2基因家族是AP[STBX]2[STBZ]/ERF超級基因家族下的1個分支，其家族蛋白包含2個重復(fù)的AP2結(jié)構(gòu)域，即AP2-R1和AP2-R2，根據(jù)2個結(jié)構(gòu)域是否含有插入序列，將AP2基因家族分為ANT組、euAP2組，此外，AP2基因家族還含有較長的識別序列GCAC（A/G）N（A/T）TCCC（A/G）ANG（C/T）元件[23-28]。AP[STBX]2/ERF超級基因家族中除AP2[STBZ]家族外，還包括ERF家族和RAV家族，ERF家族編碼蛋白含有1個AP2結(jié)構(gòu)域，RAV家族編碼蛋白含有1個AP2結(jié)構(gòu)域和1個B3結(jié)構(gòu)域，其中B3結(jié)構(gòu)域在其他植物的特異轉(zhuǎn)錄因子中是1個保守的DNA結(jié)合結(jié)構(gòu)域[21]。AP[STBX]2[STBZ]/ERF亞家族成員通過與CC、CRT/DRE、JERE或VWRE等作用元件相互作用而誘導(dǎo)下游基因的表達(dá)[29-32]。AP2參與花發(fā)育、小穗分生組織的形成、葉表皮細(xì)胞識別和胚胎發(fā)育等過程的調(diào)控[33]以及花器官的識別，如花瓣、萼片[34-38]。在AP[STBX]2[STBZ]家族編碼蛋白中，擬南芥蛋白AP213（AP2）、AP2-14（TOE1）、AP2-8（AIL6）、AP2-12（AIL7）含有EAR-like基序[17]。AP2的CAP2在鷹嘴豆（Cicer arietinum Linn.）中具有抗鹽、抗干旱脅迫的功能[39]，EgAP2-1在油棕（Elaeis guineensis Jacq.）中調(diào)控合子和體細(xì)胞胚的發(fā)育[40]。AP2/ERF轉(zhuǎn)錄因子基因家族的結(jié)構(gòu)和進(jìn)化在擬南芥、水稻和大豆中是相對保守的[41]，許多亞家族和亞類存在于這3個物種中，這說明絕大多數(shù)AP2/ERF亞家族成員存在于物種分化之前[18]。AP[STBX]2[STBZ]轉(zhuǎn)錄基因在花生種子發(fā)育的4個時期的表達(dá)存在差異[33]，其通過影響胚、胚乳和種皮的發(fā)育來控制種子的大小，同時通過對糖代謝的影響來影響種子的質(zhì)量[18，39，42]。目前尚未見關(guān)于蔓花生AP2基因家族理化特性及其進(jìn)化關(guān)系的報道，因此，本研究以蔓花生AP2蛋白序列為材料，利用生物信息學(xué)工具對蔓花生AP2基因家族的理化特性、保守基序和系統(tǒng)進(jìn)化關(guān)系進(jìn)行詳細(xì)研究，以期為進(jìn)一步研究和揭示AP2基因家族在蔓花生生長發(fā)育、基因調(diào)控中的作用提供參考依據(jù)。

1 材料與方法

1.1 材料

從PlantTFDB數(shù)據(jù)庫[43]（http：//planttfdb.cbi.pku.edu.cn/）中下載蔓花生、野花生、花生3個落花生屬和擬南芥的AP2蛋白序列，以及二穗短柄草（Brachypodium distachyon）AP2家族成員Bradi2g37800.13.p、鷹嘴豆AP2家族成員XP 004506004.1、黃花蒿（Artemisia annua）AP2家族成員Aan013559、木豆（Cajanus cajan）AP2家族成員C.cajan 21449、大豆（Glycine max）AP2家族成員Glyma.12G073300.1.p、綠豆（Vigna radiata）AP2家族成員Vradi03g03260.1、蕓豆（Phaseolus vulgaris）AP2家族成員Phvul.011G071100.1、向日葵（Helianthus annuus）AP2家族成員Han007343、芝麻（Sesamum indicum）AP2家族成員XP 011089852、蓖麻（Ricinus communis）AP2家族成員28752.m000339、玉米（Zea mays）AP2家族成員GRMZM2G700665 P03、印度水稻（Oryza sativa ssp. indica）AP2家族成員BGIOSGA018911-PA、小麥（Triticum urartu）AP2家族成員TRIUR3 14819-P1的蛋白序列。

1.2 方法

1.2.1 蔓花生AP2基因家族理化性質(zhì)的分析及亞細(xì)胞定位 利用在線程序PrtoParam[44-47]（https：//web.expasy.org/protparam/）統(tǒng)計蔓花生AP2基因家族的編碼氨基酸數(shù)量、理論等電點、分子量、帶正電荷氨基酸殘基數(shù)、帶負(fù)電荷氨基酸殘基數(shù)和脂肪族氨基酸指數(shù)及總平均親疏水性等理化性質(zhì)信息，并在PlantTFDB數(shù)據(jù)庫中以擬南芥基因序列為探針，進(jìn)行BLAST對比。利用在線軟件CELLOv2.5（http：//cello.life.nctu.edu.tw/）[48]對17條蔓花生的AP2蛋白序列進(jìn)行亞細(xì)胞定位。

1.2.2 氨基酸保守結(jié)構(gòu)域的預(yù)測分析 利用NCBI Consered Domains[49]（https：//www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi）對蔓花生AP2基因家族編碼蛋白進(jìn)行氨基酸保守結(jié)構(gòu)域的預(yù)測，對該家族各個成員的保守結(jié)構(gòu)域數(shù)量和位置進(jìn)行統(tǒng)計。

1.2.3 氨基酸信號肽及跨膜結(jié)構(gòu)的預(yù)測分析 使用SignalP 3.0 Server[44，46]（http：//www.cbs.dtu.dk/services/SignalP-3.0/）預(yù)測蔓花生AP2基因家族的17個成員的信號肽位置;在TMHMM Severv.2.0[44，46-47，49]（http：//www.cbs.dtu.dk/services/TMHMM/）在線程序中輸入AP2蛋白序列，分析該蛋白是否存在跨膜現(xiàn)象以及是否含有跨膜結(jié)構(gòu)域。

1.2.4 氨基酸親疏水性的預(yù)測分析 利用在線工具Protscale[44]（https：//web.expasy.org/protscale/）對蔓花生AP2蛋白進(jìn)行親疏水性預(yù)測分析，對其各個成員的最大親疏水性所在位點和分值進(jìn)行統(tǒng)計。

1.2.5 磷酸化位點的預(yù)測分析 使用在線工具NetPhos 3.1 Serve[44，47]（http：//www.cbs.dtu.dk/services/NetPhos/）對AP2蛋白序列的潛在磷酸化位點進(jìn)行預(yù)測分析，找出各個成員最可能的磷酸化位點及其磷酸化評分值，以及可能的磷酸化位點數(shù)量與占比。

1.2.6 蛋白質(zhì)二級結(jié)構(gòu)和三級結(jié)構(gòu)的預(yù)測分析 使用SOPMA[44]（https：//npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=npsa_sopma.html）程序?qū)βㄉ鶤P2蛋白的二級結(jié)構(gòu)進(jìn)行預(yù)測分析，對該家族各個成員的二級結(jié)構(gòu)類型在各成員中的占比和分布進(jìn)行統(tǒng)計;借助Phyre2網(wǎng)站[44，46-47，50]（http：//www.sbg.bio.ic.ac.uk/phyre2）對AP2蛋白進(jìn)行三級結(jié)構(gòu)的預(yù)測。

1.2.7 蔓花生AP2基因家族的蛋白基序分析 利用MEME在線工具[50-51]對蔓花生AP2基因家族成員進(jìn)行motif查找，設(shè)置查找的motif數(shù)量為10個，其余參數(shù)設(shè)為默認(rèn)值。

1.2.8 系統(tǒng)進(jìn)化樹的構(gòu)建 用MEGA軟件[44-46，49-50]分別構(gòu)建蔓花生AP2基因家族的系統(tǒng)發(fā)育進(jìn)化樹，落花生屬蔓花生（Arachis duranensis）、野花生（Arachis ipaensis）和花生（Arachis hypogaea L.）的系統(tǒng)發(fā)育進(jìn)化樹與蔓花生AP2基因家族、擬南芥AP2基因家族的系統(tǒng)發(fā)育進(jìn)化樹，以及蔓花生AP2家族成員Aradu.387PF、二穗短柄草AP2家族成員Bradi2g37800.13.p、擬南芥AP2家族成員AT2G28550.3、鷹嘴豆AP2家族成員XP 004506004.1、黃花蒿AP2家族成員Aan013559、木豆AP2家族成員C.cajan 21449、大豆AP2家族成員Glyma.12G073300.1.p、綠豆AP2家族成員Vradi03g03260.1、蕓豆AP2家族成員Phvul.011G071100.1、向日葵AP2家族成員Han007343、芝麻AP2家族成員XP 011089852、蓖麻AP2家族成員28752.m000339、玉米AP2家族成員GRMZM2G700665 P03、印度水稻AP2家族成員BGIOSGA018911-PA、小麥AP2家族成員TRIUR3_14819-P1的系統(tǒng)發(fā)育進(jìn)化樹。

2 結(jié)果與分析

2.1 蔓花生AP2基因家族編碼蛋白理化性質(zhì)的分析及亞細(xì)胞定位

對蔓花生AP2基因家族成員編碼蛋白的理化性質(zhì)進(jìn)行分析發(fā)現(xiàn)，該家族各蛋白的氨基酸數(shù)量范圍為237～752個，氨基酸數(shù)量最多的是 Aradu.M2I8I 家族成員（752個），氨基酸數(shù)量最少的是Aradu.A4KG4家族成員（237個）;蛋白質(zhì)的理論等電點除Aradu.T7UUH成員尚存在爭議外，其余16個成員的等電點都能分析得出，其中，除 Aradu.387PF、Aradu.RG31X、Aradu.0KI78和Aradu.GIP2Q成員的等電點大于8，Aradu.EJI6Q成員的等電點為7.70外，其余11個成員的等電點均小于7，說明該家族含有大量酸性氨基酸;蛋白不穩(wěn)定指數(shù)分析結(jié)果顯示，除Aradu.M2I8I、Aradu.A4KG4 這2個成員的不穩(wěn)定指數(shù)小于40外，其余15個成員的不穩(wěn)定指數(shù)均大于40，說明該家族成員多為不穩(wěn)定蛋白;該家族的脂肪族氨基酸指數(shù)分布范圍為4961（Aradu.RG31X）～7258（Aradu.S9CWT），脂肪族氨基酸指數(shù)反映了蛋白質(zhì)的熱穩(wěn)定性[52]，由此說明該家族蛋白質(zhì)間的熱穩(wěn)定性差異較大。由表1中AP2基因的亞細(xì)胞定位結(jié)果看出，蔓花生AP2蛋白成員均分布在細(xì)胞核中。在蔓花生AP2基因家族成員中，可按對應(yīng)氨基酸組成分為Ⅲ類，其中第Ⅰ類以絲氨酸占比最高，第Ⅱ類的絲氨酸和天冬酰胺占比一樣，第Ⅲ類以天冬酰胺占比最高（圖1）。

2.2 氨基酸保守結(jié)構(gòu)域分析

對蔓花生AP2基因家族編碼蛋白的氨基酸保守結(jié)構(gòu)域進(jìn)行預(yù)測分析，結(jié)果顯示，Aradu.387PF、Aradu.42K79、Aradu.EJI6Q、Aradu.G1CDJ、Aradu.M2I8I、Aradu.RG31X、Aradu.GIP2Q、Aradu.SE6Q0這8個成員均存在2個保守結(jié)構(gòu)域，其余9個成員均含有1個保守結(jié)構(gòu)域。由表2、圖2可以看出，蔓花生AP2基因家族編碼蛋白成員保守結(jié)構(gòu)域所含氨基酸數(shù)量大于100個的有7個;所含氨基酸數(shù)量最多的保守結(jié)構(gòu)域在Aradu.M2I8I成員上，共有135個氨基酸;所含氨基酸數(shù)量最少的保守結(jié)構(gòu)域在Aradu.A4KG4成員上，共有59個氨基酸。

2.3 氨基酸信號肽及跨膜結(jié)構(gòu)分析

對蔓花生AP2基因家族編碼蛋白各氨基酸信號肽進(jìn)行預(yù)測，分析發(fā)現(xiàn)，蔓花生AP2基因家族編碼蛋白的17個成員均無信號肽，因此推測蔓花生AP2基因家族編碼蛋白成員均為非分泌蛋白。

對蔓花生AP2基因家族編碼蛋白成員進(jìn)行跨膜結(jié)構(gòu)域預(yù)測分析發(fā)現(xiàn)，只有Aradu.SE6Q0和Aradu.42K79這2個成員存在跨膜現(xiàn)象。成員Aradu.SE6Q0在第236個氨基酸與第258個氨基酸之間存在跨膜現(xiàn)象，但并不形成跨膜結(jié)構(gòu)域;成員Aradu.42K79在第349個氨基酸與第370個氨基酸之間存在跨膜現(xiàn)象，也并不存在跨膜結(jié)構(gòu)域。經(jīng)SignalP3.0 Server程序分析可知，蔓花生AP2基因家族編碼蛋白成員均不存在跨膜結(jié)構(gòu)區(qū)域，因此推斷該家族蛋白為非跨膜蛋白。

2.4 氨基酸親疏水性預(yù)測分析

對蔓花生AP2基因家族編碼蛋白的17個成員進(jìn)行親疏水性預(yù)測分析發(fā)現(xiàn)，該家族成員Aradu.387PF的第17位精氨酸的親水性最強，為-3.756;成員Aradu.42K79的第262位甘氨酸的疏水性最強，為3.789;總平均親疏水性均小于0，說明AP2基因家族編碼蛋白成員均為親水性蛋白。最高親水性值為-1.004，在該基因家族編碼蛋白成員Aradu.A4KG4上;最低親水性值為-0.449，在成員Aradu.T7UUH上（圖3）。

2.5 磷酸化位點分析

對蔓花生AP2基因家族編碼蛋白的潛在磷酸化位點進(jìn)行預(yù)測分析發(fā)現(xiàn)，在該家族同一成員之間， 絲氨酸含量總是最高的，其次是蘇氨酸，最低的是酪氨酸;絲氨酸最可能的磷酸化位點數(shù)最多，且磷酸化評分值往往在0.997～0.998之間，只有Aradu.A4KG4、Aradu.S9CWT的磷酸化評分值低于0.997，分別為0.925、0.994;蘇氨酸、酪氨酸最可能的磷酸化位點數(shù)都只有1個，且磷酸化評分值都低于0.990;在蔓花生AP2蛋白家族成員中，Aradu.387PF具有的潛在磷酸化位點數(shù)最多，磷酸化位點數(shù)占該家族氨基酸總數(shù)的14.99%;成員Aradu.A4KG4含有的潛在磷酸化位點數(shù)最少，磷酸化位點數(shù)占該家族氨基酸總數(shù)的8.02%;在AP2蛋白家族成員中，成員Aradu.SE6Q0的絲氨酸、蘇氨酸含量都是較高的，占比分別為8.65%、5.19%;成員Aradu.A4KG4的絲氨酸、蘇氨酸含量都是較低的，占比分別為3.38%、2.53%，酪氨酸含量最高，占比為2.11%;成員Aradu.42K79的酪氨酸含量最低，占比為0.19%（圖4）。

2.6 蛋白質(zhì)二級結(jié)構(gòu)和三級結(jié)構(gòu)的預(yù)測分析

對蔓花生AP2蛋白家族成員進(jìn)行二級結(jié)構(gòu)分析，結(jié)果顯示，17個蔓花生AP2基因家族編碼蛋白成員的二級結(jié)構(gòu)均表現(xiàn)為無規(guī)則卷曲最多，β-轉(zhuǎn)角最少，有14個成員的α-螺旋數(shù)>延伸鏈數(shù)，Aradu.GG3LG、Aradu.QYV19和Aradu.T7UUH這3個成員的延伸鏈數(shù)>α-螺旋數(shù)?？傮w而言，蔓花生AP2基因家族編碼蛋白的二級結(jié)構(gòu)以無規(guī)則卷曲和α-螺旋為主，而延伸鏈和β-轉(zhuǎn)角則分散在各個氨基酸鏈中（圖5）。

對蔓花生AP2基因家族編碼蛋白的17個成員進(jìn)行三級結(jié)構(gòu)預(yù)測并對預(yù)測結(jié)果進(jìn)行比對，根據(jù)各成員三級結(jié)構(gòu)是否相似，可將該家族分為Ⅶ大類，分別用Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ、Ⅶ類表示。由圖6可以看出，第Ⅳ類所含成員數(shù)最多，一共有5個，分別是Aradu.RG31X、Aradu.M2I8I、Aradu.SE6Q0、Aradu.G1CDJ和Aradu.42K79;其次是第Ⅰ類，含有Aradu.0KI78、Aradu.387PF和Aradu.VJK1R這3個成員;第Ⅱ、Ⅲ、Ⅴ和Ⅵ類均含有2個成員，Aradu.GIP2Q、Aradu.QYV19屬于第Ⅱ類成員，第Ⅲ類成員有Aradu.A4KG4和Aradu.GG3LG，第Ⅴ類成員有Aradu.C5HAC和Aradu.UA79E，第Ⅵ類成員有Aradu.T7UUH和Aradu.S9CWT;第Ⅶ類只有1個成員，為Aradu.GIP2Q。

2.7 蔓花生AP2基因家族編碼蛋白基序分析

對蔓花生AP2基因家族編碼蛋白成員進(jìn)行Motif分析發(fā)現(xiàn)，Motif3、Motif4都含有1個共同的結(jié)構(gòu)域，即YRG;Motif5含有RAYD結(jié)構(gòu)域（圖7）。YRG結(jié)構(gòu)域有利于DNA的結(jié)合，RAYD結(jié)構(gòu)域主要調(diào)節(jié)轉(zhuǎn)錄因子與DRE順式元件結(jié)合[53]，說明蔓花生AP2基因家族編碼蛋白的Motif3、Motif4和Motif5屬于AP2結(jié)構(gòu)域。根據(jù)蔓花生AP2蛋白序列所含Motif類型的不同，將該家族分為Ⅴ大類，第Ⅰ類主要含有Motif1、Motif4和Motif8，其成員有Aradu.GG3LG、Aradu.UA79E、Aradu.S9CWT、Aradu.C5HAC、Aradu.T7UUH 和Aradu.VJK1R;第Ⅱ類包括成員Aradu.M2I8I、Aradu.QYV19，都含有Motif1、Motif3、Motif4、Motif6、Motif7、Motif8;第Ⅲ類主要含有Motif1、Motif2、Motif3和Motif5，成員有Aradu.EJI6Q、Aradu.387PF和Aradu.G1CDJ;第Ⅳ類主要含有Motif1和Motif2，成員有Aradu.A4KG4和Aradu.SE6Q0;第Ⅴ類主要含有Motif2、Motif3、Motif4、Motif5、Motif6，成員有Aradu.42K79和Aradu.GIP2Q。在蔓花生AP2蛋白序列中，含有2個AP2保守結(jié)構(gòu)域的蛋白質(zhì)序列均含有Motif5，在蔓花生AP2基因家族編碼蛋白的17個成員中，除Aradu.GIP2Q、Aradu.A4KG4不含有Motif1外，其他15個成員均含有Motif1;Motif2在該家族的15個成員中均存在，只有Aradu.QYV19、Aradu.GG3LG中不存在;有9個成員含有Motif3，14個成員含有Motif4，8個成員含有Motif5，說明Motif1、Motif2和Motif4是蔓花生AP2基因家族編碼蛋白的主要保守基序（圖8）。

2.8 系統(tǒng)進(jìn)化關(guān)系分析

通過對蔓花生AP2基因家族編碼蛋白系統(tǒng)發(fā)育進(jìn)化樹的分組比對發(fā)現(xiàn)，該家族成員Aradu.GG3LG、Aradu.UA79E、Aradu.S9CWT、Aradu.C5HAC、Aradu.T7UUH和Aradu.VJK1R的親緣關(guān)系比該家族其他成員的親緣關(guān)系近。其中Aradu.GG3LG和Aradu.UA79E、Aradu.387PF和Aradu.G1CDJ、Aradu.A4KG4和Aradu.SE6Q0、Aradu.42K79和Aradu.GIP2Q這4組的組內(nèi)成員親緣關(guān)系更近（圖8），這與以擬南芥為探針的BLAST比對結(jié)果基本吻合（表1）。

落花生屬3個物種的系統(tǒng)進(jìn)化關(guān)系可以分為Ⅶ大類，第Ⅰ類有9個成員，主要是蔓花生AP2基因家族編碼蛋白的5個成員和野花生的4個成員;第Ⅱ類有7個成員，主要是野花生AP2基因家族編碼蛋白的3個成員，蔓花生AP2基因家族編碼蛋白的3個成員及花生AP2基因家族編碼蛋白的成員Ahy016292，其中Araip.XB73L與Ahy016292的親緣關(guān)系相對較近;第Ⅲ類有7個成員，野花生有3個， 蔓花生有2個，花生有2個，其中Araip.QQN2T與Ahy022503的親緣關(guān)系相對較近，Aradu.A4KG4與Araip.80H4B的親緣關(guān)系相對較近，Aradu.0KI78與Ahy008861的親緣關(guān)系相對較近;第Ⅳ類有4個成員，分別是蔓花生的Aradu.42K79、Aradu.GIP2Q，野花生的Araip.L8K4F，花生的Ahy022725，其中Araip.L8K4F與Ahy022725的親緣關(guān)系相對較近;第Ⅴ類有7個成員，分別是蔓花生的Aradu.M2I8I、Aradu.SE6Q0、Aradu.C5HAC，野花生的Araip.AXK3N、Araip.RB7C2，花生的Ahy015793、Ahy009871;第Ⅵ類有5個成員，分別是蔓花生的Aradu.EJI6Q、Aradu.QYV19，野花生的Araip.SZ63C，花生的Ahy013042、Ahy017508;第Ⅶ類與其他類別的親緣關(guān)系相對較遠(yuǎn)，有2個成員都是野花生AP2蛋白家族成員，分別是Araip.N6C0B、Araip.E29R3（圖9-a）。

蔓花生AP2基因家族編碼蛋白和擬南芥AP2基因家族編碼蛋白的系統(tǒng)進(jìn)化關(guān)系可分為Ⅴ類，第Ⅰ類有11個成員的親緣關(guān)系相對較近，其中有7個是蔓花生的成員，4個是擬南芥的成員;第Ⅱ類全是蔓花生的4個成員; 第Ⅲ類只有1個成員，是蔓花生的Aradu.SE6Q0;第Ⅳ類有14個成員，其中蔓花生有3個，擬南芥有11個;第Ⅴ類有17個成員，其中蔓花生有9個，擬南芥有8個，成員Aradu.77UUH與成員AT5G57390.1的親緣關(guān)系相對較近（圖9-b）。

對蔓花生與二穗短柄草等15個物種AP2蛋白的系統(tǒng)進(jìn)化關(guān)系分析發(fā)現(xiàn)，蔓花生的Aradu.42K79與大豆的Glyma.12G073300.1.p親緣關(guān)系最近，與木豆的C.cajan.21449、蕓豆的Phvul.011G071100.1親緣關(guān)系較近，這與其同屬于豆科有關(guān)，但同為豆科的綠豆Vradi03g03260.1與其親緣關(guān)系較遠(yuǎn)，說明綠豆該蛋白的進(jìn)化速度更快;二穗短柄草的Bradi2g37800.13.p與蔓花生的Aradu.42K79親緣關(guān)系最遠(yuǎn)，可能與其分別為單雙子葉植物有關(guān)（圖10）。

3 討論與結(jié)論

目前，關(guān)于AP2轉(zhuǎn)錄因子的相關(guān)研究顯示，AP2亞家族主要參與植物的生長發(fā)育過程，對花發(fā)育的調(diào)控研究得最多[37，54-58]，對根的生長、果實和種子發(fā)育的調(diào)控作用也有研究[25，59]。吳彥慶等通過研究芍藥（Paeonia lactiflora Pall.）分生組織決定基因AP[STBX]2[STBZ]發(fā)現(xiàn)，plAP2蛋白為親水性不穩(wěn)定蛋白，無跨膜結(jié)構(gòu)和信號肽，非分泌蛋白[60]。阮先樂等對甘薯[Dioscorea esculenta （Lour.） Burkill]AP2蛋白的研究發(fā)現(xiàn)，AP2蛋白為不穩(wěn)定的、無跨膜結(jié)構(gòu)域的親水性非分泌蛋白[61]。鄧?yán)说韧ㄟ^對‘紅陽獼猴桃（Actinidia chinensis）AP2/EREB轉(zhuǎn)錄因子的研究發(fā)現(xiàn)，該家族成員多為酸性氨基酸，是無信號肽的非分泌親水性蛋白[62]，這與本研究中蔓花生AP2蛋白的性質(zhì)基本吻合。芍藥的AP2蛋白主要在胞質(zhì)中起著生物學(xué)功能[60]，甘薯的大多數(shù)AP2成員亞細(xì)胞定位在細(xì)胞核[61]。甘藍(lán)（Brassica oleracea L.）AP2蛋白的亞細(xì)胞定位結(jié)果顯示，其AP2蛋白多出現(xiàn)于細(xì)胞核中[63];喬永剛等通過對金銀花（Lonicera japonica Thunb.）AP2蛋白的亞細(xì)胞定位分析發(fā)現(xiàn)，其蛋白多在細(xì)胞核中[64];此外，新疆紫草[Arnebia euchroma （Royle） Johnst]的AP2蛋白也定位于細(xì)胞核中[65]。而本研究結(jié)果顯示，蔓花生AP2蛋白均定位在細(xì)胞核中，說明AP2蛋白在不同植物中具有的作用和承擔(dān)的功能可能并不相同，但其最可能在細(xì)胞核中起著生物學(xué)作用。蓖麻的16個AP2成員以α-螺旋、無規(guī)則卷曲為其主要的二級結(jié)構(gòu)[66]。郭慧等通過對甘藍(lán)AP2/ERF轉(zhuǎn)錄因子的研究發(fā)現(xiàn)，無規(guī)則卷曲、α-螺旋和延伸鏈?zhǔn)歉仕{(lán)AP2蛋白的主要二級結(jié)構(gòu)[63]。文冠果（Xanthoceras sorbifolium Bunge）AP2蛋白的二級結(jié)構(gòu)主要由無規(guī)則卷曲、α-螺旋和延伸鏈組成[67]。楊樹（Populus trichocarpa）AP2蛋白以無規(guī)則卷曲為主[68]。在本研究中，蔓花生AP2蛋白也是以無規(guī)則卷曲、α-螺旋和延伸鏈為主，說明這3種二級結(jié)構(gòu)是AP2蛋白的主要二級結(jié)構(gòu)。

本研究結(jié)果表明，蔓花生AP2蛋白為不穩(wěn)定的親水性非分泌蛋白，亞細(xì)胞定位在細(xì)胞核中，二級、三級結(jié)構(gòu)以無規(guī)則卷曲、α-螺旋和延伸鏈為主，含有YRG、RAYD結(jié)構(gòu)域，花生屬的Araip.XB73L和Ahy016292、Araip.QQN2T和Ahy022503、Aradu.A4KG4和Araip.80H4B、Aradu.0KI78和Ahy008861、Araip.L8K4F和Ahy022725這5組組內(nèi)成員的親緣關(guān)系相對較近，蔓花生AP2成員Aradu.77UUH與擬南芥AP2成員AT5G57390.1的親緣關(guān)系最為接近。物種進(jìn)化關(guān)系分析結(jié)果顯示，蔓花生的Aradu.42K79與大豆的Glyma.12G073300.1.p親緣關(guān)系最近，與二穗短柄草的Bradi2g37800.13.p親緣關(guān)系最遠(yuǎn)。這些研究結(jié)果為進(jìn)一步研究和揭示AP2基因家族在蔓花生生長發(fā)育、基因調(diào)控中的作用提供了一定的參考依據(jù)。

參考文獻(xiàn)：

[1]周 立.蔓花生（Arachis duranensis Krap.et Greg）在低溫脅迫下的生理反應(yīng)研究[D]. 重慶：西南大學(xué)，2008.

[2]張 利. Pb2+、Hg2+脅迫下蔓花生逆境生理及田間不定根發(fā)生研究[D]. 重慶：西南大學(xué)，2008.

[3]張振霞，鄭玉忠. 兩種蔓花生坪用性狀的比較[J]. 湖北農(nóng)業(yè)科學(xué)，2008，47（9）：1053-1055.

[4]鄭 荷. 地被植物蔓花生再生體系建立及多倍體誘導(dǎo)[D]. 重慶：西南大學(xué)，2008.

[5]蔚冬紅，喬善義，趙毅民.中藥合歡皮研究概況[J]. 中國中藥雜志，2004，29（7）：619-621.

[6]唐湘梧. 新的牧草栽培品種[J]. 草原與牧草，1991（3）：41-42.

[7]羅 丹，陳紅躍，劉 乾. 5種屋頂綠化植物抗旱性研究[J]. 廣東林業(yè)科技，2009，25（6）：81-85.

[8]應(yīng)朝陽，呂亮雪，劉國道，等. 干旱脅迫對多年生落花生生長和活性氧代謝的影響[J]. 熱帶作物學(xué)報，2006，27（1）：17-21.

[9]黃毅斌，應(yīng)朝陽，柯碧南，等. 兩種多年生花生品種低溫脅迫下的質(zhì)膜透性研究[J]. 四川草原，2005（10）：4-5，12.

[10]余高鏡，應(yīng)朝陽，黃毅斌，等. 2個多年生花生品種低溫脅迫下活性氧代謝特性研究[J]. 草業(yè)科學(xué)，2006（4）：29-32.

[11]周 立，鄭 荷，張 利，等. 低溫脅迫下磷酸二氫鉀對蔓花生的生理影響研究[J]. 現(xiàn)代農(nóng)業(yè)科學(xué)，2008（2）：27-28，31.

[12]彭建宗，陳兆平. 遮陰對多年生花生（Arachis pintoi）生長的影響[J]. 華南師范大學(xué)學(xué)報（自然科學(xué)版），1999（2）：92-94.

[13]彭建宗. 優(yōu)良地被植物多年生花生[J]. 廣東園林，2002（4）：41-42.

[14]陳偉俊. 一種優(yōu)良的園林地被植物——金花生[J]. 熱帶農(nóng)業(yè)科技，2006，29（1）：39.

[15]譚光營，潘志權(quán). 優(yōu)良植被——多年生花生的栽培及其應(yīng)用[J]. 廣西熱帶農(nóng)業(yè)，2004（1）：46-47.

[16]張 妍，劉 瀛，孫豐賓，等. 白樺APETALA2（AP2）轉(zhuǎn)錄因子基因的分離及其表達(dá)[J]. 林業(yè)科學(xué)研究，2012，25（2）：254-260.

[17]高春艷，吳 芮，袁 玉，等. 植物AP2/ERF轉(zhuǎn)錄因子及其在非生物脅迫應(yīng)答中的作用[J]. 江漢大學(xué)學(xué)報（自然科學(xué)版），2017，45（3）：236-240.

[18]張計育，王慶菊，郭忠仁. 植物AP2/ERF類轉(zhuǎn)錄因子研究進(jìn)展[J]. 遺傳，2012，34（7）：835-847.

[19]季愛加，羅紅梅，徐志超，等. 藥用植物轉(zhuǎn)錄因子AP2/ERF研究與展望[J]. 科學(xué)通報，2015，60（14）：1272-1284.

[20]Jofuku K D，Boer B，Montagu M V，et al. Control of Arabidopsis flower and seed development by the homeotic gene APETALA[STBX]2[STBZ][J]. Plant Cell，1994（6）：1211-1225.

[21]Toshitsugu N，Kaoru S，Tatsuhito F，et al. Genome-wide analysis of the ERF gene family in Arabidopsis and rice[J]. Plant Physiol，2006，140（2）：411-432.

[22]聞可心，劉雪梅. AP[STBX]2[STBZ]功能基因在植物花發(fā)育中的重要作用[J]. 生物技術(shù)通報，2010（2）：1-7.

[23]Houston K，McKim S M，Comadran J，et al. Variation in the interaction between alleles of HvAPETALA2 and microRNA172 determines the density of grains on the barley inflorescence[J]. Proceedings of the National Academy of Sciences of the United States of America，2013，110（41）：16675-16680.

[24]劉 強，張貴友，陳受宜. 植物轉(zhuǎn)錄因子的結(jié)構(gòu)與調(diào)控作用[J]. 科學(xué)通報，2000，45（14）：1465-1474.

[25]Yan X H，Zhang L，Chen B，et al. Functional identification and characterization of the Brassica napus transcription factor gene BnAP[STBX]2[STBZ]，the ortholog of Arabidopsis thaliana APETALA[STBX]2[STBZ][J]. PLoS One，2012，7（3）：e33890.

[26]Shigyo M，Ito M. Analysis of gymnosperm two-AP2-domain-containing genes[J]. Dev Genes Evol，2004，214（3）：105-114.

[27]Kim S，Soltis P S，Wall K，et al. Phylogeny and domain evolution in the APETALA[STBX]2[STBZ]-like gene family[J]. Mol Biol Evol，2006，23（1）：107-120.

[28]Nole-Wilson S，Krizek B A. DNA binding properties of the Arabidopsis floral development protein AINTEGUMENTA[J]. Nucleic Acids Res，2000，28：4076-4082.

[29]Ohme-Takagi M，Shinshi H. Ethylene-inducible DNA binding proteins that interact with anethylene-responsive element[J]. Plant Cell，1995，7（2）：173-182.

[30]Fits L V D，Memelink J. The jasmonate-inducible AP2/ERF-domain transcription factor ORCA3 activates gene expression via interaction with a jasmonate-responsive promoter element[J]. The Plant Journal，2001，25（1）：43-53.

[31]Gu Y Q，Wildermuth M C，Chakravarthy S，et al. Tomato transcription factors pti4，pti5，and pti6 activate defense responses when expressed in Arabidopsis[J]. Plant Cell，2002，14（4）：817-831.

[32]Sasaki K，Mitsuhara I，Seo S，et al. Two novel AP2/ERF domain proteins interact with cis-element VWRE for wound-induced expression of the Tobacco tpoxN[STBX]1[STBZ] gene[J]. The Plant Journal，2007，50（6）：1079-1092.

[33]宋佳靜. 花生突變體庫的構(gòu)建及APETALA[STBX]2[STBZ]基因的表達(dá)分析[D]. 鄭州：河南農(nóng)業(yè)大學(xué)，2013.

[34]Haughn G W，Somerville C R. Genetic control of morphogenesis in Arabidopsis[J]. Dev Genet，1988，9（2）：73-89.

[35]Bowman J L，Smyth D R，Meyerowitz E M. Genes directing flower development in Arabidopsis[J]. Plant Cell，1989，1（1）：37-52.

[36]Bowman J L，Drews G N，Meyerowitz E M. Expression of the Arabidopsis floral homeotic gene agamous is restricted to specific cell types late in flower development[J]. Plant Cell，1991，3（8）：749-758.

[37]Kunst L，Klenz J E，Martinezzapater J，et al. Ap[STBX]2[STBZ] gene determines the identity of perianth organs in flowers of Arabidopsis thaliana[J]. Plant Cell，1989，1（12）：1195-1208.

[38]Huala E，Sussex I M. Leafy interacts with floral homeotic genes to regulate Arabidopsis floral development[J]. Plant Cell，1992，4（8）：901-913.

[39]Ohto M A，F(xiàn)ischer R L，Goldberg R B，et al. Control of seed mass by APETALA[STBX]2[STBZ][J]. Proceedings of the National Academy of Sciences，2005，102（8）：3123-3128.

[40]Morcillo F，Gallard A，PilloT M，et al. EgAP2-1，an AINTEGUMENTA-like （AIL） gene expressed in meristematic and proliferating tissues of embryos in oil palm[J]. Planta，2007，226：1353-1362.

[41]Zhuang J，Anyia A，Vidmar J，et al. Discovery and expression assessment of the AP2-like genes in Hordeum vulgare[J]. Acta Physiol Plant，2011，33（5）：1639-1649.

[42]Jofuku K D，Omidyar P K，Gee Z，et al. Control of seed mass and seed yield by the floral homeotic gene APETALA[STBX]2[STBZ][J]. Proceedings of the National Academy of Sciences，2005，102（8）：3117-3122.

[43]李 寧，高升華，王 飛，等. 辣椒CONSTANS-like轉(zhuǎn)錄因子家族的生物信息學(xué)分析[J]. 北方園藝，2019（15）：1-7.

[44]王 雙，李 躍，李孟建，等. 蓖麻APs基因克隆及生物信息學(xué)分析[J]. 分子植物育種，2019，17（22）：7344-7349.

[45]閆 艷，王 希，周珂輝，等. 西洋參bHLH轉(zhuǎn)錄因子家族生物信息學(xué)分析[J]. 吉林農(nóng)業(yè)大學(xué)學(xué)報，2019，41（3）：316-323.

[46]李孟建，李國瑞，黃鳳蘭，等. 蓖麻RcGASA[STBX]9[STBZ]基因克隆及生物信息學(xué)分析[J]. 分子植物育種，2019，17（21）：6981-6986.

[47]狄建軍，孫佳欣，楊智慧，等. 蓖麻PLA[STBX]2[STBZ]α基因啟動子的克隆及生物信息學(xué)分析[J]. 分子植物育種，2019，17（15）：4961-4966.

[48]Yu C S，Chen Y C，Lu C H，et al. Prediction of protein subcellular localization[J]. Proteins：Structure，F(xiàn)unction，and Bioinformatics，2006，64（3）：643-651.

[49]劉承圓，何其光，林春花，等. 巴西橡膠樹（Hevea brasiliensis） HbMlo2基因克隆與表達(dá)分析[J]. 分子植物育種，2019，17（5）：1512-1518.

[50]李惠敏，李璐璐，陳玉梅，等. 羅漢果幾丁質(zhì)酶基因家族的生物信息學(xué)分析[J]. 生物信息學(xué)，2019（3）：167-174.

[51]陳湘瑜，徐日榮，陳 昊，等. 杜蘭落花生球蛋白基因家族的生物信息學(xué)分析[J]. 中國油料作物學(xué)報，2019，41（3）：391.

[52]徐志軍，劉 洋，徐 磊，等. 玉米轉(zhuǎn)錄因子NF-YB基因家族的生物信息學(xué)分析[J]. 分子植物育種，2019，17（12）：3807-3816.

[53]周 凡，李 丹，周 軍，等. 梨DREB基因家族全基因組鑒定及分析[J]. 基因組學(xué)與應(yīng)用生物學(xué)，2018，37（9）：3889-3902.

[54]Drews G N，Bowman J L，Meyerowitz E M. Negative regulation of the Arabidopsis homeotic gene AGAMOUS by the APETALA2 product[J]. Cell，1991，65（6）：991-1002.

[55]Shannon S，Meeks-Wagner D R. Genetic interactions that regulate inflorescence development in Arabidopsis[J]. Plant Cell，1993，5（6）：639-655.

[56]Houston K，McKim S M，Comadran J，et al. Variation in the interaction between alleles of HvAPETALA[STBX]2[STBZ] and microRNA172 determines the density of grains on the barley inflorescence[J]. Proc Natl Acad Sci USA，2013，110（41）：16675-16680.

[57]Jofuku K D，den Boer B G，van Montagu M，et al. Control of Arabidopsis flower and seed development by the homeotic gene APETALA2[J]. Plant Cell，1994，6：1211-1225.

[58]Zhang P，Yang P Z，Zhang Z Q，et al. Isolation and characterization of a buffalograss （Buchloe dactyloides） dehydration responsive element binding transcription factor，BdDREB2[J]. Gene，2014，536：123-128.

[59]Aya K，Hobo T，Sato-Izawa K，et al. A novel AP2-type transcription factor，SMALL ORGAN SIZE1，controls organ size downstream of an auxin signaling pathway[J]. Plant Cell Physiol，2014，55（5）：897-912.

[60]吳彥慶，成夢琳，趙大球，等. 芍藥分生組織決定基因APETALA[STBX]2（AP2[STBZ]）的克隆及生物信息學(xué)分析[J]. 華北農(nóng)學(xué)報，2017，32（3）：58-64.

[61]阮先樂，張 杰，張福麗，等. 甘薯AP2基因家族的生物信息學(xué)分析[J]. 分子植物育種，2017，15（6）：2066-2072.

[62]鄧 浪，沈兵琪，王連春，等. ‘紅陽獼猴桃全基因組AP[STBX]2[STBZ]/EREBP轉(zhuǎn)錄因子生物信息學(xué)分析[J]. 果樹學(xué)報，2017，34（7）：790-805.

[63]郭 慧，金司陽，劉 寒，等. 甘藍(lán)AP2/ERF轉(zhuǎn)錄因子的克隆和生物信息學(xué)分析[J]. 中國藥師，2017，20（1）：6-10.

[64]喬永剛，陳 亮，崔芬芬，等. 基于轉(zhuǎn)錄組金銀花AP2基因家族的生物信息學(xué)及表達(dá)分析[J]. 核農(nóng)學(xué)報，2019，33（9）：1698-1706.

[65]謝 騰，王 升，周良云，等. 新疆紫草AP2/ERF轉(zhuǎn)錄因子的電子克隆和生物信息學(xué)分析[J]. 中國中藥雜志，2014，39（12）：2251-2257.

[66]閆星伊，風(fēng) 蘭，韓雯毓，等. 蓖麻AP2/ERF轉(zhuǎn)錄因子家族的生物信息學(xué)分析[J]. 分子植物育種，2020，18（4）：1088-1102.

[67]劉曉輝，馬靜怡，羅 媛，等. 木本油料文冠果APETALA[STBX]2[STBZ]基因全長cDNA序列與生物信息學(xué)分析[J]. 現(xiàn)代農(nóng)業(yè)科技，2018（2）：145-146.

[68]趙金玲，姚文靜，王升級，等. 楊樹AP2/ERF轉(zhuǎn)錄因子家族生物信息學(xué)分析[J]. 東北林業(yè)大學(xué)學(xué)報，2015，43（10）：21-29.

收稿日期：2019-09-15

基金項目：四川省科技廳項目（編號：17ZZ019）;樂山師范學(xué)院人才啟動項目（編號：XJR17005）。

作者簡介：王 玲（1999—），女，四川南充人，主要從事植物分子生物學(xué)研究。E-mail：1981887804@qq.com。

通信作者：付 春，博士，講師，主要從事經(jīng)濟植物的生物信息學(xué)研究。E-mail：fuchun421@aliyun.com。