Tumor heterogeneity assessed by sequencing and fluorescence in situ hybridization (FISH) data

Computational reconstruction of clonal evolution in cancers has become a crucial tool for understanding how tumors initiate and progress and how this process varies across patients. The field still struggles, however, with special challenges of applying phylogenetic methods to cancers, such as the prevalence and importance of copy number alteration (CNA) and structural variation (SV) events in tumor evolution, which are difficult to profile accurately by prevailing sequencing methods in such a way that subsequent reconstruction by phylogenetic inference algorithms is accurate. In the present work, we develop computational methods to combine sequencing with multiplex interphase fluorescence in situ hybridization (miFISH) to exploit the complementary advantages of each technology in inferring accurate models of clonal CNA evolution accounting for both focal changes and aneuploidy at whole-genome scales. We demonstrate on simulated data that incorporation of FISH data substantially improves accurate inference of focal CNA and ploidy changes in clonal evolution from deconvolving bulk sequence data. Analysis of real glioblastoma data for which FISH, bulk sequence, and single cell sequence are all available confirms the power of FISH to enhance accurate reconstruction of clonal copy number evolution in conjunction with bulk and optionally single-cell sequence data. Availability github.com/CMUSchwartzLab/FISH_deconvolution Contact russells@andrew.cmu.edu

[1]  Jack Kuipers,et al.  Integrative inference of subclonal tumour evolution from single-cell and bulk sequencing data , 2019, Nature Communications.

[2]  A. Schäffer,et al.  Single-cell genetic analysis of ductal carcinoma in situ and invasive breast cancer reveals enormous tumor heterogeneity yet conserved genomic imbalances and gain of MYC during progression. , 2012, The American journal of pathology.

[3]  Russell Schwartz,et al.  Inferring models of multiscale copy number evolution for single-tumor phylogenetics , 2015, Bioinform..

[4]  Benjamin J. Raphael,et al.  Inferring the Mutational History of a Tumor Using Multi-state Perfect Phylogeny Mixtures. , 2016, Cell systems.

[5]  Ron Shamir,et al.  Copy-Number Evolution Problems: Complexity and Algorithms , 2016, WABI.

[6]  Icgc,et al.  Pan-cancer analysis of whole genomes , 2017, bioRxiv.

[7]  Jonathan W Heusel,et al.  Beyond sequence variation: assessment of copy number variation in adult glioblastoma through targeted tumor somatic profiling. , 2019, Human pathology.

[8]  A. Schäffer,et al.  Aneuploidy, TP53 mutation, and amplification of MYC correlate with increased intratumor heterogeneity and poor prognosis of breast cancer patients , 2018, Genes, chromosomes & cancer.

[9]  Benjamin J. Raphael,et al.  Accurate quantification of copy-number aberrations and whole-genome duplications in multi-sample tumor sequencing data , 2018, Nature Communications.

[10]  Russell Schwartz,et al.  Tumor Copy Number Deconvolution Integrating Bulk and Single-Cell Sequencing Data. , 2020, Journal of computational biology : a journal of computational molecular cell biology.

[11]  Russell Schwartz,et al.  Phylogenies Derived from Matched Transcriptome Reveal the Evolution of Cell Populations and Temporal Order of Perturbed Pathways in Breast Cancer Brain Metastases , 2019, ISMCO.

[12]  Thomas Lengauer,et al.  Mtreemix: a software package for learning and using mixture models of mutagenetic trees , 2005, Bioinform..

[13]  A. Schäffer,et al.  The evolution of tumour phylogenetics: principles and practice , 2017, Nature Reviews Genetics.

[14]  Joshua M. Korn,et al.  Comprehensive genomic characterization defines human glioblastoma genes and core pathways , 2008, Nature.

[15]  Russell Schwartz,et al.  Reconstructing Tumor phylogenies from Heterogeneous Single-Cell Data , 2007, J. Bioinform. Comput. Biol..

[16]  Russell Schwartz,et al.  Deconvolution and phylogeny inference of structural variations in tumor genomic samples , 2018, bioRxiv.

[17]  E. Schröck,et al.  Comprehensive molecular characterization of multifocal glioblastoma proves its monoclonal origin and reveals novel insights into clonal evolution and heterogeneity of glioblastomas , 2017, Neuro-oncology.

[18]  Zoltan Szallasi,et al.  Tolerance of whole-genome doubling propagates chromosomal instability and accelerates cancer genome evolution. , 2014, Cancer discovery.

[19]  Ron Shamir,et al.  Complexity and algorithms for copy-number evolution problems , 2017, Algorithms for Molecular Biology.

[20]  A. McKenna,et al.  Absolute quantification of somatic DNA alterations in human cancer , 2012, Nature Biotechnology.

[21]  P. Nowell The clonal evolution of tumor cell populations. , 1976, Science.

[22]  Niko Beerenwinkel,et al.  Computational Cancer Biology: An Evolutionary Perspective , 2016, PLoS Comput. Biol..

[23]  Russell Schwartz,et al.  Algorithms to Model Single Gene, Single Chromosome, and Whole Genome Copy Number Changes Jointly in Tumor Phylogenetics , 2014, PLoS Comput. Biol..

[24]  Jun Zhou,et al.  Analysis of gene copy number changes in tumor phylogenetics , 2016, Algorithms for Molecular Biology.

[25]  B. Taylor,et al.  Genome doubling shapes the evolution and prognosis of advanced cancers , 2018, Nature Genetics.

[26]  B. Boisselier,et al.  Whole genome duplication is an early event leading to aneuploidy in IDH-wild type glioblastoma , 2018, Oncotarget.

[27]  Russell Schwartz,et al.  Tumor Copy Number Deconvolution Integrating Bulk and Single-Cell Sequencing Data , 2018, bioRxiv.

[28]  Nancy R. Zhang,et al.  Assessing intratumor heterogeneity and tracking longitudinal and spatial clonal evolutionary history by next-generation sequencing , 2016, Proceedings of the National Academy of Sciences.

[29]  Dirk P. Kroese,et al.  Chromosome arm aneuploidies shape tumour evolution and drug response , 2020, Nature Communications.

[30]  A. Schäffer,et al.  Single-cell genetic analysis reveals insights into clonal development of prostate cancers and indicates loss of PTEN as a marker of poor prognosis. , 2014, The American journal of pathology.

[31]  Benjamin J. Raphael,et al.  Characterizing the allele- and haplotype-specific copy number landscape of cancer genomes at single-cell resolution with CHISEL , 2019 .

[32]  James D. Brenton,et al.  Phylogenetic Quantification of Intra-tumour Heterogeneity , 2013, PLoS Comput. Biol..

[33]  M. Blough,et al.  Comparative genomic and genetic analysis of glioblastoma-derived brain tumor-initiating cells and their parent tumors. , 2016, Neuro-oncology.

[34]  Russell Schwartz,et al.  Phylogenetic analysis of multiprobe fluorescence in situ hybridization data from tumor cell populations , 2013, Bioinform..

[35]  S. Gabriel,et al.  Pan-cancer patterns of somatic copy-number alteration , 2013, Nature Genetics.

[36]  Russell Schwartz,et al.  Robust unmixing of tumor states in array comparative genomic hybridization data , 2010, Bioinform..

[37]  Alberto Orfao,et al.  Detailed characterization of alterations of chromosomes 7, 9, and 10 in glioblastomas as assessed by single-nucleotide polymorphism arrays. , 2011, The Journal of molecular diagnostics : JMD.

[38]  R. Wilson,et al.  Cancer genome sequencing: a review. , 2009, Human molecular genetics.

[39]  Niko Beerenwinkel,et al.  BitPhylogeny: a probabilistic framework for reconstructing intra-tumor phylogenies , 2015, Genome Biology.