Abstract
Objectives
To determine how the intrinsic severity of successively dominant SARS-CoV-2 variants changed over the course of the pandemic.Methods
A retrospective cohort analysis in the NHS Greater Glasgow and Clyde (NHS GGC) Health Board. All sequenced non-nosocomial adult COVID-19 cases in NHS GGC with relevant SARS-CoV-2 lineages (B.1.177/Alpha, Alpha/Delta, AY.4.2 Delta/non-AY.4.2 Delta, non-AY.4.2 Delta/Omicron, and BA.1 Omicron/BA.2 Omicron) during analysis periods were included. Outcome measures were hospital admission, ICU admission, or death within 28 days of positive COVID-19 test. We report the cumulative odds ratio; the ratio of the odds that an individual experiences a severity event of a given level vs all lower severity levels for the resident and the replacement variant after adjustment.Results
After adjustment for covariates, the cumulative odds ratio was 1.51 (95% CI: 1.08-2.11) for Alpha versus B.1.177, 2.09 (95% CI: 1.42-3.08) for Delta versus Alpha, 0.99 (95% CI: 0.76-1.27) for AY.4.2 Delta versus non-AY.4.2 Delta, 0.49 (95% CI: 0.22-1.06) for Omicron versus non-AY.4.2 Delta, and 0.86 (95% CI: 0.68-1.09) for BA.2 Omicron versus BA.1 Omicron.Conclusions
The direction of change in intrinsic severity between successively emerging SARS-CoV-2 variants was inconsistent, reminding us that the intrinsic severity of future SARS-CoV-2 variants remains uncertain.Free full text
![Logo of pheelsevier](/http/europepmc.org/corehtml/pmc/pmcgifs/logo-pheelsevier.png)
Directions of change in intrinsic case severity across successive SARS-CoV-2 variant waves have been inconsistent
Associated Data
Abstract
Objectives
To determine how the intrinsic severity of successively dominant SARS-CoV-2 variants changed over the course of the pandemic.
Methods
A retrospective cohort analysis in the NHS Greater Glasgow and Clyde (NHS GGC) Health Board. All sequenced non-nosocomial adult COVID-19 cases in NHS GGC with relevant SARS-CoV-2 lineages (B.1.177/Alpha, Alpha/Delta, AY.4.2 Delta/non-AY.4.2 Delta, non-AY.4.2 Delta/Omicron, and BA.1 Omicron/BA.2 Omicron) during analysis periods were included. Outcome measures were hospital admission, ICU admission, or death within 28 days of positive COVID-19 test. We report the cumulative odds ratio; the ratio of the odds that an individual experiences a severity event of a given level vs all lower severity levels for the resident and the replacement variant after adjustment.
Results
After adjustment for covariates, the cumulative odds ratio was 1.51 (95% CI: 1.08–2.11) for Alpha versus B.1.177, 2.09 (95% CI: 1.42–3.08) for Delta versus Alpha, 0.99 (95% CI: 0.76–1.27) for AY.4.2 Delta versus non-AY.4.2 Delta, 0.49 (95% CI: 0.22–1.06) for Omicron versus non-AY.4.2 Delta, and 0.86 (95% CI: 0.68–1.09) for BA.2 Omicron versus BA.1 Omicron.
Conclusions
The direction of change in intrinsic severity between successively emerging SARS-CoV-2 variants was inconsistent, reminding us that the intrinsic severity of future SARS-CoV-2 variants remains uncertain.
Introduction
Since the SARS-CoV-2 pandemic started in late 2019, a succession of variants have achieved dominance, each replacing the previous dominant variant. From late 2020 these were designated variants of concern (VOCs); variants that exhibit increased transmission rates, antigenic differences, and/or case severity.1 The three VOCs that most impacted both the pandemic and epidemic in Scotland were Alpha (Pango lineage B.1.1.7), which emerged in the UK in September 2020,2 Delta (Pango lineage B.1.617.2), which emerged in India prior to October 20201 and spread globally in May 2021, with> 1000 introductions to the UK,3 and most recently Omicron (Pango lineage B.1.1.529), which emerged in Africa in November 2021,1 and spread very rapidly around the globe. Before the Omicron variant emerged, the Delta sublineage AY.4.2 was on course to replace the other Delta lineages, with growth rate estimates4 implying that it would become dominant in the UK in early 2022. This spread was arrested by the emergence of the more transmissible and immune-evading Omicron variant, which supplanted nearly all non-Omicron diversity5 and has itself diversified into several subvariants.
Understanding any change in disease severity associated with infection by new variants of a virus (especially one that has newly entered the human population) is critical from a clinical, public health, and basic science perspective. For example, knowledge of the severity of a new variant is vitally important in the decision-making process regarding the stringency of control measures and the roll out of vaccination and other treatments. It is expected that upon entry to a new host species from a zoonotic reservoir, the consequences of infection will be unpredictable with the severity of disease caused by the pathogen likely to be far from its evolutionary optima.6 SARS-CoV-2 is on average associated with low virulence in younger age groups whilst severe outcomes manifest in older age groups and those with comorbidities. Evolution of virulence may therefore not be strongly constrained, rather, factors governing transmission and immune evasion are likely to determine the direction of evolution.7
Our aim was to determine how the intrinsic severity of successively dominant SARS-CoV-2 variants has changed over the course of the Coronavirus disease 2019 (COVID-19) pandemic in Scotland, by making use of data from the largest health board in the country, Greater Glasgow and Clyde. To explore this, we used detailed clinical metadata and viral genomic data from National Health Service Greater Glasgow and Clyde (NHS GGC) to analyse relative case severity within 28 days of diagnosis between successive dominant lineages ( Fig. 1); B.1.177 versus Alpha, Alpha versus Delta, non-AY.4.2 versus AY.4.2 Delta, non-AY.4.2 Delta versus Omicron, and BA.1 Omicron versus BA.2 Omicron. This series of comparisons allowed us to assess trends in severity, removing this bias and accounting for other critical variables including detailed comorbidity and changes in the availability of new treatments and vaccination.
Patients and methods
This is a cohort study comparing the clinical outcomes of adults infected with a resident dominant SARS-CoV-2 lineage, with those infected with the SARS-CoV-2 lineage that would subsequently displace it.
Genome sequencing
Sequences were generated using the ARTIC Network protocol, originally developed for Oxford Nanopore-based sequencing.8 and derived versions adapted to Illumina and ARTIC-unrelated amplicon-based protocols. The COVID-19 Genomics UK Consortium (COG-UK) pipeline was used for alignment and Pango lineage assignment.9
Data inclusion/exclusion
This work was originally performed in real time for the COVID-19 response, as such, the time windows begin at the start of the first month where we believed we had enough sequences of both the resident and replacement lineages to perform the analysis, and end at the end of last month meeting the same criteria.
For the B.1.177/Alpha analysis, we included all sequenced samples with full data available on all adjustment variables from within the NHS GGC health board between 1st November 2020 and 30th January 2021 (B.1.177: n=
807; Alpha: n
=
833). A full demographic breakdown of samples is shown in Table S1. All sequences assigned as B.1.177 and associated sublineages were merged into a single category for the analysis.
For the Alpha/Delta analysis we included all sequenced samples with full metadata between 1st April 2021 and 30th June 2021 (Alpha: n=
2104; Delta: n
=
3448). All sequences assigned into B.1.617.2 and associated sublineages were merged into the Delta category for the analysis. A full demographic breakdown of samples is shown in Table S2.
For the Delta/AY.4.2 analysis we included all sequenced samples with full metadata between 1st July 2021 and 31st October 2021 (non-AY.4.2 Delta: n=
8644; AY.4.2 Delta: n
=
969). The AY.4.2 category was defined as all sequences assigned as either AY.4.2 or sublineages thereof, all remaining Delta sublineages were combined into the comparison category. A full demographic breakdown of samples is shown in Table S3.
For the non-AY.4.2 Delta/Omicron analysis, we included all sequenced samples with full metadata between 1st December 2021 and 31st December 2021 (non-AY.4.2 Delta: 1164; Omicron: 2694). The Omicron category was defined as all sequences assigned as BA.1 or B.1.1.529, the Delta category was defined as for the Delta/AY.4.2 comparison. Full demographic breakdown of samples is shown in Table S4.
For the BA.1 Omicron/BA.2 Omicron analysis, we included all sequenced samples with full metadata between 1st January 2022 and 31st March 2022 (BA.1 Omicron: 9949; BA.2 Omicron: 12994). The BA.1 Omicron category was defined as all sequences assigned as BA.1 or sublineages thereof, the BA.2 Omicron category was defined as all sequences assigned as BA.2 or sublineages thereof. Full demographic breakdown of samples is shown in Table S5.
Cases with hospital acquired COVID-19, defined as a first positive polymerase chain reaction (PCR) test occurring more than 48h following admission to hospital, and all cases younger than 18 were excluded.
Data linkage
Cohorts and de-identified linked data for the entire NHS GGC health board (1.2 million people) were prepared by the West of Scotland Safe Haven at NHS GGC. Data used in the analysis included maximum clinical severity at 28 days after the first positive test via a 4-point ordinal scale (1. No hospitalisation; 2. Hospitalisation (excluding elective surgery); 3. Admission to High Dependency Unit (HDU)/Intensive Care Unit (ICU); 4. Death), age at diagnosis, date of positive test, sex, partial postcode, number of vaccine doses, number of relevant comorbidities or risks of ill health (chronic cardiac disease, chronic respiratory disease, chronic renal disease, liver disease, dementia, chronic neurological conditions, connective tissue disease, diabetes, HIV infection, malignant tumours, clinician defined obesity, case shielding, immunosuppressive drugs, chemotherapy) and reinfection. Vaccination, comorbidities, reinfection, age and biological sex at birth were included as they are known to impact clinical outcomes. Note that, in Scotland, vaccination schedules prioritised the elderly and vulnerable. Therefore, throughout, the number of vaccine doses is correlated with age with older individuals having had more doses on average than younger individuals. Date of positive test was included to attempt to adjust for any time varying effects across the study window (e.g., changing clinical practice, hospital load, etc.). Partial postcode was included to attempt to account for spatial clustering of risk mediated though unmeasured variables (e.g., deprivation, ethnicity, etc.).
Severity was scored twice. For severity for “with” analyses, cases were assigned the most severe event that occurred within 28 days of their positive test. For the “of” analyses, events were only counted when they were explicitly linked to COVID-19 infection in the electronic patient records. Different datasets used in the linkage recorded “cause of” event differently. In datasets where ICD-10 codes were used (SMR01, accident and emergency, and deaths), a COVID-19 related ICD-10 code was required (specifically, any code starting U07, U04.9 which corresponds to an incorrect usage of the SARS ICD-10 code, and U10). In datasets where ICD-10 codes were not used (Scottish Intensive Care Society Audit Group data), the string “covid” was searched for case-insensitively in the free text entry. When the number of vaccine doses an individual had been given was calculated, if the last dose had been received less than 14 days before the date of the positive PCR test, it was ignored. Individuals with multiple confirmed episodes of infection (defined as separate positive PCR results more than 90 days apart) were marked as reinfected for any episodes after the first.
Statistical analysis of clinical data
The four-level patient outcome data were analysed using cumulative generalised additive mixed models (GAMMs) with logit links10 fit using Bayesian inference. These GAMMs included lineage, reinfection, patient sex and number of vaccine doses as categorical fixed effects and number of 4C Mortality Score comorbidities identified as a continuous fixed effect, with partial postcode included as a random effect. We included age and date of positive test as non-linear penalised regression splines. The basis dimension of the penalised regression splines was set to the number of unique dates of positive tests minus one and the number of unique ages (rounded to year) minus one respectively, with the intention that regularisation occur through the prior. Given that the pandemic was in its early stages during the first comparison (B.1.177/Alpha) and the vaccination campaign had not yet started, both reinfection and number vaccine doses received were excluded from the first model.
The same classes of parameters received the same priors in each model. The intercepts of the models were given t-distribution (location = 0, scale = 2.5, df = 3) priors, fixed effects were given normal (mean = 0, standard deviation = 2.5) priors, random effects and spline standard deviations were given exponential (mean = 2.5) priors. All severity models were fitted using the brms (v. 2.14.4) R package.11 All presented models had no divergent transitions and effective sample sizes of over 200 for all parameters.
We tested the robustness of these estimates to epidemic phase bias. Epidemic phase bias is caused by patient outcomes being correlated with the time from infection to positive test. This correlation makes the estimated odds ratio adjusting for time of positive test a biased estimator of the odds ratio adjusting for time of infection, which is the desired target of estimation.12 When lineages differ in their incidence, this effect can lead to incorrectly concluding that one lineage is associated with more severe disease than the other, when the estimated difference is driven almost entirely by the epidemic phase bias. Sensitivity to this bias was assessed using the method of Seaman et al. 12 We added four days to the population who experienced more extreme outcomes (i.e., hospitalisation, admission to ICU/HCU or death), generating a modified test time for each individual, where if the patient was not hospitalised, it was their original test time and otherwise it was this new test time. We refit the model using the modified times and using individuals whose modified times lie within the inclusion window. The resulting estimate of the cumulate odds ratio was then compared with the cumulate odds ratio estimated from the original model.
Ethical approval
Anonymised clinical data were available using the (NHS GGC SafeHaven platform and included vaccination status (dates and product names for each dose), demographic data (age and sex) comorbidities, hospital records and dates of positive and negative PCR tests, for 1.2 million inhabitants of the (NHS GGC) area over 18 years of age. Data were matched by CHI number and pseudonymised before analysis. Derogated ethical approval for the study was granted by the NHS GGC SafeHaven committee (GSH/21/IM/001).
Results
B.1.177/Alpha
Our first comparison was the severity of the SARS-CoV-2 Alpha VOC versus the previous dominant non-VOC lineage, Pango designation B.1.177 (B.1.177: n=
807; Alpha: n
=
833). The replacement of the lineages over time can be seen in Fig. 1. We found that confirmed Alpha cases were associated with more severe infection (“of” - median cumulative odds ratio: 1.60; 95% central interval: 1.10–2.30; probability of association with increased severity:>0.99; “with” – median cumulative odds ratio: 1.51; 95% central interval: 1.08–2.11; probability of association with increased severity: 0.99; breakdown of number events in each group Table S1). The breakdown of “with” severity score by age and lineage can be seen in
Fig. 2a. Parameter estimates for all parameters can be found in Table S6 and Fig. S1. We found that, in our epidemic phase bias sensitivity analysis, Alpha remained associated with increased severity, but that the magnitude was reduced, as would be expected, and that in both analyses the probability of the effect being positive was reduced below 95% (“of” - median cumulative odds ratio: 1.31; 95% central interval: 0.91–1.88; probability of association with increased severity: 0.93; “with” – median cumulative odds ratio: 1.24; 95% central interval: 0.88–1.79; probability of association with increased severity: 0.88).
![Fig. 2 Fig. 2](/http/europepmc.org/articles/PMC10234362/bin/gr2_lrg.jpg)
Comparison of case severity by age between pairs of co-circulating SARS-CoV-2 lineages Clinical severity was measured on a four-level ordinal scale based on observed outcomes within 28 days of a positive test: no hospitalisation, hospitalisation, admission to Intensive Care Unit/High Dependency Unit (ICU/HDU), death. The first column compares B.1.177 and Alpha between 1st November 2020 and 30th January 2021. The second column compares Alpha and Delta between 1st April 2021 and 30th June 2021. The third column compares Delta lineages against AY.4.2 between 1st July 2021 and 30th September 2021. The fourth column compares non-AY.4.2 Delta and Omicron between 1st December 2021 and 31st December 2021. The final column compares BA.1 Omicron and BA.2 Omicron between 1st January 2022 and 31st March 2022.
Alpha/Delta
Our second comparison was the severity of the SARS-CoV-2 Delta VOC versus the previous dominant Alpha VOC lineage (Alpha: n=
2104; Delta: n
=
3448). We estimate a substantial increase in case severity associated with Delta infections relative to Alpha (“of” - median cumulative odds ratio: 2.19; 95% central interval: 1.48–3.32; probability of association with increased severity:>0.99; “with” - median cumulative odds ratio: 2.09; 95% central interval: 1.42–3.08; probability of association with increased severity:>0.99; breakdown of number events in each group Table S2). The breakdown of “with” severity score by age and lineage can be seen in Fig. 2b. Parameter estimates for all parameters can be found in Table S7 and Fig. S1. This effect is large enough that in epidemic phase bias sensitivity bias model, we still estimate a positive effect with a probability of positivity of over 0.95 (“of” - median cumulative odds ratio: 1.45; 95% central interval: 0.98–2.16; probability of association with increased severity: 0.96; “with” - median cumulative odds ratio: 1.40; 95% central interval: 0.97–2.02; probability of association with increased severity: 0.97).
Non-AY.4.2 Delta/AY.4.2 Delta
Our third comparison was the severity of the SARS-CoV-2 Delta sublineage AY.4.2 versus the other Delta sublineages it was replacing (non-AY.4.2 Delta: n=
8644; AY.4.2 Delta: n
=
969). We estimate that the AY.4.2 lineage infections are associated with similar case severity to that seen in other Delta sublineage infections (“of” - median cumulative odds ratio: 1.05; 95% central interval: 0.79–1.39; probability of association with increased severity: 0.64; “with” - median cumulative odds ratio: 0.99; 95% central interval: 0.76–1.27; probability of association with increased severity: 0.46; breakdown of number events in each group Table S3). The breakdown of “with” severity score by age and lineage can be seen in Fig. 2c. Parameter estimates for all parameters can be found Table S8 and Fig. S1. In this case, there is no noticeable effect of epidemic phase bias, likely as the growth rate difference between the two variants was small (“of” - median cumulative odds ratio: 1.02; 95% central interval: 0.76–1.35; probability of association with increased severity: 0.57; “with” - median cumulative odds ratio: 0.98; 95% central interval: 0.75–1.27; probability of association with increased severity: 0.44).
Non-AY.4.2 Delta/Omicron
Our fourth comparison was the severity of the SARS-CoV-2 Omicron VOC versus the non-AY.4.2 Delta sublineages (non-AY.4.2 Delta: 1164; Omicron: 2694). We find that Omicron (BA.1 sublineage) infection is associated with substantially less severe disease (“of” - median cumulative odds ratio: 0.15; 95% central interval: 0.01–1.48; probability of association with increased severity: 0.06; “with” - median cumulative odds ratio: 0.49; 95% central interval: 0.22–1.06; probability of association with increased severity: 0.04; breakdown of number events in each group Table S4). The breakdown of “with” severity score by age and lineage can be seen in Fig. 2d. Parameter estimates for all parameters can be found in Table S9 and Fig. S1. Both the “with” and “of” analyses show a strong impact of epidemic phase bias. Omicron was the faster growing lineage, so the estimates of Omicron are driven more negative (“of” - median cumulative odds ratio: 0.07; 95% central interval:<0.01–0.68; probability of association with increased severity: 0.01; “with” - median cumulative odds ratio: 0.19; 95% central interval: 0.08–0.42; probability of association with increased severity:<0.01).
BA.1 Omicron/BA.2 Omicron
Our final comparison was the severity of the BA.1 Omicron sublineage versus the BA.2 Omicron sublineage (BA.1 Omicron: 9949; BA.2 Omicron: 12994). We find that BA.2 Omicron infection is associated with less severe disease than BA.1 Omicron (“of” - median cumulative odds ratio: 0.83; 95% central interval: 0.63–1.12; probability of association with increased severity: 0.11; “with” - median cumulative odds ratio: 0.86; 95% central interval: 0.68–1.09; probability of association with increased severity: 0.11; breakdown of number events in each group Table S5). The breakdown of “with” severity score by age and lineage can be seen in Fig. 2e. Parameter estimates for all parameters can be found in Table S10 and Fig. S1. This association is uncertain, but both the “with” and “of” analyses show a strong impact of epidemic phase bias. As BA.2 Omicron was the faster growing lineage, so the estimates of BA.2 Omicron are driven more negative, and our certainty of an association with a reduction of severity becomes higher (“of” - median cumulative odds ratio: 0.51; 95% central interval: 0.38–0.69; probability of association with increased severity:<0.01; “with” - median cumulative odds ratio: 0.57; 95% central interval: 0.49–0.72; probability of association with increased severity:<0.01).
Discussion
The principal findings of this study of the relative severity of COVID-19 cases caused by successive SARS-CoV-2 variant waves in Scotland was that Alpha was associated with more severe disease than B.1.177, Delta was associated with more severe disease than Alpha, non-AY.4.2 Delta and AY.4.2 Delta were associated with similar disease severity, Omicron was associated with much less severe disease than non-AY.4.2 Delta, and that BA.2 Omicron was associated with less severe disease than BA.1 Omicron. These conclusions were after accounting for comorbidities, date of positive test (to adjust for time varying effects across the study window), previous infection and vaccine availability, and robust to epidemic phase bias and the possibility of coincidental SARS-CoV-2 infection at admission to hospital. The successive replacements that we studied were not consistent in the direction of change in case severity.
Our study design has several strengths. It is the first study to our knowledge to analyse the sequential replacement of variants throughout the pandemic with respect to the progression of severity attributable to virus evolution, and to use a consistent analytical approach across sequential SARS-CoV-2 lineages. Lauring et al.13 is, in spirit, similar to our work, but their primary focus is not on the trajectory of severity, and they do not consider multiple time matched comparisons.
Secondly, we test the robustness of the severity analysis to epidemic phase bias and the impact of differences in the definition of severity outcomes (“of” analysis versus “with” analysis) across sequential variants, showing that these different definitions lead to qualitatively similar results.
Also, our approach takes a broad view of the definition of disease severity, by including community- and hospital-based cases and by considering a wider variety of clinical outcomes than those considered in most previous severity analyses. This is both a strength and limitation of our study as it provides a more accurate assessment of disease severity in the adult population than studies that have concentrated solely on hospitalised patients, or on the risk of hospitalisation, however due to our sample size the numbers of deaths and admissions to ICU in our cohort were low.
Our study does, however, have some limitations. Firstly, we only include cases with sequenced genomes, and thus our sample is biased towards cases with lower cycle threshold (Ct), because these cases are more likely to have been sequenced. This is likely to be particularly important when the Ct distribution of infections differs between variants, such as the higher Ct values seen with Omicron infection compared to Delta infection.14, 15 The direction of the sampling bias on our results depends on which of the variants is associated with the lower Ct infections. Requiring sequences for our samples also limits us to the set of individuals who have been tested by PCR, likely to represent hospitalised patients more than those in the community and assumes that this these proportions are stable across our sub-periods, an assumption which is untestable from the data we have available. We did not have access to biomarker data, and as such could not adjust for correlates of immunity (e.g., T cell counts, antibody titre) directly. Additionally, our sample size was not large enough to adjust for all the factors we would have liked to (e.g., fitting a dose by vaccine brand interaction, given that different brands are known to provide differential protection against different variants,7, 16, 17 or accounting for the differential protection induced by infection prior to vaccination).
There have been a series of other studies investigating each of the comparisons in our study individually. For Alpha versus previous variants, most previous studies in most regions have also estimated an increase in severity over extant diversity. A wide variety of end points have been used, 28-day mortality, hospitalisation, and an ordinal scale based around supplemental oxygen.18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30 Our estimates are consistent with the majority of these studies.
When considering the Delta variant, two UK community analyses found that Delta (or an S-gene proxy) infections were associated with a higher risk of admission to hospital than with Alpha.31, 32 Comparable results were observed in Danish, US, and Canadian populations with a study in a Norwegian population being the exception.33, 34, 35, 36 The US and Canadian studies also found that Delta was associated with increased risk of ICU admission and death.34, 35 Most of these studies are therefore consistent with our results.
Our estimate that AY.4.2 is associated with approximately the sample severity as other Delta sublineages is inconsistent with an analysis of the English population which found that confirmed AY.4.2 cases are associated with lower hospitalisation risk than cases associated with non-AY.4.2 Delta.37 This inconsistency between our study and others may be explained by differences in the adjustment variables used, or because the larger sample size in Nyberg et al. 37 allowed precise isolation of a small negative effect, with their effect estimate falling within our credible interval.
Our results are consistent with studies from England, Scotland, Denmark, Norway, South Africa, Canada, and the US suggesting that Omicron infections are less severe than infections with Delta.13, 38, 39, 40, 41, 42, 43, 44 This reduction in case severity resulted in increased numbers of patients being admitted to hospital with a coincidental positive SARS-CoV-2 test rather than due to COVID-19, but this did not seem to overly impact our estimate of the severity of the variant relative to Delta. Contrary to some reporting,45, 46 we find that BA.2 is infection is likely associated with a small further decrease in severity relative to BA.1 infection, potentially because we had access to genomically confirmed infections and did not have to use proxies of lineage.
It is important to appreciate measures of disease severity are highly context dependent. The clinical situation of the pandemic has shifted dramatically in Scotland during the study period, from a time with very little prior immunity to one with widespread vaccine and prior infection mediated immunity. Treatment availability with steroids, antivirals and antithrombotic agents has had a huge impact on reducing length of hospital stay and mortality.47, 48, 49 Testing patterns have also changed dramatically across the study, with periods of higher and lower rates of testing, as well as changes in the proportion of positives being sequenced. All these factors may impact the relative severity of variants, as well as the meanings of estimates of the association between lineage and severity, as the populations being compared are not the same between comparisons. For this reason, our results cannot be used to compare the intrinsic case severity of variants that did not co-circulate.
Our results demonstrate that successive variants of SARS-CoV-2 are associated with inconsistent directions of change in intrinsic disease severity after other factors are accounted for in the analysis, including comorbidities, vaccination, previous infection, and date of positive test. In keeping with this finding, we have recently demonstrated fundamental changes in the life cycle of the Omicron variants BA.1 and BA.2 that may provide a biological explanation for the substantial drop in severity associated with this variant7; Omicron is associated with less cell-to-cell fusion and tropism for nasal epithelial cells rather than cells present in the lungs as well as an endosomal rather than a direct cell entry pathway.
Given that the progression of SARS-CoV-2 virulence has not been consistent over time and that the entry pathway of SARS-CoV-2 may be changed by single amino acid changes within the S2 domain of the spike protein,50 historical trends in severity should not be used to predict the severity of future variants. This is especially true as there is no guarantee a newly emerging lineage will descend from the current Omicron diversity, and thus may have very different biological characteristics. However, once a variant has emerged, the likelihood of immune evasion and the method of cell entry may be estimated from the genome sequence. The relative reduction in severity seen with the Omicron variant should not make us complacent to the potential risks of future SARS-CoV-2 variants. Any increase to intrinsic severity in a variant with similar transmissibility to Omicron could be devastating to health systems and communities.
Results from this study were made available in a timely manner to national government agencies, and informed decisions on public health restrictions, vaccination strategies, and healthcare service planning. This study provides an important baseline for future research monitoring the relative severity of new variants and highlights the importance of ongoing genomic “early-warning” surveillance to detect new variants of concern promptly as they emerge.
Declaration of Competing Interest
The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Thomas C Williams is Principal Investigator for the BronchStart project, which is funded by the Respiratory Syncytial Virus Consortium in Europe (RESCEU), with data collection supported by the National Institute for Health Research.
Acknowledgements
The authors would like to acknowledge the work of the West of Scotland Safe Haven team in supporting extractions and linkage to de-identified NHS patient datasets. COG-UK is supported by funding from the Medical Research Council (MRC) part of UK Research & Innovation (UKRI), the National Institute of Health Research (NIHR) and Genome Research Limited, operating as the Wellcome Sanger Institute.
This work was supported by UKRI through the JUNIPER consortium (grant number MR/V038613/1, https://maths.org/juniper/); the Medical Research Council (MRC) (grant numbers MC UU 12014/12, MC UU 00002/11, https://www.ukri.org/councils/mrc/); the G2P-UK National Virology Consortium (grant number MR/W005611/1, https://www.ukri.org). The funders had no input on the study design and analysis.
Footnotes
Appendix ASupplementary data associated with this article can be found in the online version at 10.1016/j.jinf.2023.05.019.
Appendix A. Supplementary material
Supplementary material
.
Supplementary material
.
References
![[left angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/lang.gif)
![[right angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/rang.gif)
![[left angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/lang.gif)
![[right angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/rang.gif)
![[left angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/lang.gif)
![[right angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/rang.gif)
![[left angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/lang.gif)
![[right angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/rang.gif)
![[left angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/lang.gif)
![[right angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/rang.gif)
![[left angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/lang.gif)
![[right angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/rang.gif)
![[left angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/lang.gif)
![[right angle bracket]](/http/europepmc.org/corehtml/pmc/pmcents/rang.gif)
Full text links
Read article at publisher's site: https://doi.org/10.1016/j.jinf.2023.05.019
Read article for free, from open access legal sources, via Unpaywall:
http://www.journalofinfection.com/article/S0163445323002979/pdf
Citations & impact
Impact metrics
Citations of article over time
Alternative metrics
![Altmetric item for https://www.altmetric.com/details/149420993](https://api.altmetric.com/v1/donut/149420993_64.png)
Discover the attention surrounding your research
https://www.altmetric.com/details/149420993
Smart citations by scite.ai
Explore citation contexts and check if this article has been
supported or disputed.
https://scite.ai/reports/10.1016/j.jinf.2023.05.019
Article citations
Phenotyping the virulence of SARS-CoV-2 variants in hamsters by digital pathology and machine learning.
PLoS Pathog, 19(11):e1011589, 07 Nov 2023
Cited by: 1 article | PMID: 37934791 | PMCID: PMC10656012
Vaccination, immunity, and the changing impact of COVID-19 on infant health.
Proc Natl Acad Sci U S A, 120(49):e2311573120, 27 Nov 2023
Cited by: 5 articles | PMID: 38011548 | PMCID: PMC10710065
Key Considerations during the Transition from the Acute Phase of the COVID-19 Pandemic: A Narrative Review.
Vaccines (Basel), 11(9):1502, 19 Sep 2023
Cited by: 8 articles | PMID: 37766178 | PMCID: PMC10537111
Review Free full text in Europe PMC
Unraveling the Dynamics of Omicron (BA.1, BA.2, and BA.5) Waves and Emergence of the Deltacton Variant: Genomic Epidemiology of the SARS-CoV-2 Epidemic in Cyprus (Oct 2021-Oct 2022).
Viruses, 15(9):1933, 15 Sep 2023
Cited by: 2 articles | PMID: 37766339 | PMCID: PMC10535466
Change in the Clinical Picture of Hospitalized Patients with COVID-19 between the Early and Late Period of Dominance of the Omicron SARS-CoV-2 Variant.
J Clin Med, 12(17):5572, 26 Aug 2023
Cited by: 4 articles | PMID: 37685639 | PMCID: PMC10488127
Go to all (7) article citations
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Replacement dynamics and the pathogenesis of the Alpha, Delta and Omicron variants of SARS-CoV-2.
Epidemiol Infect, 151:e32, 20 Dec 2022
Cited by: 2 articles | PMID: 36535802 | PMCID: PMC9990386
Estimates of SARS-CoV-2 Omicron BA.2 Subvariant Severity in New England.
JAMA Netw Open, 5(10):e2238354, 03 Oct 2022
Cited by: 35 articles | PMID: 36282501 | PMCID: PMC9597387
Comparative analysis of the risks of hospitalisation and death associated with SARS-CoV-2 omicron (B.1.1.529) and delta (B.1.617.2) variants in England: a cohort study.
Lancet, 399(10332):1303-1312, 16 Mar 2022
Cited by: 677 articles | PMID: 35305296 | PMCID: PMC8926413
Clinical outcomes of the severe acute respiratory syndrome coronavirus 2 Omicron and Delta variant: systematic review and meta-analysis of 33 studies covering 6 037 144 coronavirus disease 2019-positive patients.
Clin Microbiol Infect, 29(7):835-844, 18 Mar 2023
Cited by: 22 articles | PMID: 36934872 | PMCID: PMC10023211
Review Free full text in Europe PMC
Rapid review and meta-analysis of serial intervals for SARS-CoV-2 Delta and Omicron variants.
BMC Infect Dis, 23(1):429, 26 Jun 2023
Cited by: 4 articles | PMID: 37365505 | PMCID: PMC10291789
Review Free full text in Europe PMC
Funding
Funders who supported this work.
Medical Research Council (10)
COVID-19 Modelling Consortium: quantitative epidemiological predictions in response to an evolving pandemic
Professor Matthew Keeling, University of Warwick
Grant ID: MR/V038613/1
Viral Genomics and Bioinformatics
Professor Andrew Davison, University of Glasgow
Grant ID: MC_UU_12014/12
Grant ID: MC UU 00002/11
Grant ID: HDRUK2021.0155
COVID-19 Genomics UK Consortium (COG-UK)
Professor Sharon Peacock, Public Health England
Grant ID: MC_PC_19027
Integrative viral genomics and bioinformatics platform
Professor David Robertson, University of Glasgow
Grant ID: MC_UU_00034/5
Preparedness Platform
Dr Emma Thomson, University of Glasgow
Grant ID: MC_UU_00034/6
G2P-UK; A National Virology Consortium to address phenotypic consequences of SARSCoV-2 genomic variation
Professor Wendy Barclay, Imperial College London
Grant ID: MR/W005611/1
Grant ID: MC UU 12014/12
Evidence Synthesis to inform health related decision making
Professor Daniela De Angelis, University of Cambridge
Grant ID: MC_UU_00002/11