Heterogeneous Tumor-Immune Microenvironments among Differentially Growing Metastases in an Ovarian Cancer Patient - PubMed (original) (raw)
Case Reports
. 2017 Aug 24;170(5):927-938.e20.
doi: 10.1016/j.cell.2017.07.025.
Danish Memon 2, Stephane Pourpe 3, Harini Veeraraghavan 4, Yanyun Li 5, Hebert Alberto Vargas 6, Michael B Gill 1, Kay J Park 7, Oliver Zivanovic 8, Jason Konner 9, Jacob Ricca 5, Dmitriy Zamarin 10, Tyler Walther 3, Carol Aghajanian 9, Jedd D Wolchok 11, Evis Sala 6, Taha Merghoub 5, Alexandra Snyder 12, Martin L Miller 13
Affiliations
- PMID: 28841418
- PMCID: PMC5589211
- DOI: 10.1016/j.cell.2017.07.025
Case Reports
Heterogeneous Tumor-Immune Microenvironments among Differentially Growing Metastases in an Ovarian Cancer Patient
Alejandro Jiménez-Sánchez et al. Cell. 2017.
Abstract
We present an exceptional case of a patient with high-grade serous ovarian cancer, treated with multiple chemotherapy regimens, who exhibited regression of some metastatic lesions with concomitant progression of other lesions during a treatment-free period. Using immunogenomic approaches, we found that progressing metastases were characterized by immune cell exclusion, whereas regressing and stable metastases were infiltrated by CD8+ and CD4+ T cells and exhibited oligoclonal expansion of specific T cell subsets. We also detected CD8+ T cell reactivity against predicted neoepitopes after isolation of cells from a blood sample taken almost 3 years after the tumors were resected. These findings suggest that multiple distinct tumor immune microenvironments co-exist within a single individual and may explain in part the heterogeneous fates of metastatic lesions often observed in the clinic post-therapy. VIDEO ABSTRACT.
Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.
Figures
Graphical abstract
Figure 1
Metastatic Tumors Exhibit Heterogeneous Growth and Somatic Mutation Patterns after Multi-line Chemotherapy (A) Representative CT scans showing concomitant progression/regression of the different resected metastatic tumors. RUQ = right upper quadrant. “Spleen” refers to the tumor deposit adjacent to the spleen. (B) CT-based volume of metastatic lesions represented with the solid vertical lines and dynamics of quantified CA125 levels with the red line indicating the CA125 upper limit of normal (35 units/ml). The x axis at the bottom shows a timeline of therapeutic interventions and clinical follow up. (C) Number of missense, silent, and nonsense mutations. (D) The phylogenetic tree represents the relationship of the samples based on binary calls of non-silent point mutations (Table S1A). Length of the branches is proportional to the number of mutations. Potential driver mutations are indicated. (E) Hierarchical cluster analysis (Euclidean distance metric and “average” linkage method) of the cellular prevalence of point mutations (n = 299) estimated with PyClone (Roth et al., 2014) (Table S1B).
Figure S1
Non-Silent Somatic Mutations and Copy-Number Alterations, Related to Figures 1 and S2 and Table S1 (A) Binary matrix of present/absent non-silent point mutations (n = 188) used for the phylogeny tree reconstruction in Figure 1D (Table S1A). (B) Relative copy-number alterations inferred from WES data of the primary and metastatic samples using CopywriteR (Kuilman et al., 2015). (C) Relative copy number profiles and tumor purity inferred after ABSOLUTE (Carter et al., 2012) analysis. Amplified and deep deleted DNA segments were defined as copy number alterations with at least ± 2 median absolute deviations for each sample. MAD = median absolute deviation.
Figure 2
Differential Expression of Immune-Related Pathways in Heterogeneously Growing Tumors (A) Expression levels and genetic alterations of genes associated with chemotherapy resistance in HGSOC (Patch et al., 2015) and multidrug resistance. Amplification and deep deletion were defined as at least ± 2 median absolute deviations of copy-number alterations for each sample (Figure S1C). (B) Single-sample gene set enrichment analysis (Barbie et al., 2009, Subramanian et al., 2005) of upregulated pathways using the KEGG (Kanehisa and Goto, 2000, Kanehisa et al., 2016) and REACTOME (Fabregat et al., 2016) databases (Tables S2D and S2E). Significantly enriched pathways (q < 0.05) with at least ± 1 log2 change relative to the median of the other samples are colored (Table S2G). False-discovery rate adjusted p value (q value) was calculated using the Benjamini-Hochberg method.
Figure S2
Gene Set Analysis of Transcript Abundance and Somatic Alteration Patterns across Samples, Related to Figure 2 and Table S2 (A–C) Gene-expression levels and genetic alterations of the DNA damage, apoptosis pathways, and caspases. (D) Expression levels of the 50 most variant genes according to their coefficient of variation (Table S2A). (E) Differential enrichment scores and enrichment q values of downregulated pathways between tumor samples (Tables S2D and S2E). No significantly enriched pathways (q < 0.05) with at least ± 1 log2 change relative to the median of the other samples were detected (Table S2G). False-discovery rate adjusted p value (q value) was calculated using the Benjamini-Hochberg method.
Figure 3
Immune Infiltration Status Shows Heterogeneous Microenvironments across Tumor Samples (A) Tumor purity and immune component estimated by analyzing Affymetrix-based transcriptomics (Table S3A) (Yoshihara et al., 2013). (B) Fractions of immune cell subsets in tumor samples inferred from gene-expression data using CIBERSORT (Newman et al., 2015). Width of bars is proportional to the −log10 p value of the deconvolution (Table S3B). CIBERSORT empirical p value, ∗p < 0.05. (C) Representative images of hematoxylin and eosin staining of tumor samples and immunofluorescence staining for DAPI, cytotoxic T cells (CD8+), helper T cells (CD4+FOXP3−), T cells (CD3+), T-regs (CD4+FOXP3+), macrophages (CD68+), and immune-checkpoint PD-L1. Complete slides are shown in Figure S3. (D) Image-based cell quantification of whole slides (Table S3C).
Figure S3
Complete Slide Hematoxylin and Eosin and Immunofluorescent Staining, Related to Figure 3 and Table S3 Hematoxylin and eosin staining of tumor samples. Immunofluorescence staining for cytotoxic T cells (CD8+), helper T cells (CD4+FOXP3−), and regulatory T cells (CD4+FOXP3+).
Figure 4
Higher HLA Expression and T Cell Oligoclonal Expansion Detected in Regressing Tumors (A) HLA-I and II gene differential expression across samples (Table S2A). (B) Number of predicted neoepitopes per sample (Tables S4B–S4D). (C) TCR sequencing of FFPE tumor samples and blood. The most prevalent TCR clonotypes (top 5 for each sample and blood) are shown (Table S5A). The blood sample was collected from the patient 550 days after secondary debulking (Figure S6A). Inset shows detection of the most frequent TCR rearrangement (CASSNDEYRGPTYEQYF) and its abundance comparison between samples (two-sided binomial tests with Benjamini-Hochberg multiple test correction, ∗∗∗ q < 0.001).
Figure S4
Neoepitope Distributions and HLA-I Neoepitope Depletion Analysis, Related to Figure 4 and Table S4 (A) Number of unique and overlapping expressed missense mutations, HLA-I and II neoepitopes between samples (Table S4D). (B) Correlations between expressed missense mutations and predicted HLA-I neoepitopes using NetMHC applied to TCGA ovarian samples (n = 150) and the primary and metastatic tumors (Tables S4E–S4G). KDE = kernel density estimate. (C) Top: Estimated neoepitope deviation from expected in the five tumor samples compared to TCGA ovarian cancer samples (n = 150). The expected number of neoepitopes was calculated by taking into account the expected number of missense mutations and the number of silent mutations according to Rooney et al., 2015 (see STAR Methods). Bottom: Neoepitope depletion analysis of 150 random unique permutations of the patient’s tumors (primary, spleen, RUQ, liver, and vaginal cuff) and their mutations. Each sample was compared against its own 150 unique permutations to control for the number of mutations (Tables S4H and S4I). Two-sided empirical p values were calculated from each distribution.
Figure S5
Predicted Immunogenicity of HLA Class I Neoepitopes, Related to Figures 1 and 4 and Tables S1 and S4 (A) Predicted immunogenic properties of trunk (clonal) and private HLA-I neoepitopes. Positive immunogenicity scores have biochemical properties associated with higher immunogenicity that outweigh properties associated with lower immunogenicity, and vice versa for negative scores (Calis et al., 2013). Horizontal lines within violin plots show the median and interquartile range of the data distribution. (B and C) Comparison between clonal and sub-clonal (including shared between two or more samples but not all) predicted immunogenicity of predicted binders and non-binders (two-sided Mann-Whitney rank test). Horizontal lines within violin plots show the median and interquartile range of the data distribution. (D–F) Probability of an HLA-I neoepitope having immunogenic properties considering its clonality and HLA-I binding affinity using the neoepitope data in (A), (B), and (C), respectively. Clonal neoepitopes have a lower probability of having immunogenic properties than sub-clonal predicted binders (chi-square test, p = 0.02). For non-binders (NetMHC score > HLA-I specific cutoff), clonal neoepitopes have a higher probability of having immunogenic properties (chi-square test, p = 0.003), as well as peptides with higher HLA-I affinities (chi-square test, p = 0.0001), although the absolute differences are minor. No significant interaction between clonality and predicted HLA-I binding affinity was detected for either binders or non-binders. GLM = generalized linear model.
Figure S6
PBMCs Sample Timeline and T Cell-Neoepitope Recognition Assay, Related to Figures 4 and 5 and Table S5 (A) Blood samples obtained from the patient 550 and 978 days after resection were used for TCR sequencing and T cell – neoepitope recognition assays respectively. (B) Experimental setup and flow cytometry gating strategy for the T cell –neoepitope recognition assays (intracellular cytokine staining assay) with surface staining of CD3, CD4, CD8, CD45, and intracellular staining of IL-4, IFN-γ, TNF-α. PBMC = peripheral blood mononuclear cells.
Figure 5
Predicted Neoepitopes with Higher Mutant than Wild-Type HLA-I Binding Affinity Elicit a T Cell Response (A) Representative scatterplots of TNF-α and IFN-γ intracellular cytokine staining of CD8+ T cells after 21 days of culture with CEF peptides or DMSO as positive and negative controls or the predicted mutant peptides (Figure S6B). CEF = Cytomegalovirus, Epstein-Barr virus, Influenza virus. (B) Percentage of CD8+ T cells with double-positive intracellular staining (TNF-α and IFN-γ) after incubation with each of the 43 predicted HLA-I neoepitopes, and HLA-I predicted binding affinity wild-type to mutant ratio (Table S5B). Mutation in gene FLG2 E1608K (P12) was found to be clonal after manual inspection in IGV (Table S1A).
Figure S7
Overall Associations between Tumor Fates and Tumor-Immune Microenvironmental Features, Related to Figures 1, 2, 3, and 4 and Tables S1, S2, S3, and S4 Cellular and molecular associations with change in tumor growth. Change in tumor growth (y axis) was calculated by dividing the tumor volume at CT scan 11 by the tumor volume at CT scan 10 (Figure 1B). Fitted curves are 2nd order polynomial regression lines plotted for trend visualization rather than prediction purposes. Capase 1 and 4 are considered inflammatory caspases involved in a type of apoptosis related to immune response called pyroptosis. The enrichment score x axis and the _q_-values come from the ssGSEA analysis. HLA-I genes include HLA-A, B, C, E, and F. HLA-II genes include HLA-DPA1, DMA, DRA, DQA1, DMB, DPB1, DQB2, DRB5, DRB1, DQB1, and DOA.
Comment in
- Cancer Evolution Constrained by the Immune Microenvironment.
McGranahan N, Swanton C. McGranahan N, et al. Cell. 2017 Aug 24;170(5):825-827. doi: 10.1016/j.cell.2017.08.012. Cell. 2017. PMID: 28841415
Similar articles
- Unraveling tumor-immune heterogeneity in advanced ovarian cancer uncovers immunogenic effect of chemotherapy.
Jiménez-Sánchez A, Cybulska P, Mager KL, Koplev S, Cast O, Couturier DL, Memon D, Selenica P, Nikolovski I, Mazaheri Y, Bykov Y, Geyer FC, Macintyre G, Gavarró LM, Drews RM, Gill MB, Papanastasiou AD, Sosa RE, Soslow RA, Walther T, Shen R, Chi DS, Park KJ, Hollmann T, Reis-Filho JS, Markowetz F, Beltrao P, Vargas HA, Zamarin D, Brenton JD, Snyder A, Weigelt B, Sala E, Miller ML. Jiménez-Sánchez A, et al. Nat Genet. 2020 Jun;52(6):582-593. doi: 10.1038/s41588-020-0630-5. Epub 2020 Jun 1. Nat Genet. 2020. PMID: 32483290 Free PMC article. - GATA3 as a regulator for naughty cancer-associated fibroblasts in the microenvironment of high-grade serous ovarian cancer.
El-Arabey AA, Abdalla M. El-Arabey AA, et al. Hum Cell. 2021 Nov;34(6):1934-1936. doi: 10.1007/s13577-021-00598-w. Epub 2021 Aug 25. Hum Cell. 2021. PMID: 34432264 No abstract available. - Association and prognostic significance of BRCA1/2-mutation status with neoantigen load, number of tumor-infiltrating lymphocytes and expression of PD-1/PD-L1 in high grade serous ovarian cancer.
Strickland KC, Howitt BE, Shukla SA, Rodig S, Ritterhouse LL, Liu JF, Garber JE, Chowdhury D, Wu CJ, D'Andrea AD, Matulonis UA, Konstantinopoulos PA. Strickland KC, et al. Oncotarget. 2016 Mar 22;7(12):13587-98. doi: 10.18632/oncotarget.7277. Oncotarget. 2016. PMID: 26871470 Free PMC article. - Low-grade serous ovarian carcinoma: an evolution toward targeted therapy.
Voutsadakis IA. Voutsadakis IA. Int J Gynecol Cancer. 2020 Oct;30(10):1619-1626. doi: 10.1136/ijgc-2019-000832. Epub 2019 Nov 27. Int J Gynecol Cancer. 2020. PMID: 31780569 Review. - Multi-omics in high-grade serous ovarian cancer: Biomarkers from genome to the immunome.
Clifford C, Vitkin N, Nersesian S, Reid-Schachter G, Francis JA, Koti M. Clifford C, et al. Am J Reprod Immunol. 2018 Aug;80(2):e12975. doi: 10.1111/aji.12975. Epub 2018 May 21. Am J Reprod Immunol. 2018. PMID: 29781548 Review.
Cited by
- Immunotherapy for ovarian cancer: towards a tailored immunophenotype-based approach.
Ghisoni E, Morotti M, Sarivalasis A, Grimm AJ, Kandalaft L, Laniti DD, Coukos G. Ghisoni E, et al. Nat Rev Clin Oncol. 2024 Sep 4. doi: 10.1038/s41571-024-00937-4. Online ahead of print. Nat Rev Clin Oncol. 2024. PMID: 39232212 Review. - Tumor battlefield within inflamed, excluded or desert immune phenotypes: the mechanisms and strategies.
Zheng S, Wang W, Shen L, Yao Y, Xia W, Ni C. Zheng S, et al. Exp Hematol Oncol. 2024 Aug 6;13(1):80. doi: 10.1186/s40164-024-00543-1. Exp Hematol Oncol. 2024. PMID: 39107856 Free PMC article. Review. - Total-Body Dynamic Imaging and Kinetic Modeling of [18F]F-AraG in Healthy Individuals and a Non-Small Cell Lung Cancer Patient Undergoing Anti-PD-1 Immunotherapy.
Omidvari N, Levi J, Abdelhafez YG, Wang Y, Nardo L, Daly ME, Wang G, Cherry SR. Omidvari N, et al. J Nucl Med. 2024 Sep 3;65(9):1481-1488. doi: 10.2967/jnumed.123.267003. J Nucl Med. 2024. PMID: 39089813 Free PMC article. - Radiology and multi-scale data integration for precision oncology.
Paverd H, Zormpas-Petridis K, Clayton H, Burge S, Crispin-Ortuzar M. Paverd H, et al. NPJ Precis Oncol. 2024 Jul 26;8(1):158. doi: 10.1038/s41698-024-00656-0. NPJ Precis Oncol. 2024. PMID: 39060351 Free PMC article. Review. - Poor efficacy of immune checkpoint inhibitor treatment in advanced thymic carcinoma patients with liver metastases.
Hao Y, Xu M, Zeng X, Wang Y, Wang W, Lin G, Li B, Huang J, Xu C, Zhang Y, Song Z. Hao Y, et al. Ther Adv Med Oncol. 2024 May 28;16:17588359241253127. doi: 10.1177/17588359241253127. eCollection 2024. Ther Adv Med Oncol. 2024. PMID: 38812990 Free PMC article.
References
- Bauer S., Groh V., Wu J., Steinle A., Phillips J.H., Lanier L.L., Spies T. Activation of NK cells and T cells by NKG2D, a receptor for stress-inducible MICA. Science. 1999;285:727–729. - PubMed
Publication types
MeSH terms
Substances
Grants and funding
- P30 CA008748/CA/NCI NIH HHS/United States
- 15601/CRUK_/Cancer Research UK/United Kingdom
- 20879/CRUK_/Cancer Research UK/United Kingdom
- 15871/CRUK_/Cancer Research UK/United Kingdom
- 15973/CRUK_/Cancer Research UK/United Kingdom
- 21141/CRUK_/Cancer Research UK/United Kingdom
- 22905/CRUK_/Cancer Research UK/United Kingdom
LinkOut - more resources
Full Text Sources
Other Literature Sources
Medical
Research Materials
Miscellaneous