Metapopulation-level adaptation of insect host plant preference and extinction-colonization dynamics in heterogeneous landscapes

Species living in highly fragmented landscapes typically occur as metapopulations with frequent turnover of local populations. The turnover rate depends on population sizes and connectivities, but it may also depend on the phenotypic and genotypic composition of populations. The Glanville fritillary butterfly (Melitaea cinxia) in Finland uses two host plant species, which show variation in their relative abundances at two spatial scales: locally among individual habitat patches and regionally among networks of patches. Female butterflies in turn exhibit spatial variation in genetically determined host plant preference within and among patch networks. Emigration, immigration and establishment of new populations have all been shown to be strongly influenced by the match between the host plant composition of otherwise suitable habitat patches and the host plant preference of migrating butterflies. The evolutionary consequences of such biased migration and colonization with respect to butterfly phenotypes might differ depending on spatial configuration and plant species composition of the patches in heterogeneous patch networks. Using a spatially realistic individual-based model we show that the model-predicted evolution of host plant preference due to biased migration explains a significant amount of the observed variation in host plant use among metapopulations living in dissimilar networks. This example illustrates how the ecological extinction-colonization dynamics may be linked with the evolutionary dynamics of life history traits in metapopulations.     

Heterogeneity of plant phenotypes caused by herbivore-specific induced responses influences the spatial distribution of herbivores

Abstract.  1. Herbivory can induce resistance in a plant and the induced phenotype may be disfavoured by subsequent herbivores. Yet, as the distance between plants in a population increases, limited mobility may make a herbivore more likely to feed and oviposit on host plants in its immediate surroundings.
2. The present study tested whether a herbivore's preference and distribution across plants with different induced phenotypes was influenced by the spatial distribution of plants. A fragmented population of Solanum dulcamara plants was created. This consisted of discrete, spatially separated patches with different histories of damage, either herbivory from adult flea beetles ( Psylliodes affinis ), tortoise beetles ( Plagiometriona clavata ), or mechanical damage. Each patch was separated by 7 m and consisted of 12 plants that were spaced 30 cm apart. Then a fixed number of adult tortoise beetles were introduced to each patch, and movement and oviposition within and between spatially separate homogeneous patches (receiving one type of damage) were compared with movement and oviposition within heterogeneous patches (containing all three types of damage) over the growing season.
3. Flea beetle and tortoise beetle herbivory consistently induced different phytochemical responses in S. dulcamara (polyphenol oxidase and peroxidase), and adult tortoise beetles avoided oviposition on the flea beetle induced plants within heterogeneous patches. However, between homogeneous patches, plant phenotype did not influence oviposition. Colonisation by naturally occurring flea beetle adults followed a similar pattern.
4. These results suggest that the heterogeneity of plant phenotypes can influence herbivore choice and distribution at small but not large spatial scales.     

Modelling the spatial dynamics and persistence of the leaf beetle Gonioctena olivacea in dynamic habitats

In dynamic landscapes natural and anthropogenic disturbance as well as succession are responsible for the emergence and subsequent disappearance of suitable habitat patches. Species inhabiting such landscapes are faced with varying number and spatial configuration of patches. A stochastic, spatially explicit simulation model was developed in order to analyse the persistence of the leaf beetle Gonioctena olivacea in a system of dynamic patches of its host plant Cytisus scoparius . The model was parameterized with data from a three-year field study on the spatial configuration, distribution, and turnover of the host plant patches as well as the patch occupancy, extinction, and colonization rates of the beetle. The simulations showed large fluctuations in the occurrence of the beetle in the patches. High levels of occupancy were related to high aggregation of the patches within the landscape. The velocity of patch turnover was found to have a severe effect on the persistence of the beetle metapopulation. Enhancing the turnover rate by only a few patches, the mean time to extinction decreases rapidly. Moreover, the results revealed that not necessarily an effect of connectivity can be detected in the analysis of occupancy patterns in dynamic landscapes, although the colonization of patches is clearly connectivity-dependent. In general, this modelling study demonstrates the importance of detailed information on patch turnover. The amount and spatial distribution of suitable habitat is a major driver of metapopulation dynamics of species in dynamic landscapes.     

Spatiotemporal Model of Barley and Cereal Yellow Dwarf Virus Transmission Dynamics with Seasonality and Plant Competition

Many generalist pathogens are influenced by the spatial distributions and relative abundances of susceptible host species. The spatial structure of host populations can influence patterns of infection incidence (or disease outbreaks), and the effects of a generalist pathogen on host community dynamics in a spatially heterogeneous community may differ from predictions derived via simple models. In this paper, we model the transmission of a generalist pathogen within a patch framework that incorporates the movement of vectors between discrete host patches to investigate the effects of local host community composition and vector movement rates on disease dynamics.     

Spatially structured metapopulation models: global and local assessment of metapopulation capacity

We model metapopulation dynamics in finite networks of discrete habitat patches with given areas and spatial locations. We define and analyze two simple and ecologically intuitive measures of the capacity of the habitat patch network to support a viable metapopulation. Metapopulation persistence capacity lambda(M) defines the threshold condition for long-term metapopulation persistence as lambda(M)>delta, where delta is defined by the extinction and colonization rate parameters of the focal species. Metapopulation invasion capacity lambda(I) sets the condition for successful invasion of an empty network from one small local population as lambda(I)>delta. The metapopulation capacities lambda(M) and lambda(I) are defined as the leading eigenvalue or a comparable quantity of an appropriate "landscape" matrix. Based on these definitions, we present a classification of a very general class of deterministic, continuous-time and discrete-time metapopulation models. Two specific models are analyzed in greater detail: a spatially realistic version of the continuous-time Levins model and the discrete-time incidence function model with propagule size-dependent colonization rate and a rescue effect. In both models we assume that the extinction rate increases with decreasing patch area and that the colonization rate increases with patch connectivity. In the spatially realistic Levins model, the two types of metapopulation capacities coincide, whereas the incidence function model possesses a strong Allee effect characterized by lambda(I)=0. For these two models, we show that the metapopulation capacities can be considered as simple sums of contributions from individual habitat patches, given by the elements of the leading eigenvector or comparable quantities. We may therefore assess the significance of particular habitat patches, including new patches that might be added to the network, for the metapopulation capacities of the network as a whole. We derive useful approximations for both the threshold conditions and the equilibrium states in the two models. The metapopulation capacities and the measures of the dynamic significance of particular patches can be calculated for real patch networks for applications in metapopulation ecology, landscape ecology, and conservation biology.     

Host plant density and patch isolation drive occupancy and abundance at a butterfly's northern range margin

Marginal populations are usually small, fragmented, and vulnerable to extinction, which makes them particularly interesting from a conservation point of view. They are also the starting point of range shifts that result from climate change, through a process involving colonization of newly suitable sites at the cool margin of species distributions. Hence, understanding the processes that drive demography and distribution at high‐latitude populations is essential to forecast the response of species to global changes. We investigated the relative importance of solar irradiance (as a proxy for microclimate), habitat quality, and connectivity on occupancy, abundance, and population stability at the northern range margin of the Oberthür's grizzled skipper butterfly Pyrgus armoricanus. For this purpose, butterfly abundance was surveyed in a habitat network consisting of 50 habitat patches over 12 years. We found that occupancy and abundance (average and variability) were mostly influenced by the density of host plants and the spatial isolation of patches, while solar irradiance and grazing frequency had only an effect on patch occupancy. Knowing that the distribution of host plants extends further north, we hypothesize that the actual variable limiting the northern distribution of P. armoricanus might be its dispersal capacity that prevents it from reaching more northern habitat patches. The persistence of this metapopulation in the face of global changes will thus be fundamentally linked to the maintenance of an efficient network of habitats.     

Fine-Grained Distribution of a Non-Native Resource Can Alter the Population Dynamics of a Native Consumer

New interactions with non-native species can alter selection pressures on native species. Here, we examined the effect of the spatial distribution of a non-native species, a factor that determines ecological and evolutionary outcomes but that is poorly understood, particularly on a fine scale. Specifically, we explored a native butterfly population and a non-native plant on which the butterfly oviposits despite the plant’s toxicity to larvae. We developed an individual-based model to describe movement and oviposition behaviors of each butterfly, which were determined by plant distribution and the butterfly''s host preference genotype. We estimated the parameter values of the model from rich field data. We simulated various patterns of plant distributions and compared the rates of butterfly population growth and changes in the allele frequency of oviposition preference. Neither the number nor mean area of patches of non-native species affected the butterfly population, whereas plant abundance, patch shape, and distance to the nearest native and non-native patches altered both the population dynamics and genetics. Furthermore, we found a dramatic decrease in population growth rates when we reduced the distance to the nearest native patch from 147 m to 136 m. Thus changes in the non-native resource distribution that are critical to the fate of the native herbivore could only be detected at a fine-grained scale that matched the scale of a female butterfly’s movement. In addition, we found that the native butterfly population was unlikely to be rescued by the exclusion of the allele for acceptance of the non-native plant as a host. This study thus highlights the importance of including both ecological and evolutionary dynamics in analyses of the outcome of species interactions and provides insights into habitat management for non-native species.     

Effects of Culling on Mesopredator Population Dynamics

Anthropogenic changes in land use and the extirpation of apex predators have facilitated explosive growth of mesopredator populations. Consequently, many species have been subjected to extensive control throughout portions of their range due to their integral role as generalist predators and reservoirs of zoonotic disease. Yet, few studies have monitored the effects of landscape composition or configuration on the demographic or behavioral response of mesopredators to population manipulation. During 2007 we removed 382 raccoons (Procyon lotor) from 30 forest patches throughout a fragmented agricultural ecosystem to test hypotheses regarding the effects of habitat isolation on population recovery and role of range expansion and dispersal in patch colonization of mesopredators in heterogeneous landscapes. Patches were allowed to recolonize naturally and demographic restructuring of patches was monitored from 2008–2010 using mark-recapture. An additional 25 control patches were monitored as a baseline measure of demography. After 3 years only 40% of experimental patches had returned to pre-removal densities. This stagnant recovery was driven by low colonization rates of females, resulting in little to no within-patch recruitment. Colonizing raccoons were predominantly young males, suggesting that dispersal, rather than range expansion, was the primary mechanism driving population recovery. Contrary to our prediction, neither landscape connectivity nor measured local habitat attributes influenced colonization rates, likely due to the high dispersal capability of raccoons and limited role of range expansion in patch colonization. Although culling is commonly used to control local populations of many mesopredators, we demonstrate that such practices create severe disruptions in population demography that may be counterproductive to disease management in fragmented landscapes due to an influx of dispersing males into depopulated areas. However, given the slow repopulation rates observed in our study, localized depopulation may be effective at reducing negative ecological impacts of mesopredators in fragmented landscapes at limited spatial and temporal scales.     

How Do Landscape Structure,Management and Habitat Quality Drive the Colonization of Habitat Patches by the Dryad Butterfly (Lepidoptera: Satyrinae) in Fragmented Grassland?

Most studies dealing with species distribution patterns on fragmented landscapes focus on the characteristics of habitat patches that influence local occurrence and abundance, but they tend to neglect the question of what drives colonization of previously unoccupied patches. In a study of the dryad butterfly, we combined classical approaches derived from metapopulation theory and landscape ecology to investigate the factors driving colonization from a recent refugium. In three consecutive transect surveys, we recorded the presence and numbers of imagos in 27 patches of xerothermic grassland and 26 patches of wet meadow. Among the predictors affecting the occurrence and abundance of the dryad, we considered environmental variables reflecting (i) habitat patch quality (e.g., goldenrod cover, shrub density, vegetation height); (ii) factors associated with habitat spatial structure (patch size, patch isolation and fragmentation); and (iii) features of patch surroundings (100-m buffers around patches) that potentially pose barriers or provide corridors. Patch colonization by the dryad was strongly limited by the distance from the species refugium in the region; there was a slight positive effect of shrub density in this respect. Butterfly abundance increased in smaller and more fragmented habitat patches; it was negatively impacted by invasive goldenrod cover, and positively influenced by the density of watercourses in patch surroundings. Nectar plant availability was positively related to species abundance in xerothermic grassland, while in wet meadow the effect was the reverse. We conclude that dryad colonization of our study area is very recent, since the most important factor limiting colonization was distance from the refugium, while the habitat quality of target patches had less relevance. In order to preserve the species, conservation managers should focus on enhancing the quality of large patches and should also direct their efforts on smaller and more fragmented ones, including those with relatively low resource availability, because such habitat fragments have an important role to play for specialist species.     

Patch occupancy of two hemipterans sharing a common host plant

Aim Two hemipteran species were chosen as a study system for the comparative analysis of patch occupancy and spatial population structure of insects sharing a common host plant. This study tested whether (1) the incidence in the host plant patches differed between the two species, and (2) the two species exhibited a different spatial population structure, i.e. were they affected differentially by isolation and area of the host plant patches. Location The porphyry landscape north of Halle (Saale) in Germany comprising 506 patches of the host plant Brachypodium pinnatum. Methods The host plant patches were surveyed for the two hemipterans. To assess the influence of patch quality on species occurrence the patches were characterized by mean cover abundance of B. pinnatum, type of subsoil, slope, exposure, and shading. The spatial configuration of the patches was considered by patch area and isolation. The influence of the habitat factors and the spatial configuration on the occupancy of the two species was analysed by logistic regression. Results Adarrus multinotatus was found in 441 patches, while Neophilaenus albipennis was found in only 90 patches. While A. multinotatus showed virtually no relationship to the habitat factors, the occupancy of N. albipennis was influenced by subsoil type, cover abundance, and shading. The effects of area and isolation on occupancy of the patches also differed between the two species. The occupancy of N. albipennis was determined largely by area and isolation, whereas in A. multinotatus no considerable effect of spatial configuration was found. Main conclusions The study revealed a marked difference between the two hemipteran species in respect of spatial population structure. Adarrus multinotatus built up a ‘patchy population’, whereas N. albipennis showed a ‘metapopulation’ structure within the same set of patches in the same landscape. Spatial population structure was found to be not only a function of spatial configuration of habitat patches, but population structure differed between the habitat generalist A. multinotatus and the habitat specialist N. albipennis.     

Host patch traits have scale-dependent effects on diversity in a stickleback parasite metacommunity

Many metacommunities are distributed across habitat patches that are themselves aggregated into groups. Perhaps the clearest example of this nested metacommunity structure comes from multi-species parasite assemblages, which occupy individual hosts that are aggregated into host populations. At both spatial scales, we expect parasite community diversity in a given patch (either individual host or population) to depend on patch characteristics that affect colonization rates and species sorting. But, are these patch effects consistent across spatial scales? Or, do different processes govern the distribution of parasite community diversity among individual hosts, versus among host patches? To answer these questions, we document the distribution of parasite richness among host individuals and among populations in a metapopulation of threespine stickleback Gasterosteus aculeatus. We find some host traits (host size, gape width) are associated with increased parasite richness at both spatial scales. Other patch characteristics affect parasite richness only among individuals (sex), or among populations (lake size, lake area, elevation and population mean heterozygosity). These results demonstrate that some rules governing parasite richness in this metacommunity are shared across scales, while others are scale-specific.     

The effective size of a metapopulation living in a heterogeneous patch network

I analyze stochastic patch occupancy models (SPOMs), which record habitat patches as empty or occupied. A problem with SPOMs has been that if the spatial structure of a heterogeneous habitat patch network is taken into account, the computational effort needed to analyze a SPOM grows as a power of 2n, where n is the number of habitat patches. I propose a computationally feasible approximation method, which approximates the behavior of a heterogeneous SPOM by an "ideal" metapopulation inhabiting a network of identical and equally connected habitat patches. The transformation to the ideal metapopulation is based on weighting the individual patch occupancies by the dynamic values of the habitat patches, which may be calculated from the deterministic mean-field approximation of the original SPOM. Conceptually, the method resembles the calculation of the effective size of a population in the context of population genetics. I demonstrate how the method may be applied to SPOMs with flexible structural assumptions and with spatially correlated and temporally varying parameter values. I apply the method to a real habitat patch network inhabited by the Glanville fritillary butterfly, illustrating that the metapopulation dynamics of this species are essentially driven by temporal variability in the environmental conditions.     

Incidence and population dynamics of the leaf beetle Gonioctena olivacea in dynamic habitats

In natural as well as in cultural landscapes, disturbance and succession are responsible for the emergence and subsequent disappearance of suitable habitat patches. The dynamics of habitat patches has important consequences for the spatial structure and dynamics of regional populations. However, there are only few studies quantifying both patch dynamics and incidence of insect species in a dynamic landscape over several years. I studied the incidence and population dynamics of the leaf beetle Gonioctena olivacea in a system of dynamic patches of the host plant Scotch broom Cytisus scoparius . The incidence of the beetle was most strongly affected by patch area, whereas connectivity, patch quality, patch age, and landscape context had no or only a minor effect when analysed with logistic regression. The size of local beetle populations was highly fluctuating between the years; however, the population dynamics of the local populations was not synchronous. Adjacent patches did not show higher degrees of synchrony than patches separated by large distances. In the three years of study, local populations became extinct through demographic or environmental stochasticity and patch destruction. Each year >10% of the patches disappeared. The extinction rate of beetles in persistent patches was decreasing with increasing patch area. On the other hand, patches newly emerged and were rapidly colonized by the beetle. The colonization rate depended on patch connectivity. Obviously, Gonioctena olivacea was capable of persisting in this system with high turnover of patches owing to its high dispersal power.     

Population level correlation between pre-alighting and post-alighting host plant preference in the Glanville fritillary butterfly

1. In many populations of the Glanville fritillary butterfly Melitaea cinxia, ovipositing females exhibit a post‐alighting preference for one of the potential host plant species available. The work reported here aimed to establish whether females with different post‐alighting preferences can discriminate between their host plant species prior to alighting, and whether pre‐alighting and post‐alighting preferences are correlated at the population level. 2. Alighting and oviposition events were recorded for groups of females from six populations in greenhouse and field experiments. 3. Landing frequencies did not change with experience, indicating that M. cinxia females did not learn from previous encounters with host plants. 4. Females from populations exhibiting post‐alighting preference searched efficiently for their host plants in the sense that they landed mainly on the species on which they oviposited predominantly. Pre‐alighting and post‐alighting preferences were correlated at the population level. 5. The correlation between pre‐alighting and post‐alighting preferences helps to explain why in nature, where the host plants often occur in distinct patches, females are more likely to colonise habitat patches in which their preferred host plant is abundant.     

Using connectivity metrics in conservation planning – when does habitat quality matter?

Aim The objective of conservation planning is often to prioritize patches based on their estimated contribution to metapopulation or metacommunity viability. The contribution that an individual patch makes will depend on its intrinsic characteristics, such as habitat quality, as well as its location relative to other patches, its connectivity. Here we systematically evaluate five patch value metrics to determine the importance of including an estimate of habitat quality into the metrics. Location We tested the metrics in landscapes designed to represent different degrees of variability in patch quality and different levels of patch aggregation. Methods In each landscape, we simulated population dynamics using a spatially explicit, continuous time metapopulation model linked to within patch logistic growth models. We tested five metrics that are used to estimate the contribution that a patch makes to metapopulation viability: two versions of the probability of connectivity index, two versions of patch centrality (a graph theory metric) and the metapopulation capacity metric. Results All metrics performed best in environments where patch quality was very variable and high quality patches were aggregated. Metrics that incorporated some measure of patch quality did better in all environments, but did particularly well in environments with high variance of patch quality and spatial aggregation of good quality patches. Main conclusions Including an estimate of patch quality significantly increased the ability of a given connectivity metric to rank correctly habitat patches according to their contribution to metapopulation viability. Incorporating patch quality is particularly important in landscapes where habitat quality is highly variable and good quality patches are spatially aggregated. However, caution should be used when applying patch metrics to homogeneous landscapes, even if good estimates of patch quality are available. Our results demonstrate that landscape structure and the degree of variability in patch quality need to be assessed prior to selecting a suitable method for estimating patch value.     

Homogenization techniques for population dynamics in strongly heterogeneous landscapes

An important problem in spatial ecology is to understand how population-scale patterns emerge from individual-level birth, death, and movement processes. These processes, which depend on local landscape characteristics, vary spatially and may exhibit sharp transitions through behavioural responses to habitat edges, leading to discontinuous population densities. Such systems can be modelled using reaction–diffusion equations with interface conditions that capture local behaviour at patch boundaries. In this work we develop a novel homogenization technique to approximate the large-scale dynamics of the system. We illustrate our approach, which also generalizes to multiple species, with an example of logistic growth within a periodic environment. We find that population persistence and the large-scale population carrying capacity is influenced by patch residence times that depend on patch preference, as well as movement rates in adjacent patches. The forms of the homogenized coefficients yield key theoretical insights into how large-scale dynamics arise from the small-scale features.     

Is the population turnover of patchy‐distributed annuals determined by dormancy dynamics or dispersal processes?

In species with fragmented distribution, regional turnover dynamics is given by the processes of local population extinction and patch (re)colonization by migrants spreading from neighboring occupied patches. In plants with dormant stages (e.g. seeds) and limited dispersal capacity, regional dynamics based on dispersal processes can be overridden by pseudo-turnover determined by signals inducing or breaking dormancy (e.g. due to changes in habitat quality) resulting in a low importance of habitat configuration and size.
In this study, I investigated the turnover dynamics of 5 annual plant species growing on ant mounds of Lasius flavus over three years. I analyzed whether the grassland-scale dynamics of these annuals is influenced by dispersal processes, or alternatively, by pseudo-turnover of soil seed populations. For that purpose I 1) searched for populations formed from soil seeds only, 2) compared the relative contribution of the soil seed bank and seed rain for population restoration after disappearance from the vegetation and 3) investigated whether colonization and extinction events are affected by patch isolation. I assumed if population turnover was rather a result of the soil seed bank dynamics then spatial effects would be hard to detect.
In spite of the presence of populations formed from soil seed and the relatively more important soil seed bank for potential population reestablishment, turnover dynamics followed the predictions of metapopulation theory. Population appearance was more probable in larger and less isolated patches. Probability of disappearance increased with decrease of population size that was negatively influenced by the patch size and its isolation. These findings indicate dispersal processes to be important in the turnover dynamics and only limited contribution of soil seed populations. Their small effectiveness is probably related to the low chance of recurrent disturbance on the mound surface.     

Habitat destruction in mutualistic metacommunities

We investigate a mutualistic metacommunity where the strength of the mutualistic interaction between species is measured by the extent to which the presence of one species on a patch either reduces the extinction rate of the others present on the same patch or increases their ability to colonize other patches. In both cases, a strong enough mutualism enables all species to persist at habitat densities where they would all be extinct in the absence of the interaction. However, a mutualistic interaction that enhances colonization enables the species to persist at lower habitat density than one that suppresses extinction. All species abruptly go extinct (catastrophe) when the habitat density is decreased infinitesimally below a critical value. A comparison of the mean field or spatially implicit case with unrestricted dispersal and colonization to all patches in the system with a spatially explicit case where dispersal is restricted to the immediate neighbours of the original patch leads to the intriguing conclusion that restricted dispersal can be favourable for species that have a beneficial effect on each other when habitat conditions are adverse. When the mutualistic interaction is strong enough, the extinction threshold or critical amount of habitat required for the persistence of all species is lower when the dispersal is locally restricted than when unrestricted ! The persistence advantage for all species created by the mutualistic interaction increases substantially with the number of species in the metacommunity, as does the advantage for restricted dispersal over global dispersal.     

The effects of habitat fragmentation via forestry plantation establishment on spatial genotypic structure in the small marsupial carnivore, Antechinus agilis

Dispersal is an important influence on species' distributions, patch colonization and population persistence in fragmented habitat. We studied the impacts of habitat fragmentation resulting from establishment of an exotic pine plantation on dispersal of the marsupial carnivore, Antechinus agilis. We applied spatial analyses of individual multilocus microsatellite genotypes and mitochondrial haplotypes to study patterns of gene flow in fragmented habitat and natural habitat 'control' areas, and how this is affected by the spatial dispersion of habitat patches, the presence of corridors and a 'mainland' source of migrants. Spatial analysis of molecular variance and partial Mantel tests confirmed the absence of cryptic barriers to gene flow in continuous habitat, which if present would confound the comparison of genetic structures in fragmented vs. unfragmented habitats. Spatial genotypic structure suggested that although dispersal was male-biased in both habitat types, fragmentation restricted dispersal of males more than that of females and the degree of restriction of male dispersal was dependent on the geographical isolation of the patch. The scale of positive genotypic structure in fragmented habitat was restricted to the two closest patches for females and the three closest patches for males. Our results provide evidence for significantly increased gene flow through habitat corridors relative to that across the matrix and for significantly lower gene flow between 'mainland' unfragmented habitat and habitat patches relative to that within either habitat type, suggesting a behavioural barrier to crossing habitat interfaces.     

Disease hotspots or hot species? Infection dynamics in multi‐host metacommunities controlled by species identity,not source location

Pathogen persistence in host communities is influenced by processes operating at the individual host to landscape‐level scale, but isolating the relative contributions of these processes is challenging. We developed theory to partition the influence of host species, habitat patches and landscape connectivity on pathogen persistence within metacommunities of hosts and pathogens. We used this framework to quantify the contributions of host species composition and habitat patch identity on the persistence of an amphibian pathogen across the landscape. By sampling over 11 000 hosts of six amphibian species, we found that a single host species could maintain the pathogen in 91% of observed metacommunities. Moreover, this dominant maintenance species contributed, on average, twice as much to landscape‐level pathogen persistence compared to the most influential source patch in a metacommunity. Our analysis demonstrates substantial inequality in how species and patches contribute to pathogen persistence, with important implications for targeted disease management.