Skip to main content

Construction and investigation of β3GNT2-associated regulatory network in esophageal carcinoma

Abstract

Background

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.

Methods

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.

Results

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.

Conclusions

Our study, for the first time, reveals the importance of β3GNT2 in ESCA progression and offers a potential therapeutic target for ESCA.

Background

Esophageal carcinoma (ESCA) is one of the most common cancers around the world [1]. China is a high incidence area for ESCA and approximately 50% of ESCA cases worldwide occur in China [2]. 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% [3]. 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 [8]. 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 [11]. 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/).

Tissue collection

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 [11]. 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).

Immunohistochemistry (IHC)

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 [12]. β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.

Luciferase assay

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 [13]. Detection was performed according to the manufacturer’s protocol. The binding of CREB1 to the β3GNT2 promoter region was evaluated by qPCR.

Animal experiments

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.

Statistical analysis

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.

Results

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.

Fig. 1
figure 1

Screening of significantly dysregulated glycosyltransferases in ESCA based on the TCGA and GTEx databases. A Flowchart of screening and identification process. B Volcano plot showing the distribution of differentially expressed glycosyltransferases. C Heatmap showing differential expression patterns of glycosyltransferases. D Survival plots for β3GNT2, ST3GAL4, β4GALT2, PIGA, and GYG1. E Pan-cancer analysis of β3GNT2 expression in paired tumor tissues and normal tissues. F Pan-cancer analysis of β3GNT2 expression in unpaired tumor and normal tissues. Data are expressed as mean ± SD. NS no significance; *P < 0.05; **P < 0.01; ***P < 0.001

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.

Table 1 Univariate and multivariate analysis of overall survival (Cox regression model)
Fig. 2
figure 2

High expression of β3GNT2 in ESCA predicts poor prognosis. A Forest plot based on multivariate Cox analysis. B A nomogram predicting 1-, 3-, 5-year survival rate. C The diagnostic value analysis of β3GNT2 by ROC curve. D Analysis of β3GNT2 mRNA expression by qPCR. E Analysis of β3GNT2 protein expression by IHC (bar value = 100 μm). F Kaplan–Meier analysis of overall survival according to β3GNT2 expression. Data are expressed as mean ± SD. **P < 0.01

Table 2 Relationship between β3GNT2 expression and clinicopathological features of ESCA patients

β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.

Fig. 3
figure 3

β3GNT2 promotes ESCA cell growth, migration, and invasion in vitro. A Analysis of β3GNT2 mRNA expression by qPCR. B Analysis of β3GNT2 protein expression by Western blot. C Cell proliferation was determined by the CCK-8 assay. D Colony formation ability was measured by colony formation assay. E Cell apoptosis was assessed by flow cytometry. F Cell migratory potential was analyzed by Transwell migration assay. G Cell invasion capacity was detected by Matrigel invasion assay. Cells were transfected with the indicated plasmids. shNC, negative control shRNA; shRNAs, β3GNT2 shRNA lentiviral vector; Mock, empty plasmid; OV, pcDNA3.1/β3GNT2 plasmid. Data are expressed as mean ± SD. *P < 0.05; **P < 0.01; ***P < 0.001

β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.

Fig. 4
figure 4

β3GNT2 facilitates ESCA growth in vivo. A Representative images of the tumors (n = 5 per group). B Measurement of tumor weight in the xenograft model. C Quantitative analysis of tumor volume. Cells were transfected with the indicated plasmids. shNC, negative control shRNA; shRNA, β3GNT2 shRNA2 lentiviral vector; Mock, empty plasmid; OV, pcDNA3.1/β3GNT2 plasmid. Data are expressed as mean ± SD. *P < 0.05

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.

Fig. 5
figure 5

EGFR can be glycosylated by β3GNT2. A GSEA plot showing a positive correlation between β3GNT2 expression and EGFR signaling. B EGFR pulled down by LEL was analyzed with an anti-EGFR antibody. IB, immunoblot; PD, pull-down. C Protein–protein interaction networks were constructed by the online GeneMANIA tool. D Protein–protein interaction networks were predicted by the String database. Cells were transfected with the indicated plasmids. shNC, negative control shRNA; shRNA, β3GNT2 shRNA2 lentiviral vector; Mock, empty plasmid; OV, pcDNA3.1/β3GNT2 plasmid

β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.

Fig. 6
figure 6

The JAK-STAT pathway is regulated by β3GNT2 in ESCA. A GSEA plot showing a positive correlation between β3GNT2 expression and the JAK-STAT pathway. B Heatmap showing the Spearman correlation coefficient between β3GNT2 and the JAK-STAT pathway-related genes in ESCA. Asterisks represent the significance of the Spearman correlation. C Analysis of the JAK-STAT pathway-related protein expression by Western blot. DG The effect of ruxolitinib on cell proliferation, colony formation, migration, and invasion was determined by CCK-8 assay (D), colony formation assay (E), Transwell migration assay (F), and Matrigel invasion assay (G). Cells were transfected with the indicated plasmids. shNC, negative control shRNA; shRNA, β3GNT2 shRNA2 lentiviral vector; Mock, empty plasmid; OV, pcDNA3.1/β3GNT2 plasmid. Data are expressed as mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001

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.

Fig. 7
figure 7

β3GNT2 is a downstream target of CREB1. A Transcription factors (TFs) of β3GNT2 predicted by KnockTF. B Correlation analysis between transcriptional activators and β3GNT2 based on the TCGA and GTEx databases. C Analysis of β3GNT2 expression through the TCGA and GTEx databases. D Correlation analysis of CREB1 mRNA expression and β3GNT2 mRNA expression in clinical ESCA samples. E Analysis of β3GNT2 mRNA expression by qPCR. F Analysis of β3GNT2 protein expression by Western blot. G CREB1 binding site in the β3GNT2 promoter region. H ChIP-qPCR analysis of CREB1 binding to the β3GNT2 promoter. I Dual-luciferase reporter activity assay. JM The proliferation, colony formation, migration, and invasion of cells co-transfected with β3GNT2 overexpression plasmid and CREB1 shRNA were measured by CCK-8 assay (J), colony formation assay (K), Transwell migration assay (L), and Matrigel invasion assay (M). Cells were transfected with the indicated plasmids. WT, wild-type; MUT, mutant; OV, pcDNA3.1/β3GNT2 plasmid; sh-Ctrl, control shRNA; sh-CREB1, CREB1 shRNA lentiviral vector. Data are expressed as mean ± SD. NA, Data not available; *P < 0.05; **P < 0.01; ***P < 0.001

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 [14]. 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.

Fig. 8
figure 8

β3GNT2 expression is controlled by miR-133b. A Prediction of regulatory miRNAs of β3GNT2. B Heatmap showing the top 50 miRNAs negatively correlated with β3GNT2 in ESCA. C Analysis of the correlation between miR-133b and β3GNT2 expression based on the TCGA database. D Correlation analysis of miR-133b and β3GNT2 expression in clinical ESCA samples. E Analysis of β3GNT2 mRNA expression by qPCR. F Analysis of β3GNT2 protein expression by Western blot. G Putative miR-133b binding site in the β3GNT2 3′-UTR. H Dual-luciferase reporter activity assay. IL The proliferation, colony formation, migration, and invasion of cells co-transfected with β3GNT2 overexpression plasmid and miR-133b mimics were measured by CCK-8 assay (I), colony formation assay (J), Transwell migration assay (K), and Matrigel invasion assay (L). M Diagram of the proposed mechanism showing how β3GNT2 regulates ESCA progression. OV, cells transfected with pcDNA3.1/β3GNT2 plasmid. Data are expressed as mean ± SD. #P > 0.05, *P < 0.05; **P < 0.01

Discussion

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 [15]. Abnormal expression of β3GNTs alters the length of polylactosamine chains, which play a critical role in carcinogenesis and cancer progression [16]. β3GNT2, as a member of the β3GNTs, has strong polylactosamine synthesizing activity and is mainly responsible for the elongation of polylactosamine [17]. 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 [23]. 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 [24]. EGFR could be glycosylated with various glycans catalyzed by different glycosyltransferases. For example, OST inhibition disrupted EGFR N-linked glycosylation [25]. Downregulation of GALNT2 suppressed the O-glycosylation of EGFR [26]. The sialylation status of EGFR was related to ST6Gal-I [27]. 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.

Conclusions

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.

Abbreviations

ESCA:

Esophageal carcinoma

β3GNT2:

β1, 3-N-Acetylglucosaminyltransferase 2

TCGA:

The Cancer Genome Atlas

GTEx:

Genotype-tissue expression

ROC:

Receiver operating characteristic

GSEA:

Gene set enrichment analysis

LEL:

Lycopersicon esculentum lectin

IHC:

Immunohistochemistry

qPCR:

Quantitative real-time PCR

CCK-8:

Cell counting Kit-8

ChIP:

Chromatin immunoprecipitation

References

  1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70:7–30.

    Article  Google Scholar 

  2. 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.

    Article  Google Scholar 

  3. 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.

    Article  CAS  Google Scholar 

  4. 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.

    Article  CAS  Google Scholar 

  5. 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.

    Article  CAS  Google Scholar 

  6. Silsirivanit A. Glycosylation markers in cancer. Adv Clin Chem. 2019;89:189–213.

    Article  CAS  Google Scholar 

  7. 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.

    Article  Google Scholar 

  8. 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.

    Article  CAS  Google Scholar 

  9. 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.

    Article  CAS  Google Scholar 

  10. 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.

    Article  Google Scholar 

  11. 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.

    Article  CAS  Google Scholar 

  12. 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.

    Article  Google Scholar 

  13. 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.

    Article  CAS  Google Scholar 

  14. Homayoonfal M, Asemi Z, Yousefi B. Targeting microRNAs with thymoquinone: a new approach for cancer therapy. Cell Mol Biol Lett. 2021;26:43.

    Article  Google Scholar 

  15. Narimatsu H. Human glycogene cloning: focus on beta 3-glycosyltransferase and beta 4-glycosyltransferase families. Curr Opin Struct Biol. 2006;16:567–75.

    Article  CAS  Google Scholar 

  16. 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.

    Article  CAS  Google Scholar 

  17. 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.

    Article  CAS  Google Scholar 

  18. 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.

    Article  CAS  Google Scholar 

  19. Henion TR, Madany PA, Faden AA, Schwarting GA. β3GnT2 null mice exhibit defective accessory olfactory bulb innervation. Mol Cell Neurosci. 2013;52:73–86.

    Article  CAS  Google Scholar 

  20. 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.

    Article  CAS  Google Scholar 

  21. 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.

    Article  CAS  Google Scholar 

  22. Owen KL, Brockwell NK, Parker BS. JAK-STAT signaling: a double-edged sword of immune regulation and cancer progression. Cancers (Basel). 2019;11:2002.

    Article  CAS  Google Scholar 

  23. 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.

    Article  CAS  Google Scholar 

  24. Azimzadeh Irani M, Kannan S, Verma C. Role of N-glycosylation in EGFR ectodomain ligand binding. Proteins. 2017;85:1529–49.

    Article  CAS  Google Scholar 

  25. 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.

    Article  Google Scholar 

  26. 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.

    Article  CAS  Google Scholar 

  27. 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.

    Article  Google Scholar 

  28. 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.

    Article  Google Scholar 

  29. 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.

    Article  CAS  Google Scholar 

  30. 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.

    Article  Google Scholar 

  31. 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.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

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).

Author information

Authors and Affiliations

Authors

Contributions

LS and QH designed the experiments. YT, YL, TT, and ST conducted the functional experiments. CH and AL analyzed the data. XD wrote the manuscript. ZL supervised the research. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Xinzhou Deng or Li Shen.

Ethics declarations

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

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

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/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s11658-022-00306-y

Keywords