- Research
- Open access
- Published:
The m6A modification mediated-lncRNA POU6F2-AS1 reprograms fatty acid metabolism and facilitates the growth of colorectal cancer via upregulation of FASN
Molecular Cancer volume 23, Article number: 55 (2024)
Abstract
Background
Long noncoding RNAs (lncRNAs) have emerged as key players in tumorigenesis and tumour progression. However, the biological functions and potential mechanisms of lncRNAs in colorectal cancer (CRC) are unclear.
Methods
The novel lncRNA POU6F2-AS1 was identified through bioinformatics analysis, and its expression in CRC patients was verified via qRT–PCR and FISH. In vitro and in vivo experiments, such as BODIPY staining, Oil Red O staining, triglyceride (TAG) assays, and liquid chromatography mass spectrometry (LC-MS) were subsequently performed with CRC specimens and cells to determine the clinical significance, and functional roles of POU6F2-AS1. Biotinylated RNA pull-down, RIP, Me-RIP, ChIP, and patient-derived organoid (PDO) culture assays were performed to confirm the underlying mechanism of POU6F2-AS1.
Results
The lncRNA POU6F2-AS1 is markedly upregulated in CRC and associated with adverse clinicopathological features and poor overall survival in CRC patients. Functionally, POU6F2-AS1 promotes the growth and lipogenesis of CRC cells both in vitro and in vivo. Mechanistically, METTL3-induced m6A modification is involved in the upregulation of POU6F2-AS1. Furthermore, upregulated POU6F2-AS1 could tether YBX1 to the FASN promoter to induce transcriptional activation, thus facilitating the growth and lipogenesis of CRC cells.
Conclusions
Our data revealed that the upregulation of POU6F2-AS1 plays a critical role in CRC fatty acid metabolism and might provide a novel promising biomarker and therapeutic target for CRC.
Background
Colorectal cancer (CRC) is the most common malignancy of the digestive system, and the second leading cause of cancer-related mortality worldwide [1]. Although significant progress has been made in surgical resection, chemotherapy, targeted therapy and immunotherapy for CRC, the prognosis of CRC patient is still unfavourable [2,3,4]. CRC hallmarks include rapid proliferation, which requires nutrients and facilitates the removal of cellular waste in nutrient-poor environments [5]. Research progress in the last decade has demonstrated that metabolic reprogramming is a pivotal general characteristic of cancer, that allow CRC cells to produce energy sufficient for rapid proliferation [6]. Accumulating evidence has confirmed aberrant metabolic reprogramming as a critical contributor to CRC initiation, proliferation and metastasis [7]. Hence, elucidating the regulatory mechanisms linking colorectal carcinogenesis factors and metabolic enzymes may provide promising avenues for improving CRC prognosis and developing effective therapeutic strategies.
As a critical metabolic feature of cancer cells, lipid metabolism has been widely reported to participate in cancer initiation, progression and drug resistance [8]. Fatty acid synthase (FASN) is a central enzyme of lipid metabolism that can catalyse de novo fatty acid (FA) biosynthesis to promote cell growth and survival [9]. For instance, loss of FBXW7β function facilitates FASN protein stability and promotes lipogenesis and growth in CRC [10]. However, the mechanisms underlying the dysregulation of lipogenesis and aberrant FASN expression in CRC still need to be clarified.
Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nucleotides in length that have limited or no ability to encode proteins, but they can regulate protein-coding gene expression by modulating mRNA transcription, translation, and post- translational modification [11, 12]. Over the past decade, lncRNAs have been reported to play vital roles in various cellular functions in CRC, including proliferation, apoptosis, migration and invasion [13]. Previously, we reported that lncRNAs played important roles in CRC tumorigenesis and progression [14,15,16]. Furthermore, several lncRNAs were reported to be critical players in lipid metabolism modulation by CRC cells. For example, ZFAS1 binds with PABP2 to stabilize SREBP1 mRNA, thus facilitating lipid accumulation in CRC cells [17]. SNHG16 acts as a “sponge” for miRNAs to promote lipid metabolism by upregulating SCD [18]. However, additional crosstalk between lncRNAs and metabolic enzymes or processes still needs to be characterized.
In the present study, we identified a novel lncRNA, POU6F2-AS1, whose expression is significantly upregulated in CRC; this upregulation phenotype is positively associated with poor clinical outcomes in CRC patients. Moreover, we found that POU6F2-AS1 drives the growth and lipogenesis of CRC cells in vitro and in vivo, and these effects partially dependent on YBX1-mediated transcriptional activation of FASN. We also demonstrated that METTL3-induced m6A modification maintains POU6F2-AS1 stability and increases its expression. Taken together, these findings revealed that POU6F2-AS1 may serve as a promising biomarker and therapeutic target for CRC.
Materials and methods
Clinical specimens
This study was reviewed and approved by the Ethics Committee of the Affiliated Hospital of Xuzhou Medical University (XYFY2022-KL473-01), and consent was obtained from all participants. A cohort comprising 84 pairs of CRC tissues and corresponding adjacent normal tissues (ANTs) was collected from CRC patients who underwent radical resection or palliative resection at the Department of Gastrointestinal Surgery of the Affiliated Hospital of Xuzhou Medical University between January 2023 and April 2023 and was used for extracting RNA to detect the expression level of POU6F2-AS1. In addition, a cohort of 60 pairs consisting of paraffin-embedded CRC tissues and corresponding ANTs was collected for tissue microarray (TMA) analysis. All CRC patients had a pathologically confirmed diagnosis and did not receive preoperative radiotherapy or chemotherapy. Clinicopathological characteristics, such as patient age and sex and tumour size and depth of invasion, etc., were also obtained.
Cell lines and culture
The human normal colorectal epithelial cell line FHC was obtained from the American Type Culture Collection (Manassas, VA, USA), while several human CRC cell lines (HCT116, DLD1, LoVo, SW620, HT29 and SW480) were obtained from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). HCT116 cells were cultured in McCoy´s 5 A medium supplemented with 10% foetal bovine serum (FBS) (Gibco, MA, USA), while SW620 cells were cultured in L-15 medium (Gibco, MA, USA). DLD1, HT-29, and LoVo cells were cultured in RPMI 1640 medium (Gibco, MA, USA), and FHC and SW480 cells were cultured in DMEM/high-glucose medium (Gibco, MA, USA). All the cells were cultured in a constant temperature incubator at 37 °C containing 5% CO2 and subjected to mycoplasma testing to confirm that they were free of contamination before use.
BODIPY staining
The lipophilic fluorescence dye BODIPY (Molecular Probes, MA, USA) was used to monitor the content of neutral lipids in CRC cells. The transfected CRC cells were inoculated into 24-well plates containing coverslips and fixed with 4% paraformaldehyde for 30 min after 24 h. The cells were subsequently incubated with 2 µM BODIPY working solution in the dark for 30 min at 37 °C. The staining procedure was stopped by washing the cells with PBS, after which the nuclei were counterstained with Hoechst 33,342 (Beyotime, Shanghai, China) for 10 min. Representative images were taken with a confocal laser scanning microscope (Carl Zeiss, Germany).
Oil Red O staining
Oil Red O staining was performed using Oil Red O Saturated Solution (Solarbio Life Sciences, Beijing, China) according to the manufacturer’s instructions. Briefly, an Oil Red O working solution was prepared by mixing 3 parts of Oil Red O saturated solution with 2 parts deionized water. The transfected CRC cells were inoculated on 14 mm round coverslips in 24-well plates, washed twice with PBS after 24 h, and fixed in 4% paraformaldehyde for 30 min before staining. For frozen tissue sectioning, optimal cutting temperature (OCT)-embedded tissue was cut into 8 μm sections, which were fixed in 4% paraformaldehyde for 30 min before staining. The cells or tissue sections were incubated in 60% isopropanol for 30 s and then in Oil Red O working solution for 20 min at room temperature. Then sections were subsequently washed with 60% isopropanol and deionized water again, counterstained with haematoxylin (Solarbio Life Sciences, Beijing, China), and photographed by an Olympus microscope (Tokyo, Japan). The Oil Red O staining area was calculated as Oil Red O+ area divided by the area of hematoxylin staining and normalized to the average of the control.
Triglyceride (TAG) assay
For the quantitative estimation of triglycerides in cells, a Triglyceride Assay Kit (Promega, WI, USA) was used in accordance with the manufacturer’s protocols. Transfected cells were inoculated into 96-well plates and subjected to a triglyceride assay after 24 h. The cells were washed twice with PBS, and glycerol detection reagent was added, followed by incubation at room temperature for 1 h to detect luminescence. The luminescence signal was measured using a Tecan Spark zymography (Zurich, CH).
Biotinylated RNA pull-down assay
In brief, T7 promoter-containing DNA was obtained by PCR amplification and purified using the TIANgel Purification Kit (TIANGEN BIOTECH, Beijing, China). The purified T7 promoter-containing DNA was incubated with a biotin RNA labelling mixture (Promega Corporation, MA, USA), T7 RNA polymerase (Thermo Fisher Scientific, MA, USA) and an RNase inhibitor for 2 h at a constant temperature of 37 °C in a PCR instrument and purified using an RNAclean Kit (TIANGEN BIOTECH, Beijing, China) to obtain the biotin-labelled RNA probe. Biotinylated RNAs were incubated with protein extracted from CRC cells and mixed with precleared streptavidin agarose resin (Thermo Fisher Scientific, MA, USA). The RNA-protein complexes were obtained by collecting the agarose resin by centrifugation, and the proteins bound therein were eluted and denatured, followed by electrophoresis on SDS-PAGE gels. Afterwards, the proteins were silver stained, and the differential bands were excised for mass spectrometry analysis (Shanghai Applied Protein Technology Co. Ltd., China). To identify which specific region of POU6F2-AS1 interacts with YBX1, we truncated POU6F2-AS1.
Luciferase reporter assay
The luciferase plasmid was constructed by Gene Create Biologicals (Wuhan, China). To investigate the transcriptional regulation of FASN, two different fragments of the FASN promoter region were PCR amplified and subsequently inserted into the KpnI/XhoI site upstream of firefly luciferase in the pGL3-Basic vector. Cells with or without knockdown or overexpression of YBX1/POU6F2-AS1 were transfected with the luciferase plasmid. Firefly luciferase activity was measured 48 h later using a Dual Luciferase Reporter Gene Assay Kit (Beyotime, Shanghai, China) according to the manufacturer’s instructions, with Renilla luciferase serving as a transfection control.
Chromatin immunoprecipitation (ChIP) assay
The ChIP assay was performed using the BeyoChIP™ ChIP Assay Kit (Beyotime, Shanghai, China) according to the manufacturer’s instructions. Briefly, 1 × 107 CRC cells were prepared and then cross-linked with 1% formaldehyde for 15 min at room temperature. The cross-linking process was quenched by adding 0.125 M glycine. Chromatin was isolated with the lysis buffer provided in the kit. Afterwards, sonication was performed to shear the DNA to 200–1000 bp. Immunoprecipitation of cross-linked protein/DNA was performed by incubating anti-YBX1 antibody (Proteintech, IL, USA) or normal rabbit IgG (Proteintech, IL, USA) with sheared chromatin overnight at 4 °C followed by incubation with Protein A/G Magnetic Beads for 2 h. Subsequently, cross-links of protein/DNA were released with 5 M NaCl, and the DNA was purified using a PCR Clean Up Kit (Beyotime, Shanghai, China). Finally, the purified DNA was subjected to qRT–PCR using the P1-P9 primer sequences listed in Additional file 1: Table S4.
Patient-derived organoid (PDO) culture model
In this study, a PDO model was established using tumour tissues from CRC patients at the Department of Gastrointestinal Surgery of the Affiliated Hospital of Xuzhou Medical University. Briefly, freshly removed tumour tissues were cut into 1–3 mm3 pieces, washed several times with antibiotic-containing PBS, and digested with 200 U/ml collagenase (Sigma, MO, USA) and 100 U/ml hyaluronidase (Sigma, MO, USA) for 30 min at 37℃. The digested tissues were filtered through a 100 μm cell filter and then centrifuged to discard the supernatant, after which the cells were resuspended in Matrigel. Then, 50 µL of Matrigel resuspension solution was added to each well of a 48-well plate, and 450 µL of CRC organic medium (bioGenous, Suzhou, China) was added to each well after solidification. The medium was changed every 2–3 days, and the organoids were passaged every 10–15 days. Lentiviral transfection was performed 24 h after organoid passaging. Organoid images were obtained with an Olympus FSX100 microscope (Olympus, Tokyo, Japan).
Animal experiments
Female BALB/c nude mice ranging from 4 to 6 weeks old were obtained from Beijing Vital River Laboratory Animal Technology Co., Ltd. (Beijing, China), and housed under specific pathogen-free conditions. To construct a xenograft tumorigenesis model, transfected HCT116 cells were injected subcutaneously into the axillary region of mice (5 × 106 cells). Tumour growth was measured every 3 days, and the tumour volume (V) was calculated according to the following formula: V= (length × width2)/2. When the maximum tumour length diameter reached 15 mm, the mice were sacrificed, and the subcutaneous tumours were removed to calculate the mass and then subjected to HE and IHC staining. The Animal Care and Use Committee of Xuzhou Medical University provided a statement of ethical approval for all animal experiments in this study.
Statistical analysis
Statistical analyses were performed using SPSS 19.0 software (IBM, Armonk, NY, USA) and GraphPad Prism version 8.0 (La Jolla, CA, USA). All the data in the present study is presented as the means ± standard deviations (SD), and all the tests were two-sided. A P value < 0.05 was considered wo indicate statistical significance. Student’s t-test or one-way analysis of variance (ANOVA) was used to evaluate the significance of differences between groups. Correlations were detected by Spearman’s correlation coefficient. The relationship between POU6F2-AS1 expression and the clinicopathological parameters of CRC patients was calculated by the chi-square test or Fisher’s exact test. Overall survival (OS) was assessed by the Kaplan–Meier method and log-rank test. Univariate and multivariate Cox proportional hazard regression models were used to evaluate the effects of POU6F2-AS1 expression or other clinicopathological parameters on survival and the hazard ratio (HR).
Additional methods
Additional experiments, including cell transfection, RNA extraction, qRT–PCR, immunohistochemistry (IHC), RNA stability assays, fluorescence in situ hybridization (FISH) and immunofluorescence (IF) staining, Western blotting, Cell Counting Kit-8 (CCK–8) assays, 5-ethynyl-2’-deoxyuridine (EdU) assays, colony formation assays, methylated RNA Immunoprecipitation (MeRIP)-qPCR, RNA immunoprecipitation (RIP)-qPCR, RNA sequencing and bioinformatics analysis, application of palmitic acid (PA) and orlistat, and liquid chromatography mass spectrometry (LC-MS) based FA analysis are described in Additional file 2: Supplemental Materials and Methods.
Results
Identification of POU6F2-AS1 and analysis of its clinical significance in CRC
To identity novel oncogenic lncRNAs that are involved in the tumorigenesis and progression of CRC, four publicly available datasets (GSE126092, GSE134525, GSE109454, and GSE84983) were downloaded from the Gene Expression Omnibus (GEO) database for analysis (Fig. S1A-D). Overall, differential expression analysis of these datasets was performed with thresholds of |log FC| ≥2 and P value < 0.05, and a Venn diagram of the differentially expressed genes (DEGs) revealed that only the lncRNA POU6F2-AS1 was significantly upregulated in all four datasets (Fig. 1A). Next, the expression level of POU6F2-AS1 was detected in The Cancer Genome Atlas (TCGA) database. Given that POU6F2-AS1 is a novel lncRNA that has not been well characterized, we first analyzed the expression of POU6F2-AS1 in various types of cancer using the TCGA database. According to the results of the pan-cancer analysis, compared with those in corresponding normal tissues, the expression levels of POU6F2-AS1 were significantly dysregulated in many solid cancers, including colon adenocarcinoma (COAD) and rectum adenocarcinoma (READ) (Fig. 1B). Next, we narrowed our focus and performed on the comprehensive analysis of the expression and clinical significance of POU6F2-AS1 in CRC (COAD and READ in TCGA databases). The POU6F2-AS1 expression level in CRC tissues was significantly greater than that in normal tissues and was positively associated with advanced TNM stage (Fig. 1C-E). Moreover, analysis of data from the Cancer Cell Line Encyclopedia (CCLE) database revealed that the POU6F2-AS1 expression level in bowel cells was greater than that in most of the other cell lines from different tissues (Fig. S1E). Similar results were observed in CRC cell lines, POU6F2-AS1 was significantly upregulated in HCT116, SW480, DLD1, SW620, HT29 and LoVo CRC cells compared with FHC cells (human normal colorectal mucosal cells) (Fig. 1F).
To verify the POU6F2-AS1 expression pattern in CRC tissues, we first evaluated its expression in 84 paired fresh frozen tissue specimens using quantitative real–time polymerase chain reaction (qRT–PCR). qRT–PCR analysis revealed that the POU6F2-AS1 expression level was significantly greater in CRC tissues than in adjacent normal tissues (ANTs) (Fig. 1G). Moreover, fluorescence in situ hybridization (FISH) assays of a tissue microarray (TMA) of 60 pairs of CRC tissues and ANTs stained with the POU6F2-AS1 probe also showed that the expression of POU6F2-AS1 was significantly upregulated in CRC tissues (Fig. 1H, Fig. S1F, G). The results of qRT–PCR and FISH analyses were consistent with those of the bioinformatic analysis of the GEO and TCGA databases.
To further characterize the clinical significance of POU6F2-AS1 in the TMA cohort, a chi-square test was performed to analyse the relationship between POU6F2-AS1 expression and clinicopathological features. The results showed that POU6F2-AS1 expression was positively associated with tumour diameter (P = 0.038), depth of invasion (P < 0.001), lymph node metastasis (P = 0.010), and TNM stage (P = 0.010) (Additional file 1: Table S1). In addition, Kaplan–Meier (K–M) survival curve analysis and Cox regression analysis were performed to analyse the prognostic value of POU6F2-AS1 in CRC. K–M survival curves revealed that a higher level of POU6F2-AS1 expression was associated with poor overall survival (OS) in CRC patients in the TMA cohort (Fig. 1I). A similar result was verified in another cohort via analysis of the corresponding the GSE16125 dataset (Fig. 1J). Univariate Cox regression analysis suggested that POU6F2-AS1 expression, TNM stage, depth of invasion, lymph node metastasis, and distant metastasis were significant risk factors in terms of OS (Fig. 1K), while multivariate Cox analysis showed that POU6F2-AS1 expression and TNM stage could serve as independent prognostic factors in terms of OS in CRC patients (Fig. 1L).
POU6F2-AS1 promotes the growth of CRC cells and PDO in vitro
The previous analysis of endogenous POU6F2-AS1 expression in CRC cell lines revealed significant upregulation in HCT116 and SW480 cells and relatively low expression in LoVo cells (Fig. 1F). To determine the biological function of POU6F2-AS1 in CRC cells, we constructed two siRNAs (si-POU6F2-AS1#1 and si-POU6F2-AS1#2) and one shRNA (sh-POU6F2-AS1) targeting POU6F2-AS1 and one overexpression plasmid (POU6F2-AS1). Subsequently, we knocked down POU6F2-AS1 in HCT116 and SW480 cells and overexpressed it in LoVo cells (Fig. S1H–J). CCK–8, EdU and colony formation assays revealed that POU6F2-AS1 knockdown significantly suppressed the growth of CRC cells, whereas POU6F2-AS1 overexpression had the opposite effect (Fig. 2A–F). Furthermore, we generated patient-derived organoid (PDO) from CRC patient tissues and found that POU6F2-AS1 knockdown suppressed PDO formation, while POU6F2-AS1 overexpression facilitated PDO formation (Fig. 2G, H). These results clearly demonstrated that POU6F2-AS1 promotes the proliferation of CRC cells and PDO.
POU6F2-AS1-mediated de novo lipogenesis is essential for CRC cell and PDO growth in vitro
To understand the underlying biological pathways by which POU6F2-AS1 facilitates CRC tumorigenesis and growth, we performed RNA-seq on CRC cells transfected with the POU6F2-AS1 overexpression plasmid or empty vector and found a total of 2030 differentially expressed genes (DEGs) (P < 0.05, |log FC| ≥ 1) were detected (Fig. 3A). KEGG enrichment analysis revealed that the significantly enriched pathways included metabolic pathways, pathways in cancer, cell cycle, colorectal cancer, central carbon metabolism in cancer (Fig. 3B). Gene ontology biological process (GO-BP) enrichment analysis also revealed that the DEGs between POU6F2-AS1 high and low groups were associated metabolic processes, biological regulation, cell proliferation, growth, etc. (Fig. S1K). To gain further insight into the molecular mechanism of POU6F2-AS1 involved in the regulation of metabolic pathways, Gene Ontology (GO) analysis of the DEGs enriched in metabolic processes revealed that the top 20 significantly enriched terms associated with biological processes included lipid metabolic process and FA metabolic process (Fig. 3C).
Abnormal lipid metabolism reprogramming can provide abundant nutrients for tumorigenesis and cancer progression [19]. Therefore, we sought to explore the role of POU6F2-AS1 in the aberrant lipid metabolism of CRC cells. Intracellular staining with the lipophilic fluorescent dye BODIPY 493/503 and Oil Red O and an intracellular triglyceride (TAG) assay indicated that POU6F2-AS1 knockdown decreased lipid stores and TAG content, while POU6F2-AS1 overexpression increased lipid stores and TAG contents in CRC cells (Fig. 3D–H; Fig. S2A–C). More importantly, liquid chromatography–mass spectrometry (LC-MS) revealed that knockdown of POU6F2-AS1 significantly decreased FAs levels in CRC cells, such as palmitic acid (C16:0), stearic acid (C18:0), palmitoleic acid (C16:1), and oleic acid (C18:1) (Fig. 3I). These data demonstrate that that POU6F2-AS1 may facilitate the accumulation of lipids by positively regulating FA metabolism.
The increased FA metabolism of cancer cells could be achieved through de novo FA synthesis, FA uptake and decreased FA catabolism [20, 21]. To clarify how POU6F2-AS1 regulate FA metabolism, we detected the effect of POU6F2-AS1 on the expression levels of key enzymes involved in FA synthesis (FASN, ACC1, and SCD1), FA β-oxidation (CPT1A) and FA uptake catabolism (CD36). The results showed that the knockdown of POU6F2-AS1 significantly reduced the relative expression of FASN, ACC1, and SCD1 at both the mRNA and protein levels, while had no effect on the expression of CPT1A and CD36 (Fig. 3J-L). Consistently, overexpression of POU6F2-AS1 remarkably increased the mRNA and protein levels of the FA synthetic enzymes in CRC cells (Fig. 3M, N), demonstrating that POU6F2-AS1 promoted de novo FA synthesis of CRC cells.
We further evaluated whether the POU6F2-AS1-mediated modulation of FA metabolism is essential for CRC cell growth. We used orlistat, a FA synthase inhibitor, to decrease the lipid content in CRC cells overexpressing POU6F2-AS1; on the other hand, we used palmitic acid (PA) to increase the lipid content in CRC cells with POU6F2-AS1 knockdown. The results showed that orlistat treatment attenuated the effects of POU6F2-AS1 overexpression on cell growth, whereas PA treatment rescued the effects exerted by POU6F2-AS1 knockdown on cell growth (Fig. S2D-K). These findings confirmed the fundamental roles of POU6F2-AS1-mediated FA accumulation in driving CRC cell growth.
POU6F2-AS1 interacts with YBX1 to promote CRC cell growth and lipogenesis
The functions of lncRNAs are tightly related to their subcellular location. To further clarify the molecular mechanism of POU6F2-AS1 in CRC, FISH was performed, and the results showed that POU6F2-AS1 was located in both the nucleus and cytoplasm (Fig. S3A). Accumulating evidence has demonstrated that interactions of lncRNAs with proteins play vital roles in cancer progression [22,23,24]. Intriguingly, our GO molecular function (GO-MF) analysis revealed that the top enriched term was binding (Fig. S1K), which indicated that POU6F2-AS1 may exert its effects by binding to specific protein(s). We subsequently screened potential POU6F2-AS1-binding protein(s) via biotin-labelled RNA pull-down and subsequent mass spectrometry (MS) analyses, and approximately 64 potential interacting proteins were pulled down with POU6F2-AS1 in HCT116 cells (Fig. 4A and Additional file 1: Table S2). Overlapping analysis of RNA-binding proteins (RBPs) predicted by the online RBPDB database and RBPmap database revealed that only YBX1 was the potential partner of POU6F2-AS1 (Fig. 4B, C).
YBX1 is a DNA/RNA binding protein that dysregulates a wide range of genes involved in cancer apoptosis, proliferation, differentiation, and drug resistance [25]. Recent studies have indicated that YBX1 is upregulated in CRC and drives CRC progression [26,27,28]. Correlation analysis according to data from the TCGA database indicated that the expression of POU6F2-A1 was positively associated with that of YBX1 in CRC (Fig. S3B). Next, biotin-labelled RNA pull-down and RNA immunoprecipitation (RIP) and assays were used to validate the physical interaction between POU6F2-AS1 and YBX1 (Fig. 4D, E). FISH and immunofluorescence (IF) colocalization assays showed that POU6F2-AS1 and YBX1 colocalized in the nucleus and cytoplasm of CRC cells (Fig. 4F). Furthermore, the catRAPID analysis and representative YBX1 peptide detected by MS revealed a potential interaction region between POU6F2-AS1 and YBX1 (Fig. 4G, C). Moreover, the secondary structure of POU6F2-AS1 was predicted via RNAfold software (Fig. S3C). To further elucidate the binding regions between POU6F2-AS1 and YBX1, we constructed one biotinylated full-length POU6F2-AS1 construct and four biotinylated fragments of POU6F2-AS1 (F1: full-length sense of POU6F2-AS1; F2: full-length antisense of POU6F2-AS1; F3: Exon1 + 2, 1–116 nt; F4: Exon2 + 3, 14–195 nt; F5: Exon4, 196–381 nt) for RNA pull-down experiments with CRC cell lysates (Fig. 4H). Deletion mapping analysis revealed that the 117–195 nt region of POU6F2-AS1 was indispensable for its interaction with YBX1 (Fig. 4I). As previously reported, YBX1 consists of three domains, the N-terminal domain (NTD), the cold shock domain (CSD), and the large C-terminal domain (CTD) (Fig. 4J) [29]. To identify the domain involved in the binding of YBX1 to POU6F2-AS1, we cloned a series of Flag-tagged truncated YBX1 proteins (P1: 1–129 aa; P2: 130–205 aa; P3: 206–324 aa; and P4: 1–324 aa) (Fig. 4K). RIP assays were conducted with an anti-Flag antibody, and the results revealed that both the P3 domain and the P4 domain specifically interact with POU6F2-AS1 (Fig. 4L), suggesting that amino acids 206–324 in the CTD domain of YBX1 are required for the interaction with POU6F2-AS1.
Next, we characterized the role of POU6F2-AS1/YBX1 regulatory axis in CRC. First, we confirmed that YBX1 expression was markedly upregulated in CRC samples based on analysis of the TCGA database (Fig. S3D, E). Afterwards, knockdown of YBX1 by specific siRNA was validated by qRT–PCR and Western blotting (Fig. S3F, G). We knockdown YBX1 in LoVo cells with stable overexpression of POU6F2-AS1 and observed that POU6F2-AS1 overexpression promoted CRC cell growth and lipogenesis; these effects were attenuated by YBX1 knockdown (Fig. 4M-O). Consistently, the decreases of size of PDO and suppression of CRC cell growth and lipogenesis caused by POU6F2-AS1 knockdown were restored by YBX1 overexpression (Fig. S3H-J). Overall, these findings suggested that YBX1 plays a carcinogenic role in CRC cells and that POU6F2-AS1 may drive the growth and lipogenesis of CRC cells by directly binding to YBX1.
POU6F2-AS1 tethers YBX1 to the FASN promoter to transcriptionally activate FASN
To further elucidate the mechanism by which the POU6F2-AS1/YBX1 axis drives oncogenesis and lipogenesis in CRC and to identify the downstream target of the axis, we assessed the lipid metabolism-related genes (LMGs) coregulated by POU6F2-AS1 and YBX1. A Venn diagram revealed that four genes (FASN, INSIG1, OLAH, and SCD) were potential targets of the POU6F2-AS1/YBX1 axis (Fig. 5A). FASN is a key enzyme in lipid metabolism, especially in the de novo synthesis of FAs. Upregulated FASN expression has been reported to be associated with cell proliferation in multiple cancer types, including CRC [10, 30, 31]. In addition, Wang et al. reported that YBX1 depletion repressed the expression of FASN in glioblastoma cells, but the involved regulatory mechanism is still unknown [32]. Hence, we focused on the role of FASN in POU6F2-AS1/YBX1 axis-dependent CRC progression.
Correlation analysis according to the TCGA database demonstrated that FASN expression was positively associated with POU6F2-AS1 and YBX1 expression in CRC (Fig. S4A, B). Next, we determined whether POU6F2-AS1 could regulate the expression of FASN by binding to YBX1. Western blotting and qRT–PCR assays revealed that the decrease in FASN levels caused by POU6F2-AS1 knockdown could be reversed by YBX1 overexpression in HCT116 cells, while the increase in FASN levels induced by POU6F2-AS1 overexpression was attenuated by YBX1 knockdown in LoVo cells (Fig. 5B).
Next, we further investigated the regulatory mechanism by which the POU6F2-AS1/YBX1 axis enhances the FASN expression. GO analysis revealed that the most significantly enriched molecular function included transcription regulator activity (Fig. S1K). YBX1 is a well-established transcription factor that has been reported to induce the transcription of oncogenic genes, such as FOXA1 [33], FZD7 [34], and HOXC8 [35] in cancers. Intriguingly, previous studies have reported that transcriptional regulation is one of the most important mechanisms of FASN overexpression [36]. We used the JASPAR website to predict the underlying binding site (s) for YBX1 in the FASN promoter region and found 6 high-confidence response elements located within the FASN promoter region, these were named YBS1-6 (Fig. 5C). These data indicated that POU6F2-AS1 might exert its effects by recruiting YBX1 to the FASN promoter to induce transcriptional activation.
To evaluate whether YBX1 transcriptionally activates FASN, two different fragments of the FASN promoter were cloned and inserted into luciferase reporter vectors. The results of luciferase reporter assays showed that the luciferase activity of the reporter containing the − 447/-136 fragment upstream of the FASN transcription start site (TSS) was significantly decreased in HCT116 and SW480 cells with YBX1 knockdown but increased in cells with YBX1 overexpression (Fig. 5D, E; Fig. S4C). Moreover, neither knockdown nor overexpression of YBX1 altered the luciferase activity of the reported containing the − 1951/-952 fragment (Fig. 5D, E; Fig. S4C). These results confirmed that the YBX1 binding site (s) were located within the − 447/-136 region of the FASN promoter. To further determine whether POU6F2-AS1 activates FASN by recruiting YBX1 to its promoter, we detected the transcriptional activity of the FASN promoter after altering the expression of POU6F2-AS1. Luciferase assays showed that overexpressing YBX1 in POU6F2-AS1 knockdown cells markedly restored FASN promoter activation. On the other hand, knockdown of YBX1 in POU6F2-AS1 overexpressing cells significantly attenuated FASN promoter activation (Fig. 5F). To further clarify which specific site (s) YBX1 binds within the FASN promoter, 9 pairs of primers (denoted P1–9) were designed for different regions of the FASN promoter (Fig. 5G). Chromatin immunoprecipitation (ChIP) assays revealed that YBX1 binds to P2–4 regions of the FASN promoter (Fig. 5H). More importantly, knockdown of POU6F2-AS1 decreased YBX1 enrichment in the P2–4 regions of the FASN promoter, indicating that deletion of POU6F2-AS1 significantly decreased the binding of YBX1 to the FASN promoter (Fig. 5I). Taken together, these results show that POU6F2-AS1 tethers YBX1 to the FASN promoter to induce transcriptional activation of FASN.
FASN is involved in POU6F2-AS1-induced growth and lipogenesis in CRC
To further validate the role of FASN in POU6F2-AS1-mediated CRC progression, bioinformatics analysis was performed. FASN expression was markedly upregulated in CRC tissues vs. normal tissues in the TCGA databases, and these results were confirmed by IHC staining of CRC tissue and ANTs (Fig. S4D, E; Fig. 5J). Next, we constructed siRNA (si-FASN) targeting FASN and one overexpression plasmid (FASN) (Fig. S4F-H). Moreover, rescue experiments showed that the inhibitory effects of POU6F2-AS1 knockdown on the proliferation and lipogenesis of CRC cells and the growth of PDO (indicated by size) were significantly restored by ectopic overexpression of FASN (Fig. 5K-N. Fig. S4I-K). On the other hand, FASN knockdown efficiently attenuated promotion of proliferation and lipogenesis of CRC cells induced by POU6F2-AS1 overexpression (Fig. S4L-N). These results suggested that the promotion of CRC cell growth and lipogenesis by POU6F2-AS1 is partially dependent on FASN.
POU6F2-AS1 drives the growth and lipogenesis of CRC cells in vivo
To further verify the in vivo contribution of the POU6F2-AS1/FASN axis to promoting CRC cell growth and lipogenesis, we performed HCT116 cell growth assays and mouse tumour xenograft model experiments. Knockdown of POU6F2-AS1 markedly inhibited cell growth and decreased tumour volume and weight, but FASN overexpression reversed these effects induced by POU6F2-AS1 knockdown, as demonstrated by the increases in cell growth and tumour volume and weight (Fig. 6A-C). Furthermore, IHC and Oil Red O staining revealed that POU6F2-AS1 knockdown decreased FASN protein levels, suppressed tumour cell growth (as measured by Ki-67 staining) and decreased lipid content, while FASN overexpression in these cells restored FASN expression and increased lipogenesis, and elevated Ki-67 staining (Fig. 6D).
m6A modification is involved in the upregulation of POU6F2-AS1 in CRC
Nevertheless, the underlying mechanism leading to the increased expression of POU6F2-AS1 is still unclear. Recent studies have shown that the N6-methyladenosine (m6A) modification is involved in the stabilization and expression of lncRNAs in cancer [37, 38]. SRAMP website predictions indicated that several m6A sites are dispersed in the POU6F2-AS1 sequence (Fig. 7A). Therefore, we wondered whether the upregulation of POU6F2-AS1 in CRC is dependent on m6A modification. Given that methyltransferase-like 3 (METTL3) and insulin like growth factor 2 mRNA binding protein 2 (IGF2BP2) are the key m6A methyltransferase “writer” and “reader”, respectively, in CRC cells, we evaluated METTL3 and IGF2BP2 expression and their correlation with POU6F2-AS1 in the TCGA database. Notably, the expression of METTL3 and IGF2BP2 was significantly upregulated in CRC tissues (Fig. 7B, C) and was positively correlated with POU6F2-AS1 expression (R = 0.6 and 0.51; Fig. 7D, E). IHC staining of TMAs further confirmed that POU6F2-AS1 expression was positively correlated with METTL3 and FASN expression (Fig. 7F, G). In addition, METTL3 shRNAs and overexpression plasmids were constructed and transfected into cells (Fig. S4O-Q). POU6F2-AS1 expression was apparently reduced with METTL3 knockdown, but increased with METTL3 overexpression (Fig. 7H; Fig. S4Q). MeRIP assays also revealed that METTL3 knockdown reduced the m6A modification level of POU6F2-AS1 in SW480 cells (Fig. 7I). Furthermore, we detected the stability of POU6F2-AS1 in SW480 cells upon METTL3 knockdown or overexpression. In cells treated with actinomycin D to block the de novo synthesis of RNA, METTL3 knockdown decreased the stability of POU6F2-AS1, whereas METTL3 overexpression had the opposite effect (Fig. 7J). Next, we elucidated to the mechanism involved in the METTL3-induced stabilization of POU6F2-AS1 in CRC cells. IGF2BP2, a well-known “reader” of m6A, is also an RNA binding protein that plays a vital role in controlling the stability of mRNAs [39]. Knockdown of IGF2BP2 in SW480 cells markedly decreased POU6F2-AS1 expression (Fig. 7K). Importantly, knockdown of IGF2BP2 decreased the half-life of POU6F2-AS1 in SW480 cells treated with actinomycin D (Fig. 7L). In addition, the interaction between IGF2BP2 and POU6F2-AS1 was inhibited after METTL3 knockdown (Fig. 7M). Overall, our data indicated that METTL3-mediated m6A modification promotes POU6F2-AS1 expression in CRC via IGF2BP2-dependent stabilization of RNA.
Discussion
Transcriptome studies have demonstrated that protein-coding RNAs account for less than 2% of the transcribed human genome, and most of the remaining RNAs are noncoding RNAs (ncRNAs) [40]. Currently, lncRNAs, which are functional ncRNAs, have attracted increasing amounts of attention in cancer research. Several phase I/II clinical trials targeting lncRNAs have been conducted to treat patients with cancer [41]. The development of aberrantly expressed lncRNAs as therapeutic tools has yielded exciting prospects. Moreover, with the widespread use and improvement of high-throughput sequencing, additional opportunities have been provided for bioinformatics-based research and clinical validation of promising lncRNAs. In the present study, we identified the novel upregulated lncRNA POU6F2-AS1 in CRC via comprehensive analysis of publicly available datasets and validated the upregulation of POU6F2-AS1 in both CRC tissues and cell lines. Moreover, we found that high expression of POU6F2-AS1 was positively associated with advanced clinicopathological stage, large tumour size, and poor OS in patients with CRC, suggesting that POU6F2-AS1 might play an essential role in CRC tumorigenesis and progression. Previous study demonstrated that POU6F2-AS1 are upregulated in lung cancer and CRC and is an oncogenic factor [42, 43]. However, the function and mechanism of POU6F2-AS1 in lipid metabolism has not been explored in CRC.
Gain- and loss-of-function experiments demonstrated that POU6F2-AS1 promotes the growth of CRC cells in vitro and in vivo. In addition, RNA sequencing analysis, BODIPY and Oil Red O staining and LC/MS data revealed that POU6F2-AS1 was able to affect lipid metabolism and FA metabolism. Reprogramming of lipid metabolism is a newly well-documented hallmark of cancer cells, and represents an attractive vulnerability that might be exploited as s therapeutic target [6]. Proliferating cancer cells are likely to obtain energy via lipid reprogramming particularly through the synthesis of FAs; overall, FA metabolism involves FA uptake, de novo synthesis, and β‐oxidation to support unlimited cellular proliferation [19, 20, 44]. In addition, cancer cells can utilize FAs to provide ingredients for synthesizing cellular organelles and membranes during rapid replication [45].
Multiple metabolic pathways contribute to lipid accumulation in cancer cells, and most of these pathways are related to the rate-limiting lipogenic enzymes, whose upregulation induces FA synthesis [46]. Accumulating evidence has shown that multiple lipogenic enzymes, such as FASN, ACC, SCD1, and ACLY, are aberrantly upregulated in a wide variety of cancers, driving FA synthesis and contributing to lipogenesis and the progression of cancer cells [47,48,49,50]. Consistently, the expression of lipogenic enzymes is also markedly increased in CRC [51]. Recently, lncRNAs were reported to participate in the regulation of lipogenic enzymes, thereby promoting the progression of cancers. For example, Zheng et al. revealed that the lncRNA TINCR could bind to ACLY and maintain its protein stability by protecting it from ubiquitin-mediated degradation, thereby facilitating de novo lipid biosynthesis and progression of nasopharyngeal carcinoma [52]. Jia et al. reported that that the lncRNA NEAT1 facilitated FA metabolism in gastric cancer via the c‑Jun/c‑Fos/SREBP1 axis [53]. Peng et al. demonstrated that the lncRNA FASRL could bind to ACACA and inhibit its phosphorylation, thus increasing FA synthesis in hepatocellular carcinoma [54]. However, the crosstalk between lncRNAs and metabolic enzymes or processes in CRC is poorly understood. In the present study, our findings revealed that POU6F2-AS1 overexpression induces lipogenesis and proliferation in CRC cells by upregulating FASN, which is the key rate‐limiting lipogenic enzyme responsible for the terminal catalytic step in the synthesis of FAs. FASN is often upregulated in cancer cells and its depletion results in antitumour effects; thus, FASN has been defined as a promising therapeutic target [55]. Consistent with this conclusion, we found that FASN knockdown could efficiently attenuated CRC cell lipogenesis and proliferation, while these processed were promoted by POU6F2-AS1 overexpression. These findings indicated cotreatment with POU6F2-AS1 siRNAs may improve the therapeutic efficacy of FASN inhibitors. However, the antitumour effects of combination therapies based on POU6F2-AS1 need to be further studied in the future.
LncRNAs typically exert their functions by acting as microRNA (miRNA) sponges and interacting with RBPs [56, 57]. In the present study, by overlapping the MS analysis and two online RBP prediction results, we determined that POU6F2-AS1 specifically interacts with YBX1. Deletion-mapping analysis revealed that the 117–195 nt region of POU6F2-AS1 was responsible for the interaction between these two genes, which is consistent with the regions predicted by catRAPID. YBX1 has been shown to play vital roles in the tumorigenesis and progression of malignancies, including CRC. Accumulating evidence has indicated that YBX1 is a vital functional partner of multiple lncRNAs in cancer. For instance, YBX1 is an interactive partner of lincNMR that regulates the expression of downstream RRM2, TYMS and TK1 by binding to their promoters, governing nucleotide metabolism and cell proliferation in cancer [58]. The lncRNA MILIP bound with YBX1 to increase the translation of Snail, leading to the enhanced metastasis of clear cell renal cell carcinoma cells [59]. LncRNA RMRP recruits YBX1 to upregulate the transcription of TGFBR1, thereby promoting the proliferation and progression of non-small cell lung cancer [60]. Previous studies demonstrated that YBX1, a well-established transcription factor, may bind to the promoters of certain genes, thereby initiating their transcriptional activation. For example, NOTCH3 [61], CBX3 [62], and ANXA8 [63] are regulated by YBX1 via transcriptional activation. In this study, we found that POU6F2-AS1 tethers YBX1 to the promoter of FASN to promote FASN activation. Furthermore, ChIP assays revealed that the region spanning from − 583 to -10 nt in the FASN promoter contains the binding sites for YBX1. Coincidently, YBX1 binding site-2 (YBX-2), YBX-3, YBX-5, and YBX-6, which were predicted by the JASPAR website, were also located in this region. These findings are consistent with conclusions reported in the literature that transcriptional regulation is one of the most important mechanisms of FASN overexpression in cancer cells [36]. However, we did not detect whether histone proteins or other transcription factors in this region cooperatively interact with YBX1 to initiate FASN transcription. Additionally, as a multifunctional DNA/RNA binding protein RNA, YBX1 also participates in the regulation of multiple DNA/RNA-dependent events, including pre-mRNA splicing, mRNA packaging, translation and stabilization, etc [29, 64]. However, whether YBX1 cooperates with POU6F2-AS1 to regulate FASN via these mechanisms needs to be explored in the future.
In addition, our work addressed the mechanism of POU6F2-AS1 upregulation in CRC. The m6A modification, a type of reversible epigenetic regulatory mechanism, has been widely reported to be involved in the stabilization and expression of RNAs, including lncRNAs, in cancer [65, 66]. For example, METTL3-mediated m6A modification increases the stabilization and expression of THAP7-AS1 in an IGF2BP1-dependent manner [37]. METTL3 catalyzes m6A modification, and YTHDC1 recognizes and stabilizes m6A to increase expression of TERRA [67]. Since the bioinformatics analysis predicted high-confidence m6A sites in the sequence of POU6F2-AS1, we further assessed whether POU6F2-AS1 was an m6A-incucible lncRNA. Coincidently, we observed that knockdown of METTL3 suppressed the expression and stability of POU6F2-AS1 in CRC cells and decreased m6A level of POU6F2-AS1. IGF2BP2 has been shown to mediate the stability of m6A-modified RNAs, such as CREB1 [68], HMGA1 [69], MSX1 and JARID2 [70] in CRC. Here, we also found that IGF2BP2 could recognize m6A sites and maintain the stability of POU6F2-AS1. Overall, these results suggested that METTL3-catalyzed and IGF2BP2- stabilized m6A modification may be involved in the upregulation of POU6F2-AS1 in CRC.
Conclusions
In summary, POU6F2-AS1 is a novel upregulated lncRNA in CRC that is stabilized by METTL3-mediated m6A modification. Moreover, POU6F2-AS1 specifically tethers YBX1 to the FASN promoter to activate its transcription, leading to the enhanced expression of this key enzyme in the de novo synthesis of FAs, thereby facilitating the progression of CRC. Thus, POU6F2-AS1 plays a pivotal role in the oncogenic METTL3/POU6F2-AS1/YBX1/FASN axis to promote CRC cell lipogenesis and growth, which provides new insight into lipid accumulation in CRC and highlights, POU6F2-AS1 as a promising therapeutic target for patients with CRC.
Data availability
The RAN-seq datasets of dysregulated lncRNAs in CRC were obtained from GEO database (https://www.ncbi.nlm.nih.gov/geo/) under the correspond accession number. The public available datasets involved in this study can be obtained from online repositories. All datasets involved in the present study are available from the corresponding author on reasonable request
Abbreviations
- CRC:
-
Colorectal cancer
- TCGA:
-
The Cancer Genome Atlas
- LncRNA:
-
long non-coding RNA
- FA:
-
Fatty acid
- YBX1:
-
Y-box binding protein 1
- FASN:
-
Fatty Acid Synthase
- ROC:
-
Receiver operating characteristic
- HE:
-
Hematoxylin eosin
- OS:
-
overall survival
- KEGG:
-
Kyoto encyclopedia of genes and genomes
- shRNA:
-
short hairpin RNA
- siRNA:
-
small interfering RNA
- GEPIA:
-
The Gene Expression Profiling Interactive Analysis
- ANTs:
-
adjacent noncancerous tissues
- qRT–PCR:
-
Real-time quantitative polymerase chain reaction
- CCK–8:
-
Cell Counting Kit–8
- FISH:
-
RNA fluorescence in situ hybridization
- IHC:
-
Immunohistochemical staining
- Edu:
-
5-Ethynyl-2’-deoxyuridine
- TAG:
-
Triglyceride
- RIP:
-
RNA immunoprecipitation
- Me-RIP:
-
Methylated RNA Immunoprecipitation
- ChIP:
-
Chromatin immunoprecipitation
- PDO:
-
Patient-derived organoid
- GEO:
-
Gene-Expression Omnibus
- LMGs:
-
Lipid metabolism-related genes
- PA:
-
Palmitic acid
- LC-MS:
-
Liquid chromatography mass spectrometry
References
Siegel RL, Miller KD, Wagle NS, Jemal A. Cancer statistics, 2023. CA Cancer J Clin. 2023;73:17–48.
Siegel RL, Wagle NS, Cercek A, Smith RA, Jemal A. Colorectal cancer statistics, 2023. CA Cancer J Clin. 2023;73:233–54.
Andrei P, Battuello P, Grasso G, Rovera E, Tesio N, Bardelli A. Integrated approaches for precision oncology in colorectal cancer: the more you know, the better. Semin Cancer Biol. 2022;84:199–213.
Guzman Y, Rios J, Paredes J, Dominguez P, Maurel J, Gonzalez-Abos C, Otero-Pineiro A, Almenara R, Ladra M, Prada B, et al. Time interval between the end of Neoadjuvant Therapy and Elective Resection of locally advanced rectal Cancer in the CRONOS Study. JAMA Surg. 2023;158:910–9.
Sedlak JC, Yilmaz OH, Roper J. Metabolism and colorectal Cancer. Annu Rev Pathol. 2023;18:467–92.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646–74.
La Vecchia S, Sebastian C. Metabolic pathways regulating colorectal cancer initiation and progression. Semin Cell Dev Biol. 2020;98:63–70.
Chen D, Zhou X, Yan P, Yang C, Li Y, Han L, Ren X. Lipid metabolism reprogramming in colorectal cancer. J Cell Biochem. 2023;124:3–16.
Fhu CW, Ali A. Fatty acid synthase: an emerging target in Cancer. Molecules 2020, 25.
Wei W, Qin B, Wen W, Zhang B, Luo H, Wang Y, Xu H, Xie X, Liu S, Jiang X, et al. FBXW7beta loss-of-function enhances FASN-mediated lipogenesis and promotes colorectal cancer growth. Signal Transduct Target Ther. 2023;8:187.
Quinn JJ, Chang HY. Unique features of long non-coding RNA biogenesis and function. Nat Rev Genet. 2016;17:47–62.
Choi SW, Kim HW, Nam JW. The small peptide world in long noncoding RNAs. Brief Bioinform. 2019;20:1853–64.
Yang Y, Yan X, Li X, Ma Y, Goel A. Long non-coding RNAs in colorectal cancer: novel oncogenic mechanisms and promising clinical applications. Cancer Lett. 2021;504:67–80.
Liu B, Song A, Gui P, Wang J, Pan Y, Li C, Li S, Zhang Y, Jiang T, Xu Y, et al. Long noncoding RNA LINC01594 inhibits the CELF6-mediated splicing of oncogenic CD44 variants to promote colorectal cancer metastasis. Cell Death Dis. 2023;14:427.
Yang T, Chen WC, Shi PC, Liu MR, Jiang T, Song H, Wang JQ, Fan RZ, Pei DS, Song J. Long noncoding RNA MAPKAPK5-AS1 promotes colorectal cancer progression by cis-regulating the nearby gene MK5 and acting as a let-7f-1-3p sponge. J Exp Clin Cancer Res. 2020;39:139.
Cheng Z, Wang J, Xu Y, Jiang T, Xue Z, Li S, Zhao Y, Song H, Song J. N7-methylguanosine-related lncRNAs: distinction between hot and cold tumors and construction of predictive models in colon adenocarcinoma. Front Oncol. 2022;12:951452.
Wang H, Chen Y, Liu Y, Li Q, Luo J, Wang L, Chen Y, Sang C, Zhang W, Ge X, et al. The lncRNA ZFAS1 regulates lipogenesis in colorectal cancer by binding polyadenylate-binding protein 2 to stabilize SREBP1 mRNA. Mol Ther Nucleic Acids. 2022;27:363–74.
Christensen LL, True K, Hamilton MP, Nielsen MM, Damas ND, Damgaard CK, Ongen H, Dermitzakis E, Bramsen JB, Pedersen JS, et al. SNHG16 is regulated by the wnt pathway in colorectal cancer and affects genes involved in lipid metabolism. Mol Oncol. 2016;10:1266–82.
Cheng C, Geng F, Cheng X, Guo D. Lipid metabolism reprogramming and its potential targets in cancer. Cancer Commun (Lond). 2018;38:27.
Rohrig F, Schulze A. The multifaceted roles of fatty acid synthesis in cancer. Nat Rev Cancer. 2016;16:732–49.
Cao D, Yang J, Deng Y, Su M, Wang Y, Feng X, Xiong Y, Bai E, Duan Y, Huang Y. Discovery of a mammalian FASN inhibitor against xenografts of non-small cell lung cancer and melanoma. Signal Transduct Target Ther. 2022;7:273.
Yao ZT, Yang YM, Sun MM, He Y, Liao L, Chen KS, Li B. New insights into the interplay between long non-coding RNAs and RNA-binding proteins in cancer. Cancer Commun (Lond). 2022;42:117–40.
Xia A, Yuan W, Wang Q, Xu J, Gu Y, Zhang L, Chen C, Wang Z, Wu D, He Q, et al. The cancer-testis lncRNA lnc-CTHCC promotes hepatocellular carcinogenesis by binding hnRNP K and activating YAP1 transcription. Nat Cancer. 2022;3:203–18.
Ren L, Fang X, Shrestha SM, Ji Q, Ye H, Liang Y, Liu Y, Feng Y, Dong J, Shi R. LncRNA SNHG16 promotes development of oesophageal squamous cell carcinoma by interacting with EIF4A3 and modulating RhoU mRNA stability. Cell Mol Biol Lett. 2022;27:89.
Kuwano M, Shibata T, Watari K, Ono M. Oncogenic Y-box binding protein-1 as an effective therapeutic target in drug-resistant cancer. Cancer Sci. 2019;110:1536–43.
Wu QN, Luo XJ, Liu J, Lu YX, Wang Y, Qi J, Liu ZX, Huang QT, Liu ZK, Lu JB, et al. MYC-Activated LncRNA MNX1-AS1 promotes the progression of Colorectal Cancer by stabilizing YB1. Cancer Res. 2021;81:2636–50.
Chen S, Li K, Guo J, Chen HN, Ming Y, Jin Y, Xu F, Zhang T, Yang Y, Ye Z, et al. circNEIL3 inhibits tumor metastasis through recruiting the E3 ubiquitin ligase Nedd4L to degrade YBX1. Proc Natl Acad Sci U S A. 2023;120:e2215132120.
Kim A, Shim S, Kim YH, Kim MJ, Park S, Myung JK. Inhibition of Y Box binding protein 1 suppresses cell growth and motility in Colorectal Cancer. Mol Cancer Ther. 2020;19:479–89.
Lyabin DN, Eliseeva IA, Ovchinnikov LP. YB-1 protein: functions and regulation. Wiley Interdiscip Rev RNA. 2014;5:95–110.
Li H, Chen Z, Zhang Y, Yuan P, Liu J, Ding L, Ye Q. MiR-4310 regulates hepatocellular carcinoma growth and metastasis through lipid synthesis. Cancer Lett. 2021;519:161–71.
Hu N, Li Y, Zhao Y, Wang Q, You JC, Zhang XD, Ye LH. A novel positive feedback loop involving FASN/p-ERK1/2/5-LOX/LTB4/FASN sustains high growth of breast cancer cells. Acta Pharmacol Sin. 2011;32:921–9.
Wang JZ, Zhu H, You P, Liu H, Wang WK, Fan X, Yang Y, Xu K, Zhu Y, Li Q et al. Upregulated YB-1 protein promotes glioblastoma growth through a YB-1/CCT4/mLST8/mTOR pathway. J Clin Invest 2022, 132.
Zhang Y, Huang YX, Wang DL, Yang B, Yan HY, Lin LH, Li Y, Chen J, Xie LM, Huang YS, et al. LncRNA DSCAM-AS1 interacts with YBX1 to promote cancer progression by forming a positive feedback loop that activates FOXA1 transcription network. Theranostics. 2020;10:10823–37.
Chen Q, Wang H, Li Z, Li F, Liang L, Zou Y, Shen H, Li J, Xia Y, Cheng Z, et al. Circular RNA ACTN4 promotes intrahepatic cholangiocarcinoma progression by recruiting YBX1 to initiate FZD7 transcription. J Hepatol. 2022;76:135–47.
Su H, Fan G, Huang J, Qiu X. LncRNA HOXC-AS3 promotes non-small-cell lung cancer growth and metastasis through upregulation of YBX1. Cell Death Dis. 2022;13:307.
Mashima T, Seimiya H, Tsuruo T. De novo fatty-acid synthesis and related pathways as molecular targets for cancer therapy. Br J Cancer. 2009;100:1369–72.
Liu HT, Zou YX, Zhu WJ, Sen L, Zhang GH, Ma RR, Guo XY, Gao P. lncRNA THAP7-AS1, transcriptionally activated by SP1 and post-transcriptionally stabilized by METTL3-mediated m6A modification, exerts oncogenic properties by improving CUL4B entry into the nucleus. Cell Death Differ. 2022;29:627–41.
Xie J, Zhang H, Wang K, Ni J, Ma X, Khoury CJ, Prifti V, Hoard B, Cerenzia EG, Yin L, et al. M6A-mediated-upregulation of lncRNA BLACAT3 promotes bladder cancer angiogenesis and hematogenous metastasis through YBX3 nuclear shuttling and enhancing NCF2 transcription. Oncogene. 2023;42:2956–70.
Zhu TY, Hong LL, Ling ZQ. Oncofetal protein IGF2BPs in human cancer: functions, mechanisms and therapeutic potential. Biomark Res. 2023;11:62.
Mattick JS, Rinn JL. Discovery and annotation of long noncoding RNAs. Nat Struct Mol Biol. 2015;22:5–7.
Adams BD, Parsons C, Walker L, Zhang WC, Slack FJ. Targeting noncoding RNAs in disease. J Clin Invest. 2017;127:761–71.
Wu XY, Xie Y, Zhou LY, Zhao YY, Zhang J, Zhang XF, Guo S, Yu XY. Long noncoding RNA POU6F2-AS1 regulates lung cancer aggressiveness through sponging miR-34c-5p to modulate KCNJ4 expression. Genet Mol Biol. 2021;44:e20200050.
Lu D, Chen A. lncRNA POU6F2-AS1 regulated by KIAA1429 contributes to Colorectal Cancer Progression in an m(6)a modification manner. Mol Biotechnol 2023.
Tang Z, Xu Z, Zhu X, Zhang J. New insights into molecules and pathways of cancer metabolism and therapeutic implications. Cancer Commun (Lond). 2021;41:16–36.
Broadfield LA, Pane AA, Talebi A, Swinnen JV, Fendt SM. Lipid metabolism in cancer: new perspectives and emerging mechanisms. Dev Cell. 2021;56:1363–93.
Koundouros N, Poulogiannis G. Reprogramming of fatty acid metabolism in cancer. Br J Cancer. 2020;122:4–22.
Zhang J, Song Y, Shi Q, Fu L. Research progress on FASN and MGLL in the regulation of abnormal lipid metabolism and the relationship between tumor invasion and metastasis. Front Med. 2021;15:649–56.
Guo Z, Huo X, Li X, Jiang C, Xue L. Advances in regulation and function of stearoyl-CoA desaturase 1 in cancer, from bench to bed. Sci China Life Sci 2023.
Icard P, Wu Z, Fournel L, Coquerel A, Lincet H, Alifano M. ATP citrate lyase: a central metabolic enzyme in cancer. Cancer Lett. 2020;471:125–34.
Mounier C, Bouraoui L, Rassart E. Lipogenesis in cancer progression (review). Int J Oncol. 2014;45:485–92.
Ecker J, Benedetti E, Kindt ASD, Horing M, Perl M, Machmuller AC, Sichler A, Plagge J, Wang Y, Zeissig S, et al. The Colorectal Cancer Lipidome: identification of a robust tumor-specific lipid species signature. Gastroenterology. 2021;161:910–923e919.
Zheng ZQ, Li ZX, Guan JL, Liu X, Li JY, Chen Y, Lin L, Kou J, Lv JW, Zhang LL, et al. Long noncoding RNA TINCR-Mediated regulation of Acetyl-CoA metabolism promotes nasopharyngeal carcinoma progression and Chemoresistance. Cancer Res. 2020;80:5174–88.
Jia Y, Yan Q, Zheng Y, Li L, Zhang B, Chang Z, Wang Z, Tang H, Qin Y, Guan XY. Long non-coding RNA NEAT1 mediated RPRD1B stability facilitates fatty acid metabolism and lymph node metastasis via c-Jun/c-Fos/SREBP1 axis in gastric cancer. J Exp Clin Cancer Res. 2022;41:287.
Peng JY, Cai DK, Zeng RL, Zhang CY, Li GC, Chen SF, Yuan XQ, Peng L. Upregulation of Superenhancer-Driven LncRNA FASRL by USF1 promotes De Novo fatty acid biosynthesis to Exacerbate Hepatocellular Carcinoma. Adv Sci (Weinh). 2022;10:e2204711.
Menendez JA, Lupu R. Fatty acid synthase (FASN) as a therapeutic target in breast cancer. Expert Opin Ther Targets. 2017;21:1001–16.
Ferre F, Colantoni A, Helmer-Citterich M. Revealing protein-lncRNA interaction. Brief Bioinform. 2016;17:106–16.
Zhou Y, Shao Y, Hu W, Zhang J, Shi Y, Kong X, Jiang J. A novel long noncoding RNA SP100-AS1 induces radioresistance of colorectal cancer via sponging miR-622 and stabilizing ATG3. Cell Death Differ. 2023;30:111–24.
Gandhi M, Gross M, Holler JM, Coggins SA, Patil N, Leupold JH, Munschauer M, Schenone M, Hartigan CR, Allgayer H, et al. The lncRNA lincNMR regulates nucleotide metabolism via a YBX1 - RRM2 axis in cancer. Nat Commun. 2020;11:3214.
Wang Y, Feng YC, Gan Y, Teng L, Wang L, La T, Wang P, Gu Y, Yan L, Li N, et al. LncRNA MILIP links YBX1 to translational activation of Snai1 and promotes metastasis in clear cell renal cell carcinoma. J Exp Clin Cancer Res. 2022;41:260.
Yin H, Chen L, Piao S, Wang Y, Li Z, Lin Y, Tang X, Zhang H, Zhang H, Wang X. M6A RNA methylation-mediated RMRP stability renders proliferation and progression of non-small cell lung cancer through regulating TGFBR1/SMAD2/SMAD3 pathway. Cell Death Differ. 2023;30:605–17.
Li Y, Gao Y, Niu X, Tang M, Li J, Song B, Guan X. LncRNA BASP1-AS1 interacts with YBX1 to regulate notch transcription and drives the malignancy of melanoma. Cancer Sci. 2021;112:4526–42.
Zhang H, Yu H, Ren D, Sun Y, Guo F, Cai H, Zhou C, Zhou Y, Jin X, Wu H. CBX3 regulated by YBX1 promotes smoking-induced pancreatic Cancer progression via inhibiting SMURF2 expression. Int J Biol Sci. 2022;18:3484–97.
Lu X, Wang J, Wang W, Lu C, Qu T, He X, Liu X, Guo R, Zhang E. Copy number amplification and SP1-activated lncRNA MELTF-AS1 regulates tumorigenesis by driving phase separation of YBX1 to activate ANXA8 in non-small cell lung cancer. Oncogene. 2022;41:3222–38.
Mordovkina D, Lyabin DN, Smolin EA, Sogorina EM, Ovchinnikov LP, Eliseeva I. Y-Box binding proteins in mRNP Assembly, translation, and Stability Control. Biomolecules 2020, 10.
Wang Z, Zhou J, Zhang H, Ge L, Li J, Wang H: RNA m(6) A methylation in cancer. Mol Oncol 2023;17:195–229.
Fu Y, Dominissini D, Rechavi G, He C. Gene expression regulation mediated through reversible m(6)a RNA methylation. Nat Rev Genet. 2014;15:293–306.
Chen L, Zhang C, Ma W, Huang J, Zhao Y, Liu H. METTL3-mediated m6A modification stabilizes TERRA and maintains telomere stability. Nucleic Acids Res. 2022;50:11619–34.
Yao B, Zhang Q, Yang Z, An F, Nie H, Wang H, Yang C, Sun J, Chen K, Zhou J, et al. CircEZH2/miR-133b/IGF2BP2 aggravates colorectal cancer progression via enhancing the stability of m(6)A-modified CREB1 mRNA. Mol Cancer. 2022;21:140.
Hou P, Meng S, Li M, Lin T, Chu S, Li Z, Zheng J, Gu Y, Bai J. LINC00460/DHX9/IGF2BP2 complex promotes colorectal cancer proliferation and metastasis by mediating HMGA1 mRNA stability depending on m6A modification. J Exp Clin Cancer Res. 2021;40:52.
Huizhe, Wu Xiangyu, Ding Xiaoyun, Hu Qing, Zhao Qiuchen, Chen Tong, Sun Yalun, Li Hao, Guo Meng, Li Ziming, Gao Weifan, Yao Lin, Zhao Kai, Li Minjie, Wei (2022) LINC01021 maintains tumorigenicity by enhancing N6-methyladenosine reader IMP2 dependent stabilization of MSX1 and JARID2: implication in colorectal cancer. Oncogene. 2022;41(13):1959–1973. https://doi.org/10.1038/s41388-022-02189-x
Acknowledgements
We would like to all individuals all individuals, including the patients, who participated in this work. We also want to express our appreciation to the team who provided the public availably GEO database, making scientific research boundaryless.
Funding
Our study was funded by grants from the National Natural Science Foundation of China (No. 82203486, 82073133), the Natural Science Foundation of Jiangsu Province (BK20231159), the Scientific Research of Jiangsu Health Committee (ZDA2020005), the ‘Six Talents Peak’ High-level Talent Project of Jiangsu Province (WSW-050), and the Xuzhou Medical Leading Talents Training Project (XWRCHT20210034).
Author information
Authors and Affiliations
Contributions
J.S. and T.J. contributed to the conception and design of the study. T.J., W.J.Q., Y.Z.X. and W.B.L. prepared the figures and drafted the manuscript. J.S. reviewed and revised the manuscript. T.J., W.J.Q., Y.Z.X., W.B.L., H.Q.H., and Q.J.L. performed the in vitro cell experiments. R.J. and X.Y.X. generated the PDOs. WJQ, H.Q.H. and Q.J.L. performed the in vivo experiments. W.J.Q, Q.Y.L and T.J. conducted the bioinformatics analysis. Z.R.F., T.X., X.Y.X. and HS participated in the collection, follow-up, detection and data analysis of clinical specimens. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was conducted in compliance with the declaration of Helsinki and received informed consent from all human subjects. The Ethnic Committee of the Affiliated Hospital of Xuzhou Medical University provided the ethics approval statements for human subjects in this study. All animal experiments were carried out according to the Health Guide for the Care and Use of Laboratory Animals approved by the Animal Experimental Research Ethics Committee of Xuzhou Medical University.
Consent for publication
All authors agree to publication of the article.
Competing interest
There are no competing interests between all authors.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Supplementary Material 1
: Original western blot images
Supplementary Material 2
: Fig. S1. (A-D) Volcano plots of the differentially expressed genes for GSE126092, GSE134525, GSE109454, and GSE84983. (E) Analysis of POU6F2-AS1 expression in cell lines of different tissue origins in the Cancer Cell Line Encyclopedia (CCLE) database. (F) FISH staining intensities of POU6F2-AS1 in CRC tissues compared with paired ANTs. T, tumour tissues; N, paired adjacent non-cancerous tissues. (G) The FISH score of POU6F2-AS1 in CRC tissues was remarkably higher than that in paired ANTs. (H, I) Validation of POU6F2-AS1 expression in HCT116 and SW480 cells transfected with siRNAs. (J) Validation of POU6F2-AS1 expression in LoVo cells transfected with overexpression plasmids. (K) Gene Ontology (GO) analysis including biological process (BP), cellular component (CC) and molecular function (MF) exhibited the significantly enriched pathways after POU6F2-AS1 overexpression. *P < 0.05, **P < 0.01, ***P < 0.001
Supplementary Material 3
: Fig. S2. (A-C) Relative triglyceride levels in CRC cells with POU6F2-AS1 knockdown or overexpression. (D-H) CCK–8, EdU and colony formation analysis of PA in HCT116 and SW480 cells with POU6F2-AS1 knockdown. (I-K) CCK–8, EdU and colony formation analysis of orlistat in LoVo cells with POU6F2-AS1 overexpression. The data are presented as the mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001.
Supplementary Material 4
: Fig. S3. (A) Representative RNA-FISH images showing that the localization of POU6F2-AS1 in CRC cells. POU6F2-AS1 probes are red, and nuclei are stained with DAPI. (B) Correlation of POU6F2-AS1 and YBX1 expression in CRC was analyzed by GEPIA 2 website. (C) The secondary structure of POU6F2-AS1 predicted by RNA fold website. (D) The YBX1 expression between CRC tissues and normal samples in TCGA database. (E) The YBX1 expression between CRC tissues and paired normal samples in TCGA database. (F, G) Validation of YBX1 expression by qRT–PCR and western blotting in LoVo cells transfected with siRNAs. (H) CRC patient-derived organoid infected with POU6F2-AS1 knockdown or YBX1 overexpression lentivirus. (I) Colony formation assay in POU6F2-AS1 knockdown of HCT116 cells with YBX1 overexpression. (J) Oil Red O staining assay in POU6F2-AS1 knockdown of HCT116 cells with YBX1 overexpression. The data are presented as the mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001.
Supplementary Material 5
: Fig. S4. (A, B) Correlation of FASN and POU6F2-AS1 or YBX1 expression in CRC was analyzed by GEPIA 2 website. (C) Luciferase reporter assays showed that the luciferase activity driven by the -447/-136 fragment of the FASN promoter region was significantly reduced in SW480 cells with YBX1 knockdown. (D, E) The FASN expression between CRC tissues and normal samples in TCGA database. (F, G) Validation of FASN expression by qRT–PCR and western blotting in LoVo cells transfected with siRNAs. (H) Validation of FASN expression by qRT–PCR in HCT116 cells transfected with FASN overexpression plasmid. (I-K) CCK–8, colony formation and Oil Red O staining assays with POU6F2-AS1 knockdown in SW480 cells with FASN overexpression. (L-N) Colony formation, CCK–8 and Oil Red O staining assays in POU6F2-AS1 overexpression LoVo cells with FASN knockdown. (O, P) Western blotting validation of METTL3 protein expression after transfection of shRNAs and overexpression plasmids into CRC cells, respectively. (Q) After transfecting overexpression plasmid of METTL3 into SW480 cells, the mRNA expression of METTL3 and POU6F2-AS1 was verified by qRT–PCR. The data are presented as the mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001
Supplementary Material 6
: Additional file 1: Table S4. The sequences of primers used for qRT–PCR
Supplementary Material 7
: Additional file 1: Table S1. Relationship between POU6F2-AS1 expression and clinicopathological features in TMAs cohort
Supplementary Material 8
: Additional file 1: Table S2. Proteins identified by mass spectrometry
Supplementary Material 9
: Additional file 1: Table S3. The sequences of oligonucleotides and probes used in this study
Supplementary Material 10
: Additional file 2: Supplemental Materials and Methods
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jiang, T., Qi, J., Xue, Z. et al. The m6A modification mediated-lncRNA POU6F2-AS1 reprograms fatty acid metabolism and facilitates the growth of colorectal cancer via upregulation of FASN. Mol Cancer 23, 55 (2024). https://doi.org/10.1186/s12943-024-01962-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12943-024-01962-8