key: cord-266022-aco5kpaj
authors: Matusiak, Magdalena; Schürch, Christian M.
title: Expression of SARS-CoV-2 entry receptors in the respiratory tract of healthy individuals, smokers and asthmatics
date: 2020-09-29
journal: Respir Res
DOI: 10.1186/s12931-020-01521-x
sha: 
doc_id: 266022
cord_uid: aco5kpaj

SARS-CoV-2 is causing a pandemic with currently > 29 million confirmed cases and > 900,000 deaths worldwide. The locations and mechanisms of virus entry into the human respiratory tract are incompletely characterized. We analyzed publicly available RNA microarray datasets for SARS-CoV-2 entry receptors and cofactors ACE2, TMPRSS2, BSG (CD147) and FURIN. We found that ACE2 and TMPRSS2 are upregulated in the airways of smokers. In asthmatics, ACE2 tended to be downregulated in nasal epithelium, and TMPRSS2 was upregulated in the bronchi. Furthermore, respiratory epithelia were negative for ACE-2 and TMPRSS2 protein expression while positive for BSG and furin, suggesting a possible alternative entry route for SARS-CoV-2.

expression measure method using affy and limma packages in R [17, 18] . For genes represented by multiple probes, the probe with the maximum average expression values in all samples was selected to represent that gene's expression. First, by plotting the first 2 principal components computed on ACE2, TMPRSS2, BSG and FURIN expression across smokers' and asthmatics' datasets, we verified that there were no detectable batch effects within each of the six microarray datasets we sought to analyze (Figs. S1A-D and S2A-D). Next, differences in gene expression in smokers vs. non-smokers and asthmatics vs. healthy individuals were modeled using linear regression, including proband age, sex, sample type and dataset as covariates in the model. When modeling the age effect with linear regression, the age of probands for whom age information was not available was set to the average age of all other probands (separately for smokers' datasets GSE63127 and asthmatics' dataset GSE4302; for age distribution, see Figs. S1E and S2E). P values were next corrected for multiple hypothesis testing using the Benjamini-Hochberg correction. In post hoc analysis, log 2 transformed data not corrected for age and sex were plotted with ggpubr [19] , and two sided Mann-Whitney U tests were performed using the Wilcox.test function in R.

Older age is an important risk factor for adverse COVID-19 outcomes [20] . Another risk factor is male sex [21] . Since there is a very high smoking rate in Chinese males compared to females-66.1% vs. 3.2% according to Ma et al. [22] -it has been suggested that smoking could be a risk factor for the more severe COVID-19 disease course observed in males [20, 23] ; however, this topic is controversially discussed [24] [25] [26] . In our analysis, we therefore removed possible confounding effects of age and sex on the status of receptor mRNA expression by regressing the linear effects of age and sex and testing our hypothesis on model residuals (Table 1 ). For samples for which sex information was not available, sex was predicted based on the expression of both X inactive specific transcript (XIST; high expression in females) and ribosomal protein S4 Y-linked 1 (RPS4Y1; high expression in males) simultaneously. In post hoc analysis, we also plotted log 2 transformed expression values not corrected for age and sex (Fig. 1) . Consistent with previous reports [7, 27] , we found significantly higher ACE2 expression in airway epithelia from healthy smokers vs. healthy non-smokers (Fig. 1a , Table 1 ). Similarly, we found significantly higher TMPR SS2 expression in smokers in one out of two datasets analyzed ( Fig. 1b, Table 1 ), whereas BSG and FURIN expression did not significantly differ between smokers and non-smokers ( Fig. 1c -d, Table 1 ). Taken together, these results indicate that, independently of sex and age, ACE2 and TMPRSS2 are upregulated in the airway epithelia of smokers. In addition, our analyses establish that BSG and FURIN, two alternate potential SARS-CoV-2 receptors, are expressed in the human respiratory tract.

We next examined four RNA microarray datasets for ACE2, TMPRSS2, BSG and FURIN expression in airway epithelia from patients with a common respiratory disease, asthma. Patients with chronic respiratory disorders including asthma are considered a COVID-19 high-risk category [3] . Interestingly, we found that ACE2 expression tended to be downregulated in nasal epithelium, whereas TMPRSS2 was significantly upregulated in bronchi and central airways of asthmatics ( Fig. 1e -f, Table 1 ). ACE2 was proposed to be an interferonstimulated gene [28] ; therefore, a potential explanation for ACE2 downregulation in asthmatics could be corticosteroid use. However, more recent data indicate that a novel, primate-specific ACE2 isoform exists that is incapable of binding SARS-CoV-2, and that this isoform is interferon-stimulated, whereas the canonical ACE2 is not [29, 30] . Further research is needed to address these highly interesting developments in more detail.

Additionally, we did not find any difference in BSG or FURIN expression between healthy and asthmatic individuals ( Fig. 1g-h, Table 1 ). These findings point towards a possible differential regulation of ACE2 and TMPRSS2 expression in airway epithelia and warrant further investigation into the underlying mechanism.

We next aimed to compare RNA expression of these receptors and cofactors to protein expression. Therefore, we examined immunohistochemistry (IHC) images from respiratory and other tissues on The Human Protein Atlas [31] . ACE-2 IHC staining was strong in epithelial cells of the duodenum and was found in other organs including testis and kidney. In contrast, human respiratory epithelial cells in samples from the nasopharynx, bronchi and lungs, as well as squamous epithelial cells from the oral mucosa, were completely negative for ACE-2 staining by IHC with two different antibodies (Fig. 2a-b) . Similar negative IHC staining results were also observed for TMPRSS2 protein (Fig. 2c) . In contrast, basigin protein was widely expressed in human tissues including heart muscle, brain, liver and kidney, and, importantly, was positive in respiratory epithelial cells from the nasopharynx and bronchi (Fig. 2d-e) . Similar to ACE-2 and TMPRSS2, basigin was negative in alveolar epithelial cells but showed multifocal positivity in cells morphologically consistent with alveolar macrophages (Fig. 2d-e) . Furthermore, basigin was strongly expressed in the oral mucosa, with a gradient from basal to apical epithelial cells (Fig. 2d-e) . Finally, furin IHC stainings were positive in nasopharynx and bronchial epithelial cells, and negative in alveolar epithelial cells, with focal positivity in alveolar macrophages (Fig. 2f) . In summary, basigin and furin protein were expressed in airway epithelia of nasopharynx and bronchi, whereas ACE-2 and TMPRSS2 protein stainings were negative. Our findings are in line with and extend recent studies addressing SARS-CoV-2 receptor and cofactor expression in the respiratory tract [32, 33] , but are in stark contrast to a 2004 study by Hamming et al. [34] , who found strong and widespread ACE-2 expression in alveolar epithelial cells and basal epithelial cells of the nasopharynx and oral mucosa. Our study highlights the discrepancies between RNA and protein expression of these receptors and cofactors, and points towards potential issues with IHC staining reproducibility and antibody specificity, important factors that need to be addressed in future investigations. One limitation of our study was that in The Human Protein Atlas, only small numbers of IHC stained samples for each tissue and molecule analyzed were available. Therefore, further studies exploring the protein expression and cellular localization of SARS-CoV-2 receptors and cofactors in each of these tissue types, ideally using tissue microarrays with large numbers of samples from multiple donor cohorts, and using multiple different antibodies, are warranted. In addition, these studies should use the recently developed multiplexed microscopy technologies [35] to address protein co-expression patterns and better delineate the cellular subsets expressing these SARS-CoV-2 entry receptors and cofactors.

Supplementary information accompanies this paper at https://doi.org/10. 1186/s12931-020-01521-x.

Abbreviations SARS-CoV-2: Severe acute respiratory syndrome coronavirus 2; COVID-19: Coronavirus disease 2019; ACE-2: Angiotensin I converting enzyme 2; TMPRSS2: Transmembrane serine protease 2; BSG: Basigin; XIST: X inactive specific transcript; RPS4Y1: Ribosomal protein S4 Y-linked 1

Datasets for receptor protein immunohistochemistry

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

Coronavirus COVID-19 Global Cases by the Center for Systems Science and Engineering (CSSE)

Preliminary estimates of the prevalence of selected underlying health conditions among patients with coronavirus disease 2019 -United States

Post-mortem examination of COVID19 patients reveals diffuse alveolar damage with severe capillary congestion and variegated findings of lungs and other organs suggesting vascular dysfunction

SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor

Immunology of COVID-19: current state of the science

ACE-2 Expression in the Small Airway Epithelia of Smokers and COPD Patients: Implications for COVID-19

SARS-CoV-2 entry factors are highly expressed in nasal epithelial cells together with innate immune genes

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

Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein

Genome-wide profiling identifies epithelial cell genes associated with asthma and with treatment response to corticosteroids

Smoking-mediated up-regulation of GAD67 expression in the human airway epithelium

Upper airway gene expression is an effective surrogate biomarker for Th2-driven inflammation in the lower airway

Asthma-COPD overlap. Clinical relevance of genomic signatures of type 2 inflammation in chronic obstructive pulmonary disease

Cigarette smoking induces changes in airway epithelial expression of genes associated with monogenic lung disorders

Altered Epithelial Gene Expression in Peripheral Airways of Severe Asthma

Affy-analysis of Affymetrix GeneChip data at the probe level

Limma powers differential expression analyses for RNA-sequencing and microarray studies

ggpubr: 'ggplot2' Based Publication Ready Plots

Covid-19: risk factors for severe disease and death

Clinical characteristics of refractory COVID-19 pneumonia in Wuhan, China

Prevalence of Cigarette Smoking and Nicotine Dependence in Men and Women Residing in Two Provinces in China

Sex difference and smoking predisposition in patients with COVID-19

Active smoking is not associated with severity of coronavirus disease 2019 (COVID-19)

COVID-19 and smoking: A systematic review of the evidence

Current Smoking is Not Associated with COVID-19

Bulk and single-cell transcriptomics identify tobacco-use disparity in lung gene expression of ACE2, the receptor of 2019-nCov

SARS-CoV-2 Receptor ACE2 Is an Interferon-Stimulated Gene in Human Airway Epithelial Cells and Is Detected in Specific Cell Subsets across Tissues

Interferons and viruses induce a novel primate-specific isoform dACE2 and not the SARS-CoV-2 receptor ACE2

Tissue-specific and interferon-inducible expression of non-functional ACE2 through endogenous retrovirus co-option

Towards a knowledge-based human protein atlas

The protein expression profile of ACE2 in human tissues

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

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

Coordinated cellular neighborhoods orchestrate Antitumoral immunity at the colorectal Cancer invasive front

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Not applicable. 

All data presented are from publicly available datasets, as detailed in the manuscript and supplemental information. Datasets for smokers vs. non-smokers:

The authors declare that they have no competing interests.Received: 25 June 2020 Accepted: 23 September 2020