Next Article in Journal
Blood Pressure and Tooth Loss: A Large Cross-Sectional Study with Age Mediation Analysis
Previous Article in Journal
Toxocariasis in a Child with Autism Spectrum Disorder
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 18
Issue 1
10.3390/ijerph18010284
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The scRNA-seq Expression Profiling of the Receptor ACE2 and the Cellular Protease TMPRSS2 Reveals Human Organs Susceptible to SARS-CoV-2 Infection

by
Jing Qi
†,
Yang Zhou
†,
Jiao Hua
†,
Liying Zhang
,
Jialin Bian
,
Beibei Liu
,
Zicen Zhao
and
Shuilin Jin
*
School of Mathematics, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to the work.
Int. J. Environ. Res. Public Health 2021, 18(1), 284; https://doi.org/10.3390/ijerph18010284
Submission received: 29 October 2020 / Revised: 21 December 2020 / Accepted: 22 December 2020 / Published: 2 January 2021
Browse Figures
Versions Notes

Abstract

:
COVID-19 patients always develop multiple organ dysfunction syndromes other than lungs, suggesting the novel virus SARS-CoV-2 also invades other organs. Therefore, studying the viral susceptibility of other organs is important for a deeper understanding of viral pathogenesis. Angiotensin-converting enzyme II (ACE2) is the receptor protein of SARS-CoV-2, and TMPRSS2 promotes virus proliferation and transmission. We investigated the ACE2 and TMPRSS2 expression levels of cell types from 31 organs to evaluate the risk of viral infection using single-cell RNA sequencing (scRNA-seq) data. For the first time, we found that the gall bladder and fallopian tube are vulnerable to SARS-CoV-2 infection. Besides, the nose, heart, small intestine, large intestine, esophagus, brain, testis, and kidney are also identified to be high-risk organs with high expression levels of ACE2 and TMPRSS2. Moreover, the susceptible organs are grouped into three risk levels based on the ACE2 and TMPRSS2 expression. As a result, the respiratory system, digestive system, and urinary system are at the top-risk level for SARS-CoV-2 infection. This study provides evidence for SARS-CoV-2 infection in the human nervous system, digestive system, reproductive system, respiratory system, circulatory system, and urinary system using scRNA-seq data, which helps in the clinical diagnosis and treatment of patients.

1. Introduction

In December 2019, a novel coronavirus pneumonia recently named Coronavirus Disease 2019 (COVID-19) by the World Health Organization (WHO) [1], was first reported in several patients in Wuhan, China [2,3,4,5]. As of December 2020, COVID-19 has spread widely in more than 200 countries; more than 71.43 million people have suffered from the disease, and over 1.6 million people have died, posing a major threat to global public health.
COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [1,6,7], which seriously damages the respiratory system. Patients develop acute respiratory infection symptoms, and even acute respiratory distress syndrome (ARDS), acute respiratory failure, and other severe symptoms [4,5,8,9]. On the other side, complications that occurred outside the lungs, including acute kidney injury, liver function damage, and multiple organ failure, were commonly described in COVID-19 cases and even led to death, suggesting that the virus invades other organs at the same time [4,5,10,11,12,13]. SARS-CoV-2 enters the cell via the binding of spike (S) protein and angiotensin-converting enzyme II (ACE2) [14,15,16,17], the receptor protein of SARS-CoV, and NL63 [18,19,20,21]. The anti-ACE2 antibody blocks the viral entry driven by SARS-CoV-2 S protein [17]. In addition, human adenocarcinoma cells (A549) were found to be incompatible with SARS-CoV-2 infection because of the low expression of ACE2 [22,23]. The cellular protease TMPRSS2 primes the S protein and promotes the transmission during the viral infection, and the inhibitor of TMPRSS2 blocks the SARS-CoV-2 infection of lung cells [17,24]. Thus, the expression of ACE2 and TMPRSS2 is indispensable for viral infection in cells. The distribution and expression of ACE2 and TMPRSS2 are strongly associated with the target organ of the SARS-CoV-2 infection [25,26,27,28]. Previous studies have shown that smoking and aging induce an increase in ACE2 expression in the human respiratory tissue, while smokers and people with old age are vulnerable to COVID-19 [25,29,30,31]. It has been shown that the androgens can up-regulate the expression of ACE2 and TMPRSS2, which is related to the higher mortality and morbidity of males relative to females in the COVID-19 pandemic [31,32]. Moreover, several studies suggested that the blocking of androgen signaling and the down-regulation of ACE2 and TMPRSS2 are protective against COVID-19 [26,33]. Targeting the transcriptional regulation of TMPRSS2 and ACE2 has been an important strategy to prevent SARS-CoV-2 infection in clinical treatment [32,34].
Single-cell RNA sequencing (scRNA-seq) can obtain the gene expression profile of a single cell, which better reveals the heterogeneity of cells, and allows better understanding of the functions of an individual cell in its microenvironment [35,36,37,38,39]. Therefore, scRNA-seq technology provides a tool to study the pathogenic mechanism of SARS-CoV-2 on human cells from the expression of genes in cell resolution. The pathological inference and gross anatomical observations indicated that the lesions caused by the novel coronavirus were mainly in the lung [9,40,41]. ACE2 is mainly expressed in type II alveolar cells (AT2) in the lung [42,43,44], which implies that the AT2 cell is vulnerable to SARS-CoV-2 infection. Researchers have obtained the susceptibility of other organs from different systems using ACE2 expression in AT2 cells as a baseline. Following this way, many respiratory districts were considered, and ACE2 was reported to be highly expressed in the nasal tissue, mouth, airway, and lung [43,45,46]. The esophagus, large intestine (ileum and colon), and pancreas were identified as high-risk organs in the digestive system [43,46,47,48]. The kidney and bladder, as major organs of the urinary system, were also indicated to be high ACE2-expressed [43,49,50,51]. Besides, the testes and uterus were manifested to be susceptible organs, implying that the reproductive system was a potential route of viral infection [52,53]. Several studies have utilized the expression of ACE2 and TMPRSS2 to predict organ susceptibility. Zhou et al. analyzed human post-mortem eyes and surgical specimens for the expression of ACE2 and TMPRSS2, and results show that the ocular surface cells including conjunctiva are at high risk for SARS-CoV-2 infection [54]. Seow et al. reported the co-expression of ACE2 and TMPRSS2 in a TROP2+ liver progenitor population and identified a potentially high-risk liver cell-type for viral ingress [55]. Lukassen et al. investigated the expression levels and distributions of ACE2 and TMPRSS2 across cell types in lung tissue and bronchial branches, respectively [56]. Furthermore, ACE2 and TMPRSS2 are highly expressed in a transient secretory cell type of bronchial branches, suggesting the increased vulnerability for SARS-CoV-2 infection in this cell type [56].
In this paper, the expression level of ACE2 and TMPRSS2, in different cell types of organs from nine major systems (including the respiratory system, digestive system, nervous system, endocrine system, reproductive system, circulatory system, urinary system, and motor system), were obtained using the scRNA-seq data. For the first time, we found that the fallopian tube and gall bladder are vulnerable to SARS-CoV-2 infection. Besides, the nose (nasal brushing epithelial cells, nasal turbinate epithelial cells, and nasal airway epithelial cells), heart, small intestine (jejunum, ileum, and duodenum), large intestine (rectum and colon), esophagus, brain (substantia nigra and cortex), testes, and kidney are predicted as high-risk organs under a more rigorous standard. Moreover, as the spike (S) protein initiated by TMPRSS2 is essential for the entry of the virus into the target cells and the transmission of the virus in the infected host, we combined the expression level of ACE2 with the expression level of TMPRSS2 to predict the risk level of each susceptible organ, and found that the respiratory system, digestive system, and urinary system are at the highest level of vulnerability to SARS-CoV-2 infection.

2. Materials and Methods

The available scRNA-seq data of healthy humans were collected for the analysis, including 31 organs from nine major human systems (Table 1). The data were downloaded from the Gene Expression Omnibus (GEO) database and the Tissue Stability Cell Atlas, and the details of the data resources are in Supplementary File Table S1.
In our paper, analysis of scRNA-seq data was conducted in R environment (version 3.6.1; The R Foundation for Statistical Computing, Topeka, KS, USA) using Seurat package (https://github.com/satijalab/seurat). First, according to the number of expressed genes, counts and cells with mitochondrial content, the low-quality cells were filtered, and the cells were reserved within the range of μ − σ and μ + σ (where μ is the mean and σ is the standard deviation of the numbers). Then, the data were normalized by the logarithmic transformation and the downstream analysis was carried out with the top 2000 most variable genes. The ScaleData function were used for linear regression and principal component analysis (PCA) was performed on the scaled data for linear dimensional reduction. Finally, we performed the cell cluster analysis, gene differential expression analysis by t-test, and annotated the cell subtypes.
Here, the gene-expressed cell means that the gene expression value UMI (unique molecular identifier) is greater than or equal to 1, or the TPM (transcripts per million) is greater than 0, or the FPKM (fragments per kilobase million) is greater than 0. The ratio of ACE2-expressed cells and the ratio of TMPRSS2-expressed cells of the cell cluster are denoted as R A C E 2 and R T M P R S S 2 , respectively. The average R A C E 2 over total AT2 cells across 8 samples is 0.79%. If R A C E 2 ≥ 0.79 % and R T M P R S S 2 > 0 , the cells in the cluster are identified as highly susceptible cells, and the corresponding organs are inferred as vulnerable to COVID-19. Moreover, we used the geometric mean of the ratios of TMPRSS2-expressed cells and the ratios of ACE2-expressed cells as the risk parameter to predict the risk for SARS-CoV-2 infection, and the risk parameter is denoted as:
R = R A C E 2 · R T M P R S S 2
According to the risk parameters, the susceptible organs were sorted into three groups. Specifically, the cells with a risk parameter greater than 0.1 are defined as level 1, which is the highest level of risk; the cells with a risk parameter of more than 0.05 and less than 0.1 are level 2, which denotes a higher risk; the rest are level 3, meaning there is an existing risk of infection.

3. Results

3.1. Respiratory System

The scRNA-seq data of the lung, nose, trachea, and bronchus in the respiratory system were collected for analysis. In the lung, AT2 cells contain an average of 0.79% ACE2-expressed cells and 21.20% TMPRSS2-expressed cells across eight samples (Figure 1A,B, Supplementary File: Figures S1–S8, and Table 2), and the expression levels of ACE2 and TMPRSS2 are high in AT2 cells (Figure 1C,D). The data of nose (nasal brushing epithelial cells, nasal turbinate epithelial cells, and nasal airway epithelial cells) contain ACE2-expressed and TMPRSS2-expressed cell clusters (Supplementary File: Figure S9–S11), and the ratios of ACE2-expressed cells of these cell clusters are all above 0.79% (Table 2); thus, the nose is identified as the high-risk organ. The low ratio of ACE2-expressed cells in the bronchus and trachea means that they are low-risk infection organs (Supplementary File: Figures S12–S13).

3.2. Digestive System

The scRNA-seq data of the jejunum, ileum, duodenum, rectum, colon, esophagus, gall bladder, pancreatic islets, liver, and stomach from the digestive system were collected for the analysis. The primordium cells from the gall bladder contain 2.6% TMPRSS2-expressed cells and 2.2% ACE2-expressed cells (Figure 2, Table 2), which means the gall bladder is vulnerable to the SARS-CoV-2 infection. Moreover, the small intestine (jejunum, ileum, and duodenum), large intestine (rectum and colon), and the esophagus are identified as high-risk organs (Supplementary File: Figures S14–S19, Table 2). However, no cell clusters from the liver, stomach, and pancreatic islets data show high ACE2 and TMPRSS2 expression levels (Supplementary File: Figures S20–S22), which demonstrates a low infection risk.

3.3. Reproductive System

The scRNA-seq data of the testis, fallopian tube, ovary, and uterus from the reproductive system were collected for analysis. The ratios of the TMPRSS2-expressed cells and the ACE2-expressed cells in the epithelial cells of the fallopian tube are 26.5% and 1.4%, respectively (Figure 3, Table 2), and the ovarian somatic cells contain 1% TMPRSS2-expressed cells and 1% ACE2-expressed cells, so the fallopian tube is identified as a high-risk organ. The testis is also identified as a high-risk organ because of the high expression level of TMPRSS2 and ACE2 (Supplementary File: Figure S23, Table 2). Low ratios of ACE2-expressed cells in the ovary and uterus mean the ovary and uterus are low infection risk organs (Supplementary File: Figures S24–S25).

3.4. Nervous System

The scRNA-seq of the substantia nigra and cortex, hippocampus, cerebellum, spinal cord, and neuronal epithelium from the nervous system were collected to infer the susceptibility of the organs. The analysis results show that ACE2 is expressed in the oligodendrocyte precursor cells and astrocytes of the substantia nigra and cortex with a high level, and TMPRSS2 is expressed as well. More specifically, astrocytes contain 1.9% ACE2-expressed cells, and oligodendrocyte precursor cells contain 1.6% ACE2-expressed cells (Supplementary File: Figure S26, Table 2). Therefore, the substantia nigra and cortex are predicted as high-risk tissues, and the brain is identified as a high-risk organ. For the analysis of other districts, cells from the hippocampus have low expression levels of TMPRSS2 and ACE2, and the cerebellum, spinal cord, and neuronal epithelium data show zero expression of TMPRSS2 and ACE2 (Supplementary File: Figures S27–S30), which demonstrates a low infection risk of these districts.

3.5. Circulatory System

The data of the heart, spleen, and artery from the circulatory system were collected to infer the susceptibility of these organs. The cardiomyocytes and cardiovascular progenitor cells from the heart contain 6.6% and 12.5% ACE2-expressed cells, respectively, and the TMPRSS2 is expressed in both cell clusters as well. Consequently, the heart is considered a high-risk organ (Supplementary File: Figure S31, Table 2). Nevertheless, almost no cells of the spleen, artery, and peripheral blood data show high TMPRSS2 and ACE2 expression levels, which means they are low-risk organs (Supplementary File: Figures S32–S34).

3.6. Urinary System

The scRNA-seq data of the kidney, ureter, and prostate from the urinary system were utilized for the data analysis. The analysis results of the kidney scRNA-seq data show high ACE2 and TMPRSS2 expression levels in the nephron epithelial cells, epithelial cells, endothelial cells, and mesangial cells. Particularly, the ratios of TMPRSS2-expressed are 10.7%, 9.6%, 12.8%, and 14.5%, respectively, and the ratios of ACE2-expressed are 2.7%, 2.7%, 2.7%, and 3.0%, respectively (Supplementary File: Figure S35, Table 2). Therefore, the kidney is considered a high-risk organ. In addition, the ACE2 is not expressed in the ureter and prostate cells (Supplementary File: Figures S36 and S37), and they are predicted to be low-risk infection organs.

3.7. Endocrine System, Motor System and Immune System

The susceptibility of organs from the endocrine system, immune system, and motor system was also considered. In the endocrine system, almost no cells from the thyroid gland data show high ACE2 and TMPRSS2 expression levels, and the thymus gland data contains no ACE2-expressed cells (Supplementary File: Figures S38 and S39). Hence, they are not high-risk organs. Likewise, few ACE2-expressed cells are found in the muscle from the motor system (Supplementary File: Figure S40). There is no ACE2 expression in the lymph nodes, tonsil (tonsil dendritic cells), and bone marrow data from the immune system (Supplementary File: Figures S41–S43), which means they are low-risk organs.

3.8. The Risk Levels of Susceptible Organs

Based on TMPRSS2 and ACE2 expression levels, we grouped the susceptible organs into three risk levels. According to the clinical implication, the lung should be the highest risk organ to SARS-CoV-2 infection. Interestingly, the large intestine (colon and rectum), esophagus, and nose (nasal airway epithelium) are the most susceptible organs, and the result indicates SARS-CoV-2 mainly attacks the respiratory system and the digestive system (Table 3, Figure 4). The kidney, small intestine (duodenum and jejunum), and fallopian tube are susceptible organs with moderate risk (Table 3, Figure 4). In addition, the testis, gall bladder, brain (substantia nigra and cortex), and heart are identified to be potentially susceptible organs (Table 3, Figure 4).

4. Discussion

ACE2 has been reported to show a significant overexpression in COVID-19 patients and is positively correlated with the expression of some other SARS-COV-2 host invasion genes [27]. Compared with current research on organ susceptibility to COVID-19 using scRNA-seq data, the advantages of our research are mainly reflected in the following aspects. Firstly, we presented a more comprehensive and rigorous analysis of 31 human organs from nine systems based on ACE2 and TMPRSS2 expression. Then, we obtained breakthrough results that the gall bladder and fallopian tube are vulnerable to SARS-CoV-2 infection. Finally, we classified susceptible organs to different risk levels and found that the respiratory system, digestive system, and urinary system are at the top-risk level for SARS-CoV-2 infection.
At present, nucleic acid testing is the most extensive testing technology for mild patients and suspected cases. The screening samples for nucleic acid detection are mostly from deep cough sputum, oropharyngeal swabs, or nasopharyngeal swabs. However, this detection technology may cause false-negative results. In some COVID-19 cases, patients have gastrointestinal (GI) symptoms, such as anorexia, vomiting and diarrhea, without any respiratory symptoms in the initial stages of the disease [70,71]. Our results verify that the esophagus and intestines (rectum, ileum, jejunum, and duodenum) are also the main sites of SARS-CoV-2 infection in the human body. Previous studies have reported that SARS-CoV-2 can be detected in the feces of infected persons, suggesting that there may be food transmission and fecal–oral transmission during the spread of SARS-CoV-2 [72,73]. Therefore, epidemic prevention experts recommend paying attention to good personal hygiene and food handling practices. In addition, neurological symptoms such as headache, loss of sense of smell, amblyopia, weakness, and myalgia were described in some COVID-19 patients [5,74,75]. Our results show that the expression ratio of ACE2 in the cortex and substantia nigra of healthy people is higher than that in AT2 cells, which indicates that SARS-CoV-2 can also invade the nervous system. Therefore, the symptoms of some organs of the digestive system and nervous system can be used as the basis for preliminary judgment of the patient’s condition.
Previous studies have reported PCR fragments of the coronavirus detected in the blood and urine of COVID-19 patients, as well as nephritis histological analysis of the autopsy results, indicating that SARS-CoV-2 causes acute tubular disease [76,77]. The renal tropism of SARS-CoV-2 revealed by the high expression of ACE2 and TMPRSS2 in the kidney may explain the renal injury in COVID-19 patients. The results show that the heart is also a susceptible organ, indicating that SARS-CoV-2 may affect the normal operation of the cardiovascular system. In fact, SARS-CoV-2 infection has been associated with cardiovascular complications, such as acute myocardial injury, and myocarditis, which increases mortality [76,78]. Besides, patients with cardiovascular and cerebrovascular diseases are prone to COVID-19-related complications due to their relatively fragile cardiovascular system [76]. Therefore, it is necessary to pay attention to the pathological changes of the urinary system and circulatory system organs when treating COVID-19 patients. Moreover, it should be noted that, according to the results of this study, SARS-CoV-2 may affect the reproductive organs. We suggest that patients should contraception during antiviral treatment, and it is recommended to continue contraception after a period of treatment.
However, the cell entry mechanisms of SARS-CoV-2 are not fully understood, and the invasion process may also be highly related to other genes than ACE2 and TMPRSS2. Therefore, our results need to be further confirmed by clinical observations and biological experiments. Due to the limitation of data collection, the susceptibility of some human organs remains to be analyzed. Moreover, whether some chronic diseases or cancer patients are more susceptible to SARS-CoV-2 infection is also worth exploring.

5. Conclusions

The clinical symptoms of patients infected by SARS-CoV-2 are mainly manifested in the respiratory system and digestive system, including cough, shortness of breath, dyspnea, and diarrhea. However, some patients also developed symptoms such as heart damage and kidney failure, indicating that the virus affected the normal function of the circulatory and urinary systems. We investigated the susceptibility of the organs and tissues in various human systems based on the scRNA-seq data analysis. In detail, 31 organs from nine major human systems were considered, out of which 11 organs were identified to be susceptible to the virus. Moreover, we classified these susceptible organs into three levels in terms of their risk, which provide novel ideas for the follow-up detection of the virus, treatment, and the monitoring of recrudescence.
Through the assessment of the susceptibility of human organs to SARS-CoV-2 infection and the accurate judgment of the viral invasion to organs, we can further understand the pathophysiology of the disease and help determine the prognostic utility of clinical and laboratory parameters related to COVID-19, which will lead to the development of more precise and effective management methods for the disease.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/18/1/284/s1, Figure S1: The lung scRNA-seq data analysis results (donor 1). (A) UMAP visualization of clustering results for the lung cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S2: The lung scRNA-seq data analysis results (donor 2). (A) UMAP visualization of clustering results for the lung cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S3: The lung scRNA-seq data analysis results (donor 3). (A) UMAP visualization of clustering results for the lung cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S4: The lung scRNA-seq data analysis results (donor 4). (A) UMAP visualization of clustering results for the lung cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S5: The lung scRNA-seq data analysis results (donor 5). (A) UMAP visualization of clustering results for the lung cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S6: The lung scRNA-seq data analysis results (donor 6). (A) UMAP visualization of clustering results for the lung cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S7: The lung scRNA-seq data analysis results (donor 7). (A) UMAP visualization of clustering results for the lung cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S8: The lung scRNA-seq data analysis results (donor 8). (A) UMAP visualization of clustering results for the lung cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S9: High ACE2 and TMPRSS2 expression levels of the mesenchymal stromal cells, plasma cells in the nasal turbinate epithelial cells. (A) UMAP visualization of clustering results for the nasal turbinate epithelial cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S10: High ACE2 and TMPRSS2 expression levels of the mesenchymal stromal cells, plasma cells in the nasal brushing epithelial cells. (A) UMAP visualization of clustering results for the nasal brushing epithelial cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S11: High ACE2 and TMPRSS2 expression levels in the nasal airway epithelial cells. (A) UMAP visualization of clustering results for the airway epithelial cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S12: The bronchus scRNA-seq data analysis results. (A) UMAP visualization of clustering results for the bronchus cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S13: The trachea scRNA-seq data analysis results. (A) UMAP visualization of clustering results for the trachea cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S14: High ACE2 and TMPRSS2 expression level of enterocyte progenitor cells and goblet cells in the jejunum. (A)UMAP visualization of clustering results for jejunum cells. (B)The ratio of ACE2-expressed cells in each cell cluster. (C)The ratio of TMPRSS2-expressed cells in each cell cluster. (D)ACE2 expression level in each cell cluster on the UMAP plot. (E)TMPRSS2 expression level in each cell cluster on the UMAP plot. (F)The expression distribution of ACE2 across each cell cluster. (G)The expression distribution of TMPRSS2 across each cell cluster, Figure S15: High ACE2 and TMPRSS2 expression level of the intestinal epithelial stem cells and enterocyte progenitor cells in the ileum. (A) UMAP visualization of clustering results for ileum cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S16: High ACE2 and TMPRSS2 expression level of the intestinal LGR5+ stem cells, epithelial stem cells, enterocyte progenitor cells, tuft progenitor cells, and enteroendocrine cells in the duodenum. (A) UMAP visualization of clustering results for duodenum cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S17: High ACE2 and TMPRSS2 expression level of goblet progenitor cells, MKI67+ progenitor cells, enterocytes, and goblet cells in the rectum. (A) UMAP visualization of clustering results for rectum cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S18: High ACE2 and TMPRSS2 expression level of the enterocytes and goblet cells in the colon. (A) UMAP visualization of clustering results for colon cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S19: High ACE2 and TMPRSS2 expression level of the secretory progenitor cells in the esophagus. (A) UMAP visualization of clustering results for esophagus cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S20: The liver scRNA-seq data analysis results. (A) UMAP visualization of clustering results for liver cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S21: The stomach scRNA-seq data analysis results. (A) UMAP visualization of clustering results for stomach cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S22: The pancreatic islets scRNA-seq data analysis results. (A) UMAP visualization of clustering results for pancreatic islets cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S23: High ACE2 and TMPRSS2 expression level of spermatogonium, peritubular myoid cells, testis somatic cells, and spermatogonial stem cells in the testis. (A) UMAP visualization of clustering results for testis cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S24: The ovary scRNA-seq data analysis results. (A) UMAP visualization of clustering results for ovary cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) ACE2 expression level in each cell cluster on the UMAP plot. (D) The expression distribution of ACE2 across each cell cluster, Figure S25: The uterus scRNA-seq data analysis results. (A) UMAP visualization of clustering results for uterus cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) ACE2 expression level in each cell cluster on the UMAP plot. (D) The expression distribution of ACE2 across each cell cluster, Figure S26: High ACE2 and TMPRSS2 expression level of oligodendrocyte precursor cells and astrocytes in the substantia nigra and cortex. (A) UMAP visualization of clustering results for the substantia nigra and cortex cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S27: The hippocampus scRNA-seq data analysis results. (A) UMAP visualization of clustering results for hippocampus cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S28: UMAP visualization of clustering results for the cerebellum cells, Figure S29: UMAP visualization of clustering results for the spinal cord cells, Figure S30: UMAP visualization of clustering results for the neuronal epithelium cells, Figure S31: High ACE2 and TMPRSS2 expression level of the cardiomyocytes and cardiovascular progenitor cells in the heart. (A) UMAP visualization of clustering results for the heart cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S32: The spleen scRNA-seq data analysis results. (A) UMAP visualization of clustering results for the spleen cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S33: The artery scRNA-seq data analysis results. (A) UMAP visualization of clustering results for the spleen cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) ACE2 expression level in each cell cluster on the UMAP plot. (D) The expression distribution of ACE2 across each cell cluster, Figure S34: UMAP visualization of clustering results for the peripheral blood cells, Figure S35: High ACE2 and TMPRSS2 expression level of the nephron epithelial cells, epithelial cells, endothelial cells, and mesangial cells in the kidney. (A) UMAP visualization of clustering results for the kidney cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S36: UMAP visualization of clustering results for the ureter cells, Figure S37: UMAP visualization of clustering results for the prostate cells, Figure S38: The thyroid gland scRNA-seq data analysis results. (A) UMAP visualization of clustering results for the thyroid gland cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) The ratio of TMPRSS2-expressed cells in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster, Figure S39: UMAP visualization of clustering results for the thymus gland cells, Figure S40: The muscle scRNA-seq data analysis results. (A) UMAP visualization of clustering results for the muscle cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) ACE2 expression level in each cell cluster on the UMAP plot. (D) The expression distribution of ACE2 across each cell cluster, Figure S41: UMAP visualization of clustering results for the lymph nodes cells, Figure S42: UMAP visualization of clustering results for the tonsil dendritic cells, Figure S43: The bone marrow scRNA-seq data analysis results. (A) UMAP visualization of clustering results for the bone marrow cells. (B) The ratio of ACE2-expressed cells in each cell cluster. (C) ACE2 expression level in each cell cluster on the UMAP plot. (D) The expression distribution of ACE2 across each cell cluster, Table S1: The data sources of organs and tissues.

Author Contributions

Conceptualization, J.Q., Y.Z., J.H. and S.J.; data collection, J.Q., Y.Z., J.H., L.Z., J.B. and B.L.; statistical analysis, J.H. and L.Z.; writing—review and editing, J.Q., Y.Z., J.H., S.J., L.Z., J.B., B.L. and Z.Z.; writing—original draft, J.Q., Y.Z. and J.H.; formal analysis, J.Q., Y.Z. and J.H.; funding acquisition, S.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (grant number 11971130).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data availability statement

The data were downloaded from the Tissue Stability Cell Atlas, GEO database, and the details of the data resource can be seen in Supplementary File Table S1.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. Naming the Coronavirus Disease (COVID-19) and the Virus that Causes It. Available online: https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/naming-the-coronavirus-disease-(covid-2019)-and-the-virus-that-causes-it (accessed on 25 December 2020).
  2. Li, Q.; Guan, X.; Wu, P.; Wang, X.; Zhou, L.; Tong, Y.; Ren, R.; Leung, K.S.M.; Lau, E.H.Y.; Wong, J.Y.; et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N. Engl. J. Med. 2020, 382, 1199–1207. [Google Scholar] [CrossRef] [PubMed]
  3. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef] [PubMed]
  4. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X.; et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 2020, 395, 497–506. [Google Scholar] [CrossRef] [Green Version]
  5. Chen, N.; Zhou, M.; Dong, X.; Qu, J.; Gong, F.; Han, Y.; Qiu, Y.; Wang, J.; Liu, Y.; Wei, Y.; et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study. Lancet 2020, 395, 507–513. [Google Scholar] [CrossRef] [Green Version]
  6. Coronaviridae Study Group of the International Committee on Taxonomy of Viruses. The species Severe acute respiratory syndrome-related coronavirus: Classifying 2019-nCoV and naming it SARS-CoV-2. Nat. Microbiol. 2020, 5, 536–544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Lu, R.; Zhao, X.; Li, J.; Niu, P.; Yang, B.; Wu, H.; Wang, W.; Song, H.; Huang, B.; Zhu, N.; et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding. Lancet 2020, 395, 565–574. [Google Scholar] [CrossRef] [Green Version]
  8. Wölfel, R.; Corman, V.M.; Guggemos, W.; Seilmaier, M.; Zange, S.; Müller, M.A.; Niemeyer, D.; Jones, T.C.; Vollmar, P.; Rothe, C.; et al. Virological assessment of hospitalized patients with COVID-2019. Nature 2020, 581, 465–469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Xu, Z.; Shi, L.; Wang, Y.; Zhang, J.; Huang, L.; Zhang, C.; Liu, S.; Zhao, P.; Liu, H.; Zhu, L.; et al. Pathological findings of COVID-19 associated with acute respiratory distress syndrome. Lancet Respir. Med. 2020, 8, 420–422. [Google Scholar] [CrossRef]
  10. Braun, F.; Lütgehetmann, M.; Pfefferle, S.; Wong, M.N.; Carsten, A.; Lindenmeyer, M.T.; Nörz, D.; Heinrich, F.; Meißner, K.; Wichmann, D.; et al. SARS-CoV-2 renal tropism associates with acute kidney injury. Lancet 2020, 396, 597–598. [Google Scholar] [CrossRef]
  11. Fanelli, V.; Fiorentino, M.; Cantaluppi, V.; Gesualdo, L.; Stallone, G.; Ronco, C.; Castellano, G. Acute kidney injury in SARS-CoV-2 infected patients. Crit. Care. 2020, 24, 155. [Google Scholar] [CrossRef] [Green Version]
  12. Cheng, Y.; Luo, R.; Wang, K.; Zhang, M.; Wang, Z.; Dong, L.; Li, J.; Yao, Y.; Ge, S.; Xu, G. Kidney disease is associated with in-hospital death of patients with COVID-19. Kidney Int. 2020, 97, 829–838. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, C.; Shi, L.; Wang, F.-S. Liver injury in COVID-19: Management and challenges. Lancet Gastroenterol. Hepatol. 2020, 5, 428–430. [Google Scholar] [CrossRef]
  14. Wrapp, D.; Wang, N.; Corbett, K.; Goldsmith, J.A.; Hsieh, C.-L.; Abiona, O.; Graham, B.S.; McLellan, J.S. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 2020, 367, 1260–1263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ou, X.; Liu, Y.; Lei, X.; Li, P.; Mi, D.; Ren, L.; Guo, L.; Guo, R.; Chen, T.; Hu, J.; et al. Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV. Nat. Commun. 2020, 11, 1620. [Google Scholar] [CrossRef] [Green Version]
  16. Wan, Y.; Shang, J.; Graham, R.; Baric, R.S.; Li, F. Receptor Recognition by the Novel Coronavirus from Wuhan: An Analysis Based on Decade-Long Structural Studies of SARS Coronavirus. J. Virol. 2020, 94, e00127-20. [Google Scholar] [CrossRef] [Green Version]
  17. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.-H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 271–280.e8. [Google Scholar] [CrossRef]
  18. Gallagher, T.M.; Buchmeier, M.J. Coronavirus spike proteins in viral entry and pathogenesis. Virology 2001, 279, 371–374. [Google Scholar] [CrossRef] [Green Version]
  19. He, L.; Ding, Y.; Zhang, Q.; Che, X.; He, Y.; Shen, H.; Wang, H.; Li, Z.; Zhao, L.; Geng, J.; et al. Expression of elevated levels of pro-inflammatory cytokines in SARS-CoV-infected ACE2+ cells in SARS patients: Relation to the acute lung injury and pathogenesis of SARS. J. Pathol. 2006, 210, 288–297. [Google Scholar] [CrossRef]
  20. Li, W.; Sui, J.; Huang, I.-C.; Kuhn, J.H.; Radoshitzky, S.R.; Marasco, W.A.; Choe, H.; Farzan, M. The S proteins of human coronavirus NL63 and severe acute respiratory syndrome coronavirus bind overlapping regions of ACE2. Virology 2007, 367, 367–374. [Google Scholar] [CrossRef] [Green Version]
  21. Wu, K.; Li, W.; Peng, G.; Li, F. Crystal structure of NL63 respiratory coronavirus receptor-binding domain complexed with its human receptor. Proc. Natl. Acad. Sci. USA 2009, 106, 19970–19974. [Google Scholar] [CrossRef] [Green Version]
  22. Blanco-Melo, D.; Nilsson-Payant, B.E.; Liu, W.C.; Uhl, S.; Hoagland, D.; Møller, R.; Oishi, K.; Panis, M.; Sachs, D.; Wang, T.T.; et al. Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19. Cell 2020, 181, 1036–1045.e9. [Google Scholar] [CrossRef] [PubMed]
  23. Harcourt, J.; Tamin, A.; Lu, X.; Kamili, S.; Sakthivel, S.K.; Murray, J.; Queen, K.; Tao, Y.; Paden, C.R.; Zhang, J.; et al. Severe Acute Respiratory Syndrome Coronavirus 2 from Patient with Coronavirus Disease, United States. Emerg. Infect. Dis. 2020, 26, 1266–1273. [Google Scholar] [CrossRef] [PubMed]
  24. Meng, T.; Cao, H.; Zhang, H.; Kang, Z.; Xu, D.; Gong, H.; Wang, J.; Li, Z.; Cui, X.; Xu, H.; et al. The insert sequence in SARS-CoV-2 enhances spike protein cleavage by TMPRSS. bioRxiv 2020. [Google Scholar] [CrossRef]
  25. Matsuyama, S.; Nao, N.; Shirato, K.; Kawase, M.; Saito, S.; Takayama, I.; Nagata, N.; Sekizuka, T.; Katoh, H.; Kato, F.; et al. Enhanced isolation of SARS-CoV-2 by TMPRSS2-expressing cells. Proc. Natl. Acad. Sci. USA 2020, 117, 7001–7003. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Stopsack, K.H.; Mucci, L.A.; Antonarakis, E.S.; Nelson, P.S.; Kantoff, P.W. TMPRSS2 and COVID-19: Serendipity or Opportunity for Intervention? Cancer Discov. 2020, 10, 779–782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Amati, F.; Vancheri, C.; Latini, A.; Colona, V.L.; Grelli, S.; D’Apice, M.R.; Balestrieri, E.; Passarelli, C.; Minutolo, A.; Loddo, S.; et al. Expression profiles of the SARS-CoV-2 host invasion genes in nasopharyngeal and oropharyngeal swabs of COVID-19 patients. Heliyon 2020, 6, e05143. [Google Scholar] [CrossRef]
  28. Li, Y.; Zhou, W.; Yang, L.; You, R. Physiological and pathological regulation of ACE2, the SARS-CoV-2 receptor. Pharmacol. Res. 2020, 157, 104833. [Google Scholar] [CrossRef]
  29. Russo, P.; Bonassi, S.; Giacconi, R.; Malavolta, M.; Tomino, C.; Maggi, F. COVID-19 and smoking: Is nicotine the hidden link? Eur. Respir. J. 2020, 55, 2001116. [Google Scholar] [CrossRef]
  30. Leung, J.M.; Yang, C.X.; Tam, A.; Shaipanich, T.; Hackett, T.L.; Singhera, G.K.; Dorscheid, D.R.; Sin, D.D. ACE-2 expression in the small airway epithelia of smokers and COPD patients: Implications for COVID-19. Eur. Respir. J. 2020, 55, 2000688. [Google Scholar] [CrossRef] [Green Version]
  31. Chakladar, J.; Shende, N.; Li, W.T.; Rajasekaran, M.; Chang, E.Y.; Ongkeko, W.M. Smoking-Mediated Upregulation of the Androgen Pathway Leads to Increased SARS-CoV-2 Susceptibility. Int. J. Mol. Sci. 2020, 21, 3627. [Google Scholar] [CrossRef]
  32. Qiao, Y.; Wang, X.-M.; Mannan, R.; Pitchiaya, S.; Zhang, Y.; Wotring, J.W.; Xiao, L.; Robinson, D.R.; Wu, Y.-M.; Tien, J.C.-Y.; et al. Targeting transcriptional regulation of SARS-CoV-2 entry factors ACE2 and TMPRSS2. Proc. Natl. Acad. Sci. USA 2020. [Google Scholar] [CrossRef]
  33. Goren, A.; Wambier, C.G.; Herrera, S.; McCoy, J.; Vañó-Galván, S.; Gioia, F.; Comeche, B.; Ron, R.; Serrano-Villar, S.; Ramos, P.M.; et al. Anti-androgens may protect against severe COVID-19 outcomes: Results from a prospective cohort study of 77 hospitalized men. J. Eur. Acad. Dermatol. Venereol. 2020. [Google Scholar] [CrossRef] [PubMed]
  34. Ragia, G.; Manolopoulos, V.G. Inhibition of SARS-CoV-2 entry through the ACE2/TMPRSS2 pathway: A promising approach for uncovering early COVID-19 drug therapies. Eur. J. Clin. Pharmacol. 2020, 76, 1623–1630. [Google Scholar] [CrossRef] [PubMed]
  35. McDavid, A.; Finak, G.; Chattopadyay, P.K.; Dominguez, M.; Lamoreaux, L.; Ma, S.S.; Roederer, M.; Gottardo, R. Data exploration, quality control and testing in single-cell qPCR-based gene expression experiments. Bioinformatics 2013, 29, 461–467. [Google Scholar] [CrossRef] [PubMed]
  36. Eberwine, J.; Sul, J.Y.; Bartfai, T.; Kim, J. The promise of single-cell sequencing. Nat. Methods 2014, 11, 25–27. [Google Scholar] [CrossRef] [PubMed]
  37. Wagner, A.; Regev, A.; Yosef, N. Revealing the vectors of cellular identity with single-cell genomics. Nat. Biotechnol. 2016, 34, 1145–1160. [Google Scholar] [CrossRef] [Green Version]
  38. Rizzetto, S.; Eltahla, A.A.; Lin, P.; Bull, R.; Lloyd, A.R.; Ho, J.W.K.; Venturi, V.; Luciani, F. Impact of sequencing depth and read length on single cell RNA sequencing data of T cells. Sci. Rep. 2017, 7, 12781. [Google Scholar] [CrossRef]
  39. Kalisky, T.; Oriel, S.; Bar-Lev, T.H.; Ben-Haim, N.; Trink, A.; Wineberg, Y.; Kanter, I.; Gilad, S.; Pyne, S. A brief review of single-cell transcriptomic technologies. Brief. Funct. Genom. 2018, 17, 64–76. [Google Scholar] [CrossRef]
  40. Bradley, B.T.; Maioli, H.; Johnston, R.; Chaudhry, I.; Fink, S.L.; Xu, H.; Najafian, B.; Marshall, D.; Lacy, J.M.; Williams, T.; et al. Histopathology and ultrastructural findings of fatal COVID-19 infections in Washington State: A case series. Lancet 2020, 396, 320–332. [Google Scholar] [CrossRef]
  41. Liu, Q.; Wang, R.; Qu, G.; Wang, Y.; Liu, P.; Zhu, Y.; Fei, G.; Ren, L.; Zhou, Y.; Liu, L. Gross examination report of a COVID-19 death autopsy. J. Forensic Med. 2020, 36, 21–23. [Google Scholar]
  42. Zhao, Y.; Zhao, Z.; Wang, Y.; Zhou, Y.; Ma, Y.; Zuo, W. Single-Cell RNA Expression Profiling of ACE2, the Receptor of SARS-CoV-2. Am. J. Respir. Crit. Care Med. 2020, 202, 756–759. [Google Scholar] [CrossRef] [PubMed]
  43. Zou, X.; Chen, K.; Zou, J.; Han, P.; Hao, J.; Han, Z. Single-cell RNA-seq data analysis on the receptor ACE2 expression reveals the potential risk of different human organs vulnerable to 2019-nCoV infection. Front. Med. 2020, 14, 185–192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Qi, F.; Qian, S.; Zhang, S.; Zhang, Z. Single cell RNA sequencing of 13 human tissues identify cell types and receptors of human coronaviruses. Biochem. Biophys. Res. Commun. 2020, 526, 135–140. [Google Scholar] [CrossRef] [PubMed]
  45. Wu, C.; Zheng, S.; Chen, Y.; Zheng, M. Single-cell RNA expression profiling of ACE2, the putative receptor of Wuhan 2019-nCoV, in the nasal tissue. medRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  46. Wu, C.; Zheng, M. Single-cell RNA expression profiling shows that ACE2, the putative receptor of Wuhan 2019-nCoV, has significant expression in the nasal, mouth, lung and colon tissues, and tends to be co-expressed with HLA-DRB1 in the four tissues. Preprints 2020. Available online: https://www.preprints.org/manuscript/202002.0247/v1 (accessed on 25 December 2020).
  47. Zhang, H.; Kang, Z.; Gong, H.; Xu, D.; Wang, J.; Li, Z.; Cui, X.; Xiao, J.; Meng, T.; Zhou, W.; et al. The digestive system is a potential route of 2019-nCov infection: A bioinformatics analysis based on single-cell transcriptomes. bioRxiv 2020. [Google Scholar] [CrossRef]
  48. Liu, F.; Long, X.; Zhang, B.; Zhang, W.; Chen, X.; Zhang, Z. ACE2 Expression in Pancreas May Cause Pancreatic Damage After SARS-CoV-2 Infection. Clin. Gastroenterol. Hepatol. 2020, 18, 2128–2130.e2. [Google Scholar] [CrossRef]
  49. Lin, W.; Hu, L.; Zhang, Y.; Ooi, J.D.; Meng, T.; Jin, P.; Ding, X.; Peng, L.; Song, L.; Xiao, Z.; et al. Single-cell analysis of ACE2 expression in human kidneys and bladders reveals a potential route of 2019-nCoV infection. bioRxiv 2020. [Google Scholar] [CrossRef]
  50. He, Q.; Mok, T.N.; Yun, L.; He, C.; Li, J.; Pan, J. Single-cell RNA sequencing analysis of human kidney reveals the presence of ACE2 receptor: A potential pathway of COVID-19 infection. Mol. Genet. Genom. Med. 2020, 8, e1442. [Google Scholar] [CrossRef]
  51. Pan, X.W.; Xu, D.; Zhang, H.; Zhou, W.; Wang, L.H.; Cui, X.G. Identification of a potential mechanism of acute kidney injury during the COVID-19 outbreak: A study based on single-cell transcriptome analysis. Intensive Care Med. 2020, 46, 1114–1116. [Google Scholar] [CrossRef] [Green Version]
  52. Wang, Z.; Xu, X. 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. Cells 2020, 9, 920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zhang, J.; Wu, Y.; Wang, R.; Lu, K.; Tu, M.; Guo, H.; Xie, W.; Qin, Z.; Li, S.; Zhu, P.; et al. Bioinformatic analysis reveals that the reproductive system is potentially at risk from SARS-CoV-2. Preprints 2020. [Google Scholar] [CrossRef]
  54. Zhou, L.; Xu, Z.; Castiglione, G.M.; Soiberman, U.S.; Eberhart, C.G.; Duh, E.J. ACE2 and TMPRSS2 are expressed on the human ocular surface, suggesting susceptibility to SARS-CoV-2 infection. Ocul. Surf. 2020, 18, 537–544. [Google Scholar] [CrossRef] [PubMed]
  55. Seow, J.J.W.; Pai, R.; Mishra, A.; Shepherdson, E.; Lim, T.K.H.; Goh, B.K.P.; Chan, J.K.Y.; Chow, P.K.H.; Ginhoux, F.; DasGupta, R.; et al. scRNA-seq reveals ACE2 and TMPRSS2 expression in TROP2+ Liver Progenitor Cells: Implications in COVID-19 associated Liver Dysfunction. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  56. Lukassen, S.; Chua, R.L.; Trefzer, T.; Kahn, N.C.; Schneider, M.A.; Muley, T.; Winter, H.; Meister, M.; Veith, C.; Boots, A.W.; et al. SARS-CoV-2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells. EMBO J. 2020, 39, e105114. [Google Scholar] [CrossRef]
  57. Madissoon, E.; Wilbrey-Clark, A.; Miragaia, R.J.; Saeb-Parsy, K.; Mahbubani, K.T.; Georgakopoulos, N.; Harding, P.; Polanski, K.; Huang, N.; Nowicki-Osuch, K.; et al. scRNA-seq assessment of the human lung, spleen, and esophagus tissue stability after cold preservation. Genome Biol. 2019, 21, 1. [Google Scholar] [CrossRef] [Green Version]
  58. Han, X.; Zhou, Z.; Fei, L.; Sun, H.; Wang, R.; Chen, Y.; Chen, H.; Wang, J.; Tang, H.; Ge, W.; et al. Construction of a human cell landscape at single-cell level. Nature 2020, 581, 303–309. [Google Scholar] [CrossRef]
  59. Gutierrez, G.D.; Kim, J.; Lee, A.H.; Tong, J.; Niu, J.; Gray, S.M.; Wei, Y.; Ding, Y.; Ni, M.; Adler, C.; et al. Gene Signature of the Human Pancreatic ε Cell. Endocrinology 2018, 159, 4023–4032. [Google Scholar] [CrossRef]
  60. Onorati, M.; Li, Z.; Liu, F.; Sousa, A.M.; Nakagawa, N.; Li, M.; Dell’Anno, M.T.; Gulden, F.O.; Pochareddy, S.; Tebbenkamp, A.T.; et al. Zika Virus Disrupts Phospho-TBK1 Localization and Mitosis in Human Neuroepithelial Stem Cells and Radial Glia. Cell Rep. 2016, 16, 2576–2592. [Google Scholar] [CrossRef] [Green Version]
  61. Zhong, S.; Ding, W.; Sun, L.; Lu, Y.; Dong, H.; Fan, X.; Liu, Z.; Chen, R.; Zhang, S.; Ma, Q.; et al. Decoding the development of the human hippocampus. Nature 2020, 577, 531–536. [Google Scholar] [CrossRef]
  62. Fan, X.; Bialecka, M.; Moustakas, I.; Lam, E.; Torrens-Juaneda, V.; Borggreven, N.V.; Trouw, L.; Louwe, L.A.; Pilgram, G.S.K.; Mei, H.; et al. Single-cell reconstruction of follicular remodeling in the human adult ovary. Nat. Commun. 2019, 10, 3164. [Google Scholar] [CrossRef] [PubMed]
  63. Hu, Z.; Artibani, M.; Alsaadi, A.; Wietek, N.; Morotti, M.; Shi, T.; Zhong, Z.; Gonzalez, L.S.; El-Sahhar, S.; Mallett, G.; et al. The Repertoire of Serous Ovarian Cancer Non-genetic Heterogeneity Revealed by Single-Cell Sequencing of Normal Fallopian Tube Epithelial Cells. Cancer Cell 2020, 37, 226–242.e7. [Google Scholar] [CrossRef] [PubMed]
  64. Guo, J.; Grow, E.J.; Mlcochova, H.; Maher, G.J.; Lindskog, C.; Nie, X.; Guo, Y.; Takei, Y.; Yun, J.; Cai, L.; et al. The adult human testis transcriptional cell atlas. Cell Res. 2018, 28, 1141–1157. [Google Scholar] [CrossRef] [PubMed]
  65. García, S.R.; Deprez, M.; Lebrigand, K.; Cavard, A.; Paquet, A.; Arguel, M.-J.; Magnone, V.; Truchi, M.; Caballero, I.; Leroy, S.; et al. Novel dynamics of human mucociliary differentiation revealed by single-cell RNA sequencing of nasal epithelial cultures. Development 2019, 146, dev177428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Reyfman, P.A.; Walter, J.M.; Joshi, N.; Anekalla, K.R.; McQuattie-Pimentel, A.C.; Chiu, S.; Fernandez, R.; Akbarpour, M.; Chen, C.I.; Ren, Z.; et al. Single-Cell Transcriptomic Analysis of Human Lung Provides Insights into the Pathobiology of Pulmonary Fibrosis. Am. J. Respir. Crit. Care Med. 2019, 199, 1517–1536. [Google Scholar] [CrossRef] [PubMed]
  67. Cui, Y.; Zheng, Y.; Liu, X.; Yan, L.; Fan, X.; Yong, J.; Hu, Y.; Dong, J.; Li, Q.; Wu, X.; et al. Single-Cell Transcriptome Analysis Maps the Developmental Track of the Human Heart. Cell Rep. 2019, 26, 1934–1950.e5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Tang-Huau, T.-L.; Gueguen, P.; Goudot, C.; Durand, M.; Bohec, M.; Baulande, S.; Pasquier, B.; Amigorena, S.; Segura, E. Human in vivo-generated monocyte-derived dendritic cells and macrophages cross-present antigens through a vacuolar pathway. Nat. Commun. 2018, 9, 2570. [Google Scholar] [CrossRef]
  69. Takeda, A.; Hollmén, M.; Dermadi, D.; Pan, J.; Brulois, K.F.; Kaukonen, R.; Lönnberg, T.; Boström, P.; Koskivuo, I.; Irjala, H.; et al. Single-Cell Survey of Human Lymphatics Unveils Marked Endothelial Cell Heterogeneity and Mechanisms of Homing for Neutrophils. Immunity 2019, 51, 561–572.e5. [Google Scholar] [CrossRef]
  70. Han, C.; Duan, C.; Zhang, S.; Spiegel, B.; Shi, H.; Wang, W.; Zhang, L.; Lin, R.; Liu, J.; Ding, Z.; et al. Digestive Symptoms in COVID-19 Patients with Mild Disease Severity: Clinical Presentation, Stool Viral RNA Testing, and Outcomes. Am. J. Gastroenterol. 2020, 115, 916–923. [Google Scholar] [CrossRef]
  71. Lee, I.C.; Huo, T.I.; Huang, Y.H. Gastrointestinal and liver manifestations in patients with COVID-19. J. Chin. Med. Assoc. 2020, 83, 521–523. [Google Scholar] [CrossRef]
  72. Hindson, J. COVID-19: Faecal-oral transmission? Nat. Rev. Gastroenterol. Hepatol. 2020, 17, 259. [Google Scholar] [CrossRef] [PubMed]
  73. Amirian, E.S. Potential fecal transmission of SARS-CoV-2: Current evidence and implications for public health. Int. J. Infect. Dis. 2020, 95, 363–370. [Google Scholar] [CrossRef] [PubMed]
  74. Pennisi, M.; Lanza, G.; Falzone, L.; Fisicaro, F.; Ferri, R.; Bella, R. SARS-CoV-2 and the Nervous System: From Clinical Features to Molecular Mechanisms. Int. J. Mol. Sci. 2020, 21, 5475. [Google Scholar] [CrossRef] [PubMed]
  75. Toptan, T.; Aktan, Ç.; Başarı, A.; Bolay, H. Case Series of Headache Characteristics in COVID-19: Headache Can Be an Isolated Symptom. Headache 2020, 60, 1788–1792. [Google Scholar] [CrossRef] [PubMed]
  76. Tsatsakis, A.; Calina, D.; Falzone, L.; Petrakis, D.; Mitrut, R.; Siokas, V.; Pennisi, M.; Lanza, G.; Libra, M.; Doukas, S.G.; et al. SARS-CoV-2 pathophysiology and its clinical implications: An integrative overview of the pharmacotherapeutic management of COVID-19. Food Chem. Toxicol. 2020, 146, 111769. [Google Scholar] [CrossRef] [PubMed]
  77. Santoriello, D.; Khairallah, P.; Bomback, A.S.; Xu, K.; Kudose, S.; Batal, I.; Barasch, J.; Radhakrishnan, J.; D’Agati, V.; Markowitz, G. Postmortem Kidney Pathology Findings in Patients with COVID-19. J. Am. Soc. Nephrol. 2020, 31, 2158–2167. [Google Scholar] [CrossRef] [PubMed]
  78. Driggin, E.; Madhavan, M.V.; Bikdeli, B.; Chuich, T.; Laracy, J.; Biondi-Zoccai, G.; Brown, T.S.; der Nigoghossian, C.; Zidar, D.A.; Haythe, J.; et al. Cardiovascular Considerations for Patients, Health Care Workers, and Health Systems During the COVID-19 Pandemic. J. Am. Coll. Cardiol. 2020, 75, 2352–2371. [Google Scholar] [CrossRef]
Ijerph 18 00284 g001 550
Figure 1. High angiotensin-converting enzyme II (ACE2) and transmembrane serine protease 2 (TMPRSS2) expression levels of type II alveolar (AT2) cells in the lung. (A), (B) R A C E 2 and R T M P R S S 2 in AT2 cells of 8 samples, respectively. The black dotted lines represent the corresponding average ratios across 8 samples. (C), (D) The expression distribution of ACE2 and TMPRSS2 in AT2 cells across 8 samples, respectively. The red dotted lines represent the corresponding average expression values of gene-expressed cells across 8 samples.
Figure 1. High angiotensin-converting enzyme II (ACE2) and transmembrane serine protease 2 (TMPRSS2) expression levels of type II alveolar (AT2) cells in the lung. (A), (B) R A C E 2 and R T M P R S S 2 in AT2 cells of 8 samples, respectively. The black dotted lines represent the corresponding average ratios across 8 samples. (C), (D) The expression distribution of ACE2 and TMPRSS2 in AT2 cells across 8 samples, respectively. The red dotted lines represent the corresponding average expression values of gene-expressed cells across 8 samples.
Ijerph 18 00284 g001
Ijerph 18 00284 g002 550
Figure 2. High ACE2 and TMPRSS2 expression levels of the primordium cells in the gall bladder. (A) Uniform manifold approximation and projection (UMAP) visualization of clustering results for gall bladder cells. (B) R A C E 2 in each cell cluster. (C) R T M P R S S 2 in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster.
Figure 2. High ACE2 and TMPRSS2 expression levels of the primordium cells in the gall bladder. (A) Uniform manifold approximation and projection (UMAP) visualization of clustering results for gall bladder cells. (B) R A C E 2 in each cell cluster. (C) R T M P R S S 2 in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster.
Ijerph 18 00284 g002
Ijerph 18 00284 g003 550
Figure 3. High ACE2 and TMPRSS2 expression levels of epithelial cells and ovarian somatic cells in the fallopian tube. (A) UMAP visualization of clustering results for the fallopian tube cells. (B) R A C E 2 in each cell cluster. (C) R T M P R S S 2 in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster.
Figure 3. High ACE2 and TMPRSS2 expression levels of epithelial cells and ovarian somatic cells in the fallopian tube. (A) UMAP visualization of clustering results for the fallopian tube cells. (B) R A C E 2 in each cell cluster. (C) R T M P R S S 2 in each cell cluster. (D) ACE2 expression level in each cell cluster on the UMAP plot. (E) TMPRSS2 expression level in each cell cluster on the UMAP plot. (F) The expression distribution of ACE2 across each cell cluster. (G) The expression distribution of TMPRSS2 across each cell cluster.
Ijerph 18 00284 g003
Ijerph 18 00284 g004 550
Figure 4. The risk levels of susceptible organs. The dot plot illustrates the rank of R A C E 2 , R T M P R S S 2 and R of 16 cell types. The sizes of the dot indicate the value of the ratio, and the colors indicate the rank of the ratio.
Figure 4. The risk levels of susceptible organs. The dot plot illustrates the rank of R A C E 2 , R T M P R S S 2 and R of 16 cell types. The sizes of the dot indicate the value of the ratio, and the colors indicate the rank of the ratio.
Ijerph 18 00284 g004
Table
Table 1. Human organs available for single-cell RNA sequencing (scRNA-seq) data analysis.
Table 1. Human organs available for single-cell RNA sequencing (scRNA-seq) data analysis.
Human SystemsHuman Organs
Digestive systemEsophagus [57]
Small intestine (jejunum, ileum, andduodenum) [58]
Large intestine (rectum, colon)
Stomach [58]
Liver [58]
Gall bladder [58]
Pancreatic islets [59]
Nervous systemBrain (substantia nigra and cortex, neuronal epithelium, and hippocampus) [60,61]
Cerebellum [58]
Spinal cord [58]
Reproductive systemOvary [62]
Fallopian tube [63]
Uterus [58]
Testis [64]
Motor systemMuscle [58]
Respiratory systemNose (nasal brushing epithelial cells, nasal turbinate epithelial cells, and nasal airway epithelium) [65]
Bronchus [65]
Lung [66]
Trachea [58]
Circulatory systemPeripheral blood [58]
Heart [67]
Artery [58]
Spleen [57]
Urinary systemKidney [58]
Ureter [58]
Prostate [58]
Immune systemTonsil (tonsil dendritic cells) [68]
Bone marrow [58]
Lymph nodes [69]
Endocrine systemThyroid [58]
Thymus [58]
Table
Table 2. Cell types with high expression levels of ACE2 and TMPRSS2.
Table 2. Cell types with high expression levels of ACE2 and TMPRSS2.
SystemsOrgansCell Types R A C E 2 R T M P R S S 2
Respiratory systemLungAT2 cells0.79%21.20%
Nose (nasal turbinate epithelial cells)Mesenchymal stromal cells1.70%8.50%
Plasma cells2.00%8.80%
Nose (nasal brushing epithelial cells)Plasma cells4.00%13.80%
Nose (nasal airway epithelial cells)Nasal airway epithelial cells8.40%23.10%
Digestive systemGall bladderPrimordium cell2.20%2.60%
Small intestine (jejunum)Enterocyte progenitor cell14.00%6.60%
Goblet cell9.10%2.60%
Small intestine (ileum)Intestinal epithelial stem cell1.50%4.20%
Enterocyte progenitor cell2.40%3.60%
Small intestine (duodenum)LGR5+ stem cell5.20%4.50%
Intestinal epithelial stem cell3.90%5.80%
Enterocyte progenitor cell4.40%6.80%
Tuft progenitor cell7.40%6.60%
Enteroendocrine cell9.40%3.80%
Large intestine (rectum)Goblet progenitor cell2.80%62.50%
MKI67+ progenitor cell10.30%75.40%
Enterocyte13.20%77.50%
Goblet cell18.20%87.30%
Large intestine (colon)Enterocyte5.70%52.10%
Goblet cell11.00%78.00%
EsophagusSecretory progenitor cell9.40%18.80%
Nervous systemBrain (substantia nigra and cortex)Oligodendrocyte precursor cell1.60%0.90%
Astrocyte1.90%1.00%
Reproductive systemFallopian tubeEpithelial cells1.40%26.50%
Ovarian somatic cell1.00%1.00%
TestisSpermatogonium1.70%4.50%
Peritubular myoid cell1.70%6.30%
Testis somatic cell2.10%4.30%
Spermatogonial stem cell1.40%2.10%
Circulatory systemHeartCardiomyocyte6.60%0.80%
Cardiovascular progenitor cell12.50%0.40%
Urinary systemKidneyNephron epithelial cell2.70%10.70%
Epithelial cell2.70%9.60%
Endothelial cell2.70%12.80%
Mesangial cell3.00%14.50%
Table
Table 3. Risk level of organs to SARS-CoV-2 infection.
Table 3. Risk level of organs to SARS-CoV-2 infection.
OrgansSystemsCell TypesRRisk Levels
Large intestine (rectum)Digestive systemGoblet cell0.3981
Large intestine (colon)Digestive systemGoblet cell0.2901
Nose (nasal airway epithelium)Respiratory systemAll cells0.1401
EsophagusDigestive systemSecretory progenitor cell0.1331
Small intestine (jejunum)Digestive systemEnterocyte progenitor cell0.0962
KidneyUrinary systemMesangial cell0.0662
Fallopian tubeReproductive systemEpithelial cell0.0602
Small intestine (duodenum)Digestive systemEnterocyte progenitor cell0.0552
Nose (nasal brushing epithelial cells)Respiratory systemPlasma cell0.0522
Nose (nasal turbinate epithelial cells)Respiratory systemPlasma cell0.0423
LungRespiratory systemAT20.0413
TestisReproductive systemPeritubular myoid cell0.0333
Small intestine (ileum)Digestive systemIntestinal epithelial stem cell0.0253
Gall bladderDigestive systemPrimordium cell0.0243
HeartCirculatory systemCardiomyocyte0.0233
Brain (substantia nigra and cortex)Nervous systemAstrocyte0.0143
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qi, J.; Zhou, Y.; Hua, J.; Zhang, L.; Bian, J.; Liu, B.; Zhao, Z.; Jin, S. The scRNA-seq Expression Profiling of the Receptor ACE2 and the Cellular Protease TMPRSS2 Reveals Human Organs Susceptible to SARS-CoV-2 Infection. Int. J. Environ. Res. Public Health 2021, 18, 284. https://doi.org/10.3390/ijerph18010284

AMA Style

Qi J, Zhou Y, Hua J, Zhang L, Bian J, Liu B, Zhao Z, Jin S. The scRNA-seq Expression Profiling of the Receptor ACE2 and the Cellular Protease TMPRSS2 Reveals Human Organs Susceptible to SARS-CoV-2 Infection. International Journal of Environmental Research and Public Health. 2021; 18(1):284. https://doi.org/10.3390/ijerph18010284

Chicago/Turabian Style

Qi, Jing, Yang Zhou, Jiao Hua, Liying Zhang, Jialin Bian, Beibei Liu, Zicen Zhao, and Shuilin Jin. 2021. "The scRNA-seq Expression Profiling of the Receptor ACE2 and the Cellular Protease TMPRSS2 Reveals Human Organs Susceptible to SARS-CoV-2 Infection" International Journal of Environmental Research and Public Health 18, no. 1: 284. https://doi.org/10.3390/ijerph18010284

APA Style

Qi, J., Zhou, Y., Hua, J., Zhang, L., Bian, J., Liu, B., Zhao, Z., & Jin, S. (2021). The scRNA-seq Expression Profiling of the Receptor ACE2 and the Cellular Protease TMPRSS2 Reveals Human Organs Susceptible to SARS-CoV-2 Infection. International Journal of Environmental Research and Public Health, 18(1), 284. https://doi.org/10.3390/ijerph18010284

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop