Differential effects of selective frankincense (Ru Xiang) essential oil versus non-selective sandalwood (Tan Xiang) essential oil on cultured bladder cancer cells: a microarray and bioinformatics study
Chinese Medicine volume 9, Article number: 18 (2014)
Frankincense (Boswellia carterii, known as Ru Xiang in Chinese) and sandalwood (Santalum album, known as Tan Xiang in Chinese) are cancer preventive and therapeutic agents in Chinese medicine. Their biologically active ingredients are usually extracted from frankincense by hydrodistillation and sandalwood by distillation. This study aims to investigate the anti-proliferative and pro-apoptotic activities of frankincense and sandalwood essential oils in cultured human bladder cancer cells.
The effects of frankincense (1,400–600 dilutions) (v/v) and sandalwood (16,000–7,000 dilutions) (v/v) essential oils on cell viability were studied in established human bladder cancer J82 cells and immortalized normal human bladder urothelial UROtsa cells using a colorimetric XTT cell viability assay. Genes that responded to essential oil treatments in human bladder cancer J82 cells were identified using the Illumina Expression BeadChip platform and analyzed for enriched functions and pathways. The chemical compositions of the essential oils were determined by gas chromatography–mass spectrometry.
Human bladder cancer J82 cells were more sensitive to the pro-apoptotic effects of frankincense essential oil than the immortalized normal bladder UROtsa cells. In contrast, sandalwood essential oil exhibited a similar potency in suppressing the viability of both J82 and UROtsa cells. Although frankincense and sandalwood essential oils activated common pathways such as inflammatory interleukins (IL-6 signaling), each essential oil had a unique molecular action on the bladder cancer cells. Heat shock proteins and histone core proteins were activated by frankincense essential oil, whereas negative regulation of protein kinase activity and G protein-coupled receptors were activated by sandalwood essential oil treatment.
The effects of frankincense and sandalwood essential oils on J82 cells and UROtsa cells involved different mechanisms leading to cancer cell death. While frankincense essential oil elicited selective cancer cell death via NRF-2-mediated oxidative stress, sandalwood essential oil induced non-selective cell death via DNA damage and cell cycle arrest.
Frankincense gum resins (known as Ru Xiang in Chinese) are obtained from Boswellia trees (family Burseraceae) and have been used for the treatment of rheumatoid arthritis and other inflammatory diseases  such as Crohn's disease . Extracts from Boswellia species resins exhibit anti-proliferative and pro-apoptotic activities in rat astrocytoma cell lines , human leukemia cell lines , and chemically induced mouse skin cancer models . The frankincense essential oil possesses anti-proliferative and pro-apoptotic activities against multiple human cancer cell lines in vitro and in vivo[6–8]. Boswellic acids were found to be the major components in frankincense extracts, with anti-tumor activity owing to their cytostatic and pro-apoptotic properties in multiple human cancer cell lines including meningioma cells , leukemia cells , hepatoma cells , melanoma cells, fibrosarcoma cells , colon cancer cells , and prostate cancer cells [14–16]. Some of the effects of frankincense essential oil were found to be related to the activities of sesquiterpenes and diterpenes .
Sandalwood (known as Tan Xiang in Chinese) belongs to the genus Santalum. The sandalwood essential oil has been used to treat skin diseases, acne, dysentery, and gonorrhea , and is considered an excellent sedating agent . Sandalwood essential oil also exhibits anti-bactericidal activity  and chemoprevention in chemically induced skin papillomas and skin cancer in CD1 mice [21, 22].
The genus Boswellia consists of four major species, while the genus Santalum consists of approximately 25 species. Each species of Boswellia and Santalum produces a slightly different type of aroma as a result of the soil and climate diversity. Frankincense and sandalwood essential oils are prepared from hydro- or steam-distillation of the plants . Medical applications of essential oils ranging from treatments for skin conditions to remedies for cancer have been based on the historical uses of these plant products. However, the active chemical compositions and mechanisms of action remain largely unclear. We previously showed that frankincense essential oil prepared from Boswellia species with different temperatures and durations of hydrodistillation possessed different chemical constituents and biological activities [7, 8].
This study aims to investigate the anti-proliferative and pro-apoptotic activities of frankincense and sandalwood essential oils on cultured human bladder cancer cells using microarrays and bioinformatics. We also intended to relate the cellular activity and gene expression profile to the chemical differences between frankincense and sandalwood essential oils.
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 Life Technologies (Grand Island, NY, USA). Frankincense (B. carterii from Somalia) and sandalwood (S. album from Sri Lanka) essential oils were obtained from Young Living Essential Oils (Lehi, UT, USA) and prepared based on previously reported procedures . An XTT cell proliferation assay kit was obtained from Roche Applied Science (Indianapolis, IN, USA). An RNeasy® Mini Kit was obtained from Qiagen (Valencia, CA, USA).
Chemical compositions of essential oils
The essential oil components were identified by gas chromatography–mass spectrometry (GC-MS), using an Agilent 7890A GC system (Agilent Technologies, Santa Clara, CA, USA) equipped with an Agilent 5975C mass selective detector (Agilent Technologies). We used either an HP-1 50 m × 0.32 mm ID × 0.5 μm film column (Agilent Technologies) or a DB-WAX 60 m × 0.32 mm ID × 0.5 μm film column (Agilent Technologies) for sandalwood essential oil, and an Rxi-5 ms 60 m × 0.25 mm ID × 0.25 μm column (Restek, Bellefonte, PA, USA) for frankincense essential oil. Retention indices were calculated by performing injections of C7-C30 normal alkanes (Supelco, Bellefonte, PA, USA) to confirm the MS identification. The detailed procedures for chemical analysis were reported previously .
As natural products, chemical constituents of frankincense and sandalwood essential oils varied from batch to batch and from season to season. A batch of distillation was used to establish a baseline relationship between chemical compositions and molecular pathways activation throughout the study.
Human bladder cell lines
Human bladder cancer J82 cells, derived from a poorly differentiated, invasive stage 3 transitional cell carcinoma , were obtained from ATCC (HTB-1; Manassas, VA, USA). 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, and 100 units/mL penicillin-100 μg/mL streptomycin. The UROtsa cell line was isolated from a primary culture of normal human urothelium and immortalized with the SV40 large T antigen . UROtsa cells were provided by Dr. Ricardo Saban (Department of Physiology, University of Oklahoma Health Sciences Center) and cultured in DMEM/F12 supplemented with 10 ng/mL EGF, 1 × ITS media supplement, and penicillin-streptomycin. Both cell lines 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 assay
Bladder cancer J82 cells and immortalized normal UROtsa cells were seeded in 96-well tissue culture plates at a density of 5 × 103 cells/well in 100 μL of their growth media to determine cell viability following treatment with frankincense and sandalwood essential oils. Following overnight incubation and adherence, an additional aliquot of growth media (100 μL) or varying dilutions of frankincense (1,400–600) (v/v) or sandalwood (16,000–7,000) (v/v) essential oil in the growth media was added to each well in triplicate. Measurement of cell viability following exposure to essential oils was performed using the XTT assay at time 0 and 24 h. During the assay, the growth medium (100 μL) was removed from each well and the XTT labeling mixture (50 μL) was added to each well. The reactions were carried out at 37°C for 4 h and the absorbance of the reaction product was recorded at 450 nm by a μQuant microplate reader (Bio-Tek; Winooski, VT, USA). Absorbance values obtained at 24 h following treatment were normalized to the values obtained at time 0 to calculate fold changes in cell viability .
RNA extraction and quality evaluation
Total RNA was isolated from J82 cells for microarray analysis to determine frankincense and sandalwood essential oil-regulated molecular responses at the mRNA level. Briefly, 5 × 105 J82 cells were seeded in 60-mm tissue culture plates, cultured overnight for adherence, and either left untreated or treated with a 1:1,100 dilution (v/v) of frankincense essential oil or a 1:11,000 (v/v) dilution of sandalwood essential oil in growth medium. Total RNA was isolated at 0 h (no treatment) and at 0.5, 1, 2, and 3 h after treatment using the RNeasy® Mini total RNA isolation kit. Total RNA concentration was determined using a nanodrop scanning spectrophotometer, and RNA was qualitatively assessed for degradation based on the ratio of 28S:18S rRNA using a capillary gel electrophoresis system (Agilent 2100 Bioanalyzer, Agilent Technologies).
Probe synthesis and hybridization
Fluorescence-labeled probes were prepared from an aliquot (250 ng) of total RNA from each time point following essential oil treatments. First strand cDNA was reverse transcribed from total RNA in T7-oligo-dT; and cRNA was synthesized by in vitro transcription from the T7 promoter and labeled with biotinylated UTP by the Illumina Total Prep RNA Amplification Kit (Ambion; Austin. TX, USA). The biotinylatedlabeled cRNA probes were hybridized overnight to Illumina human Ref-8 version 3 BeadChips containing probes representating a total of 24,526 transcripts (Illumina, San Diego, CA, USA). Following stringent washes, microarray chips were incubated with streptavidin-Cy3 (Amersham Biosciences; Piscataway, NJ, USA) and scanned on an Illumina BeadArray Reader (Illumina).
Microarray data normalization and transformation
Normalization of the microarray datasets was performed for each array by plotting a frequency histogram of the raw expression data for all genes [27, 28]. The histograms showed a right-skewed unimodal distribution curve with the mode around zero. A normal distribution curve representing the genes with low levels of expression was then fitted around the mode, mirroring the Gaussian profile of the left part of the histogram. Two parameters mean and standard deviation (SD), were defined for this distribution. The expression values of the rest of genes were then normalized to the SD of the normal distribution curve after subtraction of the mean, and log10-transformed. The arrays were then adjusted to each other by robust linear regression [29, 30]. Genes with expression value < 3.0 (or 0.477 on the log10 scale) were considered to be not expressed; this was equivalent to setting a threshold at three SD above the mean of low levels expression genes.
Identification and clustering of hypervariable (HV) genes across the time course
We defined genes with statistically significant expression levels across the time points when comparing to a reference group  as HV genes. The reference group consisted of genes expressed above background levels but with low variability as determined by an F-test, and represented the technical variability of microarray data [27, 32]. Genes whose expression levels varied significantly (P < 1/N, where N is the number of genes expressed above the background) were identified by comparing an individual gene’s variability to that of the reference group using an F-test . The threshold level at P < 1/N is a modification of the Bonferroni correction for multiple hypothesis tests. These genes were further filtered to remove those with variability arising from experimental errors by comparing the variability of the residuals in the replicated group of samples with the same variability obtained after excluding the maximum and minimum, one at a time, using a modification of the ‘leave one out’ method . A statistically significant decrease in variability after excluding one replicate provided evidence of possible error in that particular replicate. Once filtered, the remaining genes were denoted as HV genes , and were considered to be a snapshot of the dynamic biological responses following treatment with frankincense and sandalwood essential oils. The HV genes were mean-centered and clustered by hierarchical clustering using the ‘Correlation (uncentered)’ similarity metric. All statistical analyses and clustering were performed in Matlab (Natick, MA, USA).
Ontology and canonical pathway analysis
Significantly over-represented gene ontologies for the identified HV genes were analyzed using the Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov/) . The genes of interest were uploaded as Illumina Probe IDs to DAVID; and the Illumina HUMANREF-8_V3_0_R1_11282963_A array was selected as a background. Functional annotation clustering was performed with default settings and medium stringency. Statistically significant over-representation of gene ontologies was reflected by the “enrichment score”; details of its calculation can be found on the DAVID website. Gene lists were analyzed for over-represented canonical pathways using Ingenuity Pathway Analysis (IPA; Ingenuity Systems, Redwood City, CA). Genes of interest were bound into a network using the “path explorer” tool, and the network was edited using the “path designer” tool in IPA.
Semi-quantitative and quantitative reverse transcription (RT)-polymerase chain reaction (PCR) analysis
Total RNA samples were isolated from J82 cells treated with frankincense and sandalwood essential oils for RT-PCR. For semi-quantitative analysis, first-strand cDNA was synthesized from 2.5 μg of the total RNA with oligo (dT)12–18 primers (1 μg) (Life Technologies) and MMLV reverse transcriptase (200 U) (Life Technologies) in a total volume of 50 μL; the reaction was performed at 42°C for 2 h. Target mRNA species were amplified using a three-temperature cycle protocol: 94°C for 1 min, 55°C for 1 min, and 72°C for 1 min. PCR amplified products were separated on 1.5% agarose gels (Lonza, Walkersville, MD, USA) and stained with 0.5 μg/mL ethidium bromide (Life Technologies). Ethidium bromide-stained images were captured using a Gel Doc 1000 imaging system equipped with Quantity One® image analysis software (Bio-Rad Laboratories; Hercules, CA, USA).
SYBR Green-based quantitative PCR was used for quantifying the expression levels of target genes that are either commonly or differentially regulated by frankincense and sandalwood essential oils. Isolated total RNA was subjected to RNase-free DNase I (Qiagen, Valencia, CA, USA) (1 unit) digestion. PCR was performed using the above described three-temperature cycle with 30 s for each temperature reaction using a Bio-Rad iCycler real-time PCR detection system (Bio-Rad Laboratories) Melt-curve analysis was conducted after the final cycle to ensure the amplification of single species of PCR products in each sample; and agarose gel electrophoresis was used to confirm the presence of a single PCR band with the expected size in each reaction. We used β-actin as a reference gene to correct for sample-to-sample variations. Relative changes in JUN, DUSP1, PMAIP1, and ZNF311 expression after frankincense and sandalwood essential oils treatments were calculated based on Livak et al..
The results were expressed as means (SD) from at least 3 repeats. Comparisons of cell viability following frankincense oil treatment were made using one-way analysis of variance followed by a post hoc Dunnett's test. The Student’s t-test was used to compare levels of target gene expression between frankincense and sandalwood essential oil-treated J82 cells. P < 0.05 was considered statistically significant.
Identification of frankincense and sandalwood essential oil chemical components
Frankincense and sandalwood essential oils were characterized by polarimetry and GC-MS. One of the best chemical signatures for Boswellia carterii is alpha-pinene ranging between 33.9% and 62.5% with a distinctive negative optical rotation (−13.3° ± 4.9°) as we previously reported .
Sandalwood essential oil (n = 23) exhibited unique chemical compounds cis-alpha-santalol and cis-beta-santalol at high concentrations, 41.1–47.9% and 15.9–21.3%, respectively. GC-MS with nonpolar HP-1 and polar DB-WAX qualified the purity of sandalwood essential oil as the method separated these two chemical components from similar compounds . It is extremely difficult to establish standards for natural products. Even though standard operating procedures are established for harvest, storage, and processing, variable including weather, soil condition, and water quality all have significant impacts on the chemical compositions of all herbal products (i.e., essential oils in this case). Therefore, this industry and scholarly publications present the chemical compositions of essential oils using a percentage range. However, with the introduction of analytical chemistry and molecular biology, we would like to relate and translate such descriptions of chemical compositions into descriptions of molecular responses, starting with the present article. Any deviation in molecular responses from the current observations will be related to changes in chemical compositions in the future. We hope that we will be able to create a “standard” for medical applications of essential oils or for all herbal products in the future.
Tumor cell-selective versus non-selective cytotoxicity
Both frankincense and sandalwood essential oils suppressed the viability of human bladder cancer J82 cells but with different morphologies (Figure 1A). At high concentrations of both essential oils (low dilution factors), no viable cells were detected, while frankincense and sandalwood essential oils significantly suppressed J82 cell viability at 1:1,100 (P = 0.008) and 1:11,000 (P = 0.012) dilutions (v/v), respectively (Figure 1B).
Based on the results of cell viability assays, J82 cells were more sensitive than UROtsa cells to frankincense essential oil-suppressed cell viability, similar to the findings in our previous study ; the IC50 values were 1:1,250 and 1:600 dilutions (v/v) for J82 cells and UROtsa cells, respectively. In contrast, both J82 and UROtsa cells responded similarly to sandalwood essential oil treatment, with an IC50 value around 11,000 dilutions (v/v) for both cell lines (Figure 2).
Identification of HV genes in frankincense and sandalwood essential oil-treated J82 cells
Temporal changes in gene expression patterns were analyzed in total RNA samples isolated from untreated J82 cells and cells that received either frankincense or sandalwood essential oil for between 0.5 and 3 h. Full microarray data have been deposited in Gene Expression Omnibus (GEO) with the accession number GSE53171, and are accessible on the GEO web-site (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE53171).
A total of 139 genes were identified to be HV genes in frankincense essential oil-treated cells, whereas 206 HV genes were identified to be regulated by sandalwood essential oil (Additional file 1: Table S1). Among the HV genes, 51 genes were commonly regulated by both essential oils with a gradual increase in their expression levels over the experimental period (Figure 3A). Although the remaining HV genes also showed a gradual increase following stimulation with both essential oils, these genes responded differently between frankincense and sandalwood essential oils. A total of 88 HV genes were specifically regulated by frankincense essential oil, but their expression levels remained relatively constant following sandalwood essential oil treatment (Figure 3B). Another 155 HV genes were classified as sandalwood essential oil-responsive genes (Figure 3C).
Ontological comparison of essential oil-regulated genes
The 139 and 206 HV genes regulated by frankincense and sandalwood essential oils, respectively, were subjected to DAVID analysis with Gene Ontologies (GO) of the over-represented HV genes. Four groups of genes were categorized into the GO terms of transcription factor activity, histone genes, negative regulation of biological process, and apoptosis (Table 1). Transcription factor activity was the most significantly over-represented class of genes regulated in frankincense and sandalwood essential oil-treated J82 cells. The expression of some transcription factors, such as ATF3, DDIT3, EGR1, FOSB, JUN, JUNB, and MYC, was modulated by both essential oils. Other transcription factors such as FOS (Finkel–Biskis–Jinkins murine osteogenic sarcoma virus) were specifically modulated by frankincense essential oil, whereas sandalwood essential oil modulated the expression of several inhibitors of DNA binding (ID1, ID2, and ID3), along with members of the zinc finger family. In addition, all identified histone core genes, except for HIST1H3D, were commonly modulated by both essential oils.
Negative regulation of biological processes categorized by IPA was more pronounced in sandalwood essential oil-treated J82 cells. Among the most notable HV genes, growth arrest genes (GADD45A, GADD45B, and PPP1R15A) and pro-apoptotic genes (CASP9 and ING5) were specifically up-regulated by sandalwood essential oil. As part of the negative regulation of biological processes, several pro-inflammatory interleukins (IL1A, IL6, and IL8) were up-regulated by both essential oils.
Specific genes modulated by frankincense or sandalwood essential oil
The 51 genes commonly regulated by both essential oils were excluded to identify over-represented genes specifically modulated by either frankincense or sandalwood essential oil in J82 cells. The remaining 88 and 155 HV genes specifically modulated by frankincense and sandalwood essential oils, respectively (Figure 3B and C), were analyzed by DAVID. Heat shock proteins (HSP), including three members of DNAJ family (HSP40 homolog), were the main over-represented call of HV genes following frankincense essential oil treatment (Table 2). In addition to the histone core family genes (Table 1), three other members of this family were significantly over-represented among the HV genes identified in frankincense essential oil-treated J82 cells.
Transcriptional regulation was identified as the most over-represented class of genes that were specifically modulated by sandalwood essential oil in J82 cells. Among these transcription factors, members of the zinc finger family of proteins were predominantly identified (Table 2). In addition, cell death-related genes were over-represented in sandalwood essential oil-treated cells, specifically epithelial membrane protein 1, immediate early response 3, and myeloid-associated differentiation marker.
Canonical pathway comparison of genes modulated by frankincense and sandalwood essential oils
The HV genes modulated by frankincense and sandalwood essential oil treatments in J82 cells were analyzed by IPA to identify canonical pathways through connecting both essential oil-modulated genes in a single network (Figure 4). Consistently, cell death genes (67/84 genes, P = 8.27E-33) and apoptotic genes (62/84 genes, P =1.69E32) represented significantly enriched ontological categories.
Pathways including the p38 MAPK, p53, IL6, and HMGB1 signaling pathways were commonly represented by HV genes modulated by both essential oils (Table 3). IPA also identified canonical pathways that were specific for either frankincense or sandalwood essential oil treatment. For example, aryl hydrocarbon receptor (AhR) signaling and NRF2-mediated oxidative stress response genes were the two main canonical pathways over-represented by the HV genes regulated by frankincense essential oil. In contrast, ATM signaling was uniquely over-represented by the HV genes identified in J82 cells treated with sandalwood essential oil.
Analysis of transcript expression by RT-PCR
Microarray and RT-PCR analysis provided consistent results for these genes. Although both JUN and DUSP1 were up-regulated by both essential oils in J82 cells (Figure 5A), the levels of DUSP1 were significantly elevated in frankincense essential oil-treated cells (P = 0.011), whereas there was no statistical significance for JUN (P = 0.282) at 2 h following treatment (Figure 5B). DNAJB4 was up-regulated by frankincense essential oil (P = 0.043), whereas PMAIP1 (P = 0.032) and ZNF311 (P = 0.035) transcripts were specifically up-regulated in sandalwood essential oil-treated J82 cells at 2 h after stimulation.
In this study, we generated GC-MS profiles for hydrodistilled frankincense (B. carterii from Somalia) and steam-distilled sandalwood (S. album from Sri Lanka) essential oils and compared their anti-cancer activities in human bladder cancer J82 and immortalized normal human urothelial UROtsa cell lines.
A genome-wide gene expression analysis differentiated the genes and pathways that were modulated by frankincense and sandalwood essential oils in J82 cells. We compared comprehensive transcriptome expression patterns between frankincense and sandalwood essential oils based on their distinctive chemical compositions. Sandalwood essential oil contains the highest percentages of hydroxylated sesquiterpenes, such as santalols (47–56%) ; frankincense essential oil consists of high percentages of monoterpenes, such as alpha-pinene (60–80%) [7, 24].
These two essential oils had different chemical profiles [24, 39], but elicited similar cellular responses leading to activation of activator protein 1 (AP-1). AP-1 activation can be induced by environmental stresses, growth factors, and cytokines [40–42], and is implicated in cell proliferation, differentiation, apoptosis, and transformation [43, 44]. Several pro-inflammatory cytokines, including IL1A, IL6, and IL8, commonly regulated by both essential oils, have been implicated in the p38 MAPK pathway and might be responsible for subsequent AP-1 activation . However, differential expression of AP-1 components following essential oil exposure was observed. Several AP-1 subunits, FOSB, JUN, and JUNB, were regulated by both frankincense and sandalwood essential oils, while FOS expression was specifically affected by frankincense essential oil. The composition of the AP-1 heterodimer might determine the AP-1-regulated proliferative or apoptotic activities . Unique FOS regulation by frankincense essential oil might cause differential effects of AP-1 activation. FOS forms stable heterodimers with JUN and enhances AP-1 transcriptional capacity , but is not required for apoptosis . The AP-1 stability owing to FOS expression following frankincense essential oil treatment might protect normal cells against the essential oil’s cytotoxicity.
In addition to AP-1 transcription factor, EGR1 is another transcription factor that is involved in the negative regulation of biological processes . Other transcription factors, such as CEBPB, ENO3, ID2, ID3, and IRF1, specifically regulated by sandalwood essential oil, are negative regulators of biological processes, suggesting that sandalwood essential oil suppresses cancer cell viability via an array of multiple transcription factors activities. A large number of histone core genes were regulated by both frankincense and sandalwood essential oils, suggesting that both essential oils alter chromosome structure and the accessibility of DNA to transcription factors.
MYC expression was up-regulated by both frankincense and sandalwood essential oils in J82 cells. Although MYC sensitized cells to undergo apoptosis , frankincense essential oil up-regulated a MYC partner MXD1 (MAX dimerization protein 1, or MAD). MXD1 might antagonize MYC transcriptional activity by forming a DNA-binding complex with MAX in the core sequence 5'-CAC[GA]TG-3’ .
Aryl hydrocarbon signaling was specifically overrepresented in J82 cells treated with frankincense essential oil, as evidenced by the up-regulation of AhR mRNA. Activated AhR regulates downstream gene expression, including xenobiotic metabolizing enzymes and phase II metabolizing enzymes, as well as growth factors and p53 target genes such as p21CIP1 (CDKN1A). Activation of AhR signaling may be a consequence of an imbalance between the production of reactive oxygen and the detoxification of reactive intermediates in frankincense essential oil-treated J82 cells, as reflected by possible activation of the NRF2-mediated oxidative stress response pathway. Two regulators of NRF2 transcription factor, FOS and JUN, as well as NRF2-regulated stress response genes (DNAJ heat shock proteins), were specifically up-regulated by frankincense essential oil. The selective cancer cell death induced by frankincense essential oil could be a result of oxidative stress.
The effect on ataxia telangiectasia mutated (ATM) signaling was unique to sandalwood essential oil-treated J82 cells. ATM protein is a key regulator of multiple signaling cascades following damage and subsequent DNA repair. The expression of RAD50 was up-regulated in this study, and RAD50 might activate ATM. RAD50 protein is involved in DNA double-strand break repair; and its expression might be a result of DNA damage induced by sandalwood essential oil. The ATM activation was evidenced by up-regulations of multiple p53 substrates, including the GADD45 (growth arrest and DNA damage) complex and CDKN1A. Up-regulation of GADD45 and CDKN1A, which can lead to cell cycle arrest at the G2/M checkpoint with subsequent apoptosis [50, 51], was observed in sandalwood essential oil-treated J82 cells. Additionally, the involvement of p53 activation in sandalwood essential oil-treated J82 cells was evidenced by CASP6 and PMAIP1 (NOXA) expression in the apoptosis branch of the p53 pathway. The sandalwood essential oil-induced cytotoxicity might be a result of DNA instability, chromatin remodeling, and DNA double stranded breaks.
High mobility group box-1 (HMGB1) signaling, which is an integral component of oxidative stress and downstream cell survival or death, including apoptosis and autophagy [52, 53], was identified as a common canonical pathway activated by both frankincense and sandalwood essential oils. HMGB1 is involved in the assembly of nucleoprotein complexes to maintain nucleosome structure and regulate gene transcription, and is secreted by cells after stimulation with endotoxins and cytokines, including IL1, IL6, or IL8 . HMGB1 might act synergistically with oxidative stress induced by frankincense essential oil. The expression of the receptor for advanced glycation end products (RAGE) that binds to HMGB1  was specifically regulated in sandalwood essential oil-treated cells. HMGB1 has been demonstrated to induce apoptosis via the RAGE-p38 MAPK/ERK signaling pathway . Thus, sandalwood essential oil-induced HMGB1 signaling might occur synergistically with p38 MAPK signaling through the RAGE receptor and result in non-selective cell death.
Frankincense essential oil may represent a candidate on a growing list of natural compounds selectively eradicating cancer cells via oxidative stress [57–59]. To our knowledge, this is the first report suggesting that the NRF2-mediated oxidative stress response pathway might be involved in the tumor cell-specific anti-proliferative and pro-apoptotic activities of frankincense essential oil. In contrast, sandalwood essential oil appears to be a non-tumor cell-specific agent inducing DNA damage and cell cycle arrest. Further studies are required to confirm the involvement of these pathways in relevant biological systems, e.g., whether the redox status of HMGB1 induced by frankincense essential oil is involved in cancer cell-specific activity.
The effects of frankincense and sandalwood essential oils on J82 cells and UROtsa cells involved different mechanisms leading to cancer cell death. While frankincense essential oil elicited selective cancer cell death via NRF-2-mediated oxidative stress, sandalwood essential oil induced non-selective cell death via DNA damage and cell cycle arrest.
Aryl hydrocarbon receptor
Activation protein 1
Ataxia telangiectasia mutated
Database for annotation, Visualization and Integrated discovery
Epidermal growth factor
Growth arrest and DNA damage
High mobility group box-1
Heat shock protein
Immediate early response
Ingenuity pathway analysis
Transitional cell carcinoma
Finkel–Biskis–Jinkins murine osteogenic sarcoma virus.
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.
Langmead L, Rampton DS: Review article: complementary and alternative therapies for inflammatory bowel disease. Aliment Pharmacol Ther. 2006, 23: 341-349.
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.
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.
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.
Frank MB, Yang Q, Osban J, Azzarello JT, Saban MR, Saban R, Ashley RA, Welter JC, Fung KM, Lin HK: Frankincense oil derived from Boswellia carteri induces tumor cell specific cytotoxicity. BMC Complement Altern Med. 2009, 9: 6-
Suhail MM, Wu W, Cao A, Mondalek FG, Fung KM, Shih PT, Fang YT, Woolley C, Young G, Lin HK: Boswellia sacra essential oil induces tumor cell-specific apoptosis and suppresses tumor aggressiveness in cultured human breast cancer cells. BMC Complement Altern Med. 2011, 11: 129-
Ni X, Suhail MM, Yang Q, Cao A, Fung KM, Postier RG, Woolley C, Young G, Zhang J, Lin HK: Frankincense essential oil prepared from hydrodistillation of Boswellia sacra gum resins induces human pancreatic cancer cell death in cultures and in a xenograft murine model. BMC Complement Altern Med. 2012, 12: 253-
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 Med. 2002, 68: 397-401.
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 Med. 1998, 64: 328-331.
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.
Zhao W, Entschladen F, Liu H, Niggemann B, Fang Q, Zaenker KS, Han R: Boswellic acid acetate induces differentiation and apoptosis in highly metastatic melanoma and fibrosarcoma cells. Cancer Detect Prev. 2003, 27: 67-75.
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.
Pang X, Yi Z, Zhang X, Sung B, Qu W, Lian X, Aggarwal BB, Liu M: Acetyl-11-keto-β-boswellic acid inhibits prostate tumor growth by suppressing vascular endothelial growth factor receptor 2-mediated angiogenesis. Cancer Res. 2009, 69: 5893-5900.
Lu M, Xia L, Hua H, Jing Y: Acetyl-keto-β-boswellic acid induces apoptosis through a death receptor 5-mediated pathway in prostate cancer cells. Cancer Res. 2008, 68: 1180-1186.
Syrovets T, Gschwend JE, Buchele B, Laumonnier Y, Zugmaier W, Genze F, Simmet T: Inhibition of IκB kinase activity by acetyl-boswellic acids promotes apoptosis in androgen-independent PC-3 prostate cancer cells in vitro and in vivo. J Biol Chem. 2005, 280: 6170-6180.
Lin HK, Suhail MM, Fung KM, Woolley CL, Young DG: Extraction of biologically active compounds by hydrodistillation of species gum resins for anti-cancer therapy. OA Altern Med. 2013, 1: 4-
Moy RL, Levenson C, So JJ, Rock JA: Single-center, open-label study of a proprietary topical 0.5% salicylic acid-based treatment regimen containing sandalwood oil in adolescents and adults with mild to moderate acne. J Drugs Dermatol. 2012, 11: 1403-1408.
Ohmori A, Shinomiya K, Utsu Y, Tokunaga S, Hasegawa Y, Kamei C: Effect of santalol on the sleep-wake cycle in sleep-disturbed rats. Nihon Shinkei Seishin Yakurigaku Zasshi. 2007, 27: 167-171.
Misra BB, Dey S: Comparative phytochemical analysis and antibacterial efficacy of in vitro and in vivo extracts from East Indian sandalwood tree (Santalum album L.). Lett Appl Microbiol. 2012, 55: 476-486.
Dwivedi C, Abu-Ghazaleh A: Chemopreventive effects of sandalwood oil on skin papillomas in mice. Eur J Cancer Prev. 1997, 6: 399-401.
Dwivedi C, Zhang Y: Sandalwood oil prevent skin tumour development in CD1 mice. Eur J Cancer Prev. 1999, 8: 449-455.
Stichlmair J, Fair JR: Distillation : principles and practices. 1998, New York: Wiley
Woolley CL, Suhail MM, Smith BL, Boren KE, Taylor LC, Schreuder MF, Chai JK, Casabianca H, Haq S, Lin HK, Al-Shahri AA, Al-Hatmi S, Young DG: Chemical differentiation of Boswellia sacra and Boswellia carterii essential oils by gas chromatography and chiral gas chromatography–mass spectrometry. J Chromatogr A. 2012, 1261: 158-163.
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.
Petzoldt JL, Leigh IM, Duffy PG, Sexton C, Masters JR: Immortalisation of human urothelial cells. Urol Res. 1995, 23: 377-380.
Dozmorov I, Knowlton N, Tang Y, Centola M: Statistical monitoring of weak spots for improvement of normalization and ratio estimates in microarrays. BMC Bioinformatics. 2004, 5: 53-
Knowlton N, Dozmorov IM, Centola M: Microarray Data Analysis Toolbox (MDAT): for normalization, adjustment and analysis of gene expression data. Bioinformatics. 2004, 20: 3687-3690.
Hua J, Balagurunathan Y, Chen Y, Lowey J, Bittner ML, Xiong Z, Suh E, Dougherty ER: Normalization benefits microarray-based classification. EURASIP J Bioinform Syst Biol. 2006, 43056-doi:10.1155/BSB/2006/43056
Gusnanto A, Calza S, Pawitan Y: Identification of differentially expressed genes and false discovery rate in microarray studies. Curr Opin Lipidol. 2007, 18: 187-193.
Dozmorov IM, Centola M, Knowlton N, Tang Y: Mobile classification in microarray experiments. Scand J Immunol. 2005, 62 (Suppl 1): 84-91.
Dozmorov I, Centola M: An associative analysis of gene expression array data. Bioinformatics. 2003, 19: 204-211.
Varma S, Simon R: Bias in error estimation when using cross-validation for model selection. BMC Bioinformatics. 2006, 7: 91-
Dozmorov I, Knowlton N, Tang Y, Shields A, Pathipvanich P, Jarvis JN, Centola M: Hypervariable genes–experimental error or hidden dynamics. Nucleic Acids Res. 2004, 32: e147-
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: 3-
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT Method. Methods. 2001, 25: 402-408.
Howes MJ, Simmonds MS, Kite GC: Evaluation of the quality of sandalwood essential oils by gas chromatography–mass spectrometry. J Chromatogr A. 2004, 1028: 307-312.
Subasinghe U, Gamage M, Hettiarachchi DS: Essential oil content and composition of Indian sandalwood (Santalum album) in Sri Lanka. J Forest Res. 2013, 24: 127-130.
Shellie R, Marriott P, Morrison P: Comprehensive two-dimensional gas chromatography with flame ionization and time-of-flight mass spectrometry detection: qualitative and quantitative analysis of West Australian sandalwood oil. J Chromatogr Sci. 2004, 42: 417-422.
Devary Y, Gottlieb RA, Smeal T, Karin M: The mammalian ultraviolet response is triggered by activation of Src tyrosine kinases. Cell. 1992, 71: 1081-1091.
Quantin B, Breathnach R: Epidermal growth factor stimulates transcription of the c-jun proto-oncogene in rat fibroblasts. Nature. 1988, 334: 538-539.
Korber B, Mermod N, Hood L, Stroynowski I: Regulation of gene expression by interferons: control of H-2 promoter responses. Science. 1988, 239: 1302-1306.
Angel P, Karin M: The role of Jun, Fos and the AP-1 complex in cell-proliferation and transformation. Biochim Biophys Acta. 1991, 1072: 129-157.
Shaulian E, Karin M: AP-1 as a regulator of cell life and death. Nat Cell Biol. 2002, 4: E131-E136.
Eliopoulos AG, Gallagher NJ, Blake SM, Dawson CW, Young LS: Activation of the p38 mitogen-activated protein kinase pathway by Epstein-Barr virus-encoded latent membrane protein 1 coregulates interleukin-6 and interleukin-8 production. J Biol Chem. 1999, 274: 16085-16096.
Halazonetis TD, Georgopoulos K, Greenberg ME, Leder P: c-Jun dimerizes with itself and with c-Fos, forming complexes of different DNA binding affinities. Cell. 1988, 55: 917-924.
Gajate C, Alonso MT, Schimmang T, Mollinedo F: C-Fos is not essential for apoptosis. Biochem Biophys Res Commun. 1996, 218: 267-272.
Jamieson GA, Mayforth RD, Villereal ML, Sukhatme VP: Multiple intracellular pathways induce expression of a zinc-finger encoding gene (EGR1): relationship to activation of the Na/H exchanger. J Cell Physiol. 1989, 139: 262-268.
Hooker CW, Hurlin PJ: Of Myc and Mnt. J Cell Sci. 2006, 119: 208-216.
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.
Vogelstein B, Lane D, Levine AJ: Surfing the p53 network. Nature. 2000, 408: 307-310.
Tsung A, Klune JR, Zhang X, Jeyabalan G, Cao Z, Peng X, Stolz DB, Geller DA, Rosengart MR, Billiar TR: HMGB1 release induced by liver ischemia involves Toll-like receptor 4 dependent reactive oxygen species production and calcium-mediated signaling. J Exp Med. 2007, 204: 2913-2923.
Tang D, Kang R, Cheh CW, Livesey KM, Liang X, Schapiro NE, Benschop R, Sparvero LJ, Amoscato AA, Tracey KJ, Zeh HJ, Lotze MT: HMGB1 release and redox regulates autophagy and apoptosis in cancer cells. Oncogene. 2010, 29: 5299-5310.
Wang H, Bloom O, Zhang M, Vishnubhakat JM, Ombrellino M, Che J, Frazier A, Yang H, Ivanova S, Borovikova L, Manogue KR, Faist E, Abraham E, Andersson J, Andersson U, Molina PE, Abumrad NN, Sama A, Tracey KJ: HMG-1 as a late mediator of endotoxin lethality in mice. Science. 1999, 285: 248-251.
Hori O, Brett J, Slattery T, Cao R, Zhang J, Chen JX, Nagashima M, Lundh ER, Vijay S, Nitecki D, Morser J, Stern D, Schmidt AM: The receptor for advanced glycation end products (RAGE) is a cellular binding site for amphoterin. Mediation of neurite outgrowth and co-expression of rage and amphoterin in the developing nervous system. J Biol Chem. 1995, 270: 25752-25761.
Kim SW, Lim CM, Kim JB, Shin JH, Lee S, Lee M, Lee JK: Extracellular HMGB1 released by NMDA treatment confers neuronal apoptosis via RAGE-p38 MAPK/ERK signaling pathway. Neurotox Res. 2011, 20: 159-169.
Murias M, Luczak MW, Niepsuj A, Krajka-Kuzniak V, Zielinska-Przyjemska M, Jagodzinski PP, Jager W, Szekeres T, Jodynis-Liebert J: Cytotoxic activity of 3,3’,4,4’,5,5’-hexahydroxystilbene against breast cancer cells is mediated by induction of p53 and downregulation of mitochondrial superoxide dismutase. Toxicol In Vitro. 2008, 22: 1361-1370.
Verrax J, Taper H, Buc Calderon P: Targeting cancer cells by an oxidant-based therapy. Curr Mol Pharmacol. 2008, 1: 80-92.
Verrax J, Pedrosa RC, Beck R, Dejeans N, Taper H, Calderon PB: In situ modulation of oxidative stress: a novel and efficient strategy to kill cancer cells. Curr Med Chem. 2009, 16: 1821-1830.
This work was supported in part by NIH grant 5P20RR020143-06 to JDW and MD.
CW and GY are employed by the company Young Living Essential Oils. Other authors declare that they have no competing interests.
MS, CW, GY, KMF, and HKL conceived the study and designed the experiments. QY, WW, HYL, MD and JW performed the experiments. MS, CW, GY, KMF, HKL, MD and JW interpreted the experimental results. All authors prepared the manuscript preparations and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Dozmorov, M.G., Yang, Q., Wu, W. et al. Differential effects of selective frankincense (Ru Xiang) essential oil versus non-selective sandalwood (Tan Xiang) essential oil on cultured bladder cancer cells: a microarray and bioinformatics study. Chin Med 9, 18 (2014). https://doi.org/10.1186/1749-8546-9-18