- Open Access
The miR-1185-2-3p—GOLPH3L pathway promotes glucose metabolism in breast cancer by stabilizing p53-induced SERPINE1
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 47 (2021)
Phosphatidylinositol-4-phosphate-binding protein GOLPH3L is overexpressed in human ductal carcinoma of the breast, and its expression levels correlate with the prognosis of breast cancer patients. However, the roles of GOLPH3L in breast tumorigenesis remain unclear.
We assessed the expression and biological function of GOLPH3L in breast cancer by combining bioinformatic prediction, metabolomics analysis and RNA-seq to determine the GOLPH3L-related pathways involved in tumorigenesis. Dual-luciferase reporter assay and coimmunoprecipitation (Co-IP) were used to explore the expression regulation mechanism of GOLPH3L.
We demonstrated that knockdown of GOLPH3L in human breast cancer cells significantly suppressed their proliferation, survival, and migration and suppressed tumor growth in vivo, while overexpression of GOLPH3L promoted aggressive tumorigenic activities. We found that miRNA-1185-2-3p, the expression of which is decreased in human breast cancers and is inversely correlated with the prognosis of breast cancer patients, is directly involved in suppressing the expression of GOLPH3L. Metabolomics microarray analysis and transcriptome sequencing analysis revealed that GOLPH3L promotes central carbon metabolism in breast cancer by stabilizing the p53 suppressor SERPINE1.
In summary, we discovered a miRNA-GOLPH3L-SERPINE1 pathway that plays important roles in the metabolism of breast cancer and provides new therapeutic targets for human breast cancer.
Breast cancer is the most commonly diagnosed and most deadly cancer affecting women worldwide . Tumor cells present a number of characteristics, such as self-proliferation ability , apoptosis resistance , unlimited replication potential , insensitivity to growth inhibition , continuous angiogenesis , tissue invasion  and metastasis . Mammary tumorigenesis is a multistep process involving activation of oncogenes or inactivation of tumor suppressors, abnormal expression of noncoding RNA , loss of genome stability  and various genetic and epigenetic alterations . Tumor cells prefer glycolysis for glucose metabolism under aerobic conditions (Warburg effect) rather than mitochondrial oxidative phosphorylation, which is more efficient for ATP production so that tumor cells can produce more energy and various metabolites in a short period of time . ATP produced by glycolysis can satisfy high energy demand for rapid tumor growth. The intermediate products of glycolysis, such as glucose 6-phosphate and pyruvic acid, which can synthesize fatty acids and nucleic acids, regulate cell metabolism and biosynthesis. Moreover, glycolysis products acidify the microenvironment around the tumor, which is not conducive to the growth of normal cells but promotes tumor invasion and metastasis. The tumor suppressor protein p53 is the “guardian of the genome”, which plays critical roles in cell cycle regulation, DNA repair, cell differentiation and apoptosis . Tumor cells with p53 inactivation often show increasing glycolytic activity. Accumulating data show that p53 may confer tumor suppression by inhibiting the cancer metabolic switch from oxidative phosphorylation to glycolysis.
The Golgi apparatus is an important part of the cell membrane system involved in protein glycosylation, proteolytic activation and cellular secretory activity. Protein glycosylation is one of the most common posttranslational modifications, which regulates the location, function, activity, life span and diversity of proteins in tissues and cells and participates in various important life activities, such as cell recognition, differentiation, development, signal transduction and immune response. Glycosylation labels different proteins and changes the conformation of polypeptides to increase the stability of proteins. Mammalian proteins have three types of glycosylation: N-glycosylation (N-GlcNAc), O-glycosylation (O-GlcNAc) and glycosylphosphatidylinositol (GPI) anchor. In N-glycosylation, the sugar chain is covalently linked to the free NH2 group of aspartic acid of the protein. The synthesis of the N-linked sugar chain starts from the endoplasmic reticulum (ER) and is completed at the Golgi. The glycoproteins formed by the ER have similar sugar chains. After entering the Golgi apparatus from the cis surface, most of the mannose on the original sugar chain is removed during the transport process between the membrane capsules. However, different types of sugar molecules are added to different glycosyltransferases to form oligosaccharide chains with different structures. N-glycan biosynthesis coordinates the cellular response of tumor cells, determining growth, invasion and drug sensitivity ; for instance, N-acetyllactosamine N-glycans mediate PD-L1 and PD-1, affecting the efficacy of anti-PD-L1 therapies . Altering the N-glycosylation of integrins affects cis-interactions with important membrane receptors, such as EGFR, contributing to tumor cell motility and migration . The N-glycosylation product Fut8 is involved in the expression of cancer biomarkers as well as in the treatment of cancer, and GnT-V is highly associated with cancer metastasis, whereas GnT-III is associated with cancer suppression . O-linked glycosylation takes place in the Golgi apparatus. The sugar chain is covalently linked to the free OH radical of serine/threonine residues in proteins. O-GlcNAc is responsible for cancer progression. O-GlcNAc modifications regulate the activities of FoxM1 and cyclin D1, which are involved in cell cycle progression and critical to cell proliferation . O-GlcNAc has been implicated in cancer cell survival through the effect of hyper-O-GlcNAc via activation of κB-mediated signaling . Moreover, O-GlcNAc participates in cancer cell invasion and metastasis by regulating E-cadherin trafficking and function . Increased levels of O-GlcNAc transferase (OGT) have been found in breast cancer . O-GlcNAcylation serves as a nutrient sensor to modulate crosstalk with phosphorylation, such as p53 . The GPI glycosylphosphatidylinositol anchor is the only way for proteins to combine with the cell membrane, which is different from the general lipid modifying components, and its structure is extremely complex. Mutation in the X-linked PIGA gene increases the risk of developing leukemia . The spatial structure of glycoprotein determines unique glycosyltransferases to initiate certain glycosylation modifications.
Current studies show that Golgi apparatus dysfunction is associated with tumor development, but more than 80% of Golgi-related proteins have not been reported to play a role in this process. Golgi phosphoprotein 3-like (GOLPH3L) is an important protein involved in the formation of Golgi vesicles and their anterograde transport to the plasma membrane. GOLPH3L is highly expressed in various tumor tissues and promotes the proliferation of rhabdomyosarcoma cells  and is involved in the regulation of proliferation, apoptosis and the cell cycle in cervical cancer cells . Mechanistically, GOLPH3L has been reported as an activator of the NF-κB signaling pathway in ovarian cancer . Considering the critical roles of glycolysis in tumorigenesis, we demonstrate here that GOLPH3L contributes to tumorigenesis by promoting glucose metabolism in breast cancer by stabilizing certain downstream proteins of p53.
BRC patient samples
All patient-related studies were approved by the Institutional Review Board of Taishan People’s Hospital and Nanfang Hospital of Southern Medical University.
Human Cancer cell Xenograft model
All animal work was approved by the Institutional Animal Care and Use Committee (IACUC) of Southern Medical University. A total of 5 × 106 breast cancer cells were implanted into the skeletal muscle of the hind limbs of 3 ~ 4-week-old BALB/c nude mice (nu/nu). One week after transplantation, the diameter of tumors was measured every 3 days. Tumors were recovered and weighed after 3 weeks.
Human normal mammary epithelial cell line (MCF-10A) and human breast cancer cell lines (T47D, BT474, MCF-7, MDA-MB-231, SK-BR-3) were purchased from American Type Culture Collection (ATCC, Manassas, VA, USA). Human breast cancer cell lines (T47D, BT474, MCF-7, MDA-MB-231, SK-BR-3) were cultured in Roswell Park Memorial Institute 1640 medium (Gibco, USA) supplemented with 10% fetal bovine serum (FBS, HyClone, Utah, USA) and 1% penicillin-streptomycin (Pen/Strep) (Gibco) at 37 °C with 5% CO2. The human normal mammary epithelial cell line MCF-10A was cultured in media and supplements from the MEGM kit (Lonza/Clonetics, CC-3150) and 10% FBS - supplemented with 100 ng/ml cholera toxin (Sigma-Aldrich, C8052) at 37 °C with 5% CO2.
Establishment of transfected cell lines
The vectors expressing GOLPH3L-specific siRNA (RIBOBIO, Cat# siG000055204A-C) and SERPINE1-specific siRNA (ThermoFisher, Cat# 4390771). Vectors expressing human GOLPH3L cDNA and SERPINE1 cDNA were transfected into cells as previously described . Cells were selected with puromycin (2 μg/ml, GeneChem) for 3 days for stable transfection.
Western blot analysis
Cells were extracted for total protein analysis using lysis buffer, and samples were separated on 8–15% SDS-PAGE and transferred to nitrocellulose membranes, which were blocked with blocking buffer (5% skim milk in PBS with 0.05% Tween 20) and incubated with primary antibody in the blocking buffer. After being washed three times with blocking buffer, the membrane was probed with secondary antibody and developed with Supersignal West Pico or Dura (Thermo Fisher Scientific).
Quantitative PCR analysis
Real-time PCR analysis was performed using the StepOnePlus Real-Time PCR System (Applied Biosystems) with FastStart Universal SYBR Green Master Mix (Roche) as previously reported. The primer sets used were as follows:
miRNA-1185-2-3p mimic-F: AUAUACAGGGGGAGACUCUCAU;
miRNA-1185-2-3p mimic-R: GAGAGUCUCCCCCUGUAUAUUU;
miRNA-1185-2-3p inhibitor: AUGAGAGUCUCCCCCUGUAUAU.
The PCR conditions were as follows: 10 min at 95 °C, 40 cycles of 15 s at 95 °C, and 1 min at 60 °C. The average Ct value for each gene was determined from triplicate reactions and normalized to that of β-actin for genes and U6 for microRNAs (miRNAs).
Cell proliferation, apoptosis, migration assay and cell cycle assay
For EdU-high content screening of the cellular proliferation assay, cells were prepared by trypsinization and seeded onto 96-well plates at a density of 1 × 104 cells per well. After incubation for 48 h at 37 °C, the old medium was discarded, and 100 μl of medium containing EdU was added to each well and incubated for 6 h. Cells were fixed with paraformaldehyde at room temperature for 20 min. Then, 100 μl of 2 mg/ml glycine solution and 100 μl of 0.5% Triton X-100 solution were added to the wells separately, and the cells were washed twice with PBS. One hundred microliters of Apollo-643 staining solution was added in each well for 30 min and then discarded. After destaining and rinsing, DAPI reaction mixture was added to each well and incubated for 30 min. A GE IN CELL Analyzer 6500HS Confocal High Content Imaging Analyzer was used for collecting images.
Consistent treatments were applied to prepared cells in 96-well plates: 100 μl 2 mg/ml glycine solution and 100 μl 0.5% Triton X-100 solution were added into wells separately, and then the cells were washed twice with PBS. Fluorescein-dUTP solution (50 μl) mixed with TdT enzyme was added to each well and incubated for 2 h at room temperature. Then, the reaction solution was discarded, and 100 μl 2 × SSC buffer was immediately added into each well. A 1× DAPI reaction mixture was added to the plates and incubated for 30 min. A GE IN CELL Analyzer 6500HS Confocal High Content Imaging Analyzer was used for collecting images after the cells were washed three times with PBS.
Consistent treatments were applied to prepared cells in 24-well plates at a density of 2.5 × 104 cells per well. After adjusting the cell concentration to 1 × 105 cells/ml, the cells were seeded into the upper chamber at 1 × 104 cells/μl, and the lower chamber was immediately filled with 150 μl complete medium with 10% FBS as a chemoattractant and then incubated for 48 h. After the medium was discarded from the lower chamber, 150 μl PBS was added to each well for 5 min. Calcein AM cell staining solution at a working concentration of 2.5 μM was prepared with complete medium, and 150 μl staining solution was added to each well for incubation at room temperature for 30 min after the PBS was discarded. After incubation, 150 μl trypsin was added to the lower chamber and incubated at 37 °C for 15 min. Serum-containing medium was added to the chamber to stop digestion. Cells in the subchamber were aspirated for cell counting.
For the cell cycle assay, cells were prepared in 96-well plates at a density of 1 × 104 cells per well. After being washed with PBS one time, the cells were fixed with 4% paraformaldehyde for 20 min; the paraformaldehyde was removed, and 100 μl of 2 mg/ml glycine solution was added. The glycine solution was subsequently removed, and the cells were washed with PBS one time. Then, 100 μl of 0.5% Triton X-100 PBS solution was added to each well, which were then washed twice with PBS. One hundred microliters of 1× DAPI reaction mixture was added to each well, and the plates were incubated for 30 min. After being washed with PBS three times, a GE IN CELL Analyzer 6500HS Confocal High Content Imaging Analyzer was used to collect the images.
Dual-luciferase reporter assay
Cells were seeded in triplicate onto 6-well plates at a density of 4 × 105 cells/well for 48 h and transfected with 0.3 μg of REPOTM-AP-1-luc plasmid and control-luciferase plasmid separately together with 30 ng of pGMR TK renilla plasmid (GenomeDitech, Shanghai, China) using Lipofectamine™ 3000 reagent (Invitrogen, Carlsbad, USA). Firefly and Renilla luciferase activities were measured using the Dual-Luciferase Reporter Assay Kit (Promega, Madison, USA) after 48 h of transfection.
Cells were collected in centrifuge tubes filled with 1 ml methanol:acetonitrile:water (2:2:1, V/V) at a density of 2 × 107 cells/tube after washing with PBS and physiological saline solution. The liquid nitrogen is stored at − 80 °C after quick freezing. To analyze metabolomics, samples were analyzed at Applied Protein Technology (APT, Shanghai). Partial least squares discrimination analysis was applied to reveal the relationship between the expression of metabolites and sample types by calculating the cast variable importance for the projection (VIP) to measure the expression pattern of each metabolite for each group. All differentially expressed metabolites were selected for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses (VIP score > 1.0). GO was performed with KOBAS3.0 software. The differentially expressed metabolites and enriched pathways were mapped using the KEGG pathways with KOBAS3.0 software (http://www.genome.jp/kegg).
Total RNA was isolated from cells/tissues using TRIzol (Invitrogen) according to the manufacturer’s protocol. RNA purity was assessed using the ND-1000 Nanodrop. Each RNA sample had an A260:A280 ratio above 1.8 and an A260:A230 ratio above 2.0. RNA integrity was evaluated using the Agilent 2200 TapeStation (Agilent Technologies, USA), and each sample had an RNA integrity number (RIN) above 7.0. Briefly, rRNAs were removed from total RNA using the EpicentreRibo-Zero rRNA Removal Kit (Illumina, USA) and fragmented to approximately 200 bp. Subsequently, the purified RNAs were subjected to first strand and second strand cDNA synthesis followed by adaptor ligation and enrichment with a low-cycle according to the instructions of the NEBNext® Ultra™ RNA Library Prep Kit for Illumina (NEB, USA). The purified l library products were evaluated using the Agilent 2200 TapeStation and Qubit®2.0 (Life Technologies, USA) and then diluted to 10 pM for cluster generation in situ on the pair-end flow cell followed by sequencing (2 × 150 bp) with a HiSeq3000. Clean reads were obtained after removal of reads containing adaptor, poly-N and low quality from raw data. HISAT2 was used to align the clean reads to the mouse reference genome mm10 with default parameters. HTSeq was subsequently employed to convert aligned short reads into read counts for each gene model. Differential expression was assessed by DEseq using read counts as input. The Benjamini-Hochberg multiple test correction method was enabled. Differentially expressed genes were chosen according to the criteria of fold change > 2 and adjusted p-value < 0.05. All the differentially expressed genes were used for heat map analysis and KEGG ontology enrichment analyses. For KEGG enrichment analysis, a p-value < 0.05 was used as the threshold to determine significant enrichment of the gene sets.
To measure the levels of glycolytic ATP production, 1 × 104 cells were seeded into each well of black 96-well plates. To normalize the levels of protein, the same number of cells were seeded into clear bottom 96-well plates. Cells were incubated in medium containing 1 μM oligomycin (Sigma-Aldrich) to inhibit mitochondrial oxidative ATP production or 25 mM 2-deoxy-D-glucose (2-DG) to inhibit glycolytic ATP production. After washing the cells with PBS, ATP levels were measured using a kit according to the manufacturer’s protocol (PerkinElmer). Total ATP production was calculated by subtracting the amount of ATP in cells treated with both oligomycin and 2-DG from the amount of ATP in cells without treatment. To normalize the number of cells, the protein concentration was measured using the Bicinchoninic Acid Protein Assay Kit (BCA, Sigma-Aldrich) according to the manufacturer’s protocols.
Protein stability analysis
Cells were transfected with siRNAs, and cyclohexamide (CHX; 1:1000) was used to treat cells, which were harvested at various time points (0 h, 0.5 h, 1 h, 2 h, 4 h, 6 h and 12 h). Levels of various proteins were determined by Western blot analysis and quantified with ImageJ software.
Immunoprecipitation assays were performed as previously described . Briefly, cells were lysed in RIPA duffer containing protease and phosphatase inhibitors, and cells were collected after centrifugation at 12,000×g for 10 min at 4 °C. Supernatants were immunoprecipitated with antibodies followed by incubation with magnetic protein A/G beads (Pierce) for 2 h at 4 °C. The immune complexes were washed three times with PBS, resuspended in SDS-PAGE buffer and assessed by Western blot analysis.
Data were analyzed using SPSS 20.0. and two-tailed independent Student’s t-test; P < 0.05 was considered significant. Two patient cohorts were compared by Kaplan-Meier survival plot, and log-rank p-values were calculated.
GOLPH3L is highly expressed in breast Cancer and promotes the tumorigenesis of breast Cancer cells
Using Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn/index.html) and profiling the expression data of GOLPH3L in breast tumors and tumor-adjacent normal tissues in the database (GSE93601), we confirmed that GOLPH3L was dramatically overexpressed in breast cancer samples compared to normal control samples (Fig. 1a). Moreover, the expression levels of GOLPH3L were inversely correlated with the prognosis of human breast cancer patients (Fig. 1b). In this context, GOLPH3L promotes breast tumorigenesis. To explore the explicit roles of GOLPH3L in breast tumorigenesis, we investigated the mRNA and protein expression levels of GOLPH3L in breast cancer tissues and paired adjacent normal tissues as well as in various breast cancer cell lines and normal breast epithelial cell lines. The relative mRNA and protein expression of GOLPH3L in cancer tissues was much higher than that in adjacent normal tissues (Fig. 1c). Consistently, compared to the normal breast epithelial cell line MCF-10A, the expression levels of GOLPH3L mRNA and protein were increased in breast cancer cell lines, especially in T47D and BT474 cells (Fig. 1d).
Therefore, we examined the roles of GOLPH3L in T47D and BT474 cells by altering the expression of GOLPH3L through knockdown or overexpression (Figure S1a and b). Restraining GOLPH3L expression in breast cells decreased cellular proliferation and survival, while upregulating the expression of GOLPH3L in T47D cells promoted proliferation and survival (Fig. 2a and b). Transwell assays revealed that silencing GOLPH3L inhibited the migration of T47D cells, and in contrast, overexpressing GOLPH3L in cells obviously facilitated cancer cell migration (Fig. 2c). Moreover, knockdown of GOLPH3L expression significantly blocked the cell cycle at the G0/G1 phase (Fig. 2d). Further evidence of the role of GOLPH3L in tumorigenesis is that knockdown of GOLPH3L in T47D cells significantly reduced the growth of tumors formed by T47D cells in immunodeficient mice (Fig. 2e). Consistent data were obtained using BT474 breast cancer cells (Figure S1c-f). These findings demonstrate that GOLPH3L promotes the tumorigenesis of breast cancer cells in various ways.
miRNA-1185-2-3p negatively regulates GOLPH3L and positively correlates with the prognosis of human breast Cancer patients
miRNAs are 19–25 nucleotide-long noncoding RNAs that regulate gene expression via messenger RNA degradation and translation. miRNAs are involved in every aspect of biological processes, and it is worth noting that alterations in miRNA expression induce the initiation, progression and metastasis of human tumors. Considering the consequence of GOLPH3L in promoting breast tumorigenesis, we utilized the miRDB prediction program (http://mirdb.org/). Based on the predicted results (Figure S2a), we mutated the target sites of the top five predicted miRNAs in the 3′-UTR of GOLPH3L mRNA. A dual-luciferase reporter assay confirmed that miRNA-1185-2-3p regulated the levels of GOLPH3L mRNA by targeting the predicted sites within the 3′-UTR of GOLPH3L (Fig. 3a). Based on the GEO dataset GSE45666, we found that miRNA-1185-2-3p was expressed at a higher level in tumor-adjacent normal tissues than in breast tumor tissues and was inversely correlated with the prognosis of breast cancer patients (Fig. 3b, left panel). Consistently, the expression of miRNA-1185-2-3p was higher in various breast cancer cell lines than in the normal breast epithelial cell line MCF-10A (Fig. 3b, right panel). To determine the impact of miRNA-1185-2-3p on the expression of GOLPH3L, we used miRNA mimics and inhibitors to demonstrate that the expression levels of miRNA-1185-2-3p were inversely correlated with the expression of GOLPH3L in breast cancer lines (Fig. 3c and S2b). These results support the notion that miRNA-1185-2-3p directly targets GOLPH3L mRNA.
Since miRNA-1185-2-3p suppresses the expression of GOLPH3L in breast cancer cells, we induced it in the breast cancer cell lines T47D and BT474 to determine whether miRNA-1185-2-3p could have the same tumor suppressive effects as GOLPH3L knockdown. The overexpression of miRNA-1185-2-3p also inhibited proliferation, survival, and migration and arrested the cell cycle of breast cancer cells (Fig. 3d-g and S2c-f). In addition, the induction of miRNA-1185-2-3p suppressed the tumorigenesis induced by GOLPH3L overexpression (Figure S3). Therefore, these data indicate that miRNA-1185-2-3p could inhibit breast tumorigenesis by suppressing the expression of GOLPH3L.
GOLPH3L affects glucose metabolism in breast Cancer cells
To understand the mechanism by which GOLPH3L promotes tumorigenesis, we combined metabolomic analysis and transcriptome sequencing in GOLPH3L-silenced T47D breast cancer cells. We identified metabolites that were altered significantly after GOLPH3L silencing. Notably, glycolytic intermediates such as acetyl coenzyme A (acetyl-CoA), adenosine monophosphate (AMP), adenosine 5′-diphosphate (ADP) and adenosine 5′-triphosphate (ATP) were reduced when GOLPH3L expression was downregulated (Fig. 4a, S4a and Table 1). KEGG pathway analysis showed that GOLPH3L expression could markedly influence central carbon metabolism in cancer, glycolysis/gluconeogenesis and 61 other pathways (Fig. 4b and Table 2). Correlation analysis of different metabolites revealed that acetyl-CoA, AMP, ADP and ATP were positively correlated in both positive ion mode (Fig. 4c, top panel) and negative ion mode (Fig. 4c, bottom panel). Since GOLPH3L knockdown decreased the proliferation of breast cancer cells, we evaluated the impact of GOLPH3L on metabolism and found that the suppression of GOLPH3L expression decreased glycolytic activity in breast cell lines (Fig. 4d and e).
Correspondingly, high-throughput transcriptome sequencing revealed that GOLPH3L expression could markedly alter the expression of 205 genes in breast cancer cells (Fig. 5a and Table 3), and the central carbon metabolism pathway was one of the downregulated pathways in GOLPH3L-knockdown cancer cells (Fig. 5b and Table 4). Because GOLPH3L is located in the Golgi apparatus and may participate in mediating recruitment to Golgi membranes, we predicted that GOLPH3L likely functions by regulating the stability of other proteins involved in tumorigenesis via certain types of glycosylation. Therefore, we chose SERPINE1 as a candidate target through comprehensively analyzing the upregulated tumor suppressor genes and downregulated oncogenes that were related to the central carbon metabolism pathway and had glycosylation as a posttranslational modification (Fig. 5c and d).
GOLPH3L stabilizes p53-induced SERPINE1 expression in breast Cancer cells is positively correlated with increased glycolysis
Published data  and our KEGG analysis results (Figure S4c) indicated that the p53 signaling pathway regulates SERPINE1 and therefore suppresses tumor cell proliferation, invasion and migration, whereas it promotes cell apoptosis. To assess the role of SERPINE1 in breast cancer, we examined the expression level of SERPINE1. When compared to those in the tumor-adjacent normal tissues and normal breast epithelial cell line MCF-10A, the expression levels of SERPINE1 were increased in breast cancer tissues and breast cancer cell lines (Fig. 5e). In addition, the protein levels of SERPINE1 were correlated with the protein levels of GOLPH3L in breast cancer cells, supporting the notion that GOLPH3L may regulate the expression of SERPINE1 (Fig. 5f). Moreover, the expression of SERPINE1 promoted breast tumorigenesis (Fig. 5g). Using a coimmunoprecipitation (Co-IP) assay, we confirmed the interaction between GOLPH3L and SERPINE1 in breast cancer cells (Fig. 6a). To test whether GOLPH3L can stabilize SERPINE1, the alteration of the expression of GOLPH3L significantly affected the half-life of SERPINE1 in a breast cancer cell line, indicating that GOLPH3L contributes to the stabilization of SERPINE1 (Fig. 6b). Subsequently, we altered the expression of GOLPH3L in T47D cells, and the IP results showed increased ubiquitination of SERPINE1 in the presence of GOLPH3L (Fig. 6c). Furthermore, the overexpression of SERPINE1 reversed the antitumor activities induced by the suppression of GOLPH3L in breast cancer cell lines (Fig. 6d and e). Therefore, GOLPH3L promotes glucose metabolism in breast cancer and is conducive to SERPINE1 stabilization.
Breast cancer cells are produced by mutations in the DNA of normal breast cells, genetics and other factors, such as diet and exercise. Lifestyle habits can also be the cause of DNA mutations. Some DNA mutations are heritable, but the vast majority of DNA changes are acquired, and such acquired mutations occur only in breast cancer cells . Mutated DNA causes corresponding changes in genes that regulate cell proliferation and apoptosis, which can lead to uncontrolled proliferation and apoptosis, thus promoting the development of tumors. According to the New York Cancer Institute , targeted therapeutics can be categorized according to their targets: those that target the regulatory mechanism of tumor formation, those that target the tumor microenvironment, those that target the tumor immune system, those that target tumor markers and those that target tumor stem cells. Dysregulation of targets that regulate tumor formation can stimulate a series of downstream signaling pathways, resulting in abnormal division, growth, metabolism, cell environment and angiogenesis of tumor cells compared with normal cells.
Golgi phosphoprotein 3-like (GOLPH3L) is a prognostic biomarker of cervical cancer  and ovarian cancer  and may be a mitochondrial biogenesis marker in breast cancer metabolism . However, the roles of GOLPH3L in breast tumorigenesis remain unclear. SERPINE1 promotes breast cancer cell metastasis  and glycolytic metabolism in triple-negative breast cancer (TNBC)  and participates in EGFR signaling . We discovered that GOLPH3L promotes glucose metabolism in breast cancer cells by stabilizing SERPINE1, which is regulated by p53 transcription . In this context, we demonstrate that silencing the GOLPH3L gene suppresses breast tumorigenesis and glucose metabolism, while overexpression of GOLPH3L promotes breast tumorigenesis. Therefore, GOLPH3L-mediated stabilization of SERPINE1 could represent an important oncogenic pathway in the glucose metabolism of breast cancer. Our findings suggest that GOLPH3L functions as an oncogene in breast cancer by promoting cellular proliferation and migration and suppressing apoptosis. Considering the tumorigenic role of GOLPH3L, we also investigated the upstream regulators of GOLPH3L. Fortunately, we identified that miR-1185-2-3p directly targets the sites within the 3′-UTR of GOLPH3L and functionally suppresses the expression of GOLPH3L. In support of this notion, we demonstrated that miR-1185-2-3p inhibits breast tumorigenesis by suppressing GOLPH3L expression. This is the first time that miR-1185-2-3p has been identified as a regulator in cancer development. Therefore, it will be important to identify the pathway related to miR-1185-2-3p.
The treatment and prognosis of breast cancer differs according to the expression of different molecular makers such as estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor 2 (HER2). HER2+ breast cancers tend to grow and spread more aggressively. Different types of drugs such as monoclonal antibodies and kinase inhibitors have been developed to target the HER2 protein [36, 37]. Approximately 2 of 3 breast cancers are hormone receptor-positive (ER+ or PR+). Treatment with hormone therapy is often helpful in these cases, and certain targeted therapy drugs, such as CDK4/6 inhibitors , mTOR inhibitor (Everolimus) and PI3K inhibitor (Alpelisib) , can increase the efficacy of hormone therapy. Among all types of breast cancer, triple-negative breast cancer (TNBC) has the worst prognosis due to the lack of an effective target. Moreover, over 30% of breast cancer patients may suffer recurrence. In terms of metastatic breast cancer, systemic therapy usually shows unsatisfactory treatment effects. Therefore, it is of great clinical significance to discover new therapeutic targets for breast cancer. Oncogenic proteins such as hypoxia-inducible factor, Myc, p53 and PI3K/Akt/mTOR pathway proteins, which are involved in metabolic reprogramming , may serve as candidate anticancer agents. Several preclinical trials have shown effectiveness of an inhibitor targeting the glycolytic pathway . Our discovery provides a mechanistic link between miR-1185-2-3p, GOLPH3L and SERPINE1; this pathway plays a significant role in glucose metabolism in breast cancer and may serve as a novel therapeutic target for breast cancer.
We discovered a functional pathway linking miR-1185-2-3p, GOLPH3L and SERPINE1 that plays a significant role in glucose metabolism in breast cancer and provides new therapeutic targets for breast cancer treatment.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Golgi phosphoprotein 3-like
Serpin family E member 1
Epidermal growth factor receptor
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.
Ren R, Sun H, Ma C, Liu J, Wang H. Colon cancer cells secrete exosomes to promote self-proliferation by shortening mitosis duration and activation of STAT3 in a hypoxic environment. Cell Biosci. 2019;9:62.
Janzen DM, Tiourin E, Salehi JA, Paik DY, Lu J, Pellegrini M, et al. Retraction note: an apoptosis-enhancing drug overcomes platinum resistance in a tumour-initiating subpopulation of ovarian cancer. Nat Commun. 2020;11:2218.
Nyquist MD, Corella A, Coleman I, De Sarkar N, Kaipainen A, Ha G, et al. Combined TP53 and RB1 loss promotes prostate cancer resistance to a spectrum of therapeutics and confers vulnerability to replication stress. Cell Rep. 2020;31:107669.
Overdijk MB, Verploegen S, van den Brakel JH, van Bueren JJL, Vink T, van de Winkel JG, et al. Epidermal growth factor receptor (EGFR) antibody-induced antibody-dependent cellular cytotoxicity plays a prominent role in inhibiting tumorigenesis, even of tumor cells insensitive to EGFR signaling inhibition. J Immunol. 2011;187:3383–90.
Niu Y, Bao L, Chen Y, Wang C, Luo M, Zhang B, et al. HIF2-induced long noncoding RNA RAB11B-AS1 promotes hypoxia-mediated angiogenesis and breast cancer metastasis. Cancer Res. 2020;80:964–75.
Shaker H, Bundred NJ, Landberg G, Pritchard SA, Albadry H, Nicholson SL, et al. Breast cancer stromal clotting activation (tissue factor and thrombin): a pre-invasive phenomena that is prognostic in invasion. Cancer Med. 2020;9:1768–78.
Zimmer AS, Van Swearingen AED, Anders CK. HER2-positive breast cancer brain metastasis: a new and exciting landscape. Cancer Rep (Hoboken). 2020:e1274. https://doi.org/10.1002/cnr2.1274.
Cantile M, Di Bonito M, Cerrone M, Collina F, De Laurentiis M, Botti G. Long non-coding RNA HOTAIR in breast cancer therapy. Cancers (Basel). 2020;12:1197.
Di LJ, Byun JS, Wong MM, Wakano C, Taylor T, Bilke S, et al. Genome-wide profiles of CtBP link metabolism with genome stability and epithelial reprogramming in breast cancer. Nat Commun. 2013;4:1449.
Salas LA, Johnson KC, Koestler DC, O'Sullivan DE, Christensen BC. Integrative epigenetic and genetic pan-cancer somatic alteration portraits. Epigenetics. 2017;12:561–74.
Tahtouh R, Wardi L, Sarkis R, Hachem R, Raad I, El Zein N, et al. Glucose restriction reverses the Warburg effect and modulates PKM2 and mTOR expression in breast cancer cell lines. Cell Mol Biol (Noisy-le-grand). 2019;65:26–33.
Cao J, Liu X, Yang Y, Wei B, Li Q, Mao G, et al. Decylubiquinone suppresses breast cancer growth and metastasis by inhibiting angiogenesis via the ROS/p53/ BAI1 signaling pathway. Angiogenesis. 2020;23:325–38.
Dennis JW, Nabi IR, Demetriou M. Metabolism, cell surface organization, and disease. Cell. 2009;139:1229–41.
Li CW, Lim SO, Chung EM, Kim YS, Park AH, Yao J, et al. Eradication of triple-negative breast cancer cells by targeting glycosylated PD-L1. Cancer Cell. 2018;33:187–201.e10.
Vlodavsky I, Friedmann Y. Molecular properties and involvement of heparanase in cancer metastasis and angiogenesis. J Clin Invest. 2001;108:341–7.
Taniguchi N, Kizuka Y. Glycans and cancer: role of N-glycans in cancer biomarker, progression and metastasis, and therapeutics. Adv Cancer Res. 2015;126:11–51.
Caldwell SA, Jackson SR, Shahriari KS, Lynch TP, Sethi G, Walker S, et al. Nutrient sensor O-GlcNAc transferase regulates breast cancer tumorigenesis through targeting of the oncogenic transcription factor FoxM1. Oncogene. 2010;29:2831–42.
Ma Z, Vosseller K. Cancer metabolism and elevated O-GlcNAc in oncogenic signaling. J Biol Chem. 2014;289:34457–65.
Zhu W, Leber B, Andrews DW. Cytoplasmic O-glycosylation prevents cell surface transport of E-cadherin during apoptosis. EMBO J. 2001;20:5999–6007.
Yang WH, Kim JE, Nam HW, Ju JW, Kim HS, Kim YS, et al. Modification of p53 with O-linked N-acetylglucosamine regulates p53 activity and stability. Nat Cell Biol. 2006;8:1074–83.
Tanaka H, Imamura N, Oguma N, Shintani T, Tanaka K, Hyodo H, et al. Acute myelogenous leukemia with PIG-A gene mutation evolved from aplastic anemia-paroxysmal nocturnal hemoglobinuria syndrome. Int J Hematol. 2001;73:206–12.
Kunigou O, Nagao H, Kawabata N, Ishidou Y, Nagano S, Maeda S, et al. Role of GOLPH3 and GOLPH3L in the proliferation of human rhabdomyosarcoma. Oncol Rep. 2011;26:1337–42.
Feng Y, He F, Yan S, Huang H, Huang Q, Deng T, et al. The role of GOLPH3L in the prognosis and NACT response in cervical cancer. J Cancer. 2017;8:443–54.
He S, Niu G, Shang J, Deng Y, Wan Z, Zhang C, et al. The oncogenic Golgi phosphoprotein 3 like overexpression is associated with cisplatin resistance in ovarian carcinoma and activating the NF-kappaB signaling pathway. J Exp Clin Cancer Res. 2017;36:137.
Kim J, Yu L, Chen W, Xu Y, Wu M, Todorova D, et al. Wild-type p53 promotes cancer metabolic switch by inducing PUMA-dependent suppression of oxidative phosphorylation. Cancer Cell. 2019;35:191–203.e8.
Kim J, Xu S, Xiong L, Yu L, Fu X, Xu Y. SALL4 promotes glycolysis and chromatin remodeling via modulating HP1alpha-Glut1 pathway. Oncogene. 2017;36:6472–9.
Wu DM, Wang S, Wen X, Han XR, Wang YJ, Fan SH, et al. MircoRNA-1275 promotes proliferation, invasion and migration of glioma cells via SERPINE1. J Cell Mol Med. 2018;22:4963–74.
Martinez-Ramos D, Simon-Monterde L, Queralt-Martin R, Suelves-Piqueres C, Menor-Duran P, Escrig-Sos J. Breast cancer in octogenarian. Are we doing our best? A population-registry based study. Breast. 2018;38:81–5.
Tang J, Pearce L, O'Donnell-Tormey J, Hubbard-Lucey VM. Trends in the global immuno-oncology landscape. Nat Rev Drug Discov. 2018;17:922.
Sotgia F, Whitaker-Menezes D, Martinez-Outschoorn UE, Salem AF, Tsirigos A, Lamb R, et al. Mitochondria "fuel" breast cancer metabolism: fifteen markers of mitochondrial biogenesis label epithelial cancer cells, but are excluded from adjacent stromal cells. Cell Cycle. 2012;11:4390–401.
Zhang W, Xu J, Fang H, Tang L, Chen W, Sun Q, et al. Endothelial cells promote triple-negative breast cancer cell metastasis via PAI-1 and CCL5 signaling. FASEB J. 2018;32:276–88.
Humphries BA, Buschhaus JM, Chen YC, Haley HR, Qyli T, Chiang B, et al. Plasminogen activator inhibitor 1 (PAI1) promotes actin cytoskeleton reorganization and glycolytic metabolism in triple-negative breast cancer. Mol Cancer Res. 2019;17:1142–54.
Azimi I, Petersen RM, Thompson EW, Roberts-Thomson SJ, Monteith GR. Hypoxia-induced reactive oxygen species mediate N-cadherin and SERPINE1 expression, EGFR signalling and motility in MDA-MB-468 breast cancer cells. Sci Rep. 2017;7:15140.
Szoltysek K, Janus P, Zajac G, Stokowy T, Walaszczyk A, Widlak W, et al. RRAD, IL4I1, CDKN1A, and SERPINE1 genes are potentially co-regulated by NF-kappaB and p53 transcription factors in cells exposed to high doses of ionizing radiation. BMC Genomics. 2018;19:813.
Chakraborty AK, Liang K, DiGiovanna MP. Co-targeting insulin-like growth factor I receptor and HER2: dramatic effects of HER2 inhibitors on nonoverexpressing breast cancer. Cancer Res. 2008;68:1538.
Lategahn J, Hardick J, Grabe T, Niggenaber J, Jeyakumar K, Keul M, et al. Targeting Her2-insYVMA with covalent inhibitors-a focused compound screening and structure-based design approach. J Med Chem. 2020;63:11725–55.
Kennedy SP, O'Neill M, Cunningham D, Morris PG, Toomey S, Blanco-Aparicio C, et al. Preclinical evaluation of a novel triple-acting PIM/PI3K/mTOR inhibitor, IBL-302, in breast cancer. Oncogene. 2020;39:3028–40.
Wang S, Wang N, Zheng Y, Yang B, Liu P, Zhang F, et al. Caveolin-1 inhibits breast cancer stem cells via c-Myc-mediated metabolic reprogramming. Cell Death Dis. 2020;11:450.
Gandhi N, Das GM. Metabolic reprogramming in breast cancer and its therapeutic implications. Cells. 2019;8:89.
This study is supported by the National Natural Science Foundation of China (Grant No. 81902865 and Grant No. 61427807), the China Postdoctoral Science Foundation (2019 M652974), Sanming Project of Medicine in Shenzhen (SZSM201612019), the Science and Technology Foundation of Jiangmen (JK  No. 58-2020YLI010), the Medical Science and Technology Research Foundation of Guangdong Province (YWJKH  No. 15-A2020510), and the President’s Fund of Nanfang Hospital, Southern Medical University (Grant No. 2019C017).
Ethics approval and consent to participate
All institutional and national guidelines for the care and use of laboratory animals were followed.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Figure 1. GOLPH3L regulated the tumorigenic activities of BT474 cells. (A) Knockdown of GOLPH3L in MCF-10A, BT474 and T47D cells with five distinct siRNAs. (B) The overexpression (oe) of GOLPH3L in MCF-10A, BT474 and T47D cells. (C) The expression of GOLPH3L in BT474 cells promotes their proliferation, n = 3, * p < 0.05, ** p < 0.01. (D) GOLPH3L expression suppresses the apoptosis of BT474 cells, n = 3, ** p < 0.01. (E) Knockdown of GOLPH3L inhibits the migration of BT474 cells using a transwell assay, n = 3, * p < 0.05, ** p < 0.01. (F) The suppression of GOLPH3L significantly blocks the BT474 cell cycle at G0/ G1 phase, n = 3, * p < 0.05, ** p < 0.01. Supplementary Figure 2. miR-1185-2-3p negatively regulates the tumorigenesis of BT474 cells. (A) The predicted miRNAs which were most likely to regulate with GOLPH3L. (B) The overexpression of miR-1185-2-3p was achieved with miRNA mimics and inhibition miR-1185-2-3p was achieved with miRNA inhibitor in T47D and BT474, n = 3, ** p < 0.01. (C) The overexpression of miR-1185-2-3p inversely correlated with the proliferation of BT474 cells, n = 3, ** p < 0.01. (D) The overexpression of miR-1185-2-3p promoted the apoptosis of BT474 cells, n = 3, * p < 0.05. (E) miR-1185-2-3p overexpression could significantly inhibited the migration of BT474 cells, n = 3, * p < 0.05. (F) The overexpression of miR-1185-2-3p could block the cell cycle at G0/ G1 phase, n = 3, * p < 0.05. Supplementary Figure 3. miR-1185-2-3p overexpression partially reversed the tumorigenesis induced by GOLPH3L overexpression. (A) The overexpression of miR-1185-2-3p inversely correlated with the up-regulation of proliferation induced by GOLPH3L overexpression in T47D cells, n = 3, * p < 0.05, ** p < 0.01. (B) miR-1185-2-3p partially reversed the apoptosis induced by GOLPH3L overexpression in T47D cells, n = 3. * p < 0.05. (C) Overexpressed miR-1185-2-3p inhibited the up-regulation of migration caused by GOLPH3L up-regulation using transwell assay, n = 3. * p < 0.05. (D) The overexpression of miR-1185-2-3p affected the cell cycle regulation by GOLPH3L, n = 3, * p < 0.05, ** p < 0.01.Supplementary Figure 4. p53 pathway promotes the transcription of SERPINE1. (A) Volcano plots of fold changes and p-values of altered metabolites after GOLPH3L suppression in the Positive Ion Mode (top panel) and Negative Ion Mode (bottom panel). (B) Volcano plots of fold changes and p-values of RNA-seq data after silencing GOLPH3L expression in T47D cells. (C) STRING analysis predicted protein-protein network indicating that SERPINE1 interacted with p53 (top panel) and the role of SERPINE1 in p53 signaling pathway after GOLPH3L knockdown in T47D cells (bottom panel).
About this article
Cite this article
Xu, Y., Chen, W., Liang, J. et al. The miR-1185-2-3p—GOLPH3L pathway promotes glucose metabolism in breast cancer by stabilizing p53-induced SERPINE1. J Exp Clin Cancer Res 40, 47 (2021). https://doi.org/10.1186/s13046-020-01767-9
- Glucose metabolism
- p53-induced transcription