Skip to main content
Download PDF

CCNB1 is a novel prognostic biomarker and promotes proliferation, migration and invasion in Wilms tumor

BMC Medical Genomics volume 16, Article number: 189 (2023) Cite this article

Abstract

Background

Wilms tumour (WT) is a mixed type of embryonal tumour that usually occurs in early childhood. However, our knowledge of the pathogenesis or progression mechanism of WT is inadequate, and there is a scarcity of beneficial therapeutic strategies.

Methods

High-throughput RNA sequencing was employed in this study to identify differentially expressed genes (DEGs) in clinical tumor samples and matching normal tissues. The STRING database was utilized to build a protein-protein interaction (PPI) network, and the Cytohubba method was used to identify the top 10 highly related HUB genes. Then, the key genes were further screened by univariate COX survival analysis. Subsequently, the XCELL algorithm was used to evaluate the tumour immune infiltration. RT-PCR, WB, and IF were used to verify the expression level of key genes in clinical tissues and tumour cell lines. Finally, the function of the key gene was further verified by loss-of-function experiments.

Results

We initially screened 1612 DEGs, of which 1030 were up-regulated and 582 were down-regulated. The GO and KEGG enrichment analysis suggested these genes were associated with ‘cell cycle’, ‘DNA replication’. Subsequently, we identified 10 key HUB genes, among them CCNB1 was strongly related to WT patients’ overall survival. Multiple survival analyses showed that CCNB1 was an independent indicator of WT prognosis. Thus, we constructed a nomogram of CCNB1 combined with other clinical indicators. Single gene GSEA and immune infiltration analysis revealed that CCNB1 was associated with the degree of infiltration or activation status of multiple immune cells. TIDE analysis indicated that this gene was correlated with multiple key immune checkpoint molecules and TIDE scores. Finally, we validated the differential expression level of CCNB1 in an external gene set, the pan-cancer, clinical samples, and cell lines. CCNB1 silencing significantly inhibited the proliferation, migration, and invasive capabilities of WIT-49 cells, also, promoted apoptosis, and in turn induced G2 phase cell cycle arrest in loss-of-function assays.

Conclusion

Our study suggests that CCNB1 is closely related to WT progression and prognosis, and serves as a potential target.

Peer Review reports

Introduction

Wilms tumour (WT) is a mixed type of embryonal tumour that accounts for more than 7% of childhood malignancies and 90% of pediatric renal tumours [1]. WT usually occurs in early childhood, and most cases are sporadic, with only 1-2% of familial cases [2]. It is well known that the tumour grows rapidly with the main clinical manifestation of a huge abdominal mass. The recurrence rate of WT is nearly 15%, and the long-term survival rate is only 50% [3]. As medical treatments have advanced, radical resection, radiotherapy, and chemotherapy have improved the 5-year survival rate to over 90% [4]. In addition, the side effects of high-intensity radiotherapy and chemotherapy seriously affect the quality of life of WT patients [5, 6]. However, our understanding of the pathogenesis and progression mechanisms of Wilms tumor is still insufficient, and there is a lack of effective therapeutic targets [7,8,9,10]. Hence, it is critical to investigate the molecular processes behind WT development and to establish efficient biomarkers for better WT therapy and overall prognosis.

Since its first appearance in 2005 [11], high-throughput sequencing has become widely available with the popularization of sequencing technologies, and RNA-seq has become one of the most common applications of sequencing technology [12]. Differential gene expression analysis is one of the primary applications of RNA-seq, which can reveal the functional and potential molecular mechanisms of differentially expressed genes identified [13]. Importantly, differential gene expression analysis is useful for identifying potential biomarkers for cancer [14]. Bioinformatics has been widely applied to various cancer studies and has been proven effective and reliable for diagnosing cancer and targeting new tumor markers [15]. Many cancers, including Wilms’ tumor, occur due to genetic mutations. Using bioinformatics methods, we can identify and validate key genes related to cancer and provide new targets for diagnosis and treatment [16]. In this study, we used high-throughput sequencing combined with bioinformatics methods to screen key genes associated with the prognosis of Wilms’ tumor and verified their specific expression and impact on the malignant biological phenotype of Wilms’ tumor through a series of basic experiments. Figure 1 displays our main steps.

Fig. 1
figure 1

Flow chart of analysis for the present study

Full size image

Materials and methods

Clinical specimens

All programs were approved by the Ethics Committee of the affiliated Children’s Hospital of Chongqing Medical University and all patients diagnosed with WT had not received radiotherapy or chemotherapy before surgery. We collected WT tissue and normal kidney tissue from 37 patients who underwent surgery from January 2019 to September 2022 (Table 1). Each sample was immediately frozen in liquid nitrogen and preserved. All samples were identified as WT by pathology tissue department specialists, and all WT samples were classified by pathologists following the Children Oncology Group (COG) staging system [17]. We randomly selected 8 pairs of tumor samples and paired normal tissue samples for mRNA sequencing to obtain the expression data of mRNA in WT. The rest of the samples were used for gene expression validation.

Table 1 The clinicopathological features of WT patients. WT, Wilms tumor; TNM, tumor-node-metastasis; FH, favorable histology; UH, unfavorable histology
Full size table

RNA sequencing (RNA-seq) of WT tissue

The Illumina NovaSeq 6000 was used for library construction and sequencing at Shanghai Sinomics Corporation (Illumina, USA). Detailed procedures can be referred to the previously published paper [18].

Identification of DEGs

To measure mRNA levels, we utilized fragments per kilobase million (FPKM) to reflect the expression of various genes, which could be determined as follow:

$$\text{F}\text{P}\text{K}\text{M}=\frac{\text{t}\text{o}\text{t}\text{a}\text{l} \text{e}\text{x}\text{o}\text{n} \text{F}\text{r}\text{a}\text{g}\text{m}\text{e}\text{n}\text{t}\text{s}}{\text{m}\text{a}\text{p}\text{p}\text{e}\text{d} \text{r}\text{e}\text{a}\text{d}\text{s} \left(\text{M}\text{i}\text{l}\text{l}\text{i}\text{o}\text{n}\text{s}\right)\times \text{e}\text{x}\text{o}\text{n} \text{l}\text{e}\text{n}\text{g}\text{t}\text{h} \left(\text{K}\text{B}\right)}$$

The fragments within each gene interval were first counted using Stringtie software [19]. Finally, FPKM values were then calculated for each gene, and the software edgeR [20] was used to analyze the differences in gene expression between groups. The final P-value was calculated for differential expression between tumour and normal tissues. Meanwhile, the differential expression fold change (FC) was calculated based on the FPKM values. Ultimately, DEGs were evaluated using both |logFC| > 1 and Q-value 0.05 as filtering principles.

Enrichment Analysis and hub genes identification

To investigate the underlying mechanics of DEGs in WT, we used the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genome (KEGG) [21] to explore potential biological processes and pathways using the R package ‘clusterProfiler’ [22]. In addition, DEGs were uploaded to the STRING database [23] to construct a protein-protein interaction (PPI) network. Cytoscape [24] based on the Cytohubba algorithm [25] was adopted to identify the top ten highly correlated HUB genes.

Gene Set Enrichment Analysis (GSEA)

GSEA is a functional annotation tool for understanding the underlying biological pathways of different biological phenotypes or states [26]. In this study, we performed GSEA analysis using gene sets based on GO, KEGG, and ImmuneSigDB. We grouped the samples based on the median expression of CCNB1 and identified the biological processes and pathways that were significantly associated with high and low expression groups.

Survival analysis

The TARGET database was used to acquire gene transcription and patient information on WT individuals. To compare survival differences across groups, log-rank was used to assess the Kaplan-Meier (KM) survival analysis, and time-dependent receiver operating characteristic (ROC) curve analysis was done to predict accuracy. The survival outcome metric used was overall survival (OS). We conducted a multivariate Cox survival analysis of CCNB1 and other clinical indicators, resulting in three independent prognostic factors: stage, gender, and CCNB1. We then constructed a risk score (Riskscore) based on the impact of each indicator on prognosis in the COX regression model, and divided the subjects into high-risk and low-risk groups based on the median of the risk score for further survival analysis. Subsequently, the independent prognostic value of CCNB1 was constructed and evaluated in conjunction with clinical indicators.

Assessment of Immune Infiltration

The Stromalscore and Immunescore, which reflect the extent of TME-related cell infiltration in WT tumor tissue, were calculated using the “IOBR” R package, employing the ESTIMATE algorithm[27, 28]. The ESTIMATE algorithm was developed based on single-sample gene set enrichment analysis (ssGSEA). The TME consists of various types of immune cells. To assess the heterogeneous cellular landscape of the TME, cell type enrichment scores were evaluated. In the TARGET and GSE31403 datasets, the XCELL method [29] was used to calculate the degree of immune cell infiltration and immune score. In addition, in the TARGET cohort, we employed the Tumour Immune Dysfunction and Exclusion (TIDE) score to measure treatment responsiveness to immune checkpoint inhibitors [30].

Expression Verification

We obtained and downloaded the WT expression dataset (GSE66405) from the GEO platform for expression validation of CCNB1. Subsequently, a uniformly normalized pan-cancer dataset was downloaded from the Sangerbox platform (http://sangerbox.com/home.html), from which the expression data of the CCNB1 (ENSG00000134057) gene in each sample were extracted. Meanwhile, expression differences in clinical tissues and cell lines were calculated using an unpaired t-test.

RT-PCR

Tissue RNA was extracted from previously collected WT tissue and normal kidney tissue in 37 patients. Cellular RNA was extracted from WT cells (WIT-49) and normal renal epithelial cells (293T). RT-PCR methodology was performed according to previous literature [18]. The expression of CCNB1 and GAPDH was based on the formula 2^-ΔΔCt. Table 2 lists the primers that were utilized.

Table 2 Primers used for quantitative real time PCR.
Full size table

Cell transfection

Table 3 shows the CCNB1 targeting siRNA and negative control (NC) siRNA sequences (TSINGKE, Beijing, China). The WIT-49 cells were plated onto 6-well or 12-well plates and divided into four groups: experimental group for CCNB1siRNA-1, CCNB1siRNA-2, CCNB1siRNA-3, and negative control group with NC siRNA. Lipofectamine™ RNAiMAX (Invitrogen, USA) was used for transient transfection according to the manufacturer’s instructions.

Table 3 CCNB1 siRNA sequences
Full size table

Western blotting and immunofluorescence

The previous article [31] described experimental methods for western blot (Prior to antibody hybridization, we cut the membrane at corresponding band positions.) and immunofluorescence. The following antibodies were used in the experiments: primary antibody (Anti-CCNB1 Rabbit pAb, GB112098; Anti-GAPDH Rabbit pAb, GB11002, Servicebio, China), secondary antibody (HRP − conjugated secondary antibody, G1213-100UL, Servicebio, China).

Cell proliferation, Migration, and Invasion Assay

The CCK-8 counting method (MCE, HY-K0301) was used to determine cell viability. Approximately 1 × 104 WIT-49 cells were inoculated into 96-well plates and transfected with CCNB1siRNA-1 and NCsiRNA. Transwell invasion assay to determine cell invasion ability. A scratch wound healing motility assay was performed to determine cell migration using WIT-49 cells. A methodology was referenced from the previous literature of our team [32].

Flow Cytometry

Flow cytometry was used to determine the distribution of cell cycle phases and the amount of apoptosis using the BD detection kit. FlowJo software was used to handle all data.

Statistical analysis

The experimental data were analyzed and graphed using the statistical program GraphPad Prism 8. Sangerbox was used for all bioinformatics analyses. Every experiment was carried out three times. A p-value less than 0.05 was considered statistically significant.

Results

Identification of DEGs in WT tissue and functional Enrichment Analysis

We conducted mRNA sequencing analysis on 8 pairs of tumor samples and matched normal tissue samples, identifying distinct mRNA expression profiles in WT. According to the method’s guideline, a total of 1612 differentially expressed genes were revealed, with 1030 being up-regulated and 582 being down-regulated (Fig. 2B). Cluster heat map and PCA analysis showed that differential genes distinguished the tumour from normal tissue (Fig. 2A, C). Following GO and KEGG enrichment analysis, the DEGs were shown to be significantly enriched in biochemical processes such as ‘DNA replication’, ‘Cell cycle’. (Fig. 2D-H).

Fig. 2
figure 2

Identification of DEGs in WT Tissue and Functional Enrichment Analysis. (A) Heatmaps and cluster analysis of DEGs. Each row represents a gene, and each column represents a sample. Red indicates high expression, while blue indicates low expression. Letters represent tumor samples, and numbers represent normal samples (paired as follows: A-1, B-2, C-3, D-4, E-5, F-6, G-7, H-8). (B) Volcano plot of DEGs. Red represents upregulated genes (1030), blue represents downregulated genes (582). The horizontal dashed line represents an adjusted P-value of 0.05, and the vertical dashed line represents a logFC of 2. (C) PCA analysis shows a remarkable difference in transcriptomes between the tumors (red) and normal tissues (blue). (D-H) GO and KEGG enrichment analysis showed that DEGs were mainly enriched in biological pathways or processes such as ‘DNA replication’, ‘cell cycle’

Full size image

PPI Network Construction and hub genes identification

We constructed a PPI network of 1612 DEGs based on the STRING database. The PPI network includes 693 nodes and 6656 edges (Fig. 3A). We used cytoHubba to obtain the top 10 most connected Hub genes in the PPI network, including CCNB1, CDK1, CCNB2, CCNA2, BUB1, KIF11, AURKB, NCAPG, CDC45, and BUB1B (Fig. 3B).

Fig. 3
figure 3

Identification of prognosis-related HUB genes and construction of a nomogram. (A) A PPI network of 1612 DEGs based on the STRING database. Red represents upregulated genes, blue represents downregulated genes. (B) Identification of top ten Hub genes using cytoHubba algorithm. (C) Univariate Cox analysis of the 10 Hub genes showed that CCNB1 was significantly correlated with OS in WT patients. (D) KM survival curves show that high expression of CCNB1 was associated with poorer prognosis in WT. (E) Multiple COX regression analyses revealed that CCNB1 was a critical risk biomarker independent of clinical indicators. (F) Combining CCNB1 expression and clinical indicators to develop a nomogram that could predict patient survival (1-, 3-, and 5-year survival rates). (G) Combined CCNB1 expression and clinical indicators to assess patient risk scores. (H and I) ROC curve and correction curve showed that the prediction model had good prediction performance (1-, 3-, and 5-year survival rates)

Full size image

Survival and clinical relevance analysis

Univariate Cox analysis of 10 Hub genes showed that CCNB1 was significantly associated with OS of WT patients (Fig. 3C, D). Clinical correlation showed that the expression of CCNB1 was significantly associated with the pathological type and tumour progression (Table 4), suggesting that the gene may be involved in the malignant progression of WT. Multiple survival analyses showed that CCNB1 was an independent predictor of WT patients. Therefore, we constructed a nomogram (the C-index of the nomogram was 0.72) by combining other clinical indicators for predicting the near and long-term survival of patients (Fig. 3E, F). Subsequently, we calculated risk scores for each patient in a composite model with multiple indicators (Fig. 3G). The ROC curves and calibration curves showed good performance of the prediction model. The AUC values predicting 1-, 3-, and 5-year survivals were increased to 0.80, 0.75, and 0.73, respectively (Fig. 3H, I).

Table 4 Correlation between the expression of CCNB1 and the clinicopathological features of WT patients
Full size table

Analysis of the CCNB1-Associated Immune Infiltration

Single gene GSEA enrichment analysis revealed that this gene may be involved in regulating the activation and functional levels of multiple immune cells such as CD8+ T cells, CD4+ T cells, and NKT cells (Fig. 4A). Therefore, heat maps were generated by Single-sample gene set enrichment analysis (ssGSEA) to show the immune scores and relative abundance of immune cell subpopulations in the TARGET and GSE31403 datasets (Fig. 4B). The XCELL algorithm was used to calculate the ImmuneScore, StromaScore, MicroenvironmentScore, and immune infiltration level between high and low-expression subgroups of CCNB1 in the TARGET and GSE31403 cohorts (Fig. 4C). Compared to the low expression group in the TARGET cohort, the high CCNB1 expression group had a lower StromaScore. Furthermore, in the TARGET cohort, there was a significant increase in CD4 + memory T cells, CD4 + T cells, M2 macrophages, pro-B cells, γδ T cells, and Th2 cells in the CCNB1 high expression group. However, CD4 + naive T cells, CD4 + central memory (Tcm), and natural killer T (NKT) cells showed a significant increase in the CCNB1 low expression group. To explain the survival differences found in patients with good prognosis from an immune perspective, we further analyzed the differences in immune cell infiltration between CCNB1 high and low expression groups of WT patients in the GSE31403 cohort. In the GSE31403 cohort, the high expression group had lower ImmuneScore and MicroenvironmentScore compared to the low expression group. At the same time, the CCNB1 high expression group showed a significant increase in CD4 + memory T cells, CD8 + Tcm, pro-B cells, γδ T cells, and Th2 cells. The main lymphocyte subgroups involved in anti-tumor immunity, including CD4 + Tcm, CD4 + Tem, and CD8 + T cells, showed a significant increase in the CCNB1 low expression group. Taken together, these results suggest that the favorable prognosis associated with low CCNB1 expression may be partially attributed to anti-tumor immune activity. Notably, this gene showed a significant correlation with TIDE score and several immune checkpoint molecules (HMGB1, CD80, IL13, ENTPD1, BTN3A1, TGFB1, VEGFA, ICOSLG, ICAM1, BTN3A2, TNFRSF14, CX3CL1, TLR4, ADORA2A, TNFRSF18) (Fig. 4D, E), indicating that CCNB1 might be an immunological regulatory target.

Fig. 4
figure 4

Associations between CCNB1 expression and immune infiltration in TME. (A) single-gene GSEA enrichment analysis showed that CCNB1 was involved in regulating the activation and function of immune cells. (B) Single-sample gene set enrichment analysis to generate heat maps to display immune scores and relative abundance of immune cell subsets in TARGET and GSE31403 datasets. (C) XCELL algorithm was used to calculate the ImmuneScore, StromaScore, MicroenvironmentScore, and immune cell infiltration levels between CCNB1 high expression subsets and low expression subsets in TARGET and GSE31403 cohorts. (D and E) CCNBI was correlated with the TIDE score and several immune checkpoint molecules. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001

Full size image

Expression analysis in multiple datasets and clinical samples

To verify the key role of this gene, we evaluated CCNB1 expression levels in two different datasets. The results showed that this gene showed tumour-specific expression levels in two independent datasets (Fig. 5B). Pan-cancer analysis showed consistent results (Fig. 5A), suggesting this gene may be a key oncogene. Next, we performed RT-PCR in 37 WT cancer tissues and paired adjacent normal kidney tissues to verify the expression of CCNB1. The results showed that CCNB1 was significantly overexpressed in WT tissues compared with paraneoplastic tissues (Fig. 5C). Similarly, the expression of CCNB1 was significantly higher in the WT cell line WIT-49 than in the renal epithelial normal cell line 293T (Fig. 5D). Moreover, WB and IF results also showed that CCNB1 was highly expressed in WT cell lines and cancer tissues, respectively (Fig. 5E, F).

Fig. 5
figure 5

CCNB1 exhibits a tumor-specific expression pattern. (A) Pan-cancer analysis showed CCNB1 was highly expressed in a variety of tumors. (B) Two independent databases verified that CCNB1 was significantly overexpressed in WT. (C and D) RT-PCR verified that CCNB1 was highly expressed in WT tissues and cell lines. (E and F) WB and immunofluorescence were used to demonstrate the expression of CCNB1 at the protein level. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001

Full size image

Silencing of CCNB1 significantly inhibited the Proliferation, Migration, and Invasion ability of WT cells

Single-gene GSEA enrichment analysis further suggested that this gene may be involved in biological processes such as regulation of ‘cell cycle’, ‘DNA replication’ (Fig. 6A, B). To demonstrate the role of CCNB1 in WT cells, three siRNA segments targeting CCNB1 were introduced into WIT-49 cells. After transfection for 48 h, we found that the knockdown efficiency of the first segment of si-CCNB1 was the highest (Fig. 6C). When compared to the si-NC group, si-CCNB1 protein expression was considerably lower in WIT-49 cells (Fig. 6D). Furthermore, cell function studies demonstrated that silencing this gene greatly reduced WIT-49 cell proliferation, migration, and invasive capacity (Fig. 6E-G).

Fig. 6
figure 6

CCNB1 promotes the malignant biological behavior of Wilms tumor cells via multiple phenotypes. (A-B) single-gene GSEA enrichment analysis showed that CCNB1 was mainly involved in regulating cell cycle-related processes. (C) RT-PCR analysis of CCNB1 expression in WIT-49 cells after transfection with siRNAs. (D) Western blot analysis of CCNB1 expression in WIT-49 cells after transfection with si-CCNB1-1. (E-G) CCK-8 assay, cell Migration assay, and Transwell assay were used to demonstrate that silencing CCNB1 could significantly inhibit the proliferation, migration, and invasion ability of WIT-49 cells. (H and I) Flow Cytometry identified that silencing CCNB1 caused G2 phase arrest and promoted WIT-49 cell apoptosis. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001

Full size image

In addition, this gene belongs to the family of cell cycle proteins. Therefore, we further verified whether this gene is involved in regulating the cell cycle and apoptosis levels of WIT-49 cells. Relatively perspective, silencing this gene caused G2 phase arrest (Fig. 6H) and promoted the apoptosis of WIT-49 cells (Fig. 6I) by flow cytometry.

Discussion

Wilms tumor is the most prevalent malignant solid tumor in children’s urology, and the exact pathophysiology is unknown [33]. Despite advances in clinical therapy for WT, overall survival has not increased considerably, which can be related to a paucity of molecular markers for successful diagnosis and treatment. As a result, it is critical to investigate WT molecular markers in order to increase patient survival. Researchers may now integrate several bioinformatics methodologies to extensively investigate the main pathophysiology and clinical diagnosis or prognosis of many diseases at the molecular level, thanks to the fast growth of diverse bioinformatics databases and high-throughput studies [34]. Therefore, the integration of RNA sequencing data and mining in the database has become an important means to explore the pathogenesis of WT and speculate the possible markers for diagnosis and treatment.

In this study, we integrated and analyzed the high-throughput sequencing data of WT tissues and adjacent normal tissues, and found that most of the differential genes were up-regulated. GO and KEGG pathway analysis showed that DEGs were enriched in multiple pathways related to malignant biology, including ‘DNA replication’, ‘cell cycle’. To further identify key genes that play important roles, we established a PPI network based on the STRING database and screened 10 HUB genes. CCNB1 was found to be a risk factor for poor prognosis in WT patients by univariate COX analysis. Multifactorial analysis similarly demonstrated that this gene was an independent risk factor. This suggests a potential biological role for CCNB1 in WT progression. It has been shown that CCNB1 is a prognostic factor for overall survival and metastasis-free survival in breast cancer [35]. Several studies have confirmed that CCNB1 can be used as a diagnostic or prognostic biomarker for rhabdomyosarcoma, hepatocellular carcinoma, and meningioma [36,37,38]. Critically, our results of CCNB1 validation from transcriptome and protein levels and combined with clinicopathological information showed that CCNB1 was significantly highly expressed in WT. Functional experiments confirmed that silencing CCNB1 could inhibit the proliferation, invasion, and migration of WT cells, suggesting a crucial role of this gene in WT. Interestingly, there have been reports confirming the carcinogenic effect of high expression of CCNB1 in various cancers, including renal cancer, breast cancer, pancreatic cancer, hepatocellular carcinoma, and cervical cancer [39,40,41,42]. Several studies have shown that CCNB1 is a potential target for tumour intervention [43, 44]. CCNB1 has been reported to be considered other potentially useful genes for targeting hepatocellular carcinoma [45, 46]. Therefore, it is speculated that this gene may be a potential key therapeutic target for WT.

The cell cycle-related factor CCNB1 belongs to the family of cell cycle proteins [47]. It is well known that one of the distinguishing features of cancer is cell cycle dysregulation, leading to the unrestricted proliferation of cancer cells [48, 49]. In this study, GSEA enrichment analysis suggested that the gene promoted tumour progression by regulating the cell cycle, and further experiments proved that silencing the gene caused cell cycle G2 arrest and apoptosis. These findings encourage us to speculate that CCNB1 is essential for cell cycle progression and proliferation. To sum up, CCNB1 may regulate WT tumour progression through the cell cycle pathway.

Cancer patients usually have a large number of T cells, but most of them have lost their function [50]. One study found that CCNB1 caused T cell-dependent antibody responses in patients with cancer and precancerous lesions, suggesting that this gene is an important player in the immune control of tumour growth [51]. In addition, CCNB1, which is aberrantly expressed in patients with breast, lung, head, and neck cancers, can be recognized by antibodies and T cells as tumour antigens [52]. CCNB1 was discovered by Kao et al. as a common human epithelial tumor-associated antigen recognized by T lymphocytes [53]. Latner et al. [54] elucidated the enhanced expression of CCNB1 in virus-specific memory CD8+ T cells. It is worth noting that another important finding of our study is that CCNB1 is associated with immune scores and multiple immune cell infiltration levels. Single gene GSEA analysis suggests that the gene may be involved in the regulation of a variety of key immune cells, which may affect tumour immunity and lead to poor prognosis and can be used as an important indicator of cancer prognosis.

We demonstrated that CCNB1 is a promising prognostic marker and potential therapeutic target for WT by high-throughput sequencing combined with bioinformatics analysis and experimental validation. Despite the large sample bioinformatics analysis and clinical sample validation performed in this study, certain limitations remain: First, our prognostic analysis and model are based on a public cohort, and retrospective innate characteristics hinder clinical applicability. Therefore, prospective studies are needed for further validation before clinical application. Second, the downstream mechanism by which CCNB1 exerts its oncogene function remains unclear requiring further in vivo and in vitro experiments for validation. Third, the immune regulatory function of this gene in the tumour microenvironment is based on bioinformatics analysis, and the exact mechanism needs to be further explored.

In conclusion, our findings imply that CCNB1 is a significant prognostic biomarker and a possible therapeutic target for Wilms tumor.

Data Availability

The datasets generated during and/or analyses during the current study are available in the TARGET and GEO Database. The datasets presented in this study can be found in online repositories. The name of the repository and accession numbers can be found below: NCBI; GSE31403 and GSE66405. We intend to share individual deidentified participant data. RNA sequencing data from eight tumor samples and paired normal tissue samples were uploaded to the gene expression omnibus database (https://www.ncbi.nlm.nih.gov/geo/, under the accession: GSE197047). you can access it now and it’s valid forever.

Abbreviations

WT:

Wilms tumour

COG:

Children Oncology Group

TIDE:

Tumour Immune Dysfunction and Exclusion

References

  1. Treger TD, Chowdhury T, Pritchard-Jones K, Behjati S. The genetic changes of Wilms tumour. Nat Rev Nephrol. 2019;15(4):240–51.

    PubMed  Google Scholar 

  2. Vujanić GM, Gessler M, Ooms A, Collini P, Coulomb-l’Hermine A, D’Hooghe E, de Krijger RR, Perotti D, Pritchard-Jones K, Vokuhl C, et al. The UMBRELLA SIOP-RTSG 2016 Wilms tumour pathology and molecular biology protocol. Nat reviews Urol. 2018;15(11):693–701.

    Google Scholar 

  3. Maschietto M, Piccoli FS, Costa CM, Camargo LP, Neves JI, Grundy PE, Brentani H, Soares FA, Camargo B, Carraro DM. Gene expression analysis of blastemal component reveals genes associated with relapse mechanism in Wilms tumour. Eur J cancer (Oxford England: 1990). 2011;47(18):2715–22.

    CAS  Google Scholar 

  4. Davidoff AM. Wilms tumor. Adv Pediatr. 2012;59(1):247–67.

    PubMed  PubMed Central  Google Scholar 

  5. Szychot E, Apps J, Pritchard-Jones K. Wilms’ tumor: biology, diagnosis and treatment. Translational Pediatr. 2014;3(1):12–24.

    Google Scholar 

  6. Breslow NE, Lange JM, Friedman DL, Green DM, Hawkins MM, Murphy MF, Neglia JP, Olsen JH, Peterson SM, Stiller CA, et al. Secondary malignant neoplasms after Wilms tumor: an international collaborative study. Int J Cancer. 2010;127(3):657–66.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Pelosi A, Fiore PF, Di Matteo S, Veneziani I, Caruana I, Ebert S, Munari E, Moretta L, Maggi E, Azzarone B. Pediatric Tumors-Mediated Inhibitory Effect on NK Cells: The Case of Neuroblastoma and Wilms’ Tumors. Cancers 2021, 13(10).

  8. Jain J, Sutton KS, Hong AL. Progress Update in Pediatric Renal Tumors. Curr Oncol Rep. 2021;23(3):33.

    PubMed  Google Scholar 

  9. Brok J, Mavinkurve-Groothuis AMC, Drost J, Perotti D, Geller JI, Walz AL, Geoerger B, Pasqualini C, Verschuur A, Polanco A, et al. Unmet needs for relapsed or refractory Wilms tumour: mapping the molecular features, exploring organoids and designing early phase trials - a collaborative SIOP-RTSG, COG and ITCC session at the first SIOPE meeting. Eur J cancer (Oxford England: 1990). 2021;144:113–22.

    Google Scholar 

  10. van den Heuvel-Eibrink MM, Hol JA, Pritchard-Jones K, van Tinteren H, Furtwängler R, Verschuur AC, Vujanic GM, Leuschner I, Brok J, Rübe C, et al. Position paper: Rationale for the treatment of Wilms tumour in the UMBRELLA SIOP-RTSG 2016 protocol. Nat reviews Urol. 2017;14(12):743–52.

    Google Scholar 

  11. Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen Y-J, Chen Z, et al. Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005;437(7057):376–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Moon M, Nakai K. Stable feature selection based on the ensemble L 1 -norm support vector machine for biomarker discovery. BMC Genomics. 2016;17(Suppl 13):1026.

    PubMed  PubMed Central  Google Scholar 

  13. Oshlack A, Robinson MD, Young MD. From RNA-seq reads to differential expression results. Genome Biol. 2010;11(12):220.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Govindarajan M, Wohlmuth C, Waas M, Bernardini MQ, Kislinger T. High-throughput approaches for precision medicine in high-grade serous ovarian cancer. J Hematol Oncol. 2020;13(1):134.

    PubMed  PubMed Central  Google Scholar 

  15. Pucker B, Schilbert HM, Schumacher SF. Integrating Molecular Biology and Bioinformatics Education. J Integr Bioinform 2019, 16(3).

  16. Cui C, Li L, Zhen J. Bioinformatic analysis reveals the key pathways and genes in early-onset breast cancer. Med Oncol. 2018;35(5):67.

    PubMed  Google Scholar 

  17. Dome JS, Perlman EJ, Graf N. Risk stratification for wilms tumor: current approach and future directions. American Society of Clinical Oncology educational book American Society of Clinical Oncology Annual Meeting 2014:215–223.

  18. Tian XM, Xiang B, Zhang ZX, Li YP, Shi QL, Li MJ, Li Q, Yu YH, Lu P, Liu F, et al. The Regulatory Network and Role of the circRNA-miRNA-mRNA ceRNA Network in the progression and the Immune response of Wilms Tumor based on RNA-Seq. Front Genet. 2022;13:849941.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.

    CAS  PubMed  Google Scholar 

  20. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinf (Oxford England). 2010;26(1):139–40.

    CAS  Google Scholar 

  21. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–d613.

    CAS  PubMed  Google Scholar 

  24. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):11.

    Google Scholar 

  26. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Zeng D, Ye Z, Shen R, Yu G, Wu J, Xiong Y, Zhou R, Qiu W, Huang N, Sun L, et al. IOBR: Multi-Omics Immuno-Oncology Biological Research to Decode Tumor Microenvironment and Signatures. Front Immunol. 2021;12:687975.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, Treviño V, Shen H, Laird PW, Levine DA, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.

    PubMed  Google Scholar 

  29. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017;18(1):220.

    PubMed  PubMed Central  Google Scholar 

  30. Jiang P, Gu S, Pan D, Fu J, Sahu A, Hu X, Li Z, Traugh N, Bu X, Li B, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. 2018;24(10):1550–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Mi T, Jin L, Zhang Z, Wang J, Li M, Zhanghuang C, Tan X, Wang Z, Tian X, Xiang B, et al. DNA hypermethylation-regulated CX3CL1 reducing T cell infiltration indicates poor prognosis in Wilms Tumour. Front Oncol. 2022;12:882714.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Luo X, Dong J, He X, Shen L, Long C, Liu F, Liu X, Lin T, He D, Wei G. MiR-155-5p exerts tumor-suppressing functions in Wilms tumor by targeting IGF2 via the PI3K signaling pathway. Biomed pharmacotherapy = Biomedecine pharmacotherapie. 2020;125:109880.

    CAS  Google Scholar 

  33. Luo X, Deng C, Liu F, Liu X, Lin T, He D, Wei G. HnRNPL promotes Wilms tumor progression by regulating the p53 and Bcl2 pathways. OncoTargets and therapy. 2019;12:4269–79.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Buzdin A, Sorokin M, Garazha A, Glusker A, Aleshin A, Poddubskaya E, Sekacheva M, Kim E, Gaifullin N, Giese A, et al. RNA sequencing for research and diagnostics in clinical oncology. Sem Cancer Biol. 2020;60:311–23.

    CAS  Google Scholar 

  35. Niméus-Malmström E, Koliadi A, Ahlin C, Holmqvist M, Holmberg L, Amini RM, Jirström K, Wärnberg F, Blomqvist C, Fernö M, et al. Cyclin B1 is a prognostic proliferation marker with a high reproducibility in a population-based lymph node negative breast cancer cohort. Int J Cancer. 2010;127(4):961–7.

    PubMed  Google Scholar 

  36. Li Q, Zhang L, Jiang J, Zhang Y, Wang X, Zhang Q, Wang Y, Liu C, Li F. CDK1 and CCNB1 as potential diagnostic markers of rhabdomyosarcoma: validation following bioinformatics analysis. BMC Med Genom. 2019;12(1):198.

    CAS  Google Scholar 

  37. Li N, Li L, Chen Y. The Identification of Core Gene Expression Signature in Hepatocellular Carcinoma. Oxidative medicine and cellular longevity 2018, 2018:3478305.

  38. Bie L, Zhao G, Ju Y, Zhang B. Integrative genomic analysis identifies CCNB1 and CDC2 as candidate genes associated with meningioma recurrence. Cancer Genet. 2011;204(10):536–40.

    CAS  PubMed  Google Scholar 

  39. Ding K, Li W, Zou Z, Zou X, Wang C. CCNB1 is a prognostic biomarker for ER + breast cancer. Med Hypotheses. 2014;83(3):359–64.

    CAS  PubMed  Google Scholar 

  40. Ikuerowo SO, Kuczyk MA, Mengel M, van der Heyde E, Shittu OB, Vaske B, Jonas U, Machtens S, Serth J. Alteration of subcellular and cellular expression patterns of cyclin B1 in renal cell carcinoma is significantly related to clinical progression and survival of patients. Int J Cancer. 2006;119(4):867–74.

    CAS  PubMed  Google Scholar 

  41. Zhou L, Li J, Zhao YP, Cui QC, Zhou WX, Guo JC, You L, Wu WM, Zhang TP. The prognostic value of cyclin B1 in pancreatic cancer. Med Oncol (Northwood Lond Engl). 2014;31(9):107.

    Google Scholar 

  42. Bai X, Wang W, Zhao P, Wen J, Guo X, Shen T, Shen J, Yang X. LncRNA CRNDE acts as an oncogene in cervical cancer through sponging miR-183 to regulate CCNB1 expression. Carcinogenesis. 2020;41(1):111–21.

    PubMed  Google Scholar 

  43. Li M, Shang H, Wang T, Yang SQ, Li L. Huanglian decoction suppresses the growth of hepatocellular carcinoma cells by reducing CCNB1 expression. World J Gastroenterol. 2021;27(10):939–58.

    PubMed  PubMed Central  Google Scholar 

  44. Zeng H, Ji J, Song X, Huang Y, Li H, Huang J, Ma X. Stemness related genes revealed by Network Analysis Associated with Tumor Immune Microenvironment and the clinical outcome in Lung Adenocarcinoma. Front Genet. 2020;11:549213.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Lin J, Hou Y, Huang S, Wang Z, Sun C, Wang Z, He X, Tam NL, Wu C, Wu L. Exportin-T promotes tumor proliferation and invasion in hepatocellular carcinoma. Mol Carcinog. 2019;58(2):293–304.

    CAS  PubMed  Google Scholar 

  46. Zhang Q, Su R, Shan C, Gao C, Wu P. Non-SMC Condensin I Complex, Subunit G (NCAPG) is a Novel Mitotic Gene required for Hepatocellular Cancer Cell Proliferation and Migration. Oncol Res. 2018;26(2):269–76.

    PubMed  PubMed Central  Google Scholar 

  47. Miyazaki T, Arai S. Two distinct controls of mitotic cdk1/cyclin B1 activity requisite for cell growth prior to cell division. Cell cycle (Georgetown Tex). 2007;6(12):1419–25.

    CAS  PubMed  Google Scholar 

  48. Engeland K. Cell cycle arrest through indirect transcriptional repression by p53: I have a DREAM. Cell Death Differ. 2018;25(1):114–32.

    CAS  PubMed  Google Scholar 

  49. Massagué J. G1 cell-cycle control and cancer. Nature. 2004;432(7015):298–306.

    PubMed  Google Scholar 

  50. Jiang W, He Y, He W, Wu G, Zhou X, Sheng Q, Zhong W, Lu Y, Ding Y, Lu Q, et al. Exhausted CD8 + T cells in the Tumor Immune Microenvironment: New Pathways to Therapy. Front Immunol. 2020;11:622509.

    CAS  PubMed  Google Scholar 

  51. Suzuki H, Graziano DF, McKolanis J, Finn OJ. T cell-dependent antibody responses against aberrantly expressed cyclin B1 protein in patients with cancer and premalignant disease. Clin cancer research: official J Am Association Cancer Res. 2005;11(4):1521–6.

    CAS  Google Scholar 

  52. Vella LA, Yu M, Fuhrmann SR, El-Amine M, Epperson DE, Finn OJ. Healthy individuals have T-cell and antibody responses to the tumor antigen cyclin B1 that when elicited in mice protect from cancer. Proc Natl Acad Sci USA. 2009;106(33):14010–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Kao H, Marto JA, Hoffmann TK, Shabanowitz J, Finkelstein SD, Whiteside TL, Hunt DF, Finn OJ. Identification of cyclin B1 as a shared human epithelial tumor-associated antigen recognized by T cells. J Exp Med. 2001;194(9):1313–23.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Latner DR, Kaech SM, Ahmed R. Enhanced expression of cell cycle regulatory genes in virus-specific memory CD8 + T cells. J Virol. 2004;78(20):10953–9.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We appreciate the free use of TARGET and GEO databases.

Funding

This work was supported by the general project of clinical medicine research of Children’s Hospital of Chongqing Medical University (NCRC-2019-GP-08), the Medical Scientific Research Project of Chongqing (NO.2022GDRC009), the Natural Science Foundation of Chongqing (cstc2021jcyj-msxmX0345), and the general project of basic research of Key Laboratory of children’s developmental diseases of the Ministry of Education (GBRP-202105).

Author information

Authors and Affiliations

  1. Ministry of Education Key Laboratory of Child Development and Disorders, Chongqing Key Laboratory of Pediatrics, National Clinical Research Center for Child Health and Disorders, International Science and Technology Cooperation Base of Child Development and Critical Disorders, Children’s Hospital of Chongqing Medical University, Chongqing, P.R. China

    Bin Xiang, Mei-Lin Chen, Zhi-Qiang Gao, Tao Mi, Qin-Lin Shi, Jun-Jun Dong, Xiao-Mao Tian, Feng Liu & Guang-Hui Wei

  2. Department of Urology, Chongqing Key Laboratory of Children Urogenital Development and Tissue Engineering, Children’s Hospital of Chongqing Medical University, Chongqing, P.R. China

    Bin Xiang, Mei-Lin Chen, Zhi-Qiang Gao, Tao Mi, Qin-Lin Shi, Jun-Jun Dong, Xiao-Mao Tian, Feng Liu & Guang-Hui Wei

Authors
  1. Bin Xiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Mei-Lin Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Zhi-Qiang Gao
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Tao Mi
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Qin-Lin Shi
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Jun-Jun Dong
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Xiao-Mao Tian
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Feng Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Guang-Hui Wei
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

FL, XMT, and GHW contributed to creation and planning; XMT, BX, JJD, and TM contributed to data analysis; XMT, ZQG, QLS, and BX contributed to bioinformatics analysis; BX, and MLC participated in the experimental procedure. XMT, BX, and TM contributed to the manuscript writing. All authors reviewed the manuscript.

Corresponding authors

Correspondence to Xiao-Mao Tian or Feng Liu.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in strict accordance with the rules of the Declaration of Helsinki. The studies involving human participants were reviewed and approved by the Ethics Committee of the Childrens Hospital Affiliated to Chongqing Medical University, Chongqing, China (File No. 2022-50). All parents of the children who participated in this study have provided the written informed consent.

Consent for publication

Not applicable.

Informed consent

Informed consent was obtained from all individuals included in this study.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Rights and permissions

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/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xiang, B., Chen, ML., Gao, ZQ. et al. CCNB1 is a novel prognostic biomarker and promotes proliferation, migration and invasion in Wilms tumor. BMC Med Genomics 16, 189 (2023). https://doi.org/10.1186/s12920-023-01627-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-023-01627-3

Keywords

BMC Medical Genomics

ISSN: 1755-8794

Contact us