Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2021 Jul 14;12(1):4314.
doi: 10.1038/s41467-021-24467-0.

Chronic lung diseases are associated with gene expression programs favoring SARS-CoV-2 entry and severity

Linh T Bui #  1 , Nichelle I Winters #  2 , Mei-I Chung  1 , Chitra Joseph  3 , Austin J Gutierrez  1 , Arun C Habermann  2 , Taylor S Adams  4 , Jonas C Schupp  4 , Sergio Poli  5 , Lance M Peter  1 , Chase J Taylor  2 , Jessica B Blackburn  2 , Bradley W Richmond  2   6 , Andrew G Nicholson  7   8 , Doris Rassl  9 , William A Wallace  10   11 , Ivan O Rosas  12 , R Gisli Jenkins  3 , Naftali Kaminski  4 , Jonathan A Kropski  2   6   13 , Nicholas E Banovich  14 ; Human Cell Atlas Lung Biological Network
Collaborators, Affiliations

Chronic lung diseases are associated with gene expression programs favoring SARS-CoV-2 entry and severity

Linh T Bui et al. Nat Commun. .

Abstract

Patients with chronic lung disease (CLD) have an increased risk for severe coronavirus disease-19 (COVID-19) and poor outcomes. Here, we analyze the transcriptomes of 611,398 single cells isolated from healthy and CLD lungs to identify molecular characteristics of lung cells that may account for worse COVID-19 outcomes in patients with chronic lung diseases. We observe a similar cellular distribution and relative expression of SARS-CoV-2 entry factors in control and CLD lungs. CLD AT2 cells express higher levels of genes linked directly to the efficiency of viral replication and the innate immune response. Additionally, we identify basal differences in inflammatory gene expression programs that highlight how CLD alters the inflammatory microenvironment encountered upon viral exposure to the peripheral lung. Our study indicates that CLD is accompanied by changes in cell-type-specific gene expression programs that prime the lung epithelium for and influence the innate and adaptive immune responses to SARS-CoV-2 infection.

PubMed Disclaimer

Conflict of interest statement

J.A.K. has received advisory board fees from Boehringer Ingelheim, Inc, Janssen Pharmaceuticals, is on the scientific advisory board of APIE Therapeutics, and has research contracts with Genentech. In the last 36 months, N.K. reported personal fees from Biogen Idec, Boehringer Ingelheim, Third Rock, Samumed, Numedii, AstraZeneca, Life Max, Teravance, RohBar, and Pliant and Equity in Pliant; collaboration with MiRagen, AstraZeneca; Grant from Veracyte, all outside the submitted work. In addition, N.K. has a patent for New Therapies in Pulmonary Fibrosis, and Peripheral Blood Gene Expression licensed to Biotech. A.G.N. has received advisory board fees from Boehringer Ingelheim, Galapagos, Medical Quantitative Image Analysis and personal fees for educational material from Up to Date and Boehringer Ingelheim. RGJ reports grants from AstraZeneca, grants from Biogen, personal fees from Boehringer Ingelheim, personal fees from Chiesi, personal fees Daewoong, personal fees from Galapagos, grants from Galecto, grants from GlaxoSmithKline, personal fees from Heptares, nonfinancial support from NuMedii, grants and personal fees from Pliant, personal fees from Promedior, non-financial support from Redx, personal fees from Roche, other from Action for Pulmonary Fibrosis, outside the submitted work. All other authors declare no competing interests.

Figures

Fig. 1
Fig. 1. Percentage of cells expressing SARS-CoV-2 receptor genes in lung cell types in different diagnosis subgroups.
a Percentage of cells expressing ACE2 and TMPRSS2 in all cell types. Numbers are the total number of ACE2 + or TMPRSS2 + cells in each cell type in the dataset. b, c Percentage of cells expressing ACE2 (b) and TMPRSS2 (c) in each diagnosis group in the epithelial cell types. d Venn diagram shows overlapping of cells co-expressing the proposed receptors (ACE2, BSG and NRP1) and the protease TMPRSS2. e, f Percentage of cells co-expressing receptors and TMPRSS2 split by cell type and diagnosis group. Plots were generated with mean values of percentage of cells per individual samples, and data are presented as mean values ± SEM. Significant differences between diagnosis groups were calculated using Tukey_HSD test, p value < 0.05: *p value < 0.01: **p value < 0.001: ***p value < 0.0001: ****.
Fig. 2
Fig. 2. Expression profile of SARS-CoV-2 mediators and response genes in the epithelial cell population.
a Binary heatmap representing a manually curated list of genes associated with SARS-CoV-2. Orange elements indicate genes with increased expression and white elements indicate genes with decreased expression in CLD samples; Not detected: gene expression was not detected in either of the two tested populations (CLD vs. Control). Differential expressed genes (DEGs) between CLD and control samples (FDR ≤ 0.1) are outlined in black. b Violin plot depicts gene expression level in CLD and control of the two SARS-CoV-2 proteases TMPRSS2 and CTSL. c SARS-CoV-2 entry module score in different cell types, SARS-CoV-2 mediators included ACE2, BSG (CD147), NPR1, HSPA5 (GRP78), TMPRSS2, CTSL, ADAM17, FURIN. The outliers were removed in this plot, please see Supplementary Fig. 8a with outliers included. Boxes: interquartile range, lower and upper hinges correspond to the first and third quantiles, upper and lower whisker extends from the hinge to the largest values or smallest values of 1.5 x interquartile range; *p value < 0.05, **p value < 0.01, ***p value < 0.001, ****p-value < 0.0001, Tukey_HSD post-hoc test. ACE2 and ITGB6 protein expression in IPF lung sections. IPF lung sections stained for ACE2: d small airway, e large airway and f lung parenchyma. IPF lung sections stained for αvβ6: g small airway, h large airway and i lung parenchyma. j Semi-quantitative evaluation of ACE2 scoring among control (n = 12 for each tissue) and IPF (n = 62 for each tissue) sections (both the percentage of staining and staining intensity of ACE2 expression; 0-Negative; 1–0–⩽10%; 2-11–⩽25%; 3-⩽26%), data are presented as mean values ± SEM. Significant differences between IPF and control were calculated using Tukey HSD test, p value < 0.05 *. Scale bar = 100 µm. For d–j: a total of 12 normal lung samples and 62 IPF samples were used.
Fig. 3
Fig. 3. CLD AT2 cells exhibit baseline differences in gene expression profile coping with viral infection.
a Significant gene expression correlation in AT2 cells between TMPRSS2 and ACE2, BSG (CD147) and NPR1 in COPD and IPF samples, each dot represents the average expression level of the genes of interest per sample, pairwise gene correlation analysis was done using a fitting linear model and p value was calculated using Anova. b Boxplot shows differences in gene expression of selected SARS-CoV-2 response genes in the AT2 cell types among different diagnosis groups, Boxes: interquartile range, lower and upper hinges correspond to the first and third quantiles, upper and lower whisker extends from the hinge to the largest values or smallest values of 1.5 × interquartile range; **p value-adj ≤ 0.05 (negative binomial test, corrected for Age, Ethnicity, Smoking_status and Dataset). c Upset plot shows shared differential expression genes (DEGs) between different comparisons: ACE2− CLD vs. Control, ACE2 + CLD vs. Control, CLD ACE2 + vs. ACE2-, Control ACE2 + vs. ACE2- and ACE2 correlated genes in the AT2 cells. d Upregulation of two genes uniquely differentially expressed in the CLD ACE2 + vs. ACE2−. e Spearman gene correlation analysis identified genes correlated with ACE2 expression in AT2 ACE2 + cells in different diagnosis groups, p-value was adjusted using Benjamini-Hochberg corrections, dashed lines indicate the 99th percentile of Spearman rho values.
Fig. 4
Fig. 4. Analysis of SARS-CoV-2 candidate immune response genes in immune cells.
a Quantification of cell types as a percent of all immune cells in control and diseased lungs, numbers represent the total numbers of the cell type per individual samples, data are presented as mean values ± SEM. b Binary heatmap representing a manually curated list of genes associated with SARS-CoV-2. Orange elements indicate genes with increased expression and white elements indicate genes with decreased expression; Not detected: gene expression was not detected in either of the two tested populations (CLD vs. Control); DEGs with FDR ≤ 0.1 are outlined in black. c Differential expression analysis for SARS-CoV-2 immune candidate genes in cDCs, Macrophages and Monocytes. *p-adjusted value < 0.1, **p-adjusted value < 0.05, p-adjusted value was Bonferroni adjusted from Seurat FindMarkers differential expression analysis using a negative binomial test and corrected for Age, Ethnicity, Smoking_status and Dataset. d Compared to the healthy control samples, HLA type II score is higher in all disease groups (especially Other-ILD). In the T cell population, cytotoxicity scores (e) and exhaustion scores (f) are higher in the disease samples than in control samples. In a, d, e, and f: Boxes: interquartile range, lower and upper hinges correspond to the first and third quantiles, upper and lower whisker extends from the hinge to the largest values or smallest values of 1.5 x interquartile range; Tukey_HSD post-hoc test: *p value < 0.05, **p value < 0.01, ***p value < 0.001, ****p value < 0.0001. See Supplementary Fig. 10 for plots with outliers included for d–f.
Fig. 5
Fig. 5. Model of alterations in the diseased lung related to SARS-CoV2 pathogenesis.
(1) In the IPF lung, there is a proximalization of the distal airway. ACE2 + epithelial cells cluster in the small airways though total ACE2 + cell numbers are similar to control. (2) The viral entry score (accounting for all described putative receptors and proteases) is increased in diseased lungs. (3) Diseased epithelial cells have alterations in key SARS-CoV-2 response genes/pathways. (4) In the CLD lung, there is increased expression of cytotoxicity and exhaustion genes in immune cell populations and alterations in viral response pathways (interferon, antigen presentation). Figure created in Biorender.com.
See this image and copyright information in PMC

Update of

Similar articles

  • Chronic lung diseases are associated with gene expression programs favoring SARS-CoV-2 entry and severity.
    Bui LT, Winters NI, Chung MI, Joseph C, Gutierrez AJ, Habermann AC, Adams TS, Schupp JC, Poli S, Peter LM, Taylor CJ, Blackburn JB, Richmond BW, Nicholson AG, Rassl D, Wallace WA, Rosas IO, Jenkins RG, Kaminski N, Kropski JA, Banovich NE; Human Cell Atlas Lung Biological Network. Bui LT, et al. bioRxiv [Preprint]. 2021 Jan 4:2020.10.20.347187. doi: 10.1101/2020.10.20.347187. bioRxiv. 2021. Update in: Nat Commun. 2021 Jul 14;12(1):4314. doi: 10.1038/s41467-021-24467-0. PMID: 33106805 Free PMC article. Updated. Preprint.
  • Dysbalance of ACE2 levels - a possible cause for severe COVID-19 outcome in COPD.
    Fließer E, Birnhuber A, Marsh LM, Gschwandtner E, Klepetko W, Olschewski H, Kwapiszewska G. Fließer E, et al. J Pathol Clin Res. 2021 Sep;7(5):446-458. doi: 10.1002/cjp2.224. Epub 2021 May 12. J Pathol Clin Res. 2021. PMID: 33978304 Free PMC article.
  • Single-cell meta-analysis of SARS-CoV-2 entry genes across tissues and demographics.
    Muus C, Luecken MD, Eraslan G, Sikkema L, Waghray A, Heimberg G, Kobayashi Y, Vaishnav ED, Subramanian A, Smillie C, Jagadeesh KA, Duong ET, Fiskin E, Torlai Triglia E, Ansari M, Cai P, Lin B, Buchanan J, Chen S, Shu J, Haber AL, Chung H, Montoro DT, Adams T, Aliee H, Allon SJ, Andrusivova Z, Angelidis I, Ashenberg O, Bassler K, Bécavin C, Benhar I, Bergenstråhle J, Bergenstråhle L, Bolt L, Braun E, Bui LT, Callori S, Chaffin M, Chichelnitskiy E, Chiou J, Conlon TM, Cuoco MS, Cuomo ASE, Deprez M, Duclos G, Fine D, Fischer DS, Ghazanfar S, Gillich A, Giotti B, Gould J, Guo M, Gutierrez AJ, Habermann AC, Harvey T, He P, Hou X, Hu L, Hu Y, Jaiswal A, Ji L, Jiang P, Kapellos TS, Kuo CS, Larsson L, Leney-Greene MA, Lim K, Litviňuková M, Ludwig LS, Lukassen S, Luo W, Maatz H, Madissoon E, Mamanova L, Manakongtreecheep K, Leroy S, Mayr CH, Mbano IM, McAdams AM, Nabhan AN, Nyquist SK, Penland L, Poirion OB, Poli S, Qi C, Queen R, Reichart D, Rosas I, Schupp JC, Shea CV, Shi X, Sinha R, Sit RV, Slowikowski K, Slyper M, Smith NP, Sountoulidis A, Strunz M, Sullivan TB, Sun D, Talavera-López C, Tan P, Tantivit J, Travaglini KJ, Tucker NR, Vernon KA, Wadsworth MH, Waldman J, Wang X, Xu K, Yan … See abstract for full author list ➔ Muus C, et al. Nat Med. 2021 Mar;27(3):546-559. doi: 10.1038/s41591-020-01227-z. Epub 2021 Mar 2. Nat Med. 2021. PMID: 33654293 Free PMC article.
  • Classical and alternative receptors for SARS-CoV-2 therapeutic strategy.
    Masre SF, Jufri NF, Ibrahim FW, Abdul Raub SH. Masre SF, et al. Rev Med Virol. 2021 Sep;31(5):1-9. doi: 10.1002/rmv.2207. Epub 2020 Dec 26. Rev Med Virol. 2021. PMID: 33368788 Free PMC article. Review.
  • Modulation of Host Immune Response Is an Alternative Strategy to Combat SARS-CoV-2 Pathogenesis.
    Singh L, Bajaj S, Gadewar M, Verma N, Ansari MN, Saeedan AS, Kaithwas G, Singh M. Singh L, et al. Front Immunol. 2021 Jul 8;12:660632. doi: 10.3389/fimmu.2021.660632. eCollection 2021. Front Immunol. 2021. PMID: 34305892 Free PMC article. Review.
See all similar articles

Cited by

See all "Cited by" articles

References

    1. Huang C, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395:497–506. doi: 10.1016/S0140-6736(20)30183-5. - DOI - PMC - PubMed
    1. Zhu N, et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N. Engl. J. Med. 2020;382:727–733. doi: 10.1056/NEJMoa2001017. - DOI - PMC - PubMed
    1. Wang C, Horby PW, Hayden FG, Gao GF. A novel coronavirus outbreak of global health concern. Lancet. 2020;395:470–473. doi: 10.1016/S0140-6736(20)30185-9. - DOI - PMC - PubMed
    1. Zhou P, et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020;579:270–273. doi: 10.1038/s41586-020-2012-7. - DOI - PMC - PubMed
    1. Home - Johns Hopkins Coronavirus resource center. https://coronavirus.jhu.edu/.

Publication types

MeSH terms

Substances

-