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.

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.

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.

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.

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.

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.

Dynamics of breast-cancer relapse reveal late-recurring ER-positive genomic subgroups.

The rates and routes of lethal systemic spread in breast cancer are poorly understood owing to a lack of molecularly characterized patient cohorts with long-term, detailed follow-up data. Long-term follow-up is especially important for those with oestrogen-receptor (ER)-positive breast cancers, which can recur up to two decades after initial diagnosis1-6. It is therefore essential to identify patients who have a high risk of late relapse7-9. Here we present a statistical framework that models distinct disease stages (locoregional recurrence, distant recurrence, breast-cancer-related death and death from other causes) and competing risks of mortality from breast cancer, while yielding individual risk-of-recurrence predictions. We apply this model to 3,240 patients with breast cancer, including 1,980 for whom molecular data are available, and delineate spatiotemporal patterns of relapse across different categories of molecular information (namely immunohistochemical subtypes; PAM50 subtypes, which are based on gene-expression patterns10,11; and integrative or IntClust subtypes, which are based on patterns of genomic copy-number alterations and gene expression12,13). We identify four late-recurring integrative subtypes, comprising about one quarter (26%) of tumours that are both positive for ER and negative for human epidermal growth factor receptor 2, each with characteristic tumour-driving alterations in genomic copy number and a high risk of recurrence (mean 47-62%) up to 20 years after diagnosis. We also define a subgroup of triple-negative breast cancers in which cancer rarely recurs after five years, and a separate subgroup in which patients remain at risk. Use of the integrative subtypes improves the prediction of late, distant relapse beyond what is possible with clinical covariates (nodal status, tumour size, tumour grade and immunohistochemical subtype). These findings highlight opportunities for improved patient stratification and biomarker-driven clinical trials.

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.

Genome co-amplification upregulates a mitotic gene network activity that predicts outcome and response to mitotic protein inhibitors in breast cancer.

BACKGROUND: High mitotic activity is associated with the genesis and progression of many cancers. Small molecule inhibitors of mitotic apparatus proteins are now being developed and evaluated clinically as anticancer agents. With clinical trials of several of these experimental compounds underway, it is important to understand the molecular mechanisms that determine high mitotic activity, identify tumor subtypes that carry molecular aberrations that confer high mitotic activity, and to develop molecular markers that distinguish which tumors will be most responsive to mitotic apparatus inhibitors. METHODS: We identified a coordinately regulated mitotic apparatus network by analyzing gene expression profiles for 53 malignant and non-malignant human breast cancer cell lines and two separate primary breast tumor datasets. We defined the mitotic network activity index (MNAI) as the sum of the transcriptional levels of the 54 coordinately regulated mitotic apparatus genes. The effect of those genes on cell growth was evaluated by small interfering RNA (siRNA). RESULTS: High MNAI was enriched in basal-like breast tumors and was associated with reduced survival duration and preferential sensitivity to inhibitors of the mitotic apparatus proteins, polo-like kinase, centromere associated protein E and aurora kinase designated GSK462364, GSK923295 and GSK1070916, respectively. Co-amplification of regions of chromosomes 8q24, 10p15-p12, 12p13, and 17q24-q25 was associated with the transcriptional upregulation of this network of 54 mitotic apparatus genes, and we identify transcription factors that localize to these regions and putatively regulate mitotic activity. Knockdown of the mitotic network by siRNA identified 22 genes that might be considered as additional therapeutic targets for this clinically relevant patient subgroup. CONCLUSIONS: We define a molecular signature which may guide therapeutic approaches for tumors with high mitotic network activity.

A pathway-based data integration framework for prediction of disease progression.

MOTIVATION: Within medical research there is an increasing trend toward deriving multiple types of data from the same individual. The most effective prognostic prediction methods should use all available data, as this maximizes the amount of information used. In this article, we consider a variety of learning strategies to boost prediction performance based on the use of all available data. IMPLEMENTATION: We consider data integration via the use of multiple kernel learning supervised learning methods. We propose a scheme in which feature selection by statistical score is performed separately per data type and by pathway membership. We further consider the introduction of a confidence measure for the class assignment, both to remove some ambiguously labeled datapoints from the training data and to implement a cautious classifier that only makes predictions when the associated confidence is high. RESULTS: We use the METABRIC dataset for breast cancer, with prediction of survival at 2000 days from diagnosis. Predictive accuracy is improved by using kernels that exclusively use those genes, as features, which are known members of particular pathways. We show that yet further improvements can be made by using a range of additional kernels based on clinical covariates such as Estrogen Receptor (ER) status. Using this range of measures to improve prediction performance, we show that the test accuracy on new instances is nearly 80%, though predictions are only made on 69.2% of the patient cohort. AVAILABILITY: https://github.com/jseoane/FSMKL CONTACT: J.Seoane@bristol.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Classification of signaling proteins based on molecular star graph descriptors using Machine Learning models.

Signaling proteins are an important topic in drug development due to the increased importance of finding fast, accurate and cheap methods to evaluate new molecular targets involved in specific diseases. The complexity of the protein structure hinders the direct association of the signaling activity with the molecular structure. Therefore, the proposed solution involves the use of protein star graphs for the peptide sequence information encoding into specific topological indices calculated with S2SNet tool. The Quantitative Structure-Activity Relationship classification model obtained with Machine Learning techniques is able to predict new signaling peptides. The best classification model is the first signaling prediction model, which is based on eleven descriptors and it was obtained using the Support Vector Machines-Recursive Feature Elimination (SVM-RFE) technique with the Laplacian kernel (RFE-LAP) and an AUROC of 0.961. Testing a set of 3114 proteins of unknown function from the PDB database assessed the prediction performance of the model. Important signaling pathways are presented for three UniprotIDs (34 PDBs) with a signaling prediction greater than 98.0%.

Canonical correlation analysis for gene-based pleiotropy discovery.

Genome-wide association studies have identified a wealth of genetic variants involved in complex traits and multifactorial diseases. There is now considerable interest in testing variants for association with multiple phenotypes (pleiotropy) and for testing multiple variants for association with a single phenotype (gene-based association tests). Such approaches can increase statistical power by combining evidence for association over multiple phenotypes or genetic variants respectively. Canonical Correlation Analysis (CCA) measures the correlation between two sets of multidimensional variables, and thus offers the potential to combine these two approaches. To apply CCA, we must restrict the number of attributes relative to the number of samples. Hence we consider modules of genetic variation that can comprise a gene, a pathway or another biologically relevant grouping, and/or a set of phenotypes. In order to do this, we use an attribute selection strategy based on a binary genetic algorithm. Applied to a UK-based prospective cohort study of 4286 women (the British Women's Heart and Health Study), we find improved statistical power in the detection of previously reported genetic associations, and identify a number of novel pleiotropic associations between genetic variants and phenotypes. New discoveries include gene-based association of NSF with triglyceride levels and several genes (ACSM3, ERI2, IL18RAP, IL23RAP and NRG1) with left ventricular hypertrophy phenotypes. In multiple-phenotype analyses we find association of NRG1 with left ventricular hypertrophy phenotypes, fibrinogen and urea and pleiotropic relationships of F7 and F10 with Factor VII, Factor IX and cholesterol levels.

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.

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.

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.

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.

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.