- Original article
- Open Access
Personalized treatment based on mini patient-derived xenografts and WES/RNA sequencing in a patient with metastatic duodenal adenocarcinoma
Cancer Communications volume 38, Article number: 54 (2018)
Treatment guidelines for a variety of cancers have been increasingly used in clinical practice, and have resulted in major improvement in patient outcomes. However, recommended regimens (even first-line treatments) are clearly not ideal for every patients. In the present study, we used mini patient-derived xenograft (mini-PDX) and next-generation sequencing to develop personalized treatment in a patient with metastatic duodenal adenocarcinoma.
Resected metachronous metastatic tumor tissues were implanted into SCID mice to determine the sensitivity to a variety of drug regimens. Mutation profiles were assessed using both DNA whole-exome sequencing (DNA–WES) and RNA sequencing. The results of the analyses were used to select optimal treatment for the patient with metastatic duodenal adenocarcinoma.
Assessment with mini-PDX models took only 7 days. The results showed high sensitivity to S-1 plus cisplatin, gemcitabine plus cisplatin and everolimus alone. The patient received gemcitabine plus cisplatin initially, but the treatment was terminated due to toxicity. The patient was then switched to treatment with S-1 alone. The overall disease-free survival was 34 months. DNA–WES and RNA sequencing identified KRAS mutation (A146T), TP53 (C229Yfs*10) and RICTOR amplification in the metastatic duodenal adenocarcinoma. These findings provided further support to the results of the mini-PDX, and suggest mTOR inhibitors should be used if and when relapse eventually occurs in this patient.
Mini-PDX model combined with WES/RNA sequencing can rapidly assess drug sensitivity in cancer patients and reveal key genetic alterations. Further research on this technology for personalized therapy in patients with refractory malignant tumors is warranted.
Duodenal adenocarcinoma is a rare tumor that accounts for less than 1% of all gastrointestinal cancers [1, 2]. Complete surgical resection is the optimal therapeutic modality for localized duodenal adenocarcinoma. Approximately 30–40% of the patients who have undergone curative resection of the primary duodenal adenocarcinoma eventually relapse . In lymph node-positive patients, adjuvant chemotherapy could improve patient outcomes . The addition of radiotherapy to adjuvant therapy does not confer a survival benefit in high-risk patients . Due to a lack of randomized studies, optimal adjuvant chemotherapy remains a matter of debate.
In recent years, clinical oncology has greatly benefited from rapid progress in molecular testing, particularly next-generation sequencing (NGS). NGS provides higher analytical sensitivity and could simultaneously analyze numerous target genes and pathways in cancer . Identification of driver mutations in tumors provides an opportunity to guide the selection of anticancer drugs and allows more precise targeted cancer therapy . RNA sequencing quantifies the level of gene expression at the transcriptional level, and could detect sequence variants , and has been used as a new tool to prospectively test somatic mutations and germline variants. Integrating information obtained with DNA whole exome-sequencing (DNA–WES) and RNA sequencing could improve detection of cancer driver mutations .
Next-generation sequencing in combination with patient-derived xenograft (PDX) models (developed by LIDE Biotech) has been used to develop personalized treatment for melanoma . PDXs allow preclinical evaluation of treatment effects , which, unlike cell line–derived tumor models, retain mutationally heterogeneous tumor cell populations similar to that in individual patients . PDX models may also effectively mimic the treatment response of the parental tumor and can provide clues for selecting therapeutic target and regimen [13, 14]. Major obstacles in using conventional PDXs to develop personalized treatment plan for individual patients include low tumor engraftment rate and long period of time. PDX modeling is also expensive, technically cumbersome, and requires a large amount of tumor tissue . Mini-PDX is a platform with reduced complexity and faster result turn-around. It is a promising tool where tumor cells from patients maintain tumorigenicity and are then inoculated into immunocompromised mice via a special capsule to establish tumor xenografts.
In the present study, we evaluated mini-PDX models combined with NGS (DNA–WES and RNA sequencing) analysis of metastatic tissue samples to guide personalized treatment in a patient with metastatic duodenal adenocarcinoma.
DNA was extracted from fresh metastatic tumor tissues using a QIAamp DNeasy blood and tissue kit (Qiagen, Valencia, CA, USA). Raw WES reads were aligned to the reference genome (hg19) using the BWA-MEM aligner and further re-mapped for error correction and re-calibration using the Assembly Based ReAligner (ABRA) . Polymerase chain reaction (PCR) duplications were removed by Picard software. The quality control metrics included coverage distribution, sequencing error and insert size estimation. Somatic mutations and genomic amplification/deletion were then detected from BAM files. Somatic single nucleotide variants (SNVs) were identified by MuTect (v1.17) with peripheral blood. Insertion-deletion polymorphisms (INDELs) were identified using PINDEL (V0.2.4). Additional filtering steps were used to reduce background noise. The functional impact of passed mutations was annotated by SnpEff3.0. The copy number variation (CNV) regions were identified by Control-FREEC (v9.4)  with the following parameters: window = 50,000 and step = 10,000. Gene fusions were detected through an in-house pipeline.
Total RNA was isolated using the RNeasy® Mini Kit (Qiagen), and the quality of RNA sequencing paired-end reads was assessed with FastQC. Reads were then mapped to the hg19 genome with the STAR RNA sequencing aligner . Gene fusions at the transcript level were detected by STAR-fusion. Gene expression level was estimated with RSEM.
Mini-PDX assay was carried out using the OncoVee® Mini PDX kit (LIDE Biotech Co., Ltd, Shanghai, China) (Fig. 1). Briefly, metastatic tumor samples were harvested, and then washed with Hanks’ balanced salt solution (HBSS) to remove non-tumor tissue and necrotic tumor tissue in a biosafety cabinet. The tumor samples were minced, followed by incubation with collagenase at 37 °C for 1–2 h. The cells were collected and blood cells and fibroblasts were then removed. The cell suspension was transferred to HBSS-washed capsules made of hollow fiber membrane with a pore size allowing passage of molecules less than 500 kDa. The fiber system delivered media to the cells in a manner similar to the delivery of blood through the capillary networks in vivo.
SCID mice (6–8 weeks of age) were housed and monitored in a SPF animal facility. For subcutaneous implantation, a small skin incision was made and the special capsule was inserted through the subcutaneous tissue. One day after tumor cell inoculation, tumor-bearing mice were randomized to receive vehicle, S-1 alone [10 mg/kg, orally administered (po), daily], oxaliplatin alone [5 mg/kg, intravenously (iv), every 4 days], irinotecan alone (50 mg/kg, iv, every 4 days), capecitabine alone (400 mg/kg, po, daily), nab-paclitaxel alone (20 mg/kg, iv, every 4 days), everolimus alone (10 mg/kg, po, daily), cisplatin (5 mg/kg, iv, every 4 days) in combination with S-1 (10 mg/kg, po, daily) or gemcitabine (60 mg/kg, iv, every 4 days). Seven days later, the capsules with tumor cells were removed. Cell viability was evaluated using a CTG assay.
Statistical analyses were performed using SPSS version 16.0 (SPSS Inc., Chicago, IL, USA). One-way ANOVA followed by posthoc analysis was used for comparisons. P < 0.05 was considered statistically significant.
A 57-year-old man underwent a pancreaticoduodenectomy for duodenal adenocarcinoma in February 2014. The tumor penetrated the visceral peritoneum and invaded the pancreas with nodal involvement (pT4N1M0, stage IIIA, Fig. 2a). Starting from 1 month after the surgery, he received empirical treatment with 6 cycles of gemcitabine plus S-1 (GS) (gemcitabine 1000 mg/m2 on days 1 and 8, and oral S-1 60 mg twice daily on days 1–14, repeated every 3 weeks). A repeat CT scan in December 2016 revealed a lesion 2.2 cm in diameter in the root of the superior mesenteric artery and below the transverse mesocolon (Fig. 3a). PET/CT scan showed enhanced FDG uptake (Fig. 3b). The metastatic lesion was surgically removed and pathologically confirmed as metastatic duodenal adenocarcinoma in December 2016 (Fig. 2b).
Based on the results of the mini-PDX assessment, the patient was given a chemotherapy regimen that consisted of intravenous gemcitabine (1000 mg/m2) and cisplatin (25 mg/m2) on day 1 and 8, given every 3 weeks. The treatment was terminated after 4 cycles due to drug-related toxicities, including grade 3 nausea and fatigue. In May 2017, the patient was treated with 3 cycles of S-1 treatment (60 mg, po, twice daily for 14 consecutive days every 3 weeks). After the second surgery, the patient had completed adjuvant chemotherapy in February 2017. The last follow-up (9 months later) showed did not reveal tumor burden. The treatments and responses are summarized in Fig. 4.
Mini-PDX assays showed that compared with the vehicle control on day 0, cell viability was decreased by 12.0% with S-1 plus cisplatin, and by 0.8% with gemcitabine plus cisplatin. The anti-proliferation rate with oxaliplatin, S-1, and everolimus was 67.6%, 80.3%, and 90.7%, respectively on day 7 (Fig. 5). Cell viability in the nab-paclitaxel, irinotecan and capecitabine group was not significantly different from the vehicle control (P > 0.05). No significant differences in body weight were observed among different treatment groups.
DNA–WES and RNA sequencing analysis
DNA–WES and RNA sequencing were applied to a panel of 450 relevant cancer genes in metastatic duodenal adenocarcinoma tissues. Targeted gene sequencing was achieved in both tumor and normal samples with a mean coverage of 1728X and 431X, respectively. More than 85% of the reads were mapped to regions that the panel targeted or flanked in both tumor and normal samples. We identified the following somatic genetic changes in the metastatic tumor: KRAS mutation (A146T), APC (N1455Kfs*2), TP53 (C229Yfs*10), CXCR4 (V160I), TERT amplification, MCL1 amplification and RICTOR amplification. Furthermore, we detected the following pharmacogenomic biomarkers: cisplatin (GSTP1 rs1695 A/A), fluorouracil (TPMT rs1142345 A/A), and irinotecan (UGT1A1 rs8175347 6/6; SLCO1B1 rs2306283 G/A).
To confirm these mutations, we performed paired-end RNA sequencing of the metastatic sample to a depth of 52 million reads. Ninety-five percent of the reads (49 M) were aligned to the human genome (Additional file 1: Tables S1, S2). Of the aligned reads, 22.82% (12 M) were protein coding (Additional file 1: Figures S1, S2). Consistent with the results of WES, we identified the following somatic mutations in the transcriptomes of the metastatic tumor tissue: KRAS mutation (A146T), TP53 (C229Yfs*10) and RICTOR amplification.
Currently, little data is available to guide rational choice of chemotherapeutic regimen for duodenal adenocarcinoma. The current practice at many centers is to treat duodenal adenocarcinoma patients with oxaliplatin-based chemotherapy similar to that for colorectal adenocarcinoma . Gemcitabine-based adjuvant chemotherapy was offered to eligible patients with peri-ampullary duodenal adenocarcinoma in the ESPAC-3 trial . In vivo models closely mimicking the biology of duodenal adenocarcinoma in patients are urgently needed to reliably determine optimal drug sensitivities of individual cancers for personalized therapy. PDX models of certain cancers closely mirror the drug response in patients compared with other models . However, PDX models have several drawbacks that limit their clinical application. Generation of PDXs is only modestly successful and requires a large amount of tumor tissue which can usually only be obtained from surgically resected tumors. Furthermore, PDX models require 4–8 months to identify active therapies for a particular cancer patient. Elapse of time between tumor engraftment in mice and initiating treatment for patients is a limiting factor for real-time personalized medicine applications [12, 21]. To overcome these limitations, a new technology known as the mini-PDX model has emerged to guide selection of chemotherapeutic drugs. Mini-PDX models only require a small number of tumor cells and allow rapid analysis of drug sensitivity with a median latency of 7 days; thus, duodenal adenocarcinoma patients can receive personalized chemotherapy in a clinically relevant time frame. Previous studies have demonstrated a potentially strong agreement between mini-PDX- and PDXs-based prediction of drug sensitivity (overall response agreement: 89%) in a variety of solid tumors including lung cancer, pancreatic cancer, and gastric cancer. Importantly, the drug sensitivity pattern of mini-PDX can recapitulate responses in patients from which they were derived.
The median disease-free survival of duodenal adenocarcinoma patients with lymph node metastasis who underwent radical resection is 18–25 months [22, 23]. If recurrence develops, the survival time is only 5.4 months . However, the disease-free survival of our patient who received GS as an adjuvant chemotherapeutic regimen was 34 months. Sensitivity of the cancer tissues to gemcitabine and S-1 in the mini-PDX experiments seems to be consistent with clinical response to the GS regimen in this patient. Following a second resection, the patient received gemcitabine and S-1-based chemotherapy according to the mini-PDX results. Therefore, mini-PDX is suitable for selecting optimal regimens for personalized therapy clinically. However, further studies with a large patient cohort are required.
DNA–WES and RNA sequencing identified patients with KRAS mutation and RICTOR-amplified tumors. KRAS mutation correlates significantly with late stage and poor tumor differentiation in duodenal adenocarcinoma and KRAS-mutated persons are more likely to experience distant relapse . To date, there are no effective anti-Ras strategies. Activated KRAS induces a complex signaling network, including the PI3K/AKT/mTOR pathway. Due to difficulty in targeting KRAS directly, targeting downstream RAS effector signaling pathways are currently under clinical evaluation [25, 26]. RICTOR plays an important role in the PI3K/AKT/mTOR signaling pathway, which is one of the most commonly activated pathways in human cancer . RICTOR is an essential component of the mTORC2 complex and upregulation of mTORC2 activity in turn regulates cell growth, metabolism, and survival . It is also an important component of feedback loops and cross-talk in PI3K–AKT–mTOR signaling. RICTOR amplification is a rare, but recurrent somatic alteration in solid tumors . Targeting mTOR to inhibit RICTOR may be a potential therapeutic strategy . In the mini-PDX model, RICTOR amplified duodenal adenocarcinoma demonstrated sensitivity to mTORC1 inhibitor everolimus. However, inhibition of mTORC1 leads to an increase in AKT phosphorylation due to feedback loops that allow continued activity of mTORC2 . The dual mTORC1/2 inhibitor, AZD2014, showed a more potent inhibitory effect on RICTOR-amplified tumors than everolimus .
Pro-survival protein myeloid cell leukemia 1 (MCL1) amplification was identified by DNA sequencing in this patient. MCL1, one of the anti-apoptotic BCL-2 family members, is overexpressed in many human cancers , and in several key oncogenic pathways to sustain cancer cell survival and resist cancer treatment. Repression of MCL1 exerts cytotoxic effects in tumor cells , indicating that MCL1 is a promising target for the treatment of a wide range of tumors. S63845, a small molecule that specifically binds with high affinity to the BH3-binding groove of MCL1, exerts anti-tumor activity in MCL1-dependent cancer cells in vitro and in vivo, while sparing normal tissues at efficacious doses, which represents a potential breakthrough in cancer therapy .
Mini-PDX models in combination with DNA–WES/RNA analysis as an effective platform could be used to optimize the clinical management of metastatic duodenal adenocarcinoma.
Overman MJ, Hu CY, Kopetz S, Abbruzzese JL, Wolff RA, Chang GJ. A population-based comparison of adenocarcinoma of the large and small intestine: insights into a rare disease. Ann Surg Oncol. 2012;19(5):1439–45.
Cloyd JM, George E, Visser BC. Duodenal adenocarcinoma: advances in diagnosis and surgical management. World J Gastrointest Surg. 2016;8(3):212–21.
Lee HG, You DD, Paik KY, Heo JS, Choi SH, Choi DW. Prognostic factors for primary duodenal adenocarcinoma. World J Surg. 2008;32(10):2246–52.
Ecker BL, McMillan MT, Datta J, Mamtani R, Giantonio BJ, Dempsey DT, et al. Efficacy of adjuvant chemotherapy for small bowel adenocarcinoma: a propensity score-matched analysis. Cancer. 2016;122(5):693–701.
Ecker BL, McMillan MT, Datta J, Lee MK, Karakousis GC, Vollmer CM Jr, et al. Adjuvant chemotherapy versus chemoradiotherapy in the management of patients with surgically resected duodenal adenocarcinoma: a propensity score-matched analysis of a nationwide clinical oncology database. Cancer. 2017;123(6):967–76.
Sheikine Y, Kuo FC, Lindeman NI. Clinical and technical aspects of genomic diagnostics for precision oncology. J Clin Oncol. 2017;35(9):929–33.
Garnett MJ, Edelman EJ, Heidorn SJ, Greenman CD, Dastur A, Lau KW, et al. Systematic identification of genomic markers of drug sensitivity in cancer cells. Nature. 2012;483(7391):570–5.
Cirulli ET, Singh A, Shianna KV, Ge D, Smith JP, Maia JM, et al. Center for HIVAVI. Screening the human exome: a comparison of whole genome and whole transcriptome sequencing. Genome Biol. 2010;11(5):R57.
Wilkerson MD, Cabanski CR, Sun W, Hoadley KA, Walter V, Mose LE, et al. Integrated rna and DNA sequencing improves mutation detection in low purity tumors. Nucleic Acids Res. 2014;42(13):e107.
Girotti MR, Gremel G, Lee R, Galvani E, Rothwell D, Viros A, et al. Application of sequencing, liquid biopsies, and patient-derived xenografts for personalized medicine in melanoma. Cancer Discov. 2016;6(3):286–99.
Aparicio S, Hidalgo M, Kung AL. Examining the utility of patient-derived xenograft mouse models. Nat Rev Cancer. 2015;15(5):311–6.
Hidalgo M, Amant F, Biankin AV, Budinska E, Byrne AT, Caldas C, et al. Patient-derived xenograft models: an emerging platform for translational cancer research. Cancer Discov. 2014;4(9):998–1013.
Izumchenko E, Meir J, Bedi A, Wysocki PT, Hoque MO, Sidransky D. Patient-derived xenografts as tools in pharmaceutical development. Clin Pharmacol Ther. 2016;99(6):612–21.
Zhang L, Nomie K, Zhang H, Bell T, Pham L, Kadri S, et al. B-cell lymphoma patient-derived xenograft models enable drug discovery and are a platform for personalized therapy. Clin Cancer Res. 2017;23(15):4212–23.
Byrne AT, Alferez DG, Amant F, Annibali D, Arribas J, Biankin AV, Trusolino L, et al. Interrogating open issues in cancer precision medicine with patient-derived xenografts. Nat Rev Cancer. 2017;17(4):254–68.
Mose LE, Wilkerson MD, Hayes DN, Perou CM, Parker JS. Abra: improved coding indel detection via assembly-based realignment. Bioinformatics. 2014;30(19):2813–5.
Boeva V, Popova T, Bleakley K, Chiche P, Cappo J, Schleiermacher G, et al. Control-freec: a tool for assessing copy number and allelic content using next-generation sequencing data. Bioinformatics. 2012;28(3):423–5.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Gingeras TR, et al. Star: ultrafast universal rna-seq aligner. Bioinformatics. 2013;29(1):15–21.
Neoptolemos JP, Moore MJ, Cox TF, Valle JW, Palmer DH, McDonald AC, European Study Group for Pancreatic C, et al. Effect of adjuvant chemotherapy with fluorouracil plus folinic acid or gemcitabine vs observation on survival in patients with resected periampullary adenocarcinoma: the espac-3 periampullary cancer randomized trial. JAMA. 2012;308(2):147–56.
Garber K. From human to mouse and back: ‘Tumorgraft’ models surge in popularity. J Natl Cancer Inst. 2009;101(1):6–8.
Hwang CI, Boj SF, Clevers H, Tuveson DA. Preclinical models of pancreatic ductal adenocarcinoma. J Pathol. 2016;238(2):197–204.
Solaini L, Jamieson NB, Metcalfe M, Abu Hilal M, Soonawalla Z, Davidson BR, Group UKDCS., et al. Outcome after surgical resection for duodenal adenocarcinoma in the uk. Br J Surg. 2015;102(6):676–81.
Jiang QL, Huang XH, Zhang JW, Zhang SS, Wang CF. Evaluation of disease free survival after radical resection for primary duodenal adenocarcinoma. Zhonghua Yi Xue Za Zhi. 2016;96(42):3393–7.
Fu T, Guzzetta AA, Jeschke J, Vatapalli R, Dave P, Hooker CM, et al. Kras g>a mutation favors poor tumor differentiation but may not be associated with prognosis in patients with curatively resected duodenal adenocarcinoma. Int J Cancer. 2013;132(11):2502–9.
Baines AT, Xu D, Der CJ. Inhibition of ras for cancer treatment: the search continues. Future Med Chem. 2011;3(14):1787–808.
Wang Y, Kaiser CE, Frett B, Li HY. Targeting mutant kras for anticancer therapeutics: a review of novel small molecule modulators. J Med Chem. 2013;56(13):5219–30.
Rodrik-Outmezguine VS, Okaniwa M, Yao Z, Novotny CJ, McWhirter C, Banaji A, et al. Overcoming mtor resistance mutations with a new-generation mtor inhibitor. Nature. 2016;534(7606):272–6.
Guertin DA, Stevens DM, Thoreen CC, Burds AA, Kalaany NY, Moffat J, et al. Ablation in mice of the mtorc components raptor, rictor, or mlst8 reveals that mtorc2 is required for signaling to akt-foxo and pkcalpha, but not s6k1. Dev Cell. 2006;11(6):859–71.
Kim ST, Kim SY, Klempner SJ, Yoon J, Kim N, Ahn S, et al. Rapamycin-insensitive companion of mtor (rictor) amplification defines a subset of advanced gastric cancer and is sensitive to azd2014-mediated mtorc1/2 inhibition. Ann Oncol. 2017;28(3):547–54.
Driscoll DR, Karim SA, Sano M, Gay DM, Jacob W, Yu J, et al. Mtorc2 signaling drives the development and progression of pancreatic cancer. Cancer Res. 2016;76(23):6911–23.
Liu P, Gan W, Chin YR, Ogura K, Guo J, Zhang J, et al. Ptdins(3,4,5)p3-dependent activation of the mtorc2 kinase complex. Cancer Discov. 2015;5(11):1194–209.
Beroukhim R, Mermel CH, Porter D, Wei G, Raychaudhuri S, Donovan J, Meyerson M, et al. The landscape of somatic copy-number alteration across human cancers. Nature. 2010;463(7283):899–905.
Wei G, Margolin AA, Haery L, Brown E, Cucolo L, Julian B, Golub TR, et al. Chemical genomics identifies small-molecule mcl1 repressors and bcl-xl as a predictor of mcl1 dependency. Cancer Cell. 2012;21(4):547–62.
Kotschy A, Szlavik Z, Murray J, Davidson J, Maragno AL, Le Toumelin-Braizat G, et al. The mcl1 inhibitor s63845 is tolerable and effective in diverse cancer models. Nature. 2016;538(7626):477–82.
ZS and PZ designed the research. PZ and CY drafted the manuscript. YW and FZ performed the mini-PDX study. SM performed NGS. HL analyzed the data. All authors read and approved the final manuscript.
We gratefully acknowledged the help provided by staff at the Department of General Surgery, Zhejiang Hospital for assistance with the data.
The authors declare that they have no competing interests.
Availability of data and materials
Consent for publication
The patient has given his consent for this case report to be published.
Ethics approval and consent to participate
Publication of this manuscript was approved by the Ethics Committee on Scientific Research of the First Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, China; informed consent was obtained from the patient.
This study was funded by the National Natural Science Foundation of China (No. 81472346) and the National Key Research and Development Program of China (No. 2017ZX10203205).
Additional file 1: Table S1. Quality control metrics of targeted sequencing data. Table S2. Quality control metrics of RNA-seq data. Figure S1. Average coverage (Log10 transformed) of targeted genes with data deduplication. Coverage of mutant genes were shown in grey bar, error bars indicted standard deviation of sequencing depths for each gene. Coverage of all genes sequenced were shown in red line. Figure S2. Fraction of reads that mapped to coding, intergenic, ribosomal, intronic, UTR and unmapped reads for metastasis RNA-seq.
About this article
Cite this article
Zhao, P., Chen, H., Wen, D. et al. Personalized treatment based on mini patient-derived xenografts and WES/RNA sequencing in a patient with metastatic duodenal adenocarcinoma. Cancer Commun 38, 54 (2018). https://doi.org/10.1186/s40880-018-0323-y
- Duodenal adenocarcinoma
- Mini patient-derived xenograft
- Whole-exome sequencing
- RNA sequencing
- Somatic mutation
- Personalized therapy