Front. Oncol.Frontiers in OncologyFront. Oncol.2234-943XFrontiers Media S.A.10.3389/fonc.2022.876424OncologyOriginal ResearchGenetic variants in autophagy-related gene ATG2B predict the prognosis of colorectal cancer patients receiving chemotherapyYuTing1†BenShuai23†MaLing1JiangLu1ChenSilu23LinYu1ChenTao4*LiShuwei23*ZhuLingjun15*1Department of Oncology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China2Department of Genetic Toxicology, The Key Laboratory of Modern Toxicology of Ministry of Education, Center for Global Health, School of Public Health, Nanjing Medical University, Nanjing, China3Department of Environmental Genomics, Jiangsu Key Laboratory of Cancer Biomarkers, Prevention and Treatment, Collaborative Innovation Center for Cancer Personalized Medicine, Nanjing Medical University, Nanjing, China4Department of Gastrointestinal Surgery, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China5Department of Oncology, Sir Run Run Hospital, Nanjing Medical University, Nanjing, China
Edited by: Zhang Honghe, Zhejiang University, China
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Autophagy-related genes have a vital effect on colorectal cancer (CRC) by affecting genomic stability and regulating immune responses. However, the associations between genetic variants in autophagy-related genes and CRC outcomes for chemotherapy therapy remain unclear. The Cox regression model was used to evaluate the associations between single-nucleotide polymorphisms (SNPs) in autophagy-related genes and overall survival (OS) and progression-free survival (PFS) of CRC patients. The results were corrected by the false discovery rate (FDR) correction. We used the logistic regression model to investigate the associations of SNPs with the disease control rate (DCR) of patients. Gene expression analysis was explored based on an in-house dataset and other databases. The associations between gene expression and infiltrating immune cells were evaluated using the Tumor Immune Estimation Resource (TIMER) database. We observed that ATG2B rs17094017 A > T was significantly associated with increased OS (HR = 0.65, 95% CI = 0.50-0.86, P = 2.54×10-3), PFS (HR = 0.76, 95% CI = 0.62-0.93, P = 7.34×10-3), and DCR (OR = 0.60, 95% CI = 0.37-0.96, P = 3.31×10-2) of CRC patients after chemotherapy. The expression of ATG2B was down-expressed in CRC tissues than in adjacent normal tissues. Moreover, ATG2B expression influenced the infiltration of CD8+ T cells, CD4+ T cells, B cells, and T cell receptor signaling pathways, which may inhibit the occurrence of CRC by affecting the immune system. This study suggests that genetic variants in the autophagy-related gene ATG2B play a critical role in predicting the prognosis of CRC prognosis undergoing chemotherapy.
colorectal cancerATG2Bgenetic variantschemotherapyprognosisNatural Science Foundation of Jiangsu Province10.13039/501100004608Jiangsu Provincial Medical Youth Talent10.13039/501100013059Introduction
Colorectal cancer (CRC) is regarded as the third-leading cause of morbidity, as well as the third-leading cause of fatality rates in the United States (1). By 2030, the mortality of individuals with CRC will increase by more than 1.1 million worldwide, while the number of CRC cases will exceed 2.2 million (2). In China, CRC ranks the fifth in cancer‐related mortality and the second in cancer incidence (3). Several risk factors for CRC were identified by epidemiology, such as sex, age, smoking, BMI, heavy drinking, and red and processed meat (4–6). Additionally, the relationships between genetic factors and CRC risk were also verified (7). Chemotherapy, as an effective therapy for CRC, is typically used for treatment after tumor resection. Based on previous studies, a combination of leucovorin and 5-fluorouracil with either irinotecan or oxaliplatin, the first-line treatment regimen, is widely accepted therapy to treat metastatic colorectal cancer (mCRC) (8).
Autophagy is a cellular process in which macromolecules, metabolites, and damaged organelles are transported into lysosomes, where they are degraded and produced as energy materials for cell reuse. Autophagy is commonly subdivided into three types based on its mechanism: macroautophagy, microautophagy, and chaperone-mediated autophagy (9). Autophagy is of great importance in regulating cancer cell metabolism, proliferation, and survival (10, 11). Numerous studies have suggested that autophagy-related genetic variants could accelerate cancer progression (12). The effects of single nucleotide polymorphisms (SNPs) in autophagy-related genes and the outcomes and risk of breast cancer, bladder cancer, non-small cell lung cancer, CRC, and esophageal squamous cell carcinoma were reported (13–18). Recently, autophagy-related genes influence the immune response and the efficacy of immunotherapy by modulating immune system components (19). Autophagy-related genes also have an effect on radiotherapy in CRC cells (20). Furthermore, a survey conducted by Berger et al. (21) demonstrated that SNPs in autophagy-related genes contributed to the occurrence of adverse effects on chemotherapy drugs for CRC.
However, there are no related studies regarding the relationships between genetic variants in autophagy-related genes and the response to CRC chemotherapy in the Chinese population. In view of the evidence that autophagy-related genes can affect cancer prognosis, we assumed that autophagy-related gene variants were related to CRC survival after chemotherapy.
Materials and methodsStudy populations
A cohort of 344 CRC patients who underwent histopathological examination was included in the Affiliated Nanjing First Hospital and the First Affiliated Hospital of Nanjing Medical University in September 2010. According to our follow-up data, 19 patients who did not receive oxaliplatin or irinotecan-based therapy were excluded, and 325 patients were retained. Briefly, 188 CRC patients were treated with oxaliplatin-based regimens, and 137 patients underwent irinotecan-based regimens. A total of 325 CRC patients were followed up using telephone interviewing methods. The deadline for follow-up was April 2, 2016. Details of the study population have been demonstrated (22). Overall survival (OS) after chemotherapy as the primary endpoint was calculated from the time of the first chemotherapy until death or last follow-up for living patients. Additionally, progression-free survival (PFS) refers to the day elapsed from the day of chemotherapy initiation to the day of objective disease progression, death, or last follow-up. OS, PFS, and responses to chemotherapy were considered as outcomes. Peripheral venous blood was collected with written consent for each sample (5 mL). The Institutional Review Board of Nanjing Medical University authorized our research.
Clinical assessment of CRC patients
To assess CRC prognosis before therapy and after a minimum of two cycles of treatment, we used computed tomography as the detection method. Tumor responses to chemotherapy were regarded as the primary endpoint using Response Evaluation Criteria in Solid Tumors (RECIST 1.1). The prognosis and responses to chemotherapy were assessed by the complete response (CR), partial response (PR), progressive disease (PD), and stable disease (SD). The disease control rate (DCR) consisted of CR, PR, and SD.
Autophagy-related genes and SNP selection
Reactome, and Kyoto Encyclopedia of Genes and Genomes were applied to select autophagy-related genes. To find autophagy-related genes more comprehensively, we systematically searched the keywords ‘cancer’, ‘carcinoma’, ‘tumor’, ‘autophagy’, and ‘autophagy-related genes’ in PubMed. In our study, to exclude the effect of gender, we eliminated the genes located on the X chromosome. Briefly, we selected 16 candidate genes for further analysis. We compared the gene expressions between CRC tissues and normal tissues, and genes were selected by these gene screening conditions: (a) fold change > 1.2, (b) P < 0.05 and (c) call rate > 95%. Finally, 8 differentially expressed autophagy-related genes were found in CRC tissues and normal tissues for further study. The schematic diagram of SNP selection was presented in Figure 1. The flow diagram of SNP selection was similar to a previous study (23). Firstly, we extracted SNPs within 2 kb up- and down-stream regions of 8 differentially expressed autophagy-related genes using the Han Chinese in Beijing (CHB) data from the 1000 Genomes Project (March 2012) based on these selection conditions: (a) minor allele frequency (MAF) in population ≥ 0.05, (b) Hardy-Weinberg equilibrium (HWE) ≥ 0.05, (c) call rate > 95%. Secondly, SNPinfo Web Server (24), HaploReg (25), and RegulomeDB (26) were used to predict functional SNPs. SNPs were not included when the RegulomeDB score > 6. Thirdly, we selected tagged SNPs after linkage disequilibrium (LD) (r2 ≥ 0.8) by PLINK 1.09 (27). A total of 19 SNPs were retained in our study. Fourthly, the associations between 19 SNPs and the OS of CRC were assessed in the additive model. The results were corrected by the false discovery rate (FDR), and we selected SNPs of which adjusted PFDR (OS) value was < 0.05. Finally, the effects of the remaining SNPs on PFS and DCR were analyzed in the additive model, and SNPs that were statistically related to PFS and DCR were selected.
Flow diagram of SNP selection in the autophagy-related genes. First, 16 candidate genes were selected in this study. Then, 8 differentially expressed autophagy-related genes were selected by these conditions: (A) fold change > 1.2, (B)P < 0.05 and (C) call rate > 95%. SNPs within 2 kb up- and down-stream regions of 8 candidate genes were extracted using the Han Chinese in Beijing (CHB) data from the 1000 Genomes Project (March 2012) based on these selection conditions: (a) minor allele frequency (MAF) in population ≥ 0.05, (b) Hardy-Weinberg equilibrium (HWE) ≥ 0.05, (c) call rate > 95%. Tagging SNPs were recognized after linkage disequilibrium (LD) (r2 ≥ 0.8). Next, using SNPinfo Web Server, HaploReg, and RegulomeDB to predict functional SNPs. The associations between 19 SNPs and OS of CRC were assessed in the additive model, and SNPs of which adjusted PFDR (OS) value was < 0.05 were selected. Finally, the effects of the remaining SNPs on PFS and DCR were analyzed in the additive model, and SNPs that were statistically related to PFS and DCR were selected. TCGA, The Cancer Genome Atlas; MAF, minor allele frequency; HWE, Hardy‐Weinberg equilibrium; LD, linkage disequilibrium; CHB, the Han Chinese in Beijing; SNP, single nucleotide polymorphism; PFS, progression-free survival; DCR, disease control rate; OS, Overall survival; FDR, false discovery rate.
SNP genotyping
To extract genomic DNA from blood samples of CRC patients, we used the Qiagen Blood Kit (Qiagen). The method of extracting genomic DNA was the same as in a previous study (23). For genotyping in this study, Illumina Human Omni Zhonghua Bead Chips were employed. The microarray is a genome-wide SNP genotyping microarray designed specifically for Chinese people. The microarray covers optimized 900,000 labeled common and rare variants found in the Chinese population, which are strategically selected to explore novel disease and trait associations in the Chinese population. A uniform quality control protocol was used to filter samples and SNPs.
Expression analyses
We compared the expression of candidate genes in CRC tissues and noncancerous tissues using the mRNA expression data from the Gene Expression Omnibus (GEO) dataset (GSE87211), the Cancer Genome Atlas (TCGA) database and in-house RNA-Seq data. To assess targeted gene expression among various subtypes of CRC, we applied the Tumor and Immune System Interaction Database (TISIDB) (http://cis.hku.hk/TISIDB/index.php). We evaluated candidate gene expression in pan‐cancer with the Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/). To confirm the protein level of ATG2B, we used the Human Protein Atlas (HPA) database (https://www.proteinatlas.org/) and the Clinical proteomic tumor analysis consortium (CPTAC) dataset.
Functional analyses
The protein-protein interaction network with autophagy-related genes was explored using the STRING database (https://string-db.org/). The volcano plot was performed to demonstrate differentially expressed genes (DEGs) in high- and low-expressed ATG2B groups. Biological pathways and functions within the subgroups of down-regulated and up-regulated ATG2B expression were assessed using Gene Set Enrichment Analysis (GSEA). Using the Tumor Immune Estimation Resource (TIMER) database (https://cistrome.shinyapps.io/timer/) to investigate the roles of ATG2B expression in immune infiltration.
Statistical analyses
Statistical methods with univariate Cox regression analysis were utilized to assess associations between the clinical characteristics and OS of CRC patients (24). The results were corrected by the FDR, which was employed to mitigate against false-positive results. The associations between targeted SNPs and CRC prognosis were assessed depending on multivariate and univariate COX regression models (24). Logistic regression models were employed to identify the effects of SNPs on DCR. The odds ratios (ORs), hazard ratios (HRs), and 95% confidence intervals (CIs) were calculated for the genetic variants. A two-sided Student’s t test was utilized to analyze the differential expression of genes. A Kaplan-Meier analysis was applied to estimate the survival probability. P values below 0.05 were considered significant. All statistical computations were achieved by R 3.2.3 and PLINK 1.09.
ResultsSubject characteristics
Clinical characteristics of 325 CRC patients and their associations with OS were presented in Supplementary Table 1. However, no significant difference was discovered between these characteristics and OS (P > 0.05). In our study, 205 (63.08%) CRC patients were males, and 120 (36.92%) were females; 59.69% of patients were diagnosed with colon cancer. Patients with poorly differentiated tumors accounted for 21.54% of the total, and 78.46% of CRC patients were in intermediate or advanced tumor stages. Moreover, 7.08% of patients had Dukes stage C disease, and 92.92% were in Dukes stage D.
SNP selection and the genetic effect of SNPs in <italic>ATG2B</italic> on CRC
A total of 16 genes were ultimately selected in our research (Supplementary Table 2). The interactions among these proteins were presented in Supplementary Figure 1. Based on the TCGA database, we selected 8 differentially expressed autophagy-related genes in CRC tissues and normal tissues for further study (Supplementary Table 3). After quality control and LD analysis, a total of 61 SNPs remained. Then, using RegulomeDB, HaploReg, and SNPinfo Web Server to annotate the potential functions of SNPs, we selected 19 SNPs for further analysis (Supplementary Table 4).
The association between rs17094017 in <italic>ATG2B</italic> and CRC prognosis
A total of 19 candidate SNPs were evaluated for associations with the OS of CRC in the additive model. The results illustrated significant associations between rs17094017 and rs11658979 and the OS of CRC (P < 0.05) ( Table 1). After FDR correction, only ATG2B rs17094017 was significantly associated with a favorable OS (PFDR = 4.83×10-2) of CRC. Then, we analyzed whether candidate SNPs were relevant to PFS and DCR in CRC patients. In agreement with previous results, rs17094017 prolonged the PFS (HR = 0.76, 95% CI = 0.62-0.93, P = 7.34×10-3) and increased the DCR (OR = 0.60, 95% CI = 0.37-0.96, P = 3.31×10-2) of CRC (Supplementary Table 5).
Association between selected 19 SNPs and overall survival of clinical patients with CRC.
Minor allele frequencies were calculated by the in-house data.
Effect allele/reference allele.
P value for additive model adjusted for sex, age, smoking status and drinking status in Cox regression model.
The false discovery rate correction of P value.
To evaluate the effects of rs17094017 on ATG2B and CRC survival, we performed Cox regression analyses and logistic regression analyses of four models (dominant model, additive model, codominant and recessive model). Patients with the T allele had an improved prognosis (OS: HR = 0.57, 95% CI = 0.41-0.80, P = 1.03×10-3; PFS: HR = 0.74, 95% CI = 0.56-0.96, P = 2.40×10-2) compared with those harboring the A allele in the dominant model (Table 2). Moreover, a similar result was found in DCR (OR = 0.55, 95% CI = 0.31-0.99, P = 4.43×10-2) (Table 3). In addition, TT genotype carriers had a longer PFS of CRC than AA genotype carriers in the codominant model (HR = 0.53, 95% CI = 0.33-0.87, P = 1.27×10-2) (Table 2). Under the recessive model, we also found that rs17094017 with improvement of PFS was significant (HR = 0.60, 95% CI = 0.37-0.96, P = 3.34×10-2) (Table 2). However, in the recessive model, there were no differences between rs17094017 and OS or DCR (P > 0.05).
Association between ATG2B rs17094017 and the survival of clinical patients with CRC.
Genotypes
Deaths (%)
OS (overall survival)
Progress (%)
PFS (progression-free survival)
HR (95% CI)
P
HR (95% CI)a
Pa
HR (95% CI)
P
HR (95% CI)a
Pa
AA
85 (57.43)
1.00
1.00
116 (51.10)
1.00
1.00
AT
52 (35.14)
0.58 (0.41-0.82)
1.85×10-3
0.58 (0.40-0.82)
2.26×10-3
92 (40.53)
0.77 (0.59-1.02)
6.77×10-2
0.80 (0.60-1.05)
1.10×10-1
TT
11 (7.43)
0.54 (0.28-1.01)
5.27×10-2
0.54 (0.29-1.03)
6.08×10-2
19 (8.37)
0.53 (0.32-0.86)
1.07×10-2
0.53 (0.33-0.87)
1.27×10-2
Additive model
0.65 (0.50-0.85)
1.82×10-3
0.65 (0.50-0.86)
2.54×10-3
0.75 (0.61-0.91)
4.53×10-3
0.76 (0.62-0.93)
7.34×10-3
Dominant model
0.57 (0.41-0.79)
7.57×10-4
0.57 (0.41-0.80)
1.03×10-3
0.72 (0.55-0.93)
1.36×10-2
0.74 (0.56-0.96)
2.40×10-2
Recessive model
0.69 (0.37-1.28)
2.37×10-1
0.70 (0.38-1.31)
2.70×10-1
0.60 (0.37-0.96)
3.36×10-2
0.60 (0.37-0.96)
3.34×10-2
HR, hazard ratio; CI, confidence interval.
aAdjusted for sex, age, smoking and drinking status in Cox regression model.
Association between ATG2B rs17094017 and responses to chemotherapy of clinical patients with CRC.
Genotypes
PD (progress disease) (%)
DCR (disease control rate)
OR (95% CI)
P
OR (95% CI)a
Pa
AA
38 (59.37)
1.00
1.00
AT
23 (35.94)
0.59 (0.33-1.06)
8.00×10-2
0.60 (0.33-1.09)
9.19×10-2
TT
3 (4.69)
0.36 (0.10-1.28)
1.15×10-1
0.35 (0.10-1.25)
1.07×10-1
Additive model
0.60 (0.38-0.95)
3.10×10-2
0.60 (0.37-0.96)
3.31×10-2
Dominant model
0.55 (0.32-0.97)
3.88×10-2
0.55 (0.31-0.99)
4.43×10-2
Recessive model
0.46 (0.13-1.58)
2.15×10-1
0.45 (0.13-1.55)
2.03×10-1
OR, odds ratio; CI, confidence interval.
Adjusted for sex, age, smoking and drinking status in logistic regression model
Moreover, Kaplan-Meier curves were employed to verify the effects of genetic variants on ATG2B rs17094017 and CRC survival in the dominant model. The results revealed that ATG2B rs17094017 A > T might lead to a favorable prognosis (OS: HR = 0.57, 95% CI = 0.41-0.80, P = 1.03×10-3; PFS: HR = 0.74, 95% CI = 0.56-0.96, P = 2.40×10-2) (Figures 2A, B). Using the TCGA database, we compared the OS between ATG2B low-expressed group and high-expressed group. However, no significant difference was verified (P > 0.05) (Figure 2C).
Kaplan-Meier curves of CRC patients. (A–B) Kaplan-Meier curves for overall survival (A) and progression-free survival (B) for rs17094017 in clinical patients with CRC by the Cox regression model. (C) Kaplan-Meier curve for ATG2B expression levels in the TCGA database by the Cox regression model.
Stratified analyses of <italic>ATG2B</italic> rs17094017 in CRC prognosis
To compare the associations between clinical variables and CRC survival after therapy and ATG2B rs17094017, we performed stratification analyses in the dominant model, which included sex, age, cigarette smoking, alcohol intake, tumor site and grade, Dukes stage, metastases, and chemotherapy regimens.
As Figure 3 and Table 4 shown, the AT/TT genotypes had significant associations with a longer OS time than genotype CC in the subgroups of well and moderate tumor grade (HR = 0.56, 95% CI = 0.37-0.82, d P = 3.32×10-3), Dukes stage D (HR = 0.60, 95% CI = 0.42-0.86, P = 5.00×10-3), younger age (HR = 0.44, 95% CI = 0.26-0.73, P = 1.41×10-3), non-smokers (HR = 0.51, 95% CI = 0.33-0.77, P = 1.52×10-3), and non-drinkers (HR = 0.57, 95% CI = 0.38-0.84, P = 4.73×10-3). Furthermore, the AT/TT genotype carriers had improved PFS and increased DCR in non-drinkers, colon cancer, and oxaliplatin-based chemotherapy subgroups (P < 0.05). In addition, rs17094017 was non-significant with PFS and DCR of CRC for subjects below 60. Therefore, there is no sufficient reason for age to be an influential factor in the prognosis of CRC.
Stratified analyses of associations between rs17094017 and survival of clinical patients with CRC under the dominant model. (A–B) Stratified analyses of rs17094017 on overall survival (A) and progression-free survival (B). (C) Stratified analyses of rs17094017 on disease control rate.
Stratified analysis of the association between rs17094017 and the survival of clinical patients with CRC in dominant model.
Adjusted for sex, age, smoking and drinking status in Cox regression model.
Adjusted for sex, age, smoking and drinking status in logistic regression model.
Moreover, CRC prognosis was different among various chemotherapy regimens. Hence, we further used stratification analyses to evaluate whether rs17094017 had an effect on OS, PFS, and DCR of CRC patients with different chemotherapy regimens. The results documented that the AT genotype was relevant to improving prognosis (OS: HR = 0.55, 95% CI = 0.33-0.90, P = 1.70×10-2; PFS: HR = 0.63,95% CI = 0.43-0.93, P = 2.13×10-2) and increasing DCR (OR = 0.40, 95% CI = 0.16-0.97, P = 4.34×10-2) compared with AA genotype carriers receiving oxaliplatin-based chemotherapy (Supplementary Table 6). In the oxaliplatin-based chemotherapy subgroup, rs17094017 T allele prolonged PFS (HR = 0.72, 95% CI = 0.54-0.95, P = 2.15×10-2), OS (HR = 0.64, 95% CI = 0.44-0.94, P = 2.19×10-2), and increased DCR (OR = 0.39, 95% CI = 0.19-0.82, P = 1.26×10-2) under the additive model (Supplementary Table 6). Similarly, the results were discovered in the dominant model (OS: HR = 0.54, 95% CI = 0.34-0.87, P = 1.09×10-2; PFS: HR = 0.62, 95% CI = 0.43-0.90, P = 1.18×10-2; DCR; OR = 0.35, 95% CI = 0.15-0.82, P = 1.60×10-2) (Supplementary Table 6 and Supplementary Figure 2). In the irinotecan-based chemotherapy subgroup of the dominant model, patients with AT/TT genotypes had a longer OS time (HR = 0.61, 95% CI = 0.37-1.00, P = 4.92×10-2) (Supplementary Table 6 and Supplementary Figure 2).
<italic>ATG2B</italic> expression analysis in CRC tissue
We then assessed the ATG2B mRNA expression in 17 paired clinical samples and validated the results in the TCGA and GEO datasets. As presented in Figures 4A–D, ATG2B was lower-expressed in tumor tissues than in noncancerous tissues of the subjects TCGA database, P < 1.00×10-3; TCGA paired data, P < 1.00×10-3; GSE87211, P < 1.00×10-3 and in-house RNA-Seq data, P = 1.10×10-2. The association between ATG2B mRNA expression and various subtypes of CRC was assessed in the TISIB database (colon adenocarcinoma, P = 3.11×10-2; rectal adenocarcinoma, P = 4.95×10-2) (Figures 4 E, F). A similar result was discovered at the protein expression level of ATG2B (Figures 4 G, H). Likewise, we compared the ATG2B mRNA expression in other tumors to adjacent normal tissues using the GEPIA online tool. ATG2B was lowly expressed in most tumor tissues. (Supplementary Figure 3A). Based on the Cancer Cell Line Encyclopedia, ATG2B expression seemed lower in CRC cells compared to other human cancer cell lines (Supplementary Figure 3B).
ATG2B expression in CRC and adjacent normal tissues. The mRNA expression of ATG2B in (A) unpaired tissues from the TCGA database, (B) paired tissues from the TCGA database, (C) GSE87211 database, and (D) in-house RNA-Seq data. ATG2B expression among various subtypes of (E) COAD and (F) READ in the TISIB database. (G) Images of ATG2B expression in normal and CRC tissues visualized by IHC in the HPA database. (H) The protein level of ATG2B in CRC tissues and normal tissues in the CPTAC database. CIN, chromosome instability; GS, Genome Stable; HM-SNV, Hypermutated - single nucleotide variants; HM-Indel, Hypermutated - insertion-deletion; COAD, colon adenocarcinoma; READ, rectal adenocarcinoma.
To further explore the relationships between clinical variables and mRNA expression of ATG2B in CRC tissues, we chose the characteristics of tumor stages, metastasis, sex, age, family history, tumor site, BMI, and KRAS mutation status for analysis in TCGA datasets. The results revealed that ATG2B mRNA expression at each stage in tumor tissues was significantly different from those in noncancerous tissues of the subjects (P < 1.00×10-4) (Supplementary Figure 4A). The mRNA expression of ATG2B in adjacent normal tissues was higher than in CRC tissues with or without metastasis (Supplementary Figure 4B). However, no differences were found in these clinical variables (P > 0.05) (Supplementary Figures 4 C–H).
Functional prediction of <italic>ATG2B</italic>
To validate the relationship between the global gene-expression profile and ATG2B, we compared DEGs with low- and high-expressed ATG2B groups in the TCGA database. The results suggested that there was an up-regulation of 308 genes and a down-regulation of 369 genes (Supplementary Figure 5A). To further explore ATG2B-associated pathways, we analyzed significant differences between ATG2B-high and ATG2B-low expression groups using GSEA analysis in the TCGA datasets. ATG2B was significantly associated with inflammation and immune pathways, such as the Wnt signing pathway, tumor necrosis factor signing pathway, and T cell receptor signing pathway (Supplementary Figure 5B). As shown in Supplementary Figure 6A, ATG2B expression had a positive effect on dendritic cells, B cells, macrophages, CD8+ T cells, neutrophils, and CD4+ T cells in CRC by the TIMER database (P < 0.05). Furthermore, there were significant associations that were found in immune cell infiltration levels under various copy numbers of ATG2B both in colon adenocarcinoma and rectal adenocarcinoma (Supplementary Figure 6B).
Discussion
In our study, the association between SNPs in autophagy-related genes and CRC outcome for chemotherapy treatment were assessed. ATG2B rs17094017 A > T had a significant effect on the prolonged OS and PFS of CRC with oxaliplatin-based chemotherapy. Interestingly, ATG2B rs17094017 T allele was associated with DCR in CRC patients. It was concluded that rs17094017 A > T in ATG2B might predict a better prognosis for CRC patients with oxaliplatin-based chemotherapy.
Chemotherapy is one of the primary treatments for CRC. However, there is a significant difference in individual responses to the efficacy of chemotherapy. Growing evidence revealed that autophagy played a great role in mediating resistance to chemotherapy. Melanoma patients who were resistant to the BRAF inhibitor showed higher levels of autophagy through the endoplasmic reticulum stress response (25). A study revealed that inhibition of intestinal epithelial autophagy through intestinal flora could improve CRC patients’ responses to chemotherapy and alter outcomes (26). In ovarian cancer, the resistance of the cytotoxic drug paclitaxel has been attributed to autophagy induction (27).
Autophagy-related gene variants were believed to be closely related to the development of cancers, and core genes were considered to affect the functions of the cells in metabolism, proliferation, apoptosis, and immunity (28, 29). Recently, the relationships between autophagy-related gene variants and cancer prognosis have been evaluated. For example, the association between rs473543 in ATG5 and disease-free survival (DFS) of breast cancer patients undergoing chemotherapy was reported (30). Recent evidence also suggested a significant association between ATG2B rs17784271 and poor local recurrence-free survival and PFS in non-small cell lung cancer after radiotherapy (31). ATG2B rs3759601 has been shown to have therapeutic effects on bladder cancer treated with Bacillus Calmette-Guerin (16). A truncated variant of UVRAG was related to the transformation and tumor metastasis of CRC (32). ATG16L1 T300A has an effect on a good prognosis in CRC (17). Moreover, FIP200 rs1129660 played a crucial role in bevacizumab-mediated toxicity of mCRC (21). However, in the Chinese Han population, no study involving the relationships of SNPs in autophagy-related genes with CRC survival undergoing chemotherapy was found. Our study is the first to explore the association between ATG2B rs17094017 and the prognosis of CRC patients receiving chemotherapy among the Chinese population.
ATG2B, located on chromosome 14q32.2, is necessary for forming autophagosomes (33). The roles of ATG2A, ATG2B, and WIPI proteins are crucial for the membrane extension of the PI3P formation site in autophagy (12). It has been previously demonstrated that frameshift mutations of ATG2B with mononucleotide repeats occur in both CRC and gastric cancer (34). Additionally, our study revealed that ATG2B was lowly expressed in CRC tissues compared to noncancerous tissues. Reported research indicated that ATG2B displayed low expression in breast cancer (35). ATG2B was down-regulated in inflamed tissues compared to adjacent noninflamed tissues in Crohn’s disease (36). The downregulation of ATG2B activated cancer-associated fibroblasts (CAFs) by inhibiting autophagy in P53-deficient status, which accelerated the proliferation of CRC cells (37). Based on our results, the inhibition of autophagy promoted the development of cancer cells and affected chemosensitivity through suppressing ATG2B expression.
Moreover, autophagy plays an essential role in anti-tumor by regulating the immune system (38). Autophagy gene deficiency interferes with the survival, development, and differentiation of T and B cells (39, 40). Additionally, autophagy-mediated lipolysis inhibited the mitochondrial oxidative respiration pathway of neutrophil differentiation, which further affected immune defense mechanisms leading to cancer (41). Therefore, we further used GSEA analysis to predict the enrichment pathway of ATG2B. The results revealed that ATG2B was significantly related to the inflammatory and immune pathways. Based on TIMER databases, a positive correlation was discovered between ATG2B and tumor-infiltrating lymphocytes, suggesting that ATG2B might influence the immunotherapy of CRC modulating the tumor-infiltrating immune cells, but functional studies of ATG2B still need further verification.
CRC, a complex disease, which is ascribed to lifestyle and genetic factors (7). Hence, stratification analysis was used to obtain the association between clinical characteristics and CRC prognosis in this study. The results illustrated that rs17094017 was related to improved prognosis and increased DCR in colon cancer patients without drinking. A prospective cohort study suggested that patients with oxaliplatin-based chemotherapy had a longer OS time than those receiving irinotecan-based chemotherapy (42). In our study, we discovered that rs17094017 was meaningfully related to prolonged PFS, OS, and increased DCR in oxaliplatin-based chemotherapy rather than in irinotecan-based chemotherapy.
There are still several limitations to the research. Firstly, the sample size of the CRC population was relatively small. Thus, a larger population with more complete survival data is required to verify these results. Secondly, evidence concerning biological experiments with ATG2B is lacking. Hence, more basic biological research is required to further illustrate the effects of ATG2B in CRC.
Taken together, this study indicated that rs17094017 in ATG2B was related to a better outcome in CRC patients receiving chemotherapy. Furthermore, our study revealed that ATG2B rs17094017 was associated with increased DCR in CRC patients after treatment. SNP rs17094017 could serve as a novel biomarker to predict a CRC patient’s prognosis undergoing chemotherapy, providing a theoretical basis for individualized therapy in CRC. In summary, the association of genetic effects in ATG2B and the survival of CRC patients receiving chemotherapy was first explored among the Chinese population.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by The Institutional Review Board of Nanjing Medical University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
LZ, SL, and TC designed the study. TY and SB wrote the manuscript. SL critically revised the paper. SC assisted in the data analysis. LM, LJ, and YL collected the samples. All authors contributed to the article and approved the submitted version.
Funding
Natural Science Foundation of Jiangsu Province (BK20201495), Jiangsu Provincial Medical Talent (ZDRCA2016089), and National Natural Science Foundation of China (No. 82102981) supported our research.
Acknowledgments
We wish to thank Prof. Meilin Wang for providing data.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.876424/full#supplementary-material
ReferencesSiegelRLMillerKDFuchsHEJemalA. Cancer statistics, 2021. CA Cancer J Clin (2021) 71(1):7–33. doi: 10.3322/caac.21654ArnoldMSierraMSLaversanneMSoerjomataramIJemalABrayF. Global patterns and trends in colorectal cancer incidence and mortality. Gut (2017) 66(4):683–91. doi: 10.1136/gutjnl-2015-310912CaoWChenHDYuYWLiNChenWQ. Changing profiles of cancer burden worldwide and in china: a secondary analysis of the global cancer statistics 2020. Chin Med J (Engl) (2021) 134(7):783–91. doi: 10.1097/CM9.0000000000001474WolfAMDFonthamETHChurchTRFlowersCRGuerraCELaMonteSJ. Colorectal cancer screening for average-risk adults: 2018 guideline update from the american cancer society. CA Cancer J Clin (2018) 68(4):250–81. doi: 10.3322/caac.21457CaiSLiYDingYChenKJinM. Alcohol drinking and the risk of colorectal cancer death: a meta-analysis. Eur J Cancer Prev (2014) 23(6):532–9. doi: 10.1097/CEJ.0000000000000076DekkerETanisPJVleugelsJLAKasiPMWallaceMB. Colorectal cancer. Lancet (2019) 394(10207):1467–80. doi: 10.1016/s0140-6736(19)32319-0KastrinosFSamadderNJBurtRW. Use of family history and genetic testing to determine risk of colorectal cancer. Gastroenterology (2020) 158(2):389–403. doi: 10.1053/j.gastro.2019.11.029BennounaJHiretSBertautABoucheODeplanqueGBorelC. Continuation of bevacizumab vs cetuximab plus chemotherapy after first progression in kras wild-type metastatic colorectal cancer: the unicancer prodige18 randomized clinical trial. JAMA Oncol (2019) 5(1):83–90. doi: 10.1001/jamaoncol.2018.4465MizushimaNKomatsuM. Autophagy: renovation of cells and tissues. Cell (2011) 147(4):728–41. doi: 10.1016/j.cell.2011.10.026DowerCMWillsCAFrischSMWangHG. Mechanisms and context underlying the role of autophagy in cancer metastasis. Autophagy (2018) 14(7):1110–28. doi: 10.1080/15548627.2018.1450020KimmelmanACWhiteE. Autophagy and tumor metabolism. Cell Metab (2017) 25(5):1037–43. doi: 10.1016/j.cmet.2017.04.004LevineBKroemerG. Biological functions of autophagy genes: a disease perspective. Cell (2019) 176(1-2):11–42. doi: 10.1016/j.cell.2018.09.048WenJLiuHWangLWangXGuNLiuZ. Potentially functional variants of atg16l2 predict radiation pneumonitis and outcomes in patients with non–small cell lung cancer after definitive radiotherapy. J Thorac Oncol (2018) 13(5):660–75. doi: 10.1016/j.jtho.2018.01.028XieKLiangCLiQYanCWangCGuY. Role of atg10 expression quantitative trait loci in non-small cell lung cancer survival. Int J Cancer (2016) 139(7):1564–73. doi: 10.1002/ijc.30205Vera-RamirezLVodnalaSKNiniRHunterKWGreenJE. Autophagy promotes the survival of dormant breast cancer cells and metastatic tumour recurrence. Nat Commun (2018) 9(1):1944. doi: 10.1038/s41467-018-04070-6BuffenKOostingMQuintinJNgAKleinnijenhuisJKumarV. Autophagy controls bcg-induced trained immunity and the response to intravesical bcg therapy for bladder cancer. PloS Pathog (2014) 10(10):e1004485. doi: 10.1371/journal.ppat.1004485GrimmWAMesserJSMurphySFNeroTLodolceJPWeberCR. The thr300ala variant in atg16l1 is associated with improved survival in human colorectal cancer and enhanced production of type i interferon. Gut (2016) 65(3):456–64. doi: 10.1136/gutjnl-2014-308735ZhuMLYuHShiTYHeJWangMYLiQX. Polymorphisms in mtorc1 genes modulate risk of esophageal squamous cell carcinoma in eastern chinese populations. J Thorac Oncol (2013) 8(6):788–95. doi: 10.1097/JTO.0b013e31828916c6YoungTMReyesCPasnikowskiECastanaroCWongCDeckerCE. Autophagy protects tumors from t cell-mediated cytotoxicity via inhibition of tnfalpha-induced apoptosis. Sci Immunol (2020) 5(54). doi: 10.1126/sciimmunol.abb9561YukJMShinDMSongKSLimKKimKHLeeSH. Bacillus calmette-guerin cell wall cytoskeleton enhances colon cancer radiosensitivity through autophagy. Autophagy (2010) 6(1):46–60. doi: 10.4161/auto.6.1.10325BergerMDYamauchiSCaoSHannaDLSunakawaYSchirripaM. Autophagy-related polymorphisms predict hypertension in patients with metastatic colorectal cancer treated with folfiri and bevacizumab: results from tribe and fire-3 trials. Eur J Cancer (2017) 77:13–20. doi: 10.1016/j.ejca.2017.02.020ZhangKMengYCaoXXuYDuMWuY. Genetic variants in P53 signaling pathway genes predict chemotherapy efficacy in colorectal cancer. Cancer Med (2019) 8(7):3428–36. doi: 10.1002/cam4.2215WangMGuDDuMXuZZhangSZhuL. Common genetic variation in etv6 is associated with colorectal cancer susceptibility. Nat Commun (2016) 7:11478. doi: 10.1038/ncomms11478KumarSPatnaikSDixitA. Predictive models for stage and risk classification in head and neck squamous cell carcinoma (Hnscc). PeerJ (2020) 8:e9656. doi: 10.7717/peerj.9656MaXHPiaoSFDeySMcAfeeQKarakousisGVillanuevaJ. Targeting er stress-induced autophagy overcomes braf inhibitor resistance in melanoma. J Clin Invest (2014) 124(3):1406–17. doi: 10.1172/JCI70454YuTGuoFYuYSunTMaDHanJ. Fusobacterium nucleatum promotes chemoresistance to colorectal cancer by modulating autophagy. Cell (2017) 170(3):548–63.e16. doi: 10.1016/j.cell.2017.07.008ZhangSFWangXYFuZQPengQHZhangJYYeF. Txndc17 promotes paclitaxel resistance via inducing autophagy in ovarian cancer. Autophagy (2015) 11(2):225–38. doi: 10.1080/15548627.2014.998931XieXKohJYPriceSWhiteEMehnertJM. Atg7 overcomes senescence and promotes growth of brafv600e-driven melanoma. Cancer Discovery (2015) 5(4):410–23. doi: 10.1158/2159-8290.CD-14-1473CadwellK. Crosstalk between autophagy and inflammatory signalling pathways: balancing defence and homeostasis. Nat Rev Immunol (2016) 16(11):661–75. doi: 10.1038/nri.2016.100LiMMaFWangJLiQZhangPYuanP. Genetic polymorphisms of autophagy-related gene 5 (atg5) rs473543 predict different disease-free survivals of triple-negative breast cancer patients receiving anthracycline- and/or taxane-based adjuvant chemotherapy. Chin J Cancer (2018) 37(1):4. doi: 10.1186/s40880-018-0268-1WenJLiuHWangLWangXGuNLiuZ. Potentially functional variants of atg16l2 predict radiation pneumonitis and outcomes in patients with non-small cell lung cancer after definitive radiotherapy. J Thorac Oncol (2018) 13(5):660–75. doi: 10.1016/j.jtho.2018.01.028HeSZhaoZYangYO'ConnellDZhangXOhS. Truncating mutation in the autophagy gene uvrag confers oncogenic properties and chemosensitivity in colorectal cancers. Nat Commun (2015) 6:7839. doi: 10.1038/ncomms8839ZhengJ-XLiYDingY-HLiuJ-JZhangM-JDongM-Q. Architecture of the atg2b-wdr45 complex and an aromatic y/hf motif crucial for complex formation. Autophagy (2017) 13(11):1870–83. doi: 10.1080/15548627.2017.1359381KangMRKimMSOhJEKimYRSongSYKimSS. Frameshift mutations of autophagy-related genes atg2b, atg5, atg9b and atg12 in gastric and colorectal cancers with microsatellite instability. J Pathol (2009) 217(5):702–6. doi: 10.1002/path.2509ZhangXLiCWangDChenQLiC-LLiH-J. Aberrant methylation of atg2b , atg4d , atg9a and atg9b cpg island promoter is associated with decreased mrna expression in sporadic breast carcinoma. Gene (2016) 590(2):285–92. doi: 10.1016/j.gene.2016.05.036LinXTZhengXBFanDJYaoQQHuJCLianL. Microrna-143 targets atg2b to inhibit autophagy and increase inflammatory responses in crohn's disease. Inflammation Bowel Dis (2018) 24(4):781–91. doi: 10.1093/ibd/izx075InoueTHayashiYTsujiiYYoshiiSSakataniAKimuraK. Suppression of autophagy promotes fibroblast activation in p53-deficient colorectal cancer cells. Sci Rep (2021) 11(1):19524. doi: 10.1038/s41598-021-98865-1JiangGMTanYWangHPengLChenHTMengXJ. The relationship between autophagy and the immune system and its applications for tumor immunotherapy. Mol Cancer (2019) 18(1):17. doi: 10.1186/s12943-019-0944-zPuaHHGuoJKomatsuMHeYW. Autophagy is essential for mitochondrial clearance in mature t lymphocytes. J Immunol (2009) 182(7):4046–55. doi: 10.4049/jimmunol.0801143ArsovIAdebayoAKucerova-LevisohnMHayeJMacNeilMPapavasiliouFN. A role for autophagic protein beclin 1 early in lymphocyte development. J Immunol (2011) 186(4):2201–9. doi: 10.4049/jimmunol.1002223RiffelmacherTClarkeARichterFCStranksAPandeySDanielliS. Autophagy-dependent generation of free fatty acids is critical for normal neutrophil differentiation. Immunity (2017) 47(3):466–80.e5. doi: 10.1016/j.immuni.2017.08.005MarschnerNArnoldDEngelEHutzschenreuterURauhJFreierW. Oxaliplatin-based first-line chemotherapy is associated with improved overall survival compared to first-line treatment with irinotecan-based chemotherapy in patients with metastatic colorectal cancer - results from a prospective cohort study. Clin Epidemiol (2015) 7:295–303. doi: 10.2147/CLEP.S73857