key: cord-0869155-uukdrm6y
authors: Chen, Jennifer S.; Alfajaro, Mia Madel; Chow, Ryan D.; Wei, Jin; Filler, Renata B.; Eisenbarth, Stephanie C.; Wilen, Craig B.
title: Nonsteroidal Anti-inflammatory Drugs Dampen the Cytokine and Antibody Response to SARS-CoV-2 Infection
date: 2021-03-10
journal: J Virol
DOI: 10.1128/jvi.00014-21
sha: b3ad9c398f6588072746b5fcd50e8f8db0b127e9
doc_id: 869155
cord_uid: uukdrm6y

Identifying drugs that regulate severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and its symptoms has been a pressing area of investigation during the coronavirus disease 2019 (COVID-19) pandemic. Nonsteroidal anti-inflammatory drugs (NSAIDs), which are frequently used for the relief of pain and inflammation, could modulate both SARS-CoV-2 infection and the host response to the virus. NSAIDs inhibit the enzymes cyclooxygenase-1 (COX-1) and cyclooxygenase-2 (COX-2), which mediate the production of prostaglandins (PGs). Since PGs play diverse biological roles in homeostasis and inflammatory responses, inhibiting PG production with NSAIDs could affect COVID-19 pathogenesis in multiple ways, including (i) altering susceptibility to infection by modifying expression of angiotensin-converting enzyme 2 (ACE2), the cell entry receptor for SARS-CoV-2; (ii) regulating replication of SARS-CoV-2 in host cells; and (iii) modulating the immune response to SARS-CoV-2. Here, we investigate these potential roles. We demonstrate that SARS-CoV-2 infection upregulates COX-2 in diverse human cell culture and mouse systems. However, suppression of COX-2 by two commonly used NSAIDs, ibuprofen and meloxicam, had no effect on ACE2 expression, viral entry, or viral replication. In contrast, in a mouse model of SARS-CoV-2 infection, NSAID treatment reduced production of proinflammatory cytokines and impaired the humoral immune response to SARS-CoV-2, as demonstrated by reduced neutralizing antibody titers. Our findings indicate that NSAID treatment may influence COVID-19 outcomes by dampening the inflammatory response and production of protective antibodies rather than modifying susceptibility to infection or viral replication. IMPORTANCE Public health officials have raised concerns about the use of nonsteroidal anti-inflammatory drugs (NSAIDs) for treating symptoms of coronavirus disease 2019 (COVID-19). NSAIDs inhibit the enzymes cyclooxygenase-1 (COX-1) and cyclooxygenase-2 (COX-2), which are critical for the generation of prostaglandins—lipid molecules with diverse roles in homeostasis and inflammation. Inhibition of prostaglandin production by NSAIDs could therefore have multiple effects on COVID-19 pathogenesis. Here, we demonstrate that NSAID treatment reduced both the antibody and proinflammatory cytokine response to SARS-CoV-2 infection. The ability of NSAIDs to modulate the immune response to SARS-CoV-2 infection has important implications for COVID-19 pathogenesis in patients. Whether this occurs in humans and whether it is beneficial or detrimental to the host remains an important area of future investigation. This also raises the possibility that NSAIDs may alter the immune response to SARS-CoV-2 vaccination.

during murine coronavirus-induced encephalitis (22) . Therefore, reducing PGD 2 and PGE 2 levels with NSAID treatment could improve the induction of antiviral immunity and yet also promote hyperinflammatory responses. Conversely, NSAID treatment could have detrimental effects on resolution of infection by inhibiting production of PGI 2 , which is antiviral in respiratory syncytial virus infection (23) . Furthermore, NSAIDs may inhibit antibody production to SARS-CoV-2, which has been observed for other viruses, but the effect of this on disease severity is unclear as antibodies can be protective or pathogenic (24, 25) . Altogether, given the complex and sometimes conflicting roles of PGs, it is difficult to predict the overall effect of NSAIDs on the immune response to SARS-CoV-2 and, ultimately, disease outcome.

We therefore systematically assessed the effect of NSAIDs on SARS-CoV-2 infection and the immune response to SARS-CoV-2. We found that SARS-CoV-2 infection induced COX-2 expression in human cells and mice. However, suppression of COX-2 by two commonly used NSAIDs, ibuprofen and meloxicam, had no effect on ACE2 expression, viral entry, or viral replication. In a mouse model of SARS-CoV-2 infection, NSAID treatment impaired the production of proinflammatory cytokines and neutralizing antibodies but did not affect weight loss, viral burden, or activation of innate and adaptive immune cells in the lung. These results indicate that NSAID use in humans may affect COVID-19 pathogenesis by mitigating the inflammatory response and the production of protective antibodies rather than by directly influencing viral replication.

To determine the role of the COX-2 pathway in SARS-CoV-2 infection, we evaluated induction of PTGS2 (encoding COX-2) in human cells and mice. We found that SARS-CoV-2 infection of human lung cancer cell line Calu-3 led to significant upregulation of PTGS2 (Fig. 1A) . This is consistent with RNA sequencing (RNA-seq) data sets of SARS-CoV-2infected Calu-3 cells and ACE2-overexpressing A549 cells, another lung cancer cell line ( Fig. 1B and C) (26) . However, infection of human liver cancer cell line Huh7.5 did not lead to significant PTGS2 induction, demonstrating cell type specificity of PTGS2 induction by SARS-CoV-2 (Fig. 1D) .

We next assessed whether SARS-CoV-2 induces PTGS2 in a more physiologically relevant cell culture system. We cultured primary human bronchial epithelial cells (HBECs) for 28 days at an air-liquid interface, which supports pseudostratified mucociliated differentiation providing an in vitro model of airway epithelium (27) . We infected HBECs with SARS-CoV-2 at the apical surface of the culture and then performed singlecell RNA sequencing at 1, 2, and 3 days postinfection (dpi) (28) . As we previously reported that ciliated cells in air-liquid interface cultures are the major target of infection (28), we looked for PTGS2 induction in this cell type. Aggregating ciliated cells across the three time points, we found that infected ciliated cells expressed higher levels of PTGS2 compared to uninfected bystander ciliated cells (Fig. 1E) , indicating that PTGS2 is also induced by SARS-CoV-2 in a cell-intrinsic manner in ciliated cells, a physiologically relevant target cell.

To determine the relevance of these findings in vivo, we utilized transgenic mice expressing human ACE2 driven by the epithelial cell keratin 18 promoter (K18-hACE2) (29) . As SARS-CoV-2 does not efficiently interact with mouse ACE2 (4), human ACE2expressing mice are required to support SARS-CoV-2 infection (30) (31) (32) (33) (34) (35) . K18-hACE2 mice were initially developed as a model of SARS-CoV infection and have recently been demonstrated as a model of severe SARS-CoV-2 infection in the lung (29, 36) . We found that intranasal infection of K18-hACE2 mice with SARS-CoV-2 led to significant upregulation of Ptgs2 in the lung at multiple time points postinfection (Fig. 1F) , consistent with recent SARS-CoV-2-infected K18-hACE2 lung RNA-seq data (Fig. 1G) (36) . Taken together, these results demonstrate that SARS-CoV-2 infection induces PTGS2 in diverse in vitro and in vivo airway and lung systems, across multiple independent studies. These findings therefore suggest that COX-2 signaling may be a relevant pathway for regulating SARS-CoV-2 infection or the immune response to the virus.

NSAID treatment does not affect ACE2 expression in human cells and mice. We next explored whether inhibition of COX-2 could affect viral infection by regulating ACE2 expression, as has been reported in studies of diabetic rats and intestinal organoids (7, 8) . We utilized two NSAIDs, the nonselective COX-1/COX-2 inhibitor ibuprofen and the selective COX-2 inhibitor meloxicam, which are commonly used clinically. We determined the maximum nontoxic doses of ibuprofen and meloxicam on Calu-3 and Huh7.5 cells ( Fig. 2A and B ) and demonstrated that both drugs are bioactive in Calu-3 cells, as measured by reduced PGE 2 levels (Fig. 2C) . Treatment of Calu-3 or Huh7.5 cells with ibuprofen or meloxicam did not significantly affect ACE2 expression ( Fig. 2D and E). To test whether NSAID treatment affects Ace2 expression in diverse tissues in vivo, we treated C57BL/6 mice with therapeutic doses of ibuprofen and meloxicam (37) (38) (39) (40) , which did not lead to changes in Ace2 expression in the lung, heart, kidney, or ileum ( Fig. 2F to I). These data indicate that inhibition of the COX-2 pathway by NSAIDs does not affect ACE2 expression in multiple cell and tissue types in vitro or in vivo.

NSAID treatment does not affect SARS-CoV-2 entry or replication in vitro. To assess whether NSAID treatment functionally affects SARS-CoV-2 entry, we used a vesicular stomatitis virus (VSV) core expressing Renilla luciferase pseudotyped with the SARS-CoV-2 spike protein (SARS2-VSVpp). We used VSV glycoprotein (G) pseudovirus (G-VSVpp) as a control (41, 42) . Quantification of luciferase activity showed that PTGS2 expression in Calu-3 (B) and ACE2-overexpressing A549 (A549-ACE2) (C) cells following SARS-CoV-2 infection. The data are from GSE147507 (26) . (D) Huh7.5 cells were infected with SARS-CoV-2 at an MOI of 0.05. PTGS2 expression was measured at 2 dpi, normalized to ACTB. (E) HBECs were cultured at an air-liquid interface and then infected at the apical surface with 10 4 PFU of SARS-CoV-2. Cells were collected at 1, 2, and 3 dpi for single-cell RNA sequencing (scRNA-seq) (28) . A volcano plot of differentially expressed genes in infected versus bystander ciliated cells pooled from all time points is shown. PTGS2 is highlighted. (F) K18-hACE2 mice were infected intranasally with 1.2 Â 10 6 PFU of SARS-CoV-2. Ptgs2 expression in the lung was measured at 0, 2, 4, and 7 dpi. (G) Ptgs2 expression in the lung of K18-hACE2 mice following intranasal SARS-CoV-2 infection. The data are from GSE154104 (36) . All data points in this figure are presented as means 6 the standard errors of the mean (SEM). Data were analyzed by Welch's two-tailed, unpaired t test (A, D, and F); Student two-tailed, unpaired t test (B, C, and G); and two-sided Mann-Whitney U test with continuity and Benjamini-Hochberg correction (E). *, P , 0.05; ***, P , 0.001; ****, P , 0.0001. Data in panels A and D are representative of two independent experiments with three replicates per condition. ACE2 expression was measured and normalized to ACTB. (F to I) C57BL/6 mice were treated intraperitoneally with DMSO, 30 mg/kg ibuprofen, or 1 mg/kg meloxicam daily for 4 days. Ace2 expression was measured in the lung (F), heart (G), kidney (H), and ileum (I), normalized to Actb. All data points in this figure are presented as means 6 the SEM. Data were analyzed by Welch's two-tailed, unpaired t test (C to I). **, P , 0.01; ns, not significant. Data in panels A to E are representative of two independent experiments with three replicates per condition; data in panels F to I are pooled from two independent experiments with a total of four to six mice per condition. pretreatment of Calu-3 or Huh7.5 cells with ibuprofen or meloxicam did not significantly affect SARS2-VSVpp or G-VSVpp entry ( Fig. 3A and B) , confirming that NSAID inhibition of COX-2 does not impact susceptibility to infection.

Next, we studied whether COX-2 inhibition affects SARS-CoV-2 replication. Viruses from several different families have been shown to induce COX-2 signaling in host cells, which can have either proviral or antiviral functions (9, 43, 44) . To this end, we utilized a replication-competent SARS-CoV-2 expressing a mNeonGreen reporter (icSARS-CoV-2-mNG) to study the effect of COX-2 inhibition by NSAIDs on viral replication (45) . We assessed icSARS-CoV-2-mNG replication in Calu-3 cells, which upregulate PTGS2 in response to SARS-CoV-2 infection ( Fig. 1A and B) , and Huh7.5 cells, which do not (Fig.  1D ). By quantifying the percentage of mNeonGreen-expressing cells, we found that treatment of Calu-3 or Huh7.5 cells with ibuprofen or meloxicam did not impact icSARS-CoV-2-mNG replication ( Fig. 3C and D) . These results indicate that SARS-CoV-2 induction of the COX-2 pathway in Calu-3 human lung cells does not regulate viral replication.

NSAID treatment does not affect SARS-CoV-2-induced weight loss or lung viral burden in mice. Since NSAIDs do not directly affect SARS-CoV-2 entry or replication in vitro, we next investigated whether they regulate SARS-CoV-2 replication and immunity in vivo. We treated K18-hACE2 mice with meloxicam 1 day prior to SARS-CoV-2 infection and continued daily meloxicam treatment throughout the course of the infection. Dimethyl sulfoxide (DMSO) control-and meloxicam-treated K18-hACE2 mice experienced similar weight loss beginning at 5 dpi (Fig. 4A ). In addition, viral lung burden at 6 dpi was similar between infected DMSO-and meloxicamtreated mice (Fig. 4B) , suggesting that meloxicam does not affect viral replication in vitro or in vivo. NSAID treatment does not affect innate or adaptive immune cell activation in the lungs of SARS-CoV-2-infected mice. As NSAIDs regulate inflammation, we next assessed whether NSAIDs perturb the host immune response to SARS-CoV-2 in mice. First, we characterized the abundance and activation status of innate and adaptive immune cell types in the lung at 6 dpi (gating strategy is shown in Fig. S1A and B in the supplemental material). SARS-CoV-2 infection led to a slight decrease in alveolar macrophage abundance and upregulation of activation marker CD86. However, this was not modified by meloxicam treatment (Fig. 5A and B). Neutrophil counts were not significantly altered by infection or meloxicam treatment at this time point (Fig. 5C ). Activated CD69 1 NK cells increased with infection in both DMSO and meloxicamtreated mice (Fig. 5D ). Ly6C 1 monocyte/macrophage numbers and expression of activation markers CD86, MHCII, and CD64 all increased following SARS-CoV-2 infection but were not affected by meloxicam treatment (Fig. 5E to H). Furthermore, activated (CD44 1 CD69 1 ) T cell populations, including CD4 1 T cells, CD8 1 T cells, and gd T cells, increased with infection but were not significantly affected by concomitant meloxicam treatment ( Fig. 6A to C). In contrast, B cells were decreased in the lungs of infected mice, independent of meloxicam treatment (Fig. 6D ). Together, these results indicate that meloxicam treatment does not affect the overall numbers and activation status of innate or adaptive immune cells in the lung during SARS-CoV-2 infection.

NSAID treatment impairs systemic neutralizing antibody responses to SARS-CoV-2. While meloxicam did not alter the abundance of immune cell populations in the lung, NSAID treatment might regulate the function of these cell types. To assess the effect of meloxicam treatment on antibody production in infected mice, we measured the serum levels of spike-specific IgM and IgG antibodies at 6 dpi. DMSO-treated mice demonstrated detectable levels of both spike-specific IgM and IgG antibodies, which were significantly reduced in meloxicam-treated mice ( Fig. 6E and F ). Corroborating the difference in antibody titers, serum from DMSO-treated mice exhibited greater neutralization capacity compared to serum from meloxicamtreated mice (Fig. 6G) .

NSAID treatment dampens the induction of proinflammatory cytokines that are upregulated by SARS-CoV-2 infection in mice. Finally, we studied whether NSAID treatment modulates the cytokine response to SARS-CoV-2 infection. SARS-CoV-2 infection led to increased production of proinflammatory cytokines (IL-1b, IL-6, IFN-g, TNF-a, and granulocyte-macrophage colony-stimulating factor [GM-CSF]); T cell growth factors (IL-2); and chemokines (CCL2, CCL4, CXCL9, and CXCL10) that have been associated with disease severity and mortality in COVID-19 patients (Fig. 7A) (16, 17) . Among uninfected mice, meloxicam treatment led to minimal changes in cytokine 

Given the concerns about NSAID use in patients with COVID-19, we studied whether NSAIDs affect SARS-CoV-2 infection and the immune response to the virus. We found that SARS-CoV-2 infection induced PTGS2 upregulation in diverse systems, including Calu-3 and A549 lung cancer cell lines, primary HBEC air-liquid interface cultures, and the lungs of K18-hACE2 mice. Inhibition of COX-2 with the commonly used NSAIDs ibuprofen and meloxicam did not affect ACE2 expression in multiple cell and tissue types in vitro or in vivo, nor did it affect SARS-CoV-2 entry or replication. In K18-hACE2 mice infected with SARS-CoV-2, NSAID treatment impaired the production of neutralizing antibodies and proinflammatory cytokines but did not influence weight loss, viral burden, or activation state of innate and adaptive immune cells in the lung. Our findings therefore rule out a direct effect of NSAIDs on SARS-CoV-2 infection but indicate that NSAIDs could modulate COVID-19 severity by dampening production of neutralizing antibodies and inflammatory cytokines.

An important question arising from our findings is how SARS-CoV-2 infection induces COX-2 expression. One possibility is that the pattern recognition receptor retinoic acid inducible gene-I (RIG-I), which can recognize double-stranded RNA generated during viral genome replication and transcription (46) , may drive this response. Indeed, (47), and we showed here that Huh7.5 cells, which are defective in RIG-I signaling (48), do not upregulate PTGS2 in response to SARS-CoV-2. Alternatively, SARS-CoV-2 proteins may mediate the induction of COX-2 through their complex effects on host cells. In the case of SARS-CoV, transfection of plasmids encoding either the spike or the nucleocapsid genes is sufficient to stimulate COX-2 expression (12, 13) . SARS-CoV spike protein induces COX-2 expression through both calcium-dependent PKCa/ERK/NF-κB and calcium-independent PI3K/ PKC« /JNK/CREB pathways (13) , while the nucleocapsid protein directly binds to the COX-2 promoter to regulate its expression (12) . Any of these potential mechanisms are consistent with our HBEC scRNA-seq results demonstrating that SARS-CoV-2 increases PTGS2 expression in a cell-intrinsic manner.

One of the effects of NSAID treatment on the immune response to SARS-CoV-2 was impairment of early, neutralizing spike-specific antibodies. These early humoral responses are mediated by short-lived plasmablasts, and their requirement for T cell help is unclear (49) . NSAID treatment could therefore act by inhibiting activation (T cell-dependent or T cell-independent), proliferation, differentiation, or antibody-secreting capacity of spike-specific B cells and plasmablasts. B cells have been found to upregulate COX-2 expression following activation with T cell-dependent and T cell-independent stimuli (50, 51) , and treatment of purified B cell cultures with NSAIDs reduces IgM and IgG production (52). COX-2 inhibition reduces expression of BLIMP-1 and XBP-1 (53), which are essential transcription factors for plasmablast differentiation, providing a potential mechanism by which NSAIDs impair antibody production. In addition, in mice infected with vaccinia virus, antiviral antibody production is impaired by chronic but not acute COX-2 inhibition (24) , indicating that the duration of NSAID treatment may regulate the impact on antibody production. Given that K18-hACE2 mice succumb to lethal SARS-CoV-2 disease within 7 days (36, 54), we could not assess the effect of NSAID treatment on long-term antibody production. This could be explored in nonlethal models of SARS-CoV-2 infection, as well as in the setting of vaccination.

Understanding the effect of NSAID treatment on cytokine production is also critical, as cytokines may be protective early in COVID-19 but potentially pathological at later stages, thus informing the timing of immunomodulatory drugs like NSAIDs. We observed that NSAID treatment decreased the production of a subset of cytokines that were induced by infection, including IL-6, CCL2, GM-CSF, CXCL10, IL-2, and TNF-a. This is consistent with prior reports of COX-2 inhibitors decreasing the production of these cytokines in various inflammatory settings (55) (56) (57) (58) (59) . Mechanistically, NSAIDs may reduce production of these cytokines through inhibition of the PGE 2 /NF-κB positive-feedback loop, in which NF-κB and COX-2 can reciprocally activate their respective signaling pathways and amplify inflammatory responses (18, 19, 60) .

However, it is less clear whether this dampened cytokine response is beneficial, detrimental, or neutral in the setting of COVID-19 given the many roles that cytokines can play in controlling infection or driving immunopathology. GM-CSF is a myelopoietic growth factor, as well as proinflammatory cytokine, with pathogenic GM-CSF-producing Th1 cells being reported in patients with severe COVID-19 (61) . Reduction of GM-CSF production by meloxicam could therefore indicate a beneficial effect of NSAIDs on restraining hyperinflammatory responses. IL-2, IFN-g, TNF-a can be coproduced by polyfunctional T cells, which play important roles in the control of viral infections, and yet these cytokines can also promote lethal cytokine shock and tissue damage (62) (63) (64) . In addition, while IL-6 is correlated with disease severity in COVID-19, clinical trials blocking IL-6 signaling have not shown clear evidence of benefit for patients (65) . While we observed decreased cytokine responses with NSAID treatment in K18-hACE2 mice, these changes did not translate into differences in weight loss or viral burden, suggesting that these features of the disease may involve other cytokines or pathology at other sites (e.g., the brain) not affected by NSAID treatment (36, 54) . The timing, duration, and dosing of NSAID treatment may also matter for COVID-19 pathogenesis, potentially with early treatment impacting the initiation of antiviral immune responses and later treatment suppressing immune-driven pathology (14) . In the present study, we treated K18-hACE2 mice continuously throughout the course of infection, but it would be interesting to explore in future work whether limiting NSAID treatment to particular phases of disease has differential effects on pathogenesis.

Here, we demonstrated that SARS-CoV-2 infection induces COX-2 expression in cell lines, primary airway epithelial cells, and mice. Inhibition of COX-2 by NSAIDs did not affect viral entry or replication in vitro or in vivo. However, NSAID treatment impaired the production of proinflammatory cytokines and neutralizing antibodies in response to SARS-CoV-2 infection in mice. NSAIDs could therefore have complex effects on the host response to SARS-CoV-2. While studies thus far have not observed worse clinical outcomes in COVID-19 patients taking NSAIDs (66) (67) (68) (69) (70) (71) , evaluation of the breadth, potency, and durability of the humoral immune response is warranted in patients on NSAIDs in response to both natural infection and vaccination.

Cell lines. Calu-3 and Huh7.5 were from ATCC. Calu-3 cells were cultured in Eagle minimum essential medium with 10% heat-inactivated fetal bovine serum (FBS), 1% GlutaMAX (Gibco), and 1% penicillin/ streptomycin. Huh7.5 cells were cultured in Dulbecco modified Eagle medium (DMEM) with 10% heatinactivated FBS and 1% penicillin/streptomycin. All cell lines tested negative for Mycoplasma spp.

Generation of SARS-CoV-2 stocks. As previously described (42) , SARS-CoV-2 P1 stock was generated by inoculating Huh7.5 cells with SARS-CoV-2 isolate USA-WA1/2020 (BEI Resources, NR-52281) at a multiplicity of infection (MOI) of 0.01 for 3 days. The P1 stock was then used to inoculate Vero-E6 cells, and after 3 days, the supernatant was harvested and clarified by centrifugation (450 Â g for 5 min), filtered through a 0.45-mm filter, and stored in aliquots at -80°C. Virus titer was determined by plaque assay using Vero-E6 cells (42) .

To generate icSARS-CoV-2-mNG stocks (45) , lyophilized icSARS-CoV-2-mNG was reconstituted in 0.5 ml of deionized water. Then, 50 ml of virus was diluted in 5 ml of medium and added to 10 7 Vero-E6 cells. After 3 days, the supernatant was harvested and clarified by centrifugation (450 Â g for 5 min), filtered through a 0.45-mm filter, and stored in aliquots at -80°C.

All work with SARS-CoV-2 or icSARS-CoV-2-mNG was performed in a biosafety level 3 facility with approval from the office of Environmental Health and Safety and the Institutional Animal Care and Use Committee at Yale University.

Preparation of NSAIDs. Ibuprofen (I4883) and meloxicam (M3935) were purchased from Sigma-Aldrich. For cell culture experiments, ibuprofen and meloxicam were solubilized in DMSO at a stock concentration of 10 mM and then diluted in medium to make working solutions. For mouse experiments, stock solutions of ibuprofen (300 mg/ml) and meloxicam (10 mg/ml) were prepared in DMSO and then diluted 100-fold in phosphate-buffered saline (PBS) to make working solutions. To determine the maximum nontoxic dose of NSAIDs to use for cell culture experiments, cells were treated with different concentrations of ibuprofen or meloxicam for 48 h, and cell viability was measured by a CellTiter-Glo luminescent cell viability assay (Promega) according to the manufacturer's instructions.

Mice. C57BL/6J and K18-hACE2 [B6.Cg-Tg(K18-ACE2)2Prlmn/J (29)] were purchased from Jackson Laboratory. Mice were bred in-house using mating trios to enable utilization of littermates for experiments. Mice of both sexes between 6 and 8 weeks old were used for this study. C57BL/6J and K18-hACE2 mice were anesthetized using 30% (vol/vol) isoflurane diluted in propylene glycol (30% isoflurane) and administered 30 mg/kg ibuprofen, 1 mg/kg meloxicam, or an equivalent amount of DMSO intraperitoneally in a volume of 10 ml/kg daily for 4 or 7 days as indicated in the figure legends. K18-hACE2 mice were anesthetized using 30% isoflurane and administered 1.2 Â 10 6 PFU or 10 3 PFU of SARS-CoV-2 intranasally as indicated in the figure legends. Mice were monitored daily for weight and survival. Animal use and care was approved in agreement with the Yale Animal Resource Center and Institutional Animal Care and Use Committee (no. 2018-20198) according to the standards set by the Animal Welfare Act.

Analysis of RNA-seq data. We utilized RNA-seq data from recent published studies to assess the impact of SARS-CoV-2 infection on PTGS2 expression. From GSE147507 (26), we reanalyzed the raw count data from Calu-3 and A549-ACE2 cells, comparing SARS-CoV-2 infection to matched mock-treated controls. We performed differential expression analysis using the Wald test from DESeq2 (72), using a Benjamini-Hochberg adjusted P , 0.05 as the cutoff for statistical significance. For visualization of PTGS2 expression, the DESeq2-normalized counts were exported and plotted in GraphPad Prism. Statistical significance was assessed using a Student two-tailed, unpaired t test.

For analysis of HBEC air-liquid interface cultures infected with SARS-CoV-2, we utilized a previously generated catalog of differentially expressed genes that our group recently described in a preprint study (28) . The differential expression table is publicly available (https://github.com/vandijklab/HBEC _SARS-CoV-2_scRNA-seq). Here, we specifically investigated PTGS2 expression in ciliated cells, comparing infected cells to bystander cells (cells aggregated across the 1, 2, and 3 dpi time points). The cutoff for statistical significance was set at adjusted P , 0.05, and the results were visualized as a volcano plot in R.

From GSE154104 (36), we reanalyzed the raw count data from the lungs of K18-hACE2 mice infected with SARS-CoV-2, performing pairwise comparisons of mice at 2 dpi, 4 dpi, and 7 to 0 dpi controls (prior to infection). For visualization of Ptgs2 expression, the DESeq2-normalized counts were exported and plotted in GraphPad Prism. Statistical significance was assessed using a Student two-tailed, unpaired t test.

PGE 2 ELISA. Levels of PGE 2 in cell culture supernatants were measured using the Prostaglandin E 2 enzyme-linked immunosorbent assay (ELISA) kit (Cayman Chemical) according to the manufacturer's instructions. Absorbance was measured at 410 nm on a microplate reader (Molecular Devices), and PGE 2 concentrations were calculated using a standard curve.

Quantitative PCR. Cells or tissues were lysed in TRIzol (Life Technologies), and total RNA was extracted using the Direct-zol RNA Miniprep Plus kit (Zymo Research) according to the manufacturer's instructions. cDNA synthesis was performed using random hexamers and ImProm-II reverse transcriptase (Promega). qPCR was performed with Power SYBR Green (Thermo Fisher) and run on the QuantStudio3 (Applied Biosystems). Target mRNA levels were normalized to those of ACTB or Actb. The qPCR primer sequences were as follows: ACTB (human), GAGCACAGAGCCTCGCCTTT (forward) and ATCATCATCCATGGTGAGCTGG (reverse); PTGS2 (human), AGAAAACTGCTCAACACCGGAA (forward) and GCACTGTGTTTGGAGTGGGT (reverse); ACE2 (human), GGGATCAGAGATCGGAAGAAGAAAA (forward) and AAGGAGGTCTGAACA TCATCAGTG (reverse); Actb (mouse), ACTGTCGAGTCGCGTCCA (forward) and ATCCATGGCGAACTG GTGG (reverse); Ptgs2 (mouse), CTCCCATGGGTGTGAAGGGAAA (forward) and TGGGGGTCAGGGATG AACTC (reverse); and Ace2 (mouse), ACCTTCGCAGAGATCAAGCC (forward) and CCAGTGGGGCTGATGT AGGA (reverse).

Pseudovirus production. VSV-based pseudotyped viruses were produced as previously described (41, 42) . Vector pCAGGS containing the SARS-related coronavirus 2, Wuhan-Hu-1 Spike glycoprotein gene, NR-52310, was produced under HHSN272201400008C and obtained through BEI Resources, NIAID, NIH. 293T cells were transfected with the pCAGGS vector expressing the SARS-CoV-2 spike glycoprotein and then incubated with replication-deficient VSV expressing Renilla luciferase for 1 h at 37°C (41). The virus inoculum was then removed, and the cells were washed with PBS before adding media with anti-VSV-G clone I4 to neutralize residual inoculum. No antibody was added to cells expressing VSV-G. Supernatant containing pseudoviruses was collected 24 h postinoculation, clarified by centrifugation, and stored in aliquots at -80°C.

Pseudovirus entry and neutralization assays. We plated 3 Â 10 4 Calu-3 or 1 Â 10 4 Huh7.5 cells in a 100-ml volume in each well of a 96-well plate. The following day, the medium was replaced with 50 mM ibuprofen, 50 mM meloxicam, or an equivalent amount of DMSO. One day later, 10 ml of SARS-CoV-2 spike protein-pseudotyped (SARS2-VSVpp) or VSV glycoprotein-typed virus was added. The luciferase activity was measured at 24 hpi using the Renilla luciferase assay system (Promega). Each well of cells was lysed with 50 ml of lysis buffer, and 15 ml of cell lysate was then mixed with 15 ml of luciferase assay reagent. The luminescence was measured on a microplate reader (BioTek Synergy).

For neutralization assays, 2 Â 10 4 Huh7.5 cells were plated in a 100-ml volume in each well of a 96well plate. The following day, serial dilutions of serum were incubated with SARS2-VSVpp pseudovirus for 1 h at 37°C. The growth medium was then aspirated from the cells and replaced with 50 ml of the serum/virus mixture. Luciferase activity was measured at 24 hpi as detailed above. Half-maximal inhibitory concentrations (IC 50 s) were calculated as previously described (73) .

icSARS-CoV-2-mNG assay. We plated 6.5 Â 10 3 Calu-3 or 2.5 Â 10 3 Huh7.5 cells in 20 ml of phenol red-free medium containing 50 mM ibuprofen, 50 mM meloxicam, or an equivalent amount of DMSO in each well of a black-walled, clear-bottom 384-well plate. The following day, icSARS-CoV-2-mNG was added at an MOI of 1 in a 5-ml volume. The frequency of infected cells was measured by mNeonGreen expression at 1, 2, and 3 dpi by high content imaging (BioTek Cytation 5) configured with brightfield and GFP cubes. The total cell numbers were quantified by Gen5 software for brightfield images. Object analysis was used to determine the number of mNeonGreen-positive cells. The percentage of infection was calculated as the ratio of the number of mNeonGreen-positive cells to the total number of cells in brightfield.

Measurement of lung viral burden and cytokines. Lungs were perfused with 3 ml of sterile PBS. The left lobe was collected and homogenized in 1 ml of DMEM supplemented with 2% heat-inactivated FBS and 1% antibiotic-antimycotic. The viral burden was measured in lung homogenates by plaque assay on Vero-E6 cells as previously described (42) . To measure cytokines, lung homogenates were incubated with a final concentration of 1% Triton X-100 for 1 h at room temperature to inactivate SARS-CoV-2. Cytokine analysis was performed by Eve Technologies using their Mouse Cytokine Array/Chemokine Array 31-Plex (MD31) platform.

Flow cytometry analysis of lung immune cells. Lungs were perfused with 3 ml of sterile PBS. The right inferior lobe was collected and digested with 0.5 mg/ml collagenase IV (Sigma-Aldrich) and 100 U/ ml DNase I (Sigma-Aldrich) in complete RPMI for 45 min at 37°C. Single-cell suspensions of digested lung tissue were preincubated with Fc block (clone 2.4G2) for 5 min at room temperature before staining. The cells were stained with the following antibodies or viability dyes for 30 min at 4°C: PE anti-CD64 (clone X54-5/7.1), PE/Cy7 anti-Ly6C (clone HK1.4), PerCP/Cy5.5 anti-CD45.2 (clone 104), APC anti-CD86 (clone GL1), AF700 anti-CD19 (clone 6D5), DAPI (49,69-diamidino-2-phenylindole; Thermo Fisher), BV510 anti-I-A/I-E (clone M5/114.15.2), BV605 anti-CD11c (clone N418), Live/Dead Fixable Green (Thermo Fisher), PE anti-TCRg/d (clone GL3), PE/Cy7 anti-CD69 (clone H1.2F3), PerCP/Cy5.5 anti-CD8 (clone 53.6-7), APC anti-CD45.2 (clone 104), APC/Cy7 anti-TCRb (clone H57-597), AF700 anti-CD4 (clone RM4-5), PB anti-Ly6G (clone 1A8), BV510 anti-NK1.1 (clone PK136), and BV605 anti-CD44 (clone IM7). After being washed, the cells were fixed with 4% paraformaldehyde for 30 min at room temperature to inactivate SARS-CoV-2. Samples were acquired on a CytoFLEX S (Beckman Coulter) and analyzed using FlowJo software (BD).

Spike-specific ELISAs. Serum was incubated with a final concentration of 0.5% Triton X-100 and 0.5 mg/ml RNase A to inactivate any potential SARS-CoV-2. SARS-CoV-2 stabilized spike glycoprotein (BEI Resources, NR-53524) was coated at a concentration of 2 mg/ml in carbonate buffer on 96-well MaxiSorp plates (Thermo Fisher) overnight at 4°C. Plates were blocked with 1% BSA in PBS for 1 h at room temperature. Serum samples were serially diluted in 1% BSA in PBS and incubated in plates for 2 h at room temperature. Antibody isotypes were detected with anti-mouse IgM-HRP or anti-mouse IgG Fc-HRP (Southern Biotech) by incubation for 1 h at room temperature. The plates were developed with TMB stabilized chromogen (Thermo Fisher), stopped with 3 N hydrochloric acid, and read at 450 nm on a microplate reader.

Statistical analysis. Data analysis was performed using GraphPad Prism 8 unless otherwise indicated. Data were analyzed using Welch's two-tailed, unpaired t test; Student two-tailed, unpaired t test; two-tailed Mann-Whitney test; or two-way ANOVA, as indicated. P , 0.05 was considered statistically significant.

Data availability. All previously published data are available as described above. Cytokine data generated in this study are available in Table S1 in the supplemental material.

Supplemental material is available online only. SUPPLEMENTAL FILE 1, XLSX file, 0.02 MB. SUPPLEMENTAL FILE 2, PDF file, 0.3 MB.

Covid-19: ibuprofen should not be used for managing symptoms, say doctors and scientists

Novel coronavirus: anti-inflammatory medications. Medicines and Healthcare Products Regulatory Agency

Prostaglandins and inflammation

A pneumonia outbreak associated with a new coronavirus of probable bat origin

Are patients with hypertension and diabetes mellitus at increased risk for COVID-19 infection?

Cyclooxygenase activity is important for efficient replication of mouse hepatitis virus at an early stage of infection

Ibuprofen attenuates cardiac fibrosis in streptozotocininduced diabetic rats

Prostaglandin E 2 promotes intestinal repair through an adaptive cellular response of the epithelium

Activation of COX-2/PGE2 promotes sapovirus replication via the inhibition of nitric oxide production

Feline calicivirus-and murine norovirus-induced COX-2/PGE 2 signaling pathway has proviral effects

Genomic characterization and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding

Nucleocapsid protein of SARS-CoV activates the expression of cyclooxygenase-2 by binding directly to regulatory elements for nuclear factor-κB and CCAAT/ enhancer binding protein

Spike protein of SARS-CoV stimulates cyclooxygenase-2 expression via both calcium-dependent and calcium-independent protein kinase C pathways

Cytokine storm

RECOVERY Collaborative Group. 2020. Dexamethasone in hospitalized patients with Covid-19: preliminary report

Longitudinal analyses reveal immunological misfiring in severe COVID-19. 7821

Longitudinal immune profiling reveals key myeloid signatures associated with COVID-19

Non-steroidal anti-inflammatory drugs, prostaglandins, and COVID-19

Prostaglandin-cytokine crosstalk in chronic inflammation

Age-related increases in PGD 2 expression impair respiratory DC migration, resulting in diminished T cell responses upon respiratory virus infection in mice

Targeted prostaglandin E2 inhibition enhances antiviral immunity through induction of type I interferon and apoptosis in macrophages

Virus-induced inflammasome activation is suppressed by prostaglandin D2/DP1 signaling

Signaling through the prostaglandin I2 receptor IP protects against respiratory syncytial virus-induced illness

Chronic inhibition of cyclooxygenase-2 attenuates antibody responses against vaccinia infection

A perspective on potential antibody-dependent enhancement of SARS-CoV-2. 7821

Imbalanced host response to SARS-CoV-2 drives development of COVID-19

Structural and functional variations in human bronchial epithelial cells cultured in air-liquid interface using different growth media

Single-cell longitudinal analysis of SARS-CoV-2 infection in human airway epithelium

Lethal infection of K18-hACE2 mice infected with severe acute respiratory syndrome coronavirus

The pathogenicity of SARS-CoV-2 in hACE2 transgenic mice

Pathogenesis of SARS-CoV-2 in transgenic mice expressing human angiotensin-converting enzyme 2

A SARS-CoV-2 infection model in mice demonstrates protection by neutralizing antibodies

Generation of a broadly useful model for COVID-19 pathogenesis, vaccination, and treatment

Mouse model of SARS-CoV-2 reveals inflammatory role of type I interferon signaling

A mouse model of SARS-CoV-2 infection and pathogenesis

SARS-CoV-2 infection of human ACE2-transgenic mice causes severe lung inflammation and impaired function

Co-administered low doses of ibuprofen and dexamethasone produce synergistic antinociceptive effects on neuropathic mechanical allodynia in rats

Effects of buprenorphine, meloxicam, and flunixin meglumine as postoperative analgesia in mice

Antinociceptive effect of meloxicam, in neurogenic and inflammatory nociceptive models in mice

Analgesic activity of the novel COX-2 preferring NSAID, meloxicam in mono-arthritic rats: central and peripheral components

A structural basis for antibody-mediated neutralization of Nipah virus reveals a site of vulnerability at the fusion glycoprotein apex

Genomewide CRISPR screens reveal host factors critical for SARS-CoV-2 infection

Prostaglandin E2 as a modulator of viral infections

Inhibition of multiplication of parainfluenza 3 virus in prostaglandin-treated WISH cells

An infectious cDNA clone of SARS-CoV-2

Interaction of SARS and MERS coronaviruses with the antiviral interferon response

Influenza A viruses suppress cyclooxygenase-2 expression by affecting its mRNA stability

Regulating intracellular antiviral defense and permissiveness to hepatitis C virus RNA replication through a cellular RNA helicase, RIG-I

Origin and function of circulating plasmablasts during acute viral infections

Activated human B lymphocytes express cyclooxygenase-2 and cyclooxygenase inhibitors attenuate antibody production

CpG oligodeoxynucleotides induce cyclooxygenase-2 in human B lymphocytes: implications for adjuvant activity and antibody production

Ibuprofen and other widely used non-steroidal anti-inflammatory drugs inhibit antibody Effect of NSAIDs on SARS-CoV-2 Infection Journal of Virology production in human cells

Inhibition of cyclooxygenase-2 impairs the expression of essential plasma cell transcription factors and human B-lymphocyte differentiation

COVID-19 treatments and pathogenesis including anosmia in K18-hACE2 mice

Cyclooxygenase-2 inhibitors decrease interleukin-1b-stimulated prostaglandin E2 and IL-6 production by human gingival fibroblasts

Cyclooxygenase-2 promotes hepatocellular apoptosis by interacting with TNF-a and IL-6 in the pathogenesis of nonalcoholic steatohepatitis in rats

Diclofenac, a selective COX-2 inhibitor, inhibits DMH-induced colon tumorigenesis through suppression of MCP-1, MIP-1a, and VEGF

Nonsteroidal anti-inflammatory drugs may affect cytokine response and benefit healing of combatrelated extremity wounds

Selective inhibition of cyclooxygenase (COX)-2 reverses inflammation and expression of COX-2 and interleukin 6 in rat adjuvant arthritis

Prostaglandin E2-EP2-NF-κB signaling in macrophages as a potential therapeutic target for intracranial aneurysms

2020. Pathogenic T cells and inflammatory monocytes incite inflammatory storms in severe COVID-19 patients

Emergence of CD4 1 and CD8 1 polyfunctional T cell responses against immunodominant lytic and latent EBV antigens in children with primary EBV infection

Tumor necrosis factor alpha exerts powerful anti-influenza virus effects in lung epithelial cells

Synergism of TNF-a and IFN-g triggers inflammatory cell death, tissue damage, and mortality in SARS-CoV-2 infection and cytokine shock syndromes

Tocilizumab for Covid-19: the ongoing search for effective therapies

Older age and comorbidity are independent mortality predictors in a large cohort of 1305 COVID-19 patients in Michigan, United States

Bar-Haim A, Youngster I. 2020. Ibuprofen use and clinical outcomes in COVID-19 patients

OpenSAFELY: do adults prescribed Non-steroidal anti-inflammatory drugs have an increased risk of death from COVID-19

Ibuprofen and NSAID use in COVID-19 infected patients is not associated with worse outcomes: a prospective cohort study

Prior routine use of non-steroidal anti-inflammatory drugs (NSAIDs) and important outcomes in hospitalized patients with COVID-19

Adverse outcomes and mortality in users of non-steroidal anti-inflammatory drugs who tested positive for SARS-CoV-2: a Danish nationwide cohort study

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

Antigenic conservation and immunogenicity of the HIV coreceptor binding site

We acknowledge Benhur Lee, Pei-Yong Shi, the World Reference Center for Emerging Viruses and Arboviruses (WRCEVA), Paulina Pawlica, Joan Steitz, Michael Diamond, Adam Bailey, Emma Winkler, and BEI Resources for providing critical reagents and expertise. We thank all members of the Wilen and Eisenbarth labs for helpful discussions. We thank Yale Environmental Health and Safety for providing necessary training and support for SARS-CoV-2 research.This study was supported by NIH Yale University (C.B.W.) has a patent pending related to this work entitled "Compounds and Compositions for Treating, Ameliorating, and/or Preventing SARS-CoV-2 Infection and/or Complications Thereof." Yale University has committed to rapidly executable nonexclusive royalty-free licenses to intellectual property rights for the purpose of making and distributing products to prevent, diagnose, and treat COVID-19 infection during the pandemic and for a short period thereafter.