- Open Access
Artesunate shows potent anti-tumor activity in B-cell lymphoma
Journal of Hematology & Oncologyvolume 11, Article number: 23 (2018)
Although chemo-immunotherapy has led to an improved overall survival for most B-cell lymphoma types, relapsed and refractory disease remains a challenge. The malaria drug artesunate has previously been identified as a growth suppressor in some cancer types and was tested as a new treatment option in B-cell lymphoma.
We included artesunate in a cancer sensitivity drug screen in B lymphoma cell lines. The preclinical properties of artesunate was tested as single agent in vitro in 18 B-cell lymphoma cell lines representing different histologies and in vivo in an aggressive B-cell lymphoma xenograft model, using NSG mice. Artesunate-treated B lymphoma cell lines were analyzed by functional assays, gene expression profiling, and protein expression to identify the mechanism of action.
Drug screening identified artesunate as a highly potent anti-lymphoma drug. Artesunate induced potent growth suppression in most B lymphoma cells with an IC50 comparable to concentrations measured in serum from artesunate-treated malaria patients, while leaving normal B-cells unaffected. Artesunate markedly inhibited highly aggressive tumor growth in a xenograft model. Gene expression analysis identified endoplasmic reticulum (ER) stress and the unfolded protein response as the most affected pathways and artesunate-induced expression of the ER stress markers ATF-4 and DDIT3 was specifically upregulated in malignant B-cells, but not in normal B-cells. In addition, artesunate significantly suppressed the overall cell metabolism, affecting both respiration and glycolysis.
Artesunate demonstrated potent apoptosis-inducing effects across a broad range of B-cell lymphoma cell lines in vitro, and a prominent anti-lymphoma activity in vivo, suggesting it to be a relevant drug for treatment of B-cell lymphoma.
More than 60 different B-cell lymphoma types have been described, with different clinical courses and outcomes . Diffuse large B-cell lymphoma (DLBCL) and follicular lymphoma (FL) are the most common forms of aggressive non-Hodgkin lymphoma (NHL) and indolent NHL, respectively. Current standard of care often involves chemotherapy regimens combined with monoclonal antibodies targeting the B-cell surface molecule CD20 (rituximab) and has led to an improved overall survival for most lymphoma types . A number of small molecule inhibitors are currently under clinical testing, where the majority represents inhibitors of the B-cell receptor signaling (BCR) pathway and downstream signaling molecules .
Artesunate is a semisynthetic analogue of artemisinin, extracted from the sweet wormwood plant, Artemisia annua, and is the first-line treatment of severe malaria . Artemisinins have demonstrated cytotoxic effects in vitro against a range of cell lines and also to be effective against drug-resistant cancer cell lines [5, 6], although the underpinning mechanisms remain unclear [7, 8]. Few studies have been reported regarding clinical use of artesunate in cancer, but artesunate has demonstrated anti-tumor effects in a small randomized clinical trial in colon cancer and to transiently decrease tumor size and prostate-specific antigen levels in a patient with advanced prostate cancer [9, 10]. The supply of artesunate has been limited, but semisynthetic production since 2014 has increased the availability  and enables additional clinical use of artesunate beyond malaria treatment.
B lymphoma cell lines are representative models with > 80% match with the corresponding primary cancer cells [12, 13] (O. Kallioniemi, pers. com.). Here, we demonstrate a potent induction of apoptosis by artesunate in a broad range of cell lines, representing the most common types of B-cell lymphoma. Artesunate also showed potent anti-tumor efficacy in vivo in a xenograft model, providing a rationale for clinical testing as B-cell lymphoma therapy.
Small molecule inhibitors are as follows: everolimus, dasatinib, duvelisib, alisertib, ibrutinib, idelalisib, sorafenib, metformin, and entospletinib (Selleckchem, Houston, TX; further information in Additional file 1: Table S1). All chemicals were from Sigma-Aldrich (St. Louis, MD) unless otherwise noted. Primary antibodies are as follows: rabbit anti-ATF-4 (#11815), -ATF-6 (#65880) and –β-tubulin (#2128) and mouse-anti-DDIT3/CHOP (#2895) (Cell Signaling Technology, Danvers, MA), rabbit anti-GAPDH (#100118; GeneTex, Irvine, CA), mouse anti-HSP70 (#648001; BioLegend, San Diego, CA), goat anti-Lamin-B (#SC-6217; Santa Cruz Biotechnology, Dallas, TX). Secondary antibodies are as follows: HRP-conjugated goat anti-rabbit (#111-035), goat anti-mouse (#115-035), and donkey anti-goat (#705-035) IgG antibodies (Jackson Immunoresearch, West Grove, PA).
Cell lines, human samples, and culture conditions
Cell lines representing the most common B-cell lymphoma subtypes are the following: Burkitt lymphoma: BL-41, Raji, Ramos, Rec-1 (Leibniz-Institut-Deutche Sammlung von Mikroorganismen und Zellkulturen (DSMZ)) and Namalwa (gift from J. Delabie); Germinal Center B-like (GCB) DLBCL: ULA, SU-DHL-6, Oci-Ly-18 (DSMZ), Oci-Ly-2, SU-DHL-4 (gift from L. Staudt); activated B-cell like (ABC) DLBCL: U-2932 (DSMZ); DLBCL-double hit: Will-2 (DSMZ); mantle cell lymphoma: Mino, JeKo-1, Granta-519 (DSMZ); FL: Ros-50 (DSMZ), K-422 (gift from J. Delabie); immunoblastic B-cell lymphoma: DoHH-2 (DSMZ). K-422, Oci-Ly-2, Namalwa and Ros-50 were authenticated by STR-profile in 2011, BL-41, SU-DHL-4, SU-DHL-6, Ramos, Raji, Granta, WILL-2, Rec-1. Mino and DoHH-2 were DNA fingerprinted (December 2017, Forensik, Eurofins Medigenomix, Ebersberg, Germany). The cell lines obtained as gifts were second generation passages. For routine use of the cell lines, second generation splits were thawed and the experiments were performed in the time frame of 1 week after thawing and to maximum 20 passages. Mycoplasma testing is always performed after 4 weeks for all cultures using PCR-based methods (Venor GeM Mycoplasma Detection Kit), and always prior to injection of cells into mice. Cell lines were cultured in RPMI-1640 supplemented with 10% fetal calf serum (FCS) or 10% human serum (HS; TCS Biosciences, Buckingham, UK), penicillin, and streptomycin (complete media) and maintained at 37 °C in 5% CO2. B-cells were isolated from buffy coats from peripheral blood from anonymous, healthy donors at the Blood Bank (Oslo University Hospital, Norway), after informed consent and approval from regional authorities (REK S-03280). B-cells were isolated using anti-CD19 Dynabeads (Thermo Fisher Scientific, Waltham, MA), followed by the corresponding DETACHaBEADs . B-cells were maintained in complete media with rhCD40 ligand (CD40L) (0.25 μg/ml; Enzo Life Sciences, Farmingdale, NY) preincubated for 30 min with an enhancer for ligands (1 μg/ml; Enzo Life Sciences) and recombinant human interleukin (IL)21 (50 ng/ml; Thermo Fisher Scientific).
Luciferase construct and BL-41-luc cell line
The gene encoding for luciferase-TropC-GFP  was subcloned into pENTR vector and further recombined into pMP71 retroviral vector . Retroviral particles were prepared, and BL-41 cells were transduced by double spinoculation . Positive cells were sorted for GFP.
Viability and apoptosis assays
Cells were grown in 96-well plates (10,000 cells/well) for 72 h with or without artesunate (0.125–8 μM) and small molecular inhibitors (Additional file 1: Table S1). The CellTiterGlo Assay (Promega, Madison, WI) was used for measuring relative cell growth (ATP levels), using Modulus microplate reader (Turner Biosystems, Promega). Measurements were related to control and reported as relative luminescence units (RLU). Viability was measured by propidium iodide (PI) assay (1 μg/ml; Thermo Fisher Scientific) after treatment with artesunate (5–10 μM) for 72 h. Apoptosis was detected by active caspase-3 assay (Alexa-647-rabbit anti-active caspase-3 clone C92-605; BD Biosciences, CA) after 24 h treatment with artesunate (1 and 5 μM) with and without Z-VAD-FMK (1 μl/ml (ab65613), Abcam, Cambridge, UK). Terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) was determined by the In Situ Cell Death Detection Kit after 72 h treatment with artesunate (10 μM). For both apoptosis assays, cells were fixed in PFA for 5 min and permeabilized with ice cold methanol. All assays were analyzed using FACS Canto II or LSR II flow cytometer (BD). Flow cytometry data were analyzed using the online Cytobank flow cytometry software (www.cytobank.org) . DL-buthionine sulfoximine (BSO) was used at 50 μM for 18–24 h for reducing the cellular glutathione levels by inhibiting γ-glutamylcysteine synthetase. Glutathione levels were measured by GSH-Glo Glutathione assay (Promega).
Gene expression profiling
Total RNA was isolated using MiRNeasy (Qiagen, Hilden, Germany) after 4 and 12 h treatment with artesunate (5 μM). The microarray analyses were performed on Illumina’s HumanHT-12 v4 Expression BeadChip platform (Illumina, San Diego, CA). The differential gene expression analysis was done using the limma package  in R (version 3.3.1). Two technical replicates were run per cell line and time point, and the results were averaged together. The differentially expressed genes were selected based on a log fold change larger than the absolute value of 0.5, and an adjusted p value (FDR) of less than 0.01. The probes were collapsed according to gene symbol, using the annotation file for the Illumina’s HumanHT-12 v4 Expression BeadChip platform. When several probes mapped to the same gene, the probe with lowest log fold change was selected. The pathways and networks most enriched for the differential expressed genes were identified by Ingenuity Pathway Analysis (IPA) software (Qiagen) with default settings. Microarray data is available at NCBI’s Gene Expression Omnibus with accession number GSE94553 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE94553).
Cells were lysed and processed for SDS-PAGE . Miniprotean or Criterion TGX precast gels were used for SDS-PAGE (Bio-Rad Laboratories, Hercules, CA). SuperSignal West Pico or Dura (Thermo Fisher Scientific) or Clarity (Bio-Rad) was used for detection. Chemidoc MP (Bio-Rad) was applied for imaging, and image processing was performed in ImageLab (Bio-Rad), Adobe Photoshop, and Adobe Illustrator (Adobe Systems, San Jose, CA).
The care and handling of animals for the present study were in conformity with the Norwegian Food Safety Authority in compliance with the European Convention of the Protection of Vertebrates Used for Scientific Purposes (Project ID 7729). NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ (NSG) mice were bred in-house. Pilot experiments were performed with three mice in each group. Based on the results, n = 10 for treatment and control group was chosen. The mice (6–10 weeks old) were injected subcutaneously with 2 × 106 BL-41 cells expressing firefly luciferase (BL-41-luc). Tumor take was measured by IVIS at day 4 before mice were divided into control and treatment groups. The mice were divided according to tumor size within each cage for either treatment or control group to have non-biased, comparable groups. Artesunate was dissolved in EtOH/DMSO (1:1) to 370 mg/ml and diluted 1:10 in 5% Na3CO3 before injections. Mice were injected daily from day 4 with 200 mg/kg artesunate intraperitoneally or control (5% EtOH and 5% DMSO in Na3CO3). Treatment was given for 12 days, then 2 days off, followed by treatment every day until day 19 (17 injections total). Tumor growth was monitored by bioluminescent imaging (IVIS spectrum in vivo imaging system) at regular intervals. Mice were injected intraperitoneally with 150 mg/kg d-luciferin (Caliper Life Sciences, Hopkinton, MA) and imaged after 10 min. Inhalation anesthetic sevofluran (Baxter) supplemented with O2 and N2O was used during imaging. Caliper measurement was used to determine if the tumor size had reached maximum 2 cm in one direction or 2 cm3 that was the limit for euthanasia.
A Seahorse Extracellular Flux (XF96e) Analyzer (Agilent, Santa Clara, CA) was used to measure the oxygen consumption rate (OCR) which relates to mitochondrial respiration and the extracellular acidification rate (ECAR) reflecting the cellular glycolytic activity of live intact B-cell lymphoma cells, in real time with and without artesunate pretreatment (5 μM, 12 h). Prior to the assay, BL-41, SU-DHL-6, and Mino cells were attached onto Cell-Tak (Corning Inc., Corning, NY) coated 96-well XF-PS plates at a density of 1 × 105 cells/well in DMEM XF unbuffered assay media, supplemented with 2 mM sodium pyruvate, 10 mM glucose, 2 mM l-glutamine and adjusted to physiological pH (7.6) and incubated in the absence of CO2 for 1 h prior to the seahorse measurements (n = 21 technical replicates). Initially, the cell basal respiration was measured for all groups while after injection of oligomycin, a potent F1F0 ATPase inhibitor (1 μM), the OCR dropped by the amount required for ATP production by respiration. Subsequent addition of 1 μM carbonyl cyanide p-trifluoromethoxyphenylhydrazon (FCCP) uncoupled the mitochondrial electron transport from ATP synthesis and hence revealed the maximal respiration capacity in each case. Finally, addition of 1 μM antimycin A and rotenone completely inhibited the electron transport and hence respiration. Initially, the basal glycolytic rates of the treatment groups were measured while upon addition of oligomycin to the cells, being unable to produce ATP by respiration switched to their highest glycolytic capacity to compensate for their ATP requirements.
Statistical significance was determined by two-tailed unpaired Student’s t test and two-tailed heteroscedastic Student’s t test for Seahorse assay. One-sided t test was used for animal studies, in addition to log-rank test for the Kaplan-Meier plots. GraphPad Prism and Igor Pro were used for calculations. Differences were considered to be statistically significant if P < .05.
Artesunate has potent anti-lymphoma activity in vitro
To characterize the efficacy of artesunate, inhibitors currently in clinical testing for treatment of B-cell lymphoma were compared with artesunate in a drug sensitivity screen. Artesunate was identified as a promising candidate, due to its potent reduction of cell growth across all B lymphoma cell lines tested, which contrasted more variable efficacy of the other drugs (Fig. 1). To broaden the testing of artesunate, a dose-response experiment was performed in 18 different B-cell lymphoma cell lines, representing various histological types. This revealed that artesunate exhibited broad activity across the lymphoma cell lines and potently affected cell growth, whereas it had limited effect in normal B-cells activated with CD40L alone or in combination with IL21(Fig. 2a). Artesunate reduced cell growth with IC50 values from 0.189 to 3.72 μM, and a high sensitivity to artesunate (IC50 < 1 μM) was observed in 11 out of 18 cell lines tested (Fig. 2a, Additional file 2: Figure S1). A pronounced cell death was detected with PI staining after 72 h, with only minor effects in normal B-cells (Additional file 2: Figure S1B). Detection of DNA fragmentation by TUNEL assay demonstrated potent artesunate-induced apoptosis with percentage apoptotic cells ranging from 55 to 96 after 72 h exposure to artesunate, whereas normal B-cells were not affected (Fig. 2b). Induction of apoptosis was an early event, with prominent increase of active caspase 3+ cells after 24 h of artesunate treatment (Additional file 2: Figure S1C). The increase of active caspase 3+ cells by artesunate was counteracted by the presence of the pan caspase inhibitor Z-VAD-FMK (Fig. 2c), indicating that artesunate induces apoptosis in tumor cells via a caspase-dependent pathway. Cell cycle analysis did not reveal any overall changes induced by artesunate (Additional file 3: Figure S2). Together, these results indicated that induction of apoptosis represents the main mechanism of the anti-lymphoma activity of artesunate.
Artesunate has potent anti-lymphoma activity in a lymphoma xenograft model
To test the efficacy of artesunate as a potent anti-lymphoma drug in vivo, NSG mice were subcutaneously inoculated with 2 × 106 BL-41-luc lymphoma cells and subsequently treated with artesunate (200 mg/kg/day, n = 10) or left untreated (control, n = 10). Mice with varying tumor loads were evenly distributed among the groups (Additional file 4: Figure S3A and B), and treatment was initiated at day 4 after inoculation. Tumor growth was significantly reduced in mice treated with artesunate, and statistically significant differences were reached at day 12 and day 16 after start of treatment (P = .0045, day 12 and P < .0050, day 16; Fig. 3a). IVIS images of the mice at day 16 showed a distinct reduction in tumor growth in the artesunate-treated mice compared to the control group (Fig. 3b). In one of the artesunate-treated mice, the tumor was undetectable (Fig. 3b). At day 16, half of the control mice were euthanized due to maximum tumor size limitation (2 cm3). In comparison, the first artesunate-treated mouse was euthanized due to its tumor size 6 days later. All remaining mice were monitored after end of treatment (day 19), and a delayed tumor growth was observed in artesunate-treated mice. Survival analysis showed that the median time to reach maximum tumor size criteria was 17.5 versus 30.5 days for control and artesunate-treated mice, respectively (p < .0001; Fig. 3c). The artesunate dose (200 mg/kg/day) was well tolerated with no body weight loss during treatment (Additional file 4: Figure S3C). Taken together, artesunate showed potent anti-tumor response in an aggressive B-cell lymphoma xenograft model.
Artesunate induces unfolded protein response/endoplasmic reticulum stress
To obtain a global view of the transcriptional changes after artesunate treatment, gene expression profiling was performed on three different lymphoma cell lines BL-41, SU-DHL-4, and Oci-Ly-2 treated with artesunate or vehicle for 4 and 12 h. After 4 h, only few genes showed significantly different expression, whereas a large number of genes were affected after 12 h (521, 509, and 1029 genes had significantly altered expression levels in BL-41, SU-DHL-4, and Oci-Ly-2, respectively). Across the cell lines, 131 genes were commonly differentially expressed after 12 h and these genes had a very similar expression pattern (Fig. 4a). Ingenuity Pathway Analysis (IPA) uncovered “Unfolded protein response (UPR)” as the most significant pathway induced by artesunate in all three cell lines (P < 7.26 × 10−7). Other significant pathways identified by IPA were also related to stress response, including endoplasmic reticulum (ER) stress, tRNA charging, amino acid biosynthesis, and the protein ubiquitination pathway (Additional file 5: Table S2). For all three cell lines, DDIT3/CHOP was a central molecule in the top network (Additional file 6: Figure S4A-C). Other UPR-related genes, including activation transcription factor (ATF)-4, ATF-6, and GADD34, were also significantly upregulated, whereas a number of heat shock protein family genes were downregulated. IPA was repeated for a combined list of genes altered by artesunate in all three cell lines and still identified the UPR as the most significant pathway (Fig. 4b).
MYC was not a component of the pathways or networks affected by artesunate, when the three cell lines were analyzed separately in IPA, but was listed as an inhibited upstream regulator in Oci-Ly-2 (P = 1.5 × 10−7) (data not shown). However, when all three cell lines were analyzed together, MYC was a central molecule in the top network (Fig. 4c). MYC was also identified as number ten in the “upstream regulated network” analysis of the combined data (Additional file 7: Figure S5).
To test whether the artesunate-induced changes in gene expression also translated into altered expression of typical proteins involved in the UPR pathway, expression of ATF-4, ATF-6, DDIT3, and HSP70 were analyzed. Artesunate induced a strong upregulation of both ATF-4 and DDIT3 proteins after 18 h in BL-41 and SU-DHL-6 (Fig. 5a). ATF-6 is an ER bound transmembrane transcription factor that is proteolytically cleaved after accumulation of misfolded proteins in the ER . The active cytosolic fragment of ATF-6 translocates to the nucleus and acts as a transcription factor . A reduced expression of ATF-6 full-length protein was observed concomitant with an increase in the level of cleaved ATF-6 protein (Fig. 5a). To determine whether the artesunate-induced UPR response was limited to malignant lymphoma cells, normal B-cells from healthy PBMC donors were used. Whereas tunicamycin, an inducer of ER stress, strongly induced expression of DDIT3 and ATF-4 in normal B-cells, DDIT3 and ATF4 were not induced by artesunate in these cells, in concordance with the lack of induction of apoptosis (Figs. 2b and 5b). Taken together, these results show that artesunate induced ER stress and an UPR response leading to apoptosis selectively in malignant B-cells.
It has previously been suggested multiple cellular targets of artemisinins with involvement of reactive oxygen species (ROS) . We therefore studied the redox regulation of artesunate-induced growth suppression by blocking a key antioxidant defense pathway. Depletion of GSH using buthionine sulfoximine (BSO) triggered an increased sensitivity to artesunate in BL41 and SU-DHL-6 cells. In contrast, normal B-cells did not show a similar increased sensitivity to artesunate (Fig. 5c and Additional file 8: Figure S6). Taken together, malignant B-cells are prone to be more sensitive to artesunate when GSH levels are depleted.
Artesunate treatment decreases metabolic capacity in B-cell lymphoma cells
Given the well-established close interaction between ER and mitochondria, stress signals induced in the ER may negatively influence a variety of mitochondrial metabolic functions [24,25,26]. It has previously been reported that induction of ER stress leads to fragmentation of the mitochondrial structure and a decrease in mitochondrial membrane potential, oxygen production, and reduced respiration . The impact of artesunate on the metabolism was assessed using a seahorse metabolic analyzer (Fig. 6). Oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) of three cell lines (BL-41, SU-DHL-6 and Mino) pretreated or not with artesunate were measured over time and modulated by metabolic drugs in order to extract parameters of mitochondrial respiration and glycolysis (Fig. 6a, b). These experiments showed that artesunate had a significant and systematic diminution of all the mitochondrial respiration parameters measured, including basal respiration, ATP production, maximal respiration capacity, and the glycolysis capacity (basal glycolysis) in all cell lines tested (Fig. 6c, P < .05). The broad effects of artesunate on metabolism are likely contributing to its potent anti-lymphoma activity in vivo.
In this study, we tested small molecule inhibitors for potential tumor cell suppressive effects in B-cell lymphoma cell lines and compared with the anti-malaria drug artesunate. The screen identified artesunate as an attractive novel drug with potent activity across different cell lines, representative for various B-cell lymphoma histologies. Artesunate selectively and potently induced apoptosis in the B lymphoma cell lines but had minimal effect in normal B-cells. Of note, artesunate also showed potent anti-tumor effect in a highly aggressive lymphoma xenograft model.
The standard treatment of malaria is artemisinin based combinatorial treatment for 3 days , and the experience with longer treatment regimens is limited . In serum, artesunate is converted to dihydroartemisinin within minutes. In malaria patients receiving artesunate, the reported concentration of artesunate and dihydroartemisinin in serum has been measured to approximately 1 and 3 μM, respectively . In healthy volunteers, concentrations of dihydroartemisinin as high as 8 μg/ml (21 μM) have been measured after 10 min exposure to artesunate and were well tolerated [30, 31]. We found IC50 values for artesunate varying from 0.19 to 3.72 μM in B lymphoma cell lines. The strong apoptosis-inducing effect of artesunate at clinical relevant concentrations of artesunate is therefore promising.
Artesunate has earlier been shown to overcome bortezomib resistance in myeloma cells in vitro  and to have growth suppressive properties in primary myeloma cells and in myeloma and lymphoma cell lines by inducing apoptosis concomitant with downregulation of MYC . Others have also suggested MYC as a potential target for artesunate [33, 34], and the cell lines used for microarray analysis were chosen based on the MYC status (translocated, non-translocated, and mutated MYC in BL-41, SU-DHL-4, and Oci-Ly-2, respectively). MYC was identified as a central molecule in the network analysis with all three cell lines combined, suggesting MYC involvement, and demonstrated that artesunate had potent effects independent of MYC translocation and mutational status. Furthermore, artesunate also potently induced apoptosis in WILL-2 and Oci-Ly-18 cells, representing “double hit lymphoma,” having aberrant overexpression of MYC and BCL2, and also in U2932, with a subclone with “double hit” aberrations . This is an important observation as double hit lymphomas have dismal outcome .
The UPR was identified as the most deregulated pathway in response to artesunate. Additional top pathways activated by artesunate included tRNA charging, protein ubiquitination and amino acid biosynthesis, all connected to changes caused by UPR and ER stress [37, 38]. This suggests that the underpinning mechanism for artesunate-induced apoptosis is induction of ER stress. The UPR is a cellular adaptive response important for re-establishing protein-folding homeostasis by decreasing protein synthesis through phosphorylation of eIF2α and by increasing the ER protein-folding and degradation capacities through transcriptional activation by XBP1 and ATF6α [39,40,41]. The UPR is a sensor for ER stress and is activated upon environmental stress or other conditions resulting in accumulation of unfolded proteins, a key step in readjustment of the ER protein folding capacity to meet cellular needs [39, 42]. Importantly, the functional outcome of ER stress depends on intensity and duration, as the UPR is either pro-survival to preserve ER homeostasis or pro-death if the ER stress cannot be resolved [43, 44]. Therefore, in B lymphoma cells, artesunate might directly or indirectly increase the level of ER stress, which ultimately drives the cells into apoptosis. Here, we found artesunate to induce transcriptional upregulation of ATF-4 and DDIT3/CHOP, which also translated into increased protein expression. It has previously been shown that forced expression of ATF-4 in mouse embryo fibroblasts decreased the survival, whereas forced expression of DDIT3 had no effect, suggesting ATF-4 to be the key signal and DDIT3 secondary to ER stress induced cell death . Multiple pathways are involved in ER stress-induced apoptosis, and in most cases, these converge at the mitochondrial level . Artesunate greatly affected metabolism in the B-cell lymphoma cell lines as determined by a systematic diminution of their basal respiration, ATP-linked respiration, maximal respiration capacity, and glycolysis, although at varying degrees. The profound collapse of the maximal respiratory capacity in all cell lines upon artesunate treatment corresponds well with their decreased basal respiration and points to significant and irreversible damage to cell respiration. This is in line with previously reported mitochondrial fragmentation as a consequence of ER stress . Further evidence for this hypothesis comes from the observation that artesunate affected increased cytotoxicity in GSH-depleted malignant cells, but not in normal B-cells. The profound mitochondrial damage which is indicated by the metabolic analysis results could be associated with leaky electron transport chain which in turn may promote the aberrant generation of reactive oxygen species. Since GSH is a key intracellular antioxidant regulating ROS both in the cytosol and the mitochondria, the depletion of intracellular GSH could lead to increased vulnerability to oxidative stress .
In general, monotherapy in cancer treatment is unlikely to produce deep and lasting remission . Therefore, combining artesunate with other therapies may suppress multiple oncogenic pathways simultaneously and offer therapeutic synergy. Synergic effects with a new dimeric artemisinin compound, ART-838, with anti-leukemic drugs have recently been reported [48, 49]. ART-838 was found to have an improved bioavailability and increased half-life compared to artesunate. This new compound may therefore give new possibilities in cancer therapy. As induction of ER stress can be a trigger of immunogenic cell death as shown for some of the cytotoxic anti-cancer therapies such as anthracyclines and irradiation , artesunate might be used prior to immunotherapeutic maneuvers, such as checkpoint inhibition or T cell-based therapies . Furthermore, immunogenic signals released by dying tumor cells can induce antigen uptake as well as antigen processing and presentation by the antigen presenting cells. In fact, the presence of intratumor T cells may be an important participant in the continued anti-tumor response after the initial treatment with artesunate .
Our results indicate that artesunate has potent growth suppressive effects restricted to malignant B-cells that warrant clinical testing in B-cell lymphoma patients.
Diffuse large B-cell lymphoma
Deutche Sammlung von Mikroorganismen und Zellkulturen
Extracellular acidification rate
Carbonyl cyanide p-trifluoromethoxyphenylhydrazon
Fetal calf serum
False discovery rate
Green fluorescence protein
Ingenuity pathway analysis
In vivo imaging system
Mantle cell lymphoma
Nod SCID gamma
Oxygen consumption rate
Peripheral blood mononuclear cells
Polymerase chain reaction
Reactive oxygen species
Terminal deoxynucleotidyl transferase dUTP nick end labeling
Unfolded protein response
Swerdlow SH, Campo E, Pileri SA, Harris NL, Stein H, Siebert R, Advani R, Ghielmini M, Salles GA, Zelenetz AD, et al. The 2016 revision of the World Health Organization classification of lymphoid neoplasms. Blood. 2016;127(20):2375–90.
Sehn LH, Donaldson J, Chhanabhai M, Fitzgerald C, Gill K, Klasa R, MacPherson N, O'Reilly S, Spinelli JJ, Sutherland J, et al. Introduction of combined CHOP plus rituximab therapy dramatically improved outcome of diffuse large B-cell lymphoma in British Columbia. J Clin Oncol. 2005;23(22):5027–33.
Younes A, Ansell S, Fowler N, Wilson W, de Vos S, Seymour J, Advani R, Forero A, Morschhauser F, Kersten MJ, et al. The landscape of new drugs in lymphoma. Nat Rev Clin Oncol. 2017;14(6):335–46.
WHO: WHO model list of essential medicines. 2015.
Efferth T, Dunstan H, Sauerbrey A, Miyachi H, Chitambar CR. The anti-malarial artesunate is also active against cancer. Int J Oncol. 2001;18(4):767–73.
Efferth T, Sauerbrey A, Olbrich A, Gebhart E, Rauch P, Weber HO, Hengstler JG, Halatsch ME, Volm M, Tew KD, et al. Molecular modes of action of artesunate in tumor cell lines. Mol Pharmacol. 2003;64(2):382–94.
Huang C, Ba Q, Yue Q, Li J, Li J, Chu R, Wang H. Artemisinin rewires the protein interaction network in cancer cells: network analysis, pathway identification, and target prediction. Mol BioSyst. 2013;9(12):3091–100.
Ho WE, Peh HY, Chan TK, Wong WS. Artemisinins: pharmacological actions beyond anti-malarial. Pharmacol Ther. 2014;142(1):126–39.
Krishna S, Ganapathi S, Ster IC, Saeed ME, Cowan M, Finlayson C, Kovacsevics H, Jansen H, Kremsner PG, Efferth T, et al. A randomised, double blind, placebo-controlled pilot study of oral artesunate therapy for colorectal cancer. EBioMedicine. 2015;2(1):82–90.
Michaelsen FW, Saeed ME, Schwarzkopf J, Efferth T. Activity of Artemisia annua and artemisinin derivatives, in prostate carcinoma. Phytomedicine. 2015;22(14):1223–31.
Westfall PJ, Pitera DJ, Lenihan JR, Eng D, Woolard FX, Regentin R, Horning T, Tsuruta H, Melis DJ, Owens A, et al. Production of amorphadiene in yeast, and its conversion to dihydroartemisinic acid, precursor to the antimalarial agent artemisinin. Proc Natl Acad Sci U S A. 2012;109(3):111–8.
Kilpinen S, Ojala K, Kallioniemi O. A new method to interpret gene expression microarray data for individual tumors by profile alignment against a reference dataset of 10,000 samples. Cancer Res. 2011;71(8).
Ojala K, Kilpinen S, Kallioniemi O. A comprehensive bioinformatic analysis of the representativity of 460 cancer cell lines as model systems for 88 different clinical cancer types. Cancer Res. 2012;72(8).
Rasmussen AM, Smeland EB, Erikstein BK, Caignault L, Funderud S. A new method for detachment of Dynabeads from positively selected B lymphocytes. J Immunol Methods. 1992;146(2):195–202.
Loew R, Heinz N, Hampf M, Bujard H, Gossen M. Improved Tet-responsive promoters with minimized background expression. BMC Biotechnol. 2010;10:81.
Walchli S, Loset GA, Kumari S, Johansen JN, Yang W, Sandlie I, Olweus J. A practical approach to T-cell receptor cloning and expression. PLoS One. 2011;6(11):e27930.
Huse K, Bakkebo M, Walchli S, Oksvold MP, Hilden VI, Forfang L, Bredahl ML, Liestol K, Alizadeh AA, Smeland EB, et al. Role of Smad proteins in resistance to BMP-induced growth inhibition in B-cell lymphoma. PLoS One. 2012;7(10):e46117.
Kotecha N, Krutzik PO, Irish JM. Web-based analysis and publication of flow cytometry experiments. In: Current protocols in cytometry/editorial board, J Paul Robinson, managing editor [et al]; 2010. Chapter 10:Unit10 17.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Bakkebo M, Huse K, Hilden VI, Forfang L, Myklebust JH, Smeland EB, Oksvold MP. SARA is dispensable for functional TGF-beta signaling. FEBS Lett. 2012;586(19):3367–72.
Yoshida H, Haze K, Yanagi H, Yura T, Mori K. Identification of the cis-acting endoplasmic reticulum stress response element responsible for transcriptional induction of mammalian glucose-regulated proteins. Involvement of basic leucine zipper transcription factors. J Biol Chem. 1998;273(50):33741–9.
Haze K, Yoshida H, Yanagi H, Yura T, Mori K. Mammalian transcription factor ATF6 is synthesized as a transmembrane protein and activated by proteolysis in response to endoplasmic reticulum stress. Mol Biol Cell. 1999;10(11):3787–99.
Kavishe RA, Koenderink JB, Alifrangis M. Oxidative stress in malaria and artemisinin combination therapy: pros and cons. FEBS J. 2017;284(16):2579–91.
Bravo R, Vicencio JM, Parra V, Troncoso R, Munoz JP, Bui M, Quiroga C, Rodriguez AE, Verdejo HE, Ferreira J, et al. Increased ER-mitochondrial coupling promotes mitochondrial respiration and bioenergetics during early phases of ER stress. J Cell Sci. 2011;124(Pt 13):2143–52.
Ohta K, Okayama S, Togo A, Nakamura K. Three-dimensional organization of the endoplasmic reticulum membrane around the mitochondrial constriction site in mammalian cells revealed by using focused-ion beam tomography. Microscopy. 2014;63(Suppl 1):i34.
Friedman JR, Lackner LL, West M, DiBenedetto JR, Nunnari J, Voeltz GK. ER tubules mark sites of mitochondrial division. Science. 2011;334(6054):358–62.
Vannuvel K, Van Steenbrugge M, Demazy C, Ninane N, Fattaccioli A, Fransolet M, Renard P, Raes M, Arnould T. Effects of a sublethal and transient stress of the endoplasmic reticulum on the mitochondrial population. J Cell Physiol. 2016;231(9):1913–31.
WHO. Guidelines for the treatment of malaria. 3rd ed; 2015.
Davis TM, Phuong HL, Ilett KF, Hung NC, Batty KT, Phuong VD, Powell SM, Thien HV, Binh TQ. Pharmacokinetics and pharmacodynamics of intravenous artesunate in severe falciparum malaria. Antimicrob Agents Chemother. 2001;45(1):181–6.
Miller RS, Li Q, Cantilena LR, Leary KJ, Saviolakis GA, Melendez V, Smith B, Weina PJ. Pharmacokinetic profiles of artesunate following multiple intravenous doses of 2, 4, and 8 mg/kg in healthy volunteers: phase 1b study. Malar J. 2012;11:255.
Morris CA, Duparc S, Borghini-Fuhrer I, Jung D, Shin CS, Fleckenstein L. Review of the clinical pharmacokinetics of artesunate and its active metabolite dihydroartemisinin following intravenous, intramuscular, oral or rectal administration. Malar J. 2011;10:263.
Papanikolaou X, Johnson S, Garg T, Tian E, Tytarenko R, Zhang Q, Stein C, Barlogie B, Epstein J, Heuck C. Artesunate overcomes drug resistance in multiple myeloma by inducing mitochondrial stress and non-caspase apoptosis. Oncotarget. 2014;5(12):4118–28.
Holien T, Olsen OE, Misund K, Hella H, Waage A, Ro TB, Sundan A. Lymphoma and myeloma cells are highly sensitive to growth arrest and apoptosis induced by artesunate. Eur J Haematol. 2013;91(4):339–46.
Sertel S, Eichhorn T, Simon CH, Plinkert PK, Johnson SW, Efferth T. Pharmacogenomic identification of c-Myc/Max-regulated genes associated with cytotoxicity of artesunate towards human colon, ovarian and lung cancer cell lines. Molecules. 2010;15(4):2886–910.
Amini RM, Berglund M, Rosenquist R, Von Heideman A, Lagercrantz S, Thunberg U, Bergh J, Sundstrom C, Glimelius B, Enblad G. A novel B-cell line (U-2932) established from a patient with diffuse large B-cell lymphoma following Hodgkin lymphoma. Leuk Lymphoma. 2002;43(11):2179–89.
Sesques P, Johnson NA. Approach to the diagnosis and treatment of high-grade B-cell lymphomas with MYC and BCL2 and/or BCL6 rearrangements. Blood. 2017;129(3):280–8.
Guan BJ, Krokowski D, Majumder M, Schmotzer CL, Kimball SR, Merrick WC, Koromilas AE, Hatzoglou M. Translational control during endoplasmic reticulum stress beyond phosphorylation of the translation initiation factor eIF2alpha. J Biol Chem. 2014;289(18):12593–611.
Han J, Back SH, Hur J, Lin YH, Gildersleeve R, Shan J, Yuan CL, Krokowski D, Wang S, Hatzoglou M, et al. ER-stress-induced transcriptional regulation increases protein synthesis leading to cell death. Nat Cell Biol. 2013;15(5):481–90.
Hetz C. The unfolded protein response: controlling cell fate decisions under ER stress and beyond. Nat Rev Mol Cell Biol. 2012;13(2):89–102.
Walter P, Ron D. The unfolded protein response: from stress pathway to homeostatic regulation. Science. 2011;334(6059):1081–6.
Wang S, Kaufman RJ. The impact of the unfolded protein response on human disease. J Cell Biol. 2012;197(7):857–67.
Shaffer AL, Shapiro-Shelef M, Iwakoshi NN, Lee AH, Qian SB, Zhao H, Yu X, Yang L, Tan BK, Rosenwald A, et al. XBP1, downstream of Blimp-1, expands the secretory apparatus and other organelles, and increases protein synthesis in plasma cell differentiation. Immunity. 2004;21(1):81–93.
Szegezdi E, Logue SE, Gorman AM, Samali A. Mediators of endoplasmic reticulum stress-induced apoptosis. EMBO Rep. 2006;7(9):880–5.
Tabas I, Ron D. Integrating the mechanisms of apoptosis induced by endoplasmic reticulum stress. Nat Cell Biol. 2011;13(3):184–90.
Shore GC, Papa FR, Oakes SA. Signaling cell death from the endoplasmic reticulum stress response. Curr Opin Cell Biol. 2011;23(2):143–9.
Kramer RA, Zakher J, Kim G. Role of the glutathione redox cycle in acquired and de novo multidrug resistance. Science. 1988;241(4866):694–7.
Rodon J, Dienstmann R, Serra V, Tabernero J. Development of PI3K inhibitors: lessons learned from early clinical trials. Nat Rev Clin Oncol. 2013;10(3):143–53.
Moses BSFJ, Chen X, McCullough S, Kirkpatrick A, Kumar K, Tyner JW, Posner GH, Civin CI. Artemisinins synergize with BCL2 inhibitors by reducing MCL1 protein levels in leukemia cells. In: ASH abstract: Dec 11th 2017. Atlanta: Blood; 2017.
Fox JM, Moynihan JR, Mott BT, Mazzone JR, Anders NM, Brown PA, Rudek MA, Liu JO, Arav-Boger R, Posner GH, et al. Artemisinin-derived dimer ART-838 potently inhibited human acute leukemias, persisted in vivo, and synergized with antileukemic drugs. Oncotarget. 2016;7(6):7268–79.
Obeid M, Tesniere A, Ghiringhelli F, Fimia GM, Apetoh L, Perfettini JL, Castedo M, Mignot G, Panaretakis T, Casares N, et al. Calreticulin exposure dictates the immunogenicity of cancer cell death. Nat Med. 2007;13(1):54–61.
Kroemer G, Galluzzi L, Kepp O, Zitvogel L. Immunogenic cell death in cancer therapy. Annu Rev Immunol. 2013;31:51–72.
The authors thank Vera Hilden for help with the cell lines, and the core facilities at Oslo University Hospital: flow cytometry core facility for cell sorting, the Genomics and bioinformatics core facility for gene expression profiling, and the Department of comparative medicine for assistance with the daily care for the mice.
This work was supported by grants from the Research Council of Norway (Centres of Excellence program), the South-Eastern Norway Regional Health Authority (TKV), and the Norwegian Cancer Society (EBS).
Availability of data and materials
Microarray data is available at NCBI’s Gene Expression Omnibus with accession number GSE94553 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE94553).
The care and handling of animals for the present study were in conformity with the Norwegian Food Safety Authority in compliance with the European Convention of the Protection of Vertebrates Used for Scientific Purposes (Project ID 7729).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Overview of the additional drugs used in the first line screen. (DOCX 36 kb)
Figure S1. Artesunate induces apoptosis in B-cell lymphoma cell lines. (A) Dose response curves for the four most responsive cell lines and healthy CD19+ B-cells from Fig. 2a. Relative viability indicates relative luminescence units [RLU] from the CellTiter glo assay (n = 3-5, 0.125-8 μM artesunate, 72 h). (B) Propidium iodide (PI) staining detects a high level of dead cells in artesunate-exposed B-cell lymphoma cell lines (5 μM, P < .046, 10 μM, P < 0.05) compared to healthy B-cells (not significant (n.s)). Cells were treated with artesunate (5 and 10 μM, 72 h, n = 3) and stained with PI. The side scatter (SSC)- PI plot shows gating of the PI positive cells in untreated and artesunate treated BL-41 cells. (C) Artesunate induces apoptosis in B-cell lymphoma cell lines (P < 0.04), but not in healthy B-cells (n.s). Staining for active caspase 3 in lymphoma cell lines was done after artesunate treatment (5 μM art, 24 h, n = 3), fixation and permeabilization (PFA/methanol). The plots show example of gating for the positive active caspase 3 populations in both untreated and artesunate treated BL-41 cells. The cells are first gated on single cells, then living cells in forward/side scatter. (EPS 1906 kb)
Figure S2. No prominent cell cycle arrest was detected in B-cell lymphoma cell lines exposed to artesunate. Cells were incubated with artesunate (Art) (5 μM, 24 h, 0.5 mill/ml cells), fixed and permabilized (PFA/MeOH), and stored at − 20 °C over night. Near IR-dye (Life technology) was used as live/dead stain. Cell cycle analysis was performed using Hoechst 33,258 (2 μg/ml) incubated for at least 15 min on ice. An LSRII with UV laser was used for measurement and Dean-Jett-Fox model in Flow Jo v 7 for data analysis. A stack plot of (n = 3) cell cycle phases is shown in three different cell lines treated with artesunate vs control. T-test with correction for multiple testing with Holm-Sidak correction a statistical difference (art vs control) was only found in Mino cells in G1/S phase (P = 0.001) (Stack plot was done in Excel, and statistical analysis performed in GraphPad prism). (EPS 1049 kb)
Figure S3. The tumor load is evenly distributed among the treatment and control groups. (A) Tumor growth was monitored with IVIS luminescence measurements (physical units of surface radiance, photons/s/cm2/sr). Image was taken at day 0 when treatment was started (day 4 after injection of BL-41-luc). (B) Graphical distribution of luminescence measurements at treatment start. (C) Artesunate treatment does not affect the weight of the mice. An increased weight was observed for both groups during the treatment. Shown is weight (g) after treatment start (n = 10 for each group). Weight at start was 16.0-22.8 g and the mice were 6-10 weeks old. (EPS 3428 kb)
Table S2. The top five pathways from the ingenuity pathway analysis (IPA) show UPR as the top regulated pathway in all three cell lines with corresponding p-value and percentage overlap of genes. (DOCX 232 kb)
Figure S4. Top network for all three cell lines in IPA. These networks show interactions between proteins that have an upregulation (red) or downregulation (green) in gene expression in the three cell lines stimulated with artesunate (5 μM) for 12 h compared to control (untreated sample). (A) BL-41, (B) Oci-Ly-2 and (C) SU-DHL-4. (EPS 15844 kb)
Figure S5. IPA upstream regulators analysis predicts MYC to be inhibited. (A) In the top list of the upstream regulators analyzed by IPA, ATF4 and TRIB3 are among the top genes. MYC on the other hand is listed as number ten. The activation z-score predicts if a gene is an upstream activator (red) or an upstream inhibitor (blue). This means that the genes are predicted to regulate other genes, which are shown in B for MYC. (B) Figure shows upstream regulated network for MYC in BL-41, Oci-Ly-2 and SU-DHL-4 merged and show an interaction map between MYC and other regulated genes in the arrays. The dark blue color of MYC means IPA predicts inhibiton of MYC with “more confidence”. The goal of the Regulator Effects analysis in IPA is to provide insight into the causes and effects of differentially expressed genes or proteins in a dataset, and how they might cause differences in downstream outcomes. (EPS 1096 kb)
Figure S6. BSO treatment decreases the levels of GSH and attenuates artesunate in lymphoma cell lines. (A) GSH-Glo glutathione assay was used to measure GSH levels after treatment with BSO (50 μM, 24 h) in BL-41, SU-DHL-6 and normal CD19+ B-cells with or without IL21. (B) BSO does not decrease cell viability in lymphoma cell lines after 24 h treatment (50 μM), a small decrease in cell viability was seen in normal B-cells (n.s). (C) Cells were treated with BSO (50 μM) for 24 h before co-treatment with increasing concentration of artesunate for 24 h. Pre-treatment with BSO induced an increased sensitivity to artesunate in malignant cells, but not in normal B-cells (n = 3 (cell lines), n = 2 CD19+ cells; two different donors). (EPS 1102 kb)