Deciphering the Evolutionary History of Complex Rearrangements in Head and Neck Cancer Patients Using Multi-Omic Approach

Despite the large efforts in international cancer genome consortium studies, there are still a large proportion of tumors with complex genomic rearrangement often remained without a clinically relevant molecular characterization. Integration of multi-omic data helps elucidating evolutionary history of such cases and identifying predictive molecular markers. Here we present the findings of our proof-of-principle study that investigated the evolutionary history of complex rearrangements in primary head and neck tumor genomes integrating long-read whole-genome, Hi-C, and RNA sequencing. We report a HPV-positive case with development of complex genomic rearrangements tracing back to HPV-mediated genomic instability and a HPV-negative case with an enhancer hi-jacking in a region of chromothripsis predicted to co-occur with a neoloop and a super-enhancer. These structural alterations resulted in overexpression of the oncogenes CCND1 and ALK, respectively, validated with immunohistochemistry assay. Furthermore, we introduce a novel analytic approach utilizing long-read whole-genome data distinguishing somatic mutations before and after structural variants. Our findings highlight the need for multi-modal sequencing strategies to increase our understanding of cancer evolution and rare biomarkers in poorly understood cancers.

[1]  Aaron M. Streets,et al.  The complete sequence of a human genome , 2021, bioRxiv.

[2]  Rachel M. Sherman,et al.  Comprehensive analysis of structural variants in breast cancer genomes using single-molecule sequencing , 2019, bioRxiv.

[3]  Sergey Koren,et al.  Accurate circular consensus long-read sequencing improves variant detection and assembly of a human genome , 2019, Nature Biotechnology.

[4]  S. Galiacy,et al.  Targeting ALK in Cancer: Therapeutic Potential of Proapoptotic Peptides , 2019, Cancers.

[5]  Thomas Colthurst,et al.  A universal SNP and small-indel variant caller using deep neural networks , 2018, Nature Biotechnology.

[6]  S. McWeeney,et al.  Induction of anaplastic lymphoma kinase (ALK) as a novel mechanism of EGFR inhibitor resistance in head and neck squamous cell carcinoma patient-derived models , 2018, Cancer biology & therapy.

[7]  S. Mundlos,et al.  Structural variation in the 3D genome , 2018, Nature Reviews Genetics.

[8]  P. A. Futreal,et al.  Timing the Landmark Events in the Evolution of Clear Cell Renal Cell Cancer: TRACERx Renal , 2018, Cell.

[9]  Kin Chung Lam,et al.  High-resolution TADs reveal DNA sequences underlying genome organization in flies , 2017, Nature Communications.

[10]  Heng Li,et al.  Minimap2: pairwise alignment for nucleotide sequences , 2017, Bioinform..

[11]  Benjamin J. Raphael,et al.  The evolutionary history of 2,658 cancers , 2017, Nature.

[12]  Michael C. Schatz,et al.  Accurate detection of complex structural variations using single molecule sequencing , 2017, Nature Methods.

[13]  Nicolai J. Birkbak,et al.  Tracking the Evolution of Non‐Small‐Cell Lung Cancer , 2017, The New England journal of medicine.

[14]  Jacob M. Luber,et al.  HiGlass: web-based visual exploration and analysis of genome interaction maps , 2017, Genome Biology.

[15]  A. Shaw,et al.  Targeting ALK: Precision Medicine Takes on Drug Resistance. , 2017, Cancer discovery.

[16]  Md Abu Shufean,et al.  ALK: a tyrosine kinase target for cancer therapy , 2017, Cold Spring Harbor molecular case studies.

[17]  K. Trinkaus,et al.  Phase I trial of palbociclib, a selective cyclin dependent kinase 4/6 inhibitor, in combination with cetuximab in patients with recurrent/metastatic head and neck squamous cell carcinoma. , 2016, Oral oncology.

[18]  Neva C. Durand,et al.  A 3D Map of the Human Genome at Kilobase Resolution Reveals Principles of Chromatin Looping , 2014, Cell.

[19]  Heng Li Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM , 2013, 1303.3997.

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

[21]  Yuh-Seog Jung,et al.  Human papillomavirus in oropharyngeal squamous cell carcinomas in korea: Use of G1 cycle markers as new prognosticators , 2012, Head & neck.

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

[23]  A. Børresen-Dale,et al.  Mutational Processes Molding the Genomes of 21 Breast Cancers , 2012, Cell.

[24]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[25]  Colin N. Dewey,et al.  RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome , 2011, BMC Bioinformatics.

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

[27]  Cary A Moody,et al.  Human papillomavirus oncoproteins: pathways to transformation , 2010, Nature Reviews Cancer.

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

[29]  H. Katayama,et al.  A genetic variant of Aurora kinase A promotes genomic instability leading to highly malignant skin tumors. , 2009, Cancer research.

[30]  F. Gergely,et al.  Aurora-A: the maker and breaker of spindle poles , 2007, Journal of Cell Science.