Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Pericytes may facilitate SARS-CoV-2 entry into the nervous system

Pericytes may facilitate SARS-CoV-2 entry into the nervous system Objective: Although the neurological and olfactory symptoms of coronavirus disease 2019 have been identified, the neurotropic properties of the causative virus, severe acute respiratory syndrome-associated coronavirus 2 (SARS-CoV-2), remain unknown. We sought to identify the susceptible cell types and potential routes of SARS-CoV-2 entry into the central nervous system, olfactory system, and respiratory system. Methods: We collected single-cell RNA data from normal brain and nasal epithelium specimens, along with bronchial, tracheal, and lung specimens in public datasets. The susceptible cell types that express SARS-CoV-2 entry genes were identified using single-cell RNA sequencing and the expression of the key genes at protein levels was verified by immunohistochemistry. We compared the coexpression patterns of the entry receptor angiotensin-converting enzyme 2 (ACE2) and the spike protein priming enzyme transmembrane serine protease (TMPRSS)/cathepsin L among the specimens. Results: The SARS-CoV-2 entry receptor ACE2 and the spike protein priming enzyme TMPRSS/cathepsin L were coexpressed by pericytes in brain tissue; this coexpression was confirmed by immunohistochemistry. In the nasal epithelium, ciliated cells and sustentacular cells exhibited strong coexpression of ACE2 and TMPRSS. Neurons and glia in the brain and nasal epithelium did not exhibit coexpression of ACE2 and TMPRSS. However, coexpression was present in ciliated cells, vascular smooth muscle cells, and fibroblasts in tracheal tissue; ciliated cells and goblet cells in bronchial tissue; and alveolar epithelium type 1 cells, AT2 cells, and ciliated cells in lung tissue. Conclusion: Neurological symptoms in patients with coronavirus disease 2019 could be associated with SARS-CoV-2 invasion across the blood–brain barrier via pericytes. Additionally, SARS-CoV-2–induced olfactory disorders could be the result of localized cell damage in the nasal epithelium. Keywords: ACE2, coronavirus disease 2019, SARS-CoV-2, single-cell RNA sequencing, TMPRSS2 ZK, JW, and TM contributed equally to the article. Introduction Department of Rheumatology and Immunology, Changzheng Hospital, The main clinical manifestations of coronavirus disease 2019 Department of Neurosurgery, Changhai Hospital, Second Military Medical (COVID-19) resulting from severe acute respiratory syndrome University, Division of Spine, Department of Orthopedics, Shanghai General coronavirus 2 (SARS-CoV-2) infection are fever, cough, and dys- Hospital, School of Medicine, Shanghai Jiaotong University, Tongji University [1–5] e pnea. However, numerous neurological abnormalities have Cancer Center, School of Medicine, Tongji University, Department of Orthopaedic been described in patients with COVID-19; these abnormalities Oncology, Changzheng Hospital, Department of Urology, The Third Affiliated Hospital, Second Military Medical University, Department of Pathology, Shanghai include olfactory and gustatory disorders (anosmia, hyposmia, Tenth People’s Hospital, Tongji University School of Medicine, Shanghai, China, ageusia, and dysgeusia), headache, epilepsy, encephalitis, and Zeenat Qureshi Stroke Institutes and University of Missouri, Columbia, MO, [4–8] acute myelitis. SARS-CoV-2 RNA has also been identified in i j USA, Peking-Tsinghua Center for Life Sciences, School of Medicine, Tsinghua [9,10] cerebrospinal fluid via genome sequencing. These findings University, Beijing, China suggest that SARS-CoV-2 can directly or indirectly enter the *Corresponding authors: Huji Xu, Department of Rheumatology and Immunology, central nervous system (CNS). However, the underlying mecha- Changzheng Hospital, Second Military Medical University, Shanghai 200003, nism remains unknown. China. E-mail: xuhuji@smmu.edu.cn; Jianmin Liu, Department of Neurosurgery, Changhai Hospital, Second Military Medical University, Shanghai 200003, China. The CNS can be invaded through the hematogenous route or [11] E-mail: chstroke@163.com by axonal transport ; these routes involve infections of endo- thelial cells in the blood–brain barrier (BBB) or epithelial cells Copyright © 2023 The Chinese Medical Association, Published by Wolters Kluwer Health, Inc. This is an open-access article distributed under the terms of in the blood–cerebrospinal fluid barrier, as well as infections [12] the Creative Commons Attribution-Non Commercial-No Derivatives License 4.0 of the olfactory bulb through the olfactory nerve. The tro- (CCBY-NC-ND), where it is permissible to download and share the work provided pism of coronaviruses is primarily dependent on the presence of it is properly cited. The work cannot be changed in any way or used commercially attachment and entry receptors, along with protease availability without permission from the journal. and the activities of internalization-related and trafficking-re- Journal of Bio-X Research (2023) 6:23–36 [13,14] lated pathways. SARS-CoV-2 uses angiotensin-convert- [15,16] Received: 4 May 2022; Accepted: 7 December 2022 ing enzyme 2 (ACE2) as an entry receptor. Other entry receptor molecules associated with SARS-CoV-2 may include http://dx.doi.org/10.1097/JBR.0000000000000137 23 ReseaRch aR ticle Journal of Bio-X Research [17,18] CD147 and CD209. Upon receptor binding, spike protein Cell type identification and gene expression analysis priming enzymes (eg, cell surface transmembrane serine pro- Clusters were annotated based on the expression patterns of tease [TMPRSS] and endosomal cathepsins) are necessary for known cellular markers; the clustering information was visu- [19,20] some coronaviruses to achieve membrane fusion. Thus, the alized by uniform manifold approximation and projection detection of coexpression involving key genes in the viral entry (UMAP) plots. Then, lost values in the gene expression matrix process may help to identify cells that are susceptible to SARS- were imputed by the “RunALRA” function. The imputed gene CoV-2 infection; this will contribute to further investigations of expression values are displayed in feature and violin plots. the mechanisms by which SARS-CoV-2 enters the CNS. The mean expression levels of gene sets are displayed in box In this study, we used single-cell RNA sequencing (scRNA- plots. For comparisons of expression levels among datasets, seq) of normal brain and nasal epithelium specimens, as well as the quantile method was used to correct expression data; this tracheal, bronchial, and lung specimens, to identify susceptible [29,30] global normalization method is widely used. To calculate cell types and potential routes of SARS-CoV-2 entry into the the endocytosis signaling score, we obtained endocytosis-as- CNS, olfactory system, and respiratory system. We sought to sociated genes from “GOBP_ENDOCYTOSIS” gen set in the clarify how SARS-CoV-2 invades various tissues in the human MigDB dataset (https://www.gsea-msigdb.org). The average body, with a particular interest in elucidating the olfactory and expression level of all the genes in this gene set was considered CNS symptoms in patients with COVID-19. as the expression level of the endocytosis signaling. Materials and methods Immunohistochemical staining Data sources Human brain tissue specimens were obtained from five patients We analyzed single-cell transcriptomes of the brain, nasal olfac- who had undergone surgery for treatment of traumatic brain tory epithelium (OE), nasal respiratory epithelium (RE), lung, injury, including three males and two females. Human lung tis- tracheal epithelium, and bronchial epithelium (Additional Table sue specimens were obtained during routine surgical intervention 1, http://links.lww.com/JR9/A46). Nasal RE/OE and trachea from three lung cancer patients, including two males and one transcriptome data were obtained from the Gene Expression female. Ethical approval was granted by the Ethics Committee Omnibus (GEO; https://www.ncbi.nlm.nih.gov/), using acces- of Shanghai Changhai Hospital (Approval No. CHEC2020- [21] [22] sion codes GSE139522 and GSE134355. Brain transcrip- 164). All patients were fully informed of the research and signed tome data were obtained from the Sequence Read Archive the informed consent forms. Formalin-fixed paraffin-embedded (SRA; https://www.ncbi.nlm.nih.gov/sra), using accession code brain tissue sections on glass slides were incubated at 37°C over- [23] SRA667466. Lung transcriptome data were obtained from night, deparaffinized in xylene, and rehydrated in a descending Cell Atlas Data Portal (https://data.humancellatlas.org), using ethanol gradient. Then, tissue sections were incubated in antigen [24] accession code PRJEB31843. Bronchial transcriptome data retrieval solution at 95°C for 30 minutes. Next, tissue sections [25] were obtained from research data published by Lukassen et al. were incubated in 3% hydrogen peroxide for 10 minutes, fol- The single-cell data within each dataset had similar sequencing lowed by serum-free protein blocking solution (Cat# 9048-46- depth to ensure the quality of consensus sequence. Because the 8, Sigma, St. Louis, MO, USA) for 20 minutes; they were then data were obtained from public datasets, approval for our study incubated with rabbit anti-human ACE2 primary antibody (Cat# by the ethics committee and patient consent were unnecessary. ab108252, Abcam, Cambridge, UK) at a dilution ratio of 1:200 overnight at 4°C temperature and goat anti-rabbit IgG second- ary antibody (Cat# ab6721, Abcam) at a dilution ratio of 1:1000 Quality control for 1 hour at room temperature. After incubation with primary Barcodes with outlier gene counts may be indicative of dying and secondary antibodies, color development was performed cells, cells with broken membranes, or doublets. A high frac- using 3-amino-9-ethylcarbazole as the chromogenic substrate tion of mitochondria is indicative of leaky cytoplasmic mem- and tissue was counterstained with hematoxylin as described [26] [31] branes allowing mRNA to escape. We removed low-quality previously. The microscopic observation was performed by data from cells with <200 or >5000 expressed genes, as well as the optical microscope (Eclipse 80i, Nikon, Tokyo, Japan). data from cells in which >10% of unique molecular identifiers mapped to mitochondrial genes. Results Single-cell maps of mouse brain, healthy human Data integration, dimension reduction, and cell olfactory system, and healthy human respiratory system clustering [27] [28] To identify susceptible cell types in the brain, olfactory system, The R package Seurat and R package Harmony were and respiratory system, we collected scRNA-seq data for six used to process scRNA-seq data. Unique molecular identifier types of tissue from published papers (Additional Table 1, http:// counts were normalized using the “NormalizeData” function. The top 1000 highly variable genes were identified using the links.lww.com/JR9/A46). In the brain, 11 cell types were identi- “FindVariableGenes” function. Principal component analy- fied based on the analysis of 29,678 cells, after the application sis was performed on single-cell expression matrices using the of quality control criteria (Fig. 1A). The main cell types were oli- “RunPCA” function. Subsequently, “Runharmony” was used for godendrocytes, interneurons, astrocytes, microglia, and neurons batch-effect correction in each dataset. Cell clustering analysis 1 and 2. The other identified cell types included ependymal cells, was performed after the construction of a K-nearest-neighbor endothelial cells, pericytes, dividing cells, and myeloid cells. graph based on the adjusted principal component analysis The nasal epithelium was divided into two types: OE and RE. embeddings by harmony algorithm. In total, 19,903 cells and 14 cell types were identified in the OE 24 Journal of Bio-X Research ReseaRch aR ticle Figure 1. Single-cell maps of healthy human respiratory system and mouse brain tissue. UMAP plots showing the distributions of brain cells (A), nasal olfactory epithelial cells (B), nasal respiratory epithelial cells (C), tracheal epithelial cells (D), bronchial epithelial cells (E), and lung cells (F). Single-cell map of respiratory system (G). All the single-cell data in the human respiratory system were combined for another round of reduction and clustering analysis. UMAP plots show the distributions of the cell clusters in the human respiratory system. DCs=dendritic cells, HBCs=horizontal basal cells, NK=natural killer, UMAP=uniform manifold approximation and projection, VSMC=vascular smooth muscle cells. (Fig. 1B); 10,657 cells and 12 cell types were identified in the RE as well as chondrocytes, VSMCs, endothelial cells, fibroblasts, (Fig. 1C). Epithelial cell types in the RE were horizontal basal and immune cells (ie, myeloid and B cells). Additionally, 17,228 cells (HBCs) and Bowman’s gland cells, as well as sustentac- cells and seven cell types were identified in bronchial specimens ular, ciliated, and endodermal progenitor cells. Non-epithelial after the application of quality control criteria (Fig.  1E). The cell types in the RE were microvillar and endothelial cells, as main cell types identified were ciliated, basal, and goblet cells. well as fibroblasts and immune cells (ie, mast, myeloid, T, and B The other identified cell types consisted of brush cells, fibro- cells). Epithelial cell types in the OE were HBCs, pericytes, and blasts, dividing cells, and ionocytes. Bowman’s gland cells, as well as gland progenitor, sustentacular, After the application of quality control criteria, 24,710 cells and ciliated cells. Non-epithelial cell types in the OE were neu- and 14 cell types were identified in lung specimens (Fig.  1F). rons, olfactory ensheathing glia cells (OECs), fibroblasts, vascu- The cell types identified were alveolar epithelium type 1 (AT1), lar smooth muscle cells (VSMCs), and immune cells (ie, mast, alveolar epithelium type 2 (AT2), and ciliated cells, as well as myeloid, T, and B cells). fibroblasts, lymph vessel cells, VSMCs, and endothelial cells. After the application of quality control criteria, 9614 cells and The following immune cell types were identified: mast, T, B, and 10 cell types were identified in tracheal specimens (Fig. 1D). The natural killer (NK) cells, along with monocytes, macrophages, cell types were goblet, brush, ciliated, basal, and granular cells, and dendritic cells. 25 ReseaRch aR ticle Journal of Bio-X Research Because of the similar cells present among tissues in the respi- progenitor, ciliated, and sustentacular cells, as well as neu- ratory system, we constructed a map of the respiratory system rons; its expression was minimal in fibroblasts, HBCs, and based on olfactory epithelum data, nasal respiratory epithelium Bowman’s gland cells (Fig.  3C and D). CTSL was expressed data, trachea data, bronchus data, and lung data from 83,126 in all non-immune subclusters, but it was only expressed in cells. In this map, we identified 17 subclusters: HBCs, endo- myeloid cells among the immune subclusters (Fig. 3C and D). thelial cells, Bowman’s gland cells, sustentacular cells, alveolar Additionally, endocytosis signaling was strongly expressed cells, ciliated cells, neurons, OECs, fibroblasts, microvillar cells, in non-immune subclusters and weakly expressed in immune VSMCs, chondrocytes, basal cells, mast cells, myeloid cells, B subclusters (Fig.  3C and D). PIKFYVE and TPCN2 were cells, and T/NK cells (Fig. 1G). expressed in most cell types (Fig.  3D). The viral RNA repli- cation target gene MTHFD1 was detected in all subclusters in the OE, except T and B cells (Fig.  3D). Overall, the data Coexpression patterns of key entry genes in mouse brain indicated that multiple cell types in the OE could serve as tar- tissue get cells for SARS-CoV-2. The gland progenitor, ciliated, and To identify susceptible cell types in brain tissue and elucidate sustentacular cells may constitute the most susceptible cells potential routes of SARS-CoV-2 entry, we collected scRNA-seq in the OE because they exhibited the strongest expression of data from mouse brains. In total, 11 subclusters were identified ACE2 and TMPRSS2 (Fig. 3C and D). Thus, we speculate that in mouse brain cells (Fig. 2A); marker genes for each subclus- non-neuronal cell types, rather than sensory or bulb neurons, ter are shown in Figure  2B. ACE2 was strongly expressed in are responsible for anosmia and associated odor perception pericytes and weakly expressed in other subclusters (Fig. 2C). deficits in patients with COVID-19. To confirm the expression of ACE2 at the protein level, we performed immunohistochemical analysis of normal human Coexpression patterns of key entry genes in healthy brain tissues; the immunohistochemistry results were consistent human nasal respiratory epithelial tissue with the scRNA-seq findings. ACE2 was strongly expressed in vascular pericytes, as well as endothelial cells (Fig. 2D). Lung RE is a major component of the nasal epithelium. We identified alveolar cells also exhibited strong expression of ACE2 at the 12 subclusters (Fig.  4A) and the marker genes for each sub- protein level (Fig. 2E); they are regarded as the main target of cluster are shown in Figure 4B. ACE2 was strongly expressed [32,33] SARS-CoV-2 lung tissue. Therefore, pericytes may be the in sustentacular, ciliated, and microvillar cells; it was weakly main target of SARS-CoV-2 in brain tissue; entry through these expressed in endodermal progenitor cells (Fig.  4C and D). cells might be involved in indirect damage to neurons. CD147, CD147 was expressed in almost all subclusters, whereas a potential entry receptor, was expressed in almost all sub- CD209 was expressed in myeloid cells (Fig.  4D). TMPRSS clusters (Fig. 2F). The spike protein priming enzyme TMPRSS and CTSL had similar expression patterns among subclusters was nearly absent from all subclusters in brain tissue, whereas in the RE; they did not exhibit robust expression in immune cathepsin L (CTSL) was expressed in all subclusters (Fig.  2C cells (Fig. 4C and D). Thus, multiple cell types (eg, microvillar, and F); these findings suggested that spike protein hydrolysis ciliated, and sustentacular) expressed ACE2, TMPRSS/CTSL, in brain tissue is mediated by CTSL. Because endocytosis is an endocytosis signaling, PIKFYVE, TPCN2, and MTHFD1; important aspect of coronavirus infection, we examined the these cells were identified as potential targets of SARS-CoV-2 expression of the endocytosis signaling molecules phospho- (Fig. 4D). inositide kinase, FYVE-type zinc finger containing (PIKFYVE) and two pore segment channel 2 (TPCN2). We found that Coexpression patterns of key entry genes in healthy PIKFYVE was strongly expressed in ependymal cells, peri- human tracheal and bronchial tissues cytes, and neurons 1 and 2 in brain tissue (Fig.  2F), whereas the expression of TPCN2 was minimal and primarily observed Among the 11 cell types identified in tracheal tissue (Fig. 5A in myeloid cells. Furthermore, methylenetetrahydrofolate dehy- and B), ACE2 was expressed in fibroblasts, VSMCs, and cili- drogenase 1 (MTHFD1), which is an essential gene for RNA ated cells (Fig.  5C and D). CD147 was expressed in all sub- [32] virus replication, was expressed in all subclusters in brain clusters, whereas CD209 was only expressed in myeloid cells tissue (Fig.  2F). These results suggested that neurons do not (Fig.  5D). TMPRSS was mainly expressed in goblet, granule, constitute a primary site of infection; however, pericytes are the brush, basal, and ciliated cells; in contrast, CTSL was expressed most susceptible cells in brain tissue because they exhibit strong in chondrocytes, endothelial cells, VSMCs, fibroblasts, myeloid expression of ACE2 and coexpression of CTSL, which facilitate cells, and B cells (Fig.  5C and D). Endocytosis signaling, viral entry and replication. PIKFYVE, TPCN2, and MTHFD1 were strongly expressed in all cell types in tracheal tissue (Fig. 5D). Among the seven sub- clusters in bronchial tissue (Fig. 6A and B), ACE2 was mainly Coexpression patterns of key entry genes in healthy expressed in goblet and ciliated cells, whereas CD147 was human nasal olfactory epithelial tissue expressed in almost all subclusters (Fig. 6C). TMPRSS, CTSL, Next, we explored whether SARS-CoV-2 could enter the CNS endocytosis signaling, PIKFYVE, TPCN2, and MTHFD1 were through the olfactory system. We identified 14 subclusters in expressed in all cell types in bronchial tissue (Fig. 6C and D). the OE (Fig.  3A); the marker genes for each subcluster are Taken altogether, the results indicated that ciliated cells are shown in Figure  3B. ACE2 was strongly expressed in suste- the most susceptible cells in both tracheal and bronchial tissue ntacular, ciliated, and gland progenitor cells; its expression because they generally exhibit strong expression of ACE2 and was weaker in HBCs and Bowman’s gland cells (Fig. 3C and TMPRSS (Figs.  5D and 6C). Goblet cells also demonstrated D). CD147 was expressed in almost all subclusters, whereas susceptibility to SARS-CoV-2 through its coexpression of CD209 was only expressed in myeloid cells (Fig. 3D). TMPRSS ACE2, TMPRSS/CTSL, endocytosis signaling, and MTHFD1 was expressed in many subclusters in the OE, including gland (Figs. 5D and 6C). 26 Journal of Bio-X Research ReseaRch aR ticle Figure 2. Single-cell analysis of mouse brain cells. UMAP plot showing the distribution of mouse brain cells (A). Dot plot showing the marker genes for brain clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots show- ing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. Representative immunohistochemical staining of ACE2 in cerebral vascular pericytes (D) and lung tissue of humans (E). Arrows indicate ACE2-postive cells. Violin and box plots for the expression levels of ACE2, CD147, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (F). Gene expression levels are measured as log2 (TP10K+1) values and the expression levels of gene sets are measured as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projection. 27 ReseaRch aR ticle Journal of Bio-X Research Figure 3. Single-cell analysis of healthy human nasal olfactory epithelial cells. UMAP plot showing the distribution of nasal olfactory epithelial cells (A). Dot plot showing the marker genes for nasal olfactory epithelial clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indi- cates high expression of TMPRSS or CTSL. Violin and box plots for the expression levels of ACE2, CD147, CD209, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (D). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, HBCs=horizontal basal cells, OEC=olfactory ensheathing glia cells, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projection, VSMC=vascular smooth muscle cells. Coexpression patterns of key entry genes in healthy PIKFYVE, TPCN2, and MTHFD1 were strongly expressed in lung cells (Fig. 7D). human lung tissue In total, 14 subclusters (Fig.  7A) in lung tissue were identi- fied on the basis of cell type-specific marker genes (Fig.  7B). Expression patterns of SARS-CoV-2 and neurotropic ACE2 was mainly expressed in AT2 cells; it was also expressed virus receptor genes in healthy human respiratory and in AT1 and ciliated cells (Fig. 7C and D). CD209 was mainly mouse brain tissues expressed in monocytes and macrophages, whereas CLEC4M For comprehensive analysis of susceptible cell types in the respi- was expressed in endothelial cells. CD147 was expressed in ratory system and comparison of similarities among cell types, all subclusters in lung tissue (Fig.  7D). TMPRSS and CTSL we constructed a map of the respiratory system and identified 17 were expressed in most cell types; they were not expressed in subclusters; the main subclusters comprised myeloid, alveolar, mast, NK, T, or B cells (Fig. 7C and D). Endocytosis signaling, 28 Journal of Bio-X Research ReseaRch aR ticle Figure 4. Single-cell analysis of healthy human nasal respiratory epithelial cells. UMAP plot showing the distribution of nasal respiratory epithelial cells (A). Dot plot showing the marker genes for nasal respiratory epithelial clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. (D) Violin and box plots for the expression levels of ACE2, CD147, CD209, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters. Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, HBCs=horizontal basal cells, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projection. T/NK, endothelial, ciliated, neuronal, and B cells (Additional types in the respiratory system (Fig.  8A). The endocytosis-re- Fig. 1A and B, http://links.lww.com/JR9/A45). We found that lated gene PIKFYVE was expressed mainly in myeloid cells and ACE2 was strongly expressed in ciliated cells in the OE, RE, tra- fibroblasts in the RE, as well as endothelial cells in lung tissue cheal tissue, and bronchial tissue; sustentacular cells in the OE (Fig. 8A). These results showed that factors involved in SARS- and RE; gland progenitor cells in the RE; and AT2 cells in lung CoV-2 entry are strongly expressed in the OE and RE, highlight- tissue (Fig. 8A). CD209 was strongly expressed in myeloid cells ing the potential role of such cells to facilitate viral infection in the OE, RE, and tracheal tissue. TMPRSS, CTSL, endocytosis and spread. The most susceptible cells in the respiratory system signaling, TPCN2, and MTHFD1 were expressed in most cell were ciliated, gland progenitor, and sustentacular cells, as well 29 ReseaRch aR ticle Journal of Bio-X Research Figure 5. Single-cell analysis of healthy human trachea cells. UMAP plot showing the distribution of tracheal cells (A). Dot plot showing the marker genes for tracheal clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. Violin and box plots for the expression levels of ACE2, CD147, CD209, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (D). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projec- tion, VSMC=vascular smooth muscle cells. as HBCs; AT2 cells are the known vulnerable cells for SARS- A45). Notably, many viral receptors (eg, AXL, ISG15, ITGA4, [32,33] COV2, their susceptibility could be greater than the suscep- and PECAM1) are expressed in pericytes, suggesting that peri- tibility of AT2 cells, indicating that these cells could serve as the cytes mediate CNS infection (Additional Fig. 3, http://links. portal of entry for SARS-CoV-2. lww.com/JR9/A45). Furthermore, because exosome biogenesis To explore the potential mechanism of SARS-CoV-2 infec- and release are also associated with viral transport, we explored tion in the CNS, we compared the expression patterns of entry the brain-specific and OE-specific expression patterns of genes receptors for various neurotropic viruses between maps of the involved in exosome biogenesis and release. The results showed respiratory system and brain. We found that some neurotropic that neurons, ciliated cells, and OECs had strong capacities for exosome biogenesis and release, which may serve as alterna- virus entry receptors were not expressed in neurons, indicating tive pathways for CNS infection (Fig. 8B). Thus, we speculated that other mechanisms (eg, other cells) were involved in CNS the potential mechanisms involved in CNS infection (Fig. 8C). infection (Additional Figs. 2 and 3, http://links.lww.com/JR9/ 30 Journal of Bio-X Research ReseaRch aR ticle Figure 6. Single-cell analysis of healthy human bronchial epithelial cells. UMAP plot showing the distribution of bronchial epithelial cells (A). Dot plot showing the marker genes for bronchial epithelial clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. Violin and box plots for the expression levels of ACE2 CD147, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (D). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, TMPRSS=transmembrane serine protease, UMAP=uniform mani- fold approximation and projection. [33,34] Coexpression of entry receptor ACE2 and spike protein prim- coexpression patterns of viral entry genes. Although most ing enzymes in BBB pericytes enables SARS-CoV-2 entry into analyses are based on transcriptional levels, there is clear con- the cerebrospinal fluid and brain. SARS-CoV-2 may also infect sistency at the cellular level: the expression patterns of entry the telencephalon through the olfactory bulb. genes at the transcriptional level in a particular cell type can accurately predict the protein expression levels of entry recep- [35] tors. Our results were consistent with the coexpression Discussion patterns (ACE2 and TMPTSS2) previously identified among [33] scRNA-seq is a powerful technique that can be used to iden- susceptible cell types in multiple tissues. Furthermore, we tify susceptible cells through single-cell-level analyses of the obtained novel information concerning tissues such as the OE, 31 ReseaRch aR ticle Journal of Bio-X Research Figure 7. Single-cell analysis of healthy human lung cells. UMAP plot showing the distribution of lung cells (A). Dot plot showing the marker genes for lung clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. Violin and box plots for the expression levels of ACE2, CD147, CD209, CLEC4M, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (D). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. CTSL=cathepsin L, DCs=dendritic cells, NK=natural killer, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projection, VSMC=vascular smooth muscle cells. RE, and brain; we also identified coexpression patterns involv- The coronavirus infection process begins with viral bind- ing TMPRSS2 and viral replication genes. These results sug- ing to the cellular membrane through heparan sulfate pro- gest that neurological symptoms and CNS infection in patients teoglycans, which facilitate the interaction between the spike [16,36] with COVID-19 are secondary to SARS-CoV-2 invasion of the protein and the entry receptor (ie, ACE2). Recently, BBB via pericytes. Additionally, SARS-CoV-2–induced olfac- CD147 and CD209 were also identified as potential recep- tory disorders could be the result of localized damage to suste- tors for SARS-CoV-2; interactions with these proteins [17,18] ntacular cells and ciliated cells in the nasal RE and OE. also facilitate viral entry into host cells. In our study, 32 Journal of Bio-X Research ReseaRch aR ticle Figure 8. Expression levels of SARS-CoV-2 binding receptors, spike protein hydrolysis, endocytosis-related genes, and virus RNA replication genes across tissues. Violin and box plots for the expression levels of SARS-CoV-2 binding receptors, spike protein hydrolysis, endocytosis-related genes, and virus RNA replication genes across tissues (A). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. Box plots showing the expression levels of exosome biogenesis and release-related genes in the brain and olfactory epithelium (B). The expression levels of gene sets are shown as mean log2 (TP10K+1) values. Schematic depicting potential mechanisms involved in CNS infection (C). Coexpression of entry receptor ACE2 and spike protein priming enzyme cathepsin L in BBB pericytes enables SARS-CoV-2 entry into the cerebrospinal fluid and brain. SARS- CoV-2 may also infect the telencephalon through the olfactory bulb. ACE2=angiotensin-converting enzyme 2, BBB=blood–brain barrier, CNS=central nervous system, HBCs=horizontal basal cells, OEC=olfactory ensheathing glia cells, SARS-CoV-2=severe acute respiratory syndrome coronavirus 2. CD147 was expressed in most cell types; however, these in various cells. Thus, cells that include both the entry receptor cells were not consistent with published reports of infected (ACE2) and the spike protein priming enzyme (TMPRSS and/ [32] cells. Thus, the role of CD147 in SARS-CoV-2 pathoge- or CTSL) may be most susceptible to infection by SARS-CoV-2. nicity requires further investigation. In contrast, CD209 was Headache, impaired consciousness, epilepsy, encephalitis, mainly expressed in macrophages and monocytes in most and acute myelitis have been reported in patients with COVID- [7] tissues. The specific location of CD209 expression implies 19. Notably, a case of SARS-CoV-2 encephalitis was reported a role in SARS-CoV-2 entry into macrophages, which were by physicians in Beijing Ditan Hospital; genomic sequencing [32] also the target cells for SARS-CoV-2 infection, with the revealed the presence of SARS-CoV-2 in the patient’s cerebrospi- potential for initiation of a cytokine storm. nal fluid. Subsequently, a similar case was reported by physicians [10] Successful conformational changes and activation of spike in Japan. These reports indicated that SARS-CoV-2 is able to proteins result in membrane fusion; these processes require infect the human CNS in some instances. Intracranial infections [12] receptor binding, as well as appropriate protease activation. occur after damage to the BBB, which allows viral invasion. Endocytosis and cathepsin-mediated endosomal acidification In the present study, we found that the entry receptor ACE2 are also important processes during viral fusion with a target and spike protein priming enzyme CTSL were coexpressed in [37,38] [40] cell. In the present study, we found that TMPRSS and CTSL pericytes, which are specialized endothelial cells in the brain. were strongly expressed and exhibited a broad distribution. Furthermore, endothelial cells generally expressed lower levels Either TMPRSS or CTSL was expressed in ACE2-expressing of ACE2, which suggests that they can serve as an alternative cells, consistent with previous reports that ACE2 constitutes a route of infection. Notably, some vascular pericytes (eg, in [39] limiting factor for initial viral entry. Genes that encode pro- the OE and RE) did not express ACE2, indicating that brain teins involved in endocytosis and viral replication are expressed pericytes have a distinct role in blood pressure regulation and 33 ReseaRch aR ticle Journal of Bio-X Research [41] BBB maintenance. Endothelial cells in the brain are bound which were higher than the levels in AT2 cells. This finding together by tight junctions to ensure the integrity of the BBB. indicates that such ciliated cells may be the “portal” for SARS- The use of endothelial cells and pericytes as host cells could CoV-2 entry because of their anatomical location; they may also allow viral proliferation; the eventual cell death would destroy serve as a viral reservoir for persistent SARS-CoV-2 infection. the host cells and BBB, thereby releasing viral particles into the These results have important implications for understanding brain. Considering the extensive similarities described above in the symptoms of anosmia and potential mechanisms involved in terms of expression patterns for ACE2, TMPRSS2, and CTSL SARS-CoV-2 infection of the CNS. in mice and humans, neurons presumably do not constitute a primary site of infection; however, vascular pericytes may be Limitations sensitive to SARS-CoV-2, thus allowing entry through the BBB and subsequent CNS infection. In this study, the coexpression pattern of virus entry genes only Respiratory viruses may enter the CNS through neurons in indicates the susceptibility of the cell types. It does not really the olfactory and trigeminal nerves (within the nasal cavity), represent the cell types the SARS-CoV-2 infects. Furthermore, or through neurons in the vagus nerve (within the trachea and we analyzed susceptible cell types based on only normal sam- [42] lungs). Human coronavirus OC43 can induce CNS neu- ples. It would be more convincing if we used samples from ropathology through interneuron propagation and axonal COVID-19 patients to study the viral load in different cell types. [43] transport. We found that the entry receptor for human coro- So, the results were limited. navirus OC43 (HLA class I molecule) was expressed in neural cells, which enables the virus to infect neurons (Additional Fig. Conclusion 2, http://links.lww.com/JR9/A45). However, the entry receptor ACE2 and the spike protein priming enzymes were not coex- We identified cells susceptible to SARS-CoV-2 infection in the pressed in neural cells within the CNS, olfactory system, or CNS, olfactory system, and respiratory system. Our findings respiratory system. An earlier epidemic strain of coronavirus, provide insights into the mechanisms that underlie neurological SARS-CoV, also binds ACE2; SARS-CoV particles have been symptoms and CNS infection in patients with COVID-19. detected in brain tissue, where they were located almost exclu- [44] sively in neurons. Considering the strong capacities for exo- Acknowledgments some biogenesis and release in neurons, ciliated cells, and OECs, None. exosome biogenesis and release could be involved in SARS- CoV-2 infection of the CNS. However, this hypothesis requires additional investigation. Author contributions The OE serves as the main odor detection and conduction HX, JL: study design. ZK, JW, HZ, DX, HG, and WZ: data tissue; it includes ciliated olfactory sensory neurons, glia-like [45] analysis. ZK, DX, and HG: data collection and generation. HZ, sustentacular cells, and apical microvilli. Olfactory sensory JW, ZL, XC, and JX: data interpretation. TM, ZK, AQ, and HX: neurons are surrounded by sustentacular cells and the axons of manuscript drafting. AQ, HZ, TM, HX, and JL: overall super- OECs. Sustentacular cells serve as structural support for sensory vising and organizing the study. All authors approved the final neurons, whereas OECs guide olfactory sensory neuron axons version for publication. that cross the cribriform plate at the skull base and terminate in [46–48] the olfactory bulb. Additionally, olfactory sensory neurons [49] express odor receptors on dendritic cilia. Thus, SARS-CoV-2 Financial support invasion into ciliated cells and sustentacular cells may lead to This work was supported by the National Natural Ascience olfactory disorders. Foundation of China (No. 31821003 to HX) and the China SARS-CoV-2 primarily affects the respiratory system, with Ministry of Science and Technology (No. 2018AAA0100300 to [50] serious effects on the lungs. In this study, we observed coex- HX). The funders had no role in manuscript design, data collec- pression of ACE2 and TMPRSS/CTSL in the AT1, AT2, and cil- tion and analysis, preparation of the manuscript, or decision to iated cells in lung tissue. In lung alveoli, AT1 cells (ie, epithelial publish. cells) are responsible for gas exchange; AT2 cells (ie, alveolar stem-like cells) are responsible for surfactant biosynthesis and [51] Institutional review board statement immunoregulation. Thus, SARS-CoV-2 infection may cause damage to alveolar infrastructure; the resulting increase in There is no direct involvement of human subjects in this project. alveolar surface tension can induce dyspnea and hypoxia. The All the data use existing biological samples and data from prior [52] involvement of AT2 cells can lead to extensive inflammation. studies. Therefore, ethical oversight and patient consent were The host immune response produces a proinflammatory cyto- not handled in this study. kine storm, which manifests as acute lung injury and acute [53,54] respiratory distress syndrome. We found that CD209 was Data availability statement mainly expressed in macrophages and monocytes in lung tissue. Notably, CD209 may be responsible for viral antigen activation The single cell data are available at the Gene Expression Omnibus of specific T cells, as well as the secretion of perforin, interferon, under the accession code GSE139522 and GSE134355, Sequence and granzyme B; these processes can cause further injury to the ReadArchive dataset (SRA https://www.ncbi.nlm.nih.gov/sra) at [55] lungs and remote organs, such as the brain. accession code SRA667466 and Atlas Data Protal (https://data. Finally, we compared the expression levels of viral entry fac- humancellatlas.org) at accession code PRJEB31843. Brochus tors across olfactory and respiratory systems; we found gener- data will be made available Mendeley Data (https://data.men- ally high expression levels in ciliated cells in the nasal mucosa, deley.com/datasets/7r2cwbw44m/1). 34 Journal of Bio-X Research ReseaRch aR ticle [23] Zeisel A, Muñoz-Manchado AB, Codeluppi S, et al. Brain structure. Conflicts of interest Cell types in the mouse cortex and hippocampus revealed by single-cell There are no conflicts of interest. RNA-seq. Science 2015;347:1138–1142. doi:10.1126/science.aaa1934. [24] Madissoon E, Wilbrey-Clark A, Miragaia RJ, et al. scRNA-seq assessment of the human lung, spleen, and esophagus tissue stabil- References ity after cold preservation. Genome Biol 2019;21:1. doi:10.1186/ [1] Guan WJ, Ni ZY, Hu Y, et al. Clinical characteristics of coronavirus dis- s13059-019-1906-x. ease 2019 in China. N Engl J Med 2020;382:1708–1720. doi:10.1056/ [25] Lukassen S, Chua RL, Trefzer T, et al. SARS-CoV-2 receptor ACE2 NEJMoa2002032. and TMPRSS2 are primarily expressed in bronchial transient secretory [2] Chen N, Zhou M, Dong X, et al. Epidemiological and clinical charac- cells. EMBO J 2020;39:e105114. doi:10.15252/embj.20105114. teristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, [26] Luecken MD, Theis FJ. Current best practices in single-cell RNA- China: a descriptive study. Lancet 2020;395:507–513. doi:10.1016/ seq analysis: a tutorial. Mol Syst Biol 2019;15:e8746. doi:10.15252/ S0140-6736(20)30211-7. msb.20188746. [3] Huang C, Wang Y, Li X, et al. Clinical features of patients infected with [27] Stuart T, Butler A, Hoffman P, et al. Comprehensive integration 2019 novel coronavirus in Wuhan, China. Lancet 2020;395:497–506. of single-cell data. Cell 2019;177:1888–1902.e21. doi:10.1016/j. doi:10.1016/S0140-6736(20)30183-5. cell.2019.05.031. [4] Uginet M, Breville G, Assal F, et al. COVID-19 encephalopathy: clin- [28] Korsunsky I, Millard N, Fan J, et al. Fast, sensitive and accurate integra- ical and neurobiological features. J Med Virol 2021;93:4374–4381. tion of single-cell data with Harmony. Nat Methods 2019;16:1289–1296. doi:10.1002/jmv.26973. doi:10.1038/s41592-019-0619-0. [5] Sampaio Rocha-Filho PA, Albuquerque PM, Carvalho LCLS, et al. [29] Zhao Y, Wong L, Goh WWB. How to do quantile normalization Headache, anosmia, ageusia and other neurological symptoms in correctly for gene expression data analyses. Sci Rep 2020;10:15534. COVID-19: a cross-sectional study. J Headache Pain 2022;23:1–11. doi:10.1038/s41598-020-72664-6. doi:10.1186/s10194-021-01367-8. [30] Townes FW, Irizarry RA. Quantile normalization of single-cell RNA- [6] Dell’Era V, Farri F, Garzaro G, et al. Smell and taste disorders during seq read counts without unique molecular identifiers. Genome Biol COVID-19 outbreak: cross-sectional study on 355 patients. Head Neck 2020;21:160. doi:10.1186/s13059-020-02078-0. 2020;42:1591–1596. doi:10.1002/hed.26288. [31] Lukenda A, Dotlic S, Vukojevic N, et al. Expression and prognostic [7] Wu Y, Xu X, Chen Z, et al. Nervous system involvement after infec- value of putative cancer stem cell markers CD117 and CD15 in cho- tion with COVID-19 and other coronaviruses. Brain Behav Immun roidal and ciliary body melanoma. J Clin Pathol 2016;69:234–239. 2020;87:18–22. doi:10.1016/j.bbi.2020.03.031. doi:10.1136/jclinpath-2015-203130. [8] Helms J, Kremer S, Merdji H, et al. Neurologic features in severe SARS- [32] Bost P, Giladi A, Liu Y, et al. Host-viral infection maps reveal signa- CoV-2 Infection. N Engl J Med 2020;382:2268–2270. doi:10.1056/ tures of severe COVID-19 patients. Cell 2020;181:1475–1488.e12. NEJMc2008597. doi:10.1016/j.cell.2020.05.006. [9] Xu Z, Shi L, Wang Y, et al. Pathological findings of COVID-19 asso- [33] Sungnak W, Huang N, Bécavin C, et al. SARS-CoV-2 entry factors are ciated with acute respiratory distress syndrome. Lancet Respir Med highly expressed in nasal epithelial cells together with innate immune 2020;8:420–422. doi:10.1016/S2213-2600(20)30076-X. genes. Nat Med 2020;26:681–687. doi:10.1038/s41591-020-0868-6. [10] Moriguchi T, Harii N, Goto J, et al. A first case of meningitis/encepha- [34] Zhang H, Kang Z, Gong H, et al. Digestive system is a potential route litis associated with SARS-Coronavirus-2. Int J Infect Dis 2020;94:55– of COVID-19: an analysis of single-cell coexpression pattern of key 58. doi:10.1016/j.ijid.2020.03.062. proteins in viral entry process. Gut 2020;69:1010–1018. doi:10.1136/ [11] Shehata GA, Lord KC, Grudzinski MC, et al. Neurological complica- gutjnl-2020-320953. tions of covid-19: Underlying mechanisms and management. Int J Mol [35] Brann DH, Tsukahara T, Weinreb C, et al. Non-neuronal expression of Sci 2021;22:4081. doi:10.3390/ijms22084081. SARS-CoV-2 entry genes in the olfactory system suggests mechanisms [12] Li Z, Liu T, Yang N, et al. Neurological manifestations of patients with underlying COVID-19-associated anosmia. Sci Adv 2020;6:eabc5801. COVID-19: potential routes of SARS-CoV-2 neuroinvasion from the doi:10.1126/sciadv.abc5801. periphery to the brain. Front Med 2020;14:533–541. doi:10.1007/ [36] Wan Y, Shang J, Graham R, et al. Receptor recognition by the novel s11684-020-0786-5. coronavirus from Wuhan: an analysis based on decade-long struc- [13] Shang J, Wan Y, Luo C, et al. Cell entry mechanisms of SARS-CoV-2. tural studies of SARS coronavirus. J Virol 2020;94:e00127–e00120. Proc Natl Acad Sci U S A 2020;117:11727–11734. doi:10.1073/ doi:10.1128/JVI.00127-20. pnas.2003138117. [37] Knyazev E, Nersisyan S, Tonevitsky A. Endocytosis and transcytosis of [14] Chen M, Shen W, Rowan NR, et al. Elevated ACE-2 expression in the SARS-CoV-2 across the intestinal epithelium and other tissue barriers. olfactory neuroepithelium: implications for anosmia and upper respira- Front Immunol 2021;12:636966. doi:10.3389/fimmu.2021.636966. tory SARS-CoV-2 entry and replication. Eur Respir J 2020;56:2001948. [38] Zhao MM, Yang WL, Yang FY, et al. Cathepsin L plays a key role in doi:10.1183/13993003.01948-2020. SARS-CoV-2 infection in humans and humanized mice and is a prom- [15] Zhou P, Yang XL, Wang XG, et al. A pneumonia outbreak associated with ising target for new drug development. Signal Transduct Target Ther a new coronavirus of probable bat origin. Nature 2020;579:270–273. 2021;6:134. doi:10.1038/s41392-021-00558-8. doi:10.1038/s41586-020-2012-7. [39] Ou X, Liu Y, Lei X, et al. Characterization of spike glycoprotein of [16] Wrapp D, Wang N, Corbett KS, et al. Cryo-EM structure of the 2019- SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS- nCoV spike in the prefusion conformation. Science 2020;367:1260–1263. CoV. Nat Commun 2020;11:1620. doi:10.1038/s41467-020-15562-9. doi:10.1126/science.abb2507. [40] Campisi M, Shin Y, Osaki T, et al. 3D self-organized microvascular [17] Wang K, Chen W, Zhang Z, et al. CD147-spike protein is a novel route model of the human blood-brain barrier with endothelial cells, peri- for SARS-CoV-2 infection to host cells. Signal Transduct Target Ther cytes and astrocytes. Biomaterials 2018;180:117–129. doi:10.1016/j. 2020;5:283. doi:10.1038/s41392-020-00426-x. biomaterials.2018.07.014. [18] Amraei R, Yin W, Napoleon MA, et al. CD209L/L-SIGN and CD209/ [41] Hirunpattarasilp C, Attwell D, Freitas F. The role of pericytes in brain DC-SIGN act as receptors for SARS-CoV-2. ACS Cent Sci 2021;7:1156– disorders: from the periphery to the brain. J Neurochem 2019;150:648– 1165. doi:10.1021/acscentsci.0c01537. 665. doi:10.1111/jnc.14725. [19] Hoffmann M, Kleine-Weber H, Schroeder S, et al. SARS-CoV-2 cell [42] Mori I. Transolfactory neuroinvasion by viruses threatens the human entry depends on ACE2 and TMPRSS2 and is blocked by a clinically brain. Acta Virol 2015;59:338–349. doi:10.4149/av_2015_04_338. proven protease inhibitor. Cell 2020;181:271–280.e8. doi:10.1016/j. [43] Dubé M, Le Coupanec A, Wong AHM, et al. Axonal transport enables cell.2020.02.052. neuron-to-neuron propagation of human coronavirus OC43. J Virol [20] Meng T, Cao H, Zhang H, et al. The insert sequence in SARS- 2018;92:e00404–e00418. doi:10.1128/JVI.00404-18. CoV-2 enhances spike protein cleavage by TMPRSS. bioRxiv 2020. [44] Xu J, Zhong S, Liu J, et al. Detection of severe acute respiratory syndrome doi:10.1101/2020.02.08.926006. coronavirus in the brain: Potential role of the chemokine mig in patho- [21] Durante MA, Kurtenbach S, Sargi ZB, et al. Single-cell analysis of olfac- genesis. Clin Infect Dis 2005;41:1089–1096. doi:10.1086/444461. tory neurogenesis and differentiation in adult humans. Nat Neurosci [45] Villar PS, Delgado R, Vergara C, et al. Energy requirements of odor 2020;23:323–326. doi:10.1038/s41593-020-0587-9. transduction in the chemosensory cilia of olfactory sensory neu- [22] Han X, Zhou Z, Fei L, et al. Construction of a human cell land- rons rely on oxidative phosphorylation and glycolytic processing of scape at single-cell level. Nature 2020;581:303–309. doi:10.1038/ extracellular glucose. J Neurosci 2017;37:5736–5743. doi:10.1523/ s41586-020-2157-4. JNEUROSCI.2640-16.2017. 35 ReseaRch aR ticle Journal of Bio-X Research [46] Kondoh K, Lu Z, Ye X, et al. A specific area of olfactory cortex involved [51] Nabhan AN, Brownfield DG, Harbury PB, et al. Single-cell Wnt sig- in stress hormone responses to predator odours. Nature 2016;532:103– naling niches maintain stemness of alveolar type 2 cells. Science 106. doi:10.1038/nature17156. 2018;359:1118–1123. doi:10.1126/science.aam6603. [47] Le Bourhis M, Rimbaud S, Grebert D, et al. Endothelin uncouples gap [52] Kroetz DN, Allen RM, Schaller MA, et al. Type I interferon induced epi- junctions in sustentacular cells and olfactory ensheathing cells of the genetic regulation of macrophages suppresses innate and adaptive immu- olfactory mucosa. Eur J Neurosci 2014;40:2878–2887. doi:10.1111/ nity in acute respiratory viral infection. PLoS Pathog 2015;11:e1005338. ejn.12665. doi:10.1371/journal.ppat.1005338. [48] Beites CL, Kawauchi S, Crocker CE, et al. Identification and molecular [53] Desforges M, Le Coupanec A, Dubeau P, et al. Human coronaviruses and regulation of neural stem cells in the olfactory epithelium. Exp Cell Res other respiratory viruses: underestimated opportunistic pathogens of the 2005;306:309–316. doi:10.1016/j.yexcr.2005.03.027. central nervous system? Viruses 2019;12:14. doi:10.3390/v12010014. [49] Savya SP, Kunkhyen T, Cheetham CEJ. Low survival rate of young [54] Mehta P, McAuley DF, Brown M, et al. COVID-19: consider cytokine adult-born olfactory sensory neurons in the undamaged mouse olfac- storm syndromes and immunosuppression. Lancet 2020;395:1033– tory epithelium. J Bioenerg Biomembr 2019;51:41–51. doi:10.1007/ 1034. doi:10.1016/S0140-6736(20)30628-0. s10863-018-9774-8. [55] Channappanavar R, Perlman S. Pathogenic human coronavirus [50] Guan WJ, Liang WH, Zhao Y, et al. Comorbidity and its impact on infections: causes and consequences of cytokine storm and immuno- 1590 patients with COVID-19 in China: a nationwide analysis. Eur pathology. Semin Immunopathol 2017;39:529–539. doi:10.1007/ Respir J 2020;55:2000547. doi:10.1183/13993003.00547-2020. s00281-017-0629-x. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Journal of Bio-X Research Wolters Kluwer Health

Loading next page...
 
/lp/wolters-kluwer-health/pericytes-may-facilitate-sars-cov-2-entry-into-the-nervous-system-mUCHdOINIU

References

References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.

Publisher
Wolters Kluwer Health
Copyright
Copyright © 2023 The Chinese Medical Association, Published by Wolters Kluwer Health, Inc.
ISSN
2577-3585
eISSN
2096-5672
DOI
10.1097/jbr.0000000000000137
Publisher site
See Article on Publisher Site

Abstract

Objective: Although the neurological and olfactory symptoms of coronavirus disease 2019 have been identified, the neurotropic properties of the causative virus, severe acute respiratory syndrome-associated coronavirus 2 (SARS-CoV-2), remain unknown. We sought to identify the susceptible cell types and potential routes of SARS-CoV-2 entry into the central nervous system, olfactory system, and respiratory system. Methods: We collected single-cell RNA data from normal brain and nasal epithelium specimens, along with bronchial, tracheal, and lung specimens in public datasets. The susceptible cell types that express SARS-CoV-2 entry genes were identified using single-cell RNA sequencing and the expression of the key genes at protein levels was verified by immunohistochemistry. We compared the coexpression patterns of the entry receptor angiotensin-converting enzyme 2 (ACE2) and the spike protein priming enzyme transmembrane serine protease (TMPRSS)/cathepsin L among the specimens. Results: The SARS-CoV-2 entry receptor ACE2 and the spike protein priming enzyme TMPRSS/cathepsin L were coexpressed by pericytes in brain tissue; this coexpression was confirmed by immunohistochemistry. In the nasal epithelium, ciliated cells and sustentacular cells exhibited strong coexpression of ACE2 and TMPRSS. Neurons and glia in the brain and nasal epithelium did not exhibit coexpression of ACE2 and TMPRSS. However, coexpression was present in ciliated cells, vascular smooth muscle cells, and fibroblasts in tracheal tissue; ciliated cells and goblet cells in bronchial tissue; and alveolar epithelium type 1 cells, AT2 cells, and ciliated cells in lung tissue. Conclusion: Neurological symptoms in patients with coronavirus disease 2019 could be associated with SARS-CoV-2 invasion across the blood–brain barrier via pericytes. Additionally, SARS-CoV-2–induced olfactory disorders could be the result of localized cell damage in the nasal epithelium. Keywords: ACE2, coronavirus disease 2019, SARS-CoV-2, single-cell RNA sequencing, TMPRSS2 ZK, JW, and TM contributed equally to the article. Introduction Department of Rheumatology and Immunology, Changzheng Hospital, The main clinical manifestations of coronavirus disease 2019 Department of Neurosurgery, Changhai Hospital, Second Military Medical (COVID-19) resulting from severe acute respiratory syndrome University, Division of Spine, Department of Orthopedics, Shanghai General coronavirus 2 (SARS-CoV-2) infection are fever, cough, and dys- Hospital, School of Medicine, Shanghai Jiaotong University, Tongji University [1–5] e pnea. However, numerous neurological abnormalities have Cancer Center, School of Medicine, Tongji University, Department of Orthopaedic been described in patients with COVID-19; these abnormalities Oncology, Changzheng Hospital, Department of Urology, The Third Affiliated Hospital, Second Military Medical University, Department of Pathology, Shanghai include olfactory and gustatory disorders (anosmia, hyposmia, Tenth People’s Hospital, Tongji University School of Medicine, Shanghai, China, ageusia, and dysgeusia), headache, epilepsy, encephalitis, and Zeenat Qureshi Stroke Institutes and University of Missouri, Columbia, MO, [4–8] acute myelitis. SARS-CoV-2 RNA has also been identified in i j USA, Peking-Tsinghua Center for Life Sciences, School of Medicine, Tsinghua [9,10] cerebrospinal fluid via genome sequencing. These findings University, Beijing, China suggest that SARS-CoV-2 can directly or indirectly enter the *Corresponding authors: Huji Xu, Department of Rheumatology and Immunology, central nervous system (CNS). However, the underlying mecha- Changzheng Hospital, Second Military Medical University, Shanghai 200003, nism remains unknown. China. E-mail: xuhuji@smmu.edu.cn; Jianmin Liu, Department of Neurosurgery, Changhai Hospital, Second Military Medical University, Shanghai 200003, China. The CNS can be invaded through the hematogenous route or [11] E-mail: chstroke@163.com by axonal transport ; these routes involve infections of endo- thelial cells in the blood–brain barrier (BBB) or epithelial cells Copyright © 2023 The Chinese Medical Association, Published by Wolters Kluwer Health, Inc. This is an open-access article distributed under the terms of in the blood–cerebrospinal fluid barrier, as well as infections [12] the Creative Commons Attribution-Non Commercial-No Derivatives License 4.0 of the olfactory bulb through the olfactory nerve. The tro- (CCBY-NC-ND), where it is permissible to download and share the work provided pism of coronaviruses is primarily dependent on the presence of it is properly cited. The work cannot be changed in any way or used commercially attachment and entry receptors, along with protease availability without permission from the journal. and the activities of internalization-related and trafficking-re- Journal of Bio-X Research (2023) 6:23–36 [13,14] lated pathways. SARS-CoV-2 uses angiotensin-convert- [15,16] Received: 4 May 2022; Accepted: 7 December 2022 ing enzyme 2 (ACE2) as an entry receptor. Other entry receptor molecules associated with SARS-CoV-2 may include http://dx.doi.org/10.1097/JBR.0000000000000137 23 ReseaRch aR ticle Journal of Bio-X Research [17,18] CD147 and CD209. Upon receptor binding, spike protein Cell type identification and gene expression analysis priming enzymes (eg, cell surface transmembrane serine pro- Clusters were annotated based on the expression patterns of tease [TMPRSS] and endosomal cathepsins) are necessary for known cellular markers; the clustering information was visu- [19,20] some coronaviruses to achieve membrane fusion. Thus, the alized by uniform manifold approximation and projection detection of coexpression involving key genes in the viral entry (UMAP) plots. Then, lost values in the gene expression matrix process may help to identify cells that are susceptible to SARS- were imputed by the “RunALRA” function. The imputed gene CoV-2 infection; this will contribute to further investigations of expression values are displayed in feature and violin plots. the mechanisms by which SARS-CoV-2 enters the CNS. The mean expression levels of gene sets are displayed in box In this study, we used single-cell RNA sequencing (scRNA- plots. For comparisons of expression levels among datasets, seq) of normal brain and nasal epithelium specimens, as well as the quantile method was used to correct expression data; this tracheal, bronchial, and lung specimens, to identify susceptible [29,30] global normalization method is widely used. To calculate cell types and potential routes of SARS-CoV-2 entry into the the endocytosis signaling score, we obtained endocytosis-as- CNS, olfactory system, and respiratory system. We sought to sociated genes from “GOBP_ENDOCYTOSIS” gen set in the clarify how SARS-CoV-2 invades various tissues in the human MigDB dataset (https://www.gsea-msigdb.org). The average body, with a particular interest in elucidating the olfactory and expression level of all the genes in this gene set was considered CNS symptoms in patients with COVID-19. as the expression level of the endocytosis signaling. Materials and methods Immunohistochemical staining Data sources Human brain tissue specimens were obtained from five patients We analyzed single-cell transcriptomes of the brain, nasal olfac- who had undergone surgery for treatment of traumatic brain tory epithelium (OE), nasal respiratory epithelium (RE), lung, injury, including three males and two females. Human lung tis- tracheal epithelium, and bronchial epithelium (Additional Table sue specimens were obtained during routine surgical intervention 1, http://links.lww.com/JR9/A46). Nasal RE/OE and trachea from three lung cancer patients, including two males and one transcriptome data were obtained from the Gene Expression female. Ethical approval was granted by the Ethics Committee Omnibus (GEO; https://www.ncbi.nlm.nih.gov/), using acces- of Shanghai Changhai Hospital (Approval No. CHEC2020- [21] [22] sion codes GSE139522 and GSE134355. Brain transcrip- 164). All patients were fully informed of the research and signed tome data were obtained from the Sequence Read Archive the informed consent forms. Formalin-fixed paraffin-embedded (SRA; https://www.ncbi.nlm.nih.gov/sra), using accession code brain tissue sections on glass slides were incubated at 37°C over- [23] SRA667466. Lung transcriptome data were obtained from night, deparaffinized in xylene, and rehydrated in a descending Cell Atlas Data Portal (https://data.humancellatlas.org), using ethanol gradient. Then, tissue sections were incubated in antigen [24] accession code PRJEB31843. Bronchial transcriptome data retrieval solution at 95°C for 30 minutes. Next, tissue sections [25] were obtained from research data published by Lukassen et al. were incubated in 3% hydrogen peroxide for 10 minutes, fol- The single-cell data within each dataset had similar sequencing lowed by serum-free protein blocking solution (Cat# 9048-46- depth to ensure the quality of consensus sequence. Because the 8, Sigma, St. Louis, MO, USA) for 20 minutes; they were then data were obtained from public datasets, approval for our study incubated with rabbit anti-human ACE2 primary antibody (Cat# by the ethics committee and patient consent were unnecessary. ab108252, Abcam, Cambridge, UK) at a dilution ratio of 1:200 overnight at 4°C temperature and goat anti-rabbit IgG second- ary antibody (Cat# ab6721, Abcam) at a dilution ratio of 1:1000 Quality control for 1 hour at room temperature. After incubation with primary Barcodes with outlier gene counts may be indicative of dying and secondary antibodies, color development was performed cells, cells with broken membranes, or doublets. A high frac- using 3-amino-9-ethylcarbazole as the chromogenic substrate tion of mitochondria is indicative of leaky cytoplasmic mem- and tissue was counterstained with hematoxylin as described [26] [31] branes allowing mRNA to escape. We removed low-quality previously. The microscopic observation was performed by data from cells with <200 or >5000 expressed genes, as well as the optical microscope (Eclipse 80i, Nikon, Tokyo, Japan). data from cells in which >10% of unique molecular identifiers mapped to mitochondrial genes. Results Single-cell maps of mouse brain, healthy human Data integration, dimension reduction, and cell olfactory system, and healthy human respiratory system clustering [27] [28] To identify susceptible cell types in the brain, olfactory system, The R package Seurat and R package Harmony were and respiratory system, we collected scRNA-seq data for six used to process scRNA-seq data. Unique molecular identifier types of tissue from published papers (Additional Table 1, http:// counts were normalized using the “NormalizeData” function. The top 1000 highly variable genes were identified using the links.lww.com/JR9/A46). In the brain, 11 cell types were identi- “FindVariableGenes” function. Principal component analy- fied based on the analysis of 29,678 cells, after the application sis was performed on single-cell expression matrices using the of quality control criteria (Fig. 1A). The main cell types were oli- “RunPCA” function. Subsequently, “Runharmony” was used for godendrocytes, interneurons, astrocytes, microglia, and neurons batch-effect correction in each dataset. Cell clustering analysis 1 and 2. The other identified cell types included ependymal cells, was performed after the construction of a K-nearest-neighbor endothelial cells, pericytes, dividing cells, and myeloid cells. graph based on the adjusted principal component analysis The nasal epithelium was divided into two types: OE and RE. embeddings by harmony algorithm. In total, 19,903 cells and 14 cell types were identified in the OE 24 Journal of Bio-X Research ReseaRch aR ticle Figure 1. Single-cell maps of healthy human respiratory system and mouse brain tissue. UMAP plots showing the distributions of brain cells (A), nasal olfactory epithelial cells (B), nasal respiratory epithelial cells (C), tracheal epithelial cells (D), bronchial epithelial cells (E), and lung cells (F). Single-cell map of respiratory system (G). All the single-cell data in the human respiratory system were combined for another round of reduction and clustering analysis. UMAP plots show the distributions of the cell clusters in the human respiratory system. DCs=dendritic cells, HBCs=horizontal basal cells, NK=natural killer, UMAP=uniform manifold approximation and projection, VSMC=vascular smooth muscle cells. (Fig. 1B); 10,657 cells and 12 cell types were identified in the RE as well as chondrocytes, VSMCs, endothelial cells, fibroblasts, (Fig. 1C). Epithelial cell types in the RE were horizontal basal and immune cells (ie, myeloid and B cells). Additionally, 17,228 cells (HBCs) and Bowman’s gland cells, as well as sustentac- cells and seven cell types were identified in bronchial specimens ular, ciliated, and endodermal progenitor cells. Non-epithelial after the application of quality control criteria (Fig.  1E). The cell types in the RE were microvillar and endothelial cells, as main cell types identified were ciliated, basal, and goblet cells. well as fibroblasts and immune cells (ie, mast, myeloid, T, and B The other identified cell types consisted of brush cells, fibro- cells). Epithelial cell types in the OE were HBCs, pericytes, and blasts, dividing cells, and ionocytes. Bowman’s gland cells, as well as gland progenitor, sustentacular, After the application of quality control criteria, 24,710 cells and ciliated cells. Non-epithelial cell types in the OE were neu- and 14 cell types were identified in lung specimens (Fig.  1F). rons, olfactory ensheathing glia cells (OECs), fibroblasts, vascu- The cell types identified were alveolar epithelium type 1 (AT1), lar smooth muscle cells (VSMCs), and immune cells (ie, mast, alveolar epithelium type 2 (AT2), and ciliated cells, as well as myeloid, T, and B cells). fibroblasts, lymph vessel cells, VSMCs, and endothelial cells. After the application of quality control criteria, 9614 cells and The following immune cell types were identified: mast, T, B, and 10 cell types were identified in tracheal specimens (Fig. 1D). The natural killer (NK) cells, along with monocytes, macrophages, cell types were goblet, brush, ciliated, basal, and granular cells, and dendritic cells. 25 ReseaRch aR ticle Journal of Bio-X Research Because of the similar cells present among tissues in the respi- progenitor, ciliated, and sustentacular cells, as well as neu- ratory system, we constructed a map of the respiratory system rons; its expression was minimal in fibroblasts, HBCs, and based on olfactory epithelum data, nasal respiratory epithelium Bowman’s gland cells (Fig.  3C and D). CTSL was expressed data, trachea data, bronchus data, and lung data from 83,126 in all non-immune subclusters, but it was only expressed in cells. In this map, we identified 17 subclusters: HBCs, endo- myeloid cells among the immune subclusters (Fig. 3C and D). thelial cells, Bowman’s gland cells, sustentacular cells, alveolar Additionally, endocytosis signaling was strongly expressed cells, ciliated cells, neurons, OECs, fibroblasts, microvillar cells, in non-immune subclusters and weakly expressed in immune VSMCs, chondrocytes, basal cells, mast cells, myeloid cells, B subclusters (Fig.  3C and D). PIKFYVE and TPCN2 were cells, and T/NK cells (Fig. 1G). expressed in most cell types (Fig.  3D). The viral RNA repli- cation target gene MTHFD1 was detected in all subclusters in the OE, except T and B cells (Fig.  3D). Overall, the data Coexpression patterns of key entry genes in mouse brain indicated that multiple cell types in the OE could serve as tar- tissue get cells for SARS-CoV-2. The gland progenitor, ciliated, and To identify susceptible cell types in brain tissue and elucidate sustentacular cells may constitute the most susceptible cells potential routes of SARS-CoV-2 entry, we collected scRNA-seq in the OE because they exhibited the strongest expression of data from mouse brains. In total, 11 subclusters were identified ACE2 and TMPRSS2 (Fig. 3C and D). Thus, we speculate that in mouse brain cells (Fig. 2A); marker genes for each subclus- non-neuronal cell types, rather than sensory or bulb neurons, ter are shown in Figure  2B. ACE2 was strongly expressed in are responsible for anosmia and associated odor perception pericytes and weakly expressed in other subclusters (Fig. 2C). deficits in patients with COVID-19. To confirm the expression of ACE2 at the protein level, we performed immunohistochemical analysis of normal human Coexpression patterns of key entry genes in healthy brain tissues; the immunohistochemistry results were consistent human nasal respiratory epithelial tissue with the scRNA-seq findings. ACE2 was strongly expressed in vascular pericytes, as well as endothelial cells (Fig. 2D). Lung RE is a major component of the nasal epithelium. We identified alveolar cells also exhibited strong expression of ACE2 at the 12 subclusters (Fig.  4A) and the marker genes for each sub- protein level (Fig. 2E); they are regarded as the main target of cluster are shown in Figure 4B. ACE2 was strongly expressed [32,33] SARS-CoV-2 lung tissue. Therefore, pericytes may be the in sustentacular, ciliated, and microvillar cells; it was weakly main target of SARS-CoV-2 in brain tissue; entry through these expressed in endodermal progenitor cells (Fig.  4C and D). cells might be involved in indirect damage to neurons. CD147, CD147 was expressed in almost all subclusters, whereas a potential entry receptor, was expressed in almost all sub- CD209 was expressed in myeloid cells (Fig.  4D). TMPRSS clusters (Fig. 2F). The spike protein priming enzyme TMPRSS and CTSL had similar expression patterns among subclusters was nearly absent from all subclusters in brain tissue, whereas in the RE; they did not exhibit robust expression in immune cathepsin L (CTSL) was expressed in all subclusters (Fig.  2C cells (Fig. 4C and D). Thus, multiple cell types (eg, microvillar, and F); these findings suggested that spike protein hydrolysis ciliated, and sustentacular) expressed ACE2, TMPRSS/CTSL, in brain tissue is mediated by CTSL. Because endocytosis is an endocytosis signaling, PIKFYVE, TPCN2, and MTHFD1; important aspect of coronavirus infection, we examined the these cells were identified as potential targets of SARS-CoV-2 expression of the endocytosis signaling molecules phospho- (Fig. 4D). inositide kinase, FYVE-type zinc finger containing (PIKFYVE) and two pore segment channel 2 (TPCN2). We found that Coexpression patterns of key entry genes in healthy PIKFYVE was strongly expressed in ependymal cells, peri- human tracheal and bronchial tissues cytes, and neurons 1 and 2 in brain tissue (Fig.  2F), whereas the expression of TPCN2 was minimal and primarily observed Among the 11 cell types identified in tracheal tissue (Fig. 5A in myeloid cells. Furthermore, methylenetetrahydrofolate dehy- and B), ACE2 was expressed in fibroblasts, VSMCs, and cili- drogenase 1 (MTHFD1), which is an essential gene for RNA ated cells (Fig.  5C and D). CD147 was expressed in all sub- [32] virus replication, was expressed in all subclusters in brain clusters, whereas CD209 was only expressed in myeloid cells tissue (Fig.  2F). These results suggested that neurons do not (Fig.  5D). TMPRSS was mainly expressed in goblet, granule, constitute a primary site of infection; however, pericytes are the brush, basal, and ciliated cells; in contrast, CTSL was expressed most susceptible cells in brain tissue because they exhibit strong in chondrocytes, endothelial cells, VSMCs, fibroblasts, myeloid expression of ACE2 and coexpression of CTSL, which facilitate cells, and B cells (Fig.  5C and D). Endocytosis signaling, viral entry and replication. PIKFYVE, TPCN2, and MTHFD1 were strongly expressed in all cell types in tracheal tissue (Fig. 5D). Among the seven sub- clusters in bronchial tissue (Fig. 6A and B), ACE2 was mainly Coexpression patterns of key entry genes in healthy expressed in goblet and ciliated cells, whereas CD147 was human nasal olfactory epithelial tissue expressed in almost all subclusters (Fig. 6C). TMPRSS, CTSL, Next, we explored whether SARS-CoV-2 could enter the CNS endocytosis signaling, PIKFYVE, TPCN2, and MTHFD1 were through the olfactory system. We identified 14 subclusters in expressed in all cell types in bronchial tissue (Fig. 6C and D). the OE (Fig.  3A); the marker genes for each subcluster are Taken altogether, the results indicated that ciliated cells are shown in Figure  3B. ACE2 was strongly expressed in suste- the most susceptible cells in both tracheal and bronchial tissue ntacular, ciliated, and gland progenitor cells; its expression because they generally exhibit strong expression of ACE2 and was weaker in HBCs and Bowman’s gland cells (Fig. 3C and TMPRSS (Figs.  5D and 6C). Goblet cells also demonstrated D). CD147 was expressed in almost all subclusters, whereas susceptibility to SARS-CoV-2 through its coexpression of CD209 was only expressed in myeloid cells (Fig. 3D). TMPRSS ACE2, TMPRSS/CTSL, endocytosis signaling, and MTHFD1 was expressed in many subclusters in the OE, including gland (Figs. 5D and 6C). 26 Journal of Bio-X Research ReseaRch aR ticle Figure 2. Single-cell analysis of mouse brain cells. UMAP plot showing the distribution of mouse brain cells (A). Dot plot showing the marker genes for brain clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots show- ing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. Representative immunohistochemical staining of ACE2 in cerebral vascular pericytes (D) and lung tissue of humans (E). Arrows indicate ACE2-postive cells. Violin and box plots for the expression levels of ACE2, CD147, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (F). Gene expression levels are measured as log2 (TP10K+1) values and the expression levels of gene sets are measured as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projection. 27 ReseaRch aR ticle Journal of Bio-X Research Figure 3. Single-cell analysis of healthy human nasal olfactory epithelial cells. UMAP plot showing the distribution of nasal olfactory epithelial cells (A). Dot plot showing the marker genes for nasal olfactory epithelial clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indi- cates high expression of TMPRSS or CTSL. Violin and box plots for the expression levels of ACE2, CD147, CD209, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (D). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, HBCs=horizontal basal cells, OEC=olfactory ensheathing glia cells, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projection, VSMC=vascular smooth muscle cells. Coexpression patterns of key entry genes in healthy PIKFYVE, TPCN2, and MTHFD1 were strongly expressed in lung cells (Fig. 7D). human lung tissue In total, 14 subclusters (Fig.  7A) in lung tissue were identi- fied on the basis of cell type-specific marker genes (Fig.  7B). Expression patterns of SARS-CoV-2 and neurotropic ACE2 was mainly expressed in AT2 cells; it was also expressed virus receptor genes in healthy human respiratory and in AT1 and ciliated cells (Fig. 7C and D). CD209 was mainly mouse brain tissues expressed in monocytes and macrophages, whereas CLEC4M For comprehensive analysis of susceptible cell types in the respi- was expressed in endothelial cells. CD147 was expressed in ratory system and comparison of similarities among cell types, all subclusters in lung tissue (Fig.  7D). TMPRSS and CTSL we constructed a map of the respiratory system and identified 17 were expressed in most cell types; they were not expressed in subclusters; the main subclusters comprised myeloid, alveolar, mast, NK, T, or B cells (Fig. 7C and D). Endocytosis signaling, 28 Journal of Bio-X Research ReseaRch aR ticle Figure 4. Single-cell analysis of healthy human nasal respiratory epithelial cells. UMAP plot showing the distribution of nasal respiratory epithelial cells (A). Dot plot showing the marker genes for nasal respiratory epithelial clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. (D) Violin and box plots for the expression levels of ACE2, CD147, CD209, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters. Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, HBCs=horizontal basal cells, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projection. T/NK, endothelial, ciliated, neuronal, and B cells (Additional types in the respiratory system (Fig.  8A). The endocytosis-re- Fig. 1A and B, http://links.lww.com/JR9/A45). We found that lated gene PIKFYVE was expressed mainly in myeloid cells and ACE2 was strongly expressed in ciliated cells in the OE, RE, tra- fibroblasts in the RE, as well as endothelial cells in lung tissue cheal tissue, and bronchial tissue; sustentacular cells in the OE (Fig. 8A). These results showed that factors involved in SARS- and RE; gland progenitor cells in the RE; and AT2 cells in lung CoV-2 entry are strongly expressed in the OE and RE, highlight- tissue (Fig. 8A). CD209 was strongly expressed in myeloid cells ing the potential role of such cells to facilitate viral infection in the OE, RE, and tracheal tissue. TMPRSS, CTSL, endocytosis and spread. The most susceptible cells in the respiratory system signaling, TPCN2, and MTHFD1 were expressed in most cell were ciliated, gland progenitor, and sustentacular cells, as well 29 ReseaRch aR ticle Journal of Bio-X Research Figure 5. Single-cell analysis of healthy human trachea cells. UMAP plot showing the distribution of tracheal cells (A). Dot plot showing the marker genes for tracheal clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. Violin and box plots for the expression levels of ACE2, CD147, CD209, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (D). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projec- tion, VSMC=vascular smooth muscle cells. as HBCs; AT2 cells are the known vulnerable cells for SARS- A45). Notably, many viral receptors (eg, AXL, ISG15, ITGA4, [32,33] COV2, their susceptibility could be greater than the suscep- and PECAM1) are expressed in pericytes, suggesting that peri- tibility of AT2 cells, indicating that these cells could serve as the cytes mediate CNS infection (Additional Fig. 3, http://links. portal of entry for SARS-CoV-2. lww.com/JR9/A45). Furthermore, because exosome biogenesis To explore the potential mechanism of SARS-CoV-2 infec- and release are also associated with viral transport, we explored tion in the CNS, we compared the expression patterns of entry the brain-specific and OE-specific expression patterns of genes receptors for various neurotropic viruses between maps of the involved in exosome biogenesis and release. The results showed respiratory system and brain. We found that some neurotropic that neurons, ciliated cells, and OECs had strong capacities for exosome biogenesis and release, which may serve as alterna- virus entry receptors were not expressed in neurons, indicating tive pathways for CNS infection (Fig. 8B). Thus, we speculated that other mechanisms (eg, other cells) were involved in CNS the potential mechanisms involved in CNS infection (Fig. 8C). infection (Additional Figs. 2 and 3, http://links.lww.com/JR9/ 30 Journal of Bio-X Research ReseaRch aR ticle Figure 6. Single-cell analysis of healthy human bronchial epithelial cells. UMAP plot showing the distribution of bronchial epithelial cells (A). Dot plot showing the marker genes for bronchial epithelial clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. Violin and box plots for the expression levels of ACE2 CD147, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (D). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. ACE2=angiotensin-converting enzyme 2, CTSL=cathepsin L, TMPRSS=transmembrane serine protease, UMAP=uniform mani- fold approximation and projection. [33,34] Coexpression of entry receptor ACE2 and spike protein prim- coexpression patterns of viral entry genes. Although most ing enzymes in BBB pericytes enables SARS-CoV-2 entry into analyses are based on transcriptional levels, there is clear con- the cerebrospinal fluid and brain. SARS-CoV-2 may also infect sistency at the cellular level: the expression patterns of entry the telencephalon through the olfactory bulb. genes at the transcriptional level in a particular cell type can accurately predict the protein expression levels of entry recep- [35] tors. Our results were consistent with the coexpression Discussion patterns (ACE2 and TMPTSS2) previously identified among [33] scRNA-seq is a powerful technique that can be used to iden- susceptible cell types in multiple tissues. Furthermore, we tify susceptible cells through single-cell-level analyses of the obtained novel information concerning tissues such as the OE, 31 ReseaRch aR ticle Journal of Bio-X Research Figure 7. Single-cell analysis of healthy human lung cells. UMAP plot showing the distribution of lung cells (A). Dot plot showing the marker genes for lung clusters (B). Dot scale indicates the fractions of expressing cells and dot color indicates the expression level of the marker genes. UMAP plots showing the expression patterns of ACE2, TMPRSS, and CTSL (C). Red indicates high expression of ACE2 and blue indicates high expression of TMPRSSs or CTSL. Violin and box plots for the expression levels of ACE2, CD147, CD209, CLEC4M, TMPRSS, CTSL, endocytosis-related genes, PIKFYVE, TPCN2, and MTHFD1 across clusters (D). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. CTSL=cathepsin L, DCs=dendritic cells, NK=natural killer, TMPRSS=transmembrane serine protease, UMAP=uniform manifold approximation and projection, VSMC=vascular smooth muscle cells. RE, and brain; we also identified coexpression patterns involv- The coronavirus infection process begins with viral bind- ing TMPRSS2 and viral replication genes. These results sug- ing to the cellular membrane through heparan sulfate pro- gest that neurological symptoms and CNS infection in patients teoglycans, which facilitate the interaction between the spike [16,36] with COVID-19 are secondary to SARS-CoV-2 invasion of the protein and the entry receptor (ie, ACE2). Recently, BBB via pericytes. Additionally, SARS-CoV-2–induced olfac- CD147 and CD209 were also identified as potential recep- tory disorders could be the result of localized damage to suste- tors for SARS-CoV-2; interactions with these proteins [17,18] ntacular cells and ciliated cells in the nasal RE and OE. also facilitate viral entry into host cells. In our study, 32 Journal of Bio-X Research ReseaRch aR ticle Figure 8. Expression levels of SARS-CoV-2 binding receptors, spike protein hydrolysis, endocytosis-related genes, and virus RNA replication genes across tissues. Violin and box plots for the expression levels of SARS-CoV-2 binding receptors, spike protein hydrolysis, endocytosis-related genes, and virus RNA replication genes across tissues (A). Gene expression levels are shown as log2 (TP10K+1) values, and the expression levels of gene sets are shown as mean log2 (TP10K+1) values. Box plots showing the expression levels of exosome biogenesis and release-related genes in the brain and olfactory epithelium (B). The expression levels of gene sets are shown as mean log2 (TP10K+1) values. Schematic depicting potential mechanisms involved in CNS infection (C). Coexpression of entry receptor ACE2 and spike protein priming enzyme cathepsin L in BBB pericytes enables SARS-CoV-2 entry into the cerebrospinal fluid and brain. SARS- CoV-2 may also infect the telencephalon through the olfactory bulb. ACE2=angiotensin-converting enzyme 2, BBB=blood–brain barrier, CNS=central nervous system, HBCs=horizontal basal cells, OEC=olfactory ensheathing glia cells, SARS-CoV-2=severe acute respiratory syndrome coronavirus 2. CD147 was expressed in most cell types; however, these in various cells. Thus, cells that include both the entry receptor cells were not consistent with published reports of infected (ACE2) and the spike protein priming enzyme (TMPRSS and/ [32] cells. Thus, the role of CD147 in SARS-CoV-2 pathoge- or CTSL) may be most susceptible to infection by SARS-CoV-2. nicity requires further investigation. In contrast, CD209 was Headache, impaired consciousness, epilepsy, encephalitis, mainly expressed in macrophages and monocytes in most and acute myelitis have been reported in patients with COVID- [7] tissues. The specific location of CD209 expression implies 19. Notably, a case of SARS-CoV-2 encephalitis was reported a role in SARS-CoV-2 entry into macrophages, which were by physicians in Beijing Ditan Hospital; genomic sequencing [32] also the target cells for SARS-CoV-2 infection, with the revealed the presence of SARS-CoV-2 in the patient’s cerebrospi- potential for initiation of a cytokine storm. nal fluid. Subsequently, a similar case was reported by physicians [10] Successful conformational changes and activation of spike in Japan. These reports indicated that SARS-CoV-2 is able to proteins result in membrane fusion; these processes require infect the human CNS in some instances. Intracranial infections [12] receptor binding, as well as appropriate protease activation. occur after damage to the BBB, which allows viral invasion. Endocytosis and cathepsin-mediated endosomal acidification In the present study, we found that the entry receptor ACE2 are also important processes during viral fusion with a target and spike protein priming enzyme CTSL were coexpressed in [37,38] [40] cell. In the present study, we found that TMPRSS and CTSL pericytes, which are specialized endothelial cells in the brain. were strongly expressed and exhibited a broad distribution. Furthermore, endothelial cells generally expressed lower levels Either TMPRSS or CTSL was expressed in ACE2-expressing of ACE2, which suggests that they can serve as an alternative cells, consistent with previous reports that ACE2 constitutes a route of infection. Notably, some vascular pericytes (eg, in [39] limiting factor for initial viral entry. Genes that encode pro- the OE and RE) did not express ACE2, indicating that brain teins involved in endocytosis and viral replication are expressed pericytes have a distinct role in blood pressure regulation and 33 ReseaRch aR ticle Journal of Bio-X Research [41] BBB maintenance. Endothelial cells in the brain are bound which were higher than the levels in AT2 cells. This finding together by tight junctions to ensure the integrity of the BBB. indicates that such ciliated cells may be the “portal” for SARS- The use of endothelial cells and pericytes as host cells could CoV-2 entry because of their anatomical location; they may also allow viral proliferation; the eventual cell death would destroy serve as a viral reservoir for persistent SARS-CoV-2 infection. the host cells and BBB, thereby releasing viral particles into the These results have important implications for understanding brain. Considering the extensive similarities described above in the symptoms of anosmia and potential mechanisms involved in terms of expression patterns for ACE2, TMPRSS2, and CTSL SARS-CoV-2 infection of the CNS. in mice and humans, neurons presumably do not constitute a primary site of infection; however, vascular pericytes may be Limitations sensitive to SARS-CoV-2, thus allowing entry through the BBB and subsequent CNS infection. In this study, the coexpression pattern of virus entry genes only Respiratory viruses may enter the CNS through neurons in indicates the susceptibility of the cell types. It does not really the olfactory and trigeminal nerves (within the nasal cavity), represent the cell types the SARS-CoV-2 infects. Furthermore, or through neurons in the vagus nerve (within the trachea and we analyzed susceptible cell types based on only normal sam- [42] lungs). Human coronavirus OC43 can induce CNS neu- ples. It would be more convincing if we used samples from ropathology through interneuron propagation and axonal COVID-19 patients to study the viral load in different cell types. [43] transport. We found that the entry receptor for human coro- So, the results were limited. navirus OC43 (HLA class I molecule) was expressed in neural cells, which enables the virus to infect neurons (Additional Fig. Conclusion 2, http://links.lww.com/JR9/A45). However, the entry receptor ACE2 and the spike protein priming enzymes were not coex- We identified cells susceptible to SARS-CoV-2 infection in the pressed in neural cells within the CNS, olfactory system, or CNS, olfactory system, and respiratory system. Our findings respiratory system. An earlier epidemic strain of coronavirus, provide insights into the mechanisms that underlie neurological SARS-CoV, also binds ACE2; SARS-CoV particles have been symptoms and CNS infection in patients with COVID-19. detected in brain tissue, where they were located almost exclu- [44] sively in neurons. Considering the strong capacities for exo- Acknowledgments some biogenesis and release in neurons, ciliated cells, and OECs, None. exosome biogenesis and release could be involved in SARS- CoV-2 infection of the CNS. However, this hypothesis requires additional investigation. Author contributions The OE serves as the main odor detection and conduction HX, JL: study design. ZK, JW, HZ, DX, HG, and WZ: data tissue; it includes ciliated olfactory sensory neurons, glia-like [45] analysis. ZK, DX, and HG: data collection and generation. HZ, sustentacular cells, and apical microvilli. Olfactory sensory JW, ZL, XC, and JX: data interpretation. TM, ZK, AQ, and HX: neurons are surrounded by sustentacular cells and the axons of manuscript drafting. AQ, HZ, TM, HX, and JL: overall super- OECs. Sustentacular cells serve as structural support for sensory vising and organizing the study. All authors approved the final neurons, whereas OECs guide olfactory sensory neuron axons version for publication. that cross the cribriform plate at the skull base and terminate in [46–48] the olfactory bulb. Additionally, olfactory sensory neurons [49] express odor receptors on dendritic cilia. Thus, SARS-CoV-2 Financial support invasion into ciliated cells and sustentacular cells may lead to This work was supported by the National Natural Ascience olfactory disorders. Foundation of China (No. 31821003 to HX) and the China SARS-CoV-2 primarily affects the respiratory system, with Ministry of Science and Technology (No. 2018AAA0100300 to [50] serious effects on the lungs. In this study, we observed coex- HX). The funders had no role in manuscript design, data collec- pression of ACE2 and TMPRSS/CTSL in the AT1, AT2, and cil- tion and analysis, preparation of the manuscript, or decision to iated cells in lung tissue. In lung alveoli, AT1 cells (ie, epithelial publish. cells) are responsible for gas exchange; AT2 cells (ie, alveolar stem-like cells) are responsible for surfactant biosynthesis and [51] Institutional review board statement immunoregulation. Thus, SARS-CoV-2 infection may cause damage to alveolar infrastructure; the resulting increase in There is no direct involvement of human subjects in this project. alveolar surface tension can induce dyspnea and hypoxia. The All the data use existing biological samples and data from prior [52] involvement of AT2 cells can lead to extensive inflammation. studies. Therefore, ethical oversight and patient consent were The host immune response produces a proinflammatory cyto- not handled in this study. kine storm, which manifests as acute lung injury and acute [53,54] respiratory distress syndrome. We found that CD209 was Data availability statement mainly expressed in macrophages and monocytes in lung tissue. Notably, CD209 may be responsible for viral antigen activation The single cell data are available at the Gene Expression Omnibus of specific T cells, as well as the secretion of perforin, interferon, under the accession code GSE139522 and GSE134355, Sequence and granzyme B; these processes can cause further injury to the ReadArchive dataset (SRA https://www.ncbi.nlm.nih.gov/sra) at [55] lungs and remote organs, such as the brain. accession code SRA667466 and Atlas Data Protal (https://data. Finally, we compared the expression levels of viral entry fac- humancellatlas.org) at accession code PRJEB31843. Brochus tors across olfactory and respiratory systems; we found gener- data will be made available Mendeley Data (https://data.men- ally high expression levels in ciliated cells in the nasal mucosa, deley.com/datasets/7r2cwbw44m/1). 34 Journal of Bio-X Research ReseaRch aR ticle [23] Zeisel A, Muñoz-Manchado AB, Codeluppi S, et al. Brain structure. Conflicts of interest Cell types in the mouse cortex and hippocampus revealed by single-cell There are no conflicts of interest. RNA-seq. Science 2015;347:1138–1142. doi:10.1126/science.aaa1934. [24] Madissoon E, Wilbrey-Clark A, Miragaia RJ, et al. scRNA-seq assessment of the human lung, spleen, and esophagus tissue stabil- References ity after cold preservation. Genome Biol 2019;21:1. doi:10.1186/ [1] Guan WJ, Ni ZY, Hu Y, et al. Clinical characteristics of coronavirus dis- s13059-019-1906-x. ease 2019 in China. N Engl J Med 2020;382:1708–1720. doi:10.1056/ [25] Lukassen S, Chua RL, Trefzer T, et al. SARS-CoV-2 receptor ACE2 NEJMoa2002032. and TMPRSS2 are primarily expressed in bronchial transient secretory [2] Chen N, Zhou M, Dong X, et al. Epidemiological and clinical charac- cells. EMBO J 2020;39:e105114. doi:10.15252/embj.20105114. teristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, [26] Luecken MD, Theis FJ. Current best practices in single-cell RNA- China: a descriptive study. Lancet 2020;395:507–513. doi:10.1016/ seq analysis: a tutorial. Mol Syst Biol 2019;15:e8746. doi:10.15252/ S0140-6736(20)30211-7. msb.20188746. [3] Huang C, Wang Y, Li X, et al. Clinical features of patients infected with [27] Stuart T, Butler A, Hoffman P, et al. Comprehensive integration 2019 novel coronavirus in Wuhan, China. Lancet 2020;395:497–506. of single-cell data. Cell 2019;177:1888–1902.e21. doi:10.1016/j. doi:10.1016/S0140-6736(20)30183-5. cell.2019.05.031. [4] Uginet M, Breville G, Assal F, et al. COVID-19 encephalopathy: clin- [28] Korsunsky I, Millard N, Fan J, et al. Fast, sensitive and accurate integra- ical and neurobiological features. J Med Virol 2021;93:4374–4381. tion of single-cell data with Harmony. Nat Methods 2019;16:1289–1296. doi:10.1002/jmv.26973. doi:10.1038/s41592-019-0619-0. [5] Sampaio Rocha-Filho PA, Albuquerque PM, Carvalho LCLS, et al. [29] Zhao Y, Wong L, Goh WWB. How to do quantile normalization Headache, anosmia, ageusia and other neurological symptoms in correctly for gene expression data analyses. Sci Rep 2020;10:15534. COVID-19: a cross-sectional study. J Headache Pain 2022;23:1–11. doi:10.1038/s41598-020-72664-6. doi:10.1186/s10194-021-01367-8. [30] Townes FW, Irizarry RA. Quantile normalization of single-cell RNA- [6] Dell’Era V, Farri F, Garzaro G, et al. Smell and taste disorders during seq read counts without unique molecular identifiers. Genome Biol COVID-19 outbreak: cross-sectional study on 355 patients. Head Neck 2020;21:160. doi:10.1186/s13059-020-02078-0. 2020;42:1591–1596. doi:10.1002/hed.26288. [31] Lukenda A, Dotlic S, Vukojevic N, et al. Expression and prognostic [7] Wu Y, Xu X, Chen Z, et al. Nervous system involvement after infec- value of putative cancer stem cell markers CD117 and CD15 in cho- tion with COVID-19 and other coronaviruses. Brain Behav Immun roidal and ciliary body melanoma. J Clin Pathol 2016;69:234–239. 2020;87:18–22. doi:10.1016/j.bbi.2020.03.031. doi:10.1136/jclinpath-2015-203130. [8] Helms J, Kremer S, Merdji H, et al. Neurologic features in severe SARS- [32] Bost P, Giladi A, Liu Y, et al. Host-viral infection maps reveal signa- CoV-2 Infection. N Engl J Med 2020;382:2268–2270. doi:10.1056/ tures of severe COVID-19 patients. Cell 2020;181:1475–1488.e12. NEJMc2008597. doi:10.1016/j.cell.2020.05.006. [9] Xu Z, Shi L, Wang Y, et al. Pathological findings of COVID-19 asso- [33] Sungnak W, Huang N, Bécavin C, et al. SARS-CoV-2 entry factors are ciated with acute respiratory distress syndrome. Lancet Respir Med highly expressed in nasal epithelial cells together with innate immune 2020;8:420–422. doi:10.1016/S2213-2600(20)30076-X. genes. Nat Med 2020;26:681–687. doi:10.1038/s41591-020-0868-6. [10] Moriguchi T, Harii N, Goto J, et al. A first case of meningitis/encepha- [34] Zhang H, Kang Z, Gong H, et al. Digestive system is a potential route litis associated with SARS-Coronavirus-2. Int J Infect Dis 2020;94:55– of COVID-19: an analysis of single-cell coexpression pattern of key 58. doi:10.1016/j.ijid.2020.03.062. proteins in viral entry process. Gut 2020;69:1010–1018. doi:10.1136/ [11] Shehata GA, Lord KC, Grudzinski MC, et al. Neurological complica- gutjnl-2020-320953. tions of covid-19: Underlying mechanisms and management. Int J Mol [35] Brann DH, Tsukahara T, Weinreb C, et al. Non-neuronal expression of Sci 2021;22:4081. doi:10.3390/ijms22084081. SARS-CoV-2 entry genes in the olfactory system suggests mechanisms [12] Li Z, Liu T, Yang N, et al. Neurological manifestations of patients with underlying COVID-19-associated anosmia. Sci Adv 2020;6:eabc5801. COVID-19: potential routes of SARS-CoV-2 neuroinvasion from the doi:10.1126/sciadv.abc5801. periphery to the brain. Front Med 2020;14:533–541. doi:10.1007/ [36] Wan Y, Shang J, Graham R, et al. Receptor recognition by the novel s11684-020-0786-5. coronavirus from Wuhan: an analysis based on decade-long struc- [13] Shang J, Wan Y, Luo C, et al. Cell entry mechanisms of SARS-CoV-2. tural studies of SARS coronavirus. J Virol 2020;94:e00127–e00120. Proc Natl Acad Sci U S A 2020;117:11727–11734. doi:10.1073/ doi:10.1128/JVI.00127-20. pnas.2003138117. [37] Knyazev E, Nersisyan S, Tonevitsky A. Endocytosis and transcytosis of [14] Chen M, Shen W, Rowan NR, et al. Elevated ACE-2 expression in the SARS-CoV-2 across the intestinal epithelium and other tissue barriers. olfactory neuroepithelium: implications for anosmia and upper respira- Front Immunol 2021;12:636966. doi:10.3389/fimmu.2021.636966. tory SARS-CoV-2 entry and replication. Eur Respir J 2020;56:2001948. [38] Zhao MM, Yang WL, Yang FY, et al. Cathepsin L plays a key role in doi:10.1183/13993003.01948-2020. SARS-CoV-2 infection in humans and humanized mice and is a prom- [15] Zhou P, Yang XL, Wang XG, et al. A pneumonia outbreak associated with ising target for new drug development. Signal Transduct Target Ther a new coronavirus of probable bat origin. Nature 2020;579:270–273. 2021;6:134. doi:10.1038/s41392-021-00558-8. doi:10.1038/s41586-020-2012-7. [39] Ou X, Liu Y, Lei X, et al. Characterization of spike glycoprotein of [16] Wrapp D, Wang N, Corbett KS, et al. Cryo-EM structure of the 2019- SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS- nCoV spike in the prefusion conformation. Science 2020;367:1260–1263. CoV. Nat Commun 2020;11:1620. doi:10.1038/s41467-020-15562-9. doi:10.1126/science.abb2507. [40] Campisi M, Shin Y, Osaki T, et al. 3D self-organized microvascular [17] Wang K, Chen W, Zhang Z, et al. CD147-spike protein is a novel route model of the human blood-brain barrier with endothelial cells, peri- for SARS-CoV-2 infection to host cells. Signal Transduct Target Ther cytes and astrocytes. Biomaterials 2018;180:117–129. doi:10.1016/j. 2020;5:283. doi:10.1038/s41392-020-00426-x. biomaterials.2018.07.014. [18] Amraei R, Yin W, Napoleon MA, et al. CD209L/L-SIGN and CD209/ [41] Hirunpattarasilp C, Attwell D, Freitas F. The role of pericytes in brain DC-SIGN act as receptors for SARS-CoV-2. ACS Cent Sci 2021;7:1156– disorders: from the periphery to the brain. J Neurochem 2019;150:648– 1165. doi:10.1021/acscentsci.0c01537. 665. doi:10.1111/jnc.14725. [19] Hoffmann M, Kleine-Weber H, Schroeder S, et al. SARS-CoV-2 cell [42] Mori I. Transolfactory neuroinvasion by viruses threatens the human entry depends on ACE2 and TMPRSS2 and is blocked by a clinically brain. Acta Virol 2015;59:338–349. doi:10.4149/av_2015_04_338. proven protease inhibitor. Cell 2020;181:271–280.e8. doi:10.1016/j. [43] Dubé M, Le Coupanec A, Wong AHM, et al. Axonal transport enables cell.2020.02.052. neuron-to-neuron propagation of human coronavirus OC43. J Virol [20] Meng T, Cao H, Zhang H, et al. The insert sequence in SARS- 2018;92:e00404–e00418. doi:10.1128/JVI.00404-18. CoV-2 enhances spike protein cleavage by TMPRSS. bioRxiv 2020. [44] Xu J, Zhong S, Liu J, et al. Detection of severe acute respiratory syndrome doi:10.1101/2020.02.08.926006. coronavirus in the brain: Potential role of the chemokine mig in patho- [21] Durante MA, Kurtenbach S, Sargi ZB, et al. Single-cell analysis of olfac- genesis. Clin Infect Dis 2005;41:1089–1096. doi:10.1086/444461. tory neurogenesis and differentiation in adult humans. Nat Neurosci [45] Villar PS, Delgado R, Vergara C, et al. Energy requirements of odor 2020;23:323–326. doi:10.1038/s41593-020-0587-9. transduction in the chemosensory cilia of olfactory sensory neu- [22] Han X, Zhou Z, Fei L, et al. Construction of a human cell land- rons rely on oxidative phosphorylation and glycolytic processing of scape at single-cell level. Nature 2020;581:303–309. doi:10.1038/ extracellular glucose. J Neurosci 2017;37:5736–5743. doi:10.1523/ s41586-020-2157-4. JNEUROSCI.2640-16.2017. 35 ReseaRch aR ticle Journal of Bio-X Research [46] Kondoh K, Lu Z, Ye X, et al. A specific area of olfactory cortex involved [51] Nabhan AN, Brownfield DG, Harbury PB, et al. Single-cell Wnt sig- in stress hormone responses to predator odours. Nature 2016;532:103– naling niches maintain stemness of alveolar type 2 cells. Science 106. doi:10.1038/nature17156. 2018;359:1118–1123. doi:10.1126/science.aam6603. [47] Le Bourhis M, Rimbaud S, Grebert D, et al. Endothelin uncouples gap [52] Kroetz DN, Allen RM, Schaller MA, et al. Type I interferon induced epi- junctions in sustentacular cells and olfactory ensheathing cells of the genetic regulation of macrophages suppresses innate and adaptive immu- olfactory mucosa. Eur J Neurosci 2014;40:2878–2887. doi:10.1111/ nity in acute respiratory viral infection. PLoS Pathog 2015;11:e1005338. ejn.12665. doi:10.1371/journal.ppat.1005338. [48] Beites CL, Kawauchi S, Crocker CE, et al. Identification and molecular [53] Desforges M, Le Coupanec A, Dubeau P, et al. Human coronaviruses and regulation of neural stem cells in the olfactory epithelium. Exp Cell Res other respiratory viruses: underestimated opportunistic pathogens of the 2005;306:309–316. doi:10.1016/j.yexcr.2005.03.027. central nervous system? Viruses 2019;12:14. doi:10.3390/v12010014. [49] Savya SP, Kunkhyen T, Cheetham CEJ. Low survival rate of young [54] Mehta P, McAuley DF, Brown M, et al. COVID-19: consider cytokine adult-born olfactory sensory neurons in the undamaged mouse olfac- storm syndromes and immunosuppression. Lancet 2020;395:1033– tory epithelium. J Bioenerg Biomembr 2019;51:41–51. doi:10.1007/ 1034. doi:10.1016/S0140-6736(20)30628-0. s10863-018-9774-8. [55] Channappanavar R, Perlman S. Pathogenic human coronavirus [50] Guan WJ, Liang WH, Zhao Y, et al. Comorbidity and its impact on infections: causes and consequences of cytokine storm and immuno- 1590 patients with COVID-19 in China: a nationwide analysis. Eur pathology. Semin Immunopathol 2017;39:529–539. doi:10.1007/ Respir J 2020;55:2000547. doi:10.1183/13993003.00547-2020. s00281-017-0629-x.

Journal

Journal of Bio-X ResearchWolters Kluwer Health

Published: Mar 22, 2023

References