Immune landscape, evolution, hypoxia-mediated viral mimicry pathways and therapeutic potential in molecular subtypes of pancreatic neuroendocrine tumours.

OBJECTIVE: 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. DESIGN: 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. RESULTS: 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. CONCLUSION: 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.

Inter- and intra-tumoural heterogeneity in cancer-associated fibroblasts of human pancreatic ductal adenocarcinoma.

Cancer-associated fibroblasts (CAF) are orchestrators of the pancreatic ductal adenocarcinoma (PDAC) microenvironment. Stromal heterogeneity may explain differential pathophysiological roles of the stroma (pro- versus anti-tumoural) in PDAC. We hypothesised that multiple CAF functional subtypes exist in PDAC, that contribute to stromal heterogeneity through interactions with cancer cells. Using molecular and functional analysis of patient-derived CAF primary cultures, we demonstrated that human PDAC-derived CAFs display a high level of inter- and intra-tumour heterogeneity. We identified at least four subtypes of CAFs based on transcriptomic analysis, and propose a classification for human PDAC-derived CAFs (pCAFassigner). Multiple CAF subtypes co-existed in individual patient samples. The presence of these CAF subtypes in bulk tumours was confirmed using publicly available gene expression profiles, and immunostainings of CAF subtype markers. Each subtype displayed specific phenotypic features (matrix- and immune-related signatures, vimentin and α-smooth muscle actin expression, proliferation rate), and was associated with an assessable prognostic impact. A prolonged exposure of non-tumoural pancreatic stellate cells to conditioned media from cancer cell lines (cancer education experiment) induced a CAF-like phenotype, including loss of capacity to revert to quiescence and an increase in the expression of genes related to CAF subtypes B and C. This classification demonstrates molecular and functional inter- and intra-tumoural heterogeneity of CAFs in human PDAC. Our subtypes overlap with those identified from single-cell analyses in other cancers, and pave the way for the development of therapies targeting specific CAF subpopulations in PDAC. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Identifying signatures of positive selection in pigmentation genes in two South Asian populations.

OBJECTIVES: Skin pigmentation is a polygenic trait showing wide phenotypic variations among global populations. While numerous pigmentation genes have been identified to be under positive selection among European and East populations, genes contributing to phenotypic variation in skin pigmentation within and among South Asian populations are still poorly understood. The present study uses data from the Phase 3 of the 1000 genomes project focusing on two South Asian populations-GIH (Gujarati Indian from Houston, Texas) and ITU (Indian Telugu from UK), so as to decode the genetic architecture involved in adaptation to ultraviolet radiation in South Asian populations. METHODS: Statistical tests included were (1) tests to identify deviations of the Site Frequency Spectrum (SFS) from neutral expectations (Tajima's D, Fay and Wu's H and Fu and Li's D* and F*), (2) tests focused on the identification of high-frequency haplotypes with extended linkage disequilibrium (iHS and Rsb), and (3) tests based on genetic differentiation between populations (LSBL). RESULTS: Twenty-two pigmentation genes fall in the top 1% for at least one statistic in the GIH population, 5 of which (LYST, OCA2, SLC24A5, SLC45A2, and TYR) have been previously associated with normal variation in skin, hair, or eye color. In comparison, 17 genes fall in the top 1% for at least one statistic in the ITU population. Twelve loci which are identified as outliers in the ITU scan were also identified in the GIH population. CONCLUSIONS: These results suggest that selection may have affected these loci broadly across the region.

Comprehensive characterization of immune landscape of Indian and Western triple negative breast cancers.

PURPOSE: 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. EXPERIMENTAL DESIGN: 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. RESULTS: 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. CONCLUSION: 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.

Differential and longitudinal immune gene patterns associated with reprogrammed microenvironment and viral mimicry in response to neoadjuvant radiotherapy in rectal cancer.

BACKGROUND: Rectal 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. METHODS: We 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. CONCLUSION: This 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.

Heterocellular gene signatures reveal luminal-A breast cancer heterogeneity and differential therapeutic responses.

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.

Analytical Validation of Multiplex Biomarker Assay to Stratify Colorectal Cancer into Molecular Subtypes.

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.

Genomic and Transcriptomic Determinants of Therapy Resistance and Immune Landscape Evolution during Anti-EGFR Treatment in Colorectal Cancer.

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.

A Novel Statistical Method to Diagnose, Quantify and Correct Batch Effects in Genomic Studies.

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 ).

IL-11/IL11RA receptor mediated signaling: a web accessible knowledgebase.

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.