- Open Access
Construction and investigation of β3GNT2-associated regulatory network in esophageal carcinoma
Cellular & Molecular Biology Letters volume 27, Article number: 8 (2022)
Glycosyltransferases play a crucial role in various cancers. β1, 3-N-acetylglucosaminyltransferase 2, a polylactosamine synthase, is an important member of the glycosyltransferase family. However, the biological function and regulatory mechanism of β3GNT2 in esophageal carcinoma (ESCA) is still poorly understood.
The Cancer Genome Atlas and Genotype-Tissue Expression databases were used for gene expression and prognosis analysis. Quantitative real-time PCR, Western blot, and immunohistochemistry were performed to detect the expression of β3GNT2 in ESCA cell lines and tissues. In vitro assays and xenograft tumor models were utilized to evaluate the impact of β3GNT2 on ESCA progression. The downstream effectors and upstream regulators of β3GNT2 were predicted by online software and verified by functional experiments.
We found that β3GNT2 was highly expressed in ESCA tissues and positively correlated with poor prognosis in ESCA patients. β3GNT2 expression was closely associated with the tumor size, TNM stage, and overall survival of ESCA patients. Functionally, β3GNT2 promoted ESCA cell growth, migration, and invasion in vitro, as well as tumorigenesis in vivo. Mechanistically, β3GNT2 knockdown decreased the expression of the polylactosamine on EGFR. Knockdown of β3GNT2 also inhibited the JAK/STAT signaling pathway. Meanwhile, the JAK/STAT inhibitor could partly reverse the biological effects caused by β3GNT2 overexpression. Moreover, β3GNT2 expression was positively regulated by CREB1 and negatively regulated by miR-133b. Both CREB1 and miR-133b was involved in the β3GNT2-mediated ESCA progression.
Our study, for the first time, reveals the importance of β3GNT2 in ESCA progression and offers a potential therapeutic target for ESCA.
Esophageal carcinoma (ESCA) is one of the most common cancers around the world . China is a high incidence area for ESCA and approximately 50% of ESCA cases worldwide occur in China . Although great progress has been made in the treatment of ESCA, such as surgery, radiotherapy, and chemotherapy, the five-year survival rate of ESCA patients is still less than 20% . This is mainly due to the lack of effective diagnostic and prognostic biomarkers. Moreover, the mechanism involved in the occurrence and development of ESCA has not been fully elucidated. Thus, it is of great importance to explore the pathogenesis of ESCA and develop effective therapies for ESCA patients.
Glycosylation is a ubiquitous post-translational modification of proteins [4, 5]. Alteration of glycosylation, as a hallmark of cancer, is linked to tumor initiation, progression, and metastasis [6, 7]. Glycosylation is generated by a variety of glycosyltransferases through complex biosynthetic pathways. More than 200 glycosyltransferases have been reported so far . Dysregulated expression of glycosyltransferases has been implicated in ESCA. For example, high expression of C1GALT1 and OGT was associated with lymph node metastasis of ESCA patients [9, 10]. FUT8 was identified as a key driver for radioresistance in ESCA . However, a comprehensive understanding of glycosyltransferases in ESCA is still lacking.
In this study, The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases were used to identify differentially expressed glycosyltransferases between ESCA tissues and normal tissues. We found that β1, 3-N-acetylglucosaminyltransferase 2 (β3GNT2) was overexpressed in ESCA tissues and correlated with the poor prognosis of ESCA patients. We further investigated the biological function and regulatory mechanism of β3GNT2 in ESCA progression. Our results suggest that β3GNT2 may serve as a potential therapeutic target for ESCA.
Materials and methods
Data acquisition and processing
RNA-seq data, miRNA-seq data, and the corresponding clinical information were obtained from the TCGA (https://tcga-data.nci.nih.gov/tcga/) and GTEx (http://gtexportal.org/home/) platforms. Differentially expressed glycosyltransferases, transcription factors, miRNAs were screened using the DESeq2 R package. An adjusted P-value < 0.05 and |log2 fold-change (FC)|> 0.5 were set as the thresholds. Volcano plots and heatmaps were generated using the R package ggplot and pheatmap, respectively. A prognostic nomogram was constructed by R software with the rms package. The Receiver Operating Characteristic (ROC) curves were plotted and the area under the ROC curve was calculated using the R package pROC. Survival prediction was performed using the GEPIA database (http://gepia.cancer-pku.cn/). Gene Set Enrichment Analysis (GSEA) analysis was carried out using the GSEA software (http://www.broadinstitute.org/gsea). Co-expression analysis was determined by using the LinkedOmics platform (http://www.linkedomics.org/).
Sixty-five tumor tissues and paired paracancerous tissues were collected from patients who were diagnosed with ESCA and received surgery at the Taihe Hospital, Hubei University of Medicine (Shiyan, China) between 2008 and 2012. The patients, including 35 males and 30 females, were aged 35–75 years with an average age of 54.62 ± 9.14. None of the patients underwent any other treatment before surgery. Samples were stored at − 80 °C until further processing. The studies involving human participants were approved by the Ethics Committee of the Hubei University of Medicine. All patients provided informed written consent.
Cell culture and transfection
Human ESCA cell lines TE-1, KYSE150, and KYSE410 were purchased from Procell (Wuhan, China). These cells were cultured in RPMI-1640 medium (Gibco, Carlsbad, CA, USA) containing 10% FBS at 37 °C, 5% CO2. The pcDNA3.1/β3GNT2 plasmid, pcDNA3.1/CREB1 plasmid, empty pcDNA3.1 plasmid, a lentiviral vector expressing β3GNT2 shRNA, a lentiviral vector expressing CREB1 shRNA, control lentiviral vector, miR-133b inhibitor, inhibitor control, miR-133b mimics, and mimics control were constructed by GenePharma (Shanghai, China). All transfections were conducted with Lipofectamine 3000 reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol. Stable clones were selected by puromycin or G418.
Quantitative real-time PCR (qPCR) and Western blot
Total RNA was extracted using RNAiso Plus (Takara, Dalian, China). For mRNA detection, cDNA was synthesized using the PrimeScript RT Reagent Kit (Takara) and qPCR was carried out employing the SYBR-Green PCR master mix (Takara). For miRNA analysis, the One-Step PrimeScript miRNA cDNA synthesis kit (Takara) and the SYBR PrimeScript miRNA RT-PCR Kit (Takara) were used. Data were normalized to GAPDH or U6 based on the 2−ΔΔCt method. All primers were synthesized by GenePharma. Total protein was isolated using RIPA lysis buffer (Beyotime, Shanghai, China). Western blot was conducted as previously described . Primary antibodies were as follows: β3GNT2 (ab236291), CREB1 (ab32515), JAK1 (ab133666), STAT3 (ab68153), p-STAT3 (ab267373), EGFR (ab52894), and GAPDH (ab8245). All antibodies were acquired from Abcam (Shanghai, China).
Tissues were fixed in formalin and embedded in paraffin. Afterward, 5 µm sections were prepared for IHC staining. PV-9000 two-step immunohistochemical kit (ZSGB-BIO, Beijing China) was used to detect the expression of β3GNT2 protein. The positive signal was visualized by DAB detection Kit (ZSGB-BIO). Slides were read by two pathologists using a Nikon Eclipse E600 microscope (Nikon, Tokyo, Japan). IHC score = intensity score × percentage score . β3GNT2 expression was defined as follows: low expression (score 0–4) and high expression (score ≥ 5).
Cell proliferation and colony formation assays
Cell proliferation was detected by Cell Counting Kit-8 (CCK-8) assay (Dojindo, Kumamoto, Japan). Cells (1 × 103 per well) were seeded into a 96-well plate. Absorbance was measured at 450 nm. For colony formation assay, cells were added into a 6-well plate at a density of 5 × 103 cells/well. After 2 weeks, colonies (> 50 cells) were stained with 0.5% crystal violet (Beyotime).
Cell migration and invasion assays
Cells (5 × 104 in 200 μL serum-free medium) were seeded into the upper chamber of the Transwell apparatus (Corning Incorporated, Corning, NY, USA). Transwell chambers precoated with or without Matrigel (BD Biosciences, San Jose, CA, USA) were used to assess cell invasion and migration, respectively. The lower chamber was filled with 600μL of complete culture medium. After incubation for 36 h, migrated or invaded cells were stained with 0.5% crystal violet.
Cell apoptosis assay
Cell apoptosis was detected by flow cytometry using the Apoptosis Detection Kit (Absin, Shanghai, China). The percentage of apoptotic cells was defined as the sum of early apoptotic cells and late apoptotic cells. All experimental steps were conducted according to kit manual instructions.
Lectin pull-down assay
Lycopersicon esculentum lectin (LEL) agarose beads (#AL-1113, Vector Labs, Burlingame, CA, USA) were used to analyze polylactosamine chains on glycoproteins. Cell lysates (0.5 mg) were incubated with 30 μL of LEL-conjugated beads to capture the lectin-glycoprotein complexes. Precipitated glycoproteins were eluted and subjected to Western blot analysis with an anti-EGFR antibody. The EGFR in total lysates was used as the internal loading control.
The wild-type or mutant β3GNT2 3′UTR sequences containing the binding sites for miR-133b were cloned into the pmirGLO vector. Besides, the fragments of β3GNT2 promoter with wild-type or mutant CERB1 binding sites were inserted into a pGL3-basic vector. All the above plasmids were constructed by GenePharma. Cells were harvested 48 h after transfection. Luciferase activities were determined using the dual-luciferase reporter assay system (Promega, Madison, WI, USA).
Chromatin immunoprecipitation (ChIP) assay
The EZ-Magna ChIP™ Kit (Millipore, Billerica, MA, USA) was used for ChIP assay . Detection was performed according to the manufacturer’s protocol. The binding of CREB1 to the β3GNT2 promoter region was evaluated by qPCR.
Female BALB/c nude mice, 4–5 weeks old, were provided by the Animal Center of Hubei University of Medicine. Transfected cells (1 × 107) were subcutaneously injected into nude mice. Tumor volume (0.5 × length × width2) was measured every week. Four weeks after injection, all mice were sacrificed.
Data are presented as mean ± SD. Results were analyzed with Graphpad Prism (version 7) software (GraphPad, San Diego, CA, USA). Student’s t-test, Chi-square test, and one-way ANOVA were used to determine statistical significance. A correlation was analyzed using the Pearson correlation coefficient test. Overall survival curves were plotted according to the Kaplan–Meier method. A P-value < 0.05 was considered significant.
Identification of differentially expressed glycosyltransferases in ESCA
Using the RNA-seq data from the TCGA and GTEx datasets, we screened differentially expressed glycosyltransferases between ESCA tissues (182 cases) and normal esophageal tissues (666 cases) (Fig. 1A). In total, there were 63 upregulated and 57 downregulated glycosyltransferases in ESCA tissues compared with normal tissues (Fig. 1B, Additional file 1: Table S1). Aberrant expression of glycosyltransferases was presented as a heatmap plot (Fig. 1C). We further investigated the relationship between glycosyltransferases and the overall survival of ESCA patients. Kaplan–Meier analysis showed that only five glycosyltransferases, including β3GNT2, ST3GAL4, β4GALT2, PIGA, and GYG1 were related to overall survival (Fig. 1D, Additional file 1: Table S1). Among them, we observed that β3GNT2 was the most significantly upregulated glycosyltransferase in ESCA (log2|FC|= 0.86591, P < 0.001). Thus, β3GNT2 was expected to be an oncogenic factor. To elucidate the gene expression characteristics of β3GNT2 in human malignancies, pan-cancer analysis was performed using the TCGA and GTEx databases. We found that β3GNT2 expression was increased in most types of cancers (Fig. 1E, F). Based on the above analysis, β3GNT2 was selected for further study.
Clinical importance of β3GNT2 in ESCA
To evaluate the prognostic potential of β3GNT2, we downloaded the clinical information of 182 ESCA patients from the TCGA data portal. Analysis was performed using the univariate or multivariate Cox proportional hazards regression (Table 1). The univariate analysis indicated that N stage, M stage, pathologic stage, and β3GNT2 expression could influence overall survival. The multivariate analysis revealed that β3GNT2 was an independent prognostic factor for overall survival (Fig. 2A). Based on the multivariate analysis, a nomogram prognostic evaluation model was successfully constructed, which could predict the probability of 1-, 3- and 5-year overall survival of patients with ESCA (Fig. 2B). The ROC curve showed that β3GNT2 had a relatively high diagnostic accuracy, representing its strong potential for distinguishing normal samples from ESCA samples (Fig. 2C). To validate the findings obtained from public databases, we collected 65 ESCA specimens. The results of qPCR and IHC demonstrated that β3GNT2 was upregulated in ESCA tissues (Fig. 2D, E). According to Kaplan–Meier analysis, high β3GNT2 expression was linked to unsatisfactory survival in patients with ESCA (Fig. 2F). By analyzing the pathological features, we found that β3GNT2 expression was associated with the tumor size and TNM stage (Table 2). Thus, β3GNT2 could be a valuable prognostic marker for ESCA, and upregulation of β3GNT2 might contribute to ESCA development.
β3GNT2 is required for ESCA cell growth, migration, and invasion in vitro
To understand the role of β3GNT2 in ESCA, we first analyzed the protein expression of β3GNT2 in three human ESCA cell lines. We found that β3GNT2 expression was highest in TE-1 cells, whereas KYSE150 and KYSE410 cells showed moderate or lower levels of β3GNT2, respectively (Additional file 2: Figure S1). To determine whether β3GNT2 affected cell growth, we conducted loss-of-function studies in TE-1 and KYSE150 cells. The efficiency of β3GNT2 knockdown was verified by qPCR and western blot (Fig. 3A, B). CCK-8 and colony formation assays showed that β3GNT2 knockdown significantly inhibited the proliferation and colony-forming abilities of TE-1 and KYSE150 cells (Fig. 3C, D). To test whether the β3GNT2 knockdown-mediated cell growth inhibition was associated with apoptosis, flow cytometry was carried out. We observed that the proportion of apoptosis in TE-1 and KYSE150 cells was increased following β3GNT2 knockdown (Fig. 3E). Cell migration and invasion are critical steps in tumor metastasis. Transwell migration and Matrigel invasion assays revealed that the migrative and invasive capacities of TE-1 and KYSE150 cells were remarkably decreased after β3GNT2 knockdown (Fig. 3F, G). Accordingly, β3GNT2 shRNA2 was used for subsequent experiments. To validate the effect of β3GNT2 in ESCA progression, we employed the gain-of-function strategy in KYSE410 cells (Fig. 3A, B). Functionally, β3GNT2 overexpression promoted proliferation (Fig. 3C), colony formation (Fig. 3D), migration (Fig. 3F), and invasion (Fig. 3G) of KYSE410 cells. These findings suggested that β3GNT2 exhibited strong oncogenic properties in ESCA.
β3GNT2 confers in vivo tumorigenicity
Subcutaneous xenograft models were used to explore the function of β3GNT2 in vivo. TE-1 cells with or without the knockdown of β3GNT2 were subcutaneously injected into nude mice. We discovered that β3GNT2 knockdown suppressed tumor growth, as manifested by reduced tumor size (Fig. 4A), tumor weight (Fig. 4B), and tumor volume (Fig. 4C). In contrast, the opposite results were observed in mice injected with KYSE410 cells overexpressing β3GNT2 (Fig. 4A–C). These data implicated that β3GNT2 played an essential role in ESCA tumorigenesis.
EGFR glycosylation is regulated by β3GNT2
Our pathway analysis from GSEA showed that there was a positive correlation between β3GNT2 expression and EGFR signaling in ESCA (Fig. 5A). β3GNT2 has polylactosamine synthesizing activity and possesses the ability to synthesize polylactosamine structures. To investigate whether EGFR was decorated with polylactosamine and whether C1GALT1 could influence EGFR glycosylation, a lectin pull-down assay was performed. We found that EGFR could be pulled down by LEL and β3GNT2 knockdown reduced the LEL binding to EGFR. In contrast, the overexpression of β3GNT2 had the opposite effects, indicating that polylactosamine chains on EGFR were modulated by β3GNT2 (Fig. 5B). Furthermore, a gene regulatory network was constructed using the GeneMANIA website (http://genemania.org/) and the String database (https://www.string-db.org/) to determine the interactive relationship between β3GNT2 and EGFR (Fig. 5C, D). Hence, EGFR was a major downstream effector of β3GNT2 in ESCA.
β3GNT2 promotes ESCA progression through the JAK/STAT pathway
To explore the mechanism underlying the role of β3GNT2 in ESCA progression, genes co-expressed with β3GNT2 in ESCA were analyzed by using the LinkedOmics platform (Additional file 3: Table S2). These genes were then subjected to GSEA analysis. We found that β3GNT2 was closely related to the JAK/STAT pathway (Fig. 6A). The correlation between β3GNT2 and the JAK-STAT pathway-related genes was displayed by the heatmap (Fig. 6B). We discovered that the major proteins of the JAK-STAT pathway such as JAK1 and p-STAT3 were decreased after β3GNT2 knockdown but increased upon β3GNT2 overexpression (Fig. 6C). Then β3GNT2-overexpressing cells were pretreated with the JAK inhibitor, ruxolitinib (1 μM). The CCK-8, colony formation, Transwell migration, and Matrigel invasion assays showed that the promotive effects of β3GNT2 on cell proliferation, colony formation, migration, and invasion were abrogated by ruxolitinib (Fig. 6D–G). These results indicated that the oncogenic potential of β3GNT2 in ESCA was linked to the JAK/STAT pathway.
CREB1 activates β3GNT2 transcription in ESCA
To clarify the mechanism leading to β3GNT2 upregulation in ESCA, the potential transcription factors involved in the regulation of β3GNT2 were predicted using the online tool KnockTF (http://www.licpathway.net/KnockTF/index.html) (Fig. 7A). The association between transcriptional activators and β3GNT2 in ESCA was analyzed based on the TCGA and GTEx databases (Fig. 7B). Notably, CREB1 was most strongly correlated with β3GNT2 (r = 0.315). Moreover, CREB1 was overexpressed in ESCA (Fig. 7C). A positive correlation was also observed between CREB1 and β3GNT2 in our clinical samples (Fig. 7D). Thus, we speculated that CREB1 could act as an upstream regulator of β3GNT2. To test this hypothesis, we investigated whether β3GNT2 expression was regulated by CREB1. We demonstrated that silencing of CREB1 downregulated β3GNT2 mRNA and protein expression in TE-1 and KYSE150 cells, whereas CREB1 overexpression upregulated the mRNA and protein expression of β3GNT2 in KYSE410 cells (Fig. 7E, F). To identify whether CREB1 could bind to the promoter region of β3GNT2, the JASPAR database (http://jaspar.genereg.net/) was used. Through bioinformatics analysis, one putative CREB1 binding site was identified in the β3GNT2 promoter region (Fig. 7G). The enrichment of CREB1 in the β3GNT2 promoter region was validated by ChIP-qPCR (Fig. 7H). Dual-luciferase reporter assay showed that β3GNT2 transcriptional activity was increased after transfection with CREB1 plasmid. However, mutation of the putative CREB1 binding site attenuated the β3GNT2 promoter activity (Fig. 7I). To explore whether CREB1 participated in the β3GNT2-mediated ESCA progression, we conducted functional rescued experiments. We found that CREB1 silencing inhibited the proliferation, colony formation, migration, and invasion of TE-1 and KYSE150 cells, but this suppression was partially reversed by β3GNT2 overexpression (Fig. 7J–M). In summary, these results suggested that CREB1 transcriptionally upregulated β3GNT2 expression in ESCA.
miR-133b is another regulatory factor of β3GNT2
Since miRNAs are important post-transcriptional modulators of gene expression, we then investigated whether β3GNT2 was regulated by specific miRNAs . Three publicly available algorithms (TargetScan, miRWalK, and miRDB) were used to predict the potential miRNAs targeting β3GNT2. We found that β3GNT2 was the predicted target of miR-200c-3p and miR-133b (Fig. 8A). In general, miRNA expression is inversely correlated with the expression of target genes. Subsequently, the top 50 negatively related miRNAs for β3GNT2 were identified in the dataset of TCGA-ESCA (Fig. 8B). Remarkably, only miR-133b was downregulated in ESCA, and its expression was negatively associated with β3GNT2 (Fig. 8C). A negative correlation between miR-133b and β3GNT2 was also observed in our clinical samples (Fig. 8D). We further determined whether β3GNT2 expression could be modulated by miR-133b. As expected, miR-133b mimics suppressed β3GNT2 expression both at the mRNA and protein levels in TE-1 and KYSE150 cells, while miR-133b inhibitor enhanced β3GNT2 mRNA and protein expression in KYSE410 cells (Fig. 8E, F). Using the TargetScan algorithm, we discovered that the 3′-UTR of the β3GNT2 gene contained a potential miR-133b binding site (Fig. 8G). Moreover, miR-133b mimics reduced the luciferase activities in TE-1 and KYSE150 cells transfected with wild-type β3GNT2 3′-UTR but not with mutant β3GNT2 3′-UTR (Fig. 8H). Functional experiments demonstrated that the inhibitory effect of miR-133b mimics on cell proliferation, colony formation, migration, and invasion could be restored by β3GNT2 overexpression (Fig. 8I–L). These results suggested that miR-133b targeted β3GNT2 and inhibited its expression in ESCA.
The prognosis of patients with ESCA remains unsatisfactory, despite the use of multiple therapeutic approaches. Further exploration of the diagnostic and prognostic biomarkers for ESCA is warranted. Glycosyltransferases have attracted increasing interest owing to their importance in many biological processes and diseases, including cancer. In the current study, we first found that β3GNT2 expression was upregulated in ESCA. β3GNT2 upregulation was closely associated with the poor prognosis of ESCA patients. Subsequent experiments proved that β3GNT2 contributed to ESCA progression by promoting EGFR glycosylation and activating the JAK-STAT pathway. In addition, CREB1 and miR-133b were critical regulators of β3GNT2 (Fig. 8M). Thus, β3GNT2 may function as an oncogene in ESCA and may be a promising prognostic biomarker for ESCA.
β3GNTs are a family of crucial glycosyltransferases that synthesize a unique glycan structure known as polylactosamine . Abnormal expression of β3GNTs alters the length of polylactosamine chains, which play a critical role in carcinogenesis and cancer progression . β3GNT2, as a member of the β3GNTs, has strong polylactosamine synthesizing activity and is mainly responsible for the elongation of polylactosamine . Dysregulated β3GNT2 expression has been reported in several diseases [18,19,20]. Our analysis of β3GNT2 expression in clinical samples as well as from public databases both supported that β3GNT2 played an important role in the occurrence and development of ESCA. β3GNT2 knockdown inhibited ESCA growth, migration, and invasion in vitro, as well as tumor formation in vivo. Additionally, ectopic expression of β3GNT2 had the opposite biological function. The present study investigated, for the first time, the oncogenic activity of β3GNT2 in ESCA.
We explored the underlying mechanism of β3GNT2 in promoting ESCA progression. Bioinformatics analysis showed that β3GNT2 expression was linked to the JAK-STAT pathway. Numerous studies revealed that the JAK-STAT pathway was involved in cell proliferation, differentiation, apoptosis, angiogenesis, and immune regulation [21, 22]. The JAK-STAT pathway was also strongly associated with the occurrence of ESCA . In this study, we demonstrated that β3GNT2 promoted ESCA cell proliferation, invasion, and migration via the activation of the JAK-STAT pathway. More importantly, our study further identified EGFR as a downstream effector of β3GNT2 in ESCA. EGFR is a transmembrane glycoprotein and belongs to the tyrosine kinase receptor family . EGFR could be glycosylated with various glycans catalyzed by different glycosyltransferases. For example, OST inhibition disrupted EGFR N-linked glycosylation . Downregulation of GALNT2 suppressed the O-glycosylation of EGFR . The sialylation status of EGFR was related to ST6Gal-I . We confirmed that EGFR carried polylactosamine-type glycans and was a substrate of β3GNT2. These results suggested that β3GNT2 exerted its biological function in ESCA progression mainly through modulating the JAK-STAT pathway and EGFR glycosylation.
To explore the potential factors leading to β3GNT2 upregulation in ESCA, the transcription factors of β3GNT2 were predicted by bioinformatics methods. CREB1 was proposed as a candidate because it had a strong positive correlation with β3GNT2 expression in ESCA. CREB1 is a transcription factor that influences multiple cellular processes, such as cell survival, differentiation, and proliferation [28, 29]. In the present study, we confirmed that CREB1 could bind to the promoter region of β3GNT2 and then enhance its expression. Functional rescued experiments demonstrated that CREB1 silencing suppressed ESCA cell growth, migration, and invasion, but this suppression was reversed by β3GNT2 overexpression. Hence, CREB1 was a key regulator of β3GNT2 in ESCA. To further clarify the upstream regulation mechanism of β3GNT2, the regulatory miRNAs of β3GNT2 were analyzed using the online prediction tools. We discovered that miR-133b could interact with β3GNT2 by binding to its 3′-UTR and there was a negative association between miR-133b and β3GNT2 expression in ESCA. The expression of β3GNT2 was decreased by miR-133b mimics but increased by miR-133b inhibitor. Moreover, the inhibitory effect of miR-133b mimics on ESCA cell growth, migration, and invasion was recovered by β3GNT2 overexpression. It has been reported that miR-133b participates in tumor proliferation and metastasis [30, 31]. Therefore, miR-133b was another critical regulatory factor of β3GNT2 in ESCA.
There are, however, several limitations to our study. First, the expression of polylactosamine chains in ESCA tissues is not examined. Second, the association between the JAK-STAT pathway and EGFR signaling has not been elucidated. The specific glycosylation sites of EGFR remain unknown. The effect of EGFR glycosylation on ESCA progression has not yet been elucidated. Additionally, the mechanism of synergy between CREB1 and miR-133b is not clearly understood. It will be interesting to address these questions in future research.
Our results confirm that β3GNT2 predicts poor prognosis and acts as an oncogene in ESCA. β3GNT2 can enhance EGFR glycosylation and activate the JAK-STAT pathway. β3GNT2 is transcriptionally upregulated by CREB1. Meanwhile, miR-133b directly targets β3GNT2 and negatively regulates its expression. To our knowledge, this is the first study to explore the biological function, molecular mechanism, and regulatory network of β3GNT2 in ESCA. Thus, β3GNT2 may be a promising prognosis biomarker and a potential therapeutic target for ESCA.
Availability of data and materials
All remaining data are availability within the article and Additional files, or available from the authors upon request.
β1, 3-N-Acetylglucosaminyltransferase 2
The Cancer Genome Atlas
Receiver operating characteristic
Gene set enrichment analysis
Lycopersicon esculentum lectin
Quantitative real-time PCR
Cell counting Kit-8
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70:7–30.
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.
Ishikawa T, Yasuda T, Okayama T, Dohi O, Yoshida N, Kamada K, et al. Early administration of pegfilgrastim for esophageal cancer treated with docetaxel, cisplatin, and fluorouracil: a phase II study. Cancer Sci. 2019;110:3754–60.
Flynn RA, Pedram K, Malaker SA, Batista PJ, Smith BAH, Johnson AG, et al. Small RNAs are modified with N-glycans and displayed on the surface of living cells. Cell. 2021;184:3109–24.
Gupta R, Leon F, Thompson CM, Nimmakayala R, Karmakar S, Nallasamy P, et al. Global analysis of human glycosyltransferases reveals novel targets for pancreatic cancer pathogenesis. Br J Cancer. 2020;122:1661–72.
Silsirivanit A. Glycosylation markers in cancer. Adv Clin Chem. 2019;89:189–213.
Wang M, Zhu J, Lubman DM, Gao C. Aberrant glycosylation and cancer biomarker discovery: a promising and thorny journey. Clin Chem Lab Med. 2019;57:407–16.
Schjoldager KT, Narimatsu Y, Joshi HJ, Clausen H. Global view of human protein glycosylation pathways and functions. Nat Rev Mol Cell Biol. 2020;21:729–49.
Wang Y, Liao X, Ye Q, Huang L. Clinic implication of MUC1 O-glycosylation and C1GALT1 in esophagus squamous cell carcinoma. Sci China Life Sci. 2018;61:1389–95.
Qiao Z, Dang C, Zhou B, Li S, Zhang W, Jiang J, et al. O-linked N-acetylglucosamine transferase (OGT) is overexpressed and promotes O-linked protein glycosylation in esophageal squamous cell carcinoma. J Biomed Res. 2012;26:268–73.
Shen L, Xia M, Deng X, Ke Q, Zhang C, Peng F, et al. A lectin-based glycomic approach identifies FUT8 as a driver of radioresistance in oesophageal squamous cell carcinoma. Cell Oncol (Dordr). 2020;43:695–707.
Shen A, Liu L, Chen H, Qi F, Huang Y, Lin J, et al. Cell division cycle associated 5 promotes colorectal cancer progression by activating the ERK signaling pathway. Oncogenesis. 2019;8:19.
Tian T, Chen ZH, Zheng Z, Liu Y, Zhao Q, Qiu H, et al. Investigation of the role and mechanism of ARHGAP5-mediated colorectal cancer metastasis. Theranostics. 2020;10:5998–6010.
Homayoonfal M, Asemi Z, Yousefi B. Targeting microRNAs with thymoquinone: a new approach for cancer therapy. Cell Mol Biol Lett. 2021;26:43.
Narimatsu H. Human glycogene cloning: focus on beta 3-glycosyltransferase and beta 4-glycosyltransferase families. Curr Opin Struct Biol. 2006;16:567–75.
Lee PL, Kohler JJ, Pfeffer SR. Association of beta-1,3-N-acetylglucosaminyltransferase 1 and beta-1,4-galactosyltransferase 1, trans-Golgi enzymes involved in coupled poly-N-acetyllactosamine synthesis. Glycobiology. 2009;19:655–64.
Togayachi A, Kozono Y, Ishida H, Abe S, Suzuki N, Tsunoda Y, et al. Polylactosamine on glycoproteins influences basal levels of lymphocyte and macrophage activation. Proc Natl Acad Sci USA. 2007;104:15829–34.
Kataoka K, Huh NH. A novel beta1,3-N-acetylglucosaminyltransferase involved in invasion of cancer cells as assayed in vitro. Biochem Biophys Res Commun. 2002;294:843–8.
Henion TR, Madany PA, Faden AA, Schwarting GA. β3GnT2 null mice exhibit defective accessory olfactory bulb innervation. Mol Cell Neurosci. 2013;52:73–86.
Qiu H, Wu SL, Guo XH, Shen HJ, Zhang HP, Chen HL. Expression of β1,3-N-acetylglucosaminyltransferases during differentiation of human acute myeloid leukemia cells. Mol Cell Biochem. 2011;358:131–9.
Xin P, Xu X, Deng C, Liu S, Wang Y, Zhou X, et al. The role of JAK/STAT signaling pathway and its inhibitors in diseases. Int Immunopharmacol. 2020;80:106210.
Owen KL, Brockwell NK, Parker BS. JAK-STAT signaling: a double-edged sword of immune regulation and cancer progression. Cancers (Basel). 2019;11:2002.
You Z, Xu D, Ji J, Guo W, Zhu W, He J. JAK/STAT signal pathway activation promotes progression and survival of human oesophageal squamous cell carcinoma. Clin Transl Oncol. 2012;14:143–9.
Azimzadeh Irani M, Kannan S, Verma C. Role of N-glycosylation in EGFR ectodomain ligand binding. Proteins. 2017;85:1529–49.
Lopez Sambrooks C, Baro M, Quijano A, Narayan A, Cui W, Greninger P, et al. Oligosaccharyltransferase inhibition overcomes therapeutic resistance to EGFR tyrosine kinase inhibitors. Cancer Res. 2018;78:5094–106.
Sun Z, Xue H, Wei Y, Wang C, Yu R, Wang S, et al. Mucin O-glycosylating enzyme GALNT2 facilitates the malignant character of glioma by activating the EGFR/PI3K/Akt/mTOR axis. Clin Sci (Lond). 2019;133:1167–84.
Britain CM, Holdbrooks AT, Anderson JC, Willey CD, Bellis SL. Sialylation of EGFR by the ST6Gal-I sialyltransferase promotes EGFR activation and resistance to gefitinib-mediated cell death. J Ovarian Res. 2018;11:12.
Kant S, Craige SM, Chen K, Reif MM, Learnard H, Kelly M, et al. Neural JNK3 regulates blood flow recovery after hindlimb ischemia in mice via an Egr1/Creb1 axis. Nat Commun. 2019;10:4223.
Dai W, Xu Y, Mo S, Li Q, Yu J, Wang R, et al. GLUT3 induced by AMPK/CREB1 axis is key for withstanding energy stress and augments the efficacy of current colorectal cancer therapies. Signal Transduct Target Ther. 2020;5:177.
Wang QY, Zhou CX, Zhan MN, Tang J, Wang CL, Ma CN, et al. MiR-133b targets Sox9 to control pathogenesis and metastasis of breast cancer. Cell Death Dis. 2018;9:752.
Coccia E, Masanas M, López-Soriano J, Segura MF, Comella JX, Pérez-García MJ. FAIM is regulated by MiR-206, MiR-1–3p and MiR-133b. Front Cell Dev Biol. 2020;8:584606.
This work was supported by the Hubei Provincial Department of Science and Technology Innovation Group program (Grant No. 2019CFA034), the Free Exploration Foundation of the Hubei University of Medicine (Grant No. FDFR201802), the Research Project of the Health Commission of Hubei Province (Grant No. WJ2021M048), and the National Undergraduate Training Program for Innovation and Entrepreneurship (Grant No. 202113249002).
Ethics approval and consent to participate
The studies involving human participants were performed in accordance with the Declaration of Helsinki and were approved by the Ethics Committee of the Hubei University of Medicine (Approval No. 2018TH016, date: 2018.1.12). All patients provided informed written consent. The animal experiment protocol complied with the international guidelines and was approved by the Institutional Animal Care and Use Committee of the Hubei University of Medicine (Approval No. 2021015, date: 2021.3.16).
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.
Additional file 1: Table S1.
Screening of dysregulated glycosyltransferases in ESCA based on the TCGA and GTEx databases.
Additional file 2: Fig. S1.
Analysis of β3GNT2 protein expression in different ESCA cells by Western blot.
Additional file 3: Table S2.
Differentially expressed genes correlated with β3GNT2 in ESCA based on the LinkedOmics platform.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Luo, Z., Hu, Q., Tang, Y. et al. Construction and investigation of β3GNT2-associated regulatory network in esophageal carcinoma. Cell Mol Biol Lett 27, 8 (2022). https://doi.org/10.1186/s11658-022-00306-y
- Esophageal carcinoma