乘客共缺失混淆了以PTEN缺失为锚点的谷氨酰胺分解特征：位置感知特征设计的警示案例
核心概要
该研究利用TCGA泛癌图谱的GISTIC拷贝数将肿瘤分为PTEN完整、半合子缺失和纯合缺失三档，在十四种肿瘤类型中评估一个五基因谷氨酰胺分解特征（GLS、SLC1A5、GOT1、GLUD1、GPT2），发现该特征随PTEN缺失而下降（十四种全部下降、十二种显著），但这一下降并非由MYC介导，而是由染色体位置决定：GLUD1与GOT1位于10q上与PTEN相邻，三十七个无谷氨酰胺分解注释的邻近基因与PTEN拷贝数的相关性与这两个基因相当（平均rho 0.843对0.842），共缺失保真度随与PTEN距离单调下降（rho = -0.993），校正全基因组非整倍体后关联不变，而以各基因自身拷贝数作条件则使
Fig 2. Glutaminolysis signature by PTEN dosage in the six tumor types with the largest decline
bioRxiv · 第 8 页深度剖析
在十四种肿瘤类型中，谷氨酰胺分解特征随PTEN拷贝数缺失加重而下降，十二种经Benjamini-Hochberg校正后显著，最强效应见于前列腺癌（Spearman rho = -0.54，校正P = 5.3 x 10^-37），且没有任何肿瘤类型呈现信号模型所预测的正向趋势。 此前将PTEN缺失与谷氨酰胺分解联系起来的证据主要来自单个细胞系模型和前列腺癌，且从未在多肿瘤类型中、以基因剂量为关注点检验谷氨酰胺分解转录程序是否真的随PTEN拷贝数缺失而变化。 基于TCGA泛癌图谱十四种肿瘤类型、以预先设定的纳入规则（每个剂量档至少十例肿瘤）选定的队列，采用两-tailed Spearman相关并作多重校正，效应量以Cliff's delta及2,000次重抽样的bootstrap置信区间报告；效应量差异很大，前列腺解释29.1%方差，而乳腺和胃仅0.5%和0.4%。
该关联由染色体位置而非通路归属决定：GLUD1和GOT1位于10q23-q24，与PTEN共缺失，三十七个无谷氨酰胺分解注释的邻近基因与PTEN拷贝数的平均相关性为0.843，与GLUD1和GOT1的0.842相当，而位于其他染色体的三个特征基因平均仅0.034；共缺失保真度随与PTEN转录起始位点距离单调下降（rho = -0.993）。 此前没有研究把特征基因的染色体位置作为直接判别标准，来区分调控性关联与结构性乘客共缺失。 以Ensembl（GRCh38）位置对10q23-q24区间三十七个无注释基因逐一计算其拷贝数与PTEN拷贝数的相关性，并因邻近基因共享拷贝数谱而不报告名义P值；共缺失事件为局灶性，纯合缺失肿瘤中共缺失序列中位跨度0.72 Mb，90%局限于PTEN 2 Mb以内。
两种校正给出不同答案：校正全基因组非整倍体后10号染色体关联基本不变（如GLUD1在前列腺的偏相关rho -0.59对未校正-0.59），而以各基因自身拷贝数作条件则使10号染色体基因的平均系数从-0.321降至-0.022、关联消失，非10号染色体基因则始终不受影响。 全基因组非整倍体评分无法分离单条染色体臂的丢失，因此这一常用校正在此处会给出误导性的阴性结果，而基因水平拷贝数校正才是直接检验。 以Spearman偏相关分别校正每个样本的非整倍体评分、10q臂水平代理（距PTEN超过2 Mb的区间基因拷贝数中位数）以及各基因自身拷贝数，三种校正并列比较。
基于突变的自然实验把PTEN功能丧失与染色体片段丢失分开：合并所有队列时，PTEN点突变且拷贝中性的肿瘤中GLUD1低于拷贝中性野生型（435对3,561例，差异-0.24，P = 1.0 x 10^-8），但该效应由子宫内膜癌驱动，排除该谱系后无可检测差异（P = 0.22，等效性P = 0.015），而半合子缺失（均值-0.44）和纯合缺失（均值-0.89）在每次分析中都剂量依赖性地降低GLUD1；DepMap中GLS基因效应评分不随PTEN剂量分层变化（Kruskal-Wallis P = 0.74）。 乘客基因应对缺失有反应、而对单纯突变无反应，这一自然实验此前未被用于检验PTEN-谷氨酰胺分解关联。 采用两-tailed Mann-Whitney U检验，按肿瘤类型分层并在队列内z标准化后合并，等效性以预先设定的0.25个队列内标准差为界用双单侧检验评估；子宫内膜癌的例外在微卫星稳定肿瘤中效应更大（差异-0.53，P = 1.2 x 10^-6），排除了超突变作为原因。
启示与展望
该结果适用于以拷贝数定义PTEN剂量的实体瘤转录组分析，其核心结论建立在肿瘤DNA中测得的基因水平GISTIC拷贝数之上，因此不受肿瘤纯度与基质稀释影响；所提出的四步筛查流程被作者定位为筛查工作流而非经过验证的工具，前三步只需已有数据，第四步需要独立的CRISPR依赖等功能数据，且依赖本身具有情境依赖性。对希望测量PTEN缺失背景下谷氨酰胺分解转录调控的研究，作者建议使用排除10q23-q24乘客基因GLUD1和GOT1、仅依赖非10号染色体成员（GLS、SLC1A5、GPT2）的染色体感知特征，并可补充GLS2或SLC38A1等非10号染色体谷氨酰胺分解基因。
以拷贝数定义的剂量无法捕捉点突变、启动子甲基化或结构变异，因此部分功能上PTEN缺陷的肿瘤会被归为完整，使结果偏向零假设；表观遗传沉默可带来不同于拷贝数丢失的代谢后果。阈值化GISTIC判定可能在广泛低水平丢失中把样本归为纯合，而低肿瘤纯度则反向掩盖真正的纯合缺失，两者都是剂量分档的误分类来源。胶质母细胞瘤仅有十三例PTEN完整肿瘤，基线不稳定。GLS在多数细胞系中本就是弱依赖，因此DepMap结果只是限定而非证明不存在功能依赖。批量RNA-seq受肿瘤纯度与基质稀释影响，GLS和SLC1A5在免疫细胞中表达，因此非10号染色体亚特征在表达层面的解读需谨慎。前列腺癌中残留的负向关联基础尚未确定，且该分析是在泛癌结果已知后进行的探索性分析。子宫内膜癌中突变单独降低GLUD1的机制未明。该测量混淆是否出现在其他易缺失位点，本文数据无法回答。