METHODS & RESULTS: We extend and apply an environmental niche modelling framework as implemented by a previous mapping study of P. knowlesi transmission risk which included data up to 2015. We reviewed the literature from October 2015 through to March 2020 and identified 264 new records of P. knowlesi, with a total of 524 occurrences included in the current study following consolidation with the 2015 study. The modelling framework used in the 2015 study was extended, with changes including the addition of new covariates to capture the effect of deforestation and urbanisation on P. knowlesi transmission.
DISCUSSION: Our map of P. knowlesi relative transmission suitability estimates that the risk posed by the pathogen is highest in Malaysia and Indonesia, with localised areas of high risk also predicted in the Greater Mekong Subregion, The Philippines and Northeast India. These results highlight areas of priority for P. knowlesi surveillance and prospective sampling to address the challenge the disease poses to malaria elimination planning.
METHODS & RESULTS: We extend and apply an environmental niche modelling framework as implemented by a previous mapping study of P. knowlesi transmission risk which included data up to 2015. We reviewed the literature from October 2015 through to March 2020 and identified 264 new records of P. knowlesi, with a total of 524 occurrences included in the current study following consolidation with the 2015 study. The modelling framework used in the 2015 study was extended, with changes including the addition of new covariates to capture the effect of deforestation and urbanisation on P. knowlesi transmission.
DISCUSSION: Our map of P. knowlesi relative transmission suitability estimates that the risk posed by the pathogen is highest in Malaysia and Indonesia, with localised areas of high risk also predicted in the Greater Mekong Subregion, The Philippines and Northeast India. These results highlight areas of priority for P. knowlesi surveillance and prospective sampling to address the challenge the disease poses to malaria elimination planning.
METHODS: We assembled 1155 geographical records of yellow fever virus infection in people from 1970 to 2016. We used a Poisson point process boosted regression tree model that explicitly incorporated environmental and biological explanatory covariates, vaccination coverage, and spatial variability in disease reporting rates to predict the relative risk of apparent yellow fever virus infection at a 5 × 5 km resolution across all risk zones (47 countries across the Americas and Africa). We also used the fitted model to predict the receptivity of areas outside at-risk zones to the introduction or reintroduction of yellow fever transmission. By use of previously published estimates of annual national case numbers, we used the model to map subnational variation in incidence of yellow fever across at-risk countries and to estimate the number of cases averted by vaccination worldwide.
FINDINGS: Substantial international and subnational spatial variation exists in relative risk and incidence of yellow fever as well as varied success of vaccination in reducing incidence in several high-risk regions, including Brazil, Cameroon, and Togo. Areas with the highest predicted average annual case numbers include large parts of Nigeria, the Democratic Republic of the Congo, and South Sudan, where vaccination coverage in 2016 was estimated to be substantially less than the recommended threshold to prevent outbreaks. Overall, we estimated that vaccination coverage levels achieved by 2016 avert between 94 336 and 118 500 cases of yellow fever annually within risk zones, on the basis of conservative and optimistic vaccination scenarios. The areas outside at-risk regions with predicted high receptivity to yellow fever transmission (eg, parts of Malaysia, Indonesia, and Thailand) were less extensive than the distribution of the main urban vector, A aegypti, with low receptivity to yellow fever transmission in southern China, where A aegypti is known to occur.
INTERPRETATION: Our results provide the evidence base for targeting vaccination campaigns within risk zones, as well as emphasising their high effectiveness. Our study highlights areas where public health authorities should be most vigilant for potential spread or importation events.
FUNDING: Bill & Melinda Gates Foundation.
METHODS: We used macaque and mosquito species presence data, background data that captured sampling bias in the presence data, a boosted regression tree model and environmental datasets, including annual data for land classes, to predict the distributions of each vector and host species. We then compared the predicted distribution of each species with cover of each land class.
RESULTS: Fine-scale distribution maps were generated for three macaque host species (Macaca fascicularis, M. nemestrina and M. leonina) and two mosquito vector complexes (the Dirus Complex and the Leucosphyrus Complex). The Leucosphyrus Complex was predicted to occur in areas with disturbed, but not intact, forest cover (> 60% tree cover) whereas the Dirus Complex was predicted to occur in areas with 10-100% tree cover as well as vegetation mosaics and cropland. Of the macaque species, M. nemestrina was mainly predicted to occur in forested areas whereas M. fascicularis was predicted to occur in vegetation mosaics, cropland, wetland and urban areas in addition to forested areas.
CONCLUSIONS: The predicted M. fascicularis distribution encompassed a wide range of habitats where humans are found. This is of most significance in the northern part of its range where members of the Dirus Complex are the main P. knowlesi vectors because these mosquitoes were also predicted to occur in a wider range of habitats. Our results support the hypothesis that conversion of intact forest into disturbed forest (for example plantations or timber concessions), or the creation of vegetation mosaics, will increase the probability that members of the Leucosphyrus Complex occur at these locations, as well as bringing humans into these areas. An explicit analysis of disease risk itself using infection data is required to explore this further. The species distributions generated here can now be included in future analyses of P. knowlesi infection risk.
METHODOLOGY/PRINCIPAL FINDINGS: A total of 439 records of P. knowlesi infections in humans, macaque reservoir and vector species were collated. To predict spatial variation in disease risk, a model was fitted using records from countries where the infection data coverage is high. Predictions were then made throughout Southeast Asia, including regions where infection data are sparse. The resulting map predicts areas of high risk for P. knowlesi infection in a number of countries that are forecast to be malaria-free by 2025 (Malaysia, Cambodia, Thailand and Vietnam) as well as countries projected to be eliminating malaria (Myanmar, Laos, Indonesia and the Philippines).
CONCLUSIONS/SIGNIFICANCE: We have produced the first map of P. knowlesi malaria risk, at a fine-scale resolution, to identify priority areas for surveillance based on regions with sparse data and high estimated risk. Our map provides an initial evidence base to better understand the spatial distribution of this disease and its potential wider contribution to malaria incidence. Considering malaria elimination goals, areas for prioritised surveillance are identified.