Overexpression of lipid metabolism genes and PBX1 in the contralateral breasts of women with estrogen receptor-negative breast cancer

Specific differences and novel key regulatory genes of sex in influencing exceptional longevity phenotypes

BACKGROUND AND AIMS

Although the life expectancy of women systematically and robustly exceeds that of men, specific differences and molecular mechanisms of sex in influencing longevity phenotypes remain largely unknown. Therefore, we performed transcriptome sequencing of peripheral blood samples to explore regulatory mechanisms of healthy longevity by incorporating sex data.

METHODS

We selected 34 exceptional longevity (age: 98.26 ± 2.45 years) and 16 controls (age: 52.81 ± 9.78) without advanced outcomes from 1363 longevity and 692 controls recruited from Nanning of Guangxi for RNA sequencing 1. The transcriptome sequencing 1 data of 50 samples were compared by longevity and sex to screen differentially expressed genes (DEGs). Then, 121 aging samples (40-110 years old) without advanced outcomes from 355 longevity and 294 controls recruited from Dongxing of Guangxi were selected for RNA sequencing 2. The genes associated with aging from the transcriptome sequencing 2 of 121 aging samples were filtered out. Finally, the gender-related longevity candidate genes and their possible metabolic pathways were verified by cell model of aging and a real-time polymerase chain reaction (RT-PCR).

RESULTS

Metabolism differs between male and female and plays a key role in longevity. Moreover, the principal findings of this study revealed a novel key gene, UGT2B11, that plays an important role in regulating lipid metabolism through the peroxisome proliferator activated receptor gamma (PPARG) signalling pathway and ultimately improving lifespan, particularly in females.

CONCLUSION

The findings suggest specific differences in metabolism affecting exceptional longevity phenotypes between the sexes and offer novel therapeutic targets to extend lifespan by regulating lipid homeostasis.

Identification and validation of a glycosyltransferase gene signature as a novel prognostic model for lung adenocarcinoma

Background

The role of glycosyltransferase (GT) genes in lung adenocarcinoma (LUAD) needs further elucidation. Thus, our study aims to identify the prognostic gene signature of LUAD and explore its molecular functions.

Methods

We initially extracted GT gene sets from the database, and obtained mRNA expression levels and clinical data from The Cancer Genome Atlas (TCGA) database. For constructing a prognostic model for GT genes, we utilized univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses. Using the model, patients were categorized into high- and low-risk groups. Additionally, we evaluated differences in tumor immune infiltration between these groups and identified potential therapeutic drugs. Finally, we experimentally validated the expression levels of these crucial prognostic genes.

Results

We developed a risk score comprising nine GT genes (C1GALT1, FUT1, GALNT2, PLOD2, POMK, PYGB, ST3GAL6, UGT2B11, UGT3A1). Patients were then categorized into low- and high-risk groups based on this score. The low-risk group showed superior overall survival (OS) compared to the high-risk group. There were significantly distinct tumor immune microenvironment statuses observed between the two groups. We identified potential therapeutic drugs, including the MEK inhibitor (PD-184352). Finally, we verified the expression of these nine GT genes through immunohistochemistry (IHC) staining and quantitative real-time PCR (qPCR).

Conclusion

We identified a distinct LUAD GT gene signature, and these differentially expressed mRNAs could serve as valuable prognostic biomarkers and therapeutic targets. Furthermore, we experimentally validated their expression levels and identified potential therapeutic agents.

PBX1 as a novel master regulator in cancer: Its regulation, molecular biology, and therapeutic applications

PBX1 is a critical transcription factor at the top of various cell fate-determining pathways. In cancer, PBX1 stands at the crossroads of multiple oncogenic signaling pathways and mediates responses by recruiting a broad repertoire of downstream targets. Research thus far has corroborated the involvement of PBX1 in cancer proliferation, resisting apoptosis, tumor-associated neoangiogenesis, epithelial-mesenchymal transition (EMT) and metastasis, immune evasion, genome instability, and dysregulating cellular metabolism. Recently, our understanding of the functional regulation of the PBX1 protein has advanced, as increasing evidence has depicted a regulatory network consisting of transcriptional, post-transcriptional, and post-translational levels of control mechanisms. Furthermore, accumulating studies have supported the clinical utilization of PBX1 as a prognostic or therapeutic target in cancer. Preliminary results showed that PBX1 entails vast potential as a targetable master regulator in the treatment of cancer, particularly in those with high-risk features and resistance to other therapeutic strategies. In this review, we will explore the regulation, protein-protein interactions, molecular pathways, clinical application, and future challenges of PBX1.

PBX1: a TALE of two seasons-key roles during development and in cancer

Pre-B cell leukemia factor 1 (PBX1) is a Three Aminoacid Loop Extension (TALE) homeodomain-containing transcription factor playing crucial roles in organ pattering during embryogenesis, through the formation of nuclear complexes with other TALE class and/or homeobox proteins to regulate target genes. Its contribution to the development of several organs has been elucidated mainly through the study of murine knockout models. A crucial role for human development has been recently highlighted through the discovery of different de novo pathogenic PBX1 variants in children affected by developmental defects. In the adult, PBX1 is expressed in selected tissues such as in the brain, in the gastro-intestinal and urinary systems, or in hematopoietic stem and progenitor cells, while in other organs is barely detectable. When involved in the t(1;19) chromosomal translocation it acts as an oncogene, since the resulting fusion protein drives pre-B cell leukemia, due to the induction of target genes not normally targeted by the native protein. Its aberrant expression has been associated to tumor development, progression, or therapy-resistance as in breast cancer, ovarian cancer or myeloproliferative neoplasm (MPN). On the other hand, in colorectal cancer PBX1 functions as a tumor suppressor, highlighting its context-dependent role. We here discuss differences and analogies of PBX1 roles during embryonic development and in cancer, focusing mainly on the most recent discoveries.

Exploring tumor-normal cross-talk with TranNet: Role of the environment in tumor progression

There is a growing awareness that tumor-adjacent normal tissues used as control samples in cancer studies do not represent fully healthy tissues. Instead, they are intermediates between healthy tissues and tumors. The factors that contribute to the deviation of such control samples from healthy state include exposure to the tumor-promoting factors, tumor-related immune response, and other aspects of tumor microenvironment. Characterizing the relation between gene expression of tumor-adjacent control samples and tumors is fundamental for understanding roles of microenvironment in tumor initiation and progression, as well as for identification of diagnostic and prognostic biomarkers for cancers. To address the demand, we developed and validated TranNet, a computational approach that utilizes gene expression in matched control and tumor samples to study the relation between their gene expression profiles. TranNet infers a sparse weighted bipartite graph from gene expression profiles of matched control samples to tumors. The results allow us to identify predictors (potential regulators) of this transition. To our knowledge, TranNet is the first computational method to infer such dependencies. We applied TranNet to the data of several cancer types and their matched control samples from The Cancer Genome Atlas (TCGA). Many predictors identified by TranNet are genes associated with regulation by the tumor microenvironment as they are enriched in G-protein coupled receptor signaling, cell-to-cell communication, immune processes, and cell adhesion. Correspondingly, targets of inferred predictors are enriched in pathways related to tissue remodelling (including the epithelial-mesenchymal Transition (EMT)), immune response, and cell proliferation. This implies that the predictors are markers and potential stromal facilitators of tumor progression. Our results provide new insights into the relationships between tumor adjacent control sample, tumor and the tumor environment. Moreover, the set of predictors identified by TranNet will provide a valuable resource for future investigations.

Decreased Expression of GLYATL1 Predicts Poor Prognosis in Patients with Clear Cell Renal Cell Carcinoma

Background

GLYATL1 is a member of the glycine-N-acyltransferase family, which catalyses acyl group transfer. The role of GLYATL1 in cancer is largely unknown; therefore, the potential value of GLYATL1 in clear cell renal cell carcinoma (ccRCC) was explored.

Methods

The ccRCC gene expression profiles and clinical data were obtained from the University of California Santa Cruz Xena platform. Differential expression and survival analysis were performed using R software. Samples from the TIMER public database and real-world were used for validation. The potential molecular mechanism of GLYATL1 in ccRCC was explored using gene set enrichment analysis (GSEA).

Results

GLYATL1 was downregulated, indicating a poor prognosis in ccRCC. Low expression of GLYATL1 was significantly associated with advanced stage and higher histological grade ccRCC. The differential expression of GLYATL1 was validated at the protein level using clinical samples and tissue microarray. The results of GSEA showed that multiple crucial signalling pathways including fatty acid metabolism, adipogenesis, oxidative phosphorylation and epithelial-mesenchymal transition were enriched.

Conclusion

This study demonstrated that GLYATL1 downregulation has an unfavourable impact on the survival of patients with ccRCC. The resulting data indicated that GLYATL1 could be a potential new target for ccRCC therapy, which may be helpful for the personalized treatment of such individuals.

Exploring tumor-normal cross-talk with TranNet: role of the environment in tumor progression

There is a growing awareness that tumor-adjacent normal tissues used as control samples in cancer studies do not represent fully healthy tissues. Instead, they are intermediates between healthy tissues and tumors. The factors that contribute to the deviation of such control samples from healthy state include exposure to the tumor-promoting factors, tumor-related immune response, and other aspects of tumor microenvironment. Characterizing the relation between gene expression of tumor-adjacent control samples and tumors is fundamental for understanding roles of microenvironment in tumor initiation and progression, as well as for identification of diagnostic and prognostic biomarkers for cancers. To address the demand, we developed and validated TranNet, a computational approach that utilizes gene expression in matched control and tumor samples to study the relation between their gene expression profiles. TranNet infers a sparse weighted bipartite graph from gene expression profiles of matched control samples to tumors. The results allow us to identify predictors (potential regulators) of this transition. To our knowledge, TranNet is the first computational method to infer such regulation. We applied TranNet to the data of several cancer types and their matched control samples from The Cancer Genome Atlas (TCGA). Many predictors identified by TranNet are genes associated with regulation by the tumor microenvironment as they are enriched in G-protein coupled receptor signaling, cell-to-cell communication, immune processes, and cell adhesion. Correspondingly, targets of inferred predictors are enriched in pathways related to tissue remodelling (including the epithelial-mesenchymal Transition (EMT)), immune response, and cell proliferation. This implies that the predictors are markers and potential stromal facilitators of tumor progression. Our results provide new insights for the relationships between tumor adjacent control sample, tumor and the tumor environment. Moreover, the set of predictors identified by TranNet will provide a valuable resource for future investigations. The TranNet method was implemented in python, source codes and the data sets used for and generated during this study are available at the Github site https://github.com/ncbi/TranNet .

Comparative Analysis of Transcriptomic Changes including mRNA and microRNA Expression Induced by the Xenoestrogens Zearalenone and Bisphenol A in Human Ovarian Cells

Xenoestrogens are natural or synthetic compounds that mimic the effect of endogenous estrogens and might cause cancer. We aimed to compare the global transcriptomic response to zearalenone (ZEA; mycotoxin) and bisphenol A (BPA; plastic additive) with the effect of physiological estradiol (E2) in the PEO1 human ovarian cell line by mRNA and microRNA sequencing. Estrogen exposure induced remarkable transcriptomic changes: 308, 288 and 63 genes were upregulated (log₂FC > 1); 292, 260 and 45 genes were downregulated (log₂FC < -1) in response to E2 (10 nM), ZEA (10 nM) and BPA (100 nM), respectively. Furthermore, the expression of 13, 11 and 10 miRNAs changed significantly (log₂FC > 1, or log₂FC < -1) after exposure to E2, ZEA and BPA, respectively. Functional enrichment analysis of the significantly differentially expressed genes and miRNAs revealed several pathways related to the regulation of cell proliferation and migration. The effect of E2 and ZEA was highly comparable: 407 genes were coregulated by these molecules. We could identify 83 genes that were regulated by all three treatments that might have a significant role in the estrogen response of ovarian cells. Furthermore, the downregulation of several miRNAs (miR-501-5p, let-7a-2-3p, miR-26a-2-3p, miR-197-5p and miR-582-3p) was confirmed by qPCR, which might support the proliferative effect of estrogens in ovarian cells.

Spindle function and Wnt pathway inhibition by PBX1 to suppress tumor progression via downregulating DCDC2 in colorectal cancer

PBX1 is a transcription factor that regulates a variety of genes, involved in intracellular lipid metabolism, cell proliferation, and other pathways. The promoting and inhibiting function of PBX1 in various cancer types was extensively discussed, however, there have been no studies on PBX1 proteins in colorectal cancer (CRC). This study aimed to reveal the anti-tumor function of PBX1 in CRC and the underlying molecular mechanism. Bioinformatics analysis revealed that PBX1 is downregulated in CRC, indicating that is a potential antioncogene in CRC. Overexpression of PBX1 suppresses tumor growth and metastasis in vitro and in vivo. Mechanistically, we found that PBX1 acted as a transcription factor that suppressed DCDC2 expression and inhibited spindle function. Moreover, the PBX1-DCDC2 axis controlled the Wnt pathway in CRC cells. Overexpression of DCDC2 restored CRC proliferation, metastasis abilities and Wnt pathway. In conclusion, this study suggests that PBX1 acts as a transcription factor to suppress DCDC2 expression and inhibit cell proliferation and metastasis by disrupting spindle function and the Wnt pathway in CRC.

HMGCS2 Mediation of Ketone Levels Affects Sorafenib Treatment Efficacy in Liver Cancer Cells

Primary liver cancer is the fifth leading death of cancers in men, and hepatocellular carcinoma (HCC) accounts for approximately 90% of all primary liver cancer cases. Sorafenib is a first-line drug for advanced-stage HCC patients. Sorafenib is a multi-target kinase inhibitor that blocks tumor cell proliferation and angiogenesis. Despite sorafenib treatment extending survival, some patients experience side effects, and sorafenib resistance does occur. 3-Hydroxymethyl glutaryl-CoA synthase 2 (HMGCS2) is the rate-limiting enzyme for ketogenesis, which synthesizes the ketone bodies, β-hydroxybutyrate (β-HB) and acetoacetate (AcAc). β-HB is the most abundant ketone body which is present in a 4:1 ratio compared to AcAc. Recently, ketone body treatment was found to have therapeutic effects against many cancers by causing metabolic alternations and cancer cell apoptosis. Our previous publication showed that HMGCS2 downregulation-mediated ketone body reduction promoted HCC clinicopathological progression through regulating c-Myc/cyclin D1 and caspase-dependent signaling. However, whether HMGCS2-regulated ketone body production alters the sensitivity of human HCC to sorafenib treatment remains unclear. In this study, we showed that HMGCS2 downregulation enhanced the proliferative ability and attenuated the cytotoxic effects of sorafenib by activating expressions of phosphorylated (p)-extracellular signal-regulated kinase (ERK), p-P38, and p-AKT. In contrast, HMGCS2 overexpression decreased cell proliferation and enhanced the cytotoxic effects of sorafenib in HCC cells by inhibiting ERK activation. Furthermore, we showed that knockdown HMGCS2 exhibited the potential migratory ability, as well as decreasing zonula occludens protein (ZO)-1 and increasing c-Myc expression in both sorafenib-treated Huh7 and HepG2 cells. Although HMGCS2 overexpression did not alter the migratory effect, expressions of ZO-1, c-Myc, and N-cadherin decreased in sorafenib-treated HMGCS2-overexpressing HCC cells. Finally, we investigated whether ketone treatment influences sorafenib sensitivity. We showed that β-HB pretreatment decreased cell proliferation and enhanced antiproliferative effect of sorafenib in both Huh7 and HepG2 cells. In conclusion, this study defined the impacts of HMGCS2 expression and ketone body treatment on influencing the sorafenib sensitivity of liver cancer cells.

Arachidonate 15-lipoxygenase type B: Regulation, function, and its role in pathophysiology

As a lipoxygenase (LOX), arachidonate 15-lipoxygenase type B (ALOX15B) peroxidizes polyenoic fatty acids (PUFAs) including arachidonic acid (AA), eicosapentaenoic acid (EPA), docosahexaenoic acid (DHA), and linoleic acid (LA) to their corresponding fatty acid hydroperoxides. Distinctive to ALOX15B, fatty acid oxygenation occurs with positional specificity, catalyzed by the non-heme iron containing active site, and in addition to free PUFAs, membrane-esterified fatty acids serve as substrates for ALOX15B. Like other LOX enzymes, ALOX15B is linked to the formation of specialized pro-resolving lipid mediators (SPMs), and altered expression is apparent in various inflammatory diseases such as asthma, psoriasis, and atherosclerosis. In primary human macrophages, ALOX15B expression is associated with cellular cholesterol homeostasis and is induced by hypoxia. Like in inflammation, the role of ALOX15B in cancer is inconclusive. In prostate and breast carcinomas, ALOX15B is attributed a tumor-suppressive role, whereas in colorectal cancer, ALOX15B expression is associated with a poorer prognosis. As the biological function of ALOX15B remains an open question, this review aims to provide a comprehensive overview of the current state of research related to ALOX15B.

Development and Validation of a Prognostic Classifier Based on Lipid Metabolism-Related Genes for Breast Cancer

Background

The changes of lipid metabolism have been implicated in the development of many tumors, but its role in breast invasive carcinoma (BRCA) remains to be fully established. Here, we attempted to ascertain the prognostic value of lipid metabolism-related genes in BRCA.

Methods

We obtained RNA expression data and clinical information for BRCA and normal samples from public databases and downloaded a lipid metabolism-related gene set. Ingenuity Pathway Analysis (IPA) was applied to identify the potential pathways and functions of Differentially Expressed Genes (DEGs) related to lipid metabolism. Subsequently, univariate and multivariate Cox regression analyses were utilized to construct the prognostic gene signature. Functional enrichment analysis of prognostic genes was achieved by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Kaplan-Meier analysis, Receiver Operating Characteristic (ROC) curves, clinical follow-up results were employed to assess the prognostic potency. Potential compounds targeting prognostic genes were screened by Connectivity Map (CMap) database and a prognostic gene-drug interaction network was constructed using Comparative Toxicogenomics Database (CTD). Furthermore, we separately validated the selected marker genes in BRCA samples and human breast cancer cell lines (MCF-7, MDA-MB-231).

Results

IPA and functional enrichment analysis demonstrated that the 162 lipid metabolism-related DEGs we obtained were involved in many lipid metabolism and BRCA pathological signatures. The prognostic classifier we constructed comprising SDC1 and SORBS1 can serve as an independent prognostic marker for BRCA. CMap filtered 37 potential compounds against prognostic genes, of which 16 compounds could target both two prognostic genes were identified by CTD. The functions of the two prognostic genes in breast cancer cells were verified by cell function experiments.

Conclusion

Within this study, we identified a novel prognostic classifier based on two lipid metabolism-related genes: SDC1 and SORBS1. This result highlighted a new perspective on the metabolic exploration of BRCA.

Lipid exposure activates gene expression changes associated with estrogen receptor negative breast cancer

Improved understanding of local breast biology that favors the development of estrogen receptor negative (ER-) breast cancer (BC) would foster better prevention strategies. We have previously shown that overexpression of specific lipid metabolism genes is associated with the development of ER- BC. We now report results of exposure of MCF-10A and MCF-12A cells, and mammary organoids to representative medium- and long-chain polyunsaturated fatty acids. This exposure caused a dynamic and profound change in gene expression, accompanied by changes in chromatin packing density, chromatin accessibility, and histone posttranslational modifications (PTMs). We identified 38 metabolic reactions that showed significantly increased activity, including reactions related to one-carbon metabolism. Among these reactions are those that produce S-adenosyl-L-methionine for histone PTMs. Utilizing both an in-vitro model and samples from women at high risk for ER- BC, we show that lipid exposure engenders gene expression, signaling pathway activation, and histone marks associated with the development of ER- BC.

Epigenetic Determinants of Racial Disparity in Breast Cancer: Looking beyond Genetic Alterations

Simple Summary

A substantial disparity in breast cancer incidence and mortality exists between African American (AA) and European American (EA) women. However, the basis for these disparities is poorly understood. In this article, we describe that gene–environment interactions mediated through epigenetic modifications may play a significant role in racial disparities in BC incidence and outcomes. Our in silico analyses and an in-depth literature survey suggest that there exists a significant difference in epigenetic patterns between AA and EA women with breast cancer. Herein, we describe the environmental factors that contribute to these epigenetic changes, which may underlie the disparate racial burden in patients with breast cancer. We suggest that AA women with higher basal epigenetic changes, may have higher pre-disposition to cancer onset, and an aggressive disease course. Pre-existing racial differences in epigenetic profiles of breast tissues raises the possibility of examining these profiles for early diagnosis.

Abstract

Breast cancer (BC) is the most commonly diagnosed cancer in women. Despite advancements in BC screening, prevention, and treatment, BC incidence and mortality remain high among African American (AA) women. Compared with European American (EA) women, AA women tend to be diagnosed with more advanced and aggressive tumors and exhibit worse survival outcomes. Most studies investigating the determinants of racial disparities in BC have focused on genetic factors associated with African ancestry. However, various environmental and social stressors over an individual’s life course can also shape racial stratification in BC. These social and environmental exposures result in long-term changes in gene expression mediated by epigenetic mechanisms. Epigenetics is often portrayed as an intersection of socially patterned stress and genetic expression. The enduring nature of epigenetic changes makes them suitable for studying the effects of different environmental exposures over an individual’s life course on gene expression. The role of differential social and environmental exposures in racial disparities in BC suggests varied epigenetic profiles or signatures associated with specific BC subtypes in AA and EA women. These epigenetic profiles in EA and AA women could be used as biomarkers for early BC diagnosis and disease prognosis and may prove valuable for the development of targeted therapies for BC. This review article discusses the current state of knowledge regarding epigenetic differences between AA and EA women with BC. We also discuss the role of socio-environmental factors, including psychosocial stress, environmental toxicants, and dietary factors, in delineating the different epigenetic profiles in AA and EA patients with BC.

m6 A-mediated regulation of PBX1-GCH1 axis promotes gastric cancer proliferation and metastasis by elevating tetrahydrobiopterin levels

BACKGROUND

Methyltransferase 3 (METTL3)-mediated N6-methyladenosine (m6 A) RNA modification has been demonstrated to be a potential factor in promoting gastric cancer (GC). METTL3 regulates a series of signaling pathways by modifying various mRNAs. This study aimed to identify novel METTL3-mediated signaling pathways and explored possible targets for use in the clinical setting of gastric cancer.

METHODS

To investigate the proliferation and metastatic capacity of GC cell lines with METTL3 knockdown, a xenograft, lung metastasis, and popliteal lymph node metastasis model was used. The m6 A-modified RNA immunoprecipitation (Me-RIP) sequence was utilized to explore the target mRNAs of METTL3. Cell counting kit 8 and transwell assays were performed to investigate the promoting function of pre-B cell leukemia homeobox 1 (PBX1) and GTP cyclohydrolase 1 (GCH1). Western blotting and chromatin immunoprecipitation were employed to confirm the involvement of the METTL3-PBX1-GCH1 axis. ELISA and liquid chromatography-mass spectrometry were used to explore the biological function of tetrahydrobiopterin (BH4 ).

RESULTS

Knockdown of METTL3 suppressed xenograft tumor growth and lung/lymph node metastasis in vivo. Mechanistically, we found that METTL3 combined with and stabilized PBX1 mRNAs. Chromatin immunoprecipitation (ChIP) and further experiments suggested that PBX1 acted as a transcription factor inducing GCH1 expression. Moreover, the METTL3-PBX1-GCH1 axis increased BH4 levels in GC cells, thereby promoting tumor progression.

CONCLUSIONS

This study suggested that METTL3 enzymes promote tumor growth and lung/lymph node metastasis via METTL3-PBX1-GCH1 axis increasing BH4 levels in GC.

Systems medicine dissection of chr1q-amp reveals a novel PBX1-FOXM1 axis for targeted therapy in multiple myeloma

Understanding the biological and clinical impact of copy number aberrations (CNAs) on the development of precision therapies in cancer remains an unmet challenge. Genetic amplification of chromosome 1q (chr1q-amp) is a major CNA conferring an adverse prognosis in several types of cancer, including in the blood cancer multiple myeloma (MM). Although several genes across chromosome 1 (chr1q) portend high-risk MM disease, the underpinning molecular etiology remains elusive. Here, with reference to the 3-dimensional (3D) chromatin structure, we integrate multi-omics data sets from patients with MM with genetic variables to obtain an associated clinical risk map across chr1q and to identify 103 adverse prognosis genes in chr1q-amp MM. Prominent among these genes, the transcription factor PBX1 is ectopically expressed by genetic amplification and epigenetic activation of its own preserved 3D regulatory domain. By binding to reprogrammed superenhancers, PBX1 directly regulates critical oncogenic pathways and a FOXM1-dependent transcriptional program. Together, PBX1 and FOXM1 activate a proliferative gene signature that predicts adverse prognosis across multiple types of cancer. Notably, pharmacological disruption of the PBX1-FOXM1 axis with existing agents (thiostrepton) and a novel PBX1 small molecule inhibitor (T417) is selectively toxic against chr1q-amp myeloma and solid tumor cells. Overall, our systems medicine approach successfully identifies CNA-driven oncogenic circuitries, links them to clinical phenotypes, and proposes novel CNA-targeted therapy strategies in MM and other types of cancer.

The mechanism of action of Fangji Huangqi Decoction on epithelial-mesenchymal transition in breast cancer using high-throughput next-generation sequencing and network pharmacology

ETHNOPHARMACOLOGICAL RELEVANCE

Fangji Huangqi Decoction (FHD) is widely used in traditional Chinese medicine (TCM). FHD has been hypothesized to inhibit the epithelial-mesenchymal transition (EMT) process, which may positively impact breast cancer prevention and treatment. However, its exact mechanism of action is still unknown.

AIM OF THE STUDY

This study aimed to screen potential targets of FHD for the treatment of EMT in breast cancer through network pharmacology, and to verify their therapeutic effects in vitro experiments and high-throughput second-generation sequencing.

MATERIALS AND METHODS

The data sets of effective components and targets of FHD were established through the Traditional Chinese Medicine Systems Pharmacology database. The GeneCards and OMIM databases were used to establish breast cancer-related target datasets, which were then matched with the TCM target data. The interaction between key target proteins was analyzed using the STRING database; the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were used to identify the associated biological processes and enriched signal pathways, respectively. The active ingredient disease target network was analyzed using Cytoscape. Finally, next generation sequencing was used to verify the related pathways of FHD intervention in EMT in breast cancer. High-content screening was used to identify the genes/pathways affected by FHD. MDA-MB-231 and HCC-1937 breast cancer cell lines were used to evaluate the impact of FHD on migration, invasion, and EMT.

RESULTS

Eighty possible significant targets were identified for the treatment of breast cancer EMT with FHD; GO and KEGG were used to identify 173 cell biological processes associated with breast cancer (P < 0.05), including the NF-κB and PI3K-Akt signaling pathways. The high-throughput sequencing and network pharmacology results were highly consistent. The migration and invasion ability of MDA-MB-231 cells was reduced and their EMT status could be reversed by DSHR2 knockdown. The results of morphology and scratch assays showed that FHD could improve the EMT status of HCC-1973.

CONCLUSIONS

This study provides more evidence to support the clinical application of FHD, which has reliable interventional effects on breast cancer EMT. Its therapeutic effects may involve a multi-target, multi-pathway, and multi-mechanism effect.

Transposons: Unexpected players in cancer

Transposons are repetitive DNA sequences encompassing about half of the human genome. They play a vital role in genome stability maintenance and contribute to genomic diversity and evolution. Their activity is regulated by various mechanisms considering the deleterious effects of these mobile elements. Various genetic risk factors and environmental stress conditions affect the regulatory pathways causing alteration of transposon expression. Our knowledge of the biological role of transposons is limited especially in various types of cancers. Retrotransposons of different types (LTR-retrotransposons, LINEs and SINEs) regulate a plethora of genes that have a role in cell reprogramming, tumor suppression, cell cycle, apoptosis, cell adhesion and migration, and DNA repair. The regulatory mechanisms of transposons, their deregulation and different mechanisms underlying transposon-mediated carcinogenesis in humans focusing on the three most prevalent types, lung, breast and colorectal cancers, were reviewed. The modes of regulation employed include alternative splicing, deletion, insertion, duplication in genes and promoters resulting in upregulation, downregulation or silencing of genes.

PBX1: a key character of the hallmarks of cancer

Pre-B-cell leukemia homeobox transcription factor 1 (PBX1) was first identified as part of a fusion protein resulting from the chromosomal translocation t(1;19) in pre-B cell acute lymphoblastic leukemias. Since then, PBX1 has been associated with important developmental programs, and its expression dysregulation has been related to multifactorial disorders, including cancer. As PBX1 overexpression in many cancers is correlated to poor prognosis, we sought to understand how this transcription factor contributes to carcinogenesis, and to organize PBX1's roles in the hallmarks of cancer. There is enough evidence to associate PBX1 with at least five hallmarks: sustaining proliferative signaling, activating invasion and metastasis, inducing angiogenesis, resisting cell death, and deregulating cellular energetics. The lack of studies investigating a possible role for PBX1 on the remaining hallmarks made it impossible to defend or refute its contribution on them. However, the functions of some of the PBX1's transcription targets indicate a potential engagement of PBX1 in the avoidance of immune destruction and in the tumor-promoting inflammation hallmarks. Interestingly, PBX1 might be a player in tumor suppression by activating the transcription of some DNA damage response genes. This is the first review organizing PBX1 roles into the hallmarks of cancer. Thus, we encourage future studies to uncover the PBX1's underlying mechanisms to promote carcinogenesis, for it is a promising diagnostic and prognostic biomarker, as well as a potential target in cancer treatment.

HMGCS2 in metabolic pathways was associated with overall survival in hepatocellular carcinoma: A LASSO-derived study

This integrated bioinformatic study aimed to investigate potential prognostic candidates in hepatocellular carcinoma (HCC). In the GSE14520, GSE101685, and The Cancer Genome Atlas (TCGA) datasets, differentially expressed genes (DEGs) were identified and functional pathways of common DEGs were enriched. The least absolute shrinkage and selection operator (LASSO) model was used to screen the potential parameters associated with overall survival (OS) in HCC patients. Metabolic pathways were the most significantly enriched functional pathways of common DEGs in these three datasets. After LASSO model analysis, HMGCS2, UGP2, BCLC staging and TNM staging were screened as potential prognostic candidates for OS in HCC patients in GSE14520. HMGCS2 in the metabolic pathway was significantly downregulated in tumor tissues and peripheral blood mononuclear cells in HCC patients (all p < 0.05). Cox regression model indicated that HMGCS2 might be associate with OS in HCC patients in GSE14520 and in the TCGA (p = 0.029 and p = 0.05, respectively). Kaplan-Meier analysis demonstrated that HMGCS2 downregulation in tumors contributed to an unfavorable OS in HCC patients, both in GSE14520 and in the TCGA (p = 0.0001 and p = 0.0002, respectively). Additionally, HMGCS2 was significantly downregulated in HCC patients with high alpha-fetoprotein (AFP), main tumor size >5 cm, multinodular, advanced tumor staging including BCLC, TNM and CLIP (all p < 0.05). HMGCS2 was involved in metabolic pathways, and downregulated HMGCS2 in tumors was associated with unfavorable OS in HCC patients.

Identification of Ten Core Hub Genes as Potential Biomarkers and Treatment Target for Hepatoblastoma

Background

This study aimed to systematically investigate gene signatures for hepatoblastoma (HB) and identify potential biomarkers for its diagnosis and treatment.

Materials and Methods

GSE131329 and GSE81928 were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between hepatoblastoma and normal samples were identified using the Limma package in R. Then, the similarity of network traits between two sets of genes was analyzed by weighted gene correlation network analysis (WGCNA). Cytoscape was used to visualize and select hub genes. PPI network of hub genes was construed by Cytoscape. GO enrichment and KEGG pathway analyses of hub genes were carried out using ClueGO. The random forest classifier was constructed based on the hub genes using the GSE131329 dataset as the training set, and its reliability was validated using the GSE81928 dataset. The resulting core hub genes were combined with the InnateDB database to identify the innate core genes.

Results

A total of 4244 DEGs in HB were identified. WGCNA identified four modules that were significantly correlated with the disease status. A total of 114 hub genes were obtained within the top 20 genes of each node rank. 6982 relation pairs and 3700 nodes were contained in the PPI network of 114 hub genes. GO enrichment and KEGG pathway analyses of hub genes were focused on MAPK, cell cycle, p53, and other crucial pathways involved in HB. A random forest classifier was constructed using the 114 hub genes as feature genes, resulting in a 95.5% true positive rate when classifying HB and normal samples. A total of 35 core hub genes were obtained through the mean decrease in accuracy and mean decrease Gini of the random forest model. The classification efficiency of the random forest model was 81.4%. Finally, CDK1, TOP2A, ADRA1A, FANCI, XRCC1, TPX2, CCNB2, CDK4, GLYATL1, and CFHR3 were identified by cross-comparison with the InnateDB database.

Conclusion

Our study established a random forest classifier that identified 10 core genes in HB. These findings may be beneficial for the diagnosis, prediction, and targeted therapy of HB.

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.

Emerging role of pioneer transcription factors in targeted ERα positive breast cancer

Transcription factors (TFs) are modular protein groups that preferably bind to DNA sequences and guide genomic expression through transcription. Among these key regulators, “pioneer factors” are an emerging class of TFs that specifically interact with nucleosomal DNA and facilitate accessible genomic binding sites for the additional TFs. There is growing evidence of these specialized modulators in particular malignancies, as highlighted by agents’ clinical efficacy, specifically targeting nuclear hormone receptors. They have been implicated in multiple cancers more recently, with a high proportion inculpating on hormone influential cancers. Moreover, extended crosstalk and cooperation between ERα pioneering factors in estrogen-dependent breast cancer (BC) remain elucidated. This review discusses on the recent advances in our understanding of pioneer TFs in cancer, especially highlighting its potentiality to modulate chromatin condensation to permit ERα recruitment in BC cells. Through the study it was concluded that the highly prospected pioneer TFs in BC, including FOXA1, TLE1, PBX1, and GATA3, possess the potential therapeutic significance and further innovations in the field could yield targeted therapy in cancer treatment.

The expression and prognostic value of GLYATL1 and its potential role in hepatocellular carcinoma

BACKGROUND

Glycine-N-acyltransferase-like 1 (GLYATL1), which is involved in the detoxification of endogenous and exogenous acyl-CoA, promotes glutamine metabolism in xenobiotic metabolism. Recent evidence suggests an association between GLYATL1 and tumors. However, there are few comprehensive analyses of GLYATL1 in cancers. We evaluated the expression and prognostic value of GLYATL1 and explored the mechanism underlying the association between GLYATL1 and cancers.

METHODS

GLYATL1 mRNA expression across cancers was investigated in the Oncomine database and confirmed in the UALCAN and Gene Expression Profiling Interactive Analysis (GEPIA) databases. Next, its prognostic value in different cancers was revealed by PrognoScan and Kaplan-Meier plotter. According to clinicopathologic features, we conducted a subgroup analysis of the prognosis of GLYATL1 in a cohort of hepatocellular carcinoma (HCC) patients from The Cancer Genome Atlas (TCGA) and the GSE116174 dataset. We further investigated the GLYATL1 promoter methylation profile in HCC. Next, a protein-protein interaction (PPI) network was constructed via the Search Tool for the Retrieval of Interacting Genes (STRING) database. Finally, we utilized gene set enrichment analysis (GSEA) to identify significantly enriched pathways and confirmed their associations using the Tumor Immune Estimation Resource (TIMER) and GEPIA databases.

RESULTS

GLYATL1 is downregulated in many cancers and indicates a poor prognosis. Specifically, low GLYATL1 expression was associated with short overall survival (OS) in HCC patients. Interestingly, GLYATL1 expression was associated with poor OS in stage I-II HCC patients and was revealed as an independent prognostic factor. The promoter methylation level of GLYATL1 in HCC tissue was significantly higher than that in normal liver tissue. The PPI network suggested that GLYATL1 is co-expressed with ten genes, including CNGA3 and GNB5. GSEA revealed that GLYATL1 is predominantly negatively enriched in xenobiotic metabolism, and the gene association analysis in TIMER and GEPIA showed a significantly negative association between the expression of GLYATL1 and the expression of most genes involved in mitochondrial glutamine metabolism, including SLC1A5 and SLC1A11.

CONCLUSIONS

Our study is the first to shed light on the expression and prognostic value of GLYATL1 in cancers and provide a potential regulatory mechanism underlying HCC development.

PBX1 Increases the Radiosensitivity of Oesophageal Squamous Cancer by Targeting of STAT3

The radioresistance of oesophageal squamous cell carcinoma (OSCC) is a critical factor leading to a poor prognosis among patients. The expression of PBX1 is abnormally high in a broad range of human tissues, and this gene plays a key role in tumour proliferation. This research intended to explore the radiosensitization of OSCC by silencing PBX1. The OSCC cell lines KYSE450 and KYSE150 were subjected to PBX1 silencing and/or irradiation (IR). Cell proliferation, colony formation, and apoptosis were tested to evaluate the radiosensitization ability of PBX1 silencing. The levels of STAT3 and p-STAT3 in the OSCC cells were tested by Western blotting. Furthermore, KYSE150 cells with or without PBX1 silencing were xenografted into nude mice with or without radiation exposure. Concomitant PBX1 silencing and IR can obviously suppress growth and enhance radiosensitivity in OSCC cells and xenografts. Moreover, the downregulation of PBX1 inhibits the expression of STAT3 and p-STAT3. The downregulation of PBX1 may increase radiosensitivity in OSCC cells and xenografts via the PBX1/STAT3 pathway. Our findings demonstrate that PBX1 may be a potential target for promoting the effect of radiation therapy in OSCC patients.

Coexpression of PBX1 and EMP2 as Prognostic Biomarkers in Estrogen Receptor-Negative Breast Cancer via Data Mining

Previous studies revealed that PBX1 ranked the third in the differentially expressed genes about development and progression of breast cancer (BC). Nevertheless, the role of PBX1 contributing to progression of BC has been unevaluated. Here, on the basis of ONCOMINE and GOBO databases, we compared BC samples with normal controls about the expression of PBX1 in various types of cancers, as well as their related expression levels in cancer cell lines by Cancer Cell Line Encyclopedia (CCLE) analysis. It was also found that, when compared with normal controls, PBX1 was markedly higher expressed not only in BC samples but also in BC cell lines, and coexpressed with EMP2 by ONCOMINE and CCLE coexpression analysis, which was also expressed higher in BC samples and BC cell lines similarly. According to Kaplan-Meier plotter, we further explored the prognostic functions of PBX1 and EMP2 in different molecular subtypes of BC, respectively. We demonstrated that overexpression of PBX1 mRNA was correlated with worse survival in luminal B subtype BC, whereas increased EMP2 expression was associated with shorter relapse-free survival in estrogen receptor (ER)-negative patients. Combining with previous studies, we could make a conclusion that coexpression of PBX1 and EMP2 predicts poor prognosis in ER-negative BC, which could be effective biomarkers for BC.

Development and Clinical Validation of a Novel 4-Gene Prognostic Signature Predicting Survival in Colorectal Cancer

In this study, we collected genes related to energy metabolism, used gene expression data from public databases to classify molecular subtypes of colon cancer (COAD) based on the genes related to energy metabolism, and further evaluated the relationships between the molecular subtypes and prognosis and clinical characteristics. Differential expression analysis of the molecular subtypes yielded 1948 differentially expressed genes (DEGs), whose functions were closely related to the occurrence and development of cancer. Based on the DEGs, we constructed a 4-gene prognostic risk model and identified the high expression of FOXD4, ENPEP, HOXC6, and ALOX15B as a risk factor associated with a high risk of developing COAD. The 4-gene signature has strong robustness and a stable predictive performance in datasets from different platforms not only in patients with early COAD but also in all patients with colon cancer. The enriched pathways of the 4-gene signature in the high- and low-risk groups obtained by GSEA were significantly related to the occurrence and development of colon cancer. Moreover, the results of qPCR, immunohistochemistry staining and Western blot assay revealed that FOXD4, ENPEP, HOXC6, and ALOX15B are over expressed in CRC tissues and cells. These results suggesting that the signature could potentially be used as a prognostic marker for clinical diagnosis.

Computational models applied to metabolomics data hints at the relevance of glutamine metabolism in breast cancer

Background

Metabolomics has a great potential in the development of new biomarkers in cancer and it has experiment recent technical advances.

Methods

In this study, metabolomics and gene expression data from 67 localized (stage I to IIIB) breast cancer tumor samples were analyzed, using (1) probabilistic graphical models to define associations using quantitative data without other a priori information; and (2) Flux Balance Analysis and flux activities to characterize differences in metabolic pathways.

Results

On the one hand, both analyses highlighted the importance of glutamine in breast cancer. Moreover, cell experiments showed that treating breast cancer cells with drugs targeting glutamine metabolism significantly affects cell viability. On the other hand, these computational methods suggested some hypotheses and have demonstrated their utility in the analysis of metabolomics data and in associating metabolomics with patient’s clinical outcome.

Conclusions

Computational analyses applied to metabolomics data suggested that glutamine metabolism is a relevant process in breast cancer. Cell experiments confirmed this hypothesis. In addition, these computational analyses allow associating metabolomics data with patient prognosis.

HMGCS2 Mediates Ketone Production and Regulates the Proliferation and Metastasis of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the most common primary malignant tumor worldwide; however, the traditional therapeutic approaches and survival rates are still limited. To improve current therapies, it is necessary to investigate the molecular mechanisms underlying liver cancer and to identify potential therapeutic targets. The aims of this study were to verify the mechanisms and therapeutic potential of the ketogenesis rate-limiting enzyme 3-Hydroxymethylglutaryl-CoA synthase 2 (HMGCS2) in HCC. Immunohistochemical staining of human liver disease tissue arrays showed that HMGCS2 is abundantly expressed in normal liver tissues but is downregulated in cirrhosis and HCC tissues. In HCC patients, lower HMGCS2 expression was correlated with higher pathological grades and clinical stages. In our investigation of the molecular mechanisms of HMGCS2 in HCC, we showed that knockdown of HMGCS2 decreased ketone production, which promoted cell proliferation, cell migration, and xenograft tumorigenesis by enhancing c-Myc/cyclinD1 and EMT signaling and by suppressing the caspase-dependent apoptosis pathway. Ketone body treatment reduced the proliferation- and migration-promoting effects of HMGCS2 knockdown in cells. In contrast, HMGCS2 overexpression increased the intracellular ketone level and inhibited cell proliferation, cell migration, and xenograft tumorigenesis. Finally, ketogenic diet administration significantly inhibited liver cancer cell growth in mice. Our studies highlight the potential therapeutic strategy of targeting HMGCS2-mediated ketogenesis in liver cancer.

Characterization of glycine-N-acyltransferase like 1 (GLYATL1) in prostate cancer

BACKGROUND

Recent microarray and sequencing studies of prostate cancer showed multiple molecular alterations during cancer progression. It is critical to evaluate these molecular changes to identify new biomarkers and targets. We performed analysis of glycine-N-acyltransferase like 1 (GLYATL1) expression in various stages of prostate cancer in this study and evaluated the regulation of GLYATL1 by androgen.

METHOD

We performed in silico analysis of cancer gene expression profiling and transcriptome sequencing to evaluate GLYATL1 expression in prostate cancer. Furthermore, we performed immunohistochemistry using specific GLYATL1 antibody using high-density prostate cancer tissue microarray containing primary and metastatic prostate cancer. We also tested the regulation of GLYATL1 expression by androgen and ETS transcription factor ETV1. In addition, we performed RNA-sequencing of GLYATL1 modulated prostate cancer cells to evaluate the gene expression and changes in molecular pathways.

RESULTS

Our in silico analysis of cancer gene expression profiling and transcriptome sequencing we revealed an overexpression of GLYATL1 in primary prostate cancer. Confirming these findings by immunohistochemistry, we show that GLYATL1 is overexpressed in primary prostate cancer compared with metastatic prostate cancer and benign prostatic tissue. Low-grade cancers had higher GLYATL1 expression compared to high-grade prostate tumors. Our studies showed that GLYATL1 is upregulated upon androgen treatment in LNCaP prostate cancer cells which harbors ETV1 gene rearrangement. Furthermore, ETV1 knockdown in LNCaP cells showed downregulation of GLYATL1 suggesting potential regulation of GLYATL1 by ETS transcription factor ETV1. Transcriptome sequencing using the GLYATL1 knockdown prostate cancer cell lines LNCaP showed regulation of multiple metabolic pathways.

CONCLUSIONS

In summary, our study characterizes the expression of GLYATL1 in prostate cancer and explores the regulation of its regulation in prostate cancer showing role for androgen and ETS transcription factor ETV1. Future studies are needed to decipher the biological significance of these findings.

KRAS mutation and epithelial-macrophage interplay in pancreatic neoplastic transformation

Pancreatic ductal adenocarcinoma (PDA) is characterized by epithelial mutations in KRAS and prominent tumor-associated inflammation, including macrophage infiltration. But knowledge of early interactions between neoplastic epithelium and macrophages in PDA carcinogenesis is limited. Using a pancreatic organoid model, we found that the expression of mutant KRAS in organoids increased (i) ductal to acinar gene expression ratios, (ii) epithelial cells proliferation and (iii) colony formation capacity in vitro, and endowed pancreatic cells with the ability to generate neoplastic tumors in vivo. KRAS mutations induced a protumorigenic phenotype in macrophages. Altered macrophages decreased epithelial pigment epithelial derived factor (PEDF) expression and induced a cancerous phenotype. We validated our findings using annotated patient samples from The Cancer Genome Atlas (TCGA) and in our human PDA specimens. Epithelium-macrophage cross-talk occurs early in pancreatic carcinogenesis where KRAS directly induces cancer-related phenotypes in epithelium, and also promotes a protumorigenic phenotype in macrophages, in turn augmenting neoplastic growth.

Emerging roles and mechanisms of FOXC2 in cancer

Forkhead box protein C2 (FOXC2), a transcription factor of the forkhead/winged-helix family, is required for embryonic and prenatal development. FOXC2 acts as a crucial modulator during both angiogenesis and lymphangiogenesis via multiple angiogenic and lymphangiogenic pathways, respectively. Although recent studies have shed light on the emerging role of FOXC2 in cancer, very little is known about the precise underlying mechanisms. The purpose of this review is to summarize the current understanding of FOXC2 and provide potential mechanistic explanations of the relationship between FOXC2 and cancer, as well as discuss the prospect for future research in the promising prognostic value of FOXC2 in cancer.

DHRS2 inhibits cell growth and motility in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is highly prevailing in Asia and it is ranked in the most aggressive squamous cell carcinomas. High-frequency loss of heterozygosity occurred in chromosome 14q11.2 in many tumors including ESCC, suggesting that one or more tumor-suppressor genes might exist within this region. In this study, we identified the tumor-suppressing role of DHRS2 (short-chain dehydrogenase/reductase family, member 2) at 14q11.2 in ESCCs. Downregulation of DHRS2 occurred in 30.8% of primary ESCC tumor tissues vs paired non-tumorous tissues. DHRS2 downregulation was associated significantly with ESCC invasion, lymph nodes metastasis and clinical staging (P<0.001). Survival analysis revealed that DHRS2 downregulation was significantly associated with worse outcome of patients with ESCC. In vitro and in vivo studies indicated that both DHRS2 variants could suppress cell proliferation and cell motility. Moreover, we demonstrated that DHRS2 could reduce reactive oxygen species and decrease nicotinamide adenine dinucleotide phosphate (oxidized/reduced), increase p53 stability and decrease Rb phosphorylation; it also decreased p38 mitogen-activated protein kinase phosphorylation and matrix metalloproteinase 2. In summary, these findings demonstrated that DHRS2 had an important part in ESCC development and progression.

miR-107-mediated decrease of HMGCS2 indicates poor outcomes and promotes cell migration in hepatocellular carcinoma

Mitochondrial 3-hydroxy-3-methylglutaryl-CoA synthase (HMGCS2) has been implicated in human cancers, but its role and clinical significance in hepatocellular carcinoma (HCC) remain unknown. Here, we show that HMGCS2 is downregulated and exhibits antimetastatic potential in HCC. Low expression of HMGCS2 was associated with poor tumor differentiation, vascular invasion and worse overall and disease-free survivals in two independent cohorts consisting of 743 cases. In vitro data demonstrated HMGCS2 overexpression suppressed, whereas HMGCS2 silence promoted HCC cell migration via Epithelial-Mesenchymal Transition (EMT) process and the activation of ERK/c-Jun signaling pathway. Inhibition of ERK phosphorylation by PD098059 markedly attenuated the malignant phenotypes mediated by HMGCS2 siRNA. Furthermore, miR-107 was identified as an upstream regulator of HMGCS2 via directly targeting the 3'-UTR of HMGCS2 mRNA. Collectively, our findings suggest HMGCS2 serve as a promising prognostic biomarker and exert anti-tumor activity towards HCC, and therefore provide a potential target for HCC clinical intervention.