We develop and analyze an explicit multilocus genetic model of coevolution. We assume that interactions between two species (mutualists, competitors, or victim and exploiter) are mediated by a pair of additive quantitative traits that are also subject to direct stabilizing selection toward intermediate optima. Using a weak-selection approximation, we derive analytical results for a symmetric case with equal locus effects and no mutation, and we complement these results by numerical simulations of more general cases. We show that mutualistic and competitive interactions always result in coevolution toward a stable equilibrium with no more than one polymorphic locus per species. Victim-exploiter interactions can lead to different dynamic regimes including evolution toward stable equilibria, cycles, and chaos. At equilibrium, the victim is often characterized by a very large genetic variance, whereas the exploiter is polymorphic in no more than one locus. Compared to related one-locus or quantitative genetic models, the multilocus model exhibits two major new properties. First, the equilibrium structure is considerably more complex. We derive detailed conditions for the existence and stability of various classes of equilibria and demonstrate the possibility of multiple simultaneously stable states. Second, the genetic variances change dynamically, which in turn significantly affects the dynamics of the mean trait values. In particular, the dynamics tend to be destabilized by an increase in the number of loci.

To understand selection on recombination, we need to consider how linkage disequilibria develop and how recombination alters these disequilibria. Any factor that affects the development of disequilibria, including nonrandom mating, can potentially change selection on recombination. Assortative mating is known to affect linkage disequilibria but its effects on the evolution of recombination have not been previously studied. Given that assortative mating for fitness can arise indirectly via a number of biologically realistic scenarios, it is plausible that weak assortative mating occurs across a diverse set of taxa. Using a modifier model, we examine how assortative mating for fitness affects the evolution of recombination under two evolutionary scenarios: selective sweeps and mutation-selection balance. We find there is no net effect of assortative mating during a selective sweep. In contrast, assortative mating could have a large effect on recombination when deleterious alleles are maintained at mutation-selection balance but only if assortative mating is sufficiently strong. Upon considering reasonable values for the number of loci affecting fitness components, the strength of selection, and the mutation rate, we conclude that the correlation in fitness between mates is unlikely to be sufficiently high for assortative mating to affect the evolution of recombination in most species.

Character displacement, in which coevolution of similar species alters their phenotypes, can be difficult to identify on the basis of observational data alone. In two-species systems, the most commonly identified (i.e., classic) resulting pattern is greater phenotypic difference between species in sympatry than allopatry. We show that restricting studies to this pattern may exclude many instances of character displacement, particularly in the presence of spatial environmental gradients. We present four spatial models of character displacement in quantitative traits affecting competition and hybridization between the species. Our models highlight the connections between range limits and character displacement in continuous space. We conclude that the classic pattern is less likely to occur for a trait affecting resource acquisition than for a trait affecting mate choice. We also show that interspecific hybridization (when hybrids are inviable), even in very small amounts, has marked effects on the shape and stability of borders between species and the nature of character displacement. A survey of the empirical literature shows that character displacement studies often lack analysis of spatial phenotype and abundance data. We recommend more careful spatial sampling in character displacement studies, and we illustrate how comparison of clines in mean phenotype in sympatry and allopatry can be used to suggest the action of character displacement.

Humans and animals often become coinfected with pathogen strains that differ in virulence. The ensuing interaction between these strains can, in theory, be a major determinant of the direction of selection on virulence genes in pathogen populations. Many mathematical analyses of this assume that virulent pathogen lineages have a competitive advantage within coinfected hosts and thus predict that pathogens will evolve to become more virulent where genetically diverse infections are common. Although the implications of these studies are relevant to both fundamental biology and medical science, direct empirical tests for relationships between virulence and competitive ability are lacking. Here we use newly developed strain-specific real-time quantitative polymerase chain reaction protocols to determine the pairwise competitiveness of genetically divergent Plasmodium chabaudi clones that represent a wide range of innate virulences in their rodent host. We found that even against their background of widely varying genotypic and antigenic properties, virulent clones had a competitive advantage in the acute phase of mixed infections. The more virulent a clone was relative to its competitor, the less it suffered from competition. This result confirms our earlier work with parasite lines derived from a single clonal lineage by serial passage and supports the virulence-competitive ability assumption of many theoretical models. To the extent that our rodent model captures the essence of the natural history of malaria parasites, public health interventions which reduce the incidence of mixed malaria infections should have beneficial consequences by reducing the selection for high virulence.

Multilocus interactions (also known as Dobzhansky-Muller incompatibilities) are thought to be the major source of hybrid inviability and sterility. Because cytoplasmic and nuclear genomes have conflicting evolutionary interests and are often highly coevolved, cytonuclear incompatibilities may be among the first to develop in incipient species. Here, we report the discovery of cytoplasm-dependent anther sterility in hybrids between closely related Mimulus species, outcrossing M. guttatus and selfing M. nasutus. A novel pollenless anther phenotype was observed in F₂ hybrids with the M. guttatus cytoplasm (F_2G) but not in the reciprocal F_2N hybrids, F₁ hybrids or parental genotypes. The pattern of phenotypic segregation in the F_2G hybrids and two backcross populations fit a Mendelian single-locus recessive model, allowing us to map the underlying nuclear locus to a small region on LG7 of the Mimulus linkage map. Anther sterility was associated with a 20% reduction in flower size in backcross hybrids and we mapped a major cytoplasm-dependent corolla width QTL with its peak at the anther sterility locus. We argue that the cytonuclear anther sterility seen in hybrids reflects the presence of a cryptic cytoplasmic male sterility (CMS) and restorer system within the hermaphroditic M. guttatus population and therefore name the anther sterility locus restorer-of-male-fertility (RMF). The genetic mapping of RMF is a first step toward testing hypotheses about the molecular basis, individual fitness consequences, and ecological context of CMS and restoration in a system without stable CMS-restorer polymorphism (i.e., gynodioecy). The discovery of cryptic CMS in a hermaphroditic wildflower further suggests that selfish cytoplasmic evolution may play an important, but often undetected, role in shaping patterns of hybrid incompatibility and interspecific introgression in plants.

Electron transport system (ETS) function in mitochondria is essential for the aerobic production of energy. Because ETS function requires extensive interactions between mitochondrial and nuclear gene products, coadaptation between mitochondrial and nuclear genomes may evolve within populations. Hybridization between allopatric populations may then expose functional incompatibilities between genomes that have not coevolved. The intertidal copepod Tigriopus californicus has high levels of nucleotide divergence among populations at mitochondrial loci and suffers F₂ hybrid breakdown in interpopulation hybrids. We hypothesize that hybridization results in incompatibilities among subunits in ETS enzyme complexes and that these incompatibilities result in diminished mitochondrial function and fitness. To test this hypothesis, we measured fitness, mitochondrial function, and ETS enzyme activity in inbred recombinant hybrid lines of Tigriopus californicus. We found that (1) both fitness and mitochondrial function are reduced in hybrid lines, (2) only those ETS enzymes with both nuclear and mitochondrial subunits show a loss of activity in hybrid lines, and (3) positive relationships exist between ETS enzyme activity and mitochondrial function and between mitochondrial function and fitness. We also present evidence that hybrid lines harboring mitochondrial DNA (mtDNA) and mitochondrial RNA polymerase (mtRPOL) from the same parental source population have higher fitness than those with mtDNA and mtRPOL from different populations, suggesting that mitochondrial gene regulation may play a role in disruption of mitochondrial performance and fitness of hybrids. These results suggest that disruption of coadaptation between nuclear and mitochondrial genes contributes to the phenomenon of hybrid breakdown.

Selective history is thought to constrain the extent and direction of future adaptation by limiting access to genotypes that are advantageous in a novel environment. Populations of Chlamydomonas previously selected at high CO₂ were either backselected at ambient levels of CO₂, or selected at levels of CO₂ that last occurred during glaciation in the Pleistocene. There was no effect of selective history on adaptation to either level of CO₂, and the high CO₂ phenotypes were evolutionarily reversible such that fitness in ambient CO₂ returned to values seen in controls. CO₂ uptake affinity improved relative to the ancestor in both ambient and glacial CO₂, although wild-type regulation of CO₂ uptake, which deteriorated during previous selection at high CO₂, was not restored by selection at lower levels of CO₂. Trade-offs in both CO₂ uptake affinity and growth were seen after selection at any given level of CO₂. Adaptation to ambient and glacial-era levels of CO₂ produced a range of phenotypes, suggesting that chance rather than selective history contributes to the divergence of replicate populations in this system.

While evolutionary trends have long received much attention and have been widely disputed, new methods are now allowing the testing of directional hypotheses with increased rigor. Here, we test a general hypothesis about the way many kinds of discrete characters are thought to evolve, termed oligomerization. This is the tendency for serial structures (such as arthropod body and appendage segments) or armature (such as spines) to evolve primarily through loss and fusion. Focusing on the Crustacea, we use maximum likelihood methods to test for directional evolution in a large sample (> 500) of discrete traits, analyzed against molecular-based phylogenies. We find evidence for a significant trend toward trait loss, in accordance with the reduction principle. However, this trend is far from ubiquitous, with many characters exhibiting a reconstructed bias toward gains. These results suggest that caution must be used before drawing conclusions about which taxa are “primitive” or about the directionality of morphological shifts in the absence of phylogenetic analysis. Nevertheless, oligomerization—as a trend rather than a law—may be an important process that influences evolutionary trajectories from both morphological and functional perspectives.

Asexual reproduction is thought to doom organisms to extinction due to mutation accumulation and parasite exploitation. Theoretical models suggest that parthenogens may escape the negative effects of conspecifics and biological enemies through escape in space. Through intensive sequencing of a mitochondrial DNA (mtDNA) and a nuclear intron locus in sexual and parthenogenetic freshwater snails (Campeloma), I examine three questions: (1) Are sexual mtDNA lineages more restricted geographically than parthenogenetic mtDNA lineages? (2) Are independent parthenogenetic lineages shorter lived than sexual lineages? and (3) Do parthenogens have higher intraindividual nuclear sequence diversity and form well-differentiated monophyletic groups as expected under the Meselson effect? Geographic ranges of parthenogenetic lineages are significantly larger than geographic ranges of sexual lineages. Based on coalescence times under different demographic assumptions, asexual lineages are short lived, but there is variation in clonal ages. Although alternative explanations exist, these results suggest that asexual lineages may persist in the short term through dispersal, and that various constraints may cause geographic restriction of sexual lineages. Both allotriploid and diploid Campeloma parthenogens have significantly higher allelic divergence within individuals, but show limited nuclear sequence divergence from sexual ancestors. In contrast to previous allozyme evidence for nonhybrid origins of diploid Campeloma parthenogens, cryptic hybridization may account for elevated heterozygosity.

Much of sexual selection theory depends on assumptions about the genetic basis of variation in male mating success and sperm competitive ability. Despite intense interest in this topic, few genes have been identified that contribute to variation in these traits. Here we report the results of quantitative trait locus (QTL) analyses of mating success of male Drosophila melanogaster when exposed to virgin females, remating success of males with previously mated females, and both defense and offense components of sperm competition. We found two to four significant QTLs for remating success, but no QTLs for mating success, even though mating success was more genetically variable than remating success in the recombinant inbred lines used in this study. By combining these results with data from previous gene-expression experiments, we were able to identify three X-linked candidate genes for variation in remating ability. For two of these genes, QTL and expression data were completely concordant with respect to directionality of effects: high mating success was associated with high levels of gene expression and with beneficial QTL effects on the trait. We found equivocal evidence for genetic variation in sperm offense and defense in the recombinant inbred lines, and we did not find any significant QTLs for either sperm competition trait.

Starvation resistance is a trait often associated with longevity. Animals with increased longevity frequently show elevated starvation resistance and vice versa. Consequently, both life-history traits are thought to share genetic and physiological mechanisms, such as increased fat content and lowered metabolic rate. Here, we present results from 20 generations of selection on Drosophila melanogaster for increased starvation resistance at the time of adult eclosion. We observe that starvation resistance can be the result of more than one mechanism, all associated with an increase in fat resources. In general, metabolic rate is lowered under starved conditions relative to fed conditions. Metabolic rate in the starvation resistant lines is generally higher than in control lines under starved conditions. Starvation resistant flies are able to sustain a higher metabolic rate for a longer period of time when food is unavailable. This implies depletion of the increased fat reserves. However, longevity was not consistently affected by selection for increased starvation resistance. Similarly, paraquat resistance differed between selection lines and did not associate with starvation resistance, but rather with longevity. The results are discussed in relation to previous reported results on starvation resistance and its relation with mechanisms of aging and longevity.

The role of balancing selection in maintaining genetic variation for fitness is largely unresolved. This reflects the inherent difficulty in distinguishing between models of recurrent mutation versus selection, which produce similar patterns of inbreeding depression, as well as the limitations of testing such hypotheses when fitness variation is averaged across the genome. Signatures of X-linked overdominant selection are less likely to be obscured by mutational variation because X-linked mutations are rapidly eliminated by purifying selection in males. Although models maintaining genetic variation for fitness are not necessarily mutually exclusive, a series of predictions for identifying X-linked overdominant selection can be used to separate its contribution from other underlying processes. We consider the role of overdominant selection in maintaining fitness variation in a sample of 12 X chromosomes from a population of Drosophila melanogaster. Substantial variation was observed for male reproductive success and female fecundity, with heterozygous-X genotypes exhibiting the greatest degree of variance, a finding that agrees well with predictions of the overdominance model. The importance of X-linked overdominant selection is discussed along with models of recurrent mutation and sexually antagonistic selection.

Species level phylogenetic hypotheses can be used to explore patterns of divergence and speciation. In the tropics, speciation is commonly attributed to either vicariance, perhaps within climate-induced forest refugia, or ecological speciation caused by niche adaptation. Mimetic butterflies have been used to identify forest refugia as well as in studies of ecological speciation, so they are ideal for discriminating between these two models. The genus Ithomia contains 24 species of warningly colored mimetic butterflies found in South and Central America, and here we use a phylogenetic hypothesis based on seven genes for 23 species to investigate speciation in this group. The history of wing color pattern evolution in the genus was reconstructed using both parsimony and likelihood. The ancestral pattern for the group was almost certainly a transparent butterfly, and there is strong evidence for convergent evolution due to mimicry. A punctuationist model of pattern evolution was a significantly better fit to the data than a gradualist model, demonstrating that pattern changes above the species level were associated with cladogenesis and supporting a model of ecological speciation driven by mimicry adaptation. However, there was only one case of sister species unambiguously differing in pattern, suggesting that some recent speciation events have occurred without pattern shifts. The pattern of geographic overlap between clades over time shows that closely related species are mostly sympatric or, in one case, parapatric. This is consistent with modes of speciation with ongoing gene flow, although rapid range changes following allopatric speciation could give a similar pattern. Patterns of lineage accumulation through time differed significantly from that expected at random, and show that most of the extant species were present by the beginning of the Pleistocene at the latest. Hence Pleistocene refugia are unlikely to have played a major role in Ithomia diversification.

The Branchiosauridae was a clade of small amphibians from the Permo-Carboniferous with an overall salamander-like appearance. The clade is distinguished by an extraordinary fossil record that comprises hundreds of well-preserved specimens, representing a wide range of ontogenetic stages. Branchiosaurids had external gills and weakly ossified skeletons, and due to this larval appearance their status as neotenic (perennibranchiate) forms has long been accepted. Despite their extensive fossil record large specimens with an adult morphology appeared to be lacking altogether, but recently two adult specimens were identified in a rich sample of Apateon gracilis collected in the 19th century from a locality near Dresden, Saxony. These specimens are unique among branchiosaurids in showing a high level of ossification, including bones that have never been reported in a branchiosaur. These highlight the successive formation of features believed to indicate terrestrial locomotion, as well as feeding on larger prey items. Moreover, these transformations occurred in a small time window (whereas the degree of size increase is used as a proxy of time) and the degree of concentration of developmental events in branchiosaurids is unique among tetrapods outside the lissamphibians. These specimens are compared with large adults of the neotenic branchiosaurid Apateon caducus from the Saar-Nahe Basin, which despite their larger body size lack the features found in the adult A. gracilis specimens. These specimens give new insight into patterns of metamorphosis (morphological transformation) in branchiosaurids that are believed to be correlated to a change of habitat, and clearly show that different life-history pathways comparable to those of modern salamanders were already established in this Paleozoic clade.

Hybridization can be an evolutionary creative force by forming new polyploid species, creating novel genetic variation or acting as conduits of potentially advantageous traits between hybridizing forms. Evidence for the latter is often difficult to find because alleles under positive selection can spread rapidly across a hybrid zone and sweep to fixation. In Western Panama, an avian hybrid zone between two species of manakins in the genus Manacus exists where the unidirectional introgression of bright, yellow plumage into a white population provides evidence for the importance of hybrid zones as conduits of advantageous traits. Several lines of indirect evidence suggest that sexual selection favoring yellow plumage drives this asymmetrical spread, but more direct evidence is lacking. Along the edge of the hybrid zone, both yellow- and white-collared manakins are found in the same mating arenas or leks and compete for the same females (“mixed leks”), providing us with a unique opportunity to understand the dynamics of yellow plumage introgression. We studied these mixed leks to determine whether yellow males have a mating advantage over white males and, if so, whether the mating advantage is driven by male-male interactions, female choice, or both. We found that yellow males mated more than white males, suggesting that sexual selection favoring yellow males can, indeed, explain the spread of yellow plumage. However, we found that this advantage occurred only in mixed leks where the frequency of yellow males is greater than white males. This suggests that the advantage of yellow males may depend on the presence of other yellow males, which may slow the rate of introgression in leks where yellow frequency is low such as in areas where yellow males are beginning to colonize the white population. This, along with the geographic barrier posed by major rivers in the hybrid zone, may initially limit or slow the spread of yellow plumage. Finally, we found that yellow and white males were similar in aggression and body size, and held comparable positions within leks. Because these traits or factors are often important in or dictated by aggressive male-male interactions, these comparisons indicate that male-male interaction is not the primary mechanism for the spread of yellow plumage. However, white and yellow males received similar numbers of courtship visits from females but differed in the number of matings, suggesting that females actively rejected white in favor of yellow males. Our results indicate that sexual selection by female choice has driven the unidirectional introgression of yellow plumage into the white population, providing a mechanism for how hybrid zones act as conduits of novel and advantageous traits.

Where sperm competition occurs, the number and quality of sperm males inseminate relative to rival males influences fertilization success. The number of sperm males produce, however, is limited, and theoretically males should allocate sperm according to the probability of gaining future reproductive opportunities and the reproductive benefits associated with copulations. However, the reproductive opportunities and value of copulations males obtain can change over their lifetime, but whether individuals respond to such changes by adjusting the way they allocate sperm is unclear. Here we show that, in the fowl, Gallus gallus, dominant males, which have preferential access to females, modulate the number of sperm they ejaculate according to the availability of females. When presented with two females, dominant males allocated more sperm to higher quality females, whereas when females were on their own, only copulation order had an affect on their sperm numbers. In contrast, subordinate males, whose mating activity is restricted by dominant males, allocated high numbers of sperm to initial copulations, irrespective of female availability. We further show, by manipulating male social status, that sperm allocation is both phenotypically plastic, with males adjusting their patterns of sperm allocation according to their dominance rank, and intrinsic, with males being consistently different in the way they allocate sperm, once the effects of social status are taken into account. This study suggests that males have evolved sophisticated patterns of sperm allocation to respond to frequent fluctuations in the value and frequency of reproductive opportunities.

We tested whether selective breeding for early-age high voluntary exercise behavior over 16 generations caused the evolution of lifelong exercise behavior, life expectancy, and age-specific mortality in house mice (Mus domesticus). Sixteenth-generation mice from four replicate selection lines and four replicate random-bred control lines were individually housed from weaning through death and divided between two activity treatments (either with or without running wheels). Thus, there were four treatment groups: selection versus control crossed with active versus sedentary. The effects of selective breeding on life expectancy and age-specific mortality differed between females and males. In females, sedentary selection mice had early and high initial adult mortality and thus the lowest increases in mortality with age. Active selection females had the lowest early adult mortality, had limited mortality during midlife, and exhibited rapid increases in mortality rates at the very end of life; thus, they had deferred senescence. Median life expectancy was greater for both groups of selection females than for the two complementary groups of control females. Like females, sedentary selection males had the highest early adult mortality, and slow but steadily increasing mortality over the entire lifetime. Unlike the active selection females, active control males had the lowest mortality across the lifespan (until the end of life). Interestingly, the males with the lowest median life expectancy were those in the active selection treatment group. In both sexes, running (km/week) decreased over the lifetime to very low and virtually equivalent levels at the end of life in control and selection mice. Overall, these results demonstrate an evolutionary cost of selective breeding for males, regardless of exercise level, but a benefit for females when they have an outlet for the up-selected behavior. We conclude that correlated evolution of senescence occurs in mice selectively bred for high voluntary wheel running; exercise per se is beneficial for control mice of both sexes, but the impact on the effect of selection depends on sex; and the behavioral effect of exercise selection at an early age declines throughout the life span, which demonstrates decreasing genetic correlations over age for the genes involved in increased exercise.

Statistical methods are now commonly used to take into account the expected lack of independence of observations across different species (due to their phylogenetic relatedness) when computing correlations or regressions among traits. The methods are often interpreted as removing that part of the regression or correlation that is an artifact due to phylogeny and there is an expectation that the corrected regression or correlation coefficients will usually be closer to zero. It is shown here that this is not an accurate way to interpret these methods. The effect of taking phylogeny into account is to reduce the variance of the estimated regression or correlation coefficients. Their means are not because since estimates of regression coefficients are unbiased whether or not the correct phylogeny is taken into account. Estimates of correlations are only slightly biased (and in the opposite direction that many expect).

Body size is widely believed to affect the occurrence of sexual maturation. Recent studies have used changes in the age-specific body size at which the probability of maturing is 50%, a feature of probabilistic reaction norms, to quantify purported evolution of life histories. However, body size results from a combination of growth rates during successive developmental stages. Therefore, to understand the evolution of the maturation schedule, it is necessary to comprehend the relationships among body size, growth history, and maturation schedule. We examined the relationships among body size, previous growth history, and maturation probability in chum salmon (Oncorhynchus keta). In this study, previous growth history was estimated from yearly specific growth increments that provide information describing body size. Previous growth history was found to be more closely linked to maturation probability than body size. The most recent growth condition was the most important factor affecting whether a fish matured during the subsequent breeding season. Because individuals of similar body size and same age can have different growth histories, the relationship between body size and maturation probability could be plastically modified by growth history. This may violate an assumption required to infer evolution, namely that size-related maturation trends in probabilistic reaction norms are immune to growth history.