Selecting molecular therapeutic drug targets based on the expression profiles of intrahepatic cholangiocarcinomas and miRNA-mRNA regulatory networks

  • Authors:
    • Boshi Sun
    • Changming Xie
    • Tongsen Zheng
    • Dalong Yin
    • Jiabei Wang
    • Yingjian Liang
    • Yuejin Li
    • Guangchao Yang
    • Huawen Shi
    • Tiemin Pei
    • Jihua Han
    • Lianxin Liu
  • View Affiliations

  • Published online on: October 15, 2015     https://doi.org/10.3892/or.2015.4330
  • Pages: 382-390
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

The incidence of intrahepatic cholangiocarcinoma (ICC) is increasing yearly, making it the second most common carcinoma after hepatocellular carcinoma among primary malignant liver tumors. Integrated miRNA and mRNA analysis is becoming more frequently used in antitumor ICC treatment. However, this approach generates vast amounts of data, which leads to difficulties performing comprehensive analyses to identify specific therapeutic drug targets. In this study, we provide an in-depth analysis of ICC function, identifying potential highly potent antitumor drugs for antitumor therapy. Two sets of whole genome expression profiles were obtained from the Gene Expression Omnibus (GEO) database. Using modular bioinformatic analysis, six core functional modules were identified for ICC. Based on a Fisher's test of the Cmap small molecule drug database, 65 drug components were identified that regulated the genes of these six core modules. Literature mining was then used to identify 15 new potential antitumor drugs.

Introduction

Intrahepatic cholangiocarcinoma (ICC) is a malignant tumor that originates from the epithelial cells of the secondary bile duct and higher order intrahepatic bile duct branches (1). Among primary malignant liver tumors, the incidence of ICC is second to only hepatocellular carcinoma, with increasing incidence yearly (2,3). Because the pathological mechanism and tumor characteristics of ICC differ from those of hepatocellular carcinoma and extrahepatic cholangiocarcinoma, the American Joint Committee on Cancer (AJCC) has formally listed ICC as an independent hepatobiliary system tumor for staging and analysis in the 7th edition of its TNM staging system published in 2010 (4). Although ICC is a rare malignant tumor (accounting for 3% of malignant digestive system tumors and 15% of malignant liver tumors), the global incidence of ICC has increased from 0.32/100,000 to 0.85/100,000 over the past three decades (an increase of 165%). At present, ICC risk factors include primary sclerosing cholangitis, congenital structural abnormalities of the bile duct, inflammatory bowel disease, intrahepatic bile duct stones and parasitic infection. In a recent study, 8% of primary sclerosing cholangitis patients developed cholangiocarcinoma during a 5-year follow-up period. For patients with choledochal cysts or Caroli's disease and other congenital anatomical abnormalities of the biliary system, the risk of malignant transformation in patients 20 years of age or older is 15% (5). In addition, chronic viral hepatitis infection is a high risk factor for ICC. A clinical report on 317 ICC patients by Zhou et al (6) showed that the HBV infection rate of ICC patients is significantly higher than that of the healthy population. However, although many ICC risk factors have been identified, most ICC patients are not exposed to these risk factors, which is an issue for early diagnosis and treatment of ICC.

Currently, the treatment for ICC involves surgical excision as the major measure. For ICC patients who have undergone surgical excision for early to middle stage tumors, the post-operation 3-year survival rate is only 40–50% (7). For patients with middle to late stage tumors who have lost the opportunity for surgery, the average life expectancy without treatment is only 5–8 months (8). Other treatment measures include local ablation, hepatic artery embolization, systemic chemotherapy and targeted therapy. Unfortunately, ICC has a low sensitivity to chemotherapeutic drugs; therefore, there is at present no generalized chemotherapy plan. However, there is preliminary evidence that cisplatin, used alone or in combination with gemcitabine and fluoropyrimidine (S-1), can extend the life expectancy of patients who cannot have surgery. There is currently no specific targeted therapeutic drug for ICC; however, VEGR, EGFR, RAF kinase and multiple receptors involved in the onset of ICC are emerging as targets.

In recent years, modular bioinformatics analysis has been frequently applied in disease and drug studies. This approach not only elucidates the characteristics of gene regulatory networks but also comprehensively and systematically analyzes signaling pathways. Here, we used a modular bioinformatics approach to analyze the expression of ICC genes and determine their roles in ICC progression. We then used a small molecule drug database to identify drugs involved in the regulation of ICC core functional groups. We identified many drugs that regulated the core functional groups that were not previously identified for the treatment of ICC or other cancers. This approach represents a new drug discovery concept for cancer research.

Materials and methods

Screening for differentially expressed mRNAs and miRNAs in cholangiocarcinoma

We analyzed two sets of expression data for cholangiocarcinoma. These data included disease and control samples, with paired mRNAs and miRNAs. Data were acquired from the Gene Expression Omnibus (GEO) database. The series accession numbers for the mRNA and miRNA expression profiles were GSE32879 and GSE32957, respectively (Table I). We organized the data in matrix format to maintain the cholangiocarcinoma samples and normal samples from the two datasets for subsequent analysis. Of these samples, 16 were cholangiocarcinomas and 5 were normal controls. To increase the reliability of target pair prediction and reduce the false-positive rate, we used miRNA-target pairs that were supported by clip-seq data from the Starbase database. To obtain additional potential targets, we selected target pairs predicted by at least one of the five prediction databases (a total of 423,976 target pairs). We pretreated the mRNA and miRNA expression profiles with the k-nearest neighbors algorithm to fill in missing values (define, k=5). We then statistically tested the fold-change values by t-test to screen differentially expressed mRNAs and miRNAs. Fold-change values >1.5 or <1/1.5 (FDR value <0.5) were used as thresholds (Fig. 1).

Table I

Gene expression for intrahepatic cholangiocarcinoma from the GEO database.

Table I

Gene expression for intrahepatic cholangiocarcinoma from the GEO database.

SeriesArray platformSamplesExperimental designGroup
GSE32879GPL6244 Affymetrix Human Gene 1.0 ST ArrayGSM13712Gene expression profiling of 16 ICC and 7 NT tissues were perfomed.Tumor
GSm13713Tumor
GSM13714Tumor
GSM13715Tumor
GSM13716Tumor
GSM13717Tumor
GSM13718Tumor
GSM13719Tumor
GSM13720Tumor
GSM13721Tumor
GSM13722Tumor
GSM13723Tumor
GSM13724Tumor
GSM13725Tumor
GSM13726Tumor
GSM13727Tumor
GSM13728Tumor
GSM13742Control
GSM13743Control
GSM13744Control
GSM13745Control
GSM13746Control
GSM13747Control
GSM13748Control
GSE32957GPL14732 Nanostring nCounter Human microRNA Expression PlatformGSM15966Gene expression profiling of 16 ICC and 5 NT liver tissuesTumor
GSm15967Tumor
GSM15968Tumor
GSM15969Tumor
GSM15970Tumor
GSM15971Tumor
GSM15972Tumor
GSM15973Tumor
GSM15974Tumor
GSM15975Tumor
GSM15976Tumor
GSM15977Tumor
GSM15978Tumor
GSM15979Tumor
GSM15980Tumor
GSM15981Tumor
GSM15996Control
GSM15997Control
GSM15998Control
GSM15999Control
GSM16000Control

[i] ICC, intrahepatic cholangiocarcinoma; NT, non-tumor.

Construction and analysis of the miRNA regulatory networks

The miRNA-mRNA target pairs were required to meet two conditions: opposite differential expression and a predicted target regulation relation. Differential expression was mapped to the network and interacting genes were labeled. Online target gene prediction software was used to identify miRNAs that were differentially expressed relative to the predicted target genes. The predicted target genes were compared with downregulated (or upregulated) mRNAs from the whole-genome microarray data (Fig. 2A). The network conforms to the in-degree distribution (Fig. 2B). Target genes without a signal were eliminated. SPSS 15.0 statistical analysis software was used to analyze the expression values of the miRNA and the corresponding target genes in the normal control group and the tumor group. Pearson's correlation analysis, with P<0.05, was used as a determination standard to screen for target genes that were negatively correlated with miRNA expression. The screening results from each group were integrated into a regulatory network of differentially expressed miRNAs and their ICC target genes. Cytoscape 2.8.2 was used for integrative network analysis of miRNAs and target genes.

Generation of differentially expressed gene modules and analysis of module function

The differentially expressed genes were mapped to the network. Pairs in which both interacting genes were differentially expressed formed the sub-network. The MCL method was used to identify functional modules (Fig. 3). To functionally analyze the miRNAs and mRNAs that were differentially expressed in ICC, the upregulated and downregulated target genes regulated by differentially expressed miRNAs from each group were run through DAVID Bioinformatics Resources version 6.7 (http://david.Abcc.ncifcrf.gov/) for functional analysis of modules in the network.

Drug screening

The connectivity map (CMap) database contains the genome-wide transcriptional expression profile for human cells exposed to active small molecules. The CMap database includes 6,100 groups of small molecule intervention experiments (small molecule intervention group and normal control group), providing a total of 7,056 expression profiles. Of these profiles, 1,309 small molecules are represented. With these 1,309 small molecules from the CMap database, we performed differential gene analysis (drug administration vs. no drug) to obtain a gene set for each drug and its differential expression pattern (that is, how the genes were affected by each small molecule drug). We then compared the genes in each module with the differential genes corresponding to the abovementioned small molecules to determine the correspondence between the genes in each module and the small molecules. A Fisher's test was used to evaluate modular genes corresponding to the small molecules. Small molecules with p-values in the top 10 were considered the most valuable.

Literature review and drug components that connect regulatory network modules

To identify novel antitumor drugs, we analyzed 65 drug components related to the modules. Perl was used to write a program for literature mining. ActivePerl 5.16.2 was then used to mine literature information from Pubmed (NCBI). The mining scope was set to include titles and abstracts, along with drug names, modular gene names and 'cancer' as key words. This search provided articles related to each drug component and cancer. Integrating the drug component screening results, literature mining results and the modular genes, we identified the drugs with the greatest potential antitumor activity for the modular core function groups.

Results

Gene module screening

From comparison of the two sets of ICC and normal tissue data, the fold-change values were used to identify differentially expressed mRNA and miRNAs (detailed information in the Materials and methods). As shown in Table I, we identified 1,925 and 191 upregulated mRNAs and miRNAs, respectively. There were also 1,062 and 28 downregulated mRNAs and miRNAs, respectively. In addition, we performed cluster analysis on differentially expressed genes and miRNAs.

Screening of miRNA and mRNA regulatory networks and ICC modules

Although the mechanism of target gene mRNA degradation or translation inhibition caused by miRNAs is not completely clear, the complimentary pairing of a miRNA seed sequence with the 3′UTR sequence of the target gene is generally considered to be the main condition for a miRNA to act on a target gene. Accordingly, miRNA expression should be opposite to mRNA expression. We selected the upregulated (or downregulated) miRNAs from the miRNA microarray for each group. Targetscan 5.1 was used to predict the target genes from the differentially expressed miRNAs for each group. The target genes were then compared with the downregulated (or upregulated) mRNAs from the mRNA microarray to screen for target genes with differentially expressed miRNAs in each group. Subsequently, the differentially expressed miRNAs from each group and the screening results of their target genes were integrated to establish the regulatory network (Fig. 2A) of differentially expressed miRNAs for ICC and their target genes. The MCL method was used to explore the modules in the network. A total of six modules were identified (termed modules 1–6) (Fig. 3).

Functional analysis of gene modules

We analyzed the six modules in detail using Go_BP, with the arrangement based on the Go score [−log (t-test p-values)]. Each of the p-values of the 10 highest-scored Go analyses was <0.01. These modules were considered the most valuable and were used for subsequent analysis. Go analysis revealed the following: module 1 responded to hypoxia, module 2 responded to cell fraction, module 3 responded to organ development, module 4 responded to glycerophospholipid metabolism, module 5 responded to cytoplasm and module 6 responded to cytoplasm and cellular ketone metabolic process (Fig. 4).

Drug component regulation of gene modules

Based on the network data, we calculated that there were 1,309 small molecules from the different gene expression experiments; however, the final results showed that a total of 1,108 small molecules induced the expression of the different genes. For the six modules identified by bioinformatic analysis, we compared the genes in each module with the differentially expressed genes that corresponded with the small molecules obtained from the CMap database. We determined the correspondence between the genes in each module and the small molecules. Using a Fisher's exact test on the raw data, we calculated the p-values of each drug component for each module. The p-values were arranged in descending order, with the first 10 drug components with the smallest p-values considered to be the most valuable. These drug components were selected for subsequent analysis. A total of 65 drug components regulated the 6 modules (Fig. 5).

Discovery of novel potential anticancer drugs

To identify new anticancer drugs, we conducted a literature review and a retrospective study. The drug data were all from recent anticancer therapy studies. Interestingly, anticancer drug components that had not been well-studied showed potential as potent anticancer drugs (Fig. 6). Our review also identified many drug components that have been widely used in anticancer treatment, such as daunorubicin, irinotecan and mercaptopurine. However, many of the drug components we identified have not yet been well-recognized for anticancer therapy, for example, benfotiamine, telenzepine and chlorcyclizine. In addition, pimethixene, etamivan and mafenide have never been studied as anticancer drugs.

Discussion

The pathogenesis and progression of ICC is a complex, multifaceted process. During this process, various factors interact with each other. While ICC pathogenesis is an important area of study, the discovery of specific therapeutic drug targets for ICC is a more urgent issue that needs to be addressed. An in-depth understanding of drugs that can target and regulate ICC, specifically molecular-targeted therapeutic drugs, is very important for antitumor treatment.

Recently, many ICC tissues and cells from various stages have been collected as samples for whole genome microarray testing. In the present study, we integrated whole genome microarray information from ICC tissues of different origins using a modular bioinformatics approach. Modular analysis is an invaluable bioinformatics technology for functional analysis of single proteins, large-scale protein interactions and gene interactions. It has also been used to study the pathogenesis of many diseases and assess their treatments.

Among the six functional modules identified by modular analysis, module 1 contained the greatest number of genes, with the most complex interaction between the genes. Module 1 was involved in the response to hypoxia, response to drugs, protein tetramerization, response to oxidative stress and glucose metabolism. PTM (post-translational modification) is a key step in protein synthesis (9), involving the phosphorylation, glycosylation and methylation of proteins. One protein may have multiple different post-translational modifications (10). Although PTM has recently gained more appreciation, at present, PTM is rarely investigated in the context of antitumor drugs. The role of PTM has been appreciated in multiple myeloma, breast cancer and lung adenocarcinoma (1113); however, important PTMs have not been reported for ICC or even CCA. The physiology of tumor cells is very different compared to normal cells, as highlighted by differences in cell growth and differentiation, cell metabolism, biochemical composition and gene regulation. Inhibition of tumor growth is an important facet of anticancer treatment. The change in SSTR1 (somatostatin receptor type 1) in module 2 further confirmed that the modules had the core functions of ICC. Interestingly, we identified significant differences in the ADI1 gene. Because ADI1 has only been previously implicated in prostate cancer, this provides a new direction for future studies (14). Due to its important function in many biological activities, the insulin-like growth factor (IGF) signaling system has attracted wide attention. IGF signaling is associated with the onset of many diseases, such as cretinism, diabetes and cirrhosis of the liver. In module 3, IGFBP5, which is closely associated with the pathogenesis of many cancers, was highly expressed. IGFBP5 suppresses cancer in human melanoma cells (15) and plays important roles in human uterine leiomyoma cells, breast cancer, urothelial carcinoma and papillary thyroid carcinoma (1618). The role of tumor metabolism has also been appreciated in recent years. The canonical PI3K/AKT/mTOR pathway has been widely studied, as has BNIP3 (19,20). We hypothesize that ICC may inhibit the PI3K/AKT/mTOR pathway by upregulating BnIP3. However, this assertion must be verified by future experiments. The functions of modules 5 and 6 were related to cytoplasm regulation. The survival and prognosis of tumor patients are associated not only with tumor cells but also with the microenvironment of tumor cells. Thus, the pathogenesis of a tumor may be the result of the abnormal microenvironment of tumor cells. In these modules, PTMs, SSTR1, IGFBP5 and BNIP3 were most closely associated with the pathogenesis and progression of the cancer, indicating that these modules are tumor-related gene clusters, with these key genes at their cores. Moreover, in each gene cluster, multiple genes work synergistically to play an even more powerful role in tumor progression. Therefore, in antitumor therapy, the effect of a targeted drug treatment on a gene cluster is very likely better than the effect of targeting individual genes.

There has been a long history of traditional chemotherapy for ICC, with more and more drugs emerging with anticancer effects. However, because ICC itself is not sensitive to chemotherapy, there is no single effective drug, such as sorafenib for liver cancer. To screen for drug components with potential anticancer effects, we used the CMap small molecule database to identify drug components that regulated the core modules of ICC. Among the 65 drug components we identified, some drugs were found to have anticancer effects, as reported in extensive studies. For example, daunorubicin is effective for treating breast cancer (21), while irinotecan is a new drug for treating colon cancer that can strongly inhibit tumor growth and metastasis. Several studies have shown that irinotecan also has significant effects on colorectal liver metastases and metastatic colorectal cancers (22,23). Thus, many of the drug components identified through our screening have known anticancer effects, confirming the accuracy of the screening approach. Apart from the drugs that have been confirmed to have anticancer effects, literature mining identified 12 drug components with limited cancer-related data. For example, pimethixene, etamivan and mafenide have not been reported in any cancer-related studies. In our study, at least one of the drugs that regulated each module was from one of these 12 drug components. Through the Fisher's test, we determined the value of these drug components. For example, mercaptopurine, which regulated module 3, was scored first with daunorubicin, which regulated module 6, listed second. The Fisher's test indicated that the 12 identified drug components not previously reported in the anticancer literature had a strong ICC-regulating effect. This indicates that these drugs may have potent anticancer effects.

In summary, we analyzed two sets of mRNA and miRNA data for ICC. Using modular analysis, we identified six modules that were associated with ICC. Furthermore, based on the CMap database, we found 65 effective small molecule drug components that regulated the six modules. We accurately identified the most meaningful components using the Fisher's test. We then performed a detailed retrospective analysis to select 12 components from the small molecule drug components as new anticancer drug options. Our detailed target gene prediction via differential miRNA analysis and miRNA-mRNA correlation analysis revealed the regulatory roles of miRNAs in the progression of ICC. Thus, we provide a new approach for the selection of anticancer drug. Future experiments will verify the drugs we identified.

Acknowledgments

This study was supported by the Graduate Student Research Innovation fund in Heilongjiang Province of China (grant no. YJSCX2012-201HLJ), Changjiang Scholars and Innovative Research Team in University (grant no. IRT1122), the national natural Scientific Foundation of China (grant nos. 81272705, 81472322 and 81201878), China Postdoctoral Science Foundation Grant (2014M560271) and the natural Science Foundation of Heilongjiang Province (LC201437/H1617). The funders had no role in the study design, data collection and analysis, decision to publish, and preparation of the manuscript.

References

1 

Francis H, Alpini G and DeMorrow S: Recent advances in the regulation of cholangiocarcinoma growth. Am J Physiol Gastrointest Liver Physiol. 299:G1–G9. 2010. View Article : Google Scholar : PubMed/NCBI

2 

Tyson GL and El-Serag HB: Risk factors for cholangiocarcinoma. Hepatology. 54:173–184. 2011. View Article : Google Scholar : PubMed/NCBI

3 

Gatto M, Bragazzi MC, Semeraro R, Napoli C, Gentile R, Torrice A, Gaudio E and Alvaro D: Cholangiocarcinoma: Update and future perspectives. Dig Liver Dis. 42:253–260. 2010. View Article : Google Scholar : PubMed/NCBI

4 

Edge SB and Compton CC: The American Joint Committee on Cancer: The 7th edition of the AJCC cancer staging manual and the future of TNM. Ann Surg Oncol. 17:1471–1474. 2010. View Article : Google Scholar : PubMed/NCBI

5 

Mavros MN, Economopoulos KP, Alexiou VG and Pawlik TM: Treatment and prognosis for patients with intrahepatic cholangiocarcinoma: Systematic review and meta-analysis. JAMA Surg. Apr 9–2014.Epub ahead of print. View Article : Google Scholar : PubMed/NCBI

6 

Zhou H, Wang H, Zhou D, Wang H, Wang Q, Zou S, Tu Q, Wu M and Hu H: Hepatitis B virus-associated intrahepatic cholangiocarcinoma and hepatocellular carcinoma may hold common disease process for carcinogenesis. Eur J Cancer. 46:1056–1061. 2010. View Article : Google Scholar : PubMed/NCBI

7 

Sulpice L, Rayar M, Boucher E, Pele F, Pracht M, Meunier B and Boudjema K: Intrahepatic cholangiocarcinoma: Impact of genetic hemochromatosis on outcome and overall survival after surgical resection. J Surg Res. 180:56–61. 2013. View Article : Google Scholar

8 

Dodson RM, Weiss MJ, Cosgrove D, Herman JM, Kamel I, Anders R, Geschwind JF and Pawlik TM: Intrahepatic cholangiocarcinoma: Management options and emerging therapies. J Am Coll Surg. 217:736–750.e4. 2013. View Article : Google Scholar : PubMed/NCBI

9 

Olzmann JA, Brown K, Wilkinson KD, Rees HD, Huai Q, Ke H, Levey AI, Li L and Chin LS: Familial Parkinson's disease-associated L166P mutation disrupts DJ-1 protein folding and function. J Biol Chem. 279:8506–8515. 2004. View Article : Google Scholar

10 

Maita C, Maita H, Iguchi-Ariga SM and Ariga H: Monomer DJ-1 and its N-terminal sequence are necessary for mitochondrial localization of DJ-1 mutants. PLoS One. 8:e540872013. View Article : Google Scholar : PubMed/NCBI

11 

Shah SP, Lonial S and Boise LH: When cancer fights back: Multiple myeloma, proteasome inhibition, and the heat-shock response. Mol Cancer Res. 13:1163–1173. 2015. View Article : Google Scholar : PubMed/NCBI

12 

Yoo HM, Park JH, Jeon YJ and Chung CH: Ubiquitin-fold modifier 1 acts as a positive regulator of breast cancer. Front Endocrinol (Lausanne). 6:362015.

13 

Wan J, Zhan J, Li S, Ma J, Xu W, Liu C, Xue X, Xie Y, Fang W, Chin YE, et al: PCAF-primed EZH2 acetylation regulates its stability and promotes lung adenocarcinoma progression. Nucleic Acids Res. 43:3591–3604. 2015. View Article : Google Scholar : PubMed/NCBI

14 

Oram SW, Ai J, Pagani GM, Hitchens MR, Stern JA, Eggener S, Pins M, Xiao W, Cai X, Haleem R, et al: Expression and function of the human androgen-responsive gene ADI1 in prostate cancer. Neoplasia. 9:643–651. 2007. View Article : Google Scholar : PubMed/NCBI

15 

Wang J, Ding N, Li Y, Cheng H, Wang D, Yang Q, Deng Y, Yang Y, Li Y, Ruan X, et al: Insulin-like growth factor binding protein 5 (IGFBP5) functions as a tumor suppressor in human melanoma cells. Oncotarget. 6:20636–20649. 2015. View Article : Google Scholar : PubMed/NCBI

16 

Ling J, Jiang L, Zhang C, Dai J, Wu Q and Tan J: Upregulation of miR-197 inhibits cell proliferation by directly targeting IGFBP5 in human uterine leiomyoma cells. In Vitro Cell Dev Biol Anim. 51:835–842. 2015. View Article : Google Scholar : PubMed/NCBI

17 

Güllü G, Peker I, Haholu A, Eren F, Küçükodaci Z, Güleç B, Baloglu H, Erzik C, Özer A and Akkiprik M: Clinical significance of miR-140-5p and miR-193b expression in patients with breast cancer and relationship to IGFBP5. Genet Mol Biol. 38:21–29. 2015. View Article : Google Scholar : PubMed/NCBI

18 

Kavalieris L, O'Sullivan PJ, Suttie JM, Pownall BK, Gilling PJ, Chemasle C and Darling DG: A segregation index combining phenotypic (clinical characteristics) and genotypic (gene expression) biomarkers from a urine sample to triage out patients presenting with hematuria who have a low probability of urothelial carcinoma. BMC Urol. 15:232015. View Article : Google Scholar : PubMed/NCBI

19 

Huang Y, Shen P, Chen X, Chen Z, Zhao T, Chen N, Gong J, Nie L, Xu M, Li X, et al: Transcriptional regulation of BNIP3 by Sp3 in prostate cancer. Prostate. 75:1556–1567. 2015. View Article : Google Scholar : PubMed/NCBI

20 

Jiang K, Wang W, Jin X, Wang Z, Ji Z and Meng G: Silibinin, a natural flavonoid, induces autophagy via ROS-dependent mitochondrial dysfunction and loss of ATP involving BNIP3 in human MCF7 breast cancer cells. Oncol Rep. 33:2711–2718. 2015.PubMed/NCBI

21 

Cherigo L, Lopez D and Martinez-Luis S: Marine natura products as breast cancer resistance protein inhibitors. Mar Drugs. 13:2010–2029. 2015. View Article : Google Scholar : PubMed/NCBI

22 

Van Bael K, Jansen Y, Seremet T, Engels B, Delvaux G and Neyns B: A case report of long-term survival following hepatic arterial infusion of L-folinic acid modulated 5-fluorouracil combined with intravenous irinotecan and cetuximab followed by hepatectomy in a patient with initially unresectable colorectal liver metastases. Case Rep Oncol Med. 2015:4720372015.

23 

Hayashi H, Beppu T, Sakamoto Y, Miyamoto Y, Yokoyama N, Higashi T, Nitta H, Hashimoto D, Chikamoto A and Baba H: Prognostic value of Ki-67 expression in conversion therapy for colorectal liver-limited metastases. Am J Cancer Res. 5:1225–1233. 2015.PubMed/NCBI

Related Articles

Journal Cover

January-2016
Volume 35 Issue 1

Print ISSN: 1021-335X
Online ISSN:1791-2431

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Sun B, Xie C, Zheng T, Yin D, Wang J, Liang Y, Li Y, Yang G, Shi H, Pei T, Pei T, et al: Selecting molecular therapeutic drug targets based on the expression profiles of intrahepatic cholangiocarcinomas and miRNA-mRNA regulatory networks. Oncol Rep 35: 382-390, 2016.
APA
Sun, B., Xie, C., Zheng, T., Yin, D., Wang, J., Liang, Y. ... Liu, L. (2016). Selecting molecular therapeutic drug targets based on the expression profiles of intrahepatic cholangiocarcinomas and miRNA-mRNA regulatory networks. Oncology Reports, 35, 382-390. https://doi.org/10.3892/or.2015.4330
MLA
Sun, B., Xie, C., Zheng, T., Yin, D., Wang, J., Liang, Y., Li, Y., Yang, G., Shi, H., Pei, T., Han, J., Liu, L."Selecting molecular therapeutic drug targets based on the expression profiles of intrahepatic cholangiocarcinomas and miRNA-mRNA regulatory networks". Oncology Reports 35.1 (2016): 382-390.
Chicago
Sun, B., Xie, C., Zheng, T., Yin, D., Wang, J., Liang, Y., Li, Y., Yang, G., Shi, H., Pei, T., Han, J., Liu, L."Selecting molecular therapeutic drug targets based on the expression profiles of intrahepatic cholangiocarcinomas and miRNA-mRNA regulatory networks". Oncology Reports 35, no. 1 (2016): 382-390. https://doi.org/10.3892/or.2015.4330