The somatic genomic landscape of glioblastoma.

We describe the landscape of somatic genomic alterations based on multidimensional and comprehensive characterization of more than 500 glioblastoma tumors (GBMs). We identify several novel mutated genes as well as complex rearrangements of signature receptors, including EGFR and PDGFRA. TERT promoter mutations are shown to correlate with elevated mRNA expression, supporting a role in telomerase reactivation. Correlative analyses confirm that the survival advantage of the proneural subtype is conferred by the G-CIMP phenotype, and MGMT DNA methylation may be a predictive biomarker for treatment response only in classical subtype GBM. Integrative analysis of genomic and proteomic profiles challenges the notion of therapeutic inhibition of a pathway as an alternative to inhibition of the target itself. These data will facilitate the discovery of therapeutic and diagnostic target candidates, the validation of research and clinical observations and the generation of unanticipated hypotheses that can advance our molecular understanding of this lethal cancer.

Genomic and molecular characterization of preterm birth.

Preterm birth (PTB) complications are the leading cause of long-term morbidity and mortality in children. By using whole blood samples, we integrated whole-genome sequencing (WGS), RNA sequencing (RNA-seq), and DNA methylation data for 270 PTB and 521 control families. We analyzed this combined dataset to identify genomic variants associated with PTB and secondary analyses to identify variants associated with very early PTB (VEPTB) as well as other subcategories of disease that may contribute to PTB. We identified differentially expressed genes (DEGs) and methylated genomic loci and performed expression and methylation quantitative trait loci analyses to link genomic variants to these expression and methylation changes. We performed enrichment tests to identify overlaps between new and known PTB candidate gene systems. We identified 160 significant genomic variants associated with PTB-related phenotypes. The most significant variants, DEGs, and differentially methylated loci were associated with VEPTB. Integration of all data types identified a set of 72 candidate biomarker genes for VEPTB, encompassing genes and those previously associated with PTB. Notably, PTB-associated genes RAB31 and RBPJ were identified by all three data types (WGS, RNA-seq, and methylation). Pathways associated with VEPTB include EGFR and prolactin signaling pathways, inflammation- and immunity-related pathways, chemokine signaling, IFN-γ signaling, and Notch1 signaling. Progress in identifying molecular components of a complex disease is aided by integrated analyses of multiple molecular data types and clinical data. With these data, and by stratifying PTB by subphenotype, we have identified associations between VEPTB and the underlying biology.

Multiplex immunofluorescence to measure dynamic changes in tumor-infiltrating lymphocytes and PD-L1 in early-stage breast cancer.

BACKGROUND: The H&E stromal tumor-infiltrating lymphocyte (sTIL) score and programmed death ligand 1 (PD-L1) SP142 immunohistochemistry assay are prognostic and predictive in early-stage breast cancer, but are operator-dependent and may have insufficient precision to characterize dynamic changes in sTILs/PD-L1 in the context of clinical research. We illustrate how multiplex immunofluorescence (mIF) combined with statistical modeling can be used to precisely estimate dynamic changes in sTIL score, PD-L1 expression, and other immune variables from a single paraffin-embedded slide, thus enabling comprehensive characterization of activity of novel immunotherapy agents. METHODS: Serial tissue was obtained from a recent clinical trial evaluating loco-regional cytokine delivery as a strategy to promote immune cell infiltration and activation in breast tumors. Pre-treatment biopsies and post-treatment tumor resections were analyzed by mIF (PerkinElmer Vectra) using an antibody panel that characterized tumor cells (cytokeratin-positive), immune cells (CD3, CD8, CD163, FoxP3), and PD-L1 expression. mIF estimates of sTIL score and PD-L1 expression were compared to the H&E/SP142 clinical assays. Hierarchical linear modeling was utilized to compare pre- and post-treatment immune cell expression, account for correlation of time-dependent measurement, variation across high-powered magnification views within each subject, and variation between subjects. Simulation methods (Monte Carlo, bootstrapping) were used to evaluate the impact of model and tissue sample size on statistical power. RESULTS: mIF estimates of sTIL and PD-L1 expression were strongly correlated with their respective clinical assays (p < .001). Hierarchical linear modeling resulted in more precise estimates of treatment-related increases in sTIL, PD-L1, and other metrics such as CD8+ tumor nest infiltration. Statistical precision was dependent on adequate tissue sampling, with at least 15 high-powered fields recommended per specimen. Compared to conventional t-testing of means, hierarchical linear modeling was associated with substantial reductions in enrollment size required (n = 25➔n = 13) to detect the observed increases in sTIL/PD-L1. CONCLUSION: mIF is useful for quantifying treatment-related dynamic changes in sTILs/PD-L1 and is concordant with clinical assays, but with greater precision. Hierarchical linear modeling can mitigate the effects of intratumoral heterogeneity on immune cell count estimations, allowing for more efficient detection of treatment-related pharmocodynamic effects in the context of clinical trials. TRIAL REGISTRATION: NCT02950259 .

Multiscale mutation clustering algorithm identifies pan-cancer mutational clusters associated with pathway-level changes in gene expression.

Cancer researchers have long recognized that somatic mutations are not uniformly distributed within genes. However, most approaches for identifying cancer mutations focus on either the entire-gene or single amino-acid level. We have bridged these two methodologies with a multiscale mutation clustering algorithm that identifies variable length mutation clusters in cancer genes. We ran our algorithm on 539 genes using the combined mutation data in 23 cancer types from The Cancer Genome Atlas (TCGA) and identified 1295 mutation clusters. The resulting mutation clusters cover a wide range of scales and often overlap with many kinds of protein features including structured domains, phosphorylation sites, and known single nucleotide variants. We statistically associated these multiscale clusters with gene expression and drug response data to illuminate the functional and clinical consequences of mutations in our clusters. Interestingly, we find multiple clusters within individual genes that have differential functional associations: these include PTEN, FUBP1, and CDH1. This methodology has potential implications in identifying protein regions for drug targets, understanding the biological underpinnings of cancer, and personalizing cancer treatments. Toward this end, we have made the mutation clusters and the clustering algorithm available to the public. Clusters and pathway associations can be interactively browsed at m2c.systemsbiology.net. The multiscale mutation clustering algorithm is available at https://github.com/IlyaLab/M2C.

Genomically amplified Akt3 activates DNA repair pathway and promotes glioma progression.

Akt is a robust oncogene that plays key roles in the development and progression of many cancers, including glioma. We evaluated the differential propensities of the Akt isoforms toward progression in the well-characterized RCAS/Ntv-a mouse model of PDGFB-driven low grade glioma. A constitutively active myristoylated form of Akt1 did not induce high-grade glioma (HGG). In stark contrast, Akt2 and Akt3 showed strong progression potential with 78% and 97% of tumors diagnosed as HGG, respectively. We further revealed that significant variations in polarity and hydropathy values among the Akt isoforms in both the pleckstrin homology domain (P domain) and regulatory domain (R domain) were critical in mediating glioma progression. Gene expression profiles from representative Akt-derived tumors indicated dominant and distinct roles for Akt3, consisting primarily of DNA repair pathways. TCGA data from human GBM closely reflected the DNA repair function, as Akt3 was significantly correlated with a 76-gene signature DNA repair panel. Consistently, compared with Akt1 and Akt2 overexpression models, Akt3-expressing human GBM cells had enhanced activation of DNA repair proteins, leading to increased DNA repair and subsequent resistance to radiation and temozolomide. Given the wide range of Akt3-amplified cancers, Akt3 may represent a key resistance factor.

Combining dependent P-values with an empirical adaptation of Brown's method.

MOTIVATION: Combining P-values from multiple statistical tests is a common exercise in bioinformatics. However, this procedure is non-trivial for dependent P-values. Here, we discuss an empirical adaptation of Brown's method (an extension of Fisher's method) for combining dependent P-values which is appropriate for the large and correlated datasets found in high-throughput biology. RESULTS: We show that the Empirical Brown's method (EBM) outperforms Fisher's method as well as alternative approaches for combining dependent P-values using both noisy simulated data and gene expression data from The Cancer Genome Atlas. AVAILABILITY AND IMPLEMENTATION: The Empirical Brown's method is available in Python, R, and MATLAB and can be obtained from https://github.com/IlyaLab/CombiningDependentPvalues UsingEBM The R code is also available as a Bioconductor package from https://www.bioconductor.org/packages/devel/bioc/html/EmpiricalBrownsMethod.html CONTACT: Theo.Knijnenburg@systemsbiology.org SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structure-based predictions broadly link transcription factor mutations to gene expression changes in cancers.

Thousands of unique mutations in transcription factors (TFs) arise in cancers, and the functional and biological roles of relatively few of these have been characterized. Here, we used structure-based methods developed specifically for DNA-binding proteins to systematically predict the consequences of mutations in several TFs that are frequently mutated in cancers. The explicit consideration of protein-DNA interactions was crucial to explain the roles and prevalence of mutations in TP53 and RUNX1 in cancers, and resulted in a higher specificity of detection for known p53-regulated genes among genetic associations between TP53 genotypes and genome-wide expression in The Cancer Genome Atlas, compared to existing methods of mutation assessment. Biophysical predictions also indicated that the relative prevalence of TP53 missense mutations in cancer is proportional to their thermodynamic impacts on protein stability and DNA binding, which is consistent with the selection for the loss of p53 transcriptional function in cancers. Structure and thermodynamics-based predictions of the impacts of missense mutations that focus on specific molecular functions may be increasingly useful for the precise and large-scale inference of aberrant molecular phenotypes in cancer and other complex diseases.

Uterine Inflammatory Myofibroblastic Tumors: Proposed Risk Stratification Model Using Integrated Clinicopathologic and Molecular Analysis.

Inflammatory myofibroblastic tumor (IMT) of the uterus is a rare mesenchymal tumor with largely benign behavior; however, a small subset demonstrate aggressive behavior. While clinicopathologic features have been previously associated with aggressive behavior, these reports are based on small series, and these features are imperfect predictors of clinical behavior. IMTs are most commonly driven by ALK fusions, with additional pathogenic molecular alterations being reported only in rare examples of extrauterine IMTs. In this study, a series of 11 uterine IMTs, 5 of which demonstrated aggressive behavior, were evaluated for clinicopathologic variables and additionally subjected to capture-based next-generation sequencing with or without whole-transcriptome RNA sequencing. In the 6 IMTs without aggressive behavior, ALK fusions were the sole pathogenic alteration. In contrast, all 5 aggressive IMTs harbored pathogenic molecular alterations and numerous copy number changes in addition to ALK fusions, with the majority of the additional alterations present in the primary tumors. We combined our series with cases previously reported in the literature and performed statistical analyses to propose a novel clinicopathologic risk stratification score assigning 1 point each for: age above 45 years, size≥5 cm,≥4 mitotic figures per 10 high-power field, and infiltrative borders. No tumors with 0 points had an aggressive outcome, while 21% of tumors with 1 to 2 points and all tumors with ≥3 points had aggressive outcomes. We propose a 2-step classification model that first uses the clinicopathologic risk stratification score to identify low-risk and high-risk tumors, and recommend molecular testing to further classify intermediate-risk tumors.

Case report: ex vivo tumor organoid drug testing identifies therapeutic options for stage IV ovarian carcinoma.

Patients presenting with stage 4 ovarian carcinoma, including low-grade serous disease, have a poor prognosis. Although platinum-based therapies can offer some response, these therapies are associated with many side effects, and treatment resistance often develops. Toxic side effects along with disease progression render patients unable to receive additional lines of treatment and limit their options to hospice or palliative care. In this case report, we describe a patient with an unusual case of metastatic low-grade serous ovarian cancer with some features of high-grade disease who had received four previous lines of treatment and was suffering from atelectasis, pulmonary embolism, and hydronephrosis. A CLIA-certified drug sensitivity assay of an organoid culture derived from the patient's tumor (PARIS® test) identified several therapeutic options, including the combination of fulvestrant with everolimus. On this treatment regimen, the patient experienced 7 months of stable disease and survived nearly 11 months before succumbing to her disease. This case emphasizes the clinical utility of ex vivo drug testing as a new functional precision medicine approach to identify, in real-time, personalized treatment options for patients, especially those who are not benefiting from standard of care treatments.

Extraordinary clinical response to ibrutinib in low-grade ovarian cancer guided by organoid drug testing.

Low-grade serous ovarian cancer (LGSOC) typically responds poorly to standard platinum-based chemotherapy and new therapeutic approaches are needed. We describe a remarkable response to targeted therapy in a patient with platinum-resistant, advanced LGSOC who had failed standard-of-care chemotherapy and two surgeries. The patient was in rapid decline and entering hospice care on home intravenous (i.v.) opioid analgesics and a malignant bowel obstruction requiring a G-tube. Genomic analysis of the patient's tumor did not indicate obvious therapeutic options. In contrast, a CLIA-certified drug sensitivity assay of an organoid culture derived from the patient's tumor identified several therapeutic choices, including Bruton's tyrosine kinase (BTK) inhibitor ibrutinib, as well as the EGFR inhibitors afatinib and erlotinib. Following off-label administration of daily ibrutinib as monotherapy, the patient had an exceptional clinical turnaround over the following 65 weeks with normalization of CA-125 levels, resolution of the malignant bowel obstruction, halting of pain medications, and improvement of performance status from ECOG 3 to ECOG 1. After 65 weeks of stable disease, the patient's CA-125 levels began to rise, at which point the patient discontinued ibrutinib and began taking afatinib as monotherapy. The patient's CA-125 levels remained stable for an additional 38 weeks but due to anemia and rising CA-125 levels, the patient switched to erlotinib and is currently being monitored. This case highlights the clinical utility of ex vivo drug testing of patient-derived tumor organoids as a new functional precision medicine approach to identify effective personalized therapies for patients who have failed standard-of-care treatments.

A phase Ib trial of pembrolizumab plus paclitaxel or flat-dose capecitabine in 1st/2nd line metastatic triple-negative breast cancer.

Chemoimmunotherapy with anti-programmed cell death 1/ligand 1 and cytotoxic chemotherapy is a promising therapeutic modality for women with triple-negative breast cancer, but questions remain regarding optimal chemotherapy backbone and biomarkers for patient selection. We report final outcomes from a phase Ib trial evaluating pembrolizumab (200 mg IV every 3 weeks) with either weekly paclitaxel (80 mg/m2 weekly) or flat-dose capecitabine (2000 mg orally twice daily for 7 days of every 14-day cycle) in the 1st/2nd line setting. The primary endpoint is safety (receipt of 2 cycles without grade III/IV toxicities requiring discontinuation or ≥21-day delays). The secondary endpoint is efficacy (week 12 objective response). Exploratory aims are to characterize immunologic effects of treatment over time, and to evaluate novel biomarkers. The trial demonstrates that both regimens meet the pre-specified safety endpoint (paclitaxel: 87%; capecitabine: 100%). Objective response rate is 29% for pembrolizumab/paclitaxel (n = 4/13, 95% CI: 10-61%) and 43% for pembrolizumab/capecitabine (n = 6/14, 95% CI: 18-71%). Partial responses are observed in two subjects with chemo-refractory metaplastic carcinoma (both in capecitabine arm). Both regimens are associated with significant peripheral leukocyte contraction over time. Response is associated with clinical PD-L1 score, non-receipt of prior chemotherapy, and the H&E stromal tumor-infiltrating lymphocyte score, but also by a novel 27 gene IO score and spatial biomarkers (lymphocyte spatial skewness). In conclusion, pembrolizumab with paclitaxel or capecitabine is safe and clinically active. Both regimens are lymphodepleting, highlighting the competing immunostimulatory versus lymphotoxic effects of cytotoxic chemotherapy. Further exploration of the IO score and spatial TIL biomarkers is warranted. The clinical trial registration is NCT02734290.

PD-1 and ICOS coexpression identifies tumor-reactive CD4+ T cells in human solid tumors.

CD4+ Th cells play a key role in orchestrating immune responses, but the identity of the CD4+ Th cells involved in the antitumor immune response remains to be defined. We analyzed the immune cell infiltrates of head and neck squamous cell carcinoma and colorectal cancers and identified a subset of CD4+ Th cells distinct from FOXP3+ Tregs that coexpressed programmed cell death 1 (PD-1) and ICOS. These tumor-infiltrating lymphocyte CD4+ Th cells (CD4+ Th TILs) had a tissue-resident memory phenotype, were present in MHC class II-rich areas, and proliferated in the tumor, suggesting local antigen recognition. The T cell receptor repertoire of the PD-1+ICOS+ CD4+ Th TILs was oligoclonal, with T cell clones expanded in the tumor, but present at low frequencies in the periphery. Finally, these PD-1+ICOS+ CD4+ Th TILs were shown to recognize both tumor-associated antigens and tumor-specific neoantigens. Our findings provide an approach for isolating tumor-reactive CD4+ Th TILs directly ex vivo that will help define their role in the antitumor immune response and potentially improve future adoptive T cell therapy approaches.

Changes in T-cell subsets and clonal repertoire during chemoimmunotherapy with pembrolizumab and paclitaxel or capecitabine for metastatic triple-negative breast cancer.

BACKGROUND: Chemoimmunotherapy is a standard treatment for triple-negative breast cancer (TNBC), however, the impacts of different chemotherapies on T-cell populations, which could correlate with clinical activity, are not known. Quantifying T-cell populations with flow cytometry and T-cell receptor (TCR) immunosequencing may improve our understanding of how chemoimmunotherapy affects T-cell subsets, and to what extent clonal shifts occur during treatment. TCR immunosequencing of intratumoral T cells may facilitate the identification and monitoring of putatively tumor-reactive T-cell clones within the blood. METHODS: Blood and tumor biopsies were collected from patients with metastatic TNBC enrolled in a phase Ib clinical trial of first or second-line pembrolizumab with paclitaxel or capecitabine. Using identical biospecimen processing protocols, blood samples from a cohort of patients treated for early-stage breast cancer were obtained for comparison. Treatment-related immunological changes in peripheral blood and intratumoral T cells were characterized using flow cytometry and TCR immunosequencing. Clonal proliferation rates of T cells were compared based on intratumoral enrichment. RESULTS: When combined with pembrolizumab, paclitaxel and capecitabine resulted in similar time-dependent lymphodepletions across measured peripheral T-cell subsets. Their effects were more modest than that observed following curative-intent dose-dense anthracycline and cyclophosphamide (ddAC) (average fold-change in CD3+ cells, capecitabine: -0.42, paclitaxel: -0.56, ddAC: -1.21). No differences in T-cell clonality or richness were observed following capecitabine or paclitaxel-based treatments. Regression modeling identified differences in the emergence of novel T-cell clones that were not detected at baseline (odds compared with ddAC, capecitabine: 0.292, paclitaxel: 0.652). Pembrolizumab with paclitaxel or capecitabine expanded T-cell clones within tumors; however, these clones did not always expand within the blood. Proliferation rates within the blood were similar between clones that were enriched and those that were not enriched within tumors. CONCLUSION: Chemoimmunotherapy for metastatic TNBC with pembrolizumab and capecitabine or paclitaxel resulted in similar peripheral T-cell subset lymphodepletion without altering T-cell clonal diversity. Regression modeling methods are applicable in immune monitoring studies, such as this to identify the odds of novel T-cell clones emerging during treatment, and proliferation rates of tumor-enriched T-cell clones.

Transcriptomic profiles of neoantigen-reactive T cells in human gastrointestinal cancers.

Tumor-infiltrating neoantigen-reactive T cells can mediate regression of metastatic gastrointestinal cancers yet remain poorly characterized. We performed immunological screening against personalized neoantigens in combination with single-cell RNA sequencing on tumor-infiltrating lymphocytes from bile duct and pancreatic cancer patients to characterize the transcriptomic landscape of neoantigen-reactive T cells. We found that most neoantigen-reactive CD8+ T cells displayed an exhausted state with significant CXCL13 and GZMA co-expression compared with non-neoantigen-reactive bystander cells. Most neoantigen-reactive CD4+ T cells from a patient with bile duct cancer also exhibited an exhausted phenotype but with overexpression of HOPX or ADGRG1 while lacking IL7R expression. Thus, neoantigen-reactive T cells infiltrating gastrointestinal cancers harbor distinct transcriptomic signatures, which may provide new opportunities for harnessing these cells for therapy.

AACR Project GENIE: 100,000 Cases and Beyond.

The American Association for Cancer Research (AACR) Project Genomics Evidence Neoplasia Information Exchange (GENIE) is an international pan-cancer registry with the goal to inform cancer research and clinical care worldwide. Founded in late 2015, the milestone GENIE 9.1-public release contains data from >110,000 tumors from >100,000 people treated at 19 cancer centers from the United States, Canada, the United Kingdom, France, the Netherlands, and Spain. Here, we demonstrate the use of these real-world data, harmonized through a centralized data resource, to accurately predict enrollment on genome-guided trials, discover driver alterations in rare tumors, and identify cancer types without actionable mutations that could benefit from comprehensive genomic analysis. The extensible data infrastructure and governance framework support additional deep patient phenotyping through biopharmaceutical collaborations and expansion to include new data types such as cell-free DNA sequencing. AACR Project GENIE continues to serve a global precision medicine knowledge base of increasing impact to inform clinical decision-making and bring together cancer researchers internationally. SIGNIFICANCE: AACR Project GENIE has now accrued data from >110,000 tumors, placing it among the largest repository of publicly available, clinically annotated genomic data in the world. GENIE has emerged as a powerful resource to evaluate genome-guided clinical trial design, uncover drivers of cancer subtypes, and inform real-world use of genomic data. This article is highlighted in the In This Issue feature, p. 2007.

Developing an enhanced 7-color multiplex IHC protocol to dissect immune infiltration in human cancers.

The TSA Opal multiplex immunohistochemistry (mIHC) protocol (PerkinElmer) has been used to characterize immune infiltration in human cancers. This technique allows multiple biomarkers to be simultaneously stained in a single tissue section, which helps to elucidate the spatial relationship among individual cell types. We developed and optimized two improved mIHC protocols for a 7-color panel containing 6 biomarkers (CD3, CD8, CD163, PD-L1, FoxP3, and cytokeratin (CK)) and DAPI. The only difference between these two protocols was the staining sequence of those 6 biomarkers as the first sequence is PD-L1/CD163/CD8/CK/CD3/FoxP3/DAPI and the second sequence is FoxP3/CD163/CD8/CK/CD3/PD-L1/DAPI. By comparing PD-L1/FoxP3 staining in mIHC and singleplex PD-L1/FoxP3 staining on the adjacent slide, we demonstrated that the staining sequence does not affect the staining intensity of individual biomarkers as long as a proper antigen retrieval method was used. Our study suggests that use of an antigen retrieval buffer with higher pH value (such as Tris-EDTA pH9.0) than that of the stripping buffers (such as citrate buffer pH6.0) is helpful when using this advanced mIHC method to develop panels with multiple biomarkers. Otherwise, individual biomarkers may exhibit different intensities when the staining sequence is changed. By using this protocol, we characterized immune infiltration and PD-L1 expression in head and neck squamous cell carcinoma (HNSCC), breast cancer (BCa), and non-small cell lung cancer (NSCLC) specimens. We observed a statistically significant increase in CD3+ cell populations within the stroma of NSCLC as compared to BCa and increased PD-L1+ tumor cells in HNSCC as opposed to BCa.

Robust Antitumor Immunity in a Patient with Metastatic Colorectal Cancer Treated with Cytotoxic Regimens.

Microsatellite-stable (MSS) colorectal cancers are characterized by low mutation burden and limited immune-cell infiltration and thereby respond poorly to immunotherapy. Here, we report a case of metastatic MSS colorectal cancer with a robust anticancer immune response. The primary tumor was resected in 2012, and the patient received several cycles of chemotherapy until 2017. In 2018, the patient underwent a left hepatectomy to remove a new metastasis. Analysis of the metastatic tumor revealed a strong CD8+ T-cell response. A high frequency of CD8+ T cells coexpressed CD39 and CD103, a phenotype characteristic of tumor-reactive cells. Using whole-exome sequencing, we identified somatic mutations that generated peptides recognized by CD39+CD103+CD8+ T cells. The observed reactivity against the tumor was dominated by the response to a single mutation that emerged in the metastasis. Somatic mutations that were not immunogenic in the primary tumor led to robust CD8+ T-cell expansion later during disease progression. Our data suggest that the cytotoxic treatment regimen received by the patient might be responsible for this effect. Hence, the capacity of cytotoxic regimens to prime the immune system in colorectal cancer patients should be investigated further and might provide a rationale for combination with immunotherapy.

A Case Series of Metastatic Metaplastic Breast Carcinoma Treated With Anti-PD-1 Therapy.

Metaplastic breast cancer is a rare and often chemo-refractory subtype of breast cancer with poor prognosis and limited treatment options. Recent studies have reported overexpression of programmed death ligand 1 (PD-L1) in metaplastic breast cancers, and there are several reports of anti-PD-1/L1 being potentially active in this disease. In this case series, we present 5 patients with metastatic metaplastic breast cancer treated with anti-PD-1-based therapy at a single center, with 3 of 5 cases demonstrating a response to therapy, and one of the responding cases being a metaplastic lobular carcinoma with low-level hormone receptor expression. Cases were evaluated for PD-L1 expression, tumor infiltrating lymphocytes (TILs), DNA mutations, RNA sequencing, and T-cell receptor sequencing. Duration of the response in these cases was limited, in contrast to the more durable responses noted in other recently published reports.

Neoadjuvant anti-OX40 (MEDI6469) therapy in patients with head and neck squamous cell carcinoma activates and expands antigen-specific tumor-infiltrating T cells.

Despite the success of checkpoint blockade in some cancer patients, there is an unmet need to improve outcomes. Targeting alternative pathways, such as costimulatory molecules (e.g. OX40, GITR, and 4-1BB), can enhance T cell immunity in tumor-bearing hosts. Here we describe the results from a phase Ib clinical trial (NCT02274155) in which 17 patients with locally advanced head and neck squamous cell carcinoma (HNSCC) received a murine anti-human OX40 agonist antibody (MEDI6469) prior to definitive surgical resection. The primary endpoint was to determine safety and feasibility of the anti-OX40 neoadjuvant treatment. The secondary objective was to assess the effect of anti-OX40 on lymphocyte subsets in the tumor and blood. Neoadjuvant anti-OX40 was well tolerated and did not delay surgery, thus meeting the primary endpoint. Peripheral blood phenotyping data show increases in CD4+ and CD8+ T cell proliferation two weeks after anti-OX40 administration. Comparison of tumor biopsies before and after treatment reveals an increase of activated, conventional CD4+ tumor-infiltrating lymphocytes (TIL) in most patients and higher clonality by TCRß sequencing. Analyses of CD8+ TIL show increases in tumor-antigen reactive, proliferating CD103+ CD39+ cells in 25% of patients with evaluable tumor tissue (N = 4/16), all of whom remain disease-free. These data provide evidence that anti-OX40 prior to surgery is safe and can increase activation and proliferation of CD4+ and CD8+ T cells in blood and tumor. Our work suggests that increases in the tumor-reactive CD103+ CD39+ CD8+ TIL could serve as a potential biomarker of anti-OX40 clinical activity.

Germinal center reactions in tertiary lymphoid structures associate with neoantigen burden, humoral immunity and long-term survivorship in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) has traditionally been thought of as an immunologically quiescent tumor type presumably because of a relatively low tumor mutational burden (TMB) and poor responses to checkpoint blockade therapy. However, many PDAC tumors exhibit T cell inflamed phenotypes. The presence of tertiary lymphoid structures (TLS) has recently been shown to be predictive of checkpoint blockade response in melanomas and sarcomas, and are prognostic for survival in PDAC. In order to more comprehensively understand tumor immunity in PDAC patients with TLS, we performed RNA-seq, single and multiplex IHC, flow cytometry and predictive genomic analysis on treatment naïve, PDAC surgical specimens. Forty-six percent of tumors contained distinct T and B cell aggregates reflective of "early-stage TLS" (ES-TLS), which correlated with longer overall and progression-free survival. These tumors had greater CD8+ T cell infiltration but were not defined by previously published TLS gene-expression signatures. ES-TLS+ tumors were enriched for IgG1 class-switched memory B cells and memory CD4+ T cells, suggesting durable immunological memory persisted in these patients. We also observed the presence of active germinal centers (mature-TLS) in 31% of tumors with lymphocyte clusters, whose patients had long-term survival (median 56 months). M-TLS-positive tumors had equivalent overall T cell infiltration to ES-TLS, but were enriched for activated CD4+ memory cells, naive B cells and NK cells. Finally, using a TCGA-PDAC dataset, ES-TLS+ tumors harbored a decreased TMB, but M-TLS with germinal centers expressed significantly more MHCI-restricted neoantigens as determined by an in silico neoantigen prediction method. Interestingly, M-TLS+ tumors also had evidence of increased rates of B cell somatic hypermutation, suggesting that germinal centers form in the presence of high-quality tumor neoantigens leading to increased humoral immunity that confers improved survival for PDAC patients. AbbreviationsTLS: tertiary lymphoid structures; GC: germinal center(s); PDAC: pancreatic ductal adenocarcinoma; RNA-seq: RNA sequencing; BCRseq: B cell receptor sequencing; HEV: high endothelial venule; PNAd: peripheral node addressin; TMB: tumor mutational burden; TCGA: the cancer genome atlas; PAAD: pancreatic adenocarcinoma; FFPE: formalin fixed paraffin embedded; TIME: tumor immune microenvironment.