An integrative gene expression signature analysis identifies CMS4 KRAS-mutated colorectal cancers sensitive to combined MEK and SRC targeted therapy.

BACKGROUND: Over half of colorectal cancers (CRCs) are hard-wired to RAS/RAF/MEK/ERK pathway oncogenic signaling. However, the promise of targeted therapeutic inhibitors, has been tempered by disappointing clinical activity, likely due to complex resistance mechanisms that are not well understood. This study aims to investigate MEK inhibitor-associated resistance signaling and identify subpopulation(s) of CRC patients who may be sensitive to biomarker-driven drug combination(s). METHODS: We classified 2250 primary and metastatic human CRC tumors by consensus molecular subtypes (CMS). For each tumor, we generated multiple gene expression signature scores measuring MEK pathway activation, MEKi "bypass" resistance, SRC activation, dasatinib sensitivity, EMT, PC1, Hu-Lgr5-ISC, Hu-EphB2-ISC, Hu-Late TA, Hu-Proliferation, and WNT activity. We carried out correlation, survival and other bioinformatic analyses. Validation analyses were performed in two independent publicly available CRC tumor datasets (n = 585 and n = 677) and a CRC cell line dataset (n = 154). RESULTS: Here we report a central role of SRC in mediating "bypass"-resistance to MEK inhibition (MEKi), primarily in cancer stem cells (CSCs). Our integrated and comprehensive gene expression signature analyses in 2250 CRC tumors reveal that MEKi-resistance is strikingly-correlated with SRC activation (Spearman P < 10-320), which is similarly associated with EMT (epithelial to mesenchymal transition), regional metastasis and disease recurrence with poor prognosis. Deeper analysis shows that both MEKi-resistance and SRC activation are preferentially associated with a mesenchymal CSC phenotype. This association is validated in additional independent CRC tumor and cell lines datasets. The CMS classification analysis demonstrates the strikingly-distinct associations of CMS1-4 subtypes with the MEKi-resistance and SRC activation. Importantly, MEKi + SRCi sensitivities are predicted to occur predominantly in the KRAS mutant, mesenchymal CSC-like CMS4 CRCs. CONCLUSIONS: Large human tumor gene expression datasets representing CRC heterogeneity can provide deep biological insights heretofore not possible with cell line models, suggesting novel repurposed drug combinations. We identified SRC as a common targetable node--an Achilles' heel--in MEKi-targeted therapy-associated resistance in mesenchymal stem-like CRCs, which may help development of a biomarker-driven drug combination (MEKi + SRCi) to treat problematic subpopulations of CRC.

Dalotuzumab in chemorefractory KRAS exon 2 mutant colorectal cancer: Results from a randomised phase II/III trial.

Limited data are available on the efficacy of anti-IGF-1R agents in KRAS mutant colorectal cancer (CRC). We analysed the outcome of 69 chemorefractory, KRAS exon 2 mutant CRC patients who were enrolled in a double-blind, randomised, phase II/III study of irinotecan and cetuximab plus dalotuzumab 10 mg/kg once weekly (arm A), dalotuzumab 7.5 mg/kg every second week (arm B) or placebo (arm C). Objective response rate (5.6% vs. 3.1% vs. 4.8%), median progression-free survival (2.7 vs. 2.6 vs. 1.4 months) and overall survival (7.8 vs. 10.3 vs. 7.8 months) were not statistically significantly different between treatment groups. Most common grade ≥3 treatment-related toxicities included neutropenia, diarrhoea, hyperglycaemia, fatigue and dermatitis acneiform. Expression of IGF-1R, IGF-1, IGF-2 and EREG by quantitative real-time polymerase chain reaction was assessed in 351 patients from the same study with available data on KRAS exon 2 mutational status. Median cycle threshold values for all biomarkers were significantly lower (i.e., higher expression, p < 0.05) among patients with KRAS wild-type compared to those with KRAS exon 2 mutant tumours. No significant changes were found according to location of the primary tumour with only a trend towards lower expression of IGF-1 in colon compared to rectal cancers (p = 0.06). Albeit limited by the small sample size, this study does not appear to support a potential role for anti-IGF-1R agents in KRAS exon 2 mutant CRC. Data on IGF-1R, IGF-1 and IGF-2 expression here reported may be useful for patient stratification in future trials with inhibitors of the IGF pathway.

Whole-genome sequencing identifies recurrent mutations in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the most deadly cancers worldwide and has no effective treatment, yet the molecular basis of hepatocarcinogenesis remains largely unknown. Here we report findings from a whole-genome sequencing (WGS) study of 88 matched HCC tumor/normal pairs, 81 of which are Hepatitis B virus (HBV) positive, seeking to identify genetically altered genes and pathways implicated in HBV-associated HCC. We find beta-catenin to be the most frequently mutated oncogene (15.9%) and TP53 the most frequently mutated tumor suppressor (35.2%). The Wnt/beta-catenin and JAK/STAT pathways, altered in 62.5% and 45.5% of cases, respectively, are likely to act as two major oncogenic drivers in HCC. This study also identifies several prevalent and potentially actionable mutations, including activating mutations of Janus kinase 1 (JAK1), in 9.1% of patients and provides a path toward therapeutic intervention of the disease.

Functional genomics identifies therapeutic targets for MYC-driven cancer.

MYC oncogene family members are broadly implicated in human cancers, yet are considered "undruggable" as they encode transcription factors. MYC also carries out essential functions in proliferative tissues, suggesting that its inhibition could cause severe side effects. We elected to identify synthetic lethal interactions with c-MYC overexpression (MYC-SL) in a collection of ~3,300 druggable genes, using high-throughput siRNA screening. Of 49 genes selected for follow-up, 48 were confirmed by independent retesting and approximately one-third selectively induced accumulation of DNA damage, consistent with enrichment in DNA-repair genes by functional annotation. In addition, genes involved in histone acetylation and transcriptional elongation, such as TRRAP and BRD4, were identified, indicating that the screen revealed known MYC-associated pathways. For in vivo validation we selected CSNK1e, a kinase whose expression correlated with MYCN amplification in neuroblastoma (an established MYC-driven cancer). Using RNAi and available small-molecule inhibitors, we confirmed that inhibition of CSNK1e halted growth of MYCN-amplified neuroblastoma xenografts. CSNK1e had previously been implicated in the regulation of developmental pathways and circadian rhythms, whereas our data provide a previously unknown link with oncogenic MYC. Furthermore, expression of CSNK1e correlated with c-MYC and its transcriptional signature in other human cancers, indicating potential broad therapeutic implications of targeting CSNK1e function. In summary, through a functional genomics approach, pathways essential in the context of oncogenic MYC but not to normal cells were identified, thus revealing a rich therapeutic space linked to a previously "undruggable" oncogene.

Systems biology approach for new target and biomarker identification.

The pharmaceutical industry is spending increasingly large amounts of money on the discovery and development of novel medicines, but this investment is not adequately paying off in an increased rate of newly approved drugs by the FDA. The post-genomic era has provided a wealth of novel approaches for generating large, high-dimensional genetic and transcriptomic data sets from large cohorts of preclinical species as well as normal and diseased individuals. This systems biology approach to understanding disease-related biology is revolutionizing our understanding of the cellular pathways and gene networks underlying the onset of disease, and the mechanisms of pharmacological treatments that ameliorate disease phenotypes. In this article, we review a number of approaches being used by pharmaceutical and biotechnology companies, e.g., high-throughput DNA genotyping, sequencing, and genome-wide gene expression profiling, to enable drug discovery and development through the identification of new drug targets and biomarkers of disease progression, drug pharmacodynamics, and predictive markers for selecting the patients most likely to respond to therapy.

Evaluation of potential biomarkers for lenvatinib plus pembrolizumab among patients with advanced endometrial cancer: results from Study 111/KEYNOTE-146.

BACKGROUND: Lenvatinib plus pembrolizumab demonstrated clinically meaningful benefit in patients with previously treated advanced endometrial carcinoma in Study 111/KEYNOTE-146 (NCT02501096). In these exploratory analyses from this study, we evaluated the associations between clinical outcomes and gene expression signature scores and descriptively summarized response in biomarker subpopulations defined by tumor mutational burden (TMB) and DNA variants for individual genes of interest. METHODS: Patients with histologically confirmed metastatic endometrial carcinoma received oral lenvatinib 20 mg once daily plus intravenous pembrolizumab 200 mg every 3 weeks for 35 cycles. Archived formalin-fixed paraffin-embedded tissue was obtained from all patients. T-cell-inflamed gene expression profile (TcellinfGEP) and 11 other gene signatures were evaluated by RNA sequencing. TMB, hotspot mutations in PIK3CA (oncogene), and deleterious mutations in PTEN and TP53 (tumor suppressor genes) were evaluated by whole-exome sequencing (WES). RESULTS: 93 and 79 patients were included in the RNA-sequencing-evaluable and WES-evaluable populations, respectively. No statistically significant associations were observed between any of the RNA-sequencing signature scores and objective response rate or progression-free survival. Area under the receiver operating characteristic curve values for response ranged from 0.39 to 0.54; all 95% CIs included 0.50. Responses were seen regardless of TMB (≥175 or <175 mutations/exome) and mutation status. There were no correlations between TcellinfGEP and TMB, TcellinfGEP and microvessel density (MVD), or MVD and TMB. CONCLUSIONS: This analysis demonstrated efficacy for lenvatinib plus pembrolizumab regardless of biomarker status. Results from this study do not support clinical utility of the evaluated biomarkers. Further investigation of biomarkers for this regimen is warranted. TRIAL REGISTRATION NUMBER: NCT02501096.

Single-Cell Analysis Identifies NOTCH3-Mediated Interactions between Stromal Cells That Promote Microenvironment Remodeling and Invasion in Lung Adenocarcinoma.

Cancer immunotherapy has revolutionized the treatment of lung adenocarcinoma (LUAD); however, a significant proportion of patients do not respond. Recent transcriptomic studies to understand determinants of immunotherapy response have pinpointed stromal-mediated resistance mechanisms. To gain a better understanding of stromal biology at the cellular and molecular level in LUAD, we performed single-cell RNA sequencing of 256,379 cells, including 13,857 mesenchymal cells, from 9 treatment-naïve patients. Among the mesenchymal cell subsets, FAP+PDPN+ cancer-associated fibroblasts (CAF) and ACTA2+MCAM+ pericytes were enriched in tumors and differentiated from lung-resident fibroblasts. Imaging mass cytometry revealed that both subsets were topographically adjacent to the perivascular niche and had close spatial interactions with endothelial cells (EC). Modeling of ligand and receptor interactomes between mesenchymal and ECs identified that NOTCH signaling drives these cell-to-cell interactions in tumors, with pericytes and CAFs as the signal receivers and arterial and PLVAPhigh immature neovascular ECs as the signal senders. Either pharmacologically blocking NOTCH signaling or genetically depleting NOTCH3 levels in mesenchymal cells significantly reduced collagen production and suppressed cell invasion. Bulk RNA sequencing data demonstrated that NOTCH3 expression correlated with poor survival in stroma-rich patients and that a T cell-inflamed gene signature only predicted survival in patients with low NOTCH3. Collectively, this study provides valuable insights into the role of NOTCH3 in regulating tumor stroma biology, warranting further studies to elucidate the clinical implications of targeting NOTCH3 signaling. SIGNIFICANCE: NOTCH3 signaling activates tumor-associated mesenchymal cells, increases collagen production, and augments cell invasion in lung adenocarcinoma, suggesting its critical role in remodeling tumor stroma.

Associations of Tissue Tumor Mutational Burden and Mutational Status With Clinical Outcomes With Pembrolizumab Plus Chemotherapy Versus Chemotherapy For Metastatic NSCLC.

Introduction: We evaluated tissue tumor mutational burden (tTMB) and mutations in STK11, KEAP1, and KRAS as biomarkers for outcomes with pembrolizumab plus platinum-based chemotherapy (pembrolizumab-combination) for NSCLC among patients in the phase 3 KEYNOTE-189 (ClinicalTrials.gov, NCT02578680; nonsquamous) and KEYNOTE-407 (ClinicalTrials.gov, NCT02775435; squamous) trials. Methods: This retrospective exploratory analysis evaluated prevalence of high tTMB and STK11, KEAP1, and KRAS mutations in patients enrolled in KEYNOTE-189 and KEYNOTE-407 and the relationship between these potential biomarkers and clinical outcomes. tTMB and STK11, KEAP1, and KRAS mutation status was assessed using whole-exome sequencing in patients with available tumor and matched normal DNA. The clinical utility of tTMB was assessed using a prespecified cutpoint of 175 mutations/exome. Results: Among patients with evaluable data from whole-exome sequencing for evaluation of tTMB (KEYNOTE-189, n = 293; KEYNOTE-407, n = 312) and matched normal DNA, no association was found between continuous tTMB score and overall survival (OS) or progression-free survival for pembrolizumab-combination (Wald test, one-sided p > 0.05) or placebo-combination (Wald test, two-sided p > 0.05) in patients with squamous or nonsquamous histology. Pembrolizumab-combination improved outcomes for patients with tTMB greater than or equal to 175 compared with tTMB less than 175 mutations/exome in KEYNOTE-189 (OS, hazard ratio = 0.64 [95% confidence interval (CI): 0.38‒1.07] and 0.64 [95% CI: 0.42‒0.97], respectively) and KEYNOTE-407 (OS, hazard ratio = 0.74 [95% CI: 0.50‒1.08 and 0.86 [95% CI: 0.57‒1.28], respectively) versus placebo-combination. Treatment outcomes were similar regardless of KEAP1, STK11, or KRAS mutation status. Conclusions: These findings support pembrolizumab-combination as first-line treatment in patients with metastatic NSCLC and do not suggest the utility of tTMB, STK11, KEAP1, or KRAS mutation status as a biomarker for this regimen.

Association between gene expression signatures and clinical outcomes of pembrolizumab versus paclitaxel in advanced gastric cancer: exploratory analysis from the randomized, controlled, phase III KEYNOTE-061 trial.

BACKGROUND: In the randomized, controlled, phase III KEYNOTE-061 trial, second-line pembrolizumab did not significantly prolong overall survival (OS) versus paclitaxel in patients with PD-L1-positive (combined positive score ≥1) advanced gastric/gastroesophageal junction (G/GEJ) cancer but did elicit a longer duration of response and offered a favorable safety profile. This prespecified exploratory analysis was conducted to evaluate associations between tumor gene expression signatures and clinical outcomes in the phase III KEYNOTE-061 trial. METHODS: Using RNA sequencing data obtained from formalin-fixed, paraffin-embedded baseline tumor tissue samples, we evaluated the 18-gene T-cell-inflamed gene expression profile (TcellinfGEP) and 10 non-TcellinfGEP signatures (angiogenesis, glycolysis, granulocytic myeloid-derived suppressor cell (gMDSC), hypoxia, monocytic MDSC (mMDSC), MYC, proliferation, RAS, stroma/epithelial-to-mesenchymal transition/transforming growth factor-ß, WNT). The association between each signature on a continuous scale and outcomes was analyzed using logistic (objective response rate (ORR)) and Cox proportional hazards regression (progression-free survival (PFS) and OS). One-sided (pembrolizumab) and two-sided (paclitaxel) p values were calculated for TcellinfGEP (prespecified α=0.05) and the 10 non-TcellinfGEP signatures (multiplicity-adjusted; prespecified α=0.10). RESULTS: RNA sequencing data were available for 137 patients in each treatment group. TcellinfGEP was positively associated with ORR (p=0.041) and PFS (p=0.026) for pembrolizumab but not paclitaxel (p>0.05). The TcellinfGEP-adjusted mMDSC signature was negatively associated with ORR (p=0.077), PFS (p=0.057), and OS (p=0.033) for pembrolizumab, while the TcellinfGEP-adjusted glycolysis (p=0.018), MYC (p=0.057), and proliferation (p=0.002) signatures were negatively associated with OS for paclitaxel. CONCLUSIONS: This exploratory analysis of tumor TcellinfGEP showed associations with ORR and PFS for pembrolizumab but not for paclitaxel. TcellinfGEP-adjusted mMDSC signature was negatively associated with ORR, PFS, and OS for pembrolizumab but not paclitaxel. These data suggest myeloid-driven suppression may play a role in resistance to PD-1 inhibition in G/GEJ cancer and support a strategy of considering immunotherapy combinations which target this myeloid axis. TRIAL REGISTRATION NUMBER: NCT02370498.

Unique ectopic lymph node-like structures present in human primary colorectal carcinoma are identified by immune gene array profiling.

We hypothesized that immune gene-related signatures would predict the presence of unique histological features of lymphoid cell infiltrates in colorectal carcinoma (CRC) that correlate with clinical parameters. Metagene analysis with gene chip technology was performed on 326 CRCs, which were then sorted by low versus high gene scores. Microscopically, CRCs with a high gene score revealed a marked host immune response organized, remarkably, as lymphoid follicles. Proliferation involved both B and T cells. In every case, the presence of CD79a(+) B-cell precursors was identified, suggesting that the lymphoid follicles represent newly formed, ectopic lymph node-like structures. CD21(+) dendritic cells were present within the follicular germinal centers, and CD3(+) T cells were localized mainly in the parafollicular cortex zone surrounding the B-cell area of the follicles. A strong correlation between a 12-chemokine gene subset of the molecular profile and the presence of ectopic lymph node-like structures was associated with better patient survival independent of tumor staging, site location, microsatellite instability or stability, and patient treatment. These findings suggest beneficial, intratumoral immune cell priming and raise the possibility of immunotherapy intervention decisions based on molecular signatures that can identify the presence of tumor-localized, ectopic lymph node-like structures.

Tumor mutational burden predicts the efficacy of pembrolizumab monotherapy: a pan-tumor retrospective analysis of participants with advanced solid tumors.

BACKGROUND: Several studies have evaluated the relationship between tumor mutational burden (TMB) and outcomes of immune checkpoint inhibitors. In the phase II KEYNOTE-158 study of pembrolizumab monotherapy for previously treated recurrent or metastatic cancer, high TMB as assessed by the FoundationOne CDx was associated with an improved objective response rate (ORR). METHODS: We retrospectively assessed the relationship between TMB and efficacy in participants with previously treated advanced solid tumors enrolled in 12 trials that evaluated pembrolizumab monotherapy, including 3 randomized trials that compared pembrolizumab with chemotherapy. TMB was assessed in formalin-fixed, paraffin-embedded pretreatment tumor samples by whole-exome sequencing. High TMB was defined as ≥175 mutations/exome. Microsatellite instability (MSI) phenotype was based on whole-exome sequencing results. Programmed death ligand 1 (PD-L1) expression was assessed by immunohistochemistry. The primary end point was ORR assessed per RECIST V.1.1 by independent central review. Other end points included progression-free survival (PFS) assessed per RECIST V.1.1 by independent central review and overall survival (OS). RESULTS: Of the 2234 participants in the analysis, 1772 received pembrolizumab monotherapy and 462 received chemotherapy. Among the pembrolizumab-treated participants, ORR was 31.4% (95% CI 27.1 to 36.0) in the 433 participants with TMB ≥175 mutations/exome and 9.5% (95% CI 8.0 to 11.2) in the 1339 participants with TMB <175 mutations/exome. The association of TMB with ORR was observed regardless of PD-L1 expression and not driven by specific tumor types or participants with very high TMB or high MSI. In the 3 randomized controlled trials, TMB was associated with ORR (p≤0.016), PFS (p≤0.005), and OS (p≤0.029) of pembrolizumab but not of chemotherapy (p≥0.340, p≥0.643, and p≥0.174, respectively), and pembrolizumab improved efficacy versus chemotherapy in participants with TMB ≥175 mutations/exome. CONCLUSIONS: TMB ≥175 mutations/exome is associated with clinically meaningful improvement in the efficacy of pembrolizumab monotherapy and improved outcomes for pembrolizumab versus chemotherapy across a wide range of previously treated advanced solid tumor types. These data suggest TMB has broad clinical utility irrespective of tumor type, PD-L1 expression, or MSI status and support its use as a predictive biomarker for pembrolizumab monotherapy in participants with previously treated advanced solid tumors.

Transcriptomic Determinants of Response to Pembrolizumab Monotherapy across Solid Tumor Types.

PURPOSE: To explore relationships between biological gene expression signatures and pembrolizumab response. EXPERIMENTAL DESIGN: RNA-sequencing data on baseline tumor tissue from 1,188 patients across seven tumor types treated with pembrolizumab monotherapy in nine clinical trials were used. A total of 11 prespecified gene expression signatures [18-gene T-cell-inflamed gene expression profile (TcellinfGEP), angiogenesis, hypoxia, glycolysis, proliferation, MYC, RAS, granulocytic myeloid-derived suppressor cell (gMDSC), monocytic myeloid-derived suppressor cell (mMDSC), stroma/epithelial-to-mesenchymal transition (EMT)/TGFß, and WNT] were evaluated for their relationship to objective response rate (per RECIST, version 1.1). Logistic regression analysis of response for consensus signatures was adjusted for tumor type, Eastern Cooperative Oncology Group performance status, and TcellinfGEP, an approach equivalent to evaluating the association between response and the residuals of consensus signatures after detrending them for their relationship with the TcellinfGEP (previously identified as a determinant of pembrolizumab response) and tumor type. Testing of the 10 prespecified non-TcellinfGEP consensus signatures for negative association [except proliferation (hypothesized positive association)] with response was adjusted for multiplicity. RESULTS: Covariance patterns of the 11 signatures (including TcellinfGEP) identified in Merck-Moffitt and The Cancer Genome Atlas datasets showed highly concordant coexpression patterns in the RNA-sequencing data from pembrolizumab trials. TcellinfGEP was positively associated with response; signatures for angiogenesis, mMDSC, and stroma/EMT/TGFß were negatively associated with response to pembrolizumab monotherapy. CONCLUSIONS: These findings suggest that features beyond IFNγ-related T-cell inflammation may be relevant to anti-programmed death 1 monotherapy response and may define other axes of tumor biology as candidates for pembrolizumab combinations. See related commentary by Cho et al., p. 1479.

STUB1 is an intracellular checkpoint for interferon gamma sensing.

Immune checkpoint blockade (ICB) leads to durable and complete tumour regression in some patients but in others gives temporary, partial or no response. Accordingly, significant efforts are underway to identify tumour-intrinsic mechanisms underlying ICB resistance. Results from a published CRISPR screen in a mouse model suggested that targeting STUB1, an E3 ligase involved in protein homeostasis, may overcome ICB resistance but the molecular basis of this effect remains unclear. Herein, we report an under-appreciated role of STUB1 to dampen the interferon gamma (IFNγ) response. Genetic deletion of STUB1 increased IFNGR1 abundance on the cell surface and thus enhanced the downstream IFNγ response as showed by multiple approaches including Western blotting, flow cytometry, qPCR, phospho-STAT1 assay, immunopeptidomics, proteomics, and gene expression profiling. Human prostate and breast cancer cells with STUB1 deletion were also susceptible to cytokine-induced growth inhibition. Furthermore, blockade of STUB1 protein function recapitulated the STUB1-null phenotypes. Despite these encouraging in vitro data and positive implications from clinical datasets, we did not observe in vivo benefits of inactivating Stub1 in mouse syngeneic tumour models-with or without combination with anti-PD-1 therapy. However, our findings elucidate STUB1 as a barrier to IFNγ sensing, prompting further investigations to assess if broader inactivation of human STUB1 in both tumors and immune cells could overcome ICB resistance.

Putative Biomarkers of Clinical Benefit With Pembrolizumab in Advanced Urothelial Cancer: Results from the KEYNOTE-045 and KEYNOTE-052 Landmark Trials.

PURPOSE: In an exploratory analysis, we investigated the association between programmed death ligand 1 (PD-L1), tumor mutational burden (TMB), T-cell-inflamed gene expression profile (TcellinfGEP), and stromal signature with outcomes of pembrolizumab in urothelial carcinoma (UC). PATIENTS AND METHODS: Patients with advanced UC received first-line pembrolizumab 200 mg every 3 weeks in the single-arm phase II KEYNOTE-052 trial (NCT02335424) and salvage pembrolizumab 200 mg every 3 weeks or chemotherapy (paclitaxel/docetaxel/vinflunine) in the randomized phase III KEYNOTE-045 trial (NCT02256436). The association of each biomarker (continuous variable) with objective response rate (ORR), progression-free survival (PFS), and overall survival (OS) was evaluated using logistic regression (ORR) and Cox PH (PFS, OS), adjusted for ECOG PS; nominal P values were calculated without multiplicity adjustment (one-sided, pembrolizumab; two-sided, chemotherapy). Significance was prespecified at α = 0.05. RESULTS: In KEYNOTE-052, PD-L1, TMB, and TcellinfGEP were significantly associated with improved outcomes; stromal signature was significantly associated with worse outcomes. In KEYNOTE-045, although findings for TMB and TcellinfGEP with pembrolizumab were consistent with those of KEYNOTE-052, PD-L1 was not significantly associated with improved outcomes, nor was stromal signature associated with worse outcomes with pembrolizumab; chemotherapy was not associated with outcomes in a consistent manner for any of the biomarkers. Hazard ratio (HR) estimates at prespecified cutoffs showed an advantage for pembrolizumab versus chemotherapy regardless of PD-L1 or TMB, with a trend toward lower HRs in the combined positive score ≥10 and the TMB ≥175 mutation/exome subgroup. For TcellinfGEP, PFS and OS HRs were lower in the TcellinfGEP-nonlow subgroup regardless of treatment. CONCLUSIONS: Multiple biomarkers characterizing the tumor microenvironment may help predict response to pembrolizumab monotherapy in UC, and potential clinical utility of these biomarkers may be context-dependent.

Reverse Translating Molecular Determinants of Anti-Programmed Death 1 Immunotherapy Response in Mouse Syngeneic Tumor Models.

Targeting the programmed death 1/programmed death ligand 1 (PD-1/PD-L1) pathway with immunotherapy has revolutionized the treatment of many cancers. Somatic tumor mutational burden (TMB) and T-cell-inflamed gene expression profile (GEP) are clinically validated pan-tumor genomic biomarkers that can predict responsiveness to anti-PD-1/PD-L1 monotherapy in many tumor types. We analyzed the association between these biomarkers and the efficacy of PD-1 inhibitor in 11 commonly used preclinical syngeneic tumor mouse models using murinized rat anti-mouse PD-1 DX400 antibody muDX400, a surrogate for pembrolizumab. Response to muDX400 treatment was broadly classified into three categories: highly responsive, partially responsive, and intrinsically resistant to therapy. Molecular and cellular profiling validated differences in immune cell infiltration and activation in the tumor microenvironment of muDX400-responsive tumors. Baseline and on-treatment genomic analysis showed an association between TMB, murine T-cell-inflamed gene expression profile (murine-GEP), and response to muDX400 treatment. We extended our analysis to investigate a canonical set of cancer and immune biology-related gene signatures, including signatures of angiogenesis, myeloid-derived suppressor cells, and stromal/epithelial-to-mesenchymal transition/TGFß biology previously shown to be inversely associated with the clinical efficacy of immune checkpoint blockade. Finally, we evaluated the association between murine-GEP and preclinical efficacy with standard-of-care chemotherapy or antiangiogenic agents that previously demonstrated promising clinical activity, in combination with muDX400. Our profiling studies begin to elucidate the underlying biological mechanisms of response and resistance to PD-1/PD-L1 blockade represented by these models, thereby providing insight into which models are most appropriate for the evaluation of orthogonal combination strategies.

ILT3 (LILRB4) Promotes the Immunosuppressive Function of Tumor-Educated Human Monocytic Myeloid-Derived Suppressor Cells.

Myeloid-derived suppressor cells (MDSC) are immature myeloid cells that accumulate in the tumor microenvironment (TME). MDSCs have been shown to dampen antitumor immune responses and promote tumor growth; however, the mechanisms of MDSC induction and their role in promoting immune suppression in cancer remain poorly understood. Here, we characterized the phenotype and function of monocytic MDSCs (M-MDSC) generated by coculture of human peripheral blood mononuclear cells with SK-MEL-5 cancer cells in vitro. We selected the SK-MEL-5 human melanoma cell line to generate M-MDSCs because these cells form subcutaneous tumors rich in myeloid cells in humanized mice. M-MDSCs generated via SK-MEL-5 coculture expressed low levels of human leukocyte antigen (HLA)-DR, high levels of CD33 and CD11b, and suppressed both CD8+ T-cell proliferation and IFNγ secretion. M-MDSCs also expressed higher levels of immunoglobulin-like transcript 3 (ILT3, also known as LILRB4) and immunoglobulin-like transcript 4 (ILT4, also known as LILRB2) on the cell surface compared with monocytes. Therefore, we investigated how ILT3 targeting could modulate M-MDSC cell function. Treatment with an anti-ILT3 antibody impaired the acquisition of the M-MDSC suppressor phenotype and reduced the capacity of M-MDSCs to cause T-cell suppression. Finally, in combination with anti-programmed cell death protein 1 (PD1), ILT3 blockade enhanced T-cell activation as assessed by IFNγ secretion. IMPLICATIONS: These results suggest that ILT3 expressed on M-MDSCs has a role in inducing immunosuppression in cancer and that antagonism of ILT3 may be useful to reverse the immunosuppressive function of M-MDSCs and enhance the efficacy of immune checkpoint inhibitors.

Analysis of classical neutrophils and polymorphonuclear myeloid-derived suppressor cells in cancer patients and tumor-bearing mice.

In this study, using single-cell RNA-seq, cell mass spectrometry, flow cytometry, and functional analysis, we characterized the heterogeneity of polymorphonuclear neutrophils (PMNs) in cancer. We describe three populations of PMNs in tumor-bearing mice: classical PMNs, polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs), and activated PMN-MDSCs with potent immune suppressive activity. In spleens of mice, PMN-MDSCs gradually replaced PMNs during tumor progression. Activated PMN-MDSCs were found only in tumors, where they were present at the very early stages of the disease. These populations of PMNs in mice could be separated based on the expression of CD14. In peripheral blood of cancer patients, we identified two distinct populations of PMNs with characteristics of classical PMNs and PMN-MDSCs. The gene signature of tumor PMN-MDSCs was similar to that in mouse activated PMN-MDSCs and was closely associated with negative clinical outcome in cancer patients. Thus, we provide evidence that PMN-MDSCs are a distinct population of PMNs with unique features and potential for selective targeting opportunities.

Combination of EP4 antagonist MF-766 and anti-PD-1 promotes anti-tumor efficacy by modulating both lymphocytes and myeloid cells.

Prostaglandin E2 (PGE2), an arachidonic acid pathway metabolite produced by cyclooxygenase (COX)-1/2, has been shown to impair anti-tumor immunity through engagement with one or more E-type prostanoid receptors (EP1-4). Specific targeting of EP receptors, as opposed to COX-1/2 inhibition, has been proposed to achieve preferential antagonism of PGE2-mediated immune suppression. Here we describe the anti-tumor activity of MF-766, a potent and highly selective small-molecule inhibitor of the EP4 receptor. EP4 inhibition by MF-766 synergistically improved the efficacy of anti-programmed cell death protein 1 (PD-1) therapy in CT26 and EMT6 syngeneic tumor mouse models. Multiparameter flow cytometry analysis revealed that treatment with MF-766 promoted the infiltration of CD8+ T cells, natural killer (NK) cells and conventional dendritic cells (cDCs), induced M1-like macrophage reprogramming, and reduced granulocytic myeloid-derived suppressor cells (MDSC) in the tumor microenvironment (TME). In vitro experiments demonstrated that MF-766 restored PGE2-mediated inhibition of lipopolysaccharide (LPS)-induced tumor necrosis factor (TNF)-α production in THP-1 cells and human blood, and PGE2-mediated inhibition of interleukin (IL)-2-induced interferon (IFN)-γ production in human NK cells. MF-766 reversed the inhibition of IFN-γ in CD8+ T-cells by PGE2 and impaired suppression of CD8+ T-cells induced by myeloid-derived suppressor cells (MDSC)/PGE2. In translational studies using primary human tumors, MF-766 enhanced anti-CD3-stimulated IFN-γ, IL-2, and TNF-α production in primary histoculture and synergized with pembrolizumab in a PGE2 high TME. Our studies demonstrate that the combination of EP4 blockade with anti-PD-1 therapy enhances antitumor activity by differentially modulating myeloid cell, NK cell, cDC and T-cell infiltration profiles.

Classification and prediction of survival in patients with the leukemic phase of cutaneous T cell lymphoma.

We have used cDNA arrays to investigate gene expression patterns in peripheral blood mononuclear cells from patients with leukemic forms of cutaneous T cell lymphoma, primarily Sezary syndrome (SS). When expression data for patients with high blood tumor burden (Sezary cells >60% of the lymphocytes) and healthy controls are compared by Student's t test, at P < 0.01, we find 385 genes to be differentially expressed. Highly overexpressed genes include Th2 cells-specific transcription factors Gata-3 and Jun B, as well as integrin beta1, proteoglycan 2, the RhoB oncogene, and dual specificity phosphatase 1. Highly underexpressed genes include CD26, Stat-4, and the IL-1 receptors. Message for plastin-T, not normally expressed in lymphoid tissue, is detected only in patient samples and may provide a new marker for diagnosis. Using penalized discriminant analysis, we have identified a panel of eight genes that can distinguish SS in patients with as few as 5% circulating tumor cells. This suggests that, even in early disease, Sezary cells produce chemokines and cytokines that induce an expression profile in the peripheral blood distinctive to SS. Finally, we show that using 10 genes, we can identify a class of patients who will succumb within six months of sampling regardless of their tumor burden.

cDNA targets improve whole blood gene expression profiling and enhance detection of pharmocodynamic biomarkers: a quantitative platform analysis.

BACKGROUND: Genome-wide gene expression profiling of whole blood is an attractive method for discovery of biomarkers due to its non-invasiveness, simple clinical site processing and rich biological content. Except for a few successes, this technology has not yet matured enough to reach its full potential of identifying biomarkers useful for clinical prognostic and diagnostic applications or in monitoring patient response to therapeutic intervention. A variety of technical problems have hampered efforts to utilize this technology for identification of biomarkers. One significant hurdle has been the high and variable concentrations of globin transcripts in whole blood total RNA potentially resulting in non-specific probe binding and high background. In this study, we investigated and quantified the power of three whole blood profiling approaches to detect meaningful biological expression patterns. METHODS: To compare and quantify the impact of different mitigation technologies, we used a globin transcript spike-in strategy to synthetically generate a globin-induced signature and then mitigate it with the three different technologies. Biological differences, in globin transcript spiked samples, were modeled by supplementing with either 1% of liver or 1% brain total RNA. In order to demonstrate the biological utility of a robust globin artifact mitigation strategy in biomarker discovery, we treated whole blood ex vivo with suberoylanilide hydroxamic acid (SAHA) and compared the overlap between the obtained signatures and signatures of a known biomarker derived from SAHA-treated cell lines and PBMCs of SAHA-treated patients. RESULTS: We found cDNA hybridization targets detect at least 20 times more specific differentially expressed signatures (2597) between 1% liver and 1% brain in globin-supplemented samples than the PNA (117) or no treatment (97) method at FDR = 10% and p-value < 3x10-3. In addition, we found that the ex vivo derived gene expression profile was highly concordant with that of the previously identified SAHA pharmacodynamic biomarkers. CONCLUSIONS: We conclude that an amplification method for gene expression profiling employing cDNA targets effectively mitigates the negative impact on data of abundant globin transcripts and greatly improves the ability to identify relevant gene expression based pharmacodynamic biomarkers from whole blood.