A pan-cancer analysis of the microbiome in metastatic cancer.

Microbial communities are resident to multiple niches of the human body and are important modulators of the host immune system and responses to anticancer therapies. Recent studies have shown that complex microbial communities are present within primary tumors. To investigate the presence and relevance of the microbiome in metastases, we integrated mapping and assembly-based metagenomics, genomics, transcriptomics, and clinical data of 4,160 metastatic tumor biopsies. We identified organ-specific tropisms of microbes, enrichments of anaerobic bacteria in hypoxic tumors, associations between microbial diversity and tumor-infiltrating neutrophils, and the association of Fusobacterium with resistance to immune checkpoint blockade (ICB) in lung cancer. Furthermore, longitudinal tumor sampling revealed temporal evolution of the microbial communities and identified bacteria depleted upon ICB. Together, we generated a pan-cancer resource of the metastatic tumor microbiome that may contribute to advancing treatment strategies.

Identifying Epistasis in Cancer Genomes: A Delicate Affair.

Recent studies of the tumor genome seek to identify cancer pathways as groups of genes in which mutations are epistatic with one another or, specifically, "mutually exclusive." Here, we show that most mutations are mutually exclusive not due to pathway structure but to interactions with disease subtype and tumor mutation load. In particular, many cancer driver genes are mutated preferentially in tumors with few mutations overall, causing mutations in these cancer genes to appear mutually exclusive with numerous others. Researchers should view current epistasis maps with caution until we better understand the multiple cause-and-effect relationships among factors such as tumor subtype, positive selection for mutations, and gross tumor characteristics including mutational signatures and load.

A Living Biobank of Breast Cancer Organoids Captures Disease Heterogeneity.

Breast cancer (BC) comprises multiple distinct subtypes that differ genetically, pathologically, and clinically. Here, we describe a robust protocol for long-term culturing of human mammary epithelial organoids. Using this protocol, >100 primary and metastatic BC organoid lines were generated, broadly recapitulating the diversity of the disease. BC organoid morphologies typically matched the histopathology, hormone receptor status, and HER2 status of the original tumor. DNA copy number variations as well as sequence changes were consistent within tumor-organoid pairs and largely retained even after extended passaging. BC organoids furthermore populated all major gene-expression-based classification groups and allowed in vitro drug screens that were consistent with in vivo xeno-transplantations and patient response. This study describes a representative collection of well-characterized BC organoids available for cancer research and drug development, as well as a strategy to assess in vitro drug response in a personalized fashion.

A Vulnerability of a Subset of Colon Cancers with Potential Clinical Utility.

BRAF(V600E) mutant colon cancers (CCs) have a characteristic gene expression signature that is also found in some tumors lacking this mutation. Collectively, they are referred to as "BRAF-like" tumors and represent some 20% of CCs. We used a shRNA-based genetic screen focused on genes upregulated in BRAF(V600E) CCs to identify vulnerabilities of this tumor subtype that might be exploited therapeutically. Here, we identify RANBP2 (also known as NUP358) as essential for survival of BRAF-like, but not for non-BRAF-like, CC cells. Suppression of RANBP2 results in mitotic defects only in BRAF-like CC cells, leading to cell death. Mechanistically, RANBP2 silencing reduces microtubule outgrowth from the kinetochores, thereby inducing spindle perturbations, providing an explanation for the observed mitotic defects. We find that BRAF-like CCs display far greater sensitivity to the microtubule poison vinorelbine both in vitro and in vivo, suggesting that vinorelbine is a potential tailored treatment for BRAF-like CCs.

A Landscape of Pharmacogenomic Interactions in Cancer.

Systematic studies of cancer genomes have provided unprecedented insights into the molecular nature of cancer. Using this information to guide the development and application of therapies in the clinic is challenging. Here, we report how cancer-driven alterations identified in 11,289 tumors from 29 tissues (integrating somatic mutations, copy number alterations, DNA methylation, and gene expression) can be mapped onto 1,001 molecularly annotated human cancer cell lines and correlated with sensitivity to 265 drugs. We find that cell lines faithfully recapitulate oncogenic alterations identified in tumors, find that many of these associate with drug sensitivity/resistance, and highlight the importance of tissue lineage in mediating drug response. Logic-based modeling uncovers combinations of alterations that sensitize to drugs, while machine learning demonstrates the relative importance of different data types in predicting drug response. Our analysis and datasets are rich resources to link genotypes with cellular phenotypes and to identify therapeutic options for selected cancer sub-populations.

Prospective derivation of a living organoid biobank of colorectal cancer patients.

In Rspondin-based 3D cultures, Lgr5 stem cells from multiple organs form ever-expanding epithelial organoids that retain their tissue identity. We report the establishment of tumor organoid cultures from 20 consecutive colorectal carcinoma (CRC) patients. For most, organoids were also generated from adjacent normal tissue. Organoids closely recapitulate several properties of the original tumor. The spectrum of genetic changes within the "living biobank" agrees well with previous large-scale mutational analyses of CRC. Gene expression analysis indicates that the major CRC molecular subtypes are represented. Tumor organoids are amenable to high-throughput drug screens allowing detection of gene-drug associations. As an example, a single organoid culture was exquisitely sensitive to Wnt secretion (porcupine) inhibitors and carried a mutation in the negative Wnt feedback regulator RNF43, rather than in APC. Organoid technology may fill the gap between cancer genetics and patient trials, complement cell-line- and xenograft-based drug studies, and allow personalized therapy design. PAPERCLIP.

γδ T cells are effectors of immunotherapy in cancers with HLA class I defects.

DNA mismatch repair-deficient (MMR-d) cancers present an abundance of neoantigens that is thought to explain their exceptional responsiveness to immune checkpoint blockade (ICB)1,2. Here, in contrast to other cancer types3-5, we observed that 20 out of 21 (95%) MMR-d cancers with genomic inactivation of ß2-microglobulin (encoded by B2M) retained responsiveness to ICB, suggesting the involvement of immune effector cells other than CD8+ T cells in this context. We next identified a strong association between B2M inactivation and increased infiltration by γδ T cells in MMR-d cancers. These γδ T cells mainly comprised the Vδ1 and Vδ3 subsets, and expressed high levels of PD-1, other activation markers, including cytotoxic molecules, and a broad repertoire of killer-cell immunoglobulin-like receptors. In vitro, PD-1+ γδ T cells that were isolated from MMR-d colon cancers exhibited enhanced reactivity to human leukocyte antigen (HLA)-class-I-negative MMR-d colon cancer cell lines and B2M-knockout patient-derived tumour organoids compared with antigen-presentation-proficient cells. By comparing paired tumour samples from patients with MMR-d colon cancer that were obtained before and after dual PD-1 and CTLA-4 blockade, we found that immune checkpoint blockade substantially increased the frequency of γδ T cells in B2M-deficient cancers. Taken together, these data indicate that γδ T cells contribute to the response to immune checkpoint blockade in patients with HLA-class-I-negative MMR-d colon cancers, and underline the potential of γδ T cells in cancer immunotherapy.

Chromatin position effects assayed by thousands of reporters integrated in parallel.

Reporter genes integrated into the genome are a powerful tool to reveal effects of regulatory elements and local chromatin context on gene expression. However, so far such reporter assays have been of low throughput. Here, we describe a multiplexing approach for the parallel monitoring of transcriptional activity of thousands of randomly integrated reporters. More than 27,000 distinct reporter integrations in mouse embryonic stem cells, obtained with two different promoters, show ∼1,000-fold variation in expression levels. Data analysis indicates that lamina-associated domains act as attenuators of transcription, likely by reducing access of transcription factors to binding sites. Furthermore, chromatin compaction is predictive of reporter activity. We also found evidence for crosstalk between neighboring genes and estimate that enhancers can influence gene expression on average over ∼20 kb. The multiplexed reporter assay is highly flexible in design and can be modified to query a wide range of aspects of gene regulation.

Effective drug combinations in breast, colon and pancreatic cancer cells.

Combinations of anti-cancer drugs can overcome resistance and provide new treatments1,2. The number of possible drug combinations vastly exceeds what could be tested clinically. Efforts to systematically identify active combinations and the tissues and molecular contexts in which they are most effective could accelerate the development of combination treatments. Here we evaluate the potency and efficacy of 2,025 clinically relevant two-drug combinations, generating a dataset encompassing 125 molecularly characterized breast, colorectal and pancreatic cancer cell lines. We show that synergy between drugs is rare and highly context-dependent, and that combinations of targeted agents are most likely to be synergistic. We incorporate multi-omic molecular features to identify combination biomarkers and specify synergistic drug combinations and their active contexts, including in basal-like breast cancer, and microsatellite-stable or KRAS-mutant colon cancer. Our results show that irinotecan and CHEK1 inhibition have synergistic effects in microsatellite-stable or KRAS-TP53 double-mutant colon cancer cells, leading to apoptosis and suppression of tumour xenograft growth. This study identifies clinically relevant effective drug combinations in distinct molecular subpopulations and is a resource to guide rational efforts to develop combinatorial drug treatments.

Truncated FGFR2 is a clinically actionable oncogene in multiple cancers.

Somatic hotspot mutations and structural amplifications and fusions that affect fibroblast growth factor receptor 2 (encoded by FGFR2) occur in multiple types of cancer1. However, clinical responses to FGFR inhibitors have remained variable1-9, emphasizing the need to better understand which FGFR2 alterations are oncogenic and therapeutically targetable. Here we apply transposon-based screening10,11 and tumour modelling in mice12,13, and find that the truncation of exon 18 (E18) of Fgfr2 is a potent driver mutation. Human oncogenomic datasets revealed a diverse set of FGFR2 alterations, including rearrangements, E1-E17 partial amplifications, and E18 nonsense and frameshift mutations, each causing the transcription of E18-truncated FGFR2 (FGFR2ΔE18). Functional in vitro and in vivo examination of a compendium of FGFR2ΔE18 and full-length variants pinpointed FGFR2-E18 truncation as single-driver alteration in cancer. By contrast, the oncogenic competence of FGFR2 full-length amplifications depended on a distinct landscape of cooperating driver genes. This suggests that genomic alterations that generate stable FGFR2ΔE18 variants are actionable therapeutic targets, which we confirmed in preclinical mouse and human tumour models, and in a clinical trial. We propose that cancers containing any FGFR2 variant with a truncated E18 should be considered for FGFR-targeted therapies.

Neoadjuvant Immunotherapy in Locally Advanced Mismatch Repair-Deficient Colon Cancer.

BACKGROUND: Mismatch repair-deficient (dMMR) tumors can be found in 10 to 15% of patients with nonmetastatic colon cancer. In these patients, the efficacy of chemotherapy is limited. The use of neoadjuvant immunotherapy has shown promising results, but data from studies of this approach are limited. METHODS: We conducted a phase 2 study in which patients with nonmetastatic, locally advanced, previously untreated dMMR colon cancer were treated with neoadjuvant nivolumab plus ipilimumab. The two primary end points were safety, defined by timely surgery (i.e., ≤2-week delay of planned surgery owing to treatment-related toxic events), and 3-year disease-free survival. Secondary end points included pathological response and results of genomic analyses. RESULTS: Of 115 enrolled patients, 113 (98%; 97.5% confidence interval [CI], 93 to 100) underwent timely surgery; 2 patients had surgery delayed by more than 2 weeks. Grade 3 or 4 immune-related adverse events occurred in 5 patients (4%), and none of the patients discontinued treatment because of adverse events. Among the 111 patients included in the efficacy analysis, a pathological response was observed in 109 (98%; 95% CI, 94 to 100), including 105 (95%) with a major pathological response (defined as ≤10% residual viable tumor) and 75 (68%) with a pathological complete response (0% residual viable tumor). With a median follow-up of 26 months (range, 9 to 65), no patients have had recurrence of disease. CONCLUSIONS: In patients with locally advanced dMMR colon cancer, neoadjuvant nivolumab plus ipilimumab had an acceptable safety profile and led to a pathological response in a high proportion of patients. (Funded by Bristol Myers Squibb; NICHE-2 ClinicalTrials.gov number, NCT03026140.).

MED12 controls the response to multiple cancer drugs through regulation of TGF-ß receptor signaling.

Inhibitors of the ALK and EGF receptor tyrosine kinases provoke dramatic but short-lived responses in lung cancers harboring EML4-ALK translocations or activating mutations of EGFR, respectively. We used a large-scale RNAi screen to identify MED12, a component of the transcriptional MEDIATOR complex that is mutated in cancers, as a determinant of response to ALK and EGFR inhibitors. MED12 is in part cytoplasmic where it negatively regulates TGF-ßR2 through physical interaction. MED12 suppression therefore results in activation of TGF-ßR signaling, which is both necessary and sufficient for drug resistance. TGF-ß signaling causes MEK/ERK activation, and consequently MED12 suppression also confers resistance to MEK and BRAF inhibitors in other cancers. MED12 loss induces an EMT-like phenotype, which is associated with chemotherapy resistance in colon cancer patients and to gefitinib in lung cancer. Inhibition of TGF-ßR signaling restores drug responsiveness in MED12(KD) cells, suggesting a strategy to treat drug-resistant tumors that have lost MED12.

Ductal carcinoma in situ develops within clonal fields of mutant cells in morphologically normal ducts.

Mutations are abundantly present in tissues of healthy individuals, including the breast epithelium. Yet it remains unknown whether mutant cells directly induce lesion formation or first spread, leading to a field of mutant cells that is predisposed towards lesion formation. To study the clonal and spatial relationships between morphologically normal breast epithelium adjacent to pre-cancerous lesions, we developed a three-dimensional (3D) imaging pipeline combined with spatially resolved genomics on archival, formalin-fixed breast tissue with the non-obligate breast cancer precursor ductal carcinoma in situ (DCIS). Using this 3D image-guided characterization method, we built high-resolution spatial maps of DNA copy number aberration (CNA) profiles within the DCIS lesion and the surrounding normal mammary ducts. We show that the local heterogeneity within a DCIS lesion is limited. However, by mapping the CNA profiles back onto the 3D reconstructed ductal subtree, we find that in eight out of 16 cases the healthy epithelium adjacent to the DCIS lesions has overlapping structural variations with the CNA profile of the DCIS. Together, our study indicates that pre-malignant breast transformations frequently develop within mutant clonal fields of morphologically normal-looking ducts. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Loss of p53 triggers WNT-dependent systemic inflammation to drive breast cancer metastasis.

Cancer-associated systemic inflammation is strongly linked to poor disease outcome in patients with cancer1,2. For most human epithelial tumour types, high systemic neutrophil-to-lymphocyte ratios are associated with poor overall survival3, and experimental studies have demonstrated a causal relationship between neutrophils and metastasis4,5. However, the cancer-cell-intrinsic mechanisms that dictate the substantial heterogeneity in systemic neutrophilic inflammation between tumour-bearing hosts are largely unresolved. Here, using a panel of 16 distinct genetically engineered mouse models for breast cancer, we uncover a role for cancer-cell-intrinsic p53 as a key regulator of pro-metastatic neutrophils. Mechanistically, loss of p53 in cancer cells induced the secretion of WNT ligands that stimulate tumour-associated macrophages to produce IL-1ß, thus driving systemic inflammation. Pharmacological and genetic blockade of WNT secretion in p53-null cancer cells reverses macrophage production of IL-1ß and subsequent neutrophilic inflammation, resulting in reduced metastasis formation. Collectively, we demonstrate a mechanistic link between the loss of p53 in cancer cells, secretion of WNT ligands and systemic neutrophilia that potentiates metastatic progression. These insights illustrate the importance of the genetic makeup of breast tumours in dictating pro-metastatic systemic inflammation, and set the stage for personalized immune intervention strategies for patients with cancer.

In situ CRISPR-Cas9 base editing for the development of genetically engineered mouse models of breast cancer.

Genetically engineered mouse models (GEMMs) of cancer have proven to be of great value for basic and translational research. Although CRISPR-based gene disruption offers a fast-track approach for perturbing gene function and circumvents certain limitations of standard GEMM development, it does not provide a flexible platform for recapitulating clinically relevant missense mutations in vivo. To this end, we generated knock-in mice with Cre-conditional expression of a cytidine base editor and tested their utility for precise somatic engineering of missense mutations in key cancer drivers. Upon intraductal delivery of sgRNA-encoding vectors, we could install point mutations with high efficiency in one or multiple endogenous genes in situ and assess the effect of defined allelic variants on mammary tumorigenesis. While the system also produces bystander insertions and deletions that can stochastically be selected for when targeting a tumor suppressor gene, we could effectively recapitulate oncogenic nonsense mutations. We successfully applied this system in a model of triple-negative breast cancer, providing the proof of concept for extending this flexible somatic base editing platform to other tissues and tumor types.

Evaluating the effectiveness of pre-operative diagnosis of ovarian cancer using minimally invasive liquid biopsies by combining serum human epididymis protein 4 and cell-free DNA in patients with an ovarian mass.

OBJECTIVE: To assess the feasibility of scalable, objective, and minimally invasive liquid biopsy-derived biomarkers such as cell-free DNA copy number profiles, human epididymis protein 4 (HE4), and cancer antigen 125 (CA125) for pre-operative risk assessment of early-stage ovarian cancer in a clinically representative and diagnostically challenging population and to compare the performance of these biomarkers with the Risk of Malignancy Index (RMI). METHODS: In this case-control study, we included 100 patients with an ovarian mass clinically suspected to be early-stage ovarian cancer. Of these 100 patients, 50 were confirmed to have a malignant mass (cases) and 50 had a benign mass (controls). Using WisecondorX, an algorithm used extensively in non-invasive prenatal testing, we calculated the benign-calibrated copy number profile abnormality score. This score represents how different a sample is from benign controls based on copy number profiles. We combined this score with HE4 serum concentration to separate cases and controls. RESULTS: Combining the benign-calibrated copy number profile abnormality score with HE4, we obtained a model with a significantly higher sensitivity (42% vs 0%; p<0.002) at 99% specificity as compared with the RMI that is currently employed in clinical practice. Investigating performance in subgroups, we observed especially large differences in the advanced stage and non-high-grade serous ovarian cancer groups. CONCLUSION: This study demonstrates that cell-free DNA can be successfully employed to perform pre-operative risk of malignancy assessment for ovarian masses; however, results warrant validation in a more extensive clinical study.

Predicting patient response with models trained on cell lines and patient-derived xenografts by nonlinear transfer learning.

Preclinical models have been the workhorse of cancer research, producing massive amounts of drug response data. Unfortunately, translating response biomarkers derived from these datasets to human tumors has proven to be particularly challenging. To address this challenge, we developed TRANSACT, a computational framework that builds a consensus space to capture biological processes common to preclinical models and human tumors and exploits this space to construct drug response predictors that robustly transfer from preclinical models to human tumors. TRANSACT performs favorably compared to four competing approaches, including two deep learning approaches, on a set of 23 drug prediction challenges on The Cancer Genome Atlas and 226 metastatic tumors from the Hartwig Medical Foundation. We demonstrate that response predictions deliver a robust performance for a number of therapies of high clinical importance: platinum-based chemotherapies, gemcitabine, and paclitaxel. In contrast to other approaches, we demonstrate the interpretability of the TRANSACT predictors by correctly identifying known biomarkers of targeted therapies, and we propose potential mechanisms that mediate the resistance to two chemotherapeutic agents.

ecpc: an R-package for generic co-data models for high-dimensional prediction.

BACKGROUND: High-dimensional prediction considers data with more variables than samples. Generic research goals are to find the best predictor or to select variables. Results may be improved by exploiting prior information in the form of co-data, providing complementary data not on the samples, but on the variables. We consider adaptive ridge penalised generalised linear and Cox models, in which the variable-specific ridge penalties are adapted to the co-data to give a priori more weight to more important variables. The R-package ecpc originally accommodated various and possibly multiple co-data sources, including categorical co-data, i.e. groups of variables, and continuous co-data. Continuous co-data, however, were handled by adaptive discretisation, potentially inefficiently modelling and losing information. As continuous co-data such as external p values or correlations often arise in practice, more generic co-data models are needed. RESULTS: Here, we present an extension to the method and software for generic co-data models, particularly for continuous co-data. At the basis lies a classical linear regression model, regressing prior variance weights on the co-data. Co-data variables are then estimated with empirical Bayes moment estimation. After placing the estimation procedure in the classical regression framework, extension to generalised additive and shape constrained co-data models is straightforward. Besides, we show how ridge penalties may be transformed to elastic net penalties. In simulation studies we first compare various co-data models for continuous co-data from the extension to the original method. Secondly, we compare variable selection performance to other variable selection methods. The extension is faster than the original method and shows improved prediction and variable selection performance for non-linear co-data relations. Moreover, we demonstrate use of the package in several genomics examples throughout the paper. CONCLUSIONS: The R-package ecpc accommodates linear, generalised additive and shape constrained additive co-data models for the purpose of improved high-dimensional prediction and variable selection. The extended version of the package as presented here (version number 3.1.1 and higher) is available on ( https://cran.r-project.org/web/packages/ecpc/ ).

Tumour-educated platelets for breast cancer detection: biological and technical insights.

BACKGROUND: Studies have shown that blood platelets contain tumour-specific mRNA profiles tumour-educated platelets (TEPs). Here, we aim to train a TEP-based breast cancer detection classifier. METHODS: Platelet mRNA was sequenced from 266 women with stage I-IV breast cancer and 212 female controls from 6 hospitals. A particle swarm optimised support vector machine (PSO-SVM) and an elastic net-based classifier (EN) were trained on 71% of the study population. Classifier performance was evaluated in the remainder (29%) of the population, followed by validation in an independent set (37 cases and 36 controls). Potential confounding was assessed in post hoc analyses. RESULTS: Both classifiers reached an area under the curve (AUC) of 0.85 upon internal validation. Reproducibility in the independent validation set was poor with an AUC of 0.55 and 0.54 for the PSO-SVM and EN classifier, respectively. Post hoc analyses indicated that 19% of the variance in gene expression was associated with hospital. Genes related to platelet activity were differentially expressed between hospitals. CONCLUSIONS: We could not validate two TEP-based breast cancer classifiers in an independent validation cohort. The TEP protocol is sensitive to within-protocol variation and revision might be necessary before TEPs can be reconsidered for breast cancer detection.

Large-scale mutagenesis in p19(ARF)- and p53-deficient mice identifies cancer genes and their collaborative networks.

p53 and p19(ARF) are tumor suppressors frequently mutated in human tumors. In a high-throughput screen in mice for mutations collaborating with either p53 or p19(ARF) deficiency, we identified 10,806 retroviral insertion sites, implicating over 300 loci in tumorigenesis. This dataset reveals 20 genes that are specifically mutated in either p19(ARF)-deficient, p53-deficient or wild-type mice (including Flt3, mmu-mir-106a-363, Smg6, and Ccnd3), as well as networks of significant collaborative and mutually exclusive interactions between cancer genes. Furthermore, we found candidate tumor suppressor genes, as well as distinct clusters of insertions within genes like Flt3 and Notch1 that induce mutants with different spectra of genetic interactions. Cross species comparative analysis with aCGH data of human cancer cell lines revealed known and candidate oncogenes (Mmp13, Slamf6, and Rreb1) and tumor suppressors (Wwox and Arfrp2). This dataset should prove to be a rich resource for the study of genetic interactions that underlie tumorigenesis.