Association of respiratory allergy, asthma, and expression of the SARS-CoV-2 receptor ACE2
Authors: Daniel J. Jackson, William W. Busse, Leonard B. Bacharier, Meyer Kattan, George T. O’Connor, Robert A. Wood, Cynthia M. Visness, Stephen R. Durham, David Larson, Stephane Esnault, Carole Ober, Peter J. Gergen, Patrice Becker, Alkis Togias, James E. Gern, Mathew C. Altman
Journal: Journal of Allergy and Clinical Immunology (2020)

The novel coronavirus SARS-CoV-2 (COVID-19) was recognized in December 2019 as a cause of severe pneumonia and has now led to a global pandemic.1Zhou P. Yang X.L. Wang X.G. Hu B. Zhang L. Zhang W. et al.A pneumonia outbreak associated with a new coronavirus of probable bat origin.Nature. 2020; 579: 270-273Crossref PubMed Scopus (14613) Google Scholar Respiratory illnesses caused by COVID-19 cover a range of severity. The identification of risk and protective factors for disease severity from COVID-19 is critical to direct development of new treatments and infection prevention strategies. Early large case series have identified a number of risk factors for severe disease, including older age, hypertension, diabetes, cardiovascular disease, tobacco exposure, and chronic obstructive pulmonary disease.2Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72314 cases from the Chinese Center for Disease Control and Prevention [e-pub ahead of print]. JAMA https://doi.org/10.1001/jama.2020.2648. Accessed March 25, 2020.Google Scholar The US Centers for Disease Control and Prevention lists asthma as a risk factor for severe COVID-19 illness, which is logical given that many respiratory viruses have been well established to cause more serious illnesses in those with chronic airway diseases such as asthma. However, asthma and respiratory allergy have not been identified as significant risk factors for severe COVID-19 illness in case series from China.2Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72314 cases from the Chinese Center for Disease Control and Prevention [e-pub ahead of print]. JAMA https://doi.org/10.1001/jama.2020.2648. Accessed March 25, 2020.Google Scholar These preliminary reports led us to question whether we could identify features of allergy and/or asthma that could be associated with potential for reduced severity of COVID-19 illness. SARS-CoV-2 uses angiotensin-converting enzyme 2 (ACE2) as its cellular receptor, as do SARS-CoV and the coronavirus NL63.1Zhou P. Yang X.L. Wang X.G. Hu B. Zhang L. Zhang W. et al.A pneumonia outbreak associated with a new coronavirus of probable bat origin.Nature. 2020; 579: 270-273Crossref PubMed Scopus (14613) Google Scholar Higher ACE2 expression increases in vitro susceptibility to SARS-CoV,3Jia H.P. Look D.C. Shi L. Hickey M. Pewe L. Netland J. et al.ACE2 receptor expression and severe acute respiratory syndrome coronavirus infection depend on differentiation of human airway epithelia.J Virol. 2005; 79: 14614-14621Crossref PubMed Scopus (663) Google Scholar and studies examining factors that affect ACE2 gene expression have revealed that its upregulation is associated with smoking, diabetes, and hypertension, all of which are associated with increased severity of COVID-19 illness.4Brake S.J. Barnsley K. Lu W. McAlinden K.D. Eapen M.S. Sohal S.S. Smoking upregulates angiotensin-converting enzyme-2 receptor: a potential adhesion site for novel coronavirus SARS-CoV-2 (Covid-19).J Clin Med. 2020; 9: 841Crossref PubMed Scopus (368) Google Scholar We hypothesized that 1 potential explanation for the unexpected observation that asthma and other allergic diseases may not be a risk factor for severe COVID-19 disease is a reduced ACE2 gene expression in airway cells and thus decreased susceptibility to infection. To test this hypothesis, we examined whether asthma and respiratory allergy are associated with reduced ACE2 expression in airway cells from 3 different cohorts of children and adults. In all 3 studies, total RNA was extracted from nasal or lower airway epithelial brush samples, with RNA sequencing performed independently for each study as previously described and provided in detail in the Supplementary Information (available in this article's Online Repository at www.jacionline.org).5Larson D, Patel P, Salapatek AM, Couroux P, Whitehouse D, Pina A, et al. Nasal allergen challenge and environmental exposure chamber challenge: a randomized trial comparing clinical and biological responses to cat allergen. J Allergy Clin Immunol https://doi.org/10.1016/j.jaci.2020.02.024. Accessed March 25, 2020.Google Scholar Differential expression of ACE2 was assessed by using a weighted linear mixed effects model (limma) appropriate for RNA sequencing data and an empiric Bayes method. Children at high risk for asthma based on parental histories and living in urban neighborhoods were enrolled prenatally and followed prospectively in the Urban Environment and Childhood Asthma (URECA) cohort; 318 of them had nasal epithelial brushes obtained at 11 years of age. Prevalence of asthma was assessed at 10 years of age, and atopic status was defined by allergic sensitization trajectories (no or minimal, low, medium, and high) as previously described.6Bacharier L.B. Beigelman A. Calatroni A. Jackson D.J. Gergen P.J. O'Connor G.T. et al.Longitudinal phenotypes of respiratory health in a high-risk urban birth cohort.Am J Respir Crit Care Med. 2019; 199: 71-82Crossref PubMed Scopus (67) Google Scholar Additional type 2 biomarkers, including fractional exhaled nitric oxide, peripheral blood eosinophil level, and total IgE level, were evaluated by using standard methods. In URECA, allergic sensitization was inversely related to ACE2 expression in the nasal epithelium regardless of asthma status (Fig 1, A). In children with asthma, moderate allergic sensitization (fold change [FC] = 0.70; P = 4.2E–3) and high allergic sensitization (FC = 0.54; P = 6.4E–5) were associated with progressively greater reductions in ACE2 expression compared with in children with asthma but no/minimal allergic sensitization (Fig 1, B). ACE2 expression was also significantly inversely associated with type 2 biomarkers (see Table E1 in this article's Online Repository at www.jacionline.org), including the number of positive allergen-specific IgE test results (β coefficient –0.089; P = 3.1E–5), total IgE level (β coefficient –0.31; P = 5.1E–6), fractional exhaled nitric oxide (β coefficient –0.45; P = 3.4E–3), and nasal epithelial expression of IL13 (β coefficient –0.123; P = 8.6E–5). ACE2 expression was not significantly correlated with peripheral blood eosinophil level (β coefficient –0.13; P = .07). Although male sex has been associated with increased severity of COVID-19 illness,2Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72314 cases from the Chinese Center for Disease Control and Prevention [e-pub ahead of print]. JAMA https://doi.org/10.1001/jama.2020.2648. Accessed March 25, 2020.Google Scholar no sex-based differences in ACE2 expression were found in URECA. Of note, 10 participants reported nasal corticosteroid use at the time of nasal sampling, and it was not associated with alterations in ACE2 expression. We also evaluated 24 adult participants with allergic rhinitis to cat who had no asthma symptoms in the prior year, were enrolled in a study in which they underwent nasal cat allergen challenge (NAC), and had been exposed to cat allergen through an environmental exposure chamber (EEC), as previously described.5Larson D, Patel P, Salapatek AM, Couroux P, Whitehouse D, Pina A, et al. Nasal allergen challenge and environmental exposure chamber challenge: a randomized trial comparing clinical and biological responses to cat allergen. J Allergy Clin Immunol https://doi.org/10.1016/j.jaci.2020.02.024. Accessed March 25, 2020.Google Scholar Pre–allergen challenge and post–allergen challenge nasal brush samples were obtained. Allergen exposure by both NAC and EEC led to significant reductions in ACE2 expression (Fig 2, A) (with NAC, FC = 0.81 and P = 2.4E–3; with exposure through an EEC, FC = 0.79 and P = 1.6E–3). An additional cohort of 23 adult participants with mild asthma that was not treated with asthma controller therapy underwent segmental allergen bronchoprovocation to dust mite, ragweed, or cat, as previously described.7Kelly E.A. Esnault S. Liu L.Y. Evans M.D. Johansson M.W. Mathur S. et al.Mepolizumab attenuates airway eosinophil numbers, but not their functional phenotype, in asthma.Am J Respir Crit Care Med. 2017; 196: 1385-1395Crossref PubMed Scopus (99) Google Scholar Pre–allergen challenge and post–allergen challenge bronchial brushings were obtained and demonstrated significantly reduced ACE2 expression in lower airway epithelium in the post–allergen challenge samples (Fig 2, B) (FC 0.64; P = .01). From in vitro models obtained from the Gene Expression Omnibus, we assessed the effects of IL-13, a type 2 cytokine strongly related to allergic asthma, on ACE2 expression in differentiated airway epithelial cells. IL-13 significantly reduced ACE2 expression (see Fig E1 in this article's Online Repository at www.jacionline.org) in both nasal (FC = 0.44; P = 5.8E–4) and bronchial (FC = 0.80; P = 5.1E–3) epithelium. Viral respiratory infections are the most common trigger of severe asthma exacerbations in children and adults. Unexpectedly, large epidemiologic studies of the COVID-19 pandemic in China did not identify asthma as a risk factor of severe COVID-19–related illnesses.2Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72314 cases from the Chinese Center for Disease Control and Prevention [e-pub ahead of print]. JAMA https://doi.org/10.1001/jama.2020.2648. Accessed March 25, 2020.Google Scholar Here, we report that respiratory allergy and controlled allergen exposures are each associated with significant reductions in ACE2 expression. ACE2 expression was lowest in those with both high levels of allergic sensitization and asthma. Importantly, nonatopic asthma was not associated with reduced ACE2 expression. Given that ACE2 serves as the receptor for SARS-CoV-2, our findings suggest a potential mechanism of reduced COVID-19 severity in patients with respiratory allergies. However, it is likely that additional factors beyond ACE2 expression modulate the response to COVID-19 in individuals with allergy, and elucidation of these factors may also provide important insights into COVID-19 disease pathogenesis. The strengths of our study include carefully phenotyped cohorts of children and adults. Further, the allergen challenge studies included both upper and lower airway samples, with each demonstrating a consistent impact on ACE2 expression. The limitations include lack of clinical information to directly link ACE2 expression to SARS-CoV-2 infection and illness severity in our study populations. In addition, we do not have data on the ACE2 protein levels to confirm the gene expression data, although previous work suggests a direct association between ACE2 mRNA levels and ACE2 protein levels in the lung.8Li Y. Zeng Z. Cao Y. Liu Y. Ping F. Liang M. et al.Angiotensin-converting enzyme 2 prevents lipopolysaccharide-induced rat acute lung injury via suppressing the ERK1/2 and NF-kappaB signaling pathways.Sci Rep. 2016; 6: 27911Crossref PubMed Scopus (134) Google Scholar It is important to note that early data in the United States suggest a higher rate of asthma in patients hospitalized for severe COVID-19 illness, but the data do not specify whether the asthma was allergic, which is an important differentiation that relates to our findings. Nor do the data mention the potential presence of other comorbidities, such as obesity, that have been identified as risk factors for COVID-19 illness.9Centers for Disease Control and PreventionHospitalization rates and characteristics of patients hospitalized with laboratory-confirmed coronavirus disease 2019 - COVID-NET, 14 States, March 1-30, 2020.MMWR. 2020; 69Google Scholar Future studies focused on respiratory allergy, asthma, and perhaps other allergic disorders are needed to provide greater understanding of the impact of underlying allergy on COVID-19 susceptibility and disease severity. The modulation of ACE2 expression by type 2 inflammatory processes suggests the need to comprehensively evaluate the role of type 2 immune regulation in COVID-19 pathogenesis. Further elucidation of these relationships could identify novel therapeutic strategies to more effectively control this pandemic. In all 3 studies, total RNA was extracted from epithelial brush samples preserved in RLT buffer (Qiagen, Germantown, Md). Samples were thawed, vortexed, and quick-spun, after which the supernatant was transferred to fresh tubes. The samples were then spun through a Qiashredder column (Qiagen) and extracted by using RNeasy mini kits (Qiagen) with 25-μL elution volumes following the manufacturer's protocol. In the cat allergy upper airway challenge study, sequencing libraries were constructed from total RNA by using TruSeq RNA Sample Preparation Kits v2 (Illumina). In the URECA and adult asthma studies, sequencing libraries were constructed from total RNA by using SMART-Seq v4 Ultra Low Input RNA Kit (Takara). For each study, libraries were clustered onto a flow cell by using a cBOT amplification system with a HiSeq SR v4 Cluster Kit (Illumina). Single-read sequencing was carried out on a HiSeq2500 sequencer (Illumina), using a HiSeq SBS v4 Kit to generate 58-bp reads, with a target of approximately 10 million reads per sample. Samples for each study were processed and sequenced independently. Reads were processed by using workflows managed on the Galaxy platform. Reads were trimmed by 1 base at the 3' end and then trimmed from both ends until base calls had a minimum quality score of at least 30 (Galaxy FASTQ Trimmer tool [version 1.0.0]). FastqMcf (version 1.1.2) was used to remove any remaining adapter sequence. To align the trimmed reads, we used the STAR aligner with the GRCh38 reference genome and gene annotations from ensembl, release 91. Gene counts were generated by using HTSeq-count (version 0.4.1). For quality control, the samples kept were those that had counts greater than 1 million, more than 80% of reads aligned, and a median coefficient of variation (CV) coverage less than 1. Genes were filtered to include those that had a trimmed mean of M-value (TMM) normalization count of at least 1 in at least 10% of libraries and were classified as protein coding by using BioMart.E1Smedley D. Haider S. Durinck S. Pandini L. Provero P. Allen J. et al.The BioMart community portal: an innovative alternative to large, centralized data repositories.Nucleic Acids Res. 2015; 43: W589-W598Crossref PubMed Scopus (531) Google Scholar Counts were transformed to log2 counts per million along with observation-level weights by using voomWithQualityWeights from the limma R packageE2Liu R. Holik A.Z. Su S. Jansz N. Chen K. Leong H.S. et al.Why weight? Modelling sample and observational level variability improves power in RNA-seq analyses.Nucleic Acids Res. 2015; 43: e97Crossref PubMed Scopus (304) Google Scholar to create a weighted gene expression matrix suitable for downstream analyses. Differential expression of ACE2 was assessed independently in each data set by using a weighted linear mixed effects model (limma) appropriate for RNA sequencing data and an empiric Bayes method.E2Liu R. Holik A.Z. Su S. Jansz N. Chen K. Leong H.S. et al.Why weight? Modelling sample and observational level variability improves power in RNA-seq analyses.Nucleic Acids Res. 2015; 43: e97Crossref PubMed Scopus (304) Google Scholar, E3Ritchie M.E. Phipson B. Wu D. Hu Y. Law C.W. Shi W. et al.Limma powers differential expression analyses for RNA-sequencing and microarray studies.Nucleic Acids Res. 2015; 43: e47Crossref PubMed Scopus (20650) Google Scholar Mixed effects linear regression models were used; they included relevant clinical or technical variables (for URECA, cytologically determined cell percentages in the brush and the clinical site were used; for the upper airway challenge study, processing batch was used; and for the adult asthma study, no fixed effects were included) and a random effect of participant in both of the airway challenge studies. P values less than .05 were considered statistically significant. We searched the National Center for Biotechnology Information Gene Expression Omnibus for the terms IL13 and epithelial subset to the organism Homo sapiens.E4Barrett T. Wilhite S.E. Ledoux P. Evangelista C. Kim I.F. Tomashevsky M. et al.NCBI GEO: archive for functional genomics data sets--update.Nucleic Acids Res. 2013; 41: D991-D995Crossref PubMed Scopus (6441) Google Scholar From this we identified 2 studies investigating the effects of IL-13 stimulation on human airway epithelial cells grown at an air-liquid interface that had repeated measures in the IL-13 stimulation and unstimulated groups. In the GSE110799 study, human nasal epithelial cells isolated from nasal turbinates were cultured in an air-liquid interface until the full differentiation was complete, after which the differentiated cells at ALI-D47 were incubated with 100 ng/mL of IL-13 for 3 days. In GSE37693, RNA was isolated from primary culture airway epithelial cells grown at an air-liquid interface treated with or without IL-13 for 21 days.E5Alevy Y.G. Patel A.C. Romero A.G. Patel D.A. Tucker J. Roswit W.T. et al.IL-13-induced airway mucus production is attenuated by MAPK13 inhibition.J Clin Invest. 2012; 122: 4555-4568Crossref PubMed Scopus (152) Google Scholar Differential expression analysis was performed by using GEO2R, which performs voom and limmaE2Liu R. Holik A.Z. Su S. Jansz N. Chen K. Leong H.S. et al.Why weight? Modelling sample and observational level variability improves power in RNA-seq analyses.Nucleic Acids Res. 2015; 43: e97Crossref PubMed Scopus (304) Google Scholar, E3Ritchie M.E. Phipson B. Wu D. Hu Y. Law C.W. Shi W. et al.Limma powers differential expression analyses for RNA-sequencing and microarray studies.Nucleic Acids Res. 2015; 43: e47Crossref PubMed Scopus (20650) Google Scholar in the National Center for Biotechnology Information's Gene Expression Omnibus browser.Table E1Association of T2 biomarkers and nasal brush ACE2 expression in the URECA cohortBiomarkerAssociation with ACE2 expression (β coefficient)P valuePositive allergen-specific IgE–0.0893.1E–5Total IgE level–0.315.1E–6Fractional exhaled nitric oxide–0.453.4E–3Blood eosinophils–0.13.07Nasal epithelial IL-13 expression–0.1238.6E–5 Open table in a new tab COVID-19: Start with the noseJournal of Allergy and Clinical ImmunologyVol. 146Issue 5PreviewWe read with interest the data from Jackson et al1 that revealed lower angiotensin-converting enzyme 2 (ACE2) expression in nasal brushings from children with allergic sensitization, along with a progressive decline in ACE2 expression in relation to increasing IgE sensitization in those with asthma. Moreover, in nasal brushings from adults with allergic rhinitis, ACE2 expression was lower after exposure to cat allergen. Full-Text PDF