System analysis of VEGFA in renal cell carcinoma: The expression, prognosis, gene regulation network and regulation targets
Authors: Situ Yongli , Xu Qinying , Deng Li , Zhu Yan , Gao Ruxiu , Lei Lei , Shao Zheng 
Published: December 2021

Background

VEGFA is one of the most important regulators of angiogenesis and plays a crucial role in cancer angiogenesis and progression. Recent studies have highlighted a relationship between VEGFA expression and renal cell carcinoma occurrence. However, the expression level, gene regulation network, prognostic value, and target prediction of VEGFA in renal cell carcinoma remain unclear. Therefore, system analysis of the expression, gene regulation network, prognostic value, and target prediction of VEGFA in patients with renal cell carcinoma is of great theoretical significance as there is a clinical demand for the discovery of new renal cell carcinoma treatment targets and strategies to further improve renal cell carcinoma treatment efficacy.

Methods

This study used multiple free online databases, including cBioPortal, TRRUST, GeneMANIA, GEPIA, Metascape, UALCAN, LinkedOmics, Metascape, and TIMER for the abovementioned analysis.

Results

VEGFA was upregulated in patients with kidney renal clear cell carcinoma (KIRC) and kidney chromophobe (KICH), and downregulated in patients with kidney renal papillary cell carcinoma (KIRP). Moreover, genetic alterations of VEGFA were found in patients with renal cell carcinoma as follows: 4% (KIRC), 8% (KICH), and 4% (KIRP). The promoter methylation of VEGFA was lower and higher in patients with clinical stages of KIRC and stage 1 KIRP, respectively. VEGFA expression significantly correlated with KIRC and KIRP pathological stages. Furthermore, patients with KICH and KIRP having low VEGFA expression levels had a longer survival than those having high VEGFA expression levels. VEGFA and its neighboring genes functioned in the regulation of protein methylation and glycosylation, as well as muscle fiber growth and differentiation in patients with renal cell carcinoma. Gene Ontology enrichment analysis revealed that the functions of VEGFA and its neighboring genes in patients with renal cell carcinoma are mainly related to cell adhesion molecule binding, catalytic activity, acting on RNA, ATPase activity, actin filament binding, protease binding, transcription coactivator activity, cysteine-type peptidase activity, and calmodulin binding. Transcription factor targets of VEGFA and its neighboring genes in patients with renal cell carcinoma were found: HIF1A, TFAP2A, and ESR1 in KIRC; STAT3, NFKB1, and HIPK2 in KICH; and FOXO3, TFAP2A, and ETS1 in KIRP. We further explored the VEGFA-associated kinase (ATM in KICH as well as CDK1 and AURKB in KIRP) and VEGFA-associated microRNA (miRNA) targets (MIR-21 in KICH as well as MIR-213, MIR-383, and MIR-492 in KIRP). Furthermore, the following genes had the strongest correlation with VEGFA expression in patients with renal cell carcinoma: NOTCH4, GPR4, and TRIB2 in KIRC; CKMT2, RRAGD, and PPARGC1A in KICH; and FLT1, C6orf223, and ESM1 in KIRP. VEGFA expression in patients with renal cell carcinoma was positively associated with immune cell infiltration, including CD8+T cells, CD4+T cells, macrophages, neutrophils, and dendritic cells.

Conclusions

This study revealed VEGFA expression and potential gene regulatory network in patients with renal cell carcinoma, thereby laying a foundation for further research on the role of VEGFA in renal cell carcinoma occurrence. Moreover, the study provides new renal cell carcinoma therapeutic targets and prognostic biomarkers as a reference for fundamental and clinical research.