
==== Front
Mol Biomed
Mol Biomed
Molecular Biomedicine
2662-8651
Springer Nature Singapore Singapore

39227479
198
10.1186/s43556-024-00198-8
Research
SOX7 inhibits the malignant progression of bladder cancer via the DNMT3B/CYGB axis
Zhang Jingcheng 12
Zhang Wentao 12
Liu Ji 12
Liu Yuchao 12
Jiang Yufeng 23
Ainiwaer Ailiyaer 124
Chen Hanyang 12
Gu Zhuoran 12
Chen Haotian 12
Mao Shiyu 12
Guo Yadong 12
Xu Tianyuan 12
Xu Yunfei xuyunfeibb@sina.com

12
Wu Yuan 13956051264@163.com

5
Yao Xudong yxd@tongji.edu.cn

12
Yan Yang 1500069@tongji.edu.cn

12
1 grid.412538.9 0000 0004 0527 0050 Department of Urology, School of Medicine, Shanghai Tenth People’s Hospital, Tongji University, Shanghai, China
2 https://ror.org/03rc6as71 grid.24516.34 0000 0001 2370 4535 Urologic Cancer Institute, School of Medicine, Tongji University, Shanghai, China
3 grid.412538.9 0000 0004 0527 0050 Department of Urology, Chongming Branch, School of Medicine, Shanghai Tenth People’s Hospital, Tongji University, Shanghai, China
4 Department of Urology, Xinjiang Uygur Autonomous Region, Kashgar Prefecture Second People’s Hospital, Kashgar, China
5 https://ror.org/034t30j35 grid.9227.e 0000 0001 1957 3309 Department of Urology, Hefei Cancer Hospital, Chinese Academy of Sciences, Hefei, China
4 9 2024
4 9 2024
12 2024
5 3627 3 2024
22 7 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Bladder cancer (BCa) stands out as a highly prevalent malignant tumor affecting the urinary system. The Sex determining region Y-box protein family is recognized for its crucial role in BCa progression. However, the effect of Sex determining region Y-box 7 (SOX7) on BCa progression has not been fully elucidated. Herein, RNA-sequencing, western blot (WB), immunohistochemistry (IHC), immunofluorescence (IF) and tissue microarray were utilized to assess SOX7 expression in vitro and in vivo. Additionally, SOX7 expression, prognosis, and SOX7 + cytoglobin (CYGB) score were analyzed using R software. In vitro and vivo experiments were performed with BCa cell lines to validate the effect of SOX7 knockdown and overexpression on the malignant progression of BCa. The results showed that SOX7 exhibits low expression in BCa. It functions in diverse capacities, inhibiting the proliferative, migratory, and invasive capabilities of BCa. In addition, the experimental database demonstrated that SOX7 binds to the promoter of DNA methyltransferase 3 beta (DNMT3B), leading to the transcriptional inhibition of DNMT3B. This subsequently results in a reduced methylation of CYGB promoter, ultimately inhibiting the tumor progression of BCa. SOX7 + CYGB scores were significantly linked to patient prognosis. In conclusion, SOX7 inhibits the malignant progression of BCa via the DNMT3B/CYGB axis. Additionally, the SOX7 + CYGB score is capable of predicting the prognostic outcomes of BCa patients. Therefore, SOX7 and CYGB may play an important role in the progression of bladder cancer, and they can be used as prognostic markers of bladder cancer patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43556-024-00198-8.

Keywords

Bladder cancer (BCa)
Sex determining region Y-box 7 (SOX7)
Methylation
Cytoglobin (CYGB)
issue-copyright-statement© Sichuan International Medical Exchange & Promotion Association 2024
==== Body
pmcIntroduction

Bladder cancer (BCa) ranks among the most prevalent malignancies of the urinary system, with around 580,000 newly diagnosed cases worldwide in 2020 [1]. According to estimations by the World Health Organization, this number is expected to double by 2040 [2]. About 75% of BCa is non-muscle invasive BCa (NMIBC) [3]. The primary treatment approach involves bladder-preserving transurethral resection of bladder tumors combined with postoperative bladder infusion chemotherapy [4, 5]. Despite these treatments, the 5-year recurrence rate of such individuals remains as high as 50%-70%, with 10%-15% of patients progressing to muscle invasive BCa (MIBC). Every year, about 25% of new cases of BCa are categorized as MIBC, with a 5-year tumor-specific survival rate of around 50%. Platinum chemotherapy stands as the primary treatment strategy for advanced BCa; however, more than 50% of patients do not respond favorably to chemotherapy [6]. Following chemotherapy, advanced BCa patients who continue to progress or relapse typically face an average survival period of merely 12–14 months [7]. In addition, BCa has a high rate of progression and recurrence [8], underscoring the urgent need to elucidate the molecular mechanisms governing the malignant progression of BCa.

Various SOX (Sex determining region Y-box) proteins have been acknowledged for their role in disease development. More recently, they have gained recognition as important contributors to the onset and progression of human cancer [9]. Sex determining region Y-box 7 (SOX7) has been recognized as a member of the SOX family. Initially cloned by Shiozawa M in African Xenopus and mice [10, 11], it exerts a crucial function in development and endodermal differentiation. The human SOX7 gene is located on chromosome 8p23.1 and encodes 377 amino acids [12], comprising an intron and two exons. SOX7, SOX17, and SOX18 belong to the SOXF subfamily [13, 14]. Like other SOX factors, the N-terminus of SOXF factors contains highly conserved HMG DNA domains and transcriptional activation domains, which can bind to specific fragments on DNA, leading to alterations in DNA structure and transcription status [15–17]; C-end has β-catenin protein binding domain, thereby inhibiting the activity of the β-catenin/TCF transcription complex suppresses the Wnt signaling pathway [18]. SOX7 is involved in diverse developmental processes, including hematopoiesis [19], cardiovascular generation [15], endodermal differentiation [20], and myosatellite cell generation [21]. Recent studies have unveiled its tumor-suppressive role as a transcription factor in renal cell carcinoma [22], lung [23], breast [24], and prostate cancers [25]. However, the mechanism underlying the involvement of SOX7 in BCa remains unelucidated.

Globins, a protein family acknowledged for its ability to bind and transport oxygen, include prominent members like hemoglobin and myoglobin, which are pivotal for tissue oxygenation [26]. Recently, cytoglobin (CYGB), a focal point of current research within this family, has been identified as playing a unique role under certain conditions [27]. Globins are predominantly found in erythrocytes, muscle, and neural tissues; however, the expression of CYGB has been observed across diverse tissues and organs. Like other globins, CYGB possesses a compact helical structure and has the ability to reversibly bind diatomic gas molecules like NO, O2, and CO [28]. CYGB expression is known to increase in response to hypoxia, oxidative stress, or fibrotic stimuli. This is thought to be a mechanism that protects cells from oxidative stress damage. The abnormal expression of CYGB has been linked to various human cancers [29–31]. Investigations on Cygb−/− mice have revealed a significant increase in tumors in various organs such as lungs and livers in comparison to control mice, which confirms the tumor-suppressor effect of Cygb in vivo. Furthermore, the promotor region of CYGB exhibits heightened methylation in various types of cancer cell lines and solid tumors in comparison to normal tissues and non-tumor cell lines [32]. Nonetheless, the mechanism of how the CYGB promoter is methylated and how it exerts its tumor-suppressive effects has not been fully elucidated.

DNA methylation is a crucial type of epigenetic regulation involving the transfer of reactive methyl groups to a target site without altering the composition of the DNA sequence catalyzed by methyltransferases [33]. Both the dynamics and the biological outcome of methylation are the result of the activity of a complex protein mechanism consisting of writers, erasers, and readers. DNA methylation is present in numerous tissues, and DNA that is methylated has its transcriptional processes repressed, decreasing mRNA and protein expression of genes. Many diseases are associated with DNA methylation, including the development of human cancers and their progression [34].

In this research, we reveal that SOX7 is diminished in BCa, suggesting its potential role as a tumor suppressor in this context. As an oncogene, it can lead to altered CYGB expression by downregulating the transcription of DNMT3B (DNA methyltransferase 3 beta) and the methylation level of CYGB promoter. Ultimately, the altered CYGB can affect the progression of BCa. Our findings propose SOX7 and CYGB as potential biomarkers for BCa, with implications for both its prevention and diagnosis.

Result

SOX7 exhibits low expression in BCa and significantly correlates with poor prognosis

The RNA-sequencing of GSE32894 and 111 BCa tissues from Shanghai Tenth People's Hospital (STPH) revealed that SOX7 exhibited a significantly reduced expression in BCa tissues at T2-T4 stage than Tis-T1 (Fig. 1a and b). Tissue microarray staining of BCa tissues acquired from 61 patients was performed, and the results were scored based on histochemical staining (Fig. 1c). The acquired data revealed that individuals exhibiting low expression of SOX7 displayed heightened pathological and T-stage (T2-T4) BCa grades (Table. S1). Afterward, the acquired results underwent validation via the the Cancer Genome Atlas (TCGA) database, SOX7 exhibited significantly diminished expression in BCa tissues (Fig. 1d). Analysis of the GSE32894 dataset, often utilized to examine the link between genes and the prognosis of individuals with BCa, depicted a considerable association between the differential expression of SOX7 and the pelvic lymph nodes metastasis and distant metastasis in patients (Fig. 1e and f). Meanwhile, the acquired data implied that patients with higher SOX7 expression exhibited better prognosis (Fig. 1g).Fig. 1 SOX7 exhibits low expression in BCa and significantly correlates with poor prognosis. a Flowchart of RNA-sequencing of BCa. b Relationship between RNA expression of SOX7 in Tis-T1 and T2-T4 at STPH and GSE32894. c Tissue microarray staining of BCa tissues from 61 patients about NMIBC and MIBC. Scale bars: 5 mm/400 μm. d RNA expression of SOX7 in normal and BCa tissues in the TCGA database. Statistical test: Paired t test. e and f RNA expression of SOX7 and its correlation with pelvic lymph nodes metastasis (e) and distant metastasis (f) in GSE32894. Statistical test: Unpaired t test. g The link between RNA expression of SOX7 and overall survival of individuals in the GSE32894. Statistical test: Log-rank

The differential expression of SOX7 in BCa and normal bladder tissues was further verified through Immunohistochemistry (IHC), immunofluorescence (IF), and western blot (WB). The resulting data revealed that the expression of SOX7 in BCa tissues was markedly lower in comparison to the normal bladder tissues (Fig. 2a-c). Afterward, the variance in SOX7 transcription level in 25 pairs of cancer and normal bladder tissues was verified by qPCR, revealing a statistically significant reduction in the expression of SOX7 mRNA in BCa tissues (Fig. 2d). Additionally, this result remained consistent across the BCa cell lines. SV-HUC-1, a normal urothelial cell, served as the control. Compared with SV-HUC-1, RNA and protein levels of SOX7 in BCa cell lines T24, 5637, UMUC3, J82, and RT4 were significantly downregulated (Fig. 2e and f). The above findings are indicative of the consistent downregulation of SOX7 in BCa tissues and cells.Fig. 2 Downregulation of SOX7 in BCa tissues and cell lines. a IHC images illustrating the SOX7 expression in the tumor and normal tissues (n = 2). Scale bars: 50 μm. b IF images detailing the SOX7 expression in the tumor and normal tissues (n = 2). Scale bars: 50 μm. c WB describing the SOX7 expression in tumor and normal tissues (n = 5). d qPCR detection showing relative SOX7 mRNA expression in tumor and normal tissues (n = 25). Statistical test: Paired t test. e qPCR analysis displaying relative SOX7 mRNA expression in BCa cell lines. Statistical test: Unpaired t test. f WB illustrating SOX7 expression in BCa cell lines

SOX7 suppresses the proliferation of BCa cells in vivo

To verify whether SOX7 inhibits tumor progression, SOX7 was knocked down in the T24 and UMUC3 cell lines, while SOX7 overexpression was induced in UMUC3 and T24 cell lines. The efficiency of knockdown and overexpression of SOX7 was verified at the RNA and protein levels, respectively (Fig. 3a-c, Fig. S1a-b and S2a-b). After verification, the UMUC3 cell line with overexpression and T24 cell line with stable knockdown of SOX7 (SOX7-sh) were selected for subcutaneous tumor-bearing in BALB/c-nude mice (Fig. 3d). The volume of the tumor was assessed once a week and recorded (Fig. 3e). The mice were euthanized at the fourth week, and the tumor was dissected and weighed (Fig. 3f and g). The resulting data implied that heightened expression of SOX7 was capable of considerably suppressing the proliferation of BCa in vivo. Ki67 and PCNA are markers for cell proliferation. IHC performed on the extracted tumor body highlighted that the expression of Ki67 and PCNA was negatively linked to SOX7 expression (Fig. 3h and i). The above findings underscore the inhibitory role of SOX7 in the proliferation of BCa cells in vivo.Fig. 3 SOX7 suppresses the proliferation of BCa cells in vivo. a-b qPCR and WB detection of knockdown and SOX7 overexpression in T24 and UMUC3 cell lines. Statistical test: Unpaired t test. c IF images depicting knockdown and SOX7 overexpression in T23 and UMUC3 cell lines. Scale bars: 50 μm. d Schematic diagram of UMUC3 overexpressed with SOX7 injected into nude mice, with tumor size measured weekly. After four weeks, nude mice were euthanized, and tumors were removed, weighed, and recorded. Tissues were subsequently sectioned and stained. e Tumor volume growth curve. Statistical test: Unpaired t test. f Image of tumor size in nude mice after subcutaneous tumor loading for four weeks. g Histogram depicting the difference in tumor weight between the NC group and SOX7-sh group, Vector group and SOX7-oe group. Statistical test: Unpaired t test. h-i Representative IHC images of Ki67 and PCNA in mouse subcutaneous tumor tissue. Scale bars: 100 μm

Overexpression and knockdown of SOX7 affect the malignant progression of BCa in vitro

Subsequently, a series of experiments were executed to verify the impact of SOX7 on the progression of BCa. The cell counting kit-8 (CCK-8) experiment affirmed that in comparison to the vector group, cell proliferation was markedly enhanced in the SOX7 knockdown group and diminished in the SOX7 overexpression group (Fig. 4a, Fig S1c and S2c). Similarly, EdU and colony formation experiments also showed the same outcomes (Fig. 4b-c and Fig S1d-e and S2d-e). In addition, wound healing and transwell experiments revealed that the capacity of cells to migrate and invade was enhanced after knocking down SOX7 but weakened after overexpressing SOX7 (Fig. 4d-e, Fig S1f-g and S2f-g). The acquired data are indicative of the inhibitory potential of SOX7 in the malignant progression of BCa in vitro.Fig. 4 Overexpression and knockdown of SOX7 affect the malignant progression of BCa in vitro. a-b CCK-8 and colony formation assays were utilized to assess the viability of the T24 and UMUC3 cell lines following SOX7 knockdown and overexpression, respectively. Statistical test: Unpaired t test. c EdU assay was employed to estimate the viability of the T24 and UMUC3 cell lines following SOX7 knockdown and overexpression, respectively. Scale bars: 100 μm. Statistical test: Unpaired t test. d Wound-healing assay was conducted to analyze the migration abilities of T24 and UMUC3 cell lines following SOX7 knockdown and overexpression, respectively. Scale bars: 400 μm. Statistical test: Unpaired t test. e Transwell assay was executed to analyze the migratory and invasive capabilities of T24 and UMUC3 cell lines following SOX7 knockdown and overexpression, respectively. Scale bars: 100 μm. Statistical test: Unpaired t test

SOX7 upregulates CYGB expression in BCa

To investigate the impact of SOX7 on the malignant progression of BCa, RNA-sequencing was conducted on the T24 cell line overexpressing SOX7. Heatmap analysis (Fig. 5a) and Volcano plots (Fig. 5b) demonstrated that the mRNA levels of CYGB were elevated alongside the increased expression of SOX7. Consistent with the findings, KEGG pathway enrichment analysis of the RNA-sequencing results revealed enrichment in pathways in cancer (Fig. 5c). Subsequently, qPCR experiments were conducted to validate these results, confirming a significant increase in CYGB mRNA levels in cells overexpressing SOX7 (Fig. 5d).Fig. 5 SOX7 upregulates CYGB expression in BCa. a Heatmaps depicting RNA-sequencing results. b Volcano plots illustrating RNA-sequencing data. c KEGG enrichment analysis of RNA-sequencing. d mRNA expression of CYGB in T24 and UMUC3 cell lines overexpressing SOX7. Statistical test: Unpaired t test. e RNA expression of CYGB in normal and BCa tissues and disease-free survival in the TCGA database. Statistical test: Wilcoxon and Log-rank. f Link between the RNA expression level of SOX7 and tumor grade and progression in GSE13507. Statistical test: Unpaired t test

Next, exploration of CYGB expression and disease-free survival in BCa was performed using the TCGA database, revealing significantly lower expression of CYGB in BCa tissues and poorer disease-free survival (Fig. 5e). Additionally, analysis of the GSE13507 dataset indicated a correlation between CYGB expression and tumor pathological grade as well as progression (Fig. 5f). These outcomes highlight a potential oncogenic role for CYGB in BCa.

SOX7 promotes the methylation level of the CYGB promoter by inhibiting the transcription of DNMT3B

In breast cancer, it has been noted that CYGB exhibits low expression attributed to promoter methylation [27], thereby exerting a tumor-suppressive function. Similarly, a parallel phenomenon was observed in BCa. Methylation analysis of CYGB in BCa using MethHC [35] and cBioPortal [36] revealed elevated methylation levels in BCa tissues(Fig. 6a and b). Correlation analysis from the GEPIA database unveiled a negative correlation between DNMT3B and CYGB (Fig. 6c), while SOX7 and the methyltransferase DNMT3B also displayed a significant negative correlation (Fig. 6d). Given that SOX7, as a transcription factor, can bind to the DNA promoter and may inhibit the transcription of DNMT3B. Subsequent qPCR and WB analysis confirmed reduced DNMT3B expression upon overexpression of SOX7 in the T24 cell line (Fig. 6e and f).Fig. 6 SOX7 promotes the methylation level of the CYGB promoter by inhibiting the transcription of DNMT3B. a Analysis results of MethHC on methylation of CYGB. Statistical test: Unpaired t test. b Analysis results of cBioPortal on methylation of CYGB. Statistical test: Spearman. c Analysis results of GEPIA on the correlation between SOX7 and DNMT3B. Statistical test: Spearman. d Analysis results of GEPIA on the correlation between DNMT3B and CYGB. Statistical test: Spearman. e Reduced expression of DNMT3B mRNA after overexpression of SOX7 in T24 cell line. Statistical test: Unpaired t test. f Reduced expression of DNMT3B protein after overexpression of SOX7 in T24 cell line. g Motif schematic of SOX7. h Schematic diagram of the combination of SOX7 and 5’-AACAAT-3’ in the promoter region of DNMT3B. i Schematic diagram of SOX7 binding to DNA: In red is the amino acid sequence 45–113 of SOX7 and the region where it binds DNA. j ChIP-qPCR result of DNA fragments pulled-down by SOX7 antibody in T24 cell line. Non-specific IgG was used as control. Primer was specific to region of DNMT3B promoter. Statistical test: Unpaired t test. k ChIP-qPCR result of DNA fragments pulled-down by DNMT3B antibody in SOX7 overexpression T24 cell line. Non-specific IgG was used as control. Primer was specific to region of CYGB promoter. Statistical test: Unpaired t test. l Mechanism of SOX7 affecting the malignant progression of BCa

Motif of SOX7 was examined (Fig. 6g) [37] and protein-DNA binding prediction for the promoter region of SOX7 and DNMT3B was conducted. Remarkably, a binding site for SOX7 with the sequence 5’-AACAAT-3’ was identified in the DNMT3B promoter (Fig. 6h and i), consistent with the binding site for SOX7 reported in the UniPort database. Based on this, we designed primers targeting the DNMT3B promoter region, and the primers covered the binding site of SOX7. The results of chromatin immunoprecipitation-qPCR (ChIP-qPCR) showed that the primers were significantly enriched in SOX7 (Fig. 6j). To verify that DNMT3B regulates the methylation level of CYGB under the regulation of SOX7, ChIP- qPCR experiment on the CYGB promoter using DNMT3B antibody was conducted. The results showed that DNMT3B can bind to the CYGB promoter and its expression is regulated by SOX7 (Fig. 6k). Therefore, SOX7 reduces the methylation level of CYGB promoter by downregulating the expression of DNMT3B, thereby playing a role in inhibiting the progression of BCa following the increase in expression of CYGB (Fig. 6l).

Construction a nomogram to prediction of prognosis of BCa patients

To further evaluate the clinical value of SOX7 and CYGB using the TCGA database, a combined SOX7 + CYGB scoring system was constructed through bioinformatics analysis. This combined score was identified as an independent risk factor for assessing overall survival (OS) in individuals with BCa (Fig. 7a and b). Subsequently, a nomogram was developed based on this scoring system to predict patient prognosis (Fig. 7c), and the final calibration curve demonstrated that the nomogram exhibited good predictive performance for patient OS at 1, 2, and 3 years (Fig. 7d-f).Fig. 7 Combined prediction of prognosis of BCa patients by SOX7 and CYGB. a Univariate Cox regression analyses were carried out to examine variables significantly linked to OS. b Multivariate Cox regression analyses were carried out to examine variables significantly linked to OS. c Nomogram to predict the probability of 1-, 2-, and 3-year OS. d-f Calibration plots of the nomogram were utilized to predict the probability of 1-, 2-, and 3-year OS. Statistical test: Cox regression analysis (univariate and multivariate) is used to predict the predictive impact of genes and clinical pathological features

Discussion

In this investigation, a considerable downregulation of SOX7 expression in BCa was observed. Identified as a transcription factor, SOX7 demonstrates inhibitory effects on DNMT3B expression, thereby facilitating the upregulation of the downstream gene CYGB. Through this pathway, SOX7 plays a role in suppressing the migratory, proliferative, and invasive capacities of BCa. The findings of this research clarify the mechanism of BCa progression and can predict the prognosis of individuals with BCa through the combined score of SOX7 and CYGB.

More and more evidences suggest that SOX7 has a tumor-suppressive effect [38–40]. In many human cancers, the expression of SOX7 is downregulated [41], and overexpression of SOX7 protein can inhibit the proliferation of colon and lung cancer cells [23, 42]. When delving into the mechanism underlying SOX7 and cancer development, ample evidence supports the direct interaction of SOX7 with β-catenin binding, thereby facilitating the negative regulation of Wnt/β-catenin signaling to modulate cell growth [12, 43–45]. Nevertheless, subsequent studies have demonstrated that the anti-tumor impact of SOX7 is not solely reliant on the regulation of Wnt signaling. It has been observed that overexpression of mutant SOX7, lacking the β-catenin binding domain, still maintains the inhibitory activity of SOX7, as confirmed in colony formation experiments [28]. This indicates that in addition to Wnt signaling transduction, SOX7 may also participate in the regulation of other pathways related to cancer. Through experiments, it was found that SOX7 exhibits reduced expression in BCa. Following the overexpression of SOX7, a considerable decrease was observed in the proliferative, invasive, and migratory abilities of BCa. This underscores the tumor-suppressive role of SOX7 in BCa, akin to its effects observed in breast and lung cancers. Additionally, RNA-sequencing revealed a concurrent upregulation of CYGB upon SOX7 overexpression.

Recent evidence indicates that CYGB exerts anti-tumor effects across various human malignancies. In this research, a reduction in CYGB expression was observed in tissues linked to ovarian, head and neck, liver, lung, esophageal, pancreatic, and breast cancers [29, 46–53]. Moreover, diminished CYGB expression has been reported in individuals with glioma, correlating with elevated histological grade and increased risk of tumor recurrence [54]. Moreover, investigations into liver cancer development revealed that CYGB-deficient mice exhibited heightened susceptibility to liver cancer induction in both chemically induced and high-fat diet models [31]. In breast cancer, researchers have identified that CYGB inhibits breast cancer by suppressing glucose metabolism [27]. Collectively, these discoveries imply an anti-tumor role for CYGB. Similarly, this investigation highlights that CYGB, as a downstream target of SOX7, plays a tumor inhibitory function with SOX7 in BCa. Analyzing from a transcriptomic perspective, CYGB is highly expressed in normal BCa tissues, and patients with BCa with lower CYGB expression have poorer DFS, higher tumor grade, and poorer progression. Furthermore, our results indicate hypermethylation of the CYGB promoter in BCa.

Many genes are hypermethylated at their promoters and function as oncogenes in tumor development. However, upon restoration of their expression, these genes have the potential to impede tumor initiation and/or progression [55, 56]. Consistent with evidence demonstrating the epigenetic downregulation of CYGB in several types of malignancies [57, 58], demethylation treatment effectively restored CYGB expression. This confirms that promoter methylation participates in the suppression of CYGB expression across a wide range of cancers. As a member of the DNA methyltransferase family, DNMT3B has the capacity to enhance methylation on promoters of downstream genes, consequently suppressing the transcription process of these genes. Our analysis of the BCa database suggests a potential association between hypermethylation of CYGB and the malignant progression of BCa, which could be attributed to DNMT3B-mediated hypermethylation of the CYGB promoter.

In summary, this study reveals that SOX7 expression is downregulated in BCa, which has implications for regulating CYGB expression through the suppression of DNMT3B transcription and subsequent reduction methylation levels of CYGB promoter. Consequently, these changes in CYGB expression may impact the progression of BCa. Additionally, the combined SOX7 + CYGB score holds promise for predicting the prognosis of BCa patients. Nonetheless, our study has certain limitations. Such as a larger sample size is required to supplement and validate the developed scoring system. Although we analyzed the role of CYGB in BCa from the perspective of transcriptomics, the functional correlation of CYGB in BCa cells still needs further experimental verification.

Materials and methods

Data retrieval

The BCa dataset accessed at STPH, as well as The Cancer Genome Atlas (cancergenome.nih.gov/TCGA) and Gene Expression General Database (https://www.ncbi.nlm.nih.gov/GEO/), were used in this study. The STPH dataset comprised 111 patients who fulfilled the following three inclusion criteria: (1) The histological diagnosis confirmed by pathologists was BCa; (2) The prognostic information was complete and could be used for analytical studies; (3) RNA expression data met the available criteria. From November 2019 to April 2023, 111 patients with BCa were recruited from STPH data. All studies, including human participants, underwent review and approval by the STPH Ethics Committee (approval number 2021KN108). Since this investigation is retrospective in nature, informed consent of patients was not essential. Total RNA extraction, bilateral library generation, and RNA sequencing protocols for STPH were the same as established within our previously published studies [59].

Cell lines and cell culture

The human BCa cell lines utilized in this study included T24, 5637, UMUC3, J82, and RT4. Immortalized human normal bladder epithelial cell line SV-HUC-1 from the Chinese Academy of Sciences (Shanghai, China) was also used. SV-HUC-1 was cultured in F12k medium (Sigma-Aldrich, St. Louis, MO, USA), whereas RPMI-1640 medium (Thermo Fisher Science, Inc. USA) was utilized to culture T24, 5637, and UMUC3 cell lines. RT4 cell lines were grown in McCoy's 5A medium (Thermo Fisher Science, Inc. USA). The culture media utilized in these processes was enriched with 1% penicillin/streptomycin (Hyclone, Logan, UT, USA) and 10% fetal calf serum (FCS; Thermo Fisher Scientific, Inc. USA). The culturing of these cell lines was implemented in an incubator under 5% CO2 at a constant temperature of 37° C.

Western blot

Cells or tissues underwent lysis with RIPA (Beyotime, China) for 30 min, with subsequent centrifugation at 12,000 rpm for 10 min. The supernatant was taken, and total protein concentration was analyzed using Bicinchoninic Acid protein assay (BCA; Beyotime, China). Prior to Western blot analysis, protein samples (40 μg each) were subjected to electrophoresis utilizing 10% sodium dodecyl Sulfate–polyacrylamide gel (SDS-PAGE) with subsequent transfer to nitrocellulose membranes (Sigma-Aldrich; Merck KGaA). For blocking these membranes, a solution comprising PBS + 5% skim milk was utilized for 1 h, with subsequent overnight exposure to primary antibodies at 4° C. The membranes then underwent washing with PBST (thrice) and were exposed to horseradish peroxidase-coupled secondary antibodies at room temperature for 1 h. Detection of labeled proteins and determination of their concentration was carried out via a chemiluminescence imaging system (Tanon 5200 system, Tanon, Shanghai, China). Antibodies were shown at Supplementary information.

RNA isolation and quantitative real-time PCR (qPCR)

The extraction of total RNA from human tissues or cells was carried out utilizing the TRIzol reagent. The first cDNA strand was then generated utilizing a reverse transcription system kit (Vazyme Biotech Co., Ltd., China). Afterward, qPCR was performed utilizing the ABI Prism 7500 sequence detection system (Applied Biosystems, CA) and ChamQ Universal SYBR qPCR Master Mix (Vazyme, China).

The following parameters were established for qPCR: 95° C 5 min, 40 cycles (95° C 10 s, 60° C 30 s). The endogenous control was established as GAPDH. The relative folding changes were then determined in triplicate using the 2-ΔΔCT technique. Primers were shown at Supplementary information.

Lentivirus infection and cell transfection

Stable overexpression and knockdown of SOX7 in BCa cells was achieved through lentivirus infection. Initially, the plasmid was cloned into a lentiviral vector with subsequent transfection into HEK-293 T cells. The BCa cells underwent infection with the virus, and stable cell lines were established through selection with puromycin.

The cell transfection process was carried out utilizing 1 × 105 BCa cells/well inoculated into a 6-well plate. Then, siRNA or plasmid was transfected into cells utilizing lipofectamine 3000 (Invitrogen, USA) following the provided instructions. The efficiency of the transfection process was assessed via qPCR and Western blot techniques.

Nude mice xenograft assay

This assay involved subcutaneous transplantation of the tumor. Specifically, subcutaneous administration of 5 × 107 BCa cells was carried out into the right axilla of 4-week-old BALB/c-nude mice. After a 4-week period, the animals were euthanized, and the excised tumor was utilized for further follow-up experiments. Throughout the experimental period, the volume and diameter of the tumor, as well as the body weight of the mice, were monitored after a 7-day interval. Tumor volume was determined via the given formula: Tumor volume = π/6 × length × width2 [59].

Immunohistochemistry and immunofluorescence analysis

Fresh tissue was removed and underwent fixation in cold 4% paraformaldehyde (PFA), with subsequent embedding in paraffin blocks 48 h later. Following this process, the tissues were then sectioned. The sections were then de-paraffinized, graded alcohol dehydrated, antigen repaired, and blocked with 3% BSA in xylene for 1 h. After the completion of this process, incubation of these tissue sections was carried out with antibodies at 4° C for 18 h. The sections then underwent exposure to biotinylated rabbit anti-goat IgG for 20 min at 25° C, with subsequent exposure to streptavidin–horseradish peroxidase or fluorescent peroxidase for 30 min. Diaminobenzidine-H2O2 and hematoxylin were then used for tissue staining (this step is not required for IF).

The tissues were selected and grouped into MIBC and NMIBC for microarray. We regularly conduct clinical follow-up on patients and results have been published [60]. Followed by statistical analysis for group comparison using Pearson's chi-squared test.

Cell counting kit-8 assay

The degree of cell proliferation was detected utilizing the CCK-8 assay (Yeasen, Shanghai, China). For this purpose, 96-well plates were utilized for seeding 500 cells/wells. The next step involved introducing 10 μl of CCK-8 reagent into each well on days 0, 1, 2, 3, 4, and 5 of inoculation, respectively, with subsequent incubation of all samples at 37° C for 1.5 h. A microplate spectrophotometer (BioTek instrument, Winooski, VT) was utilized to examine the absorbance of the samples at 450 nm in order to determine cell proliferation.

Colony formation assay

For colony formation assays, 1 × 103 cells/wells were seeded into a 6-well plate and left for nine days. The 6-well plate was then rinsed with PBS (thrice), followed by fixation via 4% PFA. The cell community was then stained with 0.1% crystal violet (Yeasen, Shanghai, China) for 20 min, rinsed with PBS, and dried for 30 min. Finally, the cell colony was photographed with a digital camera.

EdU incorporation assay

The cells were seeded into 96-well plates at a density of 1 × 103 cells/well and evaluated for cell proliferation after 24 h using the EdU assay kit (K1077, APExBIO Corporation, US). The stained cells were then imaged using fluorescence microscopy (Nikon, Tokyo, Japan).

Wound healing assay

For the wound-healing assay, cells were cultured in 6-well plates, and cell layers were scratched by a sterile plastic suction pipette. Subsequent imaging was conducted using electron microscopy (Nikon, Tokyo, Japan) at 0, 12, and 24 h, respectively. The migratory capacity of cells was evaluated by assessing changes in the size of the injured area.

Migration and invasion assays

Transwell film with or without Matrix (Yeasen, Shanghai, China) was utilized for evaluating the invasive or migratory ability of BCa cells. Specifically, the chamber was placed into a 24-well plate, and 2.5 × 104 cells were inoculated into the upper chamber, containing 200 μl FCS-deficient medium. At the same time, 600 μl of culture medium with 10% FCS was introduced to the lower chamber. Following incubating at 37 ℃ for 24 h, the culture chamber was rinsed with PBS and fixed in 4% paraformaldehyde for about 15 min. The next step involved scraping off the cells on the upper side of the membrane with a cotton swab and staining via crystal violet at room temperature for about 15 min. Subsequently, the film was rinsed with PBS and photographed after drying out.

RNA-sequencing analysis

In this experiment, cells overexpressing SOX7 and negative control cells were extracted with subsequent isolation of the total RNA. An Illumina HiSeq 2000 (Shanghai OE Biotechnology Co., Ltd.) was utilized to conduct the RNA-seq procedure. The STAR aligner v2.5 was then utilized to comparatively assess the retrieved data with the UCSC human genome (hg19), and the relevant hits underwent quantification using feature counting software. The RNA-seq data were then subjected to analysis via "DESeq2" and "ClusterProfiler" software. The cut-off values were established as a multiplicity of change of 2 and an adjusted p-value of 0.01.

Protein-DNA molecular docking

The AlphaFold-predicted protein 3D structure file SOX7_HUMA (UniProt ID: Q9BT81) was downloaded from UniProt, and the DNA (5'-AACAAT-3') 3D structure was constructed using Discovery Studio 2019. HDOCK is designed to predict binding complexes between two molecules, including proteins and nucleic acids, by employing a hybrid docking strategy. Protein-DNA molecular docking was applied by setting appropriate docking parameters using protein as receptor and DNA as ligand.

Nomogram establishment and validation

Cox regression analyses (univariate and multivariate) were performed to predict the predictive impact of SOX7 + CYGB and other clinicopathological attributes. Factors that exhibited significance in multivariate Cox regression analyses were subsequently used to create nomograms, which underwent validation via the "Rms" R package v5.1 (https://cran.r-project.org/web/packages/rms/index.html). Validation was conducted by means of the consistency index (C-index) and the calibration graphs.

Chromatin immunoprecipitation assay

ChIP assays were conducted using the SimpleChIP Enzymatic Chromatin IP Kit (Cell Signaling Technology, MA, USA). The cells were crosslinked with 4% polyformaldehyde (Biosharp, Anhui, China) for 10 min, followed by treatment with glycine. After washing with cold PBS, the cells were lysed and the DNA was digested with micrococcal nuclease into 150–900 bp fragments. These fragments were immunoprecipitated overnight at 4 °C using antibodies against the specific protein or IgG control. The ChIP-enriched DNAs were then analyzed by qRT-PCR, and the data were normalized to the input.

Statistical analysis

All statistical analyses were conducted using SPSS software, and the results were plotted using GraphPad Prism 8 software. Shapiro–Wilk test was used to verify the normal distribution. Pairs of groups were compared using Student's t-test. Pearson's correlation analysis was employed to calculate the correlation between continuous variables. χ2 tests were used to analyze count data, and Fisher's exact test was applied when the sample size was less than 40. Cox regression analysis (univariate and multivariate) is used to predict the predictive impact of genes and clinical pathological features, in order to reduce the interference and bias of confounding variables on the analysis results. Error bar means standard deviation (SD).

Supplementary Information

Supplementary Material 1.

Acknowledgements

The authors would like to thank Dr. Shuwen Sun, Dr. Hong Wang, Dr. Bing Xiong and Dr. Yajuan Hao for excellent technical support and Bullet Edits Limited for the linguistic editing of the manuscript. The images (Figs. 1a, 3d and 6l) in this article are created on BioRender (https://www.biorender.com/).

Authors’ contributions

ZJC, ZWT wrote the manuscript, LJ generated the figures. LYC, JYF, AA, CHY GZR and CHT contributed to the clinical and experimental data, MSY, GYD and XTY analyzed the data, XYF made significant contributions in the process of revising articles, WY, YXD and YY contributed to editing the manuscript. All authors read and approved the final manuscript.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by Experimental Animal Fund of Shanghai Science and Technology Commission (No. 201409004000, No.22140903800).

Availability of data and materials

All data are available from the corresponding authors upon reasonable request.

Declarations

Ethics approval and consent to participate

The present study was approved by the ethical committee of Shanghai Tenth People’s Hospital of Tongji University. Approval number for animal experiment in ethic approval: SHDSYY-2024–1533. Approval number for clinical experiment in ethic approval: SHSY-IEC-4.1/19–120/01.

Consent for publication

Informed consent was obtained from all participants included in the study.

Competing interests

The authors have declared that no competing interest exists.

Publisher’s Note

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

Jingcheng Zhang, Wentao Zhang and Ji Liu contributed equally to this work and share the first authorship.
==== Refs
References

1. Sung H Ferlay J Siegel RL Laversanne M Soerjomataram I Jemal A Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries CA Cancer J Clin 2021 71 3 209 249 10.3322/caac.21660 33538338
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49. 10.3322/caac.21660.33538338 10.3322/caac.21660
2. Dyrskjøt L Hansel DE Efstathiou JA Knowles MA Galsky MD Teoh J Bladder cancer Nat Rev Dis Primers 2023 9 1 58 10.1038/s41572-023-00468-9 37884563
Dyrskjøt L, Hansel DE, Efstathiou JA, Knowles MA, Galsky MD, Teoh J, et al. Bladder cancer Nat Rev Dis Primers. 2023;9(1):58. 10.1038/s41572-023-00468-9.37884563 10.1038/s41572-023-00468-9
3. Lenis AT Lec PM Chamie K Mshs MD Bladder Cancer: A Review JAMA 2020 324 19 1980 1991 10.1001/jama.2020.17598 33201207
Lenis AT, Lec PM, Chamie K, Mshs MD. Bladder Cancer: A Review. JAMA. 2020;324(19):1980–91. 10.1001/jama.2020.17598.33201207 10.1001/jama.2020.17598
4. Compérat E Amin MB Cathomas R Choudhury A De Santis M Kamat A Current best practice for bladder cancer: a narrative review of diagnostics and treatments Lancet 2022 400 10364 1712 1721 10.1016/s0140-6736(22)01188-6 36174585
Compérat E, Amin MB, Cathomas R, Choudhury A, De Santis M, Kamat A, et al. Current best practice for bladder cancer: a narrative review of diagnostics and treatments. Lancet. 2022;400(10364):1712–21. 10.1016/s0140-6736(22)01188-6.36174585 10.1016/s0140-6736(22)01188-6
5. Babjuk M Burger M Capoun O Cohen D Compérat EM Dominguez Escrig JL European Association of Urology Guidelines on Non-muscle-invasive Bladder Cancer (Ta, T1, and Carcinoma in Situ) Eur Urol 2022 81 1 75 94 10.1016/j.eururo.2021.08.010 34511303
Babjuk M, Burger M, Capoun O, Cohen D, Compérat EM, Dominguez Escrig JL, et al. European Association of Urology Guidelines on Non-muscle-invasive Bladder Cancer (Ta, T1, and Carcinoma in Situ). Eur Urol. 2022;81(1):75–94. 10.1016/j.eururo.2021.08.010.34511303 10.1016/j.eururo.2021.08.010
6. Patel VG Oh WK Galsky MD Treatment of muscle-invasive and advanced bladder cancer in 2020 CA Cancer J Clin 2020 70 5 404 423 10.3322/caac.21631 32767764
Patel VG, Oh WK, Galsky MD. Treatment of muscle-invasive and advanced bladder cancer in 2020. CA Cancer J Clin. 2020;70(5):404–23. 10.3322/caac.21631.32767764 10.3322/caac.21631
7. Seidl C Targets for Therapy of Bladder Cancer Semin Nucl Med 2020 50 2 162 170 10.1053/j.semnuclmed.2020.02.006 32172801
Seidl C. Targets for Therapy of Bladder Cancer. Semin Nucl Med. 2020;50(2):162–70. 10.1053/j.semnuclmed.2020.02.006.32172801 10.1053/j.semnuclmed.2020.02.006
8. Flaig TW Spiess PE Abern M Agarwal N Bangs R Buyyounouski MK Bladder Cancer, Version 3.2024 J Natl Compr Canc Netw. 2024 22 4 216 25 10.6004/jnccn.2024.0024 38754471
Flaig TW, Spiess PE, Abern M, Agarwal N, Bangs R, Buyyounouski MK, et al. Bladder Cancer, Version 3.2024. J Natl Compr Canc Netw. 2024;22(4):216–25. 10.6004/jnccn.2024.0024.38754471 10.6004/jnccn.2024.0024
9. Grimm D Bauer J Wise P Krüger M Simonsen U Wehland M The role of SOX family members in solid tumours and metastasis Semin Cancer Biol 2020 67 Pt 1 122 153 10.1016/j.semcancer.2019.03.004 30914279
Grimm D, Bauer J, Wise P, Krüger M, Simonsen U, Wehland M, et al. The role of SOX family members in solid tumours and metastasis. Semin Cancer Biol. 2020;67(Pt 1):122–53. 10.1016/j.semcancer.2019.03.004.30914279 10.1016/j.semcancer.2019.03.004
10. Shiozawa M Hiraoka Y Komatsu N Ogawa M Sakai Y Aiso S Cloning and characterization of Xenopus laevis xSox7 cDNA Biochim Biophys Acta 1996 1309 1–2 73 76 10.1016/s0167-4781(96)00145-5 8950180
Shiozawa M, Hiraoka Y, Komatsu N, Ogawa M, Sakai Y, Aiso S. Cloning and characterization of Xenopus laevis xSox7 cDNA. Biochim Biophys Acta. 1996;1309(1–2):73–6. 10.1016/s0167-4781(96)00145-5.8950180 10.1016/s0167-4781(96)00145-5
11. Taniguchi K Hiraoka Y Ogawa M Sakai Y Kido S Aiso S Isolation and characterization of a mouse SRY-related cDNA, mSox7 Biochim Biophys Acta 1999 1445 2 225 231 10.1016/s0167-4781(99)00047-0 10320775
Taniguchi K, Hiraoka Y, Ogawa M, Sakai Y, Kido S, Aiso S. Isolation and characterization of a mouse SRY-related cDNA, mSox7. Biochim Biophys Acta. 1999;1445(2):225–31. 10.1016/s0167-4781(99)00047-0.10320775 10.1016/s0167-4781(99)00047-0
12. Takash W Cañizares J Bonneaud N Poulat F Mattéi MG Jay P SOX7 transcription factor: sequence, chromosomal localisation, expression, transactivation and interference with Wnt signalling Nucleic Acids Res 2001 29 21 4274 4283 10.1093/nar/29.21.4274 11691915
Takash W, Cañizares J, Bonneaud N, Poulat F, Mattéi MG, Jay P, et al. SOX7 transcription factor: sequence, chromosomal localisation, expression, transactivation and interference with Wnt signalling. Nucleic Acids Res. 2001;29(21):4274–83. 10.1093/nar/29.21.4274.11691915 10.1093/nar/29.21.4274
13. Bowles J Schepers G Koopman P Phylogeny of the SOX family of developmental transcription factors based on sequence and structural indicators Dev Biol 2000 227 2 239 255 10.1006/dbio.2000.9883 11071752
Bowles J, Schepers G, Koopman P. Phylogeny of the SOX family of developmental transcription factors based on sequence and structural indicators. Dev Biol. 2000;227(2):239–55. 10.1006/dbio.2000.9883.11071752 10.1006/dbio.2000.9883
14. Higashijima Y Kanki Y Molecular mechanistic insights: The emerging role of SOXF transcription factors in tumorigenesis and development Semin Cancer Biol 2020 67 Pt 1 39 48 10.1016/j.semcancer.2019.09.008 31536760
Higashijima Y, Kanki Y. Molecular mechanistic insights: The emerging role of SOXF transcription factors in tumorigenesis and development. Semin Cancer Biol. 2020;67(Pt 1):39–48. 10.1016/j.semcancer.2019.09.008.31536760 10.1016/j.semcancer.2019.09.008
15. Lilly AJ Lacaud G Kouskoff V SOXF transcription factors in cardiovascular development Semin Cell Dev Biol 2017 63 50 57 10.1016/j.semcdb.2016.07.021 27470491
Lilly AJ, Lacaud G, Kouskoff V. SOXF transcription factors in cardiovascular development. Semin Cell Dev Biol. 2017;63:50–7. 10.1016/j.semcdb.2016.07.021.27470491 10.1016/j.semcdb.2016.07.021
16. Pevny LH Lovell-Badge R Sox genes find their feet Curr Opin Genet Dev 1997 7 3 338 344 10.1016/s0959-437x(97)80147-5 9229109
Pevny LH, Lovell-Badge R. Sox genes find their feet. Curr Opin Genet Dev. 1997;7(3):338–44. 10.1016/s0959-437x(97)80147-5.9229109 10.1016/s0959-437x(97)80147-5
17. Reményi A Lins K Nissen LJ Reinbold R Schöler HR Wilmanns M Crystal structure of a POU/HMG/DNA ternary complex suggests differential assembly of Oct4 and Sox2 on two enhancers Genes Dev 2003 17 16 2048 2059 10.1101/gad.269303 12923055
Reményi A, Lins K, Nissen LJ, Reinbold R, Schöler HR, Wilmanns M. Crystal structure of a POU/HMG/DNA ternary complex suggests differential assembly of Oct4 and Sox2 on two enhancers. Genes Dev. 2003;17(16):2048–59. 10.1101/gad.269303.12923055 10.1101/gad.269303
18. Kormish JD Sinner D Zorn AM Interactions between SOX factors and Wnt/beta-catenin signaling in development and disease Dev Dyn 2010 239 1 56 68 10.1002/dvdy.22046 19655378
Kormish JD, Sinner D, Zorn AM. Interactions between SOX factors and Wnt/beta-catenin signaling in development and disease. Dev Dyn. 2010;239(1):56–68. 10.1002/dvdy.22046.19655378 10.1002/dvdy.22046
19. Lilly AJ Costa G Largeot A Fadlullah MZ Lie ALM Lacaud G Interplay between SOX7 and RUNX1 regulates hemogenic endothelial fate in the yolk sac Development 2016 143 23 4341 4351 10.1242/dev.140970 27802172
Lilly AJ, Costa G, Largeot A, Fadlullah MZ, Lie ALM, Lacaud G, et al. Interplay between SOX7 and RUNX1 regulates hemogenic endothelial fate in the yolk sac. Development. 2016;143(23):4341–51. 10.1242/dev.140970.27802172 10.1242/dev.140970
20. Butler AM, Owens DA, Wang L, King ML. A novel role for sox7 in Xenopus early primordial germ cell development: mining the PGC transcriptome. Development. 2018;145(1). 10.1242/dev.155978.
21. Rajgara RF Lala-Tabbert N Marchildon F Lamarche É MacDonald JK Scott DA SOX7 Is Required for Muscle Satellite Cell Development and Maintenance Stem Cell Reports 2017 9 4 1139 1151 10.1016/j.stemcr.2017.08.014 28943254
Rajgara RF, Lala-Tabbert N, Marchildon F, Lamarche É, MacDonald JK, Scott DA, et al. SOX7 Is Required for Muscle Satellite Cell Development and Maintenance. Stem Cell Reports. 2017;9(4):1139–51. 10.1016/j.stemcr.2017.08.014.28943254 10.1016/j.stemcr.2017.08.014
22. Wang L Fan Y Zhang L Li L Kuang G Luo C Classic SRY-box protein SOX7 functions as a tumor suppressor regulating WNT signaling and is methylated in renal cell carcinoma Faseb j 2019 33 1 254 263 10.1096/fj.201701453RR 29957056
Wang L, Fan Y, Zhang L, Li L, Kuang G, Luo C, et al. Classic SRY-box protein SOX7 functions as a tumor suppressor regulating WNT signaling and is methylated in renal cell carcinoma. Faseb j. 2019;33(1):254–63. 10.1096/fj.201701453RR.29957056 10.1096/fj.201701453RR
23. Hayano T Garg M Yin D Sudo M Kawamata N Shi S SOX7 is down-regulated in lung cancer J Exp Clin Cancer Res 2013 32 1 17 10.1186/1756-9966-32-17 23557216
Hayano T, Garg M, Yin D, Sudo M, Kawamata N, Shi S, et al. SOX7 is down-regulated in lung cancer. J Exp Clin Cancer Res. 2013;32(1):17. 10.1186/1756-9966-32-17.23557216 10.1186/1756-9966-32-17
24. Sun QY Ding LW Johnson K Zhou S Tyner JW Yang H SOX7 regulates MAPK/ERK-BIM mediated apoptosis in cancer cells Oncogene 2019 38 34 6196 6210 10.1038/s41388-019-0865-8 31332289
Sun QY, Ding LW, Johnson K, Zhou S, Tyner JW, Yang H, et al. SOX7 regulates MAPK/ERK-BIM mediated apoptosis in cancer cells. Oncogene. 2019;38(34):6196–210. 10.1038/s41388-019-0865-8.31332289 10.1038/s41388-019-0865-8
25. Peng W Guo L Tang R Liu X Jin R Dong JT Sox7 negatively regulates prostate-specific membrane antigen (PSMA) expression through PSMA-enhancer Prostate 2019 79 4 370 378 10.1002/pros.23743 30488457
Peng W, Guo L, Tang R, Liu X, Jin R, Dong JT, et al. Sox7 negatively regulates prostate-specific membrane antigen (PSMA) expression through PSMA-enhancer. Prostate. 2019;79(4):370–8. 10.1002/pros.23743.30488457 10.1002/pros.23743
26. Keller TCSt Lechauve C Keller AS Brooks S Weiss MJ Columbus L The role of globins in cardiovascular physiology Physiol Rev. 2022 102 2 859 92 10.1152/physrev.00037.2020 34486392
Keller TCSt, Lechauve C, Keller AS, Brooks S, Weiss MJ, Columbus L, et al. The role of globins in cardiovascular physiology. Physiol Rev. 2022;102(2):859–92. 10.1152/physrev.00037.2020.34486392 10.1152/physrev.00037.2020
27. Feng Y Wu M Li S He X Tang J Peng W The epigenetically downregulated factor CYGB suppresses breast cancer through inhibition of glucose metabolism J Exp Clin Cancer Res 2018 37 1 313 10.1186/s13046-018-0979-9 30545372
Feng Y, Wu M, Li S, He X, Tang J, Peng W, et al. The epigenetically downregulated factor CYGB suppresses breast cancer through inhibition of glucose metabolism. J Exp Clin Cancer Res. 2018;37(1):313. 10.1186/s13046-018-0979-9.30545372 10.1186/s13046-018-0979-9
28. Zweier JL, Hemann C, Kundu T, Ewees MG, Khaleel SA, Samouilov A, et al. Cytoglobin has potent superoxide dismutase function. Proc Natl Acad Sci U S A. 2021;118(52). 10.1073/pnas.2105053118.
29. Shaw RJ Omar MM Rokadiya S Kogera FA Lowe D Hall GL Cytoglobin is upregulated by tumour hypoxia and silenced by promoter hypermethylation in head and neck cancer Br J Cancer 2009 101 1 139 144 10.1038/sj.bjc.6605121 19568272
Shaw RJ, Omar MM, Rokadiya S, Kogera FA, Lowe D, Hall GL, et al. Cytoglobin is upregulated by tumour hypoxia and silenced by promoter hypermethylation in head and neck cancer. Br J Cancer. 2009;101(1):139–44. 10.1038/sj.bjc.6605121.19568272 10.1038/sj.bjc.6605121
30. Shivapurkar N Stastny V Okumura N Girard L Xie Y Prinsen C Cytoglobin, the newest member of the globin family, functions as a tumor suppressor gene Cancer Res 2008 68 18 7448 7456 10.1158/0008-5472.Can-08-0565 18794132
Shivapurkar N, Stastny V, Okumura N, Girard L, Xie Y, Prinsen C, et al. Cytoglobin, the newest member of the globin family, functions as a tumor suppressor gene. Cancer Res. 2008;68(18):7448–56. 10.1158/0008-5472.Can-08-0565.18794132 10.1158/0008-5472.Can-08-0565
31. le Thuy TT Morita T Yoshida K Wakasa K Iizuka M Ogawa T Promotion of liver and lung tumorigenesis in DEN-treated cytoglobin-deficient mice Am J Pathol 2011 179 2 1050 1060 10.1016/j.ajpath.2011.05.006 21684245
le Thuy TT, Morita T, Yoshida K, Wakasa K, Iizuka M, Ogawa T, et al. Promotion of liver and lung tumorigenesis in DEN-treated cytoglobin-deficient mice. Am J Pathol. 2011;179(2):1050–60. 10.1016/j.ajpath.2011.05.006.21684245 10.1016/j.ajpath.2011.05.006
32. Oleksiewicz U Liloglou T Tasopoulou KM Daskoulidou N Bryan J Gosney JR Cytoglobin has bimodal: tumour suppressor and oncogene functions in lung cancer cell lines Hum Mol Genet 2013 22 16 3207 3217 10.1093/hmg/ddt174 23591990
Oleksiewicz U, Liloglou T, Tasopoulou KM, Daskoulidou N, Bryan J, Gosney JR, et al. Cytoglobin has bimodal: tumour suppressor and oncogene functions in lung cancer cell lines. Hum Mol Genet. 2013;22(16):3207–17. 10.1093/hmg/ddt174.23591990 10.1093/hmg/ddt174
33. Dai X Ren T Zhang Y Nan N Methylation multiplicity and its clinical values in cancer Expert Rev Mol Med 2021 23 e2 10.1017/erm.2021.4 33787478
Dai X, Ren T, Zhang Y, Nan N. Methylation multiplicity and its clinical values in cancer. Expert Rev Mol Med. 2021;23: e2. 10.1017/erm.2021.4.33787478 10.1017/erm.2021.4
34. Robertson KD Jones PA DNA methylation: past, present and future directions Carcinogenesis 2000 21 3 461 467 10.1093/carcin/21.3.461 10688866
Robertson KD, Jones PA. DNA methylation: past, present and future directions. Carcinogenesis. 2000;21(3):461–7. 10.1093/carcin/21.3.461.10688866 10.1093/carcin/21.3.461
35. Huang HY Li J Tang Y Huang YX Chen YG Xie YY MethHC 2.0: information repository of DNA methylation and gene expression in human cancer Nucleic Acids Res. 2021 49 D1 D1268 d75 10.1093/nar/gkaa1104 33270889
Huang HY, Li J, Tang Y, Huang YX, Chen YG, Xie YY, et al. MethHC 2.0: information repository of DNA methylation and gene expression in human cancer. Nucleic Acids Res. 2021;49(D1):D1268-d75. 10.1093/nar/gkaa1104.33270889 10.1093/nar/gkaa1104
36. Cerami E Gao J Dogrusoz U Gross BE Sumer SO Aksoy BA The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data Cancer Discov 2012 2 5 401 404 10.1158/2159-8290.Cd-12-0095 22588877
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4. 10.1158/2159-8290.Cd-12-0095.22588877 10.1158/2159-8290.Cd-12-0095
37. Lambert SA Jolma A Campitelli LF Das PK Yin Y Albu M The Human Transcription Factors Cell 2018 172 4 650 665 10.1016/j.cell.2018.01.029 29425488
Lambert SA, Jolma A, Campitelli LF, Das PK, Yin Y, Albu M, et al. The Human Transcription Factors. Cell. 2018;172(4):650–65. 10.1016/j.cell.2018.01.029.29425488 10.1016/j.cell.2018.01.029
38. Stovall DB Wan M Miller LD Cao P Maglic D Zhang Q The regulation of SOX7 and its tumor suppressive role in breast cancer Am J Pathol 2013 183 5 1645 1653 10.1016/j.ajpath.2013.07.025 24012678
Stovall DB, Wan M, Miller LD, Cao P, Maglic D, Zhang Q, et al. The regulation of SOX7 and its tumor suppressive role in breast cancer. Am J Pathol. 2013;183(5):1645–53. 10.1016/j.ajpath.2013.07.025.24012678 10.1016/j.ajpath.2013.07.025
39. Oh KY Hong KO Huh YS Lee JI Hong SD reased expression of SOX7 induces cell proliferation and invasion and correlates with poor prognosis in oral squamous cell carcinoma J Oral Pathol Med 2017 46 9 752 758 10.1111/jop.12566 28266739
Oh KY, Hong KO, Huh YS, Lee JI, Hong SD. reased expression of SOX7 induces cell proliferation and invasion and correlates with poor prognosis in oral squamous cell carcinoma. J Oral Pathol Med. 2017;46(9):752–8. 10.1111/jop.12566.28266739 10.1111/jop.12566
40. Han L Wang W Ding W Zhang L MiR-9 is involved in TGF-β1-induced lung cancer cell invasion and adhesion by targeting SOX7 J Cell Mol Med 2017 21 9 2000 2008 10.1111/jcmm.13120 28266181
Han L, Wang W, Ding W, Zhang L. MiR-9 is involved in TGF-β1-induced lung cancer cell invasion and adhesion by targeting SOX7. J Cell Mol Med. 2017;21(9):2000–8. 10.1111/jcmm.13120.28266181 10.1111/jcmm.13120
41. Einolf HJ Lin W Won CS Wang L Gu H Chun DY Physiologically Based Pharmacokinetic Model Predictions of Panobinostat (LBH589) as a Victim and Perpetrator of Drug-Drug Interactions Drug Metab Dispos 2017 45 12 1304 1316 10.1124/dmd.117.076851 28912253
Einolf HJ, Lin W, Won CS, Wang L, Gu H, Chun DY, et al. Physiologically Based Pharmacokinetic Model Predictions of Panobinostat (LBH589) as a Victim and Perpetrator of Drug-Drug Interactions. Drug Metab Dispos. 2017;45(12):1304–16. 10.1124/dmd.117.076851.28912253 10.1124/dmd.117.076851
42. Zhang Y Huang S Dong W Li L Feng Y Pan L SOX7, down-regulated in colorectal cancer, induces apoptosis and inhibits proliferation of colorectal cancer cells Cancer Lett 2009 277 1 29 37 10.1016/j.canlet.2008.11.014 19108950
Zhang Y, Huang S, Dong W, Li L, Feng Y, Pan L, et al. SOX7, down-regulated in colorectal cancer, induces apoptosis and inhibits proliferation of colorectal cancer cells. Cancer Lett. 2009;277(1):29–37. 10.1016/j.canlet.2008.11.014.19108950 10.1016/j.canlet.2008.11.014
43. Stovall DB Cao P Sui G SOX7: from a developmental regulator to an emerging tumor suppressor Histol Histopathol. 2014 29 4 439 45 10.14670/hh-29.10.439 24288056
Stovall DB, Cao P, Sui G. SOX7: from a developmental regulator to an emerging tumor suppressor. Histol Histopathol. 2014;29(4):439–45. 10.14670/hh-29.10.439.24288056 10.14670/hh-29.10.439
44. Liu H Mastriani E Yan ZQ Yin SY Zeng Z Wang H SOX7 co-regulates Wnt/β-catenin signaling with Axin-2: both expressed at low levels in breast cancer Sci Rep 2016 6 26136 10.1038/srep26136 27188720
Liu H, Mastriani E, Yan ZQ, Yin SY, Zeng Z, Wang H, et al. SOX7 co-regulates Wnt/β-catenin signaling with Axin-2: both expressed at low levels in breast cancer. Sci Rep. 2016;6:26136. 10.1038/srep26136.27188720 10.1038/srep26136
45. Fan R He H Yao W Zhu Y Zhou X Gui M SOX7 Suppresses Wnt Signaling by Disrupting β-Catenin/BCL9 Interaction DNA Cell Biol 2018 37 2 126 132 10.1089/dna.2017.3866 29271667
Fan R, He H, Yao W, Zhu Y, Zhou X, Gui M, et al. SOX7 Suppresses Wnt Signaling by Disrupting β-Catenin/BCL9 Interaction. DNA Cell Biol. 2018;37(2):126–32. 10.1089/dna.2017.3866.29271667 10.1089/dna.2017.3866
46. McRonald FE Liloglou T Xinarianos G Hill L Rowbottom L Langan JE Down-regulation of the cytoglobin gene, located on 17q25, in tylosis with oesophageal cancer (TOC): evidence for trans-allele repression Hum Mol Genet 2006 15 8 1271 1277 10.1093/hmg/ddl042 16510494
McRonald FE, Liloglou T, Xinarianos G, Hill L, Rowbottom L, Langan JE, et al. Down-regulation of the cytoglobin gene, located on 17q25, in tylosis with oesophageal cancer (TOC): evidence for trans-allele repression. Hum Mol Genet. 2006;15(8):1271–7. 10.1093/hmg/ddl042.16510494 10.1093/hmg/ddl042
47. Xinarianos G McRonald FE Risk JM Bowers NL Nikolaidis G Field JK Frequent genetic and epigenetic abnormalities contribute to the deregulation of cytoglobin in non-small cell lung cancer Hum Mol Genet 2006 15 13 2038 2044 10.1093/hmg/ddl128 16698880
Xinarianos G, McRonald FE, Risk JM, Bowers NL, Nikolaidis G, Field JK, et al. Frequent genetic and epigenetic abnormalities contribute to the deregulation of cytoglobin in non-small cell lung cancer. Hum Mol Genet. 2006;15(13):2038–44. 10.1093/hmg/ddl128.16698880 10.1093/hmg/ddl128
48. Presneau N Dewar K Forgetta V Provencher D Mes-Masson AM Tonin PN Loss of heterozygosity and transcriptome analyses of a 1.2 Mb candidate ovarian cancer tumor suppressor locus region at 17q25.1-q25.2 Mol Carcinog. 2005 43 3 141 54 10.1002/mc.20096 15937959
Presneau N, Dewar K, Forgetta V, Provencher D, Mes-Masson AM, Tonin PN. Loss of heterozygosity and transcriptome analyses of a 1.2 Mb candidate ovarian cancer tumor suppressor locus region at 17q25.1-q25.2. Mol Carcinog. 2005;43(3):141–54. 10.1002/mc.20096.15937959 10.1002/mc.20096
49. Chua PJ Yip GW Bay BH Cell cycle arrest induced by hydrogen peroxide is associated with modulation of oxidative stress related genes in breast cancer cells Exp Biol Med (Maywood) 2009 234 9 1086 1094 10.3181/0903-rm-98 19596828
Chua PJ, Yip GW, Bay BH. Cell cycle arrest induced by hydrogen peroxide is associated with modulation of oxidative stress related genes in breast cancer cells. Exp Biol Med (Maywood). 2009;234(9):1086–94. 10.3181/0903-rm-98.19596828 10.3181/0903-rm-98
50. Wojnarowicz PM Provencher DM Mes-Masson AM Tonin PN Chromosome 17q25 genes, RHBDF2 and CYGB, in ovarian cancer Int J Oncol 2012 40 6 1865 1880 10.3892/ijo.2012.1371 22344671
Wojnarowicz PM, Provencher DM, Mes-Masson AM, Tonin PN. Chromosome 17q25 genes, RHBDF2 and CYGB, in ovarian cancer. Int J Oncol. 2012;40(6):1865–80. 10.3892/ijo.2012.1371.22344671 10.3892/ijo.2012.1371
51. Chen H Zhao X Meng T Expression and biological role of cytoglobin in human ovarian cancer Tumour Biol 2014 35 7 6933 6939 10.1007/s13277-014-1941-x 24737588
Chen H, Zhao X, Meng T. Expression and biological role of cytoglobin in human ovarian cancer. Tumour Biol. 2014;35(7):6933–9. 10.1007/s13277-014-1941-x.24737588 10.1007/s13277-014-1941-x
52. Hubers AJ Heideman DA Burgers SA Herder GJ Sterk PJ Rhodius RJ DNA hypermethylation analysis in sputum for the diagnosis of lung cancer: training validation set approach Br J Cancer 2015 112 6 1105 1113 10.1038/bjc.2014.636 25719833
Hubers AJ, Heideman DA, Burgers SA, Herder GJ, Sterk PJ, Rhodius RJ, et al. DNA hypermethylation analysis in sputum for the diagnosis of lung cancer: training validation set approach. Br J Cancer. 2015;112(6):1105–13. 10.1038/bjc.2014.636.25719833 10.1038/bjc.2014.636
53. Latina A Viticchiè G Lena AM Piro MC Annicchiarico-Petruzzelli M Melino G ΔNp63 targets cytoglobin to inhibit oxidative stress-induced apoptosis in keratinocytes and lung cancer Oncogene 2016 35 12 1493 1503 10.1038/onc.2015.222 26096935
Latina A, Viticchiè G, Lena AM, Piro MC, Annicchiarico-Petruzzelli M, Melino G, et al. ΔNp63 targets cytoglobin to inhibit oxidative stress-induced apoptosis in keratinocytes and lung cancer. Oncogene. 2016;35(12):1493–503. 10.1038/onc.2015.222.26096935 10.1038/onc.2015.222
54. Xu HW Huang YJ Xie ZY Lin L Guo YC Zhuang ZR The expression of cytoglobin as a prognostic factor in gliomas: a retrospective analysis of 88 patients BMC Cancer 2013 13 247 10.1186/1471-2407-13-247 23688241
Xu HW, Huang YJ, Xie ZY, Lin L, Guo YC, Zhuang ZR, et al. The expression of cytoglobin as a prognostic factor in gliomas: a retrospective analysis of 88 patients. BMC Cancer. 2013;13:247. 10.1186/1471-2407-13-247.23688241 10.1186/1471-2407-13-247
55. Qi YN Liu Z Hong LL Li P Ling ZQ Methyltransferase-like proteins in cancer biology and potential therapeutic targeting J Hematol Oncol 2023 16 1 89 10.1186/s13045-023-01477-7 37533128
Qi YN, Liu Z, Hong LL, Li P, Ling ZQ. Methyltransferase-like proteins in cancer biology and potential therapeutic targeting. J Hematol Oncol. 2023;16(1):89. 10.1186/s13045-023-01477-7.37533128 10.1186/s13045-023-01477-7
56. Park IG Jeon M Kim H Lee JM Coordinated methyl readers: Functional communications in cancer Semin Cancer Biol 2022 83 88 99 10.1016/j.semcancer.2021.03.015 33753223
Park IG, Jeon M, Kim H, Lee JM. Coordinated methyl readers: Functional communications in cancer. Semin Cancer Biol. 2022;83:88–99. 10.1016/j.semcancer.2021.03.015.33753223 10.1016/j.semcancer.2021.03.015
57. Daskalos A Oleksiewicz U Filia A Nikolaidis G Xinarianos G Gosney JR UHRF1-mediated tumor suppressor gene inactivation in nonsmall cell lung cancer Cancer 2011 117 5 1027 1037 10.1002/cncr.25531 21351083
Daskalos A, Oleksiewicz U, Filia A, Nikolaidis G, Xinarianos G, Gosney JR, et al. UHRF1-mediated tumor suppressor gene inactivation in nonsmall cell lung cancer. Cancer. 2011;117(5):1027–37. 10.1002/cncr.25531.21351083 10.1002/cncr.25531
58. le Thuy TT Matsumoto Y Thuy TT Hai H Suoh M Urahara Y Cytoglobin deficiency promotes liver cancer development from hepatosteatosis through activation of the oxidative stress pathway Am J Pathol 2015 185 4 1045 1060 10.1016/j.ajpath.2014.12.017 25665792
le Thuy TT, Matsumoto Y, Thuy TT, Hai H, Suoh M, Urahara Y, et al. Cytoglobin deficiency promotes liver cancer development from hepatosteatosis through activation of the oxidative stress pathway. Am J Pathol. 2015;185(4):1045–60. 10.1016/j.ajpath.2014.12.017.25665792 10.1016/j.ajpath.2014.12.017
59. Shen L Zhang J Zheng Z Yang F Liu S Wu Y PHGDH Inhibits Ferroptosis and Promotes Malignant Progression by Upregulating SLC7A11 in Bladder Cancer Int J Biol Sci 2022 18 14 5459 5474 10.7150/ijbs.74546 36147463
Shen L, Zhang J, Zheng Z, Yang F, Liu S, Wu Y, et al. PHGDH Inhibits Ferroptosis and Promotes Malignant Progression by Upregulating SLC7A11 in Bladder Cancer. Int J Biol Sci. 2022;18(14):5459–74. 10.7150/ijbs.74546.36147463 10.7150/ijbs.74546
60. Chen H Wen J Zhang W Ma W Guo Y Shen L circKDM1A suppresses bladder cancer progression by sponging miR-889–3p/CPEB3 and stabilizing p53 mRNA iScience. 2024 27 4 109624 10.1016/j.isci.2024.109624 38632984
Chen H, Wen J, Zhang W, Ma W, Guo Y, Shen L, et al. circKDM1A suppresses bladder cancer progression by sponging miR-889–3p/CPEB3 and stabilizing p53 mRNA. iScience. 2024;27(4):109624. 10.1016/j.isci.2024.109624.38632984 10.1016/j.isci.2024.109624
