Single-cell RNA Sequencing: A New Perspective in the Study of Chronic Obstructive Pulmonary Disease
-
摘要: 慢性阻塞性肺疾病(chronic obstructive pulmonary disease,COPD)是一种受遗传和环境因素交互作用影响的多基因疾病,其特征为持续的气流受限和慢性呼吸道症状。传统的基因测序方法虽然揭示了COPD的遗传学基础,但受限于细胞群体平均数据,难以完全展现疾病的复杂性。单细胞RNA测序技术在个体细胞层面分析基因表达,揭示细胞异质性,对理解COPD的病理过程、发现生物标志物及个性化医疗具有重要意义。综述了单细胞测序技术的原理以及其在COPD研究中的应用。Abstract: Chronic obstructive pulmonary disease (COPD) is a polygenic disease influenced by the interaction of genetic and environmental factors, characterized by persistent airflow limitation and chronic respiratory symptoms. Traditional gene sequencing methods have revealed the genetic basis of COPD but are limited by the average data of cell populations, making it difficult to fully display the complexity of the disease. Single-cell RNA sequencing technology provides a new perspective, allowing researchers to analyze gene expression at the individual cell level, revealing cellular heterogeneity, which is of great significance for understanding the pathological process of COPD, discovering biomarkers, and personalized medicine. This article reviews the principles of single-cell sequencing technology and its application in COPD research.
-
动脉粥样硬化(atherosclerosis,AS)是动脉粥样硬化性脑梗死(atherosclerotic cerebral infarction,ASCI)的病理基础,包括斑块形成、血管狭窄等,前者是AS典型标志,后者则是AS进展的表现[1−2]。数字减影血管造影术是血管类疾病的首选诊断手段,诊断准确性较高,但引起有创性质,可能损伤患者神经系统,引发并发症,临床应用受限[3]。颈内动脉、椎-基底动脉是前循环、后循环的主要供血动脉,经颅多普勒超声(transcranial doppler ultrasound,TCD)及颈动脉超声可反映椎-基底动脉、颈内动脉血供状态,对ASCI患者血管狭窄诊断有一定作用[4−5]。研究表明,甲状腺功能减退经损伤血管内皮功能,介导氧化应激、炎症反应,参与并诱导AS[5−6]。TCD、颈动脉超声及甲状腺功能指标(甲功3项)单一评估ASCI颈内血管狭窄程度的价值有限,基于上述,本研究创新性联合TCD、颈动脉超声及甲功3项评估ASCI颈内血管狭窄程度,为临床针对性制定方案提供参考,现报告如下。
1. 资料与方法
1.1 一般资料
本研究回顾性选取2021年10月至2023年9月邯郸市第一医院收治的200例ASCI患者临床资料。纳入标准:(1)ASCI符合《中国急性缺血性脑卒中诊治指南2018》[7]诊断;(2)经磁共振成像确诊;(3)血管狭窄经数字减影血管造影术明确诊断;(4)入院后1周内接受TCD、颈动脉超声及甲功3项检查;(5)两侧颞窗透声良好;(6)近期无抗凝药物使用史。排除标准:(1)脑出血、脑外伤及脑肿瘤等;(2)伴有甲状腺疾病史;(3)存在自身免疫性疾病;(4)伴有凝血功能异常;(5)存在心肺等脏器功能障碍;(6)近期有手术史;(7)蒙特利尔认知评估异常。最终纳入182例ASCI患者。本研究经邯郸市第一医院伦理委员会审批通过(2024-L-32)。患者均签署知情同意书。
1.2 研究方法
1.2.1 TCD检查
采用彩色多普勒超声诊断仪(德国DWLDoppler-Box),探头频率2.0 MHz,取仰卧位,经颞窗、枕窗检查大脑中动脉、大脑前动脉、大脑后动脉和椎-基底动脉,记录流向、平均血流速度、声频信号、频谱形态等信息,最终获取椎动脉血流搏动指数(pulse index,PI)、阻力指数(resistance index,RI)参数。
1.2.2 颈动脉超声检查
采用彩色多普勒超声诊断仪(Vivid E9),线阵探头频率:3~9 MHz。取仰卧位,颈后垫枕后仰暴露颈部,沿颈动脉走行探查颈内动脉外颅外段,记录内膜回声、颈内动脉距起始部20 mm处管径(diameter,D),测定颈内动脉收缩期最大流速(peak systolic velocity,PSV)、舒张末期流速(end-diastolic volume,EDV)参数。
1.2.3 甲功3项水平检测
采集受检者4 mL空腹静脉血,离心(
3000 r/min,15 min)分离血清,使用全自动生化分析仪(罗氏-601)检测血清促甲状腺激素(thyroid stimulating hormone,thyrotropin,TSH)、游离甲状腺素(free thyroxine,FT4)、游离三碘甲状腺原氨酸(free triiodothyronine,FT3)水平。1.2.4 血管狭窄程度
通过数字减影血管造影(digital subtraction angiography,DSA)检查分级[8]:直径狭窄率 < 50%为轻度狭窄,直径狭窄率50%~69%为中度狭窄,直径狭窄率70%~99%为重度狭窄,直径狭窄率99%~100%为完全闭塞。
1.3 统计学分析
采用SPSS27.0对数据进行分析,计量资料采取Bartlett方差齐性检验与夏皮罗-威尔克正态性检验,均确认具备方差齐性且服从正态分布,以均数±标准差($ \bar x \pm s $)表示,多组间比较以单因素方差进行分析,组间两两对比采用LSD-t检验,计数资料以n(%)表示,差异行χ2检验,以pearson相关系数法分析超声相关参数与甲功3项的关联性,采用受试者工作特征(receiver operating characteristic,ROC)曲线分析TCD、颈动脉超声、甲功3项单独及联合评估ASCI血管狭窄的价值,联合评估任一阳性即为阳性。检验水准α = 0.05。
2. 结果
2.1 DSA检查结果及不同血管狭窄程度ASCI患者临床资料
182例ASCI患者中男92例,女90例,年龄54~72岁,平均(63.08±4.29)岁,体质量指数20~27 kg/m2,平均(23.16±1.45)kg/m2;发病至就诊时间1~6 h,平均(3.80±0.79)h,基础疾病:高血压101例,高脂血症89例,糖尿病60例,饮酒49例,吸烟36例。
经DSA证实动脉血管轻度狭窄64例(35.16%),中度狭窄69例(37.91%),重度狭窄29例(15.93%),完全闭塞20例(10.99%)。不同血管狭窄程度各组临床资料比较,差异无统计学意义(P > 0.05),见表1。
表 1 不同血管狭窄程度ASCI患者临床资料比较[($ \bar x \pm s $)/n(%)]Table 1. Comparison of clinical data of ASCI patients with different degrees of vascular stenosis [($ \bar x \pm s $)/n(%)]项目 轻度狭窄组(n = 64) 中度狭窄组(n = 69) 重度狭窄组(n = 29) 完全闭塞组(n = 20) χ2/F P 性别 1.506 0.681 男 29(45.31) 36(52.17) 15(51.72) 12(60.00) 女 35(54.69) 33(47.83) 14(48.28) 8(40.00) 年龄(岁) 63.58 ± 4.09 62.29 ± 4.16 64.11 ± 3.92 62.72 ± 4.02 1.848 0.140 体质量指数(kg/m2) 23.18 ± 1.17 22.96 ± 1.43 23.48 ± 1.21 23.32 ± 1.29 1.248 0.294 发病至就诊时间(h) 3.72 ± 0.63 3.86 ± 0.72 3.75 ± 0.69 3.91 ± 0.75 0.685 0.563 高血压 27(50.00) 36(57.97) 22(72.41) 16(75.00) 6.455 0.092 高脂血症 25(39.06) 33(47.83) 17(58.62) 14(60.00) 4.512 0.211 糖尿病 16(25.00) 20(28.99) 13(44.83) 11(45.00) 5.486 0.140 饮酒 13(20.31) 16(23.19) 11(37.93) 9(45.00) 7.109 0.071 吸烟 10(15.63) 12(17.39) 7(24.14) 7(35.00) 4.211 0.240 2.2 不同血管狭窄程度ASCI患者超声相关参数、甲功3项比较
不同血管狭窄程度各组RI、PI、D、PSV、EDV、TSH、FT4、FT3比较,差异有统计学意义(P < 0.05),PSV、EDV、TSH、FT4随着颈内血管狭窄程度加重呈升高趋势,RI、PI、D、FT3随着颈内血管狭窄程度加重呈降低趋势(P < 0.05),见表2。
表 2 不同血管狭窄程度ASCI患者超声相关参数、甲功3项比较($ \bar x \pm s $)Table 2. Comparison of ultrasound related parameters and thyroid function test in ASCI patients with different degrees of vascular stenosis ($ \bar x \pm s $)项目 轻度狭窄组(n = 64) 中度狭窄组(n = 69) 重度狭窄组(n = 29) 完全闭塞组(n = 20) χ2/F P TCD参数 RI 0.60 ± 0.02 0.59 ± 0.03△ 0.56 ± 0.03△# 0.51 ± 0.02△#▲ 70.732 < 0.001* PI 0.92 ± 0.03 0.86 ± 0.05△ 0.74 ± 0.04△# 0.70 ± 0.03△#▲ 229.374 < 0.001* 颈内动脉超声参数 D(mm) 4.83 ± 1.06 3.79 ± 0.97△ 3.15 ± 0.82△# 2.63 ± 0.74△#▲ 37.855 < 0.001* PSV(cm/s) 171.68 ± 51.94 226.45 ± 63.67△ 323.86 ± 78.92△# 0.00 ± 0.00△#▲ 128.582 < 0.001* EDV(cm/s) 56.76 ± 8.64 85.38 ± 20.21△ 147.85 ± 43.10△# 0.00 ± 0.00△#▲ 207.157 < 0.001* 甲功3项 TSH(mIU) 7.22 ± 0.61 9.24 ± 0.98△ 12.85 ± 1.36△# 13.91 ± 1.68△#▲ 320.630 < 0.001* FT4(pmol/L) 15.41 ± 1.03 17.10 ± 1.99△ 18.13 ± 2.45△# 19.84 ± 2.72△#▲ 33.191 < 0.001* FT3(pmol/L) 2.25 ± 0.58 1.97 ± 0.43△ 1.67 ± 0.26△# 1.45 ± 0.23△#▲ 21.083 < 0.001* *P < 0.05;与轻度狭窄组比较,△P < 0.05;与中度狭窄组比较,#P < 0.05;与重度狭窄组比较,▲P < 0.05。 2.3 超声相关参数与甲功3项的相关性分析
TSH与PSV、EDV呈正相关,与RI、PI、D呈负相关,FT4与PSV、EDV呈正相关,与RI、PI、D呈负相关,FT3与RI、PI、D呈正相关,与PSV、EDV呈负相关(P < 0.05),见表3。
表 3 超声相关参数与甲功3项的相关性分析Table 3. Correlation analysis of ultrasound parameters and thyroid function tests超声相关参数 TSH FT4 FT3 r P r P r P TCD参数 RI −0.637 < 0.01* −0.606 0.01* −0.613 < 0.01* PI −0.704 < 0.01* −0.582 0.012* 0.642 < 0.01* 颈动脉超声参数 D −0.643 < 0.01* −0.624 < 0.01* 0.588 0.007* PSV 0.716 < 0.01* 0.701 < 0.01* −0.631 < 0.01* EDV 0.683 < 0.01* 0.634 < 0.01* −0.592 0.004* *P < 0.05。 2.4 TCD、颈动脉超声、甲功3项单独及联合评估ASCI血管狭窄程度的效能
将ASCI患者TCD、颈动脉超声、甲功3项参数作为检验变量,颈内血管狭窄程度作为状态变量,轻中度为阴性、重度与完全闭塞为阳性,绘制ROC曲线结果显示,TCD、颈动脉超声、甲功3项联合评估ASCI血管狭窄程度的AUC值为0.941,大于TCD、颈动脉超声、甲功3项单独评估的AUC(P < 0.05),见表4~表6、图1。
表 4 轻中度、重度与完全闭塞ASCI患者超声相关参数、甲功3项比较($ \bar x \pm s $)Table 4. Comparison of ultrasound parameters and thyroid function tests in ASCI patients with mild to moderate,severe,and complete occlusion ($ \bar x \pm s $)项目 轻中度狭窄(n=133) 重度狭窄与完全闭塞(n=49) t P TCD参数 RI 0.59 ± 0.02 0.54 ± 0.03 12.956 < 0.001* PI 0.89 ± 0.06 0.72 ± 0.05 17.691 < 0.001* 颈内动脉超声参数 D(mm) 4.29 ± 1.15 2.94 ± 0.89 7.433 < 0.001* PSV(cm/s) 200.09 ± 65.18 323.86 ± 78.92 10.716 < 0.001* EDV(cm/s) 71.61 ± 22.64 147.85 ± 43.10 15.456 < 0.001* 甲功3项 TSH(mIU) 8.27 ± 1.21 13.28 ± 1.80 21.537 < 0.001* FT4(pmol/L) 16.29 ± 2.38 18.83 ± 2.94 5.981 < 0.001* FT3(pmol/L) 2.10 ± 0.62 1.58 ± 0.33 5.580 < 0.001* *P < 0.05。 表 5 TCD、颈动脉超声、甲功3项单独及联合评估ASCI血管狭窄程度的效能Table 5. Efficacy of TCD,carotid ultrasound,and thyroid function tests in assessing ASCI vascular stenosis individually and in combination指标 AUC 95%CI 截断值 敏感度(%) 特异度(%) P RI 0.755 0.686~0.816 0.58 77.55 68.42 < 0.01* PI 0.750 0.681~0.812 0.78 75.51 65.41 < 0.01* TCD联合 0.853 0.793~0.901 12.13 71.43 85.71 < 0.01* D 0.784 0.717~0.841 3.99 mm 75.51 74.44 < 0.01* PSV 0.764 0.695~0.823 267.8 cm/s 77.55 61.65 < 0.01* EDV 0.743 0.673~0.804 158.48 cm/s 63.27 79.70 < 0.01* 颈动脉超声联合 0.873 0.815~0.917 12.18 83.67 81.20 < 0.01* TSH 0.810 0.745~0.864 12.69 mIU 79.59 76.69 < 0.01* FT4 0.787 0.720~0.844 16.51 pmol/L 79.59 66.17 < 0.01* FT3 0.741 0.671~0.803 2.12 pmol/L 71.43 66.92 < 0.01* 甲功3项联合 0.871 0.813~0.916 14.03 79.59 87.22 < 0.01* TCD、颈动脉超声、甲功3项联合 0.941 0.896~0.970 12.73 89.80 87.22 < 0.01* *P < 0.05。 表 6 ROC曲线比较Table 6. Comparison of ROC curve成对对比 AUC差异 标准误差 95%CI Z P TCD联合VS TCD、颈动脉超声、甲功3项联合 0.088 0.035 0.019~0.157 2.512 0.012* 颈动脉超声联合VS TCD、颈动脉超声、甲功3项联合 0.068 0.034 0.002~0.135 2.020 0.043* 甲功3项联合VS TCD、颈动脉超声、甲功3项联合 0.073 0.037 0.002~0.145 2.002 0.045* *P < 0.05。 3. 讨论
3.1 TCD及颈动脉超声与ASCI血管狭窄的关联
AS、血管狭窄或闭塞是脑梗死的病理基础,ASCI发病急骤、病情危重,发生后颈内血管狭窄闭塞可引起脑组织缺氧、缺血,脑神经损伤,进而造成脑组织坏死[11]。超声技术是临床诊断血管病变、评价AS的重要方法,因其无创安全、可重复、价格低廉等优点,可为临床早期治疗提供影像学依据[12]。TCD可利用超声脉冲效应提供脑动脉实时血流状况,及时监测诊断脑血管病变;频谱分析参数PI是反映血管搏动性与血管顺应性的重要指标,PI值降低提示血管顺应性变差,血管腔扩张难度增加[13−15]。RI主要体现血流灌注阻力,具有双向变化,一般在缺氧代偿期呈下降趋势,而缺氧失代偿期呈升高趋势[16]。本研究结果显示随着颈内血管狭窄程度加重,RI、PI呈降低趋势,说明可利用RI、PI参数辅助评估ASCI血管狭窄。分析其原因为ASCI患者颈内血管狭窄程度越严重,对血管内皮细胞的损害越大,血流灌注阻力越强,进一步降低血管顺应性,外周动脉明显阻塞,RI、PI显著下降[17]。椎-基底动脉、颈内动脉系统是大脑的主要供血来源,椎-基底动脉为后循环血供,颈内动脉为前循环血供,其血流供应情况与ASCI患者血管狭窄程度相关[18]。颈动脉超声可用超声波特性无创、直观观察个体颈部血管血流动力学,评估血管狭窄程度,获取颈动脉D、EDV、PSV等参数[19−20]。
3.2 甲功3项与ASCI血管狭窄的关联
临床对ASCI深入研究发现,脂质代谢障碍是AS的病变基础,且甲状腺功能减退与脂质代谢障碍有关,参与ASCI、心血管疾病进展[21−23]。FT4、FT3具有生物活性,可直接反映甲状腺功能状态,其水平不受甲状腺结合蛋白影响[24]。TSH有促垂体合成作用,参与甲状腺球蛋白合成、酪氨酸碘化,维持机体甲状腺功能[25]。以往研究证实,TSH水平升高与AS有关[26−27]。本研究显示,不同血管狭窄程度各组TSH、FT4随着颈内血管狭窄程度加重呈升高趋势,FT3随着颈内血管狭窄程度加重呈降低趋势。分析主要由于ASCI患者脑结构与功能异常改变,血管狭窄导致脑循环障碍,加重下丘脑-垂体-甲状腺轴代谢紊乱,导致TSH、FT4升高。TSH通过激活巨噬细胞、诱导线粒体呼吸链功能障碍,促进细胞炎症因子表达、自由基生成,加重血管斑块炎症,促进AS进展。随着血清FT4升高、FT3降低,可减少低密度脂蛋白胆固醇受体mRNA的表达,减弱腺苷环化酶作用,抑制脂肪分解,进而加重脂代谢异常,参与AS的发生[28−29]。本研究相关性分析显示甲功3项与血管狭窄程度显著相关。李阔等[30]研究显示,ASCI患者FT4、TSH高表达,FT3低表达,与血管狭窄程度、有无斑块有关,可佐证本研究。另外本研究发现甲功3项与TCD及颈动脉超声参数存在明显相关性,说明ASCI患者血管狭窄。陈园园等[8]研究显示,TCD、颈动脉超声联合检查评估老年脑梗死血管病变的准确度为88.33%,高于两种超声单独诊断。本研究在TCD及颈动脉超声基础上补充甲功3项常规检测,发现RI、PI、D、PSV、EDV、TSH、FT4、FT3联合评估ASCI血管狭窄程度的AUC值最大,为0.941,说明从超声血流动力学、甲状腺功能等不同方面评估ASCI血管狭窄可进一步提高诊断价值。
综上所述,ASCI血管狭窄与甲功3项、甲功3项与TCD及颈动脉超声参数之间显著相关,RI、PI、D、PSV、EDV、TSH、FT4、FT3联合评估ASCI血管狭窄的价值较高,可为不同血管狭窄程度ASCI患者临床治疗提供参考。
-
表 1 单细胞测序在COPD研究中的部分应用案例
Table 1. Some Application Cases of Single-cell Sequencing in COPD
研究内容 研究对象 主要发现 细胞亚群差异分析 COPD患者VS健康对照组 在上皮细胞、内皮细胞、免疫细胞中
均表现不同的基因表达及分子功能差异。COPD内型识别 COPD患者 COPD不同内型对环境暴露、临床表现、
疾病进展、治疗反应性等方面均存在明显异质性。COPD严重程度评估 COPD肺组织 IL-18R、CD69及穿孔素表达上调代表疾病的进展,
可评估COPD严重程度。发现治疗靶点 COPD患者及小鼠模型 SOCS1和HSPB1可能是COPD治疗新靶点;肠道微生物能通过
调节肺组织炎症反应稳定COPD肺免疫力。 -
[1] Singh D,Agusti A,Anzueto A,et al. Global Strategy for the Diagnosis,Management,and Prevention of Chronic Obstructive Lung Disease: the GOLD science committee report 2019[J]. Eur Respir J,2019,53(5): 1900164 [pii]. [2] Ishii T,Nishimura M,Akimoto A,et al. Understanding low COPD exacerbation rates in Japan: a review and comparison with other countries[J]. Int J Chron Obstruct Pulmon Dis,2018,13:3459-3471. doi: 10.2147/COPD.S165187 [3] Yang T,Cai B,Cao B,et al. Severity distribution and treatment of chronic obstructive pulmonary disease in China: baseline results of an observational study[J]. Respir Res,2022,23(1):106. doi: 10.1186/s12931-022-02021-w [4] Cazzola M,Calzetta L,Rinaldi B,et al. Management of Chronic Obstructive Pulmonary Disease in Patients with Cardiovascular Diseases[J]. Drugs,2017,77(7):721-732. doi: 10.1007/s40265-017-0731-3 [5] 周超平,李鑫辉. 单细胞转录组研究进展[J]. 生物化学与生物物理进展,2013,40(12):1193-1200. [6] Tang F,Barbacioru C,Wang Y,et al. mRNA-Seq whole-transcriptome analysis of a single cell[J]. Nat Methods,2009,6(5):377-382. doi: 10.1038/nmeth.1315 [7] Wang J,Song Y. Single cell sequencing: a distinct new field[J]. Clin Transl Med,2017,6(1):10. doi: 10.1186/s40169-017-0139-4 [8] 刘芳远,苏秀兰. 单细胞测序技术在干细胞领域的研究进展[J]. 中国细胞生物学学报,2020,42(2):328-335. [9] Hwang B,Lee J H,Bang D. Single-cell RNA sequencing technologies and bioinformatics pipelines[J]. Exp Mol Med,2018,50(8):1-14. [10] 孙楠,曲守方,陈样宜,等. 高通量基因测序仪行业标准的验证[J]. 分子诊断与治疗杂志,2021,13(7):1027-1030. doi: 10.3969/j.issn.1674-6929.2021.07.003 [11] 张成鹏. 单细胞测序技术及其应用综述[J]. 河南科学,2022,40(9):1390-1397. [12] Wang C,Xu J,Yang L,et al. Prevalence and risk factors of chronic obstructive pulmonary disease in China (the China Pulmonary Health [CPH] study): A national cross-sectional study[J]. Lancet,2018,391(10131): 1706-1717 [13] 郑烯,胡雪峰. 慢性阻塞性肺疾病发病机制的研究进展[J]. 中国细胞生物学学报,2019,41(2):304-311. [14] Liao S,Chen Y. The Role of Bioactive Small Molecules in COPD Pathogenesis[J]. COPD,2024,21(1):2307618. doi: 10.1080/15412555.2024.2307618 [15] 盛梅梅,罗专波. 慢性阻塞性肺疾病发病机制的研究现状及进展[J]. 现代医学与健康研究电子杂志,2023,7(1):38-41. [16] Ragland M F,Benway C J,Lutz S M,et al. Genetic advances in chronic obstructive pulmonary disease. Insights from COPD gene[J]. Am J Respir Crit Care Med,2019,200(6):677-690. doi: 10.1164/rccm.201808-1455SO [17] Sakornsakolpat P,Prokopenko D,Lamontagne M,et al. Genetic landscape of chronic obstructive pulmonary disease identifies heterogeneous cell-type and phenotype associations[J]. Nat Genet,2019,51(3):494-505. doi: 10.1038/s41588-018-0342-2 [18] Zhang X, Guo Y, Yang J, et al. A functional variant alters binding of activating protein 1 regulating expression of FGFZ gene associated with chronic obstructive pulmonary disease[J]. BMC Med Genet,2019,20(1):33. [19] Castaldi P J,Guo F,Qiao D,et al. Identification of functional variants in the FAM13A Chronic Obstructive Pulmonary Disease genome-wide association study locus by massively parallel reporter assays[J]. Am J Respir Crit Care Med,2019,199(1):52-61. doi: 10.1164/rccm.201802-0337OC [20] Xiao X,Cai W,Ding Z,et al. A2M Serves as promising biomarker for Chronic Obstructive Pulmonary Disease[J]. Int J Chron Obstruct Pulmon Dis,2023,18:683-692. doi: 10.2147/COPD.S407491 [21] Li Y N,Shi G C,Guan J. Advances in chest imaging in early chronic obstructive pulmonary disease[J]. Zhonghua Jie He He Hu Xi Za Zhi,2023,46(12):1266-1271. [22] Li X,Noell G,Tabib T,et al. Single cell RNA sequencing identifies IGFBP5 and QKI as ciliated epithelial cell genes associated with severe COPD[J]. Respir Res,2021,22(1):100. doi: 10.1186/s12931-021-01675-2 [23] An L,Xia H,Zheng W,et al. Comparison of cell subsets in chronic obstructive pulmonary disease and controls based on single-cell transcriptome sequencing[J]. Technol Health Care,2023,31(S1):9-24. [24] Sauler M,McDonough J E,Adams T S,et al. Characterization of the COPD alveolar niche using single-cell RNA sequencing[J]. Nat Commun,2022,13(1):494. doi: 10.1038/s41467-022-28062-9 [25] Bhatt S P,Agusti A,Bafadhel M,et al. Phenotypes,Etiotypes,and Endotypes of Exacerbations of Chronic Obstructive Pulmonary Disease[J]. Am J Respir Crit Care Med,2023,208(10):1026-1041. doi: 10.1164/rccm.202209-1748SO [26] Barnes P J. Inflammatory endotypes in COPD[J]. Allergy,2019,74(7):1249-1256. doi: 10.1111/all.13760 [27] Singh D,Kolsum U,Brightling C E,et al. Eosinophilic inflammation in COPD: Prevalence and clinical characteristics[J]. Eur Respir J,2014,44(6):1697-700 doi: 10.1183/09031936.00162414 [28] Huang Y,Niu Y,Wang X,et al. Identification of novel biomarkers related to neutrophilic inflammation in COPD[J]. Front Immunol,2024,15:1410158. doi: 10.3389/fimmu.2024.1410158 [29] Wohnhaas C T,Gindele J A,Kiechle T,et al. Cigarette smoke specifically affects small airway epithelial cell populations and triggers the expansion of inflammatory and squamous differentiation associated basal cells[J]. Int J Mol Sci,2021,22(14):7646. doi: 10.3390/ijms22147646 [30] Huang Q,Wang Y,Zhang L,et al. Single-cell transcriptomics highlights immunological dysregulations of monocytes in the pathobiology of COPD[J]. Respir Res,2022,23(1):367. doi: 10.1186/s12931-022-02293-2 [31] Villaseñor-Altamirano A B,Jain D,Jeong Y,et al. Activation of CD8(+) T cells in chronic obstructive pulmonary disease lung[J]. Am J Respir Crit Care Med,2023,208(11):1177-1195. doi: 10.1164/rccm.202305-0924OC [32] Hu Y,Shao X,Xing L,et al. Single-cell sequencing of lung macrophages and monocytes reveals novel therapeutic targets in COPD[J]. Cells,2023,12(24):2771. doi: 10.3390/cells12242771 [33] 薛廷,郑乐婷,董菲,等. 单细胞RNA测序联合实验验证树突状细胞在慢性阻塞性肺疾病中的核心基因[J]. 细胞与分子免疫学杂志,2024,40(2):97-105. [34] Fan P,Zhang Y,Ding S,et al. Integrating RNA-seq and scRNA-seq to explore the mechanism of macrophage ferroptosis associated with COPD[J]. Front Pharmacol,2023,14:1139137. doi: 10.3389/fphar.2023.1139137 [35] Lai H C,Lin T L,Chen T W,et al. Gut microbiota modulates COPD pathogenesis: Role of anti-inflammatory Parabacteroides goldsteinii lipopolysaccharide[J]. Gut,2022,71(2):309-321. [36] 王聪,刘富强,王军奎. 吸烟对肠道菌群影响的研究进展[J]. 实用心脑肺血管病杂志,2022,30(9):28-31. [37] Hicks S C,Townes F W,Teng M,Irizarry R A. Missing data and technical variability in single-cell RNA-sequencing experiments[J]. Biostatistics,2018,19(4):562-578. [38] Zhang Y,Xia R,Lv M,et al. Machine-learning algorithm-based prediction of diagnostic gene biomarkers related to immune infiltration in patients with chronic obstructive pulmonary disease[J]. Front Immunol,2022,13:740513. -