Epigenetic Control of Cancer Cell Dormancy and Awakening in Endocrine Therapy Resistance.

SUMMARY: Rosano, Sofyali, Dhiman, and colleagues show that epigenetic-related changes occur in endocrine therapy (ET)-induced dormancy in estrogen receptor positive (ER+) breast cancer, as well as in its reawakening. Targeting these epigenetic changes blocks the entrance to dormancy and reduces the persister cancer cell population, enhancing the cytotoxic effects of ET in vitro. See related article by Rosano et al., p. 866 (9).

Epigenetic Control of Cancer Cell Dormancy and Awakening in Endocrine Therapy Resistance.

SUMMARY: Rosano, Sofyali, Dhiman, and colleagues show that epigenetic-related changes occur in endocrine therapy (ET)-induced dormancy in estrogen receptor positive (ER+) breast cancer, as well as in its reawakening. Targeting these epigenetic changes blocks the entrance to dormancy and reduces the persister cancer cell population, enhancing the cytotoxic effects of ET in vitro. See related article by Rosano et al. (9).

Functional screening of amplification outlier oncogenes in organoid models of early tumorigenesis.

Somatic copy number gains are pervasive across cancer types, yet their roles in oncogenesis are insufficiently evaluated. This inadequacy is partly due to copy gains spanning large chromosomal regions, obscuring causal loci. Here, we employed organoid modeling to evaluate candidate oncogenic loci identified via integrative computational analysis of extreme copy gains overlapping with extreme expression dysregulation in The Cancer Genome Atlas. Subsets of "outlier" candidates were contextually screened as tissue-specific cDNA lentiviral libraries within cognate esophagus, oral cavity, colon, stomach, pancreas, and lung organoids bearing initial oncogenic mutations. Iterative analysis nominated the kinase DYRK2 at 12q15 as an amplified head and neck squamous carcinoma oncogene in p53-/- oral mucosal organoids. Similarly, FGF3, amplified at 11q13 in 41% of esophageal squamous carcinomas, promoted p53-/- esophageal organoid growth reversible by small molecule and soluble receptor antagonism of FGFRs. Our studies establish organoid-based contextual screening of candidate genomic drivers, enabling functional evaluation during early tumorigenesis.

p53 governs an AT1 differentiation programme in lung cancer suppression.

Lung cancer is the leading cause of cancer deaths worldwide1. Mutations in the tumour suppressor gene TP53 occur in 50% of lung adenocarcinomas (LUADs) and are linked to poor prognosis1-4, but how p53 suppresses LUAD development remains enigmatic. We show here that p53 suppresses LUAD by governing cell state, specifically by promoting alveolar type 1 (AT1) differentiation. Using mice that express oncogenic Kras and null, wild-type or hypermorphic Trp53 alleles in alveolar type 2 (AT2) cells, we observed graded effects of p53 on LUAD initiation and progression. RNA sequencing and ATAC sequencing of LUAD cells uncovered a p53-induced AT1 differentiation programme during tumour suppression in vivo through direct DNA binding, chromatin remodelling and induction of genes characteristic of AT1 cells. Single-cell transcriptomics analyses revealed that during LUAD evolution, p53 promotes AT1 differentiation through action in a transitional cell state analogous to a transient intermediary seen during AT2-to-AT1 cell differentiation in alveolar injury repair. Notably, p53 inactivation results in the inappropriate persistence of these transitional cancer cells accompanied by upregulated growth signalling and divergence from lung lineage identity, characteristics associated with LUAD progression. Analysis of Trp53 wild-type and Trp53-null mice showed that p53 also directs alveolar regeneration after injury by regulating AT2 cell self-renewal and promoting transitional cell differentiation into AT1 cells. Collectively, these findings illuminate mechanisms of p53-mediated LUAD suppression, in which p53 governs alveolar differentiation, and suggest that tumour suppression reflects a fundamental role of p53 in orchestrating tissue repair after injury.

Molecular classification and biomarkers of clinical outcome in breast ductal carcinoma in situ: Analysis of TBCRC 038 and RAHBT cohorts.

Ductal carcinoma in situ (DCIS) is the most common precursor of invasive breast cancer (IBC), with variable propensity for progression. We perform multiscale, integrated molecular profiling of DCIS with clinical outcomes by analyzing 774 DCIS samples from 542 patients with 7.3 years median follow-up from the Translational Breast Cancer Research Consortium 038 study and the Resource of Archival Breast Tissue cohorts. We identify 812 genes associated with ipsilateral recurrence within 5 years from treatment and develop a classifier that predicts DCIS or IBC recurrence in both cohorts. Pathways associated with recurrence include proliferation, immune response, and metabolism. Distinct stromal expression patterns and immune cell compositions are identified. Our multiscale approach employed in situ methods to generate a spatially resolved atlas of breast precancers, where complementary modalities can be directly compared and correlated with conventional pathology findings, disease states, and clinical outcome.

ZFP281 drives a mesenchymal-like dormancy program in early disseminated breast cancer cells that prevents metastatic outgrowth in the lung.

Increasing evidence shows that cancer cells can disseminate from early evolved primary lesions much earlier than the classical metastasis models predicted. Here, we reveal at a single-cell resolution that mesenchymal-like (M-like) and pluripotency-like programs coordinate dissemination and a long-lived dormancy program of early disseminated cancer cells (DCCs). The transcription factor ZFP281 induces a permissive state for heterogeneous M-like transcriptional programs, which associate with a dormancy signature and phenotype in vivo. Downregulation of ZFP281 leads to a loss of an invasive, M-like dormancy phenotype and a switch to lung metastatic outgrowth. We also show that FGF2 and TWIST1 induce ZFP281 expression to induce the M-like state, which is linked to CDH1 downregulation and upregulation of CDH11. We found that ZFP281 not only controls the early dissemination of cancer cells but also locks early DCCs in a dormant state by preventing the acquisition of an epithelial-like proliferative program and consequent metastases outgrowth.

Machine Learning Based Microbiome Signature to Predict Inflammatory Bowel Disease Subtypes.

Inflammatory bowel disease (IBD) is a chronic disease with unknown pathophysiological mechanisms. There is evidence of the role of microorganims in this disease development. Thanks to the open access to multiple omics data, it is possible to develop predictive models that are able to prognosticate the course and development of the disease. The interpretability of these models, and the study of the variables used, allows the identification of biological aspects of great importance in the development of the disease. In this work we generated a metagenomic signature with predictive capacity to identify IBD from fecal samples. Different Machine Learning models were trained, obtaining high performance measures. The predictive capacity of the identified signature was validated in two external cohorts. More precisely a cohort containing samples from patients suffering Ulcerative Colitis and another from patients suffering Crohn's Disease, the two major subtypes of IBD. The results obtained in this validation (AUC 0.74 and AUC = 0.76, respectively) show that our signature presents a generalization capacity in both subtypes. The study of the variables within the model, and a correlation study based on text mining, identified different genera that play an important and common role in the development of these two subtypes.

The Mettl3 epitranscriptomic writer amplifies p53 stress responses.

The p53 transcription factor drives anti-proliferative gene expression programs in response to diverse stressors, including DNA damage and oncogenic signaling. Here, we seek to uncover new mechanisms through which p53 regulates gene expression using tandem affinity purification/mass spectrometry to identify p53-interacting proteins. This approach identified METTL3, an m6A RNA-methyltransferase complex (MTC) constituent, as a p53 interactor. We find that METTL3 promotes p53 protein stabilization and target gene expression in response to DNA damage and oncogenic signals, by both catalytic activity-dependent and independent mechanisms. METTL3 also enhances p53 tumor suppressor activity in in vivo mouse cancer models and human cancer cells. Notably, METTL3 only promotes tumor suppression in the context of intact p53. Analysis of human cancer genome data further supports the notion that the MTC reinforces p53 function in human cancer. Together, these studies reveal a fundamental role for METTL3 in amplifying p53 signaling in response to cellular stress.

Inter-cellular CRISPR screens reveal regulators of cancer cell phagocytosis.

Monoclonal antibody therapies targeting tumour antigens drive cancer cell elimination in large part by triggering macrophage phagocytosis of cancer cells1-7. However, cancer cells evade phagocytosis using mechanisms that are incompletely understood. Here we develop a platform for unbiased identification of factors that impede antibody-dependent cellular phagocytosis (ADCP) using complementary genome-wide CRISPR knockout and overexpression screens in both cancer cells and macrophages. In cancer cells, beyond known factors such as CD47, we identify many regulators of susceptibility to ADCP, including the poorly characterized enzyme adipocyte plasma membrane-associated protein (APMAP). We find that loss of APMAP synergizes with tumour antigen-targeting monoclonal antibodies and/or CD47-blocking monoclonal antibodies to drive markedly increased phagocytosis across a wide range of cancer cell types, including those that are otherwise resistant to ADCP. Additionally, we show that APMAP loss synergizes with several different tumour-targeting monoclonal antibodies to inhibit tumour growth in mice. Using genome-wide counterscreens in macrophages, we find that the G-protein-coupled receptor GPR84 mediates enhanced phagocytosis of APMAP-deficient cancer cells. This work reveals a cancer-intrinsic regulator of susceptibility to antibody-driven phagocytosis and, more broadly, expands our knowledge of the mechanisms governing cancer resistance to macrophage phagocytosis.

Machine Learning Algorithms Reveals Country-Specific Metagenomic Taxa from American Gut Project Data.

In recent years, microbiota has become an increasingly relevant factor for the understanding and potential treatment of diseases. In this work, based on the data reported by the largest study of microbioma in the world, a classification model has been developed based on Machine Learning (ML) capable of predicting the country of origin (United Kingdom vs United States) according to metagenomic data. The data were used for the training of a glmnet algorithm and a Random Forest algorithm. Both algorithms obtained similar results (0.698 and 0.672 in AUC, respectively). Furthermore, thanks to the application of a multivariate feature selection algorithm, eleven metagenomic genres highly correlated with the country of origin were obtained. An in-depth study of the variables used in each model is shown in the present work.

Cell of Origin Influences Pancreatic Cancer Subtype.

Pancreatic ductal adenocarcinoma (PDAC) is a deadly disease with a 5-year survival rate of approximately 9%. An improved understanding of PDAC initiation and progression is paramount for discovering strategies to better detect and combat this disease. Although transcriptomic analyses have uncovered distinct molecular subtypes of human PDAC, the factors that influence subtype development remain unclear. Here, we interrogate the impact of cell of origin and different Trp53 alleles on tumor evolution, using a panel of tractable genetically engineered mouse models. Oncogenic KRAS expression, coupled with Trp53 deletion or point mutation, drives PDAC from both acinar and ductal cells. Gene-expression analysis reveals further that ductal cell-derived and acinar cell-derived tumor signatures are enriched in basal-like and classical subtypes of human PDAC, respectively. These findings highlight cell of origin as one factor that influences PDAC molecular subtypes and provide insight into the fundamental impact that the very earliest events in carcinogenesis can have on cancer evolution. SIGNIFICANCE: Although human PDAC has been classified into different molecular subtypes, the etiology of these distinct subtypes remains unclear. Using mouse genetics, we reveal that cell of origin is an important determinant of PDAC molecular subtype. Deciphering the biology underlying pancreatic cancer subtypes may reveal meaningful distinctions that could improve clinical intervention.This article is highlighted in the In This Issue feature, p. 521.

The AMBRA1 E3 ligase adaptor regulates the stability of cyclin D.

The initiation of cell division integrates a large number of intra- and extracellular inputs. D-type cyclins (hereafter, cyclin D) couple these inputs to the initiation of DNA replication1. Increased levels of cyclin D promote cell division by activating cyclin-dependent kinases 4 and 6 (hereafter, CDK4/6), which in turn phosphorylate and inactivate the retinoblastoma tumour suppressor. Accordingly, increased levels and activity of cyclin D-CDK4/6 complexes are strongly linked to unchecked cell proliferation and cancer2,3. However, the mechanisms that regulate levels of cyclin D are incompletely understood4,5. Here we show that autophagy and beclin 1 regulator 1 (AMBRA1) is the main regulator of the degradation of cyclin D. We identified AMBRA1 in a genome-wide screen to investigate the genetic basis of  the response to CDK4/6 inhibition. Loss of AMBRA1 results in high levels of cyclin D in cells and in mice, which promotes proliferation and decreases sensitivity to CDK4/6 inhibition. Mechanistically, AMBRA1 mediates ubiquitylation and proteasomal degradation of cyclin D as a substrate receptor for the cullin 4 E3 ligase complex. Loss of AMBRA1 enhances the growth of lung adenocarcinoma in a mouse model, and low levels of AMBRA1 correlate with worse survival in patients with lung adenocarcinoma. Thus, AMBRA1 regulates cellular levels of cyclin D, and contributes to cancer development and the response of cancer cells to CDK4/6 inhibitors.

A CRISPR/Cas9-Engineered ARID1A-Deficient Human Gastric Cancer Organoid Model Reveals Essential and Nonessential Modes of Oncogenic Transformation.

Mutations in ARID1A rank among the most common molecular aberrations in human cancer. However, oncogenic consequences of ARID1A mutation in human cells remain poorly defined due to lack of forward genetic models. Here, CRISPR/Cas9-mediated ARID1A knockout (KO) in primary TP53-/- human gastric organoids induced morphologic dysplasia, tumorigenicity, and mucinous differentiation. Genetic WNT/ß-catenin activation rescued mucinous differentiation, but not hyperproliferation, suggesting alternative pathways of ARID1A KO-mediated transformation. ARID1A mutation induced transcriptional regulatory modules characteristic of microsatellite instability and Epstein-Barr virus-associated subtype human gastric cancer, including FOXM1-associated mitotic genes and BIRC5/survivin. Convergently, high-throughput compound screening indicated selective vulnerability of ARID1A-deficient organoids to inhibition of BIRC5/survivin, functionally implicating this pathway as an essential mediator of ARID1A KO-dependent early-stage gastric tumorigenesis. Overall, we define distinct pathways downstream of oncogenic ARID1A mutation, with nonessential WNT-inhibited mucinous differentiation in parallel with essential transcriptional FOXM1/BIRC5-stimulated proliferation, illustrating the general utility of organoid-based forward genetic cancer analysis in human cells. SIGNIFICANCE: We establish the first human forward genetic modeling of a commonly mutated tumor suppressor gene, ARID1A. Our study integrates diverse modalities including CRISPR/Cas9 genome editing, organoid culture, systems biology, and small-molecule screening to derive novel insights into early transformation mechanisms of ARID1A-deficient gastric cancers.See related commentary by Zafra and Dow, p. 1327.This article is highlighted in the In This Issue feature, p. 1307.

Zmat3 Is a Key Splicing Regulator in the p53 Tumor Suppression Program.

Although TP53 is the most commonly mutated gene in human cancers, the p53-dependent transcriptional programs mediating tumor suppression remain incompletely understood. Here, to uncover critical components downstream of p53 in tumor suppression, we perform unbiased RNAi and CRISPR-Cas9-based genetic screens in vivo. These screens converge upon the p53-inducible gene Zmat3, encoding an RNA-binding protein, and we demonstrate that ZMAT3 is an important tumor suppressor downstream of p53 in mouse KrasG12D-driven lung and liver cancers and human carcinomas. Integrative analysis of the ZMAT3 RNA-binding landscape and transcriptomic profiling reveals that ZMAT3 directly modulates exon inclusion in transcripts encoding proteins of diverse functions, including the p53 inhibitors MDM4 and MDM2, splicing regulators, and components of varied cellular processes. Interestingly, these exons are enriched in NMD signals, and, accordingly, ZMAT3 broadly affects target transcript stability. Collectively, these studies reveal ZMAT3 as a novel RNA-splicing and homeostasis regulator and a key component of p53-mediated tumor suppression.

The m6A RNA demethylase FTO is a HIF-independent synthetic lethal partner with the VHL tumor suppressor.

Loss of the von Hippel-Lindau (VHL) tumor suppressor is a hallmark feature of renal clear cell carcinoma. VHL inactivation results in the constitutive activation of the hypoxia-inducible factors (HIFs) HIF-1 and HIF-2 and their downstream targets, including the proangiogenic factors VEGF and PDGF. However, antiangiogenic agents and HIF-2 inhibitors have limited efficacy in cancer therapy due to the development of resistance. Here we employed an innovative computational platform, Mining of Synthetic Lethals (MiSL), to identify synthetic lethal interactions with the loss of VHL through analysis of primary tumor genomic and transcriptomic data. Using this approach, we identified a synthetic lethal interaction between VHL and the m6A RNA demethylase FTO in renal cell carcinoma. MiSL identified FTO as a synthetic lethal partner of VHL because deletions of FTO are mutually exclusive with VHL loss in pan cancer datasets. Moreover, FTO expression is increased in VHL-deficient ccRCC tumors compared to normal adjacent tissue. Genetic inactivation of FTO using multiple orthogonal approaches revealed that FTO inhibition selectively reduces the growth and survival of VHL-deficient cells in vitro and in vivo. Notably, FTO inhibition reduced the survival of both HIF wild type and HIF-deficient tumors, identifying FTO as an HIF-independent vulnerability of VHL-deficient cancers. Integrated analysis of transcriptome-wide m6A-seq and mRNA-seq analysis identified the glutamine transporter SLC1A5 as an FTO target that promotes metabolic reprogramming and survival of VHL-deficient ccRCC cells. These findings identify FTO as a potential HIF-independent therapeutic target for the treatment of VHL-deficient renal cell carcinoma.

CRISPR screens in cancer spheroids identify 3D growth-specific vulnerabilities.

Cancer genomics studies have identified thousands of putative cancer driver genes1. Development of high-throughput and accurate models to define the functions of these genes is a major challenge. Here we devised a scalable cancer-spheroid model and performed genome-wide CRISPR screens in 2D monolayers and 3D lung-cancer spheroids. CRISPR phenotypes in 3D more accurately recapitulated those of in vivo tumours, and genes with differential sensitivities between 2D and 3D conditions were highly enriched for genes that are mutated in lung cancers. These analyses also revealed drivers that are essential for cancer growth in 3D and in vivo, but not in 2D. Notably, we found that carboxypeptidase D is responsible for removal of a C-terminal RKRR motif2 from the α-chain of the insulin-like growth factor 1 receptor that is critical for receptor activity. Carboxypeptidase D expression correlates with patient outcomes in patients with lung cancer, and loss of carboxypeptidase D reduced tumour growth. Our results reveal key differences between 2D and 3D cancer models, and establish a generalizable strategy for performing CRISPR screens in spheroids to reveal cancer vulnerabilities.

Metabolic Profiling Reveals a Dependency of Human Metastatic Breast Cancer on Mitochondrial Serine and One-Carbon Unit Metabolism.

Breast cancer is the most common cancer among American women and a major cause of mortality. To identify metabolic pathways as potential targets to treat metastatic breast cancer, we performed metabolomics profiling on the breast cancer cell line MDA-MB-231 and its tissue-tropic metastatic subclones. Here, we report that these subclones with increased metastatic potential display an altered metabolic profile compared with the parental population. In particular, the mitochondrial serine and one-carbon (1C) unit pathway is upregulated in metastatic subclones. Mechanistically, the mitochondrial serine and 1C unit pathway drives the faster proliferation of subclones through enhanced de novo purine biosynthesis. Inhibition of the first rate-limiting enzyme of the mitochondrial serine and 1C unit pathway, serine hydroxymethyltransferase (SHMT2), potently suppresses proliferation of metastatic subclones in culture and impairs growth of lung metastatic subclones at both primary and metastatic sites in mice. Some human breast cancers exhibit a significant association between the expression of genes in the mitochondrial serine and 1C unit pathway with disease outcome and higher expression of SHMT2 in metastatic tumor tissue compared with primary tumors. In addition to breast cancer, a few other cancer types, such as adrenocortical carcinoma and kidney chromophobe cell carcinoma, also display increased SHMT2 expression during disease progression. Together, these results suggest that mitochondrial serine and 1C unit metabolism plays an important role in promoting cancer progression, particularly in late-stage cancer. IMPLICATIONS: This study identifies mitochondrial serine and 1C unit metabolism as an important pathway during the progression of a subset of human breast cancers.

Chromatin regulators mediate anthracycline sensitivity in breast cancer.

Anthracyclines are a highly effective component of curative breast cancer chemotherapy but are associated with substantial morbidity1,2. Because anthracyclines work in part by inhibiting topoisomerase-II (TOP2) on accessible DNA3,4, we hypothesized that chromatin regulatory genes (CRGs) that mediate DNA accessibility might predict anthracycline response. We studied the role of CRGs in anthracycline sensitivity in breast cancer through integrative analysis of patient and cell line data. We identified a consensus set of 38 CRGs associated with anthracycline response across ten cell line datasets. By evaluating the interaction between expression and treatment in predicting survival in a metacohort of 1006 patients with early-stage breast cancer, we identified 54 CRGs whose expression levels dictate anthracycline benefit across the clinical subgroups; of these CRGs, 12 overlapped with those identified in vitro. CRGs that promote DNA accessibility, including Trithorax complex members, were associated with anthracycline sensitivity when highly expressed, whereas CRGs that reduce accessibility, such as Polycomb complex proteins, were associated with decreased anthracycline sensitivity. We show that KDM4B modulates TOP2 accessibility to chromatin, elucidating a mechanism of TOP2 inhibitor sensitivity. These findings indicate that CRGs mediate anthracycline benefit by altering DNA accessibility, with implications for the stratification of patients with breast cancer and treatment decision making.

Quantitative evidence for early metastatic seeding in colorectal cancer.

Both the timing and molecular determinants of metastasis are unknown, hindering treatment and prevention efforts. Here we characterize the evolutionary dynamics of this lethal process by analyzing exome-sequencing data from 118 biopsies from 23 patients with colorectal cancer with metastases to the liver or brain. The data show that the genomic divergence between the primary tumor and metastasis is low and that canonical driver genes were acquired early. Analysis within a spatial tumor growth model and statistical inference framework indicates that early disseminated cells commonly (81%, 17 out of 21 evaluable patients) seed metastases while the carcinoma is clinically undetectable (typically, less than 0.01 cm3). We validated the association between early drivers and metastasis in an independent cohort of 2,751 colorectal cancers, demonstrating their utility as biomarkers of metastasis. This conceptual and analytical framework provides quantitative in vivo evidence that systemic spread can occur early in colorectal cancer and illuminates strategies for patient stratification and therapeutic targeting of the canonical drivers of tumorigenesis.