Landscape drives zoonotic malaria prevalence in non-human primates
Abstract
Zoonotic disease dynamics in wildlife hosts are rarely quantified at macroecological scales due to the lack of systematic surveys. Non-human primates (NHPs) host Plasmodium knowlesi, a zoonotic malaria of public health concern and the main barrier to malaria elimination in Southeast Asia. Understanding of regional P. knowlesi infection dynamics in wildlife is limited. Here, we systematically assemble reports of NHP P. knowlesi and investigate geographic determinants of prevalence in reservoir species. Meta-analysis of 6322 NHPs from 148 sites reveals that prevalence is heterogeneous across Southeast Asia, with low overall prevalence and high estimates for Malaysian Borneo. We find that regions exhibiting higher prevalence in NHPs overlap with human infection hotspots. In wildlife and humans, parasite transmission is linked to land conversion and fragmentation. By assembling remote sensing data and fitting statistical models to prevalence at multiple spatial scales, we identify novel relationships between P. knowlesi in NHPs and forest fragmentation. This suggests that higher prevalence may be contingent on habitat complexity, which would begin to explain observed geographic variation in parasite burden. These findings address critical gaps in understanding regional P. knowlesi epidemiology and indicate that prevalence in simian reservoirs may be a key spatial driver of human spillover risk.
eLife assessment
This useful study presents findings regarding the impact of forest cover and fragmentation on the prevalence of malaria in non-human primates. The evidence supporting the claims of the authors is solid.
https://doi.org/10.7554/eLife.88616.4.sa0eLife digest
Zoonotic diseases are infectious diseases that are transmitted from animals to humans. For example, the malaria-causing parasite Plasmodium knowlesi can be transmitted from monkeys to humans through mosquitos that have previously fed on infected monkeys. In Malaysia, progress towards eliminating malaria is being undermined by the rise of human incidences of ‘monkey malaria’, which has been declared a public health threat by The World Health Organisation.
In humans, cases of monkey malaria are higher in areas of recent deforestation. Changes in habitat may affect how monkeys, insects and humans interact, making it easier for diseases like malaria to pass between them. Deforestation could also change the behaviour of wildlife, which could lead to an increase in infection rates. For example, reduced living space increases contact between monkeys, or it may prevent behaviours that help animals to avoid parasites.
Johnson et al. wanted to investigate how the prevalence of malaria in monkeys varies across Southeast Asia to see whether an increase of Plasmodium knowlesi in primates is linked to changes in the landscape. They merged the results of 23 existing studies, including data from 148 sites and 6322 monkeys to see how environmental factors like deforestation influenced the amount of disease in different places.
Many previous studies have assumed that disease prevalence is high across all macaques, monkey species that are considered pests, and in all places. But Johnson et al. found that disease rates vary widely across different regions. Overall disease rates in monkeys are lower than expected (only 12%), but in regions with less forest or more ‘fragmented’ forest areas, malaria rates are higher. Areas with a high disease rate in monkeys tend to further coincide with infection hotspots for humans. This suggests that deforestation may be driving malaria infection in monkeys, which could be part of the reason for increased human infection rates.
Johnsons et al.’s study has provided an important step towards better understanding the link between deforestation and the levels of monkey malaria in humans living nearby. Their study provides important insights into how we might find ways of managing the landscape better to reduce health risks from wildlife infection.
Introduction
Zoonotic infectious diseases arise from the spillover of pathogens into human populations, typically from a reservoir in wildlife hosts. Anthropogenic land use and land cover change have now been widely linked to infectious disease outbreaks (Brock et al., 2019; Davidson et al., 2019a; Loh et al., 2016). Such practices, including deforestation, logging, clearing for cash-crop plantations or conversion of intact forest into arable land, are accelerating across tropical forests of Southeast Asia (Fornace et al., 2021; Imai et al., 2018; Fornace et al., 2021; Imai et al., 2018). Mechanisms that underly the association between habitat disturbance and spillover risk from wildlife hosts are complex and occur over multiple spatial scales (Brock et al., 2019). In Brazil, re-emergence of Yellow Fever Virus in both NHPs and humans has been linked to areas with highly fragmented forest (Ilacqua et al., 2021). In part, an increase in ‘edge’ habitat in fragmented or mosaic landscapes can facilitate spatial overlap and altered contact patterns between wildlife, vectors, and humans (Lehman et al., 2006). Such ecological interfaces are also thought to contribute to parasite spillover in other vector-borne diseases including Zika (Li et al., 2021a), Babesiosis and Lyme disease (Simon et al., 2014), Trypanosoma cruzi (Vaz et al., 2007) and zoonotic malaria (Brock et al., 2019; Grigg et al., 2017). At the same time, habitat fragmentation can have detrimental impact on wildlife population viability, with reduced host species occupancy and reduced disease burden in highly disturbed habitats (Hanski and Ovaskainen, 2000). Disentangling this interplay is essential to inform ecological strategies for surveillance and mitigation of diseases in regions undergoing landscape change (Fornace et al., 2021).
Zoonotic P. knowlesi is a public health threat of increasing importance across Southeast Asia, following the identification of a prominent infection foci in Borneo in 2004 (Singh et al., 2004). P. knowlesi is a zoonosis, with a sylvatic cycle circulating in non-human primates (NHPs). Human cases currently occur only from spillover events (Ruiz Cuenca et al., 2022; Fornace et al., 2022; Fornace et al., 2023; Lee et al., 2011). Human transmission requires bites from infective mosquitos, primarily anopheline mosquitos of the Leucosphyrus Complex (Anopheles balabacensis, An. latens, An. introlactus) and Dirus Complex (An. dirus, An. cracens) (Moyes et al., 2016; Vythilingam et al., 2006; Wong et al., 2015). Natural hosts for P. knowlesi are typically Long-tailed macaques (Macaca fascicularis) and Southern Pig-tailed macaques (M. nemestrina) (Moyes et al., 2016), both occurring widely across Southeast Asia. Currently, distribution of P. knowlesi cases is thought to be restricted to the predicted ranges of known vector and host species (Davidson et al., 2019b), although recent studies have also identified other NHPs found to be harbouring P. knowlesi. This includes Stump-tailed macaques (M. arctoides), which are now considered to be another natural reservoir (Fungfuang et al., 2020).
Progress towards malaria elimination in Malaysia has been stymied by a recent rise in human incidence of P. knowlesi malaria. Even after accounting for increases in surveillance and diagnostic improvements it is now recognised as the most common cause of clinical malaria in Malaysia (Cooper et al., 2020). Indeed, Malaysia was the first country not to qualify for malaria elimination due to ongoing presence of zoonotic malaria and the WHO updated the guidelines to reflect zoonotic malaria as a public health threat (World Health Organization, 2021). Emergence of Plasmodium knowlesi infections has been linked to changes in land cover and land use (Fornace et al., 2021). While sporadic cases have been reported across Southeast Asia, including in Indonesia (Setiadi et al., 2016), the Philippines (Fornace et al., 2018), Vietnam (Maeno et al., 2015), Brunei (Koh et al., 2019), and Myanmar (Ghinai et al., 2017), the majority of P. knowlesi cases are found in East Malaysia (Borneo) with hotspots in the states of Sabah and Sarawak (Jeyaprakasam et al., 2020), areas that have seen extensive deforestation and landscape modification. In Sabah, human prevalence of P. knowlesi infection has recently been shown to be specifically associated with recent loss of intact forest, agricultural activities, and fragmentation across multiple localised spatial scales (Brock et al., 2019; Fornace et al., 2019b; Fornace et al., 2016).
Prevalence of the pathogen in reservoir hosts is one of three crucial factors determining the force of infection in zoonotic spillover events (Murray and Daszak, 2013). Despite this, very little is known of the impact of rapid landscape change on the distribution of P. knowlesi in NHPs. Literature on the impacts of fragmentation on primates tends to focus on primate density and abundance (Link et al., 2010; Zunino et al., 2007). What is known is that effects of land cover changes on primate-pathogen dynamics are highly variable and context-specific. Although the vector species responsible for sylvatic transmission remain unknown, the Anopheles leucospryphus group, the only vector group implicated in P. knowlesi transmission, is widely associated with secondary, disturbed forest (Brant, 2011; Hawkes et al., 2019; Wong et al., 2015). Macaques have been known to preferentially rely on fringe habitat, a behaviour that may be exaggerated in response to habitat fragmentation and facilitate exposure to vectors (Lehman et al., 2006; Stark et al., 2019). Changes to land composition can also create the biosocial conditions for higher rates of parasitism in primates. Under conditions of limited resources and reduction in viable habitat, conspecific primate density may increase as troops compete for available space. In turn, this can favour transmission via intra-species contact or allow the exchange of pathogens between troops dwelling in interior forest versus edge habitat (Faust et al., 2018; Stark et al., 2019). Habitat use may also become more intensive, preventing parasite avoidance behaviours (Nunn and Dokey, 2006). Land cover change is also known to favour more adaptable, synanthropic species such as M. fascicularis (McFarlane et al., 2012). Considering the spillover risk posed by wildlife reservoirs of P. knowlesi, clarifying any relationships between environmental factors and parasitaemia in key host species may contribute to a more comprehensive understanding of P. knowlesi transmission patterns.
Earth Observation (EO) data provides novel opportunities to investigate epidemiological patterns of diseases which are linked to environmental drivers (Kalluri et al., 2007). In relation to P. knowlesi, utility of fine-scale remote-sensing data has been demonstrated: examples include satellite-derived data used to examine household-level exposure risk in relation to proximate land configuration (Fornace et al., 2019b), UAV-imagery used to link real-time deforestation to macaque host behavioural change (Stark et al., 2019), and remote-sensing data used to interrogate risk factors for vector breeding sites (Byrne et al., 2021). Although macroecological studies that utilise geospatial data are often confounded by issues of matching temporal and spatial scales, as well as by the quality and accuracy of available georeferencing, measures can be taken to account for this when examining the role of environmental factors in modulating disease outcomes. Furthermore, ecological processes occur and interact over a range of distances, or ‘spatial scales’ (Brock et al., 2019; Fornace et al., 2016; Loh et al., 2016). This applies to determinants of vector-borne disease ecology, from larval breeding microclimate to wildlife host foraging behaviour. As multiple influential variables are rarely captured by a single scale (Cohen et al., 2016), data-driven methods can be applied to examine risk factors over multiple scales and identify covariates at their most influential extent (Byrne et al., 2021).
We hypothesise that prevalence of P. knowlesi in primate host species is spatially heterogeneous and that higher prevalence is partially driven by forest loss and fragmentation, contributing to the strong associations described between land use, land cover and human P. knowlesi risk. This study is the first to systematically assess P. knowlesi prevalence in NHPs at a regional scale, and across a wide range of habitats. In conceptual frameworks and transmission models, it is often assumed that P. knowlesi infections in NHPs are chronic (low level, persistent infection) and ubiquitous (uniformly distributed across populations; Brock et al., 2016; Jeyaprakasam et al., 2020). No studies have systematically assessed the extent and quality of all available data on P. knowlesi in NHPs. Independent studies investigating P. knowlesi in primates are typically constrained by small sample sizes and confined geographic areas, limiting inference that can be made about relationships between infection dynamics and landscape characteristics. Systematic tools developed for epidemiological studies of disease prevalence in human populations are rarely applied to the study of wildlife disease prevalence; however, such tools can be used to capture the scale and contrast required in macroecological studies to quantify disease burdens regionally. Furthermore, while recent research has shown the impact of deforestation on the distribution of macaques in the context of P. knowlesi (Moyes et al., 2016; Stark et al., 2019), associations between landscape and variation in the prevalence of simian Plasmodium spp. in primates have not been explored. We aimed to (1) assemble a georeferenced dataset of P. knowlesi in NHPs; (2) evaluate variation in NHP P. knowlesi prevalence by geographic region; and (3) assess environmental and spatial risk factors for P. knowlesi prevalence in NHPs across Southeast Asia.
Results
A systematic literature review was conducted in Medline, Embase, and Web of Science to identify articles reporting prevalence of naturally acquired Plasmodium knowlesi in NHPs. Twenty-three research articles were identified (Akter et al., 2015; Amir et al., 2020; Chang et al., 2011; Fungfuang et al., 2020; Gamalo et al., 2019; Ho et al., 2010; Zamzuri et al., 2020; Jeslyn et al., 2011; Kaewchot et al., 2022; Lee et al., 2011; Li et al., 2021b; Muehlenbein et al., 2015; Nada-Raja et al., 2022; Putaporntip et al., 2010; Saleh Huddin et al., 2019; Salwati et al., 2017; Seethamchai et al., 2008; Shahar, 2019; Gilhooly et al., 2016; personal correspondence, 2013 and 2015 ; Vythilingam et al., 2008; Yusuf et al., 2022; Zhang et al., 2016), containing 148 unique primate survey records to form the dataset for analyses (see Appendix 2 for details of JBI Critical Assessment) (Munn et al., 2015). Year of sampling ranges from 2004 to 2019. No primatological studies were identified from Vietnam, Brunei, or Timor-Leste. Full characteristics of the articles and individual study methodologies are reported in Appendix 1—table 2. Spatial resolution of the survey sites varied from GPS point coordinates to country-level administrative boundaries (Appendix 5—table 1). Geographic distribution of sampling is illustrated in Figure 1.
Overall, records report on a total of 6322 primates, with the largest proportion sampled from Peninsular Malaysia (48.5%, n=3069/6322). Primate surveys were primarily conducted on Long-tailed macaques (Macaca fascicularis) (90.5%, n=5720/6322) followed by Pig-tailed macaques (M. nemestrina; n=532/6322; Amir et al., 2020; Lee et al., 2011; Muehlenbein et al., 2015; Putaporntip et al., 2010; Appendix 1—table 3). Reported prevalence of Plasmodium knowlesi in NHPs ranged from 0% to 100%. Only 87 of the surveys (58.8%, n=87/148) reported a positive diagnosis, with the remaining 61 sites finding no molecular evidence of P. knowlesi infection (41.2%) in any primates tested. A full breakdown of P. knowlesi infection rates according to reported primate characteristics can be found in Appendix 1.
Meta-analysis of P. knowlesi prevalence
To quantify regional heterogeneity in simian cases of P. knowlesi, a one-stage meta-analysis of prevalence (number positive out of the number sampled) was conducted on primate malaria survey data. Overall pooled estimate for P. knowlesi prevalence was 11.99% (CI95% 9.35–15.26). Overall heterogeneity was assessed using the I2 statistic. Substantial between-study heterogeneity (I2 ≥75%) was found across all prevalence records (I2=80.5%; CI95% 77.3–83.1). In the sub-group analysis by region, pooled prevalence estimates are consistently low for Thailand (2.0%, CI95% 1.1–3.5%), moderate in Peninsular Malaysia (14.3%, CI95% 11.1–18.2) and elevated in Singapore (23.3%, CI95% 11.0–42.8) and Malaysian Borneo (41.1%, CI95% 20.8–64.9) (Figure 2). Sub-group heterogeneity was assessed using prediction intervals, derived from τ 2 statistic used to describe between-study variability. Prediction intervals indicate high heterogeneity of estimates within regions, consistent with expectations of high variability of prevalence across individual study sites. Detailed forest plots for individual prevalence estimates can be found in Appendix 3—figure 2.
Risk factor analysis
Covariate data and P. knowlesi prevalence data were used to fit additional models to explore the relationships between localised landscape configuration and NHP malaria prevalence. Environmental covariates were extracted from satellite-derived remote sensing datasets (Table 1) at either true sampling sites (GPS coordinates) or 10 random pseudo-sampling sites to account for geographic uncertainty in prevalence data. Host species was grouped as ‘Macaca fascicularis’ or ‘Other’ due to sample counts of <10 for certain primate species. Only 57.4% (n=85/148 records) of data included year of sampling, deemed to be insufficient to assess temporal patterns in prevalence. Tree canopy cover ranged from negligible to near total cover (100%) within buffer radii (Appendix 4—table 2). Details of covariate data processing is illustrated in Appendix 4.
Following a two-stage approach for selection of explanatory variables, tree cover and fragmentation (measured by perimeter: area ratio, PARA) were retained at 5 km as linear terms, human population density was retained at both 5 km and 20 km and primate species was retained as a categorical variable. Spearman’s rank tests for residual correlation between final variables at selected scales indicates a strong negative correlation between tree cover and fragmentation index (PARA; ρ = –0.75; Appendix 6—figure 2).
Adjusting for all other covariates in the model, we identified strong evidence of an effect between increasing tree canopy cover and higher prevalence of P. knowlesi in NHPs within a 5 km radius (aOR = 1.38, CI95% 1.19–1.60; p<0.0001). Evidence was also found for an association between likelihood of P. knowlesi and higher degrees of habitat fragmentation (PARA) within 5 km (aOR = 1.17, CI95% 1.02–1.34, p<0.0281). Evidence suggests that human population density within a 5 km radius is associated with risk of P. knowlesi in NHP (aOR = 1.36, CI95% 1.16–1.58, p=0.0001) whilst human density within 20 km has an inverse effect on likelihood of P. knowlesi (aOR = 0.56, CI95% 0.46–0.67, p<0.0001). M. fascicularis is also associated with higher prevalence relative to all other non-human primate species (aOR = 2.50, CI95% 1.31–4.85; p=0.0051). Additional complexity did not improve optimal model fit and effect modification was not pursued. In sensitivity analyses removing data points with excessive spatial uncertainty or restricting data points only to areas with high probability of macaque occurrence, evidence was consistently found that tree canopy cover (5 km) and host species exhibit a strong positive association with prevalence of P. knowlesi in NHP (Appendix 6). Final adjusted OR for the full multivariable model can be visualised in Figure 3.
Discussion
Land use and land cover change is widely linked to spillover of zoonotic pathogens from sylvatic reservoirs into human populations, and pathogen prevalence in wildlife host species is key in driving the force of infection in spillover events. Our initial analyses found that for Plasmodium knowlesi, there is substantial spatial heterogeneity and prevalence in non-human primates varies markedly between regions of Southeast Asia (Zhang et al., 2016). Consistent with our hypothesis that parasite density in primate hosts would be higher in areas experiencing habitat disturbance, we identified strong links between P. knowlesi in NHPs and measures of contemporaneous tree cover and habitat fragmentation. To our knowledge, this is the first systematic study to find evidence of landscape influencing the distribution of P. knowlesi prevalence in NHPs. Results offer evidence that P. knowlesi infection rates in NHPs are linked to changes in landscape across broad spatial scales, and that prevalence of P. knowlesi in reservoir species may be driving spillover risk across Southeast Asia. These findings could provide insight to improving surveillance of P. knowlesi and to the development of ecologically targeted interventions.
While previous studies have estimated that P. knowlesi infection would be chronic in all macaques, or as high as 50–90% for modelling P. knowlesi transmission in Malaysia (Brock et al., 2016), this data strongly suggests that this is not the case. Overall prevalence of P. knowlesi infection in all NHPs is markedly lower than usual estimates, emphasising the importance of accounting for absence data in estimations of prevalence. Considerable heterogeneity was identified between and within regional estimates for P. knowlesi across Southeast Asia, which likely reflects genuine differences according to distinct climates and habitats (Shearer et al., 2016). Malaysian Borneo was found to have an estimated prevalence over five-fold higher than West Malaysia. Crucially, such extreme prevalence estimates for NHPs in Borneo align with the known hotspot for human incidence of P. knowlesi (Cooper et al., 2020). By comparison, for Peninsular Malaysia, estimated prevalence is far lower than anticipated. Cases of human P. knowlesi do occur in West Malaysia, although transmission has been found to exhibit spatial clustering (Phang et al., 2020) which may correspond to pockets of high risk within the wider context of low prevalence of P. knowlesi in macaque populations. Regional trends in P. knowlesi also mask differences in infection rates between sample locations, driven by more localised factors. Multiple studies reported finding P. knowlesi infections in wild macaques to be low or absent in peri-domestic or urbanised areas, attributed to the absence of vector species typically found in forest fringes (Brant et al., 2016; Chua et al., 2019; Manin et al., 2016). This pattern is seen in reports from Peninsular Malaysia (Saleh Huddin et al., 2019; Vythilingam et al., 2008), Singapore (Jeslyn et al., 2011; Li et al., 2021b), and Thailand (Fungfuang et al., 2020; Putaporntip et al., 2010). The high heterogeneity of reports here suggests that the picture is even more complex. P. knowlesi infections may even vary between troops within a single study site, as was seen in the Philippines (Gamalo et al., 2019). Fine-scale interactions are unlikely to be captured by the scale of this study.
Ecological processes determining P. knowlesi infection are influenced by dynamic variables over multiple spatial scales (Cohen et al., 2016). We utilised a data-driven methodology to select variables at distances that capture maximum impact on P. knowlesi prevalence (Byrne et al., 2021; Fornace et al., 2019b), with tree cover and fragmentation influential at localised scales and human population density also exerting influence within wider radii. Contrary to previous studies on risk factors for human incidence of P. knowlesi (Fornace et al., 2019b; Fornace et al., 2016), elevation was not found to be associated with P. knowlesi in NHPs at any scale. Vector and host species composition vary substantially across tropical ecotones, and it is likely that the study extent encompasses a range of putative vectors across different landscapes, such as those of the Minimus Complex in northern regions (Parker et al., 2015) or the recently incriminated An.-collessi and An.-roperi from the Umbrosus Group (De Ang et al., 2021). Given that the vector species driving sylvatic transmission remain elusive, it is conceivable that the elevation range covers multiple vector and host species niches and explains the lack of observed relationship between elevation and P. knowlesi in NHPs. Human population density was found to be significant at multiple distances, with contrasting effects on parasite prevalence in NHP. Previous studies have found a negative association between human density and vector density and biting rates in forested landscapes (Fornace et al., 2019a). Across wide spatial scales, increased vector density in less populated, more forested areas could generate higher parasite prevalence in NHPs. At the same time Long-tailed macaques, a species shown here to have higher prevalence rates, are notorious as nuisance animals and many of the available samples were collected opportunistically in urban areas, which might underly the observed positive association between localised high human density and higher prevalence in NHP. Whilst more data would be needed to understand this interaction, this further demonstrates the importance of using approaches to identify disease dynamics across multiple spatial scales (Brock et al., 2019).
A key finding is the link between high prevalence of P. knowlesi in primate host species with high degrees of habitat fragmentation. Habitat fragmentation is a key aspect of landscape modification, where large contiguous areas of habitat (for example, forests) are broken into a mosaic of smaller patches. This disturbs the ecological structure by increasing the density of fringes or ‘edges’, dynamic habitat often at the boundaries between natural ecosystems and human-modified landscapes (Borremans et al., 2019). Other studies have linked habitat fragmentation to increased generalist parasite density in primates. In Uganda, a higher prevalence and infection risk of protozoal parasites was observed in wild populations of red colobus primates (Procolobus rufomitratus) inhabiting fragmented forests compared to those in undisturbed habitat (Gillespie and Chapman, 2008). For P. knowlesi, creation of edge habitat is thought to favour vectors of the Leucosphyrus Complex (Davidson et al., 2019a; Hawkes et al., 2019). Anopheles spp. presence can be predicted by indices of fragmentation in Sabah, Borneo, with land cover changes creating more suitable micro-climate for larval habitats (Byrne et al., 2021), and an increased abundance of An. balabacensis found in forest fringes (Hawkes et al., 2019; Wong et al., 2015). Increasing landscape complexity results in increased density of edge habitat, with conceivably higher density of vectors in forest fringes. Therefore, preferential use of fringe habitat and high exposure to vectors in forest fringes may contribute to higher conspecific transmission of P. knowlesi between primates in increasingly fragmented habitats. This finding also lends clarity to landscape fragmentation as a risk factor for human exposure to P. knowlesi in Malaysian Borneo (Brock et al., 2019; Fornace et al., 2019b), with changes in relative host density, vector density and wildlife parasite prevalence in nascent forest fringes potentially enhancing the spillover of this disease system into human populations in fragmented habitats.
Conversely, we saw a strong association between high parasite prevalence and high tree canopy coverage. Given that a strong inverse relationship with fragmentation was observed, with high tree density correlating to low fragmentation indices and vice versa, this speaks to a trade-off between dense canopy cover and high habitat complexity and suggests an ‘ideal’ amount of habitat fragmentation that facilitates prevalence in primate hosts. For animals with larger home ranges, individual-based disease models combined with movement ecology approaches have shown that the most highly fragmented areas are less favourable for maintaining parasite transmission (White et al., 2018). In Sabah, individual macaques were shown to increase ranging behaviour in response to deforestation (Stark et al., 2019). Forest edge density also peaks at intermediate levels of land conversion (Borremans et al., 2019). With smaller habitat patches in maximally fragmented landscapes potentially insufficient to support macaque troops, this interplay between disease ecology and metapopulation theory may explain why both tree density and habitat fragmentation appear to pose a greater risk for simian P. knowlesi. Likewise, this may relate to the finding that in Borneo, larger forest patches (lower fragmentation indices) were associated with P. knowlesi spillover in Borneo (Fornace et al., 2019b). Overall, this finding offers an insight to mechanisms that underpin the increased force of infection of P. knowlesi that is associated with landscape change.
There are limitations to consider in the available data and interpretation of these findings. ‘Small-study effects’ were observed in the dataset, suggestive of a bias toward positive effect estimates (Stewart et al., 2012). This may be a result of data disaggregation and small studies creating artefactually higher estimates or may reflect true bias in data collection toward areas known to be endemic for P. knowlesi and convenience sampling of macaques. Assumptions have also been made that sample site equates to habitat, which may not reflect actual habitat use, and even accurate georeferenced data points are unlikely to entirely reflect surrounding habitat within the macaque home range. Variability in study designs and data reporting also impacted geospatial accuracy. Steps were taken to account for spatial bias by extracting covariates at randomly generated pseudo-sampling points. Whilst uncertainty cannot be eliminated, we demonstrate a robust methodology to accommodate for geographical uncertainty in ecological studies. Future investigations should prioritise systematic, georeferenced sampling across a range of landscape scenarios.
Results show important regional ecological trends, but broad geographic patterns may not be generalisable at individual levels, or to all putative host species in all geographic contexts (Zhang et al., 2016). Follow-up studies should be conducted at higher spatial and temporal resolution to characterise the effect of local landscape configuration on wildlife P. knowlesi prevalence. Effects of fragmentation are likely to be dependent on land conversion type, species composition and surrounding matrix habitat (Fornace et al., 2019b). Use of perimeter: area ratio (PARA) as a fragmentation index was justified given high canopy coverage in study sites (Wang et al., 2014), although Edge Density (ED) or normalised Landscape Shape Index (nLSI) might be more appropriate in future analyses to account for variation in forest abundance. Specific land configurations have previously been linked to P. knowlesi exposure in Borneo (Fornace et al., 2019b), notably in areas where palm oil plantation is a dominant industry. Given this, broad forest classifications used here may mask important differences in P. knowlesi prevalence between land classes. As it was not possible to include contemporary land cover classifications in this analysis, future studies would also benefit from looking at specific habitat type (e.g., primary forest, agroforest, plantation).
Concluding remarks
Strong links have been identified between land use and land cover change and ecosystem perturbation that favours the transmission of vector-borne diseases (Loh et al., 2016). Prevalence of P. knowlesi in macaques is likely to be a crucial determinant of human infection risk, and more representative estimates of P. knowlesi prevalence derived here can better inform regional transmission risk models. This study also characterises landscape risk factors for heightened prevalence of P. knowlesi in NHPs. Findings provide evidence that P. knowlesi in primate hosts is partly driven by landscape modification across Southeast Asia. While the full complexity is not captured by the covariates used, it is clear that P. knowlesi infection in NHPs is not restricted to densely forested areas. This study also demonstrates the utility of systematic meta-analysis tools and remote-sensing datasets in the investigation of macroecological disease trends, in conjunction with methods to standardise a spatially heterogeneous dataset and data-driven selection of spatial scales. Gaps identified in data reporting should inform more systematic and localised primatological surveys to disentangle precise mechanisms. Notwithstanding limitations, this study highlights the marked spatial heterogeneity and role of landscape complexity in driving P. knowlesi infection rates in NHPs. Given the clear intersection between human epidemiology and wildlife ecology, it is essential that infection dynamics within wildlife reservoirs are considered in future public health interventions.
Methods
Study site
This study focused on the simian malaria Plasmodium knowlesi across Southeast Asia, within 28°30'00.0"N, 92°12'00.0"E and 11°00'00.0"S, 141°00'00.0"E. Climate mainly corresponds to the equatorial tropical zone, with high temperatures and high humidity.
Data assembly
A systematic literature review was conducted under the CoCoPop framework (Condition, Context, Population) (Ruiz Cuenca et al., 2022; Munn et al., 2015). All studies identified in the literature review were screened for data on NHPs with a confirmed P. knowlesi diagnosis or absence data (zero counts of P. knowlesi with appropriate diagnostic methods). Exclusion criteria included (a) studies exclusively relying on microscopy (Antinori et al., 2013) (b) laboratory, animal model or experimental infection studies (c) data from outside of Southeast Asia. No limit was set on the temporal range for primate survey records. Duplicate records reporting results from the same surveys were removed, with one record per survey retained. Critical appraisal of the studies was conducted using the Joanna Briggs Institute (JBI) checklist for prevalence studies (Munn et al., 2015; see Appendix 1 for details and criteria). A flowchart of the selection process is illustrated in Appendix 1—figure 3, with a full list of articles included provided in Appendix 1—table 2.
Primary outcome was defined as P. knowlesi prevalence (p, proportion positive for P. knowlesi infection from n sampled NHPs). For each independent primate study, the following variables were extracted: year of data collection, primate species sampled, primate status (wild/captive), diagnostic test (PCR/sequencing) and target gene(s), sampling method (routine/purposive), number of P. knowlesi positive samples, number of Plasmodium spp. positive samples, total number of primates tested and geographical information.
In most studies identified, study site was only geolocated to a geographic area or descriptive location. Geolocation was assigned at the lowest available level of administrative polygon (i.e. district/state/country) by cross-referencing reported sampling location with GADM (v3.6) administrative boundaries. If specific location was given, GPS coordinates were assigned via Google Maps. For data visualisation, point coordinates were plotted in QGIS (3.10.14) and R (4.1.0) software.
Meta-analysis of P. knowlesi prevalence
Meta-analysis was conducted using methods that are standard in the analysis of human disease prevalence for individual participant datasets (IDP) (Liberati et al., 2009; Stewart et al., 2012). Data were disaggregated by geographic location (site) and primate species, to illustrate variance in prevalence by survey unit (Stewart et al., 2012). One-stage meta-analysis is considered appropriate for studies where the outcome may be infrequent, so data was included in a single model under the ‘DerSimonian and Laird’ variance estimator (Munn et al., 2015). Sensitivity analyses were conducted to compare methods for the back-transformation of prevalence estimates. For studies where prevalence estimates tend towards 0% or 100%, variance tends towards 0. To stabilise the variance and enable back-transformation of zero prevalence records, logit method was selected for the transformation of prevalence, with the inverse variance method used for individual study weights (see Appendix 3 for details).
Overall heterogeneity of prevalence records was assessed using the I2 statistic (von Hippel, 2015), a relative estimate of true between-study variance. Sub-group analysis was conducted according to geographic region, with the heterogeneity of reported prevalence within regional sub-groups assessed using prediction intervals derived from the τ 2 statistic. Small-study effects, including selection and publication biases, were assessed by examining funnel plots and imputing ‘missing’ estimates using the trim-and-fill method (Lin and Chu, 2018). Full rationale and details of small-study effect assessments can be found in Appendix 3.
Remote sensing data
Satellite-derived remote sensing datasets were used to assemble local environmental and anthropogenic covariates. Gridded UN-adjusted human population estimates were assembled at 1 km resolution from WorldPop, 2018. Elevation data was obtained from NASA SRTM 90 m Digital Elevation Database v4.1 (CGIAR-CSI) (Jarvis et al., 2008) with a spatial resolution of 1 km. Contemporaneous tree cover was derived from Hanson’s Global Forest Watch (30 m) (Hansen et al., 2013), extracted for every year between 2006 and 2020.Tree cover was classified as ≥50% crown density, and then matched to primate data by sample site geolocation and by year of sample collection to account for rapid forest loss (Appendix 4—figure 1). Where a broad timeframe of sampling was provided (≥3 years), median year was used. Full details for variable selection and processing can be found in Appendix 4.
Perimeter: area ratio (PARA, ratio of patch perimeter length to patch surface area) of given land class is a key metric for habitat conversion, where a higher PARA provides a measure of boundary complexity and indicates a more fragmented landscape (McGarigal et al., 2021). Mean PARA was extracted from canopy cover within circular buffers. Habitat fragmentation has been shown to correlate with disease transmission parameters (Borremans et al., 2019; Faust et al., 2018), but definitions often lack precision and can be considered with respect to ‘separation effects’ (division and isolation of patches) and ‘geometric effects’ (changes to ratios of perimeter and core habitat; Wilkinson et al., 2018). PARA provides a measure of edge density within the buffer area (PARA >0) and has been shown to provide a good index of fragmentation and good discrimination of spatial aggregation across areas where habitat abundance (tree canopy cover) is high (Wang et al., 2014; Appendix 4—table 2, Appendix 4—figure 4).
Covariate assembly
For studies with exact GPS coordinates, precise environmental data at a single site could be obtained. For surveys published without GPS coordinates, there is considerable geographic uncertainty in the exact sampling location (Appendix 5). Uncertainty in the spatial and environmental determinants of prevalence generates a sampling bias, with the precision of covariates correlated to certain studies. Use of a single centroid proxy site is standard procedure, but often generates erroneous estimates in large or heterogenous sampling units (Cheng et al., 2021). Alternative strategies were employed to account for and mitigate the effect of spatial uncertainty and spatial bias. Each prevalence observation was replicated and assigned a random sample of environmental realisations. 10 random sampling points were generated within the sampling area provided by the study, and covariates were extracted at each proxy sampling site (Appendix 5—figure 1). Selection of random points was validated by visual inspection of the stability of model coefficients with the inclusion of an increasing number of points. Number of points was selected conservatively at the point where coefficients stabilised (n=10).
For every georeferenced sampling point, mean values for all selected covariates were extracted within buffer radii at 5 km, 10 km, and 20 km (Appendix 4). Buffer area sizes were selected to investigate multiple spatial scales over which associations between risk factors and P. knowlesi prevalence might occur. A minimum radius of 5 km was chosen to approximate the maximum ranging distance for M. fascicularis (Waxman et al., 2014), with wider radii (10–20 km) included to account for the geographic uncertainties in areal data. Flowchart of data processing chain can be found in Appendix 4—figure 2.
Analysis of environmental risk factors
Generalised linear mixed-effect regression models (GLMM) were fitted to NHP prevalence data using a binomial distribution with a logit link. To account for within-study correlation in reported average prevalence, a unique identifier combining author and study was included as a random intercept in all models. Artificial inflation of sample size in the replicated data (10 pseudo-sampling sites for data geolocated to administrative areas) was accommodated by reducing individual observation weights to 1/10th within the model.
Each covariate at each spatial scale was assessed for inclusion in the multivariable model based on bivariable analysis and a criterion of p>0.2 under likelihood ratio tests (LRT; Appendix 6—table 1). A quadratic term for the fragmentation index ‘PARA’ was included to account for possible nonlinearity. Multicollinearity among independent predictors at multiple scales was examined via variance inflation factors (VIF). The VIF of each predictor variable was examined following a stepwise procedure, starting with a saturated model and sequentially excluding the variable with the highest VIF score from the model. Stepwise selection continued in this manner until the entire subset of explanatory variables in the global model satisfied a moderately conservative threshold of VIF ≤6 (Rogerson, 2001). Qualifying variables obtained were then assessed for model inclusion using a backward stepwise strategy, removing variables with the highest p value (LRT) until a pre-defined threshold of α<0.05. Spearman’s rank tests were conducted on the selected variables to observe residual correlation, plotted as a correlation matrix (Appendix 6—figure 2).
Fully adjusted odds ratios (OR) for associations between environmental covariates and P. knowlesi prevalence were derived from the final multivariable GLMM with p values derived from LRT. Spatial sensitivity analyses were conducted by excluding data points from administrative boundaries outside a reasonable size or above a reasonable threshold of environmental certainty, according to the standard deviation (SD) of the covariate values within each set of 10 environmental realisations (Appendix 5—figure 2, Appendix 6—tables 4–6). Ecological sensitivity analyses were conducted by removing data points that fall outside areas with high predicted probability of occurrence for Macaca fascicularis, Macaca nemestrina, and Macaca leonina and running regression analyses on the constrained dataset (Moyes et al., 2016; Appendix 6—figures 3–5, Appendix 6—tables 7–9).
Appendix 1
Data assembly
Prior to conducting the study, a review of current literature was constructed to find articles related to ‘Plasmodium knowlesi’ or to both ‘malaria’ and ‘primate’, including synonyms and sub-headings. The search was elaborated to specify environmental factors (Appendix 1—figure 1). The following databases were searched:
Medline
Embase
Web of Science
Provisional data were extracted using a standardised form (Appendix 1—figure 2) using standardised definitions (Appendix 1—table 1), from which an initial set of studies were identified for this investigation. Duplicate records (confirmed/suspected to be the same specimens) were removed, with one record retained (Appendix 1—figure 3)
Of the 87 records reporting presence of P. knowlesi, only 22 records (containing 248 P. knowlesi positive macaques) report whether P. knowlesi infection was a mono-infection or mixed infection with other simian Plasmodium spp. With a low proportion of data represented, this was deemed insufficient to conduct any further investigations.
Macaca fascicularis is the predominant species tested. However, reports also include M. nemestrina (6.1%, n=9/148; 527 macaques) (Amir et al., 2020; Lee et al., 2011; Putaporntip et al., 2010; Muehlenbein et al., 2015), M. arctoides (1.4%, n=2/148; 36 macaques) (Fungfuang et al., 2020; Putaporntip et al., 2010), M.-leonina (n=1/148; 25 macaques) (Fungfuang et al., 2020), Trachypithecus obscurus (Dusky leaf monkey) (n=1/148; 7 tested) and unspecified species from the Presbytis genus (n=1/148; 2 tested) (Appendix 1—table 3). One study additionally sampled 1 Presbytis melalophos (Black-crested Sumatran langur) (Vythilingam et al., 2008), but species-specific P. knowlesi was not reported (Appendix 1—table 4).
Appendix 2
Quality appraisal
Quality was assessed using the Joanna Briggs Institute (JBI) Critical Appraisal tool for prevalence studies (Munn et al., 2015). Studies were assessed on nine standardised criteria used to inform inclusion in the meta-analysis. Full criteria and examples of scoring are given in Appendix 2—tables 1 and 2. Studies assessed to be of lower quality were those that omitted key information about the sampling method . Two studies were considered to be of higher quality owing to robust sampling and completeness of evidence (Lee et al., 2011; Saleh Huddin et al., 2019). Given the objective to assess variation in reported prevalence, and in considering the limited usefulness of criteria designed for human participants, the reliable diagnostic methods identified in all studies and the appreciable limitations in surveying wild animals, data from all studies (n=148 estimates) were included for further analyses.
Appendix 3
Meta analysis
Sub-group analysis by region was conducted under a random-effects model. Pooled estimates are then back-transformed for interpretation. The Freeman-Tukey double arcsine method is recommended in the transformation of prevalence (Munn et al., 2015). However, recent studies have found that back-transformation of the Freeman-Turkey method can generate misleading results, owing to the requirement for a global sample size for inversion (Schwarzer et al., 2019). A sensitivity analysis conducted using the logit transformation and untransformed proportions (Appendix 3—table 1) revealed a deficit in the back-transformation of the pooled prevalence estimate for Thailand under the Freeman-Tukey transformation, generating a null point estimate. To avoid this error, meta-analysis was conducted under the logit transformation with the inverse variance estimator to account for individual study weighting.
Small study effects, including selection and publication bias, were assessed by examining funnel plots and imputing ‘missing’ estimates using the trim-and-fill method (Lin and Chu, 2018). Funnel plots were generated using both SE and study size as metrics of variance, as study size has been shown to be more accurate for meta-analyses of proportions where raw estimates tend towards 0 or 1 (Munn et al., 2015). Funnel plots for the disaggregated dataset are shown using SE as the variance estimate in Appendix 3—figure 1. Asymmetry is highlighted by the trim-and-fill interpolation method, which provides an estimate of missing data and added an additional 54 imputed points to the plot.
Meta-analysis was conducted with data disaggregated by survey location and primate species (k=148). Forest plot of individual study prevalence, presented with pooled regional prevalence estimates and relative sampling effort for the disaggregated data can be visualised in Appendix 3—figure 2.
Appendix 4
Remote sensing data and covariate assembly
Environmental covariates were extracted from satellite-derived remote-sensing datasets. Elevation can be used as a proxy for vector range, with malaria transmission patterns often correspond to altitudinal ranges of vector species (Hawkes et al., 2019). Estimates of human population density provide a measure of urbanisation, used to examine risks related to human settlement proximity. Metrics of deforestation including canopy cover (cumulative loss) and degree of fragmentation, which are key determinants of macaque habitat selection (McGarigal et al., 2021) and of mosquito vector breeding sites and were derived from Hansen’s Global Forest Watch (Hansen et al., 2013) Land use classification maps derived from the Intact Forest Landscape (IFL) (Potapov et al., 2008) and Copernicus Global Land Service (100 m) (Buchhorn et al., 2020) were also explored to provide more detailed information on specific landscape composition. However, as 77.0% of NHP records were collected before the earliest land classification in 2015 (N=114/148), datasets were of limited utility and not pursued further.
Gridded UN-adjusted human population estimates were assembled at 1 km resolution from WorldPop, 2018 for multiple timepoints as a measure of urbanisation, a proxy for risks related to human settlement proximity. As minimal variation between timepoints was observed, only 2012 (median year of primate data collection) was retained.
Elevation, population density and forest cover all varied markedly across surveyed sites. Forest cover ranges from negligible to near total cover within 5 km, and up to 99.96% and 99.64% within 10 km and 20 km respectively (Appendix 4—table 2, Appendix 4—figure 4). Within a 20 km buffer, 46.1% of sites have dense forest cover ≥50% (n=683/1480) and 83.85% have moderate or high forest cover (≥20%), with similar distributions over 5 km and 10 km. Example buffers over forest cover data can be visualised in Appendix 4—figure 5.
Appendix 5
Spatial uncertainty
Available spatial resolution of the survey sites varied. 14 records (9.5%, n=14/148) could be geolocated to a point, using geographic coordinates provided or inferred. The remaining 134 were geolocated to the lowest administrative polygon according to GADM boundary definitions (Appendix 5—table 1).
Crude sensitivity analyses were initially conducted to evaluate use of centroids vs random points to approximate macaque survey site. Using GADM classifications, the largest polygon containing NHP data was identified at each administrative level. 10 points were randomly generated within each polygon in QGIS, with buffers at 5/10/20 km. Proportion of forest cover per buffer was extracted, categorised and compared to the forest cover for the centroid (Appendix 5—figure 1).
Results indicate that at district level (GID2), minimal change in forest cover between the points was observable (5 km: 0.68–0.96). However, for both the state of Southern Sumatra (GID1 5 km: 0.26–0.97) and for Myanmar (GID0 5 km: 0.16–0.99) variation was observed, with several points classifying as ‘moderate’ (20–50% cover) rather than ‘high’ (>50%) as suggested by the centroid. Overall, results show a disparity between covariates obtained at a central point compared with random points within larger administrative polygons, indicating inadequate sensitivity of centroids as a proxy for local landscape variables where there is spatial uncertainty (Cheng et al., 2021).
Appendix 6
Regression analysis
For observations with high geographic uncertainty, random pseudo-sampling of 10 sites (as described in the manuscript and Appendix 4—figure 2) was used to avoid overgeneralisations and biases typical when centroids are used as proxy sampling sites35, with data weighted accordingly. However, this has the potential to introduce extreme or unrealistic values by generating points in landscapes that are outside reasonable estimations of primate study site. Given this, further sensitivity analyses were conducted to validate the results against geographic precision of observations. Data were first truncated to include only data geolocated to administrative boundaries for relatively small area size (see Appendix 6—table 3) and exclude highly variable data from country-level boundaries (GID0). Singapore was retained as small administrative unit. GLMM regression models were fit to the truncated dataset.
Administrative boundaries are arbitrary categories and vary considerably in size and landscape consistency. To better evaluate environmental uncertainty associated with each observation, standard deviation (SD) of the covariate values within each set of 10 environmental realisations was calculated (resulting in a single standard deviation value for each covariate at each scale for every prevalence data point). Studies for which the uncertainty (SD) of covariates exceeded half of the maximum standard deviation were censored to avoid spurious associations derived from unreliable/extreme values for both forest cover (5 km) and fragmentation (5 km). Regression models were fit to the winsorized dataset and compared to results from the full dataset to ensure that associations are robust.
Sensitivity analysis has shown that the trends are robust when the data is constrained according to small administrative boundaries or by measures of spatial uncertainty in the environmental variables. However, given that a proportion of points are randomly generated, questions remain about how suitable the resulting sites are for macaque species and consequently whether the associations observed are a realistic indication of ecological trends. To address this, points were subset according to macaque species habitat suitability maps, derived from Moyes et al., 2016.
Predicted occurrence maps were combined for three species Macaca fascicularis, Macaca nemestrina and Macaca leonina to create a joint macaque extent for Southeast Asia. Binary maps of predicted habitat extent were then generated using thresholds of moderate (predictions of 0.5 and above), high (>0.75) and very high predicted probability of occurrence (>0.9) (Moyes et al., 2016; Appendix 6—figures 3–5). Datapoints from the main analysis were then overlayed with each map, and any points (±5 km buffer) that occurred outside of predicted habitat extent for macaque species were removed. Regression models were then fit to the reduced datasets to assess whether associations observed are plausible according to macaque ecology (Appendix 6—tables 7–9).
Data availability
All data of P. knowlesi prevalence in non-human primates were sourced from systematic literature review of publicly available research articles. Published primate datasets are available in Ruiz Cuenca et al., 2022 and Shahar, 2019. Full reference list for all publicly available primate infection data including details of individual articles and data types are available in Appendix 1. Earth Observation (EO) data were derived from open-source environmental datasets. References, links and details of scale and resolution and manipulation are provided in manuscript (Table 1 in Results and details in Methods) with further information provided in the Appendix 4. Example code for manipulation of environmental data and data analysis are available on an associated GitHub repository (copy archived at Johnson, 2024).
References
-
Cross-species pathogen spillover across ecosystem boundaries: mechanisms and theoryPhilosophical Transactions of the Royal Society B 374:20180344.https://doi.org/10.1098/rstb.2018.0344
-
BookChanges in Abundance, Diversity and Community Composition of Mosquitoes Based on Different Land Use in Sabah, MalaysiaImperial College London.
-
Predictive analysis across spatial scales links zoonotic malaria to deforestationProceedings. Biological Sciences 286:20182351.https://doi.org/10.1098/rspb.2018.2351
-
Identification of Hepatocystis species in a macaque monkey in northern MyanmarResearch and Reports in Tropical Medicine 2:141–146.https://doi.org/10.2147/RRTM.S27182
-
Pathogen spillover during land conversionEcology Letters 21:471–483.https://doi.org/10.1111/ele.12904
-
Association between Landscape Factors and Spatial Patterns of Plasmodium knowlesi Infections in Sabah, MalaysiaEmerging Infectious Diseases 22:201–208.https://doi.org/10.3201/eid2202.150656
-
Exposure and infection to Plasmodium knowlesi in case study communities in Northern Sabah, Malaysia and Palawan, The PhilippinesPLOS Neglected Tropical Diseases 12:e0006432.https://doi.org/10.1371/journal.pntd.0006432
-
The Monkeybar Project: population density of long-tailed macaques (Macaca fascicularis) in two different forest types in Kudat District, Sabah, MalaysiaAmerican Journal of Physical Anthropology 159:155.
-
ConferenceAssociation of Institutions for Tropical Veterinary Medicine (AITVM)14th International Conference of the Association of Institutions for Tropical Veterinary Medicine.
-
Reemergence of yellow fever in Brazil: The role of distinct landscape fragmentation thresholdsJournal of Environmental and Public Health 2021:8230789.https://doi.org/10.1155/2021/8230789
-
Molecular epidemiological investigation of Plasmodium knowlesi in humans and macaques in SingaporeVector Borne and Zoonotic Diseases 11:131–135.https://doi.org/10.1089/vbz.2010.0024
-
Plasmodium knowlesi infecting humans in Southeast Asia: What’s next?PLOS Neglected Tropical Diseases 14:e0008900.https://doi.org/10.1371/journal.pntd.0008900
-
SoftwareKnowlesiFrag, version swh:1:rev:6eecb0460cd61c26a79413bf91654557eff1380cSoftware Heritage.
-
Zoonotic pathogens survey in free-living long-tailed macaques in ThailandInternational Journal of Veterinary Science and Medicine 10:11–18.https://doi.org/10.1080/23144599.2022.2040176
-
Occupationally Acquired Plasmodium knowlesi Malaria in Brunei DarussalamSafety and Health at Work 10:122–124.https://doi.org/10.1016/j.shaw.2018.09.002
-
Edge effects and their influence on lemur density and distribution in Southeast MadagascarAmerican Journal of Physical Anthropology 129:232–241.https://doi.org/10.1002/ajpa.20241
-
Spatiotemporal distribution of zika virus and its spatially heterogeneous relationship with the environmentInternational Journal of Environmental Research and Public Health 18:290.https://doi.org/10.3390/ijerph18010290
-
Prevalence of simian malaria parasites in macaques of SingaporePLOS Neglected Tropical Diseases 15:e0009110.https://doi.org/10.1371/journal.pntd.0009110
-
Quantifying publication bias in meta-analysisBiometrics 74:785–794.https://doi.org/10.1111/biom.12817
-
Evaluating the links between biodiversity, land-use change, and infectious disease emergenceTropical Conservation pp. 79–88.
-
Investigating the contribution of peri-domestic transmission to risk of zoonotic malaria infection in humansPLOS Neglected Tropical Diseases 10:e0005064.https://doi.org/10.1371/journal.pntd.0005064
-
Accelerated diversification of nonhuman primate malarias in Southeast Asia: Adaptive radiation or geographic speciation?Molecular Biology and Evolution 32:422–439.https://doi.org/10.1093/molbev/msu310
-
Methodological guidance for systematic reviews of observational epidemiological studies reporting prevalence and cumulative incidence dataInternational Journal of Evidence-Based Healthcare 13:147–153.https://doi.org/10.1097/XEB.0000000000000054
-
Ranging patterns and parasitism in primatesBiology Letters 2:351–354.https://doi.org/10.1098/rsbl.2006.0485
-
Spatial and Temporal Analysis of Plasmodium knowlesi Infection in Peninsular Malaysia, 2011 to 2018International Journal of Environmental Research and Public Health 17:9271.https://doi.org/10.3390/ijerph17249271
-
Kasus Baru Plasmodium knowlesi pada Manusia di JambiJurnal Biotek Medisiana Indonesia 6:39–51.https://doi.org/10.22435/JBMI.V6I1.1684
-
Malaria and Hepatocystis Species in Wild Macaques, Southern ThailandThe American Journal of Tropical Medicine and Hygiene 78:646–653.https://doi.org/10.4269/ajtmh.2008.78.646
-
ThesisMolecular prevalence, spatial distribution and epidemiological risk factors of Plasmodium knowlesi infection among Macaca fascicularis raffles in Peninsular MalaysiaUniversiti Putra Malaysia.
-
Estimating Geographical Variation in the Risk of Zoonotic Plasmodium knowlesi Infection in Countries Eliminating MalariaPLOS Neglected Tropical Diseases 10:e0004915.https://doi.org/10.1371/journal.pntd.0004915
-
The heterogeneity statistic I(2) can be biased in small meta-analysesBMC Medical Research Methodology 15:35.https://doi.org/10.1186/s12874-015-0024-z
-
Natural transmission of Plasmodium knowlesi to humans by Anopheles latens in Sarawak, MalaysiaTransactions of the Royal Society of Tropical Medicine and Hygiene 100:1087–1088.https://doi.org/10.1016/j.trstmh.2006.02.006
-
Measuring habitat fragmentation: An evaluation of landscape pattern metricsMethods in Ecology and Evolution 5:634–646.https://doi.org/10.1111/2041-210X.12198
-
Inferring host range dynamics from comparative data: the protozoan parasites of new world monkeysThe American Naturalist 184:65–74.https://doi.org/10.1086/676589
-
Habitat fragmentation, biodiversity loss and the risk of novel infectious disease emergenceJournal of the Royal Society, Interface 15:20180403.https://doi.org/10.1098/rsif.2018.0403
-
Seasonal and Spatial Dynamics of the Primary Vector of Plasmodium knowlesi within a Major Transmission Focus in Sabah, MalaysiaPLOS Neglected Tropical Diseases 9:e0004135.https://doi.org/10.1371/journal.pntd.0004135
Article and author information
Author details
Funding
Wellcome Trust (221963/Z/20/Z)
- Kimberly Fornace
Royal Society (221963/Z/20/Z)
- Kimberly Fornace
Ministry of Higher Education Malaysia (LRGS/1/2018/UM/01/1)
- Reuben Sunil Kumar Sharma
World Health Organization
- Kimberly Fornace
- Chris Drakeley
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. For the purpose of Open Access, the authors have applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.
Acknowledgements
Research was supported by the Sir Henry Dale Fellowship, jointly funded by the Wellcome Trust and the Royal Society (Grant Number 221963/Z/20/Z). Data collected in Peninsular Malaysia by the Faculty of Veterinary Medicine, Universiti Putra Malaysia, was funded by the Ministry of Higher Education Malaysia (Grant LRGS/1/2018/UM/01/1/). Additional funding was supported by the World Health Organization.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Reviewed Preprint version 3:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.88616. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Johnson et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 943
- views
-
- 111
- downloads
-
- 1
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Ecology
- Microbiology and Infectious Disease
Predicting how species diversity changes along environmental gradients is an enduring problem in ecology. In microbes current theories tend to invoke energy availability and enzyme kinetics as the main drivers of temperature-richness relationships. Here we derive a general empirically-grounded theory that can explain this phenomenon by linking microbial species richness in competitive communities to variation in the temperature-dependence of their interaction and growth rates. Specifically, the shape of the microbial community temperature-richness relationship depends on how rapidly the strength of effective competition between species pairs changes with temperature relative to the variance of their growth rates. Furthermore, it predicts that a thermal specialist-generalist tradeoff in growth rates alters coexistence by shifting this balance, causing richness to peak at relatively higher temperatures. Finally, we show that the observed patterns of variation in thermal performance curves of metabolic traits across extant bacterial taxa is indeed sufficient to generate the variety of community-level temperature-richness responses observed in the real world. Our results provide a new and general mechanism that can help explain temperature-diversity gradients in microbial communities, and provide a quantitative framework for interlinking variation in the thermal physiology of microbial species to their community-level diversity.
-
- Ecology
Bumblebees (Bombus terrestris) have been shown to engage in string-pulling behavior to access rewards. The objective of this study was to elucidate whether bumblebees display means-end comprehension in a string-pulling task. We presented bumblebees with two options: one where a string was connected to an artificial flower containing a reward and the other presenting an interrupted string. Bumblebees displayed a consistent preference for pulling connected strings over interrupted ones after training with a stepwise pulling technique. When exposed to novel string colors, bees continued to exhibit a bias towards pulling the connected string. This suggests that bumblebees engage in featural generalization of the visual display of the string connected to the flower in this task. If the view of the string connected to the flower was restricted during the training phase, the proportion of bumblebees choosing the connected strings significantly decreased. Similarly, when the bumblebees were confronted with coiled connected strings during the testing phase, they failed to identify and reject the interrupted strings. This finding underscores the significance of visual consistency in enabling the bumblebees to perform the task successfully. Our results suggest that bumblebees’ ability to distinguish between continuous strings and interrupted strings relies on a combination of image matching and associative learning, rather than means-end understanding. These insights contribute to a deeper understanding of the cognitive processes employed by bumblebees when tackling complex spatial tasks.