高脂饮食小鼠肝脏转录组大规模荟萃分析：可复现核心、暴露时长相关变异及与MASLD相关人肝队列的方向一致性

《International Journal of Molecular Sciences》：Large-Scale Transcriptomic Meta-Analysis of Mouse Liver Responses to High-Fat Diet Identifies a Reproducible Core and Duration-Associated Variation

　　

饲喂高脂饮食（high-fat diet, HFD）的小鼠被广泛用于研究代谢功能障碍相关脂肪性肝病（metabolic dysfunction-associated steatotic liver disease, MASLD）相关的肝脏反应，但其跨实验可复现性仍不清楚。研究人员对64个小鼠肝脏RNA测序（RNA-sequencing）GEO系列（GEO Series, GSEs）进行了多水平荟萃分析（multilevel meta-analysis），其中包含105个HFD与对照的对比，并将结果与3个MASLD相关人肝队列进行比较。在18,234个小鼠基因中，6,920个显示出显著的合并HFD效应，而仅59个满足全部3项研究定义的可复现性标准：合并错误发现率（false discovery rate, FDR）< 0.05、95%预测区间（prediction interval）不包含零，以及≥80%的可评估GSEs与合并效应方向一致（GSE均衡方向一致性，GSE-balanced directional consistency）。这59个基因构成高置信度可复现亚群。25个Hallmark通路显示出显著的合并富集，HFD时长与2,744个基因水平效应及多个Hallmark通路相关。小鼠—人对应关系较弱，但在10,658个一对一直系同源基因（orthologs）中，小鼠方向一致性与人同向效应的比值比（odds ratio, OR）相关（每增加10个百分点，OR 1.075；95%置信区间（confidence interval, CI）1.035–1.115；p = 1.5 × 10^?4）。因此，肝脏HFD反应包括一个位于更广泛的时长敏感性转录反应中的小型可复现核心；同时，小鼠中较高的跨实验重复性与MASLD相关人肝队列中适度更高的方向一致性相关，但不构成直接模型验证。

研究背景与问题：MASLD是全球慢性肝病的重要原因，并具有代谢和临床异质性，包括性别相关易感性和进展差异。HFD小鼠被广泛用于研究MASLD相关肝脏代谢与分子反应，但实验在暴露时长、膳食组成、对照饲料配方以及动物性别、品系、年龄和基因型等方面差异较大，可影响肝脏表型及模型与人类疾病的相似程度。单一HFD实验中的转录改变不能自动视为HFD的普遍反应；合并效应显著也可能伴随效应大小甚至方向在实验间变化。因此，研究人员将平均HFD效应与跨实验可复现性作为互补特征，并进一步询问跨实验更一致重复的小鼠反应是否更可能在人MASLD相关肝病中同向变化。研究人员整合64个GSEs、105个HFD—对照对比、837个GEO样本记录（GEO Sample, GSM）和951个序列读取存档（Sequence Read Archive, SRR）accessions，标准化处理后对至少10个独立GSE可估计的18,234个小鼠基因进行多水平荟萃分析，并以预测区间和GSE均衡方向一致性评价可复现性。结论是：6,920个基因有显著合并HFD效应，但仅59个满足三项可复现性标准并构成高置信度可复现核心；25个Hallmark通路显著富集；暴露时长与2,744个基因水平效应及多个Hallmark通路相关；鼠—人全局对应弱，但小鼠方向一致性与人同向效应相关（每10个百分点OR 1.075，95% CI 1.035–1.115，p = 1.5 × 10^?4）。其意义在于区分小型可复现核心与更广泛的时长敏感转录反应，并提示跨实验重复性仅与MASLD相关人肝队列方向一致性适度相关，不构成直接模型验证。

技术方法概述：研究人员检索NCBI GEO DataSets（NCBI EDirect 12.0），纳入野生型或可提取野生型对照的全肝bulk RNA-seq、常规HFD与chow/标准/正常/低脂对照、每组≥3生物学重复且有公开原始数据的实验；排除并发疾病诱导、工程疾病基因型或不兼容测序配置。小鼠reads以STAR 2.7.10b比对GRCm39/GENCODE M39，featureCounts（Subread 2.0.6）计数，DESeq2 1.42.1进行对比水平差异表达；以未收缩log₂折叠变化及标准误平方为抽样方差。用metafor 4.8.0的rma.mv进行多水平随机效应荟萃分析（限制性最大似然，REML），随机效应含GSE和对比标识；计算I²、τ²、95%预测区间及GSE均衡方向一致性，并做留一GSE重拟合。调节因素分析用同结构；Hallmark分析用MSigDB 2025.1.Mm和fgsea 1.28.0。人肝队列GSE126848（31例/26对照）、GSE130970（72例/6对照）、GSE135251（206例/10对照）用edgeR 4.0.16、limma 3.58.1、TMM、voom及REML随机效应合并；以Ensembl Compara release 116一对一orthologs（置信度1）连接鼠人。

研究结果：
2.1 独立小鼠肝脏HFD转录组实验的大规模整合：通过预定义纳入框架、样本水平核对和最终accession映射，研究汇总64个GSEs、105个对比、837个GSM和951个SRR；元数据可得性有限，注释GSE中位暴露时长14周，约60%脂肪饮食为主，C57BL/6背景常见，对照包括chow、低脂饮食和正常饮食，性别以雄性为主，年龄信息不一致。标准化RNA-seq处理和对比水平差异表达后，18,234个基因至少见于10个独立GSE并纳入主荟萃分析。
2.2 广泛的平均HFD效应与更有限的跨实验可复现性形成对比：多水平荟萃分析显示18,234个基因中6,920个在FDR < 0.05下有显著合并效应，但跨基因中位总I²为65.5%；1,166个基因GSE均衡方向一致性≥0.80，其中1,042个同时FDR < 0.05；175个基因95%预测区间排除零，其中143个同时FDR < 0.05；63个同时满足两种可复现性指标（不论合并显著性）。这些指标仅部分重叠，说明统计显著性本身不能反映HFD反应跨实验重复程度。
2.3 59基因可复现核心与协调的通路水平HFD反应：59个基因同时满足FDR < 0.05、95%预测区间排除零和GSE均衡方向一致性≥0.80；留一GSE重拟合均保留合并效应方向；其中29个为正效应、30个为负效应。功能注释涉及脂质代谢（Ces1a）、细胞毒免疫反应（Gzmm）、染色质组织（H2ac8、H2bc8、H4c12）及DNA复制与维护（Zranb3、Stn1），但许多核心基因功能注释有限。核心基因中位总I²约0%，显著非核心基因为67.5%，部分由预测区间标准所致；GSE聚合单层REML与主估计高度相关（Pearson r = 0.991），59个核心基因方向与FDR显著性均保留。对全合并荟萃排名行基因集富集分析，25/50个Hallmark通路FDR < 0.05（19个正富集、6个负富集）；干扰素-γ反应标准化富集分数（normalized enrichment score, NES）为2.405（FDR = 9.31 × 10^?13），脂肪酸代谢NES为2.066（FDR = 3.00 × 10^?7）。
2.4 HFD暴露时长与数千基因水平HFD效应变异相关：评估暴露时长、膳食脂肪kcal%和对照饮食类型三个调节因素。14,421个可建模基因中2,744个显示显著时长相关效应（FDR < 0.05）；脂肪每增加10个百分点仅6/14,588个基因显著；对照饮食类型对206/13,300个基因有显著omnibus效应。加入时长使全部时长建模基因总异质性中位降低约2.7%，在显著时长相关基因中降低约21.6%。联合模型14,195个可评估基因中，时长对2,659个仍显著，脂肪对19个显著。探索性GSE水平分析中，时长与脂肪比例弱负相关（Spearman ρ = ?0.185；33个GSEs），性别与品系间Cramér’s V = 0.622（21个GSEs），但为描述性且类别不均。通路水平29/50个Hallmark通路呈显著时长相关富集；上皮—间质转化（epithelial–mesenchymal transition, EMT）正相关NES最高（NES = 2.734，FDR = 5.84 × 10^?19），胆汁酸代谢（NES = ?2.514，FDR = 3.76 × 10^?10）和外源物代谢（NES = ?2.333，FDR = 1.85 × 10^?11）负相关；脂肪校正后48/50个通路NES方向保留，主与校正NES相关（Pearson r = 0.965；Spearman ρ = 0.989），25个通路仍FDR < 0.05。
2.5 更可复现的小鼠反应在MASLD相关人肝队列中更常方向一致：分析GSE126848、GSE130970和GSE135251三个独立人肝队列；12,477个基因在人荟萃分析中可估计效应，中位τ² = 0.084、中位I² = 83.2%，人队列间异质性明显；高置信度一对一ortholog映射得到10,658个鼠人共同分析基因。全局水平上，合并小鼠HFD效应与合并人MASLD效应仅弱相关（Pearson r = 0.068；Spearman ρ = 0.045）。主logistic回归含10,658个orthologs，5,484个（51.5%）显示鼠人合并效应方向一致；小鼠GSE均衡方向一致性每增加10个百分点，方向一致性的OR增加7.5%（OR = 1.075，95% CI 1.035–1.115；p = 1.5 × 10^?4）；20个百分点差异对应拟合模型估计OR为1.16，属跨基因关联而非个体基因行为预测。校正小鼠效应量、研究覆盖、限制于≥20个小鼠GSEs及使用人固定效应估计后关联仍显著；留一人队列分析OR为1.069–1.116，均p < 0.001。GSE均衡方向一致性≥0.80的528个基因中306个（58.0%）显示人同向效应，低于阈值组为10,130个中5,178个（51.1%）；更严格定义下，≥0.80组更常出现达到人FDR < 0.05的同向效应（7.0% vs. 5.3%）及三个队列均同向（24.8% vs. 16.5%）。

讨论与结论翻译：讨论强调平均统计显著性与跨实验可复现性存在明显区别；59基因核心不代表HFD全部生物反应，核心外基因可能因实验异质性、GSE覆盖有限、技术变异或统计功效不足而未达标准。通路水平仍可见协调的免疫和代谢反应。暴露时长在所考察调节因素中与最多基因水平变异相关，且并未均匀降低所有基因异质性，其效应集中于部分HFD反应基因；时长相关通路提示更长暴露并非简单增强同一HFD反应，而是改变相对通路组成。脂肪kcal%的调节有限，但受常用约60%脂肪配方集中、未涵盖脂肪来源、碳水化合物和蛋白质组成、胆固醇添加、喂养方案、对照配方及实际能量摄入等限制，不能排除膳食混杂；对照饮食类型亦贡献系统变异但类别不均。鼠—人全局对应弱，可能与人类MASLD涉及肝细胞群和肝外器官互作、疾病异质性、遗传背景、合并症、用药及肝细胞组成差异，以及小鼠HFD方案差异有关；原文未直接评估各因素贡献。尽管全局相关弱，更可复现的小鼠反应更常在人肝队列中同向，但该关联不建立预测效用或直接验证。优势为大规模、实验多样的小鼠肝脏HFD转录组数据统一处理和多水平荟萃分析，可分别评价平均效应、新实验预期效应和跨GSE方向重复，并结合调节因素、通路和独立人肝队列。局限包括仅纳入公开GEO且有可及原始测序数据的研究，可能存在发表和数据可得性偏倚；元数据不完整，调节因素分析为跨研究关联而非因果；性别、品系、年龄和基因型因元数据不完整或类别不均未能作为额外全基因组调节因素可靠评估；雄性、C57BL/6和约60%脂肪饮食占优限制外推；异质性具体来源未能完全识别；低表达基因未收缩log₂折叠变化可能不稳定；59基因核心为研究定义的操作性标准；bulk肝脏转录组反映转录调控和细胞组成变化；人比较为三个异质MASLD相关肝队列，旨在评价方向一致性而非验证小鼠HFD反应为完整人类疾病模型，且部分人队列对照样本较小。研究结论翻译：在64个小鼠肝脏HFD转录组数据集中，59个基因满足全部三项研究定义的可复现性标准，定义了可复现的肝脏HFD反应核心。核心之外，暴露时长与数千个基因水平HFD效应变异及多个Hallmark通路变化相关。尽管鼠—人全局对应弱，跨实验更一致重复的小鼠反应在MASLD相关人肝队列中略更可能显示同向效应；该关联不构成对小鼠模型的直接验证。总体而言，这些结果区分了一小组重复出现的肝脏HFD反应与更广泛、跨实验变异的转录反应，其中暴露时长在所考察实验调节因素中关联最广。