关于克罗恩病因果关系洞察：氧化应激与遗传变异的作用

《Journal of Cellular and Molecular Medicine》：Causal Insights Into Crohn's Disease: The Role of Oxidative Stress and Genetic Variants

　　

克罗恩病（CD）是一种胃肠道慢性炎症性疾病，其中氧化应激是其发病机制中的一个重要促发因素。本研究利用了多组学数据，包括来自GSE216447数据集的RNA测序和三个全基因组关联研究（GWAS）数据集，以研究CD中与氧化应激相关的分子网络。使用DESeq2进行差异表达分析，随后使用clusterProfiler进行通路富集分析。使用STRING数据库构建蛋白质-蛋白质相互作用（PPI）网络。使用TwoSampleMR和MRMix进行孟德尔随机化（MR）分析，以确定遗传变异与CD之间的因果关系。进一步应用定量实时聚合酶链反应在来自静止期和活动期CD患者的肠黏膜样本中验证关键差异表达基因，包括FASN、HMGCR、ASCC3、CD101、ELOVL6、PHLDA2、PHLDA3和SCPEP1。分析在对H₂O₂干预的反应中鉴定出64个上调和46个下调的差异表达基因。与氧化应激相关的关键通路，包括p53信号通路和类固醇生物合成，被显著富集。与转录组学数据一致，qPCR证实FASN和HMGCR在静止期CD中显著上调，而ASCC3、CD101、ELOVL6、PHLDA2、PHLDA3和SCPEP1在活动期CD中显著增加（均p < 0.05）。MR分析显示，在数据集中，使用TwoSampleMR鉴定出ABCB9和OSGIN1与CD具有显著因果关系，而在MRMix中仅OSGIN1显著。在数据集中，发现ARL4C、CD101、HMGCR和IL24与CD显著相关，且TwoSampleMR和MRMix的结果存在重叠。对于数据集，ACTA2和CD101被一致鉴定为显著，表明它们在CD发病机制中的潜在作用。研究结果突出了氧化应激相关分子网络在CD中的关键参与，并强调了孟德尔随机化在阐明因果遗传因素中的实用性。qPCR验证证实了静止期CD中脂质代谢基因的持续上调以及活动期CD中炎症相关基因的显著升高，强化了氧化应激与疾病活动度之间的联系。

克罗恩病是一种以复发和缓解交替为特征的慢性炎症性肠病，其发病机制复杂，涉及遗传易感性、免疫反应、环境因素和氧化应激的多方面相互作用。氧化应激被定义为活性氧（ROS）产生与机体中和能力之间的失衡，其通过破坏肠道屏障功能、促进炎症反应和改变肠道菌群，在CD的发生发展中发挥核心作用。临床缓解期CD患者体内仍存在持续的氧化应激标志物升高现象，表明分子水平的炎症持续存在。鉴于氧化应激的关键作用，理解其内在分子网络对于识别潜在治疗靶点至关重要。多组学方法可全面揭示驱动疾病进程的分子网络，而结合孟德尔随机化（MR）分析能够推断遗传变异与疾病结局之间的因果关系。本研究旨在通过多组学方法探索CD中氧化应激潜在的分子网络，鉴定参与氧化应激反应的关键基因和通路。研究发现已发表于《Journal of Cellular and Molecular Medicine》。研究对于阐明分子机制及开发调节氧化应激的靶向疗法具有重要意义。

在研究方法上，研究人员主要整合了GSE216447数据集中H₂O₂干预的RNA测序数据及三个不同人群的GWAS数据集。其中ieu-a-10数据来自欧洲人群，ieu-a-11和ieu-a-13数据分别来自东亚和南亚人群。研究利用DESeq2进行差异表达分析，并用clusterProfiler进行GO和KEGG通路富集分析。通过基因集变异分析（GSVA）评估通路活性变化，并使用STRING数据库构建蛋白质-蛋白质相互作用（PPI）网络。基于TwoSampleMR和MRMix方法进行MR分析以推断因果遗传关联，工具变量均选自GWAS数据。为验证测序结果，研究纳入了来自20名健康对照、22名静止期以及24名活动期CD患者的肠黏膜样本进行定量实时聚合酶链反应验证。

DEGs Analysis of H2O2 Intervention
通过对H₂O₂干预后的样本进行测序分析，研究人员共鉴定出64个上调和46个下调的DEGs。其中上调最显著的基因包括GDF15、CDKN1A和MDM2等，而下调最显著的基因包括DCT、HMGCS1和LDLR等。聚类分析能够有效区分干预组与对照组，表明H₂O₂干预引起了基因表达模式的显著改变。

GO and KEGG Pathway Analysis
GO分析鉴定出109个显著富集的生物过程通路，主要涉及固醇生物合成和代谢过程。KEGG分析显示，p53信号通路、类固醇生物合成以及FoxO信号通路被显著富集。这表明氧化应激不仅引起细胞周期和凋亡相关通路的改变，还深度干扰了固醇代谢与脂质平衡。

GSVA Analysis of H2O2 Intervention
GSVA分析揭示了11条显著变化的GO通路。其中DNA导向的DNA聚合酶活性调节、蛋白酶体调节颗粒和蛋白酶体核心复合物等通路得分呈正值，表明其受到上调。该结果突显了在氧化应激环境中维持基因组稳定和加速受损蛋白降解的生物学反应机制。

PPI Network Analysis
基于鉴定出的DEGs，研究人员构建了蛋白质-蛋白质相互作用网络，并将DEGs划分为四个相互连接的聚类。这一分析有助于深入理解在H₂O₂处理后细胞内不同基因协同作用的分子机制。

Differential Expression Analysis on of H2O2 Intervention DEGs in CD Patients
研究人员在静止期CD的测序数据中重点考察了受H₂O₂干预影响的前100个DEGs，并确定FASN和HMGCR在静止期CD患者中显著上调。通过qPCR对样本队列进行验证，证实FASN和HMGCR的mRNA水平确实显著升高。对于活动期CD样本，研究发现ASCC3、CD101、ELOVL6、PHLDA2、PHLDA3和SCPEP1存在显著差异表达，qPCR验证也进一步确认了这些基因在活动期CD患者肠黏膜样本中的显著上调。

Causal Relationships Identified in CD Using TwoSampleMR and MRMix
通过MR分析识别因果关联。利用TwoSampleMR在ieu-a-10中识别出ABCB9和OSGIN1与CD存在显著因果关系，而MRMix仅识别出OSGIN1。在ieu-a-11中，ARL4C、CD101、HMGCR和IL24被TwoSampleMR确认为显著，其中ARL4C、CD101被两种方法共同验证。在ieu-a-13中，ACTA2和CD101均被两种方法鉴定为具有显著的因果效应，凸显了这些基因在CD发病机制中的潜在核心作用。

研究讨论部分指出，氧化应激在CD的发病机制中扮演了不容忽视的角色。GDF15和CDKN1A等上调基因在应对DNA损伤与抑制细胞周期中发挥保护作用，而MDM2的升高可能通过负调节p53影响肠道上皮细胞凋亡。相反，DCT和HMGCS1等下调基因反映了黑色素生成与胆固醇生物合成途径的抑制。GO与KEGG分析揭示固醇代谢异常与激酶信号传导在应对氧化应激中的关键地位。基于GSVA结果，蛋白酶体降解途径活性的增加可能是对氧化损伤的补偿性 survival 机制。针对CD患者的基因表达研究证实，FASN和HMGCR等脂质代谢基因在静止期CD中的持续上调与代谢失调密切相关；而CD101等炎症相关基因在活动期CD中的显著增加可能驱动了致病性T细胞的扩增与免疫失衡。MR分析进一步从遗传变异角度强化了ARL4C、CD101、HMGCR和IL24等基因与CD的因果关系。研究结论部分表明，本研究提供了CD中氧化应激相关分子网络的重要见解。通过结合多组学与MR分析，强调了FASN、CD101和HMGCR等特定DEGs在CD发病机制中的关键作用。qPCR验证确认FASN和HMGCR在静止期CD中持续上调，而ASCC3、CD101、ELOVL6、PHLDA2、PHLDA3和SCPEP1在活动期CD中显著升高，从而强化了氧化应激与疾病活动度之间的联系，为靶向治疗的关键途径提供了分子靶点依据。