Measuring Nonlinear Selection

Inversions Can Accumulate Balanced Sexual Antagonism: Evidence from Simulations and Drosophila Experiments

Chromosomal inversion polymorphisms can be common, but the causes of their persistence are often unclear. We propose a model for the maintenance of inversion polymorphism, which requires that some variants contribute antagonistically to two phenotypes, one of which has negative frequency-dependent fitness. These conditions yield a form of frequency-dependent disruptive selection, favoring two predominant haplotypes segregating alleles that favor opposing antagonistic phenotypes. An inversion associated with one haplotype can reduce the fitness load incurred by generating recombinant offspring, reinforcing its linkage to the haplotype and enabling both haplotypes to accumulate more antagonistic variants than expected otherwise. We develop and apply a forward simulator to examine these dynamics under a tradeoff between survival and male display. These simulations indeed generate inversion-associated haplotypes with opposing sex-specific fitness effects. Antagonism strengthens with time, and can ultimately yield karyotypes at surprisingly predictable frequencies, with striking genotype frequency differences between sexes and between developmental stages. To test whether this model may contribute to well-studied yet enigmatic inversion polymorphisms in Drosophila melanogaster, we track inversion frequencies in laboratory crosses to test whether they influence male reproductive success or survival. We find that two of the four tested inversions show significant evidence for the tradeoff examined, with In(3R)K favoring survival and In(3L)Ok favoring male reproduction. In line with the apparent sex-specific fitness effects implied for both of those inversions, In(3L)Ok was also found to be less costly to the viability and/or longevity of males than females, whereas In(3R)K was more beneficial to female survival. Based on this work, we expect that balancing selection on antagonistically pleiotropic traits may provide a significant and underappreciated contribution to the maintenance of natural inversion polymorphism.

Beauty or function? The opposing effects of natural and sexual selection on cuticular hydrocarbons in male black field crickets

Although many theoretical models of male sexual trait evolution assume that sexual selection is countered by natural selection, direct empirical tests of this assumption are relatively uncommon. Cuticular hydrocarbons (CHCs) are known to play an important role not only in restricting evaporative water loss but also in sexual signalling in most terrestrial arthropods. Insects adjusting their CHC layer for optimal desiccation resistance is often thought to come at the expense of successful sexual attraction, suggesting that natural and sexual selection are in opposition for this trait. In this study, we sampled the CHCs of male black field crickets (Teleogryllus commodus) using solid-phase microextraction and then either measured their evaporative water loss or mating success. We then used multivariate selection analysis to quantify the strength and form of natural and sexual selection targeting male CHCs. Both natural and sexual selection imposed significant linear and stabilizing selection on male CHCs, although for very different combinations. Natural selection largely favoured an increase in the total abundance of CHCs, especially those with a longer chain length. In contrast, mating success peaked at a lower total abundance of CHCs and declined as CHC abundance increased. However, mating success did improve with an increase in a number of specific CHC components that also increased evaporative water loss. Importantly, this resulted in the combination of male CHCs favoured by natural selection and sexual selection being strongly opposing. Our findings suggest that the balance between natural and sexual selection is likely to play an important role in the evolution of male CHCs in T. commodus and may help explain why CHCs are so divergent across populations and species.

Strong selection is poorly aligned with genetic variation in Ipomoea hederacea: implications for divergence and constraint

The multivariate evolution of populations is the result of the interactions between natural selection, drift, and the underlying genetic structure of the traits involved. Covariances among traits bias responses to selection, and the multivariate axis which describes the greatest genetic variation is expected to be aligned with patterns of divergence across populations. An exception to this expectation is when selection acts on trait combinations lacking genetic variance, which limits evolutionary change. Here we used a common garden field experiment of individuals from 57 populations of Ipomoea hederacea to characterize linear and nonlinear selection on 5 quantitative traits in the field. We then formally compare patterns of selection to previous estimates of within population genetic covariance structure (the G-matrix) and population divergence in these traits. We found that selection is poorly aligned with previous estimates of genetic covariance structure and population divergence. In addition, the trait combinations favored by selection were generally lacking genetic variation, possessing approximately 15%-30% as much genetic variation as the most variable combination of traits. Our results suggest that patterns of population divergence are likely the result of the interplay between adaptive responses, correlated responses, and selection favoring traits lacking genetic variation.

Using landscape genomics to delineate future adaptive potential for climate change in the Yosemite toad (Anaxyrus canorus)

An essential goal in conservation biology is delineating population units that maximize the probability of species persisting into the future and adapting to future environmental change. However, future-facing conservation concerns are often addressed using retrospective patterns that could be irrelevant. We recommend a novel landscape genomics framework for delineating future "Geminate Evolutionary Units" (GEUs) in a focal species: (1) identify loci under environmental selection, (2) model and map adaptive conservation units that may spawn future lineages, (3) forecast relative selection pressures on each future lineage, and (4) estimate their fitness and likelihood of persistence using geo-genomic simulations. Using this process, we delineated conservation units for the Yosemite toad (Anaxyrus canorus), a U.S. federally threatened species that is highly vulnerable to climate change. We used a genome-wide dataset, redundancy analysis, and Bayesian association methods to identify 24 candidate loci responding to climatic selection (R 2 ranging from 0.09 to 0.52), after controlling for demographic structure. Candidate loci included genes such as MAP3K5, involved in cellular response to environmental change. We then forecasted future genomic response to climate change using the multivariate machine learning algorithm Gradient Forests. Based on all available evidence, we found three GEUs in Yosemite National Park, reflecting contrasting adaptive optima: YF-North (high winter snowpack with moderate summer rainfall), YF-East (low to moderate snowpack with high summer rainfall), and YF-Low-Elevation (low snowpack and rainfall). Simulations under the RCP 8.5 climate change scenario suggest that the species will decline by 29% over 90 years, but the highly diverse YF-East lineage will be least impacted for two reasons: (1) geographically it will be sheltered from the largest climatic selection pressures, and (2) its standing genetic diversity will promote a faster adaptive response. Our approach provides a comprehensive strategy for protecting imperiled non-model species with genomic data alone and has wide applicability to other declining species.

An analytical pipeline to support robust research on the ecology, evolution, and function of floral volatiles

Research on floral volatiles has grown substantially in the last 20 years, which has generated insights into their diversity and prevalence. These studies have paved the way for new research that explores the evolutionary origins and ecological consequences of different types of variation in floral scent, including community-level, functional, and environmentally induced variation. However, to address these types of questions, novel approaches are needed that can handle large sample sizes, provide quality control measures, and make volatile research more transparent and accessible, particularly for scientists without prior experience in this field. Drawing upon a literature review and our own experiences, we present a set of best practices for next-generation research in floral scent. We outline methods for data collection (experimental designs, methods for conducting field collections, analytical chemistry, compound identification) and data analysis (statistical analysis, database integration) that will facilitate the generation and interpretation of quality data. For the intermediate step of data processing, we created the R package bouquet, which provides a data analysis pipeline. The package contains functions that enable users to convert chromatographic peak integrations to a filtered data table that can be used in subsequent statistical analyses. This package includes default settings for filtering out non-floral compounds, including background contamination, based on our best-practice guidelines, but functions and workflows can be easily customized as necessary. Next-generation research into the ecology and evolution of floral scent has the potential to generate broadly relevant insights into how complex traits evolve, their genomic architecture, and their consequences for ecological interactions. In order to fulfill this potential, the methodology of floral scent studies needs to become more transparent and reproducible. By outlining best practices throughout the lifecycle of a project, from experimental design to statistical analysis, and providing an R package that standardizes the data processing pipeline, we provide a resource for new and seasoned researchers in this field and in adjacent fields, where high-throughput and multi-dimensional datasets are common.

Immediate predation risk alters the relationship between potential and realised selection on male traits in the Trinidad guppy Poecilia reticulata

Imminent predation risk affects mating behaviours in prey individuals in a multitude of ways that can theoretically impact the strength of sexual selection, as well as its operation on traits. However, empirical studies of the effects of imminent predation risk on sexual selection dynamics are still scarce. Here we explore how perceived predation affects: (1) the relationship between the opportunity for selection and the actual strength of selection on male traits; and (2) which traits contribute to male fitness and the shape of selection on these traits. We simulate two consecutive reproductive episodes, under control conditions and perceived predation risk using experimental populations of Trinidad guppies. The opportunity for selection is higher under predation risk compared to the control condition, but realised selection on traits remains unaffected. Pre- and postcopulatory traits follow complex patterns of nonlinear selection in both conditions. Differences in selection gradients deviate from predictions based on evolutionary and non-lethal effects of predation, the most notable being strong disruptive selection on courtship rate under predation risk. Our results demonstrate that sexual selection is sensitive to imminent predation risk perception and reinforce the notion that both trait-based and variance-based metrics should be employed for an informative quantification.

Density-dependent patterns of multivariate selection on sperm motility and morphology in a broadcast spawning mussel

Sperm cells exhibit extraordinary phenotypic variation, both among taxa and within individual species, yet our understanding of the adaptive value of sperm trait variation across multiple contexts is incomplete. For species without the opportunity to choose mating partners, such as sessile broadcast spawning invertebrates, fertilization depends on gamete interactions, which in turn can be strongly influenced by local environmental conditions that alter the concentration of sperm and eggs. However, the way in which such environmental factors impact phenotypic selection on functional gamete traits remains unclear in most systems. Here, we analyze patterns of linear and nonlinear multivariate selection under experimentally altered local sperm densities (densities within the capture zone of eggs) on a range of functionally important sperm traits in the broadcast spawning marine mussel, Mytilus galloprovincialis. Specifically, we assay components of sperm motility and morphology across two fertilization environments that simulate either sperm limitation (when there are too few sperm to fertilize all available eggs), or sperm saturation (when there are many more sperm than required for fertilization, and the risk of polyspermy and embryonic failure is heightened). Our findings reveal that the strength, form, and targets of selection on sperm depend on the prevailing fertilization environment. In particular, our analyses revealed multiple significant axes of nonlinear selection on sperm motility traits under sperm limitation, but only significant negative directional selection on flagellum length under sperm saturation. These findings highlight the importance of local sperm densities in driving the adaptation of sperm phenotypes, particularly those related to sperm motility, in broadcast spawning invertebrates.

Morphological Evolution: Bioinspired Methods for Analyzing Bioinspired Robots

To fully understand the evolution of complex morphologies, analyses cannot stop at selection: It is essential to investigate the roles and interactions of multiple processes that drive evolutionary outcomes. The challenges of undertaking such analyses have affected both evolutionary biologists and evolutionary roboticists, with their common interests in complex morphologies. In this paper, we present analytical techniques from evolutionary biology, selection gradient analysis and morphospace walks, and we demonstrate their applicability to robot morphologies in analyses of three evolutionary mechanisms: randomness (genetic mutation), development (an explicitly implemented genotype-to-phenotype map), and selection. In particular, we applied these analytical techniques to evolved populations of simulated biorobots—embodied robots designed specifically as models of biological systems, for the testing of biological hypotheses—and we present a variety of results, including analyses that do all of the following: illuminate different evolutionary dynamics for different classes of morphological traits; illustrate how the traits targeted by selection can vary based on the likelihood of random genetic mutation; demonstrate that selection on two selected sets of morphological traits only partially explains the variance in fitness in our biorobots; and suggest that biases in developmental processes could partially explain evolutionary dynamics of morphology. When combined, the complementary analytical approaches discussed in this paper can enable insight into evolutionary processes beyond selection and thereby deepen our understanding of the evolution of robotic morphologies.

Detecting patterns of correlational selection with sampling error: A simulation study

The adoption of a multivariate perspective of selection implies the existence of multivariate adaptive peaks and pervasive correlational selection that promotes co-adaptation between traits. However, to test for the ubiquity of correlational selection in nature, we must first have a sense of how well can we estimate multivariate nonlinear selection (i.e., the γ-matrix) in the face of sampling error. To explore the sampling properties of estimated γ-matrices, we simulated inidividual traits and fitness under a wide range of sample sizes, using different strengths of correlational selection and of stabilizing selection, combined with different number of traits under selection, different amounts of residual variance in fitness, and distinct patterns of selection. We then ran nonlinear regressions with these simulated datasets to simulate γ-matrices after adding random error to individual fitness. To test how well could we detect the imposed pattern of correlational selection at different sample sizes, we measured the similarity between simulated and imposed γ-matrices. We show that detection of the pattern of correlational selection is highly dependent on the total strength of selection on traits and on the amount of residual variance in fitness. Minimum sample size needs to be at least 500 to precisely estimate the pattern of correlational selection. Furthermore, a pattern of selection in which different sets of traits contribute to different functions is the easiest to diagnose, even when using a large number of traits (10 traits), but with sample sizes in the order of 1000 individuals. Consequently, we recommend working with sets of traits from distinct functional complexes and fitness proxies less prone to effects of environmental and demographic stochasticity to test for correlational selection with lower sample sizes.

Shell colour diversification induced by ecological release: A shift in natural selection after a migration event

Ecological release is often attributed to the rapid adaptive diversification of phenotypic traits. However, it is not well understood how natural selection changes its strength and direction through the process of ecological release. Herein, we demonstrated how shell colour of the Japanese land snail Euhadra peliomphala simodae has diversified via a shift in natural selection due to ecological release after migration from the mainland to an island. This snail's shell colour diversified on the island due to disruptive selection after migration from the mainland. We used trail camera traps to identify the cause of natural selection on both the mainland and the island. We then conducted a mark-recapture experiment while collecting microhabitat use data. In total, we captured and marked around 1,700 snails on the mainland, some of which were preyed upon by an unknown predator. The trail camera traps showed that the predator is the large Japanese field mouse Apodemus speciosus, and the predatory frequency was higher on the mainland than on the island. However, this predation did not correlate with shell colour. Microhabitat use on the island was more extensive than on the mainland, with snails on the island using both ground and arboreal microhabitats. A Bayesian estimation showed that the stabilizing selection on shell colour came from factors other than predation. Our results suggest that the course of natural selection was modified due to ecological release after migration from the mainland, explaining one cause of the phenotypic diversification.

Ecological causes of fluctuating natural selection on habitat choice in an amphibian

We estimated natural selection targeting three traits related to habitat choice in a frog (Pseudacris maculata) breeding in pools on the rocky shores of Isle Royale, Michigan, over 16 years. Our aim was to identify the form and ecological causes of annual variation in directional and correlational selection as expressed in the survival and growth of tadpoles. We found directional selection favoring early breeding, but pool choice was under weak stabilizing selection. However, the form of stabilizing selection and the position of the optimum trait value shifted among years with the severity of disturbance and the intensity of biotic interactions. In years when wave wash and pool desiccation were severe, selection shifted to favor tadpoles in habitats where these risks were less pronounced. If predatory dragonfly larvae were abundant, selection favored tadpoles in small pools where dragonflies did not occur. When intraspecific competition was strong, selection favored early broods within a broader range of pool types. The agents of selection in this study—biotic interactions and disturbance—are common to many ecological systems and frequently exhibit temporal variation; this suggests that fluctuating selection may be widespread in natural populations.

Correlational selection in the age of genomics

Ecologists and evolutionary biologists are well aware that natural and sexual selection do not operate on traits in isolation, but instead act on combinations of traits. This long-recognized and pervasive phenomenon is known as multivariate selection, or-in the particular case where it favours correlations between interacting traits-correlational selection. Despite broad acknowledgement of correlational selection, the relevant theory has often been overlooked in genomic research. Here, we discuss theory and empirical findings from ecological, quantitative genetic and genomic research, linking key insights from different fields. Correlational selection can operate on both discrete trait combinations and quantitative characters, with profound implications for genomic architecture, linkage, pleiotropy, evolvability, modularity, phenotypic integration and phenotypic plasticity. We synthesize current knowledge and discuss promising research approaches that will enable us to understand how correlational selection shapes genomic architecture, thereby linking quantitative genetic approaches with emerging genomic methods. We suggest that research on correlational selection has great potential to integrate multiple fields in evolutionary biology, including developmental and functional biology, ecology, quantitative genetics, phenotypic polymorphisms, hybrid zones and speciation processes.

Multiple modes of selection can influence the role of phenotypic plasticity in species' invasions: Evidence from a manipulative field experiment

In exploring the roles of phenotypic plasticity in the establishment and early evolution of invading species, little empirical attention has been given to the importance of correlational selection acting upon suites of functionally related plastic traits in nature. We illustrate how this lack of attention has limited our ability to evaluate plasticity's role during invasion and also, the costs and benefits of plasticity. We addressed these issues by transplanting clones of European-derived Plantago lanceolata L. genotypes into two temporally variable habitats in the species' introduced range in North America. Phenotypic selection analyses were performed for each habitat to estimate linear, quadratic, and correlational selection on phenotypic trait values and plasticities in the reproductive traits: flowering onset and spike and scape lengths. Also, we measured pairwise genetic correlations for our "colonists." Results showed that (a) correlational selection acted on trait plasticity after transplantation, (b) selection favored certain combinations of genetically correlated and uncorrelated trait values and plasticities, and (c) using signed, instead of absolute, values of plasticity in analyses facilitated the detection of correlational selection on trait value-plasticity combinations and their adaptive value. Based on our results, we urge future studies on species invasions to (a) measure correlational selection and (b) retain signed values of plasticity in order to better discriminate between adaptive and maladaptive plasticity.

Sexual selection maintains a female-specific character in a species with dynamic sex roles

The effects of sexual selection are more conspicuous among male animals, and, as a result, the majority of sexual selection research focuses on males. However, burgeoning evidence suggests that sexual selection also acts on females, and there have been calls for an increased focus on females. Here, we used a multivariate approach to analyze sexual selection in Kawanaphila nartee, a spermatophore gift-giving bushcricket with dynamic sex roles. Early in the breeding season, females compete for males, and, later, when environmental food resources are more abundant, sex roles revert to Darwinian convention. Ear size, which is much greater in females than in males, has been suggested to affect female fitness as females with larger ears are more likely to reach calling males first under sex-role-reversed conditions. We tested this suggestion and found evidence of positive linear and nonlinear correlational selection acting on female ear size early in the breeding season (under reversed sex roles) but not later in the breeding season (under Darwinian sex roles). Interestingly, there was no correlation between mating success and reproductive success (Bateman gradient) at any time during the season. Together, our results indicate that even brief and circumscribed periods of intrasexual competition among females can lead to sexual selection on morphological characters and that this selection may not depend on multiple mating. Considering the wealth of reports in the literature of brief episodes of intrasexual competition among female animals, we recommend increased study of sexual selection acting on females.

Sexual selection on the genital lobes of male Drosophila simulans

Sexual selection is thought to be responsible for the rapid divergent evolution of male genitalia with several studies detecting multivariate sexual selection on genital form. However, in most cases, selection is only estimated during a single episode of selection, which provides an incomplete view of net selection on genital traits. Here, we estimate the strength and form of multivariate selection on the genitalia arch of Drosophila simulans when mating occurs in the absence of a competitor and during sperm competition, in both sperm defence and offense roles (i.e., when mating first and last). We found that the strength of sexual selection on the genital arch was strongest during noncompetitive mating and weakest during sperm offense. However, the direction of selection was similar across selection episodes with no evidence for antagonistic selection. Overall, selection was not particularly strong despite genitals clearly evolving rapidly in this species.

Predator presence and recent climatic warming raise body temperatures of island lizards

In ectothermic predator-prey relationships, evasion of predation by prey depends on physiological and behavioural responses relating to the thermal biology of both predator and prey. On Japan's Izu Islands, we investigated a prey lizard's physiological and thermal responses to the presence of a snake predator over geologic time in addition to recent climatic warming. Foraging lizard body temperatures increased by 1.3 °C from 1981 to 2019 overall, yet were 2.9 °C warmer on snake islands relative to snake-free islands. We also detected snake predator-induced selection on hind leg length, which in turn is a major determinant for sprint speed only in lizard populations exposed to predation by snakes. Accordingly, we found that warmer prey body temperatures result in faster sprint speeds by the prey at temperatures suboptimal for the snake predator, and therefore contribute to escaping predation. Given recent climatic change, further warming could irrevocably alter this and other ectothermic predator-prey relationships.

Trait differences among discrete morphs of a color polymorphic lizard, Podarcis erhardii

Color polymorphism defies evolutionary expectations as striking phenotypic variation is maintained within a single species. Color and other traits mediate social interactions, and stable polymorphism within a population is hypothesized to be related to correlational selection of other phenotypic traits among color morphs. Here, we report on a previously unknown throat color polymorphism in the Aegean Wall Lizard (Podarcis erhardii) and examine morph-correlated differences in traits important to social behavior and communication: maximum bite force capacity and chemical signal profile. We find that both sexes of P. erhardii have three color morphs: orange, yellow, and white. Moreover, orange males are significantly larger and tend to bite harder than yellow and white males. Although the established color polymorphism only partially matches the observed intraspecific variation in chemical signal signatures, the chemical profile of the secretions of orange males is significantly divergent from that of white males. Our findings suggest that morph colors are related to differences in traits that are crucial for social interactions and competitive ability, illustrating the need to look beyond color when studying polymorphism evolution.

Harsh environments promote alloparental care across human societies

Alloparental care is central to human life history, which integrates exceptionally short interbirth intervals and large birth size with an extended period of juvenile dependency and increased longevity. Formal models, previous comparative research, and palaeoanthropological evidence suggest that humans evolved higher levels of cooperative childcare in response to increasingly harsh environments. Although this hypothesis remains difficult to test directly, the relative importance of alloparental care varies across human societies, providing an opportunity to assess how local social and ecological factors influence the expression of this behaviour. We therefore, investigated associations between alloparental infant care and socioecology across 141 non-industrialized societies. We predicted increased alloparental care in harsher environments, due to the fitness benefits of cooperation in response to shared ecological challenges. We also predicted that starvation would decrease alloparental care, due to prohibitive energetic costs. Using Bayesian phylogenetic multilevel models, we tested these predictions while accounting for potential confounds as well as for population history. Consistent with our hypotheses, we found increased alloparental infant care in regions characterized by both reduced climate predictability and relatively lower average temperatures and precipitation. We also observed reduced alloparental care under conditions of high starvation. These results provide evidence of plasticity in human alloparenting in response to ecological contexts, comparable to previously observed patterns across avian and mammalian cooperative breeders. This suggests convergent social evolutionary processes may underlie both inter- and intraspecific variation in alloparental care.

Combined effects of temperature and macronutrient balance on life-history traits in Drosophila melanogaster: implications for life-history trade-offs and fundamental niche

Temperature and nutrition are amongst the most influential environmental determinants of Darwinian fitness in ectotherms. Since the ongoing climate warming is known to alter nutritional environments encountered by ectotherms, a precise understanding of the integrated effects of these two factors on ectotherm performance is essential for improving the accuracy of predictions regarding how ectotherms will respond to climate warming. Here we employed response surface methodology to examine how multiple life-history traits were expressed across a grid of environmental conditions representing full combinations of six ambient temperatures (13, 18, 23, 28, 31, 33 °C) and eight dietary protein:carbohydrate ratios (P:C = 1:16, 1:8, 1:4, 1:2, 1:1, 2:1, 4:1, 8:1) in Drosophila melanogaster. Different life-history traits were maximized in different regions in the two-dimensional temperature-nutrient space. The optimal temperature and P:C ratio identified for adult lifespan (13 °C and 1:16) were lower than those for early-life female fecundity (28 °C and 4:1). Similar divergence in thermal and nutritional optima was found between body mass at adult emergence (18 °C and P:C 1:1) and the rate of pre-adult development (28 °C and P:C 4:1). Pre-adult survival was maximized over a broad range of temperature (18-28 °C) and P:C ratio (1:8-8:1). These results indicate that the occurrence of life-history trade-offs is regulated by both temperature and dietary P:C ratio. The estimated measure of fitness was maximized at 23 °C and P:C 2:1. Based on the shape of the response surface constructed for this estimated fitness, we characterized the fundamental thermal and nutritional niche for D. melanogaster with unprecedented detail.

Patterns of environmental variance across environments and traits in domestic cattle

The variance in phenotypic trait values is a product of environmental and genetic variation. The sensitivity of traits to environmental variation has a genetic component and is likely to be under selection. However, there are few studies investigating the evolution of this sensitivity, in part due to the challenges of estimating the environmental variance. The livestock literature provides a wealth of studies that accurately partition components of phenotypic variance, including the environmental variance, in well-defined environments. These studies involve breeds that have been under strong selection on mean phenotype in optimal environments for many generations, and therefore represent an opportunity to study the potential evolution of trait sensitivity to environmental conditions. Here, we use literature on domestic cattle to examine the evolution of micro-environmental variance (CVR-the coefficient of residual variance) by testing for differences in expression of CVR in animals from the same breed reared in different environments. Traits that have been under strong selection did not follow a null expectation of an increase in CVR in heterogenous environments (e.g., grazing), a pattern that may reflect evolution of increased uniformity in heterogeneous environments. When comparing CVR across environments of different levels of optimality, here measured by trait mean, we found a reduction in CVR in the more optimal environments for both life history and growth traits. Selection aimed at increasing trait means in livestock breeds typically occurs in the more optimal environments, and we therefore suspect that the decreased CVR is a consequence of evolution of the expression of micro-environmental variance in this environment. Our results highlight the heterogeneity in micro-environmental variance across environments and point to possible connections to the intensity of selection on trait means.

Why does allometry evolve so slowly?

Morphological allometry is striking due to its evolutionary conservatism, making it an example of a certain sort of evolutionary stasis. Organisms that vary in size, whether for developmental, environmental, or evolutionary reasons, adopt shapes that are predictable from that size alone. There are two major hypotheses to explain this. It may be that natural selection strongly favors each allometric pattern, or that organisms lack the development and genetic capacity to produce variant shapes for selection to act on. Using a high-throughput system for measuring the size and shape of Drosophila wings, we documented an allometric pattern that has been virtually unchanged for 40 million years. We performed an artificial selection experiment on the static allometric slope within one species. In just 26 generations, we were able to increase the slope from 1.1 to 1.4, and decrease it to 0.8. Once artificial selection was suspended, the slope rapidly evolved back to a value near the initial static slope. This result decisively rules out the hypothesis that allometry is preserved due to a lack of genetic variation, and provides evidence that natural selection acts to maintain allometric relationships. On the other hand, it seems implausible that selection on allometry in the wing alone could be sufficiently strong to maintain static allometries over millions of years. This suggests that a potential explanation for stasis is selection on a potentially large number of pleiotropic effects. This seems likely in the case of allometry, as the sizes of all parts of the body may be altered when the allometric slope of one body part is changed. Unfortunately, hypotheses about pleiotropy have been very difficult to test. We lay out an approach to begin the systematic study of pleiotropic effects using genetic manipulations and high-throughput phenotyping.

Egg-induced changes to sperm phenotypes shape patterns of multivariate selection on ejaculates

Sperm cells exhibit extraordinary phenotypic diversity and rapid rates of evolution, yet the adaptive value of most sperm traits remains equivocal. Recent findings suggest that to understand how selection targets ejaculates, we must recognize that female-imposed physiological conditions often alter sperm phenotypes. These phenotypic changes may influence the relationships among sperm traits and their association with fitness. Here, we show that chemical substances released by eggs (known to modify sperm physiology and behaviour) alter patterns of selection on a suite of sperm traits in the mussel Mytilus galloprovincialis. We use multivariate selection analyses to characterize linear and nonlinear selection acting on sperm traits in (a) seawater alone and (b) seawater containing egg-derived chemicals (egg water). Our analyses revealed that nonlinear selection on canonical axes of multiple traits (notably sperm velocity, sperm linearity and percentage of motile sperm) was the most important form of selection overall, but importantly these patterns were only evident when sperm phenotypes were measured in egg water. These findings reveal the subtle way that females can alter patterns of selection, with the implication that overlooking environmentally moderated changes to sperm, may result in erroneous interpretations of how selection targets phenotypic (co)variation in sperm traits.

Multivariate stabilizing sexual selection and the evolution of male and female genital morphology in the red flour beetle

Male genitals are highly divergent in animals with internal fertilization. Most studies attempting to explain this diversity have focused on testing the major hypotheses of genital evolution (the lock‐and‐key, pleiotropy, and sexual selection hypotheses), and quantifying the form of selection targeting male genitals has played an important role in this endeavor. However, we currently know far less about selection targeting female genitals or how male and female genitals interact during mating. Here, we use formal selection analysis to show that genital size and shape is subject to strong multivariate stabilizing sexual selection in both sexes of the red flour beetle, Tribolium castaneum. Moreover, we show significant sexual selection on the covariance between the sexes for specific aspects of genital shape suggesting that male and female genitalia also interact to determine the successful transfer of a spermatophore during mating. Our work therefore highlights the important role that both male and female genital morphologies play in determining mating success and that these effects can occur independently, as well as through their interaction. Moreover, it cautions against the overly simplistic view that the sexual selection targeting genital morphology will always be directional in form and restricted primarily to males.

Stabilizing selection and adaptive evolution in a combination of two traits in an arctic ungulate

Stabilizing selection is thought to be common in wild populations and act as one of the main evolutionary mechanisms, which constrain phenotypic variation. When multiple traits interact to create a combined phenotype, correlational selection may be an important process driving adaptive evolution. Here, we report on phenotypic selection and evolutionary changes in two natal traits in a semidomestic population of reindeer (Rangifer tarandus) in northern Finland. The population has been closely monitored since 1969, and detailed data have been collected on individuals since they were born. Over the length of the study period (1969-2015), we found directional and stabilizing selection toward a combination of earlier birth date and heavier birth mass with an intermediate optimum along the major axis of the selection surface. In addition, we demonstrate significant changes in mean traits toward earlier birth date and heavier birth mass, with corresponding genetic changes in breeding values during the study period. Our results demonstrate evolutionary changes in a combination of two traits, which agree closely with estimated patterns of phenotypic selection. Knowledge of the selective surface for combinations of genetically correlated traits are vital to predict how population mean phenotypes and fitness are affected when environments change.

Status anxiety mediates the positive relationship between income inequality and sexualization

Income inequality generates and amplifies incentives, particularly incentives for individuals to elevate or maintain their status, with important consequences for the individuals involved and aggregate outcomes for their societies [R. G. Wilkinson, K. E. Pickett, Annu. Rev. Sociol. 35, 493-511 (2009)]. Economically unequal environments intensify men's competition for status, respect, and, ultimately, mating opportunities, thus elevating aggregate rates of violent crime and homicide [M. Daly, M. Wilson, Evolutionary Psychology and Motivation (2001)]. Recent evidence shows that women are more likely to post "sexy selfies" on social media and that they spend more on beautification in places where inequality is high rather than low [K. R. Blake, B. Bastian, T. F. Denson, et al., Proc. Natl. Acad. Sci. U.S.A. 115, 8722-8727 (2018)]. Here we test experimentally for causal links between income inequality and individual self-sexualization and status-related competition. We show that manipulating income inequality in a role-playing task indirectly increases women's intentions to wear revealing clothing and that it does so by increasing women's anxiety about their place in the social hierarchy. The effects are not better accounted for by wealth/poverty than by inequality or by modeling anxiety about same-sex competitors in place of status anxiety. The results indicate that women's appearance enhancement is partly driven by status-related goals.

Disruptive selection of shell colour in land snails: a mark–recapture study of Euhadra peliomphala simodae

Many theoretical studies have suggested that disruptive selection plays an important role in phenotypic divergence, but few studies have determined the action of disruptive selection on phenotypic divergence via field studies. This study investigated the effect of disruptive selection on shell colour polymorphism in the Japanese land snail Euhadra peliomphala simodae to determine whether extreme phenotypes of snail shell colour are favoured over intermediate phenotypes. We conducted field surveys on an oceanic island with black, yellow and intermediate-coloured E. p. simodae snails. We captured and marked ~1800 individual snails and monitored their survival over 18 months. We quantified shell colours against images and examined the frequency distribution of shell colour variation. The variation exhibited a bimodal distribution with a far lower frequency of intermediate-coloured snails than of black or yellow snails. The population sizes of the three snail groups fluctuated synchronously with the changing seasons. Bayesian estimates showed lower survival rates for juvenile intermediate-coloured snails than for juvenile black and yellow snails, implying there was disruptive selection associated with shell colour. We suggest this disruptive selection may have resulted in the evolutionary divergence of the snail’s shell colour within the lineage having high shell colour variation.

Comparative analysis of the multivariate genetic architecture of morphological traits in three species of Gomphocerine grasshoppers

Evolutionary change is the change in trait values across generations, and usually occurs in multidimensional trait space rather than along isolated traits. Genetic covariation influences the magnitude and direction of evolutionary change and can be statistically summarized by the additive genetic (co)variance matrix, G. While G can affect the response to selection, it is exposed to evolutionary change by selection and genetic drift, but the magnitude and speed of these changes are poorly understood. We use comparative G matrix analyses to assess evolution of the shape and orientation of G over longer timescales in three species of Gomphocerine grasshoppers. We estimate 10 × 10 G matrices for five morphological traits expressed in both sexes. We find low-to-moderate heritabilities (average 0.36), mostly large cross-sex correlations (average 0.54) and moderate between-trait correlations (average 0.34). G matrices differ significantly among species with wing length contributing most to these differences. Wing length is the trait that is most divergent among species, suggesting it has been under selection during species divergence. The more distantly related species, Pseudochorthippus parallelus, was the most different in the shape of G. Projection of contemporary genetic variation into the divergence space D illustrates that the major axis of genetic variation in Gomphocerippus rufus is aligned with divergence from Chorthippus biguttulus, while the major axis of genetic variation in neither of the species is aligned with the divergence between Pseudochorthippus parallelus and the other two species. Our results demonstrate significant differences in G matrices with a phylogenetic signal in the differentiation.

An evolutionary quantitative genetics model for phenotypic (co)variances under limited dispersal, with an application to socially synergistic traits

Darwinian evolution consists of the gradual transformation of heritable traits due to natural selection and the input of random variation by mutation. Here, we use a quantitative genetics approach to investigate the coevolution of multiple quantitative traits under selection, mutation, and limited dispersal. We track the dynamics of trait means and of variance-covariances between traits that experience frequency-dependent selection. Assuming a multivariate-normal trait distribution, we recover classical dynamics of quantitative genetics, as well as stability and evolutionary branching conditions of invasion analyses, except that due to limited dispersal, selection depends on indirect fitness effects and relatedness. In particular, correlational selection that associates different traits within-individuals depends on the fitness effects of such associations between-individuals. We find that these kin selection effects can be as relevant as pleiotropy for the evolution of correlation between traits. We illustrate this with an example of the coevolution of two social traits whose association within-individuals is costly but synergistically beneficial between-individuals. As dispersal becomes limited and relatedness increases, associations between-traits between-individuals become increasingly targeted by correlational selection. Consequently, the trait distribution goes from being bimodal with a negative correlation under panmixia to unimodal with a positive correlation under limited dispersal.

Predation shapes sperm performance surfaces in guppies

Sperm velocity is a key determinant of competitive fertilization success in many species. Selection is therefore expected to favour the evolution of faster sperm when the level of sperm competition is high. However, several aspects can determine the direction and strength of selection acting on this key performance trait, including ecological factors that influence both sperm competition and the strength of selection acting on correlated traits that may constrain evolutionary responses in sperm velocity. Here, we determine how a key ecological variable, the level of predation, shapes sperm swimming speed across 18 Trinidadian populations of guppies ( Poecilia reticulata). We use performance analysis, a statistical tool akin to the familiar methods of multivariate selection analyses, to determine how the level of predation influences sperm velocity (modelled as a performance trait) when accounting for correlated pre- and postcopulatory traits that are also impacted by predation. We show that predation affects the combination of pre- and postcopulatory traits that ultimately predict sperm performance. Overall, we report evidence for disruptive relationships between sperm performance and combinations of ornaments and sperm morphology, but the specific combinations of traits that predict sperm velocity depended on the level of predation. These analyses underscore the complex nonlinear interrelationships among pre- and postcopulatory traits and the importance of considering ecological factors that may ultimately change the way in which multiple traits interact to determine a trait's performance value. As such, our results are likely to be broadly applicable across systems where selection is influenced by ecological conditions.

Bite performance surfaces of three ecologically divergent Iguanidae lizards: relationships with lower jaw bones

Traits that interact to perform an ecologically relevant function are expected to be under multivariate non-linear selection. Using the lower jaw morphology as a biomechanical model, we test the hypothesis that lower jaw bones of lizards are subjected to stabilizing and correlational selection, associated with mechanical advantage and maximum bite force. We used three closely related tropidurine species that differ in size, head shape and microhabitat: Eurolophosaurus nanuzae, Tropidurus hispidus and Tropidurus semitaeniatus. We predicted a common pattern of correlational selection on bones that are part of in-levers or part of the out-lever of the lower jaw. The predicted pattern was found in E. nanuzae and T. hispidus, but this could not be shown to be statistically significant. For T. semitaeniatus, we found significant disruptive selection on a contrast involving the surangular, and also significant directional selection on linear combinations of traits in all species. The results indicate that the non-linear selection on lower jaw bones does not reflect an optimum to enhance mechanical advantage in all species. Divergent functional demands and specific ecological contexts of species seem relevant in shaping patterns of selection on morphology.

Genetic and environmental variation in transcriptional expression of seminal fluid proteins

Seminal fluid proteins (SFPs) are crucial mediators of sexual selection and sexual conflict. Recent studies have chiefly focused on environmentally induced plasticity as one source of variation in SFP expression, particularly in response to differing sperm competition levels. However, understanding the evolution of a trait in heterogenous environments requires estimates of both environmental and genetic sources of variation, as well as their interaction. Therefore, we investigated how environment (specifically mating group size, a good predictor of sperm competition intensity), genotype and genotype-by-environment interactions affect seminal fluid expression. To do so, we reared 12 inbred lines of a simultaneously hermaphroditic flatworm Macrostomum lignano in groups of either two or eight worms and measured the expression levels of 58 putative SFP transcripts. We then examined the source of variation in the expression of each transcript individually and for multivariate axes extracted from a principal component analysis. We found that mating group size did not affect expression levels according to the single transcript analyses, nor did it affect the first principal component (presumably representing overall investment in seminal fluid production). However, mating group size did affect the relative expression of different transcripts captured by the second principal component (presumably reflecting variation in seminal fluid composition). Most transcripts were genetically variable in their expression level and several exhibited genotype-by-environment interactions; relative composition also showed high genetic variation. Collectively, our results reveal the tightly integrated nature of the seminal fluid transcriptome and provide new insights into the quantitative genetic basis of seminal fluid investment and composition.

Phenotypic selection on ponderosa pine seed and seedling traits in the field under three experimentally manipulated drought treatments

Drought-related selection during seedling emergence and early development may play a strong role in adaptation. Yet this process is poorly understood and particularly so in relation to ongoing climate change. To evaluate drought-induced differences in selection during early life stages, a total of 50 maternal families sampled from three climatically disparate ponderosa pine (Pinus ponderosa Doug.) populations were grown from seed in two common garden field experiments at a location that was warmer and drier than seed origins. Three drought treatments were imposed experimentally. Phenotypic selection was assessed by relating plant fitness measured as survival or unconditional expected height at age 3 to seed density (mass per unit volume), date of emergence, and timing of shoot elongation. In the year of emergence from seed, differential mortality was particularly strong and clearly indicated selection. In contrast, selection in subsequent years was far less pronounced. Phenotypes with high seed density, an intermediate but relatively early emergence date, and high 2nd-year early-season shoot elongation exhibited the greatest estimated fitness under drought. The form of selection varied among seed sources in relation to drought treatment. Selection was generally more acute in the cases of greatest difference between drought treatment and climatic patterns of precipitation at the site of seed origin. These results suggest that populations of ponderosa pine are differentially adapted to drought patterns associated with the climate of their origin. To the extent that the phenotypic traits examined are heritable or correlated with heritable traits, our results provide insight into how tree populations may evolve in response to drought.

Dominance reversals and the maintenance of genetic variation for fitness

Antagonistic selection between different fitness components (e.g., survival versus fertility) or different types of individuals in a population (e.g., females versus males) can potentially maintain genetic diversity and thereby account for the high levels of fitness variation observed in natural populations. However, the degree to which antagonistic selection can maintain genetic variation critically depends on the dominance relations between antagonistically selected alleles in diploid individuals. Conditions for stable polymorphism of antagonistically selected alleles are narrow, particularly when selection is weak, unless the alleles exhibit "dominance reversals"-in which each allele is partially or completely dominant in selective contexts in which it is favored and recessive in contexts in which it is harmful. Although theory predicts that dominance reversals should emerge under biologically plausible conditions, evidence for dominance reversals is sparse. In this primer, we review theoretical arguments and data supporting a role for dominance reversals in the maintenance of genetic variation. We then highlight an illuminating new study by Grieshop and Arnqvist, which reports a genome-wide signal of dominance reversals between male and female fitness in seed beetles.

Prevalence of disruptive selection predicts extent of species differentiation in Lake Victoria cichlids

Theory suggests that speciation with gene flow is most likely when both sexual and ecological selection are divergent or disruptive. Divergent sexual and natural selection on the visual system have been demonstrated before in sympatric, morphologically similar sister species of Lake Victoria cichlids, but this does not explain the subtle morphological differences between them. To investigate the significance of natural selection on morphology during speciation, we here ask whether the prevalence of disruptive ecological selection differs between sympatric sister species that are at different stages of speciation. Some of our species pairs do (Pundamilia) and others do not (Neochromis) differ distinctively in sexually selected male nuptial coloration. We find that (i) evidence for disruptive selection, and for evolutionary response to it, is prevalent in traits that are differentiated between sister species; (ii) prevalence of both predicts the extent of genetic differentiation; and (iii) genetic differentiation is weaker in species pairs with conserved male nuptial coloration. Our results speak to the existence of two different mechanisms of speciation with gene flow: speciation mainly by sexual selection tightly followed by ecological character displacement in some cases and speciation mainly by divergent ecological selection in others.

Phenotypic selection on growth rhythm in whitebark pine under climatic conditions warmer than seed origins

Growth rhythm that is well synchronized with seasonal changes in local climatic conditions is understood to enhance fitness; however, rapid ongoing climate change threatens to disrupt this synchrony. To evaluate phenotypic selection on growth rhythm under expected warmer and drier future climate, seedlings from 49 populations of whitebark pine (Pinus albicaulis Engelm.) were grown and measured over more than 10 years in two common garden field experiments on sites that approximate the projected future climate of the seed origins. Selection on growth rhythm was assessed by relating individual plant fitness to timing and rate of shoot elongation. Differential survival clearly evidenced selection on growth rhythm. We detected directional and stabilizing selection that varied in magnitude between experimental sites and among years. The observed phenotypic selection supports the interpretation of clinal variation among populations within tree species as reflecting adaptive variation in response to past natural selection mediated by climate. To the extent that growth rhythm is heritable, results of the present study suggest evolution of whitebark pine toward a more distinct timing of shoot elongation and generally more rapid elongation in the immediate next generation under ongoing climate change in environments similar to the study sites.

Dynamic changes in display architecture and function across environments revealed by a systems approach to animal communication

Animal communication is often structurally complex and dynamic, with signaler and receiver behavior varying in response to multiple environmental factors. To date, studies assessing signal dynamics have mostly focused on the relationships between select signaling traits and receiver responses in a single environment. We use the wolf spider Schizocosa floridana to explore the relationships between courtship display form and function across two social contexts (female presence vs absence) and two light environments (light vs dark). We use traditional analytical methods to determine predictors of copulation success (i.e., signal function) and examine these predictors in a structural context by overlaying them on signal phenotype networks (Wilkins et al. 2015). This allows us to explore system design principles (degeneracy, redundancy, pluripotentiality), providing insight into hypotheses regarding complex signal evolution. We found that both social context and light environment affect courtship structure, although the predictors of mating success remain similar across light environments, suggesting system degeneracy. Contrastingly, the same display traits may serve different functions across social environments, suggesting pluripotentiality. Ultimately, our network approach uncovers a complexity in display structure and function that is missed by functional analyses alone, highlighting the importance of systems-based methodologies for understanding the dynamic nature of complex signals.

Unravelling anisogamy: egg size and ejaculate size mediate selection on morphology in free-swimming sperm

Gamete dimorphism (anisogamy) defines the sexes in most multicellular organisms. Theoretical explanations for its maintenance usually emphasize the size-related selection pressures of sperm competition and zygote survival, assuming that fertilization of all eggs precludes selection for phenotypes that enhance fertility. In external fertilizers, however, fertilization is often incomplete due to sperm limitation, and the risk of polyspermy weakens the advantage of high sperm numbers that is predicted to limit sperm size, allowing alternative selection pressures to target free-swimming sperm. We asked whether egg size and ejaculate size mediate selection on the free-swimming sperm of Galeolaria caespitosa, a marine tubeworm with external fertilization, by comparing relationships between sperm morphology and male fertility across manipulations of egg size and sperm density. Our results suggest that selection pressures exerted by these factors may aid the maintenance of anisogamy in external fertilizers by limiting the adaptive value of larger sperm in the absence of competition. In doing so, our study offers a more complete explanation for the stability of anisogamy across the range of sperm environments typical of this mating system and identifies new potential for the sexes to coevolve via mutual selection pressures exerted by gametes at fertilization.

Artificial Selection to Increase the Phenotypic Variance in gmax Fails

Stabilizing selection is important in evolutionary theories of the maintenance of genetic variance and has been invoked as the key process determining macroevolutionary patterns of trait evolution. However, manipulative evidence for the extent of stabilizing selection, particularly on multivariate traits, is lacking. We used artificial disruptive selection in Drosophila serrata as a tool to determine the relative strength of stabilizing selection experienced by multivariate trait combinations with contrasting levels of genetic and mutational variance. Contrary to expectation, when disruptive selection was applied to the major axis of standing genetic variance, g_max, we observed a significant and repeatable decrease in its phenotypic variance. In contrast, the multivariate trait combination predicted to be under strong stabilizing selection showed a significant and repeatable increase in its phenotypic variance. Correlated responses were observed in all selection treatments, and viability selection operating on extreme phenotypes of traits genetically correlated with those directly selected on limited our ability to increase their phenotypic range. Our manipulation revealed that multivariate trait combinations were subject to stabilizing selection; however, we did not observe a direct relationship between the strength of stabilizing selection and the levels of standing genetic variance in multivariate trait combinations. Contrasting patterns of allele frequencies underlying traits with high versus low levels of standing genetic variance may be implicated in determining the response to artificial selection in multivariate trait combinations.

The sex-specific effects of diet quality versus quantity on morphology in Drosophila melanogaster

Variation in the quality and quantity of nutrition is a major contributor to phenotypic variation in animal populations. Although we know much of how dietary restriction impacts phenotype, and of the molecular-genetic and physiological mechanisms that underlie this response, we know much less of the effects of dietary imbalance. Specifically, although dietary imbalance and restriction both reduce overall body size, it is unclear whether both have the same effect on the size of individual traits. Here, we use the fruit fly Drosophila melanogaster to explore the effect of dietary food versus protein-to-carbohydrate ratio on body proportion and trait size. Our results indicate that body proportion and trait size respond differently to changes in diet quantity (food concentration) versus diet quality (protein-to-carbohydrate ratio), and that these effects are sex specific. While these differences suggest that Drosophila use at least partially distinct developmental mechanisms to respond to diet quality versus quantity, further analysis indicates that the responses can be largely explained by the independent and contrasting effects of protein and carbohydrate concentration on trait size. Our data highlight the importance of considering macronutrient composition when elucidating the effect of nutrition on trait size, at the levels of both morphology and developmental physiology.