How does natural selection lead to cooperation between competing individuals? The Prisoner's Dilemma captures the essence of this problem. Two players can either cooperate or defect. The payoff for mutual cooperation, R, is greater than the payoff for mutual defection, P. But a defector versus a cooperator receives the highest payoff, T, where as the cooperator obtains the lowest payoff, S. Hence, the Prisoner's Dilemma is defined by the payoff ranking T > R > P > S. In a well-mixed population, defectors always have a higher expected payoff than cooperators, and therefore natural selection favors defectors. The evolution of cooperation requires specific mechanisms. Here we discuss five mechanisms for the evolution of cooperation: direct reciprocity, indirect reciprocity, kin selection, group selection, and network reciprocity (or graph selection). Each mechanism leads to a transformation of the Prisoner's Dilemma payoff matrix. From the transformed matrices, we derive the fundamental conditions for the evolution of cooperation. The transformed matrices can be used in standard frameworks of evolutionary dynamics such as the replicator equation or stochastic processes of game dynamics in finite populations.

Many molecular phylogenies show longer root-to-tip path lengths in species-rich groups, encouraging hypotheses linking cladogenesis with accelerated molecular evolution. However, the pattern can also be caused by an artifact called the node density effect (NDE): this effect occurs when the method used to reconstruct a tree underestimates multiple hits that would have been revealed by extra nodes, leading to longer root-to-tip path lengths in clades with more terminal taxa. Here we use a twofold approach to demonstrate that maximum likelihood and Bayesian methods also suffer from the NDE known to affect parsimony. First, simulations deliberately mismatching the simulation and reconstruction models show that the greater the model disparity, the greater the gap between actual and reconstructed tree lengths, and the greater the NDE. Second, taxon sampling manipulation with empirical data shows that NDE can still be present when using optimized models: across 12 datasets, 70 out of 109 sister path comparisons showed significant evidence of NDE. Unless the model fairly accurately reconstructs the real tree length—and given the complexity of real sequence evolution this may be uncommon—it will consistently produce a node density artifact. At commonly encountered divergence levels, a 10% underestimation of tree length results in ≥ 80% of simulated phylogenies showing a positive NDE. Bayesian trees have a slight but consistently stronger effect. This pervasive methodological artifact increases apparent rate heterogeneity, and can compromise investigations of factors influencing molecular evolutionary rate that use path lengths in topologically asymmetric trees.

The formation of hybrid zones between nascent species is a widespread phenomenon. The evolutionary consequences of hybridization are influenced by numerous factors, including the action of natural selection on quantitative trait variation. Here we examine how the genetic basis of floral traits of two species of Louisiana Irises affects the extent of quantitative trait variation in their hybrids. Quantitative trait locus (QTL) mapping was used to assess the size (magnitude) of phenotypic effects of individual QTL, the degree to which QTL for different floral traits are colocalized, and the occurrence of mixed QTL effects. These aspects of quantitative genetic variation would be expected to influence (1) the number of genetic steps (in terms of QTL substitutions) separating the parental species phenotypes; (2) trait correlations; and (3) the potential for transgressive segregation in hybrid populations. Results indicate that some Louisiana Iris floral trait QTL have large effects and QTL for different traits tend to colocalize. Transgressive variation was observed for six of nine traits, despite the fact that mixed QTL effects influence few traits. Overall, our QTL results imply that the genetic basis of floral morphology and color traits might facilitate the maintenance of phenotypic divergence between Iris fulva and Iris brevicaulis, although a great deal of phenotypic variation was observed among hybrids.

Kin-recognition mechanisms allow helping behaviors to be directed preferentially toward related individuals, and could be expected to evolve in many cases. However, genetic kin recognition requires a genetic polymorphism on which recognition is based, and kin discriminating behaviors will affect the evolution of such polymorphism. It is unclear whether genetic polymorphisms used in kin recognition should be maintained by extrinsic selection pressures or not, as opposite conclusions have been reached by analytical one-locus models and simulations exploring different population structures. We analyze a two-locus model in a spatially subdivided population following the island model of dispersal between demes of finite size. We find that in the absence of mutation, selection eliminates polymorphism in most cases, except with extreme spatial structure and low recombination. With mutation, the population may reach a stable limit cycle over which both loci are polymorphic; however, the average frequency of conditional helping can be high only under strong structure and low recombination. Finally, we review evidence for extrinsic selection maintaining polymorphism on which kin recognition is based.

Species with elaborate parental care often also show intense sibling competition over resources provided by parents, suggesting joint evolution of these two traits. Despite this, the evolution of elaborate parental care and the evolution of intense sibling competition are often studied separately. Here, we examine the interaction between parental food provisioning and sibling competition for resources through the joint manipulation of the presence or absence of parents and brood size in a species with facultative parental care: the burying beetle Nicrophorus vespilloides. The effect of the interaction between the presence or absence of parents and brood size was strong; brood size had a strong effect on growth when parents provided care, but no effect when parents were absent. As in previous studies, offspring grew faster when parents were present than when parents were absent, and offspring grew faster in smaller broods than in larger broods. Our behavioral observations showed that brood size had a negative effect on both the amount of time parents spent providing resources to individual offspring and the offspring's effectiveness of begging, confirming that the level of sibling competition increased with brood size. Furthermore, offspring in larger broods shifted more from begging toward self-feeding as they grew older compared to offspring in small broods. Our study provides novel insights into the joint evolution of parental care and sibling competition, and the evolution of offspring begging signals. We discuss the implications of our results in light of recent theoretical work on the evolution of parental care, sibling competition, and offspring begging signals.

The level of defense against great spotted cuckoo (Clamator glandarius) parasitism in different European populations of magpie (Pica pica) depends on selection pressures due to parasitism and gene flow between populations, which suggests the existence of coevolutionary hot spots within a European metapopulation. A mosaic of coevolution is theoretically possible at small geographical scales and with strong gene flow, because, among other reasons, plots may differ in productivity (i.e., reproductive success of hosts in the absence of parasitism) and defensive genotypes theoretically should be more common in plots of high productivity. Here, we tested this prediction by exploring the relationship between parasitism rate, level of defense against parasitism (estimated as both rejection rate and the frequency of the 457bp microsatellite allele associated with foreign egg rejection in magpies), and some variables related to the productivity (average laying date, clutch size, and number of hatchlings per nest) of magpies breeding in different subpopulations. We found that both estimates of defensive ability (egg rejection rate and frequency of the 457bp allele) covaried significantly with between-plot differences in probability of parasitism, laying date, and number of hatchlings per nest. Moreover, the parasitism rate was larger in more productive plots. These results confirm the existence of a mosaic of coevolution at a very local geographical scale, and the association between laying date and number of hatchlings with variables related to defensive ability and the selection pressure arising from parasitism supports the prediction of coevolutionary gradients in relation to host productivity.

Models of mating-system evolution emphasize the importance of frequency-dependent interactions among mating partners. It is also known that outcross siring success and the selfing rate in self-compatible hermaphrodites can be density dependent. Here, we use array experiments to show that the mating system (i.e., the outcrossing rate) and the siring success of morphs with divergent sex allocation strategies are both density dependent and frequency dependent in androdioecious populations of the wind-pollinated, annual plant Mercurialis annua. In particular, the outcrossing rate is a decreasing function of the mean interplant distance, regulated by a negative exponential pollen fall-off curve. Our results indicate that pollen dispersed from a male inflorescence are over 60% more likely to sire outcrossed progeny than equivalent pollen dispersed from hermaphrodites, likely due to the fact that males, but not hermaphrodites, disperse their pollen from erect inflorescence stalks. Because of this difference, and because males of M. annua produce much more pollen than hermaphrodites, the presence of males in the experimental arrays reduced both the selfing rate and the outcross siring success of hermaphrodites. We use our results to infer a density threshold below which males are unable to persist with hermaphrodites but above which they can invade hermaphroditic populations. We discuss our findings in the context of a metapopulation model, in which males can only persist in well-established populations but are excluded from small, sparse populations, for example, in the early stages of colonization.

Two fundamental questions dealing with simultaneous hermaphrodites are how resources are optimally allocated to the male and female function and what conditions determine shifts in optimal sex allocation with age or size. In this study, I explored multiple factors that theoretically affect fitness gain curves (that depict the relationship between sex-specific investment and fitness gains) to predict and test the overall and size-dependent sex allocation in a simultaneously hermaphroditic brooding shrimp with an early male phase. In Lysmata wurdemanni, sperm competition is absent as hermaphrodites reproducing in the female role invariably mated only once with a single other shrimp. Shrimps acting as females preferred small over large shrimps as male mating partners, male mating ability was greater for small compared to large hermaphrodites, and adolescent males were predominant in the population during the breeding season. In addition, brooding constraints were not severe and varied linearly with body size whereas the ability to acquire resources increased markedly with body size. Using sex allocation theory as a framework, the findings above permitted to infer the shape of the male and female fitness gain curves for the hermaphrodites. The absence of sperm competition and the almost unconstrained brooding capacity imply that both curves saturate, however the male curve levels off much more quickly than the female curve with increasing level of investment. In turn, the predominance of adolescent males in the population implies that the absolute gain of the female curve is greater than that of the male curve. Last, the size-dependent female preference and male mating ability of hermaphrodites determines that the absolute gain of the male curve is greater for small than for large hermaphrodites. Taking into consideration the inferred shape of the fitness gain curves, two predictions with respect to the optimal sex allocation were formulated. First, overall sex allocation should be female biased; it permits hermaphrodites to profit from the female function that provides a greater fitness return than the male function. Second, sex allocation should be size-dependent with smaller hermaphrodites allocating more than proportionally resources to male reproduction than larger ones. This size-dependent sex allocation permits hermaphrodites to profit from male mating opportunities that are the greatest at small body sizes. Size-dependent sex allocation is also expected because the male fitness gain curve decelerates more quickly than the female gain curve and experiments indicated that resources are greater for large than small hermaphrodites. These two predictions were tested when determining the sex allocation of hermaphrodites by dissecting their gonad and quantifying ovaries versus testes mass. Supporting the predictions above, hermaphrodites allocated, on average, 118 times more to the female than to the male gonad and the proportion of resources devoted to male function was higher in small than in large hermaphrodites. A trade-off between male and female allocation is assumed by theory but no negative correlation between male and female reproductive investment was observed. In L. wurdemanni, the relationship between sex-specific investment and fitness changes during ontogeny in a way that is consistent with an adjustment of sex allocation to improve size-specific reproductive success.

Genetic and developmental constraints have often been invoked to explain patterns of existing morphologies. Yet, empirical tests addressing this issue directly are still scarce. We here set out to investigate the importance of maternal body size as an evolutionary constraint on egg size in the tropical butterfly Bicyclus anynana, employing an artificial two-trait selection experiment on simultaneous changes in body and egg size (synergistic and antagonistic selection). Selection on maternal body size and egg size was successful in both the synergistic and the antagonistic selection direction. Yet, responses to selection and realized heritabilities varied across selection regimes: the most extreme values for pupal mass were found in the synergistic selection directions, whereas in the antagonistic selection direction realized heritabilities were low and nonsignificant in three of four cases. In contrast, for egg size the highest values were obtained in the lines selected for low pupal mass. Thus, selection on body size yielded a stronger correlated response in egg size than vice versa, which is likely to bias (i.e., constrain), if weakly, evolutionary change in body size. However, it seems questionable whether this will prevent evolution toward novel phenotypes, given enough time and that natural selection is strong. Correlated responses to selection were overall weak. Egg and larval development times tended to be associated with changes in maternal size, whereas variation in pupal development times weakly tended to follow variation in egg size. Lifetime fecundity was similar across selection regimes, except for females simultaneously selected for large body mass and small egg size, exhibiting increased fecundity. Multiple regressions showed that lifetime fecundity and concomitantly reproductive investment were primarily determined by longevity, as expected for an income breeder, whereas egg size was primarily determined by pupal mass. Evidence for a phenotypic trade-off between egg size and number was weak.

We study the evolution of the dispersal rate in a metapopulation model with extinction and colonization dynamics, akin to the model as originally described by Levins. To do so we extend the metapopulation model with a description of the within patch dynamics. By means of a separation of time scales we analytically derive a fitness expression from first principles for this model. The fitness function can be written as an inclusive fitness equation (Hamilton's rule). By recasting this equation in a form that emphasizes the effects of competition we show the effect of the local competition and the local population size on the evolution of dispersal. We find that the evolution of dispersal cannot be easily interpreted in terms of avoidance of kin competition, but rather that increased dispersal reduces the competitive ability. Our model also yields a testable prediction in term of relatedness and life-history parameters.

In 1996, Schluter showed that the direction of morphological divergence of closely related species is biased toward the line of least genetic resistance, represented by g_max, the leading eigenvector of the matrix of genetic variance–covariance (the G-matrix). G is used to predict the direction of evolutionary change in natural populations. However, this usage requires that G is sufficiently constant over time to have enough predictive significance. Here, we explore the alternative explanation that G can evolve due to gene flow to conform to the direction of divergence between incipient species. We use computer simulations in a mainland–island migration model with stabilizing selection on two quantitative traits. We show that a high level of gene flow from a mainland population is required to significantly affect the orientation of the G-matrix in an island population. The changes caused by the introgression of the mainland alleles into the island population affect all aspects of the shape of G (size, eccentricity, and orientation) and lead to the alignment of g_max with the line of divergence between the two populations' phenotypic optima. Those changes decrease with increased correlation in mutational effects and with a correlated selection. Our results suggest that high migration rates, such as those often seen at the intraspecific level, will substantially affect the shape and orientation of G, whereas low migration (e.g., at the interspecific level) is unlikely to substantially affect the evolution of G.

Phylogeographic studies designed to estimate rates and patterns of genetic differentiation within species often reveal unexpected and graphically striking cases of allele or haplotype sharing between species (introgression) via hybridization and backcrossing. Does introgression between species significantly influence population genetic structure relative to more conventional sources of differentiation (drift) and similarity (dispersal) among populations within species? Here we use mtDNA sequences from four species in two genera of sea urchins and sea stars to quantify the relative magnitude of gene flow across oceans and across species boundaries in the context of the trans-Arctic interchange of marine organisms between the Pacific and Atlantic oceans. In spite of the much smaller distances between sympatric congeners, rates of gene flow between sympatric species via heterospecific gamete interactions were small and significantly lower than gene flow across oceans via dispersal of planktonic larvae. We conclude that, in these cases at least, larvae are more effective than gametes as vectors of gene flow.

Plants can achieve an appropriate phenotype in particular conditions either constitutively or plastically, depending in part on the grain size of the environmental conditions being considered. Coarse-grained environmental variation should result in selection for local adaptation and no selection on plasticity to novel levels of the coarse-grained environmental factors. We tested the hypotheses that natural populations of the well-studied model system Arabidopsis thaliana are locally adapted to spatially coarse-grained environmental variation, and that the photoperiodic regime per se is at least partially responsible for that local adaptation, by exposing natural populations to photoperiodic regimes characteristic of their native and foreign (novel) environments. We also tested the hypothesis that plasticity to novel photoperiodic regimes should appear random. We found that populations showed evidence of local adaptation at a spatially coarse grain, although not to photoperiodic regime per se. We also found that the plasticities to novel photoperiodic regimes appeared random and did not generally show evidence of adaptive divergence. Our study highlights the need for caution in extrapolating from the finding of local adaptation to the causes of local adaptation.

Ecological character displacement takes place when two closely related species co-occur in only part of their geographical range, and selection to minimize competition between them promotes divergence in resource-use traits in sympatry but not in allopatry. Because populations sympatric with the heterospecific competitor will experience a different competitive environment than conspecific populations in allopatry, conspecific populations from these two competitive environments will also diverge in resource traits as an indirect consequence of interspecific ecological character displacement. Ultimately, ecologically dependent postmating isolation may arise between conspecific populations from these divergent competitive environments if offspring produced by matings between them are competitively inferior in either type of competitive environment. Yet, there are no direct tests of character displacement's role in initiating such postmating isolation. Here, we present a test by comparing the phenotypes and performances of spadefoot toad tadpoles produced from between-competitive-environment (BCE) matings versus those produced from within-competitive-environment (WCE) matings. When raised with naturally occurring competitors, BCE offspring grew significantly less and were significantly smaller than WCE offspring. BCE offspring generally performed worse even when raised alone, suggesting that they may have harbored intrinsic genetic incompatibilities. Moreover, the difference in growth and body size of BCE versus WCE offspring was significantly greater when each was raised with competitors than when each was raised alone, suggesting that BCE tadpoles were competitively inferior to WCE tadpoles. Presumably, this enhanced difference arose because BCE tadpoles produced an intermediate resource-use phenotype that is less well adapted to either competitive environment. Because larval size is under strong, positive, directional selection, reduced growth and size of BCE offspring may diminish gene flow between populations in divergent competitive environments, thereby generating postmating isolation. Thus, postmating isolation between conspecific populations, and possibly even speciation, may arise as a by-product of interactions between species.

Given the substantial costs of plant defenses against pathogens, there should be corresponding benefits that prevent resistance from being lost in natural plant populations. Here, we present evidence that systemic acquired resistance (SAR) benefits plants attacked by pathogenic bacteria in nature. In a large field experiment, we found that Arabidopsis thaliana treated with salicylic acid exhibited reduced titers of bacteria in their leaves and elevated fitness relative to controls. Most common members of the culturable bacterial community suffered this decrease, consistent with the role of SAR as a broad spectrum defense. We found no evidence of negative interactions between SAR and jasmonate-dependent resistance. Plants treated with jasmonic acid received significantly lower insect damage to their siliques, but exhibited no differences in bacterial growth or fitness relative to controls. Collectively, these data suggest a likely role of pathogenic bacteria in the maintenance of SAR, but not jasmonate-dependent resistance, in nature.