The transcription factor Foxp1 regulates aerobic glycolysis in adipocytes and myocytes.

In recent years, lactate has been recognized as an important circulating energy substrate rather than only a dead-end metabolic waste product generated during glucose oxidation at low levels of oxygen. The term "aerobic glycolysis" has been coined to denote increased glucose uptake and lactate production despite normal oxygen levels and functional mitochondria. Hence, in "aerobic glycolysis," lactate production is a metabolic choice, whereas in "anaerobic glycolysis," it is a metabolic necessity based on inadequate levels of oxygen. Interestingly, lactate can be taken up by cells and oxidized to pyruvate and thus constitutes a source of pyruvate that is independent of insulin. Here, we show that the transcription factor Foxp1 regulates glucose uptake and lactate production in adipocytes and myocytes. Overexpression of Foxp1 leads to increased glucose uptake and lactate production. In addition, protein levels of several enzymes in the glycolytic pathway are upregulated, such as hexokinase 2, phosphofructokinase, aldolase, and lactate dehydrogenase. Using chromatin immunoprecipitation and real-time quantitative PCR assays, we demonstrate that Foxp1 directly interacts with promoter consensus cis-elements that regulate expression of several of these target genes. Conversely, knockdown of Foxp1 suppresses these enzyme levels and lowers glucose uptake and lactate production. Moreover, mice with a targeted deletion of Foxp1 in muscle display systemic glucose intolerance with decreased muscle glucose uptake. In primary human adipocytes with induced expression of Foxp1, we find increased glycolysis and glycolytic capacity. Our results indicate Foxp1 may play an important role as a regulator of aerobic glycolysis in adipose tissue and muscle.

Long non-coding RNAs: definitions, functions, challenges and recommendations.

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

Long Non-Coding RNAs: Tools for Understanding and Targeting Cancer Pathways.

The regulatory nature of long non-coding RNAs (lncRNAs) has been well established in various processes of cellular growth, development, and differentiation. Therefore, it is vital to examine their contribution to cancer development. There are ample examples of lncRNAs whose cellular levels are significantly associated with clinical outcomes. However, whether these non-coding molecules can work as either key drivers or barriers to cancer development remains unknown. The current review aims to discuss some well-characterised lncRNAs in the process of oncogenesis and extrapolate the extent of their decisive contribution to tumour development. We ask if these lncRNAs can independently initiate neoplastic lesions or they always need the modulation of well characterized oncogenes or tumour suppressors to exert their functional properties. Finally, we discuss the emerging genetic approaches and appropriate animal and humanised models that can significantly contribute to the functional dissection of lncRNAs in cancer development and progression.

HnRNPK maintains single strand RNA through controlling double-strand RNA in mammalian cells.

Although antisense transcription is a widespread event in the mammalian genome, double-stranded RNA (dsRNA) formation between sense and antisense transcripts is very rare and mechanisms that control dsRNA remain unknown. By characterizing the FGF-2 regulated transcriptome in normal and cancer cells, we identified sense and antisense transcripts IER3 and IER3-AS1 that play a critical role in FGF-2 controlled oncogenic pathways. We show that IER3 and IER3-AS1 regulate each other's transcription through HnRNPK-mediated post-transcriptional regulation. HnRNPK controls the mRNA stability and colocalization of IER3 and IER3-AS1. HnRNPK interaction with IER3 and IER3-AS1 determines their oncogenic functions by maintaining them in a single-stranded form. hnRNPK depletion neutralizes their oncogenic functions through promoting dsRNA formation and cytoplasmic accumulation. Intriguingly, hnRNPK loss-of-function and gain-of-function experiments reveal its role in maintaining global single- and double-stranded RNA. Thus, our data unveil the critical role of HnRNPK in maintaining single-stranded RNAs and their physiological functions by blocking RNA-RNA interactions.

Filamin A increases aggressiveness of human neuroblastoma.

Background: The actin-binding protein filamin A (FLNA) regulates oncogenic signal transduction important for tumor growth, but the role of FLNA in the progression of neuroblastoma (NB) has not been explored. Methods: We analyzed FLNA mRNA expression in the R2 NB-database and FLNA protein expression in human NB tumors. We then silenced FLNA expression in human SKNBE2 and IMR32 NB cells by lentiviral vector encoding shRNA FLNA and assayed the cells for proliferation, migration, colony, spheroid formation, and apoptosis. SKNBE2 xenografts expressing or lacking FLNA in BALB/c nude mice were analyzed by both routine histopathology and immunohistochemistry. Results: We observed shorter patient survival with higher expression of FLNA mRNA than patients with lower FLNA mRNA expression, and high-risk NB tumors expressed higher FLNA levels. Overexpression of FLNA increased proliferation of SH-SY5 NB cells. NB cell lines transfected with siRNA FLNA proliferated and migrated less, expressed lower levels of phosphorylated AKT and ERK1/2, formed smaller colonies and spheroids, as well as increased apoptosis. After inoculation of SKNBE2 cells infected with lentivirus expressing shRNA FLNA, size of NB tumors and number of proliferating cells were decreased. Furthermore, we identified STAT3 as an interacting partner of FLNA. Silencing FLNA mRNA reduced levels of NF-κB, STAT3 and MYCN, and increased levels of p53 and cleaved caspase 3. Conclusion: Inhibition of FLNA impaired NB cell signaling and function and reduced NB tumor size in vivo, suggesting that drugs targeting either FLNA or its interaction with STAT3 may be useful in the treatment of NB.

The DNA damage inducible lncRNA SCAT7 regulates genomic integrity and topoisomerase 1 turnover in lung adenocarcinoma.

Despite the rapid improvements in unveiling the importance of lncRNAs in all aspects of cancer biology, there is still a void in mechanistic understanding of their role in the DNA damage response. Here we explored the potential role of the oncogenic lncRNA SCAT7 (ELF3-AS1) in the maintenance of genome integrity. We show that SCAT7 is upregulated in response to DNA-damaging drugs like cisplatin and camptothecin, where SCAT7 expression is required to promote cell survival. SCAT7 silencing leads to decreased proliferation of cisplatin-resistant cells in vitro and in vivo through interfering with cell cycle checkpoints and DNA repair molecular pathways. SCAT7 regulates ATR signaling, promoting homologous recombination. Importantly, SCAT7 also takes part in proteasome-mediated topoisomerase I (TOP1) degradation, and its depletion causes an accumulation of TOP1-cc structures responsible for the high levels of intrinsic DNA damage. Thus, our data demonstrate that SCAT7 is an important constituent of the DNA damage response pathway and serves as a potential therapeutic target for hard-to-treat drug resistant cancers.

NBAT1/CASC15-003/USP36 control MYCN expression and its downstream pathway genes in neuroblastoma.

BACKGROUND: MYCN has been an attractive therapeutic target in neuroblastoma considering the widespread amplification of the MYCN locus in neuroblastoma, and its established role in neuroblastoma development and progression. Thus, understanding neuroblastoma-specific control of MYCN expression at the transcriptional and post-transcriptional level would lead to identification of novel MYCN-dependent oncogenic pathways and potential therapeutic strategies. METHODS: By performing loss- and gain-of-function experiments of the neuroblastoma hotspot locus 6p22.3 derived lncRNAs CASC15-003 and NBAT1, together with coimmunoprecipitation and immunoblotting of MYCN, we have shown that both lncRNAs post-translationally control the expression of MYCN through regulating a deubiquitinase enzyme USP36. USP36 oncogenic properties were investigated using cancer cell lines and in vivo models. RNA-seq analysis of loss-of-function experiments of CASC15-003/NBAT1/MYCN/USP36 and JQ1-treated neuroblastoma cells uncovered MYCN-dependent oncogenic pathways. RESULTS: We show that NBAT1/CASC15-003 control the stability of MYCN protein through their common interacting protein partner USP36. USP36 harbors oncogenic properties and its higher expression in neuroblastoma patients correlates with poor prognosis, and its downregulation significantly reduces tumor growth in neuroblastoma cell lines and xenograft models. Unbiased integration of RNA-seq data from CASC15-003, NBAT1, USP36, and MYCN knockdowns and neuroblastoma cells treated with MYCN inhibitor JQ1, identified genes that are jointly regulated by the NBAT1/CASC15-003/USP36/MYCN pathway. Functional experiments on one of the target genes, COL18A1, revealed its role in the NBAT1/CASC15-003-dependent cell adhesion feature in neuroblastoma cells. CONCLUSION: Our data show post-translational regulation of MYCN by NBAT1/CASC15-003/USP36, which represents a new regulatory layer in the complex multilayered gene regulatory network that controls MYCN expression.

LY6K-AS lncRNA is a lung adenocarcinoma prognostic biomarker and regulator of mitotic progression.

Recent advances in genomics unraveled several actionable mutational drivers in lung cancer, leading to promising therapies such as tyrosine kinase inhibitors and immune checkpoint inhibitors. However, the tumors' acquired resistance to the newly-developed as well as existing therapies restricts life quality improvements. Therefore, we investigated the noncoding portion of the human transcriptome in search of alternative actionable targets. We identified an antisense transcript, LY6K-AS, with elevated expression in lung adenocarcinoma (LUAD) patients, and its higher expression in LUAD patients predicts poor survival outcomes. LY6K-AS abrogation interfered with the mitotic progression of lung cancer cells resulting in unfaithful chromosomal segregation. LY6K-AS interacts with and stabilizes 14-3-3 proteins to regulate the transcription of kinetochore and mitotic checkpoint proteins. We also show that LY6K-AS regulates the levels of histone H3 lysine 4 trimethylation (H3K4me3) at the promoters of kinetochore members. Cisplatin treatment and LY6K-AS silencing affect many common pathways enriched in cell cycle-related functions. LY6K-AS silencing affects the growth of xenografts derived from wildtype and cisplatin-resistant lung cancer cells. Collectively, these data indicate that LY6K-AS silencing is a promising therapeutic option for LUAD that inhibits oncogenic mitotic progression.

The long noncoding RNA TUNAR modulates Wnt signaling and regulates human ß-cell proliferation.

Many long noncoding RNAs (lncRNAs) are enriched in pancreatic islets and several lncRNAs are linked to type 2 diabetes (T2D). Although they have emerged as potential players in ß-cell biology and T2D, little is known about their functions and mechanisms in human ß-cells. We identified an islet-enriched lncRNA, TUNAR (TCL1 upstream neural differentiation-associated RNA), which was upregulated in ß-cells of patients with T2D and promoted human ß-cell proliferation via fine-tuning of the Wnt pathway. TUNAR was upregulated following Wnt agonism by a glycogen synthase kinase-3 (GSK3) inhibitor in human ß-cells. Reciprocally, TUNAR repressed a Wnt antagonist Dickkopf-related protein 3 (DKK3) and stimulated Wnt pathway signaling. DKK3 was aberrantly expressed in ß-cells of patients with T2D and displayed a synchronized regulatory pattern with TUNAR at the single cell level. Mechanistically, DKK3 expression was suppressed by the repressive histone modifier enhancer of zeste homolog 2 (EZH2). TUNAR interacted with EZH2 in ß-cells and facilitated EZH2-mediated suppression of DKK3. These findings reveal a novel cell-specific epigenetic mechanism via islet-enriched lncRNA that fine-tunes the Wnt pathway and subsequently human ß-cell proliferation.NEW & NOTEWORTHY The discovery that long noncoding RNA TUNAR regulates ß-cell proliferation may be important in designing new treatments for diabetes.

Labeling and Purification of Temporally Expressed RNAs During the S-Phase of the Cell Cycle in Living Cells.

From high-throughput DNA and RNA sequencing technologies, it is evident that more than two-thirds of the mammalian genome is transcribed and nearly 98% of the transcriptional output in humans constitute noncoding RNA, comprising tens of thousands of small and long noncoding RNAs. These observations have put the study of RNA expression levels at the center of molecular biology research. The transcriptional output of cells changes temporally throughout different cell cycle phases, or in response to a large panel of stimuli. In such instances, the measure of induced RNA transcripts might be obscured by the presence of steady-state RNA levels in the total transcriptome. With this protocol, we provide a method for labeling and purification of the nascent RNAs transcribed over short periods of time in cultured cells. The supplementation of cell culture medium with a chemically modified analog of uridine, ethynyl-uridine, allows for the subsequent biotinylation of ethynyl-uridine residues with a click-chemistry reaction. The labeled RNA is then purified on streptavidin beads and eluted. The purified RNA is suitable for use in RT-qPCR assays as well as in deep sequencing applications.

In Vivo Administration of Therapeutic Antisense Oligonucleotides.

With the rapid revolution in RNA/DNA sequencing technologies, it is evident that mammalian genomes express tens of thousands of long noncoding RNAs (lncRNAs). Since a large majority of lncRNAs have been functionally implicated in cancer development and progression, there is an increasing appreciation for the use of antisense oligonucleotide (ASO)-based therapies targeting lncRNAs in several cancers. Despite their great potential in therapeutic applications, their use is still limited due to cellular toxicity and shortcomings in achieving required stability in biological fluids and tissue uptake. To overcome these limitations, major changes in ASO chemistry have been introduced to generate second and third generation ASOs, including locked nucleic acids (LNA) technology. Here we describe two different LNA-ASO delivery approaches, a peritumoral administration and a systemic delivery in xenograft models of lung adenocarcinoma, that significantly reduced tumor growth without inducing toxicity.

Subcellular Distribution of p53 by the p53-Responsive lncRNA NBAT1 Determines Chemotherapeutic Response in Neuroblastoma.

Neuroblastoma has a low mutation rate for the p53 gene. Alternative ways of p53 inactivation have been proposed in neuroblastoma, such as abnormal cytoplasmic accumulation of wild-type p53. However, mechanisms leading to p53 inactivation via cytoplasmic accumulation are not well investigated. Here we show that the neuroblastoma risk-associated locus 6p22.3-derived tumor suppressor NBAT1 is a p53-responsive lncRNA that regulates p53 subcellular levels. Low expression of NBAT1 provided resistance to genotoxic drugs by promoting p53 accumulation in cytoplasm and loss from mitochondrial and nuclear compartments. Depletion of NBAT1 altered CRM1 function and contributed to the loss of p53-dependent nuclear gene expression during genotoxic drug treatment. CRM1 inhibition rescued p53-dependent nuclear functions and sensitized NBAT1-depleted cells to genotoxic drugs. Combined inhibition of CRM1 and MDM2 was even more effective in sensitizing aggressive neuroblastoma cells with p53 cytoplasmic accumulation. Thus, our mechanistic studies uncover an NBAT1-dependent CRM1/MDM2-based potential combination therapy for patients with high-risk neuroblastoma. SIGNIFICANCE: This study shows how a p53-responsive lncRNA mediates chemotherapeutic response by modulating nuclear p53 pathways and identifies a potential treatment strategy for patients with high-risk neuroblastoma.

Analysis of copy number alterations reveals the lncRNA ALAL-1 as a regulator of lung cancer immune evasion.

Cancer is characterized by genomic instability leading to deletion or amplification of oncogenes or tumor suppressors. However, most of the altered regions are devoid of known cancer drivers. Here, we identify lncRNAs frequently lost or amplified in cancer. Among them, we found amplified lncRNA associated with lung cancer-1 (ALAL-1) as frequently amplified in lung adenocarcinomas. ALAL-1 is also overexpressed in additional tumor types, such as lung squamous carcinoma. The RNA product of ALAL-1 is able to promote the proliferation and tumorigenicity of lung cancer cells. ALAL-1 is a TNFα- and NF-κB-induced cytoplasmic lncRNA that specifically interacts with SART3, regulating the subcellular localization of the protein deubiquitinase USP4 and, in turn, its function in the cell. Interestingly, ALAL-1 expression inversely correlates with the immune infiltration of lung squamous tumors, while tumors with ALAL-1 amplification show lower infiltration of several types of immune cells. We have thus unveiled a pro-oncogenic lncRNA that mediates cancer immune evasion, pointing to a new target for immune potentiation.

Sperm Originated Chromatin Imprints and LincRNAs in Organismal Development and Cancer.

Importance of sperm-derived transcripts and chromatin imprints in organismal development is poorly investigated. Here using an integrative approach, we show that human sperm transcripts are equally important as oocyte. Sperm-specific and sperm-oocyte common transcripts carry distinct chromatin structures at their promoters correlating with corresponding transcript levels in sperm. Interestingly, sperm-specific H3K4me3 patterns at the lincRNA promoters are not maintained in the germ layers and somatic tissues. However, bivalent chromatin at the sperm-specific protein-coding gene promoters is maintained throughout the development. Sperm-specific transcripts reach their peak expression during zygotic genome activation, whereas sperm-oocyte common transcripts are present during early preimplantation development but decline at the onset of zygotic genome activation. Additionally, there is an inverse correlation between sperm-specific and sperm-oocyte lincRNAs throughout the development. Sperm-lincRNAs also show aberrant activation in tumors. Overall, our observations indicate that sperm transcripts carrying chromatin imprints may play an important role in human development and cancer.

Understanding Long Noncoding RNA and Chromatin Interactions: What We Know So Far.

With the evolution of technologies that deal with global detection of RNAs to probing of lncRNA-chromatin interactions and lncRNA-chromatin structure regulation, we have been updated with a comprehensive repertoire of chromatin interacting lncRNAs, their genome-wide chromatin binding regions and mode of action. Evidence from these new technologies emphasize that chromatin targeting of lncRNAs is a prominent mechanism and that these chromatin targeted lncRNAs exert their functionality by fine tuning chromatin architecture resulting in an altered transcriptional readout. Currently, there are no unifying principles that define chromatin association of lncRNAs, however, evidence from a few chromatin-associated lncRNAs show presence of a short common sequence for chromatin targeting. In this article, we review how technological advancements contributed in characterizing chromatin associated lncRNAs, and discuss the potential mechanisms by which chromatin associated lncRNAs execute their functions.

Transcriptome-wide Profiling of Cerebral Cavernous Malformations Patients Reveal Important Long noncoding RNA molecular signatures.

Cerebral cavernous malformations (CCMs) are low-flow vascular malformations in the brain associated with recurrent hemorrhage and seizures. The current treatment of CCMs relies solely on surgical intervention. Henceforth, alternative non-invasive therapies are urgently needed to help prevent subsequent hemorrhagic episodes. Long non-coding RNAs (lncRNAs) belong to the class of non-coding RNAs and are known to regulate gene transcription and involved in chromatin remodeling via various mechanism. Despite accumulating evidence demonstrating the role of lncRNAs in cerebrovascular disorders, their identification in CCMs pathology remains unknown. The objective of the current study was to identify lncRNAs associated with CCMs pathogenesis using patient cohorts having 10 CCM patients and 4 controls from brain. Executing next generation sequencing, we performed whole transcriptome sequencing (RNA-seq) analysis and identified 1,967 lncRNAs and 4,928 protein coding genes (PCGs) to be differentially expressed in CCMs patients. Among these, we selected top 6 differentially expressed lncRNAs each having significant correlative expression with more than 100 differentially expressed PCGs. The differential expression status of the top lncRNAs, SMIM25 and LBX2-AS1 in CCMs was further confirmed by qRT-PCR analysis. Additionally, gene set enrichment analysis of correlated PCGs revealed critical pathways related to vascular signaling and important biological processes relevant to CCMs pathophysiology. Here, by transcriptome-wide approach we demonstrate that lncRNAs are prevalent in CCMs disease and are likely to play critical roles in regulating important signaling pathways involved in the disease progression. We believe, that detailed future investigations on this set of identified lncRNAs can provide useful insights into the biology and, ultimately, contribute in preventing this debilitating disease.

Author Correction: MEG3 long noncoding RNA regulates the TGF-ß pathway genes through formation of RNA-DNA triplex structures.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

BACH1 Stabilization by Antioxidants Stimulates Lung Cancer Metastasis.

For tumors to progress efficiently, cancer cells must overcome barriers of oxidative stress. Although dietary antioxidant supplementation or activation of endogenous antioxidants by NRF2 reduces oxidative stress and promotes early lung tumor progression, little is known about its effect on lung cancer metastasis. Here, we show that long-term supplementation with the antioxidants N-acetylcysteine and vitamin E promotes KRAS-driven lung cancer metastasis. The antioxidants stimulate metastasis by reducing levels of free heme and stabilizing the transcription factor BACH1. BACH1 activates transcription of Hexokinase 2 and Gapdh and increases glucose uptake, glycolysis rates, and lactate secretion, thereby stimulating glycolysis-dependent metastasis of mouse and human lung cancer cells. Targeting BACH1 normalized glycolysis and prevented antioxidant-induced metastasis, while increasing endogenous BACH1 expression stimulated glycolysis and promoted metastasis, also in the absence of antioxidants. We conclude that BACH1 stimulates glycolysis-dependent lung cancer metastasis and that BACH1 is activated under conditions of reduced oxidative stress.

EZH2 upregulates the PI3K/AKT pathway through IGF1R and MYC in clinically aggressive chronic lymphocytic leukaemia.

EZH2 is overexpressed in poor-prognostic chronic lymphocytic leukaemia (CLL) cases, acting as an oncogene; however, thus far, the EZH2 target genes in CLL have not been disclosed. In this study, using ChIP-sequencing, we identified EZH2 and H3K27me3 target genes in two prognostic subgroups of CLL with distinct prognosis and outcome, i.e., cases with unmutated (U-CLL, n = 6) or mutated IGHV genes (M-CLL, n = 6). While the majority of oncogenic pathways were equally enriched for EZH2 target genes in both prognostic subgroups, PI3K pathway genes were differentially bound by EZH2 in U-CLL versus M-CLL. The occupancy of EZH2 for selected PI3K pathway target genes was validated in additional CLL samples (n = 16) and CLL cell lines using siRNA-mediated EZH2 downregulation and ChIP assays. Intriguingly, we found that EZH2 directly binds to the IGF1R promoter along with MYC and upregulates IGF1R expression in U-CLL, leading to downstream PI3K activation. By investigating an independent CLL cohort (n = 96), a positive correlation was observed between EZH2 and IGF1R expression with higher levels in U-CLL compared to M-CLL. Accordingly, siRNA-mediated downregulation of either EZH2, MYC or IGF1R and treatment with EZH2 and MYC pharmacological inhibitors in the HG3 CLL cell line induced a significant reduction in PI3K pathway activation. In conclusion, we characterize for the first time EZH2 target genes in CLL revealing a hitherto unknown implication of EZH2 in modulating the PI3K pathway in a non-canonical, PRC2-independent way, with potential therapeutic implications considering that PI3K inhibitors are effective therapeutic agents for CLL.

FOXK1 and FOXK2 regulate aerobic glycolysis.

Adaptation to the environment and extraction of energy are essential for survival. Some species have found niches and specialized in using a particular source of energy, whereas others-including humans and several other mammals-have developed a high degree of flexibility1. A lot is known about the general metabolic fates of different substrates but we still lack a detailed mechanistic understanding of how cells adapt in their use of basic nutrients2. Here we show that the closely related fasting/starvation-induced forkhead transcription factors FOXK1 and FOXK2 induce aerobic glycolysis by upregulating the enzymatic machinery required for this (for example, hexokinase-2, phosphofructokinase, pyruvate kinase, and lactate dehydrogenase), while at the same time suppressing further oxidation of pyruvate in the mitochondria by increasing the activity of pyruvate dehydrogenase kinases 1 and 4. Together with suppression of the catalytic subunit of pyruvate dehydrogenase phosphatase 1 this leads to increased phosphorylation of the E1α regulatory subunit of the pyruvate dehydrogenase complex, which in turn inhibits further oxidation of pyruvate in the mitochondria-instead, pyruvate is reduced to lactate. Suppression of FOXK1 and FOXK2 induce the opposite phenotype. Both in vitro and in vivo experiments, including studies of primary human cells, show how FOXK1 and/or FOXK2 are likely to act as important regulators that reprogram cellular metabolism to induce aerobic glycolysis.