周启涛生物技术工作室
基因组与遗传2025/11/03· 编译自 Nature Genetics

IF31.7|线粒体rRNA热点突变在癌症中的功能显性效应

约4%肿瘤携带线粒体核糖体RNA热点突变,低至10%异质性即可引发线粒体功能障碍。

约4%肿瘤携带线粒体核糖体RNA热点突变,低至10%异质性即可引发线粒体功能障碍。

这篇按生信解读号的读法拆:先看背景和问题,再看作者怎么做、做出了哪几条硬发现,最后把数据和代码收走。

研究背景

肿瘤中反复出现的体细胞突变(热点)主要影响核基因组蛋白编码基因,但线粒体DNA(mtDNA)编码自身翻译机器,包括两个线粒体核糖体RNA(mt-rRNA)和22个线粒体tRNA。mt-tRNA突变可导致氧化磷酸化缺陷,且特定mt-tRNA热点已在癌症中被发现。然而,mt-rRNA是否在癌症中受到正选择尚不清楚。此前最大规模的肿瘤mtDNA研究样本量仅为本研究的百分之一到十分之一,统计功效不足。mtDNA是多拷贝基因组,突变表型取决于异质性水平。传统观点认为,致病性mtDNA突变通常为功能隐性,需超过60%异质性才表现临床表型。

科学问题

卡在哪里: 此前缺乏大规模肿瘤mtDNA测序数据,无法系统识别mt-rRNA热点突变;且mt-rRNA突变在癌症中的功能后果和显隐性特征完全未知。

本文要回答: 利用14,106例肿瘤全基因组测序数据,系统描绘mtDNA热点突变图谱,鉴定mt-rRNA热点,并通过碱基编辑和单细胞多组学解析其功能显性效应。

技术路线

作者主要用了:WGS、热点检测算法、DdCBE碱基编辑、单细胞多组学、代谢组、蛋白质组

作者首先对14,106例肿瘤全基因组测序数据进行体细胞mtDNA变异检测,应用基于背景突变率的热点检测算法识别显著复发位点。随后通过跨队列验证、种系约束分析和结构注释,评估mt-rRNA热点的选择压力和功能潜力。为验证功能,利用DdCBE碱基编辑在143B细胞中引入m.1227G>A突变,结合代谢组、蛋白质组和单细胞多组学,定量解析异质性依赖的分子表型。

Fig. 1 · 原文图,按文末许可署名使用
Fig. 1 · 原文图,按文末许可署名使用

图 1 Fig. 1。图1展示14,106例肿瘤的mtDNA热点全景:a为组织类型分布;b显示热点沿mtDNA位置分布,rRNA基因区域富集;c按突变类型分类患者比例,约4%携带rRNA热点;d显示各基因热点数量,MT-RNR1/2仅次于MT-ND5;e将热点分为蛋白编码、rRNA、tRNA等类别。

核心亮点

亮点 1|mt-rRNA是癌症热点突变的新靶点

在14,106例肿瘤中,共鉴定出138个显著SNV热点,其中34个位于MT-RNR1和MT-RNR2基因,仅次于MT-ND5。约4%的肿瘤携带至少一个rRNA热点突变。这些热点在TCGA和PCAWG队列中高度重现,排除了队列特异性假象。与蛋白编码热点不同,rRNA热点缺乏组织特异性,且异质性水平不显著升高。

Fig. 2 · 原文图,按文末许可署名使用
Fig. 2 · 原文图,按文末许可署名使用

图 2 Fig. 2。图2聚焦rRNA热点特征:a显示前四大组织谱系中各热点突变比例;b为异质性分布,rRNA热点异质性不显著升高;c展示种系同质性变异极少,表明强纯化选择;d显示热点更易破坏Watson-Crick配对;e为种系约束比较;f-g为mitoribosome结构上热点位置;h-i为tRNA/mRNA接触位点聚类。

亮点 2|热点位于种系强纯化选择位点

rRNA热点优先发生在种系受限位置(OR=2.77, P=0.049),且进化保守性显著高于非热点rRNA突变(P=0.0097)。结构分析显示,热点突变聚集于mitoribosome的tRNA接触位点(A位点、P位点)和mRNA通道入口,可能破坏翻译功能。

Fig. 3 · 原文图,按文末许可署名使用
Fig. 3 · 原文图,按文末许可署名使用

图 3 Fig. 3。图3展示m.1227G>A功能验证:a-b为m.1227G在mtSSU中的位置和氢键网络;c-e为DdCBE筛选和编辑流程;f显示mtDNA拷贝数不变;g为呼吸链亚基蛋白表达下降;h为核糖体沉降分析显示组装异常;i为耗氧率降低;j-m为代谢组变化;n-q为同位素示踪显示代谢流改变;r-s为蛋白质组显示呼吸链复合物下调。

亮点 3|m.1227G>A在低异质性下引发线粒体功能障碍

通过DdCBE碱基编辑引入m.1227G>A,在10-40%异质性下观察到呼吸链亚基NDUFB8和MT-CO2蛋白表达下降、基础耗氧率降低、代谢重编程(乳酸、TCA中间产物变化)以及核糖体组装异常。蛋白质组显示呼吸链复合物亚基广泛下调。

Fig. 4 · 原文图,按文末许可署名使用
Fig. 4 · 原文图,按文末许可署名使用

图 4 Fig. 4。图4为单细胞多组学分析:a为工作流程;b为三个文库的异质性分布;c显示MT-CO1表达随异质性增加而下降;d-e为基因-异质性相关性排序;f为GSEA火山图,显示氧化磷酸化和核糖体相关基因下调。

亮点 4|单细胞多组学揭示功能显性阈值效应

单细胞多组学分析显示,MT-CO1表达随m.1227G>A异质性增加而进行性下降,即使异质性低于50%也明显。GSEA发现高异质性细胞中氧化磷酸化和核糖体相关基因下调。阈值分析表明表型在20-30%异质性时发生最大变化,支持功能显性。

生信可带走

这一块是给做分析的人用的,不是科普点缀。

  • 方法栈: WGS、热点检测算法、DdCBE碱基编辑、单细胞多组学、代谢组、蛋白质组
  • 公开数据: 原文未给出公开组学登录号
  • 代码: 原文未给出公开 GitHub/GitLab 仓库
  • 谁该点开原文: 研究线粒体遗传学、肿瘤基因组学和代谢重编程的生信与实验人员。

带走一句

分析mtDNA突变时,不要默认功能隐性;rRNA热点在低异质性下即可产生显性表型,应纳入肿瘤进化模型。

本页供学习交流,不替代原文,也不代表原作者、期刊或资助方立场。引用科学结论请以正式发表版本为准。 图表按 CC BY(可转载,需署名) 使用。 版权声明

微信二维码

需要同类分析?扫码加微信,一对一沟通需求,交付后终身售后。

添加微信

需要同类分析?扫码加微信,一对一沟通需求,交付后终身售后。