Search

File containing Supplementary methods and figure legends

Understanding the intratumoral heterogeneity of hepatocellular carcinoma is instructive for developing personalized therapy and identifying molecular biomarkers. Here we applied whole-exome sequencing to 69 samples from 11 patients to resolve the genetic architecture of subclonal diversification. Spatial genomic diversity was found in all 11 hepatocellular carcinoma cases, with 29% of driver mutations being heterogeneous, including TERT, ARID1A, NOTCH2, and STAG2. Similar with other cancer types, TP53 mutations were always shared between all tumor regions, that is, located on the “trunk” of the evolutionary tree. In addition, we found that variants within several drug targets such as KIT, SYK, and PIK3CA were mutated in a fully clonal manner, indicating their therapeutic potentials for hepatocellular carcinoma. Temporal dissection of mutational signatures suggested that mutagenic processes associated with exposure to aristolochic acid and aflatoxin might play a more important role in early, as opposed to late, stages of hepatocellular carcinoma development. Moreover, we observed extensive intratumoral epigenetic heterogeneity in hepatocellular carcinoma based on multiple independent analytical methods and showed that intratumoral methylation heterogeneity might play important roles in the biology of hepatocellular carcinoma cells. Our results also demonstrated prominent heterogeneity of intratumoral methylation even in a stable hepatocellular carcinoma genome. Together, these findings highlight widespread intratumoral heterogeneity at both the genomic and epigenomic levels in hepatocellular carcinoma and provide an important molecular foundation for better understanding the pathogenesis of this malignancy. Cancer Res; 77(9); 2255–65. ©2017 AACR.

Files containing results related to methylation ITH

File containing all somatic mutations

Files containing resutls related to genomic and methylation ITH , detailed mutational signatures, and sample pathology

Files containing sequencing details and sample information

MEX3A is an RNA-binding protein that mediates mRNA decay through binding to 3′ untranslated regions. However, its role and mechanism in clear cell renal cell carcinoma remain unknown. In this study, we found that MEX3A expression was transcriptionally activated by ETS1 and upregulated in clear cell renal cell carcinoma. Silencing MEX3A markedly reduced clear cell renal cell carcinoma cell proliferation in vitro and in vivo. Inhibiting MEX3A induced G1/S cell-cycle arrest. Gene set enrichment analysis revealed that E2F targets are the central downstream pathways of MEX3A. To identify MEX3A targets, systematic screening using enhanced cross-linking and immunoprecipitation sequencing, and RNA-immunoprecipitation sequencing assays were performed. A network of 4,000 genes was identified as potential targets of MEX3A. Gene ontology analysis of upregulated genes bound by MEX3A indicated that negative regulation of the cell proliferation pathway was highly enriched. Further assays indicated that MEX3A bound to the CDKN2B 3′ untranslated region, promoting its mRNA degradation. This leads to decreased levels of CDKN2B and an uncontrolled cell cycle in clear cell renal cell carcinoma, which was confirmed by rescue experiments. Our findings revealed that MEX3A acts as a post-transcriptional regulator of abnormal cell-cycle progression in clear cell renal cell carcinoma., Graphical abstract, We described the comprehensive downstream targets of MEX3A using RIP-seq and eCLIP-seq. The data revealed that oncogene MEX3A was transcriptionally activated and mainly regulated G1/S transition in ccRCC, thus providing a novel molecular basis for clinical diagnosis and treatment.

To reconstruct systematically hyperactive transcription factor (TF)-dependent transcription networks in squamous cell carcinomas (SCCs), a computational method (ELMER) was applied to 1293 pan-SCC patient samples, and 44 hyperactive SCC TFs were identified. As a top candidate, DLX5 exhibits a notable bifurcate re-configuration of its bivalent promoter in cancer. Specifically, DLX5 maintains a bivalent state in normal tissues; its promoter is hypermethylation, leading to DLX5 transcriptional silencing in esophageal adenocarcinoma (EAC). In stark contrast, DLX5 promoter gains active histone marks and becomes transcriptionally activated in ESCC, which is directly mediated by SOX2. Functionally, silencing of DLX5 substantially inhibits SCC viability both in vitro and in vivo. Mechanistically, DLX5 cooperates with TP63 in regulating ∼2000 enhancers and promoters, which converge on activating cancer-promoting pathways. Together, our data establish a novel and strong SCC-promoting factor and elucidate a new epigenomic mechanism - bifurcate chromatin re-configuration - during cancer development.

Rationale: The forkhead box A1 (FOXA1) is a crucial transcription factor in initiation and development of breast, lung and prostate cancer. Previous studies about the FOXA1 transcriptional network were mainly focused on protein-coding genes. Its regulatory network of long non-coding RNAs (lncRNAs) and their role in FOXA1 oncogenic activity remains unknown. Methods: The Cancer Genome Atlas (TCGA) data, RNA-seq and ChIP-seq data were used to analyze FOXA1 regulated lncRNAs. RT-qPCR was used to detect the expression of DSCAM-AS1, RT-qPCR and Western blotting were used to determine the expression of FOXA1, estrogen receptor α (ERα) and Y box binding protein 1 (YBX1). RNA pull-down and RIP-qPCR were employed to investigate the interaction between DSCAM-AS1 and YBX1. The effect of DSCAM-AS1 on malignant phenotypes was examined through in vitro and in vivo assays. Results: In this study, we conducted a global analysis of FOXA1 regulated lncRNAs. For detailed analysis, we chose lncRNA DSCAM-AS1, which is specifically expressed in lung adenocarcinoma, breast and prostate cancer. The expression level of DSCAM-AS1 is regulated by two super-enhancers (SEs) driven by FOXA1. High expression levels of DSCAM-AS1 was associated with poor prognosis. Knockout experiments showed DSCAM-AS1 was essential for the growth of xenograft tumors. Moreover, we demonstrated DSCAM-AS1 can regulate the expression of the master transcriptional factor FOXA1. In breast cancer, DSCAM-AS1 was also found to regulate ERα. Mechanistically, DSCAM-AS1 interacts with YBX1 and influences the recruitment of YBX1 in the promoter regions of FOXA1 and ERα. Conclusion: Our study demonstrated that lncRNA DSCAM-AS1 was transcriptionally activated by super-enhancers driven by FOXA1 and exhibited lineage-specific expression pattern. DSCAM-AS1 can promote cancer progression by interacting with YBX1 and regulating expression of FOXA1 and ERα.

Hepatocellular carcinoma (HCC) is one of the most common malignancies worldwide. SET and MYND domain-containing protein 3 (SMYD3) has been shown to promote the progression of various types of human cancers, including liver cancer; however, the detailed molecular mechanism is still largely unknown. Here, we report that SMYD3 expression in HCC is an independent prognostic factor for survival and promotes the proliferation and migration of HCC cells. We observed that SMYD3 upregulated sphingosine-1-phosphate receptor 1 (S1PR1) promoter activity by methylating histone 3 (H3K4me3). S1PR1 was expressed at high levels in HCC samples, and high S1PR1 expression was associated with shorter survival. S1PR1 expression was also positively correlated with SMYD3 expression in HCC samples. We confirmed that SMYD3 promotes HCC cell growth and migration in vitro and in vivo by upregulating S1PR1 expression. Further investigations revealed that SMYD3 affects critical signaling pathways associated with the progression of HCC through S1PR1. These findings strongly suggest that SMYD3 has a crucial function in HCC progression that is partially mediated by histone methylation at the downstream gene S1PR1, which affects key signaling pathways associated with carcinogenesis and the progression of HCC.

Hepatocellular carcinoma (HCC) is one of the most common malignancies worldwide. SET and MYND domain-containing protein 3 (SMYD3) has been shown to promote the progression of various types of human cancers, including liver cancer; however, the detailed molecular mechanism is still largely unknown. Here, we report that SMYD3 expression in HCC is an independent prognostic factor for survival and promotes the proliferation and migration of HCC cells. We observed that SMYD3 upregulated sphingosine-1-phosphate receptor 1 (S1PR1) promoter activity by methylating histone 3 (H3K4me3). S1PR1 was expressed at high levels in HCC samples, and high S1PR1 expression was associated with shorter survival. S1PR1 expression was also positively correlated with SMYD3 expression in HCC samples. We confirmed that SMYD3 promotes HCC cell growth and migration in vitro and in vivo by upregulating S1PR1 expression. Further investigations revealed that SMYD3 affects critical signaling pathways associated with the progression of HCC through S1PR1. These findings strongly suggest that SMYD3 has a crucial function in HCC progression that is partially mediated by histone methylation at the downstream gene S1PR1, which affects key signaling pathways associated with carcinogenesis and the progression of HCC.

CCAAT/enhancer binding proteins (CEBPs, including CEBPA, CEBPB, CEBPD, CEBPE, CEBPG, and CEBPZ) play critical roles in a variety of physiological and pathological processes. However, the molecular characteristics and biological significance of CEBPs in esophageal squamous cell carcinoma (ESCC) have rarely been reported. Here, we show that most of the CEBPs are upregulated and accompanied with copy number amplifications in ESCC. Of note, high CEBPG expression is regulated by the ESCC specific transcription factor TP63 and serves as a prognostic factor for poor survival in ESCC patients. Functionally, CEBPG significantly promotes the proliferation and migration of ESCC cells both in vitro and in vivo. Mechanistically, CEBPG activates the PI3K-AKT signaling pathway through directly binding to distal enhancers and/or promoters of genes involved in this pathway, including genes of CCND1, MYC, CDK2, etc. These findings provide new insights into CEBPs dysregulation in ESCC and elucidate a crucial role for CEBPG in the progression of ESCC, highlighting its potential therapeutic value for ESCC treatment.

Liver fibrosis is an increasing health problem worldwide, for which no effective antifibrosis drugs are available. Although the involvement of aerobic glycolysis in hepatic stellate cell (HSC) activation has been reported, the role of pyruvate kinase M2 (PKM2) in liver fibrogenesis still remains unknown. We examined PKM2 expression and location in liver tissues and primary hepatic cells. The in vitro and in vivo effects of a PKM2 antagonist (shikonin) and its allosteric agent (TEPP-46) on liver fibrosis were investigated in HSCs and liver fibrosis mouse model. Chromatin immunoprecipitation sequencing and immunoprecipitation were performed to identify the relevant molecular mechanisms. PKM2 expression was significantly up-regulated in both mouse and human fibrotic livers compared with normal livers, and mainly detected in activated, rather than quiescent, HSCs. PKM2 knockdown markedly inhibited the activation and proliferation of HSCs in vitro. Interestingly, the PKM2 dimer, rather than the tetramer, induced HSC activation. PKM2 tetramerization induced by TEPP-46 effectively inhibited HSC activation, reduced aerobic glycolysis, and decreased MYC and CCND1 expression via regulating histone H3K9 acetylation in activated HSCs. TEPP-46 and shikonin dramatically attenuated liver fibrosis in vivo. Our findings demonstrate a nonmetabolic role of PKM2 in liver fibrosis. PKM2 tetramerization or suppression could prevent HSC activation and protects against liver fibrosis.

As the second most common malignant bone tumor in children and adolescents, Ewing sarcoma is initiated and exacerbated by a chimeric oncoprotein, most commonly, EWS-FLI1. In this study, we apply epigenomic analysis to characterize the transcription dysregulation in this cancer, focusing on the investigation of super-enhancer and its associated transcriptional regulatory mechanisms. We demonstrate that super-enhancer-associated transcripts are significantly enriched in EWS-FLI1 target genes, contribute to the aberrant transcriptional network of the disease, and mediate the exceptional sensitivity of Ewing sarcoma to transcriptional inhibition. Through integrative analysis, we identify MEIS1 as a super-enhancer-driven oncogene, which co-operates with EWS-FLI1 in transcriptional regulation, and plays a key pro-survival role in Ewing sarcoma. Moreover, APCDD1, another super-enhancer-associated gene, acting as a downstream target of both MEIS1 and EWS-FLI1, is also characterized as a novel tumor-promoting factor in this malignancy. These data delineate super-enhancer-mediated transcriptional deregulation in Ewing sarcoma, and uncover numerous candidate oncogenes which can be exploited for further understanding of the molecular pathogenesis for this disease.

Background Exosomes from cancer cells or immune cells, carrying bio-macromolecules or microRNAs (miRNAs), participate in tumor pathogenesis and progression by modulating microenvironment. Our study aims to investigate the role of these microRNA-501-3p (miR-501-3p) containing exosomes derived from tumor-associated macrophage (TAM) in the progression of pancreatic ductal adenocarcinoma (PDAC). Methods Firstly, the function of TAM recruitment in PDAC tissues was assessed, followed by identification of the effects of M2 macrophage-derived exosomes on PDAC cell activities and tumor formation and metastasis in mice. In silico analysis was conducted to predict differentially expressed genes and regulatory miRNAs related to PDAC treated with macrophages, which determined miR-501-3p and TGFBR3 for subsequent experiments. Next, gain- and loss-of-function experiments were performed to examine their role in PDAC progression with the involvement of the TGF-β signaling pathway. Results TAM recruitment in PDAC tissues was associated with metastasis. Highly expressed miR-501-3p was observed in PDAC tissues and TAM-derived exosomes. Both M2 macrophage-derived exosomes and miR-501-3p promoted PDAC cell migration and invasion, as well as tumor formation and metastasis in nude mice. MiR-501-3p was verified to target TGFBR3. PDAC cells presented with down-regulated TGFBR3, which was further decreased in response to M2 macrophage treatment. TGF-β signaling pathway activation was implicated in the promotion of miR-501-3p in PDAC development. The suppression of macrophage-derived exosomal miR-501-3p resulted in the inhibition of tumor formation and metastasis in vivo. Conclusion M2 macrophage-derived exosomal miR-501-3p inhibits tumor suppressor TGFBR3 gene and facilitates the development of PDAC by activating the TGF-β signaling pathway, which provides novel targets for the molecular treatment of PDAC. Electronic supplementary material The online version of this article (10.1186/s13046-019-1313-x) contains supplementary material, which is available to authorized users.

The bromodomain-containing protein 7 (BRD7) is a tumour suppressor protein with critical roles in cell cycle transition and transcriptional regulation. Whether BRD7 is regulated by post-translational modifications remains poorly understood. Here, we find that chemotherapy-induced DNA damage leads to the rapid degradation of BRD7 in various cancer cell lines. PARP-1 binds and poly(ADP)ribosylates BRD7, which enhances its ubiquitination and degradation through the PAR-binding E3 ubiquitin ligase RNF146. Moreover, the PARP1 inhibitor Olaparib significantly enhances the sensitivity of BRD7-positive cancer cells to chemotherapeutic drugs, while it has little effect on cells with low BRD7 expression. Taken together, our findings show that PARP1 induces the degradation of BRD7 resulting in cancer cell resistance to DNA-damaging agents. BRD7 might thus serve as potential biomarker in clinical trial for the prediction of synergistic effects between chemotherapeutic drugs and PARP inhibitors.

Esophageal squamous cell carcinoma (ESCC) is a common and aggressive malignancy, with hitherto dismal clinical outcome. Genomic analyses of patient samples reveal a complex heterogeneous landscape for ESCC, which presents in both intertumor and intratumor forms, manifests at both genomic and epigenomic levels, and contributes significantly to tumor evolution, drug resistance, and metastasis. Here, we review the important molecular characteristics underlying ESCC heterogeneity, with an emphasis on genomic aberrations and their functional contribution to cancer evolutionary trajectories. We further discuss how novel experimental tools, including single-cell sequencing and three-dimensional organoids, may advance our understanding of tumor heterogeneity. Lastly, we suggest that deciphering the mechanisms governing tumor heterogeneity holds the potential to developing precision therapeutics for ESCC patients.

Understanding the intratumoral heterogeneity of hepatocellular carcinoma is instructive for developing personalized therapy and identifying molecular biomarkers. Here we applied whole-exome sequencing to 69 samples from 11 patients to resolve the genetic architecture of subclonal diversification. Spatial genomic diversity was found in all 11 hepatocellular carcinoma cases, with 29% of driver mutations being heterogeneous, including TERT, ARID1A, NOTCH2, and STAG2. Similar with other cancer types, TP53 mutations were always shared between all tumor regions, that is, located on the “trunk” of the evolutionary tree. In addition, we found that variants within several drug targets such as KIT, SYK, and PIK3CA were mutated in a fully clonal manner, indicating their therapeutic potentials for hepatocellular carcinoma. Temporal dissection of mutational signatures suggested that mutagenic processes associated with exposure to aristolochic acid and aflatoxin might play a more important role in early, as opposed to late, stages of hepatocellular carcinoma development. Moreover, we observed extensive intratumoral epigenetic heterogeneity in hepatocellular carcinoma based on multiple independent analytical methods and showed that intratumoral methylation heterogeneity might play important roles in the biology of hepatocellular carcinoma cells. Our results also demonstrated prominent heterogeneity of intratumoral methylation even in a stable hepatocellular carcinoma genome. Together, these findings highlight widespread intratumoral heterogeneity at both the genomic and epigenomic levels in hepatocellular carcinoma and provide an important molecular foundation for better understanding the pathogenesis of this malignancy. Cancer Res; 77(9); 2255–65. ©2017 AACR.

Studies have reported that the CCN family of proteins plays an important role in stimulating tumorigenesis. However, the relationship between the CCN protein family members and the features of hepatocellular carcinoma (HCC) remains unclear. The objective of this study was to determine the relationship between the expression levels of CCN protein family members and the features of HCC. Expression levels of the CCN family of proteins in 80-paired primary HCC samples and 11 normal liver samples were determined by a quantitative real-time PCR assay. Enhanced expression of nephroblastoma overexpressed protein (NOV) and decreased expression of Wnt-induced secreted protein 1 (WISP1), cysteine-rich protein 61 (CYR61) and connective tissue growth factor (CTGF) were found in HCC samples when compared to levels in matched non-cancerous tissues. No significant difference in WISP2 was found between matched-pair samples; only a few samples showed WISP3 expression. Furthermore, the expression levels of NOV, WISP1 and CYR61 were closely correlated with certain clinical features, including venous invasion, cellular differentiation, pTNM stage, disease-free survival and overall survival. Our results suggest that HCC progression may be enhanced by NOV and suppressed by WISP1 and CYR61. Our statistical analysis suggests that these proteins may be valuable in determining the prognosis of this deadly disease and directs attention to modulating the levels of these proteins as a potential mode of therapy.

Proteasome inhibitors have been proven to be effective anticancer compounds in many tumor models, including glioblastoma multiforme (GBM). In this study, we found that the proteasome inhibitor Velcade (PS-341/bortezomib) caused GBM cell death while simultaneously activating the PI3K/Akt pathway. Therefore, we sought to investigate if the PI3K inhibitor ZSTK474 would enhance the effectiveness of Velcade in anticancer therapy. Two GBM cell lines were used to detect the effects of Velcade and ZSTK474 alone or in combination in vitro. The combination of Velcade and ZSTK474 synergistically inhibited the proliferation of GBM cell lines. Cell apoptosis was increased when exposed to Velcade and ZSTK474 in combination as shown by Annexin V analysis. Treatment with both drugs led to downregulation of the p-Akt, p-4EBP1 and p-mTOR proteins as determined by western blot analysis. The anticancer ability of Velcade for glioblastoma multiforme was, therefore, enhanced by combination with the PI3K pathway inhibitor ZSTK474 in glioblastoma multiforme.