Next Article in Journal
Identification of Cancer Stem Cell Molecular Markers and Effects of hsa-miR-21-3p on Stemness in Esophageal Squamous Cell Carcinoma
Next Article in Special Issue
From a Philosophical Framework to a Valid Prognostic Staging System of the New “Comprehensive Assessment” for Transplantable Hepatocellular Carcinoma
Previous Article in Journal
Cancer Vaccines Co-Targeting HER2/Neu and IGF1R
Previous Article in Special Issue
Assessment of a High Sensitivity Method for Identification of IDH1 R132x Mutations in Tumors and Plasma of Intrahepatic Cholangiocarcinoma Patients
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 11
Issue 4
10.3390/cancers11040519
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Establishment and Characterization of a New Intrahepatic Cholangiocarcinoma Cell Line Resistant to Gemcitabine

by
Chiara Varamo
1,2,†,
Caterina Peraldo-Neia
3,†,
Paola Ostano
3,
Marco Basiricò
4,
Chiara Raggi
5,6,
Paola Bernabei
7,
Tiziana Venesio
8,
Enrico Berrino
8,9,
Massimo Aglietta
1,4,
Francesco Leone
1,4,‡ and
Giuliana Cavalloni
4,*,‡
1
Department of Oncology, University of Turin, 10100 Torino, Italy
2
Laboratory of Tumor Inflammation and Angiogenesis, Department of Oncology, Center for Cancer Biology, KU Leuven, B3000 Leuven, Belgium
3
Cancer Genomics Lab, Fondazione Edo ed Elvo Tempia, 13900 Biella, Italy
4
Division of Medical Oncology, Candiolo Cancer Institute, FPO-IRCCS, 10060 Candiolo, Torino, Italy
5
Center for Autoimmune Liver Diseases, Humanitas Clinical and Research Center, 20089 Rozzano, Italy
6
Dept. Medicina Sperimentale e Clinica, Università di Firenze, 50100 Florence, Italy
7
Flow Cytometry Center, Candiolo Cancer Institute FPO-IRCCS, 10060 Candiolo, Torino, Italy
8
Molecular Pathology Lab, Unit of Pathology, Candiolo Cancer Institute FPO-IRCCS, 10060 Candiolo, Torino, Italy
9
Department of Medical Sciences, University of Turin, Corso Dogliotti 14, 10126 Turin, Italy
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work as co-first authors.
‡
These authors contributed equally to this work as co-last authors.
Cancers 2019, 11(4), 519; https://0-doi-org.brum.beds.ac.uk/10.3390/cancers11040519
Submission received: 7 March 2019 / Revised: 5 April 2019 / Accepted: 7 April 2019 / Published: 11 April 2019
(This article belongs to the Special Issue Models of Experimental Liver Cancer)
Browse Figures
Versions Notes

Abstract

:
Intrahepatic cholangiocarcinoma (ICC) is one of the most lethal liver cancers. Late diagnosis and chemotherapy resistance contribute to the scarce outfit and poor survival. Resistance mechanisms are still poorly understood. Here, we established a Gemcitabine (GEM) resistant model, the MT-CHC01R1.5 cell line, obtained by a GEM gradual exposure (up to 1.5 µM) of the sensitive counterpart, MT-CHC01. GEM resistance was irreversible, even at high doses. The in vitro and in vivo growth was slower than MT-CHC01, and no differences were highlighted in terms of migration and invasion. Drug prediction analysis suggested that Paclitaxel and Doxycycline might overcome GEM resistance. Indeed, in vitro MT-CHC01R1.5 growth was reduced by Paclitaxel and Doxycycline. Importantly, Doxycycline pretreatment at very low doses restored GEM sensitivity. To assess molecular mechanisms underlying the acquisition of GEM resistance, a detailed analysis of the transcriptome in MT-CHC01R1.5 cells versus the corresponding parental counterpart was performed. Transcriptomic analysis showed that most up-regulated genes were involved in cell cycle regulation and in the DNA related process, while most down-regulated genes were involved in the response to stimuli, xenobiotic metabolism, and angiogenesis. Furthermore, additional panels of drug resistance and epithelial to mesenchymal transition genes (n = 168) were tested by qRT-PCR and the expression of 20 genes was affected. Next, based on a comparison between qRT-PCR and microarray data, a list of up-regulated genes in MT-CHC01R1.5 was selected and further confirmed in a primary cell culture obtained from an ICC patient resistant to GEM. In conclusion, we characterized a new GEM resistance ICC model that could be exploited either to study alternative mechanisms of resistance or to explore new therapies.

1. Introduction

Cholangiocarcinoma (CCA) is a malignant tumor that originates from the epithelial cells of the biliary tract. It is classified as extrahepatic (ECC) and intrahepatic (ICC) cholangiocarcinoma depending on the anatomical onset of the tumor [1]. Due to its aggressiveness, CCA is characterized by a high mortality rate and poor prognosis. It is a relatively rare cancer whose diagnosis is difficult; most CCA patients are diagnosed at a state of unresectable disease and their overall 5-year survival is about 5%. The therapeutic options are limited and the only curative treatment is surgical resection with a 5-year survival reported around 20% to 50% [2]. In unresectable locally advanced and metastatic disease, no significant improvement in treatment has been achieved over the last decade. Cisplatin and Gemcitabine (GEM) combination is considered the standard of treatment, according to a phase III trial [3]. Other chemotherapy combinations have been investigated (Oxaliplatin, Capecitabine or Fluoropyrimidine, Epirubicin, Taxanes, Mytomicin), but mainly in phase II of retrospective trials, with small and heterogeneous cohorts of patients. Unfortunately, the expected overall survival at 12 months from diagnosis is about 50% [4,5]. The phenomenon of resistance to the current therapies is inevitable; indeed, CCA is able to escape several known tolerance mechanisms [6]. These mechanisms include the acquisition of genetic mutations, changes in the metabolism responsible for drug inhibition and degradation, activation of survival signaling pathways, and alterations of drug uptake and export [7]. Further, the tumor microenvironment plays an important role in chemoresistance by releasing specific factors, altering the sensitivity to chemotherapy [8].
The development of chemoresistance, in particular to GEM, represents the major point of failure in CCA treatment, responsible for progression and mortality of patients.
GEM is a deoxycytidine nucleoside analog, which is introduced in the cells by nucleoside transporters and then phosphorylated by deoxicitydine kinase (dCK), and acts by directly interfering with DNA synthesis, ribonucleotide reductase (RRMs) inhibition, and blocking cell cycle progression at the G1/S-phase transition [9]. GEM resistance has been well studied in pancreatic cancer. In particular, the involvement of several factors has been demonstrated, such as deficiencies in drug uptake, activation of alternative pathways involved in DNA repair, resistance to apoptosis, or in epithelial-mesenchymal transition (EMT) [10]. However, few information concerning GEM resistance is available in CCA. Therefore, an extensive study of the mechanisms underlying CCA GEM resistance may help to identify key molecules involved in the onset of this phenomenon.
In the present study, we established a GEM resistant CCA cell line, the MT-CHC01R1.5 and provide a biological and molecular characterization. Our data offer a useful model of chemoresistance either to study new mechanisms of resistance or to explore alternative therapeutic strategies in CCA.

2. Results

2.1. Biological Characterization of a New Established Gemcitabine Resistant Cholangiocarcinoma Cell Line

First, the GEM resistant clone, named MT-CHC01R1.5, was obtained by exposing parental cells intermittently to escalating doses of the drug, starting from 10 nM of GEM to 1.5 μM for 9 months. After 3 months of being continuously cultured in the presence of 1.5 μM, the MT-CHC01R1.5 cells were considered stable.
Next, in order to confirm GEM resistance, both MT-CHC01 and MT-CHC01R1.5 cells were treated with escalating doses of the drug. Figure 1 shows the in vitro growth inhibition of sensitive and resistant cells after GEM treatment; MT-CHC01R1.5 cells presented a remarkable resistance at high doses (until 10 μM) compared to the parental cell line. GEM growth inhibition in MT-CHC01 cells was not dose-dependent. Interestingly, MT-CHC01R1.5 resistance was irreversible; indeed, after GEM deprivation for 10 days, the resistant clone was refractory when re-exposed to the drug for 72 h (Figure 1).
To determine the population doubling time of the GEM resistant clone, we compared the cell growth curves of parental MT-CHC01 and resistant MT-CHC01R1.5 cells. As shown in Figure 2A, the doubling time of MT-CHC01 and MT-CHC01R1.5 cells was 25.3 and 49.1 h, respectively. The difference in terms of cell growth was statistically significant at 72 h (p = 0.01). We confirmed the number of viable cells by quantifying ATP, a key indicator of metabolically active cells (CellTiter GLO assay). Indeed, the MT-CHC01 proliferation rate was 1.73 ± 0.25 compared to its resistant counterpart (Figure S1).
The difference was highly appreciable in the colony formation assay (Figure 2B); indeed, the number of colonies with more than 10 cells was significantly different between the GEM sensitive and resistant clones (82.66 ± 6.42 for MT-CHC01 vs. 23.66 ± 1.5 for MT-CHC01R1.5; p = 0.0001).
Cell cycle distribution in parental and GEM resistant cells was determined by using flow cytometric analysis. As shown in Figure 3, a significant S phase expansion was observed in MT-CHC01R1.5 cells compared to MT-CHC01 parental cells after 24 h of culture (77.24% vs. 58.82% of cells, respectively, p < 0.001), reflecting the enhanced doubling time of resistant clone. At the same time, a significant decrease of G0/G1 (25.85% vs. 18.70%, p = 0.001) and G2/M (13.63% vs. 4.62%, p = 0.001) phases was observed in resistant cells. No differences were found after 48 h of culture.
From the morphological point of view, no appreciable difference was observed between parental and resistant cells. In order to verify whether GEM resistance conferred a more invasive phenotype, migration and invasion tests were carried out. MT-CHC01 and MT-CHC01R1.5 cells were added in an invasion chamber under serum free conditions in the absence or in the presence of extracellular matrix. After 24 h of incubation, no differences in both the number of migrated or invaded cells was observed, thus suggesting that in our model, the chemotherapy resistance was not linked to a more aggressive phenotype in terms of mobility potential.
Furthermore, in order to investigate the in vivo tumor growth, both MT-CHC01 and MT-CHC01R1.5 cells were subcutaneously injected in mice. As shown in Figure 4, both cell lines were able to give rise to tumors and no statistically significant difference in terms of tumorigenic potential was found in the first three weeks (p = 0.4). Indeed, after 3 weeks, tumors reached a volume ranging from 372 to 1242 mm3 in MT-CHC01 xenografts and from 282 to 1007 mm3 in MT-CHC01R1.5 xenografts. Interestingly, at the fourth week, the growth curves clearly separated in favor of mice inoculated with GEM sensitive cells (p = 0.0001; mean volume 1477.9 ± 599.9 and 707.4 ± 376.4 mm3 for MT-CHC01 and MT-CHC01R1.5 xenografts, respectively).

2.2. Gemcitabine Resistant Cells Acquired Cross-Resistance to Other Chemotherapeutic Agents

Susceptibility to Carboplatin, 5-fluorouracil (5-FU), Oxaliplatin, and Trabectedin on MT-CHC01 and MT-CHC01R1.5 were analyzed by Cell-Titer-Glo assays. As shown in Figure 5A, all the cells were resistant to Carboplatin, while in MT-CHC01R1.5, GEM resistance conferred cross-resistance to 5-FU and Trabectedin agents (Figure 5B,D, respectively). In contrast, no cross-resistance to Oxaliplatin (Figure 5C) was revealed.
The same test on the primary GEM resistant 82.3 cells demonstrated a complete cross-resistance to all evaluated drugs (Figure 5).

2.3. Gemcitabine Resistant Cells Exhibited Specific Expression Profiles with Important Effects on Biological Processes and Molecular Pathways

To assess the molecular mechanisms underlying the acquisition of GEM resistance, a detailed analysis of the transcriptome in MT-CHC01R1.5 cells versus their corresponding parental counterpart was performed. Among 736 differentially expressed genes, obtained by stringent thresholds, MT-CHC01R1.5 presented 354 down-regulated genes and 382 up-regulated ones (Table S1 and Figure S2). The gene ontology was investigated considering up- and down-regulated genes separately. Among the up-regulated genes, the most significant biological processes were functionally linked to DNA related mechanisms (elongation, replication, and metabolism) and cell cycles (Table S2). This was further confirmed as shown in Table S3 by the analysis of the cellular components and molecular functions. Additionally, process network examination by Metacore software revealed several processes involved in the cell cycle (Core, S phase, G1/S) and DNA damage (DBS, MMR, BER/NER repairs) (Table 1).
Moreover, a pathway map investigation (Table 2) showed an enrichment of different key regulators of the cell cycle, metabolism of purine and pirymidine, metabolism of folic acid, negative regulation of HIF1-alpha, activation of oxygen reactives upon oxidative stress, cytoskeleton remodeling, immune response, and DNA damage.
On the contrary, down-regulated genes were globally involved in protein folding, angiogenesis, blood vessel development, xenobiotic metabolism, and, more importantly, in mechanisms related to nutritional changes (i.e., starvation), calcium ion trafficking, hormonal regulation, and hypoxia (Table S4). Table S5 summarizes the significant cellular components and molecular functions enriched for down-regulated genes. Notably, an exploration of the process networks and pathway maps enriched in down-regulated genes showed an augmentation of proteolysis, reproduction, signal transduction mediated by ESR-1 and WNT signaling as well as apoptosis processes, as shown in Table 3 and Table 4, respectively.
Next, a drug prediction analysis based on the first 200 up-regulated genes was performed. Interestingly, eight of these genes (TUBB2, MMP1, ADH1, TYSY, POLA2, CyclineE, RRM1, DDC) were potentially druggable by chemotherapeutic or antibiotic agents, such as Paclitaxel, Doxycycline, Fomepizole, Pemetrexed, Dacarbazine, Silibinin, Triapine, and Carbidopa, respectively (Table S6). Thereafter, based on this analysis and on literature data [11,12], Paclitaxel and Doxycycline efficacy was tested. Strikingly, both drugs had a statistically significant growth inhibitory effect on MT-CHC01R1.5 cells as shown in Figure 6.

2.4. Doxycycline Restored Gemcitabine Sensitivity in Resistant MT-CHC01R1.5 Cells

As mentioned before, drug prediction studies indicated that MMP1, up-regulated in MT-CHC01R1.5 cells, may be potential druggable by the antibiotic, Doxycycline (DOXY). In order to verify the hypothesis, the sensitivity of MT-CHC01R1.5 cells to GEM in the presence of different DOXY concentrations was assessed. Since that DOXY IC50 value was 17.0 ± 2.2 μg/mL, cells were treated with non-toxic (5–10 μg/mL) and toxic doses (30 μg/mL) (Figure 7). Interestingly, pretreated cells with non-toxic DOXY doses reestablished GEM sensitivity even at very low doses. In contrast, the toxic DOXY dose (30 μg/mL) had a slight effect in restoring GEM susceptibility. Nevertheless, it was able alone to reduce cells’ viability compared to those cells routinely grown in the presence of GEM at a 1.5 μM dose (Figure 7).

2.5. Gemcitabine Resistant Cells Displayed a Modulation of Drug Resistance and Epithelial to Mesenchymal Transitions Related Genes

Next, to corroborate our previous in vitro findings, an additional molecular characterization of MT-CHC01R1.5 cells compared to their parental counterpart was achieved by RT2 profiler PCR arrays specific for drug resistance and EMT related pathways. Interestingly, 25% of drug resistance genes in MT-CHC01R1.5 cells were upregulated and these include key molecules, such as ESR2 (estrogen receptor 2), FGF2 (fibroblast growth factor 2), and SULT1E1 (sulfotransferase family 1E member) genes. Only 7 out of 84 (8.3%) were down-regulated (Figure 8A). The same analysis was performed using an EMT related genes panel. Among the 20 up-regulated genes (23.8%) (Figure 8B), the most overexpressed were COL5A2 (collagen type V alpha 2 chain), RGS2 (regulator of G-protein signaling 2), SNAI2 (snail family transcriptional repressor 2), TFPI2 (TFPI2 tissue factor pathway inhibitor 2), TGFB1 (transforming growth factor beta 1), and VIM (vimentin), thus suggesting an aggressive molecular phenotype of MT-CHC01R1.5 cells. Only 2 out of the 84 genes (2.4%) were down-regulated.
To technically validate the gene expression data, we selected nine genes of the up-regulated panel through microarray and qRT-PCR experiments; some of them overlapped between the two analyses, and the expression of all the selected genes was confirmed. Simultaneously, we evaluated the expression of the same group of genes in a primary cell culture isolated in our laboratory. This culture, named 82.3, was established from a patient primarily resistant to GEM. Interestingly, the same trend in gene expression was found (Table S7).

2.6. Methylation Assay

To evaluate the potential modification in the methylation pattern after the establishment of GEM resistance in the cell line, the well-known MGMT promoter methylation marker [13,14] was used to assess the status of both sensitive MT-CHC01 and resistant MT-CHC01R1.5 clones. As shown in Figure 9, neither the sensitive nor the resistant clone showed any methylation in all seven promoter CpG sites (mean percentage of the thymine representing the unmethylated cytosine of 100%), thus indicating that the GEM resistance of MT-CHC01R1.5 cells was not influenced by DNA methylation changes.

3. Discussion

In the present study, we investigated the important issue of acquired therapy resistance in human CCA. By continuous long-time exposure to Gemcitabine of sensitive MT-CHC01 cells, final surviving drug-resistant cells, MT-CHC01R1.5 cells, were further functionally and molecularly characterized. Importantly, acquired GEM-resistance was an irreversible phenomenon, as shown in MT-CHC01R1.5 cells when deprived of GEM for 10 days and lately re-exposed at high doses to the drug (up to 10 μM). This is also confirmed by other described CCA multidrug resistant models [15]. In contrast to the “on/off drug exposure system” described in several pancreatic and CCA models [15,16], the irreversible GEM-resistance of MT-CHC01R1.5 cells represents the key feature to ensure a reliable in vitro model for subsequent drug-resistance related analyses.
Biologically, the GEM resistant clone presented an increase in doubling time (25.3 vs. 49.1 h for MT-CHC01 and MT-CHC01R1.5, respectively) as already observed in other GEM-resistant models [15,17,18,19], including CCA. This difference could be due to a cell cycle perturbation, as S phase expansion of MT-CHC01R1.5, which may protect cells from chemotherapeutic activity, but in the meantime, there was a longer phase of DNA repair after GEM damage induction [19]. Concordantly, a statistically significant decrease in the tumor growth of MT-CHC01R1.5 was observed after four weeks from the in vivo injection compared to the parental counterpart.
Despite an early CCA response to GEM/platinum-based combination chemotherapy, most patients relapse as a result of GEM resistance onset that is often associated with a cross-resistance to other chemotherapeutic agents commonly used in the second line of treatment. Indeed, Wattanawongdon and colleagues established two different GEM resistant CCA cell lines, which also showed cross-resistance to 5-FU, Doxorubicin, and Paclitaxel [15]. In our study, we were able to confirm that acquired GEM resistance in MTCHC01R1.5 cells also conferred cross-resistance to other drugs, such as Carboplatin, 5-FU, and Trabectedin, except for Oxaliplatin. Furthermore, we expanded the analysis to primary GEM resistant 82.3 cells that were resistant to all of the already mentioned chemotherapeutic agents, including Oxaliplatin. This development of multidrug resistance suggests that other therapeutic options should be explored.
Recent advances in the treatment of several tumors have implicated the introduction of target therapies, including monoclonal antibodies and tyrosine kinase inhibitors [20]. However, in CCA, these novel therapeutic approaches have failed to provide clinical benefit [4,21,22,23].
The drug prediction analysis applied to the first 200 over-expressed genes in MT-CHC01R1.5 cells highlighted eight genes as potential therapeutic targets of chemotherapy or antibiotics. Based on this evidence, we tested the efficacy of Paclitaxel and Doxycycline, which have as targets TUBB2 and MMP1, respectively. These drugs were able to overcome GEM resistance even at low doses, comparable to those already tested in other cancer models, such as CCA, gallbladder, and ovarian cancers [24,25,26,27]. In particular, pretreatment of cells with non-toxic concentrations of Doxycycline was able to restore sensitivity even at low doses of Gemcitabine on MT-CHC01R1.5 cells.
In a recent study, it was demonstrated that albumin-bound Paclitaxel particles (nab-paclitaxel) exhibited antitumor activity in monotherapy and synergistic activity in combination with GEM treatment in metastatic pancreatic cancer [28]. Nab-paclitaxel was also FDA-approved in advanced breast cancer after failure of previous chemotherapy, and in combination with carboplatin in non-small cell lung cancer [29,30]. In a restricted series of CCA patients that failed first-line treatment, nab-paclitaxel activity was evaluated in combination with GEM or with fluoropyrimidine. The results demonstrated that patients treated with nab-paclitaxel in combination with chemotherapy showed a median overall survival of 9 months and a progression free survival of 6 months; the median time of survival after diagnosis of advanced disease was 21.5 months [31]. Furthermore, in preclinical investigation, Paclitaxel was able to inhibit EMT in CCA cells [24].
Despite the common use of Doxycycline as an antibiotic, there is some evidence of its antitumor activity, in particular, an inhibition of EMT and induction of apoptosis in different cancer types [12,32,33,34]. To date, there are no other data available on the possibility of overcoming GEM resistance with Doxycycline in other tumor models. However, it has been shown in breast cancer models that Doxycycline targets a sub-population of hypoxia-induced cancer stem cells that are Paclitaxel-resistant, overcoming hypoxia-induced drug-resistance [35]. Similarly, Doxycycline was shown to sensitize cells to tumor necrosis factor and Paclitaxel-induced apoptosis in a pancreatic in vitro model [36]. Thus, both Paclitaxel and Doxycycline are promising therapeutic agents in treatment after the first line GEM.
Several scientific studies show that the promotion of drug resistance mechanisms involves alterations in cell membrane transporters and cell distribution within the cell cycle, activation of mechanisms of DNA repair and/or anti-apoptotic processes, as well as modifications in drug metabolism, and EMT. These mechanisms can act independently or in combination through different signal transduction pathways [8].
As a matter of fact, the transcriptomic analysis revealed a perturbation of genes involved in cell cycle, in particular, in phase G1/S (i.e., MCMs, histone related genes, Rad51, PCNA, POLA2, DNA ligase I, E2F1/2), activation of DNA repair mechanisms (histone related genes, XRCC2, Rad51, PCNA, FEN1, RAD54L, FANCD2, FANCA, DNA ligase I, ChAF1 subunit A/B, RFC5, EXO1, RFC2, POLD cat (p125)), and nucleoside metabolism (RRM1, RRM2). Importantly the most significant up-regulated gene was RRM1, a well-known biomarker of GEM resistance [37]. Taken together, these data suggest that MT-CHC01R1.5 cells were able to protect themselves from the drug’s cytotoxic effects. Simultaneously, a down-regulation of genes related to apoptosis (SHIP, ActRIIB, c-Fos), regulation of the angiogenetic process (Ephrin-A, Ephrin-A1, IL-15, Angiotensin II/III, HB-EGF, Angiogenin, VEGF-A), and xenobiotic metabolism (CYPs, AKR7A3, GGT1, RORA) was evidenced. The reductions at the mRNA level of the latter two processes were not in line with those described in the literature. Besides the already well-established markers of GEM resistance, such as RRM1, RRM2, and RAD51 [17,38], which were also confirmed in our analysis, we also explored two additional panels of genes (n = 168) related to drug resistance and EMT. Among the first panel, the most perturbed genes were ESR2, FGF2, and SULT1E1 (up-regulated) and AR, BCL2L1, CYP2B6, CYP2E1, ELK1, NFKB2, NFKBIE (down-regulated), while within the EMT gene list, COL5A2, RGS2, SNAI2, TFPI2, TGFB1, VIM (up-regulated), and BMP4 and MMP9 (down-regulated) were the best represented in MT-CHC01R1.5 cells. Next, based on a comparison between the qRT-PCR and microarray data, a list of nine up-regulated genes in MT-CHC01R1.5 (BMP2, BRCA1, DHFR, JAG1, MMP1, MSH2, SPOCK2, TYMS, VIM) were selected and further confirmed in another in vitro model established in our group, obtained from an ICC patient resistant to GEM. Notably, some of these genes were already described to be linked to GEM resistance, such as VIM, a mesenchymal marker that is overexpressed in pancreatic GEM resistant cell lines [39,40]. TYMS played an important role in GEM metabolism; in sensitive cells, its expression was down-regulated by the drug, leading to apoptosis [41]. An up-regulation of this gene could cause an escape of cells from death [42]. High levels of TYMS correlated with poor responsiveness to 5-FU [43], and this correlated with MT-CHC01R1.5 unresponsiveness to this drug. JAG1, one of the ligands of NOTCH, has a proven role in cancer progression and in particular in EMT [44]. In the work of Wang and colleagues, JAG1 was overexpressed in a GEM resistant pancreatic cell line compared to its parental counterpart [40]. Moreover, in a recent work, Kim and colleagues demonstrated that in a panel of CCA cell lines, GEM resistance was associated to high levels of simultaneous expression of both APEX1 and Jagged 1 proteins, suggesting these molecules are potential therapeutic targets for CCA chemoresistance [45]. The tumor suppressor BRCA gene plays a key role in the repair of fatal DNA double strand breaks and is involved in the development of several tumors [46]. Interestingly, it was demonstrated that the BRCA/FA pathway genes and proteins were upregulated in a GEM resistant model of CCA [47].
Other genes are known to be responsible for the drug resistance phenomenon, but are not directly linked to GEM: For instance, BMP2, one of the BMP isoforms, is overexpressed in many cancer types [48,49,50] and its role is controversial. Indeed, it is a tumor suppressor in colon cancer [51], but in ovarian cancer, it promotes in vivo chemotherapy resistance by inhibition of in vitro proliferation [52]. MMP1 is found to be overexpressed in breast cancer cells resistant to Doxorubicin and is involved in the development of multidrug resistance [53,54]. MSH2, a member of the mismatch repair machinery, is found to be up-regulated in Paclitaxel resistance cell lines and, if silenced, its sensitivity is restored [55]. The role of DHFR in drug resistance is not well elucidated; Chauhan and colleagues hypothesized that its up-regulation could be implicated in Methotrexate resistance in acute lymphoblastic leukemia [56].
Additionally, our analysis showed an up-regulation of SPOCK2, a member of the SPARC family, in resistant cells. This data is in contrast with the work of Sherman-Baust and colleagues describing the down-regulation of mRNA expression in ovarian cancer cell lines resistant to Cisplatin and Doxorubicin [57].
Until now, cell lines remain the most useful and simple tool to study cancer development and drug response/resistance mechanisms.

4. Materials and Methods

4.1. Drugs

Gemcitabine hydrochloride was obtained from SANDOZ (Novartis Division, Italia), Trabectedin from PharmaMar (Pharma Mar, S.A., Madrid, Spain), 5-FU and Carboplatin from TEVA (Teva Italia srl, Milano, Italy), and Oxaliplatin from SUN Ranbaxy (Sun Pharmaceutical Industries Ltd., Goregaon, Mumbai, India). Paclitaxel and Doxycycline were obtained from Sigma Aldrich (Sigma–Aldrich, St. Louis, MO, USA). All drugs were dissolved in water for injection and aliquoted in different working solutions to avoid degradation.

4.2. Establishment of Gemcitabine Resistant Cell Line

The previously characterized CCA cell line, MT-CHC01 [58], was cultured in knockout/DMEM/F-12 complete medium containing 10% fetal bovine serum (FBS) (all from Sigma–Aldrich-Merk, Darmstadt, Germany), 100 U/mL penicillin, and 100 μg/mL streptomycin (Life Technologies Gathersburg, MD, USA). The GEM resistant clone, named MT-CHC01R1.5, was obtained by exposing parental cells intermittently to escalating doses of the drug, starting from 10 nM of GEM to 1.5 μM for 9 months. After another 3 months of being continuously cultured in the presence of 1.5 μM, the MT-CHC01R1.5 cells were considered stable.

4.3. Gemcitabine Resistant ICC Primary Cell Culture

The primary cell culture, 82.3, was obtained from a tumor specimen derived from the fourth generation of a patient-derived xenograft (PDX). The PDX was obtained from a tumor sample of a patient who underwent surgical resection for ICC. The biological material was obtained from a patient who has signed the informed consent, following institutional review board-approved protocols (001-IRCC-00IIS-10 “FPO-IRCCS, l’Istituto di Ricovero e Cura a Carattere Scientifico Candiolo (TO), Italy”). Briefly, the tumor was enzymatically digested with collagenase (200 U/mL; Sigma–Aldrich) for 3 h at 37 °C. Single cell suspension was obtained by filtering the supernatant through a 70-μm cell strainer (BD Biosciences, San Jose, CA, USA). Cells were finally re-suspended in a knockout/DMEM/F-12 complete medium in the presence of 20% FBS. After 20 passages, the immunophenotype of 82.3 cells was determined by flow cytometric analysis, evaluating CK19, CK7, and EpCAM. As shown in Figure S3, the characteristic epithelial markers are expressed. Moreover, cells are intrinsically resistant to GEM (Figure S4).

4.4. Flow Cytometry Analysis

Cells were washed in 1X PBS containing 0.1% bovine serum albumin (BSA, Sigma Aldrich, Saint Louis, MO, USA) and 0.01% sodium-azide. For cell permeabilization, when requested, the Fix and Perm reagent (BD Bioscience Milan Italy) was used following the manufacturer’s instructions. The following antibodies were used: FITC (Fluorescein isothiocyanate) conjugated mouse anti-CK7, CK19 (Abcam Cambridge, UK) and EpCAM (BD). The analysis was performed using Summit software (Beckman Coulter, Indianapolis, IN, USA). For cell cycle analysis, 1 × 106 cells were plated in 10 cm dishes for 24 and 48 h and fixed with 70% ethanol, on ice. After 24 h at −20 °C, ethanol was removed by centrifugation for 5 min at 300× g and cells were washed with PBS. The cells were then stained with 1 mL of Propidium Iodide (PI) staining solution (0.1% (v/v) Triton X-100 (Sigma-Aldrich), 10 μg/mL PI, and 100 μg/mL DNase-free RNase A (Sigma-Aldrich) in phosphate buffer saline (PBS)) and incubated in the dark at 37 °C for 15 min. Flow cytometry was performed using a CyAn™ ADP Flow Cytometer (Beckman Coulter, Indianapolis, IN, USA) and data were analyzed using FlowJo software, Version 7.6.3 (Treestar, Ashland, OR, USA). For each experiment, conducted in triplicate, 15,000 events were counted.

4.5. Cell Growth Curve

MT-CHC01 and MT-CHC01R1.5 cells (150,000/well × 2 mL) were seeded in 6-well culture plates (three-wells for each cell line) in the optimal medium. After 24, 48, and 72 h, the cells were detached and viable cells counted using an inverted microscope with 0.2% trypan blue dye. To calculate the population doubling time (DT), we used the following formula: DT = T ln2/ln(Xe/Xb), in which T is the time duration of the culture, Xe is the number of cells at the end of the incubation time, and Xb is the number of cells at the beginning of the incubation time.

4.6. Colony Formation Assay

Colonies derived from single cells were obtained by plating an equal number of MT-CHC01 and MT-CHC01R1.5 cells (500/well) in 24-well culture plates in optimal medium and then were incubated at 37 °C for five days. At the end of the period time, the medium was discarded and cells were stained with 0.1% crystal violet (Sigma-Aldrich). The number of colonies was counted for each well. The experiments were conducted three times with three points for each experiment.

4.7. Cell Growth Assay

The cell viability test was set up by plating 2500 cells per well into 96-well plates in the appropriate culture medium for 24 h. After 24 h, the cells were treated with escalating doses of GEM (from 10 μM to 0.078 μM). After 72 h, cell viability was evaluated with the Cell Titer-Glo® commercial kit (Promega Italia, Milan, Italy) following the manufacturer’s protocol. The measurement of luminescence was performed through the Glomax microplate reader (Glomax-Multi Detection System, Promega). The IC50 of the drug (concentration of drug needed to inhibit the 50% of growth) was calculated using the Calcusyn program (Biosoft, Cambridge, UK) based on the Chou-Talalay method.

4.8. In Vitro Migration and Invasion Assays

To test the motility, the transwell chambers assay was used (1 cm2/well, BD Falcon, BD Bioscience, Milan, Italy) either in the absence or in the presence of matrigel (Basement Membrane Matrix; BD MatrigelTM (BD Biosciences, San Jose, CA, USA) to evaluate the migration and invasion ability, respectively. The upper and lower cultures, in serum free conditions, were separated by 8-µm pore size poly-vinyl-pyrrolidone-free polycarbonate filters (BD Falcon). The experiments were carried out in triplicate.
After the incubation period, filters were fixed with methanol and stained with 0.5% crystal violet in 25% methanol; cells on the upper surface of the filters were removed using cotton swabs. Cells invading the lower surface were counted in five random fields and expressed as the number of invading cells per well.

4.9. In Vivo Tumor Growth Assay in NOD/SCID Mice

For in vivo studies, NOD (non-obese diabetic)/SCID (severe combined immunodeficient) female mice (4–6 weeks old) of about 20 to 25 gr/each (Charles River Laboratory, Rozzano, Milan, Italy) were maintained under sterile conditions in micro-isolator cages at the animal facilities of the IRCCS-Candiolo. All animal procedures were approved by the Institutional Ethical Committee for Animal Experimentation (Fondazione Piemontese per la Ricerca sul Cancro; 177/2015-PR and 178/2015-PR 24 March 2015) and by the Italian Ministry of Health. In three independent experiments, 6 mice were subcutaneously (s.c.) injected into the right flank under anesthesia (mixture of isoflurane and nitrous oxide) with either 3.0 × 106 of MT-CHC01 or MT-CHC01R1.5 cells in 50% growth factor-reduced BD Matrigel. Tumor diameters were measured weekly after cell injection up to tumor engraftment until 4 weeks after injection.

4.10. Gene Expression Analysis (GEP)

Total RNA was extracted from MT-CHC01 parental cells and from the corresponding GEM resistant MT-CHC01R1.5. Briefly, cells at 70% of confluence were lyzed using TRIZOL and scraped in ice. RNA extraction and purification was performed using the Absolutely RNA miRNA kit (Agilent Technologies, Santa Clara, CA 95051, USA), following the manufacturer’s instructions. RNA quantity was evaluated by Nanodrop, while the RNA integrity was evaluated by BioAnalyzer. One hundred ng of total RNA were processed using the Low Input Quick Amp Labeling Kit, one-color kit (Agilent Technologies). Six hundred ng of labeled RNA were hybridized on Whole Human Genome Microarray 4 × 44 K glass arrays. Arrays were scanned and images analyzed by the Feature Extraction Software from Agilent Technologies (version 10.7, Santa Clara, CA 95051, USA); raw data were then analyzed using the LIMMA (LInear Models for Microarray Analysis) package from Bioconductor. A threshold of |LogFC| greater than 2 and adjusted p-value < 0.001 were used to select modulated transcripts. Gene ontology analysis was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID); processes with a p-value < 0.01 being considered as statistically significant. Drug prediction, process network, and pathway map analyses were conducted using the Metacore tool (https://portal.genego.com). Microarray data were deposited in the Gene Expression Omnibus (GSE116118).

4.11. RT2 Profiler PCR Array

The Cancer Drug Resistance and the Epithelial to Mesenchymal Transition (EMT) arrays were used (Qiagen, Hilden, Germany); each 96-plate contained 84 genes and controls. Briefly, the RT2 First Strand Kit (Qiagen) was used for the reverse transcription of 0.5 µg of total RNA for each cell line. SYBR Green ROX qPCR Mastermix was used to perform the quantitative real time reaction. Data analysis was conducted using the ∆∆CT methods, as described in the protocol. Genes with |2−∆∆CT| values > 0.5 were considered modulated.

4.12. qRT-PCR Validation Assay

RNA was reverse-transcribed to cDNA with the high capacity cDNA reverse transcription kit (Applied Biosystem, Waltham, MA 02451, USA). Syber Green Mastermix was used for the round of amplification; in particular, cDNA was used for the detection of VIM, BMP2, MMP1, JAG1, DHFR, BRCA1, MSH2, SPOCK2, and TYMS deregulated genes and the PGK housekeeping gene with specific primers. Data analysis was conducted by comparing the gene expression values of MT-CHC01R1.5 with those of the parental cells, after normalization with the housekeeping PGK. Primer sequences and product lengths are summarized in Table S8.

4.13. MGMT Promoter Methylation Analyses

The methylation status of the MGMT promoter was evaluated by bisulfite-PCR and pyrosequencing. DNA bisulfite conversion was conducted using the commercial MGMT plus® (Diatech Pharmacogenetics, Jesi, Italy) kit. Briefly, 500 ng of genomic DNA were treated with sodium bisulte, inducing the deamination of the unmethylated cytosine into uracil without modifying the methylated ones. Three hundred ng of converted DNA was than amplified on a RotorGene Q (Qiagen, Hilden, Germany) instrument with a couple of primers encompassing 7 CpG sites of the MGMT promoter according to the manufacturer’s instructions; all conversion-derived uracil were replaced by thymine. The PCR condition for the MGMT gene was 95 °C for 5 min, 45 cycles of 95 °C for 30 s, 53 °C for 30 s, and 72 °C for 20 s, 72 °C for 5 min, and then green signal acquisition at 60 °C for 20 s. To evaluate the percentage of cytosine (e.g., methylated) or thymine (e.g., unmethylated) in each CpG site, Pyrosequencig analysis was carried out using the PyroMark Q96 ID system (Qiagen). Each sample was loaded two times for pyrosequencing (PSQ) and fully methylated and unmethylated DNA included in the kit were used as positive and negative controls in each experiment. PSQ analysis was performed with PyroMarker CpG software 1.0.11 (Qiagen). The software returns a mean methylation value for each 10 CpG site and the total mean of all 10 CpG sites. The samples were considered methylated when the mean of percentage of the cytosine in the 10 CpG sites was higher than 5%.

4.14. Statistical Analysis

The two ways Anova test was used to analyze cell and tumor growth as well as the differences in terms of drug response; a p-value less than 0.05 (confidence interval 95%) was considered as statistically significant. The one-way ANOVA was used to compare the response at different doses of chemotherapeutic agents within the same cell line, with a p-value less than 0.05 considered as statistically significant (confidence interval 95%). Data are expressed as the mean and standard deviation (SD) (bars) of values from at least triplicate assays. Statistical analysis was performed with Student t test. Asterisks indicate a significant p-value (* p < 0.01, ** p < 0.001, *** p < 0.0001, **** p< 0.00001).

5. Conclusions

Here, we provide a molecular and biologically functional characterization of our CCA resistant model by comparison with its parental sensitive counterpart. Further analyses focused on the methylation, microRNA, and proteome profiles as well as tailored functional and in vivo studies will be warranted to explore any possible effect of GEM and the potential targets to overcome this resistance.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6694/11/4/519/s1, Figure S1: ATP production of MT-CHC01 and MT-CHC01R1.5 cells, Figure S2: Volcano plot obtained with the deregulated genes in MT-CHC01R1.5 compared to MT-CHC01 parental cells. In blue, down-regulated genes, in red up-regulated genes, Figure S3: Immunophenotyping of 82.3 ICC primary cell culture, Figure S4: Percentage of viable cells after 72 h of gemcitabine exposure in 82.3 primary cell culture, Table S1: Differentially expressed genes in MT-CHC01R1.5 vs. MT-CHC01 parental cells, Table S2: Biological processes (level 5) enriched for up-regulated genes in MTCHC01R1.5 vs. MTCHC01 cells, Table S3: Cellular component and molecular functions enriched by up-regulated genes in MTCHC01R1.5 vs. MTCHC01 cells, Table S4: Pathways enriched for up-regulated genes in MTCHC01R1.5 vs. MTCHC01 cells, Table S5: Pathways enriched for down-regulated genes in MTCHC01R1.5 vs. MTCHC01 cells, Table S6: Drug prediction analysis: Eight potentially druggable genes from chemotherapeutic or antibiotic agents, Table S7: Expression data validation by qRT-PCR in MT-CHC01R1.5 and in 82.3 primary cell culture, Table S8: Primer sequences used for the qRT-PCR experiment of a selected panel of genes.

Author Contributions

C.V. conceived the study, performed the experiments and wrote the first article draft; C.P.-N. and G.C. conceived the study, performed the experiments and wrote the final article draft; P.O. analyzed gene expression data; M.B. performed the in vivo experiments; C.R. performed the statistical analyses; P.B. performed cytofluorimetric analysis, E.B. and T.V. performed the methylation analysis; M.A. and F.L. conceived the study and wrote the final article draft. All the authors read and approved the final work.

Funding

This work was supported by grant from Fondazione Piemontese per la Ricerca sul Cancro—Onlus (FPRC) 5X1000—“Identification of new druggable pathways in intrahepatic cholangiocarcinoma” Ministero della Salute 2010; Fondo per la ricerca locale (Linea B), LEOF_RIC_LOC_14_01 project title: “Transcriptomic and genetic analysis of paired primary and recurrent intrahepatic cholangiocarcinoma”; Progetto di Rete–NET-2011-02352137 “Genomic base triage for target therapy in colorectal cancer” Ricerca Sanitaria Finalizzata 2011; The Italian Ministry of Education, University and Research under the program FPRC “Dipartimenti di Eccellenza 2018–2022” (Project no. D15D18000410001) to the Department of Medical Sciences of the University of Turin.

Acknowledgments

We would like to thank Daniela Furlan for MGMT promoter methylation assay.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rizvi, S.; Khan, S.A.; Hallemeier, C.L.; Kelley, R.K.; Gores, G.J. Cholangiocarcinoma—Evolving concepts and therapeutic strategies. Nat. Rev. Clin. Oncol. 2018, 15, 95–111. [Google Scholar] [CrossRef]
  2. Valle, J.W.; Lamarca, A.; Goyal, L.; Barriuso, J.; Zhu, A.X. New Horizons for Precision Medicine in Biliary Tract Cancers. Cancer Discov. 2017, 7, 943–962. [Google Scholar] [CrossRef]
  3. Valle, J.W.; Furuse, J.; Jitlal, M.; Beare, S.; Mizuno, N.; Wasan, H.; Bridgewater, J.; Okusaka, T. Cisplatin and gemcitabine for advanced biliary tract cancer: A meta-analysis of two randomised trials. Ann. Oncol. 2014, 25, 391–398. [Google Scholar] [CrossRef]
  4. Eckel, F.; Schmid, R.M. Chemotherapy in advanced biliary tract carcinoma: A pooled analysis of clinical trials. Br. J. Cancer 2007, 96, 896–902. [Google Scholar] [CrossRef]
  5. Yang, R.; Wang, B.; Chen, Y.J.; Li, H.B.; Hu, J.B.; Zou, S.Q. Efficacy of gemcitabine plus platinum agents for biliary tract cancers: A meta-analysis. Anticancer Drugs 2013, 24, 871–877. [Google Scholar] [CrossRef] [PubMed]
  6. Nakamura, H.; Arai, Y.; Totoki, Y.; Shirota, T.; Elzawahry, A.; Kato, M.; Hama, N.; Hosoda, F.; Urushidate, T.; Ohashi, S.; et al. Genomic spectra of biliary tract cancer. Nat. Genet. 2015, 47, 1003–1010. [Google Scholar] [CrossRef]
  7. Pan, S.T.; Li, Z.L.; He, Z.X.; Qiu, J.X.; Zhou, S.F. Molecular mechanisms for tumour resistance to chemotherapy. Clin. Exp. Pharmacol. Physiol. 2016, 43, 723–737. [Google Scholar] [CrossRef] [PubMed]
  8. Hasan, S.; Taha, R.; Omri, H.E. Current Opinions on Chemoresistance: An Overview. Bioinformation 2018, 14, 80–85. [Google Scholar] [CrossRef] [Green Version]
  9. Huang, P.; Chubb, S.; Hertel, L.W.; Grindey, G.B.; Plunkett, W. Action of 2′,2′-difluorodeoxycytidine on DNA synthesis. Cancer Res. 1991, 51, 6110–6117. [Google Scholar]
  10. Adamska, A.; Elaskalani, O.; Emmanouilidi, A.; Kim, M.; Abdol Razak, N.B.; Metharom, P.; Falasca, M. Molecular and cellular mechanisms of chemoresistance in pancreatic cancer. Adv. Biol. Regul. 2018, 68, 77–87. [Google Scholar] [CrossRef] [PubMed]
  11. Fraveto, A.; Cardinale, V.; Bragazzi, M.C.; Giuliante, F.; De Rose, A.M.; Grazi, G.L.; Napoletano, C.; Semeraro, R.; Lustri, A.M.; Costantini, D.; et al. Sensitivity of Human Intrahepatic Cholangiocarcinoma Subtypes to Chemotherapeutics and Molecular Targeted Agents: A Study on Primary Cell Cultures. PLoS ONE 2015, 10, e0142124. [Google Scholar] [CrossRef]
  12. Chen, Z.; Wang, Y.; Liu, W.; Zhao, G.; Lee, S.; Balogh, A.; Zou, Y.; Guo, Y.; Zhang, Z.; Gu, W.; et al. Doxycycline inducible Krüppel-like factor 4 lentiviral vector mediates mesenchymal to epithelial transition in ovarian cancer cells. PLoS ONE 2014, 9, e105331. [Google Scholar] [CrossRef]
  13. Maroni, L.; Pierantonelli, I.; Banales, J.M.; Benedetti, A.; Marzioni, M. The significance of genetics for cholangiocarcinoma development. Ann. Transl. Med. 2013, 1, 28. [Google Scholar] [CrossRef]
  14. Yang, B.; House, M.G.; Guo, M.; Herman, J.G.; Clark, D.P. Promoter methylation profiles of tumor suppressor genes in intrahepatic and extrahepatic cholangiocarcinoma. Mod. Pathol. 2005, 18, 412–420. [Google Scholar] [CrossRef]
  15. Wattanawongdon, W.; Hahnvajanawong, C.; Namwat, N.; Kanchanawat, S.; Boonmars, T.; Jearanaikoon, P.; Leelayuwat, C.; Techasen, A.; Seubwai, W. Establishment and characterization of gemcitabine-resistant human cholangiocarcinoma cell lines with multidrug resistance and enhanced invasiveness. Int. J. Oncol. 2015, 47, 398–410. [Google Scholar] [CrossRef]
  16. Ioannou, N.; Seddon, A.M.; Dalgleish, A.; Mackintosh, D.; Solca, F.; Modjtahedi, H. Acquired resistance of pancreatic cancer cells to treatment with gemcitabine and HER-inhibitors is accompanied by increased sensitivity to STAT3 inhibition. Int. J. Oncol. 2016, 48, 908–918. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Wang, C.; Zhang, W.; Fu, M.; Yang, A.; Huang, H.; Xie, J. Establishment of human pancreatic cancer gemcitabine-resistant cell line with ribonucleotide reductase overexpression. Oncol. Rep. 2015, 33, 383–390. [Google Scholar] [CrossRef] [PubMed]
  18. Chen, S.Y.; Hu, S.S.; Dong, Q.; Cai, J.X.; Zhang, W.P.; Sun, J.Y.; Wang, T.T.; Xie, J.; He, H.R.; Xing, J.F.; et al. Establishment of paclitaxel-resistant breast cancer cell line and nude mice models, and underlying multidrug resistance mechanisms in vitro and in vivo. Asian Pac. J. Cancer Prev. 2013, 14, 6135–6140. [Google Scholar] [CrossRef] [PubMed]
  19. Guo, X.; Goessl, E.; Jin, G.; Collie-Duguid, E.S.; Cassidy, J.; Wang, W.; O’Brien, V. Cell cycle perturbation and acquired 5-fluorouracil chemoresistance. Anticancer Res. 2008, 28, 9–14. [Google Scholar] [PubMed]
  20. Bizama, C.; García, P.; Espinoza, J.A.; Weber, H.; Leal, P.; Nervi, B.; Roa, J.C. Targeting specific molecular pathways holds promise for advanced gallbladder cancer therapy. Cancer Treat. Rev. 2015, 41, 222–234. [Google Scholar] [CrossRef]
  21. Chiorean, E.G.; Ramasubbaiah, R.; Yu, M.; Picus, J.; Bufill, J.A.; Tong, Y.; Coleman, N.; Johnston, E.L.; Currie, C.; Loehrer, P.J. Phase II trial of erlotinib and docetaxel in advanced and refractory hepatocellular and biliary cancers: Hoosier Oncology Group GI06-101. Oncologist 2012, 17, 13. [Google Scholar] [CrossRef] [PubMed]
  22. El-Khoueiry, A.B.; Rankin, C.; Siegel, A.B.; Iqbal, S.; Gong, I.Y.; Micetich, K.C.; Kayaleh, O.R.; Lenz, H.J.; Blanke, C.D. S0941: A phase 2 SWOG study of sorafenib and erlotinib in patients with advanced gallbladder carcinoma or cholangiocarcinoma. Br. J. Cancer 2014, 110, 882–887. [Google Scholar] [CrossRef]
  23. Malka, D.; Cervera, P.; Foulon, S.; Trarbach, T.; de la Fouchardière, C.; Boucher, E.; Fartoux, L.; Faivre, S.; Blanc, J.F.; Viret, F.; et al. Gemcitabine and oxaliplatin with or without cetuximab in advanced biliary-tract cancer (BINGO): A randomised, open-label, non-comparative phase 2 trial. Lancet Oncol. 2014, 15, 819–828. [Google Scholar] [CrossRef]
  24. Hirose, A.; Tajima, H.; Ohta, T.; Tsukada, T.; Okamoto, K.; Nakanuma, S.; Sakai, S.; Kinoshita, J.; Makino, I.; Furukawa, H.; et al. Low-dose paclitaxel inhibits the induction of epidermal-mesenchymal transition in the human cholangiocarcinoma CCKS-1 cell line. Oncol. Lett. 2013, 6, 915–920. [Google Scholar] [CrossRef] [Green Version]
  25. Cadamuro, M.; Spagnuolo, G.; Sambado, L.; Indraccolo, S.; Nardo, G.; Rosato, A.; Brivio, S.; Caslini, C.; Stecca, T.; Massani, M.; et al. Low-Dose Paclitaxel Reduces S100A4 Nuclear Import to Inhibit Invasion and Hematogenous Metastasis of Cholangiocarcinoma. Cancer Res. 2016, 76, 4775–4784. [Google Scholar] [CrossRef]
  26. Wu, W.; Yu, L.H.; Ma, B.; Xu, M.J. The inhibitory effect of doxycycline on cisplatin-sensitive and -resistant epithelial ovarian cancer. PLoS ONE 2014, 9, e89841. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, S.; Liu, C.; Liu, X.; He, Y.; Shen, D.; Luo, Q.; Dong, Y.; Dong, H.; Pang, Z. Effects of matrix metalloproteinase inhibitor doxycycline and CD147 antagonist peptide-9 on gallbladder carcinoma cell lines. Tumour Biol. 2017, 39, 1010428317718192. [Google Scholar] [CrossRef] [PubMed]
  28. Von Hoff, D.D.; Ervin, T.; Arena, F.P.; Chiorean, E.G.; Infante, J.; Moore, M.; Seay, T.; Tjulandin, S.A.; Ma, W.W.; Saleh, M.N.; et al. Increased survival in pancreatic cancer with nab-paclitaxel plus gemcitabine. N. Engl. J. Med. 2013, 369, 1691–1703. [Google Scholar] [CrossRef] [PubMed]
  29. Cortazar, P.; Justice, R.; Johnson, J.; Sridhara, R.; Keegan, P.; Pazdur, R. US Food and Drug Administration approval overview in metastatic breast cancer. J. Clin. Oncol. 2012, 30, 1705–1711. [Google Scholar] [CrossRef] [PubMed]
  30. Satouchi, M.; Okamoto, I.; Sakai, H.; Yamamoto, N.; Ichinose, Y.; Ohmatsu, H.; Nogami, N.; Takeda, K.; Mitsudomi, T.; Kasahara, K.; et al. Efficacy and safety of weekly nab-paclitaxel plus carboplatin in patients with advanced non-small cell lung cancer. Lung Cancer 2013, 81, 97–101. [Google Scholar] [CrossRef]
  31. Unseld, M.; Scheithauer, W.; Weigl, R.; Kornek, G.; Stranzl, N.; Bianconi, D.; Brunauer, G.; Steger, G.; Zielinski, C.C.; Prager, G.W. Nab-paclitaxel as alternative treatment regimen in advanced cholangiocellular carcinoma. J. Gastrointest. Oncol. 2016, 7, 588–594. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. de Jong, M.C.; Hong, S.M.; Augustine, M.M.; Goggins, M.G.; Wolfgang, C.L.; Hirose, K.; Schulick, R.D.; Choti, M.A.; Anders, R.A.; Pawlik, T.M. Hilar cholangiocarcinoma: Tumor depth as a predictor of outcome. Arch. Surg. 2011, 146, 697–703. [Google Scholar] [CrossRef] [PubMed]
  33. Liu, S.; Liu, X.; Wang, H.; Zhou, Q.; Liang, Y.; Sui, A.; Yao, R.; Zhao, B.; Sun, M. Lentiviral vector-mediated doxycycline-inducible USP39 shRNA or cDNA expression in triple-negative breast cancer cells. Oncol. Rep. 2015, 33, 2477–2483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Son, K.; Fujioka, S.; Iida, T.; Furukawa, K.; Fujita, T.; Yamada, H.; Chiao, P.J.; Yanaga, K. Doxycycline induces apoptosis in PANC-1 pancreatic cancer cells. Anticancer Res. 2009, 29, 3995–4003. [Google Scholar] [PubMed]
  35. De Francesco, E.M.; Maggiolini, M.; Tanowitz, H.B.; Sotgia, F.; Lisanti, M.P. Targeting hypoxic cancer stem cells (CSCs) with Doxycycline: Implications for optimizing anti-angiogenic therapy. Oncotarget 2017, 8, 56126–56142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Fujioka, S.; Son, K.; Onda, S.; Schmidt, C.; Scrabas, G.M.; Okamoto, T.; Fujita, T.; Chiao, P.J.; Yanaga, K. Desensitization of NFκB for overcoming chemoresistance of pancreatic cancer cells to TNF-α or paclitaxel. Anticancer Res. 2012, 32, 4813–4821. [Google Scholar] [PubMed]
  37. Deng, T.; Pan, H.; Han, R.; Huang, D.; Li, H.; Zhou, L.; Wang, X.; Bai, M.; Li, X.; Liu, R.; et al. Gemcitabine sensitivity factors, hENT1 and RRM1 as potential prognostic biomarker for advanced biliary tract cancer. Int. J. Clin. Exp. Med. 2014, 7, 5041–5049. [Google Scholar]
  38. Tsai, M.S.; Kuo, Y.H.; Chiu, Y.F.; Su, Y.C.; Lin, Y.W. Down-regulation of Rad51 expression overcomes drug resistance to gemcitabine in human non-small-cell lung cancer cells. J. Pharmacol. Exp. Ther. 2010, 335, 830–840. [Google Scholar] [CrossRef] [PubMed]
  39. Shah, A.N.; Summy, J.M.; Zhang, J.; Park, S.I.; Parikh, N.U.; Gallick, G.E. Development and characterization of gemcitabine-resistant pancreatic tumor cells. Ann. Surg. Oncol. 2007, 14, 3629–3637. [Google Scholar] [CrossRef] [PubMed]
  40. Wang, Z.; Li, Y.; Kong, D.; Banerjee, S.; Ahmad, A.; Azmi, A.S.; Ali, S.; Abbruzzese, J.L.; Gallick, G.E.; Sarkar, F.H. Acquisition of epithelial-mesenchymal transition phenotype of gemcitabine-resistant pancreatic cancer cells is linked with activation of the notch signaling pathway. Cancer Res. 2009, 69, 2400–2407. [Google Scholar] [CrossRef]
  41. Ueno, H.; Kiyosawa, K.; Kaniwa, N. Pharmacogenomics of gemcitabine: Can genetic studies lead to tailor-made therapy? Br. J. Cancer 2007, 97, 145–151. [Google Scholar] [CrossRef]
  42. Colavito, D.; Cartei, G.; Dal Bianco, M.; Stecca, A.; Zustovich, F.; Dalle Carbonare, M.; Ragazzi, E.; Farina, M.; Colombrino, E.; Leon, A. Thymidylate synthetase allelic imbalance in clear cell renal carcinoma. Cancer Chemother. Pharmacol. 2009, 64, 1195–1200. [Google Scholar] [CrossRef]
  43. Mehrmohamadi, M.; Jeong, S.H.; Locasale, J.W. Molecular features that predict the response to antimetabolite chemotherapies. Cancer Metab. 2017, 5, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ntziachristos, P.; Lim, J.S.; Sage, J.; Aifantis, I. From fly wings to targeted cancer therapies: A centennial for notch signaling. Cancer Cell 2014, 25, 318–334. [Google Scholar] [CrossRef]
  45. Kim, H.B.; Cho, W.J.; Choi, N.G.; Kim, S.S.; Park, J.H.; Lee, H.J.; Park, S.G. Clinical implications of APEX1 and Jagged1 as chemoresistance factors in biliary tract cancer. Ann. Surg. Treat. Res. 2017, 92, 15–22. [Google Scholar] [CrossRef] [PubMed]
  46. D’Andrea, A.D.; Grompe, M. The Fanconi anaemia/BRCA pathway. Nat. Rev. Cancer 2003, 3, 23–34. [Google Scholar] [CrossRef] [PubMed]
  47. Nakashima, S.; Kobayashi, S.; Nagano, H.; Tomokuni, A.; Tomimaru, Y.; Asaoka, T.; Hama, N.; Wada, H.; Kawamoto, K.; Marubashi, S.; et al. BRCA/Fanconi anemia pathway implicates chemoresistance to gemcitabine in biliary tract cancer. Cancer Sci. 2015, 106, 584–591. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Zuo, W.H.; Zeng, P.; Chen, X.; Lu, Y.J.; Li, A.; Wu, J.B. Promotive effects of bone morphogenetic protein 2 on angiogenesis in hepatocarcinoma via multiple signal pathways. Sci. Rep. 2016, 6, 37499. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Chu, H.; Luo, H.; Wang, H.; Chen, X.; Li, P.; Bai, Y.; Zhang, F.; Cheng, R.; Chen, S.; Wang, Y.; et al. Silencing BMP-2 expression inhibits A549 and H460 cell proliferation and migration. Diagn. Pathol. 2014, 9, 123. [Google Scholar] [CrossRef] [PubMed]
  50. Katsuno, Y.; Hanyu, A.; Kanda, H.; Ishikawa, Y.; Akiyama, F.; Iwase, T.; Ogata, E.; Ehata, S.; Miyazono, K.; Imamura, T. Bone morphogenetic protein signaling enhances invasion and bone metastasis of breast cancer cells through Smad pathway. Oncogene 2008, 27, 6322–6333. [Google Scholar] [CrossRef] [Green Version]
  51. Vishnubalaji, R.; Yue, S.; Alfayez, M.; Kassem, M.; Liu, F.F.; Aldahmash, A.; Alajez, N.M. Bone morphogenetic protein 2 (BMP2) induces growth suppression and enhances chemosensitivity of human colon cancer cells. Cancer Cell Int. 2016, 16, 77. [Google Scholar] [CrossRef]
  52. Choi, Y.J.; Ingram, P.N.; Yang, K.; Coffman, L.; Iyengar, M.; Bai, S.; Thomas, D.G.; Yoon, E.; Buckanovich, R.J. Identifying an ovarian cancer cell hierarchy regulated by bone morphogenetic protein 2. Proc. Natl. Acad. Sci. USA 2015, 112, E6882–E6888. [Google Scholar] [CrossRef] [Green Version]
  53. Shen, C.J.; Kuo, Y.L.; Chen, C.C.; Chen, M.J.; Cheng, Y.M. MMP1 expression is activated by Slug and enhances multi-drug resistance (MDR) in breast cancer. PLoS ONE 2017, 12, e0174487. [Google Scholar] [CrossRef]
  54. Yang, J.M.; Xu, Z.; Wu, H.; Zhu, H.; Wu, X.; Hait, W.N. Overexpression of extracellular matrix metalloproteinase inducer in multidrug resistant cancer cells. Mol. Cancer Res. 2003, 1, 420–427. [Google Scholar]
  55. Zhang, J.; Yin, D.; Li, H. hMSH2 expression is associated with paclitaxel resistance in ovarian carcinoma, and inhibition of hMSH2 expression in vitro restores paclitaxel sensitivity. Oncol. Rep. 2014, 32, 2199–2206. [Google Scholar] [CrossRef]
  56. Chauhan, P.S.; Bhushan, B.; Singh, L.C.; Mishra, A.K.; Saluja, S.; Mittal, V.; Gupta, D.K.; Kapur, S. Expression of genes related to multiple drug resistance and apoptosis in acute leukemia: Response to induction chemotherapy. Exp. Mol. Pathol. 2012, 92, 44–49. [Google Scholar] [CrossRef]
  57. Sherman-Baust, C.A.; Becker, K.G.; Wood Iii, W.H.; Zhang, Y.; Morin, P.J. Gene expression and pathway analysis of ovarian cancer cells selected for resistance to cisplatin, paclitaxel, or doxorubicin. J. Ovarian Res. 2011, 4, 21. [Google Scholar] [CrossRef]
  58. Cavalloni, G.; Peraldo-Neia, C.; Varamo, C.; Casorzo, L.; Dell’Aglio, C.; Bernabei, P.; Chiorino, G.; Aglietta, M.; Leone, F. Establishment and characterization of a human intrahepatic cholangiocarcinoma cell line derived from an Italian patient. Tumour Biol. 2015, 37, 4041–4052. [Google Scholar] [CrossRef] [Green Version]
Cancers 11 00519 g001 550
Figure 1. Gemcitabine (GEM) effects on survival: Treatment was performed at 72 h using different GEM doses. MT-CHC01: parental cells; MT-CHC01R1.5: GEM resistant cells; MT-CHC01R1.5 no Gem: MT-CHC01R1.5 resistant clone deprived from GEM for 10 days and re-exposed to escalating doses of the drug for 72 h. All the experiments were conducted in three independent experiments in quadruplicate. Error bars represent mean and SD. **** p = 0.00001.
Figure 1. Gemcitabine (GEM) effects on survival: Treatment was performed at 72 h using different GEM doses. MT-CHC01: parental cells; MT-CHC01R1.5: GEM resistant cells; MT-CHC01R1.5 no Gem: MT-CHC01R1.5 resistant clone deprived from GEM for 10 days and re-exposed to escalating doses of the drug for 72 h. All the experiments were conducted in three independent experiments in quadruplicate. Error bars represent mean and SD. **** p = 0.00001.
Cancers 11 00519 g001
Cancers 11 00519 g002 550
Figure 2. Differences in cell growth and colony formation. (A) Growth curves of MT-CHC01 and MT-CHC01R1.5 cells: 1.5 × 105 cells were plated in 6-well plates in triplicate in three different experiments in optimal medium. Viable cells were counted at 24, 48, and 72 h after seeding. (B) Colony formation assay on MT-CHC01 and MT-CHC01R1.5 cells: Representative images of colony formation after 10 days of seeding. (C) Quantification of colony formation. Colonies formed by more than 10 cells were counted in triplicate in three different experiments. There is a statistically significant difference in terms of the number of colonies between MT-CHC01 and MT-CHC01R1.5 (p = 0.0001). Error bars represent the mean and SD. * p = 0.01, *** p = 0.0001.
Figure 2. Differences in cell growth and colony formation. (A) Growth curves of MT-CHC01 and MT-CHC01R1.5 cells: 1.5 × 105 cells were plated in 6-well plates in triplicate in three different experiments in optimal medium. Viable cells were counted at 24, 48, and 72 h after seeding. (B) Colony formation assay on MT-CHC01 and MT-CHC01R1.5 cells: Representative images of colony formation after 10 days of seeding. (C) Quantification of colony formation. Colonies formed by more than 10 cells were counted in triplicate in three different experiments. There is a statistically significant difference in terms of the number of colonies between MT-CHC01 and MT-CHC01R1.5 (p = 0.0001). Error bars represent the mean and SD. * p = 0.01, *** p = 0.0001.
Cancers 11 00519 g002
Cancers 11 00519 g003 550
Figure 3. Cell cycle analysis: (A) Representative flow cytometry histograms of cell populations in MT-CHC01 and MT-CHC01R1.5 cells assigned to different cell cycle phases (G0/G1, S, and G2) based on the intensity of PI staining (reflecting DNA content), gated on a single cell population. Cells were plated in optimal culture conditions (MT-CHC01R1.5 also in the presence of GEM) for 24 and 48 h and flow cytometry was performed, and data analyzed using FlowJo software. (B) Statistical analysis of the cell cycle distribution was conducted for three independent experiments. A significant increase of phase S was found in MT-CHC01R1.5 (p = 0.0001), while a significant decrease of phases G0/G1 and G2/M was revealed (p = 0.001) after 24 h. No significant differences were found after 48 h. ** p = 0.001 and **** p = 0.00001.
Figure 3. Cell cycle analysis: (A) Representative flow cytometry histograms of cell populations in MT-CHC01 and MT-CHC01R1.5 cells assigned to different cell cycle phases (G0/G1, S, and G2) based on the intensity of PI staining (reflecting DNA content), gated on a single cell population. Cells were plated in optimal culture conditions (MT-CHC01R1.5 also in the presence of GEM) for 24 and 48 h and flow cytometry was performed, and data analyzed using FlowJo software. (B) Statistical analysis of the cell cycle distribution was conducted for three independent experiments. A significant increase of phase S was found in MT-CHC01R1.5 (p = 0.0001), while a significant decrease of phases G0/G1 and G2/M was revealed (p = 0.001) after 24 h. No significant differences were found after 48 h. ** p = 0.001 and **** p = 0.00001.
Cancers 11 00519 g003
Cancers 11 00519 g004 550
Figure 4. In vivo tumor growth of MT-CHC01 and MT-CHC01R1.5 in NOD/SCID mice. The graph indicates the mean tumor volume (mm3) measured weekly (error bars: SD). Six mice for each cell lines were used and three independent experiments were conducted. **** p = 0.00001.
Figure 4. In vivo tumor growth of MT-CHC01 and MT-CHC01R1.5 in NOD/SCID mice. The graph indicates the mean tumor volume (mm3) measured weekly (error bars: SD). Six mice for each cell lines were used and three independent experiments were conducted. **** p = 0.00001.
Cancers 11 00519 g004
Cancers 11 00519 g005 550
Figure 5. Effects of several chemotherapeutic agents on MT-CHC01, MT-CHC01R1.5, and 82.3 cells. The different drugs ((A) Carboplatin; (B) 5-FU; (C) Oxaliplatin; (D) Trabectedin) at the indicated concentrations (for Carboplatin, from 40 to 0.156 μg/mL; for 5-FU, from 384 to 0.75 μM; for Oxaliplatin, from 20 to 0.078 μM; and for Trabectedin, from 10 to 0.039 nM) were added to each cell line. The effect was evaluated by the Cell Titer-Glo assay after 72 h. The values obtained are the mean with SD (bars) of three different experiments. Statistical analysis was performed with Student t test. Asterisks indicate a significant p value (* p < 0.01, ** p < 0.001, *** p < 0.0001, **** p < 0.00001). Statistical analysis was performed comparing resistant cells to MT-CHC01 cells.
Figure 5. Effects of several chemotherapeutic agents on MT-CHC01, MT-CHC01R1.5, and 82.3 cells. The different drugs ((A) Carboplatin; (B) 5-FU; (C) Oxaliplatin; (D) Trabectedin) at the indicated concentrations (for Carboplatin, from 40 to 0.156 μg/mL; for 5-FU, from 384 to 0.75 μM; for Oxaliplatin, from 20 to 0.078 μM; and for Trabectedin, from 10 to 0.039 nM) were added to each cell line. The effect was evaluated by the Cell Titer-Glo assay after 72 h. The values obtained are the mean with SD (bars) of three different experiments. Statistical analysis was performed with Student t test. Asterisks indicate a significant p value (* p < 0.01, ** p < 0.001, *** p < 0.0001, **** p < 0.00001). Statistical analysis was performed comparing resistant cells to MT-CHC01 cells.
Cancers 11 00519 g005
Cancers 11 00519 g006 550
Figure 6. Effects of Paclitaxel (A) and Doxycycline (B) on MT-CHC01R1.5 cells: MT-CHC01R1.5 cells were sensible at low doses both to Paclitaxel and Doxycycline. The viability was evaluated by Cell Titer Glo assay. All experiments were conducted three times in quadruplicate. ** p = 0.001, *** p = 0.0001, **** p = 0.00001.
Figure 6. Effects of Paclitaxel (A) and Doxycycline (B) on MT-CHC01R1.5 cells: MT-CHC01R1.5 cells were sensible at low doses both to Paclitaxel and Doxycycline. The viability was evaluated by Cell Titer Glo assay. All experiments were conducted three times in quadruplicate. ** p = 0.001, *** p = 0.0001, **** p = 0.00001.
Cancers 11 00519 g006
Cancers 11 00519 g007 550
Figure 7. Doxycycline restores the sensitivity of MT-CHC01R1.5 cells to Gemcitabine. Viability of Gemcitabine-resistant MT-CHC01 cells was tested in different culture conditions. GEM 1.5: cells routinely growing in the presence of 1.5 µM Gemcitabine for 72 h; DOXY 30-10-5: cells treated with 30 to 10 or 5 µg/mL of Doxycycline, respectively, for 72 h; DOXY 30-10-5 GEM 1.5-0.75-0.373: cells pre-treated with 30, 10, or 5 µg/mL of Doxycycline for 24 h, followed by 1.5, 0.75, or 0.375 µM Gemcitabine for 72 h; NT: not treated cells. The experiment was repeated three times with four replicates for each dose. Error bars represented means and SD. * p = 0.01, ** p = 0.001, **** p = 0.00001.
Figure 7. Doxycycline restores the sensitivity of MT-CHC01R1.5 cells to Gemcitabine. Viability of Gemcitabine-resistant MT-CHC01 cells was tested in different culture conditions. GEM 1.5: cells routinely growing in the presence of 1.5 µM Gemcitabine for 72 h; DOXY 30-10-5: cells treated with 30 to 10 or 5 µg/mL of Doxycycline, respectively, for 72 h; DOXY 30-10-5 GEM 1.5-0.75-0.373: cells pre-treated with 30, 10, or 5 µg/mL of Doxycycline for 24 h, followed by 1.5, 0.75, or 0.375 µM Gemcitabine for 72 h; NT: not treated cells. The experiment was repeated three times with four replicates for each dose. Error bars represented means and SD. * p = 0.01, ** p = 0.001, **** p = 0.00001.
Cancers 11 00519 g007
Cancers 11 00519 g008 550
Figure 8. Differentially expressed genes obtained using panels of cancer drug resistance (A) and EMT (B) related genes. Control: MT-CHC01 parental cells, which represent the reference.
Figure 8. Differentially expressed genes obtained using panels of cancer drug resistance (A) and EMT (B) related genes. Control: MT-CHC01 parental cells, which represent the reference.
Cancers 11 00519 g008
Cancers 11 00519 g009 550
Figure 9. MGMT promoter methylation level in (A) sensitive and (B) resistant cells. The highlighted regions represent the analyzed CpG site of the MGMT promoter with the percentage of C showing the level of methylation each site. The experiment was conducted in duplicate.
Figure 9. MGMT promoter methylation level in (A) sensitive and (B) resistant cells. The highlighted regions represent the analyzed CpG site of the MGMT promoter with the percentage of C showing the level of methylation each site. The experiment was conducted in duplicate.
Cancers 11 00519 g009
Table
Table 1. Enrichment by process networks of up-regulated genes.
Table 1. Enrichment by process networks of up-regulated genes.
Networksp-ValueGenes Involved
Cell cycle_S phase1.2 × 10−24MCM4, Histone H1.5, Histone H4, RFC5, MCM3, CHTF18, Rad51, Cdt1, MCM2, PCNA, FEN1, ORC1L, RFC2, MCM7, ChAF1 subunit A, MCM6, Histone H1, POLD cat (p125), POLA2, DNA ligase I, PRIM1, MCM5, MCM8, E2F1, CDC45L, MCM4/6/7 complex, ChAF1 subunit B
Cell cycle_Core6.9 × 10−13E2F2, MCM4, MCM3, Cdt1, MCM2, FEN1, Kinase MYT1, ORC1L, MCM7, MCM6, DNA ligase I, MCM5, MCM8, E2F1, CDC45L, MCM4/6/7 complex
DNA damage_DBS repair2.7 × 10−7Histone H4, XRCC2, Rad51, PCNA, FEN1, ChAF1 subunit A, RAD54L, FANCD2, FANCA, DNA ligase I, ChAF1 subunit B
DNA damage_MMR repair8.8 × 10−4RFC5, EXO1, PCNA, RFC2, POLD cat (p125)
DNA damage_Checkpoint9.8 × 10−4RFC5, CIA/ASF1, PCNA, Kinase MYT1, RFC2, FANCD2, E2F1
Cell cycle_G1-S0.001Rad51, TYSY, PCNA, b-Myb, UHRF1, Ceb1, FANCD2, E2F1
DNA damage_BER-NER repair0.002RFC5, PCNA, FEN1, RFC2, POLD cat (p125), DNA ligase I
Apoptosis_Apoptotic nucleus0.004Histone H1.5, Rad51, c-Myb, b-Myb, Histone H1, FANCD2, E2F1
Reproduction_Spermatogenesis, motility and copulation0.008Histone H1.1, PKC-beta, Histone H2, Histone H2A, PKC, DDC, Histone H1, FANCA
Cytoskeleton_Regulation of cytoskeleton rearrangement0.008Tubulin beta, Vimentin, Talin, PKC, Tubulin beta 1, Tubulin beta 2, Band 4.1-like protein 2
Table
Table 2. Enrichment by pathway maps of up-regulated genes.
Table 2. Enrichment by pathway maps of up-regulated genes.
Mapsp-ValueGenes Involved
Cell cycle_Start of DNA replication in early S phase2.68 × 10−14MCM4, MCM3, Cdt1, MCM2, ORC1L, Histone H1, MCM5, E2F1, CDC45L, MCM4/6/7 complex
Cell cycle_Transition and termination of DNA replication1.86 × 10−6MCM2, PCNA, FEN1, POLD cat (p125), DNA ligase I
NETosis in SLE3.7 × 10−6Histone H4, Histone H2, Histone H2A, PKC, Histone H1
dCTP/dUTP metabolism2.29 × 10−5Ribonucleotide reductase, RRM2, Small RR subunit, RRM1, POLD cat (p125), POLA2
dATP/dITP metabolism4.8 × 10−5Ribonucleotide reductase, RRM2, Small RR subunit, RRM1, POLD cat (p125), POLA2
Cell cycle_Sister chromatid cohesion6.3 × 10−4CHTF18, PCNA, Histone H1
Immune response_IFN-alpha/beta signaling via PI3K and NF-kB pathways6.8 × 10−4DHFR, PCNA, b-Myb, E2F1, ISG15
Cytoskeleton remodeling_Neurofilaments9.3 × 10−4Tubulin beta, Vimentin, NEFL
Transcription_Effect of Folic acid on genome stability0.001DHFR, HCP1, TYSY, ADH1
Immune response_IL-4-induced regulators of cell growth, survival, differentiation and metabolism0.001MCM4, MCM6, MMP-13, MCM5
Folic acid metabolism0.001C1TC, DHFR, HCP1, SLC19A1
DNA damage_Role of Brca1 and Brca2 in DNA repair0.001Rad51, PCNA, FANCD2
Development_Transcription regulation of granulocyte development0.001c-Myb, PKC, E2F1
Transcription_Negative regulation of HIF1A function0.002MCM3, MCM2, MCM7, MCM5
Inhibition of neutrophil migration by proresolving lipid mediators in COPD0.002TLN2, Talin, PKC, cPKC
Signal transduction_HTR2A signaling outside the nervous system0.002PKC-beta, PKC, cPKC (conventional), E2F1
DNA damage_Nucleotide excision repair0.003PCNA, POLD cat (p125), DNA ligase I
Cell adhesion_Cell-matrix glycoconjugates0.003CCL5, MMP-1, ECM1
Development_Gastrin in differentiation of the gastric mucosa0.003PKC-beta, PKC, cPKC
Complement pathway disruption in thrombotic microangiopathy0.003CCL5, PKC, cPKC
Immune response_IL-1 signaling pathway0.003CCL5, SPHK1, MMP-1, MMP-13
Oxidative stress_Activation of NOX1, NOX5, DUOX1 and DUOX2 NADPH Oxidases0.004PKC-beta, PKC, cPKC
Development_VEGF signaling and activation0.004PKC-beta, SPHK1, PKC
Glomerular injury in Lupus Nephritis0.005CCL5, PKC-beta1, MMP-1, PKC-beta2
Apoptosis and survival_DNA-damage-induced apoptosis0.005FANCD2, E2F1
Transcription_Assembly of RNA Polymerase II preinitiation complex on TATA-less promoters0.008DHFR, PCNA
Immune response_IL-16 signaling pathway0.009CCL5, PKC, cPKC
DNA damage_Mismatch repair0.009EXO1, PCNA
Table
Table 3. Enrichment by process networks of down-regulated genes.
Table 3. Enrichment by process networks of down-regulated genes.
Networksp-ValueGenes Involved
Reproduction_Gonadotropin regulation0.003ATF-3, EGR1, HB-EGF, GABA-A receptor epsilon subunit, Adenylate cyclase, FosB, AKR1C3, c-Fos
Proteolysis_Connective tissue degradation0.004Matrilysin (MMP-7), Trypsin II, Trypsin, ADAM8, SERPINA3 (ACT), Protein C inhibitor
Reproduction_FSH-beta signaling pathway0.004EGR1, TGM2, Adenylate cyclase, VEGF-A, ActRIIB, c-Fos, IBP
Signal transduction_ESR1-membrane pathway0.005HB-EGF, Adenylate cyclase, Adenylate cyclase type IV, c-Fos, Caveolin-1
Development_Regulation of angiogenesis0.007Ephrin-A, Ephrin-A1, IL-15, Angiotensin II, HB-EGF, Angiogenin, VEGF-A, Angiotensin III
Signal transduction_WNT signaling0.007Matrilysin (MMP-7), HES1, WISP2, HB-EGF, Adenylate cyclase, VEGF-A, c-Fos
Table
Table 4. Enrichment by pathway maps of down-regulated genes.
Table 4. Enrichment by pathway maps of down-regulated genes.
Mapsp-ValueGenes Involved
Protein folding and maturation_Angiotensin system maturation2.223 ×10−8Angiotensin (2–10), Angiotensin II, Angiotensin IV, Angiotensin I, Angiotensin (1–7), Angiotensinogen, Angiotensin III, Angiotensin (1–9)
Protein folding and maturation_Posttranslational processing of neuroendocrine peptides1.017 × 10−7NT, LargeNT, Trypsin, NN, LargeNN, NT/NN
Transcription_HIF-1 targets4.662 × 10−5TGM2, MCT4, IBP1, Stanniocalcin 2, VEGF-A, Cyclin G2, REDD1
Signal transduction_HTR2A signaling outside the nervous system9.750 × 10−5EGR1, HB-EGF, Adenylate cyclase, PLD1, HB-EGF (mature), Caveolin-1
Signal transduction_mTORC1 downstream signaling3.624 × 10−4eIF4A, VEGF-A, PDCD4, eIF4B, ULK1
Immune response_CD16 signaling in NK cells6.423 × 10−4PLA2, PLD1, SHIP, c-Fos, VAV-1
Reproduction_Gonadotropin-releasing hormone (GnRH) signaling7.808 × 10−4ATF-3, EGR1, Adenylate cyclase, FosB, c-Fos
L-Tryptophan metabolism (part 1)0.001SLC43A1, ASCT2 (SLC1A5), SLC7A8, SLC38A2, SLC6A14
Development_Leptin signaling via JAK/STAT and MAPK cascades0.001EGR1, VEGF-A, c-Fos
Apoptosis and survival_Apoptotic Activin A signaling0.001SHIP, ActRIIB, c-Fos
Putative pathways of hormone action in neurofibromatosis type 10.001EGR1, EGR2 (Krox20), FosB
Membrane-bound ESR1: interaction with G-proteins signaling0.002HB-EGF, Adenylate cyclase, c-Fos, Caveolin-1
Development_Non-genomic action of Retinoic acid in cell differentiation0.002TGM2, VEGF-A, c-Fos, VAV-1
IGF family signaling in colorectal cancer0.003FosB, VEGF-A, c-Fos, IBP
Regulation of lipid metabolism_FXR-dependent negative-feedback regulation of bile acids concentration0.003FGF19, CYP3A4, CYP2B6
Cytoskeleton remodeling_Fibronectin-binding integrins in cell motility0.004ITGAV, VAV-1, Caveolin-1
Development_Role of IL-8 in angiogenesis0.004HB-EGF, VEGF-A, c-Fos, Caveolin-1
Muscle contraction_Regulation of eNOS activity in endothelial cells0.004KLF2, FosB, VEGF-A, Caveolin-1
Development_Angiotensin activation of ERK0.004Angiotensin II, HB-EGF, c-Fos
Immune response_IL-6-induced acute-phase response in hepatocytes0.005Angiotensinogen, IBP1, c-Fos
Retinol metabolism0.005CYP3A5, CYP3A4, CYP2B6, CYP3A7
Immune response_IL-6 signaling pathway via JAK/STAT0.005Rac2, VEGF-A, c-Fos, VAV-1
Transport_ACM3 signaling in lacrimal glands0.006Mucin 5B, PLD1, Aquaporin 5, Lysozyme
Transcription_P53 signaling pathway0.007HSP27, VEGF-A, c-Fos
Role of tumor microenvironment in plexiform neurofibroma formation in neurofibromatosis type 10.007Rac2, VEGF-A, VAV-1
Development_Role of Activin A in cell differentiation and proliferation0.007p15, Adenylate cyclase, ActRIIB
Main pathways of Schwann cells transformation in neurofibromatosis type 1 0.008Rac2, VEGF-A, c-Fos, VAV-1
Development_Angiotensin signaling via PYK20.009Angiotensin II, c-Fos, VAV-1

Share and Cite

MDPI and ACS Style

Varamo, C.; Peraldo-Neia, C.; Ostano, P.; Basiricò, M.; Raggi, C.; Bernabei, P.; Venesio, T.; Berrino, E.; Aglietta, M.; Leone, F.; et al. Establishment and Characterization of a New Intrahepatic Cholangiocarcinoma Cell Line Resistant to Gemcitabine. Cancers 2019, 11, 519. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers11040519

AMA Style

Varamo C, Peraldo-Neia C, Ostano P, Basiricò M, Raggi C, Bernabei P, Venesio T, Berrino E, Aglietta M, Leone F, et al. Establishment and Characterization of a New Intrahepatic Cholangiocarcinoma Cell Line Resistant to Gemcitabine. Cancers. 2019; 11(4):519. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers11040519

Chicago/Turabian Style

Varamo, Chiara, Caterina Peraldo-Neia, Paola Ostano, Marco Basiricò, Chiara Raggi, Paola Bernabei, Tiziana Venesio, Enrico Berrino, Massimo Aglietta, Francesco Leone, and et al. 2019. "Establishment and Characterization of a New Intrahepatic Cholangiocarcinoma Cell Line Resistant to Gemcitabine" Cancers 11, no. 4: 519. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers11040519

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop