构建人类肠道微生物参考基因组目录
收集并组装大量公共和新增的汉中国人群粪便宏基因组数据,以构建全面的人类肠道微生物参考基因组目录。
整合9,320个公共样本和3,584个新测序汉中国人群样本的宏基因组数据,通过组装、分箱、聚类和去冗余,得到5,785个物种级基因组簇(271,480个非冗余基因组),并评估基因组质量和分类学注释。
Exploring differences in the human gut microbiome between Han Chinese and non-Chinese populations
构建大规模参考基因组目录并整合公共数据,进行队列验证,结合丰度与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
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
收集并组装大量公共和新增的汉中国人群粪便宏基因组数据,以构建全面的人类肠道微生物参考基因组目录。
整合9,320个公共样本和3,584个新测序汉中国人群样本的宏基因组数据,通过组装、分箱、聚类和去冗余,得到5,785个物种级基因组簇(271,480个非冗余基因组),并评估基因组质量和分类学注释。
量化地理、年龄、性别和BMI对肠道微生物丰度和SNP变异的贡献,以确定主要驱动因素。
基于PERMANOVA分析,分别以'地理-人群'(HC/NC)和'地理-国家'作为地理因素,对发现和复制队列的丰度及SNP距离矩阵进行变异分解。
评估HC与NC人群间微生物丰度和SNP组成的差异程度,并通过FST识别受选择压力的基因组。
计算基于Bray-Curtis和等位基因共享评分的Beta多样性,并比较人群间与人群内差异;计算每个代表基因组的FST值,并进一步计算基因的FST值。
识别仅存在于HC或NC人群中的基因组簇,并分析其功能差异。
根据基因组地理来源注释,将基因组簇分为HC特异和NC特异;进一步对Collinsella属的689个代表基因组进行KEGG模块和代谢表型富集分析。
发现流行率在HC与NC人群间显著差异的物种,并探讨其功能潜力。
使用Fisher精确检验确定流行率显著差异的物种,定义HC-PSPRGs和NC-PSPRGs,然后进行KEGG模块和代谢表型富集分析。
比较HC与NC人群内的微生物丰度和SNP多样性,揭示不同层面的多样性模式。
计算丰度Shannon指数和SNP位点核苷酸多样性π,分别作为丰度和SNP的alpha多样性,并在基因组水平比较特定菌种的SNP多样性。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 样本采集 | 粪便采集试剂盒(MCK-01 KMHD) | KMHD | MCK-01 |
| DNA提取 | 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 |