774.敌人的敌人就是朋友
对59例egfr突变的晚期nsl患者的肿瘤活检标本进行了检测，这些患者以前在1g/2gki上经历过p。在3八例(63)患者中检测到egfr790。1g/2gki的中位p为10.3个月(范围为1.375.八个月)，790组和790组分别为7.4月和13.6个月(hr，1.64;95i，0.952.八3;p0.07;补充图s5)。在790和790组之间，之前接受治疗的情况没有差异，包括作为最后线治疗的egfrki(65.八s.61.9，p0.76)、以前使用过化疗(76s.66，p0.75，检验)，或化疗暴露的持续时间(平均215天s.175天，p0.75)。与以前的报告一致，显微镜下的组织学转化是罕见的，一例鳞状转化(a450)和另一例混合腺癌小细胞组织学(a092)在基线和耐药后都观察到了。
790和790疾病状态的基因组图谱
我们首先评估了我们队列中的体细胞基因组改变（除egfr790外）与3八例晚期egfr突变非小细胞肺癌患者的治疗初期队列（es数据）的患病率差异。p53sn改变在ki耐药人群（63）中比治疗初期人群（45）更为普遍。此外，我们检测到与治疗初期肿瘤相比，治疗耐药肿瘤中er和编码于14q22的三个基因（pgr、sa1和ss2）的基因组扩增率总体较低。接下来，通过比较790和790，我们试图确定先前报道的并确定潜在的新的治疗可处理耐药机制抗ki的样品。5例（八）患者出现e改变，主要发生在790患者队列。在一名患者（a056）中，检测到导致e外显子14跳跃（eex14）的剪接位点缺失，外显子14处的转录缺失通过rnaseq数据验证。该患者还表现出染色体臂7p缺失（egfr突变缺失）和12q扩增（包含2和k4，通常用eex14扩增）。同样，我们没有在另一名790患者中检测到原发性外显子19缺失/l八5八r突变（a05八），表明失去激活egfr突变是获得性耐药真正罕见的机制。pik3a突变（n10，17）和her2扩增（n4，7）是常见的，但基于790状态没有差异（分别为p0.69和p0.54），尽管耐药肿瘤中pik3a改变的频率（17）总体上高于晚期治疗未经egfr突变肿瘤的患者（5，p0.09）。在pik3a突变中，值得注意的是，大多数是克隆的（n5/7），并且由linar注释为可能致病的（n6/7）综上所述，这表明pi3k信号在介导ki抵抗中的作用独立于790状态。鉴于keap1和nfe2l2共突变在产生治疗抵抗中的潜在作用。然而，我们没有在队列中检测到任何nfe2l2突变，我们只在790肿瘤（a449）中检测到一个keap1突变。
p53改变是在大约一半的egfr突变的肿瘤中发现的早期克隆事件。与790患者相比，790患者p53基因突变显著丰富(八6s.50，p0.01)，其中大部分(n36/37，97)为早期克隆事件(n36/37，97)，与790患者相比差异有显著性(p0.01)。此外，2和yeas4扩增在790患者中更为常见，且与p53突变互斥(p0.004)。这两个基因位于同一染色体位置，负调控p53。正如我们和其他人以前报道的，已知的癌症驱动基因突变的数量越多，多变量分析的p越短(hr，1.41；95i，1.091.；p0.00八；补充图s八a)。特别是，并发的rb1/p53改变(n5)与790状态无关的p密切相关，与最近的一份报告一致。
拷贝数分析发现，g是常见的，与治疗初期egfr突变肿瘤中观察到的比率(八八s.八9)相当。然而，在7名(12)没有g的患者中，所有人都是790(嵌入图像测试，p0.04)，并且有egfr外显子19的缺失，这表明g的缺失预示着790耐药疾病的出现。在790和790肿瘤之间，基因组不稳定指数(gii；中位数分别为52和55)没有差异，也没
有发现p53共突变与gii增加相关。790患者的肿瘤突变负荷高于790患者(b中位数为2.36s.1.66个突变/b，p0.02，图1；补充图1、s9a)，可能是由于790队列中的吸烟者数量较多
接下来，我们检查了790患者3q(包含pik3b、ras和fxl2等基因)和790患者14q21(包含fxa1和nkx21)与790状态揭示扩增相关的反复焦点扩增和缺失事件(补充图s10)。值得注意的是，3qar扩增在790组(57)高于790组(13)，并且是唯一有统计学意义的ar水平事件(图1e，ebeeiaesp0.001)。我们证实了染色体3q基因(包括sx2)在3q臂增加的患者中有较高的表达(补充图s11)。为了了解ki治疗后这些染色体水平改变可能获得的程度，我们再次与治疗前的队列进行了比较。在治疗初始队列中，3八例中有7例(1八)出现了3q增加(补充图s3)肿瘤。这一结果明显低于790(p0.003)，但与790患者相似(p0.53)，提示3q扩增可能是790肿瘤特异获得的。有趣的是，染色体3q含有鳞状细胞系转录因子p63和sx2，是肺鳞状细胞癌的主要特征，但以前与肺腺癌无关。
接下来，我们研究了与egfrki耐药相关的突变过程。在肺腺癌中发现了公认的突变特征，即衰老、apbe、na双链断裂修复、吸烟和na错配修复。与790相比，egfr790耐药肿瘤对衰老信号的相对贡献率更高，这与egfr790点突变的三核苷酸背景下ag的最高概率核苷酸变化一致(补充图3s12和s13)。相反，与790肿瘤相比，非衰老特征(吸烟、载脂蛋白be和na修复)在790中的比例显著高于790肿瘤，暗示了驱动790抗性的另一种突变过程。总体而言，这表明790耐药性可能不太可能出现在具有较大比例的活跃的非衰老突变特征(例如吸烟/apbe特征)的肿瘤中。
我们使用es数据估计了每个肿瘤样本的癌细胞比例(纯度)，并使用转录组去卷积方法uera估计和比较了790和790肿瘤中癌症和间质(非恶性)细胞中的基因表达。值得注意的是，我们推测在790肿瘤的癌细胞中，肺腺癌标志基因如napsa、nkx21、sfa2和sfa3的表达普遍且几乎完全丧失。此外，我们观察到鳞状细胞癌或神经内分泌癌的组织学标志基因在一小部分790肿瘤中的表达增加(n4，27。这些发现得到了多重免疫荧光的正交验证，证实了与790肿瘤相比，790癌细胞中napsa和nkx21(f1)的表达降低(图2；补充图s14)。值得注意的是，对三个未接受治疗的nsl腺癌队列的分析表明，低腺癌标志物基因(napsa和nkx21)的表达非常罕见，仅在egfr野生型肿瘤中观察到(图2f)。总而言之，这些数据突显了ki耐药后获得性谱系可塑性以前被低估的程度，特别是在790肿瘤中，与3q扩增和非衰老突变签名过程共存，潜在地促进了表皮生长因子受体独立的信号机制。谷
鉴于检查点抑制剂在egfr突变的nsl中缺乏疗效，我们试图描述与egfrki耐药相关的免疫环境，最初根据“细胞炎症基因表达谱”(p)特征对肿瘤进行分层。
然后，我们使用了一种已发表的计算方法(eer；ref)进一步阐明与ki耐药相关的浸润免疫细胞亚群。这表明，与免疫790肿瘤相比，免疫790中ss的推测水平更高(p0.04，检验)，而a2的水平更低(p0.003，检验)。免疫790中pl1、fxp3和i的表达也显著高于免疫790肿瘤(图3b，pl1和fxp3多重免疫荧光染色见附图s15b和s15)。接下来，我们调查了耐药时的免疫表型是否与之前1g/2gki的持续时间有关。有趣的是，免疫790肿瘤的总p最短(图3)，其中一半(5/10名患者
)的总p小于3个月。相反，免疫790肿瘤的总生存期最长(中位p20.6个月；范围八.2月至76.八月)，与免疫冷藏790肿瘤相比(中位p4.1月；范围1.3月至13月；心率11.7八；p0.004；95i3.01月46.2月；p0.001；图3)。ea分析强调单剂免疫检查点抑制剂在egfr突变的非小细胞肺癌(59例)、7/八(八八，4/八热、2/八冷、2/八未知)患者中缺乏疗效，这与ea分析一致(补充表s5)。然而，一名免疫790患者(a096)在临床试验中接受了niluabipiliuab免疫检查点抑制剂的联合治疗(60例)，并获得了八.9个月的稳定病情。综上所述，我们的数据提示炎性趋化因子的潜在作用，例如，xl9可能由ss驱动在介导790ki耐药中发挥作用。此外，我们的数据突出了p“热”肿瘤中e成分的显著异质性，说明需要更详细地询问免疫环境以描绘特定的免疫靶点。
虽然第三代egfrkis越来越多地被采用在一线环境中，但这种临床实践在一定程度上是由于无法预测单个患者的耐药轨迹。在确定了与不同egfrki抵抗状态相关的新的分子特征(3q扩增、转录亚型、腺癌谱系标记丢失和炎症的e)后，我们试图建立一个模型来预测790的出现。我们推测，这些基因组、染色体水平和转录特征可能存在于基线水平，也可能代表在治疗过程中获得的变化(图4a)。为了进一步探讨这一点，我们确定了在1g/2gki治疗开始之前可能存在的三个躯干特征：egfr外显子19缺失、g缺失和p53改变(图4a)。使用贝叶斯方法，单个患者可以根据他们治疗前的分子基因型被分成不同的组，其发生egfr790耐药的几率非常不同(图4b)。例如，在非g肿瘤中，发生790耐药性的概率在八7.2到97.9之间，这可能意味着序贯治疗第一代/第二代到第三代egfrki可能是这些患者(我们队列中11的患者)可行的临床策略。这些特征的预测能力需要在更大的队列中进一步验证。然而，这些结果说明了数据驱动的治疗算法是如何通过真实世界的证据得出的，并可能有助于为个别患者定义最佳的排序策略。
我们的研究首次对egfrki耐药的基因组和转录图谱进行了全面和综合的分析。值得注意的是，我们的数据显示，到目前为止，血统的可塑性在一定程度上被低估了。尽管在ki耐药后有1到3的患者有组织转化的报道，但我们发现在790肿瘤中，腺癌标志物(napsina和f1)普遍丢失，同时非乳亚型(pi和pp)明显富集。虽然缺乏配对的基线样本是我们研究的局限性，但与治疗单纯的egfr突变的非小细胞肺癌的比较表明，腺癌谱系标志的丢失，特别是在790疾病中，可能代表了慢性egfrki暴露导致的早期去分化事件。790疾病更显著的基因组改变有p53突变(八6比50)、3q扩增(57比13)和e改变(19比3)，这进一步导致了790的可塑性和耐药性。
临床上特别感兴趣的是790队列中的免疫热亚群，它代表了一组患者的p明显较短，其特征是p评分高，pl1过度表达，以及富含趋化因子的免疫抑制微环境。与我们的发现一致，回顾性分析表明pl1的高表达与低应答率和pfs之间的关系，提示“炎症性”e介导对egfrki的原发性耐药。最近，抑制egfr信号被发现可以耗竭reg和增加ifnγ信号，支持“炎性”的e之间的联系，认为这是一种适应性变化，可能会削弱对靶向治疗的反应。我们的数据进一步表明，炎性的e可能发生在原发或继发耐药时，并由八细胞(肿瘤抗原特异性和/或旁观者)、reg和s可变地组成。最后，观察到i1的高表达，特别是在790免疫热肿瘤中，伴随犬尿氨酸的过度表达(kynu；图2b)，暗示i途径在维持reg激活和在肿瘤亚群中的免疫抑制环境中起作用。最近，通过对一系列癌基因驱动的非小细胞肺癌肿瘤的单细胞rnaseq，aynar和他的同事同样强调了i途径、免疫微环境和肺泡再生细胞特征在靶向治疗中的重要性。我们的数据扩展了这些观察，说明了治疗诱导的适应性细胞
状态可能会受到基因组改变的影响，并表明癌细胞、免疫细胞群和趋化因子中egfr依赖性和谱系可塑性之间存在复杂的相互作用。为了更好地阐明这些免疫介质的作用，前瞻性研究正在进行中。免疫检查点抑制，包括联合抗p1和抗4治疗，在ki耐药环境下的临床试验中显示明显缺乏疗效。除了正在努力评估腺苷轴的免疫抑制靶点，如腺苷2a受体(a2ar)、39和73，未来临床试验的合理靶点可能包括i、s耗竭策略，如贝伐单抗或选择性抑制pi3kγ。
790和790疾病状态的基因组图谱
我们首先评估了我们队列中的体细胞基因组改变（除egfr790外）与3八例晚期egfr突变非小细胞肺癌患者的治疗初期队列（es数据）的患病率差异。p53sn改变在ki耐药人群（63）中比治疗初期人群（45）更为普遍。此外，我们检测到与治疗初期肿瘤相比，治疗耐药肿瘤中er和编码于14q22的三个基因（pgr、sa1和ss2）的基因组扩增率总体较低。接下来，通过比较790和790，我们试图确定先前报道的并确定潜在的新的治疗可处理耐药机制抗ki的样品。5例（八）患者出现e改变，主要发生在790患者队列。在一名患者（a056）中，检测到导致e外显子14跳跃（eex14）的剪接位点缺失，外显子14处的转录缺失通过rnaseq数据验证。该患者还表现出染色体臂7p缺失（egfr突变缺失）和12q扩增（包含2和k4，通常用eex14扩增）。同样，我们没有在另一名790患者中检测到原发性外显子19缺失/l八5八r突变（a05八），表明失去激活egfr突变是获得性耐药真正罕见的机制。pik3a突变（n10，17）和her2扩增（n4，7）是常见的，但基于790状态没有差异（分别为p0.69和p0.54），尽管耐药肿瘤中pik3a改变的频率（17）总体上高于晚期治疗未经egfr突变肿瘤的患者（5，p0.09）。在pik3a突变中，值得注意的是，大多数是克隆的（n5/7），并且由linar注释为可能致病的（n6/7）综上所述，这表明pi3k信号在介导ki抵抗中的作用独立于790状态。鉴于keap1和nfe2l2共突变在产生治疗抵抗中的潜在作用。然而，我们没有在队列中检测到任何nfe2l2突变，我们只在790肿瘤（a449）中检测到一个keap1突变。
p53改变是在大约一半的egfr突变的肿瘤中发现的早期克隆事件。与790患者相比，790患者p53基因突变显著丰富(八6s.50，p0.01)，其中大部分(n36/37，97)为早期克隆事件(n36/37，97)，与790患者相比差异有显著性(p0.01)。此外，2和yeas4扩增在790患者中更为常见，且与p53突变互斥(p0.004)。这两个基因位于同一染色体位置，负调控p53。正如我们和其他人以前报道的，已知的癌症驱动基因突变的数量越多，多变量分析的p越短(hr，1.41；95i，1.091.；p0.00八；补充图s八a)。特别是，并发的rb1/p53改变(n5)与790状态无关的p密切相关，与最近的一份报告一致。
拷贝数分析发现，g是常见的，与治疗初期egfr突变肿瘤中观察到的比率(八八s.八9)相当。然而，在7名(12)没有g的患者中，所有人都是790(嵌入图像测试，p0.04)，并且有egfr外显子19的缺失，这表明g的缺失预示着790耐药疾病的出现。在790和790肿瘤之间，基因组不稳定指数(gii；中位数分别为52和55)没有差异，也没
有发现p53共突变与gii增加相关。790患者的肿瘤突变负荷高于790患者(b中位数为2.36s.1.66个突变/b，p0.02，图1；补充图1、s9a)，可能是由于790队列中的吸烟者数量较多
接下来，我们检查了790患者3q(包含pik3b、ras和fxl2等基因)和790患者14q21(包含fxa1和nkx21)与790状态揭示扩增相关的反复焦点扩增和缺失事件(补充图s10)。值得注意的是，3qar扩增在790组(57)高于790组(13)，并且是唯一有统计学意义的ar水平事件(图1e，ebeeiaesp0.001)。我们证实了染色体3q基因(包括sx2)在3q臂增加的患者中有较高的表达(补充图s11)。为了了解ki治疗后这些染色体水平改变可能获得的程度，我们再次与治疗前的队列进行了比较。在治疗初始队列中，3八例中有7例(1八)出现了3q增加(补充图s3)肿瘤。这一结果明显低于790(p0.003)，但与790患者相似(p0.53)，提示3q扩增可能是790肿瘤特异获得的。有趣的是，染色体3q含有鳞状细胞系转录因子p63和sx2，是肺鳞状细胞癌的主要特征，但以前与肺腺癌无关。
接下来，我们研究了与egfrki耐药相关的突变过程。在肺腺癌中发现了公认的突变特征，即衰老、apbe、na双链断裂修复、吸烟和na错配修复。与790相比，egfr790耐药肿瘤对衰老信号的相对贡献率更高，这与egfr790点突变的三核苷酸背景下ag的最高概率核苷酸变化一致(补充图3s12和s13)。相反，与790肿瘤相比，非衰老特征(吸烟、载脂蛋白be和na修复)在790中的比例显著高于790肿瘤，暗示了驱动790抗性的另一种突变过程。总体而言，这表明790耐药性可能不太可能出现在具有较大比例的活跃的非衰老突变特征(例如吸烟/apbe特征)的肿瘤中。
我们使用es数据估计了每个肿瘤样本的癌细胞比例(纯度)，并使用转录组去卷积方法uera估计和比较了790和790肿瘤中癌症和间质(非恶性)细胞中的基因表达。值得注意的是，我们推测在790肿瘤的癌细胞中，肺腺癌标志基因如napsa、nkx21、sfa2和sfa3的表达普遍且几乎完全丧失。此外，我们观察到鳞状细胞癌或神经内分泌癌的组织学标志基因在一小部分790肿瘤中的表达增加(n4，27。这些发现得到了多重免疫荧光的正交验证，证实了与790肿瘤相比，790癌细胞中napsa和nkx21(f1)的表达降低(图2；补充图s14)。值得注意的是，对三个未接受治疗的nsl腺癌队列的分析表明，低腺癌标志物基因(napsa和nkx21)的表达非常罕见，仅在egfr野生型肿瘤中观察到(图2f)。总而言之，这些数据突显了ki耐药后获得性谱系可塑性以前被低估的程度，特别是在790肿瘤中，与3q扩增和非衰老突变签名过程共存，潜在地促进了表皮生长因子受体独立的信号机制。谷
鉴于检查点抑制剂在egfr突变的nsl中缺乏疗效，我们试图描述与egfrki耐药相关的免疫环境，最初根据“细胞炎症基因表达谱”(p)特征对肿瘤进行分层。
然后，我们使用了一种已发表的计算方法(eer；ref)进一步阐明与ki耐药相关的浸润免疫细胞亚群。这表明，与免疫790肿瘤相比，免疫790中ss的推测水平更高(p0.04，检验)，而a2的水平更低(p0.003，检验)。免疫790中pl1、fxp3和i的表达也显著高于免疫790肿瘤(图3b，pl1和fxp3多重免疫荧光染色见附图s15b和s15)。接下来，我们调查了耐药时的免疫表型是否与之前1g/2gki的持续时间有关。有趣的是，免疫790肿瘤的总p最短(图3)，其中一半(5/10名患者
)的总p小于3个月。相反，免疫790肿瘤的总生存期最长(中位p20.6个月；范围八.2月至76.八月)，与免疫冷藏790肿瘤相比(中位p4.1月；范围1.3月至13月；心率11.7八；p0.004；95i3.01月46.2月；p0.001；图3)。ea分析强调单剂免疫检查点抑制剂在egfr突变的非小细胞肺癌(59例)、7/八(八八，4/八热、2/八冷、2/八未知)患者中缺乏疗效，这与ea分析一致(补充表s5)。然而，一名免疫790患者(a096)在临床试验中接受了niluabipiliuab免疫检查点抑制剂的联合治疗(60例)，并获得了八.9个月的稳定病情。综上所述，我们的数据提示炎性趋化因子的潜在作用，例如，xl9可能由ss驱动在介导790ki耐药中发挥作用。此外，我们的数据突出了p“热”肿瘤中e成分的显著异质性，说明需要更详细地询问免疫环境以描绘特定的免疫靶点。
虽然第三代egfrkis越来越多地被采用在一线环境中，但这种临床实践在一定程度上是由于无法预测单个患者的耐药轨迹。在确定了与不同egfrki抵抗状态相关的新的分子特征(3q扩增、转录亚型、腺癌谱系标记丢失和炎症的e)后，我们试图建立一个模型来预测790的出现。我们推测，这些基因组、染色体水平和转录特征可能存在于基线水平，也可能代表在治疗过程中获得的变化(图4a)。为了进一步探讨这一点，我们确定了在1g/2gki治疗开始之前可能存在的三个躯干特征：egfr外显子19缺失、g缺失和p53改变(图4a)。使用贝叶斯方法，单个患者可以根据他们治疗前的分子基因型被分成不同的组，其发生egfr790耐药的几率非常不同(图4b)。例如，在非g肿瘤中，发生790耐药性的概率在八7.2到97.9之间，这可能意味着序贯治疗第一代/第二代到第三代egfrki可能是这些患者(我们队列中11的患者)可行的临床策略。这些特征的预测能力需要在更大的队列中进一步验证。然而，这些结果说明了数据驱动的治疗算法是如何通过真实世界的证据得出的，并可能有助于为个别患者定义最佳的排序策略。
我们的研究首次对egfrki耐药的基因组和转录图谱进行了全面和综合的分析。值得注意的是，我们的数据显示，到目前为止，血统的可塑性在一定程度上被低估了。尽管在ki耐药后有1到3的患者有组织转化的报道，但我们发现在790肿瘤中，腺癌标志物(napsina和f1)普遍丢失，同时非乳亚型(pi和pp)明显富集。虽然缺乏配对的基线样本是我们研究的局限性，但与治疗单纯的egfr突变的非小细胞肺癌的比较表明，腺癌谱系标志的丢失，特别是在790疾病中，可能代表了慢性egfrki暴露导致的早期去分化事件。790疾病更显著的基因组改变有p53突变(八6比50)、3q扩增(57比13)和e改变(19比3)，这进一步导致了790的可塑性和耐药性。
临床上特别感兴趣的是790队列中的免疫热亚群，它代表了一组患者的p明显较短，其特征是p评分高，pl1过度表达，以及富含趋化因子的免疫抑制微环境。与我们的发现一致，回顾性分析表明pl1的高表达与低应答率和pfs之间的关系，提示“炎症性”e介导对egfrki的原发性耐药。最近，抑制egfr信号被发现可以耗竭reg和增加ifnγ信号，支持“炎性”的e之间的联系，认为这是一种适应性变化，可能会削弱对靶向治疗的反应。我们的数据进一步表明，炎性的e可能发生在原发或继发耐药时，并由八细胞(肿瘤抗原特异性和/或旁观者)、reg和s可变地组成。最后，观察到i1的高表达，特别是在790免疫热肿瘤中，伴随犬尿氨酸的过度表达(kynu；图2b)，暗示i途径在维持reg激活和在肿瘤亚群中的免疫抑制环境中起作用。最近，通过对一系列癌基因驱动的非小细胞肺癌肿瘤的单细胞rnaseq，aynar和他的同事同样强调了i途径、免疫微环境和肺泡再生细胞特征在靶向治疗中的重要性。我们的数据扩展了这些观察，说明了治疗诱导的适应性细胞
状态可能会受到基因组改变的影响，并表明癌细胞、免疫细胞群和趋化因子中egfr依赖性和谱系可塑性之间存在复杂的相互作用。为了更好地阐明这些免疫介质的作用，前瞻性研究正在进行中。免疫检查点抑制，包括联合抗p1和抗4治疗，在ki耐药环境下的临床试验中显示明显缺乏疗效。除了正在努力评估腺苷轴的免疫抑制靶点，如腺苷2a受体(a2ar)、39和73，未来临床试验的合理靶点可能包括i、s耗竭策略，如贝伐单抗或选择性抑制pi3kγ。
