Human lung cancer harbors spatially organized stem-immunity hubs associated with response to immunotherapy

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens and found an association with beneficial response to PD-1 blockade. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcome. This hub is distinct from mature tertiary lymphoid structures and is enriched for stem-like TCF7+PD-1+CD8+ T cells, activated CCR7+LAMP3+ dendritic cells and CCL19+ fibroblasts as well as chemokines that organize these cells. Within the stem-immunity hub, we find preferential interactions between CXCL10+ macrophages and TCF7-CD8+ T cells as well as between mature regulatory dendritic cells and TCF7+CD4+ and regulatory T cells. These results provide a picture of the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Identification of transcriptional programs using dense vector representations defined by mutual information with GeneVector

Deciphering individual cell phenotypes from cell-specific transcriptional processes requires high dimensional single cell RNA sequencing. However, current dimensionality reduction methods aggregate sparse gene information across cells, without directly measuring the relationships that exist between genes. By performing dimensionality reduction with respect to gene co-expression, low-dimensional features can model these gene-specific relationships and leverage shared signal to overcome sparsity. We describe GeneVector, a scalable framework for dimensionality reduction implemented as a vector space model using mutual information between gene expression. Unlike other methods, including principal component analysis and variational autoencoders, GeneVector uses latent space arithmetic in a lower dimensional gene embedding to identify transcriptional programs and classify cell types. In this work, we show in four single cell RNA-seq datasets that GeneVector was able to capture phenotype-specific pathways, perform batch effect correction, interactively annotate cell types, and identify pathway variation with treatment over time.

Genomic and transcriptomic analysis of checkpoint blockade response in advanced non-small cell lung cancer

Anti-PD-1/PD-L1 agents have transformed the treatment landscape of advanced non-small cell lung cancer (NSCLC). To expand our understanding of the molecular features underlying response to checkpoint inhibitors in NSCLC, we describe here the first joint analysis of the Stand Up To Cancer-Mark Foundation cohort, a resource of whole exome and/or RNA sequencing from 393 patients with NSCLC treated with anti-PD-(L)1 therapy, along with matched clinical response annotation. We identify a number of associations between molecular features and outcome, including (1) favorable (for example, ATM altered) and unfavorable (for example, TERT amplified) genomic subgroups, (2) a prominent association between expression of inducible components of the immunoproteasome and response and (3) a dedifferentiated tumor-intrinsic subtype with enhanced response to checkpoint blockade. Taken together, results from this cohort demonstrate the complexity of biological determinants underlying immunotherapy outcomes and reinforce the discovery potential of integrative analysis within large, well-curated, cancer-specific cohorts.

Spatial analysis of human lung cancer reveals organized immune hubs enriched for stem-like CD8 T cells and associated with immunotherapy response

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially-localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens, and found that they were associated with beneficial responses to PD-1-blockade. Immunity hubs were enriched for many interferon-stimulated genes, T cells in multiple differentiation states, and CXCL9/10/11 + macrophages that preferentially interact with CD8 T cells. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcomes, distinct from mature tertiary lymphoid structures, and enriched for stem-like TCF7+PD-1+ CD8 T cells and activated CCR7 + LAMP3 + dendritic cells, as well as chemokines that organize these cells. These results elucidate the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Protocol for bulk RNA sequencing of enriched human neutrophils from whole blood and estimation of sample purity

Although neutrophils are the most abundant leukocyte in healthy individuals and impact outcomes of diseases ranging from sepsis to cancer, they remain understudied due to technical constraints of isolation, preservation, and sequencing. We present a modified Smart-Seq2 protocol for bulk RNA sequencing of neutrophils enriched from whole blood. We describe steps for neutrophil isolation, cDNA generation, library preparation, and sample purity estimation via a bioinformatic approach. Our approach permits the collection of large cohorts and enables detection of neutrophil transcriptomic subtypes. For complete details on the use and execution of this protocol, please refer to LaSalle et al. (2022)1 and Boribong et al. (2022).2.

Personalised progression prediction in patients with monoclonal gammopathy of undetermined significance or smouldering multiple myeloma (PANGEA): a retrospective, multicohort study

BACKGROUND

Patients with precursors to multiple myeloma are dichotomised as having monoclonal gammopathy of undetermined significance or smouldering multiple myeloma on the basis of monoclonal protein concentrations or bone marrow plasma cell percentage. Current risk stratifications use laboratory measurements at diagnosis and do not incorporate time-varying biomarkers. Our goal was to develop a monoclonal gammopathy of undetermined significance and smouldering multiple myeloma stratification algorithm that utilised accessible, time-varying biomarkers to model risk of progression to multiple myeloma.

METHODS

In this retrospective, multicohort study, we included patients who were 18 years or older with monoclonal gammopathy of undetermined significance or smouldering multiple myeloma. We evaluated several modelling approaches for predicting disease progression to multiple myeloma using a training cohort (with patients at Dana-Farber Cancer Institute, Boston, MA, USA; annotated from Nov, 13, 2019, to April, 13, 2022). We created the PANGEA models, which used data on biomarkers (monoclonal protein concentration, free light chain ratio, age, creatinine concentration, and bone marrow plasma cell percentage) and haemoglobin trajectories from medical records to predict progression from precursor disease to multiple myeloma. The models were validated in two independent validation cohorts from National and Kapodistrian University of Athens (Athens, Greece; from Jan 26, 2020, to Feb 7, 2022; validation cohort 1), University College London (London, UK; from June 9, 2020, to April 10, 2022; validation cohort 1), and Registry of Monoclonal Gammopathies (Czech Republic, Czech Republic; Jan 5, 2004, to March 10, 2022; validation cohort 2). We compared the PANGEA models (with bone marrow [BM] data and without bone marrow [no BM] data) to current criteria (International Myeloma Working Group [IMWG] monoclonal gammopathy of undetermined significance and 20/2/20 smouldering multiple myeloma risk criteria).

FINDINGS

We included 6441 patients, 4931 (77%) with monoclonal gammopathy of undetermined significance and 1510 (23%) with smouldering multiple myeloma. 3430 (53%) of 6441 participants were female. The PANGEA model (BM) improved prediction of progression from smouldering multiple myeloma to multiple myeloma compared with the 20/2/20 model, with a C-statistic increase from 0·533 (0·480-0·709) to 0·756 (0·629-0·785) at patient visit 1 to the clinic, 0·613 (0·504-0·704) to 0·720 (0·592-0·775) at visit 2, and 0·637 (0·386-0·841) to 0·756 (0·547-0·830) at visit three in validation cohort 1. The PANGEA model (no BM) improved prediction of smouldering multiple myeloma progression to multiple myeloma compared with the 20/2/20 model with a C-statistic increase from 0·534 (0·501-0·672) to 0·692 (0·614-0·736) at visit 1, 0·573 (0·518-0·647) to 0·693 (0·605-0·734) at visit 2, and 0·560 (0·497-0·645) to 0·692 (0·570-0·708) at visit 3 in validation cohort 1. The PANGEA models improved prediction of monoclonal gammopathy of undetermined significance progression to multiple myeloma compared with the IMWG rolling model at visit 1 in validation cohort 2, with C-statistics increases from 0·640 (0·518-0·718) to 0·729 (0·643-0·941) for the PANGEA model (BM) and 0·670 (0·523-0·729) to 0·879 (0·586-0·938) for the PANGEA model (no BM).

INTERPRETATION

Use of the PANGEA models in clinical practice will allow patients with precursor disease to receive more accurate measures of their risk of progression to multiple myeloma, thus prompting for more appropriate treatment strategies.

FUNDING

SU2C Dream Team and Cancer Research UK.

Longitudinal characterization of circulating neutrophils uncovers phenotypes associated with severity in hospitalized COVID-19 patients

Mechanisms of neutrophil involvement in severe coronavirus disease 2019 (COVID-19) remain incompletely understood. Here, we collect longitudinal blood samples from 306 hospitalized COVID-19+ patients and 86 controls and perform bulk RNA sequencing of enriched neutrophils, plasma proteomics, and high-throughput antibody profiling to investigate relationships between neutrophil states and disease severity. We identify dynamic switches between six distinct neutrophil subtypes. At days 3 and 7 post-hospitalization, patients with severe disease display a granulocytic myeloid-derived suppressor cell-like gene expression signature, while patients with resolving disease show a neutrophil progenitor-like signature. Humoral responses are identified as potential drivers of neutrophil effector functions, with elevated severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific immunoglobulin G1 (IgG1)-to-IgA1 ratios in plasma of severe patients who survived. In vitro experiments confirm that while patient-derived IgG antibodies induce phagocytosis in healthy donor neutrophils, IgA antibodies predominantly induce neutrophil cell death. Overall, our study demonstrates a dysregulated myelopoietic response in severe COVID-19 and a potential role for IgA-dominant responses contributing to mortality.

In vivo CRISPR screens reveal the landscape of immune evasion pathways across cancer

The immune system can eliminate tumors, but checkpoints enable immune escape. Here, we identify immune evasion mechanisms using genome-scale in vivo CRISPR screens across cancer models treated with immune checkpoint blockade (ICB). We identify immune evasion genes and important immune inhibitory checkpoints conserved across cancers, including the non-classical major histocompatibility complex class I (MHC class I) molecule Qa-1b/HLA-E. Surprisingly, loss of tumor interferon-γ (IFNγ) signaling sensitizes many models to immunity. The immune inhibitory effects of tumor IFN sensing are mediated through two mechanisms. First, tumor upregulation of classical MHC class I inhibits natural killer cells. Second, IFN-induced expression of Qa-1b inhibits CD8+ T cells via the NKG2A/CD94 receptor, which is induced by ICB. Finally, we show that strong IFN signatures are associated with poor response to ICB in individuals with renal cell carcinoma or melanoma. This study reveals that IFN-mediated upregulation of classical and non-classical MHC class I inhibitory checkpoints can facilitate immune escape.

Abstract 2259: Predictive modeling of smoldering multiple myeloma progression to multiple myeloma by continuous variable analysis

Introduction: Multiple myeloma (MM) is consistently preceded by two precursor conditions, monoclonal gammopathy of undetermined significance (MGUS) and smoldering multiple myeloma (SMM). The distinctions between MGUS, SMM, and MM rely on clinical values with inherent variation relating to tumor burden quantified from bone marrow biopsies, a measure subject to inconsistencies in location, timing, and pathologist interpretation. These challenges limit the potential of current standardized risk criteria and advocate for models that examine precursor disease kinetics. We thus developed a risk model that leverages dynamic changes in markers of precursor disease to improve clinical ability to predict time to disease progression.

Methods: To model the evolution of progression risk, we built PANGEA, an international retrospective cohort of precursor patients with baseline and serial time points of clinical and biological variables. This cohort comprises 1095 SMM patients, 254 (23%) of which progressed to MM. Using this cohort, we modeled progression to MM with Cox regression using time-dependent and continuous clinical variables. The model was trained on a subset of data restricted to patients of the Dana-Farber Cancer Institute (DFCI) and validated its performance by computing the c-statistic in a sub-cohort independent from the DFCI training cohort.

Results: The PANGEA cohort was first used to validate current models of SMM disease progression. We validated the 20/2/20 International Myeloma Working Group criteria for SMM patients using binary cutoffs of initial measurements (baseline model), and then extended this model, allowing for re-stratification by the 20/2/20 criteria over time (dynamic model). We then assessed whether rates of change in a set of myeloma-specific clinical variables unrestricted to those of the 20/2/20 criteria improved the predictive ability of the model. This improved our progression prediction as indicated by a c-statistic increase of more than 10% with respect to both 20/2/20 models (baseline and dynamic). Specifically, changes in disease indicators such as age and creatinine are highly predictive of imminent disease progression (p-value < 0.01). Finally, we clustered patients based on latent trajectories of these time-varying clinical variables and included the trajectory classes in the Cox regression. The resulting multivariable, dynamic algorithm is a dramatic improvement over current clinical standards in predicting progression from SMM to MM disease.

Conclusion: The PANGEA multivariable algorithm’s use of continuous clinical variables enhances progression risk predictions in SMM. These findings demonstrate that disease progression from SMM to MM, which likely occurs by the acquisition of sequential changes to the plasma cell clone, can be tracked by trends in clinical values, thus improving prognostication for precursor patients.

Citation Format: Annie Cowan, Habib El-Khoury, Federico Ferrari, Samuel S. Freeman, Robert Redd, Jacqueline Perry, Vidhi Patel, Priya Kaur, Hadley Barr, Katelyn Downey, David Argyelan, Anna V. Justis, David J. Lee, Elizabeth D. Lightbody, Foteini Theodorakakou, Despina Fotiou, Nikolaos Kanellias, Christine Liacos, Gad Getz, Lorenzo Trippa, Catherine Marinac, Efstathios Kastritis, Dimopoulos Meletios, Irene Ghobrial. Predictive modeling of smoldering multiple myeloma progression to multiple myeloma by continuous variable analysis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2259.

Abstract 3610: In vivo CRISPR screens reveal the landscape of immune evasion pathways across cancer

The immune system can eliminate tumors, but checkpoints enable tumors to escape immune destruction. Here, we report the systematic identification of immune evasion mechanisms using genome-scale in vivo CRISPR screens in eight murine cancer models treated with immune checkpoint blockade (ICB). We identify and validate previously unreported immune evasion genes and identify key immune inhibitory checkpoints that have a conserved role across several cancer models, such as the non-classical MHC-I molecule Qa-1b/HLA-E, which scores as the top overall sensitizing hit across all screens. Surprisingly, we find that loss of IFNγ signaling by tumor cells sensitizes 6 of 8 cancer models to ICB. While IFN-mediated inflammation has been associated with response to ICB, there have also been reports of ICB-resistance driven by IFN sensing. However, several divergent mechanisms have been proposed to explain the inhibitory effect of tumor IFN sensing, leading to uncertainty about how this key immune signaling pathway is regulating anti-tumor immunity in different contexts. Using in vivo screening data, transcriptional profiling, and genetic interaction studies, we reveal that the immune-inhibitory effects of tumor IFN sensing are the direct result of tumor upregulation of classical and non-classical MHC-I genes. The interferon-MHC-I axis can inhibit anti-tumor immunity through two mechanisms: first, upregulation of classical MHC-I inhibits the cytotoxicity of natural killer cells, which are activated by ICB. Second, IFN-mediated upregulation of Qa-1b directly inhibits cytotoxicity by effector CD8+ T cells via the NKG2A/CD94 receptor, which is induced on CD8+ T cells by ICB. Finally, we show that high interferon-stimulated gene expression in patients is associated with decreased survival in RCC and poor response to ICB in melanoma. Our study establishes a unifying mechanism to explain the inhibitory role of tumor IFN sensing, revealing that IFN-mediated upregulation of classical and non-classical MHC-I inhibitory checkpoints can facilitate immune escape.

Citation Format: Juan Dubrot, Peter P. Du, Sarah Kate Lane-Reticker, Emily A. Kessler, Audrey J. Muscato, Arnav Mehta, Samuel S. Freeman, Peter M. Allen, Kira E. Olander, Kyle M. Ockerman, Clara H. Wolfe, Fabius Wiesmann, Nelson H. Knudsen, Hsiao-Wei Tsao, Arvin Iracheta-Vellve, Emily M. Schneider, Andrea N. Rivera-Rosario, Ian C. Kohnle, Hans W. Pope, Austin Ayer, Gargi Mishra, Margaret D. Zimmer, Sarah Y. Kim, Animesh Mahapatra, Hakimeh Ebrahimi-Nik, Dennie T. Frederick, Genevieve M. Boland, W. Nicholas Haining, David E. Root, John G. Doench, Nir Hacohen, Kathleen B. Yates, Robert T. Manguso. In vivo CRISPR screens reveal the landscape of immune evasion pathways across cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3610.

Landscape of helper and regulatory antitumour CD4+ T cells in melanoma

Within the tumour microenvironment, CD4+ T cells can promote or suppress antitumour responses through the recognition of antigens presented by human leukocyte antigen (HLA) class II molecules1,2, but how cancers co-opt these physiologic processes to achieve immune evasion remains incompletely understood. Here we performed in-depth analysis of the phenotype and tumour specificity of CD4+ T cells infiltrating human melanoma specimens, finding that exhausted cytotoxic CD4+ T cells could be directly induced by melanoma cells through recognition of HLA class II-restricted neoantigens, and also HLA class I-restricted tumour-associated antigens. CD4+ T regulatory (TReg) cells could be indirectly elicited through presentation of tumour antigens via antigen-presenting cells. Notably, numerous tumour-reactive CD4+ TReg clones were stimulated directly by HLA class II-positive melanoma and demonstrated specificity for melanoma neoantigens. This phenomenon was observed in the presence of an extremely high tumour neoantigen load, which we confirmed to be associated with HLA class II positivity through the analysis of 116 melanoma specimens. Our data reveal the landscape of infiltrating CD4+ T cells in melanoma and point to the presentation of HLA class II-restricted neoantigens and direct engagement of immunosuppressive CD4+ TReg cells as a mechanism of immune evasion that is favoured in HLA class II-positive melanoma.

Implementation of a prostate cancer-specific targeted sequencing panel for credentialing of patient-derived cell lines and genomic characterization of patient samples

BACKGROUND

Primary and metastatic prostate cancers have low mutation rates and recurrent alterations in a small set of genes, enabling targeted sequencing of prostate cancer-associated genes as an efficient approach to characterizing patient samples (compared to whole-exome and whole-genome sequencing). For example, targeted sequencing provides a flexible, rapid, and cost-effective method for genomic assessment of patient-derived cell lines to evaluate fidelity to initial patient tumor samples.

METHODS

We developed a prostate cancer-specific targeted next-generation sequencing (NGS) panel to detect alterations in 62 prostate cancer-associated genes as well as recurring gene fusions with ETS family members, representing the majority of common alterations in prostate cancer. We tested this panel on primary prostate cancer tissues and blood biopsies from patients with metastatic prostate cancer. We generated patient-derived cell lines from primary prostate cancers using conditional reprogramming methods and applied targeted sequencing to evaluate the fidelity of these cell lines to the original patient tumors.

RESULTS

The prostate cancer-specific panel identified biologically and clinically relevant alterations, including point mutations in driver oncogenes and ETS family fusion genes, in tumor tissues from 29 radical prostatectomy samples. The targeted panel also identified genomic alterations in cell-free DNA and circulating tumor cells (CTCs) from patients with metastatic prostate cancer, and in standard prostate cancer cell lines. We used the targeted panel to sequence our set of patient-derived cell lines; however, no prostate cancer-specific mutations were identified in the tumor-derived cell lines, suggesting preferential outgrowth of normal prostate epithelial cells.

CONCLUSIONS

We evaluated a prostate cancer-specific targeted NGS panel to detect common and clinically relevant alterations (including ETS family gene fusions) in prostate cancer. The panel detected driver mutations in a diverse set of clinical samples of prostate cancer, including fresh-frozen tumors, cell-free DNA, CTCs, and cell lines. Targeted sequencing of patient-derived cell lines highlights the challenge of deriving cell lines from primary prostate cancers and the importance of genomic characterization to credential candidate cell lines. Our study supports that a prostate cancer-specific targeted sequencing panel provides an efficient, clinically feasible approach to identify genetic alterations across a spectrum of prostate cancer samples and cell lines.

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma

Immune checkpoint blockade (CPB) improves melanoma outcomes, but many patients still do not respond. Tumor mutational burden (TMB) and tumor-infiltrating T cells are associated with response, and integrative models improve survival prediction. However, integrating immune/tumor-intrinsic features using data from a single assay (DNA/RNA) remains underexplored. Here, we analyze whole-exome and bulk RNA sequencing of tumors from new and published cohorts of 189 and 178 patients with melanoma receiving CPB, respectively. Using DNA, we calculate T cell and B cell burdens (TCB/BCB) from rearranged TCR/Ig sequences and find that patients with TMBhigh and TCBhigh or BCBhigh have improved outcomes compared to other patients. By combining pairs of immune- and tumor-expressed genes, we identify three gene pairs associated with response and survival, which validate in independent cohorts. The top model includes lymphocyte-expressed MAP4K1 and tumor-expressed TBX3. Overall, RNA or DNA-based models combining immune and tumor measures improve predictions of melanoma CPB outcomes.

Longitudinal characterization of circulating neutrophils uncovers distinct phenotypes associated with disease severity in hospitalized COVID-19 patients

Multiple studies have identified an association between neutrophils and COVID-19 disease severity; however, the mechanistic basis of this association remains incompletely understood. Here we collected 781 longitudinal blood samples from 306 hospitalized COVID-19 ⁺ patients, 78 COVID-19 ^âˆ' acute respiratory distress syndrome patients, and 8 healthy controls, and performed bulk RNA-sequencing of enriched neutrophils, plasma proteomics, cfDNA measurements and high throughput antibody profiling assays to investigate the relationship between neutrophil states and disease severity or death. We identified dynamic switches between six distinct neutrophil subtypes using non-negative matrix factorization (NMF) clustering. At days 3 and 7 post-hospitalization, patients with severe disease had an enrichment of a granulocytic myeloid derived suppressor cell-like state gene expression signature, while non-severe patients with resolved disease were enriched for a progenitor-like immature neutrophil state signature. Severe disease was associated with gene sets related to neutrophil degranulation, neutrophil extracellular trap (NET) signatures, distinct metabolic signatures, and enhanced neutrophil activation and generation of reactive oxygen species (ROS). We found that the majority of patients had a transient interferon-stimulated gene signature upon presentation to the emergency department (ED) defined here as Day 0, regardless of disease severity, which persisted only in patients who subsequently died. Humoral responses were identified as potential drivers of neutrophil effector functions, as enhanced antibody-dependent neutrophil phagocytosis and reduced NETosis was associated with elevated SARS-CoV-2-specific IgG1-to-IgA1 ratios in plasma of severe patients who survived. In vitro experiments confirmed that while patient-derived IgG antibodies mostly drove neutrophil phagocytosis and ROS production in healthy donor neutrophils, patient-derived IgA antibodies induced a predominant NETosis response. Overall, our study demonstrates neutrophil dysregulation in severe COVID-19 and a potential role for IgA-dominant responses in driving neutrophil effector functions in severe disease and mortality.

Genome-scale screens identify factors regulating tumor cell responses to natural killer cells

To systematically define molecular features in human tumor cells that determine their degree of sensitivity to human allogeneic natural killer (NK) cells, we quantified the NK cell responsiveness of hundreds of molecularly annotated 'DNA-barcoded' solid tumor cell lines in multiplexed format and applied genome-scale CRISPR-based gene-editing screens in several solid tumor cell lines, to functionally interrogate which genes in tumor cells regulate the response to NK cells. In these orthogonal studies, NK cell-sensitive tumor cells tend to exhibit 'mesenchymal-like' transcriptional programs; high transcriptional signature for chromatin remodeling complexes; high levels of B7-H6 (NCR3LG1); and low levels of HLA-E/antigen presentation genes. Importantly, transcriptional signatures of NK cell-sensitive tumor cells correlate with immune checkpoint inhibitor (ICI) resistance in clinical samples. This study provides a comprehensive map of mechanisms regulating tumor cell responses to NK cells, with implications for future biomarker-driven applications of NK cell immunotherapies.

KIR3DL3 Is an Inhibitory Receptor for HHLA2 that Mediates an Alternative Immunoinhibitory Pathway to PD1

Blockade of the PD1 pathway is a broadly effective cancer therapy, but additional immune-inhibitory pathways contribute to tumor immune evasion. HERV-H LTR-associating 2 (HHLA2; also known as B7H5 and B7H7) is a member of the B7 family of immunoregulatory ligands that mediates costimulatory effects through its interaction with the CD28 family member transmembrane and immunoglobulin domain containing 2 (TMIGD2). However, HHLA2 has also been known to have inhibitory effects on T cells. Here, we report that we have identified killer cell immunoglobulin-like receptor, three immunoglobulin domains and long cytoplasmic tail 3 (KIR3DL3) as an inhibitory receptor for HHLA2 in T cells and natural killer (NK) cells and have generated HHLA2 and KIR3DL3 antibodies that block the immune-inhibitory activity of HHLA2, preserving the costimulatory signal. It is known that HHLA2 is frequently expressed in several tumor types, including clear cell renal cell carcinoma (ccRCC). We found that HHLA2 expression was nonoverlapping with PDL1 expression in ccRCC, suggesting that HHLA2 mediates a mechanism of tumor immune evasion that is independent from PDL1. Blockade of both the PD1 and KIR3DL3 pathways may be a more effective way to reverse tumor immune evasion.See related Spotlight on p. 128.

The Lipogenic Regulator SREBP2 Induces Transferrin in Circulating Melanoma Cells and Suppresses Ferroptosis

Circulating tumor cells (CTC) are shed by cancer into the bloodstream, where a viable subset overcomes oxidative stress to initiate metastasis. We show that single CTCs from patients with melanoma coordinately upregulate lipogenesis and iron homeostasis pathways. These are correlated with both intrinsic and acquired resistance to BRAF inhibitors across clonal cultures of BRAF-mutant CTCs. The lipogenesis regulator SREBP2 directly induces transcription of the iron carrier Transferrin (TF), reducing intracellular iron pools, reactive oxygen species, and lipid peroxidation, thereby conferring resistance to inducers of ferroptosis. Knockdown of endogenous TF impairs tumor formation by melanoma CTCs, and their tumorigenic defects are partially rescued by the lipophilic antioxidants ferrostatin-1 and vitamin E. In a prospective melanoma cohort, presence of CTCs with high lipogenic and iron metabolic RNA signatures is correlated with adverse clinical outcome, irrespective of treatment regimen. Thus, SREBP2-driven iron homeostatic pathways contribute to cancer progression, drug resistance, and metastasis. SIGNIFICANCE: Through single-cell analysis of primary and cultured melanoma CTCs, we have uncovered intrinsic cancer cell heterogeneity within lipogenic and iron homeostatic pathways that modulates resistance to BRAF inhibitors and to ferroptosis inducers. Activation of these pathways within CTCs is correlated with adverse clinical outcome, pointing to therapeutic opportunities.This article is highlighted in the In This Issue feature, p. 521.

Abstract 6073: The lipogenic regulator SREBP induces Transferrin in circulating melanoma cells, suppressing their susceptibility to ferroptosis

Circulating tumor cells (CTCs) are shed by cancers into the bloodstream, where a viable subset overcomes oxidative stress to initiate metastatic outgrowth. Clonally derived cultured CTCs from patients with BRAF-mutant melanoma reveal upregulation of lipogenesis and iron homeostasis pathways, correlated with their baseline and acquired drug resistance. In CTCs, the lipogenesis regulator SREBP directly induces transcription of the iron carrier Transferrin (TF), thereby reducing intracellular reactive oxygen species (ROS) and lipid peroxidation, and conferring resistance to BRAF inhibitors and inducers of ferroptosis. Knockdown of endogenous TF impairs tumorigenesis by melanoma CTCs; their associated soft agar clonogenic defect is rescued by the lipophilic anti-oxidants Ferrostatin-1 or Vitamin E, and by cholesterol. Single cell RNA-seq of patient-derived melanoma CTCs identifies a subset with high lipogenic, iron metabolic and proliferative signatures, correlated with adverse clinical outcome. Together, the coordinated regulation of these SREBP-driven pathways contributes to cancer progression, drug resistance and metastasis.

Citation Format: Xin Hong, Whijae Roh, Ryan J. Sullivan, Keith H. Wong, Ben S. Wittner, HongShan Guo, Taronish D. Dubash, Moshe Sade-Feldman, Ben K. Wesley, Genevieve M. Boland, Dieuwke L. Marvin, Todd Bonesteel, Chenyue Lu, Elad Horwitz, François Aguet, Samuel S. Freeman, Katherine Calhoun, Michelle K. Jewett, Linda T. Nieman, Nir Hacohen, Anders M. Näär, David T. Ting, Mehmet Toner, Shannon L. Stott, Gad Getz, Shyamala Maheswaran, Daniel A. Haber. The lipogenic regulator SREBP induces Transferrin in circulating melanoma cells, suppressing their susceptibility to ferroptosis [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 6073.

Abstract 6670: Combined signals from tumor and immune cells predict outcomes of checkpoint inhibition in melanoma

Cancer immunotherapy with checkpoint blockade has improved survival and outcomes in melanoma, but still a majority of patients do not respond. Both high tumor mutation burden (TMB) and high T cell infiltration have been associated with response, but integrative models based on DNA or RNA assays have not been comprehensively explored and validated. Focusing on melanomas from patients receiving checkpoint blockade, we generated new and aggregated existing datasets of whole exome sequencing (WES) (n = 189 total) and bulk RNA sequencing (n = 154 total) to derive genomic and transcriptomic factors that predict survival and response to immunotherapy in melanoma.

We quantified T and B cell infiltrates using rearranged T cell receptor (TCR) and immunoglobulin (Ig) sequences, respectively, from DNA or RNA sequencing. High levels of rearranged TCR reads or rearranged Ig reads in RNA-seq were associated with survival (P = 0.0046, P = 0.015) and response (P = 0.0034, P = 0.047). We created RNA-based metrics of T and B cell burden (TCBRNA or BCBRNA) by normalizing the number of rearranged TCR reads by the total number of mapped reads. When we analyzed WES data in patients for whom DNA and RNA were extracted from the same region, we found that the TCBDNA correlated with TCBRNA (rho = 0.73) and BCBDNA with BCBRNA (rho = 0.41), demonstrating that the level of lymphocyte infiltration can be estimated using rearranged TCR or Ig reads from tumor WES alone.

We found that TCBDNA and BCBDNA both associated with survival (P = 0.0023 and 0.0089). In a combined model, patients with high TMB and high TCB DNA survived longer (P = 2.4e-4, HR = 2.68) and had a higher response rate (Fisher P = 0.028). This combined model was superior to models with TMB or TCBDNA alone. Similarly, patients with high TMB and high BCBDNA had longer survival and higher response rates (log-rank P = 0.0029, HR = 2.64, Fisher P = 0.015). We reanalyzed stage III/IV melanomas from TCGA and found that the TMB high, TCBDNA high subgroup had increased survival (P = 0.007).

Next, clustering of tumor transcriptomes identified 5 tumor subtypes based on melanocyte differentiation, immune infiltration and keratin levels. These melanoma subtypes were associated with survival outcomes after immunotherapy (P = 0.019). We found that TBX3, a tumor-expressed transcription factor enriched in poorly differentiated melanomas, was over-expressed among non-responders within the immune-infiltrated subtype and among all patients (P = 3.9e-4, P = 8.7e-5). Patients whose tumors had high immune infiltrate and low expression of TBX3 had longer survival (P = 1.6e-5, HR = 3.39), however this subgroup did not have longer survival in an independent cohort (n = 73, P = 0.10, HR = 2.63). In conclusion, we demonstrate both RNA-based (immune infiltrate and tumor subtype) and DNA-based metrics (TMB/TCB or TMB/BCB) can be used as pre-treatment predictors of survival after checkpoint blockade in melanoma.

Citation Format: Samuel S. Freeman, Moshe Sade-Feldman, Jaegil Kim, Chip Stewart, Arvind Ravi, Monica Arniella, Keren Yizhak, Ignaty Leshchiner, Liudmila Elagina, Oliver Spiro, Dimitri Livitz, Daniel Rosebrock, François Aguet, Jian Carrot-Zhang, Anna Gonye, Gavin Ha, Ziao Lin, Jonathan H. Chen, Dennie T. Frederick, Michal Barzily-Rokni, Marc R. Hammond, Hans Vitzthum, Shauna M. Blackmon, Yunxin J. Jiao, Donald P. Lawrence, Lyn M. Duncan, Anat Stemmer-Rachamimov, Jennifer A. Wargo, Keith T. Flaherty, Genevieve M. Boland, Ryan J. Sullivan, Matthew Meyerson, Gad Getz, Nir Hacohen. Combined signals from tumor and immune cells predict outcomes of checkpoint inhibition in melanoma [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 6670.

Sensitive Detection of Minimal Residual Disease in Patients Treated for Early-Stage Breast Cancer

PURPOSE

Existing cell-free DNA (cfDNA) methods lack the sensitivity needed for detecting minimal residual disease (MRD) following therapy. We developed a test for tracking hundreds of patient-specific mutations to detect MRD with a 1,000-fold lower error rate than conventional sequencing.

EXPERIMENTAL DESIGN

We compared the sensitivity of our approach to digital droplet PCR (ddPCR) in a dilution series, then retrospectively identified two cohorts of patients who had undergone prospective plasma sampling and clinical data collection: 16 patients with ER+/HER2- metastatic breast cancer (MBC) sampled within 6 months following metastatic diagnosis and 142 patients with stage 0 to III breast cancer who received curative-intent treatment with most sampled at surgery and 1 year postoperative. We performed whole-exome sequencing of tumors and designed individualized MRD tests, which we applied to serial cfDNA samples.

RESULTS

Our approach was 100-fold more sensitive than ddPCR when tracking 488 mutations, but most patients had fewer identifiable tumor mutations to track in cfDNA (median = 57; range = 2-346). Clinical sensitivity was 81% (n = 13/16) in newly diagnosed MBC, 23% (n = 7/30) at postoperative and 19% (n = 6/32) at 1 year in early-stage disease, and highest in patients with the most tumor mutations available to track. MRD detection at 1 year was strongly associated with distant recurrence [HR = 20.8; 95% confidence interval, 7.3-58.9]. Median lead time from first positive sample to recurrence was 18.9 months (range = 3.4-39.2 months).

CONCLUSIONS

Tracking large numbers of individualized tumor mutations in cfDNA can improve MRD detection, but its sensitivity is driven by the number of tumor mutations available to track.

Reversion and non-reversion mechanisms of resistance to PARP inhibitor or platinum chemotherapy in BRCA1/2-mutant metastatic breast cancer

BACKGROUND

Little is known about mechanisms of resistance to poly(adenosine diphosphate-ribose) polymerase inhibitors (PARPi) and platinum chemotherapy in patients with metastatic breast cancer and BRCA1/2 mutations. Further investigation of resistance in clinical cohorts may point to strategies to prevent or overcome treatment failure.

PATIENTS AND METHODS

We obtained tumor biopsies from metastatic breast cancer patients with BRCA1/2 deficiency before and after acquired resistance to PARPi or platinum chemotherapy. Whole exome sequencing was carried out on each tumor, germline DNA, and circulating tumor DNA. Tumors underwent RNA sequencing, and immunohistochemical staining for RAD51 foci on tumor sections was carried out for functional assessment of intact homologous recombination (HR).

RESULTS

Pre- and post-resistance tumor samples were sequenced from eight patients (four with BRCA1 and four with BRCA2 mutation; four treated with PARPi and four with platinum). Following disease progression on DNA-damaging therapy, four patients (50%) acquired at least one somatic reversion alteration likely to result in functional BRCA1/2 protein detected by tumor or circulating tumor DNA sequencing. Two patients with germline BRCA1 deficiency acquired genomic alterations anticipated to restore HR through increased DNA end resection: loss of TP53BP1 in one patient and amplification of MRE11A in another. RAD51 foci were acquired post-resistance in all patients with genomic reversion, consistent with reconstitution of HR. All patients whose tumors demonstrated RAD51 foci post-resistance were intrinsically resistant to subsequent lines of DNA-damaging therapy.

CONCLUSIONS

Genomic reversion in BRCA1/2 was the most commonly observed mechanism of resistance, occurring in four of eight patients. Novel sequence alterations leading to increased DNA end resection were seen in two patients, and may be targetable for therapeutic benefit. The presence of RAD51 foci by immunohistochemistry was consistent with BRCA1/2 protein functional status from genomic data and predicted response to later DNA-damaging therapy, supporting RAD51 focus formation as a clinically useful biomarker.

Genomic and immune profiling of pre-invasive lung adenocarcinoma

Adenocarcinoma in situ and minimally invasive adenocarcinoma are the pre-invasive forms of lung adenocarcinoma. The genomic and immune profiles of these lesions are poorly understood. Here we report exome and transcriptome sequencing of 98 lung adenocarcinoma precursor lesions and 99 invasive adenocarcinomas. We have identified EGFR, RBM10, BRAF, ERBB2, TP53, KRAS, MAP2K1 and MET as significantly mutated genes in the pre/minimally invasive group. Classes of genome alterations that increase in frequency during the progression to malignancy are revealed. These include mutations in TP53, arm-level copy number alterations, and HLA loss of heterozygosity. Immune infiltration is correlated with copy number alterations of chromosome arm 6p, suggesting a link between arm-level events and the tumor immune environment.

Identification and characterization of an alternative cancer-derived PD-L1 splice variant

Therapeutic blockade of the PD-1/PD-L1 axis is recognized as an effective treatment for numerous cancer types. However, only a subset of patients respond to this treatment, warranting a greater understanding of the biological mechanisms driving immune evasion via PD-1/PD-L1 signaling and other T-cell suppressive pathways. We previously identified a head and neck squamous cell carcinoma with human papillomavirus integration in the PD-L1 locus upstream of the transmembrane domain-encoding region, suggesting expression of a truncated form of PD-L1 (Parfenov et al., Proc Natl Acad Sci USA 111(43):15544-15549, 2014). In this study, we extended this observation by performing a computational analysis of 33 other cancer types as well as human cancer cell lines, and identified additional PD-L1 isoforms with an exon 4 enrichment expressed in 20 cancers and human cancer cell lines. We demonstrate that cancer cell lines with high expression levels of exon 4-enriched PD-L1 generate a secreted form of PD-L1. Further biochemical studies of exon 4-enriched PD-L1 demonstrated that this form is secreted and maintains the capacity to bind PD-1 as well as to serve as a negative regulator on T cell function, as measured by inhibition of IL-2 and IFNg secretion. Overall, we have demonstrated that truncated forms of PD-L1 exist in numerous cancer types, and have validated that truncated PD-L1 can be secreted and negatively regulate T cell function.

Tumor fraction in cell-free DNA as a biomarker in prostate cancer

BACKGROUND

Tumor content in circulating cell-free DNA (cfDNA) is a promising biomarker, but longitudinal dynamics of tumor-derived and non-tumor-derived cfDNA through multiple courses of therapy have not been well described.

METHODS

CfDNA from 663 plasma samples from 140 patients with castration-resistant prostate cancer (CRPC) was subject to sparse whole genome sequencing. Tumor fraction (TFx) estimated using the computational tool ichorCNA was correlated with clinical features and responses to therapy.

RESULTS

TFx associated with the number of bone metastases (median TFx = 0.014 with no bone metastases, 0.047 with 1-3 bone metastases, 0.190 for 4+ bone metastases; P < 0.0001) and with visceral metastases (P < 0.0001). In multivariable analysis, TFx remained associated with metastasis location (P = 0.042); TFx was positively correlated with alkaline phosphatase (P = 0.0227) and negatively correlated with hemoglobin (Hgb) (P < 0.001), but it was not correlated with prostate specific antigen (PSA) (P = 0.75). Tumor-derived and non-tumor-derived cfDNA track together and do not increase with generalized tissue damage from chemotherapy or radiation at the time scales examined. All new treatments that led to ≥30% PSA decline at 6 weeks were associated with TFx decline when baseline TFx was >7%; however, TFx in patients being subsequently maintained on secondary hormonal therapy was quite dynamic.

CONCLUSION

TFx correlates with clinical features associated with overall survival in CRPC, and TFx decline is a promising biomarker for initial therapeutic response.

TRIAL REGISTRATION

Dana-Farber/Harvard Cancer Center (DF/HCC) protocol no. 18-135.

FUNDING

Wong Family Award in Translational Oncology, Dana Farber Cancer Institute Medical Oncology grant, Gerstner Family Foundation, Janssen Pharmaceuticals Inc., and Koch Institute Support (core) grant P30-CA14051 from the National Cancer Institute (NCI).

Defining T Cell States Associated with Response to Checkpoint Immunotherapy in Melanoma

Treatment of cancer has been revolutionized by immune checkpoint blockade therapies. Despite the high rate of response in advanced melanoma, the majority of patients succumb to disease. To identify factors associated with success or failure of checkpoint therapy, we profiled transcriptomes of 16,291 individual immune cells from 48 tumor samples of melanoma patients treated with checkpoint inhibitors. Two distinct states of CD8+ T cells were defined by clustering and associated with patient tumor regression or progression. A single transcription factor, TCF7, was visualized within CD8+ T cells in fixed tumor samples and predicted positive clinical outcome in an independent cohort of checkpoint-treated patients. We delineated the epigenetic landscape and clonality of these T cell states and demonstrated enhanced antitumor immunity by targeting novel combinations of factors in exhausted cells. Our study of immune cell transcriptomes from tumors demonstrates a strategy for identifying predictors, mechanisms, and targets for enhancing checkpoint immunotherapy.

Structural Alterations Driving Castration-Resistant Prostate Cancer Revealed by Linked-Read Genome Sequencing

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

Abstract 3001: Broad/IBM Project: Discovery of treatment resistance mechanisms through use of liquid biopsy genomics services

The Broad/IBM Cancer Resistance Project has partnered with Broad Genomics to pilot the use of cutting edge sequencing technology for the analysis of cell free DNA in blood biopsies. Working closely with the Broad's Cancer Program, Broad Genomics has developed a suite of liquid biopsy sequencing products designed to provide optimal flexibility in conducting research studies with a broad range of applications including; biomarker discovery, treatment resistance monitoring, and detection of minimal residual disease (MRD) post-surgery. Cell-free DNA is extracted from the blood, and a dual unique-molecular-indexed library is created. From this library, low coverage whole genome (ultra-low-pass 0.1x coverage) data is generated to survey sample quality and evaluate the tumor fraction in the liquid specimen. Utilizing the same library, additional assays can be selected for processing based on the research aim (Targeted Panel Assays, MRD Detection or Whole Exomes). Since our approach utilizes the same genomic material for whole genome and targeted sequencing assays, it is possible to maximize the information learned from each valuable and limited liquid biopsy specimen. Our study design takes advantage of the discovery potential of combined tissue-based sequencing and serial liquid biopsy analysis to elucidate mechanisms of cancer resistance by tracking the evolution of clonal and subclonal populations in patients samples over time. This collaboration will utilize the ultra-low-pass sequencing and whole exome sequencing together with custom analysis pipelines to correlate the genomic events with patient clinical data. We aim to process 3,000 samples from 1,000 patients over the next 3 years. To date we have processed close to 500 samples through the ultra-low-pass pipeline and 100 samples through the whole exome sequencing pipeline (results to be provided).The ability to successfully investigate treatment resistant cancers from non-invasive liquid biopsies presents new opportunities for identifying markers, understanding dynamics and monitoring tumor dissemination and clonal evolution.

Citation Format: Gad Getz, Carrie Cibulskis, Ignaty Leshchiner, Megan Hanna, Dimitri Livitz, Kara Slowik, Chaya Levovitz, Filippo Utro, Kahn Rhrissorrakrai, Denisse Rotem, Gregory Gydush, Sarah C. Reed, Justin Rhoades, Gavin Ha, Samuel S. Freeman, Christopher Lo, Mark Fleharty, Justin Abreu, Katie Larkin, Michelle Cipicchio, Brendan Blumenstiel, Matt DeFelice, Jonna Grimsby, Susanna Hamilton, Niall Lennon, Viktor A. Adalsteinsson, Laxmi Parida. Broad/IBM Project: Discovery of treatment resistance mechanisms through use of liquid biopsy genomics services [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3001.

Abstract LB-225: Liquid biopsies identify trunk mutations and reflect multiple tumors in a patient

Introduction: Precision medicine approaches to guide therapy selection require routine sampling of tumors. However, tumor biopsies are not always accessible and may be confounded by spatial heterogeneity. Liquid biopsies, including analysis of cell-free DNA (cfDNA), present a non-invasive alternative which may reflect multiple tumors in the body. Previous studies have demonstrated exome-wide concordance between single-site tumor biopsies and cfDNA, but little is known about how cfDNA reflects multiple lesions within a patient. Here we sought to determine how cfDNA reflects the body-wide tumor phylogeny, which will inform the use of cfDNA for cancer precision medicine.

Methods: We identified 20 patients with pancreatic cancer who had undergone rapid autopsy. We then screened cfDNA tumor fraction and performed whole-exome sequencing of cfDNA and multiple tumor biopsies for 3 patients with cfDNA tumor fraction &gt;10%. We inferred the tumor phylogeny and then developed a statistical approach to deconvolute the contributions to cfDNA from tumor phylogenetic nodes. Finally, we determined whether shared trunk mutations could be detected in cfDNA and tumor biopsies.

Results: For each patient, we found mutations shared between all sites and cfDNA, including putative driver mutations. We found mutations which were clonal in multiple regions were detectable in cfDNA, whereas mutations private to individual sites were never clonal in cfDNA. Through our deconvolution analysis, we found that cfDNA could not be modeled as a simple linear combination of individual sites, but rather that cfDNA represented multiple nodes in the inferred phylogeny. For two pancreatic adenocarcinoma patients, the inferred ancestor of the metastases had high estimated contribution (&gt;70%) to cfDNA, while the ancestors of the primaries had lower contributions (&lt;10%). Next, we considered trunk mutations, which originate earliest in the tumor phylogenetic tree. When we analyzed precision for detection of trunk mutations, we found on average, 71% of clonal mutations in metastases were truncal, while only 55% of clonal mutations in primary tumors were truncal. Due to copy number deletions, not all trunk mutations were detected in metastases. Finally, on average, cfDNA had equal or better precision than 83% of primaries and 88% of metastases, suggesting cfDNA may provide more accurate trunk SSNV calls than tumor biopsies.

Conclusions: Through analyzing cfDNA and synchronous tumor biopsies from the same patient, we find trunk mutations are enriched in cfDNA as compared to the average single-site biopsy. We also predict that cfDNA represents multiple nodes in the inferred phylogeny. In cases where tumor biopsies are inaccessible, we demonstrate that cfDNA might be a promising alternative to detect trunk SSNVs. These results suggest that cfDNA may be complementary to tumor biopsies for disease monitoring and treatment selection in personalized medicine.

Citation Format: Samuel S. Freeman, Ziao Lin, Gavin Ha, Ignaty Leshchiner, Justin Rhoades, Dimitri Livitz, Daniel Rosebrock, Sarah C. Reed, Gregory Gydush, Christopher Lo, Denisse Rotem, Atish D. Choudhury, Daniel G. Stover, Heather A. Parsons, Jesse S. Boehm, J Christopher Love, Matthew Meyerson, Paul Grandgenett, Michael A. Hollingsworth, Viktor A. Adalsteinsson, Gad Getz. Liquid biopsies identify trunk mutations and reflect multiple tumors in a patient [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr LB-225.

Whole-exome sequencing of cell-free DNA and circulating tumor cells in multiple myeloma

Liquid biopsies including circulating tumor cells (CTCs) and cell-free DNA (cfDNA) have enabled minimally invasive characterization of many cancers, but are rarely analyzed together. Understanding the detectability and genomic concordance of CTCs and cfDNA may inform their use in guiding cancer precision medicine. Here, we report the detectability of cfDNA and CTCs in blood samples from 107 and 56 patients with multiple myeloma (MM), respectively. Using ultra-low pass whole-genome sequencing, we find both tumor fractions correlate with disease progression. Applying whole-exome sequencing (WES) to cfDNA, CTCs, and matched tumor biopsies, we find concordance in clonal somatic mutations (~99%) and copy number alterations (~81%) between liquid and tumor biopsies. Importantly, analyzing CTCs and cfDNA together enables cross-validation of mutations, uncovers mutations exclusive to either CTCs or cfDNA, and allows blood-based tumor profiling in a greater fraction of patients. Our study demonstrates the utility of analyzing both CTCs and cfDNA in MM.

Circulating tumor cells (CTCs) and cell-free DNA (cfDNA) enables characterization of a patient’s cancer. Here, the authors analyse CTCs, cfDNA, and tumor biopsies from multiple myeloma patients to show these approaches are complementary for mutation detection, together enabling a greater fraction of patient tumors to be profiled.

Association of Cell-Free DNA Tumor Fraction and Somatic Copy Number Alterations With Survival in Metastatic Triple-Negative Breast Cancer

Purpose Cell-free DNA (cfDNA) offers the potential for minimally invasive genome-wide profiling of tumor alterations without tumor biopsy and may be associated with patient prognosis. Triple-negative breast cancer (TNBC) is characterized by few mutations but extensive somatic copy number alterations (SCNAs), yet little is known regarding SCNAs in metastatic TNBC. We sought to evaluate SCNAs in metastatic TNBC exclusively via cfDNA and determine if cfDNA tumor fraction is associated with overall survival in metastatic TNBC. Patients and Methods In this retrospective cohort study, we identified 164 patients with biopsy-proven metastatic TNBC at a single tertiary care institution who received prior chemotherapy in the (neo)adjuvant or metastatic setting. We performed low-coverage genome-wide sequencing of cfDNA from plasma. Results Without prior knowledge of tumor mutations, we determined tumor fraction of cfDNA for 96.3% of patients and SCNAs for 63.9% of patients. Copy number profiles and percent genome altered were remarkably similar between metastatic and primary TNBCs. Certain SCNAs were more frequent in metastatic TNBCs relative to paired primary tumors and primary TNBCs in publicly available data sets The Cancer Genome Atlas and METABRIC, including chromosomal gains in drivers NOTCH2, AKT2, and AKT3. Prespecified cfDNA tumor fraction threshold of ≥ 10% was associated with significantly worse metastatic survival (median, 6.4 v 15.9 months) and remained significant independent of clinicopathologic factors (hazard ratio, 2.14; 95% CI, 1.4 to 3.8; P < .001). Conclusion We present the largest genomic characterization of metastatic TNBC to our knowledge, exclusively from cfDNA. Evaluation of cfDNA tumor fraction was feasible for nearly all patients, and tumor fraction ≥ 10% is associated with significantly worse survival in this large metastatic TNBC cohort. Specific SCNAs are enriched and prognostic in metastatic TNBC, with implications for metastasis, resistance, and novel therapeutic approaches.

Genomic Heterogeneity as a Barrier to Precision Medicine in Gastroesophageal Adenocarcinoma

Gastroesophageal adenocarcinoma (GEA) is a lethal disease where targeted therapies, even when guided by genomic biomarkers, have had limited efficacy. A potential reason for the failure of such therapies is that genomic profiling results could commonly differ between the primary and metastatic tumors. To evaluate genomic heterogeneity, we sequenced paired primary GEA and synchronous metastatic lesions across multiple cohorts, finding extensive differences in genomic alterations, including discrepancies in potentially clinically relevant alterations. Multiregion sequencing showed significant discrepancy within the primary tumor (PT) and between the PT and disseminated disease, with oncogene amplification profiles commonly discordant. In addition, a pilot analysis of cell-free DNA (cfDNA) sequencing demonstrated the feasibility of detecting genomic amplifications not detected in PT sampling. Lastly, we profiled paired primary tumors, metastatic tumors, and cfDNA from patients enrolled in the personalized antibodies for GEA (PANGEA) trial of targeted therapies in GEA and found that genomic biomarkers were recurrently discrepant between the PT and untreated metastases. Divergent primary and metastatic tissue profiling led to treatment reassignment in 32% (9/28) of patients. In discordant primary and metastatic lesions, we found 87.5% concordance for targetable alterations in metastatic tissue and cfDNA, suggesting the potential for cfDNA profiling to enhance selection of therapy.Significance: We demonstrate frequent baseline heterogeneity in targetable genomic alterations in GEA, indicating that current tissue sampling practices for biomarker testing do not effectively guide precision medicine in this disease and that routine profiling of metastatic lesions and/or cfDNA should be systematically evaluated. Cancer Discov; 8(1); 37-48. ©2017 AACR.See related commentary by Sundar and Tan, p. 14See related article by Janjigian et al., p. 49This article is highlighted in the In This Issue feature, p. 1.

Transcriptomic Characterization of SF3B1 Mutation Reveals Its Pleiotropic Effects in Chronic Lymphocytic Leukemia

Mutations in SF3B1, which encodes a spliceosome component, are associated with poor outcome in chronic lymphocytic leukemia (CLL), but how these contribute to CLL progression remains poorly understood. We undertook a transcriptomic characterization of primary human CLL cells to identify transcripts and pathways affected by SF3B1 mutation. Splicing alterations, identified in the analysis of bulk cells, were confirmed in single SF3B1-mutated CLL cells and also found in cell lines ectopically expressing mutant SF3B1. SF3B1 mutation was found to dysregulate multiple cellular functions including DNA damage response, telomere maintenance, and Notch signaling (mediated through KLF8 upregulation, increased TERC and TERT expression, or altered splicing of DVL2 transcript, respectively). SF3B1 mutation leads to diverse changes in CLL-related pathways.

Metagenomic Characterization of Microbial Communities In Situ Within the Deeper Layers of the Ileum in Crohn's Disease

BACKGROUND & AIMS

Microbial dysbiosis and aberrant host-microbe interactions in the gut are believed to contribute to the development and progression of Crohn's disease (CD). Microbiome studies in CD typically have focused on microbiota in feces or superficial mucosal layers of the colon because accessing DNA from deeper layers of the bowel is challenging. In this study, we analyzed the deep tissue microbiome in patients who underwent surgical resection of the small intestine.

METHODS

Paraffin blocks were obtained from 12 CD patients undergoing ileocecal resection, and healthy ileum samples (inflammatory bowel disease-free controls) were obtained from 12 patients undergoing surgery for right-sided colon cancer. Diseased and healthy-appearing ileum was identified using microscopy, and paraffin blocks were macrodissected using a core needle to specifically isolate DNA. Illumina Whole Genome Sequencing was used for microbial sequence identification and subsequent taxonomic classification using the PathSeq tool.

RESULTS

We observed significant differences between the microbiome of CD samples vs inflammatory bowel disease-free controls, including depletion of Bacteroidetes and Clostridia. Notably, microbial composition at the phyla level did not differ markedly between healthy and diseased areas of CD patients. However, we observed enrichment of potentially pathogenic organisms at the species level.

CONCLUSIONS

Our study showed dysbiosis within deeper layers of the ileum of CD patients, specifically enrichment of enterotoxigenic Staphylococcus aureus and an environmental Mycobacterium species not described previously. Future studies with larger cohort sizes are warranted to confirm these findings. Studies would benefit from effective microbial DNA extraction methods from paraffin sections and host nucleic acid depletion approaches to increase microbial read coverage.

Abstract B2-21: Identification of somatic RNA splicing alterations in human cancers

Recent whole-exome sequencing studies have found that somatic mutations frequently occur in splicing factors across multiple cancer types, supporting the need to systematically and globally characterize splicing alterations across human cancers. Somatic alterations that affect RNA processing of cancer genes may be caused by multiple mechanisms, including mutations in splicing factors, mutations in cis-acting splice sites and other regulatory sequences, alterations in methylation, or alterations in chromatin structure. For a global look at somatic RNA splicing alterations, we are applying the Cancer Outlier Profile Analysis (COPA) approach on ~7,000 cancer transcriptomes from 12 tumor types in The Cancer Genome Atlas (TCGA) and investigating the underlying somatic mutations that cause these splicing alterations using matched whole-exome, whole-genome, and methylation data from these samples.

To identify and quantify alternative splicing in RNA-Seq data, including unannotated splicing events, we have further developed a computational pipeline called JuncBASE. To distinguish between cancer-specific splicing alterations and normal transcriptome variation, we are utilizing ~700 RNA-Seq libraries from healthy individuals from the Genotype-Tissue Expression (GTEx) project.

COPA analysis of outlier splicing events in lung adenocarcinomas and glioblastomas successfully identified known altered splicing events in MET and EGFR, respectively, that are caused by DNA level somatic mutations. The genomic mechanisms for novel somatic splicing events in known cancer genes such as FGFR3, as well as genes previously uncharacterized in cancer genomes, are currently being investigated.

To identify splicing events that may be novel somatic driver alterations, these events have been profiled using RNA-Seq data from the Cancer Cell Line Encyclopedia and are being used as biomarkers to identify genetic vulnerabilities in high-throughput shRNA screens. This work will have a significant impact on the identification of somatic splicing events that contribute to cancer pathogenesis.

Citation Format: Angela N. Brooks, Yawei Ge, Kevin Chau, Samuel S. Freeman, Gordon Saksena, Chandra Sekhar Pedamallu, Matthew Meyerson. Identification of somatic RNA splicing alterations in human cancers. [abstract]. In: Proceedings of the AACR Special Conference on Computational and Systems Biology of Cancer; Feb 8-11 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 2):Abstract nr B2-21.

Genomic aberrations in cervical adenocarcinomas in Hong Kong Chinese women

Although the rates of cervical squamous cell carcinoma have been declining, the rates of cervical adenocarcinoma are increasing in some countries. Outcomes for advanced cervical adenocarcinoma remain poor. Precision mapping of genetic alterations in cervical adenocarcinoma may enable better selection of therapies and deliver improved outcomes when combined with new sequencing diagnostics. We present whole-exome sequencing results from 15 cervical adenocarcinomas and paired normal samples from Hong Kong Chinese women. These data revealed a heterogeneous mutation spectrum and identified several frequently altered genes including FAT1, ARID1A, ERBB2 and PIK3CA. Exome sequencing identified human papillomavirus (HPV) sequences in 13 tumors in which the HPV genome might have integrated into and hence disrupted the functions of certain exons, raising the possibility that HPV integration can alter pathways other than p53 and pRb. Together, these provisionary data suggest the potential for individualized therapies for cervical adenocarcinoma based on genomic information.

Abstract 4692: Relationships between somatic genomic alterations, tumor stage and progression-free survival in cervical cancer

Background:

Cervical cancer is a major public health problem worldwide. We have recently identified novel significantly recurrent somatic mutations in HLA-B, ERBB2 and MAPK1 in cervical squamous cell carcinomas. However, the significance of somatic mutations and copy number alterations for clinical phenotype in cervical cancer is not well understood. This study seeks to identify relationships between somatic genomic alterations, epidemiological exposures, tumor stage (FIGO) and progression-free survival in cervical carcinomas.

Methods:

Cervical tumors were surgically resected or biopsied from 100 Norwegian women with tumor stages I - IV. Patients were subsequently followed for 0-109 months (mean=25 months). Whole exome sequencing (WES) was performed on DNA extracted from tumors and corresponding normal blood. Somatic single nucleotide variants and small insertion/deletions were identified by the MuTect and Indelocator algorithms respectively. The ABSOLUTE algorithm was used to classify mutations as clonal or subclonal. Somatic copy number (CN) data were derived from WES data using the CapSeg algorithm, and significantly recurrent CN alterations were identified by GISTIC2.0 analysis (q&lt;0.25). HPV typing was done by the multiplex fluorescent f-HPV DNA and MassARRAY assays. The log-rank test was used to compare survival curves.

Results:

Non-localized tumors (FIGO stages ≥ II) were associated with focal amplification of the FGFR2 gene on chromosomal cytoband 10q26 (GISTIC q = 0.18531). Indeed, 6 of 8 (75%) tumors with FGFR2 amplification were non-localized, in contrast to 16 of 92 (17%) tumors without FGFR2 amplification (p = 0.001). In addition, patients with somatic ERBB2 mutations and/or amplifications (p = 0.04), somatic TP53 mutations and/or deletions (p = 0.04), or infection with multiple HPV types (p = 0.02) had poorer prognosis (progression-free survival). We also observed a trend for higher frequency of subclonal driver events in patients with poorer survival (p = 0.07).

Conclusion:

We have identified novel relationships between somatic genomic alterations, tumor stage and patient prognosis in cervical cancer. Our data suggest a potential for exploring FGFR2 inhibition in non-localized cervical carcinomas with FGFR2 alterations in a clinical trial context.

Citation Format: Akinyemi I. Ojesina, Bradley Murray, Line Bjorge, Kathrine Woie, Camilla Krakstad, Lee Lichtenstein, Chandra Sekhar Pedamallu, Amaro Taylor-Weiner, Samuel S. Freeman, Andrew D. Cherniack, Michael S. Lawrence, Kristian Cibulskis, Scott L. Carter, Heather Walline, Thomas E. Carey, Olav K. Vintermyr, Bjorn Bertelsen, Christopher P. Crum, Gad Getz, Matthew Meyerson, Helga B. Salvesen. Relationships between somatic genomic alterations, tumor stage and progression-free survival in cervical cancer. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4692. doi:10.1158/1538-7445.AM2014-4692

Landscape of genomic alterations in cervical carcinomas

Cervical cancer is responsible for 10-15% of cancer-related deaths in women worldwide. The aetiological role of infection with high-risk human papilloma viruses (HPVs) in cervical carcinomas is well established. Previous studies have also implicated somatic mutations in PIK3CA, PTEN, TP53, STK11 and KRAS as well as several copy-number alterations in the pathogenesis of cervical carcinomas. Here we report whole-exome sequencing analysis of 115 cervical carcinoma-normal paired samples, transcriptome sequencing of 79 cases and whole-genome sequencing of 14 tumour-normal pairs. Previously unknown somatic mutations in 79 primary squamous cell carcinomas include recurrent E322K substitutions in the MAPK1 gene (8%), inactivating mutations in the HLA-B gene (9%), and mutations in EP300 (16%), FBXW7 (15%), NFE2L2 (4%), TP53 (5%) and ERBB2 (6%). We also observe somatic ELF3 (13%) and CBFB (8%) mutations in 24 adenocarcinomas. Squamous cell carcinomas have higher frequencies of somatic nucleotide substitutions occurring at cytosines preceded by thymines (Tp*C sites) than adenocarcinomas. Gene expression levels at HPV integration sites were statistically significantly higher in tumours with HPV integration compared with expression of the same genes in tumours without viral integration at the same site. These data demonstrate several recurrent genomic alterations in cervical carcinomas that suggest new strategies to combat this disease.

Sequence-based discovery of Bradyrhizobium enterica in cord colitis syndrome

BACKGROUND

Immunosuppression is associated with a variety of idiopathic clinical syndromes that may have infectious causes. It has been hypothesized that the cord colitis syndrome, a complication of umbilical-cord hematopoietic stem-cell transplantation, is infectious in origin.

METHODS

We performed shotgun DNA sequencing on four archived, paraffin-embedded endoscopic colon-biopsy specimens obtained from two patients with cord colitis. Computational subtraction of human and known microbial sequences and assembly of residual sequences into a bacterial draft genome were performed. We used polymerase-chain-reaction (PCR) assays and fluorescence in situ hybridization to determine whether the corresponding bacterium was present in additional patients and controls.

RESULTS

DNA sequencing of the biopsy specimens revealed more than 2.5 million sequencing reads that did not match known organisms. These sequences were computationally assembled into a 7.65-Mb draft genome showing a high degree of homology with genomes of bacteria in the bradyrhizobium genus. The corresponding newly discovered bacterium was provisionally named Bradyrhizobium enterica. PCR identified B. enterica nucleotide sequences in biopsy specimens from all three additional patients with cord colitis whose samples were tested, whereas B. enterica sequences were absent in samples obtained from healthy controls and patients with colon cancer or graft-versus-host disease.

CONCLUSIONS

We assembled a novel bacterial draft genome from the direct sequencing of tissue specimens from patients with cord colitis. Association of these sequences with cord colitis suggests that B. enterica may be an opportunistic human pathogen. (Funded by the National Cancer Institute and others.)

Computer-assisted analysis of particulate axoplasmic flow in organized CNS tissue cultures

Using organized nervous system tissue cultures of embryonic mouse spinal cord, time-lapse motion pictures of intra-axonal particle movements were made using Nomarski optics. The paths of 272 particles were followed and analyzed using digital computer methods. Particles had a bidirectional saltatory motility with net anterograde and retrograde velocities which were shown to be the same (about 1 micron per second). When anterograde moving and retrograde moving particles were analyzed separately for anterograde and retrograde saltations, the velocities of saltation in the direction of net movement were found to be 1.5 times greater than saltations against the direction of net movement. These differences were statistically significant. In addition, it was shown that there were regions within the axon where variations in particle motion were similar from particle to particle as they passed, and that more disturbed motion was found to occur near intra-axonal objects which were judged to be mitochondria. No decrease of particle velocity was noted in regions away from mitochondria nor was in increase in velocity noted near them. Computer-drawn particle plots illustrate the various facets of particle motion.

Experimental acetylcholine blockade of the neuromuscular junction. Effects on end plate and muscle fiber ultrastructure

To investigate the effect of acetylcholine on the formation and maintenance of the end plate, presynaptic (hemicholinium-3 and botulinum toxin) and postsynaptic (curare and alpha-bungarotoxin) blocking agents were injected into the amniotic sac or chorioallantoic circulation of chick embryos during different stages of their development (days 3 to 11, 3 to 18, and 12 to 18). Tntercostal or anterior and posterior latissimus dorsi muscles were examined by electron microscopy. Regardless of the duration of treatment, the type of agent, or the type of muscle, nerve ingrowth into muscle and end-plate formation were not arrested and the ultrastructure of nerve terminals and postsynaptic regions was normal. By contrast, muscle fiber ultrastructure was affected by acetylcholine blockade if continued beyond day 11. By day 11 both treated and control muscles were composed primarily of myotubes. Muscles treated beyond day 11 showed impaired fiber differentiation and persistence of many myotubes. In addition, many fiber regions contained unassembled fine filaments intermingled with other organelles. In hemicholinium-3 treated muscles ther were also mitochondrial abnormalities and an increased number of lipid droplets. The findings suggest that the altered structure of the neuromusclar junction in myasthenia gravis is not due to the lack of a trophic effect of acetylcholine. The impaired fiber maturation in the chick embryos may be related to the lack of acetylcholine mediated fiber activity.