Interspecific Competition, Environmental Gradients, Gene Flow, and the Coevolution of Species' Borders

Niche partitioning and individual specialisation in resources and space use of sympatric fur seals at their range margin

Ecological theory predicts niche partitioning between high-level predators living in sympatry as a mechanism to minimise the selective pressure of competition. Accordingly, male Australian fur seals Arctocephalus pusillus doriferus and New Zealand fur seals A. forsteri that live in sympatry should exhibit partitioning in their broad niches (in habitat and trophic dimensions) in order to coexist. However, at the northern end of their distributions in Australia, both are recolonising their historic range after a long absence due to over-exploitation, and their small population sizes suggest competition should be weak and may allow overlap in niche space. We found some niche overlap, yet clear partitioning in diet trophic level (δ15N values from vibrissae), spatial niche space (horizontal and vertical telemetry data) and circadian activity patterns (timing of dives) between males of each species, suggesting competition may remain an active driver of niche partitioning amongst individuals even in small, peripheral populations. Consistent with individual specialisation theory, broad niches of populations were associated with high levels of individual specialisation for both species, despite putative low competition. Specialists in isotopic space were not necessarily specialists in spatial niche space, further emphasising their diverse individual strategies for niche partitioning. Males of each species displayed distinct foraging modes, with Australian fur seals primarily benthic and New Zealand fur seals primarily epipelagic, though unexpectedly high individual specialisation for New Zealand fur seals might suggest marginal populations provide exceptions to the pattern generally observed amongst other fur seals.

Fitness surfaces and local thermal adaptation in Drosophila along a latitudinal gradient

Local adaptation is commonly cited to explain species distribution, but how fitness varies along continuous geographical gradients is not well understood. Here, we combine thermal biology and life-history theory to demonstrate that Drosophila populations along a 2500 km latitudinal cline are adapted to local conditions. We measured how heat tolerance and viability rate across eight populations varied with temperature in the laboratory and then simulated their expected cumulative Darwinian fitness employing high-resolution temperature data from their eight collection sites. Simulations indicate a trade-off between annual survival and cumulative viability, as both mortality and the recruitment of new flies are predicted to increase in warmer regions. Importantly, populations are locally adapted and exhibit the optimal combination of both traits to maximize fitness where they live. In conclusion, our method is able to reconstruct fitness surfaces employing empirical life-history estimates and reconstructs peaks representing locally adapted populations, allowing us to study geographic adaptation in silico.

A Community-Based Framework Integrates Interspecific Interactions into Forest Genetic Conservation

Forest genetic conservation is typically species-specific and does not integrate interspecific interaction and community structure. It mainly focuses on the theories of population and quantitative genetics. This approach depicts the intraspecific patterns of population genetic structure derived from genetic markers and the genetic differentiation of adaptive quantitative traits in provenance trials. However, it neglects possible interspecific interaction in natural forests and overlooks natural hybridization or subspeciation. We propose that the genetic diversity of a given species in a forest community is shaped by both intraspecific population and interspecific community evolutionary processes, and expand the traditional forest genetic conservation concept under the community ecology framework. We show that a community-specific phylogeny derived from molecular markers would allow us to explore the genetic mechanisms of a tree species interacting with other resident species. It would also facilitate the exploration of a species' ecological role in forest community assembly and the taxonomic relationship of the species with other species specific to its resident forest community. Phylogenetic β-diversity would assess the similarities and differences of a tree species across communities regarding ecological function, the strength of selection pressure, and the nature and extent of its interaction with other species. Our forest genetic conservation proposal that integrates intraspecific population and interspecific community genetic variations is suitable for conserving a taxonomic species complex and maintaining its evolutionary potential in natural forests. This provides complementary information to conventional population and quantitative genetics-based conservation strategies.

Understanding the joint evolution of dispersal and host specialisation using phytophagous arthropods as a model group

Theory generally predicts that host specialisation and dispersal should evolve jointly. Indeed, many models predict that specialists should be poor dispersers to avoid landing on unsuitable hosts while generalists will have high dispersal abilities. Phytophagous arthropods are an excellent group to test this prediction, given extensive variation in their host range and dispersal abilities. Here, we explore the degree to which the empirical literature on this group is in accordance with theoretical predictions. We first briefly outline the theoretical reasons to expect such a correlation. We then report empirical studies that measured both dispersal and the degree of specialisation in phytophagous arthropods. We find a correlation between dispersal and levels of specialisation in some studies, but with wide variation in this result. We then review theoretical attributes of species and environment that may blur this correlation, namely environmental grain, temporal heterogeneity, habitat selection, genetic architecture, and coevolution between plants and herbivores. We argue that theoretical models fail to account for important aspects, such as phenotypic plasticity and the impact of selective forces stemming from other biotic interactions, on both dispersal and specialisation. Next, we review empirical caveats in the study of this interplay. We find that studies use different measures of both dispersal and specialisation, hampering comparisons. Moreover, several studies do not provide independent measures of these two traits. Finally, variation in these traits may occur at scales that are not being considered. We conclude that this correlation is likely not to be expected from large-scale comparative analyses as it is highly context dependent and should not be considered in isolation from the factors that modulate it, such as environmental scale and heterogeneity, intrinsic traits or biotic interactions. A stronger crosstalk between theoretical and empirical studies is needed to understand better the prevalence and basis of the correlation between dispersal and specialisation.

Niche differentiation, reproductive interference, and range expansion

Understanding species distributions and predicting future range shifts requires considering all relevant abiotic factors and biotic interactions. Resource competition has received the most attention, but reproductive interference is another widespread biotic interaction that could influence species ranges. Rubyspot damselflies (Hetaerina spp.) exhibit a biogeographic pattern consistent with the hypothesis that reproductive interference has limited range expansion. Here, we use ecological niche models to evaluate whether this pattern could have instead been caused by niche differentiation. We found evidence for climatic niche differentiation, but the species that encounters the least reproductive interference has one of the narrowest and most peripheral niches. These findings strengthen the case that reproductive interference has limited range expansion and also provide a counterexample to the idea that release from negative species interactions triggers niche expansion. We propose that release from reproductive interference enables species to expand in range while specializing on the habitats most suitable for breeding.

Coral Reef Population Genomics in an Age of Global Change

Coral reefs are both exceptionally biodiverse and threatened by climate change and other human activities. Here, we review population genomic processes in coral reef taxa and their importance for understanding responses to global change. Many taxa on coral reefs are characterized by weak genetic drift, extensive gene flow, and strong selection from complex biotic and abiotic environments, which together present a fascinating test of microevolutionary theory. Selection, gene flow, and hybridization have played and will continue to play an important role in the adaptation or extinction of coral reef taxa in the face of rapid environmental change, but research remains exceptionally limited compared to the urgent needs. Critical areas for future investigation include understanding evolutionary potential and the mechanisms of local adaptation, developing historical baselines, and building greater research capacity in the countries where most reef diversity is concentrated.

A Case for the "Competitive Exclusion-Tolerance Rule" as a General Cause of Species Turnover along Environmental Gradients

AbstractClosely related, ecologically similar species often segregate their distributions along environmental gradients of time, space, and resources, but previous research suggests diverse underlying causes. Here, we review reciprocal removal studies in nature that experimentally test the role of interactions among species in determining their turnover along environmental gradients. We find consistent evidence for asymmetric exclusion coupled with differences in environmental tolerance causing the segregation of species pairs, where a dominant species excludes a subordinate from benign regions of the gradient but is unable to tolerate challenging regions to which the subordinate species is adapted. Subordinate species were consistently smaller and performed better in regions of the gradient typically occupied by the dominant species compared with their native distribution. These results extend previous ideas contrasting competitive ability with adaptation to abiotic stress to include a broader diversity of species interactions (intraguild predation, reproductive interference) and environmental gradients, including gradients of biotic challenge. Collectively, these findings suggest that adaptation to environmental challenge compromises performance in antagonistic interactions with ecologically similar species. The consistency of this pattern across diverse organisms, environments, and biomes suggests generalizable processes structuring the segregation of ecologically similar species along disparate environmental gradients, a phenomenon that we propose should be named the competitive exclusion-tolerance rule.

Genome-environment associations along elevation gradients in two snowbed species of the North-Eastern Calcareous Alps

BACKGROUND

Anthropogenic climate change leads to increasing temperatures and altered precipitation and snowmelt patterns, especially in alpine ecosystems. To understand species' responses to climate change, assessment of genetic structure and diversity is crucial as the basis for the evaluation of migration patterns, genetic adaptation potential as well as the identification of adaptive alleles.

RESULTS

We studied genetic structure, diversity and genome-environment associations of two snowbed species endemic to the Eastern Alps with a large elevational range, Achillea clusiana Tausch and Campanula pulla L. Genotyping-by-sequencing was employed to assemble loci de novo, call variants and perform population genetic analyses. Populations of either species were distinguishable by mountain, and to some extent by elevation. We found evidence for gene flow between elevations. Results of genome-environment associations suggested similar selective pressures acting on both species, emanating mainly from precipitation and exposition rather than temperature.

CONCLUSIONS

Given their genetic structure and amount of gene flow among populations the two study species are suitable to serve as a model for genetic monitoring of climate change adaptation along an elevation gradient. Consequences of climate change will predominantly manifest via changes in precipitation and, thus, duration of snow cover in the snowbeds and indirectly via shrub encroachment accompanied by increasing shading of snowbeds at lower range margins. Assembling genomes of the study species and studying larger sample sizes and time series will be necessary to functionally characterize and validate the herein identified genomic loci putatively involved in adaptive processes.

Plant Species' Capacity for Range Shifts at the Habitat and Geographic Scales: A Trade-Off-Based Framework

Climate change is causing rapid shifts in the abiotic and biotic environmental conditions experienced by plant populations, but we lack generalizable frameworks for predicting the consequences for species. These changes may cause individuals to become poorly matched to their environments, potentially inducing shifts in the distributions of populations and altering species' habitat and geographic ranges. We present a trade-off-based framework for understanding and predicting whether plant species may undergo range shifts, based on ecological strategies defined by functional trait variation. We define a species' capacity for undergoing range shifts as the product of its colonization ability and the ability to express a phenotype well-suited to the environment across life stages (phenotype-environment matching), which are both strongly influenced by a species' ecological strategy and unavoidable trade-offs in function. While numerous strategies may be successful in an environment, severe phenotype-environment mismatches result in habitat filtering: propagules reach a site but cannot establish there. Operating within individuals and populations, these processes will affect species' habitat ranges at small scales, and aggregated across populations, will determine whether species track climatic changes and undergo geographic range shifts. This trade-off-based framework can provide a conceptual basis for species distribution models that are generalizable across plant species, aiding in the prediction of shifts in plant species' ranges in response to climate change.

Adaptive, maladaptive, neutral, or absent plasticity: Hidden caveats of reaction norms

Adaptive phenotypic plasticity may improve the response of individuals when faced with new environmental conditions. Typically, empirical evidence for plasticity is based on phenotypic reaction norms obtained in reciprocal transplant experiments. In such experiments, individuals from their native environment are transplanted into a different environment, and a number of trait values, potentially implicated in individuals' response to the new environment, are measured. However, the interpretations of reaction norms may differ depending on the nature of the assessed traits, which may not be known beforehand. For example, for traits that contribute to local adaptation, adaptive plasticity implies nonzero slopes of reaction norms. By contrast, for traits that are correlated to fitness, high tolerance to different environments (possibly due to adaptive plasticity in traits that contribute to adaptation) may, instead, result in flat reaction norms. Here we investigate reaction norms for adaptive versus fitness-correlated traits and how they may affect the conclusions regarding the contribution of plasticity. To this end, we first simulate range expansion along an environmental gradient where plasticity evolves to different values locally and then perform reciprocal transplant experiments in silico. We show that reaction norms alone cannot inform us whether the assessed trait exhibits locally adaptive, maladaptive, neutral, or no plasticity, without any additional knowledge of the traits assessed and species' biology. We use the insights from the model to analyse and interpret empirical data from reciprocal transplant experiments involving the marine isopod Idotea balthica sampled from two geographical locations with different salinities, concluding that the low-salinity population likely has reduced adaptive plasticity relative to the high-salinity population. Overall, we conclude that, when interpreting results from reciprocal transplant experiments, it is necessary to consider whether traits assessed are locally adaptive with respect to the environmental variable accounted for in the experiments or correlated to fitness.

Expanded distribution and predicted suitable habitat for the critically endangered yellow-tailed woolly monkey (Lagothrix flavicauda) in Perú

The Tropical Andes Biodiversity Hotspot holds a remarkable number of species at risk of extinction due to anthropogenic habitat loss, hunting, and climate change. One of these species, the critically endangered yellow-tailed woolly monkey (Lagothrix flavicauda), was recently observed in the region Junín, 206 km south of its previously known distribution. This range extension, combined with continued habitat loss, calls for a reevaluation of the species distribution, and available suitable habitat. Here, we present novel data from surveys at 53 sites in the regions of Junín, Cerro de Pasco, Ayacucho, and Cusco. We encountered L. flavicauda at 9 sites, all in Junín, and the congeneric Lagothrix lagotricha tschudii at 20 sites, but never in sympatry. Using these new localities along with all previous geographic localities for the species, we made predictive species distribution models based on ecological niche modeling using a generalized linear model and maximum entropy. Each model incorporated bioclimatic variables, forest cover, vegetation measurements, and elevation as predictor variables. The model evaluation showed >80% accuracy for all measures. Precipitation was the strongest predictor of species presence. Habitat suitability maps illustrate potential corridors for gene flow between the southern and northern populations, although much of this area is inhabited by L. l. tschudii whereas L. flavicauda has yet to be officially confirmed in these areas, by these or any other scientific surveys. An analysis of the current protected area (PA) network showed that ~75% of remaining suitable habitat is unprotected. With this, we suggest priority areas for new PAs or expansions to existing reserves that would conserve potential corridors between L. flavicauda populations. Further surveys and characterization of the distribution in intermediate areas, combined with studies on gene flow through these areas, are still needed to protect this species.

The Resilience of Polar Collembola (Springtails) in a Changing Climate

Assessing the resilience of polar biota to climate change is essential for predicting the effects of changing environmental conditions for ecosystems. Collembola are abundant in terrestrial polar ecosystems and are integral to food-webs and soil nutrient cycling. Using available literature, we consider resistance (genetic diversity; behavioural avoidance and physiological tolerances; biotic interactions) and recovery potential for polar Collembola. Polar Collembola have high levels of genetic diversity, considerable capacity for behavioural avoidance, wide thermal tolerance ranges, physiological plasticity, generalist-opportunistic feeding habits and broad ecological niches. The biggest threats to the ongoing resistance of polar Collembola are increasing levels of dispersal (gene flow), increased mean and extreme temperatures, drought, changing biotic interactions, and the arrival and spread of invasive species. If resistance capacities are insufficient, numerous studies have highlighted that while some species can recover from disturbances quickly, complete community-level recovery is exceedingly slow. Species dwelling deeper in the soil profile may be less able to resist climate change and may not recover in ecologically realistic timescales given the current rate of climate change. Ultimately, diverse communities are more likely to have species or populations that are able to resist or recover from disturbances. While much of the Arctic has comparatively high levels of diversity and phenotypic plasticity; areas of Antarctica have extremely low levels of diversity and are potentially much more vulnerable to climate change.

Population demography maintains biogeographic boundaries

Global biodiversity is organised into biogeographic regions that comprise distinct biotas. The contemporary factors maintaining differences in species composition between regions are poorly understood. Given evidence that populations with sufficient genetic variation can adapt to fill new habitats, it is surprising that more homogenisation of species assemblages across regions has not occurred. Theory suggests that expansion across biogeographic regions could be limited by reduced adaptive capacity due to demographic variation along environmental gradients, but this possibility has not been empirically explored. Using three independently curated data sets describing continental patterns of mammalian demography and population genetics, we show that populations near biogeographic boundaries have lower effective population sizes and genetic diversity, and are more genetically differentiated. These patterns are consistent with reduced adaptive capacity in areas where one biogeographic region transitions into the next. That these patterns are replicated across mammals suggests they are stable and generalisable in their contribution to long-term limits on biodiversity homogenisation. Understanding the contemporary processes that maintain compositional differences among regional biotas is crucial for our understanding of the current and future organisation of global biodiversity.

Increasing temporal variance leads to stable species range limits

What prevents populations of a species from adapting to the novel environments outside the species' geographic distribution? Previous models highlighted how gene flow across spatial environmental gradients determines species expansion versus extinction and the location of species range limits. However, space is only one of two axes of environmental variation-environments also vary in time, and we know temporal environmental variation has important consequences for population demography and evolution. We used analytical and individual-based evolutionary models to explore how temporal variation in environmental conditions influences the spread of populations across a spatial environmental gradient. We find that temporal variation greatly alters our predictions for range dynamics compared to temporally static environments. When temporal variance is equal across the landscape, the fate of species (expansion versus extinction) is determined by the interaction between the degree of temporal autocorrelation in environmental fluctuations and the steepness of the spatial environmental gradient. When the magnitude of temporal variance changes across the landscape, stable range limits form where this variance increases maladaptation sufficiently to prevent local persistence. These results illustrate the pivotal influence of temporal variation on the likelihood of populations colonizing novel habitats and the location of species range limits.

Effects of species interactions on the potential for evolution at species' range limits

Species’ ranges are limited by both ecological and evolutionary constraints. While there is a growing appreciation that ecological constraints include interactions among species, like competition, we know relatively little about how interactions contribute to evolutionary constraints at species' niche and range limits. Building on concepts from community ecology and evolutionary biology, we review how biotic interactions can influence adaptation at range limits by impeding the demographic conditions that facilitate evolution (which we term a ‘demographic pathway to adaptation’), and/or by imposing evolutionary trade-offs with the abiotic environment (a ‘trade-offs pathway’). While theory for the former is well-developed, theory for the trade-offs pathway is not, and empirical evidence is scarce for both. Therefore, we develop a model to illustrate how fitness trade-offs along biotic and abiotic gradients could affect the potential for range expansion and niche evolution following ecological release. The model shows that which genotypes are favoured at species' range edges can depend strongly on the biotic context and the nature of fitness trade-offs. Experiments that characterize trade-offs and properly account for biotic context are needed to predict which species will expand their niche or range in response to environmental change.

This article is part of the theme issue ‘Species’ ranges in the face of changing environments (Part II)’.

The response of a metapopulation to a changing environment

A species distributed across diverse environments may adapt to local conditions. We ask how quickly such a species changes its range in response to changed conditions. Szép et al. (Szép E, Sachdeva H, Barton NH. 2021 Polygenic local adaptation in metapopulations: a stochastic eco-evolutionary model. Evolution 75, 1030-1045 (doi:10.1111/evo.14210)) used the infinite island model to find the stationary distribution of allele frequencies and deme sizes. We extend this to find how a metapopulation responds to changes in carrying capacity, selection strength, or migration rate when deme sizes are fixed. We further develop a 'fixed-state' approximation. Under this approximation, polymorphism is only possible for a narrow range of habitat proportions when selection is weak compared to drift, but for a much wider range otherwise. When rates of selection or migration relative to drift change in a single deme of the metapopulation, the population takes a time of order m^-1 to reach the new equilibrium. However, even with many loci, there can be substantial fluctuations in net adaptation, because at each locus, alleles randomly get lost or fixed. Thus, in a finite metapopulation, variation may gradually be lost by chance, even if it would persist in an infinite metapopulation. When conditions change across the whole metapopulation, there can be rapid change, which is predicted well by the fixed-state approximation. This work helps towards an understanding of how metapopulations extend their range across diverse environments. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Environmental variation and biotic interactions limit adaptation at ecological margins: lessons from rainforest Drosophila and European butterflies

Models of local adaptation to spatially varying selection predict that maximum rates of evolution are determined by the interaction between increased adaptive potential owing to increased genetic variation, and the cost genetic variation brings by reducing population fitness. We discuss existing and new results from our laboratory assays and field transplants of rainforest Drosophila and UK butterflies along environmental gradients, which try to test these predictions in natural populations. Our data suggest that: (i) local adaptation along ecological gradients is not consistently observed in time and space, especially where biotic and abiotic interactions affect both gradient steepness and genetic variation in fitness; (ii) genetic variation in fitness observed in the laboratory is only sometimes visible to selection in the field, suggesting that demographic costs can remain high without increasing adaptive potential; and (iii) antagonistic interactions between species reduce local productivity, especially at ecological margins. Such antagonistic interactions steepen gradients and may increase the cost of adaptation by increasing its dimensionality. However, where biotic interactions do evolve, rapid range expansion can follow. Future research should test how the environmental sensitivity of genotypes determines their ecological exposure, and its effects on genetic variation in fitness, to predict the probability of evolutionary rescue at ecological margins. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Understanding the biology of species' ranges: when and how does evolution change the rules of ecological engagement?

Understanding processes that limit species' ranges has been a core issue in ecology and evolutionary biology for many decades, and has become increasingly important given the need to predict the responses of biological communities to rapid environmental change. However, we still have a poor understanding of evolution at range limits and its capacity to change the ecological 'rules of engagement' that define these communities, as well as the time frame over which this occurs. Here we link papers in the current volume to some key concepts involved in the interactions between evolutionary and ecological processes at species' margins. In particular, we separate hypotheses about species' margins that focus on hard evolutionary limits, which determine how genotypes interact with their environment, from those concerned with soft evolutionary limits, which determine where and when local adaptation can persist in space and time. We show how theoretical models and empirical studies highlight conditions under which gene flow can expand local limits as well as contain them. In doing so, we emphasize the complex interplay between selection, demography and population structure throughout a species' geographical and ecological range that determines its persistence in biological communities. However, despite some impressively detailed studies on range limits, particularly in invertebrates and plants, few generalizations have emerged that can predict evolutionary responses at ecological margins. We outline some directions for future work such as considering the impact of structural genetic variants and metapopulation structure on limits, and the interaction between range limits and the evolution of mating systems and non-random dispersal. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

A needle in the haystack? Applying species co-occurrence frameworks with fish assemblage data to identify species associations and sharpen ecological hypotheses

Different species can associate or interact in many ways, and methods exist for inferring associations and underlying mechanisms from incidence data (e.g., co-occurrence frameworks). These methods have received criticism despite their recent resurgence in the literature. However, co-occurrence frameworks for identifying nonrandomly associated species pairs (e.g., aggregated or segregated pairs) have value as heuristic tools for sharpening hypotheses concerning fish ecology. This paper provides a case study examining species co-occurrence across 33 stream fish assemblages in southeastern Oklahoma, USA, which were sampled twice (1974 and 2014). This study sought to determine (a) which species were nonrandomly associated, (b) what processes might have driven these associations and (c) how consistent patterns were across time. Associations among most pairs of species (24 species, 276 unique pairs) were not significantly different from random (>80%). Among all significant, nonrandomly associated species pairs (54 unique pairs), 78% (42 pairs) were aggregated and 22% (12 pairs) segregated. Most of these (28 pairs, 52%) were hypothesized to be driven by nonbiotic mechanisms: habitat filtering (20 pairs, 37%), dispersal limitation (two pairs, 0.4%) or both (six pairs, 11%). The remaining 26 nonrandomly associated pairs (48%) had no detectable signal of spatial or environmental factors involved with the association, therefore the potential for biotic interaction was not refuted. Only five species pairs were consistently associated across both sampling periods: stonerollers Campostoma spp. and orangebelly darter Etheostoma radiosum; red shiner Cyprinella lutrensis and bullhead minnow Pimephales vigilax; bluegill sunfish Lepomis macrochirus and redear sunfish Lepomis microlophus; redfin shiner Lythrurus umbratilis and bluntnose minnow Pimephales notatus; and bigeye shiner Notropis boops and golden shiner Notemigonus crysoleucas. Frameworks for identifying nonrandomly associated species pairs can provide insight into broader mechanisms of species assembly and point to potentially interesting species interactions (out of many possible pairs). However, this approach is best applied as a tool for sharpening hypotheses to be investigated further. Rather than a weakness, the heuristic nature is the strength of such methods, and can help guide biologists toward better questions by employing relatively cheap diversity survey data, which are often already in hand, to reduce complex interaction networks down to their nonstochastic parts which warrant further investigation.

Competition, Trait Variance Dynamics, and the Evolution of a Species' Range

Geographic ranges of communities of species evolve in response to environmental, ecological, and evolutionary forces. Understanding the effects of these forces on species' range dynamics is a major goal of spatial ecology. Previous mathematical models have jointly captured the dynamic changes in species' population distributions and the selective evolution of fitness-related phenotypic traits in the presence of an environmental gradient. These models inevitably include some unrealistic assumptions, and biologically reasonable ranges of values for their parameters are not easy to specify. As a result, simulations of the seminal models of this type can lead to markedly different conclusions about the behavior of such populations, including the possibility of maladaptation setting stable range boundaries. Here, we harmonize such results by developing and simulating a continuum model of range evolution in a community of species that interact competitively while diffusing over an environmental gradient. Our model extends existing models by incorporating both competition and freely changing intraspecific trait variance. Simulations of this model predict a spatial profile of species' trait variance that is consistent with experimental measurements available in the literature. Moreover, they reaffirm interspecific competition as an effective factor in limiting species' ranges, even when trait variance is not artificially constrained. These theoretical results can inform the design of, as yet rare, empirical studies to clarify the evolutionary causes of range stabilization.

Mate choice errors may contribute to slow spread of an invasive Eurasian longhorn beetle in North America

Tetropium fuscum (Coleoptera: Cerambycidae) is a Eurasian longhorn beetle and forest pest that first became invasive to Nova Scotia, Canada around 1990. In the time since its introduction, T. fuscum has spread only about 150 km from its point of introduction. In its invasive range, T. fuscum co-exists with its congener Tetropium cinnamopterum. Although they are ecologically similar species, T. fuscum tends to infest healthier trees and has a smaller host range than T. cinnamopterum. If they successfully interbreed, this could lead to hybrid individuals that are more problematic than either parent species. On the other hand, if T. fuscum can make mating errors in the field, but is not producing hybrid offspring, then this waste of mating resources could help explain the slow spread of T. fuscum in North America. We conducted no-choice and choice mating experiments between T. fuscum and T. cinnamopterum males and females and determined that both T. fuscum and T. cinnamopterum males make mate-choice errors with heterospecific females in a laboratory setting. Our results suggest that mating errors may play a role in the slow spread of T. fuscum in North America.

Coexistence patterns and diversity in a trait-based metacommunity on an environmental gradient

The dynamics of trait-based metacommunities have attracted much attention, but not much is known about how dispersal and spatial environmental variability mutually interact with each other to drive coexistence patterns and diversity. Here, we present a spatially explicit model of competition for two essential resources in a metacommunity on a one-dimensional environmental gradient. We find that both the strength of dispersal and the range of spatial environmental variability affect coexistence patterns, spatial structure, trait distribution, and local and regional diversity. Without dispersal, species are sorted according to their optimal growth conditions on the gradient. With the onset of dispersal, source-sink effects are initiated, which increases the effects of environmental filtering and interspecific competition and generates trait lumping, so that only a few species from an environment-defined trait range can survive. Interestingly, for very large dispersal rates, species distributions become spatially homogeneous, but nevertheless two species at the extreme ends of the trade-off curve can coexist for large environmental variability. Local species richness follows a classic hump-shaped dependence on dispersal rate, while local and regional diversity exhibit a pronounced peak for intermediate values of the environmental variability. Our findings provide important insights into the factors that shape the structure of trait-based metacommunities.

Of mutualism and migration: will interactions with novel ericoid mycorrhizal communities help or hinder northward Rhododendron range shifts?

Rapid climate change imperils many small-ranged endemic species as the climate envelopes of their native ranges shift poleward. In addition to abiotic changes, biotic interactions are expected to play a critical role in plant species' responses. Below-ground interactions are of particular interest given increasing evidence of microbial effects on plant performance and the prevalence of mycorrhizal mutualisms. We used greenhouse mesocosm experiments to investigate how natural northward migration/assisted colonization of Rhododendron catawbiense, a small-ranged endemic eastern U.S. shrub, might be influenced by novel below-ground biotic interactions from soils north of its native range, particularly with ericoid mycorrhizal fungi (ERM). We compared germination, leaf size, survival, and ERM colonization rates of endemic R. catawbiense and widespread R. maximum when sown on different soil inoculum treatments: a sterilized control; a non-ERM biotic control; ERM communities from northern R. maximum populations; and ERM communities collected from the native range of R. catawbiense. Germination rates for both species when inoculated with congeners' novel soils were significantly higher than when inoculated with conspecific soils, or non-mycorrhizal controls. Mortality rates were unaffected by treatment, suggesting that the unexpected reciprocal effect of each species' increased establishment in association with heterospecific ERM could have lasting demographic effects. Our results suggest that seedling establishment of R. catawbiense in northern regions outside its native range could be facilitated by the presence of extant congeners like R. maximum and their associated soil microbiota. These findings have direct relevance to the potential for successful poleward migration or future assisted colonization efforts.

Population structure and the influence of microenvironment and genetic similarity on individual growth at Alaskan white spruce treelines

Knowledge on the adaptation of trees to rapid environmental changes is essential to preserve forests and their ecosystem services under climate change. Treeline populations are particularly suitable for studying adaptation processes in trees, as environmental stress together with reduced gene flow can enhance local adaptation. We investigated white spruce (Picea glauca) populations in Alaska on one moisture-limited and two cold-limited treeline sites with a paired plot design of one forest and one treeline population each, resulting in six plots. Additionally, one forest plot in the middle of the distribution range complements the study design. We combined spatial, climatic and dendrochronological data with neutral genetic marker of 2203 trees to investigate population genetic structure and drivers of tree growth. We used several individual-based approaches including random slope mixed-effects models to test the influence of genetic similarity and microenvironment on growth performance. A high degree of genetic diversity was found within each of the seven plots associated with high rates of gene flow. We discovered a low genetic differentiation between the three sites which was better explained by geographic distances than by environmental differences, indicating genetic drift as the main driver of population differentiation. Our findings indicated that microenvironmental features had an overall larger influence on growth performances than genetic similarity among individuals. The effects of climate on growth differed between sites but were smaller than the effect of tree size. Overall, our results suggest that the high genetic diversity of white spruce may result in a wider range of phenotypes which enhances the efficiency of selection when the species is facing rapid climatic changes. In addition, the large intra-individual variability in growth responses may indicate the high phenotypic plasticity of white spruce which can buffer short-term environmental changes and, thus, allow enduring the present changing climate conditions.

High turn-over rates at the upper range limit and elevational source-sink dynamics in a widespread songbird

The formation of an upper distributional range limit for species breeding along mountain slopes is often based on environmental gradients resulting in changing demographic rates towards high elevations. However, we still lack an empirical understanding of how the interplay of demographic parameters forms the upper range limit in highly mobile species. Here, we study apparent survival and within-study area dispersal over a 700 m elevational gradient in barn swallows (Hirundo rustica) by using 15 years of capture-mark-recapture data. Annual apparent survival of adult breeding birds decreased while breeding dispersal probability of adult females, but not males increased towards the upper range limit. Individuals at high elevations dispersed to farms situated at elevations lower than would be expected by random dispersal. These results suggest higher turn-over rates of breeding individuals at high elevations, an elevational increase in immigration and thus, within-population source-sink dynamics between low and high elevations. The formation of the upper range limit therefore is based on preference for low-elevation breeding sites and immigration to high elevations. Thus, shifts of the upper range limit are not only affected by changes in the quality of high-elevation habitats but also by factors affecting the number of immigrants produced at low elevations.

An Evolutionary and Ecological Community Model for Distribution of Phenotypes and Abundances among Competing Species

AbstractHere, we propose a theory for the structure of communities of competing species. We include ecologically realistic assumptions, such as density dependence and stochastic fluctuations in the environment, and analyze how evolution caused by r- and K-selection will affect the packing of species in the phenotypic space as well as the species abundance distribution. Species-specific traits have the same matrix G of additive genetic variances and covariances, and evolution of mean traits is affected by fluctuations in population size of all species. In general, the model produces a shape of the distributions of log abundances that is skewed to the left, which is typical of most natural communities. Mean phenotypes of the species in the community are distributed approximately uniformly on the surface of a multidimensional sphere. However, environmental stochasticity generates selection that deviates species slightly from this surface; nonetheless, phenotypic distribution will be different from a random packing of species. This model of community evolution provides a theoretical framework that predicts a relationship between the structure of the phenotypic space and the form of species abundance distributions that can be compared against time series of variation in community structure.

Environmental association modelling with loci under divergent selection predicts the distribution range of a lizard

During the historical building of a species range, individual colonizers have to confront different ecological challenges, and the capacity of the species to broaden its range may depend on the total amount of adaptive genetic variation supplied by evolution. We set out to increase our understanding of what defines a range and the role of underlying genetics by trying to predict an entire species' range from the geographical distribution of its genetic diversity under selection. We sampled five populations of the western Mediterranean lizard Psammodromus algirus that inhabit a noticeable environmental gradient of temperature and precipitation. We correlated the genotypes of 95 individuals (18-20 individuals per population) for 21 SNPs putatively under selection with environmental scores on a bioclimatic gradient, using 1 × 1 km²  grid cells as sampling units. By extrapolating the resulting model to all possible combinations of alleles, we inferred all the geographic cells that were theoretically suitable for a given amount of genetic variance under selection. The inferred distribution range overlapped to a large extent with the realized range of the species (77.46% of overlap), including an accurate prediction of internal gaps and range borders. Our results suggest an adaptability threshold determined by the amount of genetic variation available that would be required to warrant adaptation beyond a certain limit of environmental variation. These results support the idea that the expansion of a species' range can be ultimately linked to the arising of new variants under selection (either newly selected variants from standing genetic variation or innovative mutations under selection).

Emergent Shapes of Trait-Based Competition Functions from Resource-Based Models: A Gaussian Is Not Normal in Plant Communities

AbstractIn community ecology, it is widely assumed that organisms with similar traits compete more intensely with one another for resources. This assumption is often encoded into theory and empirical tests via a unimodal competition function, which predicts that per capita competitive effect declines with separation in traits. Yet it remains unknown how well this function represents the true effect of traits on competitive outcomes, especially for long-lived plant communities, where lifetime fitness is difficult to estimate. Here, we evaluate the shape of competition functions embedded in two resource-based (RB) models, wherein plants compete for shared, essential resources. In the first RB model individuals compete for two essential nutrients, and in the second they compete for light in a size-based successional setting. We compared the shapes of the competition functions that emerged from interactions within these RB models to the unimodal function and others shapes commonly applied. In few instances did the trait-based competition function emerging from the RB model even vaguely resemble any of the shapes previously used. The mismatch between these two approaches suggests that theory derived using fixed competition functions based on trait separation may not apply well to plant systems, where individuals compete for shared resources. The more promising path will be to model depletion of resources by populations in relation to their traits, with its consequences for fitness landscapes and competitive exclusion.

Worldwide Genetic Structure Elucidates the Eurasian Origin and Invasion Pathways of Dothistroma septosporum, Causal Agent of Dothistroma Needle Blight

Dothistroma septosporum, the primary causal agent of Dothistroma needle blight, is one of the most significant foliar pathogens of pine worldwide. Its wide host and environmental ranges have led to its global success as a pathogen and severe economic damage to pine forests in many regions. This comprehensive global population study elucidated the historical migration pathways of the pathogen to reveal the Eurasian origin of the fungus. When over 3800 isolates were examined, three major population clusters were revealed: North America, Western Europe, and Eastern Europe, with distinct subclusters in the highly diverse Eastern European cluster. Modeling of historical scenarios using approximate Bayesian computation revealed the North American cluster was derived from an ancestral population in Eurasia. The Northeastern European subcluster was shown to be ancestral to all other European clusters and subclusters. The Turkish subcluster diverged first, followed by the Central European subcluster, then the Western European cluster, which has subsequently spread to much of the Southern Hemisphere. All clusters and subclusters contained both mating-types of the fungus, indicating the potential for sexual reproduction, although asexual reproduction remained the primary mode of reproduction. The study strongly suggests the native range of D. septosporum to be in Eastern Europe (i.e., the Baltic and Western Russia) and Western Asia.

Sharp boundary formation and invasion between spatially adjacent periodical cicada broods

Periodical cicadas, Magicicada spp., are a useful model system for understanding the population processes that influence range boundaries. Unlike most insects, these species typically exist at very high densities (occasionally >1000/ m²) and have unusually long life-spans (13 or 17 years). They spend most of their lives underground feeding on plant roots. After the underground period, adults emerge from the ground to mate and oviposit over a period of just a few days. Collections of populations that are developmentally synchronized across large areas are known as "broods". There are usually sharp boundaries between spatially adjacent broods and regions of brood overlap are generally small. The exact mechanism behind this developmental synchronization and the sharp boundary between broods remain unknown: previous studies have focused on the impacts of predator-driven Allee-effects, competition among nymphs, and their impacts on the persistence of off-synchronized emergence events. Here, we present a nonlinear Leslie-type matrix model to additionally consider cicada movement between spatially separated broods, and examine its role in maintaining brood boundaries and within-brood developmental synchrony that is seen in nature. We successfully identify ranges of competition and dispersal that lead to stable coexistence of broods that differ between spatial patches.

Resolving spatial complexities of hybridization in the context of the gray zone of speciation in North American ratsnakes (Pantherophis obsoletus complex)

Inferring the history of divergence between species in a framework that permits the presence of gene flow has been crucial for characterizing the "gray zone" of speciation, which is the period of time where lineages have diverged but have not yet achieved strict reproductive isolation. However, estimates of both divergence times and rates of gene flow often ignore spatial information, for example when considering the location and width of hybrid zones with respect to changes in the environment between lineages. Using population genomic data from the North American ratsnake complex (Pantherophis obsoletus), we connected phylogeographic estimates of lineage structure, migration, historical demography, and timing of divergence with hybrid zone dynamics. We examined the spatial context of diversification by linking migration and timing of divergence to the location and widths of hybrid zones. Artificial neural network approaches were applied to understand how landscape features and past climate have influenced population genetic structure among these lineages. We found that rates of migration between lineages were associated with the overall width of hybrid zones. Timing of divergence was not related to migration rate or hybrid zone width across species pairs but may be related to the number of alleles weakly introgressing through hybrid zones. This research underscores how incomplete reproductive isolation can be better understood by considering differential allelic introgression and the effects of historical and contemporary landscape features on the formation of lineages as well as overall genomic estimates of migration rates through time.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Mixed company: a framework for understanding the composition and organization of mixed-species animal groups

Mixed-species animal groups (MSGs) are widely acknowledged to increase predator avoidance and foraging efficiency, among other benefits, and thereby increase participants' fitness. Diversity in MSG composition ranges from two to 70 species of very similar or completely different phenotypes. Yet consistency in organization is also observable in that one or a few species usually have disproportionate importance for MSG formation and/or maintenance. We propose a two-dimensional framework for understanding this diversity and consistency, concentrating on the types of interactions possible between two individuals, usually of different species. One axis represents the similarity of benefit types traded between the individuals, while the second axis expresses asymmetry in the relative amount of benefits/costs accrued. Considering benefit types, one extreme represents the case of single-species groups wherein all individuals obtain the same supplementary, group-size-related benefits, and the other extreme comprises associations of very different, but complementary species (e.g. one partner creates access to food while the other provides vigilance). The relevance of social information and the matching of activities (e.g. speed of movement) are highest for relationships on the supplementary side of this axis, but so is competition; relationships between species will occur at points along this gradient where the benefits outweigh the costs. Considering benefit amounts given or received, extreme asymmetry occurs when one species is exclusively a benefit provider and the other a benefit user. Within this parameter space, some MSG systems are constrained to one kind of interaction, such as shoals of fish of similar species or leader-follower interactions in fish and other taxa. Other MSGs, such as terrestrial bird flocks, can simultaneously include a variety of supplementary and complementary interactions. We review the benefits that species obtain across the diversity of MSG types, and argue that the degree and nature of asymmetry between benefit providers and users should be measured and not just assumed. We then discuss evolutionary shifts in MSG types, focusing on drivers towards similarity in group composition, and selection on benefit providers to enhance the benefits they can receive from other species. Finally, we conclude by considering how individual and collective behaviour in MSGs may influence both the structure and processes of communities.

Time spent together and time spent apart affect song, feather colour and range overlap in tinkerbirds

Most studies on the processes driving evolutionary diversification highlight the importance of genetic drift in geographical isolation and natural selection across ecological gradients. Direct interactions among related species have received much less attention, but they can lead to character displacement, with recent research identifying patterns of displacement attributed to either ecological or reproductive processes. Together, these processes could explain complex, trait-specific patterns of diversification. Few studies, however, have examined the possible effects of these processes together or compared the divergence in multiple traits between interacting species among contact zones. Here, we show how traits of two Pogoniulus tinkerbird species vary among regions across sub-Saharan Africa. However, in addition to variation between regions consistent with divergence in refugial isolation, both song and morphology diverge between the species where they coexist. In West Africa, where the species are more similar in plumage, there is possible competitive or reproductive exclusion. In Central and East Africa, patterns of variation are consistent with agonistic character displacement. Molecular analyses support the hypothesis that differences in the age of interaction among regions can explain why species have evolved phenotypic differences and coexist in some regions but not others. Our findings suggest that competitive interactions between species and the time spent interacting, in addition to the time spent in refugial isolation, play important roles in explaining patterns of species diversification.

Maladaptation beyond a geographic range limit driven by antagonistic and mutualistic biotic interactions across an abiotic gradient

Species' geographic range limits often result from maladaptation to the novel environments beyond the range margin. However, we rarely know which aspects of the n-dimensional environment are driving this maladaptation. Especially of interest is the influence of abiotic versus biotic factors in delimiting species' distributions. We conducted a 2-year reciprocal transplant experiment involving manipulations of the biotic environment to explore how spatiotemporal gradients in precipitation, fatal mammalian herbivory, and pollination affected lifetime fitness within and beyond the range of the California annual plant, Clarkia xantiana ssp. xantiana. In the first, drier year of the experiment, fitness outside the range edge was limited mainly by low precipitation, and there was some evidence for local adaptation within the range. In the second, wetter year, we did not observe abiotic limitations to plant fitness outside the range; instead biotic interactions, especially herbivory, limited fitness outside the range. Together, protection from herbivory and supplementation of pollen resulted in three- to sevenfold increases in lifetime fitness outside the range margin in the abiotically benign year. Overall, our work demonstrates the importance of biotic interactions, particularly as they interact with the abiotic environment, in determining fitness beyond geographic range boundaries.

Maladapted Prey Subsidize Predators and Facilitate Range Expansion

Dispersal of prey from predator-free patches frequently supplies a trophic subsidy to predators by providing more prey than are produced locally. Prey arriving from predator-free patches might also have evolved weaker defenses against predators and thus enhance trophic subsidies by providing easily captured prey. Using local models assuming a linear or accelerating trade-off between defense and population growth rate, we demonstrate that immigration of undefended prey increased predator abundances and decreased defended prey through eco-evolutionary apparent competition. In individual-based models with spatial structure, explicit genetics, and gene flow along an environmental gradient, prey became maladapted to predators at the predator's range edge, and greater gene flow enhanced this maladaptation. The predator gained a subsidy from these easily captured prey, which enhanced its abundance, facilitated its persistence in marginal habitats, extended its range extent, and enhanced range shifts during environmental changes, such as climate change. Once the predator expanded, prey adapted to it and the advantage disappeared, resulting in an elastic predator range margin driven by eco-evolutionary dynamics. Overall, the results indicate a need to consider gene flow-induced maladaptation and species interactions as mutual forces that frequently determine ecological and evolutionary dynamics and patterns in nature.

Implications of existing local (mal)adaptations for ecological forecasting under environmental change

Standing genetic variation represents a genetic load on population fitness but can also support a rapid response to short-term environmental change, and the greatest potential source of such standing genetic variation typically exists among locally adapted populations living along an environmental gradient. Here, we develop a spatially explicit simulation model to quantify the contribution of existing genetic variation arising from migration-mutation-selection-drift balance to time to extinction under environmental change. Simulations reveal that local adaptation across a species range associated with an underlying environmental gradient could extend time to extinction by nearly threefold irrespective of the rate of environmental change. The potential for preadapted alleles to increase the rate of adaptation changes the relative importance of established extinction risk factors; in particular, it reduced the importance of the breadth of environmental tolerance and it increased the relative importance of fecundity. Although migration of preadapted alleles generally increased persistence time, it decreased it at rates of environmental change close to the critical rate of change by creating a population bottleneck, which ultimately limited the rate at which de novo mutations could arise. An analysis of the extinction dynamics further revealed that one consequence of gene flow is the potential to maximize population growth rate in at least part of the species range, which is likely to have consequences for forecasting the consequences of ecological interactions. Our study shows that predictions of persistence time change fundamentally when existing local adaptations are explicitly taken into account, underscoring the need to preserve and manage genetic diversity.

Contrasting forms of competition set elevational range limits of species

How abiotic and biotic factors constrain distribution limits at the harsh and benign edges of species ranges is hotly debated, partly because macroecological experiments testing the proximate causes of distribution limits are scarce. It has long been recognized - at least since Darwin's On the Origin of Species - that a harsh climate strengthens competition and thus sets species range limits. Using thorough field manipulations along a large elevation gradient, we show the mechanisms by which temperature determines competition type, resulting in a transition from interference to exploitative competition from the lower to the upper elevation limits in burying beetles (Nicrophorus nepalensis). This transition is an example of Darwin's classic hypothesis that benign climates favor direct competition for highly accessible resources while harsh climates result in competition through resources of high rivalry. We propose that identifying the properties of these key resources will provide a more predictive framework to understand the interplay between biotic and abiotic factors in determining geographic range limits.

Fight or flight? Potential tradeoffs between drought defense and reproduction in conifers

Plants frequently exhibit tradeoffs between reproduction and growth when resources are limited, and often change these allocation patterns in response to stress. Shorter-lived plants such as annuals tend to allocate relatively more resources toward reproduction when stressed, while longer-lived plants tend to invest more heavily in survival and stress defense. However, severe stress may affect the fitness implications of allocating relatively more resources to reproduction versus stress defense. Increased drought intensity and duration have led to widespread mortality events in coniferous forests. In this review, we ask how potential tradeoffs between reproduction and survival influence the likelihood of drought-induced mortality and species persistence. We propose that trees may exhibit what we call 'fight or flight' behaviors under stress. 'Fight' behaviors involve greater resource allocation toward survival (e.g., growth, drought-resistant xylem and pest defense). 'Flight' consists of higher relative allocation of resources to reproduction, potentially increasing both offspring production and mortality risk for the adult. We hypothesize that flight behaviors increase as drought stress escalates the likelihood of mortality in a given location.

Understanding the margin squeeze: Differentiation in fitness-related traits between central and trailing edge populations of Corallina officinalis

ABSTRACT

Assessing population responses to climate-related environmental change is key to understanding the adaptive potential of the species as a whole. Coralline algae are critical components of marine shallow water ecosystems where they function as important ecosystem engineers. Populations of the calcifying algae Corallina officinalis from the center (southern UK) and periphery (northern Spain) of the North Atlantic species natural distribution were selected to test for functional differentiation in thermal stress response. Physiological measurements of calcification, photosynthesis, respiration, growth rates, oxygen, and calcification evolution curves were performed using closed cell respirometry methods. Species identity was genetically confirmed via DNA barcoding. Through a common garden approach, we identified distinct vulnerability to thermal stress of central and peripheral populations. Southern populations showed a decrease in photosynthetic rate under environmental conditions of central locations, and central populations showed a decline in calcification rates under southern conditions. This shows that the two processes of calcification and photosynthesis are not as tightly coupled as previously assumed. How the species as whole will react to future climatic changes will be determined by the interplay of local environmental conditions and these distinct population adaptive traits.

OPEN RESEARCH BADGES

This article has earned an Open Materials Badge for making publicly available the components of the research methodology needed to reproduce the reported procedure and analysis. All materials are available at https://doi.pangaea.de/10.1594/PANGAEA.899568.

Growth-dependent drug susceptibility can prevent or enhance spatial expansion of a bacterial population

As a population wave expands, organisms at the tip typically experience plentiful nutrients while those behind the front become nutrient-depleted. If the environment also contains a gradient of some inhibitor (e.g. a toxic drug), a tradeoff exists: the nutrient-rich tip is more exposed to the inhibitor, while the nutrient-starved region behind the front is less exposed. Here we show that this can lead to complex dynamics when the organism's response to the inhibitory substance is coupled to nutrient availability. We model a bacterial population which expands in a spatial gradient of antibiotic, under conditions where either fast-growing bacteria at the wave's tip, or slow-growing, resource-limited bacteria behind the front are more susceptible to the antibiotic. We find that growth-rate dependent susceptibility can have strong effects on the dynamics of the expanding population. If slow-growing bacteria are more susceptible, the population wave advances far into the inhibitory zone, leaving a trail of dead bacteria in its wake. In contrast, if fast-growing bacteria are more susceptible, the wave is blocked at a much lower concentration of antibiotic, but a large population of live bacteria remains behind the front. Our results may contribute to understanding the efficacy of different antimicrobials for spatially structured microbial populations such as biofilms, as well as the dynamics of ecological population expansions more generally.