skip to main content


Title: Counterintuitive scaling between population abundance and local density: Implications for modelling transmission of infectious diseases in bat populations
Abstract

Models of host–pathogen interactions help to explain infection dynamics in wildlife populations and to predict and mitigate the risk of zoonotic spillover. Insights from models inherently depend on the way contacts between hosts are modelled, and crucially, how transmission scales with animal density.

Bats are important reservoirs of zoonotic disease and are among the most gregarious of all mammals. Their population structures can be highly heterogeneous, underpinned by ecological processes across different scales, complicating assumptions regarding the nature of contacts and transmission. Although models commonly parameterise transmission using metrics of total abundance, whether this is an ecologically representative approximation of host–pathogen interactions is not routinely evaluated.

We collected a 13‐month dataset of tree‐roostingPteropusspp. from 2,522 spatially referenced trees across eight roosts to empirically evaluate the relationship between total roost abundance and tree‐level measures of abundance and density—the scale most likely to be relevant for virus transmission. We also evaluate whether roost features at different scales (roost level, subplot level, tree level) are predictive of these local density dynamics.

Roost‐level features were not representative of tree‐level abundance (bats per tree) or tree‐level density (bats per m2or m3), with roost‐level models explaining minimal variation in tree‐level measures. Total roost abundance itself was either not a significant predictor (tree‐level 3D density) or only weakly predictive (tree‐level abundance).

This indicates that basic measures, such as total abundance of bats in a roost, may not provide adequate approximations for population dynamics at scales relevant for transmission, and that alternative measures are needed to compare transmission potential between roosts. From the best candidate models, the strongest predictor of local population structure was tree density within roosts, where roosts with low tree density had a higher abundance but lower density of bats (more spacing between bats) per tree.

Together, these data highlight unpredictable and counterintuitive relationships between total abundance and local density. More nuanced modelling of transmission, spread and spillover from bats likely requires alternative approaches to integrating contact structure in host–pathogen models, rather than simply modifying the transmission function.

 
more » « less
Award ID(s):
1716698
NSF-PAR ID:
10446964
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Journal of Animal Ecology
Volume:
91
Issue:
5
ISSN:
0021-8790
Page Range / eLocation ID:
p. 916-932
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The spatial organization of populations determines their pathogen dynamics. This is particularly important for communally roosting species, whose aggregations are often driven by the spatial structure of their environment.

    We develop a spatially explicit model for virus transmission within roosts of Australian tree‐dwelling bats (Pteropusspp.), parameterized to reflect Hendra virus. The spatial structure of roosts mirrors three study sites, and viral transmission between groups of bats in trees was modelled as a function of distance between roost trees. Using three levels of tree density to reflect anthropogenic changes in bat habitats, we investigate the potential effects of recent ecological shifts in Australia on the dynamics of zoonotic viruses in reservoir hosts.

    We show that simulated infection dynamics in spatially structured roosts differ from that of mean‐field models for equivalently sized populations, highlighting the importance of spatial structure in disease models of gregarious taxa. Under contrasting scenarios of flying‐fox roosting structures, sparse stand structures (with fewer trees but more bats per tree) generate higher probabilities of successful outbreaks, larger and faster epidemics, and shorter virus extinction times, compared to intermediate and dense stand structures with more trees but fewer bats per tree. These observations are consistent with the greater force of infection generated by structured populations with less numerous but larger infected groups, and may flag an increased risk of pathogen spillover from these increasingly abundant roost types.

    Outputs from our models contribute insights into the spread of viruses in structured animal populations, like communally roosting species, as well as specific insights into Hendra virus infection dynamics and spillover risk in a situation of changing host ecology. These insights will be relevant for modelling other zoonotic viruses in wildlife reservoir hosts in response to habitat modification and changing populations, including coronaviruses like SARS‐CoV‐2.

     
    more » « less
  2. Abstract

    Changes to migration routes and phenology create novel contact patterns among hosts and pathogens. These novel contact patterns can lead to pathogens spilling over between resident and migrant populations. Predicting the consequences of such pathogen spillover events requires understanding how pathogen evolution depends on host movement behaviour. Following spillover, pathogens may evolve changes in their transmission rate and virulence phenotypes because different strategies are favoured by resident and migrant host populations. There is conflict in current theoretical predictions about what those differences might be. Some theory predicts lower pathogen virulence and transmission rates in migrant populations because migrants have lower tolerance to infection. Other theoretical work predicts higher pathogen virulence and transmission rates in migrants because migrants have more contacts with susceptible hosts.

    We aim to understand how differences in tolerance to infection and host pace of life act together to determine the direction of pathogen evolution following pathogen spillover from a resident to a migrant population.

    We constructed a spatially implicit model in which we investigate how pathogen strategy changes following the addition of a migrant population. We investigate how differences in tolerance to infection and pace of life between residents and migrants determine the effect of spillover on pathogen evolution and host population size.

    When the paces of life of the migrant and resident hosts are equal, larger costs of infection in the migrants lead to lower pathogen transmission rate and virulence following spillover. When the tolerance to infection in migrant and resident populations is equal, faster migrant paces of life lead to increased transmission rate and virulence following spillover. However, the opposite can also occur: when the migrant population has lower tolerance to infection, faster migrant paces of life can lead to decreases in transmission rate and virulence.

    Predicting the outcomes of pathogen spillover requires accounting for both differences in tolerance to infection and pace of life between populations. It is also important to consider how movement patterns of populations affect host contact opportunities for pathogens. These results have implications for wildlife conservation, agriculture and human health.

     
    more » « less
  3. Abstract

    Pathogens can spill over and infect new host species by overcoming a series of ecological and biological barriers. Hendra virus (HeV) circulates in Australian flying foxes and provides a data‐rich study system for identifying environmental drivers underlying spillover events. The frequency of spillover events to horses has varied interannually since the virus was first discovered in 1994. These observations suggest that HeV spillover events are driven, in part, by environmental factors, including loss of flying fox habitat and climate variability.

    We explicitly examine the impact of environmental variation on the risk of HeV spillover at three spatial scales relevant to this system. We use a dataset of 60 spillover events and boosted regression tree methods to identify environmental features (including concurrent and lagged temperature, rainfall, vegetation indices, land cover, and climate indices) at three spatial scales (1‐km, 20‐km, 100‐km radii) associated with horse contacts and reservoir species ecology.

    We find that temperature, local (1‐km radius) human population density, and landscape (100‐km radius) forest cover and pasture are the most influential environmental features associated with HeV spillover risk. By including multiple spatial scales and temporal lags in environmental features, we can more accurately quantify risk across space and time than with models that use a single scale. For example, high quality vegetation at the local scale and within a foraging radius (20‐km) in the concurrent month and previous years, combined with poorer quality vegetation at the landscape scale in the concurrent month increase risk of HeV spillover. These and other environmental associations likely influence the dynamic foraging behaviour of reservoir flying foxes and drive contacts that facilitate spillover into horse populations.

    Synthesis and application: Current management of HeV spillover focuses on local‐scale interventions – primarily through vaccination and detection of infected horses. Our study finds that HeV spillover risk is also driven by environmental changes over much larger scales and demonstrates management practices would benefit from incorporating landscape interventions alongside local interventions.

     
    more » « less
  4. Abstract

    Host ecological factors and external environmental factors are known to influence the structure of gut microbial communities, but few studies have examined the impacts of environmental changes on microbiotas in free‐ranging animals. Rapid land‐use change has the potential to shift gut microbial communities in wildlife through exposure to novel bacteria and/or by changing the availability or quality of local food resources. The consequences of such changes to host health and fitness remain unknown and may have important implications for pathogen spillover between humans and wildlife. To better understand the consequences of land‐use change on wildlife microbiotas, we analyzed long‐term dietary trends, gut microbiota composition, and innate immune function in common vampire bats (Desmodus rotundus) in two nearby sites in Belize that vary in landscape structure. We found that vampire bats living in a small forest fragment had more homogenous diets indicative of feeding on livestock and shifts in microbiota heterogeneity, but not overall composition, compared to those living in an intact forest reserve. We also found that irrespective of sampling site, vampire bats which consumed relatively more livestock showed shifts in some core bacteria compared with vampire bats which consumed relatively less livestock. The relative abundance of some core microbiota members was associated with innate immune function, suggesting that future research should consider the role of the host microbiota in immune defense and its relationship to zoonotic infection dynamics. We suggest that subsequent homogenization of diet and habitat loss through livestock rearing in the Neotropics may lead to disruption to the microbiota that could have downstream impacts on host immunity and cross‐species pathogen transmission.

     
    more » « less
  5. Abstract

    The exodus of flying animals from their roosting locations is often visible as expanding ring‐shaped patterns in weather radar data. The NEXRAD network, for example, archives more than 25 years of data across 143 contiguous US radar stations, providing opportunities to study roosting locations and times and the ecosystems of birds and bats. However, access to this information is limited by the cost of manually annotating millions of radar scans. We develop and deploy an AI‐assisted system to annotate roosts in radar data. We build datasets with roost annotations to support the training and evaluation of automated detection models. Roosts are detected, tracked, and incorporated into our developed web‐based interface for human screening to produce research‐grade annotations. We deploy the system to collect swallow and martin roost information from 12 radar stations around the Great Lakes spanning 21 years. After verifying the practical value of the system, we propose to improve the detector by incorporating both spatial and temporal channels from volumetric radar scans. The deployment on Great Lakes radar scans allows accelerated annotation of 15 628 roost signatures in 612 786 radar scans with 183.6 human screening hours, or 1.08 s per radar scan. We estimate that the deployed system reduces human annotation time by ~7×. The temporal detector model improves the average precision at intersection‐over‐union threshold 0.5 (APIoU = .50) by 8% over the previous model (48%→56%), further reducing human screening time by 2.3× in its pilot deployment. These data contain critical information about phenology and population trends of swallows and martins, aerial insectivore species experiencing acute declines, and have enabled novel research. We present error analyses, lay the groundwork for continent‐scale historical investigation about these species, and provide a starting point for automating the detection of other family‐specific phenomena in radar data, such as bat roosts and mayfly hatches.

     
    more » « less