Correlational selection and the evolution of genomic architecture

Diversity of Intraspecific Patterns of Brain Region Size Covariation in Fish

Traits often do not evolve in isolation or vary independently of other traits. Instead, they can be affected by covariation, both within and across species. However, the importance of within-species trait covariation and, critically, the degree to which it varies between species has yet to be thoroughly studied. Brain morphology is a trait of great ecological and behavioral importance, with regions that are hypothesized to vary in size based on behavioral and cognitive demands. Sizes of brain regions have also been shown to covary with each other across various taxa. Here, we test the degree to which covariation in brain region sizes within species has been conserved across 10 teleost fish species. These 10 species span five orders, allowing us to examine how phylogenetic proximity influences similarities in intraspecific trait covariation. Our results showed a trend that similar patterns of brain region size covariation occur in more closely related species. Interestingly, there were certain brain region pairs that showed similar levels of covariation across all species regardless of phylogenetic distance, such as the telencephalon and optic tectum, while others, such as the olfactory bulb and the hypothalamus, varied more independently. Ultimately, the patterns of brain region covariation shown here suggest that evolutionary mechanisms or constraints can act on specific brain regions independently, and that these constraints can change over evolutionary time.

Swimming performance, but not metabolism, is related to a boldness-activity syndrome in schoolmaster snapper (Lutjanus apodus)

Commercial overexploitation and climate change can alter the physiology and behavior of marine organisms, although intraspecific phenotypic responses to such changes can vary greatly depending on the environment, species, and severity of the stressor. Under the pace-of-life syndrome (POLS) hypothesis, behavior, physiology, and life-history traits are linked, and thus, affected by selection targeting any aspect of organismal biology. However, these links are understudied in tropical marine fishes, and further work is needed to better understand the impacts of fisheries and climate change on wild stocks. Moreover, tropical regions have a greater reliance on fisheries; thus investigations should focus on species with substantial socioeconomic value to ensure benefits at the local level. This study aimed to address this need by measuring the behavior (boldness and activity), metabolism, and swimming performance (using a critical swim speed [Ucrit] test) of schoolmaster snapper Lutjanus apodus in Eleuthera, the Bahamas. We report a strong positive correlation between boldness and activity, high repeatability of these behavioral metrics, and two groupings that were consistent with "proactive" and "reactive" behavioral types. These behavioral types differed significantly in their swimming performance, with reactive individuals having a 13.1% higher mean Ucrit. In contrast, no significant differences were found in the measured metabolic parameters between behavioral types. This study is the first to investigate the intraspecific links between behavior and physiology in a snapper species, using the novel and ecologically relevant comparison of Ucrit with behavioral syndrome types. These data suggest that additional research is needed to better predict the success of proactive/reactive tropical fish if overexploited and as influenced by climate change.

Mining expressed sequence tag (EST) microsatellite markers to assess the genetic differentiation of five Hynobius species endemic to Taiwan

Taiwan harbors five endemic species of salamanders (Hynobius spp.) that inhabit distinct alpine regions, contributing to population fragmentation across isolated "sky islands". With an evolutionary history spanning multiple glacial-interglacial cycles, these species represent an exceptional paradigm for exploring biogeography and speciation. However, a lack of suitable genetic markers applicable across species has limited research efforts. Thus, developing cross-amplifying markers is imperative. Expressed sequence-tag simple-sequence repeats (EST-SSRs) that amplify across divergent lineages are ideal for species identification in instances where phenotypic differentiation is challenging. Here, we report a suite of cross-amplifying EST-SSRs from the transcriptomes of the five Hynobius species that exhibit an interspecies transferability rate of 67.67%. To identify individual markers exhibiting cross-species polymorphism and to assess interspecies genetic diversity, we assayed 140 individuals from the five species across 84 sampling sites. A set of EST-SSRs with a high interspecies polymorphic information content (PIC = 0.63) effectively classified these individuals into five distinct clusters, as supported by discriminant analysis of principal components (DAPC), STRUCTURE assignment tests, and Neighbor-joining trees. Moreover, pair-wise FST values > 0.15 indicate notable between-cluster genetic divergence. Our set of 20 polymorphic EST-SSRs is suitable for assessing population structure within and among Hynobius species, as well as for long-term monitoring of their genetic composition.

A meta-analysis of whole-body and heart mass effect sizes from a long-term artificial selection experiment for high voluntary exercise

Selection experiments play an increasingly important role in comparative and evolutionary physiology. However, selection experiments can be limited by relatively low statistical power, in part because replicate line is the experimental unit for analyses of direct or correlated responses (rather than number of individuals measured). One way to increase the ability to detect correlated responses is through a meta-analysis of studies for a given trait across multiple generations. To demonstrate this, we applied meta-analytic techniques to two traits (body mass and heart ventricle mass, with body mass as a covariate) from a long-term artificial selection experiment for high voluntary wheel-running behavior. In this experiment, all four replicate High Runner (HR) lines reached apparent selection limits around generations 17-27, running approximately 2.5- to 3-fold more revolutions per day than the four non-selected Control (C) lines. Although both traits would also be expected to change in HR lines (relative heart size expected to increase, expected direction for body mass is less clear), their statistical significance has varied, despite repeated measurements. We compiled information from 33 unique studies and calculated a measure of effect size (Pearson's R). Our results indicate that, despite a lack of statistical significance in most generations, HR mice have evolved larger hearts and smaller bodies relative to controls. Moreover, plateaus in effect sizes for both traits coincide with the generational range during which the selection limit for wheel-running behavior was reached. Finally, since reaching the selection limit, absolute effect sizes for body mass and heart ventricle mass have become smaller (i.e. closer to 0).

The molding of intraspecific trait variation by selection under ecological inheritance

Organisms continuously modify their environment, often impacting the fitness of future conspecifics due to ecological inheritance. When this inheritance is biased toward kin, selection favors modifications that increase the fitness of downstream individuals. How such selection shapes trait variation within populations remains poorly understood. Using mathematical modelling, we investigate the coevolution of multiple traits in a group-structured population when these traits affect the group environment, which is then bequeathed to future generations. We examine when such coevolution favors polymorphism as well as the resulting associations among traits. We find in particular that two traits become associated when one trait affects the environment while the other influences the likelihood that future kin experience this environment. To illustrate this, we model the coevolution of (a) the attack rate on a local renewable resource, which deteriorates environmental conditions, with (b) dispersal between groups, which reduces the likelihood that kin suffers from such deterioration. We show this often leads to the emergence of two highly differentiated morphs: one that readily disperses and depletes local resources, and another that maintains these resources and tends to remain philopatric. More broadly, we suggest that ecological inheritance can contribute to phenotypic diversity and lead to complex polymorphism.

Activity-aggression behavioural syndromes exist in males but not in females of the field cricket Teleogryllus emma

Previous studies on sex differences in behaviour have largely focused on differences in average behaviours between sexes. However, males and females can diverge not only in average behaviours but also in the direction of behavioural correlations at the individual level (i.e. behavioural syndromes). Behavioural syndromes, with their potential to constrain the independent evolution of behaviours, may play a role in shaping sex-specific responses to selection and contributing to the development of sex differences in behaviour. Despite the pivotal role of behavioural syndromes in the evolution of sexual dimorphism in behaviour, robust empirical evidence of sex differences in behavioural syndromes based on repeated measurements of behaviours is scarce. In this study, we conducted repeated measurements of activity and aggression in male and female field crickets Teleogryllus emma, providing evidence of sex differences in the existence of behavioural syndromes. Males exhibited a significantly positive behavioural syndrome between activity and aggression, whereas females, in contrast, did not show any aggressive behaviour, resulting in the absence of such a syndrome. The sex differences in the existence of the activity-aggression behavioural syndromes in this species could be attributed to differences in selection. Selection favouring more active and aggressive males may have shaped a positive activity-aggression behavioural syndrome in males, whereas the absence of selection favouring female aggression may have resulted in the absence of aggression and the related behavioural syndrome in females. However, given the plasticity of behaviour with changes in age or the environment, further research is needed to explore how sex differences in the existence of activity-aggression behavioural syndromes change across contexts. Furthermore, understanding the genetic underpinning of sex differences in a behavioural syndrome would be pivotal to assess the role of behavioural syndromes in the evolution of sexual dimorphism in behaviours.

Detecting natural selection in trait-trait coevolution

No phenotypic trait evolves independently of all other traits, but the cause of trait-trait coevolution is poorly understood. While the coevolution could arise simply from pleiotropic mutations that simultaneously affect the traits concerned, it could also result from multivariate natural selection favoring certain trait relationships. To gain a general mechanistic understanding of trait-trait coevolution, we examine the evolution of 220 cell morphology traits across 16 natural strains of the yeast Saccharomyces cerevisiae and the evolution of 24 wing morphology traits across 110 fly species of the family Drosophilidae, along with the variations of these traits among gene deletion or mutation accumulation lines (a.k.a. mutants). For numerous trait pairs, the phenotypic correlation among evolutionary lineages differs significantly from that among mutants. Specifically, we find hundreds of cases where the evolutionary correlation between traits is strengthened or reversed relative to the mutational correlation, which, according to our population genetic simulation, is likely caused by multivariate selection. Furthermore, we detect selection for enhanced modularity of the yeast traits analyzed. Together, these results demonstrate that trait-trait coevolution is shaped by natural selection and suggest that the pleiotropic structure of mutation is not optimal. Because the morphological traits analyzed here are chosen largely because of their measurability and thereby are not expected to be biased with regard to natural selection, our conclusion is likely general.

ebony underpins Batesian mimicry in melanic stoneflies

The evolution of Batesian mimicry - whereby harmless species avoid predation through their resemblance to harmful species - has long intrigued biologists. In rare cases, Batesian mimicry is linked to intraspecific colour variation, in which only some individuals within a population resemble a noxious 'model'. Here, we assess intraspecific colour variation within a widespread New Zealand stonefly, wherein highly melanized individuals of Zelandoperla closely resemble a chemically defended aposematic stonefly, Austroperla cyrene. We assess convergence in the colour pattern of these two species, compare their relative palatability to predators, and use genome-wide association mapping to assess the genetic basis of this resemblance. Our analysis reveals that melanized Zelandoperla overlap significantly with Austroperla in colour space but are significantly more palatable to predators, implying that they are indeed Batesian mimics. Analysis of 194,773 genome-wide SNPs reveals an outlier locus (ebony) strongly differentiating melanic versus non-melanic Zelandoperla. Genotyping of 338 specimens from a single Zelandoperla population indicates that ebony explains nearly 70% of the observed variance in melanism. As ebony has a well-documented role in insect melanin biosynthesis, our findings indicate this locus has a conserved function across deeply divergent hexapod lineages. Distributional records suggest a link between the occurrence of melanic Zelandoperla and the forested ecosystems where the model Austroperla is abundant, suggesting the potential for adaptive shifts in this system underpinned by environmental change.

Are animal personality, body condition, physiology and structural size integrated? A comparison of species, populations and sexes, and the value of study replication

A major question in behavioural ecology is why behaviour, physiology and morphology are often integrated into syndromes. In great tits, Parus major, for example, explorative males are larger (vs. smaller) and leaner (vs. heavier) compared to less explorative individuals. Unfortunately, considerable debate exists on whether patterns found in specific studies are replicable. This debate calls for study replication among species, populations and sexes. We measured behavioural (exploration), physiological (breathing rate) and morphological traits (body mass, tarsus length, wing length, bill length) in two species (great vs. blue tits Cyanistes caeruleus), two populations (Forstenrieder Park vs. Starnberg) and two sexes (males vs. females). We then tested whether the same pattern of integration characterized all unique combinations of these three biological categories (hereafter called datasets). We used a multi-year repeated measures set-up to estimate among-individual trait correlation matrices for each dataset. We then used structural equation modelling to test for size-dependent behaviour and physiology, size-corrected (i.e. size-independent) behaviour-physiology correlations and size-corrected body mass-dependent behaviour and physiology. Finally, we used meta-analyses to test which structural paths were generally (vs. conditionally) supported (vs. unsupported). We found general and consistent support for size-dependent physiology and size-corrected body mass-dependent physiology across datasets: faster breathers were smaller but heavier for their size. Unexpectedly, condition-dependent behaviour was not supported: explorative birds were neither leaner, nor was this relationship heterogeneous across datasets. All other hypothesized patterns were dataset-specific: the covariance between size and behaviour, and between behaviour and physiology differed in sign between datasets, and both were, on average, not supported. This heterogeneity was not explained by any of our moderators: species, population or sex. The specific pattern of size- and condition-dependent physiology reported for a unique combination of species, population, and sex, thus generally predicted those in others. Patterns of size- or condition-dependent behaviour (i.e. 'personality'), or behaviour-physiology syndromes reported in specific datasets, by contrast, did not. These findings call for studies revealing the ecological background of this variation and highlight the value of study replication to help understand whether patterns of phenotypic integration reported in one study can be generalized.

Frank Beach Award Winner: The centrality of the hypothalamic-pituitary-adrenal axis in dealing with environmental change across temporal scales

Understanding if and how individuals and populations cope with environmental change is an enduring question in evolutionary ecology that has renewed importance given the pace of change in the Anthropocene. Two evolutionary strategies of coping with environmental change may be particularly important in rapidly changing environments: adaptive phenotypic plasticity and/or bet hedging. Adaptive plasticity could enable individuals to match their phenotypes to the expected environment if there is an accurate cue predicting the selective environment. Diversifying bet hedging involves the production of seemingly random phenotypes in an unpredictable environment, some of which may be adaptive. Here, I review the central role of the hypothalamic-pituitary-adrenal (HPA) axis and glucocorticoids (GCs) in enabling vertebrates to cope with environmental change through adaptive plasticity and bet hedging. I first describe how the HPA axis mediates three types of adaptive plasticity to cope with environmental change (evasion, tolerance, recovery) over short timescales (e.g., 1-3 generations) before discussing how the implications of GCs on phenotype integration may depend upon the timescale under consideration. GCs can promote adaptive phenotypic integration, but their effects on phenotypic co-variation could also limit the dimensions of phenotypic space explored by animals over longer timescales. Finally, I discuss how organismal responses to environmental stressors can act as a bet hedging mechanism and therefore enhance evolvability by increasing genetic or phenotypic variability or reducing patterns of genetic and phenotypic co-variance. Together, this emphasizes the crucial role of the HPA axis in understanding fundamental questions in evolutionary ecology.

Burning questions for a warming and changing world: 15 unknowns in plant abiotic stress

We present unresolved questions in plant abiotic stress biology as posed by 15 research groups with expertise spanning eco-physiology to cell and molecular biology. Common themes of these questions include the need to better understand how plants detect water availability, temperature, salinity, and rising carbon dioxide (CO2) levels; how environmental signals interface with endogenous signaling and development (e.g. circadian clock and flowering time); and how this integrated signaling controls downstream responses (e.g. stomatal regulation, proline metabolism, and growth versus defense balance). The plasma membrane comes up frequently as a site of key signaling and transport events (e.g. mechanosensing and lipid-derived signaling, aquaporins). Adaptation to water extremes and rising CO2 affects hydraulic architecture and transpiration, as well as root and shoot growth and morphology, in ways not fully understood. Environmental adaptation involves tradeoffs that limit ecological distribution and crop resilience in the face of changing and increasingly unpredictable environments. Exploration of plant diversity within and among species can help us know which of these tradeoffs represent fundamental limits and which ones can be circumvented by bringing new trait combinations together. Better defining what constitutes beneficial stress resistance in different contexts and making connections between genes and phenotypes, and between laboratory and field observations, are overarching challenges.

Supergenes on steroids

At the birth of supergenes, the genomic landscape is dramatically re-organized leading to pronounced differences in phenotypes and increased intrasexual diversity. Two of the best-studied supergenes in vertebrates are arguably the inversion polymorphisms on chromosomes 2 and 11 in the white-throated sparrow (Zonotrichia albicollis) and the ruff (Calidris pugnax), respectively. In both species, regions of suppressed recombination determine plumage coloration and social behavioural phenotypes. Despite the apparent lack of gene overlap between these two supergenes, in both cases the alternative phenotypes seem to be driven largely by alterations in steroid hormone pathways. Here, we explore the interplay between genomic architecture and steroid-related genes. Due to the highly pleiotropic effects of steroid-related genes and their universal involvement in social behaviour and transcriptomic regulation, processes favouring their linkage are likely to have substantial effects on the evolution of behavioural phenotypes, individual fitness, and life-history strategies. We propose that inversion-related differentiation and regulatory changes in steroid-related genes lie at the core of phenotypic differentiation in both of these interesting species. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Phenotypic and genomic signatures of interspecies cooperation and conflict in naturally occurring isolates of a model plant symbiont

Given the need to predict the outcomes of (co)evolution in host-associated microbiomes, whether microbial and host fitnesses tend to trade-off, generating conflict, remains a pressing question. Examining the relationships between host and microbe fitness proxies at both the phenotypic and genomic levels can illuminate the mechanisms underlying interspecies cooperation and conflict. We examined naturally occurring genetic variation in 191 strains of the model microbial symbiont Sinorhizobium meliloti, paired with each of two host Medicago truncatula genotypes in single- or multi-strain experiments to determine how multiple proxies of microbial and host fitness were related to one another and test key predictions about mutualism evolution at the genomic scale, while also addressing the challenge of measuring microbial fitness. We found little evidence for interspecies fitness conflict; loci tended to have concordant effects on both microbe and host fitnesses, even in environments with multiple co-occurring strains. Our results emphasize the importance of quantifying microbial relative fitness for understanding microbiome evolution and thus harnessing microbiomes to improve host fitness. Additionally, we find that mutualistic coevolution between hosts and microbes acts to maintain, rather than erode, genetic diversity, potentially explaining why variation in mutualism traits persists in nature.

Computer vision for assessing species color pattern variation from web-based community science images

Openly available community science digital vouchers provide a wealth of data to study phenotypic change across space and time. However, extracting phenotypic data from these resources requires significant human effort. Here, we demonstrate a workflow and computer vision model for automatically categorizing species color pattern from community science images. Our work is focused on documenting the striped/unstriped color polymorphism in the Eastern Red-backed Salamander (Plethodon cinereus). We used an ensemble convolutional neural network model to analyze this polymorphism in 20,318 iNaturalist images. Our model was highly accurate (∼98%) despite image heterogeneity. We used the resulting annotations to document extensive niche overlap between morphs, but wider niche breadth for striped morphs at the range-wide scale. Our work showcases key design principles for using machine learning with heterogeneous community science image data to address questions at an unprecedented scale.

•

We built a deep learning model to group color morphs from community science images

•

Our model achieved 98% accuracy for classifying striped and unstriped salamanders

•

We used our model to classify >20,000 images and built morph-specific niche models

•

We then determined if Red-backed salamanders niche partition at a range-wide scale

Genomic changes underlying repeated niche shifts in an adaptive radiation

In adaptive radiations, single lineages rapidly diversify by adapting to many new niches. Little is known yet about the genomic mechanisms involved, that is, the source of genetic variation or genomic architecture facilitating or constraining adaptive radiation. Here, we investigate genomic changes associated with repeated invasion of many different freshwater niches by threespine stickleback in the Haida Gwaii archipelago, Canada, by resequencing single genomes from one marine and 28 freshwater populations. We find 89 likely targets of parallel selection in the genome that are enriched for old standing genetic variation. In contrast to theoretical expectations, their genomic architecture is highly dispersed with little clustering. Candidate genes and genotype-environment correlations match the three major environmental axes predation regime, light environment, and ecosystem size. In a niche space with these three dimensions, we find that the more divergent a new niche from the ancestral marine habitat, the more loci show signatures of parallel selection. Our findings suggest that the genomic architecture of parallel adaptation in adaptive radiation depends on the steepness of ecological gradients and the dimensionality of the niche space.

Multispecies colour polymorphisms associated with contrasting microhabitats in two Mediterranean wrasse radiations

Intraspecific colour polymorphisms (CPs) present unique opportunities to study fundamental evolutionary questions, such as the link between ecology and phenotype, mechanisms maintaining genetic diversity and their putative role in speciation. Wrasses are highly diverse in ecology and morphology and harbour a variety of colour‐polymorphic species. In the Mediterranean Sea, wrasses of the tribe Labrini evolved two species radiations each harbouring several species with a brown and a green morph. The colour morphs occur in complete sympatry in mosaic habitats with rocky outcrops and Neptune grass patches. Morph‐specific differences had not been characterized yet and the evolutionary forces maintaining them remained unknown. With genome‐wide data for almost all Labrini species, we show that species with CPs are distributed across the phylogeny, but show evidence of hybridization. This suggests that the colour morphs are either ancient and have been lost repeatedly, that they have evolved repeatedly or have been shared via hybridization. Focusing on two polymorphic species, we find that each colour morph is more common in the microhabitat providing the best colour match and that the morphs exhibit additional behavioural and morphological differences further improving crypsis in their respective microhabitats. We find little evidence for genetic differentiation between the morphs in either species. Therefore, we propose that these colour morphs represent a multi‐niche polymorphism as an adaptation to the highly heterogeneous habitat. Our study highlights how colour polymorphism (CP) can be advantageous in mosaic habitats and that Mediterranean wrasses are an ideal system to study trans‐species polymorphisms, i.e. polymorphisms maintained across several species, in adaptive radiations.

Intraspecific diversity alters the relationship between climate change and parasitism in a polymorphic ectotherm

Climate-modulated parasitism is driven by a range of factors, yet the spatial and temporal variability of this relationship has received scant attention in wild vertebrate hosts. Moreover, most prior studies overlooked the intraspecific differences across host morphotypes, which impedes a full understanding of the climate-parasitism relationship. In the common lizard (Zootoca vivipara), females exhibit three colour morphs: yellow (Y-females), orange (O-females) and mixed (mixture of yellow and orange, M-females). Zootoca vivipara is also infested with an ectoparasite (Ophionyssus mites). We therefore used this model system to examine the intraspecific response of hosts to parasitism under climate change. We found infestation probability to differ across colour morphs at both spatial (10 sites) and temporal (20 years) scales: M-females had lower parasite infestations than Y- and O-females at lower temperatures, but became more susceptible to parasites as temperature increased. The advantage of M-females at low temperatures was counterbalanced by their higher mortality rates thereafter, which suggests a morph-dependent trade-off between resistance to parasites and host survival. Furthermore, significant interactions between colour morphs and temperature indicate that the relationship between parasite infestations and climate warming was contingent on host morphotypes. Parasite infestations increased with temperature for most morphs, but displayed morph-specific rates. Finally, infested M-females had higher reductions in survival rates than infested Y- or O-females, which implies a potential loss of intraspecific diversity within populations as parasitism and temperatures rise. Overall, we found parasitism increases with warming temperatures, but this relationship is modulated by host morphotypes and an interaction with temperature. We suggest that epidemiological models incorporate intraspecific diversity within species for better understanding the dynamics of wildlife diseases under climate warming.

Viability, behavior, and color expression in the offspring of matings between common wall lizard Podarcis muralis color morphs

Color polymorphisms are widely studied to identify the mechanisms responsible for the origin and maintenance of phenotypic variability in nature. Two of the mechanisms of balancing selection currently thought to explain the long-term persistence of polymorphisms are the evolution of alternative phenotypic optima through correlational selection on suites of traits including color and heterosis. Both of these mechanisms can generate differences in offspring viability and fitness arising from different morph combinations. Here, we examined the effect of parental morph combination on fertilization success, embryonic viability, newborn quality, antipredator, and foraging behavior, as well as inter-annual survival by conducting controlled matings in a polymorphic lacertid Podarcis muralis, where color morphs are frequently assumed to reflect alternative phenotypic optima (e.g., alternative reproductive strategies). Juveniles were kept in outdoor tubs for a year in order to study inter-annual growth, survival, and morph inheritance. In agreement with a previous genome-wide association analysis, morph frequencies in the year-old juveniles matched the frequencies expected if orange and yellow expressions depended on recessive homozygosity at 2 separate loci. Our findings also agree with previous literature reporting higher reproductive output of heavy females and the higher overall viability of heavy newborn lizards, but we found no evidence for the existence of alternative breeding investment strategies in female morphs, or morph-combination effects on offspring viability and behavior. We conclude that inter-morph breeding remains entirely viable and genetic incompatibilities are of little significance for the maintenance of discrete color morphs in P. muralis from the Pyrenees.

Color Polymorphism is a Driver of Diversification in the Lizard Family Lacertidae

Color polymorphism-two or more heritable color phenotypes maintained within a single breeding population-is an extreme type of intraspecific diversity widespread across the tree of life. Color polymorphism is hypothesized to be an engine for speciation, where morph loss or divergence between distinct color morphs within a species results in the rapid evolution of new lineages, and thus, color polymorphic lineages are expected to display elevated diversification rates. Multiple species in the lizard family Lacertidae are color polymorphic, making them an ideal group to investigate the evolutionary history of this trait and its influence on macroevolution. Here, we produce a comprehensive species-level phylogeny of the lizard family Lacertidae to reconstruct the evolutionary history of color polymorphism and test if color polymorphism has been a driver of diversification. Accounting for phylogenetic uncertainty with multiple phylogenies and simulation studies, we estimate an ancient origin of color polymorphism (111 Ma) within the Lacertini tribe (subfamily Lacertinae). Color polymorphism most likely evolved few times in the Lacertidae and has been lost at a much faster rate than gained. Evolutionary transitions to color polymorphism are associated with shifts in increased net diversification rate in this family of lizards. Taken together, our empirical results support long-standing theoretical expectations that color polymorphism is a driver of diversification.[Color polymorphism; Lacertidae; state-dependent speciation extinction models; trait-dependent diversification.].

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.

Rapid phenotypic change in a polymorphic salamander over 43 years

Color polymorphic animals offer a unique system for studying intraspecific phenotypic responses to climate change. Discrete color morphs are easy to identify, and correlated trait responses of morphs can indicate how climate warming may facilitate long-term maintenance of polymorphisms. We use a historical dataset spanning 43 years to examine temporal shifts in color morph frequency and body size in response to climate in the Eastern Red-backed Salamander, Plethodon cinereus, which contains a widespread striped/unstriped color polymorphism. We created a pipeline to extract high-throughput trait data from fluid-preserved museum specimens where we batch-photographed salamanders, de-aggregated individual specimens from photographs, and solicited help of community scientists to score color morphs. We used a linear modeling framework that includes information about spatial population structure to demonstrate that color morph frequency and body size vary in response to climate, elevation, and over time, with an overall trend of higher frequency and decreased body size of the striped morph, but increased size of the unstriped morph. These surprising results suggest that morphs may be responding to multiple climate and geographic drivers through co-adapted morphological changes. This work highlights new practices of extracting trait data from museum specimens to demonstrate species phenotypes response to climate change.

Multivariate analysis of morphology, behaviour, growth and developmental timing in hybrids brings new insights into the divergence of sympatric Arctic charr morphs

BACKGROUND

Studying the development of fitness related traits in hybrids from populations diverging in sympatry is a fundamental approach to understand the processes of speciation. However, such traits are often affected by covariance structures that complicate the comprehension of these processes, especially because the interactive relationships between traits of different nature (e.g. morphology, behaviour, life-history) remain largely unknown in this context. In a common garden setup, we conducted an extensive examination of a large suit of traits putatively involved in the divergence of two morphs of Arctic charr (Salvelinus alpinus), and investigated the consequences of potential patterns of trait covariance on the phenotype of their hybrids. These traits were measured along ontogeny and involved growth, yolk sac resorption, developmental timing (hatching and the onset of exogeneous feeding), head morphology and feeding behaviour.

RESULTS

Growth trajectories provided the strongest signal of phenotypic divergence between the two charr. Strikingly, the first-generation hybrids did not show intermediate nor delayed growth but were similar to the smallest morph, suggesting parental biases in the inheritance of growth patterns. However, we did not observe extensive multivariate trait differences between the two morphs and their hybrids. Growth was linked to head morphology (suggesting that morphological variations in early juveniles relate to simple allometric effects) but this was the only strong signal of covariance observed between all the measured traits. Furthermore, we did not report evidence for differences in overall phenotypic variance between morphs, nor for enhanced phenotypic variability in their hybrids.

CONCLUSION

Our study shed light on the multivariate aspect of development in a context of adaptive divergence. The lack of evidence for the integration of most traits into a single covariance structure suggested that phenotypic constraints may not always favour nor impede divergence toward ecological niches differing in numerous physical and ecological variables, as observed in the respective habitats of the two charr. Likewise, the role of hybridization as a disruptive agent of trait covariance may not necessarily be significant in the evolution of populations undergoing resource polymorphism.

Epistatic interactions between pterin and carotenoid genes modulate intra-morph color variation in a lizard

Color polymorphisms have become a major topic in evolutionary biology and substantial efforts have been devoted to the understanding of the mechanisms responsible for originating such colorful systems. Within-morph continuous variation, on the other hand, has been neglected in most of the studies. Here, we combine spectrophotometric/visual modeling and genetic data to study the mechanisms promoting continuous variation within categorical color morphs of Podarcis muralis. Our results suggest that intra-morph variability in the pterin-based orange morph is greater compared to white and yellow morphs. We also show that continuous variation within the orange morph is partially discriminable by conspecifics. Genotyping results indicate that allelic variants at the BCO2 locus (responsible for deposition of yellow carotenoids) contribute to generate continuous variation in orange individuals. However, other intrinsic and/or extrinsic mechanisms, such as body size, might be involved, opening a new avenue for future research on the drivers of continuous variation within-morphs.

Social selection is density dependent but makes little contribution to total selection in New Zealand giraffe weevils

Social selection occurs when traits of interaction partners influence an individual's fitness and can alter total selection strength. However, we have little idea of what factors influence social selection's strength. Further, social selection only contributes to overall selection when there is phenotypic assortment, but simultaneous estimates of social selection and phenotypic assortment are rare. Here, we estimated social selection on body size in a wild population of New Zealand giraffe weevils (Lasiorhynchus barbicornis). We measured phenotypic assortment by body size and tested whether social selection varied with sex ratio, density and interacted with the body size of the focal individual. Social selection was limited and unaffected by sex ratio or the size of the focal individual. However, at high densities social selection was negative for both sexes, consistent with size-based competitive interactions for access to mates. Phenotypic assortment was always close to zero, indicating negative social selection at high densities will not impede the evolution of larger body sizes. Despite its predicted importance, social selection may only influence evolutionary change in specific contexts, leaving direct selection to drive evolutionary change.

Do colour morphs of wall lizards express different personalities?

Colour morphs sometimes have different behavioural strategies which may be maintained by frequency or density dependence mechanisms. We investigated temporal changes in behavioural reaction to a novel environment among colour morphs (yellow, orange, white) of the European wall lizard (Podarcis muralis). Adult males were given two 15 min experimental trials, and their locomotion was highly consistent between the two trials. Boldness, freezing and escape behaviour were less repeatable. Colour morphs differed in their locomotion and freezing behaviour. Boldness was similar among the morphs, whereas escape behaviour was lowest in yellow morph. Consequently, yellow morph males tended to explore novel environments quickly and thus were more likely to move to potentially safe areas. Orange and white males showed more fear when exposed to a novel environment. Whether such alternative behavioural strategies can contribute to the maintenance of variable fitness optima among the morphs and ultimately to the maintenance of polymorphism remains open to further investigation.

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.

Computer Vision, Machine Learning, and the Promise of Phenomics in Ecology and Evolutionary Biology

For centuries, ecologists and evolutionary biologists have used images such as drawings, paintings and photographs to record and quantify the shapes and patterns of life. With the advent of digital imaging, biologists continue to collect image data at an ever-increasing rate. This immense body of data provides insight into a wide range of biological phenomena, including phenotypic diversity, population dynamics, mechanisms of divergence and adaptation, and evolutionary change. However, the rate of image acquisition frequently outpaces our capacity to manually extract meaningful information from images. Moreover, manual image analysis is low-throughput, difficult to reproduce, and typically measures only a few traits at a time. This has proven to be an impediment to the growing field of phenomics – the study of many phenotypic dimensions together. Computer vision (CV), the automated extraction and processing of information from digital images, provides the opportunity to alleviate this longstanding analytical bottleneck. In this review, we illustrate the capabilities of CV as an efficient and comprehensive method to collect phenomic data in ecological and evolutionary research. First, we briefly review phenomics, arguing that ecologists and evolutionary biologists can effectively capture phenomic-level data by taking pictures and analyzing them using CV. Next we describe the primary types of image-based data, review CV approaches for extracting them (including techniques that entail machine learning and others that do not), and identify the most common hurdles and pitfalls. Finally, we highlight recent successful implementations and promising future applications of CV in the study of phenotypes. In anticipation that CV will become a basic component of the biologist’s toolkit, our review is intended as an entry point for ecologists and evolutionary biologists that are interested in extracting phenotypic information from digital images.

Human protection drives the emergence of a new coping style in animals

Wild animals face novel environmental threats from human activities that may occur along a gradient of interactions with humans. Recent work has shown that merely living close to humans has major implications for a variety of antipredator traits and physiological responses. Here, we hypothesize that when human presence protects prey from their genuine predators (as sometimes seen in urban areas and at some tourist sites), this predator shield, followed by a process of habituation to humans, decouples commonly associated traits related to coping styles, which results in a new range of phenotypes. Such individuals are characterized by low aggressiveness and physiological stress responses, but have enhanced behavioral plasticity, boldness, and cognitive abilities. We refer to these individuals as "preactive," because their physiological and behavioral coping style falls outside the classical proactive/reactive coping styles. While there is some support for this new coping style, formal multivariate studies are required to investigate behavioral and physiological responses to anthropogenic activities.

Genetic mechanisms and correlational selection structure trait variation in a coral snake mimic

Covariation among traits shapes both phenotypic evolution and ecological interactions across space and time. However, rampant geographical variation in the strength and direction of such correlations can be particularly difficult to explain through generalized mechanisms. By integrating population genomics, surveys of natural history collections and spatially explicit analyses, we tested multiple drivers of trait correlations in a coral snake mimic that exhibits remarkable polymorphism in mimetic and non-mimetic colour traits. We found that although such traits co-occur extensively across space, correlations were best explained by a mixture of genetic architecture and correlational selection, rather than by any single mechanism. Our findings suggest that spatially complex trait distributions may be driven more by the simple interaction between multiple processes than by complex variation in one mechanism alone. These interactions are particularly important in mimicry systems, which frequently generate striking geographical variation and genetic correlations among colour pattern traits.

Most published selection gradients are underestimated: Why this is and how to fix it

Ecologists and evolutionary biologists routinely estimate selection gradients. Most researchers seek to quantify selection on individual phenotypes, regardless of whether fixed or repeatedly expressed traits are studied. Selection gradients estimated to address such questions are attenuated unless analyses account for measurement error and biological sources of within-individual variation. Estimates of standardized selection gradients published in Evolution between 2010 and 2019 were primarily based on traits measured once (59% of 325 estimates). We show that those are attenuated: bias increases with decreasing repeatability but differently for linear versus nonlinear gradients. Others derived individual-mean trait values prior to analyses (41%), typically using few repeats per individual, which does not remove bias. We evaluated three solutions, all requiring repeated measures: (i) correcting gradients derived from classic models using estimates of trait correlations and repeatabilities, (ii) multivariate mixed-effects models, previously used for estimating linear gradients (seven estimates, 2%), which we expand to nonlinear analyses, and (iii) errors-in-variables models that account for within-individual variance, and are rarely used in selection studies. All approaches produced accurate estimates regardless of repeatability and type of gradient, however, errors-in-variables models produced more precise estimates and may thus be preferable.

Maternal reproductive output and F1 hybrid fitness may influence contact zone dynamics

The outcome of secondary contact between divergent lineages or species may be influenced by both the reproductive traits of parental species and the fitness of offspring; however, their relative contributions have rarely been evaluated, particularly in longer-lived vertebrate species. We performed pure and reciprocal laboratory crosses between Ctenophorus decresii (tawny dragon) and C. modestus (swift dragon) to examine how parental reproductive traits and ecologically relevant offspring fitness traits may explain contact zone dynamics in the wild. The two species meet in a contact zone of post-F1 hybrids with asymmetric backcrossing and predominantly C. modestus mtDNA haplotypes. We found no evidence for reduced parental fecundity or offspring fitness for F1 hybrid crosses. However, maternal reproductive strategy differed between species, irrespective of the species of their mate. Ctenophorus modestus females had higher fecundity and produced more and larger clutches with lower embryonic mortality. Parental species also influenced sex ratios and offspring traits, with C. modestus ♀ × C. decresii ♂ hybrids exhibiting higher trait values for more fitness measures (growth rate, sprint speed, bite force) than offspring from all other pairings. Together, these patterns are consistent with the prevalence of C. modestus mtDNA in the contact zone, and asymmetric backcrossing likely reflects fitness effects that manifest in the F2 generation. Our results highlight how parental species can influence multiple offspring traits in different ways, which together may combine to influence offspring fitness and shape contact zone dynamics.

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism's potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species' evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module "leaf size gradient" was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under "global regulation" in vegetative bud establishment and development. Lamina basal shape and intercalary veins' number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits' ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

The Biology of Polymorphic Melanic Side-Spotting Patterns in Poeciliid Fishes

Melanin-based color patterns are an emerging model for studying molecular and evolutionary mechanisms driving phenotypic correlations. Extensive literature exists on color patterns and their correlated traits in the family Poeciliidae, indicating that these fishes are tractable models. We review the biology of polymorphic melanic side-spotting patterns characterized by macromelanophores forming irregular spotted patterns across fishes’ flanks. These patterns are present in the generaGambusia, Limia, Phalloceros, Poecilia, andXiphophorus. Their presence is controlled by dominant genes on autosomes or sex chromosomes. Variation in expression is under polygenic control; however, these genes’ identities are still largely unknown. In someGambusia holbrookiandPoecilia latipinna, expression is dependent on low temperature exposure, but underlying molecular mechanisms are unknown. Spotted fish develop melanoma in rare cases and are a well-developed model for melanoma research. Little is known about other physiological correlates except that spottedG. holbrookimales exhibit higher basal cortisol levels than unspotted males and that metabolic rate does not differ between morphs in someXiphophorusspecies. Behavioral differences between morphs are widespread, but specific to population, species, and social context. SpottedG. holbrookimales appear to be more social and more dominant. Juvenile spottedG. holbrookihave lower behavioral flexibility, and spottedX. variatusexhibit greater stress resistance. Findings conflict on whether morphs differ in sexual behavior and in sexual selection by females. Melanic side-spotting patterns are uncommon (<30%) in populations, although extreme high-frequency populations exist. This low frequency is surprising for dominant genes, indicating that a variety of selective pressures influence both these patterns and their correlated traits. Little is known about reproductive life history traits. SpottedG. holbrookiare larger and have higher survival when uncommon, but underlying mechanisms remain unknown. Spotted morphs appear to have a strong selective advantage during predation. Predators prefer to attack and consume unspotted morphs; however, this preference disappears when spottedG. holbrookimales are common, indicating negative frequency-dependent selection. Spotted morphs are preferred socially under turbid conditions, but other environmental factors that shape phenotypic correlations and morph fitness have not been studied. Finally, we present questions for future studies on melanic side-spotting patterns.

Sexual selection on performance traits in an Australian lizard with alternative reproductive tactics

Sexual selection shapes the adaptive landscape in complex ways that lead to trait integration. Much of our understanding of selection comes from studies of morphological traits. However, few studies attempt to quantify the form and direction of selection on performance even though it is predicted to be a more direct target of selection in nature. We measured sexual selection on performance traits (bite force, sprint speed and endurance) in an Australian lizard, the Eastern water skink (Eulamprus quoyii). We first staged 123 contests between size-matched males to investigate whether performance traits were important in determining contest outcome. In a second experiment, we established six breeding populations in large replicate semi-natural enclosures to estimate whether performance traits predicted reproductive success. Our results show that none of the performance measures were important in predicting contest outcome and were not generally strong predictors of reproductive success. However, our analyses suggest a complex fitness landscape driven by males adopting different alternative reproductive tactics (ARTs). We provide a rare test of the role performance plays in sexual selection and highlight the need to test common assumptions regarding the link between maximal performance and fitness. Our results suggest that performance traits may not necessarily be direct targets of sexual selection, but rather indirect targets through their integration with morphological and/or behavioural traits, highlighting a need for more explicit tests of the predicted links between performance and fitness.

Ecology shapes epistasis in a genotype-phenotype-fitness map for stick insect colour

Genetic interactions such as epistasis are widespread in nature and can shape evolutionary dynamics. Epistasis occurs due to nonlinearity in biological systems, which can arise via cellular processes that convert genotype to phenotype and via selective processes that connect phenotype to fitness. Few studies in nature have connected genotype to phenotype to fitness for multiple potentially interacting genetic variants. Thus, the causes of epistasis in the wild remain poorly understood. Here, we show that epistasis for fitness is an emergent and predictable property of nonlinear selective processes. We do so by measuring the genetic basis of cryptic colouration and survival in a field experiment with stick insects. We find that colouration shows a largely additive genetic basis but with some effects of epistasis that enhance differentiation between colour morphs. In terms of fitness, different combinations of loci affecting colouration confer high survival in one host-plant treatment. Specifically, nonlinear correlational selection for specific combinations of colour traits in this treatment drives the emergence of pairwise and higher-order epistasis for fitness at loci underlying colour. In turn, this results in a rugged fitness landscape for genotypes. In contrast, fitness epistasis was dampened in another treatment, where selection was weaker. Patterns of epistasis that are shaped by ecologically based selection could be common and central to understanding fitness landscapes, the dynamics of evolution and potentially other complex systems.

Quantifying selection on standard metabolic rate and body mass in Drosophila melanogaster

Standard metabolic rate (SMR), defined as the minimal energy expenditure required for self-maintenance, is a key physiological trait. Few studies have estimated its relationship with fitness, most notably in insects. This is presumably due to the difficulty of measuring SMR in a large number of very small individuals. Using high-throughput flow-through respirometry and a Drosophila melanogaster laboratory population adapted to a life cycle that facilitates fitness measures, we quantified SMR, body mass, and fitness in 515 female and 522 male adults. We used a novel multivariate approach to estimate linear and nonlinear selection differentials and gradients from the variance-covariance matrix of fitness, SMR, and body mass, allowing traits specific covariates to be accommodated within a single model. In males, linear selection differentials for mass and SMR were positive and individually significant. Selection gradients were also positive but, despite substantial sample sizes, were nonsignificant due to increased uncertainty given strong SMR-mass collinearity. In females, only nonlinear selection was detected and it appeared to act primarily on body size, although the individual gradients were again nonsignificant. Selection did not differ significantly between sexes although differences in the fitness surfaces suggest sex-specific selection as an important topic for further study.

Inside the supergene of the bird with four sexes

The white-throated sparrow (Zonotrichia albicollis) offers unique opportunities to understand the adaptive value of supergenes, particularly their role in alternative phenotypes. In this species, alternative plumage morphs segregate with a nonrecombining segment of chromosome 2, which has been called a 'supergene'. The species mates disassortatively with respect to the supergene; that is, each breeding pair consists of one individual with it and one without it. This species has therefore been called the "bird with four sexes". The supergene segregates with a behavioral phenotype; birds with it are more aggressive and less parental than birds without it. Here, we review our efforts to identify the genes inside the supergene that are responsible for the behavioral polymorphism. The gene ESR1, which encodes estrogen receptor α, differs between the morphs and predicts both territorial and parental behavior. Variation in the regulatory regions of ESR1 causes an imbalance in expression of the two alleles, and the degree to which this imbalance favors the supergene allele predicts territorial singing. In heterozygotes, knockdown of ESR1 causes a phenotypic switch, from more aggressive to less aggressive. We recently showed that another gene important for social behavior, vasoactive intestinal peptide (VIP), is differentially expressed between the morphs and predicts territorial singing. We hypothesize that ESR1 and VIP contribute to behavior in a coordinated way and could represent co-adapted alleles. Because the supergene contains more than 1000 individual genes, this species provides rich possibilities for discovering alleles that work together to mediate life-history trade-offs and maximize the fitness of alternative complex phenotypes.

Adaptation to drought is coupled with slow growth, but independent from phenology in marginal silver fir (Abies alba Mill.) populations

Drought is one of the most important selection pressures for forest trees in the context of climate change. Yet, the different evolutionary mechanisms, and their environmental drivers, by which certain populations become more drought tolerant than others is still little understood. We studied adaptation to drought in 16 silver fir (Abies alba Mill.) populations from the French Mediterranean Alps by combining observations on seedlings from a greenhouse experiment (N = 8,199) and on adult tress in situ (N = 315). In the greenhouse, we followed half-sib families for four growing seasons for growth and phenology traits, and tested their water stress response in a "drought until death" experiment. Adult trees in the field were assessed for δ 13C, a proxy for water use efficiency, and genotyped at 357 SNP loci. SNP data was used to generate a null expectation for seedling trait divergence between populations in order to detect the signature of selection, and 31 environmental variables were used to identify the selective environment. We found that seedlings originating from populations with low soil water capacity grew more slowly, attained a smaller stature, and resisted water stress for a longer period of time in the greenhouse. Additionally, adult trees of these populations exhibited a higher water use efficiency as evidenced by their δ 13C. These results suggest a correlated evolution of the growth-drought tolerance trait complex. Population divergence in bud break phenology was adaptive only in the second growing season, and evolved independently from the growth-drought tolerance trait complex. Adaptive divergence in bud break phenology was principally driven by the inter- and intra-annual variation in temperature at the geographic origin of the population. Our results illustrate the different evolutionary strategies used by populations to cope with drought stress at the range limits across a highly heterogeneous landscape, and can be used to inform assisted migration programs.

No evidence for differential sociosexual behavior and space use in the color morphs of the European common wall lizard (Podarcis muralis)

Explaining the evolutionary origin and maintenance of color polymorphisms is a major challenge in evolutionary biology. Such polymorphisms are commonly thought to reflect the existence of alternative behavioral or life-history strategies under negative frequency-dependent selection. The European common wall lizard Podarcis muralis exhibits a striking ventral color polymorphism that has been intensely studied and is often assumed to reflect alternative reproductive strategies, similar to the iconic "rock-paper-scissors" system described in the North American lizard Uta stansburiana. However, available studies so far have ignored central aspects in the behavioral ecology of this species that are crucial to assess the existence of alternative reproductive strategies. Here, we try to fill this gap by studying the social behavior, space use, and reproductive performance of lizards showing different color morphs, both in a free-ranging population from the eastern Pyrenees and in ten experimental mesocosm enclosures. In the natural population, we found no differences between morphs in site fidelity, space use, or male-female spatial overlap. Likewise, color morph was irrelevant to sociosexual behavior, space use, and reproductive success within experimental enclosures. Our results contradict the commonly held hypothesis that P. muralis morphs reflect alternative behavioral strategies, and suggest that we should instead turn our attention to alternative functional explanations.