Effectiveness of reactive focal mass drug administration and reactive focal vector control to reduce malaria transmission in the low malaria-endemic setting of Namibia: a cluster-randomised controlled, open-label, two-by-two factorial design trial.

BACKGROUND: In low malaria-endemic settings, screening and treatment of individuals in close proximity to index cases, also known as reactive case detection (RACD), is practised for surveillance and response. However, other approaches could be more effective for reducing transmission. We aimed to evaluate the effectiveness of reactive focal mass drug administration (rfMDA) and reactive focal vector control (RAVC) in the low malaria-endemic setting of Zambezi (Namibia). METHODS: We did a cluster-randomised controlled, open-label trial using a two-by-two factorial design of 56 enumeration area clusters in the low malaria-endemic setting of Zambezi (Namibia). We randomly assigned these clusters using restricted randomisation to four groups: RACD only, rfMDA only, RAVC plus RACD, or rfMDA plus RAVC. RACD involved rapid diagnostic testing and treatment with artemether-lumefantrine and single-dose primaquine, rfMDA involved presumptive treatment with artemether-lumefantrine, and RAVC involved indoor residual spraying with pirimiphos-methyl. Interventions were administered within 500 m of index cases. To evaluate the effectiveness of interventions targeting the parasite reservoir in humans (rfMDA vs RACD), in mosquitoes (RAVC vs no RAVC), and in both humans and mosquitoes (rfMDA plus RAVC vs RACD only), an intention-to-treat analysis was done. For each of the three comparisons, the primary outcome was the cumulative incidence of locally acquired malaria cases. This trial is registered with ClinicalTrials.gov, number NCT02610400. FINDINGS: Between Jan 1, 2017, and Dec 31, 2017, 55 enumeration area clusters had 1118 eligible index cases that led to 342 interventions covering 8948 individuals. The cumulative incidence of locally acquired malaria was 30·8 per 1000 person-years (95% CI 12·8-48·7) in the clusters that received rfMDA versus 38·3 per 1000 person-years (23·0-53·6) in the clusters that received RACD; 30·2 per 1000 person-years (15·0-45·5) in the clusters that received RAVC versus 38·9 per 1000 person-years (20·7-57·1) in the clusters that did not receive RAVC; and 25·0 per 1000 person-years (5·2-44·7) in the clusters that received rfMDA plus RAVC versus 41·4 per 1000 person-years (21·5-61·2) in the clusters that received RACD only. After adjusting for imbalances in baseline and implementation factors, the incidence of malaria was lower in clusters receiving rfMDA than in those receiving RACD (adjusted incidence rate ratio 0·52 [95% CI 0·16-0·88], p=0·009), lower in clusters receiving RAVC than in those that did not (0·48 [0·16-0·80], p=0·002), and lower in clusters that received rfMDA plus RAVC than in those receiving RACD only (0·26 [0·10-0·68], p=0·006). No serious adverse events were reported. INTERPRETATION: In a low malaria-endemic setting, rfMDA and RAVC, implemented alone and in combination, reduced malaria transmission and should be considered as alternatives to RACD for elimination of malaria. FUNDING: Novartis Foundation, Bill & Melinda Gates Foundation, and Horchow Family Fund.

Unnatural selection of salmon life histories in a modified riverscape.

Altered river flows and fragmented habitats often simplify riverine communities and favor non-native fishes, but their influence on life-history expression and survival is less clear. Here, we quantified the expression and ultimate success of diverse salmon emigration behaviors in an anthropogenically altered California river system. We analyzed two decades of Chinook salmon monitoring data to explore the influence of regulated flows on juvenile emigration phenology, abundance, and recruitment. We then followed seven cohorts into adulthood using otolith (ear stone) chemical archives to identify patterns in time- and size-selective mortality along the migratory corridor. Suppressed winter flow cues were associated with delayed emigration timing, particularly in warm, dry years, which was also when selection against late migrants was the most extreme. Lower, less variable flows were also associated with reduced juvenile and adult production, highlighting the importance of streamflow for cohort success in these southernmost populations. While most juveniles emigrated from the natal stream as fry or smolts, the survivors were dominated by the rare few that left at intermediate sizes and times, coinciding with managed flows released before extreme summer temperatures. The consistent selection against early (small) and late (large) migrants counters prevailing ecological theory that predicts different traits to be favored under varying environmental conditions. Yet, even with this weakened portfolio, maintaining a broad distribution in migration traits still increased adult production and reduced variance. In years exhibiting large fry pulses, even marginal increases in their survival would have significantly boosted recruitment. However, management actions favoring any single phenotype could have negative evolutionary and demographic consequences, potentially reducing adaptability and population stability. To recover fish populations and support viable fisheries in a warming and increasingly unpredictable climate, coordinating flow and habitat management within and among watersheds will be critical to balance trait optimization versus diversification.

Malaria risk in young male travellers but local transmission persists: a case-control study in low transmission Namibia.

BACKGROUND: A key component of malaria elimination campaigns is the identification and targeting of high risk populations. To characterize high risk populations in north central Namibia, a prospective health facility-based case-control study was conducted from December 2012-July 2014. Cases (n = 107) were all patients presenting to any of the 46 health clinics located in the study districts with a confirmed Plasmodium infection by multi-species rapid diagnostic test (RDT). Population controls (n = 679) for each district were RDT negative individuals residing within a household that was randomly selected from a census listing using a two-stage sampling procedure. Demographic, travel, socio-economic, behavioural, climate and vegetation data were also collected. Spatial patterns of malaria risk were analysed. Multivariate logistic regression was used to identify risk factors for malaria. RESULTS: Malaria risk was observed to cluster along the border with Angola, and travel patterns among cases were comparatively restricted to northern Namibia and Angola. Travel to Angola was associated with excessive risk of malaria in males (OR 43.58 95% CI 2.12-896), but there was no corresponding risk associated with travel by females. This is the first study to reveal that gender can modify the effect of travel on risk of malaria. Amongst non-travellers, male gender was also associated with a higher risk of malaria compared with females (OR 1.95 95% CI 1.25-3.04). Other strong risk factors were sleeping away from the household the previous night, lower socioeconomic status, living in an area with moderate vegetation around their house, experiencing moderate rainfall in the month prior to diagnosis and living <15 km from the Angolan border. CONCLUSIONS: These findings highlight the critical need to target malaria interventions to young male travellers, who have a disproportionate risk of malaria in northern Namibia, to coordinate cross-border regional malaria prevention initiatives and to scale up coverage of prevention measures such as indoor residual spraying and long-lasting insecticide nets in high risk areas if malaria elimination is to be realized.

The changing epidemiology of malaria elimination: new strategies for new challenges.

Malaria-eliminating countries achieved remarkable success in reducing their malaria burdens between 2000 and 2010. As a result, the epidemiology of malaria in these settings has become more complex. Malaria is increasingly imported, caused by Plasmodium vivax in settings outside sub-Saharan Africa, and clustered in small geographical areas or clustered demographically into subpopulations, which are often predominantly adult men, with shared social, behavioural, and geographical risk characteristics. The shift in the populations most at risk of malaria raises important questions for malaria-eliminating countries, since traditional control interventions are likely to be less effective. Approaches to elimination need to be aligned with these changes through the development and adoption of novel strategies and methods. Knowledge of the changing epidemiological trends of malaria in the eliminating countries will ensure improved targeting of interventions to continue to shrink the malaria map.

Fine-scale malaria risk mapping from routine aggregated case data.

BACKGROUND: Mapping malaria risk is an integral component of efficient resource allocation. Routine health facility data are convenient to collect, but without information on the locations at which transmission occurred, their utility for predicting variation in risk at a sub-catchment level is presently unclear. METHODS: Using routinely collected health facility level case data in Swaziland between 2011-2013, and fine scale environmental and ecological variables, this study explores the use of a hierarchical Bayesian modelling framework for downscaling risk maps from health facility catchment level to a fine scale (1 km x 1 km). Fine scale predictions were validated using known household locations of cases and a random sample of points to act as pseudo-controls. RESULTS: Results show that fine-scale predictions were able to discriminate between cases and pseudo-controls with an AUC value of 0.84. When scaled up to catchment level, predicted numbers of cases per health facility showed broad correspondence with observed numbers of cases with little bias, with 84 of the 101 health facilities with zero cases correctly predicted as having zero cases. CONCLUSIONS: This method holds promise for helping countries in pre-elimination and elimination stages use health facility level data to produce accurate risk maps at finer scales. Further validation in other transmission settings and an evaluation of the operational value of the approach is necessary.

The independent effect of living in malaria hotspots on future malaria infection: an observational study from Misungwi, Tanzania.

BACKGROUND: As malaria transmission declines, continued improvements of prevention and control interventions will increasingly rely on accurate knowledge of risk factors and an ability to define high-risk areas and populations at risk for focal targeting of interventions. This paper explores the independent association between living in a hotspot and prospective risk of malaria infection. METHODS: Malaria infection status defined by nPCR and AMA-1 status in year 1 were used to define geographic hotspots using two geospatial statistical methods (SaTScan and Kernel density smoothing). Other malaria risk factors for malaria infection were explored by fitting a multivariable model. RESULTS: This study demonstrated that residing in infection hotspot of malaria transmission is an independent predictor of malaria infection in the future. CONCLUSION: It is likely that targeting such hotspots with better coverage and improved malaria control strategies will result in more cost-efficient uses of resources to move towards malaria elimination.

Hot spot or not: a comparison of spatial statistical methods to predict prospective malaria infections.

BACKGROUND: Within affected communities, Plasmodium falciparum infections may be skewed in distribution such that single or small clusters of households consistently harbour a disproportionate number of infected individuals throughout the year. Identifying these hotspots of malaria transmission would permit targeting of interventions and a more rapid reduction in malaria burden across the whole community. This study set out to compare different statistical methods of hotspot detection (SaTScan, kernel smoothing, weighted local prevalence) using different indicators (PCR positivity, AMA-1 and MSP-1 antibodies) for prediction of infection the following year. METHODS: Two full surveys of four villages in Mwanza, Tanzania were completed over consecutive years, 2010-2011. In both surveys, infection was assessed using nested polymerase chain reaction (nPCR). In addition in 2010, serologic markers (AMA-1 and MSP-119 antibodies) of exposure were assessed. Baseline clustering of infection and serological markers were assessed using three geospatial methods: spatial scan statistics, kernel analysis and weighted local prevalence analysis. Methods were compared in their ability to predict infection in the second year of the study using random effects logistic regression models, and comparisons of the area under the receiver operating curve (AUC) for each model. Sensitivity analysis was conducted to explore the effect of varying radius size for the kernel and weighted local prevalence methods and maximum population size for the spatial scan statistic. RESULTS: Guided by AUC values, the kernel method and spatial scan statistics appeared to be more predictive of infection in the following year. Hotspots of PCR-detected infection and seropositivity to AMA-1 were predictive of subsequent infection. For the kernel method, a 1 km window was optimal. Similarly, allowing hotspots to contain up to 50% of the population was a better predictor of infection in the second year using spatial scan statistics than smaller maximum population sizes. CONCLUSIONS: Clusters of AMA-1 seroprevalence or parasite prevalence that are predictive of infection a year later can be identified using geospatial models. Kernel smoothing using a 1 km window and spatial scan statistics both provided accurate prediction of future infection.

Epidemiology of subpatent Plasmodium falciparum infection: implications for detection of hotspots with imperfect diagnostics.

BACKGROUND: At the local level, malaria transmission clusters in hotspots, which may be a group of households that experience higher than average exposure to infectious mosquitoes. Active case detection often relying on rapid diagnostic tests for mass screen and treat campaigns has been proposed as a method to detect and treat individuals in hotspots. Data from a cross-sectional survey conducted in north-western Tanzania were used to examine the spatial distribution of Plasmodium falciparum and the relationship between household exposure and parasite density. METHODS: Dried blood spots were collected from consenting individuals from four villages during a survey conducted in 2010. These were analysed by PCR for the presence of P. falciparum, with the parasite density of positive samples being estimated by quantitative PCR. Household exposure was estimated using the distance-weighted PCR prevalence of infection. Parasite density simulations were used to estimate the proportion of infections that would be treated using a screen and treat approach with rapid diagnostic tests (RDT) compared to targeted mass drug administration (tMDA) and Mass Drug Administration (MDA). RESULTS: Polymerase chain reaction PCR analysis revealed that of the 3,057 blood samples analysed, 1,078 were positive. Mean distance-weighted PCR prevalence per household was 34.5%. Parasite density was negatively associated with transmission intensity with the odds of an infection being subpatent increasing with household exposure (OR 1.09 per 1% increase in exposure). Parasite density was also related to age, being highest in children five to ten years old and lowest in those > 40 years. Simulations of different tMDA strategies showed that treating all individuals in households where RDT prevalence was above 20% increased the number of infections that would have been treated from 43 to 55%. However, even with this strategy, 45% of infections remained untreated. CONCLUSION: The negative relationship between household exposure and parasite density suggests that DNA-based detection of parasites is needed to provide adequate sensitivity in hotspots. Targeting MDA only to households with RDT-positive individuals may allow a larger fraction of infections to be treated. These results suggest that community-wide MDA, instead of screen and treat strategies, may be needed to successfully treat the asymptomatic, subpatent parasite reservoir and reduce transmission in similar settings.

Spatial parasite ecology and epidemiology: a review of methods and applications.

The distributions of parasitic diseases are determined by complex factors, including many that are distributed in space. A variety of statistical methods are now readily accessible to researchers providing opportunities for describing and ultimately understanding and predicting spatial distributions. This review provides an overview of the spatial statistical methods available to parasitologists, ecologists and epidemiologists and discusses how such methods have yielded new insights into the ecology and epidemiology of infection and disease. The review is structured according to the three major branches of spatial statistics: continuous spatial variation; discrete spatial variation; and spatial point processes.

Predicting future community-level ocular Chlamydia trachomatis infection prevalence using serological, clinical, molecular, and geospatial data.

Trachoma is an infectious disease characterized by repeated exposures to Chlamydia trachomatis (Ct) that may ultimately lead to blindness. Efficient identification of communities with high infection burden could help target more intensive control efforts. We hypothesized that IgG seroprevalence in combination with geospatial layers, machine learning, and model-based geostatistics would be able to accurately predict future community-level ocular Ct infections detected by PCR. We used measurements from 40 communities in the hyperendemic Amhara region of Ethiopia to assess this hypothesis. Median Ct infection prevalence among children 0-5 years old increased from 6% at enrollment, in the context of recent mass drug administration (MDA), to 29% by month 36, following three years without MDA. At baseline, correlation between seroprevalence and Ct infection was stronger among children 0-5 years old (ρ = 0.77) than children 6-9 years old (ρ = 0.48), and stronger than the correlation between active trachoma and Ct infection (0-5y ρ = 0.56; 6-9y ρ = 0.40). Seroprevalence was the strongest concurrent predictor of infection prevalence at month 36 among children 0-5 years old (cross-validated R2 = 0.75, 95% CI: 0.58-0.85), though predictive performance declined substantially with increasing temporal lag between predictor and outcome measurements. Geospatial variables, a spatial Gaussian process, and stacked ensemble machine learning did not meaningfully improve predictions. Serological markers among children 0-5 years old may be an objective tool for identifying communities with high levels of ocular Ct infections, but accurate, future prediction in the context of changing transmission remains an open challenge.

Association of Travel Distance to Nearest Abortion Facility With Rates of Abortion.

Importance: Travel distance to abortion services varies widely in the US. Some evidence shows travel distance affects use of abortion care, but there is no national analysis of how abortion rate changes with travel distance. Objective: To examine the association between travel distance to the nearest abortion care facility and the abortion rate and to model the effect of reduced travel distance. Design, Setting, and Participants: This cross-sectional geographic analysis used 2015 data on abortions by county of residence from 1948 counties in 27 states. Abortion rates were modeled using a spatial Poisson model adjusted for age, race/ethnicity, marital status, educational attainment, household poverty, nativity, and state abortion policies. Abortion rates for 3107 counties in the 48 contiguous states that were home to 62.5 million female residents of reproductive age (15-44 years) and changes under travel distance scenarios, including integration into primary care (<30 miles) and availability of telemedicine care (<5 miles), were estimated. Data were collected from April 2018 to October 2019 and analyzed from December 2019 to July 2020. Exposures: Median travel distance by car to the nearest abortion facility. Main Outcomes and Measures: US county abortion rate per 1000 female residents of reproductive age. Results: Among the 1948 counties included in the analysis, greater travel distances were associated with lower abortion rates in a dose-response manner. Compared with a median travel distance of less than 5 miles (median rate, 21.1 [range, 1.2-63.6] per 1000 female residents of reproductive age), distances of 5 to 15 miles (median rate, 12.2 [range, 0.5-23.4] per 1000 female residents of reproductive age; adjusted coefficient, -0.05 [95% CI, -0.07 to -0.03]) and 120 miles or more (median rate, 3.9 [range, 0-12.9] per 1000 female residents of reproductive age; coefficient, -0.73 [95% CI, -0.80 to -0.65]) were associated with lower rates. In a model of 3107 counties with 62.5 million female residents of reproductive age, 696 760 abortions were estimated (mean rate, 11.1 [range, 1.0-45.5] per 1000 female residents of reproductive age). If abortion were integrated into primary care, an additional 18 190 abortions (mean rate, 11.4 [range, 1.1-45.5] per 1000 female residents of reproductive age) were estimated. If telemedicine were widely available, an additional 70 920 abortions were estimated (mean rate, 12.3 [range, 1.4-45.5] per 1000 female residents of reproductive age). Conclusions and Relevance: These findings suggest that greater travel distances to abortion services are associated with lower abortion rates. The results indicate which geographic areas have insufficient access to abortion care. Modeling suggests that integrating abortion into primary care or making medication abortion care available by telemedicine may decrease unmet need.

Population size estimation of seasonal forest-going populations in southern Lao PDR.

Forest-going populations are key to malaria transmission in the Greater Mekong Sub-region (GMS) and are therefore targeted for elimination efforts. Estimating the size of this population is essential for programs to assess, track and achieve their elimination goals. Leveraging data from three cross-sectional household surveys and one survey among forest-goers, the size of this high-risk population in a southern province of Lao PDR between December 2017 and November 2018 was estimated by two methods: population-based household surveys and capture-recapture. During the first month of the dry season, the first month of the rainy season, and the last month of the rainy season, respectively, 16.2% [14.7; 17.7], 9.3% [7.2; 11.3], and 5.3% [4.4; 6.1] of the adult population were estimated to have engaged in forest-going activities. The capture-recapture method estimated a total population size of 18,426 [16,529; 20,669] forest-goers, meaning 61.0% [54.2; 67.9] of the adult population had engaged in forest-going activities over the 12-month study period. This study demonstrates two methods for population size estimation to inform malaria research and programming. The seasonality and turnover within this forest-going population provide unique opportunities and challenges for control programs across the GMS as they work towards malaria elimination.

Finding hotspots: development of an adaptive spatial sampling approach.

The identification of disease hotspots is an increasingly important public health problem. While geospatial modeling offers an opportunity to predict the locations of hotspots using suitable environmental and climatological data, little attention has been paid to optimizing the design of surveys used to inform such models. Here we introduce an adaptive sampling scheme optimized to identify hotspot locations where prevalence exceeds a relevant threshold. Our approach incorporates ideas from Bayesian optimization theory to adaptively select sample batches. We present an experimental simulation study based on survey data of schistosomiasis and lymphatic filariasis across four countries. Results across all scenarios explored show that adaptive sampling produces superior results and suggest that similar performance to random sampling can be achieved with a fraction of the sample size.

Household electricity access in Africa (2000-2013): Closing information gaps with model-based geostatistics.

Household electricity access data in Africa are scarce, particularly at the subnational level. We followed a model-based Geostatistics approach to produce maps of electricity access between 2000 and 2013 at a 5 km resolution. We collated data from 69 nationally representative household surveys conducted in Africa and incorporated nighttime lights imagery as well as land use and land cover data to produce maps of electricity access between 2000 and 2013. The information produced here can be an aid for understanding of how electricity access has changed in the region during this 14 year period. The resolution and the continental scale makes it possible to combine these data with other sources in applications in the socio-economic field, both at a local or regional level.

Predicting residential structures from open source remotely enumerated data using machine learning.

Having accurate maps depicting the locations of residential buildings across a region benefits a range of sectors. This is particularly true for public health programs focused on delivering services at the household level, such as indoor residual spraying with insecticide to help prevent malaria. While open source data from OpenStreetMap (OSM) depicting the locations and shapes of buildings is rapidly improving in terms of quality and completeness globally, even in settings where all buildings have been mapped, information on whether these buildings are residential, commercial or another type is often only available for a small subset. Using OSM building data from Botswana and Swaziland, we identified buildings for which 'type' was indicated, generated via on the ground observations, and classified these into two classes, "sprayable" and "not-sprayable". Ensemble machine learning, using building characteristics such as size, shape and proximity to neighbouring features, was then used to form a model to predict which of these 2 classes every building in these two countries fell into. Results show that an ensemble machine learning approach performed marginally, but statistically, better than the best individual model and that using this ensemble model we were able to correctly classify >86% (using independent test data) of structures correctly as sprayable and not-sprayable across both countries.

Spatial clustering of patent and sub-patent malaria infections in northern Namibia: Implications for surveillance and response strategies for elimination.

BACKGROUND: Reactive case detection (RACD) around passively detected malaria cases is a strategy to identify and treat hotspots of malaria transmission. This study investigated the unproven assumption on which this approach is based, that in low transmission settings, infections cluster over small scales. METHODS: A prospective case-control study was conducted between January 2013 and August 2014 in Ohangwena and Omusati regions in north central Namibia. Patients attending health facilities who tested positive by malaria rapid diagnostic test (RDT) (index cases) were traced back to their home. All occupants of index case households (n = 116 households) and surrounding households (n = 225) were screened for Plasmodium infection with a rapid diagnostic test (RDT) and loop mediated isothermal amplification (LAMP) and interviewed to identify risk factors. A comparison group of 286 randomly-selected control households was also screened, to compare infection levels of RACD and non-RACD households and their neighbours. Logistic regression was used to investigate spatial clustering of patent and sub-patent infections around index cases and to identify potential risk factors that would inform screening approaches and identify risk groups. Estimates of the impact of RACD on onward transmission to mosquitoes was made using previously published figures of infection rates. RESULTS: Prevalence of Plasmodium falciparum infection by LAMP was 3.4%, 1.4% and 0.4% in index-case households, neighbors of index case households and control households respectively; adjusted odds ratio 6.1 [95%CI 1.9-19.5] comparing case households versus control households. Using data from Engela, neighbors of cases had higher odds of infection [adjusted OR 5.0 95%CI 1.3-18.9] compared to control households. All infections identified by RDTs were afebrile and RDTs identified only a small proportion of infections in case (n = 7; 17%) and control (0%) neighborhoods. Based on published estimates of patent and sub-patent infectiousness, these results suggest that infections missed by RDTs during RACD would allow 50-71% of infections to mosquitoes to occur in this setting. CONCLUSION: Malaria infections cluster around passively detected cases. The majority of infections are asymptomatic and of densities below the limit of detection of current RDTs. RACD using standard RDTs are unlikely to detect enough malaria infections to dramatically reduce transmission. In low transmission settings such as Namibia more sensitive field diagnostics or forms of focal presumptive treatment should be tested as strategies to reduce malaria transmission.

Mapping land cover change over continental Africa using Landsat and Google Earth Engine cloud computing.

Quantifying and monitoring the spatial and temporal dynamics of the global land cover is critical for better understanding many of the Earth's land surface processes. However, the lack of regularly updated, continental-scale, and high spatial resolution (30 m) land cover data limit our ability to better understand the spatial extent and the temporal dynamics of land surface changes. Despite the free availability of high spatial resolution Landsat satellite data, continental-scale land cover mapping using high resolution Landsat satellite data was not feasible until now due to the need for high-performance computing to store, process, and analyze this large volume of high resolution satellite data. In this study, we present an approach to quantify continental land cover and impervious surface changes over a long period of time (15 years) using high resolution Landsat satellite observations and Google Earth Engine cloud computing platform. The approach applied here to overcome the computational challenges of handling big earth observation data by using cloud computing can help scientists and practitioners who lack high-performance computational resources.