孫琦，黃薇，劉芳，金玉環(huán)，肖建旺，黃先忠

(石河子大學生命科學學院,特色植物基因組學實驗室，新疆 石河子832003)

植物在生長過程中不斷面臨著各種生物和非生物因素的威脅，鹽脅迫已成為影響全球農(nóng)作物生產(chǎn)的主要非生物因素之一[1-5]。鈉離子的過量積累造成細胞膜透性增加、蛋白功能紊亂、有毒物質(zhì)積累及氧化損傷等，嚴重影響植物的生長發(fā)育[6-7]。植物通過不斷適應(yīng)高鹽的環(huán)境而進化出了一些有效的策略以保證生存[8-10]。在細胞水平上，通過將Na+有效分隔到液泡和其他細胞器官中以及將Na+從根部排出是避免細胞質(zhì)Na+濃度過高和維持細胞質(zhì)K+/Na+比值的有效機制[11-13]。

轉(zhuǎn)運蛋白(transport protein)是膜蛋白的一大類，介導(dǎo)生物膜內(nèi)外的化學物質(zhì)及信號的交換。因脂質(zhì)雙分子層在細胞或細胞器周圍形成了一道疏水屏障，導(dǎo)致大部分親水性化合物，如糖、氨基酸、離子等都需要特定轉(zhuǎn)運蛋白的幫助來通過疏水屏障。因此，轉(zhuǎn)運蛋白在營養(yǎng)物質(zhì)攝取、代謝產(chǎn)物釋放以及信號轉(zhuǎn)導(dǎo)等廣泛的細胞活動中起著重要的作用[14]。已有研究[15]表明植物體內(nèi)存在多個與Na+轉(zhuǎn)運相關(guān)的蛋白，如液泡膜上的Na+/H+逆向轉(zhuǎn)運蛋白(NHX)和細胞質(zhì)膜上的Na+/H+逆向轉(zhuǎn)運蛋白(SOS)。液泡膜NHX主要行使Na+在液泡中區(qū)域化的功能，在提高植物耐鹽性方面發(fā)揮了重要作用[16]。He等[17]將擬南芥AtNHX1轉(zhuǎn)入棉花后發(fā)現(xiàn)，轉(zhuǎn)基因能在高濃度NaCl條件下正常生長。Wu等[18]將水稻液泡型Na+/H+逆向轉(zhuǎn)運蛋白OsNHX1基因轉(zhuǎn)入獐毛中，轉(zhuǎn)基因獐毛植株能夠在350 mM NaCl處理10周后正常存活。另外，過量表達NHX基因在轉(zhuǎn)基因植物擬南芥(Arabidopsisthaliana)[19]，番茄(Lycopersiconesculentum)[20]，茄果(Solanumpennellii)[21]和水稻(Oryzasativa)[2]中也證實了耐鹽性的改善。

十字花科短命植物小擬南芥(Arabidopsispumila)，是十字花科擬南芥的近源物種，主要生長在干旱少雨、土壤鹽堿化嚴重的新疆的北部，具有耐鹽堿、耐干旱等特點[22]。近年來，我們從小擬南芥中克隆了一些抗逆基因[23-29]，但目前對于小擬南芥NHX2基因的克隆還沒有研究，其響應(yīng)鹽堿脅迫的機制還不清楚。

近年來利用轉(zhuǎn)錄組測序技術(shù)鑒定耐鹽相關(guān)的基因及功能分析，更好的解析了植物的耐鹽機制[30-31]。前期我們利用二代和三代測序技術(shù)測序了小擬南芥響應(yīng)高鹽脅迫的全長轉(zhuǎn)錄組，經(jīng)分析發(fā)現(xiàn)了8075個差異表達基因，初步分析了小擬南芥的鹽適應(yīng)機制[32]。本研究中，我們基于小擬南芥轉(zhuǎn)錄組數(shù)據(jù)篩選出了1157個鹽脅迫下差異表達的轉(zhuǎn)運蛋白基因，并通過K均值(K-means)方法對轉(zhuǎn)運蛋白基因進行聚類，分析了轉(zhuǎn)運蛋白基因在20個子聚類中的表達特征，從上調(diào)表達的子聚類中篩選并克隆了一個Na+/H+逆向轉(zhuǎn)運蛋白基因ApNHX2。本研究所獲得了差異表達的轉(zhuǎn)運蛋白基因，為今后深入研究小擬南芥轉(zhuǎn)運蛋白基因在耐鹽性方面的作用提供了理論依據(jù)，為探索短命植物適應(yīng)環(huán)境的分子機理奠定了基礎(chǔ)。

1 材料與方法

1.1 植物材料

新疆小擬南芥種子為本實驗室保存。種子萌發(fā)使用1/2 MS固體培養(yǎng)基。種子消毒參照文獻[23]所述方法進行。

1.2 實驗方法

1.2.1 小擬南芥脅迫處理、總RNA的提取和cDNA的合成

小擬南芥種子經(jīng)滅菌后均勻播種在1/2 MS培養(yǎng)基上，4 ℃春化處理3 d后，放置于長日照光照培養(yǎng)箱中培養(yǎng)(16 h光照/8 h黑暗,溫度為22 ℃)。將生長2周后的小擬南芥幼苗轉(zhuǎn)移至含250 mM NaCl的1/2 MS培養(yǎng)基上分別處理0、0.5、3、6、12、24和48 h后收集材料。另外收集小擬南芥成熟植株的根、莖、葉、花以及果莢組織材料。采用TIANGEN公司的RNA prep pure Plant Kit試劑盒按照說明書提取所有樣品的RNA。根據(jù)北京百泰克公司的M-MLV Reverse Transcriptase合成cDNA模板第一鏈。

1.2.2 小擬南芥響應(yīng)鹽脅迫的轉(zhuǎn)運蛋白基因分析

根據(jù)七大數(shù)據(jù)庫(Nr，Nt，Pfam，KOG/COG，Swiss-prot，KEGG，GO)的基因功能注釋結(jié)合TransportDB2.0(http://www.membranetransport.org/transportDB2/index.html/)數(shù)據(jù)庫，以transport為關(guān)鍵詞在小擬南芥轉(zhuǎn)錄組測序得到的8075個差異表達基因功能注釋信息文件[32]中手動查找轉(zhuǎn)運蛋白基因，按功能進行細化分類。將查找到的所有轉(zhuǎn)運蛋白基因在基因表達量注釋信息文件中進行搜素，記錄每個基因在鹽脅迫處理不同時間點的表達量，以FPKM(expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced，F(xiàn)PKM：每百萬堿基對每千堿基轉(zhuǎn)錄物序列的預(yù)期片段數(shù))值表示。利用R studio(https://www.rstudio.com/)繪制轉(zhuǎn)運蛋白基因的表達熱圖。

對篩選得到的轉(zhuǎn)運蛋白基因以基因的相對表達水平值log2(ratios)進行K-means聚類分析并統(tǒng)計分布在每個subcluster中的轉(zhuǎn)運蛋白基因的數(shù)量。對于每個基因來說ratios為該基因在樣品2的FPKM與該基因在樣品1的FPKM值的比值。采用相應(yīng)的距離算法計算每個基因之間的距離，根據(jù)基因的相對距離遠近分成不同的subcluster，從而實現(xiàn)聚類。

1.2.3 Na+/H+逆向轉(zhuǎn)運蛋白基因的篩選

根據(jù)K-means聚類分析得到的20個subcluster中，發(fā)現(xiàn)有一類整體趨勢表現(xiàn)為上調(diào)的K subcluster，查找分布于該K subcluster中的基因表達量log10(FPKM+1)值，利用R studio(https://www.rstudio.com/)繪制熱圖，以直觀分析各個基因在不同時間點的表達特征。將分布于該上調(diào)聚類中的Na+轉(zhuǎn)運蛋白基因序列在CDS預(yù)測信息文件中查找編碼區(qū)序列，對每1條基因在tair數(shù)據(jù)庫(http://www.arabidopsis.org/)中進行Blast比對，并記錄比對后的結(jié)果。根據(jù)比對結(jié)果發(fā)現(xiàn)Na+轉(zhuǎn)運蛋白基因中有1條與擬南芥AtNHX2基因序列相似性極高的同源基因ApNHX2，在鹽脅迫處理后的各個時間點的表達量均顯著上調(diào)。

1.2.4ApNHX2基因的克隆和測序

使用DNAStar分析ApNHX2完整開放閱讀框(ORF)，根據(jù)ORF利用Primer5.0軟件設(shè)計引物ApNHX2-F和ApNHX2-R，在兩端分別加上KpnI和XbaI酶切位點(表1)。以250 mM NaCl脅迫處理12 h后的2周大的小擬南芥幼苗cDNA第一鏈為模板進行RT-PCR擴增，在20 μL的反應(yīng)體系中加入10×PCR Buffer 2 μL，dNTP 1.5 μL，上下游引物各0.5 μL，模板cDNA 2 μL，Ex Taq (TaKaRa)酶0.3 μL，ddH2O 13.2 μL。反應(yīng)程序為94 ℃預(yù)變性2 min，94 ℃變性30 s，60 ℃退火30 s，72 ℃延伸2 min，共38個循環(huán)，終延伸為72 ℃ 10 min。PCR產(chǎn)物在1% (w/v)的瓊脂糖凝膠上檢測，將PCR產(chǎn)物利用凝膠回收試劑盒回收純化后連接到pMD19-T(TaKaRa)載體上，轉(zhuǎn)化大腸桿菌DH5α感受態(tài)細胞，經(jīng)藍白斑篩選，挑取單克隆于液體LB培養(yǎng)基中37 ℃、200 rpm下振蕩培養(yǎng)過夜，提取質(zhì)粒經(jīng)雙酶切鑒定，將鑒定正確的質(zhì)粒交由北京華大基因公司測序。

1.2.5 ApNHX2蛋白的跨膜結(jié)構(gòu)分析與進化樹構(gòu)建

利用TMHMM server v.2.0(http://www.cbs.dtu.dk/services/TMHMM/)分析該蛋白的跨膜區(qū)。利用NCBI的Blastp程序檢索ApNHX2的同源蛋白，氨基酸序列的多重比對采用Clustal X2.0程序完成，用MEGA5.0軟件構(gòu)建系統(tǒng)進化樹，其中用Neighbor-Joining方法進行1000次Bootstrap分析[33]。

1.2.6ApNHX2基因的組織表達分析

采用qRT-PCR進行基因的組織表達分析。ApNHX2基因引物為qRT-ApNHX2-F和qRT-ApNHX2-R，選擇GAPDH作為內(nèi)參基因[35]，引物為GAPDH-F和GAPDH-R (表1)。qRT-PCR采用Fast SYBR Mixture試劑盒(康為世紀)，利用美國ABI 7500 Fast實時熒光定量PCR儀(Life Technologies,Foster City,CA,USA)進行擴增。在10 μL的反應(yīng)體系中加入cDNA模板2 μL，UltraSYBR Mixture 5 μL，ROX 0.2 μL，正反向引物(10 μM)各0.2 μL，ddH2O 2.6 μL。反應(yīng)程序為95 ℃預(yù)變性10 min；95 ℃變性15 s，60 ℃退火1 min，共40個循環(huán)；60 ℃讀取熒光值。檢測每份樣品的ApNHX2基因和GAPDH內(nèi)參基因的Ct值，每份樣品3次重復(fù)，實驗數(shù)據(jù)采用2-ΔΔCt法進行相對定量分析[23]。采用SPSS17.0軟件對實驗結(jié)果進行單因素方差分析，利用Duncan法進行多組樣本之間的差異顯著性分析(P<0.05)[32]。

表1 本研究所用的PCR引物序列Tab.1 PCR primers used in this study

注：下劃線代表限制性內(nèi)切酶的酶切位點。

1.2.7 過量表達載體構(gòu)建

將測序正確的質(zhì)粒pMD19-T∶ApNHX2以及載體pCAMBIA2300-CaMV35S-OCS使用KpnI和XbaI進行雙酶切，電泳檢測后回收目的片段。利用T4連接酶4 ℃下過夜連接，將連接后的產(chǎn)物轉(zhuǎn)化DH5α感受態(tài)細胞，37 ℃下倒置培養(yǎng)14～16 h。挑取單菌落于液體LB培養(yǎng)基中搖瓶過夜培養(yǎng)，通過抽提質(zhì)粒和進行酶切驗證，利用1.2% (w/v)的瓊脂糖電泳檢測，構(gòu)建35S∶ApNHX2載體。將鑒定正確的質(zhì)粒，用凍融法轉(zhuǎn)化農(nóng)桿菌GV3101，挑取單克隆進行菌液PCR驗證，將鑒定正確的單克隆菌液甘油保存于-80 ℃?zhèn)溆谩?/p>

2 結(jié)果與分析

2.1 篩選小擬南芥中耐鹽相關(guān)的轉(zhuǎn)運蛋白基因

根據(jù)七大數(shù)據(jù)庫功能注釋從8075個差異表達基因中找到1157個耐鹽相關(guān)的轉(zhuǎn)運蛋白基因，基因表達熱圖分析表明在鹽脅迫的早期階段(0.5～6 h)，一些基因被快速激活上調(diào)表達，一些基因在脅迫的后期(12～48 h)被激活上調(diào)表達，而另外一些基因則表現(xiàn)為隨脅迫時間的增加表達量持續(xù)上升(圖1 A)。將1157個轉(zhuǎn)運蛋白基因按功能分為ABC轉(zhuǎn)運蛋白(ABC transporter)；Ca2+轉(zhuǎn)運蛋白(Calcium-transporter)；K+轉(zhuǎn)運蛋白(Potassium transporter)；Na+轉(zhuǎn)運蛋白(Sodium- transporter)以及糖轉(zhuǎn)運蛋白(Sugar transporter)等(圖1B)。統(tǒng)計主要轉(zhuǎn)運蛋白基因的數(shù)量包括24個Na+轉(zhuǎn)運蛋白基因，36個K+轉(zhuǎn)運蛋白基因，30個Ca2+轉(zhuǎn)運蛋白基因，9個Cl-轉(zhuǎn)運蛋白基因以及ATPase相關(guān)轉(zhuǎn)運蛋白家族H+-PPase、F-ATPase和P-ATPase中所含基因的數(shù)量分別為70、2和80個，從中得出P-ATPase家族中所包含的基因數(shù)量最多，而Cl-轉(zhuǎn)運蛋白相關(guān)基因數(shù)量最少(圖1B)。

A：1157個轉(zhuǎn)運蛋白基因在250 mM NaCl脅迫不同時間的表達熱圖；B：主要轉(zhuǎn)運蛋白基因的數(shù)量分布情況圖1 1157個鹽脅迫差異表達的轉(zhuǎn)運蛋白基因的分析Fig.1 Analysis of the 1157 differently expressed transporter genes under salt stress

2.2 K-means聚類方法分析轉(zhuǎn)運蛋白基因

K-means聚類分析結(jié)果顯示，1157個轉(zhuǎn)運蛋白基因分布于20個K subcluster中，其中在K6，K9以及K15子聚類中的基因數(shù)量分布較多，分別為172，193以及190個，而在K5子聚類中基因分布最少，數(shù)量為5個(圖2 A)。根據(jù)log2(ratio)值繪制的差異基因表達量聚類圖表明分布于K6子聚類中的基因相對表達水平整體表現(xiàn)出持續(xù)上調(diào)的趨勢(圖2B)。通過查找分布于K6子聚類中的基因表達量log10(FPKM+1)值所繪制的熱圖結(jié)果顯示，經(jīng)鹽脅迫處理后的基因一部分表現(xiàn)為早期(0.5～12 h)持續(xù)上調(diào)表達，一部分表現(xiàn)為后期(12～48 h)持續(xù)上調(diào)表達，而有一部分基因持續(xù)上調(diào)表達(圖2C)。有一個差異表達明顯的Na+轉(zhuǎn)運蛋白基因ApNHX2分布在K6子聚類中，隨著處理時間的增加，ApNHX2基因的表達量逐漸升高且在處理12 h時表達量達到最高，之后則表現(xiàn)出下降的趨勢，但ApNHX2基因在6個脅迫時間點的表達量均高于對照條件下(0 h)的表達量(圖2D)，暗示該基因在小擬南芥耐鹽反應(yīng)中起著重要作用。

A：1157個轉(zhuǎn)運蛋白基因在20個子聚類中的數(shù)量分布：B：K6差異基因表達量聚類圖；C：K6子聚類中172個轉(zhuǎn)運蛋白基因 響應(yīng)250 mM NaCl脅迫的表達熱圖，箭頭所指表示差異表達顯著的Na+/H+逆向轉(zhuǎn)運蛋白基因ApNHX2； D：ApNHX2基因在250 mM NaCl脅迫不同時間的表達圖2 K-means聚類分析Fig.2 K-means clustering analysis

2.3 ApNHX2基因的克隆及生物信息學分析

利用引物ApNHX2-F和ApNHX2-R擴增得到約1.6 kb的擴增產(chǎn)物。測序結(jié)果結(jié)合全長轉(zhuǎn)錄組數(shù)據(jù)庫分析表明ApNHX2全長為2299 bp，ORF大小為1626 bp，推測其編碼541個氨基酸(圖3)。ApNHX2蛋白的跨膜結(jié)構(gòu)預(yù)測結(jié)果顯示ApNHX2是一個典型的跨膜轉(zhuǎn)運蛋白，具有12個跨膜結(jié)構(gòu)區(qū)，其中包括2個未完全跨膜區(qū)(圖4)。

2.4 ApNHX2蛋白的系統(tǒng)進化分析

通過NCBI數(shù)據(jù)庫檢索ApNHX2的同源蛋白序列，利用MEGA5.0軟件構(gòu)建了小擬南芥ApNHX2與擬南芥(Arabidopsisthaliana)，野甘藍(Brassicaoleraceavar.oleracea)，歐洲油菜(Brassicanapus)，鹽角草(Salicorniaeuropaea)，木薯(Manihotesculenta)，玫瑰(Rosarugosa)，醉蝶花(Tarenayahassleriana)，榴蓮(Duriozibethinus)，濱藜(Atriplexhalimus)，橡膠樹(Heveabrasiliensis)，白樺(Betulaplatyphylla)以及胡楊(Populuspruinosa)中NHX2蛋白的系統(tǒng)進化樹(圖5)。結(jié)果表明，ApNHX2與擬南芥以及歐洲油菜等的NHX2蛋白進化關(guān)系均較近，與擬南芥NHX2屬于同一進化分支。

2.5 ApNHX2基因的組織表達分析

組織表達分析結(jié)果顯示，ApNHX2基因在小擬南芥各組織中均有表達，在花中表達量最高，在其它組織中表達量均較低且無明顯差異(圖6)。

2.6 植物過量表達載體35S∶ApNHX2的構(gòu)建

經(jīng)KpnI和XbaI酶切的pMD19-T∶ApNHX2質(zhì)粒以及pCAMBIA2300-p35s-OCS載體用T4連接酶連接后轉(zhuǎn)化大腸桿菌DH5α，抽提質(zhì)粒并進行雙酶切驗證，瓊脂糖凝膠電泳檢測結(jié)果顯示，質(zhì)粒1和質(zhì)粒2酶切后條帶均與目的基因大小一致，均在接近1600 bp處，說明目的片段已連接到表達載體上，成功構(gòu)建植物過量表達載體35S∶ApNHX2(圖7 A)。最終載體35S∶ApNHX2含有花椰菜花葉病毒CaMV35S啟動子，章魚堿合成酶基因OCS終止子及植物篩選標記新霉素磷酸轉(zhuǎn)移酶基因NPTⅡ (圖7B)。將35S∶ApNHX2載體轉(zhuǎn)化農(nóng)桿菌GV3101感受態(tài)細胞，通過菌液PCR驗證發(fā)現(xiàn)擴增出約1600 bp長度的條帶，表明該載體已成功轉(zhuǎn)化農(nóng)化菌GV3101(圖7C)。

圖3 ApNHX2的核苷酸序列及推測氨基酸序列Fig.3 Nucleotide sequence and putative amino acid sequence of ApNHX2

紅色框代表2個未完全跨膜區(qū)；綠色框代表10個強跨膜區(qū)圖4 ApNHX2蛋白的跨膜結(jié)構(gòu)預(yù)測Fig.4 Prediction of transmembrane structure of ApNHX2

▲代表本研究克隆的ApNHX2蛋白；圖中分支點的數(shù)字表示Bootstrap驗證中基于1000次 重復(fù)該節(jié)點可信度的百分比；標尺代表每一位點堿基替換率圖5 植物Na+/H+逆向轉(zhuǎn)運蛋白NHX2的系統(tǒng)進化樹分析Fig.5 Phylogenetic tree of Na+/H+ antiporter NHX2 from different plants

小寫字母表示差異顯著(P<0.05)圖6 qRT-PCR驗證的ApNHX2基因在在 根、莖、葉、花和果莢中的表達Fig.6 Expression pattern of ApNHX2 in different tissues of Arabidopsis pumila

A：ApNHX2基因的質(zhì)粒酶切檢測，M為marker Ⅲ，1,2為35S∶ApNHX2；B：35S∶ApNHX2載體模式圖；C：ApNHX2基因的 農(nóng)桿菌PCR檢測，M為marker Ⅲ，1為ApNHX2陽性對照，2-5為ApNHX2，6為ApNHX2陰性對照圖7 過表達載體35S∶ApNHX2的構(gòu)建Fig.7 Construction of plant overexpression vector 35S∶ApNHX2

3 討論

轉(zhuǎn)錄組分析在闡明各種生長和發(fā)育過程以及應(yīng)激條件下基因表達調(diào)控的復(fù)雜性中起著重要作用[35]。迄今為止的大多數(shù)研究[36-38]都報道了由二代測序(second-generation sequencing,SGS)技術(shù)產(chǎn)生的短讀序列用于分析植物組織的轉(zhuǎn)錄組，但是這些序列阻止了全長轉(zhuǎn)錄本組裝，降低了序列組裝的準確性，使生物信息學分析變得困難。單分子實時(single-molecule real-time,SMRT)測序平臺提供的第三代測序(third-generation sequencing,TGS)技術(shù)可以克服上述限制，由于其長讀取而被廣泛用于基因組測序[39-40]。盡管TGS具有相對較高的錯誤率，但通過使用短而高精度的SGS讀數(shù)進行校正可以克服這一缺點[37-38]。在本實驗室的前期研究中，我們結(jié)合SGS和TGS技術(shù)測序以產(chǎn)生更完整的A.pumila轉(zhuǎn)錄組[32]，作為本研究的基礎(chǔ)。

通過在小擬南芥轉(zhuǎn)錄組數(shù)據(jù)庫中手動查找和篩選，從8075個差異表達基因中共整理出1157個功能注釋為轉(zhuǎn)運蛋白的基因，包括24個小擬南芥Na+轉(zhuǎn)運蛋白基因(圖1)。目前隨著擬南芥和楊樹(Populuspruinosa)的轉(zhuǎn)錄組測序工作的完成，分析顯示擬南芥中存在10個Na+轉(zhuǎn)運蛋白基因[41]，楊樹中存在11個Na+轉(zhuǎn)運蛋白基因[42]，研究者們對植物中的Na+轉(zhuǎn)運蛋白基因有了進一步了解。本研究篩選出的小擬南芥Na+轉(zhuǎn)運蛋白基因數(shù)量已遠遠超過擬南芥和楊樹，因此我們推測Na+轉(zhuǎn)運蛋白在小擬南芥適應(yīng)高鹽環(huán)境中起著重要的作用。

邱生平等[43]采用RT-PCR的方法從水稻中克隆了液泡膜Na+/H+逆向轉(zhuǎn)運蛋白基因OsNHX2的全長cDNA，包含1635 bp的開放閱讀框，編碼544個氨基酸。李中虎等[44]在轉(zhuǎn)錄組測序的基礎(chǔ)上，得到了濱麥(Leymusmollis)液泡膜Na+/H+逆向轉(zhuǎn)運蛋白基因LmNHX2的DNA序列，包含1617 bp的開放閱讀框，編碼538個氨基酸。本研究基于轉(zhuǎn)錄組數(shù)據(jù)得到小擬南芥Na+/H+逆向轉(zhuǎn)運蛋白基因ApNHX2的cDNA序列，包含1626 bp的全長開放閱讀框，編碼541個氨基酸(圖3)。水稻OsNHX2和濱麥LmNHX2的跨膜結(jié)構(gòu)域分析發(fā)現(xiàn)二者均含有12個跨膜結(jié)構(gòu)區(qū)，其中包括2個未完全跨膜區(qū)[43-44]。本研究對小擬南芥ApNHX2跨膜結(jié)構(gòu)域分析發(fā)現(xiàn)與水稻OsNHX2和濱麥LmNHX2具有相似的結(jié)果，含有12個跨膜結(jié)構(gòu)區(qū)，其中包括2個未完全跨膜區(qū)(圖4)。另外，邱生平等[43]通過系統(tǒng)發(fā)育分析表明OsNHX2與液泡膜Na+/H+逆向轉(zhuǎn)運蛋白的同源關(guān)系較近而與質(zhì)膜Na+/H+逆向轉(zhuǎn)運蛋白的同源關(guān)系較遠，這些結(jié)果說明OsNHX2可能是液泡膜Na+/H+逆向轉(zhuǎn)運蛋白。而基于LmNHX2和其他NHX蛋白氨基酸序列構(gòu)建的系統(tǒng)進化樹[44]顯示，濱麥LmNHX2與小麥(Triticumaestivum) TaNHX2二者氨基酸序列相似度高達99%，這與濱麥是小麥的野生近緣種的關(guān)系相符合。本研究通過系統(tǒng)發(fā)育分析表明ApNHX2與擬南芥NHX2蛋白進化關(guān)系較近，屬于同一進化分支，這與小擬南芥是擬南芥的野生近緣種的關(guān)系相符合，可推斷出ApNHX2屬于液泡膜Na+/H+逆向轉(zhuǎn)運蛋白(圖5)。

基因表達模式的分析對于功能研究必不可少，組織表達分析表明ApNHX2在小擬南芥各組織中均有表達，但在花中表達量最高且明顯高于其他組織，說明ApNHX2很可能參與到小擬南芥花器官發(fā)育的過程(圖6)。已有研究[15]表明洋芋(Helianthustuberosus) Na+/H+逆向轉(zhuǎn)運蛋白基因HtNHX2的表達提高了水稻對鹽度的耐受性。在K+有限的鹽脅迫條件下或一般營養(yǎng)缺乏的條件下，HtNHX2的表達增加了水稻籽粒產(chǎn)量，收獲指數(shù)以及總養(yǎng)分吸收[15]。小麥Na+/H+逆向轉(zhuǎn)運蛋白基因TaNHX2的表達也表現(xiàn)為能提高苜蓿(MedicagoSativa)[45]和大豆(Glycinemax)[46]的耐鹽性。我們的研究發(fā)現(xiàn)，隨著高鹽脅迫時間的增加，ApNHX2基因的表達量逐漸升高且在處理12 h時表達量達到最高，之后表現(xiàn)出下降趨勢，但ApNHX2基因在6個脅迫時間點的表達量均高于對照條件下(0 h)的表達量(圖2D)，暗示該基因在小擬南芥耐鹽反應(yīng)中起著重要作用。本研究構(gòu)建了植物過量表達載體35S∶ApNHX2，并成功轉(zhuǎn)化農(nóng)桿菌GV3101，我們將進一步驗證ApNHX2基因功能，為今后深入研究小擬南芥轉(zhuǎn)運蛋白基因在耐鹽性方面的作用提供理論依據(jù)。