Skip to main content
  • Research article
  • Open access
  • Published:

Frankincense and myrrh and their bioactive compounds ameliorate the multiple myeloma through regulation of metabolome profiling and JAK/STAT signaling pathway based on U266 cells



Frankincense and myrrh are used as traditional anti-inflammatory and analgesic medicines in China. It has been reported that frankincense and myrrh have significant anti-tumor activities. The present study was designed to investigate the inhibitory efficacy of frankincense ethanol extracts (RXC), myrrh ethanol extracts (MYC), frankincense -myrrh ethanol extracts (YDC), frankincense -myrrh water extracts (YDS) and their main compounds on U266 human multiple myeloma cell line.


The inhibition effects of cell proliferation was evaluated by MTT assays. Cell culture supernatant was collected for estimation of cytokines. Western blot analysis was designed to investigate the regulatory of JAK/STAT signal pathway. In addition, cell metabolomics based on the ultra-performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry (UPLC/Q-TOF-MS) had been established to investigate the holistic efficacy of frankincense and myrrh on U266 cells. Acquired data were processed by partial least-squares discriminant analysis (PLS-DA) and orthogonal projection to latent structures squares-discriminant analysis (OPLS-DA) to identify potential biomarkers.


RXC, MYC significantly inhibited the proliferation of U266 cells at dose of 25–400 μg/mL, YDC and YDS at the dose of 12.5–400 μg/mL. 3-O-acetyl-α-boswellic acid, 3-acetyl-11 keto-boswellic acid and 11-keto-boswellic acid had the most significant anti- multiple myeloma activities in the 10 compounds investigated, therefore these 3 compounds were selected as representatives for Elisa assay and western blotting experiments. All the extracts and active compounds ameliorated the secretion of cytokines and down-regulated the expression of JAK/STAT signaling pathway-related proteins. Comparing RXC, MYC, YDC and YDS-treated U266 cells with vehicle control (DMSO), 13, 8, 7, 7 distinct metabolites and 2, 2, 3, 0 metabolic target pathways involved in amino acid metabolism, lipid metabolism, vitamin metabolism, arachidonic acid were identified, respectively.


Taken together our results suggest that the frankincense and myrrh and their bioactive compounds inhibit proliferation of U266 multiple myeloma cells by regulating JAK/STAT signaling pathway and cellular metabolic profile.

Peer Review reports


Today, multiple myeloma (MM) is the second most common cancer in hematological malignant disease. More than 12 thousand deaths were estimated in 2018, accounted for 2.1% of all cancer deaths and the percent surviving 5 years was 50.7% in 2008–2014, according to epidemiology, surveillance, and end result program ( [1]. The overall efficacy of current conventional therapies for MM, including combined chemotherapy, radiation therapy, stem cell transplantation, new drug treatment, etc., are not ideal. However, the experimental researches on therapeutic effect and mechanism of traditional Chinese medicine on MM have been made a breakthrough.

Frankincense and myrrh, classical traditional Chinese medicine, are commonly used in Chinese medicine as anti-inflammatory and analgesia drugs. Vitro studies have shown that frankincense, which has been commonly believed the efficacy substance is boswellic acid, have anti-inflammatory activity [2], analgesic [3], immunosuppressive and antitumor activities [4]. The main components in myrrh are volatile oils, including monoterpenoid and sesquiterpenoids, having the same effect as frankincense [5,6,7,8,9]. Previous studies have exhibited that the compatibility of frankincense and myrrh have certain impacts against tumors, including Ehrlich ascites [10], prostate cancer cells [11], and breast cancer cells [6], et al. However, there are few studies about the evaluation of the effects on MM and mechanisms.

Owing to the limitation, such as poor repeatability, less quantity of data, and time-consuming of conventional anti-tumor drugs research models, including MTT reduction assay for cellular proliferation, flow cytometry for cellular apoptosis, RT-PCR for gene expression, and Western blot for protein expression [12], the metabolomics has been more and more applied in anti-tumor mechanisms. Metabolomics is one of the newest biology approaches focused on the global metabolic status of the entire organism through non-targeted analysis of metabolites in biological samples and based on the analytical techniques including mass spectroscopy, 1H-NMR spectroscopy and liquid chromatography-mass spectroscopy [13]. Notably, cancer is a disease with mitochondrial energy metabolic disorder, which is well recognized as Warburg effect. Therefore, metabolic profiling of cells underlying treatment with anticancer drug candidates is greatly helpful for understanding their action mechanisms.

U266 is a human multiple myeloma cell, and belongs to a suspension cell line, which is the first choice for the study of multiple myeloma diseases in vitro. In this study, we designed to investigate the inhibitory efficacy of Frankincense, Myrrh and their bioactive based on U266 multiple myeloma cells and explain the action mechanism of the anti-multiple myeloma of frankincense and myrrh combining conventional molecular biology with metabolomics.


Chemicals and instruments

U266 human multiple myeloma cells were purchased from the American Type Culture Collection (VA, Catalog No.: CC-Y1527, USA). Fetal bovine serum and RPMI-1640 medium were purchased from Gibco Life Technologies (NY, USA). Methylene blue and dimethyl sulfoxide (DMSO) were purchased from Sigma-Aldrich (MO, USA). IL-6, VEGF ELISA kits were purchased from Nanjing Jiancheng Biotechnology Co., Ltd. (Nanjing, China). JAK1, p-JAK1 STAT3, p-STAT3 and GAPDH antibodies were purchased from Cell Signaling Technology (BSN, USA).

Waters Acquity TM Ultra Performance LC system (Waters, Milford, MA, USA) equipped with a Quattro Micro MS spectrometer and a Waters Xevo TM G2 QTof MS (Waters MS Technologies, Manchester, NH, USA). Deionized water was purified on a Milli-Q system (Millipore, Bedford, MA, USA). Mass Lynx v4.1 workstation was adopted to analyze the data, and Ultra-highspeed centrifuge at low temperature (Thermo Scientific, Waltham, MA, USA); DMI3000M microscope (Leica, München, Germany) were used.

Frankincense (lot number 171010) and myrrh (lot number 171108) were purchased from Suzhou Tianling Chinese Medicine Pieces Co., Ltd. (Suzhou, China). Standards 3-O-acetyl--boswellic acid (lot number 20150515), 3α-acetyl-20(29)-lupene-24-oic acid (lot number 20150513), β-boswellic acid (lot number 20150504), acetyl 11α-methoxy-β-boswellic acid (lot number 20150519), 3α-acetoxy-tirucall-7,24-dien-21-oic acid (lot number 20150523), 11-keto-β-boswellic acid (lot number 20150525), 3-acetyl-β-boswellic acid (lot number 20150410) ware purchased from Baoji Chenguang Biotechnology Co., Ltd. (Baoji, China) and the purity was greater than 98%. Abietic acid and 2-methoxy-8, 11-epoxygermacra-1(10)-7, 11-trien-6-one were laboratory-made and the purity are greater than 98% determined by UPLC-MS.

Preparation of drug solution

Ethanol extracts preparation: frankincense, myrrh and frankincense-myrrh pair (1:1 compatibility) were extracted twice with 12 times amount of 95% ethanol, 2 h each time. Combined the two supernatants, then evaporated the solvent to obtain each of the prepared samples. Each sample was separately dissolved in 1% DMSO to prepare 4 mg/mL of the RXC, MYC, and YDC technical liquid concentrates for storage.

Water extract preparation: frankincense-myrrh pair (1:1 compatibility) was extracted twice with 12 times the amount of water, 2 h each time. Combined the two supernatants, then evaporated the solvent to obtain each of the prepared samples. The sample was dissolved in water to prepare 4 mg/mL YDS technical liquid concentrates for storage.

Compounds monomer solution preparation: Each compound monomer was dissolved in 1%DMSO to prepare 1 mmol/L compound monomer technical liquid concentrates for storage.

Anti-proliferative activity

U266 cells were cultured in RPMI 1640 medium containing 10% FBS and 1% antibiotic-antimycotic and maintained in an incubator at 5% CO2 and 37 °C as ATCC described. A total of 2 × 104 cancer cells in growth media were placed in each well of a 96-well plates. After that, various concentrations of RXC, MYC, YDC, YDS (25, 50, 100, 200 and 400 μg/L) and compounds monomer (5, 10, 25, 50 and 100 μmol/L) were added in each well, then the cells were cultured for 48 h. Each of the drug groups was evaluated in triplicate. After incubation, 10 μL of 0.5 mg/mL MTT solution was added to each well in the dark and cells were incubated for another 4 h. The supernatant was discarded and 150 μL of DMSO was added. Then shaking slightly for 30 min until blue-violet crystals completely dissolved. The absorbance (A) value of each well was measured at 570 nm by a microplate reader and the cell proliferation inhibition rate was calculated. Each experiment was repeated at least three times. Proliferation inhibition rate (%) = (1 - experimental group A value / control group A value) × 100%.

ELISA assay

IL-6, VEGF kits assay were used to evaluate the change of cytokine. According to the manufacturer’s protocols, 2 × 104 cells were placed in each well of a 96-well plates. After treated with 200, 100, 50 μg/mL RXC, MYC, YDC, YDS and 100, 50, 25 μmol/L BC, AKBA, KBA for 48 h respectively, cell supernatant was collected for Elisa detection.

Western blot analysis

The myeloma cells of each group were collected after U266 cells were treated with 200 μg/mL of each extracts and 100 μmol/L of monomers for 48 h. The cells were washed with PBS, and extracted with 300 μL RIPA protein extraction lysate for 30 min in an ice bath. The supernatant was collected by centrifuged at 13,000 rpm at 4 °C for 10 min. Protein extracts were separated with SDS-PAGE and transferred to PVDF membrane. After blocking, p-JAK1, JAK1, p-STAT3, STAT3 and GAPDH antibodies were added and incubated. Membranes were washed and incubated in TBST and HRP-conjugated secondary antibody, respectively. The signals were detected using the Bio-Rad ChemiDoc XRS+.

UPLC/Q-TOF/MS analysis of metabolic Profling

Sample preparation for LC-MS

1 × 107 U266 cells were seeded in a 10 cm dish and exposed to 100 μg/mL of RXC, MYC, YDC or YDS or an equal amount of DMSO as a control. Six replicates in separate dishes for each group were analyzed. After 48 h incubation, cells were washed with PBS, then centrifuged and harvested. The cell pellets were rapid quenched with liquid nitrogen, then immediately dissolved in 1.0 mL mixture of methanol/water in ratio of 4:1 (v/v) at − 20 °C and ultra-sonicated for 30 min in an ice bath ultra-sonicator and subsequently centrifuged to collect the supernatant. The supernatant was dried, then resuspended with 100 μL of 80% methanol and filtered through 0.22 mm mesh millipore filters. Drawing 80 μL of supernatant for LC-MS analysis. In parallel a quality control (QC) sample was prepared by mixing the leftover 20 μL of supernatant from each of the 30 samples.

UPLC-Q TOF/MS conditions

Chromatographic conditions: Acquity™ UPLC BEH C18 column (2.1 mm × 50 mm, 1.7 μm); column temperature: 35 °C; flow rate 0.4 mL/min; injection volume: 3 μL; mobile phase: 0.1% formic acid water (A)-acetonitrile (B). Gradient elution conditions: 0~1 min: 5% B, 1~3 min: 5 ~ 50% B, 3~13 min: 50 ~ 85% B, 13~14 min: 85 ~ 95% B.

Mass spectrometry conditions: samples were ionized using an ESI ion source (ESI+/ESI) with a mass scan range of 100 to 1000 m/z; capillary voltage: 3.0 kV; cone voltage: 30 V; extraction voltage of 2.0 V; and ion source temperature of 120 °C; desolvation gas temperature 350 °C; cone flow gas flow 50 L/h; collision energy: 20 ~ 50 eV; desolvation gas flow 600 L/h; activation time: 30 ms; collision gas: high purity nitrogen. A solution of 200 pg/mL leucine- enkephalin (ESI+: 556.2771 m/z, ESI: 555.2615 m/z) was used as the locking mass solution to ensure the accuracy during MS analysis.

Metabolomics data processing and analysis

The raw data acquired are processed using Mass Lynx software. After chromatographic peak identification, alignment and normalization, a data matrix was obtained: the mass scan range is 100–1000 Da, the retention time range is 0–15 min, the mass deviation threshold is 0.01 Da, the mass window is 0.05 Da, the retention time window is 0.2 min, and the relative peak intensity threshold is 5%, the noise cancellation level is 6.0.

In addition, based on the impurity peaks appearing in the blank solvent control sample, a list of excluded ions was prepared and introduced into the method. The final process will result in a list of data consisting of retention time, m/z values and normalized peak areas. The data was imported into the EZ info 2.0 module for supervised partial least squares discriminant analysis (PLS-DA) and orthogonal partial least squares discriminant analysis (OPLS-DA). The Score-Plot map obtained by OPLS-DA analysis judges and displays the degree of dispersion of the metabolite profiles of samples between groups and extracts potential biomarkers. A point in OPLS-DA represents a variable. The variable importance in the projection (VIP) is measured by the value, and the variable is screened according to the VIP value. The difference between the groups was generally the normal control and the model group. This study directly compared the control group with each drug treated group.

Potential biomarkers identification and metabolic pathway analysis

Metabolites with VIP > 1, and those with P < 0.05 were compared between the control group compared with each drug treated group as potential biomarkers. They were input into the HMDB (, firstly, the molecular weight (deviation 0.5) was selected, the non-endogenous was excluded, and then the relationship with the human cell metabolism was selected, and the unrelated to multiple myeloma was further deleted to obtain the difference metabolite. After the metabolites were identified, information such as the metabolite index number, molecular formula, molecular weight, retention time, VIP value, and relative expression amount were compiled into a table, and introduced into the Metabo database ( for metabolic pathway construction.

Statistical analysis

Data were analyzed using SPSS version 19.0. All values were expressed as the means ± SD of at least three independently performed experiments. Statistical analyses were conducted with Student’s t-test. Differences with P < 0.05 were considered to be statistically significant. and *P < 0.05, **P < 0.01, ***P < 0.001 were compared between the treatment group and the control.


Determination of chemical components in frankincense-myrrh

Before investigating the anti-tumor activities of frankincense and myrrh, we analyzed the chemical components in the two drugs with the YDC and YDS as representatives. (n = 3). Ten components with abundant contents were detected out, and the contents in YDC was much higher in that in YDS (Table 1).

Table 1 The contents of 10 components in YDC and YDS

Anti-proliferative activities of frankincense and myrrh on U266 cells

RXC, MYC, at the doses of 25–400 μg/mL, and YDC, YDS, at the doses of 12.5–400 μg/mL, showed dose-dependent anti-proliferative effects on U266 cell growth (Fig. 1). The IC50 values of RXC, MYC, YDC, YDS in inhibiting the growth of U266 were 194.81, 229.07, 140.68, 316.75 μg/mL, respectively.

Fig. 1
figure 1

Anti-proliferative effects of RXC, MYC, YDC and YDS on U266 cells. Cell viability of U266 cells analyzed by MTT assay. P < 0.05 compared to the control group cells when the concentration is higher than 25 μg/mL

Evaluation of chemical components bioactivity

Ten components had different degrees of inhibition effects on U266 cells proliferation. Compared with other components, 3-O-acetyl-α-boswellic acid (BC), 3-acetyl-11 keto-boswellic acid (AKBA) and 11-keto-boswellic acid (KBA) had the most significant activities, therefore the three components were selected as representatives for Elisa assay and western blotting experiments (Table 2). The IC50 values of BC, AKBA, KBA were 35.02, 67.79, 86.56 μmol/L, respectively.

Table 2 Inhibition of the proliferation of U266 cells by the 10 chemical components

Effect of frankincense and myrrh on IL-6 and VEGF secretion of U266 cells

A significant reduction in the secretion levels of IL-6 was observed in MYC group compared to the control group, while no significant changes of the secretion levels of IL-6 was observed in other extract groups (Fig. 2a). The results were consistent with the results of BC, AKBA, KBA (Fig. 2b). All the four extracts and active components significantly inhibited the secretion of VEGF in U266 cells (Fig. 3).

Fig. 2
figure 2

The level of IL-6 secreted by U266 after be treated with extracts or active components. a RXC, MYC, YDC, YDS (50, 100 and 200 μg/mL); b BC, AKBA, KBA (25, 50 and 100 μmol/L) for 48 h.*P < 0.05, **P < 0.01, vs control group(x ± s, n = 3)

Fig. 3
figure 3

The level of VEGF secreted by U266 after be treated with extracts or active components. a RXC, MYC, YDC, YDS (50, 100 and 200 μg/mL); b BC, AKBA, KBA (25, 50 and 100 μmol/L) for 48 h.*P < 0.05, **P < 0.01, vs control group(x ± s, n = 3)

Expression of JAK1, P-JAK1, STAT3 and P-STAT3

JAK/STAT signal pathway is sustained activation status in control group cells because of the high levels expression of p-JAK1, JAK1, p-STAT3 and STAT3 proteins. Compared with control group, the expression levels of each protein were reduced to varying degrees in and U266 cells after be treated with 200 μg/mL RXC, MYC, YDC, YDS or 100 μmol/L BC, AKBA, KBA. Among them, AKBA had the best inhibitory effect on p-JAK1, KBA had the best therapeutic effect on JAK1, YDS had the most significant effect on p-STAT3, and MYC had a better effect on STAT3 than other treatment groups. Their relative reductions were 0.87, 0.60, 0.82, and 0.40, respectively. (Fig. 4).

Fig. 4
figure 4

The effect of RXC, MYC, YDC, YDS and BC, AKBA, KBA on JAK/STAT signaling pathway.*P < 0.05, **P < 0.01, vs control group(x ± s, n = 3)

Metabolomics results

Metabolomics analysis of RXC, MYC, YDC and YDS-treated cells

As shown in Fig. 5a, the QC samples are tightly clustered together in PCA scores plot, indicating that there was column stability in the whole run. While there had no obvious variation of QC samples over all observations with respect to run order in the trend plot (Fig. 5b). Ten ions chromatographic peaks were selected to evaluated the repeatability of method through six replicates of QC samples. The relative standard deviations (RSD) of ion intensity indicated that the repeatability of method established were well (Table 3).

Fig. 5
figure 5

Assessment of QC samples. a PCA score plot (PC1 versus PC2) of test samples and QC samples; b Trend plot showing the variation of t [1] over all observations. QC samples were colored as red boxes and test samples were colored as black triangle. X axis numbers represented sample number (35 injections); Y axis was arbitrary

Table 3 Variation condition of ion intensity of selected 10 ions present in the QC samples

Control and RXC, YDC treated groups were obviously separated along the first principal component, while control and MYC treated groups were significantly separated along the second principal component, and YDS treated group were not significantly different from the control group (Fig. 6).

Fig. 6
figure 6

PCA scores plot resulting from UPLC/MS spectra of U266 cells for control and treatment groups. a PCA scores plot in positive ion mode. b PCA scores plot in negative ion mode

Identification of potential metabolic biomarkers

The differential metabolites were screened according to the S-plot, VIP value, and inter-group peak area t-test between the control group and the treated groups obtained by the UPLC-QTOF/MS analysis platform. These differential metabolites identification were conducted with molecular weight, distribution and polarity of metabolite species by the Human Metabolome Database ( In this study, thirteen, eight, seven and seven endogenous metabolites involving in multiple metabolic pathways, such as amino acid metabolism, glucose metabolism, lipid metabolism, and vitamin metabolism etc., were ultimately identified which showed significantly changes between control and RXC, MYC, YDC, YDS groups, respectively (P* < 0.05 or P** < 0.01 or P*** < 0.001) (Table 4).

Table 4 Potential metabolites selected and identified between treated groups and control group

Metabolic pathway analysis

A list of metabolite identifications was imported into the MetPA database, and KEGG database was searched to explore potential metabolic pathways affected by the treatment of RXC, MYC, YDC and YDS. The pathways that with an impact value above 0.10 were screened out as the potential target pathway. RXC group mainly affected the retinol metabolic pathway and the arachidonic acid metabolic pathway (Fig. 7a). MYC group mainly affected the riboflavin metabolic pathway and the phenylalanine metabolic pathway (Fig. 7b). The metabolites affected by YDC group were involved in riboflavin metabolism pathway, sphingolipid metabolism pathway, phenylalanine metabolism pathway (Fig. 7c), however there was no potential target pathway revealed in the YDS group (Fig. 7d).

Fig. 7
figure 7

Metabolic pathways constructed by differential metabolites in the RXC, MYC, YDC and YDS group. a) - d represented RXC, MYC, YDC and YDS groups, respectively. a Arachidonic acid metabolism, b Retinol metabolism, c Phenylalanine metabolism, d Riboflavin metabolism, e Sphingolipid metabolism


The aim of the present study was to assess the inhibitory effects of frankincense and myrrh, and try to elucidate the anti- myeloma mechanism through metabolic profiling. In this study, RXC, MYC, YDC and YDS showed dose-dependent anti-proliferative effects on U266 cell growth at the doses of 25–400 μg/mL. BC, AKBA and KBA had the most significant activities, therefore the three components were selected as representatives for Elisa assay and western blotting experiments. Studies have shown that these three components have significant antitumor activity. Park B et al. reported that AKBA concentrations exceeding 50 μmol/L in U266 cells down-regulated CXCR4 expression, which is associated with the invasion and metastasis of cancer cells [14]. Zhao W et al. showed that BC-4 blocked the invasion and metastasis of B16F10 mouse melanoma cells by inducing differentiation and blocking the cell population in G1 phase and inhibiting topoisomerase II activity [15].

MM is a tumor that occurs in the terminal stage of B cell differentiation. Some cytokine abnormalities can affect the activation, development and differentiation of normal B cells, and play a key role in the occurrence and development of MM [16]. IL-6 is the most important cytokine to maintain MM cell survival and promote its proliferation. Tsuyama N et al. confirmed that MM cells expressing IL-6 have higher malignancy, faster proliferation, and are prone to drug resistance [17]. IL-6 regulates the functional status of tumor cells mainly through two signaling pathways, namely, the Ras-dependent MAPK pathway, that may be associated with IL-6 proliferation of MM cells, and the Ras-independent STATs pathway, that may be associated with IL-6 inhibition of apoptosis in MM cells [18]. VEGF is a highly specific pro-vascular endothelial growth factor that promotes increased vascular permeability, extracellular matrix degeneration, vascular endothelial cell migration, proliferation, and angiogenesis [19]. MM cells promote the expression of IL-6 in BMSCs of bone marrow stromal cells by expressing VEGF and its receptors Flt-1 and KDR, thereby promoting the proliferation of MM cells and inhibiting apoptosis [20]. In this study, MYC significant reduced in the secretion levels of IL-6 and all the four extracts and active components significantly inhibited the secretion of VEGF in U266 cells.

We further investigated the expression levels of JAK1, P-JAK1, STAT3 and P-STAT3 to prove the cell apoptosis triggered by RXC, MYC, YDC and YDS. JAK/STAT pathway is one of the most common tumor cells apoptotic signaling pathways involving numerous cytokines, growth factors and hormones [21]. Sustained activation of JAK/STAT signaling is frequently linked to cancer cell proliferation, survival, metastasis, tumor immunosuppression, and angiogenesis, and IL-6-mediated apoptosis inhibitory [22,23,24,25]. Kunnumakkara AB et al. evaluated the regulation of AKBA in U266 and MM.1S human multiple myeloma cells, SCC4 human oral squamous cells and A293 cells in vitro. The data showed that AKBA suppresses IL-6-induced STAT3 activation on those cells and phosphorylation of both Jak2 and Src. In addition, the significant inhibition of STAT3 activation by AKBA led to the suppression of gene products involved in proliferation (cyclin D1), survival (Bcl-2, Bcl-xL and Mcl-1), and angiogenesis (VEGF) [26]. Using p-STAT3, the final link of the pathway, as a standard to evaluate the regulation effects of JAK/STAT signaling pathway by treated groups, the results showed there were significant decrease in each treated groups (P < 0.01), however, there was no significant difference between the extract groups.

Meanwhile, we investigated the metabolic profile of U266 cells regulated by RXC, MYC, YDC and YDS. Control and RXC, YDC treated groups were obviously separated along the first principal component, that indicated the cellular metabolic profile were significantly altered under the treatment of RXC and YDC. The results were consistent with anti-proliferative activity. The differential metabolites were screened by the UPLC-QTOF/MS analysis platform. In this study, thirteen, eight, seven and seven endogenous metabolites showed significantly changes between control and RXC, MYC, YDC, YDS groups, respectively. Because the metabolites in RXC and YDC treated groups were the most significantly changed, further studies were focused on the metabolites L-homoserine, γ-glutamine putrescine and LysoPC (18:0), and arachidonic acid metabolism, retinol metabolism, phenylalanine metabolism, riboflavin metabolism, sphingolipid metabolism and other pathways. Which altered obviously in both RXC and YDC treated groups.

Lyso-PC, derived from glycerophospholipid metabolic pathway, has been confirmed to increase oxidative stress by inducing antibody formation and stimulating macrophages, leading to lipid peroxidation and vascular endothelial dysfunction in inflammatory state [27]. Moreover, LysoPC is an intermediate product of phospholipid choline (PC) metabolism, and PC is the main lipid component of biofilm, so the change of LysoPC reflects the state of cell membrane. Compared with the control group, the content of LysoPC (18:0) was significantly increased in RXC and YDC groups, indicating that RXC and YDC can induce oxidative stress damage and destroy the stability of the cell membrane [28].

γ-glutamine putrescine is involved in the degradation pathway of putrescine II, which is formed by putrescine and L-glutamic acid driven by ATP. Putrescine is important for maintaining the basic function and viability of cells. Recent studies have found that excessive accumulation of putrescine in cells can lead to increased uptake of the medium or increased intracellular synthesis, leading to apoptosis [29, 30]. In our study, the levels of γ-glutamine putrescine in RXC and YDC groups were significantly lower than those in the control group, indicating that RXC and YDC may block the degradation pathway of putrescine in U266 cells, and then the accumulation of γ-glutamine putrescine induce apoptosis of tumor cells.

L-homoserine is an intermediate product of methionine biosynthesis, involved in the methionine metabolic pathway. Methionine is an essential amino acid that plays an important role in protein synthesis and other biochemical processes. Studies have shown that methionine can be participate in the methylation of DNA and proteins as a methylation donor, thereby regulating the expression of genes and proteins [31]. In addition, L-methionine can inhibit the proliferation of cancer cells by inhibiting post-translational modification of the gene p53 [32]. In the present study, the L-homoserine in RXC and YDC groups were at low levels, indicating that these two may increase the activity of methionine synthesis-related enzymes, promote the conversion of L-homoserine to methionine, and further cause methylation of U266 cell DNA, thereby inhibiting cell proliferation.

Epoxyeicosatrienoic acids (EETs) are epoxy derivative of arachidonic acid formed by cytochrome P-450 (CYP) cyclooxygenase. Cheranov et al. [33] reported that EETs induced endothelial angiogenesis by stimulating VEGF expression through Src-dependent STAT-3. Research shows 17-ODYA is an effective cycysase inhibitor, which can inhibit the expression of EETs and thus reduce the production of VEGF [34]. So, RXC and YDC may inhibition of EETS secretion on cytochrome P-450 pathway promotes apoptosis of U266 cells.

Through the above pharmacodynamic research and preliminary discussion on signal pathway and metabolic pathway, we speculate that frankincense and myrrh first inhibit the secretion of IL-6 in myeloma cells, that may affect the amount of IL-6 in the bone marrow microenvironment. The decrease in the amount of IL-6 causes an activation of JAK/STAT signaling pathway. Meanwhile, the stimulation of frankincense and myrrh result in the changes of metabolic pathway, thereby inhibiting tumor angiogenesis, inducing apoptosis, causing cell necrosis. The potential mechanism of frankincense and myrrh against multiple myeloma is shown in Fig. 8.

Fig. 8
figure 8

The potential mechanism of frankincense and myrrh against MM


In summary, all of the treatment with frankincense ethanol extracts, myrrh ethanol extracts, frankincense-myrrh ethanol extracts and frankincense-myrrh water extracts showed significant anti-proliferation effect on U266 cells at the dose of 25–400 μg/mL. The inhibition effects on U266 cells proliferation of all extracts and active components were associated with the suppressed activation of JAK/STAT signaling pathways and regulation of the bone marrow microenvironment. In addition, we speculate that the proliferation inhibitory effects of U266 cells may also owe to the regulation of vitamin metabolism, amino acid metabolism and lipid metabolism. This study provides a basis for studying the mechanism of the compatibility of frankincense myrrh from the perspective of endogenous metabolites. However, the specific mechanism of action remains to be determined, and there is a lack of normal control of PBMC cells and validation of metabolites in current studies. We believe that PBMC/ in vivo toxicity and metabolite validation are necessary, and we will conduct validation tests for those in the future.

Availability of data and materials

Data are all contained within the manuscript.



3-acetyl-11 keto-boswellic acid


3-O-acetyl-α-boswellic acid




Janus kinase/signal transducer and activator of transcription


11-keto-boswellic acid


multiple myeloma


3-(4,5-dimethylthizaol-2-yl)-2,5-diphenyl tetrazolium bromide


Myrrh ethanol extracts


Frankincense ethanol extracts


Vascular endothelial growth factor


Frankincense-myrrh ethanol extracts


Frankincense-myrrh water extracts


  1. Szudy-Szczyrek A, Szczyrek M, Soroka-Wojtaszko M, Hus M. New prognostic biomarkers in multiple myeloma. Postepy Hig Med Dosw. 2016;70:811–9.

    Article  Google Scholar 

  2. Su S, Duan J, Chen T, Huang X, Shang E, Yu L, Wei K, Zhu Y, Guo J, Guo S. Frankincense and myrrh suppress inflammation via regulation of the metabolic profiling and the MAPK signaling pathway. Sci Rep. 2015;5(1):13668.

    Article  Google Scholar 

  3. Hu D, Wang C, Li F, Su S, Yang N, Yang Y, Zhu C, Shi H, Yu L, Geng X. A combined water extract of frankincense and myrrh alleviates neuropathic pain in mice via modulation of TRPV1. Neural Plast. 2017;2017:1–11.

    Article  CAS  Google Scholar 

  4. Ruiqi W, Yan W, Zuhua G, Xianjun Q. The comparative study of acetyl-11-keto-beta-boswellic acid (AKBA) and aspirin in the prevention of intestinal adenomatous polyposis in APC (min/+) mice. Drug Discov Ther. 2014;8(1):25–32.

    Article  Google Scholar 

  5. Mostafa Abbas S, Ashraf Abd-Elkhalik H. Analgesic, anti-inflammatory and anti-hyperlipidemic activities of Commiphora molmol extract (myrrh). J Intercult Ethnopharmacol. 2014;3(2):56–62.

    Article  Google Scholar 

  6. Chen Y, Zhou C, Ge Z, Liu Y, Liu Y, Feng W, Li S, Chen G, Wei T. Composition and potential anticancer activities of essential oils obtained from myrrh and frankincense. Oncol Lett. 2013;6(4):1140–6.

    Article  CAS  Google Scholar 

  7. Su S, Duan J, Tang Y, Zhang X, Yu L, Jiang F, Zhou W, Luo D, Ding A. Isolation and biological activities of neomyrrhaol and other terpenes from the resin of Commiphora myrrha. Planta Med. 2009;75(04):351–5..

  8. Rapper S, De VSF, Van KGPP, Viljoen AM, Dagne E. The additive and synergistic antimicrobial effects of select frankincense and myrrh oils--a combination from the pharaonic pharmacopoeia. Lett Appl Microbiol. 2012;54(4):352–8.

    Article  Google Scholar 

  9. Baral S, Cho DH, Pariyar R, Yoon CS, Chang BY, Kim DS, Cho HK, Kim SY, Oh H, Kim YC. The ameliorating effect of myrrh on scopolamine-induced memory impairments in mice. Evid Based Complement Alternat Med. 2015;2015:1–9.

  10. Agrawal SS, Sarita S, Rajani M, Maneesha P. Antitumor properties of Boswellic acid against Ehrlich ascites cells bearing mouse. Food Chem Toxicol. 2011;49(9):1924–34.

    Article  CAS  Google Scholar 

  11. Xiao-Ling W, Feng K, Tao S, Young CYF, Hong-Xiang L, Hui-Qing Y. Sesquiterpenoids from myrrh inhibit androgen receptor expression and function in human prostate cancer cells. Acta Pharmacol Sin. 2011;32(3):338–44.

    Article  Google Scholar 

  12. Yini W, Dan G, Zhe C, Shangfu L, Chunmei G, Deliang C, Feng L, Hongxia L, Yuyang J. Acridone derivative 8a induces oxidative stress-mediated apoptosis in CCRF-CEM leukemia cells: application of metabolomics in mechanistic studies of antitumor agents. PLoS One. 2013;8(5):e63572.

    Article  Google Scholar 

  13. Nicholson JK, Lindon JC, Holmes E. 'Metabonomics': understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data. Xenobiotica. 1999;29(11):1181–9.

    Article  CAS  Google Scholar 

  14. Byoungduck P, Bokyung S, Yadav VR, Sung-Gook C, Mingyao L, Aggarwal BB. Acetyl-11-keto-β-boswellic acid suppresses invasion of pancreatic cancer cells through the downregulation of CXCR4 chemokine receptor expression. Int J Cancer. 2011;129(1):23–33.

    Article  Google Scholar 

  15. Zhao W, Entschladen F, Liu H, Niggemann B, Fang Q, Zaenker KS, Han R. Boswellic acid acetate induces differentiation and apoptosis in highly metastatic melanoma and fibrosarcoma cells. Cancer Detect Prev. 2003;27(1):67–75.

    Article  CAS  Google Scholar 

  16. None: PP2–038 A review of the cytokine network in multiple myeloma: Diagnostic, prognostic and therapeutic implications. 2009, 48(1–2):1. (DOI:

  17. Abroun S, Liu S, Tsuyama N, Otsuyama KI, Ishikawa H, Zheng X, Obata M, Taniguchi O, Kawano MM, Li FJ. The Regulatory Mechanism of IL-6-dependent Proliferation of Human Myeloma Cells[J]. Nephron Clinical Practice. 2003;8(6):409–11.

  18. Puthier D, Bataille R, Amiot M. IL-6 up-regulates mcl-1 in human myeloma cells through JAK / STAT rather than ras / MAP kinase pathway. Eur J Immunol. 2015;29(12):3945–50.

    Article  Google Scholar 

  19. Anderson K. The vascular endothelial growth factor receptor tyrosine kinase inhibitor PTK787/ZK222584 inhibits growth and migration of multiple myeloma cells in the bone marrow microenvironment. Eur J Haematol. 2003;70(4):266–7.

    Article  Google Scholar 

  20. Klaus P, Anderson KC. The pathophysiologic role of VEGF in hematologic malignancies: therapeutic implications. Blood. 2005;105(4):1383–95.

    Article  Google Scholar 

  21. Kisseleva T, Bhattacharya S, Braunstein J, Schindler CW. Signaling through the JAK/STAT pathway, recent advances and future challenges. Gene. 2002;285(1):1–24.

    Article  CAS  Google Scholar 

  22. Taniguchi K, Karin M. IL-6 and related cytokines as the critical lynchpins between inflammation and cancer. Semin Immunol. 2014;26(1):54–74.

    Article  CAS  Google Scholar 

  23. Hua Y, Drew P, Richard J. STATs in cancer inflammation and immunity: a leading role for STAT3. Nat Rev Cancer. 2009;9(11):798–809.

    Article  Google Scholar 

  24. Qing C, Eirini B, Pasquale S, Marjan B, Sizhi Paul G, Laura D, Jared W, Till T, Selena G, Xinmin Z. The IL-6/JAK/Stat3 feed-forward loop drives tumorigenesis and metastasis. Neoplasia. 2013;15(7):848–IN45.

    Article  CAS  Google Scholar 

  25. Al-Zaid-Siddiquee T. STAT3 as a target for inducing apoptosis in solid and hematological tumors. Cell Res. 2008;18(2):254–67.

    Article  CAS  Google Scholar 

  26. Kunnumakkara AB, Nair AS, Sung B, Pandey MK, Aggarwal BB. Boswellic acid blocks signal transducers and activators of transcription 3 signaling, proliferation, and survival of multiple myeloma via the protein tyrosine phosphatase SHP-1. Mol Cancer Res. 2009;7(1):118.

    Article  CAS  Google Scholar 

  27. Catala A. Lipid peroxidation of membrane phospholipids generates hydroxy-alkenals and oxidized phospholipids active in physiological and/or pathological conditions. Chem Phys Lipids. 2009;157(1):1–11.

    Article  CAS  Google Scholar 

  28. Wright MM, Howe AG, Vanina Z. Cell membranes and apoptosis: role of cardiolipin, phosphatidylcholine, and anticancer lipid analogues. Biochem Cell Biol. 2004;82(1):18–26.

    Article  CAS  PubMed  Google Scholar 

  29. Xie X, Tome ME, Gerner EW. Loss of intracellular Putrescine Pool-size regulation induces apoptosis. Exp Cell Res. 1997;230(2):386–92.

    Article  CAS  PubMed  Google Scholar 

  30. Tome EM, Fiser MS, Payne MC, Gerner WE. Excess putrescine accumulation inhibits the formation of modified eukaryotic initiation factor 5A (eIF-5A) and induces apoptosis. Biochem J. 1997;328(Pt 3):847.

    Article  CAS  Google Scholar 

  31. Laird PW. The power and the promise of DNA methylation markers. Nat Rev Cancer. 2003;3(4):253–66.

    Article  CAS  Google Scholar 

  32. Benavides MA, Oelschlager DK, Huang-Ge Z, Stockard CR, Vital-Reyes VS, Katkoori VR, Upender M, Wenquan W, Bland KI, Grizzle WE. Methionine inhibits cellular growth dependent on the p53 status of cells. Am J Surg 2007, 193(2):0–283. (DOI:

  33. Cheranov SY, Karpurapu M, Wang D, Zhang B, Venema RC, Rao GNL. An essential role for SRC-activated STAT-3 in 14, 15-EET-induced VEGF expression and angiogenesis. Blood. 2008;111(12):5581–91.

    Article  CAS  Google Scholar 

  34. Shao J, Li Q, Wang H, Liu F, Jiang J, Zhu X, Chen Z, Zou P. P-450-dependent epoxygenase pathway of arachidonic acid is involved in myeloma-induced angiogenesis of endothelial cells. J Huazhong Univ Sci Technol. 2011;31(5):596–601.

    Article  CAS  Google Scholar 

Download references


The authors would like to thank all of the colleagues and students who contributed to this study.


This work was supported by the National Natural Science Foundation of China (NO. 30973885; 81373889; 81973708). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations



GRM and MXD performed the experiments, analyzed the data and drafted the manuscript. SCJ and ZY assisted in carrying out experiments. SLS and JAD participated in the conception and design of the study. QDW and OYZ performed the statistical analysis and revised the manuscript. All authors read and approved the final version of the manuscript.

Corresponding authors

Correspondence to Shulan Su or Jinao Duan.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Figure S1-S15 Untreated GAPDH, p-JAK1, JAK1, p-STAT3 and STAT3 protein bands images (n=3).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gao, R., Miao, X., Sun, C. et al. Frankincense and myrrh and their bioactive compounds ameliorate the multiple myeloma through regulation of metabolome profiling and JAK/STAT signaling pathway based on U266 cells. BMC Complement Med Ther 20, 96 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: