Abstract
Radiotherapy and cisplatin lead to cell killing in head and neck squamous cell carcinoma patients, but adverse events and response to treatment are not the same in patients with similar clinicopathological aspects. The aim of this prospective study was to evaluate the roles of TP53 c.215G > C, FAS c.-671A > G, FAS c.-1378G > A, FASL c.-844 C > T, CASP3 c.-1191A > G, and CASP3 c.-182-247G > T single nucleotide variants in toxicity, response rate, and survival of cisplatin chemoradiation-treated head and neck squamous cell carcinoma patients. Genomic DNA was analyzed by polymerase chain reaction for genotyping. Differences between groups of patients were analyzed by chi-square test or Fisher’s exact test, multiple logistic regression analysis, and Cox hazards model. One hundred nine patients with head and neck squamous cell carcinoma were enrolled in study. All patients were smokers and/or alcoholics. Patients with FAS c.-671GG genotype, FAS c.-671AG or GG genotype, and FASL c.-844CC genotype had 5.52 (95% confidence interval (CI): 1.42–21.43), 4.03 (95% CI: 1.51–10.79), and 5.77 (95% CI: 1.23–27.04) more chances of presenting chemoradiation-related anemia of grades 2–4, lymphopenia of grade 3 or 4, and ototoxicity of all grades, respectively, than those with the remaining genotypes. FAS c.-671GG genotype was also seen as an independent predictor of shorter event-free survival (hazard ratio (HR): 2.05; P = 0.007) and overall survival (HR: 1.83; P = 0.02) in our head and neck squamous cell carcinoma patients. These findings present, for the first time, preliminary evidence that inherited abnormalities in apoptosis pathway, related to FAS c.-671A > G and FASL c.-844 C > T single nucleotide variants, can alter toxicity and survival of tobacco- and alcohol-related head and neck squamous cell carcinoma patients homogeneously treated with cisplatin chemoradiation.
Introduction
Head and neck (HN) squamous cell carcinoma (SCC) accounts for only 5% of all new cancer cases reported worldwide, 1 but more than half of patients present advanced locoregional disease at diagnosis 2 and have 5 years overall survival (OS) rates lower than 50%. 3 HNSCC has tobacco smoking and alcohol consumption or human papillomavirus type 16 (HPV16) infection as risk factors. 4
Standard treatment for patients with locally advanced HNSCC, whose surgical resection is not amenable, is radiotherapy (RT) with cisplatin (CDDP). 5 Ionizing radiation (IR) leads to DNA single- and double-strand breaks promoting cancer cell death. 6 The cytotoxic effects of CDDP are mediated by binds with DNA forming interstrand and intrastrand cross-links. 7 Response to IR 8 and CDDP 9 damages involves P53 activation that stimulates the transcription of FAS receptor, which binds to the FAS ligand (FASL), and this interaction results in formation of the death-inducing signaling complex (DISC) and triggers apoptosis by activation of effector caspase 3 (CASP3).
The variability in adverse events and tumor response induced by RT 10 and CDDP 11 among patients receiving the same treatment can be attributed to genetic changes causing modifications in cellular phenotype. Functional studies revealed that variant alleles of TP53 c.215G > C, 12 FAS c.-671A > G, 13 FAS c.-1378G > A, 13 and FASL c.-844 C > T 14 single nucleotide variants (SNVs) reduce apoptosis capacity when compared with wild-type alleles. Reduced apoptosis was also associated with the variant alleles of CASP3 c.-1191A > G 15 and CASP3 c.-182-247G > T 16 SNVs, respectively, by in silico analysis.
TP53 c.215G > C was previously associated with myelotoxicity in ovarian cancer patients treated with CDDP-cyclophosphamide, 17 nephrotoxicity in lung cancer patients treated with CDDP or carboplatin, 18 gastrointestinal toxicity in elderly lung cancer patients treated with platinum plus navelbine, paclitaxel, docetaxel, or gemcitabine, 19 response to therapy, disease-free survival (DFS) and OS in lung cancer patients treated with platinum plus vinorelbine, paclitaxel, docetaxel, or gemcitabine, 20 and DFS in stage I or II radiation-treated HNSCC patients. 21 FAS c.-671A > G and FASL c.-844 C > T modulated DFS in patients with SCC of oropharynx (SCCO) treated with RT, chemotherapy, and/or surgery, 22 and FAS c.-1378G > A altered OS of chemoradiation-treated HPV16-positive SCCO patients; 23 it is important to note that the chemotherapeutic agents used in both studies were not specified.
The aim of this study was to evaluate the roles of TP53 c.215G > C, FAS c.-671A > G, FAS c.-1378G >A, FASL c.-844 C > T, CASP3 c.-1191A > G, and CASP3 c.-182-247G > T SNVs in toxicity, response rate, and survival of HNSCC patients treated with a homogeneous protocol consisting of RT and CDDP chemotherapy.
Materials and methods
Study population
In this prospective study, we evaluated 109 patients with HNSCC diagnosed at the Clinical Oncology Service of the General Hospital of University of Campinas between June 2011 and February 2014. The study was conducted according to the Declaration of Helsinki and was approved by the local Ethics Committee (number 274/2011), where all subjects provided written informed consent.
Patients were classified as smokers or non-smokers and drinkers or abstainers. 24 HNSCC was diagnosed 25 and staged 26 according to standard criteria.
The presence of HPV16 in tumor fragments embedded in paraffin was analyzed by immunohistochemical staining and in situ hybridization as previously described.27,28
Patients were treated with chemoradiation as definitive treatment due to locoregional unresectable tumor, refusal of surgery related to expected functional or anatomic sequels, or an organ preservation protocol. The treatment consisted of 35 sessions of RT being 2 Gy per session and intravenous CDDP at dose of 80–100 mg/m2 on days 1, 22, and 43 (100 mg/m2 of CDDP was administered to patients with Karnofsky Performance Scale (KPS) 80%–100% and without comorbidities, and patients with KPS 60%–70% and without comorbidities or KPS higher than 70% with comorbidities received CDDP at dose of 80 mg/m2). Patients who presented adverse side effects of grades 3 or 4 received CDDP at lower dose in further administrations or had suspension of CDDP.29,30 As antiemetic premedication, mannitol (125 ml of 20%), ondansetron (24 mg), and dexamethasone (20 mg) were dissolved in 3000 ml of saline 0.9% and infused before CDDP infusion; and oral dexamethasone (8 mg, every 12 h for 3 days) and metoclopramide (10 mg, every 6 h for 3 days) after each CDDP infusion.31,32 Adherence to prescriptions was considered in study. 33
Vomiting was accessed after CDDP infusion and during the four following days. Cytopenias were evaluated by hematological exams performed after CDDP infusion. Nephrotoxicity and ototoxicity were analyzed using 51Cr-EDTA glomerular filtration rate and audiometric exam, respectively, performed pre and after chemoradiation. Toxicities were graded using the National Cancer Institute (NCI) criteria version 4.0, 34 and the worst toxicity grade during treatment was considered for analysis.
The response to treatment was quantified by Response Evaluation Criteria in Solid Tumors (RECIST) 1.1. 35 Patients who failed to respond to their initial treatment regimen or relapsed received intravenous methotrexate as palliative chemotherapy. 36 Patients were followed at 3-month intervals and the latest follow-up was recorded in December 2018.
Genetic polymorphism analysis
TP53 c.215G > C (rs1042522), 37 FAS c.-671A > G (rs1800682), 38 FAS c.-1378G > A (rs2234767) and FASL c.-844 C > T (rs763110), 39 CASP3 c.-1191A > G (rs12108497), 15 and CASP3 c.-182-247G > T (rs4647601) 40 genotypes were identified in genomic DNA by polymerase chain reaction and enzymatic digestion. Briefly, the regions of interest for each gene were amplified using the mix of sterile water, 10× buffer with (NH4)2SO4 (Thermo Scientific™, USA), 25 mM MgCl2 (Thermo Scientific™, USA), 40 mM dNTP (Thermo Scientific™, USA), 10 pmoles/µL of each primer (Integrated DNA Technologies™, USA), 2.5 U Taq DNA polymerase (Thermo Scientific™, USA), and 200 ng DNA. The amplification products were subjected to fragmentation with 5 U of each enzyme: BstUI (rs1042522 and rs2234767) (Thermo Scientific™, USA), MvaI (rs1800682) (Thermo Scientific™, USA), BsrDI (rs763110) (Thermo Scientific™, USA), StuI (rs12108497) (Thermo Scientific™, USA), and Hpych4 V (rs4647601) (New England Biolabs®, USA). Subsequently, the fragments were visualized on agarose gel stained with ethidium bromide. Positive and negative controls were used in reactions. The amount of 15% of genotype determinations was carried out twice in independent experiments with 100% of concordance.
Statistical analysis
The Haploview 4.2 software (www.broad.mit.edu/mpg/haploview) was used to perform pairwise LD analyses of FAS and CASP3 haplotypes, and those with frequency higher than 10% were selected to analyses. Differences between groups were analyzed by chi-square or Fisher’s exact test. To obtain adjusted odds ratio (OR) values in associations among SNVs, toxicity, and response to treatment, the logistic regression model was used. Bonferroni correction was applied to correct for bias arising out of multiple comparisons. The statistical power of analysis (PA) was calculated using the Researcher’s Toolkit. 41
Event-free survival (EFS) and OS were calculated from date of diagnosis to date of progression, relapse, death from disease or last follow-up, and date of diagnosis until date of death by any cause or lost to follow-up, respectively. EFS and OS were calculated using the Kaplan–Meier method. The prognostic impact of clinic pathological and genetic variants was examined using univariate Cox proportional hazard ratio (HR) regression. All variables with P-value ≤ 0.15 were included in multivariate Cox regression analysis. Significant results were internally validated conducting on the final model via the bootstrap method using 1000 replications.
Significance results were two-sided and achieved when P-value was < 0.05. All tests were done using the SPSS 24.0 software (SPSS Incorporation, IL, USA).
Results
The patients’ median age was 56 years old. Most patients were males, smokers, and drinkers; the three non-smoker patients referred accentuated alcohol consumption and the nine non-drinker patients were important smokers. About 60% of patients had tumor in oral cavity or oropharynx and more than 80% of patients showed well or moderately differentiated tumor and tumors of III or IV stage. All analyzed cases tested negative for HPV16. Vomiting was reported in about 60% of cases, anemia was identified in about 97% of cases, leukopenia in 72% of cases, neutropenia in 60% of cases, lymphopenia in 97% of cases, and almost 30% of patients presented thrombocytopenia during or after treatment. About one-quarter of the available patients obtained a complete response with chemoradiation, while the remainder presented partial or stable disease (Table 1).
Clinical and tumor aspects of 109 patients with head and neck squamous cell carcinoma.
N: number of patients.
The number of patients differed from the total quoted in the study (N = 109) because it was not possible to obtain consistent information in some cases. bGrades 1–4 for myelotoxicity and ototoxicity and 1–5 for nephrotoxicity.
All patients received the total dose of RT (70 Gy). Eighty-five patients received three CDDP infusions and 24 patients were treated with two CDDP infusions due to the occurrence of myelotoxicity and/or nephrotoxicity. Median cumulative CDDP amount administered to patients was 260 mg/m2 (range: 100–300 mg/m2). Medium or high adherence (97.7%) to antiemetics was observed among patients.
The median follow-up time of the 109 patients was 22 months (range: 3–89). The estimated probabilities of 24-month EFS and OS were 41.7% and 47.2%, respectively. At the study end, 34 patients were alive (6 of them were with HNSCC and 28 without HNSCC) and 75 patients died (67 of them by disease and 8 by others causes).
Association between genetic variants, toxicity, and response to therapy
When genotypes of SNVs were analyzed in patients stratified by bone marrow toxicity to chemoradiation, it was found an excess of FAS c.-671GG genotype in patients with anemia grades 2– 4 than in those without anemia or with anemia grade 1 (30.4% versus 14.3%, P = 0.01; PA: 45.8%) after Bonferroni correction (Figure 1(a)); patients with FAS c.-671GG genotype had 5.52 more chance of presenting anemia grades 2–4 than those with FAS c.-671AA or AG genotype. FAS c.-671AG or GG genotype was also more common in patients with lymphopenia grade 3 or 4 than in those without lymphopenia or with lymphopenia grade 1 or 2 (80.0% versus 56.3%; P = 0.005; PA: 71.9%) after Bonferroni correction (Figure 1(b)); patients with FAS c.-671AG or GG genotype had 4.03 more chance of presenting lymphopenia grade 3 or 4 than those with FAS c.-671AA genotype (Table 2).

Bars of genotypes frequencies in (a) patients without anemia or with anemia grade 1 (gray bars) versus patients with anemia grades 2–4 (black bars), (b) patients without lymphopenia or with lymphopenia grade 1 or 2 (gray bars) versus patients with lymphopenia grade 3 or 4 (black bars), and (c) patients without ototoxicity (gray bars) versus with patients grades (1–4) ototoxicity (black bars). *FAS c.-671GG genotype was more common in patients with anemia grades 2–4 than in those without anemia or with anemia grade 1 (30.4% versus 14.3%, P = 0.01; PA: 45.8%), FAS c.-671AG or GG genotype was more common in patients with lymphopenia grade 3 or 4 than in those without lymphopenia or with lymphopenia grade 1 or 2 (80.0% versus 56.3%; P = 0.005; PA: 71.9%), and FASL c.-844CC genotype was more common in patients with ototoxicity (38.2% versus 9.5%, P = 0.02; power of analysis: 75.4%); patients with the respective genotypes had 5.52 more chance of presenting anemia grades 2–4, 4.03 more chance of presenting lymphopenia grade 3 or 4, and 5.77 more chance of presenting ototoxicity compared with patients with other genotypes. Multivariate Cox regression curve for (d) event-free and (e) overall survival in patients with FAS c.-671AA or AG genotype (N = 83) versus GG genotype (N = 26), (f) event-free and (g) overall survival in patients with FAS c.-1378 GG or GA genotype (N = 102) versus AA genotype (N = 7), and (h) event-free and (i) overall survival in patients with FASL c.-844CC or CT genotype (N = 82) versus TT genotype (N = 27), adjusted by gender and tumor stage.
TP53 c.215G > C, FAS c.-671A > G and c.-1378G > A, FASL c.-844 C > T, CASP3 c.-1191A > G and c.-182-247G > T single nucleotide variants and FAS and CASP3 haplotypes in 109 patients with head and neck squamous cell carcinoma stratified by myelotoxicity with treatment with chemoradiotherapy.
G: grade of toxicity; N: number of patients; OR: odds ratio; CI: confidence interval; NE: not evaluated.
ORs were adjusted by gender to lymphopenia, and by age and tumor stage to anemia. The total number of patients differed from the total quoted in the study (N = 109) because it was not possible to obtain consistent information about toxicities in some cases.
Power of analysis (PA): 45.8%; bPA: 72.4%; cPA: 56.1%; dPA: 56.0%; ePA: 71.9%; fPA: 51.7%; gPA: 51.7%. a and eSignificant even after Bonferroni correction for multiple comparisons (corrected P-value = 0.01).
The SNVs analyzed in this study did not alter the occurrence of vomiting, response to CDDP chemoradiation, and nephrotoxicity (Supplementary Table 1). However, an excess of FASL c.-844CC genotype was found in patients with ototoxicity (all grades) than in those without ototoxicity (38.2% versus 9.5%, P = 0.02; PA = 75.4%) (Figure 1(c)); patients with FASL c.-844CC genotype had 5.77 more chance of presenting ototoxicity after therapy than those with FASL c.-844CT or TT genotype.
Association between genetic variants and survival
At 24 months of follow-up, lower EFS and OS were observed in patients with stage III or IV tumor (38.8% versus 83.3%; P = 0.02) and (44.1% versus 83.3%; P = 0.01), respectively, compared with those with tumor at stage I or II (Kaplan–Meier estimates). The variable mentioned above remained predictors of shorter EFS and OS in HNSCC patients in univariate analysis. Multivariate analysis showed that tumor at III or IV stage (HR: 9.77, P = 0.02; Pbootstrap = 0.02) and (HR: 9.74, P = 0.02; Pbootstrap = 0.01), and FAS c.-671GG genotype (HR: 2.05, P = 0.007, Pbootstrap = 0.02) and (HR: 1.83, P = 0.02, Pbootstrap = 0.045) were independent predictors of shorter EFS and OS, respectively (Table 3, Figure 1(d) and (e)). EFS and OS were not altered by TP53 c.215G > C, FAS c.-1378G > A (Figure 1(f) and (g)), FASL c.-844 C > T (Figure 1(h) and (i)), CASP3 c.-1191A > G, and CASP3 c.-182-247G > T SNVs.
Association of clinicopathological aspects and P53 c.215G > C, FAS c.-671A > G and c.-1378G > A, FASL c.-844 C > T, CASP3 c.-1191A > G and c.-182-247G > T genotypes and FAS and CASP3 haplotypes of 109 head and neck squamous cell carcinoma patients with event-free survival and overall survival.
N: number, HR: hazard ratio, CI: confidence interval, NE: not evaluated.
P-values ≤ 0.15 are presented in bold letters for univariate Cox analysis.
In multivariate Cox analysis adjusted by gender and tumor stage: aPbootstrap = 0.02; bPbootstrap = 0.02; cPbootstrap = 0.01; dPbootstrap = 0.04. The total number of patients differed from the total quoted in the study (N = 109) because it was not possible to obtain consistent information about histological grade in some cases.
Discussion
We investigated, in this prospective study, the roles of TP53 c.215G > C, FAS c.671A > G and c.1378G > A, FASL c.-844 C > T, CASP3 c.-1191A > G and c.-182-247G > T SNVs, and FAS and CASP3 haplotypes in adverse events, response to therapy, and survival of HNSCC patients treated with RT and concomitant CDDP.
We initially observed that clinicopathological aspects of our patients were in general like those seen in patients from other parts of the world,29,42–45 which indicate that our sample was representative of HNSCC in medical practice settings. The more consistent difference between our cases and others seemed to occur in HNSCC carcinogenesis. Although tobacco and alcohol consumption and HPV16 infection were significant risk factors for HNSCC worldwide, 4 HPV was not identified in our patients and was uncommon in those from other regions of the country.45,46
Second, we found that patients with FAS c.-671GG genotype and FAS c.-671AG or GG genotype had 5.52 and 4.30 more chances of presenting anemia grades 2–4 and lymphopenia grade 3 or 4 with CDDP chemoradiation than those with the remaining genotypes, respectively. Myelotoxicity was not altered by TP53 c.215G > C, FAS c.1378G > A, FASL c.-844 C > T, CASP3 c.-1191A > G and c.-182-247G > T SNVs. Khrunin et al. 17 found association of TP53 c.215CC genotype with severe neutropenia in ovarian cancer patients treated with CDDP-cyclophosphamide. The difference in results obtained in our study and Khrunin’s study may be attributed to distinct types of tumors and treatments. To the best of our knowledge, the role of FAS c.-671A > G SNV in cytopenias induced by CDDP chemoradiation was not previously analyzed.
It is already known that CDDP may induce a deep transient erythropoiesis alteration leading to anemia.47,48 CDDP inhibits bone marrow cell proliferation by DNA damage, 49 possibly due to apoptosis induced by FAS activation. 9 Because the allele G of FAS c.671A > G was previously associated with lower apoptosis of damaged cells, 13 these results were not expected by us. However, activation of FAS/FASL was described as an important mechanism in erythroid maturation 50 and regulation of red blood cell homeostasis, 51 and Rehman et al. 52 observed that G allele of FAS c.-671A > G increased aplastic anemia risk. Moreover, deletion of FAS in lymphocytes caused lymphopenia in mice, 53 suggesting that FAS expression in lymphocytes is required for their survival and/or proliferation. 54 So, we hypothesize that abnormalities in erythroid maturation and lymphocytes survival and/or proliferation may have contributed to the moderate/severe anemia and lymphocitopenia seen in our patients with FAS c.-671GG genotype under CDDP chemoradiation.
We also observed that patients with FASL c.-844CC genotype had 5.77 more chance of presenting ototoxicity after CDDP chemoradiation than those with the remaining genotypes. So far, our knowledge reaches, there are no previous studies focusing on the role of FASL c.-844 C > T in ototoxicity induced by RT and CDDP.
Such association is biologically plausible because FASL c.-844CC genotype was related to higher apoptosis capacity, 14 CDDP induces apoptosis in cells of the cochlea by generating reactive oxygen species (ROS), 55 and ROS can promote programmed cell death from activation of the FAS/FASL pathway in auditory hair cell. 56
Third, no association of TP53 c.215G > C, FAS c.-671A > G and c.1378G > A, FASL c.-844 C > T, CASP3 c.-1191A > G and c.-182-247G > T with response to treatment was seen in this study. Shiraishi et al. 20 analyzed association of TP53 c.215G > C SNV with response to platinum plus vinorelbine, paclitaxel, docetaxel, or gemcitabine in lung cancer patients, and observed that patients with the CC genotype obtained a better response than others. Again, the distinct types of tumor and treatments may explain differences of results found in both studies.
Fourth, we observed that HNSCC patients with FAS c.-671GG genotype had lower EFS and OS, respectively, than those with the remaining genotypes, but TP53 c.215G > C, FAS c.1378G > A, FASL c.-844 C > T, CASP3 c.-1191A > G and c.-182-247G > T SNVs did not alter patients’ survival. Zhang et al. 22 evidenced associations of FAS c.-671AG or GG and FASL c.-844CT or TT genotypes with lower DFS in SCCO patients treated with RT, chemotherapy, and/or surgery, and the associations were particularly seen in patients with HPV16-positive tumors. No significant associations between FAS c.1378G > A SNP and patients’ survival were found in their study. Shiraishi et al. 20 found longer DFS and OS in lung cancer patients with CC genotype of TP53 c.215G > C treated with platinum plus vinorelbine, paclitaxel, docetaxel, or gemcitabine, and Azad et al. 21 observed lower DFS in early-stage HNSCC patients radiation-treated with the CC genotype of TP53 c.215G > C compared with those with the remaining genotypes. FAS c.1378AA genotype was associated with higher OS in chemoradiation-treated HPV16-positive SCCO patients analyzed by Feng et al. 23 Again, the opposed results obtained in our study and previously reported studies may be attributed to the different types of tumors, carcinogenesis, or treatments.
Since the FAS c.-671GG genotype was related to lower apoptosis of cancer cells, 13 the association of genotype with lower EFS and OS in HNSCC patients is also biologically plausible, and probably results from tumor cells survival and proliferation.
At this time, it is important to comment that although Zhang et al. 22 had seen association of FAS c.-671AG or GG genotype with lower DFS particularly in HPV16-positive SCCO patients, we found associations of the GG genotype with DFS and OS of HNSCC, which were strongly related to tobacco smoking and alcohol consumption.
Finally, we postulate that although we have not evaluated combinations of genotypes of distinct genes, we cannot exclude co-participation of TP53 c.215G > C, CASP3 c.-1191A > G, and CASP3 c.-182-247G > T SNVS in associations between FAS and myelotoxicity and patients’ survival and FASL and ototoxicity found in our study.
We are aware that our study has some inherent limitations, such as the relatively small number of patients and the lack of investigation of HPV16 in all tumor samples, which may have limited the obtaining and interpretation of the findings.
In conclusion, our data present, for the first time, associations of FAS c.-671A > G and FASL c.-844 C > T with hematological toxicity, ototoxicity, and survival of tobacco- and alcohol-related HNSCC patients homogeneously treated with CDDP chemoradiation. Nevertheless, large studies with evaluation of HPV-16 in tumor fragments are needed to validate our results and further explore the molecular mechanisms that underlie the observed associations, which may have future utility as clinical prognostic biomarkers.
Supplemental Material
Supplementary_Material – Supplemental material for FAS and FASL variations in outcomes of tobacco- and alcohol-related head and neck squamous cell carcinoma patients
Supplemental material, Supplementary_Material for FAS and FASL variations in outcomes of tobacco- and alcohol-related head and neck squamous cell carcinoma patients by Ericka Francislaine Dias Costa, Tathiane Regine Penna Lima, Leisa Lopes-Aguiar, Guilherme Augusto Silva Nogueira, Marília Berlofa Visacri, Julia Coelho França Quintanilha, Eder Carvalho Pincinato, Luciane Calonga, Fernanda Viviane Mariano, Albina Messias de Almeida Milani Altemani, João Maurício Carrasco Altemani, Patrícia Moriel, Carlos Takahiro Chone, Celso Dario Ramos and Carmen Silvia Passos Lima in Tumor Biology
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by State of São Paulo Research Foundation (FAPESP) (grant numbers 2012/01807-2, 2012/01418-6).
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
