楊丹 廉小平 周燕 張賀翠 杜丹 高啟國 任雪松 朱利泉

（西南大學(xué) 農(nóng)學(xué)部植物生理生化實驗室，重慶 400715）

結(jié)球甘藍ARC1與Exo70A1編碼區(qū)cDNA的克隆及進化分析

楊丹 廉小平 周燕 張賀翠 杜丹 高啟國 任雪松 朱利泉

（西南大學(xué) 農(nóng)學(xué)部植物生理生化實驗室，重慶 400715）

在甘藍自交不親和信號傳導(dǎo)過程中，ARC1與Exo70A1的相互作用起著承上啟下的關(guān)鍵作用。為了研究10種自交不親和結(jié)球甘藍材料ARC1和Exo70A1的進化關(guān)系，對其編碼區(qū)cDNA進行了克隆和測序；在此基礎(chǔ)上采用生物信息學(xué)方法，對這10種材料和前人已發(fā)表的83條ARC1和73條Exo70A1的相關(guān)序列的ARC1-Exo70A1的互作區(qū)和非互作區(qū)編碼序列進行了變異和進化分析。結(jié)果發(fā)現(xiàn)：（1）所選育的10種結(jié)球甘藍材料的ARC1與Exo70A1的編碼區(qū)核酸序列存在平行進化關(guān)系，且均聚類于蕓薹屬分支下；（2）ARC1比Exo70A1編碼區(qū)序列進化速率快；（3）ARC1和Exo70A1互作編碼區(qū)核酸序列存在明顯的協(xié)同進化，而ARC1和Exo70A1非互作編碼區(qū)核酸序列則在進化上呈現(xiàn)明顯分異；（4）ARC1和Exo70A1互作編碼區(qū)cDNA的進化均快于非互作編碼區(qū)。

ARC1；Exo70A1；自交不親和；進化關(guān)系；協(xié)同進化

蕓薹屬植物自交不親和反應(yīng)（self-incompatibility，SI）是沿著SCR（S cysteine-rich，S-富含半胱氨酸蛋白［1］）→SRK（S receptor kinase，S-受體激酶［2］） → ARC1（armadillo-repeat containing protein1）→Exo70A1→MOD［3］的途徑進行的。SRK作為柱頭S基因，由胞外S結(jié)構(gòu)域、跨膜結(jié)構(gòu)域和激酶結(jié)構(gòu)域的編碼區(qū)構(gòu)成，其同SCR識別的胞外S結(jié)構(gòu)域編碼區(qū)變異性較大［4-6］，而SRK蛋白完整的胞外域包括B-Lection結(jié)構(gòu)域，SLG結(jié)構(gòu)域，PAN-APPLE結(jié)構(gòu)域。在SI反應(yīng)途徑中，ARC1與Exo70A1的相互作用起著承上啟下的關(guān)鍵作用。研究表明，ARC1的C端含有一個由6-7個臂重復(fù)區(qū)組成的ARM功能域，N端含有一個UND功能域，中間有一個U-box功能域［7］。ARC1的C端結(jié)構(gòu)域是與SRK激酶域發(fā)生相互作用的區(qū)域［7］，而U-box保守結(jié)構(gòu)域是70多個氨基酸殘基構(gòu)成的UFD2（泛素融合蛋白質(zhì)），其主要作為泛素連接酶-E3起作用，在泛素化系統(tǒng)中特異性識別底物蛋白質(zhì)［8］。Exo70A1是SI信號傳導(dǎo)途徑中ARC1的下游底物，是Samuel等［9］運用酵母雙雜交文庫篩選技術(shù)從甘藍型油菜的柱頭細胞中分離出的一種新型SI信號傳導(dǎo)元件。Exo70A1蛋白的C端含有一個典型的Exo70結(jié)構(gòu)域［10］。油菜的Exo70A1基因含有12個外顯子及11個內(nèi)含子，研究表明Exo70Al可能是促進花粉水合、萌發(fā)及花粉管生長的親和因子，在SI反應(yīng)途徑中被ARC1參與的泛素/26S蛋白酶體途徑降解，從而導(dǎo)致柱頭對自花花粉的拒絕［11］。

Stone等的研究表明［7］，在ARC1和Exo70A1發(fā)生互作時，ARC1蛋白質(zhì)N端的U-box是不可缺少的結(jié)構(gòu)域；而楊佳等［11］的研究進一步表明羽衣甘藍ARC1與Exo70A1相互作用區(qū)段是ARC1的N端（UND區(qū)段）和Exo70A1的N端（此蛋白質(zhì)的1-270個氨基酸）。而目前根據(jù)S位點基因序列構(gòu)建基因樹做進化分析的文章也日益增多［12-14］。在這些研究的基礎(chǔ)上，本研究以高度自交不親和的結(jié)球甘藍為研究對象，并利用其編碼區(qū)序列分別對ARC1-Exo70A1互作區(qū)和非互作區(qū)編碼DNA進行進化分析，以期為選育的甘藍自交不親和材料在整個植物中初步建立起進化位置和進化關(guān)系。

1 材料與方法

1.1 材料

近年來選育的高度自交不親和結(jié)球甘藍（表1），由西南大學(xué)園藝園林學(xué)院蔬菜實驗室提供；大腸桿菌感受態(tài)DH5α、EasyPfu Taq DNA polymerase、膠回收試劑盒和TranScript First-strand cDNA Synthesis SuperMix、Mark III（ 目 錄 號：MD103） 購 自 北京全式金生物技術(shù)有限公司；T4 DNA 連接酶、pMDTM19-T Vector、Sac I與Sma I購自TaKaRa大連公司；RNA提取試劑盒RNAprep Pure Plant Kit（目錄號：DP432）購自天根生化科技（北京）有限公司；其他常規(guī)試劑均為國產(chǎn)分析純；引物合成與序列測定由北京華大基因公司完成。

表1 甘藍材料及其親和指數(shù)

1.2 方法

1.2.1 結(jié) 球 甘 藍 ARC1與 Exo70A1 CDS（Coding sequence，編碼區(qū)序列）的克隆及測序 運用Primer Premier 6.0軟件，依照Genbank中ARC1和Exo70A1已有的序列分別設(shè)計結(jié)球甘藍ARC1引物和Exo-70A1引物。所選擇的ARC1上游引物是（F-primer）5'-CCGGAATTCCCGATGGCCACTGATTCAGCAATGT TCGCAT-3'，下游引物是（L-primer）5'-ACGCGTCG ACGTTATCTCTGTGTTTTCTGGTCGCAC-3'；Exo70A1上游引物是（F-primer）5'-GGGAACATATGTACATG GCCGTCGATAGCGGAAT-3'，下游引物是（L-primer）5'-ACGCGTCGACCTTTACCGTCGTGGTTCATTCATAG ACT-3'。

收集開花前1 d結(jié)球甘藍花蕾的雌蕊，提取總RNA。以TranScript First-strand cDNA Synthesis SuperMix反轉(zhuǎn)錄合成進cDNA第一鏈；再選用引物對，以所得cDNA為模板，對應(yīng)的擴增ARC1和Exo70A1基因。擴增體系均為50 μL：模板4 μL，2 mmol/L dNTPs 5 μL，上下游引物各2 μL，10×Easypruf Buffer 5 μL，ddH2O 31 μL，TransTaq酶 1 μL。PCR擴增程序參數(shù)為：94℃預(yù)變性2 min；94℃預(yù)變性30 s，退火30 s，72℃延伸90 s，35個循環(huán)；于72℃反應(yīng)5 min后再4℃ 5 min，反應(yīng)終止（T：ARC1為60℃，Exo70A1為62℃），反應(yīng)結(jié)束后，置于4℃冰箱。PCR產(chǎn)物經(jīng)1%（W/V）瓊脂糖凝膠電泳，凝膠成像儀拍照，切膠，試劑盒回收純化產(chǎn)物。膠回收產(chǎn)物（ARC1、Exo70A1）分別進行TA克隆，連接到pMDTM19-T載體；重組載體轉(zhuǎn)化感受態(tài)細胞DH5α，LB固體培養(yǎng)基培養(yǎng)過夜，藍白斑篩選陽性克隆子，挑取陽性克隆并經(jīng)PCR鑒定后，分別抽提質(zhì)粒pMDTM19-T-ARC1、pMDTM19-T-Exo70A1，抽提質(zhì)粒送北京華大科技公司測序。將測序的結(jié)果進行拼接獲取結(jié)球甘藍材料的序列。

1.2.2 ARC1-Exo70A1互作區(qū)與非互作區(qū)編碼DNA的確定 不包括U-box區(qū)域的N端編碼序列（UND結(jié)構(gòu)域）是ARC1-Exo70A1的互作編碼區(qū)段［8，10，11］。結(jié)合NCBI在線網(wǎng)站查詢序列的結(jié)構(gòu)域，利用MEGA5.05軟件對ARC1序列和Exo70A1序列進行剪切，分別保留所需的互作編碼區(qū)和非互作編碼區(qū)于本地文件夾中。

1.2.3 結(jié)球甘藍ARC1和Exo70A1序列生物信息學(xué)分析 利用CLUATAL-W，Vector NTI Advance 10 進行多重序列比對分析；dnasp5軟件進行序列多態(tài)性分析。采用MEGA5.05構(gòu)建N-J模型進化樹，僅在構(gòu)建進化樹時將the bootstrap probabilities此項參數(shù)修改為1 000，其他步驟均采用默認設(shè)置。

1.2.4 ARC1和Exo70A1相似序列的獲取 NCBI及Brassica Database（http://brassicadb.org/）等網(wǎng)站分別獲取與結(jié)球甘藍ARC1和Exo70A1核酸或氨基酸序列相似性（Similarity）大于70%的ARC1（表2），Exo70A1核酸序列（表3）。結(jié)球甘藍所得材料序列命名為序列名稱加基因名，如E1-ARC1，E1-Exo70A1。對于同一種屬的不同序列在名稱后加上大寫字母A，B，C等加以區(qū)分，如序列1命名1Arabidopsis thaliana A，簡稱為1At-A。根據(jù)聚類情況來反映樹的可信度。

2 結(jié)果

2.1 結(jié)球甘藍ARC1與Exo70A1編碼區(qū)DNA的克隆、測序及序列分析

提取RNA后，經(jīng)反轉(zhuǎn)錄后獲得10種材料的cDNA。以cDNA為模板，分別擴增獲得相應(yīng)10種材料的ARC1和Exo70A1的CDS（圖1）；分別連接到pMDTM19-T載體進行測序，拼接后獲得所需序列數(shù)據(jù)，經(jīng)過Vector NTI分析可知，ARC1均為1 992 bp和 Exo70A1均為1 918 bp，與圖1中電泳結(jié)果條帶約在2 000 bp處相符。

圖1 ARC1和Exo70A1基因PCR電泳圖

10種材料的ARC1核酸序列相似度均高達97%以上，存在84個差異位點，氨基酸差異位點37個，且多為同義突變；10種材料的Exo70A1核酸序列相似度則均高于98%，有62個差異位點，氨基酸差異位點為23，也多為同義突變。圖2為ARC1與Exo70A1編碼區(qū)核酸序列與氨基酸序列局部比對圖（此部分為該序列差異位點最多的100 bp堿基區(qū)段），可以看出ARC1的核酸序列變異位點更多。綜上說明在自交不親和的結(jié)球甘藍中ARC1和Exo70A1序列高度相似，且Exo70A1比ARC1更為保守。

在線網(wǎng)站共收集獲得現(xiàn)有的83條ARC1和73

條Exo70A1本身及類似序列（表2，表3）。ARC1序列來自21個不同種屬，如蕓薹屬，擬南芥等；而Exo70A1則來自24個不同種屬。兩種基因的序列均可分為3大類：單子葉（水稻、小麥和玉米等），雙子葉植物（甘藍、白菜、擬南芥和大豆等）和低等植物（小立碗蘚和江南卷柏），但是未收集到江南卷柏的Exo701序列。序列比對結(jié)果表明：83條ARC1核酸序列相似度在41%-100%之間，存在677個差異位點；73條Exo70A1核酸序列相似度在43%-100%之間且大部分高于ARC1的相似度，有529個差異位點。

表2 ARC1序列GI號

表3 Exo70A1序列ACCESSION號

圖2 局部序列比對圖

2.2 ARC1-Exo70A1互作、非互作編碼區(qū)DNA的確定

在所獲得的73條Exo70A1的CDS序列中，其C-端均含有一個保守的Exo70結(jié)構(gòu)域。由文獻［11］可知：用ARC1 N端的UND功能域與Exo70A1的1-85氨基酸或1-270氨基酸互作均可激活報告基因，但與1-270氨基酸互作時其相互作用強度更大。據(jù)此認為所獲得的Exo70A1序列CDS的前840 bp為 Exo70A1蛋白與ARC1蛋白的主要互作編碼區(qū)。此外根據(jù)楊紅等［5］研究結(jié)果可知：RING-finger type結(jié)構(gòu)與U-box蛋白質(zhì)相似，所以認為其RING-finger結(jié)構(gòu)域功能類似于U-box結(jié)構(gòu)域；在NCBI網(wǎng)站預(yù)測結(jié)球甘藍ARC1結(jié)構(gòu)域，可知其存在U-box結(jié)構(gòu)域，且位于844-1 041 bp。

綜上可知ARC1-Exo70A1的互作區(qū)段（圖3），對序列進行剪切。ARC1以U-box保守域的第一個上游核苷酸為剪切點，對于無U-box的序列以RING-finger結(jié)構(gòu)域（表2中編號為43、45、47、48、52、53、58和61）的第一個上游核苷酸為剪切點，獲得93條互作區(qū)和非互作區(qū)ARC1序列；Exo70A1以序列840 bp處為剪切點，得到83條互作區(qū)和81條非互作區(qū)Exo70A1序列（其中編號為42和65的序列，因其非互作區(qū)不完整未使用）。

2.3 蕓薹屬與擬南芥屬序列的相對進化速率分析

10種材料結(jié)球甘藍序列相似性極高，利用其共有序列（命名為ARC1-C、Exo70A-C）進行相關(guān)分析。蕓薹屬及擬南芥屬序列進行多態(tài)性分析的結(jié)果，如表4和5所示。

非同義突變率（Ka）與同義突變率（Ks）的比值是選擇效應(yīng)和選擇方向的代表［16］。Ka/Ks大于1則說明基因受到正向選擇影響，為快速進化基因；Ka/Ks等于1則說明同義突變與非同義突變頻率相同，未受到選擇的影響；Ka/Ks小于1則說明基因受負向選擇影響，為相對保守基因。蕓薹屬與擬南芥屬的ARC1和Exo70A1序列的Ka與Ks及其比值均小于1（表4和表5），說明在蕓薹屬與擬南芥屬中，ARC1與Exo70A1這兩個基因都受負向選擇影響，是相對保守的基因。

表4 ARC1多態(tài)性分析

表5 Exo7A1多態(tài)性分析

比較不同物種同源或同類大分子的一級結(jié)構(gòu)可以得出其相對進化速率［17］，其公式為：

其中，K代表相對進化速率，2t代表進化時間，d代表替換的大分子數(shù)，N代表總的大分子數(shù)。

運用同義替換率估算分歧時間。序列間的同義突變（Ks）率如表4和5所示，且已知蕓薹屬與擬南芥的分歧發(fā)生在200萬年前（million years ago，MYA）［18］，根據(jù)上述公式可以得出K值：ARC1中Bo與At的同義突變率為1.478 7，得出KARC1為0.037 0［12］，同理可知KExo70A1為0.008 8，說明ARC1的相對進化速率比Exo70A1快，其保守性要低。而根據(jù)編碼區(qū)內(nèi)的互作區(qū)與非互作區(qū)不同區(qū)段序列的同義突變率可知，KARC1互作為0.043 1、KARC1非互作為0.031 1、KExo70A1互作為0.011 4、KExo70A1非互作為0.000 7；說明在ARC1-Exo70A1互作中，核酸序列互作區(qū)的進化速率比非互作區(qū)要快。

圖4 互作區(qū)ARC1的進化樹

2.4 10種材料的ARC1與Exo70A1的平行進化關(guān)系進化樹圖4-圖7中大部分同一種屬的序列位置均較近，說明此進化樹具有一定的可靠性。但個別同一種屬的序列分布在不同位置，，這可能是因為基因分化先于了物種分化。

2.4.1 ARC1-Exo70A1互作編碼區(qū)DNA的進化分析 互作編碼區(qū)ARC1進化樹（圖4-A）中單子葉植物與雙子葉植物均分為不同分支，而小立碗蘚和江南卷柏這兩個低等植物序列則在雙子葉部分獨立于一個分支下，說明ARC1的N端在進化過程中在雙子葉植物與單子葉植物中已經(jīng)發(fā)生了分化。結(jié)球甘藍序列與羽衣甘藍序列3聚類于同一樹枝下，說明其親緣關(guān)系很近。

Exo70A1互作區(qū)進化樹（圖5-A）的分支較短、大部分聚類到一處，說明Exo70A1保守性較高。該進化樹可大體分為3大部分：上部的雙子葉分支，中間的單子葉分支，下部的雙子葉分支。結(jié)球甘藍序列與甘藍型油菜序列9聚類到蕓薹屬的這一樹枝下，與實際情況植物分類相符合。從進化樹的分支情況可以看到蕓薹屬分支與擬南芥屬分支位置臨近，說明蕓薹屬與擬南芥屬的Exo70A1的N端在進化中是很保守的，具有很近的親緣關(guān)系。序列66即小立碗蘚這一低等植物聚類于雙子葉富集區(qū)。總體上來說Exo70A1在蕓薹屬植物的種間及屬間高度保守。

圖5 互作區(qū)Exo70A1的進化樹

根據(jù)ARC1互作區(qū)蕓薹屬部分放大圖（圖4-B）可知：蕓薹屬序列均聚類到同一樹枝下，說明蕓薹屬中ARC1的N-端保守性較好；擬南芥靠近，說明與擬南芥屬的親緣關(guān)系較為接近；蕓薹屬序列分為3部分，其中結(jié)球甘藍測序序列也分為3部分，但未發(fā)現(xiàn)其分支情況與親和指數(shù)之間的相關(guān)性。而觀察Exo70A1互作區(qū)的蕓薹屬部分（圖5-B）可發(fā)現(xiàn)：序列9甘藍型油菜與其他結(jié)球甘藍序列聚類在一起，這與ARC1互作區(qū)的分支情況存在差異；材料278，E1與序列7羽衣甘藍的親緣關(guān)系更近，并與其他結(jié)球甘藍材料聚類于不同分支下。

比較ARC1與Exo70A1互作編碼區(qū)DNA的進化樹，可發(fā)現(xiàn)3個主要的共同點：第一，單子葉植物與雙子葉植物的序列均分別聚類于不同分支之下，說明ARC1與Exo70A1的互作編碼區(qū)DNA均在單子葉植物與雙子葉植物之間已經(jīng)發(fā)生了分化；第二，蕓薹屬與擬南芥屬的序列位置都十分靠近，說明ARC1與Exo70A1的互作編碼區(qū)DNA在進化過程中至今仍然有著很近的親緣關(guān)系；第三，蕓薹屬序列分支情況類似；第四，低等植物序列均聚類于雙子葉富集區(qū)。綜上說明ARC1與Exo70A1的互作區(qū)存在協(xié)同進化。

2.4.2 ARC1-Exo70A1非互作編碼區(qū)DNA的進化分析 在ARC1與Exo70A1非互作編碼區(qū)進化樹（圖6、圖7）中，不論ARC1還是Exo70A1，蕓薹屬的非互作編碼區(qū)序列都在同一分支下；非互作區(qū)進化樹中雙子葉植物與單子葉植物也分別聚類于不同分支下，這與互作區(qū)進化樹情況相同。

ARC1非互作編碼區(qū)構(gòu)建的進化樹中（圖6）：與蕓薹屬親緣關(guān)系最近的是擬南芥（Arabidopsis），而大部分同一亞種的序列均聚類在同一分支下，表現(xiàn)出相同的亞種間同源性更高的情況；其蕓薹屬序列的分支情況與互作區(qū)部分相一致，而低等植物的序列則獨立于一個分支下，說明ARC1非互作區(qū)序列差異性更大；分支長度短于互作區(qū)，即進化速率慢。Exo70A1非互作編碼區(qū)構(gòu)建的進化樹中（圖7）：整體分支情況類似于互作區(qū)；分支長度也短于互作區(qū)；蕓薹屬序列未出現(xiàn)亞枝，說明其保守性非常好。

非互作區(qū)進化樹表明：蕓薹屬序列分支情況存在一定的差異性；其進化速率均低于互作區(qū)，與2.3中結(jié)果一致；低等植物的分支存在差異性。以上分析說明ARC1-Exo70A1互作中非互作區(qū)的DNA編碼序列比互作區(qū)比較而言更為保守。

圖6 ARC1非互作區(qū)進化樹

3 討論

本研究以自交不親和結(jié)球甘藍為研究材料，結(jié)合生物信息分析發(fā)現(xiàn)，結(jié)球甘藍ARC1與其他蕓薹屬ARC1序列相似性極高；結(jié)球甘藍Exo70A1與蕓薹屬Exo70A1序列也存在很高相似性及類似的結(jié)構(gòu)域，由此可知結(jié)球甘藍存在ARC1和Exo70A1的同源序列。

圖7 Exo70A1非互作區(qū)進化樹

所選用的10種結(jié)球甘藍的親和指數(shù)在0-1.28之間，其ARC1和Exo70A1的核酸序列無明顯差異存在，親緣關(guān)系很近。說明在高度自交不親和系的材料中ARC1和Exo70A1序列的保守性較好，具有保守的功能來維持自交不親和信號通路。因此，在實際SI系選育過程中，對于親和指數(shù)略高于1的品系不應(yīng)立即舍棄，可多觀察幾代自交不親和情況后再決定是否保留該品系。與蕓薹屬親緣關(guān)系較近的擬南芥屬中，擬南芥A. thaliana生態(tài)型均為自交親和型，而琴葉擬南芥A. lyrata卻有天然的自交不親和型［20］；但是在本文進化分析中發(fā)現(xiàn)ARC1和Exo70A1在蕓薹屬和擬南芥屬中卻具有很近的親緣關(guān)系和較高的相似度，但與禾本科等單子葉植物則進化差異較大，考慮到自交不親和機制出現(xiàn)在20 MYA之后，故可以推則ARC1在不同生態(tài)型擬南芥中對SI的貢獻不同，ARC1在擬南芥A. thaliana中較完整，在擬南芥A. lyrata中則變異加大。在禾本科等單子葉植物中則變異更大，可能承擔(dān)著其他功能。因此，所選育甘藍品系的親和指數(shù)高低可能反映著ARC1參與SI反應(yīng)的強弱程度，而Exo70A1似乎沒有這些差異。

本研究所用材料是經(jīng)過多代人工選育所得，10個品種的不結(jié)球甘藍因受人為因素控制，導(dǎo)致其變異趨同，如體現(xiàn)在10種材料的ARC1和Exo70A1的基因編碼序列具有極高的相似度，說明在同一選擇壓下基因的變異及進化會向著類似的方向進行。在Vekemans等［14］文獻中曾提到在環(huán)境條件的限制下，人工選育SI系，可能導(dǎo)致其自交親和性的加強。如干旱等不利條件或有利于其向SC系過渡，而順境中SI雜種優(yōu)勢可大大抵消其向SC系的轉(zhuǎn)變。在逆境條件下，人工選育的SI系其S位點多態(tài)性偏少，這可能是材料更為集中，基因差異較低的原因之一。因此人工選育SI系過程中應(yīng)該盡量保持良好環(huán)境，在獲得SI系純種的基礎(chǔ)上盡量減少自交次數(shù)。本研究所用材料雖是不同來源的親本，但其選擇范圍、選育條件等都可能在所測定的序列中起到了選擇作用，從而被固定為某種變異序列，導(dǎo)致10種自交不親和結(jié)球甘藍其ARC1和Exo70A1序列相似度高，存在平行進化關(guān)系，均聚類于蕓薹屬分支下。這說明短時期的逆境選育并未造成ARC1和Exo70A1編碼區(qū)DNA序列的改變，因此這種轉(zhuǎn)向是生理性的。但是選育自交不親和系最好在順境條件下進行，因為只有在順境條件下，SI所表現(xiàn)的異交優(yōu)勢，才不會被逆境掩蓋，自交不親和性也才能得以延續(xù)下去，也才具有選育SI系的遺傳和生理基礎(chǔ)。

ARC1和Exo70A1互作編碼區(qū)存在協(xié)同進化與其功能相適應(yīng)。本研究認為ARC1與Exo70A1蛋白質(zhì)的相互作用可能促成了ARC1-Exo70A1互作編碼區(qū)DNA的協(xié)同進化，即二者在進化過程中相互影響，從而具有相類似的進化方向。說明由于蛋白質(zhì)之間的相互作用需要不變的構(gòu)象，這可能在一定程度上約束了序列的變異程度，也就是說，當一方發(fā)生改變時，另一方也相應(yīng)會發(fā)生改變，這樣才能保持蛋白質(zhì)之間仍舊存在互作所需的構(gòu)象，最終形成協(xié)同進化這一趨勢。這與Goh等［21，22］提出的互作蛋白質(zhì)通常有相似的進化關(guān)系且表現(xiàn)為互作區(qū)域有相似的進化樹這一觀點相切合。此外，根據(jù)相對進化速率分析得知ARC1比Exo701的進化速度快，這與進化樹分析結(jié)果相一致，說明ARC1基因的變化性更大，其在植物中可能具有多樣性的功能。

ARC1-Exo70A1互作編碼區(qū)的變異速率比非互作編碼區(qū)更快。首先，在上游的SCR-SRK互作中，SRK的胞外域（S domain）起了主要作用［4-6］，然而S domain卻具有多態(tài)性，其氨基酸序列表現(xiàn)出一定的差異［12］。根據(jù)Kitashiba的文章［23］，可知不同胞外結(jié)構(gòu)域的SRK屬于不同S單倍型，且S單倍型的分化先于物種分化，那么SRK胞外結(jié)構(gòu)域與SCR的進化都應(yīng)受到單倍型約束，其分化也應(yīng)該先于物種分化；對于上游的蛋白質(zhì)而言，其可變性雖然較大，但還存在其他機制來保證反應(yīng)的進行，如SLG（S-glycoprotein，S-糖蛋白［24］）與同一單倍型SRK的S domain存在90%的氨基酸相似性［23］，它與SRK之間可能存在一個機制以保證SI反應(yīng)的發(fā)生。其次，在SRK-ARC1互作中ARC1的臂重復(fù)（ARM）區(qū)域（與Exo70A1的非互作區(qū)）是互作區(qū)，由本研究結(jié)果可知這個區(qū)段進化較慢，說明在SI途徑中SRKARC1互作占了比較重要的位置，需要較為保守的結(jié)構(gòu)域來維持反應(yīng)的進行。最后，在下游互作蛋白ARC1-Exo70A1中互作區(qū)段序列也具有很大的可變性，說明其可能只是SI下游反應(yīng)的一個分支，或存在其他機制來保障SI反應(yīng)的順利進行，而在楊佳等［10］的文章中提到破壞ARC1-Exo70A1途徑并不能完全打破SI，認為SI信號網(wǎng)絡(luò)可能存在其他的分支。此外，Exo70A1的Exo結(jié)構(gòu)域（與ARCI的非互作區(qū)）其進化慢，可能是與其他蛋白質(zhì)互作的區(qū)段，猜測其應(yīng)具有極其重要的作用；而目前已經(jīng)明確在酵母中Exo70結(jié)構(gòu)域能夠通過與Rho3 GTPase發(fā)生相互作用，來介導(dǎo)細胞的極性運輸［25］。Vekemans等［14］認為SI向SC的過渡是由許多獨立事件造成的，那么這些變異可能是分布于SI途經(jīng)中SCR、SRK、ARC1等相關(guān)蛋白的。短期內(nèi)，這些變異應(yīng)該是體現(xiàn)為隨機的序列變異，結(jié)合本研究，可知這些變異可能更加顯著地積累于重要蛋白的互作區(qū)核酸序列。

本研究分析了SI信號傳導(dǎo)途徑中蛋白質(zhì)ARC1和Exo-70A1，如果加入其他相關(guān)的互作蛋白進行分析，或許能發(fā)現(xiàn)SI途徑相關(guān)蛋白質(zhì)編碼DNA的進化規(guī)律，以期為其功能研究及SI信號傳導(dǎo)途徑研究提供參考。而目前已經(jīng)有許多軟件可以通過分子測年的方法來構(gòu)建系統(tǒng)發(fā)生樹的時序圖，那么下一步通過構(gòu)建時序圖便可獲得節(jié)點分歧時間，進而得到序列在各種物種中的相對進化速率，從而定量的分析不同大分子（核酸，蛋白質(zhì)等）的進化速率［19］，來說明SI途徑相關(guān)大分子在不同選擇壓力下的進化快慢，為SI的相關(guān)進化研究提供新的內(nèi)容。

4 結(jié)論

10種結(jié)球甘藍自交不親和材料之間存在平行進化關(guān)系；ARC1和Exo70A1在蕓薹屬和擬南芥屬中是相對較保守基因；在整個植物中ARC1和Exo70A1互作編碼區(qū)存在協(xié)同進化，但變異速率卻比非互作編碼區(qū)快。

［1］Schopfer CR, Nasrallah ME, Nasrallah JB. The male determinant of self-incompatibility in Brassica［J］. Science, 1999, 286（5445）：1697-1700.

［2］Stein JC, Howler B, Boyes DC, et al. Molecular cloning of a putative receptor protein kinase gene encoded at the self-incompatibility locus of Brassica oleracea［J］. Proc Natl Acad Sci USA, 1991, 88（19）：8816-8820.

［3］吳志剛. 甘藍自交不親和相關(guān)基因MOD的克隆及其與ARC1的FISH定位研究［D］. 重慶：西南大學(xué), 2011.

［4］Shimosato H, Yokota NH, Iwano M, et al. Characterization of the SP11/SCR high-affinity binding site involved in self/nonself recognition in Brassica self-incompatibility［J］. Plant Cell, 2007, 19（1）：107-117.

［5］楊紅, 朱利泉, 張賀翠, 等. 利用酵母雙雜交系統(tǒng)鑒定甘藍SCR與SRK胞外域片段間的相互作用［J］. 中國農(nóng)業(yè)科學(xué), 2011, 44（9）：1953-1962.

［6］Tantikanjana T, Nasrallah ME, Nasrallah JB. Complex networks of self-incompatibility signaling in the Brassicaceae［J］. Current Opinion in Plant Biology, 2010, 13（5）：520-526.

［7］Stone SL, Erson EM, Mullen RT, et al. ARC1 is an E3 ubiquitin ligase and promotes the ubiquitination of proteins during the rejection of self-incompatible Brassica pollen［J］. Plant Cell, 2003, 15（4）：885-898.

［8］楊昆, 張賀翠, Richard C, 等. 甘藍自交不親和信號轉(zhuǎn)導(dǎo)元件ARC1與Exo70A1的相互作用［J］. 作物學(xué)報, 2011, 37（12）：2136-2144.

［9］Samuel MA, Chong YT, Haasen KE, et al. Cellular pathways regulating responses to compatible and self-incompatible pollen in Brassica and Arabidopsis stigmas intersect at Exo70A1, a putative component of the exocyst complex［J］. Plant Cell, 2009, 21（7）：2655-2671.

［10］楊佳, 李玉花, 藍興國. 蕓苔屬植物自交不親和性S-受體激酶的內(nèi)吞作用及信號傳遞網(wǎng)絡(luò)［J］. 植物生理學(xué)報, 2012, 48（3）：211-216.

［11］楊佳. 羽衣甘藍ARC1相互作用蛋白及授粉過程中差異表達蛋白的研究［D］. 哈爾濱：東北林業(yè)大學(xué), 2012.

［12］Ryuichi I, Takeshi N. Phylogenetic analysis of Brassiceae based on the nucleotide sequences of the S-locus related gene, SLR1［J］. Theoretical and Applied Genetics, 2002, 105（8）：1159-1165.

［13］Edh K, Widén B, Ceplitis A. The evolution and diversification of S-locus haplotypes in the Brassicaceae family［J］. Genetics, 2009, 181（7）：977-984.

［14］Vekemans X, Poux C, Goubet PM, et al. The evolution of selfing from outcrossing ancestors in Brassicaceae：what have we learned from variation at the S-locus?［J］. Evol Biol, 2014, 27（7）：1372-1385.

［15］ 張桂玲. 甘藍自交不親和性快速鑒定方法的研究［D］. 哈爾濱：東北農(nóng)業(yè)大學(xué), 2003.

［16］Nielsen R. Molecular signatures of natural selection［J］. Annual Review of Genetics, 2005, 39：197-218.

［17］張昀. 生物進化［M］. 北京：北京大學(xué)出版社, 1998：207.

［18］Koch M, Haubold B, Mitchell-Olds T. Molecular systematic of the Brassicaceae：evidence from coding plastidic MatK and nuclear Chs sequences［J］. Am J Bot, 2001, 88（3）：534-544.

［19］ 黃原. 分子系統(tǒng)發(fā)生學(xué)（第1版）［M］. 北京：科學(xué)出版社, 2012：1-21.

［20］郝艾馨, 藍興國, 王宇, 等. SRK-SCR轉(zhuǎn)基因擬南芥自交不親和性的研究進展［J］. 植物生理學(xué)報, 2013, 11：1113-1120.

［21］Goh CS, Bogan AA, Joachimiak M, et al. Co-evolution of proteins with their interaction partners［J］. Mol Biol, 2000, 299（2）：283-293.

［22］Goh CS, Cohen FE. Co-evolutionary analysis reveals insights into protein-protein interactions［J］. Mol Biol, 2002, 1：177-192.

［23］Kitashiba H, Nasrallah JB. Self-incompatibility in Brassicaceae crops：lessons for interspecific incompatibility［J］. Breed Sci, 2014, 64（1）：23-37.

［24］Watanabe M, Takasaki T, Toriyama K, et al. A high degree of homology exists between the protein encoded by SLG and the S receptor domain encoded by SRK in self-incompatible Brassica campestris L.［J］. Plant Cell Physiol, 1994, 8：1221-1229.

［25］Robinson NG, Guo L, Imai J, et al. Rho3 of Saccharomyces cerevisiae, which regulates the actin cytoskeleton and exocytosis, is a GTPase which interacts with Myo2 and Exo70［J］. Mol Cell Biol, 1999, 19（5）：3580-3587.

（責(zé)任編輯 李楠）

Cloning and Phylogenetic Analysis of cDNA Sequences Coding for ARC1 and Exo70A1 in Brassica oleraces var. capitata L.

YANG Dan LIAN Xiao-ping ZHOU Yan ZHANG He-cui DU Dan GAO Qi-guo REN Xue-song ZHU Li-quan
（ Laboratory of Physiology and Biochemistry，Southwest University，Chongqing 400715）

Interaction between ARC1 and Exo70A1 in Brassica is a key point of self-incompatibility（SI）signal transduction. In order to study the evolutionary relationship of ARC1 and Exo70A1 from 10 cultivars of SI head cabbages developed in recent years，we successfully cloned and sequenced cDNAs of ARC1 and Exo70A1. Furthermore，by bioinformatics，we had mutation and evolutionary analysis of those coding region sequences，including interaction region and non-interaction regions of ARC1-Exo70A1 from the 10 cabbages and prior published 83 ARC1 and 73 Exo70A1 sequences. The results revealed that：1）The coding nucleic acid sequences of ARC1 and Exo70A1 from the 10 cabbage materials presented parallel evolutionary relationships，and belonged to the Brassica branch. 2）The evolution rate of coding sequence of ARC1 was faster than Exo70A1. 3）The interaction coding region of ARC1-Exo70A1 was obviously in the progress of co-evolution，while the non-interaction coding area of ARC1-Exo70A1 showed significant differentiation rate in evolution. 4）The evolutionary rate of interaction coding area was faster than non-interaction coding region.

ARC1；Exo70A1；self-incompatibility；evolutionary relationship；co-evolution

10.13560/j.cnki.biotech.bull.1985.2016.06.012

2015-10-07

國家自然科學(xué)基金項目（31572127），重慶市自然基金重點項目（cstc2012jjB80010）

楊丹，女，碩士研究生，研究方向：分子作物育種；E-mail：yangdan042561@163.com

朱利泉，教授，研究方向：植物分子生物學(xué)；E-mail：zhuliquan@swu.edu.cn