Molecular genetic analysis of phenotypic variation has revealed many examples of evolutionary change in the developmental pathways that control plant and animal morphology. A major challenge is to integrate the information from diverse organisms and traits to understand the general patterns of developmental evolution. This integration can be facilitated by evolutionary metamodels—traits that have undergone multiple independent changes in different species and whose development is controlled by well-studied regulatory pathways. The metamodel approach provides the comparative equivalent of experimental replication, allowing us to test whether the evolution of each developmental pathway follows a consistent pattern, and whether different pathways are predisposed to different modes of evolution by their intrinsic organization. A review of several metamodels suggests that the structure of developmental pathways may bias the genetic basis of phenotypic evolution, and highlights phylogenetic replication as a value-added approach that produces deeper insights into the mechanisms of evolution than single-species analyses.

Inbreeding depression is one of the main forces opposing the evolution of self-fertilization. Of central importance is the hypothesis that inbreeding depression and selfing coevolve antagonistically, generating either low selfing rate and high inbreeding depression or vice versa. However, there is limited evidence for this coevolution within species. We investigated this topic in the hermaphroditic snail Physa acuta. In this species, isolated individuals delay the onset of egg laying compared to individuals having access to mates. Longer delays (“waiting times”) indicate more intense selfing avoidance. We measured inbreeding depression and waiting time in a large quantitative-genetic experiment (281 outbred families derived from 26 natural populations). We observed large genetic variance for both traits and a strong positive genetic covariance between them, most of which resided within rather than among populations. It means that, within populations, individuals with higher mutation load avoided selfing more strongly on average. This genetic covariance may result from pleiotropy and/or linkage disequilibrium. Whatever its genetic architecture, the fact it emerges specifically when individuals are deprived of mates suggests it is not fortuitous and rather reflects the action of natural selection. We conclude that a diversity of mating strategies can arise within populations subjected to variation in inbreeding depression.

The evolution of matrotrophy introduces the potential for genomic conflicts between mothers and embryos. These conflicts are hypothesized to accelerate the evolution of reproductive isolation and to influence the evolution of life-history traits, reproductive structures, and genomic imprinting. These hypotheses assume offspring can influence the amount of maternal investment they receive and that there is a trade-off between maternal investment into individual offspring and maternal survival or fecundity. We used field data and laboratory crosses to test whether these assumptions are met in the matrotrophic poeciliid fish Heterandria formosa. Comparisons of life histories between two natural populations demonstrated a trade-off between the level of maternal investment into individual embryos and maternal fecundity. Laboratory crosses between individuals from these populations revealed that offspring genotype exerts an influence on the level of maternal investment and affects maternal fecundity through higher rates of embryo abortion and lower numbers of full-term offspring. Our results show that the prerequisites for parent-offspring conflict to be a potent evolutionary force in poeciliid fish are present in H. formosa. However, determining whether this conflict has shaped maternal investment in nature will require disentangling any effects of conflict from those of several ecological factors that are themselves correlated with the expected intensity of conflict.

Trade-offs between selected and nonselected environments are often assumed to exist during adaptation. This phenomenon is prevalent in microbial metabolism, where many organisms have come to specialize on a narrow breadth of substrates. One wellstudied example is methylotrophic bacteria that can use single-carbon (C₁) compounds as their sole source of carbon and energy, but generally use few, if any, multi-C compounds. Here, we use adaptation of experimental populations of the model methylotroph, Methylobacterium extorquens AM1, to C₁ (methanol) or multi-C (succinate) compounds to investigate specialization and trade-offs between these two metabolic lifestyles. We found a general trend toward trade-offs during adaptation to succinate, but this was neither universal nor showed a quantitative relationship with the extent of adaptation. After 1500 generations, succinate-evolved strains had a remarkably bimodal distribution of fitness values on methanol: either an improvement comparable to the strains adapted on methanol or the complete loss of the ability to grow on C₁ compounds. In contrast, adaptation to methanol resulted in no such trade-offs. Based on the substantial, asymmetric loss of C₁ growth during growth on succinate, we suggest that the long-term maintenance of C₁ metabolism across the genus Methylobacterium requires relatively frequent use of C₁ compounds to prevent rapid loss.

Adaptive radiation occurs when divergent natural selection in different environments leads to phenotypic differentiation. The pleiotropic effects of underlying genes can either promote or constrain this diversification. Identifying the pleiotropic effects of genes responsible for divergent traits, and testing how the environment influences these effects, can therefore help to provide an understanding of how ecology drives evolutionary change between populations. Positive selection on low-armor alleles at the Ectodysplasin (Eda) locus in threespine stickleback has led to the repeated evolution of reduced armor in populations following freshwater colonization by fully armored marine sticklebacks. Here, we demonstrate that Eda has environmentally determined pleiotropic effects on armor and growth. When raised in freshwater, reduced armor sticklebacks carrying “low” alleles at Eda had increased growth rate relative to fully armored sticklebacks carrying “complete” alleles. In saltwater treatments this growth advantage was present during juvenile growth but lost during adult growth, suggesting that in this environment stickleback are able to develop full armor plates without sacrificing overall growth rate. The environment specific pleiotropic effects of Eda demonstrate that ecological factors can mediate the influence of genetic architecture in driving phenotypic evolution. Furthermore, because size is important for mate choice in stickleback, the growth rate differences influenced by Eda may have effects on reproductive isolation between marine and freshwater populations.

The maintenance of genetic variation in traits closely associated with fitness remains a key unresolved issue in evolutionary genetics. One important qualification on the observation of genetic variation in fitness-related traits is that such traits respond asymmetrically to selection, evolving to a greater extent in the direction of lower fitness. Here we test the hypothesis that standing genetic variation in fitness-related traits is principally maintained for unfit phenotypes. Male Drosophila bunnanda vary in mating success (the primary determinant of male fitness) due to female mate choice. We used competitive mating success to partitioning males into two groups: successful (high fitness) and unsuccessful (low fitness). Relative to successful males, unsuccessful males harbored considerably greater levels of additive genetic variation for sexual signaling traits. This genetic asymmetry was detected for a multivariate trait that we demonstrated was not directly under stabilizing sexual selection, leading us to conclude the trait was under apparent stabilizing selection. Consequently, our results suggest genetic variance might be biased toward low fitness even for traits that are not themselves the direct targets of selection. Simple metrics of genetic variance are unlikely to be adequate descriptors of the complex nature of the genetic basis of traits under selection.

The genetic consequences of climate-driven range fluctuation during the Pleistocene have been well studied for temperate species, but cold-adapted (e.g., alpine, arctic) species that may have responded uniquely to past climatic events have received less attention. In particular, we have no a priori expectation for long-term evolutionary consequences of elevation shifts into and out of sky islands by species adapted to alpine habitats. Here, we examined the influence of elevation shifts on genetic differentiation and historical demography in an alpine specialist, the American pika (Ochotona princeps). Pika populations are divided into five genetic lineages that evolved in association with separate mountain systems, rather than lineages that reflect individual sky islands. This suggests a role for glacial-period elevation shifts in promoting gene flow among high-elevation populations and maintaining regional cohesion of genetic lineages. We detected a signature of recent demographic decline in all lineages, consistent with the expectation that Holocene climate warming has driven range retraction in southern lineages, but unexpected for northern populations that presumably represent postglacial expansion. An ecological niche model of past and future pika distributions highlights the influence of climate on species range and indicates that the distribution of genetic diversity may change dramatically with continued climate warming.

The phenotypic effects of genetic and environmental manipulations have been rarely investigated simultaneously. In addition to phenotypic plasticity, their effect on the amount and directions of genetic and phenotypic variation is of particular evolutionary importance because these constitute the material for natural selection. Here, we used heterozygous insertional mutations of 16 genes involved in the formation of the Drosophila wing. The flies were raised at two developmental temperatures (18°C and 28°C). Landmark-based geometric morphometries was used to analyze the variation of the wing size and shape at different hierarchical levels: among genotypes and temperatures; among individuals within group; and fluctuating asymmetry (FA). Our results show that (1) the phenotypic effects of the mutations depend on temperature; (2) reciprocally, most mutations affect wing plasticity; (3) both temperature and mutations modify the levels of FA and of among individuals variation within lines. Remarkably, the patterns of shape FA seem unaffected by temperature whereas those associated with individual variation are systematically altered. By modifying the direction of available phenotypic variation, temperature might thus directly affect the potential for further evolution. It suggests as well that the developmental processes responsible for developmental stability and environmental canalization might be partially distinct.

The distribution of a phenotypic state is often discontinuous and dispersed. An example of such a distribution can be found in the shell shapes of terrestrial gastropods, which exhibit a bimodal distribution whereby species possess either a tall shell or a flat shell. Here we propose a simple model to test the hypothesis that the bimodal distribution relates to the optimum shape for shell balance on the substrates. This model calculates the theoretical shell balance by moment and obtains empirical distribution of shell shape by compiling published data and performing a new analysis. The solution of the model supports one part of the hypothesis, showing that a low-spired shell is the best balanced and is better suited for locomotion on horizontal surface. Additionally, the model shows that both high- and low-spired shells are well balanced and suited on vertical surfaces. The shell with a spire index (shell height divided by diameter) of 1.4 is the least well balanced as a whole. Thus, spire index is expected to show a bimodal distribution with a valley at 1.4. This expectation was supported by empirical distribution of a spire index, suggesting that the bimodality of shell shape in terrestrial gastropods is related to shell balance.

An intralocus genetic conflict occurs when a locus is selected in opposing directions in different subsets of a population. Populations with two sexes have the potential to host a pair of distinct intralocus conflicts: sexual antagonism and parental antagonism. In this article, we examine the population genetic consequences of these conflicts for X-linked genes. Both conflicts are capable of maintaining genetic variation in a population, but to different degrees. For weak sexual antagonism, the X chromosome has a higher opportunity for polymorphism than the autosomes. For parental antagonism, there is a very limited opportunity for polymorphism on the X chromosome relative to autosomes or to sexual antagonism. X-linkage introduces an asymmetry in the inheritance and expression of sexually and parentally antagonistic genes that leads to a biased fixation of alleles with certain effects. We find little support for the commonly held intuition that the X chromosome should be biased toward fixing female-beneficial alleles. Contrary to this intuition, we find that the X chromosome is biased toward fixation of male-beneficial alleles for much of the range of dominance. Additionally, we find that the X chromosome is more favorable to the fixation of alleles that are beneficial when maternally derived.

Natural selection may favor two very different types of social behaviors that have costs in vital rates (fecundity and/or survival) to the actor: helping behaviors, which increase the vital rates of recipients, and harming behaviors, which reduce the vital rates of recipients. Although social evolutionary theory has mainly dealt with helping behaviors, competition for limited resources creates ecological conditions in which an actor may benefit from expressing behaviors that reduce the vital rates of neighbors. This may occur if the reduction in vital rates decreases the intensity of competition experienced by the actor or that experienced by its offspring. Here, we explore the joint evolution of neutral recognition markers and marker-based costly conditional harming whereby actors express harming, conditional on actor and recipient bearing different conspicuous markers. We do so for two complementary demographic scenarios: finite panmictic and infinite structured populations. We find that marker-based conditional harming can evolve under a large range of recombination rates and group sizes under both finite panmictic and infinite structured populations. A direct comparison with results for the evolution of marker-based conditional helping reveals that, if everything else is equal, marker-based conditional harming is often more likely to evolve than marker-based conditional helping.

Estimating dispersal distances from population genetic data provides an important alternative to logistically taxing methods for directly observing dispersal. Although methods for estimating dispersal rates between a modest number of discrete demes are well developed, methods of inference applicable to “isolation-by-distance” models are much less established. Here, we present a method for estimating ρσ2, the product of population density (ρ) and the variance of the dispersal displacement distribution (σ2). The method is based on the assumption that low-frequency alleles are identical by descent. Hence, the extent of geographic clustering of such alleles, relative to their frequency in the population, provides information about ρσ2. We show that a novel likelihood-based method can infer this composite parameter with a modest bias in a lattice model of isolation-by-distance. For calculating the likelihood, we use an importance sampling approach to average over the unobserved intraallelic genealogies, where the intraallelic genealogies are modeled as a pure birth process. The approach also leads to a likelihood-ratio test of isotropy of dispersal, that is, whether dispersal distances on two axes are different. We test the performance of our methods using simulations of new mutations in a lattice model and illustrate its use with a dataset from Arabidopsis thaliana.

Quantitative genetics models have been used to predict the constraints on local adaptation caused by gene flow between populations under migration—selection balance. One key assumption of this approach is that genetic values within a population are normally distributed. Gene flow, however, may generate distributions that are skewed toward the immigrant's mean value. If the response to selection from a skewed distribution is different from that expected under the assumption of normality, models may result in inaccurate predictions. We use individual-based computer simulations to explore this problem, comparing our results to a recent model developed by Hendry et al. (2001). We show that this model underestimates the equilibrium divergence between populations at migration—selection balance. The extent of this underestimation is correlated with the amount of genetic skew generated by migration and is partly explained by the fact that the analytical model ignores direct selection against hybrid phenotypes. We also show that all else being equal, response to selection in a population with a skewed distribution of genotypes is greater than in a population with normally distributed genotypes. The production of skew under migration—selection balance, however, is itself dependent upon the genetic architecture, with greater deviations from normality produced when alleles contributing to population differentiation have very different effect sizes. We find that both the skew and discrepancies between the models are greatest at intermediate migration rates and moderate to strong selection, which is exactly the region of parameter space that is most empirically relevant.

The spatial structuring of intraspecific genetic diversity is the result of random genetic drift, natural selection, migration, mutation, and their interaction with historical processes. The contribution of each has been typically difficult to estimate, but recent advances in statistical genetics have provided valuable new investigative tools to tackle such complexity. Using a combination of such methods, we examined the roles of environment (i.e., natural selection), random genetic processes (i.e., drift), and demography and life histories (e.g., feeding migrations) on population structure of a widely distributed and abundant marine pelagic fish of economic importance, Atlantic herring (Clupea harengus). Individuals were collected during peak spawning time from 19 spawning locations spanning the region from the western North Sea to the eastern Baltic Sea (N = 1859, eight microsatellite loci). We carried out separate analyses of neutral and selected genetic variation, which allowed us to establish that the two most important factors affecting population structure were selection due to salinity at spawning sites and feeding migrations. The genetic signal left by the demographic history of herring, on the other hand, seems to have been largely eroded, which is not surprising given the large reproductive potential and presumed enormous local effective population sizes of pelagic fish that constrain the effect of stochastic processes. The approach we used can in principle be applied to any abundant and widely distributed aquatic or terrestrial species.

We investigate how different rates of environmental change affect adaptive outcomes and dynamics by selecting Chlamydomonas populations for over 200 generations in environments where the rate of change varies. We find that slower rates of environmental change result in end populations that grow faster and pay a lower cost of adaptation than populations that adapt to a sudden change of the same magnitude. We detected partial selective sweeps in adapting populations by monitoring changes in marker frequency in each population. Although populations adapting to a sudden environmental change showed evidence of mutations of large effect segregating early on, populations adapting to slow rates of change showed patterns that were consistent with mutations of relatively small effect occurring at less predictable times. This work suggests that rates of environmental change may fundamentally alter adaptive dynamics and outcomes of adaptation by changing the size and timing of fitness increases. We suggest that using mutations of smaller effect during adaptation may result in lower levels of pleiotropy and historical constraints, which could in turn result in higher fitness by the end of the experiment.

Although several studies have demonstrated that frequency-dependent effects can promote the maintenance of cooperative behavior in microbes, experimental evidence of frequency-dependent effects in cooperative animal societies is rare. We staged mixed phenotype feeding bouts in the spider Anelosimus studiosus, which shows a within-population social polymorphism, to determine how phenotype frequency affects the foraging success of the social (cooperative) and asocial (cheater) phenotypes. Foraging performance was inferred from average change in percent mass for the respective phenotypes after staged group foraging events. We then performed a field census of multifemale colonies of A. studiosus to determine the phenotypic composition of naturally occurring colonies. Our data indicate that asocial (i.e., cheater) individuals experience negative frequency-dependent foraging success in staged foraging contests. Asocial individuals outperform social individuals when their representation is low, but lose this competitive advantage as their relative numbers increase. Naturally occurring colonies, on average, contained 58.33% social and 41.67% asocial individuals.

In proovigenic parasitoids such as Leptopilina boulardi, the female emerges with a limited egg load and no further eggs are produced during its adult life. A female thus runs the risk of exhausting this limited supply of eggs before the end of her life. Given that the production of an egg is costly, what is the evolutionarily stable egg load at emergence? This question has attracted a lot of attention in the last decade. Here, we analyze a model that allows us to track both the evolution and the population dynamics of a solitary, proovigenic parasitoid. First, we show how host-parasitoid dynamics feedbacks on the evolution of parasitoid egg load. Second, we use this model to consider the situation in which the parasitoid can be infected by a virus that manipulates the oviposition behavior of the females. In particular, we model the effect of the LbFV virus in L. boulardi, a virus that is known to enhance its horizontal transmission by increasing superparasitism (i.e., the laying of eggs in a host already parasitized). Specifically, we model (1) the effect of the virus on parasitoid egg load strategies, and (2) the evolution of egg load manipulation by the virus. This analysis yields two alternative, yet not mutually exclusive, adaptive explanations for the observation that females infected by the virus harbor higher egg loads than uninfected females. Infected females could either respond plastically to the infection status, or be manipulated by the virus. Further experimental work is required to distinguish between these two hypotheses. In a broader context, we present a general theoretical framework that allows us to study the epidemiology, the evolution, the coevolution, and the evolution of manipulation of various reproductive strategies of parasitoids.

Avian plumage colors are frequently used in studies of sexual selection, yet surprisingly little is known about how these traits evolve under different mating systems. We compared historical rates of divergence in male color patterns among the oropendolas and caciques (genera Cacicus, Gymnostinops, Ocyalus, and Psarocolius), a group with both polygynous and monogamous representatives. Reconstructing the evolution of individual color patches on a molecular phylogeny showed that overall color patterns have changed much more rapidly in oropendolas, which comprise two groups that evolved polygyny independently, than in caciques, which are predominantly monogamous. None of these taxa are notably sexually dichromatic, however, suggesting that higher rates of plumage evolution occurred in both sexes rather than just males. Despite high rates of change, color patterns show few examples of convergence among taxa, similar to the lack of homoplasy in male song among oropendolas but in a stark contrast to the repeated convergence in both plumage and song patterns found in a closely related, monogamous clade, the New World orioles (Icterus). Our results support previous suggestions that display traits evolve more rapidly and with less homoplasy in polygynous mating systems, and we provide surprising evidence that these patterns may occur in both sexes.

In recently diverged species, ancestral polymorphism and introgression can cause incongruence between gene and species trees. In the face of hybridization, few genomic regions may exhibit reciprocal monophyly, and these regions, usually evolving rapidly under selection, may be important for the maintenance of species boundaries. In animals with internal fertilization, genes encoding seminal protein are candidate barrier genes. Recently diverged hybridizing species such as the field crickets Gryllus firmus and G. pennsylvanicus, offer excellent opportunities to investigate the origins of barriers to gene exchange. These recently diverged species form a well-characterized hybrid zone, and share ancestral polymorphisms across the genome. We analyzed DNA sequence divergence for seminal protein loci, housekeeping loci, and mtDNA, using a combination of analytical approaches and extensive sampling across both species and the hybrid zone. We report discordant genealogical patterns and differential introgression rates across the genome. The most dramatic outliers, showing near-zero introgression and more structured species trees, are also the only two seminal protein loci under selection. These are candidate barrier genes with possible reproductive functions. We also use genealogical data to examine the demographic history of the field crickets and the current structure of the hybrid zone.

Diverse animals have evolved an ability to collect antimicrobial compounds from the environment as a means of reducing infection risk. Honey bees battle an extensive assemblage of pathogens with both individual and “social” defenses. We determined if the collection of resins, complex plant secretions with diverse antimicrobial properties, acts as a colony-level immune defense by honey bees. Exposure to extracts from two sources of honey bee propolis (a mixture of resins and wax) led to a significantly lowered expression of two honey bee immune-related genes (hymenoptaecin and Am Eater in Brazilian and Minnesota propolis, respectively) and to lowered bacterial loads in the Minnesota (MN) propolis treated colonies. Differences in immune expression were also found across age groups (third-instar larvae, 1-day-old and 7-day-old adults) irrespective of resin treatment. The finding that resins within the nest decrease investment in immune function of 7-day-old bees may have implications for colony health and productivity. This is the first direct evidence that the honey bee nest environment affects immune-gene expression.

The availability of increasingly comprehensive phylogenies has provided unprecedented opportunities to assess macroevolutionary patterns, yet studies on invertebrate diversification are few. In particular, despite the ecological and evolutionary importance of ants, little is known about their tempo and mode of diversification. Recent advances in ant phylogenetics can now provide a basis for rigorous analyses of the diversification of ant lineages. The goals of the present study are threefold. First, we demonstrate that a hypothesized disproportionate increase in ant diversification during the angiosperm radiation is largely artifactual. Rather, current evidence points to a fairly constant rate of lineage growth during its history. Moreover, an analysis of diversification patterns across the ant phylogeny indicates considerable rate heterogeneity among lineages. Indeed, and contrary to the expectation if lineages had experienced a single rate of lineage increase, we found no correspondence between genus age and diversity. Finally, we demonstrate a statistically significant phylogenetic signal in ant diversification: closely related genera have diversities that are more similar to one another than one would expect by chance. This suggests that the capacity for diversification may be itself a biological trait that evolved during the radiation of the family Formicidae.

Disturbance is thought to be a major factor influencing patterns of biodiversity. In addition, disturbance can modify community composition if there are species specific trade-offs between fitness and disturbance tolerance. Here, we examine the role of disturbance on the evolution of coexisting biofilm-forming morphotypes of Pseudomonas fluorescens maintained in spatially structured laboratory microcosms. We identified four heritably stable ecotypes that varied significantly in their competitiveness under different disturbance treatments. Furthermore, we identified significant trade-offs in competitiveness across disturbance treatments for three of four of these ecotypes. These trade-offs modified dominance relationships between strains and thus altered community composition, with a peak of ecotype diversity occurring at intermediate disturbance frequencies.