Skip to main content
Download PDF

Frankincense oil derived from Boswellia carteri induces tumor cell specific cytotoxicity

BMC Complementary and Alternative Medicine volume 9, Article number: 6 (2009) Cite this article

Abstract

Background

Originating from Africa, India, and the Middle East, frankincense oil has been important both socially and economically as an ingredient in incense and perfumes for thousands of years. Frankincense oil is prepared from aromatic hardened gum resins obtained by tapping Boswellia trees. One of the main components of frankincense oil is boswellic acid, a component known to have anti-neoplastic properties. The goal of this study was to evaluate frankincense oil for its anti-tumor activity and signaling pathways in bladder cancer cells.

Methods

Frankincense oil-induced cell viability was investigated in human bladder cancer J82 cells and immortalized normal bladder urothelial UROtsa cells. Temporal regulation of frankincense oil-activated gene expression in bladder cancer cells was identified by microarray and bioinformatics analysis.

Results

Within a range of concentration, frankincense oil suppressed cell viability in bladder transitional carcinoma J82 cells but not in UROtsa cells. Comprehensive gene expression analysis confirmed that frankincense oil activates genes that are responsible for cell cycle arrest, cell growth suppression, and apoptosis in J82 cells. However, frankincense oil-induced cell death in J82 cells did not result in DNA fragmentation, a hallmark of apoptosis.

Conclusion

Frankincense oil appears to distinguish cancerous from normal bladder cells and suppress cancer cell viability. Microarray and bioinformatics analysis proposed multiple pathways that can be activated by frankincense oil to induce bladder cancer cell death. Frankincense oil might represent an alternative intravesical agent for bladder cancer treatment.

Peer Review reports

Background

Frankincense resin is obtained from trees of the genus Boswellia (family Burseraceae). Incisions are made in the trunks of the trees to produce exuded gum, which appears as milk like resin. The resin hardens into orange-brown gum resin known as frankincense. There are numerous species and varieties of frankincense trees, including Boswellia serrata in India, Boswellia carteri in East Africa and China, Boswellia frereana in Somalia, and Boswellia sacra in Arabia, each producing a slightly different type of resin. Differences in soil and climate create more diversity in the resins, even within the same species. The aroma from these resins is valued for its presumed healing properties and superior qualities for religious rituals since the time of the ancient Egyptians [1], and has been used in incense, fumigants, and as a fixative in perfumes.

Frankincense resin has been considered throughout the ages to have a wealth of health supporting properties. The resins of Boswellia carteri and Boswellia serrata have been used for the treatment of rheumatoid arthritis and other inflammatory diseases [2] such as Crohn's disease [3] in traditional medicine of many countries. The anti-inflammatory activity has been attribute to the resin's ability in regulating immune cytokines production [4] and leukocyte infiltration [5, 6]. Boswellia serrata extract also exhibits anti-bacterial and anti-fungal activities [7]. Additionally, extracts from Boswellia species gum resins might possess anti-cancer activities, based on their anti-proliferative and pro-apoptotic activities in rat astrocytoma cell lines [8] and in human leukemia cell lines [9], as well as their anti-carcinogenic activity in chemically induced mouse skin cancer models [10]. Clinically, extract from the resin reduces the peritumoral edema in glioblastoma patients [8] and reverses multiple brain metastases in a breast cancer patient [11]. These results suggest that frankincense resin contains active ingredients that modulate important biological activities.

In search of the active medicinal ingredients of frankincense resins, Chevrier et al. reported that ethanol extract of Boswellia carteri resin comprises 7 boswellic acids [4]. Akihisa et al. reported that methanol extract of Boswellia carteri resin consists of 15 triterpene acids, including boswellic acids, and 2 cembrane-type diterpenes [12]. 11-keto-β-boswellic acid, the most potent anti-inflammatory component of the resin, selectively blocks leukotriene biosynthesis through inhibiting 5-lipoxygenase activity in rat neutrophilic granulocytes [13] and provides protective effects in a chemically induced mouse ulcerative colitis model [14]. Boswellic acids also prevent endotoxin/galactosamine-induced hepatitis in mice [15]. In addition, boswellic acids have been shown to possess anti-cancer activities through their cytostatic and apoptotic effects in multiple human cancer cell lines including meningioma cells [16], leukemia cells [17], hepatoma cells [18], melanoma cells, fibrosarcoma cells [19], and colon cancer cells [20].

Frankincense oil, an extract prepared by steam distillation from frankincense gum resin, is one of the most commonly used oils in aromatherapy practices. There has been considerable work done on the composition of frankincense oil from different species and commercial brands; and the constituents of frankincense oil differ according to the climate, harvest conditions, and geographical sources of frankincense resins [21]. Due to the contribution of boswellic acids, it is possible that frankincense oil also holds anti-cancer and anti-neoplastic properties. In this study, we demonstrated that a commercial source of frankincense oil can discriminate bladder cancer J82 cells from normal bladder urothelial UROtsa cells and suppress cancer cell viability. Based on gene expression analysis, frankincense oil activated several anti-proliferative and pro-apoptotic pathways that might be responsible for frankincense oil-induced cell death in J82 cells.

Methods

Reagents and chemicals

Cell culture medium [MEM and DMEM/F-12 (1:1)], fetal bovine serum (FBS), MEM vitamin solution, non-essential amino acids, epidermal growth factor (EGF), insulin-transferrin-sodium selenite (ITS) media supplement, sodium pyruvate, and penicillin-streptomycin were purchased from Invitrogen (Grand Island, NY). Frankincense oil containing 1,200 mg/ml frankincense gum resin was obtained from Young Living Essential Oils (Lehi, UT). XTT cell proliferation assay and in situ cell death detection kits were obtained from Roche (Indianapolis, IN). Trypan blue was purchased from Sigma (St. Louis, MO). RNeasy® Mini Kit was obtained from Qiagen (Valencia, CA).

Human bladder cell lines

Bladder transitional cell carcinoma J82 was obtained from ATCC (HTB-1; Manassas, VA). The J82 cell line was derived from a poorly differentiated, invasive human transitional cell bladder carcinoma (stage 3) [22]. J82 cells were maintained in growth medium consisting of MEM supplemented with 10% FBS, 0.1 mM non-essential amino acids, 1 mM sodium pyruvate, 2% MEM vitamin solution, 100 units/ml penicillin, and 100 μg/ml streptomycin. The UROtsa cell line was originally isolated from a primary culture of normal human urothelium and immortalized with a construct containing the SV40 large T antigen [23]. UROtsa cells were cultured in DMEM/F12 supplemented with 10 ng/ml EGF, 1× ITS media supplement, 100 units/ml penicillin, and 100 μg/ml streptomycin. Cells were maintained in a humidified cell incubator at 37°C and 5% CO2 and passaged every 3–4 days or when cells reached about 80% confluence.

Cell viability analysis

To determine number of viable cells following frankincense oil treatment, J82 and UROtsa cells were seeded in 96-well tissue culture plates at the density of 1 × 104 cells/mm2 in 100 μl growth medium. Following overnight incubation for adherence, 100 μl cell growth media or varying dilutions of frankincense oil (at1:600 to 1:4,000 final concentration) in their growth media were added to each well in triplicate to make a total of 200 μl. Cell viability was determined at the time of treatment and at 24 hours following frankincense oil exposure using the XTT cell proliferation assay kit. Briefly, at the time of assay, 100 μl growth media were removed from each well, and an aliquot of 50 μl XTT labeling mixture was added back to each well. Reactions were performed at 37°C for 4 hours. Absorbance was obtained by reading the plates at 450 nm wavelength using μQuant microplate reader (Bio-Tek; Winooski, VT). Absorbance values obtained at 24 hours for untreated and frankincense oil-treated cells were normalized to the values obtained at the time of treatment to calculate fold changes in cell survival.

Trypan blue exclusion was also included to determine cell viability following frankincense oil treatment. Briefly, J82 and UROtsa cells were seeded in 24-well tissue culture plates at the same density as used in XTT assay in 500 μl growth media. Following adherence, cells received either 500 μl of growth medium or varying dilutions of frankincense oil in each well. At 3 hours after frankincense oil treatment, the culture medium was collected to save the non-adherent cells; and the remaining cells were trypsinized and combined with the cells harvested from the culture medium. The cells were collected by centrifugation, and re-suspended with 200 μl phosphate buffered saline (PBS). Then, an aliquot (20 μl) of the cell suspension was mixed with the same volume of 0.4% (w/v) trypan blue solution. The cells were counted using a hemocytometer to determine the numbers of blue cells (non-viable) and bright cells (viable). Cell viability was expressed as the percentage of trypan blue positive cells compared to the total number of cells.

RNA extraction and quality evaluation

Total RNA was isolated from J82 cells for microarray analysis. Briefly, 2 × 105 J82 cells were seeded in 60 mm tissue culture plates, cultured overnight for adherence, and either left untreated or treated with 1:1,000 dilutions of frankincense oil in growth medium. Total RNA was isolated at 0 hours (no treatment) and at 0.5, 1, 2, and 3 hours after stimulation using the RNeasy® Mini total RNA isolation kit based on manufacture's recommendations (Qiagen; Valencia, CA). Total RNA concentration was determined with a nanodrop scanning spectrophotometer, and then qualitatively assessed for degradation using the ratio of 28:18s rRNA by a capillary gel electrophoresis system (Agilent 2100 Bioanalyzer, Agilent Technologies; Santa Clara CA).

RNA labeling, microarray hybridization, and scanning

A total of 250 ng of RNA from each time point was labeled using the Illumina Total Prep RNA Amplification Kit following manufacturer's directions (Ambion; Austin. TX). Briefly, cDNA was reverse transcribed from RNA after priming with T7-oligo-dT, and cRNA was synthesized in vitro from the T7 promoter while incorporating biotinylated UTP. cRNA was hybridized overnight to Illumina human Ref-8 version 3 BeadChips containing probes for a total of 24,526 transcripts. Microarray chips were washed to high stringency and labeled with streptavidin -Cy3 (Amersham Biosciences; Piscataway, NJ) prior to scanning on an Illumina BeadArray Reader.

Bioinformatics data analysis

Non-normalized fluorescent intensity of each probe on the microarray slide was obtained using the DirectHyb gene expression package in BeadStudio software (Illumina, version 3.1.3). Fluorescent intensity filtering was performed to remove genes that lacked a minimum relative fluorescence of 64 units in at least one time point. Data from the remaining probes were log transformed and quantile normalized (Matlab). A final filtering was performed to identify genes with a minimum two-fold change in normalized expression values between adjacent time points. Expression data are available on Gene Expression Omnibus (GEO) with accession number GSE14002.

TUNEL (terminal deoxynucleotidyl transferase dUTP nick end labeling) Analysis

TUNEL analysis was performed in J82 cells using an immunohistochemical (IHC)-like staining procedure as we reported [24]. Briefly, adherent J83 cells were either left untreated or treated with 1:1,000 dilution of frankincense oil. At 3 hours after treatment, both non-adherent and adherent and cells were collected following centrifugation. Cell pellets were fixed in 10% formalin, immersed in 2% agarose, and subjected to paraffin embedding. The embedded cell blocks were sectioned, dewaxed, and rehydrated. Apoptotic cells were detected using the in situ cell death detection kit. Following the terminal deoxynucleotidyl transferase reaction, fast red substrate was added for color development. Slides were then washed and sealed with an aqueous mounting medium.

DNA fragmentation analysis

To determine whether J82 cells undergo DNA fragmentation following frankincense oil treatment, 2 × 105 J82 and UROtsa cells were seeded in 60 mm tissue culture plates in their growth media, incubated overnight for adherence, and treated with a 1:1,000 dilution of frankincense oil in growth media. Cells were harvested at 0 (untreated control), 1, 3, and 6 hours following treatment; and genomic DNA was prepared and precipitated based on reported procedures [25]. Quantities of the genomic DNA were determined spectrophotometrically. Aliquots (10 μg) of the genomic DNA were separated on a 2% agarose gel; images of ethidium bromide stained gels were captured by the Gel Doc 100 system (Bio-Rad, Hercules, CA).

Statistics

The results are expressed as mean ± SEM from four experiments. Comparisons of J82 and UROtsa cell survival following frankincense oil treatment were made using the one-way analysis of variance (ANOVA) followed by post hoc Dunnett's test. P < 0.05 was considered statistically significant.

Results

Frankincense oil-suppressed bladder cell viability

The bladder carcinoma J82 cells presented a density-independent growth and grew in soft agar, but were not tumorigenic in nude mice [26]. The immortalized bladder urothelial UROtsa cells expressed SV40 large T antigen, but did not acquire characteristics of neoplastic transformation, including growth in soft agar or the development of tumors in nude mice [23]. To determine if frankincense oil suppresses bladder cell viability, both J82 and UROtsa cells were subjected to morphological evaluation and cell viability assessment. J82 cells underwent significant morphological changes, such as detaching from tissue culture plates and shrinking beginning within 3 hours following frankincense oil exposure. At 24 hours after treatment, J82 cells completely detached from tissue culture plates whereas untreated controls remained adherent to the plates (Figure 1A and 1B). In contrast, UROtsa cells remained attached to the bottom of plates and did not show noticeable morphological alterations (Figure 1C and 1D).

Figure 1
figure1

Morphological changes of bladder carcinoma J82 and bladder urothelial UROtsa cells following frankincense oil stimulation. Bladder J82 and UROtsa cells were seeded in 96-well tissue culture plates at the concentration of 1 × 104 cells/mm2, cultured overnight for adherence, and either left untreated or subjected 1:1,000 dilution of frankincense oil stimulation. Images were taken at 24 hours following treatments for (A) untreated J82 cells, (B) J82 cells treated with frankincense oil, (C) untreated UROtsa cells, and (D) UROtsa cells treated with frankincense oil using Olympus IX51 inverted microscope. Notice cell shrinkage observed in J82 cells following frankincense oil treatment. In contrast, UROtsa cells did not experience noticeable morphological alteration following the same concentration of frankincense oil exposure.

Full size image

To determine whether frankincense oil affects J82 and UROtsa cell viability, the number of viable J82 and UROtsa cells was determined following various dilutions (1:600 to 1:1,400) of frankincense oil exposure. In untreated controls, number of viable J82 cells and UROtsa cells increased 1.62 ± 0.31 and 2.72 ± 0.85 fold at 24 hours following cell seeding, respectively (Figure 2). Both J82 and UROtsa cells responded to frankincense oil treatment in a dose-dependent manner. J82 cell viability decreased when cells were treated with increasing concentrations of frankincense oil. No viable J82 cells remained at 24 hours after treatment with 1:1,100 dilution of frankincense oil (0.47 ± 0.43). In contrast, UROtsa cell viability was not significantly affected by the increasing concentrations of frankincense oil until 1:600 dilution was applied to the cells. When UROtsa cells were treated with 1:600 dilution of frankincense oil, cell viability decreased to 1.29 ± 0.77 fold as compared to untreated cells. No viable UROtsa cells were detected when the concentration of frankincense oil concentration increased to 1:400 (data not shown). Based on the XTT assay, IC50 values (the 50% inhibitory concentrations of frankincense oil) for J82 and UROtsa cells were 1:600 and 1:1,250, respectively. Trypan blue exclusion produced results similar to the XTT assay, except that J82 and UROtsa cells seem to be more sensitive to frankincense oil treatment at 1:1,300 and 1:600 dilutions, respectively (Figure 2B).

Figure 2
figure2

Bladder cell survival in response to frankincense oil exposure. Cell viability was determined using (A) a colometric XTT assay at 24 hours and (B) trypan blue exclusion at 3 hours after frankincense oil stimulation. All experiments were prepared in triplicate for XTT assay and duplicate for trypan blue exclusion. Data were presented as mean ± standard error of mean (SEM) from at least 3 independent experiments. * indicates statistical difference between frankincense oil-treated J82 cells and UROtsa cells (P < 0.05).

Full size image

Identification of frankincense oil-activated gene expression

To determine the nature of J82 cell death, microarray analysis was performed. Of the 24,526 gene probes on the microarray, 8,430 probes had a fluorescent intensity value of at least twice the background intensity for one or more time points under evaluation. A total of 122 genes in J82 cells were increased above two-fold in at least two adjacent time points by frankincense oil (see Additional file 1). Only 3 of these genes increased within the first 30 min [zinc finger protein 57, the small nucleolar RNA C/D box 48, and early growth response 1 gene (EGR1)]. Levels of EGR1 mRNA increased 5.87-fold within the first 30 min, and another 2.86-fold over the next 30 min. Another 15 genes increased at least two-fold between 30 and 60 min after frankincense oil stimulation; and 11 of which continued to show elevated expression beyond the first hour. A much larger number of genes increased in expression between 1 and 2 and between 2 and 3 hours following frankincense oil exposure.

A total of 47 genes were down-regulated in J82 cells by frankincense oil (see Additional file 2). Three genes [tubulin gamma 1, vacuolar protein sorting 11 homolog, and RNA polymerase II (DNA directed) polypeptide K] were the first to decrease greater than two-fold between 30 and 60 min after frankincense exposure. Another 12 genes decreased between 1 and 2 hours, and 32 other genes decreased between 2 and 3 hours. An additional 12 genes were identified whose levels of expression changed at least two-fold between two adjacent time points and then changed in the opposite direction at least two-fold between the next adjacent time points. These 12 genes were ankyrin repeat domain 27, chromosome 5 open reading frame 34, calcineurin binding protein 1, dodecenoyl-coenzyme A delta isomerase, dynein (axonemal, intermediate polypeptide 2), ATG2 [autophagy related 2 homolog A, N-deacetylase/N-sulfotransferase (heparan glucosaminyl) 2], oviductal glycoprotein 1, plexin A3, somatostatin receptor 1, transcriptional variant 1 of rinucleotide repeat containing 5, and zinc finger and BTB domain containing 11.

Functional grouping of frankincense oil-regulated genes

The gene products that were altered in frankincense oil-treated bladder carcinoma J82 cells were functionally grouped according to Gene Ontology classification. Based on the biological functions, gene products that function as cytokines, membrane receptors, enzymes (including kinases, peptidases, and phosphatases), and molecular transport were identified and listed in Table 1. A complete list of genes under each classification is provided in Additional file 3. Frankincense oil-regulated gene products that function as transcription factors, cell cycle arrest and cell proliferation, as well as apoptotic factors showed that frankincense oil induces cell cycle arrest and apoptosis in J82 cells.

Table 1 Functional groups of frankincense oil-regulated genes in bladder cancer J82 cells
Full size table

Transcription regulators

Two transcription factors, LOC12629 and EGR1, were immediately (within 30 min) up-regulated by frankincense oil (Table 2). Another 5 transcription factors, including ATF3, FOS, FOSB, KLF2, and ZNF234 were up-regulated within 1 hour and sustained for at least 2 hours following frankincense oil treatment. Three transcription factors, KLF4, KLF5, and ZBTB11 were up-regulated by frankincense oil between 1 and 2 hours post-treatment. The remaining 11 transcription factors (DDIT3, DEDD2, DENR, HES1, ID1, JUN, JUNB, SNAPC1, TSC22T1, UBTF, ZFP36) were considered to be late responders because their expression was altered after 2 hours of frankincense oil exposure.

Table 2 Frankincense oil-regulated transcription factors in J82 cells
Full size table

Cell cycle arrest and cell proliferation

Several gene products identified as frankincense oil responsive genes were negatively associated with regulation of cell proliferation and positively associated with cell cycle arrest (Table 3). Genes that have been identified as anti-proliferative genes, including IL8, CLK1, DLG1, KLF4, NEDD9, CDKN1A, IL1A, IL6, and SNFILK were up-regulated in frankincense oil-treated J82 cells. In addition, DDIT3 along with IL8 and CDKNIA being identified to be responsible for cell cycle arrest were also up-regulated by frankincense oil. In contrast, H2AFX and HDAC4, genes that are responsible for DNA repair and cell cycle progression, were suppressed by frankincense oil. Other anti-proliferative genes, including SSTR1, IL1A, and IL6, were also up-regulated between 0.5 and 2 hours upon frankincense oil stimulation.

Table 3 Frankincense oil-regulated growth inhibitory genes in J82 cells
Full size table

Apoptosis

Levels of a large number of genes that are responsible for apoptosis were found to be modulated by frankincense oil (Figure 3). These genes included CDKN1A, DEDD2, IER3, IL6, SGK, and TNFAIP3 (up-regulated between 1 and 2 hours and remained up-regulated) as well as GAD45B, NUDT2, and others (up-regulated between 2 and 3 hours). In addition, the cell survival gene, AXL, was down-regulated by frankincense oil. However, two anti-apoptotic genes, GSTP1 and IL1A, were up-regulated. A similar contradiction was seen with a pro-apoptotic gene, ING4, being down-regulated.

Figure 3
figure3

Hierarchical clustering of frankincense oil-regulated apoptosis-related genes in J82 cells. The map was obtained using Biometric Research Branch (BRB) ArrayTools version 3.4.0 – Beta_2 software http://linus.nci.nih.gov/BRB-ArrayTools.html after log2 transformation of fluorescence intensities. Each column represents time intervals following frankincense oil exposure, and each row represents a gene probe set. The expression levels for individual genes are indicated by green/red color indicating an elevated/suppressed level of expression, respectively.

Full size image

Frankincense oil-induced cell death

TUNEL analysis was performed to determine whether frankincense oil treated J82 cells undergo apoptosis. Frankincense oil treatment resulted in an increased number of bright red colored TUNEL positive cells as compared to untreated cells (Figure 4A). Genomic DNA fragmentation was determined between hours 1 and 6 in J82 cells following frankincense oil treatment. Agarose gel electrophoresis results showed that all genomic DNA remained as large pieces of DNA without forming a small DNA ladder (Figure 4B). There was no detectable genomic DNA for J82 cells harvested at 12 hours following frankincense oil treatment (data not shown).

Figure 4
figure4

Frankincense oil-induced J82 cell death. To determine whether frankincense oil-induced apoptosis in bladder cancer cells, J82 cells were seeded in 60 mm tissue culture plates at the concentration of 2 × 105 cells per plate, cultured overnight for adherence, and either left untreated or treated with 1:1,000 dilution of frankincense oil. (A) TUNEL analysis was performed at 3 hours following treatment. Apoptotic cells with damaged DNA were stained positive with a bright red color (inserts). (B) DNA fragmentation was determined by separating genomic DNA on a 2% agaorse gel; and the gel image was captured using Gel Doc 100 system (Bio-Rad, Hercules, CA).

Full size image

Discussion

Ranging from herbs to acupuncture, alternative medicine is becoming increasingly popular for managing health-related issues. In this brief report, we described that frankincense oil, with a window of concentration, specifically suppressed cell viability in human bladder carcinoma J82 cells, but did not affect cell viability in immortalized normal urothelial UROtsa cells. Frankincense oil suppressed J82 cell viability may be attributed to the activation of growth arrest and pro-apoptotic genes. The possibility that the witnessed differences in cell survival could be due to the presence of EGF in the media was considered. Despite this concern, we expected that UROtsa cells' resistance to frankincense oil as compared to J82 cells may not result from the presence of EGF in their growth medium based on the following observations: first, a higher concentrations of EGF (i.e. 1 μg/ml) had been reported for apoptosis protection [27], second, UROtsa cells cultured in EGF-free medium were less sensitivity to frankincense oil-suppressed cell viability as compared to J82 cells, although the overall UROtsa cell viability was reduced (data not shown), and third, four other oils, including sandalwood oil (Santalum album), balsam fir oil (Abies balsamea), palo santo oil (Bursera graveolens), and tsuga oil (Tsuga canadensis) (Young Living Essential Oils), induced nearly identical cytotoxicity in both J82 and UROtsa cells (data not shown).

Commercial frankincense oil was directly applied in our experimental cell models without modification. It was not our intention in this preliminary study to dissect the specific chemical composition of frankincense oil nor determine its efficacious dosage, since some reports indicated that total frankincense extract is more potent than pure, specific boswellic acids [9]. However, a standard assessment between chemical composition of frankincense oil and its efficacy in tumor suppression will be required in future clinical trials. In addition, frankincense oil was directed added to cell culture media in this study without the inclusion of a carrier; and dose-dependent suppression in cell viability was observed in both J82 and UROtsa cells in the absence of any oil carrier. The absence of carrier eliminated carrier-dependent effects of frankincense resin extract as reported by Chevrier et al. [4].

Gene expression analysis was terminated within 3 hours following frankincense oil treatment, since isolated RNA quality and quantity were not sufficient for microarray analysis beyond this time point. We reported a total of 122 up- and 47 down-regulated genes with greater than 2-fold induction or suppression over the period of 3 hours. These findings suggested very specific actions of frankincense oil. The genome-wide gene expression analysis supports patterns of stress, activation of cell cycle arrest, suppression of cell proliferation, and activation of apoptotic signaling in frankincense oil-treated J82 cells within 1 hours of stimulation, and some of these processes were sustained for 3 hours.

Based on the temporal regulation of the genes identified by microarray and bioinformatics analysis, we proposed that frankincense oil induces various death pathways in J82 cells. Waves of transcription factors were regulated by frankincense oil from between 30 min and 3 hours. EGR1 was one of the few genes that were rapidly up-regulated within the first 30 min. Although EGR1 has been shown to be an early gene that is immediately up-regulated in other systems and is correlated with DNA repair [28], the mechanism for elevated EGR1 expression by frankincense oil is unclear. EGR1 has been reported to increase transcription of another transcription factor ATF3 [29]; and levels of ATF3 were increased in our system between 30 and 60 min after frankincense exposure. ATF3 is induced by stresses and can bind to DDIT3 [30], which was up-regulated between 2 and 3 hours in our system. In addition, DDIT3 responds to DNA damage [31], and is responsible for cell cycle arrest [32]. Moreover, using the PAINT webtool http://www.dbi.tju.edu/dbi/tools/paint/[33] to search the TRANSFAC database, we identified EGR1 and ATF3 binding sequences upstream of multiple apoptosis-related genes identified in our system. EGR1 can bind 5'-flanking regions of 5 of the identified apoptosis-related genes (FGFR1, GADD45B, HES1, RHOB, and TRIB1) that were up-regulated between 2 and 3 hours. Binding sites for ATF3 were found in the 5'-flanking regions of 3 apoptosis-related genes (FOSB, GEM, and LAMA5) that were up-regulated between 1 and 2 hours, and 3 additional genes (HSPA1A, ID1, and JUN) between 2 and 3 hours. By similar inference, ATF3 may also account for the expression of 2 down-regulated genes: ING4 (1–2 hours) and ATG5 (2–3 hours). The sequential expression of these identified transcription factors may be ultimately responsible for cell cycle arrest, suppressed cell proliferation, and apoptosis in frankincense oil-treated J82 cells.

Suppressed cell viability and proliferation in frankincense oil-treated J82 cells was also confirmed by elevated expression of genes that are responsible for cell cycle arrest and anti-proliferation. Up-regulated IL8 [34] and CDKN1A [35] have been shown to be responsible for cell cycle arrest and suppressed cell proliferation. IL1A is a negative regulator for cell cycle progression and cell proliferation [36]; and IL6 has been shown to induce growth arrest [37]. However, levels of ING4, a molecule that is a negative regulation of cell proliferation in a human hepatocellular liver carcinoma cell line (HepG2) [38], were up-regulated in frankincense oil-treated cells. ING4 may function differently between J82 cells and HepG2 cells; or ING4 activity is suppressed by a large number of pro-apoptotic molecules in response to frankincense oil.

Frankincense oil up-regulated several pro-apoptotic genes, including CDKN1A [39], DEDD2 [40], NUDT2 [41], SGK, TNFAIP3, and IER3 [42]. Elevated expression (between 0.5 and 2 hours) followed by suppressed expression (2–3 hours) of a cell survival membrane receptor, AXL [43], suggests that cells may try to prolong cell survival following frankincense oil exposure. We also propose that frankincense oil might activate both extrinsic and intrinsic death signaling in J82 cells through stress and death receptors, respectively, to execute apoptosis. Intrinsic death signaling was suggested by up-regulated expression of SSRT1, GADD45B, DDIT3, CDKN1A that have been shown to be required for DNA damage-induced cell cycle arrest [35, 39]. Extrinsic death signaling was implicated by the up-regulated expression of DEDD2, which is a death domain receptors and induces apoptosis [40].

Although the bioinformatics and TUNEL analyses reported here suggested that frankincense oil induced apoptosis, rather than necrosis, in J82 cells, frankincense oil did not cause DNA fragmentation, a hallmark of apoptosis, in this bladder cancer cell line. It is possible that DNA fragmentation occurred between 6 and 12 hours post frankincense oil treatment. Alternatively, apoptosis without DNA fragmentation has been reported in several occasions [44–46]; and frankincense oil-induced J82 cell death may fit in this category. The detailed molecular and biological pathways utilized by frankincense oil in inducing bladder cancer cell specific cell death require further investigation.

This study helps to show that frankincense oil may be appropriate as an alternative therapy for bladder cancer. This is the first report demonstrating that frankincense oil can discriminate between bladder cancer cells and normal urothelial cells in a cell culture system, and utilizing microarray technology to identify potential biological pathways activated by frankincense oil. Our results are consistent with a news report that frankincense oil specifically targets malignant melanoma but not normal skin cells in horses http://www.purepeace.com/press/press_frankincense.pdf. Future studies are required to determine whether frankincense oil has similar effects on other bladder cancer cell lines of varying severity such as RT4, T24, and 5637, followed by in vivo studies using bladder cancer animal models. In addition, a standard manufacturing and indication needs to be applied before the commercial frankincense oil can be used as an alternative or complementary therapy for treating bladder cancer.

Conclusion

Frankincense oil can discriminate bladder cancer cells and normal urothelial cells in culture. The oil suppresses cell survival and induces apoptosis in cultured bladder cancer cells. Based on this preliminary observation, frankincense oil may represent an alternative intravesical therapy for bladder cancer, although more bladder cancer cell lines and animal models need to be used to confirm current observations.

References

  1. 1.

    Maloney GA: Gold, frankincense, and myrrh : an introduction to Eastern Christian spirituality. 1997, New York: Crossroads Pub. Co

    Google Scholar 

  2. 2.

    Banno N, Akihisa T, Yasukawa K, Tokuda H, Tabata K, Nakamura Y, Nishimura R, Kimura Y, Suzuki T: Anti-inflammatory activities of the triterpene acids from the resin of Boswellia carteri. J Ethnopharmacol. 2006, 107: 249-253. 10.1016/j.jep.2006.03.006.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Langmead L, Rampton DS: Review article: complementary and alternative therapies for inflammatory bowel disease. Aliment Pharmacol Ther. 2006, 23: 341-349. 10.1111/j.1365-2036.2006.02761.x.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Chevrier MR, Ryan AE, Lee DY, Zhongze M, Wu-Yan Z, Via CS: Boswellia carterii extract inhibits TH1 cytokines and promotes TH2 cytokines in vitro. Clin Diag Lab Immunol. 2005, 12: 575-580. 10.1128/CDLI.12.5.575-580.2005.

    CAS  Google Scholar 

  5. 5.

    Sharma ML, Khajuria A, Kaul A, Singh S, Singh GB, Atal CK: Effect of salai guggal ex-Boswellia serrata on cellular and humoral immune responses and leucocyte migration. Agents Actions. 1988, 24: 161-164. 10.1007/BF01968095.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Singh GB, Atal CK: Pharmacology of an extract of salai guggal ex-Boswellia serrata, a new non-steroidal anti-inflammatory agent. Agents Actions. 1986, 18: 407-412. 10.1007/BF01965005.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Weckesser S, Engel K, Simon-Haarhaus B, Wittmer A, Pelz K, Schempp CM: Screening of plant extracts for antimicrobial activity against bacteria and yeasts with dermatological relevance. Phytomedicine. 2007, 14: 508-516. 10.1016/j.phymed.2006.12.013.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Winking M, Sarikaya S, Rahmanian A, Jodicke A, Boker DK: Boswellic acids inhibit glioma growth: a new treatment option?. J Neurooncol. 2000, 46: 97-103. 10.1023/A:1006387010528.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Hostanska K, Daum G, Saller R: Cytostatic and apoptosis-inducing activity of boswellic acids toward malignant cell lines in vitro. Anticancer Res. 2002, 22: 2853-2862.

    CAS  PubMed  Google Scholar 

  10. 10.

    Huang MT, Badmaev V, Ding Y, Liu Y, Xie JG, Ho CT: Anti-tumor and anti-carcinogenic activities of triterpenoid, β-boswellic acid. BioFactors. 2000, 13: 225-230. 10.1002/biof.5520130135.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Flavin DF: A lipoxygenase inhibitor in breast cancer brain metastases. J Neurooncol. 2007, 82: 91-93. 10.1007/s11060-006-9248-4.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Akihisa T, Tabata K, Banno N, Tokuda H, Nishimura R, Nakamura Y, Kimura Y, Yasukawa K, Suzuki T: Cancer chemopreventive effects and cytotoxic activities of the triterpene acids from the resin of Boswellia carteri. Biol Pharm Bull. 2006, 29: 1976-1979. 10.1248/bpb.29.1976.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Safayhi H, Sailer ER, Ammon HP: Mechanism of 5-lipoxygenase inhibition by acetyl-11-keto-β-boswellic acid. Mol Pharmacol. 1995, 47: 1212-1216.

    CAS  PubMed  Google Scholar 

  14. 14.

    Anthoni C, Laukoetter MG, Rijcken E, Vowinkel T, Mennigen R, Muller S, Senninger N, Russell J, Jauch J, Bergmann J: Mechanisms underlying the anti-inflammatory actions of boswellic acid derivatives in experimental colitis. Am J Physiol Gastrointest Liver Physiol. 2006, 290: G1131-1137. 10.1152/ajpgi.00562.2005.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Safayhi H, Mack T, Ammon HP: Protection by boswellic acids against galactosamine/endotoxin-induced hepatitis in mice. Biochem Pharmacol. 1991, 41: 1536-1537. 10.1016/0006-2952(91)90575-P.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Park YS, Lee JH, Bondar J, Harwalkar JA, Safayhi H, Golubic M: Cytotoxic action of acetyl-11-keto-β-boswellic acid (AKBA) on meningioma cells. Planta Medica. 2002, 68: 397-401. 10.1055/s-2002-32090.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Shao Y, Ho CT, Chin CK, Badmaev V, Ma W, Huang MT: Inhibitory activity of boswellic acids from Boswellia serrata against human leukemia HL-60 cells in culture. Planta Medica. 1998, 64: 328-331. 10.1055/s-2006-957444.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Liu JJ, Nilsson A, Oredsson S, Badmaev V, Duan RD: Keto- and acetyl-keto-boswellic acids inhibit proliferation and induce apoptosis in Hep G2 cells via a caspase-8 dependent pathway. Int J Mol Med. 2002, 10: 501-505.

    CAS  PubMed  Google Scholar 

  19. 19.

    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 Detec Prev. 2003, 27: 67-75. 10.1016/S0361-090X(02)00170-8.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Liu JJ, Nilsson A, Oredsson S, Badmaev V, Zhao WZ, Duan RD: Boswellic acids trigger apoptosis via a pathway dependent on caspase-8 activation but independent on Fas/Fas ligand interaction in colon cancer HT-29 cells. Carcinogenesis. 2002, 23: 2087-2093. 10.1093/carcin/23.12.2087.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Mikhaeil BR, Maatooq GT, Badria FA, Amer MM: Chemistry and immunomodulatory activity of frankincense oil. Z Naturforsch [C]. 2003, 58 (3-4): 230-238.

    CAS  Article  Google Scholar 

  22. 22.

    O'Toole C, Price ZH, Ohnuki Y, Unsgaard B: Ultrastructure, karyology and immunology of a cell line originated from a human transitional-cell carcinoma. Br J Cancer. 1978, 38: 64-76.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Petzoldt JL, Leigh IM, Duffy PG, Sexton C, Masters JR: Immortalisation of human urothelial cells. Urol Res. 1995, 23: 377-380. 10.1007/BF00698738.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Yang Q, Titus M, Fung KM, Lin HK: 5α-Androstane-3α, 17β-diol supports human prostate cancer cell survival and proliferation through androgen receptor-independent signaling pathways: Implication of androgen-independent prostate cancer progression. J Cell Biochem. 2008, 104: 1612-1624. 10.1002/jcb.21731.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Mondalek FG, Lawrence BJ, Kropp BP, Grady BP, Fung KM, Madihally SV, Lin HK: The incorporation of poly(lactic-co-glycolic) acid nanoparticles into porcine small intestinal submucosa biomaterials. Biomaterials. 2008, 29: 1159-1166. 10.1016/j.biomaterials.2007.11.020.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Marshall CJ, Franks LM, Carbonell AW: Markers of neoplastic transformation in epithelial cell lines derived from human carcinomas. J Natl Cancer Inst. 1977, 58: 1743-1751.

    CAS  PubMed  Google Scholar 

  27. 27.

    Gibson EM, Henson ES, Haney N, Villanueva J, Gibson SB: Epidermal growth factor protects epithelial-derived cells from tumor necrosis factor-related apoptosis-inducing ligand-induced apoptosis by inhibiting cytochrome c release. Cancer Res. 2002, 62: 488-496.

    CAS  PubMed  Google Scholar 

  28. 28.

    Franken NA, Ten Cate R, Van Bree C, Haveman J: Induction of the early response protein EGR-1 in human tumour cells after ionizing radiation is correlated with a reduction of repair of lethal lesions and an increase of repair of sublethal lesions. Int J Oncol. 2004, 24: 1027-1031.

    CAS  PubMed  Google Scholar 

  29. 29.

    Bottone FG, Moon Y, Alston-Mills B, Eling TE: Transcriptional regulation of activating transcription factor 3 involves the early growth response-1 gene. J Pharmacol Exp Ther. 2005, 315: 668-677. 10.1124/jpet.105.089607.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Chen BP, Wolfgang CD, Hai T: Analysis of ATF3, a transcription factor induced by physiological stresses and modulated by gadd153/Chop10. Mol Cell Biol. 1996, 16 (3): 1157-1168.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Park JS, Luethy JD, Wang MG, Fargnoli J, Fornace AJ, McBride OW, Holbrook NJ: Isolation, characterization and chromosomal localization of the human GADD153 gene. Gene. 1992, 116: 259-267. 10.1016/0378-1119(92)90523-R.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Manthey KC, Rodriguez-Melendez R, Hoi JT, Zempleni J: Riboflavin deficiency causes protein and DNA damage in HepG2 cells, triggering arrest in G1 phase of the cell cycle. J Nutr Biochem. 2006, 17: 250-256. 10.1016/j.jnutbio.2005.05.004.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Vadigepalli R, Chakravarthula P, Zak DE, Schwaber JS, Gonye GE: PAINT: a promoter analysis and interaction network generation tool for gene regulatory network identification. OMICS. 2003, 7: 235-252. 10.1089/153623103322452378.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Thirumangalakudi L, Yin L, Rao HV, Grammas P: IL-8 induces expression of matrix metalloproteinases, cell cycle and pro-apoptotic proteins, and cell death in cultured neurons. J Alzheimers Dis. 2007, 11: 305-311.

    CAS  PubMed  Google Scholar 

  35. 35.

    Bunz F, Dutriaux A, Lengauer C, Waldman T, Zhou S, Brown JP, Sedivy JM, Kinzler KW, Vogelstein B: Requirement for p53 and p21 to sustain G2 arrest after DNA damage. Science. 1998, 282: 1497-1501. 10.1126/science.282.5393.1497.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Zeki K, Morimoto I, Arao T, Eto S, Yamashita U: Interleukin-1α regulates G1 cell cycle progression and arrest in thyroid carcinoma cell lines NIM1 and NPA. J Endocrinol. 1999, 160: 67-73. 10.1677/joe.0.1600067.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Sanford DC, DeWille JW: C/EBPdelta is a downstream mediator of IL-6 induced growth inhibition of prostate cancer cells. Prostate. 2005, 63: 143-154. 10.1002/pros.20159.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Zhang X, Xu LS, Wang ZQ, Wang KS, Li N, Cheng ZH, Huang SZ, Wei DZ, Han ZG: ING4 induces G2/M cell cycle arrest and enhances the chemosensitivity to DNA-damage agents in HepG2 cells. FEBS Lett. 2004, 570: 7-12. 10.1016/j.febslet.2004.06.010.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Levkau B, Koyama H, Raines EW, Clurman BE, Herren B, Orth K, Roberts JM, Ross R: Cleavage of p21Cip1/Waf1 and p27Kip1 mediates apoptosis in endothelial cells through activation of Cdk2: role of a caspase cascade. Mol Cell. 1998, 1: 553-563. 10.1016/S1097-2765(00)80055-6.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Roth W, Stenner-Liewen F, Pawlowski K, Godzik A, Reed JC: Identification and characterization of DEDD2, a death effector domain-containing protein. J Biol Chem. 2002, 277: 7501-7508. 10.1074/jbc.M110749200.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Vartanian A, Alexandrov I, Prudowski I, McLennan A, Kisselev L: Ap4A induces apoptosis in human cultured cells. FEBS Lett. 1999, 456: 175-180. 10.1016/S0014-5793(99)00956-4.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Pietzsch A, Buchler C, Aslanidis C, Schmitz G: Identification and characterization of a novel monocyte/macrophage differentiation-dependent gene that is responsive to lipopolysaccharide, ceramide, and lysophosphatidylcholine. Biochem Biophys Res Commun. 1997, 235: 4-9. 10.1006/bbrc.1997.6715.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Goruppi S, Ruaro E, Varnum B, Schneider C: Requirement of phosphatidylinositol 3-kinase-dependent pathway and Src for Gas6-Axl mitogenic and survival activities in NIH 3T3 fibroblasts. Mol Cell Biol. 1997, 17: 4442-4453.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Yamaguchi K, Uzzo R, Dulin N, Finke JH, Kolenko V: Renal carcinoma cells undergo apoptosis without oligonucleosomal DNA fragmentation. Biochem Biophys Res Commun. 2004, 318: 710-713. 10.1016/j.bbrc.2004.04.086.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Hirata H, Hibasami H, Yoshida T, Morita A, Ohkaya S, Matsumoto M, Sasaki H, Uchida A: Differentiation and apoptosis without DNA fragmentation in cultured Schwann cells derived from wallerian-degenerated nerve. Apoptosis. 1998, 3: 353-360. 10.1023/A:1009633205444.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Yuste VJ, Bayascas JR, Llecha N, Sanchez-Lopez I, Boix J, Comella JX: The absence of oligonucleosomal DNA fragmentation during apoptosis of IMR-5 neuroblastoma cells: disappearance of the caspase-activated DNase. J Biol Chem. 2001, 276: 22323-22331. 10.1074/jbc.M100072200.

    CAS  Article  PubMed  Google Scholar 

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6882/9/6/prepub

Download references

Acknowledgements

This work was supported in part by NIH grants RR16478, RR03025, RR15577, 01700172, and U2+AI062629.

Author information

Affiliations

  1. Arthritis and Immunology Research Program, Oklahoma Medical Research Foundation Microarray Research Facility, Oklahoma City, OK, 73104, USA

    Mark Barton Frank & Jeanette Osban

  2. Department of Urology, University of Oklahoma Health Sciences Center, Oklahoma City, OK, 73104, USA

    Qing Yang, Joseph T Azzarello, Richard A Ashley & Hsueh-Kung Lin

  3. Department of Physiology, University of Oklahoma Health Sciences Center, Oklahoma City, OK, 73104, USA

    Joseph T Azzarello, Marcia R Saban, Ricardo Saban & Hsueh-Kung Lin

  4. Department of Comparative Medicine, University of Oklahoma Health Sciences Center, Oklahoma City, OK, 73104, USA

    Jan C Welter

  5. Department of Pathology, University of Oklahoma Health Sciences Center, Oklahoma City, OK, 73104, USA

    Kar-Ming Fung

  6. Oklahoma City Veterans Affairs Medical Center, Oklahoma City, OK, 73104, USA

    Hsueh-Kung Lin

Authors
  1. Mark Barton Frank
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Qing Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jeanette Osban
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Joseph T Azzarello
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Marcia R Saban
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Ricardo Saban
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Richard A Ashley
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Jan C Welter
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Kar-Ming Fung
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Hsueh-Kung Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Hsueh-Kung Lin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MBF and JO conducted microarray and bioinformatics analysis. QY, JTA, and MRS performed cell biology and apoptosis analysis of cultured bladder cells. RS, RAA, JCW, KMF, and HKL conceived the idea, designed the experiments, and interpreted the experimental results. All authors contributed to manuscript preparations and approved the final manuscript.

Electronic supplementary material

12906_2008_231_MOESM1_ESM.doc

Additional file 1: Genes with minimum two-fold increase in adjacent time points. The data provided a list of all genes whose levels of expression are elevated at least two folds from one time point to the next time point. (DOC 83 KB)

12906_2008_231_MOESM2_ESM.doc

Additional file 2: Genes with minimum two-fold decrease in adjacent time points. The data listed all genes whose levels of expression are suppressed at least two folds from one time point to the next time point. (DOC 44 KB)

12906_2008_231_MOESM3_ESM.doc

Additional file 3: Functional groups of frankincense oil-regulated genes in bladder cancer J82 cells. The data provided the gene ontology classification for all genes that are regulated by frankincense oil. (DOC 50 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Frank, M.B., Yang, Q., Osban, J. et al. Frankincense oil derived from Boswellia carteri induces tumor cell specific cytotoxicity. BMC Complement Altern Med 9, 6 (2009). https://doi.org/10.1186/1472-6882-9-6

Download citation

Keywords

  • Bladder Cancer Cell
  • Bladder Cancer Cell Line
  • Boswellic Acid
  • Transitional Cell Bladder Carcinoma
  • UROtsa Cell

BMC Complementary Medicine and Therapies

ISSN: 2662-7671

Contact us