A behaviour and disease transmission model: incorporating the Health Belief Model for human behaviour into a simple transmission model.

The health and economic impacts of infectious diseases such as COVID-19 affect all levels of a community from the individual to the governing bodies. However, the spread of an infectious disease is intricately linked to the behaviour of the people within a community since crowd behaviour affects individual human behaviour, while human behaviour affects infection spread, and infection spread affects human behaviour. Capturing these feedback loops of behaviour and infection is a well-known challenge in infectious disease modelling. Here, we investigate the interface of behavioural science theory and infectious disease modelling to explore behaviour and disease (BaD) transmission models. Specifically, we incorporate a visible protective behaviour into the susceptible-infectious-recovered-susceptible (SIRS) transmission model using the socio-psychological Health Belief Model to motivate behavioural uptake and abandonment. We characterize the mathematical thresholds for BaD emergence in the BaD SIRS model and the feasible steady states. We also explore, under different infectious disease scenarios, the effects of a fully protective behaviour on long-term disease prevalence in a community, and describe how BaD modelling can investigate non-pharmaceutical interventions that target-specific components of the Health Belief Model. This transdisciplinary BaD modelling approach may reduce the health and economic impacts of future epidemics.

Mathematical models of Plasmodium vivax transmission: A scoping review.

Plasmodium vivax is one of the most geographically widespread malaria parasites in the world, primarily found across South-East Asia, Latin America, and parts of Africa. One of the significant characteristics of the P. vivax parasite is its ability to remain dormant in the human liver as hypnozoites and subsequently reactivate after the initial infection (i.e. relapse infections). Mathematical modelling approaches have been widely applied to understand P. vivax dynamics and predict the impact of intervention outcomes. Models that capture P. vivax dynamics differ from those that capture P. falciparum dynamics, as they must account for relapses caused by the activation of hypnozoites. In this article, we provide a scoping review of mathematical models that capture P. vivax transmission dynamics published between January 1988 and May 2023. The primary objective of this work is to provide a comprehensive summary of the mathematical models and techniques used to model P. vivax dynamics. In doing so, we aim to assist researchers working on mathematical epidemiology, disease transmission, and other aspects of P. vivax malaria by highlighting best practices in currently published models and highlighting where further model development is required. We categorise P. vivax models according to whether a deterministic or agent-based approach was used. We provide an overview of the different strategies used to incorporate the parasite's biology, use of multiple scales (within-host and population-level), superinfection, immunity, and treatment interventions. In most of the published literature, the rationale for different modelling approaches was driven by the research question at hand. Some models focus on the parasites' complicated biology, while others incorporate simplified assumptions to avoid model complexity. Overall, the existing literature on mathematical models for P. vivax encompasses various aspects of the parasite's dynamics. We recommend that future research should focus on refining how key aspects of P. vivax dynamics are modelled, including spatial heterogeneity in exposure risk and heterogeneity in susceptibility to infection, the accumulation of hypnozoite variation, the interaction between P. falciparum and P. vivax, acquisition of immunity, and recovery under superinfection.

Reproducibility of COVID-era infectious disease models.

Infectious disease modelling has been prominent throughout the COVID-19 pandemic, helping to understand the virus' transmission dynamics and inform response policies. Given their potential importance and translational impact, we evaluated the computational reproducibility of infectious disease modelling articles from the COVID era. We found that four out of 100 randomly sampled studies released between January 2020 and August 2022 could be completely computationally reproduced using the resources provided (e.g., code, data, instructions) whilst a further eight were partially reproducible. For the 100 most highly cited articles from the same period we found that 11 were completely reproducible with a further 22 partially reproducible. Reflecting on our experience, we discuss common issues affecting computational reproducibility and how these might be addressed.

Optimal Interruption of P. vivax Malaria Transmission Using Mass Drug Administration.

Plasmodium vivax is the most geographically widespread malaria-causing parasite resulting in significant associated global morbidity and mortality. One of the factors driving this widespread phenomenon is the ability of the parasites to remain dormant in the liver. Known as 'hypnozoites', they reside in the liver following an initial exposure, before activating later to cause further infections, referred to as 'relapses'. As around 79-96% of infections are attributed to relapses from activating hypnozoites, we expect it will be highly impactful to apply treatment to target the hypnozoite reservoir (i.e. the collection of dormant parasites) to eliminate P. vivax. Treatment with radical cure, for example tafenoquine or primaquine, to target the hypnozoite reservoir is a potential tool to control and/or eliminate P. vivax. We have developed a deterministic multiscale mathematical model as a system of integro-differential equations that captures the complex dynamics of P. vivax hypnozoites and the effect of hypnozoite relapse on disease transmission. Here, we use our multiscale model to study the anticipated effect of radical cure treatment administered via a mass drug administration (MDA) program. We implement multiple rounds of MDA with a fixed interval between rounds, starting from different steady-state disease prevalences. We then construct an optimisation model with three different objective functions motivated on a public health basis to obtain the optimal MDA interval. We also incorporate mosquito seasonality in our model to study its effect on the optimal treatment regime. We find that the effect of MDA interventions is temporary and depends on the pre-intervention disease prevalence (and choice of model parameters) as well as the number of MDA rounds under consideration. The optimal interval between MDA rounds also depends on the objective (combinations of expected intervention outcomes). We find radical cure alone may not be enough to lead to P. vivax elimination under our mathematical model (and choice of model parameters) since the prevalence of infection eventually returns to pre-MDA levels.

Japanese Encephalitis Enzootic and Epidemic Risks across Australia.

Japanese encephalitis virus (JEV) is an arboviral, encephalitogenic, zoonotic flavivirus characterized by its complex epidemiology whose transmission cycle involves reservoir and amplifying hosts, competent vector species and optimal environmental conditions. Although typically endemic in Asia and parts of the Pacific Islands, unprecedented outbreaks in both humans and domestic pigs in southeastern Australia emphasize the virus' expanding geographical range. To estimate areas at highest risk of JEV transmission in Australia, ecological niche models of vectors and waterbirds, a sample of piggery coordinates and feral pig population density models were combined using mathematical and geospatial mapping techniques. These results highlight that both coastal and inland regions across the continent are estimated to have varying risks of enzootic and/or epidemic JEV transmission. We recommend increased surveillance of waterbirds, feral pigs and mosquito populations in areas where domestic pigs and human populations are present.

Estimating the Distribution of Japanese Encephalitis Vectors in Australia Using Ecological Niche Modelling.

Recent Japanese encephalitis virus (JEV) outbreaks in southeastern Australia have sparked interest into epidemiological factors surrounding the virus' novel emergence in this region. Here, the geographic distribution of mosquito species known to be competent JEV vectors in the country was estimated by combining known mosquito occurrences and ecological drivers of distribution to reveal insights into communities at highest risk of infectious disease transmission. Species distribution models predicted that Culex annulirostris and Culex sitiens presence was mostly likely along Australia's eastern and northern coastline, while Culex quinquefasciatus presence was estimated to be most likely near inland regions of southern Australia as well as coastal regions of Western Australia. While Culex annulirostris is considered the dominant JEV vector in Australia, our ecological niche models emphasise the need for further entomological surveillance and JEV research within Australia.

A Multiscale Mathematical Model of Plasmodium Vivax Transmission.

Malaria is caused by Plasmodium parasites which are transmitted to humans by the bite of an infected Anopheles mosquito. Plasmodium vivax is distinct from other malaria species in its ability to remain dormant in the liver (as hypnozoites) and activate later to cause further infections (referred to as relapses). Mathematical models to describe the transmission dynamics of P. vivax have been developed, but most of them fail to capture realistic dynamics of hypnozoites. Models that do capture the complexity tend to involve many governing equations, making them difficult to extend to incorporate other important factors for P. vivax, such as treatment status, age and pregnancy. In this paper, we have developed a multiscale model (a system of integro-differential equations) that involves a minimal set of equations at the population scale, with an embedded within-host model that can capture the dynamics of the hypnozoite reservoir. In this way, we can gain key insights into dynamics of P. vivax transmission with a minimum number of equations at the population scale, making this framework readily scalable to incorporate more complexity. We performed a sensitivity analysis of our multiscale model over key parameters and found that prevalence of P. vivax blood-stage infection increases with both bite rate and number of mosquitoes but decreases with hypnozoite death rate. Since our mathematical model captures the complex dynamics of P. vivax and the hypnozoite reservoir, it has the potential to become a key tool to inform elimination strategies for P. vivax.

Risk prediction system for dengue transmission based on high resolution weather data.

BACKGROUND: Dengue is the fastest spreading vector-borne viral disease, resulting in an estimated 390 million infections annually. Precise prediction of many attributes related to dengue is still a challenge due to the complex dynamics of the disease. Important attributes to predict include: the risk of and risk factors for an infection; infection severity; and the timing and magnitude of outbreaks. In this work, we build a model for predicting the risk of dengue transmission using high-resolution weather data. The level of dengue transmission risk depends on the vector density, hence we predict risk via vector prediction. METHODS AND FINDINGS: We make use of surveillance data on Aedes aegypti larvae collected by the Taiwan Centers for Disease Control as part of the national routine entomological surveillance of dengue, and weather data simulated using the IBM's Containerized Forecasting Workflow, a high spatial- and temporal-resolution forecasting system. We propose a two stage risk prediction system for assessing dengue transmission via Aedes aegypti mosquitoes. In stage one, we perform a logistic regression to determine whether larvae are present or absent at the locations of interest using weather attributes as the explanatory variables. The results are then aggregated to an administrative division, with presence in the division determined by a threshold percentage of larvae positive locations resulting from a bootstrap approach. In stage two, larvae counts are estimated for the predicted larvae positive divisions from stage one, using a zero-inflated negative binomial model. This model identifies the larvae positive locations with 71% accuracy and predicts the larvae numbers producing a coverage probability of 98% over 95% nominal prediction intervals. This two-stage model improves the overall accuracy of identifying larvae positive locations by 29%, and the mean squared error of predicted larvae numbers by 9.6%, against a single-stage approach which uses a zero-inflated binomial regression approach. CONCLUSIONS: We demonstrate a risk prediction system using high resolution weather data can provide valuable insight to the distribution of risk over a geographical region. The work also shows that a two-stage approach is beneficial in predicting risk in non-homogeneous regions, where the risk is localised.

The distribution of the time taken for an epidemic to spread between two communities.

Assessing the risk of disease spread between communities is important in our highly connected modern world. However, the impact of disease- and population-specific factors on the time taken for an epidemic to spread between communities, as well as the impact of stochastic disease dynamics on this spreading time, are not well understood. In this study, we model the spread of an acute infection between two communities ('patches') using a susceptible-infectious-removed (SIR) metapopulation model. We develop approximations to efficiently evaluate the probability of a major outbreak in a second patch given disease introduction in a source patch, and the distribution of the time taken for this to occur. We use these approximations to assess how interventions, which either control disease spread within a patch or decrease the travel rate between patches, change the spreading probability and median spreading time. We find that decreasing the basic reproduction number in the source patch is the most effective way of both decreasing the spreading probability, and delaying epidemic spread to the second patch should this occur. Moreover, we show that the qualitative effects of interventions are the same regardless of the approximations used to evaluate the spreading time distribution, but for some regions in parameter space, quantitative findings depend upon the approximations used. Importantly, if we neglect the possibility that an intervention prevents a large outbreak in the source patch altogether, then intervention effectiveness is not estimated accurately.

Model-Informed Risk Assessment and Decision Making for an Emerging Infectious Disease in the Asia-Pacific Region.

BACKGROUND: Effective response to emerging infectious disease (EID) threats relies on health care systems that can detect and contain localised outbreaks before they reach a national or international scale. The Asia-Pacific region contains low and middle income countries in which the risk of EID outbreaks is elevated and whose health care systems may require international support to effectively detect and respond to such events. The absence of comprehensive data on populations, health care systems and disease characteristics in this region makes risk assessment and decisions about the provision of such support challenging. METHODOLOGY/PRINCIPAL FINDINGS: We describe a mathematical modelling framework that can inform this process by integrating available data sources, systematically explore the effects of uncertainty, and provide estimates of outbreak risk under a range of intervention scenarios. We illustrate the use of this framework in the context of a potential importation of Ebola Virus Disease into the Asia-Pacific region. Results suggest that, across a wide range of plausible scenarios, preemptive interventions supporting the timely detection of early cases provide substantially greater reductions in the probability of large outbreaks than interventions that support health care system capacity after an outbreak has commenced. CONCLUSIONS/SIGNIFICANCE: Our study demonstrates how, in the presence of substantial uncertainty about health care system infrastructure and other relevant aspects of disease control, mathematical models can be used to assess the constraints that limited resources place upon the ability of local health care systems to detect and respond to EID outbreaks in a timely and effective fashion. Our framework can help evaluate the relative impact of these constraints to identify resourcing priorities for health care system support, in order to inform principled and quantifiable decision making.

Strengthening tuberculosis control overseas: who benefits?

BACKGROUND: Although tuberculosis is a major cause of morbidity and mortality worldwide, available funding falls far short of that required for effective control. Economic and spillover consequences of investments in the treatment of tuberculosis are unclear, particularly when steep gradients in the disease and response are linked by population movements, such as that between Papua New Guinea (PNG) and the Australian cross-border region. OBJECTIVE: To undertake an economic evaluation of Australian support for the expansion of basic Directly Observed Treatment, Short Course in the PNG border area of the South Fly from the current level of 14% coverage. METHODS: Both cost-utility analysis and cost-benefit analysis were applied to models that allow for population movement across regions with different characteristics of tuberculosis burden, transmission, and access to treatment. Cost-benefit data were drawn primarily from estimates published by the World Health Organization, and disease transmission data were drawn from a previously published model. RESULTS: Investing $16 million to increase basic Directly Observed Treatment, Short Course coverage in the South Fly generates a net present value of roughly $74 million for Australia (discounted 2005 dollars). The cost per disability-adjusted life-year averted and quality-adjusted life-year saved for PNG is $7 and $4.6, respectively. CONCLUSIONS: Where regions with major disparities in tuberculosis burden and health system resourcing are connected through population movements, investments in tuberculosis control are of mutual benefit, resulting in net health and economic gains on both sides of the border. These findings are likely to inform the case for appropriate investment in tuberculosis control globally.

A metapopulation model of tuberculosis transmission with a case study from high to low burden areas.

Tuberculosis (TB) is a growing problem worldwide, especially with the emergence and high prevalence of multidrug-resistant strains. We develop a metapopulation model for TB spread, which is particularly suited to investigating transmission between areas of high and low prevalence. A case study of cross-border transmission in the Torres Strait region of Australia and Papua New Guinea (PNG) is considered and a sensitivity analysis is conducted. We find that only 6 of the 50 parameters analysed are important to the cumulative number of clinically active TB patients in the entire region. Of these, only the detection rate in PNG is found to be an important intervention parameter. We therefore give insight into the extent the area with the high burden of TB (PNG in the case study) is dominating the TB dynamics of the entire region. Furthermore, the sensitivity analysis results give insight into the data that most important to collect and refine, which is found to be data relating to the PNG parameters.