《Acta Neuropathologica》：The association between MGMT expression and outcome in H3 G34-mutant diffuse hemispheric glioma is robust to adjustment for pre-analytic and technical variables
编辑推荐：
研究人员感谢Shen和Sheng的仔细再分析，并回应其提出的技术批量效应（batch effects）可能混淆H3 G34突变型弥漫性大脑半球胶质瘤（diffuse hemispheric glioma，DHG）中MGMT表达与预后关联的观点。研究人员分析了1
研究人员感谢Shen和Sheng的仔细再分析，并回应其提出的技术批量效应（batch effects）可能混淆H3 G34突变型弥漫性大脑半球胶质瘤（diffuse hemispheric glioma，DHG）中MGMT表达与预后关联的观点。研究人员分析了110例H3 G34突变型DHG的RNA测序（RNA-seq）数据，所有样本均在St. Jude中心使用NovaSeq 6000平台测序，并纳入原研究报告中的全部RNA测序肿瘤。原分析中4例SJC-DS-1001样本因在Shen和Sheng主成分分析（principal component analysis，PCA）中最极端且来自较旧平台，已被排除。计数采用TMM标准化（edgeR），转换为log2每百万计数（counts per million，CPM），并使用limma的removeBatchEffect校正处理批次。研究人员考察了4个前分析和技术变量：收集年份、保存方法（福尔马林固定石蜡包埋（formalinfixed paraffinembedded，FFPE）或新鲜/冷冻）、数据发布年份、肿瘤含量（由H3 G34变异等位基因频率和H&E评估）。原始计数中前两个主成分随文库大小分布，并按发布年份分离；校正后发布年份不再分离。4例低复杂度FFPE样本和1例新鲜/冷冻样本为离群值，排除后结果不变。Shen和Sheng报道MGMT与CD34（r=0.76）、ACTB（r=0.75）在原始特征计数中强相关；研究人员队列中原始计数相关较弱（r=0.29和0.42），校正后消失（r=-0.04和0.01）。原始ACTB计数在新鲜/冷冻中比FFPE高4.6倍，归一化后消除。4个变量在单因素和多因素模型中均与MGMT表达、无进展生存（progressionfree survival，PFS）或总生存（overall survival，OS）无显著关联；MGMT表达按保存方法无差异（p=0.71）。在25例原报告患者中，连续MGMT表达校正4个变量后仍与PFS（每翻倍校正风险比（hazard ratio，HR）1.64，95% CI 1.20–2.24，p=0.002）和OS（HR 3.00，95% CI 1.73–5.21，p<0.001）相关。Shen和Sheng认为共享伪影可能夸大MGMT表达与基因体/内含子CpG甲基化相关；由于甲基化来自独立DNA芯片平台，测序深度伪影不能单独解释。34例双数据类型肿瘤中，MGMT表达与cg15441483、cg18453665、cg15782766甲基化Spearman ρ为0.74、0.75、0.71，校正4变量后为0.58、0.78、0.75（均p<0.001）；MGMT-STP27启动子CpG（cg12434587、cg12981137）与表达负相关（ρ=-0.48、-0.36）。研究人员结论：原结论成立，回顾性多中心RNA-seq队列的批量效应可被缓解；建议记录前分析变量、归一化并校正批次、排除低复杂度文库、用独立检测验证；MGMT表达仍需更大前瞻队列验证方可指导治疗。
研究背景方面，H3 G34突变型弥漫性大脑半球胶质瘤（diffuse hemispheric glioma，DHG）属于罕见儿童脑肿瘤，相关队列通常回顾性、多机构、跨年度并混合保存方法，容易产生高通量数据中的批量效应（batch effects）。Shen和Sheng提出，基于bulk RNA测序（RNA-seq）的MGMT表达谱可能受技术批次效应混淆，并质疑MGMT表达与基因体/内含子CpG甲基化相关性是否由共享伪影导致。为回应这一方法学争议，研究人员在《Acta Neuropathologica》发表通讯/再分析类论文，检验MGMT表达与结局的关联是否在前分析变量和技术变量校正后仍然稳健。
主要关键技术方法方面，研究人员使用110例H3 G34突变型DHG样本，均来自St. Jude中心NovaSeq 6000平台测序；排除旧平台SJC-DS-1001极端样本；采用edgeR进行TMM标准化，转换为log2每百万计数（counts per million，CPM），用limma removeBatchEffect校正处理批次；评估收集年份、保存方法（福尔马林固定石蜡包埋（FFPE）或新鲜/冷冻）、数据发布年份、肿瘤含量四类前分析/技术变量；对25例原队列做Cox回归，对34例有甲基化数据的肿瘤做Spearman相关与偏相关分析。
研究结果方面，Pre-analytic and technical variables and MGMT expression in 110 H3 G34-mutant diffuse hemispheric gliomas：通过PCA与相关性分析发现，原始计数中样本按文库大小和发布年份分离，MGMT与CD34、ACTB相关性来自测序深度缩放；TMM归一化加批次校正后，发布年份不再分离，MGMT与看家/内皮标志基因相关性消失，低复杂度FFPE样本和个别新鲜/冷冻样本为离群值但未进入原模型，排除后结论不变。
A confounder must be associated with both MGMT expression and outcome：通过单因素与多因素线性/Cox回归发现，收集年份、保存方法、发布年份、肿瘤含量均不与MGMT表达、无进展生存（PFS）或总生存（OS）显著相关，MGMT表达按保存方法无差异，说明这些变量不构成混淆因子。
continuous MGMT expression adjusted for four variables remained associated with PFS and OS：在25例原报告患者中，连续MGMT表达校正四个变量后，仍与PFS（每翻倍HR 1.64，95% CI 1.20–2.24，p=0.002）和OS（HR 3.00，95% CI 1.73–5.21，p<0.001）显著相关，证明MGMT表达与结局关联非由所检技术/前分析变量驱动。
methylation measured from DNA on a separate array platform：在34例同时有RNA-seq与DNA甲基化阵列数据的肿瘤中，MGMT表达与cg15441483、cg18453665、cg15782776等基因体/内含子CpG甲基化强相关，校正四变量后仍强；MGMT-STP27启动子CpG（cg12434587、cg12981137）与表达负相关，模式跨癌种一致且有胶质母细胞瘤细胞去甲基化实验支持，说明DNA甲基化相关不是RNA-seq深度伪影单独造成。
讨论与结论翻译方面，研究人员指出结果支持原结论，并表明回顾性、多机构RNA-seq队列中的批量效应可被缓解。研究人员同意Shen和Sheng所言此类队列需谨慎，并提出：记录前分析变量；在样本比较前归一化并校正批次；排除低复杂度文库，因为校正无法挽救此类数据；检验对这些变量调整后标志物与结局关联是否改变，并用独立检测确认关键发现。研究人员亦同意单纯相关不能证明基因体甲基化调控MGMT，且MGMT表达在指导治疗前仍须在更大、前瞻性队列中验证。结论可译为：这些结果支持原结论，并显示回顾性多机构RNA-seq队列中的批量效应可以被减轻；此类队列需要谨慎处理，研究应记录前分析变量、归一化并校正批次、排除低复杂度文库、检验变量调整后标志物—结局关联是否改变，并以独立检测确认关键发现；MGMT表达与DNA甲基化的相关性不能单独证明因果，MGMT表达作为治疗标志物仍待前瞻性大队列验证。
