key: cord-0725943-pzu6qo4d
authors: Song, Jukun; Han, Jing; Liu, Feng; Chen, Xianlin; Qian, Shenqi; Wang, Yadong; Jia, Zhenyu; Duan, Xiaofeng; Zhang, Xiangyan; Zhu, Jianguo
title: Systematic Analysis of Coronavirus Disease 2019 (COVID-19) Receptor ACE2 in Malignant Tumors: Pan-Cancer Analysis
date: 2020-10-23
journal: Front Mol Biosci
DOI: 10.3389/fmolb.2020.569414
sha: 544a6acc9dcee353dd677fd15430104dae8d58d9
doc_id: 725943
cord_uid: pzu6qo4d

BACKGROUND: Coronavirus disease 2019 (COVID-19) was first detected in patients with pneumonia in December 2019 in China and it spread rapidly to the rest of the world becoming a global pandemic. Several observational studies have reported that cancer is a risk factor for COVID-19. On the other hand, ACE2, a receptor for the SARS-CoV-2 virus, was found to be aberrantly expressed in many tumors. However, the characterization of aberrant ACE2 expression in malignant tumors has not been elucidated. Here, we conducted a systematic analysis of the ACE2 expression profile across 31 types of tumors. METHODS: Distribution of ACE2 expression was analyzed using the GTEx, CCLE, TCGA pan-cancer databases. We evaluated the effect of ACE2 on clinical prognosis using the Kaplan-Meier survival plot and COX regression analysis. Correlation between ACE2 and immune infiltration levels was investigated in various cancer types. Additionally, the correlation between ACE2 and immune neoantigen, TMB, microsatellite instability, Mismatch Repair Genes (MMRs), HLA gene members, and DNA Methyltransferase (DNMT) was investigated. The frequency of ACE2 gene mutation in various tumors was analyzed. Functional enrichment analysis was conducted in various cancer types using the GSEA method. RESULTS: In normal tissues, ACE2 was highly expressed in almost all 31 organs tested. In cancer cell lines, the expression level of ACE2 was low to medium. Although aberrant expression was observed in most cancer types, high expression of ACE2 was not linked to OS, DFS, RFS, and DFI in most tumors in TCGA pan-cancer data. We found that ACE2 expression was significantly correlated with the infiltrating levels of macrophages and dendritic cells, CD4+ T cells, CD8+ T cells, and B cells in multiple tumors. A positive correlation between ACE2 expression and immune neoantigen, TMB, and microsatellite instability was found in multiple cancers. GSEA analysis which was carried out to determine the effect of ACE2 on tumors indicated that several cancer-associated pathways and immune-related pathways were hyperactivated in the high ACE2 expression group of most tumors. CONCLUSION: These findings suggest that ACE2 is not correlated with prognosis in most cancer types. However, elevated ACE2 is significantly correlated with immune infiltrating levels, including those of CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and DCs in multiple cancers, especially in lung and breast cancer patients. These findings suggest that ACE2 may affect the tumor environment in cancer patients with COVID-19.

Accumulating evidence demonstrates that the Spike (S) protein of SARS-CoV-2 binds Angiotensin-Converting Enzyme 2 (ACE2) to trigger COVID-19 (Hoffmann et al., 2020) . Besides, the expression of ACE2 is a key determinant for the entry of SARS-CoV-2 into host cells (Hamming et al., 2004; Kuba et al., 2005) . Several studies report that the expression and distribution of ACE2 are enriched in the lungs (Lukassen et al., 2020) , esophagus (Qi et al., 2020) , kidneys (Pan et al., 2020) , bladder (Zou et al., 2020) , testis , stomach (Qi et al., 2020) , ileum (Zou et al., 2020) , and oral mucosa using a single-cell RNA sequencing technique. These findings suggested that organs with high ACE2 expression are potentially high risk for SARS-CoV-2 infection. Cancer patients in a state of systemic immunosuppression, are considered to be highly vulnerable to the COVID-19 epidemic . Current clinical studies show that COVID-19 patients with LUSC tumors have a higher risk of clinical complications and death than patients without cancer (Liang et al., 2020; Yu et al., 2020; Zhang et al., 2020) . However, a systematic analysis of the aberrant expression of ACE2 in human cancer has not yet been conducted. Therefore, in this study a bioinformatics approach was used to evaluate the prognosis of the ACE2 in TCGA Pan-cancer data. Furthermore, the association between ACE2 expression and immune neoantigen, TMB, microsatellite instability, and HLA family members was investigated.

Gene expression profiles were downloaded from the three publicly available datasets, the TCGA Pan-cancer cohort 1 , Genotype-Tissue Expression (GTEx) project 2 , and Broad 1 https://portal.gdc.cancer.gov/ 2 https://commonfund.nih.gov/GTEx/ Institute Cancer Cell Line Encyclopedia (CCLE) 3 . mRNA data in normal tissues were obtained from the GTEx project, which included 31 tissues. Distribution of the expression levels of cancer cell lines in 21 organizations was also conducted. The mRNA data of adjacent tumor tissues and tissue in 31 types of tumors were obtained from the TCGA dataset. Clinical information was also obtained (Liu et al., 2018) and the Kruskal-Wallis test was done to determine the differences among organs.

To determine the difference between normal healthy tissues, adjacent tumor samples, and tumor samples, TCGA Pan-cancer and GTEx datasets were downloaded from the UCSC XENA dataset 4 . The tumor samples and adjacent samples were obtained from the TCGA dataset while the normal samples were obtained from the GTEx dataset. The difference in ACE expression was 

To evaluate whether the ACE2 expression level was associated with tumor prognosis in various cancers, univariate COX regression analysis was performed for Overall Survival (OS), Disease-Free Survival (DFS), Disease-Specific Survival (DSS), Disease-Free Interval (DFI), and Progression-Free Interval (PFI). The threshold was adjusted to a Cox P < 0.05. Based on the median of ACE2 expression levels, samples were divided into two groups: the high expression and the low expression group. Kaplan-Meier survival analysis was used to compare the differences in 31 types of tumors. A log-rank test was used to calculate the significance of survival differences with a threshold of P < 0.05. We also evaluated the association between ACE2 expression level and other clinical features, including the TNM stage.

Tumor-infiltrating lymphocytes are independent predictors of cancer sentinel lymph node status and survival (Jin and Hu, 2020; Koletsa et al., 2020) . We investigated whether the expression of ACE2 is related to the levels of immune infiltration in different types of cancers. The CIBERSORT method was employed to evaluate the relative proportion of 22 immune cell infiltrations across multiple cancers (Newman et al., 2015) . Meanwhile, the ESTIMATE method was used to assess the immune cell infiltration levels including the immune score, tumor purity, and stromal score for each tumor sample in the TCGA pan-cohort (Yoshihara et al., 2013) . The correlation between ACE2 and scores of these immune cells in 31 types of cancers was analyzed using the Spearman correlation method. Under normal circumstances, the immune system can recognize and remove tumor cells in the tumor microenvironment. Tumor immunotherapy is a treatment method to control and eliminate tumor cells by restarting and maintaining the tumor's immune cycle thus restoring the normal anti-tumor immune response in the body. The immune checkpoint genes include monoclonal antibody class immune checkpoint inhibitors, therapeutic antibodies, cancer vaccines, cell therapy, and small molecule inhibitors (Danilova et al., 2019) . We collected more than forty immune checkpoint genes and analyzed the relationship between ACE2 gene expression and the expression of immune checkpoint genes using the Spearman correlation analysis.

Neoantigen is a neonatal antigen that is encoded by a mutant gene of a tumor cell. Different from the proteins expressed by normal cells, it is a new abnormal protein, mainly resulting from a point mutation, deletion mutation, or gene fusion (Srivastava et al., 2020) . The peptide fragments formed from these proteins after enzymatic hydrolysis are presented as antigens by DC cells to T cells. DC cells promote the activation of T cells to mature T cells that specifically recognize new tumor antigens and also promote the proliferation of the activated T cells (Rooney et al., 2015) .

Neoantigen vaccines can be developed using the immune activity of tumor neoantigens, according to a specific mutation in the tumor cells, then administered to patients to achieve the required therapeutic effect. Therefore, the number of neonatal antigens in each tumor sample was imputed and analyzed. The correlation between ACE2 expression and the gene markers of immune neoantigen was analyzed using the Spearman correlation method.

The Tumor Mutational Burden (TMB) is usually measured by the number of somatic mutations that occur within an average of 1 Mb in the coding region (exon region) of the tumor cell genome (non-synonymous mutations) (Krieger et al., 2020) . The total number of synonymous mutations indicate that the mutation pattern includes single nucleotide mutation (SNV) and small fragment insertion/deletion (Indel) and other forms of mutation. TMB is used to reflect the number of mutations in tumor cells and is a quantifiable biomarker. Here, we counted the TMB of each tumor sample separately and analyzed the relationship between gene expression and TMB using the Spearman rank correlation coefficient.

Microsatellite Instability (MSI) refers to any change in the length of a microsatellite due to the insertion or deletion of a repeating unit in a tumor compared to normal tissue. The appearance of a new microsatellite allele is a genetic phenomenon (Bonneville et al., 2017) . We analyzed the correlation between gene expression and MSI using the Spearman rank correlation coefficient. Disparities in the expression of HLA (Human Leukocyte Antigen) family members on the surface of tumor cells are part of the early and common activities that promote carcinogenesis. This is because HLA is essential for the immune recognition of tumor cells and signal transmission between tumor and immune cells (Campoli and Ferrone, 2008; Sim and Sun, 2017) . Thus, we also investigated the correlation between ACE2 expression and HLA genes.

We downloaded mutation data of multiple malignancies from the TCGA database and analyzed the alterations of the ACE2 gene within these tumors. We visualized the tumor with the most ACE2 mutations using an R data package, maftools (Mayakonda and Koeffler, 2016) .

Mismatch Repair Genes (MMRs) are the key players in the intracellular mismatch repair mechanism (Meiser et al., 2020) .

Loss of a key gene function in this mechanism can lead to DNA replication errors that cannot be repaired and in turn, lead to more somatic mutations. Therefore, TCGA expression profiling was used to evaluate the relationship between ACE2 and five MMRs genes: MLH1, MSH2, MSH6, PMS2, and EPCAM mutations.

DNA methylation is a chemical modification of DNA, which changes the genetic activity without changing the DNA sequence. DNA methylation causes changes in the chromatin structure, DNA conformation, DNA stability, and the proteins-DNA interaction, thereby regulating gene expression. In DNA methylation, a methyl group is covalently bound to the 5 carbon of cytosine, of the genomic CpG dinucleotide, under the catalysis of DNA methyltransferase . Here, we analyzed the correlation between gene expression and the expression of four methyltransferases (DNMT1, DNMT2, DNMT3A, and DNMT3B).

To determine the effect of ACE2 expression on tumors, we divided the samples into high and low expression groups according to their gene expression level and used Gene Set Enrichment Analysis (GSEA) (Barber et al., 2006) to analyze the enrichment of KEGG and HALLMARK pathways in the high and low expression groups. The gene sets using the c5 curated signatures were downloaded from the Molecular Signature Database (MSigDB) of the Broad Institute. The KEGG terms were defined between the high ACE2 expression group and the low ACE2 expression group. The significant enrichment of pathways was determined based on FDR < 0.05.

We analyzed the distribution of ACE2 in the GTEx dataset, and the results indicated that ACE2 is highly expressed in most of the organs tested. The highest expression levels were noted in the small intestine, testis, and thyroid while the lowest expression levels were in the blood, brain, and breast ( Figure 1A) . The Kruskal-Wallis test indicated that there were significant differences in the ACE2 expression levels among the organs.

The distribution of ACE2 in cancer cell lines was analyzed, and our findings showed that the expression level of ACE2 in cancer cell lines was low to medium in different tissues. The highest expression levels were noted in the intestine, salivary gland, and upper aerodigestive tract while the lowest expression levels were in the bone, pleura, and soft tissue ( Figure 1B) . The Kruskal-Wallis test indicated that there were significant differences in the ACE2 expression levels among the organs.

The distribution of ACE2 expression in the TCGA pan-cancer dataset is shown in Figure 2A . Differences in the ACE2 expression level between the primary tumor and adjacent tumor samples were observed in multiple tumor types. ACE2 expression was upregulated in seven tumors types (CHOL, GBM, KIRP,

LGG, LUAD, LUSC, READ, and UCEC) and downregulated in six tumors types (BRCA, KICH, LIHC, PRAD, STAD, and THCA). However, no significant difference was observed in ACE2 expression in BLCA, COAD, ESCA, HNSC, KIRC, LIHC, PAAD, and STAD.

Considering that there are relatively few tumor samples in the TCGA database, we combined the data of normal tissues from the GTEx database with the data of TCGA tumor tissues to analyze the ACE2 expression differences in 31 types of tumors. Our results indicate that there were significant differences in ACE2 expression in multiple tumors ( Figure 2B ). This analysis showed that ACE2 expression was higher in colorectal, gastric, kidney, lung, esophagus, pancreatic cancers, and lymphoma tumors compared to the normal tissues. Lower ACE2 expression was also observed in the breast, testis, thyroid, skin, ovary, and prostate tumors in some data sets.

We investigated whether ACE2 expression was linked with Overall Survival (OS) in cancer patients. Using univariate COX regression analysis, the forest plot in multiple cancer types is shown in Figure 3A and Table 1 . The results demonstrated that ACE2 has less prognostic influence on most cancers, except for KIRC (OS HR = 0.83, CI = 0.77-0.89, P < 0.001). Kaplan-Meier survival analysis was performed to screen for prognostic tumor types. The prognostic Kaplan-Meier survival curve is shown in Figure 3B . Higher ACE2 expression was associated with improved OS in KIRC. The results were linked with the COX results, which also indicated that ACE2 expression was not related to prognosis in most tumors.

Considering that there may be non-tumor factors leading to death during the follow-up period, we further analyzed the relationship between ACE2 expression and prognosis Disease-Specific Survival (DSS) among 31 types of tumors. The results indicated that ACE2 was not correlated with DSS in most tumors, except for KIRC (DDS HR = 0.65, CI = 0.65-0.78, P < 0.001), LGG (DDS HR = 11.97, CI = 11.97-153.64, Table 2 ). Kaplan-Meier survival analysis showed that there were significant differences in ACE2 expression in eight types of cancers (KIRC, KIRP, LGG, LIHC, LUSC, OV, UCS, and UVM) ( Figure 4B ). Among these cancers, higher ACE2 was associated with better DSS, except for LGG. Furthermore, we analyzed the correlation between ACE2 expression and prognosis Disease-Free Interval (DFI) in 31 types of tumors. The forest plot is shown in Figure 5A and Table 3 . In most tumors except for OV (DFI HR = 0.53, CI = 0.29-0.95, P < 0.001), ACE2 was not associated with DFI. Kaplan-Meier survival analysis showed that there were significant differences in ACE2 expression in ACC, LIHC, and OV ( Figure 5B) .

Additionally, we analyzed the relationship between ACE2 expression and Progression-Free Interval (PFI) in multiple cancer types. The summarized forest plot is shown in Figure 6A and Table 4 . In most tumors ACE2 was not associated with PFI except for KIRC (PFI HR = 0.81, CI = 0.75-0.87, P < 0.001),

LGG (PFI HR = 6.64, CI = 1.74-25.28, P < 0.001), SARC (PFI HR = 0.54, CI = 0.31-0.92, P < 0.001), LUSC (PFI HR = 0.77, CI = 0.61-0.97, P < 0.001), and MESO (PFI HR = 0.05, CI = 0.003-0.873, P < 0.001). Kaplan-Meier survival analysis showed that there were significant differences in ACE2 expression in COAD, KIRC, LIHI, LUSC, OV, and USC ( Figure 6B) .

The association of ACE2 with the TNM stage across 31 types of tumors was performed. Results showed no correlation between ACE2 and the TNM stage in all tumors (Supplementary Figure 1) .

We employed the CIBERSORT algorithm to determine the relative fractions of immune cell infiltration in various tumors. We found that ACE2 expression was correlated with immune infiltration levels in different types of cancers. ACE2 expression was significantly correlated with: macrophages infiltration levels in 13 types of cancers; dendritic cells in 12 types of cancers; CD4+ T cells in eight types of cancers; CD8+ T cells in three types of cancers; mast cells in eight types of cancers; B cells in five types of cancers; and NK cells in two types of cancers (Figure 7) . For example, ACE2 expression level in LUSC was linked with infiltration levels of dendritic cells (r = 0.27, P = 1.4e-09), macrophages M0 (r = -0.18, P = 6.3e-05), mast cells resting (r = 0.16, P = 0.00027), and neutrophils (r = 0.16, P = 0.00032).

Given the correlation of ACE2 expression with immune infiltration levels in different types of cancers, we analyzed the association between ACE2 and tumor purity. We found that ACE2 expression positively correlates with immune infiltration in BLCA and BRCA, LGG, SARC, TGCT, THYM, while a negative association was observed in COAD and KIRC, LIHC, SKCM FIGURE 9 | A heatmap showing the relationship between ACE2 and the expression of immune checkpoint genes. *P < 0.05, **P < 0.01; ***P < 0.001. (Figure 8) . These findings strongly suggest that ACE2 plays a vital role in immune infiltration in several tumors.

Tumor immunotherapy is a recent approach to treating malignant tumors. We collected more than forty immune checkpoint genes and analyzed their relationship with ACE2 expression. A correlation heatmap showing the association of the immune checkpoint genes and ACE2 expression is shown in Figure 9 . We found that ACE2 expression correlates with the immune checkpoint gene in PRAD, COAD, and UCEC.

We regularly collect and synthesize neoantigen vaccines, according to mutations in specific tumor cells, and using the immune activity of tumor neoantigens we immunize patients to achieve the required therapeutic effect. Therefore, we counted the number of new antigens in each tumor sample and analyzed the relationship between ACE2 expression and the number of antigens. The ACE2 expression level was positively associated with immune neoantigen in GBM, LGG, BRCA, KIRC, and KIRP while in the LUAD, COAD, LUSC, and UCEC, a negative relationship was observed (Figure 10) .

TMB reflects the number of mutations in tumor cells and is a quantifiable biomarker. Therefore, the TMB of each tumor sample was counted separately and the relationship between ACE2 expression and TMB was analyzed using the Spearman rank correlation coefficient. Significant correlations were found between ACE2 expression and TMB in BRCA, STAD, SKCM, COAD, SARC, HNSC, KIRC, KIRP, LAML, LGG, LUAD, and PAAD ( Figure 11A) . These results indicate that there is a strong relationship between ACE2 and TMB in most tumors.

Microsatellite instability (MSI) is a pattern of hypermutation that occurs in genomic microsatellites and is caused by defects in the mismatch repair system. The correlation between ACE2 and MSI was analyzed using the Spearman rank correlation coefficient. Results showed a significant correlation between ACE2 expression and MSI in BRCA, CESC, DLBC, THCA, SKCM, HNSC, and PRAD ( Figure 11B) . These results revealed a strong relationship between ACE2 and MSI in most of the tumors. Notably, many HLA genes showed a significant correlation with ACE2 expression in most tumors ( Figure 11C ).

The mutation data of 31 tumors were downloaded from TCGA, and we analyzed the ACE2 mutation in these tumors. A schematic diagram of the tumor with the highest number of mutations is shown in Figure 12 . Notably, ACE2 mutations were only observed in BLCA, BRCA, COAD, HNSC, LAML, LGG, LUAD, LUSC, OV, SKCM, STAD, and UCEC. The highest mutation was in UCEC. These findings indicate that ACE2 mutations seldom occur in most tumors.

Mismatch Repair genes are intracellular mismatch repair mechanisms. Mismatch repair deficiency is known to occur in some tumors. TCGA expression profiling was used to evaluate the association between the five key MMRs genes MLH1, MSH2, MSH6, PMS2, EPCAM, and ACE2. The heatmap is shown in Figure 13A . We found that ACE2 was significantly associated with five key MMRs in COAD, KIRC, LGG, TGCT, and THCA.

DNA methyltransferase maintains DNA methylation by transferring methyl groups to the cytosine of CpG dinucleotide islands. DNA hypermethylation in the promoter region of a tumor suppressor gene leads to gene silencing, which in turn leads to dysregulation of multiple signaling pathways associated with human malignancy. The correlation between ACE2 and the expression of four key DNA methyltransferases is shown in Figure 13B . It can be inferred that ACE2 was significantly associated with four key DNMT in THCA, THYM, STAD, PRAD, LUSC, LGG, KIRP, and COAD.

To determine the effect of ACE2 on tumors, we divided the samples into high and low groups according to the median of ACE2 expression levels and used GSEA to analyze the enrichment of the KEGG pathways in the high and low expression groups. The first three most significant pathways are shown in Figure 14 . The GSEA analysis demonstrated that several cancer-associated pathways and immune-related pathways were hyperactivated in the high ACE2 expression group in most tumors. RIG I LIKE RECEPTOR SIGNALING PATHWAY and the TOLL LIKE RECEPTOR SIGNALING PATHWAY were found in BLCA, ACC, LUAD, and PCPG. REGULATION_OF_AUTOPHAGY was found in ACC, BRCA, GBM, LUSC, and PAAD. NATURAL KILLER CELL MEDIATED FIGURE 11 | Correlation between ACE2 expression level and TMB (A) and Microsatellite instability (B), HLA genes (C). *P < 0.05, **P < 0.01; ***P < 0.001.

CYTOTOXICIT, and T CELL RECEPTOR SIGNALING PATHWAY was found in THCA, GBM, and OV.

The first public reports of SARS-CoV-2 infection, which causes a severe acute respiratory syndrome, were in December 2019. Subsequently, evidence of human-to-human transmission was released in January 2020 . COVID-19 has since caused a global pandemic and the WHO has declared the novel coronavirus disease an international public health emergency. Cancer patients are considered to be highly vulnerable to COVID-19 infection (Liang et al., 2020; Onder et al., 2020) . The patient populations with ongoing and previous cancer history represent a unique subset that is more susceptible to SARS-CoV-2 infection, and many cancer patients currently undergoing treatment are greatly concerned about SARS-CoV-2 infection. Recent studies with small cohorts (n = 12-28) of SARS-CoV-2 infected cancer patients, have proposed a potentially higher susceptibility and a poorer prognosis for COVID-19 in cancer patients compared to the general population (Liang et al., 2020; FIGURE 12 | Identification of gene mutation patterns in several cancers. Yu et al., 2020; Zhang et al., 2020) . A multicenter retrospective cohort study also reported that patients with malignant tumors have a higher risk of developing severe COVID-19 disease (Tian et al., 2020) . Moreover, another retrospective study reported that women who received an active treatment for gynecological cancer are at higher risk of severe COVID-19 (Bogani et al., 2020) . Yet another study demonstrated that the risk of COVID-19 might outweigh the urgency for cancer treatment (Indini et al., 2020) . These findings have indicated that COVID-19 may affect the prognosis of cancer patients. However, there is currently no systematic analysis framework to characterize abnormal ACE2 expression in human cancers.

The Renin-angiotensin system (RAS) is recognized as one of the most important regulators of normal and pathophysiological processes in the brain, kidney, heart, and blood vessels. The activation of the Angiotensin-Converting Enzyme 2/Angiotensin-(1-7)/Mitochondrial Assembly Receptor [ACE2/Ang-(1-7)/MasR] axis is a component of RAS and has recently been identified as an important part of gastric mucosa and cancer. ACE2 appears to play a central role in the formation of Ang-(1-7). ACE2 was linked with fibrosis, inflammation, cancer, angiogenesis (Passos-Silva et al., 2013) , and also plays an important role in COVID-19 infection. During viral infection, a key step in the process is the virus entry into the host cells and ACE2-expressing cells are target cells that are susceptible to SARS-CoV-2 infection (Parodi and Liu, 2020 ). ACE2 appears to play a central role in the ACE2/Ang-(1-7)/MasR axis. A previous meta-analysis has suggested that the ACE2/Ang-(1-7)/MasR axis plays an important role in the development of tumors (de Paula Gonzaga et al., 2020) . It is worth noting that ACE2 has both positive and negative roles in tumorigenesis (Xu et al., 2017) . On the one hand, ACE2 suppresses breast cancer angiogenesis by inhibiting the VEGF/VEGFR2/ERK pathway (Zhang et al., 2019) and on the other hand, downregulation of the ACE2/Ang-(1-7)/MasR axis enhances breast cancer metastasis by enhancing store-manipulated calcium entry (Yu et al., 2016) . Thus, the role of the ACE2/Ang-(1-7)/MasR axis in cancer is complicated.

In this work, a pan-cancer analysis of ACE2 in malignancies in the TCGA dataset was performed. The expression levels of ACE2 and systematic prognostic landscape in different types of cancers were examined using TCGA data. We observed the differential expression of ACE2 between cancer and normal tissues. We found that, compared to normal tissues, ACE2 was highly expressed in colorectal, gastric, kidney, lung, pancreatic cancers, and lymphoma tumors. Some data sets also showed that ACE2 had a lower level of expression in the breast, uterus, brain, skin, ovarian, and prostate cancers. In this work, ACE2 overexpression was found in lung cancer, hepatocellular carcinoma, and gallbladder cancer, which is in line with a previous study (Xu et al., 2017) . However, these findings contradict previous reports on ACE2 expression in breast cancer, pancreatic cancer, and oral squamous cell carcinoma (de Carvalho Fraga et al., 2017; Xu et al., 2017) . The factors that contribute to these contradictory effects of the ACE2/Ang-(1-7)/MasR axis on cancer require additional investigation.

Although ACE2 was aberrantly expressed in most tumors, its high expression was not associated with OS, DFI, PFI, and DSS, but was associated with KIRC, KIRP, LGG, LIHC, LUSC, OV, UCS, UVM, ACC, COAD, and OV. Interestingly, the high expression level of ACE2 correlated with improved prognosis of KIRC, KIRP, LIHC, LUSC, OV, UCS, UVM, ACC, COAD, and OV, but not LGG. Some reports show that ACE2 may have both positive and negative effects on the development of cancer. It has been demonstrated that overexpressed ACE2 may inhibit cell growth and vascular endothelial growth factor production in lung cancer (Cheng et al., 2016) , breast cancer (Yu et al., 2016) , colon cancer (Bernardi et al., 2012) , and pancreatic cancer (Zhou et al., 2011) . On the other hand, overexpressed ACE2 may promote the migration and invasion of human renal carcinoma cells (Zheng et al., 2015; Errarte et al., 2017) . The role of ACE2 in cancer development is, therefore, complicated. It is well-known that the interaction of SARS-CoV-2 with RAS through ACE2 is a key factor in infection. Once SARS-CoV-2 enters the target cell, the host response will be the main determinant of the severity of the subsequent pathogenesis (Gosain et al., 2020; Vaduganathan et al., 2020) . Innate lymphocytes including the invariant T (MAIT) cells and γδ Tcells respond quickly to the pathogen invasion and trigger a cytokine response that is essential for the killing of the microorganisms (Hotchkiss and Opal, 2020) . After the initial innate response, a specific adaptive immune response is required to eliminate SARS-CoV-2 (Extermann, 2000) . Lymphopenia, an independent indicator of poor prognosis in COVID-19 patients, is common and compromises the required immune response (Tan et al., 2020) . Long-lasting cytokine release, mediated by more white blood cells than T lymphocytes, may cause a "cytokine storm" and cause severe lung damage (Shi et al., 2020) . Hence, another important aspect of this study is that ACE2 expression is correlated with diverse immune infiltration levels in cancer, especially in gastric and colon cancers. Our analyses further demonstrate that ACE2 expression has a significant correlation with the infiltrating levels of dendritic cells, CD4+ T cells, CD8+ T cells, mast cells, B cells, and NK cells in most cancers. ACE2 expression also has a significantly positive correlation with the infiltrating levels of dendritic cells, macrophages M0, mast cells resting, and neutrophils in multiple cancers. ACE2 expression was linked with immune neoantigen, TMB, and microsatellite instability, especially in BRCA and SKCM. Significant correlations were found between ACE2 expression and mismatch repair genes and also with DNA methyltransferases. Thus, our study provides significant insight into the potential role of ACE2 in tumor immunology.

We found that increased ACE2 expression correlates with immune infiltration levels in most tumors and that higher ACE2 expression was related to immune neoantigen, TMB, and microsatellite instability. These findings indicate that although ACE2 is not associated with prognosis, in most cancers, it is significantly associated with immune penetration levels in a variety of tumors (especially patients with lung and breast cancer) including CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and DC. Our findings confirm that ACE2 affects the immune infiltration of cancer patients with COVID-19. Since this study is only a database analysis, it is necessary to verify the results in a larger clinical cohort.

The datasets generated for this study can be found in the online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.

JS, JH, FL, XC, JZ, and XD wrote the main manuscript text. JS and WD prepared Figures 1-14 . SQ and JS contributed to data analysis. All authors reviewed the manuscript. 

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb. 2020.569414/full#supplementary-material

Restoring function in exhausted CD8 T cells during chronic viral infection

Characterization and significance of ACE2 and Mas receptor in human colon adenocarcinoma

Cancer patients affected by COVID-19: experience from milan, lombardy

Landscape of microsatellite instability across 39 cancer types

HLA antigen changes in malignant cells: epigenetic mechanisms and biologic significance

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

ACE2 overexpression inhibits acquired platinum resistance-induced tumor angiogenesis in NSCLC

Programmed Cell Death Ligand-1 (PD-L1)

Expression Profiling Identify an Immunologic Subtype of Pancreatic Ductal Adenocarcinomas with Favorable Survival

Angiotensin-converting enzymes (ACE and ACE2) as potential targets for malignant epithelial neoplasia: review and bioinformatics analyses focused in oral squamous cell carcinoma

ACE2/Angiotensin-(1-7)/Mas receptor axis in human cancer: potential role for pediatric tumors

Expression and activity of angiotensin-regulating enzymes is associated with prognostic outcome in clear cell renal cell carcinoma patients

Note from the editors: world Health Organization declares novel coronavirus (2019-nCoV) sixth public health emergency of international concern

Measuring comorbidity in older cancer patients

COVID-19 and cancer: a comprehensive review

Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis

SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor

Activating immunity to fight a foe -a new path

Developing a risk assessment score for patients with cancer during the coronavirus disease 2019 pandemic

Tumor-infiltrating CD8 T cells predict clinical breast cancer outcomes in young women

Prognostic impact of stromal and intratumoral CD3, CD8 and FOXP3 in adjuvantly treated breast cancer: do they add information over stromal tumor-infiltrating lymphocyte density?

Targeted literature review on use of tumor mutational burden status and programmed cell death ligand 1 expression to predict outcomes of checkpoint inhibitor treatment

A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury

Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia

Cancer patients in SARS-CoV-2 infection: a nationwide analysis in China

An Integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics

SARS-CoV-2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells

Maftools: efficient analysis, visualization and summarization of MAF files from large-scale cohort based cancer studies. bioRxiv

Evaluation of implementation of risk management guidelines for carriers of pathogenic variants in mismatch repair genes: a nationwide audit of familial cancer clinics

Robust enumeration of cell subsets from tissue expression profiles

Case-Fatality rate and characteristics of patients dying in relation to COVID-19 in Italy

Identification of a potential mechanism of acute kidney injury during the COVID-19 outbreak: a study based on single-cell transcriptome analysis

From containment to mitigation of COVID-19 in the US

Angiotensin-(1-7): beyond the cardio-renal actions

Single cell RNA sequencing of 13 human tissues identify cell types and receptors of human coronaviruses

Molecular and genetic properties of tumors associated with local immune cytolytic activity

COVID-19 infection: the perspectives on immune responses

HLA-F: a new kid licensed for peptide presentation

Diverse neoantigens and the development of cancer therapies

Lymphopenia predicts disease severity of COVID-19: a descriptive and predictive study

Clinical characteristics and risk factors associated with COVID-19 disease severity in patients with cancer in Wuhan, China: a multicentre, retrospective, cohort study

Renin-angiotensin-aldosterone system inhibitors in patients with Covid-19

scRNA-seq profiling of human testes reveals the presence of the ACE2 Receptor, a target for SARS-CoV-2 infection in spermatogonia, leydig and sertoli cells

High expression of ACE2 receptor of 2019-nCoV on the epithelial cells of oral mucosa

Clinical findings in a group of patients infected with the 2019 novel coronavirus (SARS-Cov-2) outside of Wuhan, China: retrospective case series

The ACE2/Angiotensin-(1-7)/Mas receptor axis: pleiotropic roles in Cancer

Clinical characteristics of COVID-19 after gynecologic oncology surgery in three women: a retrospective review of medical records

Inferring tumour purity and stromal and immune cell admixture from expression data

Downregulation of ACE2/Ang-(1-7)/Mas axis promotes breast cancer metastasis by enhancing store-operated calcium entry

DNA methyltransferases: emerging targets for the discovery of inhibitors as potent anticancer drugs

SARS-CoV-2 transmission in patients with cancer at a tertiary care hospital in Wuhan, China

Clinical characteristics of COVID-19-infected cancer patients: a retrospective case study in three hospitals within Wuhan

ACE2 inhibits breast cancer angiogenesis via suppressing the VEGFa/VEGFR2/ERK pathway

Ang-(1-7) promotes the migration and invasion of human renal cell carcinoma cells via Mas-mediated AKT signaling pathway

Angiotensin-converting enzyme 2 acts as a potential molecular target for pancreatic cancer therapy

Single-cell RNAseq data analysis on the receptor ACE2 expression reveals the potential risk of different human organs vulnerable to 2019-nCoV infection