蛋白質組學方法評估轉基因抗蟲玉米非預期效應

郝文媛，李飛武，閆偉，李蔥蔥，郝東云，郭長虹

?

蛋白質組學方法評估轉基因抗蟲玉米非預期效應

郝文媛1,2，李飛武2，閆偉2，李蔥蔥2，郝東云1,2，郭長虹1

（1哈爾濱師范大學生命科學與技術學院，哈爾濱 150025；2吉林省農業(yè)科學院，長春130033）

【】轉基因育種是將對人類有益的外源基因通過生物技術整合進受體植物的基因組中，使其獲得通過自然進化無法獲得優(yōu)異性狀，這是轉入外源基因的“預期效應”。然而，外源基因插入受體植物還可能產生無法控制和預期的細胞、代謝或表型等方面的改變，即“非預期效應”。非預期效應是轉基因植物安全評價的核心內容之一，也是近年來研究的熱點和難點。本研究以抗蟲轉基因玉米為材料，利用差異蛋白質組學技術比較研究轉基因玉米與其非轉基因玉米之間的非預期效應，為蛋白質組學技術解析轉基因植物非預期效應提供理論參考?！尽窟x取4份已經(jīng)進入中國農業(yè)部安全評價階段轉蘇云芽孢桿菌（，Bt）內毒素基因的抗蟲玉米材料，即SK12-5zd、IE034z、Bt799z、Bt799zd，以及它們的對照玉米材料鄭單958和鄭58，分別種植于可控人工溫室中。待玉米幼苗生長至5葉1心期時，每份材料取長勢一致的6個單株最上部完全展開葉片，混合作為一個樣本進行取樣，并提取葉片蛋白質。雙向電泳技術分別分離4種轉基因玉米和它們各自對應的非轉基因材料高豐度蛋白質組，利用Image Scanner掃描儀掃描脫色后的蛋白質凝膠，PD Quest 8.0軟件（Biorad，USA）分析蛋白質凝膠圖譜，以蛋白點相對體積（% Vol）來表述每一個匹配蛋白質的表達豐度，將電泳凝膠上特異性蛋白質和相對體積大于2倍的差異蛋白質點取樣，依據(jù)相應蛋白質數(shù)據(jù)庫進行質譜鑒定，并對鑒定出的特異性和差異性蛋白質進行細胞功能富集分析（GO）和代謝途徑富集分析（KEGG），以評價轉基因抗蟲玉米的可能非預期效應?！尽客ㄟ^比較4種不同轉基因抗蟲玉米及其非轉基因對照的高豐度蛋白質組，除目標抗蟲基因外，共鑒定出61個蛋白質，其功能主要是富集于光合作用、碳固定、能量轉運等基礎細胞功能相關酶類，如核酮糖二磷酸羧化酶、ATP合成酶、丙酮酸磷酸雙激酶等。僅有少數(shù)幾個與光合作用、碳固定和ATP合成途徑等基礎代謝過程相關基因出現(xiàn)上調表達。KEGG分析表明，與對照ZD958相比，SK12-5zd的差異蛋白在光合生物碳固定途徑中顯著富集；Bt799zd的差異蛋白則分別在光合生物碳固定途徑、光合作用和碳代謝途徑中顯著富集。與對照Z58相比，IE034z和Bt799z的差異蛋白均在光合作用代謝途徑中顯著富集。證明4種轉基因材料與對照蛋白質組相似性較高，未見基因的異常表達?！尽繙y試的4種轉基因材料在蛋白質水平未發(fā)生影響其安全性的非預期效應；當前結果表明蛋白質組學技術可以用來驗證轉基因植物的非預期效應。

轉基因玉米；蛋白質組學；非預期效應；雙向電泳；差異表達蛋白

0 引言

【研究意義】植物轉基因技術的目標是將有利于植物自身或滿足人類需求的特定外源基因整合進入受體植物基因組中，使其獲得外源基因所表現(xiàn)的特定性狀并能夠穩(wěn)定遺傳。這些性狀就是轉基因植物的目標性狀，也稱為預期效應。由于外源基因整合位點不確定，加之外源基因產物本身也可能使受體植物的能量代謝和物質代謝發(fā)生變化，從而影響其固有的生理代謝過程，乃至改變受體植物的固有性狀。這些性狀改變均是轉基因育種設計過程中無法控制和預期的，這些無法控制和預期的性狀改變均稱為非預期效應[1-2]。可見，轉基因植物實際上是外源基因發(fā)揮預期的目標效應和受體植物基因組因外源基因的導入，而產生的非預期效應的綜合體。轉基因植物是否會產生非預期效應，以及對非預期效應的安全性評估一直是科學家和公眾對轉基因植物廣泛關注的話題。長期以來，各國相關部門都是基于1993年歐洲經(jīng)濟發(fā)展與合作組織OECD（Organization for Economic Co-operation and Development）提出的“實質等同性原則”對轉基因植物進行安全性評價的，評價的核心是比較轉基因植物與其非轉基因受體植物間在農藝性狀、環(huán)境適應性、營養(yǎng)及物質組成變化等的差異，以確認或排除轉基因植物是否存在非預期效應。傳統(tǒng)的檢測方法不僅難以全面涵蓋轉基因植物可能出現(xiàn)的非預期效應，且檢測指標有限、檢測費用高、費時、費力。因此，建立全面、快速、準確、靈敏的檢測技術是分析鑒定轉基因生物安全性的關鍵所在[2-3]?！厩叭搜芯窟M展】近年來，蛋白質組分析技術迅速發(fā)展，雙向電泳技術、iTRAQ技術等已成為當今蛋白質組學研究的核心技術[4]。研究人員已經(jīng)在轉基因玉米、水稻、大豆等材料中，對蛋白質組水平上的非預期效應的進行了評估，取得了較好的效果[5-11]。轉基因玉米商業(yè)化生產以來，已占全球轉基因作物種植總面積的30%，居轉基因作物的第二位，涉及的農藝性狀主要為玉米螟抗性，轉入的基因為種類繁多的來自蘇云芽孢桿菌（，Bt）的昆蟲內毒素基因[12-15]（本文統(tǒng)稱“Bt基因”）。Coll等[5]利用2D凝膠電泳、質譜技術分析比較，發(fā)現(xiàn)在同一地塊種植的2種轉Bt基因玉米MON810與其對照株系的籽粒蛋白質組幾乎一致，只有幾個蛋白質點有1—1.8倍的變化，結果表明，轉基因玉米MON810與其對照非轉基因材料具有實質等同性；Albo等[6]也發(fā)現(xiàn)MON810與其非轉基因對照的籽粒蛋白質組之間非預期差異基本等同；Zolla等[7]發(fā)現(xiàn)同一田間生長的MON810及其對照同一品系的不同世代材料的蛋白質組存在顯著差異，且在蛋白質組學水平上，環(huán)境因素對非預期變異效應的貢獻比轉基因插入更大。上述研究結果表明，蛋白質組學技術已逐漸成為檢測轉基因作物非預期效應的有效方法[16]。【本研究切入點】不同轉Bt基因的抗蟲玉米材料不但在目標基因上明顯不同，而且由于受體材料差異和回交轉育過程的不同，造成了最終的抗蟲轉基因玉米材料間存在著巨大的差異。這些遺傳背景上的差異是否會產生非預期效應仍研究甚少，目前的研究還主要集中在Bt基因的抗蟲目標預期效應上[17-22]。轉基因抗蟲玉米SK12-5、IE034、Bt799z和 Bt799zd是中國自主研發(fā)的具有潛在產業(yè)化前景的轉基因新材料，圍繞上述新材料已開展了大量育種研發(fā)工作及所獲相應育種材料的環(huán)境安全[23-25]和食用安全評價研究，但是，對這些材料的非預期效應研究未見報道。【擬解決的關鍵問題】利用雙向電泳技術對上述轉基因抗蟲玉米的4種不同遺傳背景材料的轉基因材料進行了比較分析，一方面在蛋白質組水平比較不同轉基因材料的差異；另一方面從組學的水平初步評估了實驗材料蛋白水平上是否發(fā)生了影響轉基因材料安全性的非預期效應，為非預期效應評價提供數(shù)據(jù)支撐。

1 材料與方法

1.1 材料

供試材料：編號1#為具有轉并有ZD958遺傳背景的雜交種SK12-5zd，由浙江大學提供。編號2#為轉并有Z58遺傳背景的自交系IE034z，由中國農業(yè)科學院提供。編號3#為轉并有Z58遺傳背景的自交系Bt799z，編號4#為轉并有ZD958遺傳背景的雜交種Bt799zd，均由中國農業(yè)大學提供。選取非轉基因的自交系Z58和雜交種ZD958作為對照材料，編號分別為5#和6#。編號2#和3#的轉基因自交系是目標遺傳轉化體經(jīng)2代自交后，再與Z58材料進行6代回交轉育后獲得了具有純合、穩(wěn)定遺傳背景的轉基因自交系。上述轉基因材料均已采用國家發(fā)布的轉基因檢測標準進行了外源基因分子特征鑒定，證明均為獨立轉化事件，并已進入環(huán)境釋放以上的轉基因安全評價試驗階段。

2014年6月，將試驗材料種植在吉林省公主嶺市的吉林省農業(yè)科學院的人工溫室中，幼苗長至5葉1心期進行取樣，每種試驗材料取長勢基本一致的6個單株的最上部完全展開葉片，混合作為一個樣本。液氮冷凍后，迅速放入干冰中運輸至實驗室，-80℃保存。

1.2 玉米葉片總蛋白提取

將葉片去主脈后剪成1 cm的小段，稱取0.5 g葉片加入50 mg PVPP進行充分研磨；加入7.5 mL經(jīng)-20℃預冷的蛋白提取緩沖液Ⅰ，冰浴研磨5 min，轉入預冷的10 mL離心管中，-20℃放置1 h；4℃，13 000 r/min離心45 min，棄上清；加入7.5 mL經(jīng)-20℃預冷的蛋白質提取液Ⅱ，懸浮振蕩，-20℃放置1 h；4℃，13 000 r/min離心45 min，棄上清，重復2次；加入7.5 mL 80%預冷的丙酮，充分懸浮后，離心去上清；沉淀真空干燥并稱重。每種材料各重復3次，獲得的蛋白用Bradford法[26]測定濃度后準備進行蛋白質雙向電泳分析。

1.3 蛋白質等點聚焦

采用長度24 cm，pH 3—10的線性干膠條（Biorad，美國），取蛋白提取液與水化緩沖液（8 mol·L-1Urea、2 mol·L-1硫脲、4% CHAPS、0.002%溴酚藍，用前加入65 mmol·L-1DTT和0.2% IPG Buffer）混合至總體積450 μL（含蛋白質約為1 200 μg），采用被動水化的方法處理18 h。將水化后的膠條轉移到等電聚焦槽上（膠條酸性端置聚焦盤的“+”端，堿性端置聚焦盤“-”端）。每根膠條加入3 mL礦物油，將潤濕的濾紙1/3處放在膠條的兩端夾上電極夾，扣緊聚焦盤，等點聚焦程序和參數(shù)設置參照Biorad產品說明書。等點聚焦結束后，將膠條置于平衡緩沖液（6 mol·L-1Urea、75 mmol·L-1Tris-HCl（pH 8.8）、2% SDS、25%甘油、0.002%溴酚藍）振蕩平衡2次，每次15 min。第一次平衡在緩沖液中加入1% DDT；第二次平衡在緩沖液加入2.5%碘乙酰胺（IAM）。平衡完畢后取出膠條輕輕用平衡緩沖液潤洗，然后去除多余的平衡緩沖液，準備SDS-PAGE電泳。

1.4 SDS-PAGE

將平衡好的膠條置于第二向凝膠上方，排除氣泡使二者緊密接觸，用1%瓊脂糖凝膠封固膠條，接通電源。第一階段1 W/gel運行1 h；第二階段13 W/gel運行直至溴酚藍到達凝膠底端。雙向電泳結束后，SDS-PAGE凝膠置于固定液（10%冰醋酸、40%甲醇、50%超純水）中固定3 h；凝膠固定后，超純水沖洗3次，每次15 min；加入考馬斯亮藍染色液染色（0.12% G-250、10%硫酸銨、10%磷酸、20%甲醇）振蕩過夜；超純水脫色至背景滿意為止。

1.5 質譜鑒定與數(shù)據(jù)庫分析

Image Scanner掃描儀掃描脫色后的凝膠，用Biorad PDQuest 8.0軟件分析不同玉米葉片材料的凝膠圖譜，采用Yang等[4]的方法，使用蛋白相對表達體積（% Vol）來表述每一個匹配蛋白質點的表達豐度。差異蛋白質的表達豐度分析分別以Z58和ZD958材料的凝膠作為參考膠，只有蛋白質點在至少2個重復內被發(fā)現(xiàn)才被認為是真實的蛋白質點。定義平均相對表達豐度上調（出現(xiàn)）或下調（消失）差異大于2倍（ratio＞2），并且在統(tǒng)計學T檢驗上＜0.05的蛋白質為差異蛋白質點。差異蛋白質點使用液質聯(lián)用質譜儀LC-Q-TOF-MS 6520（Agilent Technologies，USA）進行分析。使用NCBI玉米數(shù)據(jù)庫（NCBI，Taxonomy：，Entry：212235，2015/01/20）對蛋白質進行檢索。蛋白檢索結果的顯著性判別標準：一般認為PMF（Peptide Mass Fingerprinting，肽指紋圖譜）在MS/MS檢索中離子分數(shù)≥88分說明該蛋白的鑒定結果有效，其他一些輔助標準包括至少10%的序列覆蓋率以及需要有3個肽段的匹配[27]。

1.6 GO和KEGG分析

采用BLAST（https://blast.ncbi.nlm.nih.gov）方法對質譜鑒定出的特異性蛋白（E value≤10?5）與NCBI non-redundant（Nr）蛋白質數(shù)據(jù)庫進行比較，確定相應的蛋白質并獲得對應的基因信息，用DAVID方法對這些表達有差異的基因進行基因功能注釋（gene ontology，GO）分析[28]，對差異蛋白功能及代謝途徑進行KEGG通路分析[29]。

2 結果

2.1 差異表達蛋白的分離及鑒定

為了比較不同轉基因抗蟲玉米蛋白質組差異，對供試材料的葉片總蛋白進行了雙向電泳分析。結果顯示，各受試材料均可檢測到超過1 200個清晰蛋白點（圖1）。

將電泳膠圖掃描后使用Biorad PDQuest8.0軟件對每個轉基因材料及其對照（1#和6#、4#和6#、2#和5#、3#和5#）進行比較，對電泳凝膠上相對表達體積超過2倍的差異蛋白質點進行質譜儀分析，成功鑒定玉米基因組蛋白61個（圖2、表1）。進一步分析特異蛋白點（只在轉基因與對照的一個材料中出現(xiàn)）和差異蛋白點發(fā)現(xiàn)，這些蛋白點中核酮糖二磷酸羧化酶相關蛋白點15個，ATP合酶類蛋白點4個，放氧增強蛋白點13個，丙酮酸磷酸雙激酶相關蛋白點3個，GADPH酶相關蛋白點2個，類萌發(fā)素蛋白相關蛋白點4個，草銨膦乙酰轉移酶相關蛋白點2個（3#和4#），其他蛋白點6個，未知功能蛋白點10個。這些蛋白所涉及的生化反應過程主要集中在光合作用、碳固定和能量代謝等方面。根據(jù)不同蛋白的功能以及質譜檢測結果，對各類蛋白在不同轉基因材料中變化的趨勢進行了匯總分析（表2）。結果顯示，在4種轉基因玉米材料中發(fā)現(xiàn)被認為參與植物對逆境脅迫響應的類萌發(fā)素類蛋白均有顯著的表達升高；而只有質體藍素鈉一類蛋白在4#材料中發(fā)生了表達顯著降低；涉及光合作用及碳固定的相關酶類，如核酮糖二磷酸羧化酶、ATP合成酶、丙酮酸磷酸雙激酶等均發(fā)生了表達顯著上調。在同為ZD958背景的1#和4#材料中，包括類萌發(fā)素蛋白、核酮糖二磷酸羧化酶、PSⅡ的放氧增強蛋白的表達均發(fā)生了顯著的升高；而在同為Z58背景的2#和3#材料中，包括ATP合成酶、放氧增強蛋白、類萌發(fā)素蛋白、肽基脯氨酰順反式異構酶等蛋白的表達發(fā)生了顯著升高?？梢钥闯鲞z傳背景相同的轉基因玉米材料（1#和4#、2#和3#）在轉入不同抗蟲Bt基因后，其所引起的蛋白水平上的變化也存在著一定的一致性。而不同遺傳背景的轉基因材料（3#和4#）轉入相同抗蟲Bt基因后，其蛋白水平上的表達變化則差異較大。

圖2 轉基因玉米及其對照葉片蛋白質的雙向電泳對比分析

表1 轉基因玉米葉片與對照的差異表達蛋白

此外，對4份轉基因材料獲得的特異蛋白點與轉基因材料中所轉入的Bt蛋白序列進行比對，發(fā)現(xiàn)了特異蛋白點#3802、#4719、#3707、#3708分別對應Bt蛋白cry1Ab、cry1Ie、mcry1Ac、mcry1Ac（圖1）。說明轉基因材料中轉入的Bt基因在玉米中能夠正常表達，并且這些Bt蛋白能夠通過雙向電泳被分離和鑒定，最終體現(xiàn)在蛋白質組分析所獲得的特異蛋白點中。

2.2 差異蛋白的GO分析和KEGG分析

雖然通過BLASTx對試驗所得的差異蛋白點進行了初步鑒定，獲得了蛋白信息，但這些蛋白在是否在細胞組成、生物途徑和分子功能上具有共性，是否富集于某一代謝途徑，仍無法準確推斷。通過對這些差異蛋白質的生物學功能和代謝通路富集分別進行GO（圖3）和KEGG分析（圖4）。發(fā)現(xiàn)轉基因材料與對照材料相比雖有差異，但也主要集中在光合作用相關蛋白，并未發(fā)現(xiàn)特殊蛋白或特殊代謝途徑的明顯改變。SK12-5zd與ZD958比較發(fā)現(xiàn)，共有5個不同生物途徑的相關蛋白發(fā)生了差異表達，分別為光合作用（photosynthesis）相關蛋白3個、還原型戊糖磷酸循環(huán)途徑（reductive pentose-phosphate cycle）相關蛋白2個、鈣離子結合（calcium ion binding）蛋白2個、不同細胞組成相關蛋白4個（具體為膜外蛋白（extrinsic component of membrane）2個和光合系統(tǒng)Ⅱ放氧復合體（photosystem Ⅱ oxygen evolving complex）相關蛋白2個。

Bt799zd與ZD958比較發(fā)現(xiàn)，有7個不同生物途徑相關蛋白發(fā)生了差異表達，分別為光合作用相關蛋白4個，光呼吸（photorespiration）相關蛋白質2個和碳固定途徑（carbon fixation）相關蛋白2個；發(fā)現(xiàn)4個不同分子功能的蛋白發(fā)生了差異表達，分別為2個核酮糖-二磷酸羧化酶的活性蛋白（ribulose- bisphosphate carboxylase activity）和2個單加氧酶活性蛋白（monooxygenase activity）；按照細胞組成分析也發(fā)現(xiàn)了4個蛋白發(fā)生了差異表達，分別為膜外蛋白（extrinsic component of membrane）2個和光合系統(tǒng)Ⅱ放氧復合體（photosystem II oxygen evolving complex）相關蛋白2個。

IE034z與Z58比較發(fā)現(xiàn)，有4個不同生物途徑相關蛋白發(fā)生了差異表達，均屬于光合途徑相關蛋白；有5個細胞組成相關蛋白發(fā)生了差異表達，分別為膜外蛋白3個和光合系統(tǒng)Ⅱ放氧復合體相關蛋白2個；按照分子功能分析則發(fā)現(xiàn)了鈣離子結合的蛋白發(fā)生了差異表達3個。

表2 轉基因玉米差異表達蛋白匯總

Bt799z與Z58比較發(fā)現(xiàn)，有4個生物途徑相關蛋白發(fā)生了差異表達，均屬于光合途徑相關蛋白；按照細胞組成分析有6個蛋白發(fā)生了差異表達，分別為膜外蛋白3個和光合系統(tǒng)Ⅱ放氧復合體相關蛋白3個；還發(fā)現(xiàn)有3個分子功能為鈣離子結合的蛋白發(fā)生了差異表達。

通過2倍上下調控表達水平來篩選轉基因材料與其非轉基因對照材料的差異蛋白質進行KEGG分析則發(fā)現(xiàn)，在SK12-5zd與ZD958的比較中，差異蛋白在光合生物碳固定途徑（carbon fixation in photosynthetic organisms）中顯著富集；而在Bt799zd和ZD958的比較中，差異蛋白則分別在光合生物碳固定途徑、光合作用和碳代謝（carbon metabolism）幾個代謝途徑中顯著富集。在IE034z和Z58以及Bt799z和Z58的比較中均為光合作用代謝途徑中的差異蛋白顯著富集（圖4）。

圖3 轉基因玉米葉片與其對照的差異蛋白GO分析

圖4 轉基因玉米葉片與其對照的差異蛋白KEGG分析

3 討論

隨著轉基因植物的研發(fā)和推廣，轉基因生物及其產品的分析鑒定，以及轉入基因的直接表達產物和間接調控產物對轉基因生物生長發(fā)育、生態(tài)環(huán)境、食物鏈及人類食品安全的影響越來越受到關注。由于現(xiàn)有轉基因生物中多數(shù)外源基因的表達或調控產物都是蛋白質，因而在蛋白質水平上對轉基因生物進行評價顯得尤為重要，蛋白質組學的興起和發(fā)展為在蛋白質水平上對轉基因生物進行評價提供了可能[30-31]。

本研究對4種轉基因抗蟲玉米材料及其非轉基因對照玉米材料進行了蛋白質組分析。每份材料經(jīng)雙向電泳凝膠分離后均可檢測到1 200個蛋白質點，利用PDQuest8.0軟件對蛋白質組表達譜進行分析，已檢測到4份轉基因材料與對照之間差異表達的Bt殺蟲蛋白，還鑒定出61個表達量發(fā)生顯著變化的蛋白質點，主要集中于光合作用、同化物的積累及能量代謝過程。不同轉基因材料中均檢測到光合生物碳固定途徑中的差異蛋白，如光合作用中決定碳同化速率的關鍵酶——核酮糖1,5-二磷酸羧化酶、C4植物光合作用途徑中二氧化碳固定的關鍵酶——丙酮酸磷酸雙激酶、放氧增強蛋白；植物中對逆境響應的類萌發(fā)素蛋白、以及ATP合成酶類等在4個轉基因材料中均有顯著的上調表達，這些差異蛋白不僅是遺傳背景的差異，可能也與轉基因玉米具有抗蟲性相關，而GO分析和KEGG分析的結果也證實了上述差異的存在。理論上光合作用的增強，生物體抗逆能力增強，生長會更健壯，植物的光合能力也會提高，這有待于后期對材料的觀察和分析；同時光合作用和碳固定過程增強，均涉及大量的能量物質流動，這些能量物質的去向與外源蛋白的表達是否存在關聯(lián)，有待深入分析。

本研究發(fā)現(xiàn)類萌發(fā)素蛋白在所有轉基因材料中均有顯著的上調表達，該蛋白是在植物中普遍存在的位于胞外基質的可溶性糖蛋白，絕大多數(shù)為穩(wěn)定的低聚物，其作用為參與植物對逆境脅迫的響應。這類蛋白的上調表達可能預示著抗蟲轉基因玉米材料中，除依賴表達針對靶標生物的毒蛋白已達到抗蟲的目的外，還可能存在某種共有的變化機制，對其他逆境脅迫響應產生作用，而具體作用機制需通過后續(xù)的逆境脅迫響應評估及分析深入研究。

轉基因作物的蛋白質組差異一方面是由于外源基因的插入造成的，另一方面還有可能由環(huán)境因素造成[7]。本研究材料均在溫室內培養(yǎng)至5葉1心期，極大程度地避免了環(huán)境因素的影響，其結果所反應出的差異應主要為外源基因插入造成的蛋白質組差異。在此基礎上，發(fā)現(xiàn)受體材料的遺傳背景可能影響轉基因材料的蛋白表達差異。本研究中所用的轉基因自交系材料均是經(jīng)過高世代選育獲得的穩(wěn)定材料，其基因組差異極小，因此，在與其回交親本材料比較中所發(fā)現(xiàn)的蛋白質差異，應主要來源于外源基因的插入及表達所引起的蛋白質組變化。研究結果也證明在相同遺傳背景的轉基因材料中，其蛋白表達變化的一致性較高，而經(jīng)過雜交配制的轉基因雜交種材料，由于雜交過程中的不確定性，外源基因的表達和雜交過程均會對蛋白質組產生影響，因此導致檢測出了更多的差異表達蛋白，本研究結果也再次證實了不同遺傳背景但轉入相同基因的材料間（3#和4#）其蛋白水平上的表達差異較大。這與抗蟲玉米Mon810的相關蛋白質組的比較分析結果基本一致，即在轉基因玉米個體的蛋白質表達譜之間沒有顯著的個體差異，其主要差異源于基因組遺傳背景的差異[5-7]。

抗蟲轉基因玉米生物安全評價過程中非常重視其食用和飼用安全性。綜合分析本研究結果，在4份轉基因玉米材料中并沒有發(fā)現(xiàn)影響其食用或飼用安全的新蛋白，未發(fā)現(xiàn)植物生理、生化和代謝過程中發(fā)揮重要功能的蛋白質或致敏和毒素蛋白質產生或發(fā)生明顯的表達變化。由于蛋白質是具有生物功能的特殊大分子物質，轉基因植物的食用安全性變化會直接體現(xiàn)在蛋白質水平上。本研究表明抗蟲轉基因玉米在蛋白質水平上并未產生可檢出的影響其安全性的蛋白成分。轉基因材料中發(fā)生變化的蛋白多集中于光合作用、同化物的積累及能量代謝過程，這些過程也與轉基因材料安全性關系小，并非影響轉基因材料安全性的非預期效應。

4 結論

SK12-5zd、IE034z、Bt799z、Bt799zd等4種轉基因抗蟲玉米材料與對照蛋白質組相似性較高，僅有少數(shù)幾個與光合作用、碳固定和ATP合成途徑等基礎代謝過程相關基因出現(xiàn)上調表達，但未見基因的異常表達。測試的4種轉基因材料在蛋白質水平未發(fā)生影響其安全性的非預期效應。

[1] Cellini F, Chesson A, Colquhoun I, Constable A, Davies H V, Engel K H, Gatehouse A M, K?renlampi S, Kok E J, Leguay J J. Unintended effects and their detection in genetically modified crops., 2004, 42(7): 1089-1125.

[2] 趙艷, 李燕燕. 組學技術評價轉基因農作物的非預期效應. 遺傳, 2013, 35(12): 1360-1367.

Zhao Y, Li Y Y. Unintended effects assessment of genetically modified crops using omics techniques., 2013, 35(12): 1360-1367. (in Chinese)

[3] Arthur J W, Wilkins M R. Using proteomics to mine genome sequences., 2016, 3(3): 393-402.

[4] Yang P, Li F, X J, Wang X, Chen H, Chen F, Shen S H. Proteomic analysis of rice () seeds during germination., 2007, 7(18): 3358-3368.

[5] Coll A, Nadal A, Rossignol M, Puigdomènech P, Pla M. Proteomic analysis of MON810 and comparable non-GM maize varieties grown in agricultural fields., 2011, 20(4): 939-949.

[6] Albo A G, Mila S, Digilio A G, Motto M, Aime S, Corpillo D. Proteomic analysis of a genetically modified maize flour carrying Cry1Ab gene and comparison to the corresponding wild-type., 2007, 52(4): 443-455.

[7] Zolla L, Rinalducci S, Antonioli P, Righetti P G. Proteomics as a complementary tool for identifying unintended side effects occurring in transgenic maize seeds as a result of genetic modifications., 2008, 7(5): 1850-1861.

[8] Wang Y, Xu W, Zhao W, Hao J, Luo Y, Tang X, Zhang Y, Huang K. Comparative analysis of the proteomic and nutritional composition of transgenic rice seeds with Cry1ab/ac genes and their non-transgenic counterparts., 2012, 55(2): 226-233.

[9] Xue K, Yang J, Liu B, Xue D. The integrated risk assessment of transgenic rice: A comparative proteomics approach., 2012, 135(1): 314-318.

[10] Barbosa H S, Arruda S C, Azevedo R A, Arruda M A. New insights on proteomics of transgenic soybean seeds: evaluation of differential expressions of enzymes and proteins., 2012, 402(1): 299-314.

[11] Storer N P, Thompson G D, Head G P, Hutchison W, Gassmann A. Application of pyramided traits against Lepidoptera in insect resistance management for Bt crops., 2012, 3(3): 154-162.

[12] Xu L, Wang Z, Zhang J, He K, Ferry N, Gatehouse A M R. Cross-resistance of Cry1Ab-selected Asian maize borer to other Cry toxins., 2010, 134(5): 429-438.

[13] Li H, Olson M, Lin G, Hey T, Tan S Y, Narva K E.Bacillus thuringiensis Cry34Ab1/Cry35Ab1 interactions with western maize rootworm midgut membrane binding sites., 2013, 8(1): e53079.

[14] 周宗良, 林智敏, 耿麗麗, 蘇軍, 束長龍, 王鋒, 朱彥明, 張杰. 水稻中cry1Ah1基因密碼子優(yōu)化方案的比較. 生物工程學報, 2012, 28(10): 1184-1194.

Zhou Z L, Lin Z M, Geng L L, Su J, Shu C L, Wang F, Zhu Y M, Zhang J. Comparison of codon optimizations of cry1Ah1 gene in rice., 2012, 28(10): 1184-1194. (in Chinese)

[15] Guo Y Y, Tian J C, Shi W P, Dong X H, Romeis J, Naranjo S E, Hellmich R L, Shelton A M. The interaction of two-spotted spider mites, Tetranychus urticae Koch, with Cry protein production and predation by Amblyseius andersoni (Chant) in Cry1Ac/Cry2Ab cotton and Cry1F maize., 2016, 25(1): 33-44.

[16] 李欣, 黃昆侖, 朱本忠, 唐茂芝, 羅云波.利用“組學”技術檢測轉基因作物非期望效應的潛在性. 農業(yè)生物技術學報, 2005, 13(6): 802-807.

Li X, Huang K l, Zhu B z, Tang M z, Luo Y b. Potentiality of “Omics” techniques for the detection of unintended effectsin genetically modified crops., 2005, 13(6): 802-807. (in Chinese)

[17] Mathur C, Kathuria P C, Dahiya P, Singh A B. Lack of detectable allergenicity in genetically modified maize containing “Cry” proteins as compared to native maize based on in silico & in vitro analysis., 2015, 10(2): e0117340.

[18] Agapito-Tenfen S Z, Vilperte V, Benevenuto R F, Rover C M, Traavik T I, Nodari R O.Effect of stacking insecticidal cry and herbicide tolerance epsps transgenes on transgenic maize proteome., 2014, 14: 346.

[19] Hoss S, Menzel R, Gessler F, Nguyen H T, Jehle J A, Traunspurger W. Effects of insecticidal crystal proteins (Cry proteins) produced by genetically modified maize (Bt maize) on the nematode., 2013, 178: 147-151.

[20] Hendriksma H P, Kuting M, Hartel S, Nather A, Dohrmann A B, Steffan-Dewenter I, Tebbe C C. Effect of stacked insecticidal Cry proteins from maize pollen on nurse bees () and their gut bacteria., 2013, 8(3): e59589.

[21] Zhang L, Huang F, Rogers Leonard B, Chen M, Clark T, Zhu Y C, Wangila D S, Yang F, Niu Y. Susceptibility of Cry1Ab maize-resistant and susceptible strains of sugarcane borer (Lepidoptera: Crambidae) to four individual Cry proteins., 2013, 112(3): 267-272.

[22] Li Y, Meissle M, Romeis J. Use of maize pollen by adult(Neuroptera: Chrysopidae) and fate of Cry proteins in Bt-transgenic varieties., 2010, 56(2): 157-164.

[23] 王江, 武奉慈, 劉新穎, 馮樹丹, 宋新元. 轉基因抗蟲耐除草劑復合性狀玉米‘雙抗12-5’對亞洲玉米螟的抗性及對草甘膦的耐受性研究. 植物保護, 2016, 42(1): 45-50.

Wang J, Wu F C, Liu X Y, Feng S D, Song X Y. Evaluation of transgenic maize ‘Shuangkang 12-5’ with complex traits of insect- resistance and glyphosate-resistance for the resistance toand tolerance to glyphosate., 2016, 42(1): 45-50. (in Chinese)

[24] 劉新穎, 王柏鳳, 周琳, 馮樹丹, 宋新元. 轉基因抗蟲玉米IE09S034種植對田間大型土壤動物多樣性的影響. 作物雜志, 2016(1): 62-68.

Liu X Y, Wang B F, Zhou L, Feng S D, Song X Y. Effects ofon soil macrofauna diversity in transgenic corn IE09S034 fields., 2016(1): 62-68. (in Chinese)

[25] 劉新穎, 王柏鳳, 王江, 馮樹丹, 宋新元. 轉基因抗蟲玉米葉片殘體降解對土壤動物群落結構的影響. 植物保護學報, 2016, 43(3): 384-390.

Liu X Y, Wang B F, Wang J, Feng S D, Song X Y. Effects of leaf residue decomposition of-transgenic insect-resistant maize on community structure of soil animals., 2016, 43(3): 384-390. (in Chinese)

[26] Zor T, Selinger Z. Linearization of the Bradford protein assay increases its sensitivity: theoretical and experimental studies., 1996, 236: 302-308.

[27] 于濤, 李耕, 劉鵬, 董樹亭, 張吉旺, 趙斌, 柏晗.玉米早期發(fā)育階段粒位效應的蛋白質組學分析. 中國農業(yè)科學, 2016, 49(1): 54-68.

Yu T, Li G, Liu P, Dong S T, Zhang J W, Zhao B, Bai H.Proteomics analysis of grain position effects during early developmental stages of maize.,2016, 49(1): 54-68. (in Chinese)

[28] Huang D W, Sherman B T, Lempicki R A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources., 2009, 4(1): 44-57.

[29] Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome., 2004, 32(Database issue): D277-D280.

[30] Pan Y H. Application prospects of proteomicsin detection and study of genetically modified organisms., 2010, 12(1): 31-34.

[31] Cho J I, Lee G S, Park S C.Global status of GM crop development and commercialization., 2016, 43(2): 147-150.

（責任編輯 李莉，岳梅）

assessment of the unintended effects of four genetically modified Maize varieties by Proteomic approach

Hao WenYuan1,2, Li FeiWu2, Yan Wei2, Li CongCong2, Hao DongYun1,2, Guo ChangHong1

(1College of Life Sciences and Technology, Harbin Normal University, Harbin 150025;2)

【】Unintended effects of genetically modified (GM) crops are defined as the unexpected trait changes in GM breeding and research. Over the last decade, much attention has been taken to GM crops’ intended effects that are attributed to the insertion of foreign gene, leaving the unintended effects of GM crops less investigated. In this study, a comparative proteomic study using two-dimensional electrophoresis and mass spectrometry was conducted to assess the high-abundant protein variation, so as to explore possible occurrence of their unintended effects.【】Four insect resistant GM maize varieties, i.e SK12-5zd, IE034z, Bt799z and Bt799zd, which were under approval of GM safety assessment by the Ministry of Agriculture, were planted in a condition-controlled greenhouse together with their corresponding CKs, i.e Zheng58 and Zhengdan958. When seedling growing into 5-leaf stage, the top-leaf blades of each maize variety were sampled and their leaf proteins were extracted. After profiling by two-dimensional electrophoresis, on-gel leaf protein spots of each GM maize were compared against their corresponding CKs using Imager Scanner and PD Quest 8.0 (Biored, USA). The resulted unique and differential protein (with a relative volume greater than 2-fold) spots were isolated from the gel and applied for protein identification by mass spectrometry against the relevant protein database. The identified proteins were analyzed in gene ontology (GO) and KEGG for their cellular functionalities and metabolism pathway enrichment. 【】Through comparison of the unique and differential proteins between each GM maize and its corresponding no-GM CKs among 4 insect resistant GM maize varieties, a total of 61 differential proteins were identified in maize PDB (NCBI) alongside the 4 intended insect resistant gene products. The gene ontology analyses of the identified proteins revealed that they are popular cellular enzymes of fundamental metabolism, i.e Rubisco, ATP synthase, pyruvate orthophosphate dikinase, etc. Only a few genes relating to photosynthesis, ATP synthesis and carbon fixation were observed to be up-expressed, possibly due to the genetic background difference between the GM maize and its CK. KEGG analyses uncovered that, in comparison with the corresponding CKs, the differential proteins of SK12-5zd and Bt799zd were enriched similarly into photosynthesis, carbon fixation and metabolism pathways. In comparison with the CK Z58, the differential proteins of both IE034z and Bt799z were enriched into the photosynthetic metabolism pathway.Four GM maize varieties exhibited extensive similarity over their corresponding wild-type controls. It appeared to be no unintended effects observed. 【】The preliminary results revealed no obvious unintended effects in 4 insect resistant GM maize varieties in respect to high-abundant proteomic assessment.

genetically modified maize; proteomics; unintended effects; two-dimensional electrophoresis (2-DE); differentially expressed protein

10.3864/j.issn.0578-1752.2017.19.003

2017-03-20；接受日期：2017-05-03

國家轉基因生物新品種培育科技重大專項（2015ZX08013003）

郝文媛，E-mail：wenyuan_h@163.com。通信作者郭長虹，E-mail：kaku2008@hotmail.com。通信作者郝東云，E-mail：dyhao@cjaas.com