汉中国人与非中国人群肠道微生物组差异探究

Exploring differences in the human gut microbiome between Han Chinese and non-Chinese populations

作者信息Jingjing Wang, Zhenmiao Zhang, Yang Chen, Xi Zhou, Jiaxin Xiang, Chao Yang, Dmitry A Rodionov, Andrei L Osterman, Qinwei Qiu, Yusheng Deng, Yanmin Liu, Chengrui Wang, Xiaoxiao Shang, Li Huang, Chen Sun, Jianwen Guo, Zhimin Yang, Lixiang Zhai, Zhaoxiang Bian, Aiping Lyu, Lijuan Han, Wei Jia, Xiaodong Fang, Lu Zhang
PMID41545888
发布时间2026-01-17
DOI10.1186/s13059-026-03932-3

实验完整度

高

构建大规模参考基因组目录并整合公共数据,进行队列验证,结合丰度与SNP多层面分析,包含功能富集和多样性差异验证。

主要模型

人类肠道微生物组参考基因组目录(HGMRGC) 发现队列和复制队列粪便样本 Collinsella属基因组簇 Faecalibacillus intestinalis、Escherichia coli C、Lawsonibacter asaccharolyticus等代表基因组

重点核对

地理因素(HC与NC)作为主要解释变量,并控制年龄、性别和BMI 仅使用高质量基因组(完整性>90%,污染<5%)进行SNP分析 使用发现队列(2505 HC vs 1206 NC)和复制队列(2254 HC vs 601 NC)进行验证 SNP检测采用核心基因组双等位基因SNP(MAF≥1%,位点流行率≥90%),并过滤深度<4和等位基因频率<0.01的SNP

摘要

背景:人类肠道微生物群在不同人群中表现出显著多样性,受地理、饮食和生活方式等因素影响,尤其是在汉中国人和非中国人群之间。虽然先前研究主要关注肠道微生物组的分类丰度,但单核苷酸多态性(SNP)在驱动人群特异性差异中的作用在很大程度上仍未得到充分探索。结果:在本研究中,我们利用人类肠道微生物组参考基因组目录(HGMRGC)系统地研究了汉中国人与非中国人群之间的肠道微生物差异。我们观察到地理是微生物丰度和SNP变异的主要驱动因素。我们从Collinsella属中鉴定出689个具有碳水化合物利用功能差异的人群特异性基因组簇,以及108种具有显著患病率差异的物种,这些差异与维生素生物合成、抗生素耐药性和碳水化合物代谢有关。Beta多样性分析强调了微生物丰度和SNP在人群间的显著差异,而Alpha多样性分析显示非中国人群在微生物丰度方面表现出更高的多样性,而汉中国人群在SNP方面表现出更大的多样性。结论:本研究全面分析了汉中国人与非中国人群之间的肠道微生物差异,强调了人群特异性特征对微生物多样性和功能的深远影响。我们还提供了一个全面的人类肠道微生物参考基因组目录,特别关注汉中国人群,为未来肠道微生物基因组变异研究奠定了基础。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
汉中国人与非中国人群的肠道微生物组在物种丰度和SNP水平上是否存在差异?
核心机制
地理因素(作为饮食和生活方式的代表)驱动微生物SNP和丰度变异;人群特异的碳水化合物代谢、维生素合成和抗生素耐药性功能差异体现了适应性进化。
主要证据
利用HGMRGC(包含271,480个非冗余基因组和5,785个物种簇)进行PERMANOVA、FST、alpha/beta多样性分析;在发现和复制队列中验证了地理因素的主导作用,并鉴定出689个Collinsella种群特异基因组簇和108个种群特异流行基因组。
研究意义
提供了一个包含汉中国人群的全面参考基因组目录,为研究肠道微生物基因组变异的生态和进化意义及其对宿主健康的影响奠定基础。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

构建人类肠道微生物参考基因组目录

收集并组装大量公共和新增的汉中国人群粪便宏基因组数据,以构建全面的人类肠道微生物参考基因组目录。

整合9,320个公共样本和3,584个新测序汉中国人群样本的宏基因组数据,通过组装、分箱、聚类和去冗余,得到5,785个物种级基因组簇(271,480个非冗余基因组),并评估基因组质量和分类学注释。

2

确定驱动肠道微生物变异的因素

量化地理、年龄、性别和BMI对肠道微生物丰度和SNP变异的贡献,以确定主要驱动因素。

基于PERMANOVA分析,分别以'地理-人群'(HC/NC)和'地理-国家'作为地理因素,对发现和复制队列的丰度及SNP距离矩阵进行变异分解。

3

分析人群间Beta多样性和选择压力

评估HC与NC人群间微生物丰度和SNP组成的差异程度,并通过FST识别受选择压力的基因组。

计算基于Bray-Curtis和等位基因共享评分的Beta多样性,并比较人群间与人群内差异;计算每个代表基因组的FST值,并进一步计算基因的FST值。

4

鉴定人群特异性基因组簇

识别仅存在于HC或NC人群中的基因组簇,并分析其功能差异。

根据基因组地理来源注释,将基因组簇分为HC特异和NC特异;进一步对Collinsella属的689个代表基因组进行KEGG模块和代谢表型富集分析。

5

鉴定人群特异流行基因组并分析功能

发现流行率在HC与NC人群间显著差异的物种,并探讨其功能潜力。

使用Fisher精确检验确定流行率显著差异的物种,定义HC-PSPRGs和NC-PSPRGs,然后进行KEGG模块和代谢表型富集分析。

6

比较人群内Alpha多样性

比较HC与NC人群内的微生物丰度和SNP多样性,揭示不同层面的多样性模式。

计算丰度Shannon指数和SNP位点核苷酸多样性π,分别作为丰度和SNP的alpha多样性,并在基因组水平比较特定菌种的SNP多样性。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
粪便采集试剂盒(MCK-01 KMHD)KMHDMCK-01
QIAsymphony PowerFecal Pro DNA试剂盒QIAGEN--
TruSeq DNA HT样本制备试剂盒Illumina--
HiSeq X Ten测序仪Illumina--
NovaSeq测序仪Illumina--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
参考基因组目录构建
样本组成:9,320公共样本来源和3,584新测序汉中国人群样本的地域分布
阅读提示:Methods:数据收集和宏基因组测序;Additional file 1: Table S1
基因组质量筛选
高质量基因组标准:完整性>90%且污染<5%
阅读提示:Methods:MAG基因组质量评估
发现和复制队列
队列构成:HC和NC个体数量及来源;年龄、性别、BMI分布;排除近期抗生素使用者
阅读提示:Methods:发现和复制队列;Additional file 1: Table S6
丰度计算
基因组存在标准:覆盖度>40%且RPKM>0.0001
阅读提示:Methods:计算微生物丰度和流行率
SNP calling
SNP定义:核心基因组双等位基因SNP,MAF≥1%,位点流行率≥90%;过滤条件:read depth≥4,替代等位基因频率≥0.01
阅读提示:Methods:宏基因组测序数据的SNP检测
Beta多样性分析
Beta多样性计算方法:Bray-Curtis(丰度)和1-等位基因共享评分(SNP);组间比较使用单尾Mann-Whitney U检验,随机样本对数量10,000
阅读提示:Methods:Beta多样性计算;Results:Beta多样性和人群间选择压力
FST计算
定义:FST=1-(组内核苷酸多样性/组间核苷酸多样性);分析基因组条件:流行率>10%且共享SNP位点>1000
阅读提示:Methods:FST计算
人群特异流行基因组鉴定
标准:p<0.05,以及流行率条件(一个群体<5%,另一个群体至少3倍)
阅读提示:Methods:识别HC-PSPRGs和NC-PSPRGs