Inferring the differences in incubation-period and generation-interval distributions of the Delta and Omicron variants of SARS-CoV-2.

Estimating the differences in the incubation-period, serial-interval, and generation-interval distributions of SARS-CoV-2 variants is critical to understanding their transmission. However, the impact of epidemic dynamics is often neglected in estimating the timing of infection-for example, when an epidemic is growing exponentially, a cohort of infected individuals who developed symptoms at the same time are more likely to have been infected recently. Here, we reanalyze incubation-period and serial-interval data describing transmissions of the Delta and Omicron variants from the Netherlands at the end of December 2021. Previous analysis of the same dataset reported shorter mean observed incubation period (3.2 d vs. 4.4 d) and serial interval (3.5 d vs. 4.1 d) for the Omicron variant, but the number of infections caused by the Delta variant decreased during this period as the number of Omicron infections increased. When we account for growth-rate differences of two variants during the study period, we estimate similar mean incubation periods (3.8 to 4.5 d) for both variants but a shorter mean generation interval for the Omicron variant (3.0 d; 95% CI: 2.7 to 3.2 d) than for the Delta variant (3.8 d; 95% CI: 3.7 to 4.0 d). The differences in estimated generation intervals may be driven by the "network effect"-higher effective transmissibility of the Omicron variant can cause faster susceptible depletion among contact networks, which in turn prevents late transmission (therefore shortening realized generation intervals). Using up-to-date generation-interval distributions is critical to accurately estimating the reproduction advantage of the Omicron variant.

The Fading of the Mpox Outbreak Among Men Who Have Sex With Men: A Mathematical Modelling Study.

BACKGROUND: In the Netherlands, the number of mpox cases started declining before mpox vaccination was initiated. Most cases were men who have sex with men (MSM). We investigated whether the decline in mpox could be attributed to infection-induced immunity or behavioral adaptations. METHODS: We developed a transmission model and accounted for possible behavioral adaptations: fewer casual partners and shorter time until MSM with mpox refrain from sexual contacts. RESULTS: Without behavioral adaptations, the peak in modelled cases matched observations, but the decline was less steep than observed. With behavioral adaptations in the model, we found a decline of 16%-18% in numbers of casual partners in June and 13%-22% in July 2022. Model results showed a halving of the time before refraining from sex. When mpox vaccination started, 57% of MSM with very high sexual activity in the model had been infected. Model scenarios revealed that the outbreak could have waned by November 2022 even without vaccination. CONCLUSIONS: The limited duration of the mpox outbreak in the Netherlands can be ascribed primarily to infection-induced immunity among MSM with high sexual activity levels. The decline was accelerated by behavioral adaptations. Immunity among those most sexually active is essential to impede mpox resurgence.

Dynamic landscape of mpox importation risks driven by heavy-tailed sexual contact networks among men who have sex with men in 2022.

BACKGROUND: During the 2022 global mpox outbreak, the cumulative number of countries reporting their first imported case quickly rose in the early phase, but the importation rate subsequently slowed down, leaving many countries reporting no cases by the 2022 year-end. METHODS: We developed a mathematical model of international dissemination of mpox infections incorporating sexual networks and global mobility data. We used this model to characterize the mpox importation patterns observed in 2022 and to discuss the potential of further international spread. RESULTS: Our proposed model better explained the observed importation patterns than models not assuming heterogeneity in sexual contacts. Estimated importation hazards decreased in most countries, surpassing the global case count decline, suggesting a reduced per-case risk of importation. We assessed each country's potential to export mpox cases until the end of an epidemic, identifying countries capable of contributing to the future international spread. CONCLUSIONS: The accumulation of immunity among high-risk individuals over highly heterogeneous sexual networks may have contributed to the slowdown in the rate of mpox importations. Nevertheless, the existence of countries with the potential to contribute to the global spread of mpox highlights the importance of equitable resource access to prevent the global resurgence of mpox.

Time Scales of Human Mpox Transmission in The Netherlands.

Mpox has spread rapidly to many countries in nonendemic regions. After reviewing detailed exposure histories of 109 pairs of mpox cases in the Netherlands, we identified 34 pairs where transmission was likely and the infectee reported a single potential infector with a mean serial interval of 10.1 days (95% credible interval, 6.6-14.7 days). Further investigation into pairs from 1 regional public health service revealed that presymptomatic transmission may have occurred in 5 of 18 pairs. These findings emphasize that precaution remains key, regardless of the presence of recognizable symptoms of mpox.

Estimating the risk and spatial spread of measles in populations with high MMR uptake: Using school-household networks to understand the 2013 to 2014 outbreak in the Netherlands.

BACKGROUND: Measles outbreaks are still routine, even in countries where vaccination coverage exceeds the guideline of 95%. Therefore, achieving ambitions for measles eradication will require understanding of how unvaccinated children interact with others who are unvaccinated. It is well established that schools and homes are key settings for both clustering of unvaccinated children and for transmission of infection. In this study, we evaluate the potential for contacts between unvaccinated children in these contexts to facilitate measles outbreaks with a focus on the Netherlands, where large outbreaks have been observed periodically since the introduction of mumps, measles and rubella (MMR). METHODS AND FINDINGS: We created a network of all primary and secondary schools in the Netherlands based on the total number of household pairs between each school. A household pair are siblings from the same household who attend a different school. We parameterised the network with individual level administrative school and household data provided by the Dutch Ministry for Education and estimates of school level uptake of the MMR vaccine. We analysed the network to establish the relative strength of contact between schools and found that schools associated with low vaccine uptake are highly connected, aided by a differentiated school system in the Netherlands (Coleman homophily index (CHI) = 0.63). We simulated measles outbreaks on the network and evaluated the model against empirical measles data per postcode area from a large outbreak in 2013 (2,766 cases). We found that the network-based model could reproduce the observed size and spatial distribution of the historic outbreak much more clearly than the alternative models, with a case weighted receiver operating characteristic (ROC) sensitivity of 0.94, compared to 0.17 and 0.26 for models that do not account for specific network structure or school-level vaccine uptake, respectively. The key limitation of our framework is that it neglects transmission routes outside of school and household contexts. CONCLUSIONS: Our framework indicates that clustering of unvaccinated children in primary schools connected by unvaccinated children in related secondary schools lead to large, connected clusters of unvaccinated children. Using our approach, we could explain historical outbreaks on a spatial level. Our framework could be further developed to aid future outbreak response.

Higher-valency pneumococcal conjugate vaccines in older adults, taking into account indirect effects from childhood vaccination: a cost-effectiveness study for the Netherlands.

BACKGROUND: New 15- and 20-valent pneumococcal vaccines (PCV15, PCV20) are available for both children and adults, while PCV21 for adults is in development. However, their cost-effectiveness for older adults, taking into account indirect protection and serotype replacement from a switch to PCV15 and PCV20 in childhood vaccination, remains unexamined. METHODS: We used a static model for the Netherlands to assess the cost-effectiveness of different strategies with 23-valent pneumococcal polysaccharide vaccine (PPV23), PCV15, PCV20, and PCV21 for a 65-year-old cohort from a societal perspective, over a 15-year time horizon. Childhood vaccination was varied from PCV10 to PCV13, PCV15, and PCV20. Indirect protection was assumed to reduce the incidence of vaccine serotypes in older adults by 80% (except for serotype 3, no effect), completely offset by an increase in non-vaccine serotype incidence due to serotype replacement. RESULTS: Indirect effects from childhood vaccination reduced the cost-effectiveness of vaccination of older adults, depending on the serotype overlap between the vaccines. With PCV10, PCV13, or PCV15 in children, PCV20 was more effective and less costly for older adults than PPV23 and PCV15. PCV20 costs approximately €10,000 per quality-adjusted life year (QALY) gained compared to no pneumococcal vaccination, which falls below the conventional Dutch €20,000/QALY gained threshold. However, with PCV20 in children, PCV20 was no longer considered cost-effective for older adults, costing €22,550/QALY gained. As indirect effects progressed over time, the cost-effectiveness of PCV20 for older adults further diminished for newly vaccinated cohorts. PPV23 was more cost-effective than PCV20 for cohorts vaccinated 3 years after the switch to PCV20 in children. PCV21 offered the most QALY gains, and its cost-effectiveness was minimally affected by indirect effects due to its coverage of 11 different serotypes compared to PCV20. CONCLUSIONS: For long-term cost-effectiveness in the Netherlands, the pneumococcal vaccine for older adults should either include invasive serotypes not covered by childhood vaccination or become more affordable than its current pricing for individual use.

Repetition in social contacts: implications in modelling the transmission of respiratory infectious diseases in pre-pandemic and pandemic settings.

The spread of viral respiratory infections is intricately linked to human interactions, and this relationship can be characterized and modelled using social contact data. However, many analyses tend to overlook the recurrent nature of these contacts. To bridge this gap, we undertake the task of describing individuals' contact patterns over time by characterizing the interactions made with distinct individuals during a week. Moreover, we gauge the implications of this temporal reconstruction on disease transmission by juxtaposing it with the assumption of random mixing over time. This involves the development of an age-structured individual-based model, using social contact data from a pre-pandemic scenario (the POLYMOD study) and a pandemic setting (the Belgian CoMix study), respectively. We found that accounting for the frequency of contacts impacts the number of new, distinct, contacts, revealing a lower total count than a naive approach, where contact repetition is neglected. As a consequence, failing to account for the repetition of contacts can result in an underestimation of the transmission probability given a contact, potentially leading to inaccurate conclusions when using mathematical models for disease control. We, therefore, underscore the necessity of acknowledging contact repetition when formulating effective public health strategies.

Projecting COVID-19 intensive care admissions for policy advice, the Netherlands, February 2020 to January 2021.

BackgroundModel projections of coronavirus disease 2019 (COVID-19) incidence help policymakers about decisions to implement or lift control measures. During the pandemic, policymakers in the Netherlands were informed on a weekly basis with short-term projections of COVID-19 intensive care unit (ICU) admissions.AimWe aimed at developing a model on ICU admissions and updating a procedure for informing policymakers.MethodThe projections were produced using an age-structured transmission model. A consistent, incremental update procedure integrating all new surveillance and hospital data was conducted weekly. First, up-to-date estimates for most parameter values were obtained through re-analysis of all data sources. Then, estimates were made for changes in the age-specific contact rates in response to policy changes. Finally, a piecewise constant transmission rate was estimated by fitting the model to reported daily ICU admissions, with a changepoint analysis guided by Akaike's Information Criterion.ResultsThe model and update procedure allowed us to make weekly projections. Most 3-week prediction intervals were accurate in covering the later observed numbers of ICU admissions. When projections were too high in March and August 2020 or too low in November 2020, the estimated effectiveness of the policy changes was adequately adapted in the changepoint analysis based on the natural accumulation of incoming data.ConclusionThe model incorporates basic epidemiological principles and most model parameters were estimated per data source. Therefore, it had potential to be adapted to a more complex epidemiological situation with the rise of new variants and the start of vaccination.

Wastewater-based reproduction numbers and projections of COVID-19 cases in three areas in Japan, November 2021 to December 2022.

BackgroundWastewater surveillance has expanded globally as a means to monitor spread of infectious diseases. An inherent challenge is substantial noise and bias in wastewater data because of the sampling and quantification process, limiting the applicability of wastewater surveillance as a monitoring tool.AimTo present an analytical framework for capturing the growth trend of circulating infections from wastewater data and conducting scenario analyses to guide policy decisions.MethodsWe developed a mathematical model for translating the observed SARS-CoV-2 viral load in wastewater into effective reproduction numbers. We used an extended Kalman filter to infer underlying transmissions by smoothing out observational noise. We also illustrated the impact of different countermeasures such as expanded vaccinations and non-pharmaceutical interventions on the projected number of cases using three study areas in Japan during 2021-22 as an example.ResultsObserved notified cases were matched with the range of cases estimated by our approach with wastewater data only, across different study areas and virus quantification methods, especially when the disease prevalence was high. Estimated reproduction numbers derived from wastewater data were consistent with notification-based reproduction numbers. Our projections showed that a 10-20% increase in vaccination coverage or a 10% reduction in contact rate may suffice to initiate a declining trend in study areas.ConclusionOur study demonstrates how wastewater data can be used to track reproduction numbers and perform scenario modelling to inform policy decisions. The proposed framework complements conventional clinical surveillance, especially when reliable and timely epidemiological data are not available.

Contact behaviour before, during and after the COVID-19 pandemic in the Netherlands: evidence from contact surveys, 2016 to 2017 and 2020 to 2023.

BackgroundThe first wave of the COVID-19 pandemic in 2020 was largely mitigated by limiting contacts in the general population. In early 2022, most contact-reducing measures were lifted.AimTo assess whether the population has reverted to pre-pandemic contact behaviour and how this would affect transmission potential of a newly emerging pathogen.MethodsWe compared two studies on contact behaviour in the Netherlands: the PIENTER Corona study, conducted during and after the pandemic (held every 2-6 months from April 2020) and the PIENTER3 study (2016-17, as pre-pandemic baseline). In both, participants (ages 1-85 years) reported number and age group of all face-to-face persons contacted on the previous day in a survey. Transmission potential was examined using the next-generation matrix approach.ResultsWe found an average of 15.4 (95% CI: 14.3-16.4) community contacts per person per day after the pandemic in May 2023, 13% lower than baseline (17.8; 95% CI: 17.0-18.5). Among all ages, children (5-9 years) had the highest number of contacts, both pre- and post-pandemic. Mainly adults aged 20-59 years had not reverted to pre-pandemic behaviours, possibly because they more often work from home. Although the number of contacts is lower compared to the pre-pandemic period, the effect on transmission potential of a newly emerging respiratory pathogen is limited if all age groups were equally susceptible.ConclusionContinuous monitoring of contacts can signal changes in contact patterns and can define a 'new normal' baseline. Both aspects are needed to prepare for a future pandemic.

The fading of the mpox outbreak among men who have sex with men: a mathematical modelling study.

BACKGROUND: In the Netherlands, the number of mpox cases started declining before mpox vaccination was initiated. Most cases were men who have sex with men (MSM). We investigated whether the decline in mpox could be attributed to infection-induced immunity or behavioural adaptations. METHODS: We developed a transmission model and accounted for possible behavioural adaptations: less casual partners and shorter time until MSM with mpox refrain from sexual contacts. RESULTS: Without behavioural adaptations, the peak in modelled cases matched observations, but the decline was less steep than observed. With behavioural adaptations in the model, we found a decline of 16-18% in numbers of casual partners in June and 13-22% in July 2022. Model results showed a halving of the time before refraining from sex. When mpox vaccination started, 57% of MSM with very high sexual activity in the model had been infected. Model scenarios revealed that the outbreak could have waned by November 2022 even without vaccination. CONCLUSIONS: The limited duration of the mpox outbreak in the Netherlands can be ascribed primarily to infection-induced immunity among MSM with high sexual activity levels. The decline was accelerated by behavioural adaptations. Immunity among those most sexually active is essential to impede mpox resurgence.

EpiLPS: A fast and flexible Bayesian tool for estimation of the time-varying reproduction number.

In infectious disease epidemiology, the instantaneous reproduction number [Formula: see text] is a time-varying parameter defined as the average number of secondary infections generated by an infected individual at time t. It is therefore a crucial epidemiological statistic that assists public health decision makers in the management of an epidemic. We present a new Bayesian tool (EpiLPS) for robust estimation of the time-varying reproduction number. The proposed methodology smooths the epidemic curve and allows to obtain (approximate) point estimates and credible intervals of [Formula: see text] by employing the renewal equation, using Bayesian P-splines coupled with Laplace approximations of the conditional posterior of the spline vector. Two alternative approaches for inference are presented: (1) an approach based on a maximum a posteriori argument for the model hyperparameters, delivering estimates of [Formula: see text] in only a few seconds; and (2) an approach based on a Markov chain Monte Carlo (MCMC) scheme with underlying Langevin dynamics for efficient sampling of the posterior target distribution. Case counts per unit of time are assumed to follow a negative binomial distribution to account for potential overdispersion in the data that would not be captured by a classic Poisson model. Furthermore, after smoothing the epidemic curve, a "plug-in'' estimate of the reproduction number can be obtained from the renewal equation yielding a closed form expression of [Formula: see text] as a function of the spline parameters. The approach is extremely fast and free of arbitrary smoothing assumptions. EpiLPS is applied on data of SARS-CoV-1 in Hong-Kong (2003), influenza A H1N1 (2009) in the USA and on the SARS-CoV-2 pandemic (2020-2021) for Belgium, Portugal, Denmark and France.

Contact patterns of older adults with and without frailty in the Netherlands during the COVID-19 pandemic.

BACKGROUND: During the COVID-19 pandemic, social distancing measures were imposed to protect the population from exposure, especially older adults and people with frailty, who have the highest risk for severe outcomes. These restrictions greatly reduced contacts in the general population, but little was known about behaviour changes among older adults and people with frailty themselves. Our aim was to quantify how COVID-19 measures affected the contact behaviour of older adults and how this differed between older adults with and without frailty. METHODS: In 2021, a contact survey was carried out among people aged 70 years and older in the Netherlands. A random sample of persons per age group (70-74, 75-79, 80-84, 85-89, and 90 +) and gender was invited to participate, either during a period with stringent (April 2021) or moderate (October 2021) measures. Participants provided general information on themselves, including their frailty, and they reported characteristics of all persons with whom they had face-to-face contact on a given day over the course of a full week. RESULTS: In total, 720 community-dwelling older adults were included (overall response rate of 15%), who reported 16,505 contacts. During the survey period with moderate measures, participants without frailty had significantly more contacts outside their household than participants with frailty. Especially for females, frailty was a more informative predictor of the number of contacts than age. During the survey period with stringent measures, participants with and without frailty had significantly lower numbers of contacts compared to the survey period with moderate measures. The reduction of the number of contacts was largest for the eldest participants without frailty. As they interact mostly with adults of a similar high age who are likely frail, this reduction of the number of contacts indirectly protects older adults with frailty from SARS-CoV-2 exposure. CONCLUSIONS: The results of this study reveal that social distancing measures during the COVID-19 pandemic differentially affected the contact patterns of older adults with and without frailty. The reduction of contacts may have led to the direct protection of older adults in general but also to the indirect protection of older adults with frailty.

Estimated incubation period distributions of mpox using cases from two international European festivals and outbreaks in a club in Berlin, May to June 2022.

BackgroundSince May 2022, an mpox outbreak affecting primarily men who have sex with men (MSM) has occurred in numerous non-endemic countries worldwide. As MSM frequently reported multiple sexual encounters in this outbreak, reliably determining the time of infection is difficult; consequently, estimation of the incubation period is challenging.AimWe aimed to provide valid and precise estimates of the incubation period distribution of mpox by using cases associated with early outbreak settings where infection likely occurred.MethodsColleagues in European countries were invited to provide information on exposure intervals and date of symptom onset for mpox cases who attended a fetish festival in Antwerp, Belgium, a gay pride festival in Gran Canaria, Spain or a particular club in Berlin, Germany, where early mpox outbreaks occurred. Cases of these outbreaks were pooled; doubly censored models using the log-normal, Weibull and Gamma distributions were fitted to estimate the incubation period distribution.ResultsWe included data on 122 laboratory-confirmed cases from 10 European countries. Depending on the distribution used, the median incubation period ranged between 8 and 9 days, with 5th and 95th percentiles ranging from 2 to 3 and from 20 to 23 days, respectively. The shortest interval that included 50% of incubation periods spanned 8 days (4-11 days).ConclusionCurrent public health management of close contacts should consider that in approximately 5% of cases, the incubation period exceeds the commonly used monitoring period of 21 days.

Mpox outbreak in the Netherlands, 2022: public health response, characteristics of the first 1,000 cases and protection of the first-generation smallpox vaccine.

In early May 2022, a global outbreak of mpox started among persons without travel history to regions known to be enzootic for monkeypox virus (MPXV). On 8 August 2022, the Netherlands reported its 1,000th mpox case, representing a cumulative incidence of 55 per million population, one of the highest cumulative incidences worldwide. We describe characteristics of the first 1,000 mpox cases in the Netherlands, reported between 20 May and 8 August 2022, within the context of the public health response. These cases were predominantly men who have sex with men aged 31-45 years. The vast majority of infections were acquired through sexual contact with casual partners in private or recreational settings including LGBTQIA+ venues in the Netherlands. This indicates that, although some larger upsurges occurred from point-source and/or travel-related events, the outbreak was mainly characterised by sustained transmission within the Netherlands. In addition, we estimated the protective effect of first-generation smallpox vaccine against moderate/severe mpox and found a vaccine effectiveness of 58% (95% CI: 17-78%), suggesting moderate protection against moderate/severe mpox symptoms on top of any possible protection by this vaccine against MPXV infection and disease. Communication with and supporting the at-risk population in following mitigation measures remains essential.

How the Disruption in Sexually Transmitted Infection Care Due to the COVID-19 Pandemic Could Lead to Increased Sexually Transmitted Infection Transmission Among Men Who Have Sex With Men in The Netherlands: A Mathematical Modeling Study.

BACKGROUND: During the COVID-19 pandemic, the disruption in care for sexually transmitted infections (STIs) and the social distancing measures have led to reductions in STI testing and sexual behavior. We assessed the impact of these COVID-19-related changes on transmission of Chlamydia trachomatis (CT) and Neisseria gonorrhoeae (NG) among men who have sex with men (MSM) in The Netherlands. METHODS: We developed a mathematical model for CT and NG transmission among MSM, accounting for COVID-19-related changes in sexual behavior and testing in 2020 to 2021. Changes in 2020 were estimated from data from the Dutch COVID-19, Sex, and Intimacy Survey among MSM and the National Database of STI Clinics. Because of the lack of data for 2021, we examined several scenarios covering a range of changes. RESULTS: A reduction of 10% and 40% in STI testing of symptomatic and asymptomatic, respectively, individuals with a 10% to 20% reduction in numbers of casual partners (according to partner status and activity level) during the second lockdown, resulted in a 2.4% increase in CT prevalence, but a 2.8% decline in NG prevalence in 2021. A 5% and 30% reduction in STI testing of symptomatic and asymptomatic, respectively, individuals with the same reduction in casual partners resulted in a 0.6% increase in CT prevalence and a 4.9% decrease in NG prevalence in 2021. CONCLUSIONS: The disruption in STI care due to COVID-19 might have resulted in a small increase in CT prevalence, but a decrease in NG prevalence. Scaling up STI care is imperative to prevent increases in STI transmission.

Optimal vaccine allocation for COVID-19 in the Netherlands: A data-driven prioritization.

For the control of COVID-19, vaccination programmes provide a long-term solution. The amount of available vaccines is often limited, and thus it is crucial to determine the allocation strategy. While mathematical modelling approaches have been used to find an optimal distribution of vaccines, there is an excessively large number of possible allocation schemes to be simulated. Here, we propose an algorithm to find a near-optimal allocation scheme given an intervention objective such as minimization of new infections, hospitalizations, or deaths, where multiple vaccines are available. The proposed principle for allocating vaccines is to target subgroups with the largest reduction in the outcome of interest. We use an approximation method to reconstruct the age-specific transmission intensity (the next generation matrix), and express the expected impact of vaccinating each subgroup in terms of the observed incidence of infection and force of infection. The proposed approach is firstly evaluated with a simulated epidemic and then applied to the epidemiological data on COVID-19 in the Netherlands. Our results reveal how the optimal allocation depends on the objective of infection control. In the case of COVID-19, if we wish to minimize deaths, the optimal allocation strategy is not efficient for minimizing other outcomes, such as infections. In simulated epidemics, an allocation strategy optimized for an outcome outperforms other strategies such as the allocation from young to old, from old to young, and at random. Our simulations clarify that the current policy in the Netherlands (i.e., allocation from old to young) was concordant with the allocation scheme that minimizes deaths. The proposed method provides an optimal allocation scheme, given routine surveillance data that reflect ongoing transmissions. This approach to allocation is useful for providing plausible simulation scenarios for complex models, which give a more robust basis to determine intervention strategies.

Correction: Practical considerations for measuring the effective reproductive number, Rt.

[This corrects the article DOI: 10.1371/journal.pcbi.1008409.].

Social clustering of unvaccinated children in schools in the Netherlands.

For the measles-mumps-rubella (MMR) vaccine, the World Health Organization-recommended coverage for herd protection is 95% for measles and 80% for rubella and mumps. However, a national vaccine coverage does not reflect social clustering of unvaccinated children, e.g. in schools of Orthodox Protestant or Anthroposophic identity in The Netherlands. To fully characterise this clustering, we estimated one-dose MMR vaccination coverages at all schools in the Netherlands. By combining postcode catchment areas of schools and school feeder data, each child in the Netherlands was characterised by residential postcode, primary and secondary school (referred to as school career). Postcode-level vaccination data were used to estimate vaccination coverages per school career. These were translated to coverages per school, stratified by school identity. Most schools had vaccine coverages over 99%, but major exceptions were Orthodox Protestant schools (63% in primary and 58% in secondary schools) and Anthroposophic schools (67% and 78%). School-level vaccine coverage estimates reveal strong clustering of unvaccinated children. The school feeder data reveal strongly connected Orthodox Protestant and Anthroposophic communities, but separated from one another. This suggests that even at a national one-dose MMR coverage of 97.5%, thousands of children per cohort are not protected by herd immunity.

The estimated disease burden of acute COVID-19 in the Netherlands in 2020, in disability-adjusted life-years.

The impact of COVID-19 on population health is recognised as being substantial, yet few studies have attempted to quantify to what extent infection causes mild or moderate symptoms only, requires hospital and/or ICU admission, results in prolonged and chronic illness, or leads to premature death. We aimed to quantify the total disease burden of acute COVID-19 in the Netherlands in 2020 using the disability-adjusted life-years (DALY) measure, and to investigate how burden varies between age-groups and occupations. Using standard methods and diverse data sources (mandatory notifications, population-level seroprevalence, hospital and ICU admissions, registered COVID-19 deaths, and the literature), we estimated years of life lost (YLL), years lived with disability, DALY and DALY per 100,000 population due to COVID-19, excluding post-acute sequelae, stratified by 5-year age-group and occupation category. The total disease burden due to acute COVID-19 was 286,100 (95% CI: 281,700-290,500) DALY, and the per-capita burden was 1640 (95% CI: 1620-1670) DALY/100,000, of which 99.4% consisted of YLL. The per-capita burden increased steeply with age, starting from 60 to 64 years, with relatively little burden estimated for persons under 50 years old. SARS-CoV-2 infection and associated premature mortality was responsible for a considerable direct health burden in the Netherlands, despite extensive public health measures. DALY were much higher than for other high-burden infectious diseases, but lower than estimated for coronary heart disease. These findings are valuable for informing public health decision-makers regarding the expected COVID-19 health burden among population subgroups, and the possible gains from targeted preventative interventions.