基于DIA技术探究和分析不同民族人乳中差异蛋白的方法与流程

基于dia技术探究和分析不同民族人乳中差异蛋白的方法
技术领域
1.本发明属于乳品技术领域，具体涉及一种基于数据非依赖采集的蛋白组学技术(dia)对采集的不同民族的人乳进行蛋白组成和含量差异性比较分析的方法。


背景技术：

2.母乳是婴儿最天然、最理想的食物，为新生儿提供生长发育所必需的营养，并增强抵抗力。在这种情况下，作为母乳替代物，婴幼儿配方粉的科学设计和研发具有重要的意义。
3.为使婴幼儿能够获得更接近母乳的营养成分，国内外相关研究机构长期深入探索母乳组分，试图找到最符合母乳的营养均衡型婴幼儿配方乳粉。我国在母乳研究营养成分上取得了一定的研究成果，相关机构纷纷建立了中国母乳营养成分数据库，但仍存在样本地域碎片化、样本分析常规化的不足。我国人口众多、地域辽阔、膳食多样，如何利用现代技术详尽地分析中国母乳功能组分，将常规母乳营养素研究扩展到更加细微的母乳成分研究显得至关重要。我国东北朝鲜族、西北地区等母乳研究较少，对母乳的蛋白质组学研究也较少。而相关研究表明，母乳成分受种族、饮食等多种因素影响。因此，探究不同民族人乳之间的营养成分的差异性，并开发针对性的婴配粉具有重要意义。
4.乳主要包括乳蛋白、乳糖和乳脂等多种营养成分，其中乳蛋白质是最重要的组成成分。乳蛋白质可以为新儿和哺乳动物提供氮、酶、激素及免疫物质等，对其健康成长有重要的作用。乳中含有种类多样、含量丰富的蛋白质，由于遗传变异和翻译后修饰对乳蛋白的影响，使乳蛋白的组成变得复杂、庞大。其中酪蛋白(cn)、乳清蛋白(wp)和乳脂球膜蛋白(mfgm)是乳蛋白的重要组成部分。蛋白质组学是研究蛋白组的一个新的领域，主要研究蛋白质的特性，包括蛋白质表达水平、氨基酸序列、翻译后加工和蛋白质相互作用，在蛋白质水平上了解细胞的各项功能、各种生理生化过程以及疾病的病理过程等。
5.采用蛋白质组学等技术对中国东北延边地区的朝鲜族和汉族人乳进行详细的分析和比较，探索民族差异对母乳中营养成分的影响，建立特色人群母乳数据库，为开发更加接近特定人群母乳营养成分及功效的婴幼儿配方奶粉提供依据，未来将有更多的婴幼儿受惠于针对性母乳数据库。


技术实现要素：

6.本发明所要解决的技术问题是针对蛋白种类繁多的我国不同民族的母乳，提供一种通过对样品进行高速分离得不同的乳清蛋白层、酪蛋白层和乳脂肪球膜蛋白层，并采用dia技术对蛋白进行高通量分析的方法。为了达到上述技术效果，本发明包括以下技术方案：
7.一种基于dia技术探究和分析不同民族人乳中差异蛋白的方法，有以下步骤：
8.步骤1：样品采集
9.采集不同民族新鲜母乳，每例30毫升放入50毫升的无菌离心管中，并贴上标签，获
得的样品立即放置在-80℃的冰箱中储存，同一民族随机6～10个样品混合成一个混合样品，每个民族收集至少三个平行的混合样本进行生物学重复；
10.步骤2：乳清蛋白提取和分离
11.将蛋白酶抑制剂加入母乳样品，取12ml母乳样品装入15ml离心管；4℃，3000rpm，离心30min分层，最上层为脂肪层，贴挂侧管壁，轻柔倒出下层脱脂奶液至新管，留脂肪层于原离心管中；磷酸缓冲液(pbs)冲洗脂肪层5次，4℃存放待用；脱脂奶液取1ml用10％乙酸调节溶液ph值至4.6，用于酪蛋白析出，4℃，10000rpm，离心15min，使乳清蛋白与酪蛋白分离；
12.乳清蛋白沉淀：在乳清蛋白与酪蛋白分离后，取上清液，加入4倍体积预冷丙酮，-20℃过夜，沉淀用预冷丙酮清洗3次，得到乳清蛋白沉淀，吹干备用；
13.酪蛋白沉淀：在乳清蛋白与酪蛋白分离后，除去上清液，沉淀即为酪蛋白，预冷丙酮清洗3次，吹干备用；
14.脂肪蛋白提取：在母乳样品进行离心分离后得到的脂肪层中加入3倍体积甲醇，涡旋震荡，至脂肪层破碎悬浮于甲醇，4℃，9000g，离心10s；除去上清液，加入体系2倍体积的氯仿，涡旋震荡，至沉淀悬浮于氯仿；向体系中加入体系3倍体积去离子水，混匀，4℃，9000g，离心1min，溶液分三层，下层为氯仿，中间层为蛋白层，将上清液移除，向剩余体系中加入3倍体积甲醇，混匀，4℃，9000g，离心2min；除上清，将沉淀用甲醇洗涤2次，吹干，沉淀即为脂肪蛋白；
15.步骤3：过滤辅助样品制备(fasp)酶解
16.基于过滤辅助样品制备(fasp)方法，对100μg蛋白质样品进行蛋白质酶水解，所述的蛋白质样品是步骤2中分离出的乳清蛋白、酪蛋白、脂肪蛋白的一种，向蛋白质样品中添加8m尿素，使总体积达到200μl，加入二硫苏糖醇(dtt)至体系中dtt最终浓度为10mm，并使混合物在56℃下静置30分钟，添加碘乙酸(iaa)至碘乙酸的终浓度为50mm，并在室温下无光反应40min，将样品转移到10k超滤管中，然后在室温下以12000g离心，弃掉滤液，并重复该步骤三次，添加200μl浓度为50mm的碳酸氢铵溶液，然后在室温下以12000g离心，丢弃滤液，并重复该步骤3次，将胰蛋白酶以样品:酶＝50:1的质量比添加到试管中，并在37℃下酶解16小时，通过在12000g、4℃下离心获得多肽样品，然后冷冻干燥；
17.步骤4：建库样本准备
18.酶解后的多肽样品在色谱仪中进行high ph反相色谱分级，柱子在初始条件下平衡5min，柱流速维持0.4ml/min；随后柱流速增至0.7ml/min，使用1.5ml ep管，每1分钟接收1管；
19.步骤5：dda建库数据采集
20.另取酶解后的多肽样品经由轨道阱质谱仪进行分析；共上样18μl，以60/120min梯度分离：3％溶液b至100％溶液b；柱流量控制在600nl/min；所述溶液b为98％acn，2％ddh2o，ph10；
21.轨道阱质谱仪在数据依赖采集模式下运行，质谱参数设置如下：(1)msn1:扫描范围(m/z):375-1500；分辨率：120,000；agctarget:4e5；最大注入时间：50ms；(2)msn2：分辨率：15,000；agctarget＝5e4；最大注入时间：22ms；碰撞能量：30％；将获得的质谱数据用proteome discoverer 2.1.0182进行蛋白质分析，最终获得dda参考谱图库；数据处理参数为：搜库引擎采用sequestht；蛋白质数据库用uniprottaxid：9606；设置胰蛋白酶酶解；最
大允许漏切位点设为2；肽段水平的假阳性率(fdr)设置为1％；肽和片段质量误差分别设置为
±
10ppm和
±
0.02da；
22.步骤6：dia进行生物样本定量数据采集
23.另取步骤3中的多肽样品经由轨道阱质谱仪进行分析；共上样18μl，以60/120min梯度分离：3％溶液b至100％溶液b；柱流量控制在600nl/min；
24.轨道阱质谱仪在数据依赖采集模式下运行，质谱参数设置如下：(1)msn1:扫描范围(m/z):350-1500；分辨率：60,000；agctarget:4e5；最大注入时间：50ms；(2)msn2：分辨率：30,000；agctarget＝3e5；最大注入时间：54ms；碰撞能量：33％；
25.用skyline软件进行dia分析，通过导入步骤5所得的dda参考谱图库，对dia原始数据进行多肽信号提取与定量分析，设置参数为：多肽长度选6-25；子离子m/z要求为大于母离子且last ion-3；子离子最大数目设为5；子离子最小数目设为2；子离子抽提窗口采用5min；dotp需要大于等于0.6；
26.步骤7：差异表达蛋白质筛选
27.将符合fc》1.5或《0.67，且pvalue《0.05的蛋白设定为差异表达蛋白，并进行go和kegg功能富集分析和蛋白互作分析。
28.有益效果：
29.1、本发明首次从不同民族的人乳差异性出发，系统化研究其蛋白组成和含量的差异性。
30.2、本发明采用高速分离的方法得到乳清蛋白、乳脂肪球膜蛋白以及酪蛋白三层。
31.3、本发明采用dia技术对朝鲜族和汉族母乳中的乳蛋白进行系统化、高通量的研究，并进行比较分析，鉴定出的总蛋白种类多，能全面的获得不同民族之间人乳蛋白的差异，为开发针对我国特定人群的创新性婴幼儿配方食品提供技术支持和新的研发思路。
附图说明
32.图1是乳清蛋白总鉴定蛋白及差异蛋白。
33.图2是乳清蛋白中差异蛋白的go富集分析。
34.图3是乳清蛋白中差异蛋白的kegg富集分析。
35.图4是酪蛋白总鉴定蛋白及差异蛋白。
36.图5是酪蛋白中差异蛋白的go富集分析。
37.图6是酪蛋白中差异蛋白的kegg富集分析。
38.图7是乳脂肪球膜蛋白总鉴定蛋白及差异蛋白。
39.图8是乳脂肪球膜蛋白中差异蛋白的go富集分析。
40.图9是乳脂肪球膜蛋白中差异蛋白的kegg富集分析。
具体实施方式
41.下面将通过具体的实施例对本发明做进一步的详细描述。
42.实施例1：样品采集及蛋白质分离
43.步骤1：样品采集
44.根据初步调查问卷，从吉林省延边自治州选择18名朝鲜族哺乳期母亲，在采集样
本前仔细清洁乳晕和乳头数次，然后用无菌水冲洗相关部位。轻轻地挤出母乳，丢弃第一滴，然后将大约30毫升的母乳放入50毫升的无菌离心管中，并贴上标签。获得的样品立即放置在-80℃的冰箱中，然后通过冷链运输运送到实验室。每6个样品混合成一个混合样，得到3个平行样本；用于后续做三个平行样的生物学重复测定；
45.用同样的步骤得到汉族母乳的3个平行样本。
46.将蛋白酶抑制剂加入母乳样品，取12ml母乳样品装入15ml离心管；4℃，3000rpm，离心30min分层，最上层为脂肪层，贴挂侧管壁，轻柔倒出下层脱脂奶液至新管，留脂肪层于原离心管中；磷酸缓冲液(pbs)冲洗脂肪层5次，4℃存放待用；脱脂奶液取1ml用10％乙酸调节溶液ph值至4.6，用于酪蛋白析出，4℃，10000rpm，离心15min，使乳清蛋白与酪蛋白分离；
47.乳清蛋白沉淀：在乳清蛋白与酪蛋白分离后，取上清液，加入4倍体积预冷丙酮，-20℃过夜，沉淀用预冷丙酮清洗3次，得到乳清蛋白沉淀，吹干备用；
48.酪蛋白沉淀：在乳清蛋白与酪蛋白分离后，除去上清液，沉淀即为酪蛋白，预冷丙酮清洗3次，吹干备用；
49.脂肪蛋白提取：在母乳样品进行离心分离后得到的脂肪层中加入3倍体积甲醇，涡旋震荡，至脂肪层破碎悬浮于甲醇，4℃，9000g，离心10s；除去上清液，加入体系2倍体积的氯仿，涡旋震荡，至沉淀悬浮于氯仿；向体系中加入体系3倍体积去离子水，混匀，4℃，9000g，离心1min，溶液分三层，下层为氯仿，中间层为蛋白层，将上清液移除，向剩余体系中加入3倍体积甲醇，混匀，4℃，9000g，离心2min；除上清，将沉淀用甲醇洗涤2次，吹干，沉淀即为脂肪蛋白；
50.实施例2乳清蛋白的差异性分析
51.基于过滤辅助样品制备(fasp)方法，对实施例1得到的乳清蛋白质样品(100μg)进行蛋白质酶水解，向样品中添加8m尿素，使其体积达到200μl。引入二硫苏糖醇(dtt)至最终浓度为10mm，并使混合物在56℃下静置30分钟。将碘乙酸(iaa)添加到50mm的最终水平，并在室温下无光反应40min。将样品转移到10k超滤管中，然后在室温下以12000g离心，弃掉滤液，并重复该步骤三次。添加碳酸氢铵溶液(200μl，50mm)，然后在室温下以12000g离心，丢弃滤液，并重复该步骤3次。将胰蛋白酶以样品：酶＝50:1(质量比)添加到试管中，并在37℃下进行酶解16小时。通过在12000g、4℃下离心获得肽样品，然后冷冻干燥。
52.酶解后的多肽样品在agilent 1100色谱仪中进行high ph反相色谱分级。柱子在初始条件下平衡5min，柱流速维持0.4ml/min；随后柱流速增至0.7ml/min，流动相位变化情况为5-5.10min/3％b，5.10-10min/5％b，10-35min/18％b，35-45min/34％b，45-58min/95％b，58-60min/3％b。使用1.5ml ep管，每1分钟接收1管。
53.另取酶解后的多肽样品经由orbitrap fusion lumos轨道阱质谱仪进行分析；共上样18μl，以60min梯度分离：3％溶液b至100％溶液b；柱流量控制在600nl/min；轨道阱质谱仪在数据依赖采集模式下运行，质谱参数设置如下：(1)msn1:扫描范围(m/z):375-1500；分辨率：120,000；agctarget:4e5；最大注入时间：50ms；(2)msn2：分辨率：15,000；agctarget＝5e4；最大注入时间：22ms；碰撞能量：30％；分级后的建库dda数据用proteome discoverer 2.1.0182进行蛋白质分析，最终获得dda参考谱图库；相关数据处理参数如下：搜库引擎采用sequestht；蛋白质数据库用uniprottaxid：9606)；设置胰蛋白酶酶解；最大允许漏切位点设为2；肽段水平的假阳性率(fdr)设置为1％；肽和片段质量误差分别为
±
10ppm和
±
0.02da。
54.另取酶解后的多肽样品经由orbitrap fusion lumos轨道阱质谱仪进行分析；共上样18μl，以60min梯度分离：3％溶液b至100％溶液b；柱流量控制在600nl/min；轨道阱质谱仪在数据依赖采集模式下运行，质谱参数设置如下：(1)msn1:扫描范围(m/z):350-1500；分辨率：60,000；agctarget:4e5；最大注入时间：50ms；(2)msn2：分辨率：30,000；agctarget＝3e5；最大注入时间：54ms；碰撞能量：33％；用skyline软件进行dia分析，通过导入上面所得参考谱图库，对dia原始数据进行多肽信号提取与定量分析。详细参数如下：多肽长度选6-25；子离子m/z要求为大于母离子且last ion-3；子离子最大数目设为5；子离子最小数目设为2；子离子抽提窗口采用5min；dotp需要大于等于0.6。
55.上述测定，对每个民族样本做三个平行样的生物重复，结果取平均值。
56.对两个民族的测定结果进行比较，将符合fc》1.5或《0.67，且pvalue《0.05的蛋白设定为差异表达蛋白，并进行功能注释和蛋白互作分析。
57.差异蛋白列表：
58.图1以柱状图的形式显示了总鉴定蛋白，差异蛋白的数目。相较于汉族组，朝鲜族母乳样品中共鉴定到54个差异表达蛋白，其中36个蛋白表达上调，18个蛋白表达下调。上调蛋白中，细胞间粘附分子-1具有最大的fc值5.50，p值为0.0076。它被认为在多种细胞上表达，同时也是许多关键生理过程所需的最重要的保守人类受体之一。wap四个二硫核心结构域蛋白2，则是有着第二高的fc值4.72，以及0.04的p值，这表明在两组之间存在显著差异。作为wap家族成员之一，它在肿瘤的进展、恶性和转移中起着重要作用，同时研究表明，它能促进卵巢癌的转移，并可能成为卵巢癌的潜在治疗靶点。下调蛋白中，磷酸丙糖异构酶是差异最大的蛋白质，有着最小的fc值0.03，p值为0.0027。它是糖酵解和糖异生过程中催化磷酸二羟丙酮(dhap)和甘油醛-3-磷酸(g3p)相互转化的关键糖酵解酶。除了催化作用外，它还对所有活细胞的碳水化合物代谢和能量生产起着至关重要的作用。
59.差异蛋白的go富集分析：
60.通过go富集分析分析朝鲜族和汉族母乳之间的差异表达蛋白，结果表明722个go项中共有357个显著富集(p《0.05)。最显著富集的前30个go项，包括19个bp，9个cc和2个mf。作为参与蛋白数量最多的项(16)，细胞的富集因子为0.007。其p值为最小值2.85e-08，表明它是最显著富集的go项。通过go富集分析研究了韩国和汉族母乳之间的dep，结果表明，722个go项中共有357个显著富集(p《0.05)。图2以气泡图的形式显示了前30个显著富集的go项，包括19bp、9cc和2mf。作为参与蛋白质数量最多的项(16)，细胞内的富集因子为0.007，最小p值为2.85e-08。
61.在54个dep中，rac1参与go富集分析的项数最多，共101项，其次是40s核糖体蛋白s27a(96)和β淀粉样a4蛋白(95)。
62.差异蛋白的kegg富集分析：
63.采用kegg富集分析研究了朝鲜族和汉族乳清蛋白之间的差异表达蛋白。图3结果表明，129个kegg途径中共有58个显著富集(p《0.05)。涉及6种蛋白质的病毒性心肌炎是最显著富集的途径，拥有最小p值4.14e-10，其富集因子为0.1。
64.差异蛋白的蛋白互作网络分析：
65.通过mcc法挖掘排名前10位的核心蛋白，按分数排序如下：40s核糖体蛋白s27a
(150)，60s核糖体蛋白l10a(150)，60s酸性核糖体蛋白p2(146)，40s核糖体蛋白s4(145)，40s核糖体蛋白s18(120)，60s核糖体蛋白l35(120)，蛋白酶体亚单位α1型(26)，肽基脯氨酸顺反异构酶(12)、l-乳酸脱氢酶(12)和磷酸丙糖异构酶(11)。
66.实施例3：酪蛋白的差异性分析
67.基于过滤辅助样品制备(fasp)方法，对实施例1得到的酪蛋白样品(100μg)进行蛋白质酶水解，向样品中添加8m尿素，使其体积达到200μl。引入二硫苏糖醇(dtt)至最终浓度为10mm，并使混合物在56℃下静置30分钟。将碘乙酸(iaa)添加到50mm的最终水平，并在室温下无光反应40min。将样品转移到10k超滤管中，然后在室温下以12000g离心，弃掉滤液，并重复该步骤三次。添加碳酸氢铵溶液(200μl，50mm)，然后在室温下以12000g离心，丢弃滤液，并重复该步骤3次。将胰蛋白酶以样品：酶＝50:1(质量比)添加到试管中，并在37℃下进行酶解16小时。通过在12000g、4℃下离心获得肽，然后冷冻干燥。
68.酶解后的多肽样品在agilent 1100色谱仪中进行high ph反相色谱分级。柱子在初始条件下平衡5min，柱流速维持0.4ml/min；随后柱流速增至0.7ml/min，流动相位变化情况为5-5.10min/3％b，5.10-10min/5％b，10-35min/18％b，35-45min/34％b，45-58min/95％b，58-60min/3％b。使用1.5ml ep管，每1分钟接收1管。
69.另取酶解后的多肽样品经由orbitrap fusion lumos轨道阱质谱仪进行分析；共上样18μl，以60min梯度分离：0min/3％b，0-2min/8％b，2-46min/28％b，46-55min/50％b，55-56min/100％b，56-60min/100％b。柱流量控制在600nl/min；轨道阱质谱仪在数据依赖采集模式下运行，质谱参数设置如下：(1)msn1:扫描范围(m/z):375-1500；分辨率：120,000；agctarget:4e5；最大注入时间：50ms；(2)msn2：分辨率：15,000；agctarget＝5e4；最大注入时间：22ms；碰撞能量：30％；分级后的建库dda数据用proteome discoverer 2.1.0182进行蛋白质分析，最终获得dda参考谱图库；相关数据处理参数如下：搜库引擎采用sequestht；蛋白质数据库用uniprottaxid：9606)；设置胰蛋白酶酶解；最大允许漏切位点设为2；肽段水平的假阳性率(fdr)设置为1％；肽和片段质量误差分别为
±
10ppm和
±
0.02da。
70.另取酶解后的多肽样品经由orbitrap fusion lumos轨道阱质谱仪进行分析；共上样18μl，以120min梯度分离：0min/3％b，0-2min/7％b，2-92min/22％b，92-107min/40％b，107-110min/100％b，110-120min/100％b。柱流量控制在600nl/min；轨道阱质谱仪在数据依赖采集模式下运行，质谱参数设置如下：(1)msn1:扫描范围(m/z):350-1500；分辨率：60,000；agctarget:4e5；最大注入时间：50ms；(2)msn2：分辨率：30,000；agctarget＝3e5；最大注入时间：54ms；碰撞能量：33％；用skyline软件进行dia分析，通过导入上面所得参考谱图库，对dia原始数据进行多肽信号提取与定量分析。详细参数如下：多肽长度选6-25；子离子m/z要求为大于母离子且last ion-3；子离子最大数目设为5；子离子最小数目设为2；子离子抽提窗口采用5min；dotp需要大于等于0.6。
71.上述测定，对每个民族样本做三个平行样的生物重复，结果取平均值。
72.对两个民族的测定结果进行比较，将符合fc》1.5或《0.67，且pvalue《0.05的蛋白设定为差异表达蛋白，并进行功能注释。
73.差异蛋白列表：
74.图4以柱状图的形式显示了总鉴定蛋白，差异蛋白的数目。通过比较朝鲜族和汉族
的酪蛋白中所鉴定蛋白质，观察到39种显著差异表达的蛋白质，其中有10种蛋白质上调，29种下调。在上调的蛋白质中，补体成分c9的fc最高，为5.74，p值为0.0013，这是两组之间差异最大的蛋白质。c9是母乳中补体的九种成分之一，在初乳中作为补体系统的防御因子发挥作用。而在下调的蛋白中，血清转铁蛋白的fc值最低，为0.19，p值为0.0045。作为一种铁结合糖蛋白，它除了能调节体液中的铁含量外，还与多种疾病有关，如无铁蛋白血症和心血管疾病。同时血清转铁蛋白及其受体的特性可以被用来向脑和癌细胞输送药物。
75.差异蛋白的go富集分析：
76.通过go富集分析分析朝鲜族和汉族母乳之间的差异表达蛋白，结果表示545个go项中共有125个显着富集(q《0.05)。图5以气泡图的形式显示了前30个显著富集的go项，其中bp的对刺激反应(7)，cc的细胞与外细胞器(8)，mf的作用于l-氨基酸肽的肽酶活性(3)和肽酶活性(3)分别是涉及蛋白数目最多的项。血液微粒是最富集的go项，其富集因子为0.235，具有最小q值2.58e-06。经验细胞凋亡过程的负调控(2)则具有最高的富集因子0.33。该项涉及的蛋白质为纤维蛋白原α链和纤维蛋白原γ链。纤维蛋白原是一种丰富的血浆糖蛋白，由α、β和γ链组成。纤维蛋白原是止血的关键分子之一。凝血酶介导的纤维蛋白肽从纤维蛋白原中释放，将这种可溶性蛋白转化为纤维蛋白纤维网络，形成血栓的构建块。纤维蛋白原除了参与凝血，还参与炎症、细胞迁移和肿瘤发生。
77.差异蛋白的kegg富集分析：
78.对于kegg富集分析，图6结果表明113个kegg途径中共有35个被显著富集(q《0.05)。金黄色葡萄球菌的富集因子为0.103，并且拥有最小q值为1.03e-10，是最富集的途径。而涉及蛋白质数量最多的途径是金黄色葡萄球菌感染和补体和凝血级联，这二者均涉及7种蛋白质。非洲锥虫病则是有最大富集因子0.135。
79.差异蛋白的蛋白互作网络分析：
80.差异表达蛋白的蛋白互作网络分析显示存在23个节点和47条边，表明有23个差异表达蛋白参与了47种相互作用。在这23种蛋白中，19个蛋白表达下调，4个蛋白表达上调。纤维蛋白原α链与9种差异蛋白存在相互作用，是相互作用最多的蛋白质。
81.通过mcc法挖掘排名前10位的核心蛋白，按分数排序如下：纤维蛋白原α链(1442)，补体成分c9(1442)，纤溶酶原(1441)，血清转铁蛋白(1441)，纤维蛋白原γ链(1440)，cdna flj35730 fis(1440)，补体c3(724)，α-2-巨球蛋白(721)，cd59糖蛋白(5)，聚合免疫球蛋白受体(4)。
82.实施例4：脂肪蛋白的差异性分析
83.基于过滤辅助样品制备(fasp)方法，对实施例1得到的脂肪蛋白样品(100μg)进行蛋白质酶水解，向样品中添加8m尿素，使其体积达到200μl。引入二硫苏糖醇(dtt)至最终浓度为10mm，并使混合物在56℃下静置30分钟。将碘乙酸(iaa)添加到50mm的最终水平，并在室温下无光反应40min。将样品转移到10k超滤管中，然后在室温下以12000g离心，弃掉滤液，并重复该步骤三次。添加碳酸氢铵溶液(200μl，50mm)，然后在室温下以12000g离心，丢弃滤液，并重复该步骤3次。将胰蛋白酶以样品：酶＝50:1(质量比)添加到试管中，并在37℃下进行酶解16小时。通过在12000g、4℃下离心获得肽，然后冷冻干燥。
84.酶解后的多肽样品在agilent 1100色谱仪中进行high ph反相色谱分级。柱子在初始条件下平衡5min，柱流速维持0.4ml/min；随后柱流速增至0.7ml/min，流动相位变化情
况为5-5.10min/3％b，5.10-10min/5％b，10-35min/18％b，35-45min/34％b，45-58min/95％b，58-60min/3％b。使用1.5ml ep管，每1分钟接收1管。
85.另取酶解后的多肽样品经由orbitrap fusion lumos轨道阱质谱仪进行分析；共上样18μl，以60min梯度分离，流速为0.60ml/min。梯度洗脱方法为：0min/3％b，0-2min/8％b，2-46min/28％b，46-55min/50％b，55-56min/100％b，56-60min/100％b。样品采用1.5ml ep管采集，间隔1min。轨道阱质谱仪在数据依赖采集模式下运行，质谱参数设置如下：(1)msn1:扫描范围(m/z):375-1500；分辨率：120,000；agctarget:4e5；最大注入时间：50ms；(2)msn2：分辨率：15,000；agctarget＝5e4；最大注入时间：22ms；碰撞能量：30％；分级后的建库dda数据用proteome discoverer 2.1.0182进行蛋白质分析，最终获得dda参考谱图库；相关数据处理参数如下：搜库引擎采用sequestht；蛋白质数据库用uniprottaxid：9606)；设置胰蛋白酶酶解；最大允许漏切位点设为2；肽段水平的假阳性率(fdr)设置为1％；肽和片段质量误差分别为
±
10ppm和
±
0.02da。
86.另取酶解后的多肽样品经由orbitrap fusion lumos轨道阱质谱仪进行分析；共上样18μl，以60min梯度分离：0min/3％b，0-2min/8％b，2-46min/28％b，46-55min/40％b，55-56min/100％b，56-60min/100％b。；柱流量控制在600nl/min；轨道阱质谱仪在数据依赖采集模式下运行，质谱参数设置如下：(1)msn1:扫描范围(m/z):350-1500；分辨率：60,000；agctarget:4e5；最大注入时间：50ms；(2)msn2：分辨率：30,000；agctarget＝3e5；最大注入时间：54ms；碰撞能量：33％；用skyline软件进行dia分析，通过导入上面所得参考谱图库，对dia原始数据进行多肽信号提取与定量分析。详细参数如下：多肽长度选6-25；子离子m/z要求为大于母离子且last ion-3；子离子最大数目设为5；子离子最小数目设为2；子离子抽提窗口采用5min；dotp需要大于等于0.6。
87.上述测定，对每个民族样本做三个平行样的生物重复，结果取平均值。
88.对两个民族的测定结果进行比较，将符合fc》1.5或《0.67，且pvalue《0.05的蛋白设定为差异表达蛋白，并进行功能注释。
89.差异蛋白列表:
90.图7以柱状图的形式显示了总鉴定蛋白，差异蛋白的数目。通过比较朝鲜族和汉族的酪蛋白中所鉴定蛋白质，观察到238种显著差异表达的蛋白质，其中有142种蛋白质上调，96种下调。在上调的蛋白质中，聚合酶蛋白同源样结构域，家族a，成员1，亚型cra_a(phlda1)的fc最高，log2fc值为26.18，p值为0.000776，这是两组之间差异最大的蛋白质。phlda1是一种具有多种功能的蛋白，在调节细胞凋亡、增殖、分化和细胞迁移的许多生物学过程中发挥重要作用，因此在不同类型的癌症中表达改变。而在下调的蛋白中，酪氨酸蛋白激酶受体的fc值最低，log2fc值为-7.35，p值为0.01998。酪氨酸蛋白激酶受体可以协调多种细胞功能，如细胞增殖和分化，因此在多种细胞的起始和进化中发挥重要作用
91.差异蛋白的go富集分析：
92.通过go富集分析分析朝鲜族和汉族母乳之间的差异表达蛋白，结果表示共有374个显着富集(q《0.05)。图8显示了前三十项并对其进行研究，其中bp的细胞代谢过程(30)，cc的膜有界细胞器(37)，mf的催化活性(25)分别是涉及蛋白数目最多的项。细胞外间隙是最富集的go项，其富集因子为0.045，具有最小q值2.27e-15。过氧化物还蛋白活性则具有最高的富集因子0.8。该项涉及的4个差异蛋白均为过氧化物酶，属于硫醇特异性过氧化物酶，
分别催化过氧化氢和有机氢过氧化物还原为水和醇。通过解毒过氧化物和作为过氧化氢介导的信号事件的传感器，在细胞保护氧化应激中发挥作用。
93.差异蛋白的kegg富集分析：
94.对于kegg富集分析，共有118个被显著富集(q《0.05)。图9显示了前二十项并对其进行研究，数据显示，吞噬小体拥有最小q值为2.20e-12，是富集最显著的途径。而涉及蛋白质数量最多的途径是代谢途径，涉及33个差异蛋白。哮喘则是有最大富集因子0.194。过敏性哮喘是最常见的哮喘表型，通常始于儿童时期，通常同时伴有特应性皮炎和变应性鼻炎等并发症状
95.差异蛋白的蛋白互作网络分析：
96.差异表达蛋白的蛋白互作网络分析显示存在202个节点和979条边，表明有202个差异表达蛋白参与了979种相互作用。对ppi网络与对应的前20个枢纽基因的蛋白质的互作网络进行研究，数据表明存在20个节点和143条边。
97.通过degree法挖掘排名前20位的核心蛋白，按分数排序如下：肌动蛋白，细胞质1，热休克蛋白hsp 90-alpha，热休克同源物71kda蛋白，白蛋白，延伸因子1-alpha 1，40s核糖体蛋白sa，60s核糖体蛋白l9，t复合体蛋白1亚基，40s核糖体蛋白s4，40s核糖体蛋白s6，钙联接蛋白，60s核糖体蛋白l18，60s核糖体蛋白l27a，40s核糖体蛋白s19，增殖相关蛋白2g4，60s核糖体蛋白l15，真核翻译起始因子3亚基m，真核生物起始因子4a-i，真核生物起始因子4a-i以及40s核糖体蛋白s25。