色谱  2014, Vol. 32 Issue (3): 284-289   PDF (1084KB)    
扩展功能
加入收藏夹
复制引文信息
加入引用管理器
Email Alert
RSS
本文作者相关文章
郭静
晏嘉泽
郭明
靳艳
反相液相色谱-串联质谱法鉴定油菜蜂花粉中的蛋白质及活性肽
郭静1,2, 晏嘉泽2, 郭明1, 靳艳2    
1. 大连大学环境与化学工程学院, 辽宁 大连 116622;
2. 中国科学院大连化学物理研究所, 中国科学院分离分析化学重点实验室, 辽宁 大连 116023
摘要:基于鸟枪法蛋白质组学分析方法,使用反相液相色谱-串联质谱(RPLC-MS/MS)系统分析油菜蜂花粉蛋白质的胰蛋白酶酶解产物,结合数据库检索,共鉴定到353条肽段。鉴定到的肽段所归属的蛋白质中有239个蛋白质可检索到其分子生物学功能,主要功能为结合活性、酶活性、运输活性、抑制活性等。根据血管紧张素转化酶(ACE)抑制肽活性与多肽构效之间的关系,从鉴定到的肽段中筛选并适当修饰后得到5条可能具有ACE抑制活性的肽段,化学合成肽段后进行了活性验证。结果表明5条肽段均具有良好的活性,其中肽段AELDIVLALF和LAVNLIPFP表现出较高的ACE抑制活性,半数抑制浓度(IC50)分别为(10.65±0.50)μmol/L和(23.66±1.08)μmol/L。该方法速度快,成本低,大大缩短了鉴定周期,达到了高通量筛选生物活性肽的目的。
关键词反相液相色谱-串联质谱     蛋白质     血管紧张素转化酶抑制肽     油菜蜂花粉    
Application of reversed-phase liquid chromatography-tandem mass spectrometry in the identification of protein and bioactivity peptides from rape bee pollen
GUO Jing1,2, YAN Jiaze2, GUO Ming1, JIN Yan2    
1. College of Environmental and Chemical Engineering, Dalian University, Dalian 116622, China;
2. Key Laboratory of Separation Science for Analytical Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
Abstract: Based on the shotgun proteomic method, rape bee pollen protein was prepared with ultrasonic extraction and digested by trypsin, then separated and sequenced by reversed-phase liquid chromatography-tandem mass spectrometry (RPLC-MS/MS), followed by protein database searching. After the above analysis, 353 peptides were identified and the molecular biological functions of 239 proteins have been known. The identified molecular biological functions of these proteins mainly included binding activity, enzyme activity, transporter activity, inhibitor activity and so on. Five peptides were obtained after the screening and appropriate amino acid modification among the identified 353 peptides, according to the relationship between the sequence structure of the bioactivity peptide and angiotensin converting enzyme (ACE) inhibitor activity. The five peptides were speculated to have ACE inhibitor activity and synthesized to detect ACE inhibitor activity. The results showed that all of the five peptides had good ACE inhibitor activity. The peptides of AELDIVLALF and LAVNLIPFP among the five peptides displayed especially good ACE inhibition with half maximal inhibitory concentration (IC50) of (10.65±0.50) μmol/L and (23.66±1.08) μmol/L, respectively. This method is rapid, low-cost and achieves the goal of high-throughput screening of bioactivity peptides that greatly shorten the period of identification compared with traditional methods.
Key words: reversed-phase liquid chromatography-tandem mass spectrometry (RPLC-MS/MS)     protein     angiotensin converting enzyme (ACE) inhibitory peptides     rape bee pollen    

蜂花粉含有蛋白质、糖类、脂类、维生素、微量元素等化学成分,被誉为“微型营养库”和“完全营养品”。蜂花粉作为药食同源物质在我国具有悠久的食用历史,近年来研究发现蜂花粉具有清除自由基、提高免疫力、抑制前列腺疾病、改善心血管疾病、抗肿瘤等功能[ 1 ]。邓建军等[ 2 ]分离提取了油菜蜂花粉中的不同蛋白质,并研究了其抗氧化活性。张红城等[ 3 ]对6种蜂花粉中不同酶的活性进行了研究。

以动植物蛋白质为原料,用酶解法制备具有生物活性的肽是生物活性肽的研究方向之一。蜂花粉中蛋白质含量丰富,可达7% ~30% ,为制备生物活性肽提供了丰富的资源。已有研究发现由蜂花粉制备的生物活性肽具有抗氧化[ 4 ]、抗辐射[ 5 ]、抗肿瘤[ 6 ]等活性。在生物活性肽的研究中,以抑制血管紧张素转化酶(ACE)为靶标的降血压肽尤为引人关注,最早的ACE抑制肽是从蛇毒中发现的[ 7 ]。酶法降解食物源蛋白质制备ACE抑制肽具有安全性高等特点,现已从大豆蛋白质、藻蛋白质、鱼蛋白质等蛋白质酶解产物中分离得到大量的ACE抑制肽[ 8,9,10 ]

目前常用的生物活性肽研究方法是通过色谱分离等手段得到单一肽段,再利用质谱分析等技术鉴定肽段氨基酸序列,从而获得混合肽段中活性肽段的信息。但是传统的分离方法需要经过繁琐的分离纯化步骤才可以得到极少的目标产物,耗时耗力、效率低下。高通量、低损耗是质谱技术的发展趋势,相关行业研究领域也对此需求迫切[ 11 ]。基于液相色谱-质谱技术的鸟枪法是蛋白质组学研究的策略之一,该方法通过直接对蛋白质混合物进行酶解,经液相色谱-质谱分离,鉴定肽段,继而推导出样品蛋白质[ 12 ]。于海洋等[ 13 ]基于质谱学技术对锦灯笼中的蛋白质进行了鉴定及功能归属,Liu[ 14 ]等利用鸟枪法的策略分析了鹿血浆中的蛋白质,并结合ACE抑制肽的构效关系分析得到几条活性较好的ACE抑制肽。本文以油菜蜂花粉为对象,利用鸟枪法蛋白质组学研究策略直接分析其蛋白质酶解产物,研究蛋白质组成,并在此基础上筛选具有ACE抑制活性的肽段,极大地缩短了活性肽的研究周期。

1 实验部分
1.1 仪器与试剂

RPLC-MS/MS系统由Finnigan Surveyor液相色谱泵、线性离子阱(LTQ)质谱组成(美国Thermo公司);紫外/可见分光光度计(日本JASCO公司);机械式搅拌器(德国IKA公司); Magic C18 AQ反相色谱填料(5 μ m,12 nm,美国Michrom BioResources公司);石英毛细管(75 μ m×120 mm,美国Polymicro Technologies公司); Milli-Q超纯水一体机(美国Millipore公司)。

胰蛋白酶、血管紧张素转化酶(EC3.4.14.1)、马尿酰组氨酰亮氨酸(HHL)、三氟乙酸(TFA)、牛血清白蛋白(BSA)、聚乙二醇辛基苯基醚(Triton X-100)、蛋白酶抑制剂(cocktail)均购自美国Sigma-Aldrich公司。二硫苏糖醇(DTT)、碘代乙酰胺(IAA)购自美国Bio-Rad公司。乙腈(ACN,HPLC级)购自德国Merck公司。乙二胺四乙酸(EDTA)、乙二醇双(2-氨基乙基醚)四乙酸(EGTA)、苯甲基磺酰氟(PMSF)购自美国Amresco公司。油菜蜂花粉由宽甸万钧蜂业有限公司提供。其他药品和试剂均为国产分析纯或色谱纯。

1.2 蛋白质含量的测定

委托中国科学院沈阳生态所农产品安全与环境质量检测中心根据国家标准GB/T5009.5-2003 进行测定。

1.3 样品前处理
1.3.1 蛋白质的提取

向10 mL 50 mmol/L Tris-HCl缓冲液(pH 7.4,含8 mol/L 尿素、1 mmol/L EDTA和1 mmol/L EGTA)中加入0.1 g DTT,使其终浓度为65 mmol/L;再加入100 μ L Triton X-100,冰浴5~10 min后,加入100 μ L cocktail和100 μ L 100 mmol/L 的PMSF。将1.5 g经液氮环境中研磨破壁的油菜蜂花粉分散其中,超声提取3次(每次程序为:400 W超声4 s,停4 s,80个循环)。于4 ℃ 25000 g离心力下离心1 h,吸取上清液,加入4倍体积的沉淀液(丙酮 ∶ 乙醇 ∶ 乙酸=50 ∶ 50 ∶ 0.1,v/v/v),沉淀过夜;于4 ℃ 25000 g离心力下离心30 min,弃去上清液;丙酮洗涤脱色,于4 ℃ 25000 g离心力下离心15 min,弃去上清液;丙酮反复洗涤直至上清无色,再用75%(v/v)乙醇洗去丙酮。于4 ℃ 25000 g离心力下离心10 min,弃去上清液,沉淀物冷冻干燥后于-20 ℃保存。

1.3.2 蛋白质的酶解

冻干蛋白质用100 mmol/L NH4HCO3缓冲液(pH7.4,含8 mol/L 尿素)复溶,Bradford法测定蛋白质浓度。取质量浓度为2 mg/mL 的蛋白质溶液1 mL,加入DTT使其终浓度为20 mmol/L,于60 ℃水浴1 h;加入IAA使其终浓度为40 mmol/L,室温避光放置40 min;加入100 mmol/L NH4HCO3,使尿素终浓度为1 mol/L,以酶 ∶ 蛋白质=1 ∶ 25(w/w)的比例加入胰蛋白酶,于37 ℃酶解20 h。

1.3.3 蛋白质酶解产物的净化

C18 SPE柱经3倍柱体积的乙腈活化,再用3倍柱体积的0.1%(v/v)三氟乙酸溶液洗去乙腈。蛋白质酶解产物用50%(v/v)三氟乙酸溶液调节pH值至2~3后过柱。经3倍柱体积的0.1%(v/v)三氟乙酸溶液脱盐后,用1.2 mL 80%(v/v)乙腈/0.1%(v/v)三氟乙酸溶液分3次洗脱,洗脱液冷冻干燥后于-20 ℃保存。

1.4 RPLC-MS/MS分析

将冻干的酶解样品用0.1%(v/v)甲酸溶液复溶至0.4 μ g/ μ L,取20 μ L进行RPLC-MS/MS分析,重复3次。

将毛细管(内径75 μ m)的一端拉成内径约为5 μ m的尖端,把Magic C18 AQ填料压入柱内,填充柱长度在12 cm左右,将毛细管柱与质谱相连。所用的流动相A为0.1%(v/v)的甲酸溶液,流动相B为0.1%(v/v)的甲酸/乙腈溶液,流速为200 nL/min。梯度洗脱过程如下:0~2 min,0%B~5%B;2~122 min,5%B~35%B;122~125 min,35%B~80%B;125~135 min,80%B;135~137 min,80%B~0%B;137~150 min,0%B。

设置离子传输毛细管的温度为200 ℃,电喷雾电压为1.8 kV,归一化碰撞能量为35.0% 。均使用数据依赖模式对MS和MS/MS进行图谱采集,扫描模式为全扫描。全扫描中丰度最高的6个离子峰进行MS/MS扫描,其中动态排除设置为:重复次数为2,重复容忍时间为30 s,动态排除时间为90 s。利用Xcalibur软件(Thermo)进行系统控制和数据收集。

1.5 数据库检索

采集的数据用SEQUEST软件分别在蜜蜂(Apis,蛋白质数目为 14295)和油菜(Brassica campestris L.,蛋白质数目为 41262)蛋白质数据库(http:/ /www. uniprot. org/)中进行检索。检索参数为:酶为胰蛋白酶,最高允许漏切位点为2,前体离子的质量容忍度为2 Da,碎片离子的质量容忍度为1 Da。搜库结果通过本实验室创建的Armone软件[ 15 ]进行优化,筛选标准为Xcorr (cross correlation value)不小于1.9、2.2和3.75(分别对应+1,+2和+3价离子),肽段假阳性率(FDR)控制在1%以内。

1.6 肽段合成及其ACE抑制活性的检测

委托杭州中肽公司(浙江杭州)合成肽段,经RPLC-MS检测纯度。将Cushman等[ 16 ]的方法改进后进行体外ACE抑制活性测定:将50 μ L 5 mmol/L HHL与50 μ L 50 mmol/L的硼酸盐缓冲液( pH 8.3,含300 mmol/L NaCl)混合均匀,加入50 μ L一定浓度的抑制剂溶液,于37 ℃水浴5 min;加入50 μ L 0.1 U/mL的ACE,于37 ℃水浴30 min,加入200 μ L 1 mol/L 的HCl终止反应;再加入1 mL乙酸乙酯振荡15 s,于3500 r/min 离心5 min;取0.8 mL上清液,于90 ℃水浴干燥15 min,复溶于0.8 mL水中,于228 nm处检测吸光值。

抑制率计算公式为:

其中R表示抑制率;A表示对照组(将抑制剂溶液用缓冲液替代)的吸光度值;A0表示空白对照组(将抑制剂溶液和ACE均用缓冲液替代)的吸光度值;B表示样品组的吸光度值;B0表示空白样品组(将样品组的ACE用缓冲液替代)的吸光度值。
2 结果与讨论
2.1 油菜蜂花粉蛋白质分析

油菜蜂花粉样品经测定,其蛋白质含量为29.1% 。油菜蜂花粉按照1.3节中方法处理后进行RPLC-MS/MS分析,基峰色谱图如图1所示。经检索,在蜜蜂(Apis)蛋白质库中共鉴定到62条肽段,在油菜(Brassica campestris L.)蛋白质库中共鉴定到291条肽段,油菜蜂花粉中的大部分蛋白质来源于油菜。不同物种中鉴定到的相同序列不进行重复计数,鉴定到肽段的来源分布见图2。这353条肽段所归属的蛋白质,有239个蛋白质在Universal Protein数据库中可以查找到其分子生物学功能,主要有结合活性(binding activity)、酶活性(enzyme activity)、运输活性(transporter activity)、结构组成(structural constituent)、催化活性(catalytic activity)、抑制活性(inhibitor activity)等。在所鉴定到的蛋白质中,有些同时具有多重分子生物学功能,如UvrABC system protein B同时具有催化活性和结合活性。

图 1 RPLC-MS/MS分析油菜蜂花粉蛋白质酶解产物的基峰色谱图 Fig. 1 Base peak chromatogram of protein degradation from rape bee pollen by RPLC-MS/MS analysis

图 2鉴定肽段的来源分布 Fig. 2Source distribution of the identified peptides

这些蛋白质按生物学功能分类情况如图3所示,图中数值代表具有该功能的蛋白质数量,其中61%的蛋白质具有结合活性或酶活性。结合活性包括:结合离子(如锌离子、钙离子、铁离子、镁离子等)的活性;结合脱氧核糖核酸(DNA)、核糖核酸(RNA)、三磷酸腺苷(ATP)等的活性;结合脂质、铁血红素、碳水化合物以及结合辅助因子的活性。酶活性包括氧化还原酶活性、水解酶活性、异构酶活性、合成酶活性等。

图 3油菜蜂花粉中鉴定到的蛋白质的生物学功能分类 Fig. 3Biological function distribution of the proteins identified in rape bee pollen
2.2 ACE抑制肽的筛选

ACE抑制肽的生物活性受肽段氨基酸组成、序列结构和长度的影响。大多数ACE抑制肽由2~12个氨基酸组成,C末端氨基酸多为芳香族氨基酸(如Y、 F)或疏水氨基酸(如L、 I)[ 17,18 ]。Cushman等[ 19 ]认为肽段C末端三肽强烈影响其ACE抑制活性,能够与ACE活性位点相互作用的最为理想的C末端三肽序列为FAP或者WAP,而W、Y、P、F则被认为是C末端最合适的氨基酸残基;同时当N末端为支链氨基酸如V、I时,肽段对ACE活性抑制较强。Meisel等[ 20 ]通过计算机分子建模对ACE抑制肽构效关系进行研究,发现抑制肽能与不被底物所占据的亚活性位点或不同于ACE催化部位的阴离子抑制剂结合位点相互作用,ACE抑制肽中相似的正电势几乎均位于C末端的相同区域。肽段中疏水性氨基酸的含量也会对其活性有较大的影响[ 21 ]。根据以上内容推测在本实验中鉴定到的353条肽段中,AGFAGDDAPR、AELDIVLALFK、LAVNLIPFPR、IIIALLLYK和LLICGGSAYPR这5条肽段可能具有ACE抑制活性。

C末端以P、F结束的肽段对ACE具有很强的抑制活性[ 22,23 ]。Ghassem等[ 24 ]从鱼类蛋白质中分离出以P结束的肽段NGTVVFEPP,该肽段具有较高的ACE抑制活性。Saiga等[ 25 ]从鸡胸肉中分离纯化出P4肽(GFHypGTHypGLHypGF),并对其氨基酸序列结构与ACE抑制活性之间的关系进行了研究,发现当移除C末端的F后,对ACE的抑制活性大大降低,半数抑制浓度(IC50)由46 μ mol/L 升至 25000 μ mol/L。鉴于以上报道,将筛选到的AELDIVLALFK和LAVNLIPFPR肽段C端的K、R去掉,序列变为AELDIVLALF和LAVNLIPFP,推测其会具有更高的ACE抑制活性。

2.3 ACE抑制活性测定

化学合成2.2节中推测具有ACE活性的肽段并测定其ACE抑制活性,结果见表1。其中AGFAGDDAPR来源于蜜蜂,其对ACE活性抑制的IC50值为(202.48±9.60) μ mol/L。该肽段C末端三肽中有A、P两个疏水性氨基酸,N末端为支链氨基酸A,符合ACE抑制肽的构象特征。其他4条肽段均来源于花粉,其中肽段AELDIVLALF对ACE活性的抑制效果最好,其IC50值为(10.65±0.50) μ mol/L。在ACE抑制肽数据库(http:/ /www.cftri.com/pepdb/index.php)中,已知抑制ACE活性IC50值的肽段有619条,将AELDIVLALF的IC50值与之比较后排序,结果位于第147位,抑制活性居数据库中前24% 。肽段AELDIVLALF的C末端三肽均为疏水性氨基酸,分别是A、L和F,C末端为F,且整条肽段含有较多的疏水性氨基酸,包括A、V、L、F和I。LAVNLIPFP的抑制活性位于数据库中前34% ,其IC50值为(23.66±1.08) μ mol/L。其C末端三肽含芳香族氨基酸F、疏水性氨基酸P,N末端为支链氨基酸L,且含有2/3的疏水性氨基酸,包括A、V、L、F、I和P,具有良好的ACE抑制肽的结构。LLICGGSAYPR和IIIALLLYK也表现出较好的ACE抑制活性,IC50值分别为(43.58±1.96) μ mol/L 和(28.61±1.37) μ mol/L,这两条肽段C末端三肽均含有两个疏水性氨基酸,N末端均为支链氨基酸,且疏水性氨基酸在整条肽链中比例均在50%以上。本方法基于蛋白质组学分析,结合构效关系筛选生物活性肽,通过化学合成肽段后进行生物活性验证,相比需要经繁复的分离过程才能得到极少产物的常规生物活性肽鉴定、纯化方法,大大缩短了鉴定周期。

表 1 源于油菜蜂花粉的ACE抑制肽 Table 1ACE inhibitory peptides from rape bee pollen

表2是上述ACE抑制肽所属蛋白质的分子生物学功能描述,5条肽段所属的蛋白质具有结合功能、酶催化活性以及结构组成等功能。细胞膜上Na+-K+-ATP酶和Ca2+-ATP酶活性降低引起的细胞内Ca2+浓度升高是高血压发病的关键机制[ 26 ],三磷酸鸟苷(GTP)在GTP酶的催化作用下先转化为二磷酸鸟苷(GDP),之后转化成ATP。GTP酶活性及GTP结合功能将直接或间接影响Na+-K+-ATP酶和Ca2+-ATP酶的活性,从而影响高血压发病。ABC(ATP-binding cassette transporter)家族是一组ATP结合转运蛋白质,主要从事细胞内外的物质转运,早期认为该家族主要与肿瘤细胞的多药耐药性有关,近年来发现其在各个器官及疾病(尤其是心血管疾病)中表现出多种其他重要病理和生理功能[ 27 ]。从筛选出的肽段所属蛋白质的生物功能上分析,这些肽段均与高血压存在一定的相关性。

表 2 ACE抑制肽所属蛋白质的分子生物学功能 Table 2 Molecular biological functions of the proteins that ACE inhibitory peptides belong to
3 结论

本文基于鸟枪法蛋白质组学策略,利用蛋白质组学分析平台对油菜蜂花粉蛋白质进行研究,共鉴定到353条肽段,其所属蛋白质具有结合活性、酶活性及运输活性等分子生物学功能。根据ACE抑制肽的构效关系进行筛选及修饰,得到AGFAGDD-APR、AELDIVLALF、LAVNLIPFP、IIIALLLYK和LLICGGSAYPR这5条可能具有ACE抑制活性的肽段。化学合上述成肽段并进行ACE抑制活性验证,结果显示这5条肽段均具有良好的活性,其中肽段AELDIVLALF和LAVNLIPFP表现出较高的抑制活性,IC50值分别为(10.65±0.50) μ mol/L 和(23.66±1.08) μ mol/L。以上结果表明,将蛋白质组学分析平台与活性肽构效相结合可用于快速、高效筛选ACE抑制肽。

参考文献
[1]Wang K F. Nutritional Value and Therapeutic of Pollen. Beijing: Chemical Industry Press (王开发. 花粉营养价值与食疗. 北京: 化学工业出版社), 2009
[2]Deng J J, Yang H X, Cao W, et al. Food and Fermentation Industries (邓建军, 杨海霞, 曹炜, 等. 食品与发酵工业), 2011, 37(10): 218
[3]Zhang H C, Cheng M, Dong J, et al. Food Science (张红城, 程蒙, 董捷, 等. 食品科学), 2009, 30(21): 229
[4]Zhu X L. [MS Dissertation]. Beijing: Chinese Academy of Agricultural Sciences (朱晓丽. [硕士学位论文]. 北京: 中国农业科学院), 2007
[5]Wang L, Wen C R, Wang Y E, et al. Acta Nutrimenta Sinica (汪磊, 温成荣, 王玉娥, 等. 营养学报), 2010, 32(1): 72
[6]Hu X B. [PHD Dissertation]. Wuhan: Huazhong Agriculture University (胡筱波. [博士学位论文]. 武汉: 华中农业大学), 2007
[7]Ferreira S H, Bartelt D C, Greene L J. Biochemistry, 1970, 9(13): 2583  
[8]Nakahara T, Sano A, Yamaguchi H, et al. J Agric Food Chem, 2010, 58(2): 821  
[9]Sano M, Hosokawa T, Yamaguchi T, et al. J Agric Food Chem, 2002, 50(21): 6245  
[10]Wijesekara I, Qian Z J, Ryu B M, et al. Food Res Int, 2011, 44(3): 703  
[11]Lin W X, Sun X Q, Ma J. Chinese Journal of Chromatography (林维宣, 孙兴权, 马杰. 色谱), 2013, 31(5): 410
[12]Han G H, Ye M L, Jiang X G, et al. Proteomics, 2008, 8(7): 1346  
[13]Yu H Y, Yan J Z, Guo M, et al. Chinese Journal of Chromatography (于海洋, 晏嘉泽, 郭明, 等. 色谱), 2013, 31(4): 362
[14]Liu X Y, Song C X, Chen R, et al. Chin J Chem, 2010, 28: 1665  
[15]Jiang X, Ye M, Han G, et al. Anal Chem, 2010, 82(14): 6168  
[16]Cushman D W, Cheung H S. Biochem Pharmacol, 1971, 20: 1637  
[17]Nakamura Y, Yamamoto N, Sakai K, et al. J Dairy Sci, 1995, 78(6): 1253  
[18]López-Fandio R, Otte J, Van Camp J. Int Dairy J, 2006, 16(11): 1277  
[19]Cushman D W, Ondetti M A. Nat Med, 1999, 5(10): 1110  
[20]Meisel H. Livest Prod Sci, 1997, 50(1/2): 125
[21]Kohmura M, Nio N, Ariyoshi Y. Agri Biol Chem, 1990, 54(3): 835  
[22]Nakamura Y, Yamamoto N, Sakai K, et al. Dairy Sci, 1995, 78(4): 777  
[23]Vercruysse L, Van Camp J, Smagghe G. J Agric Food Chem, 2005, 53(21): 8106  
[24]Ghassem M, Arihara K, Babji A S, et al. Food Chem, 2011, 129(4): 1770  
[25]Saiga A, Okumura T, Makihara T, et al. J Agric Food Chem, 2006, 54(3): 942  
[26]Fu Y, Wang S, Lu Z, et al. Chinese Medical Journal, 1998, 111(2): 147
[27]Li G P, Ye L. Chinese Journal of General Surgery (李高鹏, 叶露. 中国普通外科杂志), 2007, 16(6): 581