Inhibition of thioredoxin activates mitophagy and overcomes adaptive bortezomib resistance in multiple myeloma
Journal of Hematology & Oncology volume 11, Article number: 29 (2018)
Although current chemotherapy using bortezomib (Velcade) against multiple myeloma in adults has achieved significant responses and even remission, a majority of patients will develop acquired resistance to bortezomib. Increased thioredoxin level has been reported to be associated with carcinogenesis; however, the role of thioredoxin in bortezomib drug resistance of myeloma remains unclear.
We generated several bortezomib-resistant myeloma cell lines by serially passaging with increased concentrations of bortezomib over a period of 1.5 years. Thioredoxin expression was measured by real-time PCR and western blot.
The role of thioredoxin in the survival of bortezomib-resistant myeloma cells was determined by specific shRNA knockdown in vitro and in vivo. Thioredoxin inhibitor (PX12) was used to determine the effectiveness of thioredoxin inhibition in the treatment of bortezomib-resistant myeloma cells. The effect of thioredoxin inhibition on mitophagy induction was examined. The correlation of thioredoxin expression with patient overall survival was interrogated. Thioredoxin expression was significantly upregulated in bortezomib-resistant myeloma cells and the change correlated with the increase of bortezomib concentration. Thioredoxin gene knockdown using specific shRNA sensitized bortezomib-resistant myeloma cells to bortezomib efficiency in vitro and in vivo. Similarly, pharmacological inhibition with PX12 inhibited the growth of bortezomib-resistant myeloma cells and overcame bortezomib resistance in vitro and in vivo. Furthermore, inhibition of thioredoxin resulted in the activation of mitophagy and blockage of mitophagy prevented the effects of PX12 on bortezomib-resistant myeloma cells, indicating that mitophagy is the important molecular mechanism for the induction of cell death in bortezomib-resistant myeloma cells by PX12. Moreover, inhibition of thioredoxin resulted in downregulation of phosphorylated mTOR and ERK1/2. Finally, thioredoxin was overexpressed in primary myeloma cells isolated from bortezomib-resistant myeloma patients and overexpression of thioredoxin correlated with poor overall survival in patients with multiple myeloma.
Our findings demonstrated that increased thioredoxin plays a critical role in bortezomib resistance in multiple myeloma through mitophagy inactivation and increased mTOR and ERK1/2 phosphorylation. Thioredoxin provides a potential target for clinical therapeutics against multiple myeloma, particularly for bortezomib-resistant/refractory myeloma patients.
Multiple myeloma, also referred as plasma cell myeloma, plasmacytic myeloma, myelomatosis, or Kahler disease, is a neoplastic malignancy characterized by the proliferation of abnormal plasma cells derived from B cells [1, 2]. These plasma cells proliferate in the bone marrow and frequently invade into the adjacent bone, causing skeletal destruction that finally results in fractures [3, 4]. In addition, multiple myeloma can cause anemia, renal insufficiency, frequent infection, and hypercalcemia and is associated with significant mortality and morbidity. Multiple myeloma accounts for about 1% of all types of malignancy and 10% of hematologic malignancies. In 2017, there were 30,280 new diagnoses and an estimated 12,590 deaths from myeloma in the USA alone. Existing therapies for multiple myeloma, including the proteasome inhibitors bortezomib (BTZ) and carfilzomib, have the potential of extending the overall survival but are not curative . Bortezomib is one of the major drugs in the treatment of multiple myeloma . Almost all myeloma patients will receive bortezomib during their course of treatment, yet nearly all patients will develop drug resistance to bortezomib . Thus, to unmask the failure to cure, drug resistance mechanisms associated with bortezomib must be characterized.
Thioredoxin (TXN) is a 12 kDa ubiquitous oxidoreductase. Along with NADPH and thioredoxin reductase, thioredoxin constitutes the thioredoxin system, one of the major disulfide reductase antioxidant systems in mammalian cells. Thioredoxin is crucial in defense against oxidative stress and in maintaining the redox environment in the cell. In addition, TXN has multi-faceted roles . TXN regulates a variety of redox-sensitive signaling pathways as well as ROS-independent genes and exerts cytoprotective effects. Overexpression of thioredoxin has been implicated in the pathogenesis of advanced malignancies, including solid cancer and hematological malignancy . Besides, there are increasing evidences of its role in the development of resistance to several chemotherapeutic agents including cisplatin and docetaxel [10, 11].
Mitophagy, mitochondrial degradation by autophagic delivery to lysosomes, is the major degradative pathway in mitochondrial turnover. Mitochondria are the essential site of aerobic energy production in eukaryotic cells. Maintaining a healthy population of mitochondria is essential to the cellular redox homeostasis and the well-being of cell self-renewal . The role of mitophagy dysfunction in cancer pathogenesis is currently an area of active investigation but the findings varied . For instance, some studies suggested that a decrease in mitophagy led to the increase in the production of free radicals and subsequent genetic instability [14, 15], thus, favoring the development of cancer. On the other hand, studies have found that mitophagy protected cancer cells from apoptosis [16, 17]; thus, an increase in mitophagy will promote cancer cell survival and progression. It is highly likely that the double-edged roles of mitophagy dysfunction in cancer pathogenesis may change significantly depending on cancer cell types . However, the role of mitophagy in chemotherapeutics-induced drug resistance is still unknown and remains to be investigated.
We hypothesized that thioredoxin plays an important role in bortezomib drug resistance in multiple myeloma, providing a novel therapeutic target against cancer drug resistance in the treatment of multiple myeloma. Additionally, we aimed at investigating the effects of thioredoxin on mitophagy in bortezomib-resistant myeloma cells.
Cell culture and generation of bortezomib-resistant myeloma cells
The human multiple myeloma cell lines MM.1S, MM.1R, OPM1, RPMI8226/Dox, and NCIH929 were purchased from ATCC Company. Cells were grown in suspension in RMPI1640 medium supplemented with 10% fetal bovine serum, 1% (v/v) penicillin, and 100 μg/mL streptomycin. Cells were maintained at 37 °C in a 5% CO2 atmosphere with a proper humidity.
To generate bortezomib-resistant myeloma cell lines, bortezomib (BTZ) was added to the multiple myeloma cell culture medium starting at 0.03 nM. The culture medium was replaced with bortezomib containing medium twice weekly and the myeloma cells were cultured for 2–4 weeks until the cells survived and became resistant to that concentration of bortezomib. The bortezomib concentration was then increased by doubling the previous concentration until BTZ reached at a final concentration of up to 8.4 nM by the end of the second year. Myeloma cell lines were maintained on BTZ containing medium until 1 week before experiments.
Drugs and reagents
The proteasome inhibitor bortezomib (PS-341, Velcade®) was obtained from Selleckchem Chemicals LLC (Houston, TX). The thioredoxin inhibitor 2-[(1-methylpropyl) dithio]-1H-imidazole (PX12) was purchased from Tocris Bioscience (Bristol, UK). The mitophagy inhibitor bafilomycin was obtained from Sigma-Aldrich (St. Louis, Missouri, USA). 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazoliumbromide (MTT) were purchased from GIBCO BRL (Grand Island, NY). The mTOR activator (MHY1485) and ERK activator (tert-butylhydroquinone; tBHQ) were purchased from Sigma-Aldrich (St. Louis, Missouri, USA). Drugs were dissolved in dimethyl sulfoxide (DMSO; Sigma-Aldrich, St. Louis, Missouri, USA) at a concentration of 100 mM as a stock solution. Anti-PINK1 antibody (ab75487) and anti-LC3B antibody (ab51520) used for western blot analysis were purchased from Abcam (Cambridge, MA). Anti-beta-actin antibody (A2228) was obtained from Sigma-Aldrich. Anti-thioredoxin antibody (C63C6), anti-AKT antibody (9272), anti-mTOR antibody (2972), anti-phospho-Akt (Ser473, 4060) antibody, anti-phospho-mTOR (Ser2448, 5536) antibody, and anti-phospho-p44/42 MAPK (ERK1/2, Thr202/Tyr204, 9101) antibody were purchased from Cell Signaling Technology Inc. (Beverly, MA, USA). Human CD138 enrichment kit (EasySep™) was purchased from StemCell Technologies (Vancouver, BC, Canada).
Equal number of parental and BTZ-resistant multiple myeloma cells were seeded in 96-well plates. Different concentrations of BTZ with or without PX12 were added to the cells in each group (at least four replicates for each group). Cells were incubated in RPMI1640 medium at 37 °C in a 5% CO2 humidified atmosphere. After incubating for 48 or 72 h, 10 μL MTT (5 mg/mL) was added to each well and plates were incubated at 37 °C for another 4 h. Finally, 100 μL of 10% sodium dodecyl sulfate (with 0.01 N HCl) was added to dissolve the crystals and absorbance was determined at 570 nm in an EL340 microplate reader (BioTek Instruments, Winooski, VT). Ratios of the 50% inhibitory concentration (IC50) value in BTZ-resistant groups to the IC50 value of parental cells were calculated and considered to be the relative indicators of drug resistance in the experimental groups.
Gene expression analyses
For TXN mRNA expression (forward: 5′-GTAGTTGACTTCTCAGCCACGTG-3′, reverse: 5′-CTGACAGTCATCCACATCTACTTC-3′), total RNAs were extracted using TRIzol reagent (Invitrogen) according to standard procedures and reverse transcribed into complementary DNA (cDNA) using a BIO-RAD iScript ™ cDNA synthesis kit (Bio-Rad, Hercules, CA, USA). Samples were then analyzed using an Applied Biosystems Real-Time PCR (SYBR Green, Bio-Rad Laboratories, Hercules, CA, USA) in triplicate. Gene expression was normalized using 18S rRNA (forward: 5′-GTAACCCGTTGAACCCCATT-3′, reverse: 5′-CCATCCAATCGGTAGTAGCG-3′).
Western blot analysis
Western blot analysis was performed as previously described . Briefly, total protein was extracted using a RIPA buffer (50 mM Tris HCl, pH 7.4/150 mM NaCl/5 mM EDTA/1% NP-40/1% sodium deoxycholate/0.1% SDS/1% aprotinin, 50 mM NaF/0.1 mM Na3VO4), and equal amounts of proteins were separated using a SDS-PAGE electrophoresis. Separated proteins were then transferred to polyvinylidene difluoride membranes (PVDF; Millipore Corp., Bedford, MA, USA) and incubated with primary antibodies for thioredoxin (1:1000), PINK1 (1:200), LC3B (1:1000), AKT/pAKT (1:1000), mTOR/p-mTOR (1:1000), pERK1/2 (1:1000), or β-actin (1:10,000) overnight at 4 °C with gentle agitation. Membranes were washed and then incubated with HRP-conjugated secondary antibodies (1:10,000) for 2 h at room temperature before signal detection by chemiluminescence (Pierce Biotechnology, Rockford, IL, USA). Densitometric quantification was performed by Image-Pro Plus 6.0 software (Media Cybernetics, Silver Springs, MD 20910, USA).
Thioredoxin-specific shRNA knockdown
Plasmids targeting human thioredoxin (shTXN1-4, catalog number RHS4430-200171579, RHS4430-200174379, RHS4430-200273352, and RHS4430-200274993) were purchased from GE Healthcare (Piscataway, NJ, USA). Plasmid for non-targeting control (shNT) and the packing and envelope vectors psPAX2 and VSVG were obtained from Addgene (Cambridge, Massachusetts). HEK293T cells were transfected with shNT or shTXN, psPAX2, and VSV-G using Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) for 24 h. The DMEM medium was changed and collected at 24 and 48 h after transfection, respectively. To obtain thioredoxin stably knockdown cells, the transduced cells were cultured with 1 μg/ml puromycin and the GFP+ cells were sorted and expanded.
Mitochondrial network by transmission electron microscopy
Conventional transmission electron microscopy analysis was performed as described previously . Briefly, human multiple myeloma cells with or without treatment were fixed by a solution containing 4% formaldehyde and 2% glutaraldehyde. Specimens were washed following by OSO4 postfixed, alcohol dehydrated, and araldite embedded. Thin sections of samples were analyzed using a FEI Tecnai G2 Twin electron microscope (FEI, Hillsboro, OR, USA).
Determination of mitochondrial membrane potential (Δψm)
JC-1 fluorescent probe kit (Molecular Probes, Eugene, OR, USA) was used to determine Δψm with two different staining spectra, the orange aggregates dye form for normally respiring cells and green monomers for cells with respiratory dysfunction (apoptotic cells). Briefly, cells with or without treatment were incubated in RPMI1640 media containing JC-1 (2 μM final concentration) at 37 °C for 15 min. Cell pellets were resuspended in cold PBS and analyzed on a flow cytometer with 488-nm excitation emission.
In vivo tumor xenograft model
All animal experiments were approved by the Animal Care Committee of Duke University Medical Center. NOD/LtSz-scid/scid (NOD/SCID) mice were purchased from Jackson Laboratories (Bar Harbor, ME, USA) and maintained in microisolator cages on laminar flow racks under pathogen-free conditions in the Division of Laboratory Animal Resources, Duke University. BTZ-resistant MM.1R multiple myeloma cells, 3 × 106 cells in 100 μL PBS, were injected subcutaneously into the dorsal lateral flank of NOD/SCID mice that had received a total body irradiation with 1.5 Gy from a 137Cs source. Engraftment of human myeloma was monitored twice per week. When the tumor was established and palpable after 10 days of xenograft, the mice were given control buffer, PX12 (12 mg/kg, i.p., twice weekly), bortezomib (0.5 mg/kg, i.p., twice weekly) or the combination of PX12, and bortezomib. Tumor length and width were measured with a caliper and tumor volume was calculated using the formula V = (L × W × W)/2, where V is the tumor volume, L is the tumor length, and W is the tumor width. At the end of the experiments, tumors were harvested and weighed. A portion of the tumors were homogenized and used for western blot analysis.
Microarray data mining
For oncomining gene expression analysis, mRNA level of TXN in normal plasma cell and myeloma cells was queried using Oncomine database (https://www.oncomine.org). The prognostic value of TXN in multiple myeloma was assessed using Mulligan and Arkansas myeloma microarrays. Overall survival was compared between high and low expression of TXN using median gene expression value as a bifurcating point.
Bone marrow aspirates were obtained from patients with myeloma with patients’ informed consent and IRB approval. The study was performed in compliance with the guidelines of the Ethical Committee of Duke University Medical Center. A total of 13 patients with myeloma, including three newly diagnosed myeloma, five relapsed myeloma, and five treated myeloma prior to stem cell transplant were enrolled in the study. Human CD138+ and CD138− cells were isolated from the bone marrow aspirates of these patients using Histopaque-1077 (Sigma) gradient separation followed by human CD138 enrichment kit (EasySep™, StemCell Technologies). The purity of the human CD138+ cells was > 95%.
All statistical analyses were performed using the Student’s t test and represented as mean ± standard error of the mean (SEM) unless noted otherwise. For in vivo experiments with ≥ 5 per group, statistical analyses were performed with two-way ANOVA followed by multiple comparisons. The p values were designated as *p < 0.05, **p < 0.01, ***p < 0.005, ****p < 0.001; n.s. non-significant (p > 0.05).
Upregulation of thioredoxin correlates with the development of bortezomib drug resistance in multiple myeloma cells
To generate adaptive BTZ-resistant multiple myeloma cell lines, we cultured five different myeloma cell lines (MM.1S, MM.1R, OPM1, RPMI8226/Dox, and NCIH929) over a period of 1.5 years with gradually increased concentrations of BTZ. At the end of the serial cultures, the IC50 of the BTZ-resistant myeloma cell lines to BZT has increased by 10–100 fold (ranging from 7.06 to 110 fold) (Fig. 1a, b and Additional file 1: Fig. S1A–1C). Compared to their parental BTZ-sensitive cells, BTZ-resistant myeloma cells demonstrated a significant upregulation of both thioredoxin mRNA gene expression and protein expression (Fig. 1c–f and Additional file 1: Fig. S1D–1I). Importantly, thioredoxin protein expression increased in parallel to the increase in drug resistance to BTZ (Fig. 1e, f). These data suggested that thioredoxin plays an important role in the development of bortezomib drug resistance in multiple myeloma.
Genetic knockdown of thioredoxin sensitizes bortezomib-resistant myeloma cells to bortezomib
We next determined the role of thioredoxin in bortezomib drug resistance using shRNA-specific gene knockdown techniques. Four shRNAs specific for human thioredoxin were transfected into MM.1R-BTZ-resistant myeloma cells and RPMI8226/Dox-BTZ-resistant myeloma cells. Thioredoxin knockdown efficiencies were determined by western blot. As shown in Fig. 2a, d, shTXN1 and shTXN3 were very effective in downregulating thioredoxin expression in both MM.1R-BTZ and RPMI8226/Dox- BTZ cells and were chosen for subsequent experiments.
MM.1R-BTZ cells transduced with shTXN1 or shTXN3 showed increased sensitivity to bortezomib after 48 and 72 h of incubation in comparison to MM.1R-BTZ cells transduced with mock control shNT (Fig. 2b, c). A similar phenotype was observed in RPMI8226/Dox-BTZ cells transduced with shTXN1 or shTXN3 (Fig. 2e, f). We further explored the effects of thioredoxin gene knockdown on overcoming bortezomib resistance in vivo. MM.1R-BTZ cells transduced with shNT and MM.1R-BTZ cells transduced with shTXN1 were subcutaneously injected into 1.5 Gy total body irradiated NOD/SCID mice. Ten days later when the tumor was established, the mice were then treated with or without bortezomib (0.5 mg/kg, i.p., twice weekly). Tumor growth was measured. MM.1R-BTZ cells transduced with shNT formed tumor mass at day 10 post inoculation and progressed rapidly up to 1100 mm3 at day 25 after xenograft into the immunodeficient mice. Treatment with bortezomib did not affect tumor growth in MM.1R-BTZ cells transduced with shNT. On the other hand, MM.1R-BTZ cells transduced with shTXN responded to bortezomib treatment and the tumor sizes of MM.1R-BTZ-shTXN1 cells were significantly smaller in size (451.75 ± 38.82 mm3 at day 25) when compared with MM.1R-BTZ-shNT tumors (Fig. 2g, p = 0.0003). These data suggested that genetic deletion of thioredoxin could abolish the cell growth advantage in BTZ-resistant multiple myeloma cells and re-sensitize these cells to bortezomib treatment.
Treatment with thioredoxin inhibitor (PX12) inhibits the growth of bortezomib-resistant myeloma cells and re-sensitizes bortezomib-resistant myeloma cells to bortezomib
We tested the effects of thioredoxin inhibitor on the survival and drug resistance in bortezomib-resistant myeloma cells in vivo and in vitro. PX12, a commercially available thioredoxin-specific inhibitor, binds to the Cys73 residue of thioredoxin, causing it to become irreversibly thioalkylated and biologically inactive . We found that PX12 inhibited the growth of MM.1R-BTZ and RPMI8226/Dox-BTZ cells in a dose-dependent manner (Fig. 3a, c). When PX12 was combined with bortezomib, it could re-sensitize bortezomib-resistant myeloma cells to bortezomib (Fig. 3b, d). We next explored the pharmacological therapeutic potential of PX12 for the treatment of BTZ-resistant multiple myeloma in vivo. For this purpose, equal volumes of MM.1R-BTZ tumor cells were xenografted into the NOD/SCID mice. Two weeks later when the tumor was formed, mice were treated with control buffer, BTZ alone (0.5 mg/kg, twice per week injection, i.p.), PX12 (12 mg/kg, twice per week, i.p.), or the combination for additional 4 weeks. As shown in Fig. 3e–g, the combination of PX12 with BTZ treatment had marked effects against MM.1R-BTZ tumor growth in NOD/SCID mice. Taken together, our findings demonstrate the potential of targeting thioredoxin to re-sensitize BTZ-resistant multiple myeloma to bortezomib treatment.
Inhibition of thioredoxin induces mitophagy in bortezomib-resistant myeloma cells
Mitophagy, the autophagy of mitochondria, is a highly specific quality control process and is currently under extensive investigation for its role in carcinogenesis, cancer progression, and resistance to chemotherapy . Thioredoxin is an antioxidant and is important in maintaining redox homeostasis. We reasoned that inhibition of thioredoxin would affect mitochondrial function and lead to the activation of mitophagy. We further hypothesized that PX12 inhibited the growth of BTZ-resistant myeloma cells and re-sensitized BTZ-resistant myeloma cells to bortezomib through the activation of mitophagy. To test this hypothesis, we measured several mitophagy activation markers, including the conversion from LC3BI to LC3BII, the expression of PETN-induced putative kinase 1 (PINK1), and mitochondria membrane potential (Δψm by JC-1) as well as the fusion of mitochondria and autophagosome under transmission electron microscopy.
Microtubule-associated protein 1A/1B-light chain 3 (LC3) is a soluble protein that is distributed ubiquitously in cells. During autophagy, the cytosolic form of LC3 (LC3-I) is conjugated to phosphatidylethanolamine to form LC3-phosphatidylethanolamine conjugate (LC3-II). Increase in LC3-II and/or reduction in LC3-I has become a reliable marker for autophagy and autophagy-related process . MM.1R-BTZ cells and RPMI8226/Dox-BTZ cells were treated with bortezomib, PX12, or the combination of bortezomib and PX12, and LC3B-I and LC3B-II were measured by western blot analysis. As shown in Fig. 4a, b, PX12 treatment significantly increased LC3B-II expression compared to that in control group. The combination of PX12 and bortezomib had further increased LC3B-II expression in MM.1R-BTZ cells and RPMI8226/Dox-BTZ cells. PINK1 is a mitochondrial serine threonine kinase and serves as a mitophagy-specific marker . Similar to LC3B-II, PX12 treatment increased PINK1 expression and the combination of PX12, and bortezomib had augmented effect on PINK1 expression in both MM.1R-BTZ- and RPMI8226/Dox-BTZ-resistant cells (Fig. 4a, b), indicating the activation of mitophagy.
To further elucidate the role of PX12 toward autophagic response in BTZ-resistant multiple myeloma, we applied transmission electron microscopy to visualize the mitophagy process and quantitated autolysosomes per 4 μm2 area. As shown in Fig. 4c, there were fewer intracellular autophagic vacuoles in MM.1R-BTZ cells treated with bortezomib. The mitochondrial was intact and separated from autophagosome. On the other hand, there was a significant increase (by 2.3 fold) in the autophagic response in MM.1R-BTZ cells treated with the combination of PX12 and BTZ, reflected by an increase in the amount of autophagy and lysosomes and the engulfed mitochondria by autophagosome (Fig. 4d). Additionally, we measured mitochondrial membrane potential (Δψm) using MitoProbe staining in MM.1R-BTZ cells treated with BTZ alone, PX12, or the combination. When mitochondrial membrane depolarizes (i.e., the loss of Δψm), JC-1 aggregates accumulate as demonstrated by the increase in the JC-1 fluorescence intensity in Q4 area (Fig. 4e–h). PX12 treatment as well as the combination of PX12 and BTZ led to accelerated Δψm loss. The percentage of MM.1R-BTZ cells undergoing mitophagy (the cells in Q4) drastically increase from 12.6% in bortezomib alone treatment group to 30.1% in bortezomib+PX12 combination group (Fig. 4f, h).
To determine if the combination of PX12 and bortezomib induced mitophay in vivo in tumor xenograft model, the tumors as described in Fig. 3 were harvested at the end of experiments and LC3B-I and PINK1 expression were measured by western blot analysis. As shown in Fig. 4i, LC3B-I was significantly reduced in the PX12+BTZ treatment group. Furthermore, PINK1 expression was significantly increased in PX12+BTZ treatment group. These data demonstrated the activation of mitophagy in vivo with PX12+BTZ treatment, which was consistent with our in vitro findings.
To determine if mitophagy activation mediates the effects of PX12+BTZ combination in BTZ-resistant multiple myeloma cells, we used bafilomycin, a mitophagy inhibitor, to block mitophagy. MM.1R-BTZ myeloma cells and RPMI8226/Dox-BTZ myeloma cells were treated with PX12, BTZ, or PX12+BTZ with or without 5 μM of bafilomycin. As shown in Fig. 4j, k, bafilomycin reversed the inhibitory effect mediated by PX12+BTZ treatment. Taken together, our results supported our hypothesis that thioredoxin inhibition enhances bortezomib sensitivity in BTZ-resistant myeloma cells via mitophagy activation.
Similarly, genetic deletion of TXN in MM.1R- and RPMI8226/Dox-BTZ myeloma cells resulted in increased expression of PINK1 and LC3B-II, and the mitochondrial membrane depolarization as measured by JC-1 aggregates (Fig. 5), consistent with the induction of mitophagy.
Inhibition of thioredoxin decreases mTOR and ERK1/2 phosphorylation
To explore the molecular mechanisms by which thioredoxin inhibition leads to the activation of mitophagy and cytotoxicities in BTZ-resistant multiple myeloma cells, we established an in silico genetic interaction network-predicted method using String (https://string-db.org/). This network analysis predicts several key hubs associated with mitophagy and thioredoxin according to the literatures search. Figure 6a, b revealed the predicted key proteins including mTOR, ERK1/2, and PINK1 and their potential relationships with mitophagy. We subsequently performed western blot analysis to validate these interactions. Our results indicated that although the total protein levels were not affected, the phosphorylation level of mTOR (Fig. 6c, d) and ERK1/2 (Fig. 6e, f) decreased dramatically in MM.1R-BTZ- and RPMI8226/Dox-BTZ-resistant cells treated with BTZ+PX12, when compared to those in the BTZ group.
To better understand whether the reduced levels of mTOR and ERK1/2 phosphorylation are directly involved in the anti-myeloma effects of thioredoxin inhibition in bortezomib-resistant myeloma cells, we treated MM.1R-BTZ and RPMI8226/Dox-BTZ-resistant cells with the BTZ+PX12 combination in the absence and presence of various concentration of mTOR activator (MHY1485) or ERK activator (tBHQ). As shown in Fig. 6g–j, the inhibitory effects of BTZ and PX12 co-treatment were reversed by MHY1485 or tBHQ in a dose-dependent manner. Thus, our data demonstrate that mTOR and ERK1/2 phosphorylation could play a role in the effects mediated by BTZ+PX12 treatment in bortezomib-resistant myeloma cells.
High expression of TXN is associated with poor prognosis in multiple myeloma patients
To investigate the clinical significance of TNX in multiple myeloma, we interrogated publicly available microarray datasets and checked TNX mRNA level expression between normal plasma cells and multiple myeloma using the Oncomine database. As shown in Fig. 7, TNX expression was significantly increased in the CD138+ plasma cells of patients with myeloma compared to normal plasma cells and tonsillar lymphoid tissues in two different datasets [24, 25] (Fig. 7a, b).
To determine the expression of thioredoxin in bortezomib relapsed/refractory myeloma patients, we isolated CD138+ cells and CD138− cells from untreated newly diagnosed myeloma patients (n = 3), myeloma patients who have just completed bortezomib containing regimen induction therapy with residual myeloma cells (n = 5), and relapsed/refractory myeloma patients (n = 5) who had previous exposure and resistant to bortezomib. When compared to the CD138− cells, TNX mRNA was significantly increased in the CD138+ plasma cells, suggesting that TXN is upregulated in myeloma cells in comparison to non-myelomatous cells of the same patient (Fig. 7c). Furthermore, consistent with the findings in our adaptive BTZ-resistant myeloma cell lines, TXN expression demonstrated a trend for elevation in CD138+ myeloma cells after treatment with bortezomib or being refractory to bortezomib (Fig. 7d) TNX expression has no notable difference in CD138− non-myelomatous cells among new diagnosed, treated, and refractory myeloma patients (Fig. 7e).
In addition, we analyzed the correlation between thioredoxin expression and overall survival using the Mulligan et al.  and Arkansas  microarray datasets. Compared to myeloma patients with low TXN expression, patients with higher TXN expression was associated with a shorter overall survival (Fig. 7f, g). Altogether, these findings demonstrated the clinical significance for TXN in myeloma, especially in bortezomib-resistant multiple myeloma.
Multiple myeloma remains an incurable disease and nearly all myeloma patients will eventually develop resistance to currently available therapeutic agents including proteasome inhibitors. Bortezomib is the first proteasome inhibitor used in the treatment of myeloma, and it has dramatically changed the landscape of the care of patients with multiple myeloma. Bortezomib can induce deeper response, lead to higher response rate, and extend the survival of patients with newly diagnosed multiple myeloma and of patients with relapsed myeloma. Unfortunately, development of drug resistance to bortezomib is inevitable and poses a major challenge in our continuous improvement of clinical outcomes for patients with multiple myeloma. Therefore, it is critical and imperative to elucidate the molecular mechanisms associated with or driving the development of bortezomib resistance. This information is essential in our efforts to overcome or re-sensitize bortezomib resistance.
In the current study, we generated over a period of 1.5 years several adaptive bortezomib-resistant myeloma cell lines (Fig. 1). We have also analyzed and compared the expression of thioredoxin in CD138+ myeloma cells from newly diagnosed myeloma patients, myeloma patients that were treated with bortezomib containing regimen, and relapsed/refractory myeloma patients who had prior exposure to bortezomib (Fig. 7). We have demonstrated that over-expression of thioredoxin was associated with the development of drug resistance to bortezomib. Thioredoxin was previously found to promote tumor growth through inhibition of apoptosis, reduce sensitivity of the tumor to drugs [28,29,30], and be associated with poor prognosis . Raninga et al. has recently shown that compared to normal plasma cells, multiple myeloma cells had higher intrinsic oxidative stress and higher levels of thioredoxin and thioredoxin reductase expression and that thioredoxin over-expression was associated with resistance to NF-κβ inhibitors . Using comparative proteomic profiling, Dytfeld et al. showed upregulation of thioredoxin expression in bortezomib-resistant myeloma cells . Our current studies were consistent with and validated these observations. Importantly, our current studies advance significantly from those studies. We indeed showed that inhibition of thioredoxin by shRNA knockdown or pharmacological approach with PX12 could re-sensitize the bortezomib-resistant myeloma cells to bortezomib in vitro and in vivo in myeloma xenograft mouse models (Figs. 2 and 3). Mechanistically, we have demonstrated that inhibition of thioredoxin re-sensitize bortezomib-resistant myeloma cells through the activation of mitophagy. Our studies provide molecular rationale and justification for targeting thioredoxin in the treatment of bortezomib relapsed/refractory myeloma.
We found that PX12-sensitized BTZ-resistant myeloma cells to bortezomib through the activation of mitophagy. We confirmed the activation of mitophagy using several methods including LC3 western blotting, PINK1 expression analysis, measurement of mitochondrial membrane potential, and the visualization of mitochondrial-autophagosome fusion with transmission electron microscopy (Fig. 4). Using bafilomycin to inhibit mitophagy, we have demonstrated the important role of mitophagy in PX12-mediated effects. Thioredoxin-specific shRNA knockdown revealed similar findings on the induction of mitophagy (Fig. 5). Thioredoxin and thioredoxin system maintain the intracellular redox homeostasis by scavenging ROS and regulating other redox proteins [31, 32]. Inhibition of thioredoxin increases intracellular oxidative stress [33, 34]. Increased production of ROS stimulates the initiation of mitophagy . Mitophagy is a cellular self-cannibalization process that captures and digests mitochondrial in lysosomes . The role of mitophagy in tumorigenesis remains a topic of debate. On one hand, mitophagy is activated in transformed cells and is beneficial for tumor maintenance and progression. On the other hand, excessive autophagy can act as a tumor-suppressive mechanism possibly through initiation of cell death or senescence [36,37,38,39,40]. The outcome of mitophagy activation in cancer depends on the stage of the disease, cell types, oncogenic drivers, and the intensity of the activation signal [35, 41, 42]. Mitophagy is also found to be involved in chemoresistance. Our current study suggested that reduced level of mitophagy is associated with the development of drug resistance in multiple myeloma cells and that activation of mitophagy could re-sensitize myeloma cells to bortezomib killing. Our study provides evidence for future exploring mitophagy pathway for the treatment of relapsed and refractory multiple myeloma.
In the current study, we delineated the molecular events involved in the PX12 + BTZ- induced mitophagy in MM cells. We described 2 signal proteins (p-mTOR and p-ERK1/2) that may play a role in PX12 + BTZ -induced mitophagy and anti-myeloma effects in MM cells. Using chemical activator of mTOR and ERK, we demonstrated an important role of mTOR and ERK in the anti-myeloma effects induced by thioredoxin inhibition. Previous studies have demonstrated that various signal pathways are involved in autophagy, including PI3K/Akt/mTOR, ERK1/2 and NF-κB. mTOR (mechanistic target of rapamycin kinase) has been shown to control multiple cellular functions such as gene transcription, protein formation, cell proliferation and senescence and cell metabolism . mTOR pathway is the master regulator of cellular metabolism and has been shown to be the major pathway regulating mitophagy. AKT was previously reported to control mitochondrial biogenesis and autophagy. In our system, we did not observe changes in AKT expression or phosphorylation after PX12 treatment. In addition, we found that treatment with PX12 resulted in desphosphorylation of ERK1/2. Our study suggested that inhibition of thioredoxin affects several molecular pathways that are important in the regulation of mitophagy.
ERK1 (p44) and ERK2 (p42) are two isoforms of ERK that are activated downstream of Ras in response to extracellular cues. ERK has been shown to induce autophagy in response to a number of anti-tumor/cytotoxic agents [44, 45]. Inhibition of ERK was associated with a decrease in autophagy and increased cellular sensitivity to tumor necrosis factor-α (TNF) in breast cancer MCF-7 cells . Interestingly, our data showed that dephosphorylation of ERK1/2 appears to be associated with the induction of mitophagy in BTZ-resistant MM cells. Additional studies are needed to further understand the ERK1/2 pathway in BTZ resistance in multiple myeloma. The discrepancy between reduced activity of ERK1/2 in our current study and the findings in breast cancer could be due to different cancer type or to different drug sensitivity.
In summary, we showed that upregulation of thioredoxin expression correlated with the development of bortezomib drug resistance in multiple myeloma. Inhibition of thioredoxin by shRNA or thioredoxin inhibitor re-sensitized bortezomib-resistant myeloma cells to bortezomib treatment. We further showed that the combination of bortezomib and thioredoxin inhibitor suppressed the growth of bortezomib-resistant myeloma cells through mitophagy activation. Our study provides new insights in the molecular mechanisms associated with bortezomib drug resistance. Our work was the first to demonstrate the novel role of mitophagy in multiple myeloma treatment.
Extracellular signal-regulated protein kinases 1 and 2
Mechanistic target of rapamycin kinase
PTEN-induced putative kinase 1
Tumor necrosis factor-α
Katzel JA, Hari P, Vesole DH. Multiple myeloma: charging toward a bright future. CA Cancer J Clin. 2007;57(5):301–18.
Wei G, Wang L, Yang H, Han X, Zheng G, Zheng W, Sun J, Shi J, Wu W, Zhao Y, et al. Clinical implications of c-maf expression in plasma cells from patients with multiple myeloma. Exper hematol oncol. 2017;6:16.
Delgado-Calle J, Bellido T, Roodman GD. Role of osteocytes in multiple myeloma bone disease. Curr Opin Support Palliat Care. 2014;8(4):407–13.
Terpos E, Berenson J, Raje N, Roodman GD. Management of bone disease in multiple myeloma. Expert Rev Hematol. 2014;7(1):113–25.
Moreau P, Richardson PG, Cavo M, Orlowski RZ, San Miguel JF, Palumbo A, Harousseau JL. Proteasome inhibitors in multiple myeloma: 10 years later. Blood. 2012;120(5):947–59.
Richardson PG, Sonneveld P, Schuster MW, Irwin D, Stadtmauer EA, Facon T, Harousseau JL, Ben-Yehuda D, Lonial S, Goldschmidt H, et al. Bortezomib or high-dose dexamethasone for relapsed multiple myeloma. N Engl J Med. 2005;352(24):2487–98.
Rastgoo N, Abdi J, Hou J, Chang H. Role of epigenetics-microRNA axis in drug resistance of multiple myeloma. J Hematol Oncol. 2017;10(1):121.
Lu J, Holmgren A. The thioredoxin antioxidant system. Free Radic Biol Med. 2014;66:75–87.
Karlenius TC, Tonissen KF. Thioredoxin and cancer: a role for thioredoxin in all states of tumor oxygenation. Cancers (Basel). 2010;2(2):209–32.
Wangpaichitr M, Sullivan EJ, Theodoropoulos G, Wu C, You M, Feun LG, Lampidis TJ, Kuo MT, Savaraj N. The relationship of thioredoxin-1 and cisplatin resistance: its impact on ROS and oxidative metabolism in lung cancer cells. Mol Cancer Ther. 2012;11(3):604–15.
Kim SJ, Miyoshi Y, Taguchi T, Tamaki Y, Nakamura H, Yodoi J, Kato K, Noguchi S. High thioredoxin expression is associated with resistance to docetaxel in primary breast cancer. Clin cancer res off j Am Assoc Cancer Res. 2005;11(23):8425–30.
Jiang P, Mizushima N. Autophagy and human diseases. Cell Res. 2014;24(1):69–79.
Yoshida GJ. Therapeutic strategies of drug repositioning targeting autophagy to induce cancer cell death: from pathophysiology to treatment. J Hematol Oncol. 2017;10(1):67.
Grek CL, Tew KD. Redox metabolism and malignancy. Curr Opin Pharmacol. 2010;10(4):362–8.
Kryston TB, Georgiev AB, Pissis P, Georgakilas AG. Role of oxidative stress and DNA damage in human carcinogenesis. Mutat Res. 2011;711(1–2):193–201.
Wu W, Tian W, Hu Z, Chen G, Huang L, Li W, Zhang X, Xue P, Zhou C, Liu L, et al. ULK1 translocates to mitochondria and phosphorylates FUNDC1 to regulate mitophagy. EMBO Rep. 2014;15(5):566–75.
Maes H, Rubio N, Garg AD, Agostinis P. Autophagy: shaping the tumor microenvironment and therapeutic response. Trends Mol Med. 2013;19(7):428–46.
Bordi M, Nazio F, Campello S. The close interconnection between mitochondrial dynamics and mitophagy in cancer. Front Oncol. 2017;7:81.
Gu Z, Xia J, Xu H, Frech I, Tricot G, Zhan F. NEK2 promotes aerobic glycolysis in multiple myeloma through regulating splicing of pyruvate kinase. J Hematol Oncol. 2017;10(1):17.
Klionsky DJ, Abdelmohsen K, Abe A, Abedin MJ, Abeliovich H, Acevedo Arozena A, Adachi H, Adams CM, Adams PD, Adeli K, et al. Guidelines for the use and interpretation of assays for monitoring autophagy (3rd edition). Autophagy. 2016;12(1):1–222.
Kirkpatrick DL, Kuperus M, Dowdeswell M, Potier N, Donald LJ, Kunkel M, Berggren M, Angulo M, Powis G. Mechanisms of inhibition of the thioredoxin growth factor system by antitumor 2-imidazolyl disulfides. Biochem Pharmacol. 1998;55(7):987–94.
Tanida I, Ueno T, Kominami E. LC3 conjugation system in mammalian autophagy. Int J Biochem Cell Biol. 2004;36(12):2503–18.
Lazarou M, Sliter DA, Kane LA, Sarraf SA, Wang C, Burman JL, Sideris DP, Fogel AI, Youle RJ. The ubiquitin kinase PINK1 recruits autophagy receptors to induce mitophagy. Nature. 2015;524(7565):309–14.
Dickens NJ, Walker BA, Leone PE, Johnson DC, Brito JL, Zeisig A, Jenner MW, Boyd KD, Gonzalez D, Gregory WM, et al. Homozygous deletion mapping in myeloma samples identifies genes and an expression signature relevant to pathogenesis and outcome. Clin cancer res off j Am Assoc Cancer Res. 2010;16(6):1856–64.
Zhan F, Barlogie B, Arzoumanian V, Huang Y, Williams DR, Hollmig K, Pineda-Roman M, Tricot G, van Rhee F, Zangari M, et al. Gene-expression signature of benign monoclonal gammopathy evident in multiple myeloma is linked to good prognosis. Blood. 2007;109(4):1692–700.
Mulligan G, Mitsiades C, Bryant B, Zhan F, Chng WJ, Roels S, Koenig E, Fergus A, Huang Y, Richardson P, et al. Gene expression profiling and correlation with outcome in clinical trials of the proteasome inhibitor bortezomib. Blood. 2007;109(8):3177–88.
Hanamura I, Huang Y, Zhan F, Barlogie B, Shaughnessy J. Prognostic value of cyclin D2 mRNA expression in newly diagnosed multiple myeloma treated with high-dose chemotherapy and tandem autologous stem cell transplantations. Leukemia. 2006;20(7):1288–90.
Raninga PV, Di Trapani G, Vuckovic S, Bhatia M, Tonissen KF. Inhibition of thioredoxin 1 leads to apoptosis in drug-resistant multiple myeloma. Oncotarget. 2015;6(17):15410–24.
Dytfeld D, Luczak M, Wrobel T, Usnarska-Zubkiewicz L, Brzezniakiewicz K, Jamroziak K, Giannopoulos K, Przybylowicz-Chalecka A, Ratajczak B, Czerwinska-Rybak J, et al. Comparative proteomic profiling of refractory/relapsed multiple myeloma reveals biomarkers involved in resistance to bortezomib-based therapy. Oncotarget. 2016;7(35):56726–36.
Baker A, Payne CM, Briehl MM, Powis G. Thioredoxin, a gene found overexpressed in human cancer, inhibits apoptosis in vitro and in vivo. Cancer Res. 1997;57(22):5162–7.
An N, Kang Y. Thioredoxin and hematologic malignancies. Adv Cancer Res. 2014;122:245–79.
Das KC, Das CK. Thioredoxin, a singlet oxygen quencher and hydroxyl radical scavenger: redox independent functions. Biochem Biophys Res Commun. 2000;277(2):443–7.
Yamamoto M, Yang G, Hong C, Liu J, Holle E, Yu X, Wagner T, Vatner SF, Sadoshima J. Inhibition of endogenous thioredoxin in the heart increases oxidative stress and cardiac hypertrophy. J Clin Invest. 2003;112(9):1395–406.
You BR, Kim SH, Park WH. Reactive oxygen species, glutathione, and thioredoxin influence suberoyl bishydroxamic acid-induced apoptosis in A549 lung cancer cells. Tumour Biol. 2015;36(5):3429–39.
White E. Deconvoluting the context-dependent role for autophagy in cancer. Nat Rev Cancer. 2012;12(6):401–10.
Elgendy M, Sheridan C, Brumatti G, Martin SJ. Oncogenic Ras-induced expression of Noxa and Beclin-1 promotes autophagic cell death and limits clonogenic survival. Mol Cell. 2011;42(1):23–35.
Choi AM, Ryter SW, Levine B. Autophagy in human health and disease. N Engl J Med. 2013;368(7):651–62.
Zeng R, Chen Y, Zhao S, Cui GH. Autophagy counteracts apoptosis in human multiple myeloma cells exposed to oridonin in vitro via regulating intracellular ROS and SIRT1. Acta Pharmacol Sin. 2012;33(1):91–100.
Dykstra KM, Allen C, Born EJ, Tong H, Holstein SA. Mechanisms for autophagy modulation by isoprenoid biosynthetic pathway inhibitors in multiple myeloma cells. Oncotarget. 2015;6(39):41535–49.
Chen S, Zhou L, Zhang Y, Leng Y, Pei XY, Lin H, Jones R, Orlowski RZ, Dai Y, Grant S. Targeting SQSTM1/p62 induces cargo loading failure and converts autophagy to apoptosis via NBK/Bik. Mol Cell Biol. 2014;34(18):3435–49.
Song IS, Kim HK, Lee SR, Jeong SH, Kim N, Ko KS, Rhee BD, Han J. Mitochondrial modulation decreases the bortezomib-resistance in multiple myeloma cells. Int J Cancer. 2013;133(6):1357–67.
Jia L, Gopinathan G, Sukumar JT, Gribben JG. Blocking autophagy prevents bortezomib-induced NF-kappaB activation by reducing I-kappaBalpha degradation in lymphoma cells. PLoS One. 2012;7(2):e32584.
Laplante M, Sabatini DM. mTOR signaling at a glance. J Cell Sci. 2009;122(Pt 20):3589–94.
Choi CH, Jung YK. Oh SH: autophagy induction by capsaicin in malignant human breast cells is modulated by p38 and extracellular signal-regulated mitogen-activated protein kinases and retards cell death by suppressing endoplasmic reticulum stress-mediated apoptosis. Mol Pharmacol. 2010;78(1):114–25.
Ellington AA, Berhow MA, Singletary KW. Inhibition of Akt signaling and enhanced ERK1/2 activity are involved in induction of macroautophagy by triterpenoid B-group soyasaponins in colon cancer cells. Carcinogenesis. 2006;27(2):298–306.
Sivaprasad U, Basu A. Inhibition of ERK attenuates autophagy and potentiates tumour necrosis factor-alpha-induced cell death in MCF-7 cells. J Cell Mol Med. 2008;12(4):1265–71.
The authors thank Dr. Benny Chen and his lab members for the technical support. Besides, we also want to thank Gene Expression Omnibus (GEO) and Oncomine databases for making their data readily available to the scientific community.
This work is supported by Duke Cancer Institute Fund (to NC and YK), NIH 1K08HL 103780-01A1 (to YK), NIH R44CA199767 (to YK), and NIH R01CA197792 (to YK). The project is supported in part by the National and Fujian Provincial Key Clinical Specialty Discipline Construction Program, China.
Availability of data and materials
The dataset(s) supporting the findings of this study are included within the article.
Ethics approval and consent to participate
The study was approved by the Duke University Institutional Review Board Committees and written informed consents were obtained before any operation to patients.
Consent for publication
The bone marrow samples were obtained under de-identifiers. No patient personal information was involved.
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Effect of BTZ on cell survival and thioredoxin expression on the parental and resistant multiple myeloma cell lines following 48 h of BTZ treatment in vitro. (A–C) Generation of adaptive bortezomib-resistant myeloma cell lines: MM.1S cells (A), OPM1 cells (B), and NCIH929 (C) were cultured with serially increased concentrations of bortezomib over a period of 1.5 years. The parental and bortezomib-resistant myeloma cells were treated with various concentration of bortezomib for 48 h and the IC50 for parental and bortezomib-resistant myeloma cells were calculated and shown. (D–F) Thioredoxin mRNA expression is upregulated in bortezomib (BTZ)-resistant MM.1S (D), OPM1 (E), and NCIH929 (F) cells. Parental cells and BTZ-resistant MM.1S, OPM1, and NCIH929-resistant RPMI8226/Dox cells at the highest bortezomib concentration were harvested and mRNA expression measured by RT-PCR. Error bars, standard error of the mean (SEM). (G–I) Thioredoxin protein expression increased in parallel to increased bortezomib drug resistance. Serial MM.1S (G), OPM1 (H) and NCIH929 (I) cells that were cultured with increased concentrations of bortezomib were harvested and protein expression measured by western blot analysis. The intensity of expression was semi-quantitated using Image-Pro Plus 6.0 software and adjusted to β-actin. Error bars, standard error of the mean (SEM); *p < 0.05, **p < 0.01, ***p < 0.005, ****p < 0.001. (TIFF 239 kb)
About this article
Cite this article
Zheng, Z., Fan, S., Zheng, J. et al. Inhibition of thioredoxin activates mitophagy and overcomes adaptive bortezomib resistance in multiple myeloma. J Hematol Oncol 11, 29 (2018). https://doi.org/10.1186/s13045-018-0575-7