儿童肥胖谱系中的代谢重塑与肠道菌群网络重构：整合多组学分析

《Gut Microbes》：Integrated multi-omics analysis reveals distinct metabolic signatures and gut microbiome alterations across the pediatric obesity spectrum

　　

儿童肥胖驱动肠道和全身隔室广泛的代谢重塑，但覆盖完整肥胖谱系的整合多组学表征仍不完整。研究人员研究了57名墨西哥儿童，按体重指数（BMI）分为正常体重（NW，n = 15）、超重（Ow，n = 8）和肥胖（Ob，n = 34）组，采用肝酶谱、非靶向粪便和血清代谢组学、16S rRNA肠道微生物组测序及抗菌药物耐药性评估。各组人口学参数可比（p = 0.632），减少了混杂。γ-谷氨酰转移酶（GGT）与肥胖关联最强（25.5 ± 6.0 vs. NW的15.0 ± 3.0 U/L；p < 0.001），作为早期肝脏应激生物标志物优于丙氨酸氨基转移酶（ALT）和天冬氨酸氨基转移酶（AST）。代谢组学分析鉴定出1088种改变的粪便代谢物和1041种改变的血清代谢物（Ob vs. NW），富集通路涵盖芳香族氨基酸代谢、一碳代谢、糖酵解/三羧酸（TCA）循环和核苷酸生物合成。关键扰动包括犬尿氨酸和3-硝基酪氨酸上调，同时胆汁酸和嘌呤代谢物下调。16S rRNA测序显示，在Ow vs. NW、Ob vs. Ow和Ob vs. NW比较中分别有78、95和28个差异丰度分类单元，肥胖中β多样性离散度升高。共丰度网络分析揭示代谢物-代谢物和分类单元-分类单元相互作用受损，包括有益的产生短链脂肪酸（SCFA）的菌群联合体丧失、胆汁酸失调和氨基酸-微生物组耦合改变。微生物组与代谢组数据整合揭示高幅度属-代谢物相关性——包括Lachnospiraceae与乙酸盐（acetate）（ρ = 0.92）以及Tyzzerella与vulpecholate（ρ = 0.94）——映射到将SCFAs、胆汁酸、色氨酸-吲哚衍生物和芳香族氨基酸与肥胖发病机制联系起来的机制通路。这些发现确定GGT为敏感的儿童代谢生物标志物，并为早期干预提供机制见解和治疗靶点。

儿童肥胖是重大全球健康挑战，影响超过3.4亿儿童青少年。肥胖可驱动肠道与全身代谢重塑，肠道微生物组被视为代谢健康的核心介导者。已有研究揭示肥胖儿童肠道中产生SCFA的Faecalibacterium、Roseburia、Agathobacter和Bacteroides fragilis等变化，并发现氨基酸代谢、脂质生物合成、芳香族氨基酸分解和能量代谢相关微生物通路增强；代谢组学也提示氨基酸、糖异生、脂肪酸合成、β-丙氨酸和尿素循环等通路扰动，一碳代谢、支链氨基酸（BCAA）、芳香族氨基酸和犬尿氨酸通路代谢物可作为早期生物标志物。然而，多数研究聚焦孤立肥胖表型或有限代谢隔室，覆盖从正常体重到超重再到肥胖全谱系的整合多组学分析仍少，尤其在拉丁美洲人群中；基于网络的代谢物-代谢物和分类单元-分类单元共丰度分析在儿童肥胖中应用不足。为此，研究人员在墨西哥儿童中开展整合多组学分析，结合肝酶谱、非靶向粪便和血清代谢组学、16S rRNA微生物组表征、抗菌药物耐药性评估和共丰度网络分析，以刻画人口学与临床特征、识别肠道和全身代谢通路扰动、定义微生物组组成与多样性模式、绘制代谢物-代谢物和微生物组-微生物组互作网络，并建立早期代谢功能障碍的整合生物标志物特征。该论文发表于《Gut Microbes》。

关键技术方法包括：研究队列为2020—2022年墨西哥城6—18岁儿童，最终57名具有完整人口学、生化、代谢组和微生物组资料者纳入分析，按WHO生长标准BMI百分位分层为NW（n=15）、Ow（n=8）和Ob（n=34；肥胖I类18、II类11、III类5）。检测空腹血ALT、AST、GGT；采用超高效液相色谱-高分辨质谱（UPLC-HRMS）进行非靶向粪便和血清代谢组学；扩增16S rRNA基因V3-V4区并在Illumina MiSeq测序；开展抗菌药物耐药性评估；构建代谢物和属水平微生物共丰度网络，并以Spearman相关分析属-代谢物关联。

3.1 人口学特征与混杂因素评估（Demographic characteristics and confounding factor assessment biomarker）：最终纳入57名儿童，NW 15例、Ow 8例、Ob 34例（I类18、II类11、III类5）；年龄（p=0.632）和性别（p=0.54）在体重状态组间可比，提示二者不是后续代谢和微生物组分析的重要混杂因素。3.2 肝酶谱显示GGT为敏感肥胖标志（Liver enzyme profiles reveal GGT as a sensitive obesity）：GGT在NW为15.0 ± 3.0 U/L、Ob为25.5 ± 6.0 U/L（p<0.001，效应量d=2.1），Ob vs. Ow p=0.008，Ow vs. NW p=0.18；ALT在Ob为24.8 ± 8.1 U/L、NW为18.2 ± 4.5 U/L（p=0.012）；AST组间稳定（p=0.324）；GGT/ALT比随肥胖严重度升高（0.82、0.85、0.90、0.93、0.92），提示早期肝脏氧化应激先于明显转氨酶升高。3.3 代谢组学揭示广泛通路重塑（Metabolomic profiling reveals extensive pathway remodeling）。3.3.1 粪便代谢组改变（Fecal metabolome alterations）：非靶向粪便代谢组学鉴定3245个特征，1088个差异代谢物；Ob vs. NW 687个、Ob vs. Ow 412个、Ow vs. NW 289个；PLS-DA显示Ob与NW分离（R²=0.68，Q²=0.52）。23条KEGG通路显著改变，包括芳香族氨基酸代谢（p=1.2×10^?8；犬尿氨酸FC=3.2、3-羟基犬尿氨酸FC=2.8、喹啉酸FC=2.1升高）、一碳代谢（p=3.5×10^?7；5-甲基四氢叶酸、5,10-亚甲基四氢叶酸、四氢叶酸下调）、能量代谢（p=8.2×10^?6；乳酸、丙酮酸、琥珀酸升高）、核苷酸代谢（p=2.1×10^?5；次黄嘌呤、黄嘌呤、尿酸降低）和胆汁酸代谢（p=4.8×10^?5；平均FC=0.58）。3.4 血清代谢组改变（Serum metabolome alterations）：血清代谢组学鉴定2834个特征，1041个差异代谢物；Ob vs. NW 612个、Ob vs. Ow 398个、Ow vs. NW 231个；PLS-DA显示稳健分组（R²=0.71，Q²=0.58）。苯丙氨酸、酪氨酸、色氨酸升高；犬尿氨酸FC=2.9，K/T比0.042 vs. 0.028（p<0.001）；3-硝基酪氨酸FC=3.4，丙二醛和4-羟基壬烯醛升高；血清乙酸盐、丙酸盐、丁酸盐轻度升高；亮氨酸、异亮氨酸、缬氨酸升高；SAM降低、SAH升高、SAM/SAH比1.2 vs. 2.8（p<0.001）。3.5 肠道微生物组组成与多样性模式（Gut microbiome composition and diversity patterns）。3.5.1 α和β多样性（Alpha and beta diversity）：16S rRNA测序获420万条高质量reads、856个ASVs；Shannon指数无显著差异（NW 4.65 ± 0.42、Ow 4.42 ± 0.51、Ob 4.55 ± 0.48；p=0.38）；β多样性组间差异显著（PERMANOVA R²=0.12，p=0.001），Ob组β多样性离散度高于NW（PERMDISP F=8.24，p=0.002）；PCoA显示Ob与NW沿PC1部分分离（解释18.2%方差）。3.5.2 微生物分类单元差异丰度（Differential abundance of microbial taxa）：DESeq2在Ow vs. NW、Ob vs. NW、Ob vs. Ow中分别鉴定78、28、95个差异ASVs；F/B比在Ob为1.8、NW为1.2（p=0.018）；Ob中Faecalibacterium、Roseburia、Agathobacter、B. fragilis、Akkermansia耗竭，Blautia、Dorea、R. gnavus、Collinsella、Escherichia-Shigella富集；科水平Lachnospiraceae、Ruminococcaceae降低，Enterobacteriaceae、Peptostreptococcaceae升高。3.6 共丰度网络分析揭示代谢与微生物互作破坏（Co-abundance network analysis reveals disrupted metabolic and microbial interactions in co-abundance networks of metabolites and microbial taxa across weight-status groups）：粪便代谢物网络中NW为50节点142边、密度0.116、模块度0.58；Ob为50节点89边、密度0.073（较NW降37%）、模块度0.38、8个弱连接组分、聚类系数0.28（降33%）；犬尿氨酸度中心性由12降至5；胆汁酸-SCFA显著边由8降至2；NW存在而Ob缺失68条边（48%），包括叶酸-氨基酸12条、胆汁酸-SCFA 8条、TCA-核苷酸10条；Ob新增23条边主要涉及氧化应激和促炎代谢物。3.6.1 微生物组共丰度网络（Microbiome co-abundance networks）：NW为42节点156边、密度0.181、模块度0.62，SCFA产生菌群模块含28个正相关，Faecalibacterium和Bacteroides介数中心性分别为0.24和0.21；Ob为42节点98边（降37%）、密度0.114、模块度0.44，Faecalibacterium介数中心性降至0.09（降62%），出现Escherichia-Shigella-R. gnavus簇12个正相关，负相关18 vs. 8；Ow网络密度0.148、模块度0.54；NW存在而Ob丢失72条边（46%），包括SCFA菌群内18个正相关、Bacteroides与丁酸产生菌14条、Akkermansia相关10条；Ob新增28条边涉及机会性 taxa。3.6.2 跨域整合：微生物组-代谢物相关性（Cross-domain integration: microbiome–metabolite correlations）：NW中Faecalibacterium-丁酸盐（ρ=0.68，p=0.002）、Roseburia-丙酸盐（ρ=0.62，p=0.005）、Bacteroides-次级胆汁酸（ρ=0.58，p=0.008）、Akkermansia-黏蛋白来源代谢物（ρ=0.54，p=0.012）呈正相关；Ob中减弱或丧失（0.28，p=0.18，降59%；0.31，p=0.15，降50%；0.22，p=0.28，降62%）；Ob新增Escherichia-Shigella-犬尿氨酸（ρ=0.71，p=0.001）、R. gnavus-3-硝基酪氨酸（ρ=0.64，p=0.003）、Collinsella-BCAA（ρ=0.59，p=0.006）；高幅度属-代谢物相关包括Lachnospiraceae-乙酸盐（ρ=0.92）和Tyzzerella-vulpecholate（ρ=0.94）。

讨论部分指出，该研究样本量中等，超重亚组仅8例，限制部分比较的统计效能；横断面设计不能推断因果；缺乏详细饮食、体力活动和社会经济数据；16S rRNA分辨率有限且无直接功能信息。讨论强调GGT较传统转氨酶更敏感，可能反映早期氧化应激和膜损伤先于肝细胞损伤；粪便和血清代谢组显示芳香族氨基酸/犬尿氨酸通路、一碳代谢、能量代谢、核苷酸和胆汁酸通路协调重塑；微生物组表现为有益SCFA产生菌耗竭和促炎 taxa 富集，且网络分析显示代谢网络碎片化、微生物网络重连、有益微生物-代谢物耦合丧失和病理互作模块出现；这些发现支持早期筛查、靶向代谢干预、微生物组靶向干预和基于个体微生物组/代谢组特征的个性化策略，但需更大队列、纵向和功能组学验证。

研究结论翻译：这项整合多组学研究全面表征了墨西哥儿童肥胖谱系中的代谢和微生物组改变。研究人员证明：（1）人口学因素（年龄、性别）在体重状态组间可比，减少了混杂效应；（2）GGT是肝脏氧化应激的敏感早期生物标志物，早于明显转氨酶升高；（3）肥胖诱导肠道和全身隔室中芳香族氨基酸、一碳、能量、核苷酸和胆汁酸通路的广泛代谢重塑；（4）肠道微生物组显示有益SCFA产生菌群耗竭和促炎 taxa 富集；（5）共丰度网络分析揭示肥胖中代谢网络碎片化、有益微生物菌群丧失、微生物-代谢物耦合受损和病理互作模块出现；（6）在该队列中，抗菌药物耐药模式似乎与年龄的关联强于与BMI状态的关联。肠道微生物组和粪便代谢组的多组学整合表明，肥胖和超重由全面代谢重编程驱动，而非单纯热量过剩。分析报告中检测到的相关性——涵盖产乙酸Lachnospiraceae、修饰胆汁酸的Tyzzerella以及代谢氨基酸的Turicibacter和Barnesiella——得到机制文献的深度佐证。菌群失调迫使微生物群改变输出：有益SCFAs和保护性吲哚衍生物耗竭，而脂肪生成底物、 overriding 法尼醇X受体（FXR）稳态的次级胆汁酸以及炎性多胺/硫化物积累。这种肠源代谢组转变破坏肠道屏障，触发代谢性内毒素血症，并决定宿主脂质和葡萄糖病理。尽管精确的机器学习特征需要在更大、平衡且采用靶向液相色谱-串联质谱（LC-MS/MS）的队列中验证，生物学映射突出表明，靶向微生物组调节（通过益生元、后生元或靶向细菌补充如4HPAA或Barnesiella）在对抗全球代谢性疾病流行方面具有巨大治疗潜力。这些发现建立了儿童肥胖的核心代谢和微生物组特征，为肥胖相关代谢功能障碍提供机制见解，并确定早期干预的潜在治疗靶点。基于网络的方法展示了仅靠差异丰度分析无法捕捉的菌群失调涌现属性，突出系统层面思维在理解复杂代谢疾病中的价值。未来需要更大队列的纵向研究、更深多组学分析和干预试验来验证这些发现并转化为临床实践。