As a critical framework for addressing a diversity of evolutionary and ecological questions, any method that provides accurate and detailed phylogeographic inference would be embraced. What is difficult to understand is the continued use of a method that not only fails, but also has never been shown to work—nested clade analysis is applied widely even though the conditions under which the method will provide reliable results have not yet been demonstrated. This contradiction between performance and popularity is even more perplexing given the recent methodological and computational advances for making historical inferences, which include estimating population genetic parameters and testing different biogeographic scenarios. Here I briefly review the history of criticisms and rebuttals that focus specifically on the high rate of incorrect phylogeographic inference of nested-clade analysis, with the goal of understanding what drives its unfettered popularity. In this case, the appeal of what nested-clade analysis claims to do—not what the method actually achieves—appears to explain its paradoxical status as a favorite method that fails. What a method promises, as opposed to how it performs, must be considered separately when evaluating whether the method represents a valuable tool for historical inference.

Algorithms are presented to simultaneously estimate probabilities of speciation and extinction, rates of anagenetic and cladogenetic phenotypic evolution, as well as ancestral character states, from a complete ultrametric species-level phylogeny with dates assigned to all bifurcations and one or more phenotypes in three or more extant species, using Metropolis–Hastings Markov Chain Monte Carlo sampling. The algorithms also estimate missing phenotypes of extant species and numbers of speciation events that occurred on all branches of the phylogeny. The algorithms are discussed and their performance is evaluated using simulated data. That evaluation shows that precise estimation of rates of evolution of one or a few phenotypes requires large phylogenies. Estimation accuracy improves with the number of species on the phylogeny.

“Dollo's law” states that, following loss, a complex trait cannot reevolve in an identical manner. Although the law has previously fallen into disrepute, it has only recently been challenged with statistical phylogenetic methods. We employ simulation studies of an irreversible binary character to show that rejections of Dollo's law based on likelihood-ratio tests of transition rate constraints or on reconstructions of ancestral states are frequently incorrect. We identify two major causes of errors: incorrect assignment of root state frequencies, and neglect of the effect of the character state on rates of speciation and extinction. Our findings do not necessarily overturn the conclusions of phylogenetic studies claiming reversals, but we demonstrate devastating flaws in the methods that are the foundation of all such studies. Furthermore, we show that false rejections of Dollo's law can be reduced by the use of appropriate existing models and model selection procedures. More powerful tests of irreversibility require data beyond phylogenies and character states of extant taxa, and we highlight empirical work that incorporates additional information.

Aposematism is one of the great mysteries of evolutionary biology. The evolution of aposematic coloration is poorly understood, but even less understood is the evolution of polymorphism in aposematic signals. Here, we use a phylogeographic approach to investigate the evolution of color polymorphism in Dendrobates pumilio, a well-known poison-dart frog (family Dendrobatidae), which displays perhaps the most striking color variation of any aposematic species. With over a dozen color morphs, ranging from bright red to dull green, D. pumilio provides an ideal opportunity to examine the evolution of color polymorphism and evolutionary shifts to cryptic coloration in an otherwise aposematic species. We constructed a phylogenetic tree for all D. pumilio color morphs from 3051bp of mtDNA sequence data, reconstructed ancestral states using parsimony and Bayesian methods, and tested the recovered tree against constraint trees using parametric bootstrapping to determine the number of changes to each color type. We find strong evidence for nearly maximal numbers of changes in all color traits, including five independent shifts to dull dorsal coloration. Our results indicate that shifts in coloration in aposematic species may occur more regularly than predicted and that convergence in coloration may indicate that similar forces are repeatedly driving these shifts.

Many tropical forest tree species have broad geographic ranges, and fossil records indicate that population disjunctions in some species were established millions of years ago. Here we relate biogeographic history to patterns of population differentiation, mutational and demographic processes in the widespread rainforest tree Symphonia globulifera using ribosomal (ITS) and chloroplast DNA sequences and nuclear microsatellite (nSSR) loci. Fossil records document sweepstakes dispersal origins of Neotropical S. globulifera populations from Africa during the Miocene. Despite historical long-distance gene flow, nSSR differentiation across 13 populations from Costa Rica, Panama, Ecuador (east and west of Andes) and French Guiana was pronounced (F_ST = 0.14, R_ST = 0.39, P < 0.001) and allele-size mutations contributed significantly (R_ST > F_ST) to the divergences between cis- and trans-Andean populations. Both DNA sequence and nSSR data reflect contrasting demographic histories in lower Mesoamerica and Amazonia. Amazon populations show weak phylogeographic structure and deviation from drift–mutation equilibrium indicating recent population expansion. In Mesoamerica, genetic drift was strong and contributed to marked differentiation among populations. The genetic structure of S. globulifera contains fingerprints of drift-dispersal processes and phylogeographic footprints of geological uplifts and sweepstakes dispersal.

Patterns of evolution are believed to vary latitudinally, but our understanding of this variation remains limited. Here we examine how patterns of subspecific diversification vary within species of birds, specifically addressing three questions: (1) Are subspecies more numerous at lower latitudes within species, consistent with greater phenotypic differentiation at lower latitudes? (2) If there are more subspecies at lower latitudes within species, can area of breeding range explain this relationship? and (3) how do latitudinal differences in subspecies within species vary geographically across the globe? Using all species with five or more subspecies from 12 of the most diverse families of birds in the world, we found consistently more subspecies at lower latitudes across all families, both hemispheres, and all continents examined. Despite the positive influence of area on the number of subspecies within species, area did not explain the greater number of subspecies at lower latitudes within species. Global patterns of subspecies support the idea that phenotypic differentiation of populations is greater at lower latitudes within species. If subspecies density provides an index of rates of incipient speciation, then our results support evolutionary hypotheses for the latitudinal diversity gradient that invoke higher tropical speciation rates.

The study of zones of secondary contact provides insight into the maintenance of reproductive isolation. Tension zone theory supplies powerful tools for assessing how dispersal and selection shape hybrid zones. We present a multimodal analysis of phenotypic clines in conjunction with clines at molecular markers in a hybrid zone between Larus glaucescens and Larus occidentalis. We developed a new method to analyze simultaneously clines of quantitative traits and molecular data. Low linkage disequilibrium and the lack of coincidence between clines at six microsatellites, a mitochondrial DNA region, and two phenotypic traits indicated introgression. However, the hypothesis of neutral diffusion was rejected based on evidence that all of the clines were concordant and narrower than expected for neutral clines, indicating some indirect selection. The analysis of phenotypic variance gave evidence of restricted phenotypic introgression and together with the bimodal distribution of phenotypes suggested that disruptive selection is acting across the hybrid zone, especially on the coloration of bare parts. Multimodal analysis of phenotypic clines also highlighted a shift between the peak of intermediates and the cline center, left behind by hybrid zone motion. High-resolution analysis of phenotypes distribution thus proved useful for detecting hybrid zone movement even without temporal data.

Understanding incipient sexual isolation and speciation is an important pursuit in evolutionary biology. The fruit fly Drosophila melanogaster is a useful model to address questions about the early stages of sexual isolation occurring within widespread species. This species exhibits sexual isolation between cosmopolitan and African flies, especially from Zimbabwe populations. In addition, we have recently described another example of partial sexual isolation between some US and Caribbean populations. This and other phenotypic data suggest that Caribbean flies might be segregating African traits. In the present work we study the geographical variation at the pheromone locus desaturase-2, as well as morphology and courtship behavior across the US–Caribbean region. We find that US and Caribbean populations show sharp geographical clines in all traits and demonstrate that Caribbean traits are more similar to those of Africa than to US populations. Further, African traits in the Caribbean are associated with sexual isolation and best explain variation in sexual isolation when all traits are considered together. These results imply that Caribbean mating preferences are likely to be based on African traits and that even at such early stages of sexual isolation, individuals may already cue in on several traits simultaneously during mate choice.

Evolutionary theory predicts that interactions between species such as resource competition or reproductive interference will generate selection for character displacement where similar species co-occur. However, the rate and direction of character displacement will depend not only on the strength of selection for trait divergence, but also on the amount of genetic variation for selected traits and the nature of genetic correlations between them. To assess the importance of genetic constraints for the evolution of character displacement, we examined the genetic architecture of a suite of floral traits previously shown to be under selection in the annual plant Ipomoea hederacea when this species co-occurs with Ipomoea purpurea. We found that the six floral traits we measured are all positively genetically correlated. We also demonstrate, using new statistical approaches, that the predicted response to selection for four of these six traits is substantially constrained by their genetic correlation structure. Most notably, the response to selection for reduced separation of the tallest and shortest anthers, which reduces the degree of detrimental heterospecific pollen flow, is substantially constrained. Our results suggest that the rate of evolution of reproductive character displacement in I. hederacea is limited by the genetic architecture of floral traits.

Evolutionary biologists explain the maintenance of intermediate levels of defense in plant populations as being due to trade-offs, or negative genetic covariances among ecologically important traits. Attempts at detecting trade-offs as constraints on the evolution of defense have not always been successful, leading some to conclude that such trade-offs rarely explain current levels of defense in the population. Using the agricultural pest Ipomoea purpurea, we measured correlations between traits involved in defense to glyphosate, the active ingredient in Roundup, a widely used herbicide. We found significant allocation costs of tolerance, as well as trade-offs between resistance and two measures of tolerance to glyphosate. Selection on resistance and tolerance exhibited differing patterns: tolerance to leaf damage was under negative directional selection, whereas resistance was under positive directional selection. The joint pattern of selection on resistance and tolerance to leaf damage indicated the presence of alternate peaks in the fitness landscape such that a combination of either high tolerance and low resistance, or high resistance and low tolerance was favored. The widespread use of this herbicide suggests that it is likely an important selective agent on weed populations. Understanding the evolutionary dynamics of herbicide defense traits is thus of increasing importance in the context of human-mediated evolution.

Subtle left–right biases are often observed in organisms with an overall bilateral symmetry. The evolutionary significance of these directional asymmetries remains uncertain, however, and scenarios of both developmental constraints and adaptation have been suggested. Reviewing the literature on asymmetry in insect wings, we analyze patterns of directional asymmetry in wing size to evaluate the possible adaptive significance of this character. We found that directional asymmetry in wing size is widespread among insects, with left- and right-biased asymmetries commonly observed. The direction of the asymmetry does not appear to be evolutionarily conserved above the species level. Overall, we argue that the very small magnitude of directional asymmetry, 0.7% of the wing size on average, associated with an extremely imprecise expression, precludes directional asymmetry from playing any major adaptive role.

Environmental niche models, which are generated by combining species occurrence data with environmental GIS data layers, are increasingly used to answer fundamental questions about niche evolution, speciation, and the accumulation of ecological diversity within clades. The question of whether environmental niches are conserved over evolutionary time scales has attracted considerable attention, but often produced conflicting conclusions. This conflict, however, may result from differences in how niche similarity is measured and the specific null hypothesis being tested. We develop new methods for quantifying niche overlap that rely on a traditional ecological measure and a metric from mathematical statistics. We reexamine a classic study of niche conservatism between sister species in several groups of Mexican animals, and, for the first time, address alternative definitions of “niche conservatism” within a single framework using consistent methods. As expected, we find that environmental niches of sister species are more similar than expected under three distinct null hypotheses, but that they are rarely identical. We demonstrate how our measures can be used in phylogenetic comparative analyses by reexamining niche divergence in an adaptive radiation of Cuban anoles. Our results show that environmental niche overlap is closely tied to geographic overlap, but not to phylogenetic distances, suggesting that niche conservatism has not constrained local communities in this group to consist of closely related species. We suggest various randomization tests that may prove useful in other areas of ecology and evolutionary biology.

Enzyme isoforms are found in many cellular reactions, and can differ in the kind of reaction they catalyze, in their substrate affinity, or in their reaction rates. The evolutionary significance of enzyme isoforms is only partially understood. We used mathematical modeling to investigate the hypothesis that isoforms may be favored by selection because they can increase the phenotypic robustness of the system. We modify a model for circadian clock gene expression in Drosophila to incorporate the presence of isoforms in the phosphorylation pathway of the period gene. We consider the case in which different isoforms catalyze the same reaction but have different affinities for the substrate. Stability is increased if there is dynamic control of the expression of isoforms relative to each other. Thus, we show that controlling isoform proportion can be a powerful mechanism for reducing the effects of variations in the values of system parameters, increasing system robustness.

We reassess the coevolution between actinomycete bacteria and fungus-gardening (attine) ants. Actinomycete bacteria are of special interest because they are metabolic mutualists of diverse organisms (e.g., in nitrogen-fixation or antibiotic production) and because Pseudonocardia actinomycetes are thought to serve disease-suppressing functions in attine gardens. Phylogenetic information from culture-dependent and culture-independent microbial surveys reveals (1) close affinities between free-living and ant-associated Pseudonocardia, and (2) essentially no topological correspondence between ant and Pseudonocardia phylogenies, indicating frequent bacterial acquisition from environmental sources. Identity of ant-associated Pseudonocardia and isolates from soil and plants implicates these environments as sources from which attine ants acquire Pseudonocardia. Close relatives of Atta leafcutter ants have abundant Pseudonocardia, but Pseudonocardia in Atta is rare and appears at the level of environmental contamination. In contrast, actinomycete bacteria in the genera Mycobacterium and Microbacterium can be readily isolated from gardens and starter-cultures of Atta. The accumulated phylogenetic evidence is inconsistent with prevailing views of specific coevolution between Pseudonocardia, attine ants, and garden diseases. Because of frequent acquisition, current models of Pseudonocardia-disease coevolution now need to be revised. The effectiveness of Pseudonocardia antibiotics may not derive from advantages in the coevolutionary arms race with specialized garden diseases, as currently believed, but from frequent recruitment of effective microbes from environmental sources. Indeed, the exposed integumental structures that support actinomycete growth on attine ants argue for a morphological design facilitating bacterial recruitment. We review the accumulated evidence that attine ants have undergone modifications in association with actinomycete bacteria, but we find insufficient support for the reverse, modifications of the bacteria resulting from the interaction with attine ants. The defining feature of coevolution—reciprocal modification—therefore remains to be established for the attine ant-actinomycete mutualism.

Generalization is at the heart of many aspects of behavioral ecology; for foragers it can be seen as an essential feature of learning about potential prey, because natural populations of prey are unlikely to be perfectly homogenous. Aposematic signals are considered to aid predators in learning to avoid a class of defended prey. Predators do this by generalizing between the appearance of prey they have previously sampled and the appearance of prey they subsequently encounter. Mimicry arises when such generalization occurs between individuals of different species. Our aim here is to explore whether the specific shape of the generalization curve can be expected to be important for theoretical predictions relating to the evolution of aposematism and mimicry. We do this by a reanalysis and development of the models provided in two recent papers. We argue that the shape of the generalization curve, in combination with the nature of genetic and phenotypic variation in prey traits, can have evolutionary significance under certain delineated circumstances. We also demonstrate that the process of gradual evolution of Müllerian mimicry proposed by Fisher is particularly efficient in populations with a rich supply of standing genetic variation in mimetic traits.

Small populations of self-incompatible plants are assumed to be threatened by a limitation of compatible mating partners due to low genetic diversity at the self-incompatibility (S) locus. In contrast, we show by using a PCR-RFLP approach for S-genotype identification that 15 small populations (N= 8–88) of the rare wild pear (Pyrus pyraster) displayed no mate limitation. S-allele diversity within populations was high (N= 9–21) as was mate availability (92.9–100%). Although population size and S-allele diversity were strongly related, no relationship was found between population size and mate availability, gene diversity (H_e), or fixation index (F_IS), based on five neutral microsatellite loci. As we determined the principal mate availability within populations based on the S-genotypes observed, the realized mate availability under natural conditions may differ from our estimates, for example, due to spatially limited pollen dispersal. We therefore urge studies on self-incompatible plants to proceed from the simple assessment of principal mate availability to the determination of realized mate availability in natural populations.

A growing number of studies highlight the nontransitive properties of ejaculates when they are in competition to fertilize a female's eggs. Increasingly, these studies suggest that postcopulatory processes act as a filter against sperm from closely related males or those with similar genotypes, limiting the deleterious effects of inbreeding on offspring fitness. We investigated the potential for such postcopulatory mechanisms of inbreeding avoidance in the guppy (Poecilia reticulata), a promiscuous livebearing fish. We used artificial insemination as a method of delivering to a female the combined ejaculates from a first cousin (relatedness coefficient r = 0.125) and an unrelated male. This method of sperm delivery controls behavioral processes of pre- and postcopulatory female choice, which can bias paternity toward unrelated males. Our genetic analysis revealed no effect of parental relatedness on paternity outcomes. The observed mean paternity share for related males (0.47) and associated variance did not differ significantly from an expected binomial distribution that assumes no biased use of sperm with respect to relatedness (0.5). Although our data provide no evidence for postcopulatory mechanisms of inbreeding avoidance, the ability of female guppies to influence ejaculate transfer and retention offers an alternative and easily testable mechanism of inbreeding avoidance in this species.

This study investigates genes enriched for expression in the spermatheca, the long-term sperm storage organ (SSO) of female Drosophila. SSO genes are likely to play an important role in processes of sexual selection such as sperm competition and cryptic female choice. Although there is keen interest in the mechanisms of sexual selection at the molecular level, very little is known about the female genes that are involved. In the present study, a high proportion of genes enriched for expression in the spermatheca are evolving rapidly. Most of the rapidly evolving genes are proteases and genes of unknown function that could play a specialized role in the spermatheca. A high percentage of the rapidly evolving genes have secretion signals and thus could encode proteins that directly interact with ejaculate proteins and coevolve with them. In addition to identifying rapidly evolving genes, the present study documents categories of genes that could play a role in spermatheca function such as storing, maintaining, and utilizing sperm. In general, candidate genes discovered in this study could play a key role in sperm competition, cryptic female choice of sperm, and sexually antagonistic coevolution, and ultimately speciation.