Next Article in Journal
Plasminogen Activator Inhibitor-1 Secretion by Autophagy Contributes to Melanoma Resistance to Chemotherapy through Tumor Microenvironment Modulation
Next Article in Special Issue
Baseline Characteristics and Outcomes of Cancer Patients Infected with SARS-CoV-2 in the Lombardy Region, Italy (AIOM-L CORONA): A Multicenter, Observational, Ambispective, Cohort Study
Previous Article in Journal
Breast Tumor Characterization Using [18F]FDG-PET/CT Imaging Combined with Data Preprocessing and Radiomics
Previous Article in Special Issue
SARS-CoV-2 Serology Monitoring of a Cancer Center Staff in the Pandemic Most Infected Italian Region
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 6
10.3390/cancers13061250
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single-Cell Expression Landscape of SARS-CoV-2 Receptor ACE2 and Host Proteases in Normal and Malignant Lung Tissues from Pulmonary Adenocarcinoma Patients

by
Guangchun Han
1,†,
Ansam Sinjab
2,†,
Kieko Hara
2,
Warapen Treekitkarnmongkol
2,
Patrick Brennan
3,
Kyle Chang
4,
Elena Bogatenkova
3,
Beatriz Sanchez-Espiridion
2,
Carmen Behrens
5,
Luisa M. Solis
2,
Boning Gao
6,
Luc Girard
6,
Jianjun Zhang
5,
Boris Sepesi
7,
Tina Cascone
5,
Lauren A. Byers
5,
Don L. Gibbons
5,
Jichao Chen
8,
Seyed Javad Moghaddam
8,
Edwin J. Ostrin
9,
Paul Scheet
10,
Junya Fujimoto
2,
Jerry Shay
11,
John V. Heymach
5,
John D. Minna
6,
Steven Dubinett
12,
Ignacio I. Wistuba
2,
Christopher S. Stevenson
13,
Avrum E. Spira
13,14,
Linghua Wang
1,* and
Humam Kadara
2,*add Show full author list remove Hide full author list
1
Department of Genomic Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
2
Department of Translational Molecular Pathology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
3
Pathology Department, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
4
Guardant Health, Redwood City, CA 94063, USA
5
Department of Thoracic, Head and Neck Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
6
Hamon Center for Therapeutic Oncology Research, University of Texas Southwestern, Dallas, TX 75390, USA
7
Department of Cardiovascular and Thoracic Surgery, The University of Texas MD Anderson Cancer Center, Houston, TX 77005, USA
8
Department of Pulmonary Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
9
Department of General Internal Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
10
Department of Epidemiology, The University of Texas MD Anderson Cancer Center, Houston, TX 77230, USA
11
Department of Cell Biology, University of Texas Southwestern, Dallas, TX 75390, USA
12
Department of Medicine, The University of California Los Angeles, Los Angeles, CA 90095, USA
13
Lung Cancer Initiative at Johnson and Johnson, Cambridge, MA 02142, USA
14
Section of Computational Biomedicine, Boston University, Boston, MA 02215, USA
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this paper.
Cancers 2021, 13(6), 1250; https://doi.org/10.3390/cancers13061250
Submission received: 13 January 2021 / Revised: 7 March 2021 / Accepted: 9 March 2021 / Published: 12 March 2021
(This article belongs to the Collection The Impact of COVID-19 Infection in Cancer)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The coronavirus disease 2019 (COVID-19) pandemic continues to spread rapidly on a global scale. When presenting with severe respiratory complications, COVID-19 results in markedly high death rates, particularly among patients with comorbidities such as cancer. Motivated by the ongoing global health crisis, we leveraged a growing in-house cohort of pulmonary tissues from lung cancer patients to analyze, at high resolution, the expression of host proteins implicated in the entryway of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into lung epithelial cells. Our results identify key pathways in lung pathobiology and inflammation that offer the potential to identify novel markers and therapeutic targets that can be repurposed for clinical management of COVID-19, particularly among lung cancer patients, a population that represents over half a million individuals in the United States alone.

Abstract

The novel coronavirus SARS-CoV-2 is the causative agent of the COVID-19 pandemic. Severely symptomatic COVID-19 is associated with lung inflammation, pneumonia, and respiratory failure, thereby raising concerns of elevated risk of COVID-19-associated mortality among lung cancer patients. Angiotensin-converting enzyme 2 (ACE2) is the major receptor for SARS-CoV-2 entry into lung cells. The single-cell expression landscape of ACE2 and other SARS-CoV-2-related genes in pulmonary tissues of lung cancer patients remains unknown. We sought to delineate single-cell expression profiles of ACE2 and other SARS-CoV-2-related genes in pulmonary tissues of lung adenocarcinoma (LUAD) patients. We examined the expression levels and cellular distribution of ACE2 and SARS-CoV-2-priming proteases TMPRSS2 and TMPRSS4 in 5 LUADs and 14 matched normal tissues by single-cell RNA-sequencing (scRNA-seq) analysis. scRNA-seq of 186,916 cells revealed epithelial-specific expression of ACE2, TMPRSS2, and TMPRSS4. Analysis of 70,030 LUAD- and normal-derived epithelial cells showed that ACE2 levels were highest in normal alveolar type 2 (AT2) cells and that TMPRSS2 was expressed in 65% of normal AT2 cells. Conversely, the expression of TMPRSS4 was highest and most frequently detected (75%) in lung cells with malignant features. ACE2-positive cells co-expressed genes implicated in lung pathobiology, including COPD-associated HHIP, and the scavengers CD36 and DMBT1. Notably, the viral scavenger DMBT1 was significantly positively correlated with ACE2 expression in AT2 cells. We describe normal and tumor lung epithelial populations that express SARS-CoV-2 receptor and proteases, as well as major host defense genes, thus comprising potential treatment targets for COVID-19 particularly among lung cancer patients.

1. Introduction

In late December 2019, an outbreak of lung pneumonia initially of unknown cause was reported in China [1]. This emerging disease, termed coronavirus disease 2019 (COVID-19), was soon thereafter attributed to infection with the novel zoonotically-transmitted coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [2]. On 11 March 2020, COVID-19 was declared a rapidly spreading global pandemic by the World Health Organization. By March 2021, SARS-CoV-2 had infected more than 118 million people worldwide, including 28 million cases and 520,000 deaths in the United States.
The clinical presentation of COVID-19 is diverse, ranging from asymptomatic infection and mild upper respiratory illness to pneumonia, acute respiratory distress syndrome (ARDS), respiratory failure, and death [1]. Recent clinical reports have suggested that old age and comorbidities such as cardiovascular disease, COPD, and cancer are risk factors for COVID-19-associated severe pneumonia and death [1]. Notably, lung cancer was found to correlate with adverse outcomes in patients with COVID-19 [3]. This has raised key questions in the clinical management of patients with both COVID-19 and lung malignancy and warrants a deeper knowledge of the yet unknown pathological mechanisms and effects of SARS-CoV-2 infection in patients with lung cancer.
Recent studies have shown that SARS-CoV-2 infects airway cells by binding its spike protein to angiotensin-converting enzyme 2 encoded by gene ACE2, the same receptor used by SARS-CoV [4,5,6,7]. ACE2 had been shown to mediate important roles in lung function, including protection from lung injury caused by SARS-CoV infection [8] and inhibition of angiogenesis in lung cancer [9]. Despite these insights, and the supposable heightened risk of lung cancer patients for COVID-19-associated mortality, the expression of SARS-CoV-2-related genes in lung tumor and uninvolved tissue is still poorly understood.
To fill these voids, we leveraged our ongoing efforts in a single-cell transcriptomic analysis of 186,916 cells, including a large number of epithelial cells (n = 70,030) derived from 5 lung adenocarcinomas (LUADs) and 14 matching uninvolved normal lung tissues. We show epithelial-specific expression patterns for ACE2 as well as TMPRSS2 and TMPRSS4, serine proteases with emerging roles in SARS-CoV-2 priming [7,10]. Among all lung epithelial subsets, we find the highest expression of ACE2 in alveolar type 2 (AT2) and malignant-enriched subsets and of TMPRSS2 and TMPRSS4 in AT2 and malignant-enriched cell populations, respectively. ACE2-positive AT2 cells also expressed genes with important yet understudied roles in lung pathobiology. Our study provides a comprehensive overview of the expression of ACE2, SARS-CoV-2-priming proteases, as well as host defense and scavenging genes in the malignant lung and nearby epithelium that may constitute targets that can be repurposed for the clinical management of COVID-19 in LUAD patients.
Some of the results of this study have been previously reported in the form of an online preprint [11].

2. Methods

Additional descriptions of methods can be found in the Online Data Supplement.

2.1. Lung Tissue Processing and Single-Cell RNA-Sequencing (scRNA-Seq)

Five patients undergoing surgical resection for primary early-stage LUAD (I-IIIA) were carefully selected for the derivation of two to three uninvolved normal lung samples for single-cell analysis (n = 19 samples, Table S1). Samples were obtained from banked or residual specimens from patients evaluated at the University of Texas MD Anderson Cancer Center. Following tissue digestion and red blood cell removal, cells were sorted (by fluorescent-activated cell sorting) into viable singlets and, in samples from Patients 2 to 5, into separate viable epithelial (EPCAM+) and nonepithelial (EPCAM−) fractions. Single-cell gene expression libraries were generated from 35 sorted fractions using the 10× Genomics platform (Pleasanton, CA, USA) and sequenced on the Illumina NovaSeq 6000 platform (San Diego, CA, USA; Online Data Supplement).

2.2. scRNA-Seq Data Analysis

Raw scRNA-seq data were preprocessed, demultiplexed, and aligned to human GRCh38 reference and feature-barcodes generated using CellRanger (10× Genomics, Pleasanton, CA, USA; version 3.0.2). Details of quality control, including quality check, data filtering, identification and removal of cellular debris, doublets and multiplets, and batch effect evaluation and correction, are found in the Online Data Supplement. Following quality filtering, a total of 186,916 cells were retained for downstream analysis. Raw unique molecular identifier (UMI) counts were log normalized and used for principal component analysis. We applied Seurat [12] for unsupervised clustering analysis and Uniform Manifold Approximation and Projection (UMAP) [13] for dimensionality reduction and visualization. Lung and airway subcluster lineage (e.g., of EPCAM+ and AT2) was defined based on the enrichment of canonical marker genes as well as top-ranked differentially expressed genes (DEGs) in each cluster using the FindAllMarkers function in the Seurat R package. We also applied the single-cell consensus clustering (SC3) approach [14] for unsupervised clustering analysis, with default parameters independent of cell lineage annotation.

2.3. Statistical Analysis

All statistical analyses were performed using R package version 3.6.0. DEG analysis (e.g., between ACE2-positive and ACE2-negative cells) was calculated using the FindMarkers function in R. Pseudobulk gene expression values for defined cell clusters were calculated by computing the mean expression of each gene across all cells in a specific cluster. Pearson’s correlation analysis was used to identify genes significantly correlated with ACE2 expression or with an AT2 meta-score. All statistical significance testing was two-sided, and results were considered statistically significant at p-values < 0.05. The Benjamini–Hochberg method was applied to control the false discovery rate (FDR) in multiple comparisons (e.g., DEG analysis) and to calculate adjusted p-values (q-values).

3. Results

3.1. Single-Cell Decoding of ACE2 Expression and SARS-CoV-2 Priming Proteases in Lung Tissues of LUAD Patients

We performed a single-cell analysis of normal lung tissues and matched early-stage LUADs from five patients using droplet-based scRNA-seq. In the first patient, we attained a limited fraction of epithelial (EPCAM+) cells by an unbiased analysis of two normal lung tissues and one LUAD (~4%, n = 624 cells), in line with studies of other organs [15]. To better capture lung epithelial cells, we performed sequencing of cells with enriched (by sorting for EPCAM) epithelial subsets from three normal lung tissues and a matched LUAD each from four additional patients. In total, 186,916 cells, 70,030 of which were epithelial, from the 5 LUADs and 14 uninvolved normal lung tissues were retained for analysis following quality control. Based on canonical gene expression markers, cells clustered into distinct epithelial, endothelial, myeloid, lymphoid, or stromal cell clusters (Figure 1A). Prompted by the ongoing COVID-19 pandemic caused by infection with the novel coronavirus SARS-CoV-2, we leveraged our unique LUAD and normal lung tissue scRNA-seq dataset to interrogate at single-cell resolution lung expression patterns of the SARS-CoV-2 receptor ACE2, as well as TMPRSS2 and TMPRSS4, two related membrane-bound serine proteases recently shown to be crucial for SARS-CoV-2 spike protein priming upon entry [7,10]. We found that all three SARS-CoV-2-related genes were nearly restricted to the epithelial compartment (Figure 1B), including when considering only cycling cells (Figure S1). These findings suggest that SARS-CoV-2 receptor and priming proteases are restricted to epithelial cells in the ecosystem of the normal and malignant lung.

3.2. Expression Patterns of ACE2, TMPRSS2, and TMPRSS4 in Normal and Malignant-Enriched Epithelial Subsets

We then probed expression patterns of the three genes among the 70,030 lung epithelial cells (5 tumor samples, n = 13,098; 14 normal lung tissues, n = 56,932) (Figure 2A, Table S2). We found on average 9307 UMIs and 2616 genes per epithelial cell (Table S2). Clustering analysis identified multiple airway lineages, including alveolar type 1 (AT1; n = 10,775), AT2 (n = 27,235), club and secretory (n = 4624), ciliated (n = 3247), as well as basal (n = 5119), including proliferative (n = 447) cells (Figure 2A and Figure S2A, Table S2). We also identified rare or transitory subpopulations such as bronchioalveolar cells, alveolar progenitors [16,17], and the CFTR-expressing pulmonary ionocytes, which represent a novel lung epithelial cell type that was recently identified in mouse and human airways [18,19,20] (Figure 2A). Additionally, a malignant-enriched cluster that was mostly derived from LUAD specimens clustered distinctly from cells of uninvolved lung tissues and exhibited high expression of tumor markers (e.g., CEACAM5), as well as comprised mixed-lineage genes in line with previous studies [21] (Figure S2A,B). We found that the fraction of ACE2-expressing cells among all lung epithelial cells was low (n = 1242, 1.8%, Table S3), albeit being the highest in AT2, club/secretory, and malignant subsets (Figure 2A). In comparison, TMPRSS2 displayed a more ubiquitous epithelial expression pattern, whereby its fraction was the highest in alveolar subsets (Figure 2B), and TMPRSS4 was largely abundant in cells of the malignant-enriched cluster (Figure 2C). Further quantification revealed that the highest fractions of ACE2-expressing cells were found in the malignant-enriched (3.5%), AT2 (2.2%), and club/secretory (2.4%) cell clusters (Figure 2D). Among those clusters with an absolute number of cells >50 and >2% ACE2-positive cells, AT2 cells exhibited the highest expression of the SARS-CoV-2 receptor (Figure 2D). The frequency of TMPRSS2-expressing epithelial cells was the highest among alveolar subclusters, including AT2 (62.5%), AT1 (57.8%), and alveolar progenitor (37.7%) cells, while cells of the malignant-enriched cluster displayed lower levels of TMPRSS2 (Figure 2E). Interestingly, the other SARS-CoV-2-priming protease, TMPRSS4, was distinctively expressed in cells of the malignant-enriched cluster and with the greatest abundance (74.6%) and highest expression levels compared to other epithelial subclusters (Figure 2F). We further interrogated tumor-associated expression patterns of these SARS-2-CoV-related host genes using bulk expression data of 52 LUAD and normal lung pairs from the TCGA cohort [22]. In line with our scRNA-seq data, TMPRSS2 was significantly downregulated in LUADs versus matched normal lung tissues, and, in sharp contrast, TMPRSS4 was largely upregulated in the LUADs (both p < 10−7; Figure S2C). Next, we examined ACE2 expression in subclusters of AT2 cells (Figure S3A). We found that the same AT2 subclusters harbored the highest average expression (AT2_c2) as well as the highest fraction of cells (AT2_c3) positive for all ACE2, TMPRSS2, and TMPRSS4 (Figure S3B–D), suggesting co-expression patterns for these SARS-CoV-2 genes in specific AT2 subsets. We also analyzed ACE2 expression in AT2 and malignant-enriched clusters based on whether patients in our single-cell cohort were on antihypertensive treatments (Table S1), since earlier studies have shown that drugs such as losartan may impact ACE2 expression levels [23]. Interestingly, while losartan-treated patients had lower fractions of ACE2-expressing AT2 cells or cells of the malignant-enriched cluster, expression levels in these subsets were significantly higher compared to their counterparts from the untreated patients (p < 3 × 10−35; Figure S4A,B). Taken together, these data describe tumor- and lung-pertinent expression patterns of the major SARS-CoV-2 receptor ACE2 and priming proteases TMPRSS2 and TMPRSS4 in LUAD patients.

3.3. Genes Co-Expressed with ACE2 in Lung Epithelial Cells

We next sought to identify differentially expressed genes (DEGs) in AT2 cells of LUAD patients based on ACE2 expression. Cutoffs of absolute gene expression (fold-change: >1.2) and a FDR (q-value < 0.05) were applied to select DEGs between ACE2-expressing (n = 607) and ACE2-absent (n = 26,628) AT2 cells. Among genes upregulated in ACE2-expressing AT2 cells, we identified genes that are highly pertinent to lung epithelial biology and disease such as HHIP (lung branching and COPD [24,25]), FGG (fibrosis, pneumonia, and inflammation [26,27]), and C4BPA (complement system, pneumonia, and infection [28,29]) (Figure 3A,B). In addition, we found that ACE2-expressing AT2 cells exhibited significantly a higher expression of scavengers such as CD36 [30], as well as DMBT1, a pattern recognition receptor that plays crucial roles in host defense against bacterial and viral pathogens, including influenza A and human immune deficiency virus 1 (HIV-1) [31] (Figure 3A,B). Interestingly, DMBT1 abundance and expression were markedly and distinctly highest in AT2 cells relative to other lung cell subsets (Figure S5A). DMBT1 expression in AT2 cells was also higher compared to each of BSG (also known as CD147), which was recently suggested to represent a potential alternate route of entry for SARS-CoV-2 [32], and other coronavirus receptors such as the MERS-CoV receptor DPP4 [33] (Figure S5A). Notably, both DMBT1 and BSG (albeit to a lesser extent), but not DPP4, exhibited a single-cell expression pattern that resembled that of ACE2 across the different AT2 subclusters (Figure S5B).
Our findings on the significant associations between ACE2 and DMBT1 prompted us to probe the correlation of both genes in the AT2 compartment. We performed a pseudobulk analysis of the AT2 cluster (by sample) in our cohort and found that DMBT1 and ACE2 exhibited a trend of correlation (R = 0.3), albeit not reaching statistical significance (Figure 4A). We also noted nonsignificant positive and negative correlations between ACE2 on the one hand and TMPRSSS2 and TMPRSS4 on the other, respectively (Figure 4A). To further interrogate those findings in larger cohorts, we performed a deconvolution analysis of bulk RNA-seq data from the TCGA LUAD cohort to estimate the abundance of AT2 cells in normal lung tissues (n = 110). ACE2, DMBT1, and TMPRSS2 expression levels significantly and positively correlated with AT2 fractions (indicated by the AT2 meta-score, p < 0.05) (Figure 4B). In contrast, TMPRSS4 expression was negatively correlated with AT2 fractions (Figure 4B), consistent with our observations of its abundance in tumors. We observed the same trends when separately analyzing stage I LUADs or stage II and III LUADs combined (Figure S6). Interestingly, the correlation between AT2 cell fractions and expression of ACE2 and DMBT1 was statistically significant only in stage I LUADs only (p < 0.05 and p = 0.01, respectively) (Figure S6). Further, among TCGA normal lung tissues with high AT2 fractions (meta-score > 15.48; n = 27), DMBT1 expression significantly and positively correlated with ACE2 (R = 0.41, p < 0.05) (Figure 4C). We did not find any significant correlation between ACE2 and TMPRSS2 or TMPRSS4 in the TCGA lung tissues by deconvolution analysis (Figure 4C). Our findings point to specific ACE2-expressing AT2 cells in the lungs of LUAD patients that also co-express other pathogen receptors and scavengers (e.g., DMBT1), thereby possibly representing a minute subpopulation with unique host defense properties and functions.

4. Discussion

We interrogated, by single-cell RNA-sequencing, 186,916 lung cells, including a large number of epithelial cells (n = 70,030), from tumor and multiple matched normal tissues to comprehensively examine abundance and expression patterns of the SARS-CoV-2 receptor ACE2 and the pathogen’s priming proteases in the lungs of LUAD patients. While ACE2 was expressed in a low fraction of lung epithelial cells (roughly 1.8%), its levels among all lung cell subsets were highest in cells of AT2 and in malignant-enriched clusters. We also found, in alveolar and malignant-enriched subsets, the highest fractions of cells expressing SARS-CoV-2-priming proteases TMPRSS2 and TMPRSS4, respectively, as well as increased expression levels of those genes. ACE2-positive AT2 cells co-expressed genes with crucial roles in lung pathological conditions such as COPD, pneumonia, and pathogen infection, including the viral binding scavenger DMBT1. We also found that DMBT1 positively correlated with AT2 cell fractions and with ACE2 itself within the AT2 compartment. Our findings suggest that cells expressing SARS-CoV-2 receptors and coreceptors in malignant lungs and surrounding normal tissues are relatively scarce, and they exhibit unique molecular and biological features that are pertinent to antiviral and host defense functions by the lung epithelium.
Our findings on the expression of ACE2 in AT2 cells are in line with previous studies [34,35]. Hamming et al. showed the abundant immunohistochemical expression of ACE2 protein in AT2 cells in the lung [35]. Ace2 was shown to not only be expressed in AT2 cells in the mouse lung but also to regulate alveolar epithelial cell survival [36]. Of note, Mossel and colleagues demonstrated that SARS-CoV, by binding to ACE2, replicates primarily in AT2 and not AT1 cells [34]. AT2 cells were also shown to be primarily responsible for an innate immune response in the lungs against SARS-CoV infection [37]. Interestingly, a recent report by Zhao and colleagues revealed that ACE2 expression was concentrated in a small population (~1.4%) of AT2 cells [38]. Here, we further interrogated a relatively larger and more diverse repertoire of airway lineages, from both lung tumor (LUAD) tissues and multiple matched normal tissues per patient, and found that not only were expression levels of ACE2 the highest in a minute fraction of AT2 cells (~2.2%, in close agreement with Zhao et al.’s report) but also that the greatest fraction of ACE2-positive cells was in cells of the malignant-enriched cluster. We also identified high expression of the serine proteases TMPRSS2 and TMPRSS4 in cells of alveolar and malignant-enriched clusters, respectively, the former of which also positively correlated with AT2 abundance. More recently, these proteases have been highlighted as additional host proteins mediating SARS-CoV-2 cleavage and internalization [7,10]. Interestingly, TMPRSS2 and TMPRSS4 had been previously ascribed diverse roles in mucosal pathobiology, including that of the lung, such as the activation of influenza and SARS-CoV viruses [39,40,41], promotion of lung fibrosis [42], as well as the induction of lung cancer stem-cell-like properties [43], growth, metastasis, and resistance to chemotherapy [44]. It is worthwhile to mention that patients with cancer, including lung malignancy, were shown to represent a population that is extremely vulnerable to COVID-19 [3,45]. It is plausible to surmise that a better understanding of the expression of SARS-CoV-2 receptors is highly warranted in lung cancer patients, particularly in LUADs, whose pathogenesis is linked to ACE2-expressing airway lineages (e.g., AT2), and also since LUADs arise in the lung periphery, a major site for the development of pneumonia in COVID-19 patients. Our finding on the expression of ACE2 and TMPRSS4 in lung cells enriched with malignant (adenocarcinoma) features may have implications for the management of COVID-19 in lung cancer, and clinical studies are warranted to explore this conjecture. Our analysis also revealed significantly elevated ACE2 levels, among ACE2-positive AT2 and cells of the malignant-enriched cluster, in patients who were on hypertension treatment with losartan compared to those without antihypertensive therapy. These findings are in line with previous reports, suggesting that antihypertensive agents may augment ACE2 expression [23], and further, may partially explain the increased prevalence of hypertension in COVID-19 patients [46]. We thus cannot construe the implications of this finding on COVID-19 management [46], and further studies (e.g., larger cohort) are warranted to better explore the association between treatment with antihypertensive drugs and expression of SARS-CoV-2 receptors in the lungs of LUAD patients.
While our study provides novel insights into the expression of SARS-CoV-2 entry factors in LUAD tissues, its results need to be interpreted with caution. Our patient cohort is limited in size, and our study was designed with a focus on the single-cell analysis of epithelial cells from LUAD specimens. Future studies interrogating larger and more diverse cohorts will provide more power to analyze the expression of SARS-CoV-2 factors by LUAD stage or across histological subtypes, as well as in individual patients and at the single-cell level. Furthermore, it is conceivable that expression patterns of ACE2 and SARS-CoV-2 host proteases in malignancies of the upper airways may be different from the patterns observed in LUAD (i.e., a cancer of the peripheral lung), and this is supported by recent studies showing the expression of ACE2 in secretory cells (such as those in the upper airway) [47,48], which we also noted in our limited cohort. Additional studies in more diverse lung cancer cohorts are thus needed to interrogate these suppositions.
We found that ACE2-positive compared to ACE2-negative AT2 cells exhibited increased levels of genes with crucial roles and expression features in lung pathological diseases. The hedgehog interacting protein HHIP was not only shown to play important roles in airway branching during lung development [24], but also single nucleotide polymorphisms of this gene are associated with increased risk for COPD [25], a pulmonary ailment characterized by chronic inflammation [49]. FGG coding for fibrinogen-gamma was shown to be induced by proinflammatory cytokines [27] and to be elevated in lung pneumonia and infection [26]. C4BPA, coding for C4BP and part of the complement system, was found to recognize and bind pneumonia-causing streptococci in the lung epithelium [28,29]. It is noteworthy that many patients with COVID-19 (e.g., those with severe disease) commonly display the same pathological phenotypes, namely lung inflammation, fibrosis, and pneumonia, linked to those ACE2-co-expressed genes. It is intriguing to suggest that perhaps this small population of ACE2-expressing cells may underlie the pathogenesis of severe ARDS, pneumonia, and respiratory failure in COVID-19 patients and, perhaps, particularly in those with LUAD. Our findings on SARS-CoV-2 receptor expression patterns in uninvolved normal lungs of patients with LUAD support the need for future studies comparing the expression of those genes between COVID-19 patients with and without LUAD.
A notable finding in our study was the co-expression of pathogens, including viral- scavengers and receptors such as CD36 and DMBT1 in ACE2-positive AT2 cells. We also found that among all lung subsets, AT2 cells distinctly and markedly displayed the highest expression (fraction and level) of DMBT1. Further, DMBT1 correlated with AT2 fractions in independent cohorts of bulk-sequenced lung tissues and positively correlated with ACE2 in the AT2 compartment. DMBT1, also known as gp340, was shown to inhibit influenza A by binding to hemagglutinin on the virus [50]. DMBT1 was also shown to interact with surfactant protein D in alveolar cells to agglutinate and inhibit influenza A virus [31]. It has been suggested that the antiviral drug oseltamivir cooperates with innate immune proteins such as DMBT1 in the inhibition of lung epithelial cell infection by influenza A virus [50]. Interestingly soluble DMBT1 in saliva was shown to exert host defense roles (neutralization or inhibition of oral transmission) against influenza A virus [50] as well as HIV-1 [51]. DMBT1 was shown to specifically inhibit HIV-1 infectivity by binding to the virus envelope protein gp120, the same protein that binds to the CD4 receptor on host T cells [52]. In a manner similar to the ACE2-mediated entry route facilitating SARS-CoV-2 infection in the lung [7], DMBT1 was shown to aid in HIV-1 transcytosis across genital tract tissue [53]. The study by Stoddard and colleagues also demonstrated that HIV-1 transport can be inhibited by antibodies or peptides that block the interaction of DMBT1 with the HIV-1 envelope protein gp120 [53]. Given our findings on DMBT1 co-expression with ACE2 in lung AT2 cells, as well as its reported binding to multiple viruses, we suggest that targeting DMBT1 using soluble peptides [54] or by antibody-based neutralization may represent a viable strategy to counteract SARS-CoV-2 infection and ameliorate COVID-19.

5. Conclusions

Single-cell transcriptomic analysis of our lung epithelial single-cell cohort demonstrated that among all lung cell subsets, cells of AT2 and malignant-enriched clusters displayed the highest relative expression of the SARS-CoV-2 receptor ACE2, albeit at low cell fractions. Our data also highlight alveolar and malignant subsets that express the serine proteases TMPRSS2 and TMPRSS4 implicated in lung pathobiology, including cancer, which also act as SARS-CoV-2 coreceptors. We further found that the viral scavenger DMBT1 is highly expressed in AT2 cells and correlates with ACE2 in this compartment. Our study underscores epithelial cell populations in the lungs of LUAD patients that express SARS-CoV-2 receptors as well as genes involved in inflammatory lung pathological conditions and host defense, and, thus, points to potential targets underlying the development of COVID-19 in lung cancer patients.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/13/6/1250/s1, Figure S1: Single-cell expression of ACE2, TMPRSS2, and TMPRSS4 in cycling lung cells, Figure S2: Single-cell clustering of lung epithelial cells and analysis of ACE2, TMPRSS2, and TMPRSS4 expression in normal versus tumor LUAD TCGA samples, Figure S3: Single-cell expression of ACE2, TMPRSS2, and TMPRSS4 in AT2 subclusters, Figure S4: Single-cell expression of ACE2 in AT2 and malignant-enriched subsets of losartan-treated versus untreated patients, Figure S5: Single-cell expression of DMBT1 and the coronavirus receptors BSG and DPP4 in lung epithelial cells, Figure S6: Correlations between estimated AT2 cell fractions and ACE2, DMBT1, TMPRSS2, and TMPRSS4 in TCGA normal lung samples of stage I or stages II and III LUADs, Table S1: Clinicopathological variables of the LUAD patients analyzed by scRNA-seq, Table S2: scRNA-seq quality metrics and number of EPCAM+ cells analyzed in tumor and multiple normal lung tissues, Table S3: Cell numbers of total and ACE2+ cells in epithelial subclusters.

Author Contributions

Conceptualization, L.W. and H.K.; Data curation, G.H., A.S., K.H., P.B., C.B., L.M.S., J.F., I.I.W., A.E.S., and H.K.; Formal analysis, G.H., A.S., C.B., S.J.M., E.J.O., P.S., J.F., L.W., and H.K.; Funding acquisition, J.D.M., I.I.W., C.S.S., A.E.S., L.W., and H.K.; Investigation, G.H., A.S., K.H., W.T., P.B., K.C., E.B., B.S.-E., L.M.S., B.G., L.G., J.Z., B.S., T.C., L.A.B., D.L.G., J.C., S.J.M., E.J.O., P.S., J.F., J.S., J.V.H., J.D.M., S.D., I.I.W., C.S.S., A.E.S., L.W., and H.K.; Methodology, G.H., A.S., K.H., W.T., P.B., K.C., E.B., J.C., J.F., L.W., and H.K.; Project administration, A.E.S., L.W., and H.K.; Resources, B.S.-E., B.S., T.C., L.A.B., D.L.G., J.V.H., J.D.M., S.D., I.I.W., C.S.S., and A.E.S.; Software, G.H. and L.W.; Supervision, C.S.S., L.W., and H.K.; Writing—original draft, G.H., A.S., L.W., and H.K.; Writing—review and editing, G.H., A.S., K.H., W.T., P.B., K.C., E.B., C.B., L.M.S., B.G., L.G., J.Z., B.S., T.C., L.A.B., D.L.G., J.C., S.J.M., E.J.O., P.S., J.F., J.S., J.V.H., J.D.M., S.D., I.I.W., C.S.S., A.E.S., L.W., and H.K. All authors have read and agreed to the published version of the manuscript.

Funding

Supported in part by research funding from Johnson and Johnson (to H.K.), National Cancer Institute (NCI) grants R01CA205608 (to H.K.), University of Texas SPORE in Lung Cancer P50CA070907 (to J.D.M.) and P50 core grant CA016672 (ATGC), and NIH grant 1S10OD024977.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of the University of Texas MD Anderson Cancer Center (PA14-0077, approved 7 April 2020).

Informed Consent Statement

Samples from all subjects involved in the study were obtained under waiver of consent from banked or residual tissues.

Data Availability Statement

All sequencing data generated in this study are deposited in the European Genome-phenome Archive (EGAS00001005021).

Conflicts of Interest

C.S.S. and A.E.S. are employees of Johnson and Johnson. H.K. receives research support from Johnson & Johnson. T.C. reports consulting fees from MedImmune/AstraZeneca and Bristol-Myers Squibb, advisory role fees from Bristol-Myers Squibb, MedImmune/AstraZeneca and EMD Serono, and clinical research funding to MD Anderson Cancer Center from Boehringer Ingelheim, MedImmune/AstraZeneca, EMD Serono, and Bristol-Myers Squibb.

References

  1. World Health Organization. Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19), 14–20 Februray 2020; WHO: Geneva, Switzerland, 2020. [Google Scholar]
  2. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef]
  3. Garassino, M.C.; Whisenant, J.G.; Huang, L.-C.; Trama, A.; Torri, V.; Agustoni, F.; Baena, J.; Banna, G.; Berardi, R.; Bettini, A.C.; et al. COVID-19 in patients with thoracic malignancies (TERAVOLT): First results of an international, registry-based, cohort study. Lancet Oncol. 2020, 21, 914–922. [Google Scholar] [CrossRef]
  4. Zhou, P.; Yang, X.-L.; Wang, X.-G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.-R.; Zhu, Y.; Li, B.; Huang, C.-L.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [Green Version]
  5. Kuba, K.; Imai, Y.; Rao, S.; Gao, H.; Guo, F.; Guan, B.; Huan, Y.; Yang, P.; Zhang, Y.; Deng, W.; et al. A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus–induced lung injury. Nat. Med. 2005, 11, 875–879. [Google Scholar] [CrossRef] [PubMed]
  6. Walls, A.C.; Park, Y.J.; Tortorici, M.A.; Wall, A.; McGuire, A.T.; Veesler, D. Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell 2020, 181, 281–292. [Google Scholar] [CrossRef]
  7. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 271–280. [Google Scholar] [CrossRef]
  8. Imai, Y.; Kuba, K.; Rao, S.; Huan, Y.; Guo, F.; Guan, B.; Yang, P.; Sarao, R.; Wada, T.; Leong-Poi, H.; et al. Angiotensin-converting enzyme 2 protects from severe acute lung failure. Nature 2005, 436, 112–116. [Google Scholar] [CrossRef] [PubMed]
  9. Wan, H.; Feng, Y.; Liu, J.; Zhang, R.; Ma, Q.; Han, B.; Xiang, Y.; Che, J.; Cao, H.; Fei, X.; et al. The angiotensin-converting enzyme 2 in tumor growth and tumor-associated angiogenesis in non-small cell lung cancer. Oncol. Rep. 2010, 23, 941–948. [Google Scholar] [CrossRef]
  10. Zang, R.; Castro, M.F.G.; McCune, B.T.; Zeng, Q.; Rothlauf, P.W.; Sonnek, N.M.; Liu, Z.; Brulois, K.F.; Wang, X.; Greenberg, H.B.; et al. TMPRSS2 and TMPRSS4 promote SARS-CoV-2 infection of human small intestinal enterocytes. Sci. Immunol. 2020, 5, 3582. [Google Scholar] [CrossRef]
  11. Han, G.; Sinjab, A.; Treekitkarnmongkol, W.; Brennan, P.; Hara, K.; Chang, K.; Bogatenkova, E.; Sanchez-Espiridion, B.; Behrens, C.; Gao, B.; et al. Single-cell analysis of human lung epithelia reveals concomitant expression of the SARS-CoV-2 receptor ACE2 with multiple virus receptors and scavengers in alveolar type II cells. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  12. Satija, R.; Farrell, J.A.; Gennert, D.; Schier, A.F.; Regev, A. Spatial reconstruction of single-cell gene expression data. Nat. Biotechnol. 2015, 33, 495–502. [Google Scholar] [CrossRef] [Green Version]
  13. Becht, E.; McInnes, L.; Healy, J.; Dutertre, C.-A.; Kwok, I.W.H.; Ng, L.G.; Ginhoux, F.; Newell, E.W. Dimensionality reduction for visualizing single-cell data using UMAP. Nat. Biotechnol. 2019, 37, 38–44. [Google Scholar] [CrossRef]
  14. Kiselev, V.Y.; Kirschner, K.; Schaub, M.T.; Andrews, T.; Yiu, A.; Chandra, T.; Natarajan, K.N.; Reik, W.; Barahona, M.; Green, A.R.; et al. SC3: Consensus clustering of single-cell RNA-seq data. Nat. Methods 2017, 14, 483–486. [Google Scholar] [CrossRef] [Green Version]
  15. Aizarani, N.; Saviano, A.; Sagar; Mailly, L.; Durand, S.; Herman, J.S.; Pessaux, P.; Baumert, T.F.; Grün, D. A human liver cell atlas reveals heterogeneity and epithelial progenitors. Nature 2019, 572, 199–204. [Google Scholar] [CrossRef] [PubMed]
  16. Ota, C.; Ng-Blichfeldt, J.-P.; Korfei, M.; Alsafadi, H.N.; Lehmann, M.; Skronska-Wasek, W.; De Santis, M.M.; Guenther, A.; Wagner, D.E.; Königshoff, M. Dynamic expression of HOPX in alveolar epithelial cells reflects injury and repair during the progression of pulmonary fibrosis. Sci. Rep. 2018, 8, 19283. [Google Scholar] [CrossRef] [Green Version]
  17. Jain, R.; Barkauskas, C.E.; Takeda, N.; Bowie, E.J.; Aghajanian, H.; Wang, Q.; Padmanabhan, A.; Manderfield, L.J.; Gupta, M.; Li, D.; et al. Plasticity of Hopx+ type I alveolar cells to regenerate type II cells in the lung. Nat. Commun. 2015, 6, 1–11. [Google Scholar] [CrossRef] [PubMed]
  18. Plasschaert, L.W.; Žilionis, R.; Choo-Wing, R.; Savova, V.; Knehr, J.; Roma, G.; Klein, A.M.; Jaffe, A.B. A single-cell atlas of the airway epithelium reveals the CFTR-rich pulmonary ionocyte. Nature 2018, 560, 377–381. [Google Scholar] [CrossRef]
  19. Montoro, D.T.; Haber, A.L.; Biton, M.; Vinarsky, V.; Lin, B.; Birket, S.E.; Yuan, F.; Chen, S.; Leung, H.M.; Villoria, J.; et al. A revised airway epithelial hierarchy includes CFTR-expressing ionocytes. Nature 2018, 560, 319–324. [Google Scholar] [CrossRef]
  20. Travaglini, K.J.; Nabhan, A.N.; Penland, L.; Sinha, R.; Gillich, A.; Sit, R.V.; Chang, S.; Conley, S.D.; Mori, Y.; Seita, J.; et al. A molecular cell atlas of the human lung from single-cell RNA sequencing. Nature 2020, 587, 619–625. [Google Scholar] [CrossRef] [PubMed]
  21. Laughney, A.M.; Hu, J.; Campbell, N.R.; Bakhoum, S.F.; Setty, M.; Lavallée, V.-P.; Xie, Y.; Masilionis, I.; Carr, A.J.; Kottapalli, S.; et al. Regenerative lineages and immune-mediated pruning in lung cancer metastasis. Nat. Med. 2020, 26, 259–269. [Google Scholar] [CrossRef]
  22. Cancer Genome Atlas Research Network. Comprehensive molecular profiling of lung adenocarcinoma. Nature 2014, 511, 543–550. [Google Scholar] [CrossRef]
  23. Klimas, J.; Olvedy, M.; Ochodnicka-Mackovicova, K.; Kruzliak, P.; Cacanyiova, S.; Kristek, F.; Krenek, P.; Ochodnicky, P. Perinatally administered losartan augments renal ACE2 expression but not cardiac or renal Mas receptor in spontaneously hypertensive rats. J. Cell. Mol. Med. 2015, 19, 1965–1974. [Google Scholar] [CrossRef]
  24. Chuang, P.-T.; McMahon, A.P. Vertebrate Hedgehog signalling modulated by induction of a Hedgehog-binding protein. Nature 1999, 397, 617–621. [Google Scholar] [CrossRef]
  25. Zhou, X.; Baron, R.M.; Hardin, M.; Cho, M.H.; Zielinski, J.; Hawrylkiewicz, I.; Sliwinski, P.; Hersh, C.P.; Mancini, J.D.; Lu, K.; et al. Identification of a chronic obstructive pulmonary disease genetic determinant that regulates HHIP. Hum. Mol. Genet. 2011, 21, 1325–1335. [Google Scholar] [CrossRef] [Green Version]
  26. Simpson-Haidaris, P.J.; Courtney, M.A.; Wright, T.W.; Goss, R.; Harmsen, A.; Gigliotti, F. Induction of fibrinogen expression in the lung epithelium during Pneumocystis carinii pneumonia. Infect. Immun. 1998, 66, 4431–4439. [Google Scholar] [CrossRef] [PubMed]
  27. Duan, H.O.; Simpson-Haidaris, P.J. Cell type-specific differential induction of the human gamma-fibrinogen promoter by interleukin-6. J. Biol. Chem. 2006, 281, 12451–12457. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Haleem, K.S.; Ali, Y.M.; Yesilkaya, H.; Kohler, T.; Hammerschmidt, S.; Andrew, P.W.; Schwaeble, W.J.; Lynch, N.J. The Pneumococcal Surface Proteins PspA and PspC Sequester Host C4-Binding Protein to Inactivate Complement C4b on the Bacterial Surface. Infect. Immun. 2018, 87. [Google Scholar] [CrossRef] [Green Version]
  29. Buffalo, C.Z.; Bahn-Suh, A.J.; Hirakis, S.P.; Biswas, T.; Amaro, R.E.; Nizet, V.; Ghosh, P. Conserved patterns hidden within group A Streptococcus M protein hypervariability recognize human C4b-binding protein. Nat. Microbiol. 2016, 1, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Silverstein, R.L.; Febbraio, M. CD36, a Scavenger Receptor Involved in Immunity, Metabolism, Angiogenesis, and Behavior. Sci. Signal. 2009, 2, re3. [Google Scholar] [CrossRef] [Green Version]
  31. Madsen, J.; Mollenhauer, J.; Holmskov, U. Review: Gp-340/DMBT1 in mucosal innate immunity. Innate Immun. 2010, 16, 160–167. [Google Scholar] [CrossRef]
  32. Wang, K.; Chen, W.; Zhou, Y.S.; Lian, J.Q.; Zhang, Z.; Du, P.; Gong, L.; Zhang, Y.; Cui, H.Y.; Geng, J.J.; et al. SARS-CoV-2 invades host cells via a novel route: CD147-spike protein. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  33. Raj, V.S.; Mou, H.; Smits, S.L.; Dekkers, D.H.W.; Müller, M.A.; Dijkman, R.; Muth, D.; Demmers, J.A.A.; Zaki, A.; Fouchier, R.A.M.; et al. Dipeptidyl peptidase 4 is a functional receptor for the emerging human coronavirus-EMC. Nature 2013, 495, 251–254. [Google Scholar] [CrossRef] [Green Version]
  34. Mossel, E.C.; Wang, J.; Jeffers, S.; Edeen, K.E.; Wang, S.; Cosgrove, G.P.; Funk, C.J.; Manzer, R.; Miura, T.A.; Pearson, L.D.; et al. SARS-CoV replicates in primary human alveolar type II cell cultures but not in type I-like cells. Virology 2008, 372, 127–135. [Google Scholar] [CrossRef] [Green Version]
  35. Hamming, I.; Timens, W.; Bulthuis, M.L.C.; Lely, A.T.; Navis, G.V.; van Goor, H. Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis. J. Pathol. 2004, 203, 631–637. [Google Scholar] [CrossRef]
  36. Uhal, B.D.; Li, X.; Xue, A.; Gao, X.; Abdul-Hafez, A. Regulation of alveolar epithelial cell survival by the ACE-2/angiotensin 1-7/Mas axis. Am. J. Physiol. Lung Cell Mol. Physiol. 2011, 301, 269–274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Qian, Z.; Travanty, E.A.; Oko, L.; Edeen, K.; Berglund, A.; Wang, J.; Ito, Y.; Holmes, K.V.; Mason, R.J. Innate Immune Response of Human Alveolar Type II Cells Infected with Severe Acute Respiratory Syndrome–Coronavirus. Am. J. Respir. Cell Mol. Biol. 2013, 48, 742–748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Zhao, Y.; Zhao, Z.; Wang, Y.; Zhou, Y.; Ma, Y.; Zuo, W. Single-Cell RNA Expression Profiling of ACE2, the Receptor of SARS-CoV-2. Am. J. Respir. Crit. Care Med. 2020, 202, 756–759. [Google Scholar] [CrossRef]
  39. Bertram, S.; Glowacka, I.; Blazejewska, P.; Soilleux, E.; Allen, P.; Danisch, S.; Steffen, I.; Choi, S.-Y.; Park, Y.; Schneider, H.; et al. TMPRSS2 and TMPRSS4 Facilitate Trypsin-Independent Spread of Influenza Virus in Caco-2 Cells. J. Virol. 2010, 84, 10016–10025. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Boöttcher, E.; Matrosovich, T.; Beyerle, M.; Klenk, H.-D.; Garten, W.; Matrosovich, M. Proteolytic Activation of Influenza Viruses by Serine Proteases TMPRSS2 and HAT from Human Airway Epithelium. J. Virol. 2006, 80, 9896–9898. [Google Scholar] [CrossRef] [Green Version]
  41. Matsuyama, S.; Nagata, N.; Shirato, K.; Kawase, M.; Takeda, M.; Taguchi, F. Efficient Activation of the Severe Acute Respiratory Syndrome Coronavirus Spike Protein by the Transmembrane Protease TMPRSS2. J. Virol. 2010, 84, 12658–12664. [Google Scholar] [CrossRef] [Green Version]
  42. Valero, A.M.; Cisneros, J.; Ramírez, R.; Gaxiola, M.; Becerril, C.; Pardo, A.; Selman, M. TMPRSS4: A novel serine protease involved in IPF development? FASEB J. 2016, 30, 50–58. [Google Scholar]
  43. Aberasturi, A.L.; Redrado, M.; Villalba, M.; Larzabal, L.; Pajares, M.J.; Garcia, J.; Evans, S.R.; Garcia-Ros, D.; Bodegas, M.E.; Lopez, L.; et al. TMPRSS4 induces cancer stem cell-like properties in lung cancer cells and correlates with ALDH expression in NSCLC patients. Cancer Lett. 2016, 370, 165–176. [Google Scholar] [CrossRef]
  44. Exposito, F.; Villalba, M.; Redrado, M.; De Aberasturi, A.L.; Cirauqui, C.; Redin, E.; Guruceaga, E.; De Andrea, C.; Vicent, S.; Ajona, D.; et al. Targeting of TMPRSS4 sensitizes lung cancer cells to chemotherapy by impairing the proliferation machinery. Cancer Lett. 2019, 453, 21–33. [Google Scholar] [CrossRef] [PubMed]
  45. Liang, S.; Willis, J.; Dou, J.; Mohanty, V.; Huang, Y.; Vilar, E.; Chen, K. Sensei: How many samples to tell evolution in single-cell studies? bioRxiv 2020. [Google Scholar] [CrossRef]
  46. Vaduganathan, M.; Vardeny, O.; Michel, T.; McMurray, J.J.; Pfeffer, M.A.; Solomon, S.D. Renin-Angiotensin-Aldosterone System Inhibitors in Patients with Covid-19. N. Engl. J. Med. 2020, 382, 1653–1659. [Google Scholar] [CrossRef]
  47. Sungnak, W.; Huang, N.; Bécavin, C.; Berg, M.; Queen, R.; Litvinukova, M.; Talavera-López, C.; Maatz, H.; Reichart, D.; Sampaziotis, F.; et al. SARS-CoV-2 entry factors are highly expressed in nasal epithelial cells together with innate immune genes. Nat. Med. 2020, 26, 681–687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Muus, C.; Luecken, M.D.; Eraslan, G.; Sikkema, L.; Waghray, A.; Heimberg, G.; Kobayashi, Y.; Vaishnav, E.D.; Subramanian, A.; Smillie, C.; et al. Single-cell meta-analysis of SARS-CoV-2 entry genes across tissues and demographics. Nat. Med. 2021, 1–14. [Google Scholar] [CrossRef]
  49. Barnes, P.J. Inflammatory mechanisms in patients with chronic obstructive pulmonary disease. J. Allergy Clin. Immunol. 2016, 138, 16–27. [Google Scholar] [CrossRef] [Green Version]
  50. White, M.R.; Crouch, E.; Van Eijk, M.; Hartshorn, M.; Pemberton, L.; Tornoe, I.; Holmskov, U.; Hartshorn, K.L. Cooperative anti-influenza activities of respiratory innate immune proteins and neuraminidase inhibitor. Am. J. Physiol. Cell. Mol. Physiol. 2005, 288, L831–L840. [Google Scholar] [CrossRef]
  51. Wu, Z.; Van Ryk, D.; Davis, C.; Abrams, W.R.; Chaiken, I.; Magnani, J.; Malamud, D. Salivary Agglutinin Inhibits HIV Type 1 Infectivity through Interaction with Viral Glycoprotein 120. AIDS Res. Hum. Retroviruses 2003, 19, 201–209. [Google Scholar] [CrossRef] [PubMed]
  52. Yoon, V.; Fridkis-Hareli, M.; Munisamy, S.; Lee, J.; Anastasiades, D.; Stevceva, L. The GP120 molecule of HIV-1 and its interaction with T cells. Curr. Med. Chem. 2010, 17, 741–749. [Google Scholar] [CrossRef] [PubMed]
  53. Stoddard, E.; Ni, H.; Cannon, G.; Zhou, C.; Kallenbach, N.; Malamud, D.; Weissman, D. gp340 Promotes Transcytosis of Human Immunodeficiency Virus Type 1 in Genital Tract-Derived Cell Lines and Primary Endocervical Tissue. J. Virol. 2009, 83, 8596–8603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Tuttolomondo, M.; Casella, C.; Hansen, P.L.; Polo, E.; Herda, L.M.; Dawson, K.A.; Ditzel, H.J.; Mollenhauer, J. Human DMBT1-Derived Cell-Penetrating Peptides for Intracellular siRNA Delivery. Mol. Ther. Nucleic Acids 2017, 8, 264–276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 13 01250 g001 550
Figure 1. Overview of 186,916 tumor- and normal-derived lung cells analyzed by scRNA-seq. (A) Bubble plot show-ing the expression of markers indicative of major lineages in the lung ecosystem. Both the fraction of cells expressing (indicated by the size of the circle) as well as their scaled expression levels (indicated by the color of the circle) are shown. (B) Uniform manifold approximation and projection (UMAP) embedding of cells from tumor and multiple nor-mal samples per patient (19 samples from 5 LUAD patients). Cells are colored by expression level of ACE2, TMPRSS2, or TMPRSS4.
Figure 1. Overview of 186,916 tumor- and normal-derived lung cells analyzed by scRNA-seq. (A) Bubble plot show-ing the expression of markers indicative of major lineages in the lung ecosystem. Both the fraction of cells expressing (indicated by the size of the circle) as well as their scaled expression levels (indicated by the color of the circle) are shown. (B) Uniform manifold approximation and projection (UMAP) embedding of cells from tumor and multiple nor-mal samples per patient (19 samples from 5 LUAD patients). Cells are colored by expression level of ACE2, TMPRSS2, or TMPRSS4.
Cancers 13 01250 g001
Cancers 13 01250 g002 550
Figure 2. Single-cell expression analysis of ACE2, TMPRSS2, and TMPRSS4 in 70,030 lung epithelial cells. (A–C) Uniform Manifold Approximation and Projection (UMAP) plots showing epithelial subclusters and expression of ACE2 (A), TMPRSS2 (B), and TMPRSS4 (C) in epithelial subclusters. Cells are colored by the expression level of each gene. AT1: alveolar type 1 cells; AT2: alveolar type 2 cells. (D–F) Bar plots showing the fraction of cells (percentage of each subcluster) expressing ACE2 (D), TMPRSS2 (E), and TMPRSS4 (F) among airway lineage clusters with the highest fractions of cells positive for each gene (2% cutoff applied). The absolute number of cells positive for each gene and within each analyzed subcluster are indicated on top of the corresponding bars. Color indicates average expression in cells positive for the gene of interest.
Figure 2. Single-cell expression analysis of ACE2, TMPRSS2, and TMPRSS4 in 70,030 lung epithelial cells. (A–C) Uniform Manifold Approximation and Projection (UMAP) plots showing epithelial subclusters and expression of ACE2 (A), TMPRSS2 (B), and TMPRSS4 (C) in epithelial subclusters. Cells are colored by the expression level of each gene. AT1: alveolar type 1 cells; AT2: alveolar type 2 cells. (D–F) Bar plots showing the fraction of cells (percentage of each subcluster) expressing ACE2 (D), TMPRSS2 (E), and TMPRSS4 (F) among airway lineage clusters with the highest fractions of cells positive for each gene (2% cutoff applied). The absolute number of cells positive for each gene and within each analyzed subcluster are indicated on top of the corresponding bars. Color indicates average expression in cells positive for the gene of interest.
Cancers 13 01250 g002
Cancers 13 01250 g003 550
Figure 3. Differentially expressed genes between ACE2-positive and -negative AT2 cells. (A) Volcano plot showing significantly differentially expressed genes (DEGs) between ACE2-positive (n = 607) and -negative (n = 26,628) AT2 cells. A cutoff of absolute gene expression (fold-change: >1.2) and a FDR (q-value < 0.05) were applied to identify the DEGs. Blue indicates downregulation, and orange indicates upregulation. (B) Violin plots showing the significant upregulation of FGG, DMBT1, and HHIP genes in ACE2-positive compared with -negative AT2 cells. FC: Fold change.
Figure 3. Differentially expressed genes between ACE2-positive and -negative AT2 cells. (A) Volcano plot showing significantly differentially expressed genes (DEGs) between ACE2-positive (n = 607) and -negative (n = 26,628) AT2 cells. A cutoff of absolute gene expression (fold-change: >1.2) and a FDR (q-value < 0.05) were applied to identify the DEGs. Blue indicates downregulation, and orange indicates upregulation. (B) Violin plots showing the significant upregulation of FGG, DMBT1, and HHIP genes in ACE2-positive compared with -negative AT2 cells. FC: Fold change.
Cancers 13 01250 g003
Cancers 13 01250 g004 550
Figure 4. Expression of the viral scavenger DMBT1 correlates with that of ACE2 in lung AT2 cells. (A) Correlation between ACE2 and DMBT1, TMPRSS2, and TMPRSS4 expression in pseudobulk data from this study. (B) Scatter plots showing significant correlations between estimated AT2 cell fractions and ACE2, DMBT1, TMPRSS2, and TMPRSS4 in TCGA normal lung samples. (C) Significant correlation between ACE2 and DMBT1 in TCGA normal lung samples with high AT2 cell fractions (meta-score > 15.48). Correlations were statistically analyzed using Pearson’s correlation coefficient.
Figure 4. Expression of the viral scavenger DMBT1 correlates with that of ACE2 in lung AT2 cells. (A) Correlation between ACE2 and DMBT1, TMPRSS2, and TMPRSS4 expression in pseudobulk data from this study. (B) Scatter plots showing significant correlations between estimated AT2 cell fractions and ACE2, DMBT1, TMPRSS2, and TMPRSS4 in TCGA normal lung samples. (C) Significant correlation between ACE2 and DMBT1 in TCGA normal lung samples with high AT2 cell fractions (meta-score > 15.48). Correlations were statistically analyzed using Pearson’s correlation coefficient.
Cancers 13 01250 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Han, G.; Sinjab, A.; Hara, K.; Treekitkarnmongkol, W.; Brennan, P.; Chang, K.; Bogatenkova, E.; Sanchez-Espiridion, B.; Behrens, C.; Solis, L.M.; et al. Single-Cell Expression Landscape of SARS-CoV-2 Receptor ACE2 and Host Proteases in Normal and Malignant Lung Tissues from Pulmonary Adenocarcinoma Patients. Cancers 2021, 13, 1250. https://doi.org/10.3390/cancers13061250

AMA Style

Han G, Sinjab A, Hara K, Treekitkarnmongkol W, Brennan P, Chang K, Bogatenkova E, Sanchez-Espiridion B, Behrens C, Solis LM, et al. Single-Cell Expression Landscape of SARS-CoV-2 Receptor ACE2 and Host Proteases in Normal and Malignant Lung Tissues from Pulmonary Adenocarcinoma Patients. Cancers. 2021; 13(6):1250. https://doi.org/10.3390/cancers13061250

Chicago/Turabian Style

Han, Guangchun, Ansam Sinjab, Kieko Hara, Warapen Treekitkarnmongkol, Patrick Brennan, Kyle Chang, Elena Bogatenkova, Beatriz Sanchez-Espiridion, Carmen Behrens, Luisa M. Solis, and et al. 2021. "Single-Cell Expression Landscape of SARS-CoV-2 Receptor ACE2 and Host Proteases in Normal and Malignant Lung Tissues from Pulmonary Adenocarcinoma Patients" Cancers 13, no. 6: 1250. https://doi.org/10.3390/cancers13061250

APA Style

Han, G., Sinjab, A., Hara, K., Treekitkarnmongkol, W., Brennan, P., Chang, K., Bogatenkova, E., Sanchez-Espiridion, B., Behrens, C., Solis, L. M., Gao, B., Girard, L., Zhang, J., Sepesi, B., Cascone, T., Byers, L. A., Gibbons, D. L., Chen, J., Moghaddam, S. J., ... Kadara, H. (2021). Single-Cell Expression Landscape of SARS-CoV-2 Receptor ACE2 and Host Proteases in Normal and Malignant Lung Tissues from Pulmonary Adenocarcinoma Patients. Cancers, 13(6), 1250. https://doi.org/10.3390/cancers13061250

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop