Evolution of enhanced innate immune suppression by SARS-CoV-2 Omicron subvariants (News)

Abstract


Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) human adaptation resulted in distinct lineages with enhanced transmissibility called variants of concern (VOCs). Omicron is the first VOC to evolve distinct globally dominant subvariants. Here we compared their replication in human cell lines and primary airway cultures and measured host responses to infection. We discovered that subvariants BA.4 and BA.5 have improved their suppression of innate immunity when compared with earlier subvariants BA.1 and BA.2. Similarly, more recent subvariants (BA.2.75 and XBB lineages) also triggered reduced innate immune activation. This correlated with increased expression of viral innate antagonists Orf6 and nucleocapsid, reminiscent of VOCs Alpha to Delta. Increased Orf6 levels suppressed host innate responses to infection by decreasing IRF3 and STAT1 signalling measured by transcription factor phosphorylation and nuclear translocation. Our data suggest that convergent evolution of enhanced innate immune antagonist expression is a common pathway of human adaptation and link Omicron subvariant dominance to improved innate immune evasion.

Main


Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) Alpha, Delta and then Omicron became sequentially dominant globally, with each evolving independently from wave 1 early lineage SARS-CoV-2 virus. Sequential lineage replacement suggests evolution of highly advantageous characteristics that effectively improved transmission. Our previous work showed that Alpha1, and also VOCs Beta to Delta2, adapted by enhancing expression of specific innate immune antagonists including Orf6, N and Orf9b, to suppress the host innate immune response initiated on infection. Since the appearance of the Omicron lineage, it is Omicron subvariants that are co-circulating, or being replaced by each other, rather than new, wave 1-derived, VOCs. The selective forces driving SARS-CoV-2 evolution may therefore have switched from being predominantly adaptation-to-host to immune escape from vaccine- and infection-driven memory responses. In fact, the first dominant Omicron subvariants BA.1 and BA.2, BA.4 and BA.5 emerged with each displaying increasing levels of antibody escape, through mutation of spike, threatening vaccine efficacy and increasing hospitalizations3,4,5,6,7,8,9,10,11,12,13,14,15,16. However, like Alpha to Delta, Omicron subvariants are also accumulating mutations beyond spike17,18, suggesting that spike-independent adaptations may also be crucial for Omicron variant dominance. In this Article, we provide evidence that, similar to VOCs Alpha to Delta, Omicron variants also improve innate immune evasion through enhancement of viral protein expression, suggesting that regulation of host responses through adapting viral protein levels is a key feature of SARS-CoV-2 evolution.

Results


To understand phenotypic differences between Omicron subvariants, and the selective forces driving their evolution, we compared replication of, and host responses to, BA.1–BA.5 with Delta, the previously dominant VOC, in Calu-3 human airway epithelial cells (HAEs; Fig. 1). We equalized input dose of each variant by viral envelope (E) gene copies (quantitative reverse transcription polymerase chain reaction, RT–qPCR) as this ensures cells are exposed to equal starting amounts of viral RNA, which is the major viral PAMP activating defensive host innate immune responses1,19. Most importantly, this approach normalizes dose independently of variant-specific differences in cell tropism or entry routes (Fig. 1a and Extended Data Fig. 1a,b)20,21,22, which we and others have shown impact both titre determination and input equalization by cell-line infectivity measurements such as 50% tissue culture infectious dose (TCID50) or plaque assay (Extended Data Fig. 1c–e). Our approach is particularly relevant for comparing Omicron subvariants because Omicron spike mutations have been shown to alter tropism, increasing cathepsin-dependent endosomal entry and reducing dependence on cell surface TMPRSS2 (refs. 20,21,22,23), irrespective of virion spike cleavage efficiency (Extended Data Fig. 1f). Endosomal cathepsins or cell surface TMPRSS2 are required to cleave spike before ACE2-mediated entry24,25. Indeed, in line with previously published data20,21,22, we have found that Omicron, particularly BA.5, has enhanced entry (cathepsin dependent and E64d sensitive) in TMPRSS2-negative cells such as Hela-ACE2 compared with previous VOCs such as Delta, whereas entry into Calu-3 cells is largely TMPRSS2 dependent (camostat sensitive) (Extended Data Fig. 1a,b), resulting in striking cell type-specific differences between variant titres by TCID50 (Extended Data Fig. 1e).

Fig. 1: BA.5 displays enhanced innate immune antagonism during infection of airway epithelial cells.

figure 1

ag, Calu-3 infection with 2,000 E copies per cell of Delta (yellow, Ο), BA.1 (blue, Ο), BA.2 (blue, Δ), BA.4 (purple, O) and BA.5 (purple, Δ), n = 3: mean viral E copies at 2 h.p.i. across three independent experiments (a); viral replication over time measured by RT–qPCR for intracellular E copies per microgram RNA (b); infection levels measured by nucleocapsid expression (% N+ by flow cytometry) (c); expression of IFNBCXCL10IFIT1IFIT2RSAD2MX1MX2 and DDX58 in infected cells over time (d); IFNβ (e) and CXCL10 (f) secretion from infected Calu-3 cells measured by ELISA at 48 h.p.i.; rescue of viral replication by JAK1-inhibitor ruxolitinib in Calu-3 cells at 48 h.p.i., where relative infection levels are shown across three independent experiments determined by E copies per microgram RNA normalized to the median infection level of the untreated control (g). hk, Primary bronchial HAEs were infected with the indicated variants at 1,500 E copies per cell: viral replication measured by intracellular E copies at 72 h.p.i. (h) and viral release into apical washes over time (i), with three biological replicates shown; expression of IFNBCXCL10IFIT1IFIT2DDX58 and RSAD2 in HAEs at 72 h.p.i., with six biological replicates shown (j); intracellular viral E copies in HAEs in the presence or absence of 5 μM ruxolitinib at 72 h.p.i., with three biological replicates shown (k). For a, one-way analysis of variance (ANOVA) with Bonferroni post-test was used. n.s., not significant at P > 0.05 for all comparisons. For bh and j, one-way ANOVA and Dunnett’s post-test were used. For i, two-way ANOVA with a Bonferroni post-test was used. For k, one-tailed unpaired Student’s t-test was used. Replicate measurements from one of three independent experiments. Fold change over mock is shown. Mean ± s.e.m. or individual datapoints are shown. h.p.i., hours post infection. Source data

Infection of Calu-3 cells with 2,000 E gene copies per cell (Fig. 1) or 200 E copies per cell (Extended Data Fig. 1) gave comparable E RNA (RT–qPCR) at 2 h post infection (h.p.i.), consistent with equal input doses (Fig. 1a and Extended Data Fig. 1g). E gene measurements during infection revealed that Omicron isolates BA.1, BA.2, BA.4 and BA.5 replicated similarly, lagging behind Delta in Calu-3 cells (Fig. 1b and Extended Data Fig. 1h–l). BA.4 replicated most slowly initially but caught up with BA.1, BA.2 and BA.5 by 24 h.p.i. (Fig. 1b and Extended Data Fig. 1). Importantly, these replication differences were observed consistently across several experiments (Fig. 1 and Extended Data Figs. 1 and 2). As E gene measurement during infection captures genomic RNA (gRNA) as well as E, S and Orf3 subgenomic RNAs (sgRNAs), we compared the levels of intracellular E RNA with those of Nsp12 and Orf1a (compare Extended Data Fig. 1h,i with Fig. 1b and Extended Data Fig. 1k,l with Extended Data Fig. 1j), which are uniquely encoded within gRNA. Importantly, the ratio of E to Nsp12 was similar until 24 h.p.i. reflecting equivalent levels of E sgRNA synthesis between variants (Extended Data Fig. 1m). Quantification of released virions by measuring E and Nsp12 RNA copies in the supernatant mirrored viral replication (Extended Data Fig. 1n–q). Similar patterns of infection were also seen when quantified by intracellular nucleocapsid (N) staining (Fig. 1c and Extended Data Fig. 1r).

BA.4 and BA.5 trigger less innate immune activation than earliest Omicron subvariants

We next compared the host innate immune response to Omicron subvariant infection of Calu-3 cells. All viral stocks were prepared in human gastrointestinal Caco-2 cells as they are naturally permissive to SARS-CoV-2 replication but do not mount a strong innate response to this infection19,26. We confirmed that viral stocks prepared in Caco-2 cells (the highest viral inoculum for each variant was 2,000 E copies per cell) did not contain measurable interferon (IFN)β and negligible IFNλ1/IFNλ3 (enzyme-linked immunosorbent assay, ELISA) (Extended Data Fig. 2a,b), ensuring differences in innate immune activation in Calu-3 infections were not a result of IFN carryover in the viral stocks.

Strikingly, we found that infection of Calu-3 cells with BA.4 and BA.5 resulted in significantly less innate immune activation compared to BA.1/BA.2, evidenced by lower induction of IFNβ (IFNB) and interferon stimulated genes (ISGs) including inflammatory chemokine CXCL10 and RSAD2DDX58IFIT1 and IFIT2 (Fig. 1d and Extended Data Fig. 2c–g) and a trend towards reduced MX1 and MX2 expression (Fig. 1d). Reduced host responses to BA.4 and BA.5 infection were also evident at the level of IFNβ and CXCL10 secretion (Fig. 1e,f). Slower replication of BA.4 probably contributes in part to reduced innate immune activation during Calu-3 infection, but BA.5 replication was similar to BA.1 and BA.2 and nonetheless induced significantly less innate immune responses. Inhibition of IFN-mediated JAK/STAT signalling with ruxolitinib, evidenced by the absence of ISG induction (Extended Data Fig. 2e,f), rescued BA.1 and BA.2 infection in Calu-3 cells to a greater degree than BA.4 or BA.5 (Fig. 1g and Extended Data Fig. 2h–j), suggesting that the greater induction of IFNβ by BA.1 and BA.2 reduced their infectivity. BA.1 to BA.5 showed similar sensitivities to a range of IFN doses used to pre-treat Calu-3 cells (Extended Data Fig. 2k–m). We therefore conclude that the differences in ruxolitinib sensitivity reflect differences in IFN induction after Calu-3 infection and not differences in IFN sensitivity. Infecting Calu-3 cells with lower virus input doses (200 E copies per cell) recapitulated our observation that Delta replicated better than Omicron BA.1–BA.5 (Extended Data Fig. 1j–l), and we again saw reduced innate immune activation by BA.4 and BA.5 compared with BA.1 and BA.2 (Extended Data Fig. 2f,g). At this lower inoculum, BA.4 infectivity was also strongly rescued by ruxolitinib treatment consistent with its slower replication being due to IFN induction (Extended Data Fig. 2i).

We next compared Omicron subvariant replication and host responses in primary HAE cultures, which better recapitulate the heterogeneous polarized epithelial layer of the respiratory tract. We have previously reported that HAEs reveal differences in VOC replication that probably reflect host adaptation, which are not always apparent in highly permissive cell lines, such as Calu-3 (refs. 1,2). Concordantly, BA.5 viral replication was higher than BA.2 and BA.4 in differentiated primary bronchial HAEs at 72 h.p.i., while apical viral release over time was comparable (Fig. 1h,i). Despite BA.4 and BA.5 replicating similarly to BA.2 in HAEs, we consistently observed reduced innate activation, measured by ISG induction, after BA.4 and BA.5 infection (IFNBCXCL10IFIT1IFIT2DDX58 and RSAD2; Fig. 1j). Inhibiting IFN signalling with JAK-inhibitor ruxolitinib suppressed ISG induction (Fig. 1j) and rescued replication of BA.2 to a greater degree than BA.4 and BA.5 (Fig. 1k). Altogether, data in Fig. 1 suggest adaptation to reduce innate immune activation between the earliest (BA.1 and BA.2) and subsequent (BA.4 and BA.5) Omicron subvariants.

SARS-CoV-2, and other respiratory viruses, reportedly replicate more efficiently in nasal and tracheal epithelial cells27, in part due to reduced innate activation and IFN responsiveness at the lower temperatures of the upper airway28,29,30. To investigate whether lower temperatures reveal further Omicron subvariant adaptation, we compared replication at 32 °C in Calu-3 cells. We found BA.1 to BA.5 all replicated less well than at 37 °C (Extended Data Fig. 3a,b) whereas Delta replication was not as temperature sensitive. As expected29, innate immune activation in response to infection, or to RNA sensing agonist poly(I:C), was largely abolished at 32 °C (measured by IFNB and CXCL10 messenger RNA induction; Extended Data Fig. 3c–e). At 37 °C, we again observed lower innate activation for BA.4 and BA.5 compared with BA.1/BA.2. In HAE, lowering the temperature to 32 °C did not impact viral replication to the same extent as in Calu-3 cells (Extended Data Fig. 3f). However, we observed reduced virus output in apical washes from infected HAE cultures for all Omicron isolates (Extended Data Fig. 3g–i). Infected HAEs at 32 °C also expressed significantly less IFNB and CXCL10 (Extended Data Fig. 3j). Overall, our data suggest that Omicron does not replicate better at 32 °C in lung epithelial cells even in the absence of an innate immune response. However, it is possible that the intra-tissue temperature throughout the airways remains closer to 37 °C than the exhaled breath temperature of 32 °C suggests31.

BA.4 and BA.5 increase Orf6 expression and efficiently antagonize innate immune activation during infection

We next investigated the mechanism underlying differential innate immune activation by Omicron subvariants. IRF3 and STAT1 are key transcription factors responding to intracellular RNA sensing and IFN production, respectively, exemplified here by poly(I:C) treatment (Extended Data Fig. 4a–c). We and others have shown that SARS-CoV-2 activates transcription factors IRF3 and STAT1 downstream of RNA sensing19,32. Consistent with their reduced innate immune triggering, we found Omicron BA.4 and BA.5 infection activated significantly less IRF3 phosphorylation than BA.2 infection (Fig. 2a–c). A similar trend was observed for STAT1 serine 727 phosphorylation, which is essential for full STAT1 transcriptional activity33, but not upstream JAK1-dependent tyrosine 701 phosphorylation (Fig. 2a,d–f). Reduction of STAT1 phosphorylation correlated with reduced STAT1 nuclear translocation during BA.4 and BA.5 infection compared with BA.2, measured by high-content single-cell immunofluorescence imaging of infected nucleocapsid-positive Calu-3 cells (Fig. 2g). These data suggest that BA.4 and BA.5 more effectively prevent intracellular activation of innate sensing pathways. We previously reported that SARS-CoV-2 VOC Alpha evolved enhanced innate immune evasion by increasing expression of key innate antagonists Orf6, Orf9b and N (Extended Data Fig. 4d), which manipulate host cell innate immune pathways1. To investigate whether Omicron subvariants have also independently evolved enhanced innate immune suppression through similar mechanisms during human adaptation, we measured viral innate antagonist protein expression during infection. Strikingly, we found that BA.4, and particularly BA.5, expressed higher levels of Orf6 and N compared with BA.1 and BA.2 (Fig. 2h–l and Extended Data Fig. 4e–k), measured at 48 h.p.i. in Calu-3 cells when E RNA levels were equivalent (Fig. 2m). Unlike previous VOCs1,2, expression of innate immune antagonist Orf9b was not detected for any Omicron isolate, possibly due to Omicron subvariants encoding lineage-specific Orf9b mutations (P10S and ΔENA at positions 27–29) altering antibody binding and precluding detection by immunoblot (Fig. 2n and Extended Data Fig. 4d). Importantly, Orf9b remained readily detectable in Delta-infected cells (Fig. 2n). Upregulation of Orf6 and N expression by BA.5 was validated using a second independent isolate (Extended Data Fig. 4l–n), and was also evident in lysates from infected HAEs (Extended Data Fig. 4o). Blocking IFN signalling with ruxolitinib rescued replication of all Omicron isolates as before (Fig. 1 and Extended Data Fig. 2) and enhanced viral protein detection by immunoblot (Fig. 2h,n and Extended Data Fig. 4e). Importantly, higher levels of BA.4 and BA.5 Orf6 and N remained apparent after ruxolitinib treatment (Fig. 2h,k,l). We previously showed that enhanced levels of Orf6, N and Orf9b protein by Alpha were associated with increased levels of the corresponding sgRNAs1. By contrast, BA.5 Orf6 and N sgRNA levels (normalized to genomic Orf1a) were not enhanced, and were only slightly upregulated during BA.4 infection (Fig. 2o), particularly in comparison with Alpha (Extended Data Fig. 4p,q). No differences were observed in S and Orf3a sgRNAs, which served as controls to rule out a general enhancement of sgRNA synthesis (Fig. 2o). Although Omicron subvariants have synonymous and non-synonymous mutations in Orf6 and N, there are no mutations that distinguish BA.4 and BA.5 from BA.1 and BA.2 that provide a simple explanation for increased Orf6 or N protein levels, including in their transcriptional regulatory sequences (Figs. 1 and 2 and Extended Data Tables 1 and 2). Thus, we hypothesize that BA.4 and BA.5 have either evolved independent mechanisms to increase Orf6 and N protein levels, or that the increase is mediated by changes elsewhere in the genome, which may impact viral translation or protein stability. Further studies are required to pinpoint the adaptations regulating Orf6 and N expression levels.

Fig. 2: BA.5 efficiently expresses SARS-CoV-2 innate antagonists during airway epithelial cell infection.

figure 2

Calu-3 cells were infected with 2,000 E copies per cell of the indicated variants. a, Western blot of STAT1-pY701, STAT1-pS727, total STAT1, IRF3-pS396, total IRF3 and β-actin at 24 h.p.i. One of four independent western blots is shown. bf, Quantification of four independent western blots showing IRF3-pS396 (b), IRF3 (c), STAT1-pS727 (d), STAT1-pY701 (e) and STAT1 (f) over β-actin at 24 h.p.i. normalized to mock. g, Quantification of STAT1 nuclear translocation detected by single-cell fluorescence microscopy over time in Calu-3 cells infected with the indicated variants. Data from 1,500 cells per condition are shown. In infected cultures, translocation was determined in N+ cells. h, Western blot of Orf6, N, spike and β-actin at 48 h.p.i. in infected cells ± 5 μM ruxolitinib (Rux). Non-specific bands detected by polyclonal anti-spike primary antibody are indicated (see Extended Data Fig. 4e for mock). One of five independent western blots shown. il, Quantification of Orf6 and N expression from five independent western blots of Calu-3 cells in the absence (i, Orf6; j, N) or presence of 5 μM ruxolitinib (k, Orf6; l, N) at 48 h.p.i., normalized to spike over BA.2. m, Viral replication in cells from hn, Representative western blot of Calu-3 cells infected with Delta, BA.1, BA.2, BA.4 and BA.5 at 2,000 E copies per cell showing Orf9b, Orf6, N and β-actin expression at 48 h.p.i. + 5 μM ruxolitinib. o, sgRNA expression of Orf6, N, spike and Orf3a normalized to Orf1a gRNA in Calu-3 cells at 48 h.p.i.; nine measurements from three independent experiments shown. For bfi–m and o, one-way analysis of variance with Dunnett’s post-test was used. For g, box-and-whisker blots show 10th–90th percentile, and groups were compared at each timepoint as indicated using a Kruskal–Wallis test. Mean ± s.e.m. or individual datapoints are shown. Source data

Orf6 expression is a major determinant of enhanced innate immune antagonism by emerging VOCs

Orf6 is a multifunctional viral accessory protein that modulates expression of host and viral proteins34,35. Orf6 selectively inhibits host transcription factor nuclear transport to potently antagonize antiviral responses during infection. To probe Orf6 mechanisms, and its contribution to enhanced innate antagonism by the VOCs, we used reverse genetics to introduce two stop codons into the Orf6 coding sequence of both Alpha (Alpha ΔOrf6) and BA.5 (BA.5 ΔOrf6), which we confirmed abolished Orf6 expression during infection (Fig. 3a,b). While Alpha ΔOrf6 replicated similarly to parental wild-type (WT) virus up to 24 h.p.i. (Fig. 3c), we observed enhanced IFNB and CXCL10 expression (Fig. 3d) and protein secretion (Extended Data Fig. 5a) during Alpha ΔOrf6 infection of Calu-3 cells compared with WT virus. Moreover, increased IRF3 nuclear translocation was evident after Alpha ΔOrf6 infection at 24 h.p.i. using single-cell quantitative immunofluorescence microscopy (Fig. 3e and Extended Data Fig. 5b). This suggests an important role for Orf6 in innate immune antagonism during viral replication1,35,36,37 and is consistent with suppression of IRF3 nuclear transport in Orf6 overexpression studies35,36,38. The reduction in Alpha ΔOrf6 replication at 48 h.p.i., and N and spike protein expression at 24 h.p.i., that was rescued by ruxolitinib treatment, is also consistent with greater IFN-mediated suppression of the Orf6 deletion mutant (Fig. 3a and Extended Data Fig. 5c).

Fig. 3: Orf6 expression is a major determinant of enhanced innate immune antagonism by emerging VOCs.

figure 3

a,b, Western blot of Alpha (a) or BA.5 (b) reverse genetic (RG) virus infections in Calu-3 cells at 24 h.p.i. ± 5 μM ruxolitinib (Rux). c, Replication of RG viruses parental Alpha WT and ΔOrf6 in Calu-3 cells infected with 2,000 E copies per cell over time. d, Gene expression in cells from c over time. e, Quantification of IRF3 nuclear translocation detected by single-cell fluorescence microscopy over time. fi, HAEs were infected with 1,500 E copies per cell of the indicated variants ± 5 μM ruxolitinib. f, Viral release into apical washes over time. g,h, Apical release in HAEs infected with Alpha WT (g) or ΔOrf6 ± 5 μM ruxolitinib (h). i, Gene expression in cells from f. Three biological replicates shown. j, Replication of RG viruses BA.5 WT, ΔOrf6 and Orf6 D61L isolates in Calu-3 cells infected with 2,000 E copies per cell over time. kIFNB expression in cells from jl, Gene expression of Calu-3 cells at 24 h.p.i. m, Western blot of STAT1-pY701, STAT1-pS727, total STAT1, IRF3-pS396, total IRF3 and β-actin at 24 h.p.i. nr, Quantification of five independent western blots showing IRF3-pS396 (n), total IRF3 (o), STAT1-pS727 (p), STAT1-pY701 (q) and total STAT1 (r) over β-actin at 24 h.p.i. s,t, Quantification of IRF3 (s) and STAT1 (t) nuclear translocation detected by single-cell fluorescence microscopy at 24 h.p.i. u,v, Replication of BA.5 WT and ΔOrf6 in HAEs infected with 1,500 E copies per cell in the absence (u) or presence (v) of 5 μM ruxolitinib. For c and d, two-way analysis of variance (ANOVA) and Bonferroni post-test were used. For es and t, data from 1,500 cells per condition are shown as box-and-whisker blots indicating 10th–90th percentile. In infected cultures, translocation was determined in N+ cells. Groups were compared by Kruskal–Wallis test. For kl and nr, one-way ANOVA with Dunnett’s post-test was used. For fiu and v, unpaired two-tailed Student’s t-test was used. Replicate measurements from one of three independent experiments. Fold change over mock is shown. Mean ± s.e.m. or individual datapoints are shown. Source data

Alpha ΔOrf6 also replicated less well than WT in HAE cells (Fig. 3f–h and Extended Data Fig. 5d). IFNB and CXCL10 gene induction, normalized to GAPDH, were similar after Alpha ΔOrf6 and WT infection (Fig. 3i), despite lower E RNA levels for Alpha ΔOrf6, consistent with increased innate immune induction by the deletion virus. Importantly, Alpha ΔOrf6 was more sensitive to ruxolitinib treatment than WT, consistent with the notion that increased IFN induction caused reduced replication of Alpha ΔOrf6 (Fig. 3g,h). To address the role of Orf6 during BA.5 infection, we compared replication of a BA.5 ΔOrf6 mutant with parental BA.5 WT virus. We also generated a BA.5 mutant bearing the Orf6 D61L mutation found in BA.2 and BA.4 that has been proposed to reduce Orf6 function2,32,39 (Fig. 3b,j). Consistent with the SARS-CoV-2 Alpha ΔOrf6 results, BA.5 ΔOrf6 showed a replication defect at 48 h.p.i. compared with BA.5 WT, and triggered significantly enhanced innate immune responses evidenced by enhanced IFNB and ISG induction (Fig. 3k,l). Deletion of Orf6 in BA.5 also increased the degree of infection-induced IRF3 and STAT1 phosphorylation (Fig. 3m–r) and nuclear translocation (Fig. 3s,t). This demonstrates that Orf6 loss enhances IRF3 and STAT1 activation despite similar levels of infection in the first 24 h.p.i., confirming the important role of Orf6 in innate immune suppression and in distinguishing BA.5 from earlier Omicron subvariants. Infection of HAEs confirmed reduced viral replication of BA.5 ΔOrf6 compared with WT BA.5, while viral release remained comparable (Fig. 3u,v and Extended Data Fig. 5e). ISG expression in HAEs was similar between WT and mutant despite lower E RNA levels during BA.5 ΔOrf6 infection, suggesting greater induction of innate immunity in the absence of Orf6 in these cells (Extended Data Fig. 5f). Interestingly, introducing the C-terminal D61L mutation into BA.5 Orf6 resulted in an intermediate innate immune phenotype measured by increased induction of IFNBCXCL10 and IFIT1 expression by the mutant virus (Fig. 3l). IRF3 phosphorylation and nuclear translocation were equivalent between BA.5 WT and Orf6 D61L (Fig. 3n–s), whereas STAT1 translocation was not antagonized by Orf6 D61L (Fig. 3t), in line with reports of a partial loss of Orf6 function in the D61L mutation2,32,39. These data suggest complex adaptation of Orf6 manipulation of innate immunity during SARS-CoV-2 Omicron lineage adaptation.

Enhanced innate antagonism is a conserved feature of dominant Omicron subvariants

During the course of this study, SARS-CoV-2 has continued to evolve and produce new Omicron subvariants (Fig. 4a and Extended Data Fig. 6a). Omicron subvariants BA.2.75, XBB.1, XBB.1.5 and BQ.1.1 have acquired increased ACE2 binding and enhanced adaptive immune evasion40,41,42,43. To test whether enhanced innate immune antagonism is consistently associated with globally successful subvariants, we compared BA.2.75, XBB.1, XBB.1.5 and BQ.1.1 isolates with BA.2 and BA.5 (Fig. 4). We equalized virus dose by Nsp12 RNA copies (RT–qPCR), a measurement of gRNA, rather than E RNA copies, due to accumulation of mutations in the E gene of later Omicron subvariants, including in the region detected by our RT–qPCR assay. We found that all Omicron subvariants retained an enhanced dependence on cathepsin, here measured in A549 cells expressing ACE2 and TMPRSS2 (Extended Data Fig. 6b). BA.2.75, XBB.1 (two independent isolates) and XBB.1.5, derived from the parental BA.2 lineage41,43, replicated comparably to earlier BA.2 and BA.5 in Calu-3 and HAEs (Fig. 4b–e and Extended Data Fig. 6c–h). BQ.1.1, which has arisen from BA.5 (ref. 43), displayed some reduction of replication in epithelial cells (Fig. 4d,e and Extended Data Fig. 6e,h). Similar to BA.5, we found that all subsequent Omicron subvariants tested triggered significantly less IFNB and CXCL10 expression than BA.2 at 24 h.p.i. (Fig. 4f). All Omicron subvariants derived from BA.2 (BA.2.75, XBB.1 and XBB.1.5) showed reduced rescue by ruxolitinib treatment, as well as reduced induction of, or sensitivity to, IFN, similar to BA.5 (Fig. 4g and Extended Data Fig. 6i). Strikingly, like BA.5, enhanced innate immune evasion by these more recent subvariants was accompanied by increased Orf6 expression for the majority of isolates (Fig. 4h,i). Reduced BQ.1.1 replication in Calu-3 cells (Fig. 4d and Extended Data Fig. 6e) prevented Orf6 and N detection in the absence of ruxolitinib (Fig. 4h). Reduced innate activation by recent Omicron subvariants also correlated with reduced IRF3 phosphorylation compared with BA.2, and reduction of STAT1 serine phosphorylation was principally observed for XBB.1 and XBB.1.5 variants (Fig. 4j–l and Extended Data Fig. 6j–l). Together these data are consistent with a trend for ongoing Omicron evolution enhancing Orf6 expression as it adapts to the human population leading to reduced innate immune responses, detectable at the level of IFN and ISG expression, and at the level of transcription factor phosphorylation and nuclear translocation. This study considering Omicron variants is very reminiscent of our previous observation of enhanced expression of key innate immune antagonists Orf6, N and Orf9b in VOCs Alpha to Delta suggesting a common evolutionary trajectory to combatting human innate immunity to enhance transmission1,2.

Fig. 4: Innate immune phenotype of dominant Omicron subvariants.

figure 4

a, Global SARS-CoV-2 variant sequence counts over time (scaled per variant), extracted from CoV-Spectrum using genomic data from GISAID. bd, Calu-3 cells were infected with 2,000 Nsp12 copies per cell. Replication of Omicron subvariants compared with BA.2 (blue) and BA.5 (purple) measured by Nsp12 copies per microgram RNA is shown for BA.2.75 (yellow; Ο) (b), XBB subvariants (XBB.1: light red, Ο; XBB.1 (B): red, Δ; XBB.1.5: dark red, □) (c) and BQ.1.1 (BQ.1.1: light green, Ο; BQ.1.1 (B): dark green, Δ) (d) isolates. e, HAEs were infected with 1,500 Nsp12 copies per cell and intracellular Nsp12 copies measured at 72 h.p.i. Three biological replicates shown. fIFNB and CXCL10 expression in Calu-3 cells infected with 2,000 Nsp12 copies per cell of the indicated Omicron subvariants at 24 h.p.i. g, Viral replication of indicated variants in Calu-3 cells in the presence or absence of 5 μM ruxolitinib at 48 h.p.i. Numbers indicate fold change in replication in the presence of 5 μM ruxolitinib. h,i, Western blot of Orf6, N, spike and β-actin at 48 h.p.i. in cells from bd in the absence (h) or presence (i) of 5 μM ruxolitinib. j, Western blot of STAT1-pY701, STAT1-pS727, total STAT1, IRF3-pS396, total IRF3 and β-actin in Calu-3 cells at 48 h.p.i. k,l, Quantification of two independent western blots of IRF3-pS396 (k) and STAT1-pS727 (l) over β-actin at 24 h.p.i. For bd, variant replication was compared with BA.2 at each timepoint using a two-way analysis of variance (ANOVA) and Bonferroni post-test. Colours indicate comparator (BA.5, purple; BA.2.75, yellow; XBB.1, light red; XBB.1 (B), red; XBB.1.5, dark red; BQ.1.1, light green; BQ.1.1 (B), dark green). For eg, one-way ANOVA with Dunnett’s post-test was used to compare all variants with BA.2. Replicate measurements from one of three independent experiments. Fold change over mock is shown. Mean ± s.e.m. or individual datapoints are shown. For f, ***P < 0.0001. Source data

Discussion


We propose a model in which the earliest host innate immune responses make an important contribution to SARS-CoV-2 transmission by influencing whether interactions with the first few cells in the airway establish a productive infection. In this model, viruses with enhanced ability to evade or antagonize innate immunity, for example, through increased Orf6 and N expression, will transmit with greater frequency because they are better at avoiding inducing, or better at shutting down, the host responses that suppress this earliest replication. This model is supported by longitudinal nasal sampling of SARS-CoV-2-infected patients shortly after confirmation of infection, which revealed pronounced and early upregulation of an innate immune response in epithelial cells that rapidly declines after symptom onset44.

How early viral manipulation of the host innate immune response influences disease is less clear. We hypothesize that, once infection of the airway is irrevocably established, innate immune suppression that permits greater levels of viral replication may in turn lead to increased disease, simply due to greater viral burden and greater inflammatory responses. Concordantly, higher baseline antiviral gene expression and more potent innate induction in the nasal epithelium of children are associated with less severe infection outcomes compared with adults45. Like others, we assume this is explained by reduced viral loads reducing disease and early IFN protecting against transmission, with late IFN responses contributing to disease46. Similarly, inborn errors of innate antiviral mechanisms and IFN autoantibodies are associated with severe coronavirus disease 2019 (COVID-19)47,48,49,50,51, assumed to be explained by greater viral loads driving increased inflammatory disease. Furthermore, clinical trials of JAK/STAT inhibitors reduced COVID-19 mortality after hospitalization52. Considering an unrelated virus, simian immunodeficiency virus in macaques, may be relevant. Here transmission efficiency and subsequent disease are also influenced by IFN at the site of infection53. In all these examples, early IFN is beneficial, reducing transmission, but late IFN is bad, increasing symptoms. Human SARS-CoV-2 challenge studies are expected to help us understand the effect of these dynamics and innate immune contributions to transmission and disease by permitting sampling before exposure and during the earliest timepoints post infection with careful assessment of disease in a highly controlled environment54,55.

We have focused on changes in expression of N and Orf6 but we expect that other viral genes contribute to evasion of innate immunity and adaptation to humans. In contrast to common cold coronaviruses, SARS-CoV-2 and its relatives encode a broad range of accessory genes56,57 that antagonize innate immunity and probably contribute to effective transmission between species. Our data suggest that upregulation of Orf6 expression is a central feature of SARS-CoV-2 adaptation to humans. Our observations using Orf6-deletion viruses confirm Orf6 to be a potent viral innate immune antagonist, as reported by others32,34,35,39, and are consistent with a model in which, like Alpha, Omicron subvariant enhancement of Orf6 expression contributes to the reduced innate immune response to infection compared with earlier Omicron viruses. Orf6 upregulation by BA.5 may, in part, explain increased pathogenicity in vivo3,4. This notion is supported by ΔOrf6 SARS-CoV-2 infection of transgenic mice or hamsters, where the Orf6 mutant causes less severe disease and there is quicker recovery from infection, despite comparable viral loads in nose and lungs32,58. Expression of accessory and structural proteins as sgRNAs during SARS-CoV-2 replication provides an elegant mechanism to selectively regulate their abundance during adaptation to host, as the level of each sgRNA and thus protein can be independently adjusted by mutation, as we found for VOCs Alpha to Delta1,2.

The earliest Omicron subvariants BA.1 and BA.2 outcompeted Delta despite not enhancing innate immune antagonism, explained by extensive antibody escape1,2,3,4,5,6,7,8,9,10,11,12,13,14 and improved spike function/stability59,60. This suggests that adaptive immunity was the strongest selection force for Omicron emergence and global dominance. We hypothesize that the acquisition of enhanced innate immune suppression by Omicron lineage variants after their initial emergence required selection for improved transmission and dominance. Thus, innate immune escape may be the second dominant selective force the virus experiences after escape from neutralizing antibodies in a population with pre-existing immunity from prior infection and vaccination. We propose that evolving to better manage host innate immunity for improved transmission is a central feature of species-specific host adaptation for all emerging viruses. Intriguingly, SARS-CoV-2 continues to jump species barriers and has been detected infecting 34 different animal species so far61, illustrating its remarkable capacity to universally antagonize species specific innate immune responses. SARS-CoV-2 will be a fantastic model to further dissect species barriers to zoonotic spillovers and understand how viruses adapt to new species.

We propose that adaptation in spike and beyond also contributes to enhanced replication in human cells1,2. This may be important for outpacing early innate responses during transmission particularly in environments with a mix of permissive and non-permissive cells such as the upper human airways in which ACE2 is expressed only on ciliated cells62. Indeed, we have found that SARS-CoV-2 replicates more slowly in primary HAE cultures than in Calu-3 cells and that HAEs better recapitulate VOC replication advantages1,2. Primary HAEs complement more tractable monoculture models, such as Calu-3 that allow mechanistic studies. We propose that linking VOC genotype to phenotype in multiple models will be essential for effective prediction of novel variant behaviour. Moreover, understanding how adaptive changes in spike, leading to altered viral tropism, influence innate immune responses also warrants further study.

This study adds to the body of evidence for innate immunity being a key barrier that must be overcome by all pandemic zoonotic viruses, particularly in the absence of immune memory in an exposure-naive species. This has also been elegantly demonstrated recently for influenza virus where avian, but not human influenza virus, is efficiently restricted by human BTN3A3 (ref. 63), which like MX1 (ref. 64) can be overcome by adaptation to the human host. Innate immune evasion has also been linked to the single pandemic human immunodeficiency virus-1 lineage65. Our findings herein have broad implications for understanding zoonotic pathogen emergence because they reveal molecular details of how SARS-CoV-2 Omicron subvariants have achieved dominance, unexpectedly by increasing specific protein expression rather than adapting by protein coding mutation. Crucially, they suggest that improvements in innate immune evasion can continue to enhance transmission, even after establishment in humans. We hypothesize an inevitable ongoing trajectory of adaptation towards escape from the innate immune mechanisms that are the gatekeepers of transmission success.

Methods


Cell culture

Calu-3 cells were purchased from AddexBio (C0016001), Caco-2 cells were a kind gift from Dalan Bailey (Pirbright Institute), Hela-ACE2 cells were a gift from James E. Voss66 and A459 cells expressing ACE2 and TMPRSS2 were previously described22. Cell lines were cultured in Dulbecco’s modified Eagle medium (DMEM) supplemented with 10% heat-inactivated foetal bovine serum (FBS, Labtech) and 100 U ml−1 penicillin–streptomycin. Cells were passaged at 80–90% confluence. For infections, Calu-3 and Caco-2 cells were seeded at 2 × 105 cells ml−1 and Hela-ACE2 cells at 1 × 105 cells ml−1 and grown to 60–80% confluence for experiments1,19. Primary normal (healthy) bronchial epithelial (NHBE-A) cells from two independent donors were cultured for five to seven passages and differentiated at an air–liquid interface as previously described1. After 21–24 days of differentiation, cells were used in infection experiments. Experiments were performed without blinding or randomization.

Viruses

SARS-CoV-2 lineages Alpha (B.1.1.7), Delta (B.1.617.2)21 and Omicron (lineage B.1.1.529.1/BA.1, lineage B.1.1.529.2/BA.2, lineage BA.2.75 (BA.2.75.3) lineage BQ.1.1 (BQ.1.1.1), lineage XBB.1) isolates were a gift from Wendy Barclay (Imperial College London, United Kingdom). Omicron BA.4 (lineage B.1.1.529.4), BA.5 (lineage B.1.1.529.5), BQ.1.1 (B) (BQ.1.1.15) and lineage XBB.1.5 (XBB.1.5.13) were a gift from Alex Sigal and Khadija Khan (Africa Health Research Institute, Durban, South Africa)7,14. SARS-CoV-2 BA.5 (B) (SARS-CoV-2/Norway/20365/2022) was obtained from the Norwegian Institute of Public Health, Oslo, Norway. Omicron isolate identity was confirmed by full genome sequencing and assigned by Nextclade v.2.14.1 (https://clades.nextstrain.org)67,68. Alpha Orf6 deletion virus (Alpha ΔOrf6) was achieved by mutation of the first two methionines: M1L (A27216T) and M19L (A27200T). Reverse genetics-derived viruses were generated as previously described69,70. In brief, to generate the WT SARS-CoV-2 Alpha variant, a set of overlapping viral genomic complementary DNA fragments were chemically synthesized (GENEWIZ). The cDNA fragment representing the 5′ terminus of the viral genome contained the bacteriophage T7 RNA polymerase promoter and the fragment representing the 3′ terminus contained the T7 RNA polymerase termination sequences. These fragments were then assembled into a full-length Alpha cDNA genome using the transformation-associated recombination (TAR) in yeast method69. To generate the Alpha virus carrying the ATG codon changes (M1L and M19L) in its Orf6 gene (to generate Alpha ΔOrf6), the relevant cDNA fragments were chemically synthesized (ThermoFisher) and the mutant viral genome assembled using TAR in yeast as described above. We similarly generated WT BA.5, BA.5 ΔOrf6 (carrying M1L and M19L changes), and BA.5 Orf6 D61L (generated by introducing the GAT → CTC nucleotide change found in BA.2) using TAR in yeast except that the assembled cDNA genomes were placed under the control of the human cytomegalovirus promoter and the relevant termination sequences. The assembled WT and Orf6 null mutant genomes were transfected into BHK-hACE2-N cells stably expressing the SARS-CoV-2 N and the human ACE2 gene for virus rescue71. The rescued viruses were passaged once (P1 stock) in Vero.E6 cells and their full genomes sequenced using Oxford Nanopore as previously described72. For Alpha and BA.5 the RG-derived viruses are referred to as WT, ΔOrf6 or Orf6 D61L to differentiate them from the clinically isolated viruses used in all other experiments. All viruses were propagated by infecting Caco-2 cells in DMEM culture medium supplemented with 1% FBS and 100 U ml−1 penicillin–streptomycin at 37 °C as previously described1,19. Virus was collected at 72 h.p.i. and clarified by centrifugation at 2,100g for 15 min at 4 °C to remove any cellular debris. Virus stocks were aliquoted and stored at −80 °C. Virus stocks were quantified by extracting RNA from 100 μl of supernatant with 1 μg ml−1 carrier RNA using Qiagen RNeasy clean-up RNA protocol, before measuring viral E RNA copies ml−1 by RT–qPCR1,19. For experiments including Omicron subvariants XBB.1 and BQ.1.1, stocks and viral replication were quantified using Nsp12 RNA copies due to accumulation of mutations in the E gene of these variants, including in the region detected by our RT–qPCR assay. Virus titres were determined by TCID50 in Hela-ACE2 cells. A total of 104 cells were seeded in 96-well plates in 100 μl. The next day, seven tenfold serial dilutions of each virus stock or supernatant were prepared and 50 µl was added to the cells in quadruplicate. Cytopathic effect (CPE) was scored at 48–72 h.p.i. TCID50 ml−1 was calculated using the Reed and Muench method, and an Excel spreadsheet created by B. D. Lindenbach was used for calculating TCID50 ml−1 values73.

To generate SARS-CoV-2 lineage frequency plots for BA.1 (B.1.1.529.1), BA.2 (B.1.1.529.2), BA.4 (B.1.1.529.4), BA.5 (B.1.1.529.5), BA.2.75 (B.1.1.529.2.75), BQ.1.1 (B.1.1.529.5.3.1.1.1.1.1.1), XBB.1 and XBB.1.5 (Fig. 4a and Extended Data Fig. 6a), the number of samples sequenced per week worldwide over all time was extracted for each variant on 5 August 2023 from CoV-Spectrum (cov-spectrum.org)74 using genomic data from the Global Initiative on Sharing All Influenza Data (GSAID)75.

Virus culture and infection

For infections, inoculum was calculated using E copies per cell quantified by RT–qPCR. Cells were inoculated with indicated variants for 2 h at 37 °C, subsequently washed with phosphate-buffered saline (PBS), and fresh DMEM culture medium supplemented with 1% FBS and 100 U ml−1 penicillin–streptomycin was added. At the indicated timepoints, cells were collected for analysis. For primary HAE infections, virus was added to the apical side for 2–3 h at 37 °C. Supernatant was then removed, and cells were washed twice with PBS. All liquid was removed from the apical side, and basal medium was replaced with fresh Pneumacult ALI medium for the duration of the experiment. Virus release was measured at the indicated timepoints by extracting viral RNA from apical PBS washes. For poly(I:C) (Sigma) stimulations, cells were transfected with poly(I:C) using Lipofectamine2000 (InvitroGen) in Opti-Mem (Thermo) for the indicated times. For IFN-sensitivity assays, cells were pre-treated with indicated concentrations or recombinant human IFNβ (Peprotech) for 18 h before infection. Cytokines were maintained throughout the experiment. For inhibition assays, cells were pre-treated with 5 μM ruxolitinib (Cambridge Bioscience), 25 μM camostat (Apexbio), 25 μM E64d (Focus Biomolecules) or dimethyl sulfoxide control for 2–3 h before SARS-CoV-2 infection. Inhibitors were maintained throughout the infection.

RT–qPCR of host and viral gene expression in infected cells

Infected cells were lysed in RLT (Qiagen) supplemented with 0.1% β-mercaptoethanol (Sigma). RNA extractions were performed according to the manufacturer’s instructions using RNeasy Micro Kits (Qiagen) including on-column DNAse I treatment (Qiagen). cDNA was synthesized using SuperScript IV (Thermo) with random hexamer primers (Thermo). RT–qPCR was performed using Fast SYBR Green Master Mix (Thermo) for host gene expression and sgRNA expression or TaqMan Master mix (Thermo Fisher Scientific) for viral RNA quantification, and reactions were performed on the QuantStudio 5 Real-Time PCR systems (Thermo Fisher Scientific). Viral E RNA copies were determined as described previously1,19. Viral sgRNAs were detected using the same forward primer against the leader sequence paired with a sgRNA specific reverse primer1,76,77. Using the 2−ΔΔCt method, sgRNA levels were normalized to GAPDH to account for differences in RNA loading and then normalized to the level of Orf1a gRNA quantified in the same way for each variant to account for differences in the level of infection. Host gene expression was determined using the 2−ΔΔCt method and normalized to GAPDH expression. The following probes and primers were used:

GAPDH forward: 5′-ACATCGCTCAGACACCATG-3′, reverse: 5′-TGTAGTTGAGGTCAATGAAGGG-3′; IFNB forward: 5′-GCTTGGATTCCTACAAAGAAGCA-3′, reverse: 5′-ATAGATGGTCAATGCGGCGTC-3′; CXCL10 forward: 5′-TGGCATTCAAGGAGTACCTC-3′, reverse: 5′-TTGTAGCAATGATCTCAACACG-3′; IFIT1 forward: 5′-CCTCCTTGGGTTCGTCTACA-3′, reverse: 5′-GGCTGATATCTGGGTGCCTA-3′; IFIT2 forward: 5′-CAGCTGAGAATTGCACTGCAA-3′, reverse: 5′-CGTAGGCTGCTCTCCAAGGA-3′; MX1 forward: 5′-ATCCTGGGATTTTGGGGCTT-3′, reverse: 5′-CCGCTTGTCGCTGGTGTCG-3′; MX2 forward: 5′-CAGCCACCACCAGGAAAC-3′, reverse 5′-TTCTGCTCGTACTGGCTGTACAG-3′, RSAD2 forward: 5′-CTGTCCGCTGGAAAGTG-3′, reverse: 5′-GCTTCTTCTACACCAACATCC-3′; DDX58 forward: 5′-CTGGACCCTACCTACATCCTG-3′, reverse: 5′-GGCATCCAAAAAGCCACGG-3′. SARS-CoV-2 E Sarbeco forward: 5′- CGTTAATAGTTAATAGCGTACTTCTTTTTC-3′; SARS-CoV-2 E Sarbeco Probe1: 5′-FAM-ACACTAGCCATCCTTACTGCGCTTCG-TAMRA-3′; SARS-CoV-2 E Sarbeco reverse 5′-ATATTGCAGCAGTACGCACACA-3′; SARS-CoV-2 Nsp12 forward: 5′-GAGTGAAATGGTCATGTGTGG-3′; SARS-CoV-2 Nsp12 reverse: 5′-CATTGGCCGTGACAGCTTGAC-3′; SARS-CoV-2 Nsp12 Probe: 5′-CTCATCAGGAGATGCCACAACTGCTTATGCTAATAG-3′; 5′ Leader forward: 5′-ACCAACCAACTTTCGATCTCTTGT-3′; Orf1a reverse: 5′-CCTCCACGGAGTCTCCAAAG-3′; Orf6 reverse: GAGGTTTATGATGTAATCAAGATTC; N reverse: 5′-CCAGTTGAATCTGAGGGTCCAC-3′; Orf3a reverse: 5′-GCAGTAGCGCGAACAAAAT-3′; S reverse: 5′-GTCAGGGTAATAAACACCACGTG-3′.

Flow cytometry

Adherent cells were trypsinized and fixed in 4% formaldehyde before intracellular staining for SARS-CoV-2 nucleocapsid (N) protein. For N detection, cells were permeabilized for 15 min with Intracellular Staining Perm Wash Buffer (BioLegend) and subsequently incubated with 1 μg ml−1 CR3009 SARS-CoV-2 cross-reactive antibody (a gift from Laura McCoy) for 30 min at room temperature. Primary antibodies were detected by incubation with secondary AlexaFluor 488-Donkey-anti-Human IgG (Jackson Labs). All samples were acquired on a BD Fortessa X20 or LSR II using BD FACSDiva software. Data were analysed using FlowJo v10.6.2 (Tree Star). Gating strategy is shown in Extended Data Fig. 7.

Cytokine secretion

Secreted mediators were detected in cell culture supernatants by ELISA. IFNβ, IFNλ1/IFNλ3 and CXCL10 were measured using Human IFN-β Quantikine ELISA Kit, Human IL-29/IL-28B (IFNλ1/IFNλ3) DuoSet ELISA or Human CXCL10/IP-10 DuoSet ELISA reagents (Bio-Techne R&D Systems) according to the manufacturer’s instructions.

Western blotting

For detection of N, Orf6, Orf9b, spike and β-actin expression, whole-cell protein lysates were extracted with RIPA buffer, and then separated by sodium dodecyl sulfate–polyacrylamide gel electrophoresis, transferred onto nitrocellulose and blocked in PBS with 0.05% Tween 20 and 5% skimmed milk. Membranes were probed with rabbit-anti-SARS spike (Invitrogen, PA1-411-1165, 1:1,000), mouse-anti-SARS-CoV-2 spike (GeneTex, 1A9, 1:1,000), rabbit-anti-Orf6 (Abnova, PAB31757, 1:1,000), rabbit-anti-Orf9b (ProSci, 9191, 1:1,000), CR3009 SARS-CoV cross-reactive human-anti-N antibody (a gift from Laura McCoy, UCL, 1:1,000), rabbit-anti-phospho-STAT1 (Ser727; Cell Signaling, cat. no. 9177, 1:1,000), rabbit-anti-phospho-STAT1 (Tyr701; Cell Signaling, cat. no. 9167, clone 58D6, 1:1,000), rabbit-anti-STAT1 (Cell Signaling, cat. no. 9172, 1:1,000), rabbit-anti-IRF3 (Cell Signaling, cat. no. 4302, 1:1,000), rabbit-anti-phospho-IRF3 (Cell Signaling, cat. no. 29047, clone D6O1M, 1:1,000) and rabbit-anti-β-actin (A2066, Sigma, 1:2,500), followed by IRDye 800CW or 680RD secondary antibodies (Abcam, goat anti-rabbit, goat anti-mouse or goat anti-human, 1:10,000). Blots were imaged using an Odyssey Infrared Imager (LI-COR Biosciences) and analysed with Image Studio Lite software. Quantifications were performed to loading controls run on the same membrane as the protein of interest. For virion blots, live virus normalized by equal total E copies was purified across a 25% sucrose cushion and concentrated by centrifugation (2 h 16,500g, 4 °C).

Immunofluorescence staining and image analysis

Infected cells were fixed using 4% paraformaldehyde/formaldehyde for 1 h at room temperature and subsequently washed with PBS. A blocking step was carried out for 35 h at room temperature with 10% goat serum/1% bovine serum albumin/0.001 Triton X-100 in PBS. Double-stranded RNA (dsRNA) and nucleocapsid detection were performed by primary incubation with rabbit-anti-IRF3 antibody (sc-33641, Santa Cruz, 1:100), rabbit-anti-STAT1 (Cell Signaling, cat. no. 14994, clone D1K9Y, 1:100), mouse-anti-dsRNA (MABE1134, Millipore, 1:100) and CR3009 SARS-CoV cross-reactive human-anti-N antibodies (1:1,000) for 18 h and washed thoroughly in PBS. Primary antibodies detection occurred using secondary anti-rabbit-AlexaFluor-488, anti-mouse-AlexaFluor-568 and anti-human-Alexa647 conjugates (Jackson ImmunoResearch, 1:500) for 1 h. All cells were labelled with Hoechst33342 (H3570, Thermo Fisher, 1:5,000). Images were acquired using the WiScan Hermes 7-Colour High-Content Imaging System (IDEA Bio-Medical) at magnification 10×/0.4 numerical aperture. Four-channel automated acquisition was carried out sequentially. Images were acquired across a well area density resulting in 31 fields of view per well and ~20,000 cells. Images were pre-processed by applying a batch rolling ball background correction in FIJI ImageJ software package78 before quantification. IRF3 and STAT1 translocation analysis was carried out using the Athena Image analysis software (IDEA Bio-Medical) and data post-processed in Python. For dsRNA, infected cell populations were determined by thresholding of populations with more than two segmented dsRNA punctae. For transcription factor translocation analysis, infected populations were determined by presence of segmented nucleocapsid signal within the cell.

Image pre-processing was carried out using a custom macro applying a 30-pixel rolling ball background subtraction to all channels. Single-cell automated image analysis was carried out using the Athena image analysis software ‘Nuclear Translocation Assay’79 (IDEA-BioMedical). Within the Athena software, nuclei were segmented using the Hoechst33342 channel, and dsRNA/N channels were segmented as ‘cytoplasmic granules’ thresholded according to the mock infected population to identify infected cells. The cellular periphery was segmented by STAT1/IRF3 channels. The raw single-cell data were processed in a Python 3 script using the Pandas Data analysis library (https://pandas.pydata.org). In short, the mean nuclear:cytoplasmic ratio was calculated from the raw data. The data were ‘top and tail’ filtered, dropping the lowest and highest percentile for the following metrics: cell area, nuclear area, mean nuclear intensity and mean cytoplasmic intensity (STAT1/IRF3). The data were filtered into ‘infected cells’ by the presence of segmented cytoplasmic granules (dsRNA/N) or ‘bystander cells’ for their absence. The filtered data were then randomly sampled in the same Python environment. ImageJ macro and Python post-processing pipelines are available upon request.

Statistical analysis

Statistical analysis was performed using GraphPad Prism9, and details of statistical tests used are indicated. Data distribution was assumed to be normal unless stated differently, but this was not formally tested. No statistical methods were used to pre-determine sample sizes, but our sample sizes are similar to those reported in previous publications1,19. Data collection and analysis were not performed blind to the conditions of the experiments.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.


Quelle : nature

Den ersten Kommentar schreiben.

Antworten

Sie müssen eingeloggt sein, um einen Kommentar zu schreiben.

Der Antigen-Schnelltest ist ein wichtiger Bestandteil der aktuellen Gesundheitspraxis. Ob für den Einsatz zu Hause, in professionellen Umgebungen oder auf Reisen – eine verlässliche Wahl.