Geospatial joint modeling of vector and parasite serology to microstratify malaria transmission.

The World Health Organization identifies a strong surveillance system for malaria and its mosquito vector as an essential pillar of the malaria elimination agenda. Anopheles salivary antibodies are emerging biomarkers of exposure to mosquito bites that potentially overcome sensitivity and logistical constraints of traditional entomological surveys. Using samples collected by a village health volunteer network in 104 villages in Southeast Myanmar during routine surveillance, the present study employs a Bayesian geostatistical modeling framework, incorporating climatic and environmental variables together with Anopheles salivary antigen serology, to generate spatially continuous predictive maps of Anopheles biting exposure. Our maps quantify fine-scale spatial and temporal heterogeneity in Anopheles salivary antibody seroprevalence (ranging from 9 to 99%) that serves as a proxy of exposure to Anopheles bites and advances current static maps of only Anopheles occurrence. We also developed an innovative framework to perform surveillance of malaria transmission. By incorporating antibodies against the vector and the transmissible form of malaria (sporozoite) in a joint Bayesian geostatistical model, we predict several foci of ongoing transmission. In our study, we demonstrate that antibodies specific for Anopheles salivary and sporozoite antigens are a logistically feasible metric with which to quantify and characterize heterogeneity in exposure to vector bites and malaria transmission. These approaches could readily be scaled up into existing village health volunteer surveillance networks to identify foci of residual malaria transmission, which could be targeted with supplementary interventions to accelerate progress toward elimination.

Fine-scale maps of malaria incidence to inform risk stratification in Laos.

BACKGROUND: Malaria risk maps are crucial for controlling and eliminating malaria by identifying areas of varying transmission risk. In the Greater Mekong Subregion, these maps guide interventions and resource allocation. This article focuses on analysing changes in malaria transmission and developing fine-scale risk maps using five years of routine surveillance data in Laos (2017-2021). The study employed data from 1160 geolocated health facilities in Laos, along with high-resolution environmental data. METHODS: A Bayesian geostatistical framework incorporating population data and treatment-seeking propensity was developed. The models incorporated static and dynamic factors and accounted for spatial heterogeneity. RESULTS: Results showed a significant decline in malaria cases in Laos over the five-year period and a shift in transmission patterns. While the north became malaria-free, the south experienced ongoing transmission with sporadic outbreaks. CONCLUSION: The risk maps provided insights into changing transmission patterns and supported risk stratification. These risk maps are valuable tools for malaria control in Laos, aiding resource allocation, identifying intervention gaps, and raising public awareness. The study enhances understanding of malaria transmission dynamics and facilitates evidence-based decision-making for targeted interventions in high-risk areas.

Mapping malaria seasonality in Madagascar using health facility data.

BACKGROUND: Many malaria-endemic areas experience seasonal fluctuations in case incidence as Anopheles mosquito and Plasmodium parasite life cycles respond to changing environmental conditions. Identifying location-specific seasonality characteristics is useful for planning interventions. While most existing maps of malaria seasonality use fixed thresholds of rainfall, temperature, and/or vegetation indices to identify suitable transmission months, we construct a statistical modelling framework for characterising the seasonal patterns derived directly from monthly health facility data. METHODS: With data from 2669 of the 3247 health facilities in Madagascar, a spatiotemporal regression model was used to estimate seasonal patterns across the island. In the absence of catchment population estimates or the ability to aggregate to the district level, this focused on the monthly proportions of total annual cases by health facility level. The model was informed by dynamic environmental covariates known to directly influence seasonal malaria trends. To identify operationally relevant characteristics such as the transmission start months and associated uncertainty measures, an algorithm was developed and applied to model realisations. A seasonality index was used to incorporate burden information from household prevalence surveys and summarise 'how seasonal' locations are relative to their surroundings. RESULTS: Positive associations were detected between monthly case proportions and temporally lagged covariates of rainfall and temperature suitability. Consistent with the existing literature, model estimates indicate that while most parts of Madagascar experience peaks in malaria transmission near March-April, the eastern coast experiences an earlier peak around February. Transmission was estimated to start in southeast districts before southwest districts, suggesting that indoor residual spraying should be completed in the same order. In regions where the data suggested conflicting seasonal signals or two transmission seasons, estimates of seasonal features had larger deviations and therefore less certainty. CONCLUSIONS: Monthly health facility data can be used to establish seasonal patterns in malaria burden and augment the information provided by household prevalence surveys. The proposed modelling framework allows for evidence-based and cohesive inferences on location-specific seasonal characteristics. As health surveillance systems continue to improve, it is hoped that more of such data will be available to improve our understanding and planning of intervention strategies.

Spatio-temporal mapping of Madagascar's Malaria Indicator Survey results to assess Plasmodium falciparum endemicity trends between 2011 and 2016.

BACKGROUND: Reliable measures of disease burden over time are necessary to evaluate the impact of interventions and assess sub-national trends in the distribution of infection. Three Malaria Indicator Surveys (MISs) have been conducted in Madagascar since 2011. They provide a valuable resource to assess changes in burden that is complementary to the country's routine case reporting system. METHODS: A Bayesian geostatistical spatio-temporal model was developed in an integrated nested Laplace approximation framework to map the prevalence of Plasmodium falciparum malaria infection among children from 6 to 59 months in age across Madagascar for 2011, 2013 and 2016 based on the MIS datasets. The model was informed by a suite of environmental and socio-demographic covariates known to influence infection prevalence. Spatio-temporal trends were quantified across the country. RESULTS: Despite a relatively small decrease between 2013 and 2016, the prevalence of malaria infection has increased substantially in all areas of Madagascar since 2011. In 2011, almost half (42.3%) of the country's population lived in areas of very low malaria risk (<1% parasite prevalence), but by 2016, this had dropped to only 26.7% of the population. Meanwhile, the population in high transmission areas (prevalence >20%) increased from only 2.2% in 2011 to 9.2% in 2016. A comparison of the model-based estimates with the raw MIS results indicates there was an underestimation of the situation in 2016, since the raw figures likely associated with survey timings were delayed until after the peak transmission season. CONCLUSIONS: Malaria remains an important health problem in Madagascar. The monthly and annual prevalence maps developed here provide a way to evaluate the magnitude of change over time, taking into account variability in survey input data. These methods can contribute to monitoring sub-national trends of malaria prevalence in Madagascar as the country aims for geographically progressive elimination.

Bayesian adaptive design: improving the effectiveness of monitoring of the Great Barrier Reef.

Monitoring programs are essential for understanding patterns, trends, and threats in ecological and environmental systems. However, such programs are costly in terms of dollars, human resources, and technology, and complex in terms of balancing short- and long-term requirements. In this work, We develop new statistical methods for implementing cost-effective adaptive sampling and monitoring schemes for coral reef that can better utilize existing information and resources, and which can incorporate available prior information. Our research was motivated by developing efficient monitoring practices for Australia's Great Barrier Reef. We develop and implement two types of adaptive sampling schemes, static and sequential, and show that they can be more informative and cost-effective than an existing (nonadaptive) monitoring program. Our methods are developed in a Bayesian framework with a range of utility functions relevant to environmental monitoring. Our results demonstrate the considerable potential for adaptive design to support improved management outcomes in comparison to set-and-forget styles of surveillance monitoring.

System-level metabolic modeling facilitates unveiling metabolic signature in exceptional longevity.

Although it is well known that metabolic control plays a crucial role in regulating the health span and life span of various organisms, little is known for the systems metabolic profile of centenarians, the paradigm of human healthy aging and longevity. Meanwhile, how to well characterize the system-level metabolic states in an organism of interest remains to be a major challenge in systems metabolism research. To address this challenge and better understand the metabolic mechanisms of healthy aging, we developed a method of genome-wide precision metabolic modeling (GPMM) which is able to quantitatively integrate transcriptome, proteome and kinetome data in predictive modeling of metabolic networks. Benchmarking analysis showed that GPMM successfully characterized metabolic reprogramming in the NCI-60 cancer cell lines; it dramatically improved the performance of the modeling with an R2 of 0.86 between the predicted and experimental measurements over the performance of existing methods. Using this approach, we examined the metabolic networks of a Chinese centenarian cohort and identified the elevated fatty acid oxidation (FAO) as the most significant metabolic feature in these long-lived individuals. Evidence from serum metabolomics supports this observation. Given that FAO declines with normal aging and is impaired in many age-related diseases, our study suggests that the elevated FAO has potential to be a novel signature of healthy aging of humans.

Human mobility patterns and malaria importation on Bioko Island.

Malaria burden on Bioko Island has decreased significantly over the past 15 years. The impact of interventions on malaria prevalence, however, has recently stalled. Here, we use data from island-wide, annual malaria indicator surveys to investigate human movement patterns and their relationship to Plasmodium falciparum prevalence. Using geostatistical and mathematical modelling, we find that off-island travel is more prevalent in and around the capital, Malabo. The odds of malaria infection among off-island travelers are significantly higher than the rest of the population. We estimate that malaria importation rates are high enough to explain malaria prevalence in much of Malabo and its surroundings, and that local transmission is highest along the West Coast of the island. Despite uncertainty, these estimates of residual transmission and importation serve as a basis for evaluating progress towards elimination and for efficiently allocating resources as Bioko makes the transition from control to elimination.

Pareto rules for malaria super-spreaders and super-spreading.

Heterogeneity in transmission is a challenge for infectious disease dynamics and control. An 80-20 "Pareto" rule has been proposed to describe this heterogeneity whereby 80% of transmission is accounted for by 20% of individuals, herein called super-spreaders. It is unclear, however, whether super-spreading can be attributed to certain individuals or whether it is an unpredictable and unavoidable feature of epidemics. Here, we investigate heterogeneous malaria transmission at three sites in Uganda and find that super-spreading is negatively correlated with overall malaria transmission intensity. Mosquito biting among humans is 90-10 at the lowest transmission intensities declining to less than 70-30 at the highest intensities. For super-spreaders, biting ranges from 70-30 down to 60-40. The difference, approximately half the total variance, is due to environmental stochasticity. Super-spreading is thus partly due to super-spreaders, but modest gains are expected from targeting super-spreaders.

Heterogeneous exposure and hotspots for malaria vectors at three study sites in Uganda.

Background: Heterogeneity in malaria transmission has household, temporal, and spatial components. These factors are relevant for improving the efficiency of malaria control by targeting heterogeneity. To quantify variation, we analyzed mosquito counts from entomological surveillance conducted at three study sites in Uganda that varied in malaria transmission intensity. Mosquito biting or exposure is a risk factor for malaria transmission. Methods: Using a Bayesian zero-inflated negative binomial model, validated via a comprehensive simulation study, we quantified household differences in malaria vector density and examined its spatial distribution. We introduced a novel approach for identifying changes in vector abundance hotspots over time by computing the Getis-Ord statistic on ratios of household biting propensities for different scenarios. We also explored the association of household biting propensities with housing and environmental covariates. Results: In each site, there was evidence for hot and cold spots of vector abundance, and spatial patterns associated with urbanicity, elevation, or other environmental covariates. We found some differences in the hotspots in rainy vs. dry seasons or before vs. after the application of control interventions. Housing quality explained a portion of the variation among households in mosquito counts. Conclusion: This work provided an improved understanding of heterogeneity in malaria vector density at the three study sites in Uganda and offered a valuable opportunity for assessing whether interventions could be spatially targeted to be aimed at abundance hotspots which may increase malaria risk. Indoor residual spraying was shown to be a successful measure of vector control interventions in Tororo, Uganda.  Cement walls, brick floors, closed eaves, screened airbricks, and tiled roofs were features of a house that had shown reduction of household biting propensity. Improvements in house quality should be recommended as a supplementary measure for malaria control reducing risk of infection.

Making the most of spatial information in health: a tutorial in Bayesian disease mapping for areal data.

Disease maps are effective tools for explaining and predicting patterns of disease outcomes across geographical space, identifying areas of potentially elevated risk, and formulating and validating aetiological hypotheses for a disease. Bayesian models have become a standard approach to disease mapping in recent decades. This article aims to provide a basic understanding of the key concepts involved in Bayesian disease mapping methods for areal data. It is anticipated that this will help in interpretation of published maps, and provide a useful starting point for anyone interested in running disease mapping methods for areal data. The article provides detailed motivation and descriptions on disease mapping methods by explaining the concepts, defining the technical terms, and illustrating the utility of disease mapping for epidemiological research by demonstrating various ways of visualising model outputs using a case study. The target audience includes spatial scientists in health and other fields, policy or decision makers, health geographers, spatial analysts, public health professionals, and epidemiologists.

The choice of spatial scales and spatial smoothness priors for various spatial patterns.

Given the drawbacks for using geo-political areas in mapping outcomes unrelated to geo-politics, a compromise is to aggregate and analyse data at the grid level. This has the advantage of allowing spatial smoothing and modelling at a biologically or physically relevant scale. This article addresses two consequent issues: the choice of the spatial smoothness prior and the scale of the grid. Firstly, we describe several spatial smoothness priors applicable for grid data and discuss the contexts in which these priors can be employed based on different aims. Two such aims are considered, i.e., to identify regions with clustering and to model spatial dependence in the data. Secondly, the choice of the grid size is shown to depend largely on the spatial patterns. We present a guide on the selection of spatial scales and smoothness priors for various point patterns based on the two aims for spatial smoothing.

The impact of spatial scales and spatial smoothing on the outcome of bayesian spatial model.

Discretization of a geographical region is quite common in spatial analysis. There have been few studies into the impact of different geographical scales on the outcome of spatial models for different spatial patterns. This study aims to investigate the impact of spatial scales and spatial smoothing on the outcomes of modelling spatial point-based data. Given a spatial point-based dataset (such as occurrence of a disease), we study the geographical variation of residual disease risk using regular grid cells. The individual disease risk is modelled using a logistic model with the inclusion of spatially unstructured and/or spatially structured random effects. Three spatial smoothness priors for the spatially structured component are employed in modelling, namely an intrinsic Gaussian Markov random field, a second-order random walk on a lattice, and a Gaussian field with Matérn correlation function. We investigate how changes in grid cell size affect model outcomes under different spatial structures and different smoothness priors for the spatial component. A realistic example (the Humberside data) is analyzed and a simulation study is described. Bayesian computation is carried out using an integrated nested Laplace approximation. The results suggest that the performance and predictive capacity of the spatial models improve as the grid cell size decreases for certain spatial structures. It also appears that different spatial smoothness priors should be applied for different patterns of point data.