A trans-omics gene-smoking interaction study of lung cancer based on consortium data
包含了多组学数据整合、两阶段验证（发现与复制）、分子修饰评分构建及跨组学整合分析等多个独立证据层级，并进行了敏感性分析和亚组分析，体现了较高的实验完整度。
欧洲血统人群（ILCCO-OncoArray、TRICL、PLCO、UKB） DNA甲基化、基因表达、蛋白质、代谢物分子特征 LDpred2预测的分子性状模型
发现阶段：ILCCO-OncoArray（16 606例病例 vs 13 905例对照） 复制阶段：TRICL、PLCO、UKB（共11 131例病例 vs 435 005例对照） 交互作用检测模型：logistic回归，调整年龄、性别、基因分型平台和10个遗传主成分 统计显著性标准：发现阶段q-FDR≤0.05，复制阶段P≤0.05且效应方向一致 MMS构建：基于5个弱相关生物标志物的效应权重，来自4个数据集荟萃分析
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
整合个体基因型数据和xQTL数据，预测多种分子生物标志物，用于后续交互作用分析。
整合ILCCO-OncoArray、TRICL、PLCO和UKB的个体基因型数据（27 737例病例和449 910例对照），以及来自GoDMC、eQTLGen、Ferkingstad等和Shin等的xQTL汇总数据，使用LDpred2预测DNA甲基化、基因表达、蛋白质和代谢物水平，仅保留后验遗传力≥0.01的生物标志物。
在发现和复制阶段验证分子特征与吸烟（包年）交互对肺癌风险的影响。
在ILCCO-OncoArray中测试交互作用（发现阶段），在TRICL、PLCO和UKB中验证（复制阶段）。使用logistic回归模型，调整年龄、性别、基因分型平台和前10个遗传主成分，检验分子特征与吸烟包年的交互项。
评估交互作用的稳健性，并探讨在不同亚组中的特异性。
使用吸烟状态替代包年作为环境暴露进行敏感性分析；按性别和组织学类型（LUAD和LUSC）进行亚组分析，检验交互作用在不同亚组中的一致性或特异性。
量化分子特征对吸烟效应的修饰作用，并基于MMS对吸烟者进行肺癌风险分层。
从8个显著生物标志物中选取5个弱相关（Pearson ρ≤0.7）的，以荟萃分析效应为权重构建MMS。根据MMS三分位数将人群分为低、中、高三组，分别评估吸烟包年与肺癌风险的关联，并在UKB中预测不同年龄和性别的肺癌发病率。
识别不同组学层（DNA甲基化→基因表达→蛋白质→代谢物）之间相关的生物标志物路径，并评估其与吸烟的交互作用。
计算不同组学层预测生物标志物之间的相关性，筛选q-FDR≤0.05且Pearson ρ＞0.3的显著相关对，构建从CpG位点到基因表达、蛋白质和代谢物的分子路径，并检验这些相关性是否对肺癌风险与吸烟的交互有贡献。
按研究目的归类文中使用的方法，便于定位所需技术。
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 数据整合与质量控制 | 个体基因型数据来源（ILCCO-OncoArray、TRICL、PLCO、UKB）及质量控制程序；xQTL数据来源（GoDMC、eQTLGen、Ferkingstad等、Shin等）及样本量；保留欧洲血统且遗传无关的个体。 阅读提示：Methods部分'Consortium-scale individual GWAS data'和'Alliance-based summary-level quantitative trait loci data' |
| 分子特征预测 | LDpred2参数（例如'sparse'选项）；后验遗传力阈值（h²≥0.01）；使用的xQTL汇总统计和LD参考面板。 阅读提示：Methods部分'Molecular trait prediction models' |
| 两阶段交互作用分析 | 发现阶段和复制阶段的数据集分配；logistic回归模型协变量（年龄、性别、基因分型平台、前10个遗传主成分）；交互项定义（分子特征×包年）；显著性标准（q-FDR≤0.05和P≤0.05）。 阅读提示：Methods部分'Association testing for xWAS' |
| 敏感性分析 | 吸烟状态分类（当前/曾经/从不）；是否使用与主要分析相同的协变量和显著性阈值。 阅读提示：Methods部分'Association testing for xWAS'中的敏感性分析描述 |
| 分子修饰评分构建 | 选取的生物标志物（cg13714459, cg13561554, cg26989927, cg15034267, RP11-326C3.14）；相关性排除标准（Pearson ρ>0.7）；效应权重来源（荟萃分析）；三分位数分组。 阅读提示：Methods部分'Molecular modifying score' |
| 跨组学整合分析 | 相关性阈值（q-FDR≤0.05且Pearson ρ>0.3）；起始CpG位点选择标准；构建路径的步骤。 阅读提示：Methods部分'Trans-omics integrative analysis' |
