基于平行反应监测验证技术分析绵羊胚胎骨骼肌蛋白质组数据

doi:10.11843/j.issn.0366-6964.2020.07.011

摘要/Abstract

摘要： 旨在对绵羊胚胎骨骼肌蛋白质组学数据进行分析，为阐明绵羊胚胎骨骼肌生长发育机制奠定基础。前期试验应用串联质谱定量法（tandem mass tag，TMT）对成年中国美利奴绵羊妊娠第85（D85）、105（D105）和135天（D135）的胚胎背最长肌进行蛋白质定量，并对D85/D105、D105/D135和D85/D135 3个比较组进行分析。本试验在此基础上，利用KEGG和平行反应监测（parallel reaction monitoring，PRM）等数据分析方法对3个比较组绵羊胚胎骨骼肌上、下调差异丰度蛋白质进行分析和验证，并对候选蛋白质进行生物信息学分析。本研究对肌纤维成熟分化标志性蛋白肌球蛋白重链（myosin-2 isoform X2，MYH）进行PRM靶向定量验证，结果表明，其变化趋势与绵羊胚胎肌纤维成熟分化趋势相符。本研究对上调差异丰度蛋白和下调差异丰度蛋白的KEGG分析结果表明，D105/D85比较组中上调差异丰度蛋白质显著富集于胰岛素等信号通路，D135/D105和D135/D85比较组中下调差异丰度蛋白质显著富集于DNA复制和蛋白质消化吸收等信号通路。与肌肉发育相关的cAMP依赖性蛋白激酶A催化亚基α（PRKACA）和葡萄糖转运蛋白成员4（GLUT4）蛋白均显著富集于胰岛素信号通路中，生物信息学分析显示，这两个蛋白质的理论分子量分别为121.73和20.64 ku，分别有147、12个带正电荷的氨基酸残基和135、12个带负电荷的氨基酸残基；理论等电点分别为8.81和6.54；亲水性平均系数分别为-0.408和0.811；PRKACA蛋白无N端糖基化位点，有45个磷酸化位点；GLUT4蛋白有1个N端糖基化位点，25个磷酸化位点；PRKACA蛋白三级结构与CAMP-2依赖蛋白激酶a嵌合融合的晶体结构相似度为85%，GLUT4蛋白三级结构与人葡萄糖转运蛋白GLUT1的相似性为78%。本研究表明，PRM定量验证趋势与绵羊胚胎肌纤维成熟分化趋势相符，PRKACA和GLUT4蛋白具有丰富的磷酸化修饰位点，并参与胰岛素信号通路调控，是重要的候选蛋白。

关键词: 绵羊, 胚胎骨骼肌, 平行反应监测, 蛋白质组学

Abstract: The purpose of this study was to analyze the data of embryonic skeletal muscle proteomic in sheep, and provide a basis for revealing the mechanism of embryonic skeletal muscle development and growth in sheep. In previous proteomic study, the embryonic longissimus dorsi proteins were quantified by using tandem mass tag (TMT) in the adult Chinese merino ewes on the Day 85 (D85), 105 (D105) and 135 (D135) of gestation, and 3 comparable groups D85/D105, D105/D135 and D85/D135 were set up. The KEGG and parallel reaction monitoring (PRM) methods were used to analyze and verify the up-and down-regulation differential abundance proteins in sheep embryonic skeletal muscle in 3 comparable groups, and bioinformatics analysis was performed for the candidate proteins. The marker protein about muscle fibers maturation and differentiation of myosin-2 isoform X2 (MYH) was verified by PRM, the results demonstrated that the change trend of MYH was consistent with the trend of sheep embryonic muscle fiber maturation and differentiation. The up-regulation and down-regulation differential abundance proteins were analyzed by KEGG method, the results showed that the up-regulation differen-tial abundance proteins in D105/D85 comparable group were significantly enriched in insulin signaling pathway, and the down-regulation differential abundance proteins in D135/D105 and D135/D85 comparable groups were significantly enriched in DNA replication and protein digestion and absorption signaling pathways. The muscle development related proteins, cAMP-dependent protein kinase catalytic subunit alpha (PRKACA) and glucose transporter member 4 isoform X1 (GLUT4), were significantly enriched in insulin signaling pathway. The bioinformatics analysis of the two proteins revealed that the theoretical molecular weight of PRKACA and GLUT4 were 121.73 and 20.64 ku, respectively. The PRKACA contained 147 positively charged amino acid residues and 135 negatively charged amino acid residues, the theoretical isoelectric point was 8.81, the hydrophilic average coefficient was -0.408; and the GLUT4 contained 12 positively charged amino acid residues and 12 negatively charged amino acid residues, the theoretical isoelectric point was 6.54, the hydrophilic average coefficient was 0.811. There were no N-terminal glycosylation site and there was 45 phosphorylation sites in PRKACA; there was an N-terminal glycosylation site and 25 phosphorylation sites in GLUT4. The homology of PRKACA and CAMP-2 relies on the chimeric fusion crystal structure of protein kinase a was 85% in tertiary structure, the homology of the GLUT4 and human glucose transporter GLUT1 was 78% in tertiary structure. The results showed that the trend of MYH quantification verified by PRM was consistent with the trend of embryonic muscle fiber maturation and differentiation. PRKACA and GLUT4 were important candidate proteins, with abundant phosphorylation sites, involved in the regulating of insulin signaling pathway.

Key words: sheep, embryonic skeletal muscle, parallel reaction monitoring, proteomics

中图分类号:

S826.2

王欣悦, 赵志达, 石田培, 尚明玉, 张莉. 基于平行反应监测验证技术分析绵羊胚胎骨骼肌蛋白质组数据[J]. 畜牧兽医学报, 2020, 51(7): 1587-1596.

WANG Xinyue, ZHAO Zhida, SHI Tianpei, SHANG Mingyu, ZHANG Li. The Data Analysis of Embryonic Skeletal Muscle Proteomic in Sheep Based on Parallel Reaction Monitoring Technology[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(7): 1587-1596.

参考文献 42

[1]	WHITE R B,BIÉRINX A S,GNOCCHI V F,et al.Dynamics of muscle fibre growth during postnatal mouse development[J].BMC Dev Biol,2010,10(1):21.
[2]	BENTZINGER C F,WANG Y X,RUDNICKI M A.Building muscle:molecular regulation of myogenesis[J]. Cold Spring Harb Perspect Biol,2012,4(2):a008342.
[3]	ASHMORE C R,ROBINSON D W,RATTRAY P,et al.Biphasic development of muscle fibers in the fetal lamb[J].Exp Neurol,1972,37(2):241-255.
[4]	HINDI S M,TAJRISHI M M,KUMAR A.Signaling mechanisms in mammalian myoblast fusion[J].Sci Signal,2013,6(272):re2.
[5]	MURPHY M,KARDON G.Origin of vertebrate limb muscle:the role of progenitor and myoblast populations[J]. Curr Top Dev Biol,2011,96:1-32.
[6]	GALLIEN S,BOURMAUD A,KIM S Y,et al.Technical considerations for large-scale parallel reaction monitoring analysis[J].J Proteomics,2014,100:147-159.
[7]	ALEKSANDR A D R G.Targeted proteomics:bench to bedside[J].J Proteomics,2018,189:3-4.
[8]	DOERR A.Targeting with PRM[J].Nat Methods,2012,9(10):950.
[9]	PICOTTI P,AEBERSOLD R.Selected reaction monitoring-based proteomics:workflows,potential,pitfalls and future directions[J].Nat Methods,2012,9(6):555-566.
[10]	WANG X Y,SHI T P,ZHAO Z D,et al.Proteomic analyses of sheep (ovis aries) embryonic skeletal muscle[J].Sci Rep,2020,10(1):1750.
[11]	COX J,NEUHAUSER N,MICHALSKI A,et al.Andromeda:a peptide search engine integrated into the MaxQuant environment[J].J Proteome Res,2011,10(4):1794-1805.
[12]	李琳,王磊,王文静,等.基于串联质谱标签法和平行反应监测技术的氟喹诺酮耐药沙门菌蛋白质组学分[J]. 微生物学通报,2018,45(7):1535-1545.LI L,WANG L,WANG W J,et al.Fluoroquinolone resistant salmonella proteomics analysis based on tandem mass tag and parallel reaction monitoring techniques[J].Microbiology China,2018,45(7):1535-1545.(in Chinese)
[13]	GASTEIGER E,HOOGLAND C,GATTIKER A,et al.Protein identification and analysis tools on the ExPASy server[M]//WALKER J M.The Proteomics Protocols Handbook.Totowa:Humana Press,2005.
[14]	BLOM N,SICHERITZ-PONTÉN T,GUPTA R,et al.Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence[J].Proteomics,2004,4(6):1633-1649.
[15]	STEENTOFT C,VAKHRUSHEV S Y,JOSHI H J,et al.Precision mapping of the human O-GalNAc glycoproteome through simple cell technology[J].EMBO J,2013,32(10):1478-1488.
[16]	BLOM N,GAMMELTOFT S,BRUNAK S.Sequence and structure-based prediction of eukaryotic protein phosphorylation sites[J].J Mole Biol,1999,294(5):1351-1362.
[17]	BUCKINGHAM M E.Muscle:the regulation of myogenesis[J].Curr Opin Genet Dev,1994,4(5):745-751.
[18]	MARX V.Targeted proteomics[J].Nat Methods,2013,10(1):19-22.
[19]	URISMAN A,LEVIN R S,GORDAN J D,et al.An optimized chromatographic strategy for multiplexing in parallel reaction monitoring mass spectrometry:insights from quantitation of activated kinases[J].Mol Cell Proteomics,2017,16(2):265-277.
[20]	PETERSON A C,RUSSELL J D,BAILEY D J,et al.Parallel reaction monitoring for high resolution and high mass accuracy quantitative,targeted proteomics[J].Mol Cell Proteomics,2012,11(11):1475-1488.
[21]	王素兰,高华萍,张菁,等.基于稳定同位素标记和平行反应监测的蛋白质组学定量技术用于肝癌生物标志物的筛选和验证[J].色谱,2017,35(9):934-940.WANG S L,GAO H P,ZHANG J,et al.Stable isotope labeling and parallel reaction monitoring-based proteomic quantification for biomarker screening and validation of hepatocellular carcinoma[J].Chinese Journal of Chromatography,2017,35(9):934-940.(in Chinese)
[22]	TIMSON D J,TRAYER H R,SMITH K J,et al.Size and charge requirements for kinetic modulation and actin binding by alkali 1-type myosin essential light chains[J].J Biol Chem,1999,274(26):18271-18277.
[23]	HAMELIN M,SAYD T,CHAMBON C,BOUIX J,et al.Proteomic analysis of ovine muscle hypertrophy[J].J Anim Sci,2006,84(12):3266-3276.
[24]	吴智春,王浩,于华芸.与细胞增殖有关的胰岛素信号通路[J].中国老年学杂志,2009,29(15):1988-1990.WU Z H,WANG H,YU H Y.Insulin signaling pathways associated with cell proliferation[J].Chinese Journal of Gerontology,2009,29(15):1988-1990.(in Chinese)
[25]	LOWENSTEIN E J,DALY R J,BATZER A G,et al.The SH2 and SH3 domain-containing protein GRB2 links receptor tyrosine kinases to ras signaling[J].Cell,1992,70(3):431-442.
[26]	SKOLNIK E Y,LEE C H,BATZER A,et al.The SH2/SH3 domain-containing protein GRB2 interacts with tyrosine-phosphorylated IRS1 and Shc:implications for insulin control of ras signalling[J].EMBO J,1993, 12(5):1929-1936.
[27]	张勇,马勇,朱宇旌,等.p38丝裂原活化蛋白激酶信号途径调节骨骼肌生长发育的机理[J].动物营养学报,2012,24(1):14-19.ZHANG Y,MA Y,ZHU Y J,et al.Mechanism of p38 MAPK signalingpathway in regulating skeletal muscle growth and development[J].Chinese Journal of Animal Nutrition,2012,24(1):14-19.(in Chinese)
[28]	RAMOCKI M B,JOHNSON S E,WHITE M A,et al.Signaling through mitogen-activated protein kinaseand Rac/Rho does not duplicate the effects of activated Ras on skeletal myogenesis[J].Mol Cell Biol,1997, 17(7):3547-3555.
[29]	PAGE J L,WANG X,SORDILLO L M,et al.MEKK1signaling through p38 leads to transcriptionalinactivation of E47 and repression of skeletal myo-genesis[J].J Biol Chem,2004,279(30):30966-30972.
[30]	TOBE K,ASAI S,MATUOKA K,et al.Cytoskeletal reorganization induced by insulin:involvement of Grb2/Ash, Ras and phosphatidylinositol 3-kinase signalling[J].Genes Cells,2003,8(1):29-40.
[31]	JENSEN J,BRENNESVIK E O,LAI Y C,et al.GSK-3β regulation in skeletal muscles by adrenaline and insulin:evidence that PKA and PKB regulate different pools of GSK-3[J].Cell Signal,2007,19(1):204-210.
[32]	AHMAD M,FLATT P R,FURMAN B L,et al.The role of the cyclic GMP-inhibited cyclic AMP-specific phosphodiesterase (PDE3) in regulating clonal BRIN-BD11 insulin secreting cell survival[J].Cell Signal, 2000, 12(8):541-548.
[33]	VANDER HAAR E,LEE S I,BANDHAKAVI S,et al.Insulin signalling to mTOR mediated by the Akt/PKB substrate PRAS40[J].Nat Cell Biol,2007,9(3):316-323.
[34]	张小宁,陆健,任航行,等.Texel和乌珠穆沁绵羊妊娠中后期胎儿背最长肌中AKT基因的发育表达模式[J].中国畜牧兽医,2012,39(8):20-25.ZHANG X N,LU J,REN H X,et al.Developmental expression patterns of AKT gene in fetal longissimus dorsi muscle between texel and ujumqin sheep during the second half of gestation[J].China Animal Husbandry & Veterinary Medicine,2012,39(8):20-25.(in Chinese)
[35]	TURNHAM R E,SCOTT J D.Protein kinase a catalytic subunit isoform PRKACA;history,function and physiology[J].Gene,2016,577(2):101-108.
[36]	SKÅLHEGG B S,HUANG Y Z,SU T,et al.Mutation of the calpha subunit of PKA leads to growth retardation and spermdysfunction[J].Mol Endocrinol,2002,16(3):630-639.
[37]	NIU W Y,HUANG C,NAWAZ Z,et al.Maturation of the regulation of GLUT4 activity by p38 MAPK during L6 cell myogenesis[J].J Biol Chem,2003,278(20):17953-17962.
[38]	MORIYA N,MIYAZAKI M.Akt1 deficiency diminishes skeletal muscle hypertrophy by reducing satellite cell proliferation[J].Am J Physiol Regul Integr Comp Physiol,2018,314(5):R741-R751.
[39]	郑琪,睢梦华,凌英会.骨骼肌卫星细胞增殖与成肌分化过程中关键信号通路的作用[J].畜牧兽医学报,2007,48(11):2005-2014.ZHENG Q,SUI M H,LING Y H.The role of key signalingpathways in the proliferation and differentiation of skeletal muscle satellite cells[J].Acta Veterinaria et Zootechnica Sinica,2007,48(11):2005-2014.(in Chinese)
[40]	LIU S Q,HAN W P,JIANG S Y,et al.Integrativetranscriptomics and proteomics analysis of longissimus dorsi muscles of Canadian double-muscled large white pigs[J].Gene,2016,577(1):14-23.
[41]	辛向博.Myostatin基因对骨骼肌代谢机制影响的定量蛋白组学与磷酸化蛋白组学研究[D].天津:天津农学院,2018.XIN X B.Quantitative proteomics andphosphoproteomics study on the effect of Myostatin gene on skeletal muscle metabolism[D].Tianjin:Tianjin Agricultural University,2018.(in Chinese)
[42]	YANG S P,LI X F,LIU X F,et al.Parallel comparative proteomics and phosphoproteomics reveal that cattle myostatin regulates phosphorylation of key enzymesin glycogen metabolism and glycolysis pathway[J]. Oncotarget,2018,9(13):11352-11370.