Bioinformatics Analysis of ZBTB16 as a Prognostic Marker for Ewing’s Sarcoma

Ding, Ke; Qiu, Wenli; Yu, Dianbo; Ma, Huade; Xie, Kangqi; Luo, Fuqiang; Li, Shanlang; Li, Zaiyong; Wei, Jihua

doi:https://doi.org/10.1155/2021/1989917

BioMed Research International

On this page

Abstract Introduction Materials Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Accelerating Precision Medicine through Integrated Analysis of Multiscale Biological Data

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 1989917 | https://doi.org/10.1155/2021/1989917

Bioinformatics Analysis of ZBTB16 as a Prognostic Marker for Ewing’s Sarcoma

Ke Ding,¹Wenli Qiu,²Dianbo Yu,¹Huade Ma,¹Kangqi Xie,¹Fuqiang Luo,¹Shanlang Li,¹Zaiyong Li,¹and Jihua Wei¹

Academic Editor: Jianxin Shi

Received16 Aug 2021

Revised03 Sept 2021

Accepted04 Sept 2021

Published06 Oct 2021

Abstract

Objective. The purpose of this study is to identify novel biomarkers for the prognosis of Ewing’s sarcoma based on bioinformatics analysis. Methods. The GSE63157 and GSE17679 datasets contain patient and healthy control microarray data that were downloaded from the Gene Expression Omnibus (GEO) database and analyzed through R language software to obtain differentially expressed genes (DEGs). Firstly, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment, protein-protein interaction (PPI) networks, and Cytoscape Molecular Complex Detection (MCODE) plug-in were then used to compute the highest scores of the module. After survival analysis, the hub genes were lastly obtained from the two module genes. Results. A total of 1181 DEGs were identified from the two GSEs. Through MCODE and survival analysis, we obtain 53 DEGs from the module and 29 overall survival- (OS-) related genes. ZBTB16 was the only downregulated gene after Venn diagrams. Survival analysis indicates that there was a significant correlation between the high expression of ZBTB16 and the OS of Ewing’s sarcoma (ES), and the low expression group had an unfavorable OS when compared to the high expression group. Conclusions. High expression of ZBTB16 may serve as a predictor biomarker of poor prognosis in ES patients.

1. Introduction

Ewing’s sarcoma (ES) is the second most common sarcoma of bone and soft tissue, usually occurring in children and adolescents [1]. The typical genetic feature is EWSR1 gene translocation, which leads to EWS-FLI1 gene fusion in t(11;22) [2], accounting for 85% of the ES family [3, 4], and is dominated by small round cells [5]. The early data report of the National Cancer Institute showed that the incidence of ES differed by gender, age, and race. Males accounted for 60.3%, mostly under 24 years old, and the incidence of whites was nine times that of African-Americans [6].

At present, optimal treatment of ES emphasizes multimodal therapy [7–9]. Chemotherapy is an efficient treatment for ES, and it significantly improves the 5-year overall survival (OS) rate [10, 11]. However, due to the lack of early efficient diagnosis, the high recurrence rate and distant metastasis of patients with ES lead to a poor prognosis [12], and the prognosis of patients with ES is not optimistic. Patients with ES showed poor survival in years, tumor cm, receiving radiotherapy alone, and receiving no treatment [13]. Therefore, there is still an urgent need for new diagnoses and treatment strategies to improve long-term survival.

With the development of tumor genomics research such as high-throughput sequencing, it is of great significance to find potential prognostic treatment pathways and targets for ES. In this study, our purpose is to investigate the biomarkers related to the prognosis of ES. We processed and analyzed microarray data from GSE63157 and GSE17679 datasets from the Gene Expression Omnibus (GEO) database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional and pathway enrichment analysis, protein-protein interaction (PPI) network construction of pivot genes, and identification of key genes that are significantly related to OS in ES are aimed at providing therapeutic strategies for the prognosis of ES.

2. Materials and Method

2.1. Data Source and Data Processing

Based on the GEO database (https://www.ncbi.nlm.nih.gov/geo/), we obtain two microarray datasets of ES: GSE63157 and GSE17679. GSE63157 was submitted by Volchenboum SL and Andrade J et al. and GSE17679 by Savola S and Klami A et al. GSE63157 included 85 ES patient samples; GSE17679 consisted of 88 ES patient samples and 18 matched human normal skeletal muscle samples (10 ES cell line samples were deleted). The R GEOquery package (version 2.60.0) does the download annotation. GSE63157 was sampled based on the GPL5715 Affymetrix Human Exon 1.0 ST Array. GSE17679 was sampled based on the GPL570 Affymetrix Human Genome U133 Plus 2.0 Array.

2.2. Identification of DEGs

All the samples’ data of GSE63157 and GSE17679 were normalized and visualized by the R “Limma” package (version 3.48.1) and also the annotation. The screening conditions for the DEGs were an absolute value of and adj. value < 0.05.

2.3. GO and KEGG Enrichment Analysis

The R “Clusterprofile” package (version 4.0.0) was used for GO and KEGG enrichment analysis and visualization. GO analysis was used to annotate the biological process (BP), cytological component (CC), and molecular function (MF) of DEGs, and KEGG enrichment analysis was used to understand the related signaling pathways; was statistically significant.

2.4. PPI Network Construction and Module Analysis

The online tool STRING (http://string-db.org) (version 11.0) was used to construct the PPI network of DEGs. The threshold for statistically significant interaction is a confidence . Cytoscape (version 3.6.1) is software providing an open platform for the visualization of biological processes, molecular interaction networks, and integration of these networks. Molecular Complex Detection (MCODE) (version 1.5.1), a plug-in of Cytoscape, is applied to discover densely connected regions and identify the significant modules through clustering a given network. MCODE , , max , node score , and degree were set as the criteria for the selection of significant modules.

2.5. Survival Analysis of Hub Gene

The survival data of 173 ES patients samples were integrated with the expression data of DEGs. Univariate and multivariate Cox regression analyses were performed through using the R “Survival” package to identify the DEGs associated with the OS time and clinical characteristics of ES, and the hub genes that could independently guide the prognosis were found by taking the overlap with module genes. The ES patients’ samples were divided into the high expression groups and the low expression groups. The difference between the hub genes and ES was further verified by the Kaplan-Meier survival curve. A Venn diagram was used for visualization, and was statistically significant.

3. Results

3.1. Normalization of Total DEGs

Gene expression data of GSE63157 and GSE17679 were standardized and displayed in a boxplot. The black line in (Figure 1(b)) is basically at the same level, indicating a high consistency of data.

(a)

(b)

3.2. Identification of DEGs

A total of 1181 DEGs were identified from GSE63157 and GSE17679 datasets, of which 364 were upregulated genes and 817 were downregulated genes (Figure 2(a)). Volcano plot of the DEGs, respectively, suggests that the expression of identified DEGs can correctly distinguish the case and normal samples.

3.3. GO and KEGG Enrichment Analysis

The R package “Clusterprofile” was used for functional enrichment analysis to assess the biological classification of DEGs (Figure 3). GO analysis showed (Figures 3(a) and 3(b)) that among biological processes, DEGs were significantly involved in the muscle system process (110 genes), muscle contraction (97 genes), energy derivation by oxidation of organic compounds (84 genes), actin-mediated cell contraction (49 genes), and muscle filament sliding (31 genes). In the cell component, DEGs were significantly involved in contractile fiber (97 genes), myofibril (94 genes), sarcomere (87 genes), I band (58 genes), and Z disc (50 genes). In molecular function, DEGs are significantly involved in binding-related projects. The molecular function includes actin binding (82 genes), structural constituent of muscle (25 genes), electron transfer activity (42 genes), actin filament binding (43 genes), and NADH dehydrogenase activity (18 genes). In addition, KEGG analysis showed (Figures 3(c) and 3(d)) that DEGs were mainly involved in diabetic cardiomyopathy (57 genes), cardiac muscle contraction (37 genes), and pathways of neurodegeneration-multiple diseases (69 genes) (Table 1).

(a)

(b)

(c)

(d)

Figure 3

GO and KEGG pathway analysis of the DEGs. The color depth of nodes refers to the -value. The size of nodes refers to the number of genes. (a) Barplot of top 5 significantly enriched biological processes, cell component, and molecular function in DEGs. (b) Bubble plot of top 5 significantly enriched biological processes, cell component, and molecular function in DEGs. (c) Barplot of top 15 significantly enriched molecular functions in DEGs. (d) Bubble plot of top 15 significantly enriched KEGG pathways in DEGs.

3.4. PPI Network Construction and Module Analysis

With a total of 645 nodes and 3300 edges, the PPI network information of the DEGs from STRING was imported into Cytoscape, and the module of DEGs was constructed. We just take the module with the highest MCODE score. This module contains 29 nodes and 406 edges (Figure 4(a)).

(a)

(b)

(c)

3.5. Hub Gene Screening and Identification

Survival analysis showed that there was a total of 72 DEGs related to OS time of patients with ES (Figure 4(b)), among which 53 genes were related to clinical characteristics. The hub gene was obtained after the intersection with 29 genes in the module with the highest MCODE score (Figure 4(c)).

3.6. Survival Analysis of Hub Gene

A total of 173 ES patients have follow-up information for survival analysis. All the ES patients were divided into the high expression group () and the low expression group (), according to the median expression levels of ZBTB16. The expression of ZBTB16 was downregulated in the ES patients group, and the OS time of patients with low expression was poorer than the high expression group. Univariate and multivariate Cox regression analyses showed that ZBTB16 may be a prognostic factor for ES. The difference is statistically significant () (Figure 5).

(a)

(b)

(c)

(d)

4. Discussion

Ewing’s sarcoma (ES) is the second most common malignant bone tumor in children and adolescents [2]. At present, vincristine, doxorubicin, cyclophosphamide, ifosfamide, and etoposide are the standard chemotherapy regimens. However, the prognosis of patients is not optimistic [14, 15] and has little effect on patients with metastasis or recurrence [16]. Microarray technology provides a new direction for the treatment of patients with ES. Therefore, it is necessary to explore appropriate and effective measures to guide the prognosis of treatment.

The purpose of our microarray-based analysis is to identify novel biomarkers related to the prognosis of patients with ES. Based on GSE63157 and GSE17679 datasets, 1181 DEGs were screened and identified from the ES patient group and normal control group, including 364 upregulated genes and 817 downregulated genes. Then, the functional enrichment analysis was carried out, and the interaction between DEGs was discussed. GO and KEGG pathway enrichment analysis shows that DEGs are mainly involved in the muscle system process, muscle contract, structural constituent of muscle, electron transfer activity, diabetic cardiomyopathy, and cardiac muscle contraction. These biological functions are consistent with our knowledge of ES. Based on the MCODE plug-in computing the highest score module and survival analysis filtering, we chose the overlapping one gene as the hub gene according to the two modules.

Zinc-finger and BTB domain-containing 16 (ZBTB16), also known as promyelocytic leukemia zinc finger protein (PLZF), was first found in t(11;17) chromosome translocation, and fusion with retinoic acid receptor α (RARa) was identified [17]. ZBTB16 is an important mediator that regulates the expression of neuropeptides [18], ion channels [19], or various receptors [20], to influence the response to incoming signals and/or downstream neuron signals. ZBTB16 is expressed in various tissues and plays a variety of biological functions through the regulation of tissue-specific target genes, which is related to adipogenesis [21], myocardial remodeling [22], osteoblast differentiation [23], and obesity [24]. The enhanced expression of ZBTB16 can increase the number of mitochondria and the ability of respiration and uncoupling. Besides, ZBTB16 is associated with various tumors. ZBTB16 is a classic androgen receptor (AR) regulatory gene, which has antiproliferative activity in prostate cancer cells [25] and plays a role in inhibiting cancer by inhibiting the MAPK pathway [26]. Secondly, ZBTB16 may become a prognostic marker for breast cancer and hepatocellular carcinoma [27, 28].

Up to now, there is no literature report on the relationship between ZBTB16 and the prognosis of Ewing’s sarcoma. Previous studies have found that ZBTB16 can be used to distinguish leiomyosarcomas from leiomyomas [29]. In addition, pathological studies have reported that ZBTB16 has high sensitivity and specificity in various tissue types of Primitive Neurological Tumor (PNET) [30]. This study verified the correlation between ZBTB16 and ES, and we need to further explore the relationship between ZBTB16 and the prognosis of ES. In our study, the expression of the ZBTB16 was downregulated in the ES group when compared to the control group. Survival analysis indicates that the high expression of the ZBTB16 was significantly higher than that in the low expression group, which further clarified that ZBTB16 was related to the prognosis of ES.

In a word, we confirmed that ZBTB16 provides valuable information for the prognosis of ES and may become a novel biomarker for the prognosis of ES patients. However, there are some limitations to this study. Due to the limited raw data, we did not consider the associations between ZBTB16 expression and clinicopathological factors (e.g., stage and grade) of patients with ES. Further research is needed to make up for the shortcomings of this article.

5. Conclusion

We identified ZBTB16 from the GEO database through bioinformatics analysis. Moreover, higher ZBTB16 expression patients can have a longer OS time. It may be the potential prognostic biomarker of ES. In addition, biological functions and mechanisms need to be studied further in ES.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Ding Ke and Qiu Wenli contributed equally to this work.

Acknowledgments

The study was supported by the (1) Guangxi Natural Science Foundation (Contract No. 2020GXNSFAA259068) and (2) First Batch of High-level Talent Scientific Research Projects of the Affiliated Hospital of Youjiang Medical University for Nationalities in 2019 (Contract No. R20196321).

References

Y. Park, M. Chang, K.-M. Kim, and D.-H. An, “Analysis of prognostic factors in extraosseous Ewing sarcoma family of tumors: review of St. Jude Children's Research Hospital experience,” Annals of Surgical Oncology, vol. 19, no. 12, pp. 3816–3822, 2012.
View at: Publisher Site | Google Scholar
H. E. Grier, “THE EWING FAMILY OF TUMORS: Ewing's Sarcoma and Primitive Neuroectodermal Tumors,” Pediatric Clinics of North America, vol. 44, no. 4, pp. 991–1004, 1997.
View at: Publisher Site | Google Scholar
M.-C. Le Deley, O. Delattre, K.-L. Schaefer et al., “Impact of EWS-ETS fusion type on disease progression in Ewing's sarcoma/peripheral primitive neuroectodermal tumor: prospective results from the cooperative Euro-E.W.I.N.G. 99 trial,” Journal of Clinical Oncology, vol. 28, no. 12, pp. 1982–1988, 2010.
View at: Publisher Site | Google Scholar
O. Delattre, J. Zucman, B. Plougastel et al., “Gene fusion with an _ETS_ DNA-binding domain caused by chromosome translocation in human tumours,” Nature, vol. 359, no. 6391, pp. 162–165, 1992.
View at: Publisher Site | Google Scholar
M. Sbaraglia, A. Righi, M. Gambarotti, and A. P. Dei Tos, “Ewing sarcoma and Ewing-like tumors,” Virchows Archiv, vol. 476, no. 1, pp. 109–119, 2020.
View at: Publisher Site | Google Scholar
M. U. Jawad, M. C. Cheung, E. S. Min, M. M. Schneiderbauer, L. G. Koniaris, and S. P. Scully, “Ewing sarcoma demonstrates racial disparities in incidence-related and sex-related differences in outcome: an analysis of 1631 cases from the SEER database, 1973-2005,” Cancer, vol. 115, no. 15, pp. 3526–3536, 2009.
View at: Publisher Site | Google Scholar
R. W. Rödl, C. Hoffmann, G. Gosheger, B. Leidinger, H. Jürgens, and W. Winkelmann, “Ewing's sarcoma of the pelvis: combined surgery and radiotherapy treatment,” Journal of Surgical Oncology, vol. 83, no. 3, pp. 154–160, 2003.
View at: Publisher Site | Google Scholar
J. Alvarez-SanNicolas, I. Gracia-Alegria, L. Trullols-Tarrago, A. Peiro-Ibañez, and C. Lamas-Gomez, “Prognostic factors and survival in Ewing's sarcoma treated by limb salvage surgery,” Clinical & Translational Oncology, vol. 21, no. 10, pp. 1374–1382, 2019.
View at: Publisher Site | Google Scholar
K. Jamshidi, M. R. Shirazi, A. Bagherifard, and A. Mirzaei, “Curettage, phenolization, and cementation in paediatric Ewing's sarcoma with a complete radiological response to neoadjuvant chemotherapy,” International Orthopaedics, vol. 43, no. 2, pp. 467–473, 2019.
View at: Publisher Site | Google Scholar
J. L. Pretz, C. M. Barysauskas, S. George et al., “Localized adult Ewing sarcoma: favorable outcomes with alternating vincristine, doxorubicin, cyclophosphamide, and ifosfamide, etoposide (VDC/IE)-based multimodality therapy,” The Oncologist, vol. 22, no. 10, pp. 1265–1270, 2017.
View at: Publisher Site | Google Scholar
J. E. Oesterheld, D. R. Reed, B. A. Setty et al., “Phase II trial of gemcitabine and nab-paclitaxel in patients with recurrent Ewing sarcoma: a report from the National Pediatric Cancer Foundation,” Pediatric Blood & Cancer, vol. 67, no. 7, 2020.
View at: Publisher Site | Google Scholar
A. Mendoza‐Naranjo, A. El‐Naggar, D. H. Wai et al., “ERBB4 confers metastatic capacity in Ewing sarcoma,” EMBO Molecular Medicine, vol. 5, no. 7, pp. 1087–1102, 2013.
View at: Publisher Site | Google Scholar
Z. H. Wan, Z. H. Huang, and L. B. Chen, “Survival outcome among patients with Ewing's sarcoma of bones and joints: a population-based cohort study,” São Paulo Medical Journal, vol. 136, no. 2, pp. 116–122, 2018.
View at: Publisher Site | Google Scholar
M. Chin, R. Yokoyama, M. Sumi et al., “Multimodal treatment including standard chemotherapy with vincristine, doxorubicin, cyclophosphamide, ifosfamide, and etoposide for the Ewing sarcoma family of tumors in Japan: results of the Japan Ewing Sarcoma Study 04,” Pediatric Blood & Cancer, vol. 67, no. 5, article e28194, 2020.
View at: Publisher Site | Google Scholar
H. E. Grier, M. D. Krailo, N. J. Tarbell et al., “Addition of ifosfamide and etoposide to standard chemotherapy for Ewing’s sarcoma and primitive neuroectodermal tumor of bone,” The New England Journal of Medicine, vol. 348, no. 8, pp. 694–701, 2003.
View at: Publisher Site | Google Scholar
J. I. Albergo, C. L. L. Gaston, M. C. Parry et al., “Risk analysis factors for local recurrence in Ewing's sarcoma: when should adjuvant radiotherapy be administered?” Bone Joint J, vol. 100-B, no. 2, pp. 247–255, 2018.
View at: Publisher Site | Google Scholar
Z. Chen, N. J. Brand, A. Chen et al., “Fusion between a novel Krüppel-like zinc finger gene and the retinoic acid receptor-alpha locus due to a variant t(11;17) translocation associated with acute promyelocytic leukaemia,” The EMBO Journal, vol. 12, no. 3, pp. 1161–1167, 1993.
View at: Publisher Site | Google Scholar
H. Cheng, S. J. Pablico, J. Lee, J. S. Chang, and S. Yu, “Zinc finger transcription factor Zbtb16 coordinates the response to energy deficit in the mouse hypothalamus,” Frontiers in Neuroscience, vol. 14, 2020.
View at: Publisher Site | Google Scholar
A. Náray-Fejes-Tóth, C. Boyd, and G. Fejes-Tóth, “Regulation of epithelial sodium transport by promyelocytic leukemia zinc finger protein,” American Journal of Physiology. Renal Physiology, vol. 295, no. 1, pp. F18–F26, 2008.
View at: Publisher Site | Google Scholar
T. Zhang, K. Dong, W. Liang et al., “G-protein-coupled receptors regulate autophagy by ZBTB16-mediated ubiquitination and proteasomal degradation of Atg14L,” eLife, vol. 4, article e06734, 2015.
View at: Publisher Site | Google Scholar
M. A. Ambele, C. Dessels, C. Durandt, and M. S. Pepper, “Genome-wide analysis of gene expression during adipogenesis in human adipose- derived stromal cells reveals novel patterns of gene expression during adipocyte differentiation,” Stem Cell Research, vol. 16, no. 3, pp. 725–734, 2016.
View at: Publisher Site | Google Scholar
E. Farrell, A. E. Armstrong, A. C. Grimes, F. J. Naya, W. J. de Lange, and J. C. Ralphe, “Transcriptome analysis of cardiac hypertrophic growth in MYBPC3-null mice suggests early responders in hypertrophic remodeling,” Frontiers in Physiology, vol. 9, p. 1442, 2018.
View at: Publisher Site | Google Scholar
F. Marofi, G. Vahedi, S. Solali et al., “Gene expression of TWIST1 and ZBTB16 is regulated by methylation modifications during the osteoblastic differentiation of mesenchymal stem cells,” Journal of Cellular Physiology, vol. 234, no. 5, pp. 6230–6243, 2019.
View at: Publisher Site | Google Scholar
C. L. Plaisier, B. J. Bennett, A. He et al., “Zbtb16 has a role in brown adipocyte bioenergetics,” Nutrition & Diabetes, vol. 2, no. 9, article e46, 2012.
View at: Publisher Site | Google Scholar
F. Jiang and Z. Wang, “Identification and characterization of PLZF as a prostatic androgen-responsive gene,” Prostate, vol. 59, no. 4, pp. 426–435, 2004.
View at: Publisher Site | Google Scholar
C.-L. Hsieh, G. Botta, S. Gao et al., “PLZF, a tumor suppressor genetically lost in metastatic castration-resistant prostate cancer, is a mediator of resistance to androgen deprivation therapy,” Cancer Research, vol. 75, no. 10, pp. 1944–1948, 2015.
View at: Publisher Site | Google Scholar
J. He, M. Wu, L. Xiong et al., “BTB/POZ zinc finger protein ZBTB16 inhibits breast cancer proliferation and metastasis through upregulating ZBTB28 and antagonizing BCL6/ZBTB27,” Clinical Epigenetics, vol. 12, no. 1, p. 82, 2020.
View at: Publisher Site | Google Scholar
A. Hui, H. Lau, C. Cao, J. Zhou, P. Lai, and S. Tsui, “Downregulation of PLZF in human hepatocellular carcinoma and its clinical significance,” Oncology Reports, vol. 33, no. 1, pp. 397–402, 2015.
View at: Publisher Site | Google Scholar
M. Momeni, T. Kalir, S. Farag et al., “Immunohistochemical detection of promyelocytic leukemia zinc finger and histone 1.5 in uterine leiomyosarcoma and leiomyoma,” Reproductive Sciences, vol. 21, no. 9, pp. 1171–1176, 2014.
View at: Publisher Site | Google Scholar
G. Q. Xiao, A. E. Sherrod, and K. M. Hurth, “ZBTB16: A new biomarker for primitive neuroectodermal tumor element / Ewing sarcoma,” Pathology, Research and Practice, vol. 215, no. 10, p. 152536, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Ke Ding et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

349

Downloads

983

Citations