木薯（Manihot esculenta Crantz）是世界上重要的粮食作物，块根可为非洲、亚洲以及中美洲提供粮食安全保障^[1]。木薯较其它作物的优势为产量高，抗旱能力强^[2]。在热带地区，木薯是继水稻和玉米之后重要的膳食来源，是全世界8亿人的主要食物来源^[3]。2014年我国木薯种植面积为38万hm²，鲜薯总产量约780万t；单产约20 t/hm²，离木薯理论产量的90 t/hm²还有很大的增长空间。因此研究影响木薯产量的因素，为进一步提高我国木薯单产提供理论依据是必要的。

木薯基因组高度杂合，后代性状分离严重，还具有自交不亲和、种子数量少等特性，导致杂交选育种效率低，这些特性已经严重阻碍木薯选育高产品种的进程^[4]。目前仅依靠传统的杂交育种手段很难快速提高木薯的选育种周期，近年来利用高通量测序技术开展野生木薯近缘种M. esculenta ssp. flabellifolia（W14）与栽培品种M. esculenta ssp. esculenta（KU50）比较基因组学研究，完成了W14和KU50的全基因组草图，注释了碳流、淀粉积累和氢氰酸合成代谢通路关键基因的生物学功能，并成功开发出数百万个全基因组分子标记^[5]，为选育种提供全基因组的平台。但通过全基因组数据只能间接推测蛋白质的功能，因为从基因表达的mRNA水平到最终合成蛋白质水平，其中包含着蛋白质翻译的调控、糖基化、磷酸化等诸多因素，这些因素都有可能改变作为直接作用因子的蛋白质功能，因而直接研究蛋白质变化具有不可替代的意义^[4]。

目前对木薯全蛋白质水平的研究还很不全面，因此本研究选用木薯高产栽培种ZM-Seaside和低产种质花叶变种作为研究材料，从全蛋白质的角度揭示这两个种质块根产量差异的主要原因，挖掘影响木薯产量的关键蛋白质，为选育高产木薯品种提供基础数据。

1 材料与方法 1.1 材料

本试验研究材料木薯高产栽培种质ZM-Seaside和低产种质花叶变种均来自中国热带农业科学院热带作物品种资源研究所木薯种质资源圃^[6]。以种植10个月后收获的木薯块根作为研究对象。

1.2 方法 1.2.1 产量的测定

种植后10个月，测定花叶变种和木薯栽培种ZM-Seaside的单株鲜薯重。每个品种选取15株。

1.2.2 木薯块根淀粉含量的测定

淀粉测定采用旋光法^[7]；应用W ZZ-1型旋光仪测定旋光物的旋光度。按照公式SC=（a×100）×100/（L×203×m）计算木薯的淀粉含量，式中SC为鲜薯粗淀粉含量（%），a为旋光度读数（度），L为观测管长度（2 dm），m为样品质量（g），203为淀粉的比旋光度。

1.2.3 木薯块根氢氰酸含量的测定

氢氰酸含量的测定采用硝酸银滴定法^[8]。使木薯浸水过夜发酵析出氢氰酸，然后将此溶液通入蒸汽蒸馏出氢氰酸，用过量的硝酸银标准溶液吸收蒸馏出来的氢氰酸，最后以标定好的硫氰化钾滴定多余的硝酸银溶液；由硝酸银用量与剩余硝酸银之差即可算出样品中氢氰酸（HCN）含量。

式中，V1为用硫氰化钾滴定25 mL硝酸银时消耗的体积（mL）；V2为滴定剩余硝酸银时消耗的体积（mL）；c为标准硫氰酸钾的浓度（mol/L）；27为氢氰酸的摩尔质量（g/mol）；m为木薯样品质量（g）。

1.2.4 木薯块根干物质率的测定

根据国际热带农业科学院中心制定的公式计算干物率^[9]。取一定质量的鲜薯于烧杯中，放在60℃干燥箱里烘干，大约4-5 d后每天测定重量直至衡重，记录衡重时木薯的重量。干物率=薯干重/鲜薯重

1.2.5 木薯块根蛋白质表达水平的分析

利用Western blot方法^[10]对3种与淀粉积累相关的蛋白质UGPase（购自Agrisera公司，货号AS05086，分子量：51.6 kD）、AGPase（购自Agrisera公司，货号AS111739，分子量：49.4 kD）和SPS（购自Agrisera公司，货号AS03035A，分子量：120-130 kD）进行表达水平分析，用Actin（购自Agrisera公司，货号AS132640，分子量：45kD）作为对照。利用ChemiImager 4400软件计算蛋白质的相对表达含量。

1.2.6 木薯块根全蛋白质的提取、分离和鉴定

木薯ZM-Seaside和花叶变种块根全蛋白质的提取均采用Chen等^[11]苯酚提取法，溶解后用Bradford试剂盒进行定量，后参照Chen等^[11]双向电泳技术将蛋白质进行分离。以花叶变种块根全蛋白质图谱为对照，采用Delta2D软件确定ZM-Seaside块根的差异蛋白质点，对平均差异表达量在±2.0以上的蛋白质进行标记^[12]。蛋白质的鉴定参照An等^[2]方法。本研究设3个生物学重复进行数据统计分析。

1.2.7 木薯蛋白质互作网络构建

采用String软件对鉴定出的差异蛋白质构建蛋白质互作网络^[13]。

1.3 数据统计分析

数据处理及统计学分析采用Excel 2010和DPS v7.55软件，差异显著性分析采用新复极差法（Duncan）^{[14, 15]}。

2 结果 2.1 木薯块根氢氰酸含量、淀粉含量及干物率

将花叶变种以及ZM-Seaside块根的薯肉跟薯皮分开处理，分别测定它们的氢氰酸含量。结果（图 1）显示，ZM-Seaside块根皮中氢氰酸含量（70.09 mg/kg）要高于花叶变种块根皮（62.61 mg/kg）中的氢氰酸含量，两种木薯块根薯肉的氢氰酸含量均较低，且没有显著差异。同时由图 1可知木薯氢氰酸主要存在于木薯块根皮中，根据NY/T 875规定食用木薯鲜薯薯肉中氢氰酸残留量不能超过50 mg/kg，因此这两种木薯块根均可食用。木薯收获后，测定得到花叶变种粗淀粉含量（25.83%）显著低于ZM-Seaside粗淀粉含量（29.18%）（图 2）；两种木薯的干物率为分别为37.16%及40.28%，栽培种ZM-Seaside的干物率显著高于花叶变种。由此可知，栽培种ZM-Seaside淀粉含量和干物率均显著高于花叶变种，进而导致其产量显著高于花叶变种。

2.2 与淀粉合成相关蛋白质表达分析

本研究利用Western blot方法研究与淀粉积累相关的蛋白质UGPase、AGPase和SPS分别在花叶变种和ZM-Seaside块根的表达水平（图 3）。将Western blot中每一种蛋白质在两个木薯种质中的总量定为100%，计算同一种蛋白质在两个木薯种质间的相对表达量。同时以Actin作为对照，确保两个木薯种质块根的上样量一致（图 3-A）。研究结果表明ZM-Seaside块根UGPase（图 3-B）、SPS（图 3-C）和AGPase（图 3-D）的表达水平显著高于花叶变种（P < 0.05）。该研究结果与木薯块根产量分析结果相一致，因而从与淀粉积累相关蛋白质水平进一步验证2种木薯块根产量的差异。

2.3 木薯块根全蛋白的提取、分离和鉴定

通过苯酚沉淀法提取块根全蛋白质，后经定量以及双向电泳分离、染色后，得到重复性较好的花叶变种及ZM-Seaside块根蛋白质双向电泳图谱（图 4-A和4-B）。以花叶变种为对照，经过Delta 2D软件分析块根的电泳图谱，得到平均差异表达量在2.0倍^[11]以上的蛋白质点39个（图 4-C），包括上调表达23个（黑色箭头所指），下调表达16个（白色箭头所指）。通过差异蛋白质点的质谱分析及NCBI数据库搜索比对，成功匹配到其中28个蛋白质（表 1）。

2.4 差异蛋白质互作调控网络的构建

将匹配得到的11个蛋白质通过String在线软件构建蛋白质互作调控网络，由图 5可知11个蛋白质节点共有47种蛋白质互作关系，其中互作最多的蛋白质是Heat shock protein（点2），它同9种蛋白质发生互作；其次是Molecular chaperone Hsp90-1（点1），含7种蛋白质互作关系，且这两种蛋白质均下调。再次是Heat-shock protein（点33，上调），含6种蛋白质互作关系，ATP synthase（点10，上调）、Small heat-shock protein（点26，下调）、18.1 kDa class I heat shock protein（点34，上调）均含4种蛋白质互作关系。Protein disulfide-isomerase（点6，下调）和Phosphoglycerate kinase，putative-R. communis（点11，上调）含3种蛋白质互作关系。Uroporphyrinogen decarboxylase，putative-R. communis（点15，上调）、Putative inorganic pyrophosphatase-Oryza sativa Japonica Group（点17，下调）、Predicted protein -P. trichocarpa（点19，下调）均有2种互作关系。

3 讨论

在前期研究中已经证明栽培种ZM-Seaside叶片的光合作用要显著优于花叶变种^[16]，而ZM-Seaside的块根产量、淀粉含量干物质率均显著高于花叶变种。由此可知叶片的光合作用能力能显著影响木薯块根的产量。前期也已经证明花叶变种叶片中与光合作用相关的蛋白质均下调^[16]，而本研究以花叶变种为对照，证明花叶变种块根中与淀粉积累相关的蛋白点也均下调。UGPase即尿苷二磷酸葡萄糖焦磷酸化酶，1953年由Munch-Petersen、Kalckar和Cutolo在酵母细胞中发现^{[17, 18]}。它催化反应UTP+葡萄糖-1-P生成UDP-葡萄糖+ PPi，而在高等植物中，产物UDPG被发现作为主要的葡萄糖基供体参与蔗糖、糖蛋白、纤维素等许多糖类代谢^[19]。吴晓俊等^[20]早已在黄芪中发现UGPase是合成多糖的关键酶，在黄芪中总多糖和可溶性多糖积累量均与UGPase的活性呈正相关，相关系数均高于0.9。UGPase被证明与AGPase相偶联形成ADPG，参与淀粉的合成^{[21, 22]}。而AGPase在植物淀粉合成中作为关键的限速酶，催化1-磷酸葡萄糖与ATP作用生成ADPG，为淀粉合成提供葡萄糖基，进而决定作物的产量^[23]，因而AGPase控制着淀粉的合成^[24]，是影响作物产量的关键酶。1955年SPS由Leloir和Cardini在小麦胚芽中发现^[25]。SPS是以UDPG为供体，以6-磷酸果糖（F-6-P）为受体的糖转移酶，6-磷酸蔗糖在磷酸蔗糖磷酸化酶的作用下脱磷酸并水解形成蔗糖和磷酸根离子^[26]。SPS是蔗糖合成过程中的限速酶^[27]。蔗糖是高等植物光合作用的主要产物，是碳运输的主要形式，也是“库”代谢的主要基质^[28]。Hubbard等^[29]证实网纹甜瓜果实蔗糖积累与SPS活性上升相关。这3种与产量相关的主要蛋白质在花叶变种中的下调与花叶变种的低产量结果一致。

以花叶变种为对照，比较栽培种ZM-Seaside块根的全蛋白质变化，成功鉴定出的39个差异蛋白质中，23个上调表达，16个下调表达。同时本研究还利用String在线软件构建差异蛋白质互作的生物调控网络，揭示了块根代谢的调控关系。在整个互作调控网络中热激蛋白Heat shock protein和分子伴侣Molecular chaperone Hsp90-1互作关系最多，这2种蛋白质是伴侣蛋白，伴侣蛋白参与到每一个生命活动中，因此伴侣蛋白的互作关系最多。Carvallo等^[30]曾报道他在研究类胡萝卜素的时候发现具有类胡萝卜素复合物和缺失类胡萝卜素复合物的两个品种中鉴定到大量差异伴侣蛋白。以及吕亚等^[31]也曾报道研究木薯叶片与光合作用日变化相关的差异蛋白的时候鉴定到大量伴侣蛋白，指出伴侣蛋白参与到各个生命活动中，在各种生命活动中起辅助作用。辅助作用的伴侣蛋白互作关系较多，其次是能量代谢ATP synthase（点10）互作关系最多。ATP用于能量供应，在栽培种ZM-Seaside块根中ATP合成酶的上调表明ZM-Seaside块根整个生命代谢活动旺盛，能量供应充足。因而其产量较花叶变种也高。本研究从蛋白质互作水平揭示了花叶变种块根产量低于ZM-Seaside的分子机理。这几个关键蛋白质有可能成为筛选木薯高产种质的标记蛋白质。后续将进一步通过酵母双杂等技术验证对这几个蛋白的互作水平进行验证。

4 结论

本研究通过直接田间测产测定高产种质和低产种质木薯块根的产量，确定两种木薯块根产量有较大差异，进而从全蛋白质水平揭示高产栽培种质ZM-Seaside块根产量显著高于低产种质花叶变种的分子机理。通过蛋白质互作调控网络推测淀粉合成相关蛋白质参与多个代谢途径，它们使各个代谢通路的蛋白质紧密相连，构成一个相互作用的生物调控网络。本研究通过蛋白质互作调控网络筛选出的关键蛋白质有可能成为选育丰产木薯品种的标记蛋白质。