Hostname: page-component-788cddb947-jbjwg Total loading time: 0 Render date: 2024-10-08T16:46:39.666Z Has data issue: false hasContentIssue false

Estimating the share of SARS-CoV-2-immunologically naïve individuals in Germany up to June 2022

Published online by Cambridge University Press:  15 February 2023

Benjamin F. Maier*
Affiliation:
Robert Koch Institute, Berlin, Germany DTU Compute, Technical University of Denmark, Kongens Lyngby, Denmark Copenhagen Center for Social Data Science, University of Copenhagen, Copenhagen, Denmark
Annika H. Rose
Affiliation:
Robert Koch Institute, Berlin, Germany Institute for Theoretical Biology and Integrated Research Institute for the Life-Sciences, Humboldt University of Berlin, Berlin, Germany
Angelique Burdinski
Affiliation:
Robert Koch Institute, Berlin, Germany Institute for Theoretical Biology and Integrated Research Institute for the Life-Sciences, Humboldt University of Berlin, Berlin, Germany
Pascal Klamser
Affiliation:
Robert Koch Institute, Berlin, Germany Institute for Theoretical Biology and Integrated Research Institute for the Life-Sciences, Humboldt University of Berlin, Berlin, Germany
Hannelore Neuhauser
Affiliation:
Robert Koch Institute, Berlin, Germany
Ole Wichmann
Affiliation:
Robert Koch Institute, Berlin, Germany
Lars Schaade
Affiliation:
Robert Koch Institute, Berlin, Germany
Lothar H. Wieler
Affiliation:
Robert Koch Institute, Berlin, Germany
Dirk Brockmann
Affiliation:
Robert Koch Institute, Berlin, Germany Institute for Theoretical Biology and Integrated Research Institute for the Life-Sciences, Humboldt University of Berlin, Berlin, Germany
*
Author for correspondence: Benjamin F. Maier, E-mail: bfmaier@physik.hu-berlin.de
Rights & Permissions [Opens in a new window]

Abstract

After the winter of 2021/2022, the coronavirus disease 2019 (COVID-19) pandemic had reached a phase where a considerable number of people in Germany have been either infected with a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant, vaccinated or both, the full extent of which was difficult to estimate, however, because infection counts suffer from under-reporting, and the overlap between the vaccinated and recovered subpopulations is unknown. Yet, reliable estimates regarding population-wide susceptibility were of considerable interest: Since both previous infection and vaccination reduce the risk of severe disease, a low share of immunologically naïve individuals lowers the probability of further severe outbreaks, given that emerging variants do not escape the acquired susceptibility reduction. Here, we estimate the share of immunologically naïve individuals by age group for each of the sixteen German federal states by integrating an infectious-disease model based on weekly incidences of SARS-CoV-2 infections in the national surveillance system and vaccine uptake, as well as assumptions regarding under-ascertainment. We estimate a median share of 5.6% of individuals in the German population have neither been in contact with vaccine nor any variant up to 31 May 2022 (quartile range [2.5%–8.5%]). For the adult population at higher risk of severe disease, this figure is reduced to 3.8% [1.6%–5.9%] for ages 18–59 and 2.1% [1.0%–3.4%] for ages 60 and above. However, estimates vary between German states mostly due to heterogeneous vaccine uptake. Excluding Omicron infections from the analysis, 16.3% [14.1%–17.9%] of the population in Germany, across all ages, are estimated to be immunologically naïve, highlighting the large impact the first two Omicron waves had until the beginning of summer in 2022. The method developed here might be useful for similar estimations in other countries or future outbreaks of other infectious diseases.

Type
Original Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press

Introduction

The coronavirus disease 2019 (COVID-19) pandemic caused by the rapid global dissemination of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its respective variants has led to a large number of infections worldwide [1]. In Germany, around 21.4 million infections had been reported until the end of May 2022. Moreover, a large part of the population had received a primary vaccination series with one of the available COVID-19 vaccines (mRNA-vaccine by BioNTech or Moderna or a vector-based vaccine by AstraZeneca or Janssen) at the time [2]. The national COVID-19 vaccination campaign began at the end of 2020 by targeting older adults, residents of nursing homes, and healthcare workers, then shifting focus to younger adults [Reference Vygen-Bonnet3]. In August 2021, a recommendation to vaccinate adolescents aged 12–17 was issued and since December 2021, children aged 5–11 years were recommended to receive a vaccination if underlying medical conditions put them at increased risk for severe disease [Reference Vygen-Bonnet4, Reference Vygen-Bonnet5]. In Germany, recovered individuals were advised not to receive a COVID-19 vaccination until 6 months [6] or 3 months [7] have passed after infection, respectively. At the time of analysis, booster vaccinations had been recommended for all persons aged 11 years and older [Reference Koch8, Reference Vygen-Bonnet9]. A central factor that would determine how the pandemic progressed in Germany was the number of people still immunologically naïve to infection, i.e. that have neither been in contact with the virus or any of its variants nor a vaccine against them. In Germany, several serological studies have been conducted [Reference Neuhauser10, 11] but up until spring 2022 no data that extended into the time of the Omicron waves was available, particularly with respect to children [Reference Neuhauser12]. Therefore, at the time we chose a mathematical modelling approach to estimate the number of immunologically naïve individuals in order to facilitate informed decisions with regard to the future trajectory of the pandemic in Germany.

To estimate the number of people that have been in contact with either virus or vaccine, one might simply summate the number of vaccinations and the number of reported infections. However, doing so ignores the fact that (a) a considerable number of vaccinated people have suffered from additional breakthrough infections (taking into account both asymptomatic and symptomatic infections herein) [13], (b) a substantial number of previously infected people have chosen to be vaccinated in accordance with national recommendations [Reference Betsch14Reference Betsch16], (c) some individuals have suffered from multiple infections [Reference Sotoodeh Ghorbani17] and (d) the exact extent of the total number of infections as compared to the reported number of infections is unknown because (i) asymptomatic infections are less likely to be identified and reported in the national surveillance system and (ii) under-ascertainment varies regionally [Reference Santos-Hövener18, Reference Neuhauser19]. In order to estimate the overlap between the vaccinated and recovered subpopulations, one may assume that the probability of any recovered individual to be vaccinated is proportional to the probability of any individual to be vaccinated. However, this largely ignores (i) the heterogeneous dynamics of the spreading disease as well as vaccination campaigns and (ii) that vaccinated individuals are less likely to suffer from an infection than unvaccinated individuals [Reference Feikin20]. Here, we introduce modelling approaches that are devised to meet the aforementioned conditions and use them to estimate the distribution of immunologically naïve, (in the infectious-disease modelling context called ‘fully susceptible’ hereafter), recovered, and vaccinated individuals in Germany, taking into account regional and age differences. We find that although the percentage of the adult population in Germany that remained fully susceptible was expected to be in the single digits after 31 May 2022, the share of unaffected children might have been considerably larger. Due to heterogeneities in vaccine uptake across German states, these values may differ by region. Our analysis cannot answer questions regarding the quality of achieved immunity against infection or disease, because we consider neither waning of immunity nor the emergence of variants with immune evasive properties, which is difficult to predict [Reference Krueger21].

Methods

We partition the population into n G = 16 regions corresponding to the German states and n A = 5 age groups corresponding to ages ‘00–04’ (infants), ‘05–11’ (children), ‘12–17’ (adolescents), ‘18–59’ (adults), ‘60 + ’ (elderly), chosen in accordance with the population structure of publicly available vaccination data [2], i.e. into 80 subpopulations. To obtain nation-wide counts of individuals in age groups, we sum the respective results over all regions, to obtain counts of individuals for all ages, we sum over all age groups. To obtain an age-independent, nation-wide result, we sum over all ages and all regions.

As we are, first and foremost, interested in estimating the proportion of individuals S  ≡ S(t = t max) that can be considered to be fully susceptible towards infection with any SARS-CoV-2 variant per region and age group, we report a simplified model here that captures the main ideas and gives the same results for S(t) as the full model which is reported in the SM, Sec. ‘Introduction’.

We consider the population of size N (an age group in a region) to be composed of susceptible (S), infected/recovered (I), infected/recovered but eligible for reinfection or vaccination (Y), vaccinated (V) and boostered (B) individuals, assuming that the population count is constant over two years such that N = S + I + Y + V + B = const.

The central problem of estimating S is to determine the overlap between recovered and vaccinated subpopulations. Given that the cumulative number of unvaccinated infected R and the number of cumulative vaccinated individuals V is known, one may naïvely assume that the probability that an infected person that was initially unvaccinated is vaccinated later on is proportional to the probability that any person in the population is vaccinated, which is given as p = V /N. This is in agreement with results of a representative survey study that suggested that recovered individuals had the same intention to vaccinate as individuals that had not yet suffered from an infection [Reference Betsch14]. Based on this assumption, the cohort size of unvaccinated and not yet infected individuals would evaluate to S  = N − (1 − V /N)R  − V . However, doing so largely ignores the time course of infections and vaccinations, with incidence and daily vaccinations peaking at different time points, where a large number of infections occurred after the peak in weekly administered vaccines. Hence, one may assume instead that when a person becomes vaccinated at time t, the probability that this person was already infected is proportional to the number of infected/recovered individuals at time t that are eligible for vaccination as p = Y/(S + Y). With incidence rates of $a_\phi \phi ( t)$ (new unvaccinated cases per day) and vaccination rates of $a_\beta \beta _S( t)$ (new vaccinations per day) obtained from data, we assume that the count of individuals in the respective states evolves dynamically as

(1)$$\partial _tS = {-}a_\phi \phi ( t) \displaystyle{S \over {( 1-r) Y + S}}-a_\beta \beta _S( t) \displaystyle{S \over {Y + S}}$$
(2)$$\partial _tI = a_\phi \phi ( t) -\displaystyle{I \over \tau }$$
(3)$$\partial _tY = \displaystyle{I \over \tau }-a_\beta \beta _S( t) \displaystyle{Y \over {Y + S}}-a_\phi \phi ( t) \displaystyle{{( 1-r) Y} \over {( 1-r) Y + S}}$$
(4)$$\partial _tV = a_\beta \beta _S( t) -a_\beta \beta _V( t) $$
(5)$$\partial _tB = a_\beta \beta _V( t).$$

The last two equations are shown here for completeness, but note that the number of vaccinated and boostered individuals can simply be obtained from data, without integrating the dynamic equations, as their integrals can be evaluated analytically and are equal to the cumulative number of respective vaccinations. Above, $a_\phi$ and $a_\beta$ are under-ascertainment ratios that account for infections and vaccinations that have not been reported. The time scale τ is equal to the average time after which an infected/recovered individual becomes eligible for reinfection or vaccination and 1 − r is the relative probability that an unvaccinated recovered person is reinfected as compared to a fully susceptible individual (Fig. 1).

Fig. 1. Simplified model schema. On each day, $a_\beta \beta _S( t) {\rm \Delta }t$ unvaccinated people become vaccinated, with under-ascertainment ratio $a_\beta$ and ${\rm \Delta }t = 1{\kern 1pt} {\rm day}$. The probability that a newly vaccinated person has been infected before is proportional to the respective size of the subpopulation of recovered people that are eligible for vaccination Y. Furthermore, on each day, $a_\phi \phi ( t) {\rm \Delta }t$ unvaccinated people become infected, with under-ascertainment ratio $a_\phi$. The probability that a newly infected person has been infected before is proportional to the respective size of the subpopulation of recovered people that are eligible for reinfection (1 − r)Y, where 1 − r is the relative reinfection probability or ‘recovered immunity’. Recovered individuals are expected to reach eligibility for reinfection/vaccination after an average duration of τ. (Note that in the full model, breakthrough and reinfections of vaccinated individuals are possible (see SM, Sec. ‘Introduction’)).

For our analysis, we draw 1000 pairs of $a_\phi$ and $a_\beta$ from shifted Gamma distributions with means $\left\langle {a_\phi } \right\rangle = 2$, $\left\langle {a_\beta } \right\rangle = 1.03$ and standard deviations ${\rm Std}[ a_\phi ] = 1$, ${\rm Std}[ a_\beta ] = 0.02$ that are bounded below by $\min ( a_ \bullet ) = 1$. Note that this distribution yields a median under-ascertainment ratio of $Q_2[ a_\phi ] = 1.7$, which is in line with results informed by seroprevalence data for Germany in 2020, published as a preprint [Reference Neuhauser19]. Furthermore, with a 97.5th percentile of 4.7, the distribution is broad enough to account for occasional high under-ascertainment ratios that have been observed locally [Reference Neuhauser10, Reference Santos-Hövener18, Reference Neuhauser19]. For infants, ascertainment is expected to be lower than for other age groups [Reference Verity22], which is why we double under-ascertainment ratios for this age group. We did not assume a higher under-ascertainment ratio for children older than 4 years, because regular screening via rapid antigen tests was mandatory in schools across the country for a considerable amount of time [Reference Berger, Fritz and Kauermann23]. We choose an eligibility time of $\tau = 90{\kern 1pt} {\rm days}$, which is approximately of the same order as the time for antibody concentrations to decay after an infection [Reference Iyer24]. While it falls in the lower bound of officially recommended time for recovered individuals to wait before getting vaccinated, surveys indicate that people might not strictly follow the official recommendation but get vaccinated earlier. Further, people with asymptomatic courses might have no knowledge about their infection, likely leading to a bias towards shorter times between infection and vaccination in those cases. The influence of lower and higher values of τ is investigated in a sensitivity analysis (see SM, Sec. ‘Methods’). The ‘recovered immunity’ parameter r quantifies the relative efficacy against reinfection. For the Alpha variant, this efficacy was observed to be lower than the vaccine efficacy against infection by mRNA- or vector-vaccines [25], but of similar order as the vaccine efficacy against infection with Delta, taking on values of r ≈ 0.65 for both. As Omicron is considered to be a variant with partial immune escape, we set a lower default value of r = 1/2 for all variants, testing r = 0 (no protection against reinfection) and r = 1 (full immunity) in sensitivity analyses.

Note that we only consider data up to 31 May 2022. We cannot reliably extrapolate the above-listed assumptions regarding under-ascertainment into the summer of 2022 and beyond because the Omicron variant and an increasing population immunity may have altered subjective perception of disease severity and resulting test usage [Reference Kolb26].

The daily vaccination rates $\beta _ \bullet ( t)$ are obtained from data [2] and averaged over calendar weeks to remove weekly modulations. Likewise, infection rates of unvaccinated individuals ϕ(t) are obtained from reported data in the German reporting system SurvStat [Reference Faensen and Krause27], which is available in aggregated form upon request. While the vaccination status is unknown for a substantial number of infections, we assume that for every day, the proportion of cases with unknown vaccination status that are, in fact, unvaccinated, is equal to the proportion of unvaccinated cases over the last seven days for which the vaccination status is known. This imputation method is performed for age- and region-stratified data.

For analyses disregarding infections with Omicron, we obtained the nation-wide and age-independent share of randomly sequenced samples in Germany [28] that the software framework ‘scorpio’ identified as ‘Omicron’ or ‘Probable Omicron’ on a per-calendar-week basis by date of extraction (‘Entnahmedatum’) as σ(t), assuming σ(t) = 0 for dates previous to 1 August 2021 and σ(t) = 1 for dates that exceed the last available date in the data. Then, all incidence rates were scaled as ϕ S,pre-Omicron(t) = ϕ S(t)[1 − σ(t)]. Note that vaccination rates are unaffected by this procedure.

Population sizes stratified by age and state were requested from destatis [29].

Eqs. (1)–(5) are integrated using Euler's method with ${\rm \Delta }t = 1{\kern 1pt} {\rm day}$ until the last day of available incidence/vaccination data. For dates where data is unavailable, we assume the respective rates are equal to zero.

Results

We find an estimated nationwide median share of fully susceptible individuals of 5.6% (quartile range [2.5%–8.5%]) until 31 May 2022. This result is, however, biased towards higher values due to a larger share of unaffected infants (36.3% [18.8%–49.8%]) and children (16.7% [4.2%–28.6%]). For age groups that are associated with a higher probability of severe disease [Reference Barek, Aziz and Islam30], we find a lower relative frequency of 3.8% [1.6%–5.9%] (adults), and 2.1% [1.0%–3.4%] (elderly), which is of the same order as the values for adolescents (1.7% [0.1%–5.2%]).

These values are achieved largely due to the respective Omicron waves in early 2022. Ignoring infections with the Omicron variant, the nationwide age-independent share of fully susceptibles increases to 16.3% [14.1%–17.9%], i.e. Omicron infections are expected to have caused a reduction in fully susceptible individuals on the order of 10 percentage points up to June 2022, although this number differs by age group. While the change in relative frequency of fully susceptibles in the ‘adult’ and ‘elderly’ age groups was only about a few percentage points (median decreases from 10.8% to 3.8% and from 4.8% to 2.1%, respectively), the three youngest age groups were affected much more strongly, with median values of fully susceptible individuals dropping from 83.3% to 36.3%, from 62.6% to 16.7% and from 20.3% to 1.7% with increasing age (cf. Fig. 2). If all variants are considered, the median share of fully susceptible ‘adults’ and ‘elderly’ barely differs (absolute difference of 1.7% points), likely due to a larger fraction of Omicron-recovered ‘adults’ (Fig. 2).

Fig. 2. Estimated nationwide relative frequency of fully susceptible individuals by age group, considering vaccinations and infections that took place up to and including May 2022. Boxes represent the area between quartiles Q 1, Q 3 and whiskers the 2.5th and 97.5th percentiles, respectively, the median is shown as a horizontal line. (Left) Considering infections with any variant. (Right) Considering infections with any variant other than Omicron and its sublineages.

Although the relative frequency of fully susceptibles varies between federal states, certain commonalities are still shared. In all states, the frequency of fully susceptible individuals decreases with age, with a strong dependence on age for children. For ages 12–17, the frequency reaches values on the same order as those of the age groups ‘adults’ and ‘elderly’ (Fig. 3). Apart from the fact that adult and elderly age groups achieve relative frequencies of fully susceptible individuals below 10%, there are no other common patterns that stand out across all states regarding these age groups. In general, these age groups show overlapping quartile intervals, with the exception of Hamburg and Bremen, where ‘adults’ show a comparatively lower relative frequency (Fig. 3). In fact, in Bremen virtually no one aged 18 and above is expected to not have been in contact with either virus or vaccine, according to the estimations.

Fig. 3. Estimated relative frequency of fully susceptible individuals by age group and region considering infections with any variant and vaccinations up to and including the Omicron wave (as of 31 May 2022).

In general, the above observations hold for the pre-Omicron analysis as well, except for the fact that, in the majority of states, the number of adults that were still unaffected decreased dramatically during the Omicron wave due to the large number of infections caused by the variant (comparing Figs 3 and 4). When excluding Omicron infections, the relative frequency of fully susceptibles differs across states on the order of about 10 percentage points, with Brandenburg and Bremen as the states with largest (21.7%) and smallest (9.6%) respective median values of fully susceptible individuals (Fig. 4). Including infections with Omicron, the median range between states is reduced to a difference of 7.9% points (median of 9.4% in Hesse and 1.4% in Saxony-Anhalt).

Fig. 4. Estimated relative frequency of fully susceptible individuals by age group and region, disregarding infections with Omicron and its sublineages, based on data available up to and including May 2022.

Our results are robust against changes in assumed eligibility time τ and recovered immunity r, varying by a few percentage points in the nationwide average for all ages. For the most at-risk age groups, i.e. adults and the elderly, these results vary even less, indicating that the influence of these parameters decreases with age (see SM, Sec. ‘Methods’ and Supplementary Fig. S3).

Regarding the detailed distribution of individuals by vaccination/infection status, we find that the largest single compartment of the model population is the group of people that has received a booster vaccination and has never been in contact with the virus (see SM Sec. ‘Methods’ and Supplementary Fig. S2), with unvaccinated recovereds comprising the second largest group. When first excluding, then including Omicron infections, both the number of non-infected vaccinateds and non-infected booster vaccinateds decreases by about 10 percentage points, demonstrating the relative efficacy of the booster vaccination against infections with the Omicron variant. The prevalence of compartments that count infected individuals decreases with the number of (breakthrough) infections per individual, which is unsurprising given that the model probability to become infected decreases exponentially with every new infection. Note that our model cannot, however, track the number of reinfections per individual between achieving the different vaccination statuses.

As under-ascertainment is expected to be larger for infants than for other age groups, we scaled the respective under-ascertainment ratio to always assume twice the value of other age groups. Because most children below 5 years of age will remain unvaccinated as per official recommendations, only infections reduce the number of fully susceptible individuals, and, therefore, the under-ascertainment ratio has a large influence (see SM, Sec. ‘Methods’ and Supplementary Fig. S4). With the degree of under-ascertainment in this age group comparatively unclear, the results must be considered relatively uncertain for this age group.

Conclusion & discussion

As the pandemic progresses, a central quantity that will determine the upcoming dynamics is the population-wide susceptibility against infection with known or future variants of SARS-CoV-2. While protection from infection, either derived from vaccination or natural infection, wanes over time and depends on the circulating virus variant, an estimation of the respective subpopulation sizes of people that suffered from one or more infections or were vaccinated/boostered gives valuable information about the size of the population that is still fully susceptible to infection, because these individuals are more prone to infection and severe disease as compared to vaccinated or recovered individuals, given that future variants do not fully escape this immunity.

Here, we found that in Germany, a nationwide single-digit percentage of individuals have not been in contact with either a variant of SARS-CoV-2 or a vaccine against them up to June 2022, yet these results vary between regions and age groups. Despite the high number of reported infections in infants and children, a considerably high percentage of these age groups might still have been fully susceptible to infection after this time. This may become problematic if a variant emerges that causes more severe disease in these age groups than previous variants. Yet, we cannot rule out the possibility that we underestimated the extent of under-ascertainment in these age groups, as the factors we used where informed by seroprevalence studies based on blood samples donated by adults (ages 18–74), while it has been reported that under-ascertainment ratios can assume values ranging from 2 up to 6 or 8 for children [Reference Brinkmann31Reference Wachter33].

In comparison, the age groups of adults and elderly showed a relatively low share of fully susceptible individuals, considering infections with all variants on the order of 5%. Only considering infections with pre-Omicron variants, however, around 9.0%–12.4% of the adult population and 3.6%–5.8% of the elderly population were expected to still be at risk of infection with variants that have a higher probability of causing severe disease than Omicron, potentially causing large outbreaks that could put high pressure on the public health system once again (with these numbers representing quartile ranges).

Across all ages and regions, we found that about 6% of the German population were expected to have not been in contact with either virus or vaccine. While at the time of analysis, no seroprevalence data that extended into the Omicron waves were available [Reference Neuhauser12], recently an interim analysis was published in a repository [Reference Lange34], based on blood samples that were taken between May and September 2022. The report concluded that 95% of the participants showed antibodies against the S-antigen, across all ages in Germany. The results are not directly comparable because these samples were taken after our study period ended and, as discussed, antibody concentration in the blood is known to wane over time, but they are in line with our estimates and support the validity of our method.

Our results are subject to a number of limitations and biases. For instance, the reported uncertainties (quartile ranges) are heavily determined by the choice of distribution of $a_\phi$. The distribution we chose has a median value of $Q_2[ a_\phi ] = 1.7$, which is slightly lower than what was observed in 2020 [Reference Neuhauser19]. Moreover, the lower distribution bound of $\min ( a_\phi ) = 1$ might be rather low, as such a value would mean that every infection has been reported, which is unlikely. Hence, at least the upper percentiles we report for S might be overestimations. Furthermore, we assume the same distribution of under-ascertainment ratios for all German states, which might not reflect potential heterogeneities in local ascertainment particularly well.

Regarding modelling choices for the eligibility time, a short average duration after infection to be eligible for vaccination leads to larger proportions of vaccine-eligible people and, hence, to a higher overlap between the vaccinated and recovered subpopulations, thus increasing the estimated number of fully susceptibles. While we chose a comparably low value of 90 days for this parameter, lower values cannot be ruled out. However, (i) the value we chose lies below the official recommendation and (ii) changes in this parameter are not expected to change our results drastically, as was shown in a sensitivity analysis.

Likewise, shorter durations of eligibility for reinfection and lower values of long-term immunity of recovered individuals increase the likelihood that a reported infection of an unvaccinated individual was, in fact, a reinfection event, thus leading to higher values of fully susceptible individuals over all. As above, our results are robust towards variations in these parameters.

Regarding results on a regional level, reported vaccinations and infections might be skewed regionally when a large number of people live in one state but traverse to others to seek medical help. These considerations might explain the extreme results observed for Hamburg and Bremen, which are city states enclosed by others.

The last German census took place in 2011 and population sizes per age group and region have been imputed for the year 2020 based on this data, thus potentially being subject to over- or under-counting. Uncertainties in population size may introduce systematic errors on the order of a few percentage points in relative frequencies. When such a relative frequency reaches low values, these absolute errors on the order of a few percentage points can lead to high relative errors in the results.

Considering incidence rates, we imputed the total number of unvaccinated cases per day from cases with undetermined vaccination status by assigning them the ‘unvaccinated’ status with probability proportional to the share of unvaccinated cases in the set of cases with determined status. This procedure can introduce systematic errors when the ascertainment of vaccination status is biased towards any of the vaccination states, which may occur, for instance, when the probability of status ascertainment increases with severity of disease. In this case, people with breakthrough infections may be less likely to have their vaccination status reported in the reporting system, which would mean that we overestimated the number of unvaccinated cases per day, introducing a bias towards lower values of the share of fully susceptible individuals.

For analyses regarding infections with variants prior to Omicron, we relied on the nationwide share of Omicron sequences, multiplying all incidence rates (regardless of region, age or vaccine status) with this function. Since vaccines assume different efficacies against infection with different variants and will likely vary across ages and regions, this assumption is expected to introduce strong bias on a fine-grained population level, which may be expected to decrease when values are aggregated over regions or ages.

Our results cannot be used to predict the future course of the pandemic directly. In fact, since SARS-CoV-2 lacks phenotypical stability and neither infection nor vaccination elicits full long-term protective immunity, especially with respect to the prevention of infection and transmission, classical herd immunity is unlikely to be reached for COVID-19 [Reference Morens, Folkers and Fauci35]. In several studies hybrid immunity resulting from infection-acquired immunity boosted with vaccination conferred the strongest or longer-lasting protection, respectively [Reference Hall36, Reference Altarawneh37]. Similarly, Omicron breakthrough infections in previously vaccinated individuals have been shown to drive cross-variant neutralisation and memory B cell formation [Reference Quandt38], suggesting that a combination of both, natural infection and vaccination, will have more impact on the future COVID-19 epidemiology than one of the events alone.

To sum up, our study shows that, presumably, only a small part of the German population has not yet been in contact with either a variant of SARS-CoV-2 or a respective vaccine against the disease they cause, up to and including May 2022. We show important proportions of fully susceptible elderly, who on average, by their age and age-associated morbidities, have a disproportionately elevated risk of severe disease. These shares differ by region and could motivate regionally targeted protection measures at the time of writing or in case of future outbreaks.

While the immunisation campaign was successful in spring and summer 2021, in particular reaching a large proportion of vulnerable people, it thereafter had difficulties to completely close immunity gaps with vaccinations, albeit enhancing the protection of a large proportion of already vaccinated people with a large booster vaccination campaign by the end of 2021. Our results show that the Omicron wave had a high impact on naturally closing the aforementioned gaps. As mentioned above, however, having been in contact with a variant of SARS-CoV-2 is not a robust equivalent of immunity and may range from mild infection followed by rapid waning of antibodies and a highly uncertain degree of immunity, to a fully vaccinated status including a booster and a breakthrough infection, which confers a more long-lasting and robust degree of protection against severe disease. At the lower end of this spectrum of presumed immunity, our analyses show that up to June 2022, one in six persons was never vaccinated but infected once or more, in the majority of cases with Omicron. This group faces higher uncertainties for later infection waves since protection against severe disease is expected to be more short-lived and too narrowly targeted to this variant.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0950268823000195.

Acknowledgements

The authors would like to thank Maria Waize and Matthias an der Heiden for helpful discussions. While working on this study B.F.M received funding through Grant CF20-0044, HOPE: How Democracies Cope with Covid-19, from the Carlsberg Foundation and was supported as an Add-on Fellow for Interdisciplinary Life Science by the Joachim Herz Stiftung.

Author contributions

All authors read and approved the manuscript. The corresponding authors had full access to all the data in the study and had final responsibility for the decision to submit for publication.

Conflict of interest

The authors declare no conflict of interests.

Data availability statement

The code and results of the simulations can be obtained from Zenodo [Reference Maier39].

References

World Health Organisation (2021) WHO Coronavirus (COVID-19) Dashboard. Available at https://covid19.who.int (Accessed 19 November 2022).Google Scholar
Robert Koch-Institut, Fachgebiet 33 (2021) COVID-19-Impfungen in Deutschland. Zenodo. Type: dataset. Available at: https://doi.org/10.5281/zenodo.5126652.Google Scholar
Vygen-Bonnet, S et al. (2020) Beschluss und Wissenschaftliche Begründung der Ständigen Impfkommission (STIKO) für die COVID-19-Impfempfehlung. Epidemiologisches Bulletin 2, 363.Google Scholar
Vygen-Bonnet, S et al. (2021) Beschluss der STIKO zur 9. Aktualisierung der COVID-19–Impfempfehlung und die dazugehörige wissenschaftliche Begründung. Epidemiologisches Bulletin 33, 346.Google Scholar
Vygen-Bonnet, S et al. (2021) STIKO-Empfehlung zur COVID-19-Impfung von Kindern im Alter von 5–11 Jahren mit dem mRNA-Impfstoff Comirnaty und die dazugehörige wissenschaftliche Begründung. Epidemiologisches Bulletin 1, 1647.Google Scholar
Ständige Impfkommission (STIKO) (2021) Beschluss der STIKO zur 10. Aktualisierung der COVID-19–Impfempfehlung. Epidemiologisches Bulletin 38, 39.Google Scholar
Ständige Impfkommission (STIKO) (2022) Beschluss der STIKO zur 18. Aktualisierung der COVID-19–Impfempfehlung. Epidemiologisches Bulletin 7, 318.Google Scholar
Koch, J et al. (2021) Wissenschaftliche Begründung der STIKO zur Aktualisierung der Empfehlung der COVID-19-Auffrischimpfung mit einem mRNA-Impfstoff. Epidemiologisches Bulletin 48, 1541.Google Scholar
Vygen-Bonnet, S et al. (2022) STIKO-Empfehlung zur COVID-19-Auffrischimpfung mit einem mRNA-Impfstoff für 12–17-jährige Kinder und Jugendliche mit dem mRNA-Impfstoff Comirnaty und die dazugehörige wissenschaftliche Begründung. Epidemiologisches Bulletin 3, 2131.Google Scholar
Neuhauser, H et al. (2021) Seroepidemiological studies on SARS-CoV-2 in samples from the general population and blood donors in Germany – findings up to August 2021. Epidemiologisches Bulletin 37, 312. Available at https://www.rki.de/EN/Content/infections/epidemiology/outbreaks/COVID-19/EB-37-2021-Art_02_en.pdf.Google Scholar
Robert Koch Institute (2022) Seroepidemiological Studies in Germany. Available at www.rki.de/covid-19-serostudies-germany.Google Scholar
Neuhauser, H et al. (2022) Observatorium serologischer Studien zu SARS-CoV-2 in Deutschland. Zenodo. Type: dataset. Available at: https://doi.org/10.5281/zenodo.7043025.Google Scholar
Robert Koch Institute (2022) Wöchentlicher Lagebericht des RKI zur Coronavirus-Krankheit-2019 (COVID-19) – 31.03.2022. Available at https://www.rki.de/DE/Content/InfAZ/N/Neuartiges_Coronavirus/Situationsberichte/Wochenbericht/Wochenbericht_2022-03-31.pdf.Google Scholar
Betsch, C et al. (2021) COSMO – Zusammenfassung und Empfehlungen Wellen 52 und 53. Available at https://projekte.uni-erfurt.de/cosmo2020/web/summary/52-53/.Google Scholar
Betsch, C et al. (2021) COVID-19 Snapshot Monitoring (COSMO Germany) – Wave 52. PsychArchives. Type: dataset. Available at https://www.psycharchives.org/handle/20.500.12034/4564.Google Scholar
Betsch, C et al. (2021) COVID-19 Snapshot Monitoring (COSMO Germany) – Wave 53. PsychArchives. Type: dataset. Available at https://www.psycharchives.org/handle/20.500.12034/4570.Google Scholar
Sotoodeh Ghorbani, S et al. (2022) Epidemiologic characteristics of cases with reinfection, recurrence, and hospital readmission due to COVID-19: a systematic review and meta-analysis. Journal of Medical Virology 94, 4453. Available at https://onlinelibrary.wiley.com/doi/10.1002/jmv.27281.CrossRefGoogle ScholarPubMed
Santos-Hövener, C et al. (2020) Serology- and PCR-based cumulative incidence of SARS-CoV-2 infection in adults in a successfully contained early hotspot (CoMoLo study), Germany, May to June 2020. Euro Surveillance: European Communicable Disease Bulletin 25, 2001752. http://doi.org/10.2807/1560-7917.ES.2020.25.47.2001752.Google Scholar
Neuhauser, H et al. (2022) Nationally representative results on SARS-CoV-2 seroprevalence and testing in Germany at the end of 2020. Scientific Reports 12, 19492.CrossRefGoogle ScholarPubMed
Feikin, DR et al. (2022) Duration of effectiveness of vaccines against SARS-CoV-2 infection and COVID-19 disease: results of a systematic review and meta-regression. The Lancet 399, 924944. Available at https://linkinghub.elsevier.com/retrieve/pii/S0140673622001520.CrossRefGoogle ScholarPubMed
Krueger, T et al. (2022) Risk assessment of COVID-19 epidemic resurgence in relation to SARS-CoV-2 variants and vaccination passes. Communications Medicine 2, 23.CrossRefGoogle ScholarPubMed
Verity, R et al. (2020) Estimates of the severity of coronavirus disease 2019: a model-based analysis. The Lancet Infectious Diseases 20, 669677.CrossRefGoogle ScholarPubMed
Berger, U, Fritz, C and Kauermann, G (2021) Eine statistische Analyse des Effekts von verpflichtenden Tests an Schulen mit Präsenzunterricht im Vergleich zum Distanzunterricht. Department of Statistics, University of Munich, Technical Reports 238, 113. https://doi.org/10.5282/ubm/epub.76005.Google Scholar
Iyer, AS et al. (2020) Persistence and decay of human antibody responses to the receptor binding domain of SARS-CoV-2 spike protein in COVID-19 patients. Science Immunology 5, eabe0367. Available at https://www.science.org/doi/abs/10.1126/sciimmunol.abe0367.CrossRefGoogle Scholar
Office for National Statistics (UK) (2021) Coronavirus (COVID-19) Infection Survey Technical Article: Impact of vaccination of testing positive in the UK: October 2021. Available at https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/articles/coronaviruscovid19infectionsurveytechnicalarticleimpactofvaccinationontestingpositiveintheuk/october2021.Google Scholar
Kolb, JJ et al. (2023) Prevalence of positive COVID-19 test results collected by digital self-report in the US and Germany. JAMA Network Open 6, e2253800-0.CrossRefGoogle ScholarPubMed
Faensen, D and Krause, G (2004) SurvStat@RKI – a web-based solution to query surveillance data in Germany. Weekly Releases (1997–2007) Euro Surveill 8, 2477. Available at: https://www.eurosurveillance.org/content/10.2807/esw.08.22.02477-en.Google Scholar
Robert Koch-Institut (2022) SARS-CoV-2 Sequenzdaten aus Deutschland. Zenodo. Type: dataset. Available at: https://doi.org/10.5281/zenodo.6381314.Google Scholar
Statistisches Bundesamt (Destatis) (2020) GENESIS-Tabelle: 12411-0012, Bevölkerung: Bundesländer, Stichtag, Altersjahre. Statistisches Bundesamt (Destatis). Type: dataset. Available at: https://www-genesis.destatis.de/genesis/online?sequenz=statistikTabellen&selectionname=12411#abreadcrumb.Google Scholar
Barek, MA, Aziz, MA and Islam, MS (2020) Impact of age, sex, comorbidities and clinical symptoms on the severity of COVID-19 cases: a meta-analysis with 55 studies and 10014 cases. Heliyon 6, e05684. Available at https://www.sciencedirect.com/science/article/pii/S2405844020325275.CrossRefGoogle Scholar
Brinkmann, F et al. (2021) Longitudinal rise in seroprevalence of SARS-CoV-2 infections in children in Western Germany – a blind spot in epidemiology? Infectious Disease Reports 13, 957964.CrossRefGoogle Scholar
Hippich, M et al. (2021) A public health antibody screening indicates a 6-fold higher SARS-CoV-2 exposure rate than reported cases in children. Med 2, 149163.e4. Available at https://linkinghub.elsevier.com/retrieve/pii/S2666634020300209.CrossRefGoogle ScholarPubMed
Wachter, F et al. (2022) Continuous monitoring of SARS-CoV-2 seroprevalence in children using residual blood samples from routine clinical chemistry. Clinical Chemistry and Laboratory Medicine (CCLM) 60, 941951. https://doi.org/10.1515/cclm-2022-0037.CrossRefGoogle ScholarPubMed
Lange, B et al. (2022) 2. Interimsanalyse des IMMUNEBRIDGE-Projektes zur Kommunikation von vorläufigen Ergebnissen an das Modellierungsnetz für schwere Infektionskrankheiten. Zenodo. Type: dataset. Available at: https://doi.org/10.5281/zenodo.7177592.Google Scholar
Morens, DM, Folkers, GK and Fauci, AS (2022) The concept of classical herd immunity may not apply to COVID-19. The Journal of Infectious Diseases 226, 195198. https://doi.org/10.1093/infdis/jiac109.CrossRefGoogle Scholar
Hall, V et al. (2022) Protection against SARS-CoV-2 after COVID-19 vaccination and previous infection. New England Journal of Medicine 386, 12071220. Available at http://www.nejm.org/doi/10.1056/NEJMoa2118691.CrossRefGoogle ScholarPubMed
Altarawneh, HN et al. (2022) Effects of previous infection and vaccination on symptomatic omicron infections. New England Journal of Medicine 387, 2134.CrossRefGoogle ScholarPubMed
Quandt, J et al. (2022) Omicron BA.1 breakthrough infection drives cross-variant neutralization and memory B cell formation against conserved epitopes. Science Immunology 7, eabq2427. Available at https://www.science.org/doi/abs/10.1126/sciimmunol.abq2427.CrossRefGoogle ScholarPubMed
Maier, BF (2022) Serocov: Estimating the distribution of COVID-19-susceptible, -recovered, and -vaccinated individuals in Germany up to June 2022. Zenodo. Type: dataset & code. Available at: https://doi.org/10.5281/zenodo.7308182.Google Scholar
Figure 0

Fig. 1. Simplified model schema. On each day, $a_\beta \beta _S( t) {\rm \Delta }t$ unvaccinated people become vaccinated, with under-ascertainment ratio $a_\beta$ and ${\rm \Delta }t = 1{\kern 1pt} {\rm day}$. The probability that a newly vaccinated person has been infected before is proportional to the respective size of the subpopulation of recovered people that are eligible for vaccination Y. Furthermore, on each day, $a_\phi \phi ( t) {\rm \Delta }t$ unvaccinated people become infected, with under-ascertainment ratio $a_\phi$. The probability that a newly infected person has been infected before is proportional to the respective size of the subpopulation of recovered people that are eligible for reinfection (1 − r)Y, where 1 − r is the relative reinfection probability or ‘recovered immunity’. Recovered individuals are expected to reach eligibility for reinfection/vaccination after an average duration of τ. (Note that in the full model, breakthrough and reinfections of vaccinated individuals are possible (see SM, Sec. ‘Introduction’)).

Figure 1

Fig. 2. Estimated nationwide relative frequency of fully susceptible individuals by age group, considering vaccinations and infections that took place up to and including May 2022. Boxes represent the area between quartiles Q1, Q3 and whiskers the 2.5th and 97.5th percentiles, respectively, the median is shown as a horizontal line. (Left) Considering infections with any variant. (Right) Considering infections with any variant other than Omicron and its sublineages.

Figure 2

Fig. 3. Estimated relative frequency of fully susceptible individuals by age group and region considering infections with any variant and vaccinations up to and including the Omicron wave (as of 31 May 2022).

Figure 3

Fig. 4. Estimated relative frequency of fully susceptible individuals by age group and region, disregarding infections with Omicron and its sublineages, based on data available up to and including May 2022.

Supplementary material: PDF

Maier et al. supplementary material

Maier et al. supplementary material

Download Maier et al. supplementary material(PDF)
PDF 547.2 KB