Allele frequencies and selection coefficients in locally adapted populations.

Understanding the role of natural selection in driving evolutionary change requires accurate estimates of the strength of selection acting at the genetic level in the wild. This is challenging to achieve but may be easier in the case of populations in migration-selection balance. When two populations are at equilibrium under migration-selection balance, there exist loci whose alleles are selected different ways in the two populations. Such loci can be identified from genome sequencing by their high values of FST. This raises the question of what is the strength of selection on locally-adaptive alleles. To answer this question we analyse a 1-locus 2-allele model of a population distributed between two niches. We show by simulation of selected cases that the outputs from finite-population models are essentially the same as those from deterministic infinite-population models. We then derive theory for the infinite-population model showing the dependence of selection coefficients on equilibrium allele frequencies, migration rates, dominance and relative population sizes in the two niches. An Excel spreadsheet is provided for the calculation of selection coefficients and their approximate standard errors from observed values of population parameters. We illustrate our results with a worked example, with graphs showing the dependence of selection coefficients on equilibrium allele frequencies, and graphs showing how FST depends on the selection coefficients acting on the alleles at a locus. Given the extent of recent progress in ecological genomics, we hope our methods may help those studying migration-selection balance to quantify the advantages conferred by adaptive genes.

Sexual imprinting leads to speciation in locally adapted populations.

Sexual imprinting is widespread in birds and other species but its existence requires explanation. Our results suggest that sexual imprinting leads to speciation in locally-adapted populations if a neutral mating cue-e.g., novel plumage coloration-arises through mutation. Importantly, the mating cue locus is not linked to adaptation loci. Local adaptation is a necessary precursor to speciation and occurs when evolution results in stable genetic polymorphisms with one allele predominating in some areas while others predominate elsewhere. Here we use a deterministic two-niche population genetic model to map the set of migration and selection rates for which polymorphic evolutionary outcomes, i.e., local adaptations, can occur. Approximate equations for the boundaries of the set of polymorphic evolutionary outcomes were derived by Bulmer (American Naturalist, 106, 254, 1972), but our results, obtained by deterministic simulation of the evolutionary process, show that one of Bulmer's equations is inaccurate except when the level of dominance is 0.5, and fails if one of the alleles is dominant. Having an accurate map of the set of migration and selection rates for which polymorphic evolutionary outcomes can occur, we then show using the model of Sibly et al. (Ecology and Evolution, 9, 13506, 2019) that local adaptation in all analyzed cases leads to speciation if a new neutral mating cue arises by mutation. We finish by considering how genome sequencing makes possible testing our model and its predictions.

Mycorrhizal type of woody plants influences understory species richness in British broadleaved woodlands.

Mature temperate woodlands are commonly dominated by ectomycorrhizal trees, whereas understory plants predominantly form arbuscular mycorrhizal associations. Due to differences in plant-fungus compatibility between canopy and ground layer vegetation the 'mycorrhizal mediation hypothesis' predicts that herbaceous plant establishment may be limited by a lack of suitable mycorrhizal fungal inoculum. We examined plant species data for 103 woodlands across Great Britain recorded in 1971 and in 2000 to test whether herbaceous plant species richness was related to the proportion of arbuscular mycorrhizal woody plants. We compared the effect of mycorrhizal type with other important drivers of woodland plant species richness. We found a positive effect of the relative abundance of arbuscular mycorrhizal woody plants on herbaceous plant species richness. The size of the observed effect was smaller than that of pH. Moreover, the effect persisted over time, despite many woodlands undergoing marked successional change and increased understorey shading. This work supports the mycorrhizal mediation hypothesis in British woodlands and suggests that increased abundance of arbuscular mycorrhizal woody plants is associated with greater understory plant species richness.

Factors affecting fisher decisions: The case of the inshore fishery for European sea bass (Dicentrarchus labrax).

Fishery management relies on forecasts of fish abundance over time and space, on scales of months and kilometres. While much research has focussed on the drivers of fish populations, there has been less investigation of the decisions made day-to-day by fishers and their subsequent impact on fishing pressure. Studies that focus on the fisher decisions of smaller vessels may be particularly important due to the prevalence of smaller vessels in many fisheries and their potential vulnerability to bad weather and economic change. Here we outline a methodology with which to identify the factors affecting fisher decisions and success as well as quantifying their effects. We analyse first the decision of when to leave port, and then the success of the fishing trip. Fisher behaviour is here analysed in terms of the decisions taken by fishers in response to bio-physical and socio-economic changes and to illustrate our method, we describe its application to the under 10-meter fleet targeting sea bass in the UK. We document the effects of wave height and show with increasing wave height fewer vessels left port to go fishing. The decision to leave port was only substantially affected by time of high tide at one of the four ports investigated. We measured the success of fishing trips by the landings of sea bass (kg) per metre of vessel length. Fishing success was lower when wave height was greater and when fish price had increased relative to the previous trip. Fuel price was unimportant, but a large proportion of the variation in success was explained by variation between individual vessels, presumably due to variation in skipper ability or technical restrictions due to vessel characteristics. The results are discussed in the context of management of sea bass and other small-scale inshore fisheries.

SEASIM-NEAM: A Spatially-Explicit Agent-based SIMulator of North East Atlantic Mackerel population dynamics.

In 2018 we published a spatially-explicit individual-based model (IBM) that uses satellite-derived maps of food availability and temperature to predict Northeast Atlantic mackerel (Scomber scombrus, NEAM) population dynamics. Since then, to address various ecological questions, we have extended the IBM to include additional processes and data. Throughout its development, technical documents have been provided in the form of e.g. supplementary information to published articles. However, we acknowledge that it would be difficult for potential users to collate information from separate supplementary documents and gain a full understanding of the current state of the IBM. Here, we provide a full technical specification of the latest version of our IBM. The technical specification is provided in the standard ODD (Overview, Design concepts and Details) format, and supplemented by a TRACE (TRAnsparent and Comprehensive model Evaludation) document. For the first time, we give our model the acronym SEASIM-NEAM: a Spatially-Explicit Agent-based SIMulator of North East Atlantic Mackerel population dynamics. This article supersedes previous documentation. Going forward we hope that this article will stimulate development of similar models.•This article collates improvements that have been made to SEASIM-NEAM over time.

Behavior underpins the predictive power of a trait-based model of butterfly movement.

Dispersal ability is key to species persistence in times of environmental change. Assessing a species' vulnerability and response to anthropogenic changes is often performed using one of two methods: correlative approaches that infer dispersal potential based on traits, such as wingspan or an index of mobility derived from expert opinion, or a mechanistic modeling approach that extrapolates displacement rates from empirical data on short-term movements.Here, we compare and evaluate the success of the correlative and mechanistic approaches using a mechanistic random-walk model of butterfly movement that incorporates relationships between wingspan and sex-specific movement behaviors.The model was parameterized with new data collected on four species of butterfly in the south of England, and we observe how wingspan relates to flight speeds, turning angles, flight durations, and displacement rates.We show that flight speeds and turning angles correlate with wingspan but that to achieve good prediction of displacement even over 10 min the model must also include details of sex- and species-specific movement behaviors.We discuss what factors are likely to differentially motivate the sexes and how these could be included in mechanistic models of dispersal to improve their use in ecological forecasting.

Multiple environmental controls explain global patterns in soil animal communities.

Soil animals play important roles in ecosystem functioning and stability, but the environmental controls on their communities are not fully understood. In this study, we compiled a dataset of soil animal communities for which the abundance and body mass of multiple soil animal groups were recorded. The mass-abundance scaling relationships were then used to investigate multiple environmental controls on soil animal community composition. The data reveal latitudinal shifts from high abundances of small soil animals at high latitudes to greater relative abundances of large soil animals at low latitudes. A hierarchical linear mixed effects model was applied to reveal the environmental variables shaping these latitudinal trends. The final hierarchical model identified mean annual temperature, soil pH and soil organic carbon content as key environmental controls explaining global mass-abundance scaling relationships in soil animal communities (R2c = 0.828, Ngroup = 117). Such relationships between soil biota with climate and edaphic conditions have been previously identified for soil microbial, but not soil animal, communities at a global scale. More comprehensive global soil community datasets are needed to better understand the generality of these relationships over a broader range of global ecosystems and soil animal groups.

The importance of including habitat-specific behaviour in models of butterfly movement.

Dispersal is a key process affecting population persistence and major factors affecting dispersal rates are the amounts, connectedness and properties of habitats in landscapes. We present new data on the butterfly Maniola jurtina in flower-rich and flower-poor habitats that demonstrates how movement and behaviour differ between sexes and habitat types, and how this effects consequent dispersal rates. Females had higher flight speeds than males, but their total time in flight was four times less. The effect of habitat type was strong for both sexes, flight speeds were ~ 2.5 × and ~ 1.7 × faster on resource-poor habitats for males and females, respectively, and flights were approximately 50% longer. With few exceptions females oviposited in the mown grass habitat, likely because growing grass offers better food for emerging caterpillars, but they foraged in the resource-rich habitat. It seems that females faced a trade-off between ovipositing without foraging in the mown grass or foraging without ovipositing where flowers were abundant. We show that taking account of habitat-dependent differences in activity, here categorised as flight or non-flight, is crucial to obtaining good fits of an individual-based model to observed movement. An important implication of this finding is that incorporating habitat-specific activity budgets is likely necessary for predicting longer-term dispersal in heterogeneous habitats, as habitat-specific behaviour substantially influences the mean (> 30% difference) and kurtosis (1.4 × difference) of dispersal kernels. The presented IBMs provide a simple method to explicitly incorporate known activity and movement rates when predicting dispersal in changing and heterogeneous landscapes.

How phenotypic matching based on neutral mating cues enables speciation in locally adapted populations.

Maynard Smith's (American Naturalist, 1966, 100, 637) suggestion that in some cases a prerequisite for speciation is the existence of local ecological adaptations has not received much attention to date. Here, we test the hypothesis using a model like that of Maynard Smith but differing in the way animals disperse between niches. In previous studies, males disperse randomly between niches but females stay put in their natal niche. As a first step toward generalizing the model, we here analyze the case that equal proportions of the two sexes disperse between niches before breeding. Supporting Maynard Smith's (1966) hypothesis, we find that once local adaptations are established, a neutral mating cue at an independent locus can rapidly enable speciation in populations with a suitable mechanism for phenotype matching. We find that stable ecological polymorphisms are relatively insensitive to the strength of selection, but depend crucially on the extent of dispersal between niches, with a threshold of ~5% if population sizes in two niches are equal. At higher levels of dispersal, ecological differentiation is lost. These results contrast with those of earlier studies and shed light on why parapatric speciation is limited by the extent of gene flow. Our testable model provides a candidate explanation for the rapid speciation rates, diversity of appearance and occurrence of "species flocks" observed among some African cichlids and neotropical birds and may also have implications for the occurrence of punctuational change on phylogenies.

Data on the movement behaviour of four species of grassland butterfly.

This Data in Brief article describes data on the movement behaviour of four species of grassland butterflies collected over three years and at four sites in southern England. The datasets consist of the movement tracks of Maniola jurtina, Aricia agestis, Pyronia tithonus, and Melanargia galathea, recorded using standard methods and presented as steps distances and turning angles. Sites consisted of nectar-rich field margins, meadows, and mown short turf grasslands with minimal flowers. In total, 783 unique movement tracks were collected. The data were used for analysing the movement behaviour of the species and for parameterising individual-based movement models.

Integrating the influence of weather into mechanistic models of butterfly movement.

BACKGROUND: Understanding the factors influencing movement is essential to forecasting species persistence in a changing environment. Movement is often studied using mechanistic models, extrapolating short-term observations of individuals to longer-term predictions, but the role of weather variables such as air temperature and solar radiation, key determinants of ectotherm activity, are generally neglected. We aim to show how the effects of weather can be incorporated into individual-based models of butterfly movement thus allowing analysis of their effects. METHODS: We constructed a mechanistic movement model and calibrated it with high precision movement data on a widely studied species of butterfly, the meadow brown (Maniola jurtina), collected over a 21-week period at four sites in southern England. Day time temperatures during the study ranged from 14.5 to 31.5 °C and solar radiation from heavy cloud to bright sunshine. The effects of weather are integrated into the individual-based model through weather-dependent scaling of parametric distributions representing key behaviours: the durations of flight and periods of inactivity. RESULTS: Flight speed was unaffected by weather, time between successive flights increased as solar radiation decreased, and flight duration showed a unimodal response to air temperature that peaked between approximately 23 °C and 26 °C. After validation, the model demonstrated that weather alone can produce a more than two-fold difference in predicted weekly displacement. CONCLUSIONS: Individual Based models provide a useful framework for integrating the effect of weather into movement models. By including weather effects we are able to explain a two-fold difference in movement rate of M. jurtina consistent with inter-annual variation in dispersal measured in population studies. Climate change for the studied populations is expected to decrease activity and dispersal rates since these butterflies already operate close to their thermal optimum.

The shark-tuna dichotomy: why tuna lay tiny eggs but sharks produce large offspring.

Teleosts such as tunas and billfish lay millions of tiny eggs weighing on the order of 0.001 g, whereas chondrichthyes such as sharks and rays produce a few eggs or live offspring weighing about 2% of adult body mass, as much as 10 000 g in some species. Why are the strategies so extreme, and why are intermediate ones absent? Building on previous work, we show quantitatively how offspring size reflects the relationship between growth and death rates. We construct fitness contours as functions of offspring size and number, and show how these can be derived from juvenile growth and survivorship curves. Convex contours, corresponding to Pearl Type 1 and 2 survivorship curves, select for extremes, either miniscule or large offspring; concave contours select for offspring of intermediate size. Of particular interest are what we call critical straight-line fitness contours, corresponding to log-linear Pearl Type 3 survivorship curves, which separate regimes that select for opposite optimal offspring sizes.

The influence of soil communities on the temperature sensitivity of soil respiration.

Soil respiration represents a major carbon flux between terrestrial ecosystems and the atmosphere, and is expected to accelerate under climate warming. Despite its importance in climate change forecasts, however, our understanding of the effects of temperature on soil respiration (RS) is incomplete. Using a metabolic ecology approach we link soil biota metabolism, community composition and heterotrophic activity to predict RS rates across five biomes. We find that accounting for the ecological mechanisms underpinning decomposition processes predicts climatological RS variations observed in an independent dataset (n = 312). The importance of community composition is evident because without it RS is substantially underestimated. With increasing temperature, we predict a latitudinal increase in RS temperature sensitivity, with Q10 values ranging between 2.33 ± 0.01 in tropical forests to 2.72 ± 0.03 in tundra. This global trend has been widely observed, but has not previously been linked to soil communities.

Equal fitness paradigm explained by a trade-off between generation time and energy production rate.

Most plant, animal and microbial species of widely varying body size and lifestyle are nearly equally fit as evidenced by their coexistence and persistence through millions of years. All organisms compete for a limited supply of organic chemical energy, derived mostly from photosynthesis, to invest in the two components of fitness: survival and production. All organisms are mortal because molecular and cellular damage accumulates over the lifetime; life persists only because parents produce offspring. We call this the equal fitness paradigm. The equal fitness paradigm occurs because: (1) there is a trade-off between generation time and productive power, which have equal-but-opposite scalings with body size and temperature; smaller and warmer organisms have shorter lifespans but produce biomass at higher rates than larger and colder organisms; (2) the energy content of biomass is essentially constant, ~22.4 kJ g-1 dry body weight; and (3) the fraction of biomass production incorporated into surviving offspring is also roughly constant, ~10-50%. As organisms transmit approximately the same quantity of energy per gram to offspring in the next generation, no species has an inherent lasting advantage in the struggle for existence. The equal fitness paradigm emphasizes the central importance of energy, biological scaling relations and power-time trade-offs in life history, ecology and evolution.

Taking error into account when fitting models using Approximate Bayesian Computation.

Stochastic computer simulations are often the only practical way of answering questions relating to ecological management. However, due to their complexity, such models are difficult to calibrate and evaluate. Approximate Bayesian Computation (ABC) offers an increasingly popular approach to this problem, widely applied across a variety of fields. However, ensuring the accuracy of ABC's estimates has been difficult. Here, we obtain more accurate estimates by incorporating estimation of error into the ABC protocol. We show how this can be done where the data consist of repeated measures of the same quantity and errors may be assumed to be normally distributed and independent. We then derive the correct acceptance probabilities for a probabilistic ABC algorithm, and update the coverage test with which accuracy is assessed. We apply this method, which we call error-calibrated ABC, to a toy example and a realistic 14-parameter simulation model of earthworms that is used in environmental risk assessment. A comparison with exact methods and the diagnostic coverage test show that our approach improves estimation of parameter values and their credible intervals for both models.

Genetic polymorphisms between altruism and selfishness close to the Hamilton threshold rb = c.

Genes that in certain conditions make their carriers altruistic are being identified, and altruism and selfishness have shown to be heritable in man. This raises the possibility that genetic polymorphisms for altruism/selfishness exist in man and other animals. Here we characterize some of the conditions in which genetic polymorphisms may occur. We show for dominant or recessive alleles how the positions of stable equilibria depend on the benefit to the recipient, b, and the cost to the altruist, c, for diploid altruists helping half or full sibs, and haplodiploid altruists helping sisters. Stable polymorphisms always occur close to the Hamilton threshold rb = c. The position of the stable equilibrium moves away 0 or 1 with both increases in c, the cost paid by the altruist, and increasing divergence from the Hamilton threshold, and alleles for selfishness can reach frequencies around 50%. We evaluate quantitative estimates of b, c and r from field studies in the light of these predictions, but the values do not fall in the regions where genetic polymorphisms are expected. Nevertheless, it will be interesting to see as genes for altruism are discovered whether they are accompanied by alternate alleles for selfishness.

Fundamental insights into ontogenetic growth from theory and fish.

The fundamental features of growth may be universal, because growth trajectories of most animals are very similar, but a unified mechanistic theory of growth remains elusive. Still needed is a synthetic explanation for how and why growth rates vary as body size changes, both within individuals over their ontogeny and between populations and species over their evolution. Here, we use Bertalanffy growth equations to characterize growth of ray-finned fishes in terms of two parameters, the growth rate coefficient, K, and final body mass, m∞. We derive two alternative empirically testable hypotheses and test them by analyzing data from FishBase. Across 576 species, which vary in size at maturity by almost nine orders of magnitude, K scaled as [Formula: see text]. This supports our first hypothesis that growth rate scales as [Formula: see text] as predicted by metabolic scaling theory; it implies that species that grow to larger mature sizes grow faster as juveniles. Within fish species, however, K scaled as [Formula: see text]. This supports our second hypothesis, which predicts that growth rate scales as [Formula: see text] when all juveniles grow at the same rate. The unexpected disparity between across- and within-species scaling challenges existing theoretical interpretations. We suggest that the similar ontogenetic programs of closely related populations constrain growth to [Formula: see text] scaling, but as species diverge over evolutionary time they evolve the near-optimal [Formula: see text] scaling predicted by metabolic scaling theory. Our findings have important practical implications because fish supply essential protein in human diets, and sustainable yields from wild harvests and aquaculture depend on growth rates.

Metabolic theory predicts whole-ecosystem properties.

Understanding the effects of individual organisms on material cycles and energy fluxes within ecosystems is central to predicting the impacts of human-caused changes on climate, land use, and biodiversity. Here we present a theory that integrates metabolic (organism-based bottom-up) and systems (ecosystem-based top-down) approaches to characterize how the metabolism of individuals affects the flows and stores of materials and energy in ecosystems. The theory predicts how the average residence time of carbon molecules, total system throughflow (TST), and amount of recycling vary with the body size and temperature of the organisms and with trophic organization. We evaluate the theory by comparing theoretical predictions with outputs of numerical models designed to simulate diverse ecosystem types and with empirical data for real ecosystems. Although residence times within different ecosystems vary by orders of magnitude-from weeks in warm pelagic oceans with minute phytoplankton producers to centuries in cold forests with large tree producers-as predicted, all ecosystems fall along a single line: residence time increases linearly with slope = 1.0 with the ratio of whole-ecosystem biomass to primary productivity (B/P). TST was affected predominantly by primary productivity and recycling by the transfer of energy from microbial decomposers to animal consumers. The theory provides a robust basis for estimating the flux and storage of energy, carbon, and other materials in terrestrial, marine, and freshwater ecosystems and for quantifying the roles of different kinds of organisms and environments at scales from local ecosystems to the biosphere.

Patterns of maximum body size evolution in Cenozoic land mammals: eco-evolutionary processes and abiotic forcing.

There is accumulating evidence that macroevolutionary patterns of mammal evolution during the Cenozoic follow similar trajectories on different continents. This would suggest that such patterns are strongly determined by global abiotic factors, such as climate, or by basic eco-evolutionary processes such as filling of niches by specialization. The similarity of pattern would be expected to extend to the history of individual clades. Here, we investigate the temporal distribution of maximum size observed within individual orders globally and on separate continents. While the maximum size of individual orders of large land mammals show differences and comprise several families, the times at which orders reach their maximum size over time show strong congruence, peaking in the Middle Eocene, the Oligocene and the Plio-Pleistocene. The Eocene peak occurs when global temperature and land mammal diversity are high and is best explained as a result of niche expansion rather than abiotic forcing. Since the Eocene, there is a significant correlation between maximum size frequency and global temperature proxy. The Oligocene peak is not statistically significant and may in part be due to sampling issues. The peak in the Plio-Pleistocene occurs when global temperature and land mammal diversity are low, it is statistically the most robust one and it is best explained by global cooling. We conclude that the macroevolutionary patterns observed are a result of the interplay between eco-evolutionary processes and abiotic forcing.

How body mass and lifestyle affect juvenile biomass production in placental mammals.

In mammals, the mass-specific rate of biomass production during gestation and lactation, here called maternal productivity, has been shown to vary with body size and lifestyle. Metabolic theory predicts that post-weaning growth of offspring, here termed juvenile productivity, should be higher than maternal productivity, and juveniles of smaller species should be more productive than those of larger species. Furthermore because juveniles generally have similar lifestyles to their mothers, across species juvenile and maternal productivities should be correlated. We evaluated these predictions with data from 270 species of placental mammals in 14 taxonomic/lifestyle groups. All three predictions were supported. Lagomorphs, perissodactyls and artiodactyls were very productive both as juveniles and as mothers as expected from the abundance and reliability of their foods. Primates and bats were unproductive as juveniles and as mothers, as expected as an indirect consequence of their low predation risk and consequent low mortality. Our results point the way to a mechanistic explanation for the suite of correlated life-history traits that has been called the slow-fast continuum.