国际肿瘤学杂志 ›› 2023, Vol. 50 ›› Issue (10): 592-599.doi: 10.3760/cma.j.cn371439-20220916-00113

• 论著 • 上一篇    下一篇

全外显子组测序筛查食管鳞状细胞癌放疗抵抗相关基因研究

陈志明1, 陈俊杰2, 李李1, 丁倩1, 韩钰楠1, 赵洪瑜1()   

  1. 1南通大学附属医院肿瘤放疗科,南通 226001
    2南通大学附属医院临床医学研究中心,南通 226001
  • 收稿日期:2022-09-16 修回日期:2023-08-08 出版日期:2023-10-08 发布日期:2023-11-08
  • 通讯作者: 赵洪瑜 E-mail:z_hy07@126.com
  • 基金资助:
    南通市基础科学研究项目(JC2020071)

Identification of critical radioresistance genes in esophageal squamous cell carcinoma by whole exome sequencing

Chen Zhiming1, Chen Junjie2, Li Li1, Ding Qian1, Han Yunan1, Zhao Hongyu1()   

  1. 1Department of Radiotherapy & Oncology, Affiliated Hospital of Nantong University, Nantong 226001, China
    2Clinical Medical Center of Affiliated Hospital of Nantong University, Nantong 226001, China
  • Received:2022-09-16 Revised:2023-08-08 Online:2023-10-08 Published:2023-11-08
  • Contact: Zhao Hongyu E-mail:z_hy07@126.com
  • Supported by:
    Nantong Basic Science Research Project(JC2020071)

摘要:

目的 比较食管鳞状细胞癌(ESCC)术后放疗后不同预后患者的基因谱差异,筛选与放疗抵抗相关的遗传变异。方法 选择2015年1月至2019年12月在南通大学附属医院接受根治性手术及术后辅助放疗的32例ESCC患者为研究对象,根据1年内是否出现照射野内复发分为复发组(放疗抵抗组,n=16)和稳定组(放疗敏感组,n=16)。分别提取患者的基因组DNA,利用全外显子组测序(WES)技术进行高通量测序。应用Trimmomatic、BWA、Picard生物信息学分析软件处理数据,通过GATK对比获得比对文件,然后利用Vardict软件从测序数据中筛选出两组的各类遗传变异。采用Kaplan-Meier法估算患者无瘤生存期(DFS)、总生存期(OS)。Cox比例风险回归模型分析影响ESCC患者DFS和OS的独立危险因素。结果 经样本数据质量控制,本研究最终纳入26例患者进行后续分析,复发组和稳定组各13例。全组非沉默肿瘤突变负荷中位数为0.95个/Mb,突变碱基替换类型均以C>T的转换为主,其次为C>G的颠换;发生频率最高的遗传变异依次是单核苷酸多态性(SNP)(75.1%)、缺失突变(13.7%)、插入突变(10.5%);复发组中肿瘤特有突变数目较稳定组稍高(中位突变数分别为36、34),且两组突变频率前十的基因谱明显不同。复发组中检测到392个特有突变基因,前5个分别为:MUC19、NPIPA5、EPPK1、FLG和FOXG1。稳定组检测到192个特有突变基因,前5个分别为TCHH、WNK1、AIM1L、COL6A5和DPCR1。复发组中位DFS和中位OS分别为15.0个月(95%CI为10.1个月~未达到)和26.2个月(95%CI为19.8个月~未达到),稳定组患者均未出现复发或转移。单因素分析发现,GRIK2(χ2=6.81,P=0.009)、MUC4(χ2=4.25,P=0.039)、MUC5B(χ2=4.03,P=0.045)、PRRG1(χ2=5.15,P=0.023)基因突变、3p缺失(χ2=4.16,P=0.041)和14q缺失(χ2=7.09,P=0.008)与DFS相关;FLG(χ2=6.41,P=0.011)、NPIPA5(χ2=4.57,P=0.033)、PKD1L2(χ2=6.41,P=0.011)、FOXG1(χ2=4.57,P=0.033)基因突变、3p缺失(χ2=3.88,P=0.049)、14q缺失(χ2=5.66,P=0.017)和18p缺失(χ2=3.85,P=0.050)与OS相关。多因素分析发现,14q缺失(HR=3.65,95%CI为1.18~11.32,P=0.025)是影响术后辅助放疗ESCC患者DFS的独立危险因素,FLG(HR=8.94,95%CI为1.52~52.74,P=0.016)、NPIPA5(HR=6.36,95%CI为1.23~33.03,P=0.028)基因突变和14q缺失(HR=3.82,95%CI为1.18~12.31,P=0.025)是影响术后辅助放疗ESCC患者OS的独立危险因素。结论 WES结果提示,术后辅助放疗ESCC复发组与稳定组的基因突变类型和突变率基本一致,但两组的基因突变谱明显不同。FLG、NPIPA5基因突变和14q缺失可作为预测术后辅助放疗ESCC患者预后的分子标志物。

关键词: 食道鳞癌, 全外显子组测序, 遗传变异, 放疗抵抗

Abstract:

Objective To compare the genetic spectrums of esophageal squamous cell carcinoma (ESCC) patients with different prognosis after postoperative radiotherapy and to screen the genetic variants associated with radiotherapy resistance. Methods A total of 32 ESCC patients who received radical surgery and postoperative adjuvant radiotherapy in Affiliated Hospital of Nantong University from January 2015 to December 2019 were selected as the study objects. According to whether there was any recurrence in the radiation field within 1 year, they were divided into a recurrence group (radiotherapy resistance group, n=16) and a stable group (radiotherapy sensitive group, n=16). Genomic DNA was extracted from patients and high-throughput sequencing was performed using whole exome sequencing (WES) technology. Biological information analysis software Trimmomatic, BWA and Picard were used to process the data and the alignment files were obtained by GATK comparison, then Vardict software was used to screen out various genetic variants from the sequencing data. The disease free survival (DFS) and overall survival (OS) were estimated by Kaplan-Meier method. Cox proportional hazard regression model was used to analyze the independent risk factors of DFS and OS of ESCC patients. Results After quality control of the sample data, 26 patients were finally included in this study for follow-up analysis, 13 in each of the recurrence and stable groups. The median tumor mutation burden of non-silent tumors in the whole group was 0.95 mutations/Mb. The substitution types of mutant bases were mainly C>T conversion, followed by C>G transmutation. The genetic variants with the highest frequency were single nucleotide polymorphism (SNP) (75.1%), deletion mutation (13.7%) and insertion mutation (10.5%). The number of tumor-specific mutations in the recurrence group was slightly higher than that in the stable group (median mutation number was 36 and 34, respectively), and the top ten gene profiles of mutation frequency were significantly different between the two groups. In the recurrence group, 392 unique mutated genes were detected, and the top five were MUC19, NPIPA5, EPPK1, FLG and FOXG1. In the stable group, 192 unique mutation genes were detected, and the top five were TCHH, WNK1, AIM1L, COL6A5 and DPCR1. The median DFS and OS were 15.0 months (95%CI: 10.1 months-not reached) and 26.2 months (95%CI: 19.8 months-not reached) in the recurrence group respectively, and no recurrence or metastasis occurred in the stable group. Univariate analysis showed that GRIK2 (χ2=6.81, P=0.009), MUC4 (χ2=4.25, P=0.039), MUC5B (χ2=4.03, P=0.045), PRRG1 (χ2=5.15, P=0.023) gene mutations, 3p deletion (χ2=4.16, P=0.041) and 14q deletion (χ2=7.09, P=0.008) were correlated with DFS. FLG (χ2=6.41, P=0.011), NPIPA5 (χ2=4.57, P=0.033), PKD1L2 (χ2=6.41, P=0.011), FOXG1 (χ2=4.57, P=0.033) gene mutations, 3p deletion (χ2=3.88, P=0.049), 14q deletion (χ2=5.66, P=0.017) and 18p deletion (χ2=3.85, P=0.050) were associated with OS. Multivariate analysis showed that 14q deletion (HR=3.65, 95%CI: 1.18-11.32, P=0.025) was an independent risk factor for DFS of ESCC patients with postoperative adjuvant radiotherapy, and FLG (HR=8.94, 95%CI: 1.52-52.74, P=0.016), NPIPA5 (HR=6.36, 95%CI: 1.23-33.03, P=0.028) gene mutation and 14q deletion (HR=3.82, 95%CI: 1.18-12.31, P=0.025) were independent risk factors for OS of ESCC patients with postoperative adjuvant radiotherapy. Conclusion The WES results suggest that the types and rates of gene mutations of the ESCC patients with postoperative adjuvant radiotherapy in the recurrence and stable groups are basically the same, but the mutation spectrum of the two groups is significantly different. FLG, NPIPA5 gene mutations and 14q deletion can be used as molecular markers to predict the prognosis of ESCC patients treated with postoperative adjuvant radiotherapy.

Key words: Esophageal squamous cell carcinoma, Whole exome sequencing, Genetic variation, Radiation resistance