Open access
Research Article
1 September 2020

Low Baseline Pulmonary Levels of Cytotoxic Lymphocytes as a Predisposing Risk Factor for Severe COVID-19

ABSTRACT

Coronavirus disease 2019 (COVID-19) is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and currently has detrimental human health, community, and economic impacts around the world. It is unclear why some SARS-CoV-2-positive individuals remain asymptomatic, while others develop severe symptoms. Baseline pulmonary levels of antiviral leukocytes, already residing in the lung prior to infection, may orchestrate an effective early immune response and prevent severe symptoms. Here, “in silico flow cytometry” was used to deconvolute the levels of all seven types of antiviral leukocytes in 1,927 human lung tissues. Baseline levels of CD8+ T cells, resting NK cells, and activated NK cells, as well as cytokines that recruit these cells, are significantly lower in lung tissues with high expression of the SARS-CoV-2 entry receptor angiotensin-converting enzyme 2 (ACE2). This is observed in univariate analyses, in multivariate analyses, and in two independent data sets. Importantly, ACE2 mRNA and protein levels very strongly correlate in human cells and tissues. The above findings also largely apply to the SARS-CoV-2 entry protease TMPRSS2. Both SARS-CoV-2-infected lung cells and COVID-19 lung tissues show upregulation of CD8+ T cell- and NK cell-recruiting cytokines. Moreover, tissue-resident CD8+ T cells and inflammatory NK cells are significantly more abundant in bronchoalveolar lavage fluids from mildly affected COVID-19 patients compared to severe cases. This suggests that these lymphocytes are important for preventing severe symptoms. Elevated ACE2 expression increases sensitivity to coronavirus infection. Thus, the results suggest that some individuals may be exceedingly susceptible to develop severe COVID-19 due to concomitant high preexisting ACE2 and TMPRSS expression and low baseline cytotoxic lymphocyte levels in the lung.
IMPORTANCE COVID-19 is caused by the highly contagious coronavirus SARS-CoV-2 and currently has detrimental human health, community, and economic impacts around the world. It is unclear why some SARS-CoV-2-positive individuals develop severe COVID-19 symptoms, which can be fatal, while others only develop mild symptoms. In the absence of an effective and widely available vaccine, it is of paramount importance that we identify risk factors for development of severe symptoms to be able to improve treatment approaches. The ACE2 gene encodes the receptor on human cells that the virus uses to infect these cells. This study finds that if the lungs of healthy individuals have high levels of ACE2, they typically have low levels of the immune cells that eliminate viruses. Therefore, some individuals may develop severe COVID-19 due to simultaneous high levels of the virus receptor and low levels of immune cells that eradicate the virus in their lungs.

INTRODUCTION

Coronaviruses are viruses belonging to the family Coronaviridae (1). They are large, single-stranded RNA viruses that often originate from bats and commonly infect mammals. While the majority of coronavirus infections cause mild symptoms, some can cause severe symptoms, such as pneumonia, respiratory failure and sepsis, which may lead to death (2, 3).
Coronavirus zoonosis constitutes a serious health risk for humans. Indeed, in recent history, transmissions of three types of coronaviruses to humans have led to various numbers of deaths. The outbreak of the severe acute respiratory syndrome (SARS) epidemic, which is caused by the SARS coronavirus (SARS-CoV), originated in Guangdong, China in 2002 and led to nearly 800 deaths (4). The Middle East respiratory syndrome coronavirus (MERS-CoV) outbreak, which emerged in Saudi Arabia in 2012, similarly caused about 800 deaths but with more than 8,000 cases, nearly four times as many cases were reported (4). Finally, coronavirus disease 19 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently causing a pandemic. On 1 May 2020, the World Health Organization reported over 3 million confirmed cases and over 220,000 patients who have succumbed to COVID-19 around the world (5). However, the actual number of deaths is probably considerably higher (6). In addition, this figure is still soaring; on 1 May 2020, the rate exceeded 6,400 deaths per day (5).
To infect target cells, coronaviruses use their spike (S) glycoprotein to bind to receptor molecules on the host cell membrane. Angiotensin-converting enzyme 2 (ACE2) has been identified as the main SARS-CoV-2 entry receptor on human cells (7, 8), while the serine protease TMPRSS2, or potentially cathepsin B and L, are used for S-protein priming to facilitate host cell entry (7). SARS-CoV-2 S protein has a 10- to 20-fold-higher affinity to human ACE2 than SARS-CoV S protein (9). Moreover, ACE2 expression proportionally increases the susceptibility to S protein-mediated coronavirus infection (1012). Hence, increased expression of ACE2 is thought to increase susceptibility to COVID-19 (1315).
Epithelial cells of the respiratory tract, including the lung, are primary SARS-CoV-2 target cells (1618). These cells can sense viral infection via pattern recognition receptors (PRRs). PRRs, including Toll-like receptors and NOD-like receptors, recognize pathogen-associated molecular patterns (PAMPs) (19). Upon PRR activation, a range of proinflammatory cytokines and chemokines are produced and released in order to activate the host’s immune system. Interferons (IFNs), in particular type I and type III IFN, are among the principal cytokines to recruit immune cells (19, 20).
Six types of leukocytes have been implicated in detecting and responding to viral infections in the lung, a major site of SARS-CoV-2 infection, which also presents with severe COVID-19 symptoms. The cytotoxic activities of CD8+ T cells and NK cells can facilitate early control of viral infections by clearing infected cells and avoiding additional viral dissemination (21, 22). Dendritic cells specialize in sensing infections, including by viruses, and inducing an immune response (23). CD4+ T cells contribute to viral clearance by promoting production of cytokines and interactions between CD8+ T cells and dendritic cells (24). M1 macrophages interact with pulmonary epithelial cells to fight viral infections in the lung (25). Finally, neutrophils may contribute to clearance of viral infections through phagocytosis of virions and viral particles. However, their precise role is uncertain (26).
SARS-CoV-2 is considerably more efficient in infection, replication, and production of infectious virus particles in human lung tissue than SARS-CoV (17). Strikingly, despite this, SARS-CoV-2 initially does not significantly induce type I, II, or III IFNs in infected human lung cells and tissue (17, 27). When this does occur, it may in fact promote further SARS-CoV-2 infection, as IFNs directly upregulate expression of the SARS-CoV-2 receptor ACE2 (28). These observations suggest that baseline levels of leukocytes, which already reside in the lung prior to infection, may be important in mounting a rapid immune response against SARS-CoV-2 infection and prevent severe COVID-19 symptoms. As stated above, ACE2 expression level may be a predictor of increased susceptibility to COVID-19 (1015). Thus, here I investigated the relationship between ACE2 and TMPRSS2 expression and the levels of seven leukocyte types implicated in antiviral immune response in human lung tissue.

RESULTS

I used bulk transcriptome sequencing (RNAseq) gene expression data from the 578 human lung tissues present in the Genotype-Tissue Expression (GTEx) database (29, 30), because this is the largest publicly available data set with clinical information. Using an established “in silico flow cytometry” pipeline (31), I estimated the levels of CD8+ T cells, resting and activated NK cells, M1 macrophages, dendritic cells, CD4+ T cells, and neutrophils in these tissues (see Fig. S1a to c and Table S1 in the supplemental material). These levels were compared to ACE2 expression levels in these lung tissues. This revealed that ACE2 expression is negatively correlated with the levels of CD8+ T cells, resting and activated NK cells, and M1 macrophages (P < 8 × 10−6, Pearson correlations) (Fig. 1a to c). However, there are no statistically significant correlations between ACE2 expression and the levels of CD4+ T cells, dendritic cells, and neutrophils (P > 0.05) (Fig. S2a to d). Thus, the levels of a majority of leukocytes involved in antiviral immune responses are significantly lower in lung tissues with high ACE2 expression levels.
FIG 1
FIG 1 Baseline levels of cytotoxic lymphocytes inversely correlate with ACE2 expression in the lung. (a to c) The correlations between the baseline levels of CD8+ T cells, resting NK cells, and activated NK cells (x axes) and those of the SARS-CoV-2 host cell receptor ACE2 (y axes) in human lung tissue are shown. Data are from the GTEx data set (n = 578). (d to f) The same correlations are shown for lung tissues in the LUG data set (n = 1,349). Regression lines and 95% confidence intervals are shown. Pearson correlation R and P values are shown, q values are Benjamini-Hochberg-adjusted P values using a false discovery rate of 0.05, and mp values are multivariate P values (Table 1).
It is possible that some of the above observations are linked to phenotypic characteristics, such as sex, age, body mass index (BMI), race, or smoking status. To test the robustness of our findings, I applied multivariable regression analysis that includes these five covariates (Table S2), as well as the levels of the above seven leukocyte types or states. This showed that only 4 of the 12 variables significantly contribute to predicting ACE2 expression levels, specifically the levels of CD8+ T cells, resting NK cells, activated NK cells, and M1 macrophages (Fig. 1a to c and Table 1). Notably, none of the five added phenotypic covariates showed statistically significant contributions. Consistently, I found limited statistically significant correlations between these variables and ACE2 expression in univariate analyses, irrespective of whether they were analyzed as continuous data or binned into discrete ordinal categories (Fig. S3a to j). Thus, the levels of four types of leukocytes that respond to viral infection are low in lung tissue with high ACE2 expression levels independently of phenotypic covariates.
TABLE 1
TABLE 1 Multivariate analysis of ACE2 expression in human lung tissue
Variableβ valueSEat valueP valueP value symbolb
CD8+ T cells−2.7190.584−4.6564.02e−06****
NK cells, resting−3.1280.583−5.3661.18e−07****
NK cells, activated−2.2810.704−3.2410.0013**
M1 macrophages−4.9631.339−3.7070.0002***
CD4+ T cells0.4530.5130.8820.3781 
Dendritic cells−12.7936.688−1.9130.0563 
Neutrophils−0.3980.469−0.8480.3971 
Sex = male (vs female)0.0650.0561.160.2466 
Age0.0040.0021.7990.0726 
BMI−0.0040.006−0.5560.5787 
Race = Caucasian (vs non-Caucasian)−0.0960.073−1.3130.1899 
Smoking status = yes (vs no)−0.0250.058−0.4330.6650 
a
SE, standard error.
b
P value symbols: **, P < 0.01; ***, P < 0.001; ****, P < 0.0001.
Next, whether the above observations could be validated in an independent cohort of individuals was tested. For this, I used, to my knowledge, the largest publicly available lung tissue data set. The Laval University, University of British-Columbia, Groningen University (LUG) data set, including microarray gene expression data of 1,349 human lung tissues, was used. Following determination of ACE2 expression levels and estimation of the levels of CD8+ T cells, resting NK cells, activated NK cells, and M1 macrophages (Fig. S1a to c and Table S1), I found that three of the four also negatively correlated with ACE2 expression in this independent data set (P < 2 × 10−8) (Fig. 1d to f). With a correlation coefficient of R = 0.096, only M1 macrophages did not correlate with ACE2 expression in this data set (Fig. S4). Thus, our observations indicate that the baseline levels of three types of cytotoxic lymphocytes, specifically CD8+ T cells, resting NK cells, and activated NK cells, are robustly and consistently low in lung tissue with high expression of the SARS-CoV-2 receptor ACE2.
To more rigorously assess our observations, a range of additional analyses was employed. Although highly statistically significant (all P values < 4.1 × 10−6, Fig. 1), the absolute Pearson R values between baseline levels of ACE2 and the three lymphocyte types seemed to be low, as they ranged between 0.2 and 0.3 (Fig. 1a to c). To test how strong these are in relative terms, the Pearson R and P values of 1,000 randomly sampled other genes were calculated. This revealed that the ACE2 R and P values were significantly lower than expected by chance (all one-sample t test P < 2.2 × 10−16; Fig. 2a and b). In addition, these R and P values ranked in the top 0.4 to 11 percentiles of strongest and most significant correlations for each of the three leukocyte types (Fig. 2a and b). Thus, the seemingly low correlations between ACE2 mRNA and cytotoxic lymphocyte levels in the lung are not only highly statistically significant but also strong in relative terms.
FIG 2
FIG 2 Levels of ACE2 mRNA, ACE2 protein, CD8+ T cells, NK cells, and cytokines in lung cells, lung tissues, and COVID-19 patient samples. (a and b) Pearson R and −log10 P values of correlations between 1,000 randomly sampled genes and the levels of indicated lymphocytes in lung tissues were determined and plotted. The NK cells were resting and activated NK cells. ACE2 Pearson R and P values are shown as red diamonds. Blue lines indicate means with 95% confidence intervals. Percentiles for ACE2 with respect to the 1,000 random R and P values are shown. Data are from the GTEx data set (n = 578). The P values are from one-sample t tests. (c) Correlations between ACE2 mRNA and protein levels in 52 cell lines. R and P values are from Pearson correlations. (d) Pearson correlation R values between mRNA and protein levels of 12,016 genes are compared to the ACE2 R coefficient (red diamond). The line and box represent the median and interquartile ranges. The ACE2 R percentile is also shown. The P value is from one-sample t test. (e) Bar graph showing the correlation between ACE2 mRNA and protein levels in human tissues. Means plus standard errors of the means (error bars) are shown. Samples are from the Human Protein Atlas. P values were from Kruskal-Wallis test. (f) Meta-analysis scatterplot showing the correlation between ACE2 mRNA and protein levels in 17 human tissues. Data are from nine different studies, as detailed in Table S3 in the supplemental material. R and P values are from Pearson correlation. (g and h) Heatmaps showing Spearman correlations between the levels of ACE2 or indicated cytotoxic lymphocytes and eight cytokines that recruit these cells in human lung tissues from the GTEx data set (n = 578). The colors of the tiles represent Spearman R, per the color bar on the right. Spearman significance levels are shown by asterisks. See also Fig. S5. (i and j) Fold increase in expression levels of indicated cytokines in Calu-3 lung cells 24 h after SARS-CoV-2 infection compared to uninfected Calu-3 cells (i), and in postmortem COVID-19 lung tissues (n = 2) to those in healthy, uninfected lung tissues (n = 2) (j). (k) Comparison of indicated fractions of lymphocyte levels in bronchoalveolar lavage fluids from mild/moderate and severe COVID-19 patients, as determined in two separate studies (46, 47). The second study also determined a tissue resident (TR) score for CD8+ T cells. Numbers in the mild/moderate column on the left show fold increase compared to the respective severe cases on the right. Asterisks in the severe column on the right represent statistical significance levels, as determined by Mann-Whitney U tests (top panel), or t tests (middle and bottom panels) comparing mild/moderate to severe cases. P value symbols: *, P < 0.05; **, P < 0.01; ***, P < 0.001; ****, P < 0.0001; †, P < 2.2 × 10−16.
Further, how well ACE2 mRNA and protein levels correlate was tested. Using the mRNA and protein levels in 52 cell lines, ACE2 mRNA levels were found to strongly correlate with ACE2 protein levels in human cells (Pearson R = 0.8155, P = 1.8 × 10−13; Fig. 2c). In fact, ACE2 ranks in the top 1.8 percentile of over 12,000 genes with the strongest mRNA-protein level correlations (P < 2.2 × 10−16; Fig. 2d). Using two ACE2-specific antibodies, immunochemistry on 40 human tissues also shows a strong ACE2 mRNA-protein correlation (P = 0.0011, Kruskal-Wallis test; Fig. 2e), and this is additionally validated by a meta-analysis conducted using nine published studies (Pearson R = 0.7130, P = 0.0013; Fig. 2f and Table S3). Therefore, I conclude that ACE2 mRNA and protein levels very strongly correlate, both in human cells and in human tissues.
Above, I found that baseline ACE2 levels in the lung negatively correlate with CD8+ T cells and resting and activated NK cells in multivariate analyses and in an independent data set (Fig. 1a to f). Several cytokines, including C-C motif chemokine ligand 2 (CCL2) to CCL5, C-X-C motif chemokine ligand 9 (CXCL9), CXCL10, CXCL16, and interleukin 18 (IL-18), are known to chemotactically attract CD8+ T cells and NK cells (3237). Consistently, I find that the baseline levels of these chemokines in human lung tissue typically significantly correlate with the baseline levels of CD8+ T cells and resting and activated NK cells (Fig. 2g and Fig. S5). Additionally, as expected given the results above, significant negative correlations were found between the levels of ACE2 and the levels of six of these eight cytokines in the lung (Fig. 2h). These findings lend further support to my previous observations, suggesting that high levels of said cytokines in the lung establish a favorable milieu for cytotoxic lymphocytes, which correlates with low ACE2 levels.
Next, the direct consequences of SARS-CoV-2 infection were assessed. In vitro SARS-CoV-2 infection of human lung cells invariably leads to upregulation of all eight above-mentioned CD8+ T cell- and NK cell-attracting cytokines, with six of these increases showing statistical significance (Fig. 2i). Similarly, compared to control lung tissues, all eight cytokines are upregulated in lung tissues from COVID-19 patients, with five showing statistical significance (Fig. 2j). Moreover, the levels of CD8+ T cells and NK cells are higher in bronchoalveolar lavage fluids from mildly affected COVID-19 patients than from severe cases, with CD8+ T cells and a subset of NK cells, inflammatory NK cells, showing a statistically significant higher level (Fig. 2k). These findings are corroborated in a different cohort of patients, additionally showing a highly significant increase in a tissue-resident signature score for CD8+ T cells (Fig. 2k). Thus, together, these observations suggest that SARS-CoV-2 infection of lung cells stimulates CD8+ T cell- and NK cell-attracting cytokines and that these cytotoxic lymphocytes are important for preventing severe symptoms of COVID-19.
Finally, whether the levels of the SARS-CoV-2 host cell protease TMPRSS2 shows similar correlations with the levels of CD8+ T cells and NK cells in the lung was tested. In univariate analyses, baseline TMPRSS2 levels in the lung show significant negative correlations with these lymphocyte levels, although in multivariate analyses, these correlations are statistically significant only for CD8+ T cells and activated NK cells (Fig. 3a to c). The corresponding R and P values are typically also significantly lower than expected by chance (Fig. 3d and e). Furthermore, TMPRSS2 mRNA and protein levels strongly correlate (R = 0.8048, P < 2.2 × 10−16, Fig. 3f), and TMPRSS2 is in the top 2.5 percentile of genes that show the strongest mRNA-protein correlation (P < 2.2 × 10−16, Fig. 3g). Additionally, TMPRSS2 expression tends to correlate negatively with CD8+ T cell- and NK cell-attracting cytokines (Fig. 3h). Therefore, albeit typically to a lesser extent, baseline TMPRSS2 expression levels in the lung negatively correlate with the levels of CD8+ T cells and NK cells in a manner similar to ACE2.
FIG 3
FIG 3 Levels of TMPRSS2 mRNA, TMPRSS2 protein, and cytokines in lung cells and tissues. (a to c) Pearson correlations between baseline levels of indicated lymphocytes and TMPRSS2 in human lung tissue, as in Fig. 1a to c. Data are from the GTEx data set (n = 578). (d and e) Pearson R and −log10 P values of correlations between 1,000 randomly sampled genes and the levels of indicated lymphocytes in lung tissues, as in Fig. 2a and b. TMPRSS2 Pearson R and P values are shown as green diamonds. (f) Correlations between TMPRSS2 mRNA and protein levels in 124 cell lines. R and P values are from Pearson correlations. (g) Pearson correlation R values between mRNA and protein levels of 12,016 genes are compared to the TMPRSS2 R coefficient (green diamond). The line and box represent the median and interquartile ranges, respectively. The TMPRSS2 R percentile is also shown. The P value is from one-sample t test. (h) Heatmap showing Spearman correlations between the levels of TMPRSS2 and cytokines in human lung tissues from the GTEx data set (n = 578), as in Fig. 2h. (i) Individuals with high baseline levels of ACE2 and TMPRSS2 show low baseline tissue-resident levels of cytotoxic lymphocytes in the lung. I propose that this may jointly predispose these individuals to development of severe COVID-19. P value symbols: *, P < 0.05; **, P < 0.01; ****, P < 0.0001; †, P < 2.2 × 10−16.
Taken together, these observations suggest that a subgroup of individuals may be exceedingly susceptible to developing severe COVID-19 due to concomitant high preexisting ACE2 and TMPRSS2 expression and low baseline levels of CD8+ T cells and NK cells in the lung (Fig. 3i).

DISCUSSION

I investigated the baseline expression levels of the SARS-CoV-2 host cell entry receptor ACE2 and the host cell entry protease TMPRSS2 and the baseline levels of all seven types of antiviral leukocytes in 1,927 human lung tissue samples. Although SARS-CoV-2 cellular tropism is broad (1618), I focused on lung tissue. In addition to epithelial cells elsewhere in the respiratory tract, alveolar epithelial cells are thought to be a primary SARS-CoV-2 entry point (16, 28). Consistently, the SARS-CoV-2 receptor ACE2 is expressed in these cells at the mRNA and protein levels (28, 3840). Moreover, in severely affected COVID-19 patients, the lungs are among the few organs that present with the most life-threatening symptoms. “Cytokine storm”-induced acute respiratory distress syndrome (ARDS), widespread alveolar damage, pneumonia, and progressive respiratory failure have been observed (41, 42). These indications frequently require admission to intensive care units (ICUs), and mechanical ventilation and may ultimately be fatal.
Early after infection, rapid activation of the innate immune system is of paramount importance for the clearance of virus infections. Infected cells typically do so through release of proinflammatory cytokines and chemokines, in particular type I and III interferons (19, 20). Notably, however, several studies have highlighted multiple complexities related specifically to SARS-CoV2 and innate immune system activation at early stages. First, unlike SARS-CoV, SARS-CoV-2-infected lung tissue initially fails to induce a range of immune cell-recruiting molecules, including several interferons (17, 27), suggesting that leukocytes are ineffectively recruited to the infected lung shortly after infection. Second, the host cell entry receptor ACE2 has been identified as an interferon target gene (28). Thus, even when interferons are upregulated in order to recruit immune cells, concomitant upregulation of ACE2 expression may in fact exacerbate SARS-CoV-2 infection (28).
These findings suggest that the levels of immune cells that already reside in the lung prior to infection may be more critical for dampening SARS-CoV-2 infection at early stages than they are for fighting infections of other viruses. Cytotoxic lymphocytes, including CD8+ T cells and NK cells, are key early responders to virus infections, and these are the cells whose baseline levels were identified here as significantly reduced in lung tissue with elevated ACE2 and TMPRSS2 expression. That these immune cells are important in preventing severe COVID-19 is supported by the fact that their levels are significantly higher in bronchoalveolar lavage fluids from patients with mild cases than from patients with severe cases. Therefore, our results suggest that individuals with increased baseline susceptibility to SARS-CoV-2 infection in the lungs may also be less well equipped from the outset to mount a rapid antiviral cellular immune response (Fig. 3i).
Several observations indicate that these cytotoxic lymphocytes are critically important for effective control of SARS-CoV-2 infection. Recent studies showed that CD8+ T cells in peripheral blood are considerably reduced and functionally exhausted in COVID-19 patients, in particular in elderly patients and in severely affected patients that require ICU admission (4345). Reduced CD8+ T cell counts also predict poor COVID-19 patient survival (43). Additionally, CD8+ T cell- and NK cell-attracting cytokines are upregulated in SARS-CoV-2-infected human lung cells and in lung tissues from COVID-19 patients, and the levels of CD8+ T cells and NK cells are higher in bronchoalveolar lavage fluids from mildly affected COVID-19 patients than in patients with severe disease (46, 47).
The five phenotypic parameters, sex, age, BMI, race, and smoking history, did not statistically significantly contribute to variation in ACE2 expression in human lung tissue, either in univariate or in multivariate analyses. This is consistent with some studies but inconsistent with others (42, 4850). These paradoxical observations may be partially explained by differing gender, age, and race distributions within each study cohort.
Further research will be required to elucidate the precise mechanisms of SARS-CoV-2-induced activation of the innate immune system early after infection. However, the link identified here between high baseline ACE2 and TMPRSS2 expression and reduced cytotoxic lymphocyte levels in human lung tissue prior SARS-CoV-2 infection is striking. It suggests that increased susceptibility to SARS-CoV-2 infection in the lungs may be accompanied by a poorer ability to mount a rapid innate immune response at early stages. This may predict long-term outcome of individuals infected with SARS-CoV-2, given that the levels of CD8+ T cells and NK cells are significantly higher in bronchoalveolar lavage fluids from mild cases compared to severe patients (46, 47). Finally, it may contribute to the substantial variation in COVID-19 clinical presentation, ranging from asymptomatic to severe respiratory and other symptoms.

MATERIALS AND METHODS

Discovery data set and processing.

Gene expression data and corresponding phenotype data from human lung tissues (n = 578) were obtained from the Genotype-Tissue Expression (GTEx) Portal (https://gtexportal.org), managed by the National Institutes of Health (NIH). Gene expression data were publicly available. Access to phenotype data required authorization. The GTEx protocol was previously described (29, 30). Briefly, total RNA was extracted from tissue. Following mRNA isolation, cDNA synthesis, and library preparation, samples were subjected to HiSeq2000 or HiSeq2500 Illumina TrueSeq RNA sequencing. Gene expression levels were obtained using RNA-SeQC v1.1.9 (51) and expressed in transcripts per million (TPM). Reported gene-level expression levels were log2 transformed, unless indicated otherwise.

Validation data set and processing.

For validation purposes, the Laval University, University of British-Columbia, Groningen University (LUG) lung tissue data set (n = 1,349) was used. This data set was accessed via Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo), accession number GSE23546, and was previously described (52). Briefly, total RNA from human lung tissue samples was isolated, quantified, quality checked, and used to generate cDNA, which was amplified and hybridized to Affymetrix gene expression arrays. Arrays were scanned, and probe-level gene expression values were normalized using robust multichip average (RMA). These normalized values were obtained from GEO. To collapse probe-level expression data to single expression levels per gene, for each gene the probe with the highest median absolute deviation (MAD) was used. The MAD for each probe p was calculated using equation 1
MAD(p)=M(|piM(p)|)
(1)
where M is the median, pi denotes probe p's expression level in sample i, and M (p) represents the median signal of probe p.

In silico cytometry.

The levels of seven types of leukocytes involved in antiviral cellular immune response, specifically CD8+ T cells, resting NK cells, activated NK cells, M1 macrophages, CD4+ T cells, dendritic cells, and neutrophils, were estimated in the discovery and validation lung tissue samples using a previously described approach (31). Specifically, the following workflow was used. First, only non-log-transformed expression values were used. Thus, where required, expression values for all samples in the discovery and validation data sets were reverse log2 transformed using equation 2
c=2cl1
(2)
where c denotes the calculated non-log2-transformed expression counts and cl denotes the previously reported log2-transformed expression counts. Next, to compensate for potential technical differences between signatures and bulk sample gene expression values due to interplatform variation, bulk mode batch correction was applied. To ensure robustness, deconvolution was statistically analyzed using 100 permutations. Pearson correlation coefficients R, root mean squared errors (RMSA), and P values are reported on a per-sample level in Fig. S1a to c and Table S1 in the supplemental material.

Univariate statistical analyses.

Log2-transformed expression levels of ACE2 and TMPRSS2 in lung tissue samples were compared to the estimated levels of seven leukocyte types or states. Pearson correlation analyses were performed to determine Pearson correlation coefficients R and P values. P values were adjusted at a false discovery rate of 0.05 to yield q values, as previously described (53). Straight lines represent the minimized sum of squares of deviations of the data points with 95% confidence intervals shown. Continuous phenotypic covariates were analyzed in the same way and, additionally, as discrete ordinal categories after binning. Discrete and binned phenotype data were statistically evaluated using Mann-Whitney U tests. All analyses were performed in the R computing environment (R Project for Statistical Computing, Vienna, Austria).

Multivariate regression analyses.

Multivariate analyses were performed using standard ordinary least-squares regression, summarized in equation 3.
Y=β0+k=1n(βkXk)+ε
(3)
where β0 denotes the intercept, while βk represents the slope of each variable Xk in a model with n variables and ε denotes the random error component. These analyses were performed using R.

mRNA and protein levels in human cells.

Available ACE2 and TMPRSS2 mRNA and corresponding protein levels in 52 and 124 human cell lines, respectively, were obtained from references 54 and 55. The correlations between mRNA and protein levels were analyzed by linear regression analysis using Pearson correlations. Pearson coefficients R for mRNA-protein level correlations were also determined for 12,015 other genes. To test whether the ACE2 and TMPRSS2 coefficients and P values were statistically significantly lower than for other genes, one-sample t tests were used.

mRNA and protein levels in human tissues.

Two types of analyses were performed to compare mRNA and protein levels in human tissues. First, mRNA expression levels from 40 human tissues were obtained from the Human Protein Atlas (https://www.proteinatlas.org). These levels represented the consensus normalized mRNA expression levels from three sources, specifically, Human Protein Atlas in-house RNAseq data, RNAseq data from the Genotype-Tissue Expression (GTEx) project and CAGE data from FANTOM5 project. Corresponding ordinal human tissue protein expression levels (“not detected,” “low,” and “'high”) were also obtained from the Human Protein Atlas. These levels were based on immunohistochemical staining of the tissues using DAB (3,3′-diaminobenzidine)-labeled antibodies (HPA000288 and CAB026174), followed by knowledge-based annotation, as described on the website. A Kruskal-Wallis test was performed to assess whether mRNA and protein levels significantly correlated.
Second, a meta-analysis was performed. For this, mRNA and protein expression levels were obtained from nine different sources. Tissue mRNA levels were obtained from six sources, as determined by Northern blotting (56), quantitative reverse transcription-PCR (RT-PCR) (57), microarray hybridization (58), RNA sequencing (59), and cap analysis of gene expression (60). Tissue protein levels were obtained from three sources, as determined by mass spectrometry (61) and immunohistochemistry (39, 62). For each source and tissue, mRNA and protein expression levels were scored as “not detected,” “low,” “intermediate,” or “high” and these received scores of 0 to 3, respectively. For each tissue, the final mRNA and protein scores were calculated by averaging the scores from the respective six and three sources (Table S3). Strength and significance level of the correlation between the final scores were determined by linear regression analysis using Pearson correlation.

Cytokine levels in SARS-CoV-2-infected lung cells.

SARS-CoV-2-induced fold increases in the expression levels of eight cytotoxic lymphocyte-attracting cytokines, CCL2-5, CXCL9, CXCL10, CXCL16, and IL-18, were determined from reference 27. These increases represent the fold increase in expression in Calu-3 lung cells, 24 h after infection with SARS-CoV-2 at a multiplicity of infection of 2, compared to uninfected Calu-3 cells.

Cytokine levels in control and COVID-19 lung tissues.

Baseline levels of the above eight cytokines in lung tissues were obtained from the GTEx project, as described above, and compared to the baseline levels of ACE2, TMPRSS2, CD8+ T cells, resting and activated NK cells in these tissues, estimated as described above. Their levels were compared using Spearman’s rank correlations (R and P values). For comparison of cytokine levels in postmortem COVID-19 lung tissues (n = 2) to those in healthy, uninfected lung tissues (n = 2), fold increases were determined following RNAseq analyses and previously reported (27).

Lymphocyte levels in COVID-19 bronchoalveolar lavage fluid samples.

The levels of CD8+ T cells, NK cells, and inflammatory NK cells in bronchoalveolar lavages from mild (n = 2) and severe (n = 22) COVID-19 patients were reported elsewhere and determined using single-cell RNAseq (46). Statistical significance levels were assessed using Mann-Whitney U tests. The levels of T cells and NK cells, as well as a CD8+ T cell tissue-resident signature score, in bronchoalveolar lavage fluid samples from moderate (n = 3) and severe/critical (n = 5) COVID-19 patients were reported in another study (47).

Data availability.

The data sets used for the analyses described in this study were obtained from dbGaP at https://www.ncbi.nlm.nih.gov/gap through dbGaP accession number phs000424.v8.p2.

ACKNOWLEDGMENTS

I thank Marianna Datseris for critically reading the manuscript.
This work was supported by the School of Biomedical Sciences, Queensland University of Technology. This study used publicly available and controlled data from the Genotype-Tissue Expression (GTEx) Project. GTEx was supported by the Common Fund of the Office of the Director of the National Institutes of Health (commonfund.nih.gov/GTEx). Additional funds were provided by the NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. Donors were enrolled at Biospecimen Source Sites funded by NCI\Leidos Biomedical Research, Inc. subcontracts to the National Disease Research Interchange (10XS170), Roswell Park Cancer Institute (10XS171), and Science Care, Inc. (X10S172). The Laboratory, Data Analysis, and Coordinating Center (LDACC) was funded through a contract (HHSN268201000029C) to the Broad Institute, Inc., Biorepository operations were funded through a Leidos Biomedical Research, Inc. subcontract to Van Andel Research Institute (10ST1035). Additional data repository and project management were provided by Leidos Biomedical Research, Inc. (HHSN261200800001E). The Brain Bank was supported by supplements to University of Miami grant DA006227. Statistical Methods development grants were made to the University of Geneva (MH090941 and MH101814), the University of Chicago (MH090951, MH090937, MH101825, and MH101820), the University of North Carolina—Chapel Hill (MH090936), North Carolina State University (MH101819), Harvard University (MH090948), Stanford University (MH101782), Washington University (MH101810), and to the University of Pennsylvania (MH101822).

Supplemental Material

File (msystems.00741-20-s0001.docx)
File (msystems.00741-20-sf001.pdf)
File (msystems.00741-20-sf002.pdf)
File (msystems.00741-20-sf003.pdf)
File (msystems.00741-20-sf004.pdf)
File (msystems.00741-20-sf005.pdf)
File (msystems.00741-20-st001.docx)
File (msystems.00741-20-st002.docx)
File (msystems.00741-20-st003.docx)
ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.

REFERENCES

1.
Coronaviridae Study Group of the International Committee on Taxonomy of Viruses. 2020. The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2. Nat Microbiol 5:536–544.
2.
Liu Y, Yan LM, Wan L, Xiang TX, Le A, Liu JM, Peiris M, Poon LLM, Zhang W. 2020. Viral dynamics in mild and severe cases of COVID-19. Lancet Infect Dis 20:656–657.
3.
Gattinoni L, Chiumello D, Caironi P, Busana M, Romitti F, Brazzi L, Camporota L. 2020. COVID-19 pneumonia: different respiratory treatments for different phenotypes? Intensive Care Med 46:1099–1102.
4.
Xie M, Chen Q. 2020. Insight into 2019 novel coronavirus − an updated interim review and lessons from SARS-CoV and MERS-CoV. Int J Infect Dis 94:119–124.
5.
World Health Organization. 2020. Coronavirus disease 2019 (COVID-19): situation report 102. World Health Organization, Geneva, Switzerland.
6.
Li R, Pei S, Chen B, Song Y, Zhang T, Yang W, Shaman J. 2020. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science 368:489–493.
7.
Hoffmann M, Kleine-Weber H, Schroeder S, Kruger N, Herrler T, Erichsen S, Schiergens TS, Herrler G, Wu NH, Nitsche A, Muller MA, Drosten C, Pohlmann S. 2020. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 181:271–280.e8.
8.
Zhou P, Yang XL, Wang XG, Hu B, Zhang L, Zhang W, Si HR, Zhu Y, Li B, Huang CL, Chen HD, Chen J, Luo Y, Guo H, Jiang RD, Liu MQ, Chen Y, Shen XR, Wang X, Zheng XS, Zhao K, Chen QJ, Deng F, Liu LL, Yan B, Zhan FX, Wang YY, Xiao GF, Shi ZL. 2020. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579:270–273.
9.
Wrapp D, Wang N, Corbett KS, Goldsmith JA, Hsieh CL, Abiona O, Graham BS, McLellan JS. 2020. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 367:1260–1263.
10.
Hofmann H, Geier M, Marzi A, Krumbiegel M, Peipp M, Fey GH, Gramberg T, Pohlmann S. 2004. Susceptibility to SARS coronavirus S protein-driven infection correlates with expression of angiotensin converting enzyme 2 and infection can be blocked by soluble receptor. Biochem Biophys Res Commun 319:1216–1221.
11.
Mossel EC, Huang C, Narayanan K, Makino S, Tesh RB, Peters CJ. 2005. Exogenous ACE2 expression allows refractory cell lines to support severe acute respiratory syndrome coronavirus replication. J Virol 79:3846–3850.
12.
Li W, Moore MJ, Vasilieva N, Sui J, Wong SK, Berne MA, Somasundaran M, Sullivan JL, Luzuriaga K, Greenough TC, Choe H, Farzan M. 2003. Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus. Nature 426:450–454.
13.
Fang L, Karakiulakis G, Roth M. 2020. Are patients with hypertension and diabetes mellitus at increased risk for COVID-19 infection? Lancet Respir Med 8:e21.
14.
Guo J, Huang Z, Lin L, Lv J. 2020. Coronavirus disease 2019 (COVID-19) and cardiovascular disease: a viewpoint on the potential influence of angiotensin-converting enzyme inhibitors/angiotensin receptor blockers on onset and severity of severe acute respiratory syndrome coronavirus 2 infection. J Am Heart Assoc 9:e016219.
15.
Zhuang MW, Cheng Y, Zhang J, Jiang XM, Wang L, Deng J, Wang PH. 2020. Increasing host cellular receptor-angiotensin-converting enzyme 2 expression by coronavirus may facilitate 2019-nCoV (or SARS-CoV-2) infection. J Med Virol.
16.
Sungnak W, Huang N, Becavin C, Berg M, Queen R, Litvinukova M, Talavera-Lopez C, Maatz H, Reichart D, Sampaziotis F, Worlock KB, Yoshida M, Barnes JL, HCA Lung Biological Network. 2020. SARS-CoV-2 entry factors are highly expressed in nasal epithelial cells together with innate immune genes. Nat Med 26:681–687.
17.
Chu H, Chan JF, Wang Y, Yuen TT, Chai Y, Hou Y, Shuai H, Yang D, Hu B, Huang X, Zhang X, Cai JP, Zhou J, Yuan S, Kok KH, To KK, Chan IH, Zhang AJ, Sit KY, Au WK, Yuen KY. 2020. Comparative replication and immune activation profiles of SARS-CoV-2 and SARS-CoV in human lungs: an ex vivo study with implications for the pathogenesis of COVID-19. Clin Infect Dis.
18.
Chu H, Chan JFW, Yuen TTT, Shuai H, Yuan S, Wang Y, Hu B, Yip CCY, Tsang JOL, Huang X, Chai Y, Yang D, Hou Y, Chik KKH, Zhang X, Fung AYF, Tsoi HY, Cai JP, Chan WM, Ip JD, Chu AWH, Zhou J, Lung DC, Kok KH, To KKW, Tsang OTY, Chan KH, Yuen KY. 2020. Comparative tropism, replication kinetics, and cell damage profiling of SARS-CoV-2 and SARS-CoV with implications for clinical manifestations, transmissibility, and laboratory studies of COVID-19: an observational study. Lancet Microbe 1:E14–E23.
19.
Newton AH, Cardani A, Braciale TJ. 2016. The host immune response in respiratory virus infection: balancing virus clearance and immunopathology. Semin Immunopathol 38:471–482.
20.
Crouse J, Kalinke U, Oxenius A. 2015. Regulation of antiviral T cell responses by type I interferons. Nat Rev Immunol 15:231–242.
21.
Schmidt ME, Varga SM. 2018. The CD8 T cell response to respiratory virus infections. Front Immunol 9:678.
22.
Hammer Q, Ruckert T, Romagnani C. 2018. Natural killer cell specificity for viral infections. Nat Immunol 19:800–808.
23.
Cook PC, MacDonald AS. 2016. Dendritic cells in lung immunopathology. Semin Immunopathol 38:449–460.
24.
Sant AJ, McMichael A. 2012. Revealing the role of CD4(+) T cells in viral immunity. J Exp Med 209:1391–1395.
25.
Harper RW. 2017. Partners in crime: epithelial priming of macrophages during viral infection. Am J Respir Cell Mol Biol 57:145–146.
26.
Galani IE, Andreakos E. 2015. Neutrophils in viral infections: current concepts and caveats. J Leukoc Biol 98:557–564.
27.
Blanco-Melo D, Nilsson-Payant BE, Liu WC, Uhl S, Hoagland D, Moller R, Jordan TX, Oishi K, Panis M, Sachs D, Wang TT, Schwartz RE, Lim JK, Albrecht RA, tenOever BR. 2020. Imbalanced host response to SARS-CoV-2 drives development of COVID-19. Cell 181:1036–1045.e9.
28.
Ziegler CGK, Allon SJ, Nyquist SK, Mbano IM, Miao VN, Tzouanas CN, Cao Y, Yousif AS, Bals J, Hauser BM, Feldman J, Muus C, Wadsworth MH, Kazer SW, Hughes TK, Doran B, Gatter GJ, Vukovic M, Taliaferro F, Mead BE, Guo Z, Wang JP, Gras D, Plaisant M, Ansari M, Angelidis I, Adler H, Sucre JMS, Taylor CJ, Lin B, Waghray A, Mitsialis V, Dwyer DF, Buchheit KM, Boyce JA, Barrett NA, Laidlaw TM, Carroll SL, Colonna L, Tkachev V, Peterson CW, Yu A, Zheng HB, Gideon HP, Winchell CG, Lin PL, Bingle CD, Snapper SB, Kropski JA, Theis FJ, et al. 2020. SARS-CoV-2 receptor ACE2 is an interferon-stimulated gene in human airway epithelial cells and is detected in specific cell subsets across tissues. Cell 181:1016–1035.e19.
29.
Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, Hasz R, Walters G, Garcia F, Young N, Foster B, Moser M, Karasik E, Gillard B, Ramsey K, Sullivan S, Bridge J, Magazine H, Syron J, Fleming J, Siminoff L, Traino H, Mosavel M, Barker L, Jewell S, Rohrer D, Maxim D, Filkins D, Harbach P, Cortadillo E, Berghuis B, Turner L, Hudson E, Feenstra K, Sobin L, Robb J, Branton P, Korzeniewski G, Shive C, Tabor D, Qi L, Groch K, Nampally S, Buia S, Zimmerman A, Smith A, Burges R, Robinson K, Valentino K, Bradbury D, Cosentino M, et al. 2013. The Genotype-Tissue Expression (GTEx) project. Nat Genet 45:580–585.
30.
GTEx Consortium, Laboratory, Data Analysis, & Coordinating Center (LDACC)—Analysis Working Group, Statistical Methods groups—Analysis Working Group, Enhancing GTEx groups (EGTEx), NIH Common Fund, NIH/NCI, NIH/NHGRI, NIH/NIMH, NIH/NIDA, Biospecimen Collection Source Site—NDRI, Biospecimen Collection Source Site—RPCI, Biospecimen Core Resource—VARI, Brain Bank Repository—University of Miami Brain Endowment Bank, Leidos Biomedical—Project Management, ELSI Study, Genome Browser Data Integration & Visualization—EBI, Genome Browser Data Integration & Visualization—UCSC Genomics Institute, University of California Santa Cruz, Lead Analysts, Laboratory, Data Analysis, & Coordinating Center (LDACC), NIH program management, Biospecimen collection, Pathology, eQTL manuscript working group, Battle A, Brown CD, Engelhardt BE, Montgomery SB. 2017. Genetic effects on gene expression across human tissues. Nature 550:204–213.
31.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. 2015. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods 12:453–457.
32.
Brown CE, Vishwanath RP, Aguilar B, Starr R, Najbauer J, Aboody KS, Jensen MC. 2007. Tumor-derived chemokine MCP-1/CCL2 is sufficient for mediating tumor tropism of adoptively transferred T cells. J Immunol 179:3332–3341.
33.
Serody JS, Burkett SE, Panoskaltsis-Mortari A, Ng-Cashin J, McMahon E, Matsushima GK, Lira SA, Cook DN, Blazar BR. 2000. T-lymphocyte production of macrophage inflammatory protein-1alpha is critical to the recruitment of CD8(+) T cells to the liver, lung, and spleen during graft-versus-host disease. Blood 96:2973–2980.
34.
Zumwalt TJ, Arnold M, Goel A, Boland CR. 2015. Active secretion of CXCL10 and CCL5 from colorectal cancer microenvironments associates with GranzymeB+ CD8+ T-cell infiltration. Oncotarget 6:2981–2991.
35.
Spranger S, Dai D, Horton B, Gajewski TF. 2017. Tumor-residing Batf3 dendritic cells are required for effector T cell trafficking and adoptive T cell therapy. Cancer Cell 31:711–723.e4.
36.
Wendel M, Galani IE, Suri-Payer E, Cerwenka A. 2008. Natural killer cell accumulation in tumors is dependent on IFN-gamma and CXCR3 ligands. Cancer Res 68:8437–8445.
37.
Fan Z, Yu P, Wang Y, Wang Y, Fu ML, Liu W, Sun Y, Fu YX. 2006. NK-cell activation by LIGHT triggers tumor-specific CD8+ T-cell immunity to reject established tumors. Blood 107:1342–1351.
38.
Bertram S, Heurich A, Lavender H, Gierer S, Danisch S, Perin P, Lucas JM, Nelson PS, Pohlmann S, Soilleux EJ. 2012. Influenza and SARS-coronavirus activating proteases TMPRSS2 and HAT are expressed at multiple sites in human respiratory and gastrointestinal tracts. PLoS One 7:e35876.
39.
Hamming I, Timens W, Bulthuis ML, Lely AT, Navis G, van Goor H. 2004. Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis. J Pathol 203:631–637.
40.
Zou X, Chen K, Zou J, Han P, Hao J, Han Z. 2020. Single-cell RNA-seq data analysis on the receptor ACE2 expression reveals the potential risk of different human organs vulnerable to 2019-nCoV infection. Front Med 14:185–192.
41.
Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, Zhang L, Fan G, Xu J, Gu X, Cheng Z, Yu T, Xia J, Wei Y, Wu W, Xie X, Yin W, Li H, Liu M, Xiao Y, Gao H, Guo L, Xie J, Wang G, Jiang R, Gao Z, Jin Q, Wang J, Cao B. 2020. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 395:497–506.
42.
Chan JF, Yuan S, Kok KH, To KK, Chu H, Yang J, Xing F, Liu J, Yip CC, Poon RW, Tsoi HW, Lo SK, Chan KH, Poon VK, Chan WM, Ip JD, Cai JP, Cheng VC, Chen H, Hui CK, Yuen KY. 2020. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet 395:514–523.
43.
Diao B, Wang C, Tan Y, Chen X, Liu Y, Ning L, Chen L, Li M, Liu Y, Wang G, Yuan Z, Feng Z, Zhang Y, Wu Y, Chen Y. 2020. Reduction and functional exhaustion of T cells in patients with coronavirus disease 2019 (COVID-19). Front Immunol 11:827.
44.
Zheng M, Gao Y, Wang G, Song G, Liu S, Sun D, Xu Y, Tian Z. 2020. Functional exhaustion of antiviral lymphocytes in COVID-19 patients. Cell Mol Immunol 17:533–535.
45.
Liu J, Li S, Liu J, Liang B, Wang X, Wang H, Li W, Tong Q, Yi J, Zhao L, Xiong L, Guo C, Tian J, Luo J, Yao J, Pang R, Shen H, Peng C, Liu T, Zhang Q, Wu J, Xu L, Lu S, Wang B, Weng Z, Han C, Zhu H, Zhou R, Zhou H, Chen X, Ye P, Zhu B, Wang L, Zhou W, He S, He Y, Jie S, Wei P, Zhang J, Lu Y, Wang W, Zhang L, Li L, Zhou F, Wang J, Dittmer U, Lu M, Hu Y, Yang D, Zheng X, et al. 2020. Longitudinal characteristics of lymphocyte responses and cytokine profiles in the peripheral blood of SARS-CoV-2 infected patients. EBioMedicine 55:102763.
46.
Wauters E, Van Mol P, Garg AD, Jansen S, Van Herck Y, Vanderbeke L, Bassez A, Boeckx B, Malengier-Devlies B, Timmerman A, Van Brussel T, Van Buyten T, Schepers R, Heylen E, Dauwe D, Dooms C, Gunst J, Hermans G, Meersseman P, Testelmans D, Yserbyt J, Matthys P, Tejpar S, CONTAGIOUS Collaborators, Neyts J, Wauters J, Qian J, Lambrechts D. 10 July 2020. Discriminating mild from critical COVID-19 by innate and adaptive immune single-cell profiling of bronchoalveolar lavages. bioRxiv preprint 196519.
47.
Liao M, Liu Y, Yuan J, Wen Y, Xu G, Zhao J, Cheng L, Li J, Wang X, Wang F, Liu L, Amit I, Zhang S, Zhang Z. 2020. Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19. Nat Med 26:842–844.
48.
Lloyd-Sherlock P, Ebrahim S, Geffen L, McKee M. 2020. Bearing the brunt of covid-19: older people in low and middle income countries. BMJ 368:m1052.
49.
Yang X, Yu Y, Xu J, Shu H, Xia J, Liu H, Wu Y, Zhang L, Yu Z, Fang M, Yu T, Wang Y, Pan S, Zou X, Yuan S, Shang Y. 2020. Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study. Lancet Respir Med 8:475–481.
50.
Smith JC, Sausville EL, Girish V, Yuan ML, Vasudevan A, John KM, Sheltzer JM. 2020. Cigarette smoke exposure and inflammatory signaling increase the expression of the SARS-CoV-2 receptor ACE2 in the respiratory tract. Dev Cell 53:514–529.e3.
51.
DeLuca DS, Levin JZ, Sivachenko A, Fennell T, Nazaire MD, Williams C, Reich M, Winckler W, Getz G. 2012. RNA-SeQC: RNA-seq metrics for quality control and process optimization. Bioinformatics 28:1530–1532.
52.
Obeidat M, Hao K, Bosse Y, Nickle DC, Nie Y, Postma DS, Laviolette M, Sandford AJ, Daley DD, Hogg JC, Elliott WM, Fishbane N, Timens W, Hysi PG, Kaprio J, Wilson JF, Hui J, Rawal R, Schulz H, Stubbe B, Hayward C, Polasek O, Jarvelin MR, Zhao JH, Jarvis D, Kahonen M, Franceschini N, North KE, Loth DW, Brusselle GG, Smith AV, Gudnason V, Bartz TM, Wilk JB, O’Connor GT, Cassano PA, Tang W, Wain LV, Soler Artigas M, Gharib SA, Strachan DP, Sin DD, Tobin MD, London SJ, Hall IP, Pare PD. 2015. Molecular mechanisms underlying variations in lung function: a systems genetics analysis. Lancet Respir Med 3:782–795.
53.
Benjamini Y, Hochberg Y. 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Stat Methodol 57:289–300.
54.
Ghandi M, Huang FW, Jané-Valbuena J, Kryukov GV, Lo CC, McDonald ER, Barretina J, Gelfand ET, Bielski CM, Li H, Hu K, Andreev-Drakhlin AY, Kim J, Hess JM, Haas BJ, Aguet F, Weir BA, Rothberg MV, Paolella BR, Lawrence MS, Akbani R, Lu Y, Tiv HL, Gokhale PC, de Weck A, Mansour AA, Oh C, Shih J, Hadi K, Rosen Y, Bistline J, Venkatesan K, Reddy A, Sonkin D, Liu M, Lehar J, Korn JM, Porter DA, Jones MD, Golji J, Caponigro G, Taylor JE, Dunning CM, Creech AL, Warren AC, McFarland JM, Zamanighomi M, Kauffmann A, Stransky N, Imielinski M, et al. 2019. Next-generation characterization of the Cancer Cell Line Encyclopedia. Nature 569:503–508.
55.
Nusinow DP, Szpyt J, Ghandi M, Rose CM, McDonald ER, III, Kalocsay M, Jané-Valbuena J, Gelfand E, Schweppe DK, Jedrychowski M, Golji J, Porter DA, Rejtar T, Wang YK, Kryukov GV, Stegmeier F, Erickson BK, Garraway LA, Sellers WR, Gygi SP. 2020. Quantitative proteomics of the Cancer Cell Line Encyclopedia. Cell 180:387–402.e16.
56.
Tipnis SR, Hooper NM, Hyde R, Karran E, Christie G, Turner AJ. 2000. A human homolog of angiotensin-converting enzyme. Cloning and functional expression as a captopril-insensitive carboxypeptidase. J Biol Chem 275:33238–33243.
57.
Harmer D, Gilbert M, Borman R, Clark KL. 2002. Quantitative mRNA expression profiling of ACE 2, a novel homologue of angiotensin converting enzyme. FEBS Lett 532:107–110.
58.
Wu C, Jin X, Tsueng G, Afrasiabi C, Su AI. 2016. BioGPS: building your own mash-up of gene annotations and expression profiles. Nucleic Acids Res 44:D313–D316.
59.
Uhlen M, Zhang C, Lee S, Sjostedt E, Fagerberg L, Bidkhori G, Benfeitas R, Arif M, Liu Z, Edfors F, Sanli K, von Feilitzen K, Oksvold P, Lundberg E, Hober S, Nilsson P, Mattsson J, Schwenk JM, Brunnstrom H, Glimelius B, Sjoblom T, Edqvist PH, Djureinovic D, Micke P, Lindskog C, Mardinoglu A, Ponten F. 2017. A pathology atlas of the human cancer transcriptome. Science 357:eaan2507.
60.
FANTOM Consortium and the RIKEN PMI and CLST (DGT), Forrest ARR, Kawaji H, Rehli M, Baillie JK, de Hoon MJL, Haberle V, Lassmann T, Kulakovskiy IV, Lizio M, Itoh M, Andersson R, Mungall CJ, Meehan TF, Schmeier S, Bertin N, Jørgensen M, Dimont E, Arner E, Schmidl C, Schaefer U, Medvedeva YA, Plessy C, Vitezic M, Severin J, Semple CA, Ishizu Y, Young RS, Francescatto M, Alam I, Albanese D, Altschuler GM, Arakawa T, Archer JAC, Arner P, Babina M, Rennie S, Balwierz PJ, Beckhouse AG, Pradhan-Bhatt S, Blake JA, Blumenthal A, Bodega B, Bonetti A, Briggs J, Brombacher F, Burroughs AM, Califano A, Cannistraci CV, Carbajo D, Chen Y, et al. 2014. A promoter-level mammalian expression atlas. Nature 507:462–470.
61.
Kim M-S, Pinto SM, Getnet D, Nirujogi RS, Manda SS, Chaerkady R, Madugundu AK, Kelkar DS, Isserlin R, Jain S, Thomas JK, Muthusamy B, Leal-Rojas P, Kumar P, Sahasrabuddhe NA, Balakrishnan L, Advani J, George B, Renuse S, Selvan LDN, Patil AH, Nanjappa V, Radhakrishnan A, Prasad S, Subbannayya T, Raju R, Kumar M, Sreenivasamurthy SK, Marimuthu A, Sathe GJ, Chavan S, Datta KK, Subbannayya Y, Sahu A, Yelamanchi SD, Jayaram S, Rajagopalan P, Sharma J, Murthy KR, Syed N, Goel R, Khan AA, Ahmad S, Dey G, Mudgal K, Chatterjee A, Huang T-C, Zhong J, Wu X, Shaw PG, et al. 2014. A draft map of the human proteome. Nature 509:575–581.
62.
Thul PJ, Akesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H, Alm T, Asplund A, Bjork L, Breckels LM, Backstrom A, Danielsson F, Fagerberg L, Fall J, Gatto L, Gnann C, Hober S, Hjelmare M, Johansson F, Lee S, Lindskog C, Mulder J, Mulvey CM, Nilsson P, Oksvold P, Rockberg J, Schutten R, Schwenk JM, Sivertsson A, Sjostedt E, Skogs M, Stadler C, Sullivan DP, Tegel H, Winsnes C, Zhang C, Zwahlen M, Mardinoglu A, Ponten F, von Feilitzen K, Lilley KS, Uhlen M, Lundberg E. 2017. A subcellular map of the human proteome. Science 356:eaal3321.

Information & Contributors

Information

Published In

cover image mSystems
mSystems
Volume 5Number 527 October 2020
eLocator: 10.1128/msystems.00741-20
Editor: Paola Flórez de Sessions, Oxford Nanopore Technologies

History

Received: 31 July 2020
Accepted: 15 August 2020
Published online: 1 September 2020

Keywords

  1. COVID-19
  2. SARS-CoV-2
  3. ACE2
  4. TMPRSS2
  5. T cells
  6. NK cells

Contributors

Author

Institute of Health and Biomedical Innovation, Queensland University of Technology (QUT), Faculty of Health, School of Biomedical Sciences, Brisbane, QLD, Australia
Centre for Data Science, Queensland University of Technology (QUT), Brisbane, QLD, Australia
University of Queensland Diamantina Institute, The University of Queensland, Translational Research Institute, Brisbane, QLD, Australia

Editor

Paola Flórez de Sessions
Editor
Oxford Nanopore Technologies

Notes

Address correspondence to [email protected].

Metrics & Citations

Metrics

Note:

  • For recently published articles, the TOTAL download count will appear as zero until a new month starts.
  • There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
  • Citation counts come from the Crossref Cited by service.

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. For an editable text file, please select Medlars format which will download as a .txt file. Simply select your manager software from the list below and click Download.

View Options

Figures and Media

Figures

Media

Tables

Share

Share

Share the article link

Share with email

Email a colleague

Share on social media

American Society for Microbiology ("ASM") is committed to maintaining your confidence and trust with respect to the information we collect from you on websites owned and operated by ASM ("ASM Web Sites") and other sources. This Privacy Policy sets forth the information we collect about you, how we use this information and the choices you have about how we use such information.
FIND OUT MORE about the privacy policy