基于iTRAQ技术对棉花叶片响应化学打顶的差异蛋白质组学分析

侯晓梦，刘连涛，李梦，孙红春，张永江，杜欢，李存东



基于iTRAQ技术对棉花叶片响应化学打顶的差异蛋白质组学分析

侯晓梦1,2，刘连涛1,2，李梦4，孙红春1,2，张永江1,2，杜欢2,3，李存东2

（1河北农业大学农学院，河北保定 071000；2河北省作物生长调控实验室，河北保定 071000；3河北农业大学生命科学学院，河北保定 071000；4河北省农林科学院，石家庄 050051）

【】比较化学打顶和人工打顶方式下棉花植株生理变化及蛋白质的差异表达，为化学打顶的作用机理提供理论依据。【】以黄河流域大面积推广的冀棉863为试验品种，于2015—2016年设置人工打顶、化学打顶和未打顶3种处理，于7月20日进行统一打顶处理，化学打顶剂为人工喷施，用量为1.125 L·hm-2。打顶处理后定期测定各处理间棉花株高与主茎功能叶内源激素含量。株高测量为子叶节到主茎生长点的高度，使用直尺测量。使用酶联免疫法测定棉花功能叶的生长素（IAA）、赤霉素（GA3）、脱落酸（ABA）和玉米素核苷（ZR）含量。采用iTRAQ技术对人工打顶和化学打顶处理的打顶后20 d的主茎功能叶进行差异蛋白质组学分析。【】与人工打顶的棉花相比，化学打顶处理株高显著高于人工打顶处理，两年试验中分别高11.8%和14.5%，但显著低于未打顶处理，两年试验中分别低6.0%和6.5%，喷施化学打顶剂有效抑制了棉花株高的增长。不同打顶处理对棉花功能叶GA3含量影响较大，打顶后GA3含量变化为单峰曲线，处理30 d各处理之间达到显著差异，GA3含量为未打顶＞化学打顶＞人工打顶，30 d后化学打顶与未打顶处理呈下降趋势，人工打顶处理则在20 d时出现下降趋势，在处理后50 d时各处理GA3含量无显著差异。2016年IAA含量峰值出现在处理后40 d，化学打顶处理峰值显著低于其他两个处理，2015年3种打顶处理间无显著差异。ABA含量在处理后40 d时达到最大值，未打顶处理峰值显著低于其他两个处理。3种打顶处理的ZR含量无显著差异。化学打顶与人工打顶处理相比，iTRAQ标记方法检测到69个差异表达蛋白，29个上调表达，40个下调表达，其中碳水化合物和能量代谢相关的蛋白多下调表达，降低了植株的长势；与GA调节正相关蛋白多上调表达，增强GA效应。【】化学打顶能有效控制棉花株高，对棉花功能叶的GA含量影响较大，化学打顶处理含量显著高于人工打顶处理，与人工打顶处理相比，化学打顶与植物生长发育相关蛋白多下调表达，可能是植株通过降低碳水化合物合成，减少能量代谢，增加GA含量，激活GA效应来实现株高的控制。

棉花；化学打顶；激素；iTRAQ技术；差异蛋白质

0 引言

【研究意义】棉花是中国重要的经济作物，随着劳动力紧缺以及成本增加，植棉的机械化、规模化成为现代植棉业的必然趋势[1]。目前，人工打顶制约了棉花生产轻简化、规模化、精准化和机械化作业的发展[2]。发展和运用可代替人工打顶的化学打顶技术，具有广阔前景。探索化学打顶剂对棉花形态性状、激素调控及蛋白质差异表达的效应，将明确其调控的分子机制，为开发新型生长调节剂及棉花化学打顶技术的应用提供参考，具有重要的理论和实践意义。【前人研究进展】打顶是棉花栽培过程中的关键环节，化学打顶更是棉花轻简化栽培的必然趋势。前人就化学打顶与人工打顶对棉花的影响差异进行了许多研究，发现化学打顶可以起到类似人工打顶的效果，但不同打顶方式塑造的棉花株型不同[3]，同时对棉花内源激素含量[4]、物质分配[5-6]、代谢活动等产生的影响不同。袁青锋等[7]研究表明，化学打顶与人工打顶处理相比，棉花的株高、果枝数和结铃数等较高，两种处理的单铃重、衣分、籽棉和皮棉产量水平相当；杨成勋等[8]发现棉花化学打顶后，株宽显著低于人工打顶，改善了冠层中下部透光率，增加了光合面积和光合时间，利于产量的形成；康正华等[9]对棉花进行了不同化学打顶剂处理，表明部分化学打顶剂处理具有增产效果；李莉等[10]研究表明，打顶后加速植物生长的激素含量降低，促进植物成熟与衰老的激素含量升高，加速了棉花的衰老。氟节胺打顶处理后棉花倒四叶内源激素发生了不同程度的变化，与人工打顶处理相比变化趋势及作用时间存在差异，最终达到打顶的效果[11]。由上可知，激素在棉花生长发育过程中发挥着重要作用，研究植物激素的信号传导和作用机理也是蛋白质组学的重要内容之一。杨惠娟等[12]针对烟草打顶的研究发现，打顶前后有9个蛋白质差异表达，其中，促进烟草叶片光合作用和抗性相关蛋白的表达增强，与光呼吸能量代谢相关的蛋白表达下调。iTRAQ技术[13]是一种植物蛋白质应答分子机制有效研究技术。基于iTRAQ技术研究化学打顶条件下的蛋白质组学，有利于了解打顶后棉花的生理与分子反应机制。【本研究切入点】前人对化学打顶技术研究多集中于农艺性状、经济性状以及产量品质等方面，对棉花蛋白质组差异表达的影响鲜有报道。【拟解决的关键问题】通过设置人工与化学打顶两个处理，对棉花形态性状、激素调控机制及蛋白质差异表达情况开展研究，探究其生理调控的分子机制，为棉花轻简化栽培提供理论依据。

1 材料与方法

1.1 试验材料与田间设计

试验于2015—2016年在河北农业大学清苑实验农场进行，土质为壤土。2015年试验地的土壤基础地力情况：有机质含量19.79 g·kg-1、碱解氮含量70.84 mg·kg-1、速效磷含量35.42 mg·kg-1、速效钾含量149.45 mg·kg-1；2016年试验地的土壤基础肥力情况：有机质含量18.41 g·kg-1、碱解氮含量52.14 mg·kg-1、速效磷含量28.66 mg·kg-1和速效钾含量170.67 mg·kg-1。

供试品种为冀棉863，设人工打顶、化学打顶和未打顶3个处理。采用随机区组试验设计，3次重复，每个小区长12 m，宽11.5 m。宽窄行种植，大行行距100 cm，小行50 cm，密度45 000株/hm2。

所有处理于7月20日进行打顶处理。化学打顶剂为水乳剂，主要成分包括缩节胺（N,N- Dimethylpiperidinium chloride，DPC）、缓释剂、助剂等，由中国农业大学农学与生物技术学院作物化学控制研究中心和北京市农业技术推广站共同研发，新疆金棉科技有限责任公司生产。化学打顶剂用量为1.125 L·hm-2，人工喷施。其他田间管理措施采用常规高产栽培模式。

1.2 田间取样及调查

株高调查：在7月20日打顶后，每10 d选择10株代表性植株进行调查，至株高不再变化。

激素样品取样：于打顶日开始，每10 d取一次样品，共取6次，选取部位为主茎最高效叶片倒四叶，测定生长素（IAA）、赤霉素（GA3）、脱落酸（ABA）和玉米素核苷（ZR）含量。

蛋白质组学样品：在打顶后20 d（8月10日）选取人工打顶和化学打顶处理主茎倒四叶，叶片装入样品袋后，用液氮速冻后转至-80℃保存备用。

1.3 蛋白质组测定分析

1.3.1 激素含量测定 激素测定方法采用间接酶联免疫吸附法（ELISA）。试剂盒由中国农业大学化控室提供，使用美国BIO-RAD 680酶标仪测定OD值，并制作logit曲线计算样品中激素的含量。

1.3.2 蛋白质组学测定 蛋白提取与定量：样品用液氮磨成干粉，加入200 μL TEAB溶解；超声破碎15 min后，12 000 r/min离心20 min取上清；加入4倍体积的冷丙酮（含终浓度为10 mmol·L-1DTT）沉淀2 h；12 000 r/min离心20 min，收集沉淀；再加入800 μL的冷丙酮（含终浓度为10 mmol·L-1DTT）重悬沉淀；12 000 r/min离心20 min，收集沉淀，然后风干沉淀；加入100 μL TEAB溶解蛋白得到蛋白样品。选择常规的Bradford定量方法定量总蛋白质。

蛋白酶切与标记：选取100 μg蛋白质，体积整体调节到100 μL，然后加入500 μL 50 mmol·L-1NH4HCO3稀释，加入2 μg Tryspin 酶液；37℃消化过夜8—16 h；取出上述酶解液，加入等体积的0.1% FA酸化；将上述酸化后的酶解液加入到活化平衡后的Strata-X C18柱子中除盐，冷冻抽干后用20 μL 0.5 mol·L-1TEAB复溶。标记采用8-plex标记，人工打顶样品分别用115、116和117标记，化学打顶样品分别用118、119和121标记，6个样本标记完以后等量混合。

肽段分离与蛋白质鉴定：使用Thermo DINOEX Ultimate 3000 BioRS色谱仪进行分离，将混合后的样本分成12个组分。分离后的肽段进行LC-MS/MS分析，质谱仪是Thermo公司的Orbitrap质量检测器，型号为Q exactive HF。

蛋白质功能鉴定：对鉴定蛋白质进行GO、COG和KEGG分析，全面反映蛋白质功能。Gene ontology（GO）是一个国际标准化的基因功能分类体系，提供了一套动态更新的标准词汇表来全面描述生物体中基因和基因产物的属性。Cluster of orthologous groups of proteins（COG，蛋白相邻类的聚簇）是对蛋白质进行直系同源分类的数据库。Kyoto encyclopedia of genes and genomes（KEGG）是系统分析基因功能、基因组信息数据库，它有助于把基因及表达信息作为一个整体网络进行研究。

1.4 数据统计与分析

株高及激素含量数据采用Microsoft Excel 2010进行统计，通过SPSS 21.0进行方差分析检验其差异显著性。

蛋白质数据分析采用对MS/MS质谱数据经过系列优化处理后与数据库进行相似性比较打分从而进行蛋白鉴定。使用与AB Sciex 5600 plus配套的搜索引擎ProteinpilotTM V4.5。对于proteinpilot的鉴定结果进一步过滤，对于鉴定到的蛋白，认为unused score≥1.3（即可信度水平在95%以上），每个蛋白至少包含一个unique肽段的蛋白为可信蛋白；对于鉴定的肽段，以conf≥95过滤，即可信度在95%以上认为可信肽段，对于蛋白质定量，为得到更全面的关于某蛋白的定量信息，Proteinpilot软件使用了conf≥15的所有肽段，不符合该条件的肽段不包含在本报告中。

2 结果

2.1 不同打顶处理对棉花株高的影响

由图1可以看出，两年试验中3种打顶处理的株高从处理后20 d开始出现差异显著，表现为未打顶处理最高，化学打顶处理次之，人工打顶最低。

2.2 不同打顶处理对棉花功能叶激素含量的影响

由图2可以看出，2015年试验中3种打顶处理棉花功能叶中GA3含量呈先升高后降低趋势，在处理后20 d时未打顶处理含量显著高于其他两个处理，在处理后30 d时化学打顶处理显著高于人工打顶处理，3种打顶处理之间差异显著，处理后40 d时化学打顶与未打顶处理无显著差异，显著高于人工打顶处理。打顶处理后，棉花叶片IAA和ABA含量先升高，在处理后40 d时达到最大值，之后逐渐下降，在处理后30 d时未打顶处理IAA含量显著低于其他两个处理，未打顶处理ABA含量峰值显著低于其他两个处理，其余时期无显著差异。ZR含量在处理后20 d内逐渐升高，之后在250 ng·g-1值上下波动，3种打顶处理之间无显著差异。

图中处理间标有相同字母表示差异不显著（P＞0.05）。下同

图2 2015年不同处理对棉花主茎功能叶激素含量的影响

由图3可以看出，2016年试验中3种打顶处理GA3含量变化差异较大，在处理后30 d时未打顶处理显著高于其他两个处理，在处理40 d之后未打顶＞化学打顶＞人工打顶，3种处理之间达到显著差异。3种打顶处理的IAA和ABA含量变化均为单峰曲线，在处理后40 d时达到最大值，化学打顶处理IAA含量峰值显著低于其他两个处理，未打顶处理ABA含量峰值显著低于其他两处理。ZR含量呈先升高后降低的趋势，在处理后40 d达到最大值，3种打顶之间无显著差异。

综上可知，不同打顶处理对棉花功能叶GA3含量影响较大，对IAA和ABA含量有一定影响，对ZR含量基本无影响。与人工打顶处理相比，化学打顶与未打顶处理后期叶片GA3含量下降相对较晚。

图3 2016年不同处理对棉花主茎功能叶激素含量的影响

2.3 不同打顶处理对棉花蛋白质组表达影响

2.3.1 蛋白质定量分析 对人工打顶与化学打顶处理棉花主茎功能叶进行iTRAQ蛋白质组学分析，以人工打顶处理叶片为对照，比较两者之间的差异表达蛋白，当差异倍数达到1.2倍及以上（即up-regulate≥1.2和down-regulate≤0.83），且经过显著性统计检验其-value≤0.05时，视为显著差异蛋白。结果显示总鉴定蛋白为3 065个，两个处理之间检测到的总差异蛋白质数量为69个，其中29个上调，40个下调表达。

本研究中，质谱产生的二级谱图数、解析的二级谱图数分别为352 145和76 544，谱图鉴定率达到了21%以上，鉴定到的肽段和蛋白数总体分别为14 676和3 147（表1）。

2.3.2 鉴定蛋白质功能注释 通过对鉴定蛋白质进行GO、KEGG和COG注释，以利用不同来源的数据库全面反映蛋白质的功能，揭示蛋白质在各个生命活动中的生物学意义。从图4可以看出，由于背景注释库的局限性，GO功能注释到的蛋白质最多，为2 996个蛋白质，COG注释到2 159个蛋白质，KEGG注释到1 858个蛋白质。

2.3.3 鉴定蛋白质的COG注释分析 COG是对蛋白质进行直系同源分类的数据库，通过对鉴定到的蛋白质同COG数据库进行了比对，对这些蛋白质进行了功能分类（图5）。其中，1.9%蛋白质功能未知，10.5%参与能量代谢，7%参与氨基酸运输及代谢，9%涉及碳水化合物运输及代谢，其他差异蛋白可分别归入翻译后修饰、翻译、辅酶运输与代谢、信号传导等方面。

图4 不同功能注释结果统计图

表1 蛋白质鉴定信息统计表

*表示可信度至少为95%，**表示至少含有2个unique肽段的鉴定蛋白质数目

* means significant difference at 0.05 level, ** means the number of identified proteins containing 2 unique peptides at least

A：核糖核酸加工与修饰RNA processing and modification；B：染色质结构和动态变化Chromatin structure and dynamics；C：能量产物与转换Energy production and conversion；D：细胞循环调控、细胞分裂、染色体分裂Cell cycle control, cell division, chromosome partitioning；E：氨基酸转运与代谢Amino acid transport and metabolism；F：核苷酸转运与代谢Nucleotide transport and metabolism；G：碳水化合物转运与代谢Carbohydrate transport and metabolism；H：辅酶转运与代谢Coenzyme transport and metabolism；I：脂质转运与代谢Lipid transport and metabolism；J：翻译、核糖体结构和生物合成Translation, ribosomal structure and biogenesis；K：翻译Translation；L：复制、重组和修饰Replication, recombination and repair；M：细胞壁/薄膜/膜的生物合成Cell wall/membrane/envelope biogenesis；N：细胞活力Cell motility；O：翻译后修饰、蛋白质周转、伴侣Posttranslational modification, protein turnover, chaperones；P：无机离子转运与代谢Inorganic ion transport and metabolism；Q：次级代谢产物的合成、转运与分解Secondary metabolites biosynthesis, transport and catabolism；R：一般功能预测General function prediction only；S：未知功能Function unknown；T：信号转导机制Signal transduction mechanisms；U：细胞内转运、分泌和液泡运输Intracellular trafficking, secretion, and vesicular transport；V：防御机制Defense mechanisms；Z：细胞骨架Cytoskeleton

2.3.4 鉴定和差异蛋白质GO分析 为了能够更加清楚的分析差异蛋白质的功能，对差异蛋白质进行独立的功能注释分析（图6）。GO总共有3个本体，分别描述基因的分子功能、细胞组分、参与的生物过程。在生理过程中，有19个功能存在明显差异，细胞成分组织或合成、细胞过程、代谢过程、生物调节等功能涉及的差异蛋白比例较大。在细胞组分中，细胞、细胞组分、膜腔、细胞器等两处理之间差异蛋白显著。在分子功能中，两处理在抗氧化物酶活性、结合、催化酶活性和结构分子活性等方面都出现了显著差异。

A：生物附着biological adhesion；B：生物调节biological regulation；C：碳固定carbon utilization；D：细胞增殖cell proliferation；E：细胞成分组成或生物合成cellular component organization or biogenesis；F：细胞过程cellular process；G：死亡death；H：发育过程developmental process；I：建立定位establishment of localization；J：生长growth；K：免疫系统过程immune system process；L：定位localization；M：转移locomotion；N：代谢过程metabolic process；O：多机体过程multi-organismal process；P：多细胞生物过程multicellular organismal process；Q：生物过程的负调节negative regulation of biological process；R：色素沉着pigmentation；S：生物过程的正调节positive regulation of biological process；T：生物过程调节regulation of biological process；U：复制reproduction；V：繁殖过程reproductive process；W：刺激响应response to stimulus；X：周期过程rhythmic process；Y：信号signaling；Z：病毒繁殖viral reproduction；AB：细胞cell；AC：细胞组分cell part；AD：胞外区extracellular region；AE：胞外区组分extracellular region part；AF：高分子配合物macromolecular complex；AG：膜封闭腔membrane-enclosed lumen；AH：细胞器organelle；AI：细胞器组分organelle part；AJ：抗氧化活性antioxidant activity；AK：结合binding；AL：催化活性catalytic activity；AM：电子载体活性electron carrier activity；AN：酶调节剂活性enzyme regulator activity；AO：分子传导活性molecular transducer activity；AP：核酸结合转录因子活性nucleic acid binding transcription factor activity；AQ：蛋白结合转录因子活性protein binding transcription factor activity；AR：结构分子活性structural molecule activity；AS：转移载体活性transporter activity

2.3.5 鉴定差异蛋白质基本信息 表2列出了所有差异蛋白质的相关信息。差异表达蛋白质为69个，其中涉及碳水化合物和能量代谢的蛋白质有14个，涉及蛋白质和氨基酸代谢有16个，应激与防御相关蛋白有14个，还涉及其他代谢过程。主要富集的Pathway途径有精氨酸和脯氨酸代谢、光合作用碳固定过程、谷胱甘肽代谢和苯丙氨酸代谢等。能量代谢相关蛋白多下调表达，应激相关蛋白多上调表达。

3 讨论

3.1 不同打顶处理下棉花株高与激素含量的差异

由本研究结果可知，不同打顶处理后棉花株高出现显著差异，处理后期化学打顶处理显著高于人工打顶，但显著低于未打顶处理，这可能是喷施化学打顶剂后棉花需要一定的时间吸收转化，在这期间棉花继续生长，这与袁青锋等[7-8]研究结果一致。这说明化学打顶处理能在一定程度上控制棉花株高。不同打顶处理后棉花主茎功能叶激素含量也表现出一定差异，但两年结果中3种打顶处理的ZR含量在各时期差异不显著。两年试验中3种打顶处理前期各激素含量基本无差异，这与赵强[14]研究结果一致，但与叶春秀等[11]研究结果不相同；前人研究表明人工打顶导致棉花生长中后期GA合成量减少[10]，本试验后期化学打顶的GA3含量显著高于人工打顶，说明化学打顶可以缓解人工打顶后期GA含量快速下降，可能利于防止棉花早衰，有待于进一步研究证实。

表2 差异蛋白质相关信息

3.2 不同打顶处理下棉花功能叶的蛋白质组差异

在蛋白质组水平上，应用iTRAQ技术对两种打顶方式下棉花主茎功能叶进行了蛋白质组学分析，鉴定出了一定数量差异蛋白，涉及与植物生长发育密切相关的碳水化合物和能量代谢等诸多方面，一些与GA代谢相关的蛋白质表达丰度也发生了变化。

3.2.1 与能量代谢及应激代谢相关的蛋白质 打顶处理后，棉花主茎功能叶中与参与碳水化合物和能量代谢相关的许多蛋白质的表达丰度发生变化，涉及糖酵解、三羧酸循环、卡尔文循环等碳水化合物代谢途径。在三羧酸循中，乌头酸水合酶是三羧酸循环中的辅酶，催化柠檬酸和异柠檬酸相互转化。柠檬酸合酶催化来自糖酵解或其他异化反应的乙酰CoA与草酰乙酸缩合合成柠檬酸，是TCA循环的关键酶。本研究中顺乌头酸水合酶（aconitate hydratase）和ATP-柠檬酸合成酶链蛋白2（ATP-citrate synthase beta chain protein 2）表达上调，这表明化学打顶与人工打顶处理相比，促进了三羧酸循环，加快了碳水化合物代谢。V型质子ATP酶催化亚基A蛋白（V-type proton ATPase catalytic subunit A-like protein）表达上调，促进了短链脂肪酸的代谢。

在C3植物叶片中，丙酮酸磷酸双激酶（pyruvate orthophosphate dikinase，PPDK）是糖酵解途径的关键酶，参与碳固定，并为氨基酸合成提供碳骨架[15]。3-磷酸甘油醛脱氢酶在细胞的碳代谢的中心环节，催化糖酵解中唯一的氧化反应[16]。本试验中丙酮酸磷酸双激酶与3-磷酸甘油醛脱氢酶2（glyceraldehyde-3- phosphate dehydrogenase 2）均表达下调，这表明与人工打顶相比，化学打顶处理的叶片中糖酵解途径减弱，这说明人工打顶后棉花需要更多的能量来响应，棉花化学打顶有效缓解棉花因打顶造成的剧烈反应，从而维持棉花的持续稳定增长。

谷胱甘肽代谢途径对维持细胞内活性氧与自由基动态平衡、抵御胁迫具有重要作用[17]。此途径中谷胱甘肽-S-转移酶F9样蛋白（glutathione S-transferase F9-like protein）和谷胱甘肽合成酶（glutathione synthetase）表达上调，谷胱甘肽还原酶（glutathione reductase）表达下调，总体上加快了谷胱甘肽的合成以及对有害物质的清除。精氨酸和脯氨酸代谢途径中delta-1-吡咯啉-5-羧酸脱氢酶12A1（delta-1-pyrroline- 5-carboxylate dehydrogenase 12A1）和天冬氨酸转氨酶（aspartate aminotransferase）表达上调，加快了三羧酸循环，乙酰谷氨酸激酶（acetylglutamate kinase）表达下调，对尿素循环造成了一定影响，精氨酸合成多胺调节植物的生长发育，从而提高棉花植株的抗逆能力[18]。

前人研究表明菠菜受到硫化氢胁迫后，光合作用和能量代谢提高来响应胁迫[19]，本试验中化学打顶处理后棉花叶片中碳水化合物和能量代谢相关蛋白多数下调表达，抗逆相关蛋白质表达上调，这说明化学打顶对棉花的机械损伤更小，人工打顶后棉花提高能量代谢来抵御打顶造成的损伤，但化学打顶后可能产生抑制植株生长的物质，棉株通过自身抗逆反应来适应自身变化。

3.2.2 与GA代谢相关的蛋白质 DELLA蛋白在GA和IAA、ABA等其他植物激素信号传导途径中具有交互调节作用，是GA调控作用的关键，其在其他环境信号间的调控机制也被逐一发现[20]，如光的信号转导[21]。DELLA蛋白是一类被广泛定性的GA信号抑制因子[22]。GID1蛋白与GA及DELLA蛋白结合后形成赤霉素受体GA-GID1-DELLA蛋白复合体，然后SCF聚合体标记该三聚体，诱导泛素26S蛋白酶体降解DELLA蛋白，解除DELLA蛋白对植物生长的抑制作用，产生GA效应[23-26]。该试验中26S蛋白酶体的非ATP酶调节亚基6同源物（26S proteasome non-ATPase regulatory subunit 6 homolog）表达上调，加快降解DELLA蛋白，促进植株生长。因此，与人工打顶处理相比，化学打顶处理更能促进棉花的纵向生长。

ABA抑制GA诱导的一些过程，如PKABA蛋白激酶是ABA正调节相关蛋白，它过表达可抑制GA调节的GAMYB和-淀粉酶基因的表达[27-28]。GA也调节ABA的代谢和信号转导[29]，HvSPY是GA负调控因子，却在ABA信号转导中起正调节作用，但不依赖ABA途径[30]。本研究中-双加氧酶1（Alpha-dioxygenase 1）参与ABA负调控介导的信号转导通路，其表达上调，可能抑制ABA表达，延缓植物衰老。镁螯合酶亚基Chl I-1（Magnesium-chelatase subunit ChlI-1）是叶绿素生物合成的关键酶，正调节ABA信号，但不与ABA结合[31]，该蛋白表达上调，促进ABA产生效应。并且GA的变化会直接影响IAA的生物合成与降解，IAA又直接作用于细胞，导致细胞或者细胞壁的伸长变异，这又会导致植株节间长度或者节数变化，最终导致植株形态变化。此外，与光合作用相关的蛋白质光合NDH亚基的亚基B 2（photosynthetic NDH subunit of subcomplex B 2）和叶绿素a、b结合蛋白CP26（chlorophyll a-b binding protein CP26）参与控制碳水化合物代谢，影响蔗糖利用效率和细胞壁的合成，其表达下调导致化学打顶处理植株光合能力减弱，本研究中后期人工打顶处理棉花主茎叶GA含量显著低于化学打顶处理，光合能力也受到抑制，这与前人研究结果一致[32-33]。这说明打顶存在对GA调控的抑制，但化学打顶对GA的抑制比人工打顶要弱，GA减少会加速叶片的衰老，化学打顶能在一定程度上防止棉花早衰。

前人研究表明GA与类胡萝卜素合成途径有共同的前体，并且有些类胡萝卜素能直接合成ABA[34]。糖酵解产物乙酰辅酶A和丙酮酸促进前体物质IPP的生成，加快了植物MEP途径，为棉花叶片GA含量升高提供了条件。GA代谢及信号传导正相关蛋白多上调，与人工打顶相比，化学打顶促进了GA的合成与表达，植物通过减少赤霉素来减缓植株生长从而适应环境，这同样说明化学打顶与人工打顶相比，对植株造成的伤害更小，更有利于植株的生长发育，但关于不同打顶方式影响棉花生长的具体分子机制还有待于进一步研究。

4 结论

与人工打顶相比，化学打顶处理有效控制了棉花株高，提高了生育后期GA含量，IAA、ABA和ZR含量无显著差异。基于iTRAQ技术的蛋白质组学分析，鉴定得到69个差异表达蛋白质，与植物生长发育相关的碳水化合物和能量代谢相关蛋白多下调表达，与GA调控正调节相关蛋白多上调表达，人工打顶激活了植物的能量代谢与应激反应来调控棉花的生长，而化学打顶更多的是通过调控GA激素的代谢过程来适应外部环境的变化。

[1] 肖才升, 李育强, 李庠, 陈浩东, 王洪, 何叔军, 匡政成. 洞庭湖植棉区油后直播棉的发展前景探讨. 中国棉花, 2013, 40(11): 21-22.

XIAO C S, LI Y Q, LI X, CHEN H D, WANG H, HE S J, KUANG Z C. The development prospects of cotton directly seeded after rape in Dongting Lakeside region., 2013, 40(11): 21-22. (in Chinese)

[2] 毛树春. 我国棉花种植技术的现代化问题. 中国棉花, 2010, 37(3): 2-5.

MAO S C. The modernization problems of cotton planting technology in China., 2010, 37(3): 2-5. (in Chinese)

[3] 李雪, 朱昌华, 夏凯, 甘立军. 辛酸甲酯, 癸酸甲酯和6-BA对棉花去顶的影响. 棉花学报, 2009, 21(1): 70-72.

LI X, ZHU C H, XIA K, GAN L J. Effects of methyl octanoate, methy decanoate and 6-BA on topping in cotton., 2009, 21(1): 70-72. (in Chinese)

[4] 李莉, 田长彦, 陈冠文, 李中军. 棉花打顶后倒四叶内激素变化对棉花产量和品质的影响. 干旱区研究, 2014, 31(1): 111-117.

LI L, TIAN C Y, CHEN G W, LI Z J. Effect of hormone change in reciprocal fourth leaf on cotton yield and quality after topping., 2014, 31(1): 111-117. (in Chinese)

[5] 毛树春, 董和林, 裴建忠. 棉花栽培新技术. 上海: 上海科学技术出版社, 2002: 78-110.

MAO S C, DONG H L, PEI J Z.. Shanghai: Shanghai science and Technology Press, 2002: 78-110. (in Chinese)

[6] 叶伟. 打顶、化控对棉花蕾铃形成及脱落的影响. 安徽农学通报, 2009, 15(13): 93-141.

YE W. The effect of topping and chemical control on the formation and shedding of cotton buds., 2009, 15(13): 93-141. (in Chinese)

[7] 袁青锋, 张静, 管利军, 徐宇强, 朱晓平. 化学打顶对棉花生长发育和产量的影响. 江苏农业科学, 2015, 43(5): 72-74.

YUAN Q F, ZHANG J, GUAN L J, XU Y Q, ZHU X P. Effects of chemical topping on cotton growth and yield., 2015, 43(5): 72-74. (in Chinese)

[8] 杨成勋, 张旺锋, 徐守振, 随龙龙, 梁福斌, 董恒义. 喷施化学打顶剂对棉花冠层结构及群体光合生产的影响. 中国农业科学, 2016, 49(9): 1672-1684.

YANG C X, ZHANG W F, XU S Z, SUI L L, LIANG F B, DONG H Y. Effects of spraying chemical topping agents on canopy structure and canopy photosynthetic production in cotton., 2016, 49(9): 1672-1684. (in Chinese)

[9] 康正华, 赵强, 娄善伟, 葛杰, 张海祥. 不同化学打顶剂对棉花农艺及产量性状的影响. 新疆农业科学, 2015, 52(7): 1200-1208.

KANG Z H, ZHAO Q, LOU S W, GE J, ZHANG H X. Effects on the agronomic and economic characters of cotton by applying different topping chemicals., 2015, 52(7): 1200-1208. (in Chinese)

[10] 李莉, 田长彦, 黄子蔚, 王平, 陈冠文, 王强, 李中军. 打顶后棉花倒四叶内源激素的变化及与早衰的关系. 干旱地区农业研究, 2010, 28(3): 162-168.

LI L, TIAN C Y, HUANG Z W, WANG P, CHEN G W, WANG Q, LI Z J. Influence of decapitation on endogenous hormone of inverse fourth leaf and premature senility of cotton., 2010, 28(3): 162-168. (in Chinese)

[11] 叶春秀, 王志军, 王刚, 庄振刚, 谢宗铭. 化学打顶剂对北疆早熟陆地棉倒四叶内源激素含量的影响//中国农学会棉花分会2016年年会论文汇编, 2016: 7.

YE C X, WANG Z J, WANG G, ZHUANG Z G, XIE Z M. Effects of chemical topping on endogenous hormones contents in early maturing upland cotton () in Northern Xinjiang//, 2016: 7. (in Chinese)

[12] 杨惠娟, 许俐, 王景, 刘国顺, 史宏志, 崔红, 钟伯雄. 打顶前后烤烟叶片蛋白质表达差异的研究. 中国烟草学报, 2012, 18(5): 73-78.

YANG H J, XU L, WANG J, LIU G S, SHI H Z, CUI H, ZHONG B X. An investigation into differences in protein expression in flue-cured tobacco leaf before and after topping., 2012, 18(5): 73-78. (in Chinese)

[13] PIERCE A, UNWIN R D, EVANS C A, GRIFFITHS S, CARNEY L, ZHANG L Q, JAWORSKA E, LEE C F, BLINCO D, OKONIEWSKI M J, MILLER C J, BITTON D A, SPOONCER E, WHETTON A D. Eight-channel iTRAQ enables comparison of the activity of six leukemogenic tyrosine kinases., 2008, 7(5): 853-863.

[14] 赵强. 棉花化学打顶及其配套技术体系研究[D]. 北京: 中国农业大学, 2011.

ZHAO Q. Study on the characteristics and technical systems of chemical detopping in cotton (L.)[D]. Beijing: China Agricultural University, 2011. (in Chinese)

[15] DOUBNEROVÁ V, RYŠLAVÁ H. What can enzymes of C4 photosynthesis do for C3 plants under stress?, 2011, 180(4): 575-583.

[16] LAO Y M, LU Y, JIANG J G, LUO L X. Six regulatory elements lying in the promoter region imply the functional diversity of chloroplast GAPDH in., 2012, 60(36): 9211-9220.

[17] 闫慧芳, 毛培胜, 夏方山. 植物抗氧化剂谷胱甘肽研究进展. 草地学报, 2013, 21(3): 428-434.

YAN H F, MAO P S, XIA F S. Research progress in plant antioxidant glutathione., 2013, 21(3): 428-434. (in Chinese)

[18] TAKAHASHI T, KAKEHI J I. Polyamines: ubiquitous polycations with unique roles in growth and stress responses., 2010, 105(1): 1-6.

[19] 张川芳, 赵琪, 肖振, 蔡晓峰, 王全华, 戴绍军. 菠菜发育与逆境应答蛋白质组学研究进展. 园艺学报, 2016, 43(9): 1752-1764.

ZHANG C F, ZHAO Q, XIAO Z, CAI X F, WANG Q H, DAI S J. Advances in proteomics analysis ofdevelopment and stress response., 2016, 43(9): 1752-1764. (in Chinese)

[20] ACHARD P, CHENG H, GRAUWE L D, DECAT J, SCHOUTTETEN H, MORITZ T, DOMINIQUE V D S, PENG J R, HARBERD N P. Integration of plant responses to environmentally activated phytohormonal signals., 2006, 311(5757): 91-94.

[21] GAO S J, XIE X Z, YANG S G, CHEN Z P, WANG X J. The changes of GA level and signaling are involved in the regulation of me-socotyl elongation during blue light mediated de-etiolation in., 2012, 39(4): 4091-4100.

[22] DILL A, JUNG H S, SUN T P. The DELLA motif is essential for gibberellin-induced degradation of RGA., 2001, 98(24): 14162-14167.

[23] SASAKI A, ITOH H, GOMI K, UEGUCHI-TANAKA M, ISHIYAMA K, KOBAYASHI M, JEONG D H, AN G, KITANO H, ASHIKARI M, MATSUOKA M. Accumulation of phosphorylated repressor for gibberellin signaling in an F-box mutant., 2003, 299(5614): 1896-1898.

[24] MCGINNIS K M, THOMAS S G, SOULE J D, STRADER L C, ZALE J M, SUN T P, STEBER C M. TheSLEEPY1 gene encodes a putative F-box subunit of an SCF E3 ubiquitin ligase., 2003, 15(5): 1120-1130.

[25] VIERSTRA R D. The ubiquitin/26S proteasome pathway, the complex last chapter in the life of many plant proteins., 2003, 8(3): 135-142.

[26] UEGUCHI-TANAKA M, NAKAJIMA M, KATOH E, OHMIYA H, ASANO K, SAJI S, XIANG H Y, ASHIKARI M, KITANO H, YAMAGUCHI I, MATSUOKA M. Molecular interactions of a soluble gibberellin receptor, GID1, with a rice DELLA protein, SLR1, and gibberellin., 2007, 19(7): 2140-2155.

[27] ZENTELLA R, YAMAUCHI D, HO T D. Molecular dissection of the gibberellin/abscisic acid signaling pathways by transiently expressed RNA interference in barley aleurone cells., 2002, 14(9): 2289-2301.

[28] GÓMEZ-CADENAS A, ZENTALLA R, WALKER-SIMMONS M K, Ho T H D. Gibberellin/abscisic acid antagonism in barley aleurone cells: site of action of the protein kinase PKABA1 in relation to gibberellin signaling molecules., 2001, 13(3): 667-679.

[29] ITOH H, UEGUCHI-TANAKA M, SATO Y, ASHIKARI M, MATSUOKA M. The gibberellin signaling pathway is regulated by the appearance and disappearance of SLENDER RICE1 in nuclei., 2002, 14(1): 57-70.

[30] ROBERTSON M. Increased dehydrin promoter activity caused by HvSPY is independent of the ABA response pathway., 2003, 34(1): 39-46.

[31] 杜淑媛. 镁螯合酶各组分在ABA信号转导中的功能[D]. 北京: 清华大学, 2012.

DU S Y. Roles of different components of magnesium chelatase in abscisic acid signal transduction[D]. Beijing: Tsinghua University, 2012. (in Chinese)

[32] 杜尧舜. 增施CO2和生长调节剂对草莓光合作用的影响. 浙江农业学报, 2000, 12(3): 144-146.

DU Y S. Effects of CO2enrichment and plant growth regulators on photosynthesis in strawberry., 2000, 12(3): 144-146. (in Chinese)

[33] YUAN L, XU D Q. Stimulatory effect of exogenous GA3on photosynthesis and the level of endogenous GA3in soybean leaf., 2001, 28(4): 317-320.

[34] BARTLEY G E, SCOLNIK P A. Plant carotenoids: Pigments for photo protection, visual attraction, and human health., 1995, 7: 1027-1038.

（责任编辑 李莉，岳梅）

Differential Proteomics Analysis of Cotton Leaf Response to Chemical Topping Based on iTRAQ Technique

HOU XiaoMeng1,2, LIU LianTao1,2, LI Meng4, SUN HongChun1,2, ZHANG YongJiang1,2, DU Huan2,3, LI CunDong2

(1College of Agronomy, Hebei Agricultural University, Baoding 071000, Hebei;2Key Laboratory of Hebei Crop Growth Regulation, Baoding 071000, Hebei;3College of Life Sciences, Hebei Agricultural University, Baoding 071000, Hebei;4Hebei Academy of Agriculture and Forestry Sciences, Shijiazhuang 050051)

【】The objective of this study is to compare the effects of manual topping and chemical topping on physiological changes and proteins differentially expression of cotton. The results of this study would provide a theoretical basis for the mechanism of chemical topping. 【】Manual topping, chemical topping and no top pruning were practiced on cultivar GM863 which was planted in large area of the Yellow River basin in 2015-2016. Cotton was treated on July 20th simultaneously. Artificial spraying chemical topping agent and the dosage was 1.125 L·hm-2. The plant height and hormone contents of functional leaves were determined regularly after topping. The plant height was measured with the height from the cotyledon node to the growth point of the main stem by using ruler. The hormone contents of functional leaves, including IAA, GA3, ABA and ZR were measured by enzyme-linked immunosorbent assay. The differential proteomics analysis was carried out on the functional leaves of manual topping and chemical topping by using iTRAQ technique. 【】Compared with the cotton of manual topping treatment, chemical topping treatment was significantly higher than manual topping treatment, which was 11.8% and 14.5% higher than others in the two-year trial, but was significantly lower than no top pruning, which was 6.0% and 6.5% lower than others, vertical growth was significantly inhibited after spraying. Different topping treatments greatly affected the GA3content of cotton functional leaves, the change of GA3content was a single peak curve, the contents of GA3were no top pruning＞chemical topping＞manual topping, and reached significant difference. The GA3contents of chemical topping and no top pruning were decreased at 30 d after topping and manual topping was decreased at 20 d. There was no significant difference among three treatments at 50 d after topping. The maximum IAA content in chemical topping treatment was significantly lower than that in the other two treatments at 40 d after topping in 2016, and there was no significant difference among the three treatments in 2015. The maximum ABA content in no top pruning treatment was significantly lower than that in the other two treatments at 40 d after topping. There was no significant difference among the three treatments of the contents of ZR. Compared with manual topping treatment, 69 differentially expressed proteins were detected by using iTRAQ in leaves of chemical topping treatment, 29 up-regulated and 40 down-regulated expressions, in which carbohydrate and energy metabolism related proteins were most down-regulated, and it was indicated that the growth of plants was decreased. The up-regulation of most GA regulated positive correlation proteins enhanced the GA effect. 【】Chemical topping can effectively control the plant height of cotton, and which has a great effect on the GA content of functional leaves. The contents of GA in chemical topping treated leaves were significantly higher than that of manual topping treatment. Compared with manual topping treatment, the proteins related to plant growth and development in chemical topping treatment were most down-regulated expression, may be the plants reduced carbohydrate synthesis and energy metabolism, and activated GA effect by increasing GA content to control plant height.

cotton; chemical topping; hormone; iTRAQ technique; differential proteins

10.3864/j.issn.0578-1752.2017.19.004

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

国家自然科学基金（31171495，31571610）

侯晓梦，E-mail：houxiaomeng66@sina.com。通信作者李梦，E-mail：bdlimeng@163.com。通信作者李存东，E-mail：nxylcd@mail.hebau.edu.cn