1.
Introduction
Esophageal cancer is the eighth most common type of cancer and the sixth most deadly in terms of mortality rates [1,2]. Worldwide, there are about 456,000 newly diagnosed cases of esophageal cancer every year [3]. Among several histological types of esophageal cancer, esophageal squamous cell carcinoma (ESCC) accounts for about 90% of all cases [3]. Currently, the main clinical treatment methods for ESCC include endoscopic resection, surgery, radiotherapy, and chemotherapy [4,5]. Early diagnosis and treatment of ESCC are beneficial for patients in a 5-year survival rate [6]. However, although it is easy to perform endoscopic examination in the esophagus, the early symptoms of ESCC are not obvious, causing a delay in diagnosis [7,8]. Quite a few ESCC patients were diagnosed at advanced clinical stages [9]; and worse still, ESCC patients with advanced clinical stages tend to have dropped significantly 5-year survival rate [6]. Therefore, to improve the prognosis of patients, it is critical to study the molecular mechanisms of ESCC. According to previous reports, several risk factors, such as drinking, smoking, diet, nutrition, heredity and body mass index, play essential roles in the occurrence and development of ESCC [1,10,11]. However, the molecular pathogenesis of ESCC is extremely complex, involving several prospects, such as host genetics, host-microbe interactions, and epigenetic imbalance [12,13,14,15]. Therefore, it is necessary to explore the molecular mechanisms of ESCC and find new clinical markers beneficial for identifying and predicting prognosis of the disease.
Belonging to the inhibitor of the apoptosis (IAP) family, baculoviral IAP repeat containing 5 (BIRC5, also known as API4, EPR-1, survivin) protein is an inhibitor of apoptosis encoded by its homonymous gene – BIRC5. In mammalian cells, it is involved in the control of mitosis, the regulation of apoptosis, and the cell stress response. Increased BIRC5 expression can be seen in multiple types of cancer, such as non-small cell lung cancer [15] and breast cancer [16]. In ESCC, BIRC5 protein level is also upregulated, which is associated with poor clinical prognosis of ESCC [17,18]. Meanwhile, BIRC5 can reduce the sensitivity of tumor cells to radiation induction [19]. These studies indicate that BIRC5 plays an important role in ESCC.
However, the current reports on BIRC5 and ESCC are limited; most of them only focus on the protein level and are based on a small sample size. Some studies have shown that reduced BIRC5 protein level can affect the development of ESCC by inhibiting the migration and invasion of tumor cells [19], but the mechanism of BIRC5 in the occurrence and development of ESCC is not clear. Also, whether the overexpression of BIRC5 is related to the differentiation, invasion depth, stage, and metastasis of ESCC cells remains controversial [17,19], Therefore, it is necessary to further research the clinical significance and biological role of BIRC5 expression in ESCC and explore the mechanism of the occurrence and development of ESCC.
In this study, by comprehensively using ESCC clinical samples and information from in-house, literature, and multiple databases, we attempted to explore BIRC5 expression in ESCC at both mRNA and protein levels. Furthermore, we have explored the clinical significance of BIRC5 expression in the prognosis and distinction of ESCC and the underlying molecular mechanism of the disease. From this, the research contributes to understanding of the pathogenesis of ESCC, and is conducive to identify ESCC and prediction in prognosis of the disease.
2.
Materials and methods
2.1. Tissue microarray and in-house BIRC5 protein expression assay
The current study with approval number 2018-KY-(0162), was approved by the Ethics Committee of the Second Affiliated Hospital of Guangxi Medical University, China. The inclusion criteria were: (1) patients diagnosed as ESCC by pathology; (2) patients with completely preserved paraffin specimens; (3) participants aged no less than 18 years; (4) patients with oral informed consent. Patients who did not agree to participate in or did not provide sufficient clinical information for subsequent research were excluded.
To understand the difference in BIRC5 protein expression between ESCC patients and non-ESCC patients, five formalin fixed paraffin embedded tissue microarray sections (ESC1504, ESC1503, ESC242, ESC241, and ESC481), purchased from Fanpu Biotech (Guilin, China), were used to explore BIRC5 protein expression. These sections were used for immunohistochemistry. First, the tissue sections were dewaxed and repaired. After that, 0.01 M citrate buffer (pH = 6.0) was used to soak tissue sections for antigen extraction. Then, the endogenous peroxidase was inactivated by 3% H2O2. Moreover, the tissue sections were placed in rabbit anti-human BIRC5 monoclonal antibody (ab126762, Abcam, UK, dilution ratio 1:100), while the negative control slides were placed in phosphate buffer, and then both were cultured overnight at 4℃. Furthermore, a second antibody labeled horseradish peroxidase (Changdao Biotechnology Co., Ltd., Shanghai, China) was added to the tissue sections, which would be heated to 25℃ for 25 minutes and dyed with 3, 3'-diaminobenzidine. Afterwards, tissue sections were dehydrated and sealed, and then evaluated under a microscope to observe the staining of the nucleus and/or cytoplasm, in which the positive staining was brown granules, and the negative staining was blue granules. In this study, the staining intensity score and the percentage record score were evaluated for each sample. Two senior pathologists individually and randomly selected 10 areas from each staining section for evaluation and recording. The staining intensity score was as follows: integer scores 0–3 respectively represented no, light, moderate, and strong staining; the percentage score was as follows: integer scores 0–4 respectively represented < 5%, 5%–25%, 26%–50%, 51%–75% and > 75% positive cells. The total protein expression score was obtained by multiplying the intensity score and the recorded percentage score. Supplementary Appendix 1 shows the clinical parameters of samples in tissue microarray.
2.2. Collection of public microarrays and RNA sequencing data sets of BIRC5 mRNA expression
To further explore the significance of BIRC5 expression in ESCC, through public databases and published articles, we collected microarrays, RNA sequencing data sets related to BIRC5 mRNA expression, and clinical information. The databases consulted by us include Gene Expression Omnibus, Oncomine, ArrayExpress, and the GDC Data Portal. The search terms used for microarrays and RNA sequencing data sets was: "esophag* and (mRNA or gene)". The included criteria for data sets were: (1) homo sapiens-related research. (2) ESCC-related tissues or cells; (3) data with mRNA expression of BIRC5. Data sets with duplicate, incomplete, or unqualified sample numbers were excluded. The steps for selecting data sets are listed in Supplementary Appendix 2.
2.3. Statistical analysis
2.3.1. Data processing
All data processing steps were executed in R software (v4.2.0). Expression of genes was first transformed by log2 (x + 1). Then, for those platforms (e.g., GPL20795) containing more than one data set, the Surrogate Variable Analysis software package [20] was used to reduce the batch effects. An example of before and after removing batch effects can be seen in Supplementary Appendix 3, and the data was clearly clustered together after removing the batch effects. The limma [21] or edge [22] software packages were used to standardize the gene expression value of each data set.
2.3.2. BIRC5 expression in ESCC and non-ESCC
In R (v4.2.0), the t test and the standardized mean difference (SMD) were used to analyze whether there was a significant difference in BIRC5 expression between the ESCC group and the non-ESCC control group. The fixed effect model and random effect model were used to calculate SMD. When I2 value (I2 test) < 50% or p value (chi-square test) > 0.1, which indicates that the heterogeneity of SMD is not obvious, the fixed effect model should be used. Otherwise, the random effect model was used. For SMD, that the 95% confidence interval (CI) does not include zero represents its statistical significance. The funnel plot of Begg's test was applied to assess publication bias, and there was no significant publication bias when p > 0.1. Violin plots and forests plots were drawn, respectively, based on ggplot2 and meta [23] packages.
2.3.3. Clinical significance of BIRC5 expression in ESCC
The Wilcoxon test was used to explore the correlation between the expression of BIRC5 and the clinical parameters of ESCC patients, including gender, age, survival status, tumor stage, lymph node stage, clinical stage, tumor location, alcohol use, and tobacco use. Receiver operating characteristic curves (ROCs) and summary receiver operating characteristic curve (sROC) were used to identify ESCC samples and non-ESCC samples based on area under the curve (AUC, ranging from zero to one). The closer the AUC is to one, the better the discrimination effect of BIRC5 expression. Also, the range and criteria of sensitivity and specificity are the same as AUC. For likelihood ratios, the bigger the positive likelihood ratio or the smaller the negative likelihood ratio, the increasing superior the discrimination effect of BIRC5 expression. The effect of BIRC5 expression on the prognosis of ESCC patients was explored by Kaplan-Meier survival curve and restricted mean survival time (RMST). All above plots were drawn in R (v4.2.0).
2.3.4. Potential mechanisms of BIRC5 in ESCC
The limma [21] software package was used to screen differential expression genes. A gene, with the absolute value of log2 (fold change) ≥ 1 and SMD > 0 (the 95% CI of SMD does not contain zero), was identified as a differential high-expression gene (DHEG). At the same time, positive correlation genes (Pearson coefficient ≥ 3, p < 0.05) of BIRC5 expression were screened. The positive correlation genes and DHEGs were crossed to obtain the BIRC5-related DHEGs (RDHEGs). To explore the potential molecular mechanism of BIRC5 in ESCC, Gene Ontology (GO) terms, Disease Ontology (DO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and Reactome pathway analyses were carried out based on RDHEGs. A term or a pathway with an adjusted p value < 0.05 was selected for the study. Also, the identification of hub genes in ESCC was accomplished by using a protein-protein interaction (PPI) network through RDHEGs.
To explore the mechanism of regulating BIRC5 expression, we predicted the transcription factors regulating BIRC5 expression using the following methods: First, we obtained the base sequence of the underlying promoter region of BIRC5 (1 kb upstream of the BIRC5 transcription start site) from the National Center for Biotechnology Information. Second, transcription factors (TFs) that may regulate BIRC5 expression were predicted from animal transcription factor (through the base sequences) [24] and Cistrome Data Browser [25] databases. Third, the initial TFs were obtained by the intersection of predicted TFs and RDHEGs. Fourth, matched sequences between the underlying promoter region of BIRC5 and motifs of initial TFs (obtained from the JASPAR database [26,27]) were explored via MEME-Suite [28]. Finally, the potential TFs regulating BIRC5 were identified from the initial TFs with chromatin immunoprecipitation sequencing (CHIP Seq) data in Cistrome Data Browser.
In this study, when a p value < 0.05, there was statistical significance. Supplementary Appendix 4 shows the main design of this study. To avoid confusion, in this study, the BIRC5 gene was presented in an italic form; the "BIRC5 protein" just represented the protein product of the gene BIRC5, while the "BIRC5" without "protein" represented both BIRC5 gene and BIRC5 protein.
3.
Results
3.1. Samples and their information of data sets included in the study
A total of 60 data sets (e.g., GSE70409) for BIRC5 mRNA expression and one in-house data set for BIRC5 protein expression were included for further research, containing 1,119 samples of the ESCC group and 1,004 samples of the non-ESCC group (Table 1). Then, the data sets with the same platform were reorganized into a new data set, and the batch effects were removed. For example, the new data set called "GPL13497" consists of data sets GSE97558 and GSE45168. A total of 16 new data sets were used for further research (Table 1).
3.2. BIRC5 expression in ESCC
Under the microscope, it can be seen that BIRC5 protein is mainly expressed in ESCC tissues (Figure 1), and upregulated expression of BIRC5 protein is detected in the ESCC group instead of the non-ESCC group via tissue microarrays (Figure 2, p < 0.05). Much more interestingly, consistent with the expression of protein level in ESCC, the results of almost every data set (except GSE45350) consistently show that the BIRC5 expression in the ESCC group was significantly higher than that in the non-ESCC group (Figure 2, p < 0.05).
Fifteen data sets were tested for heterogeneity, and the I2 = 86% (p value < 0.01) indicated significant heterogeneity, for which a random effect model was selected for calculating SMD. It can be seen from Figure 3A that BIRC5 mRNA and protein levels are both highly expressed in the ESCC group rather than the non-ESCC group (SMD = 2.21, 95% CI: 1.84 to 2.58), and Begg's funnel plot shows no significant publication bias in the SMD results, p = 0.729 (Figure 3B). These results indicate that BIRC5 is highly expressed in ESCC.
3.3. The significance of BIRC5 expression in ESCC
3.3.1. Relationship between clinical parameters and BIRC5 expression
There exist three data sets—in-house tissue microarrays, GSE53624, and the Cancer Genome Atlas (TCGA)—with clinical parameters. Thus, the relationship between clinical parameters and BIRC5 expression was explored through the three data sets. As shown in Supplementary Appendix 5, although no clinical parameter was found to be related to BIRC5 expression at both protein and mRNA levels, the mRNA expression of BIRC5 was higher in patients aged 65 or older and in non-lower esophageal tumor locations (p < 0.05, Supplementary Appendixes 5A–5B); also, BIRC5 protein expression was increased in higher lymph node stages (N1–N3) and clinical stages (Ⅲ-Ⅳ) (p < 0.05, Supplementary Appendix 5C).
3.3.2. Relationship between BIRC5 expression and prognosis of ESCC patients
Based on prognosis information from ESCC samples in the three data sets—in-house tissue microarrays, GSE53624, and TCGA—the median of the expression value of BIRC5 was used as the cut-off value for survival analyses. Through the Kaplan–Meier curves (Figure 4A), there is a cross between the high expression group and the low expression group for each data set. Thus, the RMST instead of the log-rank test was suitable for evaluating the relationship between BIRC5 expression and the prognosis of ESCC patients [29]. For in-house tissue microarrays, although patients in the high BIRC5 protein level group tended to have a shorter RMST, the difference was not statistically significant (p > 0.05, Figure 4B). For GSE53624, ESCC patients from the high BIRC5 expression group tended to have a shorter RMST than those from the low BIRC5 expression group (p < 0.05, Figure 4B). However, this was not consistent with the TCGA data set, with RMST of the high BIRC5 expression group = 3.23 years and that of the low expression of BIRC5 group = 2.18 years (p < 0.05, Figure 4B).
3.3.3. Effect of BIRC5 expression in distinguishing ESCC
Via ROCs and sROC, BIRC5 expression shows significant identifying ability between ESCC and non-ESCC samples, with all AUCs of ROCs > 0.82 (Supplementary Appendix 6) and AUC of sROC was 0.96 (95% CI: 0.93–0.97, Figure 5). The sensitivity and specificity were 0.92 (95% CI: 0.84–0.96) and 0.91 (95% CI: 0.87–0.94), respectively. The positive likelihood ratio and negative likelihood ratio were 10.71 (95% CI: 7.11–16.14) and 0.09 (95% CI: 0.05–0.17), respectively (Supplementary Appendix 7). These indicate that BIRC5 expression has an outstanding ability to distinguish ESCC samples from non-ESCC samples.
3.4. The potential mechanism of BIRC5 in ESCC
3.4.1. Identification of RDHEGs for exploring mechanism of BIRC5 in ESCC
A total of 8,857 genes with the absolute value of log2 (fold change) ≥ 1 were screened (Supplementary Appendix 8). Among them, 4,484 were selected with SMD > 0 and the 95% CI of SMD did not contain zero (data not shown). In the genes related to BIRC5 expression, 647 were identified to be positively correlated with BIRC5 expression in at least five data sets. Then, 445 RDHEGs were obtained by crossing DHEGs with BIRC5-positive correlation genes (data not shown).
3.4.2. Enrichment analyses for exploring mechanism of BIRC5 in ESCC
In GO analysis, RDHEGs significantly clustered in GO terms related to DNA replication, including "Chromosomal region" (cell components), "Nuclear division" (biological processes), and "Catalytic activity acting on DNA" (molecular functions) (Figure 6). In DO enrichment analysis, RDHEGs significantly involved in a variety of tumors, including "Hereditary breast ovarian cancer" and "Sensory system cancer" (Figure 6B).
In the main KEGG analysis, the top pathway of RDHEGs is the "Cell cycle, " while the Reactome enrichment pathway of RDHEGs is the "Cell cycle checkpoint" (Figure 7A). Both KEGG and Reactome analyses suggest that RDHEGs play an essential role in the cell cycle. Thus, we further performed PPI analysis based on 88 RDHEGs in the cell cycle checkpoint and cell cycle pathways. Then, CDK1 (cyclin dependent kinase 1), MAD2L, and CDC20 (cell cycle 20) are identified as hub genes in these RDHEGs (Figure 7B), indicating their important role in the cell cycle.
To explore the mechanism of high expression of BIRC5 in ESCC, TFs that may regulate BIRC5 expression were predicted. 422 and 100 TFs were predicted from Animal TFDB and Cistrome Data Browse, respectively. After the intersection of these TFs and RDHEGs, four TFs—E2F1, EZH2, MAZ, and TFDP1—were screened (Supplementary Appendix 9A). The four TFs were highly expressed in ESCC (SMDs > 0, 95% CI excluding 0) (Supplementary Appendix 9B) and positively correlated with the expression of BIRC5 (Supplementary Appendix 10). Among the four TFs, E2F1 has been reported to regulate BIRC5 expression [30], while no motif of EZH2 can be found in JASPAR. Thus, MAZ and TFDP1 were further analyzed. With the MEME-Suite analysis tool, among the upstream BIRC5 transcription initiation site (chr78213252–7821425) locate the binding sequences of MAZ and TFDP1 (Supplementary Appendix 11). Moreover, there were ChIP-Seq binding peaks of MAZ and TFDP1 in the potential promoter region upstream of the transcription start site of BIRC5, suggesting that these two TFs regulate the transcription of BIRC5 and affect its expression (Figure 8). Therefore, MAZ and TFDP1 may regulate BIRC5 expression.
4.
Discussion
In this study, to promote understanding of the pathogenesis of ESCC and provide a new idea for the clinical distinction and treatment of ESCC, we pay attention to the changes in the expression level of BIRC5 in ESCC, the clinical significance of BIRC5 expression level, and the underlying molecular mechanism of it in ESCC. By analyzing numerous samples (n = 2,123) from multiple research centers, with a t test and calculating SMD, we found that BIRC5 mRNA and protein are highly expressed in the ESCC group. The high expression of BIRC5 in ESCC is related to the clinical characteristics of ESCC patients, including age, lymph node stages, and clinical stages. Upregulated BIRC5 expression makes it feasible to distinguish ESCC from non-ESCC, suggesting its important clinical significance in ESCC. Also, with several analytical methods, this study also reveals that BIRC5 plays an essential role in ESCC and that it is mainly related to the cell cycle. The two TFs, MAZ and TFDP1, contribute to the upregulated expression of BIRC5 in ESCC; and as far as we know, no relevant reports have been found before this study, demonstrating the novelty of our research. In conclusion, this study demonstrated that elevated BIRC5 expression may have important clinical significance in ESCC. Its molecular mechanism in ESCC was initially discussed, which contributes to the understanding of ESCC in pathogenesis.
ESCC is one of the most common types of malignant tumors in humans. Most ESCC patients were diagnosed at an advanced stage, and their prognosis is not ideal [31]. Therefore, exploring new markers for distinguishing ESCC from non-ESCC and predicting prognosis of ESCC patients contributes to improving the survival rate of patients [32]. As a member of the apoptosis-suppressing gene family, BIRC5 is of great significance in a variety of malignant tumors. For instance, upregulated BIRC5 expression can be seen in multiple cancers, such as non-small cell lung cancer [15] and breast cancer [16], and it may become a marker for early detection and/or treatment for these two cancers. Elevated BIRC5 expression has been revealed in ESCC based on previous reports [15,33], but some limitations are generally common in these reports, such as small sample size and focusing on a single level (protein) of BIRC5 expression. In order to further verify the expression of BIRC5 in ESCC, this study, first based on more than 2,000 samples from multiple research centers and various statistical methods, disclosed that the mRNA and protein expression levels of BIRC5 were upregulated in ESCC.
The high expression of BIRC5 may be related to the poor prognosis of ESCC patients, and it makes it possible to distinguish ESCC from non-ESCC. Regarding the clinical significance of BIRC5 in ESCC, Shang et al. [19] found that as ESCC progressed, the expression level of BIRC5 increased. However, whether BIRC5 expression is correlated with the tumor size, lymph node metastasis, and advanced clinical stage of ESCC patients is still controversial [15,17]. In this study, we have not yet observed that BIRC5 is related to several clinical parameters—gender, survival status, tumor size, distant metastasis, alcohol use, and tobacco consumption —of ESCC patients at both the mRNA and protein levels. However, high mRNA expression of BIRC5 can be detected in patients with the location of ESCC lesions (lower esophagus) and an older age (≥ 65 years old). The protein expression of BIRC5 is increased in the higher lymph node stage (N1–N3) and late clinical stages (Ⅲ-Ⅳ). These results suggest that BIRC5 may be associated with the prognosis of ESCC. However, according to previous reports, with BIRC5 expression, ESCC patients has better [34] or worse [35] prognosis is still controversial. Interestingly, such a contradiction can be seen between mRNA and protein expression levels of BIRC5 in this study. Based on samples from TCGA, ESCC patients with high expression of BIRC5 had longer RMST, which is contrary to the result of GSE53624. However, we tended to believe that BIRC5 expression was associated with a worse prognosis of ESCC patients for these reasons: (1) The sample in the GSE53624 data set (n = 119) is larger than the sample in the TCGA data set (n = 81), and the result of a larger sample size is more reliable. (2) The prognostic data of the GSE53624 data set is more orderly than that of TCGA (in the TCGA data set, before the 2-year survival time, the prognosis of patients with high expression of BIRC5 is worse, while after 2-year survival time, patients with upregulated expression of BIRC5 show better prognosis), so results based on GSE53624 are more credible. (3) The high expression of BIRC5 protein is related to advanced lymph node and clinical stages. Therefore, current studies suggest that elevated BIRC5 expression may be a risk factor for poor prognosis in ESCC patients. It is also worth noting that, based on the expression data of BIRC5, the samples of the ESCC group can be well distinguished from the samples of the non-ESCC group, suggesting the potential of BIRC5 as an ESCC discrimination indicator; a finding that, as far as we know, has not been reported before. In summary, the current study indicates that BIRC5 may have potential significance in the clinical application of ESCC, such as treatment and early discrimination of the disease.
To improve the understanding of the role of BIRC5 in ESCC, we also analyzed the molecular mechanism of BIRC5 in ESCC. The expression products of RDHEGs screened based on BIRC5 were mainly distributed in the cell chromosomal region, mainly involved in the process of cell nuclear division and DNA replication and participate in catalyzing DNA activity. The important period of DNA replication is the S phase of the interphase (interphase is one phase of the cell cycle, and the other is the division phase). Interestingly, although the BIRC5 (also known as survivin) protein encoded by BIRC5 mainly aggregates in G2 during cell division, it has been localized in the nucleus during the S phase [36], suggesting the possibility of BIRC5 participating in DNA replication. Also, BIRC5 expression may not only be related to cell proliferation (DNA replication), but also affect the progression of a variety of tumors (such as ovarian cancer and sensory system cancer) based on DO analysis, implying the key role of BIRC5 in ESCC. Further analysis revealed that RDHEGs significantly aggregated in the cell cycle pathways of KEGG and Reactome, indicating that the mechanism of BIRC5 in ESCC may be related to the cell cycle. According to previous studies, BIRC5 protein is essential for cell division and can inhibit cell death [37]. The mechanism remains unclear, and the cell cycle may be one of the factors involved. In addition to the interphase shown above, the division phase may also be the stage where BIRC5 protein participates in the cell cycle process. During the division phase, BIRC5 protein helps the chromosomes to be arranged correctly in the later stage by positioning the chromosomal passenger complex at the centromere of the ante-mid period of mitosis [37]. Thus, BIRC5 may participate in cell proliferation by maintaining the cell cycle. In short, RDHEGs screened based on BIRC5 are mainly involved in the cell cycle, especially DNA replication, suggesting that BIRC5 may participate in the occurrence and development of ESCC through these pathways.
With PPI network analysis of RDHEGs enriched in the cell cycle pathway, CDK1, CDC20, and MAD2L were found to be the hub genes of these RDHEGs in the cell cycle pathway. CDK1, a serine/threonine protein kinase that acts on the G2/M point of the cell cycle, is a necessary condition for eukaryotic cell division and thus participates in cell proliferation. Abnormal regulation of CDK1 leads to abnormal cell differentiation and cell cycles, and eventually leads to malignant tumor formation [38]. Like CDK1, CDC20 is overexpressed in a variety of human tumors and is considered to show a carcinogenic effect in human tumorigenesis [39]. Its expression was associated with poor prognosis in patients with breast cancer [40] and colorectal cancer [41]. Another hub gene, MAD2L, was thought to be associated with gastric cancer progression by participating in the cell cycle [42], although it is rarely reported in tumors. As RDHEGs screened based on BIRC5, the three genes – CDK1, CDC20, and MAD2L all play a role in promoting cancer and are also closely related to the cell cycle. Therefore, it is worth further exploring whether BIRC5 plays a role in ESCC by interacting with these three genes.
In this study, we also found that MAZ and TFDP1 may be transcription factors regulating BIRC5 expression. Studies have shown that MAZ is overexpressed in prostate cancer tissues with bone metastasis and further enhanced in metastatic bone tissues [43]. The TFDP1 gene was amplified in non-small cell lung cancer [44], and was involved in the amplification of hepatocellular carcinoma by up-regulating the expression of CCNE1 [45]. These findings suggest that MAZ and TFDP1 are involved in the expression regulation of cancer-related genes. In this study, these two TFs were not only upregulated in ESCC (similar to BIRC5), but also significantly positively correlated with the expression of BIRC5 in ESCC. Furthermore, for MAZ and TFDP1, there were ChIP-Seq binding peaks in the potential promoter region of BIRC5. These results suggest that the transcription and expression of BIRC5 in ESCC may be caused by the regulation of MAZ and TFDP1, but this needs further study.
Although we have revealed some noteworthy findings, this study still has limitations: First, all the samples we collected are from tissues or cells, and we failed to collect serum samples from ESCC patients to study the clinical significance of BIRC5 in ESCC, such as whether it is possible to quickly screen ESCC by detecting serum BIRC5 levels. At the same time, this study lacks samples with detailed prognostic parameters, so it is impossible to carry out prognostic-related research through large samples. Also, in the future, the regulation of MAZ and TFDP1 on BIRC5 should be verified through in vivo and in vitro experiments.
5.
Conclusions
In summary, the current studies have revealed that the mRNA and protein of BIRC5 are highly expressed in the ESCC group, and the highly expressed BIRC5 may be used as a potential marker for discrimination and treatment of ESCC. Also, this study also explored the potential mechanism of BIRC5 in ESCC and found that the cell cycle may be an important way for BIRC5 to participate in the occurrence and development of ESCC, and the regulation of MAZ and TFDP1 may be a mechanism for the high expression of BIRC5 in ESCC.
Acknowledgments
This work was supported by the College Student Innovation and Entrepreneurship Training Program Project of the First Clinical Medical University of Guangxi Medical University (2020YFYA05, 2020YFYA06), the College Student Innovation and Entrepreneurship Training Program Project (202010598052) and the Guangxi Medical University Future Academic Stars Project (WLXSZX20081). The results shown in the study are in part based upon data generated by the TCGA Research Network: www.cancer.gov/tcga.
Conflict of interest
The authors have no conflicts of interest to declare.