Evolutionary biologists assume that species formation requires a drastic reduction in gene exchange between populations, but the rate sufficient to prevent speciation is unknown. To study speciation, we use a new class of population genetic models that incorporate simple developmental genetic rules, likely present in all organisms, to construct the phenotype. When we allow replicate populations to evolve in parallel to a new, shared optimal phenotype, often their hybrids acquire poorly regulated phenotypes: Dobzhansky-Muller incompatibilities arise and postzygotic reproductive isolation evolves. Here we show that, although gene exchange does inhibit this process, it is the proportion of migrants exchanged (m) rather than the number of migrants (Nm) that is critical, and rates as high as 16 individuals exchanged per generation still permit the evolution of postzygotic isolation. Stronger directional selection counters the inhibitory effect of gene flow, increasing the speciation probability. We see similar results when populations in a standard two-locus, two-allele Dobzhansky-Muller model are subject to simultaneous directional selection and gene flow. However, in developmental pathway models with more than two loci, gene flow is more able to impede speciation. Genetic incompatibilities arise as frequent by-products of adaptive evolution of traits determined by regulatory pathways, something that does not occur when phenotypes are modeled using the standard, additive genetic framework. Development therefore not only constrains the microevolutionary process, it also facilitates the interactions among genes and gene products that make speciation more likely—even in the face of strong gene flow.

Seed dormancy plays an important role in germination ecology and seed plant evolution. Morphological seed dormancy is caused by an underdeveloped embryo that must mature prior to germination. It has been suggested that the presence of an underdeveloped embryo is plesiomorphic among seed plants and that parallel directional change in embryo morphology has occurred separately in gymnosperms and in angiosperms. We test these hypotheses using original data on embryo morphology of key basal taxa, a published dataset, and the generalized least squares (GLS) method of ancestral character state reconstruction. Reconstructions for embryo to seed ratio (E:S) using family means for 179 families showed that E:S has increased between the ancestral angiosperm and almost all extant angiosperm taxa. Species in the rosid clade have particularly large embryos relative to the angiosperm ancestor. Results for the gymnosperms show a similar but smaller increase. There were no statistically significant differences in E:S between basal taxa and any derived group due to extremely large standard errors produced by GLS models. However, differences between reconstructed values for the angiosperm ancestor and more highly nested nodes are large and these results are robust to topological and branch-length manipulations. Our analysis supports the idea that the underdeveloped embryo is primitive among seed plants and that there has been a directional change in E:S within both angiosperms and gymnosperms. Our analysis suggests that dormancy enforced by an underdeveloped embryo is plesiomorphic among angiosperms and that nondormancy and other dormancy types probably evolved within the angiosperms. The shift in E:S was likely a heterochronic change, and has important implications for the life history of seed plants.

Experimental and comparative evidence implies that homoploid hybrid speciation is a reproducible process, mediated in part by ecological selection. Here, molecular data from the chloroplast genome and 17 nuclear microsatellite loci were employed to determine whether a well-documented homoploid hybrid species, Helianthus paradoxus, has arisen multiple times. Helianthus paradoxus is ecologically divergent from its parental species, and has a disjunct geographic distribution consistent with multiple origins. The molecular data, however, strongly support a single hybrid origin. First, all sampled populations of H. paradoxus are fixed for a single chloroplast DNA (cpDNA) haplotype, whereas local populations of both parental species, H. annuus and H. petiolaris, have multiple cpDNA haplotypes. Second, H. paradoxus populations form a single, well-supported clade (99.8% bootstrap support) in a neighbor-joining tree based on microsatellite allele frequencies. The microsatellite data also tentatively place the origin of H. paradoxus between 75,000 years and 208,000 years before present, indicating that anthropogenic disturbance likely did not play a role in the formation of this species. Finally, the genetic structure of this species is not consistent with passive riparian dispersal, which has been suggested for other wetland plant species, but may be explained by dispersal mechanisms implicated for H. annuus, such as large migratory mammals.

The genetic basis of species differences provides insight into the mode and tempo of phenotypic divergence. We investigate the genetic basis of floral differences between two closely related plant taxa with highly divergent mating systems, Mimulus guttatus (large-flowered outcrosser) and M. nasutus (small-flowered selfer). We had previously constructed a framework genetic linkage map of the hybrid genome containing 174 markers spanning approximately 1800 cM on 14 linkage groups. In this study, we analyze the genetics of 16 floral, reproductive, and vegetative characters measured in a large segregating M. nasutus × M. guttatus F₂ population (N = 526) and in replicates of the parental lines and F₁ hybrids. Phenotypic analyses reveal strong genetic correlations among floral traits and epistatic breakdown of male and female fertility traits in the F₂ hybrids. We use multitrait composite interval mapping to jointly locate and characterize quantitative trait loci (QTLs) underlying interspecific differences in seven floral traits. We identified 24 floral QTLs, most of which affected multiple traits. The large number of QTLs affecting each trait (mean = 13, range = 11–15) indicates a strikingly polygenic basis for floral divergence in this system. In general, QTL effects are small relative to both interspecific differences and environmental variation within genotypes, ruling out QTLs of major effect as contributors to floral divergence between M. guttatus and M. nasutus. QTLs show no pattern of directional dominance. Floral characters associated with pollinator attraction (corolla width) and self-pollen deposition (stigma-anther distance) share several pleiotropic or linked QTLs, but unshared QTLs may have allowed selfing to evolve independently from flower size. We discuss the polygenic nature of divergence between M. nasutus and M. guttatus in light of theoretical work on the evolution of selfing, genetics of adaptation, and maintenance of variation within populations.

Crucial to understanding the process of natural selection is characterizing phenotypic selection. Measures of phenotypic selection can be biased by environmental variation among individuals that causes a spurious correlation between a trait and fitness. One solution is analyzing genotypic data, rather than phenotypic data. Genotypic data, however, are difficult to gather, can be gathered from few species, and typically have low statistical power. Environmental correlations may act through traits other than through fitness itself. A path analytic framework, which includes measures of such traits, may reduce environmental bias in estimates of selection coefficients. We tested the efficacy of path analysis to reduce bias by re-analyzing three experiments where both phenotypic and genotypic data were available. All three consisted of plant species (Impatiens capensis, Arabidopsis thaliana, and Raphanus sativus) grown in experimental plots or the greenhouse. We found that selection coefficients estimated by path analysis using phenotypic data were highly correlated with those based on genotypic data with little systematic bias in estimating the strength of selection. Although not a panacea, using path analysis can substantially reduce environmental biases in estimates of selection coefficients. Such confidence in phenotypic selection estimates is critical for progress in the study of natural selection.

Allocation to sexual reproduction is an important life-history trait in clonal plants. Different selection pressures between competitive and competition-free environments are likely to result in the evolution of specialized genotypes and to maintain genetic variation in reproductive allocation. Moreover, selection may also result in the evolution of plastic allocation strategies. The necessary prerequisite for evolution, heritable genetic variation, can best be studied with selection experiments. Starting from a base population of 102 replicated genotypes of the clonal herb Ranunculus reptans, we imposed selection on the proportion of flowering rosettes in the absence of competition (base population: mean = 0.391, broad-sense heritability = 0.307). We also selected on the plasticity in this trait in response to competition with a naturally coexisting grass in a parallel experiment (base population: 14% lower mean in the presence of competition, broad-sense heritability = 0.072). After two generations of bidirectional selection, the proportion of flowering rosettes was 26% higher in the high line than in the low line (realized heritability ± SE = 0.205 ± 0.017). Moreover, genotypes of the high line had 11% fewer carpels per flower, a 22% lower proportion of rooted rosettes, and a 39% smaller average distance between rosettes within a clone. In the second experiment, we found no significant responses to selection for high and low plasticity in the proportion of flowering rosettes (realized heritability ± SE = −0.002 ± 0.013). Our study indicates a high heritability and potential for further evolution of the proportion of flowering rosettes in R. reptans, but not for its plasticity, which may have been fixed by past evolution at its current level. Moreover, our results demonstrate strong genetic correlations between allocation to sexual reproduction and other clonal life-history characteristics.

Models allowing the coexistence of females and hermaphrodites in gynodioecious populations assume a simple genetic system of sex determination, a seed fitness advantage of females (compensation), and a negative pleiotropic effect of nuclear sex-determining genes on fitness (cost of restoration). In Lobelia siphilitica, sex is determined by both mitochondrial genes causing cytoplasmic male sterility (CMS) and nuclear genes that restore fertility when present with specific CMS haplotypes (nuclear restorers). I tested for a cost of restoration in L. siphilitica by measuring restored hermaphrodites for five fitness components and estimating the number of nuclear restorers by crosses with females carrying CMS1 and CMS2. A cost of restoration appears as a significant negative coefficient (B) in the regression model explaining fitness. I found that hermaphrodites carrying more nuclear restorer genes for CMS2 (or restorer genes of greater effect) have lower pollen viability (B = −1.08, P = 0.001). This pollen viability cost of restoration in L. siphilitica supports the theoretical prediction that negative pleiotropic effects of restorers will exist in populations of gynodioecious species containing females. The existence of such a cost supports the view that gynodioecy can be a stable breeding system in nature.

Genetically based variation in outcrossing rate generates lineages within populations that differ in their history of inbreeding. According to some models, mating-system modifiers in such populations will demonstrate both linkage and identity disequilibrium with fitness loci, resulting in lineage-specific inbreeding depression. Other models assert that differences among families in levels of inbreeding depression are mainly attributable to random accumulation of genetic load, unrelated to variation at mating-system loci. We measured female reproductive success of selfed and outcrossed progeny from naturally occurring lineages of Datura stramonium, a predominantly self-fertilizing annual weed that has heritable variation in stigma-anther separation, a trait that influences selfing rates. Progeny from inbred lineages (as identified by high degree of anther-stigma overlap) showed equal levels of seed production, regardless of cross type. Progeny from mixed lineages (as identified by relatively high separation between anthers and stigma) showed moderate levels of inbreeding depression. We found a significant correlation between anther-stigma separation and relative fitness of selfed and outcrossed progeny, suggesting that family-level inbreeding depression may be related to differences among lineages in inbreeding history in this population. Negative inbreeding depression in putatively inbred lineages may be due in part to additive effects or to epistatic interactions among loci.

Black mustard, Brassica nigra, is highly variable in both constitutive resistance and induction response following damage by herbivores. A focal population from Ithaca, New York, was used to test the following two predictions of optimal defense theory: (1) that allocation to resistance will reduce plant performance in the absence of herbivores; and (2) that induction response will be negatively correlated with constitutive resistance. The experiment consisted of a half-sib mating design with 47 paternal families and four dams per sire, fully crossed with a damage treatment consisting of 25% leaf removal by the cabbage white butterfly, Pieris rapae, when plants had four leaves. Leaf trichome density, sinigrin concentration, and glucobrassicin concentration were 38%, 19%, and 16% higher, respectively, for the seventh leaf of damaged plants. Paternal families did not vary significantly in their induction response. Narrow-sense heritabilities were h²_s = 0.51, 0.76, and 0.50 for constitutive leaf trichome density, sinigrin concentration, and glucobrassicin concentration, respectively. Positive genetic correlations were found between glucobrassicin concentration and days to first flower, suggesting a genetic cost of resistance. Induction responses were negatively correlated with constitutive allocation for leaf trichome density and sinigrin concentration. The results were therefore consistent with optimal defense theory, offering modest evidence for both predictions.

Adaptive phenotypic plasticity in chemical defense is thought to play a major role in plant-herbivore interactions. We investigated genetic variation for inducibility of defensive traits in wild radish plants and asked if the evolution of induction is constrained by costs of phenotypic plasticity. In a greenhouse experiment using paternal half-sibling families, we show additive genetic variation for plasticity in glucosinolate concentration. Genetic variation for glucosinolates was not detected in undamaged plants, but was significant following herbivory by a specialist herbivore, Pieris rapae. On average, damaged plants had 55% higher concentrations of glucosinolates compared to controls. In addition, we found significant narrow-sense heritabilities for leaf size, trichome number, flowering phenology, and lifetime fruit production. In a second experiment, we found evidence of genetic variation in induced plant resistance to P. rapae. Although overall there was little evidence for genetic correlations between the defensive and life-history traits we measured, we show that more plastic families had lower fitness than less plastic families in the absence of herbivory (i.e., evidence for genetic costs of plasticity). Thus, there is genetic variation for induction of defense in wild radish, and the evolution of inducibility may be constrained by costs of plasticity.

Distinguishing the historical effects of gene migration and vicariance on contemporary genetic structure is problematic without testable biogeographic hypotheses based on preexisting geological and environmental evidence. The availability of such hypotheses for North America's Sonoran Desert has contributed to our understanding of the effect of historical vicariance and dispersal events on the diversification of this region's vertebrate biota but have not yet been applied to its flora. In this paper we describe a detailed allozyme analysis of the population genetic structure and phylogeography of the Sonoran Desert columnar cactus, Lophocereus schottii (senita). Inferred phylogroup distributions reflect two historical vicariance events: (1) a middle Pliocene northward transgression of the Sea of Cortéz that is reflected in well-supported Baja California peninsular and continental phylogroups but not in current taxonomic treatments of the species; and (2) a late Pliocene transpeninsular seaway across southern Baja that is reflected in tentative support for peninsular and southern Cape Region phylogroups corresponding to taxonomic varieties L. schottii var. schottii and L. schottii var. australis, respectively. A middle Pleistocene midpeninsular seaway hypothesized to explain congruent phylogroup distributions in several vertebrate taxa is not reflected in L. s. var. schottii, nor is the distinction of a third variety, L. s. var. tenuis, from continental populations of L. s. var. schottii. Linear regression of pairwise estimates of interpopulation differentiation (M̂ and F_ST/[1 − F_ST]) on interpopulation geographic distance revealed significant evidence of isolation by distance within peninsular and continental phylogroups but not between them, consistent with historical vicariance between but not within these regions. We also found significant evidence of isolation by distance between putative L. s. var. schottii and L. s. var. australis phylogroups, suggesting that reproductive isolation between peninsular and Cape Region forms is incomplete. Within peninsular, but not continental, phylogroups, northward range expansion from southern Pleistocene refugia is reflected in significant declines in genetic variation with increasing latitude and in an area phenogram in which populations are progressively nested from south (ancestral) to north (descendant) along the Baja peninsula. Although the geographic concordance of phylogenetic topologies suggests that ancient vicariance events, and not dispersal, have primarily influenced the biogeographic distributions of Baja's vertebrate biota, the phylogeographic structure of L. schottii suggests that Sonoran Desert plant species may exhibit genetic signatures of postglacial range expansion and gene flow as well as vicariance.

A few marine cases have demonstrated morphological and genetic divergence in the absence of spatial barriers to gene flow, suggesting that the initial phase of speciation is possible without geographic isolation. In the Bocas del Toro Archipelago of the Atlantic Coast of Panama, we found two morphotypes of the scleractinian coral Favia fragum with opposing depth distributions. One morphotype fit the classical description of F. fragum and was most abundant at 3 m depth. A second morphotype was distinguished by raised corallites and was restricted to ≤ 1 m depth. The two morphotypes overlapped in distribution at 1 m depth. Multivariate analysis of polyp-level characters (shape and distribution of septa within corallites) divided samples into two groups corresponding to initial qualitative observations of colony shape and corallite relief. To determine whether reduced gene flow maintains morphological variation, we measured the frequencies of alleles at five allozyme loci in both morphotypes at three sites 1–2 km distant. While there were significant differences in allele frequencies between morphotypes within sites, there were also frequency differences among sites at most loci, with the exception of nearly fixed alleles at the PGM locus. Extremely low heterozygosity permitted us to use haplotypes to compare genetic distance between morphotypes and among sites. Comparisons between haplotype data and a null model assuming gene flow between morphotypes showed that the two morphotypes shared significantly fewer haplotypes than expected, and average genetic distance between morphotypes was significantly greater than expected. Partitioning haplotype variation with analysis of molecular variance demonstrated that 35% of the variation was explained by morphotype, whereas 28% of the variation was explained by site. Two PGM heterozygotes and several individuals homozygous for rare PGM alleles are consistent with hybridization, and perhaps introgression by selfing within morphotypes. We consider three hypotheses for this morphological and genetic divergence in F. fragum: (1) intraspecific polymorphism, (2) incipient species, (3) biological species; and discuss the role of reproductive characters in a divergence-with-gene flow mechanism of speciation.

The maintenance of two genetically distinct reproductive modes such as outcrossing and selfing within a population of animals or plants is still a matter of considerable debate. Hermaphroditic parasites often reproduce either alone by selfing or in pairs by outcrossing. They can be used as a model to study potential benefits of outcrossing. Any advantage from outcrossing may be important, especially in host-parasite coevolution, but has not, to our knowledge, been studied yet in any parasite species. We studied the potential effect of outcrossing in a tapeworm, Schistocephalus solidus, on both infection success and growth in its first intermediate host, the copepod Macrocyclops albidus. Tapeworms that had been obtained from natural populations of three-spined sticklebacks (Gasterosteus aculeatus) were allowed to reproduce either alone or in pairs, in an in vitro system that replaced the final host's gut. This resulted in either selfed or outcrossed offspring, respectively. In one part of the experiment, copepods were exposed to either selfed or outcrossed parasites, in a second part to both types simultaneously, in order to study the effect of competition between them. To discriminate parasites of either origin within the same host, a novel method for fluorescent vital labeling was used. We show here for the first time that outcrossed parasites had a higher infection success and faster development in the host. This advantage of outcrossing became apparent only in the competitive situation, in which superior abilities of parasites to extract limiting resources from the host become crucial.

Doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA) has been demonstrated in both mytilid and unionid bivalves. Under DUI, females pass on their mtDNA to both sons and daughters, whereas males pass on their mtDNA to only sons. In mytilids, the loss of an original male (or M) mitotype, with its subsequent replacement by that lineage's female (or F) mitotype, has been called a role-reversal or, more specifically, a masculinization event. Multiple masculinization events have been inferred during the evolutionary history of mytilids but not unionids. The perceived lack of role-reversal events in unionids may represent a significant difference in the evolutionary dynamics of DUI between the two bivalve taxa or simply a lack of sufficient taxon sampling in unionids. To evaluate these alternative hypotheses, six additional unionoidean bivalve genera were sampled for DUI including one genus from the sister taxon of the Unionidae, the Hyriidae. Phylogenetic analyses of 619 base pairs of cytochrome c oxidase I (COI) from eight genera (nine species) of unionoidean bivalves, plus the sister taxon to the Unionoida, Neotrigonia, revealed that the M and F unionoidean mitotypes were contained in gender-specific, topologically congruent clades. This supports the hypothesis that either role-reversal events do not occur in unionoideans or, if they do occur, their products are ephemeral in an evolutionary sense. Furthermore, the fact that the mantle-tissue-derived Neotrigonia mitotype is the sister mitotype to the unionoidean F mitotype clade suggests that DUI has been operating with high fidelity in unionoids for at least 200 million years. A relatively low incidence of interspecific hybridization in unionoideans and a possibly obligate role for the M mitotype in unionoidean gender determination are offered as potential explanations for the disparate evolutionary dynamics of DUI observed between mytilid and unionoidean bivalves.

Drosophila yakuba is a species widespread in Africa, whereas D. santomea, its newly discovered sister species, is endemic to the volcanic island of São Tomé in the Gulf of Guinea. Drosophila santomea probably formed after colonization of the island by its common ancestor with D. yakuba. The two species differ strikingly in pigmentation: D. santomea, unlike the other eight species in the D. melanogaster subgroup, almost completely lacks dark abdominal pigmentation. D. yakuba shows the sexually dimorphic pigmentation typical of the group: both sexes have melanic patterns on the abdomen, but males are much darker than females. A genetic analysis of this species difference using morphological markers shows that the X chromosome accounts for nearly 90% of the species difference in the area of abdomen that is pigmented and that at least three genes (one on each major chromosome) are involved in each sex. The order of chromosome effects on pigmentation area are the same in males and females, suggesting that loss of pigmentation in D. santomea may have involved the same genes in both sexes. Further genetic analysis of the interspecific difference between males in pigmentation area and intensity using molecular markers shows that at least five genes are responsible, with no single locus having an overwhelming effect on the trait. The species difference is thus oligogenic or polygenic. Different chromosomal regions from each of the two species influenced pigmentation in the same direction, suggesting that the species difference (at least in males) is due to natural or sexual selection and not genetic drift. Measurements of sexual isolation between the species in both light and dark conditions show no difference, suggesting that the pigmentation difference is not an important cue for interspecific mate discrimination. Using DNA sequence differences in nine noncoding regions, we estimate that D. santomea and D. yakuba diverged about 400,000 years ago, a time similar to the divergences between two other well-studied pair of species in the subgroup, both of which also involved island colonization.

Rapid changes in climate may impose strong selective pressures on organisms. Evolutionary responses to climate change have been observed in natural populations, yet no example has been documented for a metabolic enzyme locus. Furthermore, few studies have linked physiological responses to stress with allozyme genotypic variation. We quantified changes in allele frequency between 1988 and 1996 at three allozyme loci (isocitrate dehydrogenase, Idh; phosphoglucose isomerase, Pgi; and phosphoglucomutase, Pgm) for the leaf beetle Chrysomela aeneicollis in the Bishop Creek region of the Sierra Nevada of California (2900–3300 m). Beetles often experience high daytime (>32°C) and extremely low nighttime (<−5°C) temperatures during summer. Bishop Creek weather station data indicated that conditions were unusually dry before 1988, and that conditions were cool and wet during the years preceding the 1996 collection. We found directional changes in allele frequency at Pgi (11% increase in the Pgi-1 allele), but not at Idh or Pgm. We also found that physiological response to thermal extremes depended on Pgi genotype. Pgi 1–1 individuals induced expression of a 70-kD heat shock protein (HSP) at lower temperatures than 1–4 or 4–4 individuals, and 1–1 individuals expressed higher levels of HSP70 after laboratory exposure to temperatures routinely experienced in nature. Survival after nighttime laboratory exposure to subzero temperatures depended on gender, previous exposure to cold, and Pgi genotype. Females expressed higher levels of HSP70 than males after exposure to heat, and recovery by female Pgi 1–1 homozygotes after exposure to cold (−5°C) was significantly better than 1–4 or 4–4 genotypes. These data suggest that the cooler climate of the mid-1990s may have caused an increase in frequency of the Pgi-1 allele, due to a more robust physiological response to cold by Pgi 1–1 and 1–4 genotypes.

Maternally transmitted bacteria that kill male hosts early in their development are found in many insects. These parasites typically infect 1–30% of wild females, but in a few species of insects, prevalences exceed 95%. We investigated one such case in the butterfly Acraea encedon, which is infected with a male-killing Wolbachia bacterium. We measured three key parameters that affect the prevalence of the parasite: transmission efficiency, rate of survival of infected males, and the direct cost of infection. We observed that all wild females transmit the bacterium to all their offspring and that all infected males die in wild populations. We were unable to detect any physiological cost to infection in lab culture. These observations explain the high prevalence of the A. encedon male killer, as theory predicts that under these conditions the parasite will spread to fixation. This will occur provided the death of males provides some benefit to the surviving infected females. The problem therefore becomes why the bacterium has not reached fixation and driven the butterfly extinct due to the shortage of males. We therefore investigated whether males choose to mate with uninfected rather than infected females, as this would prevent the bacterium from reaching fixation. We tested this hypothesis in the “lekking swarms” of virgin females found in the most female-biased populations, and were unable to detect any evidence of mate choice. In conclusion, this male killer has spread to high prevalence because it has a high transmission efficiency and low cost, but the factors maintaining uninfected females in the population remain unknown.

The pattern of character variation within a hybrid zone, the hybrid zone structure, has been used to infer the processes that maintain hybrid zones. Unfortunately it is difficult to infer process from structure alone because many different processes can produce the same pattern of character variation. Mosaic hybrid zones may be maintained by exogenous selection in a heterogeneous environment and/or endogenous selection against hybrid individuals; habitat preference, premating isolating barriers and/or fertility selection can also contribute. The spatial scale at which a hybrid zone is sampled affects its apparent structure; a hybrid zone may appear clinal at one scale and mosaic at another. Here, we sample the mosaic hybrid zone between two field crickets, Gryllus firmus and G. pennsylvanicus, at a scale that spans the boundaries between individual soil-habitat patches. From our analysis, we find that at fine scales, the mosaic hybrid zone resolves into a set of steep clines across patch boundaries. Both morphological and molecular traits exhibit sharp and generally concordant clines. However, clines for mitochondrial DNA and one anonymous nuclear marker are clearly displaced as a result of current hybridization or past introgression (the “ghost of hybridization past”). Thus, scale is important for the structure of this and probably other hybrid zones. The extremely sharp, concordant clines across patch boundaries indicate that the cricket hybrid zone is undoubtedly structured by selection. However, the detailed mechanisms responsible for the maintenance of the hybrid zone—whether endogenous selection against hybrids, exogenous selection by the environment, and/or behavioral preferences for mates or habitats—remain to be elucidated. Determining these mechanisms will depend on closer inspection of the organisms themselves and their interactions, as is the case for all hybrid zones.

Inclusive fitness benefits have been suggested to be a major selective force behind the evolution of cooperative breeding. We investigated the fitness benefits selecting for cooperative breeding in the Seychelles warbler, Acrocephalus sechellensis. A microsatellite-based genotyping method was used to determine the relatedness of subordinates to group offspring in an isolated population of Seychelles warblers. The indirect and direct breeding benefits accruing to individual subordinates were then calculated for every successful breeding event over a three-year period. We show that female subordinates frequently gained parentage and that this, combined with high levels of extragroup paternity, resulted in low levels of relatedness between subordinates and nondescendent offspring within a territory. Direct breeding benefits were found to be significantly higher than indirect kin benefits for both female and male subordinates. As predicted, female subordinates gained significantly more direct breeding opportunities and therefore higher inclusive fitness benefits by being a subordinate within a group than did males. This may explain why most subordinates in the Seychelles warbler are female.

Studies of phenotype-environment associations in adaptive radiation have focused largely on morphological traits related to resource-based phenotypic differences. The genetic basis of adaptive behaviors implicated in population divergence remains poorly understood, as few studies have tested the hypothesis of behavioral phenotype-environment associations. We provide evidence of a phenotype-environment association for differential adaptive swimming behaviors through experiments conducted on dwarf, normal, and hybrid lake whitefish (Coregonus clupeaformis). Highly significant differences were observed for depth selection, directional changes, and burst swimming, implicating a genetic basis for these behaviors. Hybrid crosses revealed that depth selection is under additive genetic control, while dominance effects were suggested for directional changes and burst swimming. Estimates for the genetic basis of behavioral differentiation from an animal model were consistent with these observations. Comparative estimates of behavioral differentiation (Q_ST) against neutral expectations (F_ST) revealed pronounced departures from neutral expectations in all three behavioral phenotypes, consistent with the hypothesis that directional selection has driven the divergence of behavior in dwarf and normal lake whitefish ecotypes.

The relationship between body size and metabolic rate is a crucial issue in organismal biology and evolution. There has been considerable debate over whether the scaling exponent of the relationship is 0.75 (Kleiber's Law) or 0.67. Here we show that determination of this exponent for mammals depends on both the evolutionary tree and the regression model used in the comparative analysis. For example, more recent molecular-based phylogenies tend to support a 0.67 exponent, whereas older phylogenies, mostly based on morphological data, suggest a 0.75 exponent. However, molecular phylogenies yield more variable results than morphological phylogenies and thus are not currently helping to resolve the issue.