Targeting LxCxE cleft pocket of retinoblastoma protein in M2 macrophages inhibits ovarian cancer progression.

Ovarian cancer remains a major health threat with limited treatment options available. It is characterized by immunosuppressive tumor microenvironment (TME) maintained by tumor- associated macrophages (TAMs) hindering anti-tumor responses and immunotherapy efficacy. Here we show that targeting retinoblastoma protein (Rb) by disruption of its LxCxE cleft pocket, causes cell death in TAMs by induction of ER stress, p53 and mitochondria-related cell death pathways. A reduction of pro-tumor Rb high M2-type macrophages from TME in vivo enhanced T cell infiltration and inhibited cancer progression. We demonstrate an increased Rb expression in TAMs in women with ovarian cancer is associated with poorer prognosis. Ex vivo, we show analogous cell death induction by therapeutic Rb targeting in TAMs in post-surgery ascites from ovarian cancer patients. Overall, our data elucidates therapeutic targeting of the Rb LxCxE cleft pocket as a novel promising approach for ovarian cancer treatment through depletion of TAMs and re-shaping TME immune landscape. Statement of significance: Currently, targeting immunosuppressive myeloid cells in ovarian cancer microenvironment is the first priority need to enable successful immunotherapy, but no effective solutions are clinically available. We show that targeting LxCxE cleft pocket of Retinoblastoma protein unexpectedly induces preferential cell death in M2 tumor-associated macrophages. Depletion of immunosuppressive M2 tumor-associated macrophages reshapes tumor microenvironment, enhances anti-tumor T cell responses, and inhibits ovarian cancer. Thus, we identify a novel paradoxical function of Retinoblastoma protein in regulating macrophage viability as well as a promising target to enhance immunotherapy efficacy in ovarian cancer.

Multiple sclerosis patient-derived spontaneous B cells have distinct EBV and host gene expression profiles in active disease.

Epstein-Barr virus (EBV) is an aetiologic risk factor for the development of multiple sclerosis (MS). However, the role of EBV-infected B cells in the immunopathology of MS is not well understood. Here we characterized spontaneous lymphoblastoid cell lines (SLCLs) isolated from MS patients and healthy controls (HC) ex vivo to study EBV and host gene expression in the context of an individual's endogenous EBV. SLCLs derived from MS patient B cells during active disease had higher EBV lytic gene expression than SLCLs from MS patients with stable disease or HCs. Host gene expression analysis revealed activation of pathways associated with hypercytokinemia and interferon signalling in MS SLCLs and upregulation of forkhead box protein 1 (FOXP1), which contributes to EBV lytic gene expression. We demonstrate that antiviral approaches targeting EBV replication decreased cytokine production and autologous CD4+ T cell responses in this ex vivo model. These data suggest that dysregulation of intrinsic B cell control of EBV gene expression drives a pro-inflammatory, pathogenic B cell phenotype that can be attenuated by suppressing EBV lytic gene expression.

Microbial gene expression analysis of healthy and cancerous esophagus uncovers bacterial biomarkers of clinical outcomes.

Local microbiome shifts are implicated in the development and progression of gastrointestinal cancers, and in particular, esophageal carcinoma (ESCA), which is among the most aggressive malignancies. Short-read RNA sequencing (RNAseq) is currently the leading technology to study gene expression changes in cancer. However, using RNAseq to study microbial gene expression is challenging. Here, we establish a new tool to efficiently detect viral and bacterial expression in human tissues through RNAseq. This approach employs a neural network to predict reads of likely microbial origin, which are targeted for assembly into longer contigs, improving identification of microbial species and genes. This approach is applied to perform a systematic comparison of bacterial expression in ESCA and healthy esophagi. We uncover bacterial genera that are over or underabundant in ESCA vs healthy esophagi both before and after correction for possible covariates, including patient metadata. However, we find that bacterial taxonomies are not significantly associated with clinical outcomes. Strikingly, in contrast, dozens of microbial proteins were significantly associated with poor patient outcomes and in particular, proteins that perform mitochondrial functions and iron-sulfur coordination. We further demonstrate associations between these microbial proteins and dysregulated host pathways in ESCA patients. Overall, these results suggest possible influences of bacteria on the development of ESCA and uncover new prognostic biomarkers based on microbial genes. In addition, this study provides a framework for the analysis of other human malignancies whose development may be driven by pathogens.

An African-Specific Variant of TP53 Reveals PADI4 as a Regulator of p53-Mediated Tumor Suppression.

TP53 is the most frequently mutated gene in cancer, yet key target genes for p53-mediated tumor suppression remain unidentified. Here, we characterize a rare, African-specific germline variant of TP53 in the DNA-binding domain Tyr107His (Y107H). Nuclear magnetic resonance and crystal structures reveal that Y107H is structurally similar to wild-type p53. Consistent with this, we find that Y107H can suppress tumor colony formation and is impaired for the transactivation of only a small subset of p53 target genes; this includes the epigenetic modifier PADI4, which deiminates arginine to the nonnatural amino acid citrulline. Surprisingly, we show that Y107H mice develop spontaneous cancers and metastases and that Y107H shows impaired tumor suppression in two other models. We show that PADI4 is itself tumor suppressive and that it requires an intact immune system for tumor suppression. We identify a p53-PADI4 gene signature that is predictive of survival and the efficacy of immune-checkpoint inhibitors. SIGNIFICANCE: We analyze the African-centric Y107H hypomorphic variant and show that it confers increased cancer risk; we use Y107H in order to identify PADI4 as a key tumor-suppressive p53 target gene that contributes to an immune modulation signature and that is predictive of cancer survival and the success of immunotherapy. See related commentary by Bhatta and Cooks, p. 1518. This article is highlighted in the In This Issue feature, p. 1501.

Computational Methods Summarizing Mutational Patterns in Cancer: Promise and Limitations for Clinical Applications.

Since the rise of next-generation sequencing technologies, the catalogue of mutations in cancer has been continuously expanding. To address the complexity of the cancer-genomic landscape and extract meaningful insights, numerous computational approaches have been developed over the last two decades. In this review, we survey the current leading computational methods to derive intricate mutational patterns in the context of clinical relevance. We begin with mutation signatures, explaining first how mutation signatures were developed and then examining the utility of studies using mutation signatures to correlate environmental effects on the cancer genome. Next, we examine current clinical research that employs mutation signatures and discuss the potential use cases and challenges of mutation signatures in clinical decision-making. We then examine computational studies developing tools to investigate complex patterns of mutations beyond the context of mutational signatures. We survey methods to identify cancer-driver genes, from single-driver studies to pathway and network analyses. In addition, we review methods inferring complex combinations of mutations for clinical tasks and using mutations integrated with multi-omics data to better predict cancer phenotypes. We examine the use of these tools for either discovery or prediction, including prediction of tumor origin, treatment outcomes, prognosis, and cancer typing. We further discuss the main limitations preventing widespread clinical integration of computational tools for the diagnosis and treatment of cancer. We end by proposing solutions to address these challenges using recent advances in machine learning.

A deep learning approach reveals unexplored landscape of viral expression in cancer.

About 15% of human cancer cases are attributed to viral infections. To date, virus expression in tumor tissues has been mostly studied by aligning tumor RNA sequencing reads to databases of known viruses. To allow identification of divergent viruses and rapid characterization of the tumor virome, we develop viRNAtrap, an alignment-free pipeline to identify viral reads and assemble viral contigs. We utilize viRNAtrap, which is based on a deep learning model trained to discriminate viral RNAseq reads, to explore viral expression in cancers and apply it to 14 cancer types from The Cancer Genome Atlas (TCGA). Using viRNAtrap, we uncover expression of unexpected and divergent viruses that have not previously been implicated in cancer and disclose human endogenous viruses whose expression is associated with poor overall survival. The viRNAtrap pipeline provides a way forward to study viral infections associated with different clinical conditions.

The microbiome-derived metabolite TMAO drives immune activation and boosts responses to immune checkpoint blockade in pancreatic cancer.

The composition of the gut microbiome can control innate and adaptive immunity and has emerged as a key regulator of tumor growth, especially in the context of immune checkpoint blockade (ICB) therapy. However, the underlying mechanisms for how the microbiome affects tumor growth remain unclear. Pancreatic ductal adenocarcinoma (PDAC) tends to be refractory to therapy, including ICB. Using a nontargeted, liquid chromatography-tandem mass spectrometry-based metabolomic screen, we identified the gut microbe-derived metabolite trimethylamine N-oxide (TMAO), which enhanced antitumor immunity to PDAC. Delivery of TMAO intraperitoneally or via a dietary choline supplement to orthotopic PDAC-bearing mice reduced tumor growth, associated with an immunostimulatory tumor-associated macrophage (TAM) phenotype, and activated effector T cell response in the tumor microenvironment. Mechanistically, TMAO potentiated the type I interferon (IFN) pathway and conferred antitumor effects in a type I IFN-dependent manner. Delivering TMAO-primed macrophages intravenously produced similar antitumor effects. Combining TMAO with ICB (anti-PD1 and/or anti-Tim3) in a mouse model of PDAC significantly reduced tumor burden and improved survival beyond TMAO or ICB alone. Last, the levels of bacteria containing CutC (an enzyme that generates trimethylamine, the TMAO precursor) correlated with long-term survival in patients with PDAC and improved response to anti-PD1 in patients with melanoma. Together, our study identifies the gut microbial metabolite TMAO as a driver of antitumor immunity and lays the groundwork for potential therapeutic strategies targeting TMAO.

Mutated processes predict immune checkpoint inhibitor therapy benefit in metastatic melanoma.

Immune Checkpoint Inhibitor (ICI) therapy has revolutionized treatment for advanced melanoma; however, only a subset of patients benefit from this treatment. Despite considerable efforts, the Tumor Mutation Burden (TMB) is the only FDA-approved biomarker in melanoma. However, the mechanisms underlying TMB association with prolonged ICI survival are not entirely understood and may depend on numerous confounding factors. To identify more interpretable ICI response biomarkers based on tumor mutations, we train classifiers using mutations within distinct biological processes. We evaluate a variety of feature selection and classification methods and identify key mutated biological processes that provide improved predictive capability compared to the TMB. The top mutated processes we identify are leukocyte and T-cell proliferation regulation, which demonstrate stable predictive performance across different data cohorts of melanoma patients treated with ICI. This study provides biologically interpretable genomic predictors of ICI response with substantially improved predictive performance over the TMB.

Gene-expression profiles of pretreatment biopsies predict complete response of rectal cancer patients to preoperative chemoradiotherapy.

PURPOSE: Preoperative (neoadjuvant) chemoradiotherapy (CRT) and total mesorectal excision is the standard treatment for rectal cancer patients (UICC stage II/III). Up to one-third of patients treated with CRT achieve a pathological complete response (pCR). These patients could be spared from surgery and its associated morbidity and mortality, and assigned to a "watch and wait" strategy. However, reliably identifying pCR based on clinical or imaging parameters remains challenging. EXPERIMENTAL DESIGN: We generated gene-expression profiles of 175 patients with locally advanced rectal cancer enrolled in the CAO/ARO/AIO-94 and -04 trials. One hundred and sixty-one samples were used for building, training and validating a predictor of pCR using a machine learning algorithm. The performance of the classifier was validated in three independent cohorts, comprising 76 patients from (i) the CAO/ARO/AIO-94 and -04 trials (n = 14), (ii) a publicly available dataset (n = 38) and (iii) in 24 prospectively collected samples from the TransValid A trial. RESULTS: A 21-transcript signature yielded the best classification of pCR in 161 patients (Sensitivity: 0.31; AUC: 0.81), when not allowing misclassification of non-complete-responders (False-positive rate = 0). The classifier remained robust when applied to three independent datasets (n = 76). CONCLUSION: The classifier can identify >1/3 of rectal cancer patients with a pCR while never classifying patients with an incomplete response as having pCR. Importantly, we could validate this finding in three independent datasets, including a prospectively collected cohort. Therefore, this classifier could help select rectal cancer patients for a "watch and wait" strategy. TRANSLATIONAL RELEVANCE: Forgoing surgery with its associated side effects could be an option for rectal cancer patients if the prediction of a pathological complete response (pCR) after preoperative chemoradiotherapy would be possible. Based on gene-expression profiles of 161 patients a classifier was developed and validated in three independent datasets (n = 76), identifying over 1/3 of patients with pCR, while never misclassifying a non-complete-responder. Therefore, the classifier can identify patients suited for "watch and wait".

Identification of combinations of somatic mutations that predict cancer survival and immunotherapy benefit.

Cancer evolves through the accumulation of somatic mutations over time. Although several methods have been developed to characterize mutational processes in cancers, these have not been specifically designed to identify mutational patterns that predict patient prognosis. Here we present CLICnet, a method that utilizes mutational data to cluster patients by survival rate. CLICnet employs Restricted Boltzmann Machines, a type of generative neural network, which allows for the capture of complex mutational patterns associated with patient survival in different cancer types. For some cancer types, clustering produced by CLICnet also predicts benefit from anti-PD1 immune checkpoint blockade therapy, whereas for other cancer types, the mutational processes associated with survival are different from those associated with the improved anti-PD1 survival benefit. Thus, CLICnet has the ability to systematically identify and catalogue combinations of mutations that predict cancer survival, unveiling intricate associations between mutations, survival, and immunotherapy benefit.

Hard wiring of normal tissue-specific chromosome-wide gene expression levels is an additional factor driving cancer type-specific aneuploidies.

BACKGROUND: Many carcinomas have recurrent chromosomal aneuploidies specific to the tissue of tumor origin. The reason for this specificity is not completely understood. METHODS: In this study, we looked at the frequency of chromosomal arm gains and losses in different cancer types from the The Cancer Genome Atlas (TCGA) and compared them to the mean gene expression of each chromosome arm in corresponding normal tissues of origin from the Genotype-Tissue Expression (GTEx) database, in addition to the distribution of tissue-specific oncogenes and tumor suppressors on different chromosome arms. RESULTS: This analysis revealed a complex picture of factors driving tumor karyotype evolution in which some recurrent chromosomal copy number reflect the chromosome arm-wide gene expression levels of the their normal tissue of tumor origin. CONCLUSIONS: We conclude that the cancer type-specific distribution of chromosomal arm gains and losses is potentially "hardwiring" gene expression levels characteristic of the normal tissue of tumor origin, in addition to broadly modulating the expression of tissue-specific tumor driver genes.

Pyrvinium Pamoate Induces Death of Triple-Negative Breast Cancer Stem-Like Cells and Reduces Metastases through Effects on Lipid Anabolism.

Cancer stem-like cells (CSC) induce aggressive tumor phenotypes such as metastasis formation, which is associated with poor prognosis in triple-negative breast cancer (TNBC). Repurposing of FDA-approved drugs that can eradicate the CSC subcompartment in primary tumors may prevent metastatic disease, thus representing an effective strategy to improve the prognosis of TNBC. Here, we investigated spheroid-forming cells in a metastatic TNBC model. This strategy enabled us to specifically study a population of long-lived tumor cells enriched in CSCs, which show stem-like characteristics and induce metastases. To repurpose FDA-approved drugs potentially toxic for CSCs, we focused on pyrvinium pamoate (PP), an anthelmintic drug with documented anticancer activity in preclinical models. PP induced cytotoxic effects in CSCs and prevented metastasis formation. Mechanistically, the cell killing effects of PP were a result of inhibition of lipid anabolism and, more specifically, the impairment of anabolic flux from glucose to cholesterol and fatty acids. CSCs were strongly dependent upon activation of lipid biosynthetic pathways; activation of these pathways exhibited an unfavorable prognostic value in a cohort of breast cancer patients, where it predicted high probability of metastatic dissemination and tumor relapse. Overall, this work describes a new approach to target aggressive CSCs that may substantially improve clinical outcomes for patients with TNBC, who currently lack effective targeted therapeutic options. SIGNIFICANCE: These findings provide preclinical evidence that a drug repurposing approach to prevent metastatic disease in TNBC exploits lipid anabolism as a metabolic vulnerability against CSCs in primary tumors.

Single-Cell-Derived Primary Rectal Carcinoma Cell Lines Reflect Intratumor Heterogeneity Associated with Treatment Response.

PURPOSE: The standard treatment of patients with locally advanced rectal cancer consists of preoperative chemoradiotherapy (CRT) followed by surgery. However, the response of individual tumors to CRT is extremely diverse, presenting a clinical dilemma. This broad variability in treatment response is likely attributable to intratumor heterogeneity (ITH). EXPERIMENTAL DESIGN: We addressed the impact of ITH on response to CRT by establishing single-cell-derived cell lines (SCDCL) from a treatment-naïve rectal cancer biopsy after xenografting. RESULTS: Individual SCDCLs derived from the same tumor responded profoundly different to CRT in vitro. Clonal reconstruction of the tumor and derived cell lines based on whole-exome sequencing revealed nine separate clusters with distinct proportions in the SCDCLs. Missense mutations in SV2A and ZWINT were clonal in the resistant SCDCL, but not detected in the sensitive SCDCL. Single-cell genetic analysis by multiplex FISH revealed the expansion of a clone with a loss of PIK3CA in the resistant SCDCL. Gene expression profiling by tRNA-sequencing identified the activation of the Wnt, Akt, and Hedgehog signaling pathways in the resistant SCDCLs. Wnt pathway activation in the resistant SCDCLs was confirmed using a reporter assay. CONCLUSIONS: Our model system of patient-derived SCDCLs provides evidence for the critical role of ITH for treatment response in patients with rectal cancer and shows that distinct genetic aberration profiles are associated with treatment response. We identified specific pathways as the molecular basis of treatment response of individual clones, which could be targeted in resistant subclones of a heterogenous tumor.

Interplay between DNA damage repair and apoptosis shapes cancer evolution through aneuploidy and microsatellite instability.

Driver mutations and chromosomal aneuploidy are major determinants of tumorigenesis that exhibit complex relationships. Here, we identify associations between driver mutations and chromosomal aberrations that define two tumor clusters, with distinct regimes of tumor evolution underpinned by unique sets of mutations in different components of DNA damage response. Gastrointestinal and endometrial tumors comprise a separate cluster for which chromosomal-arm aneuploidy and driver mutations are mutually exclusive. The landscape of driver mutations in these tumors is dominated by mutations in DNA repair genes that are further linked to microsatellite instability. The rest of the cancer types show a positive association between driver mutations and aneuploidy, and a characteristic set of mutations that involves primarily genes for components of the apoptotic machinery. The distinct sets of mutated genes derived here show substantial prognostic power and suggest specific vulnerabilities of different cancers that might have therapeutic potential.

The GENDULF algorithm: mining transcriptomics to uncover modifier genes for monogenic diseases.

Modifier genes are believed to account for the clinical variability observed in many Mendelian disorders, but their identification remains challenging due to the limited availability of genomics data from large patient cohorts. Here, we present GENDULF (GENetic moDULators identiFication), one of the first methods to facilitate prediction of disease modifiers using healthy and diseased tissue gene expression data. GENDULF is designed for monogenic diseases in which the mechanism is loss of function leading to reduced expression of the mutated gene. When applied to cystic fibrosis, GENDULF successfully identifies multiple, previously established disease modifiers, including EHF, SLC6A14, and CLCA1. It is then utilized in spinal muscular atrophy (SMA) and predicts U2AF1 as a modifier whose low expression correlates with higher SMN2 pre-mRNA exon 7 retention. Indeed, knockdown of U2AF1 in SMA patient-derived cells leads to increased full-length SMN2 transcript and SMN protein expression. Taking advantage of the increasing availability of transcriptomic data, GENDULF is a novel addition to existing strategies for prediction of genetic disease modifiers, providing insights into disease pathogenesis and uncovering novel therapeutic targets.

Functional Genomic Complexity Defines Intratumor Heterogeneity and Tumor Aggressiveness in Liver Cancer.

Chronic inflammation and chromosome aneuploidy are major traits of primary liver cancer (PLC), which represent the second most common cause of cancer-related death worldwide. Increased cancer fitness and aggressiveness of PLC may be achieved by enhancing tumoral genomic complexity that alters tumor biology. Here, we developed a scoring method, namely functional genomic complexity (FGC), to determine the degree of molecular heterogeneity among 580 liver tumors with diverse ethnicities and etiologies by assessing integrated genomic and transcriptomic data. We found that tumors with higher FGC scores are associated with chromosome instability and TP53 mutations, and a worse prognosis, while tumors with lower FGC scores have elevated infiltrating lymphocytes and a better prognosis. These results indicate that FGC scores may serve as a surrogate to define genomic heterogeneity of PLC linked to chromosomal instability and evasion of immune surveillance. Our findings demonstrate an ability to define genomic heterogeneity and corresponding tumor biology of liver cancer based only on bulk genomic and transcriptomic data. Our data also provide a rationale for applying this approach to survey liver tumor immunity and to stratify patients for immune-based therapy.

A unique insert in the genomes of high-risk human papillomaviruses with a predicted dual role in conferring oncogenic risk.

The differences between high risk and low risk human papillomaviruses (HR-HPV and LR-HPV, respectively) that contribute to the tumorigenic potential of HR-HPV are not well understood but can be expected to involve the HPV oncoproteins, E6 and E7. We combine genome comparison and machine learning techniques to identify a previously unnoticed insert near the 3'-end of the E6 oncoprotein gene that is unique to HR-HPV. Analysis of the insert sequence suggests that it exerts a dual effect, by creating a PDZ domain-binding motif at the C-terminus of E6, as well as eliminating the overlap between the E6 and E7 coding regions in HR-HPV. We show that, as a result, the insert might enable coupled termination-reinitiation of the E6 and E7 genes, supported by motifs complementary to the human 18S rRNA. We hypothesize that the added functionality of E6 and positive regulation of E7 expression jointly account for the tumorigenic potential of HR-HPV.

In silico learning of tumor evolution through mutational time series.

Cancer arises through the accumulation of somatic mutations over time. Understanding the sequence of mutation occurrence during cancer progression can assist early and accurate diagnosis and improve clinical decision-making. Here we employ long short-term memory (LSTM) networks, a class of recurrent neural network, to learn the evolution of a tumor through an ordered sequence of mutations. We demonstrate the capacity of LSTMs to learn complex dynamics of the mutational time series governing tumor progression, allowing accurate prediction of the mutational burden and the occurrence of mutations in the sequence. Using the probabilities learned by the LSTM, we simulate mutational data and show that the simulation results are statistically indistinguishable from the empirical data. We identify passenger mutations that are significantly associated with established cancer drivers in the sequence and demonstrate that the genes carrying these mutations are substantially enriched in interactions with the corresponding driver genes. Breaking the network into modules consisting of driver genes and their interactors, we show that these interactions are associated with poor patient prognosis, thus likely conferring growth advantage for tumor progression. Thus, application of LSTM provides for prediction of numerous additional conditional drivers and reveals hitherto unknown aspects of cancer evolution.