Integration of TP53, DREAM, MMB-FOXM1 and RB-E2F target gene analyses identifies cell cycle gene regulatory networks

At a Crossroads to Cancer: How p53-Induced Cell Fate Decisions Secure Genome Integrity

P53 is known as the most critical tumor suppressor and is often referred to as the guardian of our genome. More than 40 years after its discovery, we are still struggling to understand all molecular details on how this transcription factor prevents oncogenesis or how to leverage current knowledge about its function to improve cancer treatment. Multiple cues, including DNA-damage or mitotic errors, can lead to the stabilization and nuclear translocation of p53, initiating the expression of multiple target genes. These transcriptional programs may be cell-type- and stimulus-specific, as is their outcome that ultimately imposes a barrier to cellular transformation. Cell cycle arrest and cell death are two well-studied consequences of p53 activation, but, while being considered critical, they do not fully explain the consequences of p53 loss-of-function phenotypes in cancer. Here, we discuss how mitotic errors alert the p53 network and give an overview of multiple ways that p53 can trigger cell death. We argue that a comparative analysis of different types of p53 responses, elicited by different triggers in a time-resolved manner in well-defined model systems, is critical to understand the cell-type-specific cell fate induced by p53 upon its activation in order to resolve the remaining mystery of its tumor-suppressive function.

CXCR4 hyperactivation cooperates with TCL1 in CLL development and aggressiveness

Aberrant CXCR4 activity has been implicated in lymphoma pathogenesis, disease progression, and resistance to therapies. Using a mouse model with a gain-of-function CXCR4 mutation (CXCR4C1013G) that hyperactivates CXCR4 signaling, we identified CXCR4 as a crucial activator of multiple key oncogenic pathways. CXCR4 hyperactivation resulted in an expansion of transitional B1 lymphocytes, which represent the precursors of chronic lymphocytic leukemia (CLL). Indeed, CXCR4 hyperactivation led to a significant acceleration of disease onset and a more aggressive phenotype in the murine Eµ-TCL1 CLL model. Hyperactivated CXCR4 signaling cooperated with TCL1 to cause a distinct oncogenic transcriptional program in B cells, characterized by PLK1/FOXM1-associated pathways. In accordance, Eµ-TCL1;CXCR4C1013G B cells enriched a transcriptional signature from patients with Richter's syndrome, an aggressive transformation of CLL. Notably, MYC activation in aggressive lymphoma was associated with increased CXCR4 expression. In line with this finding, additional hyperactive CXCR4 signaling in the Eµ-Myc mouse, a model of aggressive B-cell cancer, did not impact survival. In summary, we here identify CXCR4 hyperactivation as a co-driver of an aggressive lymphoma phenotype.

DREAM On: Cell Cycle Control in Development and Disease

Perfectly orchestrated periodic gene expression during cell cycle progression is essential for maintaining genome integrity and ensuring that cell proliferation can be stopped by environmental signals. Genetic and proteomic studies during the past two decades revealed remarkable evolutionary conservation of the key mechanisms that control cell cycle-regulated gene expression, including multisubunit DNA-binding DREAM complexes. DREAM complexes containing a retinoblastoma family member, an E2F transcription factor and its dimerization partner, and five proteins related to products of Caenorhabditis elegans multivulva (Muv) class B genes lin-9, lin-37, lin-52, lin-53, and lin-54 (comprising the MuvB core) have been described in diverse organisms, from worms to humans. This review summarizes the current knowledge of the structure, function, and regulation of DREAM complexes in different organisms, as well as the role of DREAM in human disease. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The chromatin-binding domain of Ki-67 together with p53 protects human chromosomes from mitotic damage

Vertebrate mammals express a protein called Ki-67 which is most widely known as a clinically useful marker of highly proliferative cells. Previous studies of human cells indicated that acute depletion of Ki-67 can elicit a delay at the G1/S boundary of the cell cycle, dependent on induction of the checkpoint protein p21. Consistent with those observations, we show here that acute Ki-67 depletion causes hallmarks of DNA damage, and the damage occurs even in the absence of checkpoint signaling. This damage is not observed in cells traversing S phase but is instead robustly detected in mitotic cells. The C-terminal chromatin-binding domain of Ki-67 is necessary and sufficient to protect cells from this damage. We also observe synergistic effects when Ki-67 and p53 are simultaneously depleted, resulting in increased levels of chromosome bridges at anaphase, followed by the appearance of micronuclei. Therefore, these studies identify the C terminus of Ki-67 as an important module for genome stability.

Transcription factor RFX7 governs a tumor suppressor network in response to p53 and stress

Despite its prominence, the mechanisms through which the tumor suppressor p53 regulates most genes remain unclear. Recently, the regulatory factor X 7 (RFX7) emerged as a suppressor of lymphoid neoplasms, but its regulation and target genes mediating tumor suppression remain unknown. Here, we identify a novel p53-RFX7 signaling axis. Integrative analysis of the RFX7 DNA binding landscape and the RFX7-regulated transcriptome in three distinct cell systems reveals that RFX7 directly controls multiple established tumor suppressors, including PDCD4, PIK3IP1, MXD4, and PNRC1, across cell types and is the missing link for their activation in response to p53 and stress. RFX7 target gene expression correlates with cell differentiation and better prognosis in numerous cancer types. Interestingly, we find that RFX7 sensitizes cells to Doxorubicin by promoting apoptosis. Together, our work establishes RFX7’s role as a ubiquitous regulator of cell growth and fate determination and a key node in the p53 transcriptional program.

Tumor suppressor p53: Biology, signaling pathways, and therapeutic targeting

TP53 is the most commonly mutated gene in human cancer with over 100,000 literature citations in PubMed. This is a heavily studied pathway in cancer biology and oncology with a history that dates back to 1979 when p53 was discovered. The p53 pathway is a complex cellular stress response network with multiple diverse inputs and downstream outputs relevant to its role as a tumor suppressor pathway. While inroads have been made in understanding the biology and signaling in the p53 pathway, the p53 family, transcriptional readouts, and effects of an array of mutants, the pathway remains challenging in the realm of clinical translation. While the role of mutant p53 as a prognostic factor is recognized, the therapeutic modulation of its wild-type or mutant activities remain a work-in-progress. This review covers current knowledge about the biology, signaling mechanisms in the p53 pathway and summarizes advances in therapeutic development.

KBoost: a new method to infer gene regulatory networks from gene expression data

Reconstructing gene regulatory networks is crucial to understand biological processes and holds potential for developing personalized treatment. Yet, it is still an open problem as state-of-the-art algorithms are often not able to process large amounts of data within reasonable time. Furthermore, many of the existing methods predict numerous false positives and have limited capabilities to integrate other sources of information, such as previously known interactions. Here we introduce KBoost, an algorithm that uses kernel PCA regression, boosting and Bayesian model averaging for fast and accurate reconstruction of gene regulatory networks. We have benchmarked KBoost against other high performing algorithms using three different datasets. The results show that our method compares favorably to other methods across datasets. We have also applied KBoost to a large cohort of close to 2000 breast cancer patients and 24,000 genes in less than 2 h on standard hardware. Our results show that molecularly defined breast cancer subtypes also feature differences in their GRNs. An implementation of KBoost in the form of an R package is available at: https://github.com/Luisiglm/KBoost and as a Bioconductor software package.

A subset of CB002 xanthine analogs bypass p53-signaling to restore a p53 transcriptome and target an S-phase cell cycle checkpoint in tumors with mutated-p53

Mutations in TP53 occur commonly in the majority of human tumors and confer aggressive tumor phenotypes, including metastasis and therapy resistance. CB002 and structural-analogs restore p53 signaling in tumors with mutant-p53 but we find that unlike other xanthines such as caffeine, pentoxifylline, and theophylline, they do not deregulate the G2 checkpoint. Novel CB002-analogs induce pro-apoptotic Noxa protein in an ATF3/4-dependent manner, whereas caffeine, pentoxifylline, and theophylline do not. By contrast to caffeine, CB002-analogs target an S-phase checkpoint associated with increased p-RPA/RPA2, p-ATR, decreased Cyclin A, p-histone H3 expression, and downregulation of essential proteins in DNA-synthesis and DNA-repair. CB002-analog #4 enhances cell death, and decreases Ki-67 in patient-derived tumor-organoids without toxicity to normal human cells. Preliminary in vivo studies demonstrate anti-tumor efficacy in mice. Thus, a novel class of anti-cancer drugs shows the activation of p53 pathway signaling in tumors with mutated p53, and targets an S-phase checkpoint.

Replication-dependent histone biosynthesis is coupled to cell-cycle commitment

The current model of replication-dependent (RD) histone biosynthesis posits that RD histone gene expression is coupled to DNA replication, occurring only in S phase of the cell cycle once DNA synthesis has begun. However, several key factors in the RD histone biosynthesis pathway are up-regulated by E2F or phosphorylated by CDK2, suggesting these processes may instead begin much earlier, at the point of cell-cycle commitment. In this study, we use both fixed- and live-cell imaging of human cells to address this question, revealing a hybrid model in which RD histone biosynthesis is first initiated in G1, followed by a strong increase in histone production in S phase of the cell cycle. This suggests a mechanism by which cells that have committed to the cell cycle build up an initial small pool of RD histones to be available for the start of DNA replication, before producing most of the necessary histones required in S phase. Thus, a clear distinction exists at completion of mitosis between cells that are born with the intention of proceeding through the cell cycle and replicating their DNA and cells that have chosen to exit the cell cycle and have no immediate need for histone synthesis.

A subset of CB002 xanthine analogs bypass p53-signaling to restore a p53 transcriptome and target an S-phase cell cycle checkpoint in tumors with mutated-p53

Mutations in TP53 occur commonly in the majority of human tumors and confer aggressive tumor phenotypes, including metastasis and therapy resistance. CB002 and structural-analogs restore p53 signaling in tumors with mutant-p53 but we find that unlike other xanthines such as caffeine, pentoxifylline, and theophylline, they do not deregulate the G2 checkpoint. Novel CB002-analogs induce pro-apoptotic Noxa protein in an ATF3/4-dependent manner, whereas caffeine, pentoxifylline, and theophylline do not. By contrast to caffeine, CB002-analogs target an S-phase checkpoint associated with increased p-RPA/RPA2, p-ATR, decreased Cyclin A, p-histone H3 expression, and downregulation of essential proteins in DNA-synthesis and DNA-repair. CB002-analog #4 enhances cell death, and decreases Ki-67 in patient-derived tumor-organoids without toxicity to normal human cells. Preliminary in vivo studies demonstrate anti-tumor efficacy in mice. Thus, a novel class of anti-cancer drugs shows the activation of p53 pathway signaling in tumors with mutated p53, and targets an S-phase checkpoint.

SLC7A11 Is a Superior Determinant of APR-246 (Eprenetapopt) Response than TP53 Mutation Status

APR-246 (eprenetapopt) is in clinical development with a focus on hematologic malignancies and is promoted as a mutant-p53 reactivation therapy. Currently, the detection of at least one TP53 mutation is an inclusion criterion for patient selection into most APR-246 clinical trials. Preliminary results from our phase Ib/II clinical trial investigating APR-246 combined with doublet chemotherapy [cisplatin and 5-fluorouracil (5-FU)] in metastatic esophageal cancer, together with previous preclinical studies, indicate that TP53 mutation status alone may not be a sufficient biomarker for APR-246 response. This study aims to identify a robust biomarker for response to APR-246. Correlation analysis of the PRIMA-1 activity (lead compound to APR-246) with mutational status, gene expression, protein expression, and metabolite abundance across over 700 cancer cell lines (CCL) was performed. Functional validation and a boutique siRNA screen of over 850 redox-related genes were also conducted. TP53 mutation status was not consistently predictive of response to APR-246. The expression of SLC7A11, the cystine/glutamate transporter, was identified as a superior determinant of response to APR-246. Genetic regulators of SLC7A11, including ATF4, MDM2, wild-type p53, and c-Myc, were confirmed to also regulate cancer-cell sensitivity to APR-246. In conclusion, SLC7A11 expression is a broadly applicable determinant of sensitivity to APR-246 across cancer and should be utilized as the key predictive biomarker to stratify patients for future clinical investigation of APR-246.

Homeostatic apoptosis prevents competition-induced atrophy in follicular B cells

While the intrinsic apoptosis pathway is thought to play a central role in shaping the B cell lineage, its precise role in mature B cell homeostasis remains elusive. Using mice in which mature B cells are unable to undergo apoptotic cell death, we show that apoptosis constrains follicular B (FoB) cell lifespan but plays no role in marginal zone B (MZB) cell homeostasis. In these mice, FoB cells accumulate abnormally. This intensifies intercellular competition for BAFF, resulting in a contraction of the MZB cell compartment, and reducing the growth, trafficking, and fitness of FoB cells. Diminished BAFF signaling dampens the non-canonical NF-κB pathway, undermining FoB cell growth despite the concurrent triggering of a protective p53 response. Thus, MZB and FoB cells exhibit a differential requirement for the intrinsic apoptosis pathway. Homeostatic apoptosis constrains the size of the FoB cell compartment, thereby preventing competition-induced FoB cell atrophy.

p63 and p53: Collaborative Partners or Dueling Rivals?

The tumor suppressor p53 and its oncogenic sibling p63 (ΔNp63) direct opposing fates in tumor development. These paralog proteins are transcription factors that elicit their tumor suppressive and oncogenic capacity through the regulation of both shared and unique target genes. Both proteins predominantly function as activators of transcription, leading to a paradigm shift away from ΔNp63 as a dominant negative to p53 activity. The discovery of p53 and p63 as pioneer transcription factors regulating chromatin structure revealed new insights into how these paralogs can both positively and negatively influence each other to direct cell fate. The previous view of a strict rivalry between the siblings needs to be revisited, as p53 and p63 can also work together toward a common goal.

ETV4 and ETV5 drive synovial sarcoma through cell cycle and DUX4 embryonic pathway control

Synovial sarcoma is an aggressive malignancy with no effective treatments for patients with metastasis. The synovial sarcoma fusion SS18-SSX, which recruits the SWI/SNF-BAF chromatin remodeling and polycomb repressive complexes, results in epigenetic activation of FGF receptor (FGFR) signaling. In genetic FGFR-knockout models, culture, and xenograft synovial sarcoma models treated with the FGFR inhibitor BGJ398, we show that FGFR1, FGFR2, and FGFR3 were crucial for tumor growth. Transcriptome analyses of BGJ398-treated cells and histological and expression analyses of mouse and human synovial sarcoma tumors revealed prevalent expression of two ETS factors and FGFR targets, ETV4 and ETV5. We further demonstrate that ETV4 and ETV5 acted as drivers of synovial sarcoma growth, most likely through control of the cell cycle. Upon ETV4 and ETV5 knockdown, we observed a striking upregulation of DUX4 and its transcriptional targets that activate the zygotic genome and drive the atrophy program in facioscapulohumeral dystrophy patients. In addition to demonstrating the importance of inhibiting all three FGFRs, the current findings reveal potential nodes of attack for the cancer with the discovery of ETV4 and ETV5 as appropriate biomarkers and molecular targets, and activation of the embryonic DUX4 pathway as a promising approach to block synovial sarcoma tumors.

Dying to Survive-The p53 Paradox

The p53 tumour suppressor is best known for its canonical role as "guardian of the genome", activating cell cycle arrest and DNA repair in response to DNA damage which, if irreparable or sustained, triggers activation of cell death. However, despite an enormous amount of work identifying the breadth of the gene regulatory networks activated directly and indirectly in response to p53 activation, how p53 activation results in different cell fates in response to different stress signals in homeostasis and in response to p53 activating anti-cancer treatments remains relatively poorly understood. This is likely due to the complex interaction between cell death mechanisms in which p53 has been activated, their neighbouring stressed or unstressed cells and the local stromal and immune microenvironment in which they reside. In this review, we evaluate our understanding of the burgeoning number of cell death pathways affected by p53 activation and how these may paradoxically suppress cell death to ensure tissue integrity and organismal survival. We also discuss how these functions may be advantageous to tumours that maintain wild-type p53, the understanding of which may provide novel opportunity to enhance treatment efficacy.

Progesterone Receptors Promote Quiescence and Ovarian Cancer Cell Phenotypes via DREAM in p53-Mutant Fallopian Tube Models

CONTEXT

The ability of ovarian steroids to modify ovarian cancer (OC) risk remains controversial. Progesterone is considered to be protective; recent studies indicate no effect or enhanced OC risk. Knowledge of progesterone receptor (PR) signaling during altered physiology that typifies OC development is limited.

OBJECTIVE

This study defines PR-driven oncogenic signaling mechanisms in p53-mutant human fallopian tube epithelia (hFTE), a precursor of the most aggressive OC subtype.

METHODS

PR expression in clinical samples of serous tubal intraepithelial carcinoma (STIC) lesions and high-grade serous OC (HGSC) tumors was analyzed. Novel PR-A and PR-B isoform-expressing hFTE models were characterized for gene expression and cell cycle progression, emboli formation, and invasion. PR regulation of the DREAM quiescence complex and DYRK1 kinases was established.

RESULTS

STICs and HGSC express abundant activated phospho-PR. Progestin promoted reversible hFTE cell cycle arrest, spheroid formation, and invasion. RNAseq/biochemical studies revealed potent ligand-independent/-dependent PR actions, progestin-induced regulation of the DREAM quiescence complex, and cell cycle target genes through enhanced complex formation and chromatin recruitment. Disruption of DREAM/DYRK1s by pharmacological inhibition, HPV E6/E7 expression, or DYRK1A/B depletion blocked progestin-induced cell arrest and attenuated PR-driven gene expression and associated OC phenotypes.

CONCLUSION

Activated PRs support quiescence and pro-survival/pro-dissemination cell behaviors that may contribute to early HGSC progression. Our data support an alternative perspective on the tenet that progesterone always confers protection against OC. STICs can reside undetected for decades prior to invasive disease; our studies reveal clinical opportunities to prevent the ultimate development of HGSC by targeting PRs, DREAM, and/or DYRKs.

p53 Frameshift Mutations Couple Loss-of-Function with Unique Neomorphic Activities

p53 mutations that result in loss of transcriptional activity are commonly found in numerous types of cancer. While the majority of these are missense mutations that map within the central DNA-binding domain, truncations and/or frameshift mutations can also occur due to various nucleotide substitutions, insertions, or deletions. These changes result in mRNAs containing premature stop codons that are translated into a diverse group of C-terminally truncated proteins. Here we characterized three p53 frameshift mutant proteins expressed from the endogenous TP53 locus in U2OS osteosarcoma and HCT116 colorectal cancer cell lines. These mutants retain intact DNA-binding domains but display altered oligomerization properties. Despite their abnormally high expression levels, they are mostly transcriptionally inactive and unable to initiate a stimuli-induced transcriptional program characteristic of wild-type p53. However, one of these variant p53 proteins, I332fs*14, which resembles naturally expressed TAp53 isoforms β and γ, retains some residual antiproliferative activity and can induce cellular senescence in HCT116 cells. Cells expressing this mutant also display decreased motility in migration assays. Hence, this p53 variant exhibits a combination of loss-of-gain and gain-of-function characteristics, distinguishing it from both wild type p53 and p53 loss. IMPLICATIONS: p53 frameshift mutants display a mixture of residual antiproliferative and neomorphic functions that may be differentially exploited for targeted therapy.

Dual proteome-scale networks reveal cell-specific remodeling of the human interactome

Thousands of interactions assemble proteins into modules that impart spatial and functional organization to the cellular proteome. Through affinity-purification mass spectrometry, we have created two proteome-scale, cell-line-specific interaction networks. The first, BioPlex 3.0, results from affinity purification of 10,128 human proteins-half the proteome-in 293T cells and includes 118,162 interactions among 14,586 proteins. The second results from 5,522 immunoprecipitations in HCT116 cells. These networks model the interactome whose structure encodes protein function, localization, and complex membership. Comparison across cell lines validates thousands of interactions and reveals extensive customization. Whereas shared interactions reside in core complexes and involve essential proteins, cell-specific interactions link these complexes, "rewiring" subnetworks within each cell's interactome. Interactions covary among proteins of shared function as the proteome remodels to produce each cell's phenotype. Viewable interactively online through BioPlexExplorer, these networks define principles of proteome organization and enable unknown protein characterization.

CK1-mediated phosphorylation of FAM110A promotes its interaction with mitotic spindle and controls chromosomal alignment

Progression through the cell cycle is driven by cyclin-dependent kinases that control gene expression, orchestration of mitotic spindle, and cell division. To identify new regulators of the cell cycle, we performed transcriptomic analysis of human non-transformed cells expressing a fluorescent ubiquitination-based cell cycle indicator and identified 701 transcripts differentially expressed in G1 and G2 cells. Family with sequence similarity 110 member A (FAM110A) protein is highly expressed in G2 cells and localized at mitotic spindle and spindle poles during mitosis. Depletion of FAM110A impairs chromosomal alignment, delays metaphase-to-anaphase transition, and affects spindle positioning. Using mass spectrometry and immunoprecipitation, we identified casein kinase I (CK1) in complex with FAM110A during mitosis. CK1 phosphorylates the C-terminal domain of FAM110A in vitro, and inhibition of CK1 reduces phosphorylation of mitotic FAM110A. Wild-type FAM110A, but not the FAM110A-S252-S255A mutant deficient in CK1 phosphorylation, rescues the chromosomal alignment, duration of mitosis, and orientation of the mitotic spindle after depletion of endogenous FAM110A. We propose that CK1 regulates chromosomal alignment by phosphorylating FAM110A and promoting its interaction with mitotic spindle.

Quantitative imaging of RAD51 expression as a marker of platinum resistance in ovarian cancer

Early relapse after platinum chemotherapy in epithelial ovarian cancer (EOC) portends poor survival. A-priori identification of platinum resistance is therefore crucial to improve on standard first-line carboplatin-paclitaxel treatment. The DNA repair pathway homologous recombination (HR) repairs platinum-induced damage, and the HR recombinase RAD51 is overexpressed in cancer. We therefore designed a REMARK-compliant study of pre-treatment RAD51 expression in EOC, using fluorescent quantitative immunohistochemistry (qIHC) to overcome challenges in quantitation of protein expression in situ. In a discovery cohort (n = 284), RAD51-High tumours had shorter progression-free and overall survival compared to RAD51-Low cases in univariate and multivariate analyses. The association of RAD51 with relapse/survival was validated in a carboplatin monotherapy SCOTROC4 clinical trial cohort (n = 264) and was predominantly noted in HR-proficient cancers (Myriad HRDscore < 42). Interestingly, overexpression of RAD51 modified expression of immune-regulatory pathways in vitro, while RAD51-High tumours showed exclusion of cytotoxic T cells in situ. Our findings highlight RAD51 expression as a determinant of platinum resistance and suggest possible roles for therapy to overcome immune exclusion in RAD51-High EOC. The qIHC approach is generalizable to other proteins with a continuum instead of discrete/bimodal expression.

Identification and functional characterization of new missense SNPs in the coding region of the TP53 gene

Infrequent and rare genetic variants in the human population vastly outnumber common ones. Although they may contribute significantly to the genetic basis of a disease, these seldom-encountered variants may also be miss-identified as pathogenic if no correct references are available. Somatic and germline TP53 variants are associated with multiple neoplastic diseases, and thus have come to serve as a paradigm for genetic analyses in this setting. We searched 14 independent, globally distributed datasets and recovered TP53 SNPs from 202,767 cancer-free individuals. In our analyses, 19 new missense TP53 SNPs, including five novel variants specific to the Asian population, were recurrently identified in multiple datasets. Using a combination of in silico, functional, structural, and genetic approaches, we showed that none of these variants displayed loss of function compared to the normal TP53 gene. In addition, classification using ACMG criteria suggested that they are all benign. Considered together, our data reveal that the TP53 coding region shows far more polymorphism than previously thought and present high ethnic diversity. They furthermore underline the importance of correctly assessing novel variants in all variant-calling pipelines associated with genetic diagnoses for cancer.

E2F1 Maintains Gastric Cancer Stemness Properties by Regulating Stemness-Associated Genes

Purpose

To determine the regulatory role of E2F1 in maintaining gastric cancer stemness properties and the clinical significance of E2F1 in gastric cancer.

Materials and Methods

We conducted a tumor spheroid formation assay to enrich gastric cancer stem-like cells. The protein and mRNA expression levels of genes were measured using Western Blot and qRT-PCR. Lentivirus-mediated overexpression and downregulation of E2F1 were performed to evaluate the effect of E2F1 on the stemness properties of gastric cancer cells. The effect of E2F1 on gastric cancer cell sensitivity of 5-Fu was evaluated using cell viability assay and TdT-mediated dUTP Nick-End Labeling staining. We also analyzed the association between E2F1 expression and clinical characteristics in gastric cancer patients. The KM plotter database was used to analyze the relationship between E2F1 and overall survival in GC patients.

Results

We found that E2F1 expression was significantly higher in gastric cancer tissues than in the paired adjacent normal tissues (p < 0.05) and was positively correlated with tumor size (p < 0.05), T stage (p < 0.05), and differentiation degree (p < 0.05). KM plotter database demonstrated a close association between higher E2F1 expression level and worse overall survival of gastric cancer patients (p < 0.05). In vitro assay illustrated that E2F1 could regulate the expression of stemness-associated genes, such as BMI1, OCT4, Nanog, and CD44, and maintain the tumor spheroid formation ability of gastric cancer cells. E2F1 enhanced 5-Fu resistance in gastric cancer cells, and the E2F1 expression level was correlated with the prognosis of gastric cancer patients receiving 5-Fu therapy. The expression levels of stemness-associated genes were also significantly higher in gastric cancer tissues than the paired adjacent normal tissues (p < 0.05). A positive correlation was observed between E2F1 and BMI1 (r = 0.422, p < 0.05), CD44 (r = 0.634, p < 0.05), OCT4 (r = 0.456, p < 0.05), and Nanog (r = 0.337, p < 0.05) in gastric cancer tissues. The co-overexpression of E2F1 and stemness-associated genes was associated with worse overall survival.

Conclusion

E2F1 plays a significant role in gastric cancer progression by maintaining gastric cancer stemness properties through the regulation of stemness-associated genes. The close association between E2F1 and poor prognosis of patients suggests that E2F1 could serve as a prognostic biomarker and a therapeutic target in gastric cancer patients.

Exploration of immune-related genes in high and low tumor mutation burden groups of chromophobe renal cell carcinoma

Renal cell carcinoma (RCC) is one of most common cancers with gradually increasing incidence and high mortality. Chromogenic RCC (chRCC) is the third most common histological subtype of RCC, accounting for approximately 5–7% of RCC. In our study, the transcriptome expression profile data (n=89) of chRCC, corresponding clinical data (n=113) and the somatic mutation data (n=66) were obtained from the TCGA database. We first analyzed the mutation data of chRCC patients and divided chRCC patients into high and low tumor mutation burden (TMB) groups based on the median TMB. We found that high TMB was significantly associated with worse prognosis and could promote tumor metastasis and development. Moreover, four different immune-related genes (BIRC5, PDGFRL, INHBE, IL20RB) were also identified. We found that BIRC5 was significantly overexpressed in the high TMB group and correlated with worse prognosis. The results of univariate and multivariate COX analyses demonstrated that BIRC5 (hazard ratio (HR) = 2.094) may serve as a prognostic indicator for patients with chRCC with high TMB. In addition, we identified the possible functional pathways of BIRC5 through gene set enrichment analysis (GSEA) enrichment. A positive correlation was obtained between BIRC5 and the abundance of CD4⁺ T cells. The results of our study revealed their correlation between the immune-related genes and clinicopathologic features as well as potential functional pathways as well as immune infiltrating cells, which may provide more data about the development of chRCC immunotherapy.

Mining TCGA database for tumor mutation burden and their clinical significance in bladder cancer

BACKGROUND

Bladder cancer is the ninth most-common cancer worldwide and it is associated with high morbidity and mortality. Tumor mutational burden (TMB) is an emerging biomarker in cancer characterized by microsatellite instability. TMB has been described as a powerful predictor of tumor behavior and response to immunotherapy.

METHODS

A total of 443 bladder cancer samples obtained from The Cancer Genome Atlas (TCGA) were analyzed for mutation types, TMB values, and prognostic value of TMB. Differentially expressed genes (DEGs) were identified from the TMB groupings. Functional analysis was performed to assess the prognostic value of the first 30 core genes. CIBERSORT algorithm was used to determine the correlation between the immune cells and TMB subtypes.

RESULTS

Single nucleotide polymorphism (SNP) and C>T were reported as the most common missense mutations and we also identified a high rate of mutations in TP53, TTN, KMT2D. Bladder cancer patients with high TMB showed a better prognosis. Enrichment analysis of the DEGs revealed that they were involved in the regulation of the P13K-Akt signaling pathway, cytokine-cytokine receptor interaction, and Ras signaling pathway. The high expression of hub genes ADRA2A, CXCL12, S1PR1, ADAMTS9, F13A1, and SPON1 was correlated with poor overall survival. Besides, significant differences in the composition of the immune cells of T cells CD8, T cells CD4 memory activated, NK cells resting and Mast cells resting were observed.

CONCLUSIONS

The present study provides a comprehensive and systematic analysis of the prediction of TMB in bladder cancer and its clinical significance. Also, the study provides additional prognostic information and opportunities for immunotherapy in bladder cancer.

Mechanisms of Cellular Senescence: Cell Cycle Arrest and Senescence Associated Secretory Phenotype

Cellular senescence is a stable cell cycle arrest that can be triggered in normal cells in response to various intrinsic and extrinsic stimuli, as well as developmental signals. Senescence is considered to be a highly dynamic, multi-step process, during which the properties of senescent cells continuously evolve and diversify in a context dependent manner. It is associated with multiple cellular and molecular changes and distinct phenotypic alterations, including a stable proliferation arrest unresponsive to mitogenic stimuli. Senescent cells remain viable, have alterations in metabolic activity and undergo dramatic changes in gene expression and develop a complex senescence-associated secretory phenotype. Cellular senescence can compromise tissue repair and regeneration, thereby contributing toward aging. Removal of senescent cells can attenuate age-related tissue dysfunction and extend health span. Senescence can also act as a potent anti-tumor mechanism, by preventing proliferation of potentially cancerous cells. It is a cellular program which acts as a double-edged sword, with both beneficial and detrimental effects on the health of the organism, and considered to be an example of evolutionary antagonistic pleiotropy. Activation of the p53/p21^WAF1/CIP1 and p16^INK4A/pRB tumor suppressor pathways play a central role in regulating senescence. Several other pathways have recently been implicated in mediating senescence and the senescent phenotype. Herein we review the molecular mechanisms that underlie cellular senescence and the senescence associated growth arrest with a particular focus on why cells stop dividing, the stability of the growth arrest, the hypersecretory phenotype and how the different pathways are all integrated.

Integrative analysis regarding the correlation between GAS2 family genes and human glioma prognosis

Background

Emerging oncogenes were reportedly linked to the complicated subtypes and pathogenesis of clinical gliomas. Herein, we first comprehensively explored the potential correlation between growth‐arrest‐specific two family genes (GAS2, GAS2L1, GAS2L2, GAS2L3) and gliomas by bioinformatics analysis and cellular experiments.

Methods

Based on the available datasets of TCGA (The Cancer Genome Atlas), CGGA (Chinese Glioma Genome Atlas), and Oncomine databases, we performed a series of analyses, such as gene expression, survival prognosis, DNA methylation, immune infiltration, and partner enrichment. We also utilized two glioma cell lines to conduct the colony formation and wound‐healing assay.

Results

GAS2L3 gene was highly expressed in glioma tissues compared to normal brain tissues (p < 0.05). We further observed the relationship between the high expressed GAS2L3 and poor clinical prognosis of brain low‐grade glioma (LGG) cases in our Cox proportional hazard model (hazard ratio [HR] = 0.1715, p < 0.001). Moreover, DNA hypomethylation status of GAS2L3 was correlated with the high expression of GAS2L3 in LGG tissues and the poor clinical prognosis of primary glioma cases (p < 0.05). We also found that the high expression of GAS2L3 was associated with the infiltration level of immune cells, especially the T cells (p < 0.0001). Functional enrichment analysis of GAS2L3‐correlated genes and interaction partners further indicated that GAS2L3 might take part in the occurrence of glioma by influencing a series of biological behaviors, such as cell division, cytoskeleton binding, and cell adhesion. Additionally, our cellular experiment data suggested that a highly expressed GAS2L3 gene contributes to the enhanced proliferation and migration of glioma cells.

Conclusion

This study first analyzed the potential role of GAS2 family genes, especially GAS2L3, in the clinical prognosis and possible functional mechanisms of glioma, which gives a novel insight into the relationship between GAS2L3 and LGG.

Interpreting Sequence Variation in PDAC-Predisposing Genes Using a Multi-Tier Annotation Approach Performed at the Gene, Patient, and Cohort Level

We investigated germline variation in pancreatic ductal adenocarcinoma (PDAC) predisposition genes in 535 patients, using a custom-built panel and a new complementary bioinformatic approach. Our panel assessed genes belonging to DNA repair, cell cycle checkpoints, migration, and preneoplastic pancreatic conditions. Our bioinformatics approach integrated annotations of variants by using data derived from both germline and somatic references. This integrated approach with expanded evidence enabled us to consider patterns even among private mutations, supporting a functional role for certain alleles, which we believe enhances individualized medicine beyond classic gene-centric approaches. Concurrent evaluation of three levels of evidence, at the gene, sample, and cohort level, has not been previously done. Overall, we identified in PDAC patient germline samples, 12% with mutations previously observed in pancreatic cancers, 23% with mutations previously discovered by sequencing other human tumors, and 46% with mutations with germline associations to cancer. Non-polymorphic protein-coding pathogenic variants were found in 18.4% of patient samples. Moreover, among patients with metastatic PDAC, 16% carried at least one pathogenic variant, and this subgroup was found to have an improved overall survival (22.0 months versus 9.8; p=0.008) despite a higher pre-treatment CA19-9 level (p=0.02). Genetic alterations in DNA damage repair genes were associated with longer overall survival among patients who underwent resection surgery (92 months vs. 46; p=0.06). ATM alterations were associated with more frequent metastatic stage (p = 0.04) while patients with BRCA1 or BRCA2 alterations had improved overall survival (79 months vs. 39; p=0.05). We found that mutations in genes associated with chronic pancreatitis were more common in non-white patients (p<0.001) and associated with longer overall survival (52 months vs. 26; p=0.004), indicating the need for greater study of the relationship among these factors. More than 90% of patients were found to have variants of uncertain significance, which is higher than previously reported. Furthermore, we generated 3D models for selected mutant proteins, which suggested distinct mechanisms underlying their dysfunction, likely caused by genetic alterations. Notably, this type of information is not predictable from sequence alone, underscoring the value of structural bioinformatics to improve genomic interpretation. In conclusion, the variation in PDAC predisposition genes appears to be more extensive than anticipated. This information adds to the growing body of literature on the genomic landscape of PDAC and brings us closer to a more widespread use of precision medicine for this challenging disease.

Oncogenic B-Myb Is Associated With Deregulation of the DREAM-Mediated Cell Cycle Gene Expression Program in High Grade Serous Ovarian Carcinoma Clinical Tumor Samples

Cell cycle control drives cancer progression and treatment response in high grade serous ovarian carcinoma (HGSOC). MYBL2 (encoding B-Myb), an oncogene with prognostic significance in several cancers, is highly expressed in most HGSOC cases; however, the clinical significance of B-Myb in this disease has not been well-characterized. B-Myb is associated with cell proliferation through formation of the MMB (Myb and MuvB core) protein complex required for transcription of mitotic genes. High B-Myb expression disrupts the formation of another transcriptional cell cycle regulatory complex involving the MuvB core, DREAM (DP, RB-like, E2F, and MuvB), in human cell lines. DREAM coordinates cell cycle dependent gene expression by repressing over 800 cell cycle genes in G0/G1. Here, we take a bioinformatics approach to further evaluate the effect of B-Myb expression on DREAM target genes in HGSOC and validate our cellular model with clinical specimens. We show that MYBL2 is highly expressed in HGSOC and correlates with expression of DREAM and MMB target genes in both The Cancer Genome Atlas (TCGA) as well as independent analyses of HGSOC primary tumors (N = 52). High B-Myb expression was also associated with poor overall survival in the TCGA cohort and analysis by a DREAM target gene expression signature yielded a negative impact on survival. Together, our data support the conclusion that high expression of MYBL2 is associated with deregulation of DREAM/MMB-mediated cell cycle gene expression programs in HGSOC and may serve as a prognostic factor independent of its cell cycle role. This provides rationale for further, larger scale studies aimed to determine the clinical predictive value of the B-Myb gene expression signature for treatment response as well as patient outcomes.

FAM72, Glioblastoma Multiforme (GBM) and Beyond

Simple Summary

Glioblastoma multiforme (GBM) is a serious and aggressive cancer disease that has not allowed scientists to rest for decades. In this review, we consider the new gene pair |-SRGAP2–FAM72-| and discuss its role in the cell cycle and the possibility of defining new therapeutic approaches for the treatment of GBM and other cancers via this gene pair |-SRGAP2–FAM72-|.

Abstract

Neural stem cells (NSCs) offer great potential for regenerative medicine due to their excellent ability to differentiate into various specialized cell types of the brain. In the central nervous system (CNS), NSC renewal and differentiation are under strict control by the regulation of the pivotal SLIT-ROBO Rho GTPase activating protein 2 (SRGAP2)—Family with sequence similarity 72 (FAM72) master gene (i.e., |-SRGAP2–FAM72-|) via a divergent gene transcription activation mechanism. If the gene transcription control unit (i.e., the intergenic region of the two sub-gene units, SRGAP2 and FAM72) gets out of control, NSCs may transform into cancer stem cells and generate brain tumor cells responsible for brain cancer such as glioblastoma multiforme (GBM). Here, we discuss the surveillance of this |-SRGAP2–FAM72-| master gene and its role in GBM, and also in light of FAM72 for diagnosing various types of cancers outside of the CNS.

MMB-FOXM1-driven premature mitosis is required for CHK1 inhibitor sensitivity

To identify genes whose loss confers resistance to CHK1 inhibitors, we perform genome-wide CRISPR-Cas9 screens in non-small-cell lung cancer (NSCLC) cell lines treated with the CHK1 inhibitor prexasertib (CHK1i). Five of the top six hits of the screens, MYBL2 (B-MYB), LIN54, FOXM1, cyclin A2 (CCNA2), and CDC25B, are cell-cycle-regulated genes that contribute to entry into mitosis. Knockout of MMB-FOXM1 complex components LIN54 and FOXM1 reduce CHK1i-induced DNA replication stress markers and premature mitosis during Late S phase. Activation of a feedback loop between the MMB-FOXM1 complex and CDK1 is required for CHK1i-induced premature mitosis in Late S phase and subsequent replication catastrophe, indicating that dysregulation of the S to M transition is necessary for CHK1 inhibitor sensitivity. These findings provide mechanistic insights into small molecule inhibitors currently studied in clinical trials and provide rationale for combination therapies.

Branigan et al., by using genome-wide CRISPR screens, identify the MMB-FOXM1 complex as being required for CHK1 inhibitor (CHK1i) sensitivity. Their study shows that CHK1i-induced premature activation of the G2/M transcriptional program by this complex triggers a breakdown in the separation of DNA synthesis and mitosis, leading to replication catastrophe.

PAF remodels the DREAM complex to bypass cell quiescence and promote lung tumorigenesis

The DREAM complex orchestrates cell quiescence and the cell cycle. However, how the DREAM complex is deregulated in cancer remains elusive. Here, we report that PAF (PCLAF/KIAA0101) drives cell quiescence exit to promote lung tumorigenesis by remodeling the DREAM complex. PAF is highly expressed in lung adenocarcinoma (LUAD) and is associated with poor prognosis. Importantly, Paf knockout markedly suppressed LUAD development in mouse models. PAF depletion induced LUAD cell quiescence and growth arrest. PAF is required for the global expression of cell-cycle genes controlled by the repressive DREAM complex. Mechanistically, PAF inhibits DREAM complex formation by binding to RBBP4, a core DREAM subunit, leading to transactivation of DREAM target genes. Furthermore, pharmacological mimicking of PAF-depleted transcriptomes inhibited LUAD tumor growth. Our results unveil how the PAF-remodeled DREAM complex bypasses cell quiescence to promote lung tumorigenesis and suggest that the PAF-DREAM axis may be a therapeutic vulnerability in lung cancer.

Identification of an immune-related RNA-binding protein signature to predict survival and targeted therapy responses in liver cancer

RNA-binding proteins (RBPs) play crucial roles in multiple cancers. However, very few RBPs and their association with immune genes have been systematically studied in liver cancer (LC). We aimed to identify an immune-related RBP signature to predict the survival of LC patients. Bioinformatics methods were used to identify differentially expressed, immune-related, and prognostic RBPs and to develop an immune-related RBP signature based on data from the Cancer Genome Atlas (TCGA) cohort. We obtained eight differentially expressed, immune-related, and prognostic RBPs to construct a risk signature. The signature could effectively distinguish between high- and low-risk patients, and its predictive capacity was validated in the International Cancer Genomics Consortium (ICGC) cohort. We speculated that the high-risk group was more sensitive to targeted therapy. The immune-related RBP signature is an independent prognostic biomarker for LC patients and can expand the application of targeted therapy through patient stratification.

Restoring the DREAM Complex Inhibits the Proliferation of High-Risk HPV Positive Human Cells

High-risk (HR) human papillomaviruses are known causative agents in 5% of human cancers including cervical, ano-genital and head and neck carcinomas. In part, HR-HPV causes cancer by targeting host-cell tumor suppressors including retinoblastoma protein (pRb) and RB-like proteins p107 and p130. HR-HPV E7 uses a LxCxE motif to bind RB proteins, impairing their ability to control cell-cycle dependent transcription. E7 disrupts DREAM (Dimerization partner, RB-like, E2F and MuvB), a transcriptional repressor complex that can include p130 or p107, but not pRb, which regulates genes required for cell cycle progression. However, it is not known whether disruption of DREAM plays a significant role in HPV-driven tumorigenesis. In the DREAM complex, LIN52 is an adaptor that binds directly to p130 via an E7-like LxSxE motif. Replacement of the LxSxE sequence in LIN52 with LxCxE (LIN52-S20C) increases p130 binding and partially restores DREAM assembly in HPV-positive keratinocytes and human cervical cancer cells, inhibiting proliferation. Our findings demonstrate that disruption of the DREAM complex by E7 is an important process promoting cellular proliferation by HR-HPV. Restoration of the DREAM complex in HR-HPV positive cells may therefore have therapeutic benefits in HR-HPV positive cancers.

Sirtuin 5 Is Regulated by the SCFCyclin F Ubiquitin Ligase and Is Involved in Cell Cycle Control

The ubiquitin-proteasome system is essential for cell cycle progression. Cyclin F is a cell cycle-regulated substrate adapter F-box protein for the Skp1, CUL1, and F-box protein (SCF) family of E3 ubiquitin ligases. Despite its importance in cell cycle progression, identifying cyclin F-bound SCF complex (SCF^{Cyclin F}) substrates has remained challenging. Since cyclin F overexpression rescues a yeast mutant in the cdc4 gene, we considered the possibility that other genes that genetically modify cdc4 mutant lethality could also encode cyclin F substrates. We identified the mitochondrial and cytosolic deacylating enzyme sirtuin 5 (SIRT5) as a novel cyclin F substrate. SIRT5 has been implicated in metabolic processes, but its connection to the cell cycle is not known. We show that cyclin F interacts with and controls the ubiquitination, abundance, and stability of SIRT5. We show SIRT5 knockout results in a diminished G₁ population and a subsequent increase in both S and G₂/M. Global proteomic analyses reveal cyclin-dependent kinase (CDK) signaling changes congruent with the cell cycle changes in SIRT5 knockout cells. Together, these data demonstrate that SIRT5 is regulated by cyclin F and suggest a connection between SIRT5, cell cycle regulation, and metabolism.

MYBL2-Driven Transcriptional Programs Link Replication Stress and Error-prone DNA Repair With Genomic Instability in Lung Adenocarcinoma

It has long been recognized that defects in cell cycle checkpoint and DNA repair pathways give rise to genomic instability, tumor heterogeneity, and metastasis. Despite this knowledge, the transcription factor-mediated gene expression programs that enable survival and proliferation in the face of enormous replication stress and DNA damage have remained elusive. Using robust omics data from two independent studies, we provide evidence that a large cohort of lung adenocarcinomas exhibit significant genome instability and overexpress the DNA damage responsive transcription factor MYB proto-oncogene like 2 (MYBL2). Across two studies, elevated MYBL2 expression was a robust marker of poor overall survival and disease-free survival outcomes, regardless of disease stage. Clinically, elevated MYBL2 expression identified patients with aggressive early onset disease, increased lymph node involvement, and increased incidence of distant metastases. Analysis of genomic sequencing data demonstrated that MYBL2 High lung adenocarcinomas had elevated somatic mutation burden, widespread chromosomal alterations, and alterations in single-strand DNA break repair pathways. In this study, we provide evidence that impaired single-strand break repair, combined with a loss of cell cycle regulators TP53 and RB1, give rise to MYBL2-mediated transcriptional programs. Omics data supports a model wherein tumors with significant genomic instability upregulate MYBL2 to drive genes that control replication stress responses, promote error-prone DNA repair, and antagonize faithful homologous recombination repair. Our study supports the use of checkpoint kinase 1 (CHK1) pharmacological inhibitors, in targeted MYBL2 High patient cohorts, as a future therapy to improve lung adenocarcinoma patient outcomes.

Tumor Mutation Burden, Immune Cell Infiltration, and Construction of Immune-Related Genes Prognostic Model in Head and Neck Cancer

Background: Head and neck squamous cell carcinoma (HNSCC) is the sixth most common malignancy worldwide, and the prognosis of HNSCC remains bleak. Numerous studies revealed that the tumor mutation burden (TMB) could predict the survival outcomes of a variety of tumors. Objectives: This study aimed to investigate the TMB and immune cell infiltration in these patients and construct an immune-related genes (IRGs) prognostic model. Methods: The expression data of 546 HNSCC patients were obtained from The Cancer Genome Atlas (TCGA) database. All patients were divided into high- and low- TMB groups, and the relationship between TMB and clinical relevance was further analyzed. The differentially expressed genes (DEGs) were identified using the R software package, limma. Functional enrichment analyses were conducted to identify the significantly enriched pathways between two groups. CIBERSORT algorithm was adopted to calculate the abundance of 22 leukocyte subtypes. The IRGs prognostic model was constructed via the multivariate Cox regression analysis. Results: Missense mutation and single nucleotide variants (SNV) were the most predominant mutation types in HNSCC. TP53, TTN, and FAT1 were the most frequently mutated genes. Patients with high TMB were observed with worse survival outcomes. The functional analysis of TMB associated DEGs showed that the identified DEGs mainly involved in spliceosome, RNA degradation, proteasome, and RNA polymerase pathways. We observed that macrophages, T cells CD8, and T cells CD4 memory were the most commonly infiltrated subtypes of immune cells in HNSCC. Finally, an IRGs prognostic model was constructed, and the AUC of the ROC curve was 0.635. Conclusions: Our results suggest that high TMB is associated with poor prognosis in HNSCC patients. The constructed model has potential prognostic value for the prognosis of these individuals, and it needs to be further validated in large-scale and prospective studies.

Benzo[a]pyrene represses DNA repair through altered E2F1/E2F4 function marking an early event in DNA damage-induced cellular senescence

Transcriptional regulation of DNA repair is of outmost importance for the restoration of DNA integrity upon genotoxic stress. Here we report that the potent environmental carcinogen benzo[a]pyrene (B[a]P) activates a cellular DNA damage response resulting in transcriptional repression of mismatch repair (MMR) genes (MSH2, MSH6, EXO1) and of RAD51, the central homologous recombination repair (HR) component, ultimately leading to downregulation of MMR and HR. B[a]P-induced gene repression is caused by abrogated E2F1 signalling. This occurs through proteasomal degradation of E2F1 in G2-arrested cells and downregulation of E2F1 mRNA expression in G1-arrested cells. Repression of E2F1-mediated transcription and silencing of repair genes is further mediated by the p21-dependent E2F4/DREAM complex. Notably, repression of DNA repair is also observed following exposure to the active B[a]P metabolite BPDE and upon ionizing radiation and occurs in response to a p53/p21-triggered, irreversible cell cycle arrest marking the onset of cellular senescence. Overall, our results suggest that repression of MMR and HR is an early event during genotoxic-stress induced senescence. We propose that persistent downregulation of DNA repair might play a role in the maintenance of the senescence phenotype, which is associated with an accumulation of unrepairable DNA lesions.

High Risk α-HPV E6 Impairs Translesion Synthesis by Blocking POLη Induction

Simple Summary

Cervical cancers (CaCx) are caused by the expression of human papillomavirus oncogenes (HPV E6 and E7). Here, in vitro assays, computational approaches and immunohistochemical analysis of cervical biopsies show that HPV oncogenes impair translesion synthesis (TLS). This limits the pathway’s ability to prevent replication stress from causing fork collapse and DNA damage. As a result, HPV oncogenes make cells more sensitive to replication stressing agents, such as Cisplatin. Mechanistically, HPV E6 prevents replication stress from triggering the accumulation of a TLS-specific polymerase (POLη). Supplying exogenous POLη to CaCx cells rescues TLS and lowers Cisplatin toxicity.

Abstract

High risk genus α human papillomaviruses (α-HPVs) express two versatile oncogenes (α-HPV E6 and E7) that cause cervical cancer (CaCx) by degrading tumor suppressor proteins (p53 and RB). α-HPV E7 also promotes replication stress and alters DNA damage responses (DDR). The translesion synthesis pathway (TLS) mitigates DNA damage by preventing replication stress from causing replication fork collapse. Computational analysis of gene expression in CaCx transcriptomic datasets identified a frequent increased expression of TLS genes. However, the essential TLS polymerases did not follow this pattern. These data were confirmed with in vitro and ex vivo systems. Further interrogation of TLS, using POLη as a representative TLS polymerase, demonstrated that α-HPV16 E6 blocks TLS polymerase induction by degrading p53. This doomed the pathway, leading to increased replication fork collapse and sensitivity to treatments that cause replication stress (e.g., UV and Cisplatin). This sensitivity could be overcome by the addition of exogenous POLη.

Dissecting the DNA binding landscape and gene regulatory network of p63 and p53

The transcription factor p53 is the best-known tumor suppressor, but its sibling p63 is a master regulator of epidermis development and a key oncogenic driver in squamous cell carcinomas (SCC). Despite multiple gene expression studies becoming available, the limited overlap of reported p63-dependent genes has made it difficult to decipher the p63 gene regulatory network. Particularly, analyses of p63 response elements differed substantially among the studies. To address this intricate data situation, we provide an integrated resource that enables assessing the p63-dependent regulation of any human gene of interest. We use a novel iterative de novo motif search approach in conjunction with extensive ChIP-seq data to achieve a precise global distinction between p53-and p63-binding sites, recognition motifs, and potential co-factors. We integrate these data with enhancer:gene associations to predict p63 target genes and identify those that are commonly de-regulated in SCC representing candidates for prognosis and therapeutic interventions.

FOXM1 drives HPV+ HNSCC sensitivity to WEE1 inhibition

Head and neck squamous cell carcinoma (HNSCC) associated with high-risk human papilloma virus (HPV) infection is a growing clinical problem. The WEE1 kinase inhibitor AZD1775 (WEE1i) overrides cell cycle checkpoints and is being studied in HNSCC regimens. We show that the HPV16 E6/E7 oncoproteins sensitize HNSCC cells to single-agent WEE1i treatment through activation of a FOXM1-CDK1 circuit that drives mitotic gene expression and DNA damage. An isogenic cell system indicated that E6 largely accounts for these phenotypes in ways that extend beyond p53 inactivation. A targeted genomic analysis implicated FOXM1 signaling downstream of E6/E7 expression and analyses of primary tumors and The Cancer Genome Atlas (TCGA) data revealed an activated FOXM1-directed promitotic transcriptional signature in HPV+ versus HPV- HNSCCs. Finally, we demonstrate the causality of FOXM1 in driving WEE1i sensitivity. These data suggest that elevated basal FOXM1 activity predisposes HPV+ HNSCC to WEE1i-induced toxicity and provide mechanistic insights into WEE1i and HPV+ HNSCC therapies.

P63 modulates the expression of the WDFY2 gene which is implicated in cancer regulation and limb development

TP63 is a member of the TP53 gene family, sharing a common gene structure that produces two groups of mRNAs' encoding proteins with different N-terminal regions (ΔN and TA isoforms); both transcripts are also subjected to alternative splicing mechanisms at C-terminus, generating a variety of isoforms. p63 is a master regulator of epidermal development and homoeostasis as well as an important player in tumorigenesis and cancer progression with both oncogenic and tumour suppressive roles. A number of studies have aimed at the identification of p63 target genes, allowing the dissection of the molecular pathways orchestrated by the different isoforms. In the present study we investigated in more detail the p63 responsiveness of the WDFY2 (WD repeat and FYVE domain containing 2) gene, encoding for an endosomal protein identified as a binding partner of the PI-3K/AKT signalling pathway. We showed that overexpression of different p63 isoforms was able to induce WDFY2 expression in TP53-null cells. The p63-dependent transcriptional activation was associated with specific response elements (REs) that have been identified by a bioinformatics tool and validated by yeast- and mammal-based assays. Interestingly, to confirm that WDFY2 belongs to the p63 network of cancer regulation, we analysed the impact of WDFY2 alterations, by showing its frequent deletion in different types of tumours and suggesting its expression level as a prognostic biomarker. Lastly, we identified a chromosomal translocation involving the WDFY2 locus in a patient affected by a rare congenital limb anomaly, indicating WDFY2 as a possible susceptibility gene placed downstream p63 in the network of limb development.

Tumor suppressor p53: from engaging DNA to target gene regulation

The p53 transcription factor confers its potent tumor suppressor functions primarily through the regulation of a large network of target genes. The recent explosion of next generation sequencing protocols has enabled the study of the p53 gene regulatory network (GRN) and underlying mechanisms at an unprecedented depth and scale, helping us to understand precisely how p53 controls gene regulation. Here, we discuss our current understanding of where and how p53 binds to DNA and chromatin, its pioneer-like role, and how this affects gene regulation. We provide an overview of the p53 GRN and the direct and indirect mechanisms through which p53 affects gene regulation. In particular, we focus on delineating the ubiquitous and cell type-specific network of regulatory elements that p53 engages; reviewing our understanding of how, where, and when p53 binds to DNA and the mechanisms through which these events regulate transcription. Finally, we discuss the evolution of the p53 GRN and how recent work has revealed remarkable differences between vertebrates, which are of particular importance to cancer researchers using mouse models.

Clinical significance of expression of PPP1R105 in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is a common malignant tumor with a poor prognosis. In this study, we evaluated the expression of the PPP1R105 gene in HCC and its relationship with prognosis through bioinformatics analysis. Then, we used immunohistochemistry to verify the expression of PPP1R105 protein.

AIM

To investigate the expression of PPP1R105 (Ki-67 or MKI67) in HCC and its clinical significance.

METHODS

The relative expression data of the PPP1R105 gene in various tumor tissues were retrieved from the tumor protein database and TCGA database and compared between cancer and paired normal tissues. The protein-protein interaction network (PPI) of PPP1R105 was constructed by using the STRING database. The progression free survival (DFS) and overall survival (OS) between the high and low PPP1R105 expression groups were compared. Eighty-one HCC patients who underwent surgery were retrospectively included in the study, and PPP1R105 protein expression was detected by immunohistochemistry. The relationship between PPP1R105 protein expression and patients' features was evaluated.

RESULTS

The expression levels of PPP1R105 mRNA in different tumors were not significantly different. In HCC, the expression level of PPP1R105 mRNA was significantly higher than that of normal liver tissue, and the expression level was related with clinical stage. The aggregation index of PPI was 0.99, and the enrichment of protein interaction network was significant (P < 0.01). PPP1R105 was mainly enriched in cell cycle, cell aging, virus carcinogenesis, and p53 signaling pathway. The positive correlation between KIF18B and PPP1R105 was the most significant (r_person = 0.91, P < 0.05), while the negative correlation between MCELL and PPP1R105 was the most significant (r_person = -0.59, P < 0.05). DFS (HR = 1.90, P < 0.01) and OS (HR = 1.90, P < 0.01) in the high PPP1R105 expression group were significantly shorter than those in the low expression group. PPP1R105 protein was mainly expressed in the cytoplasm and nucleolus of the cells. Of the included 81 patients, 30 (37.0%) showed high expression of PPP1R105 protein. The proportion of patients with high expression of PPP1R105 protein was higher than that of patients with low expression of MIK67 protein (P < 0.05).

CONCLUSION

PPP1R105 is highly expressed in patients with HCC, which is related to a poor DFS and OS. PPP1R105 can be used as a poor prognostic molecular marker in HCC.

The differences of regulatory networks between papillary and anaplastic thyroid carcinoma: an integrative transcriptomics study

BACKGROUND

Unlike papillary thyroid cancer (PTC), anaplastic thyroid carcinoma (ATC) is extremely aggressive and rapidly lethal without effective therapies. However, the differences of master regulators and regulatory networks between PTC and ATC remain unclear. Methods: Three representative datasets comprising 32 ATC, 69 PTC, and 78 normal thyroid tissue samples were combined to form a large dataset. Differentially expressed genes (DEGs) were identified and enriched by limma package and gene set enrichment analysis, respectively. Subsequently, protein-protein interaction network and transcription factors (TFs) regulatory network were constructed to identify gene modules and master regulators. Further, master regulators were validated by RT-PCR and western blot. Finally, Kaplan-Meier plotter was applied to evaluate their prognostic values. Results: A total of 560 DEGs were identified as ATC-specific malignant signature. The regulatory network analysis showed that nine master regulators were significantly correlated with three gene modules and potentially regulated the expression of DEGs in three gene modules, respectively. Furthermore, CREB3L1, FOSL2, E2F1 and CAT were significantly associated with overall survival of thyroid cancer patients. FOXM1, FOSL2, MYBL2, AVEN and E2F1 were unfavorable factors of recurrence-free survival (RFS), while CAT was a favorable factor of RFS. RT-PCR and western blot confirmed that six TFs were obviously up-regulated in ATC tissues/cell line as compared with PTC and normal thyroid tissues/cell lines, respectively. In addition, 19 ATC-specific kinases were identified to illustrate the potential post-translational modification. Conclusions: Our findings provide a comprehensive insight into malignant mechanism of ATC, which may indicate their value in the future investigation of ATC.

The Myb-MuvB Complex Is Required for YAP-Dependent Transcription of Mitotic Genes

YAP and TAZ, downstream effectors of the Hippo pathway, are important regulators of proliferation. Here, we show that the ability of YAP to activate mitotic gene expression is dependent on the Myb-MuvB (MMB) complex, a master regulator of genes expressed in the G2/M phase of the cell cycle. By carrying out genome-wide expression and binding analyses, we found that YAP promotes binding of the MMB subunit B-MYB to the promoters of mitotic target genes. YAP binds to B-MYB and stimulates B-MYB chromatin association through distal enhancer elements that interact with MMB-regulated promoters through chromatin looping. The cooperation between YAP and B-MYB is critical for YAP-mediated entry into mitosis. Furthermore, the expression of genes coactivated by YAP and B-MYB is associated with poor survival of cancer patients. Our findings provide a molecular mechanism by which YAP and MMB regulate mitotic gene expression and suggest a link between two cancer-relevant signaling pathways.

Complex Cartography: Regulation of E2F Transcription Factors by Cyclin F and Ubiquitin

The E2F family of transcriptional regulators sits at the center of cell cycle gene expression and plays vital roles in normal and cancer cell cycles. Whereas control of E2Fs by the retinoblastoma family of proteins is well established, much less is known about their regulation by ubiquitin pathways. Recent studies placed the Skp1-Cul1-F-box-protein (SCF) family of E3 ubiquitin ligases with the F-box protein Cyclin F at the center of E2F regulation, demonstrating temporal proteolysis of both activator and atypical repressor E2Fs. Importantly, these E2F members, in particular activator E2F1 and repressors E2F7 and E2F8, form a feedback circuit at the crossroads of cell cycle and cell death. Moreover, Cyclin F functions in a reciprocal circuit with the cell cycle E3 ligase anaphase-promoting complex/cyclosome (APC/C), which also controls E2F7 and E2F8. This review focuses on the complex contours of feedback within this circuit, highlighting the deep crosstalk between E2F, SCF-Cyclin F, and APC/C in regulating the oscillator underlying human cell cycles.

The pseudo-caspase FLIP(L) regulates cell fate following p53 activation

p53 is the most frequently mutated, well-studied tumor-suppressor gene, yet the molecular basis of the switch from p53-induced cell-cycle arrest to apoptosis remains poorly understood. Using a combination of transcriptomics and functional genomics, we unexpectedly identified a nodal role for the caspase-8 paralog and only human pseudo-caspase, FLIP(L), in regulating this switch. Moreover, we identify FLIP(L) as a direct p53 transcriptional target gene that is rapidly up-regulated in response to Nutlin-3A, an MDM2 inhibitor that potently activates p53. Genetically or pharmacologically inhibiting expression of FLIP(L) using siRNA or entinostat (a clinically relevant class-I HDAC inhibitor) efficiently promoted apoptosis in colorectal cancer cells in response to Nutlin-3A, which otherwise predominantly induced cell-cycle arrest. Enhanced apoptosis was also observed when entinostat was combined with clinically relevant, p53-activating chemotherapy in vitro, and this translated into enhanced in vivo efficacy. Mechanistically, FLIP(L) inhibited p53-induced apoptosis by blocking activation of caspase-8 by the TRAIL-R2/DR5 death receptor; notably, this activation was not dependent on receptor engagement by its ligand, TRAIL. In the absence of caspase-8, another of its paralogs, caspase-10 (also transcriptionally up-regulated by p53), induced apoptosis in Nutlin-3A-treated, FLIP(L)-depleted cells, albeit to a lesser extent than in caspase-8-proficient cells. FLIP(L) depletion also modulated transcription of canonical p53 target genes, suppressing p53-induced expression of the cell-cycle regulator p21 and enhancing p53-induced up-regulation of proapoptotic PUMA. Thus, even in the absence of caspase-8/10, FLIP(L) silencing promoted p53-induced apoptosis by enhancing PUMA expression. Thus, we report unexpected, therapeutically relevant roles for FLIP(L) in determining cell fate following p53 activation.

Promiscuous Targeting of Cellular Proteins by Vpr Drives Systems-Level Proteomic Remodeling in HIV-1 Infection

HIV-1 encodes four “accessory proteins” (Vif, Vpr, Vpu, and Nef), dispensable for viral replication in vitro but essential for viral pathogenesis in vivo. Well characterized cellular targets have been associated with Vif, Vpu, and Nef, which counteract host restriction and promote viral replication. Conversely, although several substrates of Vpr have been described, their biological significance remains unclear. Here, we use complementary unbiased mass spectrometry-based approaches to demonstrate that Vpr is both necessary and sufficient for the DCAF1/DDB1/CUL4 E3 ubiquitin ligase-mediated degradation of at least 38 cellular proteins, causing systems-level changes to the cellular proteome. We therefore propose that promiscuous targeting of multiple host factors underpins complex Vpr-dependent cellular phenotypes and validate this in the case of G2/M cell cycle arrest. Our model explains how Vpr modulates so many cell biological processes and why the functional consequences of previously described Vpr targets, identified and studied in isolation, have proved elusive.

•

HIV-1 Vpr is responsible for almost all proteomic changes in HIV-1-infected cells

•

Vpr directly targets multiple nuclear proteins for degradation

•

Vpr cellular phenotypes (e.g., cell cycle arrest) stem from broad substrate targeting

•

Targeting of a few proteins is conserved across diverse primate lentiviral species

Exploration of gene expression profiles and immune microenvironment between high and low tumor mutation burden groups in prostate cancer

BACKGROUND

Tumor mutation burden (TMB) has been established as a biomarker for response to immune therapy and prognosis in various cancers. However, the association between TMB and prognosis of prostate cancer (PCa) remains unclear. This study aimed to investigate the impact of TMB in biochemical recurrence (BCR) and the immune microenvironment in high and low TMB groups.

METHODS

Mutation data, gene expression, clinicopathological information were downloaded from The Cancer Genome Atlas (TCGA). Mutation types and TMB values were identified. All samples were divided into high and low TMB groups with median TMB value as the cutoff point. The BCR-free survival rates, Differentially expressed genes (DEGs) and immune cells infiltrations in different TMB groups were identified.

RESULTS

The most common variant type and SNV were single nucleotide polymorphism and C > T. respectively. High TMB level was significantly associated with older age, positive lymph node, higher International Society of Urological Pathology (ISUP) grade, advanced stage and poor BCR-free survival. 132 DEGs were identified and involved in receptor ligand activity and hormone activity. High expression of six core genes UBE2C, PLK1, CDC20, BUB1, CDK1 and HJURP were associated with worse BCR-free survival. The analysis of immune cells infiltration revealed that the amount of activated CD4⁺ memory T cells was significantly different in high and low TMB groups.

CONCLUSIONS

The current study comprehensively described the summary of mutation and TMB related DEGs in PCa. TMB was associated with BCR-free survival and the infiltration of activated CD4⁺ memory T cells in the immune microenvironment.

Genomic analyses of early responses to radiation inglioblastoma reveal new alterations at transcription,splicing, and translation levels

High-dose radiation is the main component of glioblastoma therapy. Unfortunately, radio-resistance is a common problem and a major contributor to tumor relapse. Understanding the molecular mechanisms driving response to radiation is critical for identifying regulatory routes that could be targeted to improve treatment response. We conducted an integrated analysis in the U251 and U343 glioblastoma cell lines to map early alterations in the expression of genes at three levels: transcription, splicing, and translation in response to ionizing radiation. Changes at the transcriptional level were the most prevalent response. Downregulated genes are strongly associated with cell cycle and DNA replication and linked to a coordinated module of expression. Alterations in this group are likely driven by decreased expression of the transcription factor FOXM1 and members of the E2F family. Genes involved in RNA regulatory mechanisms were affected at the mRNA, splicing, and translation levels, highlighting their importance in radiation-response. We identified a number of oncogenic factors, with an increased expression upon radiation exposure, including BCL6, RRM2B, IDO1, FTH1, APIP, and LRIG2 and lncRNAs NEAT1 and FTX. Several of these targets have been previously implicated in radio-resistance. Therefore, antagonizing their effects post-radiation could increase therapeutic efficacy. Our integrated analysis provides a comprehensive view of early response to radiation in glioblastoma. We identify new biological processes involved in altered expression of various oncogenic factors and suggest new target options to increase radiation sensitivity and prevent relapse.