精准医学杂志投稿

发表期刊：Nat Commun 发表日期：2022 Mar 28 DOI:          最近的研究表明，包括 AML 在内的儿童和成人癌症的基因组图谱存在显着差异。 与成人相比，儿童 AML 在 MYC ITD 和 WT1 中往往表现出更高的突变频率，而在 DNMT3A 和 TP53 中的突变频率较低。 同时有研究表明，不同的种族背景可能对疾病发展和进展的分子驱动因素产生深远影响。根据这一证据，值得注意的是，尽管综合基因组研究揭示了儿童 AML 的关键基因组异常，但这些观察结果主要基于西方人群患者的基因组分析，中国 AML 患者仍然缺乏基因组分析。 （1）骨髓样本取自 2001-2018 年在上海儿童医学中心 (SCMC) 血液肿瘤科诊断为 AML 的 292 名患者。 （2）收集了10个具有匹配RNA seq和WGS数据的所有诊断样本 (1)转录组测序（RNA-seq）和分析、SNV与Indel分析 (2)驱动突变分析：PeCanPIE24和MutClan分析用于识别驱动突变；从 Cloud收集的983例儿童肿瘤中，共有6975733个已发表的体细胞突变用于构建突变簇背景，如果突变在PeCanPIE分析中被归类为gold或显著位于突变簇中，则该突变将被视为潜在的驱动因素 (3)RNA-seq检测驱动突变：10个具有匹配RNA seq和WGS数据的所有诊断样本；对RNA-seq数据进行SNV和INDEL分析；收集WGS检测到的SNV和INDEL作为基准。PeCanPIE应用于所有突变，本分析仅包括分类为金或银的突变。 (4)统计分析：相关性分析、生存分析、Cox回归分析         通过研究2001-2018年在上海儿童医学中心（SCMC）诊断和治疗的292名中国儿童AML患者，分析了基因组改变。中国人群的临床特征与西方人群的临床特征相当。然而，中国AML患者较年轻，中位年龄为岁，而西方队列为岁。患者接受AML-SCMC-2009-A和AML-SCMC-2009-B方案的预后无显著差异或在一段时间内接受治疗的患者之间无显著差异（图S1）。转录组测序（RNA-seq）应用于所有肿瘤样本，并分析序列突变和基因重排。         RNA-seq数据分析显示，292例患者中有200例（）发生224次重排，涉及97个基因（图1a）。与之前的报道一致，在中国AML患者中检测到的最常见融合包括RUNX1-RUNX1T1、KMT2A重排和NUP98重排。此外，在本研究队列的三名患者中发现了涉及XPO1基因的复发性帧内融合（图1b），包括XPO1-TNRC18融合和XPO1-MLLT10融合。值得注意的是，在另外两名最近在SCMC诊断的AML患者中也检测到XPO1-TNRC18融合。有趣的是，所有四个携带XPO1-TNRC18的病例都属于M7组，并且没有任何已知的AML驱动融合，支持这种融合是这些病例中AML的潜在驱动因素，并且在FAB组中定义了以前未分类的分子亚型。总体而言，XPO1-TNRC18病例占AML M7亚型的（队列中40例中有2例）。在目前的分析中，还发现了以前未观察到的其他基因重排，包括PTPRA-FUS、ZEB2-ATIC和MSI2-UBE3C。         作者优化了变体调用和处理管道，以分析来自纯肿瘤RNA序列数据的潜在体细胞和癌症相关序列突变。首先为了评估这种方法的性能，将其应用于从之前发表的10例儿科ALL病例中收集的RNA-seq数据，还对每个病例的肿瘤和缓解样本的全基因组测序（WGS）数据进行了匹配。结果表明，对纯肿瘤RNA序列数据的分析成功地识别了WGS发现的的驱动突变（21例中的18例）。同时，RNA-seq分析检测到另外9种潜在的驱动突变，包括NRAS G13D和KRAS G13D等。在这9个突变中，有8个在之前的研究中使用基因组DNA被纳入捕获验证实验。在这8个突变中，有7个验证成功。值得注意的是，所有9个突变均为亚克隆，由于覆盖范围不足而在WGS分析中被遗漏。这些结果证明了从RNA-seq分析序列突变的能力，尤其是用于检测亚克隆变异。使用这种方法，作者接下来确定了总共975个影响305个基因的非同义序列突变。这些突变包括707个单核苷酸变异（SNV）和268个插入/缺失（indels），每个病例检测到的中位数为4个突变。通过应用组合策略进一步分析潜在的驱动突变，将PECANPI24的突变致病性分析与MutClan的突变聚类分析相结合，共鉴定出572个潜在的驱动因素变异，影响73个基因。此外，CICERO25检测到24个内部串联重复（ITD），影响FLT3和MYC。总的来说，（n=239）的病例检测到了驱动序列突变，每个病例检测到三个驱动突变的中位数。        在之后的分析中只关注驱动突变。在这些驱动因素中，超过5%的患者有10个基因反复突变（图2a），包括FLT3、KIT、NRAS、KRAS、CEBPA、ASXL2、PTPN11、CSF3R、GATA2和JAK2。值得注意的是，作者发现了儿童AML的驱动基因，包括LZTR1和SPOP，之前未发现与儿童AML相关联，以及ARID2 和 SH2B3 的功能缺失突变，据报道它们在其他儿科癌症（如 ALL）中具有致病性，但在 AML 中没有。尽管这些基因中的大多数基因的基因组突变先前已经在儿童AML中报道过，但与代表西方人群的TARGET AML研究相比，中国队列中的突变发生率有显著不同。在SCMC或TARGET队列中>4%的患者中，驱动基因内的21个基因或热点反复突变，其中10个（）显示出显著不同的突变频率（图2b、c）。其中4例在中国患者中显示出较高的突变频率，包括ASXL2、JAK2、CSF3R胞质结构域和KIT外显子17（KIT-E17)。另一方面，中国患者的FLT3、FLT3 ITD、NRAS、WT1、NPM1和TET2突变频率较低。值得注意的是，中国和西方患者之间的不同突变频率主要由3-14岁的患者构成（图S4a）。此外，作者发现FLT3和NRAS突变在中国队列的年轻患者中更常见（图S4b）。         作者进一步整合了不同的变异类型以分析中国儿童AML的基因组图谱。除了上述融合和序列突变外，用RNA-seq分析发现CBL外显子8/9缺失。本研究队列中有12例（）发现了这种局部缺失，与TARGET队列相当。总的来说，在93个基因中发现的驱动突变被分为六条途径（图3a）。AML中最常见的突变途径是转录调节、表观遗传学和RAS信号传导。共有 的患者检测到激活其他信号通路的突变，包括JAK-STAT等。值得注意的是，在中国患者中，RAS信号通路的突变频率显著降低，这与在RAS途径中观察到的单个基因的低突变频率一致。         接下来，作者研究了在中国AML患者中检测到的驱动基因突变的配对关系。共发现150对基因、突变热点或FAB组（排除融合伙伴或单个基因内不同结构域之间的配对后为143对）显著同时或完全突变（图3b）。通过这项分析，在建AML中立了的多个关联，包括CSF3R和KIT、ASXL2和KIT、DHX15和KIT、DNM2和JAK3等的共突变。对这些以前未曾描述过的关联的观察，只能部分解释为这些基因在中国 AML 患者中的突变频率较高，因为在不同人群中也发现了完全不同的突变关联。         生存分析显示驱动基因组畸变与患者预后相关（图4a）。CBFB-MYH11融合与预后良好相关，而NUP98-KDM5A/NSD1、FUS-ERG和   CBFA2T2-GLIS2与预后不良相关。另一方面，作者注意到RUNX1-RUNX1T1融合和KMT2A重排患者的5年无事件生存率（EFS）相似。关于这些突变，发现CEBPA、NPM1和GATA2的突变与良好的预后相关，而RUNX1和FLT3、ITD的突变与较差的预后相关，这与之前在西方队列中的发现一致。携带上述驱动基因畸变的患者在中国和西方队列中的预后情况相似（图S6）。         除了这些先前建立的关联，还发现携带UBTF突变的患者与野生型相比预后更差（图S7a）。与单独携带FLT3变体的患者相比，携带FLT3变体并伴有UBTF突变、RUNX1突变或NUP98重排的患者的预后更差（图S7b-d）。另一方面，FLT3 ITD和NPM1突变的患者在当前的分析中显示出良好的结果（图S7e）。         研究中应用了单变量和多变量Cox回归分析。单变量Cox回归显示上述结果一致，而多变量Cox回归显示CBFA2T3-GLIS2、FUS-ERG、NUP98重排、FLT3 ITD和RUNX1突变与不良预后独立相关，而GATA2与良好预后独立相关。由于UBTF和CEBPA突变分别与FLT3 ITD和GATA2存在显著的共突变，因此不包括在多变量Cox回归模型中。此外，将上述五个与不良预后相关的基因组因素合并为高危基因型，以及CR1状态和GATA2，再次进行多变量Cox回归分析。结果显示，合并高危基因型是与不良预后显著相关的独立危险因素。         在本研究队列中，有78名患者至少有一项上述改变与良好或不良结果相关。接下来，作者检查了剩下的患者是否有任何其他与临床相关的改变，将重点放在那些缺乏任何畸变的患者身上，这些畸变与预后有着良好的关联。除上述78名患者外，还有24名患者因携带TP53或ASXL1突变、DEK-NUP214融合或确认的染色体异常而被进一步排除在外。发现，对于剩下的患者，第一周期诱导治疗后的治疗反应是与患者预后最显著相关的因素之一（图S8a）。在一个诱导周期（CR1）后未达到完全缓解的患者显示出不良结局，与携带与不良结局相关的基因组变异的患者相似（图S8b）。另一方面，虽然CR1患者的预后相对较好，但与携带与良好预后相关的变异的患者相比，这些患者的预后较差（图S8c），表明该CR1组中存在混合患者。         接下来，分析了129例CR1患者的每个融合亚型的基因组畸变与预后之间的关系，包括RUNX1-RUNX1T1，KMT2A重排，其他融合病例，融合阴性病例。发现CSF3R和KIT-E17的突变在RUNX1-RUNX1T1患者中与不良预后相关。CSF3R和KIT外显子17突变在RUNX1-RUNX1T1融合亚型中显著共同发生。事实上，与携带CSF3R和KIT外显子17突变的RUNX1-RUNX1T1融合亚型CR1患者相比，携带CSF3R和KIT外显子17突变的CR1患者的预后更差(图4b）。         接下来，作者根据本研究建立的临床相关性修改了欧洲白血病（ELN）遗传风险分类模型。提出的SCMC pAML模型（图5a）的特点是调整了几种基因畸变的风险分类。例如，FUS-ERG、CBFA2T3-GLIS2、NUP98-KDM5A和NUP98-NSD1被归为高风险组（HR），而KMT2A重排被归为中等风险组（IR）。重要的是，RUNX1-RUNX1T1融合的患者在之前的模型中与良好的预后相关，根据CR1状态和随后获得的突变（包括CSF3R和KIT-E17突变）进一步细分。与ELN模型相比，SCMC pAML模型识别出更多的HR患者和更少的低风险（LR）患者，以及相似数量的IR患者（图5b）。尽管两种模型分类的风险组在患者预后方面存在显著差异，但使用SCMC pAML分类的LR组和IR组显示出更有利的结果（图5c），5年EFS率分别为和，这些比率显著高采用ELN模型分层的LR组和IR组。另一方面，与ELNHR组相比，SCMC pAML HR组的预后更差，5年EFS发生率分别为和，无统计学差异。此外，在多变量cox模型中，将SCMC pAML模型分为HR组是一个独立的风险因素，SCMC pAML HR患者的不良事件风险显著增加。         研究展示了中国儿科AML驱动器改变的基因组景观，并发现了以前未描述的基因组畸变，包括 XPO1-TNRC18 融合。全面比较中国和西方AML队列，发现基因组改变特征明显。例如，中国AML患者在 KIT 和 CSF3R 中表现出突变，而RAS信号通路中的基因突变较少。同时提出了一种精细的预后风险分类模型，该模型更好地反映了中国AML患者的不良事件风险。本研究结果揭示了一个临床相关的突变谱，该谱在突变频率和突变共发生模式方面与西方队列不同。这些发现进一步阐明了儿童急性髓系白血病的复杂性，并强调了在精确医学时代为临床管理建立风险分层时考虑种族背景的重要性。