A two-stage genome-wide association study identifies novel germline genetic variations in CACNA2D3 associated with radiotherapy response in nasopharyngeal carcinoma
Journal of Translational Medicine volume 21, Article number: 11 (2023)
Radiotherapy (RT) is the standard treatment for nasopharyngeal carcinoma (NPC). However, due to individual differences in radiosensitivity, biomarkers are needed to tailored radiotherapy to cancer patients. However, comprehensive genome-wide radiogenomic studies on them are still lacking. The aim of this study was to identify genetic variants associated with radiotherapy response in patients with NPC.
This was a large‑scale genome-wide association analysis (GWAS) including a total of 981 patients. 319 individuals in the discovery stage were genotyped for 688,783 SNPs using whole genome-wide screening microarray. Significant loci were further genotyped using MassARRAY system and TaqMan SNP assays in the validation stages of 847 patients. This study used logistic regression analysis and multiple bioinformatics tools such as PLINK, LocusZoom, LDBlockShow, GTEx, Pancan-meQTL and FUMA to examine genetic variants associated with radiotherapy efficacy in NPC.
After genome-wide level analysis, 19 SNPs entered the validation stage (P < 1 × 10− 6), and rs11130424 ultimately showed statistical significance among these SNPs. The efficacy was better in minor allele carriers of rs11130424 than in major allele carriers. Further stratified analysis showed that the association existed in patients in the EBV-positive, smoking, and late-stage (III and IV) subgroups and in patients who underwent both concurrent chemoradiotherapy and induction/adjuvant chemotherapy.
Our study showed that rs11130424 in the CACNA2D3 gene was associated with sensitivity to radiotherapy in NPC patients.
Trial registration number: Effect of genetic polymorphism on nasopharyngeal carcinoma chemoradiotherapy reaction, ChiCTR-OPC-14005257, Registered 18 September 2014, http://www.chictr.org.cn/showproj.aspx?proj=9546.
Nasopharyngeal carcinoma is a malignant tumor arising from the endothelium of the nasopharynx  Its incidence is highly geographically specific, with a high prevalence in countries of China, India and Vietnam. China accounted for approximately half of all new cases of NPC worldwide in 2018 .
The main treatment for NPC is radiotherapy  The current radiation treatment strategy for NPC is to use simple radiotherapy for early-stage NPC and concurrent chemoradiotherapy for intermediate and advanced stages  Although radiotherapy is the major treatment for NPC, there is some individual difference in the outcome of patients receiving radiotherapy, and the exact mechanism has not yet been fully understood. There is currently a problem of over- or under-treatment in clinical care as there is no good way to accurately identify patients who may benefit from radiotherapy prior to treatment. In the treatment of NPC, people seek to improve the efficacy of radiotherapy and increase patient survival.
With regard to the development of biomarkers related to the efficacy of radiotherapy, the concept of radiogenomics has been proposed  Radiogenomics is based on the theory that the efficacy of radiotherapy is not only related to the treatment scheme, radiation dose and irradiation location, but also to the genetic characteristics of the patients. Candidate gene association and genome-wide association analysis are the main approaches in radiogenomics studies. Most current studies about biomarkers of tumor radiotherapy efficacy are based on candidate gene association studies. For example, X-ray repair cross complementing 1 (XRCC1) 399Gln was found to be significantly associated with radiotherapy resistance in patients with oesophageal carcinoma ; Epidermal growth factor receptor (EGFR) R497K can be used as a predictor of short-term efficacy and late recurrence or metastasis of radiotherapy in cervical cancer  In addition, rs3135001 is a useful biomarker for predicting for the sensitivity of head and neck squamous cell carcinoma to radiotherapy  However, there are few GWAS on NPC radiotherapy sensitivity, but it is equally important to explore the genetic variants associated with NPC radiotherapy sensitivity.
In this study, we conducted a two-stage GWAS study in 1166 NPC patients to comprehensively examine the entire genome at a higher resolution than candidate gene association studies. We selected multiple NPC susceptibility-associated loci and validated these findings in two additional cohort of samples. Our findings clearly indicate that the CACNA2D3 gene plays a role in the mechanism of radiotherapy sensitivity in NPC.
Study population and data collection
All patients in this study were recruited from Jiangxi Provincial Cancer Hospital and Hunan Provincial Cancer Hospital from 2014 to 2022 respectively, who were diagnosed with NPC by pathological histological examination, while having no distant metastases and no other concurrent tumors. They are mainly treated by Intensity Modulated Radiation Therapy (IMRT), which allows for more precise application of different doses of radiation to different target areas. The radiation dose depends on the total tumor volume (GTV), the clinical target volume (CTV) and the planned target volume (PTV). A single dose of radiation therapy is usually about 2 Gy, once a day, 5 days a week. Primary lesions and positive cervical lymph nodes received 66–77 Gy and 54–60 Gy, respectively, and patients were treated in approximately 30–33 sessions. The patients have not received surgery, targeted therapy, immunotherapy or other anti-tumor treatment prior to radiotherapy. Cisplatin was used as the base regimen for adjuvant chemotherapy, neoadjuvant chemotherapy or concurrent radiotherapy. Exclusion criteria are: (1) breastfeeding pregnant patients; (2) having serious complications.
Clinicians used the WHO Response Evaluation Criteria in Solid Tumors (RECIST) to assess the efficacy of each patient after radiotherapy based on the patient’s MRI medical records. In this study, post-radiotherapy efficacy was assessed for the primary lesion and positive lymph nodes. The results of radiotherapy efficacy assessment were classified into four categories: complete remission (CR), partial remission (PR), stable disease (SD) and progressive disease (PD). For association analysis, all patients were further divided into two groups: the group with good outcome as CR, and the group with poor outcome as non-CR (including PR, SD and PD).
The study protocol was approved by the Ethics Committee of Jiangxi Cancer Hospital and Hunan Cancer Hospital. Clinical information of all patients from both hospitals were collected with the same protocol and were not any overlap of patients. The ethical aspects of this study are in line with the Declaration of Helsinki. All subjects were provided with written informed consents. We applied this study for clinical admission in the Chinese Clinical Trial Register (registration number: ChiCTR-OPC-14005257).
Genotype detection, quality control and imputation
For discovery stage samples, genotyping was performed using the Illumina Infinium Global Screening Array-24 v1.0 (GSA) microarray, and the data generated by the Illumina microarray platform was visualized and analyzed according to GenomeStudio software. For quality control, the raw data were first selected for sample and SNP selection using PLINK (version 1.9) software. SNP quality control (QC) conditions: (1) SNP call rate ≥ 95%; (2) HWE ≥ 1 × 10− 4; (3) MAF ≥ 0.01; (4) localization on autosomes. Sample QC conditions: (1) Sample call rate ≥ 95%; (2) Correcting for Cryptic relatedness. Cryptic relatedness refers to unknown more recent relatedness including family relationships such as grandparent-grandchild and full sibling pairs; (3) The heterozygosity rate was less than three standard deviations from the mean of all participants; (4) Genotype-estimated sex was discordant with biological sex. Using 1000 Genomes stage 3 v4 as the reference panel, IMPUTE2 was used for genotype imputation. SNPs with INFO below 0.4 and MAF<0.05 were excluded. Finally, 4,112,760 unique SNPs were included in the GWAS analysis. For samples in the validation stage 1, genotyping was performed using the Sequenom MassARRAY system, and genotyping data were obtained using TYPER 4.0 software. Exclusion criteria at this stage included (1) Sample call rate ≥ 95%; (2) SNP detection rate ≥ 95%; (3) MAF ≥ 0.01; (4) HWE ≥ 0.05; (5) Sample relatedness check; (6) The heterozygosity rate was less than three standard deviations from the mean of all participants; and (7) Genotype-estimated sex was discordant with biological sex. For samples in the validation stage 2, genotyping was conducted using the TaqMan SNP assays (Applied Biosystems, Foster City, CA) with a Roche LightCycler 480 Real-Time PCR System (Roche, Basle, Switzerland) for rs11130424. The amplification reaction was carried out in 10 µL of the PCR mixture containing 5 µL of LightCycler® 480 Probes Master (Roche, Switzerland), 0.5 µM of forward primer rs11130424-F: 5′- ATCATGGACGGTCCTTGCT-3′ and 0.5 µM of the reverse primer rs11130424-R: 5′- GGACATGTCAATTTAAAGTTGCTAT − 3′, 0.4 µM of each allele-specific probes (rs11130424_G: 5′-FAM- AAACTCTGTTACTCATAGAA -MGB-3′ and rs11130424_T:5′-VIC- AACTCTGTTATTCATAGAAT -MGB-3′) and 10 ng of DNA. PCR was carried out using a LightCycler® 96 System (Roche, Switzerland) under optimized conditions (preliminary denaturation – 94 ℃–3 min; 40 cycles: 94 ℃–5 s, 60 ℃–30 s). Fluorescence measurements were carried out at the annealing step of PCR in FAM and VIC channels. The results of genotyping were analyzed using the LightCycler® 96 thermal cycler software.
Genome-wide association analysis
Principal component analysis (PCA) was conducted on all patients in the discovery stage using the SNPRelate procedure. The distribution of the samples was plotted using the R package “ggplot2”. Samples that deviated significantly from the overall were removed and then reconducted for PCA. Clinical factors that have a significantly effect on the phenotype were considered as covariates. T-tests were used for continuous variables and logistic tests for categorical variables, and clinical factors that differed markedly in subgroups were selected as covariates for correction. The significantly associated clinical factors in multivariate regression were adjusted in genome-wide association analyses. For two phenotypes of primary lesion and positive lymph node efficacy, genome-wide association analyses were performed under additive genetic effects assumption, using a logistic regression model adjusting sex, smoking and drinking as covariates. Quantitative-Quantitative plots and Manhattan plots were plotted via the R package “qqman”. A genomic inflation factor λ was calculated for population stratification using PLINK (version 1.9). The deviation of the observed versus the expected distribution of the P values was represented by the inflation factor λ. The closer λ was to 1, the less stratified the population was considered to be.
Regional plot was generated using LocusZoom (http://locuszoom.sph.umich.edu/). The LD and haplotype block estimations analyses were performed using PLINK and LDBlockShow (1.35 version) . Expression quantitative trait loci (eQTL) and DNA methylation quantitative trait loci (meQTL) were analyzed using Genotype-Tissue Expression (GTEx) project (https://gtexportal.org/) and Pancan-meQTL database (http://bioinfo.life.hust.edu.cn/Pancan-meQTL/) respectively . Functional annotation of GWAS results and chromatin interaction visualization were generated using the online tool functional mapping and annotation (FUMA) v1.3.4c (https://fuma.ctglab.nl/) [11, 12]. Transcription regulation was predicted based on the data of ENCODE database (https://www.encodeproject.org/) . 1000Genomes project (https://www.ncbi.nlm.nih.gov/variation/tools/1000genomes/) and the gnomAD database (http://gnomad-sg.org/) were used to query the racial frequencies of loci in each ethnicity, respectively.
Pearson’s chi-square test or Fisher’s exact test was used to assess the relationship between clinical characteristics and radiotherapy response. The influences of genetic variants on the radiotherapy sensitivity were examined using multivariable logistic regression by calculating odds ratios (ORs) and their 95% confidence intervals (CIs). In this study, all P values were two-sided, and P < 0.05 was considered statistically significant. All these analyses were performed by PLINK (1.9 version), SPSS (20.0 version), and R v0.4.0.3 (www.r-project.org).
As shown in Fig. 1, this was a two-stage GWAS study. The study population initially consisted of 1645 patients with NPC and 479 were excluded due to failure of sample quality control. Ultimately, 1166 patients were included in the two-stage analysis. Of these, 319 patients from the Jiangxi Provincial Cancer Hospital were recruited into the discovery phase cohort. 662 patients from Hunan Cancer Hospital and 185 patients from Jiangxi Cancer Hospital were used to validate the study results. Detailed clinical characteristics are shown in Table 1. Here we analyzed 2 indicators of primary lesions and positive lymph node efficacy after radiotherapy. The characteristics of NPC patients involved in these association analyses are summarized in Additional file 1: Tables S1 and S2. For each phenotype, the association of clinical characteristics of age, gender, BMI, smoking status, EBV infection, clinical stage, and treatment scheme between different subgroups was analyzed. All significant clinical factors were considered as covariates and adjusted in subsequent GWAS analysis.
In the discovery stage, peripheral blood DNA samples from 319 patients were subjected to genome-wide microarray analysis and 688,783 autosomal SNPs were included in the association analysis; in the validation stage, peripheral blood DNA samples from 847 patients were subjected to mass array microarray and TaqMan SNP assays analysis of candidate loci. For each sensitivity phenotype, associations between clinical characteristics such as age, sex, Body Mass Index (BMI), smoking status, EBV infection, clinical stage and treatment regimen were analyzed between different subgroups. All clinical factors that would have an impact were considered as covariates and were adjusted for in subsequent association analyses.
Q-Q plots showed that the expected P-value distributions for all phenotypes were only slightly deviated in the tail region (Additional file 1: Fig. S1). The genomic inflation factor λ for the efficacy of primary lesions and positive lymph nodes after radiotherapy were 1.013 and 1.011 respectively, which were close to 1. In addition, the PCA results showed that there were no deviations from the sample for each efficacy (Additional file 1: Fig. S2). These results suggest that the patient population is ethnically homogeneous and that the observed differences in association analysis are not caused by the internal structure of the population. The Manhattan plot in (Fig. 2a, b) shows the results for both efficacy results of the discovery stage. We further selected significant variants with P < 1 × 10− 6 for genotyping in the validation stage samples. A total of 19 loci entered the validation stage. For the SNP locus rs6414584 screened in the discovery stage, the mass array genotyping was not performed again in the validation stage as it did not pass the conditions of quality control, and the P value and OR value in the validation stage were not given in the table. 18 loci were finally subjected to association analysis in the validation stage. For the efficacy of primary lesions, no SNPs were successfully replicated as being associated with the specific efficacy in the validation cohort (Additional file 1: Table S2). For the efficacy of positive lymph nodes, rs11130424 was successfully validated by showing statistical significance (Additional file 1: Table S3). We further confirmed the effect of rs11130424 on the sensitivity of radiotherapy for NPC with another validation cohort, and statistically significant difference still existed for rs11130424. After combination of the two validation stages, rs11130424 showed a highly statistically significant difference (Table 2). We suggest that rs11130424 may have a correlation with the efficacy of positive lymph nodes after radiotherapy. Thus, it was conducted further analysis.
As previously reported, clinical factors may also influence the sensitivity of radiotherapy. Therefore, we conducted a stratified analysis of rs11130424. Patients were grouped differently according to smoking status, EBV concentration, clinicopathological grade and treatment regimen. Correlations between the rs11130424 locus and these clinical factors were calculated. Clinical characteristics showed some significant influence on the genetic correlation of rs11130424. As shown in Table 3, for smoking status, EBV infection, clinical stage and treatment scheme, the association existed in subgroups of EBV positive, smoking, late stage (III and IV) and patients received both concurrent chemoradiotherapy and induction/adjuvant chemotherapy. This result suggests that clinical factors may influence the results of radiotherapy sensitivity analysis. For example, more than half of patients who were homozygous mutant with late stage (III and IV) were effectively managed after receiving radiotherapy, and the same trend was observed in the subgroup of patients who smoked and received both chemoradiotherapy and induction/adjuvant chemotherapy (Fig. 3).
Based on these results, we concluded that rs11130424 correlates with radiotherapy response in patients of EBV positive, smoking, late stage (III and IV) and receiving both concurrent chemoradiotherapy and induction/adjuvant chemotherapy.
Gene mapping associated with radiotherapy efficacy
To learn the functional consequence of rs11130424, we performed a fine-mapping and functional annotation. As indicated in Fig. 4a, rs11130424 were located in the chromosome 3p21.1. We first performed a linkage disequilibrium analysis using LDBlockShow. Analysis of a block spanning 49.41 kb showed that rs11130424 was not in linkage with any proximal SNPs (Fig. 4b). We further analyzed the eQTL and meQTL profiles of this locus in GTEX and Pancan-meQTL databases. Although rs11130424 was located in the intron region of calcium voltage-gated channel auxiliary subunit alpha2delta 3(CACNA2D3), eQTL analysis did not show that it had an effect with the expression of CACNA2D3 (Additional file 1: Fig. S3). Moreover, the methylation level of the locus differed significantly among genotypes in endometrial carcinoma (P = 4.87 × 10− 6, Fig. 4c). Whether the corresponding methylation site affects the expression of gene needs to be further explored. With the development of three-dimensional genomics, chromatin loops can spatially bring enhancers and promoters closer together, allowing enhancers to bind to promoters and play a role in regulating genes. 3D chromatin interaction indicated that they interacted with coiled-coil domain containing 66 (CCDC66) (Fig. 4d).
Ethnic frequency analysis
Considering that it was in Chinese NPC patients that the rs11130424 locus was identified, we explored its ethnic differences by comparing MAF in five major populations. As shown in Additional file 1: Fig. S4 in the Supplement, the locus showed significant ethnic differences with a high frequency in East Asian populations, further indicating the potential and value of the rs11130424 locus as a biomarker in Chinese populations.
In our study, we explored the genetic variants associated with radiotherapy response using GWAS in a total of 1166 NPC patients. We identified a novel germline mutation of rs11130424 in the CACNA2D3 gene. The G allele carriers were more resistant to radiotherapy. Further stratified analysis showed that there was also a correlation existed in subgroups of EBV positive, smoking, late stage (III and IV) and patients received both concurrent chemoradiotherapy and induction/adjuvant chemotherapy. The efficacy was better in minor allele carriers of rs11130424 than major allele. Rs11130424 is located on the intron region of the CACNA2D3 gene, and the mapped genes were CACNA2D3 and CCDC66.
Radiotherapy is the main treatment for NPC, but there is some individual difference in the efficacy of treatment for patients receiving radiotherapy, the exact mechanism of which is not fully explained. For nasopharyngeal carcinoma, little is known about the genetic variants associated with radiotherapy sensitivity. Yang Z et al. showed an association between variants in autophagy-associated genes (ATG) and radiotherapy efficacy in nasopharyngeal carcinoma. ATG10 rs10514231, rs1864183 and rs4703533 were significantly associated with poorer radiotherapy efficacy at the primary site and positive lymph nodes, while ATG16L2 rs10898880 was positively associated with radiotherapy efficacy . A number of genetic variants have been reported to be associated with radiotherapy response for cancer, including XRCC1, X-ray repair cross complementing 2 (XRCC2), xeroderma pigmentosum group D (XPD), and EGFR [15, 16]. However, these studies have mainly been conducted by candidate gene association strategies, focusing on DNA repair genes. GWAS provides a powerful strategy to discover novel genetic variants. Therefore, we conducted a GWAS study to identify genetic variants associated with radiotherapy sensitivity in NPC patients at a genome-wide level. After a validation stage of association analysis, we eventually found that rs11130424 on chromosome 3p21.1 was associated with response to radiotherapy. This locus is a novel finding and has not been previously reported in nasopharyngeal carcinoma or other cancers.
To date, genome-wide association analysis studies of radiotherapy in head and neck cancers, including nasopharyngeal carcinoma, have focused more on radiotherapy toxicity, and they have revealed a number of genetic variants associated with multiple toxicities, including acute xerostomia , radiotherapy-induced mucositis , and radiotherapy-induced brain injury.  Previously, we also performed a GWAS study of multiple acute radiotherapy toxicities in 1084 patients with nasopharyngeal carcinoma. Rs6711678, rs4848597, rs4848598 and rs2091255 on chromosome 2q14.2 and rs5845477 are novel risk loci for skin toxicity and dysphagia . In this study, we provided the first GWAS study for radiotherapy response in NPC patients. We found that rs11130424 on chromosome 3p21.1 was correlated with radiotherapy sensitivity. Gene mapping and annotation indicated its implicated gene was CACNA2D3 and CCDC66.
CACNA2D3 encodes a subunit of the voltage-dependent calcium channel complex that regulates intracellular Ca2 + influx. It was revealed that CACNA2D3 promoted radiation-induced apoptosis by regulating intracellular calcium ion concentration; and affected the expression of several genes, correlating with the activity of multiple signaling pathways in oesophageal cancer (Fig. 5); in particular, CACNA2D3 inhibited the phosphorylation of PI3K and AKT . Studies have demonstrated that activation of PI3K signaling pathways is closely associated with the development of NPC . As the PI3K/AKT signaling pathway is highly activated in NPC and is closely associated with tumor radiotherapy sensitivity, CACNA2D3 is involved in regulating PI3K, and we speculate that this may be the mechanism by which variants in the CACNA2D3 gene are associated with the occurrence of radiotherapy resistance in patients. In addition, it has been reported that increased levels of Ca2 + in the cytoplasm activate the MAPK cascade , and Ras, Raf, MEK and ERK proteins are key factors in this pathway. Studies have shown that CACNA2D3/CA2 + induces phosphorylation of p38 MAPK in endometrial cancer, which in turn induces apoptosis in tumor cells . CACNA2D3 may influence radiotherapy resistance by regulating the PI3K/AKT and Ras/MARK pathways (Fig. 5).
CCDC66 is a member of the coiled-coil structural domain (CCDC) proteins, which act primarily as oncogenes to regulate tumor proliferation, metastasis, angiogenesis and apoptosis. Although there is no evidence that CCDC66 is directly involved in resistance to radiotherapy, another protein in the family, CCDC86, promotes the proliferation and invasive migration of NPC cells by regulating the activation of the PI3K/AKT pathway by EGFR. EGFR is a key upstream regulator of the PI3K/AKT pathway, and EGFR phosphorylation can indirectly activate PI3K . It has been implicated that the EGFR-PI3K/AKT pathway induces NPC cell senescence and suppresses the tumor stem cell phenotype [26, 27]. We hypothesized that CCDC66 may have a similar regulatory mechanism in nasopharyngeal carcinoma radiotherapy resistance (Fig. 5).
In addition, the gene expression profile and clinical characteristics of GSE48503 were obtained from the GEO database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi) for download. CACNA2D3 and CCDC66 genes were also found among the differentially expressed gene profiles of radiation-resistant and radiation-sensitive NPC cells, which also reaffirmed the reliability of our results (Additional file 1: Fig. S5).
The present study has several limitations. Firstly, the sample size of our study did not have sufficient power to conclusively demonstrate the association between rs11130424 and radiotherapy efficacy, and also reduced the possibility of identifying other genetic variants. Another limitation is that the exact mechanism by which methylation at the rs11130424 locus affects the phenotype is unclear, as it is located in the intron region of CACNA2D3 and the role of methylation in the gene has not been fully established. And in vitro functional studies are needed to validate it. Finally, more newly identified genetic variants can be combined with clinical factors to establish a comprehensive mathematical model for predicting individual efficacy in NPC in the future, thus providing a reference for individualized clinical treatment.
Availability of data and materials
The datasets and code that support the findings of this study are available from the corresponding author on reasonable request.
Genome-wide association analysis
X-ray repair cross complementing 1
Epidermal growth factor receptor
Intensity Modulated Radiation Therapy
Total tumor volume
Clinical target volume
Planned target volume
Response Evaluation Criteria in Solid Tumors
Principal component analysis
Expression quantitative trait loci
DNA methylation quantitative trait loci
Genotype-Tissue Expression project
Functional mapping and annotation
Body Mass Index
Calcium voltage-gated channel auxiliary subunit alpha2delta 3
Coiled-coil domain containing 66
X-ray repair cross complementing 2
Xeroderma pigmentosum group D
Coiled-coil structural domain
Wong KCW, Hui EP, Lo KW, Lam WKJ, Johnson D, Li L, et al. Nasopharyngeal carcinoma: an evolving paradigm. Nat Rev Clin Oncol. 2021;18(11):679–95.
Chen YP, Chan ATC, Le QT, Blanchard P, Sun Y, Ma J. Nasopharyngeal carcinoma. Lancet. 2019;394(10192):64–80.
Guigay J, Temam S, Bourhis J, Pignon JP, Armand JP. Nasopharyngeal carcinoma and therapeutic management: the place of chemotherapy. Ann Oncol. 2006;17(Suppl 10):304–7. x.
Lee AW, Ma BB, Ng WT, Chan AT. Management of nasopharyngeal carcinoma: current practice and future perspective. J Clin Oncol. 2015;33(29):3356–64.
Kerns SL, Ostrer H, Rosenstein BS. Radiogenomics: using genetics to identify cancer patients at risk for development of adverse effects following radiotherapy. Cancer Discov. 2014;4(2):155–65.
Yoon HH, Catalano PJ, Murphy KM, Skaar TC, Philips S, Powell M, et al. Genetic variation in DNA-repair pathways and response to radiochemotherapy in esophageal adenocarcinoma: a retrospective cohort study of the Eastern Cooperative Oncology Group. BMC Cancer. 2011;11:176.
Jin G, Fan XM, Li KX, Niu SH, Zhang QY. The association between epidermal growth factor receptor single nucleotide polymorphisms and Radiochemotherapy Response in Cervical Cancer. Pathol Oncol Res. 2020;26(2):1255–61.
Li Q, Liang Y, Liu Z, Yu C. Associations of GWAS-Identified risk loci with progression, efficacy and toxicity of Radiotherapy of Head and Neck squamous cell carcinoma treated with Radiotherapy. Pharmgenomics Pers Med. 2021;14:1205–10.
Dong SS, He WM, Ji JJ, Zhang C, Guo Y, Yang TL. LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files. Brief Bioinform. 2021;22:4.
Gong J, Wan H, Mei S, Ruan H, Zhang Z, Liu C, et al. Pancan-meQTL: a database to systematically evaluate the effects of genetic variants on methylation in human cancer. Nucleic Acids Res. 2019;47(D1):D1066-d1072.
Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26(18):2336–7.
Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. 2017;8(1):1826.
Davis CA, Hitz BC, Sloan CA, Chan ET, Davidson JM, Gabdank I, et al. The encyclopedia of DNA elements (ENCODE): data portal update. Nucleic Acids Res. 2018;46(D1):D794-d801.
Yang Z, Liu Z. Potentially functional variants of autophagy-related genes are associated with the efficacy and toxicity of radiotherapy in patients with nasopharyngeal carcinoma. Mol Genet Genomic Med. 2019;7(12):e1030.
Yang Z, Liu Z. Potential functional variants in DNA repair genes are associated with efficacy and toxicity of radiotherapy in patients with non-small-cell lung cancer. J Oncol. 2020;2020:3132786.
Butkiewicz D, Krześniak M, Gdowicz-Kłosok A, Giglok M, Marszałek-Zeńczak M, Suwiński R. Polymorphisms in EGFR gene predict clinical outcome in unresectable non-small cell lung cancer treated with radiotherapy and platinum-based chemoradiotherapy. Int J Mol Sci. 2021;22:11.
Naderi E, Crijns APG, Steenbakkers R, van den Hoek JGM, Boezen HM, Alizadeh BZ, et al. A two-stage genome-wide association study of radiation-induced acute toxicity in head and neck cancer. J Transl Med. 2021;19(1):481.
Schack LMH, Naderi E, Fachal L, Dorling L, Luccarini C, Dunning AM, et al. A genome-wide association study of radiotherapy induced toxicity in head and neck cancer patients identifies a susceptibility locus associated with mucositis. Br J Cancer. 2022;126(7):1082–90.
Wang TM, Shen GP, Chen MY, Zhang JB, Sun Y, He J, et al. Genome-wide association study of susceptibility loci for radiation-induced brain injury. J Natl Cancer Inst. 2019;111(6):620–8.
Wang Y, Xiao F, Zhao Y, Mao CX, Yu LL, Wang LY, et al. A two-stage genome-wide association study to identify novel genetic loci associated with acute radiotherapy toxicity in nasopharyngeal carcinoma. Mol Cancer. 2022;21(1):169.
Nie C, Qin X, Li X, Tian B, Zhao Y, Jin Y, et al. CACNA2D3 enhances the chemosensitivity of esophageal squamous cell carcinoma to cisplatin via inducing Ca(2+)-mediated apoptosis and suppressing PI3K/Akt pathways. Front Oncol. 2019;9:185.
Lin DC, Meng X, Hazawa M, Nagata Y, Varela AM, Xu L, et al. The genomic landscape of nasopharyngeal carcinoma. Nat Genet. 2014;46(8):866–71.
Chen KC, Kao PH, Lin SR, Chang LS. Upregulation of Fas and FasL in Taiwan cobra phospholipase A2-treated human neuroblastoma SK-N-SH cells through ROS- and Ca2+-mediated p38 MAPK activation. J Cell Biochem. 2009;106(1):93–102.
Kong X, Li M, Shao K, Yang Y, Wang Q, Cai M. Progesterone induces cell apoptosis via the CACNA2D3/Ca2+/p38 MAPK pathway in endometrial cancer. Oncol Rep. 2020;43(1):121–32.
Martelli AM, Evangelisti C, Chiarini F, Grimaldi C, McCubrey JA. The emerging role of the phosphatidylinositol 3-kinase/ akt/mammalian target of rapamycin signaling network in cancer stem cell biology. Cancers (Basel). 2010;2(3):1576–96.
Chan KC, Ting CM, Chan PS, Lo MC, Lo KW, Curry JE, et al. A novel Hsp90 inhibitor AT13387 induces senescence in EBV-positive nasopharyngeal carcinoma cells and suppresses tumor formation. Mol Cancer. 2013;12(1):128.
Ma L, Zhang G, Miao XB, Deng XB, Wu Y, Liu Y, et al. Cancer stem-like cell properties are regulated by EGFR/AKT/β-catenin signaling and preferentially inhibited by gefitinib in nasopharyngeal carcinoma. Febs j. 2013;280(9):2027–41.
We thank all patients who participated in the study and the participating clinic staff for their contribution to data collection.
This work was supported by the National Natural Science Foundation of China (82073943, 81773823 and 81974511), the Hunan Cancer Hospital Climb Plan (YF2020011), the Scientific and Technological Project of Changsha (kq2004147), Science Foundation of Hunan Province (NO. 2022JJ30363), the Wisdom Accumulation and Talent Cultivation Project of the Third Xiangya Hospital of Central South University (YX202110), and The Open Fund for Scientific Research of NHC Key Laboratory of Personalized Diagnosis and Treatment of Nasopharyngeal Carcinoma (2021NPCK01).
Ethics approval and consent to participate
The study protocol was approved by the Ethics Committee of Jiangxi Province Tumor Hospital and Hunan Cancer Hospital. All subjects were provided with written informed consents.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Quantile–quantile (QQ) plot of observed association P-values (y-axis) against expected P-values (x-axis) in the discovery stage. X-axis represents –log10 expected P-values and Y-axis represents –log10 observed P-values. Figure S2. Distribution of samples according to PCA analysis in discovery stage. Figure 3. eQTL violin plot from GTEx displaying significant association between rs11130424 and the expression of CACNA2D3 within the whole blood. Figure S4. The MAF of rs11130424 in different ethnic populations. SAS: South Asian, EUR: Europe, EAS: East Asian, AMR: American, AFR: African. Figure S5. Volcano plot showed differentially expressed genes in the sensitive and resistant groups. Table S1. Characteristics of NPC patients involved in primary lesion efficacy association analysis. Table S2 Characteristics of NPC patients involved in positive lymph node efficacy association analysis. Table S3. SNPs significantly associated with primary lesion efficacy in nasopharyngeal carcinoma patients after radiotherapy. Table S4. SNPs significantly associated with positive lymph node efficacy in nasopharyngeal carcinoma patients after radiotherapy.
About this article
Cite this article
Yu, LL., Hu, BW., Huang, HX. et al. A two-stage genome-wide association study identifies novel germline genetic variations in CACNA2D3 associated with radiotherapy response in nasopharyngeal carcinoma. J Transl Med 21, 11 (2023). https://doi.org/10.1186/s12967-022-03819-4