Despite biomarker stratification, the anti-EGFR antibody cetuximab is only effective against a subgroup of colorectal cancers (CRCs). This genomic and transcriptomic analysis of the cetuximab ...resistance landscape in 35 RAS wild-type CRCs identified associations of NF1 and non-canonical RAS/RAF aberrations with primary resistance and validated transcriptomic CRC subtypes as non-genetic predictors of benefit. Sixty-four percent of biopsies with acquired resistance harbored no genetic resistance drivers. Most of these had switched from a cetuximab-sensitive transcriptomic subtype at baseline to a fibroblast- and growth factor-rich subtype at progression. Fibroblast-supernatant conferred cetuximab resistance in vitro, confirming a major role for non-genetic resistance through stromal remodeling. Cetuximab treatment increased cytotoxic immune infiltrates and PD-L1 and LAG3 immune checkpoint expression, potentially providing opportunities to treat cetuximab-resistant CRCs with immunotherapy.
•NF1 and non-canonical KRAS and BRAF aberrations associate with cetuximab resistance•Genetic resistance drivers are absent in most biopsies that acquired resistance•Stromal remodeling is an alternative non-genetic mechanism of cetuximab resistance•Cetuximab-mediated immune modulation may sensitize CRCs to immunotherapy
Woolston et al. show that in metastatic colorectal cancer cetuximab resistance can be conferred by genetic mechanisms, such as NF1 loss or RAS/RAF alterations, or by transcriptomic changes that induce a stroma-rich phenotype. They also provide a rationale for combining cetuximab with immunotherapy.
Previously, we classified colorectal cancers (CRCs) into five CRCAssigner (CRCA) subtypes with different prognoses and potential treatment responses, later consolidated into four consensus molecular ...subtypes (CMS). Here we demonstrate the analytical development and validation of a custom NanoString nCounter platform-based biomarker assay (NanoCRCA) to stratify CRCs into subtypes. To reduce costs, we switched from the standard nCounter protocol to a custom modified protocol. The assay included a reduced 38-gene panel that was selected using an in-house machine-learning pipeline. We applied NanoCRCA to 413 samples from 355 CRC patients. From the fresh frozen samples (n = 237), a subset had matched microarray/RNAseq profiles (n = 47) or formalin-fixed paraffin-embedded (FFPE) samples (n = 58). We also analyzed a further 118 FFPE samples. We compared the assay results with the CMS classifier, different platforms (microarrays/RNAseq) and gene-set classifiers (38 and the original 786 genes). The standard and modified protocols showed high correlation (> 0.88) for gene expression. Technical replicates were highly correlated (> 0.96). NanoCRCA classified fresh frozen and FFPE samples into all five CRCA subtypes with consistent classification of selected matched fresh frozen/FFPE samples. We demonstrate high and significant subtype concordance across protocols (100%), gene sets (95%), platforms (87%) and with CMS subtypes (75%) when evaluated across multiple datasets. Overall, our NanoCRCA assay with further validation may facilitate prospective validation of CRC subtypes in clinical trials and beyond.
•By comprehensively profiling the immune transcriptome of triple-negative breast cancers (TNBCs) and applying machine learning methods, we identified three immunotypes (1, 2 and 3) with significant ...differences in prognosis and therapy responses.•TNBC is much more common in India than in the West. However, there was no significant difference in immune transcriptome between Indian and Western TNBCs in our study.•Our analysis shows biological, signaling and clinical features of three TNBC immunotypes.•Immunotype-1 gene signature is associated with improved prognosis and treatment responses in a cross-cancer comparison analysis of melanoma patients treated with anti-PDL1 therapy.•Our study identified a potential opportunity to stratify patients for MAGEA3 and anti-PDL1-based therapies, which warrants further validation.
Triple-negative breast cancer (TNBC) is a heterogeneous disease with a significant challenge to effectively manage in the clinic worldwide. Immunotherapy may be beneficial to TNBC patients if responders can be effectively identified. Here we sought to elucidate the immune landscape of TNBCs by stratifying patients into immune-specific subtypes (immunotypes) to decipher the molecular and cellular presentations and signaling events of this heterogeneous disease and associating them with their clinical outcomes and potential treatment options.
We profiled 730 immune genes in 88 retrospective Indian TNBC samples using the NanoString platform, established immunotypes using non-negative matrix factorization-based machine learning approach, and validated them using Western TNBCs (n=422; public datasets). Immunotype-specific gene signatures were associated with clinicopathological features, immune cell types, biological pathways, acute/chronic inflammatory responses, and immunogenic cell death processes. Responses to different immunotherapies associated with TNBC immunotypes were assessed using cross-cancer comparison to melanoma (n=504). Tumor-infiltrating lymphocytes (TILs) and pan-macrophage spatial marker expression were evaluated.
We identified three robust transcriptome-based immunotypes in both Indian and Western TNBCs in similar proportions. Immunotype-1 tumors, mainly representing well-known claudin-low and immunomodulatory subgroups, harbored dense TIL infiltrates and T-helper-1 (Th1) response profiles associated with smaller tumors, pre-menopausal status, and a better prognosis. They displayed a cascade of events, including acute inflammation, damage-associated molecular patterns, T-cell receptor-related and chemokine-specific signaling, antigen presentation, and viral-mimicry pathways. On the other hand, immunotype-2 was enriched for Th2/Th17 responses, CD4+ regulatory cells, basal-like/mesenchymal immunotypes, and an intermediate prognosis. In contrast to the two T-cell enriched immunotypes, immunotype-3 patients expressed innate immune genes/proteins, including those representing myeloid infiltrations (validated by spatial immunohistochemistry), and had poor survival. Remarkably, a cross-cancer comparison analysis revealed the association of immunotype-1 with responses to anti-PD-L1 and MAGEA3 immunotherapies.
Overall, the TNBC immunotypes identified in TNBCs reveal different prognoses, immune infiltrations, signaling, acute/chronic inflammation leading to immunogenic cell death of cancer cells, and potentially distinct responses to immunotherapies. The overlap in immune characteristics in Indian and Western TNBCs suggests similar efficiency of immunotherapy in both populations if strategies to select patients according to immunotypes can be further optimized and implemented.
BackgroundRectal cancers show a highly varied response to neoadjuvant radiotherapy/chemoradiation (RT/CRT) and the impact of the tumor immune microenvironment on this response is poorly understood. ...Current clinical tumor regression grading systems attempt to measure radiotherapy response but are subject to interobserver variation. An unbiased and unique histopathological quantification method (change in tumor cell density (ΔTCD)) may improve classification of RT/CRT response. Furthermore, immune gene expression profiling (GEP) may identify differences in expression levels of genes relevant to different radiotherapy responses: (1) at baseline between poor and good responders, and (2) longitudinally from preradiotherapy to postradiotherapy samples. Overall, this may inform novel therapeutic RT/CRT combination strategies in rectal cancer.MethodsWe generated GEPs for 53 patients from biopsies taken prior to preoperative radiotherapy. TCD was used to assess rectal tumor response to neoadjuvant RT/CRT and ΔTCD was subjected to k-means clustering to classify patients into different response categories. Differential gene expression analysis was performed using statistical analysis of microarrays, pathway enrichment analysis and immune cell type analysis using single sample gene set enrichment analysis. Immunohistochemistry was performed to validate specific results. The results were validated using 220 pretreatment samples from publicly available datasets at metalevel of pathway and survival analyses.ResultsΔTCD scores ranged from 12.4% to −47.7% and stratified patients into three response categories. At baseline, 40 genes were significantly upregulated in poor (n=12) versus good responders (n=21), including myeloid and stromal cell genes. Of several pathways showing significant enrichment at baseline in poor responders, epithelial to mesenchymal transition, coagulation, complement activation and apical junction pathways were validated in external cohorts. Unlike poor responders, good responders showed longitudinal (preradiotherapy vs postradiotherapy samples) upregulation of 198 immune genes, reflecting an increased T-cell-inflamed GEP, type-I interferon and macrophage populations. Longitudinal pathway analysis suggested viral-like pathogen responses occurred in post-treatment resected samples compared with pretreatment biopsies in good responders.ConclusionThis study suggests potentially druggable immune targets in poor responders at baseline and indicates that tumors with a good RT/CRT response reprogrammed from immune “cold” towards an immunologically “hot” phenotype on treatment with radiotherapy.
Genome projects now generate large-scale data often produced at various time points by different laboratories using multiple platforms. This increases the potential for batch effects. Currently there ...are several batch evaluation methods like principal component analysis (PCA; mostly based on visual inspection), and sometimes they fail to reveal all of the underlying batch effects. These methods can also lead to the risk of unintentionally correcting biologically interesting factors attributed to batch effects. Here we propose a novel statistical method, finding batch effect (findBATCH), to evaluate batch effect based on probabilistic principal component and covariates analysis (PPCCA). The same framework also provides a new approach to batch correction, correcting batch effect (correctBATCH), which we have shown to be a better approach to traditional PCA-based correction. We demonstrate the utility of these methods using two different examples (breast and colorectal cancers) by merging gene expression data from different studies after diagnosing and correcting for batch effects and retaining the biological effects. These methods, along with conventional visual inspection-based PCA, are available as a part of an R package exploring batch effect (exploBATCH; https://github.com/syspremed/exploBATCH ).
Abstract
Interleukin-11 (IL-11) is a pleiotropic cytokine that belongs to gp130 family. It plays a significant role in the synthesis and maturation of hematopoietic cells, inhibition of adipogenesis, ...regulation of embryo implantation, and trophoblasts invasion. Although IL-11 signaling has been described in several biological processes, a centralized resource documenting these molecular reactions induced by IL-11 is not publicly available. In the current study, we have manually annotated the molecular reactions and interactions induced by IL-11 from literature available. We have documented 40 unique molecules involved in 18 protein-protein interactions, 26 enzyme-substrate reactions, 7 translocation events, and 4 activation/ inhibition reactions. We have also annotated 23 genes reported to be differentially regulated under IL-11 stimulation. We have enabled the data availability in standard exchange formats from 'NetPath', a repository for signaling pathways. We believe that this will help in the identification of potential therapeutic targets in IL-11-associated disorders.
A comprehensive analysis of the immune landscape of pancreatic neuroendocrine tumours (PanNETs) was performed according to clinicopathological parameters and previously defined molecular subtypes to ...identify potential therapeutic vulnerabilities in this disease.
Differential expression analysis of 600 immune-related genes was performed on 207 PanNET samples, comprising a training cohort (n=72) and two validation cohorts (n=135) from multiple transcriptome profiling platforms. Different immune-related and subtype-related phenotypes, cell types and pathways were investigated using different in silico methods and were further validated using spatial multiplex immunofluorescence.
The study identified an immune signature of 132 genes segregating PanNETs (n=207) according to four previously defined molecular subtypes: metastasis-like primary (MLP)-1 and MLP-2, insulinoma-like and intermediate. The MLP-1 subtype (26%-31% samples across three cohorts) was strongly associated with elevated levels of immune-related genes, poor prognosis and a cascade of tumour evolutionary events: larger hypoxic and necroptotic tumours leading to increased damage-associated molecular patterns (viral mimicry), stimulator of interferon gene pathway, T cell-inflamed genes, immune checkpoint targets, and T cell-mediated and M1 macrophage-mediated immune escape mechanisms. Multiplex spatial profiling validated significantly increased macrophages in the MLP-1 subtype.
This study provides novel data on the immune microenvironment of PanNETs and identifies MLP-1 subtype as an immune-high phenotype featuring a broad and robust activation of immune-related genes. This study, with further refinement, paves the way for future precision immunotherapy studies in PanNETs to potentially select a subset of MLP-1 patients who may be more likely to respond.
Breast cancer is a highly heterogeneous disease. Although differences between intrinsic breast cancer subtypes have been well studied, heterogeneity within each subtype, especially luminal-A cancers, ...requires further interrogation to personalize disease management. Here, we applied well-characterized and cancer-associated heterocellular signatures representing stem, mesenchymal, stromal, immune, and epithelial cell types to breast cancer. This analysis stratified the luminal-A breast cancer samples into five subtypes with a majority of them enriched for a subtype (stem-like) that has increased stem and stromal cell gene signatures, representing potential luminal progenitor origin. The enrichment of immune checkpoint genes and other immune cell types in two (including stem-like) of the five heterocellular subtypes of luminal-A tumors suggest their potential response to immunotherapy. These immune-enriched subtypes of luminal-A tumors (containing only estrogen receptor positive samples) showed good or intermediate prognosis along with the two other differentiated subtypes as assessed using recurrence-free and distant metastasis-free patient survival outcomes. On the other hand, a partially differentiated subtype of luminal-A breast cancer with transit-amplifying colon-crypt characteristics showed poor prognosis. Furthermore, published luminal-A subtypes associated with specific somatic copy number alterations and mutations shared similar cellular and mutational characteristics to colorectal cancer subtypes where the heterocellular signatures were derived. These heterocellular subtypes reveal transcriptome and cell-type based heterogeneity of luminal-A and other breast cancer subtypes that may be useful for additional understanding of the cancer type and potential patient stratification and personalized medicine.