COPD患者肺组织microRNA表达与定量肺气肿及肺功能关联的探索性转录组分析

《International Journal of Chronic Obstructive Pulmonary Disease》：Association of Lung Tissue MicroRNA Expression with Quantitative Emphysema in COPD: An Exploratory Transcriptomic Analysis

　　

背景：microRNA（miRNA）可能参与与肺气肿相关的生物学通路，但其与人肺组织中定量肺气肿的关联仍不确定。方法：研究人员分析了GSE201121中117例慢性阻塞性肺疾病（COPD）/肺气肿患者的肺组织miRNA表达和临床表型。经表达过滤后，研究人员采用线性回归评估839个miRNA位点与肺气肿百分比（定量计算机断层扫描（CT）衍生的低于?950 HU（Hounsfield unit）的低衰减区百分比）的关联，模型校正年龄和性别，并使用HC3异方差稳健标准误。采用Benjamini-Hochberg错误发现率（BH-FDR）控制多重检验。显著位点被合并为唯一成熟miRNA。二次分析考察其与FEV₁（第一秒用力呼气容积）占预计值百分比、FEV₁/FVC（用力肺活量）和一氧化碳弥散量（DL_CO）的关联。敏感性分析排除强影响观测、使用对数转换的肺气肿百分比，并在既往吸烟者中校正包年。研究人员使用multiMiR查询经实验验证的miRNA-靶标相互作用，随后进行Gene Ontology Biological Process（GO-BP）富集分析。结果：14个miRNA在BH-FDR<0.05水平与肺气肿百分比相关。8个miRNA与肺气肿呈正相关（miR-144-5p、miR-190a、miR-20b-5p、miR-451a、miR-3912、miR-15a-5p、miR-202-5p和miR-4497），6个呈负相关（miR-489、miR-181b-5p、miR-378a-5p、miR-3614-5p、miR-197-3p和miR-877-3p）。全部14个miRNA与肺功能的关联与肺气肿结果一致。经表型特异性BH-FDR校正后，12个miRNA与FEV₁占预计值百分比相关，12个与FEV₁/FVC相关，8个与DL_CO相关。所有14个关联在敏感性分析中均保持方向与显著性。对实验支持的靶标进行富集分析，识别出与miRNA调控、凋亡、缺氧/氧反应、细胞周期调控和血管发育相关的过程。结论：14个成熟miRNA的肺组织表达与COPD患者的肺气肿百分比和肺功能相关；而基于实验支持相互作用数据库的靶标富集结果为探索性，需要独立和功能验证。关键词：microRNA；COPD；肺气肿；肺功能。

研究背景：慢性阻塞性肺疾病（COPD，chronic obstructive pulmonary disease）为异质性疾病，持续性气流受限和结构异常是其核心特征，其中肺实质肺气肿性破坏是疾病表型的主要决定因素之一。即使肺功能损害相近，患者肺气肿范围仍差异显著，且肺气肿负担越重，肺功能和功能状态越差。定量计算机断层扫描（CT，computed tomography）测量肺气肿较视觉或定性分级更客观、可重复，可实现与分子和生理数据的连续组织水平关联，而非仅按“有/无肺气肿”分组比较。既往研究提示COPD人肺组织miRNA表达改变，并涉及组织损伤、修复、细胞外基质调节和气道生物学等通路；部分miRNA与肺气肿严重度或肺气肿性肺破坏相关。然而，肺组织miRNA表达与COPD患者内连续定量肺气肿的关系仍未完全阐明。既往研究多比较不同肺病组、使用肺气肿的生理替代指标，或在包含非COPD人群的人群中评估定量肺气肿。因此，研究人员开展本研究，旨在识别COPD个体中与定量肺气肿相关的肺组织miRNA，并检验这些miRNA与肺功能指标是否呈一致关联；同时通过敏感性分析评估关联稳健性，并探索实验支持的靶基因及富集的生物学过程，为观察到的miRNA-肺气肿关联提供可能的生物学背景。研究得出14个成熟miRNA与肺气肿百分比及肺功能相关，靶标富集涉及凋亡、缺氧/氧反应、细胞周期调控和血管过程，但结论具有假设生成性，需独立和功能验证。

关键技术方法概括（不超过250字）：研究人员基于公共Gene Expression Omnibus数据集GSE201121，纳入117例有定量肺气肿资料的COPD/肺气肿患者，匹配肺组织miRNA表达与临床及定量影像表型。表达过滤后保留在至少50%样本中表达值>0的839个miRNA位点；以年龄和性别校正线性回归分析标准化miRNA表达与肺气肿百分比（低于?950 HU的CT低衰减区百分比）的关联，采用HC3异方差稳健标准误，并以BH-FDR控制多重检验。二次分析采用年龄和性别校正模型评估14个成熟miRNA与FEV₁占预计值百分比、FEV₁/FVC及DL_CO的关联。敏感性分析排除强影响观测（n=109）、采用对数转换肺气肿百分比，并在既往吸烟者中校正年龄、性别和包年。研究人员使用multiMiR查询实验验证的miRNA-靶标相互作用，并按miRNA与肺气肿关联方向分组，使用clusterProfiler进行GO-BP富集分析。

研究结果：
研究人群（Study Population）：该横断面分析纳入117例COPD患者；平均年龄65.5岁（SD，10.0），男性71例（60.7%），既往吸烟者104例（88.9%），平均包年58.2；平均FEV₁为50.1%预计值，平均DL_CO为55.9%预计值，肺气肿百分比中位数7.57%（IQR，1.95–28.32%）。通过描述性分析给出研究人群基线特征。

与定量肺气肿相关的miRNA（miRNAs Associated with Quantitative Emphysema）：在839个筛查miRNA位点中，年龄和性别校正线性回归显示15个位点（代表14个唯一成熟miRNA）达到BH-FDR<0.05。8个成熟miRNA与肺气肿百分比正相关：miR-144-5p、miR-190a、miR-20b-5p、miR-451a、miR-3912、miR-15a-5p、miR-202-5p和miR-4497；6个负相关：miR-489、miR-181b-5p、miR-378a-5p、miR-3614-5p、miR-197-3p和miR-877-3p。校正后效应估计范围为每1-SD更高表达对应?6.10个百分点（miR-378a-5p）至+6.70个百分点（miR-144-5p）。

与肺功能的关联（Associations with Pulmonary Function）：对14个肺气肿相关miRNA分别以FEV₁占预计值百分比、FEV₁/FVC和DL_CO为结局进行年龄和性别校正线性回归，所有14个miRNA与三项肺功能指标的关联方向均与肺气肿关联相反，即与更高肺气肿相关的miRNA倾向于与更差肺功能相关，反之亦然。经表型特异性BH-FDR校正后，12/14个miRNA与FEV₁占预计值百分比相关，12/14个与FEV₁/FVC相关，8/14个与DL_CO相关。

敏感性分析（Sensitivity Analyses）：排除8个强影响观测后（n=109），全部14个成熟miRNA保持相同关联方向且均保持名义显著（P<0.05），回归系数中位绝对变化约10.2%；以对数转换肺气肿百分比为结局时，14个miRNA仍保持相同方向且名义显著；在仅限既往吸烟者的敏感性分析中，校正年龄、性别和包年后，14个miRNA仍保持相同方向且名义显著，回归系数相对主模型中位绝对变化为4.7%。通过预设敏感性分析，研究人员得出主关联总体稳健的结论。

靶标分析与生物学过程富集（Target Analysis and Biological Process Enrichment）：基于实验支持的miRNA-靶标分析共识别140个唯一基因，其中89个与肺气肿正相关miRNA相连，55个与负相关miRNA相连，4个基因（BCL2、PTEN、MAPK1和TMED7）为两个方向基因集共享。GO-BP富集显示，与更高肺气肿相关miRNA的靶标富集于miRNA转录和代谢调控、细胞周期G1/S转换、凋亡信号、肌肉细胞增殖和血管发育等过程；与更低肺气肿相关miRNA的靶标富集于凋亡、缺氧和氧水平反应及JNK信号等过程。研究人员将这些富集结果视为探索性，提示观察到的miRNA-肺气肿关联可能存在相关生物学背景。

讨论总结：研究人员指出，该研究在COPD患者中评估肺组织miRNA表达与定量肺气肿的关系，发现14个唯一成熟miRNA，其中8个正相关、6个负相关；这些关联不仅限于CT定义的肺气肿，全部14个miRNA与FEV₁占预计值百分比、FEV₁/FVC和DL_CO的关联方向均与肺气肿结果一致，经表型特异性FDR校正后分别有12、12和8个miRNA与上述肺功能指标相关。实验支持靶标的功能注释涉及凋亡、缺氧和氧水平反应、细胞周期调控及miRNA调控通路。鉴于横断面设计，研究人员认为这些发现属于假设生成而非验证性，不能解释为这些miRNA因果性促进或防止肺气肿性破坏。研究扩展了既往在异质性肺疾病中研究miRNA表达的工作，聚焦COPD内定量肺气肿，系统筛查肺组织miRNA与连续CT肺气肿表型的关联，并进一步在FEV₁占预计值百分比、FEV₁/FVC和DL_CO中检验，显示结构肺气肿的分子关联伴随方向一致的生理关联。DLCO一致关联可能与肺气肿性破坏降低肺泡-毛细血管气体交换面积并伴肺微血管损伤有关；实验研究提示失调miRNA可通过Notch和VEGF信号调节肺微血管内皮存活和血管生成。与既往研究比较，miR-15a-5p曾在COPDGene全血miRNA-mRNA网络中作为与肺气肿百分比相关的枢纽，miR-202-5p在重度早发COPD患者肺组织中报告升高；其他肺模型中对miR-190a、miR-451、miR-181b-5p、miR-378和miR-197-3p的发现方向兼容。但一致性并非普遍，例如本研究 miR-144-5p较高表达与更重肺气肿相关，而实验性miR-144过表达在LPS诱导急性肺损伤中保护内皮屏障破坏和血管渗漏；本研究miR-489较高表达与更轻肺气肿相关，而抑制miR-489在新生高氧模型中改善肺泡分隔。研究人员认为这些差异提示miRNA效应具有情境依赖性，可能反映疾病过程、发育阶段、细胞类型、组织组成和成熟miRNA链特异性差异，因此不宜对单个miRNA统一赋予致病或保护角色。富集分析中，负相关miRNA靶标富集于凋亡相关通路和缺氧反应，正相关miRNA靶标富集于miRNA调控、细胞周期控制、凋亡信号和血管发育；凋亡相关通路在两个方向组均出现，可能反映miRNA复杂调控效应。局限性包括：横断面设计不能确认miRNA表达改变先于肺气肿还是继发于结构肺损伤；bulk肺组织反映多种肺细胞群信号，肺气肿相关细胞组成差异可能部分贡献关联；尽管控制多重比较且敏感性分析一致，仍可能存在吸入治疗等残余混杂、样本量小、单一数据集限制精确性和推广性，需更大独立肺组织队列重复；靶标和富集分析依赖既往实验支持的miRNA-靶标相互作用及数据库实验背景，富集不能证明特定miRNA在人肺气肿肺组织中直接调控、激活或因果参与；本研究未包括qPCR分子验证、实验扰动或独立验证。总体意义在于提供与COPD定量肺气肿及肺功能相关的候选miRNA和可能生物学过程，为后续细胞特异性表达、mRNA谱、纵向表型和功能实验提供依据。

研究结论（翻译）：总之，14个成熟miRNA的肺组织表达与COPD个体的定量肺气肿相关。这些miRNA还与肺功能呈方向一致的关联，支持其与结构和功能疾病特征的相关性。富集分析提示可能涉及凋亡、缺氧相关反应、细胞周期调控和血管过程。这些发现具有假设生成性，需要独立重复和功能验证。