Microenvironment drives cell state, plasticity, and drug response in pancreatic cancer.

Prognostically relevant RNA expression states exist in pancreatic ductal adenocarcinoma (PDAC), but our understanding of their drivers, stability, and relationship to therapeutic response is limited. To examine these attributes systematically, we profiled metastatic biopsies and matched organoid models at single-cell resolution. In vivo, we identify a new intermediate PDAC transcriptional cell state and uncover distinct site- and state-specific tumor microenvironments (TMEs). Benchmarking models against this reference map, we reveal strong culture-specific biases in cancer cell transcriptional state representation driven by altered TME signals. We restore expression state heterogeneity by adding back in vivo-relevant factors and show plasticity in culture models. Further, we prove that non-genetic modulation of cell state can strongly influence drug responses, uncovering state-specific vulnerabilities. This work provides a broadly applicable framework for aligning cell states across in vivo and ex vivo settings, identifying drivers of transcriptional plasticity and manipulating cell state to target associated vulnerabilities.

Cancer aneuploidies are shaped primarily by effects on tumour fitness.

Aneuploidies-whole-chromosome or whole-arm imbalances-are the most prevalent alteration in cancer genomes1,2. However, it is still debated whether their prevalence is due to selection or ease of generation as passenger events1,2. Here we developed a method, BISCUT, that identifies loci subject to fitness advantages or disadvantages by interrogating length distributions of telomere- or centromere-bounded copy-number events. These loci were significantly enriched for known cancer driver genes, including genes not detected through analysis of focal copy-number events, and were often lineage specific. BISCUT identified the helicase-encoding gene WRN as a haploinsufficient tumour-suppressor gene on chromosome 8p, which is supported by several lines of evidence. We also formally quantified the role of selection and mechanical biases in driving aneuploidy, finding that rates of arm-level copy-number alterations are most highly correlated with their effects on cellular fitness1,2. These results provide insight into the driving forces behind aneuploidy and its contribution to tumorigenesis.

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.

PI3K Pathway Alterations in Peritoneal Metastases are Associated with Earlier Recurrence in Patients with Colorectal Cancer Undergoing Optimal Cytoreductive Surgery.

BACKGROUND: Colorectal cancer with peritoneal metastasis (CRC-PM) represents a biologically heterogeneous disease; yet little is known regarding the impact of tumor biology on survival outcomes following optimal cytoreductive surgery (CRS). We analyzed the frequency of alterations in cancer signaling pathways in patients with CRC-PM and their impact on recurrence-free survival (RFS) following optimal CRS. METHODS: Thirty-five consecutive CRC-PM patients who underwent optimal CRS/HIPEC and next generation sequencing of peritoneal metastases were included in the study. Alterations in eight cancer-related signaling pathways were analyzed: Wnt/APC, p53, RTK-RAS, PI3K, TGF-B, Notch, Myc, and cell cycle. The association of pathway alterations with RFS and OS following optimal cytoreduction was estimated using Cox proportional hazard modeling. RESULTS: The most frequently altered pathways were Wnt/APC (63%), p53 (63%), RTK-RAS (60%), and PI3K (23%). Among optimally cytoreduced patients with CRC-PM, PI3K pathway alterations were an independent predictor of worse RFS (hazard ratio 3.2, 95% confidence interval CI 1.3-8.3, p = 0.01) with a clinically meaningful impact on median months to recurrence (5 vs. 13 months, p = 0.02). Alterations in p53, Wnt, and RTK-RAS pathways were not significantly associated with a difference in RFS following CRS. Alterations in the four pathways were not associated with differences in OS following CRS (median OS was 50 (interquartile range 23-80) months). CONCLUSIONS: In patients with CRC-PM, PI3K pathway alterations are associated with earlier recurrence following optimal CRS, which may represent a distinct molecular subtype. This novel finding can tailor clinical trials by using PIK3CA-directed interventions to reduce risk of recurrence after optimal CRS.

Transcriptional dysregulation in developing trigeminal sensory neurons in the LgDel mouse model of DiGeorge 22q11.2 deletion syndrome.

LgDel mice, which model the heterozygous deletion of genes at human chromosome 22q11.2 associated with DiGeorge/22q11.2 deletion syndrome (22q11DS), have cranial nerve and craniofacial dysfunction as well as disrupted suckling, feeding and swallowing, similar to key 22q11DS phenotypes. Divergent trigeminal nerve (CN V) differentiation and altered trigeminal ganglion (CNgV) cellular composition prefigure these disruptions in LgDel embryos. We therefore asked whether a distinct transcriptional state in a specific population of early differentiating LgDel cranial sensory neurons, those in CNgV, a major source of innervation for appropriate oropharyngeal function, underlies this departure from typical development. LgDel versus wild-type (WT) CNgV transcriptomes differ significantly at E10.5 just after the ganglion has coalesced. Some changes parallel altered proportions of cranial placode versus cranial neural crest-derived CNgV cells. Others are consistent with a shift in anterior-posterior patterning associated with divergent LgDel cranial nerve differentiation. The most robust quantitative distinction, however, is statistically verifiable increased variability of expression levels for most of the over 17 000 genes expressed in common in LgDel versus WT CNgV. Thus, quantitative expression changes of functionally relevant genes and increased stochastic variation across the entire CNgV transcriptome at the onset of CN V differentiation prefigure subsequent disruption of cranial nerve differentiation and oropharyngeal function in LgDel mice.

Quantification of aneuploidy in targeted sequencing data using ASCETS.

SUMMARY: The expansion of targeted panel sequencing efforts has created opportunities for large-scale genomic analysis, but tools for copy-number quantification on panel data are lacking. We introduce ASCETS, a method for the efficient quantitation of arm and chromosome-level copy-number changes from targeted sequencing data. AVAILABILITY AND IMPLEMENTATION: ASCETS is implemented in R and is freely available to non-commercial users on GitHub: https://github.com/beroukhim-lab/ascets, along with detailed documentation. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

scReQTL: an approach to correlate SNVs to gene expression from individual scRNA-seq datasets.

BACKGROUND: Recently, pioneering expression quantitative trait loci (eQTL) studies on single cell RNA sequencing (scRNA-seq) data have revealed new and cell-specific regulatory single nucleotide variants (SNVs). Here, we present an alternative QTL-related approach applicable to transcribed SNV loci from scRNA-seq data: scReQTL. ScReQTL uses Variant Allele Fraction (VAFRNA) at expressed biallelic loci, and corelates it to gene expression from the corresponding cell. RESULTS: Our approach employs the advantage that, when estimated from multiple cells, VAFRNA can be used to assess effects of SNVs in a single sample or individual. In this setting scReQTL operates in the context of identical genotypes, where it is likely to capture RNA-mediated genetic interactions with cell-specific and transient effects. Applying scReQTL on scRNA-seq data generated on the 10 × Genomics Chromium platform using 26,640 mesenchymal cells derived from adipose tissue obtained from three healthy female donors, we identified 1272 unique scReQTLs. ScReQTLs common between individuals or cell types were consistent in terms of the directionality of the relationship and the effect size. Comparative assessment with eQTLs from bulk sequencing data showed that scReQTL analysis identifies a distinct set of SNV-gene correlations, that are substantially enriched in known gene-gene interactions and significant genome-wide association studies (GWAS) loci. CONCLUSION: ScReQTL is relevant to the rapidly growing source of scRNA-seq data and can be applied to outline SNVs potentially contributing to cell type-specific and/or dynamic genetic interactions from an individual scRNA-seq dataset. AVAILABILITY: https://github.com/HorvathLab/NGS/tree/master/scReQTL.

ReQTL: identifying correlations between expressed SNVs and gene expression using RNA-sequencing data.

MOTIVATION: By testing for associations between DNA genotypes and gene expression levels, expression quantitative trait locus (eQTL) analyses have been instrumental in understanding how thousands of single nucleotide variants (SNVs) may affect gene expression. As compared to DNA genotypes, RNA genetic variation represents a phenotypic trait that reflects the actual allele content of the studied system. RNA genetic variation at expressed SNV loci can be estimated using the proportion of alleles bearing the variant nucleotide (variant allele fraction, VAFRNA). VAFRNA is a continuous measure which allows for precise allele quantitation in loci where the RNA alleles do not scale with the genotype count. We describe a method to correlate VAFRNA with gene expression and assess its ability to identify genetically regulated expression solely from RNA-sequencing (RNA-seq) datasets. RESULTS: We introduce ReQTL, an eQTL modification which substitutes the DNA allele count for the variant allele fraction at expressed SNV loci in the transcriptome (VAFRNA). We exemplify the method on sets of RNA-seq data from human tissues obtained though the Genotype-Tissue Expression (GTEx) project and demonstrate that ReQTL analyses are computationally feasible and can identify a subset of expressed eQTL loci. AVAILABILITY AND IMPLEMENTATION: A toolkit to perform ReQTL analyses is available at https://github.com/HorvathLab/ReQTL. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Characterization of an autotrophic bioreactor microbial consortium degrading thiocyanate.

Thiocyanate (SCN-) forms as a by-product of cyanidation during gold ore processing and can be degraded by a variety of microorganisms utilizing it as an energy, nitrogen, sulphur and/or carbon source. In complex consortia inhabiting bioreactor systems, a range of metabolisms are sustained by SCN- degradation; however, despite the addition or presence of labile carbon sources in most bioreactor designs to date, autotrophic bacteria have been found to dominate key metabolic functions. In this study, we cultured an autotrophic SCN--degrading consortium directly from gold mine tailings. In a batch-mode bioreactor experiment, this consortium degraded 22 mM SCN-, accumulating ammonium (NH4+) and sulphate (SO42-) as the major end products. The consortium consisted of a diverse microbial community comprised of chemolithoautotrophic members, and despite the absence of an added organic carbon substrate, a significant population of heterotrophic bacteria. The role of eukaryotes in bioreactor systems is often poorly understood; however, we found their 18S rRNA genes to be most closely related to sequences from bacterivorous Amoebozoa. Through combined chemical and phylogenetic analyses, we were able to infer roles for key microbial consortium members during SCN- biodegradation. This study provides a basis for understanding the behaviour of a SCN- degrading bioreactor under autotrophic conditions, an anticipated approach to remediating SCN- at contemporary gold mines.

Genomic and Immunophenotypic Landscape of Acquired Resistance to PD-(L)1 Blockade in Non-Small-Cell Lung Cancer.

PURPOSE: Although immune checkpoint inhibitors (ICI) have extended survival in patients with non-small-cell lung cancer (NSCLC), acquired resistance (AR) to ICI frequently develops after an initial benefit. However, the mechanisms of AR to ICI in NSCLC are largely unknown. METHODS: Comprehensive tumor genomic profiling, machine learning-based assessment of tumor-infiltrating lymphocytes, multiplexed immunofluorescence, and/or HLA-I immunohistochemistry (IHC) were performed on matched pre- and post-ICI tumor biopsies from patients with NSCLC treated with ICI at the Dana-Farber Cancer Institute who developed AR to ICI. Two additional cohorts of patients with intervening chemotherapy or targeted therapies between biopsies were included as controls. RESULTS: We performed comprehensive genomic profiling and immunophenotypic characterization on samples from 82 patients with NSCLC and matched pre- and post-ICI biopsies and compared findings with a control cohort of patients with non-ICI intervening therapies between biopsies (chemotherapy, N = 32; targeted therapies, N = 89; both, N = 17). Putative resistance mutations were identified in 27.8% of immunotherapy-treated cases and included acquired loss-of-function mutations in STK11, B2M, APC, MTOR, KEAP1, and JAK1/2; these acquired alterations were not observed in the control groups. Immunophenotyping of matched pre- and post-ICI samples demonstrated significant decreases in intratumoral lymphocytes, CD3e+ and CD8a+ T cells, and PD-L1-PD1 engagement, as well as increased distance between tumor cells and CD8+PD-1+ T cells. There was a significant decrease in HLA class I expression in the immunotherapy cohort at the time of AR compared with the chemotherapy (P = .005) and the targeted therapy (P = .01) cohorts. CONCLUSION: These findings highlight the genomic and immunophenotypic heterogeneity of ICI resistance in NSCLC, which will need to be considered when developing novel therapeutic strategies aimed at overcoming resistance.

Clinicogenomic characterization of inflammatory breast cancer.

Background: Inflammatory breast cancer (IBC) is a rare and poorly characterized type of breast cancer with an aggressive clinical presentation. The biological mechanisms driving the IBC phenotype are relatively undefined-partially due to a lack of comprehensive, large-scale genomic studies and limited clinical cohorts. Patients and Methods: A retrospective analysis of 2457 patients with metastatic breast cancer who underwent targeted tumor-only DNA-sequencing was performed at Dana-Farber Cancer Institute. Clinicopathologic, single nucleotide variant (SNV), copy number variant (CNV) and tumor mutational burden (TMB) comparisons were made between clinically confirmed IBC cases within a dedicated IBC center versus non-IBC cases. Results: Clinicopathologic differences between IBC and non-IBC cases were consistent with prior reports-including IBC being associated with younger age at diagnosis, higher grade, and enrichment with hormone receptor (HR)-negative and HER2-positive tumors. The most frequent somatic alterations in IBC involved TP53 (72%), ERBB2 (32%), PIK3CA (24%), CCND1 (12%), MYC (9%), FGFR1 (8%) and GATA3 (8%). A multivariate logistic regression analysis revealed a significant enrichment in TP53 SNVs in IBC; particularly in HER2-positive and HR-positive disease which was associated with worse outcomes. Tumor mutational burden (TMB) did not differ substantially between IBC and non-IBC cases and a pathway analysis revealed an enrichment in NOTCH pathway alterations in HER2-positive disease. Conclusion: Taken together, this study provides a comprehensive, clinically informed landscape of somatic alterations in a large cohort of patients with IBC. Our data support higher frequency of TP53 mutations and a potential enrichment in NOTCH pathway activation-but overall; a lack of major genomic differences. These results both reinforce the importance of TP53 alterations in IBC pathogenesis as well as their influence on clinical outcomes; but also suggest additional analyses beyond somatic DNA-level changes are warranted.

Emergency Department Preparedness to Care for Sexual Assault Survivors: A Nationwide Study.

INTRODUCTION: Emergency departments (ED) provide trauma-informed care to sexual assault (SA) survivors and connect them with comprehensive services. Through surveying SA survivor advocates, we aimed to 1) document updated trends in the quality of care and resources offered to SA survivors and 2) identify potential disparities according to geographic regions in the US, urban vs rural clinic locations, and the availability of sexual assault nurse examiners (SANE). METHODS: We conducted a cross-sectional study between June-August 2021, surveying SA advocates who were dispatched from rape crisis centers to support survivors during ED care. Survey questions addressed two major themes in quality of care: staff preparedness to provide trauma-response care; and available resources. Staff preparedness to provide trauma-informed care was assessed through observations of staff behaviors. We used Wilcoxon rank-sum and Kruskal-Wallis tests to analyze differences in responses according to geographic regions and SANE presence. RESULTS: A total of 315 advocates from 99 crisis centers completed the survey. The survey had a participation rate of 88.7% and a completion rate of 87.9%. Advocates who indicated that a higher proportion of their cases were attended by SANEs were more likely to report higher rates of trauma-informed staff behaviors. For example, the recalled rate of staff asking patients for consent at every step of the exam was significantly associated with SANE presence (P < 0.001). With respect to access to resources, 66.7% of advocates reported that hospitals often or always have evidence collection kits available; 30.6% reported that resources such as transportation and housing are often or always available, and 55.3% reported that SANEs are often or always part of the care team. The SANEs were reported to be more frequently available in the Southwest than in other US regions (P < 0.001) and in urban as opposed to rural areas (P < 0.001). CONCLUSION: Our study indicates that support from sexual assault nurse examiners is highly associated with trauma-informed staff behaviors and comprehensive resources. Urban-rural and regional disparities exist regarding access to SANEs, suggesting that elevating nationwide quality and equity in care of survivors of sexual assault requires increased investments in SANE training and coverage.

STK11 Inactivation Predicts Rapid Recurrence in Inoperable Early-Stage Non-Small-Cell Lung Cancer.

PURPOSE: Molecular factors predicting relapse in early-stage non-small-cell lung cancer (ES-NSCLC) are poorly understood, especially in inoperable patients receiving radiotherapy (RT). In this study, we compared the genomic profiles of inoperable and operable ES-NSCLC. MATERIALS AND METHODS: This retrospective study included 53 patients with nonsquamous ES-NSCLC (stage I-II) treated at a single institution (University of Chicago) with surgery (ie, operable; n = 30) or RT (ie, inoperable; n = 23) who underwent tumor genomic profiling. A second cohort of ES-NSCLC treated with RT (Stanford, n = 39) was included to power clinical analyses. Prognostic gene alterations were identified and correlated with clinical variables. The primary clinical end point was the correlation of prognostic genes with the cumulative incidence of relapse, disease-free survival, and overall survival (OS) in a pooled RT cohort from the two institutions (N = 62). RESULTS: Although the surgery cohort exhibited lower rates of relapse, the RT cohort was highly enriched for somatic STK11 mutations (43% v 6.7%). Receiving supplemental oxygen (odds ratio [OR] = 5.5), 20+ pack-years of tobacco smoking (OR = 6.1), and Black race (OR = 4.3) were associated with increased frequency of STK11 mutations. In the pooled RT cohort (N = 62), STK11 mutation was strongly associated with inferior oncologic outcomes: 2-year incidence of relapse was 62% versus 20% and 2-year OS was 52% versus 85%, remaining independently prognostic on multivariable analyses (relapse: subdistribution hazard ratio = 4.0, P = .0041; disease-free survival: hazard ratio, 6.8, P = .0002; OS: hazard ratio, 6.0, P = .022). STK11 mutations were predominantly associated with distant failure, rather than local. CONCLUSION: In this cohort of ES-NSCLC, STK11 inactivation was associated with poor oncologic outcomes after RT and demonstrated a novel association with clinical hypoxia, which may underlie its correlation with medical inoperability. Further validation in larger cohorts and investigation of effective adjuvant systemic therapies may be warranted.

Impact of Aneuploidy and Chromosome 9p Loss on Tumor Immune Microenvironment and Immune Checkpoint Inhibitor Efficacy in NSCLC.

INTRODUCTION: Although gene-level copy number alterations have been studied as a potential biomarker of immunotherapy efficacy in NSCLC, the impact of aneuploidy burden and chromosomal arm-level events on immune checkpoint inhibitor (ICI) efficacy in NSCLC is uncertain. METHODS: Patients who received programmed cell death protein 1 or programmed death-ligand 1 (PD-L1) inhibitor at two academic centers were included. Across all 22 chromosomes analyzed, an arm was considered altered if at least 70% of its territory was either gained or deleted. Among nonsquamous NSCLCs which underwent targeted next-generation sequencing, we retrospectively quantified aneuploidy using the adjusted fraction of chromosomal arm alterations (FAA), defined as the number of altered chromosome arms divided by the number of chromosome arms assessed, adjusted for tumor purity. RESULTS: Among 2293 nonsquamous NSCLCs identified, the median FAA increased with more advanced cancer stage and decreased with higher PD-L1 tumor proportion score (TPS) levels (median FAA in TPS < 1%: 0.09, TPS 1%-49%: 0.08, TPS ≥ 50%: 0.05, p < 0.0001). There was a very weak correlation between FAA and tumor mutational burden when taken as continuous variables (R: 0.07, p = 0.0005). A total of 765 advanced nonsquamous NSCLCs with available FAA values were treated with ICIs. With decreasing FAA tertiles, there was a progressive improvement in objective response rate (ORR 15.1% in upper tertile versus 23.2% in middle tertile versus 28.4% in lowest tertile, p = 0.001), median progression-free survival (mPFS 2.5 versus 3.3 versus 4.1 mo, p < 0.0001), and median overall survival (mOS 12.5 versus 13.9 versus 16.4 mo, p = 0.006), respectively. In the arm-level enrichment analysis, chromosome 9p loss (OR = 0.22, Q = 0.0002) and chromosome 1q gain (OR = 0.43, Q = 0.002) were significantly enriched in ICI nonresponders after false discovery rate adjustment. Compared with NSCLCs without chromosome 9p loss (n = 452), those with 9p loss (n = 154) had a lower ORR (28.1% versus 7.8%, p < 0.0001), a shorter mPFS (4.1 versus 2.3 mo, p < 0.0001), and a shorter mOS (18.0 versus 9.6 mo, p < 0.0001) to immunotherapy. In addition, among NSCLCs with high PD-L1 expression (TPS ≥ 50%), chromosome 9p loss was associated with lower ORR (43% versus 6%, p < 0.0001), shorter mPFS (6.4 versus 2.6 mo, p = 0.0006), and shorter mOS (30.2 versus 14.3 mo, p = 0.0008) to immunotherapy compared with NSCLCs without 9p loss. In multivariable analysis, adjusting for key variables including FAA, chromosome 9p loss, but not 1q gain, retained a significant impact on ORR (hazard ratio [HR] = 0.25, p < 0.001), mPFS (HR = 1.49, p = 0.001), and mOS (HR = 1.47, p = 0.003). Multiplexed immunofluorescence and computational deconvolution of RNA sequencing data revealed that tumors with either high FAA levels or chromosome 9p loss had significantly fewer tumor-associated cytotoxic immune cells. CONCLUSIONS: Nonsquamous NSCLCs with high aneuploidy and chromosome 9p loss have a distinct tumor immune microenvironment and less favorable outcomes to ICIs.

Clinical Targeted Next-Generation Panel Sequencing Reveals MYC Amplification Is a Poor Prognostic Factor in Osteosarcoma.

PURPOSE: Osteosarcoma risk stratification, on the basis of the presence of metastatic disease at diagnosis and histologic response to chemotherapy, has remained unchanged for four decades, does not include genomic features, and has not facilitated treatment advances. We report on the genomic features of advanced osteosarcoma and provide evidence that genomic alterations can be used for risk stratification. MATERIALS AND METHODS: In a primary analytic patient cohort, 113 tumor and 69 normal samples from 92 patients with high-grade osteosarcoma were sequenced with OncoPanel, a targeted next-generation sequencing assay. In this primary cohort, we assessed the genomic landscape of advanced disease and evaluated the correlation between recurrent genomic events and outcome. We assessed whether prognostic associations identified in the primary cohort were maintained in a validation cohort of 86 patients with localized osteosarcoma tested with MSK-IMPACT. RESULTS: In the primary cohort, 3-year overall survival (OS) was 65%. Metastatic disease, present in 33% of patients at diagnosis, was associated with poor OS (P = .04). The most frequently altered genes in the primary cohort were TP53, RB1, MYC, CCNE1, CCND3, CDKN2A/B, and ATRX. Mutational signature 3 was present in 28% of samples. MYC amplification was associated with a worse 3-year OS in both the primary cohort (P = .015) and the validation cohort (P = .012). CONCLUSION: The most frequently occurring genomic events in advanced osteosarcoma were similar to those described in prior reports. MYC amplification, detected with clinical targeted next-generation sequencing panel tests, is associated with poorer outcomes in two independent cohorts.

Tumor aneuploidy predicts survival following immunotherapy across multiple cancers.

Tumor mutational burden (TMB) has emerged as a promising biomarker of immunotherapy response across multiple cancer types; however, clinical outcomes among patients with low TMB tumors are heterogeneous. Herein, we demonstrate that tumor aneuploidy provides independent prognostic value among patients with lower TMB (<80th percentile) tumors treated with immunotherapy. A higher aneuploidy score is associated with poor prognosis following immunotherapy among tumors with low TMB, but not those with high TMB. Importantly, aneuploidy scores can be calculated from existing clinical targeted sequencing infrastructure, facilitating deployment of aneuploidy scores as a clinical biomarker.

A proliferative subtype of colorectal liver metastases exhibits hypersensitivity to cytotoxic chemotherapy.

Personalized treatment approaches for patients with limited liver metastases from colorectal cancer are critically needed. By leveraging three large, independent cohorts of patients with colorectal liver metastases (n = 336), we found that a proliferative subtype associated with elevated CIN70 scores is linked to immune exclusion, increased metastatic proclivity, and inferior overall survival in colorectal liver metastases; however, high CIN70 scores generate a therapeutic vulnerability to DNA-damaging therapies leading to improved treatment responses. We propose CIN70 as a candidate biomarker to personalize systemic treatment options for patients with metastatic colorectal cancer. These findings are potentially broadly applicable to other human cancers.

Highly aneuploid non-small cell lung cancer shows enhanced responsiveness to concurrent radiation and immune checkpoint blockade.

Over 500 clinical trials are investigating combination radiotherapy and immune checkpoint blockade (ICB) as cancer treatments; however, the majority of trials have found no positive interaction. Here we perform a comprehensive molecular analysis of a randomized phase I clinical trial of patients with non-small cell lung cancer (NSCLC) treated with concurrent or sequential ablative radiotherapy and ICB. We show that concurrent treatment is superior to sequential treatment in augmenting local and distant tumor responses and in improving overall survival in a subset of patients with immunologically cold, highly aneuploid tumors, but not in those with less aneuploid tumors. In addition, radiotherapy alone decreases intratumoral cytotoxic T cell and adaptive immune signatures, whereas radiotherapy and ICB upregulates key immune pathways. Our findings challenge the prevailing paradigm that local ablative radiotherapy beneficially stimulates the immune response. We propose the use of tumor aneuploidy as a biomarker and therapeutic target in personalizing treatment approaches for patients with NSCLC treated with radiotherapy and ICB.