Unexposed populations and potential COVID-19 burden in European countries

Determinants of remaining burden

The potential remaining COVID-19 burden in a country depends on several factors, but the key determinants are: (i) the proportion of each age group in the population with naturally-acquired immunity, (ii) the proportion of each age group with vaccine-induced immunity, and (iii) the age structure of the population, which determines the population-level risk of severe COVID-19 outcomes (hospitalisation/death). Estimating the level of naturally-acquired immunity in the population from case incidence data is difficult due to variation in reporting levels over the course of the pandemic and across countries. With reliably reported death time series and high-quality estimates of the infection fatality rate (IFR), infection time series can in principle be inferred by deconvolving the death time series and scaling it by the inverse of the IFR (3,4). However, since vaccination reduces the risk of death given SARS-CoV-2 infection and most countries have followed an age-based vaccine rollout, the IFR for COVID-19 has varied asynchronously between age groups over time, and this needs to be accounted for in the calculation. We therefore use age-stratified death (5,6) and vaccination (7) data to infer age-stratified infection time series, accounting for the impact of vaccination on the IFR for each age group over time.

Vaccination and the infection fatality rate

COVID-19 vaccines have been shown to offer strong protection against disease, hospitalisation and death, but weaker protection against infection (8,9). Their efficacy against different outcomes also varies according to vaccine type, dose and virus variant (8). Thus, as vaccination coverage has increased (Figure 1A), the overall IFR has decreased, but in a way that depends on the speed and composition of the vaccine rollout (by age and vaccine product) and the changing proportions of different SARS-CoV-2 variants. We account for the varying impacts of vaccination on different disease outcomes by estimating the overall time-varying IFR for each age group for each country as a weighted average of the IFR for unvaccinated individuals and that for vaccinated individuals, averaging the efficacies of different vaccine types and doses against different outcomes for the different circulating variants according to their relative proportions in each country (see Supplement). This gives the overall country-level IFRs shown in Figure 1B. The pattern of decrease in the IFR is similar for most countries in Europe, except England and Romania, where respectively earlier and slower vaccine rollouts have led to earlier and slower declines, and the overall pre-vaccination IFR was higher (close to 1% vs an average of 0.8%) for countries with older populations, such as Italy and Portugal.

• Download figure
• Open in new tab

Figure 1. (A) Overall first-dose vaccine coverage over time since 1st December 2020 in 19 European countries and (B) corresponding population-weighted average infection fatality rates (IFR) over time.

Shaded bands in (B) show 95% credible intervals based on a truncated normal approximation to posterior distribution of IFR in (10). See Supplementary Figure S1 for vaccine coverage and IFR over time by age group for each country.

Historical burden and naturally-acquired immunity

Using the age- and time-dependent IFR for each country (Supplementary Figure S1) we infer age-stratified infection time series from age-stratified death time series (11,12) for each country (Figure 2A), and from these calculate the cumulative proportion of (unvaccinated and vaccinated) individuals who have been infected and therefore have some degree of immunity. We make the simplifying assumption that naturally-acquired and vaccine-induced immunity do not wane, given the relatively slow waning of immunity suggested by current evidence (13), and assume complete cross-protection between different variants, given limited observation of immune escape for the Alpha and Delta variants (14). We find considerable variation in the estimated proportion of the population with naturally-acquired immunity between countries, ranging from 3.0% (95% credible interval (CI) 2.0-5.2%) in Norway over the whole population to 70% (95% CI 61-84%) in Romania.

• Download figure
• Open in new tab

Figure 2. (A)-(B) Inferred age-stratified infection time series from age-stratified death time series (Supplementary Figure 2) for countries with (A) <1 million and (B) ≥1 million estimated cumulative infections. (C) Cumulative proportion infected or vaccinated against seroprevelance estimates from serological surveys in the SeroTracker database (15).

Solid lines and shaded bands show maximum a posteriori estimates and 95% credible intervals for inferred infections; dashed lines in (B) show where vaccination has contributed to immunity; points in (B) show central estimates from seroprevalence surveys plotted at the mid-point of the survey period, with geographic scope as shown in the legend. Estimated cumulative infection counts are shifted 18 days later in (B) to account for the delay between infection and seroconversion (16).

Vaccine-induced immunity and breakthrough infections

To calculate the proportion of the population with immunity acquired through infection or vaccination, or both, it is necessary to estimate the proportion of infections that are breakthrough infections, i.e. among vaccinees, to avoid double counting these individuals when tallying the overall number with immunity. Assuming that vaccine-induced immunity, like naturally-acquired immunity, wanes slowly (13) and that vaccine rollout occurred sufficiently early and quickly in most European countries that the number of vaccinated individuals rapidly dominated the number previously infected, the proportion of infections that are breakthrough infections can be approximated via the relationship between vaccine efficacy against infection and the proportion of infections amongst vaccinated individuals from the “screening method” for calculating vaccine efficacy (17) (see Supplement). This approach suggests that in many countries there are very few or no individuals in older age groups (40+yrs) who are both unvaccinated and unexposed (Figure 2B and Figure 3A), and that the majority of unexposed unvaccinated individuals are in younger age groups (<40yrs), with the overall proportion of unvaccinated and unexposed individuals varying from 0% (95% CI 0-1%) in Hungary to 31% (95% CI 8-34%) in Slovenia.

• Download figure
• Open in new tab

Figure 3. (A) Population pyramids showing the proportions of each age group in each country that are partially/fully vaccinated (and not previously infected), previously infected, and unvaccinated and unexposed. (B) Estimated maximum remaining hospitalisations per 100,000 population by age group, country and vaccination status (assuming no waning of immunity or emergence of immune escape variants), and (C) relationship between overall remaining hospitalisations per 100,000 population and proportion who have received at least one vaccine dose across countries.

Note upper limit of horizontal axis in (B) has been truncated to ensure differences between countries remain visible and vertical axis in (C) is on log scale. See Supplementary Figure S3 for pattern of maximum remaining deaths by age and country and its relationship with vaccine coverage, and Supplementary Figure S4 for the relationship between remaining hospitalisations and deaths and population age structure.

Ethical statement

Ethical approval was not necessary for this modelling study as the analysis uses only publicly available aggregated secondary data.

Remaining hospitalisations and deaths

Having estimated the remaining number of unvaccinated and unexposed individuals in each age group in each country, we calculate the potential remaining COVID-19 burden as the number of hospitalisations/deaths that would occur per 100,000 individuals if the entire population were to be (re-)exposed now (Figure 3). In other words, the burden if all unexposed and unvaccinated individuals were infected and a proportion of vaccinated individuals and previously infected individuals were infected dependent on the protection against infection afforded by vaccination/previous infection. This corresponds to the maximum potential short-term burden of hospitalisations and deaths as it does not account for waning of immunity, population growth or potential emergence of immune escape variants, and assumes that vaccine coverage, the infection hospitalisation rate (IHR) and the IFR stay at their current levels (Figure 1 and Supplementary Figure S1).

The estimated maximum overall remaining hospitalisations per 100,000 individuals ranges widely across countries, from 62 (95% CI 58-67) in England to 840 (95% CI 760-920) in Romania (Figure 3C). This reflects the combination of high vaccine coverage and high prior levels of exposure in England and the low vaccine coverage, particularly among older individuals, in Romania than other countries (Figure 3A). In absolute terms, the estimated maximum number of remaining hospitalisations ranges from 6,300 (95% CI 5,500-7,700) for Denmark to 280,000 (95% CI 260,000-310,000) for Germany. The comparative pattern of maximum remaining COVID-19 deaths across countries is similar (Supplementary Figure S3), given the strong influence of age on the IFR as well as the IHR (10,18), with the burden of deaths per 100,000 people ranging from 19 (95% CI 18-20) for England to 360 (95% CI 330-380) for Romania, and absolute numbers of deaths ranging from 1,200 (95% CI 1,000-1,400) for Slovenia to 115,000 (95% CI 104,000-129,000) for Germany. Aside from Romania, where vaccine coverage is low, countries with a combination of lower vaccine coverage among older age groups, relatively low prior exposure and older populations (Austria, Finland, Germany, Greece, the Netherlands, Slovenia) have the highest maximum remaining burdens (Figure 3C and Supplementary Figure S4), as they have the potential to experience much higher numbers of hospitalisations and deaths among the elderly than countries with younger populations and high coverage in older age groups.

Discussion and conclusion

Our results suggest that the potential remaining burden of COVID-19 hospitalisations and deaths across the 19 European countries considered is substantial, amounting to over 900,000 hospitalisations and 300,000 deaths, and that it varies considerably between countries, with countries that have experienced less transmission so far, have lower vaccine coverage and/or have older populations having much higher potential outstanding burdens. These estimates represent an upper bound on the short-term burden that could theoretically occur from a complete lifting of non-pharmaceutical interventions and return to pre-pandemic contact levels. We do not consider longer-term burden, as this would depend on the impact of waning immunity on severe disease, which remains to be established (13), and would also be influenced by susceptible individuals being replenished through births. Continued non-pharmaceutical interventions and efforts to achieve high vaccine coverage are therefore required in the short term in these countries to limit severe COVID-19 outcomes. Exactly how transmission dynamics evolve in different countries in the longer term will depend on local policies, and this will be explored through more detailed modelling building on this analysis.

Funding

This work was funded by the National Insitute for Heatlh Research (NIHR) (NIHR200908) and the Wellcome Trust (206250/Z/17/Z) (LACC, TWR, AJK). The work was partly supported by funding from the European Union’s Horizon 2020 research and innovation programme - project EpiPose (101003688: RCB); the FCDO/Wellcome Trust (Epidemic Preparedness Coronavirus research programme 221303/Z/20/Z: KvZ); the Wellcome Trust (210758/Z/18/Z: SA); NIHR (NIHR200929) (NGD); UK MRC (MC_PC_19065) (NGD) and UKRI Research England (NGD). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Authors’ contributions

LACC and AJK conceptualised the study and devised the methodology. LACC carried out the analysis and wrote the first draft of the manuscript. RCB compiled the information on vaccine efficacies. RCB, TWR, SA, KvZ and NGD contributed to the analysis and interpretation of the data and presentation of figures, and reviewed and edited the manuscript. All authors approved the final version of the manuscript for submission.