A reciprocal feedback between N6-methyladenosine reader YTHDF3 and lncRNA DICER1-AS1 promotes glycolysis of pancreatic cancer through inhibiting maturation of miR-5586-5p
Journal of Experimental & Clinical Cancer Research volume 41, Article number: 69 (2022)
Glycolysis is a pivotal process in metabolic reprogramming of tumorigenesis. Previous research has indicated that lncRNAs might play crucial roles in glycolysis of various tumors. However, the function of lncRNAs in glycolysis of pancreatic cancer has not been fully elucidated.
Bio-information analyses were applied to reveal the potential glycolysis-associated lncRNA. RT-PCR and fluorescence in situ hybridization (FISH) assays were applied to detect the expression of antisense RNA1 of DICER1 (DICER1-AS1) in pancreatic cancer tissues and cell lines. Gain- and loss-of-function experiments were performed to evaluate the roles of DICER1-AS1 in glycolysis and tumorigenesis of PC. Mechanistic experiments including luciferase reporter assay, RNA immunoprecipitation (RIP), and chromatin immunoprecipitation (ChIP) were employed to uncover the downstream targets and regulatory mechanism of DICER1-AS1 in glycolysis of PC.
Bio-information analysis indicated that DICER1-AS1 was downregulated in PC and negatively correlated with glycolytic genes expression. Meanwhile, overexpression of DICER1-AS1 inhibited glycolysis, proliferation, and metastasis of PC cells both in vitro and in vivo. Mechanistically, DICER1-AS1 promoted transcription of its sense gene DICER1 by recruiting transcriptional factor YY1 to the DICER1 promoter. Meanwhile, DICER1 promoted maturation of miR-5586-5p which consequently inhibited glycolytic gene expression including LDHA, HK2, PGK1, and SLC2A1. Notably, enhanced interaction between N6-methyladenosine (m6A) reader YTHDF3 and DICER1-AS1 led to degradation of DICER1-AS1 in response to glucose depletion. Moreover, our data revealed that YTHDF3 was a critical target for miR-5586-5p, by which forming a negative feedback with DICER1-AS1 to regulate glycolysis of PC.
Our results implicate a negative feedback of m6A reader YTHDF3 and glycolytic lncRNA DICER1-AS1 is involved in glycolysis and tumorigenesis of PC.
To maintain the urgent demand for tumorigenesis and aggressiveness, pancreatic cancer cells have extensively reprogrammed energy metabolism . One of the most common and important alterations of metabolism is aerobic glycolysis, which is termed the “Warburg effect”. Glycolysis displays a less efficient metabolism with more lactate production, lower extracellular pH value, and more glucose consumption, even under an oxygen-enriched microenvironment . Pancreatic cancer cells show significantly activated glycolysis, which was associated with poor prognosis in pancreatic cancer (PC) patients . Thereby, targeting the Warburg effect has been regarded as a promising therapeutic direction in tumors. Recent studies revealed that disrupting SLC2A1, LDHA or other key glycolytic enzymes resulted in remarkably reduced conversion of glucose to pyruvate and ATP, and elevated mitochondrial oxidative phosphorylation . Meanwhile, inhibitors of SLC2A1 or LDHA in combination with tumor chemotherapeutics displayed significant synergistic antitumor effect in vitro, which highlighting the promising role for cancer therapeutic strategy of glycolysis [5, 6]. Nevertheless, the metabolic process and molecular mechanisms of glycolysis in PC progression remain uncertain and need to be further determined.
Long noncoding RNAs (lncRNAs) are a type of special transcript longer than 200 nucleotides and lack protein-coding potential. Although lncRNAs were regarded as merely transcriptional ‘‘noise’’ before, emerging evidence has shown that lncRNAs play pivotal roles in transcriptional regulation, epigenetic gene regulation, and tumors . And our previous study has revealed that lncRNA GLS-AS, dysregulated under glucose and glutamine deprivation, repressed pancreatic cancer progression by repressing the Myc/GLS pathway . Specifically, emerging studies have shown the key role of lncRNAs in the regulation of glycolysis in cancers. For example, lncRNA MACC1-AS1 promotes glycolysis progression and anti-oxidative capabilities under metabolic stress in gastric cancer . Besides, lncRNA LINRIS has been reported as an independent prognostic factor that participated in the regulation of aerobic glycolysis of colorectal cancer . Similarly, Wang et al. have reported that lncRNA PVT1 facilitates PC glycolysis and progression by regulating the miR-519d-3p/HIF-1A pathway . Nevertheless, the molecular mechanisms of lncRNAs for the Warburg effect and its role in PC progression have not been fully discovered. Therefore, our present study sought to find a lncRNA that is critical for glycolysis and tumorigenesis of PC.
In this study, we identified that DICER1-AS1 was down-regulated in PC and was correlated with a poor prognosis. Meanwhile, both in vitro and in vivo experiments validated that DICER1-AS1 functioned as an inhibitor of PC in glycolysis and tumorigenesis by promoting DICER1-mediated maturation of miR-5586-5p. Notably, we further disclosed DICER1-AS1 was degraded in glucose deprivation and further explored the m6A-correlated regulation.
The cell lines (SW1990, AsPC-1, PANC-1, BxPC-3) of human pancreatic cancer were obtained from the American Type Culture Collection (ATCC, USA). The normal pancreatic ductal epithelial cell line (HPDE) was obtained from Bena Culture Collection (BNCC, China). Details are provided in Additional file 1.
Patients and clinical samples
The clinical samples including tumor tissues and matched noncancerous tissues were obtained from patients with pancreatic cancer at the Pancreatic Disease Institute of Union Hospital (Wuhan, China). Details are provided in Additional file 1.
RNA pull-down assay
For RNA pull-down assay, RNA was synthesized in vitro via RT-PCR using specific primers containing a T7 RNA-polymerase promoter sequence according to the manufacturer's instructions of MAXIscript® Kit (Ambion, USA). Then the RNA was end-labeled with desthiobiotin by using Pierce RNA 3′ End Desthiobiotinylation Kit (Thermo Scientific, USA). Finally, RNA pull-down assay was performed using the Pierce™ Magnetic RNA–Protein Pull-Down Kit (Thermo Scientific, USA) according to the manufacturer's instructions. Details are provided in Additional file 1.
m6A modifications of RNAs were measured by methylated RNA immunoprecipitation (MeRIP) assay. Details are provided in Additional file 1.
All statistical analyses were performed using R-3.0.2 software (http://cran.r-project.org/bin/windows/base/old/3.0.2/). Details are provided in Additional file 1.
Further applied methods
Additional RNA isolation, reverse transcription, and quantitative real-time PCR (qRT-PCR), cell transfection, RNA fluorescence in situ hybridization (RNA-FISH) assay, MTT assay, Wound healing assay, Transwell assay, Western blot assay, xenograft assay, Glucose uptake assay, Lactate Production Assay, Luciferase reporter assay, RNA stability assay, Chromatin Immunoprecipitation (ChIP) assay, RNA pull-down assay, RNA-binding protein immunoprecipitation (RIP), immunohistochemistry (IHC), Immunoflourescence, Gene set enrichment analysis (GSEA) are further described in the Additional file 1.
DICER1-AS1 is downregulated in pancreatic cancer tissues and negatively correlated with glycolysis pathway
To explore the potential glycolysis-correlated lncRNAs in pancreatic cancer, we first analyzed the public dataset of pancreatic adenocarcinoma (PAAD) derived from The Cancer Genome Atlas (TCGA) database . We identified 7 lncRNAs that were differentially expressed in PC patients based on the varied clinical progression, and negatively associated with the glycolysis pathway by the analysis from the gene set enrichment analysis (GSEA)  (Fig. 1A-B, Additional file 3: Figure S1A). Among these lncRNAs, antisense RNA of DICER1 (DICER1-AS1) was downregulated in PC tissues derived from GSE41368 dataset (Fig. 1C), then was specifically focused. The GSEA analysis demonstrated that DICER1-AS1 was negatively correlated with gene sets including glycolysis (Fig. 1D), cell metastasis, and pancreatic cancer (Additional file 3: Figure S1B). Gene ontology (GO) analyses revealed that DICER1-AS1 was correlated carbohydrate catabolic process and production of miRNA (Additional file 3: Figure S1C). Coincidentally, GSEA analysis based on GSE41368 also revealed that DICER1-AS1 was negatively correlated with glycolysis process (Additional file 3: Figure S1D). The TCGA database revealed that DICER1-AS1 was downregulated in PC with advanced clinical stages and tumor grade, and correlated with longer overall survival (OS), disease-free interval (DFI), progression-free interval (PFI), and disease-specific survival (DSS) (Fig. 1F, Additional file 3: Figure S1A). Meanwhile, DICER1-AS1 level was also downregulated in PC tissues derived from cohort 1 (Additional file 2: Table S3), which was composed of 86 paired PC tissues and normal tissues, as well as correlated with longer overall survival (Additional file 3: Figure S1E-F). Furthermore, receiver operating characteristic (ROC) curves showed that DICER1-AS1 can identify PC patients from healthy individuals both in the GEO database (GSE41368, GSE16515) and cohort 1 (Additional file 3: Figure S1G). Subsequently, we analyzed the correlation of DICER1-AS1 expression with different clinicopathological features in our cohort1 and found that the DICER1-AS1 level was significantly associated with tumor size and TNM stage (Additional file 3: Figure S1H). Univariate and multivariate regression analyses showed that DICER1-AS1 expression was an independent predictor of PC patients with significant hazard ratios for predicting clinical outcomes (Additional file 2: Table S4).
Meanwhile, FISH analysis confirmed that DICER1-AS1 expression was lower in PC tissues compared with normal pancreatic (NP) tissues (Fig. 1G). Besides, DICER1-AS1 expression was significantly lower in different pancreatic cell lines than normal human pancreatic duct epithelial cells (HPDE) (Fig. 1H). For a good efficiency of RNA interference experiments, further study was investigated in PC cells (PANC-1, BxPC-3) with relatively high expression of DICER1-AS1. Subcellular fractionation assays showed that DICER1-AS1 is located both in the nucleus and cytoplasm (Fig. 1H). Furthermore, the online coding potential tool (CPC and CPAT) [14, 15] confirmed that DICER1-AS1 has a very low protein-coding ability (Additional file 3: Figure S2A). Besides, ribosome profiling data showed a remarkably low-level peak in the DICER1-AS1 gene indicating a low protein-coding ability of DICER1-AS1  (Additional file 3: Figure S2B). Meanwhile, we found that DICER1-AS1 has low evolutionary conservation by using PhyloP analysis across 100 vertebrates (Additional file 3: Figure S2C). After mining public datasets (32 studies) of cBioPortal for Cancer Genomics, we found that there was no amplification and deep deletion of the DICER1-AS1 gene in pancreatic cancer (Additional file 3: Figure S2D). Taken together, these results implied that DICER1-AS1 was a glycolysis-associated lncRNA that might function as an inhibitor in PC progression.
DICER1-AS1 inhibits glycolysis, proliferation and metastasis of PC cells
To further investigate the biological impact of DICER1-AS1 in pancreatic cancer cells, DICER1-AS1 was manipulated with overexpressed plasmid or siRNA (siDICER1-AS1#1–3), respectively (Fig. 2A, Additional file 3: Figure S3A). We found that the overexpression of DICER1-AS1 significantly inhibited (Fig. 2B-D), while the knockdown of DICER1-AS1 obviously increased the proliferation, migration, and invasion of BxPC-3 and PANC-1 cells (Figure S3B-D). We next identified whether altered DICER1-AS1 levels directly influence glycolytic metabolism in pancreatic cancer cells. Indeed, overexpression of DICER1-AS1 significantly decreased, while DICER1-AS1 knockdown increased the extracellular acidification rate (ECAR), glucose consumption, lactate production, and ATP level. On the contrary, oxygen consumption rate (OCR) was notably increased by DICER1-AS1 overexpression but decreased by DICER1-AS1 knockdown (Fig. 2E-I, Additional file 3: Figure S3E-I). In addition, the 2-DG, an inhibitor of glycolysis, significantly blocked the glucose uptake and lactate production induced by DICER-AS1 knockdown (Additional file 3: Figure S3J-K). These data indicate that DICER1-AS1 is an inhibitor for pancreatic cancer by suppressing glycolysis.
DICER1 is a critical target for DICER1-AS1 regulating glycolysis
Since DICER1-AS1 is an antisense RNA of DICER1 (Fig. 3A), we wonder whether DICER1 might be a critical target of DICER1-AS1. Coincidently, DICER1-AS1 was positively correlated with DICER1 expression in pancreatic cancer (Fig. 3B), as well as in a variety of cancers or normal tissues (Additional file 3: Figure S4A and 4B). Moreover, TCGA and GEO databases displayed that low expression of DICER1 correlated with a poor overall survival both in pancreatic cancer (Fig. 3C) and a variety of malignant tumors (Additional file 3: Figure S4C). Meanwhile, the copy number of the DICER1 gene was neither significantly altered in the different clinical stages nor associated with the survival of PC patients from the TCGA database (Additional file 3: Figure S5A-B), which indicates that the mRNA level but not copy number of DICER1 gene is decisive for pancreatic cancer.
Coincidently, GO and Kyoto Encyclopedia of Genes and Genomes (KEGG)  analyses further indicated that DICER1 was related to glucose homeostasis and glycolysis (Fig. 3D-E). Meanwhile, overexpression or knockdown of DICER1-AS1 could increase or decrease DICER1 expression of PC cells in both mRNA and protein level (Fig. 3F, Additional file 3: Figure S5C-D), which was further validated with RNA-FISH and immunofluorescence assays (Fig. 3G, Additional file 3: Figure S5E). After being treated with transcription inhibitor actinomycin D, the decay of DICER1 was not affected by overexpression or depletion of DICER1-AS1 both in BxPC-3 and PANC-1 cells (Fig. 3H, Figure S5F), which indicated that DICER1 would be regulated by DICER1-AS1 in transcription but not post-transcription level. Additionally, DNA fragments containing the DICER1 promoter were inserted into the promoter region of firefly luciferase reporter plasmids to perform the luciferase assay. The results showed that overexpression of DICER1-AS1 enhanced, while depletion of DICER1-AS1 decreased the transcriptional activity of the DICER1 promoter (Fig. 3I, Figure S5G). These results further confirmed that DICER1-AS1 regulated DICER1 expression at the transcriptional level.
To further investigate the importance of DICER1 in the regulation of glycolysis, DICER1-AS1/DICER1 pathway was modulated by rescue or reverse experiments. Results showed that overexpression or knockdown of DICER1-AS1 induced or inhibited the levels of DICER1, which could be reversed by inhibition or overexpression of DICER1 (Fig. 3J, Figure S5H). Moreover, depletion of DICER1 obviously rescued the glucose uptake and lactate production of PC cells which were inhibited by DICER1-AS1 overexpression (Fig. 3K-L). In contrast, the glucose uptake and lactate production promoted by DICER1-AS1 knockdown was remarkably reversed by DICER1 overexpression (Additional file 3: Figure S5I-J). Collectively, our data intensively implied that DICER1-AS1 regulates the glycolysis of PC cells by regulating DICER1 expression.
Transcription factor YIN-YANG-1 (YY1) is a DICER1-AS1 binding protein
To explore the underlying regulatory role of DICER1-AS1 on DICER1, we performed an intersection analysis of transcription factors which were correlated with DICER1 expression (P < 0.01, r > 0) in TCGA and GEO database. Further over-lapping analysis with protein-RNA pairs estimating database catRAPID  and transcription factor binding database JASPAR (https://jaspar.genereg.net/) revealed that transcription factor YY1 might participate in DICER1-AS1/DICER1 pathway (Fig. 4A). As shown in Fig. 4B, the catRAPID algorithm showed an intensive interaction between DICER1-AS1 and YY1. Meanwhile, analysis of RNA-seq data from the GEO, TCGA, and CCLE databases demonstrated that DICER1 was significantly correlated with YY1 mRNA expression both in cancer tissue and tumor cell lines (Fig. 4C, Additional file 3: Figure S6A). Moreover, low expression of YY1 was associated with poor outcomes of PC patients and other types of cancers derived from the TCGA and GEO database (Fig. 4D, Additional file 3: Figure S6B).
To verify the interaction between DICER1-AS1 and YY1, RIP assay was applied by using the anti-YY1 antibody and IgG antibody. The RIP assay confirmed a specific binding between DICER1-AS1 to YY1 both in BxPC-3 and PANC-1 cells (Fig. 4E). Furthermore, biotin-labeled RNA pull-down and western blot assays further validated the interaction of YY1 with DICER1-AS1 (Fig. 4F). Besides, deletion-mapping analyses indicated that exon 4 (374-830nt) of DICER1-AS1 was essential for its interaction with YY1 protein (Fig. 4G). Coincidently, overexpression of full or 374-830nt of DICER1-AS1, but not depletion of 374-830nt and antisense, could increase the expression of DICER1 in BxPC-3 and PANC-1 cells (Fig. 4H). To determine which functional domain of DICER1 is essential for the interaction with DICER1-AS1, RIP assay was performed by using Flag-tagged full-length and truncated mutants of YY1 protein. The results showed that the zinc-finger (295–414 aa) domain of Flag-tagged DICER1 protein was crucial for its interaction with DICER1-AS1 (Fig. 4I). The TASSER and SWISS-MODEL software [19, 20] were used to predict the 3D structure of the full-length and the zinc-finger domain of YY1 protein, respectively. And the RNA secondary structure of DICER1-AS1 was predicted via the RNA fold web server. Then the visual structure of the YY1/DICER1-AS1 binding complex was presented in Fig. 4J. Collectively, these results indicated that YY1 is a DICER1-AS1 binding protein which might be essential for DICER1 expression.
DICER1-AS1 promotes DICER1 transcription by recruiting YY1 to the DICER1 promoter
According to the analysis from the JASPAR database, the promoter area of the DICER1 gene possesses two potential binding positions for transcription factor YY1 (Fig. 5A). Importantly, the Cistrome Data Browser database  showed a significant enhancement of YY1 at the DICER1 promoter region accompanied by enriched peaks of histone modifications (H3K4me3 and H3K27ac) which indicating an active transcription (Fig. 5B). Furthermore, analysis of histone transcriptional markers of activation revealed that DICER1-AS1 overexpression increased, while DICER1-AS1 knockdown decreased the interaction of H3K4me3 and H3K27ac with DICER1 promoter (Additional file 3: Figure S7A-B). These results strongly suggested that DICER1-AS1 may have a potential role in regulating DICER1 transcription. Moreover, ChIP analysis confirmed that YY1 was accumulated in the P2 site of the DICER1 promoter (Fig. 5C). To verify the influence of YY1 on DICER1 transcription, we detected the luciferase activities of DICER1 promoter vectors containing the wild-type (WT) or mutant (MUT1-3) binding sites (Fig. 5D). The luciferase assays demonstrated that YY1 overexpression or knockdown could increase or decrease the luciferase intensity in PC cells transfected with WT plasmid and MUT1 plasmid, which further confirmed that the P2 site of the DICER1 promoter is necessary for YY1 to interact with DICER1 promoter (Fig. 5E). Coincidently, overexpression of YY1 increased, while silencing of YY1 decreased the expression of DICER1 (Fig. 5F). Furthermore, silencing of YY1 reversed the effects of DICER1-AS1 overexpression, while overexpression of YY1 rescued the effects of DICER1-AS1 knockdown on DICER1 promoter activity, YY1 accumulation on DICER1 promoter, DICER1 expression, as well as glucose uptake and lactate production in PC cells (Fig. 5G-K, Figure S7C-H). These data indicated that DICER1-AS1 promotes DICER1 transcription by recruiting YY1 to the DICER1 promoter.
miR-5586-5p is an essential target for DICER1-AS1/DICER1 pathway
Since DICER1 plays an important role in the biogenesis of miRNAs, we presumed whether specific miRNAs might be contributing to the DICER1-AS1/DICER1 pathway. To determine the putative target miRNA, we analyzed the miRNAs positively correlated with expression of DICER1-AS1 and DICER1, as well as overall survival of PC in the TCGA database (Fig. 6A). The overlapping analyses indicated that miR-5586-5p, miR-29c-5p, miR-1301-3p, and miR-1224-5p were positively correlated with both DICER1-AS1 and DICER1 (Fig. 6B-C). Additionally, these 4 miRNAs were associated with favorable overall survival and relapse-free survival (Additional file 3: Figure S8A-B), and negatively associated with advanced clinical stages and grades (Fig. 6D, Additional file 3: Figure S8C-E). Although without obvious alteration of pre-miRNAs of the four miRNAs (Additional file 3: Figure S8F-G), knockdown or overexpression of DICER1 significantly decreased or increased expression of the 4 miRNAs in PC cells (Fig. 6E, Additional file 3: Figure S8H). Among those miRNAs, miR-5586-5p displayed the most significant change, then we selected miR-5586-5p as the critical target for further investigation. Consistently, RIP assays displayed that knockdown of DICER1 significantly decreased, but overexpression of DICER1 increased the binding of DICER1 to pre-mir-5586 (Fig. 6F-G). Moreover, overexpression or knockdown of DICER1 prevented the alteration of miR-5586-5p expression induced by the knockdown or overexpression of DICER1-AS1 (Fig. 6H). And miR-5586-5p was downregulated in PC tissues compared with normal tissues. (Additional file 3: Figure S8I). Meanwhile, miR-5586-5p mimics or inhibitors successfully reversed the changes of glucose uptake and lactate production in PC cells which were induced by knockdown or overexpression of DICER1-AS1 (Fig. 6I-J, Additional file 3: Figure S8J). Subsequently, analysis of TCGA database showed that the miR-5586-5p level was negatively associated with the tumor stage (Additional file 2: Table S5). Univariate and multivariate regression analyses showed that miR-5586-5p expression was an independent predictor of PC patients with significant hazard ratios for predicting favorably clinical outcomes (Additional file 2: Table S6). Collectively, these data demonstrated that miR-5586-5p is a critical target of DICER1-AS1/DICER1 pathway and functions as an inhibitor in pancreatic cancer.
MiR-5586-5p regulates glycolysis by inducing mRNA degradation of glycolytic genes
Since GSEA analysis revealed that miR-5586-5p was negatively correlated with the gene set of pancreatic cancer and glycolysis (Fig. 7A), we further investigated whether miR-5586-5p can regulate glycolytic gene expression. To test this possibility, we first performed the bioinformatics analysis by using the TargetScan database and then performed overlapping analyses with TCGA and GEO databases (GSE16515, GSE41368). The analysis demonstrated 32 potential targets of miR-5586-5p were overexpressed in PC tissues and correlated with poor survival (Fig. 7B), which including 4 glycolytic genes SLC2A1, LDHA, HK2, and PGK1 (Additional file 3: Figure S9A-C). The miR-5586-5p binding sites of the 4 glycolytic genes were predicted by the TargetScan database (Fig. 7C). As shown in Fig. 7D, four glycolytic genes (SLC2A1, LDHA, HK2, PGK1) were negatively correlated with miR-5586-5p.
Furthermore, miR-5586-5p mimics substantially reduced, whereas miR-5586-5p inhibitors increased the mRNA and protein levels of SLC2A1, LDHA, HK2, and PGK1 (Fig. 7E-F, Additional file 3: Figure S9D-E). Since AGO2 participates in the formation of the RISC complex, we further investigated whether miR-5586-5p affects the binding of AGO2 with target mRNAs. RIP assays indicated that mimics or inhibitors of miR-5586-5p significantly facilitated or repressed the interaction of AGO2 with SLC2A1, LDHA, HK2, and PGK1 (Fig. 7G). Meanwhile, we found that mimics or inhibitors of miR-5586-5p significantly increased or reduced the degradation of SLC2A1, SLC2A1, HK2, and PGK1 mRNAs in BxPC-3 cells treated with actinomycin D (Additional file 3: Figure S9F). To further verify these glycolytic genes are the direct target of miR-5586-5p, we constructed luciferase reporter plasmids containing wild type or mutant type 3’UTR (WT, MUT) of these genes (Fig. 7H). Luciferase assays indicated that miR-5586-5p mimics decreased, but the inhibitors increased the reporter activities, which was further reversed by knockdown or overexpression of DICER1-AS1 (Fig. 7I). Moreover, qRT-PCR and western blot assays confirmed that knockdown of DICER1-AS1 relieved the inhibition of miR-5586-5p mimics on the glycolytic genes, whereas overexpression of DICER1-AS1 deprived the upregulation of miR-5586-5p inhibitors on the glycolytic genes (Fig. 7J-K). Collectively, these results demonstrated that miR-5586-5p regulates the stability of the mRNAs of glycolytic genes to participate in the DICER1-AS1/DICER1 pathway.
N6-methyladenosine reader YTHDF3 decreases DICER1-AS1 stability of PC cells in response to glucose deprivation
Numerous lncRNAs were aberrantly expressed for adaption to nutrient stress in tumor microenvironment, including hypoxia, acidity, glutamine deprivation, and glucose deprivation [22, 23]. Interestingly, DICER1-AS1 was only notably decreased during glucose deprivation, but without significant alteration under hypoxia, acidity, and glutamine deprivation in PC cells (Fig. 8A, Additional file 3: Figure S10A). To further investigate whether the glucose deprivation induced downregulation of DICER1-AS1 of PC cells, the luciferase reporter assay was performed. Basic pGL3 plasmid and pGL3 plasmid containing DICER1-AS1 promoter was transfected into BxPC-3 and PANC-1 cells. The results revealed that the luciferase intensity was not changed in both PANC-1 and BxPC-3 cells, indicating that DICER1-AS1 might be regulated at the post-transcriptional level during glucose deprivation (Fig. 8B). Nevertheless, accelerated RNA decay of DICER1-AS1 upon glucose starvation was confirmed in PC cells (Fig. 8C).
Since research had referred that m6A modification contributes to degradation of lncRNAs , we then wonder whether the enhanced degradation of DICER1-AS1 during glucose deprivation was mediated by m6A modification. According to the online tool m6A-Atlas, a high-confidence collection of reliable m6A sites , DICER1-AS1 has m6A sites in the exon region (at Chr14: 95,646,121, 95,646,197). Meanwhile, bioinformatics database m6A2Target (http://m6a2target.canceromics.org/#/) also demonstrated that DICER1-AS1 interacted with a m6A reader YTH-domain family member 3 (YTHDF3) (Fig. 8D). Meanwhile, analysis of TCGA database showed that YTHDF3 was significantly increased in pancreatic cancer tissues, as well as negatively correlated with DICER1-AS1 expression and overall survival (Fig. 8E-G). Then we found that knockdown or overexpression of YTHDF3 increased or decreased the expression of DICER1-AS1 (Fig. 8H, Additional file 3: Figure S10B). Besides, overexpression or knockdown of YTHDF3 increased or decreased the mRNA and protein levels of glycolytic genes (SLC2A1, LDHA, HK2, PGK1). Meanwhile, the expression of DICER1 was reduced or increased by YTHDF3 overexpression and knockdown, respectively (Additional file 3: Figure S10C-D).
To validate the results of m6A databases, we detected the m6A modification of DICER1-AS1 in PC cells by the MeRIP-qPCR assay. The results showed that DICER1-AS1 was enriched with m6A antibody in both PANC-1 and BxPC-3 cells (Fig. 8I). RIP assay with YTHDF3 antibody further verified that YTHDF3 bonded to DICER1-AS1 to form a m6A modification complex (Fig. 8J). In addition, with increased interaction between YTHDF3 and DICER1-AS1 induced by YTHDF3 overexpression, YTHDF3 overexpression decreased m6A enrichment and stability of DICER1-AS1. On the contrary, the m6A enrichment and stability of DICER1-AS1 were enhanced by knockdown of YTHDF3, accompanied by decreased interaction between YTHDF3 and DICER1-AS1 (Fig. 8K-M, Additional file 3: Figure S10E-G). Coincidently, the RIP and MeRIP assay further revealed that glucose deprivation notably increased interaction between YTHDF3 and DICER1-AS1, which consequently decreased m6A enriched DICER1-AS1 (Fig. 8N-O, Additional file 3: Figure S10H-I). Furthermore, the glucose deprivation-enhanced interaction between DICER1-AS1 and YTHDF3 was significantly inhibited by YTHDF3 knockdown, while the decreased stability of DICER1-AS1 during glucose deprivation was obviously increased after YTHDF3 depletion (Fig. 8P, Additional file 3: Figure S10J-K). Coincidently, MeRIP assay showed the inhibition of the m6A enrichment of DICER1-AS1 under glucose deprivation was reversed by YTHDF3 knockdown (Additional file 3: Figure S10L). Taken together, our data revealed that YTHDF3 decreases DICER1-AS1 stability via a m6A dependent manner in response to energy stress.
YHDF3 is a reciprocal target for DICER1-AS1 via miR-5586-5p-induced inhibition
Furthermore, we detected whether YTHDF3 expression was affected in response to glucose deprivation. Interestingly, both YTHDF3 mRNA and protein expression were increased during glucose deprivation (Additional file 3: Figure S11A). Meanwhile, we found that miR-5586-5p was downregulated in BxPC-3 and PANC-1 cells under glucose deprivation (Additional file 3: Figure S11B). Moreover, TargetScan predicted a binding site for miR-5586-5p in 3’UTR of YTHDF3, which suggested that YTHDF3 would be mediated by miR-5586-5p. And the reporter plasmid harboring the wild type (WT) or mutant type (MUT) 3’UTR region of YTHDF3 was constructed and transfected into BxPC-3 cells to perform luciferase reporter assay (Additional file 3: Figure S11C). The results showed that miR-5586-5p mimics or inhibitors significantly decreased or increased the luciferase reporter activity in WT cells, but without alteration in MUT cells (Additional file 3: Figure S11D-E). Moreover, miR-5586-5p mimics or inhibitors significantly reversed or enhanced the YTHDF3 expression under glucose deprivation (Additional file 3: Figure S11F-G). The luciferase reporter assay showed miR-mimics reversed the decrease of luciferase density of WT cells which was affected in response to glucose deprivation (Additional file 3: Figure S11H). Therefore, these results implicated that YTHDF3 is a critical target of miR-5586-5p.
To explore the role of YTHDF3 and miR-5586-5p on the DICER1-AS1 expression during glucose deprivation, PC cells were transfected with vector or YTHDF3, and those co-transfected with miR-5586-5p mimics or inhibitors. The results showed that YTHDF3-mediated inhibition of DICER1-AS1 level was neutralized by miR-5586-5p mimics (Additional file 3: Figure S11I). Furthermore, YTHDF3 knockdown increased the DICER1-AS1 level, which was reversed by miR-5586-5p inhibitors (Additional file 3: Figure S11J). Together, these data indicate a reciprocal feedback of YTHDF3 and DICER1A-S1/miR-5586-5p in PC cells in response to glucose deprivation.
YTHDF3/DICER1-AS1/DICER1/miR-5586-5p axis is pivotal for glycolysis and tumorigenesis of PC
To further validate the function of DICER1-AS1 and its pathway in glycolysis of PC cells, we next performed immunohistochemical staining in 86 PC patients’ tissues of cohort 1. Consistently, the samples with low expression of DICER1-AS1 showed strong staining for YTHDF3 and glycolytic genes (SLC2A1, LDHA, HK2, PGK1), as well as weak staining of DCIER1. Whereas samples with high expression of DICER1-AS1 showed strong staining for DICER1, but weak staining for YTHDF3 and glycolytic genes (Fig. 9A-B). Kaplan–Meier survival analyses from TCGA database also indicated that low expression of DICER1-AS1/DICER1, DICER1/miR-5586-5p, and DICER1-AS1/miR-5586-5p groups showed a poor overall survival of PC patients. Meanwhile, low expression of YTHDF3 accompanied by high expression of miR-5586-5p, DICER1-AS1, or DICER1 showed a long overall survival of PC patients (Fig. 9C). To further evaluate the biological function of DICER1-AS1 in vivo, BxPC-3 cells were transfected stably by using lentivirus containing negative control (LV-NC) or DICER1-AS1 sequences (LV-DICER1-AS1). These transfected cells were further implanted into nude mice to construct a subcutaneous xenograft tumor model. The results revealed that the stable expression of DICER1-AS1 into BxPC-3 cells resulted in retarded growth and lower tumor weight of subcutaneous xenograft tumors (Fig. 9D-E). To evaluate the tumor metastasis in vivo, the nude mice were received tail vein injections with transfected BxPC-3 cells. We found that the number of lung metastases in the LV-DICER1-AS1 group was decreased compared to the control group (Fig. 9F). Correspondingly, the DICER1-AS1 transfected tumors displayed upregulation of DICER1-AS1, DICER1, and miR-5586-5p, as well as downregulation of YTHDF3, SLC2A1, LDHA, HK2, and PGK1 (Additional file 3: Figure S12A-D). Collectively, these results suggested that YTHDF3/DICER1-AS1/DICER1/miR-5586-5p axis is pivotal for glycolysis and tumorigenesis of pancreatic cancer.
Recently, the non-coding RNAs (ncRNAs), including miRNAs and lncRNAs have received growing attention in human tumors , particularly in pancreatic cancer [8, 27, 28]. Although studies have evaluated the involvement of lncRNAs in the dysregulation of glycolysis metabolism in pancreatic cancers, the precise mechanism needs to be further investigated [29, 30]. In this study, we identify that DICER1-AS1 is an inhibitor for glycolysis and tumorigenesis of pancreatic cancer by inducing miRNA-mediated mRNA degradation of glycolytic genes. Moreover, the present study indicates that DICER1-AS1 is epigenetically downregulated by its target m6A reader YTHDF3, which forms a negative feedback.
The bioinformatic analyses showed that DICER1-AS1 was downregulated in PC tissues and negatively correlated with glycolytic pathway and overall survival. We further identified that overexpression of DICER1-AS1 significantly inhibited glycolysis and consequently prohibited proliferation and metastasis of pancreatic cancer cells both in vitro and in vivo. Moreover, we further displayed that DICER1-AS1 was downregulated in PC tissues in our cohort 1 and associated with poor overall survival of PC patients. Thus, our study implicates that DICER1-AS1 is an inhibitor for glycolysis and tumorigenesis of pancreatic cancer. Differently, research from Gu et al. reported that DICER1-AS1 promoted the proliferation, invasion, and autophagy of osteosarcoma cells via miR-30b/ATG5 axis . Other research also indicated that DICER1-AS1 functioned as promoter in colorectal cancer by regulating miR-296-5p/STAT3 pathway . Therefore, these researches implicate that the function of DICER1-AS1 in heterogeneous cancer, which needs further investigation.
Natural antisense transcripts are a cluster of RNAs that are transcribed from the opposite strand of a sense DNA strand. It’s reported that antisense RNA can alter the expression of sense genes at the transcriptional or post-transcriptional level [33, 34]. Similarly, our present study demonstrated the inhibition on glycolysis of DICER1-AS1 was dependent on regulating DICER1 transcription. Meanwhile, as a critical processor for maturation of miRNA, we further founded that DICER1 promoted maturation of miR-5586-5p. Research also showed that knockdown of DICER1 significantly inhibited xenograft tumor growth of BxPC-3 cells by suppressing expression of miR-22, miR-143, let-7i, and miR-29b . PICK1 inhibits the processing of pre-mir-615-3p to mature miR-615-3p via interfering with the binding of DICER1 to Smad2/3 in breast cancer cells . Thereby, our research provides a new clue for interaction between lncRNA and miRNA which involving DICER1-mediated miRNA maturation.
YY1 is a ubiquitously distributed transcription factor belonging to the GLI-Krüppel family of zinc-finger protein, associated with multiple cellular processes . YY1 suppressed the proliferation via inhibiting the transcription of lncRNA SOX2OT, thereby decreased the expression of its target gene SOX2 in PC . Other studies also showed that YY1 suppressed the proliferation, invasion, and metastatic properties of PC cells by downregulating MMP10 , as well as FER and MMP2 . Similarly, our present study validated that DICER1-AS1 promoted DCIER1 transcription by recruiting YY1 to promoter of DICER1. Meanwhile, YY1 was downregulated in PC tissues and correlated with DICER1 expression and poor prognosis of PC patients. Therefore, these data also identify that YY1 functions as an inhibitor in PC by interacting with DICER1-AS1/DICER1 pathway.
It’s reported that miRNAs play a vital role in metabolic homeostasis including glycolysis in tumors. For example, circMAT2B upregulated PKM2 by functioning as a “sponge” of miR-338-3p, meanwhile promotes the glycolysis level and hepatocellular carcinoma progression . Besides, cell-secreted miR-122 is associated with metastasis in breast cancer patients by inhibiting glycolysis-related pyruvate kinase . The present study demonstrated that miR-5586-5p induced degradation of glycolytic genes including SLC2A1, LDHA, HK2, and PGK1 by binding to 3’UTR of each mRNA, which consequently inhibited glycolysis of PC cells. Meanwhile, GEO database displayed that the 4 glycolytic genes were highly expressed in PC tissues and negatively correlated with miR-5586-5p and overall survival of PC patients. Similarly, research from Bai et al. showed that a seven-miRNA expression-based prognostic signature including miR-5586-5p was downregulated in PC which might stratify patients with PC into low- and high-risk groups . Thereby, these data intensively reminder that miR-5586 is an inhibitor of PC and provides further evidence that miRNAs play critical roles in glycolysis of tumor.
Recent studies showed that m6A modification, as the most abundant RNA modification, was responsible for the dysregulation of lncRNAs . Previous studies showed that the stability of lncRNA GAS5 was destroyed by YTHDF3 and consequently promoted progression of colorectal cancer . Similarly, our study showed glucose deprivation significantly induced degradation of DICER1-AS1, which was mediated by m6A reader YTHDF3. Interestingly, our research also displayed that YTHDF3 was a critical target for miR-5586-5p. To explore the comprehensive effects of YTHDF3/DICER1-AS1/DICER1/miR-5586-5p axis on the glycolysis and prognosis of PC, we performed IHC assay and Kaplan–Meier survival analyses. The results demonstrated that high expression of DICER1-AS1 correlated with high expression of DCIER1, but correlated with low expression of YTHDF3 and the four glycolytic genes. Moreover, high expression of DICER1-AS1 along with high expression of DICER1 or miR-5586 was associated with longer overall survival. On the contrary, high expression of YTHDF3 along with low expression of DICER1 or DICER1-AS1 demonstrated shorter overall survival. Moreover, the in vivo experiments further identified that overexpression of DICER1-AS1 significantly inhibited the growth and metastasis of PC tumors.
In summary, these results revealed an m6A-modified DICER1-AS1 is involved in glycolysis of PC by regulating glycolytic genes expression through processing maturation of miR-5586-5p and implicates that DICER1-AS1 might be a favorable biomarker and therapeutic target for treatment of PC.
Availability of data and materials
All data generated or analyzed during this study are included either in this article or in the supplementary Materials and Methods, Tables, and Figures.
Noncancerous peritumoral tissues
Long non-coding RNAs
Fluorescence in situ hybridization
Methylated RNA immunoprecipitation
Quantitative real-time PCR
American Type Culture Collection
National Comprehensive Cancer Network
Human pancreatic duct epithelial
The Cancer Genome Atlas
Kyoto Encyclopedia of Genes and Genomes
Gene set enrichment analysis
Coding-Potential Assessment Tool
Coding Potential Calculator
Cancer Cell Line Encyclopedia
Receiver operating characteristic curve
Halbrook CJ, Lyssiotis CA. Employing Metabolism to Improve the Diagnosis and Treatment of Pancreatic Cancer. Cancer Cell. 2017;31:5–19.
Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324:1029–33.
Yang J, Ren B, Yang G, Wang H, Chen G, You L, Zhang T, Zhao Y. The enhancement of glycolysis regulates pancreatic cancer metastasis. Cell Mol Life Sci. 2020;77:305–21.
Ooi AT, Gomperts BN. Molecular Pathways: Targeting Cellular Energy Metabolism in Cancer via Inhibition of SLC2A1 and LDHA. Clin Cancer Res. 2015;21:2440–4.
Liu Y, Cao Y, Zhang W, Bergmeier S, Qian Y, Akbar H, Colvin R, Ding J, Tong L, Wu S, et al. A small-molecule inhibitor of glucose transporter 1 downregulates glycolysis, induces cell-cycle arrest, and inhibits cancer cell growth in vitro and in vivo. Mol Cancer Ther. 2012;11:1672–82.
Zhao Y, Liu H, Liu Z, Ding Y, Ledoux SP, Wilson GL, Voellmy R, Lin Y, Lin W, Nahta R, et al. Overcoming trastuzumab resistance in breast cancer by targeting dysregulated glucose metabolism. Cancer Res. 2011;71:4585–97.
Ponting CP, Oliver PL, Reik W. Evolution and functions of long noncoding RNAs. Cell. 2009;136:629–41.
Deng SJ, Chen HY, Zeng Z, Deng S, Zhu S, Ye Z, He C, Liu ML, Huang K, Zhong JX, et al. Nutrient Stress-Dysregulated Antisense lncRNA GLS-AS Impairs GLS-Mediated Metabolism and Represses Pancreatic Cancer Progression. Cancer Res. 2019;79:1398–412.
Zhao Y, Liu Y, Lin L, Huang Q, He W, Zhang S, Dong S, Wen Z, Rao J, Liao W, Shi M. The lncRNA MACC1-AS1 promotes gastric cancer cell metabolic plasticity via AMPK/Lin28 mediated mRNA stability of MACC1. Mol Cancer. 2018;17:69.
Wang Y, Lu JH, Wu QN, Jin Y, Wang DS, Chen YX, Liu J, Luo XJ, Meng Q, Pu HY, et al. LncRNA LINRIS stabilizes IGF2BP2 and promotes the aerobic glycolysis in colorectal cancer. Mol Cancer. 2019;18:174.
Sun J, Zhang P, Yin T, Zhang F, Wang W. Upregulation of LncRNA PVT1 Facilitates Pancreatic Ductal Adenocarcinoma Cell Progression and Glycolysis by Regulating MiR-519d-3p and HIF-1A. J Cancer. 2020;11:2572–9.
Integrated Genomic Characterization of Pancreatic Ductal Adenocarcinoma. Cancer Cell. 2017;32:185-203.e113.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, Gao G. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007;35:W345-349.
Wang L, Park HJ, Dasari S, Wang S, Kocher JP, Li W. CPAT: Coding-Potential Assessment Tool using an alignment-free logistic regression model. Nucleic Acids Res. 2013;41:e74.
Michel AM, Fox G. A MK, De Bo C, O’Connor PB, Heaphy SM, Mullan JP, Donohue CA, Higgins DG, Baranov PV: GWIPS-viz: development of a ribo-seq genome browser. Nucleic Acids Res. 2014;42:D859-864.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, Benner C, Chanda SK. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.
Livi CM, Klus P, DelliPonti R, Tartaglia GG. catRAPID signature: identification of ribonucleoproteins and RNA-binding regions. Bioinformatics. 2016;32:773–5.
Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, Heer FT, de Beer TAP, Rempfer C, Bordoli L, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46:W296-w303.
Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat Methods. 2015;12:7–8.
Zheng R, Wan C, Mei S, Qin Q, Wu Q, Sun H, Chen CH, Brown M, Zhang X, Meyer CA, Liu XS. Cistrome Data Browser: expanded datasets and new tools for gene regulatory analysis. Nucleic Acids Res. 2019;47:D729-d735.
Denko NC. Hypoxia, HIF1 and glucose metabolism in the solid tumour. Nat Rev Cancer. 2008;8:705–13.
Pan M, Reid MA, Lowman XH, Kulkarni RP, Tran TQ, Liu X, Yang Y, Hernandez-Davies JE, Rosales KK, Li H, et al. Regional glutamine deficiency in tumours promotes dedifferentiation through inhibition of histone demethylation. Nat Cell Biol. 2016;18:1090–101.
Ni W, Yao S, Zhou Y, Liu Y, Huang P, Zhou A, Liu J, Che L, Li J. Long noncoding RNA GAS5 inhibits progression of colorectal cancer by interacting with and triggering YAP phosphorylation and degradation and is negatively regulated by the m(6)A reader YTHDF3. Mol Cancer. 2019;18:143.
Tang Y, Chen K, Song B, Ma J, Wu X, Xu Q, Wei Z, Su J, Liu G, Rong R, et al. m6A-Atlas: a comprehensive knowledgebase for unraveling the N6-methyladenosine (m6A) epitranscriptome. Nucleic Acids Res. 2021;49:D134-d143.
Anastasiadou E, Jacob LS, Slack FJ. Non-coding RNA networks in cancer. Nat Rev Cancer. 2018;18:5–18.
Zhou C, Yi C, Yi Y, Qin W, Yan Y, Dong X, Zhang X, Huang Y, Zhang R, Wei J, et al. LncRNA PVT1 promotes gemcitabine resistance of pancreatic cancer via activating Wnt/β-catenin and autophagy pathway through modulating the miR-619-5p/Pygo2 and miR-619-5p/ATG14 axes. Mol Cancer. 2020;19:118.
Franses JW, Philipp J, Missios P, Bhan I, Liu A, Yashaswini C, Tai E, Zhu H, Ligorio M, Nicholson B, et al. Pancreatic circulating tumor cell profiling identifies LIN28B as a metastasis driver and drug target. Nat Commun. 2020;11:3303.
Zhai S, Xu Z, Xie J, Zhang J, Wang X, Peng C, Li H, Chen H, Shen B, Deng X. Epigenetic silencing of LncRNA LINC00261 promotes c-myc-mediated aerobic glycolysis by regulating miR-222-3p/HIPK2/ERK axis and sequestering IGF2BP1. Oncogene. 2021;40:277–91.
Zhou X, Gao W, Hua H, Ji Z. LncRNA-BLACAT1 Facilitates Proliferation, Migration and Aerobic Glycolysis of Pancreatic Cancer Cells by Repressing CDKN1C via EZH2-Induced H3K27me3. Front Oncol. 2020;10:539805.
Gu Z, Hou Z, Zheng L, Wang X, Wu L, Zhang C. LncRNA DICER1-AS1 promotes the proliferation, invasion and autophagy of osteosarcoma cells via miR-30b/ATG5. Biomed Pharmacother. 2018;104:110–8.
Ma C, Ma N, Qin L, Miao C, Luo M, Liu S. DICER1-AS1 Promotes the Malignant Behaviors of Colorectal Cancer Cells by Regulating miR-296-5p/STAT3 Axis. Cancer Manag Res. 2020;12:10035–46.
Faghihi MA, Wahlestedt C. Regulatory roles of natural antisense transcripts. Nat Rev Mol Cell Biol. 2009;10:637–43.
Katayama S, Tomaru Y, Kasukawa T, Waki K, Nakanishi M, Nakamura M, Nishida H, Yap CC, Suzuki M, Kawai J, et al. Antisense transcription in the mammalian transcriptome. Science. 2005;309:1564–6.
Wang X, Zhao J, Huang J, Tang H, Yu S, Chen Y. The regulatory roles of miRNA and methylation on oncogene and tumor suppressor gene expression in pancreatic cancer cells. Biochem Biophys Res Commun. 2012;425:51–7.
Lei B, Wang D, Zhang M, Deng Y, Jiang H, Li Y. miR-615-3p promotes the epithelial-mesenchymal transition and metastasis of breast cancer by targeting PICK1/TGFBRI axis. J Exp Clin Cancer Res. 2020;39:71.
Shi Y, Seto E, Chang LS, Shenk T. Transcriptional repression by YY1, a human GLI-Krüppel-related protein, and relief of repression by adenovirus E1A protein. Cell. 1991;67:377–88.
Zhang JJ, Zhu Y, Zhang XF, Liu DF, Wang Y, Yang C, Shi GD, Peng YP, Zhang K, Tian L, et al. Yin Yang-1 suppresses pancreatic ductal adenocarcinoma cell proliferation and tumor growth by regulating SOX2OT-SOX2 axis. Cancer Lett. 2017;408:144–54.
Zhang JJ, Zhu Y, Xie KL, Peng YP, Tao JQ, Tang J, Li Z, Xu ZK, Dai CC, Qian ZY, et al. Yin Yang-1 suppresses invasion and metastasis of pancreatic ductal adenocarcinoma by downregulating MMP10 in a MUC4/ErbB2/p38/MEF2C-dependent mechanism. Mol Cancer. 2014;13:130.
Chen Q, Zhang JJ, Ge WL, Chen L, Yuan H, Meng LD, Huang XM, Shen P, Miao Y, Jiang KR. YY1 inhibits the migration and invasion of pancreatic ductal adenocarcinoma by downregulating the FER/STAT3/MMP2 signaling pathway. Cancer Lett. 2019;463:37–49.
Li Q, Pan X, Zhu D, Deng Z, Jiang R, Wang X. Circular RNA MAT2B Promotes Glycolysis and Malignancy of Hepatocellular Carcinoma Through the miR-338-3p/PKM2 Axis Under Hypoxic Stress. Hepatology. 2019;70:1298–316.
Fong MY, Zhou W, Liu L, Alontaga AY, Chandra M, Ashby J, Chow A, O’Connor ST, Li S, Chin AR, et al. Breast-cancer-secreted miR-122 reprograms glucose metabolism in premetastatic niche to promote metastasis. Nat Cell Biol. 2015;17:183–94.
Bai X, Lu D, Lin Y, Lv Y, He L. A seven-miRNA expression-based prognostic signature and its corresponding potential competing endogenous RNA network in early pancreatic cancer. Exp Ther Med. 2019;18:1601–8.
Zuo X, Chen Z, Gao W, Zhang Y, Wang J, Wang J, Cao M, Cai J, Wu J, Wang X. M6A-mediated upregulation of LINC00958 increases lipogenesis and acts as a nanotherapeutic target in hepatocellular carcinoma. J Hematol Oncol. 2020;13:5.
We sincerely appreciate the researchers who worked for this study.
Full list of author information is available at the end of the article
This research was supported by the Clinical Research Physician Program of Tongji Medical College, Huazhong University of Science and Technology and the National Science Foundation of China [81672406, 81872030].
The Ethics Committee at the Academic Medical Center of Huazhong University of Science and Technology approved this study. The Animal Research Committee of the Academic Medical Center at Huazhong University of Science and Technology reviewed all animal studies, which were consistent with Institutional and Animal Care and Use Committee guidelines.
Consent for publication
All authors agree with the manuscript content.
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The sequences of PCR primers. Table S2. The sequences for gene knockdown. Table S3. Clinical information of 86 cases PAAD patients (cohort 1). Table S4. Univariate and multivariate analyses of DICER1-AS1 level and overall survival (TCGA). Table S5. Correlation between miR-5586-5p expression and clinicopathological features in PAAD patients (TCGA). Table S6. Univariate and multivariate analyses of the miR-5586-5p level and overall survival (TCGA).
Figure S1. DICER1-AS1 is downregulated in pancreaticcancer tissues and negatively correlated with glycolysis pathway. Figure S2. Characteristics of lncRNA DICER1-AS1. Figure S3.Knockdown of DICER1-AS1 promotes glycolysis, proliferationand metastasis of PC cells. Figure S4. The expression correlation between DICER1-AS1and DICER1 and Kaplan-Meier survival curves of DICER1 in tumor datasets. Figure S5. DICER1 is a critical targetfor DICER1-AS1 regulating glycolysis. Figure S6. The expressioncorrelation between YY1 and DICER1 in the CCLE databaseand Kaplan-Meier survival curves of YY1 in tumordatasets. Figure S7. DICER1-AS1 regulates DICER1 transcription byrecruiting YY1 to the promoter of DICER1. Figure S8.miR-5586-5pis an essential target for the DICER1-AS1/DICER1 pathway. Figure S9. miR-5586-5p regulates glycolysis by inducingmRNA degradation of glycolytic genes.
About this article
Cite this article
Hu, Y., Tang, J., Xu, F. et al. A reciprocal feedback between N6-methyladenosine reader YTHDF3 and lncRNA DICER1-AS1 promotes glycolysis of pancreatic cancer through inhibiting maturation of miR-5586-5p. J Exp Clin Cancer Res 41, 69 (2022). https://doi.org/10.1186/s13046-022-02285-6
- Pancreatic cancer