How does the process of life-history evolution interplay with population dynamics? Almost all models that have addressed this question assume that any combination of phenotypic traits uniquely determine the ecological population state. Here we show that if multiple ecological equilibria can exist, the evolution of a trait that relates to competitive performance can undergo adaptive reversals that drive cyclic alternation between population equilibria. The occurrence of evolutionary reversals requires neither environmentally driven changes in selective forces nor the coevolution of interactions with other species. The mechanism inducing evolutionary reversals is twofold. First, there exist phenotypes near which mutants can invade and yet fail to become fixed; although these mutants are eventually eliminated, their transitory growth causes the resident population to switch to an alternative ecological equilibrium. Second, asymmetrical competition causes the direction of selection to revert between high and low density. When ecological conditions for evolutionary reversals are not satisfied, the population evolves toward a steady state of either low or high abundance, depending on the degree of competitive asymmetry and environmental parameters. A sharp evolutionary transition between evolutionary stasis and evolutionary reversals and cycling can occur in response to a smooth change in ecological parameters, and this may have implications for our understanding of size-abundance patterns.

Sexually transmitted diseases (STDs) have been shown to increase the costs of multiple mating and therefore favor relatively monogamous mating strategies. We examine another way in which STDs can influence mating systems in species in which female choice is important. Because more popular males are more likely to become infected, STDs can counteract any selective pressure that generates strong mating skews. We build two models to investigate female mate choice when the sexual behavior of females determines the prevalence of infection in the population. The first model has no explicit social structure. The second model considers the spatial distribution of matings under social monogamy, when females mated to unattractive males seek extrapair fertilizations from attractive males. In both cases, the STD has the potential to drastically reduce the mating skew. However, this reduction does not always happen. If the per contact transmission probability is low, the disease dies out and is of no consequence. In contrast, if the transmission probability is very high, males are likely to be infected regardless of their attractiveness, and mating with the most attractive males imposes again no extra cost for the female. We also show that optimal female responses to the risk of STDs can buffer the prevalence of infection to remain constant, or even decrease, with increasing per contact transmission probabilities. In all cases considered, the feedback between mate choice strategies and STD prevalence creates frequency-dependent fitness benefits for the two alternative female phenotypes considered (choosy vs. randomly mating females or faithful vs. unfaithful females). This maintains mixed evolutionarily stable strategies or polymorphisms in female behavior. In this way, a sexually transmitted disease can stabilize the populationwide proportion of females that mate with the most attractive males or that seek extrapair copulations.

The most enigmatic sexual manipulation by Wolbachia endosymbionts is cytoplasmic incompatibility (CI): infected males are reproductively incompatible with uninfected females. In this paper, we extend the theory on population dynamics and evolution of CI, with emphasis on haplodiploid species. First, we focus on the problem of the threshold to invasion of the Wolbachia infection in a population. Simulations of the dynamics of infection in small populations show that it does not suffice to assume invasion by drift alone (or demographic “accident”). We propose several promising alternatives that may facilitate invasion of Wolbachia in uninfected populations: sex-ratio effects, meta population structure, and other fitness-compensating effects. Including sex-ratio effects of Wolbachia allows invasion whenever infected females produce more infected daughters than uninfected females produce uninfected daughters. Several studies on haplodiploid species suggest the presence of such sex-ratio effects. The simple metapopulation model we analyzed predicts that, given that infecteds are better “invaders,” uninfecteds must be better “colonizers” to maintain coexistence of infected and uninfected patches. This condition seems more feasible for species that suffer local extinction due to predation (or parasitization) than for species that suffer local extinction due to overexploiting their resource(s). Finally, we analyze the evolution of CI in haplodiploids once a population has been infected. Evolution does not depend on the type of CI (female mortality or male production), but hinges solely on decreasing the fitness cost and/or increasing the transmission efficiency. Our models offer new perspectives for increasing our understanding of the population and evolutionary dynamics of CI.

Ciliates provide a powerful system to analyze the evolution of duplicated α-tubulin genes in the context of single-celled organisms. Genealogical analyses of ciliate α-tubulin sequences reveal five apparently recent gene duplications. Comparisons of paralogs in different ciliates implicate differing patterns of substitutions (e.g., ratios of replacement/synonymous nucleotides and radical/conservative amino acids) following duplication. Most substitutions between paralogs in Euplotes crassus, Halteria grandinella and Paramecium tetraurelia are synonymous. In contrast, α-tubulin paralogs within Stylonychia lemnae and Chilodonella uncinata are evolving at significantly different rates and have higher ratios of both replacement substitutions to synonymous substitutions and radical amino acid changes to conservative amino acid changes. Moreover, the amino acid substitutions in C. uncinata and S. lemnae paralogs are limited to short stretches that correspond to functionally important regions of the α-tubulin protein. The topology of ciliate α-tubulin genealogies are inconsistent with taxonomy based on morphology and other molecular markers, which may be due to taxonomic sampling, gene conversion, unequal rates of evolution, or asymmetric patterns of gene duplication and loss.

Anadenobolus excisus is a large species of millipede endemic to the Caribbean Island of Jamaica. Initial detailed morphological studies showed little or no discrete variation across this species' distribution in somatic or, in particular, genitalic morphology. However, a molecular survey based on ∼1000 base pairs of the mitochondrial (mtDNA) 16S rRNA gene that examines 242 individuals sampled from 54 localities reveals three highly divergent mtDNA lineages. A lack of discrete morphological differentiation suggests that genetic and morphological divergence may be decoupled, a pattern inconsistent with a number of evolutionary models. In contrast to minimal morphological divergence, size variation among mtDNA lineages suggests that character displacement has occurred and that these lineages are cohesive in sympatry. We conclude that A. excisus is actually a complex of three cryptic species and that morphological approaches to delineating millipede species may sometimes underestimate evolutionary diversity.

The general purpose of the paper is to test evolutionary optimality theories with experimental data on reproduction, energy consumption, and longevity in a particular Drosophila genotype. We describe the resource allocation in Drosophila females in terms of the oxygen consumption rates devoted to reproduction and to maintenance. The maximum ratio of the component spent on reproduction to the total rate of oxygen consumption, which can be realized by the female reproductive machinery, is called metabolic reproductive efficiency (MRE). We regard MRE as an evolutionary constraint. We demonstrate that MRE may be evaluated for a particular Drosophila phenotype given the fecundity pattern, the age-related pattern of oxygen consumption rate, and the longevity. We use a homeostatic model of aging to simulate a life history of a representative female fly, which describes the control strain in the long-term experiments with the Wayne State Drosophila genotype. We evaluate the theoretically optimal trade-offs in this genotype. Then we apply the Van Noordwijk-de Jong resource acquisition and allocation model, Kirkwood's disposable soma theory, and the Partridge-Barton optimality approach to test if the experimentally observed trade-offs may be regarded as close to the theoretically optimal ones. We demonstrate that the two approaches by Partridge-Barton and Kirkwood allow a positive answer to the question, whereas the Van Noordwijk-de Jong approach may be used to illustrate the optimality. We discuss the prospects of applying the proposed technique to various Drosophila experiments, in particular those including manipulations affecting fecundity.

We report an assay of egg-to-adult viability in full-sibling mutation accumulation (MA) lines derived from a completely homozygous population of Drosophila melanogaster and maintained for 210 generations. A simultaneous evaluation was also made of a large population derived from the same origin and maintained as a control for the same period. We also present computer simulations to explore the possible decline in viability of the control population due to mutation accumulation and the possible effect of selection within and between MA lines. For this purpose, we used two mutational models independent from the data analyzed and based on radically different assumptions. The first model implies a large number of mutations of small effect, whereas the second implies a much smaller number of mutations with much larger effects. The observed rate of decline in mean viability was very small but significant (0.077%). The rate of increase in among line variance (0.189 × 10⁻³) was similar to those obtained previously in the same lines. The simulation results indicated that a model of many mutations of small effect is incompatible with the evolution of the mean viability of the control and MA lines over generations, the distribution of line means after 210 generations of mutation accumulation, and the pattern of line extinction over generations. Basically, this model predicted a large drop in viability, both in the control and particularly the MA lines, that is not observed empirically. It also predicted a rate of line extinction too low in the early generations and too high in the later ones. In contrast, the model based on few mutations of large effect was generally consistent with all the observations.

Patterns of clonal diversity in parthenogenetic species can be influenced by occasional sexual reproduction, selection due to environmental heterogeneity, and genetic drift. Here we investigate geographic patterns of clonal diversity in the obligate parthenogenetic earth mite species, Penthaleus major, in southeastern Australia. Spatial autocorrelation analysis indicated a relationship between geographic distance and clonal diversity at the sites. However, there was no evidence for an association between distance and clonal composition or the frequency of the two most common clones. Regression analyses indicated that clonal diversity and the frequency of one of the common clonal groups were associated with latitude and/or climatic variables, whereas the other common group was not associated with any of those variables. In four transects (each 250–400 km in length) extending from the center of the distribution of P. major to its margins, there was a consistent and continuous decline in clonal diversity as the margin was approached but no evidence that one clonal group predominated at the margin. This argues against the importance of a general-purpose genotype being favored at margins and suggests that either selection acts to decrease local diversity at sites near the margin or that diversity declines due to drift at these sites.

I present patterns characterizing the evolution of intrinsic postzygotic isolation in Lepidoptera by analyzing data from the literature on genetic distance, strength of hybrid sterility and inviability, biogeography, and natural hybridization. Using genetic distance as a proxy for time, I investigate the time-course of the evolution of postzygotic isolation and the waiting times to particular hybrid fitness problems. The results show that postzygotic isolation increases gradually as species diverge, but that hybrid sterility evolves faster than hybrid inviability. The overwhelming preponderance of female-specific hybrid problems in Lepidoptera shows that Haldane's rule (the preferential sterility or inviability of the heterogametic sex) is well obeyed. Together the rates and patterns characterizing the accumulation of postzygotic isolation allow several tests of the composite theory of Haldane's rule. Interestingly, comparing these data with those from Drosophila reveals that Haldane's rule for sterility evolves as fast (if not faster) in Lepidoptera. Finally, I show that a substantial fraction of sympatric species hybridizes in nature and that the majority of these suffer some level of hybrid sterility or inviability.

The grasshoppers Chorthippus brunneus and Chorthippus jacobsi are easily distinguished by male calling song and the number of stridulatory pegs on the hind femur, and form a mosaic hybrid zone in northern Spain. In this paper, we fit a two-dimensional cline to variation in male calling song characters, which are of particular interest as they are likely to be involved in mate choice by females. As with variation in peg number, local habitat makes only a small contribution in explaining deviations in mean song score from clinal expectations. However, the fitted width of the cline for song characters is significantly narrower than for peg number, suggesting that mating signals may be associated with reduced hybrid fitness in the field and that recombination rates are sufficient to allow clines for different characters to diverge in width. Despite this, estimates for the overall elevation in linkage disequilibrium at the zone center, based on covariance between peg and song characters, reveal a substantial overrepresentation of parental genotypes at the cline center relative to the expectations of a tension zone of similar width. Examination of covariance at individual sites reveals that this inflated estimate of linkage disequilibrium is caused by several sites where the distribution of phenotypes is effectively bimodal. This substantial variation in linkage disequilibrium at the cline center could result from local variation in the strength of assortative mating or selection against hybrids, or may reflect the long-distance colonization of empty habitat from outside the hybrid zone, which would continually create new contacts between parental genotypes at the cline center. Hybrid zones like this, in which strong linkage disequilibrium occurs in some situations but not in others, are of particular relevance to speciation research and allow investigation of the spread of combinations of alleles through different genetic and ecological backgrounds.

We investigated the interplay between natural selection and gene flow in the adaptive divergence of threespine stickleback (Gasterosteus aculeatus) that reside parapatrically in lakes and streams. Within the Misty Lake system (Vancouver Island, British Columbia), stickleback from the inlet stream (flowing into the lake) have fewer gill rakers and deeper bodies than stickleback from the lake—differences thought to facilitate foraging (benthic macroinvertebrates in the stream vs. zooplankton in the open water of the lake). Common-garden experiments demonstrated that these differences have a genetic basis. Reciprocal transplant enclosure experiments showed that lake and inlet stickleback grow best in their home environments (although differences were subtle and often not significant). Release-recapture experiments in the inlet showed that lake fish are less well-suited than inlet fish for life in the stream (higher mortality or emigration in lake fish). Morphological divergence in the wild and under common rearing was greater between the lake and the inlet than between the lake and the outlet. Genetic divergence (mitochondrial DNA and microsatellites) was greatest between the lake and the upper inlet (1.8 km upstream from the lake), intermediate between the lake and the lower inlet (0.9 km upstream), and least between the lake and the outlet stream (1.2 km downstream). Relative levels of gene flow estimated from genetic data showed the inverse pattern. The negative association between morphological divergence and gene flow is consistent with the expectation that gene flow can constrain adaptation. Estimated absolute levels of gene flow also implied a constraint on adaptation in the outlet but not the inlet. Our results suggest that natural selection promotes the adaptive divergence of lake and stream stickleback, but that the magnitude of divergence can be constrained by gene flow.

Causal evidence linking resource competition to species divergence is scarce. In this study, we coupled field observations with experiments to ask if the degree of character displacement reflects the intensity of competition between two closely related spadefoot toads (Spea bombifrons and S. multiplicata). Tadpoles of both species develop into either a small-headed omnivorous morph, which feeds mostly on detritus, or a large-headed carnivorous morph, which specializes on and whose phenotype is induced by fairy shrimp. Previously, we found that S. multiplicata are inferior competitors for fairy shrimp and are less likely to develop into carnivores in sympatry with S. bombifrons. We compared four key trophic characters in S. multiplicata across natural ponds where the frequency of S. bombifrons varied. We found that S. multiplicata became increasingly more omnivore-like as the relative abundance of S. bombifrons increased. Moreover, in controlled laboratory populations, S. multiplicata became increasingly more omnivore-like and S. bombifrons became increasingly more carnivore-like as we increased the relative abundance of the other species. Phenotypic plasticity helped mediate this divergence: S. multiplicata became increasingly less likely to eat shrimp and develop into carnivores in the presence of S. bombifrons, a superior predator on shrimp. However, divergence also reflected differences in canalized traits: When reared under common conditions, S. multiplicata tadpoles became increasingly less likely to produce carnivores as their natal pond decreased in elevation. Presumably, this pattern reflected selection against carnivores in lower-elevation ponds, because S. bombifrons became increasingly more common with decreasing elevation. Local genetic adaptation to the presence of S. bombifrons was remarkably fine grained, with differences in carnivore production detected between populations a few kilometers apart. Our results suggest that the degree of character displacement potentially reflects the intensity of competition between interacting species and that both phenotypic plasticity and fine-scale genetic differentiation can mediate this response. Moreover, these results provide causal evidence linking resource competition to species divergence.

Understanding the fitness consequences of inbreeding (inbreeding depression) is of importance to evolutionary and conservation biology. There is ample evidence for inbreeding depression in captivity, and data from wild populations are accumulating. However, we still lack a good quantitative understanding of inbreeding depression and what influences its magnitude in natural populations. Specifically, the relationship between the magnitude of inbreeding depression and environmental severity is unclear. We quantified inbreeding depression in survival and reproduction in populations of cactus finches (Geospiza scandens) and medium ground finches (Geospiza fortis) living on Isla Daphne Major in the Galápagos Archipelago. Our analyses showed that inbreeding strongly reduced the recruitment probability (probability of breeding given that an adult is alive) in both species. Additionally, in G. scandens, first-year survival of an offspring with f = 0.25 was reduced by 21% and adults with f = 0.25 experienced a 45% reduction in their annual probability of survival. The magnitude of inbreeding depression in both adult and juvenile survival of this species was strongly modified by two environmental conditions, food availability and number of competitors. In juveniles, inbreeding depression was only present in years with low food availability, and in adults inbreeding depression was five times more severe in years with low food availability and large population sizes. The combination of relatively severe inbreeding depression in survival and the reduced recruitment probability led to the fact that very few inbred G. scandens ever succeeded in breeding. Other than recruitment probability, no other trait showed evidence of inbreeding depression in G. fortis, probably for two reasons: a relatively high rate of extrapair paternity (20%), which may lead to an underestimate of the apparent inbreeding depression, and low sample sizes of highly inbred G. fortis, which leads to low statistical power. Using data from juvenile survival, we estimated the number of lethal equivalents carried by G. scandens, G. fortis, and another congener, G. magnirostris. These results suggest that substantial inbreeding depression can exist in insular populations of birds, and that the magnitude of the inbreeding depression is a function of environmental conditions.

Despite the importance of Darwin's finches to the development of evolutionary theory, the origin of the group has only recently been examined using a rigorous, phylogenetic methodology that includes many potential outgroups. Knowing the evolutionary relationships of Darwin's finches to other birds is important for understanding the context from which this adaptive radiation arose. Here we show that analysis of mitochondrial DNA sequence data from the cytochrome b gene confirm that Darwin's finches are monophyletic. In addition, many taxa previously proposed as the sister taxon to Darwin's finches can be excluded as their closest living relative. Darwin's finches are part of a well-supported monophyletic group of species, all of which build a domed nest. All but two of the non-Darwin's finches included in this clade occur on Caribbean islands and most are Caribbean endemics. These close relatives of Darwin's finches show a diversity of bill types and feeding behaviors similar to that observed among Darwin's finches themselves. Recent studies have shown that adaptive evolution in Darwin's finches occurred relatively quickly. Our data show that among the relatives of Darwin's finches, the evolution of bill diversity was also rapid and extensive.

Efforts to evaluate the evolutionary and ecological dynamics of conspecific brood parasitism in birds and other animals have focused on the fitness costs of parasitism to hosts and fitness benefits to parasites. However, it has been speculated recently that, in species with biparental care, host males might cooperate with parasitic females by allowing access to the host nest in exchange for copulations. We develop a cost-benefit model to explore the conditions under which such host-parasite cooperation might occur. When the brood parasite does not have a nest of her own, the only benefit to the host male is siring some of the parasitic eggs (quasi-parasitism). Cooperation with the parasite is favored when the ratio of host male paternity of his own eggs relative to his paternity of parasitic eggs exceeds the cost of parasitism. When the brood parasite has a nest of her own, a host male can gain additional, potentially more important benefits by siring the high-value, low-cost eggs laid by the parasite in her own nest. Under these conditions, host males should be even more likely to accept parasitic eggs in return for copulations with the parasitic female. We tested these predictions for American coots (Fulica americana), a species with a high frequency of conspecific brood parasitism. Multilocus DNA profiling indicated that host males did not sire any of the parasitic eggs laid in host nests, nor did they sire eggs laid by the parasite in her own nest. We used field estimates of the model parameters from a four-year study of coots to predict the minimum levels of paternity required for the costs of parasitism to be offset by the benefits of mating with brood parasites. Observed levels of paternity were significantly lower than those predicted under a variety of assumptions, and we reject the hypothesis that host males cooperated with parasitic females. Our model clarifies the specific costs and benefits that influence host-parasite cooperation and, more generally, yields precise predictions about expected levels of host male paternity. These predictions will enable a more rigorous assessment of field studies designed to test adaptive hypotheses of host-parasite cooperation.

To study the correlated evolution of locomotor behavior and exercise physiology, we conducted an artificial selection experiment. From the outbred Hsd:ICR strain of Mus domesticus, we began eight separate lines, each consisting of 10 breeding pairs. In four of the lines, we used within-family selection to increase voluntary wheel running. The remaining four lines were random-bred (within lines) to serve as controls. Various traits have been monitored to test for correlated responses. Here, we report on organ masses, with emphasis on the triceps surae muscle complex, an important extensor of the ankle. Mice from the selected lines exhibit reduced total body mass, increased relative (mass-corrected) kidney mass, and reduced relative triceps surae mass. In addition, a discrete muscle-mass polymorphism was observed: some individuals had triceps surae that were almost 50% lighter than normal for their body mass. This small-muscle phenotype was observed in only three of the eight lines: in one control line, it has fluctuated in frequency between zero and 10%, whereas in two of the selected lines it has increased in frequency to approximately 50% by generation 22. Data from a set of parents and offspring (generations 23 and 24) are consistent with inheritance as a single autosomal recessive allele. Evidence for the adaptive significance of the small-muscle allele was obtained by fitting multiple-generation data to hierarchical models that include effects of genetic drift and/or selection. The small-muscle allele is estimated to have been present at low frequency (approximately 7%) in the base population, and analysis indicates that strong selection favors the allele in the selected but not control lines. We hypothesize that the small muscles possess functional characteristics and/or that the underlying allele causes pleiotropic effects (e.g., reduced total body mass; increased relative heart, liver, and kidney mass) that facilitate high levels of wheel running. Nevertheless, at generation 22, wheel running of affected individuals did not differ significantly from those with normal-sized muscles, and the magnitude of response to selection has been similar in all four selected lines, indicating that multiple genetic “solutions” are possible in response to selection for high activity levels.

Jarman (1974) proposed a series of relationships between habitat use, food dispersion, and social behavior and hypothesized a series of evolutionary steps leading to sexual dimorphism in body size through sexual selection in African antelope species. The hypothesis states that sexual size dimorphism evolved in a three-step process. Initially, ancestral monomorphic and monogamous ungulate species occupying closed habitats radiated into open grassland habitats. Polygynous mating systems then rapidly evolved in response to the aggregation of males and females, perhaps in relation to the clumped distribution of food resources in open habitats. Subsequently, size dimorphism evolved in those species occupying open habitats, but not in species that remained in closed habitats or retained monogamy. This hypothesis has played an important role in explaining the origins of sexual dimorphism in mammals. However, the temporal sequence of the events that Jarman proposed has never been demonstrated. Here we use a phylogeny of extant ungulate species, along with maximum-likelihood statistical techniques, to provide a test of Jarman's hypothesis.

There is considerable interest in comparing genetic variance-covariances matrices (G matrix). However, present methods are difficult to implement and cannot readily be extended to incorporate effects of other variables such as habitat, sex, or location. In this paper I present a method based on MANOVA that can be done using only standard statistical packages (coding for the method using SPLUS is available from the author). The crux of the approach is to use the jackknife method to estimate the pseudovalues of the estimates; these estimates can then be used as datapoints in a MANOVA. I illustrate the method using two published datasets: (1) variation in G matrices resulting from differences in rearing condition, species, and sex in the crickets Gryllus firmus and G. pennsylvanicus; and (2) variation in G matrices associated with habitat and history in the amphipod Gammarus minus.

Phylogenetic tree imbalance was originally believed to indicate differences in evolutionary rates within trees, but other sources of imbalance have been identified, such as tree incompleteness and low quality of the data. To examine the effect of data quality, I calculated Colless's index for 69 recent complete phylogenies. On average, these phylogenies were more unbalanced than phylogenies generated by the equal rates Markov (ERM) model. I tried Mooers's (1995) method to correct for tree size, but his measure appeared to become dependent on tree size when there are large trees (i.e., > 14 tips) in a collection. Instead I corrected for tree size by taking the difference between Colless's index of observed trees and the ERM model expectation for a tree of the same size. The balance measure thus obtained did not correlate significantly to consistency and retention indices as indicators of data quality. It was also independent of the factors kingdom (plants and animals) and taxon level at the tips and type of data (molecular, morphological, and combined).

Complex hybridity is a rare diploid genetic system of plants, extensively characterized in Oenothera, in which heterozygosity for one or more reciprocal translocations is maintained by means of autogamy and a balanced lethal system. It is visible at metaphase I of meiosis as rings or chains of chromosomes held together by terminalized chiasmata. Phylogenetic analysis based on 274 random amplified polymorphic DNA markers showed that in the Australian endemic Isotoma petraea (Lobeliaceae) the genetic system had a single origin, as a ring-of-six, in the Pigeon Rock population. It subsequently spread to other populations to produce hybrids incorporating additional chromosomes into the rings. Our research supports the suggestion that complex hybridity is an evolutionary response to intense inbreeding and selected because it allows masking of accumulated deleterious alleles.

An explicit assumption of studies that employ a mitochondrial DNA (mtDNA) molecular clock is that mtDNA evolves independently of morphology. Here we report a very strong correlation between egg size divergence and cytochrome c oxidase-1 (CO1) amino acid sequence divergence among sister species of bivalve molluscs separated by the Central American Isthmus (i.e., “geminate” species). Analyses of the molecular data reveal that CO1 sequences likely did not diverge as a function of time or evolve in response to positive natural selection. Given that an excess of CO1 amino acid polymorphism exists within species (as expected if most mutations are only slightly deleterious), a third hypothesis is that reductions in effective population size could simultaneously increase the fixation rate of nearly neutral mtDNA polymorphisms and in some way also facilitate egg size evolution. The remarkable strength of the relationship between egg size and CO1 amino acid sequence demonstrates that, even in the absence of an obvious functional relationship or clock-like evolution, the amounts of molecular and morphological change can be tightly correlated, and therefore may reflect common processes. Accordingly, the assumption that the evolutionary divergence of molecules and morphology are independent must always be carefully examined.