Testing whether some species traits have a significant effect on diversification rates is central in the assessment of macroevolutionary theories. However, we still lack a powerful method to tackle this objective. I present a new method for the statistical analysis of diversification with species traits. The required data are observations of the traits on recent species, the phylogenetic tree of these species, and reconstructions of ancestral values of the traits. Several traits, either continuous or discrete, and in some cases their interactions, can be analyzed simultaneously. The parameters are estimated by the method of maximum likelihood. The statistical significance of the effects in a model can be tested with likelihood ratio tests. A simulation study showed that past random extinction events do not affect the Type I error rate of the tests, whereas statistical power is decreased, though some power is still kept if the effect of the simulated trait on speciation is strong. The use of the method is illustrated by the analysis of published data on primates. The analysis of these data showed that the apparent overall positive relationship between body mass and species diversity is actually an artifact due to a clade-specific effect. Within each clade the effect of body mass on speciation rate was in fact negative. The present method allows to take both effects (clade and body mass) into account simultaneously.

We present here a stochastic two-locus, two-habitat model for the evolution of migration with local adaptation and kin selection. One locus determines the migration rate while the other causes local adaptation. We show that the opposing forces of kin competition and local adaptation can lead to the existence of one or two convergence stable migration rates, notably depending on the recombination rate between the two loci. We show that linkage between migration and local adaptation loci has two antagonist effects: when linkage is tight, cost of local adaptation increases, leading to smaller equilibrium migration rates. However, when linkage is tighter, the population structure at the migration locus tends to be very high because of the indirect selection, and thus equilibrium migration rates increases. This result, qualitatively different from results obtained with other models of migration evolution, indicates that ignoring drift or the detail of the genetic architecture may lead to incorrect conclusions.

Under the coalescent model for population divergence, lineage sorting can cause considerable variability in gene trees generated from any given species tree. In this paper, we derive a method for computing the distribution of gene tree topologies given a bifurcating species tree for trees with an arbitrary number of taxa in the case that there is one gene sampled per species. Applications for gene tree distributions include determining exact probabilities of topological equivalence between gene trees and species trees and inferring species trees from multiple datasets. In addition, we examine the shapes of gene tree distributions and their sensitivity to changes in branch lengths, species tree shape, and tree size. The method for computing gene tree distributions is implemented in the computer program COAL.

The initial evolution of aposematic and mimetic antipredator signals is thought to be paradoxical because such coloration is expected to increase the risk of predation before reaching a stage when predators associate it effectively with a defense. We propose, however, that constraints associated with the alternative strategy, cryptic coloration, may facilitate the evolution of antipredator signals and thus provide a solution for the apparent paradox. We tested this hypothesis first using an evolutionary simulation to study the effect of a constraint due to habitat heterogeneity, and second using a phylogenetic comparison of the Lepidoptera to investigate the effect of a constraint due to prey motility. In the evolutionary simulation, antipredator warning coloration had an increased probability to invade the prey population when the evolution of camouflage was constrained by visual difference between microhabitats. The comparative study was done between day-active lepidopteran taxa, in which camouflage is constrained by motility, and night-active taxa, which rest during the day and are thus able to rely on camouflage. We compared each of seven phylogenetically independent day-active groups with a closely related nocturnal group and found that antipredator signals have evolved at least once in all the diurnal groups but in none of their nocturnal matches. Both studies lend support to our idea that constraints on crypsis may favor the evolution of antipredator warning signals.

Gametophytic self-incompatibility (SI) in plants is a widespread mechanism preventing self-fertilization and the ensuing inbreeding depression, but it often evolves to self-compatibility. We analyze genetic mechanisms for the breakdown of gametophytic SI, incorporating a dynamic model for the evolution of inbreeding depression allowing for partial purging of nearly recessive lethal mutations by selfing, and accounting for pollen limitation and sheltered load linked to the S-locus. We consider two mechanisms for the breakdown of gametophytic SI: a nonfunctional S-allele and an unlinked modifier locus that inactivates the S-locus. We show that, under a wide range of conditions, self-compatible alleles can invade a self-incompatible population. Conditions for invasion are always less stringent for a nonfunctional S-allele than for a modifier locus. The spread of self-compatible genotypes is favored by extremely high or low selfing rates, a small number of S-alleles, and pollen limitation. Observed parameter values suggest that the maintenance of gametophytic SI is caused by a combination of high inbreeding depression in self-incompatible populations coupled with intermediate selfing rates of the self-compatible genotypes and sheltered load linked to the S-locus.

To test the hypothesis that genes have extended phenotypes on the community, we quantified how genetic differences among cottonwoods affect the diversity, abundance, and composition of the dependent arthropod community. Over two years, five major patterns were observed in both field and common-garden studies that focused on two species of cottonwoods and their naturally occurring F₁ and backcross hybrids (collectively referred to as four different cross types). We did not find overall significant differences in arthropod species richness or abundance among cottonwood cross types. We found significant differences in arthropod community composition among all cross types except backcross and narrowleaf cottonwoods. Thus, even though we found similar richness among cross types, the species that composed the community were significantly different. Using vector analysis, we found that the shift in arthropod community composition was correlated with percent Fremont alleles in the host plant, which suggests that the arthropod community responds to the underlying genetic differences among trees. We found 13 arthropod species representing different trophic levels that were significant indicators of the four different cross types. Even though arthropod communities changed in species composition from one year to the next, the overall patterns of community differences remained remarkably stable, suggesting that the genetic differences among cross types exert a strong organizing influence on the arthropod community. Together, these results support the extended phenotype concept. Few studies have observationally and experimentally shown that entire arthropod communities can be structured by genetic differences in their host plants. These findings contribute to the developing field of community genetics and suggest a strategy for conserving arthropod diversity by promoting genetic diversity in their host plants.

Parasite-mediated selection is potentially of great importance in modulating genetic diversity. Genetic variation for resistance, the fuel for natural selection, appears to be common in host-parasite interactions, but responses to selection are rarely observed. In the present study, we tested whether environmental variation could mediate infection and determine evolutionary outcomes. Temperature was shown to dramatically alter the potential for parasite-mediated selection in two independent laboratory infection experiments at four temperatures. The bacterial parasite, Pasteuria ramosa, was extremely virulent at 20°C and 25°C, sterilizing its host, Daphnia magna, so that females often never produced a single brood. However, at 10°C and 15°C, the host-parasite interaction was much more benign, as nearly all females produced broods before becoming sterile. This association between virulence and temperature alone could stabilize coexistence and lead to the maintenance of diversity, because it would weaken parasite-mediated selection during parts of the season. Additionally, highly significant genotype-by-environment interactions were found, with changes in clone rank order for infection rates at different temperatures. Our results clearly show that the outcome of parasite-mediated selection in this system is strongly context dependent.

Phenotypic plasticity is an important response mechanism of plants to environmental heterogeneity. Here, we explored the genetic basis of plastic responses of Arabidopsis thaliana to water deficit by experimentally mapping quantitative trait loci (QTL) in two recombinant inbred populations (Cvi × Ler and Ler × Col). We detected genetic variation and significant genotype-by-environment interactions for many traits related to water use. We also mapped 26 QTL, including six for carbon isotope composition (δ¹³C). Negative genetic correlations between fruit length and fruit production as well as between flowering time and branch production were corroborated by QTL colocalization, suggesting these correlations are due to pleiotropy or physical linkage. Water-limited plants were more apically dominant with greater root:shoot ratios and higher δ¹³C (higher water-use efficiency) when compared to well-watered plants. Many of the QTL effects for these traits interacted significantly with the irrigation treatment, suggesting that the observed phenotypic plasticity is genetically based. We specifically searched for epistatic (QTL-QTL) interactions using a two-dimensional genome scan, which allowed us to detect epistasis regardless of additive genetic effects. We found several significant QTL-QTL interactions including three that exhibited environmental dependence. These results provide preliminary evidence for proposed genetic mechanisms underlying phenotypic plasticity.

Conspecific sperm precedence occurs when females are exposed to sperm from males of multiple species, but preferentially use sperm of a conspecific. Conspecific sperm precedence and its mechanisms have been documented widely in terrestrial species, in which complex female behaviors or reproductive tract morphologies can allow many opportunities for female choice and sperm competition, however, the opportunity for conspecific sperm precedence in free spawning marine invertebrates has been largely ignored. Two sea urchin species, Echinometra oblonga and E. sp. C, have high levels of interspecific fertilization in no-choice lab crosses, but no natural hybrids have been found. We performed competitive fertilization assays to test for conspecific sperm precedence and found that eggs of both species showed a marked preference for conspecific sperm when fertilized with heterospecific sperm mixtures. Strong rejection of heterospecific sperm would not have been predicted from no-choice assays and helps explain the lack of natural hybrids. We also found significant variation in hybridization success among crosses. Conspecific sperm precedence in free spawning invertebrates shows that the simple surfaces of eggs and sperm provide ample opportunity for egg choice and sperm competition even in the absence of intricate behavior or complex reproductive morphologies.

Numerous studies have reported that females benefit from mating with multiple males (polyandry) by minimizing the probability of fertilization by genetically incompatible sperm. Few, however, have directly attributed variation in female reproductive success to the fertilizing capacity of sperm. In this study we report on two experiments that investigated the benefits of polyandry and the interacting effects of males and females at fertilization in the free-spawning Australian sea urchin Heliocidaris erythrogramma. In the first experiment we used a paired (split clutch) experimental design and compared fertilization rates within female egg clutches under polyandry (eggs exposed to the sperm from two males simultaneously) and monandry (eggs from the same female exposed to sperm from each of the same two males separately). Our analysis revealed a significant fertilization benefit of polyandry and strong interacting effects of males and females at fertilization. Further analysis of these data strongly suggested that the higher rates of fertilization in the polyandry treatment were due to an overrepresentation of fertilizations due to the most compatible male. To further explore the interacting effects of males and females at fertilization we performed a second factorial experiment in which four males were crossed with two females (in all eight combinations). In addition to confirming that fertilization success is influenced by male × female interactions, this latter experiment revealed that both sexes contributed significant variance to the observed patterns of fertilization. Taken together, these findings highlight the importance of male × female interactions at fertilization and suggest that polyandry will enable females to reduce the cost of fertilization by incompatible gametes.

The perceived wide geographic range of organisms in the sea, facilitated by ready dispersal of waterborne dispersal stages, is a challenge for hypotheses of marine speciation but a boon to efforts of marine conservation. Wide species ranges are especially striking in the reef-rich Indo-west Pacific, the largest and most diverse marine biogeographic region, extending across half the planet. The insular marine biota of the tropical Pacific is characterized by wide-ranging species and provides the most striking examples of long distance dispersal, with endemism largely confined to the most remote island groups. Here we show that the gastropod Astralium “rhodostomum” has developed endemic clades on almost every Pacific archipelago sampled, a pattern unprecedented in marine biogeography, and reminiscent of the terrestrial biota of oceanic islands. Mitochondrial DNA sequences indicate that this species-complex is comprised of at least 30 geographically isolated clades, separated by as little as 180 km. Evidence suggests that such fine scale endemism and high diversity is not exceptional, but likely characterizes a substantial fraction of the reef biota. These results imply that (1) marine speciation can regularly occur over much finer spatial scales than generally accepted, (2) the diversity of coral reefs is even higher than suggested by morphology-based estimates, and (3) conservation efforts need to focus at the archipelagic level in the sea as on land.

The extent to which natural selection shapes phenotypic variation has long been a matter of debate among those studying organic evolution. We studied the patterns of gene expression polymorphism and divergence in several datasets that ranged from comparisons between two very closely related laboratory strains of mice to comparisons across a considerably longer time scale, such as between humans and chimpanzees, two species of mice, and two species of Drosophila. The results were analyzed and interpreted in view of neutral models of phenotypic evolution. Our analyses used a number of metrics to show that most mRNA levels are evolutionary stable, changing little across the range of taxonomic distances compared. This implies that, overall, widespread stabilizing selection on transcription levels has prevented greater evolutionary changes in mRNA levels. Nevertheless, the range of rates of divergence is large with highly significant differences in the rate and patterns of transcription divergence across functional classes defined on the basis of the gene ontology annotation (primates and mice datasets) or on the basis of the pattern of sex-biased gene expression (Drosophila). Moreover, rates of divergence of sex-biased genes in the contrast between Drosophila species show a distinct pattern from that observed in the contrast between populations of D. melanogaster. Hence, we discuss the time scale of the changes observed and its consequences for the relationship between variation in gene expression within and between species. Finally, we argue that differences in mRNA levels of the magnitudes observed herein could be explained by a remarkably small number of generations of directional selection.

The hypothesis that sexual selection drives the evolution of condition dependence is not firmly supported by empirical evidence, and the process remains poorly understood. First, even though sexual competition typically involves multiple traits, studies usually compare a single sexual trait with a single “control” trait, ignoring variation among sexual traits and raising the possibility of sampling bias. Second, few studies have addressed the genetic basis of condition dependence. Third, even though condition dependence is thought to result from a form of sex-specific epistasis, the evolution of condition dependence has never been considered in relation to intralocus sexual conflict. We argue that condition dependence may weaken intersexual genetic correlations and facilitate the evolution of sexual dimorphism. To address these questions, we manipulated an environmental factor affecting condition (larval diet) and examined its effects on four sexual and four nonsexual traits in Prochyliza xanthostoma adults. As predicted by theory, the strength of condition dependence increased with degree of exaggeration among male traits. Body shape was more condition dependent in males than in females and, perhaps as a result, genetic and environmental effects on body shape were congruent in males, but not in females. However, of the four male sexual traits, only head length was significantly larger in high-condition males after controlling for body size. Strong condition dependence was associated with reduced intersexual genetic correlation. However, homologous male and female traits exhibited correlated responses to condition, suggesting an intersexual genetic correlation for condition dependence itself. Our findings support the role of sexual selection in the evolution of condition dependence, but reveal considerable variation in condition dependence among sexual traits. It is not clear whether the evolution of condition dependence has mitigated or exacerbated intralocus sexual conflict in this species.

Although female multiple mating is ubiquitous in insects, its adaptive significance remains poorly understood. Benefits to multiple mating can accrue via direct material benefits, indirect genetic benefits, or both. We investigated the effects of multiple mating in the decorated cricket, Gryllodes sigillatus, by simultaneously varying the number of times that females mated and the number of different males with which they mated, measuring aspects of female fecundity and elements of offspring performance and viability. Multiple matings resulted in enhanced female fitness relative to single matings when females mated with different partners, but not when females mated repeatedly with the same male. Specifically, polyandrous females produced significantly more offspring surviving to reproductive maturity than did monogamous females mating once or mating repeatedly with the same male. These results suggest that the benefit females gain from multiple mating is influenced primarily by genetic and not material benefits.

The social spiders are unusual among cooperatively breeding animals in being highly inbred. In contrast, most other social organisms are outbred owing to inbreeding avoidance mechanisms. The social spiders appear to originate from solitary subsocial ancestors, implying a transition from outbreeding to inbreeding mating systems. Such a transition may be constrained by inbreeding avoidance tactics or fitness loss due to inbreeding depression. We examined whether the mating system of a subsocial spider, in a genus with three social congeners, is likely to facilitate or hinder the transition to inbreeding social systems. Populations of subsocial Stegodyphus lineatus are substructured and spiders occur in patches, which may consist of kin groups. We investigated whether male mating dispersal prevents matings within kin groups in natural populations. Approximately half of the marked males that were recovered made short moves (< 5m) and mated within their natal patch. This potential for inbreeding was counterbalanced by a relatively high proportion of immigrant males. In mating experiments, we tested whether inbreeding actually results in lower offspring fitness. Two levels of inbreeding were tested: full sibling versus non-sib matings and matings of individuals within and between naturally occurring patches of spiders. Neither full siblings nor patch mates were discriminated against as mates. Sibling matings had no effect on direct fitness traits such as fecundity, hatching success, time to hatching and survival of the offspring, but negatively affected offspring growth rates and adult body size of both males and females. Neither direct nor indirect fitness measures differed significantly between within patch and between-patch pairs. We tested the relatedness between patch mates and nonpatch mates using DNA fingerprinting (TE-AFLP). Kinship explained 30% of the genetic variation among patches, confirming that patches are often composed of kin. Overall, we found limited male dispersal, lack of kin discrimination, and tolerance to low levels of inbreeding. These results suggest a history of inbreeding which may reduce the frequency of deleterious recessive alleles in the population and promote the evolution of inbreeding tolerance. It is likely that the lack of inbreeding avoidance in subsocial predecessors has facilitated the transition to regular inbreeding social systems.

Trinidad guppies (Poecilia reticulata) are distributed along an environmental gradient in carotenoid availability that limits the carotenoid content of the orange spots of males. The amount of synthetic red pteridines (drosopterins) in the orange spots covaries with the carotenoid content, such that the ratio of the two types of pigments is roughly conserved across streams. Carotenoids and drosopterins have different spectral properties and thus the ratio of the two types of pigments affects the shape of the orange spot reflectance spectrum. Geographic conservation of the carotenoid:drosopterin ratio suggests that males may be under selection to maintain a particular hue. We tested this hypothesis by comparing the pigmentation and coloration of guppies from six streams in the field to that of second-generation descendants of the same populations raised on three dietary carotenoid levels in the laboratory. The results show clearly that the geographic variation in drosopterin production is largely genetic and that the hue of the orange spots is conserved among populations in the field, relative to the laboratory diet groups. This is a countergradient pattern because genetic differences between populations in drosopterin production mask the effect of carotenoid availability on the hue of the orange spots. The potential for countergradient sexual selection to contribute to reproductive isolation between populations is discussed.

Heterozygosity-fitness correlations (HFCs) have been reported in populations of many species. We provide evidence for a positive correlation between genetic variability and growth rate at 12 allozyme loci in a catadromous marine fish species, the European eel (Anguilla anguilla L.). More heterozygous individuals show a significantly higher length and weight increase and an above average condition index in comparison with more homozygous individuals. To a lesser extent, six microsatellite loci show a similar pattern, with positive but not significant correlations between heterozygosity and growth rate. The HFCs observed could be explained by an effect of either direct allozyme overdominance or associative overdominance. Selection affecting some of the allozyme loci would explain the greater strength of the HFCs found at allozymes in comparison with microsatellites and the lack of correlation between MLH at allozymes and MLH at microsatellites. Associative overdominance (where allozyme loci are merely acting as neutral markers of closely linked fitness loci) might provide an explanation for the HFCs if we consider that allozyme loci have a higher chance than microsatellites to be in linkage disequilibrium with fitness loci.

Given that evolutionary divergence in mating signals leads to reproductive isolation in numerous animal taxa, understanding what drives signal divergence is fundamental to our understanding of speciation. Mating signals are thought to diverge via several processes, including (1) as a by-product of morphological adaptation, (2) through direct adaptation to the signaling environment, or (3) to facilitate species recognition. According to the first two hypotheses, birdsongs diversify in different foraging niches and habitats as a product of selection for optimal morphology and efficient sound transmission, respectively. According to the third hypothesis, they diversify as a result of selection against maladaptive hybridization. In this study I test all three hypotheses by examining the influence of morphology, acoustic environment, and the presence of closely related congeners on song structure in 163 species of antbird (Thamnophilidae). Unlike oscine passerines, these Neotropical suboscines make ideal subjects because they develop their songs without learning. In other words, patterns of vocal divergence are not complicated by cultural evolution. In support of the morphological adaptation hypothesis, body mass correlates with the acoustic frequency of songs, and bill size with temporal patterning. These relationships were robust, even when controlling for phylogenetic inertia using independent contrasts, suggesting that there has been correlated evolution between morphological and acoustic traits. The results also support the acoustic adaptation hypothesis: birds which habitually sing in the understory and canopy produce higher-pitched songs than those that sing in the midstory, suggesting that song structure is related to the sound transmission properties of different habitat strata. Finally, the songs of sympatric pairs of closely related species are more divergent than those of allopatric pairs, as predicted by the species recognition hypothesis. To my knowledge, these data provide the first direct evidence that species recognition and ecological adaptation operate in tandem, and that the interplay between these factors drives the evolution of mating signals in suboscine birds.

How often will natural selection drive parallel evolution at the DNA sequence level? More precisely, what is the probability that selection will cause two populations that live in identical environments to substitute the same beneficial mutation? Here I show that, under fairly general conditions, the answer is simple: if a wild-type sequence can mutate to n different beneficial mutations, replicate populations will on average fix the same mutation with probability P = 2/(n 1). This probability, which is derived using extreme value theory, is independent of most biological details, including the length of the gene in question and the precise distribution of fitness effects among alleles. I conclude that the probability of parallel evolution under natural selection is nearly twice as large as that under neutrality.

Sex allocation theory predicts that parents should bias their reproductive investments toward the offspring sex generating the greatest fitness return. When females are the heterogametic sex (e.g., ZW in butterflies, some lizards, and birds), production of daughters is associated with an increased risk of offspring inviability due to the expression of paternal, detrimental recessives on the Z chromosome. Thus, daughters should primarily be produced when mating with partners of high genetic quality. When female sand lizards (Lacerta agilis) mate with genetically superior males, exhibiting high MHC Class I polymorphism, offspring sex ratios are biased towards daughters, possibly due to recruitment of more Z-carrying oocytes when females have assessed the genetic quality of their partners. If our study has general applicability across taxa, it predicts taxon-specific sex allocation effects depending on which sex is the heterogametic one.

It is a well-known phenomenon that islands can support populations of gigantic or dwarf forms of mainland conspecifics, but the variety of explanatory hypotheses for this phenomenon have been difficult to disentangle. The highly venomous Australian tiger snakes (genus Notechis) represent a well-known and extreme example of insular body size variation. They are of special interest because there are multiple populations of dwarfs and giants and the age of the islands and thus the age of the tiger snake populations are known from detailed sea level studies. Most are 5000–7000 years old and all are less than 10,000 years old. Here we discriminate between two competing hypotheses with a molecular phylogeography dataset comprising approximately 4800 bp of mtDNA and demonstrate that populations of island dwarfs and giants have evolved five times independently. In each case the closest relatives of the giant or dwarf populations are mainland tiger snakes, and in four of the five cases, the closest relatives are also the most geographically proximate mainland tiger snakes. Moreover, these body size shifts have evolved extremely rapidly and this is reflected in the genetic divergence between island body size variants and mainland snakes. Within south eastern Australia, where populations of island giants, populations of island dwarfs, and mainland tiger snakes all occur, the maximum genetic divergence is only 0.38%. Dwarf tiger snakes are restricted to prey items that are much smaller than the prey items of mainland tiger snakes and giant tiger snakes are restricted to seasonally available prey items that are up three times larger than the prey items of mainland tiger snakes. We support the hypotheses that these body size shifts are due to strong selection imposed by the size of available prey items, rather than shared evolutionary history, and our results are consistent with the notion that adaptive plasticity also has played an important role in body size shifts. We suggest that plasticity displayed early on in the occupation of these new islands provided the flexibility necessary as the island's available prey items became more depauperate, but once the size range of available prey items was reduced, strong natural selection followed by genetic assimilation worked to optimize snake body size. The rate of body size divergence in haldanes is similar for dwarfs (h_g = 0.0010) and giants (h_g = 0.0020– 0.0025) and is in line with other studies of rapid evolution. Our data provide strong evidence for rapid and repeated morphological divergence in the wild due to similar selective pressures acting in different directions.

In a recent paper, Wright et al. (2003) argue for the hypothesis that greater biologically available energy elevates the rate of molecular evolution. However, their results are also consistent with alternative hypotheses that invoke either environmentally driven variation in effective population sizes, or natural selection, or both. The available energy gradient cited by Wright et al. is linearly correlated with temperature fluctuations, and the observed rate heterogeneity could be a consequence of this environmental variability. The distribution of phylogenetic branch lengths alone is insufficient to distinguish between the hypotheses, and complementary approaches are suggested.

A recent study of environmental effects on rates of molecular evolution in the plant subgenus Mearnsia shows that species occurring in more equatorial latitudes have higher rates of substitution in rDNA sequences as compared to their more southerly congeners (Wright et al. 2003). However, we believe that the statistical approach employed by Wright et al. (2003) insufficiently accounts for the phylogenetic nonindependence of the species examined, given that all six equatorial species of Mearnsia form a clade. To distinguish between the effect of latitude and that of phylogenetic nonindependence, we have employed a variety of comparative approaches that use independent contrasts to test for an effect of environment across this entire subgenus. We find very little evidence for an effect of latitude on rate of molecular evolution using these approaches and believe that the shared evolutionary history of the clade is a plausible explanation of the apparent rate difference between equatorial and subequatorial Mearnsia species.