- Open Access
A novel patient-derived organoids-based xenografts model for preclinical drug response testing in patients with colorectal liver metastases
Journal of Translational Medicine volume 18, Article number: 234 (2020)
Cancer-related mortality in patients with colorectal cancer (CRC) is predominantly caused by development of colorectal liver metastases (CLMs). How to screen the sensitive chemotherapy and targeted therapy is the key element to improve the prognosis of CLMs patients. The study aims to develop patient-derived organoids-based xenografted liver metastases (PDOX-LM) model of CRC, to recapitulate the clinical drug response.
We transplanted human CRC primary tumor derived organoids in murine spleen to obtain xenografted liver metastases in murine liver. Immunohistochemistry (IHC) staining, whole-exome and RNA sequencing, and drug response testing were utilized to identify the homogeneity in biological and genetic characteristics, and drug response between the PDOX-LM models and donor liver metastases.
We successfully established PDOX-LM models from patients with CLMs. IHC staining showed that positive expression of CEA, Ki67, VEGF, FGFR2 in donor liver metastases were also well preserved in matched xenografted liver metastases. Whole-exon sequencing and transcriptome analysis showed that both xenografted and donor liver metastases were highly concordant in somatic variants (≥ 0.90 frequency of concordance) and co-expression of driver genes (Pearson’s correlation coefficient reach up to 0.99, P = 0.001). Furthermore, drug response testing showed that the PDOX-LM models can closely recapitulated the clinical response to mFOLFOX6 regiments.
This PDOX-LM model provides a more convenient and informative platform for preclinical testing of individual tumors by retaining the histologic and genetic features of donor liver metastases. This technology holds great promise to predict treatment sensitivity for patients with CLMs undergoing chemotherapy.
Colorectal cancer (CRC) represents the third leading cause of cancer-related deaths worldwide . Cancer-related mortality is predominantly caused by development of colorectal liver metastases (CLMs), which occurs in approximately 50% of all patients during the course of their illness . Despite application of target agents and advanced therapeutic schedule, the median survival of patients with metastatic CRC (mCRC) is limited at approximately 30 months . Therefore, how to screen the sensitive drugs for CLMs is critical to improve survival of these patients.
The pre-clinical animal models play critical roles in precision medicine, by elucidating biomarkers that predict drug response and identify patients who are most likely to benefit from a specific treatment [4, 5]. An ideal animal model should recapitulate the molecular, histopathological, and etiological characteristics of the donor tumors . Currently, a number of animal models have evolved toward sophisticated systems that mimic the pathology of CRC. In particular, patient-derived xenografts (PDX) models and genetically engineered mouse models (GEMMs) have been widely used in metastases research [4, 7]. On one hand, PDX models have advantages in replicating the characteristics and genetic diversity of the donor tumor, and these models provide platform for pre-clinical drug response testing. However, the metastatic PDX models are time-consuming, expensive, and technically challenging . On the other hand, the GEMMs are popular in the research of carcinogenic progress and the mechanisms of specific cancer-related genes. However, they usually cannot fully reproduce the genetic complexity of human tumors in the clinical practice.
Recently, the emerging organoid culture technology allows for cancer cells to be grown in a 3-dimensional matrix and have been a critical progress in cancer research . Prior studies have demonstrated that patient-derived cancer organoids can be easily maintained and manipulated in vitro and faithfully recapitulate characteristics of the patient tissues [9, 10]. The establishment of living tumor organoid biobanks offers a platform for high-throughput drug screens . Although the patient-derived organoids (PDOs) represents a powerful resource for finding effective therapeutic strategies direct to specific tumor subtypes, it does not account for the interplay between tumor cells and the surrounding tissue microenvironment, because this interplay could not be recapitulated in a dish [9,10,11]. Many studies have pointed to the importance of tumor microenvironment in influencing tumor cell identity and behavior [12,13,14,15], emphasizing the necessity of validating the results obtained in vitro in animal model systems.
To overcome these limitations, we have recently developed a method to transplant human CRC organoids in murine spleen to obtain xenografted liver metastases model, named patient-derived organoids-based xenografted liver metastases (PDOX-LM) model. This transplantation approach allows primary tumor formation in murine spleen and the spontaneous development of metastases in murine liver. And we further verified the homogeneity between the xenografted and donor liver metastases in histological feature, tumor-specific protein, genetic and transcriptome characterization and drug response.
Colonic tissues were obtained from the Zhongshan Hospital with informed consent and the study was approved by the ethical committee. All patients were diagnosed with colorectal cancer. From the resected colorectal segment, tumor tissue was isolated. The isolation of tumor epithelium was performed essentially as described by Hans Clevers .
Organoid preparation and transplantation
Human tissue-derived organoids were obtained and cultured according to Hans Clevers . The expanded organoids were harvested with Cell Recovery Solution (BD Biosciences) and suspended in Matrigel at a concentration of 1 × 105 cells/ml. A total 200 μl (2 × 104 cells) of Matrigel-organoid suspension was injected into the spleen of each Balb/c-nu mice. The mice were euthanized 7 weeks after xenotransplantation. The spleens, the livers, the lungs and the other abdominal organs were isolated and examined. The size of the engrafted or metastasized tumors were determined on measurement by vernier caliper. Each CRC organoid line was transplanted into four independent Balb/c-nu mice. Some samples were then stored in liquid nitrogen, and the rest of samples were formalin-fixed and embedded in paraffin for subsequent IHC analysis. All animal procedures were approved by the Animal Ethics Committee of Zhongshan hospital and complied with the Guide for the Care and Use of Laboratory Animals.
Whole-exome sequencing and single nucleotide polymorphisms
For each sample, 250 ng DNA was sheared and subject to whole-exome sequencing using the Agilent v2 capture probe set and sequenced by HiSeq 2500 using 76 base pair reads, as previously described [17, 18]. A median 12.4 Gb of unique sequence was generated for each sample. Sequence data were locally realigned to improve sensitivity and reduce alignment artifacts prior to identification of mutations, insertions, and deletions as previously described [19,20,21].
Somatic copy-number analysis
Somatic copy-number analysis was performed using segmented copy-number profiles generated from whole-exome sequencing using the SegSeq algorithm . Briefly, read depth in tumor and normal pairs was calculated to provide relative copy-number ratios at each exon followed by circular binary segmentation . The GISTIC2.0 algorithm  was used to investigate focal and arm-level copy-number changes (Additional file 2: Table S1A–E).
RNA sequence data processing
RNA was extracted from xenografts and tumor tissues using QIAGEN RNA mini kit. RNA sample was hybridized on Affymetrix Human Gene 2.0 ST arrays. The raw CEL files were processed with Affymetrix Power Tools using the Hg19 genome build and NetAffx annotation dating from 09-30-2018. Between-array normalization was performed using rma-sketch, within APT. This resulted in an intensity matrix of 21,681 genes by 2 samples. For analysis of individual genes, data were analyzed using the R2 web application, which is freely available at http://r2.amc.nl.
Drug response testing
PDOX-LM models based on Patient #3, #4, #5 and #6 were randomly distributed into two groups: the control group (n = 16) for intraperitoneal injection with phosphate buffer saline [PBS, 10 mM, 0.1 ml per mouse] 2 times every weeks, and the chemotherapy group (n = 16) for intraperitoneal injection with the combination of oxaliplatin (0.4 mg/ml PBS) and 5-Fu (4 mg/ml PBS) [0.1 ml per mouse] 2 times every week. All intraperitoneal injection started from the 7th weeks since the splenic xenotransplantation of mCRC PDOs, and last for 4 weeks. In addition, the mice in chemotherapy group were further divided into PR group and PD group according to the clinical drug response of the donor patients. Supplementary methods were detail showed in the Additional file 3, and the primary antibodies used for IF or IHC staining were showed in the Table S2.
Protocol of establishment and analysis of a PDOX-LM model
Here we presented a detailed description of the procedure for the PDOX-LM model (Additional file 1: Figure S1). This approach allowed the formation of xenografted metastases in murine liver, but not in extra-hepatic sites. By immunohistochemical (IHC) staining, whole-exome and RNA sequence, and drug response testing, we verified that the xenografted liver metastases formed by organoids of primary tumor recapitulated similar mutations, gene expression and drug response of the donor liver metastases.
Establishment of patient-derived organoids from primary tumors of CRC patients
We adopted the Wnt-dependency of cancer stem cell to selectively expand tumor organoids, according to the previous reports [10, 16, 25]. In this study, surgically resected tissues were obtained from previously untreated CRC patients. Primary tumor derived organoids were all successfully generated from the samples of 6 metastatic mCRC and 6 localized CRC cases (Table 1). To develop a long-term expandable organoid culture, the combination of epithelial growth factor (EGF) and the Wnt amplifier R-spondin1 was essential to add in the culture medium [8, 14, 15].
We showed the number of CRC organoids varied among different patient samples, with some tumors rendering thousands of organoids whereas others yielding only 10–20 primary organoids, which was similar with the previous studies . Notably, the tumor organoids derived from mCRC cases demonstrated a higher proliferating ability (at least 200 organoids or 1 × 105 cells/ml Matrigel) than those derived from localized CRC cases (a median of 40 organoids or 5 × 104 cells/ml Matrigel) (Fig. 1a).
And we also identified the histologic features of CRC PDOs by hematoxylin and eosin (H&E) staining and characteristic biomarkers by immunofluorescence (IF) staining (Fig. 1b). In all cases, the cultured organoids were similar to the matched donor tumors in secondary architecture, nuclear pleomorphism, nuclear to cytoplasmic ratio, presence of prominent nucleoli, and mitotic rate. The CRC organoids commonly developed crypt-like structures reminiscent of the malignant glands seen within the patient’s cancer (Fig. 1a). In addition, these key phenotypic characteristics, such as lgr5 and E-cadherin, are maintained across passages. According to previous researches, Lgr5, an intestinal stem cell-specific membrane protein, confirmed the stemness of organoids , and positive expression of E-cadherin confirmed the epithelial source of organoids (Fig. 1b).
Establishment of PDOX-LM models
After cutting a freshly procured primary CRC tumor specimen into small tissue blocks (~ 8 mm3), we obtained millions of single cells and some cell clusters by physical and chemical separation. Then we cultured the tumor organoids in vitro according to the protocols. After expanding, the organoids (1 × 105 cells/ml Matrigel) of low-passage (≤ 6 passage) were engrafted into the murine spleens. In brief, we cut through the abdominal walls along the costal margin, and immediately injected 200 μl Matrigel-organoids suspension into the bottom of murine spleen. Following two factors may account for the successful establishment of PDOX-LM models: successful expanding organoids in vitro and immediate stanch of copious bleeding from murine spleen.
The CRC organoids derived from the primary tumors of 6 mCRC and 6 localized CRC cases were transplanted into the spleen of balb/c-nu male mice. Interestingly, we observed that the organoids of mCRC source other than localized CRC source successfully formed macrometastatic colonies (≥ 1 mm in size) in the livers . In addition, using this transplanting method, macrometastatic lesion could only be found in murine livers, but not in extra-hepatic sites (Fig. 1c).
In terms of the dynamic monitoring marker in CRC, serum Carcinoembryonic antigen (CEA) is often used as a typical tumor marker . The in vivo study showed that, after transplantation, we monitored human CEA levels in murine blood dynamically. Human CEA were detectable in murine blood 1 week after transplantation, with concentrations progressively increasing over the following weeks (Fig. 1d). We observed that the serum CEA in PDOX-LM models of mCRC source was higher than that in models of localized CRC source (Fig. 1e). And the serum CEA levels in successfully established PDOX-LM models were higher than that of the mouse models without xenografted liver metastases (Fig. 1f). The serum CEA levels were comparable to those detected in human cancer patients, reaching up to 104 ng/ml (Table 2).
In addition, to explore the index to monitor the establishment of PDOX-LM models, we monitored the body weight of mice after transplantation dynamically. We observed that 8 mice (8/24, 33.3%) had a weight loss at 5th week, and 21 mice (21/24, 87.5%) had a weight loss at 6th week (Fig. 1g). All mice with weight loss were found to bear xenografted liver metastases when sacrificed at 7th week.
All mice were sacrificed in 7th week since transplantation in murine spleens. For transplanted organoids from Patient #1, at most two macrometastatic colonies (≥ 1 mm in size) were formed in the livers, whereas transplanted organoids from Patient #2 formed at least 4 colonies. And for transplanted organoids from Patient #3, #4, #5 and #6, the numbers of macrometastatic colonies was 2–5 in the livers (Table 2). The metastatic capacity of transplanting organoids seemed like substantially diverse.
Histological feature and tumor-specific protein of xenografted and donor liver metastases
Histological features of the parental tumors were well preserved in these macrometastatic lesions, especially in secondary architecture, nuclear pleomorphism, nuclear to cytoplasmic ratio, presence of prominent nucleoli, and mitotic rate (Fig. 2). By means of IHC staining, positive staining of CEA, Ki67, VEGF, FGFR2 and negative staining of EGFR, CK7 of the parental lesions were also well preserved in matched organoids and xenografted liver metastases (Fig. 2). Notably, the xenografted and donor liver metastases showed similar patterns of expression of these CRC tumor-specific markers.
Genomic characterization of xenografted and donor liver metastases
To demonstrate the maintenance of molecular features between the xenografted and donor liver metastases, we performed DNA copy-number analyses and whole-exome sequencing platforms between the xenografted and donor liver metastases [10, 16]. Matched xenografted and donor liver metastases were sequenced for 2 mCRC patients accepting simultaneous resection of primary tumors and liver metastases. And the sequencing depth were greater than 500 × median coverage across the interrogated genes.
Genomic DNA were isolated from liver metastases and matched normal colonic mucosa for sequencing in order to identify tumor-specific somatic mutations . Genomic DNA were available for comparative analysis only for patient #1 and #2. The mutation rates per Mb varied widely for different tumor biopsy specimen (range 3.2 to 8.6), with a median value of 6.4 in the xenografted liver metastases, higher than the median rate of 3.7 in the donor liver metastases (Fig. 3a, Additional file 2: Table S1A–C). Mutations were predominantly CpG to T transitions, consistent with results from large-scale CRC sequencing [9, 16, 25]. Both xenografted liver metastases displayed non-hypermutation (≤ 10 mutations/Mb), which was consistent with the matched donor liver metastases.
Mutations known to enhance the selective growth advantage of the cancer cells, or driver mutations, were specifically examined. These driver mutations were overwhelmingly present in the xenografted liver metastases. Somatic variants within the coding regions in xenografted liver metastases were highly concordant with the corresponding donor liver metastases (Patient #1, 0.96 frequency of concordance; Patient #2, 0.90 frequency of concordance) (Fig. 3b, Additional file 2: Table S1C). Indeed, combined analysis of SCNAs and single nucleotide variants (SNVs) to infer Cancer Cell Fractions (CCF) [28, 29] in the donor and xenografted liver metastases, revealed that the common CRC driver mutations were maintained in the establishment of model (Additional file 2: Table S1D).
We showed that both xenografted and matched donor liver metastases shared many missense growth promoting mutation, such as TP53 (exon3,G > A), POLE (exon8,C > A), AXIN2 (exon 7 G > C; exon 2 G > A), TCF7L2 (exon13 C > A) and RNF43 (exon2 C > T) alterations, suggesting they arose from the same somatically altered progenitor cell (Additional file 2: Table S1B). In brief, the most commonly altered genes in donor biopsy [25, 30, 31] were well represented in the xenografted liver metastases (Fig. 3b). Inactivating alterations to the tumor suppressors APC, TP53, FBXW7 and SMAD4, as well as activating mutations in KRAS (codon 12) were observed (Additional file 2: Table S1B) . Interestingly, the successful establishment of both KRAS mutant type and KRAS wild type organoids implied that KRAS status did not influence the establishment of PDOX-LM models. On the other hand, discordant mutations were also assessed for their biological significance in cancer, based on Cancer Gene Census and data reported from the PanCancer analysis of 5000 whole exomes [31, 32]. Only 3.9% (14/356) of discordant mutations found in xenografted liver metastases affected cancer-associated genes, including a third hit to APC, which was already biallelically inactivated in Patient #1 samples, and MEGF6 mutation in Patient #2 (Fig. 3b, c, Additional file 2: Table S1E). The discordant mutations had a mean allelic frequency of 7.4% and 10.4% for the donor and xenografted biopsy, respectively.
Mutations of genes in DNA mismatch repair (MMR)-associated pathways has been reported to associate with a hypermutated phenotype . We showed that missense mutations were absent in MSH3 in both cases, consistent with their clinical classification as non-hypermutated CRC cases . However, the limited sample size did not allow a statistical analysis for somatic copy number alterations to identify significant regions of amplification and deletions.
In aggregate, our xenografted liver metastases faithfully captured the genomic features of the donor liver metastases.
Transcriptome analysis of donor and xenografted liver metastases
The murine livers provide the similar mesenchyme, blood vessels, hepatic stellate cell, etc. And donor and xenografted liver metastases were formed under closely resemble condition. Therefore, we assumed that this model allows for direct gene expression analysis of subclone with liver affinity. Transcriptomic characterization was analyzed using Affymetrix single transcript arrays. Figure 4a, b showed the volcano plots of differentially expressed genes. Both donor and xenografted liver metastases exhibited high heterogeneity. Top 10 up-regulated genes and down-regulated genes in patient #1 and #2 were also listed (Fig. 4a, b). The most of up- or down-regulated genes were not common CRC driver genes except NOTCH3 and MYH11. Collagen-associated genes, including COL1A1, COL5A1, COL18A1 and COL9A3 were examined in top 10 up- or down-regulated genes (Fig. 4a, b), which attributed to less stoma composition in xenografted liver metastases than that of matched donor. Next, we searched for genes co-expressed between donor and xenografted liver metastases. Both donor and xenografted liver metastases co-expressed driver genes including TP53, KRAS, NRAS, EGFR, ERBB2, APC, CTNNB1, AXIN2, SMAD4 and CDK2 (fold change < 4) (Table 3) . Indeed, large sample analysis are warranted to confirm whether xenografted liver metastases faithfully capture the expression of CRC driver gene in donor liver metastases.
Drug response-associated mutations and gene expression
To identify the accordance of drug response between donor and xenografted liver metastases, we compared drug response-associated mutations and gene expression. To compare the predictive value in drug response, target genes for 30 targeted agents or chemotherapy agents were available for both xenografted and matched liver metastases. The analysis included 20 genes identified as mutant, or frameshift or non-frameshift insertion as described by previous reports (Fig. 5a) [35, 36]. KRAS and NRAS status were detected the same between donor and xenografted liver metastases by WES, which is associated with the response of the anti-EGFR inhibitors cetuximab and BIBW2992 (afatinib) (Figs. 5a). The donor and xenografted liver metastases with mutant-type KRAS demonstrated a higher expression than the counterpart with wild-type KRAS (Figs. 5a, b, Additional file 2: Table S1B).
We also identified a number of discordant mutant-type target genes during the establishment of xenograft model, including IGF2R, EGFR, POLE, AXIN2, APC, CDX2, ERBB3, TP53. Loss-of-function mutations of the tumor suppressor TP53 were associated with resistance to cisplatin and nutlin-3a, an inhibitor of MDM2 (Fig. 5a, Additional file 2: Table S1B). Only xenografted liver metastases were nonsynonymous mutant-type for TP53 by WES, which predicted a different response to cisplatin or nutlin-3a between donor and xenografted liver metastases. Similarly, donor and xenografted liver metastases demonstrated distinct sensitivity to the pan-ERBB inhibitor AZD8931 and the chemotherapeutic gemcitabine.
We also performed a validation screen in transcriptome analysis and compared expression of responses-associated genes between donor and xenografted specimens (Fig. 5b). A total of 30 drug response-associated target genes were summarized in Fig. 5b. Both donor and xenografted liver metastases demonstrated identical drug response to cisplatin, olaparib, PD-173074, cetuximab, trastuzumab, afatinib, AZD8931, lapatinib, trametinib, SCH772984, EMD-1214063, deltarasin, MK-2206, nutlin-3a,PD-0332991, LGK974, LY2109761, BIRB-1532 and dasatanib. In contrast, donor and xenografted liver metastases demonstrated different response to KU55933, axitinib, Sorafenib, CH5424802, taselisib, QSI-027, vorinostat, dabrafenib, BIRB 0796 and GSK126 (Fig. 5b).
PDOX-LM models mimicking clinical drug response of chemotherapy
In order to verify if our PDOX-LM models reliably recapitulate patients’ response to clinical drugs, we looked for established models in which the corresponding donor patients had been treated with mFOLFOX6 chemotherapy (Oxaliplatin and 5-FU were both key component). And this PDOX-LM models can predict the drug response of clinical chemotherapy treatment (Fig. 6a). Patient #3, #4, #5 and #6 were used in drug response testing, showing primary resistance to chemotherapy (Patient #3, #4, #6: partial response; Patient #5: progressive disease; according to RECIST 1.1 criteria; Fig. 6b). PDOX-LM models were divided into the control group (n = 16), PR group (n = 12) and PD group (n = 4). The overall survival (OS) of PR group were significantly prolonged (Fig. 6c), which was closely recapitulated the clinical drug response of the donor patients.
In the present study, we established 6 well-characterized PDOX-LM models. We combined organoid culture and in vivo splenic injection technology to obtain a rapid formation of metastatic xenografts in murine liver. As previous reports, in vivo models mimicking metastases have proven more informative than subcutaneous grown models in CRC [5, 37, 38], with improved microenviroment-metastases interaction [12,13,14], and evaluation of drug efficacy [39, 40]. And this study firstly revealed a rapid way to establish in vivo metastatic models based on primary tumor derived organoids with a high success rate. This PDOX-LM models could faithfully recapitulate molecular profiles and metastatic behavior of the metastatic tumors in mCRC patients.
We used organoid culture to closely recapitulate the genetic and morphological heterogeneous composition of the original tumor, offering great promise as a pre-clinical cancer model. According to the reports by Van de Wetering et al. , PDOs gave rise to large numbers of different organoids (10–1000), suggesting that the heterogeneous composition of the original tumor was largely conserved. In line with the previous results, we demonstrated that the organoids derived from mCRC presented a more vigorous reproductive capacity than those from non-metastatic CRC. Specifically, our results indicated a high density of organoids (100,000 cells/ml) may be an important contributing factor for the establishment of a PDOX-LM model. This may hasten more cancer stem cells and account for a successful establishment of in vivo metastatic model based on primary tumor-derived organoids.
The PDX approach has shown its high degree of translatability to the derived patient and offer in vivo models to study drug resistance and design novel treatment schedules in pre-clinical researches previously. However, it is reported that CRC PDX model require 6–8 months of propagation to be useful for treatment, making PDX model unsuitable for timely prediction for sensitive drugs . And no more than 5% of the patient-derived liver metastases were successfully formed in murine model [5, 41], because most mice died of giant xenografts in situ before the metastases occurred. And the successful xenografted liver metastases took quite a longer period (about 6 months since transplantation). Even in excellent work by Arianna Fumagalli , their murine models also have some limitations. The liver metastases were high frequently accompanied by lung metastases in murine model. All PDOs need further editing in cancer driver gene to enhance invasiveness before transplantation, which caused high cost and are not easy to master the technology. Our PDOX-LM models have overcome those limitations, because it possessed the characteristics such as liver specific metastases, without gene editing, and simple operation for establishment. The establishment of PDOX-LM models took a mean period of 7 weeks, equivalent to the length of time to establishing CRC cell lines-based live metastases model. In addition, unlike subcutaneous xenografts, an increase of CEA levels was detectable in murine serum during the development of xenografted liver metastases . Therefore, our model can mimic the biological behavior of human CRC in murine blood, facilitating the dynamical monitor of metastatic xenografts.
The present study further confirmed that the PDOX-LM could recapitulate molecular and phenotypic heterogeneity of donor liver metastases. It has been well established that Wnt and MAPK pathways account for CRC progression and metastases [13, 25, 27], and the presence of APC, KRAS, BRAF, PIK3CA, SMAD4 and TP53 mutations closely relates to treatment decision and prognosis. The sequencing data verified that a similar genomic alteration in KRAS, BRAF, PIK3CA and SMAD4 in both xenografted and donor liver metastases. It suggested that the common maintenance of subclone comprising major driver gene mutations during tumor progression in donor tumor or mice model . Interestingly, xenografted and donor liver metastases have discordant driving alterations in APC and TP53. It is very likely that the newly established model exhibit different degrees of fitness and pathway activity because of the aggressive organoids subpopulations had greater survival rates in vivo.
Furthermore, the transcriptome analysis verified that the PDOX-LM model could replicate expression of driver genes of the donor liver metastases. High similarity RNA profiles were observed between human donor and xenografted liver metastases. Especially for driver genes, such as APC, KRAS, BRAF, PIK3CA, SMAD4 and TP53, donor and xenografted liver metastases demonstrated highly consistent in both cases. According to Fang et al. , metastatic lesions are enriched in mutations of genes affecting PI3K-Akt signaling, cell adhesion and stellate-cell activation in the liver, the predominant metastasis site in patients. In this study, donor and xenografted liver metastases were highly consistent in expression of these genes, which means metastatic xenografts may narrow the difference in some cancer driver mutations in donor liver metastases.
In view of clinical application, we hypothesized that the PDOX-LM model can help screen the sensitive drugs. Given the complex molecular makeup of the tumor, it is important to identify the sensitive regiments to achieve a better treatment outcome.
Thus, using the WES and transcriptome analysis, we explored the power of PDOX-LM model to predict the drug response. According to ESMO guideline, cetuximab is used as the first-line target agent for mCRC, and our models could mimic the drug response of donor liver metastases by examining KRAS mutation status and expression level. Furthermore, it could exhibit similar response in multiple drugs, such as trastuzumab, afatinib, AZD8931, lapatinib, trametinib and SCH772984, which was due to the similar mutation status and expression in these drug-associated gene. However, our models failed to mimic the response of axitinib, sorafenib, taselisib and dabrafenib. One reason for this failure is that donor and xenografted liver metastases may partly root from different subclones of the identical primary tumor . And another reason is that only tumor epithelial, but not mesenchymal niche, were maintained during the culture of organoids.
To verify whether our model can mimic the drug response in clinical tumor, we also performed the drug response testing. It showed that PDOX-LM models indeed recapitulated the clinical drug response to mFOLFOX6 in the donor patients. Notably, according to Ooft SN et al. , the PDO culture unable to predict the response of the combination oxaliplatin-5FU. And the absence of stromal cell in the PDO culture may partially account for this result [6, 9]. The existence of stromal cells compromise the anti-cancer effect of the cytotoxic chemotherapy [6, 45,46,47,48,49]. In our PDOX-LM models, the PDOX involved stromal cells in tumor xenografts development in vivo, so our PDOX-LM models could recapitulate the response of oxaliplatin and 5-Fu in donor tumor. Our drug testing results were similar as the published PDX model based on patient-derived colon cancer cells by Isabel Puig et al. . Moreover, the rapid establishment of PDOX-LM models could increase the possibility of translational clinical application.
Indeed, this study has some limitations. First, more tumor samples should be included to offer further insight into the homogeneity and heterogeneity between donor and xenografted liver metastases. Second, targeted agents need to be explored by this model. Moreover, more biomarker, such as circulating tumor cells or DNA, could be verified in further by testing the murine blood.
In conclusion, the PDOX-LM models, reflecting the pheno- and genotype of the original patient samples, facilitate the rational selection of appropriate metastatic models for future drug experiments. This experimental approach could reduce the high failure rate observed in early preclinical trials in mCRC, which have no regard for the biological or molecular grounds when selecting experimental drugs or the patient cohort. And it is possible that advanced genetic analysis will discover novel genetic alterations accounting for the malignant progression and drug response.
Availability of data and materials
The authors declare that all the data supporting the findings of this study are available within the article.
Metastatic colorectal cancer
Colorectal liver metastases
- PDOX-LM model:
Patient-derived organoids-based xenografted liver metastases model
- H&E staining:
Hematoxylin and eosin staining
- IF staining:
- IHC staining:
Cancer Cell Fractions
Response Evaluation Criteria in Solid Tumors
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69:7–34.
Zhang L, Zhang L, Wang H, Chen L, Sui G. Diagnostic performance of contrast-enhanced ultrasound and magnetic resonance imaging for detecting colorectal liver metastases: a systematic review and meta-analysis. Dig Liver Dis. 2019;51:1241–8.
Troussard X, Cornet E. Outline for writing an article for current treatment options in oncology: splenic lymphoma with villous lymphocytes. Curr Treat Options Oncol. 2007;8:97–108.
Brown KM, Xue A, Julovi SM, Gill AJ, Pavlakis N, Samra JS, et al. Using patient-derived xenograft models of colorectal liver metastases to predict chemosensitivity. J Surg Res. 2018;227:158–67.
Hoffman RM. Patient-derived orthotopic xenografts: better mimic of metastasis than subcutaneous xenografts. Nat Rev Cancer. 2015;15:451–2.
Romano G, Chagani S, Kwong LN. The path to metastatic mouse models of colorectal cancer. Oncogene. 2018;37:2481–9.
Martin ES, Belmont PJ, Sinnamon MJ, Richard LG, Yuan J, Coffee EM, et al. Development of a colon cancer GEMM-derived orthotopic transplant model for drug discovery and validation. Clin Cancer Res. 2013;19:2929–40.
Inoue A, Deem AK, Kopetz S, Heffernan TP, Draetta GF, Carugo A. Current and future horizons of patient-derived xenograft models in colorectal cancer translational research. Cancers (Basel). 2019;11:1321.
Ooft SN, Weeber F, Dijkstra KK, McLean CM, Kaing S, van Werkhoven E, et al. Patient-derived organoids can predict response to chemotherapy in metastatic colorectal cancer patients. Sci Transl Med. 2019;11:eaay2574.
Sato T, Stange DE, Ferrante M, Vries RG, Van Es JH, Van den Brink S, et al. Long-term expansion of epithelial organoids from human colon, adenoma, adenocarcinoma, and Barrett’s epithelium. Gastroenterology. 2011;141:1762–72.
Liu X, Cheng Y, Abraham JM, Wang Z, Wang Z, Ke X, et al. Modeling Wnt signaling by CRISPR-Cas9 genome editing recapitulates neoplasia in human Barrett epithelial organoids. Cancer Lett. 2018;436:109–18.
Jackstadt R, van Hooff SR, Leach JD, Cortes-Lavaud X, Lohuis JO, Ridgway RA, et al. Epithelial NOTCH signaling rewires the tumor microenvironment of colorectal cancer to drive poor-prognosis subtypes and metastasis. Cancer Cell. 2019;36:319–36.
Calon A, Espinet E, Palomo-Ponce S, Tauriello DV, Iglesias M, Cespedes MV, et al. Dependency of colorectal cancer on a TGF-beta-driven program in stromal cells for metastasis initiation. Cancer Cell. 2012;22:571–84.
Zomer A, Maynard C, Verweij FJ, Kamermans A, Schafer R, Beerling E, et al. In Vivo imaging reveals extracellular vesicle-mediated phenocopying of metastatic behavior. Cell. 2015;161:1046–57.
Schwitalla S, Fingerle AA, Cammareri P, Nebelsiek T, Goktuna SI, Ziegler PK, et al. Intestinal tumorigenesis initiated by dedifferentiation and acquisition of stem-cell-like properties. Cell. 2013;152:25–38.
van de Wetering M, Francies HE, Francis JM, Bounova G, Iorio F, Pronk A, et al. Prospective derivation of a living organoid biobank of colorectal cancer patients. Cell. 2015;161:933–45.
Fisher S, Barry A, Abreu J, Minie B, Nolan J, Delorey TM, et al. A scalable, fully automated process for construction of sequence-ready human exome targeted capture libraries. Genome Biol. 2011;12:R1.
Imielinski M, Berger AH, Hammerman PS, Hernandez B, Pugh TJ, Hodis E, et al. Mapping the hallmarks of lung adenocarcinoma with massively parallel sequencing. Cell. 2012;150:1107–20.
Ojesina AI, Lichtenstein L, Freeman SS, Pedamallu CS, Imaz-Rosshandler I, Pugh TJ, et al. Landscape of genomic alterations in cervical carcinomas. Nature. 2014;506:371–5.
Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013;31:213–9.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43:491–8.
Chiang DY, Getz G, Jaffe DB, O’Kelly MJ, Zhao X, Carter SL, et al. High-resolution mapping of copy-number alterations with massively parallel sequencing. Nat Methods. 2009;6:99–103.
Olshen AB, Venkatraman ES, Lucito R, Wigler M. Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics. 2004;5:557–72.
Mermel CH, Schumacher SE, Hill B, Meyerson ML, Beroukhim R, Getz G. GISTIC20 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol. 2011;12:R41.
Cancer Genome Atlas Network. Comprehensive molecular characterization of human colon and rectal cancer. Nature. 2012;487:330–7.
Matano M, Date S, Shimokawa M, Takano A, Fujii M, Ohta Y, et al. Modeling colorectal cancer using CRISPR-Cas9-mediated engineering of human intestinal organoids. Nat Med. 2015;21:256–62.
Vitiello PP, De Falco V, Giunta EF, Ciardiello D, Cardone C, Vitale P, et al. Clinical practice use of liquid biopsy to identify RAS/BRAF mutations in patients with metastatic colorectal cancer (mCRC): A single institution experience. Cancers (Basel). 2019;11:1504.
Li G, Pavlick D, Chung JH, Bauer T, Tan BA, Peguero J, et al. Genomic profiling of cell-free circulating tumor DNA in patients with colorectal cancer and its fidelity to the genomics of the tumor biopsy. J Gastrointest Oncol. 2019;10:831–40.
Han SW, Kim HP, Shin JY, Jeong EG, Lee WC, Lee KH, et al. Targeted sequencing of cancer-related genes in colorectal cancer using next-generation sequencing. PLoS ONE. 2013;8:e64271.
Petaccia DMM, Melo FM, Ribeiro H, Marques MC, Kagohara LT, Begnami MD, et al. KRAS mutation status is highly homogeneous between areas of the primary tumor and the corresponding metastasis of colorectal adenocarcinomas: one less problem in patient care. Am J Cancer Res. 2017;7:1978–89.
Lawrence MS, Stojanov P, Mermel CH, Robinson JT, Garraway LA, Golub TR, et al. Discovery and saturation analysis of cancer genes across 21 tumour types. Nature. 2014;505:495–501.
Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, et al. A census of human cancer genes. Nat Rev Cancer. 2004;4(3):177–83.
Bourdais R, Rousseau B, Pujals A, Boussion H, Joly C, Guillemin A, et al. Polymerase proofreading domain mutations: new opportunities for immunotherapy in hypermutated colorectal cancer beyond MMR deficiency. Crit Rev Oncol Hematol. 2017;113:242–8.
Bao Y, Guo Y, Yang Y, Wei X, Zhang S, Zhang Y, et al. PRSS8 suppresses colorectal carcinogenesis and metastasis. Oncogene. 2019;38(4):497–517.
Chadwick RB, Jiang GL, Bennington GA, Yuan B, Johnson CK, Stevens MW, et al. Candidate tumor suppressor RIZ is frequently involved in colorectal carcinogenesis. Proc Natl Acad Sci USA. 2000;97:2662–7.
Iglesias D, Fernandez-Peralta AM, Nejda N, Daimiel L, Azcoita MM, Oliart S, et al. RIS1, a gene with trinucleotide repeats, is a target in the mutator pathway of colorectal carcinogenesis. Cancer Genet Cytogenet. 2006;167:138–44.
Bettenworth D, Mucke MM, Schwegmann K, Faust A, Poremba C, Schafers M, et al. Endoscopy-guided orthotopic implantation of colorectal cancer cells results in metastatic colorectal cancer in mice. Clin Exp Metastasis. 2016;33:551–62.
Prieto V, Ludwig JM, Farris AB, Nagaraju GP, Lawal TO, El-Rayes B, et al. Establishment of human metastatic colorectal cancer model in rabbit liver: a pilot study. PLoS ONE. 2017;12:e177212.
Bhagavatula SK, Upadhyaya K, Miller BJ, Bursch P, Lammers A, Cima MJ, et al. An interventional image-guided microdevice implantation and retrieval method for in vivo drug response assessment. Med Phys. 2019;46:5134–43.
Thiel C, Cordes H, Conde I, Castell JV, Blank LM, Kuepfer L. Model-based contextualization of in vitro toxicity data quantitatively predicts in vivo drug response in patients. Arch Toxicol. 2017;91:865–83.
Fu XY, Besterman JM, Monosov A, Hoffman RM. Models of human metastatic colon cancer in nude mice orthotopically constructed by using histologically intact patient specimens. Proc Natl Acad Sci USA. 1991;88:9345–9.
Fumagalli A, Suijkerbuijk S, Begthel H, Beerling E, Oost KC, Snippert HJ, et al. A surgical orthotopic organoid transplantation approach in mice to visualize and study colorectal cancer progression. Nat Protoc. 2018;13:235–47.
Metildi CA, Kaushal S, Luiken GA, Talamini MA, Hoffman RM, Bouvet M. Fluorescently labeled chimeric anti-CEA antibody improves detection and resection of human colon cancer in a patient-derived orthotopic xenograft (PDOX) nude mouse model. J Surg Oncol. 2014;109:451–8.
Fang W, Radovich M, Zheng Y, Fu CY, Zhao P, Mao C, et al. ‘Druggable’ alterations detected by Ion Torrent in metastatic colorectal cancer patients. Oncol Lett. 2014;7(6):1761–6.
Mao Y, Feng Q, Zheng P, Yang L, Liu T, Xu Y, et al. Low tumor purity is associated with poor prognosis, heavy mutation burden, and intense immune phenotype in colon cancer. Cancer Manage Res. 2018;10:3569–77.
Adamska A, Elaskalani O, Emmanouilidi A, Kim M, Abdol RN, Metharom P, et al. Molecular and cellular mechanisms of chemoresistance in pancreatic cancer. Adv Biol Regul. 2018;68:77–87.
Garg AD, More S, Rufo N, Mece O, Sassano ML, Agostinis P, et al. Trial watch: immunogenic cell death induction by anticancer chemotherapeutics. Oncoimmunology. 2017;6(12):e1386829.
Bezu L, Gomes-de-Silva LC, Dewitte H, Breckpot K, Fucikova J, Spisek R, et al. Combinatorial strategies for the induction of immunogenic cell death. Front Immunol. 2015;6:187.
Zitvogel L, Apetoh L, Ghiringhelli F, Kroemer G. Immunological aspects of cancer chemotherapy. Nat Rev Immunol. 2008;8:59–73.
Puig I, Chicote I, Tenbaum SP, Arques O, Herance JR, Gispert JD, et al. A personalized preclinical model to evaluate the metastatic potential of patient-derived colon cancer initiating cells. Clin Cancer Res. 2013;19:6787–801.
This study was supported by Grants from the National Natural Science Foundation of China (81602035, 81272390), the Shanghai Engineering Research Center of Colorectal Cancer Minimally Invasive (17DZ2250600) and the Shanghai Municipal Health Commission: Shanghai Outstanding Youth Specialist Training Program (Q2017-059), Clinical Science and Technology Innovation Project of Shanghai (SHDC12016104), and Shanghai Science and Technology Committee Project (17411951300, S2019-069).
Ethical approval and consent to participate
Ethics approval and consent to participate all procedures performed in this study involving human participants and animals were approved by the Ethics Committee of Zhongshan Hospital. Written informed consent was signed by each patient.
Consent for publication
The authors report no conflicts of interest in this work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 2: Table S1. A. 356 exonic SNV matched with TCGA-COAD on dbSNP Ids in PDOX-LM models. B. 36 Recurrent mutations in PDOX-LM models. C Somatic mutation rates in PDOX-LM models. D Cancer Cell Fractions rates; E. Discordant mutation between xenografted and donor liver metastases. A–E related to Figs. 3 and 5.
About this article
Cite this article
Jian, M., Ren, L., He, G. et al. A novel patient-derived organoids-based xenografts model for preclinical drug response testing in patients with colorectal liver metastases. J Transl Med 18, 234 (2020). https://doi.org/10.1186/s12967-020-02407-8
- Colorectal cancer
- Liver metastases
- Animal model