Multi-Omics Integrated Analysis of the Protective Effect of EZH2 Inhibition in Mice with Renal Ischemia-Reperfusion Injury.

INTRODUCTION: Acute kidney injury (AKI) is a common clinical syndrome associated with high morbidity and mortality. Inhibition of the methyltransferase enhancer of zeste homolog 2 (EZH2) by its inhibitor 3-deazaneplanocin A (3-DZNeP) exerts renal benefits in acute renal ischemia-reperfusion injury (IRI). However, the underlying mechanisms are not completely known. This study aimed to elucidate the pathological mechanism of EZH2 in renal IRI by combination of multi-omics analysis and expression profiling in a public clinical cohort. METHODS: In this study, C57BL/6 J mice were used to establish the AKI model, which were treated with 3-DZNeP for 24 h. Kidney samples were collected for RNA-seq analysis, which was combined with publicly available EZH2 chromatin immunoprecipitation sequencing (ChIP-seq) data of mouse embryonic stem cell for a joint analysis to identify differentially expressed genes. Several selected differentially expressed genes were verified by quantitative PCR. Finally, single-nucleus sequencing data and expression profiling in public clinical datasets were used to confirm the negative correlation of the selected genes with EZH2 expression. RESULTS: 3-DZNeP treatment significantly improved renal pathology and function in IRI mice. Through RNA-seq analysis combined with EZH2 ChIP-seq database, 162 differentially expressed genes were found, which might be involved in EZH2-mediated pathology in IRI kidneys. Four differential expressed genes (Scd1, Cidea, Ghr, and Kl) related to lipid metabolism or cell growth were selected based on Gene Ontology and Kyoto Encyclopedia of Genes and Genome enrichment analysis, which were validated by quantitative PCR. Data from single-nucleus RNA sequencing revealed the negative correlation of these four genes with Ezh2 expression in different subpopulations of proximal tubular cells in IRI mice in a different pattern. Finally, the negative correlation of these four genes with EZH2 expression was confirmed in patients with AKI in two clinical datasets. CONCLUSIONS: Our study indicates that Scd1, Cidea, Ghr, and Kl are downstream genes regulated by EZH2 in AKI. Upregulation of EZH2 in AKI inhibits the expression of these four genes in a different population of proximal tubular cells to minimize normal physiological function and promote acute or chronic cell injuries following AKI.

A biomimetic nanoplatform for customized photothermal therapy of HNSCC evaluated on patient-derived xenograft models.

Cancer cell membrane (CCM) derived nanotechnology functionalizes nanoparticles (NPs) to recognize homologous cells, exhibiting translational potential in accurate tumor therapy. However, these nanoplatforms are majorly generated from fixed cell lines and are typically evaluated in cell line-derived subcutaneous-xenografts (CDX), ignoring the tumor heterogeneity and differentiation from inter- and intra- individuals and microenvironments between heterotopic- and orthotopic-tumors, limiting the therapeutic efficiency of such nanoplatforms. Herein, various biomimetic nanoplatforms (CCM-modified gold@Carbon, i.e., Au@C-CCM) were fabricated by coating CCMs of head and neck squamous cell carcinoma (HNSCC) cell lines and patient-derived cells on the surface of Au@C NP. The generated Au@C-CCMs were evaluated on corresponding CDX, tongue orthotopic xenograft (TOX), immune-competent primary and distant tumor models, and patient-derived xenograft (PDX) models. The Au@C-CCM generates a photothermal conversion efficiency up to 44.2% for primary HNSCC therapy and induced immunotherapy to inhibit metastasis via photothermal therapy-induced immunogenic cell death. The homologous CCM endowed the nanoplatforms with optimal targeting properties for the highest therapeutic efficiency, far above those with mismatched CCMs, resulting in distinct tumor ablation and tumor growth inhibition in all four models. This work reinforces the feasibility of biomimetic NPs combining modular designed CMs and functional cores for customized treatment of HNSCC, can be further extended to other malignant tumors therapy.

Multi-omics characterization of WNT pathway reactivation to ameliorate BET inhibitor resistance in liver cancer cells.

The Bromodomain and Extra-terminal domain (BET) proteins are promising targets in treating cancers. Although BET inhibitors have been in clinical trials, they are limited by lacking of suitable biomarkers to indicate drug responses in different cancers. Here we identify DHRS2, ETV4 and NOTUM as potential biomarkers to indicate drug resistance in liver cancer cells of a recently discovered BET inhibitor, Hjp-6-171. Furthermore, we confirm that reactivation of WNT pathway, the target of NOTUM, contributes to the drug sensitivity restoration in Hjp-6-171 resistant cells. Specially, combinations of Hjp-6-171 and a GSK3ß inhibitor CHIR-98014 show remarkable therapeutic effects in vitro and in vivo. Integrating RNA-seq and ChIP-seq data, we reveal the expression signature of ß-catenin regulated genes is contrary in sensitive cells to that in resistant cells. We propose WNT signaling molecules such as ß-catenin and ETV4 to be candidate biomarkers to indicate BET inhibitor responses in liver cancer patients.

Characterization of the genomic landscape and actionable mutations in Chinese breast cancers by clinical sequencing.

The remarkable advances in next-generation sequencing technology have enabled the wide usage of sequencing as a clinical tool. To promote the advance of precision oncology for breast cancer in China, here we report a large-scale prospective clinical sequencing program using the Fudan-BC panel, and comprehensively analyze the clinical and genomic characteristics of Chinese breast cancer. The mutational landscape of 1,134 breast cancers reveals that the most significant differences between Chinese and Western patients occurred in the hormone receptor positive, human epidermal growth factor receptor 2 negative breast cancer subtype. Mutations in p53 and Hippo signaling pathways are more prevalent, and 2 mutually exclusive and 9 co-occurring patterns exist among 9 oncogenic pathways in our cohort. Further preclinical investigation partially suggests that NF2 loss-of-function mutations can be sensitive to a Hippo-targeted strategy. We establish a public database (Fudan Portal) and a precision medicine knowledge base for data exchange and interpretation. Collectively, our study presents a leading approach to Chinese precision oncology treatment and reveals potentially actionable mutations in breast cancer.

Truncated HDAC9 identified by integrated genome-wide screen as the key modulator for paclitaxel resistance in triple-negative breast cancer.

Rationale: Paclitaxel resistance is a major concern when treating triple-negative breast cancer (TNBC) patients. We aimed to identify candidates causing paclitaxel resistance and explore their significance in TNBC therapeutics. Methods: A genome-wide CRISPR screening, integrated with transcriptome analyses, was performed to identify candidates involved in paclitaxel-resistant TNBCs. Cell proliferation, cytotoxicity, immunofluorescent staining, and xenograft assays were conducted to verify the phenotypes of paclitaxel resistance induced by candidate genes, both in vitro and in vivo. RNA sequencing, Western blotting, and chromatin immunoprecipitation assays were used to explore the underlying mechanisms. Results: MEF2-interacting transcriptional repressor (MITR), the truncated isoform of histone deacetylase 9 (HDAC9) lacking the deacetylation domain, was enriched in paclitaxel-resistant cells. Elevated MITR expression resulted in increased interleukin-11 (IL11) expression and activation of downstream JAK/STAT3 signaling. Mechanistically, MITR counteracted MEF2A-induced transcriptional suppression of IL11, ultimately causing paclitaxel resistance. By contrast, pharmacological inhibition of JAK1/2 by ruxolitinib reversed paclitaxel resistance both in vitro and in vivo. Conclusion: Our in vitro and in vivo genetic and cellular analyses elucidated the pivotal role of MITR/MEF2A/IL11 axis in paclitaxel resistance and provided a novel therapeutic strategy for TNBC patients to overcome poor chemotherapy responses.

Dissecting the heterogeneity of the alternative polyadenylation profiles in triple-negative breast cancers.

Background: Triple-negative breast cancer (TNBC) is an aggressive malignancy with high heterogeneity. However, the alternative polyadenylation (APA) profiles of TNBC remain unknown. Here, we aimed to define the characteristics of the APA events at post-transcription level among TNBCs. Methods: Using transcriptome microarray data, we analyzed APA profiles of 165 TNBC samples and 33 paired normal tissues. A pooled short hairpin RNA screen targeting 23 core cleavage and polyadenylation (C/P) genes was used to identify key C/P factors. Results: We established an unconventional APA subtyping system composed of four stable subtypes: 1) luminal androgen receptor (LAR), 2) mesenchymal-like immune-activated (MLIA), 3) basal-like (BL), 4) suppressed (S) subtypes. Patients in the S subtype had the worst disease-free survival comparing to other patients (log-rank p = 0.021). Enriched clinically actionable pathways and putative therapeutic APA events were analyzed among each APA subtype. Furthermore, CPSF1 and PABPN1 were identified as the master C/P factors in regulating APA events and TNBC proliferation. The depletion of CPSF1 or PABPN1 weakened cell proliferation, enhanced apoptosis, resulted in cell cycle redistribution and a reversion of APA events of genes associated with tumorigenesis, proliferation, metastasis and chemosensitivity in breast cancer. Conclusions: Our findings advance the understanding of tumor heterogeneity regulation in APA and yield new insights into therapeutic target identification in TNBC.

Competitive endogenous RNA is an intrinsic component of EMT regulatory circuits and modulates EMT.

The competitive endogenous RNA (ceRNA) hypothesis suggests an intrinsic mechanism to regulate biological processes. However, whether the dynamic changes of ceRNAs can modulate miRNA activities remains controversial. Here, we examine the dynamics of ceRNAs during TGF-ß-induced epithelial-to-mesenchymal transition (EMT). We observe that TGFBI, a transcript highly induced during EMT in A549 cells, acts as the ceRNA for miR-21 to modulate EMT. We further identify FN1 as the ceRNA for miR-200c in the canonical SNAIL-ZEB-miR200 circuit in MCF10A cells. Experimental assays and computational simulations demonstrate that the dynamically induced ceRNAs are directly coupled with the canonical double negative feedback loops and are critical to the induction of EMT. These results help to establish the relevance of ceRNA in cancer EMT and suggest that ceRNA is an intrinsic component of the EMT regulatory circuit and may represent a potential target to disrupt EMT during tumorigenesis.

Multi-Omics Profiling Reveals Distinct Microenvironment Characterization and Suggests Immune Escape Mechanisms of Triple-Negative Breast Cancer.

PURPOSE: The tumor microenvironment has a profound impact on prognosis and immunotherapy. However, the landscape of the triple-negative breast cancer (TNBC) microenvironment has not been fully understood. EXPERIMENTAL DESIGN: Using the largest original multi-omics dataset of TNBC (n = 386), we conducted an extensive immunogenomic analysis to explore the heterogeneity and prognostic significance of the TNBC microenvironment. We further analyzed the potential immune escape mechanisms of TNBC. RESULTS: The TNBC microenvironment phenotypes were classified into three heterogeneous clusters: cluster 1, the "immune-desert" cluster, with low microenvironment cell infiltration; cluster 2, the "innate immune-inactivated" cluster, with resting innate immune cells and nonimmune stromal cells infiltration; and cluster 3, the "immune-inflamed" cluster, with abundant adaptive and innate immune cells infiltration. The clustering result was validated internally with pathologic sections and externally with The Cancer Genome Atlas and METABRIC cohorts. The microenvironment clusters had significant prognostic efficacy. In terms of potential immune escape mechanisms, cluster 1 was characterized by an incapability to attract immune cells, and MYC amplification was correlated with low immune infiltration. In cluster 2, chemotaxis but inactivation of innate immunity and low tumor antigen burden might contribute to immune escape, and mutations in the PI3K-AKT pathway might be correlated with this effect. Cluster 3 featured high expression of immune checkpoint molecules. CONCLUSIONS: Our study represents a step toward personalized immunotherapy for patients with TNBC. Immune checkpoint inhibitors might be effective for "immune-inflamed" cluster, and the transformation of "cold tumors" into "hot tumors" should be considered for "immune-desert" and "innate immune-inactivated" clusters.

Genomic and Transcriptomic Landscape of Triple-Negative Breast Cancers: Subtypes and Treatment Strategies.

We comprehensively analyzed clinical, genomic, and transcriptomic data of a cohort of 465 primary triple-negative breast cancer (TNBC). PIK3CA mutations and copy-number gains of chromosome 22q11 were more frequent in our Chinese cohort than in The Cancer Genome Atlas. We classified TNBCs into four transcriptome-based subtypes: (1) luminal androgen receptor (LAR), (2) immunomodulatory, (3) basal-like immune-suppressed, and (4) mesenchymal-like. Putative therapeutic targets or biomarkers were identified among each subtype. Importantly, the LAR subtype showed more ERBB2 somatic mutations, infrequent mutational signature 3 and frequent CDKN2A loss. The comprehensive profile of TNBCs provided here will serve as a reference to further advance the understanding and precision treatment of TNBC.

CircRNAFisher: a systematic computational approach for de novo circular RNA identification.

Circular RNAs (circRNAs) are emerging species of mRNA splicing products with largely unknown functions. Although several computational pipelines for circRNA identification have been developed, these methods strictly rely on uniquely mapped reads overlapping back-splice junctions (BSJs) and lack approaches to model the statistical significance of the identified circRNAs. Here, we reported a systematic computational approach to identify circRNAs by simultaneously utilizing BSJ overlapping reads and discordant BSJ spanning reads to identify circRNAs. Moreover, we developed a novel procedure to estimate the P-values of the identified circRNAs. A computational cross-validation and experimental validations demonstrated that our method performed favorably compared to existing circRNA detection tools. We created a standalone tool, CircRNAFisher, to implement the method, which might be valuable to computational and experimental scientists studying circRNAs.

PHF5A Epigenetically Inhibits Apoptosis to Promote Breast Cancer Progression.

Alternative splicing (AS) and its regulation play critical roles in cancer, yet the dysregulation of AS and its molecular bases in breast cancer development have not yet been elucidated. Using an in vivo CRISPR screen targeting RNA-binding proteins, we identified PHD finger protein 5A (PHF5A) as a key splicing factor involved in tumor progression. PHF5A expression was frequently upregulated in breast cancer and correlated with poor survival, and knockdown of PHF5A significantly suppressed cell proliferation, migration, and tumor formation. PHF5A was required for SF3b spliceosome stability and linked the complex to histones, and the PHF5A-SF3b complex modulated AS changes in apoptotic signaling. In addition, expression of a short truncated FAS-activated serine/threonine kinase (FASTK) protein was increased after PHF5A ablation and facilitated Fas-mediated apoptosis. This PHF5A-modulated FASTK-AS axis was widely present in breast cancer specimens, particularly those of the triple-negative subtype. Taken together, our findings reveal that PHF5A serves as an epigenetic suppressor of apoptosis and thus provides a mechanistic basis for breast cancer progression and may be a valuable therapeutic target.Significance: This study provides an epigenetic mechanistic basis for the aggressive biology of breast cancer and identifies a translatable therapeutic target. Cancer Res; 78(12); 3190-206. ©2018 AACR.

Computational identification of mutually exclusive transcriptional drivers dysregulating metastatic microRNAs in prostate cancer.

Androgen-ablation therapies, which are the standard treatment for metastatic prostate cancer, invariably lead to acquired resistance. Hence, a systematic identification of additional drivers may provide useful insights into the development of effective therapies. Numerous microRNAs that are critical for metastasis are dysregulated in metastatic prostate cancer, but the underlying molecular mechanism is poorly understood. We perform an integrative analysis of transcription factor (TF) and microRNA expression profiles and computationally identify three master TFs, AR, HOXC6 and NKX2-2, which induce the aberrant metastatic microRNA expression in a mutually exclusive fashion. Experimental validations confirm that the three TFs co-dysregulate a large number of metastasis-associated microRNAs. Moreover, their overexpression substantially enhances cell motility and is consistently associated with a poor clinical outcome. Finally, the mutually exclusive overexpression between AR, HOXC6 and NKX2-2 is preserved across various tissues and cancers, suggesting that mutual exclusivity may represent an intrinsic characteristic of driver TFs during tumorigenesis.

Synergistic action of master transcription factors controls epithelial-to-mesenchymal transition.

Epithelial-to-mesenchymal transition (EMT) is a complex multistep process in which phenotype switches are mediated by a network of transcription factors (TFs). Systematic characterization of all dynamic TFs controlling EMT state transitions, especially for the intermediate partial-EMT state, represents a highly relevant yet largely unexplored task. Here, we performed a computational analysis that integrated time-course EMT transcriptomic data with public cistromic data and identified three synergistic master TFs (ETS2, HNF4A and JUNB) that regulate the transition through the partial-EMT state. Overexpression of these regulators predicted a poor clinical outcome, and their elimination readily abolished TGF-ß-induced EMT. Importantly, these factors utilized a clique motif, physically interact and their cumulative binding generally characterized EMT-associated genes. Furthermore, analyses of H3K27ac ChIP-seq data revealed that ETS2, HNF4A and JUNB are associated with super-enhancers and the administration of BRD4 inhibitor readily abolished TGF-ß-induced EMT. These findings have implications for systematic discovery of master EMT regulators and super-enhancers as novel targets for controlling metastasis.

3'UTR shortening identifies high-risk cancers with targeted dysregulation of the ceRNA network.

Competing endogenous RNA (ceRNA) interactions form a multilayered network that regulates gene expression in various biological pathways. Recent studies have demonstrated novel roles of ceRNA interactions in tumorigenesis, but the dynamics of the ceRNA network in cancer remain unexplored. Here, we examine ceRNA network dynamics in prostate cancer from the perspective of alternative cleavage and polyadenylation (APA) and reveal the principles of such changes. Analysis of exon array data revealed that both shortened and lengthened 3'UTRs are abundant. Consensus clustering with APA data stratified cancers into groups with differing risks of biochemical relapse and revealed that a ceRNA subnetwork enriched with cancer genes was specifically dysregulated in high-risk cancers. The novel connection between 3'UTR shortening and ceRNA network dysregulation was supported by the unusually high number of microRNA response elements (MREs) shared by the dysregulated ceRNA interactions and the significantly altered 3'UTRs. The dysregulation followed a fundamental principle in that ceRNA interactions connecting genes that show opposite trends in expression change are preferentially dysregulated. This targeted dysregulation is responsible for the majority of the observed expression changes in genes with significant ceRNA dysregulation and represents a novel mechanism underlying aberrant oncogenic expression.

Ponatinib efficiently kills imatinib-resistant chronic eosinophilic leukemia cells harboring gatekeeper mutant T674I FIP1L1-PDGFRα: roles of Mcl-1 and ß-catenin.

BACKGROUND: T674I FIP1L1-PDGFRα in a subset of chronic eosinophilic leukemia (CEL) is a gatekeeper mutation that is resistant to many tyrosine kinase inhibitors (TKIs) (e.g., imatinib, nilotinib and dasatinib), similar to T315I Bcr-Abl. Therefore, novel TKIs effective against T674I FIP1L1-PDGFRα are needed. Ponatinib (AP24534) is a novel orally bioavailable TKI against T315I Bcr-Abl, but it is not clear whether ponatinib is effective against T674I FIP1L1-PDGFRα. The purpose of this study was to examine the effect of ponatinib on T674I FIP1L1-PDGFRα. METHODS: Molecular docking analysis in silico was performed. The effects of ponatinib on PDGFRα signaling pathways, apoptosis and cell cycling were examined in EOL-1, BaF3 cells expressing either wild type (WT) or T674I FIP1L1-PDGFRα. The in vivo antitumor activity of ponatinib was evaluated with xenografted BaF3-T674I FIP1L1-PDGFRα cells in nude mice models. RESULTS: Molecular docking analysis revealed that ponatinib could bind to the DFG (Asp-Phe-Gly)-out state of T674I PDGFRα. Ponatinib potently inhibited the phosphorylation of WT and T674I FIP1L1-PDGFRα and their downstream signaling molecules (e.g., Stat3, Stat5). Ponatinib strikingly inhibited the growth of both WT and T674I FIP1L1-PDGFRα-carrying CEL cells (IC50: 0.004-2.5 nM). It induced apoptosis in CEL cells with caspase-3-dependent cleavage of Mcl-1, and inhibited tyrosine phosphorylation of ß-catenin to decrease its stability and pro-survival functions. In vivo, ponatinib abrogated the growth of xenografted BaF3-T674I FIP1L1-PDGFRα cells in nude mice. CONCLUSIONS: Ponatinib is a pan-FIP1L1-PDGFRα inhibitor, and clinical trials are warranted to investigate its efficacy in imatinib-resistant CEL.

Substituted indolin-2-ones as p90 ribosomal S6 protein kinase 2 (RSK2) inhibitors: Molecular docking simulation and structure-activity relationship analysis.

A series of novel indolin-2-ones inhibitors against p90 ribosomal S6 protein kinase 2 (RSK2) were designed and synthesized and their structure-activity relationship (SAR) was studied. The most potent inhibitor, compound 3s, exhibited potent inhibition against RSK2 with an IC50 value of 0.5 µM and presented a satisfactory selectivity against 23 kinases. The interactions of these inhibitors with RSK2 were investigated based on the proposed binding poses with molecular docking simulation. Four compounds and six compounds exhibited moderate anti-proliferation activities against PC 3 cells and MCF-7 cells, respectively.

Synthesis, activity evaluation, and docking analysis of barbituric acid aryl hydrazone derivatives as RSK2 inhibitors.

The 90 kDa ribosomal S6 kinases (RSKs), especially RSK2, have attracted attention for the development of new anticancer agents. Through structural optimization of the hit compound 1 from our previous study, a series of barbituric acid aryl hydrazone analogues were designed and synthesized as potential RSK2 inhibitors. The most potent one, compound 9, showed a higher activity against RSK2 with an IC50 value of 1.95 µM. To analyze and elucidate their structure-activity relationship, the homology model of RSK2 N-terminal kinase domain was built and molecular docking simulations were performed, which provide helpful clues to design new inhibitors with desired activities.

Insulin-like growth factor-1 receptor (IGF-1R) kinase inhibitors in cancer therapy: advances and perspectives.

The insulin-like growth factors (IGF) and their receptors play pivotal roles in cellular signaling transduction and thus regulate cell growth, differentiation, apoptosis, transformation and other important physiological progresses. The insulin-like growth factor 1 receptor (IGF-1R) mainly engages in the Ras/mitogen-activated protein kinase (MAPK) pathway and the phosphoinositide 3-kinase/protein kinase B (PI3K/AKT) pathway, and also forms cross-talk with the epidermal growth factor receptor (EGFR) pathway. Currently, it draws more attention since its overexpression has been demonstrated in various human cancers, such as colorectal cancer, breast cancer, prostate cancer and lung tumors, thus the strategy targeting the IGF-1R would be promising in treatment of these cancers. There are already dozens of agents developed for the inhibition of IGF-1R, which are categorized into monoclonal antibodies, small molecule inhibitors and so on. While in this review, small molecule inhibitors would be the focus for detailed discussion. Herein, we updated previously reported research papers and reviews in this field and summarized developments of small molecule inhibitors up to 2011. Finally, we proposed the application of network pharmacology methods to reconsider the clinical use of inhibitors with concomitant IR inhibition or other kinases inhibition, hoping that more optimal combinations would be obtained for cancer therapy.

SHAFTS: a hybrid approach for 3D molecular similarity calculation. 2. Prospective case study in the discovery of diverse p90 ribosomal S6 protein kinase 2 inhibitors to suppress cell migration.

We described a prospective application of ligand-based virtual screening program SHAFTS to discover novel inhibitors for p90 ribosomal S6 protein kinase 2 (RSK2). Taking the putative 3D conformations of two weakly binding RSK2 NTKD inhibitors as query templates, SHAFTS was used to perform 3D similarity based virtual screening because of a lack of crystal structure of RSK2 protein, thus leading to the identification of several novel scaffolds that would have been missed by conventional 2D fingerprint methods. The most potent hit compounds show low micromolar inhibitory activities against RSK2. In particular, one of the hit compounds exhibits potent antimigration activity against the MDA-MB-231 tumor cell. The results exemplified SHAFTS' application in active enrichment and scaffold hopping, which is of general interest for lead identification in drug discovery endeavors and also provides novel scaffolds that lay the foundation for uncovering new RSK2 regulatory mechanisms.