Accumulation of Immunity in Heavy-Tailed Sexual Contact Networks Shapes Mpox Outbreak Sizes.

Many countries affected by the global outbreak of mpox in 2022 have observed a decline in cases. Our mathematical model accounting for heavy-tailed sexual partnership distributions suggests that mpox epidemics can hit the infection-derived herd immunity threshold and begin to decline, with <1% of sexually active men who have sex with men infected regardless of interventions or behavioral changes. We consistently found that many countries and US states experienced an epidemic peak, with cumulative cases of around 0.1% to 0.5% among men who have sex with men. The observed decline in cases may not necessarily be attributable to interventions or behavioral changes primarily.

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.

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.

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.

Patterns of SARS-CoV-2 circulation revealed by a nationwide sewage surveillance programme, the Netherlands, August 2020 to February 2022.

BackgroundSurveillance of SARS-CoV-2 in wastewater offers a near real-time tool to track circulation of SARS-CoV-2 at a local scale. However, individual measurements of SARS-CoV-2 in sewage are noisy, inherently variable and can be left-censored.AimWe aimed to infer latent virus loads in a comprehensive sewage surveillance programme that includes all sewage treatment plants (STPs) in the Netherlands and covers 99.6% of the Dutch population.MethodsWe applied a multilevel Bayesian penalised spline model to estimate time- and STP-specific virus loads based on water flow-adjusted SARS-CoV-2 qRT-PCR data for one to four sewage samples per week for each of the more than 300 STPs.ResultsThe model captured the epidemic upsurges and downturns in the Netherlands, despite substantial day-to-day variation in the measurements. Estimated STP virus loads varied by more than two orders of magnitude, from ca 1012 virus particles per 100,000 persons per day in the epidemic trough in August 2020 to almost 1015 per 100,000 in many STPs in January 2022. The timing of epidemics at the local level was slightly shifted between STPs and municipalities, which resulted in less pronounced peaks and troughs at the national level.ConclusionAlthough substantial day-to-day variation is observed in virus load measurements, wastewater-based surveillance of SARS-CoV-2 that is performed at high sampling frequency can track long-term progression of an epidemic at a local scale in near real time.

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.

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.

Estimated incubation period for monkeypox cases confirmed in the Netherlands, May 2022.

In May 2022, monkeypox outbreaks have been reported in countries not endemic for monkeypox. We estimated the monkeypox incubation period, using reported exposure and symptom-onset times for 18 cases detected and confirmed in the Netherlands up to 31 May 2022. Mean incubation period was 9.0 [corrected] days (5th-95th percentiles: 4.2-17.3), underpinning the current recommendation to monitor or isolate/quarantine case contacts for 21 days. However, as the incubation period may differ between different transmission routes, further epidemiological investigations are needed.

A scenario modelling analysis to anticipate the impact of COVID-19 vaccination in adolescents and children on disease outcomes in the Netherlands, summer 2021.

BackgroundSince the roll-out of COVID-19 vaccines in late 2020 and throughout 2021, European governments have relied on mathematical modelling to inform policy decisions about COVID-19 vaccination.AimWe present a scenario-based modelling analysis in the Netherlands during summer 2021, to inform whether to extend vaccination to adolescents (12-17-year-olds) and children (5-11-year-olds).MethodsWe developed a deterministic, age-structured susceptible-exposed-infectious-recovered (SEIR) model and compared modelled incidences of infections, hospital and intensive care admissions, and deaths per 100,000 people across vaccination scenarios, before the emergence of the Omicron variant.ResultsOur model projections showed that, on average, upon the release of all non-pharmaceutical control measures on 1 November 2021, a large COVID-19 wave may occur in winter 2021/22, followed by a smaller, second wave in spring 2022, regardless of the vaccination scenario. The model projected reductions in infections/severe disease outcomes when vaccination was extended to adolescents and further reductions when vaccination was extended to all people over 5 years-old. When examining projected disease outcomes by age group, individuals benefitting most from extending vaccination were adolescents and children themselves. We also observed reductions in disease outcomes in older age groups, particularly of parent age (30-49 years), when children and adolescents were vaccinated, suggesting some prevention of onward transmission from younger to older age groups.ConclusionsWhile our scenarios could not anticipate the emergence/consequences of SARS-CoV-2 Omicron variant, we illustrate how our approach can assist decision making. This could be useful when considering to provide booster doses or intervening against future infection waves.

Impact of the Heterogeneity in Free Chlorine, UV254, and Ozone Susceptibilities Among Coxsackievirus B5 on the Prediction of the Overall Inactivation Efficiency.

The disinfection susceptibilities of viruses vary even among variants, yet the inactivation efficiency of a certain virus genotype, species, or genus was determined based on the susceptibility of its laboratory strain. The objectives were to evaluate the variability in susceptibilities to free chlorine, UV254, and ozone among 13 variants of coxsackievirus B5 (CVB5) and develop the model allowing for predicting the overall inactivation of heterogeneous CVB5. Our results showed that the susceptibilities differed by up to 3.4-fold, 1.3-fold, and 1.8-fold in free chlorine, UV254, and ozone, respectively. CVB5 in genogroup B exhibited significantly lower susceptibility to free chlorine and ozone than genogroup A, where the laboratory strain, Faulkner, belongs. The capsid protein in genogroup B contained a lower number of sulfur-containing amino acids, readily reactive to oxidants. We reformulated the Chick-Watson model by incorporating the probability distributions of inactivation rate constants to capture the heterogeneity. This expanded Chick-Watson model indicated that up to 4.2-fold larger free chlorine CT is required to achieve 6-log inactivation of CVB5 than the prediction by the Faulkner strain. Therefore, it is recommended to incorporate the variation in disinfection susceptibilities for predicting the overall inactivation of a certain type of viruses.

Estimating the asymptomatic proportion of SARS-CoV-2 infection in the general population: Analysis of nationwide serosurvey data in the Netherlands.

BACKGROUND: The proportion of SARS-CoV-2 positive persons who are asymptomatic-and whether this proportion is age-dependent-are still open research questions. Because an unknown proportion of reported symptoms among SARS-CoV-2 positives will be attributable to another infection or affliction, the observed, or 'crude' proportion without symptoms may underestimate the proportion of persons without symptoms that are caused by SARS-CoV-2 infection. METHODS: Based on two rounds of a large population-based serological study comprising test results on seropositivity and self-reported symptom history conducted in April/May and June/July 2020 in the Netherlands (n = 7517), we estimated the proportion of reported symptoms among those persons infected with SARS-CoV-2 that is attributable to this infection, where the set of relevant symptoms fulfills the ECDC case definition of COVID-19, using inferential methods for the attributable risk (AR). Generalised additive regression modelling was used to estimate the age-dependent relative risk (RR) of reported symptoms, and the AR and asymptomatic proportion (AP) were calculated from the fitted RR. RESULTS: Using age-aggregated data, the 'crude' AP was 37% but the model-estimated AP was 65% (95% CI 63-68%). The estimated AP varied with age, from 74% (95% CI 65-90%) for < 20 years, to 61% (95% CI 57-65%) for the 50-59 years age-group. CONCLUSION: Whereas the 'crude' AP represents a lower bound for the proportion of persons infected with SARS-CoV-2 without COVID-19 symptoms, the AP as estimated via an attributable risk approach represents an upper bound. Age-specific AP estimates can inform the implementation of public health actions such as targetted virological testing and therefore enhance containment strategies.

Estimating the Asymptomatic Ratio of Norovirus Infection During Foodborne Outbreaks With Laboratory Testing in Japan.

BACKGROUND: Foodborne norovirus outbreak data in Japan from 2005-2006, involving virological surveillance of all symptomatic and asymptomatic individuals, were reanalyzed to estimate the asymptomatic ratio of norovirus infection along with the risk of infection and the probability of virus shedding. METHODS: Employing a statistical model that is considered to capture the data-generating process of the outbreak and virus surveillance, maximum likelihood estimation of the asymptomatic ratio was implemented. RESULTS: Assuming that all norovirus outbreaks (n = 55) were the result of random sampling from an identical distribution and ignoring genogroup and genotype specificities, the asymptomatic ratio was estimated at 32.1% (95% confidence interval [CI], 27.7-36.7). Although not significant, separate estimation of the asymptomatic ratio of the GII.4 genotype appeared to be greater than other genotypes and was estimated at 40.7% (95% CI, 32.8-49.0). CONCLUSION: The present study offered the first explicit empirical estimates of the asymptomatic ratio of norovirus infection in natural infection settings. The estimate of about 30% was consistent with those derived from volunteer challenge studies. Practical difficulty in controlling GII.4 outbreaks was supported by the data, considering that a large estimate of the asymptomatic ratio was obtained for the GII.4 genotype.

Dynamics of the pneumonic plague epidemic in Madagascar, August to October 2017.

Transmission potential and severity of pneumonic plague in Madagascar were assessed. Accounting for reporting delay, the reproduction number was estimated at 1.73. The case fatality risk was estimated as 5.5%. Expected numbers of exported cases from Madagascar were estimated across the world and all estimates were below 1 person from August to October, 2017.

Comparative assessment of primary and secondary infection risks in a norovirus outbreak using a household model simulation.

Diarrheal diseases can be transmitted via both primary infection due to exposures to contaminated materials from the environment and secondary infection due to person-to-person contacts. Usually, the importance of secondary infection is empirically assessed by fitting mathematical models to the epidemic curves. However, these empirical models may not be applicable to other epidemic cases because they are developed only for the target epidemics and they don't consider the detail routes of infection. In our previous study, we developed a theoretical model taking into account the various routes of infection that commonly occur in households (e.g., shaking hands, food handling, and changing diapers). This model was made flexible and applicable to any epidemics by means of adjusting model parameters. In this study, we proposed a new index "Vulnerability indicator to secondary infection (VISI)", which expressed a ratio of secondary infection to primary infection risks and calculated this index in a simulated norovirus (NoV) epidemic that involved 10,000 households. The results demonstrated that households composed of more than three members including infant(s) had much higher levels of VISI (5-45) than two-member-households with VISI (0.1-4). These results concluded that the infants were likely to be a hub of secondary infections in highly dense families and therefore careful handling of diapers was deemed indispensible in such families to effectively control the secondary infections.

No detectable fitness cost of infection by cell-fusing agent virus in Aedes aegypti mosquitoes.

Aedes mosquitoes are well-known vectors of arthropod-borne viruses (arboviruses). Mosquitoes are more frequently infected with insect-specific viruses (ISVs) that cannot infect vertebrates. Some ISVs interfere with arbovirus replication in mosquito vectors, which has gained attention for potential use against arbovirus transmission. Cell-fusing agent virus (CFAV), a widespread ISV, can reduce arbovirus dissemination in Ae. aegypti. However, vectorial capacity is largely governed by other parameters than pathogen load, including mosquito survival and biting behaviour. Understanding how ISVs impact these mosquito fitness-related traits is critical to assess the potential risk of using ISVs as biological agents. Here, we examined the effects of CFAV infection on Ae. aegypti mosquito fitness. We found no significant reduction in mosquito survival, blood-feeding behaviour and reproduction, suggesting that Ae. aegypti is tolerant to CFAV. The only detectable effect was a slight increase in human attraction of CFAV-infected females in one out of eight trials. Viral tolerance is beneficial for introducing CFAV into natural mosquito populations, whereas the potential increase in biting activity must be further investigated. Our results provide the first insight into the link between ISVs and Aedes mosquito fitness and highlight the importance of considering all aspects of vectorial capacity for arbovirus control using ISVs.

Genotype-dependent kinetics of enterovirus inactivation by free chlorine and ultraviolet (UV) irradiation.

Inactivation kinetics of enterovirus by disinfection is often studied using a single laboratory strain of a given genotype. Environmental variants of enterovirus are genetically distinct from the corresponding laboratory strain, yet it is poorly understood how these genetic differences affect inactivation. Here we evaluated the inactivation kinetics of nine coxsackievirus B3 (CVB3), ten coxsackievirus B4 (CVB4), and two echovirus 11 (E11) variants by free chlorine and ultraviolet irradiation (UV). The inactivation kinetics by free chlorine were genotype- (i.e., susceptibility: CVB5 < CVB3 ≈ CVB4 < E11) and genogroup-dependent and exhibited up to 15-fold difference among the tested viruses. In contrast, only minor (up to 1.3-fold) differences were observed in the UV inactivation kinetics. The differences in variability between the two disinfectants could be rationalized by their respective inactivation mechanisms: inactivation by UV mainly depends on the genomic size and composition, which was similar for all viruses tested, whereas free chlorine targets the viral capsid protein, which exhibited critical differences between genogroups and genotypes. Finally, we integrated the observed variability in inactivation rate constants into an expanded Chick-Watson model to estimate the overall inactivation of an enterovirus consortium. The results highlight that the distribution of inactivation rate constants and the abundance of each genotype are essential parameters to accurately predict the overall inactivation of an enterovirus population by free chlorine. We conclude that predictions based on inactivation data of a single variant or reference pathogen alone likely overestimate the true disinfection efficiency of free chlorine.

Duration of SARS-CoV-2 viral shedding in faeces as a parameter for wastewater-based epidemiology: Re-analysis of patient data using a shedding dynamics model.

Wastewater-based epidemiology (WBE) is one of the most promising approaches to effectively monitor the spread of COVID-19. The virus concentration in faeces and its temporal variations are essential information for WBE. While some clinical studies have reported SARS-CoV-2 concentrations in faeces, the value varies amongst patients and changes over time. The present study aimed to examine how the temporal variations in the concentration of virus in faeces affect the monitoring of disease incidence. We reanalysed the experimental findings of clinical studies to estimate the duration of virus shedding and the faecal virus concentration. Available experimental data as of 23 October 2020 were collected. The viral shedding kinetics was modelled, and the dynamic model was fitted to the collected data by a Bayesian framework. Using posterior distributions, the duration of viral shedding and the concentration of virus copies in faeces over time were computed. We estimated the median concentration of SARS-CoV-2 in faeces as 3.4 (95% CrI: 0.24-6.5) log copies per gram-faeces over the shedding period, and our model implied that the duration of viral shedding was 26.0 days (95% CrI: 21.7-34.9), given the current standard quantification limit (Ct = 40). With simulated incidences, our results also indicated that a one-week delay between symptom onset and wastewater sampling increased the estimation of incidence by a factor of 17.2 (i.e., 101.24 times higher). Our results demonstrated that the temporal variation in virus concentration in faeces affects microbial monitoring systems such as WBE. The present study also implied the need for adjusting the estimates of virus concentration in faeces by incorporating the kinetics of unobserved concentrations. The method used in this study is easily implemented in further simulations; therefore, the results of this study might contribute to enhancing disease surveillance and risk assessments that require quantities of virus to be excreted into the environment.

Age-dependent association between SARS-CoV-2 cases reported by passive surveillance and viral load in wastewater.

The actual number of individuals infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is difficult to estimate using a case-reporting system (i.e., passive surveillance) alone because of asymptomatic infection. While wastewater-based epidemiology has been implemented as an alternative/additional monitoring tool to reduce reporting bias, the relationship between passive and wastewater surveillance data has not yet been explicitly examined. As there is strong age dependency in the symptomatic ratio of SARS-CoV-2 infections, here, we aimed to estimate i) an age-dependent association between the number of reported cases and viral load in wastewater and ii) the time lag between these time series. The viral load in wastewater was modeled as a combination of contributions from virus shedding by different age groups, incorporating the delay, and fitted with daily case count data collected from the Massachusetts Department of Public Health and SARS-CoV-2 RNA concentration in wastewater recorded by the Massachusetts Water Resources Authority. The estimated lag between the time series of viral loads in wastewater and of reported cases was 10.8 (95% confidence interval: 10.2-11.6) and 8.8 (8.4-9.1) days for the northern and southern areas of the wastewater treatment plant, respectively. The estimated contribution rate of a reported case to the viral load in wastewater in the 0-19 yr age group was 0.38 (0.35-0.41) and 0.40 (0.37-0.43) for the northern and southern areas, and that in the 80+ yr age group was 0.67 (0.65-0.69) and 0.51 (0.49-0.52) for the northern and southern areas, respectively. The estimated lag between these time series suggested the predictability of reported cases 10 days later using viral loads in wastewater. The contribution of a reported case in passive surveillance to the viral load in wastewater differed by age, suggesting a large variation in viral shedding kinetics among age groups.