《International Journal of Genomics》：Single Cell Transcriptomic Atlas Reveals Cellular Heterogeneity and Molecular Mechanisms in Pediatric Airway Hyperresponsiveness Epithelium
编辑推荐：
背景：气道高反应性（airway hyperresponsiveness, AHR）是儿童慢性炎性气道疾病的重要临床特征，但其细胞结构仍未完全阐明。由于批量转录组学（bulk transcriptomics）方法可能掩盖来自稀有或状态特异性细胞群体的信号，单细
背景：气道高反应性（airway hyperresponsiveness, AHR）是儿童慢性炎性气道疾病的重要临床特征，但其细胞结构仍未完全阐明。由于批量转录组学（bulk transcriptomics）方法可能掩盖来自稀有或状态特异性细胞群体的信号，单细胞RNA测序（single-cell RNA sequencing, scRNA-seq）为研究高反应性气道中的上皮与免疫细胞多样性提供了有用框架。
方法：研究人员对来源于儿童AHR患者支气管刷检样本的公共scRNA-seq数据（GEO: GSM8722221，10x Genomics Chromium平台）进行处理，使用Cell Ranger进行质控过滤后，采用PCA（15个主成分）、UMAP/tSNE嵌入及Leiden聚类（resolution=0.8）进行分析。细胞类型通过经典标志物注释并经Seurat标签转移验证。伪时间轨迹通过Monocle（DDRTree）推断，并经scVelo、Slingshot和PAGA验证。巨噬细胞亚群在resolution=1.2下解析并通过hdWGCNA分析。GO和KEGG富集分析使用clusterProfiler（Benjamini-Hochberg校正；校正p<0.05）。网络药理学整合TCM BTMAN 2.0、Swiss Target Prediction、SEA及PPI拓扑分析，以识别小青龙汤（Xiao Qing Long Tang）的靶点与草药化合物。
结果：质量控制支持该数据集的分析适用性（线粒体比例与nCount_RNA的r=-0.31；数据集间r=0.96）。Actb、Gapdh、Tuba1b、Hmga2、Vim、Col1a1、Zeb1、Snai2、Twist1和Mmp9等为高变基因，经JackStraw检验后保留前15个主成分。UMAP和tSNE解析出八个具有协调分离的气道黏膜细胞群体。TSLP、IL33、CCL26、S100A8、POSTN、CLCA1、FCGR3B、TNFRSF8、FOXJ1和SCGB1A1的细胞类型特异性表达与已知的炎症、重塑和黏液纤毛程序一致。高分辨率巨噬细胞分析鉴定出23个推测性亚群，研究人员谨慎地将其解读为计算定义的细胞状态。富集分析突出细胞质翻译、核糖体生物合成、氧化磷酸化及相关的线粒体特征。网络药理学鉴定VEGFA、PTGS2和ESR1为优先靶点，疾病相关与化合物相关靶点集合之间产生83个共享候选靶点。
结论：该研究提供了儿童AHR单细胞图谱以及上皮-免疫异质性、巨噬细胞状态多样性和候选治疗网络的整合性、假说生成性视角。研究结果应被解读为需要健康儿童对照、独立队列和功能实验验证的计算证据。
本研究以儿童气道高反应性（airway hyperresponsiveness, AHR）为研究对象，利用公共数据库中的单细胞RNA测序（scRNA-seq）数据构建细胞水平的转录组图谱，并结合网络药理学方法探索潜在治疗靶点。该论文发表于《International Journal of Genomics》。
**一、研究背景与问题提出**
儿童AHR是全球约3亿患者罹患的 prevalent 慢性呼吸系统疾病，其发病率在不同地理和社会经济背景下持续上升，临床特征表现为气道炎症、支气管高反应性及进行性结构重塑三联征，具体症状包括反复发作的呼吸困难、喘息、胸闷和持续性咳嗽。尽管研究者已在AHR病理生理学方面积累了数十年基础，但高反应性气道的固有细胞异质性、遗传易感性位点与环境损伤及免疫回路失调之间的复杂相互作用，仍阻碍了机制研究向有效个体化治疗的转化。
呼吸道黏膜是宿主免疫系统与吸入环境之间的主要界面，在物理屏障维持、黏液纤毛清除及先天性和适应性免疫调节中发挥着不可或缺的作用。在AHR患者中，该黏膜层经历深刻的病理学转变，包括上皮屏障破坏、病理性黏液高分泌、纤毛搏动频率受损及上皮修复机制失调。这些分子和细胞改变是AHR核心临床特征的关键机制，并可能涉及多个相互依赖细胞类型之间时空协调的改变，但这些细胞类型的精确贡献仍远未阐明。
传统的批量RNA测序（bulk RNA-sequencing）虽提供了AHR相关基因表达的基础转录图谱，但该方法因对异质性细胞群体进行平均而从根本上掩盖了疾病相关信号，因RNA提取方案引入系统性偏倚，且无法捕捉细胞间通讯至关重要的微环境背景。单细胞RNA测序（scRNA-seq）则变革了研究者以单细胞分辨率探究疾病过程的能力，能够识别稀有细胞群体、精确表征细胞状态转变，并从快照转录数据推断发育轨迹。在气道疾病背景下，scRNA-seq已超越批量方法的局限，解析出上皮亚型特异性转录程序，列举了离子细胞（ionocyte）和多纤毛前体细胞（deuterosomal cell）等功能相关的稀有群体，并揭示了批量平均所掩盖的细胞间信号轴。
尽管新兴的单细胞研究已开始揭示哮喘及相关气道疾病的显著细胞多样性，但全面绘制儿童AHR气道黏膜完整细胞生态系统并整合功能注释的单细胞图谱仍然明显缺乏。本研究即针对这一空白，通过对儿童AHR气道黏膜组织的详细scRNA-seq分析，旨在构建高反应性气道黏膜的综合细胞图谱，明确细胞类型特异性分子标志，阐明上皮分化轨迹和谱系关系，高分辨率表征巨噬细胞亚群多样性和功能特化，识别疾病相关转录网络及上游调控机制，并整合网络药理学进行治疗靶点优先排序。
**二、主要技术方法**
本研究所用数据来源于公共基因表达综合数据库（Gene Expression Omnibus, GEO），登录号为GSM8722221，系采用10x Genomics Chromium平台从儿童AHR患者支气管刷检样本中获取的转录组数据，刷检样本主要采集气道表面上皮，亦可回收黏膜和气道腔中的免疫细胞。数据处理采用Cell Ranger进行参考基因组比对、基因水平UMI计数和细胞条形码去多重化，经PCA降维（保留15个主成分）、UMAP/tSNE嵌入和Leiden聚类（整体数据集resolution=0.8，巨噬细胞亚群resolution=1.2）进行细胞分群。细胞类型通过经典谱系标志物注释，并经Seurat标签转移和差异表达分析验证。伪时间轨迹使用Monocle（DDRTree方法）重建，以基底细胞（basal cell）为根状态，并经scVelo、Slingshot和PAGA交叉验证。巨噬细胞亚群分析采用hdWGCNA进行基因共表达网络分析。功能富集分析使用clusterProfiler进行GO和KEGG分析（Benjamini-Hochberg校正，校正p<0.05）。网络药理学分析整合TCM BTMAN 2.0、Swiss Target Prediction和SEA数据库，构建蛋白-蛋白互作（PPI）网络，并通过拓扑学参数（度中心性、中介中心性和接近中心性）进行靶点优先排序。数据整合采用Harmony算法校正批次效应，通过k近邻批次效应检验（kBET）、局部逆辛普森指数（LISI）和轮廓系数评估校正效果。
**三、研究结果**
**从原始转录组到可解释的细胞景观：儿童AHR scRNA-seq数据的预处理与降维分析**：质量控制显示线粒体基因比例与总RNA计数呈适度负相关（r=-0.31），而检测基因数与总计数之间显示出强内部一致性（r=0.96），表明低质量条形码和明显技术异常值在降维和聚类前已得到充分控制。Actb、Gapdh、Tuba1b、Hmga2、Vim、Col1a1、Zeb1、Snai2、Twist1和Mmp9等为高变基因，经JackStraw置换检验支持前15个主成分的统计相关性（所有p<0.05）。
**儿童气道高反应性黏膜组织的单细胞图谱构建与细胞分类**：UMAP降维显示细胞分离为离散区域，tSNE呈现相似的分离模式，表明主要细胞群不依赖于单一可视化方法。基于标志物的注释鉴定出主要气道黏膜群体，包括上皮、内皮、基底、杯状、纤毛和免疫富集区室。杯状细胞和巨噬细胞富集群体的存在与AHR中的黏液失调和先天性免疫激活一致。
**空间解析的基因表达特征界定儿童AHR黏膜组织的细胞身份与炎症回路**：TSLP和IL33在限制性上皮相关区域富集，与其在2型气道炎症中作为上游警报素（alarmin）的已知作用一致；CCL26表达提示与嗜酸性粒细胞募集的可能关联；S100A8和FCGR3B标记免疫富集区域，分别与髓系和中性粒细胞相关活性相符；POSTN显示更广泛的基质样分布，与重塑相关的细胞外基质信号一致；CLCA1集中于紧凑的分泌性上皮簇；FOXJ1和SCGB1A1分别支持纤毛细胞和俱乐部细胞（club cell）的注释。
**儿童气道高反应性黏膜组织的谱系承诺与转录特化**：差异表达分析支持细胞类型注释框架，微血管内皮细胞、纤毛细胞、基底细胞和杯状细胞显示独特的标志物图谱。Monocle伪时间分析将基底细胞置于推断的根状态附近，提示从基底样状态向更特化上皮群体的逐步转录转变。层次聚类将共表达基因划分为不同模块，各代表在上皮谱系承诺特定阶段活跃的协调生物学程序或调控回路。
**解构巨噬细胞多样性：儿童气道高反应性的亚群结构与谱系特异性功能**：高分辨率单细胞分析鉴定出23个巨噬细胞相关亚群（M1-M23），经层次聚类组织为促炎型、抗炎型和组织重塑型等更广泛功能类别。UMAP可视化显示部分亚群占据紧凑区域提示共享转录程序，而其他亚群更为分散，可能代表中间或情境依赖性激活状态。这些空间模式为AHR微环境中巨噬细胞功能多样性提供线索，但各亚群的生物学独立性需通过外部数据集和实验检测确证。
**炎症与气道重塑交叉点：儿童AHR细胞的通路富集与网络水平功能分析**：GO生物学过程分析显示细胞质翻译、核糖体组装、ATP代谢过程、核糖核蛋白复合物生物合成、氧化磷酸化、有氧呼吸和核糖体生物合成等通路的主导富集；GO细胞组分分析提示胞质核糖体、核糖体亚基、黏着斑和细胞-基质连接；GO分子功能分析识别核糖体结构成分、钙黏蛋白结合、rRNA结合和NADH脱氢酶活性。KEGG通路分析鉴定核糖体、冠状病毒病-COVID-19、帕金森病、氧化磷酸化、亨廷顿病、朊病毒病、化学致癌作用-活性氧（ROS）、糖尿病心肌病、肌萎缩侧索硬化症（ALS）和阿尔茨海默病为最显著富集通路。但需指出，KEGG神经退行性疾病术语通常由共享氧化磷酸化、线粒体和核糖体基因驱动，故应主要解读为代谢和氧化应激相关转录活性的证据，而非神经退行性机制直接参与儿童AHR。
**网络药理学与儿童气道高反应性的靶点识别**：通过单细胞结果作为生物学起点，网络药理学分析汇总了来自公共数据库的疾病相关靶基因，使用TCM BTMAN 2.0及Swiss Target Prediction和SEA检索小青龙汤的化合物-靶点关联。疾病靶点与化合物靶点的交集产生83个共享治疗候选靶点。PPI网络构建后，通过拓扑测量优先排序VEGFA、PTGS2、ESR1、AKT1、TP53、TNF、IL6、STAT3、EGFR和PIK3CA等靶点。麻黄（MA HUANG）和半夏（BAN XIA）在所检测草药中显示最广的靶点覆盖度。层次聚类热图显示半夏、干姜（GAN JIANG）和麻黄具有更广泛的预测靶点覆盖，而五味子（WU WEI ZI）、炙甘草（ZHI GAN CAO）和桂枝（GUI ZHI）则更集中于选定基因群周围。UpSet图总结草药间靶点重叠，最大交集包含23个共享靶点，提示配方内可能存在互补性靶点覆盖，但这些计算关联应被视为进一步对接、扰动和体内验证的候选对象。
**四、讨论与结论**
本研究提供了儿童AHR黏膜的细胞分辨率图谱，将已知的AHR相关炎症、上皮、重塑和巨噬细胞相关信号组织到细胞解析框架中，主要价值在于以儿科单细胞背景精细化组织已知AHR相关程序。
基底细胞在图谱中的显著代表与其作为气道上皮多能祖细胞的建立角色一致，伪时间分析提示的高反应性条件下基底细胞转录状态改变，可能反映慢性炎症气道中修复程序的异常激活。TSLP和IL33的差异表达为警报素驱动的上皮-免疫交互提供分子证据，其空间协调性表达提示协同警报素信号放大下游2型免疫级联。CCL26和FCGR3B的表达模式突出嗜酸性和嗜中性粒细胞炎症的互补机制。POSTN和CLCA1在特定细胞亚群中的富集表达指向活跃的上皮下重塑和杯状细胞驱动的黏液病理。23个巨噬细胞相关亚群的鉴定提示儿童AHR气道黏膜内存在显著的先天性免疫异质性，但高分辨率聚类可能过度分割连续转录梯度，M1-M23标签最好理解为计算定义的巨噬细胞状态。
GO和KEGG富集分析显示突出的细胞质翻译、核糖体生物合成、氧化磷酸化和线粒体通路特征，表明分析气道细胞中生物合成和代谢活动增强，与炎症、上皮修复和黏液产生的能量需求一致。网络药理学鉴定VEGFA、PTGS2和ESR1为疾病相关与化合物相关网络交集内的优先靶点，但从小青龙汤靶点预测到scRNA-seq发现的转化仍属计算性质，缺乏分子对接、扰动实验、动物模型或临床验证。
研究的主要局限在于基于单一公共GEO数据集，未纳入匹配的健康儿童对照或独立验证队列，限制了区分疾病特异性变化与基线变异的能力。伪时间轨迹、巨噬细胞极化状态和通路水平解读仍为计算推断；支气管刷检样本主要代表黏膜表面，可能未完全捕获更深的气道区室。网络药理学部分仅提供靶点优先排序，未证实结合、通路调节或治疗效力。
综上所述，本研究呈现了儿童AHR黏膜的细胞分辨率图谱，突出了上皮炎症信号、巨噬细胞相关转录多样性、代谢通路活性和候选化合物-靶点网络。研究发现将在健康对照、独立队列和功能模型中的验证情况，决定其转化价值。
