Elevated temperature inhibits SARS-CoV-2

Full text

Turn on search term navigation

Introduction

The pandemic spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2; [1–3]) is an ongoing international health crisis with over 245 million infections and 5 million reported deaths worldwide to date (WHO; https://covid19.who.int; November 2021). The spectrum of SARS-CoV-2–related disease (Coronavirus Disease 2019 (COVID-19)) is highly variable, ranging from asymptomatic viral shedding to acute respiratory distress syndrome (ARDS), multiorgan failure, and death. Besides coughing, dyspnea, and fatigue, fever (also known as pyrexia) is one of the most frequently reported symptoms [4–8]. Fever is an evolutionarily conserved response to microbial infection and inflammation that can influence the regulation of multiple cellular processes, including host innate and adaptive immune responses [9,10]. Unlike hyperthermia or heat stroke, fever represents a controlled shift in body temperature regulation induced by the expression of exogenous (microbial) and endogenous (host) pyrogenic regulatory factors, including pathogen associated molecular patterns (PAMPs) and proinflammatory cytokines (e.g., interleukin 6 (IL-6)) [9,10]. Body temperature naturally varies throughout the day, with age, sex, and ethnic origin being contributing factors [9,10]. In healthy middle-aged adults, febrile temperatures range from 38 to 41°C (ΔT approximately 1 to 4°C above baseline), with low (38 to 39°C), moderate (39.1 to 40°C), high (40.1 to 41.1°C), and hyperpyrexia (>41.1°C) temperature ranges [10]. Temperature elevation is known to confer protection against a number of respiratory pathogens [9,11], with antipyretic drug treatment leading to increased mortality in intensive care unit (ICU) patients infected with influenza A virus (IAV) [12–14].

With respect to COVID-19, up to 90% of hospitalized patients show low (44%) to moderate (13% to 34%) grade fever [4–8], with ICU patients presenting a 10% higher prevalence of fever relative to non-ICU patients [6,15]. In vitro studies have shown that SARS-CoV-2 replicates more efficiently at lower temperatures associated with the upper respiratory airway (33°C), which correlates with an overall weaker interferon (IFN)-mediated immune response to infection relative to lower respiratory airway (37°C) infection [16]. These data suggest that tissue temperature could be a significant factor in the host immune response to SARS-CoV-2 infection and COVID-19 disease progression. However, it remains to be determined what effect elevated temperature (≥37°C) has on SARS-CoV-2 replication. We therefore set out to determine the net effect of temperature elevation on SARS-CoV-2 infection within respiratory epithelial tissue.

Utilizing a three-dimensional (3D) respiratory model that closely mimics the tissue physiology and cellular tropism of SARS-CoV-2 infection observed in the respiratory airway [16–24], we demonstrate temperature elevation (≥39°C) to restrict the replication and propagation of SARS-CoV-2 independently of the induction of IFN-mediated antiviral immune defenses. We show that respiratory epithelium remains permissive to SARS-CoV-2 infection at temperatures up to 40°C, but to restrict the initiation of viral transcription leading to reduced levels of viral RNA (vRNA) replication and apical shedding of infectious virus. We identify temperature to play an important role in the differential regulation of multiple epithelial host responses to SARS-CoV-2 infection, including epigenetic, long noncoding RNA (lncRNA), and immunity-related pathways. Our data identify an important role for tissue temperature in the epithelial restriction of SARS-CoV-2 replication independently of canonical IFN-mediated antiviral immune defenses previously reported to restrict SARS-CoV-2 infection.

Results

Differentiation of primary human bronchial epithelial cells into ciliated respiratory epithelium supports SARS-CoV-2 replication in vitro

To establish a model suitable to study the effect of temperature on SARS-CoV-2 replication, we differentiated primary human bronchiolar epithelial (HBEp) cells isolated from 3 independent donors into stratified respiratory epithelium. Hematoxylin and eosin (H&E) and immunohistochemistry (IHC) staining of respiratory cultures demonstrated these tissues to contain a mixture of epithelial and goblet cells, with high levels of apical ciliation and angiotensin-converting enzyme 2 (ACE2) expression, the principal surface receptor of SARS-CoV-2 (Fig 1A and 1B) [25–28]. Infection of respiratory airway cultures with SARS-CoV-2 (England/02/2020; MOI 0.05, 104 plaque-forming unit (PFU)/Tissue) at 37°C demonstrated these tissues to support infection and viral replication, with intraepithelial and apical vRNA accumulation readily detected by in situ hybridization (ISH) by 72 h postinfection (Fig 1A and 1B). Notably, we observed discrete clusters of vRNA accumulation within respiratory epithelia (Fig 1A, arrows), indicative of localized areas of intraepithelial infection, propagation, and spread [29,30]. The overall morphology of the respiratory epithelium remained largely intact, with little shedding of ciliated cells from the epithelial surface (Fig 1A, 120 h). Measurement of genome copies by reverse transcription quantitative PCR (RT-qPCR) and infectious virus by TCID50 in apical washes collected over time demonstrated the linear phase of virus release to occur between 48 and 120 h (Fig 1C and 1D). A substantial drop in apical infectious titers was observed at 144 h postinfection (Fig 1D), consistent with multiple waves of virus release over prolonged periods of incubation in primary airway cultures [31].

[Figure omitted. See PDF.]

Fig 1. Differentiation of primary bronchial epithelial airway cultures supports SARS-CoV-2 replication.

Primary bronchial epithelial cells isolated from 3 independent donors (donor 1, male Caucasian aged 63 years; donor 2, Hispanic male aged 62 years; donor 3, Caucasian female aged 16 years; all nonsmokers) were seeded onto 6.5 mm transwells and differentiated under ALI conditions for ≥35 days. Ciliated respiratory cultures were mock treated (media only) or SARS-CoV-2 (SCV2; MOI 0.05, 104 PFU/Tissue) infected at 37°C for the indicated times (h). (A, B) Representative images of H&E, ACE2 IHC (brown), or SARS-CoV-2 RNA ISH (red) stained sections. Hematoxylin was used as a counter stain. Scale bars = 20 μm. (C) Genome copies per ml of SARS-CoV-2 in apical washes harvested over time as determined RT-qPCR. N ≥ 7 tissues per condition derived from a minimum of 3 independent biological experiments. Means and SD shown. (D) TCID50 assay measuring infectious viral load in apical washes harvested from SARS-CoV-2 infected tissues over time. Means and SD shown. N = 6 tissues per condition derived from a minimum of 3 independent biological experiments. Raw values presented in S9 Data. ACE2, angiotensin-converting enzyme 2; ALI, air–liquid interface; H&E, hematoxylin and eosin; IHC, immunohistochemistry; ISH, in situ hybridization; PFU, plaque-forming unit; RT-qPCR, reverse transcription quantitative PCR; SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2; vRNA, viral RNA.

https://doi.org/10.1371/journal.pbio.3001065.g001

Respiratory airway cultures induce a heat stress response at elevated temperature

While the heat stress response has been extensively investigated in two-dimensional (2D) cell culture model systems [32,33], this pathway remains poorly characterized in 3D respiratory epithelia under air–liquid interface (ALI). We therefore investigated the influence of temperature on our 3D respiratory model. Mock-treated respiratory cultures were incubated at 37, 39, or 40°C (representative of core body temperature and low to moderate grade febrile temperatures, respectively) prior to tissue fixation or RNA extraction and RNA sequencing (RNA-Seq). H&E staining and transepithelial electrical resistance (TEER) measurements demonstrated no morphological or membrane integrity changes to the respiratory epithelium upon temperature elevation (Fig 2A and 2B). RNA-Seq analysis of respiratory epithelia demonstrated no difference in the expression level of ACE2 (Fig 2C, p = 0.6753) or a reference set of 24 genes known to be constitutively expressed across a wide range of tissues and cell types (Fig 2D, p = 0.2136; [34–36]) upon temperature elevation. These data indicate that general transcription remains largely unperturbed at elevated temperatures up to 40°C. Out of the 867 differentially expressed genes (DEGs) identified (Fig 2E; Q < 0.05), DEG enrichment was observed across multiple pathways in tissues incubated at 40 relative to 37°C, including cellular response to stress, RNA polymerase I promoter opening, and DNA methylation (Fig 2F and S1 Data). Significant DEG enrichment was also observed in cellular response to heat stress (Fig 2F [arrow], 2G; S1 Data). We conclude that our respiratory model induces a significant heat stress response at elevated temperature without visible damage to the epithelium or induction of innate immune defenses that may otherwise influence the interpretation of SARS-CoV-2 replication studies at elevated temperature.

[Figure omitted. See PDF.]

Fig 2. Respiratory airway cultures induce a heat stress response upon incubation at elevated temperature.

Ciliated respiratory cultures differentiated from HBEp cells isolated from 3 independent donors (donor 1, male Caucasian aged 63 years; donor 2, Hispanic male aged 62 years; donor 3, Caucasian female aged 16 years; all nonsmokers) were incubated at 37, 39, or 40°C for 72 h prior to fixation or RNA extraction and RNA-Seq. (A) Representative H&E stained tissue sections (donor 1) at 72 h. Scale bars = 20 μm. (B) TEER (Ohm per cm2) measurements of respiratory epithelia integrity (donor 1) at 37 or 40°C. N = 3 tissues per condition; lines, mean; all data points shown; p-values shown, unpaired two-tailed t test. (C) ACE2 expression levels (CPM) in respiratory cultures incubated at 37 or 40°C; means and SD shown; p-value shown, unpaired two-tailed t test. (D) Expression values (log2 CPM) of a reference gene set (24 genes; [34–36]) in respiratory cultures incubated at 37 or 40°C; p-value shown, paired two-tailed t test. (E) Scatter plots showing high confidence (Q < 0.05) DEG transcripts identified between respiratory cultures incubated at 37 or 40°C; up-regulated DEGs, red circles; down-regulated DEGs, blue circles. (F) Reactome pathway analysis of mapped up-regulated DEGs. Top 30 pathways (FDR < 0.05; plotted as −log10 FDR) shown. Dotted line, threshold of significance (−log10 FDR of 0.05). (G) Expression values (log2 CPM) of DEGs associated with cellular response to heat stress pathway (R-HSA-3371556; arrow in F); p-value shown, paired two-tailed t test. (C to H) RNA-Seq data derived from RNA isolated from 3 donors (donors 1 to 3) per sample condition derived from 3 independent experiments per donor. Raw values presented in S1 and S9 Data. ACE2, angiotensin-converting enzyme 2; CPM, counts per million; DEG, differentially expressed gene; FDR, false discovery rate; HBEp, human bronchiolar epithelial; H&E, hematoxylin and eosin; TEER, transepithelial electrical resistance.

https://doi.org/10.1371/journal.pbio.3001065.g002

Elevated temperature restricts the replication of SARS-CoV-2 in respiratory airway cultures

We next examined the effect of elevated temperature on SARS-CoV-2 replication. Respiratory cultures were incubated at 37, 39, or 40°C for 24 h prior to SARS-CoV-2 infection and continued incubation at their respective temperature. Measurement of genome copies and infectious virus from apical washes collected over time (24 to 72 h) demonstrated that extracellular viral titers were significantly decreased at both 39 and 40°C relative to 37°C (Fig 3A and 3B). Importantly, SARS-CoV-2 infection at 37°C for 24 h prior to temperature upshift also led to reduced genome copy titers in apical washes (Fig 3C; p = 0.0056). These data indicate that temperature elevation, either prior to or following SARS-CoV-2 infection, is inhibitory to viral replication and/or release of viral particles. These data identify febrile temperature ranges associated with low (39°C) to moderate (40°C) grade fever to restrict SARS-CoV-2 propagation in respiratory epithelia.

[Figure omitted. See PDF.]

Fig 3. Elevated temperature restricts SARS-CoV-2 replication in respiratory airway cultures.

Ciliated respiratory cultures differentiated from HBEp cells isolated from donor 1 (male Caucasian aged 63 years) were incubated at 37, 39, or 40°C for 24 h prior to SARS-CoV-2 (SCV2; MOI 0.05, 104 PFU/Tissue) infection and incubation at the indicated temperatures. Apical washes were collected over time (as indicated). (A) Genome copies per ml of SARS-CoV-2 in apical washes were determined RT-qPCR. N ≥ 11 tissues per condition derived from a minimum of 3 independent biological experiments. Black line, median; whisker, 95% confidence interval; all data points shown; p-values shown, one-way ANOVA Kruskal–Wallis test. (B) TCID50 assay measuring infectious viral load in apical washes harvested from SARS-CoV-2–infected respiratory cultures. N ≥ 11 tissues per condition derived from a minimum of 3 independent biological experiments. Black line, median; whisker, 95% confidence interval; all data points shown; p-values shown, one-way ANOVA Kruskal–Wallis test. (C) Ciliated respiratory cultures were infected with SARS-CoV-2 (MOI 0.05, 104 PFU/Tissue) at 37°C and incubated for 24 h prior to continued incubation at 37 or temperature upshift to 40°C. Genome copies per ml of SARS-CoV-2 in apical washes collected at 72 h were determined RT-qPCR. N = 9 tissues per condition derived from a minimum of 3 independent biological experiments. Black line, median; whisker, 95% confidence interval; all data points shown; p-value shown, Mann–Whitney U test. Raw values presented in S9 Data. HBEp, human bronchiolar epithelial; PFU, plaque-forming unit; RT-qPCR, reverse transcription quantitative PCR; SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2.

https://doi.org/10.1371/journal.pbio.3001065.g003

We next investigated whether this temperature-dependent restriction was donor dependent. Respiratory cultures were incubated at 37 or 40°C for 24 h prior to SARS-CoV-2 infection and continued incubation at their respective temperature. Measurement of genome copy titers in apical washes demonstrated viral titers were significantly reduced at 40°C in a donor-independent manner (Fig 4A and 4B). Notably, donor-dependent profiles of restriction were observed (Fig 4B), with respiratory tissue derived from donor 1 being the most refractory to SARS-CoV-2 replication at 40°C. We next examined if the decreased yield of infectious virus observed at elevated temperature occurred due to viral genetic mutation. Amplicon sequencing of vRNA isolated from donor 1 apical washes at 72 h identified a total of 7 unique single-nucleotide polymorphisms (SNPs) relative to input sequence with no indels being identified (Table 1 and S2 Data). SNPs were only identified in genomic sequences isolated from apical washes incubated at 37°C (Table 1). Thus, the temperature-dependent restriction observed in SARS-CoV-2 replication at 40°C is not a consequence of viral mutation. As the generation of SNPs requires vRNA replication, the absence of genomic SNPs at 40°C may reflect a temperature-dependent block in vRNA replication. To investigate this, we quantified the intracellular levels of vRNA from within infected tissues. RT-qPCR analysis using 2 independent primer-probe sets (N and ORF1a) demonstrated significantly lower levels of intracellular vRNA in tissues incubated at 40°C relative to control samples (Fig 4C). ISH staining of tissue sections confirmed abundant levels of intraepithelial vRNA accumulation within all infected tissues incubated at 37°C, with little to no staining observed at 40°C (Fig 4D). Together, these data indicate that respiratory tissue remains permissive to SARS-CoV-2 infection (viral entry) but refractory to vRNA replication to levels sufficient for robust ISH detection.

[Figure omitted. See PDF.]

Fig 4. Elevated temperature restricts SARS-CoV-2 replication in respiratory airway cultures in a donor-independent manner.

https://doi.org/10.1371/journal.pbio.3001065.g004

[Figure omitted. See PDF.]

Table 1. SARS-CoV-2 sequence summary.

SARS-CoV-2 RNA isolated from apical washes collected from donor 1 infected tissues incubated at 37 or 40°C at 72 h were subjected to amplicon sequencing. Reference (ref) sequences (Wuhan, MN908947.3; England 2, EPI_ISL_407073) were used for consensus alignment and SNP identification. Unique SNPs (not detected in input and >5% base count frequency) and frequency (%base counts) shown. Full SNP analysis presented in S2 Data.

https://doi.org/10.1371/journal.pbio.3001065.t001

Elevated temperature restricts the onset of SARS-CoV-2 transcription independently of the IFN-mediated innate immune defenses

As we could detect significant levels of intracellular SARS-CoV-2 RNA within infected tissues at elevated temperature (Fig 4C), the temperature-dependent restriction observed could reflect enhanced levels of immune activation following pattern recognition receptor (PRR) detection of viral PAMPs [37,38]. To investigate this, we compared the relative transcriptome profiles from mock-treated or SARS-CoV-2–infected tissues at 37 or 40°C derived from 3 independent donors (Figs 5 and S1–S4, S1 and S3–S5 Data). Principle component analysis (PCA) demonstrated host transcriptome signatures to be donor specific (S1 Fig), with donor 3 (Caucasian female aged 16 years) demonstrating the most divergence from either of the other 2 donors (donor 1, Caucasian male aged 63 years; donor 2, Hispanic male aged 62 years). Analysis of pooled data identified distinct profiles of DEG (Q < 0.05) expression for each paired experimental condition analyzed (Fig 5A and 5B and S6 Data), with clusters of unique and shared DEG expression (Fig 5B, purple lines). Significant differences in lncRNA expression were also identified between each paired condition (S5 Fig and S7 Data). Pathway analysis identified distinct profiles of DEG enrichment, with infection at 37°C inducing an elevated immune signature (immune system process, GO:0002376) relative to all other conditions (Fig 5C, arrow). Notably, SARS-CoV-2 infection at 40°C failed to register a significant immune response above that observed for mock treatment (Fig 5C [arrow], gray box; p > 0.05). Collectively, these data identify tissue temperature to play an important role in the differential regulation of transcriptional host responses to SARS-CoV-2 infection, with enriched levels of immune gene expression observed during infection at 37°C. Investigating further, we compared the transcriptome profiles from SARS-CoV-2–infected tissues at 37 and 40°C. Out of the 1,934 DEGs (Q < 0.05) identified (Fig 6A), pathway analysis identified immune system and cytokine-related pathways to be suppressed during SARS-CoV-2 infection at 40 relative to 37°C (Fig 6B and 6C and S8 Data). DEG analysis identified IFNB1 (IFNβ) and IFNL1-3 (IFNλ) transcript levels to only be up-regulated above 1 CPM during SARS-CoV-2 infection at 37°C (Fig 6D and S8 Data). While a small number of IFN-stimulated genes (ISGs) were identified to be up-regulated within infected tissues at 40°C relative to mock treatment (9 ISGs in total; S3 Fig), no net increase in immune DEG expression was observed between mock or SARS-CoV-2–infected tissues at 40°C relative to mock treatment at 37°C (Fig 6C). Together, these data suggest that SARS-CoV-2 infection of respiratory epithelia at 40°C fails to induce a robust IFN response. To corroborate this result, we performed RT-qPCR on RNA extracted from mock or SARS-CoV-2–infected tissues derived from multiple donors. The induction of IFNB1 and IFNL1 expression was readily observed in SARS-CoV-2–infected tissues at 37°C, with no significant change observed during infection at 40°C relative to mock-treated samples (Figs 6E and S6). Similarly, significant levels of ISG induction (Mx1, ISG15, and IFIT2) were only observed in SARS-CoV-2–infected tissues at 37°C (Figs 6F and S6). Notably, the levels of ISG induction at 37°C were relatively weak (3- to 6-fold relative to mock treatment; Fig 6F), indicative of low levels of innate immune stimulation at 72 h postinfection [20,39]. Infection of respiratory cultures treated with the Janus kinase (JAK) inhibitor Ruxolitinib (Ruxo; 5 μM), which blocks IFN-mediated JAK–STAT signaling during virus infection [40,41], demonstrated no significant change in SARS-CoV-2 apical titers relative to carrier control (DMSO) treatment (Fig 6G). Thus, in contrast to temperature elevation (Figs 4 and 5), these data demonstrate IFN-mediated innate immune defenses to have a negligible effect on viral yield up to 72 h postinfection in our model system. We next examined the influence of temperature on SARS-CoV-2 replication in VeroE6 cells, a Chlorocebus sabaeus (African green monkey) cell line known to be permissive to SARS-CoV-2 infection but defective in type-I IFN-mediated host defenses [26,42,43]. Temperature elevation to 40°C restricted SARS-CoV-2 replication, with or without preincubation at elevated temperature, relative to infection at 37°C (S7 Fig). Collectively, these data demonstrate the temperature-dependent restriction of SARS-CoV-2 replication to occur independently of IFN-mediated antiviral immune defenses.

[Figure omitted. See PDF.]

Fig 5. Respiratory airway cultures induce distinct transcriptional host responses to SARS-CoV-2 infection at elevated temperature.

https://doi.org/10.1371/journal.pbio.3001065.g005

[Figure omitted. See PDF.]

Fig 6. Elevated temperature restricts SARS-CoV-2 replication in respiratory airway cultures independently of IFN-mediated antiviral immune defenses.

Ciliated respiratory cultures differentiated from primary HBEp cells isolated from 3 independent donors (donor 1, male Caucasian aged 63 years; donor 2, Hispanic male aged 62 years; donor 3, Caucasian female aged 16 years; all nonsmokers) were incubated at 37 or 40°C for 24 h prior to mock (media) or SARS-CoV-2 (SCV2; MOI 0.05, 104 PFU/Tissue) infection. Tissues were incubated at their respective temperatures for 72 h prior to RNA extraction and RNA-Seq or RT-qPCR. (A) Scatter plots showing high confidence (Q < 0.05) DEG transcripts identified between SARS-CoV-2–infected respiratory cultures at 37 or 40°C; up-regulated DEGs, red circles; down-regulated DEGs, blue circles. (B) Reactome pathway analysis of mapped down-regulated DEGs. Top 25 down-regulated (FDR < 0.05) pathways shown (blue bars; plotted as −log10 FDR). Dotted line, threshold of significance (−log10 FDR of 0.05). (C) Expression profile (log2 CPM) of down-regulated immune system DEGs (R-HSA168256; arrow in B) relative to expression levels in mock tissue at 37 or 40°C. Black line, median; dotted lines, fifth and 95th percentile range; p-values shown, one-way ANOVA (top), paired two-tailed t test (bottom). (D) Expression levels (log2 CPM) of immune system DEGs (identified in B, black arrow) relative to mock at 37 or 40°C. Every second gene labeled. Red arrows highlight genes chosen for RT-qPCR validation. (A to D) RNA-Seq data derived from RNA isolated from 3 donors (donors 1 to 3) per sample condition derived from 3 independent experiments per donor. (E/F) RT-qPCR (ΔΔCT) quantitation of IFN (IFNB1 and IFNL1) or ISG (Mx1, ISG15, and IFIT2) transcript levels within mock or SARS-CoV-2–infected respiratory cultures derived from 3 independent donors (transcript values normalized to GAPDH per sample condition). N ≥ 3 tissues per donor per condition derived from 3 independent biological experiments. Means and SD shown; all data points shown; p-values shown, one-way ANOVA. Individual ΔCT values per donor are shown in S6 Fig for donor comparison. (G) Respiratory cultures differentiated from donor 1 HBEp cells were pretreated with Ruxo (5 μM) or carrier control (DMSO) for 16 h prior to SARS-CoV-2 infection (MOI 0.05, 104 PFU/Tissue) and continued incubation at 37°C in the presence of inhibitor or carrier control. Apical washes were collected at the indicated times (h) and genome copies per ml determined RT-qPCR. N = 3 tissues per sample condition derived from 3 independent biological experiments. Means and SD shown; p-values shown, Mann–Whitney U test. Raw values presented in S8 and S9 Data. CPM, counts per million; DEG, differentially expressed gene; FDR, false discovery rate; HBEp, human bronchiolar epithelial; IFN, interferon; ISG, IFN-stimulated gene; PFU, plaque-forming unit; RT-qPCR, reverse transcription quantitative PCR; Ruxo, Ruxolitinib; SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2.

https://doi.org/10.1371/journal.pbio.3001065.g006

We next investigated at which stage in the replication cycle SARS-CoV-2 became restricted at elevated temperature. RNA-Seq analysis of SARS-CoV-2–infected tissues demonstrated no difference in the total number of mapped reads (MRs; human + viral) between 37 and 40°C (Fig 7A, p = 0.5690; S9 Data). However, a significant decrease in the number of viral MRs spanning multiple SARS-CoV-2 ORFs was detected (Fig 7B–7D). Importantly, a substantial decrease in the expression of viral subgenomic RNAs (sgRNAs) was also detected at elevated temperature (Fig 7E, p = 0.0018). As the expression of SARS-CoV-2 sgRNAs requires the onset of viral transcription [44], these data identify tissue temperature to play an important role in the regulation of viral transcription. Correspondingly, indirect immunofluorescence staining of tissue sections demonstrated significantly fewer SARS-CoV-2 nucleocapsid (N) positive foci in tissues infected at 40 relative to 37°C, confirming reduced levels of viral gene expression (Fig 7F and 7G, p < 0.0001). Costaining of tissue sections for Mx1 expression, a well-established ISG product (Fig 6D and 6F), demonstrated Mx1 expression to be localized in proximity to SARS-CoV-2 infectious foci at 37°C, with little to no staining observed within infected respiratory epithelia at 40°C (Fig 7F). We conclude that infection of respiratory tissue at elevated temperature restricts SARS-CoV-2 replication through a mechanism that inhibits viral transcription independently of IFN-mediated antiviral immune defenses previously shown to restrict SARS-CoV-2 replication [37,38].

[Figure omitted. See PDF.]

Fig 7. Elevated temperature restricts SARS-CoV-2 transcription in respiratory airway cultures.

https://doi.org/10.1371/journal.pbio.3001065.g007

Discussion

A defining symptom of COVID-19 is the onset of fever with a febrile temperature range of 38 to 41°C [4–8]. However, the effect of elevated temperature on SARS-CoV-2 tissue tropism and replication has remained to be determined. Here, we identify a temperature-sensitive phenotype in SARS-CoV-2 replication in respiratory epithelia that occurs independently of the robust induction of IFN-mediated innate immune defenses known to restrict SARS-CoV-2 replication [37,38]. The differentiation of stratified respiratory epithelium has proven to be a valuable research tool to investigate the cellular tropism, replication kinetics, and immune regulation of SARS-CoV-2, as these tissue models mimic many aspects of infection observed in animal models and COVID-19 patients [16–24,29,45–49]. While the use of such 3D models represents an important advancement over traditional two-dimensional cell culture systems, the absence of circulating immune cells (e.g., macrophages, natural killer cells, dendritic cells, and neutrophils), which modulate fever and proinflammatory immune responses to infection is an important limiting factor [9,10,50]. Thus, we limit our conclusions to the effect of elevated temperature on SARS-CoV-2 infection and replication within respiratory epithelium.

Consistent with previous reports [16,17,19–21,30,38,49], SARS-CoV-2 infection of respiratory epithelium at 37°C induced a proinflammatory immune signature (S2D Fig; IFNβ and IFNλ1–3). While SARS-CoV-2 is known to induce a comparatively weak immune response relative to other respiratory viruses (e.g., IAV) [20,39], we observed localized areas of epithelial infection (Fig 1A, arrows) that were coincident with elevated levels of ISG expression by 72 h postinfection (Fig 7F, Mx1). These findings are consistent with reports that have shown localized areas of infection in the respiratory airway of infected ferrets and surface epithelium of organoid cultures [29–31,49]. Thus, the proportion of epithelial infection and/or rate of intraepithelial spread may account for the differences observed in immune signature reported between different respiratory pathogens [20,51]. As we observed localized areas of ISG induction, our data support a model of infection where cytokine (e.g., IFN) receptor-binding consumption kinetics may be a rate limiting factor in the protection of respiratory epithelia to SARS-CoV-2 infection [52,53]. Such observations warrant further investigation to determine whether these localized areas of epithelial immune protection are induced by specific cytokines (e.g., IFNβ or IFNλ1–3) and/or correlate with COVID-19 disease progression.

We demonstrate 3D respiratory epithelial cultures under ALI to induce a robust heat stress response upon temperature elevation without loss of epithelium integrity, disruption to general cellular transcription, or damage-associated molecular pattern (DAMP) activation of innate immune defenses (Figs 2 and S4). Temperature elevation alone was not sufficient to block SARS-CoV-2 entry (Figs 3, 4, and S7), but refractory to SARS-CoV-2 transcription leading to reduced levels of vRNA accumulation and apical shedding of infectious virus (Figs 3, 4, and 7). Thus, we identify febrile temperatures associated with COVID-19 to play an important role in limiting the epithelial replication of SARS-CoV-2 infection. Importantly, SARS-CoV-2 restriction occurred independently of the robust induction of type-I (IFNβ) or type-III (IFNλ) IFN-mediated immune defenses (Figs 5, 6, and S7), despite infected tissues having abundant levels of intracellular vRNA (Fig 4C). We posit that the lack of immune induction observed at 40°C is likely a consequence of diminished levels of vRNA replication, which has been shown to play an important role in the production of PAMPs (e.g., dsRNA intermediates) required for the activation of innate immune defenses to other respiratory pathogens [54–57]. However, PRR detection of PAMPs plays an important role in the induction of a proinflammatory cytokines (including IL-6) required to mount a fever response to microbial challenge in vivo [9,10]. Thus, additional animal studies will be required to determine whether the temperature-dependent restriction of SARS-CoV-2 occurs independently of an IFN response in vivo. While speculative, we posit that low to moderate grade fever may confer protection to respiratory tissue within SARS-CoV-2–infected individuals as a component of a homeostatically controlled non-hyperinflammatory immune response to infection.

While we identify tissue temperature to play an important role in the regulation of SARS-CoV-2 transcription and replication in respiratory epithelia (Fig 7D and 7E), the precise molecular mechanism(s) of restriction remain to be determined. The heat stress response has been shown to play both a positive (proviral) and negative (antiviral) role during virus replication [58–60]. Our transcriptomic analysis identified respiratory epithelial host responses to be differentially regulated in response to both temperature elevation and SARS-CoV-2 infection (Figs 5, 6, and S1–S5). Thus, SARS-CoV-2 infection of respiratory epithelia at elevated temperature elicits a distinct host response that may contribute to the restriction of SARS-CoV-2 replication at multiple stages of infection. For example, we identify lncRNAs to be differentially expressed in response to both SARS-CoV-2 infection and temperature (S5 Fig), which are known to influence the outcome of virus infection independently of IFN-mediated immune defenses [61]. Thus, multiple gene products and/or pathways may contribute to the sequential or accumulative restriction of SARS-CoV-2 replication within respiratory epithelia at elevated temperature.

We identify temperature elevation to correlate with lower levels of SARS-CoV-2 sgRNA expression (Fig 7E), identifying a temperature-dependent block in SARS-CoV-2 transcription. We hypothesize that this restriction may relate to an inhibition in SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) activity and/or binding affinity to vRNA, as the genomic replication activity of the IAV RdRp polymerase is known to be restricted at elevated temperature (41°C) [62]. Notably, avian IAV cold adaptation of RdRp is a host determinant of IAV zoonosis, which naturally replicates in the intestinal tract of birds at higher temperatures [63]. As circulating strains of SARS-CoV-2 are susceptible to nonsynonymous mutation within RdRp coding sequences (http://cov-glue.cvr.gla.ac.uk/#/home), molecular studies are warranted to determine if such amino acid substitutions might influence the thermal restriction of SARS-CoV-2 replication. While SARS-CoV-2 has been reported to replicate to higher viral titers at temperatures associated with the upper respiratory airway (33°C) relative to core body temperature (37°C) [16], the replication of coronaviruses at febrile body temperatures has remained poorly characterized. While temperature-sensitive (ts) coronavirus mutants have been identified for feline infectious peritonitis virus (tsFIPV) and murine hepatitis virus (MHV tsNC11), their parental wild-type derivatives were able to replicate at temperatures ≥39°C [64,65]. The replication defect for tsFIPV was attributed to a block in viral maturation that limited its replication in cats to the upper respiratory airway [64], whereas the growth defect in tsNC11 was attributed to a coding mutation within the macrodomain and papain-like protease 2 domain of the nonstructural protein 3 [65]. Thus, we present the first evidence demonstrating a circulating strain of coronavirus to be sensitive to temperature thermoregulation. Genomic analysis of SARS-CoV-2 isolated from apical washes failed to identify any unique SNPs from infected tissues at 40°C relative to input sequence. These data support the temperature-dependent block in SARS-CoV-2 propagation to occur prior to vRNA replication, a prerequisite requirement for genomic mutation. As such, heat adaptation gain of function experiments through serial passage of SARS-CoV-2 at elevated temperature (≥39°C) may shed light on whether the restriction observed is related to viral and/or cellular host factors. Importantly, however, appropriate levels of biosafety (ethical and genetic modification) should be considered prior to such experimentation.

In summary, we identify an important role for tissue temperature in the cellular restriction of SARS-CoV-2 during infection of respiratory epithelia that occurs independently of the robust induction of canonical IFN-mediated antiviral immune defenses known to restrict SARS-CoV-2. We demonstrate tissue temperature to significantly influence the differential regulation of epithelial host responses to SARS-CoV-2 infection and the progression of viral transcription upon tissue infection. Future investigation is warranted to determine the precise mechanism(s) of restriction, as this may uncover novel avenues for therapeutic intervention in the treatment of COVID-19.

Materials and methods

Virus

SARS-CoV-2 (SCV2) strain BetaCoV/England/02/2020/EPI_ISL_407073 was isolated from a COVID-19 patient (Public Health England). The virus was passaged 3 times in Vero E6 cells and genotype sequence confirmed by Illumina sequencing. All experiments were performed in a Biosafety level 3-laboratory at the MRC-University of Glasgow Centre for Virus Research (SAPO/223/2017/1a).

SARS-CoV-2 replication kinetics in VeroE6 cells

Vero E6 cells (a gift from Michelle Bouloy, Institute Pasteur, France) were propagated in Dulbecco’s Modified Eagle Medium (DMEM GlutaMAX; ThermoFisher, 31966–021) supplemented with 10% fetal calf serum (FCS; ThermoFisher, 10499044 [Lot Number 08G8293K]) at 37°C with 5% CO2. Cells were seeded the day before infection at a density of 2 × 105 cells per well in a 12-well plate and infected with 104 PFU of SARS-CoV-2 diluted in serum-free DMEM for 120 min with occasional rocking. The inoculum was removed and replaced by DMEM + 10% FCS. Supernatant was collected every 24 h and snap frozen prior to analysis. For temperature inhibition experiments, cells were either preincubated at 37 or 40°C for 24 h prior to infection with continued incubation at their respective temperatures or incubated at the indicated temperatures after infection at 37°C after the addition of overlay medium (as indicated).

Differentiation and infection of respiratory epithelium

Three independent donors of primary HBEp cells were used for this study. Donor 1, 63-year-old Caucasian male (EpithelixSarl, 02AB077201F2); donor 2, 62-year-old Hispanic male (EpithelixSarl, AB079301); donor 3, 16-year-old Caucasian female (Lonza, CC-2540S). All donors were healthy and nonsmokers. For expansion and amplification, cells were propagated at 37°C with 5% CO2 in hAEC culture medium (Epithelix, EP09AM). All differentiation experiments were performed on cells that had been passaged a total of 4 times. Cells were seeded onto transwells (Falcon, 734–0036) at a cell density of 3 × 104 cells per transwell and grown to confluency. Cells were differentiated under ALI in PneumaCult-ALI medium containing hydrocortisone and heparin supplements (STEMCELL Technologies, 05001) as specified by the manufacturer’s guidelines. Cells were differentiated for a minimum of 35 days, with media changed every 48 h and washed twice weekly in serum-free DMEM after day 20 post-ALI to remove mucus build up. ALI cultures were apically washed twice before infection, 24 h prior to infection and immediately preceding infection. Cultures were inoculated apically with 100 μl of serum-free DMEM containing 104 PFU (MOI of 0.05) of SARS-CoV-2 (based on Vero E6 titres) for 120 min. The inoculum was removed, and the apical surface washed once in serum-free DMEM (0 h time point). Unless stated otherwise, ALI cultures were incubated at 37°C with 5% CO2. For temperature inhibition experiments, ALI cultures were preincubated at 37 or 40°C for 24 h prior to infection with continued incubation at their respective temperatures. For temperature upshift experiments, ALI cultures were infected at 37°C and incubated for 24 h prior to temperature upshift to 40°C or continued incubation at 37°C. For Ruxolitinib (Ruxo) experiments, DMSO (carrier control) or Ruxo (Selleckchem, S1378) was added to basal medium at a final concentration of 5 μM for 16 h prior to infection at 37°C. Cultures were incubated at 37°C and basal media changed daily containing fresh inhibitor. Apical washes were collected every 24 h in 200 μl of serum-free DMEM incubated for 30 min to collect infectious virus. Samples were stored at −80°C until required. Tissues were fixed at the indicated time points in 8% formaldehyde (Fisher Scientific, F/1501/PB17) for 16 to 24 h prior to paraffin embedding and processing. TEER measurements were conducted using an ERS-2 Voltohmmeter (Millicell; MERS00002). PBS was added to apical and basal chambers of ALI cultures or empty transwells and allowed to stand at room temperature (RT) for 15 min. Readings were taken prior to temperature shift (day 0) and every 24 h for 3 days. TEER measurements were calculated by subtracting the average of empty transwells from ALI epithelial readings and reported as Ωcm2.

TCID50

Vero E6 cells (subclone 6F5, MESO) were seeded the day before infection into 96-well plates at a density of 1.1 × 104 cells/well in DMEM + 10% FCS at 37°C with 5% CO2. Virus supernatants were serially diluted 1:10 in DMEM + 2% FCS in a total volume of 100 μl. Cells were infected in triplicate and incubated for 3 days at 37°C with 5% CO2. Plates were incubated for 3 days at 37°C in 5% CO2. Cells were fixed in 8% formaldehyde (Fisher Scientific, F/1501/PB17) for ≥2.5 h. Cytopathic effect was scored by staining with 0.1% Coomassie Brilliant Blue (BioRad, 1610406) in 45% methanol and 10% glacial acetic acid. TCID50 was calculated according to Reed–Muench and Spearman–Kärber as described [66].

Plaque assay

Vero E6 cells (subclone 6F5, MESO) were seeded at a density of 2.5 × 105 cells/well in 12-well plates the day before infection. Ten-fold serial dilutions of virus were prepared in DMEM + 2% FCS in a final volume of 500 μl. Cells were infected for 60 min with occasional rocking at 37°C in 5% CO2. The inoculum was removed and 1 ml of overlay containing 0.6% Avicel (Avicel microcrystalline cellulose, RC-591) made up in DMEM containing 2% FCS was added to each well. Cells were incubated for 3 days at 37°C in 5% CO2 prior to fixation in 8% formaldehyde and Coomassie brilliant blue staining (as described above). Plaques were counted manually and the PFU per ml calculated.

Reverse transcription quantitative PCR (RT-qPCR)

RNA was extracted from TRIzol treated apical wash samples using a QIAGEN RNeasy Mini kit (Qiagen, 74106) or in a 96-well plate format using the KingFisher Flex purification system (Thermo Scientific, 5400610) and RNAdvance kit (Beckman Coulter, A35604) following manufacturer’s protocols. SARS-CoV-2 RNA was quantified using a NEB Luna Universal Probe One-Step RT-qPCR Kit (New England Biolabs, E3006) and 2019-nCoV CDC N1 primers and probes (IDT, 10006713) or against the Orf1a open reading frame (Thermofisher, 4448508. Forward 5′- GACATAGAAGTTACTGGCGATAG-3′. Probe: VIC- ACCCCGTGACCTTGGTGCTTGT-MGB/NFQ. Reverse: 5′- TTAATATGACGCGCACTACAG-3′). Genome copy numbers were quantified using a standard curve generated from serial dilutions of an RNA standard. The RNA standard was calibrated using a plasmid (2019-nCoV_N; IDT, 10006625) that was quantified using droplet digital PCR. Values were normalized to copies/ml for apical washes or copies/tissue for cell-associated RNA.

IFN-β/λ and ISG expression analysis

For IFN-β and IFN-λ quantification, cDNA was prepared from extracted RNA using Maxima First Strand cDNA Synthesis kit for RT-qPCR (Thermofisher, K1672).

qPCR reactions were performed using TakyonTM Low Rox Probe MasterMix dTTP Blue (Eurogentec, B0701) with primers and probe detecting GAPDH (Forward: 5′- GAAGGTGAAGGTCGGAGT -3′, Reverse:5′- GAAGATGGTGATGGGATTTC-3′, Probe: FAM- CAAGCTTCCCGTTCTCAGCC -TAMRA). IFN-β and IFN-λ were analyzed with Applied Biosystems, assay identification numbers Hs01077958_s1 and Hs00601677_g1, following manufacturer’s instructions. For ISG analysis, cDNA was generated from extracted RNA using a high-capacity cDNA reverse transcription kit (Thermofisher, 4368814). qPCR reactions were performed using a Brilliant III Ultra-fast qPCR master mix (Agilent: 600884) to quantify GAPDH, Mx1, ISG15, and IFIT2/ISG54 expression using TaqMan gene-specific primer (FAM/MGB) probe mixes (Life Technologies): assay ID Mx1 (HS00895608_m1), ISG15 (Hs01921425_s1), IFIT2 (Hs01922738_s1), and GAPDH (4333764F).

In situ hybridization and immunostaining of respiratory tissue sections

Paraffin-embedded tissues were sectioned using a microtome (approximately 2 to 3 μm thick) and mounted on glass slides. Tissue sections were stained with H&E. For immunofluorescence staining, tissue sections underwent antigen retrieval by citrate (pH 6) or EDTA (pH 8) pressure cooking (as indicated below). RNAscope was used for the detection of SARS-CoV-2 RNA using a Spike-specific probe set (Advanced Cell Diagnostics, 848561 and 322372) following the manufacturer’s protocol, which included a pretreatment with boiling in target solution and proteinase K treatment. ACE2 was detected using a rabbit polyclonal anti-hACE2 antibody (Cell Signaling, 4355S; citrate antigen retrieval) and EnVision+ anti-rabbit HRP (Agilent, K4003). Slides were scanned with a bright field slide scanner (Leica, Aperio Versa 8). Mx1 was detected using a mouse monoclonal anti-Mx1 antibody ([67]; EDTA antigen retrieval). Nucleocapsid protein was detected using a sheep polyclonal anti-SARS-CoV-2 N protein antibody (University of Dundee, DA114; EDTA antigen retrieval). Secondary antibodies for detection were rabbit anti-sheep AlexaFluor 555 (abcam, ab150182) and rabbit anti-mouse AlexaFluor 488 (Sigma-Aldrich, SAB4600056). Nuclei were stained with ProLong Gold (Life Technologies, P36941). Images were collected using a Zeiss LSM 710 confocal microscope using 40× Plan-Apochromat oil objective lens (numerical aperture 1.4) using 405, 488, and 543 nm laser lines. Zen black software (Zeiss) was used for image capture and exporting images, with minimal adjustment (image rotation) in Adobe Photoshop for presentation.

RNA sequencing (RNA-Seq)

Infected or mock-treated respiratory cultures derived from a minimum of 3 independent experiments per donor per sample condition were harvested by scraping into TRIzol Reagent and transferred into tubes containing 2.8 mm ceramic beads (Stretton Scientific, P000916-LYSK0-A.0). Samples were homogenized (two 20 s pulses with a 30-s interval at RT) using a Percellys Cryolys Evolution Super Homogenizer (Bertin Instruments, P000671-CLYS2-A) at 5,500 rounds per min. The homogenized suspension was loaded into QIAshredder tubes (Qiagen, 79654) and centrifuged (12,000 × g for 2 min at RT). A volume of 0.25 ml of chloroform (VWR Life Sciences, 0757) was added to the eluate and incubated at RT for 7 min prior to centrifugation (12,000 × g, 15 min, 4°C). The aqueous phase was transferred into a fresh tube, mixed with 250 μl of 100% ethanol, and RNA isolated using RNAeasy columns (Qiagen, 74104) following the manufacturer’s protocol, which included a 15-min DNase treatment (Qiagen, 79254) treatment at RT. Eluted RNA was quantified using Qubit Fluorometer 4 (Life Technologies, Q33238), Qubit RNA HS Assay (Life Technologies, Q32855) and dsDNA HS Assay Kits (Life Technologies, Q32854) and quality controlled on a 4200 TapeStation System (Agilent Technologies, G2991A) with a High Sensitivity RNA Screen Tape assay (Agilent Technologies, 5067–5579). All samples had a RIN score of ≥8.8. Around 220 ng of total RNA was used to prepare libraries for sequencing using an Illumina TruSeq Stranded mRNA Library Prep kit (Illumina, 20020594) and SuperScript II Reverse Transcriptase (Invitrogen, 18064014) according to the manufacturer’s instructions. Dual indexed libraries were PCR amplified, then cleaned up with Agencourt AMPure XP magnetic beads (Beckman Coulter), quantified using Qubit Fluorometer 4 (Life Technologies, Q33238) and Qubit dsDNA HS Assay Kit (Life Technologies Q32854), and the size distribution assessed using a 4200 TapeStation System (Agilent Technologies, G2991A) with a High Sensitivity D1000 Screen Tape assay (Agilent Technologies, 5067–5584). Libraries were pooled in equimolar concentrations and sequenced using an Illumina NextSeq 500/550 sequencer (Illumina, FC-404-2005). At least 95% of the reads generated presented a Q score of ≥30. RNA-Seq reads were quality assessed using FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc)). Sequence adaptors removed using TrimGalore (https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/)). RNA-Seq reads were aligned to the Homo sapiens genome (GRCh38), downloaded via Ensembl using HISAT2. HISAT2 is a fast and sensitive splice aware mapper, which aligns RNA-Seq reads to mammalian-sized genomes using FM index strategy [68]. FeatureCounts [69] was used to calculate MR counts that were normalized to counts per million (CPM; unless otherwise stated). Generalized linear models (GLMs) with multifactor designs was used for DEG analysis in EdgeR [70]. A false discovery rate (FDR) p-value <0.05 (Q < 0.05) was used as a cutoff of significant differential gene expression. Sequencing reads were aligned to SARS-CoV-2 (GenBank accession MN908947.3) using BWA [71]. Periscope was used to quantify sgRNA expression levels [72]. Values were normalized to ORF1a. Only high-confidence (Q < 0.05) DEGs were used for pathway analysis in Reactome (https://reactome.org) [73,74] or differential pathway enrichment in Metascape (https://metascape.org/gp/index.html#/main/step1) [75]. In Reactome, the gene mapping tool was used as a filter to identify pathways enriched (overrepresented) for mapped DEG entities. FDR values <0.05 were considered significant for pathway enrichment. In Metascape, all DEGs were used for differential pathway analysis. Pathway p-values <0.05 were considered significant. Differentially expressed lncRNAs (Q < 0.05) were identified using the EnsemblBioMarttool (http://www.ensembl.org/biomart/martview/05285d5f063a05a82b8ba71fe18a0f18). Heat maps were plotted in GraphPad Prism (version 9.2.0). Mean CPM values of 0 were normalized to 0.01 for log2 presentation. Venn diagrams were plotted using http://bioinformatics.psb.ugent.be/webtools/Venn/.

Viral genomic amplicon sequencing

vRNA was extracted from input or extracellular released virus collected from apical washes (72 h postinfection) using TRIzol Reagent. RNA quality control was performed using Qubit (Invitrogen) and Agilent 4200 Tapestation System (Agilent). For amplicon-based sequencing, amplicons were prepared according to the protocol developed by the ARTIC Network using primer version 3 (https://artic.network/ncov-2019). DNA fragments were cleaned using AMPure beads (Beckman Coulter). Illumina sequencing libraries were prepared using 40 ng of amplicon and the KAPA HyperPrep kit (Roche, KK8502), as per manufacturer’s protocol. Samples were indexed using the NEBNext Multiplex Oligos for Illumina Dual Index Primers Set 1 (New England Biolabs, E6440S) and libraries amplified with 8 cycles of PCR. Amplified libraries were cleaned using AMPure beads (Beckman Coulter) and quantified using Qubit dsDNA HS kit (Invitrogen) and run on an Agilent 4200 Tapestation System (Agilent). Equimolar amounts of libraries were pooled, loaded at a final concentration of 11 pM, and sequenced in an Illumina MiSeq System (Illumina) using a MiSeq v2 cartridge 500 cycle kit. Sequences have been deposited in the European Nucleotide Archive (https://www.ebi.ac.uk/ena/browser/home), accession number PRJEB41332. Reads were aligned against Wuhan-Hu-1 reference strain (GenBank accession MN908947.3) using BWA [71] followed by primer trimming and consensus calling with iVar [76]. The assembled data were parsed using DiversiTools (http://josephhughes.github.io/btctools/) to determine the frequency of nucleotides at each site and observe the mutations. COV-GLUE (http://cov-glue.cvr.gla.ac.uk/#/home) was used for identifying amino acid replacement following by mutation.

Statistical analysis

The number (N) of tissues sampled per experimental condition is shown throughout and derived from a minimum of 3 independent biological experiments. GraphPad Prism (version 9.2.0) was used for statistical analysis on all data points per sample condition. Statistical tests and p-values are shown throughout. Statistically significant differences were accepted at p < 0.05.

Supporting information

S1 Fig. PCA of SARS-CoV-2–infected respiratory cultures derived from 3 independent donors.

https://doi.org/10.1371/journal.pbio.3001065.s001

(TIF)

S2 Fig. SARS-CoV-2 infection of respiratory airway cultures induces a type-I and type-III IFN response.

https://doi.org/10.1371/journal.pbio.3001065.s002

(TIF)

S3 Fig. Identification of DEGs in mock and SARS-CoV-2–infected respiratory cultures incubated at 40°C.

https://doi.org/10.1371/journal.pbio.3001065.s003

(TIF)

S4 Fig. Identification of DEGs in mock and SARS-CoV-2–infected respiratory cultures incubated at 37 and 40°C, respectively.

https://doi.org/10.1371/journal.pbio.3001065.s004

(TIF)

S5 Fig. Respiratory airway cultures induce distinct lncRNA and miRNA transcriptional host responses to SARS-CoV-2 infection at elevated temperature.

https://doi.org/10.1371/journal.pbio.3001065.s005

(TIF)

S6 Fig. Elevated temperature restricts SARS-CoV-2 replication in respiratory airway cultures in a donor- and IFN-independent manner.

Ciliated respiratory cultures differentiated from primary HBEp cells isolated from 3 independent donors (donor 1 [D1], male Caucasian aged 63 years; donor 2 [D2], Hispanic male aged 62 years; donor 3 [D3], Caucasian female aged 16 years) were incubated at 37 or 40°C for 24 h prior to mock (media only) or SARS-CoV-2 (SCV2; MOI 0.05, 104 PFU/Tissue) infection. Tissues were incubated at their respective temperatures for 72 h prior to RNA extraction and RT-qPCR. (A/B) Quantitation of IFN (IFNB1 and IFNL1) or ISG (Mx1, ISG15, and IFIT2) transcript levels within mock of SARS-CoV-2–infected respiratory cultures. N ≥ 3 tissues per sample condition derived from a minimum of 3 independent biological experiments. Means and SD shown; all data points (normalized [GAPDH] ΔCT values) shown; p-values shown, one-way ANOVA Kruskal–Wallis test. Dotted line, threshold of linear assay detection. ΔCT values presented to aid comparison between donor samples relative to that of averaged ΔΔCT values presented in Fig 6E and 6F per experimental condition. Raw values presented in S9 Data. HBEp, human bronchiolar epithelial; IFN, interferon; ISG, IFN-stimulated gene; PFU, plaque-forming unit; RT-qPCR, reverse transcription quantitative PCR; SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2.

https://doi.org/10.1371/journal.pbio.3001065.s006

(TIF)

S7 Fig. Elevated temperature restricts SARS-CoV-2 replication in Vero E6 cells.

Vero E6 cells were infected with SARS-CoV-2 (SCV2; MOI 0.05 PFU/cell) at 37°C prior to temperature elevation and incubation at 37 or 40°C (A/B) or preincubated at 37 or 40°C for 24 h prior to infection and continued incubation at their respective temperatures (C). (A) TCID50 titers of supernatants derived from SARS-CoV-2–infected Vero E6 cells incubated at 37°C over time (h). Means and SD shown. (B/C) TCID50 SARS-CoV-2 titers at 24 and 48 h postinfection. Left-hand panel; means and SD. Right-hand panel; black line, median; whisker, 95% confidence interval; all data points shown; p-values shown, Mann–Whitney U test. (A-C) N = 3 independent biological experiments performed in triplicate. Raw values presented in S9 Data. PFU, plaque-forming unit; SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2.

https://doi.org/10.1371/journal.pbio.3001065.s007

(TIF)

S1 Data. RNA-Seq data analysis for Mck-treated respiratory tissues incubated at 37 and 40°C.

Mck, mock; RNA-Seq, RNA sequencing.

https://doi.org/10.1371/journal.pbio.3001065.s008

(XLSX)

S2 Data. Identification of SCV2 SNPs from individual infected respiratory tissues incubated at 37 (37C-1 to 3) or 40°C (40C-1 to 3).

Sequences were aligned to Wuhan (MN908947.3) as a reference sequence. SCV2, SARS-CoV-2; SNP, single-nucleotide polymorphism.

https://doi.org/10.1371/journal.pbio.3001065.s009

(XLSX)

S3 Data. RNA-Seq data analysis for Mck-treated or SCV2-infected respiratory tissues at 37°C.

Mck, mock; RNA-Seq, RNA sequencing; SCV2, SARS-CoV-2.

https://doi.org/10.1371/journal.pbio.3001065.s010

(XLSX)

S4 Data. RNA-Seq data analysis for Mck-treated or SCV2-infected respiratory tissues at 40°C.

Mck, mock; RNA-Seq, RNA sequencing; SCV2, SARS-CoV-2.

https://doi.org/10.1371/journal.pbio.3001065.s011

(XLSX)

S5 Data. RNA-Seq data analysis for Mck-treated or SCV2-infected respiratory tissues at 37°C and 40°C, respectively.

Mck, mock; RNA-Seq, RNA sequencing; SCV2, SARS-CoV-2.

https://doi.org/10.1371/journal.pbio.3001065.s012

(XLSX)

S6 Data. Comparative DEG analysis of RNA-Seq data derived from Mck-treated or SCV2-infected respiratory tissues at 37 and 40°C derived from S1, S3, S4, and S5 data analysis.

DEG, differentially expressed gene; Mck, mock; RNA-Seq, RNA sequencing; SCV2, SARS-CoV-2.

https://doi.org/10.1371/journal.pbio.3001065.s013

(XLSX)

S7 Data. Comparative lncRNA and miRNA DEG analysis of RNA-Seq data derived from Mck-treated or SCV2-infected respiratory tissues at 37 and 40°C derived from S1, S3, S4, and S5 data analysis.

DEG, differentially expressed gene; lncRNA, long noncoding RNA; Mck, mock; miRNA, microRNA; RNA-Seq, RNA sequencing; SCV2, SARS-CoV-2.

https://doi.org/10.1371/journal.pbio.3001065.s014

(XLSX)

S8 Data. RNA-Seq data analysis for SCV2-infected respiratory tissues at 37 and 40°C.

RNA-Seq, RNA sequencing; SCV2, SARS-CoV-2.

https://doi.org/10.1371/journal.pbio.3001065.s015

(XLSX)

S9 Data. Underlying data used for quantitative analysis in this study.

https://doi.org/10.1371/journal.pbio.3001065.s016

(XLSX)

Acknowledgments

The authors thank Lynn Marion Stevenson, Frazer Bell, and Lynn Oxford (College of Medical, Veterinary and Life Sciences, University of Glasgow) for their exceptional efforts during the UK lockdown. We also thank Matthew Parker (University of Sheffield) for assistance in our periscope analysis.

Citation: Herder V, Dee K, Wojtus JK, Epifano I, Goldfarb D, Rozario C, et al. (2021) Elevated temperature inhibits SARS-CoV-2 replication in respiratory epithelium independently of IFN-mediated innate immune defenses. PLoS Biol 19(12): e3001065. https://doi.org/10.1371/journal.pbio.3001065

References

1. Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N Engl J Med. 2020;382(8):727–33. Epub 2020 Jan 25. pmid:31978945; PubMed Central PMCID: PMC7092803.

2. Lu R, Zhao X, Li J, Niu P, Yang B, Wu H, et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet. 2020;395(10224):565–74. Epub 2020 Feb 3. pmid:32007145; PubMed Central PMCID: PMC7159086.

3. Wu F, Zhao S, Yu B, Chen YM, Wang W, Song ZG, et al. Author Correction: A new coronavirus associated with human respiratory disease in China. Nature. 2020;580(7803):E7. Epub 2020 Apr 17. pmid:32296181; PubMed Central PMCID: PMC7608129.

4. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395(10223):497–506. Epub 2020 Jan 28. pmid:31986264; PubMed Central PMCID: PMC7159299.

5. Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. Clinical Characteristics of Coronavirus Disease 2019 in China. N Engl J Med. 2020;382(18):1708–20. Epub 2020 Feb 29. pmid:32109013; PubMed Central PMCID: PMC7092819.

6. Fu L, Wang B, Yuan T, Chen X, Ao Y, Fitzpatrick T, et al. Clinical characteristics of coronavirus disease 2019 (COVID-19) in China: A systematic review and meta-analysis. J Infect. 2020;80(6):656–65. Epub 2020 Apr 14. pmid:32283155; PubMed Central PMCID: PMC7151416.

7. Hoang A, Chorath K, Moreira A, Evans M, Burmeister-Morton F, Burmeister F, et al. COVID-19 in 7780 pediatric patients: A systematic review. EClinicalMedicine. 2020;24:100433. Epub 2020 Aug 9. pmid:32766542; PubMed Central PMCID: PMC7318942.

8. Xu T, Chen C, Zhu Z, Cui M, Chen C, Dai H, et al. Clinical features and dynamics of viral load in imported and non-imported patients with COVID-19. Int J Infect Dis. 2020;94:68–71. Epub 2020 Mar 18. pmid:32179140; PubMed Central PMCID: PMC7270709.

9. Evans SS, Repasky EA, Fisher DT. Fever and the thermal regulation of immunity: the immune system feels the heat. Nat Rev Immunol. 2015;15(6):335–49. Epub 2015 May 16. pmid:25976513; PubMed Central PMCID: PMC4786079.

10. Ogoina D. Fever, fever patterns and diseases called ’fever’—a review. J Infect Public Health. 2011;4(3):108–24. Epub 2011 Aug 17. pmid:21843857.

11. Foxman EF, Storer JA, Fitzgerald ME, Wasik BR, Hou L, Zhao H, et al. Temperature-dependent innate defense against the common cold virus limits viral replication at warm temperature in mouse airway cells. Proc Natl Acad Sci U S A. 2015;112(3):827–32. Epub 2015 Jan 7. pmid:25561542; PubMed Central PMCID: PMC4311828.

12. Earn DJ, Andrews PW, Bolker BM. Population-level effects of suppressing fever. Proc Biol Sci. 2014;281(1778):20132570. Epub 2014 Jan 24. pmid:24452021; PubMed Central PMCID: PMC3906934.

13. Schulman CI, Namias N, Doherty J, Manning RJ, Li P, Elhaddad A, et al. The effect of antipyretic therapy upon outcomes in critically ill patients: a randomized, prospective study. Surg Infect. 2005;6(4):369–75. Epub 2006 Jan 26. pmid:16433601.

14. Ryan M, Levy MM. Clinical review: fever in intensive care unit patients. Crit Care. 2003;7(3):221–5. Epub 2003 Jun 10. pmid:12793871; PubMed Central PMCID: PMC270667.

15. Cao Z, Li T, Liang L, Wang H, Wei F, Meng S, et al. Clinical characteristics of Coronavirus Disease 2019 patients in Beijing, China. PLoS ONE. 2020;15(6):e0234764. Epub 2020 Jun 20. pmid:32555674; PubMed Central PMCID: PMC7299347.

16. V’Kovski P, Gultom M, Kelly JN, Steiner S, Russeil J, Mangeat B, et al. Disparate temperature-dependent virus-host dynamics for SARS-CoV-2 and SARS-CoV in the human respiratory epithelium. PLoS Biol. 2021;19(3):e3001158. Epub 2021 Mar 30. pmid:33780434; PubMed Central PMCID: PMC8032198.

17. Ravindra NG, Alfajaro MM, Gasque V, Huston NC, Wan H, Szigeti-Buck K, et al. Single-cell longitudinal analysis of SARS-CoV-2 infection in human airway epithelium identifies target cells, alterations in gene expression, and cell state changes. PLoS Biol. 2021;19(3):e3001143. Epub 2021 Mar 18. pmid:33730024; PubMed Central PMCID: PMC8007021 "Compounds and Compositions for Treating, Ameliorating, and/or Preventing SARS-CoV-2 Infection and/or Complications Thereof." Yale University has committed to rapidly executable non-exclusive 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.

18. Zhu N, Wang W, Liu Z, Liang C, Wang W, Ye F, et al. Morphogenesis and cytopathic effect of SARS-CoV-2 infection in human airway epithelial cells. Nat Commun. 2020;11(1):3910. Epub 2020 Aug 9. pmid:32764693; PubMed Central PMCID: PMC7413383.

19. Chua RL, Lukassen S, Trump S, Hennig BP, Wendisch D, Pott F, et al. COVID-19 severity correlates with airway epithelium-immune cell interactions identified by single-cell analysis. Nat Biotechnol. 2020;38(8):970–9. Epub 2020 Jun 28. pmid:32591762.

20. Blanco-Melo D, Nilsson-Payant BE, Liu WC, Uhl S, Hoagland D, Moller R, et al. Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19. Cell. 2020;181(5):1036–45 e9. Epub 2020 May 18. pmid:32416070; PubMed Central PMCID: PMC7227586.

21. Pizzorno A, Padey B, Julien T, Trouillet-Assant S, Traversier A, Errazuriz-Cerda E, et al. Characterization and Treatment of SARS-CoV-2 in Nasal and Bronchial Human Airway Epithelia. Cell Rep Med. 2020;1(4):100059. Epub 2020 Aug 25. pmid:32835306; PubMed Central PMCID: PMC7373044.

22. Deshmukh V, Motwani R, Kumar A, Kumari C, Raza K. Histopathological observations in COVID-19: a systematic review. J Clin Pathol. 2021;74(2):76–83. Epub 2020 Aug 21. pmid:32817204.

23. Polak SB, Van Gool IC, Cohen D, von der Thusen JH, van Paassen J. A systematic review of pathological findings in COVID-19: a pathophysiological timeline and possible mechanisms of disease progression. Mod Pathol. 2020;33(11):2128–38. Epub 2020 Jun 24. pmid:32572155; PubMed Central PMCID: PMC7306927.

24. Bradley BT, Maioli H, Johnston R, Chaudhry I, Fink SL, Xu H, et al. Histopathology and ultrastructural findings of fatal COVID-19 infections in Washington State: a case series. Lancet. 2020;396(10247):320–32. Epub 2020 Jul 20. pmid:32682491; PubMed Central PMCID: PMC7365650.

25. Li W, Moore MJ, Vasilieva N, Sui J, Wong SK, Berne MA, et al. Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus. Nature. 2003;426(6965):450–4. Epub 2003 Dec 4. pmid:14647384; PubMed Central PMCID: PMC7095016.

26. Hoffmann M, Kleine-Weber H, Schroeder S, Kruger N, Herrler T, Erichsen S, et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell. 2020;181(2):271–80 e8. Epub 2020 Mar 7. pmid:32142651; PubMed Central PMCID: PMC7102627.

27. Ziegler CGK, Allon SJ, Nyquist SK, Mbano IM, Miao VN, Tzouanas CN, 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(5):1016–35 e19. Epub 2020 May 16. pmid:32413319; PubMed Central PMCID: PMC7252096.

28. Trypsteen W, Van Cleemput J, Snippenberg WV, Gerlo S, Vandekerckhove L. On the whereabouts of SARS-CoV-2 in the human body: A systematic review. PLoS Pathog. 2020;16(10):e1009037. Epub 2020 Oct 31. pmid:33125439.

29. Zaeck LM, Scheibner D, Sehl J, Muller M, Hoffmann D, Beer M, et al. Light Sheet Microscopy-Assisted 3D Analysis of SARS-CoV-2 Infection in the Respiratory Tract of the Ferret Model. Viruses. 2021;13(3). Epub 2021 Apr 4. pmid:33807059; PubMed Central PMCID: PMC8004956.

30. Triana S, Metz-Zumaran C, Ramirez C, Kee C, Doldan P, Shahraz M, et al. Single-cell analyses reveal SARS-CoV-2 interference with intrinsic immune response in the human gut. Mol Syst Biol. 2021;17(4):e10232. Epub 2021 Apr 28. pmid:33904651; PubMed Central PMCID: PMC8077299.

31. Hao S, Ning K, Kuz CA, Vorhies K, Yan Z, Qiu J. Long-Term Modeling of SARS-CoV-2 Infection of In Vitro Cultured Polarized Human Airway Epithelium. mBio 2020;11(6). Epub 2020 Nov 8. pmid:33158999.

32. Trinklein ND, Murray JI, Hartman SJ, Botstein D, Myers RM. The role of heat shock transcription factor 1 in the genome-wide regulation of the mammalian heat shock response. Mol Biol Cell. 2004;15(3):1254–61. Epub 2003 Dec 12. pmid:14668476; PubMed Central PMCID: PMC363119.

33. Morimoto RI. Cells in stress: transcriptional activation of heat shock genes. Science. 1993;259(5100):1409–10. Epub 1993 Mar 5. pmid:8451637

34. Moein S, Javanmard SH, Abedi M, Izadpanahi MH, Gheisari Y. Identification of Appropriate Housekeeping Genes for Gene Expression Analysis in Long-term Hypoxia-treated Kidney Cells. Adv Biomed Res. 2017;6:15. Epub 2017 Mar 17. pmid:28299307; PubMed Central PMCID: PMC5343613.

35. Curina A, Termanini A, Barozzi I, Prosperini E, Simonatto M, Polletti S, et al. High constitutive activity of a broad panel of housekeeping and tissue-specific cis-regulatory elements depends on a subset of ETS proteins. Genes Dev. 2017;31(4):399–412. Epub 2017 Mar 10. pmid:28275002; PubMed Central PMCID: PMC5358759.

36. Eisenberg E, Levanon EY. Human housekeeping genes, revisited. Trends Genet. 2013;29(10):569–74. Epub 2013 Jul 3. pmid:23810203.

37. Felgenhauer U, Schoen A, Gad HH, Hartmann R, Schaubmar AR, Failing K, et al. Inhibition of SARS-CoV-2 by type I and type III interferons. J Biol Chem. 2020;295(41):13958–64. Epub 2020 Jun 27. pmid:32587093; PubMed Central PMCID: PMC7549028.

38. Stanifer ML, Kee C, Cortese M, Zumaran CM, Triana S, Mukenhirn M, et al. Critical Role of Type III Interferon in Controlling SARS-CoV-2 Infection in Human Intestinal Epithelial Cells. Cell Rep. 2020;32(1):107863. Epub 2020 Jul 2. pmid:32610043; PubMed Central PMCID: PMC7303637.

39. Hadjadj J, Yatim N, Barnabei L, Corneau A, Boussier J, Smith N, et al. Impaired type I interferon activity and inflammatory responses in severe COVID-19 patients. Science. 2020;369(6504):718–24. Epub 2020 Jul 15. pmid:32661059; PubMed Central PMCID: PMC7402632.

40. Stewart CE, Randall RE, Adamson CS. Inhibitors of the interferon response enhance virus replication in vitro. PLoS ONE. 2014;9(11):e112014. Epub 2014 Nov 13. pmid:25390891; PubMed Central PMCID: PMC4229124.

41. Neubauer A, Johow J, Mack E, Burchert A, Meyn D, Kadlubiec A, et al. The janus-kinase inhibitor ruxolitinib in SARS-CoV-2 induced acute respiratory distress syndrome (ARDS). Leukemia. 2021;35(10):2917–23. Epub 2021 Aug 14. pmid:34385593; PubMed Central PMCID: PMC8358255.

42. Desmyter J, Melnick JL, Rawls WE. Defectiveness of interferon production and of rubella virus interference in a line of African green monkey kidney cells (Vero). J Virol. 1968;2(10):955–61. Epub 1968 Oct 1. PubMed Central PMCID: PMC375423. pmid:4302013

43. Mosca JD, Pitha PM. Transcriptional and posttranscriptional regulation of exogenous human beta interferon gene in simian cells defective in interferon synthesis. Mol Cell Biol. 1986;6(6):2279–83. Epub 1986 Jun 1. pmid:3785197; PubMed Central PMCID: PMC367773.

44. Kim D, Lee JY, Yang JS, Kim JW, Kim VN, Chang H. The Architecture of SARS-CoV-2 Transcriptome. Cell. 2020;181(4):914–21 e10. Epub 2020 Apr 25. pmid:32330414; PubMed Central PMCID: PMC7179501.

45. Kim YI, Kim SG, Kim SM, Kim EH, Park SJ, Yu KM, et al. Infection and Rapid Transmission of SARS-CoV-2 in Ferrets. Cell Host Microbe. 2020;27(5):704–9 e2. Epub 2020 Apr 8. pmid:32259477; PubMed Central PMCID: PMC7144857.

46. Hassan AO, Case JB, Winkler ES, Thackray LB, Kafai NM, Bailey AL, et al. A SARS-CoV-2 Infection Model in Mice Demonstrates Protection by Neutralizing Antibodies. Cell. 2020;182(3):744–53 e4. Epub 2020 Jun 20. pmid:32553273; PubMed Central PMCID: PMC7284254.

47. Sia SF, Yan LM, Chin AWH, Fung K, Choy KT, Wong AYL, et al. Pathogenesis and transmission of SARS-CoV-2 in golden hamsters. Nature. 2020;583(7818):834–8. Epub 2020 May 15. pmid:32408338; PubMed Central PMCID: PMC7394720.

48. Milewska A, Kula-Pacurar A, Wadas J, Suder A, Szczepanski A, Dabrowska A, et al. Replication of Severe Acute Respiratory Syndrome Coronavirus 2 in Human Respiratory Epithelium. J Virol. 2020;94(15). Epub 2020 May 22. pmid:32434888; PubMed Central PMCID: PMC7375387.

49. Youk J, Kim T, Evans KV, Jeong YI, Hur Y, Hong SP, et al. Three-Dimensional Human Alveolar Stem Cell Culture Models Reveal Infection Response to SARS-CoV-2. Cell Stem Cell. 2020;27(6):905–19 e10. Epub 2020 Nov 4. pmid:33142113; PubMed Central PMCID: PMC7577700.

50. Muralidharan S, Mandrekar P. Cellular stress response and innate immune signaling: integrating pathways in host defense and inflammation. J Leukoc Biol. 2013;94(6):1167–84. Epub 2013 Aug 31. pmid:23990626; PubMed Central PMCID: PMC3828604.

51. Kindler E, Jonsdottir HR, Muth D, Hamming OJ, Hartmann R, Rodriguez R, et al. Efficient replication of the novel human betacoronavirus EMC on primary human epithelium highlights its zoonotic potential. mBio 2013;4(1):e00611–12. Epub 2013 Feb 21. pmid:23422412; PubMed Central PMCID: PMC3573664.

52. Oyler-Yaniv A, Oyler-Yaniv J, Whitlock BM, Liu Z, Germain RN, Huse M, et al. A Tunable Diffusion-Consumption Mechanism of Cytokine Propagation Enables Plasticity in Cell-to-Cell Communication in the Immune System. Immunity. 2017;46(4):609–20. Epub 2017 Apr 9. pmid:28389069; PubMed Central PMCID: PMC5442880.

53. Oyler-Yaniv J, Oyler-Yaniv A, Shakiba M, Min NK, Chen YH, Cheng SY, et al. Catch and Release of Cytokines Mediated by Tumor Phosphatidylserine Converts Transient Exposure into Long-Lived Inflammation. Mol Cell. 2017;66(5):635–47 e7. Epub 2017 Jun 3. pmid:28575659; PubMed Central PMCID: PMC6611463.

54. Tapia K, Kim WK, Sun Y, Mercado-Lopez X, Dunay E, Wise M, et al. Defective viral genomes arising in vivo provide critical danger signals for the triggering of lung antiviral immunity. PLoS Pathog. 2013;9(10):e1003703. Epub 2013 Nov 10. pmid:24204261; PubMed Central PMCID: PMC3814336.

55. Killip MJ, Jackson D, Perez-Cidoncha M, Fodor E, Randall RE. Single-cell studies of IFN-beta promoter activation by wild-type and NS1-defective influenza A viruses. J Gen Virol. 2017;98(3):357–63. Epub 2016 Dec 17. pmid:27983470; PubMed Central PMCID: PMC5721924.

56. Vasilijevic J, Zamarreno N, Oliveros JC, Rodriguez-Frandsen A, Gomez G, Rodriguez G, et al. Reduced accumulation of defective viral genomes contributes to severe outcome in influenza virus infected patients. PLoS Pathog. 2017;13(10):e1006650. Epub 2017 Oct 13. pmid:29023600; PubMed Central PMCID: PMC5638565.

57. Russell AB, Elshina E, Kowalsky JR, Te Velthuis AJW, Bloom JD. Single-Cell Virus Sequencing of Influenza Infections That Trigger Innate Immunity. J Virol. 2019;93(14). Epub 2019 May 10. pmid:31068418; PubMed Central PMCID: PMC6600203.

58. Li G, Zhang J, Tong X, Liu W, Ye X. Heat shock protein 70 inhibits the activity of Influenza A virus ribonucleoprotein and blocks the replication of virus in vitro and in vivo. PLoS ONE. 2011;6(2):e16546. Epub 2011 Mar 11. pmid:21390211; PubMed Central PMCID: PMC3044721.

59. Hirayama E, Atagi H, Hiraki A, Kim J. Heat shock protein 70 is related to thermal inhibition of nuclear export of the influenza virus ribonucleoprotein complex. J Virol. 2004;78(3):1263–70. Epub 2004 Jan 15. pmid:14722281; PubMed Central PMCID: PMC321380.

60. Wan Q, Song D, Li H, He ML. Stress proteins: the biological functions in virus infection, present and challenges for target-based antiviral drug development. Signal Transduct Target Ther. 2020;5(1):125. Epub 2020 Jul 15. pmid:32661235; PubMed Central PMCID: PMC7356129.

61. Wang P, Xu J, Wang Y, Cao X. An interferon-independent lncRNA promotes viral replication by modulating cellular metabolism. Science. 2017;358(6366):1051–5. Epub 2017 Oct 28. pmid:29074580.

62. Dalton RM, Mullin AE, Amorim MJ, Medcalf E, Tiley LS, Digard P. Temperature sensitive influenza A virus genome replication results from low thermal stability of polymerase-cRNA complexes. Virol J. 2006;3:58. Epub 2006 Aug 29. pmid:16934156; PubMed Central PMCID: PMC1569369.

63. Long JS, Mistry B, Haslam SM, Barclay WS. Host and viral determinants of influenza A virus species specificity. Nat Rev Microbiol. 2019;17(2):67–81. Epub 2018 Nov 30. pmid:30487536.

64. Christianson KK, Ingersoll JD, Landon RM, Pfeiffer NE, Gerber JD. Characterization of a temperature sensitive feline infectious peritonitis coronavirus. Arch Virol. 1989;109(3–4):185–96. Epub 1989 Jan 1. pmid:2558634; PubMed Central PMCID: PMC7086852.

65. Deng X, Mettelman RC, O’Brien A, Thompson JA, E T, O’Brien TE, et al. Analysis of Coronavirus Temperature-Sensitive Mutants Reveals an Interplay between the Macrodomain and Papain-Like Protease Impacting Replication and Pathogenesis. J Virol. 2019;93(12):e021140–18. pmid:30918076

66. Hierholzer JC, Killington RA. 2—Virus isolation and quantitation. In: Mahy BWJ, Kangro HO, editors. Virology Methods Manual. London: Academic Press; 1996. p. 25–46.

67. Flohr F, Schneider-Schaulies S, Haller O, Kochs G. The central interactive region of human MxA GTPase is involved in GTPase activation and interaction with viral target structures. FEBS Lett. 1999;463(1–2):24–8. Epub 1999 Dec 22. pmid:10601631.

68. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60. Epub 2015 Mar 10. pmid:25751142; PubMed Central PMCID: PMC4655817.

69. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30. Epub 2013 Nov 15. pmid:24227677.

70. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40. Epub 2009 Nov 17. pmid:19910308; PubMed Central PMCID: PMC2796818.

71. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60. Epub 2009 Apr 20. pmid:19451168; PubMed Central PMCID: PMC2705234.

72. Parker MD, Lindsey BB, Leary S, Gaudieri S, Chopra A, Wyles M, et al. Subgenomic RNA identification in SARS-CoV-2 genomic sequencing data. Genome Res. 2021;31(4):645–58. Epub 2021 Mar 17. pmid:33722935; PubMed Central PMCID: PMC8015849.

73. Fabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2018;46(D1):D649–D55. Epub 2017 Nov 18. pmid:29145629; PubMed Central PMCID: PMC5753187.

74. Jassal B, Matthews L, Viteri G, Gong C, Lorente P, Fabregat A, et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020;48(D1):D498–D503. Epub 2019 Nov 7. pmid:31691815; PubMed Central PMCID: PMC7145712.

75. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523. Epub 2019 Apr 5. pmid:30944313; PubMed Central PMCID: PMC6447622.

76. Grubaugh ND, Gangavarapu K, Quick J, Matteson NL, De Jesus JG, Main BJ, et al. An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar. Genome Biol. 2019;20(1):8. Epub 2019 Jan 10. pmid:30621750; PubMed Central PMCID: PMC6325816.

AuthorAffiliation

About the Authors:

Vanessa Herder

Contributed equally to this work with: Vanessa Herder, Kieran Dee

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0003-4814-1382

Kieran Dee

Contributed equally to this work with: Vanessa Herder, Kieran Dee

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-3280-6768

Joanna K. Wojtus

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0003-4345-3475

Ilaria Epifano

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-2056-0510

Daniel Goldfarb

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0001-6014-387X

Christoforos Rozario

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0003-2438-5219

Quan Gu

Roles Data curation, Formal analysis, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-1201-6734

Ana Da Silva Filipe

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-9442-2903

Kyriaki Nomikou

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

Jenna Nichols

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0003-0093-7170

Ruth F. Jarrett

Roles Conceptualization, Formal analysis, Methodology, Resources

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

Andrew Stevenson

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-5325-5360

Steven McFarlane

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

Meredith E. Stewart

Roles Methodology, Resources

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0001-6275-2740

Agnieszka M. Szemiel

Roles Methodology, Resources

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0003-3085-9994

Rute M. Pinto

Roles Methodology, Resources

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0003-3291-1397

Andreu Masdefiol Garriga

Roles Data curation, Formal analysis

Affiliations MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom, University of Glasgow School of Veterinary Medicine, Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-2094-3001

Chris Davis

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-7317-3266

Jay Allan

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-6482-7125

Sheila V. Graham

Roles Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing

* E-mail: [email protected] (SVG); [email protected] (PRM); [email protected] (CB)

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-7140-8279

Pablo R. Murcia

Roles Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing

* E-mail: [email protected] (SVG); [email protected] (PRM); [email protected] (CB)

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-4352-394X

Chris Boutell

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (SVG); [email protected] (PRM); [email protected] (CB)

Affiliation: MRC-University of Glasgow Centre for Virus Research (CVR), Glasgow, Scotland United Kingdom

ORCID logo https://orcid.org/0000-0002-2970-7785

Word count: 13415

Show less

© 2021 Herder et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.

Abstract

Translate

The pandemic spread of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the etiological agent of Coronavirus Disease 2019 (COVID-19), represents an ongoing international health crisis. A key symptom of SARS-CoV-2 infection is the onset of fever, with a hyperthermic temperature range of 38 to 41°C. Fever is an evolutionarily conserved host response to microbial infection that can influence the outcome of viral pathogenicity and regulation of host innate and adaptive immune responses. However, it remains to be determined what effect elevated temperature has on SARS-CoV-2 replication. Utilizing a three-dimensional (3D) air–liquid interface (ALI) model that closely mimics the natural tissue physiology of SARS-CoV-2 infection in the respiratory airway, we identify tissue temperature to play an important role in the regulation of SARS-CoV-2 infection. Respiratory tissue incubated at 40°C remained permissive to SARS-CoV-2 entry but refractory to viral transcription, leading to significantly reduced levels of viral RNA replication and apical shedding of infectious virus. We identify tissue temperature to play an important role in the differential regulation of epithelial host responses to SARS-CoV-2 infection that impact upon multiple pathways, including intracellular immune regulation, without disruption to general transcription or epithelium integrity. We present the first evidence that febrile temperatures associated with COVID-19 inhibit SARS-CoV-2 replication in respiratory epithelia. Our data identify an important role for tissue temperature in the epithelial restriction of SARS-CoV-2 independently of canonical interferon (IFN)-mediated antiviral immune defenses.

Details

Title

Elevated temperature inhibits SARS-CoV-2 replication in respiratory epithelium independently of IFN-mediated innate immune defenses

Author

Vanessa Herder https://orcid.org/0000-0003-4814-1382; Kieran Dee https://orcid.org/0000-0002-3280-6768; Joanna K. Wojtus https://orcid.org/0000-0003-4345-3475; Ilaria Epifano https://orcid.org/0000-0002-2056-0510; Daniel Goldfarb https://orcid.org/0000-0001-6014-387X; Christoforos Rozario https://orcid.org/0000-0003-2438-5219; Quan Gu https://orcid.org/0000-0002-1201-6734; Ana Da Silva Filipe https://orcid.org/0000-0002-9442-2903; Nomikou, Kyriaki; Jenna Nichols https://orcid.org/0000-0003-0093-7170; Jarrett, Ruth F; Andrew Stevenson https://orcid.org/0000-0002-5325-5360; McFarlane, Steven; Meredith E. Stewart https://orcid.org/0000-0001-6275-2740; Agnieszka M. Szemiel https://orcid.org/0000-0003-3085-9994; Rute M. Pinto https://orcid.org/0000-0003-3291-1397; Andreu Masdefiol Garriga https://orcid.org/0000-0002-2094-3001; Chris Davis https://orcid.org/0000-0002-7317-3266; Jay Allan https://orcid.org/0000-0002-6482-7125; Sheila V. Graham https://orcid.org/0000-0002-7140-8279; Pablo R. Murcia https://orcid.org/0000-0002-4352-394X; Chris Boutell https://orcid.org/0000-0002-2970-7785

First page

e3001065

Section

Research Article

Publication year

2021

Publication date

Dec 2021

Publisher

Public Library of Science

ISSN

15449173

e-ISSN

15457885

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1371/journal.pbio.3001065

ProQuest document ID

2620105491

Elevated temperature inhibits SARS-CoV-2 replication in respiratory epithelium independently of IFN-mediated innate immune defenses

Jump to:

Full text

Abstract

Details

Suggested sources