Interpretable systems biomarkers predict response to immune-checkpoint inhibitors

Summary Cancer cells can leverage several cell-intrinsic and -extrinsic mechanisms to escape immune system recognition. The inherent complexity of the tumor microenvironment, with its multicellular and dynamic nature, poses great challenges for the extraction of biomarkers of immune response and immunotherapy efficacy. Here, we use RNA-sequencing (RNA-seq) data combined with different sources of prior knowledge to derive system-based signatures of the tumor microenvironment, quantifying immune-cell composition and intra- and intercellular communications. We applied multi-task learning to these signatures to predict different hallmarks of immune responses and derive cancer-type-specific models based on interpretable systems biomarkers. By applying our models to independent RNA-seq data from cancer patients treated with PD-1/PD-L1 inhibitors, we demonstrated that our method to Estimate Systems Immune Response (EaSIeR) accurately predicts therapeutic outcome. We anticipate that EaSIeR will be a valuable tool to provide a holistic description of immune responses in complex and dynamic systems such as tumors using available RNA-seq data.

[1]  G. Freeman,et al.  Interferon-γ-induced activation of JAK1 and JAK2 suppresses tumor cell susceptibility to NK cells through upregulation of PD-L1 expression , 2015, Oncoimmunology.

[2]  J. Christensen,et al.  Cytotoxic CD8+ T cells in cancer and cancer immunotherapy , 2020, British Journal of Cancer.

[3]  P. Sharma,et al.  Immune Checkpoint Targeting in Cancer Therapy: Toward Combination Strategies with Curative Potential , 2015, Cell.

[4]  P. Blackshear,et al.  Oncogenic RAS Signaling Promotes Tumor Immunoresistance by Stabilizing PD-L1 mRNA , 2017, Immunity.

[5]  G. Berchem,et al.  Improving Cancer Immunotherapy by Targeting the Hypoxic Tumor Microenvironment: New Opportunities and Challenges , 2019, Cells.

[6]  Eytan Ruppin,et al.  Robust prediction of response to immune checkpoint blockade therapy in metastatic melanoma , 2018, Nature Medicine.

[7]  M. Park,et al.  Irradiation of breast cancer cells enhances CXCL16 ligand expression and induces the migration of natural killer cells expressing the CXCR6 receptor. , 2016, Cytotherapy.

[8]  Christian H. Holland,et al.  Benchmark and integration of resources for the estimation of human transcription factor activities. , 2019, Genome research.

[9]  M. Reina,et al.  Role of LFA-1 and ICAM-1 in Cancer , 2017, Cancers.

[10]  P. Hersey,et al.  Inducible but Not Constitutive Expression of PD-L1 in Human Melanoma Cells Is Dependent on Activation of NF-κB , 2015, PloS one.

[11]  Hideo Baba,et al.  CXCL9, CXCL10, CXCL11/CXCR3 axis for immune activation - A target for novel cancer therapy. , 2018, Cancer treatment reviews.

[12]  D. Olive,et al.  B and T Lymphocyte Attenuator Is Highly Expressed on CMV-Specific T Cells during Infection and Regulates Their Function , 2010, The Journal of Immunology.

[13]  Monika S. Kowalczyk,et al.  A Cancer Cell Program Promotes T Cell Exclusion and Resistance to Checkpoint Blockade , 2018, Cell.

[14]  Alexander W. Jung,et al.  Pan-cancer computational histopathology reveals mutations, tumor composition and prognosis , 2019, Nature Cancer.

[15]  O. Elemento,et al.  Cancer systems biology: embracing complexity to develop better anticancer therapeutic strategies , 2014, Oncogene.

[16]  D. Klein The Tumor Vascular Endothelium as Decision Maker in Cancer Therapy , 2018, Front. Oncol..

[17]  Federica Eduati,et al.  Toward Systems Biomarkers of Response to Immune Checkpoint Blockers , 2020, Frontiers in Oncology.

[18]  V. Sondak,et al.  12-Chemokine Gene Signature Identifies Lymph Node-like Structures in Melanoma: Potential for Patient Selection for Immunotherapy? , 2012, Scientific Reports.

[19]  H. Young,et al.  Human B Cell Activation by Autologous NK Cells Is Regulated by CD40-CD40 Ligand Interaction: Role of Memory B Cells and CD5+ B Cells , 2001, The Journal of Immunology.

[20]  Zhuo Chen,et al.  Establishment and validation of an immune-based prognostic score model in glioblastoma. , 2020, International immunopharmacology.

[21]  D. Jackson,et al.  Cooperation between Constitutive and Inducible Chemokines Enables T Cell Engraftment and Immune Attack in Solid Tumors. , 2019, Cancer cell.

[22]  B. Becher,et al.  BAFF-secreting neutrophils drive plasma cell responses during emergency granulopoiesis , 2016, The Journal of experimental medicine.

[23]  Pornpimol Charoentong,et al.  Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade , 2016, bioRxiv.

[24]  Matthew D. Hellmann,et al.  Immune‐Related Adverse Events Associated with Immune Checkpoint Blockade , 2018, The New England journal of medicine.

[25]  J. Lunceford,et al.  IFN-g – related mRNA profile predicts clinical response to PD-1 blockade , 2019 .

[26]  Mark D. M. Leiserson,et al.  Precision Oncology: The Road Ahead. , 2017, Trends in molecular medicine.

[27]  Lior Pachter,et al.  Near-optimal probabilistic RNA-seq quantification , 2016, Nature Biotechnology.

[28]  Laurence Zitvogel,et al.  The immune contexture in cancer prognosis and treatment , 2017, Nature Reviews Clinical Oncology.

[29]  C. Robert A decade of immune-checkpoint inhibitors in cancer therapy , 2020, Nature Communications.

[30]  Koichi S. Kobayashi,et al.  NLRC5/CITA: A Key Player in Cancer Immune Surveillance. , 2017, Trends in cancer.

[31]  S. Bonini,et al.  Who benefits from immunotherapy? , 1997, Allergy.

[32]  J. Carpten,et al.  Translating RNA sequencing into clinical diagnostics: opportunities and challenges , 2016, Nature Reviews Genetics.

[33]  A. Regev,et al.  Integrative molecular and clinical modeling of clinical outcomes to PD1 blockade in patients with metastatic melanoma , 2019, Nature Medicine.

[34]  J. Wargo,et al.  Primary, Adaptive, and Acquired Resistance to Cancer Immunotherapy , 2017, Cell.

[35]  D. Aldinucci,et al.  The CCL5/CCR5 Axis in Cancer Progression , 2020, Cancers.

[36]  J. Lunceford,et al.  IFN-&ggr;–related mRNA profile predicts clinical response to PD-1 blockade , 2017, The Journal of clinical investigation.

[37]  Sun Wook Cho,et al.  CXCL16 positively correlated with M2-macrophage infiltration, enhanced angiogenesis, and poor prognosis in thyroid cancer , 2019, Scientific Reports.

[38]  Vanessa M. Hubbard-Lucey,et al.  Trends in the global immuno-oncology landscape , 2018, Nature Reviews Drug Discovery.

[39]  Timothy A. Chan,et al.  The hallmarks of successful anticancer immunotherapy , 2018, Science Translational Medicine.

[40]  Inference and analysis of cell-cell communication using CellChat , 2021, Nature communications.

[41]  J. Modiano,et al.  Fas ligand based immunotherapy: A potent and effective neoadjuvant with checkpoint inhibitor properties, or a systemically toxic promoter of tumor growth? , 2016, Discovery medicine.

[42]  Joon-Oh Park,et al.  Comprehensive molecular characterization of clinical responses to PD-1 inhibition in metastatic gastric cancer , 2018, Nature Medicine.

[43]  Peter S. Linsley,et al.  Copy Number Loss of the Interferon Gene Cluster in Melanomas Is Linked to Reduced T Cell Infiltrate and Poor Patient Prognosis , 2014, PloS one.

[44]  Matthieu Texier,et al.  Safety profiles of anti-CTLA-4 and anti-PD-1 antibodies alone and in combination , 2016, Nature Reviews Clinical Oncology.

[45]  Jeffrey E. Lee,et al.  B cells and tertiary lymphoid structures promote immunotherapy response , 2020, Nature.

[46]  Annelisa M Cornel,et al.  MHC Class I Downregulation in Cancer: Underlying Mechanisms and Potential Targets for Cancer Immunotherapy , 2020, Cancers.

[47]  G. Korbutt,et al.  Mannose 6-Phosphate/Insulin-like Growth Factor II Receptor Is a Death Receptor for Granzyme B during Cytotoxic T Cell–Induced Apoptosis , 2000, Cell.

[48]  T. Chan,et al.  Tumor and Microenvironment Evolution during Immunotherapy with Nivolumab , 2017, Cell.

[49]  J. Lasota,et al.  Clinicopathologic profile, immunophenotype, and genotype of CD274 (PD-L1)-positive colorectal carcinomas , 2017, Modern Pathology.

[50]  P. Mischel,et al.  Loss of tumor suppressor PTEN function increases B7-H1 expression and immunoresistance in glioma , 2007, Nature Medicine.

[51]  J. Bergh,et al.  Genetic, transcriptional and post-translational regulation of the programmed death protein ligand 1 in cancer: biology and clinical correlations , 2018, Oncogene.

[52]  J. Hadfield,et al.  RNA sequencing: the teenage years , 2019, Nature Reviews Genetics.

[53]  Liping Wang,et al.  Roles of IFN-γ in tumor progression and regression: a review , 2020, Biomarker Research.

[54]  D. Ray,et al.  Molecular Actions of PPARα in Lipid Metabolism and Inflammation. , 2018, Endocrine reviews.

[55]  Ou Sha,et al.  Trichosanthin increases Granzyme B penetration into tumor cells by upregulation of CI-MPR on the cell surface , 2017, Oncotarget.

[56]  J. Sáez-Rodríguez,et al.  Perturbation-response genes reveal signaling footprints in cancer gene expression , 2016, Nature Communications.

[57]  M. Socinski,et al.  First‐Line Nivolumab in Stage IV or Recurrent Non–Small‐Cell Lung Cancer , 2017, The New England journal of medicine.

[58]  C. Ager,et al.  Targeted hypoxia reduction restores T cell infiltration and sensitizes prostate cancer to immunotherapy , 2018, The Journal of clinical investigation.

[59]  Nuno A. Fonseca,et al.  Transcription Factor Activities Enhance Markers of Drug Sensitivity in Cancer. , 2018, Cancer research.

[60]  Jian Su,et al.  Potential Predictive Value of TP53 and KRAS Mutation Status for Response to PD-1 Blockade Immunotherapy in Lung Adenocarcinoma , 2016, Clinical Cancer Research.

[61]  O. Lantz,et al.  Anti-NKG2A mAb Is a Checkpoint Inhibitor that Promotes Anti-tumor Immunity by Unleashing Both T and NK Cells , 2018, Cell.

[62]  J. Galon,et al.  Approaches to treat immune hot, altered and cold tumours with combination immunotherapies , 2019, Nature Reviews Drug Discovery.

[63]  J. Borst,et al.  Targeting the T-cell co-stimulatory CD27/CD70 pathway in cancer immunotherapy: rationale and potential. , 2015, Immunotherapy.

[64]  C. Sawyers,et al.  The phosphatidylinositol 3-Kinase–AKT pathway in human cancer , 2002, Nature Reviews Cancer.

[65]  G. Freeman,et al.  Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response , 2018, Nature Medicine.

[66]  Adam A. Margolin,et al.  The Cancer Cell Line Encyclopedia enables predictive modeling of anticancer drug sensitivity , 2012, Nature.

[67]  O. Harismendy,et al.  Deciphering cell–cell interactions and communication from gene expression , 2020, Nature reviews. Genetics.

[68]  Yu Wang,et al.  Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data , 2017, Genome Medicine.

[69]  Rajarsi R. Gupta,et al.  Spatial Organization and Molecular Correlation of Tumor-Infiltrating Lymphocytes Using Deep Learning on Pathology Images. , 2018, Cell reports.

[70]  F. Eduati,et al.  Multi-Omics Profiling of the Tumor Microenvironment: Paving the Way to Precision Immuno-Oncology , 2018, Front. Oncol..

[71]  Mirjana Efremova,et al.  CellPhoneDB: inferring cell–cell communication from combined expression of multi-subunit ligand–receptor complexes , 2020, Nature Protocols.

[72]  S. Elledge,et al.  Tumor aneuploidy correlates with markers of immune evasion and with reduced response to immunotherapy , 2017, Science.

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

[74]  K. Tsuneyama,et al.  High-level expression of chemokine CXCL16 by tumor cells correlates with a good prognosis and increased tumor-infiltrating lymphocytes in colorectal cancer. , 2007, Cancer research.

[75]  Lieping Chen,et al.  Interferon regulatory factor‐1 is prerequisite to the constitutive expression and IFN‐γ‐induced upregulation of B7‐H1 (CD274) , 2006, FEBS letters.

[76]  Jing-mei Liu,et al.  CXCL9: evidence and contradictions for its role in tumor progression , 2016, Cancer medicine.

[77]  Bertram Klinger,et al.  Perturbation-response genes reveal signaling footprints in cancer gene expression , 2016 .

[78]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[79]  Tae Min Kim,et al.  Genomic landscape associated with potential response to anti-CTLA-4 treatment in cancers , 2017, Nature Communications.

[80]  Jun Li,et al.  TCPA: a resource for cancer functional proteomics data , 2013, Nature Methods.

[81]  A. Mantovani,et al.  Antagonistic Inflammatory Phenotypes Dictate Tumor Fate and Response to Immune Checkpoint Blockade , 2020, Immunity.

[82]  Y. M. Lee,et al.  Emerging role of RUNX3 in the regulation of tumor microenvironment , 2018, BMB reports.

[83]  N. Hacohen,et al.  Molecular and Genetic Properties of Tumors Associated with Local Immune Cytolytic Activity , 2015, Cell.

[84]  I. Mellman,et al.  Oncology meets immunology: the cancer-immunity cycle. , 2013, Immunity.

[85]  D. Schadendorf,et al.  Tertiary lymphoid structures improve immunotherapy and survival in melanoma , 2020, Nature.

[86]  Piero Carninci,et al.  A draft network of ligand–receptor-mediated multicellular signalling in human , 2015, Nature Communications.

[87]  Delin Liu,et al.  Immunotherapy strategy of EGFR mutant lung cancer. , 2018, American journal of cancer research.

[88]  Christian H. Holland,et al.  Robustness and applicability of transcription factor and pathway analysis tools on single-cell RNA-seq data , 2020, Genome Biology.

[89]  Weiping Zou,et al.  Chemokines in the cancer microenvironment and their relevance in cancer immunotherapy , 2017, Nature Reviews Immunology.

[90]  Wei Huang,et al.  Targeting STAT3 in Cancer Immunotherapy , 2020, Molecular cancer.

[91]  P. A. Futreal,et al.  Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance , 2017, Science Translational Medicine.

[92]  Mehmet Gönen,et al.  Bayesian Efficient Multiple Kernel Learning , 2012, ICML.

[93]  Stacey D. Finley,et al.  Integrative Approaches to Cancer Immunotherapy. , 2019, Trends in cancer.

[94]  W. Gu,et al.  PD-L1 Distribution and Perspective for Cancer Immunotherapy—Blockade, Knockdown, or Inhibition , 2019, Front. Immunol..

[95]  M. Yuan,et al.  Model selection and estimation in regression with grouped variables , 2006 .

[96]  Z. Trajanoski,et al.  Quantifying tumor-infiltrating immune cells from transcriptomics data , 2018, Cancer Immunology, Immunotherapy.

[97]  Charles Schmidt The benefits of immunotherapy combinations. , 2017 .

[98]  G. Guarda,et al.  NLRC5, a promising new entry in tumor immunology , 2016, Journal of Immunotherapy for Cancer.

[99]  J. Madore,et al.  Distinct Immune Cell Populations Define Response to Anti-PD-1 Monotherapy and Anti-PD-1/Anti-CTLA-4 Combined Therapy. , 2019, Cancer cell.

[100]  Christian H. Holland,et al.  Transfer of regulatory knowledge from human to mouse for functional genomic analysis. , 2020, Biochimica et biophysica acta. Gene regulatory mechanisms.

[101]  Z. Trajanoski,et al.  Next-generation computational tools for interrogating cancer immunity , 2019, Nature Reviews Genetics.

[102]  I. Mellman,et al.  Elements of cancer immunity and the cancer–immune set point , 2017, Nature.

[103]  R. Homer,et al.  Inhibition of pulmonary fibrosis in mice by CXCL10 requires glycosaminoglycan binding and syndecan-4. , 2010, The Journal of clinical investigation.

[104]  E. Puré,et al.  Cancer-associated fibroblasts: key determinants of tumor immunity and immunotherapy , 2020, Current opinion in immunology.

[105]  T. Graeber,et al.  Interferon Receptor Signaling Pathways Regulating PD-L1 and PD-L2 Expression. , 2019, Cell reports.

[106]  F. Garrido,et al.  Adenovirus expressing β2-microglobulin recovers HLA class I expression and antitumor immunity by increasing T-cell recognition , 2014, Cancer Gene Therapy.

[107]  Yunyan Gu,et al.  A qualitative transcriptional signature for predicting microsatellite instability status of right-sided Colon Cancer , 2019, BMC genomics.

[108]  Joshua M. Stuart,et al.  The Cancer Genome Atlas Pan-Cancer analysis project , 2013, Nature Genetics.

[109]  Mark J. Smyth,et al.  The TRAIL apoptotic pathway in cancer onset, progression and therapy , 2008, Nature Reviews Cancer.

[110]  M. Shimada,et al.  High STAT4 Expression Indicates Better Disease-free Survival in Patients with Gastric Cancer. , 2017, Anticancer research.

[111]  E. Smit,et al.  The Potential of Combined Immunotherapy and Antiangiogenesis for the Synergistic Treatment of Advanced NSCLC , 2017, Journal of thoracic oncology : official publication of the International Association for the Study of Lung Cancer.

[112]  P. Wen,et al.  Neoadjuvant anti-PD-1 immunotherapy promotes a survival benefit with intratumoral and systemic immune responses in recurrent glioblastoma , 2018, Nature Medicine.

[113]  C. Caux,et al.  Cold Tumors: A Therapeutic Challenge for Immunotherapy , 2019, Front. Immunol..

[114]  Yan Guo,et al.  A Cell-Based Systems Biology Assessment of Human Blood to Monitor Immune Responses after Influenza Vaccination , 2015, PloS one.

[115]  Trevor Hastie,et al.  Regularization Paths for Generalized Linear Models via Coordinate Descent. , 2010, Journal of statistical software.

[116]  H. Zou,et al.  Regularization and variable selection via the elastic net , 2005 .

[117]  F. Pontén,et al.  CDX2: A Prognostic Marker in Metastatic Colorectal Cancer Defining a Better BRAF Mutated and a Worse KRAS Mutated Subgroup , 2020, Frontiers in Oncology.

[118]  Douglas A. Lauffenburger,et al.  Analysis of Single-Cell RNA-Seq Identifies Cell-Cell Communication Associated with Tumor Characteristics , 2018, Cell reports.

[119]  F. Marincola,et al.  BTLA marks a less-differentiated tumor-infiltrating lymphocyte subset in melanoma with enhanced survival properties , 2015, Oncoimmunology.

[120]  I. Melero,et al.  ICAM-1-LFA-1 Dependent CD8+ T-Lymphocyte Aggregation in Tumor Tissue Prevents Recirculation to Draining Lymph Nodes , 2018, Front. Immunol..

[121]  Greta Volpedo,et al.  Immune Suppression Mediated by STAT4 Deficiency Promotes Lymphatic Metastasis in HNSCC , 2020, Frontiers in Immunology.

[122]  R. Bourgon,et al.  TGF-β attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells , 2018, Nature.

[123]  J. Wargo,et al.  Loss of IFN-γ Pathway Genes in Tumor Cells as a Mechanism of Resistance to Anti-CTLA-4 Therapy , 2016, Cell.

[124]  D. Pardoll,et al.  Mechanisms regulating PD-L1 expression on tumor and immune cells , 2019, Journal of Immunotherapy for Cancer.

[125]  T. Gajewski,et al.  Impact of oncogenic pathways on evasion of antitumour immune responses , 2018, Nature Reviews Cancer.