To the Editor:
On March 11, 2020, the World Health Organization declared the coronavirus disease (COVID-19) outbreak a pandemic. As of May 28, 2020, laboratories had confirmed 5,701,337 COVID-19 cases, and 357,668 deaths had been reported in 216 countries, areas, or territories (1). COVID-19 is caused by a new type of pathogenic coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is phylogenetically similar to SARS-CoV, with approximately 80% identity between the genomes (2). SARS viruses affect the respiratory tract and cause an acute respiratory response through the same cell-entry receptor, ACE2 (angiotensin-converting enzyme 2), which is the only experimentally confirmed SARS-CoV-2 receptor. SARS-CoV-2 infection also uses activation of the spike proteins found on the surface of the virus for cellular entry. The best candidates for priming spike proteins are two host cell enzymes called Furin and TMPRSS2 (2). In the current severe global emergency, to enable effective prevention and care, it is imperative to identify potential risk factors, such as cigarette smoking, which is a substantial risk factor for various important bacterial and viral infections.
Some of the results of this study were previously published in preprint form (medRxiv, https://www.medrxiv.org/content/10.1101/2020.02.05.20020107v3).
We evaluated a comprehensive set of transcriptomic data sets to investigate the associations of smoking with ACE2, FURIN, and TMPRSS2 gene expression in lung tissues. Two data sets were generated using normal lung tissues from patients with lung adenocarcinoma: a Caucasian RNA-sequencing (RNA-seq) data set from The Cancer Genome Atlas (n = 48) (3) and an Asian RNA-seq data set from the Gene Expression Omnibus (GSE40419, n = 74) (4). We included three polyethnic microarray data sets of gene expression in healthy small airway epithelium samples (GSE63127, n = 230 [5]; GSE19667, n = 116 [6]; and GSE5058, n = 24 [7]) and large airway epithelium samples (GSE7895, n = 104 [8]). In addition, we analyzed three microarray data sets of samples derived from healthy subjects and patients with chronic obstructive pulmonary disease (COPD), including small airway epithelium samples from current smokers (from GSE5058, n = 26 [7]), bronchial airway epithelium samples from current and former smokers (GSE37147, n = 238 [9]), and lung samples from white patients (n = 438) who underwent lung cancer surgery at the Institut Universitaire de Cardiologie et de Pneumologie de Québec (10). RNA-seq data sets were generated with the Illumina HiSeq platform and microarray data sets were generated with Affymetrix arrays. All of the data were deidentified, and the study approvals were obtained in the original studies (2–10). A total of 1,286 assay results were evaluated. We considered the fragments per kilobase per million mapped reads for RNA-seq data and robust multiarray average values for microarray data to represent normalized gene expression. All data were log2 transformed to improve normality.
Smoking status (never, former, or current smoker) was identified based on self-reported smoking history. Association tests were performed using a linear model with log2 ACE2, FURIN, or TMPRSS2 gene expression as the dependent variable and smoking status or COPD status as an independent variable. A meta-analysis was performed by pooling the effect sizes and SEs estimated from each study using a random-effects model. Age and sex were included as covariates. Although we did not observe significant associations of age and sex with the expression of ACE2 and FURIN, we found a negative correlation between TMPRSS2 expression and age in some of the data sets. Data management, statistical analyses, and visualizations were performed using R 3.6.1.
A single-cell RNA-seq data set (GSE131391) (11) was also analyzed. In this analysis, bronchial epithelial cells, single ALCAM+ epithelial cells, and CD45+ white blood cells were profiled from six never-smokers and six current smokers. Sequencing read counts in single cells were downloaded, and subsequent data analyses, including data normalization, high variable feature selection, data scaling, dimension reduction, and cluster identification, were performed using the Seurat 3.0 package. We used SCANNER for data visualization and cell type identification.
We identified upregulation of pulmonary ACE2 gene expression in ever-smokers compared with nonsmokers in all data sets, irrespective of tissue subset or COPD status (Figure 1). A meta-analysis showed that ever-smoking significantly and substantially increased pulmonary ACE2 expression by 25% (P value = 1.4 × 10−16; Figure 1). Similarly, smoking status (never, former, or current smoker) was also significantly associated with ACE2 pulmonary expression in the meta-analysis (β = 0.14, P = 2.0 × 10−6; Figure 1). The significant effect of smoking on ACE2 pulmonary expression identified in this study may suggest an increased risk for viral binding and entry of SARS-CoV and SARS-CoV-2 in lungs of smokers. FURIN was also upregulated by smoking, but to a lower extent compared with ACE2. TMRPSS2 gene expression in lung was not associated with smoking (Figure 1).
We also evaluated the effect of COPD on gene expression. When we stratified the data by smoking status, we observed a trend (β = 0.08, P = 0.07) for higher ACE2 levels in patients with COPD, but the results were not consistent across data sets (Figure 1). In the Institut Universitaire de Cardiologie et de Pneumologie de Québec data, ACE2 expression was upregulated in patients with COPD (P = 0.0006), but the effect was attenuated after adjustment for smoking status (P = 0.03).
We further evaluated the effect of smoking on ACE2 pulmonary expression in single bronchial epithelial cells from six never-smokers and six current smokers. We found that smoking remodeled cells in the bronchial epithelium, with a loss of club cells and extensive hyperplasia of goblet cells. The ACE2 gene was mainly expressed in goblet cells in smokers, and in club cells in never-smokers (Figure 2). This result is consistent with a very recent study that found the highest ACE2 expression in alveolar type II cells (which derive from club cells) and in a transient secretory cell type in subsegmental bronchial branches (12). This may indicate that smokers have a risk of COVID-19 infection complications based on their ACE2 expression profiles, which could contribute to variations in infection susceptibility, disease severity, and treatment outcome.
Despite a significant increase in the prevalence of electronic cigarettes (E-cigs), to date no studies have compared single (E-cig only) and dual (E-cig and tobacco) users. The mechanisms underlying tobacco-related upregulation of ACE2 pulmonary expression, as well as the degree to which smoking affects infection susceptibility and clinical manifestations, are unknown. Further mechanistic studies are needed to address these issues. Although our knowledge is currently limited, this study indicates that smoking could be a risk factor for COVID-19 by affecting ACE2 expression, and provides valuable information for identifying and stratifying more susceptible populations.
1. | World Health Organization. Coronavirus disease (COVID-19) outbreak situation [updated 2020 May 28; accessed 2020 May 29]. Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019. |
2. | Walls AC, Park YJ, Tortorici MA, Wall A, McGuire AT, Veesler D. Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell 2020;181:281–292, e6. |
3. | Weinstein JN, Collisson EA, Mills GB, Shaw KR, Ozenberger BA, Ellrott K, et al.; Cancer Genome Atlas Research Network. The Cancer Genome Atlas Pan-Cancer Analysis Project. Nat Genet 2013;45:1113–1120. |
4. | Seo JS, Ju YS, Lee WC, Shin JY, Lee JK, Bleazard T, et al. The transcriptional landscape and mutational profile of lung adenocarcinoma. Genome Res 2012;22:2109–2119. |
5. | Tilley AE, Staudt MR, Salit J, Van de Graaf B, Strulovici-Barel Y, Kaner RJ, et al. Cigarette smoking induces changes in airway epithelial expression of genes associated with monogenic lung disorders. Am J Respir Crit Care Med 2016;193:215–217. |
6. | Strulovici-Barel Y, Omberg L, O’Mahony M, Gordon C, Hollmann C, Tilley AE, et al. Threshold of biologic responses of the small airway epithelium to low levels of tobacco smoke. Am J Respir Crit Care Med 2010;182:1524–1532. |
7. | Tilley AE, Harvey BG, Heguy A, Hackett NR, Wang R, O’Connor TP, et al. Down-regulation of the notch pathway in human airway epithelium in association with smoking and chronic obstructive pulmonary disease. Am J Respir Crit Care Med 2009;179:457–466. |
8. | Beane J, Sebastiani P, Liu G, Brody JS, Lenburg ME, Spira A. Reversible and permanent effects of tobacco smoke exposure on airway epithelial gene expression. Genome Biol 2007;8:R201. |
9. | Steiling K, van den Berge M, Hijazi K, Florido R, Campbell J, Liu G, et al. A dynamic bronchial airway gene expression signature of chronic obstructive pulmonary disease and lung function impairment. Am J Respir Crit Care Med 2013;187:933–942. |
10. | Bossé Y, Postma DS, Sin DD, Lamontagne M, Couture C, Gaudreault N, et al. Molecular signature of smoking in human lung tissues. Cancer Res 2012;72:3753–3763. |
11. | Duclos GE, Teixeira VH, Autissier P, Gesthalter YB, Reinders-Luinge MA, Terrano R, et al. Characterizing smoking-induced transcriptional heterogeneity in the human bronchial epithelium at single-cell resolution. Sci Adv 2019;5:eaaw3413. |
12. | Lukassen S, Chua RL, Trefzer T, Kahn NC, Schneider MA, Muley T, et al. SARS-CoV-2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells. EMBO J [online ahead of print] 4 Apr 2020; DOI: https://doi.org/10.15252/embj.20105114. |
C.A. was supported by Cancer Prevention Research Institute of Texas grant RR170048 and by NIH/National Cancer Institute grant U19CA203654. F.K. was supported by NIH grants R01 ES029442, R01 AI135803, and W81XWH-16-1-0361 and by VA Merit grant CX000104.
Author Contributions: Conception and design: G.C. and C.I.A. Analysis and interpretation: G.C., Y.B., and C.I.A. Drafting of the manuscript for important intellectual content: G.C., Y.B., F.X., F.K., and C.I.A.
Originally Published in Press as DOI: 10.1164/rccm.202003-0693LE on April 24, 2020
Author disclosures are available with the text of this letter at www.atsjournals.org.