- 1Department of Pharmacy, The First Affiliated Hospital of Third Military Medical University (Army Medical University), Chongqing, China
- 2Department of Pharmacy, Chongqing Emergency Medical Center, Chongqing, China
Background: ABCG2 and ABCB1 are genes related to the pharmacokinetics of sunitinib and have been associated with its toxicity and efficacy. However, the results have been controversial. This study aimed to evaluate the associations of ABCG2 and ABCB1 polymorphisms with sunitinib-induced toxicity and efficacy in renal cell carcinoma (RCC) by meta-analysis.
Methods:PubMed, EMBASE, Cochrane Library, and Web of Science were systematically searched for studies investigating the associations of the ABCG2 rs2231142 polymorphism with sunitinib-induced toxicity and the associations of the ABCB1 rs1128503 and ABCB1 rs2032582 polymorphisms with sunitinib-induced toxicity and clinical outcomes. The associations were evaluated by effect size (ES) with 95% confidence intervals (CIs).
Results: Eight and five studies were included in the toxicity and efficacy analysis, respectively, including a total of 1081 RCC patients. The ABCG2 rs2231142 A allele was associated with an increased risk of sunitinib-induced thrombocytopenia and hand-foot syndrome (HFS) in Asians (ES = 1.65, 95% CI = 1.15–2.36, p = 0.006; ES = 1.52, 95% CI = 1.02–2.27, p = 0.041). However, the ABCG2 rs2231142 polymorphism was not associated with sunitinib-induced hypertension or neutropenia (ES = 1.09, 95% CI = 0.69–1.73, p = 0.701; ES = 0.87, 95% CI = 0.57–1.31, p = 0.501). Compared with the C allele, the ABCB1 rs1128503 T allele was associated with a decreased risk of sunitinib-induced hypertension but worse progression-free survival (PFS) (ES = 0.44, 95% CI = 0.26–0.77, p = 0.004; ES = 1.36, 95% CI = 1.07–1.73, p = 0.011). There was no significant association between the T allele or C allele of ABCB1 rs1128503 and overall survival (OS) (ES = 0.82, 95% CI = 0.61–1.10, p = 0.184). The ABCB1 rs2032582 T allele was associated with worse PFS than the other alleles (ES = 1.46, 95% CI = 1.14–1.87, p = 0.003), while there was no significant association between the T allele or other alleles and sunitinib-induced hypertension, HFS, or OS (ES = 0.77, 95% CI = 0.46–1.29, p = 0.326; ES = 1.02, 95% CI = 0.65–1.62, p = 0.919; ES = 1.32, 95% CI = 0.85–2.05, p = 0.215).
Conclusion: The results indicate that the ABCG2 rs2231142 polymorphism may serve as a predictor of sunitinib-induced thrombocytopenia and HFS in Asians, while the ABCB1 rs1128503 polymorphism may serve as a predictor of sunitinib-induced hypertension, and both the ABCB1 rs1128503 and rs2032582 polymorphisms may serve as predictors of PFS in RCC. These results suggest a possible application of individualized use of sunitinib according to the genetic background of patients.
Introduction
Renal cell carcinoma (RCC) is a common cancer with high malignancy, and 20–30% of RCC patients already suffer from metastatic lesions at the time of diagnosis (Basso et al., 2010; Sun et al., 2011). Sunitinib is a first-generation tyrosine kinase inhibitor (TKI) that was approved for the treatment of metastatic RCC (mRCC) in 2006 (Diekstra et al., 2016). In patients with mRCC, sunitinib has been associated with improvements in progression-free survival (PFS) and overall survival (OS) compared with interferon-alpha (Motzer et al., 2009), and it achieves significant response rates in both Western and Japanese patients (Motzer et al., 2006; Tomita et al., 2010). The Metastatic Renal Cell Carcinoma database Consortium International (IMDC) model is the most widely used prognostic models for the prognosis of mRCC in clinical practice and clinical trials (Rini et al., 2018). Currently, sunitinib remains the first-line standard of care for IMDC favorable-risk patients, and for IMDC intermediate- and poor-risk patients, immune checkpoint inhibition, sequencing, and combined systemic therapy have been reported to have an OS benefit (de Velasco et al., 2019; Loo et al., 2019; Stuhler et al., 2020). Increased exposure to sunitinib is associated with improved PFS and OS but also an increased risk for adverse events (Houk et al., 2009; Motzer et al., 2009). Sunitinib-associated toxicities include thrombocytopenia, neutropenia, leucopenia, hand-foot syndrome (HFS), hypertension, mucositis, and diarrhea (Kalantari, 2009; Escudier et al., 2014; Atkins et al., 2015). The therapeutic efficacy and toxicity of sunitinib are very heterogeneous and have been difficult to predict before treatment initiation. Single-nucleotide polymorphisms (SNPs) in genes encoding metabolism enzymes or transporters related to the pharmacokinetics (PK) and pharmacodynamics (PD) of sunitinib have been identified to be associated with the toxicity and efficacy of sunitinib in previous studies (Garcia-Donas et al., 2015; Rodriguez-Antona and Taron, 2015), especially SNPs in genes related to the PK pathways of sunitinib in patients with RCC (Junker et al., 2013; Fiszer and Krzyzosiak, 2014; Narjoz et al., 2015).
Sunitinib is a substrate of ATP-binding cassette member B1 (ABCB1) and another efflux transporter encoded by ATP-binding cassette member G2 (ABCG2) (Low et al., 2016). ABCG2 and ABCB1 are the PK-related genes of sunitinib (Watanabe et al., 2017). The most common functional SNP in ABCG2 was reported to be rs2231142 (421C/A), and those in ABCB1 were rs2032582 (2677G/TA), rs1128503 (1236T/C), and rs1045642 (3435C/T) (Kato et al., 2017; Watanabe et al., 2017; Reustle et al., 2018). ABCG2 rs2231142 was located at Q141K, ABCB1 rs1045642 at | 1154 | , ABCB1 rs1128503 at G412G, and ABCB1 rs2032582 at A893 S/T. The associations of ABCB1 and ABCG2 polymorphisms with sunitinib-induced toxicity and clinical outcomes in patients with RCC have been investigated. However, the associations were controversial. ABCG2 rs2231142 was reported to have no association with thrombocytopenia in Whites and Asians (van Erp et al., 2009; Kato et al., 2017; Zhang et al., 2018). However, in other studies, ABCG2 rs2231142 was associated with severe thrombocytopenia (Low et al., 2016), and patients carrying the ABCG2 rs2231142 AA genotype were more likely to develop thrombocytopenia, neutropenia, and HFS even after adjustment (Kim et al., 2013). Chu et al. reported that the ABCG2 rs2231142 A allele was correlated with a 3-fold decrease in the risk of neutropenia (Chu et al., 2015). Garcia-Donas et al. reported that ABCG2 rs2231142 was not associated with hypertension but seemed to confer protection against HFS, while ABCB1 rs1128503 and rs2032582 seemed to confer protection against hypertension (Garcia-Donas et al., 2011). The ABCB1 rs2032582 TT, AT, and GT genotypes were reported to be significantly correlated with grade 2 and grade 3 HFS (Zhang et al., 2018). However, de Velasco et al. reported that ABCB1 rs2032582 was not associated with sunitinib-induced toxicity (de Velasco et al., 2016). Beuselinck et al. reported that both the PFS and OS of patients with mRCC were significantly associated with SNP rs1128503 in ABCB1 but were not associated with ABCB1 rs2032582 (Beuselinck et al., 2013), while ABCB1 rs1128503 and rs2032582 were not associated with PFS or OS in patients with mRCC in another study (Garcia-Donas et al., 2011).
Identifying the effects of ABCG2 and ABCB1 polymorphisms on sunitinib-induced toxicity and efficacy in patients with RCC could help to optimize the therapeutic management strategy and maximize the clinical benefits of sunitinib. Here we used a meta-analysis to evaluate the associations of ABCG2 rs2231142, ABCB1 rs1128503, and ABCB1 rs2032582 polymorphisms with sunitinib-induced toxicity and those of ABCB1 rs1128503 and ABCB1 rs2032582 polymorphisms with PFS and OS in patients with RCC.
Methods
Search Strategy and Selection Criteria
PubMed, EMBASE, Cochrane Library, and Web of Science were systematically searched using the keywords “sunitinib’’ and ‘‘polymorphisms” to identify eligible studies. Only articles in English were included. The last search was updated on 7 February 2020. We also retrieved eligible studies in the references of all relevant articles. Furthermore, we contact the authors to request missing information to strengthen the analysis.
The inclusion criteria for studies were as follows: 1) studies of RCC patients who received sunitinib treatment; 2) studies evaluating the associations of ABCG2 rs2231142 polymorphism with sunitinib-induced thrombocytopenia, HFS, and hypertension; 3) studies evaluating the associations of ABCB1 rs1128503 and rs2032582 polymorphisms with sunitinib-induced hypertension and HFS; 4) studies evaluating the associations of ABCB1 rs1128503 and rs2032582 polymorphisms with PFS and OS; and 5) studies with available odds ratios (ORs) between genotypes.
The exclusion criteria for studies were as follows: 1) letters, reviews, and case reports; 2) nonclinical studies; and 3) studies that were duplicate publications. In addition, we pooled complete data in the meta-analysis if multiple studies had overlapping data.
Data Extraction
Two authors independently selected the relevant articles and then extracted the following data from the articles: first author’s name, publication year, study design, sample size, treatment regimen, age, sex, region/ethnicity, ABCG2 rs2231142, ABCB1 rs1128503 and ABCB1 rs2032582 phenotype, genotyping methods, evidence of Hardy-Weinberg equilibrium in cases, adverse events and clinical outcome in ABCG2 rs2231142, ABCB1 rs1128503 and ABCB1 rs2032582 phenotypes, and the ORs of ABCG2 rs2231142 A allele vs. C allele, ABCB1 rs1128503 T allele vs. C allele, and ABCB1 rs2032582 T allele vs. other alleles. If there was any controversy, it was resolved by discussion among the authors.
Statistical Analysis
STATA software (version 12.0; Stata Corp, College Station, TX, USA) was used to perform the statistical tests. The associations of ABCG2 rs2231142, ABCB1 rs1128503, and ABCB1 rs2032582 polymorphisms with sunitinib-induced toxicity, and ABCB1 rs1128503 and ABCB1 rs2032582 polymorphisms with PFS and OS were evaluated by effect size (ES) with 95% CIs (Cheng et al., 2020). The heterogeneity between studies was assessed by the chi-square-based Q-test and I2 tests. Both no-heterogeneity criteria were required to be met (p > 0.05 AND I2 < 50%) to use the fixed-effects model; otherwise, the random-effects model was used to calculate the pooled ES. The Z test was used to investigate the significance of the ES, and PZ < 0.05 was considered statistically significant. The stability of the results was evaluated by sensitivity analysis. The publication bias among studies was determined by Egger’s and Begg’s tests, and PE < 0.05 was considered significant.
Results
Study Selection
The study selection process is shown in Figure 1. We initially retrieved 815 articles from electronic databases, including PubMed, EMBASE, Science Direct, and Web of Science. We included 34 relevant studies after removing the duplicates and reviewing the titles, abstracts, and full texts. Then, 25 articles were excluded because they did not determine the ORs between genotypes. Ultimately, 8 articles published between 2011 and 2018 assessing the relationships of ABCG2 rs2231142, ABCB1 rs1128503, and ABCB1 rs2032582 polymorphisms with sunitinib-induced toxicity, and 5 articles published between 2011 and 2017 assessing the relationships of ABCB1 rs1128503 and ABCB1 rs2032582 polymorphisms with PFS and OS were included in the current meta-analysis (Garcia-Donas et al., 2011; Beuselinck et al., 2013; Kim et al., 2013; Diekstra et al., 2015; Chu et al., 2015; de Velasco et al., 2016; Low et al., 2016; Numakura et al., 2017; Zhang et al., 2018).
Characteristics of the Included Studies
The main characteristics of all eligible studies are shown in Tables 1–3. In total, 9 studies including 1081 patients were enrolled in the pooled analysis. The total sample size is the sum of the sample in each study except for Garcia-Donas et al. study, because patients in the study by Diekstra et al. were pooled from 5 studies including study of Garcia-Donas et al., and we only used data from Garcia-Donas et al. when those from Diekstra et al. did not indicate the specific toxicity or efficacy. The included studies are all observational. Five studies were prospective cohort studies, and four studies were retrospective cohort studies. Five studies (55.6%) reported on Asian individuals, and 4 (44.4%) studies reported mainly on Caucasian individuals. Sunitinib was given as monotherapy in the included studies. The variant allele A frequency of ABCG2 rs2231142 in Asians was significantly higher than that in Caucasians (34.1% vs. 5.3%), the variant allele C frequency of ABCB1 rs1128503 in the two different ethnic groups was similar (37.6% vs. 31.1%), and the variant AT genotype frequency of ABCB1 rs2032582 in Asians was a little higher than that in Caucasians (29.9% vs. 16.5%) in the included studies. The selected results contained the necessary information according to the STrengthening the REporting of Genetic Association Studies (STREGA) (Little et al., 2009).
Relationship Between the ABCG2 rs2231142 Polymorphism and Toxicity
Thinking the consistency of the overall results, we initially analyzed the associations of the ABCG2 polymorphisms with sunitinib-induced toxicities without regard to ethnic background. The meta-analysis results are shown in Table 4. There was no heterogeneity between studies assessing the relationships of the ABCG2 rs2231142 polymorphism with sunitinib-induced thrombocytopenia (I2 = 0, p > 0.05), and the fixed-effects model was applied to the analysis. There was heterogeneity between studies assessing the relationships of ABCG2 rs2231142 polymorphism with sunitinib-induced hypertension, HFS, and neutropenia (I2 = 61.4%, 68.3%, and 68.6%, respectively), and the random-effects model was applied to the analysis. Based on the results, we found that compared with the C allele, the ABCG2 rs2231142 A allele was significantly associated with an increased risk of sunitinib-induced thrombocytopenia in Asians (ES = 1.65, 95% CI = 1.15–2.36, p = 0.006; Figure 2), while there were no significant associations between the A allele and C allele in sunitinib-induced hypertension, HFS, or neutropenia (ES = 1.09, 95% CI = 0.69–1.73, p = 0.701; ES = 1.40, 95% CI = 0.94–2.08, p = 0.094; ES = 0.87, 95% CI = 0.57–1.31, p = 0.501; Figure 2). Since the patients included in the meta-analysis of thrombocytopenia and neutropenia were all Asian (Table 2), we only conducted a subgroup analysis of HFS and hypertension in Asians. Compared with the C allele, the ABCG2 rs2231142 A allele was significantly associated with an increased risk of sunitinib-induced HFS (ES = 1.52, 95% CI = 1.02–2.27, p = 0.041; Figure 3). There was also no significant associations between the A allele and C allele in sunitinib-induced hypertension in Asians (ES = 1.17, 95% CI = 0.74–1.86, p = 0.504; Figure 3).
FIGURE 2. Forest plots estimating the associations of ABCG2 rs2231142 polymorphism with sunitinib-induced toxicity and the corresponding sensitivity analysis (A allele vs. C allele). The ABCG2 rs2231142 A allele was significantly associated with an increased risk of thrombocytopenia in Asians compared with the C allele, while there were no significant associations between the A allele and C allele in hypertension, HFS, or neutropenia. The sensitivity analysis showed that no single study qualitatively altered the pooled ES of thrombocytopenia and neutropenia; however, the study by Diekstra et al. altered the hypertension results, and the study by Garcia-Donas et al. altered the HFS results.
FIGURE 3. Forest plots estimating the associations of the ABCG2 rs2231142 polymorphism with sunitinib-induced toxicities in Asians and the corresponding sensitivity analysis (A allele vs. C allele). Compared with the C allele, the ABCG2 rs2231142 A allele was significantly associated with an increased risk of HFS in Asians. There were no significant associations between the A allele and C allele in hypertension in this population. The sensitivity analysis showed that no single study qualitatively altered the pooled ES.
Correlations of the ABCB1 rs1128503 Polymorphism With Toxicity and Clinical Outcomes
We only performed the meta-analysis of rs1128503 in ABCB1 in mixed ethnicity as the variant allele C frequency in Asians and Caucasians was similar. There was no heterogeneity between studies assessing the associations of the ABCB1 rs1128503 polymorphism with sunitinib-induced hypertension (I2 = 0, p > 0.05), and the fixed-effects model was applied to the analysis. There was heterogeneity between studies assessing the relationship of ABCB1 rs1128503 polymorphism with sunitinib-induced PFS and OS (I2 = 82.4%, 68.3%, p < 0.05), and the random-effects model was applied to the analysis. The results showed that compared with the C allele, the ABCB1 rs1128503 T allele was significantly associated with a decreased risk of sunitinib-induced hypertension and worse PFS (ES = 0.44, 95% CI = 0.26–0.77, p = 0.004; ES = 1.36, 95% CI = 1.07–1.73, p = 0.011; Figure 4). There was no significant association between the T allele and C allele in sunitinib-induced OS (ES = 0.82, 95% CI = 0.61–1.10, p = 0.184; Figure 4).
FIGURE 4. Forest plots estimating the associations of the ABCB1 rs1128503 polymorphism with sunitinib-induced hypertension, progression-free survival (PFS), and overall survival (OS), and the corresponding sensitivity analysis (T allele vs. C allele). Compared with the C allele, the ABCB1 rs1128503 T allele was significantly associated with a decreased risk of hypertension and worse PFS. There was no significant association between the T allele and C allele in OS. The sensitivity analysis showed that no single study qualitatively altered the pooled ES of hypertension; however, the study by Beuselinck et al. altered the PFS results, and the study by Diekstra et al. altered the OS results.
Correlations of the ABCB1 rs2032582 Polymorphism With Toxicity and Clinical Outcomes
There was heterogeneity between studies assessing the associations of the ABCB1 rs2032582 polymorphism with sunitinib-induced hypertension and OS (I2 = 80.6%, 63.3%), and the random-effects model was applied to the analysis. There was no heterogeneity between studies assessing the relationship of the ABCB1 rs2032582 polymorphism with sunitinib-induced HFS or PFS (I2 = 48.4%, 0, p > 0.05), and the fixed-effects model was applied to the analysis. The results showed that compared with the other alleles, the ABCB1 rs2032582 T allele was significantly associated with worse PFS (ES = 1.46, 95% CI = 1.14–1.87, p = 0.003; Figure 5). In contrast, there was no significant association between the T allele and other alleles in sunitinib-induced hypertension, HFS, or OS (ES = 0.77, 95% CI = 0.46–1.29, p = 0.326; ES = 1.02, 95% CI = 0.65–1.62, p = 0.919; ES = 1.32, 95% CI = 0.85–2.05, p = 0.215; Figure 5).
FIGURE 5. Forest plots estimating the associations of the ABCB1 rs2032582 polymorphism with sunitinib-induced hypertension, hand-foot syndrome (HFS), progression-free survival (PFS), and overall survival (OS), and the corresponding sensitivity analysis (T allele vs. other alleles). Compared with the other alleles, the ABCB1 rs2032582 T allele was significantly associated with worse PFS. In contrast, there was no significant association between the T allele and other alleles in sunitinib-induced hypertension, HFS, or OS. The sensitivity analysis showed that no single study qualitatively altered the pooled ES.
Sensitivity Analysis and Publication Bias
The sensitivity analysis in the meta-analysis of the relationships of the ABCG2 rs2231142 polymorphism with sunitinib-induced thrombocytopenia and neutropenia showed that no single study qualitatively altered the pooled ES, however, the study by Diekstra et al. may qualitatively altered the result of hypertension and that by Garcia-Donas et al. may qualitatively altered the result of HFS (Figure 2). When we omitted the article, the meta-analysis results for HFS changed but not the hypertension results, and no single study qualitatively altered the pooled ES in either analysis, which provided evidence of the stability of the meta-analysis (Figure 3). There was no publication bias for thrombocytopenia and HFS (p > 0.05, Table 4).
The sensitivity analysis in the meta-analysis of the relationships of the ABCB1 rs1128503 polymorphism with sunitinib-induced hypertension showed that no single study qualitatively altered the pooled ES; however, the study by Beuselinck et al. may qualitatively altered the result of PFS and that by Diekstra et al. may qualitatively altered the result of OS (Figure 4). When we omitted the article, the meta-analysis results for PFS (ES = 1.58, 95% CI = 1.23–2.04, p < 0.001) and OS (ES = 1.27, 95% CI = 0.77–2.11, p = 0.353) were not changed, and no single study qualitatively altered the pooled ES in either analysis (Figure 6).
FIGURE 6. Forest plots estimating the associations of ABCB1 rs1128503 polymorphism with sunitinib-induced progression-free survival (PFS) and overall survival (OS) when we omitted the article that altered the pooled ES (T allele vs. C allele). The ABCB1 rs1128503 polymorphism was associated with PFS but not OS, and no single study qualitatively altered the pooled ES in both analysis.
The sensitivity analysis in the meta-analysis of the relationships of the ABCB1 rs2032582 polymorphism with sunitinib-induced hypertension, HFS, PFS, and OS showed that no single study qualitatively altered the pooled ES (Figure 5). There was no publication bias for PFS in ABCB1 rs1128503, or ABCB1 rs2032582 polymorphisms (p > 0.05, Table 4).
Discussion
In the current meta-analysis, we initially analyzed the associations of the ABCG2 rs2231142 polymorphism with sunitinib-induced thrombocytopenia, HFS, hypertension, and neutropenia; the associations of the ABCB1 rs1128503 polymorphism with sunitinib-induced hypertension, PFS, and OS; and the associations of ABCB1 rs2032582 polymorphisms with sunitinib-induced hypertension, HFS, PFS, and OS, without regard to ethnic background. However, we found heterogeneity between studies based on the results. We further considered the origin of the heterogeneity, and we found that the heterogeneity might derive from ethnicity as the variant allele A frequency of ABCG2 rs2231142 in Asians was significantly higher than that in Caucasians, then we conducted a subgroup analysis of hypertension and HFS in Asians. Next, we did the sensitivity analysis of all meta-analysis to evaluate their stability. The results showed that one article with population Caucasians altered the results of HFS in ABCG2, which indicated that the heterogeneity might derive from ethnicity. For ABCB1 rs1128503 and ABCB1 rs2032582 results, there were no article qualitatively altered the results. The sensitivity analysis results provided evidence of the stability of the meta-analysis.
The therapeutic efficacy of sunitinib is heterogeneous among different IMDC risk groups (Rini et al., 2018). In studies included in the meta-analysis, they did not mention the IMDC stratification, and we could not analyze the effect of IMDC on toxicity and clinical outcomes. The treatment regimen of sunitinib in the included articles was not completely consistent, and five studies had three starting dose (50 or 37.5 or ≤37.5, mg/day) (Garcia-Donas et al., 2011; Diekstra et al., 2015; Chu et al., 2015; Numakura et al., 2017; Zhang et al., 2018). In Chu et al. study, the starting dose (50 or 37.5 or ≤37.5, mg/day) was not associated with toxicities, including leukopenia, neutropenia, or diarrhea but was associated with PFS and OS (Chu et al., 2015), and we used the results of multivariate analysis in the meta-analysis, which showed no significant association. Zhang et al. study showed no significant association between toxicity of grade <3 or ≥3 and each dosage (25, 37.5, and 50 mg/d) (Zhang et al., 2018). The effect of dose on toxicities and clinical outcome were not mentioned in the other three studies (Garcia-Donas et al., 2011; Diekstra et al., 2015; Numakura et al., 2017). RCC is more common in men. However, the included studies did not investigate the differences between these genetic polymorphisms in different genders. In the included studies, Chu et al. reported that the gender was not associated with the leucopenia, neutropenia, diarrhea, PFS or OS (Chu et al., 2015). In Diekstra et al. study, gender seemed to show association with thrombocytopenia and toxicity grade >2 in univariate analysis, however, there was no association in the multivariate analysis (Diekstra et al., 2015). Zhang et al. showed that gender was associated with hypertension in univariate analysis, but the results of multivariate analysis showed no significant association (Zhang et al., 2018). In Low et al. study, gender was only associated with leucopenia (Low et al., 2016), which was not included in the meta-analysis. In summary, we used ORs of multivariate analysis in the pooled ES analysis, excluding other possible factors affecting toxicity and efficacy of sunitinib, such as gender, age, sunitinib dosage, laboratory indexes, and other related gene polymorphisms in each included study, which strengthened our conclusions.
Compared with non-Asians, Asians are more likely to develop sunitinib-induced adverse effects (Kim et al., 2011; Lee et al., 2014). The incidence of sunitinib-induced hematotoxicity in Japanese patients with RCC is higher when compared to European populations (Low et al., 2016; Kato et al., 2017). The AA genotype of rs2231142 in ABCG2 has been associated with an increase in systemic exposure to sunitinib, possibly causing sunitinib-related adverse events. In the current meta-analysis, we found that compared with the C allele, the ABCG2 rs2231142 A allele was significantly associated with an increased risk of sunitinib-induced thrombocytopenia in Asians. However, we could not find associations of ABCG2 rs2231142 polymorphisms with hypertension and neutropenia in the meta-analysis.
Garcia-Donas et al. reported that the A allele of ABCG2 rs2231142 seemed to confer protection against HFS in Caucasians (Garcia-Donas et al., 2011), while in the other four studies carried out in Asians, an increased risk of development of HFS was found (Kim et al., 2013; Low et al., 2016; Numakura et al., 2017; Zhang et al., 2018). In our results, there was no significant association between the A allele and C allele in sunitinib-induced HFS. However, when we conducted a subgroup analysis in Asians, we found that compared with the C allele, the ABCG2 rs2231142 A allele was significantly associated with an increased risk of sunitinib-induced HFS. A possible reason for this difference might be the small number of patients with homozygous alleles among Caucasians. The ABCG2 rs2231142 A allele is located within the ATP-binding cassette domain that regulates the ATP binding activity of ABCG2 protein, which reduces ATPase activity (Mizuarai et al., 2004), induces the reduction of transport activity, increases drug accumulation, decreases the efflux velocity of the drug (Imai et al., 2002), and thus increases the systemic exposure to sunitinib and is subsequently more likely to lead to the development of toxicity (Mizuno et al., 2012). It has been reported that the variant allele A of ABCG2 rs2231142 in Asians is approximately 30%, in Caucasians is approximately 10%, and in African Americans is approximately 5% (Diekstra et al., 2016). Diekstra et al. reported that the ABCG2 rs2231142 polymorphism was significantly associated with sunitinib-induced dose reduction after cycle 1, 2, or 3 (Diekstra et al., 2015). Although other previous studies have reported that the ABCG2 rs2231142 polymorphism is not significantly associated with toxicity-related dose reduction (Garcia-Donas et al., 2011; Kato et al., 2017; Zhang et al., 2018) or time to dose reduction in mRCC patients (Numakura et al., 2017), further studies should verify whether dose adjustment based on early onset thrombocytopenia and HFS in Asians prolongs sunitinib treatment.
Patients carrying the variant genotypes of rs1128503 and rs2032582 in ABCB1 showed an increased clearance of sunitinib and its active metabolite SU12662, consequently leading to lower exposure to sunitinib, reduced hypertension, and decreases in PFS and OS (Garcia-Donas et al., 2011; Beuselinck et al., 2014; Diekstra et al., 2014; Teo et al., 2015). In the current meta-analysis, we found that compared with the C allele, the ABCB1 rs1128503 T allele was significantly associated with a decreased risk of sunitinib-induced hypertension but worse PFS in Asians and Caucasians; and the ABCB1 rs2032582 T allele was significantly associated with worse PFS compared with other alleles, which was consistent with previous studies (Diekstra et al., 2015; Chu et al., 2015; Numakura et al., 2017). However, there was no significant association between the ABCB1 rs2032582 T and G alleles in sunitinib-induced hypertension, HFS, or OS. Four included studies also investigated the association of haplotype ABCB1 (rs1128503, rs2032582, rs1045642) with sunitinib induced toxicities and clinical outcome. In Kim et al. study, there were no significant differences among CGC, TTT, and TGC haplotype of ABCB1 (rs1128503, rs2032582, rs1045642) in thromobocytopenia, neutropenia, anemia, or HFS (Kim et al., 2013). Diekstra et al. reported that PFS and OS were improved in the presence of CGT in haplotype ABCB1 (rs1128503, rs2032582, rs1045642) (Diekstra et al., 2015), however, there was no significant difference in median PFS between the present haplotype and absent haplotype in another study (Beuselinck et al., 2013). In Chu et al. study, ABCB1 rs1045642, 1236 rs1128503, rs2032582 TTT haplotype was correlated with a 10-fold (p = 0.03) decrease in the risk of neutropenia and inferior PFS and OS (Chu et al., 2015). We should also pay attention to patients carrying the opposite haplotype.
SNPs in sunitinib target candidate genes, including vascular endothelial growth factor receptors (VEGFRs) 1, 2, and 3; Fms-like tyrosine kinase 3 receptor (FLT3); and PK-related genes, including cytochrome P450 1A1 (CYP1A1), CYP3A5, NR1/2, and NR1/3, have also been reported to be associated with sunitinib-induced toxicity and efficacy in patients with mRCC (van Erp et al., 2009; Garcia-Donas et al., 2011; van der Veldt et al., 2011; Eechoute et al., 2012; Kim et al., 2012; Beuselinck et al., 2013; Beuselinck et al., 2014; Kato et al., 2017; Watanabe et al., 2017). SNPs rs2010963 and rs2070744 in VEGF were reported to be associated with increased chances for the occurrence and duration of hypertension (Eechoute et al., 2012; Kim et al., 2012). The FLT3 738 T/C polymorphism and the CYP3A5 A allele (CYP3A5*1) of rs776746 were reported to be associated with dose reduction due to sunitinib-induced toxicity (Garcia-Donas et al., 2011; Kato et al., 2017). Van Erp et al. observed an increased risk for leukopenia in CYP1A1 rs1048943 and FLT3 rs1933437 and an absence of CAG in the NR1/3 haplotype (rs2307424, rs2307418, rs4073054); CYP1A1 rs1048943 was associated with mucosal inflammation; and EGFR2 rs2305948 was associated with any toxicity > grade 2 (van Erp et al., 2009). Watanabe A et al. reported that the STAT3 polymorphism contributes to a risk factor for stomatitis (Watanabe et al., 2017). The rs9582036, rs9582036, and rs9554320 in VEGF-R1; rs2981582 and rs2305948 in VEGF-R2; rs307826 and rs307821 in VEGF-R3; and rs3025039 in VEGF-A were reported to be associated with PFS and OS in mRCC patients (Garcia-Donas et al., 2011; Kim et al., 2012; Beuselinck et al., 2013; Beuselinck et al., 2014). Compared with the disadvantageous genetic profile carriers, carriers of the genetic profile (at least an A allele in CYP3A5, a TCG copy in ABCB1, or a missing CAT copy in NR1/3) showed improved PFS and OS (van der Veldt et al., 2011). SNP rs2276707 in NR1/2 and SNPs rs2307424 and rs4073054 in NR1/3 were also reported to be associated with PFS and OS in patients with mRCC (Beuselinck et al., 2013).
We did not include rs1045642 in ABCB1 in the meta-analysis because we could not obtain the ORs between the genotypes or because the included article numbers were below three. ABCG2 polymorphisms have also been reported to be associated with other sunitinib-induced toxicity in mRCC patients, such as hypothyroidism (Werbrouck et al., 2018). The T alleles of ABCB1 rs1128503 and ABCB1 rs2032582 were reported to be correlated with decreases in the risk of neutropenia and diarrhea, respectively (Chu et al., 2015). We did not include all sunitinib-induced toxicities or SNPs associated with sunitinib-induced toxicity and efficacy for the above reason, or there was no controversy.
To our knowledge, the present study is the first meta-analysis to investigate the associations of ABCB1 and ABCG2 polymorphisms with sunitinib-induced toxicity and efficacy in patients with RCC. Based on our results, compared with the C allele, the ABCG2 rs2231142 allele A was significantly associated with increased risks of sunitinib-induced thrombocytopenia and HFS in Asians, while the T alleles of ABCB1 rs1128503 and ABCB1 rs2032582 were significantly associated with worse PFS than the other alleles. The results indicate that the ABCG2 rs2231142 polymorphism may serve as a predictor of sunitinib-induced thrombocytopenia and HFS in Asians, while the ABCB1 rs1128503 polymorphism may serve as a predictor of sunitinib-induced hypertension, and both the ABCB1 rs1128503 and rs2032582 polymorphisms may serve as predictors of PFS in RCC. Although our analysis showed no significant relationships of ABCB1 and ABCG2 polymorphisms with other sunitinib-induced severe toxicity or OS, we should also pay more attention to the use of sunitinib in Asian patients. The results may support the possible application of individualized use of sunitinib according to the genetic background of patients with RCC. Genotyping for ABCG2 rs2231142, ABCB1 rs1128503 and rs2032582 polymorphism could become a clinical routine practice to select the appropriate dose to decrease the risk of sunitinib-induced thrombocytopenia and HFS in Asians while ensuring efficacy. The population pharmacokinetic model based on the genotypes should be established to predict the occurrence of sunitinib-induced toxicities and the efficacy.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
LC and FS designed the study, performed the literature search and screening, performed the data analysis, and drafted the manuscript. ZC performed the literature search and screening and extracted the literature data. PY and ZL assisted in extracting the literature data. BW assisted in the figure design. All authors approved the final version of the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Atkins, M. B., Gravis, G., Drosik, K., Demkow, T., Tomczak, P., Wong, S. S., et al. (2015). Trebananib (AMG 386) in combination with sunitinib in patients with metastatic renal cell cancer: an open-label, multicenter, phase II study. J. Clin. Oncol. 33 (30), 3431–3438. doi:10.1200/JCO.2014.60.6012
Basso, M., Cassano, A., and Barone, C. (2010). A survey of therapy for advanced renal cell carcinoma. Urol. Oncol. 28 (2), 121–133. doi:10.1016/j.urolonc.2009.04.015
Beuselinck, B., Karadimou, A., Lambrechts, D., Claes, B., Wolter, P., Couchy, G., et al. (2013). Single-nucleotide polymorphisms associated with outcome in metastatic renal cell carcinoma treated with sunitinib. Br. J. Cancer 108 (4), 887–900. doi:10.1038/bjc.2012.548
Beuselinck, B., Lambrechts, D., Van Brussel, T., Wolter, P., Cardinaels, N., Joniau, S., et al. (2014). Efflux pump ABCB1 single nucleotide polymorphisms and dose reductions in patients with metastatic renal cell carcinoma treated with sunitinib. Acta Oncol. 53 (10), 1413–1422. doi:10.3109/0284186x.2014.918276
Cheng, L., Wang, Y., Li, X., Feng, W., Weng, B., Yuan, Q., et al. (2020). Meta-analysis of the associations of CYP2B6-516G>T polymorphisms with efavirenz-induced central nervous system side effects and virological outcome in HIV-infected adults. Pharmacogenomics J. 20 (2), 246–259. doi:10.1038/s41397-019-0112-2
Chu, Y. H., Li, H., Tan, H., Koh, V., Lai, J., Phyo, W. M., et al. (2015). Association of ABCB1 and FLT3 polymorphisms with toxicities and survival in asian patients receiving sunitinib for renal cell carcinoma. PLos One 10 (8), e0134102. doi:10.1371/journal.pone.0134102
de Velasco, G., Gray, K. P., Hamieh, L., Urun, Y., Carol, H. A., Fay, A. P., et al. (2016). Pharmacogenomic markers of targeted therapy toxicity in patients with metastatic renal cell carcinoma. Eur. Urol. Focus 2 (6), 633–639. doi:10.1016/j.euf.2016.03.017
de Velasco, G., Bex, A., Albiges, L., Powles, T., Rini, B. I., Motzer, R. J., et al. (2019). Sequencing and combination of systemic therapy in metastatic renal cell carcinoma. Eur. Urol. Oncol. 2 (5), 505–514. doi:10.1016/j.euo.2019.06.022
Diekstra, M. H., Klümpen, H. J., Lolkema, M. P., Yu, H., Kloth, J. S., Gelderblom, H., et al. (2014). Association analysis of genetic polymorphisms in genes related to sunitinib pharmacokinetics, specifically clearance of sunitinib and SU12662. Clin. Pharmacol. Ther. 96 (1), 81–89. doi:10.1038/clpt.2014.47
Diekstra, M. H., Swen, J. J., Boven, E., Castellano, D., Gelderblom, H., Mathijssen, R. H., et al. (2015). CYP3A5 and ABCB1 polymorphisms as predictors for sunitinib outcome in metastatic renal cell carcinoma. Eur. Urol. 68 (4), 621–629. doi:10.1016/j.eururo.2015.04.018
Diekstra, M. H., Swen, J. J., Gelderblom, H., and Guchelaar, H. J. (2016). A decade of pharmacogenomics research on tyrosine kinase inhibitors in metastatic renal cell cancer: a systematic review. Expert Rev. Mol. Diagn. 16 (5), 605–618. doi:10.1586/14737159.2016.1148601
Eechoute, K., van der Veldt, A. A., Oosting, S., Kappers, M. H., Wessels, J. A., Gelderblom, H., et al. (2012). Polymorphisms in endothelial nitric oxide synthase (eNOS) and vascular endothelial growth factor (VEGF) predict sunitinib-induced hypertension. Clin. Pharmacol. Ther. 92 (4), 503–510. doi:10.1038/clpt.2012.136
Escudier, B., Porta, C., Bono, P., Powles, T., Eisen, T., Sternberg, C. N., et al. (2014). Randomized, controlled, double-blind, cross-over trial assessing treatment preference for pazopanib versus sunitinib in patients with metastatic renal cell carcinoma: PISCES Study. J. Clin. Oncol. 32 (14), 1412–1418. doi:10.1200/JCO.2013.50.8267
Fiszer, A., and Krzyzosiak, W. J. (2014). Oligonucleotide-based strategies to combat polyglutamine diseases. Nucleic Acids Res. 42 (11), 6787–6810. doi:10.1093/nar/gku385
Garcia-Donas, J., Esteban, E., Leandro-García, L. J., Castellano, D. E., del Alba, A. G., Climent, M. A., et al. (2011). Single nucleotide polymorphism associations with response and toxic effects in patients with advanced renal-cell carcinoma treated with first-line sunitinib: a multicentre, observational, prospective study. Lancet Oncol. 12 (12), 1143–1150. doi:10.1016/s1470-2045(11)70266-2
Garcia-Donas, J., Rodriguez-Moreno, J. F., Romero-Laorden, N., and Rodriguez-Antona, C. (2015). Renal carcinoma pharmacogenomics and predictors of response: steps toward treatment individualization. Urol. Oncol. 33 (4), 179–186. doi:10.1016/j.urolonc.2013.09.015
Houk, B. E., Bello, C. L., Kang, D., and Amantea, M. (2009). A population pharmacokinetic meta-analysis of sunitinib malate (SU11248) and its primary metabolite (SU12662) in healthy volunteers and oncology patients. Clin. Cancer Res. 15 (7), 2497–2506. doi:10.1158/1078-0432.CCR-08-1893
Imai, Y., Nakane, M., Kage, K., Tsukahara, S., Ishikawa, E., Tsuruo, T., et al. (2002). C421A polymorphism in the human breast cancer resistance protein gene is associated with low expression of Q141K protein and low-level drug resistance. Mol. Cancer Ther. 1 (8), 611–616.
Junker, K., Ficarra, V., Kwon, E. D., Leibovich, B. C., Thompson, R. H., and Oosterwijk, E. (2013). Potential role of genetic markers in the management of kidney cancer. Eur. Urol. 63 (2), 333–340. doi:10.1016/j.eururo.2012.09.040
Kalantari, H. R. (2009). Predictive factors for severe toxicity of sunitinib in unselected patients with advanced renal cell cancer. Br. J. Cancer 101 (7), 1222–1224. doi:10.1038/sj.bjc.660530310.1038/sj.bjc.6605303
Kato, R., Kato, Y., Matsuura, T., Kanehira, M., Takata, R., and Obara, W. (2017). Characteristics of early-onset hematotoxicity of sunitinib in Japanese patients with renal cell carcinoma. BMC Cancer 17 (1), 214. doi:10.1186/s12885-017-3205-9
Kim, H. R., Park, H. S., Kwon, W. S., Lee, J. H., Tanigawara, Y., Lim, S. M., et al. (2013). Pharmacogenetic determinants associated with sunitinib-induced toxicity and ethnic difference in Korean metastatic renal cell carcinoma patients. Cancer Chemother. Pharmacol. 72 (4), 825–835. doi:10.1007/s00280-013-2258-y
Kim, H. S., Hong, M. H., Kim, K., Shin, S. J., Ahn, J. B., Jeung, H. C., et al. (2011). Sunitinib for Asian patients with advanced renal cell carcinoma: a comparable efficacy with different toxicity profiles. Oncology 80 (5-6), 395–405. doi:10.1159/000330361
Kim, J. J., Vaziri, S. A., Rini, B. I., Elson, P., Garcia, J. A., Wirka, R., et al. (2012). Association of VEGF and VEGFR2 single nucleotide polymorphisms with hypertension and clinical outcome in metastatic clear cell renal cell carcinoma patients treated with sunitinib. Cancer 118 (7), 1946–1954. doi:10.1002/cncr.26491
Lee, S. H., Bang, Y. J., Mainwaring, P., Ng, C., Chang, J. W., Kwong, P., et al. (2014). Sunitinib in metastatic renal cell carcinoma: an ethnic Asian subpopulation analysis for safety and efficacy. Asia Pac. J. Clin. Oncol. 10 (3), 237–245. doi:10.1111/ajco.12163
Little, J., Higgins, J. P., Ioannidis, J. P., Moher, D., Gagnon, F., von Elm, E., et al. (2009). STrengthening the REporting of genetic association studies (STREGA): an extension of the STROBE statement. PLos Med. 6 (2), e22. doi:10.1371/journal.pmed.1000022
Loo, V., Salgia, M., Bergerot, P., Philip, E. J., and Pal, S. K. (2019). First-line systemic therapy for metastatic clear-cell renal cell carcinoma: critical appraisal of emerging options. Target. Oncol. 14 (6), 639–645. doi:10.1007/s11523-019-00676-y
Low, S. K., Fukunaga, K., Takahashi, A., Matsuda, K., Hongo, F., Nakanishi, H., et al. (2016). Association study of a functional variant on ABCG2 gene with sunitinib-induced severe adverse drug reaction. PLoS One 11 (2), e0148177. doi:10.1371/journal.pone.0148177
Mizuarai, S., Aozasa, N., and Kotani, H. (2004). Single nucleotide polymorphisms result in impaired membrane localization and reduced atpase activity in multidrug transporter ABCG2. Int. J. Cancer 109 (2), 238–246. doi:10.1002/ijc.11669
Mizuno, T., Fukudo, M., Terada, T., Kamba, T., Nakamura, E., Ogawa, O., et al. (2012). Impact of genetic variation in breast cancer resistance protein (BCRP/ABCG2) on sunitinib pharmacokinetics. Drug Metab. Pharmacokinet. 27 (6), 631–639. doi:10.2133/dmpk.dmpk-12-rg-026
Motzer, R. J., Rini, B. I., Bukowski, R. M., Curti, B. D., George, D. J., Hudes, G. R., et al. (2006). Sunitinib in patients with metastatic renal cell carcinoma. JAMA 295 (21), 2516–2524. doi:10.1001/jama.295.21.2516
Motzer, R. J., Hutson, T. E., Tomczak, P., Michaelson, M. D., Bukowski, R. M., Oudard, S., et al. (2009). Overall survival and updated results for sunitinib compared with interferon alfa in patients with metastatic renal cell carcinoma. J. Clin. Oncol. 27 (22), 3584–3590. doi:10.1200/JCO.2008.20.1293
Narjoz, C., Cessot, A., Thomas-Schoemann, A., Golmard, J. L., Huillard, O., Boudou-Rouquette, P., et al. (2015). Role of the lean body mass and of pharmacogenetic variants on the pharmacokinetics and pharmacodynamics of sunitinib in cancer patients. Invest. New Drugs 33 (1), 257–268. doi:10.1007/s10637-014-0178-2
Numakura, K., Tsuchiya, N., Kagaya, H., Takahashi, M., Tsuruta, H., Inoue, T., et al. (2017). Clinical effects of single nucleotide polymorphisms on drug-related genes in Japanese metastatic renal cell carcinoma patients treated with sunitinib. Anticancer Drugs 28 (1), 97–103. doi:10.1097/CAD.0000000000000425
Reustle, A., Fisel, P., Renner, O., Büttner, F., Winter, S., Rausch, S., et al. (2018). Characterization of the breast cancer resistance protein (BCRP/ABCG2) in clear cell renal cell carcinoma. Int. J. Cancer 143 (12), 3181–3193. doi:10.1002/ijc.31741
Rini, B. I., Hutson, T. E., Figlin, R. A., Lechuga, M. J., Valota, O., Serfass, L., et al. (2018). Sunitinib in patients with metastatic renal cell carcinoma: clinical outcome according to international metastatic renal cell carcinoma database Consortium risk group. Clin. Genitourin Cancer 16 (4), 298–304. doi:10.1016/j.clgc.2018.04.005
Rodríguez-Antona, C., and Taron, M. (2015). Pharmacogenomic biomarkers for personalized cancer treatment. J. Intern. Med. 277 (2), 201–217. doi:10.1111/joim.12321
Stühler, V., Maas, J. M., Rausch, S., Stenzl, A., and Bedke, J. (2020). Immune checkpoint inhibition for the treatment of renal cell carcinoma. Expert Opin. Biol. Ther. 20 (1), 83–94. doi:10.1080/14712598.2020.1677601
Sun, M., Thuret, R., Abdollah, F., Lughezzani, G., Schmitges, J., Tian, Z., et al. (2011). Age-adjusted incidence, mortality, and survival rates of stage-specific renal cell carcinoma in North America: a trend analysis. Eur. Urol. 59 (1), 135–141. doi:10.1016/j.eururo.2010.10.029
Teo, Y. L., Wee, H. L., Chue, X. P., Chau, N. M., Tan, M. H., Kanesvaran, R., et al. (2015). Effect of the CYP3A5 and ABCB1 genotype on exposure, clinical response and manifestation of toxicities from sunitinib in Asian patients. Pharmacogenomics J. 16 (1), 47–53. doi:10.1038/tpj.2015.13
Tomita, Y., Shinohara, N., Yuasa, T., Fujimoto, H., Niwakawa, M., Mugiya, S., et al. (2010). Overall survival and updated results from a phase II study of sunitinib in Japanese patients with metastatic renal cell carcinoma. Jpn. J. Clin. Oncol. 40 (12), 1166–1172. doi:10.1093/jjco/hyq146
van der Veldt, A. A., Eechoute, K., Gelderblom, H., Gietema, J., Guchelaar, H. J., van Erp, N. P., et al. (2011). Genetic polymorphisms associated with a prolonged progression-free survival in patients with metastatic renal cell cancer treated with sunitinib. Clin. Cancer Res. 17 (3), 620–629. doi:10.1158/1078-0432.CCR-10-1828
van Erp, N. P., Eechoute, K., van der Veldt, A. A., Haanen, J. B., Reyners, A. K., Mathijssen, R. H., et al. (2009). Pharmacogenetic pathway analysis for determination of sunitinib-induced toxicity. J. Clin. Oncol. 27 (26), 4406–4412. doi:10.1200/JCO.2008.21.7679
Watanabe, A., Yamamoto, K., Ioroi, T., Hirata, S., Harada, K., Miyake, H., et al. (2017). Association of single nucleotide polymorphisms in STAT3, ABCB1, and ABCG2 with stomatitis in patients with metastatic renal cell carcinoma treated with sunitinib: a retrospective analysis in Japanese patients. Biol. Pharm. Bull. 40 (4), 458–464. doi:10.1248/bpb.b16-00875
Werbrouck, E., Bastin, J., Lambrechts, D., Verbiest, A., Van Brussel, T., Lerut, E., et al. (2019). ABCG2 Polymorphism rs2231142 and hypothyroidism in metastatic renal cell carcinoma patients treated with sunitinib. Acta Clin. Belg. 74 (3), 180–188. doi:10.1080/17843286.2018.1477229
Keywords: renal cell carcinoma, sunitinib, ABCG2, ABCB1, polymorphism, meta-analysis
Citation: Sun F, Chen Z, Yao P, Weng B, Liu Z and Cheng L (2021) Meta-Analysis of ABCG2 and ABCB1 Polymorphisms With Sunitinib-Induced Toxicity and Efficacy in Renal Cell Carcinoma. Front. Pharmacol. 12:641075. doi: 10.3389/fphar.2021.641075
Received: 13 December 2020; Accepted: 19 January 2021;
Published: 08 March 2021.
Edited by:
Moneeza Kalhan Siddiqui, University of Dundee, United KingdomReviewed by:
Salvador F. Aliño, University of Valencia, SpainJatinder K. Lamba, University of Florida, United States
Copyright © 2021 Sun, Chen, Yao, Weng, Liu and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lin Cheng, cheng7zhu@163.com
†These authors have contributed equally to this work