An SEIR network epidemic model with manual and digital contact tracing allowing delays.

We consider an SEIR epidemic model on a network also allowing random contacts, where recovered individuals could either recover naturally or be diagnosed. Upon diagnosis, manual contact tracing is triggered such that each infected network contact is reported, tested and isolated with some probability and after a random delay. Additionally, digital tracing (based on a tracing app) is triggered if the diagnosed individual is an app-user, and then all of its app-using infectees are immediately notified and isolated. The early phase of the epidemic with manual and/or digital tracing is approximated by different multi-type branching processes, and three respective reproduction numbers are derived. The effectiveness of both contact tracing mechanisms is numerically quantified through the reduction of the reproduction number. This shows that app-using fraction plays an essential role in the overall effectiveness of contact tracing. The relative effectiveness of manual tracing compared to digital tracing increases if: more of the transmission occurs on the network, when the tracing delay is shortened, and when the network degree distribution is heavy-tailed. For realistic values, the combined tracing case can reduce R0 by 20%-30%, so other preventive measures are needed to reduce the reproduction number down to 1.2-1.4 for contact tracing to make it successful in avoiding big outbreaks.

SIRS epidemics with individual heterogeneity of immunity waning.

In the current paper we analyse an extended SIRS epidemic model in which immunity at the individual level wanes gradually at exponential rate, but where the waning rate may differ between individuals, for instance as an effect of differences in immune systems. The model also includes vaccination schemes aimed to reach and maintain herd immunity. We consider both the informed situation where the individual waning parameters are known, thus allowing selection of vaccinees being based on both time since last vaccination as well as on the individual waning rate, and the more likely uninformed situation where individual waning parameters are unobserved, thus only allowing vaccination schemes to depend on time since last vaccination. The optimal vaccination policies for both the informed and uniformed heterogeneous situation are derived and compared with the homogeneous waning model (meaning all individuals have the same immunity waning rate), as well as to the classic SIRS model where immunity at the individual level drops from complete immunity to complete susceptibility in one leap. It is shown that the classic SIRS model requires least vaccines, followed by the SIRS with homogeneous gradual waning, followed by the informed situation for the model with heterogeneous gradual waning. The situation requiring most vaccines for herd immunity is the most likely scenario, that immunity wanes gradually with unobserved individual heterogeneity. For parameter values chosen to mimic COVID-19 and assuming perfect initial immunity and cumulative immunity of 12 months, the classic homogeneous SIRS epidemic suggests that vaccinating individuals every 15 months is sufficient to reach and maintain herd immunity, whereas the uninformed case for exponential waning with rate heterogeneity corresponding to a coefficient of variation being 0.5, requires that individuals instead need to be vaccinated every 4.4 months.

Management of cervical CSF-venous fistula causing acute cognitive impairment and coma.

CSF-venous fistulas (CVFs) are increasingly recognised as a cause of spontaneous intracranial hypotension. They may present atypically including with brain sagging pseudo-dementia. Cervical CVFs are rare and their management can be difficult due to associated eloquent nerve roots. We report the case of a 49-year-old woman who presented with cognitive decline progressing to coma. Brain imaging showed features of spontaneous intracranial hypotension and a right C7 CVF was identified at digital subtraction and CT myelography. Initial treatment with CT-guided injection of fibrin sealant produced temporary improvement in symptoms before surgical treatment resulted in total clinical remission and radiological resolution.

Extending susceptible-infectious-recovered-susceptible epidemics to allow for gradual waning of immunity.

Susceptible-infectious-recovered-susceptible (SIRS) epidemic models assume that individual immunity wanes in one leap, from complete immunity to complete susceptibility. For many diseases immunity on the contrary wanes gradually, something that has become even more evident during COVID-19 pandemic where also recently infected have a reinfection risk, and booster vaccines are given to increase immunity. Here, a novel mathematical model is presented allowing for the gradual decay of immunity following linear or exponential waning functions. The two new models and the SIRS model are compared assuming all three models have the same cumulative immunity. When no intervention is put in place, we find that the long-term prevalence is higher for the models with gradual waning. If aiming for herd immunity by continuous vaccination, it is shown that larger vaccine quantities are required when immunity wanes gradually compared with results obtained from the SIRS model, and this difference is the biggest for the most realistic assumption of exponentially waning of immunity. For parameter choices fitting to COVID-19, the critical amount of vaccine supply is about 50% higher if immunity wanes linearly, and more than 150% higher when immunity wanes exponentially, when compared with the classic SIRS epidemic model.

Bayesian nowcasting with leading indicators applied to COVID-19 fatalities in Sweden.

The real-time analysis of infectious disease surveillance data is essential in obtaining situational awareness about the current dynamics of a major public health event such as the COVID-19 pandemic. This analysis of e.g., time-series of reported cases or fatalities is complicated by reporting delays that lead to under-reporting of the complete number of events for the most recent time points. This can lead to misconceptions by the interpreter, for instance the media or the public, as was the case with the time-series of reported fatalities during the COVID-19 pandemic in Sweden. Nowcasting methods provide real-time estimates of the complete number of events using the incomplete time-series of currently reported events and information about the reporting delays from the past. In this paper we propose a novel Bayesian nowcasting approach applied to COVID-19-related fatalities in Sweden. We incorporate additional information in the form of time-series of number of reported cases and ICU admissions as leading signals. We demonstrate with a retrospective evaluation that the inclusion of ICU admissions as a leading signal improved the nowcasting performance of case fatalities for COVID-19 in Sweden compared to existing methods.

Monitoring real-time transmission heterogeneity from incidence data.

The transmission heterogeneity of an epidemic is associated with a complex mixture of host, pathogen and environmental factors. And it may indicate superspreading events to reduce the efficiency of population-level control measures and to sustain the epidemic over a larger scale and a longer duration. Methods have been proposed to identify significant transmission heterogeneity in historic epidemics based on several data sources, such as contact history, viral genomes and spatial information, which may not be available, and more importantly ignore the temporal trend of transmission heterogeneity. Here we attempted to establish a convenient method to estimate real-time heterogeneity over an epidemic. Within the branching process framework, we introduced an instant-individualheterogenous infectiousness model to jointly characterize the variation in infectiousness both between individuals and among different times. With this model, we could simultaneously estimate the transmission heterogeneity and the reproduction number from incidence time series. We validated the model with data of both simulated and real outbreaks. Our estimates of the overall and real-time heterogeneities of the six epidemics were consistent with those presented in the literature. Additionally, our model is robust to the ubiquitous bias of under-reporting and misspecification of serial interval. By analyzing recent data from South Africa, we found evidence that the Omicron might be of more significant transmission heterogeneity than Delta. Our model based on incidence data was proved to be reliable in estimating the real-time transmission heterogeneity.

Analysing the Effect of Test-and-Trace Strategy in an SIR Epidemic Model.

Consider a Markovian SIR epidemic model in a homogeneous community. To this model we add a rate at which individuals are tested, and once an infectious individual tests positive it is isolated and each of their contacts are traced and tested independently with some fixed probability. If such a traced individual tests positive it is isolated, and the contact tracing is iterated. This model is analysed using large population approximations, both for the early stage of the epidemic when the "to-be-traced components" of the epidemic behaves like a branching process, and for the main stage of the epidemic where the process of to-be-traced components converges to a deterministic process defined by a system of differential equations. These approximations are used to quantify the effect of testing and of contact tracing on the effective reproduction numbers (for the components as well as for the individuals), the probability of a major outbreak, and the final fraction getting infected. Using numerical illustrations when rates of infection and natural recovery are fixed, it is shown that Test-and-Trace strategy is effective in reducing the reproduction number. Surprisingly, the reproduction number for the branching process of components is not monotonically decreasing in the tracing probability, but the individual reproduction number is conjectured to be monotonic as expected. Further, in the situation where individuals also self-report for testing, the tracing probability is more influential than the screening rate (measured by the fraction infected being screened).

Modelling preventive measures and their effect on generation times in emerging epidemics.

We present a stochastic epidemic model to study the effect of various preventive measures, such as uniform reduction of contacts and transmission, vaccination, isolation, screening and contact tracing, on a disease outbreak in a homogeneously mixing community. The model is based on an infectivity process, which we define through stochastic contact and infectiousness processes, so that each individual has an independent infectivity profile. In particular, we monitor variations of the reproduction number and of the distribution of generation times. We show that some interventions, i.e. uniform reduction and vaccination, affect the former while leaving the latter unchanged, whereas other interventions, i.e. isolation, screening and contact tracing, affect both quantities. We provide a theoretical analysis of the variation of these quantities, and we show that, in practice, the variation of the generation time distribution can be significant and that it can cause biases in the estimation of reproduction numbers. The framework, because of its general nature, captures the properties of many infectious diseases, but particular emphasis is on COVID-19, for which numerical results are provided.

The risk for a new COVID-19 wave and how it depends on R 0, the current immunity level and current restrictions.

The COVID-19 pandemic has hit different regions differently. The current disease-induced immunity level î in a region approximately equals the cumulative fraction infected, which primarily depends on two factors: (i) the initial potential for COVID-19 in the region (R 0), and (ii) the preventive measures put in place. Using a mathematical model including heterogeneities owing to age, social activity and susceptibility, and allowing for time-varying preventive measures, the risk for a new epidemic wave and its doubling time are investigated. Focus lies on quantifying the minimal overall effect of preventive measures p Min needed to prevent a future outbreak. It is shown that î plays a more influential roll than when immunity is obtained from vaccination. Secondly, by comparing regions with different R 0 and î it is shown that regions with lower R 0 and low î may need higher preventive measures (p Min) compared with regions having higher R 0 but also higher î, even when such immunity levels are far from herd immunity. Our results are illustrated on different regions but these comparisons contain lots of uncertainty due to simplistic model assumptions and insufficient data fitting, and should accordingly be interpreted with caution.

Quantifying the preventive effect of wearing face masks.

An important task in combating the current Covid-19 pandemic lies in estimating the effect of different preventive measures. Here, we focus on the preventive effect of enforcing the use of face masks. Several publications study this effect, however, often using different measures such as: the relative attack rate in case-control studies, the effect on incidence growth/decline in a specific time frame and the effect on the number of infected in a given time frame. These measures all depend on community-specific features and are hence not easily transferred to other community settings. We argue that a more universal measure is the relative reduction in the reproduction number, which we call the face mask effect, E FM. It is shown how to convert the other measures to E FM. We also apply the methodology to four empirical studies using different effect-measures. When converted to estimates of E FM, all estimates lie between 15 and 40%, suggesting that mandatory face masks reduce the reproduction number by an amount in this range, when compared with no individuals wearing face masks.

Inflow restrictions can prevent epidemics when contact tracing efforts are effective but have limited capacity.

When a region tries to prevent an outbreak of an epidemic, two broad strategies are available: limiting the inflow of infected cases by using travel restrictions and quarantines or limiting the risk of local transmission from imported cases by using contact tracing and other community interventions. A number of papers have used epidemiological models to argue that inflow restrictions are unlikely to be effective. We simulate a simple epidemiological model to show that this conclusion changes if containment efforts such as contact tracing have limited capacity. In particular, our results show that moderate travel restrictions can lead to large reductions in the probability of an epidemic when contact tracing is effective but the contact tracing system is close to being overwhelmed.

Inferring transmission heterogeneity using virus genealogies: Estimation and targeted prevention.

Spread of HIV typically involves uneven transmission patterns where some individuals spread to a large number of individuals while others to only a few or none. Such transmission heterogeneity can impact how fast and how much an epidemic spreads. Further, more efficient interventions may be achieved by taking such transmission heterogeneity into account. To address these issues, we developed two phylogenetic methods based on virus sequence data: 1) to generally detect if significant transmission heterogeneity is present, and 2) to pinpoint where in a phylogeny high-level spread is occurring. We derive inference procedures to estimate model parameters, including the amount of transmission heterogeneity, in a sampled epidemic. We show that it is possible to detect transmission heterogeneity under a wide range of simulated situations, including incomplete sampling, varying levels of heterogeneity, and including within-host genetic diversity. When evaluating real HIV-1 data from different epidemic scenarios, we found a lower level of transmission heterogeneity in slowly spreading situations and a higher level of heterogeneity in data that included a rapid outbreak, while R0 and Sackin's index (overall tree shape statistic) were similar in the two scenarios, suggesting that our new method is able to detect transmission heterogeneity in real data. We then show by simulations that targeted prevention, where we pinpoint high-level spread using a coalescence measurement, is efficient when sequence data are collected in an ongoing surveillance system. Such phylogeny-guided prevention is efficient under both single-step contact tracing as well as iterative contact tracing as compared to random intervention.

Key questions for modelling COVID-19 exit strategies.

Combinations of intense non-pharmaceutical interventions (lockdowns) were introduced worldwide to reduce SARS-CoV-2 transmission. Many governments have begun to implement exit strategies that relax restrictions while attempting to control the risk of a surge in cases. Mathematical modelling has played a central role in guiding interventions, but the challenge of designing optimal exit strategies in the face of ongoing transmission is unprecedented. Here, we report discussions from the Isaac Newton Institute 'Models for an exit strategy' workshop (11-15 May 2020). A diverse community of modellers who are providing evidence to governments worldwide were asked to identify the main questions that, if answered, would allow for more accurate predictions of the effects of different exit strategies. Based on these questions, we propose a roadmap to facilitate the development of reliable models to guide exit strategies. This roadmap requires a global collaborative effort from the scientific community and policymakers, and has three parts: (i) improve estimation of key epidemiological parameters; (ii) understand sources of heterogeneity in populations; and (iii) focus on requirements for data collection, particularly in low-to-middle-income countries. This will provide important information for planning exit strategies that balance socio-economic benefits with public health.

A mathematical model reveals the influence of population heterogeneity on herd immunity to SARS-CoV-2.

Despite various levels of preventive measures, in 2020, many countries have suffered severely from the coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. Using a model, we show that population heterogeneity can affect disease-induced immunity considerably because the proportion of infected individuals in groups with the highest contact rates is greater than that in groups with low contact rates. We estimate that if R 0 = 2.5 in an age-structured community with mixing rates fitted to social activity, then the disease-induced herd immunity level can be ~43%, which is substantially less than the classical herd immunity level of 60% obtained through homogeneous immunization of the population. Our estimates should be interpreted as an illustration of how population heterogeneity affects herd immunity rather than as an exact value or even a best estimate.

Completing Linnaeus's inventory of the Swedish insect fauna: Only 5,000 species left?

Despite more than 250 years of taxonomic research, we still have only a vague idea about the true size and composition of the faunas and floras of the planet. Many biodiversity inventories provide limited insight because they focus on a small taxonomic subsample or a tiny geographic area. Here, we report on the size and composition of the Swedish insect fauna, thought to represent roughly half of the diversity of multicellular life in one of the largest European countries. Our results are based on more than a decade of data from the Swedish Taxonomy Initiative and its massive inventory of the country's insect fauna, the Swedish Malaise Trap Project The fauna is considered one of the best known in the world, but the initiative has nevertheless revealed a surprising amount of hidden diversity: more than 3,000 new species (301 new to science) have been documented so far. Here, we use three independent methods to analyze the true size and composition of the fauna at the family or subfamily level: (1) assessments by experts who have been working on the most poorly known groups in the fauna; (2) estimates based on the proportion of new species discovered in the Malaise trap inventory; and (3) extrapolations based on species abundance and incidence data from the inventory. For the last method, we develop a new estimator, the combined non-parametric estimator, which we show is less sensitive to poor coverage of the species pool than other popular estimators. The three methods converge on similar estimates of the size and composition of the fauna, suggesting that it comprises around 33,000 species. Of those, 8,600 (26%) were unknown at the start of the inventory and 5,000 (15%) still await discovery. We analyze the taxonomic and ecological composition of the estimated fauna, and show that most of the new species belong to Hymenoptera and Diptera groups that are decomposers or parasitoids. Thus, current knowledge of the Swedish insect fauna is strongly biased taxonomically and ecologically, and we show that similar but even stronger biases have distorted our understanding of the fauna in the past. We analyze latitudinal gradients in the size and composition of known European insect faunas and show that several of the patterns contradict the Swedish data, presumably due to similar knowledge biases. Addressing these biases is critical in understanding insect biomes and the ecosystem services they provide. Our results emphasize the need to broaden the taxonomic scope of current insect monitoring efforts, a task that is all the more urgent as recent studies indicate a possible worldwide decline in insect faunas.

Introducing pre-exposure prophylaxis to prevent HIV acquisition among men who have sex with men in Sweden: insights from a mathematical pair formation model.

OBJECTIVES: Since 2017, the Public Health Agency of Sweden recommends that pre-exposure prophylaxis (PrEP) for HIV should be offered to high-risk individuals, in particular to men who have sex with men (MSM). The objective of this study is to develop a mathematical model investigating the effect of introducing PrEP to MSM in Sweden. DESIGN: A pair formation model, including steady and casual sex partners, is developed to study the impact of introducing PrEP. Two groups are included in the model: sexually high active MSM and sexually low active MSM. Three mixing assumptions between the groups are considered. SETTING: A gay-friendly MSM HIV/sexually transmitted infection testing clinic in Stockholm, Sweden. This clinic started offering PrEP to MSM in October 2018. PARTICIPANTS: The model is calibrated according to detailed sexual behaviour data gathered in 2015 among 403 MSM. RESULTS: By targeting sexually high active MSM, a PrEP coverage of 3.5% of the MSM population (10% of all high actives) would result in the long-term HIV prevalence to drop considerably (close to 0%). While targeting only low actives would require a PrEP coverage of 35% for a similar reduction. The main effect of PrEP is the reduced susceptibility, whereas the increased HIV testing rate (every third month) among PrEP users plays a lesser role. CONCLUSIONS: To create a multifaceted picture of the effects of interventions against HIV, we need models that include the different stages of HIV infection and real-world data on detailed sexual behaviour to calibrate the mathematical models. Our findings conclude that targeting HIV high-risk individuals, within HIV risk populations such as MSM, with PrEP programmes could greatly decrease the long-term HIV prevalence in Sweden. Therefore, risk stratification of individuals is of importance in PrEP implementation programmes, to ensure optimising the effect and cost-effectiveness of such programmes.

Model selection and parameter estimation for dynamic epidemic models via iterated filtering: application to rotavirus in Germany.

Despite the wide application of dynamic models in infectious disease epidemiology, the particular modeling of variability in the different model components is often subjective rather than the result of a thorough model selection process. This is in part because inference for a stochastic transmission model can be difficult since the likelihood is often intractable due to partial observability. In this work, we address the question of adequate inclusion of variability by demonstrating a systematic approach for model selection and parameter inference for dynamic epidemic models. For this, we perform inference for six partially observed Markov process models, which assume the same underlying transmission dynamics, but differ with respect to the amount of variability they allow for. The inference framework for the stochastic transmission models is provided by iterated filtering methods, which are readily implemented in the R package pomp by King and others (2016, Statistical inference for partially observed Markov processes via the R package pomp. Journal of Statistical Software69, 1-43). We illustrate our approach on German rotavirus surveillance data from 2001 to 2008, discuss practical difficulties of the methods used and calculate a model based estimate for the basic reproduction number $R_0$ using these data.

Dynamic modeling of hepatitis C transmission among people who inject drugs.

To reach the WHO goal of hepatitis C elimination, it is essential to identify the number of people unaware of their hepatitis C virus (HCV) infection and to investigate the effect of interventions on the disease transmission dynamics. In many high-income countries, one of the primary routes of HCV transmission is via contaminated needles shared by people who inject drugs (PWIDs). However, substantial underreporting combined with high uncertainty regarding the size of this difficult to reach population, makes it challenging to estimate the core indicators recommended by the WHO. To support progress toward the elimination goal, we present a novel multi-layered dynamic transmission model for HCV transmission within a PWID population. The model explicitly accounts for disease stage (acute and chronic), injection drug use status (active and former PWIDs), status of diagnosis (diagnosed and undiagnosed) and country of disease acquisition (domestic or abroad). First, based on this model, and using routine surveillance data, we estimate the number of undiagnosed PWIDs, the true incidence, the average time until diagnosis, the reproduction numbers and associated uncertainties. Second, we examine the impact of two interventions on disease dynamics: (1) direct-acting antiviral drug treatment, and (2) needle exchange programs. As a proof of concept, we illustrate our results for a specific data set. In addition, we develop a web application to allow our model to be explored interactively and with different parameter values.

Getting more from heterogeneous HIV-1 surveillance data in a high immigration country: estimation of incidence and undiagnosed population size using multiple biomarkers.

BACKGROUND: Most HIV infections originate from individuals who are undiagnosed and unaware of their infection. Estimation of this quantity from surveillance data is hard because there is incomplete knowledge about (i) the time between infection and diagnosis (TI) for the general population, and (ii) the time between immigration and diagnosis for foreign-born persons. METHODS: We developed a new statistical method for estimating the incidence of HIV-1 and the number of undiagnosed people living with HIV (PLHIV), based on dynamic modelling of heterogeneous HIV-1 surveillance data. The methods consist of a Bayesian non-linear mixed effects model using multiple biomarkers to estimate TI of HIV-1-positive individuals, and a novel incidence estimator which distinguishes between endogenous and exogenous infections by modelling explicitly the probability that a foreign-born person was infected either before or after immigration. The incidence estimator allows for direct calculation of the number of undiagnosed persons. The new methodology is illustrated combining heterogeneous surveillance data from Sweden between 2003 and 2015. RESULTS: A leave-one-out cross-validation study showed that the multiple-biomarker model was more accurate than single biomarkers (mean absolute error 1.01 vs ≥1.95). We estimate that 816 [95% credible interval (CI) 775-865] PLHIV were undiagnosed in 2015, representing a proportion of 10.8% (95% CI 10.3-11.4%) of all PLHIV. CONCLUSIONS: The proposed methodology will enhance the utility of standard surveillance data streams and will be useful to monitor progress towards and compliance with the 90-90-90 UNAIDS target.

Estimation in emerging epidemics: biases and remedies.

When analysing new emerging infectious disease outbreaks, one typically has observational data over a limited period of time and several parameters to estimate, such as growth rate, the basic reproduction number R0, the case fatality rate and distributions of serial intervals, generation times, latency and incubation times and times between onset of symptoms, notification, death and recovery/discharge. These parameters form the basis for predicting a future outbreak, planning preventive measures and monitoring the progress of the disease outbreak. We study inference problems during the emerging phase of an outbreak, and point out potential sources of bias, with emphasis on: contact tracing backwards in time, replacing generation times by serial intervals, multiple potential infectors and censoring effects amplified by exponential growth. These biases directly affect the estimation of, for example, the generation time distribution and the case fatality rate, but can then propagate to other estimates such as R0 and growth rate. We propose methods to remove or at least reduce bias using statistical modelling. We illustrate the theory by numerical examples and simulations.