Species distribution modeling for disease ecology: A multi-scale case study for schistosomiasis host snails in Brazil.

Species distribution models (SDMs) are increasingly popular tools for profiling disease risk in ecology, particularly for infectious diseases of public health importance that include an obligate non-human host in their transmission cycle. SDMs can create high-resolution maps of host distribution across geographical scales, reflecting baseline risk of disease. However, as SDM computational methods have rapidly expanded, there are many outstanding methodological questions. Here we address key questions about SDM application, using schistosomiasis risk in Brazil as a case study. Schistosomiasis is transmitted to humans through contact with the free-living infectious stage of Schistosoma spp. parasites released from freshwater snails, the parasite's obligate intermediate hosts. In this study, we compared snail SDM performance across machine learning (ML) approaches (MaxEnt, Random Forest, and Boosted Regression Trees), geographic extents (national, regional, and state), types of presence data (expert-collected and publicly-available), and snail species (Biomphalaria glabrata, B. straminea, and B. tenagophila). We used high-resolution (1km) climate, hydrology, land-use/land-cover (LULC), and soil property data to describe the snails' ecological niche and evaluated models on multiple criteria. Although all ML approaches produced comparable spatially cross-validated performance metrics, their suitability maps showed major qualitative differences that required validation based on local expert knowledge. Additionally, our findings revealed varying importance of LULC and bioclimatic variables for different snail species at different spatial scales. Finally, we found that models using publicly-available data predicted snail distribution with comparable AUC values to models using expert-collected data. This work serves as an instructional guide to SDM methods that can be applied to a range of vector-borne and zoonotic diseases. In addition, it advances our understanding of the relevant environment and bioclimatic determinants of schistosomiasis risk in Brazil.

A flexible model for thermal performance curves.

Temperature responses of many biological traits-including population growth, survival, and development-are described by thermal performance curves (TPCs) with phenomenological models like the Briere function or mechanistic models related to chemical kinetics. Existing TPC models are either simple but inflexible in shape, or flexible yet difficult to interpret in biological terms. Here we present flexTPC: a model that is parameterized exclusively in terms of biologically interpretable quantities, including the thermal minimum, optimum, and maximum, and the maximum trait value. FlexTPC can describe unimodal temperature responses of any skewness and thermal breadth, enabling direct comparisons across populations, traits, or taxa with a single model. We apply flexTPC to various microbial and entomological datasets, compare results with the Briere model, and find that flexTPC often has better predictive performance. The interpretability of flexTPC makes it ideal for modeling how thermal responses change with ecological stressors or evolve over time.

Plasticity in mosquito size and thermal tolerance across a latitudinal climate gradient.

Variation in heat tolerance among populations can determine whether a species is able to cope with ongoing climate change. Such variation may be especially important for ectotherms whose body temperatures, and consequently, physiological processes, are regulated by external conditions. Additionally, differences in body size are often associated with latitudinal clines, thought to be driven by climate gradients. While studies have begun to explore variation in body size and heat tolerance within species, our understanding of these patterns across large spatial scales, particularly regarding the roles of plasticity and genetic differences, remains incomplete. Here, we examine body size, as measured by wing length, and thermal tolerance, as measured by the time to immobilisation at high temperatures ("thermal knockdown"), in populations of the mosquito Aedes sierrensis collected from across a large latitudinal climate gradient spanning 1300 km (34-44° N). We find that mosquitoes collected from lower latitudes and warmer climates were more tolerant of high temperatures than those collected from higher latitudes and colder climates. Moreover, body size increased with latitude and decreased with temperature, a pattern consistent with James' rule, which appears to be a result of plasticity rather than genetic variation. Our results suggest that warmer environments produce smaller and more thermally tolerant populations.

A Mosquito Parasite Is Locally Adapted to Its Host but Not Temperature.

AbstractClimate change will alter interactions between parasites and their hosts. Warming may affect patterns of local adaptation, shifting the environment to favor the parasite or host and thus changing the prevalence of disease. We assessed local adaptation to hosts and temperature in the facultative ciliate parasite Lambornella clarki, which infects the western tree hole mosquito Aedes sierrensis. We conducted laboratory infection experiments with mosquito larvae and parasites collected from across a climate gradient, pairing sympatric or allopatric populations across three temperatures that were either matched or mismatched to the source environment. Lambornella clarki parasites were locally adapted to their hosts, with 2.6 times higher infection rates on sympatric populations compared with allopatric populations, but they were not locally adapted to temperature. Infection peaked at the intermediate temperature of 12.5°C, notably lower than the optimum temperature for free-living L. clarki growth, suggesting that the host's immune response can play a significant role in mediating the outcome of infection. Our results highlight the importance of host selective pressure on parasites, despite the impact of temperature on infection success.

Climate and urbanization drive changes in the habitat suitability of Schistosoma mansoni competent snails in Brazil.

Schistosomiasis is a neglected tropical disease caused by Schistosoma parasites. Schistosoma are obligate parasites of freshwater Biomphalaria and Bulinus snails, thus controlling snail populations is critical to reducing transmission risk. As snails are sensitive to environmental conditions, we expect their distribution is significantly impacted by global change. Here, we used machine learning, remote sensing, and 30 years of snail occurrence records to map the historical and current distribution of forward-transmitting Biomphalaria hosts throughout Brazil. We identified key features influencing the distribution of suitable habitat and determined how Biomphalaria habitat has changed with climate and urbanization over the last three decades. Our models show that climate change has driven broad shifts in snail host range, whereas expansion of urban and peri-urban areas has driven localized increases in habitat suitability. Elucidating change in Biomphalaria distribution-while accounting for non-linearities that are difficult to detect from local case studies-can help inform schistosomiasis control strategies.

Re-assessing thermal response of schistosomiasis transmission risk: Evidence for a higher thermal optimum than previously predicted.

The geographical range of schistosomiasis is affected by the ecology of schistosome parasites and their obligate host snails, including their response to temperature. Previous models predicted schistosomiasis' thermal optimum at 21.7°C, which is not compatible with the temperature in sub-Saharan Africa (SSA) regions where schistosomiasis is hyperendemic. We performed an extensive literature search for empirical data on the effect of temperature on physiological and epidemiological parameters regulating the free-living stages of S. mansoni and S. haematobium and their obligate host snails, i.e., Biomphalaria spp. and Bulinus spp., respectively. We derived nonlinear thermal responses fitted on these data to parameterize a mechanistic, process-based model of schistosomiasis. We then re-cast the basic reproduction number and the prevalence of schistosome infection as functions of temperature. We found that the thermal optima for transmission of S. mansoni and S. haematobium range between 23.1-27.3°C and 23.6-27.9°C (95% CI) respectively. We also found that the thermal optimum shifts toward higher temperatures as the human water contact rate increases with temperature. Our findings align with an extensive dataset of schistosomiasis prevalence in SSA. The refined nonlinear thermal-response model developed here suggests a more suitable current climate and a greater risk of increased transmission with future warming for more than half of the schistosomiasis suitable regions with mean annual temperature below the thermal optimum.

Re-assessing thermal response of schistosomiasis transmission risk: evidence for a higher thermal optimum than previously predicted.

The geographical range of schistosomiasis is affected by the ecology of schistosome parasites and their obligate host snails, including their response to temperature. Previous models predicted schistosomiasis' thermal optimum at 21.7 °C, which is not compatible with the temperature in sub-Saharan Africa (SSA) regions where schistosomiasis is hyperendemic. We performed an extensive literature search for empirical data on the effect of temperature on physiological and epidemiological parameters regulating the free-living stages of S. mansoni and S. haematobium and their obligate host snails, i.e., Biomphalaria spp. and Bulinus spp., respectively. We derived nonlinear thermal responses fitted on these data to parameterize a mechanistic, process-based model of schistosomiasis. We then re-cast the basic reproduction number and the prevalence of schistosome infection as functions of temperature. We found that the thermal optima for transmission of S. mansoni and S. haematobium range between 23.1-27.3 °C and 23.6-27.9 °C (95 % CI) respectively. We also found that the thermal optimum shifts toward higher temperatures as the human water contact rate increases with temperature. Our findings align with an extensive dataset of schistosomiasis prevalence in SSA. The refined nonlinear thermal-response model developed here suggests a more suitable current climate and a greater risk of increased transmission with future warming for more than half of the schistosomiasis suitable regions with mean annual temperature below the thermal optimum.

Population-specific thermal responses contribute to regional variability in arbovirus transmission with changing climates.

Temperature is increasing globally, and vector-borne diseases are particularly responsive to such increases. While it is known that temperature influences mosquito life history traits, transmission models have not historically considered population-specific effects of temperature. We assessed the interaction between Culex pipiens population and temperature in New York State (NYS) and utilized novel empirical data to inform predictive models of West Nile virus (WNV) transmission. Genetically and regionally distinct populations from NYS were reared at various temperatures, and life history traits were monitored and used to inform trait-based models. Variation in Cx. pipiens life history traits and population-dependent thermal responses account for a predicted 2.9°C difference in peak transmission that is reflected in regional differences in WNV prevalence. We additionally identified genetic signatures that may contribute to distinct thermal responses. Together, these data demonstrate how population variation contributes to significant geographic variability in arbovirus transmission with changing climates.

Interconnecting global threats: climate change, biodiversity loss, and infectious diseases.

The concurrent pressures of rising global temperatures, rates and incidence of species decline, and emergence of infectious diseases represent an unprecedented planetary crisis. Intergovernmental reports have drawn focus to the escalating climate and biodiversity crises and the connections between them, but interactions among all three pressures have been largely overlooked. Non-linearities and dampening and reinforcing interactions among pressures make considering interconnections essential to anticipating planetary challenges. In this Review, we define and exemplify the causal pathways that link the three global pressures of climate change, biodiversity loss, and infectious disease. A literature assessment and case studies show that the mechanisms between certain pairs of pressures are better understood than others and that the full triad of interactions is rarely considered. Although challenges to evaluating these interactions-including a mismatch in scales, data availability, and methods-are substantial, current approaches would benefit from expanding scientific cultures to embrace interdisciplinarity and from integrating animal, human, and environmental perspectives. Considering the full suite of connections would be transformative for planetary health by identifying potential for co-benefits and mutually beneficial scenarios, and highlighting where a narrow focus on solutions to one pressure might aggravate another.

Local tree cover predicts mosquito species richness and disease vector presence in a tropical countryside landscape.

Context: Land use change drives both biodiversity loss and zoonotic disease transmission in tropical countryside landscapes. Developing solutions for protecting countryside biodiversity, public health, and livelihoods requires understanding the scales at which habitat characteristics such as land cover shape biodiversity, especially for arthropods that transmit pathogens. Evidence increasingly shows that species richness for many taxa correlates with local tree cover. Objectives: We investigated whether mosquito species richness, community composition, and presence of disease vector species responded to land use and tree cover - and if so, whether at spatial scales similar to other taxa. Methods: We paired a field survey of mosquito communities in agricultural, residential, and forested lands in rural southern Costa Rica with remotely sensed tree cover data. We compared mosquito community responses to tree cover surrounding survey sites measured across scales, and analyzed community responses to land use and environmental gradients. Results: Tree cover was positively correlated with mosquito species richness, and negatively correlated with the presence of the common invasive dengue vector Aedes albopictus, particularly at small spatial scales of 80 - 200m. Land use predicted community composition and Ae. albopictus presence. Environmental gradients of tree cover, temperature, and elevation explained 7% of species turnover among survey sites. Conclusions: The results suggest that preservation and expansion of tree cover at local scales can protect biodiversity for a wide range of taxa, including arthropods, and also confer protection against disease vector occurrence. The identified spatial range of tree cover benefits can inform land management for conservation and public health protection.

Perceived experts are prevalent and influential within an antivaccine community on Twitter.

Perceived experts (i.e. medical professionals and biomedical scientists) are trusted sources of medical information who are especially effective at encouraging vaccine uptake. The role of perceived experts acting as potential antivaccine influencers has not been characterized systematically. We describe the prevalence and importance of antivaccine perceived experts by constructing a coengagement network of 7,720 accounts based on a Twitter data set containing over 4.2 million posts from April 2021. The coengagement network primarily broke into two large communities that differed in their stance toward COVID-19 vaccines, and misinformation was predominantly shared by the antivaccine community. Perceived experts had a sizable presence across the coengagement network, including within the antivaccine community where they were 9.8% of individual, English-language users. Perceived experts within the antivaccine community shared low-quality (misinformation) sources at similar rates and academic sources at higher rates compared to perceived nonexperts in that community. Perceived experts occupied important network positions as central antivaccine users and bridges between the antivaccine and provaccine communities. Using propensity score matching, we found that perceived expertise brought an influence boost, as perceived experts were significantly more likely to receive likes and retweets in both the antivaccine and provaccine communities. There was no significant difference in the magnitude of the influence boost for perceived experts between the two communities. Social media platforms, scientific communications, and biomedical organizations may focus on more systemic interventions to reduce the impact of perceived experts in spreading antivaccine misinformation.

Non-household environments make a major contribution to dengue transmission: Implications for vector control.

Aedes-borne pathogens have been increasing in incidence in recent decades despite vector control activities implemented in endemic settings. Vector control for Aedes-transmitted arboviruses typically focuses on households because vectors breed in household containers and bite indoors. Yet, our recent work shows a high abundance of Aedes spp. vectors in public spaces. To investigate the impact of non-household environments on dengue transmission and control, we used field data on the number of water containers and abundance of Aedes mosquitoes in Household (HH) and Non-Household (NH) environments in two Kenyan cities, Kisumu and Ukunda, from 2019-2022. Incorporating information on human activity space, we developed an agent-based model to simulate city-wide conditions considering HH and five types of NH environments in which people move and interact with other humans and vectors during peak biting times. We additionally evaluated the outcome of vector control activities implemented in different environments in preventive (before an epidemic) and reactive (after an epidemic commences) scenarios. We estimated that over half of infections take place in NH environments, where the main spaces for transmission are workplaces, markets, and recreational locations. Accordingly, results highlight the important role of vector control activities at NH locations to reduce dengue. A greater reduction of cases is expected as control activities are implemented earlier, at higher levels of coverage, with greater effectiveness when targeting only NH as opposed to when targeting only HH. Further, local ecological factors such as the differential abundance of water containers within cities are also influential factors to consider for control. This work provides insight into the importance of vector control in both household and non-household environments in endemic settings. It highlights a specific approach to inform evidence-based decision making to target limited vector control resources for optimal control.

Climate warming is expanding dengue burden in the Americas and Asia.

Climate change poses significant threats to public health, with dengue representing a growing concern due to its high existing burden and sensitivity to climatic conditions. Yet, the quantitative impacts of temperature warming on dengue, both in the past and in the future, remain poorly understood. In this study, we quantify how dengue responds to climatic fluctuations, and use this inferred temperature response to estimate the impacts of historical warming and forecast trends under future climate change scenarios. To estimate the causal impact of temperature on the spread of dengue in the Americas and Asia, we assembled a dataset encompassing nearly 1.5 million dengue incidence records from 21 countries. Our analysis revealed a nonlinear relationship between temperature and dengue incidence with the largest marginal effects at lower temperatures (around 15°C), peak incidence at 27.8°C (95% CI: 27.3 - 28.2°C), and subsequent declines at higher temperatures. Our findings indicate that historical climate change has already increased dengue incidence 18% (12 - 25%) in the study region, and projections suggest a potential increase of 40% (17 - 76) to 57% (33 - 107%) by mid-century depending on the climate scenario, with some areas seeing up to 200% increases. Notably, our models suggest that lower emissions scenarios would substantially reduce the warming-driven increase in dengue burden. Together, these insights contribute to the broader understanding of how long-term climate patterns influence dengue, providing a valuable foundation for public health planning and the development of strategies to mitigate future risks due to climate change.

Mosquito thermal tolerance is remarkably constrained across a large climatic range.

How mosquitoes may respond to rapid climate warming remains unknown for most species, but will have major consequences for their future distributions, with cascading impacts on human well-being, biodiversity and ecosystem function. We investigated the adaptive potential of a wide-ranging mosquito species, Aedes sierrensis, across a large climatic gradient by conducting a common garden experiment measuring the thermal limits of mosquito life-history traits. Although field-collected populations originated from vastly different thermal environments that spanned over 1200 km, we found limited variation in upper thermal tolerance between populations. In particular, the upper thermal limits of all life-history traits varied by less than 3°C across the species range and, for most traits, did not differ significantly between populations. For one life-history trait-pupal development rate-we did detect significant variation in upper thermal limits between populations, and this variation was strongly correlated with source temperatures, providing evidence of local thermal adaptation for pupal development. However, we found that maximum environmental temperatures across most of the species' range already regularly exceed the highest upper thermal limits estimated under constant temperatures. This result suggests that strategies for coping with and/or avoiding thermal extremes are likely key components of current and future mosquito thermal tolerance.

Climate and urbanization drive changes in the habitat suitability of Schistosoma mansoni competent snails in Brazil.

Schistosomiasis is a neglected tropical disease caused by Schistosoma parasites. Schistosoma are obligate parasites of freshwater Biomphalaria snails, so controlling snail populations is critical to reducing transmission risk. As snails are sensitive to environmental conditions, we expect their distribution is significantly impacted by global change. Here, we leveraged machine learning, remote sensing, and 30 years of snail occurrence records to map the historical and current distribution of competent Biomphalaria throughout Brazil. We identified key features influencing the distribution of suitable habitat and determined how Biomphalaria habitat has changed with climate and urbanization over the last three decades. Our models show that climate change has driven broad shifts in snail host range, whereas expansion of urban and peri-urban areas has driven localized increases in habitat suitability. Elucidating change in Biomphalaria distribution - while accounting for non-linearities that are difficult to detect from local case studies - can help inform schistosomiasis control strategies.

Local tree cover predicts mosquito species richness and disease vector presence in a tropical countryside landscape.

Context: Land use change drives both biodiversity loss and zoonotic disease transmission in tropical countryside landscapes. Developing solutions for protecting countryside biodiversity, public health, and livelihoods requires understanding the scales at which habitat characteristics such as land cover shape biodiversity, especially for arthropods that transmit pathogens. Evidence increasingly shows that species richness for many taxa correlates with local tree cover. Objectives: We investigated whether mosquito species richness, community composition, and presence of disease vector species responded to land use and tree cover - and if so, whether at spatial scales similar to other taxa. Methods: We paired a field survey of mosquito communities in agricultural, residential, and forested lands in rural southern Costa Rica with remotely sensed tree cover data. We compared mosquito community responses to tree cover surrounding survey sites measured across scales, and analyzed community responses to land use and environmental gradients. Results: Tree cover was positively correlated with mosquito species richness, and negatively correlated with the presence of the common invasive dengue vector Aedes albopictus , particularly at small spatial scales of 80 - 200m. Land use predicted community composition and Ae. albopictus presence. Environmental gradients of tree cover, temperature, and elevation explained 7% of species turnover among survey sites. Conclusions: The results suggest that preservation and expansion of tree cover at local scales can protect biodiversity for a wide range of taxa, including arthropods, and also confer protection against disease vector occurrence. The identified spatial range of tree cover benefits can inform land management for conservation and public health protection.

Global malaria predictors at a localized scale.

Malaria is a life-threatening disease caused by Plasmodium parasites transmitted by Anopheles mosquitoes. In 2021, more than 247 million cases of malaria were reported worldwide, with an estimated 619,000 deaths. While malaria incidence has decreased globally in recent decades, some public health gains have plateaued, and many endemic hotspots still face high transmission rates. Understanding local drivers of malaria transmission is crucial but challenging due to the complex interactions between climate, entomological and human variables, and land use. This study focuses on highly climatically suitable and endemic areas in Côte d'Ivoire to assess the explanatory power of coarse climatic predictors of malaria transmission at a fine scale. Using data from 40 villages participating in a randomized controlled trial of a household malaria intervention, the study examines the effects of climate variation over time on malaria transmission. Through panel regressions and statistical modeling, the study investigates which variable (temperature, precipitation, or entomological inoculation rate) and its form (linear or unimodal) best explains seasonal malaria transmission and the factors predicting spatial variation in transmission. The results highlight the importance of temperature and rainfall, with quadratic temperature and all precipitation models performing well, but the causal influence of each driver remains unclear due to their strong correlation. Further, an independent, mechanistic temperature-dependent R 0 model based on laboratory data aligns well with observed malaria incidence rates, emphasizing the significance and predictability of temperature suitability across scales. By contrast, entomological variables, such as entomological inoculation rate, were not strong predictors of human incidence in this context. Finally, the study explores the predictors of spatial variation in malaria, considering land use, intervention, and entomological variables. The findings contribute to a better understanding of malaria transmission dynamics at local scales, aiding in the development of effective control strategies in endemic regions.

Temperature dependence of mosquitoes: comparing mechanistic and machine learning approaches.

Mosquito vectors of pathogens (e.g., Aedes , Anopheles , and Culex spp. which transmit dengue, Zika, chikungunya, West Nile, malaria, and others) are of increasing concern for global public health. These vectors are geographically shifting under climate and other anthropogenic changes. As small-bodied ectotherms, mosquitoes are strongly affected by temperature, which causes unimodal responses in mosquito life history traits (e.g., biting rate, adult mortality rate, mosquito development rate, and probability of egg-to-adult survival) that exhibit upper and lower thermal limits and intermediate thermal optima in laboratory studies. However, it remains unknown how mosquito thermal responses measured in laboratory experiments relate to the realized thermal responses of mosquitoes in the field. To address this gap, we leverage thousands of global mosquito occurrences and geospatial satellite data at high spatial resolution to construct machine-learning based species distribution models, from which vector thermal responses are estimated. We apply methods to restrict models to the relevant mosquito activity season and to conduct ecologically-plausible spatial background sampling centered around ecoregions for comparison to mosquito occurrence records. We found that thermal minima estimated from laboratory studies were highly correlated with those from the species distributions (r = 0.90). The thermal optima were less strongly correlated (r = 0.69). For most species, we did not detect thermal maxima from their observed distributions so were unable to compare to laboratory-based estimates. The results suggest that laboratory studies have the potential to be highly transportable to predicting lower thermal limits and thermal optima of mosquitoes in the field. At the same time, lab-based models likely capture physiological limits on mosquito persistence at high temperatures that are not apparent from field-based observational studies but may critically determine mosquito responses to climate warming.

Mosquito thermal tolerance is remarkably constrained across a large climatic range.

How mosquitoes may respond to rapid climate warming remains unknown for most species, but will have major consequences for their future distributions, with cascading impacts on human well-being, biodiversity, and ecosystem function. We investigated the adaptive potential of a wide-ranging mosquito species, Aedes sierrensis, across a large climatic gradient by conducting a common garden experiment measuring the thermal limits of mosquito life history traits. Although field-collected populations originated from vastly different thermal environments that spanned over 1,200 km, we found remarkably limited variation in upper thermal tolerance between populations, with the upper thermal limits of fitness varying by <1°C across the species range. For one life history trait-pupal development rate-we did detect significant variation in upper thermal limits between populations, and this variation was strongly correlated with source temperatures, providing evidence of local thermal adaptation for pupal development. However, we found environmental temperatures already regularly exceed our highest estimated upper thermal limits throughout most of the species range, suggesting limited potential for mosquito thermal tolerance to evolve on pace with warming. Strategies for avoiding high temperatures such as diapause, phenological shifts, and behavioral thermoregulation are likely important for mosquito persistence.