Introduction
Lung cancer is the most commonly diagnosed malignancy [1, 2], and more than 80% of patients with lung cancer are diagnosed with NSCLC [3]. In Yunnan Province, especially in Xuanwei City, the incidence of NSCLC is ranked the highest in world [4]. A mortality survey, conducted in Xuanwei during 2004–2005, revealed that the average mortality rate in NSCLC was 91.3 per 100, 000 in Xuanwei City [5], and the age at onset of lung cancer was 15–25 years younger than the average age in China [5]. Another survey, conducted in Xuanwei during 2011–2013, indicated that the male and female mortality rates in NSCLC were three and six times higher, respectively, in this region than in rural areas of China [6]. Due to these factors, NSCLC occurring in this area is known as Xuanwei Lung Cancer (XWLC). Even in the past few decades, despite the latest developments in NSCLC diagnosis and therapeutics, the overall survival (OS) rate of NSCLC patients still remains poor [7]. Therefore, it is essential to identify novel biomarkers with high specificity and sensitivity to increase the effectiveness of diagnostic, prognostic, and therapeutic strategies for NSCLC.
O-Glycan synthesis enzyme glucosaminyl (N-acetyl) transferase 3, mucin type (GCNT3/GNTM/C24GNT/C2GNT2/C2GNTM/C2/4GnT) is on human chromosome 15q22.2 and encodes a member of the N-acetylglucosaminyltransferase family, which catalyzes the formation of core 2 and core 4 O-glycans on O-linked glycosylation in mucin biosynthesis [8]. Cancer cells develop abnormal accumulation of mucin that covers epithelial cells and hinders drug access, which leads to cancer chemoresistance [9]. Recent research has shown that GCNT3 promotes human tumorigenesis. Rao et al. (2016) reported that GCNT3 is significantly overexpressed in human pancreatic cancer, and high GCNT3 expression is associated with the smoking and drinking habits of patients, as well as a diagnosis of diabetes [10]. GCNT3 has also been found to be under-expressed in stage II colon tumors compared with human normal colon tissues [11]. Moreover, GCNT3 gene is probably involved in the immunotoxicity, teratogenicity, and disruption of lipid metabolism caused by polycyclic aromatic hydrocarbons [12]. However, little is known about the role of GCNT3 in NSCLC, especially in XWLC, and its biological function and clinical significance has not yet been established.
MicroRNA (miRNAs) regulate a broad range of cellular processes, including cell growth, proliferation and apoptosis, through post-translational modulation of their targets [13]. We focused on miR-302b-3p in this study, because it has been confirmed – through a new method known as miRNA immunoprecipitation (miRIP) assay, established by Xiaoping Su – that it targets GCNT3 [14]. MiR-302b is a promising biomarker and candidate therapeutic target for multiple types of human cancers [15-17]. It has been found to sensitize breast cell lines in cisplatin treatment [15], suppress cell proliferation in hepatocellular carcinoma [16], and induce apoptosis in epithelial ovarian carcinoma [17]. Therefore, it is of great importance to clarify the regulatory pattern between GCNT3 and miR-302b-3p during NSCLC tumorigenesis.
In this study, GCNT3 was found via mRNA microarray to be up-regulated in NSCLC tissues compared with adjacent noncancerous tissues. qRT-PCR and western blot analysis demonstrated that GCNT3 was overexpressed in both NSCLC cancer cell lines (XWLC and NCI-H292) and primary lung cancer tissues from patients in Xuanwei, Yunnan Province, China. Knockdown of GCNT3 inhibited NSCLC cell proliferation, migration and invasion, while GCNT3 overexpression promoted these activities. miRIP and dual_luciferase reporter gene assays showed that miR-302b-3p could target GCNT3. The modulatory pattern between GCNT3 and miR-302b-3p and its role in XWLC tumorigenesis was investigated in NSCLC cells. GCNT3 was found to be a target of miR-302b-3p and enabled miR-302b-3p to participate in the pathogenesis and progress of NSCLC. Moreover, western blot experiments revealed that E-cadherin, N-cadherin, vimentin, cyclin D1 and p-Erk were downstream molecules of the miR-302b-3p/GCNT3 pathway. Our results indicated that GCNT3, negatively controlled by miR-302b-3p, regulated proliferation, migration and invasion in NSCLC cells by facilitating the activation of the ERK signaling pathway and EMT, which could be a potential biomarker for NSCLC diagnosis and therapy.
Materials and Methods
Collection of specimens and clinical data
This study was approved by the Ethics Committee of School of Medicine Yunnan University, and performed in accordance with the standards set by the Declaration of Helsinki. Written informed consent was obtained from each subject before participation in this study. All samples of NSCLC tissues and adjacent non-tumor tissues were collected from patients undergoing surgical procedures at the Third Affiliated Hospital of Yunnan Kunming Medical University. Only samples from NSCLC patients who had not undergone operative chemotherapy, radiotherapy, or targeted therapy were collected. Among them, 15 pairs of NSCLC and adjacent non-cancerous lung tissue samples obtained during surgery were snap-frozen immediately in liquid nitrogen, and stored at –80°C until further processing to avoid degradation of the RNA. Another 92 pairs of NSCLC and adjacent noncancerous lung tissues were paraffin-embedded. The clinicopathological features of the patients were collected from medical records and were summarized in Table 1. Tumor staging was based on the seventh edition of the Union for International Cancer Control (UICC) TNM staging system for lung cancer [18]. All 92 NSCLC cases underwent follow-up via telephone, and the OS rate of each patient was calculated from the time of the initial surgery to the end of followup or the day of death.
Table 1.
Correlation between GCNT3 expression and clinicopathologic parameters in NSCLC Patients (n=92). Note: TNM, tumor-node-metastasis staging system; P values were calculated using chi-square test; *p < 0.05, **p < 0.01
Cell lines
Two human NSCLC cell lines (A549 and NCI-H292) and a human normal pulmonary epithelial cell line (BEAS-2B) were purchased from American Type Culture Collection (Manassas, VA), and another two local NSCLC cell lines (XWLC and GLC) were isolated from the cancer tissues of patients lived in Xuanwei, a region in the Yunnan Province of China. Among them, A549, GLC and XWLC cell lines were cultured in Dulbecco’s modified Eagle’s medium (DMEM, HyClone, Logan, UT) medium containing 10% fetal bovine serum (HyClone) and 1% penicillin and streptomycin (Life Technologies, Carlsbad, CA), while NCI-H292 and BEAS-2B cell lines were maintained in RPMI 1640 (HyClone) medium with 10% fetal bovine serum (HyClone) and 1% penicillin and streptomycin (Life Technologies), at 37°C in a humidified air atmosphere containing 5% CO2.
Human Gene Expression Microarray assay
Total RNA from 15 pairs of NSCLC and adjacent non-cancerous lung tissues were isolated using the TRIzol agent (Invitrogen, Carlsbad, CA) following the manufacturers’ instruction. The quality of RNA samples was tested via the NanoDrop2000 (Thermo-Fisher Scientific, Waltham, MA, USA). Analysis of the expression profile of NSCLC-related genes was performed by a commercially available oligonucleotide microarray (Affymetrix GeneChip Human Transcriptome Array 2.0, Thermo Fisher Scientific) following the manufacturer’s recommendations. The array images were further analyzed by Agilent Feature Extraction software (version_10.7.1.1, Agilent Technologies, Santa Clara, CA). The raw data were normalized using the quantile algorithm. Differentially expressed genes were then identified through fold changes, and the P values were calculated with the Student’s t-test. The threshold values set for up- and down-regulated genes was a fold change ≥ 2.0 and a P value< 0.05.
Immunohistochemistry (IHC)
GCNT3 protein expression was assessed by immunohistochemistry (IHC) in 92 pairs of NSCLC and adjacent non-cancerous lung tissue sections. Tissues were fixed in 10% neutral formalin and embedded in paraffin. The paraffin-embedded NSCLC specimens were sliced into 4–μm–thick sections. Sections were then deparaffinized with xylene and rehydrated with a graded series of ethyl alcohol. Antigen retrieval was achieved by microwave treatment in 0.01 M sodium citrate buffer (pH 6.0) for 10 min. Endogenous peroxidase was inactivated by treatment with 0.3% hydrogen peroxide for 10 min. Tissue sections were heated for 5 min at 121°C and then cooled to 85°C. To prevent non-specific binding, we blocked the sections with goat serum (BioGenex, San Ramon, CA) for 30 min at room temperature. We then incubated the sections with anti-GCNT3 (TA332087, 1: 50, Origene, Rockville, MD) at 4°C overnight, followed by biotinylated goat anti rabbit IgG secondary antibody. For the negative control (NC), sections were treated with phosphate-buffered saline (PBS) instead of antibody. After washing, immunoreactivity was developed using diaminobenzidine tetrahydrochloride (DAB) solution (Vector Laboratories, Burlingame, CA). Tumor specimens were then counterstained with hematoxylin, dehydrated in alcohol, and mounted. The scores were calculated by evaluating the staining intensity and percentage of stained cells in representative areas via two investigators blinded to the clinical parameters. The expression score in IHC, defined as the staining intensity multiplied by the percentage of tumor positive area, was semi-quantitatively determined by the H-score (0–300). A median value derived from the testing set and applied to the validation set determined the cut-off point for tumoral GCNT3 high/low-level expression. The cut-off point was set at 75.
RNA extraction and quantitative real-time PCR (qRT-PCR) assays
Total RNA was extracted from NSCLC specimens or cultured cells by TRIZOL reagent (Invitrogen, Carlsbad, CA) following the manufacturer’s instructions and treated with DNase I (Invitrogen, Carlsbad, CA) to eliminate potential DNA contamination. The extracted RNA was reverse-transcribed into cDNA by using a Reverse Transcription Kit (Takara, Dalian, China) according to the manufacturer’s protocol. Quantitative real-time PCR (qRT-PCR) was performed with SYBR Premix Ex Taq (Takara, Dalian) on a CFX96 TouchTM Real-Time PCR Detection System (Bio-Rad Laboratories, Hercules, CA), and the data were collected. The mRNA transcription levels were normalized to the expression of GAPDH, and the miRNA transcription levels were normalized to the expression of U6. The sequences of the primers were as follows: GCNT3 (forward: 5′-GCCAGTAAGCTGGTTCGG-3′, reverse: 5′-GCCAGTAAGCTGGTTCGG-3′), GAPDH (forward: 5′-GACAGTCAGCCGCATCTTCT-3′, reverse: 5′-TTAAAAGCAGCCCTGGTGAC-3′), miR-302b-3p (loop: 5′-GTCGTATCCAGTCCAGGGACCGAGGACTGGATACGACCTACTA-3′, forward: 5′-GCGTAAGTGCTTCCATGTT-3′, reverse: 5′-TCCAGGGACCGAGGA-3′) and U6(loop: 5′-CGCTTCACGAATTTGCGTGTCAT-3′, forward: 5′-GCTTCGGCAGCACATATACTAAAAT-3′, reverse 5′-CGCTTCACGAATTTGCGTGTCAT-3′.). The expression levels of the genes of interest, GCNT3 or miR-302b-3p, were relative to the threshold cycle (CT) values, and converted to fold changes. All experiments were conducted in triplicate.
Plasmid generation and transfection
The full-length GCNT3 sequence was synthesized by Beijing Genomics Institute (Beijing, China) and subcloned into the pEX-3 vector (GenePharma, Shanghai, China) to generate the pEX-3-GCNT3 following the manufacturer’s instructions. An empty pEX-3 vector was used as a control. The recombinant plasmids were extracted by Pure YieldTM Plasmid Midiprep kit (Promega Corp., Madison, WI) and sequenced (Sain Biological Corporation, Shanghai, China). pEX-3-GCNT3 and empty pEX-3 plasmids (5µg plasmid DNA) were transfected into NSCLC cells using Lipofectamine® 3000 Reagent (Thermo Fisher Scientific) in a 6-well cell culture plate according to the manufacturer’s instruction, respectively. Cells were harvested 48 hours after transfection, to detect overexpression or for further functional assays.
Small interfering RNA and miRNA transfection
For RNAi-mediated knockdown of GCNT3, small interfering RNAs (siRNAs) specific to GCNT3 (GCUACUGCGAGCUGUGUAUTT) and a control siRNA (UUCUCCGAACGUGUCACGUTT) were designed and synthesized by GenePharma (shanghai, China). The miR-302b-3p mimics (UAAGUGCUUCCAUGUUUUAGUAG) and control NC mimic (UUCUCCGAACGUGUCACGUTT) were also provided by GenePharma. XWLC and NCI-H292 cells were transfected with the RNA oligoribonucleotides with Lipofectamine RNAiMAX (Invitrogen, Carlsbad, CA) and Opti-MEM (Gibco, Carlsbad, CA) according to the manufacturer’s recommendations. Cells were then incubated for 48h before further functional assays.
Cell proliferation assays
Cell proliferation was measured using Cell Counting Kit-8 (CCK-8) Assay kit (Dojindo, Kusatsu, Japan). The transfected cells were respectively seeded into 96-well plates (Corning Costar, Corning, NY) at a density of 3.0×103/ well/100 µL, and incubated in a water-saturated carbon dioxide incubator. CCK8 solution (10μl) was added to each well at 37°C for 1h, and the absorbance in each well was measured at a wavelength of 450 nm with the Thermo Scientific Multiskan FC microplate reader. Measurements of cell proliferation were performed at 24, 48, and 72 hours after transfection. All experiments were performed in triplicate. For the colony assay, cell transfection was performed 48h later as described above. Cells were then placed into 6-cm cell culture dishes at a density of 3000 cells per dish, and incubated for 10 days with medium changed every 3 days. The culture dishes were washed with phosphate buffer saline (PBS), fixed with 4% paraformaldehyde, stained with 0.25% crystal violet and counted. Colony-forming rates were then analyzed.
Cell migration and invasion assays
Cell migration and invasion assays were performed using 24-well Transwell chambers with 8-μm pores filters (Corning Costar, NY). Thirty-six hours after transfection, cells were serum starved for 12 hours, digested with trypsin, and resuspended in serum free medium. For the migration assay, a total of 10×104 cells in 200μL serum-free medium were seeded into the upper chambers of Transwell plates. For the invasion assay, the same numbers of cells were placed into the upper chamber of a well coated with Matrigel (BD Biosciences, San Jose, CA). The lower chambers were filed with 800µl of medium with 20% fetal bovine serum as a chemoattractant. After incubation at 37°C with 5% CO2 for 48h, non-migratory cells were removed from the upper surface of the filter with cottons swabs, while migrated cells on the lower surface of the filter were fixed in 4% paraformaldehyde and stained with Giemsa. The numbers of migrating or invading cells were counted in 5 randomly selected fields under a microscope (Olympus, Tokyo). All experiments were performed in triplicate.
Flow cytometric analysis for cell cycle and apoptosis
Cells were harvested and washed with PBS 48 hours after transfection. For the cell cycle analysis, cells were fixed with ice-cold 70% ethanol at 4°C overnight. The fixed cells were washed with PBS, incubated with ribonuclease A (BD Biosciences) at 37°C for 30 min, and stained with propidium iodide (PI, BD Biosciences) for 30 min. DNA content was analyzed on a flow cytometer (BD FACSJazz flow cytometer, BD Biosciences). The percentages of the cells in G0–G1, S, and G2–M phases were assessed with Flowjo7.61 (Tree Star Inc., Ashland, OR), and the difference between each group was compared. For the apoptosis assays, the transfected cells were double stained with FITC-Annexin V and PI by using FITC Annexin V Apoptosis Detection Kit (BD Biosciences) according to the protocols. The stained cells were analyzed by BD FACSJazz flow cytometer (BD Biosciences), and the data was analyzed using Cell Quest software (BD Biosciences). The percentages of apoptotic cells were counted and compared between each group. All measurements were repeated in triplicate.
miRIP assay
The following experiments were performed as described previously [14]. The probes are listed in the Table 3. Cells were transfected with the biotin-tagged specific probe or control probe and harvested 24 h after transfection. Cells were crosslinked by formaldehyde, equilibrated in glycine buffer and scraped with lysis buffer. Cell samples were sonicated and centrifuged. The supernatant was transferred to a 2 mL tube and 50 μl was separately saved for input analysis. The supernatant lysates were incubated with M-280 beads for 1h with rotating. The beads-sample mixture was washed and incubated to reverse the formaldehyde cross-links. Subsequently, RNAs were re-purified by TRIzol. Multiplex miRNA array was used to analyze the sequences of 184 miRNAs simultaneously from a single RT-PCR reaction to discover the miRNAs that would target GCNT3 in XWLC cells.
Western blot
Total protein was extracted using lysis buffer (Biomiga, Inc, San Diego) supplemented with protease inhibitor cocktail (Roche, Mannheim, Germany) and phosphatase inhibitor cocktail (Roche), quantified using the bicinchoninic acid (BCA) method, and denatured with called sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) sample buffer (Takara). Protein samples (15 µg) were separated by 10% SDS-PAGE, and transferred onto polyvinylidene difluoride (PVDF, 0.20µm) membranes (Millipore, Bedford, MA). The membranes were blocked with blocking buffer (5% skim milk or bovine serum albumin in Tris-buffered saline with Tween 20) and then incubated with appropriate primary antibodies namely, GCNT3 (ab77726, 1: 1000, Abcam, Cambridge, UK), E-cadherin (ab15148, 1: 500, Abcam, Cambridge, UK), N-cadherin (ab18203, 1: 1000, Abcam, Cambridge, UK), vimentin (ab137321, 1: 3000, Abcam, Cambridge, UK), Cyclin D1 (ab134175, 1: 1000, Abcam, Cambridge, UK), Erk1/2 (ab17942, 1: 1000, Abcam, Cambridge, UK), p-Erk (ab76299, 1: 5000, Abcam, Cambridge, UK), β-Tubulin (MA5-16308, 1: 5000, Invitrogen) overnight at 4°C. The membranes were washed for 10min three times with TBST and subsequently incubated with the corresponding secondary antibody (1: 5000, Santa Cruz Biotechnology, Santa Cruz, CA) for 1 h at room temperature. Membranes were then washed another three times for 10 min each with TBST. Protein bands were visualized using Super Signal electrochemiluminescence (#34580, Thermo Scientific) and quantitated with image Pro Plus 6.0, and the data were normalized to β-tubulin loading controls. All analyses were performed at least three times.
Dual-luciferase reporter assay
The oligonucleotides in the 3’ untranslated region (UTR) of the human GCNT3 harboring a putative miR-302b-3p target binding sequences was synthesized by GenePharm (Shanghai, China). As an NC, oligonucleotides containing a mutated miR-302b-3p targeting site, which was created by replacing target binding sequences with complementary bases, were also designed. The restriction enzyme cleavage sites of SacI and XhoI in the multiple cloning site (MCS) of the GP-miRGLO dual-luciferase miRNA target expression vector (Promega) were chosen as cloning sites. The wild-type and mutant 3’UTR oligonucleotides were subcloned into the downstream of the firefly luciferase gene using GP-miRGLO vector (Promega) to generate the GP-miRGLO-GCNT3-wild-type and GP-miRGLO-GCNT3-mutant, following the manufacturer’s instructions. The recombinant plasmids were extracted by Pure YieldTM Plasmid Midiprep kit (Promega) and sequenced (Sain Biological Corporation, Shanghai).
For the luciferase reporter assay, 293T cells were seeded in a 24-well plate the day before transfection. The cells were co-transfected with 2µg of the recombinant plasmids (wild-type or mutant GCNT3-3’UTR) and 50nM miRNA mimics (miR-302b-3p or NC) per well. At 24 and 48 h after co-transfection, cell lysates were collected, and luciferase activities were detected using a Dual-Luciferase Reporter System (Promega).
Statistical analysis
All statistical analyses were performed using GraphPad Prism 5.0 (GraphPad software, La Jolla, CA) and SPSS version 19.0 (SPSS, Inc., Chicago, IL). Continuous variables between two groups were estimated by the Student’s t-test. The relationship between GCNT3 protein expression in the IHC test and the clinicopathological factors (for categorical data) was calculated using the Chi-square test. Kaplan-Meier survival analyses were performed in 92 cases of specimens with lower or higher GCNT3 expression, and the difference was analyzed by the log-rank test. For each test, a two-sided P value < 0.05 was considered to indicate a statistically significant difference. *, **, and *** indicate statistical significance at P values of 0.05, 0.01 and 0.001, respectively.
Results
Differentially expressed mRNAs in NSCLC tissues
We performed a human gene expression microarray on 15 pairs of NSCLC and adjacent non-cancerous lung tissues. Differential expression of genes was identified by fold changes as well as P values calculated with the t-test. Since the threshold had been set for up-regulated and down-regulated genes at a fold change ≥ 2.0 and a P value <0.05, a total of 623 up-regulated and 2118 down-regulated mRNAs were identified in NSCLC samples compared with adjacent normal lung tissues by Student’s t-test method, as depicted in the volcano plot (Fig. 1A) and the heat map (Fig. 1B). Among them, SPP1 [19], GCNT3 [11], TMPRSS4[20], SPINK1 [21], HABP2 [22], and CXCL13 [23] are known oncogenes for tumor i-genesis, while two, namely TMPRSS4 and HABP2, were identified as oncogenes in lung cancer. The top 10 up-regulated mRNAs in NSCLC tissues are listed in Table 2. We selected GCNT3 for the subsequent analysis.
Fig. 1.
GCNT3 was overexpressed in NSCLC, and its high expression was involved with poor prognosis. (A) Volcano plot of P values as a function of the weighted fold change for mRNAs in NSCLC tissues. (B) Heatmap of expression profiles for the 2741 mRNAs showed significant expression changes (623 up-regulated and 2118 downregulated). (C) mRNA levels of GCNT3 were examined in 15 pairs NSCLC tissues and noncancerous tissues by qPCR, and normalized to GAPDH expression (P< 0.001). (D) GCNT3 expression in NSCLC and normal tissues based on TCGA dataset (P< 0.001). (E) The protein expression in NSCLC tissues and adjacent noncancerous tissues was assessed by immunohistochemistry assay. (F) GCNT3 was upregulated in Carcinoma tissues compared to the Par-carcinoma tissues, determined by IHC assay (P< 0.001). (G) Kaplan-Meier survival analysis showed that the overall survival of NSCLC patients with higher GCNT3 expression was significantly shorter than those with lower GCNT3 expression (P< 0.05; Scale bar: 100 μm).
GCNT3 expression was upregulated in NSCLC specimens, and high expression of GCNT3 correlated with poor prognosis in NSCLC patients.
To investigate the role of GCNT3 in the tumorigenesis of NSCLC, we firstly examined the expression of GCNT3 in NSCLC tissues. In Fig. 1C, 15 pairs of NSCLC and adjacent non-cancerous tissues had markedly increased levels of GCNT3 mRNA in NSCLC tissues compared with adjacent non-cancerous tissues by real-time PCR (P < 0.001). Similar results were observed after re-analyzing gene expression data in 226 NSCLC tissues and 20 normal tissues downloaded from The Cancer Genome Atlas (TCGA: https://tcga-data.nci.nih.gov/docs/publications/tcga/), a public dataset with clinical information and mRNA expression levels of all patients at the Oncomine website (https://www.oncomine.org) (Fig. 1D). Taken together, these results indicated that GCNT3 might play a critical role in the progression of NSCLC.
We then evaluated GCNT3 expression in 92 pairs of NSCLC and adjacent non-cancerous tissues using IHC staining. The results indicated that GCNT3 was highly increased in NSCLC tissues at the protein level (Fig. 1E-1F). Furthermore, according to IHC, we divided the 92 NSCLC patients into a high GCNT3 expression group (above average GCNT3 expression, n = 55) and a low expression group (below average GCNT3 expression, n = 37). As presented in Table 1, we explored the clinical significance of GCNT3 in NSCLC, and the correlation between GCNT3 expression levels and the patients’ clinical features were analyzed by the Chi-square test. The results implied that the expression of GCNT3 was closely associated with TNM stage (P = 0.0078), lymph node metastasis (P = 0.0016), and survival years after surgery (P = 0.0033). Nevertheless, there was no correlation between the expression level of GCNT3 and age, sex or tumor size (all P > 0.05). The correlation between the expression of GCNT3 protein and prognosis in NSCLC patients was assessed by Kaplan-Meier survival curves and a log-rank test. The results showed that higher expression of GCNT3 correlated with shorter OS among NSCLC patients (P = 0.042, Fig. 1G). The P value was assessed by the log-rank test.
Knockdown of GCNT3 inhibited NSCLC cell proliferation, migration, and invasion, and induced cell cycle arrest and apoptosis.
To further investigate the biological function of GCNT3 in NSCLC, we detected the expression of GCNT3 in four human NSCLC cell lines (A549, GLC, XWLC and NCI-H292) and a human Normal pulmonary epithelial cell line (BEAS-2B). As illustrated in Fig. 2A, XWLC and NCI-H292 cells exhibited higher expression of GCNT3 at the protein level. Thus, we knocked down endogenous GCNT3 expression by siRNA transfection in XWLC and NCI-H292 cells (Fig. 2B). After GCNT3 knockdown, growth curves detected by CCK8 assay implied that the proliferative ability of both XWLC and NCI-H292 cells was remarkably inhibited (all P < 0.05, Fig. 2C). Furthermore, the colony formation assay indicated that the colony formation rates of GCNT3-knockdown XWLC and NCI-H292 cells were much lower than those transfected with siNC (Fig. 2D). These data indicated that downregulation of GCNT3 repressed the proliferation of NSCLC cells.
Fig. 2.
GCNT3 protein levels in NSCLC cell lines, and GCNT3 knockdown inhibited NSCLC cell proliferation. (A) Western blot assay showing GCNT3 protein levels in four human NSCLC cell lines (A549, GLC, XWLC and NCI-H292) and a human Normal pulmonary epithelial cell line (BEAS-2B). (B) Western blot below suggesting significant knockdown of GCNT3 mediated by siRNA after 48hours post-transfection in two NSCLC cell lines (XWLC and NCI-H292). β-tubulin was used as the loading control. (C) CCK8 assay was performed to measure the proliferation of XWLC and NCI-H292 cell lines at 0, 24, 48 and 72 h after siRNA and siNC (control) transfection. (D) Colonyforming assay was used to detect the proliferation of si-GCNT3-transfected XWLC and NCI-H292 cell lines. Data was based on at least 3 independent experiments, and shown as mean ± SD. (*P< 0.05; **P< 0.01; ***P< 0.001).
A flow cytometry assay was conducted to identify the stages of cell cycle progression in XWLC and NCI-H292 cells transfected with GCNT3-siRNA or NC-si. The results indicated that GCNT3-siRNA transfection in XWLC and NCI-H292 cells led to a remarkable increase in the percentages of cells at G0/G1 phase and a marked decrease in cells at S phase. G2/M phase was not obviously affected after downregulation of GCNT3. (all P < 0.01, Fig. S1A - for all supplemental material see www.karger.com/10.1159/000494482/) To investigate the potential mechanisms underlying the suppression of growth after GCNT3 downregulation, annexin V-FITC/PI staining assay was performed to assess the effect on NSCLC cell apoptosis. The results revealed that the percentage of cells in early apoptosis in XWLC-siGCNT3 and NCI-H292-GCNT3 was higher than that in XWLC-NC (P < 0.01) and NCI-H292-NC (P < 0.01), and the percentage of cells in late apoptosis in XWLC-siGCNT3 and NCI-H292-GCNT3 was also higher than in XWLC-NC (P < 0.01) and NCI-H292-NC (P < 0.05). These results indicated that the inhibition of GCNT3 expression promoted apoptosis in NSCLC cells, which suggested that GCNT3 contributed to tumor progression by inhibiting apoptosis (Fig. S1B).
Since the metastasis of cancer cells is known to be a crucial indicator of tumor progression, we detected the migration and invasion capacities of NSCLC cells after transfecting them with siGCNT3 or siNC via transwell assays with non-Matrigel-coated and Matrigel coated inserts, respectively. As shown in Fig. S1D and 3E, GCNT3 knockdown clearly repressed the migration and invasion capacities of XWLC and NCI-H292 cells (all P < 0.01).
To investigate the regulatory mechanism of GCNT3, we detected certain vital proteins that might be involved in EMT, cell cycle progression and the Erk pathways. As shown in Fig. S1C, western blotting assay showed that the protein levels of cyclin D1, N-cadherin, vimentin and p-Erk (phosphorylation of Erk) were downregulated, while thoses of E-caderherin were markedly increased, after GCNT3 knockdown. However, the levels of Erk protein showed no significant change. These data suggested that GCNT3 could regulate migration, invasion, the cell cycle and Erk pathway associated proteins.
GCNT3 overexpression promoted proliferation, migration and invasion in NSCLC cells
To further unravel the function of GCNT3 in the oncogenesis of NSCLC, GCNT3 was significantly upregulated in XWLC and NCI-H292 cells by GCNT3 plasmid transfection. The colony formation assay results suggested that upregulation of GCNT3 expression enhanced colony forming abilities in XWLC and NCI-H292 cells (Fig. S2B-C, P < 0.05). Tumor migration (Fig. S2D) and invasion (Fig. S2E) were also enhanced by transfecting GCNT3 plasmid into XWLC and NCI-H292 cells. We detected a couple of proteins that might be involved with the Erk signaling pathways and EMT. Western blotting indicated that the protein levels of N-cadherin, vimentin and p-Erk were upregulated, while the level of E-cadherin was markedly decreased, after GCNT3 overexpression (Fig. S2A).
miR-302b-3p directly targeted the GCNT3, and regulated NSCLC tumor progression
miRNAs can regulate a broad range of cellular processes through post-translational modulation of their targets [13]. To understand the mechanism by which GCNT3 affected NSCLC progression risk after miRNA, the miRIP method was used to identify the miRNAs targeting GCNT3 [14]. Ten probes of GCNT3 were transfected into XWLC cells to hybridize with GCNT3 mRNA for 24 h (Fig. S3A). The associated mRNAs of GCNT3 were purified using Streptavidin Dynabeads (M-280) and isolated by TRIzol. Multiplex miRNA array was used to analyze the sequences of 184 miRNAs simultaneously from a single RT-PCR reaction to discover the associated miRNAs targeting GCNT3 in XWLC cells. miR-302b-3p, interacted with GCNT3 mRNA robustly, was selected for further investigation (Fig. S3B). We then detected the expression of miR-302b-3p in four human NSCLC cell lines (A549, GLC, XWLC and NCI-H292) and a human normal pulmonary epithelial cell line (BEAS-2B). As illustrated in Fig. S3C, XWLC and NCI-H292 cells exhibited lower expression of miR-302b-3p. These results indicated that miR-302b-3p correlated negatively with GCNT3 in NSCLC cell lines.
To further investigate the function of miR-302b-3p in NSCLC cells, we transfected a miR-302b-3p mimic into XWLC and NCI-H292 cells. As illustrated in Fig. S3D, miR-302b-3p overexpression inhibited GCNT3 expression at the protein level in XWLC and NCI-H292 cells. After transfection with miR-302b-3p, growth curves detected by the CCK-8 assay implied that the proliferative ability of XWLC and NCI-H292 cells was remarkably inhibited (all P < 0.05, Fig. S3E). Tumor migration (Fig. S3F) and invasion (Fig. S3G) activities were also suppressed by transfection with miR-302b-3p in XWLC and NCI-H292 cells. Our data suggested that miR-302b-3p directly targeted GCNT3 and was involved in tumorigenesis in NSCLC.
To confirm whether GCNT3 gene was a direct target of miR-302b-3p, a dual luciferase reporter assay was performed in 293T cells. We first performed a bioinformatics analysis using Miranda (http://34.236.212.39/microrna/home.do) to find out the putative miR-302b-3p target binding sequences in the 3’UTR of the human GCNT3. (Fig. 3A) The partial sequence in the 3’UTR of the human GCNT3 harboring putative miR-302b-3p target binding sequences was synthesized. As an NC, a fragment with the same part of GCNT3 3’UTR with the miR-302b-3p binding site replaced by complementary bases was also synthesized. These two fragments were subcloned into a dual-luciferase reportor vector to generate the GP-miRGLO-GCNT3-wild-type and GP-miRGLO-GCNT3-mutant. (Fig. 3B) The miRNA (miR-302b-3p or NC mimics, 50nM) and the recombinant plasmids (GP-miRGLO-GCNT3-wild-type or GP-miRGLO-GCNT3-mutant, 2µg) were co-transfected into 293T cells for 24 h and 48 h. The results indicated that miR-302b-3p strongly inhibited the activity of the firefly luciferase gene of the reporter vector containing the GCNT3 3’UTR with a wild-type miR-302b-3p binding site, but did not inhibit the reporter vector with mutant GCNT3 3’UTR (Fig. 3C). This result indicated that miR-302b-3p inhibited the expression of GCNT3 by directly targeting the 3’UTR.
Fig. 3.
miR-302b-3p was found to be the GCNT3-binding miRNA through dual luciferase reporter gene assays. (A) The putative miR-302b-3p target binding sequences in the 3’UTR of the human GCNT3 predicted by bioinformatics analysis was listed below. (B) The GCNT3-wild-type with putative binding site of miR-302b-3p and the mutant generated by replacing the binding region of miR-302b-3p with complementary bases were subcloned into the GP-miRGLO vector, respectively. The restriction enzyme cutting sites of the SacI and XhoI were chosen as the cloning sites. (C) For luciferase assay, the cells were co-transfected with 2ug of the recombinant plasmids (WT or mutant GCNT3-3’UTR) and 50nM miRNA mimics (miR-302b-3p or NC) per well. Luciferase activities were measured twenty-four and forty-eight hours after transfection. Firefly luciferase activities were normalized by the activity of renilla luciferase. The experiments were performed in triplicate. (*P< 0.05; **P< 0.01; ***P< 0.001).
miR-302b-3p overexpression suppressed NSCLC cell proliferation, migration and invasion in a GCNT3-dependent manner
We found that GCNT3 was the target of miR-302b-3p, and facilitated NSCLC progression. Thus, rescue experiment was conducted to identify whether GCNT3 overexpression could reverse the phenotype of miR-302b-3p overexpression. The colony formation assay results revealed that the overexpression of miR-302b-3p suppressed colony formting abilities in XWLC and NCI-H292 cells (Fig. S4B, P < 0.05), while miR-302b-3p and GCNT3 co-overexpression reversed the miR-302b-3p overexpression-induced decline in colony forming abilities (Fig. S4B). The transwell assay with or without Matrigel suggested that miR-302b-3p overexpression markedly inhibited migration and invasion in NSCLC cells, while miR-302b-3p and GCNT3 co-overexpression restored these abilities compared with miR-302b-3p overexpression. (Fig. S4C and D, P < 0.01) In addition, we examined the expression of some proteins regulated by the Erk signaling pathway and EMT. Western blot results revealed that the overexpression of miR-302b-3p downregulated the protein levels of GCNT3, N-cadherin, vimentin, and p-Erk and elevated the levels of E-cadherin. In contrast, co-transfection with miR-302b-3p and GCNT3 abrogated the effects of miR-302b-3p on the levels of these proteins (Fig. S4A).
Discussion
Epidemiological studies [4-6] indicate that XWLC has peculiar traits, such as high rates of incidence and death, low male-to-female ratio, younger age at onset, and adenocarcinoma as the dominant type. Strong evidence revealed that XWLC is a special type of adenocarcinoma caused by indoor air pollution from household burning of solid fuels, which contain carcinogenics polycyclic aromatic hydrocarbons [24, 25]. Thereby, it is essential to clarify the pathogenetic mechanisms in this kind of lung cancer in the Xuanwei area to enable early detection, and prognostic and therapeutic strategies. In this study, differential expression of mRNAs was screened through mRNA microarray between NSCLC tissues and matched adjacent non-cancerous lung tissues from the patients in Xuanwei, Yunnan Province, China. We focused on GCNT3, because when it is up-regulated, it has been found to be associated with clinical indexes in patients with NSCLC. Knockdown of GCNT3 inhibited the proliferation, migration and invasion ability of NSCLC cells, while GCNT3 overexpression promoted these activities. Furthermore, GCNT3 was a target of miR-302b-3p as confirmed in miRIP and dual luciferase reporter gene assays. Overexpression of miR-302b-3p inhibited cell proliferation, migration and invasion by suppressing GCNT3 expression, while miR-302b-3p and GCNT3 co-overexpression reversed the miR-302b-3p overexpression phenotype. In addition, western blotting revealed that E-cadherin, N-cadherin, vimentin, cyclin D1 and p-Erk were downstream molecules of the miR-302b-3p/GCNT3 pathway.
GCNT3 is a member of the N-acetylglucosaminyltransferase family involved in mucin biosynthesis [8]. Abnormal expression of GCNT3 promotes the formation of mucins to enhance the flexibility of cancer cells and allow them to adapt to adverse tumor microenvironment conditions [9-12]. Analysis of clinical data from NSCLC patients and expression data from TCGA indicated that GCNT3 was up-regulated in NSCLC tissues and associated with TNM stage, lymph node metastasis and patients’ OS (Table 1 and Fig. 1). In addition, GNCT3 was also over-expressed in NSCLC cells compared with normal cells (Fig. 2). The GCNT3 inhibition and overexpression functional assays were performed in XWLC/NCI-H292 cell lines to explore the biological roles of GCNT3 in NSCLC tumorigenesis (Fig. 2, Fig. S1 and S2). The results showed that knockdown of GCNT3 exerted tumor-suppressive effects by reducing NSCLC cell proliferation, while GCNT3 overexpression promoted tumorigenesis by facilitating NSCLC cell proliferation. Transwell assay showed that knockdown of GCNT3 in XWLC and NCI-H292 cells suppressed cell migration and invasion, while overexpression of GCNT3 in A549 cells promoted migration and invasion. These results indicated that GCNT3 could be highly related to tumorigenesis in NSCLC.
As important regulators in downstream oncogenic pathways, miRNAs have been found to contribute to oncogenesis [26]. Alterations in miRNA expression have been reported in several human tumors including those in the colon [27], prostate [28], lung [29], breast [30], liver [31], pancreas [32], and bladder [33]. The mi-RIP method was performed to precisely identify the miRNAs that can bind to the GCNT3 mRNA in NSCLC cells [14]. miR-302b-3p was found to interact robustly with GCNT3 mRNA, and was selected for the next investigation. In recent studies, miR-302b has been found to be a tumor suppressor in various tumors, including esophageal cancer [34], hepatocellular carcinoma [35], and gastric cancer [36]. Zhang et al. (2017) reported that miR-302b decreased esophageal cancer growth by targeting ERBB4, IRF2 and CXCR4 and inhibiting CRI critical pathway and expression of downstream cytokines [34]. Wang et al. (2016) identified that miR-302b inhibited human hepatocellular carcinoma cell metastasis and invasion by targeting the AKT2/NF-kB/MMP-2 pathway [35]. In addition, miR-302b was also found to target and negatively regulate the expression of CDK2 to decrease cell proliferation and arrest the cell cycle in gastric cancer tissues by suppressing the Erk signal pathway [36]. In the present study, miR-302b-3p was found to be the GCNT3-binding miRNA through miRIP and dual luciferase reporter gene assays (Fig. S3 and Fig. 3). It was significantly down-regulated in NSCLC cell lines. The results also showed that overexpression of miR-320b-3p down-regulated the expression of GCNT3 and decreased the proliferation, migration and invasion abilities of NSCLC cells. (Fig. S3) The rescue assay showed that these effects could be abrogated by re-expression of GCNT3, indicating that miR-302b-3p acted as a repressor of NSCLC cells in a GCNT3-dependent manner. In addition, miR-302b-3p overexpression remarkably inhibited NCI-H292 cell proliferation after incubation for 24h, whereas the inhibition of cell proliferation by si-GCNT3 was revealed in NCI-H292 48 h after incubation. (Fig. S4) The reason for these results is that miR-302b not only inhibited the expression of GCNT3, but also modulated the expression of its target genes, such as ERBB4, IRF2, CXCR4 [34], AKT2/NF-kB/MMP-2 pathway [35], CDK2 [36], ErbB4 et al. [37], and AKT2 [38]. As we know, one single miRNA can regulate a broad range of cellular processes including cell growth, proliferation and apoptosis, by post-translational modulation of the target genes [13]. Therefore, the overexpression of miR-302b-3p might be more effective than down-regulation of GCNT3 in inhibiting the proliferation of NSCLC cells.
Recent evidence has shown that GCNT3 promotes human tumorigenesis [9, 10], and the EMT process has been confirmed to be an essential process for tumor invasion and metastasis [39]. E-cadherin [40], N-cadherin [41] and vimentin [42] have been identified as the main factors in the EMT process. In the present study, we revealed that knockdown of GCNT3 expression or restoration of miR-302b-3p can suppress the migratory and invasive phenotype in XWLC and NCI-H292 cells, and increase the level of E-cadherin, while decreasing the expression of N-cadherin and vimentin. These results suggested that GCNT3 and miR-302b-3p affect NSCLC metastasis partly via the EMT. In addition, the Erk signaling pathway has been shown to be frequently activated in numerous kinds of human cancers [43-45]. Aberrant activation of phosphorylated Erk1/2 induced by miRNAs has been observed in lung cancer [45-47]. Recent reports point out that miR-379 acts as a tumor suppressor in NSCLC by indirectly regulation of ERK signaling pathways [45], miR-7 inhibits the activation of ERK/MAPK signaling pathway by down-regulating FAK expression [46], and miR-622 overexpression decreases invasion in lung tumor cells by inhibiting HIF-1alpha via inactivation of ERK signaling in A549 cell [47]. Consistent with these findings, we found that the p-Erk level was significantly decreased in response to GCNT3 silencing or miR-302b-3p overexpression in XWLC and NCI-H292 cells, suggesting that the effect of the miR-302b-3p/GCNT3 signaling pathway in NSCLC may depend on the activation of p-Erk. Moreover, we investigated potential targets that could be responsible for GCNT3 mediated cell cycle arrest. In this study, the expression of cyclin D1 was significantly decreased in GCNT3-knockdown XWLC and NCI-H292 cells, which was consistent with the results of cell cycle analysis. Cyclin D1 is a well-established cell cycle regulatory protein, involved in the pathogenesis of several human tumor types, including lung cancer, breast cancer and pancreatic cancer [48-50]. Cyclin D1 controls cell cycle transition from the G1 to S phase, which is essential for the dividing cell to enter the DNA synthesis phase [51]. Thus, we speculated that GCNT3 could positively regulate cell cycle progression in NSCLC by altering the expression of cyclin D1. All together, these findings indicated that E-cadherin, N-cadherin, vimentin, cyclin D1 and Erk were the downstream factors of the axis of miR-302b-3p/GCNT3.
In conclusion, GCNT3 was regarded as an oncogene in NSCLC, and a target of miR-302b-3p, which regulated cell proliferation, migration and invasion in a GCNT3-dependent manner. Moreover, E-cadherin, N-cadherin, vimentin, cyclin D1 and p-Erk appeared to be downstream molecules of the miR-302b-3p/GCNT3 pathway.
Conclusion
Taken together, the miR-302b-3p/GCNT3 pathway may participate in the initiation and progression of NSCLC, and may be an important cancer promoter and novel candidate for NSCLC therapy.
Acknowledgements
We are grateful to all the patients participating in this study. This work was financially supported by the grants from the National Natural Science Foundation of China (No.81860499 to CJ. Xiao and No.81460435, No. 81860494 to SY. Zheng), and the Academic young scholar of distinction for doctoral post graduate in Yunnan province (C6155501).
Disclosure Statement
The authors declare that they have no conflict of interest.


Get Permission
