key: cord-0761128-pcekcvwr authors: Blanco-Melo, Daniel; Nilsson-Payant, Benjamin E.; Liu, Wen-Chun; Uhl, Skyler; Hoagland, Daisy; Møller, Rasmus; Jordan, Tristan X.; Oishi, Kohei; Panis, Maryline; Sachs, David; Wang, Taia T.; Schwartz, Robert E.; Lim, Jean K.; Albrecht, Randy A.; tenOever, Benjamin R. title: Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19 date: 2020-05-15 journal: Cell DOI: 10.1016/j.cell.2020.04.026 sha: dea4b05ab91f7eb501a9205c53c2f510bdccf147 doc_id: 761128 cord_uid: pcekcvwr Viral pandemics, such as the one caused by SARS-CoV-2, pose an imminent threat to humanity. Because of its recent emergence, there is a paucity of information regarding viral behavior and host response following SARS-CoV-2 infection. Here we offer an in-depth analysis of the transcriptional response to SARS-CoV-2 compared with other respiratory viruses. Cell and animal models of SARS-CoV-2 infection, in addition to transcriptional and serum profiling of COVID-19 patients, consistently revealed a unique and inappropriate inflammatory response. This response is defined by low levels of type I and III interferons juxtaposed to elevated chemokines and high expression of IL-6. We propose that reduced innate antiviral defenses coupled with exuberant inflammatory cytokine production are the defining and driving features of COVID-19. Correspondence res2025@med.cornell.edu (R.E.S.), jean.lim@mssm.edu (J.K.L.), randy.albrecht@mssm.edu (R.A.A.), benjamin.tenoever@mssm.edu (B.R.t.) In Brief In comparison to other respiratory viruses, SARS-CoV-2 infection drives a lower antiviral transcriptional response that is marked by low IFN-I and IFN-III levels and elevated chemokine expression, which could explain the proinflammatory disease state associated with COVID-19. Coronaviruses are a diverse group of single-stranded positivesense RNA viruses with a wide range of vertebrate hosts (Cui et al., 2019) . Four common coronavirus genera (alpha, beta, gamma, and delta) circulate among vertebrates and cause mild upper respiratory tract illnesses in humans and gastroenteritis in animals (Weiss and Navas-Martin, 2005) . However, in the past two decades, three highly pathogenic human betacoronaviruses have emerged from zoonotic events (Amanat and Krammer, 2020) . In 2002 In -2003 , severe acute respiratory syndromerelated coronavirus 1 (SARS-CoV-1) infected $8,000 people worldwide with a case fatality rate of $10%, followed by Middle East respiratory syndrome-related coronavirus (MERS-CoV), which has infected $2,500 people with a case fatality rate of $36% since 2012 (de Wit et al., 2016) . At present, the world is suffering from a pandemic of SARS-CoV-2, which causes coronavirus disease 2019 and has a global mortality rate that remains to be determined Zhu et al., 2020) . SARS-CoV-2 infection is characterized by a range of symptoms including fever, cough, and general malaise in the majority of cases . More severe cases of COVID-19 develop acute respiratory distress syndrome and acute lung injury, leading to morbidity and mortality caused by damage to the alveolar lumen leading to inflammation and pneumonia (Wolfel et al., 2020; Xu et al., 2020) . The physiological response to virus infection is generally initiated at the cellular level following replication (tenOever, 2016) . After virus entry, the infected cell detects the presence of virus replication through any one of a number of pattern recognition receptors (PRRs) (Janeway and Medzhitov, 2002) . These receptors serve as sentinels for a variety of microbes inside and outside of the cell by physically engaging distinct structures that are shared among different pathogens. In the case of virus infection, cellular detection of replication is largely mediated by a family of intracellular PRRs that sense aberrant RNA structures that often form during virus replication (Janeway and Medzhitov, 2002) . Engagement of virus-specific RNA structures culminates in oligomerization of these receptors and activation of downstream transcription factors, most notably interferon regulator factors (IRFs) and nuclear factor kB (NF-kB) (Hur, 2019) . Transcriptional activation of IRFs and NF-kB results in the launch of two general antiviral programs. The first is engagement of cellular antiviral defenses, which is mediated by transcriptional induction of type I and III interferons (IFN-I and IFN-III, respectively) and subsequent upregulation of IFN-stimulated genes (ISGs) (Lazear et al., 2019) . The second arm of the antiviral response involves recruitment and coordination of specific subsets of leukocytes, which is orchestrated primarily by chemokine secretion (Proudfoot, 2002; Sokol and Luster, 2015) . This broad antiviral response puts selective pressure on viruses and has resulted in the evolution of countless viral countermeasures (García-Sastre, 2017) . Therefore, the host response to a virus is generally not uniform, and infections can inflict different degrees of morbidity and mortality. The current pandemic of COVID-19 is an acute and rapidly developing global health crisis. To better understand the molecular basis of the disease, we sought to characterize the transcriptional response to infection in a variety of model systems, including in vitro tissue culture, ex vivo infection of primary cells, and in vivo samples derived from COVID-19 patients and animals. We chose to characterize the transcriptional response to SARS-CoV-2 and determine how it compares with common respiratory viruses, including influenza A virus (IAV). These two respiratory viruses encode a variety of different antagonists to the IFN-I and -III response (Frieman and Baric, 2008; García-Sastre, 2017) . For the closely related SARS-CoV-1, IFN antagonism has been attributed to ORF3B, ORF6, and the nucleocapsid (N) gene products (Frieman et al., 2010; Kopecky-Bromberg et al., 2007) . SARS-CoV-1 also encodes nsp1, a nuclease that has been implicated in cleaving host mRNA to prevent ribosomal loading and causing host shutoff (Kamitani et al., 2006) . Similar to SARS-CoV-1, IAV also encodes the IFN-I and -III antagonist nonstructural protein 1 (NS1), which blocks initial detection by the PRR through binding and masking aberrant RNA produced during infection . Here we compare the transcriptional response of SARS-CoV-2 with other respiratory viruses to identify transcriptional signatures that may underlie COVID-19 biology. These data demonstrate that the overall transcriptional induction to SARS-CoV-2 is aberrant. Despite virus replication, the host response to SARS-CoV-2 fails to launch a robust IFN-I and -III response while simultaneously inducing high levels of chemokines needed to recruit effector cells. Because a waning immune response would enable sustained viral replication, these findings may explain why serious cases of COVID-19 are more frequently observed in individuals with comorbidities. Defining the Transcriptional Response to SARS-CoV-2 Relative to Other Respiratory Viruses To compare the transcriptional response of SARS-CoV-2 with other respiratory viruses, including MERS-CoV, SARS-CoV-1, human parainfluenza virus 3 (HPIV3), respiratory syncytial virus (RSV), and IAV, we first chose to focus on infection in a variety of respiratory cell lines (Figure 1 ). To this end, we collected poly(A) RNA from infected cells and performed RNA sequencing (RNAseq) to estimate viral load. These data show that virus infection levels ranged from 0.1% to more than 50% of total RNA reads ( Figure 1A ). In agreement with others (Harcourt et al., 2020) , we found A549 lung alveolar cells to be relatively non-permissive to SARS-CoV-2 replication, in contrast to Calu-3 cells (0.1% versus 15% total reads, respectively). The low rate of infection in A549 cells is postulated to be the result of low expression of the viral receptor ACE2 (Harcourt et al., 2020; Hoffmann et al., 2020) . To bypass this restriction, we supplemented A549 cells with a vector expressing mCherry or ACE2 ( Figures 1B-1D) . In low-MOI infections (MOI, 0.2), exogenous expression of ACE2 enabled SARS-CoV-2 to replicate and comprise $54% of the total reads mapping more than 3003 coverage across the $30-kb genome ( Figures 1A and 1B) . Western blot analyses corroborated these RNA-seq data, showing Nucleocapsid (N) expression only in cells supplemented with ACE2 ( Figure 1C ). Furthermore, qPCR analyses of these cells demonstrated that the levels of Envelope (E) and non-structural protein 14 (nsp14) were more than three orders of magnitude higher in the presence of ACE2 ( Figure 1D ). It is noteworthy that, despite this dramatic increase in viral load, we observed neither activation of TBK1, the kinase responsible for IFN-I and IFN-III expression, nor induction of STAT1 and MX1, IFN-I-stimulated genes ( Figure S1A ; Sharma et al., 2003) . The lack of IFN-I and -III engagement in ACE2-expressing A549 cells could, however, be overcome by using a 10-fold increase in virus (MOI, 2) despite the fact that total viral reads after 24 h of replication where comparable with low-MOI conditions ( Figures 1A, 1B , and 1F). To determine whether SARS-CoV-2 is sensitive to IFN-I, we next treated cells with universal IFNb and assessed viral levels at the RNA and protein levels (Figures S1B and S1C). These data demonstrate that addition of IFN-I resulted in a dramatic reduction in virus replication, in agreement with the findings of others (Lokugamage et al., 2020) . We also observed no increase in viral Spike levels (or a considerable effect on viral reads) when IFN-I signaling was blocked by addition of ruxolitinib, a JAK1 and 2 kinase inhibitor, despite significantly preventing induction of ISGs (Figures S1D-S1H). In contrast, ruxolitinib treatment had a minimal effect on induction of cytokines and chemokines, indicating that the high induction of these genes in SARS-CoV-2 infection is independent of IFN-I and -III signaling ( Figure S1G ). To next determine how each of these in vitro infections alters the host transcriptional landscape, we first performed a differential expression analysis, comparing infected cell conditions with their respective mock conditions. These analyses indicate that the transcriptional response in cells that allows high replication of SARS-CoV-2 is significantly different from the host response of all other viruses tested ( Figure 1E ). Moreover, SARS-CoV-2 infection in unmodified A549 cells shows a unique response compared with SARS-CoV-1 despite comparable levels of viral load ( Figures 1A and 1E) . Last, MERS-CoV infection, which approaches $50% of total reads at 24 hours post-infection (hpi), clusters together with SARS-CoV-1 and IAV, reflecting an overall repression of the host antiviral response ( Figures 1E and 1F ). Conversely, HPIV3 and RSV comprise a unique cluster denoted by high expression of IFNs and ISGs ( Figures 1E and 1F ). Interestingly, low-MOI SARS-CoV-2-infected A549 cells expressing ACE2 (A549-ACE2) show no significant IFN-I or IFN-III expression but instead display moderate levels of a subset of ISGs and a unique proinflammatory cytokine signature ( Figure 1F ). This signature is also present in high-MOI infections of SARS-CoV-2 in A549-ACE2 and Calu-3 cells, together with more than 6,000 other differentially expressed genes, further explaining Table S1 ). Furthermore, high-MOI infection in these cells also led to high induction of IFNs and ISGs observed for HPIV3 and RSV, despite remarkable differences in viral replication ($60% total reads in A549-ACE2 cells compared with $15% in Calu-3 cells) (Figures 1A and 1F ). The discrepancy between the levels of viral replication and IFN production/signaling suggests that, although SARS-CoV-2 is capable of engaging the IFN-I and IFN-III systems, this response is prevented by an antagonist that is rendered ineffective under high-MOI conditions. Alternatively, these data may instead indicate that high-MOI conditions in cell culture result in formation of pathogen-associated molecular patterns (PAMPs), which may or may not reflect physiological conditions in vivo. Given the disparate results of our in vitro cell culture systems, we next sought to determine how normal human bronchial epithelial (NHBE) cells respond to SARS-CoV-2 infection, in contrast to treatment with IFN-I alone or infection with wild-type (WT) IAV or a mutant IAV lacking its antiviral antagonist (IAVDNS1) (Figure 2) . Treatment of NHBE cells with IFN-I resulted in significant induction of 381 genes, most of these also differentially expressed in IAVDNS1 infection, together outlining a robust innate immune response in these cells (Figure 2A ). In contrast, and despite different levels of replication, the transcriptional response to infection with SARS-CoV-2 and WT IAV are similar in magnitude but different in nature, with only 8 shared significantly induced genes, including interleukin-6 (IL-6), IRF9, ICAM1, and tumor necrosis factor (TNF) (Figure 2A ; Figure S2A ). To further understand the global host response as it pertained to each of these conditions, we grouped these samples in a PCA space ( Figure 2B ). This analysis shows progressive transcriptional perturbations along principal component one, which accounts for more than 60% of sample variation. In this space, SARS-CoV-2 elicits the most modest transcriptional changes, followed by IAV, IFNb treatment, and IAVDNS1 ( Figure 2B ). Gene enrichment analyses of differentially expressed transcripts illustrate a diminished IFN-I signaling biology for SARS-CoV-2 and IAV infections ( Figures 2C and 2D ). In both examples, IFN-I and IFN-III are undetectable, but a very small subset of ISGs is induced ( Figures 2C, 2D , and S2B; Table S2 ). In the case of IAV, this diminished antiviral response is mediated by expression of NS1 because IAVDNS1 infections result in robust IFNB and IFNL1-3 induction ( Figures 2C, 2D , and S2B; Table S2 ). Despite a complete lack of IFN expression, the response to SARS-CoV-2 in NHBE cells still elicited a strong chemotactic and inflammatory response, indicated by expression of CCL20, CXCL1, IL-1B, IL-6, CXCL3, CXCL5, CXCL6, CXCL2, CXCL16, and TNF ( Figures 2C and 2E ; Table S2 ). In addition to the modest IFN-I response, SARS-CoV-2 in NHBE cells also triggers some unique pathways, including a response to IFN-II (which is also observed in response to IAVDNS1), and significant enrichment in chemokine signaling ( Figure 2C ). In Vitro Response to SARS-CoV-2 To determine whether the limited response to SARS-CoV-2 observed so far was a by-product of cell culture, we next pursued an in vivo longitudinal study in animals. To this end, we chose to perform SARS-CoV-2 infections in ferrets because this has been described as an appropriate animal model (Kim et al., 2020) . Ferrets were infected intranasally with SARS-CoV-2 or influenza A/California/04/2009 and monitored by nasal wash, which generates a small pellet of cells from the upper respiratory tract. RNA-seq was performed on these cells, enabling us to quantify viral load over time. Reads from nasal wash 1 day after infection revealed a low level of virus replication comprising 0.006% of total reads ( Figure 3A ). Three days after infection, virus replication peaked at 1.2% of the total sequencing reads before decreasing to 0.05% of total reads on day 7 and completely clearing the virus by day 14 ( Figure 3A ). As a comparison, a sublethal infection of IAV comprised less than 0.03% of total reads on day 7 from the same sample type ( Figure S3A ). The presence of virus in the nasal passage would further suggest that these ferrets had the potential to transmit the virus, in agreement with the findings of others (Kim et al., 2020; Varble et al., 2014) . To characterize the response to SARS-CoV-2 over time, upper respiratory cell populations were compared with mock-treated ferrets. On day 1 after infection, we observed very little transcriptional difference correlating to the amount of virus detection at this time ( Figure 3B ). By day 3, we observed the beginning of a cytokine response marked by CCL8 and CXCL9, consistent with what was observed in cell culture. By day 7, despite waning levels of virus, the cytokine response continued to expand and included CCL2, CCL8, and CXCL9, among others ( Figure 3B ; Table S3 ). Moreover, we noted evidence of mixed leukocyte infiltration with significant upregulation in CD163, CD226, CCR5, CCR6, CXCR1, CXCR2, and CXCR7 ( Figure 3C ). Overall, the magnitude of this transcriptional response in the upper respiratory tract was significantly lower compared with a comparable IAV infection ( Figure S3B ). However, although IAV induces a greater number of genes, SARS-CoV-2 generates a unique gene signature enriched for cell death and leukocyte activation, including transcripts such as IL1A and CXCL8 (GO: 0008219 and GO: 0045431; Table S3 ). In contrast, the transcriptional footprint of IAV as it pertains to the cellular antiviral response was strikingly greater in magnitude than that observed for SARS-CoV-2 and included the IFN signature genes MX1, ISG20, OASL, and Tetherin ( Figure S3B ; Table S3 ). By day 14, we detected no viral reads for SARS-CoV-2, and the observed cytokines returned to baseline, with the exception of IL-6 and IL1RN or IL1RA, which remained elevated, similar to results observed with MERS (Pascal et al., 2015; Figures 3B and 3C) . (E) Heatmap as in (D) for genes belonging to GO annotations for cytokine activity and chemokine activity (GO: 0005125, GO: 0008009). The graphs depict the log 2 (fold change) of DEGs of infected compared with mock-treated cells. Genes included have a log 2 (fold change) of more than 1 and a p-adjusted value of less than 0.05. Last, to investigate how the host response to SARS-CoV-2 and IAV affected the respiratory tract, we next performed parallel infections and examined the trachea on day 3. With both infections, we observed very low levels of virus but a robust transcriptional response (Table S3) . Gene enrichment analysis of differentially expressed transcripts implicated two populations of immune cell signatures ( Figure 3D ). The first population included common markers for monocytes and lymphocytes, and the induction of these genes was comparable between SARS-CoV-2 and IAV ( Figures 3D and S3C) . Intriguingly, unique gene signatures from SARS-CoV-2-infected trachea that were largely absent in response to IAV align with those of progenitor cells from the hematopoietic lineage, suggesting that infection may be inducing hematopoiesis ( Figures 3D and S3C ) (Lefranç ais et al., 2017; Yoshida et al., 2019) . Additional research in this area will be required to ascertain whether this is a contributing factor toward the development of COVID-19. Following the characterization of SARS-CoV-2 infection in ferrets, we next sought to correlate these results with natural human infections. To this end, we first compared post-mortem lung samples from COVID-19-positive patients with biopsied healthy lung tissue from uninfected individuals. Transcriptional profiling of these samples, all derived from males older than 60 years (n = 2 for each group), demonstrated $2,000 differen-tially expressed genes with enrichment for the innate and humoral responses (Figures 4A and S4A ). Genes significantly induced in response to SARS-CoV-2 included a subset of ISGs with no IFN-I or IFN-III detected by RNA-seq or semiquantitative PCR ( Figure S4B ; Table S4 ). In addition to genes implicated in innate antiviral immunity, SARS-CoV-2 also induced robust levels of chemokines, including CCL2, CCL8, and CCL11 ( Figure 4A ). Despite the limited number of patients analyzed, these data corroborate our findings in NHBE and ferrets ( Figure S4A ). Next we wished to further validate our findings with a larger cohort of patients through direct detection of circulating cytokines induced by SARS-CoV-2 infection. To this end, we obtained serum from two cohorts of individuals from the Kaiser Santa Clara testing facility (Santa Clara, CA). These two cohorts tested positive for SARS-CoV-2 by nasopharyngeal swabs or were admitted to the hospital for non-COVID-19-related respiratory issues (n = 24 for each group). During initial analyses, these serum samples consistently tested negative for IFNb and the IFNl family of IFNs ( Figure 4B ). Moreover, analyses of cytokines and chemokines quantified in individual serum samples revealed enhancement of generalized inflammation among COVID-19 patients, marked by a significant increase in circulating IL-6, IL1RA, CCL2, CCL8 CXCL2, CXCL8, CXCL9, and CXCL16 levels (Figure 4C ). Significant elevation of CXCL9 and CXCL16 (chemoattractants of T or natural killer (NK) cells, respectively), CCL8 and CCL2 (which recruit monocytes and/or macrophages), and CXCL8 (a classic neutrophil chemoattractant) suggest that the presence of these cells may be a primary driver of the signature pathology observed in COVID-19 patients (Proudfoot, 2002) . Although this sample size is not necessarily representative of the whole population of infected COVID-19 patients, our data are consistent with what we observe using our other model systems. Additional sampling will be required to validate these findings. In the present study, we focus on defining the host response to SARS-CoV-2 and other human respiratory viruses in cell lines, primary cell cultures, ferrets, and COVID-19 patients. In general, our data show that the overall transcriptional footprint of SARS-CoV-2 infection was distinct in comparison with other highly pathogenic coronaviruses and common respiratory viruses such as IAV, HPIV3, and RSV. It is noteworthy that, despite a reduced IFN-I and -III response to SARS-CoV-2, we observed a consistent chemokine signature. One exception to this observation is the response to high-MOI infection in A549-ACE2 and Calu-3 cells, where replication was robust and an IFN-I and -III signature could be observed. In both of these examples, cells were infected at a rate to theoretically deliver two functional virions per cell in addition to any defective interfering particles within the virus stock that were not accounted for by plaque assays. Under these conditions, the threshold for PAMP may be achieved prior to the ability of the virus to evade detection through production of a viral antagonist. Alternatively, addition of multiple genomes to a single cell may disrupt the stoichiometry of viral components, which, in turn, may itself generate A C B (A) Volcano plot depicting DEGs in post-mortem lung samples of two COVID-19 patients compared with healthy lung biopsies. DEGs (p-adjusted < 0.05) with a |log 2 (fold change)| of more than 2 are indicated in red. Non-significant DEGs with a |log 2 (fold change)| of more than 2 are indicated in green. (B and C) Cytokine profiles of COVID-19 patients. Sera of 24 COVID-19 patients and 24 SARS-CoV-2-negative controls were analyzed by ELISA for the protein levels of (B) IFN-I and IFN-III or (C) a broad panel of cytokines. The dotted line depicts the limit of detection. Statistical significance was calculated by Mann-Whitney non-parametric t test. NS, non-significant; *p < 0.05, **p < 0.005, ***p < 0.0001. PAMPs that would not form otherwise. These ideas are supported by the fact that, at a low-MOI infection in A549-ACE2 cells, high levels of replication could also be achieved, but in the absence of IFN-I and -III induction. Taken together, these data suggest that, at low MOIs, the virus is not a strong inducer of the IFN-I and -III system, as opposed to conditions where the MOI is high. These dynamics are also likely to contribute to development of COVID-19 during the course of infection (Wolfel et al., 2020) . A recurrent observation in each of our systems is a robust production of cytokines and their subsequent transcriptional response. According to our longitudinal in vivo data, this response starts as early as 3 days after infection and continues beyond clearance of the virus. A recent study analyzing severe versus mild cases of COVID-19 showed that peripherally derived macrophages predominated in the lungs of severe cases (Liao et al., 2020) . Consistent with this, we found, in all of our systems, significant induction of monocyte-associated chemokines such as CCL2 and CCL8. In addition, our data suggest that neutrophils could also contribute to the disease observed in COVID-19 patients, as demonstrated by CXCL2 and CXCL8 induction. This is consistent with data showing elevated circulating neutrophil levels among COVID-19 patients Qin et al., 2020) , which may have prognostic value for identifying individuals at risk for developing severe disease. It is also noteworthy that two of the cytokines uniquely elevated in response to SARS-CoV-2 are IL-6 and IL1RA, suggesting that there might be a parallel between COVID-19 and cytokine release syndrome (CRS), a complication commonly seen following CAR T-cell treatment (Giavridis et al., 2018) . Should this be true, drugs such as tocilizumab and anakinra may prove to be beneficial for treatment of COVID-19 (Norelli et al., 2018) . Future studies will be needed to address this formally. Like SARS-CoV-2, the clinical manifestation of SARS-CoV-1 has been proposed to stem from a dysregulated immune response in patients and delayed expression of IFN-I (Channappanavar et al., 2016; Law et al., 2005; Menachery et al., 2014) . Based on animal models, SARS-CoV-1 has been found to induce a robust cytokine response that generally shows a delay in IFN-I, culminating in improper recruitment of inflammatory monocytemacrophage populations (Channappanavar et al., 2016) . This dynamic seems to be in line with what we observed with SARS-CoV-2 because low levels of IFN-I and -III are likely produced in response to infection. Given the moderate viral replication levels observed in vivo, one explanation for the low IFN expression could be that a small subset of cells is refractory to the antagonistic mechanism of SARS-CoV-2, producing sufficient amounts of IFN-I and/or IFN-III to guide immune cell activation and ISG induction. What makes SARS-CoV-2 distinct from other viruses used in this study is the propensity to selectively induce morbidity and mortality in older populations (Novel Coronavirus Pneumonia Emergency Response Epidemiology Team, 2020). The physiological basis for this morbidity is believed to be selective death of type II pneumocytes, which results in loss of air exchange and fluid leakage into the lungs (Qian et al., 2013; Xu et al., 2020) . Although it remains to be determined whether the inappropriate inflammatory response to SARS-CoV-2 is responsible for the abnormally high lethality in older populations, it does explain why the virus is generally asymptomatic in young people with healthy and robust immune systems . Given the results here, it is tempting to speculate that an already restricted immune response in the aging population prevents successful inhibition of viral spread at early stages of infection, further exacerbating the morbidity and mortality observed for this age group (Jing et al., 2009; Montecino-Rodriguez et al., 2013) . Taken together, the data presented here suggest that the response to SARS-CoV-2 is imbalanced with regard to controlling virus replication versus activation of the adaptive immune response. Given this dynamic, treatments for COVID-19 have less to do with the IFN response and more to do with controlling inflammation. Because our data suggest that numerous chemokines and ILs are elevated in COVID-19 patients, future efforts should focus on U.S. Food and Drug Administration (FDA)approved drugs that can be rapidly deployed and have immunomodulating properties. Detailed methods are provided in the online version of this paper and include the following: PSMB9 IFITM3 HERC6 BST2 OAS3 PARP9 GBP4 LAMP3 NLRC5 GBP2 TRIM25 MYD88 TAP1 TRIM14 TRIM5 EIF2AK2 OGFR TRIM21 SLC25A28 TDRD7 IRF9 ISG20 SP100 PNPT1 C19orf66 PLSCR1 PARP12 IFITM2 STAT1 TLR3 MOV10 TMEM140 NMI CNP PSMB8 TREX1 HLA−C NUB1 ADAR TXNIP TRIM56 HLA−H RNF31 HLA−E FADD JAK1 CNOT7 LSM14A PDE12 IFNAR2 AXL TRAFD1 PTPN1 TRIM26 HLA−A PTPN2 IRF5 ABCE1 TPR CACTIN ELF1 IL15 CASP8 B2M NCOA7 GATA3 IL4R IRF2 HLA−B TYK2 SETD2 KLHL20 RIPK2 CSF1 IP6K2 TBK1 LPAR6 IRF4 STAT2 USP18 SP110 DHX58 EPSTI1 DDX60 GMPR IRF1 HLA−F IFITM1 BATF2 IFIH1 IRF7 IFI44 CMPK2 OASL IFNB1 CoV-2 IFN IAV IAV NS1 IFNA1 IFNA10 IFNA13 IFNA14 IFNA16 IFNA17 IFNA2 IFNA21 IFNA22P IFNA4 IFNA5 IFNA6 IFNA7 IFNA8 IFNAR1 IFNAR2 IFNB1 IFNE IFNG IFNGR1 IFNGR2 IFNK IFNL1 IFNL2 IFNL3 IFNL4 IFNLR1 IFNW1 0 Log2 fold change -Log10P Log2 fold change MX1 OASL Tetherin ISG20 Figure S3 . Transcriptional Response to SARS-CoV-2 and IAV in Ferrets, Related to Figure 3 (A) Read coverage along the IAV genome. Graph indicates the number of viral reads per each position of the IAV virus genome identified in RNA extracted from nasal washes of ferrets at 7 days post infection. Scaled model of the concatenated IAV segments is depicted below. (B) Volcano plots indicating differentially expressed genes of ferrets infected with SARS-CoV-2 or IAV for 7 days. Differentially expressed genes (p-adjusted value < 0.05) with a |Log 2 (Fold Change)| > 2 are indicated in red. Non-significant differentially expressed Genes with a |Log 2 (Fold Change)| > 2 are indicated in green. (C) Cellular profiling from a subset of genes selectively enriched in response to SARS-CoV-2 compared to IAV, as determined by the Immunological Genome Project. Figure S4 . Unique and Shared Biological Processes between Different Models of SARS-CoV-2 Infection, Related to Figure 4 (A) Dotplot visualization of enriched GO terms in NHBE cells, ferrets and COVID-19 patients. Gene enrichment analyses were performed using STRING against the GO dataset for biological processes. The color of the dots represents the false discovery rate (FDR) value for each enriched GO term and its size represents the percentage of genes enriched in the total gene set. (B) Semiquantitative PCR analysis of healthy and COVID-19 derived lung tissues. Image indicates the relative expression of SARS-CoV-2 nsp14, IFNB and tubulin transcripts in healthy human biopsies and biological replicates of lung tissue from COVID-19 patients. Additionally cDNA of A549 cells infected with IAVDNS1 are included as controls. ll Article SARS-CoV-2 Vaccines: Status Report An inability to maintain the ribonucleoprotein genomic structure is responsible for host detection of negative-sense RNA viruses Dysregulated Type I Interferon and Inflammatory Monocyte-Macrophage Responses Cause Lethal Pneumonia in SARS-CoV-Infected Mice Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study Molecular Diagnosis of a Novel Coronavirus (2019-nCoV) Causing an Outbreak of Pneumonia Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR Origin and evolution of pathogenic coronaviruses SARS and MERS: recent insights into emerging coronaviruses Mechanisms of severe acute respiratory syndrome pathogenesis and innate immunomodulation. Microbiol SARS-CoV pathogenesis is regulated by a STAT1 dependent but a type I, II and III interferon receptor independent mechanism Ten Strategies of Interferon Evasion by Viruses Influenza A virus lacking the NS1 gene replicates in interferon-deficient systems CAR T cell-induced cytokine release syndrome is mediated by macrophages and abated by IL-1 blockade Iduronic acid-containing glycosaminoglycans on target cells are required for efficient respiratory syncytial virus infection Severe Acute Respiratory Syndrome Coronavirus 2 from Patient with 2019 Novel Coronavirus Disease, United States SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor Double-Stranded RNA Sensors and Modulators in Innate Immunity Shared and Distinct Functions of Type I and Type III Interferons The lung is a site of platelet biogenesis and a reservoir for haematopoietic progenitors The landscape of lung bronchoalveolar immune cells in COVID-19 revealed by single-cell RNA sequencing Sequential Immunization With Live-Attenuated Chimeric Hemagglutinin-Based Vaccines Confers Heterosubtypic Immunity Against Influenza A Viruses in a SARS-CoV-2 is sensitive to type I interferon pretreatment Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 SARS-CoV-2 Infection in Children Pathogenic influenza viruses and coronaviruses utilize similar and contrasting approaches to control interferon-stimulated gene responses Causes, consequences, and reversal of immune system aging Monocytederived IL-1 and IL-6 are differentially required for cytokine-release syndrome and neurotoxicity due to CAR T cells The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) in China MERS-CoV infection. Proc. Natl. Acad. Sci. USA Chemokine receptors: multifaceted therapeutic targets Innate immune response of human alveolar type II cells Dysregulation of immune response in patients with COVID-19 in Wuhan, China Triggering the interferon antiviral response through an IKK-related pathway The chemokine system in innate immunity STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets The Evolution of Antiviral Defense Systems Influenza A virus transmission bottlenecks are defined by infection route and recipient host Coronavirus pathogenesis and the emerging pathogen severe acute respiratory syndrome coronavirus. Microbiol ggplot2: Elegant Graphics for Data Analysis A penalized matrix decomposition, with applications to sparse principal components and canonical correlation analysis Virological assessment of hospitalized patients with COVID-2019 A new coronavirus associated with human respiratory disease in China Pathological findings of COVID-19 associated with acute respiratory distress syndrome The cis-Regulatory Atlas of the Mouse Immune System A Novel Coronavirus from Patients with Pneumonia in China Materials Availability All materials and reagents will be made available upon instalment of a material transfer agreement (MTA) The original sequencing datasets for SARS-CoV-1 and MERS-CoV infections can be found on the NCBI Gene Expression Omnibus (GEO) server under the accession number GSE56192. EXPERIMENTAL MODEL AND SUBJECT AVAILABILITY Cell cultures and primary cells Human adenocarcinomic alveolar basal epithelial (A549) cells (ATCC, CCL-185), human adenocarcinomic lung epithelial (Calu-3) cells (ATCC, HTB-55), human HEp-2 cells (ATCC, CCL-23), Madin-Darby Canine Kidney (MDCK) cells (ATCC, CCL-34), MDCK-NS1 cells (García-Sastre et al., 1998) and African green monkey kidney epithelial Vero E6 cells (ATCC, CRL-1586) were maintained at 37 C and 5% CO 2 in Dulbecco's Modified Eagle Medium (DMEM, GIBCO) supplemented with 10% Fetal Bovine Serum (FBS, Corning) All animals were confirmed to be seronegative for circulating influenza A (H1N1) viruses, influenza A (H3N2) viruses and influenza B viruses prior to purchase. Ferrets were housed in cages in the enhanced BSL-3 facility of the Emerging Pathogens Institute at the Icahn School of Medicine at Mount Sinai. All animal experiments were performed according to protocols approved by the Institutional Animal Care and Use Committee (IACUC) and Institutional Biosafety Committee of the Icahn School of For semiquantitative PCR analyses, two additional lung samples were derived post-mortem from males over 60 years of age. The uninfected human lung samples (n = 2) were obtained post surgery through the Mount Sinai Institutional Biorepository and Molecular Pathology Shared Resource Facility (SRF) in the Department of Pathology. The Biorepository operates under a Mount Sinai Institutional Review Board (IRB) approved protocol and follows guidelines set by HIPAA. Sera were obtained from the Kaiser Santa Clara testing facility USA) and titers in Vero E6 cells as previously described. All work involving live SARS-CoV-2 was performed in the CDC/USDA-approved BSL-3 facility of the Icahn School of Medicine at Mount Sinai in accordance with institutional biosafety requirements. METHOD DETAILS RNA-Seq of viral infections Approximately 5 3 10 5 A549 or Calu-3 cells were infected with influenza A/Puerto Rico/8/1934 (H1N1) virus (IAV), human respiratory syncytial virus (RSV), human parainfluenza virus 3 (HPIV3) or SARS-CoV-2 as indicated. Infections with IAV were performed at a multiplicity of infection of 5 for 9 h in DMEM supplemented with 0.3% BSA, 4.5 g/L D-glucose, 4 mM L-glutamine and 1 mg/ml TPCKtrypsin. Infections with RSV and HPIV3 were performed at an MOI of 2 for 24 h in DMEM supplemented with 2% FBS, 4.5 g/L D-glucose and 4 mM L-glutamine. Infections with SARS-CoV-2 were performed at an MOI of 0.2 or 2 for 24 h in DMEM supplemented with 2% FBS, 4.5 g/L D-glucose, 4 mM L-glutamine, 10 mM Non-Essential Amino Acids, 1 mM Sodium Pyruvate and 10 mM HEPES. Approximately 5 3 10 5 NHBE cells were infected with either SARS-CoV-2 at an MOI of 2 for 24 h or influenza A/Puerto Rico/8/1934 (H1N1) virus or influenza A/Puerto Rico/8/1934 (H1N1) virus lacking the NS1 gene at an MOI of 3 for 12 h in bronchial epithelial growth media supplemented with BEGM SingleQuots. As a comparison to viral infection, NHBE cells were treated with 100 units / mL of IFNb for 4 -12 h. Total RNA from infected and mock infected cells was lysed in TRIzol (Invitrogen) and extracted and DNase I treated using Direct-zol RNA Miniprep kit (Zymo Research) according to the manufacturer's instructions. RNA-seq libraries of polyadenylated RNA were prepared using the TruSeq RNA Library Prep Kit v2 (Illumina) or TruSeq Stranded mRNA Library Prep Kit (Illumina) according to the manufacturer's instructions. Sequencing libraries were sequenced on an Illumina NextSeq 500 platform. Adenovector transductions Approximately 5 3 10 5 A549 cells were transduced with Adenovectors purchased from Vector Biolabs at an MOI of 500 to induce expression of mCherry (Ad-mCherry) and human ACE2 (Ad-GFP-h-ACE2). 48 h post-transduction, efficient gene expression of delivered fluorescent proteins was confirmed by fluorescent microscopy using an EVOS M5000 Imaging System. Drug treatments Approximately 5 3 10 5 Vero E6 cells were infected with SARS-CoV-2 at an MOI of 0.05 in DMEM supplemented with 2% FBS, 4.5 g/L D-glucose, 4 mM L-glutamine, 10 mM Non-Essential Amino Acids, 1 mM Sodium Pyruvate and 10 mM HEPES. Vero E6 cells were treated with 100 units of universal Type-I IFNB 2 h post-infection. NHBE cells were treated with 100 units of human IFNb as indicated. Cells were harvested for RNA and protein analysis 24 hpi as described below. Approximately 2.5 3 10 5 A549 cells transduced with an ACE2 Adenovector were pre-treated with 500 nM Ruxolitinib or DMSO control for 1 h in infection media before infection with SARS-CoV-2 at an MOI of 2 for 24 h. Cells were harvested for protein analysis as described below. Western blot Protein was extracted from cells in Radioimmunoprecipitation assay (RIPA) lysis buffer containing 1X cOmplete Protease Inhibitor Cocktail (Roche) and 1X Phenylmethylsulfonyl fluoride (Sigma Aldrich) prior to safe removal from the BSL-3 facility. Samples were analyzed by SDS-PAGE and transferred onto nitrocellulose membranes. Proteins were detected using mouse monoclonal anti-Actin IRDye 800CW, 926-32210) and goat anti-rabbit (IRDye 680RD Ferrets were randomly assigned to the different treatment groups (naive, n = 2; SARS-CoV-2 infection, n = 6; influenza A virus (pH1N1) infection, n = 2; influenza A virus (H3N2) infection, n = 2). Both influenza A virus and SARS-CoV-2 infections of ferrets were performed simultaneously in the BSL-3 facility. For influenza A virus infections, all naive ferrets were infected intranasally with 10 5 PFU of influenza A/California/04/2009 (pH1N1) virus or 10 6 PFU of influenza A/Texas/71/2017 (H3N2) virus. Nasal washes were collected from anesthetized ferrets challenged with influenza A/California/ 04/2009 (pH1N1) virus on day 7 post infection and preserved at À80 C. Trachea were collected from euthanized ferrets challenged with influenza A/Texas/71/2017 (H3N2) virus on day 3 post infection and preserved at À80 C. For SARS-CoV-2 virus infections, all naive ferrets were infected with 5 3 10 4 PFU of SARS-CoV-2 isolate USA-WA1/2020 To diminish the noise introduced by variables inherent to the use of different cell types and systems, our differential expression analyses were always performed by matching each experimental condition with a corresponding mock treated sample that counted for the cell type, collection time, concurrent animal controls, etc. The raw sequencing data (fastq files) for the SARS-CoV-1 and MERS infections were downloaded from GEO (GSE56192, including their corresponding mock-treated controls) and processed in the same way as the rest of our experimental conditions. In order to capture the whole breadth of the response to IFNb treatment we pooled samples from 4, 6 and 12 hr post treatment and compare them together to mock treated cells. Differentially expressed genes (DEGs) were characterized for each sample (|L2FC| > 1, padjusted-value < 0.05) and were used as query to search for enriched biological processes (Gene ontology BP) and network analysis of protein interactions using STRING All RNA-Seq data performed in this paper can be found on the NCBI Gene Expression Omnibus (GEO) under accession number GSE147507. All non-RNA-seq statistical analyses were performed as Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact, Benjamin tenOever (benjamin.tenoever@mssm.edu).