key: cord-299989-p59u6qa0
authors: Zhang, Lei; Han, Xiaohong; Shi, Yuankai
title: Comparative analysis of SARS-CoV-2 receptor ACE2 expression in multiple solid tumors and matched non-diseased tissues
date: 2020-06-18
journal: Infect Genet Evol
DOI: 10.1016/j.meegid.2020.104428
sha: 
doc_id: 299989
cord_uid: p59u6qa0

The emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a global public health emergency. SARS-CoV-2 employs the host cell receptor ACE2 for cellular entry. Nonetheless, the differences in ACE2 expression pattern in lung versus other normal and solid tumor tissues remain incompletely characterized. Here, we analyze a large data set comprising ACE2 mRNA expression for 7592 tissue samples across 22 types of primary solid tumor and 4461 samples across matched 18 non-diseased tissues. Our results unravel eight normal tissues and 10 primary solid tumors, which might be at high risk of SARS-CoV-2 infection. These findings may provide additional insight into the prevention and treatment of SARS-CoV-2 infection, in particular for patients with these 10 vulnerable cancer types.

manifestations have been observed (N. Guan et al., 2020; Huang et al., 2020; Wang et al., 2020) ,such as diarrhea, nausea or vomiting, liver abnormality, acute cardiac injury, and acute kidney injury. It is reported that cancer patients might harbor a higher risk of SARS-CoV-2 infection and inferior prognosis than those in infection without cancer (Liang et al., 2020) .

However, whether a heterogeneity of risk for infection exists among various cancer types remains unclear.

Here, we retrieved ACE2 mRNA expression data of 7592 tissue samples across 22 primary solid tumor types in The Cancer Genome Atlas (TCGA) and 4461 samples across 18 matched nondiseased tissues in Genotype-Tissue Expression (GTEx) from UCSC Xena (https://xena.ucsc.edu),

where TCGA and GTEx data were co-analyzed by the same Toil RNA-seq pipeline to eliminate computational batch effects (Vivian et al., 2017) . The expression values of ACE2 were quantified by RNA-Seq by Expectation-Maximization algorithm (Li and Dewey, 2011) and then normalized using the upper quartile method. The normalized values were log 2 -transformed after adding an offset of 1 to avoid taking log of zero before analysis. Further, virus abundances for TCGA tumors, quantified by numbers of virus-supporting reads per hundred million reads processed (RPHM), were obtained (Cao et al., 2016) . We defined a tumor sample with RPHM  100 for a given virus as virus-positive and examined ACE2 expression across seven tumor types with frequent viral presence according to the previous study (Cao et al., 2016) (Supplementary Table) . Mann-Whitney U test was used to compare the expression level between two groups. This study is exempt from ethical review because its data are publicly available and deidentified.

As shown in Fig. 1 , we observed a widespread distribution of ACE2 in these normal tissues, which is consistent with previous reports (Harmer et al., 2002) . Notably, eight normal tissues, including testis, kidney, thyroid, pancreas, breast, esophagus, liver, and ovary, had significantly higher ACE2 levels than lung as a reference (all P < 0.05), and expression levels in colon and bladder were similar to that in lung (both P > 0.05). The differences in the expression abundance between these tissues and lung indicate possible SARS-CoV-2 infection in extrapulmonary organs. For instance, J o u r n a l P r e -p r o o f Journal Pre-proof highest ACE2 abundance in testis may implicate its great possibility of SARS-CoV-2 exposure. Of note, a pathological analysis of testes from six patients who died of SARS showed that orchitis is a complication of SARS (Xu et al., 2006) . Thus, we propose strengthening follow-ups for reproductive functions of recovered SARS-CoV-2 male patients. Additionally, comparably high expression observed in kidney, liver, and colon may partially contribute to acute kidney injury, liver impairment, and diarrhea at onset of COVID-19, respectively (N. Guan et al., 2020; Huang et al., 2020; Wang et al., 2020) . Interestingly, our finding that ACE2 was highly expressed in breast appears to be in contrast to a retrospective study on nine pregnant women with COVID-19 in the third trimester, in which the colostrum from six patients tested negative for SARS-CoV-2 (H. . However, considering the small sample size and short duration of the study period, the risk of vertical transmission via breastfeeding deserves further investigations. Cancer cases analyzed in the study by Liang et al. (2020) is, however, limited in number (five lung cancers, four colorectal cancers, three breast cancers, two bladder cancers, and four other types of cancer). This restricts its ability to draw conclusions about the risk of SARS-CoV-2 infection in subgroups of specific cancer type. We hypothesize that another contributing mechanism may be increased likelihood of SARS-CoV-2 entry into certain cancer tissues due to aberrantly abundant ACE2 expression. Therefore, during this pandemic, we propose reinforcing personal protection, J o u r n a l P r e -p r o o f Journal Pre-proof such as remote medical counselling, minimizing the number of hospital visits, and appropriate isolation procedures when admitted to hospitals, for cancer patients, especially patient subgroups with these 10 solid tumor types. Interestingly, our analysis for seven types of tumors related to viruses revealed that virus-positive HNSC samples showed significantly lower ACE2 abundance than virus-negative ones (P < 0.01) (Fig. 3) , suggesting potential viral roles of ACE2 expression in HNSC.

In conclusion, we performed the first, to our knowledge, large-scale comparative analysis of ACE2 expression across multiple solid tumors and matched non-diseased tissues based on a consistently analyzed expression repository, which highlights eight normal tissues and 10 primary solid tumors with potentially similar or greater risk of SARS-CoV-2 exposure compared with lung, and identify a potential association between HNSC-related viruses and ACE2 expression. Given the large sample sizes for these 18 risky candidates, our results may be statistically robust and reliable;

notably, we did not use normal tissue samples from TCGA (adjacent to the tumor), because they are typically limited in number and their proximity to tumor may introduce signals of tumor microenvironment in their ACE2 expression profile (Aran et al., 2017) ; moreover, we used expression data unified by a standardized bioinformatic pipeline (Toil RNA-seq), enabling the direct comparison of ACE2 expression level from two sources (TCGA and GTEx). Our findings may contribute additional insight into the prevention and treatment of COVID-19, especially in patient subgroups with certain vulnerable cancer types. However, further clinical and autopsy studies are required to validate these findings. Box plots display the median and interquartile range, whiskers extend to 1.5 times the interquartile range, and outlier data are shown as dots. **, P < 0.01; ns, not significant. CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; COAD, colon adenocarcinoma; ESCA, esophageal carcinoma; HNSC, head and neck squamous cell carcinoma; LIHC, liver hepatocellular carcinoma; READ, rectum adenocarcinoma; STAD, stomach adenocarcinoma. 

Comprehensive analysis of normal adjacent to tumor transcriptomes

Divergent viral presentation among human tumors and adjacent normal tissues

Clinical characteristics and intrauterine vertical transmission potential of COVID-19 infection in nine pregnant women: a retrospective review of medical records

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

Clinical Characteristics of Coronavirus Disease 2019 in China

Quantitative mRNA expression profiling of

Toil enables reproducible, open source, big biomedical data analyses

Clinical Characteristics of 138 Hospitalized Patients With 2019 Novel Coronavirus-Infected Pneumonia in Wuhan

Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72 314 Cases From the Chinese Center for Disease Control and Prevention

Orchitis: a complication of severe acute respiratory syndrome (SARS)

A pneumonia outbreak associated with a new coronavirus of probable bat origin