Using Oncomine database, GEPIA, and UALCAN to analyze the expression of BTG1
We obtained data from the BTG1 study in 453 different types of tumors through the Oncomine database. A total of 43 studies suggested statistical differences in BTG1 mRNA levels between tumors and normal tissues, of which 21 studies showed that BTG1 expression levels in tumors were significantly increased, and another 22 studies showed BTG1 expression levels in tumors significantly reduced. Compared with normal tissues, BTG1 expression was increased in brain and CNS cancer, cervical cancer, head and neck cancer, kidney cancer, and other cancers. BTG1 showed low expression in breast cancer, colorectal cancer, leukemia, lymphoma, ovarian cancer, and sarcoma. In addition, the expression of BTG1 in lung cancer was increased in one study and decreased in two studies (Figure 1A). Through GEPIA, we obtained the relative expression of BTG1 in each tumor and normal tissues (Figure 1B-C).
Next, we conducted a subgroup analysis of UCEC cases through UALCAN, which included 546 UCEC cases. The results showed that, compared with the normal control group, BTG1 mRNA expression was low in primary endometrial cancer tissues (P = 1.91589999420927E-09; Figure 2A). The expression level of BTG1 in stage 4 endometrial carcinoma was significantly lower than that in stage 1, stage 2, stage 3 (stage1 vs stage 4, P=3.029100E-04; stage2 vs stage 4, P =2.382500E-03; stage 3 vs stage 4, P=4.759200E-03; Figure 2B). For patients of different weights, there was no difference between normal weight, extreme weight, and obese; however, BTG1 expression of extreme obese group is higher than extreme weight or obese (extreme weight vs extreme obese, P=8.015500E -03; obese vs extreme obese, P=2.182600E-03; Figure 2C). For BTG1 expression in different age groups, the expression level of BTG1 in patients between 61 and 80yearsold was lower than the expression level of BTG1 in patients between 41 and 60 years old (age (41-60 yrs)-vs-Age(61-80 yrs), P=4.821900E-02; Figure 2D). For the expression level of BTG1 in different histological subtypes, the expression level of BTG1 in endometrioid was higher than in serous and mixed serous and endometrioid subtypes, but there was no difference in expression between serous and mixed serous or endometrioid. (endometrioid vs serous, P=1.76350000047343E-07; endometrioid vs mixed serous and endometrioid, P=1.958680E-03; Figure E) Furthermore, the expression level of BTG1 in TP53-Mutant EC was significantly lower than that of TP53-NonMutant EC (TP53-Mutant-vs-TP53-NonMutant; P=1.71269998183732E-09; Figure F). There was no difference in the expression of BTG1 between patients of different races or menopause status (Figure 2G-2H).
High expression of BTG1 can improve EC prognosis
In the KM plotter online analysis tool, we set the follow-up time to 120 months. A total of 543 EC cases in the database met these criteria. The results showed that patients with high BTG1 EC had significantly higher OS than patients with low BTG1 EC (HR = 0.52, 0.32–0.85, logrank P = 0.008) (Figure 2I). At the same time, in UALCAN, patients with high BTG1 expression (136 patients) had a longer survival time (P = 0.031) compared with patients with low/moderate BTG1 expression (407 patients) (Figure 2J). However, survival analysis of patients with different menopause status, body weight, and race indicated that there was no difference in survival time (Figure 2K-M). These results show that patients of EC with high BTG1 expression have a better prognosis.
BTG1 mRNA expression in EC and relative clinicopathological analysis
qRT-PCR detection showed that compared with normal human endometrial tissue, BTG1 had low expression in human EC tissue (P≤0.001, Figure 3A). Clinicopathological analysis showed the expression of BTG1 was related to invasion depth (P=0.031) and FIGO stage (P=0.012). The relative invasion depth of BTG1 in EC with invasion depth ≥1/2 showed less expression than EC with an invasion depth<1/2 EC. BTG1 was relatively under-expressed in FIGO stage II, III and IV compared to FIGO stage I (Figure 3B-C). BTG1 and EC pathological tissue type (P=0.450), fertility history (P=0.694), lymphatic metastasis (P=0.148), menopause (P=0.206), estrogen receptor (ER) (P=0.930), progesterone receptor (PR) (P=0.163), and age of diagnosis (P=0.227) were not related (Figure 3D-J).
BTG1 expression affects the prognosis of EC
A total of 70 patients with EC were followed up until February 15, 2019. The longest and shortest survival times were 60 months and 1 month, respectively. A total of 13 of 70 patients with EC died. Univariate results indicated that the average survival time of the BTG1 high-expression group was 57.2 months, and the average survival time of the BTG1 low-expression group was 48.9 months. Low expression of BTG1 was significantly correlated with OS shortening (P = 0.041) (Figure 3K). These results are consistent with the results we obtained based on TCGA.
Functional enrichment of BTG1 in EC patients
Using the functional module in LinkedOmics, the mRNA-related genes of 176 UCEC patients in TCGA were analyzed. As shown in the volcano graph, there were 2805 genes that significantly positively correlated with BTG1 (dark red dots), and 1267 genes that significantly negatively correlated with BTG1 (dark green dots) (FDR <0.01) (Figure 4A). The heat map showed that the top 50 gene sets were positively and negatively related to BTG1 (Figure 4B–C). We queried the functions of the top 50 gene sets positively and negatively related to BTG1. The results showed that BTG1 plays an important role in regulating embryonic development, tumorigenesis, apoptosis, and cell cycle. The statistical scatter plot of each gene showed the expression of BTG1 and SH3BGRL (Pearson-Correlation: 0.6438, P=9.797e-14), PCMTD1 (Pearson-Correlation: 0.6322, P=3.642e-13), and TET2 (Pearson-Correlation: 0.6207, P=1.273e-12) had a strong correlation (Figure 4D-F).
Next, we used Metascape for GO enrichment analysis to analyze the function of BTG1 and the top 200 related differentially expressed genes. The results showed that BTG1 and its related differentially expressed genes are mainly involved in the ING2 complex, protein-containing complex disassembly, translation, GTPase regulator activity, exchange factor activity, and cellular response to epidermal growth factor stimulus (Figure 5A-B, Table S3). The biological processes of BTG1 and its related genes included protein-containing complex disassembly, translation, cellular response to epidermal growth factor stimulus, and regulation of GTPase activity (Figure 5C-D, Table S4). The molecular functions of BTG1 and its related genes included GTPase regulator activity and guanyl-nucleotide exchange factor activity (Figure 5E-F, Table S5). In order to further analyze the relationship between BTG1 and UCEC, we conducted a PPI network analysis and analysis of important genetic components in MCODE. Research suggests that its biological function includes the following aspects: major pathway of rRNA processing in the nucleolus and cytosol,rRNA processing in the nucleus and cytosol, metabolism of RNA, mRNA splicing, major pathway, mRNA splicing, processing of capped intron-containing pre-mRNA, 55S ribosome, mitochondrial,mitochondrial translation elongation, mitochondrial translation termination, resolution of sister chromatid cohesion, separation of sister chromatids, mitotic prometaphase, ING2 complex, HDACs deacetylate histones, and chromatin organization (Figure 6A-B).
BTG1 miRNA target and interacting protein network in EC
In order to further study the functional targets of BTG1 in EC, we used miRDB to obtain 58 miRNAs that can bind to BTG1 (Table 1). Using the protein interaction network constructed by GeneMANIA (Figure 6C), we obtained interacting proteins, including HOXB9, a transcription factor involved in cell proliferation and differentiation. We further revealed the mutual regulation of BTG1 and HOXB9 (Figure 6D).
BTG1 inhibits cell proliferation, migration, invasion, and promotes apoptosis in ECCs
In order to study the effect of BTG1 on the malignant biological behavior of ECCs, we first transfected Ishikawa cells and HEC-1A cells with a BTG1 overexpression plasmid, a BTG1 knockdown plasmid, or a corresponding negative control. The transfection efficiency was verified by qRT-PCR. A knockdown plasmid with the best knockdown effect (BTG1-RNAi (6008-2) was selected for subsequent experiments (Figure S1). Through EDU experiments, we found that overexpression of BTG1 could inhibit the proliferation of Ishikawa cells and HEC-1A cells, while knockdown of BTG1 could promote the proliferation of Ishikawa cells and HEC-1A cells (Figure 7A). We obtained through a wound healing assay that overexpression of BTG1 inhibited the invasion of Ishikawa cells and HEC-1A cells. In contrast, knocking down BTG1 promoted the migration of Ishikawa cells and HEC-1A cells (Figure 7B). The transwell cell invasion assay showed that BTG1 overexpression inhibited the invasion of Ishikawa cells and HEC-1A cells. In contrast, knocking down BTG1 expression promoted the invasion of Ishikawa cells and HEC-1A cells (Figure 7C). Finally, we used flow cytometry to examine the effect of BTG1 knockdown and overexpression on apoptosis in Ishikawa cells and HEC-1A cells. The results showed that BTG1 overexpression promoted apoptosis and BTG1 knockdown suppressed apoptosis (Figure 7D).
BTG1 inhibits the EMT process in endometrial cancer cells
In order to explore the mechanism of BTG1 on the EMT process of ECCs, we detected the expression of E-cadherin, vimentin, and N-cadherin by western blot. Overexpression of BTG1 increased the expression of E-cadherin, while the expression of N-cadherin and vimentin decreased. Knocking down BTG1 reduced the expression of E-cadherin and increased the expression of N-cadherin and Vimentin (Figure 8A). Our results indicate that BTG1 can inhibit the EMT process in ECCs.