《Bioinformatics and Biology Insights》：Tear MicroRNA-Based Molecular Signature for Keratoconus Stratification: A Cross-Cohort Integration Study
编辑推荐：
摘要专业翻译
圆锥角膜（Keratoconus，KC）是一种进行性角膜扩张性疾病，其特征为结构重塑和免疫相关分子改变。然而，泪液来源的微小RNA（microRNA，miRNA）能否反映角膜分子变化仍不清楚。研究人员分析了三个公开可用的角膜转录组数据集（N=11
摘要专业翻译
圆锥角膜（Keratoconus，KC）是一种进行性角膜扩张性疾病，其特征为结构重塑和免疫相关分子改变。然而，泪液来源的微小RNA（microRNA，miRNA）能否反映角膜分子变化仍不清楚。研究人员分析了三个公开可用的角膜转录组数据集（N=112），通过差异表达分析和基因集富集分析（Gene Set Enrichment Analysis，GSEA）鉴定免疫相关基因表达模式。在一项前瞻性泪液队列（N=20）中，研究人员进行了泪液miRNA谱分析，并使用最小绝对收缩和选择算子（Least Absolute Shrinkage and Selection Operator，LASSO）回归筛选候选miRNA。研究人员将miRNA模式与免疫相关转录组改变进行整合，并使用标准化hsa-miR-326表达和角膜曲率测量值（K2）构建了探索性分子-结构指数。研究人员在KC角膜中鉴定了免疫相关转录改变，并推导出五-miRNA泪液特征。该特征在独立角膜上皮miRNA数据集中表现出判别性能（受试者工作特征曲线下面积AUC=0.89，95%置信区间CI：0.69-1.00；N=16），并在另一个泪液队列中通过逆转录定量聚合酶链反应（RT-qPCR）进一步评估（N=52）。在这些候选miRNA中，hsa-miR-326在各数据集中表现出一致的下调。网络分析确定了候选miRNA与免疫相关基因（包括IL1B、CD4、CXCL10和TGFB1）之间的关联。探索性分子-结构指数与角膜曲率陡峭化呈负相关（Spearman ρ = -0.86，p < 0.001）。这些发现提示，泪液来源的miRNA可能提供与角膜免疫相关改变互补的分子信息，并值得在KC特征化中进一步研究。
论文解读
一、研究背景与意义
圆锥角膜（Keratoconus，KC）是一种进行性角膜扩张性疾病，以角膜基质变薄、前凸、瘢痕形成和不同程度视觉畸变为特征。传统上KC被视为非炎症性疾病，但越来越多的证据表明其是一种受遗传易感性、环境刺激、氧化应激和慢性揉眼行为共同影响的多因素复杂疾病。基质金属蛋白酶（Matrix Metalloproteinases，MMPs）和促炎介质的表达改变已被报道。组织学上，KC涉及胶原排列紊乱、角膜细胞密度降低和蛋白聚糖组成改变，导致组织水平生物力学性能受损和细胞外基质（Extracellular Matrix，ECM）重塑。
当前基于影像的诊断方法所依赖的临床参数（如角膜曲率测量值）主要反映已确立的结构改变，在识别疾病早期分子水平变化方面能力有限。泪液作为非侵入性生物标志物的潜在来源受到关注，其中微小RNA（microRNA，miRNA）因其在体液中良好的稳定性而成为优秀的生物标志物候选。然而，跨队列变异因患者异质性和研究设计差异而仍然显著，且泪液来源miRNA与角膜转录组改变的整合有限，阻碍了候选分子信号的生物学情境化解读。更重要的是，泪液来源分子信号与角膜结构改变之间的生物学联系仍不明确，限制了其在分层中的转化解释。为此，本研究开展了跨队列整合分析，以定义免疫相关基因特征，并通过机器学习从前瞻性收集的泪液miRNA谱中识别候选分子特征，最终将泪液来源miRNA特征与角膜结构参数整合，探索KC中的分子-结构耦合。
二、研究设计与技术方法
研究人员从NCBI数据库获取了三个独立RNA测序数据集（GEO登录号GSE77938及BioProject登录号PRJNA312169、PRJNA636666、PRJNA1184491），共112个角膜组织样本。采用共识K均值聚类进行样本分层与纯化，通过主成分分析（Principal Component Analysis，PCA）评估批次效应并使用limma包的removeBatchEffect函数进行校正。差异表达基因（Differentially Expressed Genes，DEGs）分析使用DESeq2包，筛选标准为校正后p值<0.05且|log2FoldChange|>2。基因集富集分析（Gene Set Enrichment Analysis，GSEA）使用MSigDB数据库的c2.cp.kegg.v11.0.symbols基因集，FDR阈值<0.05。
前瞻性泪液队列来自林口长庚纪念医院（IRB No. 202401867B0），包括20名参与者（10例KC，10例对照）进行miRNA测序。使用miRNeasy Mini Kit提取总RNA，在Illumina NovaSeq平台进行75bp单端测序，miRNA注释基于miRBase v22版本。特征选择采用LASSO逻辑回归模型（glmnet包），通过十折交叉验证确定最优正则化参数。miRNA-mRNA调控网络通过整合泪液来源miRNA特征与角膜基质转录组改变构建，预测靶标从miRDB数据库获取并与GSEA鉴定的174个核心mRNA取交集。蛋白质-蛋白质相互作用（Protein-Protein Interaction，PPI）网络使用STRING数据库（版本12.0）构建，通过Cytoscape的cytoHubba插件基于Degree中心性鉴定枢纽基因。独立验证队列包含52名参与者（20例KC、20例对照、8例干眼症、4例白内障），通过RT-qPCR验证hsa-miR-326表达，以U6作为内参，采用2-ΔΔCt方法进行相对定量。
三、研究结果
聚类分析与代表性组别选择
对整合队列（N=112）应用共识K均值聚类后，通过肘部法则确定最优聚类数为4。其中两个聚类具有明确生物学相关性：一个以对照样本为主（94.7%；KC：1，CON：18），另一个富集KC样本（80.3%；KC：45，CON：11）。排除两个混合组成聚类后，保留63个样本（45例KC，18例对照）用于下游分析。批次效应校正后，保留样本在PCA空间中KC与对照的分离度改善，PC1和PC2累计方差解释率从42.1%提升至47.7%。
KC与对照组差异表达分析
差异表达分析鉴定出5598个显著差异表达基因。GSEA显示KC组富集类固醇激素生物合成、黏蛋白型O-聚糖生物合成和视黄醇代谢等通路；对照组相对富集趋化因子信号、TNF信号和IL-17信号等免疫相关通路，这一模式可能反映KC组织中免疫稳态和应激反应信号的失调。
候选miRNA特征鉴定
泪液miRNA测序分析中，火山图鉴定出大量差异表达miRNA，其中hsa-miR-1260a为最显著下调标志物，hsa-miR-222-3p和hsa-miR-221-5p在KC样本中显著上调。LASSO回归筛选出7个候选miRNA：hsa-miR-1260a、hsa-miR-320b、hsa-miR-326、hsa-miR-30b-3p、hsa-miR-1291、hsa-miR-1290在KC组一致显著下调，hsa-miR-222-3p在KC样本中显著上调。
miRNA-mRNA调控网络整合与功能特征化
将7个候选miRNA与174个角膜基质来源的KEGG通路引导差异表达基因整合，鉴定出33个重叠靶基因。PPI网络分析基于Degree中心性确定前10个枢纽基因，其中IL1B、CD4、CXCL10和TGFB1为核心枢纽基因。精化网络突出5个核心miRNA（hsa-miR-326、hsa-miR-1291、hsa-miR-1290、hsa-miR-222-3p和hsa-miR-30b-3p）直接调控枢纽基因。hsa-miR-326为网络中连接度最高的节点，可能调控IL1B、TGFB1和MMP1等多个炎症和信号基因。功能富集分析显示这些枢纽基因参与细胞因子介导的信号通路、应激反应、细胞对机械刺激的反应和细胞迁移等关键生物学过程。
角膜上皮组织外部验证
在独立角膜上皮miRNA数据集（PRJNA842431，N=16）中，hsa-miR-1291（p=0.009）、hsa-miR-326（p=0.014）和hsa-miR-30b-3p（p=0.014）在KC组显著下调，与泪液发现高度一致；hsa-miR-1290（p=0.529）和hsa-miR-222-3p（p=0.894）无显著差异表达。基于泪液谱训练的五-miRNA面板在该外部验证队列中取得AUC=0.89（95% CI：0.69-1.00）的判别性能。
hsa-miR-326的RT-qPCR实验验证
在独立临床队列（N=52）中，RT-qPCR证实KC患者泪液hsa-miR-326表达较健康对照显著下调（p<0.05），且KC组表达水平最低，白内障、干眼症和健康对照组之间无显著差异。研究人员构建了探索性分子-结构指数（Molecular-Structural Index，MSI），定义为z(miR-326) ? z(K2)，该指数与Max K2呈显著负相关（Spearman ρ = -0.86，p<0.001，n=20）。
hsa-miR-326与临床特征关联
未观察到特应性、睑板腺功能障碍、瘙痒和隐形眼镜使用等眼表相关特征与hsa-miR-326表达存在显著关联。hsa-miR-326单独与Max K2、最小角膜厚度和ABCD最大评分无显著相关性，但MSI与Max K2强相关。性别分层分析未发现女性和男性参与者在hsa-miR-326表达、MSI、Max K2或最小角膜厚度方面存在显著差异。
四、讨论与结论
讨论部分指出，KC相关分子改变组织为相互依存的免疫-应激反应模块而非独立基因变化。外部验证的部分可重复性提示组织或区室依赖性变异，反映细胞内上皮调控与细胞外泪液分泌动力学之间的差异。研究结果支持KC与慢性免疫失调相关的观点，免疫激活在KC中组织为协调的信号网络，可能通过持续细胞因子信号破坏ECM合成与降解平衡，促进基质变薄和进行性生物力学减弱。单细胞转录组研究进一步支持KC中免疫相关失调存在于特定角膜细胞群体。研究强调，泪液来源miRNA可能反映与角膜结构改变相关的上游分子状态，但网络构建以泪液来源miRNA锚定至预设免疫相关角膜基因为基础，结果应视为假设生成的关联而非确定性调控机制证据。从转化角度看，泪液miRNA谱可能作为传统结构评估的补充分子信息，而非独立诊断生物标志物。MSI作为概念验证性复合指数，其与疾病严重程度的关联提示分子-结构联合表征可能比单一组分更好地捕捉疾病复杂性，但应严格视为假设生成的临床相关性指标而非经过验证的预测工具。
结论部分总结：本研究通过整合泪液来源miRNA谱与角膜转录组数据集，鉴定了KC中免疫相关分子特征。通过共识聚类降低生物学异质性、差异基因表达分析、通路富集和LASSO特征选择，研究人员鉴定了五个泪液来源miRNA和174个免疫相关mRNA，这些分子与区分KC和对照聚类的分子模式相关。IL1B、CXCL10、CXCL11、TGFB1和CD4等关键基因构成紧密互连的网络，参与细胞因子信号传导、机械应激反应和ECM重塑。多个泪液来源miRNA（包括hsa-miR-326、hsa-miR-320b和hsa-miR-30b-3p）在独立角膜上皮数据集中获得外部验证支持，提示跨数据集一致的生物学关联。其中hsa-miR-326作为整合网络中的优先候选分子，值得在未来的功能研究中进一步探索。该整合分析框架为连接泪液来源分子信号与KC角膜转录组改变提供了基础，并可能为未来研究疾病异质性的分子特征化奠定基础。
