The Fine-Scale Landscape of Immunity and Parasitism in a Wild Ungulate Population

Fine-scale spatial variation in fitness, inbreeding, and inbreeding depression in a wild ungulate

Environmental stress can exacerbate inbreeding depression by amplifying differences between inbred and outbred individuals. In wild populations, where the environment is often unpredictable and stress can be highly detrimental, the interplay between inbreeding depression and environmental variation is potentially important. Here, we investigate variation in inbreeding level, fitness and strength of inbreeding depression across a fine-scale geographic area (~12 km2) in an individually monitored population of red deer (Cervus elaphus). We show that northern regions of the study area have lower birth weights, lower juvenile survival rates, and higher inbreeding coefficients. Such fine-scale differences in inbreeding coefficients could be caused by the mating system of red deer combined with female density variation. We then tested for an inbreeding depression-by-environment interaction (ID × E) in birth weight and juvenile survival, by fitting an interaction term between the inbreeding coefficient and geographic location. We find that inbreeding depression in juvenile survival is stronger in the harsher northern regions, indicating the presence of ID × E. We also highlight that the ability to infer ID × E is affected by the variation in inbreeding within each geographic region. Therefore, for future studies on ID × E in wild populations, we recommend first assessing whether inbreeding and traits vary spatially or temporally. Overall, this is one of only a handful of studies to find evidence for ID × E in a wild population-despite its prevalence in experimental systems-likely due to intense data demands or insufficient variation in environmental stress or inbreeding coefficients.

Divergent age-related changes in parasite infection occur independently of behaviour and demography in a wild ungulate

As animals age, they exhibit a suite of phenotypic changes, often including reductions in movement and social behaviour ('behavioural ageing'). By altering an individual's exposure to parasites, behavioural ageing may influence infection status trajectories over the lifespan. However, these processes could be confounded by age-related changes in other phenotypic traits, or by selective disappearance of certain individuals owing to parasite-induced mortality. Here, we uncover contrasting age-related patterns of infection across three helminth parasites in wild adult female red deer (Cervus elaphus). Counts of strongyle nematodes (order: Strongylida) increased with age, while counts of liver fluke (Fasciola hepatica) and tissue worm (Elaphostrongylus cervi) decreased, and lungworm (Dictyocaulus) counts did not change. These relationships could not be explained by socio-spatial behaviours, spatial structuring, or selective disappearance, suggesting behavioural ageing is unlikely to be responsible for driving age trends. Instead, social connectedness and strongyle infection were positively correlated, such that direct age-infection trends were directly contrasted with the effects implied by previously documented behavioural ageing. This suggests that behavioural ageing may reduce parasite exposure, potentially countering other age-related changes. These findings demonstrate that different parasites can show contrasting age trajectories depending on diverse intrinsic and extrinsic factors, and that behaviour's role in these processes is likely to be complex and multidirectional.This article is part of the discussion meeting issue 'Understanding age and society using natural populations'.

Expanding theory, methodology and empirical systems at the spatial-social interface

All animals exhibit some combination of spatial and social behaviours. A diversity of interactions occurs between such behaviours, producing emergent phenomena at the spatial-social interface. Untangling and interrogating these complex, intertwined processes can be vital for identifying the mechanisms, causes and consequences of behavioural variation in animal ecology. Nevertheless, the integrated study of the interactions between spatial and social phenotypes and environments (at the spatial-social interface) is in its relative infancy. In this theme issue, we present a collection of papers chosen to expand the spatial-social interface along several theoretical, methodological and empirical dimensions. They detail new perspectives, methods, study systems and more, as well as offering roadmaps for applied outputs and detailing exciting new directions for the field to move in the future. In this Introduction, we outline the contents of these papers, placing them in the context of what comes before, and we synthesize a number of takeaways and future directions for the spatial-social interface. This article is part of the theme issue 'The spatial-social interface: a theoretical and empirical integration'.

Understanding spatiotemporal effects of food supplementation on host-parasite interactions using community-based science

Supplemental feeding can increase the overall health of animals but also can have variable effects on how animals defend themselves against parasites. However, the spatiotemporal effects of food supplementation on host-parasite interactions remain poorly understood, likely because large-scale, coordinated efforts to investigate them are difficult. Here, we introduce the Nest Parasite Community Science Project, which is a community-based science project that coordinates studies with bird nest box 'stewards' from the public and scientific community. This project was established to understand broad ecological patterns between hosts and their parasites. The goal of this study was to determine the effect of food supplementation on eastern bluebirds (Sialia sialis) and their nest parasite community across the geographic range of the bluebirds from 2018 to 2021. We received 674 nests from 69 stewards in 26 states in the eastern United States. Nest box stewards reported whether or not they provided mealworms or suet near nesting bluebirds, then they followed the nesting success of the birds (number of eggs laid and hatched, proportion that hatched, number and proportion of nestlings that successfully fledged). We then identified and quantified parasites in the nests. Overall, we found that food supplementation increased fledging success. The most common nest parasite taxon was the parasitic blow fly (Protocalliphora sialia), but a few nests contained fleas (Ceratophyllus idius, C. gallinae and Orchopeas leucopus) and mites (Dermanyssus spp. and Ornithonyssus spp.). Blow flies were primarily found at northern latitudes, where food supplementation affected blow fly prevalence. However, the direction of this effect varied substantially in direction and magnitude across years. More stewards fed bluebirds at southern latitudes than at northern latitudes, which contradicts the findings of other community-based science projects. Overall, food supplementation of birds was associated with increased host fitness but did not appear to play a consistent role in defence against these parasites across all years. Our study demonstrates the importance of coordinated studies across years and locations to understand the effects of environmental heterogeneity, including human-based food supplementation, on host-parasite dynamics.

Spatiotemporal-social association predicts immunological similarity in rewilded mice

Environmental influences on immune phenotypes are well-documented, but our understanding of which elements of the environment affect immune systems, and how, remains vague. Behaviors, including socializing with others, are central to an individual's interaction with its environment. We therefore tracked behavior of rewilded laboratory mice of three inbred strains in outdoor enclosures and examined contributions of behavior, including associations measured from spatiotemporal co-occurrences, to immune phenotypes. We found extensive variation in individual and social behavior among and within mouse strains upon rewilding. In addition, we found that the more associated two individuals were, the more similar their immune phenotypes were. Spatiotemporal association was particularly predictive of similar memory T and B cell profiles and was more influential than sibling relationships or shared infection status. These results highlight the importance of shared spatiotemporal activity patterns and/or social networks for immune phenotype and suggest potential immunological correlates of social life.

Age-specific impacts of vegetation functional traits on gastrointestinal nematode parasite burdens in a large herbivore

Gastrointestinal nematode (GIN) parasites play an important role in the ecological dynamics of many animal populations. Recent studies suggest that fine-scale spatial variation in GIN infection dynamics is important in wildlife systems, but the environmental drivers underlying this variation remain poorly understood. We used data from over two decades of GIN parasite egg counts, host space use, and spatial vegetation data from a long-term study of Soay sheep on St Kilda to test how spatial autocorrelation and vegetation in an individual's home range predict parasite burden across three age groups. We developed a novel approach to quantify the plant functional traits present in a home range to describe the quality of vegetation present. Effects of vegetation and space varied between age classes. In immature lambs, strongyle parasite faecal egg counts (FEC) were spatially structured, being highest in the north and south of our study area. Independent of host body weight and spatial autocorrelation, plant functional traits predicted parasite egg counts. Higher egg counts were associated with more digestible and preferred plant functional traits, suggesting the association could be driven by host density and habitat preference. In contrast, we found no evidence that parasite FEC were related to plant functional traits in the host home range in yearlings or adult sheep. Adult FEC were spatially structured, with highest burdens in the north-east of our study area, while yearling FEC showed no evidence of spatial structuring. Parasite burdens in immature individuals appear more readily influenced by fine-scale spatial variation in the environment, highlighting the importance of such heterogeneity for our understanding of wildlife epidemiology and health. Our findings support the importance of fine-scale environmental variation for wildlife disease ecology and provides new evidence that such effects may vary across demographic groups within a population.

Spatial Variation in the Inflammatory Response of House Sparrows in their Native Range

Characterizing spatial differences in wildlife immunity is the first step to identify environmental drivers of host defense and disease risks. The house sparrow (Passer domesticus) is a model system for ecoimmunology, but spatial differences in immunity have been largely restricted to the invasive range of this global species. We provide an initial test of spatial variation in immune response to phytohemagglutinin in the native range, finding that birds from Romania have greater inflammatory responses than birds from Egypt. Future broad surveys across the house sparrow native range could contextualize these differences and determine underlying drivers.

A mixed-model approach for estimating drivers of microbiota community composition and differential taxonomic abundance

Next-generation sequencing (NGS) and metabarcoding approaches are increasingly applied to wild animal populations, but there is a disconnect between the widely applied generalized linear mixed model (GLMM) approaches commonly used to study phenotypic variation and the statistical toolkit from community ecology typically applied to metabarcoding data. Here, we describe the suitability of a novel GLMM-based approach for analyzing the taxon-specific sequence read counts derived from standard metabarcoding data. This approach allows decomposition of the contribution of different drivers to variation in community composition (e.g., age, season, individual) via interaction terms in the model random-effects structure. We provide guidance to implementing this approach and show how these models can identify how responsible specific taxonomic groups are for the effects attributed to different drivers. We applied this approach to two cross-sectional data sets from the Soay sheep population of St. Kilda. GLMMs showed agreement with dissimilarity-based approaches highlighting the substantial contribution of age and minimal contribution of season to microbiota community compositions, and simultaneously estimated the contribution of other technical and biological factors. We further used model predictions to show that age effects were principally due to increases in taxa of the phylum Bacteroidetes and declines in taxa of the phylum Firmicutes. This approach offers a powerful means for understanding the influence of drivers of community structure derived from metabarcoding data. We discuss how our approach could be readily adapted to allow researchers to estimate contributions of additional factors such as host or microbe phylogeny to answer emerging questions surrounding the ecological and evolutionary roles of within-host communities. IMPORTANCE NGS and fecal metabarcoding methods have provided powerful opportunities to study the wild gut microbiome. A wealth of data is, therefore, amassing across wild systems, generating the need for analytical approaches that can appropriately investigate simultaneous factors at the host and environmental scale that determine the composition of these communities. Here, we describe a generalized linear mixed-effects model (GLMM) approach to analyze read count data from metabarcoding of the gut microbiota, allowing us to quantify the contributions of multiple host and environmental factors to within-host community structure. Our approach provides outputs that are familiar to a majority of field ecologists and can be run using any standard mixed-effects modeling packages. We illustrate this approach using two metabarcoding data sets from the Soay sheep population of St. Kilda investigating age and season effects as worked examples.

Behavioural individuality determines infection risk in clonal ant colonies

In social groups, infection risk is not distributed evenly across individuals. Individual behaviour is a key source of variation in infection risk, yet its effects are difficult to separate from other factors (e.g., age). Here, we combine epidemiological experiments with chemical, transcriptomic, and automated behavioural analyses in clonal ant colonies, where behavioural individuality emerges among identical workers. We find that: (1) Caenorhabditis-related nematodes parasitise ant heads and affect their survival and physiology, (2) differences in infection emerge from behavioural variation alone, and reflect spatially-organised division of labour, (3) infections affect colony social organisation by causing infected workers to stay in the nest. By disproportionately infecting some workers and shifting their spatial distribution, infections reduce division of labour and increase spatial overlap between hosts, which should facilitate parasite transmission. Thus, division of labour, a defining feature of societies, not only shapes infection risk and distribution but is also modulated by parasites.

Fine-Scale Landscape Epidemiology: Sarcoptic Mange in Bare-Nosed Wombats (Vombatus ursinus)

Landscape epidemiology provides a valuable framework to interpret, predict, and manage spatiotemporal patterns of disease. Yet, owing to the difficulty of detecting pathogen occurrence in free-ranging wildlife, disentangling the factors driving disease dynamics remains a considerable challenge, particularly at fine spatial scales. Here, we investigated the fine-scale landscape epidemiology of sarcoptic mange-a visually apparent disease caused by the mite Sarcoptes scabiei-in bare-nosed wombats (Vombatus ursinus), by: (1) characterizing the distribution and density of wombats within the landscape and (2) examining the effect of environmental variation on the occurrence and apparent prevalence of mange. Wombats were heterogeneously distributed over 19.4 km of transect space (0.096-1.39 wombats ha-1) and seven months of time (increasing by a factor of 1.76). Wombat density was negatively associated with distance to vegetation cover, supporting a general propensity for wombats to occur and burrow near vegetation (native and exotic, excluding pasture). The apparent prevalence of mange varied spatially (3.1-37.5%), with the probability of disease greater in wombats with minimal vegetation and low-lying pans in their estimated home range. We observed trends of increased prevalence in areas with more burrows available per wombat and in individuals occurring near vegetation cover (although not within their home range). Wombat density and active burrow density did not influence the prevalence of mange. This research emphasizes the fine scale at which spatiotemporal patterns of disease can manifest and is the first to investigate the influence of host density for any species with indirect transmission of S. scabiei. Collectively, our results suggest that individuals inhabiting less optimal habitat (pasture) may be at greater risk of disease, or that diseased wombats may be competitively excluded from more optimal habitat (vegetated areas). We discuss implications for understanding and managing mange in wombats and cross-applicability to other mange-affected species with environmental transmission.

The Unusual Value of Long-Term Studies of Individuals: The Example of the Isle of Rum Red Deer Project

Long-term studies of individuals enable incisive investigations of questions across ecology and evolution. Here, we illustrate this claim by reference to our long-term study of red deer on the Isle of Rum, Scotland. This project has established many of the characteristics of social organization, selection, and population ecology typical of large, polygynous, seasonally breeding mammals, with wider implications for our understanding of sexual selection and the evolution of sex differences, as well as for their population dynamics and population management. As molecular genetic techniques have developed, the project has pivoted to investigate evolutionary genetic questions, also breaking new ground in this field. With ongoing advances in genomics and statistical approaches and the development of increasingly sophisticated ways to assay new phenotypic traits, the questions that long-term studies such as the red deer study can answer become both broader and ever more sophisticated. They also offer powerful means of understanding the effects of ongoing climate change on wild populations.

Forecasting parasite sharing under climate change

Species are shifting their distributions in response to climate change. This geographic reshuffling may result in novel co-occurrences among species, which could lead to unseen biotic interactions, including the exchange of parasites between previously isolated hosts. Identifying potential new host-parasite interactions would improve forecasting of disease emergence and inform proactive disease surveillance. However, accurate predictions of future cross-species disease transmission have been hampered by the lack of a generalized approach and data availability. Here, we propose a framework to predict novel host-parasite interactions based on a combination of niche modelling of future host distributions and parasite sharing models. Using the North American ungulates as a proof of concept, we show this approach has high cross-validation accuracy in over 85% of modelled parasites and find that more than 34% of the host-parasite associations forecasted by our models have already been recorded in the literature. We discuss potential sources of uncertainty and bias that may affect our results and similar forecasting approaches, and propose pathways to generate increasingly accurate predictions. Our results indicate that forecasting parasite sharing in response to shifts in host geographic distributions allow for the identification of regions and taxa most susceptible to emergent pathogens under climate change. This article is part of the theme issue 'Infectious disease macroecology: parasite diversity and dynamics across the globe'.

The role of maternally transferred antibodies in maternal performance in red deer

Maternal effects are ubiquitous. Yet, the pathways through which maternal effects occur in wild mammals remain largely unknown. We hypothesise that maternal immune transfer is a key mechanism by which mothers can affect their offspring fitness, and that individual variation in maternally derived antibodies mainly depends on a mother's characteristics and the environmental conditions she experiences. To test this, we assayed six colostrum-derived antibodies in the plasma of 1447 neonates in a wild red deer population. Neonatal antibody levels were mainly affected by maternal genes, environmental variation and costs of prior reproductive investment. We found consistent heterogeneity in maternal performance across traits, with mothers producing the heaviest calves also having calves with more antibodies. Unexpectedly, antibody levels were not associated with calf survival. We provide a unique example of how evolutionary theory on maternal effects can be used to gain insight into the causes of maternal effects in wild populations.

Land use-induced spillover: a call to action to safeguard environmental, animal, and human health

The rapid global spread and human health impacts of SARS-CoV-2, the virus that causes COVID-19, show humanity's vulnerability to zoonotic disease pandemics. Although anthropogenic land use change is known to be the major driver of zoonotic pathogen spillover from wildlife to human populations, the scientific underpinnings of land use-induced zoonotic spillover have rarely been investigated from the landscape perspective. We call for interdisciplinary collaborations to advance knowledge on land use implications for zoonotic disease emergence with a view toward informing the decisions needed to protect human health. In particular, we urge a mechanistic focus on the zoonotic pathogen infect-shed-spill-spread cascade to enable protection of landscape immunity-the ecological conditions that reduce the risk of pathogen spillover from reservoir hosts-as a conservation and biosecurity priority. Results are urgently needed to formulate an integrated, holistic set of science-based policy and management measures that effectively and cost-efficiently minimise zoonotic disease risk. We consider opportunities to better institute the necessary scientific collaboration, address primary technical challenges, and advance policy and management issues that warrant particular attention to effectively address health security from local to global scales.

Fitness Costs of Parasites Explain Multiple Life-History Trade-Offs in a Wild Mammal

AbstractReproduction in wild animals can divert limited resources away from immune defense, resulting in increased parasite burdens. A long-standing prediction of life-history theory states that these parasites can harm the reproductive individual, reducing its subsequent survival and fecundity, producing reproduction-fitness trade-offs. Here, we examined associations among reproductive allocation, immunity, parasitism, and subsequent survival and fecundity in a wild population of individually identified red deer (Cervus elaphus). Using path analysis, we investigated whether costs of lactation in terms of downstream survival and fecundity were mediated by changes in strongyle nematode count and mucosal antibody levels. Lactating females exhibited increased parasite counts, which were in turn associated with substantially decreased fitness in the following year in terms of overwinter survival, fecundity, subsequent calf weight, and parturition date. This study offers observational evidence for parasite regulation of multiple life-history trade-offs, supporting the role of parasites as an important mediating factor in wild mammal populations.

Negative density-dependent parasitism in a group-living carnivore

Animals living at high population densities commonly experience greater exposure to disease, leading to increased parasite burdens. However, social animals can benefit immunologically and hygienically from cooperation, and individuals may alter their socio-spatial behaviour in response to infection, both of which could counteract density-related increases in exposure. Consequently, the costs and benefits of sociality for disease are often uncertain. Here, we use a long-term study of a wild European badger population (Meles meles) to investigate how within-population variation in host density determines infection with multiple parasites. Four out of five parasite taxa exhibited consistent spatial hotspots of infection, which peaked among badgers living in areas of low local population density. Combined movement, survival, spatial and social network analyses revealed that parasite avoidance was the likely cause of this negative density dependence, with possible roles for localized mortality, encounter-dilution effects, and micronutrient-enhanced immunity. These findings demonstrate that animals can organize their societies in space to minimize parasite infection, with important implications for badger behavioural ecology and for the control of badger-associated diseases.

Personality affects dynamics of an experimental pathogen in little brown bats

Host behaviour can affect host-pathogen dynamics and theory predicts that certain individuals disproportionately infect conspecifics during an epidemic. Consistent individual differences in behaviour, or personality, could influence this variation with the most exploratory or sociable individuals most likely to spread pathogens. We quantified exploration and sociability in little brown bats (Myotis lucifugus) and then experimentally manipulated exposure to a proxy pathogen (i.e. ultraviolet (UV) fluorescent powder) to test two related hypotheses: (i) more sociable and more exploratory individuals would be more likely to transmit infections to other individuals, and (ii) more sociable and more exploratory individuals uninfected with an invading pathogen would be more likely to acquire infections. We captured 10 groups of 16 bats at a time and held each group in an outdoor flight tent equipped with roosting-boxes. We used hole-board and Y-maze tests to quantify exploration and sociability of each bat and randomly selected one individual from each group for 'infection' with non-toxic, UV fluorescent powder. Each group of 10 bats was released into the flight tent for 24 h, which represented an experimental infection trial. After 24 h, we removed bats from the trial, photographed each individual under UV light and quantified infection intensity from digital photographs. As predicted, the exploratory behaviour of the experimentally infected individual was positively correlated with infection intensity in their group-mates, while more exploratory females had higher pathogen acquisition. Our results highlight the potential influence of host personality and sex on pathogen dynamics in wildlife populations.

Avian Health across the Landscape: Nestling Immunity Covaries with Changing Landcover

The landscape composition of an organism's home range or territory should influence aspects of its condition, including measures of immune function. Changes in immunocompetence arising from variation in landcover may provide important links between habitat changes and patterns of disease spread. To establish a baseline understanding for whether immune measures covary with changes in landcover, we examined associations between immunological parameters and landcover composition for adults and nestlings of five shrubland bird species. Specifically, we examined the bacteria-killing ability (BKA) of the blood plasma and profiles of the five avian leukocytes as our measures of immune function, and assessed the proportion of area around each bird's nest that was composed of the four major landcover types in the Midwestern USA: row crop agriculture, developed, forest, and grass/shrub. We performed landcover assessments at 100 and 1000 m radius buffers to identify whether associations between habitat and immune function differed at the two spatial scales. As part of this work, we examined age and species-related immunological variation, as well as associations among the immune parameters. There was little evidence linking variation in immune function to landcover composition for the adults at either spatial scale, but there were numerous associations for nestlings, and these were stronger at the 1000 than 100 m spatial scale. The proportion of grass/shrub around the nest had the largest impact on immune function, although the effect varied by immune parameter and species. BKA and basophils were inversely associated with grass/shrub for all species, whereas lymphocytes were positively associated with grass/shrub for all species. We also documented species-level differences among adults and nestlings for BKA and all leukocytes except monocytes. As expected, we found that nestlings had reduced levels of BKA, lymphocytes, monocytes, and elevated heterophils compared with adults (except for field sparrow-Spizella pusilla-nestlings, which had higher lymphocytes). Basophils generally did not differ by age class, and eosinophils exhibited species-specific patterns, in which they were higher for nestling American robins (Turdus migratorius) and gray catbirds (Dumetella carolinensis) compared with adults, but lower in the other nestlings. Heterophils and lymphocytes were inversely associated for all species and age classes, and basophil levels were positively associated with BKA across species and age classes. Together, these findings bolster our understanding of age and species-specific variation in immune function, and provide evidence that immune measures can covary with changes in landcover.

Multi-Scale Drivers of Immunological Variation and Consequences for Infectious Disease Dynamics

The immune system is the primary barrier to parasite infection, replication, and transmission following exposure, and variation in immunity can accordingly manifest in heterogeneity in traits that govern population-level infectious disease dynamics. While much work in ecoimmunology has focused on individual-level determinants of host immune defense (e.g., reproductive status and body condition), an ongoing challenge remains to understand the broader evolutionary and ecological contexts of this variation (e.g., phylogenetic relatedness and landscape heterogeneity) and to connect these differences into epidemiological frameworks. Ultimately, such efforts could illuminate general principles about the drivers of host defense and improve predictions and control of infectious disease. Here, we highlight recent work that synthesizes the complex drivers of immunological variation across biological scales of organization and scales these within-host differences to population-level infection outcomes. Such studies note the limitations involved in making species-level comparisons of immune phenotypes, stress the importance of spatial scale for immunology research, showcase several statistical tools for translating within-host data into epidemiological parameters, and provide theoretical frameworks for linking within- and between-host scales of infection processes. Building from these studies, we highlight several promising avenues for continued work, including the application of machine learning tools and phylogenetically controlled meta-analyses to immunology data and quantifying the joint spatial and temporal dependencies in immune defense using range expansions as model systems. We also emphasize the use of organismal traits (e.g., host tolerance, competence, and resistance) as a way to interlink various scales of analysis. Such continued collaboration and disciplinary cross-talk among ecoimmunology, disease ecology, and mathematical modeling will facilitate an improved understanding of the multi-scale drivers and consequences of variation in host defense.

Leukocyte Profiles Reflect Geographic Range Limits in a Widespread Neotropical Bat

Quantifying how the environment shapes host immune defense is important for understanding which wild populations may be more susceptible or resistant to pathogens. Spatial variation in parasite risk, food and predator abundance, and abiotic conditions can each affect immunity, and these factors can also manifest at both local and biogeographic scales. Yet identifying predictors and the spatial scale of their effects is limited by the rarity of studies that measure immunity across many populations of broadly distributed species. We analyzed leukocyte profiles from 39 wild populations of the common vampire bat (Desmodus rotundus) across its wide geographic range throughout the Neotropics. White blood cell differentials varied spatially, with proportions of neutrophils and lymphocytes varying up to six-fold across sites. Leukocyte profiles were spatially autocorrelated at small and very large distances, suggesting that local environment and large-scale biogeographic factors influence cellular immunity. Generalized additive models showed that bat populations closer to the northern and southern limits of the species range had more neutrophils, monocytes, and basophils, but fewer lymphocytes and eosinophils, than bats sampled at the core of their distribution. Habitats with access to more livestock also showed similar patterns in leukocyte profiles, but large-scale patterns were partly confounded by time between capture and sampling across sites. Our findings suggest that populations at the edge of their range experience physiologically limiting conditions that predict higher chronic stress and greater investment in cellular innate immunity. High food abundance in livestock-dense habitats may exacerbate such conditions by increasing bat density or diet homogenization, although future spatially and temporally coordinated field studies with common protocols are needed to limit sampling artifacts. Systematically assessing immune function and response over space will elucidate how environmental conditions influence traits relevant to epidemiology and help predict disease risks with anthropogenic disturbance, land conversion, and climate change.

Macroimmunology: The drivers and consequences of spatial patterns in wildlife immune defence

The prevalence and intensity of parasites in wild hosts varies across space and is a key determinant of infection risk in humans, domestic animals and threatened wildlife. Because the immune system serves as the primary barrier to infection, replication and transmission following exposure, we here consider the environmental drivers of immunity. Spatial variation in parasite pressure, abiotic and biotic conditions, and anthropogenic factors can all shape immunity across spatial scales. Identifying the most important spatial drivers of immunity could help pre-empt infectious disease risks, especially in the context of how large-scale factors such as urbanization affect defence by changing environmental conditions. We provide a synthesis of how to apply macroecological approaches to the study of ecoimmunology (i.e. macroimmunology). We first review spatial factors that could generate spatial variation in defence, highlighting the need for large-scale studies that can differentiate competing environmental predictors of immunity and detailing contexts where this approach might be favoured over small-scale experimental studies. We next conduct a systematic review of the literature to assess the frequency of spatial studies and to classify them according to taxa, immune measures, spatial replication and extent, and statistical methods. We review 210 ecoimmunology studies sampling multiple host populations. We show that whereas spatial approaches are relatively common, spatial replication is generally low and unlikely to provide sufficient environmental variation or power to differentiate competing spatial hypotheses. We also highlight statistical biases in macroimmunology, in that few studies characterize and account for spatial dependence statistically, potentially affecting inferences for the relationships between environmental conditions and immune defence. We use these findings to describe tools from geostatistics and spatial modelling that can improve inference about the associations between environmental and immunological variation. In particular, we emphasize exploratory tools that can guide spatial sampling and highlight the need for greater use of mixed-effects models that account for spatial variability while also allowing researchers to account for both individual- and habitat-level covariates. We finally discuss future research priorities for macroimmunology, including focusing on latitudinal gradients, range expansions and urbanization as being especially amenable to large-scale spatial approaches. Methodologically, we highlight critical opportunities posed by assessing spatial variation in host tolerance, using metagenomics to quantify spatial variation in parasite pressure, coupling large-scale field studies with small-scale field experiments and longitudinal approaches, and applying statistical tools from macroecology and meta-analysis to identify generalizable spatial patterns. Such work will facilitate scaling ecoimmunology from individual- to habitat-level insights about the drivers of immune defence and help predict where environmental change may most alter infectious disease risk.