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 (10–12). Hence, increased expression of ACE2 is thought to increase susceptibility to COVID-19 (13–15).
Epithelial cells of the respiratory tract, including the lung, are primary SARS-CoV-2 target cells (16–18). 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 (10–15). 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
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
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
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 (32–37). 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
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 (16–18), 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, 38–40). 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 (43–45). 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, 48–50). 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
(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
(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.
(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).
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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
Duijf PHG 2020. Low Baseline Pulmonary Levels of Cytotoxic Lymphocytes as a Predisposing Risk Factor for Severe COVID-19. mSystems5:10.1128/msystems.00741-20.
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.
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.