《Frontiers in Genetics》：Single-cell transcriptomic profiling uncovers key molecular signatures in glioma pathogenesis
编辑推荐：
背景：胶质瘤是中枢神经系统最常见且致命的原发性恶性肿瘤，其特点是显著的细胞异质性和不良预后。在单细胞分辨率下全面表征胶质瘤对于识别新的治疗靶点和改善患者预后至关重要。方法：研究人员从基因表达综合数据库（GEO）获取胶质瘤样本，进行了全面的单细胞RNA测序（sc
背景：胶质瘤是中枢神经系统最常见且致命的原发性恶性肿瘤，其特点是显著的细胞异质性和不良预后。在单细胞分辨率下全面表征胶质瘤对于识别新的治疗靶点和改善患者预后至关重要。方法：研究人员从基因表达综合数据库（GEO）获取胶质瘤样本，进行了全面的单细胞RNA测序（scRNA-seq）分析。采用了先进的计算方法，包括主成分分析（PCA）、统一流形逼近与投影（UMAP）、t-分布随机邻域嵌入（t-SNE）和差异表达分析，以表征细胞异质性。进行了基因共表达网络构建和通路富集分析以识别功能模块。候选生物标志物在GL261胶质瘤细胞、C8-D1A正常胶质细胞和椎间盘髓核细胞中通过定量实时荧光PCR（qRT-PCR）和酶联免疫吸附试验（ELISA）进行了验证。结果：质量控制分析揭示了高质量的单细胞数据，每个细胞的中位基因计数为5,437，UMI计数为14,207。共鉴定出4,753个高度可变基因，并使用Louvain算法划分了22个不同的细胞簇。PCA负载分析识别了关键贡献基因，包括细胞周期调节因子（CDK1、TOP2A、BIRC5）、免疫相关基因（C1QA、C1QB、HLA-DRB6）和神经谱系标记物（MOBP、MAG、GJB1）。细胞类型注释确定了七种主要细胞群：星形胶质细胞、少突胶质细胞、小胶质细胞、神经干细胞、OPC/不成熟神经元、周细胞和T细胞。差异表达分析发现了847个上调基因和652个下调基因（|log2FC| > 1，校正P < 0.05）。基因共表达网络分析揭示了五个以枢纽基因为中心的模块，包括CLEC12A、CLU、AQP4和S100A16。通路富集分析显示，细胞周期、Notch信号通路、MAPK通路和神经发生显著参与。实验验证证实，与正常胶质细胞相比，Plp1在GL261细胞中显著下调（0.38 ± 0.05倍，P < 0.01），而Fth1（2.15 ± 0.28倍，P < 0.001）和Gm42418（5.67 ± 0.52倍，P < 0.001）则显著上调。结论：这项全面的单细胞转录组分析成功表征了胶质瘤的细胞异质性，识别了不同的细胞群和分子特征。PLP1、FTH1和GM42418被验证为潜在的分子生物标志物，为理解胶质瘤发病机制和潜在治疗靶点提供了新的见解。
胶质瘤是中枢神经系统最常见且最具破坏性的原发性恶性肿瘤，约占所有恶性脑肿瘤的80%，其特点是显著的细胞异质性、侵袭性的生物学行为以及在综合治疗干预下仍然不良的预后。尽管过去几十年在外科手术、放射治疗和化疗方案方面取得了实质性进展，但最具侵袭性的亚型——胶质母细胞瘤的中位生存期仍约为15个月，这凸显了对新型诊断生物标志物和治疗靶点的迫切需求。近年来，转录组学分析技术的突破性进展，特别是单细胞RNA测序（scRNA-seq）技术的出现，使研究人员能够在前所未有的分辨率上全面表征肿瘤微环境内的细胞异质性，揭示了以前未被识别的细胞亚群、发育层级和驱动肿瘤发生、发展及治疗抵抗的转录程序，从而推动了我们对胶质瘤生物学的理解。胶质瘤的肿瘤微环境包含一个由肿瘤细胞、免疫细胞、神经细胞和基质成分组成的复杂生态系统。理解其细胞组成和分子特征对于开发有效的治疗策略至关重要。然而，将生物信息学发现转化为临床相关生物标志物，需要从单细胞转录组学到蛋白质水平的跨平台系统验证。为应对这些挑战，研究人员利用基因表达综合数据库（GEO，数据集编号：GSM9116755）的公开单细胞转录组数据，结合系统性的实验验证，对胶质瘤的细胞复杂性进行了全面的单细胞RNA测序分析。通过降维、差异表达谱分析、基因共表达网络构建和通路富集分析等先进的计算分析，研究人员识别了在细胞群间表现出独特表达模式的候选分子生物标志物，并在建立的胶质瘤细胞系（GL261）和正常胶质细胞（C8-D1A）中通过定量实时荧光PCR（qRT-PCR）和酶联免疫吸附试验（ELISA）对这些发现进行了验证，为识别的分子特征的生物学相关性提供了有力证据。这项研究结果发表于《Frontiers in Genetics》期刊。
研究人员为开展研究用到的主要技术方法包括：首先，从GEO数据库获取单细胞RNA测序原始数据，并进行严格的质量控制，过滤低质量细胞和潜在双联体。其次，利用Seurat软件包对数据进行标准化，识别高度可变基因，并通过主成分分析（PCA）、统一流形近似与投影（UMAP）和t-分布随机邻域嵌入（t-SNE）进行降维和可视化。接着，应用Louvain算法进行无监督聚类，鉴定出不同的细胞群，并通过经典标记基因表达模式进行细胞类型注释。此外，利用FindAllMarkers函数进行差异表达分析，构建基因共表达网络并识别功能模块，并利用clusterProfiler进行通路富集分析。最后，在GL261胶质瘤细胞、C8-D1A正常胶质细胞和椎间盘髓核细胞模型中，采用qRT-PCR和ELISA技术对候选生物标志物（Plp1、Fth1、Gm42418）的表达进行了实验验证。
研究结果部分内容总结如下：
**单细胞RNA测序质量控制与初步聚类分析**：通过小提琴图评估每个细胞的基因计数（nFeature_RNA）和UMI总计数（nCount_RNA），证实了数据的高质量，中位数分别为5,437和14,207。密度散点图显示大多数细胞线粒体基因百分比（percent.mt）低于5%，表明细胞状态健康。通过方差稳定变换（vst）方法鉴定了4,753个高度可变基因。对标准化后的表达矩阵进行主成分分析（PCA），PCA负载热图显示，前几个主成分的主要贡献基因包括细胞周期相关基因（如CDK1、TOP2A、BIRC5）、免疫相关基因（如C1QA、C1QB、HLA-DRB6）和神经谱系标记物（如MOBP、MAG、GJB1）。基于前30个主成分进行UMAP分析和Louvain聚类，成功划分出22个不同的细胞簇（Clusters 0-21）。
**通过多维基因表达分析表征细胞簇异质性**：t-SNE可视化结合密度轮廓线确认了22个细胞群的清晰空间分隔。树状图显示了细胞簇大小的分布差异。UMAP特征图、小提琴图、山脊图和标记基因表达点图展示了特定标记基因（如A2M、CLEC12A、OPALIN、EGFR等）在各细胞簇中的表达模式，揭示了簇特异性的转录特征。细胞复杂性分析显示，每个簇的平均检测到的基因数量存在异质性，簇3的转录复杂性最高，簇7相对较低。
**细胞类型注释与质量指标分析**：基于经典标记基因的表达模式进行细胞类型注释。Z-score归一化的热图显示了每个簇的独特表达特征。通过堆叠条形图量化并可视化了跨簇的细胞类型组成，最终鉴定出七种主要细胞群：星形胶质细胞、少突胶质细胞、小胶质细胞、神经干细胞、OPC/不成熟神经元、周细胞和T细胞。质量指标箱线图显示大多数簇具有一致的nFeature_RNA和nCount_RNA分布，表明聚类由真实的转录变异驱动，而非技术伪影。细胞复杂性散点图证实了UMI计数与基因检测数量之间的正相关性（Pearson r = 0.433）。
**差异表达、基因共表达网络及通路富集分析**：火山图可视化了差异表达分析结果，共鉴定出847个显著上调和652个显著下调基因（|log2FC| > 1，校正P < 0.05）。基因共表达网络分析揭示了以枢纽基因为中心的协调表达基因模块，枢纽基因包括CLEC12A、ATP13A4、CLU、AQP4、S100A16等。三维PCA可视化提供了细胞异质性的综合空间表示。通路富集点图显示，差异表达基因显著富集于细胞周期、Notch信号通路、MAPK通路和神经发生等与胶质瘤相关的通路。细胞类型表达相关性热图分析了肿瘤微环境中不同细胞群之间的关系，例如T细胞与祖细胞、神经干细胞之间呈负相关。
**候选生物标志物的实验验证**：在GL261胶质瘤细胞中，通过qRT-PCR验证了三个候选生物标志物的表达。结果显示，髓鞘蛋白脂质蛋白基因Plp1显著下调（0.38 ± 0.05倍，P < 0.01），铁蛋白重链基因Fth1显著上调（2.15 ± 0.28倍，P < 0.001），而长链非编码RNA基因Gm42418的上调幅度最大（5.67 ± 0.52倍，P < 0.001）。椎间盘髓核细胞中这三个基因的表达水平介于正常胶质细胞和胶质瘤细胞之间。ELISA实验也从蛋白质水平对这些发现进行了验证。
**细胞转录复杂性与基因调控网络架构**：通过六边形密度图进一步确认了UMI计数与基因计数之间强烈的正相关性（Pearson r = 0.89），证实了数据的技术可重复性。差异基因表达分析识别了炎症标记物（如Cxcl10）和细胞外基质重塑因子（如Mmp9）等显著上调的基因。基因共表达网络分析确定了与免疫反应、细胞周期调节、细胞外基质组织、神经元信号传导和代谢过程相关的五个主要功能模块。三维UMAP可视化整合了簇分离、表达梯度和细胞大小指标，显示某些簇之间可能存在渐进的转录状态转换，而另一些簇则保持离散的身份。
讨论与结论总结：
本研究通过全面的单细胞转录组分析，成功表征了胶质瘤的细胞异质性，识别出22个具有独特分子特征的细胞群，代表了肿瘤微环境中的七种主要细胞类型。主成分分析（PCA）负载分析揭示细胞周期调节因子、免疫标记物和神经谱系基因是细胞异质性的主要贡献者。基因共表达网络分析鉴定出以CLEC12A、CLU和AQP4等基因为中心的功能模块。通路富集分析表明细胞周期、Notch信号通路、MAPK通路和神经发生在胶质瘤发病机制中起着关键作用。实验验证证实了PLP1在胶质瘤细胞中下调，而FTH1和GM42418上调。这些发现为理解胶质瘤生物学提供了新的见解，并确定了用于诊断和治疗开发的潜在分子生物标志物。特别是，FTH1的上调可能与胶质瘤细胞增强的抗铁死亡防御机制有关，而GM42418（一种长链非编码RNA）的显著上调是一个值得深入研究的潜在新靶点。研究也存在一定局限性，如分析基于单一样本数据，未来需要在多中心独立队列中进行验证，以及需要进行功能研究以阐明这些生物标志物在胶质瘤发病机制中的具体作用。
