《International Journal of Chronic Obstructive Pulmonary Disease》：Identification and Validation of Mannose Metabolism-Related Biomarkers in COPD Through Integrated Bioinformatics and Machine Learning Analysis: A Pilot Study
编辑推荐：
背景：慢性阻塞性肺疾病（COPD）是一种进行性呼吸系统疾病，以持续气流受限和持续炎症过程为特征。近年来，甘露糖代谢已成为慢性疾病发展的重要因素。本研究探讨了甘露糖代谢相关基因（MMRGs）如何参与COPD发病机制，并评估了其作为候选诊断和治疗靶点的效用。方法：
背景：慢性阻塞性肺疾病（COPD）是一种进行性呼吸系统疾病，以持续气流受限和持续炎症过程为特征。近年来，甘露糖代谢已成为慢性疾病发展的重要因素。本研究探讨了甘露糖代谢相关基因（MMRGs）如何参与COPD发病机制，并评估了其作为候选诊断和治疗靶点的效用。方法：研究人员通过GEO数据库获取COPD患者和健康对照者的血样基因表达数据。鉴定差异表达基因（DEGs）并与MMRGs取交集以获得候选基因。应用三种机器学习算法并结合独立数据集表达验证来识别生物标志物。构建列线图预测模型，并使用受试者工作特征（ROC）曲线分析评估其诊断性能。随后进行基因集富集分析（GSEA）、免疫浸润分析、药物预测和分子对接。结果：从1685个DEGs中鉴定出19个候选基因，随后筛选出两个生物标志物：MAN1C1和MAN2B2。基于两者构建的列线图模型显示出中等区分效能（曲线下面积（AUC）=0.701）。此外，GSEA分析显示两者均共同富集于TNF靶点上调基因集等通路。免疫浸润结果显示，COPD与对照组之间共有12类免疫细胞存在显著差异（p < 0.05），例如活化B细胞。最后，药物预测分别揭示MAN1C1和MAN2B2有12种和3种潜在药物，其中曲古抑菌素A显示出潜在的结合构象。结论：本研究揭示了MMRGs在COPD中的潜在作用，并鉴定了新型生物标志物。这些发现为COPD的早期诊断、个体化治疗和药物开发提供了新的见解和研究基础。
**论文解读：甘露糖代谢相关生物标志物在慢性阻塞性肺疾病中的鉴定与验证**
**研究背景与问题**
慢性阻塞性肺疾病（chronic obstructive pulmonary disease, COPD）是一种以持续性气流受限和慢性气道炎症为特征的进行性呼吸系统疾病，全球疾病负担沉重，2021年全球疾病负担（GBD）数据显示其导致约350万人死亡，预计2020至2050年间全球经济损失达5.5万亿美元。尽管早期诊断和治疗取得进展，但COPD的病理生理机制仍复杂不清，亟需发现潜在生物标志物以促进早期诊断、靶向治疗开发及改善患者预后。近年来，甘露糖代谢（mannose metabolism）在慢性炎症和肿瘤等疾病中的关键作用被逐渐揭示，但甘露糖代谢相关基因（mannose metabolism-related genes, MMRGs）在COPD中的角色尚未明确。外周血作为全身免疫和炎症反应的关键载体，可反映肺部病理变化，且易于获取，适宜大规模探索性分析。因此，研究人员开展本研究，旨在通过整合生物信息学与机器学习方法，鉴定MMRGs中与COPD相关的生物标志物，并评估其诊断和治疗潜力。该研究发表于《International Journal of Chronic Obstructive Pulmonary Disease》。
**主要技术方法**
研究人员利用公共基因表达数据库（Gene Expression Omnibus, GEO）获取训练集GSE42057（94例COPD患者、42例健康对照）和验证集GSE94916（6例COPD患者、6例健康对照）的外周血转录组数据。主要关键技术包括：①差异表达分析（运用limma R包筛选差异表达基因，differentially expressed genes, DEGs）；②与从分子特征数据库（Molecular Signatures Database, MSigDB）提取的160个MMRGs取交集，获得候选基因；③三种机器学习算法——最小绝对收缩和选择算子（least absolute shrinkage and selection operator, LASSO）回归、随机森林（random forest, RF）和支持向量机-递归特征消除（support vector machine-recursive feature elimination, SVM-RFE）联合筛选候选生物标志物，并在独立数据集进行表达验证；④基于筛选出的生物标志物构建列线图（nomogram）预测模型，通过受试者工作特征（receiver operating characteristic, ROC）曲线和校准曲线评估诊断效能；⑤基因集富集分析（gene set enrichment analysis, GSEA）和基因集变异分析（gene set variation analysis, GSVA）探索通路富集情况；⑥单样本基因集富集分析（ssGSEA）进行免疫浸润分析；⑦通过GeneMANIA数据库构建功能交互网络，利用ENCORI和miRDB数据库构建lncRNA-miRNA-mRNA调控网络；⑧基于DSigDB数据库进行药物预测，并采用CB-Dock2进行分子对接分析。
**研究结果**
*共获得19个候选基因*：在GSE42057训练集中，共鉴定出1685个DEGs（1208个上调，477个下调）。与160个MMRGs取交集后获得19个重叠基因，定义为候选基因。基因本体论（Gene Ontology, GO）富集分析显示这些候选基因主要富集于糖脂生物合成过程、GPI锚定生物合成过程、α-甘露糖苷酶活性等通路；京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）富集分析揭示其显著参与GPI锚定生物合成、果糖和甘露糖代谢、AMPK信号通路等。蛋白-蛋白相互作用（protein-protein interaction, PPI）网络显示候选基因编码蛋白间存在复杂交互，如PGAP2与PIGO。
*MAN1C1和MAN2B2被鉴定为生物标志物*：通过LASSO、RF和SVM-RFE三种机器学习算法筛选，取交集获得5个候选生物标志物（MAN1C1、MTMR1、DPY19L1、PIGK、MAN2B2）。在GSE42057和GSE94916两个数据集中进行表达验证，仅MAN1C1和MAN2B2在两个数据集中均表现出统计学显著差异（p < 0.05）且变化方向一致，即在COPD标本中显著下调。亚细胞定位分析显示MAN1C1主要定位于高尔基体、细胞外基质和内质网，MAN2B2主要定位于溶酶体和细胞外区域。
*基于MAN1C1和MAN2B2构建的列线图模型在COPD风险预测中显示出中等诊断效能*：基于两个生物标志物构建的列线图模型，总分为85.8时对应的疾病预测概率为72.4%。校准曲线显示预测概率与实际发生概率高度一致（Hosmer-Lemeshow拟合优度检验p=0.36），ROC曲线分析显示曲线下面积（area under the curve, AUC）为0.701，表明模型具有中等区分能力。
*GSEA和GSVA分析揭示MAN1C1和MAN2B2在COPD免疫应答和代谢异常中的潜在作用*：GSEA分析显示MAN1C1和MAN2B2共同富集于38条通路，包括TNF靶点上调基因集（TNF target up）和蕈样肉芽肿CD4
+T细胞特征上调基因集等。MAN1C1还富集于血管免疫母细胞T细胞淋巴瘤下调基因集等，MAN2B2富集于肿瘤坏死因子信号通路下游靶基因、血小板特异性基因表达特征等。GSVA分析显示，与对照组相比，COPD标本中短链脂肪酸代谢和骨骼肌底物代谢相关通路显著上调，而CHSY1和CHST3功能缺陷相关通路下调。
*免疫细胞浸润分析揭示MAN1C1与COPD免疫微环境的关系*：采用ssGSEA分析28种免疫细胞，发现COPD组与对照组间12种免疫细胞浸润水平存在显著差异（p < 0.05），包括活化B细胞、活化CD4 T细胞、活化CD8 T细胞、活化树突状细胞（DC）、中央记忆CD4 T细胞、中央记忆CD8 T细胞、效应记忆CD8 T细胞、γδ T细胞、未成熟DC、巨噬细胞、中性粒细胞和调节性T细胞（Treg）。其中，活化B细胞、活化CD4 T细胞、活化CD8 T细胞和中央记忆CD4 T细胞在COPD组中显著下调；活化DC、中央记忆CD8 T细胞、γδ T细胞、未成熟DC、巨噬细胞、中性粒细胞和Treg显著上调。相关性分析显示巨噬细胞与中性粒细胞呈最强正相关（cor=0.84），活化CD8 T细胞与巨噬细胞呈最强负相关（cor=-0.66）。MAN1C1与活化B细胞（cor=0.34）、活化CD4 T细胞（cor=0.57）、中央记忆CD4 T细胞（cor=0.37）和效应记忆CD8 T细胞（cor=0.30）呈显著正相关。
*MAN1C1和MAN2B2的功能相互作用和转录后调控机制*：GeneMANIA网络分析显示，MAN1C1和MAN2B2的功能相互作用以物理交互为主（77.64%），并涉及共表达、计算预测、共定位等。富集分析显示两者与甘露糖苷酶活性、O-聚糖特异性糖苷水解酶功能、糖苷键裂解活性显著相关。构建的lncRNA-miRNA-mRNA调控网络揭示，MAN1C1受17个miRNA（如hsa-miR-376a-3p）和9个lncRNA调控，MAN2B2受5个miRNA（如hsa-miR-627-5p）调控，表明两者具有不同的转录后调控机制。
*药物预测和分子对接分析揭示潜在调控药物*：mRNA-药物相互作用网络显示，MAN1C1对应12种潜在靶向药物（如AM580、曲古抑菌素A），MAN2B2对应3种（如EINECS 250-892-2、维替泊芬）。分子对接分析显示，曲古抑菌素A与MAN1C1蛋白的结合能为-7.0 kcal/mol，与活性位点氨基酸残基GLN-535、GLU-604和SER-602形成氢键；EINECS 250-892-2与MAN2B2蛋白的结合能为-7.9 kcal/mol，同样与上述残基形成氢键，表明两者可能作为潜在调控药物。
**讨论与结论**
讨论部分指出，MAN1C1（α-1,2-甘露糖苷酶I）在COPD中下调，可能与活化B细胞和CD4 T细胞浸润正相关，通过影响糖基化调控受体聚集和信号转导，参与慢性炎症进程。MAN2B2（α-1,6-甘露糖苷酶）参与溶酶体糖蛋白降解，其功能缺陷与糖基化相关免疫失调相关。GSEA显示两者共同富集于TNF信号通路，与既往研究一致，提示其可能通过糖代谢和糖基化修饰影响细胞因子受体和免疫信号网络。GSVA发现短链脂肪酸代谢通路在COPD中上调，支持肠道菌群-肺轴在COPD免疫失调中的作用。免疫浸润分析揭示COPD患者存在免疫反应调节异常，巨噬细胞与中性粒细胞强正相关，MAN1C1与活化CD4 T细胞正相关，提示其可能在免疫微环境中发挥重要作用。药物预测和分子对接表明曲古抑菌素A（组蛋白去乙酰化酶抑制剂）和EINECS 250-892-2可能通过调节免疫细胞功能影响COPD进展。研究结论：本研究提示MMRGs可能参与COPD发病机制，并可作为潜在生物标志物，为COPD的早期诊断、个体化治疗和药物开发提供了新的见解和研究基础。但需注意，这些发现基于生物信息学分析，存在样本量小、缺乏实验验证等局限，需进一步大规模队列和实验研究确认。
