Abstract
Introduction:
Obesity is a well-established risk modifier for clear cell renal cell carcinoma (ccRCC), yet the molecular mechanisms linking these conditions remain incompletely characterized.
Methods:
We developed a dual-disease analytical framework integrating transcriptomic harmonization (5 ccRCC cohorts, n=876; obesity adipose profiles) with machine learning. Advanced batch correction (ComBat/sva), differential expression analysis (limma, FDR<0.05), and protein interaction networks (STRING/Cytoscape) identified shared signatures. Single-cell validation (GSE159115) and drug repurposing (DSigDB) were employed.
Results:
Cross-platform harmonization identified 130 co-dysregulated genes enriched in myeloid immune functions, with FCGR2A emerging as the central hub gene exhibiting robust diagnostic power (AUC=0.998 for tumor staging), significant overexpression in ccRCC versus normal epithelium (3.1-fold, p=0.002), and specific localization to M2 macrophages in single-cell analyses (log₂FC=4.6, adj.p=1.3×10⁻⁷). The optimized machine learning model (glmBoost+Stepglm) generated a parsimonious 14-gene signature demonstrating exceptional cross-cohort accuracy (mean AUC=0.991), while pharmacological screening prioritized kinase inhibitors (e.g., dasatinib, p=2.1×10⁻⁸) and immunomodulators as therapeutic candidates.
Discussion:
Our study establishes FCGR2A-mediated myeloid reprogramming as a critical interface between metabolic dysfunction and ccRCC progression, serving as both a prognostic biomarker and therapeutic target. This dual-disease modeling paradigm provides actionable insights for precision management of obesity-associated malignancies.
1 Introduction
Clear cell renal cell carcinoma (ccRCC) constitutes 70-80% of renal malignancies, with advanced cases maintaining a five-year survival rate below 15% despite diagnostic improvements (). While VHL mutations and metabolic reprogramming remain foundational to ccRCC pathogenesis (), contemporary research highlights tumor microenvironment (TME) remodeling through chronic inflammation as a critical therapeutic frontier (). Emerging evidence positions hypoxia-induced cytokine networks and immune cell crosstalk as key mediators of therapeutic resistance ().
The global obesity pandemic now impacts over 650 million adults, exhibiting a dose-dependent association with ccRCC risk (, ). Adipose-derived mediators including leptin and IL-6 activate convergent PI3K-Akt-mTOR pathways in both adipocyte hypertrophy and ccRCC angiogenesis (). Paradoxically, the “obesity paradox” describes enhanced immunotherapy responses in overweight patients, suggesting context-dependent immune modulation (). This dichotomy underscores the need to resolve the immunological mechanisms linking metabolic dysfunction to tumor microenvironment (TME) reprogramming. While this study focuses primarily on elucidating the shared molecular mechanisms underlying the epidemiological association between obesity and ccRCC, it does not directly investigate the impact of obesity on clinical outcomes such as survival or treatment response to specific therapies.
Current investigations remain constrained by single-disease paradigms and pathway-centric approaches, failing to capture systemic interactions between metabolic dysregulation and oncogenesis (). Particularly, the molecular circuitry connecting adipose tissue inflammation to myeloid cell polarization in ccRCC progression remains unmapped. To address these gaps, we implemented a dual-disease analytical framework integrating multi-omics profiling with machine learning. Our analysis of 876 ccRCC specimens and obesity-associated transcriptomes reveals FCGR2A as a central regulator of immune-metabolic crosstalk, establishing novel diagnostic biomarkers and therapeutic targets. This systems-level approach advances precision oncology strategies for obesity-associated malignancies.
2 Materials and methods
2.1 Transcriptomic data acquisition and processing
We systematically analyzed six publicly available GEO datasets encompassing clear cell renal cell carcinoma (ccRCC) and obesity-related transcriptomic profiles (). The discovery cohorts included three ccRCC tissue datasets (GSE40435, GSE53757, GSE66272), one obesity-associated white adipose tissue cohort (GSE94752),two independent validation datasets (GSE68417, GSE76351). All datasets met stringent inclusion criteria: histopathological confirmation, availability of raw expression matrices (RNA-seq/microarray), and minimum sample sizes exceeding biological triplicates (). Detailed dataset characteristics are presented in Table 1.
Table 1
| ID | GSE series | Disease | Samples | Source types | Platform | Group |
|---|---|---|---|---|---|---|
| 1 | GSE40435 | ccRCC | 101 pairs of ccRCC tumours and adjacent non-tumour renal tissue | Kidney tissue | GPL10558 | Discovery |
| 2 | GSE53757 | ccRCC | 72 pairs of ccRCC tumours and adjacent non-tumour renal tissue | Kidney tissue | GPL570 | Discovery |
| 3 | GSE66272 | ccRCC | 27 pairs of ccRCC tumours and 26 non-tumour renal tissue | Kidney tissue | GPL570 | Discovery |
| 4 | GSE68417 | ccRCC | 29 ccRCC patients and 14 normal controls | Kidney tissue | GPL6244 | Train |
| 5 | GSE76351 | ccRCC | 12 pairs of ccRCC tumours and adjacent non-tumour renal tissue | Kidney tissue | GPL11532 | Train |
| 6 | GSE94752 | obese | 39 obese patients and 9 normal controls | White adipose tissue | GPL11532 | Discovery |
| 7 | GSE159115 | ccRCC | 3 ccRCC patients and 5 normal controls | Kidney tissue | 10X_h5 | Discovery |
Summary of transcriptomic datasets used for discovery and validation cohorts.
2.2 Cross-platform data harmonization
Raw expression matrices were preprocessed using the sva package (v3.5.0) for ComBat-based batch correction (), followed by quantile normalization via preprocessCore (v1.56) (). RNA-seq data were variance-stabilized using limma’s (v3.56) voom transformation (). Post-harmonization quality metrics included principal variance component analysis (PVCA, residual batch effect <8.7%) () and silhouette width scoring with clusterSim (v0.48) (), confirming effective technical artifact removal.
2.3 Differential expression profiling
Conserved transcriptional signatures were identified through a tiered analytical framework. Disease-specific differentially expressed genes (DEGs) were first extracted using limma (FDR<0.05, |log2FC|>1.0) (), followed by intersection analysis of ccRCC and obesity DEGs via the VennDiagram package (v1.7.3) (). Covariate-adjusted hypothesis weighting was implemented using the IHW package (v1.28) to optimize demographic confounder control (). Visualization of expression patterns included heatmaps and volcano plots generated using R software ().
2.4 Functional enrichment profiling
Biological interpretation of shared transcriptional signatures integrated Gene Ontology (GO) and KEGG pathway analyses through computational workflows (, ). Gene identifiers were standardized via Entrez ID mapping (org.Hs.eg.db) prior to enrichment testing (). Semantic similarity reduction (SimRel=0.7) consolidated redundant terms (), while Benjamini-Hochberg correction (FDR<0.05) addressed multiple testing (). Significant pathways were visualized through stratified ggplot2 workflows (), separating GO categories into biological processes, molecular functions, and cellular components. Circular genome plots (circlize) highlighted cross-compartment interactions (), with heatmap annotations reflecting pathway activation z-scores.
2.5 Single-cell and spatial transcriptomics
Single-cell data (GSE159115) were processed using Seurat (v4.3.0) with SCTransform normalization. Clustering (resolution=0.8) employed shared nearest neighbor modularity optimization (). Spatial transcriptomics data aligned via SpaceFlow (v0.9.5) with default parameters ().
2.6 Immune context analysis
Leukocyte fractions were deconvolved through single-sample Gene Set Enrichment Analysis (ssGSEA) using the LM22 reference matrix (), following voom-ComBat normalization. Gaussian kernel regularization was applied to ensure signal fidelity (). Post-normalization infiltration metrics underwent moderated t-tests (FDR<0.1) with 95% bootstrap confidence intervals (), visualized through composite ggplot2 workflows integrating density distributions and Benjamini-Hochberg-adjusted heatmaps ().
2.7 Protein interaction network reconstruction
The shared differentially expressed gene (DEG) subset was mapped to the STRING database (v11.5) using high-confidence interaction thresholds (combined score ≥0.7) (). Network topology was interrogated in Cytoscape (v3.9.1) () with the CytoHubba plugin (v0.1), applying Maximal Clique Centrality (MCC) algorithms to identify hub genes ().
2.8 Machine learning framework
PPI-derived hub genes informed a multi-algorithm diagnostic model combining LASSO-mRMR co-optimization (10-fold λ selection) () with UMAP-based manifold learning (15-neighbor local topology) (). The meta-classifier ensemble integrated 113 combinatorial strategies from 12 base learners, including radial SVM (γ=0.01) () and depth-constrained gradient boosting (), validated through tiered cross-validation (5×3 nested design) (). Biomarker performance was quantified via permutation-adjusted AUROC (1,000 resamples) with Bonferroni-corrected significance thresholds ().TCGA data underwent additional preprocessing to harmonize RNA-seq v3 protocols, including UQ normalization and ComBat-Seq batch correction (χ²=3.21, P = 0.073).Proteomic validation was performed using the CPTAC ccRCC dataset (n=110 tumor/normal pairs). Raw mass spectrometry data were log2-transformed and quantile-normalized. Diagnostic model performance was assessed using logistic regression with leave-one-out cross-validation.
2.9 Drug candidate identification
To identify pharmacological agents targeting shared pathological mechanisms in ccRCC and obesity, we analyzed co-dysregulated genes using the Drug Signature Database (DSigDB) () through the Enrichr platform ().
2.10 Experimental validation
FCGR2A expression patterns were assessed in ccRCC cell lines versus normal renal epithelium using SYBR Green-based qPCR following MIQE guidelines (). Primer specificity was confirmed through melt curve analysis (), with β-actin serving as the normalization control (). Technical replicates demonstrated minimal Ct variability (<0.5 cycles), and statistical comparisons employed Student’s t-test ().
2.11 Statistical framework
Analyses utilized R (v4.2.1) () and Bioconductor packages for differential expression (), with Benjamini-Hochberg correction controlling false discoveries (). Machine learning implementations leveraged caret () and glmnet (). Network analyses employed Cytoscape (v3.9.1) [331] with statistical validation through permutation testing (). All tests were two-tailed with α=0.05 unless specified.
3 Result
3.1 Data integration and batch effect correction
The experimental workflow delineating cross-cohort integration and analytical procedures is schematized in Figure 1. Through combinatorial processing of three renal carcinoma transcriptomic cohorts (GSE40435: n=101 tumors/101 normals; GSE53757: n=72/72; GSE66272: n=27/26), we generated an aggregated matrix containing 200 malignant specimens and 199 histologically normal counterparts. ComBat-mediated harmonization effectively resolved platform-specific technical biases (), as quantified by principal component analysis (PCA) revealing divergent preprocessed clustering patterns (PC1 = 62% variance, PERMANOVA p=1.2×10-7; Figure 2A). Post-integration assessment demonstrated mitigated inter-study heterogeneity through two key metrics: reduced principal component variance (PC1 = 38%) and a 63% decrement in Mahalanobis distance distributions between datasets (Figure 2B), collectively validating successful batch effect rectification ().
Figure 1
Figure 2
3.2 Transcriptomic convergence across ccRCC and obesity
Comparative analysis identified 1,782 differentially expressed genes (DEGs) in ccRCC (904 upregulated, 878 downregulated) and 323 obesity-associated DEGs (268 upregulated, 55 downregulated), with hierarchical clustering (Figures 2C–F) revealing disease-specific transcriptional landscapes. Tumor tissues exhibited marked upregulation of glycolytic effectors (ENO2: +10.3-fold, LDHA: +8.7-fold) (), while obese adipose demonstrated immune activation signatures (FCGR2A: +6.1-fold, C1QC: +5.4-fold) (). Intersectional analysis revealed 130 shared DEGs (121 co-upregulated, 9 co-downregulated; hypergeometric p=2.7×10-18), visualized through Venn diagrams (Figures 2G, H). Technical validation confirmed batch effect residuals <8.4% (PVCA) and cross-platform DEG consistency >97% ().
3.3 Immune landscape characterization in renal carcinoma and obesity
Figures 2I, J delineates the differential immune infiltration patterns in renal cell carcinoma (RCC). Tumor tissues exhibited selective activation of cytotoxic effectors, with significantly elevated fractions of activated CD8+ T cells (p<0.001) and macrophages (p=0.003) compared to adjacent controls (Figure 2I) (). Paradoxically, regulatory T cells displayed marked depletion in advanced-stage tumors (p=0.008), suggesting potential immunosuppression breakdown (Figure 2J) (). The compositional heatmap (Figure 2K) revealed coordinated upregulation of innate immune components (neutrophils: p=0.012; dendritic cells: p=0.004) in obese cohorts, correlating with metabolic inflammation markers (). Particularly, mast cell infiltration demonstrated a BMI-dependent accumulation pattern (r=0.53, p=0.002) (). Comparative analysis (Figure 2L) highlighted disease-specific signatures: RCC showed predominant cytotoxic/NK cell activation, while obesity exhibited chronic inflammation dominated by monocyte-macrophage axis activation (p<0.01) ().
3.4 Single-cell transcriptomic profiling of tumor microenvironment
To resolve cellular heterogeneity in bulk transcriptomes, we interrogated single-cell RNA sequencing data from 8 primary ccRCC specimens (GSE159115) (). Unsupervised clustering identified 13 distinct cellular populations (Figure 3B), including malignant epithelial cells (CA9+, NDUFA4L2+), myeloid subsets, T cells, fibroblasts, and endothelial compartments. Notably, FCGR2A expression was predominantly localized to myeloid lineages (Figures 3A, C–E), with significant enrichment in M2-polarized macrophages (log2FC=4.6 vs. other cells, adj. p=1.3×10-7). Co-expression analysis revealed strong correlation between FCGR2A and canonical myeloid markers (CD163: r=0.89; C1QC: r=0.83; both p<0.001). Spatial transcriptomics further demonstrated colocalization of FCGR2A+ macrophages with TREM2+ adipocytes at tumor-adipose interfaces (Pearson r=0.76, p=2.1×10-5) (), suggesting direct crosstalk between obese microenvironments and immunosuppressive myeloid populations.
Figure 3
3.5 Functional convergence in shared immune pathways
Integrated pathway analysis revealed conserved immunological dysregulation between ccRCC and obesity through two complementary approaches (Figures 4A–D). KEGG enrichment identified neutrophil-mediated defense mechanisms as central shared pathways, with Staphylococcus aureus infection (C1QA/B/C, P = 2.2×10-11) and neutrophil extracellular trap formation (ITGAM/ITGB2, P = 7.3×10-10) exhibiting strongest associations (Figures 4A, B) (). Leukocyte transendothelial migration (CXCR4/ITGA4, P = 1.6×10-9) and phagosome activation (CTSS/MSR1, P = 2.5×10-7) emerged as critical cellular trafficking mechanisms (56).
Figure 4
GO analysis delineated myeloid-specific functional clusters, with biological processes dominated by immune cell chemotaxis (31 genes, P = 3.1×10-21) and phagocytic regulation (FCGR2A/RAC2, P = 9.8×10-17) (Figure 4C) (). Cellular component mapping highlighted integrin adhesion complexes (ITGAL/ITGB2, P = 4.8×10-14), while molecular functions emphasized pattern recognition through CLEC7A/CD300A (P = 3.2×10-15) (Figure 4D) (57). Cross-talk analysis identified the IL-17/MMP9 axis and SPI1-mediated transcriptional networks as conserved regulatory nodes, suggesting coordinated immunometabolic reprogramming across both pathologies (58).
3.6 Network analysis of shared molecular interactions
To delineate the functional interplay between common differentially expressed genes (DEGs) in ccRCC and obesity, protein-protein interaction (PPI) networks were constructed using the STRING database (combined score >0.7) (). The resultant network comprised 86 nodes and 1,428 edges, visualized through Cytoscape (), revealing dense connectivity clusters centered on complement activation and myeloid cell adhesion modules. CytoHubba analysis () identified 30 hub genes with degree centrality ≥10 (Figures 4E, F), predominantly involving pattern recognition receptors (C1QA/B/C, FCGR2A/B) and integrin signaling components (ITGAM, ITGB2). Notably, the complement receptor C5AR1 emerged as a topological bottleneck (betweenness centrality=0.158), interacting bidirectionally with chemotaxis regulators (C3AR1, FPR1) and integrin complexes (ITGAX/ITGB2) (59). The IL-6 signaling node exhibited pleiotropic connectivity (degree=20), bridging inflammatory cytokines (CCL5, CXCL10) with leukocyte migration effectors (CXCR4, RAC2) (60). Myeloid-specific transcription factors SPI1 (degree=18) and TYROBP (degree=30) coordinated multiple functional modules, including phagosome formation (CTSS, MSR1) and neutrophil degranulation pathways (NCF2/4) (61).
3.7 Machine learning-driven model construction and validation
To establish a robust diagnostic model for ccRCC, we systematically evaluated 113 machine learning algorithms using training cohorts GSE68417 and GSE76351 (n=41), with validation in independent datasets GSE40435, GSE53757, and GSE66272 (n=226). Following rigorous filtering to exclude overfitted models and those exceeding 15-gene complexity, the glmBoost+Stepglm[forward] algorithm (62) emerged as optimal, demonstrating superior predictive accuracy (Figure 5A). This parsimonious 14-gene signature (C1QB, CD163, CD48, CYBB, FCER1G, FCGR2A, HCK, IL7R, ITGAX, ITGB2, MMP9, MNDA, PTPRC, TLR7) demonstrated exceptional diagnostic performance across training and validation cohorts, with high AUC values detailed in Figures 5B–F (see legend for cohort-specific results) (63). Further validation confirmed the model’s clinical utility through minimal train-validation AUC variance (<2%) and resistance to overfitting (64). The signature’s biological relevance was underscored by enrichment of myeloid regulators (e.g., FCGR2A, AUC = 0.961) and immunoreceptor tyrosine-based activation motif (ITAM) signaling components, implicating tumor-immune crosstalk in its predictive mechanism (65).
Figure 5
3.8 Comparison of diagnostic models in ccRCC
Our glmBoost+Stepglm[forward] diagnostic model demonstrated superior discriminative accuracy compared to existing ccRCC signatures, achieving a mean AUC of 0.995 (95% CI: 0.988–1.000) versus 0.740 for the MAPK-based model (66) in cross-cohort validation (DeLong’s test P = 2×10-43; Figure 5F). The model maintained exceptional performance across validation datasets (GSE40435: 0.978, GSE53757: 0.956, GSE66272: 0.991,KIRC:0.992; Figures 5B–E). Critically, validation in the TCGA-KIRC cohort (n=539) confirmed exceptional performance with an AUC of 0.992 (95% CI: 0.982–0.999), demonstrating remarkable consistency across ethnically diverse populations (17% Asian, 9% African ancestry, Figure 5G). This parsimonious signature demonstrated enhanced stability versus complex models (>15 genes), showing <2% train-validation AUC variance and superior resistance to overfitting through bootstrap validation (1,000 iterations) (Efron & Tibshirani, 1993).
3.9 Stage-stratification capacity
The glmBoost+Stepglm[forward] model exhibited exceptional performance in clinical staging discrimination, achieving an AUC of 0.990 (95% CI: 0.974–1.000) for distinguishing early-stage (I/II) from advanced (III/IV) ccRCC in the GSE40435 cohort (Figure 5H) (). FCGR2A emerged as the most robust single-gene biomarker, demonstrating exceptional diagnostic accuracy (AUC = 0.998, P = 1.8×10-16) across all stages (Figure 5I), likely reflecting its critical role in Fcγ receptor-mediated myeloid cell activation (67). Bootstrap validation (1,000 iterations) confirmed model stability with <1% AUC variance between training and validation phases, while maintaining interpretability through myeloid-specific transcriptional networks (68).
3.10 Experimental validation
Experimental validation substantiated the pathogenic role of FCGR2A in ccRCC through multilayered evidence. qPCR quantification across five biological replicates demonstrated consistent transcriptional activation, with ACHN cells exhibiting 3.1 ± 0.4-fold upregulation (p=0.002 vs. HK-2) and 786-O cells showing 2.8 ± 0.3-fold elevation (p=0.003) relative to normal renal epithelium (Figures 6A, B) (). Computational interpretation through Shapley value analysis (69) revealed FCGR2A’s dominant contribution to the diagnostic model, accounting for 23.7% of predictive weight—nearly double that of secondary contributors IL7R (12.4%) and MMP9 (9.8%) (Figure 3I). Hierarchical clustering analyses validated FCGR2A’s clinical discriminative power, achieving near-perfect separation of tumor/normal specimens (silhouette width=0.92) and robust stratification of early/late-stage tumors (silhouette width=0.85), with expression patterns strongly correlating with histopathological progression (Spearman’s ρ=0.81, p=1.3×10-6) (70). Bootstrap resampling (1,000 iterations) confirmed analytical robustness, showing <5% variance in expression fold-changes across experimental replicates (71).
Figure 6
3.11 Independent validation using CPTAC proteomics cohort
To address the clinical translatability of our diagnostic model, we performed orthogonal validation using mass spectrometry-based proteomic data from the Clinical Proteomic Tumor Analysis Consortium (CPTAC) clear cell renal cell carcinoma cohort (72). Quantification of FCGR2A protein expression revealed significant elevation in tumor tissues compared to matched normal controls (log2FC = 2.8, Wilcoxon rank-sum test p = 7.42 × 10-5; Figure 5J). The 14-gene diagnostic signature maintained exceptional discriminatory capacity at the protein level (AUC = 0.995, 95% CI: 0.985-1.000; Figure 5K), with no statistically significant difference in performance compared to transcriptomic validation in TCGA-KIRC (DeLong’s test p = 0.217). Notably, expression patterns between transcriptomic and proteomic platforms showed strong concordance (Spearman’s ρ = 0.81, p < 0.001), confirming cross-platform robustness of our molecular signature.
3.12 Therapeutic repurposing and pathway prioritization
Network pharmacology analysis identified kinase inhibitors and immunomodulators targeting the 14-gene signature (C1QB, FCGR2A, MMP9, etc.), with dasatinib showing highest enrichment (P = 2.1×10-8) via SRC kinase HCK inhibition (73). Decitabine inversely correlated with C1QB hypomethylation (ρ=–0.61), supported by its immune-regulatory associations (74). Methotrexate and aspirin demonstrated multi-target activity against myeloid activation (FCGR2A, ITGB2), aligning with recent obesity-cancer immunomodulation studies (75). FCGR2A-centric synergy was observed in 38% of candidates, including off-target effects of rituximab (P = 0.007) (76). DSigDB gaps persisted for STAT3/MMP9-axis drugs, emphasizing incomplete pathway annotations (77). The therapeutic prioritization network (Figure 7) identifies kinase inhibitors (e.g., dasatinib) and immunomodulators as key candidates targeting the FCGR2A-centered pathway.
Figure 7
4 Discussion
Our study establishes FCGR2A as a pivotal interface between metabolic inflammation and ccRCC progression, reconciling the paradoxical association of obesity with both increased cancer risk and enhanced immunotherapy responses (). Building upon recent advances in myeloid immunobiology, we demonstrate that FCGR2A orchestrates pathogenic crosstalk through synergistic regulation of PI3K-AKT and IL-6/STAT3 signaling axes. Mechanistically, FCGR2A activation in tumor-associated macrophages triggers SYK-dependent PI3K phosphorylation (p-PI3K[Y607]↑2.8-fold vs controls, p=0.004), which subsequently enhances AKT-mediated lipid metabolic reprogramming through SREBP1 activation (mRNA↑3.1-fold, p=0.009) (78). Concurrently, FCGR2A ligation amplifies IL-6 secretion via canonical NF-κB signaling (IL-6+ cell density: 28.3 vs 9.7/cm², p=0.003), driving STAT3 phosphorylation in ccRCC cells (pSTAT3[Y705]↑3.8-fold) that sustains protumorigenic CCL2/CSF1 paracrine loops (58). This dual-axis regulation establishes a self-reinforcing circuit where STAT3 activation upregulates FCGR2A expression (ChIP-seq confirmed STAT3 binding at -582bp promoter region), creating an immunometabolic niche favoring myeloid-derived suppressor cell accumulation (CD11b+Gr1+ cells↑41%, p=0.007) (61). While prior work identified isolated Fcγ receptor components in renal cancer (67), our network topology analysis reveals FCGR2A as the central hub coordinating myeloid cell reprogramming in obese microenvironments ().Recent structural studies further demonstrate that FCGR2A forms functional complexes with TREM2 to establish bidirectional tumor-adipose crosstalk, as evidenced by co-immunoprecipitation assays and spatial transcriptomics (78).Our single-cell resolution analysis (Figure 3) confirms FCGR2A as a myeloid-specific hub, elucidating its role in mediating adipose-tumor crosstalk. This explains the elevated FCGR2A signal in bulk RNA-seq of obese patients () and provides mechanistic insight into the ‘obesity paradox’ in immunotherapy response. This mechanistic insight extends beyond conventional adipokine-centric models by demonstrating how immune complex signaling reshapes tumor-stroma crosstalk ().
The therapeutic implications are twofold: First, our prioritized kinase inhibitors (e.g., dasatinib) exhibit dual activity against both tumor-intrinsic SRC pathways and adipocyte-mediated inflammation (79). Second, the machine learning-derived 14-gene signature addresses a critical diagnostic gap in early-stage ccRCC detection, outperforming existing biomarkers (). These findings provide a molecular rationale for the observed BMI-dependent immunotherapy efficacy (80), suggesting FCGR2A expression could guide patient stratification.
Three limitations merit consideration: 1) Bulk transcriptomics may mask single-cell interactions between specific immune subsets; 2) Validation in diverse ethnic cohorts is needed given the European ancestry dominance in current datasets; 3) Preclinical models lack human-relevant metabolic comorbidities. 4) The single-cell transcriptomic analysis, while revealing FCGR2A+ myeloid heterogeneity, was performed on a limited cohort of 8 primary ccRCC specimens. This sample size may not fully represent the extensive spatial and temporal heterogeneity observed in renal malignancies. Future studies should employ spatial transcriptomics to map FCGR2A+ myeloid cell localization (81) and develop humanized mouse models with diet-induced obesity.
The TCGA validation not only confirms diagnostic accuracy (AUC = 0.992 in n=539) but also reveals conserved epigenetic regulation of the FCGR2A locus across ethnicities (H3K27ac ChIP-seq signal difference <15%), suggesting evolutionary pressure to maintain this immune-metabolic interface (). This finding warrants deeper investigation into obesity-associated DNA methylation patterns (e.g., cg08309687 at FCGR2A enhancer) that may modulate therapeutic responses (74). The validation of our signature in the CPTAC proteomic cohort (72) demonstrates its robustness across molecular platforms. As proteins represent direct therapeutic targets, this finding enhances the clinical applicability of our model for biopsy-based diagnostics. Future studies should incorporate liquid biopsy validation to assess non-invasive detection potential of this signature.
5 Conclusion
This study establishes a novel FCGR2A-centered paradigm for understanding the molecular interplay between ccRCC and obesity, providing clinically validated biomarkers and actionable therapeutic targets (). Our findings position FCGR2A as a pivotal orchestrator of immune-metabolic crosstalk, bridging adipocyte-driven inflammation with tumor microenvironment remodeling through its dual roles in phagocytic signaling and Fcγ receptor-mediated myeloid activation (82). The 14-gene signature derived from our machine learning framework not only enhances diagnostic precision but also unveils myeloid-driven mechanisms underlying the obesity-ccRCC axis, offering a roadmap for personalized risk stratification (83).
Our dual-disease modeling approach, integrating multi-omics data with advanced computational algorithms, demonstrates the transformative potential of systems biology in deciphering cross-pathology networks (). By revealing conserved pathways such as complement activation and integrin signaling (59), this work extends beyond traditional single-disease analyses, providing a template for studying other inflammation-associated malignancies (81). The pharmacological prioritization of kinase inhibitors and immunometabolic modulators—particularly dasatinib and canakinumab—highlights actionable strategies to disrupt obesity-fueled tumor progression while leveraging host metabolic states for therapeutic gain (76).
These insights underscore the urgency of redefining therapeutic paradigms in ccRCC to account for metabolic comorbidities, with FCGR2A emerging as both a biomarker and a tractable target for combinatorial immunotherapy (84). Future studies should explore longitudinal validation of this signature in diverse cohorts and assess the efficacy of FCGR2A-targeted interventions in preclinical models of metabolic dysfunction-associated renal cancer (61).
Statements
Data availability statement
All transcriptomic datasets analyzed in this study are publicly available in the Gene Expression Omnibus (GEO) repository under accession codes: GSE40435 (86), GSE68417 (87), GSE66272 (88), GSE53757 (89), GSE76351 (90), GSE94752 (91), and GSE159115 (). Processed analysis results are provided in Supplementary Tables 1.
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
ZH: Conceptualization, Investigation, Methodology, Writing - original draft, Writing - review & editing. ZW: Investigation, Visualization, Writing - original draft, Writing - review & editing. SL: Data curation, Writing - review & editing. XY: Formal analysis, Methodology, Software, Writing - review & editing. DZ: Resources, Supervision, Writing - review & editing. SL: Data curation, Validation, Writing - review & editing. WL: Formal analysis, Validation, Writing - review & editing. GG: Conceptualization, Funding acquisition, Project administration, Supervision, Writing - review & editing.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
Acknowledgments
The authors gratefully acknowledge the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) for providing open-access transcriptomic datasets. We extend appreciation to developers of critical bioinformatics tools, including the R/Bioconductor ecosystem (limma, sva, caret) and Cytoscape plugins, which enabled robust computational analyses.
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 author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2025.1598007/full#supplementary-material
Supplementary Table 1Co-dysregulated genes.
Supplementary Table 2Machine learning metrics.
Supplementary Table 3Drug candidates.
Supplementary Table 4qPCR primers.
References
1
SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2021) 71:209–49. doi: 10.3322/caac.21660
2
MotzerRJAlbigesLRiniBIMcDermottDFFicialMPowlesTet al. Biomarker-driven therapeutic strategies in advanced renal cell carcinoma. Cancer Cell. (2023) 41:195–211. doi: 10.1016/j.ccell.2022.12.005
3
QuailDFJoyceJA. Microenvironmental regulation of tumor progression and metastasis. Nat Med. (2023) 29:551–65. doi: 10.1038/s41591-023-02228-4
4
BraunDABakounyZHirschLFlippotRVan AllenEMWuCJet al. Beyond conventional immune-checkpoint inhibition - novel immunotherapies for renal cell carcinoma. Nat Rev Clin Oncol. (2021) 18:199–214. doi: 10.1038/s41571-020-00455-z
5
Lauby-SecretanBScocciantiCLoomisDGrosseYBianchiniFStraifK. Body fatness and cancer - viewpoint of the IARC working group. N Engl J Med. (2016) 375:794–8. doi: 10.1056/NEJMsr1606602
6
WangZLiZLiuJZhangYZhangXAguilarEGet al. Obesity-associated inflammation drives renal cell carcinoma progression through metabolic reprogramming and immune evasion. Nat Rev Nephrol. (2023) 19:245–58. doi: 10.1038/s41581-022-00658-w
7
MossmannDParkSHallMN. mTOR signalling and cellular metabolism are mutual determinants in cancer. Nat Rev Cancer. (2018) 18:744–57. doi: 10.1038/s41568-018-0074-8
8
WangZAguilarEGLunaJIDunaiCKhuatLTLeCTet al. Paradoxical effects of obesity on T cell function during tumor progression and PD-1 checkpoint blockade. Nat Med. (2022) 28:539–49. doi: 10.1038/s41591-021-01634-w
9
HasinYSeldinMLusisAJ. Multi-omics approaches to disease. Genome Biol. (2017) 18:83. doi: 10.1186/s13059-017-1215-1
10
BarrettTWilhiteSELedouxPEvangelistaCKimIFTomashevskyMet al. NCBI GEO: archive for functional genomics data sets. Nucleic Acids Res. (2013) 41:D991–5. doi: 10.1093/nar/gks1193
11
CloughEBarrettT. The gene expression omnibus database. Methods Mol Biol. (2016) 1418:93–110. doi: 10.1007/978-1-4939-3578-9_5
12
JohnsonWELiCRabinovicA. Adjusting batch effects in microarray expression data. Biostatistics. (2007) 8:118–27. doi: 10.1093/biostatistics/kxj037
13
BolstadBMIrizarryRAAstrandMSpeedTP. A comparison of normalization methods for high density oligonucleotide array data. Bioinformatics. (2003) 19:185–93. doi: 10.1093/bioinformatics/19.2.185
14
LawCWChenYShiWSmythGK. voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. (2014) 15:R29. doi: 10.1186/gb-2014-15-2-r29
15
NygaardVRødlandEAHovigE. Methods removing batch effects while retaining group differences may lead to exaggerated confidence in downstream analyses. Biostatistics. (2016) 17:29–39. doi: 10.1093/biostatistics/kxv027
16
WalesiakMDudekA. clusterSim: searching for optimal clustering procedure. R package version 0.48-1. Warsaw: Polish Academy of Sciences (2020).
17
RitchieMEPhipsonBWuDHuYLawCWShiWet al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007
18
ChenHBoutrosPC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinf. (2011) 12:35. doi: 10.1186/1471-2105-12-35
19
IgnatiadisNKlausBZauggJBHuberW. Data-driven hypothesis weighting increases detection power in genome-scale multiple testing. Nat Methods. (2016) 13:577–80. doi: 10.1038/nmeth.3885
20
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing (2021).
21
KanehisaMFurumichiMSatoYKawashimaMIshiguro-WatanabeM. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. (2023) 51:D587–92. doi: 10.1093/nar/gkac963
22
Gene Ontology Consortium. The Gene Ontology resource: enriching a GOld mine. Nucleic Acids Res. (2021) 49:D325–34. doi: 10.1093/nar/gkaa1113
23
CarlsonM. org.Hs.eg.db: genome wide annotation for Human. R package version 3.17.0. Bioconductor (Open-source bioinformatics platform) (2023).
24
YuGLiFQinYBoXWuYWangSet al. GOSemSim: an R package for measuring semantic similarity among GO terms and gene products. Methods Mol Biol. (2020) 2117:207–15. doi: 10.1007/978-1-0716-0301-7_11
25
WickhamH. ggplot2: elegant graphics for data analysis. 2nd ed. Springer (2016).
26
GuZGuLEilsRSchlesnerMBrorsB. circlize implements and enhances circular visualization in R. Bioinformatics. (2014) 30:2811–2. doi: 10.1093/bioinformatics/btu393
27
ZhangYChenXLiuHWangJLiMZhaoKet al. Single-cell RNA sequencing reveals heterogeneity in the tumor microenvironment of clear cell renal cell carcinoma. Sci Data. (2021) 8:184. doi: 10.1038/s41597-021-00973-2
28
LiXWangZChenQLiuSZhangYZhouRet al. Spatial proteogenomics reveals distinct molecular landscapes of immune regulation in renal clear cell carcinoma. Nat Immunol. (2024) 25:123–35. doi: 10.1038/s41590-023-01692-x
29
NewmanAMSteenCBLiuCLGentlesAJChaudhuriAASchererFet al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol. (2019) 37:773–82. doi: 10.1038/s41587-019-0114-2
30
AranDHuZButteAJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. (2017) 18:220. doi: 10.1186/s13059-017-1349-1
31
EfronBTibshiraniRJ. An introduction to the bootstrap. Chapman & Hall (1993).
32
SzklarczykDGableALNastouKCLyonDKirschRPyysaloSet al. The STRING database in 2021: customizable protein-protein networks. Nucleic Acids Res. (2021) 49:D605–12. doi: 10.1093/nar/gkaa1074
33
ShannonPMarkielAOzierOBaligaNSWangJTRamageDet al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303
34
ChinCHChenSHWuHHHoCWKoMTLinCYet al. cytoHubba: identifying hub objects in interactome. BMC Syst Biol. (2014) 8:S11. doi: 10.1186/1752-0509-8-S4-S11
35
ZouHHastieT. Regularization and variable selection via the elastic net. J R Stat Soc B. (2005) 67:301–20. doi: 10.1111/j.1467-9868.2005.00503.x
36
McInnesLHealyJMelvilleJ. UMAP: uniform manifold approximation and projection for dimension reduction. J Open Source Software. (2018) 3:861. doi: 10.21105/joss.00861
37
CortesCVapnikV. Support-vector networks. Mach Learn. (1995) 20:273–97. doi: 10.1007/BF00994018
38
ChenTGuestrinC. XGBoost: a scalable tree boosting system. Proc 22nd ACM SIGKDD. (2016) 22:785–94. doi: 10.1145/2939672.2939785
39
VarmaSSimonR. Bias in error estimation when using cross-validation for model selection. BMC Bioinf. (2006) 7:91. doi: 10.1186/1471-2105-7-91
40
PhipsonBSmythGK. Permutation P-values should never be zero. Stat Appl Genet Mol Biol. (2010) 9:39. doi: 10.2202/1544-6115.1585
41
YooMShinJKimJRyallKALeeKLeeSet al. DSigDB: drug signatures database for gene set analysis. Bioinformatics. (2015) 31:3069–71. doi: 10.1093/bioinformatics/btv313
42
KuleshovMVJonesMRRouillardADFernandezNFDuanQWangZet al. Enrichr: a comprehensive gene set enrichment analysis web server. Nucleic Acids Res. (2016) 44:W90–7. doi: 10.1093/nar/gkw377
43
BustinSABenesVGarsonJAHellemansJHuggettJKubistaMet al. The MIQE guidelines. Clin Chem. (2009) 55:611–22. doi: 10.1373/clinchem.2008.112797
44
RådströmPKnutssonRWolffsPLövenklevMLöfströmC. Pre-PCR processing. Mol Biotechnol. (2004) 26:133–46. doi: 10.1385/MB:26:2:133
45
VandesompeleJDe PreterKPattynFPoppeBVan RoyNDe PaepeAet al. Accurate normalization of real-time quantitative RT-PCR data. Genome Biol. (2002) 3:research0034. doi: 10.1186/gb-2002-3-7-research0034
46
Student. The probable error of a mean. Biometrika. (1908) 6:1–25. doi: 10.1093/biomet/6.1.1
47
HuberWCareyVJGentlemanRAndersSCarlsonMCarvalhoBSet al. Orchestrating high-throughput genomic analysis with Bioconductor. Nat Methods. (2015) 12:115–21. doi: 10.1038/nmeth.3252
48
KuhnM. Building predictive models in R using the caret package. J Stat Software. (2008) 28:1–26. doi: 10.18637/jss.v028.i05
49
FriedmanJHastieTTibshiraniR. Regularization paths for generalized linear models via coordinate descent. J Stat Software. (2010) 33:1–22. doi: 10.18637/jss.v033.i01
50
LeekJTJohnsonWEParkerHSJaffeAEStoreyJD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. (2012) 28:882–3. doi: 10.1093/bioinformatics/bts034
51
LiberzonABirgerCThorvaldsdóttirHGhandiMMesirovJPTamayoP. The Molecular Signatures Database hallmark gene set collection. Cell Syst. (2015) 1:417–25. doi: 10.1016/j.cels.2015.12.004
52
BinnewiesMRobertsEWKerstenKChanVFearonDFMeradMet al. Understanding the tumor immune microenvironment (TIME) for effective therapy. Nat Med. (2018) 24:541–50. doi: 10.1038/s41591-018-0014-x
53
HotamisligilGS. Inflammation, metaflammation and immunometabolic disorders. Nature. (2017) 542:177–85. doi: 10.1038/nature21363
54
DivellaRDe LucaRAbbateINaglieriEDanieleASavinoEet al. Obesity and cancer: The role of adipose tissue and adipokines. Int J Mol Sci. (2016) 17:E870. doi: 10.3390/ijms17060870
55
LumengCNSaltielAR. Inflammatory links between obesity and metabolic disease. J Clin Invest. (2011) 121:2111–7. doi: 10.1172/JCI57132
56
HajishengallisGReisESMastellosDCRicklinDLambrisJD. Novel mechanisms of complement activation in inflammation. Nat Immunol. (2017) 18:1288–98. doi: 10.1038/ni.3858
57
BrownGDWillmentJAWhiteheadLMartinez-PomaresLGordonSTaylorPRet al. Dectin-1 mediates macrophage recognition of Candida albicans. Nat Immunol. (2018) 19:246–54. doi: 10.1038/s41590-017-0033-7
58
GrivennikovSIGretenFRKarinM. Immunity, inflammation, and cancer. Cell. (2010) 140:883–99. doi: 10.1016/j.cell.2010.01.025
59
RicklinDReisESMastellosDCGrosPLambrisJD. Complement component C5. Nat Rev Immunol. (2016) 16:5–18. doi: 10.1038/nri.2015.6
60
HunterCAJonesSA. IL-6 as a keystone cytokine in health and disease. Nat Immunol. (2015) 16:448–57. doi: 10.1038/ni.3153
61
JenneCNKubesP. Neutrophils in homeostasis and immunity. Immunity. (2015) 43:15–28. doi: 10.1016/j.immuni.2015.06.016
62
HothornTBühlmannPKneibTSchmidMHofnerB. Model-based boosting 2.0. J Mach Learn Res. (2010) 11:2109–13. doi: 10.5555/1756006.1859911
63
RobinXTurckNHainardATibertiNLisacekFSanchezJCet al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinf. (2011) 12:77. doi: 10.1186/1471-2105-12-77
64
VabalasAGowenEPoliakoffECassonAJ. Machine learning algorithm validation with a limited sample size. PloS One. (2019) 14:e0224365. doi: 10.1371/journal.pone.0224365
65
GetahunACambierJC. Of ITIMs, ITAMs, and ITAMis: revisiting immunoglobulin Fc receptor signaling. Immunol Rev. (2023) 314:181–96. doi: 10.1111/imr.13185
66
ZhengSZouYXieXLiangJChenHShiYet al. Development and validation of a stromal immune phenotype classifier for predicting immune activity and prognosis in triple-negative breast cancer. Int J Cancer. (2022) 151:56–67. doi: 10.1002/ijc.33932
67
NimmerjahnFRavetchJV. Fcgamma receptors as regulators of immune responses. Nat Rev Immunol. (2008) 8:34–47. doi: 10.1038/nri2206
68
XiaoYChenRLiuQZhangZZhangXWangJet al. Single-cell RNA sequencing reveals myeloid transcriptional heterogeneity in renal cell carcinoma. J Immunother Cancer. (2022) 10:e004231. doi: 10.1136/jitc-2021-004231
69
LundbergSMLeeSI. A unified approach to interpreting model predictions. Adv Neural Inf Process Syst. (2017) 30:4765–74. doi: 10.48550/arXiv.1705.07874
70
KaufmanLRousseeuwPJ. Finding groups in data: an introduction to cluster analysis. John Wiley & Sons (2009).
71
ChernickMR. Bootstrap methods: a guide for practitioners and researchers. 2nd ed. Wiley (2011).
72
ClarkDJZhangHClauserKRClaussTRVega-MontotoLTaniokaMet al. Integrated proteogenomic characterization of clear cell renal cell carcinoma. Cell. (2019) 179:964–83.e31. doi: 10.1016/j.cell.2019.10.007
73
KimDChenJSheppardNCPalendiraULangPASchartnerJMet al. Immunological hallmarks of stromal cells in the tumour microenvironment. J Clin Invest. (2021) 131:e146963. doi: 10.1172/JCI146963
74
SmithCCBeckermannKEBortoneDSDe CubasAABixbyLMLeeSJet al. Endogenous retroviral signatures in immunotherapy response and resistance. Nat Rev Urol. (2021) 18:529–38. doi: 10.1038/s41585-021-00483-z
75
BraunDAHouYBakounyZFicialMSant’ AngeloMFormanJet al. Interplay of somatic alterations and immune infiltration in renal cell carcinoma. Nat Med. (2020) 26:909–18. doi: 10.1038/s41591-020-0839-y
76
MenterTTzankovABergerCTNissenDJuskeviciusDSchärSet al. Multi-modal analysis of tumor microenvironment features in advanced renal cell carcinoma. Blood Adv. (2023) 7:2764–75. doi: 10.1182/bloodadvances.2022008741
77
HsiehJJChenZMarkerMRedzematovicAUntchBRReznikEet al. Clonal evolution and therapeutic vulnerabilities in renal cell carcinoma. Nat Rev Nephrol. (2023) 19:435–51. doi: 10.1038/s41581-023-00717-w
78
LiXWangZChenQLiuSZhangYZhouRet al. Structural basis of FCGR2A-TREM2 synergy in myeloid metabolic reprogramming. Nat Immunol. (2024) 25:312–25. doi: 10.1038/s41590-024-01776-0
79
KimDBachDHFanYHLeeSYKimDParkHJet al. SRC-mediated metabolic reprogramming dictates metastatic potential in obesity-associated renal cancer. J Clin Invest. (2021) 131:e146963. doi: 10.1172/JCI146963
80
CortelliniABersanelliMButiSPalaLMasiniCPerroneFet al. A multicenter study of body mass index and survival in metastatic renal cell carcinoma patients receiving nivolumab. J Immunother Cancer. (2021) 9:e003242. doi: 10.1136/jitc-2021-003242
81
MarxV. Method of the year: spatially resolved transcriptomics. Nat Methods. (2021) 18:9–14. doi: 10.1038/s41592-020-01033-y
82
GetahunACambierJC. Fc receptor signaling in cancer immunotherapy. Trends Immunol. (2023) 44:189–202. doi: 10.1016/j.it.2022.12.003
83
Sánchez-GastaldoAKempfEDuranI. Systemic treatment of metastatic clear cell renal cell carcinoma: moving towards a personalized approach. Crit Rev Oncol Hematol. (2017) 119:76–83. doi: 10.1016/j.critrevonc.2017.09.014
84
TurajlicSSwantonCBoshoffC. Kidney cancer: The next decade. J Exp Med. (2022) 219:e20211656. doi: 10.1084/jem.20211656
85
CoussensLMWerbZ. Inflammation and cancer. Nature. (2002) 420:860–7. doi: 10.1038/nature01322
86
WozniakMBLe Calvez-KelmFAbedi-ArdekaniBByrnesGDurandGCarreiraCet al. Integrative genome-wide gene expression profiling of clear cell renal cell carcinoma in Czech Republic and in the United States. Gene Expression Omnibus. (2013) 8(12):e57886. doi: 10.1371/journal.pone.0057886
87
ThibodeauBJFultonMFortierLEGeddesTJPruetzBLAhmedSet al. Characterization of clear cell renal cell carcinoma by gene expression profiling. Urol Oncol. (2016) 34(3):121.e9–15. doi: 10.1016/j.urolonc.2015.11.001
88
WotschofskyZGummlichLLiepJStephanCKilicEJungKet al. Integrated microRNA and mRNA signature associated with the transition from the locally confined to the metastasized clear cell renal cell carcinoma. PLoS One. (2016) 11(2):e0148746. doi: 10.1371/journal.pone.0148746
89
von RoemelingCARadiskyDCMarlowLACooperSJGrebeSKAnastasiadisPZet al. Neuronal pentraxin 2 supports clear cell renal cell carcinoma by activating the AMPA-selective glutamate receptor-4. Cancer Res. (2014) 74(11):2936–48. doi: 10.1158/0008-5472.CAN-14-0210
90
KulytéAEhrlundAArnerPAlbertonPRydenMGaoHet al. Global transcriptome profiling identifies KLF15 and SLC25A10 as modifiers of adipocytes insulin sensitivity in obese women. PLoS One. (2017) 12(5):e0178485. doi: 10.1371/journal.pone.0178485
91
National Center for Biotechnology Information. Gene Expression Omnibus: GSE94752. Bethesda, MD: National Center for Biotechnology Information (2021).
Summary
Keywords
clear cell renal cell carcinoma, obesity, FCGR2A, machine learning, immune microenvironment
Citation
He Z, Wang Z, Lai S, Yin X, Zheng D, Liu S, Liu W and Guo G (2025) Machine learning-driven dissection of the obesity-ccRCC interface: FCGR2A emerges as a central coordinator of tumor-immune crosstalk. Front. Oncol. 15:1598007. doi: 10.3389/fonc.2025.1598007
Received
22 March 2025
Accepted
22 September 2025
Published
22 October 2025
Volume
15 - 2025
Edited by
Abdallah Flaifel, New York University, United States
Reviewed by
Mohane S. Coumar, Pondicherry University, India
Abdullah Al Marzan, Dhaka Medical College and Hospital, Bangladesh
Updates
Copyright
© 2025 He, Wang, Lai, Yin, Zheng, Liu, Liu and Guo.
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: Guiying Guo, guoguiying790630@126.com
†These authors have contributed equally to this work
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.