肠道微生物健康指数GMWI2

作者Jaeyun Sung在2020年和2024年的两篇文章，分别提出了计算个体肠道微生物健康指数的方法:

• A predictive index for health status using species-level gut microbiome profiling
• Gut Microbiome Wellness Index 2 enhances health status prediction from gut microbiome taxonomic profiles

研究设计

作者首先收集了70各个研究中12,532份肠道微生物样本，基于一些排除标准，最终有8,869份样本保留，其中57个健康的研究有5757份样本和18个非健康研究有3112份研究。所有样本采用相同的生物信息学流程处理。菌群丰度使用MetaPhlan3计算，根据相对丰度阈值>=0.00001确定每个样本的菌群存在/不存在。

Lasso惩罚逻辑回归提取健康和疾病相关分类群

使用Lasso惩罚逻辑回归(Lasso-penalized logistic regression)对菌群的矩阵进行分析，得到:

• 59个菌(taxa)的系数是正
• 3091个菌的系数是0
• 50个菌的系数是负

非零系数的菌包括，1个冈、3个目、7科、19属、79种

系数被定义为当样本中的相关分类群从缺失变为存在时，属于健康患者的样本的log-odds的期望变化

假设我们有一组数据，记录了不同年龄段的人群是否患病的情况。通过逻辑回归分析，我们得到了关于年龄这一自变量的Odds Ratio值，比如OR=1.1。

这个OR=1.1的含义是：在控制其他可能影响患病的因素（如性别、遗传、生活习惯等）不变的情况下，年龄每增加一岁，个体患病的概率与不患病的概率之比将增加10%。换句话说，如果一个人在某个年龄段患病的概率是某个值，那么当他年龄增加一岁时，他患病的“胜率”（即患病的概率与不患病的概率之比）将变为原来的1.1倍。

根据https://github.com/biobakery/MetaPhlAn/wiki/MetaPhlAn-3.1安装MetaPhlAn的数据库，作者使用MetaPhlAn的数据库的版本是mpa_v30_CHOCOPhlAn_201901，详情见。

下载指定版本的MetaPhlAn数据库到当前文件夹的db目录:

metaphlan --install --index mpa_v30_CHOCOPhlAn_201901 --bowtie2db db

docker 命令

docker run --rm  -it \
	-u $(id -u):$(id -g) \
	-v $PWD:$PWD \
	-w $PWD \
	wybioinfo/gmwi2:1.6 \
	metaphlan --install --index mpa_v30_CHOCOPhlAn_201901 --bowtie2db db

注意此时会比较慢，MetaPhlAn需要下载数据库并建立索引

Downloading MetaPhlAn database
Please note due to the size this might take a few minutes

Downloading http://cmprod1.cibio.unitn.it/biobakery3/metaphlan_databases/mpa_v30_CHOCOPhlAn_201901.tar
Downloading file of size: 366.62 MB
366.62 MB 100.00 %  16.09 MB/sec  0 min -0 sec         
Downloading http://cmprod1.cibio.unitn.it/biobakery3/metaphlan_databases/mpa_v30_CHOCOPhlAn_201901.md5
Downloading file of size: 0.00 MB
0.01 MB 12800.00 %  31.60 MB/sec  0 min -0 sec         

Decompressing db/mpa_v30_CHOCOPhlAn_201901.fna.bz2 into db/mpa_v30_CHOCOPhlAn_201901.fna

Building Bowtie2 indexes