key: cord-0921542-qvxxrw1a authors: Niu, Xinhao; Zhu, Zhe; Shao, Enming; Bao, Juan title: ACE2 Is a Prognostic Biomarker and Associated with Immune Infiltration in Kidney Renal Clear Cell Carcinoma: Implication for COVID-19 date: 2021-01-28 journal: J Oncol DOI: 10.1155/2021/8847307 sha: 96e970c7c51f092f42b6045387c158e30baca31a doc_id: 921542 cord_uid: qvxxrw1a BACKGROUND: KIRC is one of the most common cancers with a poor prognosis. ACE2 was involved in tumor angiogenesis and progression in many malignancies. The role of ACE2 in KIRC is still ambiguous. METHODS: Various bioinformatics analysis tools were investigated to evaluate the prognostic value of ACE2 and its association with immune infiltration in KIRC. RESULTS: ACE2 was shown to be downregulated in KIRC at the mRNA and protein level. Low expression of ACE2 protein in KIRC patients was observed in subgroup analyses based on gender, age, weight, tumor grade, and cancer stage. Upregulation of ACE2 in KIRC was associated with a favorable prognosis. ACE2 mRNA expression showed a positive correlation with the abundance of immune cells (B cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells) and the level of immune markers of different immune cells in KIRC. ACE2 expression could affect, in part, the immune infiltration and the advanced cancer stage. Moreover, enrichment analysis revealed that ACE2 in KIRC were mainly involved in translation factor activity, immunoglobulin binding, metabolic pathways, transcriptional misregulation in cancerous cells, cell cycle, and ribosomal activity. Several ACE2-associated kinases, miRNA, and transcription factor targets in KIRC were also identified. CONCLUSION: ACE2 was downregulated in KIRC and served as a prognostic biomarker. It was also shown to be associated with immune infiltration. Kidney cancer is one of the most common malignances globally, accounting for about 4.5% of all newly diagnosed malignances [1] . It is anticipated that 73,750 people would be newly diagnosed with kidney cancer and 14,830 patients are likely to die because of the disease in the USA in 2020 [2] . e most common subtype of renal cancer is kidney renal clear cell carcinoma (KIRC), which makes up over 70% of kidney cancers [3] . Surgery excision remains the primary therapy for KIRC due to the growing resistance to radiotherapy and chemotherapy [4] . Much worse, the prognosis of KIRC patients tends to be poor, especially for patients in an advanced stage. e five-year overall survival rate of stage IV patients is less than 10% [5] . Previous studies have revealed that immune infiltration is significantly linked to the survival of KIRC patients. [6, 7] . Immunotherapy has been suggested as the treatment for metastatic KIRC [8, 9] . erefore, clarifying the role of immune infiltration in KIRC and identifying immuneassociated markers for the prognosis for KIRC are particularly necessary. Angiotensin converting enzyme 2 (ACE2) is a member of the renin angiotensin system (RAS) whose open reading framework encodes an 805-amino-acid polypeptide [10] . Increasing evidence indicates a significant function of ACE2 in the tumor angiogenesis and its progression in many cancers, such as thyroid carcinoma, breast carcinoma, and lung adenocarcinoma [11] [12] [13] . ACE2 has also been suggested as a biomarker for many diseases, including squamous cell/adenosquamous carcinoma, endometrial carcinoma, and hypertension [10, 14, 15] . However, limited studies have clarified the function of ACE2 in immune infiltration and its role in the prognosis in KIRC. Coronavirus disease 2019 , caused by the novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was initially found in Wuhan of China since December 2019 [16, 17] . It is well known that the functional host receptor of SARS-CoV-2 is ACE2 [18, 19] . Over 10 million peoples were diagnosed with COVID-19 and over 520000 peoples died of this disease globally until July 1, 2020. As we have seen, the prognosis of COVID-19 patients with KIRC remains ambiguous. erefore, our study was performed to detect ACE2 levels and the prognostic value in patients with KIRC. e function of ACE2 in immune infiltration in KIRC was also clarified. Our results may provide additional evidence regarding the role of ACE2 and immune infiltration in patients with KIRC. and Human Protein Atlas. ACE2 expression in KIRC was identified in the Oncomine (https://www.oncomine.org/), UALCAN (http://ualcan.path.uab.edu/cgi-bin/ualcan-res. pl), and Human Protein Atlas (https://www.proteinatlas. org/). ACE2 mRNA levels in various malignances, including KIRC, were determined with the Oncomine database and the threshold was set to the P value � 0.05 and fold change (FC) � 2, as well as gene ranking � top 10% [20] . In order to further detect the ACE2 protein expression in KIRC, we then used UALCAN and Human Protein Atlas. Based on data from Clinical Proteomic Tumor Analysis Consortium (CPTAC), UALCAN could be also used to detect the ACE2 protein expression in various subtribes of patients with KIRC [21] . e Human Protein Atlas is a program designed to map all of the human proteins in the cells, tissues, and organs [22] . Immunohistochemical staining of ACE2 in KIRC was obtained from the Human Protein Atlas. Plotter. In order to evaluate the significance of ACE2 level in the prognosis of KIRC, GEPIA (http://gepia.cancer-pku.cn/) [23] , OSkirc (http://bioinfo.henu.edu.cn/KIRC/KIRCList. jsp) [24] , and KM plotter databases (https://kmplot.com/) were applied separately. e median value of ACE2 expression was utilized to identify high/low ACE2 expression patients and the P value was set as 0.05. In Meier plotter, subgroup prognosis analysis based on different clinicopathologic features and immune cells in KIRC was performed using TCGA KIRC dataset. Analysis. TIMER (https:// cistrome.shinyapps.io/timer/) is a comprehensive tool providing immune infiltrates analysis across TCGA tumors [25] . Immune cell infiltration and immune biomarker expression were correlated with ACE2 and were evaluated with Spearman's correlation analysis using the TCGA KIRC dataset. e immune cells included were B cells, CD4+ Tcells, CD8+ Tcells, neutrophils, macrophages, and dendritic cells. Immune biomarkers were excluded because they have already been described in previous studies [26] [27] [28] . Analysis. cBioPortal (http://www.cbioportal.org) is a TCGA visual tool used to perform genome analysis [29] . We analyzed ACE2 genetic alteration in KIRC with the threshold as ±2.0 in mRNA expression z-scores (RNASeq V2 RSEM) and protein expression z-scores (RPPA). In order to verify the ACE2-associated functions in KIRC, LinkedOmics (http://www.linkedomics.org/), a comprehensive tool for multiomics analysis, was used [30] . A Pearson correlation test was used to explore genes that are linked to ACE2 in KIRC, while gene set enrichment analysis (GSEA) was performed for the enrichment analyses (GO and KEGG pathways), and ACE2-associated targets (kinase, miRNA, and transcription factor) were obtained with GSEA. ese analyses were carried out using the TCGA KIRC dataset, with a P value < 0.05. We initially detected the mRNA and protein expression of ACE2 in KIRC in Oncomine, UALCAN, and Human Protein Atlas. According to the data from Oncomine, ACE2 mRNA expression was dramatically reduced in KIRC when compared with normal kidney tissues (Figures 1(a) -1(c)). A gene expression profile also revealed that ACE2 mRNA expression was reduced in KIRC when compared with normal kidney tissues, with an FC of −2.843 as well as a P value of 0.01 (Figure 1 (b)) [31] . Another study indicated that ACE2 mRNA is expressed 5.131 times more in renal tissues than in KIRC tissues ( Figure 1 (c), P � 1.50E − 10) [32] . In order to further verify these results, we decided to use the CPTAC dataset to observe ACE2 protein expression. As expected, the results demonstrated a downwards regulation of ACE2 protein expression in KIRC when compared with normal kidney tissues (Figure 1(d) ). ACE2 protein expression was detected with staining and the expression data from the Human Protein Atlas. Interestingly, the immunohistochemical staining map suggested a low protein expression of ACE2 in KIRC tissues with a high protein expression of ACE2 in normal kidney tissues (Figure 1(e) ). However, we evaluated ACE2 protein expression in various subtribes of patients with KIRC. e results are shown in Figure 2 . is indicates a low expression of the ACE2 protein in KIRC patients in the subtribes analyses based on gender, age, weight, tumor grade, and cancer stage. erefore, ACE2 was downregulated in KIRC and may be involved in tumor progression. A Kaplan-Meier curve was applied using TCGA KIRC and GSE29609 datasets for prognosis analysis. KIRC patients 2 Journal of Oncology Data are mean ± SE. * P < 0.05; * * P < 0.01; * * * P < 0.001. with a high level of ACE2 expression were strongly correlated with better overall survival (OS) (Figure 3 (a), logrank P � 1.1e − 05) and disease-free survival (DFS) rates ( Figure 3 (b), logrank P � 0.000034). us, ACE2 could potentially serve as a prognostic biomarker in KIRC patients. e correlation between ACE2 expression and clinical characteristics of KIRC patients in the Kaplan-Meier plot was also explored to see how ACE2 expression affects the prognosis of patients with KIRC. As shown in Table 1 , increasing levels of ACE2 were linked to better prognosis in male and female patients and high/low mutation burden patients (all P < 0.05). Moreover, an increased expression level of ACE2 was linked to better prognosis in tumor grades 2 to 4 of KIRC patients. However, there is not enough data about KIRC patients in tumor grade 1 to perform the same analysis. Specifically, the increasing level of ACE2 was linked to better prognosis in cancers in stages 2 to 4 of KIRC patients (All P < 0.05) but was not linked to better prognosis in cancer stage 1 patients (HR � 0.57, P � 0.069, Table 1 ). ese data demonstrate that ACE2 expression could affect the prognosis of KIRC patients with advanced cancer stage. Previous studies have highlighted the significance of the tumor immune infiltration in the prognosis of renal cancer [6, 33] . erefore, we evaluated the correlation between ACE2 mRNA expression and immune infiltration in KIRC using the TIMER database. Interestingly, ACE2 mRNA expression showed a positive link to the abundance of B cells (P � 9.78e − 07), CD8+ T cells (P � 0.00395), macrophages (P � 0.0275), neutrophils (P � 0.00742), and dendritic cells (P � 0.0116) (Figure 4 (a)). Conversely, the copy number alteration of ACE2 could inhibit immune infiltration ( Figure 4 (b)). We further investigated if the expression of ACE2 was associated with immune markers of different immune cells in KIRC. As expected, a significant correlation was obtained between the expression of ACE2, and most of the immune markers in KIRC after tumor purity modulation were performed (Table 2) . Specifically, ACE2 was strongly linked to CD8A and CD8B (CD8+ T cell), CD19 and CD79A (B cell), CD86, and CD115 (monocyte), as well as CCL2 and CD68 (TAM). ACE2 was also positively linked to all markers of M1 macrophage (INOS, IRF5, and PTGS2). Moreover, ACE2 levels showed a positive association with most markers of natural killer cell (KIR2DL1, KIR2DL3, KIR2DL4, KIR3DL1, and KIR3DL2), Dendritic cell (HLA-DPB1, HLA-DQB1, HLA-DRA, HLA-DPA1, CD1 C, and NRP1), and 2 (GATA3, STAT6, and STAT5A). Similarly, ACE2 in KIRC showed a positive correlation with STAT3 in 17, FOXP3, STAT5B, and TGFB1 in Treg, as well as TIM-3 in T cell exhaustion ( Table 2) . Taken together, ACE2 was associated with tumor immune infiltration in KIRC, and ACE2 may play a vital role in immune escape in the KIRC microenvironment. e abovementioned results found that the expression levels of ACE2 were associated with favorable prognoses and immune infiltration in patients with KIRC. A prognostic analysis was performed to verify if the expression of ACE2 affects prognosis and immune infiltration in KIRC. is was based on immune cells using the Kaplan-Meier plotter. As we could see in Figure 5 , high expression of ACE2 in KIRC from the cohorts of enriched/ decreased basophils ( Figure 5 . erefore, ACE2 may affect the prognosis of patients with KIRC, in part, due to immune infiltration. Genomic mutations are known to be significantly linked to tumorigenesis. In our study, genetic alteration analysis of ACE2 in KIRC patient datasets revealed that a total of 9% of genetic alterations in ACE2 in KIRC and the genetic alteration form contained missense mutations, truncating mutations, deep deletions, and low mRNA (Figure 7 (a)). Moreover, ACE2 mutations could lead to protein change, including E489 * and I21 V (Figure 7(b) ). Interestingly, we found that ACE2 alterations in KIRC predicted a worse overall survival rate (P � 0.00121, Figure 7 (c)). ese findings suggest that an ACE2 genetic alteration may regulate tumorigenesis and its progression to KIRC, thus impacting the patients' prognosis. e TCGA KIRC dataset was analyzed with LinkedOmics. Figure 8 (a) shows that 3792 genes were positively linked to ACE2, and 6892 genes were negatively linked to ACE2 (false discovery rate <0.01). e top 50 significant genes that showed a positive and negative correlation with ACE2 were also obtained (Figures 8(b) and 8(c) ). GSEA was performed to analyze GO enrichment analysis, which revealed that ACE2 in KIRC were mainly involved in extracellular structure organization, small molecule catabolic processes, cellular amino acid metabolic processes, translation factor activity, structural constituent of ribosomes, immunoglobulin binding, cytokine receptor binding, and RNA binding (Figure 8 To further clarify the underlining mechanisms of how ACE2 affects tumorigenesis and the progression of KIRC, we explored ACE2-associated kinase, miRNA, and transcription factor targets in KIRC using GSEA in LinkedOmics. As a result, the top five most significant ACE2-associated kinase targets in KIRC were Kinase_LCK, Kinase_LYN, Kinase_SYK, Kinase_JAK3, and Kinase_HCK (Table 3, Table 3 , all P < 0.05). In the transcription factor target analysis, the results demonstrated V$IRF_Q6, V$NFKB_Q6_01, V$ELF1_Q6, V$PEA3_Q6, and V$PU1_Q6 as the ACE2-associated targets in KIRC (Table 3 , all P < 0.05). ACE2, a novel identified component of RAS, could regulate the tumorigenesis and progression in cancers and serve as a biomarker for many diseases [34] [35] [36] . Moreover, increasing evidence highlights the association between ACE2, tumor microenvironment, and immune infiltration [10, 37] . However, there were limited studies that clarified the function of ACE2 in immune infiltration and the prognosis of KIRC. erefore, our study was undertaken. e expression analysis revealed that ACE2 was downregulated in KIRC patients at the mRNA and protein level, and a low expression of ACE2 protein in KIRC patients was obtained in the subgroup analysis. ese results indicate that ACE2 may play a significant role in KIRC. Further prognosis analysis indicated that high ACE2 level in KIRC patients was linked to a favorable prognosis in both the TCGA and GEO cohorts, suggesting ACE2 could be a novel prognostic biomarker for KIRC and the prediction of a favorable outcome. ACE2 has also been suggested as a biomarker for other diseases or cancers. In thyroid carcinoma, ACE2 was employed as a biomarker and was also found to regulate tumor progression [13] . Moreover, ACE2 acted as a biomarker in chronic kidney disease and associated with higher risk for silent atherosclerosis [36] . A steady accumulation of data suggests that immune-cell infiltration could regulate tumor progression and metastasis, thus affecting the patients' prognosis [38, 39] . In our study, we also clarified the correlation between ACE2 and immune infiltration. We found ACE2 to be positively associated with the abundance of immune cells, including B cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells. Moreover, a strong correlation between the expression of ACE2 and most of the immune biomarker sets were analyzed. ese immune cells or biomarkers were known to be involved in tumor progression and as a biomarker for prognosis or therapy of KIRC. Ying et al. found that tumor microenvironment B cells were associated with poor survival and reduced response to treatment [40] . Moreover, CD39+ CD8+ T cells were shown to act as prognostic biomarkers in patients with KIRC and were also used to indicate poor prognosis [41] . In our study, we also found that ACE2 may affect the prognosis of KIRC patients, in part, due to immune infiltration. However, previous studies have suggested that T-regulatory cells are correlated with the poor outcomes of patients with KIRC [6] . However, we found that ACE2 positively associated with the abundance of several immune cells. us, all correlations of ACE2 and infiltrating immune cells in KIRC may not be favorable. Further study should be performed to verify these results and observations. e enrichment analysis suggested that ACE2 in KIRC were primarily involved in translation factor activity, immunoglobulin binding, metabolic pathways, transcriptional misregulation in cancer cells, cell cycle, and ribosomes. ese findings are consistent with the previous study where ACE2 was associated with cell cycle transcription [42] . Misregulation of ribosome functions and the cell cycle has 3 ACE2 TMEM27 GLYATL1 CLCN5 LRP2 SLC13A1 PHYHIPL SLC3A1 GBA3 CUBN FUT6 SLC27A2 HAO2 PKLR AFTPH PDZK1 GLYAT AGXT2 SLC16A9 SLCO4C1 TMEM174 FMO1 TINAG DDC SLC22A11 NAT8 SLC1A1 ACMSD RAB3IP LRRC19 SATB2 SLC5A10 EHHADH BBOX1 SLC7A9 SLC17A3 SLC22A24 SLC5A1 GIPC2 C11or154 MIOX NAT8B LOC15332B GJB1 CYP4A11 AVPR1B TAL2 SLC22A13 USH1C SLC22A12 (b) Figure 8 : Continued. C1or121 TPST2 MINA TRNP1 GPRC5A TMEM158 COL6A1 PCOLCE SMARCD3 GNAS GLT8D1 BALAP2L1 KCNG1 CD55 ETV4 DBN1 PPM1J ST5 GAD1 PSD3 PODNL1 SNX21 MAGED4 CACNB3 IMPDH1 SIX4 PCBP3 GPSM1 LMO1 MAGED4B TUBB3 MPP2 CREB3L1 LEF1 ITPKA Journal of Oncology been linked to many diseases, in particular cancers [43] . Our study further confirmed the association of ACE2 with transcriptional misregulation and the cell cycle in patients with KIRC. Genomic instability and mutagenesis could result in tumor genesis and progression while kinase could stabilize and repair genomic DNA [44] . In our study, several ACE2associated kinase targets included LCK, LYN, SYK, JAK3, and HCK. Interestingly, these kinases were mainly involved in regulating genomic stability, cell cycle, ribosomal activity, and immune infiltration [45, 46] . ese are consistent with the findings that ACE2 in KIRC were primarily involved in the cell cycle, ribosomal activity, and immune infiltration. Previous studies suggest that JAK3 is used as a biomarker and is associated with immune infiltration in patients with KIRC [7] . Another study found that alterations in the JAK3 pathway were involved in CD8 T cell immune deviation in RCC [47] . erefore, ACE2 may regulate immune infiltration via JAK3 kinase. However, further research needs to be done to confirm this. Transcriptional dysregulation and cell cycle disorder could lead to constant proliferation and abnormal cell invasion, which is the fundamental feature of cancers. Transcription factors are the key regulator of transcriptional activity and the cell cycle [48] . Our study identified several ACE2-associated transcription factor targets in KIRC, such as V$IRF_Q6, V$NFKB_Q6_01, V$ELF1_Q6, V$PEA3_Q6, and V$PU1_Q6. PEA3 could block cell cycle progression in breast cancer [49] . Moreover, PEA3 could facilitate cell invasion and metastasis in colorectal carcinoma [50] . erefore, ACE2 may regulate the cell cycle and progression of KIRC via these transcription factors. Further studies should be performed to test this hypothesis. In conclusion, our results indicate that ACE2 may be developed as a prognostic biomarker. It is associated with immune infiltration in patients with KIRC. is lays a foundation for further study of the function of ACE2 in the carcinogenesis and progression of KIRC. Data Availability e analyzed datasets generated during the study are available from the corresponding author on reasonable request. e authors declare that they have no conflicts of interest. Xinhao Niu performed data analysis work and aided in writing the manuscript. Juan Bao designed the study and assisted in writing the manuscript. Zhe Zhu and Enming Shao edited the manuscript. All authors read and approved the final manuscript. Xinhao Niu and Zhe Zhu contributed equally. Association between CD8+ T-cell infiltration and breast cancer survival in 12 439 patients CA: A Identification of therapeutic targets and prognostic biomarkers among CXC chemokines in the renal cell carcinoma microenvironment NR1B2 suppress kidney renal clear cell carcinoma (KIRC) progression by regulation of LATS 1/2-YAP signaling Identification of therapeutic targets and prognostic biomarkers among CXC chemokines in the renal cell carcinoma microenvironment Immune infiltration in renal cell carcinoma JAK3 is a potential biomarker and associated with immune infiltration in kidney renal clear cell carcinoma Treatment of renal cell carcinoma: current status and future directions e society for immunotherapy of cancer consensus statement on immunotherapy for the treatment of advanced renal cell carcinoma (RCC) ACE2 correlated with immune infiltration serves as a prognostic biomarker in endometrial carcinoma and renal papillary cell carcinoma: implication for COVID-19 ACE2 overexpression inhibits acquired platinum resistanceinduced tumor angiogenesis in NSCLC Downregulation of ACE2/ Ang-(1-7)/Mas axis promotes breast cancer metastasis by enhancing store-operated calcium entry Angiotensin converting enzymes ACE and ACE2 in thyroid cancer progression ACE2 and FZD1 are prognosis markers in squamous cell/adenosquamous carcinoma and adenocarcinoma of gallbladder Hypertension and hypertensive left ventricular hypertrophy are associated with ACE2 genetic polymorphism COVID-19: learning from lessons to guide treatment and prevention interventions e emerging COVID-19 research: dynamic and regularly updated science maps and analyses Angiotensin-converting enzyme 2 (ACE2), SARS-CoV-2 and the pathophysiology of coronavirus disease 2019 (COVID-19) Coronavirus disease (COVID-19): a machine learning bibliometric analysis ONCOMINE: a cancer microarray database and integrated data-mining platform UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses Proteomics. Tissue-based map of the human proteome GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses OSkirc: a web tool for identifying prognostic biomarkers in kidney renal clear cell carcinoma TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells Genome-wide association analysis identifies genetic correlates of immune infiltrates in solid tumors Gene expression markers of tumor infiltrating leukocytes e role of tumour-associated macrophages in bone metastasis Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal Link-edOmics: analyzing multi-omics data within and across 32 cancer types Clear cell sarcoma of the kidney: up-regulation of neural markers with activation of the sonic hedgehog and Akt pathways Gene signatures of progression and metastasis in renal cell cancer Extent and characteristics of immune infiltration in clear cell renal cell carcinoma and the prognostic value Overexpression of ACE2 produces antitumor effects via inhibition of angiogenesis and tumor cell invasion in vivo and in vitro e angiotensin-converting enzyme 2 in tumor growth and tumor-associated angiogenesis in non-small cell lung cancer Circulating angiotensin converting enzyme 2 activity as a biomarker of silent atherosclerosis in patients with chronic kidney disease Angiotensinconverting enzyme 2 attenuates the metastasis of non-small cell lung cancer through inhibition of epithelial-mesenchymal transition Bioinformatic identification of renal cell carcinoma microenvironmentassociated biomarkers with therapeutic and prognostic value Mining therapeutic and prognostic significance of STATs in renal cell carcinoma with bioinformatics analysis Identification and validation of dichotomous immune subtypes based on intratumoral immune cells infiltration in clear cell renal cell carcinoma patients Tumor-infiltrating CD39 (+) CD8 (+) T cells determine poor prognosis and immune evasion in clear cell renal cell carcinoma patients, Cancer immunology Ace2 receives helping hand for cell-cycle transcription New roles for Dicer in the nucleolus and its relevance to cancer Multiple functions of p21 in cell cycle, apoptosis and transcriptional regulation after DNA damage Exploration of the therapeutic aspects of Lck: a kinase target in inflammatory mediated pathological conditions SYK expression in human breast cancer JAK3/STAT5/6 pathway alterations are associated with immune deviation in CD8 T cells in renal cell carcinoma patients Comparative cell cycle transcriptomics reveals synchronization of developmental transcription factor networks in cancer cells Antitumoral actions of the anti-obesity drug orlistat (Xenical) in breast cancer cells: blockade of cell cycle progression, promotion of apoptotic cell death and PEA3-mediated transcriptional repression of Her2/neu (erbB-2) oncogene Pea3 expression promotes the invasive and metastatic potential of colorectal carcinoma