MiR-3168, miR-6125, and miR-4718 as Potential Predictors of Cisplatin-Induced Nephrotoxicity in Patients with Head and Neck Cancer

doi:10.21203/rs.3.rs-147058/v1

Download PDF

Research Article

MiR-3168, miR-6125, and miR-4718 as Potential Predictors of Cisplatin-Induced Nephrotoxicity in Patients with Head and Neck Cancer

https://doi.org/10.21203/rs.3.rs-147058/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this older preprint version

Read the latest preprint version →

No biomarker is available for identifying cancer patients at risk of developing nephrotoxicity when treated with cisplatin. We performed microRNA (miRNA) sequencing using plasma collected five days after cisplatin treatment (D5) from twelve patients with head and neck cancer with and without nephrotoxicity (grade ≥2 increased serum creatinine). The most differentially expressed miRNAs between the two groups were selected for quantification at baseline and D5 in a larger cohort of patients. The association between miRNAs and nephrotoxicity was evaluated by calculating the odds ratio (OR) from univariate logistic regression. Receiver operating characteristic curves (ROC) were used to estimate the area under the curve (AUC), sensitivity, and specificity. MiR-3168 (p=1.98×10⁻⁸), miR-4718 (p=4.24×10⁻⁵), and miR-6125 (p=6.60×10⁻⁵) were the most differentially expressed miRNAs and were further quantified in 43, 48, and 53 patients, respectively. The baseline expression of miR-3168 (p=0.0456, OR=1.03, 95% CI: 1.00–1.06) and miR-4718 (p=0.0388, OR=1.56, 95% CI: 1.03–2.46) were associated with an increased risk of nephrotoxicity, whereas miR-6125 showed a trend (p=0.0618, OR=1.73, 95% CI: 0.98–3.29). MiR-4718 showed the highest AUC (0.77, 95% CI: 1.74–5.8) with sensitivity of 66.76 and specificity of 79.49. We therefore report baseline plasma miRNAs as potential predictors of cisplatin-induced nephrotoxicity.

Oncology

Cisplatin

nephrotoxicity

microRNAs

miR-3168

miR-6125

miR-4718

Cisplatin is a widely used chemotherapeutic agent for treating several types of tumors, including head and neck cancer[1]. The antitumor efficacy of cisplatin relies on its interaction with DNA and the formation of DNA adducts through covalent bonds, mainly at N7 of adenine and guanine. The DNA adducts arrest the cell cycle in the G2 phase and block DNA replication, leading to the activation of apoptosis in tumor cells[2].

Despite the proven efficacy of cisplatin treatment, patients frequently experience toxicity that limits the duration and efficacy of therapy. Nephrotoxicity is one of the most frequent and dose-limiting cisplatin-induced toxicities, with a prevalence that varies between 30% and 50% in patients treated with ≥80 mg/m² cisplatin[3,4].

The mechanisms of cisplatin-induced nephrotoxicity are not entirely understood. The proposed mechanisms include the following: (1) activation of apoptosis due to cisplatin-induced damage in renal tubular cells[5]; (2) induction of oxidative stress resulting from increased production of reactive oxygen species due to cisplatin-induced damage in mitochondrial DNA[3]; (3) accumulation of cisplatin in renal tubular cells due to the high affinity of cisplatin to the copper transporter CTR1 and the organic cation transporter 2 (OCT2)[6,7]; and (4) inhibition of Na⁺/K⁺ ATPase pump in renal cells by cisplatin[8].

Currently, no validated biomarkers are available for identifying patients at risk of cisplatin-induced nephrotoxicity and traditional markers used to evaluate nephrotoxicity after cisplatin administration, including serum creatinine (SCr) and blood urea nitrogen (BUN), lack specificity and sensitivity[9]. SCr is the most used indicator for detecting renal damage, but it is not an optimal renal maker, showing no significant changes in serum levels until approximately half of the nephrons are lost[10]. Thus, the need for identifying new biomarkers to evaluate cisplatin-induced nephrotoxicity.

Circulating microRNAs (miRNAs) have been investigated as potential biomarkers of cisplatin-induced nephrotoxicity because they are easily detected in biological fluids[11]. MiRNAs are small, single-stranded RNA molecules of approximately 22 non-coding nucleotides, which participate in the post-transcriptional regulation of gene expression[12]. Most studies were carried out using animal models[11]. Only one study was performed in urine samples of cisplatin-treated patients with mesothelioma to investigate specific miRNAs (miR-21, miR-200c, and miR-423) but could not validate a marker of cisplatin-induced nephrotoxicity[13]. Various circulating miRNAs are found in specific fluids; however, the same miRNAs can be found in different fluids. The plasma has the most variety of miRNAs, probably because it captures miRNAs from different cell types during circulation[14].

To our knowledge, this is the first study aimed to identify circulating plasma miRNAs as biomarkers of cisplatin-induced nephrotoxicity using the patients’ samples. We performed miRNA sequencing using plasma samples from patients with head and neck cancer, with and without cisplatin-induced nephrotoxicity, and tested three selected miRNAs for validation in a larger cohort of patients. Moreover, we performed a bioinformatics analysis to identify the mechanisms of miRNAs and their target genes to explain why only a subset of patients develop cisplatin-induced nephrotoxicity.

Patient characteristics and toxicities

A total of 60 patients were included in this study. Nephrotoxicity (grade ≥2 increased SCr) was observed in 11 (18.3%) patients and grade 3 in 3 (5.0%) patients. No case of grade 4 increased SCr was reported. Twelve out of 60 patients were selected for miRNAs sequencing and divided into two groups: 6 patients with grade ≥2 increased SCr (nephrotoxicity group) and six patients with grade = 0 increased SCr (non-nephrotoxicity group). The clinical characteristics and nephrotoxicity parameters of all patients included in the study and patients whose miRNA samples were sequenced are shown in Table 1.

MiRNA sequencing results

Thirty-three miRNAs were differentially expressed between the nephrotoxicity and non-nephrotoxicity groups with a fold regulation (FR) >2 or FR < −2, with no comments or comment “A” according to the GeneGlobe Data Analysis Center (Qiagen) (Supplementary Table S1). The most differentially expressed miRNAs were miR-3168 (p = 1.98 × 10⁻⁸, FR = 8.08), miR-6125 (p = 6.60 × 10⁻⁵, FR = 5.31), and miR-4718 (p = 4.24 × 10⁻⁵, FR = 5.12) upregulated in the nephrotoxicity group. Volcano plots for all miRNAs detected through sequencing (including those with comments “B” and “C”) are shown in Supplementary Fig. S1.

Validation of miR-3168, miR-6125, and miR-4718

The expression of miR-3168, miR-6125, and miR-4718 was verified in 50, 60, and 55 patients, respectively. However, 7 patients were excluded because cel-miR-39 expression was higher than 2 standard deviations (SD, 2 patients with grade ≥2 increased SCr and 5 patients with grade <2 increased SCr). Thus, 43, 53, and 48 patients were analyzed for the expressions of miR-3168, miR-6125, and miR-4718, respectively.

MiR-4718 expression was higher before cisplatin administration in patients with grade ≥2 increased SCr (p = 0.0433). MiR-3168 (p = 0.0570) and miR-6125 (p = 0.0658) also showed a trend toward higher expression in baseline samples of patients with grade ≥2 increased SCr (Fig. 1 and Supplementary Table S2). Thus, we focused our further analyses on evaluating the baseline expressions of miR-3168, miR-6125, and miR-4718.

MiR-4718 (p = 0.0388, odds ratio, OR = 1.56, 95% confidence interval, CI: 1.03–2.46) and miR-3168 (p = 0.0456, OR = 1.03, 95% CI: 1.00–1.06) were statistically significantly associated with a higher risk of grade ≥2 increased SCr and miR-6128 showed a trend (p = 0.0681, OR = 1.73, 95% CI: 0.98–3.29) in the univariate regression analyses. In the predictive performance analysis (Fig. 1), miR-4718 showed the highest AUC-ROC (AUC = 0.77, 95% CI: 1.74–5.8) and specificity (79.49), whereas miR-6125 showed the highest sensitivity (88.89).

The baseline expressions of miR-3168, miR-6125, and miR-4718 were highly correlated (SupplementaryTable S3), with a variance inflation factor (VIF) of 3.60, 3.82, and 4.37, respectively, indicating moderate multicollinearity. In multivariate logistic regression analysis (Supplementary Table S4), none of the miRNAs were statistically significantly associated with grade ≥2 increased SCr, probably due to multicollinearity. We combined baseline miRNA expression by summing the expression values of miR-3168, miR-6125, and miR-4718 and performed an univariate logistic regression, adjusting for age and gender. Although the combined miRNAs were statistically significantly associated with a higher risk of ≥2 increased SCr (p = 0.0409, OR = 1.03, 95% CI: 1.00–1.05) with an area under ther curve (AUC) = 0.73 (95% CI: 0.51–0.93), no advantage in prediction performance was observed compared with individually analyzed miRNAs.

Bioinformatics analysis

Differentially expressed miRNAs, with an FR > 2.5 (upregulated), including miR-3168, miR-6125, miR-4718, miR-5694, miR-203a-3p, miR-141-5p, miR-7977, and miR-1303, were selected for bioinformatics analysis. The matrix, showing the top 100 target genes predicted using the upregulated miRNAs is shown in Supplementary Fig. S2. MiR-7977 and miR-203a-3p did not present targets predicted using at least six algorithms and were excluded from further analysis

Differentially expressed miRNAs with an FR < −2.5 (downregulated), including miR-17-5p, miR-1185-1-3p, miR-766-3p, miR-151b, miR-151a-5p, and miR-485-3p, were selected for bioinformatics analysis. The matrix with the top 100 predicted target genes of downregulated miRNAs is shown in Supplementary Fig. S3.

Gene set enrichment analysis was performed with 2017 predicted target genes of upregulated miRNAs and 1677 predicted target genes of downregulated miRNAs. The top 50 canonical signaling pathways of predicted target genes of upregulated and downregulated miRNAs are shown in Fig. 2. For both upregulated and downregulated miRNAs, an enrichment of signaling pathways involved in cisplatin-induced nephrotoxicity was observed, including the nuclear factor kappa B (NF-κB), AMP-activated protein kinase (AMPK), IL-1 signaling pathways, transforming growth factor β (TGF-β), and ErbB/ErbB2-ErbB3/ErbB4 signaling pathways.

Figure 3 highlights the target genes from the miRNA matrices directly involved in cisplatin-induced nephrotoxicity, which can partially explain the nephrotoxicity in patients from the nephrotoxicity group.

We have investigated, for the first time, the role of plasma miRNAs as biomarkers of cisplatin-induced nephrotoxicity using patients’ samples. We used miRNA sequencing to identify miR-3168, miR-6125, and miR-4718 as potential biomarkers of cisplatin-induced nephrotoxicity. The evaluation of miRNAs in a larger cohort of patients provides more evidence of their potential utility. A bioinformatics analysis approach was used to build the post-transcriptional network regulated by the miRNAs and identify genes affected by the differential expression of miRNAs and explain why a subset of cancer patients develop nephrotoxicity during cisplatin treatment.

We performed a miRNA sequencing of plasma samples collected from six patients with and without nephrotoxicity, respectively, five days after cisplatin administration (D5). The most statistically significant differentially expressed miRNAs (miR-3168, miR-6125, and miR-4718) between the two groups of patients were quantified in the plasma samples of a larger cohort of patients. Interestingly, miR-3168, miR-6125, and miR-4718 had a higher expression in baseline samples of patients who developed nephrotoxicity after cisplatin treatment (miR-4718 was statistically significant, whereas miR-3168 and miR-625 showed a trend) and we focused our analysis on their utility as predictors of cisplatin-induced nephrotoxicity.

MiR-4718 was found to be the most promising marker associated with grade ≥2 increased SCr, with an AUC-ROC of 77 and sensitivity of 66.76 and specificity of 79.49 of in predicting cisplatin-induced nephrotoxicity. MiR-3168 was also statistically significant associated with grade ≥2 increased SCr; however, it showed the lowest predictive performance. MiR-6125 did not reach statistical significance, with only a trend being observed, but showed an intermediate predictive performance when compared to miR-4718 and miR-3168 and had the best sensitivity (88.89) in predicting cisplatin-induced nephrotoxicity. MiR-3168, miR-6125, and miR-4718 correlated with each other, and multicollinearity does not suggest an additive effect of the miRNAs in a multivariate model. Thus, the potential of miR-3168, miR-6125, and miR-4718 as predictors of cisplatin-induced nephrotoxicity might rely on their individual effects.

To understand how miR-3168, miR-6125, and miR-4718 are involved in cisplatin-induced nephrotoxicity, we used a bioinformatics approach for the three candidate miRNAs and additional miRNAs with FR > 2.5 or FR < −2.5 identified through sequencing in this study. Target gene enrichment analysis revealed critical signaling pathways involved in cisplatin-induced nephrotoxicity, including the ErbB, NF-κB, TGF-β, IL-1, and AMPK pathways. Upregulated miRNAs targeting genes of the ErbB signaling pathway, such as miR-3168, which targets PDK, could downregulate the pathway, leading to cisplatin-induced apoptosis in renal cells (Fig. 3). MiR-146b was shown to contribute to nephrotoxicity by reducing ErbB4/HER4 expression[15]; although miR-146b was not identified in our study, this corroborates our hypothesis that miRNAs are key regulators of the ErbB signaling pathway during cisplatin-induced nephrotoxicity.

The NF-κB signaling pathway is involved in the induction of miR-375 expression by cisplatin, suppressing the activation of hepatocyte nuclear factor-1β (HNF-1β), which acts as a nephroprotectant and prevents renal damage in tubular cells[16]. Cisplatin has been reported to induce NF-κB phosphorylation and its subsequent translocation to the nucleus, promoting the transcription of inflammatory factors, including TNF-α and IL-1, which contribute to renal damage[17,18]. TNF-α plays a central role in the inflammatory response triggered by cisplatin in the renal cell[19] and contributes to a significant increase in the expression of IL-1β and TGF-β mRNA in rat renal cells after cisplatin treatment[20]. IL-1β is also involved in the recruitment of monocytes in the inflamed renal parenchyma[17].

AMPK regulates autophagy in kidney cells. The inhibition of AMPK was reported to suppress cisplatin-induced autophagy in tubular cells, followed by increased DNA damage and activation of p53[21]. Autophagy contributes toward removing damaged mitochondria and decrease oxidative stress[22]. The inhibition of autophagy in renal cells may also be related to increased oxidative stress, hence, increased nephrotoxicity. Several upregulated miRNAs identified in patients with nephrotoxicity in our study target members of the AMPK signaling pathway, mainly PRKAB2 (by miR-4718 and miR-3168, SupplementaryFig. S2), which could explain nephrotoxicity by repression of autophagy in tubular cells.

This study also aimed to explore where the differentially expressed miRNAs act during the mechanism of action of cisplatin inside the renal cell (Fig. 3). Important genes involved in the mechanism of action of cisplatin are potentially targeted by the miRNAs identified in our sequencing. These differentially expressed miRNAs probably co-operate to modulate different pathways and could be released into the plasma upon renal damage and detected as circulating miRNAs.

The ErbB signaling pathway is an example of the contribution of different miRNAs in regulating genes that encode proteins that interact and lead to cisplatin-induced apoptosis in renal cells. The DNA repair pathway is another example. The reducing activity of the nucleotide-excision repair (NER) pathway could reduce the cell’s ability to repair cisplatin-induced DNA damage, thereby inducing cell death[23], due to the downregulation of ERCC1 (by miR-4718 and miR-6125), TOP2A (by miR-3168), and BRCA1 (by miR-1303).

The downregulation of miR-17-5p and miR-485-3p could explain the increased activity of cisplatin receptor CTR1, which leads to a higher uptake of cisplatin by renal cells, leading to nephrotoxicity[7]. The regulation of genes involved in the mitochondrial apoptosis pathway may also contribute to higher nephrotoxicity, suggested by a decrease in the activity of the anti-apoptotic protein Bcl-2 by miR-3168, miR-6125, and miR-1303. In addition, genes of the cisplatin detoxification pathway, which includes the conjugation of cisplatin with glutathione, are shown to be target of miR-3168, miR-6125 and miR-1303, which could reduce the content of glutathione S-transferase (GST) and reduced glutathione (GSH). Although the role of GST in cisplatin-induced nephrotoxicity is controversial[24], evidence suggests that the decreased GST/GSH leads to a higher concentration of nonconjugated cisplatin inside the renal cell, contributing to higher nephrotoxicity[25].

There are some limitations in the analyses presented. Our study comprised of 60 patients, of which 11 experienced grade ≥2 increased SCr; thus, the sample size was small. Furthermore, given the number of events and concerns of overfitting, we adjusted the analysis only for age and gender. Potential confounders might also be also considered, such as concomitant medication use. We have performed only in silico analysis to hypothesize how miRNAs can modulate cisplatin-induced nephrotoxicity and our hypotheses should be evaluated through functional analysis. Moreover, our findings require independent confirmation in a larger cohort of patients.

In conclusion, this study provided evidence of baseline plasmatic expression of miR-3168, miR-6125, and miR-4718 as potential predictors of cisplatin-induced nephrotoxicity. The expression of miR-3168, miR-6125, and miR-4718 could be assessed before treatment to improve the use of cisplatin in the clinical setting and selecting patients to be treated with cisplatin with an improved risk/benefit ratio.

Patient selection, treatment regimen, and toxicity

This nested case-control study included patients with primary squamous cell carcinoma of the head and neck who received high doses of cisplatin (80 or 100 mg/m²) concomitantly with conventional radiotherapy.

Blood samples were collected before and on D5 for SCr and BUN measurement. SCr clearance was estimated using the Cockcroft–Gault formula[26]. Nephrotoxicity was classified as grade ≥2 increased SCr according to Common Toxicity Criteria for Adverse Events (CTCAE) version 4. Grade ≥2 increased SCr is classified as more than two times SCr baseline values.

Sequencing of miRNAs

MiRNAs were extracted from samples collected on D5 from six patients with grade ≥2 increased SCr (nephrotoxicity group) and six patients with grade = 0 increased SCr (non-nephrotoxicity group) using the miRNeasy Serum/Plasma Kit (Qiagen, cat no. 217184). Library preparation was performed using the QIAseq^TM miRNA Library Kit (Qiagen, cat no. 331502). The quality control of the libraries was performed by analyzing their size using the 4200 TapeStation Instrument (Agilent Technologies) and the concentration by using Qubit Fluorometric Quantitation (Thermo Fisher Scientific). The molarity of each library was calculated and diluted to 4 nM for sequencing.

Sequencing was performed using the MiSeq Reagent Kit v3, 150-cycle (Illumina, MS-102-3001). Sequencing data analysis was performed using the GeneGlobe Data Analysis Center (Qiagen). DESeq2 method was used for normalization. GeneGlobe provides the fold-change (FC = miRNA expression of the nephrotoxicity group / miRNA expression of the non-nephrotoxicity group), FR (FR = FC, if FC ≥1 or FR = FC ≥1, if FC < 1) and p-value based on the Wald test for each miRNA. GeneGlobe also provides comments for different miRNAs. These comments can be “A”, if the miRNA's average expression level is relatively low (<10) in either the nephrotoxicity or the non-nephrotoxicity group and is reasonably high in the other group (>10); “B”, if the miRNA's average expression level is relatively low (<10) in both the nephrotoxicity and non-nephrotoxicity groups; or “C”, if the miRNA has an expression count of zero in both the nephrotoxicity and non-nephrotoxicity groups, meaning that its expression was undetected, making the FC result erroneous and uninterpretable.

Validation of selected miRNAs

The three most statistically significant differentially expressed miRNAs, with the highest FR identified through sequencing, were selected for validation in a larger cohort of patients. After extracting the miRNAs from all plasma samples collected on the baseline and on D5 from all patients included in the study, cDNA synthesis was performed using the TaqMan™ Advanced miRNA cDNA Synthesis Kit (Applied Biosystems, Cat No./ID: A28007) and qPCR using TaqMan™ Advanced miRNA Assays (Applied Biosystems, Cat No./ID:A25576). In addition, qPCR of the exogenous control cel-miR-39 and the endogenous control hsa-miR-16 for normalization was performed. Hsa-miR-16 was selected as an endogenous control because its expression was shown to be stable in plasma samples[27]. We normalized our analyses using cel-miR-39 due to its smaller variations of cycle threshold (CT) than has-miR-16. Patients with cel-miR-39 expression above two SDs were excluded from the analysis.

MiRNA qPCR results were analyzed using the QuantStudioTM Real Time PCR Software 6. Each miRNA had its expression evaluated and relative expressions were obtained using the 2^−∆∆CT method[28], where ∆CT = candidate CTmiRNA − CT cel-miR-39 and ∆∆CT = ∆CT – mean ∆CTs of patients with grade <2 increased SCr.

Statistical analysis

Statistical analysis was performed using the R software environment for statistical computing[29]. MiRNA expression was compared between the nephrotoxicity and non-nephrotoxicity groups at baseline and D5 using the Mann–Whitney U tests. For comparisons with a p-value of <0.10, we tested the association between miRNA expression and grade ≥2 increased SCr by calculating the OR from a univariate logistic regression analysis, considering grade ≥2 increased SCr as a binary outcome, adjusted for age and gender. Receiver operating characteristic curves (ROC) were used to estimate the AUC, optimal cutpoint, sensitivity, and specificity.

Spearman correlations and multicollinearity using VIF were calculated for miRNA expression. Multivariate logistic regression analyses were performed to evaluate the relationship between grade ≥2 increased SCr and miRNAs expression at baseline. Age and gender were included as covariates in the multivariate model.

Bioinformatics analysis

Differentially expressed miRNAs with FR > 2.5 or FR < −2.5 were selected for bioinformatics analysis. All miRNAs selected presented no comments or comment “A” and a p-value of <0.05. The identification of predicted miRNA target genes was performed using the miRWalk platform[30], which provides a list of predicted miRNA target genes according to twelve different algorithms, including TargetScan[31]. Only target genes predicted by TargetScan and at least five different databases were selected for the following analyses.

Two matrices were constructed to identify the interaction between miRNAs and their predicted target genes: one for upregulated miRNAs and one for downregulated miRNAs. The matrices were sorted according to the potential target genes of different miRNAs (stronger evidence); genes targeted by at least two different miRNAs were selected for unsupervised enrichment analysis using the Ingenuity Pathway Analysis (IPA^®, Qiagen) software to identify the main canonical signaling pathways involving differentially expressed miRNAs.

The KeggMapper[32] was used to visualize the canonical signaling pathways involving miRNA target genes. We selected the signaling pathway for platinum-resistance and adapted it to the signaling pathway of cisplatin-induced nephrotoxicity in renal cells. Thus, we could visualize how miRNAs regulate genes in cisplatin-induced nephrotoxicity.

Ethics approval and consent to participate

All participants provided written informed consent for sample collection and analysis. The study was performed in accordance with the Declaration of Helsinki and was approved by the local IRB (The Research Ethics Committee of the School of Medical Sciences/UNICAMP).

Data availability

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

The authors would like to acknowledge the São Paulo Research Foundation-FAPESP (2017/02338-0 and 2017/11329-4) and the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brazil (CAPES) (Finance Code 001) for the funding and Enago (www.enago.com) for the English language review.

Author contributions

Conceptualization: JCFQ and PM. Methodology: JCFQ and PM. Software: JCFQ and MVG. Formal Analysis: JCFQ. Investigation: JCFQ, MAC, NGT, LBB, JMO, TSC. Resources: JBB, EP, MHH, CSPL, and PM. Data Curation: JCFQ. Writing – Original Draft Preparation: JCFQ. Writing – Review & Editing: MVG and PM. Supervision: MHH, CSPL, and PM. Funding Acquisition: JCFQ, MAC, and PM.

Competing interests

The author declare no competing interests

Funding

This work was supported by the São Paulo Research Foundation-FAPESP (2017/02338-0 and 2017/11329-4) and by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brazil (CAPES) (Finance Code 001).

[1] Gold, J.M. & Raja, A. Cisplatin. In: StatPearls [Internet]. Treasure Island (FL): StatPearls Publishing. https://www.ncbi.nlm.nih.gov/books/NBK547695/ (2020).

[2] Siddik, Z.H. Cisplatin: Mode of cytotoxic action and molecular basis of resistance. Oncogene22, 7265-7279 (2003).

[3] Quintanilha, J.C.F. et al. Cisplatin-induced human peripheral blood mononuclear cells’ oxidative stress and nephrotoxicity in head and neck cancer patients: the influence of hydrogen peroxide. Mol. Cell. Biochem. 440, 139-145 (2018).

[4] Visacri, M.B. et al. Adverse drug reactions and kinetics of cisplatin excretion in urine of patients undergoing cisplatin chemotherapy and radiotherapy for head and neck cancer: A prospective study. Daru25, 12 (2017).

[5] Miller, R.P., Tadagavadi, R.K., Ramesh, G. & Reeves, W.B. Mechanisms of cisplatin nephrotoxicity. Toxins (Basel) 2, 2490-2518 (2010).

[6] Ciarimboli, G. Membrane transporters as mediators of cisplatin side-effects. Anticancer Res.34, 547-550 (2014).

[7] Pabla, N., Murphy, R.F., Liu, K. & Dong, Z. The copper transporter Ctr1 contributes to cisplatin uptake by renal tubular cells during cisplatin nephrotoxicity. Am. J. Physiol. Renal Physiol.296, F505-F511 (2009).

[8] Eljack, N.D. et al. Mechanisms of cell uptake and toxicity of the anticancer drug cisplatin. Metallomics 6, 2126-2133 (2014).

[9] Herrera-Pérez, Z., Gretz, N. & Dweep, H. A Comprehensive Review on the Genetic Regulation of Cisplatin-induced Nephrotoxicity. Curr. Genomics17, 279-293 (2016).

[10] Gautier, J.C. et al. Evaluation of novel biomarkers of nephrotoxicity in two strains of rat treated with cisplatin. Toxicol. Pathol.38, 943-956 (2010).

[11] Quintanilha, J.C.F., Saavedra, K.F., Visacri, M.B., Moriel, P. & Salazar, L.A. Role of epigenetic mechanisms in cisplatin-induced toxicity. Crit. Rev. Oncol. Hematol.137, 131-142 (2019).

[12] Lagos-Quintana, M., Rauhut, R., Lendeckel, W. & Tuschl, T. Identification of novel genes coding for small expressed RNAs. Science294, 853-858 (2001).

[13] Pavkovic, M. et al. Detection of drug-induced acute kidney injury in humans using urinary KIM-1, miR-21, -200c, and -423. Toxicol. Sci. 152, 205-213 (2016).

[14] Weber, J.A. et al. The microRNA spectrum in 12 body fluids. Clin. Chem. 56, 1733-1741 (2010).

[15] Zhu, Y. et al. MicroRNA-146b, a Sensitive Indicator of Mesenchymal Stem Cell Repair of Acute Renal Injury. Stem Cells Transl. Med. 5, 1406-1415 (2016).

[16] Hao, J. et al. MicroRNA-375 is induced in cisplatin nephrotoxicity to repress hepatocyte nuclear factor 1-β. J. Biol. Chem. 292, 4571-4582 (2017).

[17] Volarevic, V. et al. Molecular mechanisms of cisplatin-induced nephrotoxicity: A balance on the knife edge between renoprotection and tumor toxicity. J. Biomed. Sci. 26, 25 (2019).

[18] Ramesh, G., Zhang, B., Uematsu, S., Akira, S. & Reeves, W.B. Endotoxin and cisplatin synergistically induce renal dysfunction and cytokine production in mice. Am. J. Physiol. Renal Physiol.239, F325-F232 (2007).

[19] Ramesh, G. & Brian Reeves, W. TNF-α mediates chemokine and cytokine expression and renal injury in cisplatin nephrotoxicity. J. Clin. Invest. 110, 835-842 (2002).

[20] Kumar, P., Sulakhiya, K., Barua, C.C. & Mundhe, N. TNF-α, IL-6 and IL-10 expressions, responsible for disparity in action of curcumin against cisplatin-induced nephrotoxicity in rats. Mol. Cell. Biochem. 431, 113-122 (2017).

[21] Wei, L. et al. AMP-activated protein kinase regulates autophagic protection against cisplatin-induced tissue injury in the kidney. Genet. Mol .Res.14, 12006-12015 (2015).

[22] White, E. Deconvoluting the context-dependent role for autophagy in cancer. Nat. Rev. Cancer12, 401-410 (2012).

[23] Kelland, L. The resurgence of platinum-based cancer chemotherapy. Nat. Rev. Cancer 7, 573-584 (2007).

[24] Townsend, D.M., Tew, K.D., He, L., King, J.B. & Hanigan, M.H. Role of glutathione S-transferase Pi in cisplatin-induced nephrotoxicity. Biomed. Pharmacother.63, 79-85 (2009).

[25] Sadzuka, Y., Shimizu, Y. & Takino, Y. Role of glutathione S-transferase isoenzymes in cisplatin-induced nephrotoxicity in the rat. Toxicol. Lett. 70, 221-222 (1994).

[26] Cockcroft, D.W. & Gault, M.H. Prediction of creatinine clearance from serum creatinine. Nephron 16, 31-41 (1976).

[27] Poel, D., Buffart, T.E., Oosterling-Jansen, J., Verheul, H.M.W. & Voortman, J. Evaluation of several methodological challenges in circulating miRNA qPCR studies in patients with head and neck cancer. Exp. Mol. Med.50, e454 (2018).

[28] Livak, K.J. & Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2^-ΔΔCT method. Methods 25, 402-408 (2001).

[29] RCore Team. R: A language and environment for statistical computing.Vienna, Austria: R Foundation for Statistical Computing (2020).

[30] Dweep, H. & Gretz, N. MiRWalk2.0: A comprehensive atlas of microRNA-target interactions. Nat. Methods12, 697 (2015).

[31] Agarwal, V., Bell, G.W., Nam, J.W. & Bartel, D.P. Predicting effective microRNA target sites in mammalian mRNAs. eLife4, e05005 (2015).

[32] Kanehisa, M. & Sato, Y. KEGG Mapper for inferring cellular functions from protein sequences. Protein Sci.29, 28-35 (2020).

Table 1. Patients with head and neck cancer treated with cisplatin. KPS Karnofsky performance status, SD standard deviation.

Patient and clinical characteristics
		Patients whose miRNA samples were sequenced		All patients (n=60)
		Non-nephrotoxicity group (n=6)	Nephrotoxicity group (n=6)	All patients (n=60)
Age (mean ± SD, years)		60.3 ± 4.9	57.5 ± 7.4	58.4 ± 7.3
Sex (n, %)
Men		5 (83.3)	4 (66.7)	54 (90.0)
Women		1 (16.7)	2 (33.3)	6 (10.0)
Ethnicity (n, %)
Caucasian	4 (66.7)		6 (100.0)	45 (75.0)
Non-Caucasian	2 (33.3)		0 (0.0)	15 (25.0)
Smoking category (n, %)
Never smoked		1 (16.7)	1 (16.7)	6 (10.0)
Light smoker		0 (0.0)	0 (0.0)	3 (5.0)
Moderate smoker		0 (0.0)	1 (16.7)	4 (6.7)
Heavy smoker		5 (100.0)	4 (66.6)	47 (78.3)
Drinking category (n, %)
Abstainer		2 (33.3)	2 (33.3)	9 (15.0)
Light drinker		0 (0.0)	0 (0.0)	5 (8.3)
Moderate drinker		0 (0.0)	0 (0.0)	2 (3.3)
Heavy drinker		3	1 (16.7)	16 (26.7)
Very heavy drinker		1 (16.7)	3 (50.0)	28 (46.7)
KPS (n, %)
60		0 (0.0)	0 (0.0)	1 (1.7)
70		0 (0.0)	0 (0.0)	1 (1.7)
80		2 (33.3)	3 (50.0)	9 (15.0)
90		4 (66.7)	3 (50.0)	40 (66.6)
100		0 (0.0)	0 (0.0)	9 (15.0)
Tumor site (n, %)
Oropharynx		0 (0.0)	1 (16.7)	15 (25.0)
Hypopharynx		2 (33.3)	3 (50.0)	9 (15.0)
Larynx		1 (16.7)	1 (16.7)	13 (21.7)
Oral cavity		3 (50.0)	1 (16.7)	20 (33.3)
Hypopharynx + Larynx		0 (0.0)	0 (0.0)	1 (1.7)
na		0 (0.0)	0 (0.0)	2 (3.3)
Stage (n, %)
I		0 (0.0)	0 (0.0)	0 (0.0)
II		0 (0.0)	0 (0.0)	6 (10.0)
III		1 (16.7)	0 (0.0)	7 (11.6)
IV		5 (83.3)	6 (100.0)	46 (76.7)
na		0 (0.0)	0 (0.0)	1 (1.7)
Nephrotoxicity parameters
Serum creatinine, SCr (mg/dL)
Baseline		0.94 ± 0.18	0.75 ± 0.10	0.82 ± 0.20
D5		0.94 ± 0.21	3.39 ± 2.13	1.36 ± 1.01
Creatinine clearance (mL/min)
Baseline		70.20 ± 22.93	97.62 ± 40.55	87.70 ± 26.90
D5		69.72 ± 24.01	25.05 ± 11.46	61.86 ± 25.54
Blood urea nitrogen, BUN (mg/dL)
Baseline		38.50 ± 14.28	21.33 ± 9.48	28.75 ± 10.95
D5		45.67 ± 14.88	92.00 ± 40.61	53.63 ± 22.99

SupplementaryFile.pdf

Download PDF

Version 1

posted

You are reading this older preprint version

Read the latest preprint version →

MiR-3168, miR-6125, and miR-4718 as Potential Predictors of Cisplatin-Induced Nephrotoxicity in Patients with Head and Neck Cancer

Archived Versions:

Version 1

Abstract

Figures

Introduction

Results

Discussion

Methods

Declarations

References

Tables

Supplementary Files

Archived Versions:

Version 1