key: cord-0872175-pq3n18ae authors: Pascall, D. J.; Vink, E.; Mollett, G.; Blacow, R.; Bulteel, N.; Campbell, R.; Campbell, A.; Clifford, S.; Davis, C.; da Silva Filipe, A.; Fjodorova, L.; Forrest, R.; Goldstein, E.; Gunson, R.; Haughney, J.; Holden, M. T.; Honour, P.; Hughes, J.; James, E.; Lewis, T.; McHugh, M.; Onishi, Y.; Parcell, B.; Robertson, D. L.; El Sakka, N.; Shabaan, S.; Shepherd, J. G.; Smollett, K.; Templeton, K.; Wastnedge, E.; Williams, T. C.; Nyberg, T.; Seaman, S. R.; MacLean, O.; Wilkie, C.; Ray, S.; The COVID-19 Genomics UK consortium,; Thomson, E. C. title: Inconsistent directions of change in case severity across successive SARS-CoV-2 variant waves suggests an unpredictable future date: 2022-03-27 journal: nan DOI: 10.1101/2022.03.24.22272915 sha: ea134292210feee8c6f517b2a06dc61eb02c6c80 doc_id: 872175 cord_uid: pq3n18ae Objective To determine how the severity of successively dominant SARS-CoV-2 variants has changed over the course of the COVID-19 pandemic. Design Prospective cohort analysis. Setting Community- and hospital- sequenced COVID-19 cases in the NHS Greater Glasgow and Clyde (NHS GG&C) Health Board (1.2 million people). Participants All sequenced non-nosocomial adult COVID-19 cases in NHS GG&C identified to be infected with the relevant SARS-CoV-2 lineage during the following analysis periods. B.1.177/Alpha analysis: 1st November 2020 - 30th January 2021 (n = 1640). Alpha/Delta analysis: 1st April - 30th June 2021 (n = 5552). AY.4.2 Delta/non-AY.4.2 Delta analysis: 1st July - 31st October 2021 (n = 9613). Non-AY.4.2 Delta/Omicron analysis: 1st - 31st December 2021 (n = 3858). Main outcome measures Admission to hospital, admission to ICU, or death within 28 days of first positive COVID-19 test Results In the B.1.177/Alpha analysis, 300 of 807 (37.2%) B.1.177 cases were recorded as hospitalised or having a more severe outcome, compared to 232 of 833 (27.9%) Alpha cases. After adjusting for the following covariates: age, sex, time of positive test, comorbidities and partial postcode, the cumulative odds ratio was 1.51 (95% central credible interval 1.08-2.11) for Alpha versus B.1.177. In the Alpha/Delta analysis, 113 of 2104 (5.4%) Alpha cases were recorded as hospitalised or having a more severe outcome, compared to 230 of 3448 (6.7%) Delta cases. After adjusting for the above covariates plus number of vaccine doses and reinfection, the cumulative odds ratio was 2.09 (95% central credible interval 1.42-3.08) for Delta versus Alpha. In the non-AY.4.2 Delta/AY.4.2 Delta analysis, 845 of 8644 (9.8%) non-AY.4.2 Delta cases were recorded as hospitalised or having a more severe outcome, compared to 101 of 969 (10.4%) AY.4.2 Delta cases. After adjusting for the previously stated covariates, the cumulative odds ratio was 0.99 (95% central credible interval 0.76-1.27) for AY.4.2 Delta versus non-AY.4.2 Delta. In the non-AY.4.2 Delta/Omicron analysis, 30 of 1164 (2.6%) non-AY.4.2 Delta cases were recorded as hospitalised or having a more severe outcome, compared to 26 of 2694 (1.0%) Omicron cases. After adjusting for the previously listed covariates, the median cumulative odds ratio was 0.49 (95% central credible interval 0.22-1.06) for Omicron versus non-AY.4.2 Delta. Conclusions The direction of change in disease severity between successively emerging SARS-CoV-2 variants of concern was inconsistent. This heterogeneity in virulence between variants, coupled with independent evolutionary emergence, demonstrates that severity associated with future SARS-CoV-2 variants is inherently unpredictable. 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 2020 [1] 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 very rapidly spread 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 estimates [4] implying that it would become dominant in the UK in early 2022. This spread was arrested by the emergence of the more transmissible and immuneevading Omicron variant, which has supplanted nearly all non-Omicron diversity [5] . 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 a vitally important in the decision-making process for 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] . To generate empirical data to test this hypothesis, we prospectively linked detailed clinical metadata and viral genomic data from NHS GG&C 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, and non-AY.4.2 Delta versus Omicron. We test the robustness of these estimates to epidemic phase bias, a bias caused by patient outcomes being correlated with the time from infection to positive test, resulting in an estimated odds ratio adjusting for time of positive to test being a biased estimator of the odds ratio adjusting for time of infection [8] . 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. 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. Sequences were generated using the ARTIC Network protocol, originally developed for Oxford Nanopore-based sequencing [9] . and derived versions adapted to Illumina and ARTIC-unrelated amplicon-based protocols. The COG-UK pipeline was used for alignment and Pango lineage assignment [10] . For the B.1.177/Alpha analysis, we included all sequenced samples with full data available on all adjustment variables from within the NHS GG&C 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 S4 . Cases with hospital acquired COVID-19, defined as a first positive PCR test occurring more than 48 hours following admission to hospital, and all cases younger than 18 were excluded. Cohorts and de-identified linked data were prepared by the West of Scotland Safe Haven at NHS GG&C. Data used in the analysis included admission to hospital, date of hospital admission, date of hospital discharge, date of death (inside or outside hospital) and 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 HDU/ICU; 4. Death), age at diagnosis, 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. 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. The four-level patient outcome data were analysed using cumulative generalised additive mixed models (GAMMs) with logit links [11] fit using Bayesian inference. These GAMMs included lineage, reinfection, patient sex and number of vaccine doses as categorical fixed effects and number of ISARIC4C identified comorbidities 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 parameter 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 [12] . All presented models had no divergent transitions and effective sample sizes of over 200 for all parameters. Sensitivity to epidemic phase bias was assessed using the method of Seaman et al. [8] . We added four days to the population who experienced more extreme outcomes (i.e. hospitalisation, . CC-BY 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted March 27, 2022. ; 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. We took the Scottish sequences and looked at the count of non-reference nonsynonymous mutations found each day and used the nlstools package in R [13] to model the daily logistic growth rate of each lineage. We took defining mutations for each lineage to model the growth rate of the lineage. The defining mutations for the variants were chosen as those which were not present in the previously dominant lineage. The defining mutations were N501Y for Alpha, L452R for Delta, A222V for AY.4.2, and N501Y for Omicron. For AY.4.2 we calculated the growth rate of A222V relative to L452R to represent the advantage relative to basal delta rather than against the shifting-inproportion blend of Delta and Alpha variants. The period for each growth rate estimate covered the same window as the clinical analysis. The growth rate was taken from the growth parameter in the following equation of the regression. We used counts on both sides of the equation to down weight days with limited data providing noisier proportions. 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). The replacement of the lineages over time and their growth rates can be seen in Figure 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 that effect is positive: >0.99; "with" -median cumulative odds ratio: 1.51; 95% central interval: 1.08-2.11; probability that effect is positive: 0.99). 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 S5 . 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 that effect is positive: 0.93; "with" -median cumulative odds ratio: 1.24; 95% central interval: 0.88-1.79; probability that effect is positive: 0.88). 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 that effect is positive: >0.99; "with" -median cumulative odds ratio: 2.09; 95% central interval: 1.42-3.08; probability that effect is positive: >0.99). 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 S6 . 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 is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted March 27, 2022. ; interval: 0.98-2.16; probability that effect is positive: 0.96; "with" -median cumulative odds ratio: 1.40; 95% central interval: 0.97-2.02; probability that effect is positive: 0.97). 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 that effect is positive: 0.64; "with" -median cumulative odds ratio: 0.99; 95% central interval: 0.76-1.27; probability that effect is positive: 0.46). 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 S7 . 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 that effect is positive: 0.57; "with" -median cumulative odds ratio: 0.98; 95% central interval: 0.75-1.27; probability that effect is positive: 0.44). 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 that effect is positive: 0.06; "with" -median cumulative odds ratio: 0.49; 95% central interval: 0.22-1.06; probability that effect is positive: 0.04). 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 S8 . 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 that effect is positive: 0.01; "with"median cumulative odds ratio: 0.19; 95% central interval: 0.08-0.42; probability that effect is positive: <0.01). 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, and Omicron was associated with much less severe disease than non-AY.4.2 Delta. These conclusions were after accounting for comorbidities, changes in treatment 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 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 2022 [14] 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. 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. We also 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. Our study does, however, have some limitations. Firstly, we only include cases with sequenced genomes, and thus our sample predominantly includes cases with lower Ct, because these cases are more likely to have been sequenced. This is likely to be . CC-BY 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted March 27, 2022. ; particularly important when the Ct distribution of infections differs between variants. Additionally, this limits us to the set of individuals who have been tested by PCR, likely to represent hospitalised patients more than those in the community. Also, 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, 15, 16] ). 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 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 [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] . Our estimates are consistent with the majority of these studies. Our sample is smaller than was used in some of these studies, but we benefit from much higher resolution clinical data, and being able to control for comorbidities. 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 [28, 29] . Comparable results were observed in Danish, US and Canadian populations with a study in a Norwegian population being the exception [30] [31] [32] [33] . The US and Canadian studies also found that Delta was associated with increased risk of ICU admission and death [31, 32] . 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[34] . 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. 2022[33] 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, Canada and the US suggesting that Omicron infections are less severe than infections with Delta [14, . 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. 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 [38] [39] [40] . Testing patterns have also changed dramatically across the study, with periods of higher and lower rates of testing. All of these factors may impact the relative severity of variants. For this reason, our results cannot be used to compare the intrinsic case severity of variants that were not co-circulating at the same time. Our results demonstrate that successive variants of SARS-CoV-2 are associated with inconsistent differences in disease severity after other factors are accounted for in the analysis, including comorbidities, vaccination, previous infection and changes in treatment. In keeping with this finding, emerging data of fundamental changes in the life cycle of the Omicron variants BA.1 and BA.2 may provide a biological explanation for the substantial drop in severity associated with this variant. 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 direction of the evolution of SARS-CoV-2 virulence has not been consistent over time, it is historical trends in severity cannot be used to predict the severity of future variants. However, once a variant has emerged, the likelihood of immune evasion and the method of cell entry may be . CC-BY 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted March 27, 2022. ; 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 in disease severity in a variant with similar transmissibility to Omicron could be devastating to health systems and communities. 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 in a timely manner. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted March 27, 2022. ; is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted March 27, 2022. ; https://doi.org/10.1101/2022.03.24.22272915 doi: medRxiv preprint World Health Organization Tracking SARS-CoV-2 variants Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations Context-specific emergence and growth of the SARS-CoV-2 Delta variant SARS-CoV-2 variants of concern and variants under investigation in England: Technical briefing 29 Tracking of variants The three Ts of virulence evolution during zoonotic emergence The hypertransmissible SARS-CoV-2 Omicron variant exhibits significant antigenic change, vaccine escape and a switch in cell entry mechanism Adjusting for time of infection or positive test when estimating the risk of a post-infection outcome in an epidemic nCoV-2019 sequencing protocol v3 (LoCost) V.3. protocols.io CLIMB-COVID: continuous integration supporting decentralised sequencing for SARS-CoV-2 genomic surveillance Ordinal Regression Models in Psychology: A Tutorial An R Package for Bayesian Multilevel Models Using Stan A Toolbox for Nonlinear Regression in R: The Package nlstools Clinical severity of, and effectiveness of mRNA vaccines against, covid-19 from omicron, delta, and alpha SARS-CoV-2 variants in the United States: prospective observational study Omicron extensively but incompletely escapes Pfizer BNT162b2 neutralization SARS-CoV-2 variants of concern and variants under investigation in England: Technical briefing 34 Increased mortality in community-tested cases of SARS-CoV-2 lineage B.1.1.7. Nature Risk of mortality in patients infected with SARS-CoV-2 variant of concern 202012/1: matched cohort study Case fatality risk of the SARS-CoV-2 variant of concern B.1.1.7 in England Risk of hospital admission for patients with SARS-CoV-2 variant B.1.1.7: cohort analysis Assessment of Mortality and Hospital Admissions Associated with Confirmed Infection with SARS-CoV-2 Variant of Concern VOC-202012/01 (B.1.1.7) a Matched Cohort and Time-to-Event Analysis Severity of Severe Acute Respiratory System Coronavirus 2 (SARS-CoV-2) Alpha Variant (B.1.1.7) in England Descriptive comparison of admission characteristics between pandemic waves and multivariable analysis of the association of the Alpha variant (B.1.1.7 lineage) of SARS-CoV-2 with disease severity in inner London Genomic characteristics and clinical effect of the emergent SARS-CoV-2 B.1.1.7 lineage in London, UK: a whole-genome sequencing and hospital-based cohort study Clinical and Virological Features of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variants of Concern: A Retrospective Cohort Study Comparing B.1.1.7 (Alpha), B.1.351 (Beta), and B.1.617.2 (Delta) SARS-CoV-2 lineage B.1.1.7 is associated with greater disease severity among hospitalised women but not men: multicentre cohort study The SARS-CoV-2 Alpha variant caused increased clinical severity of disease in Scotland: a genomics-based prospective cohort analysis. medRxiv Public Health Scotland, EAVE II Collaborators. SARS-CoV-2 Delta VOC in Scotland: demographics, risk of hospital admission, and vaccine effectiveness Hospital admission and emergency care attendance risk for SARS-CoV-2 delta (B.1.617.2) compared with alpha (B.1.1.7) variants of concern: a cohort study Hospitalisation associated with SARS-CoV-2 delta variant in Denmark 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. 18-29 739 1242 30-39 574 887 40-49 382 633 50-59 247 426 60-69 120 167 70-79 28 61 80-89 11 28 90+ <5 <5