Application of single-cell RNA sequencing in optimizing a combinatorial therapeutic strategy in metastatic renal cell carcinoma

BackgroundIntratumoral heterogeneity hampers the success of marker-based anticancer treatment because the targeted therapy may eliminate a specific subpopulation of tumor cells while leaving others unharmed. Accordingly, a rational strategy minimizing survival of the drug-resistant subpopulation is essential to achieve long-term therapeutic efficacy.ResultsUsing single-cell RNA sequencing (RNA-seq), we examine the intratumoral heterogeneity of a pair of primary renal cell carcinoma and its lung metastasis. Activation of drug target pathways demonstrates considerable variability between the primary and metastatic sites, as well as among individual cancer cells within each site. Based on the prediction of multiple drug target pathway activation, we derive a combinatorial regimen co-targeting two mutually exclusive pathways for the metastatic cancer cells. This combinatorial strategy shows significant increase in the treatment efficacy over monotherapy in the experimental validation using patient-derived xenograft platforms in vitro and in vivo.ConclusionsOur findings demonstrate the investigational application of single-cell RNA-seq in the design of an anticancer regimen. The approach may overcome intratumoral heterogeneity which hampers the success of precision medicine.

[1]  Huanming Yang,et al.  Single-Cell Exome Sequencing Reveals Single-Nucleotide Mutation Characteristics of a Kidney Tumor , 2012, Cell.

[2]  Philippe Dessen,et al.  Characterization of a Large Panel of Patient-Derived Tumor Xenografts Representing the Clinical Heterogeneity of Human Colorectal Cancer , 2012, Clinical Cancer Research.

[3]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

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

[5]  Helga Thorvaldsdóttir,et al.  Integrative Genomics Viewer , 2011, Nature Biotechnology.

[6]  E. S. Venkatraman,et al.  A faster circular binary segmentation algorithm for the analysis of array CGH data , 2007, Bioinform..

[7]  L. Truong,et al.  PAX2 and PAX8 expression in primary and metastatic renal tumors: a comprehensive comparison. , 2012, Archives of pathology & laboratory medicine.

[8]  Jing Yang,et al.  Epithelial–mesenchymal plasticity in carcinoma metastasis , 2013, Genes & development.

[9]  G. Keating Afatinib: A Review of Its Use in the Treatment of Advanced Non-Small Cell Lung Cancer , 2014, Drugs.

[10]  Charles J. Vaske,et al.  Single-cell analyses of transcriptional heterogeneity during drug tolerance transition in cancer cells by RNA sequencing , 2014, Proceedings of the National Academy of Sciences.

[11]  N. McGranahan,et al.  The causes and consequences of genetic heterogeneity in cancer evolution , 2013, Nature.

[12]  C. Klein,et al.  Parallel progression of primary tumours and metastases , 2009, Nature Reviews Cancer.

[13]  J. Cheville,et al.  Clear Cell Renal Cell Carcinoma: Gene Expression Analyses Identify a Potential Signature for Tumor Aggressiveness , 2005, Clinical Cancer Research.

[14]  E. Lander,et al.  Erratum: Core epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudinlow and metaplastic breast cancer subtypes (Proceedings of the National Academy of Sciences of the United States of America (2010) 107:35 (15449-15454)) , 2010 .

[15]  S. Ramaswamy,et al.  Systematic identification of genomic markers of drug sensitivity in cancer cells , 2012, Nature.

[16]  Y. Okada,et al.  Invasion and metastasis of renal cell carcinoma , 2014, Medical Molecular Morphology.

[17]  R. Milo,et al.  Dynamic Proteomics of Individual Cancer Cells in Response to a Drug , 2008, Science.

[18]  Pornpimol Charoentong,et al.  ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks , 2009, Bioinform..

[19]  David McDermott,et al.  Temsirolimus, interferon alfa, or both for advanced renal-cell carcinoma. , 2007, The New England journal of medicine.

[20]  J. Montero,et al.  Inhibition of Src Family Kinases and Receptor Tyrosine Kinases by Dasatinib: Possible Combinations in Solid Tumors , 2011, Clinical Cancer Research.

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

[22]  J. Uhm Clonal selection drives genetic divergence of metastatic medulloblastoma , 2012 .

[23]  Richard Durbin,et al.  Fast and accurate long-read alignment with Burrows–Wheeler transform , 2010, Bioinform..

[24]  C. Charbonneau,et al.  Epidemiologic and socioeconomic burden of metastatic renal cell carcinoma (mRCC): a literature review. , 2008, Cancer treatment reviews.

[25]  F. Erdoğan,et al.  Prognostic significance of morphologic parameters in renal cell carcinoma , 2004, International journal of clinical practice.

[26]  M. Higashiyama,et al.  Differences in chemosensitivity between primary and paired metastatic lung cancer tissues: In vitro analysis based on the collagen gel droplet embedded culture drug test (CD-DST). , 2012, Journal of thoracic disease.

[27]  G. Hannon,et al.  Patient-derived tumor xenografts: transforming clinical samples into mouse models. , 2013, Cancer research.

[28]  N. Cox,et al.  Clinical drug response can be predicted using baseline gene expression levels and in vitro drug sensitivity in cell lines , 2014, Genome Biology.

[29]  Mark T. W. Ebbert,et al.  Tumor grafts derived from women with breast cancer authentically reflect tumor pathology, growth, metastasis and disease outcomes , 2011, Nature Medicine.

[30]  Victor Raskin,et al.  School of Mechanical Engineering , 2007 .

[31]  Pablo Cingolani,et al.  © 2012 Landes Bioscience. Do not distribute. , 2022 .

[32]  T. Choueiri,et al.  External validation and comparison with other models of the International Metastatic Renal-Cell Carcinoma Database Consortium prognostic model: a population-based study. , 2013, The Lancet. Oncology.

[33]  Y. Nagashima,et al.  Contribution of the Src family of kinases to the appearance of malignant phenotypes in renal cancer cells , 2005, Molecular carcinogenesis.

[34]  L. Gunaratnam,et al.  Silencing of epidermal growth factor receptor suppresses hypoxia-inducible factor-2-driven VHL-/- renal cancer. , 2005, Cancer research.

[35]  T. Hudson,et al.  The Genetic Basis for Cancer Treatment Decisions , 2012, Cell.

[36]  Thomas D. Wu,et al.  A comprehensive transcriptional portrait of human cancer cell lines , 2014, Nature Biotechnology.

[37]  G. Getz,et al.  Inferring tumour purity and stromal and immune cell admixture from expression data , 2013, Nature Communications.

[38]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[39]  F. Waldman,et al.  Epidermal growth factor receptor expression is associated with rapid tumor cell proliferation in renal cell carcinoma. , 1997, Human pathology.

[40]  Nci Dream Community A community effort to assess and improve drug sensitivity prediction algorithms , 2014 .

[41]  Steven J. M. Jones,et al.  Comprehensive molecular characterization of clear cell renal cell carcinoma , 2013, Nature.

[42]  Laura M. Heiser,et al.  A community effort to assess and improve drug sensitivity prediction algorithms , 2014, Nature Biotechnology.

[43]  P. Sorger,et al.  Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis , 2009, Nature.

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

[45]  Cheng Li,et al.  Adjusting batch effects in microarray expression data using empirical Bayes methods. , 2007, Biostatistics.

[46]  Marie Joseph,et al.  Gene Signatures of Progression and Metastasis in Renal Cell Cancer , 2005, Clinical Cancer Research.

[47]  Juergen Friedrich,et al.  Spheroid-based drug screen: considerations and practical approach , 2009, Nature Protocols.

[48]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

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

[50]  S. Pomer,et al.  Concomitant overexpression of the EGFR and erbB‐2 genes in renal cell carcinoma (RCC) is correlated with dedifferentiation and metastasis , 1996, International journal of cancer.

[51]  H. Seol,et al.  Patient-Derived Xenografts from Non–Small Cell Lung Cancer Brain Metastases Are Valuable Translational Platforms for the Development of Personalized Targeted Therapy , 2014, Clinical Cancer Research.

[52]  Shawn M. Gillespie,et al.  Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma , 2014, Science.

[53]  N. Navin Cancer genomics: one cell at a time , 2014, Genome Biology.

[54]  P. A. Futreal,et al.  Genomic architecture and evolution of clear cell renal cell carcinomas defined by multiregion sequencing , 2014, Nature Genetics.

[55]  J. Yokota,et al.  Molecular footprints of human lung cancer progression , 2004, Cancer science.

[56]  Michael L. Gatza,et al.  A pathway-based classification of human breast cancer , 2010, Proceedings of the National Academy of Sciences.

[57]  Sang Hyun Yi,et al.  High-throughput screening (HTS) of anticancer drug efficacy on a micropillar/microwell chip platform. , 2014, Analytical chemistry.

[58]  R. Kamm,et al.  Microfluidic assay for simultaneous culture of multiple cell types on surfaces or within hydrogels , 2012, Nature Protocols.

[59]  Justin Guinney,et al.  GSVA: gene set variation analysis for microarray and RNA-Seq data , 2013, BMC Bioinformatics.

[60]  R. Arceci Intratumor Heterogeneity and Branched Evolution Revealed by Multiregion Sequencing , 2012 .

[61]  R. Motzer,et al.  Comparative effectiveness of axitinib versus sorafenib in advanced renal cell carcinoma (AXIS): a randomised phase 3 trial , 2011, The Lancet.

[62]  W. Linehan,et al.  The genetic basis of kidney cancer: a metabolic disease , 2010, Nature Reviews Urology.

[63]  Jordan C. Ciciliano,et al.  Corrigendum: RNA sequencing of pancreatic circulating tumour cells implicates WNT signalling in metastasis , 2012, Nature.

[64]  Joshua F. McMichael,et al.  Genome Remodeling in a Basal-like Breast Cancer Metastasis and Xenograft , 2010, Nature.

[65]  Yang Zhao,et al.  EMT Phenotype is Induced by Increased Src Kinase Activity via Src-Mediated Caspase-8 Phosphorylation , 2012, Cellular Physiology and Biochemistry.

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

[67]  T. Pawson,et al.  Effects of dasatinib on EphA2 receptor tyrosine kinase activity and downstream signalling in pancreatic cancer , 2008, British Journal of Cancer.

[68]  Steven J. M. Jones,et al.  Comprehensive molecular characterization of clear cell renal cell carcinoma , 2013, Nature.

[69]  Andrew M. K. Brown,et al.  Variable Clonal Repopulation Dynamics Influence Chemotherapy Response in Colorectal Cancer , 2013, Science.

[70]  D. Trump Lapatinib versus hormone therapy in patients with advanced renal cell carcinoma: A randomized Phase III clinical trial , 2008 .

[71]  Do-Hyun Nam,et al.  Single-cell mRNA sequencing identifies subclonal heterogeneity in anti-cancer drug responses of lung adenocarcinoma cells , 2015, Genome Biology.

[72]  H. Moch,et al.  Brain metastasis in renal cancer patients: metastatic pattern, tumour-associated macrophages and chemokine/chemoreceptor expression , 2013, British Journal of Cancer.

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

[74]  Michael B Atkins,et al.  Resistance to targeted therapy in renal-cell carcinoma. , 2009, The Lancet. Oncology.