Diagnosis, Therapy, and Prognosis for Hepatocellular CarcinomaView this Special Issue
Research Article | Open Access
Gene Ontology and Expression Studies of Strigolactone Analogues on a Hepatocellular Carcinoma Cell Line
Human hepatocellular carcinoma (HCC) is the most common and recurrent type of primary adult liver cancer without any effective therapy. Plant-derived compounds acting as anticancer agents can induce apoptosis by targeting several signaling pathways. Strigolactone (SL) is a novel class of phytohormone, whose analogues have been reported to possess anticancer properties on a panel of human cancer cell lines through inducing cell cycle arrest, destabilizing microtubular integrity, reducing damaged in the DNA repair machinery, and inducing apoptosis. In our previous study, we reported that a novel SL analogue, TIT3, reduces HepG2 cell proliferation, inhibits cell migration, and induces apoptosis. To decipher the mechanisms of TIT3-induced anticancer activity in HepG2, we performed RNA sequencing and the differential expression of genes was analyzed using different tools. RNA-Seq data showed that the genes responsible for microtubule organization such as TUBB, BUB1B, TUBG2, TUBGCP6, TPX2, and MAP7 were significantly downregulated. Several epigenetic modulators such as UHRF1, HDAC7, and DNMT1 were also considerably downregulated, and this effect was associated with significant upregulation of various proapoptotic genes including CASP3, TNF-α, CASP7, and CDKN1A (p21). Likewise, damaged DNA repair genes such as RAD51, RAD52, and DDB2 were also significantly downregulated. This study indicates that TIT3-induced antiproliferative and proapoptotic activities on HCC cells could involve several signaling pathways. Our results suggest that TIT3 might be a promising drug to treat HCC.
In 2012, 0.8 million patients were diagnosed with liver cancer, the seventh highest age-related incidence rate globally . Human hepatocellular carcinoma (HCC) is the most common type of primary liver cancer in adults, as well as the most frequently recurrent malignancy without any effective therapy [2, 3]. HCC is the third most common cancer-related cause of mortality globally . Several etiological factors could lead to the development of HCC including hepatitis B virus (HBV), hepatitis C virus (HCV), alcohol, cirrhosis, and nonalcoholic fatty liver disease (NAFLD) .
There exists an inveterate history of compounds, derived from plants, serving as anticancer agents . These compounds can exert their inhibitory effects on cancer cells by targeting several pathways including cell cycle arrest, cell proliferation, and apoptosis. Strigolactones (SLs) are a novel class of phytohormones, which control the branching of shoot architecture by hindering growth and self-renewal of axillary meristem cells [5, 6]. It has previously been reported that synthetic SL analogues instigate G2/M cell cycle arrest and apoptosis by regulating the p38 and JNK1/2 MAPKs signaling pathways, causing the induction of stress in an array of solid and nonsolid human cancer cells, including prostate, colon, leukemia, osteosarcoma, and lung cancer cell lines.
It has only slight effects on the growth, survival, and viability of nontransformed human fibroblasts, healthy primary prostate cells, and mammary epithelial cells [7, 8]. SL analogues demonstrated their anticancer effects in a xenograft model of breast cancer  and have also been shown to affect the integrity of the microtubule network by impeding the migration of highly invasive breast cancer cell lines .
Recently, SL analogues have been found to destabilize the genomic DNA of cancer cells by inducing DNA double-strand breaks (DSBs) and activating the DNA damage and simultaneously hindering DNA repair, resulting in cell death. It is noteworthy that these activities of SL analogues have not been reported in nontransformed BJ fibroblast cells . Furthermore, the efficiency of the delivery of SL analogues (hence their therapeutic efficacy) to prostate cancer cells may get enhanced by encapsulation of SL analogues in glutathione/pH-responsive nanosponges .
Synthetic SL analogues have been reported to downregulate RAD51 expression through ubiquitination in a proteasome-dependent way, hence, reducing the localization of RAD51 to DSB sites . RAD51 is a crucial component of a prominent DNA repair pathway, the homology-directed repair (HDR) machinery. Overexpression of RAD51 causes an increase in DNA repair activity, which could result in resistance to DNA damage that is usually imposed by radiotherapy or chemotherapy [12, 13].
In our previous study, we found that the newly synthesized SL analogue, TIT3, inhibits proliferation and induces apoptosis of HepG2 (hepatocellular carcinoma) cells with minimal effects on healthy cells . The molecular structure of TIT3 is shown in Figure 1.
In the present study, we performed RNA sequencing and the differential expression of genes was analyzed using different tools to analyze and investigate the differential gene expression of HepG2 cells treated with TIT3 and to disclose the possible signaling pathways leading to the inhibition of cell proliferation and induction of apoptosis.
2. Materials and Methods
2.1. Cell Culture and Treatment
HepG2 cells were obtained from ATCC (Manassas, Virginia, USA). These cells were then sustained at 37°C, in a humidified incubator, at 5% CO2. DMEM (UFC Biotech, Riyadh, KSA) supplemented with 10% fetal bovine serum (FBS) (LifeTech, catalogue no: 16000-044) and 1% (100 U/ml) penicillin-streptomycin (LifeTech, catalogue no: 15140-122) was used to maintain cells.
2.2. RNA-Seq, Differentially Expressed Genes, and Bioinformatics Analysis
HepG2 cells were treated with 60 μM of SL analogue TIT3 (IC50 of 63.46 μM)  for 24 hours, in triplicate. An RNeasy kit (QIAGEN) was used to extract the total RNA, and the concentration of RNA was quantified. A bioanalyzer was used to analyze the quality of the total RNA (RIN ). Sequencing libraries were then generated using TruSeq Stranded mRNA Sample Preparation Kits (Illumina, CA) from 2500 ng of the total RNA from each of the three replicates. The Illumina HiSeq 2000 system was used to conduct 50 bp long single-end deep sequencing. The FASTX-Toolkit was used to remove the adaptor sequence and filtering of low-quality base call and low-quality reads. The short filtered sequencing reads that were acquired were mapped to the human genome by the TopHat2 and Subreads package; the featureCounts function was used to quantify the gene expression values. These gene expression values were then used to calculate the size of the library and dataset dispersion for the analysis of differentially expressed genes . Differential gene expression was examined using the R/Bioconductor package edgeR and established by log fold change (LogFC) and false discovery rate (FDR) ( or ≤-1; ).
2.3. Bioinformatics Analysis
The gene set functional analysis and pathway analysis were analyzed using the gene ontology (GO) and KEGG pathway. The gene IDs of interest were converted to EntrezID and uploaded to DAVID bioinformatics tools. GO and KEGG pathway analysis were performed by setting all the GO terms and KEGG pathway genes as background genes. Overrepresented GO terms or pathways are determined by the enrichment score (, gene ).
3.1. Gene Expression Is Regulated by TIT3
Data obtained from HepG2 cells treated with 60 μM of SL analogue of TIT3 revealed that the mRNA expression of 3240 genes was modulated, with 1473 genes being upregulated (log fold ; ) and 1767 genes being downregulated (log fold ; ). The number of altered transcripts has been organized based on the log fold change (LogFC) or the value (Tables 1 and 2). Overall, the number of transcripts being upregulated was fewer than the number of transcripts being downregulated.
3.2. Gene Enrichment Analysis of Altered Transcripts
The gene enrichment analysis of gene ontology (GO) terms () revealed that there was a significant increase in the negative regulation of transcription by the RNA polymerase II promoter and negative regulation of G1/S transition of mitosis and a substantial decrease in the damaged DNA repair genes. A summary of GO analysis with different biological processes, cell components, and molecular functions of upregulated and downregulated transcripts in HepG2 cells treated with TIT3 is shown in Figures 2 and 3, respectively.
3.3. KEGG Pathway Analysis
The KEGG pathway analysis revealed the probability of the involvement of apoptosis pathways involving TNF and PI3K-Akt (Figures 4–6). There was a significant decrease in the genes involved in the organization of microtubules such as BUB1B, TUBB, TUBG2, TUBGCP6, TPX2, and MAP7 (; ). A significant decrease was also found in the expression levels of crucial epigenetic players UHRF1, DNMT1, and HDAC7, known to inhibit the expression of several tumor suppressor genes in cancer (; ).
Additionally, genes responsible for DNA damage repair including DDB2, RAD51, and RAD52 were substantially downregulated (Table 3). The expression levels of several tumor suppressor genes such as NKX-3, FLCN, ING1, SIK1, and TP53INP1 (; ), and genes exhibiting proapoptotic activities such as TNF-α, LTA (TNF-β), CASP3, MOAP1, and CASP7 (; ) as well as genes having antiproliferative effects such as CDKN1A (p21) (; ) were significantly increased in response to TIT3 treatment (Table 4). Figure 7 shows the probable interactions of genes of the different transcriptional regulators, and Figure 8 depicts a heat map, representing the comprehensive regulation of gene expressions in terms of their values and LogFC.
Fold change treated vs control.
Fold change treated vs control.
Synthetic SL analogues have been reported to induce cell cycle arrest and apoptosis in both solid and hematological tumors by targeting several signaling pathways [7, 8]. Our previous study showed that synthetic SL analogue TIT3 inhibited the proliferation and migration and induced apoptosis of HepG2 cells with minimal toxicity towards healthy noncancerous cells . Since TIT3 impeded the migration of HepG2 cells, we suggested that such an effect is a result of the interference with the organization of the microtubular network. Data obtained from RNA-Seq showed significant downregulation of BUB1B, TUBB, TUBG2, TUBGCP6, TPX2, and MAP7 genes, which are known to be involved in the microtubular organization, suggesting that the antimetastatic and proapoptotic effects of TIT3 could be challenged by mechanisms involved in the organization of microtubules. Our results are consistent with several studies showing that SL analogues can induce apoptosis in breast cancer cells . Other cancer cell lines such as melanoma, colon, lung, prostate, and osteosarcoma were also reported to be affected by SL analogues through targeting of the microtubular network .
Furthermore, our results showed that the histone deacetylase 7 (HDAC7) was downregulated in response to TIT3 treatment. Interestingly, the unusual activity of the histone deacetylases including HDAC7 has been reported in many types of cancers . HDAC inhibitors such as TSA, SAHA, and MS-275 are useful in the chemotherapeutic regimens of many cancers including HCC and significantly inhibit cell proliferation, and migration/invasion induces cell cycle arrest and apoptosis of HCC . Knockdown of HDAC7 led to G1/S arrest in different cancer cells through the upregulation of the cell cycle inhibitor CDKN1A (p21). Interestingly, an increase in the expression levels of p21 mRNA has also been observed in TIT3-treated HepG2 cancer cell lines . All the evidences above suggest that TIT3 could act as a HDAC inhibitor causing cell cycle arrest through a p21-dependent mechanism inducing inhibition of HepG2 cell proliferation.
Moreover, UHRF1, a well-documented regulator of gene expression in cancer [19, 20], was downregulated by TIT3 suggesting that UHRF1 could be a potent target for TIT3 in HCC. In alignment with our results, UHRF1 overexpression has been demonstrated to cause tumorigenesis in different cancer types including HCC . UHRF1 inhibition by using pharmacological compounds is associated with the reactivation of various tumor suppressor genes, thus suppressing the proliferation of cancer cells by inducing apoptosis .
Double-strand breaks (DSBs) are the most notable form of DNA damage, and once the DSBs are formed, cells may undergo either of the two repair mechanisms: nonhomologous end joining (NHEJ) or homology-directed repair (HDR) [23, 24]. Previous studies have revealed that cancer cells which lack HDR are quite sensitive to DNA-damaging agents . Our data obtained from RNA-Seq showed that TIT3 induced the downregulation of damaged DNA repair genes including DDB2, RAD51, and RAD52. Therefore, TIT3 could be an inhibitor of DNA repair proteins. In support of our results, the evidence is available in the literature that SL analogues can hamper HDR and impair DSB repair .
Caspase 3 (CASP 3) is known as an executioner caspase in apoptosis resulting in the inhibition of proliferation of HepG2 cancer cells . TNF-α is also known to modulate proliferation, differentiation, and apoptosis or necrotic cell death in several different cell types including HepG2 cancer cells [27, 28]. Our results illustrated that the primary genes responsible for apoptosis including CASP 3, CASP 7, TNF-α, TNF-β, and MOAP1 were significantly upregulated by TIT3 treatment on HepG2 cancer cells.
Overall, we propose that the inhibition of HepG2 cancer cell growth was due to an interplay of genes wherein the treatment of TIT3 significantly altered their expression levels. Altered gene expressions affected cell proliferation, cell cycle, metastasis, and apoptosis. TIT3 could also be an inhibitor of HDAC and can target the organization of the microtubular network as well as affect the genes involved in DNA repair.
We provided evidence that TIT3 targets several critical pathways in HepG2 cells. Therefore, to obtain and establish a deeper understanding of the molecular mechanisms exerted by TIT3, molecular biology techniques such as Western blotting, qPCR, microarray, and proteomics must be done to reveal the specific targets.
The data used to support the findings of this study are available from the corresponding author upon request.
The data and material presented in this manuscript neither have been published before nor have been submitted for publication to another scientific journal or are being considered for publication elsewhere.
Conflicts of Interest
The authors have no relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript. This includes employment, consultancies, honoraria, stock ownership, grants, or patents received or pending and royalties.
MNH, SR, and MAH conducted experiments under the supervision of MA, TA, and HC. JR, WW, and SM analyzed the data. MA, MH, SR, CB, and MZ wrote most of the manuscript. SSM, AA, TK, and KA applied for funding and helped in the data analysis. All authors read and approved the final manuscript.
The authors acknowledge the technical support from the Deanship of Scientific Research, King Abdulaziz University, Jeddah, and KFMRC for providing the RNA-Seq facility. This work was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under grant no. RG-43-130-40. The authors, therefore, acknowledge with thanks DSR for the technical and financial support.
- Y. A. Ghouri, I. Mian, and J. H. Rowe, “Review of hepatocellular carcinoma: epidemiology, etiology, and carcinogenesis,” Journal of Carcinogenesis, vol. 16, no. 1, p. 1, 2017.
- C.-J. Liu, J.-H. Yang, F.-Z. Huang et al., “The role of MiR-99b in mediating hepatocellular carcinoma invasion and migration,” European Review for Medical and Pharmacological Sciences, vol. 22, no. 8, pp. 2273–2281, 2018.
- V. Zahir, M. N. Niki, P. Hajian, F. H. Shirazi, and H. Mirzaei, “Evaluation of silibinin effects on the viability of HepG2 (human hepatocellular liver carcinoma) and HUVEC (human umbilical vein endothelial) cell lines,” Iranian Journal of Pharmaceutical Research, vol. 17, no. 1, pp. 261–267, 2018.
- D. J. Newman and G. M. Cragg, “Advanced preclinical and clinical trials of natural products and related compounds from marine sources,” Current Medicinal Chemistry, vol. 11, no. 13, pp. 1693–1713, 2004.
- V. Gomez-Roldan, S. Fermas, P. B. Brewer et al., “Strigolactone inhibition of shoot branching,” Nature, vol. 455, no. 7210, pp. 189–194, 2008.
- M. Umehara, A. Hanada, S. Yoshida et al., “Inhibition of shoot branching by new terpenoid plant hormones,” Nature, vol. 455, no. 7210, pp. 195–200, 2008.
- C. B. Pollock, H. Koltai, Y. Kapulnik, C. Prandi, and R. I. Yarden, “Strigolactones: a novel class of phytohormones that inhibit the growth and survival of breast cancer cells and breast cancer stem-like enriched mammosphere cells,” Breast Cancer Research and Treatment, vol. 134, no. 3, pp. 1041–1055, 2012.
- C. B. Pollock, S. McDonough, V. S. Wang et al., “Strigolactone analogues induce apoptosis through activation of P38 and the stress response pathway in cancer cell lines and in conditionally reprogrammed primary prostate cancer cells,” Oncotarget, vol. 5, no. 6, pp. 1683–1698, 2014.
- E. Mayzlish-Gati, D. Laufer, C. F. Grivas et al., “Strigolactone analogs act as new anti-cancer agents in inhibition of breast cancer in xenograft model,” Cancer Biology and Therapy, vol. 16, no. 11, pp. 1682–1688, 2015.
- M. P. Croglio, J. M. Haake, C. P. Ryan et al., “Analogs of the novel phytohormone, strigolactone, trigger apoptosis and synergize with PARP inhibitors by inducing DNA damage and inhibiting DNA repair,” Oncotarget, vol. 7, no. 12, pp. 13984–14001, 2016.
- M. Argenziano, C. Lombardi, B. Ferrara et al., “Glutathione/PH-responsive nanosponges enhance strigolactone delivery to prostate cancer cells,” Oncotarget, vol. 9, no. 88, pp. 35813–35829, 2018.
- D. Schild and C. Wiese, “Overexpression of RAD51 suppresses recombination defects: a possible mechanism to reverse genomic instability,” Nucleic Acids Research, vol. 38, no. 4, pp. 1061–1070, 2010.
- A. Ward, K. K. Khanna, and A. P. Wiegmans, “Targeting homologous recombination, new pre-clinical and clinical therapeutic combinations inhibiting RAD51,” Cancer Treatment Reviews, vol. 41, no. 1, pp. 35–45, 2015.
- M. N. Hasan, H. Choudhry, S. S. Razvi et al., “Synthetic strigolactone analogues reveal anti-cancer activities on hepatocellular carcinoma cells,” Bioorganic and Medicinal Chemistry Letters, vol. 28, no. 6, pp. 1077–1083, 2018.
- Y. Lu, N. Starkey, W. Lei et al., “Inhibition of hedgehog-signaling driven genes in prostate cancer cells by Sutherlandia frutescens extract,” PLoS One, vol. 10, no. 12, article e0145507, 2015.
- W. Weichert, “HDAC expression and clinical prognosis in human malignancies,” Cancer Letters, vol. 280, no. 2, pp. 168–176, 2009.
- S.-O. Kim, B.-T. Choi, I.-W. Choi et al., “Anti-invasive activity of histone deacetylase inhibitors via the induction of Egr-1 and the modulation of tight junction-related proteins in human hepatocarcinoma cells,” BMB Reports, vol. 42, no. 10, pp. 655–660, 2009.
- K. Qu, T. Lin, J. Wei et al., “Cisplatin induces cell cycle arrest and senescence via upregulating P53 and P21 expression in HepG2 cells,” Journal of Southern Medical University, vol. 33, no. 9, pp. 1253–1259, 2013.
- P. Bashtrykov, G. Jankevicius, R. Z. Jurkowska, S. Ragozin, and A. Jeltsch, “The UHRF1 protein stimulates the activity and specificity of the maintenance DNA methyltransferase DNMT1 by an allosteric mechanism,” The Journal of Biological Chemistry, vol. 289, no. 7, pp. 4106–4115, 2014.
- C. Bianchi and R. Zangi, “UHRF1 discriminates against binding to fully-methylated CpG-sites by steric repulsion,” Biophysical Chemistry, vol. 171, pp. 38–45, 2013.
- R. Mudbhary, Y. Hoshida, Y. Chernyavskaya et al., “UHRF1 overexpression drives DNA hypomethylation and hepatocellular carcinoma,” Cancer Cell, vol. 25, no. 2, pp. 196–209, 2014.
- M. Alhosin, T. Sharif, M. Mousli et al., “Down-regulation of UHRF1, associated with re-expression of tumor suppressor genes, is a common feature of natural compounds exhibiting anti-Cancer properties,” Journal of Experimental & Clinical Cancer Research, vol. 30, no. 1, p. 41, 2011.
- F. Liang, M. Han, P. J. Romanienko, and M. Jasin, “Homology-directed repair is a major double-strand break repair pathway in mammalian cells,” Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 9, pp. 5172–5177, 1998.
- M. Zimmermann and T. de Lange, “53BP1: pro choice in DNA repair,” Trends in Cell Biology, vol. 24, no. 2, pp. 108–117, 2014.
- C. J. Lord and A. Ashworth, “The DNA damage response and cancer therapy,” Nature, vol. 481, no. 7381, pp. 287–294, 2012.
- Y. Kondoh, T. Kawada, and R. Urade, “Activation of caspase 3 in HepG2 cells by elaidic acid (T18:1),” Biochimica et Biophysica Acta, vol. 1771, no. 4, pp. 500–505, 2007.
- J. Chen, X. Sun, T. Xia, Q. Mao, and L. Zhong, “Pretreatment with dihydroquercetin, a dietary flavonoid, protected against concanavalin A-induced immunological hepatic injury in mice and TNF-α/ActD-induced apoptosis in HepG2 cells,” Food & Function, vol. 9, no. 4, pp. 2341–2352, 2018.
- L. A. Tartaglia and D. V. Goeddel, “Two TNF receptors,” Immunology Today, vol. 13, no. 5, pp. 151–153, 1992.
Copyright © 2019 Mohammed Nihal Hasan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.