DeTiN: overcoming tumor-in-normal contamination

Randomized Phase II Study Evaluating the Addition of Pembrolizumab to Radium-223 in Metastatic Castration-resistant Prostate Cancer

The checkpoint immunotherapeutic pembrolizumab induces responses in a small minority of patients with metastatic castration-resistant prostate cancer (mCRPC). Radium-223 (R223) may increase immunogenicity of bone metastases and increase pembrolizumab (P) activity. In a randomized phase II study, we assessed the effect of R223+P compared with R223 on tumor immune infiltration, safety, and clinical outcomes in patients with mCRPC. The primary endpoint was differences in CD4+ and CD8+ T-cell infiltrate in 8-week versus baseline bone metastasis biopsies; secondary endpoints were safety, radiographic progression-free survival (rPFS), and overall survival (OS). Of the 42 treated patients (29 R223+P, 13 R223), 18 R223+P and 8 R223 patients had evaluable paired tumor biopsies. Median fold-change of CD4+ T cells was -0.7 (range: -9.3 to 4.7) with R223+P and 0.1 (-11.1 to 3.7) with R223 (P = 0.66); for CD8+ T cells, median fold-change was -0.6 (-7.4 to 5.3) with R223+P and -1.3 (-3.1 to 4.8) with R223 (P = 0.66). Median rPFS and OS was 6.1 (95% confidence interval: 2.7-11.0) and 16.9 months [12.7-not reached (NR)], respectively, with R223+P and 5.7 (2.6-NR) and 16.0 (9.0-NR), respectively, with R223. Although R223+P was well tolerated with no unexpected toxicity, the combination did not improve efficacy. High-dimensional flow cytometry demonstrated minimal immune modulation with R223, whereas R223+P induced CTLA-4 expression on circulating CD4+ T cells. Clinical responders possessed lower circulating frequencies of Ki67+ T and myeloid cells at baseline and higher circulating frequencies of TIM-3+ T and myeloid cells by week 9. Although R223+P did not induce T-cell infiltration into the tumor microenvironment, exhaustion of induced peripheral T-cell immune responses may dampen the combination's clinical activity.

A comprehensive performance evaluation, comparison, and integration of computational methods for detecting and estimating cross-contamination of human samples in cancer next-generation sequencing analysis

Cross-sample contamination is one of the major issues in next-generation sequencing (NGS)-based molecular assays. This type of contamination, even at very low levels, can significantly impact the results of an analysis, especially in the detection of somatic alterations in tumor samples. Several contamination identification tools have been developed and implemented as a crucial quality-control step in the routine NGS bioinformatic pipeline. However, no study has been published to comprehensively and systematically investigate, evaluate, and compare these computational methods in the cancer NGS analysis. In this study, we comprehensively investigated nine state-of-the-art computational methods for detecting cross-sample contamination. To explore their application in cancer NGS analysis, we further compared the performance of five representative tools by qualitative and quantitative analyses using in silico and simulated experimental NGS data. The results showed that Conpair achieved the best performance for identifying contamination and predicting the level of contamination in solid tumors NGS analysis. Moreover, based on Conpair, we developed a Python script, Contamination Source Predictor (ConSPr), to identify the source of contamination. We anticipate that this comprehensive survey and the proposed tool for predicting the source of contamination will assist researchers in selecting appropriate cross-contamination detection tools in cancer NGS analysis and inspire the development of computational methods for detecting sample cross-contamination and identifying its source in the future.

Allosteric PI3Kα Inhibition Overcomes On-target Resistance to Orthosteric Inhibitors Mediated by Secondary PIK3CA Mutations

PIK3CA mutations occur in ∼8% of cancers, including ∼40% of HR-positive breast cancers, where the PI3K-alpha (PI3Kα)-selective inhibitor alpelisib is FDA approved in combination with fulvestrant. Although prior studies have identified resistance mechanisms, such as PTEN loss, clinically acquired resistance to PI3Kα inhibitors remains poorly understood. Through serial liquid biopsies and rapid autopsies in 39 patients with advanced breast cancer developing acquired resistance to PI3Kα inhibitors, we observe that 50% of patients acquire genomic alterations within the PI3K pathway, including PTEN loss and activating AKT1 mutations. Notably, although secondary PIK3CA mutations were previously reported to increase sensitivity to PI3Kα inhibitors, we identified emergent secondary resistance mutations in PIK3CA that alter the inhibitor binding pocket. Some mutations had differential effects on PI3Kα-selective versus pan-PI3K inhibitors, but resistance induced by all mutations could be overcome by the novel allosteric pan-mutant-selective PI3Kα-inhibitor RLY-2608. Together, these findings provide insights to guide strategies to overcome resistance in PIK3CA-mutated cancers.

SIGNIFICANCE

In one of the largest patient cohorts analyzed to date, this study defines the clinical landscape of acquired resistance to PI3Kα inhibitors. Genomic alterations within the PI3K pathway represent a major mode of resistance and identify a novel class of secondary PIK3CA resistance mutations that can be overcome by an allosteric PI3Kα inhibitor. See related commentary by Gong and Vanhaesebroeck, p. 204 . See related article by Varkaris et al., p. 240 . This article is featured in Selected Articles from This Issue, p. 201.

Exosomal lncCRLA is predictive for the evolvement and development of lung adenocarcinoma

Lung adenocarcinoma, the most common histological subtype of non-small cell lung cancer, exhibits heterogeneity that enables adaptability, limits therapeutic success, and remains incompletely understood. Our team uncovers that lncRNA related to chemotherapy resistance in lung adenocarcinoma (lncCRLA) is preferentially expressed in lung adenocarcinoma cells with the mesenchymal phenotype. lncCRLA can not enhance chemotherapy resistance in lung adenocarcinoma due to its binding to RIPK1 in exosomes, which is released into intercellular media and transferred by exosomes from mesenchymal-like to epithelial-like cells. However, plasmatic lncCRLA corresponding to tissue lncCRLA functions as a preferred biomarker to reflect the response to chemotherapy and disease progression of lung adenocarcinoma. Through single-cell sequencing, RNA-Mutect technique and spatial transcriptomics, a handful of hybrid EMT cells with elevated lncCRLA are characterized as the origin of lung adenocarcinoma, which are indiscriminated from hybrid EMT cells by the in-depth sequencing. Plasmatic lncCRLA is properly predictive for the preinvasive lesion of lung adenocarcinoma that would evolve to invasive lesion. That notion is confirmed by a brand-new transgenic mouse model in which EMT is tracked by Cre and Dre system. Dasatinib is potential to hinder the spontaneous progression from preinvasive to invasive lesion of lung adenocarcinoma. Together, plasmatic lncCRLA is defined as a brand-new circulating biomarker to predict the occurrence and evolvement of lung adenocarcinoma, a light for early detection of lung adenocarcinoma.

Computational validation of clonal and subclonal copy number alterations from bulk tumor sequencing using CNAqc

Copy number alterations (CNAs) are among the most important genetic events in cancer, but their detection from sequencing data is challenging because of unknown sample purity, tumor ploidy, and general intra-tumor heterogeneity. Here, we present CNAqc, an evolution-inspired method to perform the computational validation of clonal and subclonal CNAs detected from bulk DNA sequencing. CNAqc is validated using single-cell data and simulations, is applied to over 4000 TCGA and PCAWG samples, and is incorporated into the validation process for the clinically accredited bioinformatics pipeline at Genomics England. CNAqc is designed to support automated quality control procedures for tumor somatic data validation.

Clinical application of tumour-in-normal contamination assessment from whole genome sequencing

The unexpected contamination of normal samples with tumour cells reduces variant detection sensitivity, compromising downstream analyses in canonical tumour-normal analyses. Leveraging whole-genome sequencing data available at Genomics England, we develop a tool for normal sample contamination assessment, which we validate in silico and against minimal residual disease testing. From a systematic review of [Formula: see text] patients with haematological malignancies and sarcomas, we find contamination across a range of cancer clinical indications and DNA sources, with highest prevalence in saliva samples from acute myeloid leukaemia patients, and sorted CD3+ T-cells from myeloproliferative neoplasms. Further exploration reveals 108 hotspot mutations in genes associated with haematological cancers at risk of being subtracted by standard variant calling pipelines. Our work highlights the importance of contamination assessment for accurate somatic variants detection in research and clinical settings, especially with large-scale sequencing projects being utilised to deliver accurate data from which to make clinical decisions for patient care.

The role of APOBEC3B in lung tumor evolution and targeted cancer therapy resistance

In this study, the impact of the apolipoprotein B mRNA-editing catalytic subunit-like (APOBEC) enzyme APOBEC3B (A3B) on epidermal growth factor receptor (EGFR)-driven lung cancer was assessed. A3B expression in EGFR mutant (EGFRmut) non-small-cell lung cancer (NSCLC) mouse models constrained tumorigenesis, while A3B expression in tumors treated with EGFR-targeted cancer therapy was associated with treatment resistance. Analyses of human NSCLC models treated with EGFR-targeted therapy showed upregulation of A3B and revealed therapy-induced activation of nuclear factor kappa B (NF-κB) as an inducer of A3B expression. Significantly reduced viability was observed with A3B deficiency, and A3B was required for the enrichment of APOBEC mutation signatures, in targeted therapy-treated human NSCLC preclinical models. Upregulation of A3B was confirmed in patients with NSCLC treated with EGFR-targeted therapy. This study uncovers the multifaceted roles of A3B in NSCLC and identifies A3B as a potential target for more durable responses to targeted cancer therapy.

Clinical Outcomes and Evolution of Clonal Hematopoiesis in Patients with Newly Diagnosed Multiple Myeloma

Clonal hematopoiesis (CH) at time of autologous stem cell transplant (ASCT) has been shown to be associated with decreased overall survival (OS) and progression-free survival (PFS) in patients with multiple myeloma not receiving immunomodulatory drugs (IMiD). However, the significance of CH in newly diagnosed patients, including transplant ineligible patients, and its effect on clonal evolution during multiple myeloma therapy in the era of novel agents, has not been well studied. Using our new algorithm to differentiate tumor and germline mutations from CH, we detected CH in approximately 10% of 986 patients with multiple myeloma from the Clinical Outcomes in MM to Personal Assessment of Genetic Profile (CoMMpass) cohort (40/529 transplanted and 59/457 non-transplanted patients). CH was associated with increased age, risk of recurrent bacterial infections and cardiovascular disease. CH at time of multiple myeloma diagnosis was not associated with inferior OS or PFS regardless of undergoing ASCT, and all patients benefited from IMiD-based therapies, irrespective of the presence of CH. Serial sampling of 52 patients revealed the emergence of CH over a median of 3 years of treatment, increasing its prevalence to 25%, mostly with DNMT3A mutations.

SIGNIFICANCE

Using our algorithm to differentiate tumor and germline mutations from CH mutations, we detected CH in approximately 10% of patients with newly diagnosed myeloma, including both transplant eligible and ineligible patients. Receiving IMiDs improved outcomes irrespective of CH status, but the prevalence of CH significantly rose throughout myeloma-directed therapy.

Genomic Tumor Correlates of Clinical Outcomes Following Organ-Sparing Chemoradiation Therapy for Bladder Cancer

PURPOSE

There is an urgent need for biomarkers of radiation response in organ-sparing therapies. Bladder preservation with trimodality therapy (TMT), consisting of transurethral tumor resection followed by chemoradiation, is an alternative to radical cystectomy for muscle-invasive bladder cancer (MIBC), but molecular determinants of response are poorly understood.

EXPERIMENTAL DESIGN

We characterized genomic and transcriptomic features correlated with long-term response in a single institution cohort of patients with MIBC homogeneously treated with TMT. Pretreatment tumors from 76 patients with MIBC underwent whole-exome sequencing; 67 underwent matched transcriptomic profiling. Molecular features were correlated with clinical outcomes including modified bladder-intact event-free survival (mBI-EFS), a composite endpoint that reflects long-term cancer control with bladder preservation.

RESULTS

With a median follow-up of 74.6 months in alive patients, 37 patients had favorable long-term response to TMT while 39 had unfavorable long-term response. Tumor mutational burden was not associated with outcomes after TMT. DNA damage response gene alterations were associated with improved locoregional control and mBI-EFS. Of these alterations, somatic ERCC2 mutations stood out as significantly associated with favorable long-term outcomes; patients with ERCC2 mutations had significantly improved mBI-EFS [HR, 0.15; 95% confidence interval (CI), 0.06-0.37; P = 0.030] and improved BI-EFS, an endpoint that includes all-cause mortality (HR, 0.33; 95% CI, 0.15-0.68; P = 0.044). ERCC2 mutant bladder cancer cell lines were significantly more sensitive to concurrent cisplatin and radiation treatment in vitro than isogenic ERCC2 wild-type cells.

CONCLUSIONS

Our data identify ERCC2 mutation as a candidate biomarker associated with sensitivity and long-term response to chemoradiation in MIBC. These findings warrant validation in independent cohorts.

TSC2 inactivation, low mutation burden and high macrophage infiltration characterise hepatic angiomyolipomas

AIMS

Although TSC1 or TSC2 inactivating mutations that lead to mTORC1 hyperactivation have been reported in hepatic angiomyolipomas (hAML), the role of other somatic genetic events that may contribute to hAML development is unknown. There are also limited data regarding the tumour microenvironment (TME) of hAML. The aim of the present study was to identify other somatic events in genomic level and changes in TME that contribute to tumorigenesis in hAML.

METHODS AND RESULTS

In this study, we performed exome sequencing in nine sporadic hAML tumours and deep-coverage targeted sequencing for TSC2 in three additional hAML. Immunohistochemistry and multiplex immunofluorescence were carried out for 15 proteins to characterise the tumour microenvironment and assess immune cell infiltration. Inactivating somatic variants in TSC2 were identified in 10 of 12 (83%) cases, with a median allele frequency of 13.6%. Five to 18 somatic variants (median number: nine, median allele frequency 21%) not in TSC1 or TSC2 were also identified, mostly of uncertain clinical significance. Copy number changes were rare, but detection was impaired by low tumour purity. Immunohistochemistry demonstrated numerous CD68+ macrophages of distinct appearance from Küpffer cells. Multiplex immunofluorescence revealed low numbers of exhausted PD-1+/PD-L1+, FOXP3+ and CD8+ T cells.

CONCLUSION

hAML tumours have consistent inactivating mutations in TSC2 and have a low somatic mutation rate, similar to other TSC-associated tumours. Careful histological review, standard IHC and multiplex immunofluorescence demonstrated marked infiltration by non-neoplastic inflammatory cells, mostly macrophages.

Pan-cancer proteogenomics connects oncogenic drivers to functional states

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Targeted in silico characterization of fusion transcripts in tumor and normal tissues via FusionInspector

Here, we present FusionInspector for in silico characterization and interpretation of candidate fusion transcripts from RNA sequencing (RNA-seq) and exploration of their sequence and expression characteristics. We applied FusionInspector to thousands of tumor and normal transcriptomes and identified statistical and experimental features enriched among biologically impactful fusions. Through clustering and machine learning, we identified large collections of fusions potentially relevant to tumor and normal biological processes. We show that biologically relevant fusions are enriched for relatively high expression of the fusion transcript, imbalanced fusion allelic ratios, and canonical splicing patterns, and are deficient in sequence microhomologies between partner genes. We demonstrate that FusionInspector accurately validates fusion transcripts in silico and helps characterize numerous understudied fusions in tumor and normal tissue samples. FusionInspector is freely available as open source for screening, characterization, and visualization of candidate fusions via RNA-seq, and facilitates transparent explanation and interpretation of machine-learning predictions and their experimental sources.

The peptide woods are lovely, dark and deep: Hunting for novel cancer antigens

Harnessing the patient's immune system to control a tumor is a proven avenue for cancer therapy. T cell therapies as well as therapeutic vaccines, which target specific antigens of interest, are being explored as treatments in conjunction with immune checkpoint blockade. For these therapies, selecting the best suited antigens is crucial. Most of the focus has thus far been on neoantigens that arise from tumor-specific somatic mutations. Although there is clear evidence that T-cell responses against mutated neoantigens are protective, the large majority of these mutations are not immunogenic. In addition, most somatic mutations are unique to each individual patient and their targeting requires the development of individualized approaches. Therefore, novel antigen types are needed to broaden the scope of such treatments. We review high throughput approaches for discovering novel tumor antigens and some of the key challenges associated with their detection, and discuss considerations when selecting tumor antigens to target in the clinic.

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.

Inferring early genetic progression in cancers with unobtainable premalignant disease

Analysis of premalignant tissue has identified the typical order of somatic events leading to invasive tumors in several cancer types. For other cancers, premalignant tissue is unobtainable, leaving genetic progression unknown. Here, we demonstrate how to infer progression from exome sequencing of primary tumors. Our computational method, PhylogicNDT, recapitulated the previous experimentally determined genetic progression of human papillomavirus-negative (HPV^-) head and neck squamous cell carcinoma (HNSCC). We then evaluated HPV⁺ HNSCC, which lacks premalignant tissue, and uncovered its previously unknown progression, identifying early drivers. We converted relative timing estimates of driver mutations and HPV integration to years before diagnosis based on a clock-like mutational signature. We associated the timing of transitions to aneuploidy with increased intratumor genetic heterogeneity and shorter overall survival. Our approach can establish previously unknown early genetic progression of cancers with unobtainable premalignant tissue, supporting development of experimental models and methods for early detection, interception and prognostication.

Workflow enabling deepscale immunopeptidome, proteome, ubiquitylome, phosphoproteome, and acetylome analyses of sample-limited tissues

Serial multi-omic analysis of proteome, phosphoproteome, and acetylome provides insights into changes in protein expression, cell signaling, cross-talk and epigenetic pathways involved in disease pathology and treatment. However, ubiquitylome and HLA peptidome data collection used to understand protein degradation and antigen presentation have not together been serialized, and instead require separate samples for parallel processing using distinct protocols. Here we present MONTE, a highly sensitive multi-omic native tissue enrichment workflow, that enables serial, deep-scale analysis of HLA-I and HLA-II immunopeptidome, ubiquitylome, proteome, phosphoproteome, and acetylome from the same tissue sample. We demonstrate that the depth of coverage and quantitative precision of each 'ome is not compromised by serialization, and the addition of HLA immunopeptidomics enables the identification of peptides derived from cancer/testis antigens and patient specific neoantigens. We evaluate the technical feasibility of the MONTE workflow using a small cohort of patient lung adenocarcinoma tumors.

Distinct Mechanisms of Mismatch-Repair Deficiency Delineate Two Modes of Response to Anti-PD-1 Immunotherapy in Endometrial Carcinoma

Mismatch repair-deficient (MMRd) cancers have varied responses to immune-checkpoint blockade (ICB). We conducted a phase II clinical trial of the PD-1 inhibitor pembrolizumab in 24 patients with MMRd endometrial cancer (NCT02899793). Patients with mutational MMRd tumors (6 patients) had higher response rates and longer survival than those with epigenetic MMRd tumors (18 patients). Mutation burden was higher in tumors with mutational MMRd compared with epigenetic MMRd; however, within each category of MMRd, mutation burden was not correlated with ICB response. Pretreatment JAK1 mutations were not associated with primary resistance to pembrolizumab. Longitudinal single-cell RNA-seq of circulating immune cells revealed contrasting modes of antitumor immunity for mutational versus epigenetic MMRd cancers. Whereas effector CD8+ T cells correlated with regression of mutational MMRd tumors, activated CD16+ NK cells were associated with ICB-responsive epigenetic MMRd tumors. These data highlight the interplay between tumor-intrinsic and tumor-extrinsic factors that influence ICB response.

SIGNIFICANCE

The molecular mechanism of MMRd is associated with response to anti-PD-1 immunotherapy in endometrial carcinoma. Tumors with epigenetic MMRd or mutational MMRd are correlated with NK cell or CD8+ T cell-driven immunity, respectively. Classifying tumors by the mechanism of MMRd may inform clinical decision-making regarding cancer immunotherapy. This article is highlighted in the In This Issue feature, p. 247.

Integrative Analysis of Germline Rare Variants in Clear and Non-Clear Cell Renal Cell Carcinoma

IMPORTANCE

RCC encompasses a set of histologically distinct cancers with a high estimated genetic heritability, of which only a portion is currently explained. Previous rare germline variant studies in RCC have usually pooled clear and non-clear cell RCCs and have not adequately accounted for population stratification that may significantly impact the interpretation and discovery of certain candidate risk genes.

OBJECTIVE

To evaluate the enrichment of germline PVs in established cancer-predisposing genes (CPGs) in clear cell and non-clear cell RCC patients compared to cancer-free controls using approaches that account for population stratification and to identify unconventional types of germline RCC risk variants that confer an increased risk of developing RCC.

DESIGN SETTING AND PARTICIPANTS

In 1,436 unselected RCC patients with sufficient data quality, we systematically identified rare germline PVs, cryptic splice variants, and copy number variants (CNVs). From this unselected cohort, 1,356 patients were ancestry-matched with 16,512 cancer-free controls, and gene-level enrichment of rare germline PVs were assessed in 143 CPGs, followed by an investigation of somatic events in matching tumor samples.

MAIN OUTCOMES AND MEASURES

Gene-level burden of rare germline PVs, identification of secondary somatic events accompanying the germline PVs, and characterization of less-explored types of rare germline PVs in RCC patients.

RESULTS

In clear cell RCC (n = 976 patients), patients exhibited significantly higher prevalence of PVs in VHL compared to controls (OR: 39.1, 95% CI: 7.01-218.07, p-value:4.95e-05, q-value:0.00584). In non-clear cell RCC (n = 380 patients), patients carried enriched burden of PVs in FH (OR: 77.9, 95% CI: 18.68-324.97, p-value:1.55e-08, q-value: 1.83e-06) and MET (OR: 1.98e11, 95% CI: 0-inf, p-value: 2.07e-05, q-value: 3.50e-07). In a CHEK2-focused analysis with European cases and controls, clear cell RCC patients (n=906 European patients) harbored nominal enrichment of the previously reported low-penetrance CHEK2 variants, p.Ile157Thr (OR:1.84, 95% CI: 1.00-3.36, p-value:0.049) and p.Ser428Phe (OR:5.20, 95% CI: 1.00-26.40, p-value:0.045) while non-clear cell RCC patients (n=295 European patients) exhibited nominal enrichment of CHEK2 LOF germline PVs (OR: 3.51, 95% CI: 1.10-11.10, p-value: 0.033). RCC patients with germline PVs in FH, MET, and VHL exhibited significantly earlier age of cancer onset compared to patients without any germline PVs in CPGs (Mean: 46.0 vs 60.2 years old, Tukey adjusted p-value < 0.0001), and more than half had secondary somatic events affecting the same gene (n=10/15, 66.7%, 95% CI: 38.7-87.0%). Conversely, patients with rare germline PVs in CHEK2 exhibited a similar age of disease onset to patients without any identified germline PVs in CPGs (Mean: 60.1 vs 60.2 years old, Tukey adjusted p-value: 0.99), and only 30.4% of the patients carried secondary somatic events in CHEK2 (n=7/23, 95% CI: 14.1-53.0%). Finally, rare pathogenic germline cryptic splice variants underexplored in RCC were identified in SDHA and TSC1, and rare pathogenic germline CNVs were found in 18 patients, including CNVs in FH, SDHA, and VHL.

CONCLUSIONS AND RELEVANCE

This systematic analysis supports the existing link between several RCC risk genes and elevated RCC risk manifesting in earlier age of RCC onset. Our analysis calls for caution when assessing the role of germline PVs in CHEK2 due to the burden of founder variants with varying population frequency in different ancestry groups. It also broadens the definition of the RCC germline landscape of pathogenicity to incorporate previously understudied types of germline variants, such as cryptic splice variants and CNVs.

Evolutionary history of transformation from chronic lymphocytic leukemia to Richter syndrome

Richter syndrome (RS) arising from chronic lymphocytic leukemia (CLL) exemplifies an aggressive malignancy that develops from an indolent neoplasm. To decipher the genetics underlying this transformation, we computationally deconvoluted admixtures of CLL and RS cells from 52 patients with RS, evaluating paired CLL-RS whole-exome sequencing data. We discovered RS-specific somatic driver mutations (including IRF2BP2, SRSF1, B2M, DNMT3A and CCND3), recurrent copy-number alterations beyond del(9p21)(CDKN2A/B), whole-genome duplication and chromothripsis, which were confirmed in 45 independent RS cases and in an external set of RS whole genomes. Through unsupervised clustering, clonally related RS was largely distinct from diffuse large B cell lymphoma. We distinguished pathways that were dysregulated in RS versus CLL, and detected clonal evolution of transformation at single-cell resolution, identifying intermediate cell states. Our study defines distinct molecular subtypes of RS and highlights cell-free DNA analysis as a potential tool for early diagnosis and monitoring.

The dual role of CD70 in B-cell lymphomagenesis

BACKGROUND

CD70 is a costimulatory molecule that is transiently expressed on a small set of activated lymphocytes and is involved in T-cell-mediated immunity. However, the role of CD70 in B-cell malignancies remains controversial.

METHODS

We investigated the clinical relevance of CD70 genetic alterations and its protein expression in two diffuse large B-cell lymphoma (DLBCL) cohorts with different ethnic backgrounds. We also performed transcriptomic analysis to explore the role of CD70 alterations in tumour microenvironment. We further tested the blockade of CD70 in combination with PD-L1 inhibitor in a murine lymphoma model.

RESULTS

We showed that CD70 genetic aberrations occurred more frequently in the Chinese DLBCL cohort (56/233, 24.0%) than in the Swedish cohort (9/84, 10.8%), especially in those with concomitant hepatitis B virus (HBV) infection. The CD70 genetic changes in DLBCL resulted in a reduction/loss of protein expression and/or CD27 binding, which might impair T cell priming and were independently associated with poor overall survival. Paradoxically, we observed that over-expression of CD70 protein was also associated with a poor treatment response, as well as an advanced disease stage and EBV infection. More exhausted CD8+ T cells were furthermore identified in CD70 high-expression DLBCLs. Finally, in a murine lymphoma model, we demonstrated that blocking the CD70/CD27 and/or PD1/PD-L1 interactions could reduce CD70+ lymphoma growth in vivo, by directly impairing the tumour cell proliferation and rescuing the exhausted T cells.

CONCLUSIONS

Our findings suggest that CD70 can play a role in either tumour suppression or oncogenesis in DLBCL, likely via distinct immune evasion mechanisms, that is, impairing T cell priming or inducing T cell exhaustion. Characterisation of specific dysfunction of CD70 in DLBCL may thus provide opportunities for the development of novel targeted immuno-therapeutic strategies.

Clinical Efficacy and Whole-Exome Sequencing of Liquid Biopsies in a Phase IB/II Study of Bazedoxifene and Palbociclib in Advanced Hormone Receptor-Positive Breast Cancer

PURPOSE

Sensitivity to endocrine therapy (ET) is critical for the clinical benefit from the combination of palbociclib plus ET in hormone receptor-positive/HER2-negative (HR+/HER2-) advanced breast cancer. Bazedoxifene is a third-generation selective estrogen receptor (ER) modulator and selective ER degrader with activity in preclinical models of endocrine-resistant breast cancer, including models harboring ESR1 mutations. Clinical trials in healthy women showed that bazedoxifene is well tolerated.

PATIENTS AND METHODS

We conducted a phase Ib/II study of bazedoxifene plus palbociclib in patients with HR+/HER2- advanced breast cancer who progressed on prior ET (N = 36; NCT02448771).

RESULTS

The study met its primary endpoint, with a clinical benefit rate of 33.3%, and the safety profile was consistent with what has previously been seen with palbociclib monotherapy. The median progression-free survival (PFS) was 3.6 months [95% confidence interval (CI), 2.0-7.2]. An activating PIK3CA mutation at baseline was associated with a shorter PFS (HR = 4.4; 95% CI, 1.5-13; P = 0.0026), but activating ESR1 mutations did not impact the PFS. Longitudinal plasma circulating tumor DNA whole-exome sequencing (WES; N = 68 plasma samples) provided an overview of the tumor heterogeneity and the subclonal genetic evolution, and identified actionable mutations acquired during treatment.

CONCLUSIONS

The combination of palbociclib and bazedoxifene has clinical efficacy and an acceptable safety profile in a heavily pretreated patient population with advanced HR+/HER2- breast cancer. These results merit continued investigation of bazedoxifene in breast cancer.

Immune biomarkers of response to immunotherapy in patients with high-risk smoldering myeloma

Patients with smoldering multiple myeloma (SMM) are observed until progression, but early treatment may improve outcomes. We conducted a phase II trial of elotuzumab, lenalidomide, and dexamethasone (EloLenDex) in patients with high-risk SMM and performed single-cell RNA and T cell receptor (TCR) sequencing on 149 bone marrow (BM) and peripheral blood (PB) samples from patients and healthy donors (HDs). We find that early treatment with EloLenDex is safe and effective and provide a comprehensive characterization of alterations in immune cell composition and TCR repertoire diversity in patients. We show that the similarity of a patient's immune cell composition to that of HDs may have prognostic relevance at diagnosis and after treatment and that the abundance of granzyme K (GZMK)+ CD8+ effector memory T (TEM) cells may be associated with treatment response. Last, we uncover similarities between immune alterations observed in the BM and PB, suggesting that PB-based immune profiling may have diagnostic and prognostic utility.

Molecular map of chronic lymphocytic leukemia and its impact on outcome

Recent advances in cancer characterization have consistently revealed marked heterogeneity, impeding the completion of integrated molecular and clinical maps for each malignancy. Here, we focus on chronic lymphocytic leukemia (CLL), a B cell neoplasm with variable natural history that is conventionally categorized into two subtypes distinguished by extent of somatic mutations in the heavy-chain variable region of immunoglobulin genes (IGHV). To build the 'CLL map,' we integrated genomic, transcriptomic and epigenomic data from 1,148 patients. We identified 202 candidate genetic drivers of CLL (109 new) and refined the characterization of IGHV subtypes, which revealed distinct genomic landscapes and leukemogenic trajectories. Discovery of new gene expression subtypes further subcategorized this neoplasm and proved to be independent prognostic factors. Clinical outcomes were associated with a combination of genetic, epigenetic and gene expression features, further advancing our prognostic paradigm. Overall, this work reveals fresh insights into CLL oncogenesis and prognostication.

A patient-driven clinicogenomic partnership for metastatic prostate cancer

Molecular profiling studies have enabled discoveries for metastatic prostate cancer (MPC) but have predominantly occurred in academic medical institutions and involved non-representative patient populations. We established the Metastatic Prostate Cancer Project (MPCproject, mpcproject.org), a patient-partnered initiative to involve patients with MPC living anywhere in the US and Canada in molecular research. Here, we present results from our partnership with the first 706 MPCproject participants. While 41% of patient partners live in rural, physician-shortage, or medically underserved areas, the MPCproject has not yet achieved racial diversity, a disparity that demands new initiatives detailed herein. Among molecular data from 333 patient partners (572 samples), exome sequencing of 63 tumor and 19 cell-free DNA (cfDNA) samples recapitulated known findings in MPC, while inexpensive ultra-low-coverage sequencing of 318 cfDNA samples revealed clinically relevant AR amplifications. This study illustrates the power of a growing, longitudinal partnership with patients to generate a more representative understanding of MPC.

Crowdis et al. describe the MPCproject (mpcproject.org), a decentralized initiative to partner with patients with metastatic prostate cancer in the US and Canada to accelerate molecular research. The authors describe clinicogenomic results from the first 706 geographically diverse patient partners and lay the foundation for sustained and inclusive partnership in this disease.

Tissue-resident memory and circulating T cells are early responders to pre-surgical cancer immunotherapy

Neoadjuvant immune checkpoint blockade has shown promising clinical activity. Here, we characterized early kinetics in tumor-infiltrating and circulating immune cells in oral cancer patients treated with neoadjuvant anti-PD-1 or anti-PD-1/CTLA-4 in a clinical trial (NCT02919683). Tumor-infiltrating CD8 T cells that clonally expanded during immunotherapy expressed elevated tissue-resident memory and cytotoxicity programs, which were already active prior to therapy, supporting the capacity for rapid response. Systematic target discovery revealed that treatment-expanded tumor T cell clones in responding patients recognized several self-antigens, including the cancer-specific antigen MAGEA1. Treatment also induced a systemic immune response characterized by expansion of activated T cells enriched for tumor-infiltrating T cell clonotypes, including both pre-existing and emergent clonotypes undetectable prior to therapy. The frequency of activated blood CD8 T cells, notably pre-treatment PD-1-positive KLRG1-negative T cells, was strongly associated with intra-tumoral pathological response. These results demonstrate how neoadjuvant checkpoint blockade induces local and systemic tumor immunity.

Genetic subtypes of smoldering multiple myeloma are associated with distinct pathogenic phenotypes and clinical outcomes

Smoldering multiple myeloma (SMM) is a precursor condition of multiple myeloma (MM) with significant heterogeneity in disease progression. Existing clinical models of progression risk do not fully capture this heterogeneity. Here we integrate 42 genetic alterations from 214 SMM patients using unsupervised binary matrix factorization (BMF) clustering and identify six distinct genetic subtypes. These subtypes are differentially associated with established MM-related RNA signatures, oncogenic and immune transcriptional profiles, and evolving clinical biomarkers. Three genetic subtypes are associated with increased risk of progression to active MM in both the primary and validation cohorts, indicating they can be used to better predict high and low-risk patients within the currently used clinical risk stratification models.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Cancer Neoantigens: Challenges and Future Directions for Prediction, Prioritization, and Validation

Prioritization of immunogenic neoantigens is key to enhancing cancer immunotherapy through the development of personalized vaccines, adoptive T cell therapy, and the prediction of response to immune checkpoint inhibition. Neoantigens are tumor-specific proteins that allow the immune system to recognize and destroy a tumor. Cancer immunotherapies, such as personalized cancer vaccines, adoptive T cell therapy, and immune checkpoint inhibition, rely on an understanding of the patient-specific neoantigen profile in order to guide personalized therapeutic strategies. Genomic approaches to predicting and prioritizing immunogenic neoantigens are rapidly expanding, raising new opportunities to advance these tools and enhance their clinical relevance. Predicting neoantigens requires acquisition of high-quality samples and sequencing data, followed by variant calling and variant annotation. Subsequently, prioritizing which of these neoantigens may elicit a tumor-specific immune response requires application and integration of tools to predict the expression, processing, binding, and recognition potentials of the neoantigen. Finally, improvement of the computational tools is held in constant tension with the availability of datasets with validated immunogenic neoantigens. The goal of this review article is to summarize the current knowledge and limitations in neoantigen prediction, prioritization, and validation and propose future directions that will improve personalized cancer treatment.

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.

Calling small variants using universality with Bayes-factor-adjusted odds ratios

The application of next-generation sequencing in research and particularly in clinical routine requires highly accurate variant calling. Here we describe UVC, a method for calling small variants of germline or somatic origin. By unifying opposite assumptions with sublation, we discovered the following two empirical laws to improve variant calling: allele fraction at high sequencing depth is inversely proportional to the cubic root of variant-calling error rate, and odds ratios adjusted with Bayes factors can model various sequencing biases. UVC outperformed other variant callers on the GIAB germline truth sets, 192 scenarios of in silico mixtures simulating 192 combinations of tumor/normal sequencing depths and tumor/normal purities, the GIAB somatic truth sets derived from physical mixture, and the SEQC2 somatic reference sets derived from the breast-cancer cell-line HCC1395. UVC achieved 100% concordance with the manual review conducted by multiple independent researchers on a Qiagen 71-gene-panel dataset derived from 16 patients with colon adenoma. UVC outperformed other unique molecular identifier (UMI)-aware variant callers on the datasets used for publishing these variant callers. Performance was measured with sensitivity-specificity trade off for called variants. The improved variant calls generated by UVC from previously published UMI-based sequencing data provided additional insight about DNA damage repair. UVC is open-sourced under the BSD 3-Clause license at https://github.com/genetronhealth/uvc and quay.io/genetronhealth/gcc-6-3-0-uvc-0-6-0-441a694.

Subtype heterogeneity and epigenetic convergence in neuroendocrine prostate cancer

Neuroendocrine carcinomas (NEC) are tumors expressing markers of neuronal differentiation that can arise at different anatomic sites but have strong histological and clinical similarities. Here we report the chromatin landscapes of a range of human NECs and show convergence to the activation of a common epigenetic program. With a particular focus on treatment emergent neuroendocrine prostate cancer (NEPC), we analyze cell lines, patient-derived xenograft (PDX) models and human clinical samples to show the existence of two distinct NEPC subtypes based on the expression of the neuronal transcription factors ASCL1 and NEUROD1. While in cell lines and PDX models these subtypes are mutually exclusive, single-cell analysis of human clinical samples exhibits a more complex tumor structure with subtypes coexisting as separate sub-populations within the same tumor. These tumor sub-populations differ genetically and epigenetically contributing to intra- and inter-tumoral heterogeneity in human metastases. Overall, our results provide a deeper understanding of the shared clinicopathological characteristics shown by NECs. Furthermore, the intratumoral heterogeneity of human NEPCs suggests the requirement of simultaneous targeting of coexisting tumor populations as a therapeutic strategy.

Neuroendocrine carcinomas (NECs) arise from different anatomic sites, but have similar histological and clinical features. Here, the authors show that the epigenetic landscape of a range of NECs converges towards a common epigenetic state, while distinct subtypes occur within neuroendocrine prostate cancer contributing to intratumor heterogeneity in clinical samples.

Discovery and Features of an Alkylating Signature in Colorectal Cancer

Several risk factors have been established for colorectal cancer, yet their direct mutagenic effects in patients' tumors remain to be elucidated. Here, we leveraged whole-exome sequencing data from 900 colorectal cancer cases that had occurred in three U.S.-wide prospective studies with extensive dietary and lifestyle information. We found an alkylating signature that was previously undescribed in colorectal cancer and then showed the existence of a similar mutational process in normal colonic crypts. This alkylating signature is associated with high intakes of processed and unprocessed red meat prior to diagnosis. In addition, this signature was more abundant in the distal colorectum, predicted to target cancer driver mutations KRAS p.G12D, KRAS p.G13D, and PIK3CA p.E545K, and associated with poor survival. Together, these results link for the first time a colorectal mutational signature to a component of diet and further implicate the role of red meat in colorectal cancer initiation and progression. SIGNIFICANCE: Colorectal cancer has several lifestyle risk factors, but the underlying mutations for most have not been observed directly in tumors. Analysis of 900 colorectal cancers with whole-exome sequencing and epidemiologic annotations revealed an alkylating mutational signature that was associated with red meat consumption and distal tumor location, as well as predicted to target KRAS p.G12D/p.G13D.This article is highlighted in the In This Issue feature, p. 2355.

Molecular correlates of response to eribulin and pembrolizumab in hormone receptor-positive metastatic breast cancer

Immune checkpoint inhibitors (ICIs) have minimal therapeutic effect in hormone receptor-positive (HR+ ) breast cancer. We present final overall survival (OS) results (n = 88) from a randomized phase 2 trial of eribulin ± pembrolizumab for patients with metastatic HR+ breast cancer, computationally dissect genomic and/or transcriptomic data from pre-treatment tumors (n = 52) for molecular associations with efficacy, and identify cytokine changes differentiating response and ICI-related toxicity (n = 58). Despite no improvement in OS with combination therapy (hazard ratio 0.95, 95% CI 0.59-1.55, p = 0.84), immune infiltration and antigen presentation distinguished responding tumors, while tumor heterogeneity and estrogen signaling independently associated with resistance. Moreover, patients with ICI-related toxicity had lower levels of immunoregulatory cytokines. Broadly, we establish a framework for ICI response in HR+ breast cancer that warrants diagnostic and therapeutic validation. ClinicalTrials.gov Registration: NCT03051659.

Differential whole-genome doubling and homologous recombination deficiencies across breast cancer subtypes from the Taiwanese population

Whole-genome doubling (WGD) is an early macro-evolutionary event in tumorigenesis, involving the doubling of an entire chromosome complement. However, its impact on breast cancer subtypes remains unclear. Here, we performed a comprehensive and quantitative analysis of WGD and its influence on breast cancer subtypes in patients from Taiwan and consequently highlight the genomic association between WGD and homologous recombination deficiency (HRD). A higher manifestation of WGD was reported in triple-negative breast cancer, conferring high chromosomal instability (CIN), while HER2 + tumors exhibited early WGD events, with widely varied CIN levels, compared to luminal-type tumors. An association of higher activity of de novo indel signature 2 with WGD and HRD in Taiwanese breast cancer patients was reported. A control test between WGD and pseudo non-WGD samples was further employed to support this finding. The study provides a better comprehension of tumorigenesis in breast cancer subtypes, thus assisting in personalized treatment.

Wu, Hsieh et al. analyze Taiwanese breast cancer patient samples using whole-exome sequencing to examine the heterogeneity and homogeneity in the timing and dependencies of somatic aberrations across disease subtypes. The authors focus on somatic alterations and related features that correlate with whole genome doubling, including homologous recombination deficiencies.

Molecular features of exceptional response to neoadjuvant anti-androgen therapy in high-risk localized prostate cancer

High-risk localized prostate cancer (HRLPC) is associated with a substantial risk of recurrence and disease mortality. Recent clinical trials have shown that intensifying anti-androgen therapies administered before prostatectomy can induce pathologic complete responses or minimal residual disease, called exceptional response, although the molecular determinants of these clinical outcomes are largely unknown. Here, we perform whole-exome and transcriptome sequencing on pre-treatment multi-regional tumor biopsies from exceptional responders (ERs) and non-responders (NRs, pathologic T3 or lymph node-positive disease) to intensive neoadjuvant anti-androgen therapies. Clonal SPOP mutation and SPOPL copy-number loss are exclusively observed in ERs, while clonal TP53 mutation and PTEN copy-number loss are exclusively observed in NRs. Transcriptional programs involving androgen signaling and TGF-β signaling are enriched in ERs and NRs, respectively. These findings may guide prospective validation studies of these molecular features in large HRLPC clinical cohorts treated with neoadjuvant anti-androgens to improve patient stratification.

Nivolumab in combination with cabozantinib for metastatic triple-negative breast cancer: a phase II and biomarker study

This single-arm phase II study investigated the efficacy and safety of cabozantinib combined with nivolumab in metastatic triple-negative breast cancer (mTNBC). The primary endpoint was objective response rate (ORR) by RECIST 1.1. Biopsies at baseline and after cycle 1 were analyzed for tumor-infiltrating lymphocytes (TILs), PD-L1, and whole-exome and transcriptome sequencing. Only 1/18 patients achieved a partial response (ORR 6%), and the trial was stopped early. Toxicity led to cabozantinib dose reduction in 50% of patients. One patient had a PD-L1-positive tumor, and three patients had TILs > 10%. The responding patient had a PD-L1-negative tumor with low tumor mutational burden but high TILs and enriched immune gene expression. High pretreatment levels of plasma immunosuppressive cytokines, chemokines, and immune checkpoint molecules were associated with rapid progression. Although this study did not meet its primary endpoint, immunostaining, genomic, and proteomic studies indicated a high degree of tumor immunosuppression in this mTNBC cohort.

Parallel Genomic Alterations of Antigen and Payload Targets Mediate Polyclonal Acquired Clinical Resistance to Sacituzumab Govitecan in Triple-Negative Breast Cancer

Sacituzumab govitecan (SG), the first antibody-drug conjugate (ADC) approved for triple-negative breast cancer, incorporates the anti-TROP2 antibody hRS7 conjugated to a topoisomerase-1 (TOP1) inhibitor payload. We sought to identify mechanisms of SG resistance through RNA and whole-exome sequencing of pretreatment and postprogression specimens. One patient exhibiting de novo progression lacked TROP2 expression, in contrast to robust TROP2 expression and focal genomic amplification of TACSTD2/TROP2 observed in a patient with a deep, prolonged response to SG. Analysis of acquired genomic resistance in this case revealed one phylogenetic branch harboring a canonical TOP1 ^E418K resistance mutation and subsequent frameshift TOP1 mutation, whereas a distinct branch exhibited a novel TACSTD2/TROP2 ^T256R missense mutation. Reconstitution experiments demonstrated that TROP2^T256R confers SG resistance via defective plasma membrane localization and reduced cell-surface binding by hRS7. These findings highlight parallel genomic alterations in both antibody and payload targets associated with resistance to SG. SIGNIFICANCE: These findings underscore TROP2 as a response determinant and reveal acquired SG resistance mechanisms involving the direct antibody and drug payload targets in distinct metastatic subclones of an individual patient. This study highlights the specificity of SG and illustrates how such mechanisms will inform therapeutic strategies to overcome ADC resistance.This article is highlighted in the In This Issue feature, p. 2355.

Phase II Multicenter Study of Enzalutamide in Metastatic Castration-Resistant Prostate Cancer to Identify Mechanisms Driving Resistance

PURPOSE

Enzalutamide is a second-generation androgen receptor (AR) inhibitor that has improved overall survival (OS) in metastatic castration-resistant prostate cancer (CRPC). However, nearly all patients develop resistance. We designed a phase II multicenter study of enzalutamide in metastatic CRPC incorporating tissue and blood biomarkers to dissect mechanisms driving resistance.

PATIENTS AND METHODS

Eligible patients with metastatic CRPC underwent a baseline metastasis biopsy and then initiated enzalutamide 160 mg daily. A repeat metastasis biopsy was obtained at radiographic progression from the same site when possible. Blood for circulating tumor cell (CTC) analysis was collected at baseline and progression. The primary objective was to analyze mechanisms of resistance in serial biopsies. Whole-exome sequencing was performed on tissue biopsies. CTC samples underwent RNA sequencing.

RESULTS

A total of 65 patients initiated treatment, of whom 22 (33.8%) had received prior abiraterone. Baseline biopsies were enriched for alterations in AR (mutations, amplifications) and tumor suppression genes (PTEN, RB1, and TP53), which were observed in 73.1% and 92.3% of baseline biopsies, respectively. Progression biopsies revealed increased AR amplifications (64.7% at progression vs. 53.9% at baseline) and BRCA2 alterations (64.7% at progression vs. 38.5% at baseline). Genomic analysis of baseline and progression CTC samples demonstrated increased AR splice variants, AR-regulated genes, and neuroendocrine markers at progression.

CONCLUSIONS

Our results demonstrate that a large proportion of enzalutamide-treated patients have baseline and progression alterations in the AR pathway and tumor suppressor genes. We demonstrate an increased number of BRCA2 alterations post-enzalutamide, highlighting the importance of serial tumor sampling in CRPC.

Radiation-related genomic profile of papillary thyroid carcinoma after the Chernobyl accident

The 1986 Chernobyl nuclear power plant accident increased papillary thyroid carcinoma (PTC) incidence in surrounding regions, particularly for radioactive iodine (¹³¹I)-exposed children. We analyzed genomic, transcriptomic, and epigenomic characteristics of 440 PTCs from Ukraine (from 359 individuals with estimated childhood ¹³¹I exposure and 81 unexposed children born after 1986). PTCs displayed radiation dose-dependent enrichment of fusion drivers, nearly all in the mitogen-activated protein kinase pathway, and increases in small deletions and simple/balanced structural variants that were clonal and bore hallmarks of nonhomologous end-joining repair. Radiation-related genomic alterations were more pronounced for individuals who were younger at exposure. Transcriptomic and epigenomic features were strongly associated with driver events but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth after environmental radiation exposure.

Evolutionary Trajectories and Genomic Divergence in Localized Breast Cancers after Ipsilateral Breast Tumor Recurrence

Simple Summary

Ipsilateral breast tumor relapse (IBTR) occurs in 5–10% of localized breast cancers (BCs) within 10 years of incidence, despite proper treatment of the primary lesion. However, the clinical consequences of evolutionary trajectories of BC cells and their impact on IBTR remain poorly understood. Here, we conducted a longitudinal genomic analysis of 10 matched localized BC patients with IBTR. Overall, we identified the differences in homologous recombination deficiency, chromosomal instability, and somatic mutation drivers between primary and relapsed lesions. Our analyses highlighted three clonal architectures that shape by distinct mutagenic processes and subclonal diversification during relapse progression. Finally, this study provided a framework, which integrated actionable biomarkers with clonal architectures, towards improvement of future treatment decisions.

Abstract

The evolutionary trajectories that drive clinical and therapeutic consequences in localized breast cancers (BCs) with ipsilateral breast tumor relapse (IBTR) remain largely unknown. Analyses of longitudinal paired whole-exome sequencing data from 10 localized BC patients with IBTR reveal that, compared to primary breast tumors, homologous recombination (HR) deficiency, inactivation of the HR pathway, chromosomal instability, and somatic driver mutations are more frequent. Furthermore, three major models of evolution in IBTR are summarized, through which relative contributions of mutational signatures shift, and the subclonal diversity expansions are shown. Optimal treatment regimens are suggested by the clinically relevant molecular features, such as HR deficiency (20%) or specific alterations (30%) with sensitivity to available FDA-approved drugs. Finally, a rationale for the development of the therapeutic management framework is provided. This study sheds light on the complicated evolution patterns in IBTR and has significant clinical implications for future improvement of treatment decisions.

A first-generation pediatric cancer dependency map

Exciting therapeutic targets are emerging from CRISPR-based screens of high mutational-burden adult cancers. A key question, however, is whether functional genomic approaches will yield new targets in pediatric cancers, known for remarkably few mutations, which often encode proteins considered challenging drug targets. To address this, we created a first-generation pediatric cancer dependency map representing 13 pediatric solid and brain tumor types. Eighty-two pediatric cancer cell lines were subjected to genome-scale CRISPR-Cas9 loss-of-function screening to identify genes required for cell survival. In contrast to the finding that pediatric cancers harbor fewer somatic mutations, we found a similar complexity of genetic dependencies in pediatric cancer cell lines compared to that in adult models. Findings from the pediatric cancer dependency map provide preclinical support for ongoing precision medicine clinical trials. The vulnerabilities observed in pediatric cancers were often distinct from those in adult cancer, indicating that repurposing adult oncology drugs will be insufficient to address childhood cancers.

Proteogenomic Landscape of Breast Cancer Tumorigenesis and Targeted Therapy

The integration of mass spectrometry-based proteomics with next-generation DNA and RNA sequencing profiles tumors more comprehensively. Here this "proteogenomics" approach was applied to 122 treatment-naive primary breast cancers accrued to preserve post-translational modifications, including protein phosphorylation and acetylation. Proteogenomics challenged standard breast cancer diagnoses, provided detailed analysis of the ERBB2 amplicon, defined tumor subsets that could benefit from immune checkpoint therapy, and allowed more accurate assessment of Rb status for prediction of CDK4/6 inhibitor responsiveness. Phosphoproteomics profiles uncovered novel associations between tumor suppressor loss and targetable kinases. Acetylproteome analysis highlighted acetylation on key nuclear proteins involved in the DNA damage response and revealed cross-talk between cytoplasmic and mitochondrial acetylation and metabolism. Our results underscore the potential of proteogenomics for clinical investigation of breast cancer through more accurate annotation of targetable pathways and biological features of this remarkably heterogeneous malignancy.

Integrated molecular drivers coordinate biological and clinical states in melanoma

We performed harmonized molecular and clinical analysis on 1,048 melanomas and discovered markedly different global genomic properties among subtypes (BRAF, (N)RAS, NF1, Triple Wild-Type), subtype-specific preferences for secondary driver genes, and active mutational processes previously unreported in melanoma. Secondary driver genes significantly enriched in specific subtypes reflected preferential dysregulation of additional pathways, such as induction of TGF-β signaling in BRAF melanomas and inactivation of the SWI/SNF complex in (N)RAS melanomas, and select co-mutation patterns coordinated selective response to immune checkpoint blockade. We also defined the mutational landscape of Triple Wild-Type melanomas and identified enrichment of DNA repair defect signatures in this subtype, which were associated with transcriptional downregulation of key DNA repair genes and may revive previously discarded or currently unconsidered therapeutic modalities for genomically stratified melanoma patient subsets. Broadly, harmonized meta-analysis of melanoma whole-exomes identified distinct molecular drivers that may point to multiple opportunities for biological and therapeutic investigation.

Genomic analyses of flow-sorted Hodgkin Reed-Sternberg cells reveal complementary mechanisms of immune evasion

Classical Hodgkin lymphoma (cHL) is composed of rare malignant Hodgkin Reed-Sternberg (HRS) cells within an extensive, but ineffective, inflammatory/immune cell infiltrate. HRS cells exhibit near-universal somatic copy gains of chromosome 9p/9p24.1, which increase expression of the programmed cell death protein 1 (PD-1) ligands. To define genetic mechanisms of response and resistance to PD-1 blockade and identify complementary treatment targets, we performed whole-exome sequencing of flow cytometry-sorted HRS cells from 23 excisional biopsies of newly diagnosed cHLs, including 8 Epstein-Barr virus-positive (EBV+) tumors. We identified significantly mutated cancer candidate genes (CCGs) as well as somatic copy number alterations and structural variations and characterized their contribution to disease-defining immune evasion mechanisms and nuclear factor κB (NF-κB), JAK/STAT, and PI3K signaling pathways. EBV- cHLs had a higher prevalence of genetic alterations in the NF-κB and major histocompatibility complex class I antigen presentation pathways. In this young cHL cohort (median age, 26 years), we identified a predominant mutational signature of spontaneous deamination of cytosine- phosphate-guanines ("Aging"), in addition to apolipoprotein B mRNA editing catalytic polypeptide-like, activation-induced cytidine deaminase, and microsatellite instability (MSI)-associated hypermutation. In particular, the mutational burden in EBV- cHLs was among the highest reported, similar to that of carcinogen-induced tumors. Together, the overall high mutational burden, MSI-associated hypermutation, and newly identified genetic alterations represent additional potential bases for the efficacy of PD-1 blockade in cHL. Of note, recurrent cHL alterations, including B2M, TNFAIP3, STAT6, GNA13, and XPO1 mutations and 2p/2p15, 6p21.32, 6q23.3, and 9p/9p24.1 copy number alterations, were also identified in >20% of primary mediastinal B-cell lymphomas, highlighting shared pathogenetic mechanisms in these diseases.

Comprehensive characteristics of somatic mutations in the normal tissues of patients with cancer and existence of somatic mutant clones linked to cancer development

BACKGROUND

Somatic mutations are a major driver of cancer development and many have now been identified in various cancer types, but the comprehensive somatic mutation status of the normal tissues matched to tumours has not been revealed.

METHOD

We analysed the somatic mutations of whole exome sequencing data in 392 patient tumour and normal tissue pairs based on the corresponding blood samples across 10 tumour types.

RESULTS

Many of the mutations involved in oncogenic pathways such as PI3K, NOTCH and TP53, were identified in the normal tissues. The ageing-related mutational signature was the most prominent contributing signature found and the mutations in the normal tissues were frequently in genes involved in late replication time (p<0.0001). Variants were rarely overlapping across tissue types but shared variants between normal and matched tumour tissue were present. These shared variants were frequently pathogenic when compared with non-shared variants (p=0.001) and showed a higher variant-allele-fraction (p<0.0001). Normal tissue-specific mutated genes were frequently non-cancer-associated (p=0.009). PIK3CA mutations were identified in 6 normal tissues and were harboured by all of the matched cancer tissues. Multiple types of PIK3CA mutations were found in normal breast and matched cancer tissues. The PIK3CA mutations exclusively present in normal tissue may indicate clonal expansions unrelated to the tumour. In addition, PIK3CA mutation was appeared that they arose before the occurrence of the allelic imbalance.

CONCLUSION

Our current results suggest that somatic mutant clones exist in normal tissues and that their clonal expansion could be linked to cancer development.

Mechanisms and therapeutic implications of hypermutation in gliomas

A high tumour mutational burden (hypermutation) is observed in some gliomas1-5; however, the mechanisms by which hypermutation develops and whether it predicts the response to immunotherapy are poorly understood. Here we comprehensively analyse the molecular determinants of mutational burden and signatures in 10,294 gliomas. We delineate two main pathways to hypermutation: a de novo pathway associated with constitutional defects in DNA polymerase and mismatch repair (MMR) genes, and a more common post-treatment pathway, associated with acquired resistance driven by MMR defects in chemotherapy-sensitive gliomas that recur after treatment with the chemotherapy drug temozolomide. Experimentally, the mutational signature of post-treatment hypermutated gliomas was recapitulated by temozolomide-induced damage in cells with MMR deficiency. MMR-deficient gliomas were characterized by a lack of prominent T cell infiltrates, extensive intratumoral heterogeneity, poor patient survival and a low rate of response to PD-1 blockade. Moreover, although bulk analyses did not detect microsatellite instability in MMR-deficient gliomas, single-cell whole-genome sequencing analysis of post-treatment hypermutated glioma cells identified microsatellite mutations. These results show that chemotherapy can drive the acquisition of hypermutated populations without promoting a response to PD-1 blockade and supports the diagnostic use of mutational burden and signatures in cancer.

Identification of cancer driver genes based on nucleotide context

Cancer genomes contain large numbers of somatic mutations, but few of these mutations drive tumor development. Current approaches identify driver genes based on mutational recurrence, or they approximate the functional consequences of nonsynonymous mutations using bioinformatic scores. While passenger mutations are enriched in characteristic nucleotide contexts, driver mutations occur in functional positions, which are not necessarily surrounded by a particular nucleotide context. We observed that mutations in contexts that deviate from the characteristic contexts around passenger mutations provide a signal in favor of driver genes. We therefore developed a method that combines this feature with the signals traditionally used for driver gene identification. We applied our method to whole-exome sequencing data from 11,873 tumor-normal pairs and identified 460 driver genes that clustered into 21 cancer-related pathways. Our study provides a resource of driver genes across 28 tumor types with additional driver genes identified based on mutations in unusual nucleotide contexts.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.