Fundamental immune-oncogenicity trade-offs define driver mutation fitness.

Missense driver mutations in cancer are concentrated in a few hotspots1. Various mechanisms have been proposed to explain this skew, including biased mutational processes2, phenotypic differences3-6 and immunoediting of neoantigens7,8; however, to our knowledge, no existing model weighs the relative contribution of these features to tumour evolution. We propose a unified theoretical 'free fitness' framework that parsimoniously integrates multimodal genomic, epigenetic, transcriptomic and proteomic data into a biophysical model of the rate-limiting processes underlying the fitness advantage conferred on cancer cells by driver gene mutations. Focusing on TP53, the most mutated gene in cancer1, we present an inference of mutant p53 concentration and demonstrate that TP53 hotspot mutations optimally solve an evolutionary trade-off between oncogenic potential and neoantigen immunogenicity. Our model anticipates patient survival in The Cancer Genome Atlas and patients with lung cancer treated with immunotherapy as well as the age of tumour onset in germline carriers of TP53 variants. The predicted differential immunogenicity between hotspot mutations was validated experimentally in patients with cancer and in a unique large dataset of healthy individuals. Our data indicate that immune selective pressure on TP53 mutations has a smaller role in non-cancerous lesions than in tumours, suggesting that targeted immunotherapy may offer an early prophylactic opportunity for the former. Determining the relative contribution of immunogenicity and oncogenic function to the selective advantage of hotspot mutations thus has important implications for both precision immunotherapies and our understanding of tumour evolution.

Mutant p53 disrupts mammary tissue architecture via the mevalonate pathway.

p53 is a frequent target for mutation in human tumors, and mutant p53 proteins can actively contribute to tumorigenesis. We employed a three-dimensional culture model in which nonmalignant breast epithelial cells form spheroids reminiscent of acinar structures found in vivo, whereas breast cancer cells display highly disorganized morphology. We found that mutant p53 depletion is sufficient to phenotypically revert breast cancer cells to a more acinar-like morphology. Genome-wide expression analysis identified the mevalonate pathway as significantly upregulated by mutant p53. Statins and sterol biosynthesis intermediates reveal that this pathway is both necessary and sufficient for the phenotypic effects of mutant p53 on breast tissue architecture. Mutant p53 associates with sterol gene promoters at least partly via SREBP transcription factors. Finally, p53 mutation correlates with highly expressed sterol biosynthesis genes in human breast tumors. These findings implicate the mevalonate pathway as a therapeutic target for tumors bearing mutations in p53.

The interplay between epigenetic changes and the p53 protein in stem cells.

Epigenetic programs regulate the development and maintenance of organisms over a lifetime. These programs are carried out through chemical modifications of DNA and proteins such as histones and transcription factors. These epigenetic modifications are less stable than genetic alterations and even reversible under a variety of circumstances, such as developmental changes, regeneration of tissues, cell divisions, aging, and pathological conditions observed in many cancers. The p53 protein not only enforces the stability of the genome by the prevention of genetic alterations in cells but also plays a role in regulating the epigenetic changes that can occur in cells. The full-length p53 protein is largely inactive in stem cells but, when activated, helps to commit these cells to developmental lineages through a series of epigenetic changes. Just as p53 impacts epigenetic change, the enzyme activities that carry out epigenetic protein modifications act on the p53 protein and its splice variants in stem and progenitor cells to silence or activate its transcriptional activities. Thus, there is a great deal of cross-talk between the p53 protein and epigenetic programs. This review collects the diverse experimental evidence that leads to these conclusions. This in turn permits new ideas and directions for the treatment of cancers, reactivating developmental pathways for tissue regeneration and responses to the impact of aging.

The p53 family: guardians of maternal reproduction.

The p53 family of proteins consists of p53, p63 and p73, which are transcription factors that affect both cancer and development. It is now emerging that these proteins also regulate maternal reproduction. Whereas p63 is important for maturation of the egg, p73 ensures normal mitosis in the developing blastocyst. p53 subsequently regulates implantation of the embryo through transcriptional control of leukaemia inhibitory factor. Elucidating the cell biological basis of how these factors regulate female fertility may lead to new approaches to the control of human maternal reproduction.

A neoantigen fitness model predicts tumour response to checkpoint blockade immunotherapy.

Checkpoint blockade immunotherapies enable the host immune system to recognize and destroy tumour cells. Their clinical activity has been correlated with activated T-cell recognition of neoantigens, which are tumour-specific, mutated peptides presented on the surface of cancer cells. Here we present a fitness model for tumours based on immune interactions of neoantigens that predicts response to immunotherapy. Two main factors determine neoantigen fitness: the likelihood of neoantigen presentation by the major histocompatibility complex (MHC) and subsequent recognition by T cells. We estimate these components using the relative MHC binding affinity of each neoantigen to its wild type and a nonlinear dependence on sequence similarity of neoantigens to known antigens. To describe the evolution of a heterogeneous tumour, we evaluate its fitness as a weighted effect of dominant neoantigens in the subclones of the tumour. Our model predicts survival in anti-CTLA-4-treated patients with melanoma and anti-PD-1-treated patients with lung cancer. Importantly, low-fitness neoantigens identified by our method may be leveraged for developing novel immunotherapies. By using an immune fitness model to study immunotherapy, we reveal broad similarities between the evolution of tumours and rapidly evolving pathogens.

Identification of unique neoantigen qualities in long-term survivors of pancreatic cancer.

Pancreatic ductal adenocarcinoma is a lethal cancer with fewer than 7% of patients surviving past 5 years. T-cell immunity has been linked to the exceptional outcome of the few long-term survivors, yet the relevant antigens remain unknown. Here we use genetic, immunohistochemical and transcriptional immunoprofiling, computational biophysics, and functional assays to identify T-cell antigens in long-term survivors of pancreatic cancer. Using whole-exome sequencing and in silico neoantigen prediction, we found that tumours with both the highest neoantigen number and the most abundant CD8+ T-cell infiltrates, but neither alone, stratified patients with the longest survival. Investigating the specific neoantigen qualities promoting T-cell activation in long-term survivors, we discovered that these individuals were enriched in neoantigen qualities defined by a fitness model, and neoantigens in the tumour antigen MUC16 (also known as CA125). A neoantigen quality fitness model conferring greater immunogenicity to neoantigens with differential presentation and homology to infectious disease-derived peptides identified long-term survivors in two independent datasets, whereas a neoantigen quantity model ascribing greater immunogenicity to increasing neoantigen number alone did not. We detected intratumoural and lasting circulating T-cell reactivity to both high-quality and MUC16 neoantigens in long-term survivors of pancreatic cancer, including clones with specificity to both high-quality neoantigens and predicted cross-reactive microbial epitopes, consistent with neoantigen molecular mimicry. Notably, we observed selective loss of high-quality and MUC16 neoantigenic clones on metastatic progression, suggesting neoantigen immunoediting. Our results identify neoantigens with unique qualities as T-cell targets in pancreatic ductal adenocarcinoma. More broadly, we identify neoantigen quality as a biomarker for immunogenic tumours that may guide the application of immunotherapies.

MicroRNA-101 Suppresses Tumor Cell Proliferation by Acting as an Endogenous Proteasome Inhibitor via Targeting the Proteasome Assembly Factor POMP.

Proteasome inhibition represents a promising strategy of cancer pharmacotherapy, but resistant tumor cells often emerge. Here we show that the microRNA-101 (miR-101) targets the proteasome maturation protein POMP, leading to impaired proteasome assembly and activity, and resulting in accumulation of p53 and cyclin-dependent kinase inhibitors, cell cycle arrest, and apoptosis. miR-101-resistant POMP restores proper turnover of proteasome substrates and re-enables tumor cell growth. In ERα-positive breast cancers, miR-101 and POMP levels are inversely correlated, and high miR-101 expression or low POMP expression associates with prolonged survival. Mechanistically, miR-101 expression or POMP knockdown attenuated estrogen-driven transcription. Finally, suppressing POMP is sufficient to overcome tumor cell resistance to the proteasome inhibitor bortezomib. Taken together, proteasome activity can not only be manipulated through drugs, but is also subject to endogenous regulation through miR-101, which targets proteasome biogenesis to control overall protein turnover and tumor cell proliferation.

Functional network analysis reveals an immune tolerance mechanism in cancer.

We present a technique to construct a simplification of a feature network which can be used for interactive data exploration, biological hypothesis generation, and the detection of communities or modules of cofunctional features. These are modules of features that are not necessarily correlated, but nevertheless exhibit common function in their network context as measured by similarity of relationships with neighboring features. In the case of genetic networks, traditional pathway analyses tend to assume that, ideally, all genes in a module exhibit very similar function, independent of relationships with other genes. The proposed technique explicitly relaxes this assumption by employing the comparison of relational profiles. For example, two genes which always activate a third gene are grouped together even if they never do so concurrently. They have common, but not identical, function. The comparison is driven by an average of a certain computationally efficient comparison metric between Gaussian mixture models. The method has its basis in the local connection structure of the network and the collection of joint distributions of the data associated with nodal neighborhoods. It is benchmarked on networks with known community structures. As the main application, we analyzed the gene regulatory network in lung adenocarcinoma, finding a cofunctional module of genes including the pregnancy-specific glycoproteins (PSGs). About 20% of patients with lung, breast, uterus, and colon cancer in The Cancer Genome Atlas (TCGA) have an elevated PSG+ signature, with associated poor group prognosis. In conjunction with previous results relating PSGs to tolerance in the immune system, these findings implicate the PSGs in a potential immune tolerance mechanism of cancers.

The evolution of thymic lymphomas in p53 knockout mice.

Germline deletion of the p53 gene in mice gives rise to spontaneous thymic (T-cell) lymphomas. In this study, the p53 knockout mouse was employed as a model to study the mutational evolution of tumorigenesis. The clonality of the T-cell repertoire from p53 knockout and wild-type thymic cells was analyzed at various ages employing TCRß sequencing. These data demonstrate that p53 knockout thymic lymphomas arose in an oligoclonal fashion, with tumors evolving dominant clones over time. Exon sequencing of tumor DNA revealed that all of the independently derived oligoclonal mouse tumors had a deletion in the Pten gene prior to the formation of the TCRß rearrangement, produced early in development. This was followed in each independent clone of the thymic lymphoma by the amplification or overexpression of cyclin Ds and Cdk6. Alterations in the expression of Ikaros were common and blocked further development of CD-4/CD-8 T cells. While the frequency of point mutations in the genome of these lymphomas was one per megabase, there were a tremendous number of copy number variations producing the tumors' driver mutations. The initial inherited loss of p53 functions appeared to delineate an order of genetic alterations selected for during the evolution of these thymic lymphomas.

Mdm2: Open questions.

The Mdm2 oncoprotein and its association with p53 were discovered 30 years ago, and a cornucopia of activities and regulatory pathways have been associated with it. In this review, we will raise questions about Mdm2 and its cousin Mdm4 that we consider worth pursuing in future research, reaching from molecular structures and intracellular activities all the way to development, evolution, and cancer therapy. We anticipate that such research will not only close a few gaps in our knowledge but could add new dimensions to our current view. This compilation of questions contributes to the preparation for the 10th Mdm2 Workshop in Tokyo.

Insights into immune system development and function from mouse T-cell repertoires.

The ability of the adaptive immune system to respond to arbitrary pathogens stems from the broad diversity of immune cell surface receptors. This diversity originates in a stochastic DNA editing process (VDJ recombination) that acts on the surface receptor gene each time a new immune cell is created from a stem cell. By analyzing T-cell receptor (TCR) sequence repertoires taken from the blood and thymus of mice of different ages, we quantify the changes in the VDJ recombination process that occur from embryo to young adult. We find a rapid increase with age in the number of random insertions and a dramatic increase in diversity. Because the blood accumulates thymic output over time, blood repertoires are mixtures of different statistical recombination processes, and we unravel the mixture statistics to obtain a picture of the time evolution of the early immune system. Sequence repertoire analysis also allows us to detect the statistical impact of selection on the output of the VDJ recombination process. The effects we find are nearly identical between thymus and blood, suggesting that our analysis mainly detects selection for proper folding of the TCR receptor protein. We further find that selection is weaker in laboratory mice than in humans and it does not affect the diversity of the repertoire.

P53 and The Immune Response: 40 Years of Exploration-A Plan for the Future.

The p53 field was born from a marriage of the techniques of cancer virus research and immunology. Over the past 40 years, it has followed the path of cancer research. Now cancer treatments are turning to immunotherapy, and there are many hints of the role of the p53 protein in both the regulation of the innate immune system and as an antigen in adaptive immune responses. The p53 gene and protein are part of the innate immune system, and play an important role in infectious diseases, senescence, aging, and the surveillance of repetitive DNA and RNAs. The mutant form of the p53 protein in cancers elicits both a B-cell antibody response (a tumor antigen) and a CD-8 killer T-cell response (a tumor-specific transplantation antigen). The future will take the p53-immune response field of research into cancer immunotherapy, autoimmunity, inflammatory responses, neuro-degeneration, aging, and life span, and the regulation of epigenetic stability and tissue regeneration. The next 40 years will bring the p53 gene and its proteins out of a cancer focus and into an organismic and environmental focus.

Lysine methylation represses p53 activity in teratocarcinoma cancer cells.

TP53 (which encodes the p53 protein) is the most frequently mutated gene among all human cancers, whereas tumors that retain the wild-type TP53 gene often use alternative mechanisms to repress the p53 tumor-suppressive function. Testicular teratocarcinoma cells rarely contain mutations in TP53, yet the transcriptional activity of wild-type p53 is compromised, despite its high expression level. Here we report that in the teratocarcinoma cell line NTera2, p53 is subject to lysine methylation at its carboxyl terminus, which has been shown to repress p53's transcriptional activity. We show that reduction of the cognate methyltransferases reactivates p53 and promotes differentiation of the NTera2 cells. Furthermore, reconstitution of methylation-deficient p53 mutants into p53-depleted NTera2 cells results in elevated expression of p53 downstream targets and precocious loss of pluripotent gene expression compared with re-expression of wild-type p53. Our results provide evidence that lysine methylation of endogenous wild-type p53 represses its activity in cancer cells and suggest new therapeutic possibilities of targeting testicular teratocarcinoma.

The Evolution of Tumor Formation in Humans and Mice with Inherited Mutations in the p53 Gene.

While tumors are very heterogeneous in their origins, mutations in the p53 gene and inactivation of p53 gene functions are the most common feature that predispose to the formation of cancers in humans. Inherited p53 mutations lead to different tumor types at very different frequencies and at very different ages than somatic p53 mutations. The reasons for this are explored. When the first mutation arises in a stem cell (a gatekeeper mutation) it selects for a specific subset of second mutations which in turn select for mutations in a third subset of genes. The nature of the first mutation in a tumor determines, by selection, the functional types of subsequent mutations. Inherited mutations occur at different developmental times and in different orders of mutational sequences than somatic mutations. The excess risk of developing a cancer with an inherited p53 mutation is two- to three-fold in endodermal derived tissues compared with 100- to 1000-fold for ectodermal and mesenchymal derived tissues. By contrast, endodermal derived tumors with somatic p53 mutations occur at very high frequencies (70-100%). These evolutionary restrictions upon the mutational path that tumor development may take could open up new avenues for therapy and prevention.

Negative regulation of tumor suppressor p53 by microRNA miR-504.

Tumor suppressor p53 plays a central role in tumor prevention. p53 protein levels and activity are under a tight and complex regulation in cells to maintain the proper function of p53. MicroRNAs play a key role in the regulation of gene expression. Here we report the regulation of p53 through miR-504. miR-504 acts as a negative regulator of human p53 through its direct binding to two sites in the p53 3' untranslated region. Overexpression of miR-504 decreases p53 protein levels and functions in cells, including p53 transcriptional activity, p53-mediated apoptosis, and cell-cycle arrest in response to stress, and furthermore promotes tumorigenecity of cells in vivo. These results demonstrate the direct negative regulation of p53 by miR-504 as a mechanism for p53 regulation in cells, which highlights the importance of microRNAs in tumorigenesis.

P53 and the defenses against genome instability caused by transposons and repetitive elements.

The recent publication by Wylie et al. is reviewed, demonstrating that the p53 protein regulates the movement of transposons. While this work presents genetic evidence for a piRNA-mediated p53 interaction with transposons in Drosophila and zebrafish, it is herein placed in the context of a decade or so of additional work that demonstrated a role for p53 in regulating transposons and other repetitive elements. The line of thought in those studies began with the observation that transposons damage DNA and p53 regulates DNA damage. The presence of transposon movement can increase the rate of evolution in the germ line and alter genes involved in signal transduction pathways. Transposition can also play an important role in cancers where the p53 gene function is often mutated. This is particularly interesting as recent work has shown that de-repression of repetitive elements in cancer has important consequences for the immune system and tumor microenvironment.

DTIE, a novel core promoter element that directs start site selection in TATA-less genes.

The transcription start site (TSS) determines the length and composition of the 5' UTR and therefore can have a profound effect on translation. Yet, little is known about the mechanism underlying start site selection, particularly from promoters lacking conventional core elements such as TATA-box and Initiator. Here we report a novel mechanism of start site selection in the TATA- and Initiator-less promoter of miR-22, through a strictly localized downstream element termed DTIE and an upstream distal element. Changing the distance between them reduced promoter strength, altered TSS selection and diminished Pol II recruitment. Biochemical assays suggest that DTIE does not serve as a docking site for TFIID, the major core promoter-binding factor. TFIID is recruited to the promoter through DTIE but is dispensable for TSS selection. We determined DTIE consensus and found it to be remarkably prevalent, present at the same TSS downstream location in ≈20.8% of human promoters, the vast majority of which are TATA-less. Analysis of DTIE in the tumor suppressor p53 confirmed a similar function. Our findings reveal a novel mechanism of transcription initiation from TATA-less promoters.

TAp73 opposes tumor angiogenesis by promoting hypoxia-inducible factor 1α degradation.

Tumor hypoxia and hypoxia-inducible factor 1 (HIF-1) activation are associated with cancer progression. Here, we demonstrate that the transcription factor TAp73 opposes HIF-1 activity through a nontranscriptional mechanism, thus affecting tumor angiogenesis. TAp73-deficient mice have an increased incidence of spontaneous and chemically induced tumors that also display enhanced vascularization. Mechanistically, TAp73 interacts with the regulatory subunit (α) of HIF-1 and recruits mouse double minute 2 homolog into the protein complex, thus promoting HIF-1α polyubiquitination and consequent proteasomal degradation in an oxygen-independent manner. In human lung cancer datasets, TAp73 strongly predicts good patient prognosis, and its expression is associated with low HIF-1 activation and angiogenesis. Our findings, supported by in vivo and clinical evidence, demonstrate a mechanism for oxygen-independent HIF-1 regulation, which has important implications for individualizing therapies in patients with cancer.