Next Article in Journal
Andrias davidianus Ranavirus (ADRV) Genome Replicate Efficiently by Engaging Cellular Mismatch Repair Protein MSH2
Next Article in Special Issue
Two Separate Clusters of SARS-CoV-2 Delta Variant Infections in a Group of 41 Students Travelling from India: An Illustration of the Need for Rigorous Testing and Quarantine
Previous Article in Journal
Role of Microbiota in Viral Infections and Pathological Progression
Previous Article in Special Issue
Prevalence of Anti-SARS-CoV-2 Antibodies and Potential Determinants among the Belgian Adult Population: Baseline Results of a Prospective Cohort Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 5
10.3390/v14050951
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SARS-CoV-2 Virion Infectivity and Cytokine Production in Primary Human Airway Epithelial Cells

by
Thuc Nguyen Dan Do
,
Sandra Claes
,
Dominique Schols
,
Johan Neyts
* and
Dirk Jochmans
*
Department of Microbiology, Immunology and Transplantation, Rega Institute, KU Leuven, 3000 Leuven, Belgium
*
Authors to whom correspondence should be addressed.
Viruses 2022, 14(5), 951; https://doi.org/10.3390/v14050951
Submission received: 4 April 2022 / Revised: 26 April 2022 / Accepted: 28 April 2022 / Published: 2 May 2022
(This article belongs to the Special Issue SARS-CoV-2 Research in Belgium)
Browse Figures
Versions Notes

Abstract

:
The emergence of new SARS-CoV-2 variants and the replacement of preceding isolates have been observed through B.1.1.7, B.1.351, B.1.617.2, and B.1.1.529 lineages (corresponding to alpha, beta, delta, and omicron variants of concern (VoC), respectively). However, there is still a lack of biological evidence to which extent those VoC differ from the ancestral lineages. By exploiting human airway epithelial cell (HAEC) cultures, which closely resemble the human airway architecture and physiology, we report distinctive SARS-CoV-2 tropism in different respiratory tissues. In general, SARS-CoV-2 VoC predominantly infect and replicate in HAEC better than the progenitor USA-WA1 isolate or the BavPat1 isolate, which contains the D614G mutation, even though there is little to no difference between variants regarding their infectivity (i.e., virion-per-vRNA copy ratio). We also observe differential tissue-specific innate immunity activation between the upper and lower respiratory tissues in the presence of the virus. Our study provides better comprehension of the behavior of the different VoC in this physiologically relevant ex vivo model.

1. Introduction

Since SARS-CoV-2 was first identified, a new isolate containing the D614G spike (S) substitution was reported with higher infectivity and transmission in in vitro and in vivo models [1]. This 614G virus quickly dominated the global pandemic, especially in Europe. Later, in 2020, the B.1.1.7 (alpha) variant firstly discovered in the UK quickly became dominant worldwide, replacing the 614G lineage before it was usurped by the B.1.162.2 (delta) variant. Meanwhile, B.1.351 (beta) emerged in South Africa and was mainly confined there and its surrounding countries. At the end of 2021, B.1.1.529 (omicron) was first identified in South Africa and rapidly spread to other countries. Epidemiologically, not only is the human-to-human transmission of those isolates more superior than that of other lineages, but their potential escape of vaccine-acquired immunity and adverse impact on therapeutics categorizes them as variants of concern (VoC) [2,3,4,5,6]. Additionally, the delta and alpha variants have been associated with the increase in disease severity [7,8].
New mutations have been detected across VoC’s genomes, both outside and inside the S protein-coding region. As they share amino acid substitutions at the receptor-binding domain (RBD) of the S protein, such as D614G and N501Y, the B.1.1.7 and B.1.351 variants are resistant to the neutralization activity of monoclonal antibodies [9]. Other data indicates that B.1.351 reduces the protective efficacy of current vaccines more considerably than B.1.1.7 [10,11]. The K417N, E848K, and N501Y amino acid changes in key residues of the S protein seem to have greater influence on the neutralization of those variants than other mutations [10]. Overall, VoC and emergent variants [12] with similar mutated residues in the S protein may cause new resurgences of COVID-19 cases worldwide.
To address this, we assessed the susceptibility of different sources of human primary airway epithelial cells (HAEC), cultured at the air–liquid interface (ALI), to authentic alpha, beta, delta, and omicron viruses in comparison with their ancestral lineages. In addition, we analyzed the induction of cytokine production. Our study provides better understanding of the behavior of the different VoC in this physiologically relevant ex vivo model.

2. Materials and Methods

2.1. SARS-CoV-2 Isolates and Virus Stocks

The ancestral Severe Acute Respiratory Syndrome-related Coronavirus 2 (SARS-CoV-2) strains, Germany/BavPat1/2020 isolate (EPI_ISL_406862|28 January 2020, in this study referred to as BavPat1, kindly provided by Prof. Christian Drosten, University Hospital Charité, Berlin, Germany) and USA/WA-CDC-WA1/2020 isolate (EPI_ISL_404895|19 January 2020, referred to as USA-WA1, provided by BEI Resources) were retrieved from COVID-19 positive patients in early 2020, and have been described previously [13,14]. The USA-WA1 strain is recognized as lineage A and is most similar to the original Wuhan-Hu-1 strain (GeneBank accession number MN908947.3) whereas the BavPat1 is designated as lineage B and carries the D614G amino acid substitution in the spike protein.
The SARS-CoV-2 VoC B.1.1.7 (EPI_ISL_791333|21 December 2020), B.1.351 (EPI_ISL_896474|11 January 2021), B.1.617.2 (EPI_ISL_2425097|20 April 2021), and B.1.1.529 (EPI_ISL_6794907|24 November 2021) isolates were named according to their own lineages. All strains were isolated in-house from nasopharyngeal swabs of travelers returning to Belgium. The generation of virus stock was fully reported elsewhere [15,16,17]. Only early passages were used for the ALI experiment, i.e., passage 3 (P3) of BavPat1, P5 of USA-WA1, P2 of B.1.1.7, B.1.351, B.1.617.2, and B.1.1.529.
All infectious virus-related work was conducted under a biosafety level 3 (BSL-3) facility at the Rega Institute for Medical Research, KU Leuven, according to institutional guidelines.

2.2. Cell Culture

The African monkey kidney cell line VeroE6 (ATCC CRL-1586) was cultured in Dulbecco’s modified Eagle’s medium (DMEM; Gibco, catalogue no. 41965-039) supplemented with 10% v/v heat-inactivated fetal bovine serum (HI-FBS; HyClone, catalogue no. SV03160.03) and 1% v/v Penicillin-Streptomycin 10,000 U/mL (P/S; Gibco, catalogue no. 15140148) at 37 °C and 5% CO2. End-point titrations were performed with 2% HI-FBS-containing medium instead of 10%.

2.3. Airway Ex Vivo Models and Study Design

A total of 150 air–liquid interface cultures were used in this study. That included 60 human nasal airway epithelial cell cultures (HnAEC; catalogue no. EP02MP) obtained from a pool of fourteen healthy donors (n = 9), 30 human bronchial airway epithelial cell cultures (HbAEC; catalogue no. EP01MD) from one donor (n = 6), and 60 human small airway epithelial cell cultures (HsAEC; catalogue no. EP21SA) from two different donors (n = 12). HAEC were received from Epithelix (Geneva, Switzerland) in an air–liquid interphase set-up. Cultures model the architecture and cellular complexity of the conducting airway at specific regions and are ready-to-use. After arrival, each insert was washed with pre-warmed medium and maintained in corresponding MucilAir medium (Epithelix, catalogue no. EP04MM) or SmallAir medium (Epithelix, catalogue no. EP64SA) at 37 °C and 5% CO2 for at least 4 days before use. On the day of the experiment, the HAEC were exposed to 100 µL of SARS-CoV-2 inoculum at the apical side for 1.5 h at 35 °C. The first apical wash with medium was collected 24 h later (day 1 post-infection (p.i.)). Every day from day 0, subsequent apical washes were collected and stored at −80 °C for analysis.

2.4. RNA Extraction and Quantitative Reverse Transcription-PCR (RT-qPCR)

Viral RNA in the apical wash was isolated using the Cells-to-cDNATM II cell lysis buffer kit (Thermo Fisher Scientific, Merelbeke, Belgium, catalogue no. AM8723). Briefly, 5 µL wash fluid was added to 50 µL lysis buffer, incubated at room temperature (RT) for 10 min and then at 75 °C for 15 min. 150 µL nuclease-free water was additionally added to the mixture prior to RT-qPCR. In parallel, a ten-fold serial dilution of 2019_nCoV_N_positive control from IDT Technologies (catalogue no. 10006625) was used to generate a standard curve. The amount of viral RNA expressed as number of copies per mL (copies/mL) was quantified by RT-qPCR on a LighCycler96 platform (Roche, Anderlecht, Belgium) using the iTaq universal probes one-step kit (Bio-Rad, Temse, Belgium, catalogue no. 1725141) and a commercial mix of primers for the N gene (forward primer 5′-GACCCCAAAATCAGCGAAAT-3′, reverse primer 5′-TCTGGTTACTGCCAGTTGAATCTG-3′) and probes (5′-FAM-ACCCCGCATTACGTTTGGTGGACC-BHQ1-3′) manufactured at IDT Technologies (catalogue no. 10006606). The reaction was prepared according to the manufacturer’s protocol.

2.5. End-Point Titration Assay

VeroE6 cells were seeded in 96-well tissue culture plates (Corning, Amsterdam, The Netherlands, catalogue no. 353072) at a density of 1 × 104 cells per well. After 24 h, serial dilutions of ALI wash fluid were prepared in the plates. Cells were incubated for 3 days at 37 °C and evaluated microscopically for the absence or presence of the virus-induced cytopathic effect (CPE). The infectious viral titer was determined by end-point titration, expressed as 50% tissue culture infectious dose per mL (TCID50/mL). Virus titers were calculated by using the Reed and Muench method as previously reported [18].

2.6. Cell Viability Assay

At the end of the experiment, HAEC were washed once with pre-warmed 1× PBS (Thermo Fisher Scientific, catalogue no. 14190-094) before being incubated with 200 µL 0.25% trypsin-EDTA (Thermo Fisher Scientific, catalogue no. 25200056) at 37 °C for 2 h. Single cell suspension was stained with acridine orange and propidium iodide (AOPI) solution. The number of living cells and dead cells were automatically counted by a LunaTM automated cell counter (Logos Biosystems, Villeneuve d’Ascq, France)), giving the cell viability of HAEC as the percentage of living cells over total cell numbers.
In case of nasal HAEC cell growth was determined via a colorimetric method using the CYQUANTTM LDH Cytotoxicity Assay Kit (Thermo Fisher Scientific, catalogue no. C20301). After being washed apically, cultures were refreshed with new basolateral medium and lysed for 2 h. The activity of released lactate dehydrogenase (LDH) in the basal medium was proportionally quantified through the level of formazan formation at 490 nm. Cell viability was expressed as the percentage of LDH activity in infected groups compared to the non-infected one.

2.7. Profiling of Chemokine and Cytokine Secretion

Apical wash fluid and/or basolateral medium were analyzed for a panel of chemokines and cytokines—IP-10, IL-6, IL-10, IFN-ɣ, and TNF-α using customized ProcartaPlex assay kits (Thermo Fisher ScientificInvitrogen) according to the manufacturer’s protocol. Those cytokines have been reported as being linked to disease exacerbation in patients with COVID-19 [19]. The analyte concentrations were determined using the included standard curve of known concentrations coupled with the Luminex XPONENT v4.2 software.

2.8. Statistical Analysis

The number of cultures and independent experiments that were performed is indicated in each legend of the individual figure. Results are represented as median and 95% confidence intervals for cytokine and chemokine expressions or mean ± standard deviation in other experiments. All statistical comparisons in the study were performed in GraphPad Prism v9.2.0 (GraphPad Software, San Diego, CA, USA). Statistical significance was determined using Kruskal–Wallis with Dunnett’s multiple comparison test unless specifically mentioned otherwise. p-values of ≤0.05 were considered statistically significant. In the figures, “ns” indicates a p > 0.05 and an asterisks indicates a statistical significance level of * p < 0.05, ** p < 0.01, and *** p < 0.001.

3. Results

To investigate the relative infectivity and cytokine secretion profile of VoC B.1.1.7 and B.1.351, we infected HAEC from a small airway origin (HsAEC) with 2 × 103 TCID50 of each SARS-CoV-2 isolate including the ancestral isolates USA-WA1 and BavPat1. As expected, only 17% of the cultures (2 out of 12) inoculated with the USA-WA1 isolate resulted in productive infection while 42%, 92%, and 75% of samples showed a productive infection when inoculated with BavPat1, B.1.1.7, and B.1.351 isolates, respectively (n = 12) (Figure 1A). Only infected samples were presented in the next figures. B.1.1.7 and B.1.351 variants replicated more efficiently during the first 6 days of the experiment, reaching a peak of 109 RNA copies/mL (Figure 1B) and 105 TCID50/mL (Figure 1C) while only 106 and 107 RNA copies/mL and ~103 TCID50/mL were recorded for the comparators (USA-WA1 and BavPat1, respectively). Although it was noticed that infectious particles were significantly produced on day 6 in B.1.351-infected inserts (p = 0.0476; Figure 1H), no difference in relative infectivity was observed for all four VoC (p = 0.9181; Figure 1F) (p = 0.7611; Figure 1I). All HAEC-cultures maintained more than 90% cell viability (Figure 1J).
Likewise, a similar observation was made using HAEC of bronchial origin (HbAEC) (Figure 2). Whereas six out of six samples were infected by either the B.1.1.7 or the B.1.351 variant, only 50% of the cultures could be infected with either USA-WA-1 or BavPat1 isolates (n = 6) (Figure 2A). At 2 d.p.i., the viral RNA loads were determined at 2 × 105, 7 × 105, 1 × 108, and 2 × 106 RNA copies/mL (Figure 2D) and the infectious virus titers were 5 × 104, 1 × 105, 5 × 106, and 3 × 105 TCID50/mL (Figure 2E) in USA-WA1, BavPat1, B.1.1.7, and B.1.351 groups, respectively. As a result, no disparity in the virion infectivity was observed (p = 0.1169; Figure 2F). Of note, the alpha-infected cultures produced a lower infectivity ratio in comparison to the ancestral strain USA-WA1 on day 4 p.i. (p = 0.0162; Figure 2I). Cultures incubated with any isolate maintained at least 90% cell viability after the 10-day experiment (Figure 2J).
To further comprehend the immune response of the HAEC towards SARS-CoV-2 infection, the expression levels of cytokines related to COVID-19 in patients were analyzed in the apical washes of the cultures. Only for the B.1.351 and B.1.1.7 isolates we could see a difference in cytokine production compared to the uninfected controls. In the B.1.351-infected group, an upregulation of IP-10 levels of 5.5 fold (p = 0.0044; Figure 3A) to 7.5 fold (p = 0.0004; Figure 4A) compared to the uninfected group was observed at 6 days p.i. A 4.6-fold increase in IP-10 expression was induced by the B.1.1.7 variant on day 4 p.i. (p = 0.0114; Figure 4A). No upregulation in IL-6 production was observed (Figure 3B, Figure 4B). Likewise, no difference in other cytokine levels (i.e., IL-10, IFN-ɣ, TNF-α) was noted (data not shown). The overall discrepancy in cytokine production between the alpha and beta variants versus basal strains was limited.
Since the delta (B.1.617.2) and omicron (B.1.1.529) variants became dominant worldwide in mid and late 2021, replacing the preceding isolates in most countries, we performed an additional experiment exploring SARS-CoV-2 VoC infection using HAEC from nasal regions. After 24-h p.i. (h.p.i.), viral RNA of all SARS-CoV-2-infected groups were detectable with a large replication advantage in omicron-infected samples (Figure 5A). At 2 d.p.i., the omicron isolate produced the highest quantity of vRNA (2.5 × 1010 RNA copies/mL, p < 0.0001) and virions (5 × 106 TCID50/mL, p = 0.0066) when compared with USA-WA1 (6.3 × 107 RNA copies/mL and 4 × 104 TCID50/mL) and the alpha (1010 RNA copies/mL, p = 0.0020; 6.3 × 106 TCID50/mL, p = 0.0041) and delta (4 × 109 RNA copies/mL, p = 0.0366; 3 × 106 TCID50/mL, p = 0.0303) isolates (Figure 5B). Again, comparable infectivity of the virions produced by the different groups was noted (Figure 5B). There was no influence of SARS-CoV-2 infection on HnAEC cell viability (Figure 5C). Interestingly, infection with either omicron, delta, or BavPat1 led to an upregulation of IP-10 and/or IL-6 expression in apical washes and/or basal medium compared to non-infected samples at 48 h.p.i. (Figure 6). For instance, IP-10 and IL-6 levels were both induced in the B.1.1.529-infected group while only the former was upregulated in B.1.617.2-infected cultures, with a statistical significance compared to the negative control (p < 0.0001 for B.1.1.529/IP-10, p = 0.0179–0.0068 for B.1.617.2/IP-10, p = 0.0429–0.0469 for B.1.1.529-IL-6). From day 4 p.i., the production of IP-10 and IL-6 were markedly accelerated in most of SARS-CoV-2-infected cultures, but not in the alpha and delta ones. Taken together, comparable infectivity was observed between samples infected with VoC and ancestral isolates, and alpha, beta, omicron, and delta variants have different effects on cytokine secretions in the human respiratory tract (HRT) at the onset of SARS-CoV-2 infection.
By comparing the replication kinetics of SARS-CoV-2 strains in different inserts generated from distinctive regions of the HRT, we observed that relative virion infectivity was not only tissue-dependent (p < 0.0001; two-way ANOVA) but also SARS-CoV-2 variant-dependent (p = 0.0030, two-way ANOVA). As infected samples from small airway tissues did not effectively produce detectable RNA copies and/or infectious progenies after 2 d.p.i., we excluded this data from the HsAEC experiments for further analysis.
In general, bronchial tissues produced a higher virion-per-vRNA copy ratio (around 1–2 logs increase) than the nasal tissues (Figure 7A, Table 1). In addition, the beta variant generated higher infectivity than the alpha or USA-WA1 whilst the omicron isolate produced the lowest ratio of virion-per-vRNA copy in nasal cultures. The VoC (alpha and beta) were inclined to have lower virion:vRNA copy levels than the ancestral lineages, but no difference between isolates was noted in bronchial samples (Figure 7B).

4. Discussion

The human respiratory tract (HRT), the primary target of SARS-CoV-2 infection, extends from the nasal cavity to the terminal bronchioles and is fully covered by the airway epithelium. By forming a pseudostratified epithelial layer, the organotypic HAEC-model closely resembles the human airway architecture and physiology with its dynamic cellular composition including basal cells, ciliated cells, and goblet cells [20]. Although SARS-CoV-2 genetic diversity rapidly evolved, little is known about the comparison of VoC infection and virulence in a physiologically relevant human model, potentially explaining the current epidemiological observations. Here, we report distinctive SARS-CoV-2 tropism in different respiratory tissues. In general, the SARS-CoV-2 VoC predominantly infects and replicates in HAEC better than the ancestral USA-WA1 isolate or the ancestral BavPat1 isolate, which contains the D614G mutation. A possible explanation is that the heavily mutated spike (S) protein, including the P681R amino acid change in the delta variant, enhances the cleavage efficiency of S1 and S2 subunits, resulting in more productive virus entry into airway epithelial cells, thereby boosting the replication fitness of the mutant [21]. Another S substitution (i.e., N501Y seen in alpha and beta isolates) also confers a major fitness advantage in the upper airway replication and transmission via an increased receptor-binding domain (RBD)/angiotensin-converting enzyme 2 (ACE2) interaction [22]. Mutations in the nucleocapsid protein (R203K and G204R observed in omicron) is associated with increased viral replication [23]. To this end, our finding is in line with other reports [21,22,24,25]. In particular, the extremely fast replication of omicron could be due to its adapted cellular tropism (less transmembrane protease serine 2 (TMPRSS2) dependent), allowing the virus to infect more nasal cells via the endocytic cathepsin L pathway than the competitors [24,25,26].
When it comes to their infectivity (i.e., virion-per-vRNA copy ratio), there is little to no difference between variants in comparison with the ancestral USA-WA1 lineage. One possible reason for the discrepancy concerning the virion-per-vRNA copy ratio observed in alpha-infected cultures (Figure 2I) is that other undefined mutations around the furin cleavage site (e.g., T716I, N679K) modulate the proteolytic cleavage of the alpha S protein, thus, destabilizing the S trimer conformation, leading to unstable virions [26]. Our observation is in accordance with the replicative fitness of alpha lineage in HAEC previously reported [27]. Of note, since not all extracellular vRNA represents equal infectious progeny, it is of importance to compare the viral infectivity based on the virion:vRNA copy value rather than the sole parameter of vRNA copy or infectious titer. This perhaps gives better insights into the transmission dynamics of viral pathogens [28]. The lower virion:vRNA copy ratio is often attributed to incomplete genomes, lethal mutations, deficits in virion capsid/envelope structure, or host barriers which restrict early viral replication [29]. As the new SARS-CoV-2 variants evolve to more productively enter the airway epithelial cells [21,22], it is unclear why their virion-per-copy ratios are comparable to that of the progenitor USA-WA1 and to what extent that might explain the corresponding human-to-human transmission. Therefore, further experiments are required to elucidate the stability of VoC virion structure and its corresponding infectivity to explain the similar virion:vRNA copy observed in this study.
The epithelial composition of different anatomical regions in HRT has also been reported to partly influence SARS-CoV-2 replication and spread [20]. High co-expression levels of ACE2 and TMPRSS2 (two main entry receptors for SARS-CoV-2 infection) have been found in both secretory and ciliated cells [30,31], facilitating SARS-CoV-2 replication [32]. Notably, the co-expression of ACE2 and TMPRSS2 genes were more enriched in the nasal cavity than in the lower airway and parenchyma [31]. This explains why the highest percentage (100%) of infected samples in our study is from nasal tissue, followed by bronchial and small airway tissues.
In the current study, we observed the upregulation of proinflammatory cytokines and chemokines (IL-6, IP-10) in omicron-inoculated nasal tissues 48 h.p.i. This could be due to the extremely high viral replication of omicron already on day 1 p.i. Later, elevated production of IP-10 and IL-6 was also observed in other groups at a similar level, suggesting a saturated condition. A lesser influence on protein secretions was noted in the alpha and delta groups, implying variant-specific immune responses. Meanwhile, the high concentrations of IL-6 and IP-10 in patients’ sera, together with the infiltrating macrophages, have been depicted as the “immune signature” of SARS-CoV-2 and as the prototype factors driving COVID-19 progression [33]. On the other hand, it is reported that the systemic and mucosal immunity show distinct responses to SARS-CoV-2 infection, with a strong impact on the nasalpharyngeal cytokine profile [34]. This could perhaps explain the unusual cytokine release profile in nasal tissues, which is probably not the source of cytokines in the plasma of COVID-19 patients. In small airway and bronchial cultures, the induction of the epithelium-secreting IP-10 protein is not in accordance with the secretion of the IL-6 protein. Taken together, it is plausible that epithelial cells are conditioned to produce immune-associated proteins to reduce viral susceptibility. The restricted innate immune response to SARS-CoV-2 reported herein is in agreement with data from another study [35] and expands the observations in cytokine secretion in primary HAEC and in lung explants [35,36,37]. In this regard, the differential tissue-specific innate immune activation between the upper and lower respiratory tissues could also be attributed to the discrepancy in observed virion infectivity.
One of the limitations of our study is the lack of a direct competition experiment, which could expose replicative differences not detected in single growth kinetic assays, as previously reported [21,22,25,27,38].

5. Conclusions

To conclude, our findings reveal that despite the comparable virion infectivity, the different SARS-CoV-2 VoC are endowed with faster and enhanced replication competence in HAEC cultures compared with ancestral lineages, explaining their more efficient transmission in the human population. Moreover, a distinctive tropism of respiratory tissues in response to SARS-CoV-2 infection is a pivotal determinant in elucidating host–pathogen interactions or selecting a suitable model for antiviral testing.

Author Contributions

Conceptualization, J.N., D.J. and T.N.D.D.; methodology, T.N.D.D. and S.C.; formal analysis, T.N.D.D. and S.C.; investigation, T.N.D.D.; writing—original draft preparation, T.N.D.D.; writing—review and editing, J.N., D.J., D.S. and S.C.; visualization, T.N.D.D.; supervision, J.N. and D.J.; funding acquisition, J.N. and D.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by the Bill & Melinda Gates Foundation (BMGF) under grant number INV-006366 and the European Health Emergency Preparedness and Response Authority (HERA). T.N.D.D received the fellowship from the European Union’s Horizon 2020 grant OrganoVIR ‘Organoids for Virus Research—An innovative training—ITN programme’. Part of this research work was performed using the ‘Caps-It’ research infrastructure (project ZW13-02) that was financially supported by the Hercules Foundation and Rega Foundation, KU Leuven.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We thank Birgit Voeten and Niels Cremers for their excellent technical assistance. We also thank Piet Maes for kindly providing the SARS-CoV-2 isolates used in this study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Baric, R.S. Emergence of a Highly Fit SARS-CoV-2 Variant. N. Engl. J. Med. 2020, 383, 2684–2686. [Google Scholar] [CrossRef] [PubMed]
  2. Davies, N.G.; Abbott, S.; Barnard, R.C.; Jarvis, C.I.; Kucharski, A.J.; Munday, J.D.; Pearson, C.A.B.; Russell, T.W.; Tully, D.C.; Washburne, A.D.; et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science 2021, 372, eabg3055. [Google Scholar] [CrossRef] [PubMed]
  3. Tegally, H.; Wilkinson, E.; Giovanetti, M.; Iranzadeh, A.; Fonseca, V.; Giandhari, J.; Doolabh, D.; Pillay, S.; San, E.J.; Msomi, N.; et al. Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa. medRxiv 2020. [Google Scholar] [CrossRef]
  4. Washington, N.L.; Gangavarapu, K.; Zeller, M.; Bolze, A.; Cirulli, E.T.; Schiabor Barrett, K.M.; Larsen, B.B.; Anderson, C.; White, S.; Cassens, T.; et al. Emergence and rapid transmission of SARS-CoV-2 B.1.1.7 in the United States. Cell 2021, 184, 2587–2594.e7. [Google Scholar] [CrossRef]
  5. Farinholt, T.; Doddapaneni, H.; Qin, X.; Menon, V.; Meng, Q.; Metcalf, G.; Chao, H.; Gingras, M.C.; Avadhanula, V.; Farinholt, P.; et al. Transmission event of SARS-CoV-2 delta variant reveals multiple vaccine breakthrough infections. BMC Med. 2021, 19, 255. [Google Scholar] [CrossRef]
  6. Gu, H.; Krishnan, P.; Ng, D.Y.M.; Chang, L.D.J.; Liu, G.Y.Z.; Cheng, S.S.M.; Hui, M.M.Y.; Fan, M.C.Y.; Wan, J.H.L.; Lau, L.H.K.; et al. Probable Transmission of SARS-CoV-2 Omicron Variant in Quarantine Hotel, Hong Kong, China, November 2021. Emerg Infect. Dis 2022, 28, 460–462. [Google Scholar] [CrossRef]
  7. Davies, N.G.; Jarvis, C.I.; Group, C.C.-W.; Edmunds, W.J.; Jewell, N.P.; Diaz-Ordaz, K.; Keogh, R.H. Increased mortality in community-tested cases of SARS-CoV-2 lineage B.1.1.7. Nature 2021, 593, 270–274. [Google Scholar] [CrossRef]
  8. Sheikh, A.; McMenamin, J.; Taylor, B.; Robertson, C.; Public Health Scotland and the EAVE II Collaborators. SARS-CoV-2 Delta VOC in Scotland: Demographics, risk of hospital admission, and vaccine effectiveness. Lancet 2021, 397, 2461–2462. [Google Scholar] [CrossRef]
  9. Wang, P.; Nair, M.S.; Liu, L.; Iketani, S.; Luo, Y.; Guo, Y.; Wang, M.; Yu, J.; Zhang, B.; Kwong, P.D.; et al. Antibody resistance of SARS-CoV-2 variants B.1.351 and B.1.1.7. Nature 2021, 593, 130–135. [Google Scholar] [CrossRef]
  10. Liu, Y.; Liu, J.; Xia, H.; Zhang, X.; Fontes-Garfias, C.R.; Swanson, K.A.; Cai, H.; Sarkar, R.; Chen, W.; Cutler, M.; et al. Neutralizing Activity of BNT162b2-Elicited Serum. N. Engl. J. Med. 2021, 384, 1466–1468. [Google Scholar] [CrossRef]
  11. Mahase, E. Covid-19: Novavax vaccine efficacy is 86% against UK variant and 60% against South African variant. BMJ 2021, 372, n296. [Google Scholar] [CrossRef] [PubMed]
  12. Naveca, F.G.; Nascimento, V.; de Souza, V.C.; Corado, A.D.L.; Nascimento, F.; Silva, G.; Costa, Á.; Duarte, D.; Pessoa, K.; Mejía, M.; et al. COVID-19 in Amazonas, Brazil, was driven by the persistence of endemic lineages and P.1 emergence. Nat. Med. 2021, 27, 1230–1238. [Google Scholar] [CrossRef] [PubMed]
  13. Bohmer, M.M.; Buchholz, U.; Corman, V.M.; Hoch, M.; Katz, K.; Marosevic, D.V.; Bohm, S.; Woudenberg, T.; Ackermann, N.; Konrad, R.; et al. Investigation of a COVID-19 outbreak in Germany resulting from a single travel-associated primary case: A case series. Lancet Infect. Dis 2020, 20, 920–928. [Google Scholar] [CrossRef]
  14. Harcourt, J.; Tamin, A.; Lu, X.; Kamili, S.; Sakthivel, S.K.; Murray, J.; Queen, K.; Tao, Y.; Paden, C.R.; Zhang, J.; et al. Severe Acute Respiratory Syndrome Coronavirus 2 from Patient with Coronavirus Disease, United States. Emerg Infect. Dis 2020, 26, 1266–1273. [Google Scholar] [CrossRef] [PubMed]
  15. Abdelnabi, R.; Boudewijns, R.; Foo, C.S.; Seldeslachts, L.; Sanchez-Felipe, L.; Zhang, X.; Delang, L.; Maes, P.; Kaptein, S.J.F.; Weynand, B.; et al. Comparing infectivity and virulence of emerging SARS-CoV-2 variants in Syrian hamsters. EBioMedicine 2021, 68, 103403. [Google Scholar] [CrossRef]
  16. Abdelnabi, R.; Foo, C.S.; Jochmans, D.; Vangeel, L.; De Jonghe, S.; Augustijns, P.; Mols, R.; Weynand, B.; Wattanakul, T.; Hoglund, R.M.; et al. The oral protease inhibitor (PF-07321332) protects Syrian hamsters against infection with SARS-CoV-2 variants of concern. Nat. Commun 2022, 13, 719. [Google Scholar] [CrossRef]
  17. Abdelnabi, R.; Foo, C.S.; Zhang, X.; Lemmens, V.; Maes, P.; Slechten, B.; Raymenants, J.; Andre, E.; Weynand, B.; Dallmeier, K.; et al. The omicron (B.1.1.529) SARS-CoV-2 variant of concern does not readily infect Syrian hamsters. Antivir. Res. 2022, 198, 105253. [Google Scholar] [CrossRef]
  18. Reed, L.J.; Muench, H. A simple method of estimating fifty per cent endpoints 12. Am. J. Epidemiol. 1938, 27, 493–497. [Google Scholar] [CrossRef]
  19. Yang, L.; Xie, X.; Tu, Z.; Fu, J.; Xu, D.; Zhou, Y. The signal pathways and treatment of cytokine storm in COVID-19. Signal. Transduct Target. Ther. 2021, 6, 255. [Google Scholar] [CrossRef]
  20. Wang, Y.; Thaler, M.; Ninaber, D.K.; van der Does, A.M.; Ogando, N.S.; Beckert, H.; Taube, C.; Salgado-Benvindo, C.; Snijder, E.J.; Bredenbeek, P.J.; et al. Impact of human airway epithelial cellular composition on SARS-CoV-2 infection biology. bioRxiv 2021. [Google Scholar] [CrossRef]
  21. Liu, Y.; Liu, J.; Johnson, B.A.; Xia, H.; Ku, Z.; Schindewolf, C.; Widen, S.G.; An, Z.; Weaver, S.C.; Menachery, V.D.; et al. Delta spike P681R mutation enhances SARS-CoV-2 fitness over Alpha variant. bioRxiv 2021. [Google Scholar] [CrossRef]
  22. Liu, Y.; Liu, J.; Plante, K.S.; Plante, J.A.; Xie, X.; Zhang, X.; Ku, Z.; An, Z.; Scharton, D.; Schindewolf, C.; et al. The N501Y spike substitution enhances SARS-CoV-2 infection and transmission. Nature 2022, 602, 294–299. [Google Scholar] [CrossRef] [PubMed]
  23. Wu, H.; Xing, N.; Meng, K.; Fu, B.; Xue, W.; Dong, P.; Tang, W.; Xiao, Y.; Liu, G.; Luo, H.; et al. Nucleocapsid mutations R203K/G204R increase the infectivity, fitness, and virulence of SARS-CoV-2. Cell Host Microbe 2021, 29, 1788–1801.e6. [Google Scholar] [CrossRef] [PubMed]
  24. Hui, K.P.Y.; Ho, J.C.W.; Cheung, M.C.; Ng, K.C.; Ching, R.H.H.; Lai, K.L.; Kam, T.T.; Gu, H.; Sit, K.Y.; Hsin, M.K.Y.; et al. SARS-CoV-2 Omicron variant replication in human bronchus and lung ex vivo. Nature 2022, 603, 715–720. [Google Scholar] [CrossRef] [PubMed]
  25. Peacock, T.P.; Brown, J.C.; Zhou, J.; Thakur, N.; Newman, J.; Kugathasan, R.; Sukhova, K.; Kaforou, M.; Bailey, D.; Barclay, W.S. The SARS-CoV-2 variant, Omicron, shows rapid replication in human primary nasal epithelial cultures and efficiently uses the endosomal route of entry. bioRxiv 2022. [Google Scholar] [CrossRef]
  26. Du, X.; Tang, H.; Gao, L.; Wu, Z.; Meng, F.; Yan, R.; Qiao, S.; An, J.; Wang, C.; Qin, F.X. Omicron adopts a different strategy from Delta and other variants to adapt to host. Signal. Transduct. Target. Ther. 2022, 7, 45. [Google Scholar] [CrossRef] [PubMed]
  27. Touret, F.; Luciani, L.; Baronti, C.; Cochin, M.; Driouich, J.S.; Gilles, M.; Thirion, L.; Nougairede, A.; de Lamballerie, X. Replicative Fitness of a SARS-CoV-2 20I/501Y.V1 Variant from Lineage B.1.1.7 in Human Reconstituted Bronchial Epithelium. mBio 2021, 12, e0085021. [Google Scholar] [CrossRef]
  28. McCormick, W.; Mermel, L.A. The basic reproductive number and particle-to-plaque ratio: Comparison of these two parameters of viral infectivity. Virol. J. 2021, 18, 92. [Google Scholar] [CrossRef]
  29. Sanjuan, R. Collective properties of viral infectivity. Curr. Opin. Virol. 2018, 33, 1–6. [Google Scholar] [CrossRef]
  30. Lukassen, S.; Chua, R.L.; Trefzer, T.; Kahn, N.C.; Schneider, M.A.; Muley, T.; Winter, H.; Meister, M.; Veith, C.; Boots, A.W.; et al. SARS-CoV-2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells. EMBO J. 2020, 39, e105114. [Google Scholar] [CrossRef]
  31. Sungnak, W.; Huang, N.; Becavin, C.; Berg, M.; Queen, R.; Litvinukova, M.; Talavera-Lopez, C.; Maatz, H.; Reichart, D.; Sampaziotis, F.; et al. SARS-CoV-2 entry factors are highly expressed in nasal epithelial cells together with innate immune genes. Nat. Med. 2020, 26, 681–687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Zhu, N.; Wang, W.; Liu, Z.; Liang, C.; Wang, W.; Ye, F.; Huang, B.; Zhao, L.; Wang, H.; Zhou, W.; et al. Morphogenesis and cytopathic effect of SARS-CoV-2 infection in human airway epithelial cells. Nat. Commun 2020, 11, 3910. [Google Scholar] [CrossRef] [PubMed]
  33. Coperchini, F.; Chiovato, L.; Rotondi, M. Interleukin-6, CXCL10 and Infiltrating Macrophages in COVID-19-Related Cytokine Storm: Not One for All But All for One! Front. Immunol 2021, 12, 668507. [Google Scholar] [CrossRef] [PubMed]
  34. Smith, N.; Goncalves, P.; Charbit, B.; Grzelak, L.; Beretta, M.; Planchais, C.; Bruel, T.; Rouilly, V.; Bondet, V.; Hadjadj, J.; et al. Distinct systemic and mucosal immune responses during acute SARS-CoV-2 infection. Nat. Immunol. 2021, 22, 1428–1439. [Google Scholar] [CrossRef] [PubMed]
  35. Gamage, A.M.; Tan, K.S.; Chan, W.O.Y.; Liu, J.; Tan, C.W.; Ong, Y.K.; Thong, M.; Andiappan, A.K.; Anderson, D.E.; Wang, Y.; et al. Infection of human Nasal Epithelial Cells with SARS-CoV-2 and a 382-nt deletion isolate lacking ORF8 reveals similar viral kinetics and host transcriptional profiles. PLoS Pathog. 2020, 16, e1009130. [Google Scholar] [CrossRef]
  36. Chu, H.; Chan, J.F.; Wang, Y.; Yuen, T.T.; Chai, Y.; Hou, Y.; Shuai, H.; Yang, D.; Hu, B.; Huang, X.; et al. Comparative Replication and Immune Activation Profiles of SARS-CoV-2 and SARS-CoV in Human Lungs: An Ex Vivo Study With Implications for the Pathogenesis of COVID-19. Clin. Infect. Dis. 2020, 71, 1400–1409. [Google Scholar] [CrossRef] [Green Version]
  37. Hui, K.P.Y.; Cheung, M.C.; Perera, R.; Ng, K.C.; Bui, C.H.T.; Ho, J.C.W.; Ng, M.M.T.; Kuok, D.I.T.; Shih, K.C.; Tsao, S.W.; et al. Tropism, replication competence, and innate immune responses of the coronavirus SARS-CoV-2 in human respiratory tract and conjunctiva: An analysis in ex-vivo and in-vitro cultures. Lancet Respir. Med. 2020, 8, 687–695. [Google Scholar] [CrossRef]
  38. Ulrich, L.; Halwe, N.J.; Taddeo, A.; Ebert, N.; Schon, J.; Devisme, C.; Trueb, B.S.; Hoffmann, B.; Wider, M.; Fan, X.; et al. Enhanced fitness of SARS-CoV-2 variant of concern Alpha but not Beta. Nature 2022, 602, 307–313. [Google Scholar] [CrossRef]
Viruses 14 00951 g001 550
Figure 1. SARS-CoV-2 VoC replication kinetics in human small airway epithelial cell cultures. Viral replication of SARS-CoV-2 VoC at the viral input of 2 × 103 TCID50/insert in HsAEC cultures. The cultures were incubated at 35 °C and maintained for 10 days. (A) Percentages of infected cultures in total tested samples. Two independent experiments with HsAEC from two different donors were performed. Six independent inserts were used for each condition per experiment. (B,D,G) Viral RNA or (C,E,H) infectious particles in apical washes of infected cultures were determined by RT-qPCR or by end-point titration assay, respectively. (F,I) Relative virion infectivity (i.e., the ratio of the number infectious particles over the number of viral RNA). (J) Cell viability of HsAEC cultures. Infected samples were represented as closed markers whereas non-infected samples were indicated as open markerss. Dots represent individual samples. The dotted lines indicate the lower limit of detection (LLOD). ns p > 0.05, * p < 0.05.
Figure 1. SARS-CoV-2 VoC replication kinetics in human small airway epithelial cell cultures. Viral replication of SARS-CoV-2 VoC at the viral input of 2 × 103 TCID50/insert in HsAEC cultures. The cultures were incubated at 35 °C and maintained for 10 days. (A) Percentages of infected cultures in total tested samples. Two independent experiments with HsAEC from two different donors were performed. Six independent inserts were used for each condition per experiment. (B,D,G) Viral RNA or (C,E,H) infectious particles in apical washes of infected cultures were determined by RT-qPCR or by end-point titration assay, respectively. (F,I) Relative virion infectivity (i.e., the ratio of the number infectious particles over the number of viral RNA). (J) Cell viability of HsAEC cultures. Infected samples were represented as closed markers whereas non-infected samples were indicated as open markerss. Dots represent individual samples. The dotted lines indicate the lower limit of detection (LLOD). ns p > 0.05, * p < 0.05.
Viruses 14 00951 g001
Viruses 14 00951 g002 550
Figure 2. SARS-CoV-2 VoC replication kinetics in human bronchial airway epithelial cell cultures. Viral replication of SARS-CoV-2 VoC at the viral input of 2 × 103 TCID50/insert in HbAEC cultures. The cultures were incubated at 35 °C and maintained for 10 days. (A) Percentages of infected cultures in total tested samples. Six independent inserts were used for each infected condition while only four inserts were used for non-infected control group. (B,D,G) Viral RNA or (C,E,H) infectious particles in apical washes of infected cultures were determined by RT-qPCR or by end-point titrations, respectively. (F,I) Relative virion infectivity. (J) Cell viability of HbAEC cultures. Infected samples were represented as filled markers whereas non-infected samples were indicated as open markers. Dots represent individual samples. The dotted lines indicate the lower limit of detection (LLOD). ns p > 0.05, * p < 0.05, ** p < 0.01.
Figure 2. SARS-CoV-2 VoC replication kinetics in human bronchial airway epithelial cell cultures. Viral replication of SARS-CoV-2 VoC at the viral input of 2 × 103 TCID50/insert in HbAEC cultures. The cultures were incubated at 35 °C and maintained for 10 days. (A) Percentages of infected cultures in total tested samples. Six independent inserts were used for each infected condition while only four inserts were used for non-infected control group. (B,D,G) Viral RNA or (C,E,H) infectious particles in apical washes of infected cultures were determined by RT-qPCR or by end-point titrations, respectively. (F,I) Relative virion infectivity. (J) Cell viability of HbAEC cultures. Infected samples were represented as filled markers whereas non-infected samples were indicated as open markers. Dots represent individual samples. The dotted lines indicate the lower limit of detection (LLOD). ns p > 0.05, * p < 0.05, ** p < 0.01.
Viruses 14 00951 g002
Viruses 14 00951 g003 550
Figure 3. Cytokine and chemokine expression profiles after infection of HsAEC with different SARS-CoV-2 isolates. Production of (A) IP-10 or (B) IL-6 in HsAEC post-infection in the presence or absence of SARS-CoV-2. For this analysis, the infected samples from two experiment samples were pooled. The dotted lines represent the lower limit of quantification (LLOQ). Statistical significance between variants and uninfected control was calculated by Kruskal–Wallis with two-sided Dunn’s post hoc test. ** p < 0.01.
Figure 3. Cytokine and chemokine expression profiles after infection of HsAEC with different SARS-CoV-2 isolates. Production of (A) IP-10 or (B) IL-6 in HsAEC post-infection in the presence or absence of SARS-CoV-2. For this analysis, the infected samples from two experiment samples were pooled. The dotted lines represent the lower limit of quantification (LLOQ). Statistical significance between variants and uninfected control was calculated by Kruskal–Wallis with two-sided Dunn’s post hoc test. ** p < 0.01.
Viruses 14 00951 g003
Viruses 14 00951 g004 550
Figure 4. Cytokine and chemokine expression profiles after infection of HbAEC with different SARS-CoV-2 isolates. Production of (A) IP-10 or (B) IL-6 in HbAEC post-infection in the presence or absence of SARS-CoV-2. The dotted lines represent the lower limit of quantification (LLOQ). Statistical significance between the variants and uninfected control was calculated by Kruskal–Wallis with two-sided Dunn’s post hoc test. * p < 0.05, *** p < 0.001.
Figure 4. Cytokine and chemokine expression profiles after infection of HbAEC with different SARS-CoV-2 isolates. Production of (A) IP-10 or (B) IL-6 in HbAEC post-infection in the presence or absence of SARS-CoV-2. The dotted lines represent the lower limit of quantification (LLOQ). Statistical significance between the variants and uninfected control was calculated by Kruskal–Wallis with two-sided Dunn’s post hoc test. * p < 0.05, *** p < 0.001.
Viruses 14 00951 g004
Viruses 14 00951 g005 550
Figure 5. The replication profile of different SARS-CoV-2 isolates in nasal HAEC. The nasal tissues were infected with SARS-CoV-2 at 2 × 103 TCID50/culture at 35 °C for 7 days. (A) Viral loads from the apical washes are expressed as log10 SARS-CoV-2 RNA copies per mL. (B) Viral RNA copies, infectious particles, and relative virion infectivity at 2 d.p.i. are expressed as log10 RNA copies per mL, log10 TCID50 per mL, and log10 inverse ratio of TCID50-per-copy, respectively. (C) Cell viability of HnAEC after 7 days p.i. Each dot represents individual data. The lower limit of detection (LLOD) is indicated as dotted lines. Statistical significance was determined using Kruskal–Wallis with Dunnett’s multiple comparison test with ns p > 0.05, * p < 0.05, ** p < 0.01, **** p < 0.0001.
Figure 5. The replication profile of different SARS-CoV-2 isolates in nasal HAEC. The nasal tissues were infected with SARS-CoV-2 at 2 × 103 TCID50/culture at 35 °C for 7 days. (A) Viral loads from the apical washes are expressed as log10 SARS-CoV-2 RNA copies per mL. (B) Viral RNA copies, infectious particles, and relative virion infectivity at 2 d.p.i. are expressed as log10 RNA copies per mL, log10 TCID50 per mL, and log10 inverse ratio of TCID50-per-copy, respectively. (C) Cell viability of HnAEC after 7 days p.i. Each dot represents individual data. The lower limit of detection (LLOD) is indicated as dotted lines. Statistical significance was determined using Kruskal–Wallis with Dunnett’s multiple comparison test with ns p > 0.05, * p < 0.05, ** p < 0.01, **** p < 0.0001.
Viruses 14 00951 g005
Viruses 14 00951 g006 550
Figure 6. The expression of COVID-19-related chemokine and cytokine in HnAEC. Concentrations of IP-10 (A,B) and IL-6 (C,D) in the presence or absence of SARS-CoV-2 are expressed as pg per mL. Mean values with 95% CI are shown and statistical significance between the variants and non-infected control was calculated by Kruskal–Wallis with two-sided Dunn’s post hoc test. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.
Figure 6. The expression of COVID-19-related chemokine and cytokine in HnAEC. Concentrations of IP-10 (A,B) and IL-6 (C,D) in the presence or absence of SARS-CoV-2 are expressed as pg per mL. Mean values with 95% CI are shown and statistical significance between the variants and non-infected control was calculated by Kruskal–Wallis with two-sided Dunn’s post hoc test. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.
Viruses 14 00951 g006
Viruses 14 00951 g007 550
Figure 7. Comparison of SARS-CoV-2 virion infectivity in different human respiratory tissues. For this analysis, infected samples on day 2 p.i. from independent experiments were pooled. The relative infectivity is expressed as log10 inverse ratio of TCID50-per-copy. A comparison between the different cell types (A) as well as a comparison between the different VoCs (B) is shown. Each dot represents individual cultures and statistical significance was determined by two-way ANOVA with a Turkey post-hoc test. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.
Figure 7. Comparison of SARS-CoV-2 virion infectivity in different human respiratory tissues. For this analysis, infected samples on day 2 p.i. from independent experiments were pooled. The relative infectivity is expressed as log10 inverse ratio of TCID50-per-copy. A comparison between the different cell types (A) as well as a comparison between the different VoCs (B) is shown. Each dot represents individual cultures and statistical significance was determined by two-way ANOVA with a Turkey post-hoc test. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.
Viruses 14 00951 g007
Table
Table 1. Descriptive summary of log10 inverse ratio of TCID50-per-vRNA copy.
Table 1. Descriptive summary of log10 inverse ratio of TCID50-per-vRNA copy.
Insert OriginSARS-CoV-2 StrainMeanStandard DeviationN
Nasal insertsUSA-WA13.140.639
BavPat12.630.539
B.1.1.73.270.449
B.1.3512.300.349
B.1.617.23.100.419
B.1.1.5293.710.649
Bronchial insertsUSA-WA10.980.363
BavPat10.710.343
B.1.1.71.560.566
B.1.3511.380.486
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Do, T.N.D.; Claes, S.; Schols, D.; Neyts, J.; Jochmans, D. SARS-CoV-2 Virion Infectivity and Cytokine Production in Primary Human Airway Epithelial Cells. Viruses 2022, 14, 951. https://doi.org/10.3390/v14050951

AMA Style

Do TND, Claes S, Schols D, Neyts J, Jochmans D. SARS-CoV-2 Virion Infectivity and Cytokine Production in Primary Human Airway Epithelial Cells. Viruses. 2022; 14(5):951. https://doi.org/10.3390/v14050951

Chicago/Turabian Style

Do, Thuc Nguyen Dan, Sandra Claes, Dominique Schols, Johan Neyts, and Dirk Jochmans. 2022. "SARS-CoV-2 Virion Infectivity and Cytokine Production in Primary Human Airway Epithelial Cells" Viruses 14, no. 5: 951. https://doi.org/10.3390/v14050951

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop