《Ophthalmology Science》：A human pharmacogenomics approach provides insight into the pathogenesis and pathophysiology of steroid-induced ocular hypertension
编辑推荐：
目的：识别与糖皮质激素诱导的眼高压（SI-OHT）相关的单核苷酸多态性（SNP）及其潜在靶基因，并进行功能注释分析。设计：全基因组关联研究（GWAS）。参与者：在单一临床中心接受角膜移植术后的福克斯角膜内皮营养不良（FECD）患者（N=439）。术后患者自我滴
目的：识别与糖皮质激素诱导的眼高压（SI-OHT）相关的单核苷酸多态性（SNP）及其潜在靶基因，并进行功能注释分析。设计：全基因组关联研究（GWAS）。参与者：在单一临床中心接受角膜移植术后的福克斯角膜内皮营养不良（FECD）患者（N=439）。术后患者自我滴注1%醋酸泼尼松龙滴眼液。分别在术前及术后1、3、6和12个月测量眼压（IOP）。方法：收集唾液样本提取DNA并进行基因分型。以IOP的最大变化值作为GWAS的定量性状。采用HRC参考面板进行基因型填补后，利用线性回归分析确定定量性状与基因型之间的关联。复制分析纳入4个独立队列（N=49-103），这些参与者此前已进行过基因分型。主要结果指标：研究人员发现共46个达到全基因组显著性水平（P<5×10-8）的SNP，聚类于29个不同风险位点；另有623个达到提示性显著性水平（P<5×10-6）的SNP，聚类于323个风险位点。大多数SNP为罕见或低频变异，且具有较大的效应值。研究人员通过比对基因表达谱数据及注释分析，验证了441个优先靶基因。在前29个风险位点中，31%与其他高眼压表型的风险位点共定位，另有2个位点通过生物学证据与其他高眼压表型相关联，总重叠率达38%。三个提示性显著性的SNP在不同队列中得到独立验证。许多优先靶基因并不在房水流出通路的小梁网（TM）或邻管组织（JCT）中表达，这与当前的疾病范式相悖。注释分析提示了新的病理生理机制。结论：自20世纪50年代GC首次用于治疗眼病以来，为何部分患者会发生SI-OHT而另一些则不会一直未明。本研究结果表明SI-OHT具有遗传学基础，支持其与其他高眼压表型的关联性，并为后续实验室研究提供了假设依据。已发现并验证的SNP可用于SI-OHT风险预测，优先靶基因则有望成为SI-OHT管理的干预靶点。
研究背景与意义
糖皮质激素（GC）广泛应用于眼科抗炎治疗及术后管理，然而约30%至50%的成年患者会出现眼压升高，即糖皮质激素诱导的眼高压（SI-OHT），严重者可进展为青光眼。尽管这一现象自20世纪50年代被发现以来长期困扰临床，但其背后的遗传机制始终未明。既往研究多聚焦于小梁网细胞的体外基因表达改变，但受限于细胞模型的单一性及样本量不足，难以全面揭示复杂的体内调控网络。此外，原发性开角型青光眼（POAG）患者多为类固醇反应者，且类固醇反应者发生POAG的风险显著高于非反应者，提示SI-OHT与其他高眼压表型之间存在共同的遗传基础。因此，阐明SI-OHT的遗传结构，不仅有助于实现个性化的GC用药风险评估，也为开发新的治疗靶点提供了理论依据。本研究成果发表于《Ophthalmology Science》。
关键技术方法
研究人员采用全基因组关联研究（GWAS）设计，构建了包含439例欧洲裔福克斯角膜内皮营养不良（FECD）患者的发现队列，所有患者均接受了穿透性角膜内皮移植术并规范使用1%醋酸泼尼松龙滴眼液。研究以术后IOP相对于基线的最大变化值（deltaIOPmax）作为定量性状。基因分型采用Illumina Infinium Global Screening Array-24芯片，并利用HRC参考面板进行基因型填补。统计分析采用线性混合效应模型，校正了基线IOP、年龄、性别及前四个主成分，并引入遗传关系矩阵以控制隐性亲缘关系。为验证结果，研究纳入了4个独立的复制队列，分别来自美国印第安纳波利斯、佛罗里达及印度金奈，涵盖不同的GC给药方式及种族背景。研究人员进一步结合GTEx数据库、RegulomeDB及DAVID等生物信息学工具，对风险位点进行了功能注释、靶基因优先级排序及通路富集分析。
研究结果
SNP关联分析结果
研究人员分别对术后3个月和12个月的IOP变化进行了GWAS分析。3个月的分析因避开了药物减量及临床干预，显示出更强的统计效力。分析共鉴定出46个达到全基因组显著性（P<5×10-8）的SNP，分布于29个独立风险位点；此外还有623个达到提示性显著性（P<5×10-6）的SNP。与常见复杂疾病不同，这些SI-OHT相关的SNP大多为罕见或低频变异（MAF<0.05），且表现出极大的效应值，最高可达27.33 mmHg。
遗传负荷分析
针对全基因组显著位点构建的遗传风险评分显示，个体携带的风险等位基因数量与IOP升高幅度呈显著正相关。3个月结局的遗传负荷每增加一个单位，deltaIOPmax平均增加2.58 mmHg（P=4.33×10-37）；12个月结局同样呈现强相关性，证实了累积遗传效应对疾病表型的显著影响。
靶基因优先级排序与验证
研究人员假设靶基因应在房水流出通路（AOPs）中表达，并结合GTEx数据库的eQTL/sQTL数据，最终筛选出441个优先靶基因。通过与既往GC处理的转录组数据进行比对，71%的靶基因被证实受GC调控。值得注意的是，许多关键基因（如COL11A1、RBFOX1等）并不在小梁网或邻管组织中高表达，提示SI-OHT的病理机制可能涉及眼内更广泛的细胞类型及系统性调控网络。
风险位点与其他高眼压表型的共定位
在29个顶级风险位点中，31%（9个位点）与POAG、原发性闭角型青光眼（PACG）、剥脱综合征（PEXS）等其他高眼压表型的风险位点存在基因组共定位。其中，位于COL11A1基因区域的风险位点最为显著，与5种不同表型（包括II型Stickler综合征）相关联，但这些关联SNP并非相同的变异，也不处于连锁不平衡状态，提示不同疾病可能通过影响同一基因的不同调控元件而致病。
独立复制
虽然全基因组显著的稀有变异因频率过低未能在所有队列中复制，但在合并分析的3个月QT数据中，有5个提示性SNP的P值得到改善。更重要的是，在印度金奈的队列中，研究人员成功复制了3个SNP，分别位于AGTR1和RP1基因区域，且效应方向一致。
功能注释
针对COL11A1位点的深入分析发现，SI-OHT相关的稀有变异多位于基因下游或内含子区域，且多为剪接数量性状位点（sQTL），倾向于抑制第41号内含子的正确剪切，导致功能性蛋白表达下降。这一机制在不同高眼压表型中呈现出趋同性。通路富集分析揭示了三个核心生物学过程：心肌传导相关的钾离子通道功能（Muscle Contraction）、Notch信号通路（Signaling by NOTCH）以及糖胺聚糖代谢（Metabolism of Carbohydrates）。其中，Notch信号被证实与GC受体（GR）信号存在交叉对话，其抑制是GC发挥转录激活作用的前提。
讨论与结论
本研究通过大规模药物基因组学分析，首次系统描绘了SI-OHT的遗传图谱。研究发现SI-OHT主要由大量低频和罕见的功能缺失型变异驱动，这解释了为何传统针对常见变异的GWAS在此领域进展缓慢。研究挑战了“SI-OHT仅源于小梁网纤维化”的传统观念，指出COL11A1等基因的表达下调导致的细胞外基质生物力学改变，以及Notch等非经典信号通路的紊乱，共同构成了SI-OHT的复杂病理基础。此外，SI-OHT与POAG虽共享部分风险基因，但致病变异各异，支持两者为发病机制不同而终末表型相似的独立疾病。综上所述，该研究所鉴定的SNP组合可作为临床预测GC致眼压升高的有效工具，而AGTR1、NOTCH及COL11A1等相关通路则为开发防止视力丧失的新型疗法指明了方向。
