- Research article
- Open Access
- Open Peer Review
Establishment of a patient-derived intrahepatic cholangiocarcinoma xenograft model with KRAS mutation
BMC Cancer volume 16, Article number: 90 (2016)
Intrahepatic cholangiocarcinoma (ICC) is an aggressive, highly lethal tumors and lacks of effective chemo and targeted therapies. Cell lines and animal models, even partially reflecting tumor characteristics, have limits to study ICC biology and drug response. In this work, we created and characterized a novel ICC patient-derived xenograft (PDX) model of Italian origin.
Seventeen primary ICC tumors derived from Italian patients were implanted into NOD (Non-Obese Diabetic)/Shi-SCID (severe combined immunodeficient) mice. To verify if the original tumor characteristics were maintained in PDX, immunohistochemical (cytokeratin 7, 17, 19, and epithelial membrane antigen) molecular (gene and microRNA expression profiling) and genetic analyses (comparative genomic hybridization array, and mutational analysis of the kinase domain of EGFR coding sequence, from exons 18 to 21, exons 2 to 4 of K-RAS, exons 2 to 4 of N-RAS, exons 9 and 20 of PI3KCA, and exon 15 of B-RAF) were performed after tumor stabilization.
One out of 17 (5.8 %) tumors successfully engrafted in mice. A high molecular and genetic concordance between primary tumor (PR) and PDX was confirmed by the evaluation of biliary epithelial markers, tissue architecture, genetic aberrations (including K-RAS G12D mutation), and transcriptomic and microRNA profiles.
For the first time, we established a new ICC PDX model which reflects the histology and genetic characteristics of the primary tumor; this model could represent a valuable tool to understand the tumor biology and the progression of ICC as well as to develop novel therapies for ICC patients.
Cholangiocarcinoma (CCA) is the most common biliary tract neoplasm of the biliary tree, classified, according to its site of origin, as intrahepatic (ICC), perihilar or extrahepatic (ECC) cholangiocarcinoma [1, 2]. These subtypes differ in their biology, clinical-pathological characteristics and management. ICC accounts for approximately 10-15 % of CCA [3, 4], although its incidence is different worldwide with a higher incidence in Asia (96 per 100,000 in Thailand) , but is increasing also in other geographic regions . Several risk factors of CCA, include infectious and inflammatory diseases, congenital conditions, drugs, and toxins. However, recent studies identified new and emerging risk factors for ICC, occupational and environment-related [7, 8]: the chronic viral hepatitis, liver cirrhosis-alcohol-related, smoking, obesity, diabetes and asbestos [9–13]. Patients with unresectable disease (70-90 %) have a poor prognosis with a survival of less than 12 months following diagnosis.
The lack of effective therapies prompts to identify alternative approaches, based on a deepen molecular knowledge. The high throughput techniques, i.e. gene and microRNA profiling, next generation sequencing (NGS), exome sequencing, provide huge amount of data and information suitable to identify potential drug targets [14, 15]. Nowadays, pathogenesis and drug response are usually studied on preclinical models represented by cell lines, primary cultures, and xenografts.
In particular, xenografts and orthotopic models obtained by CCA cell lines, carcinogen-induced and genetically engineered mouse model for CCA has been created . In the last years, patient-derived cancer xenograft (PDX) models have been established by directly engrafting surgically resected human tumor tissues into immune compromised mice. Molecular and genetic analysis demonstrated that PDXs rely primary tumor characteristics, making them suitable models to study pathogenesis and to test anti-cancer drugs activity. PDXs are established from different cancer types, including gastric, breast, ovarian, colon, lung, prostate, and pancreatic cancers [17–23].
To date, no human CCA models derived from tumor patients have been developed.
Here, we established and characterized a patient-derived ICC model derived from a patient of Italian origin. This model will be helpful either to provide a more suitable model for preclinical studies or to test drug efficacy.
Establishment and characterization of patient derived xenograft (PDX)
Tumor samples were obtained from Italian patients subjected to surgical resection for ICC. Biological material was obtained from patients who has signed the informed consent, following institutional review board-approved protocols (“PROFILING Protocol, n° 001-IRCC-00 IIS-10” approved by Comitato Etico Interaziendale of A.O.U. San Luigi Gonzaga, Orbassano, Torino, Italy). This institutional study provides molecular genetic analysis, set up of primary cultures and the creation of PDX from tumor biological samples (primary tumor, metastasis, tumor cells taken under paracentesis or thoracentesis procedures, and blood). We have overall implanted 17 fresh tumor specimens from ICC patients, 14 primary (PR) and 3 recurrent tumors, here named from CHC001 to CHC020.
For PDX establishment, NOD (Non-Obese Diabetic)/Shi-SCID (severe combined immunodeficient) female mice (4–6 weeks old) (Charles River Laboratory) were maintained under sterile conditions in micro-isolator cages at the animal facilities of IRCCS-Candiolo. All animal procedures were approved by the Institutional Ethical Commission for Animal Experimentation (Fondazione Piemontese per la Ricerca sul Cancro) and by the Italian Ministry of Health. Mice were subcutaneously grafted with a fragment of 4x4 mm of representative tumor.
The expression of biliary markers Cytokeratin (CK) 7, 17, 19, and epithelial membrane antigen (EMA)  was evaluated by immunohistochemical analysis (IHC) to compare the characteristics of primary and engrafted tumor. Slides were incubated with primary antibodies followed by the appropriate secondary antibodies; the reaction was visualized by DAB (3,3-diaminiobenzidine) and counterstained with hematoxylin.
Comparative genomic hybridization array
Genomic DNA of PR and its PDX at fourth generation was extracted from formalin fixed, paraffin embedded (FFPE) tissues using the QiAmp FFPE DNA mini Kit (Qiagen). High-resolution oligonucleotide comparative genomic hybridization (CGH) arrays analysis was performed following standard operating procedures of Agilent Technologies. One thousand ng of DNA were digested by a double enzymatic digestion (AluIþRsa I), fragmented, amplified, and purified. After the quantification with Nanodrop, 2 μg of genomic DNA of both tumor and control from Promega (Human Genomic DNA Female N 30742202/male N 30993901) were labeled with CY5-dCTPs and CY3-dCTP, respectively, and hybridized on glass arrays (2 X105 K) at 65C° for 40 hours at 20 rpm. Slides were then washed, scanned on an Agilent 4000C dual laser scanner and images analyzed with Feature Extraction v10.5 software. Raw txt files were then loaded into Cytogenomics software for data processing and visualization.
Gene and microRNA expression analysis
For gene expression analysis (GEP), tissues were homogenized by using TissueLyser LT (Qiagen s.r.l. Milano, Italy) and total RNA (mRNA and microRNA) was extracted and purified by Absolutely RNA miRNA kit (Agilent Technologies), following manufacturers’ protocols. Quantitative and qualitative evaluation of total RNA was performed by Nanodrop and BioAnalyzer, respectively. For GEP analysis, 100 ng of total RNA were amplified and labeled using Low Input Quick Amp Labeling Kit, one-color kit (Agilent Technologies). Six hundred ng of labeled RNA were hybridized on SurePrint G3 Human Gene Expression 8x60K v2 glass arrays. Arrays were scanned and images analyzed by the Feature Extraction Software from Agilent Technologies (version 10.7); raw data were then processed using the Bioconductor package Limma (Linear models for microarray analysis). Background correction was performed with the normexp method with an offset of 50, and quantile was used for the between-array normalization. The empirical Bayes method was used to compute a moderated t-statistics.
For microRNA analysis, 100 ng of total RNA were labeled using the miRNA Complete Labeling and Hyb Kit and hybridized on Human miRNA Microarray Kit Release 16.0, 8x60K. Arrays were scanned and images analyzed by the Feature Extraction Software from Agilent Technologies (version 10.7). Raw data elaboration was carried out with Bioconductor (http://www.bioconductor.org/) , using R statistical language. Background correction was performed with the normexp method, and quantile was used for the between-array normalization. External datasets: GSE26566 and GSE47764 datasets, containing normal bile duct gene and miRNA expression profiles respectively, were downloaded from the GEO website (http://www.ncbi.nlm.nih.gov/geo/). To merge these raw data to our own, we first averaged the signal at probe level (for microRNA arrays, performed on two different versions of Agilent platform) or at gene symbol level (for gene expression arrays, performed on two different platforms). The obtained matrices were then merged and normalized with the quantile function. The LIMMA (LInear Models for Microarray Analysis) package was used to identify differentially expressed genes/microRNAs in tumor versus normal samples. The empirical Bayes method was used to compute a moderated t-statistics .
MicroRNA validation by qRT-PCR
MicroRNA of PDX and of a pool of liver normal tissues was transcribed in cDNA by using TaqMan microRNA Reverse Transcription Kit (Applied Biosystem) using specific primers for mir-21, mir-199, mir-200, mir-31, and for the housekeeping RPL-21. The TaqMan microRNA Assays (with the different fluorescent probes) and the TaqMan Universal MasterMix NO Amperase UNG were used to perform the quantitative Real-time PCR. All the experiments were carried out in triplicate in optical grade 96-well plates. Quantitative analysis was performed by the measurement of Ct values; briefly, to calculate the relative expression of the target microRNA normalized to RPL21, the average of target Ct was subtracted from the average of RPL21 Ct(ΔCt). The amount of target, normalized to an endogenous reference and relative to a calibrator (fold-change) is given by 2-ΔΔCt where the calculation of ΔΔCt involves subtraction by the ΔCt calibrator value (pool of liver normal tissues).
Genomic DNA was extracted by using QIAamp DNA FFPE Mini kit (Qiagen, Milan, Italy) following the manufactures’s instructions. For formalin fixed and paraffin embedded (FFPE) tumor the neoplastic area was obtained by laser microdissection (VSL-337ND-S, Spectra-Physics, Mountain View, CA). The kinase domain of EGFR coding sequence, from exons 18 to 21, was amplified by using primers and nested polymerase chain reaction (PCR) conditions previously described by Lynch and coll . Exons 2 to 4 of K-RAS and N-RAS, exons 9 and 20 of PI3KCA, exon 15 of B-RAF were amplified by PCR as previously described [28, 29]. PCR products were then purified using Wizard® SV Gel and PCR Clean-Up System (Promega, Italy) and sense and antisense sequences were obtained using forward and reverse internal primers, respectively. Each exon was sequenced using the BigDye Terminator Cycle sequence following the PE Applied Biosystem strategy and Applied Biosystem ABI PRISM3100 DNA Sequencer (Applied Biosystem, Forster City, CA). Mutations were confirmed performing two independent PCR amplifications.
Generation and characterization of BTC patient derived xenografts
ICC tumors obtained from surgery were subcutaneously implanted into NOD/SCID mice as described in the Materials and Methods section. Characteristics of tumor patients were summarized in Table 1. Ten patients were females and seven males and the age ranged from 44 to 82; 14 out of 17 (82.4 %) tumor specimens were primary tumors and 3 out of 17 (17.6 %) were recurrences.
Only one tumor out of 17 (5.8 %) was successfully engrafted. It was a primary tumor and was histopathologically classified as pT2b pN0, moderately differentiated (G2) ICC. Tumor sample was also evaluated for the presence of HBV or HCV markers, resulting negative. Patient had chronic colecystitis, but did not have liver cirrhosis or chronic liver disease, primary sclerosing cholangitis diabetes, obesity.
Primary tumor, named CHC001 PDX, was successfully engrafted in mice at the first generation after 4 months; after reaching a volume of 1000 mm3, tumor was explanted and re-implanted in new mice. Starting from the second generation, the latency of growth was decreased from 4 months to 1 month until the stabilization obtained at the fourth generation. If cryopreserved in DMSO 10 % and FBS 90 % in culture medium, it was able to successfully engraft in mice when re-implanted. After stabilization, immunohistochemical and molecular investigations were performed to verify if both features were retained in the PDX.
Immunohistochemistry analysis for the expression of Cytokeratin 7, 17, 19 and EMA as well as the Hematoxillin & Eosin staining  showed that PDX retained the same morphology of PR up to the fourth generation as well as the same immunoreactivity (Fig. 1).
The genomic status of PR and of its PDX was assessed by array CGH technique. As shown in Fig. 2, we found a concordance between the two samples (r = 0.64 by Pearson correlation); the number of common chromosomal alteration was 24 with 7 gained regions and 17 lost regions; the most statistically significant chromosome regions included the loss of the regions in 3p, 5q, 6p, 8p, 9p, 14q, 18q, and the gain of the regions 1p, 2q, 3q, and 12p, 15q, and 20q. Table 2 summarizes the common aberrant regions.
Further, we revealed that PDX acquired other alterations, in particular the loss of 3p, the entire 4, 6q, the entire 7, 10p, 11p, 12q, 15q, 17p, 19p, 21q and 22q, and the gained regions in 5p, 10q, 13q, 15q, and 20q.
To further characterize the PDX model, we selected genes allocated in the aberrant regions typical of PDX; considering the first 500 amplified or deleted genes, respectively, we performed Gene Ontology (GO) analysis, and GO categories are summarized in Additional file 1: Table S1.
Gene expression profiling
Gene expression analysis was performed on the primitive tumor and on the PDX at the fourth passage. Figure 3 showed the correlation plot of differentially expressed genes obtained by Pearson correlation function; this correlation is very high (r = 0.94), enforcing that PDX retained primary tumor characteristics. In order to find the peculiar characteristics of this tumor, common differentially expressed genes were compared to six normal bile duct samples, belonging to the cohort of Andersen and collaborators . Genes list was filtered on adjusted p-value (<0,00001) and the most significant 300 probes were analyzed for Gene Ontology; we found that down-regulated genes are involved in blood coagulation, inflammation response, and in lipid metabolism; on the contrary, up-regulated genes globally affected DNA biosynthesis processes, as nucleosome assembly and organization, translation, underlying that tumor cells are more active rather than normal cells. Even the high correlation of gene expression data, we found 63 up-regulated and 276 down-regulated genes altered in PDX versus primary tumor (Additional file 2: Table S2). Further, we compared differentially expressed genes in PDX with the list of genes allocated in amplified or deleted regions found in PDX; 5 up-regulated and 32 down-regulated genes were found to be overlapped (Additional file 3: Table S3).
MicroRNA expression profiling
The comparison between PR and its PDX revealed a high correlation in terms of microRNA expression (r = 0.92 by Pearson correlation), as shown in Fig. 4. Common deregulated microRNAs were compared with those obtained by normal bile duct in a work of Peng et al. . Row data were filtered with a logFC < or > 0.58 and a p-value of < 0.01. An unsupervised hierarchical cluster showed the deregulated microRNAs among primary and PDX tumors compared to normal bile duct (Fig. 5). Twenty-eight microRNAs (Table 3), of which 7 down-regulated and 21 up-regulated were selected. Nine out of 28 microRNAs are involved in the negative or positive regulation of cell cycle, apoptosis, migration and proliferation, underlying that these processes are altered in tumor cells. In order to enforce these data, we validated the expression of 4 microRNAs by qRT-PCR. As shown in Additional file 4: Table S4, the trend of expression of mir-21, mir-200, mir-199, and mir-31 is confirmed.
Furthermore, we analyzed if PDX acquired peculiar characteristics in terms of microRNA expression; Additional file 5: Table S5 showed that only let-7a-5p, miR-15b-5p, let-7d-5p, miR-200b-5p were down-regulated in PDX compared to primary tumor.
Mutational analysis of the kinase domain of EGFR coding sequence, from exons 18 to 21, exons 2,3 and 4 of K-RAS, exons 2,3 and 4 of N-RAS, exons 9 and 20 of PI3KCA, and exon 15 of B-RAF were performed on PR and on PDX. As shown in Fig. 6, only the sequence of K-RAS exon 2 is mutated (G12D mutation) in the primary tumor (panel B) and is maintained in PDX (panel C).
Intrahepatic cholangiocarcinoma constitutes the second most common primary hepatic malignancy with a very poor prognosis [32, 33]. Thus, the identification of alternative therapeutic options is an urgent step to improve the outcome of these patients. ICC PDX models could represent an useful tool either to study the disease from biological and molecular aspects or to investigate response to new therapies. Here, we established and characterized, for the first time, an Italian ICC PDX derived from fresh tumor tissue.
We subcutaneously implanted 17 ICC fresh tumor tissues into immunocompromised mice and we obtained a rate of successful engraftment of 5.8 % (1/17). The engraftment was reached after 4 months from implant, while for the subsequent generations the latency was significantly reduced to one month. The same result was obtained re-implanting archival frozen tissues.
The limited success of engraftment of these tumors is not clear. For colorectal cancer PDX the engraftment rate is 67 % ; as concerning mammary tumors, the rate is higher with metastatic tissues rather than primary tumors; moreover, graft achievement depends on other factors, as tumor histotypes, grading, and on the presence of Estrogen and HER2 receptors . We can speculate that the presence of K-RAS mutation in our PDX model could be a driver of the more aggressive phenotype, thus explaining the successful engraftment, as shown in colorectal cancer PDX model .
K-RAS mutations are one of the biological determinants of anti-EGFR target therapy resistance in colorectal cancer . Although the role of K-RAS in response to the anti-EGFR therapy in CCA is controversial [37–39], this model could be suitable for the evaluation of the effectiveness of alternative therapies in K-RAS mutated patients for whom anti- EGFR therapies are unfit.
We further compared immunophenotypical and molecular features of PR with its corresponding PDX and we found that both tissue architecture and immunoreactivity of biliary epithelial markers were maintained in PDX.
The genetic relationship between PR and its PDX was established by array-CGH analysis; some genetic alterations were found and maintained from PR to PDX. Some of these regions, in particular the loss of 3p, 6p, 8p, 9p, and 14q regions and the gain of 3q, and 20q regions, are common with the previously found by Miller and coll. in an ICC case series . Other genetic alterations were found only in PDX; this could suggest that, even if the PDX retained the main characteristics of primary tumor, i) the murine environment leads to the acquirement of further chromosomic alterations, ii) the tumor experiences progression regardless of recipient and acquires new chromosomal aberrations, as previously demonstrated by Shiraishi and collaborators , iii) the more aggressive cell subpopulation is selected in the murine model.
Comparing transcriptomic profiling of primitive tumor and its PDX, we found a high correlation in terms of gene and microRNA expression, demonstrating that the PDX retained most of primary tumor genetic characteristics. To further characterize our model, we identified a panel of deregulated genes comparing both PR and PDX tumor with published normal bile duct epithelia; to overcome the lack of normal samples in our Institution, we used external dataset of normal biliary tissues, even introducing a possible bias. We select a panel of down-regulated genes involved in blood coagulation, inflammation response, and in lipid metabolism processes and a panel of up-regulated genes involved in DNA biosynthesis processes. We also selected a panel of twenty-eight microRNAs (7 down-regulated and 21 up-regulated), most of them involved in the cell cycle, apoptosis, migration and proliferation regulation. In particular, we found an up-regulation of miR-21, already described in CCA; in fact, the mir-21 overexpression is typical of ICC compared to both normal tissues and hepatic cancer [42, 43]. Moreover, functional studies on CCA cell lines showed the potential oncogenic role of miR-21 by inhibiting PDCD4 and TIMP3, involved in apoptosis and in the inhibition of the matrix metalloproteinases, respectively .
In conclusion, in this study we firstly established an ICC PDX model and characterized it for genetic and molecular alterations; we demonstrated that this model recapitulates the histological characteristics and maintains most of the genetic features of primary tumor, providing a reliable tool to study this neoplasia and to test the efficacy of new drug.
Next generation sequencing
Patient-derived cancer xenograft
Epithelial membrane antigen
Comparative genomic hybridization array
(Formalin fixed, paraffin embedded)
Gene expression profiling
fetal bovine serum
Yamasaki S. Intrahepatic cholangiocarcinoma: macroscopic type and stage classification. J Hepatobiliary Pancreat Surg. 2003;10(4):288–91.
Welzel TM, McGlynn KA, Hsing AW, O’Brien TR, Pfeiffer RM. Impact of classification of hilar cholangiocarcinomas (Klatskin tumors) on the incidence of intra- and extrahepatic cholangiocarcinoma in the United States. J Natl Cancer Inst. 2006;98(12):873–5.
Hammill CW, Wong LL. Intrahepatic cholangiocarcinoma: a malignancy of increasing importance. J Am Coll Surg. 2008;207(4):594–603.
Patel T. Increasing incidence and mortality of primary intrahepatic cholangiocarcinoma in the United States. Hepatology. 2001;33(6):1353–7.
Dodson RM, Weiss MJ, Cosgrove D, Herman JM, Kamel I, Anders R, Geschwind JF, Pawlik TM. Intrahepatic cholangiocarcinoma: management options and emerging therapies. J Am Coll Surg. 2013;217(4):736–50. e734.
Tsai S, Nathan H, Pawlik TM. Primary liver cancer: intrahepatic cholangiocarcinoma emerges from the shadows. Updates Surg. 2010;62(1):5–9.
Kumagai S, Kurumatani N, Arimoto A, Ichihara G. Cholangiocarcinoma among offset colour proof-printing workers exposed to 1,2-dichloropropane and/or dichloromethane. Occup Environ Med. 2013;70(7):508–10.
Brandi G, Di Girolamo S, Farioli A, de Rosa F, Curti S, Pinna AD, Ercolani G, Violante FS, Biasco G, Mattioli S. Asbestos: a hidden player behind the cholangiocarcinoma increase? Findings from a case-control analysis. Cancer Causes Control. 2013;24(5):911–8.
Tyson GL, El-Serag HB. Risk factors for cholangiocarcinoma. Hepatology. 2011;54(1):173–84.
Shin HR, Lee CU, Park HJ, Seol SY, Chung JM, Choi HC, Ahn YO, Shigemastu T. Hepatitis B and C virus, Clonorchis sinensis for the risk of liver cancer: a case-control study in Pusan, Korea. Int J Epidemiol. 1996;25(5):933–40.
Donato F, Gelatti U, Tagger A, Favret M, Ribero ML, Callea F, Martelli C, Savio A, Trevisi P, Nardi G. Intrahepatic cholangiocarcinoma and hepatitis C and B virus infection, alcohol intake, and hepatolithiasis: a case-control study in Italy. Cancer Causes Control. 2001;12(10):959–64.
Zhou YM, Yin ZF, Yang JM, Li B, Shao WY, Xu F, Wang YL, Li DQ. Risk factors for intrahepatic cholangiocarcinoma: a case-control study in China. World J Gastroenterol. 2008;14(4):632–5.
Palmer WC, Patel T. Are common factors involved in the pathogenesis of primary liver cancers? A meta-analysis of risk factors for intrahepatic cholangiocarcinoma. J Hepatol. 2012;57(1):69–76.
Andersen JB, Spee B, Blechacz BR, Avital I, Komuta M, Barbour A, Conner EA, Gillen MC, Roskams T, Roberts LR, et al. Genomic and genetic characterization of cholangiocarcinoma identifies therapeutic targets for tyrosine kinase inhibitors. Gastroenterology. 2012;142(4):1021–31. e1015.
Ross JS, Wang K, Gay L, Al-Rohil R, Rand JV, Jones DM, Lee HJ, Sheehan CE, Otto GA, Palmer G, et al. New routes to targeted therapy of intrahepatic cholangiocarcinomas revealed by next-generation sequencing. Oncologist. 2014;19(3):235–42.
Ko KS, Peng J, Yang H. Animal models of cholangiocarcinoma. Curr Opin Gastroenterol. 2013;29(3):312–8.
Zhang H, Cohen AL, Krishnakumar S, Wapnir IL, Veeriah S, Deng G, Coram MA, Piskun CM, Longacre TA, Herrler M, et al. Patient-derived xenografts of triple-negative breast cancer reproduce molecular features of patient tumors and respond to mTOR inhibition. Breast Cancer Res. 2014;16(2):R36.
Scott CL, Becker MA, Haluska P, Samimi G. Patient-derived xenograft models to improve targeted therapy in epithelial ovarian cancer treatment. Front Oncol. 2013;3:295.
Zhang T, Zhang L, Fan S, Zhang M, Fu H, Liu Y, Yin X, Chen H, Xie L, Zhang J, et al. Patient-Derived Gastric Carcinoma Xenograft Mouse Models Faithfully Represent Human Tumor Molecular Diversity. PLoS One. 2015;10(7):e0134493.
Ilie M, Nunes M, Blot L, Hofman V, Long-Mira E, Butori C, Selva E, Merino-Trigo A, Vénissac N, Mouroux J, et al. Setting up a wide panel of patient-derived tumor xenografts of non-small cell lung cancer by improving the preanalytical steps. Cancer Med. 2015;4(2):201–11.
Cho YB, Hong HK, Choi YL, Oh E, Joo KM, Jin J, Nam DH, Ko YH, Lee WY, et al. Colorectal cancer patient-derived xenografted tumors maintain characteristic features of the original tumors. J Surg Res. 2014;187(2):502–9.
Lin D, Wyatt AW, Xue H, Wang Y, Dong X, Haegert A, Wu R, Brahmbhatt S, Mo F, Jong L, et al. High fidelity patient-derived xenografts for accelerating prostate cancer discovery and drug development. Cancer Res. 2014;74(4):1272–83.
Mattie M, Christensen A, Chang MS, Yeh W, Said S, Shostak Y, Capo L, Verlinsky A, An Z, Joseph I, et al. Molecular characterization of patient-derived human pancreatic tumor xenograft models for preclinical and translational development of cancer therapeutics. Neoplasia. 2013;15(10):1138–50.
Akiba J, Nakashima O, Hattori S, Tanikawa K, Takenaka M, Nakayama M, Kondo R, Nomura Y, Koura K, Ueda K, et al. Clinicopathologic analysis of combined hepatocellular-cholangiocarcinoma according to the latest WHO classification. Am J Surg Pathol. 2013;37(4):496–505.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 2004, 3:Article3.
Lynch TJ, Bell DW, Sordella R, Gurubhagavatula S, Okimoto RA, Brannigan BW, Harris PL, Haserlat SM, Supko JG, Haluska FG, et al. Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. N Engl J Med. 2004;350(21):2129–39.
Pignochino Y, Sarotto I, Peraldo-Neia C, Penachioni JY, Cavalloni G, Migliardi G, Casorzo L, Chiorino G, Risio M, Bardelli A, et al. Targeting EGFR/HER2 pathways enhances the antiproliferative effect of gemcitabine in biliary tract and gallbladder carcinomas. BMC Cancer. 2010;10:631.
Negru S, Papadopoulou E, Apessos A, Stanculeanu DL, Ciuleanu E, Volovat C, Croitoru A, Kakolyris S, Aravantinos G, Ziras N, et al. KRAS, NRAS and BRAF mutations in Greek and Romanian patients with colorectal cancer: a cohort study. BMJ Open. 2014;4(5):e004652.
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. Epub ahead of print.
Peng F, Jiang J, Yu Y, Tian R, Guo X, Li X, Shen M, Xu M, Zhu F, Shi C, et al. Direct targeting of SUZ12/ROCK2 by miR-200b/c inhibits cholangiocarcinoma tumourigenesis and metastasis. Br J Cancer. 2013;109(12):3092–104.
Zhang GW, Lin JH, Qian JP, Zhou J. Identification of risk and prognostic factors for patients with clonorchiasis-associated intrahepatic cholangiocarcinoma. Ann Surg Oncol. 2014;21(11):3628–37.
Siegel R, Ma J, Zou Z, Jemal A. Cancer statistics, 2014. CA Cancer J Clin. 2014;64(1):9–29.
Marangoni E, Vincent-Salomon A, Auger N, Degeorges A, Assayag F, de Cremoux P, de Plater L, Guyader C, De Pinieux G, Judde JG, et al. A new model of patient tumor-derived breast cancer xenografts for preclinical assays. Clin Cancer Res. 2007;13(13):3989–98.
Bertotti A, Migliardi G, Galimi F, Sassi F, Torti D, Isella C, Corà D, Di Nicolantonio F, Buscarino M, Petti C, et al. A molecularly annotated platform of patient-derived xenografts (“xenopatients”) identifies HER2 as an effective therapeutic target in cetuximab-resistant colorectal cancer. Cancer Discov. 2011;1(6):508–23.
Misale S, Yaeger R, Hobor S, Scala E, Janakiraman M, Liska D, Valtorta E, Schiavo R, Buscarino M, Siravegna G, et al. Emergence of KRAS mutations and acquired resistance to anti-EGFR therapy in colorectal cancer. Nature. 2012;486(7404):532–6.
Chen JS, Hsu C, Chiang NJ, Tsai CS, Tsou HH, Huang SF, Bai LY, Chang IC, Shiah HS, Ho CL, et al. A KRAS mutation status-stratified randomized phase II trial of gemcitabine and oxaliplatin alone or in combination with cetuximab in advanced biliary tract cancer. Ann Oncol. 2015;26(5):943–9.
Malka D, Cervera P, Foulon S, Trarbach T, de la Fouchardière C, Boucher E, Fartoux L, Faivre S, Blanc JF, 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(8):819–28.
Borbath I, Ceratti A, Verslype C, Demols A, Delaunoit T, Laurent S, Deleporte A, Vergauwe P, Van Maanen A, Sempoux C, et al. Combination of gemcitabine and cetuximab in patients with advanced cholangiocarcinoma: a phase II study of the Belgian Group of Digestive Oncology. Ann Oncol. 2013;24(11):2824–9.
Miller G, Socci ND, Dhall D, D’Angelica M, DeMatteo RP, Allen PJ, Singh B, Fong Y, Blumgart LH, Klimstra DS, et al. Genome wide analysis and clinical correlation of chromosomal and transcriptional mutations in cancers of the biliary tract. J Exp Clin Cancer Res. 2009;28:62.
Shiraishi K, Okita K, Harada T, Kusano N, Furui T, Kondoh S, Oga A, Kawauchi S, Fukumoto Y, Sasaki K, et al. Comparative genomic hybridization analysis of genetic aberrations associated with development and progression of biliary tract carcinomas. Cancer. 2001;91(3):570–7.
Karakatsanis A, Papaconstantinou I, Gazouli M, Lyberopoulou A, Polymeneas G, Voros D. Expression of microRNAs, miR-21, miR-31, miR-122, miR-145, miR-146a, miR-200c, miR-221, miR-222, and miR-223 in patients with hepatocellular carcinoma or intrahepatic cholangiocarcinoma and its prognostic significance. Mol Carcinog. 2013;52(4):297–303.
Chen L, Yan HX, Yang W, Hu L, Yu LX, Liu Q, Li L, Huang DD, Ding J, Shen F, et al. The role of microRNA expression pattern in human intrahepatic cholangiocarcinoma. J Hepatol. 2009;50(2):358–69.
Selaru FM, Olaru AV, Kan T, David S, Cheng Y, Mori Y, Yang J, Paun B, Jin Z, Agarwal R, et al. MicroRNA-21 is overexpressed in human cholangiocarcinoma and regulates programmed cell death 4 and tissue inhibitor of metalloproteinase 3. Hepatology. 2009;49(5):1595–601.
This work was supported by grant from “Associazione Italiana Ricerca sul Cancro–AIRC 5X1000 2010-Ministry of Health, FPO. Project n°16:30 “Identificazione di nuove vie di trasduzione del segnale intracellulare sensibili ai farmaci nel colangiocarcinoma intraepatico (ICC)”. Fondazione Piemontese per la Ricerca sul Cancro - Onlus - “Identification of new druggable pathways in intrahepatic cholangiocarcinoma” 5 per Mille 2010 Ministero della Salute.
Authors declare that they have no competing interests.
GC and CPN (co-authors): designed the study, conceived and carried out the experiments, FS: carried out IHC experiments; GC: interpreted microarray data; IS: performed IHC analysis;, FL and MA: designed the study, wrote and revised the article. All authors have read and approved the manuscript.
Gene Ontology categories (biological processes, cellular components, molecular functions, and pathways) altered in PDX. (XLSX 15 kb)
Differentially expressed genes in PDX versus primary tumor (PR). (XLSX 22 kb)
Overlapped deregulated genes obtained comparing differentially expressed genes (GEP) with the list of genes allocated in amplified or deleted regions (aCGH) found in PDX. (XLSX 13 kb)
Relative quantitation of deregulated microRNA using the comparative Ct method. (DOCX 18 kb)
Deregulated microRNA in PDX compared to primary tumor (PR). (XLSX 9 kb)