key: cord-0703764-n3itos4s
authors: He, Chaobin; Hua, Xin; Sun, Shuxin; Li, Shaolong; Wang, Jun; Huang, Xin
title: Integrated Bioinformatic Analysis of SARS-CoV-2 Infection Related Genes ACE2, BSG and TMPRSS2 in Aerodigestive Cancers
date: 2021-03-10
journal: J Inflamm Res
DOI: 10.2147/jir.s300127
sha: eafb2f4405a9e4425ccefba4fc4e48a3292c6f23
doc_id: 703764
cord_uid: n3itos4s

BACKGROUND: Cancer patients are more vulnerable to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection than the general population, with lung epithelial cells or enterocytes being the main targets. However, the expressions of SARS-CoV-2 entry-related genes in aerodigestive cancers have not been fully elucidated. METHODS: In this study, the expressions of SARS-CoV-2 receptors and cofactors, including angiotensin I-converting enzyme 2 (ACE2), basigin (BSG) and transmembrane serine protease 2 (TMPRSS2), were comprehensively assessed. We compared BSG and TMPRSS2 expressions between aerodigestive cancers and matched normal tissues through Gene Expression Profiling Interactive Analysis 2 (GEPIA2). Furthermore, expressions in healthy colon tissues at different anatomical locations were explored using the Genotype-Tissue Expression (GTEx) dataset. In addition, expressions among different tumor stages and the prognostic values were detected through GEPIA2. Moreover, the correlation between gene expression and immune infiltration was explored via Tumor Immune Estimation Resource (TIMER). Finally, expressions in primary colorectal cancer (CRC), lung metastasis and liver metastasis were investigated using the Gene Expression Omnibus (GEO) dataset GSE41258. RESULTS: Similar to ACE2, TMPRSS2 and BSG were also highly expressed in the digestive tracts. Intriguingly, BSG/TMPRSS2 expression in adjacent normal colon tissue or lung tissue was higher than that in corresponding healthy tissue, whereas they varied not among different tumor stages and correlated not with prognosis in aerodigestive cancers. Moreover, ACE2 was expressed at higher levels in lung metastases from CRC than in normal lung tissues. CONCLUSION: SARS-CoV-2 entry genes were highly expressed in CRC, and we reported for the first time higher expression of ACE2 in lung metastases from CRC than in normal lung, indicating that these patients may be more susceptible to extrapulmonary or pulmonary SARS-CoV-2 infection. Since our study is a bioinformatic analysis, further experimental evidences and clinical data are urgently needed.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of coronavirus disease 2019 (COVID- 19) , which has been declared to be a pandemic in 2020. As of 06 February 2021, 105,915,287 cases and 2,309,033 deaths have been reported globally (https://www.worldometers.info/coronavirus/). Angiotensin I-converting enzyme 2 (ACE2) is a functional receptor that plays a crucial role in SARS-CoV-2 infection, and transmembrane serine protease 2 (TMPRSS2) plays a main supporting role; 1-3 recently basigin (BSG) has been identified as another host receptor. 4, 5 Cancer patients are more susceptible to infections because of the immunosuppression caused by malignancy or anticancer treatments, 6 such as chemotherapy. Through bioinformatic analysis, many studies have indicated that cancer tissues expressed ACE2, BSG or TMPRSS2 at higher levels than matched healthy tissues in lung adenocarcinoma (LUAD), 2, 7, 8 stomach adenocarcinoma (STAD) and colon adenocarcinoma (COAD). 9, 10 This evidences provide the rationale for the high susceptibility of cancer patients to SARS-CoV-2. 11, 12 However, not all cancer patients are similarly at risk for COVID-19. 13, 14 SARS-CoV-2 is transmitted mainly through the respiratory or digestive route: 3 the lungs are in direct contact with the external environment, although the expression of ACE2 in lung tissues was generally lower than that in the digestive tract, pulmonary epithelial cells are the main targets of SARS-CoV-2. 3, 15, 16 Moreover, another study identified SARS-CoV-2 in resected specimens from colorectal cancer (CRC). 17 Gastrointestinal cancers mainly metastasize to the liver and to a lesser extent the lung; therefore, the susceptibility of patients with aerodigestive cancers, including lung metastases from CRC, is worth discussing.

The gene expressions in normal tissues were downloaded from the Genotype-Tissue Expression (GTEx), and the difference between men and women was explored using the R package ggpubr. 10 Expressions of BSG and TMPRSS2 in human cancer tissues or cell lines were downloaded from the Gene Expression Profiling Interactive Analysis 2 (GEPIA2; http://gepia2.cancer-pku. cn/) 18 or the Cancer Cell Line Encyclopedia (CCLE; https://portals.broadinstitute.org/ccle). 19 

Comparison of BSG and TMPRSS2 expressions between cancer and matched normal tissues was performed through GEPIA2, which has incorporated RNA-Seq data from The Cancer Genome Atlas (TCGA) and GTEx datasets. Likewise, we have compared gene expressions among COAD, adjacent normal colon and healthy colon; or lung cancer, adjacent normal lung and healthy lung. 10 As previously reported, 10,20 data were normalized in R environment using NormalizeBetweenArrays.

In addition, we explored the expression levels among different tumor stages using GEPIA2.

Tumor Immune Estimation Resource (TIMER; http:// timer.cistrome.org/) is adopted for the analysis of immune cell infiltration in cancer (including B cell, CD8 + T cell, CD4 + T cell, macrophage, neutrophil and dendritic cell), 21 as reported previously. 22 

OncoVar is an integrated database and analysis platform (https://oncovar.org/), by reanalyzing 10,769 exomes from 33 cancer types in TCGA and 1942 genomes from 18 cancer types in International Cancer Genome Consortium (ICGC), it could be used to identify driver mutations and driver genes. 23 Genes were classified into four grades from nonpathogenic to pathogenic: level 0 (non-pathogenic, score = 0), level 1 (possible pathogenic, score = 1), level 2 (likely pathogenic, 1 < score ≤ 10), level 3 (probable pathogenic, 10 < score < 20) and level 4 (pathogenic, score ≥ 20).

The relationship between BSG/TMPRSS2 expression and the prognosis, including overall survival (OS) and diseasefree survival (DFS), was analyzed using GEPIA2 by the Log-rank test.

Expression of SARS-CoV-2 Infection Related Genes in Metastases from CRC ACE2, BSG and TMPRSS2 expressions in primary CRC at different anatomical sites, lung metastases and liver metastases were explored using the TCGA or Gene Expression Omnibus (GEO) dataset GSE41258, with the exclusion of repeats and 28 metastases that exhibited high levels of stromal contamination. 24 

Comparisons were performed using Kruskal-Wallis tests or unpaired t-test (GraphPad Prism 7.0 software). The strength of the correlation was determined with the following guide: 0.00-0.19, very weak; 0.20-0.39, weak; 0.40--0.59, moderate; 0.60-0.79, strong; 0.80-1.0, very strong. 25 Two-side p < 0.05 was considered statistically significant. 

There have been reports about ACE2, BSG or TMPRSS2 expression in humans under physiological or pathological conditions. 5, 9, 10, 26 We have previously reported the high expression of ACE2 in normal or cancerous tissues of the digestive, urinary and male reproductive tracts; 10 herein we found that BSG and TMPRSS2 were also highly expressed in the normal digestive tract, and there was no big difference between men and women across different tissues ( Figure 1A and B). Likewise, gastrointestinal cancer tissues (from TCGA; Figure 1C and D) or cancer cell lines (from CCLE; Figure 1E and F) expressed BSG and TMPRSS2 at relatively high levels across different cancer types.

In the comparison between cancer tissues and adjacent normal tissues by using the Gene Expression Profiling Interactive Analysis 2 (GEPIA2) (from TCGA; Figure  S1 ), cancer tissues expressed BSG at a higher level than adjacent normal tissues in CHOL, ESCA, LIHC and STAD, but at a lower level in COAD (Figure S1A and C); the level of TMPRSS2 in COAD, READ, LUAD and LUSC was lower than that in matched adjacent normal tissues ( Figure S1B and D). After a combination of healthy tissues from the Genotype-Tissue Expression (GTEx) (Figure 2A -D), results showed that cancer tissues expressed BSG at a higher level than matched normal tissues in CHOL, ESCA, LIHC and pancreatic adenocarcinoma (PAAD); cancer tissues expressed TMPRS2 at a higher level in COAD and READ, but at a lower level in ESCA, PAAD, LUAD and LUSC. These results indicated a different expression pattern between adjacent normal tissues in TCGA and healthy tissues in GTEx, especially TMPRSS2 expression in the colon and rectum. It should be noted that the matched GTEx data for READ were from the healthy colon, therefore we compared the expression of BSG and TMPRSS2 among cancer, adjacent normal tissue and healthy tissue in the colon or lung. Results showed that adjacent normal colon expressed BSG and TMPRSS2 at different levels compared with a healthy colon, and transverse colon expressed TMPRSS2 at higher levels than sigmoid colon ( Figure 2E and F); BSG and TMPRSS2 expressions in the adjacent normal lung were generally higher than in healthy lung ( Figure 2G and H). In addition, we analyzed the correlation between BSG/TMPRSS2 expression and immune infiltration in COAD and lung cancers (LUAD and LUSC), and found a weak correlation between BSG/TMPRSS2 expression and immune cells infiltration (including B cell, CD4+ T cell, neutrophil and Dendritic cell; Figure 3A -F). Moreover, there were no statistically significant differences in the expression levels among different stages of aerodigestive cancers, except for ESCA, PAAD and LUAD for BSG ( Figure 4A -I); LIHC and PAAD for TMPRSS2 ( Figure 5A -I).

The Prognosis Value and Driver Gene Analysis of ACE2, BSG or TMPRSS2 in Aerodigestive Cancers

The prognosis of BSG and TMPRSS2 in aerodigestive cancers were explored. As shown in Figure 6 , BSG correlated negatively with OS and DFS in LIHC and LUAD ( Figure 6A-D) , whereas there was no significant correlation with prognosis (OS or PFS) in TMPRSS2 expression ( Figure 6E and F). In addition, we found that none of ACE2, BSG or TMPRSS2 was candidate cancer driver gene in the aerodigestive cancers (Table S1 ).

Since ACE2, BSG and TMPRSS2 were expressed at high levels in gastrointestinal cancers. First, we investigated their expression in CRC at different anatomical locations using TCGA data or GSE41258 dataset. The results showed that there was little difference in the expression of ACE2 ( Figure 7A and D), BSG ( Figure 7B and E) or TMPRSS2 ( Figure 7C and F). We next compared the expression levels among primary CRC, liver metastasis, normal liver tissue, lung metastasis and normal lung tissue. The comparative analysis between metastases and related normal tissues (the lung or liver) showed that ACE2 was expressed at a higher level in lung metastasis than in normal lung tissue, whereas the difference between liver metastasis and normal liver tissue did not reach statistical significance ( Figure 7G at a higher level than normal liver ( Figure 7H ); although TMPRSS2 expression in liver/lung metastases was higher than that in matched normal liver/lung, the differences did not reach statistical significance ( Figure 7I ). In addition, there were no statistically significant differences in ACE2, BSG or TMPRSS2 expression among primary CRC, liver metastasis and lung metastasis ( Figure 7G -I).

Patients with cancers are immunocompromised and vulnerable to severe infections, including SARS-CoV-2. For cancer patients with COVID-19, increased age, male sex and comorbidity were associated with mortality. 27, 28 First, SARS-CoV-2 invaded host cells through receptor ACE2 or BSG, and cofactors, such as TMPRSS2. 5, 26 

ACE2, BSG and TMPRSS2 were both highly expressed in human colon tissues, including normal and cancerous tissues; and there were no big differences between male and female ( Figure 1 ). 10 Therefore, the gender differences in morbidity and mortality may be attributed to the sex differences in immune. [29] [30] [31] Second, we found that both BSG and TMPRSS2 expression levels in adjacent normal tissue were generally higher than corresponding cancer or healthy tissue in colon or lung (Figure 2 ), which could be partially explained by the immune cells infiltration (Figure 3) , since TMPRSS2 correlated with immune cell infiltration and there were differences in immune infiltration among cancers, adjacent normal tissues and healthy tissues. [32] [33] [34] These findings were roughly consistent with that others have detected SARS-CoV-2 in enterocytes from adjacent normal colorectal tissues resected from CRC patients with COVID-19. 17 In addition, different anatomical sites of the colon expressed ACE2, BSG or TMPRSS2 at different levels ( Figure 2E and F). 10 However, they were not significantly (Figures 4 and 5) , and were not related to prognosis ( Figure 6 ) in COAD or READ, which was consistent with the negative results of OncoVar analysis (Table S1) .

Third, based on the abovementioned results, we aimed to explore the expression levels of SARS-CoV-2 entryrelated genes in lung metastases from CRC, since CRC expressed high levels of ACE2, BSG and TMPRSS2, whereas lung presented the most commonly affected site. Intriguingly, there was little difference in ACE2, BSG or TMPRSS2 expression among CRC from different anatomical sites ( Figure 7A -C, from TCGA; Figure 7D -F, from GSE41258), given that there were differences between normal transverse and sigmoid colon ( Figure 2E and F). 10 In addition, results showed that lung metastasis expressed ACE2, but not BSG or TMPRSS2 at higher levels than normal lung ( Figure 7G -I), since ACE2 was the most important gene for SARS-CoV-2 entry, our study suggested that patients with lung metastases could be more vulnerable to SARS-CoV-2 than the general population. Supportively, others have improved SARS-CoV-2 infection by transducing A549 cells (human adenocarcinomic alveolar basal epithelial) with a vector over-expressing human ACE2. 35 Moreover, the little difference among primary CRC, liver metastasis and lung metastasis indicated that ACE2, BSG or TMPRSS2 should be not relevant to the metastatic capacity, which was consistent with the finding that they were not significantly related to prognosis in COAD or READ (including OS and DFS; Figure 6 ). 10 In conjunction with our previous publication, 10 the results of this study showed that patients with aerodigestive cancers could be more susceptible to SARS-CoV-2; these findings were consistent with other clinical observations that lung cancer and CRC were the most common types in COVID-19 patients; [36] [37] [38] in addition, ACE2 expression in lung metastases from CRC was higher than that in the normal lung (Figure 7) , raising the possibility that these patients may be more vulnerable to pulmonary SARS-CoV -2 infection regardless of immunosuppression caused by malignancy or anticancer therapeutics, and lung metastases or circulating CRC cells could serve as a reservoir or transporter of SARS-CoV-2 for replication and dissemination. Moreover, the COVID-19 pandemic had a hysteresis effect on cancer diagnosis and treatment, and there was a report that SARS-CoV-2 infection could induce epithelial to mesenchymal transition in lung cancer cells, 39 which could induce cancer progression and metastasis.

Besides lung cells, SARS-CoV-2 can infect enterocytes, 15 cholangiocytes 40 and cardiomyocytes. 41 However, to this date, there is only limited evidence that SARS-CoV-2 could infect cancer cells. 39 Although studies have reported the morbidity of primary cancers to SARS-CoV-2, 11 we presented the rational of high susceptibility to SARS-CoV-2 in patients with aerodigestive cancers.

SARS-CoV-2 infection-related genes ACE2, BSG and TMPRSS2 were relatively highly expressed in aerodigestive cancers. To the best of our knowledge, this is the first study suggesting the genetic susceptibility of lung metastases to SARS-CoV-2. Patients with lung metastases of CRC should be more cautious to SARS-CoV-2, and be given more supervision. Since our data and other studies mainly depend on bioinformatics analysis, further validation with real-world or experimental data is needed. Gene Expression Omnibus; GEPIA2, Gene Expression Profiling Interactive Analysis 2; GTEx, Genotype-Tissue Expression; ICGC, International Cancer Genome Consortium; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; PAAD, pancreatic adenocarcinoma; READ, rectum adenocarcinoma; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2; STAD, stomach adenocarcinoma; TCGA, The Cancer Genome Atlas; TIMER, Tumor Immune Estimation Resource; TMPRSS2, transmembrane serine protease 2; TPM, transcripts per million.

All data are available in the manuscript or upon request to the corresponding author.

All authors report no conflicts of interest in this work.

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

Analysis of the susceptibility of lung cancer patients to SARS-CoV-2 infection

Pathophysiology, transmission, diagnosis, and treatment of Coronavirus disease 2019 (COVID-19): a review

SARS-CoV-2 invades host cells via a novel route: CD147-spike protein

The expression of SARS-CoV-2 receptor ACE2 and CD147, and protease TMPRSS2 in human and mouse brain cells and mouse brain tissues

Myeloid suppressor cells in cancer and autoimmunity

Expression of the SAR2-Cov-2 receptor ACE2 reveals the susceptibility of COVID-19 in non-small cell lung cancer

Structural variations and expression profiles of the SARS-CoV-2 host invasion genes in lung cancer

Pan-cancer analysis of transmembrane protease serine 2 and cathepsin L that mediate cellular SARS-CoV-2 infection leading to COVID-19

Bioinformatic analysis of correlation between immune infiltration and COVID-19 in cancer patients

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

Cancer associates with risk and severe events of COVID-19: a systematic review and meta-analysis

Are all patients with cancer at heightened risk for severe Coronavirus Disease

COVID-19 prevalence and mortality in patients with cancer and the effect of primary tumour subtype and patient demographics: a prospective cohort study

SARS-CoV-2 productively infects human gut enterocytes

Quantitative mRNA expression profiling of ACE 2, a novel homologue of angiotensin converting enzyme

Postoperative onset and detection of SARS-CoV-2 in surgically resected specimens from gastrointestinal cancer patients with pre/asymptomatic COVID-19

GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis

The cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity

Increased expression of TET3 predicts unfavorable prognosis in patients with ovarian cancer-a bioinformatics integrative analysis

TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells

Induced CD10 expression during monocyte-to-macrophage differentiation identifies a unique subset of macrophages in pancreatic ductal adenocarcinoma

OncoVar: an integrated database and analysis platform for oncogenic driver variants in cancers

Association of survival and disease progression with chromosomal instability: a genomic exploration of colorectal cancer

LAYN is a prognostic biomarker and correlated with immune infiltrates in gastric and colon cancers

Expression of SARS-CoV-2 entry receptors in the respiratory tract of healthy individuals, smokers and asthmatics

COVID-19 mortality in patients with cancer on chemotherapy or other anticancer treatments: a prospective cohort study

Clinical impact of COVID-19 on patients with cancer (CCC19): a cohort study

Covid-19: male disadvantage highlights the importance of sex disaggregated data

Sex differences in immune responses that underlie COVID-19 disease outcomes

Racial and gender-based differences in COVID-19. Front Public Health

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

Widespread parainflammation in human cancer

Understanding immune phenotypes in human gastric disease tissues by multiplexed immunohistochemistry

Imbalanced host response to SARS-CoV-2 drives development of COVID-19

Patients with cancer appear more vulnerable to SARS-CoV-2: a multicenter study during the COVID-19 outbreak

Outcomes of novel coronavirus disease 2019 (COVID-19) infection in 107 patients with cancer from Wuhan

Which type of cancer patients are more susceptible to the SARS-COX-2: evidence from a meta-analysis and bioinformatics analysis

SARS-CoV-2 infection induces EMT-like molecular changes, including ZEB1-mediated repression of the viral receptor ACE2, in lung cancer models

Recapitulation of SARS-CoV-2 infection and cholangiocyte damage with human liver ductal organoids

Human iPSC-derived cardiomyocytes are susceptible to SARS-CoV-2 infection

The Journal of Inflammation Research is an international, peerreviewed open-access journal that welcomes laboratory and clinical findings on the molecular basis, cell biology and pharmacology of inflammation including original research, reviews, symposium reports, hypothesis formation and commentaries on: acute/chronic inflammation; mediators of inflammation; cellular processes; molecular mechanisms; pharmacology and novel anti-inflammatory drugs; clinical conditions involving inflammation. The manuscript management system is completely online and includes a very quick and fair peerreview system. Visit http://www.dovepress.com/testimonials.php to read real quotes from published authors.