关于欧洲人群中系统性疾病与圆锥角膜因果关联的孟德尔随机化研究报告

一、 研究作者、机构与发表信息

本研究的主要作者为Xiaxue Chen和Lanlan Chen。Xiaxue Chen所属机构为中国吉林省长春市吉林大学第二医院眼科；Lanlan Chen所属机构为中国吉林省长春市吉林大学第一医院肝胆胰外科、普外科中心。该研究发表于American Journal of Ophthalmology期刊，于2024年4月30日被接受，正式发表于2024年，卷号为265，页码为189-199。

二、 学术背景与研究目的

本研究属于眼科学与遗传流行病学的交叉领域，具体聚焦于探讨多种常见系统性疾病与圆锥角膜（Keratoconus， KC）之间是否存在因果关联。圆锥角膜是一种以角膜非炎性变薄、不对称性锥形前凸为特征的角膜病变，常导致进行性视力损害，晚期可能需要进行角膜移植，给患者和社会带来沉重负担。尽管其病因被认为是遗传和环境因素复杂互作的结果，但具体机制尚不完全明确。

长期以来，临床观察和流行病学研究提示圆锥角膜可能与多种系统性疾病共存，例如花粉症/过敏性鼻炎/湿疹、糖尿病、阻塞性睡眠呼吸暂停、甲状腺功能障碍等。然而，这些关联性研究的结果常常相互矛盾，且多为横断面或病例对照研究，难以区分因果关系与混杂因素（如共同的环境因素、反向因果关系等）的影响。例如，关于糖尿病对圆锥角膜是起保护作用还是风险因素，文献中存在截然不同的观点。因此，厘清这些系统性疾病与圆锥角膜之间的因果关系，对于理解圆锥角膜的病因学、识别高危人群以及制定临床干预策略具有重要的科学和临床意义。

在此背景下，本研究旨在利用孟德尔随机化（Mendelian Randomization， MR） 这一遗传流行病学方法。MR利用与暴露因素（此处指各种系统性疾病或生理指标）强相关的遗传变异作为工具变量（Instrumental Variables， IVs），来推断暴露与结局（此处指圆锥角膜）之间的潜在因果关系。由于基因型在配子形成时随机分配，且通常不受后天环境因素影响，MR方法能够有效减少传统观察性研究中的混杂偏倚，并有助于评估因果方向。

本研究的具体目的是：在欧洲人群中，运用两样本MR分析，系统性地评估一系列常见系统性疾病（包括糖尿病、花粉症/过敏性鼻炎/湿疹、阻塞性睡眠呼吸暂停、甲状腺功能障碍、主动脉瘤、重度抑郁症、炎症性肠病、二尖瓣脱垂）与圆锥角膜发病风险之间的因果关联。

三、 详细研究流程与方法

本研究是一项基于公开数据的两样本孟德尔随机化分析，其工作流程严格遵循STROBE-MR报告规范，并基于MR的三个核心假设：1) 工具变量与暴露强相关；2) 工具变量与混杂因素无关；3) 工具变量仅通过暴露影响结局。

流程一：数据来源与工具变量筛选 研究者从多个公共全基因组关联研究（Genome-Wide Association Studies， GWAS）数据库中获取了暴露（系统性疾病）和结局（圆锥角膜）的遗传数据。所有原始GWAS均已获得相应伦理委员会批准，参与者均知情同意。 1. 结局数据：圆锥角膜的GWAS数据来源于一项包括2116名欧洲病例和24,626名欧洲对照的研究（GWAS Catalog ID: gcst90013442）。 2. 暴露数据：针对每种系统性疾病，研究者从特定GWAS中筛选与之强相关的单核苷酸多态性（Single Nucleotide Polymorphisms， SNPs）作为工具变量。 * 糖尿病：为避免单纯诊断标签的局限性，研究者采用了一项近期MR研究中的方法，选取了40个与2型糖尿病相关且同时与糖化血红蛋白（HbA1c）水平相关的SNPs，以代表长期平均血糖水平。 * 花粉症/过敏性鼻炎/湿疹：数据来自MRC-IEU，包含106,201例病例和355,812例对照，筛选出151个SNPs。 * 阻塞性睡眠呼吸暂停：基于一项多特征GWAS荟萃分析，筛选出29个SNPs。 * 甲状腺功能：使用甲状腺刺激激素（TSH）和游离甲状腺素（FT4）作为连续变量指标，分别筛选出57个和29个SNPs。 * 主动脉瘤、重度抑郁症、克罗恩病、溃疡性结肠炎、二尖瓣脱垂：均从相应的大型GWAS中筛选SNPs，显著性阈值多为p < 5×10^-8，并基于连锁不平衡（r^2 < 0.001或0.01）进行聚类。 为确保工具变量的强度，研究者计算了F统计量，并剔除了F < 10的弱工具变量。最终用于分析的SNPs数量从29到151个不等，所有SNPs的平均F统计量均大于10，表明工具变量强度足够。

流程二：孟德尔随机化分析与敏感性检验 研究采用多种MR模型进行分析，以确保结果的稳健性。 1. 主要分析方法：逆方差加权法（Inverse-Variance Weighted， IVW），特别是随机效应IVW模型，作为推断因果关联的主要方法。 2. 补充分析方法：同时使用MR-Egger回归和加权中位数法进行估计，以在不同假设下验证结果。 3. 敏感性分析：为评估结果是否受水平多效性（即SNP通过非暴露途径影响结局）或异质性的影响，研究进行了以下检验： * 异质性检验：使用Cochran‘s Q检验。 * 水平多效性检验：使用MR-Egger截距检验。 * 离群值检测与校正：使用MR-PRESSO方法检测并剔除可能具有水平多效性的离群SNPs，然后重新计算效应估计值。 * 留一法分析：依次剔除每个SNP后重复MR分析，以检查是否有单个SNP对结果产生过度影响。 4. 进一步验证：对于IVW方法显示有显著因果关联的暴露，研究者额外使用了CAUSE（Causal Analysis Using Summary Effect Estimates）模型进行验证，该模型能更好地处理复杂的多效性结构，降低假阳性结论的风险。 5. 多重检验校正：使用错误发现率（False Discovery Rate， FDR）方法对p值进行校正。 所有统计分析均使用R软件（4.3.0版）进行，主要利用TwoSampleMR、MR-PRESSO和CAUSE等软件包。

四、 主要研究结果

研究的主要发现基于随机效应IVW模型的结果，其他方法作为补充和验证。

1. 具有显著正因果关联的疾病：

   • 花粉症/过敏性鼻炎/湿疹：遗传预测的花粉症/过敏性鼻炎/湿疹与圆锥角膜风险显著正相关（比值比[OR] = 10.144；95%置信区间[CI]， 2.441-42.149；p = 0.001）。这意味着，遗传上易患这些过敏性疾病的人群，患圆锥角膜的风险显著升高。这一结果在CAUSE模型中也得到强力支持（p = 4.5×10^-9）。
   • 溃疡性结肠炎：遗传预测的溃疡性结肠炎与圆锥角膜风险显著正相关（OR = 1.147；95% CI， 1.054-1.248；p = 0.002）。CAUSE模型同样验证了这一关联（p = 4.8×10^-5）。

2. 未发现显著因果关联的疾病/指标：

   • 长期高血糖状态（以HbA1c相关遗传工具代表）：与圆锥角膜风险无显著因果关联（OR = 0.320；95% CI， 0.029-3.549；p = 0.353）。这表明，遗传预测的较高长期血糖水平并未显示出对圆锥角膜具有潜在的保护作用，推翻了之前一些观察性研究的假设。
   • 其他疾病：阻塞性睡眠呼吸暂停（OR = 0.775；p = 0.580）、甲状腺刺激激素TSH（OR = 0.900；p = 0.169）、游离甲状腺素FT4（OR = 0.682；p = 0.022，但经FDR校正后不显著）、主动脉瘤（OR = 0.995；p = 0.937）、重度抑郁症（OR = 1.227；p = 0.627）、克罗恩病（OR = 1.096；p = 0.062）、二尖瓣脱垂（OR = 0.980；p = 0.826）均未显示出与圆锥角膜有显著的因果关联（所有比较p > 0.05）。

3. 敏感性分析结果：

   • 在HbA1c的分析中检测到异质性（Cochran‘s Q检验p = 0.030），但由于采用了随机效应IVW模型，且结果本身不显著，因此不影响“无因果关联”的结论。
   • 在其他大多数分析中，未检测到显著的异质性或水平多效性，表明MR分析的核心假设得到较好满足，结果相对稳健。
   • 留一法分析未发现任何单个SNP对上述显著结果（花粉症/过敏性疾病、溃疡性结肠炎）或主要阴性结果产生决定性影响。

五、 研究结论与意义

本研究通过大规模的孟德尔随机化分析，为欧洲人群中系统性疾病与圆锥角膜的因果关联提供了新的遗传学证据。

主要结论： 1. 花粉症/过敏性鼻炎/湿疹和溃疡性结肠炎是圆锥角膜的潜在因果风险因素。遗传易感性导致个体罹患这些疾病的同时，也显著增加了患圆锥角膜的风险。 2. 长期高血糖状态（糖尿病） 并未显示出对圆锥角膜具有遗传学上的保护作用。之前一些观察性研究中发现的“保护性关联”很可能受到混杂因素或反向因果关系的影响。 3. 在本研究分析的样本量和遗传工具条件下，未发现阻塞性睡眠呼吸暂停、甲状腺功能障碍、主动脉瘤、重度抑郁症、克罗恩病和二尖瓣脱垂与圆锥角膜存在因果关联。

科学价值与应用价值： * 病因学启示：该研究首次利用MR方法系统评估了多种系统性疾病与圆锥角膜的因果关系，为圆锥角膜的病因学研究提供了更高级别的证据。特别是明确了过敏性疾病和溃疡性结肠炎的风险角色，提示免疫和炎症通路可能在圆锥角膜发病中起重要作用（例如，过敏患者可能因眼痒而增加揉眼行为，而揉眼是已知的圆锥角膜风险因素；溃疡性结肠炎相关的系统性炎症可能影响角膜稳态）。 * 临床实践指导：研究结果对临床具有指导意义。眼科医生在面对患有严重花粉症、过敏性鼻炎、湿疹或溃疡性结肠炎的患者时，应提高对圆锥角膜的警惕性，考虑进行早期筛查（如角膜地形图检查）。同时，对于合并糖尿病的圆锥角膜患者，不应因过去某些研究中所谓的“保护作用”而延误对圆锥角膜的诊断和治疗，应积极控制血糖并处理角膜问题。 * 研究方法示范：研究展示了在复杂疾病病因探索中，如何严谨地应用MR方法，包括精心选择工具变量（如使用代表连续生理指标的遗传工具）、采用多种互补的分析模型、以及进行全面的敏感性分析，为同类研究提供了范本。

六、 研究亮点

1. 研究方法的先进性与严谨性：这是首批应用大规模、系统性的两样本孟德尔随机化分析来探讨圆锥角膜与多种系统性疾病因果关系的研究之一。研究设计巧妙，通过遗传工具变量有效控制了混杂偏倚，为推断因果关系提供了比传统观察性研究更强的证据。
2. 对争议问题的澄清：特别是关于糖尿病与圆锥角膜关系的长期争议，本研究通过使用代表长期血糖水平的遗传工具，得出了“无保护作用”的明确结论，对既往相互矛盾的研究发现进行了重要整合和澄清。
3. 发现新的潜在因果风险因素：不仅确认了过敏性疾病的风险，还首次通过MR方法提示溃疡性结肠炎与圆锥角膜可能存在因果关联，为探索炎症性肠病与眼部疾病的共病机制开辟了新方向。
4. 数据规模与质量：研究利用了近年来发布的多个大型GWAS的汇总数据，样本量庞大，统计效能较高（多数分析统计效能>80%）。并且严格筛选工具变量，保证了分析的可靠性。
5. 全面的敏感性分析框架：研究并未止步于主要分析，而是进行了包括MR-Egger、加权中位数、MR-PRESSO、留一法以及CAUSE模型在内的一系列敏感性分析，极大地增强了研究结果的稳健性和可信度。

七、 其他有价值的内容

研究者在讨论部分也坦诚指出了本研究的局限性： 1. 残留多效性的可能性：尽管进行了多种检验，但孟德尔随机化分析无法完全排除未被检测到的水平多效性。 2. 疾病覆盖范围有限：由于缺乏公开可用的高质量GWAS数据，本研究未能纳入一些此前报道可能与圆锥角膜相关的疾病，如唐氏综合征、某些结缔组织病（如马凡综合征、埃勒斯-当洛综合征）。 3. 人群局限性：所有数据均来源于欧洲血统人群，因此研究结论可能无法直接推广至其他种族或族群。 4. 统计效能：对于部分暴露（如主动脉瘤、TSH），本研究的统计效能相对较低（45%和38%），未来需要更大样本量的研究来进一步确认或排除其微弱的因果效应。 这些局限性的阐明为未来研究指明了方向，例如需要收集更多种族的数据、开展针对特定疾病的更深入机制研究（如探究溃疡性结肠炎相关的炎症因子与角膜细胞生物学行为的关系）等。