Abstract
Chemoradiation-induced hearing loss (CRIHL) is one of the most devasting side effects for nasopharyngeal carcinoma (NPC) patients, which seriously affects survivors’ long-term quality of life. However, few studies have comprehensively characterized the risk factors for CRIHL. In this study, we found that age at diagnosis, tumor stage, and concurrent cisplatin dose were positively associated with chemoradiation-induced hearing loss. We performed a genome-wide association study (GWAS) in 777 NPC patients and identified rs1050851 (within the exon 2 of NFKBIA), a variant with a high deleteriousness score, to be significantly associated with hearing loss risk (HR = 5.46, 95% CI 2.93–10.18, P = 9.51 × 10–08). The risk genotype of rs1050851 was associated with higher NFKBIA expression, which was correlated with lower cellular tolerance to cisplatin. According to permutation-based enrichment analysis, the variants mapping to 149 hereditary deafness genes were significantly enriched among GWAS top signals, which indicated the genetic similarity between hereditary deafness and CRIHL. Pathway analysis suggested that synaptic signaling was involved in the development of CRIHL. Additionally, the risk score integrating genetic and clinical factors can predict the risk of hearing loss with a relatively good performance in the test set. Collectively, this study shed new light on the etiology of chemoradiation-induced hearing loss, which facilitates high-risk individuals’ identification for personalized prevention and treatment.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
Platinum compounds are the most widely used chemotherapy agents for cancer worldwide, but are also known for their serious side effects of ototoxicity. In addition, radiotherapy to head and neck tumors, such as nasopharyngeal carcinoma (NPC), often inevitably damages the adjacent organs of auditory, which may further increase the risk of hearing loss (HL). Previous studies reported that both platinum and radiation could contribute to DNA damage to auditory cells (Mujica-Mota et al. 2014), leading to apoptosis and death of hair cells, stria vascularis endothelial cells, and spiral ganglion neurons (Landier 2016; Shi et al. 2021). Platinum-based concurrent chemoradiotherapy has been recommended as the standard treatment modality for advanced NPC (Pfister et al. 2020), a malignancy particularly prevalent in the east and southeast Asia, especially in Southern China (Chen et al. 2019). Therefore, the side effect of chemoradiation-induced hearing loss (CRIHL) is particularly devasting for NPC patients (Low et al. 2006), with 21–34% developing hearing impairment (Wu et al. 2014; Tang et al. 2018; Zhang et al. 2019b). However, since the lack of effective treatment, affected patients suffered from a progressive and persistent process of hearing loss, often leading to dramatic declines in their quality of life.
Several clinical factors, including age (Kwong et al. 1996; Wang et al. 2015), cochlea radiation dose, and cisplatin dose (Chen et al. 2006; Chan et al. 2009; Wang et al. 2015), have been reported to be associated with chemotherapy or radiotherapy-related ototoxicity. Nevertheless, these clinical parameters could not fully account for inter-patient variability, suggesting there remain other potential factors playing roles in individual predisposition to treatment-induced hearing loss. Genome-wide association studies (GWASs) have reported that genetic variants were associated with cisplatin or radiation-induced ototoxicity in pediatric cancer survivors (mainly medulloblastoma, osteosarcoma, and acute lymphoblastic leukemia) or adult testicular patients (Xu et al. 2015; Wheeler et al. 2017; Meijer et al. 2021; Trendowski et al. 2021), indicating the underlying role of genetic variants in cisplatin or radiation-induced hearing loss. However, previous GWASs were mainly conducted in European populations of pediatric patients and those with certain cancer types, only focusing on the risk of hearing loss induced by cisplatin or radiotherapy alone. Researches to investigate the genetic backgrounds for hearing loss resulting from chemotherapy in combination with radiotherapy in other ancestries were limited. Especially for NPC, which is treated with a combined modality of chemoradiotherapy and has an extremely unbalanced ethnical distribution, the etiology of CRIHL has never been fully studied yet.
To comprehensively characterize the risk factors for CRIHL, we evaluated the influence of clinical factors and identified the genetic factors in NPC patients who underwent combined chemoradiotherapy. We think the current study would be helpful to identify high-risk individuals for personalized prevention and treatment for CRIHL.
Materials and methods
Patient recruitment and follow-up
A total of 1278 histologically confirmed NPC patients, without distant metastatic at diagnosis, were recruited from Sun Yat-sen University Cancer Center (SYSUCC), between 2009 and 2015. Of these 1278 patients, 21 patients younger than 18 years old at diagnosis, 265 patients with self-reported hearing loss before cancer treatment, and 215 patients without concurrent chemoradiotherapy were excluded. Then, a total of 777 eligible patients were included in the present study. All 777 patients in this study underwent concurrent chemoradiotherapy, which has been recommended as the standard treatment protocol for advanced NPC based on the National Comprehensive Cancer Network guidelines (Pfister et al. 2020). And all patients underwent a standardized intensity-modulated radiotherapy (IMRT) with a prescribed total dose of 66–72 Gy in 28–33 fractions (Wen et al. 2021). Most of patients (704/777) received platinum-based concurrent chemoradiotherapy with/without induction chemotherapy and/or adjuvant chemotherapy. Induction chemotherapy mainly consisted of cisplatin (80–100 mg/m2) plus docetaxel and/or fluorouracil every 3 weeks for 2–3 cycles (Yao et al. 2017). Concurrent chemotherapy was mainly cisplatin administered weekly (30–40 mg/m2, for 4–6 cycles) or triweekly (80–100 mg/m2, for 2–3 cycles) during radiation therapy (Yao et al. 2017).
After completion of treatment, patients were followed up every 3 months during the first 3 years and then every 6 months. In addition to completing clinical examination and laboratory testing, patient-reported adverse events were recorded during outpatient clinic follow-up visits. Patients reporting the event of hearing loss in outpatient clinic follow-up records documented by physicians were defined as cases with hearing loss, whereas those without reported hearing loss were classified as the non-hearing loss group. The details for the collection of other clinical data are described in the Supplementary file 1. The time to hearing loss event was defined as the interval between the initiation of radiotherapy or chemotherapy and the time when the event of hearing loss was first documented. The follow-up time was measured from the start of cancer treatment to the date of the last outpatient clinic follow-up. The Institutional Review Board of Sun Yat-sen University Cancer Center approved this study. A brief overview of study design is presented in Fig. S1.
Genotyping, imputation and quality control
DNA was isolated from peripheral blood samples using a commercial DNA extraction kit (Qiagen). Genotyping was performed on the Illumina Infinium Global Screening Array. A standard quality control protocol for genotypes was implemented in PLINK v1.9. The detailed information about genotyping and quality control was described in our previous study (He et al. 2022). Individuals with call rate < 95%, gender discrepancy, extreme heterozygosity rate (inbreeding coefficient F > 6 SDs), or pairwise identify by descent (IBD) > 0.25 were removed from the analysis. Besides, genetic principal component analysis (PCA) was applied to identify population stratification outliers (> 6 SDs on any one of the top 10 PCs using five iterations) based on EIGENSTRAT (Price et al. 2006). SNPs on autosomal chromosomes were included, but those with call rate < 95%, minor allele frequency (MAF) < 0.01 or deviating from Hardy-Weinberg equilibrium (HWE) (P < 1.0 × 10–12) were removed. Imputation was performed with SHAPEIT2 and IMPUTE2, using 1000 Genome Phase III integrated variant set of the entire population as a reference panel. Further quality control was conducted and SNPs with INFO < 0.8, call rate < 95%, or MAF < 0.01 were excluded. Finally, a total of 777 subjects and 5,248,879 variants were retained for subsequent analysis.
Genome-wide association study
We performed a genome-wide association analysis in 777 patients to identify potential genetic variants by assessing the associations between 5,248,879 SNPs and risk of hearing loss, adjusting for age at diagnosis, T stage, concurrent cisplatin dose, and the first two genetic principal components. Manhattan plot and quantile–quantile plot were generated for showing the overall significance of GWAS. The genomic inflation factor (λGC) was calculated to evaluate the deviation of the observed versus the expected distribution of P-values. Kaplan–Meier curve was drawn to illustrate the relationships between SNP genotypes and CRIHL after adjusting for the same covariates aforementioned. Additionally, since cisplatin is the most ototoxic platinum compound (Romano et al. 2020) and the main agent used in NPC chemotherapy, we conducted further subgroup analysis in 516 patients treated with cisplatin-based chemotherapy to explore the robustness of associations for the SNPs.
Permutation-based enrichment analysis
We performed permutation tests to explore the genetic similarity between hereditary hearing loss and chemoradiation-induced hearing loss by examining the enrichment of the genetic variants mapping to the reported hereditary hearing loss genes among GWAS top signals. 149 hereditary deafness genes have been reported to be associated with autosomal non-syndromic and/or syndromic hearing loss from the Hereditary Hearing Loss Homepage (https://hereditaryhearingloss.org/). Since approximately 90% of linkage disequilibrium (LD) blocks in the Asian population were observed within the 50 kb gene boundary extension (Wang et al. 2007), we selected 50 kb as the gene boundary for the permutation-based enrichment analysis, which is also referred to the previous studies (Hoffmann et al. 2016; Ghisdal et al. 2017; Wheeler et al. 2017). Then, 56,040 genetic variants within 50 kb of these 149 genes were included for further analysis. In each shuffle, the patient phenotypes, including age at diagnosis, T stage, concurrent cisplatin dose, and ototoxicity (as a time-to-event variable), were randomly permuted relative to the patient genotypes and the association analyses were repeated on the permuted resamples. We counted the number of deafness gene SNPs with a range of significance thresholds from P < 5.0 × 10–02 to P < 5.0 × 10–06 in association test for each replicate to examine the robustness of enrichment analysis. 1000 permutation shuffles were conducted to construct the null distribution. The distribution of the number of significant genetic variants within 50 kb of deafness genes was displayed in a histogram, compared with the observed number in the current study. The empirical P-value was obtained by calculating the probability of detecting the expected number of significant SNPs larger than or equal to the observed number.
Functional annotations
We conducted functional annotations for variants in moderate to high LD (r2 ≥ 0.6) with the lead variants by FUMA (https://fuma.ctglab.nl/). The LD between genetic variants was calculated on the basis of 1000 Genomes Phase III East Asian population. We annotated SNPs by performing ANNOVAR gene-based annotation. The details of functional prediction scores, chromatin marks, and expression quantitative trait loci (eQTL) analysis are shown in the Supplementary file 1.
Gene expression and cellular sensitivity to cisplatin or radiation
To determine whether the expression of susceptible genes exerts effects on chemosensitivity, we examined the relationship between gene expression and cellular sensitivity to cisplatin in the central nervous system (CNS) cancer cell lines (mainly glioma cell lines), as well as other cancer cell lines across 11 histology types (bladder/urinary tract, bone, bowel, breast, kidney, lung, ovary/fallopian tube, pancreas, peripheral nervous system, skin and soft tissue), which was referred to in previous studies (Wheeler et al. 2017; El Charif et al. 2019) since there is no available relevant data on auditory cell lines. Baseline gene expression data in different cell lines was generated by the Cancer Cell Line Encyclopedia (CCLE) project (Barretina et al. 2012). Cisplatin sensitivity, measured as the concentration of cisplatin inhibiting half of cellular growth (IC50), was derived from the Genomics of Drug Sensitivity in Cancer (GDSC) database (Yang et al. 2013). We obtained data on gene expression and cisplatin sensitivity from the curated pharmacogenomic datasets by R package “PharmacoGx” (Smirnov et al. 2016).
Similarly, we examined the correlations of gene expression with radiosensitivity in the cell line types mentioned above (Trendowski et al. 2021). Radiosensitivity was measured by the surviving fraction after 2 Gy radiation (SF2) from the 9-days viability assay in tumor cell lines (Yard et al. 2016). We obtained data on gene expression and radiosensitivity from the curated radiogenomic datasets by R package “RadioGx” (Manem et al. 2019). Associations of gene expression with cellular sensitivity to cisplatin or radiation were estimated by Spearman’s rank correlation test in R 4.1.1.
Pathway analysis
We conducted pathway enrichment analysis to explore the biological processes involved in the development of CRIHL by R package “clusterProfiler” (Yu et al. 2012). Genes, which had at least one SNP with P < 0.001 in the GWAS analysis within their 20 kb upstream/downstream regions, were kept, as the previous study described (Wang et al. 2019). We included 799 candidate genes for pathway enrichment analysis based on Gene Ontology (GO) database. Gene sets with a size less than 10 or more than 500 were excluded from the analysis. Benjamini and Hochberg approach was used to correct for multiple testing.
Construction of a risk score model
By combining genetic and clinical factors, we constructed a risk score model to estimate the 5-years risk of CRIHL via Cox regression model. Eight SNPs, including seven lead SNPs with P < 1.0 × 10–06 in the GWAS (rs1050851, rs117098517, rs1485149, rs17010289, rs79938362, rs73335760 and rs201061882) and one lead SNP with P < 1.0 × 10–05 within the hereditary deafness gene regions (rs2275994), were included in the risk score. The clinical factors of age at diagnosis, tumor stage, and concurrent cisplatin dose were also included in the risk score. All 777 subjects in this study were randomly dichotomized into a training set (n = 389) and a test set (n = 388). The risk score model was developed based on the samples in the training set and then validated in the test set. We generated the risk score by summing the values weighted by the coefficients of the genetic and clinical variables. The risk score for each patient in this study was calculated, and the median risk score in the training set was regarded as a cutoff value to divide subjects into low or high-risk score groups.
Statistical analysis
Univariate Cox regression analysis was carried out to test the association of clinical factors with CRIHL. Variables with a two-tailed P-value of less than 0.05 in univariate analysis were adjusted in multivariate Cox regression analysis. The statistically significant (P < 0.05) clinical factors determined by multivariate analysis were included as covariates in genome-wide association analysis. The optimal cutoff points of continuous variables were defined by receiver operating characteristic (ROC) curve analyses, and categorical variables were classified based on clinical relevance. The GWAS analysis was conducted by Cox proportional hazard model under an additive assumption using the R package “gwasurvivr” (Rizvi et al. 2019). Stratified analysis was conducted to investigate the stability property of association results of the lead variants. All subjects were dichotomized based on age at diagnosis, sex, T stage, pre-treatment EBV DNA levels, treatment modality, primary tumor dose and concurrent cisplatin dose, respectively. In each subgroup, hazard ratio (HR), 95% confidence interval (CI) and P-value were estimated via Cox regression model with adjustment for age at diagnosis, T stage and concurrent cisplatin dose (except for itself in the subgroups defined by these variables). All statistical tests were two-sided and performed in R 4.1.1.
Results
Patient characteristics
777 eligible NPC patients were recruited and followed up in this study. With a median follow-up time of 70.7 months (IQR 52.7–85.5), 132 (17.0%) patients developed chemoradiation-induced hearing loss. The median time interval from the initiation of cancer treatment to the event of hearing loss occurrence was 18.4 months. The clinical characteristics were shown in Table 1. By univariate Cox regression analyses, age at diagnosis, T stage, induction cisplatin dose, and concurrent cisplatin dose were significantly correlated with CRIHL. By multivariate analyses to identify independent clinical risk factors, we found that elder age (for ≥ 50 years vs. < 50 years at diagnosis: adjusted HR = 1.72, 95% CI 1.19–2.49, P = 0.004), advanced tumor stage (for T3–T4 vs. T1–T2: adjusted HR = 1.66, 95% CI 1.12–2.46, P = 0.012), and higher concurrent cisplatin dose (for ≥ 155 mg/m2 vs. < 155 mg/m2: adjusted HR = 1.46, 95% CI 1.01–2.12, P = 0.046) were significantly associated with higher risk of hearing loss. However, the other factors of sex, smoking, drinking, pre-treatment EBV DNA loads, radiation dose for the primary tumor, induction cisplatin dose, hypertension, and diabetes were not significantly associated with hearing loss risk.
GWAS analysis for chemoradiation-induced hearing loss
The GWAS was conducted in 777 patients to assess the associations between 5,248,879 SNPs and the risk of hearing loss by the Cox proportional hazard model under an additive assumption. The quantile–quantile plot showed a good match between the distribution of the observed P-values and the expected P-values by chance (λGC = 1.01), and an apparent deviation within the tail of the distribution, indicating true associations (Fig. S2). GWAS identified 30 loci surpassing the suggestive threshold (P < 1.0 × 10–05; Table S1). rs1050851 was the most significant SNP associated with the risk of CRIHL (HR = 5.46, 95% CI 2.93–10.18, P = 9.51 × 10–08; Fig. 1A and Table S2). Since cisplatin is the most ototoxic platinum compound (Romano et al. 2020) and the main agent used in NPC chemotherapy, we conducted a subgroup analysis in 516 patients who were treated with cisplatin-based chemotherapy and the clinical characteristics were shown in Table S3. We found rs1050851 exceeded the significant threshold of P < 5.0 × 10–08 (HR = 9.16, 95% CI 4.76–17.64, P = 3.32 × 10–11; Table S4) in the subgroup patients who received cisplatin-based chemotherapy, indicating the robustness of the association. Patients with the risk allele (A) at rs1050851 were more susceptible to hearing loss in comparison with those carrying the non-risk allele (G) (Fig. 1C). To be specific, the incidence of hearing loss among patients carrying AG genotype at rs1050851 was 55.0%, which was much higher than that of 15.8% among patients carrying GG genotype (no AA genotype was observed).rs1050851 is a synonymous coding variant within the exon 2 of NFKBIA gene (Fig. 1B). To measure the deleteriousness of this SNP, we annotated with functional prediction scores and found it has a very high CADD score of 17.10, indicating a deleterious mutation (Table S5). eQTL analysis showed that the risk allele at rs1050851 was correlated with higher expression of NFKBIA in human brain tissues (such as the cerebellar cortex; Fig. 1D and Table S6), as well as in blood tissues (Table S7). Similar eQTL associations were also seen between genotype of rs3138054 (r2 = 0.78 with rs1050851) and NFKBIA expression (Fig. S3, Tables S6 and S7). Additionally, both rs1050851 and rs3138054 were located within the regions containing enhancer or promoter histone markers (eg, H3K27ac, H3K4me1, and H3K4me3) in multiple tissues, including ganglion eminence derived neurospheres, brain cingulate gyrus and lymphoblastoid cells, as well as within the DNase I hypersensitive site (DHS) peaks in lymphoblastoid cells (Fig. 1E). These results suggested that rs1050851 and rs3138054 might modulate the expression of NFKBIA through affecting the promoter or enhancer activity, or chromatin accessibility, to influence the risk of CRIHL.
By cisplatin and radiation cytotoxicity analysis, we found that higher expression of NFKBIA gene was significantly associated with lower cisplatin IC50 (that is, lower tolerance to cisplatin) in CNS cancer cell lines (ρ = − 0.416, P = 0.005; Fig. 1F), while the expression of NFKBIA gene was not associated with radiosensitivity in the cell lines (ρ = 0.118, P = 0.500), indicative of a harmful function of NFKBIA against cisplatin-induced injury. NFKBIA gene encodes a member of NF-κB inhibitor family, IκBα protein, which is involved in inflammatory responses in the CNS (O'Neill and Kaltschmidt 1997). Interestingly, we found that the significant association of cisplatin sensitivity with NFKBIA expression was specific to cancer cell lines from CNS tissues, not from other tissue types (Fig. 1G).
By and large, carriers with the minor allele (A) at rs1050851 would have higher NFKBIA expression and lower cisplatin tolerance, being prone to cisplatin-induced injury, which may partially explain the increased hearing loss risk for patients carrying the risk allele at rs1050851 in the GWAS analysis.
Enrichment analysis of hereditary hearing loss genes
To delineate a comprehensive overview of genetic etiology for hearing loss in cancer patients undergoing chemoradiation, we further evaluated whether the known 149 autosomal hereditary deafness genes were involved in CRIHL. The quantile–quantile plot indicated that variants in known deafness genes were more likely to correlate with CRIHL than expected by chance (Fig. 2A). Of 149 autosomal hearing loss genes tested, 118 genes had at least one SNP with P-value < 0.05 (Table S8). Interestingly, of 118 deafness genes, several genes have been reported to be associated with auditory neuropathy, a hearing dysfunction characterized by impaired transmission of signal to the auditory nerve by the presynaptic inner hair cells (Starr et al. 1996), such as PJVK (Delmaghani et al. 2006), NARS2 (Simon et al. 2015), SLC17A8 (Ruel et al. 2008), DIAPH3 (Schoen et al. 2010) and DIAPH1 (Wu et al. 2020).
The permutation-based enrichment analysis showed that SNPs within 50 kb of the 149 genes were significantly enriched among the GWAS top signals (P < 1.0 × 10–04) for CRIHL with an empirical P-value of 0.034 (Fig. 2B). The results of enrichment analysis were consistent across a range of P-value thresholds for GWAS from P < 5.0 × 10–05 to P < 5.0 × 10–06 and remained significant (for threshold: PGWAS < 5.0 × 10–05, empirical P = 0.016; for threshold: PGWAS < 1.0 × 10–05, empirical P = 0.033; for threshold: PGWAS < 5.0 × 10–06, empirical P = 0.027; Fig. S4). But the results showed a non-significant enrichment (empirical P > 0.05) for the P-value thresholds from P < 5.0 × 10–02 to P < 5.0 × 10–04.rs2275994, located at 18 kb downstream of GRHL2, was the most significant SNP within the deafness loci (HR = 3.85, 95% CI 2.24–6.61, P = 1.10 × 10–06; Fig. 2C and Table S9). Patients carrying the risk allele (T) at rs2275994 were more likely to develop hearing loss than those with the non-risk allele (C) (Fig. 2D). This variant has a relatively low CADD score of 7.12 (Table S10). eQTL analysis in human brain tissues (such as the temporal cortex) showed that the risk allele at rs2275994 was correlated with lower expression of GRHL2 (Fig. 2E and Table S11). Similar eQTL associations between GRHL2 expression and another two variants (rs10091039 and rs16868134) in moderate LD with rs2275994 (r2 = 0.69 and 0.66, respectively) were also found in multiple brain tissues (Fig. S5 and Table S11). Additionally, rs2275994 resided in the DHS peaks in lymphoblastoid cells, and rs10091039 was located within the DHS peaks in fetal brain tissues, and rs10091039 and rs16868134 were located within the regions embracing enhancer histone marks (H3K4me1) in neuron cultured cells (Fig. 2F). These results suggested that rs2275994, rs10091039 and rs16868134 might modulate the expression of GRHL2 by affecting chromatin accessibility, or enhancer activity to influence the risk of CRIHL.
Using cisplatin and radiation cytotoxicity analysis, we found that lower expression levels of GRHL2 were significantly correlated with lower radiation SF2 (that is, lower tolerance to radiation) in CNS cancer cell lines (ρ = 0.341, P = 0.045; Fig. 2G), while the expression of GRHL2 was not significantly associated with cisplatin sensitivity in the cell lines (ρ = 0.045, P = 0.768), indicative of a protective effect of GRHL2 on radiation-related damage. Similarly, the significant association between radiosensitivity and GRHL2 expression was specifically exited in tumor cell lines of CNS origin instead of other tissue origins (Fig. 2H).
In short, carriers with the minor allele (T) at rs2275994 would have lower GRHL2 expression and lower radiation tolerance, being susceptible to radiation-induced damage, which may partially explain the increased hearing loss risk for patients carrying the risk allele at rs2275994 in the GWAS analysis.
Stratified analysis
To examine the robustness of the association of ototoxicity risk with the lead variants (rs1050851 and rs2275994) identified above, we conducted stratified analysis by various clinical subgroups, such as by different groups of age (< 50 and ≥ 50 years), sex (male and female), T stage (T1–T2 and T3–T4), pre-treatment EBV DNA levels (< 4000 and ≥ 4000 copies/ml), treatment modality (CCRT alone and CCRT + IC/AC), primary tumor dose (< 70.2 and ≥ 70.2 Gy) and concurrent cisplatin dose (< 155 and ≥ 155 mg/m2). We found the associations between the risk of hearing loss and rs1050851 or rs2275994 were relatively stable and still remained significant within most of clinical subgroups (Fig. S6).
Pathway analysis
To further evaluate the potential biological mechanisms responsible for the development of CRIHL, we conducted pathway enrichment analysis using 799 candidate genes that had at least one SNP with P < 0.001 in the GWAS analysis within their 20-kb upstream/downstream regions. Based on GO enrichment analysis, candidate genes were enriched in the gene sets such as “synaptic membrane”, “neuron to neuron synapse” and “modulation of chemical synaptic transmission” (Fig. 3). This result indicated that synaptic signaling and neuronal connectivity might participate in the development of CRIHL.
Predictive risk score model for chemoradiation-induced hearing loss
To predict the risk of CRIHL, we developed risk score models using clinical factors only, genetic factors only, as well as integrating both clinical and genetic factors in a training set (n = 389), and validated their predictive values in a test set (n = 388). The clinical factors consisted of age at diagnosis, tumor stage, and concurrent cisplatin dose. The genetic factors included rs1050851, rs117098517, rs1485149, rs17010289, rs79938362, rs73335760, rs201061882 and rs2275994. The area under the curve (AUC) of the risk score models derived from clinical model, genetic model, and combined model in the training set were 0.63, 0.76 and 0.80, respectively (Fig. 4A). The risk score integrating genetic and clinical factors was well-replicated in the test set (AUC = 0.78; Fig. 4B), as well as in all samples (n = 777, AUC = 0.79; Fig. 4C).
To investigate whether the risk score integrating genetic and clinical variables can predict the development of CRIHL, we regarded the median risk score of 0.96 in the training set as the cutoff point to split all subjects into low or high-risk score groups. We found that the risk of CRIHL in patients with high-risk score was significantly higher than those with low-risk score (training set: HR = 5.27, P = 3.85 × 10–09; test set: HR = 4.54, P = 2.84 × 10–07; all samples: HR = 4.88, P = 6.88 × 10–15; Fig. 4D–F). These results suggested that the risk score combining the genetic and clinical variables could predict the development of CRIHL with a relatively good performance, which could facilitate better risk stratification.
Discussion
Chemotherapy or radiotherapy-related hearing loss has been considered one of the most common and important survivorship concerns for cancer patients. We conducted this study to explore the clinical and genetic factors on the risk of CRIHL, and provided evidence that rs1050851 within NFKBIA locus was associated with CRIHL by bioinformatic functional analyses. Moreover, we found a genetic similarity between hereditary deafness and CRIHL. These findings have extended knowledge of treatment-related ototoxicity, providing opportunities to deliver more insights into the etiology of CRIHL, and new clues for preventing or alleviating the side effect of cancer treatment.
Several clinical risk factors have been reported to be associated with cisplatin or radiation-induced hearing loss (Bhandare et al. 2010; Trendowski et al. 2019), and our study further confirmed that age at diagnosis, tumor stage, and concurrent cisplatin dose were independent risk factors for CRIHL. In this study, we found that survivors who were older at diagnosis and with later T stage were more likely to experience ototoxicity after treatment, which was in accord with previous studies on NPC patients (Kwong et al. 1996; Wang et al. 2015). The survivors who were treated with a higher concurrent cisplatin dose would have a significantly higher risk of ototoxicity, while there was no significant association between induction cisplatin dose and the ototoxicity risk, demonstrating that concurrent cisplatin exerted more effects than induction or adjuvant cisplatin on hearing loss, which is also consistent with the findings from another study on hearing loss for NPC patients (Chan et al. 2009).
In addition to clinical variables, GWAS identified potential susceptible genetic variants for CRIHL. Our in-silico analysis showed that rs1050851, within the exon 2 of NFKBIA, was associated with cisplatin-induced hearing loss. NFKBIA encoding protein IκBα binds with NF-κB as an inactivated complex under resting conditions (Shih et al. 2011). When exposed to external harmful stimulation, IκBα is degraded and the activated NF-κB is released, which is involved in the inflammatory response (Vallabhapurapu and Karin 2009; Shih et al. 2011). The role of NF-κB pathway in hearing loss is complicated and multi-faceted. For example, the deficiency of NF-κB is associated with auditory nerve degeneration and increased noise-induced hearing loss in mice (Lang et al. 2006). In humans, the molecular basis of ROR1 missense mutation-caused inner ear anomalies and deafness are linked to NF-κB deficiency (Diaz-Horta et al. 2016), and this pathway also protects cochlear hair cells from aminoglycoside-induced ototoxicity (Jiang et al. 2005). However, activation of NF-κB pathway has been linked to cochlear inflammation in rats with noise-induced hearing loss (Zhang et al. 2019a). Previous studies have indicated that treatment with cisplatin in cultured auditory cells (HEI-OC1) would decrease cell viability and increase the release of proinflammatory cytokines, accompanied by the activation of the NF-κB pathway (So et al. 2007, 2008). NF-κB signaling could regulate the inflammatory processes in the central nervous system and play an important role in synaptic transmission and neuronal plasticity (O'Neill and Kaltschmidt 1997). Interestingly, our pathway enrichment analysis suggested similar pathways of synaptic signaling (such as “synaptic membrane”, “neuron to neuron synapse”, and “modulation of chemical synaptic transmission”) might participate in the development of CRIHL. Therefore, the current study provided potential clues between rs1050851 and the risk of CRIHL, and further study should be conducted to explore the underlying mechanism of NFKBIA in the development of chemoradiation-induced ototoxicity, such as through the neuroinflammatory processes and synaptic signaling pathway.
Moreover, we found a genetic similarity between hereditary deafness and CRIHL by permutation-based enrichment analysis. The SNPs within the regions of 149 hereditary deafness genes were more likely to be correlated with CRIHL than expected by chance and were significantly enriched among the GWAS top signals of CRIHL. rs2275994 (downstream of GRHL2, a known hereditary deafness gene) was the most significant SNP associated with CRIHL among the regions of 149 hereditary deafness genes. Patients carrying the risk allele (T) of rs2275994 may have lower GRHL2 expression, which exerted risk effects on radiation damage. GRHL2 gene encodes a transcription factor involved in epithelial differentiation, morphogenesis, and maintenance (Werth et al. 2010). Mutations in this gene, usually truncating variants, have been reported to cause autosomal dominant non-syndromic hearing loss (Peters et al. 2002; Vona et al. 2013; Trebusak Podkrajsek et al. 2021). Grhl2 protein deficiency in a zebrafish model directly reduced the expression of cldnb and epcam, junction proteins in otic epithelial cells, resulting in abnormal development of the inner ear and impaired hearing (Han et al. 2011). In addition, several lines of evidence support our assumption that hereditary deafness and CRIHL share genetic etiology. For example, variants in WFS1, a gene involved in both autosomal non-syndromic hearing loss and Wolfram syndrome, showed significant associations with cisplatin-induced ototoxicity in testicular cancer survivors (Wheeler et al. 2017). Mutations in hereditary deafness genes, such as TRIOBP, ILDR1, and EYA4, also exhibited significant correlations with age-related hearing impairment (Hoffmann et al. 2016; Wells et al. 2019).
Our study also has some limitations. First, the phenotype of hearing loss in this study was from the medical records, which are mainly based on self-reported symptoms from patients and subjective inquiry from physicians, rather than objective assessment by pure tone audiometry. A systematic review demonstrated that self-reported hearing loss based on questionnaires could accurately distinguish subjects with hearing difficulty from those without hearing difficulty (Chou et al. 2011). Previous studies utilizing self-reported hearing status have also successfully identified several genetic variants for radiation-related hearing loss in childhood cancer survivors (Trendowski et al. 2021), and age-related hearing impairment in UK Biobank populations (Wells et al. 2019). Therefore, we think hearing evaluation based on subjective sensation may be a feasible and cost-effective measurement for alternative, when audiometry test is absent. However, future work is needed to clearly elaborate or fully clarify the definition of hearing loss through pure tone audiometry, otoscopy, acoustic impedance testing, and so on, to generate more objective and reliable results. Second, the radiation dosage parameters of the inner ear or cochlea were not available in our study. It is reported that the radiation dose to the inner ear or cochlea can much better reflect the hearing loss induced by radiation therapy (Pan et al. 2005; Chan et al. 2009; Theunissen et al. 2015). However, we made some efforts to overcome this limitation. We collected the information on the prescribed total dose to the nasopharynx and included tumor stage as a covariate in GWAS analysis, which may reflect the dosage-volume parameter of the cochlea dose to some degree, since patients in later tumor stages are more likely to receive higher radiation doses to the gross tumor as well as adjacent organs including cochlea. Additionally, we did stratification analysis by tumor stage or radiation dose to the nasopharynx and found the genetic effects remained similar and robust within different subgroups (Fig. S6). Third, we think future study with larger sample size is warranted to replicate the current findings and identify more variants for CRIHL. In current study, we tested previously reported SNPs associated with cisplatin or radiation-induced hearing loss for replication. Of note, rs34533789 (within the intron of ATXN1), which has been reported to be associated with radiation-induced hearing loss for childhood cancer survivors (Trendowski et al. 2021), was replicated in our study (HR = 1.35, P value = 0.04).
In summary, our study has increased the knowledge of clinical and genetic risk factors of CRIHL, which may assist in the development of accurate prediction models to identify high-risk individuals, and contribute to precision prevention to improve quality of life for cancer survivors. However, a multicenter study with a larger sample size is needed to validate our findings and discover additional variants for CRIHL. Additional functional studies to uncover the molecular basis of CRIHL are also required, which will prompt advances in protective interventions and less ototoxic treatments.
Availability of data and material
The datasets used and/or analyzed during the current study are available at Research Data Deposit (RDD) public platform (http://www.researchdata.org.cn, accession number: RDDA2023873981).
References
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, Wilson CJ, Lehar J, Kryukov GV, Sonkin D et al (2012) The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature 483:603–607
Bhandare N, Jackson A, Eisbruch A, Pan CC, Flickinger JC, Antonelli P, Mendenhall WM (2010) Radiation therapy and hearing loss. Int J Radiat Oncol Biol Phys 76:S50-57
Chan SH, Ng WT, Kam KL, Lee MC, Choi CW, Yau TK, Lee AW, Chow SK (2009) Sensorineural hearing loss after treatment of nasopharyngeal carcinoma: a longitudinal analysis. Int J Radiat Oncol Biol Phys 73:1335–1342
Chen WC, Jackson A, Budnick AS, Pfister DG, Kraus DH, Hunt MA, Stambuk H, Levegrun S, Wolden SL (2006) Sensorineural hearing loss in combined modality treatment of nasopharyngeal carcinoma. Cancer 106:820–829
Chen YP, Chan ATC, Le QT, Blanchard P, Sun Y, Ma J (2019) Nasopharyngeal carcinoma. Lancet 394:64–80
Chou R, Dana T, Bougatsos C, Fleming C, Beil T (2011) Screening adults aged 50 years or older for hearing loss: a review of the evidence for the U.S. preventive services task force. Ann Intern Med 154:347–355
Delmaghani S, del Castillo FJ, Michel V, Leibovici M, Aghaie A, Ron U, Van Laer L, Ben-Tal N, Van Camp G, Weil D et al (2006) Mutations in the gene encoding pejvakin, a newly identified protein of the afferent auditory pathway, cause DFNB59 auditory neuropathy. Nat Genet 38:770–778
Diaz-Horta O, Abad C, Sennaroglu L, Foster J 2nd, DeSmidt A, Bademci G, Tokgoz-Yilmaz S, Duman D, Cengiz FB, Grati M et al (2016) ROR1 is essential for proper innervation of auditory hair cells and hearing in humans and mice. Proc Natl Acad Sci USA 113:5993–5998
El Charif O, Mapes B, Trendowski MR, Wheeler HE, Wing C, Dinh PC Jr, Frisina RD, Feldman DR, Hamilton RJ, Vaughn DJ et al (2019) Clinical and genome-wide analysis of cisplatin-induced tinnitus implicates novel ototoxic mechanisms. Clin Cancer Res 25:4104–4116
Ghisdal L, Baron C, Lebranchu Y, Viklicky O, Konarikova A, Naesens M, Kuypers D, Dinic M, Alamartine E, Touchard G et al (2017) Genome-wide association study of acute renal graft rejection. Am J Transplant 17:201–209
Han Y, Mu Y, Li X, Xu P, Tong J, Liu Z, Ma T, Zeng G, Yang S, Du J et al (2011) Grhl2 deficiency impairs otic development and hearing ability in a zebrafish model of the progressive dominant hearing loss DFNA28. Hum Mol Genet 20:3213–3226
He YQ, Wang TM, Ji M, Mai ZM, Tang M, Wang R, Zhou Y, Zheng Y, Xiao R, Yang D et al (2022) A polygenic risk score for nasopharyngeal carcinoma shows potential for risk stratification and personalized screening. Nat Commun 13:1966
Hoffmann TJ, Keats BJ, Yoshikawa N, Schaefer C, Risch N, Lustig LR (2016) A large genome-wide association study of age-related hearing impairment using electronic health records. PLoS Genet 12:e1006371
Jiang H, Sha SH, Schacht J (2005) NF-kappaB pathway protects cochlear hair cells from aminoglycoside-induced ototoxicity. J Neurosci Res 79:644–651
Kwong DL, Wei WI, Sham JS, Ho WK, Yuen PW, Chua DT, Au DK, Wu PM, Choy DT (1996) Sensorineural hearing loss in patients treated for nasopharyngeal carcinoma: a prospective study of the effect of radiation and cisplatin treatment. Int J Radiat Oncol Biol Phys 36:281–289
Landier W (2016) Ototoxicity and cancer therapy. Cancer 122:1647–1658
Lang H, Schulte BA, Zhou D, Smythe N, Spicer SS, Schmiedt RA (2006) Nuclear factor κB deficiency is associated with auditory nerve degeneration and increased noise-induced hearing loss. J Neurosci 26:3541–3550
Low WK, Toh ST, Wee J, Fook-Chong SM, Wang DY (2006) Sensorineural hearing loss after radiotherapy and chemoradiotherapy: a single, blinded, randomized study. J Clin Oncol 24:1904–1909
Manem VS, Lambie M, Smith I, Smirnov P, Kofia V, Freeman M, Koritzinsky M, Abazeed ME, Haibe-Kains B, Bratman SV (2019) Modeling cellular response in large-scale radiogenomic databases to advance precision radiotherapy. Cancer Res 79:6227–6237
Meijer AJM, Diepstraten FA, Langer T, Broer L, Domingo IK, Clemens E, Uitterlinden AG, de Vries ACH, van Grotel M, Vermeij WP et al (2021) TCERG1L allelic variation is associated with cisplatin-induced hearing loss in childhood cancer, a PanCareLIFE study. NPJ Precis Oncol 5:64
Mujica-Mota MA, Lehnert S, Devic S, Gasbarrino K, Daniel SJ (2014) Mechanisms of radiation-induced sensorineural hearing loss and radioprotection. Hear Res 312:60–68
O’Neill LA, Kaltschmidt C (1997) NF-κ B: a crucial transcription factor for glial and neuronal cell function. Trends Neurosci 20:252–258
Pan CC, Eisbruch A, Lee JS, Snorrason RM, Ten Haken RK, Kileny PR (2005) Prospective study of inner ear radiation dose and hearing loss in head-and-neck cancer patients. Int J Radiat Oncol Biol Phys 61:1393–1402
Peters LM, Anderson DW, Griffith AJ, Grundfast KM, San Agustin TB, Madeo AC, Friedman TB, Morell RJ (2002) Mutation of a transcription factor, TFCP2L3, causes progressive autosomal dominant hearing loss, DFNA28. Hum Mol Genet 11:2877–2885
Pfister DG, Spencer S, Adelstein D, Adkins D, Anzai Y, Brizel DM, Bruce JY, Busse PM, Caudell JJ, Cmelak AJ et al (2020) Head and neck cancers, version 2.2020, NCCN clinical practice guidelines in oncology. J Natl Compr Canc Netw 18:873–898
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D (2006) Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 38:904–909
Rizvi AA, Karaesmen E, Morgan M, Preus L, Wang J, Sovic M, Hahn T, Sucheston-Campbell LE (2019) gwasurvivr: an R package for genome-wide survival analysis. Bioinformatics 35:1968–1970
Romano A, Capozza MA, Mastrangelo S, Maurizi P, Triarico S, Rolesi R, Attina G, Fetoni AR, Ruggiero A (2020) Assessment and management of platinum-related ototoxicity in children treated for cancer. Cancers (basel) 12(5):1266
Ruel J, Emery S, Nouvian R, Bersot T, Amilhon B, Van Rybroek JM, Rebillard G, Lenoir M, Eybalin M, Delprat B et al (2008) Impairment of SLC17A8 encoding vesicular glutamate transporter-3, VGLUT3, underlies nonsyndromic deafness DFNA25 and inner hair cell dysfunction in null mice. Am J Hum Genet 83:278–292
Schoen CJ, Emery SB, Thorne MC, Ammana HR, Sliwerska E, Arnett J, Hortsch M, Hannan F, Burmeister M, Lesperance MM (2010) Increased activity of diaphanous homolog 3 (DIAPH3)/diaphanous causes hearing defects in humans with auditory neuropathy and in Drosophila. Proc Natl Acad Sci USA 107:13396–13401
Shi W, Hou X, Bao X, Hou W, Jiang X, Ma L, Jiang X, Dong L (2021) Mechanism and protection of radiotherapy induced sensorineural hearing loss for head and neck cancer. Biomed Res Int 2021:3548706
Shih VF, Tsui R, Caldwell A, Hoffmann A (2011) A single NFκB system for both canonical and non-canonical signaling. Cell Res 21:86–102
Simon M, Richard EM, Wang X, Shahzad M, Huang VH, Qaiser TA, Potluri P, Mahl SE, Davila A, Nazli S et al (2015) Mutations of human NARS2, encoding the mitochondrial asparaginyl-tRNA synthetase, cause nonsyndromic deafness and Leigh syndrome. PLoS Genet 11:e1005097
Smirnov P, Safikhani Z, El-Hachem N, Wang D, She A, Olsen C, Freeman M, Selby H, Gendoo DM, Grossmann P et al (2016) PharmacoGx: an R package for analysis of large pharmacogenomic datasets. Bioinformatics 32:1244–1246
So H, Kim H, Lee JH, Park C, Kim Y, Kim E, Kim JK, Yun KJ, Lee KM, Lee HY et al (2007) Cisplatin cytotoxicity of auditory cells requires secretions of proinflammatory cytokines via activation of ERK and NF-κB. J Assoc Res Otolaryngol 8:338–355
So H, Kim H, Kim Y, Kim E, Pae HO, Chung HT, Kim HJ, Kwon KB, Lee KM, Lee HY et al (2008) Evidence that cisplatin-induced auditory damage is attenuated by downregulation of pro-inflammatory cytokines via Nrf2/HO-1. J Assoc Res Otolaryngol 9:290–306
Starr A, Picton TW, Sininger Y, Hood LJ, Berlin CI (1996) Auditory neuropathy. Brain 119(Pt 3):741–753
Tang LQ, Chen DP, Guo L, Mo HY, Huang Y, Guo SS, Qi B, Tang QN, Wang P, Li XY et al (2018) Concurrent chemoradiotherapy with nedaplatin versus cisplatin in stage II-IVB nasopharyngeal carcinoma: an open-label, non-inferiority, randomised phase 3 trial. Lancet Oncol 19:461–473
Theunissen EA, Zuur CL, Jozwiak K, Lopez-Yurda M, Hauptmann M, Rasch CR, van der Baan S, de Boer JP, Dreschler WA, Balm AJ (2015) Prediction of hearing loss due to cisplatin chemoradiotherapy. JAMA Otolaryngol Head Neck Surg 141:810–815
Trebusak Podkrajsek K, Tesovnik T, Bozanic Urbancic N, Battelino S (2021) Novel GRHL2 gene variant associated with hearing loss: a case report and review of the literature. Genes (basel) 12(4):484
Trendowski MR, El Charif O, Dinh PC Jr, Travis LB, Dolan ME (2019) Genetic and modifiable risk factors contributing to cisplatin-induced toxicities. Clin Cancer Res 25:1147–1155
Trendowski MR, Baedke JL, Sapkota Y, Travis LB, Zhang X, El Charif O, Wheeler HE, Leisenring WM, Robison LL, Hudson MM et al (2021) Clinical and genetic risk factors for radiation-associated ototoxicity: a report from the Childhood Cancer Survivor Study and the St. Jude Lifetime Cohort. Cancer 127:4091–4102
Vallabhapurapu S, Karin M (2009) Regulation and function of NF-κB transcription factors in the immune system. Annu Rev Immunol 27:693–733
Vona B, Nanda I, Neuner C, Muller T, Haaf T (2013) Confirmation of GRHL2 as the gene for the DFNA28 locus. Am J Med Genet A 161A:2060–2065
Wang Y, Fong PY, Leung FC, Mak W, Sham PC (2007) Increased gene coverage and Alu frequency in large linkage disequilibrium blocks of the human genome. Genet Mol Res 6:1131–1141
Wang J, Chen YY, Tai A, Chen XL, Huang SM, Yang C, Bao Y, Li NW, Deng XW, Zhao C et al (2015) Sensorineural hearing loss after combined intensity modulated radiation therapy and cisplatin-based chemotherapy for nasopharyngeal carcinoma. Transl Oncol 8:456–462
Wang TM, Shen GP, Chen MY, Zhang JB, Sun Y, He J, Xue WQ, Li XZ, Huang SY, Zheng XH et al (2019) Genome-wide association study of susceptibility loci for radiation-induced brain injury. J Natl Cancer Inst 111:620–628
Wells HRR, Freidin MB, Zainul Abidin FN, Payton A, Dawes P, Munro KJ, Morton CC, Moore DR, Dawson SJ, Williams FMK (2019) GWAS identifies 44 independent associated genomic loci for self-reported adult hearing difficulty in UK biobank. Am J Hum Genet 105:788–802
Wen DW, Lin L, Mao YP, Chen CY, Chen FP, Wu CF, Huang XD, Li ZX, Xu SS, Kou J et al (2021) Normal tissue complication probability (NTCP) models for predicting temporal lobe injury after intensity-modulated radiotherapy in nasopharyngeal carcinoma: a large registry-based retrospective study from China. Radiother Oncol 157:99–105
Werth M, Walentin K, Aue A, Schonheit J, Wuebken A, Pode-Shakked N, Vilianovitch L, Erdmann B, Dekel B, Bader M et al (2010) The transcription factor grainyhead-like 2 regulates the molecular composition of the epithelial apical junctional complex. Development 137:3835–3845
Wheeler HE, Gamazon ER, Frisina RD, Perez-Cervantes C, El Charif O, Mapes B, Fossa SD, Feldman DR, Hamilton RJ, Vaughn DJ et al (2017) Variants in WFS1 and other mendelian deafness genes are associated with cisplatin-associated ototoxicity. Clin Cancer Res 23:3325–3333
Wu F, Wang R, Lu H, Wei B, Feng G, Li G, Liu M, Yan H, Zhu J, Zhang Y et al (2014) Concurrent chemoradiotherapy in locoregionally advanced nasopharyngeal carcinoma: treatment outcomes of a prospective, multicentric clinical study. Radiother Oncol 112:106–111
Wu K, Wang H, Guan J, Lan L, Zhao C, Zhang M, Wang D, Wang Q (2020) A novel variant in diaphanous homolog 1 (DIAPH1) as the cause of auditory neuropathy in a Chinese family. Int J Pediatr Otorhinolaryngol 133:109947
Xu H, Robinson GW, Huang J, Lim JY, Zhang H, Bass JK, Broniscer A, Chintagumpala M, Bartels U, Gururangan S et al (2015) Common variants in ACYP2 influence susceptibility to cisplatin-induced hearing loss. Nat Genet 47:263–266
Yang W, Soares J, Greninger P, Edelman EJ, Lightfoot H, Forbes S, Bindal N, Beare D, Smith JA, Thompson IR et al (2013) Genomics of Drug Sensitivity in Cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res 41:D955-961
Yao JJ, Yu XL, Zhang F, Zhang WJ, Zhou GQ, Tang LL, Mao YP, Chen L, Ma J, Sun Y (2017) Radiotherapy with neoadjuvant chemotherapy versus concurrent chemoradiotherapy for ascending-type nasopharyngeal carcinoma: a retrospective comparison of toxicity and prognosis. Chin J Cancer 36:26
Yard BD, Adams DJ, Chie EK, Tamayo P, Battaglia JS, Gopal P, Rogacki K, Pearson BE, Phillips J, Raymond DP et al (2016) A genetic basis for the variation in the vulnerability of cancer to DNA damage. Nat Commun 7:11428
Yu G, Wang LG, Han Y, He QY (2012) clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16:284–287
Zhang G, Zheng H, Pyykko I, Zou J (2019a) The TLR-4/NF-κB signaling pathway activation in cochlear inflammation of rats with noise-induced hearing loss. Hear Res 379:59–68
Zhang Y, Chen L, Hu GQ, Zhang N, Zhu XD, Yang KY, Jin F, Shi M, Chen YP, Hu WH et al (2019b) Gemcitabine and cisplatin induction chemotherapy in nasopharyngeal carcinoma. N Engl J Med 381:1124–1135
Funding
This study was funded by the National Key Research and Development Program of China (2021YFC2500400), the Basic and Applied Basic Research Foundation of Guangdong Province, China (2021B1515420007), the Special Support Program for High-level Professionals on Scientific and Technological Innovation of Guangdong Province, China (2014TX01R201), National Natural Science Foundation of China (81903395, 82003520), the Key Area Research and Development Program of Guangdong Province, China (2019B110233004).
Author information
Authors and Affiliations
Contributions
YQH performed the quality control of data, data analysis and interpretation, and manuscript editing. LTL performed data collection, data analysis and interpretation, and manuscript draft. TMW provided support for algorithms and manuscript editing. WQX, DWY, DHL, HD, RWX, CMD, WLZ, YL, YXW, QLW, TZ performed data collection and data analysis. XZL, XHZ, PFZ, SDZ, YZH performed sample collection and data acquisition. YS participated in manuscript review. WHJ conceived and designed the study, revised and reviewed the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors have declared no conflict of interest.
Ethical approval
The study protocol was approved by the Institutional Review Board of Sun Yat-sen University Cancer Center (B2019-026-01).
Consent to participate
All subjects have provided written informed consents.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
He, YQ., Luo, LT., Wang, TM. et al. Clinical and genome-wide association analysis of chemoradiation-induced hearing loss in nasopharyngeal carcinoma. Hum Genet 142, 759–772 (2023). https://doi.org/10.1007/s00439-023-02554-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00439-023-02554-0