Different forms of superspreading lead to different outcomes: heterogeneity in infectiousness and contact behavior relevant for the case of SARS-CoV-2 (preprint)

Superspreading events play an important role in the spread of SARS-CoV-2 and several other pathogens. Hence, while the basic reproduction number of the original Wuhan SARS-CoV-2 is estimated to be about 3 for Belgium, there is substantial inter-individual variation in the number of secondary cases each infected individual causes. Multiple factors contribute to the occurrence of superspreading events: heterogeneity in infectiousness and susceptibility, variations in contact behavior, and the environment in which transmission takes place. While superspreading has been included in several infectious disease transmission models, our understanding of the effect that these different forms of superspreading have on the spread of pathogens and the effectiveness of control measures remains limited. To disentangle the effects of infectiousness-related heterogeneity on the one hand and contact-related heterogeneity on the other, we implemented both forms of superspreading in an individual-based model describing the transmission and spread of SARS-CoV-2 in the Belgian population. We considered its impact on viral spread as well as on the effectiveness of social distancing. We found that the effects of superspreading driven by heterogeneity in infectiousness are very different from the effects of superspreading driven by heterogeneity in contact behavior. On the one hand, a higher level of infectiousness-related heterogeneity results in less outbreaks occurring following the introduction of one infected individual. Outbreaks were also slower, with a lower peak which occurred at a later point in time, and a lower herd immunity threshold. Finally, the risk of resurgence of an outbreak following a period of lockdown decreased. On the other hand, when contact-related heterogeneity was high, this also led to smaller final sizes, but caused outbreaks to be more explosive in regard to other aspects (such as higher peaks which occurred earlier, and a higher herd immunity threshold). Finally, the risk of resurgence of an outbreak following a period of lockdown increased. Determining the contribution of both source of heterogeneity is therefore important but left to be explored further. Author summaryTo investigate the effect of different sources of superspreading on disease dynamics, we implemented superspreading driven by heterogeneity in infectiousness and heterogeneity in contact behavior into an individual-based model for the transmission of SARS-CoV-2 in the Belgian population. We compared the impact of both forms of superspreading in a scenario without interventions as well as in a scenario in which a period of strict social distancing (i.e. a lockdown) is followed by a period of partial release. We found that both forms of superspreading have very different effects. On the one hand, increasing the level of infectiousness-related heterogeneity led to less outbreaks being observed following the introduction of one infected individual in the population. Furthermore, final outbreak sizes decreased, and outbreaks became slower, with lower and later peaks, and a lower herd immunity threshold. Finally, the risk for resurgence of an outbreak following a period of lockdown also decreased. On the other hand, when contact-related heterogeneity was high, this also led to smaller final sizes, but caused outbreaks to be more explosive regarding other aspects (such as higher peaks that occurred earlier). The herd immunity threshold also increased, as did the risk of resurgence of outbreaks.

Inferring age-specific differences in susceptibility to and infectiousness upon SARS-CoV-2 infection based on Belgian social contact data (preprint)

Several important aspects related to SARS-CoV-2 transmission are not well known due to a lack of appropriate data. However, mathematical and computational tools can be used to extract part of this information from the available data, like some hidden age-related characteristics. In this paper, we investigate age-specific differences in susceptibility to and infectiousness upon contracting SARS-CoV-2 infection. More specifically, we use panel-based social contact data from diary-based surveys conducted in Belgium combined with the next generation principle to infer the relative incidence and we compare this to real-life incidence data. Comparing these two allows for the estimation of age-specific transmission parameters. Our analysis implies the susceptibility in children to be around half of the susceptibility in adults, and even lower for very young children (preschooler). However, the probability of adults and the elderly to contract the infection is decreasing throughout the vaccination campaign, thereby modifying the picture over time.

Leveraging of SARS-CoV-2 PCR cycle thresholds values (Ct) to forecast COVID-19 trends (preprint)

IntroductionWe assessed the usefulness of SARS-CoV-2 RT-PCR cycle thresholds (Ct) values trends produced by the LHUB-ULB (a consolidated microbiology laboratory located in Brussels, Belgium) for monitoring the epidemics dynamics at local and national levels and for improving forecasting models. MethodsSARS-CoV-2 RT-PCR Ct values produced from April 1, 2020, to May 15, 2021, were compared with national COVID-19 confirmed cases notifications according to their geographical and time distribution. These Ct values were evaluated against both a phase diagram predicting the number of COVID-19 patients requiring intensive care and an age-structured model estimating COVID-19 prevalence in Belgium. ResultsOver 155,811 RT-PCR performed, 12,799 were positive and 7,910 Ct values were available for analysis. The 14-day median Ct values were negatively correlated with the 14-day mean daily positive tests with a lag of 17 days. In addition, the 14-day mean daily positive tests in LHUB-ULB were strongly correlated with the 14-day mean confirmed cases in the Brussels-Capital and in Belgium with coinciding start, peak and end of the different waves of the epidemic. Ct values decreased concurrently with the forecasted phase-shifts of the diagram. Similarly, the evolution of 14-day median Ct values was negatively correlated with daily estimated prevalence for all age-classes. ConclusionWe provide preliminary evidence that trends of Ct values can help to both follow and predict the epidemics trajectory at local and national levels, underlining that consolidated microbiology laboratories can act as epidemic sensors as they gather data that are representative of the geographical area they serve.

Covid-19 Belgium: Extended SEIR-QD model with nursery homes and long-term scenarios-based forecasts from school opening (preprint)

We model the evolution of the covid-19 epidemic in Belgium with an age-structured extended SEIR-QD epidemic model with separated consideration for nursery homes. All parameters of the model are estimated using a MCMC method, except integrated data on social contacts. The model is calibrated on hospitals' data, number of deaths, nursery homes' tests and serological tests. We present the current situation on September 2020 as well as long-term scenarios-based forecasts with the possibility of a second wave in function of new transmissions from contacts at school.

Covid-19 Belgium: Extended SEIR-QD model with nursing homes and long-term scenarios-based forecasts (preprint)

Following the spread of the COVID-19 pandemic and pending the establishment of vaccination campaigns, several non pharmaceutical interventions such as partial and full lockdown, quarantine and measures of physical distancing have been imposed in order to reduce the spread of the disease and to lift the pressure on healthcare system. Mathematical models are important tools for estimating the impact of these interventions, for monitoring the current evolution of the epidemic at a national level and for estimating the potential long-term consequences of relaxation of measures. In this paper, we model the evolution of the COVID-19 epidemic in Belgium with a deterministic age-structured extended compartmental model. Our model takes special consideration for nursing homes which are modelled as separate entities from the general population in order to capture the specific delay and dynamics within these entities. The model integrates social contact data and is fitted on hospitalisations data (admission and discharge), on the daily number of COVID-19 deaths (with a distinction between general population and nursing home related deaths) and results from serological studies, with a sensitivity analysis based on a Bayesian approach. We present the situation as in November 2020 with the estimation of some characteristics of the COVID-19 deduced from the model. We also present several mid-term and long-term projections based on scenarios of reinforcement or relaxation of social contacts for different general sectors, with a lot of uncertainties remaining.