Epigenetic factor EPC1 is a master regulator of DNA damage response by interacting with E2F1 to silence death and activate metastasis-related gene signatures

EPC1/2 regulate hematopoietic stem and progenitor cell proliferation by modulating H3 acetylation and DLST

Enhancers of polycomb 1 (EPC1) and 2 (EPC2) are involved in multiple biological processes as components of histone acetyltransferases/deacetylase complexes and transcriptional cofactors, and their dysfunction was associated with developmental defects and diseases. However, it remains unknown how their dysfunction induces hematopoietic stem and progenitor cell (HSPC) defects. Here, we show that depletion of EPC1/2 significantly reduced the number of hematopoietic stem and progenitor cells (HSPCs) in the aorta-gonad mesonephros and caudal hematopoietic tissue regions by impairing HSPC proliferation, and consistently downregulated the expression of HSPC genes in K562 cells. This study demonstrates the functions of EPC1/2 in regulating histone H3 acetylation, and in regulating DLST (dihydrolipoamide S-succinyltransferase) via H3 acetylation and cooperating with transcription factors serum response factor and FOXR2 together, and in the subsequent HSPC emergence and proliferation. Our results demonstrate the essential roles of EPC1/2 in regulating H3 acetylation, and DLST as a linkage between EPC1 and EPC2 with mitochondria metabolism, in HSPC emergence and proliferation.

Logic-based modeling and drug repurposing for the prediction of novel therapeutic targets and combination regimens against E2F1-driven melanoma progression

Melanoma presents increasing prevalence and poor outcomes. Progression to aggressive stages is characterized by overexpression of the transcription factor E2F1 and activation of downstream prometastatic gene regulatory networks (GRNs). Appropriate therapeutic manipulation of the E2F1-governed GRNs holds the potential to prevent metastasis however, these networks entail complex feedback and feedforward regulatory motifs among various regulatory layers, which make it difficult to identify druggable components. To this end, computational approaches such as mathematical modeling and virtual screening are important tools to unveil the dynamics of these signaling networks and identify critical components that could be further explored as therapeutic targets. Herein, we integrated a well-established E2F1-mediated epithelial-mesenchymal transition (EMT) map with transcriptomics data from E2F1-expressing melanoma cells to reconstruct a core regulatory network underlying aggressive melanoma. Using logic-based in silico perturbation experiments of a core regulatory network, we identified that simultaneous perturbation of Protein kinase B (AKT1) and oncoprotein murine double minute 2 (MDM2) drastically reduces EMT in melanoma. Using the structures of the two protein signatures, virtual screening strategies were performed with the FDA-approved drug library. Furthermore, by combining drug repurposing and computer-aided drug design techniques, followed by molecular dynamics simulation analysis, we identified two potent drugs (Tadalafil and Finasteride) that can efficiently inhibit AKT1 and MDM2 proteins. We propose that these two drugs could be considered for the development of therapeutic strategies for the management of aggressive melanoma.

In silico analysis prediction of HepTH1-5 as a potential therapeutic agent by targeting tumour suppressor protein networks

Many studies reported that the activation of tumour suppressor protein, p53 induced the human hepcidin expression. However, its expression decreased when p53 was silenced in human hepatoma cells. Contrary to Tilapia hepcidin TH1-5, HepTH1-5 was previously reported to trigger the p53 activation through the molecular docking approach. The INhibitor of Growth (ING) family members are also shown to directly interact with p53 and promote cell cycle arrest, senescence, apoptosis and participate in DNA replication and DNA damage responses to suppress the tumour initiation and progression. However, the interrelation between INGs and HepTH1-5 remains unknown. Therefore, this study aims to identify the mechanism and their protein interactions using in silico approaches. The finding revealed that HepTH1-5 and its ligands had interacted mostly on hotspot residues of ING proteins which involved in histone modifications via acetylation, phosphorylation, and methylation. This proves that HepTH1-5 might implicate in an apoptosis signalling pathway and preserve the protein structure and function of INGs by reducing the perturbation of histone binding upon oxidative stress response. This study would provide theoretical guidance for the design and experimental studies to decipher the role of HepTH1-5 as a potential therapeutic agent for cancer therapy. Communicated by Ramaswamy H. Sarma.

Drug Repurposing at the Interface of Melanoma Immunotherapy and Autoimmune Disease

Cancer cells have a remarkable ability to evade recognition and destruction by the immune system. At the same time, cancer has been associated with chronic inflammation, while certain autoimmune diseases predispose to the development of neoplasia. Although cancer immunotherapy has revolutionized antitumor treatment, immune-related toxicities and adverse events detract from the clinical utility of even the most advanced drugs, especially in patients with both, metastatic cancer and pre-existing autoimmune diseases. Here, the combination of multi-omics, data-driven computational approaches with the application of network concepts enables in-depth analyses of the dynamic links between cancer, autoimmune diseases, and drugs. In this review, we focus on molecular and epigenetic metastasis-related processes within cancer cells and the immune microenvironment. With melanoma as a model, we uncover vulnerabilities for drug development to control cancer progression and immune responses. Thereby, drug repurposing allows taking advantage of existing safety profiles and established pharmacokinetic properties of approved agents. These procedures promise faster access and optimal management for cancer treatment. Together, these approaches provide new disease-based and data-driven opportunities for the prediction and application of targeted and clinically used drugs at the interface of immune-mediated diseases and cancer towards next-generation immunotherapies.

KAT2A/E2F1 Promotes Cell Proliferation and Migration via Upregulating the Expression of UBE2C in Pan-Cancer

Various studies have shown that lysine acetyltransferase 2A (KAT2A), E2F transcription factor 1 (E2F1), and ubiquitin conjugating enzyme E2 C (UBE2C) genes regulated the proliferation and migration of tumor cells through regulating the cell cycle. However, there is a lack of in-depth and systematic research on their mechanisms of action. This study analyzed The Cancer Genome Atlas (TCGA) to screen potential candidate genes and the regulation network of KAT2A and E2F1 complex in pan-cancer. Quantitative real-time PCR (qRT-PCR) and Western blotting (WB), cell phenotype detection, immunofluorescence co-localization, chromatin immunoprecipitation assay (ChIP), and RNA-Seq techniques were used to explore the functional of a candidate gene, UBE2C. We found that the expression of these three genes was significantly higher in more than 10 tumor types compared to normal tissue. Moreover, UBE2C was mainly expressed in tumor cells, which highlighted the impacts of UBE2C as a specific therapeutic strategy. Moreover, KAT2A and E2F1 could promote cell proliferation and the migration of cancer cells by enhancing the expression of UBE2C. Mechanically, KAT2A was found to cooperate with E2F1 and be recruited by E2F1 to the UBE2C promoter for elevating the expression of UBE2C by increasing the acetylation level of H3K9.

Loss of Rbl2 (Retinoblastoma-Like 2) Exacerbates Myocardial Ischemia/Reperfusion Injury

Background The postmitotic state of adult cardiomyocytes, maintained by the cell cycle repressor Rbl2 (retinoblastoma-like 2), is associated with considerable resistance to apoptosis. However, whether Rbl2 regulates cardiomyocyte apoptosis remains unknown. Methods and Results Here, we show that ablation of Rbl2 increased cardiomyocyte apoptosis following acute myocardial ischemia/reperfusion injury, leading to diminished cardiac function and exaggerated ventricular remodeling in the long term. Mechanistically, ischemia/reperfusion induced expression of the proapoptotic protein BCL2 interacting protein 3 (Bnip3), which was augmented by deletion of Rbl2. Because the Bnip3 promoter contains an adenoviral early region 2 binding factor (E2F)-binding site, we further showed that loss of Rbl2 upregulated the transcriptional activator E2F1 but downregulated the transcriptional repressor E2F4. In cultured cardiomyocytes, treatment with H₂O₂ markedly increased the levels of E2F1 and Bnip3, resulting in mitochondrial depolarization and apoptosis. Depletion of Rbl2 significantly augmented H₂O₂-induced mitochondrial damage and apoptosis in vitro. Conclusions Rbl2 deficiency enhanced E2F1-mediated Bnip3 expression, resulting in aggravated cardiomyocyte apoptosis and ischemia/reperfusion injury. Our results uncover a novel antiapoptotic role for Rbl2 in cardiomyocytes, suggesting that the cell cycle machinery may directly regulate apoptosis in postmitotic cardiomyocytes. These findings may be exploited to develop new strategies to limit ischemia/reperfusion injury in the treatment of acute myocardial infarction.

Multiomics Integrative Analysis Identifying EPC1 as a Prognostic Biomarker in Head and Neck Squamous Cell Carcinoma

Background. Biomarker research in head and neck squamous cell carcinoma (HNSCC) is constantly revealing promising findings. An enhancer of polycomb homolog 1 (EPC1) was found to play a procancer role in nasopharyngeal carcinoma (NPC), but its role in HNSCC with strong heterogeneity is still unclear. Herein, we investigated the prognostic significance and related mechanisms of EPC1 in HNSCC. Methods. The Kaplan-Meier plotter was used to evaluate the prognostic significance of EPC1. Based on a range of published public databases, the multiomics expression of EPC1 in HNSCC was explored to investigate the mechanisms affecting prognosis. Results. According to the clinical data, high EPC1 expression in HNSCC was a predictor of patient prognosis (hazard ratio $\left(\text{HR}\right)=0.64$ ; 95% confidence interval (CI) 0.49-0.83; $P<0.01$ ). EPC1 expression varied among clinical subtypes and was related to key factors, such as TP53 and human papillomavirus (HPV) ( $P<0.05$ ). At the genetic level, EPC1 expression level may be associated with protein phosphorylation, cell adhesion, cancer-related pathways, etc. For the noncoding region, a competing endogenous RNA network was constructed, and 6 microRNAs and 12 long noncoding RNAs were identified. At the protein level, a protein-protein interaction (PPI) network related to EPC1 expression was constructed and found to be involved in HPV infection, endocrine resistance, and multiple cancer pathways. At the immune level, EPC1 expression was correlated with a variety of immune cells and immune molecules, which together constituted the immune microenvironments of tumors. Conclusion. High EPC1 expression may predict a better prognosis in HNSCC, as it is more frequently found in HNSCC with HPV infection. EPC1 may participate in the genomics, transcriptomics, proteomics, and immunomics of HNSCC, and the results can provide a reference for the development of targeted drugs and evaluation of patient prognosis.

Comprehensive Pan-Cancer Analysis of Senescence With Cancer Prognosis and Immunotherapy

Senescence is a double-edged sword in tumorigenesis and affects the immunotherapy response through the modulation of the host’s immune system. However, there is currently a lack of comprehensive analysis of the senescence-related genes (SRGs) in human cancers, and the predictive role of senescence in cancer immunotherapy response has not been explored. The multi-omics approaches were performed in this article to conduct a systematic pan-cancer genomic analysis of SRGs in cancer. In addition, we calculated the generic senescence score (SS) to quantify the senescence levels in cancers and explored the correlations of SS with cancer prognosis, biological processes, and tumor microenvironment (TME). The gene signatures were deregulated in multiple cancers and indicated a context-dependent correlation with prognosis, tumor-immune evasion, and response to therapy across various tumor types. Further analysis disclosed that SS was positively associated with the infiltration levels of immune suppressive cells, including induced Tregs (iTregs), central memory Ts (Tcms), and natural Tregs (nTregs), and negatively associated with immune killer cells, including natural killers (NKs) and mucosal-associated invariant Ts (MAITs). Moreover, the SS was significantly correlated with tumor-associated macrophages (TAMs), cancer-associated fibroblasts (CAFs), immune-related genes, and immune checkpoints and had a predictive value of immunotherapy response. Thus, the expression of SRGs was involved in resistance to several anticancer drugs. Our work illustrates the characterization of senescence across various malignancies and highlights the potential of senescence as a biomarker of the response to immunotherapy.

Fusion genes in gynecologic tumors: the occurrence, molecular mechanism and prospect for therapy

Gene fusions are thought to be driver mutations in multiple cancers and are an important factor for poor patient prognosis. Most of them appear in specific cancers, thus satisfactory strategies can be developed for the precise treatment of these types of cancer. Currently, there are few targeted drugs to treat gynecologic tumors, and patients with gynecologic cancer often have a poor prognosis because of tumor progression or recurrence. With the application of massively parallel sequencing, a large number of fusion genes have been discovered in gynecologic tumors, and some fusions have been confirmed to be involved in the biological process of tumor progression. To this end, the present article reviews the current research status of all confirmed fusion genes in gynecologic tumors, including their rearrangement mechanism and frequency in ovarian cancer, endometrial cancer, endometrial stromal sarcoma, and other types of uterine tumors. We also describe the mechanisms by which fusion genes are generated and their oncogenic mechanism. Finally, we discuss the prospect of fusion genes as therapeutic targets in gynecologic tumors.

The miR-378c-Samd1 circuit promotes phenotypic modulation of vascular smooth muscle cells and foam cells formation in atherosclerosis lesions

MicroRNAs have emerged as key regulators in vascular diseases and are involved in the formation of atherosclerotic lesions. However, the atherosclerotic-specific MicroRNAs and their functional roles in atherosclerosis are unclear. Here, we report that miR-378c protects against atherosclerosis by directly targeting Sterile Alpha Motif Domain Containing 1 (Samd1), a predicted transcriptional repressor. miR-378c was strikingly reduced in atherosclerotic plaques and blood of acute coronary syndrome (ACS) patients relative to healthy controls. Suppression of miR-378c promoted vascular smooth muscle cells (VSMCs) phenotypic transition during atherosclerosis. We also reported for the first time that Samd1 prolonged immobilization of LDL on the VSMCs, thus facilitated LDL oxidation and subsequently foam cell formation. Further, we found that Samd1 contains predicted DNA binding domain and directly binds to DNA regions as a transcriptional repressor. Together, we uncovered a novel mechanism whereby miR-378c-Samd1 circuit participates in two key elements of atherosclerosis, VSMCs phenotypic transition and LDL oxidation. Our results provided a better understanding of atherosclerosis pathophysiology and potential therapeutic management by targeting miR-378c-Samd1 circuit.

Transcriptional E2F1/2/3/6 as potential prognostic biomarkers in cutaneous melanoma

Although the abnormal expression of members of the E2F family has been reported to participate in carcinogenesis in many human types of cancer, the bioinformatics role of the E2F family in melanoma is unknown. This research was designed to detect the expression, methylation, prognostic value and potential effects of the E2F family in melanoma. We investigated E2F family mRNA expression from the Oncomine and GEPIA databases and their methylation status in the MethHC database. Meanwhile, we detected the relative E2F family expression levels by qPCR and immunohistochemistry. Kaplan-Meier Plotter was used to draw survival analysis charts, and gene functional enrichment analyses were applied through cBioPortal database analysis. E2F1/2/3/4/5/6 mRNA and proteins were clearly upregulated in cutaneous melanoma patients, and high expression levels of E2F1/2/3/6 were statistically related to high methylation levels. Increased mRNA expression of E2F1/2/3/6 was related to lower overall survival rates (OS) and disease-free survival (DFS) in cutaneous melanoma cases. Meanwhile, E2F1/2/3/6 carried out these effects through regulating multiple signaling pathways, including the MAPK, PI3K-Akt and p53 signaling pathways. Taking together, our findings suggest that E2F1/2/3/6 could act as potential targets for precision therapy in cutaneous melanoma patients.

Development and validation of a nomogram with an epigenetic signature for predicting survival in patients with lung adenocarcinoma

Epigenetic factors play crucial roles in carcinogenesis by modifying chromatin architecture. Here, we established an epigenetic biosignature-based model for examining survival in patients with lung adenocarcinoma (LUAD). We retrieved gene-expression profiles and clinical data from The Cancer Genome Atlas and Gene Expression Omnibus and clustered the data into training (n = 490) and Validation (n = 226) datasets, respectively. To establish an epigenetic model, we identified prognostic epigenetic regulation-related genes by LASSO and Cox regression analyses, and established a novel 11-gene signature, including EPC1, GADD45A, HCFC2, RCOR1, SMARCAL1, TLE2, TRIM28, and ZNF516, for predicting LUAD overall survival (OS). The biosignature performed optimally in both the training and validation sets according to receiver operating characteristic and calibration plots. Moreover, the biosignature classified patients into high- and low-risk clusters with distinct survival times, with Cox regression analysis revealing the biosignature as an independent LUAD prognostic index. Furthermore, the generated nomogram integrating the prognostic gene biosignature and clinical indices predicted LUAD OS with high efficiency and outperformed tumor-node-metastasis staging in LUAD survival prediction. These results demonstrated the efficacy of the epigenetic signature prognostic nomogram for reliably predicting LUAD OS and its potential application for informing clinical decision making and individualized treatment.

Modeling epistasis in mice and yeast using the proportion of two or more distinct genetic backgrounds: Evidence for "polygenic epistasis"

Background

The majority of quantitative genetic models used to map complex traits assume that alleles have similar effects across all individuals. Significant evidence suggests, however, that epistatic interactions modulate the impact of many alleles. Nevertheless, identifying epistatic interactions remains computationally and statistically challenging. In this work, we address some of these challenges by developing a statistical test for polygenic epistasis that determines whether the effect of an allele is altered by the global genetic ancestry proportion from distinct progenitors.

Results

We applied our method to data from mice and yeast. For the mice, we observed 49 significant genotype-by-ancestry interaction associations across 14 phenotypes as well as over 1,400 Bonferroni-corrected genotype-by-ancestry interaction associations for mouse gene expression data. For the yeast, we observed 92 significant genotype-by-ancestry interactions across 38 phenotypes. Given this evidence of epistasis, we test for and observe evidence of rapid selection pressure on ancestry specific polymorphisms within one of the cohorts, consistent with epistatic selection.

Conclusions

Unlike our prior work in human populations, we observe widespread evidence of ancestry-modified SNP effects, perhaps reflecting the greater divergence present in crosses using mice and yeast.

Many statistical tests which link genetic markers in the genome to differences in traits rely on the assumption that the same polymorphism will have identical effects in different individuals. However, there is substantial evidence indicating that this is not the case. Epistasis is the phenomenon in which multiple polymorphisms interact with one another to amplify or negate each other’s effects on a trait. We hypothesized that individual SNP effects could be changed in a polygenic manner, such that the proportion of as genetic ancestry, rather than specific markers, might be used to capture epistatic interactions. Motivated by this possibility, we develop a new statistical test that allowed us to examine the genome to identify polymorphisms which have different effects depending on the ancestral makeup of each individual. We use our test in two different populations of inbred mice and a yeast panel and demonstrate that these sorts of variable effect polymorphisms exist in 14 different physical traits in mice and 38 phenotypes in yeast as well as in murine gene expression. We use the term “polygenic epistasis” to distinguish these interactions from the more conventional two- or multi-locus interactions.

Identification of Novel Fusion Genes in Bone and Soft Tissue Sarcoma and Their Implication in the Generation of a Mouse Model

Fusion genes induced by chromosomal aberrations are common mutations causally associated with bone and soft tissue sarcomas (BSTS). These fusions are usually disease type-specific, and identification of the fusion genes greatly helps in making precise diagnoses and determining therapeutic directions. However, there are limitations in detecting unknown fusion genes or rare fusion variants when using standard fusion gene detection techniques, such as reverse transcription-polymerase chain reaction (RT-PCR) and fluorescence in situ hybridization (FISH). In the present study, we have identified 19 novel fusion genes using target RNA sequencing (RNA-seq) in 55 cases of round or spindle cell sarcomas in which no fusion genes were detected by RT-PCR. Subsequent analysis using Sanger sequencing confirmed that seven out of 19 novel fusion genes would produce functional fusion proteins. Seven fusion genes detected in this study affect signal transduction and are ideal targets of small molecule inhibitors. YWHAE-NTRK3 expression in mouse embryonic mesenchymal cells (eMCs) induced spindle cell sarcoma, and the tumor was sensitive to the TRK inhibitor LOXO-101 both in vitro and in vivo. The combination of target RNA-seq and generation of an ex vivo mouse model expressing novel fusions provides important information both for sarcoma biology and the appropriate diagnosis of BSTS.

LncRNA-SLC16A1-AS1 induces metabolic reprogramming during Bladder Cancer progression as target and co-activator of E2F1

Long non-coding RNAs (lncRNAs) have emerged as integral components of E2F1-regulated gene regulatory networks (GRNs), but their implication in advanced or treatment-refractory malignancy is unknown. Methods: We combined high-throughput transcriptomic approaches with bioinformatics and structure modeling to search for lncRNAs that participate in E2F1-activated prometastatic GRNs and their phenotypic targets in the highly-relevant case of E2F1-driven aggressive bladder cancer (BC). RNA immunoprecipitation was performed to verify RNA-protein interactions. Functional analyses including qRT-PCR, immunoblotting, luciferase assays and measurement of extracellular fluxes were conducted to validate expression and target gene regulation. Results: We identified E2F1-responsive lncRNA-SLC16A1-AS1 and its associated neighboring protein-coding gene, SLC16A1/MCT1, which both promote cancer invasiveness. Mechanistically, upon E2F1-mediated co-transactivation of the gene pair, SLC16A1-AS1 associates with E2F1 in a structure-dependent manner and forms an RNA-protein complex that enhances SLC16A1/MCT1 expression through binding to a composite SLC16A1-AS1:E2F1-responsive promoter element. Moreover, SLC16A1-AS1 increases aerobic glycolysis and mitochondrial respiration and fuels ATP production by fatty acid β-oxidation. These metabolic changes are accompanied by alterations in the expression of the SLC16A1-AS1:E2F1-responsive gene PPARA, a key mediator of fatty acid β-oxidation. Conclusions: Our results unveil a new gene regulatory program by which E2F1-induced lncRNA-SLC16A1-AS1 forms a complex with its transcription factor that promotes cancer metabolic reprogramming towards the acquisition of a hybrid oxidative phosphorylation/glycolysis cell phenotype favoring BC invasiveness.

RNA isoform screens uncover the essentiality and tumor-suppressor activity of ultraconserved poison exons

While RNA-seq has enabled comprehensive quantification of alternative splicing, no correspondingly high-throughput assay exists for functionally interrogating individual isoforms. We describe pgFARM (paired guide RNAs for alternative exon removal), a CRISPR-Cas9-based method to manipulate isoforms independent of gene inactivation. This approach enabled rapid suppression of exon recognition in polyclonal settings to identify functional roles for individual exons, such as an SMNDC1 cassette exon that regulates pan-cancer intron retention. We generalized this method to a pooled screen to measure the functional relevance of 'poison' cassette exons, which disrupt their host genes' reading frames yet are frequently ultraconserved. Many poison exons were essential for the growth of both cultured cells and lung adenocarcinoma xenografts, while a subset had clinically relevant tumor-suppressor activity. The essentiality and cancer relevance of poison exons are likely to contribute to their unusually high conservation and contrast with the dispensability of other ultraconserved elements for viability.

Topoisomerase inhibitors promote cancer cell motility via ROS-mediated activation of JAK2-STAT1-CXCL1 pathway

Background

Topoisomerase inhibitors (TI) can inhibit cell proliferation by preventing DNA replication, stimulating DNA damage and inducing cell cycle arrest. Although these agents have been commonly used in the chemotherapy for the anti-proliferative effect, their impacts on the metastasis of cancer cells remain obscure.

Methods

We used the transwell chamber assay to test effects of Topoisomerase inhibitors Etoposide (VP-16), Adriamycin (ADM) and Irinotecan (CPT-11) on the migration and invasion of cancer cells. Conditioned medium (CM) from TI-treated cells was subjected to Mass spectrometry screening. Gene silencing, neutralizing antibody, and specific chemical inhibitors were used to validate the roles of signaling molecules.

Results

Our studies disclosed that TI could promote the migration and invasion of a subset of cancer cells, which were dependent on chemokine (C-X-C motif) ligand 1 (CXCL1). Further studies disclosed that TI enhanced phosphorylation of Janus kinase 2 (JAK2) and Signal transducers and activators of transcription 1 (STAT1). Silencing or chemical inhibition of JAK2 or STAT1 abrogated TI-induced CXCL1 expression and cell motility. Moreover, TI increased cellular levels of reactive oxygen species (ROS) and promoted oxidation of Protein Tyrosine Phosphatase 1B (PTP1B), while reduced glutathione (GSH) reversed TI-induced JAK2-STAT1 activation, CXCL1 expression, and cell motility.

Conclusions

Our study demonstrates that TI can promote the expression and secretion of CXCL1 by elevating ROS, inactivating PTP1B, and activating JAK2-STAT1 signaling pathway, thereby promoting the motility of cancer cells.

Electronic supplementary material

The online version of this article (10.1186/s13046-019-1353-2) contains supplementary material, which is available to authorized users.

Drug Repositioning Inferred from E2F1-Coregulator Interactions Studies for the Prevention and Treatment of Metastatic Cancers

Metastasis management remains a long-standing challenge. High abundance of E2F1 triggers tumor progression by developing protein-protein interactions (PPI) with coregulators that enhance its potential to activate a network of prometastatic transcriptional targets.

Methods: To identify E2F1-coregulators, we integrated high-throughput Co-immunoprecipitation (IP)/mass spectometry, GST-pull-down assays, and structure modeling. Potential inhibitors of PPI discovered were found by bioinformatics-based pharmacophore modeling, and transcriptome profiling was conducted to screen for coregulated downstream targets. Expression and target gene regulation was validated using qRT-PCR, immunoblotting, chromatin IP, and luciferase assays. Finally, the impact of the E2F1-coregulator complex and its inhibiting drug on metastasis was investigated in vitro in different cancer entities and two mouse metastasis models.

Results: We unveiled that E2F1 forms coactivator complexes with metastasis-associated protein 1 (MTA1) which, in turn, is directly upregulated by E2F1. The E2F1:MTA1 complex potentiates hyaluronan synthase 2 (HAS2) expression, increases hyaluronan production and promotes cell motility. Disruption of this prometastatic E2F1:MTA1 interaction reduces hyaluronan synthesis and infiltration of tumor-associated macrophages in the tumor microenvironment, thereby suppressing metastasis. We further demonstrate that E2F1:MTA1 assembly is abrogated by small-molecule, FDA-approved drugs. Treatment of E2F1/MTA1-positive, highly aggressive, circulating melanoma cells and orthotopic pancreatic tumors with argatroban prevents metastasis and cancer relapses in vivo through perturbation of the E2F1:MTA1/HAS2 axis.

Conclusion: Our results propose argatroban as an innovative, E2F-coregulator-based, antimetastatic drug. Cancer patients with the infaust E2F1/MTA1/HAS2 signature will likely benefit from drug repositioning.

A Well-Controlled BioID Design for Endogenous Bait Proteins

The CRISPR/Cas9 revolution is profoundly changing the way life sciences technologies are used. Many assays now rely on engineered clonal cell lines to eliminate the overexpression of bait proteins. Control cell lines are typically nonengineered cells or engineered clones, implying a considerable risk for artifacts because of clonal variation. Genome engineering can also transform BioID, a proximity labeling method that relies on fusing a bait protein to a promiscuous biotin ligase, BirA*, resulting in the tagging of vicinal proteins. We here propose an innovative design to enable BioID for endogenous proteins wherein we introduce a T2A-BirA* module at the C-terminus of endogenous p53 by genome engineering, leading to bicistronic expression of both p53 and BirA* under control of the endogenous promoter. By targeting a Cas9-cytidine deaminase base editor to the T2A autocleavage site, we can efficiently derive an isogenic population expressing a functional p53-BirA* fusion protein. Using quantitative proteomics we show significant benefits over the classical ectopic expression of p53-BirA*, and we provide a first well-controlled view of the proximal proteins of endogenous p53 in colon carcinoma cells. This novel application for base editors expands the CRISPR/Cas9 toolbox and can be a valuable addition for synthetic biology.

Emerging functional markers for cancer stem cell-based therapies: Understanding signaling networks for targeting metastasis

Metastasis is one of the most challenging issues in cancer patient management, and effective therapies to specifically target disease progression are missing, emphasizing the urgent need for developing novel anti-metastatic therapeutics. Cancer stem cells (CSCs) gained fast attention as a minor population of highly malignant cells within liquid and solid tumors that are responsible for tumor onset, self-renewal, resistance to radio- and chemotherapies, and evasion of immune surveillance accelerating recurrence and metastasis. Recent progress in the identification of their phenotypic and molecular characteristics and interactions with the tumor microenvironment provides great potential for the development of CSC-based targeted therapies and radical improvement in metastasis prevention and cancer patient prognosis. Here, we report on newly uncovered signaling mechanisms controlling CSC's aggressiveness and treatment resistance, and CSC-specific agents and molecular therapeutics, some of which are currently under investigation in clinical trials, gearing towards decisive functional CSC intrinsic or surface markers. One special research focus rests upon subverted regulatory pathways such as insulin-like growth factor 1 receptor signaling and its interactors in metastasis-initiating cell populations directly related to the gain of stem cell- and EMT-associated properties, as well as key components of the E2F transcription factor network regulating metastatic progression, microenvironmental changes, and chemoresistance. In addition, the study provides insight into systems biology tools to establish complex molecular relationships behind the emergence of aggressive phenotypes from high-throughput data that rely on network-based analysis and their use to investigate immune escape mechanisms or predict clinical outcome-relevant CSC receptor signaling signatures. We further propose that customized vector technologies could drastically enhance systemic drug delivery to target sites, and summarize recent progress and remaining challenges. This review integrates available knowledge on CSC biology, computational modeling approaches, molecular targeting strategies, and delivery techniques to envision future clinical therapies designed to conquer metastasis-initiating cells.

RNAi-Based Identification of Gene-Specific Nuclear Cofactor Networks Regulating Interleukin-1 Target Genes

The potent proinflammatory cytokine interleukin (IL)-1 triggers gene expression through the NF-κB signaling pathway. Here, we investigated the cofactor requirements of strongly regulated IL-1 target genes whose expression is impaired in p65 NF-κB-deficient murine embryonic fibroblasts. By two independent small-hairpin (sh)RNA screens, we examined 170 genes annotated to encode nuclear cofactors for their role in Cxcl2 mRNA expression and identified 22 factors that modulated basal or IL-1-inducible Cxcl2 levels. The functions of 16 of these factors were validated for Cxcl2 and further analyzed for their role in regulation of 10 additional IL-1 target genes by RT-qPCR. These data reveal that each inducible gene has its own (quantitative) requirement of cofactors to maintain basal levels and to respond to IL-1. Twelve factors (Epc1, H2afz, Kdm2b, Kdm6a, Mbd3, Mta2, Phf21a, Ruvbl1, Sin3b, Suv420h1, Taf1, and Ube3a) have not been previously implicated in inflammatory cytokine functions. Bioinformatics analysis indicates that they are components of complex nuclear protein networks that regulate chromatin functions and gene transcription. Collectively, these data suggest that downstream from the essential NF-κB signal each cytokine-inducible target gene has further subtle requirements for individual sets of nuclear cofactors that shape its transcriptional activation profile.

Blockage of the mevalonate pathway overcomes the apoptotic resistance to MEK inhibitors with suppressing the activation of Akt in cancer cells

With increasing clinical demands for MEK inhibitors in cancer treatment, overcoming the resistance to MEK inhibitors is an urgent problem to be solved. Numerous reports have shown that MEK inhibition results in the activation of PI3K-Akt signaling, which may confer apoptotic resistance to MEK inhibitors. We here demonstrate that the blockade of the mevalonate pathway using the antilipidemic drug statins represses Akt activation following MEK inhibition and induces significant apoptosis when co-treated with CH5126766 or trametinib. These events were clearly negated by the addition of mevalonate or geranylgeranyl pyrophosphate, indicating that the protein geranylgeranylation is implicated in the apoptotic resistance to MEK inhibitors. Furthermore, mechanistically, the combined treatment of CH5126766 with statins upregulated TNF-related apoptosis-inducing ligand (TRAIL), which was dependent on inhibition of the mevalonate pathway and is involved in apoptosis induction in human breast cancer MDA-MB-231 cells. The present study not only revealed that the mevalonate pathway could be targetable to enhance the efficacy of MEK inhibitors, but also proposes that combinatorial treatment of MEK inhibitors with statins may be a promising therapeutic strategy to sensitize cancer cells to apoptosis.

A web platform for the network analysis of high-throughput data in melanoma and its use to investigate mechanisms of resistance to anti-PD1 immunotherapy

Cellular phenotypes are established and controlled by complex and precisely orchestrated molecular networks. In cancer, mutations and dysregulations of multiple molecular factors perturb the regulation of these networks and lead to malignant transformation. High-throughput technologies are a valuable source of information to establish the complex molecular relationships behind the emergence of malignancy, but full exploitation of this massive amount of data requires bioinformatics tools that rely on network-based analyses. In this report we present the Virtual Melanoma Cell, an online tool developed to facilitate the mining and interpretation of high-throughput data on melanoma by biomedical researches. The platform is based on a comprehensive, manually generated and expert-validated regulatory map composed of signaling pathways important in malignant melanoma. The Virtual Melanoma Cell is a tool designed to accept, visualize and analyze user-generated datasets. It is available at: https://www.vcells.net/melanoma. To illustrate the utilization of the web platform and the regulatory map, we have analyzed a large publicly available dataset accounting for anti-PD1 immunotherapy treatment of malignant melanoma patients.

Advances in cancer stem cell targeting: How to strike the evil at its root

Cancer progression to metastatic stages is still unmanageable and the promise of effective anti-metastatic therapy remains largely unmet, emphasizing the need to develop novel therapeutics. The special focus here is on cancer stem cells (CSC) as the seed of tumor initiation, epithelial-mesenchymal transition, chemoresistance and, as a consequence, drivers of metastatic dissemination. We report on targeted therapies gearing towards the CSC's internal and membrane-anchored markers using agents such as antibody derivatives, nucleic therapeutics, small molecules and genetic payloads. Another emphasis lies on novel proceedings envisaged to deliver current and prospective therapies to the target sites using newest viral and non-viral vector technologies. In this review, we summarize recent progress and remaining challenges in therapeutic strategies to combat CSC.

Critical genomic regulation mediated by Enhancer of Polycomb

Enhancer of Polycomb (EPC) was first identified for its contributions to development in Drosophila and was soon-thereafter purified as a subunit of the NuA4/TIP60 acetyltransferase complex. Since then, EPC has often been left in the shadows as an essential, yet non-catalytic subunit of NuA4/TIP60; however, its deep conservation and disease association make clear that it warrants additional attention. In fact, recent studies in yeast demonstrated that its Enhancer of Polycomb, Epl1, was just as important for gene expression and acetylation as is the catalytic subunit of NuA4. Despite its conservation, studies of EPC have often remained siloed between organisms. Here, our goal is to provide a cohesive view of the current state of the EPC literature as it stands among the major model organisms in which it has been studied. EPC is involved in multiple processes, beginning with its cardinal role in regulating global and targeted histone acetylation. EPC also frequently serves as an important interaction partner in these basic cellular functions, as well as in multicellular development, such as in hematopoiesis and skeletal muscle differentiation, and in human disease. Taken together, a unifying theme from these studies highlights EPC as a critical genomic regulator.

Expression level is a key determinant of E2F1-mediated cell fate

The Rb/E2F network has a critical role in regulating cell cycle progression and cell fate decisions. It is dysfunctional in virtually all human cancers, because of genetic lesions that cause overexpression of activators, inactivation of repressors, or both. Paradoxically, the downstream target of this network, E2F1, is rarely strongly overexpressed in cancer. E2F1 can induce both proliferation and apoptosis but the factors governing these critical cell fate decisions remain unclear. Previous studies have focused on qualitative mechanisms such as differential cofactors, posttranslational modification or state of other signaling pathways as modifiers of the cell fate decisions downstream of E2F1 activation. In contrast, the importance of the expression levels of E2F1 itself in dictating the downstream phenotypes has not been rigorously studied, partly due to the limited resolution of traditional population-level measurements. Here, through single-cell quantitative analysis, we demonstrate that E2F1 expression levels have a critical role in determining the fate of individual cells. Low levels of exogenous E2F1 promote proliferation, moderate levels induce G1, G2 and mitotic cell cycle arrest, and very high levels promote apoptosis. These multiple anti-proliferative mechanisms result in a strong selection pressure leading to rapid elimination of E2F1-overexpressing cells from the population. RNA-sequencing and RT-PCR revealed that low levels of E2F1 are sufficient to induce numerous cell cycle-promoting genes, intermediate levels induce growth arrest genes (i.e., p18, p19 and p27), whereas higher levels are necessary to induce key apoptotic E2F1 targets APAF1, PUMA, HRK and BIM. Finally, treatment of a lung cancer cell line with a proteasome inhibitor, MLN2238, resulted in an E2F1-dependent mitotic arrest and apoptosis, confirming the role of endogenous E2F1 levels in these phenotypes. The strong anti-proliferative activity of moderately overexpressed E2F1 in multiple cancer types suggests that targeting E2F1 for upregulation may represent an attractive therapeutic strategy in cancer.

Chromatin Regulation by the NuA4 Acetyltransferase Complex Is Mediated by Essential Interactions Between Enhancer of Polycomb (Epl1) and Esa1

Enzymes that modify and remodel chromatin act in broadly conserved macromolecular complexes. One key modification is the dynamic acetylation of histones and other chromatin proteins by opposing activities of acetyltransferase and deacetylase complexes. Among acetyltransferases, the NuA4 complex containing Tip60 or its Saccharomyces cerevisiae ortholog Esa1 is of particular significance because of its roles in crucial genomic processes including DNA damage repair and transcription. The catalytic subunit Esa1 is essential, as are five noncatalytic NuA4 subunits. We found that of the noncatalytic subunits, deletion of Enhancer of polycomb (Epl1), but not the others, can be bypassed by loss of a major deacetylase complex, a property shared by Esa1 Noncatalytic complex subunits can be critical for complex assembly, stability, genomic targeting, substrate specificity, and regulation. Understanding the essential role of Epl1 has been previously limited, a limitation now overcome by the discovery of its bypass suppression. Here, we present a comprehensive in vivo study of Epl1 using the powerful tool of suppression combined with transcriptional and mutational analyses. Our results highlight functional parallels between Epl1 and Esa1 and further illustrate that the structural role of Epl1 is important for promotion of Esa1 activity. This conclusion is strengthened by our dissection of Epl1 domains required in vivo for interaction with specific NuA4 subunits, histone acetylation, and chromatin targeting. These results provide new insights for the conserved, essential nature of Epl1 and its homologs, such as EPC1/2 in humans, which is frequently altered in cancers.

Deciphering hallmark processes of aging from interaction networks

BACKGROUND

Aging is broadly considered to be a dynamic process that accumulates unfavourable structural and functional changes in a time dependent fashion, leading to a progressive loss of physiological integrity of an organism, which eventually leads to age-related diseases and finally to death.

SCOPE OF REVIEW

The majority of aging-related studies are based on reductionist approaches, focusing on single genes/proteins or on individual pathways without considering possible interactions between them. Over the last few decades, several such genes/proteins were independently analysed and linked to a role that is affecting the longevity of an organism. However, an isolated analysis on genes and proteins largely fails to explain the mechanistic insight of a complex phenotype due to the involvement and integration of multiple factors.

MAJOR CONCLUSIONS

Technological advance makes it possible to generate high-throughput temporal and spatial data that provide an opportunity to use computer-based methods. These techniques allow us to go beyond reductionist approaches to analyse large-scale networks that provide deeper understanding of the processes that drive aging.

GENERAL SIGNIFICANCE

In this review, we focus on systems biology approaches, based on network inference methods to understand the dynamics of hallmark processes leading to aging phenotypes. We also describe computational methods for the interpretation and identification of important molecular hubs involved in the mechanistic linkage between aging related processes. This article is part of a Special Issue entitled "System Genetics" Guest Editor: Dr. Yudong Cai and Dr. Tao Huang.