摘要： 鹽芥是研究耐鹽機理的模式植物。為從蛋白質水平揭示鹽芥響應鹽脅迫的分子機制，本研究采用同位素標記相對和絕對定量(iTRAQ) 技術對不同NaCl濃度處理7 d的鹽芥葉片進行差異蛋白質組學分析。結果表明，鹽芥葉片中共鑒定到4 607個蛋白質，其中281個蛋白質的表達豐度顯著增加，95個蛋白質的表達豐度顯著降低。鹽脅迫差異表達蛋白質的KEGG代謝通路和蛋白質互作網絡分析結果表明，促進植株光合作用可幫助鹽芥適應低鹽環境;抑制葉綠素和支鏈氨基酸合成、調控應激反應基因的表達是鹽芥應對中鹽環境的重要因素;而有效清除活性氧、提高滲透物質積累量和增加能量供應可能是鹽芥耐受高鹽環境的關鍵。本研究結果為揭示鹽芥應答鹽脅迫的分子響應機制提供了理論基礎。

　　關鍵詞： 鹽芥;鹽脅迫;差異蛋白質;iTRAQ

　　鹽脅迫是自然界中主要的非生物脅迫之一，土壤中的高濃度鹽分會使植物體內離子失衡，產生滲透壓力，從而嚴重影響植物的生長和發育[1]。鹽芥是一年生草本植物，屬真鹽生植物，多生長于鹽漬化土壤中。鹽芥是擬南芥的近親[2-3]，同樣具有生長周期短、基因組構成小(約為擬南芥的2倍)、種子數目多、可利用農桿菌侵染花序法進行遺傳轉化試驗、遺傳轉化效率高等特點。因此，近幾年鹽芥被提出作為耐鹽分子機理研究的理想模式植物[2，4-5]。雖然鹽芥和擬南芥有很多相似的特性，但兩者之間的耐鹽性卻存在很大差異。通過轉錄組和蛋白質組學技術分析擬南芥和鹽芥在響應鹽脅迫的差異時，發現擬南芥中存在的鹽脅迫響應基因，不管是否處于鹽脅迫環境，在鹽芥中均能高豐度表達[6-7]，說明鹽芥在面對逆境時是“有備無患”;而且鹽芥中存在特有的新陳代謝途徑，使得其在鹽脅迫條件下能減少由植物激素所誘導的生理修飾[8]，從而能在一定程度上承受更大的脅迫壓力。同時研究發現，鹽脅迫條件下，鹽芥根系的磷酸化蛋白質組發生改變，這些變化的磷酸化蛋白質參與了信號轉導、活性氧清除、能量途徑、蛋白質合成以及蛋白質折疊等過程[9]，這些研究為將來從磷酸化水平上研究鹽芥的耐鹽機制提供了重要的理論依據。

　　目前，人們已從鹽芥中鑒定到一些與鹽脅迫耐受相關的基因和調控因子，如編碼K+/Na+轉運蛋白基因(ThHKT1)[10-11]、編碼細胞質膜和液泡膜Na+/H+逆向轉運蛋白基因(ThNHX1)和(ThSOS1)[12-14]、焦磷酸酶基因(TsVP) [15-17]、高親和性K+轉運體(ThHAK5)[18]、脯氨酸合成基因(ThP5CS)[19]、Cu/Zn超氧化物歧化酶基因(ThCSD)[20]等，且部分基因的生物學功能已得到一定的闡述。我們前期對鹽芥葉片和葉綠體響應鹽脅迫的比較蛋白質組學的研究結果表明，鹽芥可通過Na+液泡區隔化，積累滲透調節物質(如淀粉、可溶性糖、脯氨酸等)，以及維持光合效率和生長發育等方式來適應鹽脅迫環境[21-22]。本研究在已有的研究基礎上，應用先進的同位素相對標記與絕對定量技術(iTRAQ)對不同鹽濃度處理條件下的鹽芥葉片進行比較蛋白質組學研究，并重點分析不同鹽濃度處理對鹽芥葉片蛋白質表達豐度的影響，期望揭示鹽芥適應不同鹽濃度條件的蛋白質化學機制，為揭示鹽芥耐鹽分子機制提供參考。

　　1 材料與方法

　　1.1 試驗材料與試驗設計

　　將鹽芥種子直接點播在基質土(營養土∶蛭石=1∶1，體積比)中，放置在培養箱中進行培養。培養條件：溫度為白天22 ℃，晚上20 ℃;相對濕度為60%±5%;16 h光照，8 h黑暗。種子發芽后，用1/2濃度Hoagland營養液每3 d澆灌1次。2個月后，挑選長勢一致的鹽芥植株用含有0 mmol/L、200 mmol/L(低鹽)、400 mmol/L(中鹽)和600 mmol/L(高鹽)NaCl的1/2濃度Hoagland營養液進行澆灌，每天更換新的營養液，每個鹽濃度處理36株植株。收集連續澆灌7 d(根據前期的預試驗確定的鹽處理時間)的鹽芥葉片，每12株為1個生物學重復，液氮速凍，-80 ℃保存備用。

　　1.2 鹽芥葉片蛋白質的提取、酶解及iTRAQ分析

　　參照Wang等[23]的BPP法提取鹽芥葉片總蛋白質，用預冷的甲醇和丙酮清洗蛋白質沉淀，然后溶解在蛋白質裂解液(7 mol/L尿素，2 mol/L硫脲，4% CHAPS(一種非變性的兩性離子型去垢劑)，30 mmol/L Tris-HCl，pH 8.5)中。以牛血清白蛋白(BSA)為標準溶液，通過Bradford測定法測定蛋白質濃度[24]。取100 μg用不同濃度的NaCl處理的鹽芥葉片總蛋白質，用胰蛋白酶進行消化后，使用iTRAQ試劑分別標記0 mmol/L、200 mmol/L、400 mmol/L和 600 mmol/L NaCl處理的上述葉片總蛋白質酶解產物后，再進行等量混合。采用強陽離子變換(SCX)色譜柱(4.6 mm×250.0 mm， Aqua C18， 5 μm， 100 )對iTRAQ標記的肽段進行預分級，再通過液相串聯質譜對標記肽進行鑒定，獲取差異表達肽段信息。

　　1.3 質譜數據搜庫

　　使用Proteinpilot軟件(Version 4.5)對質譜結果進行搜庫，搜索數據庫為自建的鹽芥蛋白質數據庫，蛋白數據從美國國立生物技術信息中心(NCBI)網站(https：//www.ncbi.nlm.nih.gov/genome/12266)上下載。搜庫參數設置為：酶類為Trypsin，錯誤剪切位點數為1，固定修飾位點為Carbamidomethyl(C)，可變修飾位點為Oxidation(M)。將0 mmol/L NaCl處理鹽芥葉片蛋白質作為對照組，200 mmol/L NaCl、400 mmol/L NaCl、600 mmol/L NaCl作為處理組，按處理組與對照組離子的峰面積比值，選擇置信度在95%以上的結果進行報告。質譜鑒定出的可信蛋白質篩選參數應滿足：可信度在95%以上的肽置信水平，至少鑒定出2條肽段，且錯誤發現率(FDR)≤1%。

　　1.4 質譜數據生物信息學分析

　　通過KEGG途徑分析對可信鑒定的蛋白質參與的代謝通路進行富集。同時，對鹽脅迫差異表達蛋白質(變化倍數>1.50或<0.67，P<0.05)進行KEGG代謝途徑和蛋白質間相互作用網絡分析(https：//string-db.org/cgi/input.pl)。

　　1.5 部分鹽脅迫差異表達蛋白質對應基因的qRT-PCR分析

　　提取用不同鹽濃度處理的鹽芥葉片總RNA，取1 μg總RNA反轉錄成cDNA，并將反轉錄后的cDNA樣品稀釋5倍后用于qRT-PCR分析，每個樣品至少重復3次qRT-PCR試驗。取1 μl稀釋后的cDNA加入到SYBR Green PCR Master mix體系中，使用Mx3005P熒光定量PCR儀進行qRT-PCR試驗，鹽芥的actin基因(NCBI基因登錄號312283264)作為內參基因。用于qRT-PCR的引物序列詳見表1。

　　2 結果與分析

　　2.1 鹽芥葉片蛋白質種類分析及鹽脅迫差異表達蛋白質的篩選

　　應用iTRAQ技術結合AB5600+高端生物質譜儀對鹽芥葉片總蛋白質進行定量分析，3次iTRAQ重復鑒定出的可信蛋白質數分別為3 750個、3 897個和3 901個，共計4 607個蛋白質至少在1次重復中鑒定到的可信的鹽芥葉片總蛋白質，其中3 811個蛋白質至少在2次重復試驗中得到可信鑒定(圖1A)。這3 811個可信蛋白質參與了18條代謝途徑，其中碳水化合物代謝途徑參與的蛋白質種類最多(988個蛋白質)，其次是翻譯類蛋白質(895個蛋白質)與折疊、分類和降解類蛋白質(830個蛋白質)、氨基酸代謝(579個蛋白質)、信號轉導(547個蛋白質)、運輸和分解代謝(522個蛋白質)、脂質代謝(496個蛋白質)、環境適應(452種蛋白質)、能量代謝(416種蛋白質)等代謝途徑類蛋白質(圖1B)。

　　推薦閱讀：高中生物課程作業設計與實施的實踐探索