Evolution

When two populations come into secondary contact, assortative mating can act as a barrier to gene flow. However, when assortative mating is incomplete, associations between preference and cue loci can break down, eroding assortative mating and collapsing species boundaries. Here, we investigate factors that can contribute to the maintenance of differentiated reproductive populations despite gene flow, resulting in either clinal hybrid zones or overlap zones between distinct populations. By simulating secondary contact on a uniform two-dimensional landscape using an individual-based computer model, we examine how the potential outcomes of secondary contact depend on search costs and/or pleiotropy between various combinations of mating cue, mating preference, and hybrid viability traits. We find that search costs can maintain stable mating trait clines, even in the absence of hybrid inviability traits. Pleiotropy between mating cues and hybrid inviability (i.e., "magic cues") temporarily stabilizes the mating cue cline but fails to maintain mating preference differentiation, since preference is not itself under direct selection. As the preference cline collapses, assortative mating breaks down, and -- absent search costs or very strong hybrid inviability -- the mating cue cline subsequently collapses as well. Thus, without preference-cue or preference-viability pleiotropy, search costs are essential for maintaining mating trait differentiation. Changes to the genetic architecture of mating traits strongly alter the distribution of phenotypes, affecting whether parental populations can coexist spatially. Small parameter changes can abruptly shift outcomes, highlighting the sensitivity of these systems.

Understanding how divergence in reproductive traits can promote speciation remains a fundamental question in evolutionary biology. Sperm morphology and kinematics diverge rapidly across species. However, the effect of hybridization between recently diverged species on sperm traits remains unclear, limiting our understanding of how reproductive isolation evolves. Here, we evaluated sperm morphology, kinematics, and trait integration in species of a young (~10,000 years), sympatric Cyprinodon pupfish radiation from San Salvador Island, Bahamas, as well as fertile advanced-generation hybrids between two of these species. We found significant divergence in flagellum length, midpiece area, and sperm velocity among species. In contrast, hybrids displayed transgressive kinematic profiles defined by high velocities, reduced path curvature, and distinct patterns of sperm kinematic integration. Our findings suggest that hybridization between recently diverged species may reorganize the underlying control of sperm locomotor mechanisms, generating novel phenotypes that could contribute to reproductive isolation in the early stages of speciation.

A long-standing question in evolutionary biology is: Under what circumstances can speciation occur despite hybridisation or because of hybridisation? Some models of hybrid incompatibilities predict that speciation can occur even in the presence of gene flow and strong selection against hybrids, an outcome also influenced by genetic contributions from parental species and genetic architecture. On the other hand, empirical work has shown that heterosis can counteract the effect of incompatibilities, hindering the speciation process. Theoretical models that simultaneously consider the positive and negative impacts of hybridisation on fitness remain scarce, raising questions about the effects of hybrid incompatibilities in the presence of heterosis. To address this question, we study how (Bateson)-Dobzhansky-Muller incompatibilities (BDMIs) interact with overdominant mutations in a two-locus population genetics model of an isolated hybrid population. We find that the strength of overdominance relative to incompatibilities determines the long-term genetic composition of the hybrid population. We show that high recombination exposes incompatibilities to selection and reduces the frequency of derived alleles in the hybrid population, limiting the strength of BDMIs that can be maintained by the balance with overdominance. We also show how neutral variation is affected by the strength of selection and the recombination rate between incompatible loci, generating patterns that include an increase in local variation resembling associative overdominance, or a reduction resembling background selection. Such variation of neutral variation, particularly at intermediate distances from BDMI loci, can generate peaks or troughs of diversity that are explained by the recombination rate between BDMI loci, and initial proportions of admixture between parental populations. Our work demonstrates how the genomic conflict caused by the interplay of overdominance and hybrid incompatibilities, recombination, and parental contributions, shape the genome of an isolated hybrid population.

Across animals, variation in adult body size is accompanied by coordinated variation in the size of individual morphological traits. However, the same morphological trait can scale differently with body size depending on what drives the size variation. In Drosophila melanogaster, for example, wing size scales differently with body size when size varies because of developmental nutrition versus developmental temperature. Whether the genetic basis of size plasticity and scaling is shared across different environmental regulators of size remains unclear, but is central to predicting how selection acts on the developmental mechanisms that regulate trait size, plasticity and morphological scaling. Using ~200 isogenic D. melanogaster lineages, we measured wing and leg size across nutritional and thermal treatments. For each lineage, we estimated nutritional and thermal plasticity for both traits, as well as the wing-leg individual-level scaling relationship, or ILSR, generated by each environmental source of size variation. We found extensive genetic variation in both thermal and nutritional plasticity for wings and legs, and in the slope of the ILSR between them. However, a lineages thermal plasticity was genetically uncorrelated with its nutritional plasticity for either trait, and we detected no genetic correlation between the slopes of thermal and nutritional wing-leg ILSRs. We also found no genetic correlation in the slope of nutritional wing-leg ILSRs across temperatures. Thus, the slope of a lineages nutritional ILSR at 17{degrees}C was not predictive of its slope at 25{degrees}C of 28{degrees}C. Nevertheless, the overall pattern of nutritional ILSRs was conserved across temperatures. These results suggest that the genetic architecture of size plasticity and scaling depends on the environmental source of size variation. Consequently, the evolutionary response of scaling to selection in heterogeneous environments may not be predictable from genetic variation measured in any single environment.

Hybridization, the interbreeding between species or genetically distinct populations, can lead to deleterious fitness consequences, but simultaneously it can boost adaptive potential by increasing genetic variation, especially in novel environments. However, how incompatibilities and beneficial genetic combinations interplay across generations remains poorly understood. Here, we tracked how two fitness proxies, the absolute number of adult offspring and the proportion of eggs that reached adulthood, evolve in parental and hybrid populations at three different time points over 50 generations, in both novel and ancestral environments. To test this, we used two geographically distinct populations with low divergence (Dxy=0.002) of the two-spotted spider mite (Tetranychus urticae). Across the first three generations, hybrids showed significantly lower fitness than parental populations in the ancestral environment, indicating incompatibilities between the parental genomes. In contrast, hybrid and parental fitnesses were similar in the novel environment, indicating that the impact of incompatibilities was minor compared to the selection imposed by the novel environment. However, after 50 generations, hybrids displayed similar fitness relative to parental populations in all environments, suggesting resolution of the incompatibilities. Furthermore, around generation 45, hybrids temporarily outperformed parental populations in a novel environment, suggesting a transient window of higher adaptive potential, before fitness stabilized again by generation 50. In conclusion, we show that hybrid populations of T. urticae can swiftly purge incompatibilities when genetic divergence is low. These findings suggest that the dynamics of incompatibility resolution and adaptive potential of novel haplotypes play out over a long time frame, highlighting the importance of tracking hybrid fitness past the first few generations.

The evolutionary patterns of trait diversification provide insights into the function of the trait. Early burst of trait evolution is often associated with adaptive radiation, rapidly diversifying the trait in response to vacant niches followed by the slowdown of the diversification with niche filling, whereas late burst is more likely to be associated with sexual selection, possibly contributing to reproductive barriers between closely related species. Here, we studied the diversification of tail fork depth through time in hirundines to infer its function, which remains unclear due to the competing two alternative hypotheses: the sexual selection hypothesis, which is a classic explanation of deeply forked tails, proposed that this trait has evolved via sexual selection, which was then challenged by the viability selection hypothesis, which proposed that deeply forked tails have mainly evolved via viability selection for enhancing aerodynamic performance during aerial foraging on large prey. We found a late burst of tail fork depth, but not of bill length, i.e., an index of prey size. The observed pattern is consistent with the sexual selection hypothesis but not with the viability selection hypothesis.

Sensory systems mediate the interaction between organisms and their environment, but how complex sensory pathways evolve and relate to variation in perception and behavior across ecological contexts, remains poorly understood, especially for terrestrial taxa. Here, we investigate whole-visual-system adaptation in Heliconius erato butterflies. Using continent-wide sampling, we demonstrate that within H. erato, facet count significantly decreased with increasing elevation. Common-garden rearing of low-elevation H. erato populations from Ecuador and their high-elevation sister species, H. himera, showed that eye and brain morphology are heritable, and comparisons to genomic measures of divergence indicates that this variation is due to divergent selection. Parallel comparisons from Colombia involving H. chestertonii (high elevation) and H. erato venus (low elevation) further revealed that eye and brain morphology can evolve as independent, decoupled traits. For both locations, differences in visual acuity correlated with variation in facet count. We also observed parallel evolution of spectral sensitivity, with independent high-elevation populations having fewer red-reflecting lateral filtering pigments. To experimentally link visual system morphology to behavior, we assessed visual acuity in second-generation H. erato cyrbia-H. himera hybrids. Overall, acuity was influenced by facet count, and when analyzed together with brain morphology, by a positive interaction between facet count and optic lobe volume, demonstrating that structural investment in the eye and neural expansion combine to maximize visual perception. This work shows that visual adaptation is a multi-layered process whereby sensory traits can evolve independently under localized ecological pressures, but evolution across the visual pathway contributes to refinements in behavioral performance.

Recognizing patterns in functional trait evolution is a necessary step in testing macroevolutionary questions. Quantification of these patterns and interpretation of their generative processes relies on an ever-expanding suite of comparative approaches, but current methods oversimplify the process-to-pattern mapping. This simplification may promote binary classifications of patterns and their drivers, such as adaptive versus non-adaptive or constrained versus unconstrained. A potentially more robust and evolutionarily informative alternative is to fit an expanded suite of evolutionary models at different levels of taxonomic or ecological resolution to reveal how evolutionary drivers imprint simultaneously on observed patterns of trait evolution. Here, we perform a clade-wide analysis of functional trait evolution across Anura (frogs and toads). Focusing on three functionally important traits, we quantify the relative fit of random, directional, and bounded models at the intra- and inter-microhabitat levels to explore how changes in trait function across microhabitats translate into different evolutionary regimes. We recover heterogeneous regime-level distributions of model support that imply complex underlying evolutionary dynamics, while also revealing methodological biases and model identifiability issues. These findings underscore the need to develop more robust tools for evolutionary model fitting and advance beyond binary frameworks for interpreting evolutionary processes.

The evolutionary origins of asexuality remain poorly understood, despite extensive research on its ecological and evolutionary consequences. Asexuality often arises through hybridization between species with intermediate genomic divergence, implying that hybrid-induced asexuality may be partly repeatable. The Amazon molly (Poecilia formosa), the first asexual vertebrate known to science, challenges this view: repeated experimental crosses between its extant parental species have failed to recreate a stable Amazon molly-like lineage. This apparent paradox gave rise to the Rare Formation Hypothesis, which proposes that stable asexuality requires an exceptionally specific genomic combination. Here, we combine experimental crosses, molecular cytogenetics, and population genomics to test whether ancestral introgression before the hybrid speciation event set the stage for the singular origin of the Amazon molly. We show that most experimental hybrids are viable but sexual, but that a subset of F1 hybrids produce unreduced eggs through a mechanism distinct from that of the Amazon molly. Population genomic analyses reveal that introgression between parental species likely predated the formation of the Amazon molly, and shared homozygous tracts across Amazon molly genomes support inheritance from admixed progenitors. Together, our findings reconcile the repeatable and contingent views of the origin of asexuality, suggesting that ancestral introgression may be the missing mechanism assembling the rare genomic combinations required for seemingly unrepeatable evolutionary innovations, including the emergence of asexual species.

In many species with female sperm storage, ejaculates from multiple males overlap in the female reproductive tract, making sperm competitive ability a key component of male reproductive fitness and a target of rapid evolutionary change in the underlying genes. Here, we used controlled laboratory assays of Drosophila melanogaster sperm competition, with doubly-mated females and paternity assignment of offspring, to ask whether a Bradley-Terry framework can effectively summarize and predict competitive outcomes. The Bradley-Terry model is a probabilistic approach that estimates a latent "ability" score for each contestant based on outcomes of pairwise contests, and thus is naturally suited to data from sperm competition, which are intrinsically pairwise. We selected five distinct male genotypes: four carried strongly expressed RFP or GFP markers that allowed us to distinguish their heterozygous offspring under UV illumination, and the fifth was Canton-S, a standard wild-type genotype that served as our reference. Using Canton-S females, we assayed all 20 ordered pairwise combinations of first and second male, recorded successful double matings, and quantified the offspring sired by each male. We then extended the Bradley-Terry model to estimate genotype-specific competitive success separately for first-male "defense" (fertilization success following initial mating, also called "P1") and second-male "offense" (fertilization success following a remating, also called "P2"). This framework provides a flexible and efficient way to integrate results across large arrays of pairwise mating tests and to derive predictive scores for sperm competitive performance.

Understanding how populations persist in gradually deteriorating environments through evolution is a central question in ecology and evolutionary biology. Previous studies have primarily focused on identifying the critical rate of environmental change beyond which extinction is certain. However, the existence of a viable equilibrium when the rate is below the threshold does not guarantee that a population can survive the transient dynamics to reach it. Using a quantitative genetic model that explicitly incorporates feedback among population size, genetic variance, and mean trait evolution, we show that population persistence can exhibit bistability when the rate of environmental change is below the extinction threshold. Specifically, extinction still occurs if the initial population size and genetic variance fall below a critical level. The initial state also influences the eco-evolutionary dynamics, such that a temporary increase or decline in population size and/or genetic variance does not necessarily predict the ultimate fate of the population. Therefore, in addition to estimating the critical rate of environmental change for extinction, characterizing current population size, genetic variation, and the degree of maladaptation may improve predictions of extinction risk in deteriorating environments.

Evolutionary diversity in metabolic rate underlies differences in physiology, morphology, and life history across the tree of life. Cell size has been proposed as an important determinant of metabolic rate. The mechanisms underlying this proposed connection are based on the lower surface area to volume ratios in larger cells. As relative surface area decreases, the cost of maintaining ion gradients across the cell membrane through action of the Na+/K+-ATPase pump are posited to decrease, lowering overall metabolic costs. Despite strong theoretical support for this model, and its incorporation into broader models of life history evolution, empirical measurement of Na+/K+-ATPase activity in species that differ in cell size has been lacking. Here, we study nine species of salamanders of the genus Plethodon that span a large range of cell sizes approaching the animal upper limit. We compare basal cellular respiration rates, relative cost of the Na+/K+-ATPase pump, and maximal mitochondrial respiration rates in liver and heart tissue. Contrary to predictions, we find no support for a relationship between cell size and any of these mitochondrial respiratory variables. We reconcile this surprising result with broader phylogenetic studies showing a lack of correlation between cell size and metabolic rate at the organismal level.

Dominance is a widespread feature of genetic variants which affects life-history traits and fitness. Although dominance is generally thought to be an intrinsic property of genetic variants, it can sometimes evolve, as in the classic case of melanism in the peppered moth. The broader question of how likely dominance is to evolve is, however, controversial, because conditions favouring dominance evolution are often restrictive. Here, we revisit Haldanes classic hypothesis that dominance might evolve during the spread of beneficial mutations to fixation (i.e., during selective sweeps). We first confirm results of earlier models that sweeps of unconditionally beneficial mutations generate little potential for dominance to evolve, even in cases where modifier alleles segregate prior to selective sweeps. However, when sweeping beneficial alleles trade off between different environments -- which we explore with the illustrative case of sexually antagonistic selection -- the scope for dominance evolution expands. This occurs because modifier alleles can alter dominance separately in each environment, increasing the mean fitness of heterozygotes, prolonging the sojourn time of the sweep, and generating more heterozygosity upon which the modifier can act. In extreme cases, beneficial mutations that were initially destined for fixation can undergo a "dominance reversal" as a result of dominance evolution, converting them to balanced polymorphisms. We quantify how regularly dominance reversals of sweeping sexually antagonistic alleles can be expected to evolve. Overall, our results highlight conditions that allow the dominance of beneficial mutations to evolve, which we discuss in light of data on the frequency of selective sweeps, standing genetic variation for modifiers, and plasticity of modifier effects.

How broad phenotypic variation is maintained in natural populations in the face of selection is a central question in evolutionary biology. We address this question in the water strider Microvelia longipes, where males exhibit striking variation in rear leg length used in male-male contests for dominance. Using reaction norm experiments on inbred lines, we demonstrate that phenotypic plasticity contributes to expanding phenotypic variation, but requires high genetic variation to generate the broad range of trait expression observed in natural populations. Experimental evolution favouring trait exaggeration revealed that directional sexual selection not only fails to erode variation of male rear leg length, but rather amplifies it beyond the natural distribution. Additionally, male-limited selection in favour of dominance generated substantial fecundity costs in females, underscoring the role of sexual conflict driven by females in constraining exaggerated secondary sexual traits in males. Our findings show that sexually antagonistic selection and directional sexual selection jointly generate high genetic variation, which phenotypic plasticity inflates into broad phenotypic distribution of male weapon size. This provides an empirical explanation for the high variability of male exaggerated weapons in nature.

Hosts defend themselves against parasites through resistance (reducing parasite burden) and tolerance (reducing the fitness cost of infection without affecting parasites). This distinction has important evolutionary implications: resistance is predicted to maintain polymorphism while tolerance tends to fix, and only resistance is expected to provoke parasite counter-adaptation. The reaction-norm framework, which infers tolerance from the slope of host fitness regressed on parasite burden, assumes that a shallow slope reflects parasite-independent host protection. We test this assumption using a within-host model in two variants: microparasites (Model 1, with within-host replication) and macroparasites (Model 2, without). Sublethal immunity impairs the host-exploitation rate of the parasite, reducing both growth and per-parasite virulence without killing them. We show that this generates systematic slope differences among host genotypes that the framework interprets as variation in tolerance. Furthermore, the ranking of slopes across genotypes reverses between linear and sigmoidal damage functions: under linear damage, the strongest immune responder appears most tolerant; under sigmoidal damage, the weakest responder does. Decomposition of the damage reduction shows that virulence reduction accounts for the majority of the effect across both model variants. Thus, the reaction-norm slope cannot determine whether host fitness is maintained by parasite-independent tissue protection or by sublethal impairment of parasites.

The repeated evolution of certain complex traits within a given lineage is compelling and presents an opportunity for understanding selective and genetic features that promote rapid multi-trait adaptation. The North American plant genus Penstemon is one such example, with numerous evolutionary shifts from ancestral bee syndrome to hummingbird syndrome flowers. Here, we traced the evolution of multi-trait floral phenotypes within a focal clade of Penstemon using a new whole genome phylogenomic estimate and found four independent origins of the hummingbird syndrome. We found strong evolutionary convergence in floral traits across these four origins and observed that the hummingbird syndrome assembled rapidly, without leaving a signal of stepwise modification of traits. Patterns of evolutionary correlations among floral traits in bee syndrome Penstemon species likely enable rapid shifts to hummingbird syndrome. Using phylogenomic tests for introgression, we found no evidence that adaptive introgression has fueled the four repeated origins of hummingbird syndrome flowers. Patterns of allele sharing were instead consistent with substantial levels of incomplete lineage sorting, suggesting repeated complex adaptation involves de novo mutation or adaptation from ancestral variation.

Many endosymbioses in eukaryotes superficially appear to be beneficial to both participants. However, there is little direct evidence for this, and symbioses naturally set up conditions in which each member of the pair is under selection to extract resources from the other. Ultimately, the endosymbiont either evolves to be in conflict with the interests of the host or to act cooperatively with the host contrary to its own best interests. Focusing on obligate symbioses, we develop theory to clarify the population-genetic conditions favoring the alternative outcomes. The balance is usually tipped in favor of exploitation by the symbiont, particularly when the number of symbionts within host cells is high, selection is strong on symbionts relative to hosts, there is horizontal transfer of symbionts, and/or the symbionts have accelerated mutation rates or turnover times. If the symbiont conditions the host-cell biology to enhance within-host population sizes, selection for selfish symbionts will be further enhanced by the diminished level of within-host drift. Although the host evolves in parallel to exploit resources from the endosymbiont, the net result is often a stalemate in which the host is no better off than prior to host-symbiont coevolution. Strict vertical inheritance can result in an evolutionary alignment of interests of the endosymbiont and the host, as this minimizes the possibility of within-host selection, but even here there is a critical host population size below which the symbiont evolves to exploit the host. These results suggest that the evolutionary enslavement of a symbiont to benefit a host species requires a narrow mix of population-biological features of both participants.

Brain size and structure are hypothesized to be adaptively designed to satisfy the behavioral requirements of securing food and living socially. The importance of these socioecological and sociobiological selective forces in brain evolution is constantly debated. Socioecological divergence is striking in the Neotropical ant genus Neoponera: N. apicalis is a generalist solitary predator forming small colonies of ~100 whereas N. commutata colonies are approximately 10 times as large and workers pheromonally organize cooperatively raids only on Syntermes termite colonies. We interspecifically compared the size and structure of the compound eyes, size and number of antennal glomeruli, mosaic brain scaling and synaptic processing (microglomeruli-MG). Our results indicate that N. apicalis workers have a larger number of ommatidia, antennal lobe glomeruli, and allometrically larger antennal and optic lobes than N. commutata. These sensory traits were associated with differences in higher-order processing architectures in the mushroom body (MB) microglomeruli (MG). N. commutata workers had an allometrically larger MB, perhaps due to their socially complex chemical foraging communication, although MG density in N. apicalis was higher in both the MB lip and collar, regions associated with processing olfactory and visual information, respectively. The increase in MG density in N. apicalis may be associated with higher demands for navigation, learning, and memory, as well as a higher density of antennal lobe glomeruli to support prey odor discrimination. In contrast, N. commutata workers had larger ommatidia and antennal lobe glomeruli. Larger ommatidia correlate with their diurnal/nocturnal habits and a larger MB Our findings indicate that differences in behavioral performance demands associated with socioecological differentiation are reflected in variation in visual and olfactory system structure, brain size, mosaicism, and synaptic organization. Our results support both social and ecological brain hypothesis as drivers of mosaic brain evolution.

Reproductive success often depends on coordinated combinations of morphology, dispersal ability, and mating behavior. Supergenes, genomic regions of suppressed recombination, allow such combinations to be inherited by offspring as a single unit. In ants, independently evolved supergenes control colony queen number, yet few studies have investigated their joint influence on morphology, mating, and colony-founding in sexuals. Formica cinerea provides a unique opportunity to address this question because it harbors two supergenes that together produce three queen and male morphs: large monogyne, large polygyne, and small polygyne. Here we show that these supergenes jointly shape an integrated suite of traits across the reproductive cycle. Wing area was primarily associated with the chromosome 3 supergene, with monogyne individuals having larger wings than polygyne individuals. Thorax volume was associated with the chromosome 9 supergene, with small polygyne individuals having reduced thorax volume regardless of social origin. Mating was assortative for both supergenes in large morphs but random in small polygyne queens. Independent colony founding was almost exclusively performed by large monogyne queens; initial egg production was unaffected by mate genotype. These findings show that the two supergenes jointly coordinate dispersal morphology, mate choice, and colony-founding into coherent reproductive strategies, preventing maladaptive intermediate phenotypes.

Internal evolutionary conflicts arise when elements within an organism have diverging fitness interests. Examples range from meiotic drive and cytoplasmic male sterility to transposable elements and supernumerary B chromosomes. While once treated as genetic curiosities, they are now seen as widespread and major drivers of eukaryotic genome evolution. Yet their study remains fragmented, with no clear entry point not only for those who wish to gain an overview of theoretical advances, or those who wish to construct models of their own. Here, we discuss ways in which internal evolutionary conflicts have been modelled and develop a common population genetic framework for building such models. The framework provides explicit criteria for what counts as conflict, distinguishing it from fitness trade-offs, and formalises how and when internal conflicts arise. By treating different cases within the same structure, it shows that these diverse phenomena share a common logic.