Next Article in Journal
Real World Outcomes versus Clinical Trial Results of Durvalumab Maintenance in Veterans with Stage III Non-Small Cell Lung Cancer
Next Article in Special Issue
Pharmacogenomics Testing in Phase I Oncology Clinical Trials: Constructive Criticism Is Warranted
Previous Article in Journal
The Central Role of the Ubiquitin–Proteasome System in EBV-Mediated Oncogenesis
Previous Article in Special Issue
The Tumour Suppressor CYLD Is Required for Clathrin-Mediated Endocytosis of EGFR and Cetuximab-Induced Apoptosis in Head and Neck Squamous Cell Carcinoma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 14
Issue 3
10.3390/cancers14030612
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Systematic Review

Pharmacogenomics of Vincristine-Induced Peripheral Neuropathy in Children with Cancer: A Systematic Review and Meta-Analysis

by
Aniek Uittenboogaard
1,2,*,
Céline L. G. Neutel
3,
Johannes C. F. Ket
4,
Festus Njuguna
5,
Alwin D. R. Huitema
2,6,7,
Gertjan J. L. Kaspers
1,2,*,† and
Mirjam E. van de Velde
1,†
1
Emma Children’s Hospital, Amsterdam UMC, Vrije Universiteit Amsterdam, Pediatric Oncology, 1105 AZ Amsterdam, The Netherlands
2
Princess Máxima Center for Pediatric Oncology, 3584 CS Utrecht, The Netherlands
3
Department of Neurosurgery, Radboud University Medical Center, 6525 GA Nijmegen, The Netherlands
4
Medical Library, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
5
Department of Pediatric Oncology, Moi University, Eldoret 30107, Kenya
6
Department of Pharmacy & Pharmacology, The Netherlands Cancer Institute, 1066 CX Amsterdam, The Netherlands
7
Department of Clinical Pharmacy, University Medical Center Utrecht, Utrecht University, 3584 CX Utrecht, The Netherlands
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cancers 2022, 14(3), 612; https://doi.org/10.3390/cancers14030612
Submission received: 29 November 2021 / Revised: 20 January 2022 / Accepted: 24 January 2022 / Published: 26 January 2022
(This article belongs to the Special Issue Pharmacokinetics/Pharmacodynamics/Pharmacogenetics Based Personalized Medicine for Cancer Therapy)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Vincristine is a drug that is part of the treatment for many children with cancer. Its main side-effect is vincristine-induced peripheral neuropathy (VIPN), which often presents as tingling, pain, and lack of strength in the hands and feet. It is not yet possible to predict which children will suffer from VIPN. In this review, we report on all genetic variations that are associated with VIPN. We found that variations in genes related to vincristine transport, cell structure, hereditary nerve disease, and genes without a previously known connection to vincristine or VIPN are related to VIPN. Variations in genes involved in vincristine breakdown are not significantly associated with VIPN. In conclusion, genetic variations affect a child’s tendency to develop VIPN. In the future, this information might be used to predict the risk of VIPN and adapt treatment on this.

Abstract

Vincristine-induced peripheral neuropathy (VIPN) is a debilitating side-effect of vincristine. It remains a challenge to predict which patients will suffer from VIPN. Pharmacogenomics may explain an individuals’ susceptibility to side-effects. In this systematic review and meta-analysis, we describe the influence of pharmacogenomic parameters on the development of VIPN in children with cancer. PubMed, Embase and Web of Science were searched. In total, 1597 records were identified and 21 studies were included. A random-effects meta-analysis was performed for the influence of CYP3A5 expression on the development of VIPN. Single-nucleotide polymorphisms (SNPs) in transporter-, metabolism-, cytoskeleton-, and hereditary neuropathy-associated genes and SNPs in genes previously unrelated to vincristine or neuropathy were associated with VIPN. CYP3A5 expression status was not significantly associated with VIPN. The comparison and interpretation of the results of the included studies was limited due to heterogeneity in the study population, treatment protocol and assessment methods and definitions of VIPN. Independent replication is essential to validate the clinical significance of the reported associations. Future research should aim for prospective VIPN assessment in both a discovery and a replication cohort. Ultimately, the goal would be to screen patients upfront to determine optimal vincristine dosage with regards to efficacy and risk of VIPN.

1. Introduction

Vincristine is an important chemotherapeutic agent that is commonly used in treatment for pediatric cancers. It is approved by the United States Food and Drug Administration (FDA) for the treatment of acute lymphoblastic leukemia (ALL), Hodgkin and non-Hodgkin lymphoma, neuroblastoma, rhabdomyosarcoma, low-grade glioma and nephroblastoma. Furthermore, off-label uses include the treatment of Ewing sarcoma and medulloblastoma [1,2]. The main side-effect of vincristine is vincristine-induced peripheral neuropathy (VIPN), which often presents as a symmetric sensory-motoric neuropathy progressing distally to proximally [1,2]. Presenting signs include foot drop, loss of deep tendon reflexes, impaired balance, pain or tingling [1,2]. In addition, patients can suffer from autonomic symptoms such as constipation or orthostatic hypotension. The reported prevalence of VIPN varies, depending on assessment method and study population, but it is estimated that the majority of patients receiving vincristine will experience some form of VIPN during treatment [1,2,3,4]. Up to 30% of patients may suffer from severe VIPN, requiring dose reduction or cessation of treatment [3,5]. Suffering from VIPN is associated with a lower health-related quality of life, both by self- and proxy assessment and consistently when using different assessment tools for VIPN [6]. This effect of VIPN on health-related quality of life seems to persevere after treatment, as was shown in a recent study in ALL survivors in which over 16% suffered from long-term VIPN and experienced impact on both physical health and social functioning [7].
It is recognized that different populations might have an altered risk for VIPN [3]. Older age has been associated with an increased risk of VIPN, although results have been inconsistent [8,9,10,11,12]. In addition, white children appear to have a higher risk of VIPN than black children [3,9,12,13,14,15], which is corroborated by a recent study in Kenyan pediatric cancer patients in which only one out of 78 black patients developed severe VIPN and less than 5% developed clinically relevant VIPN, despite the use of sensitive assessment methods [16]. Interestingly, these children are being treated at a higher vincristine dose than what is common in Western countries (2.0 mg/m2 as opposed to 1.5 mg/m2) [1,16]. Studies assessing the relationship between VIPN and vincristine pharmacokinetics (PK) have shown inconsistent results. Some studies show a correlation between VIPN and PK parameters such as area under the curve (AUC) [17], an estimate of vincristine exposure, and intercompartmental clearance [18], whereas others do not confirm these findings [19,20,21,22]. Therefore, potential risk factors for VIPN could be genetic variations in genes involved in vincristine PK, such as variations in the cytochrome (CYP) 450 family of enzymes. Vincristine is predominantly metabolized by CYP3A4 and CYP3A5, of which the latter has a higher intrinsic clearance [23]. Genetic variants in both enzymes result in different metabolic activity [23,24]. Racial populations have different distributions of wild-type and variant CYP3A4/5 alleles [25,26,27]. Combined with the observation that black patients develop less VIPN, it has led to the hypothesis that faster clearance of vincristine in black children results in a lower risk of VIPN in comparison to white patients [14]. Indeed, several studies have described the effect of variations in CYP3A4 and CYP3A5 on the development of VIPN [8,13,14,16,20,28,29,30,31,32]. Differences in VIPN prevalence across populations may thus stem from variations in genetic background, which can be studied via the rapidly expanding field of pharmacogenomics.
Pharmacogenomics aims to assess the influence of genomics on an individuals’ treatment response and susceptibility to side-effects, such as VIPN [33,34]. Often, the effect of single nucleotide polymorphisms (SNPs) is assessed [35,36]. The frequency distribution of major and minor alleles varies across racial groups and study populations, which has been well characterized in large projects such as the 1000 Genomes Project and the genome Aggregation Database (gnomAD) [37,38]. Pharmacogenomics aims to find those SNPs or genetic variations that are biologically relevant [35,36]. Two main study designs have been used to assess this: candidate gene studies or population-based genome- or exome-wide association studies (GWAS or EWAS respectively) [39,40]. Candidate gene studies determine, a priori, a set of genes, based on available literature or mechanism of action, whose influence on a certain outcome is to be assessed [39]. Population-based GWAS or EWAS, on the other hand, assess the whole exome or genome (by whole exome sequencing (WES) or whole genome sequencing (WGS)) for genetic variation in relation to a certain outcome measure [39]. These studies may therefore result in previously unknown genotype—phenotype associations.
Pharmacogenomics can serve as a guidance tool for precision therapy in which a priori a patients’ genetic susceptibility for side-effects or therapeutic efficacy is determined. Although this has been implemented in clinical practice for some drugs, such as thiopurine methyltransferase (TPMT), this is currently not possible for vincristine [41,42]. Especially since there is a lack of understanding of what causes variability in VIPN across patients, pharmacogenomics can provide valuable insight into the pathogenesis of VIPN. If genes affecting vincristine PK are implicated, this may emphasize the potential of therapeutic drug monitoring. Moreover, since it is unlikely that VIPN is caused by differences in PK alone, variation in cellular sensitivity to vincristine and in neuronal pathways could be contributing factors. The implication of genes related to neuronal pathways, the cytoskeleton or cellular integrity with VIPN might then help guiding clinicians in deciding a priori if patients have a high chance of being developing (clinically relevant) VIPN and thus if patients should be monitored more closely than others, or even given an adapted vincristine dosage. In contrast, other patients might be identified who tolerate higher levels of vincristine and might thus not benefit from the generally applied dose capping at 1.5 mg/m2. Ultimately, the goal would be to develop a protocol for vincristine in which patients are stratified based on the presence of genetic polymorphisms and given a dosage that limits the risk of severe VIPN while maintaining the highest possible therapeutic efficacy. However, to explore this possibility, the first step is to provide a detailed overview of the effect of SNPs in all reported genes so far on VIPN. Therefore, in this systematic review, we aim to describe the influence of pharmacogenomic parameters on the development of VIPN in children with cancer. Furthermore, we performed a meta-analysis on the influence of CYP3A5 expression status on the development of VIPN.

2. Materials and Methods

2.1. Protocol and Registration

This study was performed according to the Preferred Reporting Items for Systematic Reviews and Meta-analyses guidelines (PRISMA) [43]. The study protocol was registered at the PROSPERO International prospective register of systematic reviews (registration number CRD42021210437) [44].

2.2. Eligibility Criteria

We included prospective and retrospective case-control or cohort studies assessing the relation between VIPN and pharmacogenomic parameters in five or more children with cancer. No restrictions regarding the characteristics of the children with cancer were applied. Pharmacogenomic parameters could include RNA variations such as microRNAs (miRNA) as well. To make sure we described the relationship between VIPN and pharmacogenomic parameters, we excluded descriptive reviews and studies in which no distinction could be made whether patients suffered from VIPN or neuropathy due to other causes, such as diabetes or Charcot-Marie-Tooth disease. If necessary, authors were contacted for clarification or additional data. Systematic reviews and meta-analyses were screened for additional inclusions.

2.3. Information Sources

PubMed, Embase, and Clarivate Analytics/Web of Science Core Collection were searched from inception up to 30 September 2021 (by AU and JCFK). Search terms were used as controlled vocabulary (e.g., MeSH) as well as free text terms for title, abstract and author keywords were amongst others: ‘pharmacogenetics’ and ‘child’ and ‘vincristine’ and ‘neuropathy’ or ‘constipation’. No limitations on language or date were applied. See Supplementary Materials 1 for the full search strategies in all databases.

2.4. Study Selection

Title and abstract screening were performed independently by two reviewers (AU and CLGN) based on pre-defined in- and exclusion criteria. Next, studies were screened full-text for eligibility by two reviewers (AU and CLGN). In case of disagreement, a third reviewer was consulted (MEvdV). Reasons for exclusion after full-text screening were documented. A meta-analysis was performed for studies assessing the relationship between CYP3A5 expression status and VIPN.

2.5. Risk of Bias Assessment

Risk of bias was assessed by two independent reviewers (AU and CLGN) according to a modified version of the quality assessment tool for quantitative studies by the Effective Public Health Practice Project (EPHPP), scoring each study as strong, moderate or weak in seven domains: study design, confounders, blinding, data collection methods, withdrawals and drop-outs, analysis, and selection of reported results [45]. The global rating of each study was determined as follows: no weak ratings of the domains resulted in an overall strong rating, one weak rating of the domains resulted in an overall moderate rating, and two or more weak ratings of the domains resulted in an overall weak rating.

2.6. Data Extraction and Synthesis

For the included studies, study and baseline characteristics were extracted according to a data extraction template (Supplementary Materials 2). The main outcome was the relationship between VIPN and pharmacogenomic parameters expressed as an effect size (odds ratio (OR)) or p-value. ORs were either reported by studies or calculated by the authors. If studies reported multiple ORs using different definitions of cases of VIPN and controls, the significant result with highest clinical relevance was shown (severe VIPN or any grade VIPN). If studies performed both univariate and multivariate analyses, both effect sizes were included in this study.
For the meta-analysis on CYP3A5 expression status and VIPN, dominant OR was calculated if raw data were available, or pre-calculated ORs provided by the authors were used. CYP3A5 expressers were defined as having at least one functional allele (*1) and CYP3A5 non-expressers were defined as having only non-functional variant alleles. If needed, authors were contacted for missing data.

2.7. Statistical Analysis

The meta-analysis on CYP3A5 was performed in R, version 3.6.1, using the ‘meta’ package (Rstudio Inc.) [46]. A random-effects model was applied to pool odds ratios since considerable between-study heterogeneity was suspected. The Paule–Mandel procedure was used to estimate variance τ2. Heterogeneity was estimated using I2 with the following interpretations: I2 of 25% indicated low heterogeneity, I2 of 50% indicated moderate heterogeneity, and I2 of 75% indicated high heterogeneity. The inverse-variance approach was used with the ‘metagen’ command. If study heterogeneity I2 was higher than 50%, an assessment of outliers or influential cases was performed. The risk of publication bias was assessed via evaluation of a funnel plot and Egger’s testing for asymmetry. A two-sided p-value of <0.05 was considered statistically significant.

3. Results

3.1. Study Selection

We identified 1597 reports through database searching (Figure 1). After removal of duplicates, the abstracts and titles of 1367 records were screened. Of these, 109 were selected for full-text review. One report was sought for retrieval and no full-text version was published. Since insufficient data was available in the abstract, the report was excluded. Assessment of eligibility resulted in the exclusion of 88 reports based on: consisted of a narrative review (35 reports), no data available on pharmacogenomic parameters (24 reports), no data available on VIPN (13 reports), no administration of vincristine (12 reports), and same data were used as in another report (four reports). Finally, 21 reports were included in this systematic review.

3.2. Study Characteristics

Sociodemographic and baseline characteristics of the included studies can be found in Table 1. Eighteen studies followed a candidate gene approach [8,10,11,13,14,16,20,22,28,29,30,31,32,47,48,49,50,51], whereas three studies were population-based GWAS or EWAS [9,52,53] (Table 1). Four studies included a replication cohort to confirm their findings of the discovery cohort [8,9,52,53]. In total, the number of included patients with available genotype and VIPN data ranged from 24 to 1132. The majority of included patients were diagnosed with ALL and were white (Table 1). The prevalence of moderate to severe VIPN (grade 2–4) ranged from 19.5–53.2%, with the exception of the study by Skiles et al. who reported an incidence of 2.8% in black Kenyan patients [16]. Studies used different definitions of cases (patients with VIPN) and controls (patients without VIPN) (Table 2). Different measurement tools for VIPN were used in the different studies, most often the Common Terminology Criteria for Adverse Events (CTCAE), in which a subset of items was used to score peripheral neuropathy [8,9,10,13,14,16,20,29,30,47,50,51,52,53], followed by the modified Balis scale [9,16,31,32] and World Health Organization (WHO) scale [11,48,49] were used. Finally, The Children’s Cancer Group (CCG) toxicity criteria, National Cancer Institute (NCI) common toxicity criteria, Total Neuropathy Score—Pediatric Vincristine (TNS-PV) and pediatric modified total neuropathy score (ped-mTNS) were all used in one study each [16,22,28,53]. Seven studies assessed VIPN prospectively [9,16,20,28,32,47,53], while the rest of the studies assessed VIPN retrospectively.

3.3. Risk of Bias

An overview per domain of risk of bias can be found in Table S1. Twelve [8,10,11,22,29,30,31,48,49,50,51,52] and nine studies [9,13,14,16,20,28,32,47,53] scored an overall strong and moderate rating on risk of bias, respectively. No studies received an overall weak rating on risk of bias.
Table
Table 1. Sociodemographic and clinical characteristics of studies included in the systematic review.
Table 1. Sociodemographic and clinical characteristics of studies included in the systematic review.
Author and Year of PublicationStudy DesignPatients with Genotype + VIPN Data (n)Patient CharacteristicsVincristine DosageVIPNGlobal Rating Risk of Bias Assessment
Disease StudiedAgeMale (%)Race (%)Single Dosage, (per mg/m2 and max)Cumulative Dosage (mg)Method Used for VIPN AssessmentPrevalence VIPN
Abaji et al., 2018—QcALL cohort [52]EWAS237ALL82.7% <10 y/o, 17.3% ≥10 y/o.54.9All white 1.5, max. 2.0 Not available NCI-CTCAE 3.0, retrospective
Grade 3–4 peripheral neuropathy		14.8%Strong
Abaji et al., 2018—AIEOP cohort [52]EWAS405ALL83.2% <10 y/o, 16.8% ≥10 y/o.53.1All white1.5, max. 2.0Not available NCI-CTCAE 3.0, retrospective
Grade 3–4 peripheral neuropathy		3.2%Strong
Abo–Bakr et al., 2017 [47]Candidate gene 97ALL79.4% ≤10 y/o, 20.6% >10 y/o58.8All white 1.5, max. 2.0Not available NCI-CTCAE 3.0, prospective
Foot drop, ileus, vocal cord paralysis, ptosis		Foot drop: 4.1%Moderate
Aplenc et al., 2003 [28]Candidate gene, case–control 533ALL70.0% ≤5 y/o, 30.0% >5 y/o32.55.8 black
94.2 other 		1.5, max. not available46.5–64.5CCG toxicity criteria, prospective
Grade 3 or 4 peripheral neuropathy		5.3%Moderate
Ceppi et al., 2014 [8]Candidate gene320ALL80.0% ≤10 y/o, 20.0% >10 y/o 55.3All white 1.5–2.0, max. 2.073.5–74.0NCI-CTCAE 3.0, retrospective
Peripheral neuropathy		Grade 1–2: 20.0%
Grade 3–4: 10.6%		Strong
Diouf et al., 2015—St. Jude cohort [9]GWAS St. Jude: 222. ALL68.9% ≤10 y/o, 31.1% >10 y/o 42.367.1 white, 19.8 black, 14.0 other 1.5, max. 2.0
COG: 1.5– 		54.0NCI-CTCAE 1.0, prospective
Grade 2–4 peripheral neuropathy		28.8%Moderate
Diouf et al., 2015—COG cohort [9]GWAS99Relapsed ALL47.5% ≤10 y/o, 52.5% >10 y/o 59.660.6 white, 1.0 black,
38.3 other 		1.5–2.0, max. 2.0–2.578.0–97.5Modified Balis scale, prospective
Grade 2–4 peripheral neuropathy		22.9%Moderate
Egbelakin et al., 2011 [29]Candidate gene107ALL Not availableNot available92.5 white
0.9 black
6.5 other 		1.5, max. 2.0 Not availableNCI-CTCAE 3.0, retrospective
Peripheral and autonomic neuropathy		Grade 1–4: 98.1%
Grade 3–4: 53.2%		Strong
Guilhaumou et al., 2011 [20]Candidate–gene24Solid tumors57.7% <10 y/o, 42.3% ≥10 y/o57.7All white 1.5, max 2.0Mean (SD) at time of enrolment: 7.35 (5.30)NCI-CTCAE 3.0, prospective
Pain, peripheral neuropathy, gastro–intestinal toxicity		33.3%Moderate
Gutierrez–Camino et al., 2016 [10]Candidate gene142ALL88.7% ≤10 y/o, 11.3% >10 y/o 57.0All white1.5, max 2.015.0–30.0NCI-CTCAE 1.0, retrospective
Grade 2–4 peripheral neuropathy		25.4%Strong
Gutierrez–Camino et al., 2017 [48]Candidate gene (miRNA) 155ALL Mean (SD): 5.1 (3.2) y/o 58.9Mainly white1.5, max 2.015.0–30.0WHO criteria, retrospective
Peripheral neuropathy		Grade 1–2: 16.0%
Grade 3–4: 10.1%		Strong
Kayilioğlu et al., 2017 [30]Candidate gene, case–controlCases: 115 (VCR), controls: 50 (no VCR) Cases: ALL and solid tumors.
Controls: no neurological disorders or symptoms		Mean (SD): ALL 7.0 (4.6), solid tumors 7.5 (5.0), controls 10.2 (4.6)ALL and solid tumors: 61.7
Controls: 62.0		All white 1.5, max 2.0Mean (SD) total: ALL 7.71 (0.89), solid tumors 6.5 (1.5) NCI-CTCAE 3.0, retrospective
Grade 2–5 neurotoxicity		20.8%Strong
Kishi et al., 2007 [13]Candidate gene240ALL70.4% ≤10 y/o, 29.6% >10 y/o 59.269.6 white
18.3 black
12.1 other 		1.5, max 2.054.0–97.5NCI-CTCAE 1.0, prospective/retrospective not available.
Peripheral neuropathy and constipation		Grade 3: 12.1%
Grade 4: 0.4%		Moderate
Li et al., 2019—POG cohort [53]GWAS 1069.ALLNot available52.3 All white1.5, max not available18–23 doses of 1.5 mg/m2NCI-CTCAE 2.0, prospective
Grade 3–5 peripheral neuropathy.		4.8%Moderate
Li et al., 2019—ADVANCE cohort [53]GWAS63ALLMean (SD): 8.2 (4.7) y/o46.0All white1.5, max 2.0Not availableTNS–PV, prospective.
Sensory symptoms, temperature and vibration sensibility, strength, tendon reflexes.		Mean + SD: 3.8 (2.6)Moderate
Lopez–Lopez et al., 2016 [11]Candidate gene 133ALLMean (SD): 5.5 (3.4) y/o56.6Mainly white1.5, max 2.015.0–30.0WHO criteria, retrospective
Peripheral neuropathy		Grade 1–2: 18.4%
Grade 3–4:
11.8%		Strong
Martin–Guerrero et al., 2019 [49]Candidate gene 133ALLMean (SD): 5.5 (3.4) y/o56.6Mainly white1.5, max 2.015.0–30.0WHO criteria, retrospective
Grade 2–4
peripheral neuropathy		25.4%Strong
McClain et al., 2018 [31]Candidate gene 239ALLMean (SD): 5.8 (3.9) y/o53.1All white Not availableMean (SD), at time of event: extensive metabolizers: 10.0 (5.7), intermediate: 13.4 (13.6), poor: 10.4 (8.9)Modified Balis scale, retrospective
Grade 3–4 peripheral neuropathy		Grade 3–4: 18.4%Strong
Plasschaert et al., 2004 [22]Candidate gene 52ALL73.1% < 10 y/o, 26.9% ≥ 10 y/o61.598.1 white
1.9 other 		Once 1.5, other doses 2.0, max. 2.513.5 mg/m2NCI common toxicity criteria
Constipation		Grade 1–2: 55.8%, Grade 3–4 26.9%Strong
Renbarger et al., 2008 [14]Race as surrogate for genotype, case–control Cases: 21 black
Controls: 92 white		ALL Mean (SD): black: 8.2 (4.8) y/o, white: 5.0 (3.1) y/oCases + controls: 50.481.4 white
18.6 black		Not availableMean (SD), Caucasians: 48.5 (14.3), AAs: 42.4 (11.6)NCI-CTCAE 3.0, retrospective
Neurotoxicity		Grade 1–4: 34.8% white, 4.8 black Moderate
Sims et al., 2016 [32]Candidate gene52BALL 77.4% < 10 y/o, 22.6% ≥ 10 y/o62.268.5 white
31.5 black		1.5, max. 2.0Not availableModified Balis scale, prospective
Peripheral neuropathy, constipation if grade 3–4		Grade 1–4: 80.6% white, 76.5% black Moderate
Skiles et al., 2018 [16]Candidate gene72Leukemia, lymphoma, solid tumors Mean (SD): low expressers: 6.1 (5.2), intermediate: 6.5 (4.0), high: 6.1 (4.6)53.8All black Kenyan2.0, max. 2.58.5 mg/m2NCI-CTCAE 4.0, modified Balis scale, Faces Pain Scale, Pediatric Neuropathic Pain Scale, ped–mTNS, all prospective.
Peripheral neuropathy and neuropathic pain		NCI–CTCAE: grade 2–4: 2.8%. Ped–mTNS: 4.3% 5 or higher.Moderate
Wright et al., 2019 [51]Candidate gene, case–control Cases: 167 (VIPN), controls: 57 (no VIPN)ALLMedian (IQR): cases 4.8 (3.3–9.0), controls: 5.4 (3.3–9.0)Cases: 60.4, controls: 40.4Mainly whiteNot availableMedian + IQR: cases: 61.4 (48.0–72.0), controls: 66.0 (51.0–74.8)NCI-CTCAE 4.0, retrospective
Peripheral neuropathy		Grade 2–4: 167 cases Strong
Zgheib et al., 2018 [50]Candidate gene 133ALLMean (SD): 6.7 (5.0)57.1All whiteInduction and re–induction: 1.5, max. 2.0. Continuation: 2.0, max. 2.0Mean (SD), patients without VIPN: 66.0 (6.1), with VIPN grade 2–4: 27.9 (12.1)NCI-CTCAE 4.0, retrospective
Peripheral neuropathy		Grade 2–4: 19.5%Strong
EWAS = exome-wide association study, ALL = acute lymphoblastic leukemia, NCI-CTCAE = National Cancer Institute—Common Toxicity Criteria for Adverse Events, CCG = Children’s Cancer Group, GWAS = genome-wide association study, SD = standard deviation, miRNA = microRNA, WHO = World Health Organization, TNS-PV = Total Neuropathy Score—Pediatric Vincristine, ped-mTNS = pediatric modified total neuropathy score, NCI = National Cancer Institute, IQR = interquartile range.
Table
Table 2. Single-nucleotide polymorphisms that were significantly associated with vincristine-induced peripheral neuropathy in the pediatric oncology population.
Table 2. Single-nucleotide polymorphisms that were significantly associated with vincristine-induced peripheral neuropathy in the pediatric oncology population.
GeneSNPAllele, Major/MinorAuthor and Year of PublicationMAF (%)Number of Patients (n) Method Effect Size Effect Size with 95% CI (If Applicable)Effect
Cases of VIPN *Controls *
Transport
ABCB1rs4728709C/TCeppi et al., 2014 [8]TT/TC: 17.1
CC: 82.9		63 (grade 1–2)214 (grade 0)Dominant OR 0.3 (0.1–0.9)Protective 1
rs10244266T/GLopez-Lopez et al., 2016 [11]14.346 (WHO grade 1–4)103 (WHO grade 0) Dominant OR 2.60 (1.16–5.83)Risk 2
rs10268314T/CLopez-Lopez et al., 2016 [11]14.327 (WHO grade 1–2)103 (WHO grade 0)Dominant OR 3.19 (1.23–8.25)Risk 2
rs10274587G/A Lopez-Lopez et al., 2016 [11]14.627 (WHO grade 1–2)103 (WHO grade 0) Dominant OR 3.48 (1.36–8.86)Risk 2
ABCC1rs1967120T/CLopez-Lopez et al., 2016 [11]27.318 (WHO grade 3–4)103 (WHO grade 0) Dominant OR 0.29 (0.09–0.99)Protective 2
rs3743527C/TLopez-Lopez et al., 2016 [11]19.746 (WHO grade 1–4)103 (WHO grade 0) Dominant OR 0.32 (0.13–0.79)Protective 2
rs3784867C/T Wright et al., 2019 [51]32.0170 (grade 2–4)57 (grade 0) Additive OR 4.91 (1.99–12.10)Risk 3
rs11642957T/CLopez-Lopez et al., 2016 [11]48.146 (WHO grade 1–4)103 (WHO grade 0)Dominant OR 0.43 (0.19–0.98)Protective 2
rs11864374G/ALopez-Lopez et al., 2016 [11]24.446 (WHO grade 1–4)103 (WHO grade 0)Dominant OR 0.35 (0.15–0.79)Protective 2
rs12923345T/CLopez-Lopez et al., 2016 [11]15.446 (WHO grade 1–4)103 (WHO grade 0)Dominant OR 2.39 (1.08–5.25)Risk 2
rs17501331A/GLopez-Lopez et al., 2016 [11]13.246 (WHO grade 1–4)103 (WHO grade 0)Dominant OR 2.50 (1.10–5.68)Risk 2
ABCC2rs12826G/ALopez-Lopez et al., 2016 [11]42.646 (WHO grade 1–4)103 (WHO grade 0)Dominant OR 0.24 (0.10–0.54)Protective
rs3740066G/ALopez-Lopez et al., 2016 [11]36.246 (WHO grade 1–4)103 (WHO grade 0)Dominant OR 0.23 (0.10–0.53)Protective
rs2073337A/GLopez-Lopez et al., 2016 [11]45.818 (WHO grade 3–4)103 (WHO grade 0)Dominant OR 0.35 (0.10–1.24)Protective
rs4148396C/TLopez-Lopez et al., 2016 [11]42.146 (WHO grade 1–4)103 (WHO grade 0)Dominant OR 0.36 (0.16–0.81)Protective
rs11190298G/ALopez-Lopez et al., 2016 [11]45.046 (WHO grade 1–4)103 (WHO grade 0)Recessive OR 2.44 (1.01–5.86)Risk
ABCC1/RALPB1: miR–3117rs12402181G/AGutierrez–Camino et al., 2017 [48]14.819 (WHO grade 3–4)128 (WHO grade 0)Dominant OR 0.13 (0.02–0.99)Protective 2
Vincristine metabolism
CYP3A4rs2740574A/G(*1B) Aplenc et al., 2003 [28] 8.628 (CCG grade 3–4)505 (CCG grade 0–2)Allelic OR 0 (0–0.75)Protective 2
Guilhaumou et al., 2011 [20]6.3Nr of neurotoxicity eventsChi–squarep = 1.00Not significant
Kishi et al., 2007 [13]AA: 79.6
AG + GG: 20.4		30 (grade 2–4)210 (grade 0–1)Dominant OR 1.37 (0.57–3.29)Not significant
GSTM1DeletionNon–null/nullKishi et al., 2007 [13]Non–null: 57.5
Null: 42.5		30 (grade 2–4)210 (grade 0–1)OR 0.46 (0.22–0.94)Protective2
VDR rs1544410G/AKishi et al., 2007 [13]GG: 45.8
AA and AG: 54.2		30 (grade 2–4)210 (grade 0–1)Recessive OR 2.22 (1.06–4.67)Risk
Cytoskeleton–associated
ACTG1rs1135989G/ACeppi et al., 2014 [8]36.538 (grade 3–4)214 (grade 0)Dominant OR 2.8 (1.3–6.3)Risk 1
CAPG rs2229668G/ACeppi et la. 2014 [8]12.639 (grade 3–4)214 (grade 0)Dominant OR 2.1 (1.1–3.7)Risk 1
rs3770102C/ACeppi et al., 2014 [8]41.439 (grade 3–4)214 (grade 0)Dominant OR 0.1 (0.01–0.8)Protective 1
CEP72 rs924607C/TDiouf et al., 2015—St. Jude cohort [9]36.764 (grade 2–4)158 (grade 0) Recessive OR 5.5 (2.5–12.2)Risk
Diouf et al., 2015—COG cohort [9]36.422 (grade 2–4)74 (grade 0) Recessive OR 3.8 (1.3–11.4)Risk
Gutierrez–Camino et al., 2016 [10]39.436 (WHO grade 2–4)106 (WHO grade 0–1)Recessive OR 0.7 (0.2–2.4)Not significant
Wright et al., 2019 [51] TT: 13.5
CT and CC: 86.5		156 (grade 2–4)56 (grade 0) Recessive OR 3.4 (0.9–12.6)Not significant
Zgheib et al., 2018 [50]36.923 (grade 2–4)107 (grade 0–1)Recessive OR 1.04 (0.32–3.43)Not significant
MAPT rs11867549A/GMartin–Guerrero et al., 2019 [49]22.518 (WHO grade 3–4)103 (WHO grade 0)Dominant OR 0.21 (0.04–0.96)Protective 2
SYNE2rs2781377G/AAbaji et al., 2018—QcALL cohort [52]7.835 (grade 3–4)201 (grade 0)Additive OR 2.5 (1.2–5.2)Risk
TUBB2B: miR–202rs12355840T/CMartin–Guerrero et al., 2019 [49]23.427 (WHO grade 1–2)103 (WHO grade 0)Dominant OR 2.88 (1.07–7.72)Risk
Hereditary neuropathy
SLC5A7rs1013940T/CWright et al., 2019 [51]15.2170 (grade 2–4)57 (grade 0)Additive OR 8.60 (1.68–44.15)Risk 3
Other (GWAS/EWAS studies)
BAHD1rs3803357C/AAbaji et al., 2018—QcALL cohort [52]41.735 (grade 3–4)201 (grade 0) Dominant OR 0.35 (0.2–0.7)Protective
COCH rs1045466T/GLi et al., 2020—POG cohort [53]38 Maximum neuropathy scoreDominant HR 0.27 (0.16–0.50)Protective
Li et al., 2020—ADVANCE cohort [53]33 Linear regression −3.56 (−5.45;−1.67)Protective
Chromosome 12/ chemerinrs7963521T/CLi et al., 2020—POG cohort [53]41Maximum neuropathy scoreAdditive HR 2.23 (1.49–3.35)Risk
Li et al., 2020—ADVANCE cohort [53]43 Additive HR 2.16 (0.53–3.70)Not significant
ETAA1 rs17032980A/GDiouf et al., 2015—St. Jude cohort [9]26.664 (grade 2–4)158 (grade 0)Allelic OR 3.17 (1.95–5.17)Risk
Diouf et al., 2015—COG cohort [9] 19.222 (grade 2–4)74 (grade 0)Allelic OR 10.4 (2.97–36.15)Risk
MRPL4rs10513762C/TAbaji et al., 2018—QcALL cohort [52]7.035 (grade 3–4)202 (grade 0)Dominant OR 3.3 (1.4–7.7)Risk
MTNR1Brs12786200C/TDiouf et al., 2015—St. Jude cohort [9] 22.764 (grade 2–4)158 (grade 0)Allelic OR 0.23 (0.13–0.40)Protective
Diouf et al., 2015—COG cohort [9] 20.722 (grade 2–4)74 (grade 0)Allelic OR 0.24 (0.08–0.76)Protective
Zgheib et al., 2018 [50]18.123 (grade 2–4)107 (grade 0–1)Dominant OR 0.59 (0.22–1.62)Not significant
NDUFAF6rs7818688C/ADiouf et al., 2015—St. Jude cohort [9]12.664 (grade 2–4)158 (grade 0)Allelic OR 4.26 (2.45–7.42)Risk
Diouf et al., 2015—COG cohort [9]14.122 (grade 2–4)74 (grade 0)Allelic OR 4.59 (1.35–15.59)Risk
TMEM215rs4463516C/GDiouf et al., 2015—St. Jude cohort [9]33.664 (grade 2–4)158 (grade 0)Allelic OR 3.17 (1.95–5.17)Risk
Diouf et al., 2015—COG cohort [9]24.222 (grade 2–4)74 (grade 0)Allelic OR 4.94 (1.65–14.79)Risk
miRNA
miR–4481rs7896283T/CGutierrez–Camino et al., 2017 [48]37.519 (WHO grade 3–4)128 (WHO grade 0)Dominant OR 4.69 (1.43–15.43)Risk 2
miR–6076rs35650931G/C Gutierrez–Camino et al., 2017 [48]8.747 (WHO grade 1–4)128 (WHO grade 0) Dominant OR 0.22 (0.05–0.97)Protective 2
SNP = single nucleotide polymorphism, MAF = minor allele frequency, CI = confidence interval, OR = odds ratio, ABCB1 = ATP binding cassette subfamily B member 1, ABCC1 = ATP binding cassette subfamily C member 1, ABCC2 = ATP binding cassette subfamily C member 2, RALPB1 = RalA binding protein 1, miR = microRNA, CYP3A4 = cytochrome P450 3A4, GSTM1 = glutathione S-transferase mu 1, VDR = vitamin D receptor, CAPG = capping actin protein gelsolin like, CEP72 = centrosomal protein 72, MAPT = microtubule associated protein tau, TUBB2B = tubulin beta 2B class IIB, ACTG1 = actin gamma 1, SYNE2 = spectrin repeat containing nuclear envelope protein 2, SLC5A7 = solute carrier family 5 member 7, BAHD1 = bromo adjacent homology domain containing 1, COCH = cochlin, ETAA1 = Ewing’s tumor-associated antigen 1, MRPL4 = mitochondrial ribosomal protein L4, MTNR1B = melatonin receptor 1B, NDUFAF6 = NADH: ubiquinone oxidoreductase complex assembly factor 6, TMEM215 = transmembrane protein 215. * Grades are referring to CTCAE grades unless mentioned otherwise. 1 Significance threshold not adjusted for multiple comparisons. 2 Significance threshold was not met after correcting for multiple comparisons. 3 Significance threshold was not adjusted for multiple comparisons, but associations p < 0.001 were prioritized. Odds ratios (OR) were defined as following: recessive OR meant that the risk of VIPN increased y-fold if two copies of the minor allele (genotype: aa) or genetic variation were present; dominant OR meant that the risk of VIPN increased y-fold if either one or two copies of the minor allele were present (genotypes: Aa or aa); allelic OR meant that the risk of VIPN increased y-fold with each additional copy of the minor allele or genetic variation; and the additive OR meant that the risk of VIPN increased y-fold for the heterozygous genotype (Aa) and 2y-fold for the homozygous variant genotype (aa).

3.4. Association between Pharmacogenomic Parameters and VIPN

Table 2 and Table 3 show an overview of all SNPs found to have a statistically significant and non-significant association with VIPN, respectively. Figure 2 shows a schematic overview of the function of genes associated with VIPN. Sixteen SNPs in three ATP-binding cassette transporter genes (ABCB1, ABCC1, ABCC2) and one SNP in an miRNA targeting ABCC1/RalA binding protein 1 (RALPB1) were described to be significantly associated with VIPN (Table 2). Ten SNPs were associated with a protective effect against VIPN, whereas seven SNPs were associated with an increased risk of VIPN. Of note, the strongest protective associations with high precision were reported for SNPs rs3740066 and rs12826 in ABCC2 (OR 0.23, 95% CI 0.10−0.53, and 0.24, 95% CI 0.10−0.54 respectively). The strongest risk association with acceptable precision was reported for rs3784867 in ABCC1 (OR 4.91, 95% CI 1.99−12.10).
In terms of metabolism-associated genes, a deletion in glutathione S-transferase mu 1 (GSTM1) and an SNP in vitamin D receptor (VDR) were implicated with a heightened and a decreased risk to VIPN, respectively (Table 2) [13]. Furthermore, six SNPs in cytoskeleton-associated genes or in miRNAs targeting those were associated with VIPN (microtubule associated protein tau (MAPT), targeting tubulin beta 2B class IIB (TUBB2), actin gamma 1 (ACTG1), capping actin protein gelsolin like (CAPG) and spectrin repeat containing nuclear envelope protein 2 (SYNE2)) (Table 2). Of those, two SNPs were related to microtubules (MAPT and TUBB2) and associated with a protective effect and an increased risk of VIPN, respectively (Table 2) [49]. The four other SNPs were located in cytoskeleton-associated genes (ACTG1, CAPG, and SYNE2) and associated with a CTCAE grade 3−4 VIPN (Table 2) [8,52]. The latter passed the stringent significance threshold for multiple comparisons, but the results could not be confirmed in a replication cohort [52]. The strongest protective association was noted for SNP rs3770102 in CAPG with an effect size of 0.1, although the uncertainty was high (95% CI 0.01−0.8). One SNP in a gene associated with hereditary neuropathies (solute carrier family 5 member 7 (SLC5A7)) resulted in an increased susceptibility to VIPN (Table 2) [51]. The reported effect size was large, but the size of the confidence interval indicated relatively high uncertainty (OR 8.60, 95% CI 1.68−44.15) Except for the SNP in SYNE2, all aforementioned SNPS were solely assessed in a discovery cohort and no replication studies were performed for any of those associations [52].
Four studies assessed the influence of SNP rs924607 in CEP72 on the development of CTCAE or WHO grade 2−4 VIPN [9,10,50,51]. In both their discovery and replication cohort, Diouf et al. described an increased risk of VIPN in patients with the risk genotype [9]. The strongest association was seen in the discovery cohort (OR 5.5, 95% CI 2.5−12.2); this effect size was smaller in the replication cohort (OR 3.8, 95% CI 1.3−11.4) where the prevalence of VIPN was also slightly lower (28.8 and 22.9% respectively). Three replication studies could not confirm these findings [10,50,51].
GWAS or EWAS demonstrated significant associations between VIPN and eight SNPs in genes previously not associated with neuropathy, vincristine mechanism of action or metabolism (Table 2). All studies first reporting these associations made use of both a discovery and replication cohort to validate their results [9,52,53]. SNPs in cochlin (COCH), Ewing’s tumor-associated antigen 1 (ETAA1), melatonin receptor 1B (MTNR1B), NADH: ubiquinone oxidoreductase complex assembly factor (NDUFAF6), and transmembrane protein 215 (TMEM215) were significantly associated with VIPN both in a discovery and replication cohort, whereas this relationship was only established in the discovery cohort for SNPs in bromo adjacent homology domain containing 1 (BAHD1), chromosome 12/chemerin, and mitochondrial ribosomal protein L4 (MRPL4). The described SNPS in BAHD1 and COCH were protective against VIPN. The strongest protective association with high precision was reported for the latter (OR 0.27, 95% CI 0.16−0.50). The SNPs in chromosome 12/chemerin, ETAA1, MRPL4, NDUFAF6, and TMEM215 were associated with an increased risk of VIPN. The SNP in ETAA1 showed a strong effect on risk of VIPN, especially in the replication cohort of Diouf et al., although the precision was relatively low (OR 10.4, 95% CI 2.97−36.15). Moreover, the SNPs in NDUFAF6 and TMEM215 also showed relatively large effect sizes with acceptable uncertainty in both a discovery and replication cohort. Finally, Diouf et al. described an SNP in melatonin receptor 1B (MTNR1B) as protective against VIPN both in a discovery and replication cohort with a large effect size and high precision (OR 0.23, 95% CI 0.13−0.40, and OR 0.24, 95% CI 0.08−0.76), but another study by Zgheib et al. could not confirm these results [9,50]. All significant associations passed the stringent threshold for multiple comparisons.
Gutierrez-Camino et al. found two SNPs in miRNA to be associated with VIPN, of which one miRNA could be related to the axon-guidance pathway, whereas the other could not be related to any known vincristine- or neurotoxicity-related pathway (Table 2) [48].
Several studies assessed the influence of covariates such as cumulative vincristine dosage, treatment protocol, and patient characteristics on their results, but these covariates did not have a significant influence on the reported associations (Table S2). Only the significant associations reported by Diouf et al. did not maintain their significance when corrected for genetically defined ancestry and cumulative vincristine dosage (Table S2) [9].
Table
Table 3. Single-nucleotide polymorphisms that were not significantly associated with vincristine-induced peripheral neuropathy in the pediatric oncology population.
Table 3. Single-nucleotide polymorphisms that were not significantly associated with vincristine-induced peripheral neuropathy in the pediatric oncology population.
GeneSNPAuthor and Year of Publication
ABCB1rs1045642Plasschaert et al., 2004 [22], Ceppi et al., 2014 [8], Zgheib et al., 2018 [50]
rs1128503Ceppi et al., 2014 [8], Zgheib et al., 2018 [50]
rs2032582Plasschaert et al., 2004 [22], Ceppi et al., 2014 [8]
ABCC2rs717620Zgheib et al., 2018 [50]
ACTG1rs1139405Ceppi et al., 2014 [8]
rs7406609Ceppi et al., 2014 [8]
CAPGrs6886Ceppi et al., 2014 [8]
CYP1A1rs4646903Abo-Bakr et al., 2017 1 [47]
GSTP1rs1695Kishi et al., 2007 [13], Abo-Bakr et al., 2017 1 [47]
GSTT1DeletionKishi et al., 2007 [13]
MAP4rs11268924Ceppi et al., 2014 [8]
rs1137524Ceppi et al., 2014 [8]
rs1875103Ceppi et al., 2014 [8]
rs11711953Ceppi et al., 2014 [8]
MDR1Exon 21, G > T/AKishi et al., 2007 [13]
Exon 26, C/TKishi et al., 2007 [13]
MTHFRrs1801133Kishi et al., 2007 [13]
rs1801131Kishi et al., 2007 [13]
SLC19A1rs1051266Kishi et al., 2007 [13]
TPMTCombined genotypes: 238GG, 460GG, 719AA/othersKishi et al., 2007 [13]
TUBBrs6070697Ceppi et al., 2014 [8]
rs10485828Ceppi et al., 2014 [8]
TYMSEnhancer repeat: others/3AND3Kishi et al., 2007 [13]
UGT1A1Enhancer repeat: others/7AND7Kishi et al., 2007 [13]
VDRrs2228570Kishi et al., 2007 [13]
XRCC1rs1799782Abo-Bakr et al., 2017 1 [47]
CYP1A1 = cytochrome P450 family 1 subfamily A member 1, GSTP1 = glutathione S-transferase pi 1, GSTT1 = glutathione S-transferase theta 1, MAP4 = microtubule-associated protein 4, MDR1 = multidrug resistance mutation 1, MTHFR = methylenetetrahydrofolate reductase, SLC19A1 = solute carrier family 19 member 1, TPMT = thiopurine methyltransferase, TYMS = thymidylate synthetase, UGT1A1 = uridine glucuronosyltransferase 1A1, XRCC1 = X-ray repair cross-complementing protein 1. 1 Association could not be tested due to small number of patients with VIPN.

3.5. CYP3A4 and CYP3A5

In regard to CYP3A4, Aplenc et al. found an SNP in CYP3A4 to be protective against VIPN [28], but two follow-up studies could not replicate these findings (Table 2) [13,20]. Furthermore, ten studies assessed the influence of CYP3A5 expression on the development of VIPN [8,13,14,16,20,28,29,30,31,32]. Of those studies, nine either presented a pre-calculated OR or raw data to calculate an OR and could thus be included in the meta-analysis [8,13,14,16,20,28,29,31,32]. If possible, dominant ORs were calculated based on data presented in the article or additional data provided by the authors (Tables S3 and S4). As shown in Figure 3, there was no statistically significant pooled effect between CYP3A5 expression status and the development of VIPN (pooled OR 0.69, 95% CI 0.38−1.26, I2 = 50%, τ2 = 0.33). The study by Kayilioğlu et al., which could not be included in the meta-analysis due to the unavailability of appropriate data, did not find a significant association either (Table S3). However, all studies either found no effect of CYP3A5 status or found that expression of CYP3A5 was a protective factor for VIPN; the opposite was not reported. Of note, the included studies all used different definitions for cases (patients with VIPN) and controls (patients without VIPN) (Table S3), likely contributing to the moderate heterogeneity. In addition, all studies performed genotyping to determine CYP3A5 expression status, whereas Renbarger et al., the study with the strongest association, used race as a surrogate. All studies compared expressers of CYP3A5 to non-expressers, except for Ceppi et al. who calculated allelic OR (Table S3). Evaluation of the funnel plot and Egger’s test for asymmetry were not indicative of obvious publication bias (p-value Egger’s test = 0.40) (Figure S1).

4. Discussion

This systematic review shows that pharmacogenomic parameters have a significant influence on VIPN in children with cancer and show potential for clinical relevance. Several SNPs in genes related to vincristine metabolism, hereditary neuropathy, the cytoskeleton and microtubules have been associated with VIPN. Furthermore, population-based GWAS and EWAS identified significant interactions with SNPs in genes previously unrelated to VIPN or vincristine. Our meta-analysis showed that CYP3A5 expression status was not a significant risk factor for VIPN.
Several significant associations were found between SNPs in the ABC family of genes (ABCB1, ABCC1, ABCC2). These genes code for transmembrane proteins that mediate vincristine efflux across cell membranes; variations may thus contribute to different vincristine levels and therefore VIPN (Figure 2) [54,55]. Three candidate gene studies described associations between ABCB1, ABCC1 and ABCC2 and VIPN, all in children with ALL [8,11,51]. Of note, the vast majority of the associations were reported in the same discovery cohort (Lopez-Lopez et al.) [11]. Interestingly, all associations between SNPs in ABCC2 and VIPN passed a significance threshold corrected for multiple comparisons and showed the largest effect sizes, suggesting that a stronger relationship may exist between ABCC2 and VIPN than between ABCB1/ABCC1 and VIPN. Indeed, several cell line studies have shown that ABCC2 function was associated with vincristine resistance or sensitivity [56,57,58]. However, none of the reported associations have been replicated in other cohorts and results should thus be interpreted with caution.
Furthermore, SNPs in cytoskeleton-associated genes were associated with VIPN (Figure 2). Vincristine exerts its cytostatic effect via binding to the β-subunit of tubulins, which inhibits microtubule polymerization and consequently causes arrest of mitosis in the metaphase [1,59]. During cell division, there is a well-known interaction between microtubules and the actin cytoskeleton; the latter contributes to mitotic spindle assembly and formation [60,61,62]. It is possible that SNPs in genes that affect microtubule formation or the actin cytoskeleton affect binding of vincristine to tubulins or the effect of vincristine binding to tubulins. While this can result in an altered risk of VIPN, one could also hypothesize that this influences the effect of vincristine on mitotic spindle disintegration and thus ultimately the cytotoxic effect. Should that be the case, patients with a lower risk of VIPN might also experience less antitumor effect in comparison with patients with a higher risk of VIPN, which would argue for dose individualization in which standard dose capping is not applied to every patient. Future studies assessing the relationship between VIPN incidence and long-term treatment outcome, correcting for received cumulative vincristine dosage, may provide further insight. Of note, the studies reporting these associations concerned predominantly white patients with ALL and except for one study, the reported associations have not been assessed in a replication cohort [8,49,52]. Therefore, these results regarding SNPs in microtubule- and cytoskeleton-associated genes should be interpreted with caution until independent replication is performed. An association that has been replicated in several independent studies is the association between rs924607 in CEP72 and VIPN. CEP72 encodes for a centrosomal protein that is required for adequate chromosome segregation [63,64]. Centrosomes enable correct alignment of chromosomes during mitosis by controlling the position and orientation of the microtubule spindles at the spindle poles [63,64]. A recent meta-analysis on the effect of this SNP in CEP72 on VIPN in children with ALL confirmed this finding across three studies in the continuation phase of treatment [65]. A clinical trial enrolling newly diagnosed children with ALL and lymphomas is randomizing patients with the high risk (TT) genotype between a decreased dosage (1.0 mg/m2) and conventional dosage (1.5 mg/m2) of vincristine during the continuation phase of treatment [66]. Recruitment is still ongoing. However, it is important to note that the effect of CEP72 on VIPN likely differs depending on treatment phase and genetic background, since this association was not significant in a white Spanish population in the induction phase and a white Arab population during induction and continuation phases [10,50].
In addition, we assessed the effect of CYP3A5 expression status on VIPN in a meta-analysis and found an overall pooled effect of 0.69 (95% CI 0.38−1.26) (Figure 2). Two studies reported a significant effect of CYP3A5 expression status on VIPN. Renbarger et al. found the strongest association, but it is important to note that they used race as a surrogate for CYP3A5 expression status [14]. This can be debated, since white children can express a CYP3A5 as well, albeit less often than in black children (10−20% and >55%, respectively) [67,68]. Furthermore, Kishi et al. found a significant association in a relatively large cohort (240 children) with a prevalence of 12.5% of severe VIPN [13]. However, all other studies could not replicate these findings, even those with sample sizes adequately powered to detect a difference, such as the population wide GWAS or EWAS. Nonetheless, no study reported CYP3A5 expression as a risk factor for VIPN. In conclusion, this meta-analysis shows that there is no significant effect of CYP3A5 expression status on VIPN.
The comparison and interpretation of the results of the included studies is limited due to heterogeneity in the study population, treatment protocol and varying assessment methods and definitions of VIPN. Firstly, cumulative vincristine dosage varied between 6.5 and 97.5 mg across studies. Since cumulative vincristine dosage likely is an independent risk factor for VIPN, it could be of influence when establishing a relationship between a pharmacogenomic parameter and VIPN [1,3]. However, except for Diouf et al., studies that included cumulative vincristine dosage or treatment protocol as covariates did not find an effect on outcome [8,9,11,13,49,51,52]. Another source of heterogeneity was the variety of measurement approaches to report VIPN. The majority of studies used the NCI-CTCAE for peripheral neuropathy, but other studies used the modified Balis scale, WHO scale, or other methods, to quantify VIPN. The sensitivity and specificity differ across assessment methods and their results can thus not be compared one-on-one [1,3,17]. Similarly, seven studies assessed VIPN prospectively, whereas the other studies assessed VIPN retrospectively. Retrospective VIPN assessment is less sensitive, especially when quantifying the presence of any grade or low grade VIPN [1,3,17]. Moreover, the majority of included studies had a study population of less than 150 patients and were thus limited by a relatively small sample size. This is further reinforced by the observation that most studies assessed the relationship between several genetic variations and VIPN and thus performed multiple comparisons. Therefore, it is advised to counteract the likelihood of false positives by adjusting the significance threshold with for example Bonferroni or False Discovery Rate (FDR) correction [69,70]. This systematic review shows that approximately half of the reported significant associations either did not pass a stringent significance threshold, or that it was not applied. However, associations that did not pass the stringent significance threshold could still be biologically relevant since (stringent) statistical significance should be interpreted within the context of the study design [35,39]. A lack of statistical significance could be an indication that the study was not adequately powered to detect a difference [35,39].
Independent replication is essential to validate the clinical significance of reported associations [33,34]. This could be facilitated by the uniform and reliable assessment of VIPN while employing a sensitive assessment method such as the ped-mTNS or ped-TNV [2,3]. A larger number of patients could be included if VIPN was consistently noted in patient charts. Subsequently, this would allow for reliable comparisons between studies. Furthermore, the growing availability of high-throughput techniques allows for genome- or exome wide analysis in an increasing number of studies. Interestingly, in this systematic review and meta-analysis, the majority of included studies followed a candidate gene approach, limiting the findings to pre-defined selection of genes. The studies that employed population-wide GWA or EWA analysis expanded our knowledge by uncovering genotype-phenotype associations that were not previously described in relation to VIPN. To actualize the potential of pharmacogenomic testing, future studies should apply sensitive and uniform measurement approaches to report VIPN while employing robust genotyping methods such as EWAS or GWAS. These data could be used by consortiums such as the Clinical Pharmacogenetics Implementation Consortium (CPIC) to create guidelines for clinicians to implement pharmacogenomic testing into clinical practice. The CPIC assesses pharmacogenomic studies by assigning levels of evidence and syntheses the results [71]. Another promising approach to assess the influence of pharmacogenomics on VIPN is combining the effect of different SNPs or genetic variations in one effect size, since single SNPs most likely have limited clinical relevance. Such an approach was adopted by Abaji et al., in which a combined-effect model was established to assess the additive effect of several SNPs associated with VIPN [52]. Patients were classified in risk groups according to their weighted genetic risk score; and this model could successfully predict the risk of VIPN in both their discovery and replication cohort [52].
The strength of this systematic review and meta-analysis is the inclusion of all studies assessing the relationship between pharmacogenomic parameters and VIPN. No restrictions were applied regarding patient or study characteristics. Furthermore, we followed the PRISMA guidelines and two independent authors performed screening, data extraction and risk of bias assessment. A standardized risk of bias assessment was performed with the validated EPHPP tool [45]. The weakness of this study is that the meta-analysis was limited to CYP3A5 expression status and that the other pharmacogenomic parameters could not be assessed in a meta-analysis.

5. Conclusions

From this systematic review, we can conclude that the following pharmacogenomic parameters have a significant influence on VIPN in children with cancer: SNPs in ABCB1, ABCC1, ABCC2, CYP3A4, GSTM1, VDR, ACTG1, CAPG, CEP72, MAPT, SYNE2, TUBB2B, SLC5A7, BAHD1, COCH, chromosome 12/chemerin, ETAA1, MRPL4, MTNR1B, NDUFAF6, TMEM215 and in three miRNAs. Our meta-analysis shows that CYP3A5 expression does not result in a heightened susceptibility of VIPN. To actualize the potential of pharmacogenomic testing, future research should prospectively assess VIPN with a sensitive measurement tool in both a discovery and replication cohort. Ultimately, the goal would be to develop an individualized protocol based on a patients’ genotype, taking all risk and protective genes into account, and subsequently give patients a dosage that limits the risk of VIPN while maintaining highest possible therapeutic efficacy. Dosage reductions or cessation of treatment, or for some patients even standardized dose capping, would no longer be necessary.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers14030612/s1, Supplementary Materials 1, Supplementary Materials 2, Figure S1: Funnel plot for evaluation of small study effects in the meta-analysis on the effect of CYP3A5 expression status on vincristine-induced peripheral neuropathy, Table S1: Risk of bias of the included studies, Table S2: Studies that assessed the effect of covariates on significant associations, Table S3: Effect of CYP3A5 single-nucleotide polymorphisms (SNPs) on vincristine-induced peripheral neuropathy (VIPN). Reference allele is *1 (extensive metabolizer), variant alleles are *3 (rs776746), *6 (rs10264272), and *7 (rs41303343). Patients with at least one *1 allele are considered to be expressers of CYP3A5. Patients without *1 allele are considered to be non-expressers of CYP3A5. Recessive OR: non-expresser (*3/*3) compared to expresser (other genotypes), Table S4: Raw data used to calculate odds ratios (OR) for the effect of CYP3A5 single-nucleotide polymorphisms (SNPs) on vincristine-induced peripheral neuropathy (VIPN). Reference allele is *1 (extensive metabolizer), variant alleles are *3 (rs776746), *6 (rs10264272), and *7 (rs41303343). Patients with at least one *1 allele are considered to be expressers of CYP3A5. Patients without *1 allele are considered to be non-expressers of CYP3A5. Dominant OR: non-expresser (*3/*3) compared to expresser (other genotypes).

Author Contributions

Conceptualization: A.U., M.E.v.d.V., G.J.L.K. Methodology: A.U., C.L.G.N., J.C.F.K. Validation: A.U., C.L.G.N. Formal analysis: A.U., C.L.G.N. Writing—original draft preparation: A.U. Writing—review & editing: A.U., M.E.v.d.V., F.N., A.D.R.H., G.J.L.K. Visualization: A.U. Supervision: M.E.v.d.V., F.N., A.D.R.H., G.J.L.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available in this article and in the Supplementary Materials.

Acknowledgments

The authors would like to thank Skiles et al. and Guilhaumou et al. for contributing original data to this study. Furthermore, they would like to gracefully acknowledge Sharon Remmelzwaal for advice on the statistical analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mora, E.; Smith, E.M.L.; Donohoe, C.; Hertz, D.L. Vincristine-induced peripheral neuropathy in pediatric cancer patients. Am. J. Cancer Res. 2016, 6, 2416–2430. [Google Scholar] [PubMed]
  2. Smith, E.M.L.; Kuisell, C.; Cho, Y.; Kanzawa-Lee, G.A.; Gilchrist, L.S.; Park, S.B.; Scott, M.R.; Alberti, P. Characteristics and patterns of pediatric chemotherapy-induced peripheral neuropathy: A systematic review. Cancer Treat. Res. Commun. 2021, 28, 100420. [Google Scholar] [CrossRef] [PubMed]
  3. Van de Velde, M.E.; Kaspers, G.L.; Abbink, F.C.H.; Wilhelm, A.J.; Ket, J.C.F.; van den Berg, M.H. Vincristine-induced peripheral neuropathy in children with cancer: A systematic review. Crit. Rev. Oncol. Hematol. 2017, 114, 114–130. [Google Scholar] [CrossRef] [PubMed]
  4. Madsen, M.L.; Due, H.; Ejskjær, N.; Jensen, P.; Madsen, J.; Dybkær, K. Aspects of vincristine-induced neuropathy in hematologic malignancies: A systematic review. Cancer Chemother. Pharmacol. 2019, 84, 471–485. [Google Scholar] [CrossRef] [Green Version]
  5. Lavoie Smith, E.M.; Li, L.; Chiang, C.; Thomas, K.; Hutchinson, R.J.; Wells, E.M.; Ho, R.H.; Skiles, J.; Chakraborty, A.; Bridges, C.M.; et al. Patterns and severity of vincristine-induced peripheral neuropathy in children with acute lymphoblastic leukemia. J. Peripher. Nerv. Syst. JPNS 2015, 20, 37–46. [Google Scholar] [CrossRef]
  6. Van de Velde, M.E.; van den Berg, M.H.; Kaspers, G.J.L.; Abbink, F.C.H.; Twisk, J.W.R.; van der Sluis, I.M.; van den Bos, C.; van den Heuvel-Eibrink, M.M.; Segers, H.; Chantrain, C.; et al. The association between vincristine-induced peripheral neuropathy and health-related quality of life in children with cancer. Cancer Med. 2021, 10, 8172–8181. [Google Scholar] [CrossRef]
  7. Tay, C.G.; Lee, V.W.M.; Ong, L.C.; Goh, K.J.; Ariffin, H.; Fong, C.Y. Vincristine-induced peripheral neuropathy in survivors of childhood acute lymphoblastic leukaemia. Pediatr. Blood Cancer 2017, 64, e26471. [Google Scholar] [CrossRef]
  8. Ceppi, F.; Langlois-Pelletier, C.; Gagné, V.; Rousseau, J.; Ciolino, C.; De Lorenzo, S.; Kevin, K.M.; Cijov, D.; Sallan, S.E.; Silverman, L.B.; et al. Polymorphisms of the vincristine pathway and response to treatment in children with childhood acute lymphoblastic leukemia. Pharmacogenomics 2014, 15, 1105–1116. [Google Scholar] [CrossRef] [Green Version]
  9. Diouf, B.; Crews, K.R.; Lew, G.; Pei, D.; Cheng, C.; Bao, J.; Zheng, J.J.; Yang, W.; Fan, Y.; Wheeler, H.E.; et al. Association of an inherited genetic variant with vincristine-related peripheral neuropathy in children with acute lymphoblastic leukemia. JAMA 2015, 313, 815–823. [Google Scholar] [CrossRef]
  10. Gutierrez-Camino, A.; Martin-Guerrero, I.; Lopez-Lopez, E.; Echebarria-Barona, A.; Zabalza, I.; Ruiz, I.; Guerra-Merino, I.; Garcia-Orad, A. Lack of association of the CEP72 rs924607 TT genotype with vincristine-related peripheral neuropathy during the early phase of pediatric acute lymphoblastic leukemia treatment in a Spanish population. Pharm. Genom. 2016, 26, 100–102. [Google Scholar] [CrossRef]
  11. Lopez-Lopez, E.; Gutierrez-Camino, A.; Astigarraga, I.; Navajas, A.; Echebarria-Barona, A.; Garcia-Miguel, P.; Garcia de Andoin, N.; Lobo, C.; Guerra-Merino, I.; Martin-Guerrero, I.; et al. Vincristine pharmacokinetics pathway and neurotoxicity during early phases of treatment in pediatric acute lymphoblastic leukemia. Pharmacogenomics 2016, 17, 731–741. [Google Scholar] [CrossRef]
  12. Anghelescu, D.L.; Faughnan, L.G.; Jeha, S.; Relling, M.V.; Hinds, P.S.; Sandlund, J.T.; Cheng, C.; Pei, D.; Hankins, G.; Pauley, J.L.; et al. Neuropathic pain during treatment for childhood acute lymphoblastic leukemia. Pediatr. Blood Cancer 2011, 57, 1147–1153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Kishi, S.; Cheng, C.; French, D.; Pei, D.; Das, S.; Cook, E.H.; Hijiya, N.; Rizzari, C.; Rosner, G.L.; Frudakis, T.; et al. Ancestry and pharmacogenetics of antileukemic drug toxicity. Blood 2007, 109, 4151–4157. [Google Scholar] [CrossRef] [PubMed]
  14. Renbarger, J.L.; McCammack, K.C.; Rouse, C.E.; Hall, S.D. Effect of race on vincristine-associated neurotoxicity in pediatric acute lymphoblastic leukemia patients. Pediatr. Blood Cancer 2008, 50, 769–771. [Google Scholar] [CrossRef]
  15. Smitherman, A.B.; Faircloth, C.B.; Deal, A.; Troy, M.; Gold, S.H. Vincristine toxicity with co-administration of fluconazole during induction therapy for pediatric acute lymphoblastic leukemia. Pediatr. Blood Cancer 2017, 64, e26525. [Google Scholar] [CrossRef]
  16. Skiles, J.L.; Chiang, C.; Li, C.H.; Martin, S.; Smith, E.L.; Olbara, G.; Jones, D.R.; Vik, T.A.; Mostert, S.; Abbink, F.; et al. CYP3A5 genotype and its impact on vincristine pharmacokinetics and development of neuropathy in Kenyan children with cancer. Pediatr. Blood Cancer 2018, 65, e26854. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Lavoie Smith, E.M.; Li, L.; Hutchinson, R.J.; Ho, R.; Burnette, W.B.; Wells, E.; Bridges, C.; Renbarger, J. Measuring vincristine-induced peripheral neuropathy in children with acute lymphoblastic leukemia. Cancer Nurs. 2013, 36, E49–E60. [Google Scholar] [CrossRef] [Green Version]
  18. Van de Velde, M.E.; Panetta, J.C.; Wilhelm, A.J.; van den Berg, M.H.; van der Sluis, I.M.; van den Bos, C.; Abbink, F.C.H.; van den Heuvel-Eibrink, M.M.; Segers, H.; Chantrain, C.; et al. Population pharmacokinetics of vincristine related to infusion duration and peripheral neuropathy in pediatric oncology patients. Cancers 2020, 12, 1789. [Google Scholar] [CrossRef]
  19. Crom, W.R.; de Graaf, S.S.; Synold, T.; Uges, D.R.; Bloemhof, H.; Rivera, G.; Christensen, M.L.; Mahmoud, H.; Evans, W.E. Pharmacokinetics of vincristine in children and adolescents with acute lymphocytic leukemia. J. Pediatr. 1994, 125, 642–649. [Google Scholar] [CrossRef]
  20. Guilhaumou, R.; Simon, N.; Quaranta, S.; Verschuur, A.; Lacarelle, B.; Andre, N.; Solas, C. Population pharmacokinetics and pharmacogenetics of vincristine in paediatric patients treated for solid tumour diseases. Cancer Chemother. Pharmacol. 2011, 68, 1191–1198. [Google Scholar] [CrossRef]
  21. Moore, A.S.; Norris, R.; Price, G.; Nguyen, T.; Ni, M.; George, R.; van Breda, K.; Duley, J.; Charles, B.; Pinkerton, R. Vincristine pharmacodynamics and pharmacogenetics in children with cancer: A limited-sampling, population modelling approach. J. Paediatr. Child Health 2011, 47, 875–882. [Google Scholar] [CrossRef] [PubMed]
  22. Plasschaert, S.L.; Groninger, E.; Boezen, M.; Kema, I.; de Vries, E.G.; Uges, D.; Veerman, A.J.; Kamps, W.A.; Vellenga, E.; de Graaf, S.S.; et al. Influence of functional polymorphisms of the MDR1 gene on vincristine pharmacokinetics in childhood acute lymphoblastic leukemia. Clin. Pharm. 2004, 76, 220–229. [Google Scholar]
  23. Dennison, J.B.; Jones, D.R.; Renbarger, J.L.; Hall, S.D. Effect of CYP3A5 expression on vincristine metabolism with human liver microsomes. J. Pharmacol. Exp. Ther. 2007, 321, 553–563. [Google Scholar] [CrossRef] [PubMed]
  24. Arbitrio, M.; Scionti, F.; Di Martino, M.T.; Pensabene, L.; Tassone, P.; Tagliaferri, P. Pharmacogenetics/pharmacogenomics of drug-metabolizing enzymes and transporters. In Reference Module in Biomedical Sciences; Elsevier: Amsterdam, The Netherlands, 2021. [Google Scholar]
  25. Kuehl, P.; Zhang, J.; Lin, Y.; Lamba, J.; Assem, M.; Schuetz, J.; Watkins, P.B.; Daly, A.; Wrighton, S.A.; Hall, S.D.; et al. Sequence diversity in CYP3A promoters and characterization of the genetic basis of polymorphic CYP3A5 expression. Nat. Genet. 2001, 27, 383–391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Blanco, J.G.; Edick, M.J.; Hancock, M.L.; Winick, N.J.; Dervieux, T.; Amylon, M.D.; Bash, R.O.; Behm, F.G.; Camitta, B.M.; Pui, C.-H.; et al. Genetic polymorphisms in CYP3A5, CYP3A4 and NQO1 in children who developed therapy-related myeloid malignancies. Pharm. Genom. 2002, 12, 605–611. [Google Scholar] [CrossRef]
  27. Roy, J.-N.; Lajoie, J.; Zijenah, L.S.; Barama, A.; Poirier, C.; Ward, B.J.; Roger, M. CYP3A5 genetic polymorphisms in different ethnic populations. Drug Metab. Dispos. 2005, 33, 884–887. [Google Scholar] [CrossRef]
  28. Aplenc, R.; Glatfelter, W.; Han, P.; Rappaport, E.; La, M.; Cnaan, A.; Blackwood, M.A.; Lange, B.; Rebbeck, T. CYP3A genotypes and treatment response in paediatric acute lymphoblastic leukaemia. Br. J. Haematol. 2003, 122, 240–244. [Google Scholar] [CrossRef]
  29. Egbelakin, A.; Ferguson, M.J.; MacGill, E.A.; Lehmann, A.S.; Topletz, A.R.; Quinney, S.K.; Li, L.; McCammack, K.C.; Hall, S.D.; Renbarger, J.L. Increased risk of vincristine neurotoxicity associated with low CYP3A5 expression genotype in children with acute lymphoblastic leukemia. Pediatr. Blood Cancer 2011, 56, 361–367. [Google Scholar] [CrossRef] [Green Version]
  30. Kayilioğlu, H.; Kocak, U.; Kan Karaer, D.; Percin, E.F.; Sal, E.; Tekkesin, F.; Isik, M.; Oner, N.; Belen, F.B.; Yilmaz Keskin, E.; et al. Association of CYP3A5 expression and vincristine neurotoxicity in pediatric malignancies in Turkish population. J. Pediatr. Hematol. Oncol. 2017, 39, 458–462. [Google Scholar] [CrossRef]
  31. McClain, C.A.; Bernhardt, M.B.; Berger, A.; Bernini, J.C.; Marquez-Do, D.; Winslow, R.; Scheurer, M.E.; Schafer, E.S. Pharmacogenetic association with neurotoxicity in Hispanic children with acute lymphoblastic leukaemia. Br. J. Haematol. 2018, 181, 684–687. [Google Scholar] [CrossRef] [Green Version]
  32. Sims, R.P. The effect of race on the CYP3A-mediated metabolism of vincristine in pediatric patients with acute lymphoblastic leukemia. J. Oncol. Pharm. Pract. 2016, 22, 76–81. [Google Scholar] [CrossRef] [PubMed]
  33. Lee, J.W.; Aminkeng, F.; Bhavsar, A.P.; Shaw, K.; Carleton, B.C.; Hayden, M.R.; Ross, C.J. The emerging era of pharmacogenomics: Current successes, future potential, and challenges. Clin. Genet. 2014, 86, 21–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Weinshilboum, R.M.; Wang, L. Pharmacogenomics: Precision medicine and drug response. Mayo Clin. Proc. 2017, 92, 1711–1722. [Google Scholar] [CrossRef]
  35. Alwi, Z.B. The Use of SNPs in pharmacogenomics studies. Malays. J. Med. Sci. 2005, 12, 4–12. [Google Scholar]
  36. Wang, J.; Pang, G.S.; Chong, S.S.; Lee, C.G. SNP web resources and their potential applications in personalized medicine. Curr. Drug Metab. 2012, 13, 978–990. [Google Scholar] [CrossRef] [PubMed]
  37. Auton, A.; Abecasis, G.R.; Altshuler, D.M.; Durbin, R.M.; Abecasis, G.R.; Bentley, D.R.; Chakravarti, A.; Clark, A.G.; Donnelly, P.; Eichler, E.E.; et al. A global reference for human genetic variation. Nature 2015, 526, 68–74. [Google Scholar] [CrossRef] [Green Version]
  38. Karczewski, K.J.; Francioli, L.C.; Tiao, G.; Cummings, B.B.; Alföldi, J.; Wang, Q.; Collins, R.L.; Laricchia, K.M.; Ganna, A.; Birnbaum, D.P.; et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 2020, 581, 434–443. [Google Scholar] [CrossRef]
  39. Clarke, G.M.; Anderson, C.A.; Pettersson, F.H.; Cardon, L.R.; Morris, A.P.; Zondervan, K.T. Basic statistical analysis in genetic case-control studies. Nat. Protoc. 2011, 6, 121–133. [Google Scholar] [CrossRef] [Green Version]
  40. Zondervan, K.T.; Cardon, L.R. Designing candidate gene and genome-wide case-control association studies. Nat. Protoc. 2007, 2, 2492–2501. [Google Scholar] [CrossRef] [Green Version]
  41. Lennard, L.; Cartwright, C.S.; Wade, R.; Vora, A. Thiopurine methyltransferase and treatment outcome in the UK acute lymphoblastic leukaemia trial ALL2003. Br. J. Haematol. 2015, 170, 550–558. [Google Scholar] [CrossRef] [Green Version]
  42. Stocco, G.; Cheok, M.H.; Crews, K.R.; Dervieux, T.; French, D.; Pei, D.; Yang, W.; Cheng, C.; Pui, C.H.; Relling, M.V.; et al. Genetic polymorphism of inosine triphosphate pyrophosphatase is a determinant of mercaptopurine metabolism and toxicity during treatment for acute lymphoblastic leukemia. Clin. Pharm. 2009, 85, 164–172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Liberati, A.; Altman, D.G.; Tetzlaff, J.; Mulrow, C.; Gøtzsche, P.C.; Ioannidis, J.P.A.; Clarke, M.; Devereaux, P.J.; Kleijnen, J.; Moher, D. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate healthcare interventions: Explanation and elaboration. BMJ 2009, 339, b2700. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Uittenboogaard, A.; Van de Velde, M.E.; Kaspers, G.J. Pharmacogenomics of Vincristine-Induced Peripheral Neuropathy: A Systematic Review (CRD42021210437). Available online: https://www.crd.york.ac.uk/prospero/display_record.php?RecordID=210437 (accessed on 30 September 2021).
  45. Armijo-Olivo, S.; Stiles, C.R.; Hagen, N.A.; Biondo, P.D.; Cummings, G.G. Assessment of study quality for systematic reviews: A comparison of the cochrane collaboration risk of bias tool and the effective public health practice project quality assessment tool: Methodological research. J. Eval. Clin. Pract. 2012, 18, 12–18. [Google Scholar] [CrossRef] [PubMed]
  46. Harrer, M.; Cuijpers, P.; Furukawa, T.A.; Ebert, D.D. Doing Meta-Analysis with R: A Hands-On Guide, 1st ed.; Chapman & Hall/CRC Press: Boca Raton, FL, USA; London, UK, 2021. [Google Scholar]
  47. Abo-Bakr, A.; Mossallam, G.; El Azhary, N.; Hafez, H.; Badawy, R. Impact of CYP1A1, GSTP1 and XRCC1 genes polymorphisms on toxicity and response to chemotherapy in childhood acute lymphoblastic leukemia. J. Egypt. Natl. Cancer Inst. 2017, 29, 127–133. [Google Scholar] [CrossRef] [PubMed]
  48. Gutierrez-Camino, Ã.; Umerez, M.; Martin-Guerrero, I.; García de Andoin, N.; Santos, B.; Sastre, A.; Echebarria-Barona, A.; Astigarraga, I.; Navajas, A.; Garcia-Orad, A. Mir-pharmacogenetics of Vincristine and peripheral neurotoxicity in childhood B-cell acute lymphoblastic leukemia. Pharm. J. 2017, 18, 704–712. [Google Scholar] [CrossRef] [PubMed]
  49. Martin-Guerrero, I.; Gutierrez-Camino, A.; Echebarria-Barona, A.; Astigarraga, I.; Garcia de Andoin, N.; Navajas, A.; Garcia-Orad, A. Variants in vincristine pharmacodynamic genes involved in neurotoxicity at induction phase in the therapy of pediatric acute lymphoblastic leukemia. Pharm. J. 2019, 19, 564–569. [Google Scholar] [CrossRef]
  50. Zgheib, N.K.; Ghanem, K.M.; Tamim, H.; Aridi, C.; Shahine, R.; Tarek, N.; Saab, R.; Abboud, M.R.; El-Solh, H.; Muwakkit, S.A. Genetic polymorphisms in candidate genes are not associated with increased vincristine-related peripheral neuropathy in Arab children treated for acute childhood leukemia: A single institution study. Pharm. Genom. 2018, 28, 189–195. [Google Scholar] [CrossRef]
  51. Wright, G.E.B.; Amstutz, U.; Drögemöller, B.I.; Shih, J.; Rassekh, S.R.; Hayden, M.R.; Carleton, B.C.; Ross, C.J.D. Pharmacogenomics of vincristine-induced peripheral neuropathy implicates pharmacokinetic and inherited neuropathy genes. Clin. Pharmacol. Ther. 2019, 105, 402–410. [Google Scholar] [CrossRef] [Green Version]
  52. Abaji, R.; Ceppi, F.; Patel, S.; Gagné, V.; Xu, C.J.; Spinella, J.F.; Colombini, A.; Parasole, R.; Buldini, B.; Basso, G.; et al. Genetic risk factors for VIPN in childhood acute lymphoblastic leukemia patients identified using whole-exome sequencing. Pharmacogenomics 2018, 19, 1181–1193. [Google Scholar] [CrossRef]
  53. Li, L.; Sajdyk, T.; Smith, E.M.L.; Chang, C.W.; Li, C.; Ho, R.H.; Hutchinson, R.; Wells, E.; Skiles, J.L.; Winick, N.; et al. Genetic variants associated with vincristine-induced peripheral neuropathy in two populations of children with acute lymphoblastic leukemia. Clin. Pharm. 2019, 105, 1421–1428. [Google Scholar] [CrossRef]
  54. Dean, M.; Rzhetsky, A.; Allikmets, R. The human ATP-binding cassette (ABC) transporter superfamily. Genome Res. 2001, 11, 1156–1166. [Google Scholar] [CrossRef] [PubMed]
  55. Hodges, L.M.; Markova, S.M.; Chinn, L.W.; Gow, J.M.; Kroetz, D.L.; Klein, T.E.; Altman, R.B. Very important pharmacogene summary: ABCB1 (MDR1, P-glycoprotein). Pharm. Genom. 2011, 21, 152–161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Cui, Y.; König, J.; Buchholz, J.K.; Spring, H.; Leier, I.; Keppler, D. Drug resistance and ATP-dependent conjugate transport mediated by the apical multidrug resistance protein, MRP2, permanently expressed in human and canine cells. Mol. Pharmacol. 1999, 55, 929–937. [Google Scholar] [PubMed]
  57. Folmer, Y.; Schneider, M.; Blum, H.E.; Hafkemeyer, P. Reversal of drug resistance of hepatocellular carcinoma cells by adenoviral delivery of anti-ABCC2 antisense constructs. Cancer Gene Ther. 2007, 14, 875–884. [Google Scholar] [CrossRef]
  58. Liu, Z.; Qiu, M.; Tang, Q.-L.; Liu, M.; Lang, N.; Bi, F. Establishment and biological characteristics of oxaliplatin-resistant human colon cancer cell lines. Chin. J. Cancer 2010, 29, 661–667. [Google Scholar] [CrossRef] [Green Version]
  59. Below, J.; Das, M.J. Vincristine. In StatPearls; StatPearls Publishing LLC.: Treasure Island, FL, USA, 2021. [Google Scholar]
  60. Lancaster, O.M.; Baum, B. Shaping up to divide: Coordinating actin and microtubule cytoskeletal remodelling during mitosis. Semin. Cell Dev. Biol. 2014, 34, 109–115. [Google Scholar] [CrossRef]
  61. Kunda, P.; Baum, B. The actin cytoskeleton in spindle assembly and positioning. Trends Cell Biol. 2009, 19, 174–179. [Google Scholar] [CrossRef]
  62. Wang, Y.; Stear, J.H.; Swain, A.; Xu, X.; Bryce, N.S.; Carnell, M.; Alieva, I.B.; Dugina, V.B.; Cripe, T.P.; Stehn, J.; et al. Drug targeting the actin cytoskeleton potentiates the cytotoxicity of low dose vincristine by abrogating actin-mediated repair of spindle defects. Mol. Cancer Res. 2020, 18, 1074–1087. [Google Scholar] [CrossRef] [Green Version]
  63. Meraldi, P. Centrosomes in spindle organization and chromosome segregation: A mechanistic view. Chromosome Res. 2016, 24, 19–34. [Google Scholar] [CrossRef]
  64. Oshimori, N.; Li, X.; Ohsugi, M.; Yamamoto, T. Cep72 regulates the localization of key centrosomal proteins and proper bipolar spindle formation. EMBO J. 2009, 28, 2066–2076. [Google Scholar] [CrossRef] [Green Version]
  65. Zečkanović, A.; Jazbec, J.; Kavčič, M. Centrosomal protein72 rs924607 and vincristine-induced neuropathy in pediatric acute lymphocytic leukemia: Meta-analysis. Future Sci. OA 2020, 6, FSO582. [Google Scholar] [CrossRef]
  66. Total Therapy XVII for Newly Diagnosed Patients with Acute Lymphoblastic Leukemia and Lymphoma. ClinicalTrials.gov Identifier: NCT03117751. Available online: https://clinicaltrials.gov/ct2/show/NCT03117751 (accessed on 30 September 2021).
  67. Garsa, A.A.; McLeod, H.L.; Marsh, S. CYP3A4 and CYP3A5genotyping by Pyrosequencing. BMC Med Genet. 2005, 6, 19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Lamba, J.K.; Lin, Y.S.; Schuetz, E.G.; Thummel, K.E. Genetic contribution to variable human CYP3A-mediated metabolism. Adv. Drug Deliv. Rev. 2002, 54, 1271–1294. [Google Scholar] [CrossRef]
  69. Chen, S.-Y.; Feng, Z.; Yi, X. A general introduction to adjustment for multiple comparisons. J. Thorac. Dis 2017, 9, 1725–1729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Menyhart, O.; Weltz, B.; Győrffy, B. MultipleTesting.com: A tool for life science researchers for multiple hypothesis testing correction. PLoS ONE 2021, 16, e0245824. [Google Scholar] [CrossRef]
  71. Relling, M.V.; Klein, T.E. CPIC: Clinical pharmacogenetics implementation consortium of the pharmacogenomics research network. Clin. Pharmacol. Ther. 2011, 89, 464–467. [Google Scholar] [CrossRef]
Cancers 14 00612 g001 550
Figure 1. PRISMA flow diagram identification of studies included in the systematic review. VIPN = vincristine-induced peripheral neuropathy.
Figure 1. PRISMA flow diagram identification of studies included in the systematic review. VIPN = vincristine-induced peripheral neuropathy.
Cancers 14 00612 g001
Cancers 14 00612 g002 550
Figure 2. Schematic overview of the function of genes associated with VIPN. Red: described SNPs in this gene are associated with a higher risk of VIPN; green: described SNPs in this gene are associated with a lower risk of VIPN, brown: described SNPs in this gene are associated with both a higher and lower risk of VIPN (different per SNP). Created with BioRender.com.
Figure 2. Schematic overview of the function of genes associated with VIPN. Red: described SNPs in this gene are associated with a higher risk of VIPN; green: described SNPs in this gene are associated with a lower risk of VIPN, brown: described SNPs in this gene are associated with both a higher and lower risk of VIPN (different per SNP). Created with BioRender.com.
Cancers 14 00612 g002
Cancers 14 00612 g003 550
Figure 3. Forest plot showing the effect of CYP3A5 expression status on VIPN, ORs describe the effect of expression of CYP3A5 in comparison to non-expression of CYP3A5 [8,13,14,16,20,28,29,31,32]. The functional allele is *1 and variant alleles are *3 (rs776746), *6 (rs10264272), and *7 (rs41303343). A dominant model was adopted: patients with at least one *1 allele are considered to be expressers of CYP3A5. Patients without *1 allele are considered to be non-expressers of CYP3A5.
Figure 3. Forest plot showing the effect of CYP3A5 expression status on VIPN, ORs describe the effect of expression of CYP3A5 in comparison to non-expression of CYP3A5 [8,13,14,16,20,28,29,31,32]. The functional allele is *1 and variant alleles are *3 (rs776746), *6 (rs10264272), and *7 (rs41303343). A dominant model was adopted: patients with at least one *1 allele are considered to be expressers of CYP3A5. Patients without *1 allele are considered to be non-expressers of CYP3A5.
Cancers 14 00612 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Uittenboogaard, A.; Neutel, C.L.G.; Ket, J.C.F.; Njuguna, F.; Huitema, A.D.R.; Kaspers, G.J.L.; van de Velde, M.E. Pharmacogenomics of Vincristine-Induced Peripheral Neuropathy in Children with Cancer: A Systematic Review and Meta-Analysis. Cancers 2022, 14, 612. https://doi.org/10.3390/cancers14030612

AMA Style

Uittenboogaard A, Neutel CLG, Ket JCF, Njuguna F, Huitema ADR, Kaspers GJL, van de Velde ME. Pharmacogenomics of Vincristine-Induced Peripheral Neuropathy in Children with Cancer: A Systematic Review and Meta-Analysis. Cancers. 2022; 14(3):612. https://doi.org/10.3390/cancers14030612

Chicago/Turabian Style

Uittenboogaard, Aniek, Céline L. G. Neutel, Johannes C. F. Ket, Festus Njuguna, Alwin D. R. Huitema, Gertjan J. L. Kaspers, and Mirjam E. van de Velde. 2022. "Pharmacogenomics of Vincristine-Induced Peripheral Neuropathy in Children with Cancer: A Systematic Review and Meta-Analysis" Cancers 14, no. 3: 612. https://doi.org/10.3390/cancers14030612

APA Style

Uittenboogaard, A., Neutel, C. L. G., Ket, J. C. F., Njuguna, F., Huitema, A. D. R., Kaspers, G. J. L., & van de Velde, M. E. (2022). Pharmacogenomics of Vincristine-Induced Peripheral Neuropathy in Children with Cancer: A Systematic Review and Meta-Analysis. Cancers, 14(3), 612. https://doi.org/10.3390/cancers14030612

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop