Effect of lentivirus-mediated shRNA inactivation of HK1, HK2, and HK3 genes in colorectal cancer and melanoma cells

Transcriptional regulation and post-translational modifications in the glycolytic pathway for targeted cancer therapy

Cancer cells largely rely on aerobic glycolysis or the Warburg effect to generate essential biomolecules and energy for their rapid growth. The key modulators in glycolysis including glucose transporters and enzymes, e.g. hexokinase 2, enolase 1, pyruvate kinase M2, lactate dehydrogenase A, play indispensable roles in glucose uptake, glucose consumption, ATP generation, lactate production, etc. Transcriptional regulation and post-translational modifications (PTMs) of these critical modulators are important for signal transduction and metabolic reprogramming in the glycolytic pathway, which can provide energy advantages to cancer cell growth. In this review we recapitulate the recent advances in research on glycolytic modulators of cancer cells and analyze the strategies targeting these vital modulators including small-molecule inhibitors and microRNAs (miRNAs) for targeted cancer therapy. We focus on the regulation of the glycolytic pathway at the transcription level (e.g., hypoxia-inducible factor 1, c-MYC, p53, sine oculis homeobox homolog 1, N6-methyladenosine modification) and PTMs (including phosphorylation, methylation, acetylation, ubiquitination, etc.) of the key regulators in these processes. This review will provide a comprehensive understanding of the regulation of the key modulators in the glycolytic pathway and might shed light on the targeted cancer therapy at different molecular levels.

Exploring the prognostic value of HK3 and its association with immune infiltration in glioblastoma multiforme

Background: Hexokinase 3 (HK3) is one of the key enzymes involved in glucose phosphorylation (the first step in most glucose metabolic pathways). Many studies have demonstrated the vital role of dysregulation of HK3 in several tumors. However, there is a need for in-depth characterization of the role of HK3 in glioblastoma multiforme (GBM). Methods: All data were sourced from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA). Kaplan-Meier analysis and univariate regression were applied for survival analysis. Gene set enrichment analysis (GSEA) was used for enrichment analysis. Tumor Immune Single Cell Hub (TISCH) database was applied for single-cell analysis. Tumor Immune Dysfunction and Exclusion (TIDE) analysis was applied to evaluate the immune response. Results: HK3 expression was upregulated in GBM and correlated with poor prognosis. The high HK3 expression group was primarily enriched in adaptive immune response, chemokine signaling pathway, and cytokine-cytokine receptor interaction. The high HK3 expression group showed significantly greater enrichment of the majority of immune cells and immune-related pathways. HK3 showed significant correlation with most immune cells, especially macrophages (p < .001, R = .81). TISCH analysis showed that HK3 was predominantly expressed in macrophages in most cancers. HK3 showed significant correlation with most immune-related genes, such as PD-1 (p < .001, R = .41), PDL-1 (p < .001, R = .27), and CTLA-4 (p < .001, R = .29). TIDE analysis revealed that the low HK3 expression group has a lower TIDE score and may benefit from immunotherapy. Drug sensitivity analysis showed that patients with high HK3 expression frequently showed drug resistance. Conclusion: HK3 was associated with poor prognosis and may serve as a biomarker of macrophages in GBM. HK3 was also associated with immune response and drug resistance. Our findings may provide novel insights for GBM immunotherapy.

3-Bromopyruvate Suppresses the Malignant Phenotype of Vemurafenib-Resistant Melanoma Cells

(1) BRAF mutations are associated with high mortality and are a substantial factor in therapeutic decisions. Therapies targeting BRAF-mutated tumors, such as vemurafenib (PLX), have significantly improved the overall survival of melanoma patients. However, patient relapse and low response rates remain challenging, even with contemporary therapeutic alternatives. Highly proliferative tumors often rely on glycolysis to sustain their aggressive phenotype. 3-bromopyruvate (3BP) is a promising glycolysis inhibitor reported to mitigate resistance in tumors. This study aimed to evaluate the potential of 3BP as an antineoplastic agent for PLX-resistant melanoma treatment. (2) The effect of 3BP alone or in combination with PLX on viability, proliferation, colony formation, cell death, migration, invasion, epithelial-mesenchymal marker and metabolic protein expression, extracellular glucose and lactate, and reactive species were evaluated in two PLX-resistant melanoma cell lines. (3) 3BP treatment, which was more effective as monotherapy than combined with PLX, disturbed the metabolic and epithelial-mesenchymal profile of PLX-resistant cells, impairing their proliferation, migration, and invasion and triggering cell death. (4) 3BP monotherapy is a potent metabolic-disrupting agent against PLX-resistant melanomas, supporting the suppression of the malignant phenotype in this type of neoplasia.

Long non-coding RNAs play an important regulatory role in tumorigenesis and tumor progression through aerobic glycolysis

Compared to normal cells, cancer cells generate ATP mainly through aerobic glycolysis, which promotes tumorigenesis and tumor progression. Long non-coding RNAs (LncRNAs) are a class of transcripts longer than 200 nucleotides with little or without evident protein-encoding function. LncRNAs are involved in the ten hallmarks of cancer, interestingly, they are also closely associated with aerobic glycolysis. However, the mechanism of this process is non-transparent to date. Demonstrating the mechanism of lncRNAs regulating tumorigenesis and tumor progression through aerobic glycolysis is particularly critical for cancer therapy, and may provide novel therapeutic targets or strategies in cancer treatment. In this review, we discuss the role of lncRNAs and aerobic glycolysis in tumorigenesis and tumor progression, and further explore their interaction, in hope to provide a novel therapeutic target for cancer treatment.

The Regulatory Network of MicroRNA in the Metabolism of Colorectal Cancer

Colorectal cancer (CRC) is the third most common malignant tumor in the world. During the progression of CRC, the entire metabolic network undergoes reprogramming, including marked changes in the regulation of glucose, lipid and amino acid metabolism. Although microRNAs (miRNAs) account for only 1% of the entire human genome, they play an important role in almost all physiological and pathological processes in the body. MiRNAs can react directly with key enzymes in the metabolic processes. MiRNAs also interact with other ncRNAs, as a member of non-coding RNA (ncRNA), to form their own regulatory network in various oncogenic pathways of CRC metabolism. The progression of colorectal cancer is closely related to the intestinal flora, where miRNAs act as important mediators. Understanding how miRNAs act in the regulatory network of CRC metabolism is helpful to elucidate the characteristics of tumor occurrence, proliferation, metastasis and drug resistance. This review summarizes the role of miRNAs in the metabolism of CRC and how miRNAs interact with key enzymes, ncRNA and intestinal flora to further discuss how miRNAs affect CRC and realize some new strategies for the early diagnosis and treatment of CRC.

Gene-gene interaction analysis incorporating network information via a structured Bayesian approach

Increasing evidence has shown that gene-gene interactions have important effects in biological processes of human diseases. Due to the high dimensionality of genetic measurements, interaction analysis usually suffers from a lack of sufficient information and has unsatisfactory results. Biological network information has been massively accumulated, allowing researchers to identify biomarkers while taking a system perspective, conducting network selection (of functionally related biomarkers), and accommodating network structures. In main-effect-only analysis, network information has been incorporated. However, effort has been limited in interaction analysis. Recently, link networks that describe the relationships between genetic interactions have been demonstrated as effective for revealing multiscale hierarchical organizations in networks and providing interesting findings beyond node networks. In this study, we develop a novel structured Bayesian interaction analysis approach to effectively incorporate network information. This study is among the first to identify gene-gene interactions with the assistance of network selection, while simultaneously accommodating the underlying network structures of both main effects and interactions. It innovatively respects multiple hierarchies among main effects, interactions, and networks. The Bayesian technique is adopted, which may be more informative for estimation and prediction over some other techniques. An efficient variational Bayesian expectation-maximization algorithm is developed to explore the posterior distribution. Extensive simulation studies demonstrate the practical superiority of the proposed approach. The analysis of TCGA data on melanoma and lung cancer leads to biologically sensible findings with satisfactory prediction accuracy and selection stability.

Comprehensive Analysis of Monocarboxylate Transporter 4 (MCT4) expression in breast cancer prognosis and immune infiltration via integrated bioinformatics analysis

Lactate blunts the anticancer immune response in breast cancer (BC). However, little is known about the exact effect of lactate transporters such as monocarboxylate transporter 4 (MCT4) on immunotherapy. In this study, we investigated the expression status and prognostic value of MCT4 in BC through large-scale transcriptome data. Our results showed that MCT4 was overexpressed in BC, particularly in the basal-like molecular subtype. Overexpression of MCT4 was significantly correlated with high BC lesion grade and poor prognosis. Enrichment analysis indicated that the MCT4-related genes were involved in immune- and metabolism-related bioprocesses, such as myeloid leukocyte activation, the adaptive immune system, and catabolic process. We also found that the expression of MCT4 in BC lesions was associated with immune cell infiltration and glycolytic rate-limiting enzymes like pyruvate kinase M2 (PKM2) and hexokinases-3 (HK3). Our observations indicate that MCT4 may play a pivotal role in the maintenance of the tumor immune microenvironment (TIME) through metabolic reprogramming. The enzymes of the glycolysis pathway (MCT4, PKM2, and HK3) may thus serve as new targets to modulate the TIME and enhance immunotherapy efficiency.[Figure: see text].

Genetic Mutations and Non-Coding RNA-Based Epigenetic Alterations Mediating the Warburg Effect in Colorectal Carcinogenesis

Colorectal cancer (CRC) development is a gradual process defined by the accumulation of numerous genetic mutations and epigenetic alterations leading to the adenoma-carcinoma sequence. Despite significant advances in the diagnosis and treatment of CRC, it continues to be a leading cause of cancer-related deaths worldwide. Even in the presence of oxygen, CRC cells bypass oxidative phosphorylation to produce metabolites that enable them to proliferate and survive-a phenomenon known as the "Warburg effect". Understanding the complex glucose metabolism in CRC cells may support the development of new diagnostic and therapeutic approaches. Here we discuss the most recent findings on genetic mutations and epigenetic modulations that may positively or negatively regulate the Warburg effect in CRC cells. We focus on the non-coding RNA (ncRNA)-based epigenetics, and we present a perspective on the therapeutic relevance of critical molecules and ncRNAs mediating the Warburg effect in CRC cells. All the relevant studies were identified and assessed according to the genes and enzymes mediating the Warburg effect. The findings summarized in this review should provide a better understanding of the relevance of genetic mutations and the ncRNA-based epigenetic alterations to CRC pathogenesis to help overcome chemoresistance.

Mitochondrial Voltage-Dependent Anion Channel 1-Hexokinase-II Complex-Targeted Strategy for Melanoma Inhibition Using Designed Multiblock Peptide Amphiphiles

Targeted therapies of melanoma are of urgent need considering the resistance of this aggressive type of cancer to chemotherapeutics. The voltage-dependent anion channel 1 (VDAC1)-hexokinase-II (HK-II) complex is an emerging target for novel anticancer therapies based on induced mitochondria-mediated apoptosis. The low cell membrane permeability of the anticancer 12-mer peptide N-Ter (RDVFTKGYGFGL) derived from the N-terminal fragment of the VDAC1 protein impedes the intracellular targeting. Here, novel multiblock VDAC1-derived cationic amphiphilic peptides (referred to as Pal-N-Ter-TAT, pFL-N-Ter-TAT, and Pal-pFL-N-Ter-TAT) are designed with a self-assembly propensity and cell-penetrating properties. The created multiblock amphiphilic peptides of partial α-helical conformations form nanoparticles of ellipsoid-like shapes and are characterized by enhanced cellular uptake. The amphiphilic peptides can target mitochondria and dissociate the VDAC1-HK-II complex at the outer mitochondrial membrane, which result in mitochondria-mediated apoptosis. The latter is associated with decrease of the mitochondrial membrane potential, cytochrome c release, and changes of the expression levels of the apoptotic proteins in A375 melanoma cells. Importantly, the mitochondrial VDAC1-derived amphiphilic peptides have a comparable IC₅₀ value for melanoma cells to a small-molecule drug, sorafenib, which has been previously used in clinical trials for melanoma. These results demonstrate the potential of the designed peptide constructs for efficient melanoma inhibition.

Hexokinase 3 dysfunction promotes tumorigenesis and immune escape by upregulating monocyte/macrophage infiltration into the clear cell renal cell carcinoma microenvironment

Purpose: This study aimed to identify the potential prognostic role of HK3 and provide clues about glycolysis and the microenvironmental characteristics of ccRCC. Methods: Based on the Cancer Genome Atlas (TCGA, n = 533) and Gene expression omnibus (GEO) (n = 127) databases, real-world (n = 377) ccRCC cohorts, and approximately 15,000 cancer samples, the prognostic value and immune implications of HK3 were identified. The functional effects of HK3 in ccRCC were analyzed in silico and in vitro. Results: The large-scale findings suggested a significantly higher HK3 expression in ccRCC tissues and the predictive efficacy of HK3 for tumor progression and a poor prognosis. Next, the subgroup survival and Cox regression analyses showed that HK3 serves as a promising and independent predictive marker for the prognosis and survival of patients with ccRCC from bioinformatic databases and real-world cohorts. Subsequently, we found that HK3 could be used to modulate glycolysis and the malignant behaviors of ccRCC cells. The comprehensive results suggested that HK3 is highly correlated with the abundance of immune cells, and specifically stimulates the infiltration of monocytes/macrophages presenting surface markers, regulates the immune checkpoint molecules PD-1 and CTLA-4 of exhaustive T cells, restrains the immune escape of tumor cells, and prompts the immune-rejection microenvironment of ccRCC. Conclusion: In conclusion, the large-scale data first revealed that HK3 could affect glycolysis, promote malignant biologic processes, and predict the aggressive progression of ccRCC. HK3 may stimulate the abundance of infiltrating monocytes/macrophages presenting surface markers and regulate the key molecular subgroups of immune checkpoint molecules of exhaustive T cells, thus inducing the microenvironmental characteristics of active anti-tumor immune responses.

Glycolysis-Related Genes Serve as Potential Prognostic Biomarkers in Clear Cell Renal Cell Carcinoma

Metabolic rearrangement is a marker of cancer that has been widely studied in recent years. One of the major metabolic characteristics of tumor cells is the high levels of glycolysis, even under aerobic conditions, a phenomenon that is called the "Warburg effect." We investigated the expression and copy number variation (CNV) frequency of all glycolysis-related genes in multiple cancer types and found many differentially expressed genes, particularly in clear cell renal cell carcinoma (ccRCC). Single nucleotide variants (SNVs) showed that the overall average mutation frequency for all genes was low. The purpose of this study was to establish a predictive model by studying glycolysis-related genes in ccRCC. We compared the expression of glycolysis-related genes in 539 ccRCC tissues and 72 normal renal tissues from The Cancer Genome Atlas dataset and identified 17 upregulated and 26 downregulated genes. Pathway analysis revealed that PSAT1 and SDHB could activate the cell cycle, RPIA could activate the DNA damage response, and HK3 could activate apoptosis and EMT signaling, while PDK2 could inhibit apoptosis. The results of the drug sensitivity analysis suggested that some of these differentially expressed genes were positively correlated with drug sensitivity. Thirteen genes were selected from the gene coexpression network and the LASSO regression analysis. The Kaplan-Meier overall survival curves showed that the expression of upregulated genes in ccRCC patients was associated with lower overall survival. We established a predictive model consisting of 13 genes (RPIA, G6PD, PSAT1, ENO2, HK3, IDH1, PDK4, PGM2, PGK1, FBP1, OGDH, SUCLA2, and SUCLG2). This predictive model correlated well with the development and progression of ccRCC. Thus, it is of great value in the diagnosis and prognostic evaluation of ccRCC and may aid the identification of potential prognostic biomarkers and drug targets.

NETO2 Is Deregulated in Breast, Prostate, and Colorectal Cancer and Participates in Cellular Signaling

The NETO2 gene (neuropilin and tolloid-like 2) encodes a protein that acts as an accessory subunit of kainate receptors and is predominantly expressed in the brain. Upregulation of NETO2 has been observed in several tumors; however, its role in tumorigenesis remains unclear. In this study, we investigated NETO2 expression in breast, prostate, and colorectal cancer using quantitative PCR (qPCR), as well as the effect of shRNA-mediated NETO2 silencing on transcriptome changes in colorectal cancer cells. In the investigated tumors, we observed both increased and decreased NETO2 mRNA levels, presenting no correlation with the main clinicopathological characteristics. In HCT116 cells, NETO2 knockdown resulted in the differential expression of 17 genes and 2 long non-coding RNAs (lncRNAs), associated with the upregulation of circadian rhythm and downregulation of several cancer-associated pathways, including Wnt, transforming growth factor (TGF)-β, Janus kinase (JAK)-signal transducer and activator of transcription (STAT), mitogen-activated protein kinase (MAPK), and phosphatidylinositol 3-kinase (PI3K)/protein kinase B (AKT) pathways. Furthermore, we demonstrated the possibility to utilize a novel model organism, short-lived fish Nothobranchius furzeri, for evaluating NETO2 functions. The ortholog neto2b in N. furzeri demonstrated a high similarity in nucleotide and amino acid sequences with human NETO2, as well as was characterized by stable expression in various fish tissues. Collectively, our findings demonstrate the deregulation of NETO2 in the breast, prostate, and colorectal cancer and its participation in the tumor development primarily through cellular signaling.

Hypoxia-inducible hexokinase-2 enhances anti-apoptotic function via activating autophagy in multiple myeloma

Multiple myeloma (MM) is an incurable hematopoietic neoplasm derived from plasma cells, and existing in the bone marrow. Recent developments in the field of myeloma onco-biology have enabled the use of proteasome inhibitors (PIs) as key drugs for MM. PIs can increase cell sensitivity to endoplasmic reticulum stress, leading to apoptosis of myeloma cells. PI cannot kill all myeloma cells, however; one reason of this might be activation of autophagy via hypoxic stress in the bone marrow microenvironment. Hypoxia-inducible gene(s) that regulate autophagy may be novel therapeutic target(s) for PI-resistant myeloma cells. Here, a hypoxia-inducible glycolytic enzyme hexokinase-2 (HK2) was demonstrated to contribute by autophagy activation to the acquisition of an anti-apoptotic phenotype in myeloma cells. We found that hypoxic stress led to autophagy activation accompanied by HK2 upregulation in myeloma cells. Under hypoxic conditions, HK2 knockdown inhibited glycolysis and impaired autophagy, inducing apoptosis. The cooperative effects of a PI (bortezomib) against immunodeficient mice inoculated with HK2-knocked down myeloma cells were examined and significant tumor reduction was observed. An HK2 inhibitor, 3-bromopyruvate (3-BrPA), also induced apoptosis under hypoxic rather than normoxic conditions. Further examination of the cooperative effects between 3-BrPA and bortezomib on myeloma cells revealed a significant increase in apoptotic myeloma cells. These results strongly suggested that HK2 regulates the activation of autophagy in hypoxic myeloma cells. Cooperative treatment using PI against a dominant fraction, and HK2 inhibitor against a minor fraction, adapted to the bone marrow microenvironment, may lead to deeper remission for refractory MM.

MicroRNA-513a-3p regulates colorectal cancer cell metabolism via targeting hexokinase 2

Disruption of cell metabolism is a hallmark of cancer cells. Accumulating evidence suggests that microRNAs (miRNAs/miRs) are involved in almost all physiological and pathological processes. The aberrant expression of miRNAs induces metabolic reprogramming in cancer cells and thus, promotes proliferation. In the current study, miR-513a-3p was identified as a significantly downregulated miRNA in colorectal cancer cells and tumors. Overexpression of miR-513a-3p in colorectal cancer cells inhibited proliferation and glycolysis. A well-documented metabolic regulator, hexokinase 2 (HK2), was predicted and validated HK2to be a target gene of miR-513a-3p in colorectal cancer cells. In addition, overexpression of HK2 reversed the miR-513a-3p mimic-induced inhibition of proliferation. The association between HK2 and miR-513a-3p was further observed in tumors collected from patients with colorectal cancer. The findings suggest that miR-513a-5p may inhibit glycolysis in colorectal cancer cells via repressing HK2 expression, indicating that miR-513a-5p may be a tumor suppressor in colorectal cancer.

Circular RNA circPITX1 knockdown inhibits glycolysis to enhance radiosensitivity of glioma cells by miR-329-3p/NEK2 axis

Background

Numerous circular RNAs (circRNAs) have been recognized as vital modulators of human malignancies, including glioma. Whereas, the functional role of circRNA Pituitary Homeo Box 1 (circPITX1) in the radioresistance of glioma cells remains largely uncertain.

Methods

Quantitative real-time PCR (qRT-PCR) or western blot analysis was employed to examine the expression of circPITX1, microRNA (miR)-329-3p and NIMA-related kinase 2 (NEK2). 3-(4, 5-dimethylthiazol-2-y1)-2, 5-diphenyl tetrazolium bromide (MTT) assay was used to determine cell viability. Glycolysis was assessed by commercial kits and western blot analysis. Colony formation assay was conducted to analyze cell survival and clonogenicity capacity. The relationship among circPITX1, miR-329-3p and NEK2 was confirmed via dual-luciferase reporter assay. The in vivo function of circPITX1 was evaluated by tumor xenograft assay.

Results

Expression of circPITX1 and NEK2 was up-regulated in glioma tissues and cells, while miR-329-3p exhibited reverse trend. CircPITX1 knockdown repressed viability, glycolysis and colony formation, but promoted radiosensitivity of glioma cells, as well as inhibited tumor growth in vivo. MiR-329-3p was a target miRNA of circPITX1 and miR-329-3p deficiency reversed knockdown of circPITX1-mediated glycolysis inhibition and radioresistance reduction. MiR-329-3p exerted inhibitory effects on glycolysis and radioresistance of glioma cells by targeting NEK2. CircPITX1 facilitated NEK2 expression by sponging miR-329-3p. Glycolytic inhibitor 2-deoxy-d-glucose (2-DG) disposition weakened the promoted impact on glycolysis caused by circPITX1.

Conclusion

CircPITX1 knockdown reduced glycolysis to contribute to radiosensitivity in glioma through miR-329-3p/NEK2 axis, providing a possible mechanism of circPITX1 in the development of glioma.

Bioinformatic identification of differentially expressed genes associated with prognosis of locally advanced lymph node-positive prostate cancer

Prostate cancer (PCa) is one of the primary causes of cancer-related mortality in men worldwide. Patients with locally advanced PCa with metastases in regional lymph nodes are usually marked as a high-risk group. One of the chief concerns for this group is to make an informed decision about the necessity of conducting adjuvant androgen deprivation therapy after radical surgical treatment. During the oncogenic transformation and progression of the disease, the expression of many genes is altered. Some of these genes can serve as markers for diagnosis, predicting the prognosis or effectiveness of drug therapy, as well as possible therapeutic targets. We undertook bioinformatic analysis of the RNA-seq data deposited in The Cancer Genome Atlas consortium database to identify possible prognostic markers. We compared the groups with favorable and unfavorable prognosis for the cohort of patients with PCa showing lymph node metastasis (pT2N1M0, pT3N1M0, and pT4N1M0) and for the most common molecular type carrying the fusion transcript TMPRSS2-ERG. For the entire cohort, we revealed at least six potential markers (IDO1, UGT2B15, IFNG, MUC6, CXCL11, and GBP1). Most of these genes are involved in the positive regulation of immune response. For the TMPRSS2-ERG subtype, we also identified six genes, the expression of which may be associated with prognosis: TOB1, GALNT7, INAFM1, APELA, RAC3, and NNMT. The identified genes, after additional studies and validation in the extended cohort, could serve as a prognostic marker of locally advanced lymph node-positive PCa.

Characterization of glycolysis-related gene expression in malignant melanoma

Malignant melanoma exhibits a distinct metabolic phenotype with high glycolytic activity. Previously, we have shown that glucose transporter isoform 1 (GLUT1) favors growth and metastasis of malignant melanoma. In this study, we investigated the expression of GLUT1 and the further glycolysis-related genes hexokinase 1 and 2 (HK1, HK2), lactate dehydrogenase A (LDH-A) and monocarboxylate transporters 1 and 4 (MCT1, MCT4) in eleven human melanoma cell lines under normoxic and hypoxic conditions. Furthermore, a set of 25 human malignant melanoma tissue samples was analyzed. Under hypoxic conditions, we could observe a significant upregulation of hypoxia-inducible factor 1 alpha (HIF-1a) target genes GLUT1, HK2 and LDH-A, but not MCT4. While under normoxic conditions the expression of glycolysis-related genes showed no correlation with origin or BRAF mutation status, GLUT1 expression was significantly elevated in metastatic and BRAF-V600E mutated melanoma cell lines under hypoxic conditions. Furthermore, GLUT1 expression in human melanoma tissue samples correlated significantly with HK1, LDH-A and MCT1 expression, confirming a glycolytic phenotype. Notably, Cyclin D1 expression, which is used as a prognostic marker for the outcome of melanoma patients, as it is associated with proliferation and invasiveness of melanoma, significantly correlated with GLUT1, HK1, LDH-A and MCT1 expression. In summary, our findings provide further evidence that enhanced glycolytic activity in melanoma favors disease progression and is an attractive therapeutic target for this highly aggressive tumor.

Retracted: MicroRNA-34a inhibit hepatocellular carcinoma progression by repressing hexokinase-1

Hepatocellular carcinoma (HCC) is known as a frequent type of primary cancer in the liver, and it is the third-most common cause of cancer-related death all over the world. However, the molecular mechanism in the progression of HCC is still unclear. The current study was designed to investigate the expression and function of microRNA-34a (miR-34a) in HCC. In HCC tissues and cells, the expression levels of miR-34a were analyzed by quantitative real-time polymerase chain reaction. The association between the level of miR-34a and hexokinase (HK)-1 was also investigated via luciferase reporter assay. Cell viability and proliferation were detected by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay and flow cytometry. To assess whether miR-34a can limit tumor growth in vivo, animal models and terminal deoxynucleotidyl transferase dUTP nick end labeling assay were used for examining the role of miR-34a on the development of HCC and cell apoptosis. The expression level of miR-34a was reduced in HCC samples and cells. The expression of miR-34a was associated with the viability and proliferation capacity of HCC cells, and miR-34a could inhibit HCC cells proliferation by inhibiting HK1. In the mouse model of HCC, volumes and weight of the tumors were significantly decreased by transfection with miR-34a mimic compared with the control group. Furthermore, miR-34a mimics could induce apoptosis in a greater proportion of cells compared with the control group. Taken together, the data may provide some novel insights into the molecular mechanism of miR-34a and HK1 in the progression of HCC. Thus, miR-34a/HK1 axis might be a novel promising therapeutic target for treating HCC.

The CIMP-high phenotype is associated with energy metabolism alterations in colon adenocarcinoma

Background

CpG island methylator phenotype (CIMP) is found in 15–20% of malignant colorectal tumors and is characterized by strong CpG hypermethylation over the genome. The molecular mechanisms of this phenomenon are not still fully understood. The development of CIMP is followed by global gene expression alterations and metabolic changes. In particular, CIMP-low colon adenocarcinoma (COAD), predominantly corresponded to consensus molecular subtype 3 (CMS3, “Metabolic”) subgroup according to COAD molecular classification, is associated with elevated expression of genes participating in metabolic pathways.

Methods

We performed bioinformatics analysis of RNA-Seq data from The Cancer Genome Atlas (TCGA) project for CIMP-high and non-CIMP COAD samples with DESeq2, clusterProfiler, and topGO R packages. Obtained results were validated on a set of fourteen COAD samples with matched morphologically normal tissues using quantitative PCR (qPCR).

Results

Upregulation of multiple genes involved in glycolysis and related processes (ENO2, PFKP, HK3, PKM, ENO1, HK2, PGAM1, GAPDH, ALDOA, GPI, TPI1, and HK1) was revealed in CIMP-high tumors compared to non-CIMP ones. Most remarkably, the expression of the PKLR gene, encoding for pyruvate kinase participating in gluconeogenesis, was decreased approximately 20-fold. Up to 8-fold decrease in the expression of OGDHL gene involved in tricarboxylic acid (TCA) cycle was observed in CIMP-high tumors. Using qPCR, we confirmed the increase (4-fold) in the ENO2 expression and decrease (2-fold) in the OGDHL mRNA level on a set of COAD samples.

Conclusions

We demonstrated the association between CIMP-high status and the energy metabolism changes at the transcriptomic level in colorectal adenocarcinoma against the background of immune pathway activation. Differential methylation of at least nine CpG sites in OGDHL promoter region as well as decreased OGDHL mRNA level can potentially serve as an additional biomarker of the CIMP-high status in COAD.

Electronic supplementary material

The online version of this article (10.1186/s12881-019-0771-5) contains supplementary material, which is available to authorized users.

Mutational load in carotid body tumor

Background

Carotid body tumor (CBT) is a rare neoplasm arising from paraganglion located near the bifurcation of the carotid artery. There is great intra-tumor heterogeneity, and CBT development could be associated with both germline and somatic allelic variants. Studies on the molecular genetics of CBT are limited, and the molecular mechanisms of its pathogenesis are not fully understood. This work is focused on the estimation of mutational load (ML) in CBT.

Methods

Using the NextSeq 500 platform, we performed exome sequencing of tumors with matched lymph node tissues and peripheral blood obtained from six patients with CBT. To obtain reliable results in tumors with low ML, we developed and successfully applied a complex approach for the analysis of sequencing data. ML was evaluated as the number of somatic variants per megabase (Mb) of the target regions covered by the Illumina TruSeq Exome Library Prep Kit.

Results

The ML in CBT varied in the range of 0.09–0.28/Mb. Additionally, we identified several pathogenic/likely pathogenic somatic and germline allelic variants across six patients studied (including TP53 variants).

Conclusions

Using the developed approach, we estimated the ML in CBT, which is much lower than in common malignant tumors. Identified variants in known paraganglioma/pheochromocytoma-causative genes and novel genes could be associated with the pathogenesis of CBT. The obtained results expand our knowledge of the mutation process in CBT as well as the biology of tumor development.

Nfe2l1-silenced insulinoma cells acquire aggressiveness and chemoresistance

The transcription factor nuclear factor erythroid 2-like 1 (NFE2L1 or NRF1) is involved in various critical cell processes such as maintenance of ubiquitin-proteasome system and regulation of the cellular antioxidant response. We previously determined that pancreatic β-cell-specific Nfe2l1-knockout mice had hyperinsulinemia and that silencing of Nfe2l1 in mouse islets or MIN6 insulinoma β-cells induced elevated basal insulin release and altered glucose metabolism. Hypoglycemia is a major issue with aggressive insulinomas, although a role of NFE2L1 in this pathology is not defined. In the present work, we studied the tumorigenicity of Nfe2l1-deficient insulinoma MIN6 cells (Nfe2l1-KD) and sensitivity to chemotherapy. Nfe2l1-KD cells grew faster and were more aggressive than Scramble cells in vitro In a mouse allograft transplantation model, insulinomas arising from Nfe2l1-KD cells were more aggressive and chemoresistant. The conclusion was amplified using streptozotocin (STZ) administration in an allograft transplantation model in diabetic Akita background mice. Furthermore, Nfe2l1-KD cells were resistant to damage by the chemotherapeutic drugs STZ and 5-fluorouracil, which was linked to binding of hexokinase 1 with mitochondria, enhanced mitochondrial membrane potential and closed mitochondrial potential transition pore. Overall, both in vitro and in vivo data from Nfe2l1-KD insulinoma cells provided evidence of a previously un-appreciated action of NFE2L1 in suppression of tumorigenesis. Nfe2l1 silencing desensitizes insulinoma cells and derived tumors to chemotherapeutic-induced damage, likely via metabolic reprograming. These data indicate that NFE2L1 could potentially play an important role in the carcinogenic process and impact chemosensitivity, at least within a subset of pancreatic endocrine tumors.

HK3 overexpression associated with epithelial-mesenchymal transition in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is a common cancer worldwide. The main cause of death in CRC includes tumor progression and metastasis. At molecular level, these processes may be triggered by epithelial-mesenchymal transition (EMT) and necessitates specific alterations in cell metabolism. Although several EMT-related metabolic changes have been described in CRC, the mechanism is still poorly understood.

RESULTS

Using CrossHub software, we analyzed RNA-Seq expression profile data of CRC derived from The Cancer Genome Atlas (TCGA) project. Correlation analysis between the change in the expression of genes involved in glycolysis and EMT was performed. We obtained the set of genes with significant correlation coefficients, which included 21 EMT-related genes and a single glycolytic gene, HK3. The mRNA level of these genes was measured in 78 paired colorectal cancer samples by quantitative polymerase chain reaction (qPCR). Upregulation of HK3 and deregulation of 11 genes (COL1A1, TWIST1, NFATC1, GLIPR2, SFPR1, FLNA, GREM1, SFRP2, ZEB2, SPP1, and RARRES1) involved in EMT were found. The results of correlation study showed that the expression of HK3 demonstrated a strong correlation with 7 of the 21 examined genes (ZEB2, GREM1, TGFB3, TGFB1, SNAI2, TWIST1, and COL1A1) in CRC.

CONCLUSIONS

Upregulation of HK3 is associated with EMT in CRC and may be a crucial metabolic adaptation for rapid proliferation, survival, and metastases of CRC cells.

Upregulation of NETO2 gene in colorectal cancer

BACKGROUND

Neuropilin and tolloid-like 2 (NETO2) is a single-pass transmembrane protein that has been shown primarily implicated in neuron-specific processes. Upregulation of NETO2 gene was also detected in several cancer types. In colorectal cancer (CRC), it was associated with tumor progression, invasion, and metastasis, and seems to be involved in epithelial-mesenchymal transition (EMT). However, the mechanism of NETO2 action is still poorly understood.

RESULTS

We have revealed significant increase in the expression of NETO2 gene and deregulation of eight EMT-related genes in CRC. Four of them were upregulated (TWIST1, SNAIL1, LEF1, and FOXA2); the mRNA levels of other genes (FOXA1, BMP2, BMP5, and SMAD7) were decreased. Expression of NETO2 gene was weakly correlated with that of genes involved in the EMT process.

CONCLUSIONS

We found considerable NETO2 upregulation, but no significant correlation between the expression of NETO2 and EMT-related genes in CRC. Thus, NETO2 may be involved in CRC progression, but is not directly associated with EMT.

The role of PPAR alpha in perfluorooctanoic acid induced developmental cardiotoxicity and l-carnitine mediated protection-Results of in ovo gene silencing

Perfluorooctanoic acid (PFOA) is a persistent organic pollutant. This study established an in ovo peroxisome proliferator-activated receptor alpha (PPAR alpha) silencing model in chicken embryo heart, and investigated the role of PPAR alpha in PFOA induced developmental cardiotoxicity. The in ovo silencing was achieved by introducing lentivirus expressing PPAR alpha siRNA into ED2 chicken embryo via microinjection (0.05ul/g egg weight). Transfection efficacy was confirmed by fluorescent microscopy and western blotting. To assess the developmental cardiotoxicity, cardiac function (heart rate) and morphology (right ventricular wall thickness) were measured in D1 hatchling chickens. 2mg/kg (egg weight) PFOA exposure at ED0 induced significant elevation of heart rate and thinning of right ventricular wall thickness in D1 hatchling chickens. PPAR alpha silencing did not prevent PFOA-induced elevation of heart rate; however, it did significantly increase the right ventricular wall thickness as compared to PFOA exposed animals. Meanwhile, PPAR alpha silencing did not abolish the protective effects exerted by exposure to 100mg/kg (egg weight) l-carnitine. In conclusion, PFOA-induced heart rate elevation is likely PPAR alpha independent, while the right ventricular wall thinning seems to be PPAR alpha dependent. The protective effects of l-carnitine do not require PPAR alpha.

miR-142-5p promotes development of colorectal cancer through targeting SDHB and facilitating generation of aerobic glycolysis

Aberrant expression of miRNAs contributes to the development of human malignancies. A recent study revealed that miR-142-5p is increased in the serum of colorectal cancer (CRC) patients compared to health people. Using starBase v2.0, we found that succinate dehydrogenase-B (SDHB) is a potential target of miR-142-5p, while SDHB is negatively correlated to cancer development through regulating energetic metabolism. Based on these information, this study further examined the expression profiles of miR-142-5p and SDHB in CRC tissues and cell lines using PCR and Western blotting. Transfection experiment and luciferase assay were performed to identify relationship between miR-142-5p and SDHB. Oxygen intake, glucose consumption and production of lactic acid were used to evaluate the influence on energetic metabolism. CRC growth and proliferation were assessed by in vitro and in vivo studies. Results showed that miR-142-5p was up-regulated in CRC, but SDHB was down-regulated. SDHB was confirmed as a target of miR-142-5p, and decreased SDHB in CRC was result from the abnormal up-regulation of miR-142-5p. Lose of SDHB by miR-142-5p inhibited oxygen intake by CRC cells, but increased glucose consumption and lactate production. These suggest miR-142-5p up-regulation in CRC probably facilitates generation of aerobic glycolysis by reducing SDHB. miR-142-5p promoted proliferation and colony formation of CRC, but inhibited apoptosis. SDHB overexpression abrogated these effect of miR-142-5p, which indicates that SDHB depletion mediates tumor-promoting actions of miR-142-5p. This study added novel insight into the CRC development regulated by miR-142-5p. It may be a promising therapy target in the future molecular therapy.