引言：不同生态位下的乳汁演化
乳汁是哺乳动物幼崽的主要营养来源，其成分演化以满足物种特定的生态和发育需求。在干旱和半沙漠生态系统中，骆驼奶是一种独特的营养基质，支持幼驼在极端温度波动、水源稀缺和饲料有限的环境下生长。与牛、羊等在持续干旱中生存艰难的反刍动物相比，骆驼展现出非凡的韧性，在极端条件下仍能维持长期的产奶能力。作为伪反刍动物，骆驼的胃肠道结构使其能更高效地利用粗纤维，从低质量饲料中提取更多的脂肪酸前体。因此，骆驼奶中长链脂肪酸比例高达约96.4%，而牛奶中为85.3%。然而，仅凭脂肪酸组成的差异难以全面洞察乳脂的结构组织和功能潜力。脂质成分是乳汁功能属性的核心，它们不仅是高密度的能量储备和细胞膜的结构成分，也是细胞内信号分子的前体。特定脂质种类在调节肠道稳态和神经发育中起着关键作用，这些过程对于幼崽在恶劣环境中的生存至关重要。
结果与讨论
2.1. 分析方法验证
通过质量控样品的稳定性和重现性分析，确认了超高效液相色谱-串联质谱平台的可靠性，为后续脂质组学数据的稳健性奠定了基础。
2.2. 比较脂质组学分析
非靶向脂质组学分析在骆驼奶和牛奶中共鉴定出2460种脂质化合物。根据LIPID MAPS分类系统，主要分为五大类：甘油磷脂（GPs，占50.00%）、鞘脂（SPs，20.46%）、脂肪酸（FAs，11.36%）、固醇脂（STs，9.09%）和甘油脂（GLs，9.09%）。其中，GP是乳汁中的主要脂质类别。亚类分析进一步识别出44个脂质亚类，其中甘油三酯（TG，占35.49%，873种）、磷脂酰胆碱（PC，18.25%，449种）和二酰甘油（DG，11.42%，281种）占主导地位。尽管某些亚类如Hex1Cer（0.57%）含量极低，但它们可能具有重要的生物学调节功能。
2.3. 乳汁脂质组的多元统计分析
主成分分析和偏最小二乘判别分析显示，骆驼奶与牛奶的脂质组存在显著差异。在正、负离子模式下，两组样品均能清晰分离，且PLS-DA模型显示出高预测精度和可靠性，为寻找物种特异性生物标志物建立了稳健框架。
2.4. 差异脂质物种的鉴定
通过严格的统计标准，共鉴定出498个差异脂质物种，涵盖18个亚类。骆驼奶中显著富集了18种脂质，全部属于甘油磷脂类，包括6种PC、6种PE、2种PI等。相反，牛奶中显著富集了480种脂质，其中411种是TG。对脂质亚类丰度的定量分析进一步揭示了两者的核心差异：在骆驼奶中，PC（约49%）和PI（约22%）占主导；而在牛奶中，TG占比高达97.84%，PC、PE和SM等磷脂各占不足1%。
2.5. 候选生物标志物的发现与验证
为寻找可区分两种乳品的生物标志物，研究根据差异倍数对脂质进行了排序。骆驼奶中PI(18:0/22:3)和PE(18:0/22:3)极度富集，而牛奶中TG(4:0/14:0/18:0 CHO)和TG(4:0/10:0/12:0)则呈现数量级上的优势上调。利用这四种脂质的定量浓度数据构建正交偏最小二乘判别分析模型，能清晰区分两类样品，置换检验证实模型有效且未过拟合。受试者工作特征曲线分析显示，这四种标志物组合的曲线下面积达到1.0，表现出卓越的分类性能。
2.6. 比较氧化脂质组学分析
靶向定量分析揭示了源自多不饱和脂肪酸的三类氧化脂质的差异。骆驼奶中含有显著更高水平的PUFA衍生物，特别是11,12-diHETrE和5-HEPE。这些氧化脂质是炎症、氧化应激反应和免疫调节的关键介质，其存在可能与骆驼适应极端环境压力有关。相比之下，牛奶中8-HETrE等氧化脂质水平较高。这种组成差异可能源于牛在瘤胃消化过程中发生广泛的生物氢化，降低了乳中不饱和脂肪酸的丰度；而作为伪反刍动物的骆驼则保留了更高比例的非修饰不饱和脂肪酸及其生物活性氧化衍生物。
2.7. 脂质相关性网络
相关性网络分析揭示了脂质在物种内部，尤其是亚类内部的强烈相互作用。骆驼奶中显示出比牛奶更多的正相关性，特别是在PC、PE和PI类别内部。这些紧密协调的磷脂模块可能对幼驼在干旱环境下的神经发育和肠道健康适应具有潜在贡献。
2.8. KEGG通路富集分析
通路富集分析比较了骆驼奶与牛奶差异脂质相关的代谢通路，鉴定出35条KEGG代谢通路，其中11条是显著富集的脂质代谢通路。骆驼奶脂质显著富集于甘油磷脂代谢和磷脂酰肌醇信号系统等通路，这与其高丰度的PC和PE结构相符，这些分子对膜生物发生和突触可塑性至关重要。相比之下，牛奶脂质主要富集于脂肪消化与吸收、甘油脂代谢和胰岛素抵抗等与能量相关的通路，这与其TG主导的组成一致，旨在为温带环境下出生后的快速发育提供快速的能量动员支持。
材料与方法
3.1. 样品采集
研究在中国新疆裕民县江阿布拉克村的半沙漠草原进行。样品采集自8头健康的准噶尔双峰驼和8头新疆褐牛，所有动物均为处于泌乳期（产后2-3个月）的健康雌性，在相同的自由放牧条件下饲养，主要采食耐旱植被。所有乳汁样品在清晨挤奶时段采集，经液氮速冻后于-80°C保存。
3.2. 仪器与试剂
色谱系统采用Vanquish™ UHPLC和ExionLC™ AD系统，质谱仪为Q Exactive™ HF-X和QTRAP®6500+。使用了Accucore C30和Atlantis Premier BEH C18色谱柱。试剂均为高效液相色谱级。
3.3. 样品前处理
脂质提取采用改良的Matyash方法。取100 µL乳汁，与甲醇、甲基叔丁基醚及SPLASH™ LIPIDOMIX™内标混合，涡旋、离心后收集有机相，氮吹干燥，用异丙醇复溶供LC-MS/MS分析。氧化脂质提取则采用固相萃取柱进行。
3.4. UHPLC-MS/MS分析
为确保数据质量，每批次分析均包含程序空白以监控背景污染，所有样品采用随机进样顺序，并穿插质控样品以监控仪器稳定性。非靶向脂质组学在正、负离子模式下进行，扫描范围m/z 114–1700。靶向氧化脂质组学在正离子模式下进行。
3.5. 统计分析
使用LipidSearch 5.1软件进行峰对齐和脂质定量。采用metaX平台进行主成分分析和偏最小二乘判别分析等多变量统计分析。差异脂质的筛选标准为变量重要性投影值大于1.0，p值小于0.05，且差异倍数绝对值大于1。使用R语言进行热图绘制、相关性网络构建和通路富集分析。脂质注释参考了京都基因与基因组百科全书、人类代谢组数据库和LIPID MAPS数据库。
结论
本研究综合运用非靶向脂质组学和靶向氧化脂质组学方法，系统描绘了骆驼奶与牛奶的脂质谱。研究发现两者存在深刻的组成差异：骆驼奶显著富集18种甘油磷脂，而牛奶则被480种脂质（主要是甘油三酯）主导。定量亚类分析表明，磷脂酰胆碱和磷脂酰肌醇构成了骆驼奶的主要脂质类别，与牛奶中甘油三酯占绝对主导的图谱形成鲜明对比。研究鉴定出的四种脂质——PI(18:0/22:3)、PE(18:0/22:3)、TG(4:0/14:0/18:0 CHO)和TG(4:0/10:0/12:0)——可作为区分骆驼奶和牛奶的有效候选生物标志物。通路富集分析揭示了11条显著改变的脂质代谢通路，主要与甘油磷脂代谢相关。这些结果从分子层面阐释了骆驼奶脂质组成可能具有的生态适应意义，加深了我们对极端环境下骆驼生理的理解，并为探索提升骆驼奶生产及其潜在营养与功能价值的策略提供了有价值的参考数据。
