- 1Department of Pharmacotherapy and Translational Research, College of Pharmacy, University of Florida, Gainesville, FL, United States
- 2Department of Pathology, Immunology and Laboratory Medicine, College of Medicine, University of Florida, Gainesville, FL, United States
- 3Division of Hematology and Oncology, Department of Medicine, University of Florida, Gainesville, FL, United States
- 4Division of Hematology and Oncology, Vanderbilt University Medical Center, Preston Research Building, Nashville, TN, United States
- 5Department of Malignant Hematology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, United States
- 6Cardio-Oncology Center of Excellence, Division of Cardiology, Department of Medicine, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, United States
- 7Department of Internal Medicine, College of Medicine, University of Florida, Gainesville, FL, United States
- 8UF Health Cancer Center, Gainesville, FL, United States
- 9Department of Pharmacodynamics, College of Pharmacy, University of Florida, Gainesville, FL, United States
- 10Center for Pharmacogenomics and Precision Medicine, College of Pharmacy, University of Florida, Gainesville, FL, United States
- 11Division of Cardiovascular Medicine, Department of Medicine and Center for Regenerative Medicine, University of Florida, Gainesville, FL, United States
- 12Novant Health Forsyth Medical Center, Hematology, Transplantation, and Cellular Therapy Division, Winston-Salem, NC, United States
Multiple myeloma (MM) is the second most frequent hematologic cancer in the United States. Carfilzomib (CFZ), an irreversible proteasome inhibitor being used to treat relapsed and refractory MM, has been associated with cardiotoxicity, including heart failure. We hypothesized that a multi-omics approach integrating data from different omics would provide insights into the mechanisms of CFZ-related cardiovascular adverse events (CVAEs). Plasma samples were collected from 13 MM patients treated with CFZ (including 7 with CVAEs and 6 with no CVAEs) at the University of Florida Health Cancer Center. These samples were evaluated in global metabolomic profiling, global proteomic profiling, and microRNA (miRNA) profiling. Integrative pathway analysis was performed to identify genes and pathways differentially expressed between patients with and without CVAEs. The proteomics analysis identified the up-regulation of lactate dehydrogenase B (LDHB) [fold change (FC) = 8.2, p = 0.01] in patients who experienced CVAEs. The metabolomics analysis identified lower plasma abundance of pyruvate (FC = 0.16, p = 0.0004) and higher abundance of lactate (FC = 2.4, p = 0.0001) in patients with CVAEs. Differential expression analysis of miRNAs profiling identified mir-146b to be up-regulatein (FC = 14, p = 0.046) in patients with CVAE. Pathway analysis suggested that the pyruvate fermentation to lactate pathway is associated with CFZ-CVAEs. In this pilot multi-omics integrative analysis, we observed the down-regulation of pyruvate and up-regulation of LDHB among patients who experienced CVAEs, suggesting the importance of the pyruvate oxidation pathway associated with mitochondrial dysfunction. Validation and further investigation in a larger independent cohort are warranted to better understand the mechanisms of CFZ-CVAEs.
Introduction
Multiple myeloma (MM) is a malignant plasma cell disorder and is the second most frequent hematologic cancer in the United States (1). Proteasome inhibitors (PIs) are among the most important classes of drugs to treat newly diagnosed, relapsed, and refractory MM (2). The mechanism of PI is dependent on the inhibition of proteasomal activity that disrupts the cell signaling pathways and leads to apoptosis and cell death (3, 4). So far, three PIs, bortezomib (Velcade®), carfilzomib (CFZ, Kyprolis®), and ixazomib (Ninlaro®), have been approved by the United States Food and Drug Administration (FDA) (5). CFZ is a second-generation PI that irreversibly inhibits 20S proteasome and is FDA-approved to treat relapsed or refractory MM patients who have received one to three previous treatments for MM. MM patients treated with CFZ and dexamethasone had improved survival compared with those treated with bortezomib and dexamethasone (6, 7).
Patients with MM often have cardiovascular diseases (CVDs) that may be a result of factors unrelated to MM, including older age, diabetes, obesity, or factors related to MM (such as amyloidosis), and other factors associated with the treatment of MM (such as anthracyclines, corticosteroids, alkylating agents, and PIs) (8). In an interim analysis of the phase 3 ENDEAVOR study, CFZ resulted in a significant and clinically meaningful survival benefit in relapsed MM patients compared with bortezomib, but it also was associated with significantly higher grade 3 cardiotoxicity, including heart failure (HF). In two meta-analyses, CFZ was associated with an increased incidence of cardiovascular adverse events (CVAEs) (8–18%) (2, 5) including hypertension, HF, arrhythmia, and cardiomyopathy (9, 10). CFZ accumulates in the heart, causes potent inhibition of the cardiac proteasome, and also decreases the 20S chymotrypsin-like activity to 10% in cardiomyocytes and 45% in bone marrow cells within 5 min after administration (11). Clinical trials and real-world experiences reported increasing incidence of CVAEs due to CFZ (CFZ-CVAE) in patients with MM compared with patients without MM with an adjusted hazard ratio (HR) ranging from 1.74 to 4.09 (8), which led to hospitalization, discontinuation of CFZ treatment, and also treatment-related death in ~5% of patients (12, 13). A critical clinical implication of this cardiotoxicity is treatment interruption that is associated with cancer recurrence (14). So far, there has been no human study investigating the mechanisms of CFZ-CVAE nor germline genetic risk factors for CFZ-CVAE. An integrative approach that combines multi-omics data, such as metabolomics, proteomics, transcriptomics, and microRNA (miRNA), may provide important insights into complex biological processes related to CFZ-CVAE.
This pilot study aimed to explore an integrative multi-omics approach to investigate data from different omics platforms including proteomics, metabolomics, and miRNAs to gain insights into the molecular mechanisms of CFZ-CVAE in MM patients.
Materials and Methods
Patients
The study was a retrospective cohort study of MM patients treated with CFZ at the University of Florida Health Cancer Center. Post-treatment whole blood samples were collected from MM patients who provided consent. CVAEs were evaluated in patients who had HF symptoms that led to obtaining the echocardiography and were confirmed by a cardiologist. A total of 13 MM patients were analyzed, including 7 who experienced CVAEs. This study was approved by the University of Florida Institutional Review Board (IRB) (IRB# 201702876).
Global Metabolomic Profiling
Global metabolomic profiling of plasma samples was performed on a Thermo Q-Exactive Orbitrap mass spectrometer coupled with ultra-high-performance liquid chromatography (UHPLC) at the Southeast Center for Integrated Metabolomics at the University of Florida. Details of the analytical method are described in the method section of the Supplementary Materials. All metabolites from positive and negative ion modes were separately subjected to quality control steps, data normalization, and statistical analyses using MetaboAnalyst 4.0 (15, 16). T-test and partial least square discrimination analysis (PLS-DA) were performed to identify metabolites with different abundance in those with and without CVAEs. Known metabolites with fold change (FC) of >1.5 and raw p-value of <0.001 were selected for pathway analysis using Ingenuity Pathway Analysis (IPA) (Qiagen IPA, USA).
Global Proteomics Profiling
Global proteomic profiling was performed on the plasma samples using the isobaric tags for relative and absolute quantification (iTRAQ) high-throughput proteomic method on these 13 MM patients treated with CFZ to identify proteomic biomarkers associated with CFZ-CVAEs. For quantification, only MS/MS spectra that were unique to a particular peptide and the sum of the signal-to-noise ratios for all the peak pairs >9 were used. The ratios with p < 0.05 present in at least two replicates were considered significant. Only the significant ratios from the replicates were used to calculate the average ratio for the protein. Details of the analytical method are described in the Method section of the Supplementary Materials.
Circulating miRNA Extraction and Profiling
Total RNA was extracted from the plasma samples using the miRVana microRNA Isolation Kit (Thermo Fisher Scientific, USA) and reverse transcribed using Megaplex™ RT Primers. Preamplification of cDNA with Megaplex™ PreAmp primers was performed according to the manufacturer's protocol and recommendations (Thermo Fisher Scientific, USA). miRNA expression was quantified using TaqMan® OpenArray® Human MicroRNA Panel containing 754 human miRNAs in QuantStudio™ 12K Flex Real-Time PCR.
Statistical Analysis
The characteristics of the patients were summarized with mean and standard deviation, median and range, or number and percentages as appropriate. The characteristics of patients with and without CVAEs were compared using Fisher's exact test for categorical variables and Mann–Whitney U test for continuous variables.
The high-throughput data generated from TaqMan® OpenArray® RT-qPCR, the miRNAs with Cq-values above 35, and the amplification score below 1.24 were excluded. The combination of miR-24 and miR-16 was used for normalization according to the NormFinder software suggestion with the most stable value of 0.005 (17). MiRNAs were determined to be significantly up- or down-regulated in MM samples with CVAEs and FC >1.5 or <0.5, p ≤ 0.05; the 2−ΔΔCT method was used for analysis (18). All statistical analyses were performed using SAS v.9.4 (Cary, NC, USA), Statistical Package for the Social Sciences SPSS v.26 (IBM Corp., USA), or R software.
Multi-Omics Integrative Analysis and Prediction of Biological Functions and Networks
The canonical pathways and networks associated with the differentially expressed metabolites, proteins, and miRNAs were identified using Qiagen's Ingenuity® Pathway Analysis version 6.0 (IPA® QIAGEN Bioinformatics, CA, USA; www.qiagen.com/ingenuity). IPA helps to predict the biological activity across thousands of datasets that represent disease conditions. By incorporating the top proteins, metabolites, and miRNAs, the canonical pathways associated with the multi-omics data can be identified. The pathways enriched with the top signals in any of the multi-omics data sources were identified using Fisher's exact test.
Results
A total of 13 MM patients treated with CFZ were analyzed, including 7 (53.8%) who experienced CVAEs. The overall median age was 58 years (range 46–76), and 38.5% were men. The characteristics of patients with or without CVAEs are summarized in Table 1. There were no statistical differences in demographics, comorbidities, or MM status between the patients who developed CVAEs and those who did not. There was no significant difference in baseline left ventricular ejection fraction (LVEF) between the two groups. However, there was a trend of lower LVEF in patients with CVAE (median 32.5) than in those without CVAEs (median of 57.5) (p = 0.08).
Metabolomic Profiling
A total of 8,636 metabolites were identified in the untargeted metabolomic profiling. After quality control, a total of 375 known metabolites were analyzed. The PLS-DA scores plots indicated a good separation of metabolites between patients with and without CVAEs (Figures 1A,B). Seventeen metabolites had FC >1.5 and raw p < 0.001. The top 25 metabolites that differentiated these two groups of patients in positive and negative ionization are illustrated in the heatmaps (Figures 1C,D). Among the top metabolites, lower intensities of lysophospholipids were identified in patients with CVAEs compared with patients without CVAEs. For example, the abundance of lysophosphatidylcholines [LysoPC (18:3), HMDB0010387] in patients with CVAEs was 33% of that in patients without CVAEs (p = 0.0003). The abundance of LysoPC (20:3, HMDB0010393) in patients with CVAEs was 39% of that in patients without CVAEs (p = 0.0005). The abundance of lysophosphatidylethanolamine [LysoPE (20:3), HMDB0011515] in patients with CVAEs was 44% of that as compared with patients without CVAEs (p = 0.0005) (Figures 2A–C).
Figure 1. Metabolomic profiling of MM patients treated with CFZ. Green color: patients without CVAE. Red color: patients with CVAE. (A) PLS-DA 2D scores plot of metabolites in the positive ion channel; (B) PLS-DA 2D scores plot of metabolites in the negative ion channel; (C) heatmap of the top 25 metabolites in the positive ion channel; (D) heatmap of the top 25 metabolites in the negative ion channel.
Figure 2. Figure showing the normalized abundance of metabolites with significantly different levels between patients with and without CVAEs after treatment with CFZ. (A) LysoPC 18:3, (B) LysoPC 20:3, (C) LysoPE 20:3, (D) pyruvic acid, (E) N-lactoyl methionine.
Additionally, the metabolomics data analysis identified lower abundance of pyruvic acid (HMDB0000243) in patients with CVAEs as 16% of that seen in patients without CVAEs (p = 0.0004), but an increased abundance of the lactate modified amino acid (N-lactoyl-methionine, HMDB0062182) (2.4-fold) in patients with CVAEs compared with the abundance in patients without CVAEs (p = 0.0001) (Figures 2D,E). All top metabolites that differentiated patients with and without CVAEs (with FC >1.5 or <0.5, p < 0.01) are listed in Table 2.
Table 2. Top metabolites that differentiate patients with and without CVAEs (with fold change >1.5 or <0.5 and p < 0.01).
Proteomic Profiling
From a total of 374 identified proteins, 347 were quantified with 202 excluded due to variation within the non-CVAE group. Fifty-seven proteins were differentially expressed, including 12 proteins down-regulated and 45 up-regulated in MM patients with CVAEs compared with those without CVAEs (FC >1.5 or <0.5, p < 0.05) (Table 3). The differential expression between the two groups is presented in the volcano plot (Figure 3). The top proteins that were up-regulated in patients with CVAEs compared with those without CVAEs include complement factor D (CFD) (FC = 6.15, p = 0.0006), lactate dehydrogenase B (LDHB) (FC = 8.24, p = 0.015), and insulin-like growth factor-binding protein 2 (IGFBP2) (FC = 7.67, p = 0.0017) (Table 3).
Table 3. The list of a total of 57 differentially expressed proteins at <0.5 or >1.5 of fold changes with p < 0.05.
Figure 3. Volcano plot of differentially expressed proteins between patients with and without CVAEs after treatment with CFZ. Negative log10 p-values were plotted on the y-axis, and log2 normalized fold changes expression levels on the x-axis. Significant differential expression with p-value threshold values (α = 0.05) was detected at <0.5 or >1.5-fold changes. Red and green dots indicate up- and down-regulated proteins in patients with CVAEs relative to patients without CVAEs.
MiRNA Profiling
TaqMan OpenArray analyzed a total of 750 miRNAs in 13 MM patients treated with CFZ. Twenty-one of these miRNAs were identified to be up-regulated or down-regulated in patients with CVAEs compared with patients without CVAEs (Table 4). Figure 4 shows the volcano plot of the OpenArray miRNAs data that were up-regulated or down-regulated in MM patients with CVAEs compared with those without CVAEs. Differential expression analysis of miRNAs profiling identified mir-146b to be up-regulated in those with CVAEs compared with those without CVAEs (FC = 14, p = 0.046).
Table 4. MicroRNA up-regulated or down-regulated in MM with or without CVAEs as determined by Open Array fold change <0.5 and ≥1.5.
Figure 4. Volcano plot of differential expression of miRNA between patients with and without CVAEs. MicroRNAs are differentially expressed in MM patients treated with CFZ. Volcano plot of significantly up-regulated (red dots) and down-regulated (green dots) microRNAs in patients with and without CVAEs as detected by TaqMan OpenArray profiling. Negative log10 p-values were plotted on the y-axis, and log2 normalized fold changes expression levels on the x-axis. A significant differential expression was detected with the p ≤ 0.05.
Pathway Analysis Using IPA
The top proteins and metabolites that were differentially expressed in the two groups were incorporated into the IPA analysis. The top two most significant canonical pathways associated with CFZ-CVAEs are the glutamate-dependent acid resistance pathway (p = 0.006) and pyruvate fermentation to lactate pathway (p = 0.01).
Discussion
In this pilot study, we performed integrative multi-omics analyses of metabolomic, proteomic, and miRNA data to identify differences among MM patients with or without CVAEs. Up-regulation of LDHB, lower plasma abundance of pyruvate, higher plasma abundance of lactate, and up-regulation of mir-146b were observed in patients with CVAEs compared with those without CVAEs. The most significant pathway identified was the glutamate-dependent acid resistance pathway that is related to fatty acid biotransformation activity. (19) The second top significant pathway, the pyruvate oxidation pathway, was related to mitochondrial dysfunction (20).
Mitochondria is the active organelle the produces most reactive oxygen species (ROS) that leads to peroxidation of lipids and DNA damage in cardiomyocytes (21–23). In cardiac cells, the LDHB converts lactate to pyruvate, which is considered the critical metabolite in glycolysis and glucose oxidation. The pyruvate oxidation (pyruvate decarboxylation) occurs after pyruvate is converted to acetyl-CoA in the mitochondria by the pyruvate dehydrogenase complex (PDC). Acetyl-CoA then is used in the citric acid cycle (CAC) for adenosine triphosphate (ATP) production, which is essential for cardiomyocyte contractility (24). The decarboxylation of pyruvate is dominant in the mitochondria of healthy cardiac cells (24). Mitochondrial protein complex (MPC) is a mitochondrial inner-protein complex that consists of two subunits: MPC1 and MPC2. Both subunits are essential for the normal function of MPC, and loss of any subunit would affect the mitochondrial pyruvate transport and the pyruvate supporting the oxygen consumption in cardiac cells (25). A recent study identified a correlation between LDHB and cardiac hypertrophy and suggested a molecular mechanism for regulating abnormal energy metabolism in cardiac hypertrophy (26–28). Based on our preliminary results, there was compelling evidence for the pyruvate oxidation pathway's involvement in CFZ-CVAEs.
A previous study of mice published by Efentakis and colleagues suggested that CFZ induces cardiotoxicity through increased protein phosphatase 2 (PP2A) activity and inhibition of AMP-activated protein kinase (AMPK) and its downstream autophagic targets (4). Our untargeted metabolomic analyses identified significantly lower abundance of phospholipids, such as LysoPC. In another study in non-cancer patients with Heart Failure Reduced Ejection Fraction (HFrEF), the authors also found decreasing phospholipids compared with healthy controls matched by age, sex, body mass index, and ischemic heart disease status (29). Lower lysophospholipid levels coupled with higher long-chain acylcarnitine levels may reflect a shift toward increased myocardial glucose oxidation, with down-regulation of fatty acid oxidation as the result of mitochondrial dysfunction in HF patients. Lower lysophospholipid levels result from increased PP2A activity, which is an essential and ubiquitously expressed serine–threonine phosphatase, plays a critical role in cellular processes, such as cell proliferation and apoptosis, and is responsible for more than 90% of protein dephosphorylation in cardiac cells (30, 31). Studies have shown that lower levels of lysophospholipids are likely to increase PP2A activity (32) and inhibition of AMPK (33). Our findings of metabolomic readout of MM patients treated with CFZ were consistent with a previously published study in mice.
We also discovered the CFD to be significantly up-regulated in patients with CVAEs. CFD is targeted by miR-17, miR-19b, miR-92a, miR-30b,c, and miR-106a (34). CFD also promotes the alternative pathway complement (AP) activation, and the complement activation was shown to occur in chronic HF and is associated with adverse clinical events (35). These findings suggest that CFZ may induce cardiac dysfunction via increasing CFD activity, leading to the complement system's activation. We also observed an elevation of IGFBP2 in patients who experienced CVAEs. Previous studies reported that IGFBP2 might play an essential role as a potential diagnostic and prognostic marker for suspected HF (30, 36, 37). However, there is still a lack of information about the structure and function of IGFBPs associated with CVD.
In our study, we also detected the up-regulation of miR-146b in patients with CVAEs. In a study performed on rats, the authors reported the role of miR-146b during myocardial ischemia/reperfusion injury (IRI) in myocardial I/R rat and hypoxia/reoxygenation (H/R) injured cardiomyocytes. They found that miR-146b overexpression significantly decreased the myocardial infarct size and cardiomyocytes apoptotic rates and release of creatine kinase and lactate dehydrogenase during IRI by targeting Smad4 and modulating the expression of c-fos and c-JUN in the myocardium (31). Another study showed that miR-146b inhibition increased hypoxia-induced cardiomyocyte apoptosis, which indicated the potential protective effect of miR-146b in ischemia injury (38). However, the role of miR-146b in myocardial ischemic injury is still unclear (4).
A healthy heart relies primarily on fatty acid oxidation for ATP production, with a smaller contribution from glucose oxidation. In early HF, glucose oxidation increases, and glucose is converted to pyruvate in the cytosol, then transported into the mitochondria, and converted to acetyl Co-A by pyruvate dehydrogenase for ATP production and to lactate by LDHB. Both fatty acid and glucose oxidation are decreased in advanced HF (39–41). In this cohort study, proteomics analysis identified the up-regulation of LDHB in patients who experienced CVAEs. On the other hand, the metabolomics data analysis showed lower level of pyruvate but increased level of lactate in patients with CVAEs. A recent study identified a correlation between LDHB and cardiac hypertrophy and suggested a molecular mechanism for regulating abnormal energy metabolism in cardiac hypertrophy (26). LDHB may be an appropriate target to improve abnormal cardiac energy metabolism in the setting of cardiac hypertrophy.
Based on the integrative analysis of multi-omics data, we identified the association between the pyruvate fermentation to lactate (pyruvate oxidation) pathway and CFZ-CVAEs. Our results highlight the increase in lactate level and decrease of pyruvate level after CFZ administration, which supports the idea that anaerobic glycolysis occurring in myocardial ischemia leads to mitochondrial metabolic dysfunction and reduction in ATP production through oxidative phosphorylation (42, 43). We propose a possible mechanism for the CFZ-CVAEs in Figure 5. A lower level of pyruvate and a higher level of lactate in patients with CVAEs are consistent with cardiomyocytes switching substrate preference from fatty acids to glucose because of mitochondrial dysfunction.
Figure 5. Proposed molecular mechanism of CFZ-CVAEs. Red arrows indicate the changes in this study. CFZ, carfilzomib; LDHB, lactate dehydrogenase B; PDH, pyruvate dehydrogenase; NADH, nicotinamide adenine dinucleotide hydrogen; ATP, adenosine triphosphate.
It is important to acknowledge a couple of limitations of our study. First, this was a pilot study with a small sample size selected from a single center. Therefore, the metabolomic, proteomic, and miRNA biomarkers that differentiate MM patients with and without CVAEs need to be further validated before they can be considered useful clinically. Nevertheless, these preliminary findings point to a plausible pathway that indicates the CFZ-CVAEs might be related to mitochondrial dysfunction. Second, since the plasma samples were obtained after the patients were treated with CFZ, the identified biomarkers have no predictive value. In order to identify clinically useful predictive biomarkers, pre-treatment samples need to be analyzed. However, the differences between patients with and without CVAEs provided essential insights into the mechanisms of CFZ-CVAEs.
In summary, in this pilot multi-omics integrative analysis, we observed the down-regulation of pyruvate and up-regulation of LDHB among patients with CFZ-CVAEs, suggesting the importance of pyruvate oxidation pathway and mitochondrial dysfunction in CFZ-CVAEs. Further investigation in a larger independent cohort is warranted to better understand the mechanisms of CFZ-CVAEs.
Data Availability Statement
The data presented in the study are deposited in the ProteomeXchange Consortium via the PRIDE partner repository, accession number (PXD024281).
Ethics Statement
The studies involving human participants were reviewed and approved by University of Florida Institutional Review Board. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
MT wrote the manuscript. MT, LC, YH, and YG performed the statistical analyses. TG performed the metabolomics analysis. SS assisted in the pathway analysis. CS, NW, and JM collected the patient samples. YG, RC, RB, MF, DD, LY, TL, KM, and CP provided the critical revision of the manuscript. YG secured the funding for this study. All authors contributed to the article and approved the submitted version.
Funding
This research was supported in part by the University of Florida Health Cancer Center and in part by the University of Florida Clinical and Translational Science Institute, which is supported in part by the NIH National Center for Advancing Translational Sciences under award number UL1TR001427.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.645122/full#supplementary-material
References
1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019 (US statistics). CA Cancer J Clin. (2019) 69:7–34. doi: 10.3322/caac.21551
2. Gandolfi S, Laubach JP, Hideshima T, Chauhan D, Anderson KC, Richardson PG. The proteasome and proteasome inhibitors in multiple myeloma. Cancer Meta Rev. (2017) 36:561–84. doi: 10.1007/s10555-017-9707-8
3. Hasinoff BB, Patel D, Wu X. Molecular mechanisms of the cardiotoxicity of the proteasomal-targeted drugs bortezomib and carfilzomib. Cardiovasc Toxicol. (2017) 17:237–50. doi: 10.1007/s12012-016-9378-7
4. Efentakis P, Kremastiotis G, Varela A, Nikolaou PE, Papanagnou ED, Davos CH, et al. Molecular mechanisms of carfilzomib-induced cardiotoxicity in mice and the emerging cardioprotective role of metformin. Blood. (2019) 133:710–23. doi: 10.1182/blood-2018-06-858415
5. Teicher BA, Tomaszewski JE. Competitive landscape report. Biochem Pharmacol. (2015) 96:1–9. doi: 10.1016/j.bcp.2015.04.008
6. Dimopoulos MA, Goldschmidt H, Niesvizky R, Joshua D, Chng WJ, Oriol A, et al. Carfilzomib or bortezomib in relapsed or refractory multiple myeloma (ENDEAVOR): an interim overall survival analysis of an open-label, randomised, phase 3 trial. Lancet Oncol. (2017) 18:1327–37. doi: 10.1016/S1470-2045(17)30578-8
7. Libby E, Garcia D, Quintana D, Fekrazad MH, Bauman J, Ebaid A, et al. Disease-specific survival for patients with multiple myeloma: significant improvements over time in all age groups. Leuk Lymphoma. (2014) 55:2850–7. doi: 10.3109/10428194.2014.897700
8. Kistler KD, Kalman J, Sahni G, Murphy B, Werther W, Rajangam K, et al. Incidence and risk of cardiac events in patients with previously treated multiple myeloma versus matched patients without multiple myeloma: an observational, retrospective, cohort study. Clin Lymphoma, Myeloma Leuk. (2017) 17:89–96.e3. doi: 10.1016/j.clml.2016.11.009
9. Atrash S, Tullos A, Panozzo S, Bhutani M, Van Rhee F, Barlogie B, et al. Cardiac complications in relapsed and refractory multiple myeloma patients treated with carfilzomib. Blood Cancer J. (2015) 5:e272. doi: 10.1038/bcj.2014.93
10. Siegel D, Martin T, Nooka A, Harvey RD, Vij R, Niesvizky R, et al. Integrated safety profile of single-agent carfilzomib: experience from 526 patients enrolled in 4 phase II clinical studies. Haematologica. (2013) 98:1753–61. doi: 10.3324/haematol.2013.089334
11. Heckmann MB, Doroudgar S, Katus HA, Lehmann LH. Cardiovascular adverse events in multiple myeloma patients. J Thorac Dis. (2018) 10:S4296–305. doi: 10.21037/jtd.2018.09.87
12. Muchtar E, Gatt ME, Rouvio O, Ganzel C, Chubar E, Suriu C, et al. Efficacy and safety of salvage therapy using carfilzomib for relapsed or refractory multiple myeloma patients: a multicentre retrospective observational study. Br J Haematol. (2016) 172:89–96. doi: 10.1111/bjh.13799
13. Danhof S, Schreder M, Rasche L, Strifler S, Einsele H, Knop S. ‘Real-life’ experience of preapproval carfilzomib-based therapy in myeloma - analysis of cardiac toxicity and predisposing factors. Eur J Haematol. (2016) 97:25–32. doi: 10.1111/ejh.12677
14. Yu AF, Yadav NU, Lung BY, Eaton AA, Thaler HT, Hudis CA, et al. Trastuzumab interruption and treatment-induced cardiotoxicity in early HER2-positive breast cancer. Breast Cancer Res Treat. (2015) 149:489–95. doi: 10.1007/s10549-014-3253-7
15. Chong J, Wishart DS, Xia J. Using metaboanalyst 4.0 for comprehensive and integrative metabolomics data analysis. Curr Protoc Bioinforma. (2019) 68:e86. doi: 10.1002/cpbi.86
16. Xia J, Psychogios N, Young N WD. MetaboAnalyst: a web server for metabolomic data analysis and interpretation. Nucleic Acids Res. (2009) W652–60. doi: 10.1093/nar/gkp356
17. Andersen CL, Jensen JL, Ørntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. (2004) 64:5245–50. doi: 10.1158/0008-5472.CAN-04-0496
18. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. (2001) 25:402–8. doi: 10.1006/meth.2001.1262
19. Woo JM, Kim JW, Song JW, Blank L1M, Park JB. Activation of the glutamic acid-dependent acid resistance system in escherichia coli BL21(DE3) leads to increase of the fatty acid biotransformation activity. PLoS ONE. (2016) 11:e0163265. doi: 10.1371/journal.pone.0163265
20. Brooks GA, Dubouchaud H, Brown M, Sicurello JP, Eric Butz C. Role of mitochondrial lactate dehydrogenase and lactate oxidation in the intracellular lactate shuttle. Proc Natl Acad Sci USA. (1999) 96:1129–34. doi: 10.1073/pnas.96.3.1129
21. Gewirtz DA. A critical evaluation of the mechanisms of action proposed for the antitumor effects of the anthracycline antibiotics adriamycin and daunorubicin. Biochem Pharmacol. (1999) 57:727–41. doi: 10.1016/S0006-2952(98)00307-4
22. Angsutararux P, Luanpitpong S, Issaragrisil S. Chemotherapy-induced cardiotoxicity : overview of the roles of oxidative stress. (2015) 2015:795602. doi: 10.1155/2015/795602
23. Lipshultz SE, Alvarez JA, Scully RE. Anthracycline associated cardiotoxicity in survivors of childhood cancer. Heart. (2007) 94:525–33. doi: 10.1136/hrt.2007.136093
24. Comte B, Vincent G, Bouchard B, Des Rosiers C. Probing the origin of acetyl-CoA and oxaloacetate entering the citric acid cycle from the 13C labeling of citrate released by perfused rat hearts. J Biol Chem. (1997) 272:26117–24. doi: 10.1074/jbc.272.42.26117
25. Herzig S, Raemy E, Montessuit S, Veuthey JL, Zamboni N, Westermann B, et al. Identification and functional expression of the mitochondrial pyruvate carrier. Science. (2012) 336:93–6. doi: 10.1126/science.1218530
26. Feng H, Wu J, Chen P, Wang J, Deng Y, Zhu G, et al. MicroRNA-375-3p inhibitor suppresses angiotensin II-induced cardiomyocyte hypertrophy by promoting lactate dehydrogenase B expression. J Cell Physiol. (2019) 234:14198–209. doi: 10.1002/jcp.28116
27. Frey N, Katus HA, Olson EN, Hill JA. Hypertrophy of the heart: a new therapeutic target? Circulation. (2004) 109:1580–9. doi: 10.1161/01.CIR.0000120390.68287.BB
28. Shimizu I, Minamino T. Physiological and pathological cardiac hypertrophy [Internet]. J Mol Cell Cardiol. (2016) 97:245–62. doi: 10.1016/j.yjmcc.2016.06.001
29. Marcinkiewicz-Siemion M, Ciborowski M, Ptaszynska-Kopczynska K, Szpakowicz A, Lisowska A, Jasiewicz M, et al. LC–MS-based serum fingerprinting reveals significant dysregulation of phospholipids in chronic heart failure. J Pharm Biomed Anal. (2018) 154:354–63. doi: 10.1016/j.jpba.2018.03.027
30. Hoeflich A, David R, Hjortebjerg R. Current IGFBP-related biomarker research in cardiovascular disease-We need more structural and functional information in clinical studies. Front Endocrinol. (2018) 9:388. doi: 10.3389/fendo.2018.00388
31. Di YF, Li DC, Shen YQ, Wang CL, Zhang DY, Shang AQ, et al. MiR-146b protects cardiomyocytes injury in myocardial ischemia/reperfusion by targeting Smad4. Am J Transl Res. (2017) 9:656–63.
32. Huynh K, Barlow CK, Jayawardana KS, Weir JM, Mellett NA, Cinel M, et al. High-Throughput plasma lipidomics: detailed mapping of the associations with cardiometabolic risk factors. Cell Chem Biol. (2019) 26:71–84.e4. doi: 10.1016/j.chembiol.2018.10.008
33. Zhu Y, Feng Y, Shen L, Xu D, Wang B, Ruan K, et al. Effect of metformin on the urinary metabolites of diet-induced-obese mice studied by ultra performance liquid chromatography coupled to time-of-flight mass spectrometry (UPLC-TOF/MS). J Chromatogr B Anal Technol Biomed Life Sci. (2013) 925:110–6. doi: 10.1016/j.jchromb.2013.02.040
34. Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. (2015) 4:1–38. doi: 10.7554/eLife.05005.028
35. Shahini N, Michelsen AE, Nilsson PH, Ekholt K, Gullestad L, Broch K, et al. The alternative complement pathway is dysregulated in patients with chronic heart failure. Sci Rep. (2017) 7:1–10. doi: 10.1038/srep42532
36. Berry M, Galinier M, Delmas C, Fournier P, Desmoulin F, Turkieh A, et al. Proteomics analysis reveals IGFBP2 as a candidate diagnostic biomarker for heart failure. IJC Metab Endocr. (2015) 6:5–12. doi: 10.1016/j.ijcme.2014.11.003
37. Barutaut M, Fournier P, Peacock WF, Evaristi MF, Caubère C, Turkieh A, et al. Insulin-like growth factor binding protein 2 predicts mortality risk in heart failure. Int J Cardiol. (2020) 300:245–51. doi: 10.1016/j.ijcard.2019.09.032
38. Li JW, He SY, Feng ZZ, Zhao L, Jia WK, Liu P, et al. MicroRNA-146b inhibition augments hypoxia-induced cardiomyocyte apoptosis. Mol Med Rep. (2015) 12:6903–10. doi: 10.3892/mmr.2015.4333
39. Fillmore N, Mori J, Lopaschuk GD. Mitochondrial fatty acid oxidation alterations in heart failure, ischaemic heart disease and diabetic cardiomyopathy. Br J Pharmacol. (2014) 171:2080–90. doi: 10.1111/bph.12475
40. Lionetti V, Stanley WC, Recchia FA. Modulating fatty acid oxidation in heart failure. Cardiovasc Res. (2011) 90:202–9. doi: 10.1093/cvr/cvr038
41. Li X, Liu J, Lu Q, Ren D, Sun X, Rousselle T, et al. AMPK: A therapeutic target of heart failure—not only metabolism regulation. Biosci Rep. (2019) 39:1–13. doi: 10.1042/BSR20181767
42. Barry WH. Heart physiology from cell to circulation, 4th ed. Circulation. (2004) 110:e313. doi: 10.1161/01.CIR.0000143724.99618.62
Keywords: proteasome inhibitors, Cardio-oncology, carfilzomib, metabolomcis, proteomic
Citation: Tantawy M, Chekka LM, Huang Y, Garrett TJ, Singh S, Shah CP, Cornell RF, Baz RC, Fradley MG, Waheed N, DeRemer DL, Yuan L, Langaee T, March K, Pepine CJ, Moreb JS and Gong Y (2021) Lactate Dehydrogenase B and Pyruvate Oxidation Pathway Associated With Carfilzomib-Related Cardiotoxicity in Multiple Myeloma Patients: Result of a Multi-Omics Integrative Analysis. Front. Cardiovasc. Med. 8:645122. doi: 10.3389/fcvm.2021.645122
Received: 22 December 2020; Accepted: 02 March 2021;
Published: 29 April 2021.
Edited by:
Carlo Gabriele Tocchetti, University of Naples Federico II, ItalyReviewed by:
Christian Cadeddu Dessalvi, University of Cagliari, ItalyAlessandra Cuomo, Federico II University Hospital, Italy
Copyright © 2021 Tantawy, Chekka, Huang, Garrett, Singh, Shah, Cornell, Baz, Fradley, Waheed, DeRemer, Yuan, Langaee, March, Pepine, Moreb and Gong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yan Gong, gong@cop.ufl.edu