CHARLS| 这才是真正的双因素对死亡影响的分析，交互后再亚组

欢迎参加孟德尔随机化课程(11.18-11.19)

2023年10月4日，暨南大学第二临床医学院学者在《Frontiers in Endocrinology》（二区，IF=5.2）发表题为：" The relationship between sarcopenia and mortality in Chinese community-dwelling adults: a 7-year cohort study with propensity score matching and Mendelian randomization" 的研究论文。

这项研究是基于来自中国健康与养老追踪调查（CHARLS）中的数据。研究旨在评估中国社区人群中肌肉减少症的7年死亡风险，并探索肌肉减少症的组成部分与任何死亡之间的因果关系。研究结果表明，中国人群中肌减少症及其组成对死亡风险的显著影响。特别是，低握力和通常的行走速度是死亡风险的显著因素。

这个研究主题并不新颖，同主题同一数据库的文章也有发表，那么这篇文章是如何创新的呢？

摘要与主要结果

一、摘要

背景：肌肉减少症与不良健康结果相关，包括死亡风险增加。本研究旨在评估中国社区人群中肌肉减少症的7年死亡风险，并探索肌肉减少症的组成部分与任何死亡之间的因果关系。

方法：数据来源于2011年至2018年间进行的中国健康与养老追踪调查(CHARLS)。使用亚洲少肌症工作组(AWGS) 2019标准诊断少肌症。使用逻辑回归、Kaplan-Meier (KM)生存分析和倾向评分匹配及逆概率加权。使用欧洲人口数据进行的孟德尔随机化(MR)分析被用于评估少肌症和任何死亡之间的因果关系。

结果：该研究包括9006名参与者:3892名无肌肉减少症，3570名有可能肌肉减少症，1125名有肌肉减少症，419名有严重肌肉减少症。在7年的随访中，有871例死亡，包括196例肌肉减少症和133例严重肌肉减少症。KM曲线显示肌肉减少症有较高的死亡风险。在校正混杂变量后，与肌肉减少症患者相比，肌肉减少症患者7年死亡率的优势比(ORs)为1.41 (95% CI，1.06-1.87)(p < 0.05)。严重肌肉减少症的or为2.11 (95% CI为1.51-2.95)。倾向评分匹配分析和治疗权重分析的逆概率证实了这些发现。肌肉减少症和7年死亡率的校正or在45-60岁组为2.94 (95% CI，1.6-5.39)，在60-80岁组为1.72 (95% CI，1.11-2.68)，在≥80岁组为5.03 (95% CI，0.48-52.65)。严重肌肉减少症和7年死亡率的or在45-60岁年龄组为6.92 (95% CI，1.95-24.5)，在60-80岁年龄组为2.59 (95% CI，1.61-4.17)，在≥80岁年龄组为12.52 (95% CI，1.18-133.18)。利用负方差加权(IVW)方法进行的磁共振分析揭示了60岁及以上老年人握力低、通常走路速度和死亡风险之间的因果联系。

结论：这项研究强调了中国人群中肌减少症及其组成对死亡风险的显著影响。特别是，低握力和通常的行走速度是死亡风险的显著因素。

二、研究结果

1.基线特征

该研究包括9006名参与者，其中45.9%为男性，54.1%为女性。其中，无肌肉减少症3892例，可能有肌肉减少症3570例，有肌肉减少症1125例，重度肌肉减少症419例。在7年的随访中，有871例死亡，其中147例无肌肉减少症，395例可能有肌肉减少症，196例有肌肉减少症，133例有严重肌肉减少症。

2.logistic回归分析

单因素和多因素logistic回归分析显示，7年死亡率的危险因素为肌肉减少症（OR：1.41,95%CI，1.06-1.87）、重度肌少症（OR=2.11,95%CI，1.51-2.95）、年龄（OR=1.09,95%CI，1.08-1.1）、糖尿病（OR=1.51,95%CI，1.26-1.82）、高血压（OR=1.31,95%CI，1.11-1.56）、 癌症（OR = 4.46,95% CI，2.45-8.09）、卒中（OR = 1.64,95% CI，1.13-2.38）、肺部疾病（OR = 1.82,95% CI，1.45-2.28）和记忆相关疾病（OR = 2.11,95% CI，1.37-3.25）（p < 0.05）。7年死亡率的保护因素为女性（OR=0.52,95%CI，0.43-0.63）、已婚（OR=0.33,95%CI，0.17-0.61）和离婚（OR=0.46,95%CI，0.24-0.9）（p < 0.05）。

3.肌肉减少症与全因死亡率的关系

根据 Kaplan-Meier 生存曲线，肌肉减少症和重度肌肉减少症的 7 年死亡率风险更高。在调整年龄、体重指数（BMI）、教育水平、婚姻状况、城市地区、饮酒、吸烟、糖尿病、高血压、癌症、心脏病、中风、肺病、关节病、肝病、肾病、消化系统疾病、哮喘和记忆相关疾病（模型4）后，肌少症和重度肌少症对7年死亡率的OR分别为1.41（95%CI，1.06-1.87）和2.11（95%CI， 1.51–2.95)。7年死亡率中可能肌肉减少症的未调整OR为3.17（2.61-3.85），差异有统计学意义（p < 0.05），但调整后的OR为1.11-1.23，无统计学差异（p > 0.05）

进一步按年龄分层进行亚组分析

4. PSM 和 IPTW 分析

将肌肉减少症和重度肌肉减少症作为单独的组进行倾向评分匹配。在肌肉减少症和重度肌肉减少症中，PSM 和 IPTW 分析显示 7 年死亡率的 OR 显著升高 （p < 0.05）。倾向评分匹配后，肌少症的 7 年死亡率 OR 为 1.33 （1.09-1.63），重度肌肉减少症的 OR 为 1.7 （1.25-2.32）。

经逆概率加权分析后，肌肉减少症和重度肌少症的7年死亡率OR分别为1.77（1.49-2.09）和2.28（1.69-3.07）

5. 孟德尔随机化研究

使用 MR 方法分析了 283 个 SNP，我们发现 MR-Egger 方法有显著结果 （OR = 1.261,95% CI，1.056–1.507，p = 0.011）。加权中位数法的OR为0.948（95%CI，0.850-1.058，p = 0.341），而IVW方法的OR为0.954（95%CI，0.892-1.021，p = 0.176 ）。评估了异质性 （p > 0.05） 和潜在的多效性 （p < 0.05）

设计与统计学方法

一、研究设计

P（Population）参与者：

队列研究：2011-2018年中国健康与养老追踪调查（CHARLS）中有评估肌肉减少症指标的人

孟德尔随机化研究：使用欧洲人群进行

E（exposure）暴露因素：肌肉减少症

O（outcome）结局：死亡

S（Study design）研究类型：队列研究

二、统计方法

1.队列研究

1.1.统计描述+差异性分析：对于连续变量，提供了95%的置信区间（CIs），而对于分类变量则提供了百分比频率。为了比较连续数据和分类数据，我们使用了t检验和χ2检验。

1.2.通过使用logistic回归模型，可以确定死亡风险。Kaplan-Meier 曲线以视觉方式描绘。为了减少潜在的选择偏倚，倾向评分匹配 （PSM） 用于平衡有和没有肌肉减少症或严重肌肉减少症的参与者之间的协变量。使用逻辑回归模型计算个体倾向评分后，使用卡尺宽度为倾向评分 0.2 个标准差的最近邻匹配技术将最低握力 （HGS） 组患者与其他组患者配对。

1.3.使用逆概率加权 （IPTW） 进行回归分析

1.4.使用 R 软件包和 Free Statistics 软件版本 1.7 进行统计分析。统计显著性由双侧 p 值 <0.05 确定。

2孟德尔随机化研究

2.1.采用逆方差加权（IVW）方法进行孟德尔随机化（MR）分析

2.2.SNP选择过程涉及多种方法：显著性阈值为p < 5 × 10−8用于鉴定SNP，实现全基因组意义。阈值为 r2< 0.001（聚集距离为 10,000 kb）被设置为排除处于连锁不平衡状态的 SNP。

2.3.在解决潜在的多效性问题时，PhenoScanner数据库被用于SNP鉴定和随后的排除。如F统计量<10所示，弱工具变量被系统地排除在分析之外。

2.4.在每次 MR 分析之前，采用 MR–多效性残差和异常值 （MR-PRESSO） 方法以消除潜在的异常值。

2.5.为了验证这些发现，包括MR-Egger、加权中位数和加权模式在内的替代方法与IVW一起应用。使用 Cochran 的 Q 检验评估异质性，而通过 MR-Egger 截距检验和留一法分析评估多效性。

2.6.所有统计分析均使用R软件（4.3.0版）进行MR分析和敏感性评估。

小结

这篇文章主题很早就已经发过，对于已经被发表过的主题如何进行创新：

1、扩大/缩小人群范围，在某一特定人群中研究

2、统计方法创新：例如倾向得分加权、倾向得分匹配、中介分析、交互作用分析等等

3、联合孟德尔随机化研究

4、联合多个数据库

另外，现在越来越多的文章做双因素对疾病的影响，包括遗传和环境，遗传和行为，环境和行为等，但很多都是进行了亚组分析或把两个因素两两组合生成一个新变量，但真正探讨双因素作用的还是得用交互分析，存在交互后再进行亚组分析。

一个专门做公共数据库的公众号，关注我们