Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 26 November 2024
Sec. Cancer Molecular Targets and Therapeutics
This article is part of the Research Topic Renewed Insight into Cancer Mechanism and Therapy View all articles

Comprehensive prognostic gene identification and functional characterization of GRAMD1A in Wilms tumor: development of risk prediction models and therapeutic implications

Qiang ZengQiang Zeng1Junfeng TaoJunfeng Tao1Lilu QinLilu Qin1Yong ZengYong Zeng2Zhong LiuZhong Liu1Mingxian Xu*Mingxian Xu3*Linshan Zeng*Linshan Zeng2*
  • 1Department of Pediatric Surgery, Jiangxi Maternal and Child Health Hospital, Jiangxi Children’s Medical Center, Nanchang, Jiangxi, China
  • 2Department of Pediatric Surgery, The First Affiliated Hospital of Gannan Medical University, Ganzhou, Jiangxi, China
  • 3Department of Musculoskeletal Oncology, The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, Guangdong, China

Background: Wilms tumor (WT) is the most common pediatric kidney cancer, with survival rates exceeding 90% in localized cases. However, advanced or recurrent WT remains difficult to treat due to poor prognosis and limited knowledge of its molecular mechanisms. Gene expression profiling has shown promise in identifying prognostic markers and therapeutic targets. This study aimed to identify key prognostic genes and pathways in WT, construct risk prediction models, and validate their role in tumor progression.

Methods: RNA sequencing and clinical data from 136 WT patients were obtained from the TARGET database. Differential gene expression analysis was conducted using GEO datasets GSE11024 and GSE66405 to compare WT and normal kidney tissues. Identified differentially expressed genes (DEGs) underwent Gene Ontology (GO) and KEGG pathway enrichment analysis to explore biological functions and pathways associated with WT progression. Univariate Cox regression was used to assess the association between DEGs and overall survival (OS) and progression-free survival (PFS). LASSO regression models were developed for risk stratification, and model accuracy was evaluated using time-dependent ROC curves. External validation confirmed key hub genes, while functional assays in WT cell lines (WiT-49) assessed the role of GRAMD1A in tumor behavior.

Results: A total of 3,395 DEGs were identified, with 1,564 upregulated and 1,831 downregulated genes. Enrichment analyses revealed significant pathways involved in cell cycle regulation and metabolic reprogramming. Six key genes (GRAMD1A, PLXNA3, SPR, EBAG9, RBM47, and RIDA) were associated with both OS and PFS. LASSO models demonstrated strong predictive performance, with GRAMD1A identified as a major risk factor. External validation confirmed differential expression, and functional assays showed that GRAMD1A silencing significantly inhibited WT cell viability, proliferation, migration, and invasion.

Conclusions: This study identifies novel prognostic genes and potential therapeutic targets in WT. GRAMD1A, SPR, EBAG9, RBM47, and RIDA play critical roles in WT progression, with GRAMD1A as a key oncogenic factor, offering potential for risk stratification and future therapeutic intervention.

Introduction

Wilms tumor (WT) is the most prevalent renal malignancy in infants and young children (8, 9). WT is closely linked to early nephrogenesis, exhibiting morphological and genetic similarities with the early stages of kidney development (1012). Standard treatment regimens primarily involve a multimodal approach comprising surgery, chemotherapy, and radiotherapy. Despite significant reductions in mortality through advancements in these therapies, the issues of tumor recurrence and treatment-related toxicity persist, highlighting the urgent need for more individualized therapeutic strategies (13, 14). While nephrectomy remains a widely accepted and effective treatment, it presents significant challenges in pediatric patients, particularly in those undergoing active growth or with bilateral involvement (1518).

Clonal evolution in tumor cells is shaped by both morphological and genetic progression. Numerous copy number alterations (CNAs) serve as prognostic biomarkers for patient stratification. For instance, the WT1 tumor suppressor gene, which is critical in kidney development, has been shown to play a pivotal role in Wilms tumor (WT) pathogenesis. Mutations in WT1 are associated with abnormal tumor formation and developmental defects (1921). Thus, optimizing therapeutic strategies based on clinical and biological risk factors is key to advancing personalized treatment and improving patient outcomes (5, 11, 22). A deeper understanding of tumor biology is essential for elucidating the genetic basis of WT and facilitating the development of personalized therapeutic approaches. Although significant progress has been made, the molecular mechanisms driving WT development, and the identification of potential therapeutic targets remain incompletely understood. In this study, we aim to identify potential prognostic and therapeutic targets, as well as the underlying molecular mechanisms of WT, by leveraging data from the TARGET and GEO databases. Additionally, we seek to construct a reliable clinical prognostic model based on key genes, offering promising therapeutic avenues for WT patients.

Methods

Dataset retrieval and preprocessing

The overall analytical process is depicted in Figure 1. RNA-seq data were sourced from the TARGET database, encompassing clinical data from 136 Wilms tumor (WT) cases. In addition, the GSE11024 (1), GSE66405 (2),GSE73209 (3), and GSE110696 (4) datasets were retrieved from the NCBI GEO database (http://www.ncbi.nih.gov/geo). Table 1 outlines the specific applications of these datasets in the study.

Figure 1
www.frontiersin.org

Figure 1. Schematic overview of the analytical workflow employed in this study.

Table 1
www.frontiersin.org

Table 1. Overview of principal attributes of the datasets employed.

Transcriptomic analysis of differential gene expression

Differential expression analysis was performed using the R package “limma” (version 4.0.3), which employs a negative binomial model to estimate log2 fold changes and p-values for each gene. Genes with a p-value of < 0.05 were considered significantly differentially expressed. Volcano plots were generated to illustrate the distribution of differentially expressed genes (DEGs) between Wilms tumor (WT) and normal kidney tissues, utilizing data from the GEO database. The raw data, downloaded in MINiML format, underwent differential expression analysis via the “limma” package in R software. To reduce false positives, adjusted p-values were calculated using the Benjamini-Hochberg method, with DEGs defined as those with an adjusted p-value < 0.05 and a log2 fold change > 1.3 or < -1.3. The overlapping DEGs from datasets GSE66405 and GSE11024 were compared using a Venn diagram, and the overlap coefficient was calculated to assess consistency between the datasets. Additionally, validation of key hub genes was conducted using two independent datasets, GSE73209 and GSE110696 (5).

Functional enrichment analysis: gene ontology and KEGG pathways

To characterize the biological functions and pathways associated with differentially expressed genes (DEGs) identified from GEO datasets, functional enrichment analyses were performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The “clusterProfiler” package in R was employed for this purpose, with significant pathways defined by an adjusted p-value < 0.05 using the Benjamini-Hochberg correction. Enrichment analyses were conducted across three categories: biological processes (BP), molecular functions (MF), and cellular components (CC), with a minimum of 10 genes per category included in the analysis. Gene sets were retrieved from the Molecular Signatures Database (MSigDB) (6) to ensure comprehensive coverage. Visualization of the enriched GO terms and KEGG pathways was accomplished using the “enrichplot” R pack (7).

Survival analysis

To evaluate the prognostic significance of key hub genes, survival analysis was performed using the TARGET Wilms tumor (WT) cohort. Patients were stratified into high- and low-expression groups based on the median expression level of each gene. Kaplan-Meier survival curves were generated to assess the relationship between gene expression levels and overall survival (OS) or progression-free survival (PFS). Statistical significance was determined using the log-rank test, with p-values < 0.05 considered significant. Hazard ratios (HR) and 95% confidence intervals (CI) were calculated using univariate Cox proportional hazards regression to quantify the impact of gene expression on patient survival outcomes.

Development and validation of the prognostic risk model

The prognostic risk model was developed by integrating gene expression data from the TARGET Wilms tumor cohort with patient survival outcomes. Key survival-associated genes were identified using the least absolute shrinkage and selection operator (LASSO) regression, combined with 10-fold cross-validation to optimize the model. Risk scores were subsequently calculated for each patient based on the regression coefficients and corresponding gene expression values, with the calculation method detailed in a previous study.

Cell culture

HEK 293T cells were acquired from the American Type Culture Collection (ATCC, USA), and WiT49 cells were obtained as previously reported (5). WiT49 cells were cultured in a 1:1 mixture of high-glucose Dulbecco’s Modified Eagle Medium (DMEM) and Nutrient Mixture F-12 (DMEM/F-12), supplemented with 10% fetal bovine serum (FBS), 100 U/mL penicillin, and 100 µg/mL streptomycin, as described in previous studies. HEK 293T cells were maintained in DMEM with 10% FBS, supplemented with the same antibiotics. Both cell lines were routinely tested for mycoplasma contamination and confirmed to be mycoplasma-free before use in experiments, ensuring cell line authenticity and reliability for downstream applications.

Transfection

Transfection was performed using Lipofectamine 3000 (Invitrogen, USA) according to the manufacturer’s protocol. WiT49 cells were seeded in 6-well plates at a density optimized for 70-80% confluency at the time of transfection. Small interfering RNA (siRNA) targeting GRAMD1A (si-GRAMD1A) and a negative control siRNA (si-NC) were purchased from IGEbio (Guangzhou, China). The specific sequences used for each siRNA are provided in Table 2. Transfection efficiency was confirmed 48 hours post-transfection by quantifying GRAMD1A mRNA and protein expression levels using RT-qPCR and Western blot, respectively. This ensured effective knockdown of GRAMD1A, validating the subsequent functional assays.

Table 2
www.frontiersin.org

Table 2. Oligonucleotide sequences.

Western blotting

Total proteins were extracted from WiT49 and HEK 293T cells using RIPA lysis buffer supplemented with protease inhibitors. Protein concentrations were quantified using a BCA protein assay kit (Thermo Fisher Scientific). Equal amounts of protein (25 μg per lane) were subjected to separation by SDS-PAGE and subsequently transferred onto 0.22 μm PVDF membranes (Millipore). The membranes were blocked in 5% non-fat milk diluted in Tris-buffered saline with 0.1% Tween-20 (TBST) for 1 hour at room temperature. Following blocking, membranes were incubated overnight at 4°C with primary antibodies against GRAMD1A (1:1000, NBP1-93730, Novus) or α-Tubulin (1:5000, 11224-1-AP, Proteintech), with α-Tubulin serving as the loading control. After three washes in TBST, membranes were incubated with HRP-conjugated secondary antibodies (1:5000, Thermo Fisher Scientific) for 1 hour at room temperature. Protein bands were visualized using enhanced chemiluminescence (ECL) reagent and imaged with a ChemiDoc XRS+ imaging system (Bio-Rad). The relative expression levels of proteins were quantified using ImageJ software, normalized to α-Tubulin to ensure consistent loading across samples.

RT-PCR

Total RNA was isolated from Wilms’ tumor (WiT-49) cells and normal renal epithelial (HEK 293T) cells using TRIzol reagent (Invitrogen, USA), following the manufacturer’s protocol. The extracted RNA was quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific). Complementary DNA (cDNA) was synthesized from 1 μg of total RNA using the M-MLV reverse transcription kit (M1705, Promega). Quantitative PCR (qPCR) was performed on the synthesized cDNA using SYBR Master Mix (DRR041B, TAKARA) with a LightCycler 480 II instrument (Roche). The qPCR cycling conditions included an initial denaturation at 95°C for 5 minutes, followed by 40 cycles of denaturation at 95°C for 10 seconds, annealing at 60°C for 20 seconds, and extension at 72°C for 30 seconds. Relative gene expression levels of GRAMD1A and GAPDH were calculated using the 2^-ΔΔCt method, with GAPDH as the internal control for normalization. Primer sequences used for qRT-PCR are detailed in Table 3. All experiments were conducted in triplicate to ensure reproducibility.

Table 3
www.frontiersin.org

Table 3. Primers for Quantitative Real-Time PCR (qRT-PCR).

Cell viability and colony formation assay

WiT49 cells, following transfection with either si-GRAMD1A or si-NC, were seeded into 96-well plates at a density of 3,000 cells per well. Cell viability was assessed at 24, 48, and 72 hours post-transfection using the Cell Counting Kit-8 (CCK-8; Dojindo, Japan). After adding 10 μL of CCK-8 reagent per well, cells were incubated at 37°C in a humidified 5% CO2 atmosphere for 2 hours, and absorbance was measured at 450 nm using a microplate reader.

For the colony formation assay, transfected WiT49 cells were plated at a density of 500 cells per well in 6-well plates and cultured for 12 days, with the culture medium refreshed every 4 days. Colonies formed were fixed with 4% formaldehyde for 15 minutes, stained with 0.1% crystal violet, and counted manually under a microscope.

Transwell assay

Transfected WiT49 cells were resuspended in serum-free DMEM medium and seeded at a density of 5 × 104 cells per insert into the upper chamber of transwell inserts (8 μm pore size; Corning, USA). For the migration assay, inserts without Matrigel coating were used, while for the invasion assay, the membranes were coated with 100 μL of Matrigel (1:8 dilution, BD Bioscience, USA) and allowed to solidify for 1 hour at 37°C. The lower chambers were filled with medium containing 15% FBS, serving as a chemoattractant. After 24 hours of incubation at 37°C in a humidified incubator with 5% CO2, non-migratory/invading cells on the upper side of the membrane were carefully removed using a cotton swab. Cells that migrated or invaded to the lower surface were fixed with 4% paraformaldehyde for 10 minutes, stained with 0.1% crystal violet (Solarbio, China) for 30 minutes at room temperature, washed with phosphate-buffered saline (PBS), and counted under a phase-contrast microscope in five randomly selected fields per membrane.

Statistical analysis

Data analysis and visualization were conducted in R (version 4.0.3). The Kruskal-Wallis test evaluated continuous variable differences, and the Wilcoxon test compared two groups. Hazard ratios (HR) and 95% confidence intervals (CI) were calculated via Cox regression using the survival package. Kaplan-Meier analysis and log-rank tests were used to compare survival curves. Statistical significance was set at p < 0.05 (*p < 0.05, **p < 0.01, ***p < 0.001).

Result

Identification of differentially expressed genes in Wilms’ tumors patients

RNA sequencing data from the GEO datasets GSE11024 (Wilms tumor, n=27; normal kidney, n=12) and GSE66405 (Wilms tumor, n=28; normal kidney, n=4) were analyzed to determine differentially expressed genes (DEGs). In the GSE11024 dataset, a total of 7,061 DEGs were identified (Figure 2A), while 5,525 DEGs were detected in the GSE66405 dataset (Figure 2B). To evaluate the consistency of these findings, a Venn diagram was constructed to identify intersecting genes, revealing 3,419 common DEGs between the two datasets (Figure 2C). After eliminating duplicates, 3,395 unique DEGs with consistent expression profiles were retained, comprising 1,564 upregulated genes (Figure 2D) and 1,831 downregulated genes (Figure 2E). These DEGs are likely involved in critical biological processes relevant to the pathogenesis of Wilms tumor.

Figure 2
www.frontiersin.org

Figure 2. Differential gene expression analysis. (A) Volcano plot depicting differentially expressed genes (DEGs) in dataset GSE11024. (B) Volcano plot depicting DEGs in dataset GSE66405. (C–E) Venn diagrams illustrating intersecting gene sets: (C) Total intersecting genes between datasets, (D) Overlapping upregulated genes, and (E) Overlapping downregulated genes.

Functional enrichment analyses of DEGs

To elucidate the underlying biological roles and potential mechanisms of the identified DEGs, we performed KEGG pathway and Gene Ontology (GO) enrichment analyses separately on the 1,564 upregulated and 1,831 downregulated genes.

For the upregulated DEGs, KEGG pathway enrichment analysis identified 18 significantly enriched pathways, with the most notable ones being the Cell cycle, Wnt signaling pathway, and Nucleocytoplasmic transport (Figure 3A). GO enrichment analysis revealed 511 significantly enriched Biological Process (BP) terms, primarily associated with organelle fission, nuclear division, and chromosome segregation (Figure 3B). A total of 48 enriched terms for Molecular Function (MF) were found, mainly related to tubulin binding, ATP hydrolysis activity, and catalytic activity, acting on DNA (Figure 3C). Additionally, 99 enriched terms for Cellular Component (CC) were identified, including chromosomal region, spindle, and condensed chromosome (Figure 3D).

Figure 3
www.frontiersin.org

Figure 3. Functional Enrichment Analysis of Differentially Expressed Genes (DEGs). (A, E) KEGG pathway enrichment analysis for differentially expressed genes: (A) Pathways enriched for upregulated genes, (E) Pathways enriched for downregulated genes. (B, F) Gene Ontology (GO) enrichment analysis for biological processes involving DEGs: (B) Biological processes enriched for upregulated genes, (F) Biological processes enriched for downregulated genes. (C, G) GO enrichment analysis for molecular functions of DEGs: (C) Molecular functions enriched for upregulated genes, (G) Molecular functions enriched for downregulated genes. (D, H) GO enrichment analysis for cellular components of DEGs: (D) Cellular components enriched for upregulated genes, (H) Cellular components enriched for downregulated genes.

For the downregulated DEGs, KEGG pathway enrichment analysis identified 65 significantly enriched pathways, primarily related to Thermogenesis, Chemical carcinogenesis - reactive oxygen species, and Oxidative phosphorylation (Figure 3E). GO enrichment analysis of these DEGs identified 670 enriched Biological Process (BP) terms, with key processes involving the small molecule catabolic process, generation of precursor metabolites and energy, and fatty acid metabolic process (Figure 3F). For Molecular Function (MF), 168 enriched terms were found, including active transmembrane transporter activity, ion transmembrane transporter activity, and lyase activity (Figure 3G). Moreover, 97 enriched Cellular Component (CC) terms were identified, such as mitochondrial matrix, mitochondrial inner membrane, and apical part of the cell (Figure 3H). These enrichment analyses suggest that dysregulation of the cell cycle, Wnt signaling, and metabolic reprogramming, particularly involving oxidative phosphorylation and fatty acid metabolism, play critical roles in the progression of Wilms tumor (WT), highlighting potential therapeutic targets.

Prognostic gene identification and analysis in Wilms tumor patients

To identify key prognostic genes in Wilms tumor (WT), clinical data and RNA-seq data from 136 WT patients in the TARGET database were analyzed. A univariate Cox regression analysis was performed to assess the association between the 3,395 DEGs and overall survival (OS) as well as progression-free survival (PFS). As a result, 95 genes were found to be significantly associated with OS (Table 4), and 154 genes were significantly associated with PFS (Table 5). Further analysis using a Venn diagram revealed six genes that were significantly related to both OS and PFS: GRAMD1A, PLXNA3, SPR, EBAG9, RBM47, and RIDA (Figure 4A). Univariate Cox regression analysis was further performed on the six candidate genes identified from the intersection to assess their association with overall survival (Figure 4B).

Table 4
www.frontiersin.org

Table 4. Univariate analysis of genes associated with overall survival (OS).

Table 5
www.frontiersin.org

Table 5. Univariate analysis of genes associated with progression-free survival (PFS).

Figure 4
www.frontiersin.org

Figure 4. Identification of prognosis-related genes in Wilms tumor. (A) Venn diagram illustrating the overlap of genes significantly associated with overall survival (OS) and progression-free survival (PFS). (B) Univariate Cox regression analysis of six candidate genes identified from the intersection, assessing their association with overall survival. (C-H) Kaplan-Meier survival curves depicting the relationship between the expression levels of each intersecting gene and overall survival in patients with Wilms tumor. High and low expression groups are stratified by median expression levels, demonstrating the prognostic value of each gene.

Subsequent analysis of the expression levels of these six genes in relation to OS demonstrated that high expression of GRAMD1A (Figure 4C) and PLXNA3 (Figure 4D) was associated with poor OS, while low expression of SPR (Figure 4E), EBAG9 (Figure 4F), RBM47 (Figure 4G), and RIDA (Figure 4H) was correlated with unfavorable OS outcomes. Similarly, the expression levels of these six genes were analyzed in relation to PFS, as shown in Figure 5A. The relationship between the expression of these genes and PFS revealed that high expression of GRAMD1A (Figure 5B) and PLXNA3 (Figure 5C) was linked to poor PFS, whereas low expression of SPR (Figure 5D), EBAG9 (Figure 5E), RBM47 (Figure 5F), and RIDA (Figure 5G) was associated with poor PFS outcomes. These findings suggest that both high expression of GRAMD1A and PLXNA3, and low expression of SPR, EBAG9, RBM47, and RIDA are associated with poorer OS and PFS in WT patients, highlighting their potential as prognostic biomarkers.

Figure 5
www.frontiersin.org

Figure 5. Validation of Prognostic Significance of Six Key Genes in Relation to Progression-Free Survival (PFS). (A) Univariate Cox regression analysis of six intersecting genes to determine their prognostic significance for PFS. (B–G) Kaplan-Meier survival curves illustrating the association between the expression levels of each of the intersecting genes and PFS in Wilms tumor patients. Each plot shows the impact of high versus low gene expression on patient outcomes, providing insight into their potential prognostic value.

Construction of the prognostic risk model and analysis

To further evaluate the predictive capacity of the six identified prognostic genes, two LASSO regression-based prognostic risk models (OS/PFS) were constructed to assess their effectiveness in predicting overall survival (OS) and progression-free survival (PFS) in WT patients (Figure 6). The LASSO regression analysis for OS identified four genes (GRAMD1A, SPR, EBAG9, and RBM47) that were significantly associated with OS (Figures 6A–E). A risk score formula was developed as follows:

Figure 6
www.frontiersin.org

Figure 6. Evaluation of Prognostic Risk Score Models Based on Six Key Prognostic Genes in Wilms Tumor (WT) Patients. (A, F) LASSO (Least Absolute Shrinkage and Selection Operator) regression coefficients for the selected prognostic genes, with optimization performed using the lambda parameter to minimize prediction error. (B, G) Partial likelihood deviance as a function of log(λ), demonstrating the selection of the optimal penalty parameter in the LASSO model. (C, H) Comprehensive representation of risk score distribution, patient survival status, and a heatmap illustrating the expression levels of key genes within the developed risk model, highlighting the distinction between high-risk and low-risk groups. (D, I) Kaplan-Meier survival analysis comparing overall survival (OS) or progression-free survival (PFS) between high-risk and low-risk groups, stratified by the median risk score, emphasizing the prognostic value of the risk model. (E, J) Time-dependent Receiver Operating Characteristic (ROC) curve analyses assessing the predictive accuracy of the gene signatures over time for OS or PFS, with Area Under the Curve (AUC) values indicating the robustness of the prognostic model. (A–E) represent the LASSO regression model for overall survival, whereas (F–J) depict the LASSO regression model for progression-free survival. Each sub-panel contributes to validating the model’s efficiency in predicting clinical outcomes for WT patients.

Riskscore=(0.3325)*GRAMD1A+(-0.1101)*SPR+(-0.5246)*EBAG9+(-0.053)*RBM47(lambda.min=0.0344). These results suggest that GRAMD1A is a risk factor for OS, whereas SPR, EBAG9, and RBM47 are protective factors. Figure 6C shows the distribution of risk scores, survival status, and a heatmap of the four key gene expression levels. WT patients were stratified into high- and low-risk groups based on the median risk score, with the low-risk group showing significantly better OS (p=0.0217, HR=1.895) (Figure 6D). Time-dependent ROC analysis was used to assess the model’s accuracy in predicting WT patient survival. The area under the curve (AUC) at 1, 3, 5, and 10 years was 0.547, 0.703, 0.686, and 0.568, respectively (Figure 6E), indicating the model’s potential for OS prediction. The LASSO regression analysis for PFS identified five key genes (Figures 6F–J). The following equation was used to calculate the PFS risk score:

Riskscore=(0.0657)GRAMD1A+(0.0769)SPR+(0.4489)EBAG9+(0.0582)RBM47+(0.5172)RIDA (lambda.min=0.0302)

In the PFS risk model, GRAMD1A was identified as a risk factor, while SPR, EBAG9, RBM47, and RIDA were protective factors. Patients in the low-risk group had significantly better PFS compared to the high-risk group (p < 0.001, HR = 2.678) (Figure 6I). The AUC values for PFS prediction at 1, 3, 5, and 8 years were 0.690, 0.829, 0.853, and 0.822, respectively (Figure 6J). In summary, the PFS risk model, constructed using GRAMD1A, SPR, EBAG9, RBM47, and RIDA, demonstrates better predictive accuracy than the OS risk model, making it a valuable tool for prognostic assessment in WT patients. The identified prognostic genes, particularly GRAMD1A, SPR, EBAG9, RBM47, and RIDA, offer significant potential as both prognostic biomarkers and therapeutic targets for improving personalized treatment in Wilms tumor patients.

Validation of prognostic key hub gene expression in Wilms tumor using external datasets

To validate the differential expression of the five key genes identified in the prognostic risk model, we utilized two external datasets, GSE73209 and GSE110696, to assess their expression levels in Wilms tumor (WT) and normal kidney tissues (Figure 7). In the GSE73209 validation dataset, GRAMD1A (Figure 7A) was significantly overexpressed in WT samples compared to normal kidney tissues. Conversely, RIDA (Figure 7B), RBM47 (Figure 7C), and SPR (Figure 7E) exhibited higher expression levels in normal kidney tissues, while EBAG9 (Figure 7D) showed no significant differential expression between WT and normal tissues. Consistent findings were obtained from the GSE110696 dataset, further corroborating these expression patterns (Figures 7F–J). These results reinforce the robustness and prognostic relevance of the constructed risk prediction model, reinforcing their potential as biomarkers for improving risk stratification and personalized treatment strategies in Wilms tumor patients.

Figure 7
www.frontiersin.org

Figure 7. Validation of the Expression Patterns of Five Key Hub Genes in External Datasets. (A, F) Validation of GRAMD1A gene expression in Wilms Tumor (WT) samples compared to normal kidney tissues. (B, G) Validation of RIDA gene expression levels in WT versus normal tissues. (C, H) Validation of RBM47 gene expression to assess differences between WT and normal kidney samples. (D, I) EBAG9 gene expression analysis to confirm differential expression in WT. (E, J) SPR gene expression comparison between WT and normal tissues. (A–E) Validation using the external dataset GSE73209, highlighting expression levels in WT samples versus normal kidney tissues. (F–J) Validation using the external dataset GSE110696, demonstrating the robustness and consistency of differential expression for the hub genes identified in our study. These validations confirm the expression trends of the five key hub genes identified, supporting their potential role in Wilms Tumor pathology and their utility as prognostic markers.

Cell functional verification

From the four key genes identified in the above analysis, GRAMD1A, which was highly expressed in tumor tissues, was selected for further functional investigation. GRAMD1A expression was confirmed to be upregulated at both the RNA and protein levels in WT cells (Figure 8A). To explore the functional role of GRAMD1A, transfection experiments were performed in WiT-49 cells using siRNA. Cells were transfected with either si-NC or si-GRAMD1A, and knockdown efficiency was assessed after 48 hours (Figure 8B). Functional assays revealed that silencing GRAMD1A significantly inhibited cell viability, proliferation, migration, and invasion in WiT-49 cells (Figures 8C–E).

Figure 8
www.frontiersin.org

Figure 8. Functional Characterization of GRAMD1A in WiT-49 Wilms Tumor Cells. (A) GRAMD1A expression levels in WiT-49 cells, demonstrating baseline expression prior to gene silencing. (B) Verification of GRAMD1A knockdown efficiency in WiT-49 cells transfected with GRAMD1A-specific siRNA (siGRAMD1A) or negative control siRNA (siNC). GRAMD1A mRNA levels were measured using RT-qPCR, while protein levels were assessed via Western blot analysis. (C) Colony formation assay showing the effects of GRAMD1A silencing on the proliferation capability of WiT-49 cells, highlighting reduced colony formation in GRAMD1A knockdown cells compared to siNC. (D) CCK-8 assay results depicting the effect of GRAMD1A knockdown on cell viability over time in WiT-49 cells. Cells transfected with siGRAMD1A exhibited significantly reduced viability compared to control cells, indicating GRAMD1A’s role in promoting cell proliferation. (E) Transwell migration and invasion assays used to assess the impact of GRAMD1A silencing on the migratory and invasive behavior of WiT-49 cells. GRAMD1A knockdown led to a marked reduction in both migration and invasion compared to the negative control group. These results collectively demonstrate the role of GRAMD1A in regulating the proliferation, migration, and invasion of WiT-49 cells, further supporting its potential as a therapeutic target in Wilms Tumor. (*p < 0.05, **p < 0.01, ***p < 0.001).

Discussion

Wilms tumor (WT) is the most prevalent pediatric renal malignancy (8, 23). Epidemiological studies have revealed variations in WT incidence across different geographic regions and ethnic groups (24), with the disease accounting for approximately 5% of all childhood malignancies in individuals under 15 years of age (25). WT is more common in females, with the incidence peaking in males at 17.9 cases per million person-years at one year of age (26). As an embryonal tumor, WT is thought to arise from developmental arrest or disruption during kidney formation. However, despite significant research efforts, the genetic underpinnings of WT remain incompletely understood, limiting advancements in therapeutic approaches. The identification of prognostic biomarkers and elucidation of the molecular mechanisms driving WT are critical for improving our understanding of this disease and optimizing treatment strategies. In this study, we provide novel insights into the molecular pathogenesis of WT and identify potential prognostic markers that may be instrumental in enhancing patient outcomes. By integrating clinical data and gene expression profiles from multiple datasets, we identified key differentially expressed genes (DEGs) and developed a robust prognostic risk model for WT patients. These findings not only contribute to a deeper understanding of WT biology but also highlight promising molecular targets for therapeutic intervention.

In this study, differential gene expression analysis was performed using transcriptomic data from the GEO database, identifying 3,419 differentially expressed genes (DEGs), including 1,564 upregulated and 1,831 downregulated genes. Subsequently, by integrating data from the TARGET database, six key genes significantly associated with Wilms tumor (WT) prognosis were identified. Among these, GRAMD1A and PLXNA3 were found to be overexpressed in WT and correlated with poorer prognosis, while SPR, EBAG9, RBM47, and RIDA showed reduced expression, which was linked to unfavorable outcomes. These six genes were further used to develop LASSO regression models to predict overall survival (OS) and progression-free survival (PFS) in WT patients. In the OS prediction model, GRAMD1A was identified as a risk factor, whereas SPR, EBAG9, and RBM47 acted as protective factors, with the model achieving an AUC of 0.703 for 3-year OS prediction. In the PFS prediction model, GRAMD1A was similarly identified as a risk factor, while SPR, EBAG9, RBM47, and RIDA were protective factors, resulting in AUC values of 0.690, 0.829, 0.853, and 0.822 at 1, 3, 5, and 8 years, respectively. These results demonstrate the strong prognostic value of these genes for predicting patient outcomes. Notably, the PFS model exhibited superior prognostic value compared to the OS model. Furthermore, our model demonstrated a significant improvement in prognostic accuracy over our previous model (5), as well as those developed by other researchers (27, 28). External validation using two independent datasets confirmed the differential expression of these genes in WT samples compared to normal kidney tissues. Specifically, GRAMD1A was significantly overexpressed in tumor tissues, while RIDA, RBM47, and SPR showed higher expression levels in normal kidney tissues. These findings collectively highlight the prognostic significance of these genes in WT and suggest their potential as targets for future therapeutic interventions.

GRAM Domain Containing 1A (GRAMD1A) mediates non-vesicular cholesterol transport from the plasma membrane (PM) to the endoplasmic reticulum (ER), serving as a key cholesterol transporter (29, 30). Previous studies have demonstrated that GRAMD1A is recruited to autophagosome initiation sites, and its cholesterol transport activity is inhibited by autophagy-related proteins. This inhibition is sufficient to suppress autophagosome biogenesis, indicating that GRAMD1A-mediated cholesterol transport is essential for autophagosome formation (31). In hepatocellular carcinoma (HCC), high GRAMD1A expression is associated with poor prognosis, and it has been shown to promote the self-renewal of liver cancer stem cells and drive hepatocarcinogenesis through the upregulation of STAT5 (32). Although the role of GRAMD1A in tumorigenesis has been largely unexplored in other cancer types, our study is the first to demonstrate that elevated GRAMD1A expression is linked to poor overall survival (OS) and progression-free survival (PFS) in Wilms tumor (WT) patients. Furthermore, GRAMD1A showed substantial diagnostic value in our LASSO-based prognostic model. To elucidate the functional role of GRAMD1A in WT, in vitro functional assays were conducted, revealing that GRAMD1A knockdown significantly inhibited tumor cell proliferation, invasion, and migration.

Reactive Intermediate Imine Deaminase A Homolog (RIDA), also referred to as UK114, was first identified in a 1993 study on chromatin-associated proteins, noted for its unusual solubility in perchloric acid. Early structural analyses indicated potential substrate interactions and underscored the conserved nature of this protein family (33). RidA proteins possess conserved enamine/imine deaminase activity, and the absence of RidA in bacteria, plants, and yeast leads to phenotypes such as nutrient deficiencies, mitochondrial maintenance defects, and disruptions in one-carbon (C1) metabolism (34, 35). In oncological research, RIDA/UK114 has been identified as a tumor antigen (3638). Studies have shown that RIDA/UK114 expression is tissue-specific, with the highest levels found in the liver and kidneys. In hepatocellular carcinoma (HCC), UK114 expression is significantly downregulated at both the mRNA and protein levels, and its decreased expression correlates with the Edmondson-Steiner grade of tumor differentiation, making it a potential biomarker for HCC staging (39). Despite these insights, the precise biological functions of RIDA remain incompletely understood. In this study, we demonstrated that RIDA expression is significantly downregulated in Wilms tumor (WT) compared to normal kidney tissues. Patients with lower RIDA expression had poorer prognoses, and in our prognostic model, higher RIDA expression was identified as a protective factor in WT.

RNA Binding Motif Protein 47 (RBM47) is a single-stranded RNA-binding protein that plays critical roles in various RNA processes, including alternative splicing, RNA stability, and RNA editing (4042). RBM47 is essential for embryonic endoderm development (43), and conditional expression of RBM47 alleles has been shown to cause fetal intestinal developmental defects and growth retardation (40, 44). As a novel and evolutionarily conserved RNA-binding protein (RBP) in vertebrates, RBM47 is increasingly recognized for its role as a tumor regulator (45). It has been reported to be downregulated in papillary thyroid carcinoma (PTC) tissues and cells, and its overexpression can induce autophagy and suppress PTC cell proliferation (46). In hepatocellular carcinoma (HCC), RBM47 upregulation significantly inhibits tumor progression in vitro, primarily through the upregulation of UPF1, serving as a tumor suppressor by acting as a DNA/RNA binding protein at both transcriptional and post-transcriptional levels (47). Furthermore, RBM47 has been shown to modulate intestinal injury and tumorigenesis by altering pathways related to cell proliferation, oxidative stress, and inflammation (48).

Sepiapterin reductase (SPR) is an aldo-keto reductase that catalyzes the NADPH-dependent reduction of pterin derivatives, playing a crucial role in the biosynthesis of tetrahydrobiopterin (BH4) (49). SPR is widely distributed across various tissues and is implicated in numerous diseases, including neurological dysfunction, chronic pain, cardiovascular diseases, and cancer (49). In hepatocellular carcinoma (HCC), elevated SPR expression is significantly associated with shorter patient survival, suggesting that SPR could serve as a potential prognostic marker for HCC. Moreover, SPR is involved in the progression of HCC through a non-enzymatic mechanism, controlling tumor development via the FoxO3a/Bim signaling pathway (50).

The enrichment analysis of upregulated and downregulated genes in Wilms tumor was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Upregulated DEGs were significantly enriched in the cell cycle, Wnt signaling pathway, and nucleocytoplasmic transport pathways, while GO analysis indicated involvement in organelle fission, nuclear division, and chromosome segregation. These results are consistent with the classification of WT as an embryonal malignancy resulting from developmental disruption, highlighting its connection to chromosomal segregation and cell proliferation (51, 52). In contrast, the downregulated DEGs were associated with pathways related to thermogenesis, chemical carcinogenesis, and oxidative phosphorylation. GO analysis further indicated key processes involved in small molecule catabolism, precursor metabolite generation, and fatty acid metabolism, suggesting their role in tumor cell energy metabolism. Previous studies have demonstrated that WT1 suppresses thermogenesis-related genes and modulates metabolic processes (53). Tumor cells often undergo metabolic reprogramming, with oxidative phosphorylation serving as a critical energy source. Its inhibition is emerging as a potential therapeutic target (54). Additionally, oxidative phosphorylation, mediated by macrophages and monocytes, has been closely associated with WT prognosis and clinical outcomes (55, 56).

Although the findings of this study provide valuable insights, several limitations must be addressed. First, the relatively small sample size derived from publicly available datasets, such as TARGET and GEO, may not adequately capture the genetic and clinical heterogeneity of Wilms tumor, potentially affecting the generalizability of our results. The retrospective nature of the study further introduces limitations, including possible biases related to the imbalance of clinicopathological features and treatment heterogeneity across patients. These factors limit the broader application of the findings, especially in diverse clinical settings. Future studies should focus on expanding the cohort size and incorporating more diverse populations to ensure the robustness of the results, while also utilizing prospective study designs to minimize bias. Moreover, although our in vitro assays successfully demonstrated the role of GRAMD1A in WT, additional in vivo studies are necessary to comprehensively evaluate its therapeutic potential and safety in a physiological context. The molecular mechanisms by which GRAMD1A contributes to tumor progression are not fully elucidated, underscoring the need for further mechanistic exploration. Additionally, while this study identified other prognostic genes (PLXNA3, SPR, EBAG9, RBM47, and RIDA) that could be of clinical significance, their precise roles in WT remain largely unexplored. Detailed mechanistic studies on these genes are warranted to better understand their contribution to WT pathogenesis and their interplay with known oncogenic pathways. Expanding this understanding will help validate these genes as potential biomarkers or therapeutic targets, ultimately enhancing personalized treatment approaches for Wilms tumor.

Conclusions

In summary, this study identified critical prognostic genes and molecular pathways associated with the progression of Wilms tumor, with GRAMD1A emerging as a key regulator. The prognostic risk models developed in this study provide valuable tools for enhancing the accuracy of outcome prediction and informing personalized therapeutic strategies for patients with Wilms tumor. Continued validation and investigation of these findings in clinical contexts may facilitate the development of novel targeted therapies and contribute to improved clinical outcomes for affected patients.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.

Ethics statement

Ethical approval was not required for the studies on humans in accordance with the local legislation and institutional requirements because only commercially available established cell lines were used.

Author contributions

QZ: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. JT: Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Validation, Writing – original draft. LQ: Conceptualization, Data curation, Methodology, Software, Visualization, Writing – original draft. YZ: Conceptualization, Data curation, Formal Analysis, Methodology, Validation, Visualization, Writing – original draft. ZL: Data curation, Formal Analysis, Investigation, Methodology, Validation, Writing – original draft. MX: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. LZ: Conceptualization, Formal Analysis, Investigation, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. Science and Technology Research Project of Jiangxi Provincial Department of Education, GJJ2203544. Jiangxi Children’s Medical Center (National Children’s Regional Medical Center), 2024B004.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Kort EJ, Farber L, Tretiakova M, Petillo D, Furge KA, Yang XJ, et al. The E2F3-Oncomir-1 axis is activated in Wilms’ tumor. Cancer Res. (2008) 68:4034–38. doi: 10.1158/0008-5472.CAN-08-0592

PubMed Abstract | Crossref Full Text | Google Scholar

2. Su C, Zheng J, Chen S, Tuo J, Su J, Ou X, et al. Identification of key genes associated with cancer stem cell characteristics in Wilms’ tumor based on bioinformatics analysis. Ann Transl Med. (2022) 10:1204. doi: 10.21037/atm-22-4477

PubMed Abstract | Crossref Full Text | Google Scholar

3. Karlsson J, Valind A, Jansson C, O’Sullivan MJ, Holmquist ML, Gisselsson D. Aberrant epigenetic regulation in clear cell sarcoma of the kidney featuring distinct DNA hypermethylation and EZH2 overexpression. Oncotarget. (2016) 7:11127–36. doi: 10.18632/oncotarget.7152

PubMed Abstract | Crossref Full Text | Google Scholar

4. Murphy AJ, Chen X, Pinto EM, Williams JS, Clay MR, Pounds SB, et al. Forty-five patient-derived xenografts capture the clinical and biological heterogeneity of Wilms tumor. Nat Commun. (2019) 10:5806. doi: 10.1038/s41467-019-13646-9

PubMed Abstract | Crossref Full Text | Google Scholar

5. Zeng Q, Liu T, Qin L, Wang C, Peng G, Liu Z, et al. Screening of potential hub genes involved in Kidney Wilms tumor via bioinformatics analysis and experimental validation. BMC Cancer. (2024) 24:771. doi: 10.1186/s12885-024-12541-x

PubMed Abstract | Crossref Full Text | Google Scholar

6. Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. (2015) 1:417–25. doi: 10.1016/j.cels.2015.12.004

PubMed Abstract | Crossref Full Text | Google Scholar

7. Wang N, Wei L, Liu D, Zhang Q, Xia X, Ding L, et al. Identification and validation of autophagy-related genes in diabetic retinopathy. Front Endocrinol (Lausanne). (2022) 13:867600. doi: 10.3389/fendo.2022.867600

PubMed Abstract | Crossref Full Text | Google Scholar

8. Pastore G, Znaor A, Spreafico F, Graf N, Pritchard-Jones K, Steliarova-Foucher E. Malignant renal tumors incidence and survival in European children (1978-1997): report from the Automated Childhood Cancer Information System project. Eur J Cancer. (2006) 42:2103–14. doi: 10.1016/j.ejca.2006.05.010

PubMed Abstract | Crossref Full Text | Google Scholar

9. Ooms AH, Gadd S, Gerhard DS, Smith MA, Guidry AJ, Meerzaman D, et al. Significance of TP53 mutation in wilms tumors with diffuse anaplasia: A report from the children’s oncology group. Clin Cancer Res. (2016) 22:5582–91. doi: 10.1158/1078-0432.CCR-16-0985

PubMed Abstract | Crossref Full Text | Google Scholar

10. Young MD, Mitchell TJ, Vieira BF, Tran M, Stewart BJ, Ferdinand JR, et al. Single-cell transcriptomes from human kidneys reveal the cellular identity of renal tumors. Science. (2018) 361:594–99. doi: 10.1126/science.aat1699

PubMed Abstract | Crossref Full Text | Google Scholar

11. Treger TD, Chowdhury T, Pritchard-Jones K, Behjati S. The genetic changes of Wilms tumor. Nat Rev Nephrol. (2019) 15:240–51. doi: 10.1038/s41581-019-0112-0

PubMed Abstract | Crossref Full Text | Google Scholar

12. Vujanic GM, Gessler M, Ooms A, Collini P, Coulomb-l’Hermine A, D’Hooghe E, et al. The UMBRELLA SIOP-RTSG 2016 Wilms tumor pathology and molecular biology protocol. Nat Rev Urol. (2018) 15:693–701. doi: 10.1038/s41585-018-0100-3

PubMed Abstract | Crossref Full Text | Google Scholar

13. Charles A. Phylogeny of wilms tumor? Ebiomedicine. (2016) 9:11–2. doi: 10.1016/j.ebiom.2016.06.044

PubMed Abstract | Crossref Full Text | Google Scholar

14. Green DM. The treatment of stages I-IV favorable histology Wilms’ tumor. J Clin Oncol. (2004) 22:1366–72. doi: 10.1200/JCO.2004.08.008

PubMed Abstract | Crossref Full Text | Google Scholar

15. Lopes RI, Ming J, Koyle MA, Grant R, Fonseca A, Lorenzo AJ. Zero-ischemia” Laparoscopic-assisted partial nephrectomy for the management of selected children with wilms tumor following neoadjuvant chemotherapy. Urology. (2017) 100:103–10. doi: 10.1016/j.urology.2016.08.051

PubMed Abstract | Crossref Full Text | Google Scholar

16. Godzinski J, Graf N, Audry G. Current concepts in surgery for Wilms tumor–the risk and function-adapted strategy. Eur J Pediatr Surg. (2014) 24:457–60. doi: 10.1055/s-0034-1396425

PubMed Abstract | Crossref Full Text | Google Scholar

17. Warmann SW, Godzinski J, van Tinteren H, Heij H, Powis M, Sandstedt B, et al. Minimally invasive nephrectomy for Wilms tumors in children - data from SIOP 2001. J Pediatr Surg. (2014) 49:1544–48. doi: 10.1016/j.jpedsurg.2014.06.005

PubMed Abstract | Crossref Full Text | Google Scholar

18. Schmidt A, Warmann SW, Urla C, Schaefer J, Fideler F, Fuchs J. Patient selection and technical aspects for laparoscopic nephrectomy in Wilms tumor. Surg Oncol. (2019) 29:14–9. doi: 10.1016/j.suronc.2019.02.007

PubMed Abstract | Crossref Full Text | Google Scholar

19. Royer-Pokora B, Wruck W, Adjaye J, Beier M. Gene expression studies of WT1 mutant Wilms tumor cell lines in the frame work of published kidney development data reveals their early kidney stem cell origin. PloS One. (2023) 18:e270380. doi: 10.1371/journal.pone.0270380

PubMed Abstract | Crossref Full Text | Google Scholar

20. Morrison AA, Viney RL, Saleem MA, Ladomery MR. New insights into the function of the Wilms tumor suppressor gene WT1 in podocytes. Am J Physiol Renal Physiol. (2008) 295:F12–17. doi: 10.1152/ajprenal.00597.2007

PubMed Abstract | Crossref Full Text | Google Scholar

21. Mrowka C, Schedl A. Wilms’ tumor suppressor gene WT1: from structure to renal pathophysiologic features. J Am Soc Nephrol. (2000) 11 Suppl 16:S106–15. doi: 10.1681/ASN.V11suppl_2s106

PubMed Abstract | Crossref Full Text | Google Scholar

22. Bitaraf M, Mahmanzar M, Zafari N, Mohammadpour H, Vasei M, Moradi ML, et al. The potential key genes and pathways associated with Wilms tumor in quest of proper candidates for diagnostic and therapeutic purposes. Sci Rep. (2022) 12:17906. doi: 10.1038/s41598-022-22925-3

PubMed Abstract | Crossref Full Text | Google Scholar

23. Aldrink JH, Heaton TE, Dasgupta R, Lautz TB, Malek MM, Abdessalam SF, et al. Update on wilms tumor. J Pediatr Surg. (2019) 54:390–97. doi: 10.1016/j.jpedsurg.2018.09.005

PubMed Abstract | Crossref Full Text | Google Scholar

24. Nakata K, Colombet M, Stiller CA, Pritchard-Jones K, Steliarova-Foucher E. Incidence of childhood renal tumors: An international population-based study. Int J Cancer. (2020) 147:3313–27. doi: 10.1002/ijc.33147

PubMed Abstract | Crossref Full Text | Google Scholar

25. Dome JS, Perlman EJ, Graf N. Risk stratification for Wilms tumor: current approach and future directions. Am Soc Clin Oncol Educ Book. (2014) 34:215–23. doi: 10.14694/EdBook_AM.2014.34.215

PubMed Abstract | Crossref Full Text | Google Scholar

26. Zeng Q, Xu M, Zeng L, Tao J, Qin L, Zeng Y, et al. Correction to Lancet Oncol. Lancet Oncol. (2017) 18:719–31. doi: 10.1016/S1470-2045(17)30369-8

PubMed Abstract | Crossref Full Text | Google Scholar

27. Ding C, Li G, Li Y, Gao H, Sun F. The construction and experimental verification of a 6-LncRNA model based on Lactic acid metabolism in the tumor microenvironment of Wilms tumor. Gene. (2025) 932:148898. doi: 10.1016/j.gene.2024.148898

PubMed Abstract | Crossref Full Text | Google Scholar

28. Guo Y, Lu W, Zhang Z, Liu H, Zhang A, Zhang T, et al. A novel pyroptosis-related gene signature exhibits distinct immune cells infiltration landscape in Wilms’ tumor. BMC Pediatr. (2024) 24:279. doi: 10.1186/s12887-024-04731-0

PubMed Abstract | Crossref Full Text | Google Scholar

29. Naito T, Saheki Y. GRAMD1-mediated accessible cholesterol sensing and transport. Biochim Biophys Acta Mol Cell Biol Lipids. (2021) 1866:158957. doi: 10.1016/j.bbalip.2021.158957

PubMed Abstract | Crossref Full Text | Google Scholar

30. Wu YW, Waldmann H. Toward the role of cholesterol and cholesterol transfer protein in autophagosome biogenesis. Autophagy. (2019) 15:2167–68. doi: 10.1080/15548627.2019.1666595

PubMed Abstract | Crossref Full Text | Google Scholar

31. Laraia L, Friese A, Corkery DP, Konstantinidis G, Erwin N, Hofer W, et al. The cholesterol transfer protein GRAMD1A regulates autophagosome biogenesis. Nat Chem Biol. (2019) 15:710–20. doi: 10.1038/s41589-019-0307-5

PubMed Abstract | Crossref Full Text | Google Scholar

32. Fu B, Meng W, Zhao H, Zhang B, Tang H, Zou Y, et al. GRAM domain-containing protein 1A (GRAMD1A) promotes the expansion of hepatocellular carcinoma stem cell and hepatocellular carcinoma growth through STAT5. Sci Rep. (2016) 6:31963. doi: 10.1038/srep31963

PubMed Abstract | Crossref Full Text | Google Scholar

33. Irons JL, Hodge-Hanson K, Downs DM. RidA proteins protect against metabolic damage by reactive intermediates. Microbiol Mol Biol Rev. (2020) 84:e00024-20. doi: 10.1128/MMBR.00024-20

PubMed Abstract | Crossref Full Text | Google Scholar

34. Ernst DC, Downs DM. 2-Aminoacrylate Stress Induces a Context-Dependent Glycine Requirement in ridA Strains of Salmonella enterica. J Bacteriol. (2016) 198:536–43. doi: 10.1128/JB.00804-15

PubMed Abstract | Crossref Full Text | Google Scholar

35. Irons J, Sacher JC, Szymanski CM, Downs DM. Cj1388 Is a RidA Homolog and Is Required for Flagella Biosynthesis and/or Function in Campylobacter jejuni. Front Microbiol. (2019) 10:2058. doi: 10.3389/fmicb.2019.02058

PubMed Abstract | Crossref Full Text | Google Scholar

36. Bussolati G, Geuna M, Bussolati B, Millesimo M, Botta M, Bartorelli A. Cytolytic and tumor inhibitory antibodies against UK114 protein in the sera of cancer patients. Int J Oncol. (1997) 10:779–85. doi: 10.3892/ijo.10.4.779

PubMed Abstract | Crossref Full Text | Google Scholar

37. Funaro A, Horenstein AL, Ghisolfi G, Bussolati B, Bartorelli A, Bussolati G. Identification of a 220-kDa membrane tumor-associated antigen by human anti-UK114 monoclonal antibodies selected from the immunoglobulin repertoire of a cancer patient. Exp Cell Res. (1999) 247:441–50. doi: 10.1006/excr.1998.4384

PubMed Abstract | Crossref Full Text | Google Scholar

38. Ceciliani F, Faotto L, Negri A, Colombo I, Berra B, Bartorelli A, et al. The primary structure of UK114 tumor antigen. FEBS Lett. (1996) 393:147–50. doi: 10.1016/0014-5793(96)00850-2

PubMed Abstract | Crossref Full Text | Google Scholar

39. Chong CL, Huang SF, Hu CP, Chen YL, Chou HY, Chau GY, et al. Decreased expression of UK114 is related to the differentiation status of human hepatocellular carcinoma. Cancer Epidemiol Biomarkers Prev. (2008) 17:535–42. doi: 10.1158/1055-9965.EPI-07-0506

PubMed Abstract | Crossref Full Text | Google Scholar

40. Fossat N, Tourle K, Radziewic T, Barratt K, Liebhold D, Studdert JB, et al. C to U RNA editing mediated by APOBEC1 requires RNA-binding protein RBM47. EMBO Rep. (2014) 15:903–10. doi: 10.15252/embr.201438450

PubMed Abstract | Crossref Full Text | Google Scholar

41. Wolfe AD, Arnold DB, Chen XS. Comparison of RNA editing activity of APOBEC1-A1CF and APOBEC1-RBM47 complexes reconstituted in HEK293T cells. J Mol Biol. (2019) 431:1506–17. doi: 10.1016/j.jmb.2019.02.025

PubMed Abstract | Crossref Full Text | Google Scholar

42. Cieply B, Park JW, Nakauka-Ddamba A, Bebee TW, Guo Y, Shang X, et al. Multiphasic and dynamic changes in alternative splicing during induction of pluripotency are coordinated by numerous RNA-binding proteins. Cell Rep. (2016) 15:247–55. doi: 10.1016/j.celrep.2016.03.025

PubMed Abstract | Crossref Full Text | Google Scholar

43. Shivalingappa P, Singh DK, Sharma V, Arora V, Shiras A, Bapat SA. RBM47 is a critical regulator of mouse embryonic stem cell differentiation. Stem Cell Rev Rep. (2023) 19:475–90. doi: 10.1007/s12015-022-10441-w

PubMed Abstract | Crossref Full Text | Google Scholar

44. Fossat N, Radziewic T, Jones V, Tourle K, Tam PP. Conditional restoration and inactivation of Rbm47 reveal its tissue-context requirement for viability and growth. Genesis. (2016) 54:115–22. doi: 10.1002/dvg.22920

PubMed Abstract | Crossref Full Text | Google Scholar

45. Shivalingappa P, Sharma V, Shiras A, Bapat SA. RNA binding motif 47 (RBM47): emerging roles in vertebrate development, RNA editing and cancer. Mol Cell Biochem. (2021) 476:4493–505. doi: 10.1007/s11010-021-04256-5

PubMed Abstract | Crossref Full Text | Google Scholar

46. Qin Y, Sun W, Wang Z, Dong W, He L, Zhang T, et al. RBM47/SNHG5/FOXO3 axis activates autophagy and inhibits cell proliferation in papillary thyroid carcinoma. Cell Death Dis. (2022) 13:270. doi: 10.1038/s41419-022-04728-6

PubMed Abstract | Crossref Full Text | Google Scholar

47. Guo T, You K, Chen X, Sun Y, Wu Y, Wu P, et al. RBM47 inhibits hepatocellular carcinoma progression by targeting UPF1 as a DNA/RNA regulator. Cell Death Discovery. (2022) 8:320. doi: 10.1038/s41420-022-01112-3

PubMed Abstract | Crossref Full Text | Google Scholar

48. Soleymanjahi S, Blanc V, Molitor EA, Alvarado DM, Xie Y, Gazit V, et al. RBM47 regulates intestinal injury and tumorigenesis by modifying proliferation, oxidative response, and inflammatory pathways. JCI Insight. (2023) 8(9):e161118. doi: 10.1172/jci.insight.161118

PubMed Abstract | Crossref Full Text | Google Scholar

49. Wu Y, Chen P, Sun L, Yuan S, Cheng Z, Lu L, et al. Sepiapterin reductase: Characteristics and role in diseases. J Cell Mol Med. (2020) 24:9495–506. doi: 10.1111/jcmm.15608

PubMed Abstract | Crossref Full Text | Google Scholar

50. Wu Y, Du H, Zhan M, Wang H, Chen P, Du D, et al. Sepiapterin reductase promotes hepatocellular carcinoma progression via FoxO3a/Bim signaling in a nonenzymatic manner. Cell Death Dis. (2020) 11:248. doi: 10.1038/s41419-020-2471-7

PubMed Abstract | Crossref Full Text | Google Scholar

51. Cheng C, Cai Y, Liu X, Wu Y, Cheng Q, Wu Y, et al. KHSRP modulated cell proliferation and cell cycle via regulating PPP2CA and p27 expression in Wilms tumor. Cell Signal. (2022) 100:110447. doi: 10.1016/j.cellsig.2022.110447

PubMed Abstract | Crossref Full Text | Google Scholar

52. Tang J, Wang F, Cheng G, Si S, Sun X, Han J, et al. Wilms’ tumor 1-associating protein promotes renal cell carcinoma proliferation by regulating CDK2 mRNA stability. J Exp Clin Cancer Res. (2018) 37:40. doi: 10.1186/s13046-018-0706-6

PubMed Abstract | Crossref Full Text | Google Scholar

53. Kirschner KM, Foryst-Ludwig A, Gohlke S, Li C, Flores RE, Kintscher U, et al. Wt1 haploinsufficiency induces browning of epididymal fat and alleviates metabolic dysfunction in mice on high-fat diet. Diabetologia. (2022) 65:528–40. doi: 10.1007/s00125-021-05621-1

PubMed Abstract | Crossref Full Text | Google Scholar

54. Ashton TM, McKenna WG, Kunz-Schughart LA, Higgins GS. Oxidative phosphorylation as an emerging target in cancer therapy. Clin Cancer Res. (2018) 24:2482–90. doi: 10.1158/1078-0432.CCR-17-3070

PubMed Abstract | Crossref Full Text | Google Scholar

55. Aminzadeh S, Vidali S, Sperl W, Kofler B, Feichtinger RG. Energy metabolism in neuroblastoma and Wilms tumor. Transl Pediatr. (2015) 4:20–32. doi: 10.3978/j.issn.2224-4336.2015.01.04

PubMed Abstract | Crossref Full Text | Google Scholar

56. Meng J, Chen Y, Lu X, Ge Q, Yang F, Bai S, et al. Macrophages and monocytes mediated activation of oxidative phosphorylation implicated the prognosis and clinical therapeutic strategy of Wilms tumor. Comput Struct Biotechnol J. (2022) 20:3399–408. doi: 10.1016/j.csbj.2022.06.052

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Wilms tumor, GRAMD1A, prognostic markers, LASSO regression, cell cycle

Citation: Zeng Q, Tao J, Qin L, Zeng Y, Liu Z, Xu M and Zeng L (2024) Comprehensive prognostic gene identification and functional characterization of GRAMD1A in Wilms tumor: development of risk prediction models and therapeutic implications. Front. Oncol. 14:1501718. doi: 10.3389/fonc.2024.1501718

Received: 25 September 2024; Accepted: 04 November 2024;
Published: 26 November 2024.

Edited by:

Dan Liu, Sichuan University, China

Reviewed by:

Yuxuan Song, Peking University People’s Hospital, China
Beichen Li, Institute of Science Tokyo, Japan
Alec Lee, University of Basel, Switzerland

Copyright © 2024 Zeng, Tao, Qin, Zeng, Liu, Xu and Zeng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Linshan Zeng, emxzZWQxMjE2QDE2My5jb20=; Mingxian Xu, eHVteDhAbWFpbDIuc3lzdS5lZHUuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.