小鼠血清代谢组遗传结构揭示羧酸酯酶（carboxyl esterases）作为循环脂肪酸代谢的主调控因子

《PLOS Genetics》：Genetic architecture of the murine serum metabolome reveals carboxyl esterases as master regulators of circulating fatty acid metabolism

　　

背景：循环代谢物和脂质的系统性生化多样性反映了遗传变异与环境暴露的综合效应。人类代谢物数量性状基因座（metabolite quantitative trait locus, mQTL）研究已建立基因-代谢物关联，但性别、饮食、年龄、用药和环境暴露可掩盖遗传贡献。遗传多样性模型系统为分离循环代谢组生化多样性的遗传贡献提供了强有力的补充策略。方法/主要发现：研究人员对来自多样性远交群（Diversity Outbred, DO）群体的541只小鼠血清样本进行质谱分析，在240种代谢物、561种脂质、43种氧化脂质（oxylipins）和4,465个串联质谱（MS/MS）特征中鉴定出1,933个mQTL。8、9、10和17号染色体上的QTL热点提示不同位点参与血清脂质、代谢物和肽样特征的调控。8号染色体Ces1和Ces2基因簇与脂质、酰基肉碱（acylcarnitines）和酰基辅酶A（acyl-CoAs）相关，而9号和10号染色体位点包括其他脂质和代谢物QTL。17号染色体上主要组织相容性复合体（major histocompatibility complex, MHC）和补体C3（C3）区域与未靶向肽样MS/MS特征相关。与匹配的红细胞（red blood cell）mQTL比较揭示广泛的区室特异性遗传控制。结论/意义：总体而言，这些发现提供了循环代谢组遗传调控的高分辨率图谱，提供了补充和扩展人类代谢遗传学的机制性见解。

循环代谢组和脂质组是连接遗传变异与复杂性状和疾病的动态中间表型。人类mQTL研究已识别许多调控代谢物和脂质水平的位点，并阐明与心血管疾病、糖尿病和神经退行性疾病相关的通路，包括支链氨基酸代谢、脂肪酸代谢和脂质运输。然而，人类研究受有限遗传多样性和异质性环境暴露（exposome）制约，生活方式、饮食、药物等强烈调节循环代谢组，使遗传效应与环境效应难以区分；人类GWAS候选基因的功能验证亦存在障碍。经典近交小鼠虽已提供代谢调控见解，但遗传多样性有限；多亲本群体如协作交叉（Collaborative Cross, CC）及其远交衍生DO小鼠保留实验控制并捕获八个创始株的遗传多样性，近似自然群体规模变异，支持复杂性状高分辨率映射。既往DO研究已映射影响脂质代谢、葡萄糖稳态、胰岛素敏感性和衰老的位点，但调控循环代谢物和脂质的更广泛遗传网络仍未完全理解，既往多映射单个QTL而未充分定义代谢变异的遗传结构。为此，研究人员利用DO群体高分辨率映射血清代谢物和脂质的遗传决定因素，结合靶向与非靶向代谢组学、脂质组学和氧化脂质分析识别循环生化多样性新遗传驱动因素，并利用其他DO队列的基因和蛋白表达数据定义和优先候选基因；此前该队列分子谱分析聚焦红细胞，本研究是同一批DO小鼠首个以血清为重点的代谢物和脂质谱，旨在推进代谢调控理解并为未来机制和转化研究奠基。

关键技术方法概括（不超过250字）：研究纳入来自Jackson Laboratory的DO小鼠（第45和46代；541只有血清样本，298雌、243雄；2只未通过基因型质控，539只用于遗传分析），小鼠自由采食、未标准化禁食时间。采用质谱代谢组学、脂质组学和氧化脂质分析血清；未靶向MS/MS特征用Proteome Discoverer 3.2中MASCOT和SEQUEST HT检索UniProt Mus musculus数据库。用R/qtl2进行QTL基因组扫描，LOD>7为宽松阈值，LOD>10且4 Mbp内至少8个QTL定义热点；用MOFA整合组学层；整合来自Svenson、Attie和Aging三个DO队列的13个肝脏、脂肪、胰岛、骨骼肌、心脏、肾脏和骨eQTL/pQTL数据集优先候选基因；并与同一队列已发表红细胞QTL比较。

研究结果：
1. 血清生化特征的性别差异：通过线性混合模型检测，FDR<0.01的显著性别差异见于119种代谢物、297种脂质、30种氧化脂质和1,238个未鉴定MS/MS特征；脂肪酸通路中多数显著特征在雄性更高，雌性相关例子包括trimethylamine N-oxide和kynurenine；因性别与批次部分混杂，研究人员强调应谨慎解释。
2. QTL分析揭示驱动血清生化变异的独特热点位点：窄义遗传力中位数代谢物19.9%、脂质18.6%，氧化脂质和未鉴定MS/MS特征较低；检测到163个代谢物QTL、324个脂质QTL、17个氧化脂质QTL和1,429个未鉴定MS/MS QTL，并定义7个热点：8号染色体94和105.5 Mbp、9号染色体96 Mbp、10号染色体75.5 Mbp、1号染色体171 Mbp、17号染色体36和58 Mbp。
3. MOFA揭示QTL热点跨血清生化特征的共调控模式：MOFA识别12个解释>1%总体变异的潜在因子；潜在因子1捕获脂质50.4%、未靶向MS/MS 10.1%、代谢物5.9%的变异，并映射至8号染色体94 Mbp热点、解释其14.4%变异；潜在因子6解释未靶向MS/MS 9.7%变异并映射至17号染色体36 Mbp热点、解释58.4%变异，提示部分热点跨多组学层而部分层特异。
4. 染色体8 QTL热点：羧酸酯酶基因簇调控循环脂质和肉碱水平：Ces1热点（94 Mbp）代表61个QTL（LOD>7为127个），映射31种脂质、1种氧化脂质TUDCA和29个未靶向特征；Ces2热点（105.5 Mbp）代表41个QTL（LOD>7为45个），映射21种酰基肉碱、5种酰基辅酶A和15个未靶向特征；两热点haplotype effects不同，Ces1b、Ces1c和Ces2h为优先候选，肝脏local eQTL/pQTL倒置提示其高表达与低血清水平相关，另一DO队列肝脏和血浆脂质QTL支持Ces2。
5. 染色体9 QTL热点：鉴定调控脂质代谢的新位点：96 Mbp热点代表9个QTL（LOD>7为12个），包括Gk5，仅映射磷脂酰胆碱和甘油三酯等脂质；haplotype effects近双等位，NZO、CAST和PWK携带与较低脂质水平相关等位基因；候选包括Gk5、Pls1、Atr、Xrn1、Tfdp2和Atp1b3，local eQTL/pQTL提供部分支持。
6. 染色体10 QTL热点：氨基酸代谢物QTL富集：75.5 Mbp热点代表30个QTL（LOD>7为39个），映射胆碱、生物素、甘油磷酸乙醇胺和胞苷等代谢物及未靶向特征；haplotype effects类似Ces2热点，候选基因涉及线粒体氨基酸分解代谢和氮稳态，包括Ggt、Ddt和Susd2，local eQTL/pQTL支持Ggt5、Susd2、Gm20441、Gstt3和Gstt1。
7. 未靶向MS/MS特征QTL热点揭示多态性肽和蛋白特征：未靶向特征映射至与已鉴定代谢物和脂质共享的热点（9号染色体除外），并独特映射至1号染色体171 Mbp和17号染色体36与58 Mbp；独特热点基于完整质量电离态[M + 2H]⁺²和[M + 3H]⁺³富集肽；Mascot/SEQUEST检索获得70个高置信肽谱匹配、30个肽条目和28个独特小鼠蛋白登录号；1号染色体热点中Apoa2的错义变异株分布模式与肽样信号haplotype effects高度一致。
8. 染色体17 MHC热点与肽样MS/MS特征相关：36 Mbp热点有14个未靶向MS/MS QTL（LOD>7为44个），跨越MHC class II位点，多等位haplotype effects提示高度遗传异质性；58 Mbp热点由8个QTL（LOD>7为11个）组成，其中7个强QTL LOD>20且呈肽样，区域包括C3，双等位haplotype effects区分野生来源与经典实验室株，候选错义变异为rs107650602和rs52547087；结果提示免疫基因变异可能通过补体激活或蛋白水解降解等机制影响血清肽丰度。
9. 血清-红细胞代谢物比较识别有限保守遗传调控：与同一队列已发表红细胞QTL比较，血清QTL与红细胞QTL大多不同，无热点重叠，无脂质或氧化脂质QTL共享；157个共同定量代谢物中仅N⁶-methyl-L-lysine和N-Acetyl-D-phenylalanine的4个QTL位置一致且血清与红细胞水平相关，提示少数代谢物存在稳定系统调控，而多数性状呈基质或区室特异。

讨论部分总结：染色体8 Ces1和Ces2热点均与脂质、肉碱和酰基辅酶A相关，但haplotype effects和分子谱不同；Ces1调控脂质和一种氧化脂质，Ces2调控肉碱和酰基辅酶A，TUDCA在Ces1内独特haplotype effect提示Ces1旁系同源物底物特异功能。Ces基因既往多研究异源物和前药代谢，近期研究提示Ces1d和Ces1g参与肝甘油三酯水解和脂蛋白分泌；eQTL/pQTL整合优先Ces1b、Ces1c和Ces2h，组织表达与haplotype效应一致。9号染色体Gk5为有趣候选，但主要表达于皮肤，需未来功能研究；10号染色体热点富集氨基酸代谢物，候选Ggt和Gstt涉及谷胱甘肽代谢和含氮化合物处理；17号染色体MHC热点机制未明，C3热点有两个错义变异与haplotype效应匹配。研究局限包括DO小鼠环境暴露和进化历史异于人类、未靶向工作流需正交验证、LD区域内无法完全验证因果基因、自由采食和未标准化禁食时间、性别与批次部分混杂、体重未记录且不能独立于体重解释QTL。研究结论翻译：尽管存在这些局限，研究结果提供了代谢调控多基因架构的广泛见解；通过在遗传多样性DO队列中提供参考图谱，本研究为解释受暴露组异质性混杂的人类mQTL信号提供框架；识别已建立代谢酶的新作用，提出Gk5等新候选，并强调MHC位点的免疫-代谢汇聚；这些发现支持未来机制研究，并可能为代谢和炎症疾病的治疗策略提供信息；未来应探索羧酸酯酶基因簇和Gk5等已识别遗传变异如何在分子水平影响代谢通路，以及饮食、药物或生理应激等环境因素如何与遗传易感性互作；鉴于DO小鼠运动反应存在已知遗传变异，确定代谢遗传架构在不同扰动和应激下如何变化将是重要下一步，对理解基因-环境互作和提高代谢组学遗传学的转化相关性至关重要。