Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Cancer stemness, intratumoral heterogeneity, and immune response across cancers

Alex Miranda, Phineas T. Hamilton, Allen W. Zhang, Swetansu Pattnaik, Etienne Becht, View ORCID ProfileArtur Mezheyeuski, Jarle Bruun, Patrick Micke, Aurélien de Reynies, and Brad H. Nelson
  1. aDeeley Research Centre, BC Cancer, Victoria, BC V8R 6V5, Canada;
  2. bDepartment of Molecular Oncology, BC Cancer, Vancouver, BC V5Z 4E6, Canada;
  3. cCentre for Molecular Medicine and Therapeutics, BC Children’s Hospital Research Institute, Vancouver, BC V5Z 4H4, Canada;
  4. dGraduate Bioinformatics Training Program, University of British Columbia, Vancouver, BC V6T 1Z3, Canada;
  5. eThe Kinghorn Cancer Centre and Cancer Division, Garvan Institute of Medical Research, Darlinghurst, NSW 2010, Australia;
  6. fSingapore Immunology Network, Agency for Science, Technology and Research, 138648 Singapore;
  7. gDepartment of Immunology, Genetics, and Pathology, Uppsala University, 751 85 Uppsala, Sweden;
  8. hDepartment of Molecular Oncology, Institute for Cancer Research, The Norwegian Radium Hospital, Oslo University Hospital, 0379 Oslo, Norway;
  9. iProgramme Cartes d’Identité des Tumeurs, Ligue Nationale Contre le Cancer, 75013 Paris, France;
  10. jDepartment of Biochemistry and Microbiology, University of Victoria, Victoria, BC V8P 3E6, Canada;
  11. kDepartment of Medical Genetics, University of British Columbia, Vancouver, BC V6T 1Z3, Canada

See allHide authors and affiliations

PNAS April 30, 2019 116 (18) 9020-9029; first published April 17, 2019; https://doi.org/10.1073/pnas.1818210116
Alex Miranda
aDeeley Research Centre, BC Cancer, Victoria, BC V8R 6V5, Canada;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Phineas T. Hamilton
aDeeley Research Centre, BC Cancer, Victoria, BC V8R 6V5, Canada;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Allen W. Zhang
bDepartment of Molecular Oncology, BC Cancer, Vancouver, BC V5Z 4E6, Canada;
cCentre for Molecular Medicine and Therapeutics, BC Children’s Hospital Research Institute, Vancouver, BC V5Z 4H4, Canada;
dGraduate Bioinformatics Training Program, University of British Columbia, Vancouver, BC V6T 1Z3, Canada;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Swetansu Pattnaik
eThe Kinghorn Cancer Centre and Cancer Division, Garvan Institute of Medical Research, Darlinghurst, NSW 2010, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Etienne Becht
fSingapore Immunology Network, Agency for Science, Technology and Research, 138648 Singapore;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Artur Mezheyeuski
gDepartment of Immunology, Genetics, and Pathology, Uppsala University, 751 85 Uppsala, Sweden;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Artur Mezheyeuski
Jarle Bruun
hDepartment of Molecular Oncology, Institute for Cancer Research, The Norwegian Radium Hospital, Oslo University Hospital, 0379 Oslo, Norway;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Patrick Micke
gDepartment of Immunology, Genetics, and Pathology, Uppsala University, 751 85 Uppsala, Sweden;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Aurélien de Reynies
iProgramme Cartes d’Identité des Tumeurs, Ligue Nationale Contre le Cancer, 75013 Paris, France;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brad H. Nelson
aDeeley Research Centre, BC Cancer, Victoria, BC V8R 6V5, Canada;
jDepartment of Biochemistry and Microbiology, University of Victoria, Victoria, BC V8P 3E6, Canada;
kDepartment of Medical Genetics, University of British Columbia, Vancouver, BC V6T 1Z3, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: bnelson@bccrc.ca
  1. Edited by Douglas Hanahan, Swiss Federal Institute of Technology Lausanne, Lausanne, Switzerland, and approved March 22, 2019 (received for review October 22, 2018)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

The exclusion of immune cells from the tumor microenvironment has been associated with poor prognosis in the majority of cancers. We report that when considering 21 solid cancer types, immune cell exclusion is widely associated with the presence of a stem cell-like phenotype in tumors (“stemness”). Stemness positively correlates with higher intratumoral heterogeneity, possibly by protecting antigenic clones from elimination by the immune system. The activation of a stemness program appears to limit antitumor immune responses via tumor cell-intrinsic silencing of endogenous retrovirus expression, repression of type I interferon signaling, and up-regulation of immunosuppressive checkpoints. Our work suggests that targeting the stemness phenotype in cancer will promote T cell infiltration and render tumors more responsive to immune control.

Abstract

Regulatory programs that control the function of stem cells are active in cancer and confer properties that promote progression and therapy resistance. However, the impact of a stem cell-like tumor phenotype (“stemness”) on the immunological properties of cancer has not been systematically explored. Using gene-expression–based metrics, we evaluated the association of stemness with immune cell infiltration and genomic, transcriptomic, and clinical parameters across 21 solid cancers. We found pervasive negative associations between cancer stemness and anticancer immunity. This occurred despite high stemness cancers exhibiting increased mutation load, cancer-testis antigen expression, and intratumoral heterogeneity. Stemness was also strongly associated with cell-intrinsic suppression of endogenous retroviruses and type I IFN signaling, and increased expression of multiple therapeutically accessible immunosuppressive pathways. Thus, stemness is not only a fundamental process in cancer progression but may provide a mechanistic link between antigenicity, intratumoral heterogeneity, and immune suppression across cancers.

  • cancer stemness
  • antitumor immunity
  • intratumoral heterogeneity

Tumor infiltration by T cells has been associated with improved clinical outcomes in a broad range of tumor types. Despite this, a large proportion of solid cancers appears nonpermissive to lymphocyte infiltration or nonimmunogenic (immunologically “cold”), and thus protected from cytolytic attack by lymphocytes, such as CD8+ T cells (1). The advent of immunotherapies, such as immune checkpoint inhibitors that rely on preexisting antitumor immune responses, has made an improved understanding of the mechanisms underlying the cold tumor phenotype essential.

Mounting evidence suggests that tumor cells can exhibit stem cell-like properties, ranging from characteristic gene-expression profiles to experimentally validated long-term self-renewal and repopulation capacities. The cancer stem cell (CSC) hypothesis posits that a subpopulation of tumor cells is capable of self-renewal and is responsible for the long-term maintenance of tumors (2). This hypothesis provides compelling explanations for clinical observations, such as therapeutic resistance, tumor dormancy, and metastasis (3). CSCs have been identified in a variety of human tumors, as assayed by their ability to initiate tumor growth in immunocompromised mice (4, 5). However, considerable controversy remains as to how best to define CSCs and the extent to which different tumor types exhibit a hierarchical organization. These controversies notwithstanding, there is increasing evidence that stem cell-associated molecular features, often referred to as “stemness,” are biologically important in cancer (6). It is unclear whether the stemness phenotype reflects the presence of bona fide CSCs in tumors or simply the coopting of stem cell-associated programs by non-CSC tumor cells (or both). Whatever the underlying mechanism may be, stemness has emerged as an important phenomenon due to its strong association with poor outcomes in a wide variety of cancers (7, 8). Moreover, stemness appears to be a convergent phenotype in cancer evolution (9, 10), suggesting it is a fundamentally important property of malignancy.

The evolution of transformed cells in the tumor microenvironment is shaped by diverse selective pressures, including the host immune response. Experimental work has shown that embryonic, mesenchymal, and induced pluripotent stem cells possess immune modulatory properties, while resistance to immune-mediated destruction has also recently been shown to be an intrinsic property of quiescent adult tissue stem cells (11) and CSCs (12). Similarly, immune selection has been shown to drive tumor evolution toward a stemness phenotype that inhibits cytotoxic T cell responses (13). Moreover, a recent analysis of The Cancer Genome Atlas (TCGA) revealed negative associations between stemness and some metrics of tumor leukocyte infiltration (14). Finally, CSCs have been proposed as a driver of intratumoral heterogeneity (6, 9). Consistent with this, we (15) and others (16) have reported negative associations between immune cell infiltration and intratumoral heterogeneity.

Motivated by these observations, we hypothesized that the stemness phenotype of cancer cells may confer immunosuppressive properties on tumors, resulting in immunologically cold microenvironments that both foster and maintain intratumoral heterogeneity. To address this, we performed an integrated analysis of stemness, immune response, and intratumoral heterogeneity across cancers. We recover pervasive negative associations between antitumor immunity and stemness, and strong positive associations between stemness and intratumoral heterogeneity. We further find that cancer cell lines with high stemness have cell-intrinsic immunosuppressive features, suggesting that immunologically cold microenvironments can arise due to the presence of high-stemness cancer cells. We propose that cancer stemness provides a link between tumor antigenicity, intratumoral heterogeneity, immune suppression, and the resulting evolutionary trajectories in human cancer.

Results

Derivation and Comparison of Stemness Signatures.

Recent studies have provided evidence that cancer stemness can be represented by core gene-expression programs across diverse cancer types (8, 10, 14, 17, 18). Building on this prior work, we inferred tumor stemness from cancer transcriptomes using single-sample gene set enrichment analysis (ssGSEA) with a modified version of a gene set developed by Palmer et al. (17) to measure the level of plasticity and differentiation of mesenchymal stem cells, pluripotent stem cells, terminally differentiated tissues, and human tumors across >3,200 microarray samples. Intriguingly, the authors of this gene set identified a cluster of “immune” genes that negatively loaded the principal components they used to infer stemness; however, they did not further explore this relationship. To adopt this gene set for use in ssGSEA and avoid biasing our analysis toward recovering negative associations between stemness and immunity, we omitted this immune gene cluster from our signature. We also omitted cell proliferation markers to avoid recovering a signature of proliferation rather than stemness (19).

We validated the performance of the resulting 109 gene signatures (Dataset S1A) on diverse datasets, finding that it recapitulated the expected degree of stemness in both malignant and nonmalignant cell populations (SI Appendix, Fig. S1). We further validated this signature using the stem cell-based validation dataset used by Malta et al. (14) in their recent pan-cancer stemness analysis (GSE30652), which revealed similarly high classification accuracy for our signature compared with theirs [multiclass area under the curve (AUC) 0.92 vs. 0.91, respectively] (SI Appendix, Fig. S2 A and B). We then analyzed RNA sequencing data from 8,290 primary cancers representing 21 solid cancer types (from TCGA) and found our signature showed good concordance with that of Malta et al. and two other recently published signatures (Spearman’s ρ = 0.43, 0.74, and 0.66, respectively) (SI Appendix, Fig. S2 C–F).

Stemness Varies Across Cancers and Predicts Patient Survival.

Using our signature, we found that stemness varied strongly across TCGA samples, with cancer type explaining 54% of the variation (ANOVA; adjusted R2) (Fig. 1A). Consistent with prior reports of stemness being a negative prognostic factor (7, 8), we found a strong negative relationship between median stemness and median overall survival across cancer types (Fig. 1B) (ρ = −0.60; P = 0.004; n = 21 cancers). Within cancer types, Cox regression likewise showed stemness to be significantly negatively prognostic for overall survival in the majority of cancers (Cox proportional hazards; P < 0.05), and significantly positively prognostic for none (Fig. 1C), underscoring the relevance of this signature both within and across cancers. We also noted a significant decrease in the magnitude of the hazard associated with stemness within cancers as median stemness increased (ρ = −0.59; P < 0.01) (Fig. 1C), pointing to a potential threshold effect with a saturating hazard in cancers with higher average stemness. For reference, we compared the prognostic association of our ssGSEA-based stemness with the mRNAsi signature of Malta et al. (14) derived using one-class logistic regression (OCLR), for which positive associations in some cancers were reported. Using pan-cancer Cox regressions stratified by cancer, we found ssGSEA-based stemness to be substantially more predictive of survival in this modeling framework [log hazard ratio = 0.23 ± 0.03 (coefficient ± SE); P < 10−15 vs. not significant for mRNA stemness index (mRNAsi)], demonstrating this signature uncovers the negative outcomes associated with high-stemness cancers expected from previous reports (7, 8). This relationship held when controlling for tumor purity, demonstrating that a relationship between tumor purity and patient survival is not a confounding factor (and see below).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Stemness and survival across cancers. (A) Stemness score varies widely across 21 solid cancers from TCGA. Each point represents an individual case, and cancer types are ordered by median stemness score (z-scored ssGSEA). (B) Median survival decreases with increasing median stemness (P = 0.004). Gray points represent cancers in which median overall survival times were not evaluable. (C) Stemness associates with poor outcome within cancers. Log hazard ratio (±95% CI) for the association of stemness with overall survival is shown. Hazard decreases with increasing average stemness of cancers (P = 0.008). Cox models control for patient age and tumor purity. Cancer acronyms are used as defined by TCGA (https://portal.gdc.cancer.gov).

Stemness Negatively Associates with Immune Cell Infiltration Across Solid Cancers.

To evaluate the relationship between stemness and antitumor immunity, we generated signatures of predicted immune cell infiltration for each patient sample using xCell, an ssGSEA-based tool that infers cellular content in the tumor microenvironment (20). CD8+ T cells, which have a well-established association with favorable prognosis in a majority of solid cancers (21), showed a clear negative association with stemness for most cancers (Fig. 2A). We also considered other cell types important for antitumor immunity, including NK cells and B cells, and again observed recurrent negative associations with stemness (Fig. 2A). Additional cell types, such as CD4+ T cells, Tregs, and neutrophils, showed more variable associations with stemness, indicating this relationship does not apply to all infiltrating immune cell populations in all cancers (Fig. 2A).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Stemness negatively associates with immune cell signatures. (A) Circos plot showing the association between stemness score and the presence of 8 inferred immune cell types across 21 cancer types (colored bars in outer ring). The color and height of the inner bars represent the Spearman correlation ρ values for each cell type and cancer type. (B) Volcano plot showing the association between stemness score and immune signature (sum of z-scored signatures of CD8+ T cells, NK cells, and B cells) for each cancer. The dashed line indicates Padj = 0.05. (C) Association between stemness score and immune signature in the different molecular subtypes of endometrial (UCEC) and breast (BRCA) cancer, and within primary and metastatic melanoma (SKCM) samples (Padj < 10−7). Each point represents one case. Colors indicate the different molecular subtypes of UCEC and BRCA, or sample types for SKCM. CN, copy number; MSI, microsatellite instable; POLE, polymerase epsilon.

To generate a robust single score for antitumor immunity (hereafter referred to as “immune signature”), we aggregated xCell predictions for CD8+ T cell, NK cell, and B cell infiltration, reasoning that these cells represent important anticancer effector cells across diverse cancers (22, 23). Indeed, this immune signature was significantly associated with increased survival in the majority of cancers (SI Appendix, Fig. S3A; overall survival curves for all patients stratified by median stemness and immune signature are shown in SI Appendix, Fig. S3B). A notable exception was kidney cancer (i.e., KIRP), where a negative association was observed, in accord with prior reports for this malignancy (24). Consistent with our analyses using single immune cell-type scores (Fig. 2A), the immune signature showed a negative association with stemness within nearly all cancers (Fig. 2B). This negative association was also recovered using other published stemness gene sets (18, 19, 25) (Dataset S1B), most clearly with stemness scores reflecting NANOG, SOX2, and MYC signaling, and to a lesser extent with those reflecting embryonic stem cell programs (SI Appendix, Fig. S4). Furthermore, when we used CIBERSORT infiltration scores (26) in place of our immune signature, we found that increased stemness was associated with strong polarization of infiltrating leukocyte populations toward a macrophage-dominated, CD8+ T cell-depleted composition for most cancers (linear model controlling for cancer type; P < 10−15 for both cell fractions).

Whereas the preceding analyses were performed within cancer types, we also evaluated the relationship between stemness and immune signature across cancer types. Unexpectedly, we found no significant association between median immune signature and median stemness score across cancers (P = 0.82). This suggests that factors other than stemness control the differences in immune cell infiltration across cancer types, while the association with stemness applies within individual cancer types. It also demonstrates that our stemness metric is not recovering negative associations with immunity simply due to the lower tumor purity that is inextricably associated with the presence of infiltrating immune cells.

We next examined whether the negative association between immune signature and stemness was influenced by tumor subtype or stage. Here, we focused on breast and endometrial cancers, which have well-characterized subtypes with strong prognostic associations, and melanoma, for which both primary and metastatic samples are available within TCGA. In breast cancer, stemness varied markedly across subtypes (ANOVA; adjusted R2 = 0.26), with the basal subtype having the highest stemness, as expected (19), and the luminal-A subtype the lowest (Fig. 2C). In endometrial cancer, the highest stemness score was observed in high-copy number (CN) alteration (CN-high) and polymerase-epsilon mutant tumors, and the lowest stemness score was seen in low-CN alteration tumors (CN-low) [Tukey honest significant difference (HSD); P = 0.003; CN-high vs. CN-low tumors]. Finally, we observed a substantially higher stemness score in metastatic compared with primary melanoma lesions (Fig. 2C). In all these cancers, we observed recurrent negative associations between stemness and immune signature which remained significant when controlling for cancer subtype and tumor purity (see below; linear models; P < 10−7).

To investigate in an unbiased manner whether processes apart from antitumor immunity negatively correlate with stemness, we conducted differential expression tests to identify gene-expression patterns associated with the lowest versus highest stemness quintiles (<20th versus >80th percentiles) for each analyzed cancer type. Even with this unbiased approach, nearly all of the pathways recurrently enriched in low-stemness samples within a cancer were immune-related (SI Appendix, Fig. S5). Recognizing that the presence of nonmalignant cells can confound expression analyses of bulk-sequenced tumor samples by diluting tumor-specific expression signatures, we performed additional analyses to control for such effects. First, using recently published estimates of purity across TCGA (27), we found that the association between tumor purity and stemness was negligible and failed to reach significance in a pan-cancer linear model controlling for cancer type (P = 0.23). Second, we refit the differential expression models described above to control for tumor purity and repeated the pathway enrichment analyses, still finding that immune-associated pathways were enriched in low-stemness tumors (SI Appendix, Fig. S6).

To evaluate the relative contributions of malignant versus nonmalignant cells (e.g., stromal cells) to the stemness score, we applied our stemness metric to a single-cell RNA sequencing (scRNA-seq) dataset of lung cancer (28), which comprises a comprehensive inventory of different cell types from the tumor microenvironment (>52,000 cells from five patients). The average stemness score was much higher in cancer cells than any other cell type, including fibroblasts, myeloid cells, or lymphocytes (SI Appendix, Fig. S7) (Tukey HSD; P < 10−15), strongly supporting the notion that the stemness signature largely emanates from cancer cells rather than stromal or other cell types in the tumor microenvironment.

Stemness Associates with Immunologically Cold Cancers Measured via Immunohistochemistry.

To confirm the negative association between stemness and lymphocyte infiltration, we turned to three patient cohorts with matched immunohistochemistry (IHC)-based T cell infiltration scores and gene-expression data suitable for computing a stemness score. Using a cohort of 33 colorectal cancer patients (29), we found a strong negative association between stemness and total infiltrating CD3+ T cells (ρ = −0.63; P < 0.001) (Fig. 3A). Furthermore, in this cohort the xCell-based immune signature was strongly correlated with infiltrating CD3+ cells (ρ = 0.69; P < 0.001), supporting its fidelity for measuring immune-cell infiltration. With this validation in hand, we compared the xCell-based immune signature and stemness for the total patient cohort with available microarray data (n = 585), and again observed a clear negative correlation (ρ = −0.22; P < 10−7).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Relationship between stemness and immune cell infiltrates in different cancer cohorts scored via IHC. (A) Stemness score negatively associates with tumor-infiltrating CD3+ T cells in colorectal cancer (P < 0.01; n = 33; data from ref. 29). Each point represents one patient sample. (B) Stemness negatively associates with tumor-infiltrating T cells (in this case the sum of CD4+ and CD8+ cells) in lung cancer (P = 0.07; n = 35; data from ref. 30). Colors denote adeno versus squamous cell lung cancers. (C) Stemness is significantly associated with tumor-infiltrating CD3+ T cells in a multisite dataset of high-grade ovarian cancer (P = 0.028; n = 44 samples from 12 patients; data from ref. 15). Colors represent individual patients. (D) By hotspot analysis, the fraction of tissue area occupied by colocalizing tumor and immune cells (FCI) is negatively associated with stemness in a multisite dataset from high grade ovarian cancer (P < 0.01; 44 samples from 12 patients; data from ref. 15). Colors represent individual patients.

We next evaluated this relationship in a cohort of 35 lung cancer patients with matched RNA-sequencing and IHC-based quantitation of immune cell infiltrates (30). For consistency with the above analysis, we calculated the infiltration of T cells by summing previously calculated CD4+ and CD8+ cell fractions. Although we observed a negative association between the stemness score and the percent of infiltrating T cells, this did not reach statistical significance (ρ = −0.32; P = 0.07) (Fig. 3B). Nonetheless, the additional RNA-seq data in this cohort revealed a clear negative association between the stemness and immune signature at the transcriptional level (n = 199 tumor samples; ρ = −0.19; P = 0.007).

Finally, we evaluated a small cohort of high-grade serous ovarian cancer (HGSC) cases (15) for which matched IHC-based T cell counts and microarray-based gene-expression data were available for multiple tumor sites within each patient (n = 44 samples from 12 patients). Consistent with the above findings, we found a negative association between stemness and total CD3+ T cells (negative binomial mixed effects model; P = 0.0028) (Fig. 3C). In addition, we previously subjected samples from this cohort to Getis-Ord Gi* “hotspot” analysis to quantify immune cell engagement with tumor cells (15). Intriguingly, all hotspot metrics showed a clear negative association with stemness (mixed-effects models; P < 0.01; 44 samples from 12 patients) (Fig. 3D, FCI shown) (31), suggesting that stemness negatively influences lymphocyte engagement with tumor cells.

Stemness Associates with Intratumoral Heterogeneity.

Stemness has been proposed to foster tumor clone diversity, with the replicative potential of CSCs enabling greater tumor heterogeneity (6, 9). Our hypothesis suggests that stemness could additionally promote intratumoral heterogeneity by inhibiting immune selection against new cancer clones. These predictions have not, to our knowledge, been systematically tested within or across cancers. Therefore, we compared stemness and intratumoral heterogeneity using data from two recent TCGA studies (26, 32). Using data from the first study (32), we found a dramatic positive correlation between median stemness and median number of clones in a cancer (n = 935 patients across 11 cancers; ρ = 0.75; P = 0.008) (Fig. 4A). Furthermore, we found a positive association between stemness and tumor clone count in a linear model controlling for cancer site and tumor purity, indicating that these associations are discernible both across and within cancers (P = 0.0002) (Fig. 4B). Using pan-cancer predictions of intratumoral heterogeneity from a second, larger TCGA-based study (26), we recovered a similarly strong association (ρ = 0.64; P = 0.002; n = 6,791 samples across 21 cancers) (Fig. 4C), which was again significant across all samples when controlling for cancer site and tumor purity in a linear model (P < 10−15) (Fig. 4D).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Stemness associates with intratumoral heterogeneity within and across cancers. (A) Median stemness and median clonality [inferred by Andor et al. (32)] are strongly correlated across cancers (n = 11; P = 0.008). (B) Stemness score and clonality [inferred by Andor et al. (32)] are correlated across patients while controlling for cancer type (n = 935; P = 0.0002). Colored points represent different tumor sites. (C) Median stemness score and median intratumoral heterogeneity score [inferred by Thorsson et al. (26)] are strongly correlated across cancers (n = 20; P = 0.002). (D) Stemness score and intratumoral heterogeneity [inferred by Thorsson et al. (26)] are correlated across patients while controlling for cancer type (n = 6,791; P < 10−15). Colored points represent different tumor sites. Spearman ρ values are shown.

Potential Mechanisms and Consequences of Stemness-Associated Immunosuppression.

We evaluated the association between stemness and several known genetic and environmental factors that affect antitumor immunity. Antitumor immunity involves T cell recognition of neo-antigens arising from somatic mutations (33). Therefore, we examined the association between stemness, immune signature, and nonsynonymous mutation load, analyzing TCGA samples with available mutation calls (n = 6,682). While median mutation load correlated with median immune signature across cancers (ρ = 0.49; P = 0.02, n = 21), there was generally little correlation within cancers, as has been reported in other TCGA-based analyses (34) (SI Appendix, Fig. S8A). Across cancers, median mutation load showed a positive association with median stemness (ρ = 0.52; P = 0.015, n = 21) (Fig. 5A). Similarly, within cancers, we generally found positive associations between mutation load and stemness (Fig. 5B). We recovered qualitatively similar but slightly stronger associations between stemness and neoantigen load computed with NetMHCpan in the above pan-cancer analysis (26) (SI Appendix, Fig. S8 B and C).

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Mutation load, CT antigen expression and ERV associations with stemness. (A) Median stemness and median mutation load are positively correlated across cancers (n = 21; P = 0.015). Mutation load is represented as log-transformed nonsynonymous mutations per base (log10 ns mutations per base pair). (B) Volcano plot reveals stemness score and mutation load correlate within some cancers (upper right quadrant). The x axis represents Spearman correlation ρ values, and the y axis represents −log10-adjusted P values (Padj). Dashed red line indicates the significance threshold, Padj value = 0.05. (C) Stemness score and CT antigen expression (ssGSEA of CT antigen gene set) positively correlate in most cancers. Bar plots show the Spearman ρ values for each cancer type, and asterisks denote Padj < 0.05. (D) Redundancy analysis triplot reveals stemness negatively associates with multivariate ERV expression (P < 0.001; 33 ERVs evaluated in 4,252 samples, analysis conditioned by cancer type).

Cancer-testis (CT) antigens are tumor antigens that normally are expressed in gametogenic tissue but become aberrantly expressed in a broad range of malignancies, often leading to an immune response (35). Using a set of 201 CT genes curated by the CTdatabase (Dataset S2) (36), we generated an ssGSEA CT antigen score and found strong positive associations with stemness within cancers (Fig. 5C). Accordingly, there was a generally negative association between CT antigen score and the immune signature, which was significant in 8 of 21 cancers (Padj < 0.05) (SI Appendix, Fig. S9). Thus, like neoantigens, CT antigens show a positive association with stemness and a negative association with immune signature.

Normal stem cells have been shown to suppress endogenous retrovirus (ERV) expression, presumably to prevent insertional mutagenesis in long-lived stem cell lineages (37). Conversely, ERV expression can be activated in cancer cells (38, 39), where it can potentially elicit antitumor immune responses by activating viral defense mechanisms and the type I IFN response (40) or by yielding immunogenic foreign epitopes (41). Despite these possibilities, associations between immunity and ERV expression were inconsistently observed in a recent pan-cancer analysis (42). To better understand this relationship, we investigated interactions between stemness, immune signature, and ERV expression. Because of the repetitive nature of ERVs, we used ERV-specific read-mappings (42) to evaluate ERV expression in 4,252 TCGA samples that overlapped with our stemness and immune signature analysis. Using redundancy analysis (a constrained extension of principal components analysis), we found that multivariate ERV expression was not clearly associated with the immune signature, consistent with prior reports (42); members of the ERVK family were an exception, showing moderate positive associations with immune signature (Fig. 5D). In contrast, ERV expression showed a pervasive negative association with stemness (P < 0.001) (Fig. 5D), consistent with the notion that suppression of ERV expression is a feature of the stem cell phenotype (43). Thus, both immune signature and ERV expression are negatively associated with stemness, but these appear to be largely orthogonal relationships.

Tumor Cell Intrinsic Mechanisms of Stemness-Mediated Immunosuppression.

To address whether the negative association between stemness and immune signature is attributable to cancer cell-intrinsic processes, we calculated stemness scores for 1,048 cancer cell lines using gene-expression data from the Cancer Cell Line Encyclopedia (CCLE) (44), as well as the cancer cell fraction from the scRNA-seq study of lung cancer, mentioned above (28). We first assessed associations between stemness scores and the expression of 11 mapped ERVs in CCLE transcriptomes and found that three of three ERVs with nonnegligible expression levels across cell lines were significantly negatively associated with stemness, two of which remained significant when controlling for tissue of origin (linear models; Padj < 0.05) (SI Appendix, Fig. S10). We next generated an ssGSEA score for type I IFN signaling (reactome IFN α/β pathway gene set) and generally observed negative correlations with stemness within cancers in the TCGA dataset (Fig. 6A). When evaluating this association in the CCLE dataset, we observed a clear negative relationship between the IFN signature and stemness (ρ = −0.22. P < 10−10) (SI Appendix, Fig. S11A), which remained significant when controlling for tissue of origin of the cell line (linear model; P < 0.001) and when omitting cell lines derived from hematopoietic lineages (linear model; P = 0.02).

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Cell-intrinsic stemness score associates with decreased type I IFN signaling and increased expression of CD276 and PVR. (A) Volcano plot reveals that stemness score and type I IFN signaling (reactome IFN α/β pathway ssGSEA) are negatively correlated in most cancers (upper left quadrant). The x axis represents Spearman correlation ρ values, and the y axis represents −log10 adjusted P values (Padj). Dashed red line indicates the significance threshold, Padj = 0.05. (B) Stemness score is negatively associated with type I IFN signaling (reactome IFN α/β pathway ssGSEA) in tumor cells from four of five lung cancer patients, based on scRNA-seq data from ref. 28. Each point represents a single tumor cell. (C and D) Heatmaps showing Spearman correlations for stemness and select immunosuppressive genes based on data from the CCLE (C) and TCGA (D). Spearman correlations were calculated within tissues represented in the CCLE by more than 10 cell lines. Genes are ranked according to the final column (“overall”), which represents the correlation across all samples, irrespective of cell line or tumor type. Red-blue intensities reflect the correlation ρ values. Asterisks denote Benjamini–Hochberg-corrected significant associations (Padj < 0.05).

Analysis of the aforementioned lung cancer scRNA-seq dataset (28) also revealed a striking negative association between cancer cell-intrinsic stemness and IFN α/β signaling in four of five patients (Fig. 6B). To further test this association in nonneoplastic lineages, we took advantage of the stem cell gene-expression dataset previously used to validate stemness signatures (GSE30652), which likewise showed a striking negative association between stemness and type I IFN signaling (ρ = −0.81; P < 10−15) (SI Appendix, Fig. S11B).

Finally, to examine other cell-intrinsic mechanisms of immunosuppression, we analyzed a curated list of immunosuppressive genes previously reported to be expressed in human cancer cells. Using both the CCLE and pan-cancer TCGA datasets, the expression of each of these genes was assessed in relation to our stemness signature (Fig. 6 C and D). This revealed positive associations between stemness and a number of immunosuppressive genes, including CD276 [B7-H3, shown to inhibit T cell activation and autoimmunity (45)], PVR [CD155, a member of the B7/CD28 superfamily, shown to exhibit potent inhibitory action in different subsets of immune cells (46)], and TGFB1 [a key player in the induction of immunological tolerance (47)]. Thus, the stemness phenotype is associated with expression of several gene products that could potentially serve as targets for immune modulation.

Discussion

Although cancer stemness, antitumor immunity, and intratumoral heterogeneity have all emerged as important features of cancer in recent years, their covariation across cancers has not been systematically investigated. Here we report that stemness is associated with suppressed immune response, higher intratumoral heterogeneity, and dramatically worse outcome for the majority of cancers. Although correlative analyses such as ours do not reveal causality, we propose that the stemness phenotype found in cancer cells, similar to that in normal stem cells, involves the expression of immunosuppressive factors that engender the formation of immune-privileged microenvironments in which tumor clone diversification can occur. The resulting heterogeneity may provide a substrate for the selection of treatment-resistant clones, resulting in inferior clinical outcomes. Thus, our findings implicate stemness as a shared therapeutic target to achieve the dual objectives of constraining tumor evolution and enhancing antitumor immunity.

A recent pan-cancer analysis reported inconsistent relationships between cancer stemness and immunity, recovering negative relationships between stemness and tumor-infiltrating lymphocytes for some cancers and positive relationships for others (14). While this work provided a valuable perspective on stemness across cancers, our results differ in that we recover much stronger and more pervasive negative relationships with immune infiltration and survival. In contrast to our ssGSEA approach, the OCLR approach used above was validated on a cohort that lacked any malignant samples. Moreover, when we reproduced their analysis, we found many components of their stemness score (i.e., positive OCLR model weights) were immunologically relevant genes; for example, among the 50 most positive gene weights were IDO1, LCK, KLRG2, PSMB9 (a component of the immunoproteasome), and multiple TNF-receptors. With the OCLR approach, higher expression of such immune genes contributes positively to the stemness score, precluding an unbiased assessment of the relationship between stemness and tumor immunity. As mentioned, it is also possible that an immune cell signature in the tumor microenvironment could negatively correlate with transcriptional signatures from cancer cells simply through dilution of cancer-specific transcripts by infiltrating immune cells. We took numerous precautions to ensure this was not driving our results. These included controlling for tumor purity in linear models throughout the analyses, showing a nonsignificant relationship between stemness and tumor purity, and identifying a negative relationship between stemness and cell-intrinsic IFN signaling in individual lung cancer cells (by scRNA-seq; Fig. 6B) and in cancer cell line transcriptomes (SI Appendix, Fig. S11A).

A contribution of cancer stemness to intratumoral heterogeneity has been postulated for some time (6, 9), but direct evidence has been lacking. We recovered a dramatic positive association between stemness and multiple metrics of intratumoral heterogeneity across cancers (Fig. 4), which is especially noteworthy given that these metrics were derived from different data types (i.e., mRNA vs. DNA). Given recent work from our group and others linking increased intratumoral heterogeneity with decreased immune cell infiltration (15, 16), one could speculate that stemness might contribute to intratumoral heterogeneity by both increasing the replicative capacities of individual tumor clones and by shielding antigenic clones from elimination by the immune system.

We found generally positive associations between stemness and mutation load within cancers (Fig. 5B), and clear evidence of this across cancers (Fig. 5A), consistent with studies demonstrating accumulation of mutations in normal adult stem cells (48). We also found strong positive associations between stemness and CT antigen expression (Fig. 5C), which is consistent with reports of CT antigen expression in mesenchymal (49), embryonic stem cells (50), and CSCs (51). It is also consistent with prior studies (42) and the present work finding nonsignificant or negative correlations between immune infiltration and CT antigen expression (Fig. 5D). Thus, the negative association between stemness and immune response is not readily attributable to low neoantigen or CT antigen load, strongly implicating the involvement of other mechanisms.

ERVs, which constitute ∼8% of the human genome, are known to be suppressed in pluripotent and embryonic stem cells (43) yet activated in human cancer (38, 39), leading us to ask which behavior would predominate in high stemness cancers. We found a strong negative association between stemness and ERV expression in both TCGA (Fig. 5D) and cancer cell line data (SI Appendix, Fig. S10), indicating that the stemness phenotype in human cancer retains this property of normal stem cells.

We identified a negative correlation between stemness and type I IFN signaling in the TCGA (Fig. 6A), CCLE (SI Appendix, Fig. S11A), and scRNA-seq datasets (Fig. 6B). Whether ERV suppression underlies the low intrinsic IFN signaling in high stemness cancer cells awaits experimental investigation. In support of our results, an attenuated innate immune response is a major characteristic of embryonic stem cells (52), and we find clear confirmation of this in nonneoplastic stem cell transcriptomes (SI Appendix, Fig. S11B). Collectively, these data suggest that activation of a stemness program in tumors could limit antitumor immune responses by silencing ERVs and repressing type I IFN signaling in a cell-intrinsic manner (53).

Using CCLE data, we found a clear association between stemness and the expression of several immunosuppressive genes, including CD276, PVR, and TGFB1 (Fig. 6C). These associations are especially intriguing given that CCLE-based gene expression profiles are independent of any ongoing influence of the immune system. CD276, a B7 family ligand, is now being clinically targeted due to expression on both cancer cells and tumor-infiltrating blood vessels (54). Intriguingly, CD276 is coexpressed with CD133, a marker that distinguishes cell populations enriched for CSCs in colorectal cancer (55). PVR is a key ligand in an emerging checkpoint pathway involving TIGIT, an inhibitory receptor expressed on T cells and other immune cells (46). Although no association between stem cells and PVR has been described so far in mammalian stem cells, expression of PVR can be activated by sonic hedgehog signaling (56), a pathway essential for self-renewal and cell fate determination in normal and CSCs (57). TGFB1 and other TGFB family members have well-documented roles in development (58) and CSC proliferation and maintenance (59, 60).

Recent work has demonstrated that tumor-intrinsic oncogenic signaling pathways have immune suppressive properties (61), but this too can be understood in the context of stemness. For example, molecular pathways involving WNT/B-catenin, MYC, PTEN, and LKB1 have been implicated in the inhibition of antitumor immunity (61), yet they also play important roles in stem cell maintenance (62⇓⇓–65). Thus, stemness may provide a unifying framework for understanding how various oncogenic signaling pathways engender an immunosuppressive tumor microenvironment.

Although it remains unclear if stemness metrics derived from bulk tumor samples represent rare populations of bona fide CSCs or a wholesale shift of the cancer cell population toward a higher stemness phenotype, our findings provide rationale for therapeutic targeting of the stemness phenotype itself. In particular, if stemness plays a causative role in the formation of cold tumor microenvironments, it may prove beneficial to target specific molecules or pathways that are inherent to the stemness phenotype, such as the aforementioned immunosuppressive molecules. High-stemness cancers might also be rendered more sensitive to immune control by administering drugs that induce cell differentiation to irreversibly disrupt the stemness phenotype (66). This might bring the additional benefit of constraining further tumor evolution, creating the conditions for durable clinical responses.

Materials and Methods

Data Acquisition and Processing.

All analyses were conducted with R software newer than version 3.4.2. We extracted clinical parameters and molecular subtypes for TCGA data from the pan-cancer curated clinical data of Liu et al. (67). Mutation data were downloaded from Firebrowse (www.firebrowse.org) as the number of nonsynonymous mutations per base. Tumor purity (27), mRNAsi (14), intratumoral heterogeneity (26, 32), and CIBERSORT and neoantigen scores (26) were accessed by extracting the relevant data for overlapping samples from the respective supplemental materials.

For correlations between stemness scores and IHC-based immune cell counts, immune cell infiltration data were provided by the authors of the respective studies (29). We obtained matching expression data (microarray or RNA-sequencing) from the Gene Expression Omnibus (GEO: GSE39582, GSE81089), or from the authors (15, 29, 30). Other expression data for validation was obtained from the GEO (GSE30652, GSE15192, GSE31257, GSE76009).

We accessed RNA sequencing as upper quartile-normalized fragments per kilobase of transcript per million mapped reads using the TCGAbiolinks R/Bioconductor package (68), for each cancer of interest, and expression data were merged across cancers. For genes with multiple annotated transcripts, we selected the transcript with the highest expression to represent the gene, then filtered the expression set to include only primary samples (except for melanoma, for which we included metastases), removed patients with duplicate samples, and removed any patients without a consensus purity score in Aran et al. (27) to enable purity corrections in analyses. For microarray datasets, we converted probe IDs to human gene symbols using biomaRt (69), and retained the probe with the highest expression for each gene, as above. We accessed the scRNA-seq data and tSNE embeddings in Lambrechts et al. (28) from EBI (E-MTAB-6149), while cell type annotations and anonymized patient codes were provided by the authors. Where appropriate (e.g., for linear modeling), expression data were log2(x + 1) -transformed.

We calculated stemness and other ssGSEA signatures using the GSVA package in R (gene sets in Dataset S1) (70) without normalization, and subsequently scaled values as z-scores within datasets of interest. For ssGSEA calculations on scRNA-seq, we first omitted genes below the median of average expression across samples; this left representation for 61 of 109 of the genes in the stemness signature. xCell enrichment scores were calculated in R, using the rawEnrichment analysis (20) function, which omits scaling scores to [0, 1] and correction for correlations among related cell types (20), as we sought to avoid introducing nonlinearities from these steps into analysis. To generate the immune signature, we summed z-scored signatures of cell types of interest (CD8+ T cells, NK cells, B cells). Because z-scoring of ssGSEA scores was done within each dataset, these scores should not be directly compared across datasets.

Quantification and Statistical Analysis.

We used nonparametric Spearman’s correlation to assess pairwise associations between variables of interest within cancers, or for median values across cancers, adjusting for multiple tests using the Benjamini–Hochberg method, where appropriate. For analyses across multiple cancer types or subtypes, we used linear models, controlling for site or subtype as fixed main effects, and inspecting model residuals to ensure model assumptions were reasonable. For analyses controlling for purity, purity (as consensus purity estimate from ref. 27) was included as a covariate (main effect) in linear models.

We conducted survival analyses using Cox proportional hazards models, calculating 95% CIs on log hazard ratios. We tested model assumptions using cox.zph (71). Where there were significant violations of model assumptions (P < 0.05), we inspected model Schoenfield residuals. We found that for the few instances in which model assumptions were violated, this was attributable to higher than predicted survival for some long-term survivors. For analyses across all patients and cancer types, we stratified Cox models by cancer type. For IHC data with multiple samples taken from the same patient, we modeled patient as a random effect in linear mixed-effects models implemented in lme4 in R (72) and assessed the significance of fixed effects of interest using likelihood-ratio tests on nested models.

Pathway enrichment analysis was conducted using ReactomePA (73) after testing for differential expression between quintiles using limma (74). For enrichment analysis, we selected the top 1,000 significantly down-regulated genes in high stemness cancers based on the moderated t statistic; in cancers where <1,000 genes were significantly down-regulated (Padj < 0.05), we selected all significantly down-regulated genes for downstream analysis. Significance of enrichment was evaluated at Padj < 0.05, and recurrently enriched pathways were defined as those that were significantly enriched in the greatest number of cancers. In limma analyses that included purity as a covariate, purity was log-transformed for consistency with the transformation of expression values.

To assess ERV expression, we first variance-filtered mapped ERVs to select those above the median interquartile range of expression using the genefilter R package (75). We conducted partial redundancy analysis using the vegan package (76), implementing the default distance metric (Euclidean). We conditioned the analysis by cancer type to control for cancer type-specific effects, and tested the significance of multivariate associations using permutation tests (n = 1,000).

Acknowledgments

We thank Dr. John Dick (University of Toronto) and Dr. Julian Lum (BC Cancer) for critical review of the manuscript, and Dr. Peter Watson (BC Cancer) for helpful discussions. We are grateful to The Cancer Genome Atlas and Cancer Cell Line Encyclopedia for access to data that enabled this study. Funding was provided by the BC Cancer Foundation (to B.H.N.), Canadian Cancer Society (to B.H.N.), Canada’s Networks and Centres of Excellence (to B.H.N.), Canadian Institutes of Health Research (to B.H.N.), Terry Fox Research Institute (to B.H.N.), Cancer Research Society (to B.H.N.), Vanier Canada Graduate Scholarship (to A.W.Z.), and Canadian Institutes of Health Research Postdoctoral fellowships (to P.T.H. and A. Miranda).

Footnotes

  • ↵1A. Miranda. and P.T.H. contributed equally to this work.

  • ↵2To whom correspondence should be addressed. Email: bnelson{at}bccrc.ca.
  • Author contributions: A. Miranda, P.T.H., and B.H.N. designed research; A. Miranda and P.T.H. performed research; E.B., A. Mezheyeuski, J.B., P.M., and A.d.R. contributed new reagents/analytic tools; A. Miranda, P.T.H., A.W.Z., S.P., and B.H.N. analyzed data; B.H.N. conducted study supervision; and A. Miranda, P.T.H., and B.H.N. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1818210116/-/DCSupplemental.

Published under the PNAS license.

References

  1. ↵
    1. Gajewski TF
    (2015) The next hurdle in cancer immunotherapy: Overcoming the non-T-cell-inflamed tumor microenvironment. Semin Oncol 42:663–671.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Batlle E,
    2. Clevers H
    (2017) Cancer stem cells revisited. Nat Med 23:1124–1134.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Nassar D,
    2. Blanpain C
    (2016) Cancer stem cells: Basic concepts and therapeutic implications. Annu Rev Pathol 11:47–76.
    OpenUrlCrossRef
  4. ↵
    1. O’Brien CA,
    2. Pollett A,
    3. Gallinger S,
    4. Dick JE
    (2007) A human colon cancer cell capable of initiating tumour growth in immunodeficient mice. Nature 445:106–110.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Al-Hajj M,
    2. Wicha MS,
    3. Benito-Hernandez A,
    4. Morrison SJ,
    5. Clarke MF
    (2003) Prospective identification of tumorigenic breast cancer cells. Proc Natl Acad Sci USA 100:3983–3988.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Kreso A,
    2. Dick JE
    (2014) Evolution of the cancer stem cell model. Cell Stem Cell 14:275–291.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Ng SW, et al.
    (2016) A 17-gene stemness score for rapid determination of risk in acute leukaemia. Nature 540:433–437.
    OpenUrlCrossRef
  8. ↵
    1. Smith BA, et al.
    (2018) A human adult stem cell signature marks aggressive variants across epithelial cancers. Cell Rep 24:3353–3366.e5.
    OpenUrl
  9. ↵
    1. Greaves M
    (2013) Cancer stem cells as ‘units of selection’ Evol Appl 6:102–108.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Chen H,
    2. He X
    (2016) The convergent cancer evolution toward a single cellular destination. Mol Biol Evol 33:4–12.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Agudo J, et al.
    (2018) Quiescent tissue stem cells evade immune surveillance. Immunity 48:271–285.e5.
    OpenUrlCrossRef
  12. ↵
    1. Bruttel VS,
    2. Wischhusen J
    (2014) Cancer stem cell immunology: Key to understanding tumorigenesis and tumor immune escape? Front Immunol 5:360.
    OpenUrlPubMed
  13. ↵
    1. Noh KH, et al.
    (2012) Cancer vaccination drives Nanog-dependent evolution of tumor cells toward an immune-resistant and stem-like phenotype. Cancer Res 72:1717–1727.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Malta TM, et al.
    (2018) Machine learning identifies stemness features associated with oncogenic dedifferentiation. Cell 173:338–354.e15.
    OpenUrlCrossRef
  15. ↵
    1. Zhang AW, et al.
    (2018) Interfaces of malignant and immunologic clonal dynamics in ovarian cancer. Cell 173:1755–1769.
    OpenUrl
  16. ↵
    1. Safonov A, et al.
    (2017) Immune gene expression is associated with genomic aberrations in breast cancer. Cancer Res 77:3317–3324.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Palmer NP,
    2. Schmid PR,
    3. Berger B,
    4. Kohane IS
    (2012) A gene expression profile of stem cell pluripotentiality and differentiation is conserved across diverse solid and hematopoietic cancers. Genome Biol 13:R71.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Shats I, et al.
    (2011) Using a stem cell-based signature to guide therapeutic selection in cancer. Cancer Res 71:1772–1780.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Ben-Porath I, et al.
    (2008) An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat Genet 40:499–507.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Aran D,
    2. Hu Z,
    3. Butte AJ
    (2017) xCell: Digitally portraying the tissue cellular heterogeneity landscape. Genome Biol 18:220.
    OpenUrl
  21. ↵
    1. Gooden MJ,
    2. de Bock GH,
    3. Leffers N,
    4. Daemen T,
    5. Nijman HW
    (2011) The prognostic influence of tumour-infiltrating lymphocytes in cancer: A systematic review with meta-analysis. Br J Cancer 105:93–103.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Melero I,
    2. Rouzaut A,
    3. Motz GT,
    4. Coukos G
    (2014) T-cell and NK-cell infiltration into solid tumors: A key limiting factor for efficacious cancer immunotherapy. Cancer Discov 4:522–526.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Sarvaria A,
    2. Madrigal JA,
    3. Saudemont A
    (2017) B cell regulation in cancer and anti-tumor immunity. Cell Mol Immunol 14:662–674.
    OpenUrlCrossRef
  24. ↵
    1. Jochems C,
    2. Schlom J
    (2011) Tumor-infiltrating immune cells and prognosis: The potential link between conventional cancer therapy and immunity. Exp Biol Med (Maywood) 236:567–579.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Bhattacharya B, et al.
    (2004) Gene expression in human embryonic stem cell lines: Unique molecular signature. Blood 103:2956–2964.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Thorsson V, et al., Cancer Genome Atlas Research Network
    (2018) The immune landscape of cancer. Immunity 48:812–830.e14.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Aran D,
    2. Sirota M,
    3. Butte AJ
    (2015) Systematic pan-cancer analysis of tumour purity. Nat Commun 6:8971.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Lambrechts D, et al.
    (2018) Phenotype molding of stromal cells in the lung tumor microenvironment. Nat Med 24:1277–1289.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Becht E, et al.
    (2016) Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol 17:218.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Mezheyeuski A, et al.
    (2018) Multispectral imaging for quantitative and compartment-specific immune infiltrates reveals distinct immune profiles that classify lung cancer patients. J Pathol 244:421–431.
    OpenUrl
  31. ↵
    1. Nawaz S,
    2. Heindl A,
    3. Koelble K,
    4. Yuan Y
    (2015) Beyond immune density: Critical role of spatial heterogeneity in estrogen receptor-negative breast cancer. Mod Pathol 28:766–777.
    OpenUrlCrossRefPubMed
  32. ↵
    1. Andor N, et al.
    (2016) Pan-cancer analysis of the extent and consequences of intratumor heterogeneity. Nat Med 22:105–113.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Matsushita H, et al.
    (2012) Cancer exome analysis reveals a T-cell-dependent mechanism of cancer immunoediting. Nature 482:400–404.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Şenbabaoğlu Y, et al.
    (2016) Tumor immune microenvironment characterization in clear cell renal cell carcinoma identifies prognostic and immunotherapeutically relevant messenger RNA signatures. Genome Biol 17:231.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Simpson AJ,
    2. Caballero OL,
    3. Jungbluth A,
    4. Chen YT,
    5. Old LJ
    (2005) Cancer/testis antigens, gametogenesis and cancer. Nat Rev Cancer 5:615–625.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Almeida LG, et al.
    (2009) CTdatabase: A knowledge-base of high-throughput and curated data on cancer-testis antigens. Nucleic Acids Res 37:D816–D819.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Gaudet F, et al.
    (2004) Dnmt1 expression in pre- and postimplantation embryogenesis and the maintenance of IAP silencing. Mol Cell Biol 24:1640–1648.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Strissel PL, et al.
    (2012) Reactivation of codogenic endogenous retroviral (ERV) envelope genes in human endometrial carcinoma and prestages: Emergence of new molecular targets. Oncotarget 3:1204–1219.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Schmitt K,
    2. Reichrath J,
    3. Roesch A,
    4. Meese E,
    5. Mayer J
    (2013) Transcriptional profiling of human endogenous retrovirus group HERV-K(HML-2) loci in melanoma. Genome Biol Evol 5:307–328.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Chiappinelli KB, et al.
    (2015) Inhibiting DNA methylation causes an interferon response in cancer via dsRNA including endogenous retroviruses. Cell 162:974–986.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Boller K,
    2. Janssen O,
    3. Schuldes H,
    4. Tönjes RR,
    5. Kurth R
    (1997) Characterization of the antibody response specific for the human endogenous retrovirus HTDV/HERV-K. J Virol 71:4581–4588.
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Rooney MS,
    2. Shukla SA,
    3. Wu CJ,
    4. Getz G,
    5. Hacohen N
    (2015) Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell 160:48–61.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Yang BX, et al.
    (2015) Systematic identification of factors for provirus silencing in embryonic stem cells. Cell 163:230–245.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Barretina J, et al.
    (2012) The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature 483:603–607.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Lee YH, et al.
    (2017) Inhibition of the B7-H3 immune checkpoint limits tumor growth by enhancing cytotoxic lymphocyte function. Cell Res 27:1034–1045.
    OpenUrl
  46. ↵
    1. Mahnke K,
    2. Enk AH
    (2016) TIGIT-CD155 interactions in melanoma: A novel co-inhibitory pathway with potential for clinical intervention. J Invest Dermatol 136:9–11.
    OpenUrl
  47. ↵
    1. Johnston CJ,
    2. Smyth DJ,
    3. Dresser DW,
    4. Maizels RM
    (2016) TGF-β in tolerance, development and regulation of immunity. Cell Immunol 299:14–22.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Blokzijl F, et al.
    (2016) Tissue-specific mutation accumulation in human adult stem cells during life. Nature 538:260–264.
    OpenUrlCrossRefPubMed
  49. ↵
    1. Saldanha-Araujo F, et al.
    (2010) Cancer/testis antigen expression on mesenchymal stem cells isolated from different tissues. Anticancer Res 30:5023–5027.
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Lifantseva N, et al.
    (2011) Expression patterns of cancer-testis antigens in human embryonic stem cells and their cell derivatives indicate lineage tracks. Stem Cells Int 2011:795239.
    OpenUrlPubMed
  51. ↵
    1. Yamada R, et al.
    (2013) Preferential expression of cancer/testis genes in cancer stem-like cells: Proposal of a novel sub-category, cancer/testis/stem gene. Tissue Antigens 81:428–434.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Guo YL, et al.
    (2015) Attenuated innate immunity in embryonic stem cells and its implications in developmental biology and regenerative medicine. Stem Cells 33:3165–3173.
    OpenUrlCrossRef
  53. ↵
    1. Fuertes MB,
    2. Woo SR,
    3. Burnett B,
    4. Fu YX,
    5. Gajewski TF
    (2013) Type I interferon response and innate immune sensing of cancer. Trends Immunol 34:67–73.
    OpenUrlCrossRefPubMed
  54. ↵
    1. Seaman S, et al.
    (2017) Eradication of tumors through simultaneous ablation of CD276/B7-H3-positive tumor cells and tumor vasculature. Cancer cell 31:501–515.e8.
    OpenUrl
  55. ↵
    1. Bin Z, et al.
    (2014) Overexpression of B7-H3 in CD133+ colorectal cancer cells is associated with cancer progression and survival in human patients. J Surg Res 188:396–403.
    OpenUrlCrossRefPubMed
  56. ↵
    1. Solecki DJ,
    2. Gromeier M,
    3. Mueller S,
    4. Bernhardt G,
    5. Wimmer E
    (2002) Expression of the human poliovirus receptor/CD155 gene is activated by sonic hedgehog. J Biol Chem 277:25697–25702.
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Cochrane CR,
    2. Szczepny A,
    3. Watkins DN,
    4. Cain JE
    (2015) Hedgehog signaling in the maintenance of cancer stem cells. Cancers (Basel) 7:1554–1585.
    OpenUrl
  58. ↵
    1. Watabe T,
    2. Miyazono K
    (2009) Roles of TGF-beta family signaling in stem cell renewal and differentiation. Cell Res 19:103–115.
    OpenUrlCrossRefPubMed
  59. ↵
    1. Peñuelas S, et al.
    (2009) TGF-beta increases glioma-initiating cell self-renewal through the induction of LIF in human glioblastoma. Cancer Cell 15:315–327.
    OpenUrlCrossRefPubMed
  60. ↵
    1. Ikushima H, et al.
    (2009) Autocrine TGF-beta signaling maintains tumorigenicity of glioma-initiating cells through Sry-related HMG-box factors. Cell Stem Cell 5:504–514.
    OpenUrlCrossRefPubMed
  61. ↵
    1. Spranger S,
    2. Gajewski TF
    (2018) Impact of oncogenic pathways on evasion of antitumour immune responses. Nat Rev Cancer 18:139–147.
    OpenUrlCrossRef
  62. ↵
    1. Nusse R
    (2008) Wnt signaling and stem cell control. Cell Res 18:523–527.
    OpenUrlCrossRefPubMed
  63. ↵
    1. Murphy MJ,
    2. Wilson A,
    3. Trumpp A
    (2005) More than just proliferation: Myc function in stem cells. Trends Cell Biol 15:128–137.
    OpenUrlCrossRefPubMed
  64. ↵
    1. Hill R,
    2. Wu H
    (2009) PTEN, stem cells, and cancer stem cells. J Biol Chem 284:11755–11759.
    OpenUrlAbstract/FREE Full Text
  65. ↵
    1. Gurumurthy S, et al.
    (2010) The Lkb1 metabolic sensor maintains haematopoietic stem cell survival. Nature 468:659–663.
    OpenUrlCrossRefPubMed
  66. ↵
    1. de Thé H
    (2018) Differentiation therapy revisited. Nat Rev Cancer 18:117–127.
    OpenUrl
  67. ↵
    1. Liu J, et al.
    (2018) An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics. Cell 173:400–416.e11.
    OpenUrlCrossRefPubMed
  68. ↵
    1. Colaprico A, et al.
    (2016) TCGAbiolinks: An R/bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res 44:e71.
    OpenUrlCrossRefPubMed
  69. ↵
    1. Durinck S,
    2. Spellman PT,
    3. Birney E,
    4. Huber W
    (2009) Mapping identifiers for the integration of genomic datasets with the R/bioconductor package biomaRt. Nat Protoc 4:1184–1191.
    OpenUrlCrossRefPubMed
  70. ↵
    1. Hänzelmann S,
    2. Castelo R,
    3. Guinney J
    (2013) GSVA: Gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics 14:7.
    OpenUrlCrossRefPubMed
  71. ↵
    1. Therneau TM
    (2018) survival: Survival Analysis. R Package Version 2.38. Available at https://CRAN.R-project.org/package=survival. Accessed on December 20, 2018.
  72. ↵
    1. Bates D,
    2. Mächler M,
    3. Bolker BM,
    4. Walker SC
    (2015) Fitting linear mixed-effects models using lme4. J Stat Softw 67:1–48.
    OpenUrlCrossRefPubMed
  73. ↵
    1. Yu G,
    2. He QY
    (2016) ReactomePA: An R/bioconductor package for reactome pathway analysis and visualization. Mol Biosyst 12:477–479.
    OpenUrlCrossRefPubMed
  74. ↵
    1. Ritchie ME, et al.
    (2015) limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 43:e47.
    OpenUrlCrossRefPubMed
  75. ↵
    1. Gentleman R,
    2. Carey V,
    3. Huber W,
    4. Hahne F
    (2018) genefilter: Methods for Filtering Genes from High-Throughput Experiments. R Package Version 1.64.0. Available at bioconductor.org/packages/genefilter. Accessed on December 20, 2018.
  76. ↵
    1. Oksanen J, et al.
    (2015) vegan: Community Ecology Package. R Package Version 2.5-4. Available at https://CRAN.R-project.org/package=vegan. Accessed on December 20, 2018.
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Cancer stemness, intratumoral heterogeneity, and immune response across cancers
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Cancer stemness, intratumoral heterogeneity, and immune response across cancers
Alex Miranda, Phineas T. Hamilton, Allen W. Zhang, Swetansu Pattnaik, Etienne Becht, Artur Mezheyeuski, Jarle Bruun, Patrick Micke, Aurélien de Reynies, Brad H. Nelson
Proceedings of the National Academy of Sciences Apr 2019, 116 (18) 9020-9029; DOI: 10.1073/pnas.1818210116

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Cancer stemness, intratumoral heterogeneity, and immune response across cancers
Alex Miranda, Phineas T. Hamilton, Allen W. Zhang, Swetansu Pattnaik, Etienne Becht, Artur Mezheyeuski, Jarle Bruun, Patrick Micke, Aurélien de Reynies, Brad H. Nelson
Proceedings of the National Academy of Sciences Apr 2019, 116 (18) 9020-9029; DOI: 10.1073/pnas.1818210116
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Biological Sciences
  • Medical Sciences
Proceedings of the National Academy of Sciences: 116 (18)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Model of the Amazon forest
News Feature: A sea in the Amazon
Did the Caribbean sweep into the western Amazon millions of years ago, shaping the region’s rich biodiversity?
Image credit: Tacio Cordeiro Bicudo (University of São Paulo, São Paulo, Brazil), Victor Sacek (University of São Paulo, São Paulo, Brazil), and Lucy Reading-Ikkanda (artist).
Syrian archaeological site
Journal Club: In Mesopotamia, early cities may have faltered before climate-driven collapse
Settlements 4,200 years ago may have suffered from overpopulation before drought and lower temperatures ultimately made them unsustainable.
Image credit: Andrea Ricci.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490