Skip to main content
Download PDF

Matrine inhibits the growth of natural killer/T-cell lymphoma cells by modulating CaMKIIγ-c-Myc signaling pathway

BMC Complementary Medicine and Therapies volume 20, Article number: 214 (2020) Cite this article

Abstract

Background

C-Myc overexpression is associated with poor prognosis and aggressive progression of natural killer/T-cell lymphoma (NKTCL). Matrine, a main alkaloid of the traditional Chinese herb Sophora flavescens Ait, has been shown to inhibit cellular proliferation and induce apoptosis of various cancer cells. The present study investigated the effects and possible mechanisms of matrine inhibiting the growth of natural killer/T-cell lymphoma cells.

Methods

The effects of matrine on the proliferation, apoptosis and expression of apoptotic molecules, STAT3, LMP1, RUNX3, EZH2 and activation of CaMKIIγ/c-Myc pathway were examined in cultured NKTCL cell line NK92 cells.

Results

In cultured NK92 cells, matrine inhibited the proliferation in a dose and time dependent manner. The IC50 value of matrine was 1.71 mM for 72 h post exposure in NK92 cells. Matrine induced apoptosis with decreased Bcl-2 expression and the proteasome-dependent degradation of c-Myc protein in NK92 cells. c-Myc protein half-life in NK92 was reduced from 80.7 min to 33.4 min after matrine treatment, which meant the stability of c-Myc was decreased after matrine exposure. Furthermore, we found that matrine downregulated c-Myc phosphorylation at Ser62 together with the inhibition of CaMKIIγ, a key regulator of c-Myc protein in NKTCL. The downregulation of c-Myc transcription by matrine was mediated through LMP1 inhibition. We also observed that anti-proliferative activity of matrine was irrelevant to STAT3, RUNX3 and EZH2.

Conclusions

The results of the present study indicated that matrine inhibits the growth of natural killer/T-cell lymphoma cells by modulating LMP1-c-Myc and CaMKIIγ-c-Myc signaling pathway.

Peer Review reports

Background

NK/T-cell lymphoma (NKTCL), or extranodal NK/T-cell lymphoma, nasal type as classified by the World Health Organization, is an aggressive non-Hodgkin lymphoma originated from NK cells or cytotoxic T cells with a strong association with Epstein Barr Virus (EBV) [1]. It has a predilection for the nose and upper aerodigestive tract tissues. NKTCL is prevalent in Asian and Latin American, and usually affects males than females. The standard therapy for NKTCL is radiotherapy with or without L-asparaginase based combination chemotherapy. However, outcomes are poor with a five-year overall survival of 52% at best for relapsed or refractory NKTCL patients [2,3,4]. Effective targeted therapy is urgently needed, especially in the relapsed or refractory setting.

MYC is recognized as a master regulator of numerous biological and disease processes. c-Myc is extensively expressed in the normal cells, indicating high proliferative capacity [5]. The half-lives of c-Myc mRNA and protein are very short. The degradation of c-Myc protein is strictly controlled by rapid ubiquitin-mediated proteolysis [6]. Deregulation of the c-Myc is known to play a key role in the pathogenesis of numerous cancers including lymphoma [7]. c-Myc is overexpressed in the most NKTCL where translocations of MYC gene don’t exist [8]. The overexpression of c-Myc and the anti-apoptotic protein Bcl-2 has been correlated with poor prognosis [9]. c-Myc remains to be an “undruggable target” [10]. It is intrinsically disordered, lacking globular functional domains. Direct inhibition of c-Myc is a big challenge in cancer medicine.

Matrine is a pleiotropic alkaloid isolated from Sophora flavescens Ait, which has pharmacological activities including anti-inflammatory, anti-viral, and anti-fibrotic activities [11,12,13,14]. Recently, several studies have demonstrated that matrine has antitumor activity against various types of cancers including leukemia, multiple myeloma, gastric cancer [15,16,17,18]. However, the precise mechanism underlying the antitumor functions of matrine remains unclear. Therefore, we designed this study to investigate the antitumor effect of matrine in human NKTCL cells and its related molecular mechanism.

Methods

Cell lines and reagents

The human NKTCL NK92 cell line was obtained from the DSMZ collection (Germany) and maintained in MEM alpha medium supplemented with 12.5% fetal bovine serum (Gibco), 12.5% horse serum (Gibco), 10 ng/mL IL-2 (PeproTech, USA) in a humidified 5% CO2 atmosphere at 37 °C. Matrine, purchased from Nanjing Zelang Medical Technology Co., Ltd. (China), was dissolved in MEM alpha medium. Vindesine sulfate, purchased from Hangzhou Minsheng Pharmaceutical Co., Ltd. (China), was dissolved in 0.9% NaCl. Methylthiazolyldiphenyl-tetrazolium bromide (MTT) was obtained from Amresco (USA). MG132 and cycloheximide (CHX) were purchased from Cayman Chemical (USA). Dead Cell Apoptosis Kit (Cat# V13241) and TRIzol were bought from Invitrogen (USA). HiScript II Q RT reagent kit and ChamQ™ SYBR qPCR Master Mix kit were purchased from Vazyme (Nanjing, China). The antibodies for Caspase-3 (Cat# 9662), PARP (Cat# 9532), Bcl-2 (Cat# 4223), Bax (Cat# 5023), Stat3 (Cat# 12640), phospho-STAT3 (Tyr 705) (Cat# 9145) and phospho-c-Myc (Ser62) (Cat# 13748) were obtained from Cell Signaling Technology (USA). The antibodies for c-Myc (Cat# ab32072) and EBV LMP1 (Cat# ab78113) were obtained from Abcam (USA). Ca2+/calmodulin-dependent protein kinase II γ (CaMKIIγ) antibody (Cat# AP7208a) was obtained from Abgent (Suzhou, China). GAPDH antibody (Cat# 60004–1-Ig) was purchased from Proteintech (USA). Prestained and western blot marker was bought from Haigene (Harbin, China).

PBMC isolation

Whole blood samples were collected in 6 mL ethylene diamine tetraacetic acid-K2 vacuum blood tubes from healthy donors and mixed immediately after the collection by inverting 10 times. Peripheral blood mononuclear cells (PBMCs) isolation was performed immediately by using Ficoll density gradient centrifugation. Total 10 mL whole blood was used to mix with 10 mL 0.9% NaCl. Each 4 mL diluted blood sample was carefully layered on the 4 mL Ficoll medium (Tianjin Haoyang Biological Manufacture Co., Ltd., China) in 15 mL conical tube. The samples were then centrifuged for 20 min at 600×g, 20 °C. The lymphocyte-containing band in the opaque interface was carefully transferred into a clean 50 mL conical tube. The collected PBMCs washed three times in DMEM medium were maintained in DMEM medium with 10% FBS, 4 mg/L PHA (Sigma-Aldrich, USA) and used for cell viability assay on matrine immediately.

Cell viability analysis

The cell viability was determined by MTT assay. About 40,000 NK92 cells and 200,000 PBMCs were plated into each well of 96 well plates, respectively. NK92 cells were treated by increasing concentrations of matrine (125, 250, 500, 1000, 2000, 4000 μM) for 24 h, 48 h and 72 h, and vindesine sulfate (positive control) (0.032, 0.16, 0.8, 4, 20, 100 μM) for 72 h. PBMCs were treated by the same increasing concentrations of matrine for 72 h. Negative control cells were treated with MEM alpha medium and 0.9% NaCl, respectively. MTT (5 mg/mL) was added for the viability assay. The absorbance of the solution was read by a microplate reader ELx808 (BioTek, USA), using a test wavelength of 490 nm. Results obtained were expressed as percentage inhibition rate to test agents. Half maximal inhibitory concentration (IC50) was calculated by GraphPad Prism program.

Flow cytometric analysis for apoptosis

Approximately 5 × 105 NK92 cells were seeded into each well of six well plates and then incubated with 1, 2 and 4 mM matrine for 48 h. NK92 cells mixed with MEM alpha medium only were used for the control. The cells were washed with ice-cold PBS and then resuspended in 1 × binding buffer. Each sample was incubated with 5 μL Alexa Fluor 488 annexin V and 1 μL 100 μg/mL PI for 15 min protected from light at room temperature according to the manufacturer’s instruction. The samples were analyzed using flow cytometry (Navios, Beckman Coulter, USA) within an hour.

Quantitative RT-PCR analysis

Approximately 7.5 × 105 NK92 cells were seeded into each well of six well plates and then incubated with 1.2, 2.4 and 3.6 mM matrine for 24 h. NK92 cells mixed with MEM alpha medium only were used for the control. The cells were washed with ice-cold PBS and total RNA, including small RNAs, was extracted by TRIzol. cDNAs were synthesized from total RNA using the HiScript II Q RT reagent kit with gene-specific primers. The ChamQ™ SYBR qPCR Master Mix kit was used for the thermocycling reaction according to the manufacturer’s instructions in the CFX384 Real Time PCR system (Bio-Rad, USA). The threshold cycle (Ct) was determined using default threshold settings. All experiments were done in triplicates. β-Actin and U6 snRNA were used as the controls to normalize mRNA and miRNA input, respectively. The relative gene expression quantification was calculated with the 2−ΔΔct method.

Western blotting

Cells were washed twice with PBS buffer, and total cellular proteins were extracted and subjected to SDS-PAGE, and then transferred to PVDF membranes (Merck, USA) and blocked with 5% nonfat milk in TBS–Tween 20 (TBST). The membranes were then reacted with primary antibodies overnight at 4 °C. After 3 washes with TBST, membranes were probed with a horseradish peroxidase-conjugated secondary antibody for 1 h at room temperature, and reacted with Immobilon Western Chemiluminescent HRP Substrate (Millipore, USA). Protein levels were measured with the densitometric intensity.

Cycloheximide chase analysis

NK92 cells (7.5 × 105) were treated with or without 1.96 mM matrine for 12 h. Cells were then treated with 100 μg/mL cycloheximide and harvested at indicated time points, and western blotting was performed.

Statistical analysis

Data were expressed as mean ± standard deviation. Student’s t-test was applied for comparison of the means of two groups, and one way Analysis of Variance (ANOVA) was used to assess the level of significance between the means of multiple groups. Statistical significance was defined as p < 0.05.

Results

Matrine inhibits the growth of NK92

The cell viability assay was performed to evaluate the percentage inhibition rate and IC50 of matrine in NK92 cells and PBMCs. As shown in Fig. 1a, matrine displayed anti-proliferative activity in NK92 cells in a dose dependent and time dependent manner. The IC50 values of matrine in NK92 cells were 6.32 ± 0.02, 1.96 ± 0.03 and 1.71 ± 0.05 mM, respectively, for 24 h, 48 h and 72 h while that of vindesine in NK92 cells for 72 h was 0.64 ± 0.03 μM. The effect of matrine on the normal lymphocytes was explored by using PBMCs isolated from five healthy donors for MTT assay. The percentage inhibition rate of matrine on the proliferation of normal PBMCs induced by PHA for 72 h was 0.88 ± 0.51, 2.32 ± 0.42 and 5.35 ± 0.41%, respectively, for 0.5 mM, 1 mM and 2 mM matrine treatment (Fig. 1a). Although exposure to matrine for a longer time than 48 h was found to be more potent in inhibiting NK92 cell viability, it inclined to result in increased frequency of necrotic cells in NK92 cells, so we determined the treatment with matrine for no longer than 48 h for the following experiments in NK92 cells.

Fig. 1
figure 1

Anti-proliferation and apoptosis induction of matrine in NKTCL cells. a NK92 cells and PBMCs were treated with matrine and vindesine at different concentrations for different times. The total viable cells were determined by MTT assay. b NK92 cells were exposed to matrine at different concentrations for 48 h and then determined for apoptotic cells by annexin V and PI staining using flow cytometry. c Percentage (%) of apoptotic cells induced by matrine at various concentrations. Analyses in triplicates. (*p < 0.05, **p < 0.01 compare to 0 mM group)

Full size image

Matrine induces apoptosis

To confirm whether the growth inhibition of NK92 cells induced by matrine was caused by apoptosis, NK92 cells were exposed to matrine at 0, 1, 2 and 4 mM for 48 h, and the occurrence of apoptosis was identified using the annexin V and PI staining because annexin V bound membrane phosphatidyl serine in the apoptosis and viability dye PI resolved late-stage apoptotic and necrotic cells from early-stage apoptotic cells (Fig. 1b). Overall, the result in Fig. 1c shows a significant increase in the percentage of apoptotic cells from 3.83 ± 0.78% in untreated cells to 7.9 ± 0.66, 13.33 ± 0.7 and 62.9 ± 0.66%, respectively, for 1 mM, 2 mM and 4 mM matrine treatment.

Matrine regulates the expression of apoptosis-related proteins

To explore the mechanism responsible for matrine mediated apoptosis, the apoptotic protein expressions were evaluated by western blot. Figure 2a and Supplementary Figure 1 show the results of western blot for cleaved PARP, cleaved Capase-3, Bcl-2, and Bax proteins relative to GAPDH in control cultures and cultures exposed to 1.96 mM (IC50 of 48 h) matrine for 48 h. Matrine treatment induced the significant upregulation of cleaved PARP and cleaved Capase-3 with normalized relative expression of 1.46 ± 0.27 and 1.63 ± 0.30, respectively (Fig. 2b). Significantly decreased expression of Bcl-2 was induced with normalized relative expression of 0.83 ± 0.06 while the expression of Bax was not significantly influenced by matrine. Therefore, the ratio of Bcl-2/Bax proteins was significantly decreased to 0.77 ± 0.10 after matrine treatment. The results clearly indicated that matrine induces apoptosis of NK92 cells via activation of the mitochondrial pathway.

Fig. 2
figure 2

Matrine induced apoptosis of NKTCL cells via activation of the mitochondrial pathway. a NK92 cells were treated with 1.96 mM matrine for 48 h, followed by western blot. GAPDH was used as a loading control. b The relative intensities of target proteins were normalized to those of GAPDH. Analyses in triplicates. (n.s., not significant; *p < 0.05)

Full size image

Matrine inhibits the growth of NK92 independent of JAK/STAT3 pathway

Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway was deregulated in most NKTCL cases evaluated by STAT3 phosphorylation [19]. STAT3 is activated by phosphorylation at Tyr705, which induces dimerization, nuclear translocation and DNA binding. The phosphorylation of STAT3 at Tyr705 in NK92 cells was then analyzed. NK92 cells were treated with matrine at 1.96 mM for 48 h. The cells were then collected for analysis of phospho-STAT3 (Tyr705) and STAT3 expression by western blot. As shown in Fig. 3 and Supplementary Figure 2, the level of phospho-STAT3 (Tyr705) protein in NK92 cells was significantly increased after matrine treatment. The results demonstrated that the anti-proliferative activity of matrine in NK92 is independent of JAK/STAT3 pathway. The growth inhibition of NK92 on matrine may be mediated by the suppression of the vital targets of STAT3 transcriptional activity.

Fig. 3
figure 3

Matrine inhibited NKTCL cells independent of JAK/STAT3 pathway. a NK92 cells were treated with 1.96 mM matrine for 48 h, followed by western blot for STAT3, p-STAT3 (Tyr705) antibodies. GAPDH was used as a loading control. b The relative intensities of target proteins were normalized to those of GAPDH. Analyses in triplicates. (*p < 0.05, **p < 0.01)

Full size image

Matrine downregulates c-Myc expression by decreased transcription and increased degradation

c-Myc is the principal target of STAT3 transcriptional activity and is overexpressed in the most NKTCL [20] . The expression of c-Myc protein in NK92 cells was then analyzed by western blot. NK92 cells were treated with matrine at 1.96 mM for 48 h. The levels of c-Myc protein in NK92 cells were reduced with normalized relative expression of 0.70 ± 0.03 after matrine treatment (Fig. 4a and Supplementary Figure 3). These results indicated that the growth inhibition induced by matrine in NK92 cells is associated with the downregulation of c-Myc protein.

Fig. 4
figure 4

Decreased c-Myc protein induced by matrine and rescued by proteasome inhibitor. a Effect of matrine on c-Myc protein expression in NK92 cells. NK92 cells were treated with matrine at 1.96 mM for 48 h, and c-Myc protein levels were measured by western blot. b c-Myc mRNA levels in NK92 cells were determined by quantitative RT-PCR at 24 h after matrine treatment at 0, 1.2, 2.4, and 3.6 mM (*p < 0.05, **p < 0.01 compared to 0 mM group). c CHX chase assay for the half-time of c-Myc. NK92 cells were treated with or without 1.96 mM matrine for 12 h. Cells were then treated with CHX (100 μg/mL) for the indicated minutes, and western blotting was performed. d c-Myc levels were quantified relative to GAPDH levels and graphed as percent c-Myc protein remaining after CHX treatment. Half-lives of c-Myc were calculated from exponential line equations and shown for each treat. e c-Myc protein levels were determined at 6 h post-treatment of MG132 and/or matrine, and (f) the relative intensities of target proteins were normalized to those of GAPDH. Analyses in triplicates. (*p < 0.05, **p < 0.01)

Full size image

Transcription analysis of c-Myc gene was then processed. NK92 cells were incubated with 1.2, 2.4 and 3.6 mM matrine, respectively, for 24 h. The qRT-PCR results showed that c-Myc mRNA levels were significantly reduced after matrine treatment in a dose dependent manner (Fig. 4b). Normalized mRNA expression of c-Myc was 0.84 ± 0.04, 0.65 ± 0.02 and 0.57 ± 0.05, respectively, for 1.2 mM, 2.4 mM and 3.6 mM matrine treatment (Fig. 4b).

The degradation of c-Myc protein was also analyzed by cycloheximide chase assay. The translation inhibitor cycloheximide was added into the control NK92 cells and matrine-treated NK92 cells, and c-Myc protein levels were evaluated at different time points by western blot. Half-lives of c-Myc protein were then calculated by GraphPad Prism program. As shown in Fig. 4c, d and Supplementary Figure 4, c-Myc half-life in matrine-treated NK92 cells was about 33.4 min, while that in the control NK92 cells was about 80.7 min. The results indicated that c-Myc degradation in matrine-treated NK92 is promoted, and the stability of c-Myc in matrine-treated NK92 is decreased.

To explore the possible degradation pathway of c-Myc in NK92, the proteasome inhibitor MG132 was used to block the activity of proteasome. c-Myc protein levels were determined at 6 h post-treatment of MG132 with or without matrine by western blot. As shown in Fig. 4e, f and Supplementary Figure 5, MG132 treatment prevented matrine-induced c-Myc protein degradation. The results demonstrated that matrine promotes c-Myc protein degradation in NK92 cells in a proteasome-dependent manner.

Matrine downregulates c-Myc phosphorylation at Ser62 through CaMKIIγ inhibition

c-Myc protein stability is regulated by two phosphorylation sites with opposite functions: serine 62 phosphorylation stabilizes c-Myc while threonine 58 phosphorylation promotes c-Myc degradation [21]. The phosphorylation of c-Myc at Ser62 in NK92 cells was then analyzed by western blot. NK92 cells were treated with matrine for 48 h. The levels of phospho-c-Myc (Ser62) in NK92 cells were significantly decreased after matrine treatment suggesting the stability of c-Myc is declined (Fig. 5a, b and Supplementary Figure 6).

Fig. 5
figure 5

Matrine inhibited NKTCL cells through CaMKIIγ/c-Myc pathway. a NK92 cells were treated with 1.96 mM matrine for 48 h, followed by western blot for c-Myc, p-c-Myc (Ser62), CaMKIIγ and LMP1 antibodies. GAPDH was used as loading control. b The relative intensities of target proteins were normalized to those of loading control (**p < 0.01). c LMP1, RUNX3, EZH2, miR-26a, miR-26b and miR-101 transcription levels in NK92 cells were determined by quantitative RT-PCR at 24 h after matrine treatment at 0, 1.2, 2.4, and 3.6 mM (*p < 0.05, **p < 0.01 compared to 0 mM group). Analyses in triplicates

Full size image

CAMKIIγ was recently identified to be one of 102 potential genes involved in a synthetic lethal interaction with c-Myc [22]. CaMKIIγ phosphorylates c-Myc at Ser62 directly in T cell lymphoma [23]. The regulation of CaMKIIγ in matrine treatment NK92 cells was then analyzed by western blot. As shown in Fig. 5a, b and Supplementary Figure 6, the levels of CaMKIIγ protein were also significantly decreased after matrine treatment. These results indicated that matrine downregulates c-Myc phosphorylation at Ser62 by targeting CaMKIIγ.

Matrine downregulates c-Myc expression through LMP1 inhibition

c-Myc is a transcriptional target of the EBV protein LMP1 [24], it has been proposed that c-Myc upregulation in NKTCL is mediated by EBV [8]. The regulation of LMP1 in matrine treatment NK92 cells was analyzed by western blot. NK92 cells were treated with matrine for 48 h. As shown in Fig. 5a, b and Supplementary Figure 6, significantly decreased expression of LMP1 was induced with normalized relative expression of 0.59 ± 0.13 after matrine treatment. Transcription analysis of LMP1 gene was also processed in NK92 cells. The qRT-PCR results showed that LMP1 mRNA levels were significantly increased after matrine treatment (Fig. 5c). The results showed that the downregulation of c-Myc transcription by matrine is mediated through LMP1 inhibition. The decreased protein expression of LMP1 was caused by the translation inhibition, not due to the transcription inhibition.

Matrine downregulates c-Myc expression irrelevant to RUNX3 and EZH2

RUNX3 is a known target of c-Myc transcriptional activity. It is oncogenic in NKTCL [25]. Transcription analysis of RUNX3 gene in NK92 cells showed that RUNX3 mRNA levels were significantly increased after matrine treatment (Fig. 5c), which implied that RUNX3 transcription is not decreased by downregulated c-Myc protein after matrine treatment.

c-Myc activation has been shown to lead to EZH2 overexpression by suppression of its negative regulatory miRNAs, such as miR-26 and miR-101, in NKTCL [26]. EZH2 silences the tumor suppressor genes with its histone methyltransferase. Transcription analysis of EZH2, miR-26 and miR-101 genes in NK92 cells showed that EZH2 mRNA levels were increased after matrine treatment, while the miRNA levels of miR-26a, miR-26b and miR-101 were decreased (Fig. 5c). These results indicated that miRNA transcriptions of miR-26a, miR-26b and miR-101 are not negatively regulated by decreased c-Myc protein after matrine treatment, and corresponding EZH2 transcription is not inhibited. The above mentioned results demonstrated that downregulated c-Myc expression by matrine is irrelevant to RUNX3 and EZH2 in NK92 cells.

Discussion

Matrine is a natural tetracyclo-quinolizindine alkaloid with potent antitumor functions towards various types of cancers. Ma et al. reported that matrine inhibited the growth of human chronic myeloid leukemia K562 cells with an IC50 of 2 mM for 48 h treatment [15]. Han et al. showed that matrine inhibited human myeloma cell lines RPMI8226 and U266 with IC50 at 48 h of 4.55 mM and 5.36 mM, respectively [17]. In present study, matrine inhibited the growth of NKTCL cell line NK92 cells in a dose and time dependent manner. It had anti-proliferative activity against NK92 cells with IC50 of 1.96 ± 0.03 mM and 1.71 ± 0.05 mM, respectively, for 48 h and 72 h treatment. Vindesine is an alkaloid for the treatment of non-Hodgkin’s lymphoma in clinic. It inhibited the growth of NK92 cells with an IC50 of 0.64 ± 0.03 μM for 72 h treatment. Since the IC50 of matrine is much higher than that of vindesine, the effect of matrine on the normal lymphocytes is of much concern. No significant inhibition on the proliferation of PBMCs induced by PHA was found with 0.5 mM matrine treatment. Only slight inhibition of PBMCs was seen between 1 mM and 2 mM matrine treatment for 72 h (Fig. 1a). Han et al. reported that lower concentrations of matrine (1, 2, 4, 6 mM) had no effects on the proliferation of PBMCs and higher concentrations of matrine (8, 12, 20 mM) suppressed the proliferation of PBMCs within 72 h [17]. They also showed that 2, 6, 12 and 20 mM concentrations of matrine had no effects on the induction of apoptosis of PBMCs for 48 h [17]. Our data and previous literature support that matrine with lower than 2 mM has no effects on the normal PBMCs, which will be beneficial for NKTCL patients.

Our previous study reported that matrine induced apoptosis in leukemia cells [16]. In this study, matrine induced apoptosis of NKTCL cell line NK92 cells in a dose dependent manner from 3.83 ± 0.78% to 62.9 ± 0.66% after matrine treatment at different concentrations for 48 h (Fig. 1c). It is interesting to notice that apoptosis induced by 2 mM matrine was too small (less than 14%) even though this dose was more than IC50 (1.96 mM), which suggests matrine does not preferentially induce apoptosis. The growth inhibition of NK92 cells induced by matrine was partially caused by apoptosis. Furthermore, matrine upregulated the products of cleaved Capase-3 and cleaved PARP, and downregulated the expression of Bcl-2 and reduced the ratio of Bcl-2/Bax (Fig. 2b). Our finding indicated that matrine induces apoptosis of NK92 cells through the activation of the mitochondrial pathway.

JAK/STAT pathway plays an important role in the pathogenesis of NKTCL through its pro-proliferative activity [20]. STAT3 phosphorylation activation at Tyr705 is shown to provide a growth advantage by upregulating the expression of c-Myc through binding to its promoter in NKTCL [19, 20]. In present study, matrine promoted the STAT3 phosphorylation at Tyr705 in NK92 cells while the growth of cells was inhibited (Fig. 3). We speculate that the expression of vital targets genes of pSTAT3, such as c-Myc, are inhibited by matrine. This speculation was confirmed by the expression analysis of c-Myc (Fig. 4). Our results show that JAK/STAT3 pathway is not involved in anti-proliferative activity of matrine in NK92.

c-Myc is the principal target of STAT3 transcriptional activity and is overexpressed extensively in NKTCL. Ng SB et al. reported that 45.4% NKTCLs had overexpressed c-Myc by immunohistochemistry [8]. Genome-wide gene expression profiling studies have identified the activation of c-Myc in NKTCL. c-Myc is known to cause widespread miRNA repression [27], it is suggested that c-Myc activation may be one of the mechanism in the downregulation of miRNA in NKTCL [28]. Our data showed that c-Myc protein expression was inhibited by matrine in NK92 cells. c-Myc gene transcription was decreased, and the degradation of c-Myc protein was accelerated by matrine treatment. c-Myc protein half-life was much shorter after exposure to matrine in NK92 cells, which meant the stability of c-Myc in matrine-treated NK92 is declined (Fig. 4). It is interesting to find that matrine promoted c-Myc protein degradation in NK92 cells through a proteasome-dependent manner (Fig. 4). Our data showed that matrine is an effective anti-proliferative agent that is primarily mediated through inhibition of c-Myc.

Recently, Gu et al. reported that CaMKIIγ phosphorylated Ser62 of c-Myc and facilitated the stability of c-Myc in T cell lymphoma. Inhibition of CaMKIIγ alleviated T cell lymphoma burden in mice [23]. As we know, most of T cell lymphoma has c-Myc protein overexpression. Since c-Myc protein overexpression also exists in NKTCL, we postulate that CaMKIIγ may also play an important role in NKTCL. As expected, a positive correlation between CaMKIIγ and p-c-Myc (Ser62)/c-Myc was observed in NK92 cells. The levels of phospho-c-Myc (Ser62) and CaMKIIγ in NK92 cells were together remarkably decreased after matrine treatment (Fig. 5). Our findings demonstrated that matrine inhibits the growth of NKTCL cells by modulating CaMKIIγ-c-Myc pathway. Although the precise role of CaMKIIγ in NKTCL is not known, we propose that CaMKIIγ affects the stability of c-Myc protein and may be responsible for sustaining malignant growth of NKTCL. Our data supported that CaMKIIγ inhibition is an effective strategy for targeting c-Myc-driven NKTCL.

Previous study reported that oncoprotein LMP1 of EBV directly regulated the transcription of c-Myc, and LMP1 was critical for EBV-mediated B-lymphocyte transformation [24]. It has been postulated that c-Myc activation in NKTCL is mediated primarily through LMP1 [8]. As expected, a positive correlation between LMP1 and c-Myc was observed in NK92 cells. LMP1 protein and c-Myc mRNA in NK92 cells were together remarkably decreased after matrine treatment (Fig. 5). The transcription inhibition of c-Myc was mediated by the downregulated LMP1 protein. Furthermore, we found that the decreased protein expression of LMP1 was caused by the translation inhibition, not due to the transcription inhibition, which implies that LMP1-specific miRNAs may involve this regulation with matrine treatment. Present findings indicated that matrine inhibits the growth of NKTCL cells by modulating LMP1-c-Myc pathway.

Direct inhibition of c-Myc remains an elusive goal in cancer medicine. Targeting c-Myc function via BET bromodomain inhibition has been validated by studies in blood cancer. JQ1, the first small molecule bromodomain inhibitor, caused the down-regulation of the c-Myc transcription and c-Myc-dependent target genes in acute myeloid leukemia, Burkitt lymphoma and diffuse large B-cell lymphoma [29, 30]. Alkaloids represent an important group of anticancer drugs. Several alkaloids, such as vindesine, berberine and matrine, are well-known as potent chemotherapeutic agents. Quinolino-benzo-[5,6]-dihydroisoquindolium compounds of berberine derivatives, an effective G-quadruplex ligand targeting c-Myc oncogene, resulted in the down-regulation of c-Myc gene transcription in leukemia cell line HL60 [31]. Recently, Yu et al. reported that tosyl chloride-berbamine, small molecule analog of berbamine, eliminated c-Myc-positive leukemia in vitro and in vivo by targeting CaMKIIγ/Myc axis [32]. Our current study showed that matrine inhibited the transcription of c-Myc by targeting LMP1 and promoted c-Myc protein degradation by targeting CaMKIIγ in the NK92 cells (Fig. 6). CaMKIIγ/Myc axis could be a valid target for developing small molecule-based new therapies for treating MYC-mediated NKTCL. Matrine could be beneficial for c-Myc-driven NKTCL patients.

Fig. 6
figure 6

Cartoon diagram of the mechanisms for growth suppression of NKTCL cells by matrine. Matrine inhibits the transcription of c-Myc through the downregulated LMP1 protein. Matrine downregulates c-Myc phosphorylation at Ser62 through CaMKIIγ inhibition, and then promotes the destabilization and degradation of c-Myc protein in a proteasome-dependent manner

Full size image

Limitations of this study include the focus on the specific NKTCL cell line. The precise mechanism of c-Myc downregulation by matrine in NKTCL requires further investigation. The characterization of these responses to primary NKTCL cells remains to be determined. Therefore, the effects of matrine on NKTCL in vivo will be determined in the future studies.

Conclusions

In summary, this study demonstrated that matrine possesses anti-proliferation and apoptosis induction effects in NK/T-cell lymphoma cells. Moreover, to the best of our knowledge, the present study was the first to show the mechanisms for suppression of NKTCL by matrine may involve the inhibition of LMP1-c-Myc and CaMKIIγ-c-Myc signaling pathway. There is no available therapeutic drugs targeting c-Myc in NKTCL, our data support that matrine could be a promising compound for treatment of c-Myc-driven NK/T-cell lymphoma.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CaMKIIγ:

Ca2+/calmodulin-dependent protein kinase II γ

CHX:

Cycloheximide

EBV:

Epstein barr virus

IC50 :

Half maximal inhibitory concentration

MTT:

Methylthiazolyldiphenyl-tetrazolium bromide

NKTCL:

Natural killer/T-cell lymphoma

STAT:

Signal transducer and activator of transcription

References

  1. Swerdlow SH, Campo E, Pileri SA, Harris NL, Stein H, Siebert R, Advani R, Ghielmini M, Salles GA, Zelenetz AD, Jaffe ES. The 2016 revision of the World Health Organization classification of lymphoid neoplasms. Blood. 2016;127:2375–90.

    Article  CAS  Google Scholar 

  2. Qi SN, Xu LM, Yuan ZY, Wu T, Zhu SY, Shi M, Su H, Wang Y, He X, Zhang LL, Wu G, Qu BL, Qian LT, Hou XR, Zhang FQ, Zhang YJ. Effect of primary tumor invasion on treatment and survival in extranodal nasal-type NK/T-cell lymphoma in the modern chemotherapy era: a multicenter study from the China lymphoma collaborative group (CLCG). Leuk Lymphoma. 2019. https://doi.org/10.1080/10428194.2019.1602265.

  3. Tse E, Kwong YL. The diagnosis and management of NK/T-cell lymphomas. J Hematol Oncol. 2017;10:85.

    Article  Google Scholar 

  4. Kwong YL, Kim WS, Lim ST, Kim SJ, Tang T, Tse E, Leung AY, Chim CS. SMILE for natural killer/T-cell lymphoma: analysis of safety and efficacy from the Asia lymphoma study group. Blood. 2012;120:2973–80.

    Article  CAS  Google Scholar 

  5. Pelengaris S, Khan M, Evan G. C-MYC: more than just a matter of life and death. Nat Rev Cancer. 2002;2:764–76.

    Article  CAS  Google Scholar 

  6. Salghetti SE, Kim SY, Tansey WP. Destruction of Myc by ubiquitin-mediated proteolysis: cancer-associated and transforming mutations stabilize Myc. EMBO J. 1999;18:717–26.

    Article  CAS  Google Scholar 

  7. Gabay M, Li Y, Felsher DW. MYC activation is a hallmark of cancer initiation and maintenance. Cold Spring Harb Perspect Med. 2014. https://doi.org/10.1101/cshperspect.a014241.

  8. Ng SB, Selvarajan V, Huang G, Zhou J, Feldman AL, Law M, Kwong YL, Shimizu N, Kagami Y, Aozasa K, Salto-Tellez M, Chng WJ. Activated oncogenic pathways and therapeutic targets in extranodal nasal-type NK/T cell lymphoma revealed by gene expression profiling. J Pathol. 2011;223:496–510.

    Article  CAS  Google Scholar 

  9. Wang JH, Bi XW, Li PF, Xia ZJ, Huang HQ, Jiang WQ, Zhang YJ, Wang L. Overexpression of MYC and BCL2 predicts poor prognosis in patients with Extranodal NK/T-cell lymphoma, nasal type. J Cancer. 2017;8:793–800.

    Article  CAS  Google Scholar 

  10. McKeown MR, Bradner JE. Therapeutic strategies to inhibit MYC. Cold Spring Harb Perspect Med. 2014. https://doi.org/10.1101/cshperspect.a014266.

  11. Fu S, Sun C, Tao X, Ren Y. Anti-inflammatory effects of active constituents extracted from Chinese medicinal herbs against Propionibacterium acnes. Nat Prod Res. 2012;26:1746–9.

    Article  CAS  Google Scholar 

  12. Yang Y, Xiu J, Zhang X, Zhang L, Yan K, Qin C, Liu J. Antiviral effect of matrine against human enterovirus 71. Molecules. 2012;17:10370–6.

    Article  CAS  Google Scholar 

  13. Wan XY, Luo M, Li XD, He P. Hepatoprotective and anti-hepatocarcinogenic effects of glycyrrhizin and matrine. Chem Biol Interact. 2009;181:15–9.

    Article  CAS  Google Scholar 

  14. Zheng Y, Jiang X, Gao F, Song J, Sun J, Wang L, Sun X, Lu Z, Zhang H. Identification of plant-derived natural products as potential inhibitors of the Mycobacterium tuberculosis proteasome. BMC Complement Altern Med. 2014;14:400.

    Article  Google Scholar 

  15. Ma L, Zhu Z, Jiang L, Sun X, Lu X, Zhou M, Qian S, Jianyong L. Matrine suppresses cell growth of human chronic myeloid leukemia cells via its inhibition of the interleukin-6/Janus activated kinase/signal transducer and activator of transcription 3 signaling cohort. Leuk Lymphoma. 2015;56:2923–30.

    Article  CAS  Google Scholar 

  16. Wu D, Shao K, Zhou Q, Sun J, Wang Z, Yan F, Liu T, Wu X, Ye B, Huang H, Zhou Y. Autophagy and ubiquitin-mediated Proteolytic degradation of PML/Rarα fusion protein in Matrine-induced differentiation sensitivity recovery of ATRA-resistant APL (NB4-LR1) cells: in vitro and in vivo studies. Cell Physiol Biochem. 2018;48:2286–301.

    Article  CAS  Google Scholar 

  17. Han Y, Zhang S, Wu J, Yu K, Zhang Y, Yin L, Bi L. Matrine induces apoptosis of human multiple myeloma cells via activation of the mitochondrial pathway. Leuk Lymphoma. 2010;51:1337–46.

    Article  CAS  Google Scholar 

  18. Li Y, Zhang J, Ma H, Chen X, Liu T, Jiao Z, He W, Wang F, Liu X, Zeng X. Protective role of autophagy in matrine-induced gastric cancer cell death. Int J Oncol. 2013;42:1417–26.

    Article  CAS  Google Scholar 

  19. Küçük C, Jiang B, Hu X, Zhang W, Chan JK, Xiao W, Lack N, Alkan C, Williams JC, Avery KN, Kavak P, Scuto A, Sen E, Gaulard P, Staudt L, Iqbal J, Zhang W, Cornish A, Gong Q, Yang Q, Sun H, d'Amore F, Leppä S, Liu W, Fu K, de Leval L, McKeithan T, Chan WC. Activating mutations of STAT5B and STAT3 in lymphomas derived from γδ-T or NK cells. Nat Commun. 2015;6:6025.

    Article  Google Scholar 

  20. Coppo P, Gouilleux-Gruart V, Huang Y, Bouhlal H, Bouamar H, Bouchet S, Perrot C, Vieillard V, Dartigues P, Gaulard P, Agbalika F, Douay L, Lassoued K, Gorin NC. STAT3 transcription factor is constitutively activated and is oncogenic in nasal-type NK/T-cell lymphoma. Leukemia. 2009;23:1667–78.

    Article  CAS  Google Scholar 

  21. Sears R, Nuckolls F, Haura E, Taya Y, Tamai K, Nevins JR. Multiple Ras-dependent phosphorylation pathways regulate Myc protein stability. Genes Dev. 2000;14:2501–14.

    Article  CAS  Google Scholar 

  22. Toyoshima M, Howie HL, Imakura M, Walsh RM, Annis JE, Chang AN, Frazier J, Chau BN, Loboda A, Linsley PS, Cleary MA, Park JR, Grandori C. Functional genomics identififies therapeutic targets for MYC-driven cancer. Proc Natl Acad Sci U S A. 2012;109:9545–50.

    Article  CAS  Google Scholar 

  23. Gu Y, Zhang J, Ma X, Kim BW, Wang H, Li J, Pan Y, Xu Y, Ding L, Yang L, Guo C, Wu X, Wu J, Wu K, Gan X, Li G, Li L, Forman SJ, Chan WC, Xu R, Huang W. Stabilization of the c-Myc protein by CAMKIIgamma promotes T cell lymphoma. Cancer Cell. 2017;32:115–28.

    Article  CAS  Google Scholar 

  24. Dirmeier U, Hoffmann R, Kilger E, Schultheiss U, Briseño C, Gires O, Kieser A, Eick D, Sugden B, Hammerschmidt W. Latent membrane protein 1 of Epstein–Barr virus coordinately regulates proliferation with control of apoptosis. Oncogene. 2005;24:1711–7.

    Article  CAS  Google Scholar 

  25. Selvarajan V, Osato M, Nah GSS, Yan J, Chung TH, Voon DC, Ito Y, Ham MF, Salto-Tellez M, Shimizu N, Choo SN, Fan S, Chng WJ, Ng SB. RUNX3 is oncogenic in natural killer/T-cell lymphoma and is transcriptionally regulated by MYC. Leukemia. 2017;31:2219–27.

    Article  CAS  Google Scholar 

  26. Yan J, Ng SB, Tay JL, Lin B, Koh TL, Tan J, Selvarajan V, Liu SC, Bi C, Wang S, Choo SN, Shimizu N, Huang G, Yu Q, Chng WJ. EZH2 overexpression in natural killer/T-cell lymphoma confers growth advantage independently of histone methyltransferase activity. Blood. 2013;121:4512–20.

    Article  CAS  Google Scholar 

  27. Chang TC, Yu D, Lee YS, Wentzel EA, Arking DE, West KM, Dang CV, Thomas-Tikhonenko A, Mendell JT. Widespread microRNA repression by Myc contributes to tumorigenesis. Nat Genet. 2008;40:43–50.

    Article  CAS  Google Scholar 

  28. Ng SB, Yan J, Huang G, Selvarajan V, Tay JL, Lin B, Bi C, Tan J, Kwong YL, Shimizu N, Aozasa K, Chng WJ. Dysregulated microRNAs affect pathways and targets of biologic relevance in nasal-type natural killer/T-cell lymphoma. Blood. 2011;118:4919–29.

    Article  CAS  Google Scholar 

  29. Mertz JA, Conery AR, Bryant BM, Sandy P, Balasubramanian S, Mele DA, Bergeron L, Sims RJ III. Targeting MYC dependence in cancer by inhibiting BET bromodomains. Proc Natl Acad Sci. 2011;108:16669–74.

    Article  CAS  Google Scholar 

  30. Chapuy B, McKeown MR, Lin CY, Monti S, Roemer MG, Qi J, Rahl PB, Sun HH, Yeda KT, Doench JG, Reichert E, Kung AL, Rodig SJ, Young RA, Shipp MA, Bradner JE. Discovery and characterization of super-enhancer-associated dependencies in diffuse large B cell lymphoma. Cancer Cell. 2013;24:777–90.

    Article  CAS  Google Scholar 

  31. Ma Y, Ou TM, Tan JH, Hou JQ, Huang SL, Gu LQ, Huang ZS. Quinolino-benzo-[5, 6]-dihydroisoquindolium compounds derived from berberine: a new class of highly selective ligands for G-quadruplex DNA in c-myc oncogene. Eur J Med Chem. 2011;46:1906–13.

    Article  CAS  Google Scholar 

  32. Yu Q, Wang P, Yang L, Wu Z, Li S, Xu Y, Wu B, Ma A, Gan X, Xu R. Novel synthetic tosyl chloride-berbamine regresses lethal MYC-positive leukemia by targeting CaMKIIγ/Myc axis. Biomed Pharmacother. 2019;117:109134.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by Zhejiang Provincial Natural Science Foundation of China under Grants No.LY15H290004 and No.LY16H290002; National TCM Clinical Research Bases of China under Grant No.JDZX2015116; and Health Commission of Zhejiang Province under Grant No.2015ZDA021. The funding body provided financial support and had no role in the design of the study, collection, analysis, interpretation of data, or writing the manuscript.

Author information

Authors and Affiliations

  1. The First Affiliated Hospital, Zhejiang Chinese Medical University, No. 54 Youdian Road, Zhejiang, 310006, Hangzhou, China

    Jianyou Gu, Yu Zhang, Xiao Wang, Jingjing Xiang, Shu Deng, Dijiong Wu, Junfa Chen, Lihong Yu, Yan Zhou & Jianping Shen

  2. Key Laboratory of Integrative Chinese and Western Medicine for the Diagnosis and Treatment of Circulatory Diseases of Zhejiang Province, No. 54 Youdian Road, Zhejiang, 310006, Hangzhou, China

    Jianyou Gu & Xiao Wang

  3. Shaoxing Second Hospital, No. 123 Yanan Road, Shaoxing, Zhejiang, 312000, China

    Yaokun Wang

Authors
  1. Jianyou Gu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yu Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Xiao Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jingjing Xiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Shu Deng
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Dijiong Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Junfa Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Lihong Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Yan Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Yaokun Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Jianping Shen
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

JG and JS: conception and design of research; XW, JX, LY, YZ and YW performed experiments; YZ, SD, DW and JC analyzed data and interpreted results of experiments; JG and JS drafted manuscript; JG and JS edited and revised manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jianping Shen.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of the First Affiliated Hospital of Zhejiang Chinese Medical University (approval no. 2016-KL-019-02) and conducted according to the principles expressed in the Declaration of Helsinki. All participants gave written informed consent before the study.

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: Supplementary Figure 1

. Matrine induced the expression of apoptosis-related proteins in NKTCL cells. NK92 cells (5 × 105) were treated with 1.96 mM matrine for 48 h, followed by western blot. NK92 cells treated without matrine were used as control. (A) Representative WB result of Caspase-3 and cleaved Caspase-3. (B) Representative WB result of PARP and cleaved PARP. (C) Representative WB result of Bcl-2. (D) Representative WB result of Bax. (E) WB result of GAPDH, the loading control for A, B, C and D. Supplementary Figure 2. Matrine induced the expression and phosphorylation of STAT3 in NKTCL cells. NK92 cells (5 × 105) were treated with 1.96 mM matrine for 48 h, followed by western blot. NK92 cells treated without matrine were used as control. (A) Representative WB result of phosphorylation of STAT3 at Tyr705. (B) Representative WB result of STAT3. (C) WB result of GAPDH, the loading control for A and B. Supplementary Figure 3. Matrine decreased the expression of c-Myc protein in NKTCL cells. NK92 cells (5 × 105) were treated with 1.96 mM matrine for 48 h, followed by western blot. NK92 cells treated without matrine were used as control. (A) Representative WB result of c-Myc. (B) WB result of GAPDH, the loading control for A. Supplementary Figure 4. Matrine promoted c-Myc protein degradation in NKTCL cells. Cycloheximide chase assay was used for the half-time of c-Myc protein. NK92 cells (1 × 106) were treated with or without 1.96 mM matrine for 12 h. Cells were then treated with cycloheximide (100 μg/mL) for the indicated minutes, and western blotting was performed. NK92 cells treated without matrine were used as control. (A) Representative WB result of c-Myc in matrine treated NK92 cells. (B) WB result of GAPDH, the loading control for A. (C) Representative WB result of c-Myc in the control NK92 cells. (D) WB result of GAPDH, the loading control for C. Supplementary Figure 5. MG132 prevented matrine-induced c-Myc protein degradation in NKTCL cells. NK92 cells (5 × 105) were treated with 1.96 mM matrine, 10 μM MG132 with or without 1.96 mM matrine, respectively, for 6 h, followed by western blot. NK92 cells treated without matrine and MG132 were used as control. (A) Representative WB result of c-Myc. (B) WB result of GAPDH, the loading control for A. Supplementary Figure 6. Matrine inhibited NKTCL cells through CaMKIIγ/c-Myc pathway. NK92 cells (5 × 105) were treated with 1.96 mM matrine for 48 h, followed by western blot. NK92 cells treated without matrine were used as control. (A) Representative WB result of p-c-Myc (Ser62). (B) Representative WB result of c-Myc. (C) Representative WB result of CaMKIIγ. (D) Representative WB result of LMP1. (E) WB result of GAPDH, the loading control for A, B, C and D.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Gu, J., Zhang, Y., Wang, X. et al. Matrine inhibits the growth of natural killer/T-cell lymphoma cells by modulating CaMKIIγ-c-Myc signaling pathway. BMC Complement Med Ther 20, 214 (2020). https://doi.org/10.1186/s12906-020-03006-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12906-020-03006-2

Keywords

BMC Complementary Medicine and Therapies

ISSN: 2662-7671

Contact us