1 材料与方法
1.1 试验动物及饲养管理
1.2 血液采集与血浆制备
1.3 血浆生化指标的测定
1.4 血浆激素浓度的测定
1.5 血浆代谢组的测定
1.6 数据处理与分析
2 结果与分析
2.1 不同乳蛋白产量奶牛的血浆生化指标
表1 不同乳蛋白产量奶牛的血浆生化指标Table 1 Plasma biochemical indices in dairy cows with different milk protein yields |
| 项目 Items | 组别Groups | 均值标准误 SEM | P值 P-value | |
|---|---|---|---|---|
| LMP | HMP | |||
| 谷丙转氨酶ALT/(U/L) | 29.21 | 21.28 | 1.452 | 0.002 |
| 谷草转氨酶AST/(U/L) | 68.35 | 56.43 | 5.615 | 0.078 |
| 总蛋白TP/(g/L) | 74.29 | 77.37 | 3.307 | 0.388 |
| 白蛋白ALB/(g/L) | 27.53 | 27.21 | 1.055 | 0.777 |
| 球蛋白GLB/(g/L) | 46.76 | 50.15 | 3.194 | 0.330 |
| 葡萄糖GLU/(mmol/L) | 3.84 | 3.71 | 0.104 | 0.242 |
| 乳酸脱氢酶LDH/(U/L) | 551.82 | 497.18 | 57.796 | 0.381 |
| 甘油三酯TG/(mmol/L) | 0.15 | 0.19 | 0.021 | 0.079 |
| 肌酐CR/(μmol/L) | 125.24 | 118.04 | 13.411 | 0.610 |
| 总胆固醇TC/(mmol/L) | 5.86 | 3.40 | 0.412 | 0.001 |
| 低密度脂蛋白LDL/(mmol/L) | 5.56 | 3.05 | 0.500 | 0.002 |
| 高密度脂蛋白HDL/(mmol/L) | 1.20 | 0.81 | 0.089 | 0.005 |
| 尿素氮UN/(mmol/L) | 5.44 | 4.02 | 0.452 | 0.020 |
| 总胆汁酸TBA/(μmol/L) | 10.47 | 13.90 | 2.075 | 0.150 |
| 总胆红素TBIL/(μmol/L) | 0.76 | 0.83 | 0.067 | 0.351 |
| 谷丙谷草比ALT/AST | 2.34 | 2.66 | 0.218 | 0.198 |
| 白球比ALB/GLB | 0.59 | 0.55 | 0.043 | 0.399 |
2.2 不同乳蛋白产量奶牛的血浆激素浓度
表2 不同乳蛋白产量奶牛的血浆激素浓度Table 2 Plasma hormone concentrations in dairy cows with different milk protein yields |
| 项目 Items | 组别Groups | 均值标准误 SEM | P值 P-value | |
|---|---|---|---|---|
| LMP | HMP | |||
| 生长激素GH/(ng/mL) | 14.03 | 8.76 | 1.430 | 0.010 |
| 催乳素PRL/(mIU/L) | 445.30 | 673.54 | 91.431 | 0.047 |
| 胰岛素INS/(mIU/L) | 26.94 | 26.72 | 4.774 | 0.964 |
| 胰岛素生长因子-1 IGF-1/(ng/mL) | 199.64 | 208.43 | 38.323 | 0.826 |
2.3 血浆代谢组鉴定
2.3.1 数据预处理
2.3.2 PCA和PLS-DA置换检验
图2 样本质控主成分分析散点图图A为正离子模式下分析结果,图B为负离子模式下分析结果。LMP:LMP组样本;HMP:HMP组样本;QC:所有待测溶液的等体积混合样本。下图同。 Fig.2 Sample quality control PCA scatter plot Figure A shows the analysis results under positive ion mode, while figure B shows those under negative ion mode. LMP: samples of LMP group; HMP: samples of HMP group; QC: an equal-volume mixture of all tested solutions. The same as below. |
图4 样本PLS-DA置换检验图横坐标为置换保留比例,纵坐标为R2Y 和 Q2。保留比例为 1.00时对应的点即为原始模型的 R2Y 与 Q2值。 Fig.4 Sample PLS-DA permutation test plot The X-axis represents the permutation retention ratio, and the Y-axis represents R2Y and Q2. The point corresponding to a retention ratio of 1.00 indicates the R2Y and Q2 values of the original model. |
2.3.3 不同乳蛋白产量奶牛的血浆差异代谢物
图5 差异代谢物火山图红色的点代表发生下调的差异代谢物(VIP>1且P<0.05),绿色的点代表发生上调的差异代谢物(VIP>1且P<0.05)。 Fig.5 Volcano plots of differential metabolites Red dots represent significantly downregulated differential metabolites with VIP>1 and P<0.05, while green dots denote significantly upregulated differential metabolites with VIP>1 and P<0.05. |
表3 正离子模式下血浆差异代谢物Table 3 Plasma differential metabolites under positive ion mode |
| 代谢物 Metabolites | 组别Groups | 变量重要 性投影 VIP | P值 P-value | |
|---|---|---|---|---|
| LMP | HMP | |||
| 磷脂酰胆碱18∶2 LPC18∶2 | 33 602 752.09 | 21 808 933.68 | 35.872 | 0.009 |
| 甘油磷酰胆碱Glycerophosphocholine | 5 980 027.21 | 3 124 603.14 | 17.664 | 0.006 |
| 齿阿米素Visnagin | 21 696 534.45 | 23 897 376.96 | 14.982 | 0.025 |
| 拉莫三嗪Lamotrigine | 7 457 284.35 | 8 846 136.37 | 11.875 | 0.026 |
| 2-十六烷酰硫基-1-乙基磷酰胆碱 2-hexadecanoylthio-1-ethylphosphorylcholine | 1 368 725.98 | 202 623.77 | 11.459 | 0.002 |
| 异乐灵Isopropalin | 1 849 962.26 | 773 950.31 | 11.210 | 0.001 |
| 1-硬脂酰-2-油酰-sn-甘油-3-磷酸乙醇胺 1-stearoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine | 1 873 778.87 | 805 648.50 | 10.985 | 0.002 |
| 1-硬脂酰-2-油酰-sn-甘油-3-磷酸胆碱 1-octadecanoyl-2-octadecenoyl-sn-glycero-3-phosphocholine | 2 062 582.85 | 3 143 401.72 | 10.598 | 0.014 |
| 伊立玛醌Ilimaquinone | 5 876 009.66 | 6 855 208.90 | 9.932 | 0.026 |
| 1-(1Z-十六碳烯基)-sn-甘油-3-磷酸胆碱 1-(1Z-hexadecenyl)-sn-glycero-3-phosphocholine | 1 310 960.11 | 506 206.92 | 9.703 | <0.001 |
| 2-磷酸甘油酸2-phosphoglycerate | 12 283 841.01 | 13 146 703.48 | 9.214 | 0.038 |
| 硫醚磷脂酰胆碱Thioetheramide-PC | 1 800 322.28 | 1 123 735.52 | 8.578 | 0.010 |
| 花生四烯酰硫代磷酰胆碱Arachidonoylthiophosphorylcholine | 3 512 640.88 | 2 770 775.46 | 8.499 | 0.046 |
| (3R,6R)-6-[(3R,5S,7R,8R,9S,10S,12S,13R,14S,17R)-3,7-二乙酰氧基-12-羟基-10,13-二甲基十六氢-1H-环戊烷并[a]菲-17-基]庚烷-1,3-二基二乙酸酯 (3R,6R)-6-[(3R,5S,7R,8R,9S,10S,12S,13R,14S,17R)-3,7-diacetoxy-12-hydroxy-10,13-dimethylhexadecahydro-1H-cyclopenta[a]phenanthren-17-yl]heptane-1,3-diyl diacetate | 644 795.34 | 48 678.48 | 8.378 | <0.001 |
| 诺米林Nomilin | 1 004 745.92 | 393 522.71 | 8.320 | 0.002 |
| 1-油酰-sn-甘油-3-磷酸胆碱 1-oleoyl-sn-glycero-3-phosphocholine | 1 911 789.63 | 2 596 519.74 | 8.024 | 0.049 |
| 焦谷氨酰-苯丙氨酰-精氨酸Pyroglu-Phe-Arg | 1 004 538.61 | 560 503.08 | 7.200 | <0.001 |
| 月桂酸Dodecanoic acid | 1 754 126.39 | 2 205 260.67 | 7.139 | 0.003 |
| 1-花生四酰基-2-羟基-sn-甘油-3-磷酸胆碱 1-arachidoyl-2-hydroxy-sn-glycero-3-phosphocholine | 633 073.11 | 235 372.62 | 6.285 | 0.027 |
| 氯溴隆Chlorbromuron | 3 738 299.60 | 3 364 727.53 | 6.123 | 0.040 |
| 醋酸甲地孕酮Megestrol acetate | 330 474.21 | 56 529.03 | 5.628 | <0.001 |
| 1-十五烷酰基-sn-甘油-3-磷酸胆碱 1-pentadecanoyl-sn-glycero-3-phosphocholine | 969 311.76 | 666 079.38 | 5.393 | 0.043 |
| 达托利昔布Dactolisib | 438 658.44 | 213 292.70 | 5.148 | <0.001 |
| 丙溴磷Profenofos | 1 473 122.48 | 1 233 625.45 | 5.102 | 0.008 |
表4 负离子模式下血浆差异代谢物Table 4 Plasma differential metabolites under negative ion mode |
| 代谢物 Metabolites | 组别Groups | 变量重 要性投影 VIP | P值 P-value | |
|---|---|---|---|---|
| LMP | HMP | |||
| 6-硫酸氨基葡萄糖D-glucosamine,6-sulfate | 542 361.80 | 199 584.81 | 10.713 | 0.014 |
| 香柑醇Bergaptol | 479 070.08 | 134 008.15 | 10.542 | 0.024 |
| 1-硬脂酰-sn-甘油-3-磷酰胆碱 1-stearoyl-2-hydroxy-sn-glycero-3-phosphocholine | 747 908.94 | 458 564.09 | 10.364 | 0.003 |
| 雌三醇E3 | 25 849.40 | 327 408.42 | 10.116 | 0.013 |
| 绿僵菌素A Destruxin A | 424 887.79 | 180 535.54 | 9.602 | 0.001 |
| 刃天青Resazurin | 170 697.08 | 402 599.48 | 9.362 | <0.001 |
| 1-[2,4,6-三羟基-3-(7-羟基-2-(4-羟基苯基)-3,4-二氢-2H-色烯-4-基)苯基]十二烷-1-酮 1-[2,4,6-trihydroxy-3-[7-hydroxy-2-(4-hydroxyphenyl)-3,4-dihydro-2H-chromen-4-yl]phenyl]dodecan-1-one | 23 178.82 | 253 728.23 | 8.990 | 0.007 |
| 异秦皮碱Isogentisin | 597 549.17 | 376 622.61 | 8.830 | 0.004 |
| (2-氨基乙氧基)[2-[(9Z)-十六碳-9-烯酰氧基]-3-[十八碳-1,11-二烯-1-基氧基]丙氧基]次膦酸 (2-aminoethoxy)[2-[hexadec-9-enoyloxy]-3- [octadeca-1,11-dien-1-yloxy]propoxy]phosphinic acid | 305 068.48 | 99 588.21 | 8.659 | 0.007 |
| 2-[(10E,12Z)-9-羟基十八碳-10,12-二烯酰]氧基-3-[(Z)-十八碳-1-烯氧基]丙基2-(三甲基铵)乙基磷酸酯 Pc(18∶1e/9-HODE) | 656 441.44 | 463 913.85 | 8.252 | 0.012 |
| L-氢化乳清酸L-dihydroorotate | 686 827.33 | 526 567.14 | 6.937 | 0.042 |
| 白三烯B4 Leukotriene B4 | 151 340.45 | 33 511.35 | 6.533 | 0.003 |
| 2-花生四烯酰-1-硬脂酰-sn-甘油 2-arachidonoyl-1-stearoyl-sn-glycerol | 127 504.84 | 239 985.53 | 6.338 | 0.004 |
| 光甘草内酯Glabrolide | 36 373.98 | 152 800.24 | 6.145 | 0.026 |
| 1-(9Z,12Z-十八碳二烯酰)-2-羟基-sn-甘油-3-磷酸乙醇胺 1-(9Z,12Z-octadecadienoyl)-2-hydroxy-sn-glycero-3- phosphoethanolamine | 153 032.34 | 62 695.52 | 5.749 | 0.006 |
| 3-[(Z)-十六碳-1-烯氧基]-2-[(10E,12Z)-9-羟基十八碳-10,12-二烯酰]氧基丙基2-(三甲基铵)乙基磷酸酯 Pc(16∶1e/9-HODE) | 422 008.18 | 332 753.33 | 5.531 | 0.013 |
| 苯乙酰甘氨酸Phenaceturic acid | 153 179.37 | 61 451.21 | 5.494 | 0.020 |
| 亚牛磺酸Hypotaurine | 218 269.27 | 124 739.59 | 5.331 | 0.038 |
| 1-硬脂酰基-2-羟基-sn-甘油-3-磷酸乙醇胺 1-stearoyl-2-hydroxy-sn-glycero-3-phosphoethanolamine | 211 761.64 | 136 617.20 | 5.302 | 0.001 |
| 3-[(4-氯苯甲酰基)氨基]丙酸 3-[(4-chlorobenzoyl)amino]propanoic acid | 162 225.80 | 81 735.53 | 5.150 | 0.030 |
2.3.4 不同乳蛋白产量奶牛的血浆差异代谢通路
图6 KEGG通路富集分析柱状图横坐标为通路经过超几何检验后得到的P值取-log10后的数值(-log100.05≈1.301 03,-log100.01=2),纵坐标为通路名称。 Fig.6 KEGG pathway enrichment analysis bar chart The X-axis represents the -log10(P-value) obtained from hypergeometric tests on pathways (-log100.05≈1.301 03,-log100.01=2), while the Y-axis displays pathway names. |
