PYCR in Kidney Renal Papillary Cell Carcinoma: Expression, Prognosis, Gene Regulation Network, and Regulation Targets

Background : Pyrroline-5-carboxylate reductase (PYCR) includes three human genes encoding three isozymes, PYCR1 , PYCR2 , and PYCR3 (or PYCRL ), which facilitate the final step in the conversion of glutamine to proline. These genes play important roles in regulating the cell cycle and redox homeostasis as well as promoting growth signaling pathways. Proline is abnormally upregulated in a variety of cancers, and as the last key enzyme in proline production, PYCR plays an integral role in promoting tumorigenesis and cancer progression. However, its role in patients with kidney renal papillary cell carcinoma (KIRP) has not been fully elucidated. In this study, we aimed to systematically analyze the expression, gene regulatory network, prognostic value, and target prediction of PYCR in patients with KIRP, elucidate the association between PYCR expression and KIRP, and identify potential new targets for the clinical treatment of KIRP. Methods : We systematically analyzed the expression, prognosis, gene regulatory network, and regulatory targets of PYCR1, PYCR2, and PYCRL in KIRP using multiple online databases including cBioPortal, STRING, MethSurv, GeneMANIA, Gene Expression Profiling Interactive Analysis (GEPIA), Metascape, UALCAN, LinkedOmics, and TIMER. Results : The expression levels of PYCR1 , PYCR2, and PYCRL were considerably upregulated in patients with KIRP based on sample type, sex, age, and individual cancer stage. PYCR1 and PYCR2 transcript levels were markedly upregulated in females than in males, and patients aged 21 – 40 years had higher PYCR1 and PYCR2 transcript levels than those in other age groups. Interestingly, PYCR2 transcript levels gradually decreased with age. In addition, the expressions of PYCR1 and PYCR2 were notably correlated with the pathological stage of KIRP. Patients with KIRP with low PYCR1 and PYCR2 expression had longer survival than those with high PYCR1 and PYCR2 expression. PYCR1 , PYCR2, and PYCRL were altered by 4%, 7%, and 6%, respectively, in 280 patients with KIRP. The methylation levels of cytosine-phosphate-guanine (CpG) sites in PYCR were markedly correlated with the prognosis of patients with KIRP. PYCR1 , PYCR2 , PYCRL, and their neighboring genes form a complex network of interactions. The molecular functions of the genes, as demonstrated by their corresponding Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, included calcium channel activity, phospholipid binding, RNA polymerase II-specificity, and kinase and GTPase-regulatory activities. PYCR1 , PYCR2 , and PYCRL targeted miR-21, miR-221, and miR-222, resulting in a better prognosis of KIRP. We analyzed mRNA sequencing data from 290 patients with KIRP and found that ADA , NPM3 , and TKT were positively associated with PYCR1 expression; PFDN2 , JTB , and HAX1 were positively correlated with PYCR2 expression; SHARPIN , YDJC , and NUBP2 were positively correlated with PYCRL expression; PYCR1 was positively correlated with B cell and CD8+ T-cell infiltration levels; macrophage infiltration was negatively correlated with PYCR2 expression; and PYCRL expression was negatively correlated with B-cell, CD8+ T cell, and dendritic cell infiltration levels. Conclusions : PYCR1 , PYCR2, and PYCRL may be potential therapeutic and prognostic biomarkers for patients with KIRP. The regulation of microRNAs (miRNAs), including miR-21, miR-221, and miR-222, may prove an important strategy for KIRP treatment.


Introduction
Cancer is the leading cause of human death and a major public health problem in most countries [1,2].Kidney renal papillary cell carcinoma (KIRP) accounts for 10-20% of renal cell carcinomas (RCC) and is the second most common type of RCC [3].Currently, the efficacy of therapeutic drugs remains unsatisfactory in clinical practice [4].The etiology of KIRP remains unclear.Some studies have reported that the loss or mutation of various oncogenes in patients with KIRP, which leads to the excitation or inhibition of their functions, is one of the inevitable initial steps in the occurrence and development of KIRP [5].Therefore, mining new biomarkers and potential regulatory targets is crucial to improve the morbidity and survival of patients with KIRP.
Pyrroline-5-carboxylate reductase (PYCR) includes three human genes encoding three isozymes, PYCR1, PYCR2, and PYCR3 (or PYCRL), which facilitate the final step in the conversion of glutamine to proline.PYCR is upregulated in various cancer tissues [6].High expression of PYCR is positively correlated with poor cancer prognosis.Increased expression of PYCR leads to elevated proline concentrations, which are thought to be metabolic addictions in cancer cells, making proline a trusted therapeutic target [7].Numerous studies have shown that PYCR plays a crucial role in tumor growth and progression [8][9][10].In addition, proline has been shown to act as a potent anti-oxidative stress agent capable of protecting tumor cells from various reactive oxygen species inducers.Taken together, downregulating the expression level of PYCR may be an effective cancer treatment strategy and is expected to be a target for blocking tumor progression and improving survival.
The role of PYCR in KIRP is not well understood.Therefore, this study systematically analyzed the expression, gene regulatory network, prognostic value, and target prediction of PYCR in patients with KIRP, elucidated the association between PYCR and KIRP, and identified new potential targets for KIRP therapy.

UALCAN
UALCAN (http://ualcan.path.uab.edu/analysis.html) is an online professional database for analyzing tumor gene expression and methylation levels [11][12][13].We used UAL-CAN to analyze the expression and methylation levels of PYCR1, PYCR2, and PYCRL in healthy subjects and patients with KIRP.We used the student's t-test for comparative analysis, and the difference was considered significant at a p-value less than 0.05.

GEPIA
Gene Expression Profiling (GEPIA) (http://gepia.cancer-pku.cn/index.html) is a free online platform for analyzing the correlation of gene expression levels with tumor pathological stage and prognostic value [13].We used GEPIA to analyze the pathological staging correlation and prognostic value of the expression level of PYCR1, PYCR2, and PYCRL in patients with KIRP.Student's t-test was used for comparative analysis, and differences were considered significant at a p-value less than 0.05.

cBioPortal
cBioPortal (http://cbioportal.org) is an online professional database used to analyze genetic alterations in tumors [11][12][13].We used the cBioPortal database to analyze genetic alterations in PYCR1, PYCR2, PYCRL, and their neighboring genes.A total of 280 KIRP samples were analyzed, and mRNA expression z-scores were obtained relative to all samples (log RNA Seq V2 RSEM) using a z-score threshold of ± 2.0.

STRING
STRING (https://string-db.org/cgi/input.pl) is an online professional database that analyzes protein-protein in-teractions (PPI) [13].We used STRING to build a lowconfidence level (0.150) PPI network and screening criteria for species defined as "humans".Finally, K-means cluster analysis was performed on PYCR1, PYCR2, PYCRL, and their neighboring genes.

Metascape
Metascape (https://metascape.org) is a professionalfree tool that analyzes Gene Ontology (GO) functions and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment [13].We used Metascape to analyze the GO function and KEGG pathway enrichment of PYCR1, PYCR2, PYCRL, and their altered neighboring genes in KIRP.

LinkedOmics
LinkedOmics (http://www.linkedomics.org/) is a free online platform used to analyze miRNA target enrichment and differentially expressed genes associated with tumor genes [13].The LinkedOmics database was used to analyze miRNA target enrichment and differentially expressed genes associated with PYCR1, PYCR2, and PYCRL.

MethSurv
MethSurv (https://biit.cs.ut.ee/methsurv/) is a free online platform that provides survival analyses based on DNA methylation.We used MethSurv to analyze DNA methylation and prognosis at PYCR sites in patients with KIRP.

PYCR Expression in Patients with KIRP
We first compared the expression levels of PYCR1, PYCR2, and PYCRL in patients with KIRP according to sample type, sex, age, and individual cancer stage.Our results showed that PYCR1, PYCR2, and PYCRL transcript levels were significantly upregulated in patients with KIRP (p < 0.05; Fig. 1).In patients with KIRP, PYCR1 and PYCR2 transcript levels were significantly upregulated in females compared with in males (p < 0.01; Fig. 1B,F), and patients aged 21-40 years had higher PYCR1 and PYCR2 transcript levels than patients in other age groups (p < 0.05; Fig. 1C,G).Interestingly, PYCR2 transcript levels gradually decreased with age (p < 0.05; Fig. 1G).The transcript levels of PYCR1 and PYCR2 in different cancer stages were significantly higher than those in normal individuals and gradually increased with an increase in cancer stage (p < 0.05; Fig. 1D,H).In addition, we evaluated the correlation between the differential expression of PYCR1 and PYCR2 and the pathological stage of KIRP.We found a significant correlation between the expression of PYCR1 (p = 6.02 × 10 −14 ; Fig. 2A) and PYCR2 (p = 0.00145; Fig. 2B) and the pathological stage of patients with KIRP.Finally, we used GEPIA to assess the prognostic value of PYCR1, PYCR2, and PYCRL expression in patients with KIRP.Our results showed that patients with KIRP with low PYCR1 expression had a longer overall survival than those with high PYCR1 expression (p = 0.00094; Fig. 3A).Similarly, patients with KIRP with low PYCR1 and PYCR2 expression had longer disease-free survival than those with high PYCR1 and PYCR2 expression (p = 0.03 and p = 0.017, re-spectively; Fig. 3B,D).

Genetic Alteration and Promoter Methylation of PYCR in Patients with KIRP
We further assessed the genetic alterations in PYCR1, PYCR2, and PYCRL in 280 patients with KIRP using TCGA.We found that PYCR1 was altered by 4% in patients with KIRP, with the type of genetic alteration mainly including amplification, high RNA levels, and low RNA levels (Fig. 4A).However, PYCR2 was altered by 7% in patients with KIRP, with the type of genetic alteration mainly including missense mutations, high RNA levels, and low RNA levels (Fig. 4B).PYCRL was altered by 6% in patients with KIRP, with the type of genetic alter- ation mainly including missense mutations, truncating mutations, deep deletions, high RNA levels, and low RNA levels (Fig. 4C).Next, we assessed the promoter methylation levels of PYCR1, PYCR2, and PYCRL in patients with KIRP using UALCAN.We found that PYCR1 and PY-CRL promoter methylation levels were significantly downregulated in patients with KIRP (p < 0.05; Fig. 4D,F).Among the 18 predicted CpG sites of PYCR1, cg25759517 and cg19202384 were significantly correlated with KIRP prognosis (Fig. 5A) (Table 1).Patients with high PYCR1 methylation at these CpG sites had better overall survival than those with low PYCR1 methylation (Fig. 5B,C).Furthermore, Among the nine predicted CpG sites of PYCR2, cg07049680, cg23091741, and cg06086141 were considerably correlated with KIRP prognosis (Fig. 5D) (Table 1).Patients with high PYCR2 methylation at CpG sites (cg07049680 and cg23091741) had worse overall survival than patients with low PYCR2 methylation (Fig. 5E,F).However, patients with high PYCR2 methylation at CpG sites (cg06086141) had better overall survival than patients with low PYCR2 methylation (Fig. 5G).Among the 12 predicted CpG sites of PYCRL, cg26507094 was markedly correlated with KIRP prognosis (Fig. 5H) (Table 1).Patients with high PYCRL methylation at these CpG sites had worse overall survival than those with low PYCRL methylation (Fig. 5I).

Neighboring Gene Alteration and PYCR Interaction Network in Patients with KIRP
We evaluated alterations in the neighboring genes of PYCR1, PYCR2, and PYCRL in patients with KIRP using cBioPortal.We found a gene alteration frequency of ≥25.00%, ≥10.00%, and ≥11.76% for the 50 most frequently altered neighboring genes of PYCR1, PYCR2, and PYCRL in patients with KIRP, respectively (Tables 2,3,4).The most frequently altered neighboring genes of PYCR1 in patients with KIRP were ALYREF (33.33%),ANAPC11 (33.33%), and ARHGDIA (33.33%) (Table 2).However, PBRM1 (20.00%),CDKN2A (20.00%), and ALK (15.00%) were the most frequently altered neighboring genes of PYCR2 in patients with KIRP (Table 3).The most frequently altered neighboring genes of PYCRL in patients with KIRP were MT-CO2 (23.53%),BLK (17.65%), and C2CD5 (17.65%) (Table 4).We further evaluated the potential interactions between PYCR1, PYCR2, PYCRL, and their neighboring genes.Our results showed that 43 nodes, 188 edges, and 3 clusters were obtained in the constructed PPI network of PYCR1 and its neighboring genes in patients with KIRP (Fig. 6A).Furthermore, PYCR1 and its neighboring genes were linked to a complex interaction network (70 genes and 226 edges) through co-expression, physical interactions, shared protein domains, and prediction (Fig. 6B).However, we found that 46 nodes, 158 edges, and 3 clusters were obtained in the constructed PPI network of PYCR2 and its neighboring genes in patients with KIRP (Fig. 6C).Furthermore, PYCR2 and its neighboring genes were linked to a complex interaction network (71 genes and 217 edges) through co-expression, physical interactions, genetic interactions, and prediction (Fig. 6D).Additionally, 33 nodes, 126 edges, and 3 clusters were obtained in the constructed PPI network of PYCRL and its neighboring genes in patients with KIRP (Fig. 6E).In addition, PYCRL and its neighboring genes were linked to a complex interaction network (60 genes and 188 edges) through co-expression, physical interactions, shared protein domains, and prediction (Fig. 6F).

GO and KEGG Pathway Enrichment Analyses
We further performed GO and KEGG pathway enrichment analyses of PYCR1, PYCR2, PYCRL, and their top 50 altered neighboring genes in patients with KIRP using Metascape.We found that the molecular functions related to PYCR1 and its neighboring genes mainly included voltage-gated calcium channel activity, phospholipid binding, and ubiquitin-like protein ligase binding (Fig. 7A).Their cellular components include nuclear ubiquitin ligase complexes, tertiary granules, and cation channel complexes (Fig. 7B).Their biological processes included protein localization to lysosomes, regulation of wound healing, and negative regulation of cell activation (Fig. 7C).Arrhythmogenic right ventricular cardiomyopathy and cardiac muscle contraction were enriched according to KEGG pathway analysis (Fig. 7D).However, our results showed that nuclear receptor binding, DNA-binding transcription activator activity, RNA polymerase II-specific, phosphotransferase activity, alcohol group as acceptor, and oxidoreductase activity were molecular functions related to PYCR2 and its neighboring genes (Fig. 7E).The RNA polymerase II transcription regulator complex and the cytoplasmic region were their cellular components (Fig. 7F).Cellular response to oxidative stress, regulation of the G1/S transition of the mitotic cell cycle, and the enzyme-linked receptor protein signaling pathway were the main biological processes (Fig. 7G).KEGG pathway analysis revealed that PYCR2 and its neighboring genes were associated with pathways in cancer (Fig. 7H).Furthermore, the molecular functions related to PYCRL and its neighboring genes mainly include oxidoreductase activity, kinase activity, and GTPase regulator activity (Fig. 7I).The cellular components were spindles (Fig. 7J).The mitotic cell cycle process, positive regulation of protein localization, small molecule biosynthetic process, and chemotaxis were the main biological processes (Fig. 7K).Arginine and proline metabolism and amino acid biosynthesis were found to be enriched according to KEGG pathway analysis (Fig. 7L).

Discussion
PYCR has been reported to be overexpressed in various types of tumors, such as bladder, prostate, and gastric cancers [14].However, PYCR expression has not yet been reported in patients with KIRP.Initially, we compared the expression levels of PYCR in patients with KIRP according to sample types, sex, age, and individual cancer stage and found that PYCR transcript levels were consider-ably upregulated in patients with KIRP.Interestingly, our results showed that PYCR1 and PYCR2 transcript levels were markedly upregulated in females than in males, and patients aged 21-40 years had higher PYCR1 and PYCR2 transcript levels than patients in other age groups.However, studies have confirmed that the incidence of RCC in males is higher than that in females [15].Whether PYCR1 and PYCR2 differences in expression by sex and age are associated with the onset and prognosis of KIRP has not yet been reported.In addition, the number of cases in some individual groups was small.Hence, this requires further investigation.Furthermore, we found an extensive positive correlation between PYCR1, PYCR2, and the pathological stage of patients with KIRP.Our results also showed that PYCR1, PYCR2, and PYCRL were altered by 4%, 7%, and 6%, respectively, in 280 patients with KIRP.We also found that PYCR1 and PYCRL promoter methylation levels were considerably downregulated in patients with KIRP.Genetic changes and promoter methylation often result in abnormal gene expression and function.Increased expression of PYCR1 and PY-CRL caused by genetic changes and promoter methylation may also serve as an important factor in KIRP progresion.In addition, we found that the methylation levels of CpG sites in PYCR (PYCR1: cg25759517 and cg19202384; PYCR2: cg07049680, cg23091741, and cg06086141; PY-CRL: cg26507094) were notably correlated with the prognosis of patients with KIRP.This revealed that the methyla- tion levels of PYCR act as effective prognostic biomarkers for KIRP, indicating that PYCR may play a critical role in KIRP progression.The most frequently altered neighboring genes of PYCR for PYCR1 were (ALYREF, ANAPC11, and ARHGDIA), PYCR2 were (PBRM1, CDKN2A, and ALK), and PYCRL were (MT-CO2, BLK, and C2CD5).In cancer patients with high PYCR1 expression, increased PYCR1 activity is associated with higher expression of oncogenes, such as ALYREF, ANAPC11, and ARHGDIA.PYCR1-driven oncogenes promote tumor cell proliferation and metastasis [16][17][18].PBRM1, CDKN2A, and ALK are frequently mutated in lung cancer and are associated with tumor resistance [19].MT-CO2 is a mitochondrial gene whose mutation has been widely reported in various human tumors.MT-CO2 variation may be a potential prognostic biomarker in MUTYH-associated polyposis patients [20].CDKN2A is the second most common tumor suppressor gene in cancers.CDKN2A is frequently altered in chordoma [21].In conclusion, these gene alterations may be involved in the occurrence and development of KIRP to different degrees and in different pathways.Finally, we assessed the prognostic value of PYCR expression in patients with KIRP.Our results showed that patients with KIRP with low PYCR1 and PYCR2 expression had longer survival rates than those with high PYCR1 and PYCR2 expression.Thus, PYCR1 and PYCR2 may serve as potential prognostic markers in patients with KIRP.
We further evaluated the potential interactions between PYCR and its neighboring genes.We found that PYCR and its neighboring genes are linked to a complex interaction network through co-expression, physical interactions, and prediction.Next, we evaluated the functions of PYCR and its neighboring genes.We found that the molecular functions of PYCR1 and its neighboring genes were tumor-associated voltage-gated calcium channel activity, phospholipid binding, and ubiquitin-like protein ligase binding.For instance, the voltage-gated calcium channel is closely related to cancer pain and inhibiting its activity can reduce cancer pain [22].In addition, inhibition of T-type voltage-gated calcium channels showed antiproliferative and cytotoxic effects and is an important anticancer target [23].Neuronally expressed developmentally downregulated 4 (NEDD4) is a ubiquitin-like protein ligase that accelerates tumor growth and metastasis by regulating the microRNA-132/NRF2 axis [24].NEDD4 is a potential therapeutic target for the prevention of tumor recurrence and metastasis.In addition, we found that nuclear receptor binding, DNA-binding transcription activator activity, and RNA polymerase II-specificity, which are related to tumorigenesis and development, were the molecular functions of PYCR2 and its neighboring genes.For example, overexpression of nuclear receptor-binding protein 1 inhibits colorectal cancer cell proliferation and promotes apoptosis [25].The KEGG pathway of PYCR2 and its neighboring genes is involved in cancer pathways.In addition, the molecular functions associated with PYCRL and its neighboring genes mainly include oxidoreductase, kinase, and GTPase-regulatory activities related to tumor metabolism.Their KEGG pathway analysis revealed their involvment in proline metabolism.Proline has been shown to act as a potent anti-oxidative stress agent capable of protecting cells from various reactive oxygen species (ROS) inducers.Excessive ROS production is detrimental, through oxidative damage to DNA, proteins, and lipids, to the survival of cancer cells.In summary, regulating the functions of PYCR and its neighboring genes may be an important strategy for the treatment of KIRP.
We explored the correlation between differentially expressed genes and PYCR expression in patients with KIRP.We found that the expression of 20,023 genes was correlated with PYCR1 expression.ADA, NPM3, and TKT were the top three genes positively correlated with the expression of PYCR1.Adenosine deaminase (ADA) regulates intracellular and extracellular adenosine concentrations by converting adenosine to inosine.Adenosine accumulates in high concentrations in tumors and exerts a series of malignancy-promoting effects [40].Transketolase (TKT) is a key enzyme in the non-oxidative pentose phosphate pathway.TKT can promote the development of liver cancer in a non-metabolic manner through its nuclear localization and EGFR pathway [41].In addition, we found that PFDN2, JTB, and HAX1 were the top three genes positively correlated with the expression of PYCR2.Prefoldin subunits (PFDN2) are involved in cancer progression.PFDN2 may serve as poor prognostic markers for gastric cancer [42].JTB is a tumor suppressor gene in various malignancies and can be used as a biomarker for breast cancer [43].Studies have confirmed that HS-1-associated protein-1 (HAX1) enhances non-small cell lung cancer survival and metastasis through the AKT/mammalian target of the rapamycin (MTOR) and double minute 2 protein (MDM2)/P53 signaling pathways [44].Furthermore, SHARPIN, YDJC, and NUBP2 were the top three genes positively correlated with the expression of PYCRL.SHARPIN plays an important role in tumorigenesis.Overexpression of SHARPIN promotes tumor progression in ovarian cancer [45].Chitooligosaccharide deacetylase homolog (YDJC) contributes to induction of protein phosphatase 2A (PP2A) ubiquitination by increasing epithelial-mesenchymal transition (EMT).Hence, YDJC may become a new target for antitumor therapy [46].Therefore, targeting these genes may provide adjuvant therapy for KIRP.
Immunotherapy is a novel clinical cancer treatment.Tumor immune infiltration is closely related to clinical prognosis.As expected, the expression level of PYCR in patients with KIRP was positively or negatively correlated with the level of immune cell infiltration.B-cell-mediated inflammation may play an important role in promoting the progression of aggressive malignancies.However, data from preclinical studies have shown significant differences in the role of B cells at different stages of cancer development.In mouse models of precancerous lesions, B cells appear to cause inflammation, which in turn induces precancerous lesions [47].CD8+ T cells play a central role in antitumor immunity.CD8+ T cells are activated upon stimulation by tumor antigens and then differentiate into functional cells and migrate to the tumor site.Activated CD8+ T cells can directly destroy tumor cells by releasing perforin and granzymes and by inducing death ligand/death receptormediated apoptosis.They also secrete cytokines to regulate the immune system in tumor cells.However, disturbances in mitochondrial dynamics can affect CD8+ T cell function.Overall, mitochondrial abnormalities in CD8+ T cells contribute to cancer development [48].Therefore, targeting PYCR or PYCR-related regulatory targets may be a feasible strategy for regulating immune cell infiltration in patients with KIRP.

Conclusions
This study systematically analyzed the expression, gene regulatory network, prognostic value, and target prediction of PYCR in patients with KIRP, elucidated the relationship between PYCR and KIRP, and provided new insights into KIRP treatment.

Fig. 1 .
Fig. 1.The expression of PYCR in patients with KIRP (UAL-CAN).(A) The transcription expression of PYCR1 based on sample types.(B) The transcription expression of PYCR1 based on the sex of the patient.(C) The transcription expression of PYCR1 based on the age of the patient.(D) The transcription expression of PYCR1 based on individual cancer stage.(E) The transcription expression of PYCR2 based on sample types.(F) The transcription expression of PYCR2 based on the sex of the patient.(G) The transcription expression of PYCR2 based on the age of the patient.(H) The transcription expression of PYCR2 based on individual cancer stage.(I) The transcription expression of PYCRL based on sample types.(J) The transcription expression of PYCRL based on the sex of the patient.(K) The transcription expression of PYCRL based on the age of the patient.(L) The transcription expression of PYCRL based on individual cancer stage.*p < 0.05; **p < 0.01; ***p < 0.001.

Fig. 3 .
Fig. 3.The prognostic value of PYCR in KIRP (GEPIA).(A) The overall survival curve of PYCR1 in patients with KIRP.(B) The disease-free survival cure of PYCR1 in patients with KIRP.(C) The overall survival curve of PYCR2 in patients with KIRP.(D) The disease-free survival cure of PYCR2 in patients with KIRP.(E) The overall survival curve of PYCRL.(F) The diseasefree survival cure of PYCRL in patients with KIRP.

Fig. 4 .
Fig. 4. Genetic alteration and promoter methylation of PYCR in KIRP (UALCAN).(A) Genetic alteration of PYCR1in patients with KIRP.(B) Genetic alteration of PYCR2 in patients with KIRP.(C) Genetic alteration of PYCRL in patients with KIRP.(D) Promoter methylation of PYCR1in healthy individuals and patients with KIRP.(E) Promoter methylation of PYCR2 in healthy individuals and patients with KIRP.(F) Promoter methylation of PYCRL in healthy individuals and patients with KIRP.

Fig. 5 .
Fig. 5. DNA methylation at CpG sites and prognostic value of PYCR in KIRP (MethSurv).(A) Heatmap showing the PYCR1 DNA methylation at CpG sites.(B) Patients with KIRP with higher PYCR1 methylation of cg25759517 CpG sites had a better overall survival than those with lower PYCR1 methylation (HR = 0.512, p = 0.038).(C) Patients with KIRP with higher PYCR1 methylation of cg19202384 CpG sites had a better overall survival than those with lower PYCR1 methylation (HR = 0.452, p = 0.015).(D) Heatmap showing the PYCR2 DNA methylation at CpG sites.(E) Patients with KIRP with higher PYCR2 methylation of cg07049680 CpG sites had a worse overall survival than those with lower PYCR2 methylation (HR = 1.907, p = 0.045).(F) Patients with KIRP with higher PYCR2 methylation of cg23091741 CpG sites had a worse overall survival than those with lower PYCR2 methylation (HR = 2.265, p = 0.014).(G) Patients with KIRP with higher PYCR2 methylation of cg23091741 CpG sites had a better overall survival than those with lower PYCR1 methylation (HR = 0.422, p = 0.011).(H) Heatmap showing the PYCRL DNA methylation at CpG sites.(I) Patients with KIRP with higher PYCRL methylation of cg26507094 CpG sites had a worse overall survival than those with lower PYCRL methylation (HR = 2.611, p = 0.0033).

Fig. 6 .
Fig. 6.Interaction analyses of PYCR and their neighboring genes in KIRP (STRING and GeneMANIA).(A) PPI network of PYCR1 and its neighboring genes in patients with KIRP (STRING).(B) Network analyses of PYCR1 and its neighboring genes in patients with KIRP (GeneMANIA).(C) PPI network of PYCR2 and its neighboring genes in patients with KIRP (STRING).(D) Network analyses of PYCR2 and its neighboring genes in patients with KIRP (GeneMANIA).(E) PPI network of PYCRL and its neighboring genes in patients with KIRP (STRING).(F) Network analyses of PYCRL and its neighboring genes in patients with KIRP (GeneMANIA).

Fig. 7 .
Fig. 7. GO function and KEGG pathways enrichment analyses of PYCR1, PYCR2, PYCRL, and their neighboring genes in patients with KIRP (Metascape).(A) Molecular functions of PYCR1 and its neighboring genes.(B) Cellular components of PYCR1 and its neighboring genes.(C) Biological processes of PYCR1 and its neighboring genes.(D) KEGG pathway analysis of PYCR1 and its neighboring genes.(E) Molecular functions of PYCR2 and its neighboring genes.(F) Cellular components of PYCR2 and its neighboring genes.(G) Biological processes of PYCR2 and its neighboring genes.(H) KEGG pathway analysis of PYCR2 and its neighboring genes.(I) Molecular functions of PYCRL and its neighboring genes.(J) Cellular components of PYCRL and its neighboring genes.(K) Biological processes of PYCRL and its neighboring genes.(L) KEGG pathway analysis of PYCRL and its neighboring genes.

Fig. 8 .
Fig. 8. Genes differentially expressed in correlation with PYCR expression in KIRP (LinkedOmics).(A) Pearson test was used to analyze correlations between PYCR1 expression and genes differentially expressed in patients with KIRP.(B,C) Heatmaps showing genes positively and negatively correlated, respectively, with PYCR1 in patients with KIRP (top 50 genes).(D,E,F) The scatter plots show Pearson's correlation of PYCR1 expression with expression of ADA, NPM, and TKT, respectively, in patients with KIRP.(G) Pearson test was used to analyze correlations between PYCR2 expression and genes differentially expressed in patients with KIRP.(H,I) Heatmaps showing genes positively and negatively correlated, respectively, with PYCR2 in patients with KIRP (top 50 genes).(J,K,L) The scatter plots show Pearson's correlation of PYCR2 expression with expression of PFDN2, JTB, and HAX1, respectively, in patients with KIRP.(M) Pearson test was used to analyze correlations between PY-CRL expression and genes differentially expressed in patients with KIRP.(N,O) Heatmaps showing genes positively and negatively correlated, respectively, with PYCRL in patients with KIRP (top 50 genes).(P,Q,R) The scatter plots show Pearson's correlation of PYCRL expression with expression of SHARPIN, YDJC, and NUBP2, respectively, in patients with KIRP.