Immune response to SARS-CoV-2 in the nasal mucosa in children and adults
========================================================================

* Clarissa M Koch
* Andrew D Prigge
* Kishore R Anekalla
* Avani Shukla
* Hanh Chi Do-Umehara
* Leah Setar
* Jairo Chavez
* Hiam Abdala-Valencia
* Yuliya Politanska
* Nikolay S Markov
* Grant R Hahn
* Taylor Heald-Sargent
* L Nelson Sanchez-Pinto
* William J Muller
* Alexander V Misharin
* Karen M Ridge
* Bria M Coates

## Abstract

**Rationale** Despite similar viral load and infectivity rates between children and adults infected with SARS-CoV-2, children rarely develop severe illness. Differences in the host response to the virus at the primary infection site are among the proposed mechanisms.

**Objectives** To investigate the host response to SARS-CoV-2, respiratory syncytial virus (RSV), and influenza virus (IV) in the nasal mucosa in children and adults.

**Methods** Clinical outcomes and gene expression in the nasal mucosa were analyzed in 36 children hospitalized with SARS-CoV-2 infection, 24 children with RSV infection, 9 children with IV infection, 16 adults with mild to moderate SARS-CoV-2 infection, and 7 healthy pediatric and 13 healthy adult controls.

**Results** In both children and adults, infection with SARS-CoV-2 leads to an interferon response in the nasal mucosa. The magnitude of the interferon response correlated with the abundance of viral reads and was comparable between symptomatic children and adults infected with SARS-CoV-2 and symptomatic children infected with RSV and IV. Cell type deconvolution identified an increased abundance of immune cells in the samples from children and adults with a viral infection. Expression of *ACE2* and *TMPRSS2* – key entry factors for SARS-CoV-2 – did not correlate with age or presence or absence of viral infection.

**Conclusions** Our findings support the hypothesis that differences in the immune response to SARS-CoV-2 determine disease severity, independent of viral load and interferon response at the primary infection primary site.

Keywords
*   COVID-19
*   pneumonia
*   viral infections
*   interferons

## Introduction

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the current pandemic known as coronavirus disease 2019 (COVID-19). A notable feature of SARS-CoV-2 infection is the dramatic heterogeneity in clinical phenotypes. While most infected individuals remain asymptomatic or develop a very mild disease, some develop severe SARS-CoV-2 pneumonia and acute respiratory distress syndrome. These severe infections account for most of the public health and societal impacts of the COVID-19 pandemic and disproportionately affect the elderly while mostly sparing children (1–3). This is different from other common viral respiratory infections, including respiratory syncytial virus (RSV) and influenza virus (IV), which tend to be most severe in young children and elderly adults. Of the approximately 225,000 deaths from COVID-19 reported in the first eight months of the pandemic in the United States, only 189 occurred in children less than 18 years of age (4), despite at least 277,285 documented infections in children over the same period (5).

Identifying the mechanisms driving morbidity and mortality differences between children and adults infected with SARS-CoV-2 may provide essential insights for treating patients with COVID-19. Infection begins in the nasal mucosa (6, 7), where the SARS-CoV-2 spike protein binds angiotensin-converting enzyme 2 (ACE2) on the surface of epithelial cells. In concert with host proteases, principally transmembrane serine protease 2 (TMPRSS2), the virus enters the cell to begin replication (8–10). Multiple studies demonstrated that *ACE2* expression in the airways increases with age (11, 12), leading to lower infectivity and lower viral load in children. However, subsequent studies have demonstrated that SARS-CoV-2 viral load is similar in children and adults, suggesting the relative protection from COVID-19 in children is not due to age-related differences in infectivity (13–16). As there are important differences between the pediatric and adult immune systems, the host response to SARS-CoV-2 could explain the observed differences in n children and adults’ clinical phenotypes. The primary host response to viral infection is antiviral interferon signaling to limit viral replication, recruit immune cells, and clear infected cells (17). However, some viruses have developed strategies to evade this response, and SARS- CoV-2 is known to inhibit host interferon responses *in vitro (18, 19)*. This has led investigators to speculate that differences in antiviral interferon signaling at the primary infection site might be responsible for differences in the frequency of severe COVID-19 between children and adults. To date, no studies have compared the host response to the virus at the primary site of SARS-CoV-2 infection in children and adults.

This study aimed to assess how the host response to SARS-CoV-2 at the primary site of infection compares between children and adults, and how it may differ from the host response to other respiratory viruses, namely R.S.V. and IV, that are known to cause severe disease in children. Curettage and transcriptional profiling of the nasal mucosa can provide valuable insights into host-pathogen interactions at the initial infection site that may drive disease outcomes (20–23). Therefore, we performed R.N.A. sequencing of nasal mucosa samples from children and adults with SARS-CoV-2, children with R.S.V. or IV, and healthy children and adults. Local antiviral interferon signaling did not differ between children and adults infected with SARS-CoV-2 and was similar to children with severe R.S.V. infection. The magnitude of the interferon response only correlated with viral load and not the severity of illness. Our findings suggest that the local interferon response to SARS-CoV-2 in the upper respiratory tract does not differ between children and adults, arguing that factors other than the local antiviral response drive differences in the clinical phenotypes associated with SARS-CoV-2 infection.

## Methods

### Study population

Our study included the following six groups: healthy children, children infected with SARS- CoV-2, children infected with RSV, children infected with IV, healthy adults, and adults infected with SARS-CoV-2. Approval for this study was obtained from the Institutional Review Board at Ann & Robert H. Lurie Children’s Hospital of Chicago (Chicago, IL, U.S.A.). Informed consent was obtained from adult participants and parents/guardians of pediatric participants. The medical records of all pediatric patients admitted between December 2019, and November 2020 were screened for inclusion criteria and the absence of exclusion criteria (see below). We sequenced samples from healthy children recruited for a previous study as negative controls. These samples were obtained from healthy children attending a well-child check at Ann & Robert H. Lurie Children’s Hospital of Chicago Outpatient Center in May 2018. If the patient was eligible for the study, a study team member approached for consent. For SARS-CoV-2-infected children, both the child and accompanying adult were invited to participate.

### Inclusion criteria

Children were eligible for inclusion in the SARS-CoV-2-infected group if they had a positive polymerase chain reaction (PCR) test for SARS-CoV-2 at admission to the hospital or our drive through testing center (**Supplemental Figure 1A**). Symptoms were collected based on reports by participants or their guardians.

Adults accompanying a SARS-CoV-2-infected child or who were themselves admitted and cared for in our pediatric hospital with a positive PCR test for SARS-CoV-2 were eligible for inclusion (**Supplemental Figure 1B**). Polymerase chain reaction testing was not available for the adults accompanying a child prior to enrollment. RNA sequencing identified SARS-CoV-2 RNA in adult participants categorized with SARS-CoV-2. Adults without detectable SARS-CoV-2 via RNA sequencing were classified as healthy controls. Self-reported symptoms were collected from adults at the time of enrollment and via follow-up weekly phone calls for 3 weeks. Adults were asked about 11 of the most common symptoms reported with SARS-CoV-2 infection.

Participants were eligible for inclusion in the positive viral control group if they had a PCR confirming infection with RSV or IV, were admitted to the pediatric intensive care unit (PICU), and had symptoms of pulmonary parenchymal disease, including any of the following: 1) infiltrates, areas of atelectasis, or hyperinflation on chest radiograph, 2) rales, wheezing, retractions, tachypnea, or any form of respiratory distress on clinical exam, or 3) an oxygen requirement to maintain peripheral oxygen saturation levels above 95%.

Healthy pediatric controls were eligible for inclusion if they were less than 18 years of age and had no signs or symptoms of viral infection at the time of their scheduled well-child visit with their Primary Care Provider (**Supplemental Figure 1C**). All healthy pediatric controls were clinically stable and without apparent viral or bacterial infection when the nasal epithelial scrapings were obtained. Additionally, metagenomic analysis of RNA-seq data using the OneCodex platform did not detect any viral reads for 29 common human respiratory viruses. Healthy controls were followed for 2 weeks via electronic medical record to ensure no infections occurred within that time frame.

### Exclusion criteria

Subjects were not eligible for inclusion in the pediatric healthy control group if they had any significant medical co-morbidities, including pre-existing lung disorders, pre-existing cardiac disease, immunodeficiency, malignancy, and neurologic disorders, increasing the risk of aspiration or respiratory failure. Children with a clinical concern for multisystem inflammatory syndrome in children (MIS-C) were not eligible for inclusion in the SARS- CoV-2 cohort.

### Patient Enrollment

We enrolled 163 participants in the study (**Supplemental Figure 1**). Prior to sample collection, 10 participants withdrew consent. Eight participants were discharged before sample collection. After the isolation of RNA, 26 samples were excluded due to low RNA quality. Following RNA sequencing, an additional 14 samples were excluded: 1 for low-quality RNA-seq library, 3 for sex mismatch between biological sex indicated in the research records and expression of *XIST* or *RPS4Y1* genes, 3 for negative PCR testing results, 1 for multisystem inflammatory syndrome in children (MIS-C) diagnosis, 1 for unavailable patient research records, and 5 healthy controls for detection of viral reads for common respiratory viruses. A cohort of 105 participants: 7 healthy children, 36 children infected with SARS-CoV-2, 24 children infected with RSV, 9 children infected with IV, 13 healthy adults, and 16 adults infected with SARS-CoV-2 were used for the analysis.

### Clinical Data

Clinical data were collected and managed using REDCap electronic data capture tools hosted at the Northwestern University Clinical and Translational Sciences Institute (24). Hospital length of stay was determined from admission and discharge dates. Readmissions within 24 hours of discharge were considered a continuation of hospitalization. If a patient had multiple PICU stays during a single encounter, the PICU length of stay was cumulative. Advanced respiratory support was defined by the use of high flow nasal cannula (HFNC), non-invasive positive pressure ventilation (NIPPV), or invasive mechanical ventilation (IMV). Respiratory support and oxygen duration were determined by documentation on a calendar day of non-home advanced respiratory support or oxygen to exclude chronic baseline support. Pediatric acute respiratory distress syndrome was defined using the Pediatric Acute Lung Injury Consensus Conference criteria (25). Detection of a co-pathogen was defined as a positive viral PCR or any positive bacterial or fungal cultures.

### Nasal curettage and RNA isolation

Due to infection control procedures during the COVID-19 pandemic, nasal mucosa curettage was performed during routine care by health care providers caring for the patient. Providers received training just before sample collection by trained study personnel. Nasal epithelial cells were obtained by mucosal scrape biopsy of the inferior turbinate using a sterile plastic curette (Rhino-Pro curette, Arlington Scientific) as previously described (26) to obtain a predominantly epithelial cell population. A total of 3 single-pass scrapings per nares were performed per subject. The curette tips were cut and placed into RNase-free collection tubes containing 350 μL of RLT plus buffer (Qiagen, Germantown, MD, U.S.A.) supplemented with 2-mercaptoethanol, vortexed vigorously, and stored at −80°C. RNA was isolated using Qiagen RNeasy extraction kits (Qiagen).

### RNA sequencing and analysis

RNA quality and quantity were assessed using TapeStation 4200 High Sensitivity RNA tapes (Agilent), and RNA-seq libraries were prepared from 1 ng of total RNA using SMARTer Stranded Total RNA-seq Kit v2 (Takara Bio). After QC using TapeStation 4200 High Sensitivity DNA tapes (Agilent), dual indexed libraries were pooled and sequenced on a NextSeq 500 instrument (Illumina), 75 cycles, single-end, to an average sequencing depth of 5.88M reads.

FASTQ files were generated using bcl2fastq (Illumina). Viral RNA was detected using a custom hybrid genome prepared by joining FASTA, GFF, and GTF files for GRCh37.87, SARS-CoV-2 (NC\_045512.2), Influenza A/California/07/2009 (GCF_001343785.1), and RSV/S2 ts1C (GCF_000856445.1) genomes. An additional negative-strand transcript spanning the entirety of the SARS-CoV-2 genome was then added to the GTF and GFF files to enable detection of SARS-CoV-2 replication. To facilitate reproducible analysis, samples were processed using the publicly available nf-core/RNA-seq pipeline version 1.4.2 implemented in Nextflow 19.10.0 using Singularity 3.2.1-1 with the minimal command nextflow run nf-core/rnaseq -r 1.4.2 –singleEnd -profile singularity –reverseStranded --three_prime_clip_r2 3 (27–29). Briefly, lane-level reads were trimmed using trimGalore! 0.6.4 and aligned to the hybrid genome described above using STAR 2.6.1d (30). Gene-level assignment was then performed using featureCounts 1.6.4 (31). First, putative sample swaps were identified by comparing known patient sex with sex determined by *XIST* and *RPS4Y1* expression levels. No compromised samples were not included in the analysis.

### Weighted gene co-expression network analysis (WGCNA)

WGCNA was performed manually using WGCNA version 1.69 with default settings unless otherwise noted (32). Highly variable genes among all participant samples (n=105) and pediatric participant samples only (n=76) were included in analyses. To best capture patterns of co-regulation, a signed network was used. Using the pickSoftThreshold function, we empirically determined a soft threshold of 3 to best fit the network structure. A minimum module size of 30 genes was chosen. Module eigengenes were then related to patient and sample metadata using biweight midcorrelation. Module GO enrichment was then determined using GOrilla (33).

### Deconvolution of bulk RNA-seq signatures

Deconvolution of bulk RNA-seq signatures was performed using AutoGeneS v1.0.3 (34). For details of the functions used and their parameters, see the code available at [https://github.com/NUPulmonary/2021_Koch](https://github.com/NUPulmonary/2021_Koch). Briefly, we used a single-cell RNA-seq dataset from Ordovas-Montanes et al. (35), which contains nasal epithelial and immune cells, to train the AutoGeneS model. AutoGeneS selects cell type-specific genes by both minimizing the correlation between and maximizing the distance between the cell type-specific average gene expression profiles. The model was then applied to bulk RNA-seq data to estimate specific cell types’ proportions using regression.

### Statistical Analysis

Categorical data were summarized by percentages and compared using Fisher’s exact test. Continuous, non-parametric data were summarized using medians with interquartile ranges and compared using Mann-Whitney U or Kruskal-Wallis testing. The Benjamini-Hochberg procedure was used for false discovery rate correction when performing multiple comparisons. Data were considered statistically significant at P-adjusted<0.05. Statistical analyses were performed using Graphpad Prism version 9.0.0 and R version 4.0.3. Data were visualized using ggplot2 version 3.3.1 and Graphpad Prism.

### Data and code availability

Processed RNA-seq counts are available as Online Supplemental Table 1. Raw data is in the process of being deposited to dbGaP/SRA. Code is available via [https://github.com/NUPulmonary/2021\_Koch](https://github.com/NUPulmonary/2021_Koch).

## Results

### Description of the cohort

The demographics of the pediatric and adult cohorts are summarized in **Table 1** and **Supplemental Table 1**. When considering only children, patients with RSV were younger, and a larger proportion was previously healthy than those with SARS-CoV-2 (**Table 1**). The ages and proportion of individuals with chronic medical conditions in the uninfected adult group were no different from the SARS-CoV-2 infected adult group (**Supplemental Table 1**).

View this table:
[Table 1:](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/T1)

Table 1: 
Description of Pediatric Cohort. Adjusted p-value <0.05 were considered statistically significant, calculated by Kruskal-Wallis test with Benjamini-Hochberg FDR correction.

**Figure 1A** summarizes the clinical course of the hospitalized pediatric participants infected with IV, RSV, or SARS-CoV-2. There was no difference in the time from symptom onset to nasal sample collection between the groups (**Figure 1B**). Consistent with their inclusion as positive controls for moderate and severe viral respiratory tract infections, children with RSV were hospitalized longer and received higher respiratory support levels than children with SARS-CoV-2 (**Figure 1C, D, Table 1**). Similarly, children with IV received increased respiratory support than children with SARS-CoV-2 (**Figure 1D, Table 1**). Of the 36 children with SARS-CoV-2, 35 were admitted to the hospital, and 8 required ICU level care. Of those 8 participants, seven were admitted with respiratory symptoms requiring non-invasive respiratory support, while 1 was admitted with severe diabetic ketoacidosis in the setting of fever. None of the SARS-CoV-2-positive children were intubated. However, one patient with a history of obesity (BMI 37.2) met the criteria for pediatric acute respiratory distress syndrome (ARDS) while supported on full-face bilevel positive airway pressure. The age, sex, race, and ethnicity of the children with IV and the uninfected children did not differ from those with SARS-CoV-2 (**Supplemental Figure 1D, F-I**). A larger proportion of children with SARS-CoV-2 identified as Hispanic compared to RSV (**Supplemental Figure 1H**). The sex, race, and ethnicity were not different between uninfected and infected adult groups (**Supplemental Figure 1E-I, Supplemental Table 1**).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F1)

Figure 1. Demographics of the cohort.
**A**. Clinical course of pediatric participants. Timing of symptom onset (diamond), admission (square), sample collection (triangle), discharge (circle), hospital stay (thin line), and duration of advanced (HFNC, NIPPV, IMV) respiratory support (bold line) of hospitalized children with SARS-CoV-2 (n=35), RSV (n=25), and IV (n=8). **B**. Self or parent/guardian reported time since symptom onset at the time of sample collection was similar between all participant groups. **C**. Children with RSV infection had a longer hospital length of stay than children with SARS-CoV-2 infection. Pairwise comparisons of medians performed using the Mann-Whitney U test. **D**. Higher proportions of children with IV and RSV infection required HFNC, NIPPV, and IMV at peak illness severity when compared to children with SARS-CoV-2 infection. Proportions compared using Fisher’s Exact Test. All p-values were adjusted using FDR correction. Differences were not significant (p-adjusted>0.05) unless noted. RSV=respiratory syncytial virus, IV=influenza virus, HFNC=high flow nasal cannula, NIPPV=non-invasive positive pressure ventilation, IMV=intermittent mechanical ventilation, FDR=false discovery rate.

The majority of the adult participants with SARS-CoV-2 reported mild illness (**Supplemental Table 2**). There was no difference in the time from symptom onset to nasal sample collection between adults and all pediatric groups (**Figure 1B**). Two adult participants (age 19 and 20 years) were admitted to our pediatric hospital at the time of enrollment. An additional adult participant reported admission to an outside hospital after enrollment in our study. The two patients admitted to our hospital had respiratory illnesses requiring admission to the PICU for non-invasive positive pressure ventilation with supplemental oxygen. Neither met the ARDS criteria. The participant admitted to an outside hospital reported an admission duration of 3 days and did not require intensive care unit level care.

The symptoms reported by participants or the parent/guardian for children with SARS- CoV-2 infection are summarized in **Supplemental Table 1**. Of the 16 adults with SARS- CoV-2 infection (see methods), 7 (44%) were asymptomatic, and the other 9 (56%) reported symptoms. The majority of children with SARS-CoV-2 infection experienced respiratory symptoms (57%). Other symptoms included fever only (14%), gastrointestinal symptoms (9%), non-respiratory viral syndromes (6%), and asymptomatic (14%) (**Table 1**). The asymptomatic children with SARS-CoV-2 were detected on routine screening at admission.

### SARS-CoV-2 viral load and expression of SARS-CoV-2 entry factors in the nasal mucosa does not correlate with age

To characterize our participant sample data and confirm clinical diagnoses, we queried bulk RNA-seq data obtained from curettage of nasal mucosa for the presence of SARS- CoV-2, RSV, or IV reads. Viral reads for SARS-CoV-2 were detected in 28 of 36 samples collected from children diagnosed with SARS-CoV-2 infection by PCR performed in the hospital diagnostic laboratory (**Figure 2A**). Of the 24 children diagnosed with RSV, all had detectable RSV reads in RNA-seq data. Six out of 9 children with IV had detectable Influenza A viral reads in their nasal mucosa. Of the 29 adults accompanying children with SARS-CoV-2 infection, 14 had detectable SARS-CoV-2 reads, and 2 young adults admitted to the pediatric hospital with a positive SARS-CoV-2 PCR test did not have detectable viral reads. Nine samples were obtained from healthy children. RNA-seq did not detect any SARS-CoV-2, RSV, or IV reads in these samples; however, a more in-depth analysis of these samples using the TWIST viral panel as implemented in the OneCodex platform demonstrated the presence of enterovirus and parainfluenza virus in two of the samples. These samples were excluded from subsequent analysis, leaving 7 subjects in the healthy control group (**Supplemental Figure 1C**). In agreement with previous reports (13–16), we found that SARS-CoV-2 viral load was similar between children and adults and did not correlate with *ACE2* expression (**Figure 2B, C**). Several studies have reported that the expression of *ACE2* and *TMPRSS2* – two key entry factors for SARS-CoV-2 – increases with age (11, 12). While some reports suggested that interferons produced in response to viral infection can lead to increased expression of *ACE2*, potentially contributing to the spread of infection (7, 36), other reports have found no correlation or the opposite effect (37–39). We found that the expression of *ACE2* and *TMPRSS2* was not different between healthy children and adults. Moreover, none of the viral infections led to increased expression of *ACE2* and *TMPRSS2* (**Figure 2D**). Similarly, *ACE2* and *TMPRSS2* expression did not correlate with age (**Figure 2E**). Taken together, these data demonstrate that transcriptomic data confirmed clinical diagnoses for SARS-CoV-2, IV, and RSV infections. Moreover, in contrast to previous studies, we did not find a significant correlation between SARS-CoV-2 entry factors and age in our cohort.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F2)

Figure 2. SARS-CoV-2 viral load and expression of SARS-CoV-2 entry factors in the nasal mucosa does not correlate with age
We compared the nasal transcriptome of 36 children with SARS-CoV-2 infection with 7 healthy children, 24 children with RSV infection, 9 children with IV infection, 16 adults with SARS-CoV-2 infection, and 13 healthy adults. **A**. Detection of viral reads for SARS-CoV- 2, influenza A/California/07/2009, and RSV/S2 ts1C in RNA-seq data for all participant samples (n=105). Heatmap of viral genes for SARS-CoV-2, IV, and RSV. Log-transformed pseudocounts (CPM+1) are shown. **B**. Normalized reads were averaged for all 11 SARS-CoV-2 viral genes. Average SARS-CoV-2 reads did not differ between children (yellow, circle) and adults (purple, triangle). Log-transformed pseudocounts (CPM+1) are shown. **C**. Averaged SARS-CoV-2 viral reads did not correlate with *ACE2* expression. Log-transformed pseudocounts (CPM+1) are shown. Pearson correlation coefficient (r) and p-value are shown. **D**. Normalized expression of *ACE2* and *TMPRSS2* was similar across all participant groups; differences not significant, Kruskal-Wallis testing with Benjamini-Hochberg FDR correction, all adjusted p-values >0.05. Log-transformed pseudocounts (CPM+1) are shown. **E**. *ACE2* and *TMPRSS2* expression did not correlate with age. Pearson correlation coefficient (r) and p-value are shown. Log-transformed pseudocounts (CPM+1) are shown.

### Cell-type deconvolution of transcriptomic signatures identifies an influx of immune cells into the nasal mucosa upon viral infection

Bulk RNA-seq data obtained from complex tissue samples composed of multiple cell types can obscure disease-specific signatures as RNA molecules from different cell types are averaged. Accordingly, our initial exploratory analysis did not reveal a distinct structure or pattern within our dataset (**Supplemental Figure 3A, B**). To address the potential heterogeneity in cell-type composition in our samples, we performed *in silico* cell type deconvolution of our bulk RNA-seq data using cell type-specific reference data from the nasal mucosa (35). Hierarchical clustering of predicted cell type abundance distinguished 6 clusters of samples (**Figure 3A, Online Supplemental Table 2**). Four clusters were characterized by an abundance of immune cells (clusters 1–4), while two clusters were enriched with secretory and ciliated epithelial cells (clusters 5 and 6). In line with previous reports (40, 41), deconvolution analysis demonstrated that most nasal mucosa samples from healthy participants were enriched for epithelial cells, with all of the samples from healthy children dominated by a ciliated cell signature (**Figure 3A, B**). In contrast, samples from children with RSV infection were equally divided between the secretory cell cluster and the immune cell clusters, consistent with previous reports of replacement of ciliated cells with secretory cells and the recruitment of neutrophils in RSV infection (42, 43). Increased immune cell abundance was observed in most samples from children infected with IV and SARS-CoV-2 (**Figure 3A, B, Supplemental Figure 3C**). The time from onset of symptoms to the nasal mucosa samples’ collection did not differ between epithelial and immune dominant subsets (**Figure 3C**). In children and adults infected with SARS-CoV-2, viral reads levels were increased in immune-enriched samples compared to epithelial-enriched samples, suggesting increased viral load is associated with immune cell recruitment to the nasal mucosa (**Figure 3D**). To confirm that differences in cell-type composition were not obscuring differences in *ACE2* expression in children and adults, we analyzed *ACE2* expression in samples from the epithelial cell clusters and the immune cell clusters separately. There remained no difference in *ACE2* expression between children and adults **(Figure 3E**). Interestingly, deconvolution analysis of the two samples that were excluded from the healthy children group because of the incidental detection of enterovirus and parainfluenza virus found that these samples had an increased percentage of neutrophils, thus supporting transcriptomic profiling of nasal mucosa as a sensitive tool for assessment of viral infections (**Online Supplemental Table 2**). Taken together, these data demonstrate that samples obtained from the nasal mucosa of healthy and virus-infected participants are characterized by distinct epithelial-or immune-cell abundance profiles. Specifically, samples obtained from participants infected with SARS-CoV-2, IV, or RSV are characterized by an increased abundance of immune cells, whereas healthy participant samples are enriched for epithelial cells.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F3)

Figure 3. Cell type deconvolution of the transcriptomic signatures in the nasal mucosa demonstrates an influx of immune cells after viral infection.
**A**. Heatmap demonstrating the relative abundance of cell types as determined by *in silico* deconvolution of bulk RNA-seq gene expression profiles (see Methods). Hierarchical clustering on samples from 105 participants identified 6 clusters enriched for immune or epithelial cells. Z-score normalized relative abundance values are shown. **B**. Relative abundance of immune and epithelial cell types across all six participant groups. Participants with a viral infection had a significantly higher relative abundance of immune cells than healthy participants. Adjusted p-values calculated using the Kruskal-Wallis test with Benjamini-Hochberg FDR correction are shown. Adjusted p-values <0.05 were considered statistically significant. **C**. Number of days between symptom onset and day of sample collection were similar for immune-enriched and epithelial-enriched samples. n.s.= not significant, by Mann-Whitney U test. **D**. SARS-CoV-2 viral load, defined as average normalized reads for 11 SARS-CoV-2 viral genes, is significantly higher in immune-enriched clusters in comparison to epithelium-enriched clusters. Log-transformed pseudocounts (CPM+1) are shown. The adjusted p-value calculated by the Mann-Whitney U test is shown. **E**. Normalized expression of *ACE2* in epithelial- and immune-enriched samples is similar across all participant groups; differences are not significant, Kruskal-Wallis testing with Benjamini-Hochberg FDR correction, all adjusted p-values >0.05. Log-transformed pseudocounts (CPM+1) are shown.

### Viral infection is associated with increased expression of interferon and innate immune response genes in the nasal mucosa

To investigate whether specific gene expression patterns correlate with clinical parameters, we performed a weighted gene co-expression network analysis (WGCNA). Analysis of all participant samples (n=105, both children and adults) identified 7 modules, 3 of which significantly correlated with one or more clinical parameters **(Figure 4A**). Module 7 positively correlates with age and does not correlate with any other parameter, suggesting it contains age-related genes not associated with a viral infection. Enrichment analysis revealed that this age-related module enriches for genes involved in GO biological processes such as *RNA processing* and includes several small nucleolar RNA genes (snoRNAs) such as *SNORA49* and *SNORD8* (Online Supplemental Table 3), which have been proposed as aging markers (44)(**Figure 4B**). Modules 3 and 4 positively correlated with the number of SARS-CoV-2 viral reads. Module 3, which positively correlated with SARS-CoV-2 and IV reads, enriched for genes involved in *response to virus* and *type I interferon signaling pathway* (**Figure 4B**). Module 4 also correlated with RSV reads and enriched for genes involved in GO biological processes associated with *neutrophil activation, chemotaxis*, and *leukocyte differentiation*.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F4)

Figure 4. Transcriptional signatures associated with viral infection enrich for inflammatory and interferon response genes in children and adults.
**A**. 2170 HVGs for all participant samples (n=105) were analyzed using WGCNA. Module-trait relationships were identified for 7 modules. Modules 1-7 contain 33, 199, 248, 992, 375, 32, and 51 genes, respectively. 240 genes that did not fall into any module were assigned to the null module. **B**. GO biological process enrichment for modules correlating with SARS-CoV-2 viral reads. **C**. 2024 HVGs identified among pediatric participant samples only (n=76) were analyzed using WGCNA. Module-trait relationships were identified for 7 modules. Modules 1-7 contain 43, 385, 823, 92, 200, 204 and 46 genes, respectively. 231 genes that did not fall into any module were assigned to the null module. **D**. GO biological process enrichment for modules correlating with SARS-CoV-2 reads. **E**. Scatterplot of correlation between Gene Significance (GS) for SARS-CoV-2 reads vs. Module Membership (MM) in modules 3,4 and 5. Pearson correlation coefficients (r) and p-values are shown. P-values <0.05 were considered statistically significant. Linear regression line (blue) with a 95% confidence interval (grey) are shown. HVGs - highly variable genes, WGCNA - Weighted Gene Co-Expression Network Analysis. Coll\_O2 – On supplemental oxygen at time of collection. Coll\_none – No respiratory support at time of sample collection. Coll\_cannula – Nasal cannula/simple mask at time of sample collection. Coll\_ppv – Positive pressure ventilation (PPV) at time of sample collection. Peak\_o2 – Received supplemental oxygen. Peak\_none – No respiratory support Peak\_cannula – Nasal cannula/simple mask was highest respiratory support during course of illness. Peak\_ppv – PPV was highest respiratory support during course of illness. Sx2coll – Time between symptom onset and collection (days). Icu_los – ICU length of stay. Hosp_los – Hospital length of stay. Supp_dur – Duration of support (days). Oxy\_duration – Duration of oxygen therapy.

We next focused only on pediatric samples to identify distinct modules correlated with specific viral infections and clinical outcomes. WGCNA identified 7 modules within the pediatric cohort (n=76)(**Figure 4C**). Modules 3, 4, and 5 correlated with SARS-CoV-2 reads (**Figure 4C**). Module 3 also correlated with RSV reads and enriched for genes involved in the GO biological processes *inflammatory response, regulation of IL-6 production*, and *T cell activation* (**Figure 4D**). Module 5 correlated with both SARS-CoV- 2 and IV reads and enriched for genes involved in *response to virus, type I interferon signaling*, and chemokines, including *CXCL10* (**Figure 4E)**. Module 4 is the only module unique to SARS-CoV-2 reads. The 92 genes that make up this module did not enrich for specific GO biological processes but included the antiviral transcription factor *IRF8*, the chemokines *CXCL9* and C*XCL11*, and the T cell marker CD4 (**Figure 4E**). There was no significant correlation between any of the 7 modules and severity of illness or outcomes. Taken together, these data demonstrate that in participants with SARS-CoV-2, RSV, or IV infection, the nasal mucosa is characterized by the expression of genes associated with innate immune and interferon response.

### Local interferon response is independent of age and virus

Since WGCNA identified antiviral interferon signaling as a shared feature of viral infection with SARS-CoV-2, IV, or RSV in our cohort, we next focused on this pathway. We found that gene lists for type I and type II interferon responses largely overlap. Therefore, we compiled a list of IFN-response genes using hallmark pathways for Interferon alpha and gamma responses (45, 46). Comparison of normalized expression across all participant groups revealed an increased expression of interferon response genes in patients with RSV, IV, and SARS-CoV-2 compared to healthy controls that did not correlate with age (**Figure 5A-C**). Specifically, we did not observe differences in the interferon response between children and adults infected with SARS-CoV-2 (**Figure 5B**). While neither age nor type of viral infection correlated with the interferon response, we found a significant correlation with the viral load measured by an average number of viral reads (**Figure 5D**). Consistent with the lack of correlation between viral load, the severity of illness, or the time from symptom onset, the interferon response did not correlate with clinical outcomes such as duration of respiratory support and ICU length of stay (**Supplemental Figure 5A, B**). Interestingly, although we detected an interferon response signature in the majority of virus-infected participants, we did not detect the expression of secreted type I, type II, or type III interferons at the transcriptional level in a majority of patients with RSV or SARS- CoV-2 infection (**Figure 5E**). Secreted interferons were primarily detected in participants with IV infection, consistent with this cohort having the highest average expression of interferon response genes overall (**Figure 5B, E**). This suggests that the nasal mucosa is not a source of secreted interferons in RSV and SARS-CoV-2 infection, at least when these infections are brought to medical attention.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F5)

Figure 5. Local interferon response is independent of age and virus.
**A**. Heatmap of z-score normalized expression (CPM) of 222 Type I and II IFN response genes. Hierarchical clustering on genes (rows) did not identify distinct clusters. Samples within SARS-CoV-2, IV, and RSV groups are ordered by ascending viral reads for SARS-CoV- 2, IV, or RSV, respectively. **B**. Average normalized expression of IFN response genes (combined Hallmark pathway lists for IFNA and IFNG response) by participant group. Asterisks indicate statistical significance (adjusted p-value <0.05, Kruskal-Wallis test with Benjamini-Hochberg FDR correction). **C**. Scatter plot of correlation between age (years) and average IFN response reads. Log-transformed CPM counts are shown. Pearson correlation coefficient (r) and p-value are shown. **D**. Scatter plot of correlation between average IFN response reads (log2CPM) and SARS-CoV-2, RSV, or IV viral loads, respectively. Pearson r-values and p-values are shown. **E**. Heatmap of z-score normalized expression (CPM) for 19 secreted Type I, II, and III interferons.

### Increased expression of innate immune response genes in children with high SARS-CoV-2 viral load

Although IFN response was not different between children and adults infected with SARS- CoV-2 (**Figure 5B**), we wanted to determine if the expression of genes in other pathways was different. Therefore, we next focused only on children and adults with high viral loads (average SARS-CoV-2 viral read > 6 CPM), indicating active infection. We found no difference in the magnitude of the interferon response between children and adults with a high average viral load (**Figure 6A**). The cellular composition was also similar between the groups, allowing us to compare gene expression differences without cellular composition bias (**Figure 6B**). We compared children and adults with high SARS-CoV-2 viral load and identified 576 differentially expressed genes (**Figure 6C**). Genes upregulated in children compared to adults were enriched for innate immune processes, including *cellular response to IL-1* and *inflammatory response* (**Figure 6C, D**). Inflammatory cytokines, their receptors, and chemokines, like *IL18, IL1RL2, CXCL3* and *CCL2*, were among the genes increased in children (**Figure 6C, E**). Participants with low to no SARS-CoV-2 viral load had fewer differentially expressed genes between children and adults and the enrichment analysis did not identify any significant GO biological processes (**Supplemental Figures 6A-B**). Furthermore, we did not observe difference in the IFN response or ACE2 between children and adults with low viral reads (**Supplemental Figures 6B**). There were also no differences in immune or epithelial cellular abundance and no differences in clinical outcomes between low and high viral counts groups (**Supplemental Figures 6C-D**). This data suggests that SARS-CoV-2- infected children with a high viral load have a similar interferon response but an increased innate immune response in the nasal mucosa compared to SARS-CoV-2-infected adults.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F6)

Figure 6. Increased expression of innate immune response genes in children with high SARS-CoV-2 viral load.
**A**. The average normalized IFN response gene expression was similar in children and adults with high viral load (>6 CPM average). n.s. = not statistically significant by Kruskal-Wallis test with Benjamini-Hochberg FDR correction. **B**. Relative cellular abundance. n.s. = not significant by Kruskal-Wallis test with Benjamini-Hochberg FDR correction. **C**. MA plot of differentially expressed genes (edgeR p-value <0.05 and log2FC > |1.0| between adults and children with high SARS-CoV-2 viral load (> 6 CPM average). Number of genes upregulated (red) and downregulated (blue) in adults compared to children are shown. **D**. Genes upregulated in children enriched for innate immune-related GO biological processes. **E**. Scatter plots of normalized expression for *IL18, IL1RL2*, and *CXCL3* are shown. Line indicates median. Log-transformed CPM counts are shown.

## Discussion

The majority of children infected with SARS-CoV-2 will at most experience mild symptoms similar to common respiratory viral infections. The mechanisms underlying this relative protection from a severe infection in children compared to adults are unclear. Differences in expression of the SARS-CoV-2 entry factors or response to the virus at the primary site of infection between children and adults have been proposed as potential mechanisms (47, 48); however, these hypotheses have not been directly tested to date. Here we present the first study that directly compares host response to SARS-CoV-2 at the primary site of the infection between children and adults. In our cohort of children infected with SARS-CoV-2, the local response in the nasal mucosa was largely the same as non-hospitalized adults infected with SARS-CoV-2 and children with severe RSV and IV infections. This similarity in response to the virus, particularly the interferon response, suggests that age-related differences in SARS-CoV-2 infection outcomes are independent of the primary response to the virus in the nasal mucosa.

Age-related differences in expression of *ACE2* and *TMPRSS2* have been proposed as modulators of clinical severity of SARS-CoV-2 infections. ACE2 is the cell surface receptor bound by the SARS-CoV-2 spike protein, and TMPRSS2 is a host protease that cleaves the spike protein to allow for efficient binding and cell entry (8). Previous studies have demonstrated an association between age and the expression of *ACE2* and *TMPRSS2*, suggesting that low levels of expression of *ACE2* and *TMPRSS2* in pediatric airway epithelial cells could protect from severe infection (11, 12, 36). In our cohort, we did not observe a difference in the expression of *ACE2* or *TMPRSS2* between children and adults. While some reports suggested that *ACE2* expression may increase in response to interferon stimulation (7, 36), we saw no increase in the expression of *ACE2* and *TMPRSS2* expression during active SARS-CoV-2 or other viral infections, supporting previous reports that interferon signaling does not lead to increased *ACE2* expression (37–39).

Our study’s unique strength is the inclusion of the appropriate controls: healthy children and adults and children infected with other respiratory viruses. This allowed us to report our findings on SARS-CoV-2 infection in the context of RSV and IV infections, which are known to induce severe illness requiring hospitalization and respiratory support. We found a high degree of similarity in the local response to SARS-CoV-2, RSV, and IV infections. This agrees with Mick and colleagues’ report, who performed metatranscriptomics on nasopharyngeal swabs and found substantial overlap between SARS-CoV-2 and other acute respiratory infections in adults (49). Thus, our findings suggest that other factors, such as systemic immune responses, may determine the differences in disease severity between children and adults.

Given that the expression of SARS-CoV-2 entry factors and the local response at the primary site of infection does not differ between children and adults, what factors could explain relative protection in children and unpredictable susceptibility in elderly adults infected with SARS-CoV-2? As humans have numerous encounters with betacoronaviruses, including seasonal coronaviruses, they form adaptive immunity (50). Thus, even before the COVID-19 pandemic, a substantial number of individuals in the population had antibodies and T cells that recognize betacoronaviruses and can cross-react with SARS-CoV-2. Such cross-reactive memory T cells have been reported for SARS-CoV (51) and recently for SARS-CoV-2 (52–56). This observation initially led to speculation that these cross-reactive antibodies and T cells could play a protective role in children and young adults and that vanishing adaptive immunity would make elderly individuals more susceptible to SARS-CoV-2. However, in another study, a larger number of SARS-CoV-2 cross-reactive memory T cells was observed independently in the elderly and patients with severe compared with mild COVID-19 (56). Importantly, these cross-reactive T cells exhibited lower avidity and reduced antiviral responses when stimulated with SARS-CoV-2 peptides compared with T cells from patients who recovered from COVID-19. Thus, instead of being protective, these cross-reactive low avidity T cells, which are less efficient at virus elimination, could contribute to prolonged production of inflammatory signaling. This hypothesis is attractive as it explains the differences in the epidemiology of SARS-CoV-2, which disproportionately affects elderly individuals. Children and young adults, who likely have fewer lifetime coronavirus exposures, often have only mild symptoms, despite having viral loads in the upper airways comparable to adults (3, 14). While multiple studies have investigated adaptive immune response in adults infected with COVID-19, similar systematic studies in the pediatric population are lacking.

Our study has several limitations. First, our institution utilized two different PCR platforms for the detection of SARS-CoV-2 during the study period. Due to differences in methodology, viral loads cannot be directly compared between the platforms. Moreover, PCR testing was not available for the adults accompanying children infected with SARS- CoV-2 during the study period. To overcome this limitation, we used dual RNA-seq, which allows for simultaneous host transcriptomic analysis and pathogen characterization (21). By creating custom hybrid genomes, we were able to identify and quantify viral transcripts in our samples as a surrogate for viral load, thus, enabling comparison across all groups in our study. Second, while our study includes children infected with other respiratory viruses – RSV and IV, which serve as proper controls to SARS-CoV-2-infected children, it primarily includes hospitalized children with SARS-CoV-2, representing a minority of pediatric SARS-CoV-2 infections. However, 14% of our pediatric cohort were asymptomatic with incidental SARS-CoV-2 diagnoses. The local host response in these asymptomatic children’s nasal mucosa could not be distinguished from more severely affected pediatric participants. Similarly, our study only included two hospitalized young adults with moderate SARS-CoV-2 infection and did not include any severely ill hospitalized adults requiring intubation and mechanical ventilation. Therefore, hospitalized adults with severe SARS-CoV-2 may have a different local response at the primary site of infection. However, the high degree of similarity in the local response between different respiratory viruses in our cohort suggests that other factors, such as the systemic immune response, determine disease severity in SARS-CoV-2 infection. Supporting this interpretation, recent studies using single-cell RNA-seq profiling of nasopharyngeal swabs demonstrated activation of an interferon-response signature in epithelial cells and influx of immune cells in hospitalized adults with severe COVID-19 (57, 58). Third, our transcriptomic analysis was performed on a biopsy from a complex tissue containing multiple epithelial and immune cell types (40). Thus, our results are affected by the changes in the abundance of the specific cell types, either due to age- or disease-related differences in tissue composition or operator sampling bias. We applied *in silico* tissue deconvolution (34) to account for differences in tissue composition, which allowed us to perform a more accurate comparison between samples with similar composition. This comparison demonstrated that among subjects with high SARS-CoV- 2 viral load, children exhibit increased activation of the innate immune system on a per-cell level. Ultimately, the use of techniques that enable analysis at the single-cell level, such as single-cell RNA-seq, will allow a more precise assessment of cell-type-specific changes associated with SARS-CoV-2 or other viral infections. Recently, single-cell RNA- seq has been applied to understand virus-host interactions in children infected with IV (59), and similar efforts are underway in SARS-CoV-2 (60).

In summary, we present evidence that the response to SARS-CoV-2 infection at the primary site of infection is similar between children and adults and comparable in its magnitude to other common respiratory viruses. Our findings suggest that factors outside of the primary infection site are likely drivers of the differences in severity of SARS-CoV- 2 infection between children and adults.

## Data Availability

Processed RNA-seq counts are available as Online Supplemental Table 1. Raw data is in the process of being deposited to dbGaP/SRA. Code is available via https://github.com/NUPulmonary/2021_Koch.

## Figure and Table Legends

![Supplemental Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F7.medium.gif)

[Supplemental Figure 1.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F7)

Supplemental Figure 1. Flow chart of study subject recruitment and sample analysis.
**A**. Enrollment of children with RSV, IV, and SARS-CoV-2 infection from December 2019 to November 2020. **B**. Enrollment of adults accompanying children with SARS-CoV-2 infection and young adults with SARS-CoV-2 infection admitted to our hospital from March 2020 to November 2020. **C**. Enrollment of healthy children in a previous study in May 2018. **D**. Children with RSV infection were younger than those with SARS-CoV-2 and IV infections. **E**. Distribution of ages between adults infected with SARS-CoV-2 and healthy adults was similar. **F**. Proportion of females (blue) and males (green) were similar between all groups. **G**. Proportions of self or parent/guardian reported race differed between children with SARS-CoV-2 and RSV infections. **H**. A larger proportion of children with SARS-CoV-2 infection identified as Hispanic or Latino when compared to children with RSV infection. **I**. Among children requiring advanced respiratory support, there was no difference in the duration of respiratory support. Pairwise comparisons of medians calculated using the Mann-Whitney U test. Proportions compared using Fisher’s Exact Test. All p-values were adjusted using FDR correction. Differences were not significant (p-adjusted>0.05) unless noted. RSV=respiratory syncytial virus, IV=influenza virus, HFNC=high flow nasal cannula, NIPPV=non-invasive positive pressure ventilation, IMV=intermittent mechanical ventilation, FDR=false discovery rate.

![Supplemental Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F8.medium.gif)

[Supplemental Figure 3.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F8)

Supplemental Figure 3. 
**A**. Principal component analysis (PCA) of normalized counts data (CPM) from all participant samples (n=105) demonstrates a lack of structure related to age or diagnosis. **B**. Heatmap of Pearson correlation of mean gene expression across all samples, hierarchical clustering performed on rows and columns. Participant samples do not cluster by age or diagnosis, except for the healthy children group. **C**. Relative abundance of immune and epithelial cell types across all six participant groups. Adjusted p-values were calculated by Kruskal-Wallis with Benjamini-Hochberg FDR correction.

![Supplemental Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F9.medium.gif)

[Supplemental Figure 5.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F9)

Supplemental Figure 5. 
**A**. Time from symptom onset to day of collection in children and adults with SARS-CoV-2, RSV, or IV infection. Pearson correlation coefficient (r) and p- value are shown. **B**. Scatter plot of correlation between support duration (days) or ICU length of stay compared to average IFN response. Pearson correlation coefficient (r) and p-value are shown.

![Supplemental Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/28/2021.01.26.21250269/F10.medium.gif)

[Supplemental Figure 6.](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/F10)

Supplemental Figure 6. 
**A**. MA plot of differentially expressed genes (edgeR p-value <0.05 and log2FC > |1.0| between adults and children with low to no SARS-CoV-2 viral load (< 6 CPM average). Number of genes upregulated (red) and downregulated (blue) in adults compared to children are shown. **B**. SARS-CoV-2 viral load, defined as average normalized reads for 11 SARS-CoV-2 viral genes, is similar between children and adults in the low and high viral reads groups. Log-transformed pseudocounts (CPM+1) are shown. The average normalized IFN response gene expression was similar between children and adults with low (< 6 CPM average) or high viral load (> 6 CPM average). *ACE2* expression did not differ between children and adults with low (< 6 CPM average) or high viral load (> 6 CPM average). n.s. = not statistically significant by Kruskal-Wallis test with Benjamini-Hochberg FDR correction. **C**. Relative cellular abundance. n.s. = not significant by Kruskal-Wallis test with Benjamini-Hochberg FDR correction. **D**. No difference in clinical characteristics of pediatric participants with low (< 6 CPM average) or high viral load (> 6 CPM average).

View this table:
[Supplemental Table 1:](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/T2)

Supplemental Table 1: 
Description of Adult Cohort. Adjusted p-value <0.05 were considered statistically significant, calculated by Kruskal-Wallis test with Benjamini-Hochberg FDR correction.

View this table:
[Supplemental Table 2:](http://medrxiv.org/content/early/2021/01/28/2021.01.26.21250269/T3)

Supplemental Table 2: 
Self or parent/guardian reported SARS-CoV-2 symptoms.

## Acknowledgments and funding

This research was supported in part through the computational resources and staff contributions provided by the Genomics Compute Cluster, which is jointly supported by the Feinberg School of Medicine, the Center for Genetic Medicine, and Feinberg’s Department of Biochemistry and Molecular Genetics, the Office of the Provost, the Office for Research, and Northwestern Information Technology. The Genomics Compute Cluster is part of Quest, Northwestern University’s high-performance computing facility, to advance research in genomics.

The authors thank all the participants for enrolling in this study, The authors Rogan Grant for his help with the RNA-seq pipeline and Benjamin Singer and G.R.S. Budinger for their thoughtful discussions.

Andrew D. Prigge was supported by the Gorter Family Foundation. Alexander V. Misharin was supported by N.I.H. grants U19AI135964, P01AG049665, R01HL153312, and NUCATS COVID-19 Rapid Response Grant. Karen M. Ridge was supported by N.I.H. grants P01AG049665 and P01GM096971. Bria M. Coates was supported by funds provided by the Manne Research Institute COVID-19 Springboard Exploratory Research Award and N.I.H. grant K08HL143127. Taylor Heald-Sargent and L. Nelson Sanchez-Pinto were supported by funds provided by the Manne Research Institute COVID-19 Springboard Exploratory Research Award.

The authors do not declare a conflict of interest.

## Footnotes

*   # Share last authorship

*   **Funding sources:** Manne Research Institute, NHLBI K08HL143127

*   Received January 26, 2021.
*   Revision received January 26, 2021.
*   Accepted January 28, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Castagnoli R, Votto M, Licari A, Brambilla I, Bruno R, Perlini S, Rovida F, Baldanti F, Marseglia GL. Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV- 2) Infection in Children and Adolescents: A Systematic Review. JAMA Pediatr 2020;174:882–889.
    
    

2.  2.Ludvigsson JF. Systematic review of COVID-19 in children shows milder cases and a better prognosis than adults. Acta Paediatr 2020;109:1088–1095.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/apa.15270&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

3.  3.Dong Y, Mo X, Hu Y, Qi X, Jiang F, Jiang Z, Tong S. Epidemiology of COVID-19 Among Children in China. Pediatrics 2020;145.:
    
    

4.  4.CDC. COVID-19 cases, deaths, and trends in the U.S. 2020;at <[https://covid.cdc.gov/covid-data-tracker/](https://covid.cdc.gov/covid-data-tracker/)>.
    
    

5.  5.Leeb RT, Price S, Sliwa S, Kimball A, Szucs L, Caruso E, Godfred-Cato S, Lozier M. COVID-19 Trends Among School-Aged Children - United States, March 1- September 19, 2020. MMWR Morb Mortal Wkly Rep 2020;69:1410–1415.
    
    

6.  6.Sungnak W, Huang N, Bécavin C, Berg M, Queen R, Litvinukova M, Talavera-López C, Maatz H, Reichart D, Sampaziotis F, Worlock KB, Yoshida M, Barnes JL, HCA Lung Biological Network. SARS-CoV-2 entry factors are highly expressed in nasal epithelial cells together with innate immune genes. Nat Med 2020;26:681–687.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0868-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

7.  7.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 2nd, 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, et al. 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 2020;181:1016–1035.e19.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.04.035&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

8.  8.Hoffmann M, Kleine-Weber H, Schroeder S, Krüger N, Herrler T, Erichsen S, Schiergens TS, Herrler G, Wu N-H, Nitsche A, Müller MA, Drosten C, Pöhlmann S. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020;181:271–280.e8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.02.052&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

9.  9.Walls AC, Park Y-J, Tortorici MA, Wall A, McGuire AT, Veesler D. Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell 2020;183:1735.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.11.032&link_type=DOI) 

10. 10.Wrapp D, Wang N, Corbett KS, Goldsmith JA, Hsieh C-L, Abiona O, Graham BS, McLellan JS. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 2020;367:1260–1263.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjcvNjQ4My8xMjYwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDEvMjgvMjAyMS4wMS4yNi4yMTI1MDI2OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

11. 11.Bunyavanich S, Do A, Vicencio A. Nasal Gene Expression of Angiotensin-Converting Enzyme 2 in Children and Adults. JAMA 2020;323:2427–2429.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.8707&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

12. 12.Muus C, Luecken MD, Eraslan G, Waghray A, Heimberg G, Sikkema L, Kobayashi Y, Vaishnav ED, Subramanian A, Smilie C, Others. Integrated analyses of single-cell atlases reveal age, gender, and smoking status associations with cell type-specific expression of mediators of SARS-CoV-2 viral entry and highlights inflammatory programs in putative target cells. BioRxiv 2020;at <[https://www.biorxiv.org/content/10.1101/2020.04.19.049254v1.full-text](https://www.biorxiv.org/content/10.1101/2020.04.19.049254v1.full-text)>.
    
    

13. 13.Yonker LM, Neilan AM, Bartsch Y, Patel AB, Regan J, Arya P, Gootkind E, Park G, Hardcastle M, St John A, Appleman L, Chiu ML, Fialkowski A,  De la Flor D, Lima R, Bordt EA, Yockey LJ, D’Avino P, Fischinger S, Shui JE, Lerou PH, Bonventre JV, Yu XG, Ryan ET, Bassett IV, Irimia D, Edlow AG, Alter G, Li JZ, et al. Pediatric Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2): Clinical Presentation, Infectivity, and Immune Responses. J Pediatr 2020;227:45–52.e5.
    
    

14. 14.Jones TC, Mühlemann B, Veith T, Biele G, Zuchowski M, Hoffmann J, Stein A, Edelmann A, Corman VM, Drosten C. An analysis of SARS-CoV-2 viral load by patient age. Infectious Diseases (except HIV/AIDS) 2020;doi:10.1101/2020.06.08.20125484.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNi4wOC4yMDEyNTQ4NHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDEvMjgvMjAyMS4wMS4yNi4yMTI1MDI2OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

15. 15.Heald-Sargent T, Muller WJ, Zheng X, Rippe J, Patel AB, Kociolek LK. Age-Related Differences in Nasopharyngeal Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Levels in Patients With Mild to Moderate Coronavirus Disease 2019 (COVID-19). JAMA Pediatr 2020;174:902–903.
    
    

16. 16.Baggio S, L’Huillier AG, Yerly S, Bellon M, Wagner N, Rohr M, Huttner A, Blanchard-Rohner G, Loevy N, Kaiser L, Jacquerioz F, Eckerle I. SARS-CoV-2 viral load in the upper respiratory tract of children and adults with early acute COVID-19. Clin Infect Dis 2020;doi:10.1093/cid/ciaa1157.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa1157&link_type=DOI) 

17. 17.Iwasaki A, Foxman EF, Molony RD. Early local immune defences in the respiratory tract. Nat Rev Immunol 2017;17:7–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri.2016.117&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27890913&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

18. 18.Blanco-Melo D, Nilsson-Payant BE, Liu W-C, Uhl S, Hoagland D, Møller R, Jordan TX, Oishi K, Panis M, Sachs D, Wang TT, Schwartz RE, Lim JK, Albrecht RA, tenOever BR. Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19. Cell 2020;181:1036–1045.e9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.04.026&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

19. 19.Hadjadj J, Yatim N, Barnabei L, Corneau A, Boussier J, Smith N, Péré H, Charbit B, Bondet V, Chenevier-Gobeaux C, Breillat P, Carlier N, Gauzit R, Morbieu C, Pène F, Marin N, Roche N, Szwebel T-A, Merkling SH, Treluyer J-M, Veyer D, Mouthon L, Blanc C, Tharaux P-L, Rozenberg F, Fischer A, Duffy D, Rieux-Laucat F, Kernéis S, et al. Impaired type I interferon activity and inflammatory responses in severe COVID-19 patients. Science 2020;doi:10.1126/science.abc6027.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjkvNjUwNC83MTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMS8yOC8yMDIxLjAxLjI2LjIxMjUwMjY5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

20. 20.Do LAH, Pellet J, van Doorn HR, Tran AT, Nguyen BH, Tran TTL, Tran QH, Vo QB, Tran Dac NA, Trinh HN, Nguyen TTH, Le Binh BT, Nguyen HMK, Nguyen MT, Thai QT, Vo TV, Ngo NQM, Dang TKH, Cao NH, Van Tran T, Ho LV, De Meulder B, Auffray C, Hofstra J-J, Farrar J, Bryant JE, de Jong M, Hibberd ML. Host Transcription Profile in Nasal Epithelium and Whole Blood of Hospitalized Children Under 2 Years of Age With Respiratory Syncytial Virus Infection. J Infect Dis 2017;217:134–146.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jix519&link_type=DOI) 

21. 21.Wesolowska-Andersen A, Everman JL, Davidson R, Rios C, Herrin R, Eng C, Janssen WJ, Liu AH, Oh SS, Kumar R, Fingerlin TE, Rodriguez-Santana J, Burchard EG, Seibold MA. Dual RNA-seq reveals viral infections in asthmatic children without respiratory illness which are associated with changes in the airway transcriptome. Genome Biol 2017;18:12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-016-1140-8&link_type=DOI) 

22. 22.Poole A, Urbanek C, Eng C, Schageman J, Jacobson S, O’Connor BP, Galanter JM, Gignoux CR, Roth LA, Kumar R, Lutz S, Liu AH, Fingerlin TE, Setterquist RA, Burchard EG, Rodriguez-Santana J, Seibold MA. Dissecting childhood asthma with nasal transcriptomics distinguishes subphenotypes of disease. J Allergy Clin Immunol 2014;133:670–8.e12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2013.11.025&link_type=DOI) 

23. 23.Altman MC, Gill MA, Whalen E, Babineau DC, Shao B, Liu AH, Jepson B, Gruchalla RS, O’Connor GT, Pongracic JA, Kercsmar CM, Khurana Hershey GK, Zoratti EM, Johnson CC, Teach SJ, Kattan M, Bacharier LB, Beigelman A, Sigelman SM, Presnell S, Gern JE, Gergen PJ, Wheatley LM, Togias A, Busse WW, Jackson DJ. Transcriptome networks identify mechanisms of viral and nonviral asthma exacerbations in children. Nat Immunol 2019;20:637–651.
    
    

24. 24.Harris PA, Taylor R, Thielke R, Payne J, Gonzalez N, Conde JG, Others. A metadata-driven methodology and workflow process for providing translational research informatics support. J Biomed Inform 2009;42:377–381.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jbi.2008.08.010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18929686&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264958800018&link_type=ISI) 

25. 25.Pediatric Acute Lung Injury Consensus Conference Group. Pediatric acute respiratory distress syndrome: consensus recommendations from the Pediatric Acute Lung Injury Consensus Conference. Pediatr Crit Care Med 2015;16:428–439.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/PCC.0000000000000350&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25647235&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

26. 26.Sala MA, Balderas-Martínez YI, Buendía-Roldan I, Abdala-Valencia H, Nam K, Jain M, Bhorade S, Bharat A, Reyfman PA, Ridge KM, Pardo A, Sznajder JI, Budinger GRS, Misharin AV, Selman M. Inflammatory pathways are upregulated in the nasal epithelium in patients with idiopathic pulmonary fibrosis. Respir Res 2018;19:233.
    
    

27. 27.Di Tommaso P, Chatzou M, Floden EW, Barja PP, Palumbo E, Notredame C. Nextflow enables reproducible computational workflows. Nat Biotechnol 2017;35:316–319.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.3820&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28398311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

28. 28.Ewels PA, Peltzer A, Fillinger S, Patel H, Alneberg J, Wilm A, Garcia MU, Di Tommaso P, Nahnsen S. The nf-core framework for community-curated bioinformatics pipelines. Nat Biotechnol 2020;38:276–278.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41587-020-0439-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32055031&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

29. 29.Kurtzer GM, Sochat V, Bauer MW. Singularity: Scientific containers for mobility of compute. PLoS One 2017;12:e0177459.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0177459&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28494014&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

30. 30.Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 2013;29:15–21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts635&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23104886&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312654600003&link_type=ISI) 

31. 31.Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014;30:923–930.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btt656&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24227677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000334078300005&link_type=ISI) 

32. 32.Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 2008;9:559.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-9-559&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19114008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

33. 33.Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics 2009;10:48.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-10-48&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19192299&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

34. 34.Aliee H, Theis F. AutoGeneS: Automatic gene selection using multi-objective optimization for RNA-seq deconvolution. bioRxiv 2020;at <[https://www.biorxiv.org/content/10.1101/2020.02.21.940650v1.abstract](https://www.biorxiv.org/content/10.1101/2020.02.21.940650v1.abstract)>.
    
    

35. 35.Ordovas-Montanes J, Dwyer DF, Nyquist SK, Buchheit KM, Vukovic M, Deb C, Wadsworth MH 2nd, Hughes TK, Kazer SW, Yoshimoto E, Cahill KN, Bhattacharyya N, Katz HR, Berger B, Laidlaw TM, Boyce JA, Barrett NA, Shalek AK. Allergic inflammatory memory in human respiratory epithelial progenitor cells. Nature 2018;560:649–654.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

36. 36.Sajuthi SP, DeFord P, Li Y, Jackson ND, Montgomery MT, Everman JL, Rios CL, Pruesse E, Nolin JD, Plender EG, Wechsler ME, Mak ACY, Eng C, Salazar S, Medina V, Wohlford EM, Huntsman S, Nickerson DA, Germer S, Zody MC, Abecasis G, Kang HM, Rice KM, Kumar R, Oh S, Rodriguez-Santana J, Burchard EG, Seibold MA. Type 2 and interferon inflammation regulate SARS-CoV-2 entry factor expression in the airway epithelium. Nat Commun 2020;11:5139.
    
    

37. 37.Chu H, Chan JF-W, Wang Y, Yuen TT-T, Chai Y, Hou Y, Shuai H, Yang D, Hu B, Huang X, Zhang X, Cai J-P, Zhou J, Yuan S, Kok K-H, To KK-W, Chan IH-Y, Zhang AJ, Sit K-Y, Au W-K, Yuen K-Y. 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 2020;doi:10.1093/cid/ciaa410.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa410&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32270184&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

38. 38.Onabajo OO, Banday AR, Yan W, Obajemu A, Stanifer ML, Santer DM, Florez-Vargas O, Piontkivska H, Vargas J, Kee C, Tyrrell DLJ, Mendoza JL, Boulant S, Prokunina-Olsson L. Interferons and viruses induce a novel primate-specific isoform dACE2 and not the SARS-CoV-2 receptor ACE2. bioRxiv 2020;doi:10.1101/2020.07.19.210955.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wNy4xOS4yMTA5NTV2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAxLzI4LzIwMjEuMDEuMjYuMjEyNTAyNjkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

39. 39.de Lang A, Osterhaus ADME, Haagmans BL. Interferon-gamma and interleukin-4 downregulate expression of the SARS coronavirus receptor ACE2 in Vero E6 cells. Virology 2006;353:474–481.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.virol.2006.06.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16860835&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

40. 40.Deprez M, Zaragosi L-E, Truchi M, Becavin C, Ruiz García S, Arguel M-J, Plaisant M, Magnone V, Lebrigand K, Abelanet S, Brau F, Paquet A, Pe’er D, Marquette C-H, Leroy S, Barbry P. A Single-cell Atlas of the Human Healthy Airways. Am J Respir Crit Care Med 2020;doi:10.1164/rccm.201911-2199OC.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.201911-2199OC&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32726565&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

41. 41.Vieira Braga FA, Kar G, Berg M, Carpaij OA, Polanski K, Simon LM, Brouwer S, Gomes T, Hesse L, Jiang J, Fasouli ES, Efremova M, Vento-Tormo R, Talavera-López C, Jonker MR, Affleck K, Palit S, Strzelecka PM, Firth HV, Mahbubani KT, Cvejic A, Meyer KB, Saeb-Parsy K, Luinge M, Brandsma C-A, Timens W, Angelidis I, Strunz M, Koppelman GH, et al. A cellular census of human lungs identifies novel cell states in health and in asthma. Nat Med 2019;25:1153–1163.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/S41591-019-0468-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31209336&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

42. 42.Persson BD, Jaffe AB, Fearns R, Danahay H. Respiratory syncytial virus can infect basal cells and alter human airway epithelial differentiation. PLoS One 2014;9:e102368.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0102368&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25033192&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

43. 43.Sebina I, Phipps S. The Contribution of Neutrophils to the Pathogenesis of RSV Bronchiolitis. Viruses 2020;12.:
    
    

44. 44.Steinbusch MMF, Fang Y, Milner PI, Clegg PD, Young DA, Welting TJM, Peffers MJ. Serum snoRNAs as biomarkers for joint ageing and post traumatic osteoarthritis. Sci Rep 2017;7:43558.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/srep43558&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28252005&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

45. 45.Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 2005;102:15545–15550.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTAyLzQzLzE1NTQ1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDEvMjgvMjAyMS4wMS4yNi4yMTI1MDI2OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

46. 46.Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics 2011;27:1739–1740.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr260&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21546393&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291261300036&link_type=ISI) 

47. 47.Steinman JB, Lum FM, Ho PP-K, Kaminski N, Steinman L. Reduced development of COVID-19 in children reveals molecular checkpoints gating pathogenesis illuminating potential therapeutics. Proc Natl Acad Sci U S A 2020;117:24620–24626.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzQwLzI0NjIwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDEvMjgvMjAyMS4wMS4yNi4yMTI1MDI2OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

48. 48.Zimmermann P, Curtis N. Why is COVID-19 less severe in children? A review of the proposed mechanisms underlying the age-related difference in severity of SARS-CoV-2 infections. Arch Dis Child 2020;doi:10.1136/archdischild-2020-320338.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImFyY2hkaXNjaGlsZCI7czo1OiJyZXNpZCI7czoyNjoiYXJjaGRpc2NoaWxkLTIwMjAtMzIwMzM4djEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMS8yOC8yMDIxLjAxLjI2LjIxMjUwMjY5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

49. 49.Mick E, Kamm J, Pisco AO, Ratnasiri K, Babik JM, Castañeda G, DeRisi JL, Detweiler AM, Hao SL, Kangelaris KN, Kumar GR, Li LM, Mann SA, Neff N, Prasad PA, Serpa PH, Shah SJ, Spottiswoode N, Tan M, Calfee CS, Christenson SA, Kistler A, Langelier C. Upper airway gene expression reveals suppressed immune responses to SARS-CoV-2 compared with other respiratory viruses. Nat Commun 2020;11:5854.
    
    

50. 50.Sariol A, Perlman S. Lessons for COVID-19 Immunity from Other Coronavirus Infections. Immunity 2020;53:248–263.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

51. 51.Gioia C, Horejsh D, Agrati C, Martini F, Capobianchi MR, Ippolito G, Poccia F. T-Cell response profiling to biological threat agents including the SARS coronavirus. Int J Immunopathol Pharmacol 2005;18:525–530.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/039463200501800312&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16164833&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000232655100012&link_type=ISI) 

52. 52.Grifoni A, Weiskopf D, Ramirez SI, Mateus J, Dan JM, Moderbacher CR, Rawlings SA, Sutherland A, Premkumar L, Jadi RS, Marrama D, de Silva AM, Frazier A, Carlin AF, Greenbaum JA, Peters B, Krammer F, Smith DM, Crotty S, Sette A. Targets of T Cell Responses to SARS-CoV-2 Coronavirus in Humans with COVID- 19 Disease and Unexposed Individuals. Cell 2020;181:1489–1501.e15.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=WOS:000543822100009.English&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

53. 53.Le Bert N, Tan AT, Kunasegaran K, Tham CYL, Hafezi M, Chia A, Chng MHY, Lin M, Tan N, Linster M, Chia WN, Chen MI-C, Wang L-F, Ooi EE, Kalimuddin S, Tambyah PA, Low JG-H, Tan Y-J, Bertoletti A. SARS-CoV-2-specific T cell immunity in cases of COVID-19 and SARS, and uninfected controls. Nature 2020;doi:10.1038/s41586-020-2550-z.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2550-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32668444&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

54. 54.Mateus J, Grifoni A, Tarke A, Sidney J, Ramirez SI, Dan JM, Burger ZC, Rawlings SA, Smith DM, Phillips E, Mallal S, Lammers M, Rubiro P, Quiambao L, Sutherland A, Yu ED, da Silva Antunes R, Greenbaum J, Frazier A, Markmann AJ, Premkumar L, de Silva A, Peters B, Crotty S, Sette A, Weiskopf D. Selective and cross-reactive SARS-CoV-2 T cell epitopes in unexposed humans. Science 2020;doi:10.1126/science.abd3871.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjExOiIzNzAvNjUxMi84OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAxLzI4LzIwMjEuMDEuMjYuMjEyNTAyNjkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

55. 55.Braun J, Loyal L, Frentsch M, Wendisch D, Georg P, Kurth F, Hippenstiel S, Dingeldey M, Kruse B, Fauchere F, Baysal E, Mangold M, Henze L, Lauster R, Mall MA, Beyer K, Röhmel J, Voigt S, Schmitz J, Miltenyi S, Demuth I, Müller MA, Hocke A, Witzenrath M, Suttorp N, Kern F, Reimer U, Wenschuh H, Drosten C, et al. SARS-CoV-2-reactive T cells in healthy donors and patients with COVID-19. Nature 2020;587:270–274.
    
    

56. 56.Bacher P, Rosati E, Esser D, Martini GR, Saggau C, Schiminsky E, Dargvainiene J, Schöder I, Wieters I, Khodamoradi Y, Eberhardt F, Maria J G, Neb H, Sonntagbauer M, Conrad C, Tran F, Rosenstiel P, Markewitz R, Wandinger K-P, Augustin M, Rybniker J, Kochanek M, Leypoldt F, Cornely OA, Koehler P, Franke A, Scheffold A. Low avidity CD4+ T cell responses to SARS-CoV-2 in unexposed individuals and humans with severe COVID-19. Immunity 2020;.:
    
    

57. 57.Trump S, Lukassen S, Anker MS, Chua RL, Liebig J, Thürmann L, Corman VM, Binder M, Loske J, Klasa C, Krieger T, Hennig BP, Messingschlager M, Pott F, Kazmierski J, Twardziok S, Albrecht JP, Eils J, Hadzibegovic S, Lena A, Heidecker B, Bürgel T, Steinfeldt J, Goffinet C, Kurth F, Witzenrath M, Völker MT, Müller SD, Liebert UG, et al. Hypertension delays viral clearance and exacerbates airway hyperinflammation in patients with COVID-19. Nat Biotechnol 2020;doi:10.1038/s41587-020-00796-1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41587-020-00796-1&link_type=DOI) 

58. 58.Chua RL, Lukassen S, Trump S, Hennig BP, Wendisch D, Pott F, Debnath O, Thürmann L, Kurth F, Völker MT, Kazmierski J, Timmermann B, Twardziok S, Schneider S, Machleidt F, Müller-Redetzky H, Maier M, Krannich A, Schmidt S, Balzer F, Liebig J, Loske J, Suttorp N, Eils J, Ishaque N, Liebert UG, von Kalle C, Hocke A, Witzenrath M, et al. COVID-19 severity correlates with airway epithelium-immune cell interactions identified by single-cell analysis. Nat Biotechnol 2020;doi:10.1038/s41587-020-0602-4.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41587-020-0602-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32591762&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F28%2F2021.01.26.21250269.atom) 

59. 59.Cao Y, Guo Z, Vangala P, Donnard E, Liu P, McDonel P, Ordovas-Montanes J, Shalek AK, Finberg RW, Wang JP, Garber M. Single-cell analysis of upper airway cells reveals host-viral dynamics in influenza infected adults. Cold Spring Harbor Laboratory 2020;doi:10.1101/2020.04.15.042978.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wNC4xNS4wNDI5Nzh2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAxLzI4LzIwMjEuMDEuMjYuMjEyNTAyNjkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

60. 60.Chan Zuckerberg Initiative Single-Cell COVID-19 Consortia, Ballestar E, Farber DL, Glover S, Horwitz B, Meyer K, Nikolic M, Ordovas-Montanes J, Sims P, Shalek A, Vandamme N, Vandekerckhove L, Vento-Tormo R, Villani AC. Single cell profiling of COVID-19 patients: an international data resource from multiple tissues. bioRxiv 2020;doi:10.1101/2020.11.20.20227355.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2020.11.20.20227355&link_type=DOI)