Whole genomes define concordance of matched primary, xenograft, and organoid models of pancreas cancer

Pancreatic ductal adenocarcinoma (PDAC) has the worst prognosis among solid malignancies and improved therapeutic strategies are needed to improve outcomes. Patient-derived xenografts (PDX) and patient-derived organoids (PDO) serve as promising tools to identify new drugs with therapeutic potential in PDAC. For these preclinical disease models to be effective, they should both recapitulate the molecular heterogeneity of PDAC and validate patient-specific therapeutic sensitivities. To date however, deep characterization of the molecular heterogeneity of PDAC PDX and PDO models and comparison with matched human tumour remains largely unaddressed at the whole genome level. We conducted a comprehensive assessment of the genetic landscape of 16 whole-genome pairs of tumours and matched PDX, from primary PDAC and liver metastasis, including a unique cohort of 5 'trios' of matched primary tumour, PDX, and PDO. We developed a pipeline to score concordance between PDAC models and their paired human tumours for genomic events, including mutations, structural variations, and copy number variations. Tumour-model comparisons of mutations displayed single-gene concordance across major PDAC driver genes, but relatively poor agreement across the greater mutational load. Genome-wide and chromosome-centric analysis of structural variation (SV) events highlights previously unrecognized concordance across chromosomes that demonstrate clustered SV events. We found that polyploidy presented a major challenge when assessing copy number changes; however, ploidy-corrected copy number states suggest good agreement between donor-model pairs. Collectively, our investigations highlight that while PDXs and PDOs may serve as tractable and transplantable systems for probing the molecular properties of PDAC, these models may best serve selective analyses across different levels of genomic complexity.

[1]  O. Myklebost,et al.  Generation and Characterisation of Novel Pancreatic Adenocarcinoma Xenograft Models and Corresponding Primary Cell Lines , 2014, PloS one.

[2]  Andrew H. Beck,et al.  PharmacoGx: an R package for analysis of large pharmacogenomic datasets , 2015, Bioinform..

[3]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[4]  Michael C. Rusch,et al.  CREST maps somatic structural variation in cancer genomes with base-pair resolution , 2011, Nature Methods.

[5]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[6]  Lincoln D. Stein,et al.  Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes , 2012, Nature.

[7]  A. Bouchard-Côté,et al.  PyClone: statistical inference of clonal population structure in cancer , 2014, Nature Methods.

[8]  Niklaus J Grünwald,et al.  vcfr: a package to manipulate and visualize variant call format data in R , 2017, Molecular ecology resources.

[9]  Yong-shu He,et al.  [Structural variation in the human genome]. , 2009, Yi chuan = Hereditas.

[10]  Igor Jurisica,et al.  Hypoxia predicts aggressive growth and spontaneous metastasis formation from orthotopically grown primary xenografts of human pancreatic cancer. , 2011, Cancer research.

[11]  Kenneth H. Yu,et al.  Preclinical models of pancreatic ductal adenocarcinoma. , 2017, Chinese clinical oncology.

[12]  Paul Shannon,et al.  VariantAnnotation: a Bioconductor package for exploration and annotation of genetic variants , 2014, Bioinform..

[13]  Agnieszka K Witkiewicz,et al.  Pancreatic cancer cell lines as patient-derived avatars: genetic characterisation and functional utility , 2017, Gut.

[14]  Elizabeth M. Smigielski,et al.  dbSNP: the NCBI database of genetic variation , 2001, Nucleic Acids Res..

[15]  M. DePristo,et al.  A framework for variation discovery and genotyping using next-generation DNA sequencing data , 2011, Nature Genetics.

[16]  P. Spellman,et al.  Subtypes of Pancreatic Ductal Adenocarcinoma and Their Differing Responses to Therapy , 2011, Nature Medicine.

[17]  Gordon Keller,et al.  Ductal pancreatic cancer modeling and drug screening using human pluripotent stem cell– and patient-derived tumor organoids , 2015, Nature Medicine.

[18]  Robert Gentleman,et al.  Software for Computing and Annotating Genomic Ranges , 2013, PLoS Comput. Biol..

[19]  Robert Gentleman,et al.  rtracklayer: an R package for interfacing with genome browsers , 2009, Bioinform..

[20]  Lincoln D. Stein,et al.  Association of Distinct Mutational Signatures With Correlates of Increased Immune Activity in Pancreatic Ductal Adenocarcinoma , 2017, JAMA oncology.

[21]  Thomas Zichner,et al.  DELLY: structural variant discovery by integrated paired-end and split-read analysis , 2012, Bioinform..

[22]  A. Børresen-Dale,et al.  Copynumber: Efficient algorithms for single- and multi-track copy number segmentation , 2012, BMC Genomics.

[23]  Aaron R. Quinlan,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2022 .

[24]  M. Spector,et al.  Organoid Models of Human and Mouse Ductal Pancreatic Cancer , 2015, Cell.

[25]  Tao Xie,et al.  Whole Exome Sequencing of Rapid Autopsy Tumors and Xenograft Models Reveals Possible Driver Mutations Underlying Tumor Progression , 2015, PloS one.

[26]  Gun Ho Jang,et al.  A renewed model of pancreatic cancer evolution based on genomic rearrangement patterns , 2016, Nature.

[27]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[28]  R. Gibbs,et al.  Genomic analyses identify molecular subtypes of pancreatic cancer , 2016, Nature.

[29]  Cassandra Willyard Copy number variations' effect on drug response still overlooked , 2015, Nature Medicine.

[30]  Wendy S. W. Wong,et al.  Strelka: accurate somatic small-variant calling from sequenced tumor-normal sample pairs , 2012, Bioinform..

[31]  R. Hruban,et al.  Molecular Characteristics of Pancreatic Ductal Adenocarcinoma , 2011, Pathology research international.

[32]  Hongen Zhang,et al.  RCircos: an R package for Circos 2D track plots , 2013, BMC Bioinformatics.

[33]  Trevor J Pugh,et al.  Oncotator: Cancer Variant Annotation Tool , 2015, Human mutation.

[34]  K. Lillemoe,et al.  Tumor engraftment in patient-derived xenografts of pancreatic ductal adenocarcinoma is associated with adverse clinicopathological features and poor survival , 2017, PloS one.

[35]  A. Ranga,et al.  Artificial three-dimensional niches deconstruct pancreas development in vitro , 2013, Development.

[36]  Huamin Wang,et al.  Clinical and Genetic Implications of DNA Mismatch Repair Deficiency in Patients With Pancreatic Ductal Adenocarcinoma , 2017, JAMA surgery.

[37]  H T Lynch,et al.  Review of the Lynch syndrome: history, molecular genetics, screening, differential diagnosis, and medicolegal ramifications , 2009, Clinical genetics.

[38]  Gordon B Mills,et al.  Integrated Patient-Derived Models Delineate Individualized Therapeutic Vulnerabilities of Pancreatic Cancer. , 2016, Cell reports.

[39]  Davide Prandi,et al.  Personalized In Vitro and In Vivo Cancer Models to Guide Precision Medicine. , 2017, Cancer discovery.

[40]  A. Sivachenko,et al.  Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples , 2013, Nature Biotechnology.

[41]  M Eileen Dolan,et al.  Clinically relevant genetic variations in drug metabolizing enzymes. , 2011, Current drug metabolism.

[42]  Obi L. Griffith,et al.  GenVisR: Genomic Visualizations in R , 2016, bioRxiv.

[43]  Todd R. Golub,et al.  Abstract 1028: Patient-derived xenografts undergo mouse-specific tumor evolution , 2018, Tumor Biology.

[44]  Steven J. M. Jones,et al.  Integrated Genomic Characterization of Pancreatic Ductal Adenocarcinoma. , 2017, Cancer cell.

[45]  Manuel Hidalgo,et al.  Patient-derived xenograft models: an emerging platform for translational cancer research. , 2014, Cancer discovery.

[46]  Mingming Jia,et al.  COSMIC: exploring the world's knowledge of somatic mutations in human cancer , 2014, Nucleic Acids Res..

[47]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[48]  J. Kench,et al.  Whole genomes redefine the mutational landscape of pancreatic cancer , 2015, Nature.

[49]  Jeremy Wazny,et al.  Xenome—a tool for classifying reads from xenograft samples , 2012, Bioinform..