Next Article in Journal
Effects of Aging and Methionine Restriction on Rat Kidney Metabolome
Next Article in Special Issue
ADGRL4/ELTD1 Silencing in Endothelial Cells Induces ACLY and SLC25A1 and Alters the Cellular Metabolic Profile
Previous Article in Journal / Special Issue
Historical Biobanks in Breast Cancer Metabolomics— Challenges and Opportunities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Organ NMR Metabolomics to Assess In Vivo Overall Metabolic Impact of Cisplatin in Mice

by
Tatiana J. Carneiro
1,
Rita Araújo
1,
Martin Vojtek
2,
Salomé Gonçalves-Monteiro
2,
Carmen Diniz
2,
Ana L.M. Batista de Carvalho
3,
Maria Paula M. Marques
3,4 and
Ana M. Gil
1,*
1
Department of Chemistry and CICECO–Aveiro Institute of Materials, University of Aveiro, 3810-193 Aveiro, Portugal
2
LAQV/REQUIMTE, Department of Drug Sciences, Laboratory of Pharmacology, Faculty of Pharmacy, University of Porto, 4150-755 Porto, Portugal
3
“Química-Física Molecular”, University of Coimbra, 3004-535 Coimbra, Portugal
4
Department of Life Sciences, Faculty of Science and Technology, University of Coimbra, 3000-456 Coimbra, Portugal
*
Author to whom correspondence should be addressed.
Metabolites 2019, 9(11), 279; https://doi.org/10.3390/metabo9110279
Submission received: 24 October 2019 / Revised: 8 November 2019 / Accepted: 11 November 2019 / Published: 13 November 2019
(This article belongs to the Special Issue Metabolomics in the Study of Disease)

Abstract

:
This work describes, to our knowledge, the first NMR metabolomics analysis of mice kidney, liver, and breast tissue in response to cisplatin exposure, in search of early metabolic signatures of cisplatin biotoxicity. Balb/c mice were exposed to a single 3.5 mg/kg dose of cisplatin and then euthanized; organs (kidney, liver, breast tissue) were collected at 1, 12, and 48 h. Polar tissue extracts were analyzed by NMR spectroscopy, and the resulting spectra were studied by multivariate and univariate analyses. The results enabled the identification of the most significant deviant metabolite levels at each time point, and for each tissue type, and showed that the largest metabolic impact occurs for kidney, as early as 1 h post-injection. Kidney tissue showed a marked depletion in several amino acids, comprised in an overall 13-metabolites signature. The highest number of changes in all tissues was noted at 12 h, although many of those recovered to control levels at 48 h, with the exception of some persistently deviant tissue-specific metabolites, thus enabling the identification of relatively longer-term effects of cDDP. This work reports, for the first time, early (1–48 h) concomitant effects of cDDP in kidney, liver, and breast tissue metabolism, thus contributing to the understanding of multi-organ cDDP biotoxicity.

1. Introduction

Cisplatin (cDDP) is a platinum(II)-containing drug (cis-[Pt(NH3)2Cl2]) which is well known for its anticancer properties [1,2]. The impact of cDDP in cell division was observed for the first time in 1965 (in E. coli) and the compound was approved as an antineoplastic agent in 1978 [2]. CDDP is presently used, often in combination with other drugs, for treatment of many different types of cancer, e.g., head and neck, lung, testicular, ovarian, and bladder cancer [1,3]. However, its clinical use is limited to some extent due to frequently acquired drug resistance and severe cytotoxic effects such as nephrotoxicity, hepatotoxicity, cardiotoxicity, neurotoxicity, and bone marrow toxicity [3]. Attempts to overcome these drawbacks have been made through the development of several other Pt (or other metal)-based drugs, such as carboplatin (a Pt complex with bidentate cyclobutane dicarboxylic acid replacing the chloride ions), which represents one of the most successful [1]. Furthermore, the combination of cDDP with other chemotherapeutic agents, such as taxanes (paclitaxel), anthracyclines (doxorubicin), or antimetabolites (tegafur-uracil, gemcitabine, capecitabine), has been extensively used to reduce cDDP deleterious side effects and acquired resistance, while increasing patient compliance to treatment [1,4,5].
CDDP biotoxicity is known to lead to considerable deviant effects at the intracellular level (namely DNA damage, organelle dysfunction, oxidative stress, or inflammation), as well as on the plasma membrane, which then translates into deviations in endogenous metabolic status [1,2]. This abnormal metabolic behavior may potentially unveil new early biomarkers of cDDP biotoxicity, which may enable timely recognition and management in a clinical context; metabolomics is an enticing means to identify such biomarkers. Indeed, metabolomics has emerged as a promising approach to evaluate in vitro and in vivo responses to drugs [6,7]. This strategy provides valuable insight into the metabolic perturbations produced locally (in affected tissues) and at the systemic level, as viewed through the metabolome of biofluids [8,9]. In the context of response to cDDP, most studies have addressed its role in cancer treatment, with fewer reports dedicated to cDDP toxicity alone (i.e., in healthy cell lines or animal models). In cancer treatment with cDDP, metabolomics has been applied to several human cancer cell lines in search of biomarkers of chemotherapy efficacy (lung [10,11], breast [12,13], bone [14,15], brain [16], and ovarian [17,18]) or resistance (ovarian [17,18,19], ovarian and lung [20], breast [21]). CDDP treatment efficacy and nephrotoxicity have also been assessed through kidney tissue metabolomics in a mouse model of metastatic melanoma [22], and biofluids, namely serum/plasma [23,24,25] and urine [26], in a model of lung cancer, both for cDDP sole [24,25] and combined administration [23,24].
Regarding the use of metabolomics to improve the understanding of cDDP biotoxicity mechanisms, some studies have addressed healthy liver [27], lung [10], and kidney [28] cell lines. In particular, human liver L02 cells incubated with cDDP and analyzed (upon extraction) by NMR metabolomics [27] showed drug effects on cell metabolism for variable orders of magnitude of drug concentrations (1 nM to 1 mM), noting that the highest dose induced pronounced perturbations in the levels of amino acids and nucleotides. NMR metabolomics was also used to analyze the (polar) metabolic profiles of human healthy lung cells (MRC5 cell line), compared to cancer cells (A549 cell line), upon treatment with several chemotherapeutics agents, including cDDP [10]. The heathy cell line showed decreasing trends in several metabolites (e.g., branched-chain amino acids (BCAAs), glutamine and glycine, lactate, glutathione, and taurine) as an indication of drug toxicity effects, possibly underlying the response of cancer cells to cDDP treatment. In the case of renal cell lines, a study of the human renal proximal tubular RPTEC/TERT1 cell line was carried out through ‘multi-omics’ approaches, namely transcriptomics, proteomics, and metabolomics [28]. Both cellular supernatants and lysates analyses revealed distinct perturbations, particularly in relation to choline, fatty acids, and glutathione metabolism. Furthermore, metabolomics has also focused on healthy rodent models to address tissue damage upon cDDP exposure [29] and to help assess nephrotoxicity, since this is the major dose-limiting factor related to cDDP treatment [5]. Nephrotoxicity-related metabolic changes have been observed in the biofluids (blood plasma/serum and urine) of rodent models, reflecting clear dose- and time-dependent deviant metabolic responses to cisplatin exposure [30,31,32,33,34,35,36,37,38,39,40,41]. Although biofluids metabolic profiles carry valuable information on the systemic metabolic response to cDDP, and may be an important source of non-invasive early markers of drug toxicity, tissue metabolic analyses (including but not restricted to kidney) may provide important information on local deviant metabolism and interorgan metabolic interplays. To our knowledge, only cisplatin-exposed kidney tissue has been studied through metabolomics, notably using MS-based untargeted metabolomics [3,42,43]. The first studies addressed Sprague-Dawley (SD) rats, with an early report monitoring both urine and kidney tissue after 1, 5, and 28 days post cDDP injection (0.5 mg/kg) [42]. An interesting result related to the early observation of altered levels of polyamines and amino acids (namely BCCAs) in urine, whereas tissue revealed depleted levels of several amino acids (consistently with their increased excretion) and nucleosides. The metabolic heterogeneity in kidney was addressed in two further metabolomic studies, either considering medulla and cortex measurements separately (using 2.5, 5, and 10 mg/kg cDDP) [43], or by using MS imaging (MSI) techniques in tandem with LC-MS (using 30 mg/kg cDDP) [3]. Upon 7 days post-injection with different cisplatin doses [43], metabolic changes were shown to be more significant and dose-dependent in kidney medulla, compared to cortex. Many amino acids were found to be down-regulated in both kidney sections, in tandem with up-regulation of fatty acids (FA), diacylglycerols (DG), ceramide, and sphingosine, along with a set of metabolites noted to distinguish the metabolic profiles of medulla and cortex. On a shorter time-scale (3–72 h), in a study of mice exposed to a large cDDP dose (30 mg/kg b.w.) [3], significant region-specific metabolic differences were confirmed, and early changes (at times < 24 h) were noted in several compounds, such as amino acids (up-regulated, in apparent contradiction with previous observations), nucleotides, and intermediates of folate metabolism.
In this work, a 1H NMR-based metabolomics strategy was carried out, for the first time to our knowledge, to assess the impact of cisplatin on the metabolism of the kidney, adding the novel study of the concomitant metabolic changes in the liver and breast tissues of the same healthy mice. The latter type of tissue should be particularly useful as a control matrix for studies of breast cancer animal models treated with cisplatin or other drugs. As a technique which is complementary to MS-based methods, NMR metabolomics should confirm previous observations in kidney tissue (for comparable timescales and drug doses), while unveiling novel metabolic changes, not only in kidney, but also in liver and breast tissue. Focusing on a very short timescale of exposure to cDDP (1 to 48 h), this study provides an early dynamic metabolic picture of the interplay of different organs in responding to the drug, thus unveiling potential early multi-organ metabolic markers of cDDP biotoxicity.

2. Results

2.1. Typical 1H NMR Spectra of Aqueous Extracts of Mice Kidney, Liver and Breast Tissue

Figure 1 shows the average 1H NMR spectra of polar extracts obtained from the kidney, liver, and breast tissue of healthy mice (controls), including semi-quantitative indication of predominant metabolites (compounds noted in red). Table S1 presents the overall list of assigned peaks (52 in total) and corresponding information regarding metabolite predominance in each type of tissue. In general, our results are consistent with previous reports of 1H NMR spectra of kidney [44,45,46] and liver extracts [46,47,48,49], adding new information regarding the polar metabolome of mice kidney (namely, changes in the levels of inosine monophosphate (IMP), uridine monophosphate (UMP), allantoin, threonine, 3-hydroxyisobutyrate (3-HIBA), trimethylamine N-oxide (TMAO), uridine diphosphate-glucuronate (UDP-GlcA)), liver (asparagine, glycerophosphocholine (GPC), hypoxanthine, hippurate, IMP, threonine, 3-HIBA, TMAO, and UDP-GlcA), and breast tissue, for which, to our knowledge, no previous 1H NMR metabolomics study has been reported. By visual inspection of the spectra in Figure 1, it becomes clear that alanine, lactate, and taurine predominate in all three tissue types, whereas specific characteristics are noted, as expected, for different tissues: (1) kidney comprises relatively high levels of betaine, m-inositol, and inosine, (2) liver shows high levels of acetone, glucose, glycogen, and glutathione (reduced form, GSH), and (3) breast tissue shows high levels of creatine, phosphocholine, and glycerol moieties arising from glycerolipids. The following section describes the changes in the complex metabolic profiles of each tissue type upon time-course exposure to cDDP.

2.2. Impact of Mice Exposure to Cisplatin

2.1.1. Kidney Metabolic Profiling

Unsupervised analysis with PCA of the spectra of both controls and cDDP-exposed groups showed group separation only in the case of kidney (Figure 2, left), reflecting a stronger impact of the drug on the metabolic profile of that organ compared to liver and breast tissue. The corresponding PLS-DA score plots (Figure 2, right) confirmed such separation, given the satisfactory predictive power observed (Q2 = 0.48, compared to Q2 = 0.38 and 0.30 for liver and breast tissue, respectively). Particularly for kidney samples, it is interesting to note, through examination of PLS-DA plots, that both control and exposed kidney samples seemed to follow a slight time course tendency towards positive LV2, with relatively larger sample dispersion for the drug-exposed group.
In order to systematically identify all the metabolites in kidney varying with exposure, at some stage of exposure, pairwise multivariate analysis was carried out for each timepoint (see results for t = 12 h in Figure 3). In spite of the low sample numbers, the PLS-DA model obtained for kidney tissue at 12 h post-injection showed high robustness (predictive power Q2 = 0.73), together with sample groups for 1 h (Q2 = 0.81, not shown), both showing better separation than for 48 h (Q2 = 0.54, not shown), which seems to indicate some degree of recovery from the shorter-term impact of cDDP. The underlying varying metabolites were identified through inspection of the LV1 loading plots corresponding to each time (e.g., for 12 h, as shown in Figure 3a, right), and their variation quantified by integration (Table 1). A set of 14 metabolites (including unassigned spin system U6) described the short-term impact of cDDP on kidney tissue, followed by an enlarged signature of 21 compounds at 12 h post-injection, and by a lesser number of changes (11-metabolite signature) at 48 h. The heatmap of these variations, organized by compound families (Figure 4), clearly shows that kidney is the most affected organ in terms of metabolism, compared to liver and breast tissue. Hence, its close examination will facilitate the identification of the short- and longer-term effects of cDDP (in the study range of 1–48 h), as well as metabolic recoveries and persistent deviations. Also, as this is, to our knowledge, the first NMR metabolomics study of cDDP exposure in mouse kidney, the present results provide additional information to existing MS-based reports on healthy animals [3,42,43], namely regarding cDDP-induced changes in 3-HIBA (↑↓), allantoin (↓, 12 h), adenosine monophosphate (AMP) (↑), betaine (↑,12 h), choline (↓), fumarate (↓, 48 h), m-inositol (↑), phosphocholine (PC) (↑, 12 h), taurine (↑, 1 h), trimethylamine (TMA) (↓, 48 h), and UDP-GlcA (↑, 12 h) (see Table 1).
The shorter-term (1 h) effects of cDDP on kidney comprise a general decrease in amino acids, with the exception of taurine (which increases significantly), in tandem with decreased free choline, and increased levels of ADP/AMP/UMP, 3-HIBA, and m-inositol (see Table 1, Figure 4a). Medium-term effects (12 h) involve a wider range of compounds, including additional decreases in asparagine, creatine, allantoin, 3-HIBA (inversed variation in relation to 1 h), and many unassigned resonances, as well as increased levels of PC, glucose, nicotinamide adenine dinucleotide (NAD+), betaine, and UDP-GlcA. The column corresponding to 48 h shows that many metabolites approach levels characteristic of controls, as illustrated by white or light colors without statistical relevance (Figure 4). More persistent amino acid variations relate to decreased alanine, tyrosine, and the three BCCAs. Other longer-term deviant levels comprise decreased levels of hypoxanthine, fumarate, niacinamide (NAM), TMA, unassigned resonance U1 (δ 0.89, triplet), and increased m-inositol.

2.1.2. Liver Metabolic Profiling

In the case of liver, this is, to the best of our knowledge, the first metabolomics study of liver tissue exposed to cDDP using either NMR or MS. In liver, PCA and PLS-DA scores show either no group separation (PCA, see Figure 2b, left) or weak separation in supervised analysis (see Figure 2b, right). In PLS-DA, a slight time-course tendency is hinted at in controls, but no clear time tendency is seen in drug-exposed liver. Interestingly, this group shows the least dispersion in PLS-DA scores among all three treated tissues (see Figure 2, right), thus reflecting the higher tissue homogeneity of liver, and suggesting that liver exposure to cDDP seems to be the least time-dependent. At 12 h post-injection, a separation is clear in both PCA and PLS-DA scores (Q2 = 0.70) (Figure 3b), and the profile illustrated in the corresponding loadings plot (Figure 3b, right) reflects changes in some of the specific predominant metabolites in liver, namely lactate, GSH, and ketone bodies.
The current results show that the short-term effects of cisplatin in liver are limited to decreases in lactate and threonine (as previously reported for human liver cell line L02 exposed to cDDP [27]), to which we add three decreased unassigned resonances (U7 δ 0.73 singlet, U10, δ 4.31 doublet, and U11, δ 6.18 singlet), with U7 tentatively assigned to bile acids (Table 2). The stronger effects of cDDP on the liver metabolic profile take place at 12 h with a 15-metabolite signature, followed by recovered levels for most of those metabolites at 48 h and longer-term novel variations in 7 metabolites (in addition to reversed changes in unassigned resonances U10 to U12) (see Table 2). Inspection of the corresponding heatmap (Figure 4b) shows that visible metabolite variations tend to generally comprise decreases at 1 and 12 h (note GSH and acetone marked decreases at 12 h, in particular), with few exceptions of increasing tendencies (note increased histidine and lactate at 12 h). Conversely, at 48 h, there is a clear predominance of increased metabolite levels in all families, except organic acids (namely, note the significant acetate decrease at 48 h). Hence, the dynamics of cDDP exposure seems to be slightly delayed in liver compared to kidney, with most changes occurring at 12 h post-injection and beyond. Notably, the most marked persistent changes at 48 h comprise increased levels of tyrosine, ADP/AMP/ATP/IMP/NAD+, and four unassigned resonances, in tandem with an acetate decrease.

2.1.3. Breast Tissue Metabolic Profiling

Breast tissue is shown to exhibit the least affected metabolic profile following treatment with cisplatin (Figure 2c), with no clear time course trend, which may also reflect the higher heterogeneity of this tissue. However, it is worthy of note that a slightly higher number of metabolite changes take place at 12 h post-injection, compared to 1 and 48 h (Table 3). The first indication of deviant metabolism is increased levels of the ketone body 3-hydroxibutyrate (3-HBA), which, however, quickly levels out back to control levels at 12 h, when decreased inosine and increased ADP levels are also noted (in tandem with an increase in U13, δ 1.25, singlet). All these changes are, however, restored to control values at 48 h, and only glutamine and U14 (δ 7.68, singlet) remain at lower levels, compared to controls. To our knowledge, this is the only study to date to have reported on drug-elicited, metabolic changes in breast tissue. As this type of tissue is mainly composed of lipids, the additional analysis of the lipophilic metabolome would certainly be of interest. However, the polar metabolome changes reported here already advance the early importance of 3-HBA, followed by a nucleotide deviant response, and finally, an important role played by glutamine in the later stages (see Figure 4c).

3. Discussion

Regarding the effect of cDDP on mouse organ metabolism, this work is, to the best of our knowledge, the first NMR-based metabolomics study. In the case of kidney metabolism, as a complementary strategy to MS-based metabolomics (in terms of higher holistic nature and reproducibility, simpler non-destructive sample handling, albeit lower sensitivity), this NMR work complements previous MS studies of mice kidney extracts carried out to characterize cDDP effects [3,42,43]. Furthermore, the present work considers effects as early as 1 h post-injection for all three tissues types, a time range, to the best of our knowledge, only comparable to that of one previous MS-based report on kidney exposed to cDDP ([3], 3–48 h).
In terms of the relative dynamics of the metabolic effects, it was observed that kidney suffers a large impact of the drug at only 1 h, evidencing the strongest short-term effect within the different tissues considered. At 1 h, the general amino acids decrease noted in kidney tissue is consistent with previous reports that employed cDDP doses of similar orders of magnitude (0.5 to 10 mg/mL, [42,43]), although tissue metabolic changes were measured within the timescale of days. One such study showed that a tandem increase in amino acids occurs in urine; this was interpreted as evidence of reabsorption impairment of amino acids in kidney due to nephrotoxicity and, thus, inducing increased amino acid excretion. However, the authors [43] also note some exceptions to the tissue decreasing levels tendency, e.g., for glutamine, similar to our observation regarding taurine. On the other hand, another recent report using a combination of MS imaging and LC-MS [3] indicated a general amino acid increase upon drug administration at times as short as 3 h, again with some exceptions (e.g., decreases in glutamate, glutamine, and cysteine). The latter study, however, involved the administration of a 10-fold higher cDDP concentration (30 mg/kg), which suggests that a strong dose-dependent effect on amino acid levels occurs, from decreased levels, at doses of a few mg/kg, to increased general levels, at much higher cDDP doses. Incidentally, U-shaped variations in metabolite levels have also been reported for human liver cell line, subjected to cDDP in the 1 nM–1 mM range [27], suggesting that similar behavior may be observed in vivo (namely in kidney). It is interesting to note, however, that in spite of the two previous, apparently contradicting studies in relation to amino acids levels in kidney tissue [3,43], glutamate levels behave similarly, in that they oppose the trend of most amino acids. In the present study, we noted taurine as differing from the decreased overall trend, which is consistent with a biochemical proximity between taurine and glutamate, as suggested before [3] for glutamate/hypotaurine/cysteine, all of which are involved in taurine and hypotaurine metabolism. In addition, cysteine (and methionine) metabolism is associated with glutathione metabolism, where GSH plays an important role in anti-oxidative mechanisms. Other early effects of cDDP on kidney comprise increased m-inositol, adenine, and uridine nucleotides, namely ADP/AMP/UMP, the latter in apparent contradiction to the recent high-dose cDDP report [3], again evidencing an apparent reversal effect induced by larger doses, as noted for amino acids. At 12 h post-injection, there is evidence of a deviant behavior in energy metabolism (glucose, creatine), cell membrane metabolism (choline and PC), and glycosylation metabolism (UDP-GlcA), in addition to a large number of still unassigned compounds. This strong deviant signature is mostly reversed to controls at 48 h, with the exception of a few compounds, including alanine, tyrosine, and the three BCCAs. Increased levels of the latter were particularly noted to be a strong excretory marker of cDDP renal toxicity [42], consistent with impaired reabsorption of amino acids and their subsequent decreased levels in the tissue, which we confirm for cDDP doses of a few mg/kg. Alanine decrease stands out as a possible global indicator of longer-term cDDP-induced nephrotoxicity, since it was also noted by MS metabolomics of kidney extracts in melanoma tumor-bearing animals treated with cDDP for 21 days [22]. Deviant levels of hypoxanthine, fumarate, m-inositol, NAM (vitamin B3), and TMA remain at 48 h as relatively late effects of cDDP. The latter two compounds may be indicative of an abnormal gut microflora metabolic response to cDDP, consistent with previous reports of dysbiotic intestinal microflora in relation to kidney disease [52].
The early amino acid depletion in kidney is accompanied by a qualitative (non-significant) tendency for amino acids decrease also in liver, although statistically relevant effects are few, and are related to decreased lactate levels and a few other unassigned resonances. The major impact of cDDP on liver takes place at 12 h and 48 h, but in a distinct way: at 12 h, ketone bodies (acetone and 3-HBA) seem to be used up for enhanced energy production; this may justify the enhanced levels of AMP/ATP/IMP/NAD+ noted at 48 h post-injection. The profile in unassigned resonances shows a general reversal from 12 to 48 h, which in itself is interesting as part of the characteristic signatures, thus supporting further endeavors to identify these resonances. In the case of liver, this is, to the best of our knowledge, the first NMR metabolomics study of liver tissue exposed to cDDP; the only comparison possible with previous results refers to the human liver L02 cell line [27]. In cells, however, cDDP seems to induce extensive changes in many amino acid levels, not observed here, while no changes are noted in ketone bodies. The same study does show some changes in lactate and GSH levels depending on cDDP doses, which suggests that lactate/GSH changes may be translational to in vivo.
Finally, in tandem with the early effects on kidney and liver discussed above, the sole early change in breast tissue (increased ketone body 3-HBA) may reflect a deviant shift to enhanced fatty acid oxidation which quickly leveled off subsequently. The strong nitrogenated bases decrease at 12 h, and the usage of glutamine at 48 h requires further investigation, although we propose at this point that nucleotide metabolism and glutamine usage, possibly to sustain energy production, may contribute to the metabolic breast tissue signature of cDDP response.

4. Materials and Methods

4.1. Chemicals

Cisplatin (cis-dichlorodiammine platinum (II), 99.9%) was purchased from Sigma-Aldrich (Sintra, Portugal). Euthasol® solution (400 mg/mL pentobarbital sodium) was from Le Vet (Oudewater, Netherlands). All reagents were of analytical grade.

4.2. Ethical Considerations

Handling and care of animals were consistent with Portuguese (Decreto-Lei n.°113/2013) and European legislation (Directive 2010/63/EU) on the protection of animals used for scientific purposes, and were in agreement with the recommendations in the Guide for Care and Use of Laboratory Animals of the National Institutes of Health (NIH). The study protocol was approved by the Ethics Committee for Animal Experimentation of the Faculty of Pharmacy of the University of Porto, Porto, Portugal (Permit Number: 25-10-2015), and by the Ethics Committee and the Organ Responsible for the Welfare of Animals of ICBAS-UP, Porto, Portugal (Permit number 134/2015). The ARRIVE Guidelines were followed for reporting in vivo experiments [53].

4.3. Animals

Six-week old, Specific-Pathogen-Free (SPF), female BALB/cByJ mice (30 animals in total) were purchased from Charles River Laboratories (L’Arbresle, France) and acclimatized for 1 week at ICBAS-UP Rodent Animal House Facility (Porto, Portugal). Animals were randomly distributed in groups of five per cage (in individually ventilated cages) with enrichment material (corncob bedding, paper role tube, and one large sheet of tissue paper for nesting), and were housed in a SPF environment with ad libitum access to water and standard pellet food under controlled 12 h light/dark cycles (lights on at 7.00 AM), temperature (22 ± 2 °C), and humidity (50 ± 10%).

4.4. In Vivo Experimental Procedures

After 1 week of acclimatization, animals were randomly allocated into two groups (15 animals per group): cDDP and vehicle (PBS). At the beginning and end of the experiments, the animals weighed 19.9 ± 1.5 g and 20.0 ± 1.6 g, respectively. CDDP 3.5 mg/kg solution in PBS: H2PO4 1.5 mM, Na2HPO4 4.3 mM, KCl 2.7 mM, NaCl 150 mM, pH 7.4 and PBS alone were administered in a single dose, via intraperitoneal injection (volume of injection = 200 µL) in the cDDP and vehicle groups, respectively. After drug/vehicle administration, five animals from each group were euthanized, at each time-point (1, 12, and 48 h), with pentobarbital intraperitoneal injection (120 mg/kg) followed by cardiac puncture. The left kidney, median liver lobe, and thoracic mammary glands were excised and snap frozen in liquid nitrogen and stored at −80 °C for NMR analysis. All injected solutions were sterile filtered. No biochemical measurements were performed to detect kidney damage, as the volume of blood obtained was too small (under 400 µL and no sample, in some cases). In any case, all animals showed no loss of weight (see the weight values above), no evidence of decreased food or drink intake, and no behavioral changes whatsoever, which indicated the absence of significant kidney damage. This was also supported by identical observations in a small pilot study (n = 2) of our own with animals injected with 3 cycles of cDDP (3.5 mg/kg per cycle), and by a previous report that blood urea nitrogen (BUN) in mice increased after 24 h for a 30 mg/kg b.w. cDDP dose [3].

4.5. Sample Preparation for NMR

Tissues were weighted (ca. 50, 60, and 35 mg for kidney (superior half of organ), liver (median liver lobe), and breast tissue (predominantly consisting of mammary gland), respectively) and ground by mechanical maceration with cooled liquid N2 using a pestle and mortar [54,55,56]. Tissue extracts were then prepared using the biphasic methanol/chloroform/water (2.0:2.0:1.0) procedure [57]. Briefly, samples were homogenized in 8.0 mL/g cold 80% methanol, 4.0 mL/g cold chloroform, and 2 mg/L cold water, vortexed for 60 s and kept on ice for 10 min [57]. Samples were centrifuged (8000 rpm, 5 min, 23 °C), and polar phases were removed, vacuum-dried, and stored at −80 °C until analysis. Before NMR acquisition, aqueous extracts were suspended in 650 µL of 100 mM sodium phosphate buffer (pH 7.4, in D2O containing 0.25% 3-(trimethylsilyl)-propionic-2,2,3,3-d4 acid (TSP, for chemical shift referencing)), homogenized, and 600 µL of the resulting solution was then transferred to NMR tubes.

4.6. NMR Measurements

All NMR spectra were acquired on a Bruker AVANCE III spectrometer operating at 500.13 MHz for 1H. The unidimensional (1D) proton NMR spectra of the aqueous extracts were recorded at 298 K, using the “noesypr1d” pulse sequences (Bruker library), with 2.34 s acquisition time, 2 s relaxation delay, 512 scans, 7002.801 Hz spectral width, and 32 k data points. Each free-induction decay was zero-filled to 64 k points and multiplied by a 0.3 Hz eponential function prior to Fourier transformation. Spectra were manually phased, baseline-corrected, and chemical-shift referenced to TSP. Bidimensional (2D) NMR homonuclear (TOCSY) and heteronuclear (HSQC) spectra were recorded for selected samples to aid spectral assignment. Peak assignment was based on comparison with data available on Bruker BBIOREFCODE spectral database and the human metabolome database (HMDB) [58], as well as on existing literature.

4.7. Data Processing and STATISTICS

The 1D NMR spectra were converted into data matrices (AMIX-viewer 3.9.14, BrukerBiospin, Rheinstetten, Germany) prior multivariate analysis, and exclusion of the water (δ 4.5–5.2) and methanol (singlet at δ 3.36) regions was applied. Spectra were aligned by recursive, segment-wise peak alignment (RSPA) (Matlab 8.3.0, The MathWorks Inc., Natick, Massachusetts, USA) and normalized to spectral total area, to account for sample concentration differences. Multivariate analysis of the 1D spectra was carried out using principal component analysis (PCA) and partial least-squares discriminant analysis (PLS-DA) after unit variance (UV) scaling of the spectra (SIMCA-P 11.5; Umetrics, Umeå, Sweden). PLS-DA loadings were back-transformed, multiplying each variable by its standard deviation, and colored according to variable importance to the projection (VIP) (Matlab 8.3.0, The MathWorks Inc., Natick, Massachusetts, USA). PLS-DA models were considered statistically robust for predictive power (Q2) values ≥ 0.05. The relevant resonances identified in PLS-DA loadings plots were integrated (Amix 3.9.5, Bruker BioSpin, Rheinstetten, Germany), normalized, and variations assessed by univariate analysis (Shapiro-Wilk test used to assess if data was normally distributed, Student’s t test used for normally-distributed data, or Wilcoxon test used for non-normally distributed data) (R-statistical software). Significantly changed metabolites (p < 0.05) were identified, and the corresponding effect size values calculated [50]. False discovery rate (FDR) correction based on the Benjamini and Hochberg method [51] was used to correct p-values for multiple comparisons.

5. Conclusions

This work describes, to the authors´ knowledge, the first NMR metabolomics analysis of mice kidney, liver, and breast tissue in response to cisplatin exposure (single 3.5 mg/kg injection), to search for early metabolic signatures of cDDP biotoxicity. The results obtained show that the largest metabolic impact occurs for kidney, as expected, as early as at 1 h post-injection. Kidney shows a marked depletion in amino acids, comprising a 14-metabolite signature at 1 h. The highest number of variations is noted at 12 h after drug administration, although many of these recover to control levels at 48 h, with the exception of persistent changes in 10 specific metabolites (namely, alanine, BCCAs, tyrosine, hypoxanthine, fumarate m-inositol, niacinamide, and trimethylamine). Liver and breast tissue also show more marked metabolic changes at 12 h, with liver standing out for its persistent increase in the levels of several nucleotides as a longer-term cisplatin effect. The present study, however, presents some limitations, namely regarding some still unassigned resonances included in all tissue signatures (an observation transversal to most metabolomic studies, either NMR- or MS-based), and the lack of biochemical measurements of biotoxicity in both kidney and liver tissues (although there was no behavioral evidence of significant organ damage). In spite of these, this work successfully reports the very early (1–12 h) cisplatin effects in kidney, liver, and breast tissue metabolism, thus contributing to our understanding of multi-organ cDDP biotoxicity.

Supplementary Materials

The following is available online at https://www.mdpi.com/2218-1989/9/11/279/s1, Table S1: List of compounds and corresponding spin systems identified in the 500 MHz 1H NMR spectra of aqueous extracts of kidney, liver and breast tissue of healthy BALB/c mice.

Author Contributions

The authors of this paper contributed for the work carried out as follows: conceptualization, A.M.G., M.P.M.M., A.L.M.B.d.C. and C.D.; methodology, C.D., M.V., S.G.-M., R.A., A.M.G., T.J.C.; validation, C.D., M.V., S.G.-M.; formal analysis, R.A., T.J.C. and A.M.G.; resources, C.D., M.P.M.M. and A.M.G.; data curation, T.J.C. and A.M.G.; writing—original draft preparation, T.J.C., C.D., M.V., S.G.-M. and A.M.G.; writing—review and editing, T.J.C., M.V., S.G.-M., M.P.M.M., A.L.M.B.d.C., C.D. and A.M.G.; supervision, C.D., M.P.M.M. and A.M.G.; project administration, C.D., M.P.M.M. and A.M.G.; funding acquisition, C.D., M.P.M.M., A.L.M.B.d.C. and A.M.G.

Funding

This research was developed within the scope of the CICECO—Aveiro Institute of Materials, FCT UID/CTM/50011/2019 project (AMG), financed by national funds through the FCT/MCTES and by the Portuguese Foundation for Science and Technology (FCT) through LAQV/REQUIMTE FCT UID/QUI/50006/2019 (CD), UID/MULTI/00070/2019 (MPM), POCI-01-0145-FEDER-0016786 and Centro-01-0145-FEDER-029956 (co-financed by COMPETE 2020, Portugal 2020 and European Community through FEDER). The authors acknowledge further financial support from FCT through PTDC/QEQ-MED/1890/2014, within Project 3599–to Promote Scientific Production and Technological Development as well as the formation of thematic networks (3599-PPCDT)–jointly financed by the European Community Fund FEDER. We also acknowledge the Portuguese National NMR Network (PTNMR), supported by FCT funds as the NMR spectrometer used is part of PTNMR and partially supported by Infrastructure Project Nº 022161 (co-financed by FEDER through COMPETE 2020, POCI and PORL and FCT through PIDDAC). M.V. thanks FCT and PhD Program in Medicines and Pharmaceutical Innovation (i3DU) for PhD grant PD/BD/135460/2017 and T.J.C. thanks FCT for her PhD grant SFRH/BD/145920/2019, both grants funded by the European Social Fund of the European Union and national funds FCT/MCTES.

Acknowledgments

M Spraul, Bruker BioSpin (Germany) is thanked for access to software and Bruker spectral database. The authors thank Sara Capas-Peneda for procedural support and valuable recommendations during in vivo experiments. We thank Joana B. Sousa, Maria do Céu Pereira, Mónica Caldas and João S. Marques for their help during tissue collections.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Dasari, S.; Tchounwou, P.B. Cisplatin in cancer therapy: Molecular mechanisms of action. Eur. J. Pharmacol. 2014, 740, 364–378. [Google Scholar] [CrossRef] [PubMed]
  2. Manohar, S.; Leung, N. Cisplatin nephrotoxicity: A review of the literature. J. Nephrol. 2018, 31, 15–25. [Google Scholar] [CrossRef] [PubMed]
  3. Irie, M.; Hayakawa, E.; Fujimura, Y.; Honda, Y.; Setoyama, D.; Wariishi, H.; Hyodo, F.; Miura, D. Analysis of spatiotemporal metabolomic dynamics for sensitively monitoring biological alterations in cisplatin-induced acute kidney injury. Biochem. Biophys. Res. Commun. 2018, 496, 140–146. [Google Scholar] [CrossRef] [PubMed]
  4. Achkar, I.W.; Abdulrahman, N.; Al-Sulaiti, H.; Joseph, J.M.; Uddin, S.; Mraiche, F. Cisplatin based therapy: The role of the mitogen activated protein kinase signaling pathway. J. Transl. Med. 2018, 16. [Google Scholar] [CrossRef]
  5. Perše, M.; Veceric-Haler, C. Cisplatin-Induced Rodent Model of Kidney Injury: Characteristics and Challenges. BioMed Res. Int. 2018, 2018, 1462802. [Google Scholar] [CrossRef]
  6. Khoo, S.H.G.; Al-Rubeai, M. Metabolomics as a Complementary Tool in Cell Culture. Biotechnol. Appl. Biochem. 2007, 47, 71–84. [Google Scholar] [CrossRef]
  7. Cuperlovic-Culf, M.; Barnett, A.D.; Culf, A.S.; Chute, I. Cell Culture Metabolomics: Applications and Future Directions. Drug Discov. Today 2010, 15, 610–621. [Google Scholar] [CrossRef]
  8. Coen, M.; Holmes, E.; Lindon, J.C.; Nicholson, J.K. NMR-Based Metabolic Profiling and Metabonomic Approaches to Problems in Molecular Toxicology. Chem. Res. Toxicol. 2008, 21, 9–27. [Google Scholar] [CrossRef]
  9. Robertson, D.G.; Watkins, P.B.; Reily, M.D. Metabolomics in Toxicology: Preclinical and Clinical Applications. Toxicol. Sci. 2011, 120, S146–S170. [Google Scholar] [CrossRef]
  10. Lee, S.; Kang, S.; Park, S. Comparison of Metabolic Profiles of Normal and Cancer Cells in Response to Cytotoxic Agents. J. Korean Magn. Reson. Soc. 2017, 21, 31–43. [Google Scholar] [CrossRef]
  11. Duarte, I.F.; Ladeirinha, A.F.; Lamego, I.; Gil, A.M.; Carvalho, L.; Carreira, I.M.; Melo, J.B. Potential Markers of Cisplatin Treatment Response Unveiled by NMR Metabolomics of Human Lung Cells. Mol. Pharm. 2013, 10, 4242–4251. [Google Scholar] [CrossRef] [PubMed]
  12. Maria, R.M.; Altei, W.F.; Selistre-de-Araujo, H.S.; Colnago, L.A. Impact of chemotherapy on metabolic reprogramming: Characterization of the metabolic profile of breast cancer MDA-MB-231 cells using H-1 HR-MAS NMR spectroscopy. J. Pharm. Biomed Anal. 2017, 146, 324–328. [Google Scholar] [CrossRef] [PubMed]
  13. Mandujano-Tinoco, E.A.; Gallardo-Perez, J.C.; Marin-Hernandez, A.; Moreno-Sanchez, R.; Rodriguez-Enriquez, S. Anti-mitochondrial therapy in human breast cancer multi-cellular spheroids. Biochim. Biophys. Acta 2013, 1833, 541–551. [Google Scholar] [CrossRef] [PubMed]
  14. Lamego, I.; Marques, M.P.M.; Duarte, I.F.; Martins, A.S.; Oliveira, H.; Gil, A.M. Impact of the Pd(2)Spermine Chelate on Osteosarcoma Metabolism: An NMR Metabolomics Study. J. Proteome Res. 2017, 16, 1773–1783. [Google Scholar] [CrossRef] [PubMed]
  15. Duarte, I.F.; Lamego, I.; Marques, J.; Marques, M.P.M.; Blaise, B.J.; Gil, A.M. Nuclear Magnetic Resonance (NMR) Study of the Effect of Cisplatin on the Metabolic Profile of MG-63 Osteosarcoma Cells. J. Proteome Res. 2010, 9, 5877–5886. [Google Scholar] [CrossRef] [PubMed]
  16. Pan, X.Y.; Wilson, M.; Mirbahai, L.; McConville, C.; Arvanitis, T.N.; Griffin, J.L.; Kauppinen, R.A.; Peet, A.C. In Vitro Metabonomic Study Detects Increases in UDP-GlcNAc and UDP-GalNAc, as Early Phase Markers of Cisplatin Treatment Response in Brain Tumor Cells. J. Proteome Res. 2011, 10, 3493–3500. [Google Scholar] [CrossRef]
  17. Alonezi, S.; Tusiimire, J.; Wallace, J.; Dufton, M.J.; Parkinson, J.A.; Young, L.C.; Clements, C.J.; Park, J.K.; Jeon, J.W.; Ferro, V.A.; et al. Metabolomic Profiling of the Synergistic Effects of Melittin in Combination with Cisplatin on Ovarian Cancer Cells. Metabolites 2017, 7, 14. [Google Scholar] [CrossRef]
  18. Poisson, L.M.; Munkarah, A.; Madi, H.; Datta, I.; Hensley-Alford, S.; Tebbe, C.; Buekers, T.; Giri, S.; Rattan, R. A metabolomic approach to identifying platinum resistance in ovarian cancer. J. Ovarian Res. 2015, 8, 13. [Google Scholar] [CrossRef]
  19. De Castro, F.; Benedetti, M.; Antonaci, G.; Del Coco, L.; De Pascali, S.A.; Muscella, A.; Marsigliante, S.; Fanizzi, F.P. Response of Cisplatin Resistant Skov-3 Cells to [Pt(O,O ‘-Acac)(gamma-Acac)(DMS)] Treatment Revealed by a Metabolomic H-1-NMR Study. Molecules 2018, 23, 2301. [Google Scholar] [CrossRef]
  20. Obrist, F.; Michels, J.; Durand, S.; Chery, A.; Pol, J.; Levesque, S.; Joseph, A.; Astesana, V.; Pietrocola, F.; Wu, G.S.; et al. Metabolic vulnerability of cisplatin-resistant cancers. EMBO J. 2018, 37. [Google Scholar] [CrossRef]
  21. Zong, L.; Pi, Z.F.; Liu, S.; Liu, Z.Q.; Song, F.R. Metabolomics analysis of multidrug-resistant breast cancer cells in vitro using methyl-tert-butyl ether method. RSC Adv. 2018, 8, 15831–15841. [Google Scholar] [CrossRef]
  22. Chao, W.-W.; Cheng, Y.-W.; Chen, Y.-R.; Lee, S.-H.; Cjou, C.-Y.; Shyur, L.-F. Phyto-sesquiterpene lactone deoxyelephantopin and cisplatin synergistically suppress lung metastasis of B16 melanoma in mice with reduced nephrotoxicity. Phytomedicine 2019, 56, 194–206. [Google Scholar] [CrossRef] [PubMed]
  23. Tian, Y.H.; Wang, Z.J.; Liu, X.H.; Duan, J.C.; Feng, G.S.; Yin, Y.X.; Gu, J.; Chen, Z.L.; Gao, S.G.; Bai, H.; et al. Prediction of Chemotherapeutic Efficacy in Non-Small Cell Lung Cancer by Serum Metabolomic Profiling. Clin. Cancer Res. 2018, 24, 2100–2109. [Google Scholar] [CrossRef] [PubMed]
  24. Xu, S.; Zhou, Y.P.; Geng, H.; Song, D.D.; Tang, J.; Zhu, X.M.; Yu, D.; Hu, S.; Cui, Y.F. Serum Metabolic Profile Alteration Reveals Response to Platinum-Based Combination Chemotherapy for Lung Cancer: Sensitive Patients Distinguished from Insensitive ones. Sci. Rep. 2017, 7, 17524. [Google Scholar] [CrossRef]
  25. Tandberg, D.J.; Holt, T.; Kelsey, C.R. Plasma Metabolites and Risk of Radiation-induced Esophagitis: A Secondary Analysis from a Prospective Study. Anticancer Res. 2017, 37, 719–725. [Google Scholar] [CrossRef]
  26. Doskocz, M.; Marchewka, Z.; Jez, M.; Passowicz-Muszynska, E.; Dlugosz, A. Preliminary Study on J-Resolved NMR Method Usability for Toxic Kidney’s Injury Assessment. Adv. Clin. Exp. Med. 2015, 24, 629–635. [Google Scholar] [CrossRef]
  27. Liu, S.; Wang, W.; Zhou, X.; Gu, R.; Ding, Z. Dose responsive effects of cisplatin in L02 cells using NMR-based metabolomics. Environ. Toxicol. Pharmacol. 2014, 37, 150–157. [Google Scholar] [CrossRef]
  28. Wilmes, A.; Bielow, C.; Ranninger, C.; Bellwon, P.; Aschauer, L.; Limonciel, A.; Chassaigne, H.; Kristl, T.; Aiche, S.; Huber, C.G.; et al. Mechanism of cisplatin proximal tubule toxicity revealed by integrating transcriptomics, proteomics, metabolomics and biokinetics. Toxicol. In Vitro 2015, 30, 117–127. [Google Scholar] [CrossRef]
  29. Heyman, S.N.; Rosen, S.; Rosenberger, C. Animal models of renal dysfunction: Acute kidney injury. Expert Opin. Drug Discov. 2009, 4, 629–641. [Google Scholar] [CrossRef]
  30. Portilla, D.; Li, S.; Nagothu, K.K.; Megyesi, J.; Kaissling, B.; Schnackenberg, L.; Safirstein, R.L.; Beger, R.D. Metabolomic study of cisplatin-induced nephrotoxicity. Kidney Int. 2006, 69, 2194–2204. [Google Scholar] [CrossRef]
  31. Zhang, P.; Chen, J.; Wang, Y.; Huang, Y.; Tian, Y.; Zhang, Z.; Xu, F. Discovery of Potential Biomarkers with Dose- and Time-Dependence in Cisplatin Induced Nephrotoxicity Using Metabolomics Integrated with a Principal Component-Based Area Calculation Strategy. Chem. Res. Toxicol. 2016, 29, 776–783. [Google Scholar] [CrossRef]
  32. Li, Y.; Zhang, X.; Zhou, H.; Fan, S.; Wang, Y.; Zhang, L.; Ju, L.; Wu, X.; Wu, H.; Zhang, Y. Metabonomics study on nephrotoxicity induced by intraperitoneal and intravenous cisplatin administration using rapid resolution liquid chromatography coupled with quadrupole-time-of- flight mass spectrometry (RRLC-Q-TOF-MS). RSC Adv. 2014, 4, 8260. [Google Scholar] [CrossRef]
  33. Zhang, P.; Li, W.; Chen, J.; Li, R.; Zhang, Z.; Huang, Y.; Xu, F. Branched-Chain Amino Acids as Predictors for Individual Differences of Cisplatin Nephrotoxicity in Rats: A Pharmacometabonomics Study. J. Proteome Res. 2017, 16, 1753–1762. [Google Scholar] [CrossRef] [PubMed]
  34. Uehara, T.; Horinouchid, A.; Morikawaa, Y.; Tonomuraa, Y.; Minamie, K.; Onob, A.; Yamatec, J.; Yamadab, H.; Ohnof, Y.; Urushidanib, T. Identification of metabolomic biomarkers for drug-induced acute kidney injury in rats. J. Appl. Toxicol. 2014, 34, 1087–1095. [Google Scholar] [CrossRef] [PubMed]
  35. Ezaki, T.; Nishiumi, S.; Azumaa, T.; Yoshida, M. Metabolomics for the early detection of cisplatin induced nephrotoxicity. Toxicol. Res. 2017, 6, 843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Xu, E.Y.; Perlina, A.; Vu, H.; Troth, S.P.; Brennan, R.J.; Aslamkhan, A.G.; Xu, Q. Integrated Pathway Analysis of Rat Urine Metabolic Profiles and Kidney Transcriptomic Profiles To Elucidate the Systems Toxicology of Model Nephrotoxicants. Chem. Res. Toxicol. 2008, 21, 1548–1561. [Google Scholar] [CrossRef]
  37. Ryu, S.H.; Lee, J.D.; Kim, J.W.; Kim, S.; Kim, K.B. 1H NMR toxicometabolomics following cisplatin-induced nephrotoxicity in male rats. J. Toxicol. Sci. 2019, 44, 57–71. [Google Scholar] [CrossRef]
  38. Espandiari, P.; Rosenzweig, B.; Zhang, J.; Zhou, Y.; Schnackenberg, L.; Vaidya, V.S.; Goering, P.L.; Brown, R.P.; Bonventre, V.; Mahjoob, K.; et al. Age-related differences in susceptibility to cisplatin-induced renal toxicity. J. Appl. Toxicol. 2009, 30, 172–182. [Google Scholar] [CrossRef] [Green Version]
  39. Won, A.J.; Kim, S.; Kim, Y.G.; Kim, K.-B.; Choi, W.S.; Kacew, S.; Kim, K.S.; Jung, J.H.; Lee, B.M.; Kimb, S.; et al. Discovery of urinary metabolomic biomarkers for early detection of acute kidney injury. Mol. Biosyst. 2016, 12, 133. [Google Scholar] [CrossRef]
  40. Portilla, D.; Schnackenberg, L.; Beger, R.D. Metabolomics as an Extension of Proteomic Analysis: Study of acute kidney injury. Semin. Nephrol. 2007, 27, 609–620. [Google Scholar] [CrossRef] [Green Version]
  41. Kwon, H.K.; Kim, M.; Wen, H.; Kang, S.; Yang, H.-J.; Choi, M.-J.; Lee, H.S.; Choi, D.; Park, I.S.; Suh, I.J.; et al. Predicting idiopathic toxicity of cisplatin by a pharmacometabonomic approach. Kidney Int. 2011, 79, 529–537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Boudonck, K.J.; Mitechell, M.W.; Német, L.; Keresztes, L.; Nyska, A.; Shinar, D.; Rosenstock, M. Discovery of Metabolomics Biomarkers for Early Detection of Nephrotoxicity. Toxicol. Pathol. 2009, 37, 280–292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zhang, P.; Chen, J.-Q.; Huang, W.-Q.; Li, W.; Huang, Y.; Zhang, Z.-J.; Xu, F.-G. Renal Medulla is More Sensitive to Cisplatin than Cortex Revealed by Untargeted Mass Spectrometry-Based Metabolomics in Rats. Sci. Rep. 2017, 7, 44804. [Google Scholar] [CrossRef] [PubMed]
  44. Zhong, F.; Liu, X.; Zhou, Q.; Hao, X.; Lu, Y.; Guo, S.; Wang, W.; Lin, D.; Chen, N. 1H NMR spectroscopy analysis of metabolites in the kidneys provides new insight into pathophysiological mechanisms: Applications for treatment with Cordyceps sinensis. Nephrol. Dial. Transplant. 2012, 27, 556–565. [Google Scholar] [CrossRef] [Green Version]
  45. Ling, Y.S.; Liang, H.-J.; Chung, M.-H.; Lin, M.-H.; Lin, C.-Y. NMR- and MS-based metabolomics: Various organ responses following naphthalene intervention. Mol. BioSyst. 2014, 10, 1918–1931. [Google Scholar] [CrossRef]
  46. Zheng, H.; Ni, Z.; Cai, A.; Zhang, X.; Chen, J.; Shu, Q.; Gao, H. Balancing metabolome coverage and reproducibility for untargeted NMR-based metabolic profiling in tissue samples through mixture design methods. Anal. Bioanal. Chem. 2018, 410, 7783–7792. [Google Scholar] [CrossRef]
  47. Jiang, L.; Si, Z.-H.; Li, M.-H.; Zhao, H.; Fu, Y.-H.; Xing, Y.-X.; Hong, W.; Ruan, L.Y.; Li, P.-M.; Wang, J.-S. 1H NMR-based metabolomics study of liver damage induced by ginkgolic acid (15:1) in mice. J. Pharm. Biomed. Anal. 2017, 136, 44–54. [Google Scholar] [CrossRef]
  48. Chen, M.; Zheng, H.; Xu, M.; Zhao, L.; Zhang, Q.; Song, J.; Zhao, Z.; Lu, S.; Weng, Q.; Wu, X.; et al. Changes in hepatic metabolic profile during the evolution of STZ-induced diabetic rats via an 1H NMR-based metabonomic investigation. Biosci. Rep. 2019, 39, 12. [Google Scholar] [CrossRef] [Green Version]
  49. Feng, X.; Xu, H.; Chen, J.-F.; Ruan, L.-Y.; Zhao, W.-L.; Meng, H.-H.; Liu, W.-Y.; Zhao, W.-L.; Zheng, Q.; Liu, Z.-C.; et al. Potential hepatoxicity risk of the shell of Herpetospermum caudigerum Wall in rats based on 1H-NMR metabonomicsXin. J. Pharm. Biomed. Anal. 2019, 176, 112800. [Google Scholar] [CrossRef]
  50. Berben, L.; Sereika, S.M.; Engberg, S. Effect size estimation: Methods and examples. Int. J. Nurs. Stud. 2012, 49, 1039–1047. [Google Scholar] [CrossRef]
  51. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Stat. Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  52. Guldris, S.C.; Parra, E.G.; Amenós, A.C. Gut microbiota in chronic kidney disease. Nefrologia 2017, 37, 9–19. [Google Scholar] [CrossRef]
  53. Kikenny, C.; Browne, W.J.; Cuthill, I.C.; Emerson, M.; Altman, D.G. Improving Bioscience Research Reporting: The ARRIVE Guidelines for Reporting Animal Research. PLoS Biol. 2010, 8, e1000412. [Google Scholar] [CrossRef]
  54. Beckonert, O.; Keun, H.C.; Ebbels, T.M.; Bundy, J.; Holmes, E.; Lindon, J.C.; Nicholson, J.K. Metabolic profiling, metabolomic and metabonomic procedures for NMR spectroscopy of urine, plasma, serum and tissue extracts. Nat. Protoc. 2007, 2, 2692–2703. [Google Scholar] [CrossRef]
  55. Lin, C.Y.; Wu, H.; Tjeerdema, R.S.; Viant, M.R. Evaluation of metabolite extraction strategies from tissue samples using NMR metabolomics. Metabolomics 2007, 3, 55–67. [Google Scholar] [CrossRef]
  56. Le Belle, J.E.; Harris, N.G.; Williams, S.R.; Bhakoo, K.K. A comparison of cell and tissue extraction techniques using high-resolution 1H-NMR spectroscopy. NMR Biomed. 2002, 15, 37–44. [Google Scholar] [CrossRef]
  57. Wu, H.; Southam, A.D.; Hines, A.; Viant, M.R. High-throughput tissue extraction protocol for NMR- and MS-based metabolomics. Anal. Biochem. 2008, 372, 204–212. [Google Scholar] [CrossRef]
  58. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vázquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0—The Human Metabolome Database for 2018. Nucleic Acids Res. 2018, 46, D608–D617. [Google Scholar] [CrossRef]
Figure 1. Typical 500 MHz 1H NMR spectra of aqueous extracts of (a) kidney, (b) liver, and (c) breast tissue from the control group (BALB/c mice) 1 h after injection with phosphate buffer saline. Compound names in red represent the predominant metabolites in the NMR spectra of each type of tissue. * Cut-off spectral region due to water suppression. 3-Letter codes are used for amino acids; Acet: acetate; Ado: adenosine; ADP: adenosine diphosphate; AMP: adenosine monophosphate; ATP: adenosine triphosphate; Cho: choline; Cre: creatine; Glc: glucose; GPC: glycerophosphocholine; GSH: glutathione (reduced); Hip: hippurate; HX: hypoxanthine; Ino: inosine; Lac: lactate; NAD+: nicotinamide adenine dinucleotide (reduced); NAM: niacinamide; PC: phosphocholine; Tau: taurine; TMAO: trimethylamine N-oxide; UDP-GlcA: uridine diphosphate-glucuronate; Urd: uridine.
Figure 1. Typical 500 MHz 1H NMR spectra of aqueous extracts of (a) kidney, (b) liver, and (c) breast tissue from the control group (BALB/c mice) 1 h after injection with phosphate buffer saline. Compound names in red represent the predominant metabolites in the NMR spectra of each type of tissue. * Cut-off spectral region due to water suppression. 3-Letter codes are used for amino acids; Acet: acetate; Ado: adenosine; ADP: adenosine diphosphate; AMP: adenosine monophosphate; ATP: adenosine triphosphate; Cho: choline; Cre: creatine; Glc: glucose; GPC: glycerophosphocholine; GSH: glutathione (reduced); Hip: hippurate; HX: hypoxanthine; Ino: inosine; Lac: lactate; NAD+: nicotinamide adenine dinucleotide (reduced); NAM: niacinamide; PC: phosphocholine; Tau: taurine; TMAO: trimethylamine N-oxide; UDP-GlcA: uridine diphosphate-glucuronate; Urd: uridine.
Metabolites 09 00279 g001
Figure 2. PCA (left) and PLS-DA (right) scores scatter plots for 1H NMR spectra of polar extracts of all time-course samples of (a) kidney, (b) liver, and (c) breast tissue of BALB/c mice.
Figure 2. PCA (left) and PLS-DA (right) scores scatter plots for 1H NMR spectra of polar extracts of all time-course samples of (a) kidney, (b) liver, and (c) breast tissue of BALB/c mice.
Metabolites 09 00279 g002
Figure 3. Scores scatter plots from PCA (left) and PLS-DA (middle) analysis obtained for 1H NMR spectra of (a) kidney, (b) liver, and (c) breast tissue of BALB/c mice at 12 h post-injection with cisplatin. LV1 loadings plots from PLS-DA analysis (right), colored according to variable importance to the projection (VIP) and with main peak assignment indicated. 3-Letter codes are used for amino acids; 3-HBA: 3-hydroxybutyrate; Ado: adenosine; ADP: adenosine monophosphate; AMP: adenosine monophosphate; ATP: adenosine monophosphate; Cho: choline; Cre: creatine; Glc: glucose; GPC: glycerophosphocholine; Ino: inosine; Lac: lactate; PC: phosphocholine; Tau: taurine.
Figure 3. Scores scatter plots from PCA (left) and PLS-DA (middle) analysis obtained for 1H NMR spectra of (a) kidney, (b) liver, and (c) breast tissue of BALB/c mice at 12 h post-injection with cisplatin. LV1 loadings plots from PLS-DA analysis (right), colored according to variable importance to the projection (VIP) and with main peak assignment indicated. 3-Letter codes are used for amino acids; 3-HBA: 3-hydroxybutyrate; Ado: adenosine; ADP: adenosine monophosphate; AMP: adenosine monophosphate; ATP: adenosine monophosphate; Cho: choline; Cre: creatine; Glc: glucose; GPC: glycerophosphocholine; Ino: inosine; Lac: lactate; PC: phosphocholine; Tau: taurine.
Metabolites 09 00279 g003
Figure 4. Heatmap illustrating the metabolic variations in (a) kidney, (b) liver, and (c) breast tissue of BALB/c mice at 1 h, 12 h, and 48 h post-injection with cisplatin, comparatively to the controls group. The scale is color-coded as a function of Effect Size, from minimum (dark blue) to maximum (dark red) values. Lines and columns represent metabolites and subjects, respectively. Asterisks represent the confidence interval of each variation: * p-value < 5.0 × 10−2; ** p-value < 1.0 × 10−2; *** p-value < 1.0 × 10−3; **** p-value < 1.0 × 10−4. 3-Letter codes are used for amino acids; ADP: adenosine diphosphate; AMP: adenosine monophosphate; Cho: choline; Glc: glucose; HX: hypoxanthine; NAD+: nicotinamide adenine dinucleotide (reduced); PC: phosphocholine; UMP: uridine monophosphate.
Figure 4. Heatmap illustrating the metabolic variations in (a) kidney, (b) liver, and (c) breast tissue of BALB/c mice at 1 h, 12 h, and 48 h post-injection with cisplatin, comparatively to the controls group. The scale is color-coded as a function of Effect Size, from minimum (dark blue) to maximum (dark red) values. Lines and columns represent metabolites and subjects, respectively. Asterisks represent the confidence interval of each variation: * p-value < 5.0 × 10−2; ** p-value < 1.0 × 10−2; *** p-value < 1.0 × 10−3; **** p-value < 1.0 × 10−4. 3-Letter codes are used for amino acids; ADP: adenosine diphosphate; AMP: adenosine monophosphate; Cho: choline; Glc: glucose; HX: hypoxanthine; NAD+: nicotinamide adenine dinucleotide (reduced); PC: phosphocholine; UMP: uridine monophosphate.
Metabolites 09 00279 g004
Table 1. List of metabolite variations in kidney of mice exposed to cDDP at 1, 12, and 48 h post-injection, relative to the controls. Only variations with 95% significance level (p-value < 0.05) are shown; ES, effect size and corresponding error [50].
Table 1. List of metabolite variations in kidney of mice exposed to cDDP at 1, 12, and 48 h post-injection, relative to the controls. Only variations with 95% significance level (p-value < 0.05) are shown; ES, effect size and corresponding error [50].
MetaboliteδH ppm (multiplicity)1 h12 h48 h
ES±Errorp-ValueES±Errorp-ValueES±Errorp-Value
3-HIBA1.09 (d)4.7±2.51.9 × 10−4 a−2.2±1.69.5 × 10−3----
ADP c8.28 (s)1.8±1.62.9 × 10−21.8±1.51.6 × 10−2----
Ala b1.48 (d)−3.0±1.99.5 × 10−3 a----−2.2±1.71.1 × 10−2
Allantoin5.39 (s)----−2.3±1.69.3 × 10−3----
AMP8.61 (s)1.8±1.53.3 × 10−21.9±1.52.8 × 10−2----
Asn b, c 2.96 (m)----−1.9±1.52.3 × 10−2----
Betaine3.90 (s)----1.8±1.53.8 × 10−2----
Cho3.21 (s)−2.1±1.61.6 × 10−2−3.2±1.92.7 × 10−3----
Cre c3.93 (s)----−2.5±1.61.1 × 10−2----
Fumarate6.52 (s)--------−3.5±2.12.7 × 10−3 a
Glc b5.23 (d)----1.9±1.52.7 × 10−2----
Hypoxanthine b8.21 (s)--------−2.6±1.83.3 × 10−2
Ile b, c1.01 (d)−3.1±1.94.1 × 10−3 a−1.9±1.53.3 × 10−2−2.2±1.71.5 × 10−2
Leu b, c0.96 (t)−3.5±2.11.2 × 10−3 a−1.8±1.53.0 × 10−2−1.8±1.62.5 × 10−2
m-Inositol3.62 (t)2.5±1.76.8 × 10−3 a1.5±1.44.4 × 10−21.6±1.53.2 × 10−2
Niacinamide b7.60 (dd)--------−1.7±1.53.6 × 10−2
PC3.22 (s)----1.9±1.51.7 × 10−2----
Phe b,c7.33 (d)−3.2±2.01.6 × 10−2--------
Tau3.42 (t)3.0±1.93.2 × 10−3 a--------
TMA2.89 (s)---------4.3±2.41.1 × 10−3 a
Tyr b, c 6.90 (d)−2.1±1.61.6 × 10−2-----3.1±1.92.5 × 10−3 a
NAD+ c8.43 (s)----1.7±1.41.6 × 10−2----
UDP-GlcA7.95 (d)----1.7±1.43.0 × 10−2----
UMP c5.99 (m)1.8±1.62.8 × 10−2--------
Val b, c 1.05 (d)−2.5±1.78.9 × 10−3 a−2.1±1.52.0 × 10−2−2.0±1.62.1 × 10−2
U10.89 (t)--------−2.5±1.81.6 × 10−2
U20.93 (s)----−2.1±1.61.2 × 10−2----
U31.62 (d)----−2.8±1.75.0 × 10−3----
U42.92 (s)----−1.7±1.43.5 × 10−2----
U53.15 (s)----−1.9±1.53.2 × 10−2----
U63.35 (s)2.0±1.61.8 × 10−21.8±1.52.2 × 10−2----
a Metabolic variations remain significant after False Discovery Rate correction [51]; b metabolite variations reported previously [3,42,43], c metabolite variations noted in different direction compared to previous report [3]. 3-Letter codes are used for amino acids; 3-HIBA: 3-hydroxyisobutyrate; ADP: adenosine diphosphate; AMP: adenosine monophosphate; Cho: choline; Cre: creatine; Glc: glucose: PC: phosphocholine; Tau: taurine; TMA: trimethylamine; NAD+: nicotinamide adenine dinucleotide; UDP-GlcA: uridine-diphosphate-glucoronate; Ui: unassigned i; UMP: uridine monophosphate.
Table 2. List of metabolite variations in liver of mice exposed to cDDP at 1, 12 and 48 h post-injection, relatively to the controls. Only variations with 95% significance level (p-value < 0.05) are shown; ES, effect size and corresponding error [50].
Table 2. List of metabolite variations in liver of mice exposed to cDDP at 1, 12 and 48 h post-injection, relatively to the controls. Only variations with 95% significance level (p-value < 0.05) are shown; ES, effect size and corresponding error [50].
MetaboliteδH ppm (multiplicity)1 h12 h48 h
ES±Errorp-ValueES±Errorp-ValueES±Errorp-Value
2-aminobutyrate 0.80 (t)----−1.6±1.44.4 × 10−2----
3-HBA1.20 (d)----−1.8±1.53.1 × 10−2----
Acetate b1.92 (s)--------−2.5±1.72.0 × 10−2
Acetone2.24 (s)----−5.1±2.65.2 × 10−5a----
ADP8.54 (s)--------2.6±1.81.6 × 10−2
AMP b4.51 (dd)--------2.8±1.83.8 × 10−3
ATP b8.52 (s)--------2.3±1.71.6 × 10−2
DMA2.73 (s)----−2.4±1.66.0 × 10−3 a----
Formate8.46 (s)----−1.9±1.57.9 × 10−3 a----
GSH2.55 (m)----−4.0±2.12.7 × 10−4 a----
His7.08 (s)----2.2±1.69.110−3----
IMP8.58 (s)--------1.7±1.53.2 × 10−2
Lactate b4.10 (q)−2.7±1.76.9 × 10−31.7±1.53.0 × 10−2----
NAD+8.43 (s)--------2.6±1.89.9 × 10−3
Tyr6.90 (d)--------1.9±1.62.3 × 10−2
UDP-GlcA7.95 (d)-----1.7±1.43.0 × 10−2----
U70.73 (s)−2.2±1.61.6 × 10−2----2.6±1.81.6 × 10−2
U80.85 (t)----−1.6±1.41.6 × 10−2----
U93.10 (d)----2.5±1.79.8 × 10−3 a----
U104.31 (d)−1.6±1.44.7 × 10−2−1.8±1.52.1 × 10−21.7±1.53.3 × 10−2
U116.18 (s)−1.7±1.53.6 × 10−2−2.6±1.77.8 × 10−3 a2.2±1.72.3 × 10−2
U128.28 (br)----−1.8±1.53.5 × 10−21.8±1.54.2 × 10−2
a Metabolic variations remaining significant after False Discovery Rate correction [51]; b metabolite variations observed in common with a previous study on the human liver cell line L02 [27] as, to our knowledge, no previous reports have been found in relation to liver tissue, either of murine or human nature; Tentative of assignment. 3-HBA: 3-hydroxybutyrate; ATP: adenosine triphosphate; DMA: dimethylamine; GSH: glutathione (reduced) IMP: inosine monophosphate; other abbreviations are as defined in Table 1.
Table 3. List of metabolite variations in breast tissue of mice exposed to cDDP at 1, 12, and 48 h post-injection, relatively to controls. Only variations with 95% significance level (p-value < 0.05) are shown; ES, effect size and corresponding error [50]; ª Metabolic variations remaining significant after False Discovery Rate correction [51]. Ado: adenosine; Ino: inosine; other abbreviations are as defined in the legends of Table 1 and Table 2.
Table 3. List of metabolite variations in breast tissue of mice exposed to cDDP at 1, 12, and 48 h post-injection, relatively to controls. Only variations with 95% significance level (p-value < 0.05) are shown; ES, effect size and corresponding error [50]; ª Metabolic variations remaining significant after False Discovery Rate correction [51]. Ado: adenosine; Ino: inosine; other abbreviations are as defined in the legends of Table 1 and Table 2.
MetaboliteδH ppm (multiplicity)1 h12 h48 h
ES±Errorp-valueES±Errorp-valueES±Errorp-value
3-HBA1.20 (d)1.8±1.54.3 × 10−2--------
Ado/Ino6.10 (d)----−2.4±1.66.1 × 10−3a----
ADP8.54 (s)----1.7±1.43.1 × 10−2 a----
Gln2.45 (m)--------−2.2±1.61.2 × 10−2
Ino8.35 (s)----−2.8±1.72.4 × 10−3 a----
U131.25 (s)----2.0±1.52.6 × 10−2 a----
U147.68 (s)--------−1.9±1.63.1 × 10−2

Share and Cite

MDPI and ACS Style

Carneiro, T.J.; Araújo, R.; Vojtek, M.; Gonçalves-Monteiro, S.; Diniz, C.; Batista de Carvalho, A.L.M.; Marques, M.P.M.; Gil, A.M. Multi-Organ NMR Metabolomics to Assess In Vivo Overall Metabolic Impact of Cisplatin in Mice. Metabolites 2019, 9, 279. https://doi.org/10.3390/metabo9110279

AMA Style

Carneiro TJ, Araújo R, Vojtek M, Gonçalves-Monteiro S, Diniz C, Batista de Carvalho ALM, Marques MPM, Gil AM. Multi-Organ NMR Metabolomics to Assess In Vivo Overall Metabolic Impact of Cisplatin in Mice. Metabolites. 2019; 9(11):279. https://doi.org/10.3390/metabo9110279

Chicago/Turabian Style

Carneiro, Tatiana J., Rita Araújo, Martin Vojtek, Salomé Gonçalves-Monteiro, Carmen Diniz, Ana L.M. Batista de Carvalho, Maria Paula M. Marques, and Ana M. Gil. 2019. "Multi-Organ NMR Metabolomics to Assess In Vivo Overall Metabolic Impact of Cisplatin in Mice" Metabolites 9, no. 11: 279. https://doi.org/10.3390/metabo9110279

APA Style

Carneiro, T. J., Araújo, R., Vojtek, M., Gonçalves-Monteiro, S., Diniz, C., Batista de Carvalho, A. L. M., Marques, M. P. M., & Gil, A. M. (2019). Multi-Organ NMR Metabolomics to Assess In Vivo Overall Metabolic Impact of Cisplatin in Mice. Metabolites, 9(11), 279. https://doi.org/10.3390/metabo9110279

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