Abstract
Background/Aim: The aim of the present study was to investigate the clinical significance of 7 single nucleotide polymorphisms (SNPs) of vascular endothelial growth factor A (VEGFA), endothelin receptor type A (EDNRA), nibrin (NBS1) and Fas cell surface death receptor (FAS) genes in patients with nasopharyngeal carcinoma (NPC). Patients and Methods: Genomic DNA was extracted from the peripheral blood specimens of 60 patients. Genotyping of 4 VEGFA polymorphisms, namely VEGFA -1154 G/A (rs1570360), -2578 A/C (rs699947), -1498 C/T (rs833061) and +936 C/T (rs3025039), as well as EDNRA SNP p.His323 (rs5333), NBS1 p.E185Q (rs1805794) and FAS -671 A/G (rs1800682) was performed by using Sanger sequencing. Results: The VEGFA +936 CC genotype was more frequent in tumors with bilateral infiltration of pterygoid plates compared to those with ipsilateral (76.9% vs. 69.6%, p=0.008) and no infiltration of pterygoid plates (76.9% vs. 68.8%, p=0.023). VEGFA -2578, VEGFA -1154 and VEGFA +936 were significantly associated with infiltration to the pterygoid processes (p=0.011, p=0.041 and p=0.032, respectively). EDNRA H323H TT genotype was marginally associated with infiltration to the ipsilateral medial pterygoid muscles (p=0.045). A trend towards longer overall survival was observed for VEGFA -2578 CC as compared to AC (HR=0.24, p=0.060). Conclusion: The studied VEGFA SNPs seem to be associated with the local extension of the NPC and maybe with the clinical outcome of this patient group.
Nasopharyngeal carcinoma (NPC) is a rare human malignancy with an annual incidence of 0.1-1.1/100,000 people with Caucasian origin in western countries. Particularly in Southeast Europe, the annual incidence reaches 1 per 100,000 inhabitants (1, 2). In contrast, the highest incidence rates are noted in Southern Asia, Alaska and the Mediterranean Basin (1). The unbalanced geographical incidence of this tumor may represent the role of the environment and the genetic background in the development of NPC (3, 4).
Imaging diagnosis of NPC is based on magnetic resonance imaging (MRI) along with computed tomography (CT), which accurately provide definition of the tumor extension and determine the treatment choice (5-9). The cornerstone of management of patients with early-stage NPC is radiation therapy (RT), while the combination of RT with chemotherapy is necessary for the management of locally advanced disease (10). Despite the combined imaging applications with RT and systemic therapies treatment still fails, especially when the tumor is in the advanced stage (11). Hence, prognostic factors, such as gene single-nucleotide polymorphisms (SNPs), needed to be further developed (12).
Many studies have shown that NPC pathogenesis is strongly associated with Epstein-Barr virus (EBV) infection, independently of ethnic origin, and for that reason NPC has been mainly characterized as a virus-related human malignancy (13, 14). In addition, it seems that EBV-related NPC is associated with a specific histological subtype. The undifferentiated type of NPC, strictly associated with EBV infection, is predominant in endemic areas, while the differentiated and non-keratinizing squamous cell variants are common in the USA and North European countries (11).
Except for EBV infection, multiple other factors participate in the pathogenesis of NPC, such as genetic, epigenetic and environmental factors as well as dietary habits (15). On the other hand, exposure to both radiation and chemicals can lead to lesions of genomic DNA. Double strand breaks (DSB) are a relatively dangerous form of DNA lesions, which, in combination with the failure of DNA repair mechanisms, can lead to genomic instability and probably cancer (6-8). Nibrin (NBS1) plays an important role in the cellular response to DNA damage and maintenance of genomic integrity. NBS1, as a part of a greater complex MRN (MRE11/RAD50/NBS1), modulates the DNA damage signal and promotes the repair (16). Any genetic variation of the NBS1 gene may have deleterious effects on cell capability to respond to DNA damage, predisposing to cancer and also to the genetic disease named Nijmegen Breakage Syndrome (NBS) (17).
Several studies have investigated a genetic variant of NBS1, E185Q polymorphism, in various types of cancer. A meta-analysis, a decade ago, showed that there might be an association between NBS1 8360 G/C (E185Q) polymorphism and cancer risk, especially in Caucasians (18). Zeng et al., who carried out a hospital-based case–control study to investigate the possible relationship of the E185Q polymorphism with the risk of NPC in Chinese population, found a correlation with advanced NPC stage and with progression of the disease to invasive and metastatic (18, 19). Although the NBS1 E185Q polymorphism is a commonly studied polymorphism, this SNP has not been studied regarding its role in local extension of NPC.
In addition, angiogenesis also plays a significant role in NPC development and progression (20, 21). VEGFA is a central molecule in the angiogenesis process and for that reason numerous studies have investigated the clinical significance of VEGFA variants in cancer (22-24). In the current study, VEGFA -1154 G/A (rs1570360), -2578 A/C (rs699947), -1498 C/T (rs833061) and +936 C/T (rs3025039) have been evaluated.
The axis of endothelin-1 (END1)/endothelin A receptor (EDNRA) also plays a significant role in angiogenesis. Focusing on the NPC, many publications have shown that EDNRA is overexpressed in 73.9% of NPC cases and it is also an independent prognostic factor for distant metastases (25). Furthermore, it has been shown that EDNRA SNP rs5333 has a prognostic significance for patients with locoregionally advanced NPC (26), as well as it has been associated with lower risk (grade 2-3) for radiation-induced xerostomia in NPC patients (27).
Another pivotal cellular function, a hallmark of cancer, is the apoptosis, which also exerts a central role in NPC (1). FAS is a cell surface receptor that can interact with its natural FAS ligand (FASL) to trigger apoptosis. Mutations in the FAS gene and its ligand can contribute to cell transformation by blocking the apoptotic signal transduction pathway, leading to loss of balance between proliferation and cell death (28). In addition, variations of FAS/FASL have been associated with impaired apoptosis in cancer. Particularly, in a meta-analysis of 52 studies, Xu et al. evaluated certain polymorphisms and found associations between different types of cancer. FAS -670 A/G polymorphism was associated with decreased risk of prostate cancer and melanoma, while the FAS -1377 AA genotype had an increased risk for breast, gastric and esophageal cancer. In addition, the FASL -844 C/T polymorphism was correlated with increased risk for gastric, esophageal and ovarian cancer (29). Also, a specific SNP of FAS, -670 A/G (rs1800682), has attracted the interest as it has been associated with esophageal cancer, hepatocellular carcinoma, and prostate cancer susceptibility (30), as well as with papillary thyroid cancer (31). On the contrary, rs1800682 wasn't associated with cervical cancer risk (32).
In this study, we investigated the possible association between the aforementioned SNPs with clinicopathological parameters, disease extension to parapharyngeal spaces, infiltration to adjacent structures, metastasis to regional lymph nodes as well as with the clinical outcome of NPC patients.
Patients and Methods
Patients. A cohort of 60 Greek patients with mainly locally advanced NPC, were retrieved from the Hellenic Cooperative Oncology Group (HeCOG) electronic database. Their formalin-fixed paraffin-embedded (FFPE) blocks and gDNA were also retrieved from the HeCOG's Tumor Repository and were analyzed in the Laboratory of Molecular Oncology (MOL, by the Hellenic Foundation for Cancer Research/Aristotle University of Thessaloniki, Thessaloniki, Greece). The study protocol was approved by the Bioethics Committee of the Aristotle University of Thessaloniki, School of Health Sciences, Faculty of Medicine (#2/23March2016) and was designed according to the ethical guidelines of the Helsinki Declaration (2013) (33). Written informed consents were obtained from all participants.
The clinical, CT and MRI data of these patients, who were diagnosed with NPC, were retrospectively reviewed. We evaluated the involvement or infiltration of structures adjacent to the nasopharynx (especially the base of the skull) and of the parapharyngeal space. TNM staging was applied according to the revised 8th edition of the AJCC staging system (Table I). Tumor without extension to the parapharyngeal space (line from the free edge of the medial pterygoid plate to the lateral border of the carotid artery) was categorized as degree A (Figure 1). Tumor with extension to parapharyngeal space and infiltration of the medial/lateral pterygoid, and prevertebral muscles (line from the base of the medial pterygoid plate to the styloid process) was categorized as degree B (Figure 2). Tumor invasion of the bony structures of skull base and pterygoid structures (line from the edge of the lateral pterygoid process to the posterior border of the ascending branch of the mandible) was categorized as degree C (Figure 3). Tumor with intracranial extension and involvement of cranial nerves and hypopharynx was categorized as degree C (Figure 4). Patients were treated with induction chemotherapy followed by concurrent chemoradiotherapy or chemoradiotherapy alone.
SNPs selection. Literature (PubMed, Google Scholar) was reviewed by using suitable keywords (SNPs, nasopharyngeal cancer, risk, MRI, angiogenesis, apoptosis, biomarker, cancer). Particular SNPs were found in VEGFA, EDNRA, NBS1 and FAS genes, which had been previously assessed regarding their role in cancer. Their selection was based on the functionality as well as minor allele frequency. Based on these, VEGFA -1154 G/A (rs1570360), -2578 A/C (rs699947), -1498 C/T (rs833061), located within the promoter and +936 C/T (rs3025039) in the 3’-UTR region, as well as EDNRA SNP p.His323 (rs5333), a synonymous SNP in exon 6, NBS1 p.E185Q (rs1805794), a non-synonymous SNP within exon 5, and FAS -670 A/G (rs1800682), a SNP in the 5’-end regulatory, were finally selected. In all instances, the minus symbol (-) describes nucleotides before translation start site and plus (+) nucleotides following translation end site.
SNPs genotyping. DNA was extracted from peripheral blood samples corresponding to the 60 patients under study, with a standard desalting method (34). Genotyping for the selected SNPs of VEGFA, EDNRA, FAS and NBS1 was performed with Sanger sequencing, as previously described (35).
The primers used for detection of VEGFA SNPs were (sense and antisense): -2578 A/G/C/T, 5’-AGGATGGGGCTGACTAGGTAA-3’ and 5’-CCCCCTTTTCCTCCAACTCTC-3’, 324 bp amplicon; -1498 C/T 5’-GCCCATTCCCTCTTTAGCCA-3’ and 5’-AGTGAGGTTACGTGCGGACAG-3’, 325 bp; -1154 G/A, 5’-CTGCTCCCTCCTCGCCAATGC -3’ and 5’-CCAAGCCTCCGCGATCCT-3’, 308 bp; +936 C/T, 5’-TCACCAGGAAAGACTGATACA-3’ and 5’-GGTGGGTGTGTCTACAGGA-3’, 295 bp. Furthermore, the primers used for detection of EDNRA p.His323His were 5’-ATTCTTTCTCTGGTGTCTGC-3’ and 5’-GAAAATCTGAGAAACTCCAAT-3’, 281 bp; of FAS -670A/G were 5’-TGCGATTTGGCTTAAGTTGT-3’ and 5’-GGCTTCTGCTGTAGTTCAACC-3’, 312 bp; of NBS1 p.E185Q 5’-AGAGAGATGAAAGGGAAA-3’ and 5’- ATTACATCCTGAAACAAGCAT-3’, 333 bp amplicon, respectively.
Genotypic analysis of all SNPs involved polymerase chain reaction (PCR) amplification of the target gene region with the M13-coupled primers mentioned above, followed by Sanger sequencing. All PCR reactions were performed on a GeneAmp PCR system, while sense and antisense sequencing was carried out, using M13 forward and reverse primers, with the Big Dye Teminator kit v.1.1 (AppliedBiosystems, Foster City, CA, USA), according to the manufacturer's instructions. Samples were then purified utilizing the ethanol/EDTA precipitation protocol and loaded on the ABI 3130xl Genetic Analyzer. Sequences were visualized upon capillary electrophoresis in the ABI3130 genetic analyzer, initially called with the Sequencing Analysis software and further analyzed with the SeqScape software v2.5 (Applied Biosystems/Biosolutions, Athens).
Evaluation of EBV status. For the evaluation of EBV status, FFPE tumor tissue blocks were retrieved from the HeCOG Tumor Repository. EBV status was determined based on the assessment of EBV-encoded small RNA (EBER) using either the chromogenic in situ hybridization (CISH) or the polymerase chain reaction (PCR).
In the case of EBER detection by CISH, FFPE blocks underwent histological review by an experienced pathologist, who marked tumor dense areas for Tissue Microarray (TMA) construction. Upon histological evaluation, available tumors were transferred into low-density TMAs, whereby three cores (0.6 mm in diameter) were obtained from appropriate regions of each tumor with the use of a manual arrayer (Model I, Beecher Instruments, San Prairie, WI) and included in the TMA block. Subsequently, EBER detection was performed on 4 μm sections from the TMA block with the use of a Bond Max autostainer (Leica Microsystems, Wetzlar, Germany), according to the provider's specifications, as previously described (34-37). Tumor cells containing EBER transcripts were considered positive when an intense, brown, predominantly nuclear staining was present in >1% of tumor cells (38). Only positive cases for the poly-A probe, which is indicative of the preservation of intact mRNA in cells was considered eligible for assessment of EBER status. For PCR-based EBV detection, DNA was extracted from 10 μm FFPE sections, macrodissected manually to obtain tumor dense areas as marked on respective H&E sections by the pathologist, using the QIAamp DNA Mini Kit (Qiagen GmbH, Hilden, Germany), according to standard procedures. For PCR-based EBV detection, specific primers were used targeting a segment of the BamH1W region, as described previously.
Statistical analysis. Group comparisons of categorical data were performed with the Chi-square or Fisher's exact (where appropriate) test. Overall survival (OS) was measured from the date of diagnosis to the date of patient's death or last contact, whichever occurred first. Survival curves were estimated using the Kaplan Meier method and compared across groups with the log-rank test. The associations between factors of interest and death rates were assessed with hazard ratios estimated with Cox proportional hazard regression models. Statistical significance was set at a two-sided p<0.05. Statistical analyses were performed using the Statistical Analyses Software (SAS; version 9.3, SAS Institute Inc., Cary, NC).
Results
Study population. Clinicopathological characteristics of the studied patients as well as tumor radiological parameters are shown in Table II. A total of 60 patients with NPC, aged between 14 and 78 years, were included in this study with 78.3% of the study population being males and 21.7% females. Histological subtype was confirmed from biopsies taken from the primary site or from infiltrated cervical lymph nodes in a few cases. Thirty-five patients (58.3%) had an undifferentiated carcinoma, 7 patients (11.7%) had squamous cell carcinoma (SCC) and 18 (30%) had a non-keratinising carcinoma. EBV status was available for 47 patients, with 43 of them were positive (91.5%) and 4 were negative (8.5%). The majority of patients (81.7%) had stage II and III disease, while 10 patients had stage IV disease. Lateral extension A was observed in 20 patients (33.3%), lateral extension B in 43.3% and C in 23.3% of the study cohort. It is of note that only 7 patients (11.7%) had intracranial involvement, while 50% of the patients had one wall and the other 50% had two walls extension of the tumor. During the follow-up period [median 133.34 months (95%CI=119.8-147.18)], one patient was lost from follow-up and 18 deaths (30%) were reported. Thirteen of them (72.2%) were due to NPC, 3 to other tumors, while one patient died of pneumonia and another due to old age. The median OS had not been reached yet at the time of the analysis. The 3- and 5-year OS rates were 86.2% and 80.8%, respectively.
Genotype frequencies. The frequency distributions of the examined SNPs are presented in Table III. The genotyping results showed that most patients were heterozygous for VEGFA -2578 AC (50%), VEGFA -1498 CT (45%), VEGFA -1154 GG (41.1%) and FAS -670 AG (66.1%).
VEGFA +936 was associated with infiltration to the pterygoid plates. A significant association was observed between infiltration to the pterygoid plates and VEGFA +936 (p=0.008). Since there was only one patient with infiltration of contralateral pterygoid plates, the 2×2 associations between the rest of the categories were examined to detect where the true association lies. The VEGFA +936 CC genotype was found to be more frequent in tumors with bilateral infiltration of pterygoid plates as compared to those with ipsilateral (76.9% vs. 69.6%, p=0.008) and no infiltration of pterygoid plates (76.9% vs. 68.8%, p=0.023).
VEGFA -2578, VEGFA -1154 and VEGFA +936 were significantly associated with infiltration of the pterygoid processes. Additionally, VEGFA -2578, VEGFA -1154 and VEGFA +936 were significantly associated with infiltration to the pterygoid processes (p=0.011, p=0.041 and p=0.032, respectively). Examining the 2×2 associations among negative, ipsilateral and bilateral infiltration to pterygoid process, VEGFA -2578 AC was more frequent in tumors with ipsilateral infiltration to the pterygoid process as compared to those with negative infiltration (71.4% vs. 35.3%, p=0.004) as well as compared to those with bilateral infiltration of pterygoid process (71.4% vs. 40%, p=0.035). VEGFA -1154 GG was more frequent in tumors with infiltration of the bilateral pterygoid process as opposed to those with no infiltration (73.3% vs. 21.1%, p=0.008) and VEGFA +936 CC was more frequently observed in tumors with bilateral than in tumors with ipsilateral infiltration of the pterygoid process (84.6% vs. 66.7%, p=0.016). However, the power of these results is limited by the small number of patients in the groups under comparison and interpretation should be conducted with caution until further validation in larger cohorts.
SNPs and clinical outcome. In terms of OS, only a trend towards longer survival was observed for VEGFA -2578 CC as compared to AC (HR=0.24, 95%CI=0.05-1.06, p=0.060), while the rest of the examined polymorphisms were not found to be prognostic for OS (Table IV). Upon adjustment for age, the only parameter that was found to marginally significantly affect OS in the univariate analysis, the hazard ratio for VEGFA -2578 CC was 0.28 (95%CI=0.06-1.27), however statistical significance was not reached (p=0.098). As previously, the power of the aforementioned results is reduced by the small number of patients experiencing the event of interest in the comparable groups and need to be further validated in larger cohorts of NPC patients.
Associations with clinicopathological parameters. The associations of SNPs with selected clinical features are presented in Table V. Patients with undifferentiated tumors were more frequently homozygous GG for VEGFA -1154 (p=0.030). In addition, smokers were less likely to have either VEGFA -1498 CT genotype (30.4% vs. 60%, p=0.043) or -2578 AC (44.1% vs. 61.1%, p=0.038) compared to non-smokers. No further significant associations were observed between the polymorphisms and the clinical characteristics.
No associations were observed with tumor extension. Similarly, no significant associations were observed between any of the examined polymorphisms and tumor location, lateral extension or infiltration to the prevertebral muscles (Table VI), while EDNRA H323H TT genotype was marginally significantly associated with infiltration to the ipsilateral medial pterygoid muscles (p=0.045).
Discussion
The study of SNPs and their association with certain clinicopathological features, tumor development, progression and overall survival has become a hot area of cancer research (39, 40). Over the years, various technologies have been developed for efficient SNP genotyping, including array-based hybridization, discrete PCR-based SNP assays as well as targeted Sanger sequencing, while recently the advent of next-generation sequencing technology (NGS) has allowed high-throughput but also flexible SNP genotypic analysis, proving advantage in terms of cost-effectiveness, number of SNPs and samples processed at once (41). This technological explosion has also helped the study of SNPs in the context of pharmacogenomics, contributing to the individualized management of cancer patients (42).
In the current study, seven SNPS of the VEGFA, EDNRA, NBS1 and FAS genes were evaluated in patients with NPC, assessing their correlations with clinicopathological parameters, radiologic traits and survival outcome of NPC patients. Three of these (VEGFA +936, VEGFA -2578, VEGFA -1154) seem to have a role in NPC. The fact that these SNPs are variants of the VEGFA gene, which is a significant component of angiogenesis, may reflects the important role of angiogenesis in NPC pathogenesis and development (20, 21). Our observation is supported from numerous studies, which have characterized these SNPs as functional and their status has been associated with risk, treatment response and prognosis of many malignant and non-malignant diseases (43-45). It seems that the studied SNPs exert a significant role in the pathogenesis of cancer affecting important functions, but their role needs to be clarified further in order to become clinically useful.
An interesting finding of the present study is the association of VEGFA +936 CC homozygosity with the infiltration to the pterygoid plates and processes. Although this SNP has been studied regarding its role in the NPC, this is the first report which associates VEGFA +936 with infiltration of NPC primary tumor adjacent tissues (46). Supportive to our observation is the fact the T allele of this SNP has been associated with plasma concentration (47) as well as with the expression in tumor tissues of VEGFA (48). Additionally, T allele of VEGFA +936 has been related to reduced risk for breast cancer in a cohort of Austrian individuals (49). Similarly to our findings, Liu et al. has documented that VEGFA +936 has no association with overall survival of patients with locoregionally advanced NPC treated with cisplatin- and fluorouracil-based chemoradiotherapy (50).
Our study further revealed an association of VEGFA -2578 AC genotype with increased frequency of ipsilateral infiltration of the pterygoid process. In addition, VEGFA -2578 CC genotype showed marginally favorable prognostic significance in univariate analysis, a relation which did not remain upon age adjustment. However, Wang et al. have demonstrated that NPC patients harboring the VEGFA −2578 CC genotype exhibited increased aggressiveness, larger size, poorer differentiation and advanced stage disease compared to patients harboring the VEGFA -2578 A allele (51). In addition, Nasr et al. have found a significant association between the VEGFA -2578 C allele carriers and aggressive forms of NPC, which were characterized by large tumors and advanced tumor stages (52).
Intriguing was also our finding that VEGFA -1154 GG was associated with infiltration to the bilateral pterygoid processes and to undifferentiated tumors. Interestingly, based on a meta-analysis from Hong et al., no association between VEGFA -1154G/A polymorphism and the risk of cancer has been described with the exception of non-Caucasians, who seem to have higher cancer risk (53). In addition, Masago et al. have shown that VEGFA -1154 is an independent prognostic factor for NSCLC patients (54). Similarly, Verboom et al. have reported that rs1570360 is associated with PFS in patients with advanced GIST receiving imatinib (55).
Various SNPs have also been identified on genes of EDN1 and EDNRA and have been associated with prognosis of some diseases. Colombo et al., reported that EDNRA H323H polymorphism was associated with a substantially increased risk of heart failure while Herrmann et al. have found that the EDNRA/H323H polymorphism predicts survival in patients with dilated cardiomyopathy (26, 56-59). In our group, the EDNRA H323H TT genotype was marginally significantly associated with infiltration to the ipsilateral medial pterygoid muscles. No statistically significant association was observed between the ENDRA polymorphism and tumor location or extension of the NPC.
Furthermore, it should be kept in mind that there are some limitations in the current study. First of all, a major limitation of our study is the small number of the enrolled patients, which didn't permit to show statistically significant associations. In addition, it would be more informative if a two-phase design had been followed and so, our results need further validation in a larger cohort.
In conclusion, three well described SNPs of the VEGFA gene (rs3025039, rs1570360, rs699947) seem to have an association with infiltration of adjacent tissues in patients with NPC, while rs699947 may also have prognostic value for the same patients. Further validation is needed in a larger cohort of NPC patients.
Footnotes
↵* These Authors contributed equally to this study.
Authors' Contributions
Elizabeth Psoma: resources, investigation, writing-original draft preparation. Georgia-Angeliki Koliou: formal analysis, writing-original draft preparation. Foteinos-Ioannis Dimitrakopoulos: writing-original draft preparation. Kyriaki Papadopoulou: investigation, writing-review and editing. Dimitra Rontogianni: investigation. Mattheos Bobos: investigation. Anastasios Visvikis: resources. Paris A. Kosmidis: resources. George Fountzilas: resources. John Constantinidis: resources. Anna Kalogera-Fountzila: conceptualization, investigation, supervision.
This article is freely accessible online.
Conflicts of Interest
All Authors declare no financial competing interests regarding this study.
- Received December 13, 2019.
- Revision received December 23, 2019.
- Accepted January 6, 2020.
- Copyright© 2020, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved