《Frontiers in Immunology》：A regulatory network underlying idiopathic pulmonary fibrosis
编辑推荐：
摘要背景：特发性肺纤维化（IPF）是一种进行性间质性肺病，其中遗传易感性与上皮、免疫和间充质重塑相互作用。尽管染色体11p15.5位点包含位于MUC5B和TOLLIP附近的已确认IPF易感信号，但该区域更广泛的调控结构仍未完全解析。方法：研究人员整合了IPF全
摘要背景：特发性肺纤维化（IPF）是一种进行性间质性肺病，其中遗传易感性与上皮、免疫和间充质重塑相互作用。尽管染色体11p15.5位点包含位于MUC5B和TOLLIP附近的已确认IPF易感信号，但该区域更广泛的调控结构仍未完全解析。方法：研究人员整合了IPF全基因组关联研究汇总统计数据与甲基化、表达和蛋白质数量性状位点，使用基于汇总数据的孟德尔随机化（SMR）。SMR优先排序的候选基因在独立的转录组学和甲基化队列中进行了评估，并通过microRNA、转录因子、蛋白质相互作用、机器学习、单细胞和空间转录组学分析进一步情境化。在博来霉素诱导的肺纤维化大鼠模型中评估了纤维化相关的表达模式。结果：分析恢复了已确认的MUC5B和TOLLIP信号，并将BRSK2优先排序为一个相对未被充分探索的候选基因，该基因得到基于eQTL的SMR和独立分子证据的支持。BRSK2 pQTL关联未通过HEIDI检验，因此未被解释为趋同的蛋白质水平遗传证据。网络分析将BRSK2与细胞周期、代谢应激和衰老相关程序联系起来，而跨队列机器学习将FOXA2、CDC25B和NFE2优先排序为信息性网络特征。单细胞和空间分析将BRSK2优先定位于成纤维细胞和肌成纤维细胞区室，以及组织学纤维化严重程度较高的区域。在纤维化大鼠肺中，BRSK2表达增加，而FOXA2和CDC25B在转录和蛋白质水平上均降低。结论：这些发现细化了染色体11p15.5 IPF易感位点的分子图谱，并将BRSK2优先排序为IPF相关促纤维化成纤维细胞状态的候选成分。其因果贡献、直接调控关系和治疗可行性需要有针对性的机制验证。
**论文解读：特发性肺纤维化染色体11p15.5易感位点的调控网络解析**
**研究背景**
特发性肺纤维化（IPF）是一种进行性且最终致命的间质性肺病，其特征是持续性成纤维细胞活化、细胞外基质过度沉积和肺泡结构不可逆破坏。当前IPF发病机制模型强调受损上皮细胞、免疫细胞群、成纤维细胞和肌成纤维细胞在空间异质性肺微环境中的适应不良相互作用。然而，遗传易感性如何转化为细胞类型特异性转录程序和持续性间充质重塑的分子机制仍不完全清楚。遗传学研究已显著推进了对IPF易感性的理解，并反复涉及染色体11p15.5区域，包括粘蛋白5B（MUC5B）和Toll相互作用蛋白（TOLLIP）附近的变异。然而，将疾病相关变异分配给其功能性效应基因仍具挑战性，因为大多数风险变异位于非编码区，可能通过组织和细胞状态特异性调控机制起作用，且不一定调控最近的基因。因此，11p15.5位点更广泛的调控结构，特别是超出已确认的MUC5B和TOLLIP信号的部分，仍未被充分解析。整合全基因组关联研究数据与分子数量性状位点，包括甲基化、表达和蛋白质QTL，为优先排序其分子性状与疾病风险共享遗传关联的基因提供了途径。在此背景下，基于汇总数据的孟德尔随机化（SMR）有其用途，但其解释仍对连锁不平衡、异质性、组织背景和工具变量有效性敏感，因此需要来自互补分子和实验证据的佐证。本研究旨在通过多层整合框架系统解析染色体11p15.5 IPF易感位点的调控结构，区分已确定的易感信号与额外候选调控因子，定义其在纤维化肺中的细胞和空间背景，并评估所得关联在不同分子层面的一致性。
**研究方法**
研究人员开展了一项多层面整合研究，主要关键技术方法包括：1）基于汇总数据的孟德尔随机化（SMR）分析，整合IPF GWAS汇总统计数据与来自DNA甲基化（mQTL）、基因表达（eQTL）和蛋白质丰度（pQTL）的数量性状位点数据集，并使用异质性检验（HEIDI）区分共享遗传信号与潜在连锁不平衡；2）在独立转录组学和甲基化队列（GSE213001, GSE134692, GSE231693, GSE175458）中对SMR优先排序的候选基因进行差异表达和甲基化验证；3）构建整合microRNA（miRNA）、转录因子（TF）和蛋白质-蛋白质相互作用（PPI）的调控网络，并进行功能富集分析；4）应用包含多种算法的机器学习框架（如弹性网络、岭回归、套索等）对BRSK2相关网络进行特征优先排序，并使用SHAP分析量化特征贡献；5）整合单细胞RNA测序数据（Seurat v5.1.0, Harmony v1.2.3）和空间转录组学数据（从E-MTAB-14121获取）分析候选基因的细胞和空间分布。样本队列来源包括公开数据库GWAS Catalog、Gene Expression Omnibus (GEO)和ArrayExpress，以及一个由4名对照和4名IPF患者组成的空间转录组学队列（共25个Visium切片）。此外，研究人员在博来霉素诱导的肺纤维化大鼠模型中进行实验验证，通过qRT-PCR和免疫荧光检测候选基因的表达变化。
**研究结果**
**3.1 多组学优先排序与IPF易感性相关的候选基因**
通过整合IPF GWAS与mQTL、eQTL和pQTL数据集的SMR分析，在染色体11p15.5易感位点识别出多个与IPF风险相关的分子性状。eQTL整合分析将MUC5B、TOLLIP和BRSK2优先排序。其中，MUC5B和TOLLIP是已确认的IPF易感基因，而BRSK2作为一个相对未被探索的候选基因，获得了显著的表达水平SMR证据支持。然而，BRSK2的pQTL关联未通过HEIDI检验（P
HEIDI = 1.76e-06），因此未被解释为独立的蛋白质水平遗传证据。
**3.2 候选基因表达和DNA甲基化的独立验证**
在独立IPF转录组数据集中，MUC5B和BRSK2的表达显著增加，TOLLIP表达显著降低，与eQTL-based SMR估计方向一致。差异甲基化分析在BRSK2位点内识别出7个CpG位点（如cg07652408, cg14306344等）的显著改变，以及MUC5B和TOLLIP相关CpG位点的改变。
**3.3 候选IPF基因的整合调控网络**
研究人员构建了整合miRNA、TF和PPI的调控网络。miRNA分析显示BRSK2和MUC5B共享多个候选miRNA（如hsa-miR-214-3p）。TF网络分析发现EGR1、MYC、TFAP2A在IPF中增加，而FOXA2、NFE2减少。PPI网络分析将BRSK2与CDC25B、WEE1、PKD1等细胞周期和应激反应相关蛋白连接。KEGG富集分析显示，BRSK2相关网络富集于细胞周期调节、细胞衰老和AMPK信号通路。
**3.4 机器学习识别与IPF相关的关键调控特征**
多算法机器学习框架应用在三个独立转录组队列中，FOXA2、CDC25B和NFE2被识别为最具影响力的三个特征，其低表达将模型预测向IPF类别偏移。单基因ROC分析显示CDC25B、FOXA2和TOLLIP具有高区分能力。
**3.5 调控网络在IPF中的细胞和空间背景**
整合的单细胞RNA测序分析显示，BRSK2优先在成纤维细胞和肌成纤维细胞中表达，而MUC5B主要在上皮细胞（特别是AT2细胞）中表达，TOLLIP分布更广泛。空间转录组学分析显示，BRSK2在轻度和中度纤维化区域表达相对有限，但在中度和重度纤维化区域更为突出，与其在成纤维细胞和肌成纤维细胞中的优先表达一致。
**3.6 对BRSK2相关网络成药性的计算评估**
通过数据库挖掘和文献整理，研究人员识别出多个候选化合物。分子对接显示GW296115和星形孢菌素对BRSK2具有最高的预测结合评分。随后进行的100纳秒分子动力学模拟评估了BRSK2-喜树碱、BRSK2-木犀草素等复合物的稳定性。
**3.7 肺纤维化中BRSK2相关表达模式的实验验证**
在博来霉素诱导的肺纤维化大鼠模型中，通过qRT-PCR证实BRSK2、MUC5B、MYC、PKD1和WEE1的mRNA水平显著增加，而CDC25B、FOXA2和TOLLIP降低。免疫荧光染色在蛋白质水平上验证了BRSK2增加，FOXA2和CDC25B减少。
**结论与讨论**
研究结论总结为：这些发现细化了染色体11p15.5 IPF易感位点的分子图谱，并将BRSK2优先排序为IPF相关促纤维化成纤维细胞状态的候选成分。其因果贡献、直接调控关系和治疗可行性需要有针对性的机制验证。
讨论部分指出，该研究将MUC5B、TOLLIP和BRSK2分别置于上皮、免疫相关和间充质重塑背景下。BRSK2的优先排序基于其通过HEIDI检验的eQTL-based SMR关联、一致的表达和甲基化改变，以及其在成纤维细胞/肌成纤维细胞中的富集和与纤维化严重程度的关联。然而，其pQTL关联的异质性、缺乏精细定位和染色质相互作用数据、部分网络分析基于数据库预测以及博来霉素模型与人类IPF的差异均是本研究的局限性。研究人员强调，BRSK2的因果作用、与TGF-β1响应性染色质调控的关系及药理靶向性仍需针对性机制验证。
