协调的线粒体基因组表达程序标志低级别胶质瘤中分化神经谱系和突触转录状态

《International Journal of Molecular Sciences》：A Coordinated Mitochondrial Genome Expression Program Marks a Differentiated Neural-Lineage and Synaptic Transcriptional State in Lower-Grade Glioma

　　

线粒体调控细胞代谢、信号传导和细胞命运，然而，在低级别胶质瘤（LGG）中，协调的线粒体基因组转录的生物学意义仍不清楚。研究人员整合了来自TCGA和CGGA队列的转录组学和临床数据，以研究13个线粒体DNA编码的蛋白编码基因（mtPCGs）的表达。所有mtPCGs的表达升高始终与有利的分子亚型和延长的总生存期（OS）相关。全转录组相关性分析和基因本体论（Gene Ontology）分析揭示了一个富集于突触信号传导、神经传递、突触后组织和神经元分化的协调程序，同时伴有与神经元和突触标记物的正相关，以及与干性和免疫相关标记物的负相关。尽管线粒体转录评分在调整已建立的临床分子变量后并不与生存独立相关，但它始终标记出生物学上侵袭性较低的肿瘤。协调的线粒体基因组表达识别出一种与有利LGG生物学相关的神经元样转录状态，揭示了线粒体转录、突触身份和胶质瘤分化之间的意外联系。

**论文解读：协调的线粒体基因组表达程序标志低级别胶质瘤中分化神经谱系和突触转录状态**

**研究背景与问题**

低级别胶质瘤（LGG）是一组弥漫性浸润性胶质瘤，包括世界卫生组织（WHO）2级和3级肿瘤。现代分类整合了组织病理学特征与分子标志物，特别是异柠檬酸脱氢酶（IDH）突变状态和1p/19q联合缺失状态。这些分子改变定义了生物学上不同的胶质瘤亚型，其临床行为、治疗反应和患者预后显著不同。线粒体除了产生ATP外，还参与细胞信号传导、分化和细胞死亡调控，但其在LGG中的生物学意义尚不完全清楚。线粒体DNA（mtDNA）编码13个蛋白编码基因（mtPCGs），这些基因是氧化磷酸化（OXPHOS）系统的核心组分。已有研究表明，mtDNA拷贝数、突变以及线粒体基因转录抑制可影响癌症中的关键信号通路，但协调的线粒体基因组转录在LGG中的生物学意义仍不清楚。因此，研究人员旨在通过整合TCGA（癌症基因组图谱）和CGGA（中国胶质瘤基因组图谱）两大独立队列的转录组学和临床数据，系统表征13个mtPCGs在LGG中的表达模式、预后意义及其全转录组相关性，以揭示其潜在生物学功能。

**研究内容与结论**

研究人员发现，所有13个mtPCGs的表达升高在TCGA和CGGA队列中均一致地与有利的分子亚型（IDH突变且1p/19q共缺失型）和延长的总生存期（OS）显著相关。尽管在调整年龄、肿瘤分级、IDH状态和1p/19q状态后，线粒体转录评分与生存的关联不再具有独立统计学意义，但其始终标记出生物学上侵袭性较低的肿瘤。全转录组相关性分析和基因本体论（GO）富集分析揭示了一个与线粒体评分高度相关的协调转录程序，该程序显著富集于化学突触传递、神经元突触可塑性调节、电压门控钠通道簇集和逆行神经元致密核心囊泡运输等突触和神经元过程。此外，线粒体评分与神经元和突触标记物（如VAMP2、SHANK1、SYP）呈正相关，与少突胶质细胞谱系标记物（如MBP、PLP1）呈弱正相关，而与免疫细胞标记物（如C1QA、PTPRC）和干性标记物（如NES、SOX9）呈负相关。这些发现表明，协调的线粒体基因组表达与一种分化神经谱系和突触转录状态相关，这种状态与更有利的LGG生物学行为相关联。

**重要意义**

该研究发表在《International Journal of Molecular Sciences》。研究首次聚焦于线粒体基因组自身编码的13个mtPCGs，而非核编码的线粒体相关基因，揭示了协调的线粒体基因组表达程序与LGG中神经元样分化状态之间的意外联系。这一发现为理解线粒体转录在胶质瘤分化中的作用提供了新视角，并提示线粒体基因表达可能作为识别生物学上惰性肿瘤的转录状态标志物，而非独立的预后因子。研究结果为后续探索线粒体功能与神经谱系分化之间的机制关联奠定了基础。

**关键技术与方法**

研究人员使用了来自TCGA（https://www.cancer.gov/tcga）和CGGA（http://www.cgga.org.cn）两个独立队列的RNA-seq数据及临床信息，包括WHO 2级和3级胶质瘤样本。主要技术方法包括：1）利用DESeq2进行方差稳定化变换和基因-wise z-score标准化，计算每个样本的13个mtPCGs表达的算术平均值作为线粒体转录评分；2）采用Kaplan-Meier生存分析和log-rank检验评估mtPCGs表达与总生存期（OS）的关联，并使用survminer包确定最佳截断值；3）通过多变量Cox比例风险回归模型调整年龄、肿瘤分级、IDH突变状态和1p/19q共缺失状态，评估线粒体评分的独立预后价值；4）使用Spearman秩相关分析进行全转录组基因与线粒体评分的相关性分析，并选取前200个正相关基因进行基因本体论（GO）富集分析（Enrichr平台）；5）计算线粒体评分与代表性细胞谱系标志物（神经元、突触、少突胶质细胞、免疫、干性等）的Spearman相关系数。

**研究结果**

**2.1 mtPCGs在不同LGG类型中的表达**
研究人员首先根据组织学分类和分子亚型分析mtPCGs表达模式。在TCGA和CGGA数据集中，按组织学类型分层时，未观察到显著总体差异；但按分子亚型分析时，发现明确的表达模式：IDH突变且1p/19q共缺失的LGG（LGG-IDH-mut-codel）中mtPCGs表达最高，IDH野生型（LGG-IDH-wt）表达最低，而IDH突变但无1p/19q共缺失（LGG-IDH-mut-non-codel）表达居中。这一模式在两个队列中均一致。

**2.2 高mtPCG表达与LGG患者更长OS相关**
Kaplan-Meier分析显示，所有13个mtPCGs的高表达在TCGA和CGGA队列中均与更长的总生存期（OS）显著相关。将13个mtPCGs表达合并为线粒体基因评分后，该评分的高值也与更长生存期相关。但多变量Cox回归分析表明，在调整年龄、肿瘤分级、IDH状态和1p/19q共缺失状态后，线粒体评分与OS的关联不再具有统计学显著性，提示其预后价值依赖于已知的分子预后因素。

**2.3 mtPCGs表达之间的相关性**
通过计算13个mtPCGs之间的两两相关性，发现大多数基因之间存在中等到强的正相关（Spearman或Pearson相关系数），表明存在一个协调的基因表达程序。层次聚类热图直观展示了这一共表达结构。

**2.4 线粒体转录程序的功能表征**
利用全转录组相关性分析，研究人员识别出与线粒体评分最正相关的200个基因。GO富集分析显示，这些基因显著富集于化学突触传递、神经元突触可塑性调节、电压门控钠通道簇集和逆行神经元致密核心囊泡运输等神经元和突触过程。分子功能分析进一步突出通道活性和可溶性N-乙基马来酰亚胺敏感因子附着蛋白受体（SNARE）结合，支持与神经元信号和囊泡运输的关联。代表性正相关基因包括ATP1A3、VAMP2、SHANK1等。

**2.5 线粒体评分与细胞谱系标志物的关联**
线粒体评分与神经元和突触标志物（VAMP2 r=0.44, SHANK1 r=0.44, SYP r=0.39, SNAP25 r=0.36, RBFOX3 r=0.35）以及少突胶质细胞谱系标志物（MBP r=0.25, PLP1 r=0.21, MAG r=0.19, MOG r=0.17）呈正相关；而与免疫标志物（C1QA r=-0.26, PTPRC r=-0.20, AIF1 r=-0.17）和干性标志物（NES r=-0.27, SOX9 r=-0.15）呈负相关。

**讨论与结论**

讨论部分指出，本研究利用大规模癌症基因组学数据集，发现13个mtPCGs的协调表达与LGG的有利分子亚型和更长生存期密切相关，但并非独立预后因子，而是反映了一种与肿瘤分化和低侵袭性相关的更广泛生物学表型。该线粒体转录程序与突触和神经元过程富集，以及神经元和少突胶质细胞谱系标志物正相关、免疫和干性标志物负相关，表明协调的线粒体基因表达可能标志一种分化神经谱系状态。研究人员提出，增加的mtPCG表达可能为高度耗能的突触传递和神经可塑性提供能量支持，从而赋予肿瘤更成熟的神经元样表型。研究局限性包括：转录组关联不能直接推断因果关系；mRNA丰度与蛋白质水平可能不一致；未能评估单个mtPCG在不同患者或WHO分级中的过表达频率；以及数据来自bulk肿瘤组织，不能明确归因于特定细胞类型。未来需结合单细胞RNA测序、空间转录组学或功能实验进一步验证。

**结论部分翻译**：我们的研究结果表明，13个mtPCGs的协调表达在LGG中始终与有利的分子亚型和更长的生存期相关。除了其预后关联外，线粒体转录程序还与涉及神经元和突触生物学过程的基因相关联，并与神经元、突触和少突胶质细胞谱系标志物呈正相关，而与免疫、血管和干性相关标志物呈负相关。这些发现共同提示，升高的线粒体基因表达发生于一种更广泛的分化神经谱系和突触转录状态中，该状态与更有利的LGG生物学行为相关。尽管线粒体评分在调整已确定的分子预后因素后并非独立预后因子，但其与有利临床分子特征的可重复关联凸显了协调的线粒体基因表达作为LGG转录异质性中一个具有生物学信息组成部分的作用。