Theoretical studies of speciation have been dominated by numerical simulations aiming to demonstrate that speciation in a certain scenario may occur. What is needed now is a shift in focus to identifying more general rules and patterns in the dynamics of speciation. The crucial step in achieving this goal is the development of simple and general dynamical models that can be studied not only numerically but analytically as well. I review some of the existing analytical results on speciation. I first show why the classical theories of speciation by peak shifts across adaptive valleys driven by random genetic drift run into trouble (and into what kind of trouble). Then I describe the Bateson-Dobzhansky-Muller (BDM) model of speciation that does not require overcoming selection. I describe exactly how the probability of speciation, the average waiting time to speciation, and the average duration of speciation depend on the mutation and migration rates, population size, and selection for local adaptation. The BDM model postulates a rather specific genetic architecture of reproductive isolation. I then show exactly why the genetic architecture required by the BDM model should be common in general. Next I consider the multilocus generalizations of the BDM model again concentrating on the qualitative characteristics of speciation such as the average waiting time to speciation and the average duration of speciation. Finally, I consider two models of sympatric speciation in which the conditions for sympatric speciation were found analytically. A number of important conclusions have emerged from analytical studies. Unless the population size is small and the adaptive valley is shallow, the waiting time to a stochastic transition between the adaptive peaks is extremely long. However, if transition does happen, it is very quick. Speciation can occur by mutation and random drift alone with no contribution from selection as different populations accumulate incompatible genes. The importance of mutations and drift in speciation is augmented by the general structure of adaptive landscapes. Speciation can be understood as the divergence along nearly neutral networks and holey adaptive landscapes (driven by mutation, drift, and selection for adaptation to a local biotic and/or abiotic environment) accompanied by the accumulation of reproductive isolation as a by-product. The waiting time to speciation driven by mutation and drift is typically very long. Selection for local adaptation (either acting directly on the loci underlying reproductive isolation via their pleiotropic effects or acting indirectly via establishing a genetic barrier to gene flow) can significantly decrease the waiting time to speciation. In the parapatric case the average actual duration of speciation is much shorter than the average waiting time to speciation. Speciation is expected to be triggered by changes in the environment. Once genetic changes underlying speciation start, they go to completion very rapidly. Sympatric speciation is possible if disruptive selection and/or assortativeness in mating are strong enough. Sympatric speciation is promoted if costs of being choosy are small (or absent) and if linkage between the loci experiencing disruptive selection and those controlling assortative mating is strong.

We present moments and likelihood methods that estimate a DNA substitution rate from a group of closely related sister species pairs separated at an assumed time, and we test these methods with simulations. The methods also estimate ancestral population size and can test whether there is a significant difference among the ancestral population sizes of the sister species pairs. Estimates presented in the literature often ignore the ancestral coalescent prior to speciation and therefore should be biased upward. The simulations show that both methods yield accurate estimates given sample sizes of five or more species pairs and that better likelihood estimates are obtained if there is no significant difference among ancestral population sizes. The model presented here indicates that the larger than expected variation found in multitaxa datasets can be explained by variation in the ancestral coalescence and the Poisson mutation process. In this context, observed variation can often be accounted for by variation in ancestral population sizes rather than invoking variation in other parameters, such as divergence time or mutation rate. The methods are applied to data from two groups of species pairs (sea urchins and Alpheus snapping shrimp) that are thought to have separated by the rise of Panama three million years ago.

Pleurocarpous mosses, characterized by lateral female gametangia and highly branched, interwoven stems, comprise three orders and some 5000 species, or almost half of all moss diversity. Recent phylogenetic analyses resolve the Ptychomniales as sister to the Hypnales plus Hookeriales. Species richness is highly asymmetric with approximately 100 Ptychomniales, 750 Hookeriales, and 4400 Hypnales. Chloroplast DNA (cpDNA) sequences were obtained to compare partitioning of molecular diversity among the orders with estimates of species richness, and to test the hypothesis that either the Hookeriales or Hypnales underwent a period (or periods) of exceptionally rapid diversification. Levels of biodiversity were quantified using explicitly historical “phylogenetic diversity” and nonhistorical estimates of standing sequence diversity. Diversification rates were visualized using lineage-through-time (LTT) plots, and statistical tests of alternative diversification models were performed using the methods of Paradis (1997). The effects of incomplete sampling on the shape of LTT plots and performance of statistical tests were investigated using simulated phylogenies with incomplete sampling. Despite a much larger number of accepted species, the Hypnales contain lower levels of (cpDNA) biodiversity than their sister group, the Hookeriales, based on all molecular measures. Simulations confirm previous results that incomplete sampling yields diversification patterns that appear to reflect a decreasing rate through time, even when the true phylogenies were simulated with constant rates. Comparisons between simulated results and empirical data indicate that a constant rate of diversification cannot be rejected for the Hookeriales. The Hypnales, however, appear to have undergone a period of exceptionally rapid diversification for the earliest 20% of their history.

Theory predicts a significant relationship between the size of a population and the magnitude and composition of its genetic load, but few natural populations have been investigated. We examined the magnitude of genetic load due to recessive deleterious alleles (GL) both segregating and fixed within Gentianella germanica populations of varying size by selfing and reciprocally crossing plants within and between natural populations according to a partial diallel design and by comparing the performance of the experimental progeny in a common-garden experiment. The results show that GL for total fitness in small populations (fewer than 200 plants) was mainly due to fixed recessive deleterious alleles, whereas GL for total fitness in larger populations (more than 200 plants) appeared to be mainly due to segregating deleterious recessive alleles. The total fitness of selfed plants increased with decreasing population size, indicating some purging of deleterious alleles associated with declining population sizes. The magnitudes of GL due to fixed deleterious alleles in small populations and segregating deleterious alleles in large populations, however, were overall similar, suggesting that purging selection was an insignificant force when compared to genetic drift in determining the magnitude of GL in small natural populations in this species. The results of this study highlight the importance of population size in determining the dynamics of genetic loads of natural populations and are overall in line with a large body of theoretical work indicating that small populations may face higher extinction risks due to the fixation and accumulation of deleterious alleles of small effect.

Orchids of the genus Ophrys (Orchidaceae) are pollinated by male bees and wasps through sexual deception. The Ophrys sphegodes group encompasses several closely related species that differ slightly in floral morphology and are pollinated by different solitary bee species. Populations representing different species of the O. sphegodes group often flower simultaneously in sympatry. To test whether gene flow across the species boundaries occurs in these sympatric populations, or whether they are reproductively isolated, we examined the distribution of genetic variation within and among populations and species of this group. We collected at each of five different localities in southern France and Italy two sympatric, co-flowering Ophrys populations, representing six Ophrys species in total. The six microsatellite loci surveyed were highly variable. Genetic differentiation among geographically distant populations of the same species was lower than differentiation among sympatric populations of different species. However, the strength of genetic differentiation among species was among the lowest reported for orchids. Genotype assignment tests and marker-based estimates of gene flow revealed that gene flow across species boundaries occurred and may account for the low observed differentiation among species. These results suggest that sexual deceit pollination in Ophrys may be less specific than thought, or that rare mistakes occur.

Genetic variation in tolerance for seaweed chemical defenses among populations of marine herbivores is poorly described, either because few marine studies have addressed the issue or because the feeding preferences of populations of marine herbivores rarely differentiate. We address this using the polyphagous amphipod Ampithoe longimana (Crustacea), which commonly consumes, and lives on, terpene-rich brown seaweeds in the tropical genus Dictyota. Previous work demonstrated that in areas where the amphipod and seaweeds overlap geographically, these chemically deterrent seaweeds provide the amphipod with spatial refuges from predation by omnivorous fishes that avoid consuming Dictyota. However, the amphipod's distribution extends northward of the seaweed's distribution, making this benefit of association unavailable to more northerly populations of amphipods. On average, populations sympatric with Dictyota have stronger feeding preference for Dictyota species and greater fitness when raised on Dictyota than do populations that are outside Dictyota's geographic endpoint. These results are consistent with the hypothesis that in areas where the amphipods co-exist with Dictyota, selection favors amphipods that tolerate Dictyota's chemical defenses and thereby access its enemy-free space. Amphipods allopatric with Dictyota are unable to tolerate Dictyota's chemical defenses, either because of selection, drift, or some mix of both forces. A trade-off between preference for Dictyota and for the chemically distinct seaweed Hypnea musciformis is suggested by the finding that populations with high preference for Dictyota tended to have lower preference for Hypnea. However, this population-level pattern appears to be the result of independent evolution and not a genetic trade-off: an analysis of full-sib families within a single population detected heritable variation in preferences for Dictyota and suggested this for Hypnea (although the latter was not statistically significant), yet no family-level covariation was detected. Phylogeographic analysis of mitocondrial DNA and nuclear sequences indicates a strong historical break between populations sympatric with Dictyota and populations more than 500 km beyond Dictyota's geographic endpoint. The historical separation of these most northerly populations could have allowed a random accumulation of alleles to effectively degrade tolerance for Dictyota, although selection may also be responsible for the degradation.

While molecular and quantitative trait variation may be theoretically correlated, empirical studies using both approaches frequently reveal discordant patterns, and these discrepancies can contribute to our understanding of evolutionary processes. Here, we assessed genetic variation in six populations of the copepod Tigriopus californicus. Molecular variation was estimated using five polymorphic microsatellite loci, and quantitative variation was measured using 22-life history and morphometric characters. Within populations, no correlation was found between the levels of molecular variation (heterozygosity) and quantitative variation (heritability). Between populations, quantitative subdivision (Q_ST) was correlated with molecular subdivision when measured as F_ST but not when measured as R_ST. Unlike most taxa studied to date, the overall level of molecular subdivision exceeded the level of quantitative subdivision (F_ST = 0.80, R_ST = 0.89, Q_ST = 0.30). Factors that could contribute to this pattern include stabilizing or fluctuating selection on quantitative traits or accelerated rates of molecular evolution.

The genetic incompatibilities that underlie F₂ hybrid breakdown and reproductive isolation between allopatric populations may be susceptible to environmental interactions. Here we show that epistatic interactions between cytochrome c (CYC) alleles and mitochondrial DNA (mtDNA) variation are dramatically influenced by environmental temperature in interpopulation hybrids of the copepod Tigriopus californicus. CYC is a nuclear-encoded gene that functionally interacts with electron transport system (ETS) complexes composed in part of mtDNA-encoded proteins. Previous studies have provided evidence for functional coadaptation between CYC and ETS complex IV (cytochrome c oxidase) and for cytoplasmic effects on the fitness of CYC genotype in copepod hybrids. In this study, selection on CYC genotype is shown to continue into advanced generation hybrids (F₂–F₈) increasing the likelihood that CYC itself is involved in the interaction (and not a linked factor). Relative viabilities varied markedly between copepods raised in two different temperature/light regimes. These results suggest that both intrinsic coadaptation and extrinsic selection will influence the outcome of natural hybridizations between populations. Furthermore, the results indicate that the fitness of particular hybrid genotypes depends on additional non-mtDNA encoded genes that interact with CYC.

Sea urchins are widely used to study both fertilization and development. In this study we combine the two fields to examine the evolution of reproductive isolation in the genus Heliocidaris. Heliocidaris tuberculata develops indirectly via a feeding larva, whereas the only other species in the genus, H. erythrogramma, has evolved direct development through a nonfeeding larva. We estimated the time of divergence between H. erythrogramma and H. tuberculata from mitochondrial DNA divergence, quantified levels of gametic compatibility between the two species in cross-fertilization assays, and examined the mode of evolution of the sperm protein bindin by sequencing multiple alleles of the two species. Bindin is the major component of the sea urchin sperm acrosomal vesicle, and is involved in sperm-egg attachment and fusion. Based on our analyses, we conclude that: the two species of Heliocidaris diverged less than five million years ago, indicating that direct development can evolve rapidly in sea urchins; since their divergence, the two species have become gametically incompatible; Heliocidaris bindin has evolved under positive selection; and this positive selection is concentrated on the branch leading to H. erythrogramma. Three hypotheses can explain the observed pattern of selection on bindin: (1) it is a correlated response to the evolution of direct development in H. erythrogramma; (2) it is the result of an intraspecific process acting in H. erythrogramma but not in H. tuberculata; or (3) it is the product of reinforcement on the species that invests more energy into each egg to avoid hybridization.

Self-compatible hermaphroditic organisms that mix self-fertilization and outcrossing are of great interest for investigating the evolution of mating systems. We investigate the evolution of selfing in Lymnaea truncatula, a self-compatible hermaphroditic freshwater snail. We first analyze the consequences of selfing in terms of genetic variability within and among populations and then investigate how these consequences along with the species ecology (harshness of the habitat and parasitism) might govern the evolution of selfing. Snails from 13 localities (classified as temporary or permanent depending on their water availability) were sampled in western Switzerland and genotyped for seven microsatellite loci. F_IS (estimated on adults) and progeny array analyses (on hatchlings) provided similar selfing rate estimates of 80%. Populations presented a low polymorphism and were highly differentiated (F_ST = 0.58). Although the reproductive assurance hypothesis would predict higher selfing rate in temporary populations, no difference in selfing level was observed between temporary and permanent populations. However, allelic richness and gene diversity declined in temporary habitats, presumably reflecting drift. Infection levels varied but were not simply related to either estimated population selfing rate or to differences in heterozygosity. These findings and the similar selfing rates estimated for hatchlings and adults suggest that within-population inbreeding depression is low in L. truncatula.

The cytochrome c oxidase enzyme (COX) is comprised of 10 nuclear-encoded subunits and three mitochondrial-encoded subunits in close physical association in the inner mitochondrial membrane. COX passes electrons from cytochrome c to molecular oxygen and pumps protons into the inner mitochondrial space for ATP production. Selection on nuclear-mitochondrial interactions within species should lead to coadaptation of the proteins comprising this important enzyme. Under this model, there should be relatively little disruption of COX activity when mitochondrial genomes are crossed among strains within species. A more pronounced disruption of activity is expected when the mitochondrial genome is expressed in the nuclear background of a different species. We test these hypotheses in Drosophila using hybridization and backcrossing among lines of D. simulans and D. mauritiana. Disrupted cytonuclear genotypes were constructed using backcrosses between two lines of D. simulans (siI and siII) that introduced each divergent mitochondrial DNA (mtDNA) into each nuclear background due to maternal inheritance of mtDNA. Similar crosses were used to introduce each D. simulans mtDNA into the D. mauritiana maI nuclear background. Reconstituted cytonuclear control genotypes were constructed by backcrossing the initial F₁ females to males of the maternal genotype. COX enzyme activities were compared among these disrupted and reconstituted backcross genotypes within and between species. The disruption effect on COX activity was restricted to males of interspecific genotypes. These data support the coadaptation hypothesis and are consistent with predictions that the evolution of modifiers of male mitochondrial dysfunction is hindered by the maternal inheritance of mtDNA. New sequence data for nuclear encoded subunits of COX identified amino acids that may play a role in the disruption effect.

The evolution of sexual dimorphism may occur when natural and sexual selection result in different optimum trait values for males and females. Perhaps the most prominent examples of sexual dimorphism occur in sexually selected traits, for which males usually display exaggerated trait levels, while females may show reduced expression of the trait. In some species, females also exhibit secondary sexual traits that may either be a consequence of a correlated response to sexual selection on males or direct sexual selection for female secondary sexual traits. In this experiment, we simultaneously measure the intersex genetic correlations and the relative strength of sexual selection on males and females for a set of cuticular hydrocarbons in Drosophila serrata. There was significant directional sexual selection on both male and female cuticular hydrocarbons: the strength of sexual selection did not differ among the sexes but males and females preferred different cuticular hydrocarbons. In contrast with many previous studies of sexual dimorphism, intersex genetic correlations were low. The evolution of sexual dimorphism in D. serrata appears to have been achieved by sex-limited expression of traits controlled by genes on the X chromosome and is likely to be in its final stages.

Many host-specific parasites are restricted to a limited range of host species by ecological barriers that impede dispersal and successful establishment. In some cases, microevolutionary differentiation is apparent on top of host specificity, as evidenced by significant parasite population genetic structure among host populations. Ecological barriers responsible for specificity and genetic structure can, in principle, reinforce macroevolutionary processes that generate congruent host-parasite phylogenies. However, few studies have explored both the micro- and macroevolutionary ramifications of close association in a single host-parasite system. Here we compare the macroevolutionary histories of two genera of feather lice (Phthiraptera: Ischnocera) that both parasitize New World pigeons and doves (Aves: Columbiformes). Earlier work has shown that dove body lice (genus Physconelloides) are more host specific and have greater population genetic structure than dove wing lice (Columbicola). We reconstructed phylogenies for representatives of the two genera of lice and their hosts, using nuclear and mitochondrial DNA sequences. The phylogenies were well resolved and generally well supported. We compared the phylogenies of body lice and wing lice to the host phylogeny using reconciliation analyses. We found that dove body lice show strong evidence of cospeciation whereas dove wing lice do not. Although the ecology of body and wing lice is very similar, differences in their dispersal ability may underlie these joint differences in host specificity, population genetic structure, and coevolutionary history.

Photoperiodism, the ability to assess the length of day or night, enables a diverse array of plants, birds, mammals, and arthropods to organize their development and reproduction in concert with the changing seasons in temperate climatic zones. For more than 60 years, the mechanism controlling photoperiodic response has been debated. Photoperiodism may be a simple interval timer, that is, an hourglasslike mechanism that literally measures the length of day or night or, alternatively, may be an overt expression of an underlying circadian oscillator. Herein, we test experimentally whether the rhythmic response in Wyeomyia smithii indicates a causal, necessary relationship between circadian rhythmicity and the evolutionary modification of photoperiodic response over the climatic gradient of North America, or may be explained by a simple interval timer. We show that a day-interval timer is sufficient to predict the photoperiodic response of W. smithii over this broad geographic range and conclude that rhythmic responses observed in classical circadian-based experiments alone cannot be used to infer a causal role for circadian rhythmicity in the evolution of photoperiodic time measurement. More importantly, we argue that the pursuit of circadian rhythmicity as the central mechanism that measures the duration of night or day has distracted researchers from consideration of the interval-timing processes that may actually be the target of natural selection linking internal photoperiodic time measurement to the external seasonal environment.

We examined the possible role of hybridization in the invasion process of the African honeybee by testing two hypotheses regarding fluctuating asymmetry (FA), a measure of developmental stability, in wing characteristics: (1) FA should be higher in hybrid versus parental genotypes of African and European races; (2) FA should be lower in African bees compared to hybrid and European workers. Parental and reciprocal hybrid worker genotypes were cross fostered in common-hive rearing environments. We did not find greater FA for wing size and shape in the hybrids compared to both parental types. However, we did find significantly lower FA of shape in the African workers compared to the European and hybrid workers, suggesting that European bees and their hybrids may have compromised fitness relative to African bees. We also found that the two hybrid genotypes significantly differed in overall wing size and shape. If these differences affect wing aerodynamics, then the paternity of hybrids may influence worker performance and could potentially contribute to the loss of European matrilines. Hybridization had few consistent effects on directional asymmetry for wing size and shape. Genotypic factors played a far greater role in determining the effect of hybridization on wing morphology than did differences in rearing environment. Thus, African bees may have lower FA for wing shape (and by inference greater developmental stability) relative to European and hybrid workers, which may contribute to the ability of African bees to displace European honeybee races in invaded regions.

In a colony headed by a single monandrous foundress, theories predict that conflicts between a queen and her workers over both sex ratio and male production should be intense. If production of males by workers is a function of colony size, this should affect sex ratios, but few studies have examined how queens and workers resolve both conflicts simultaneously. We conducted field and laboratory studies to test whether sex-ratio variation can be explained by conflict over male production between queen and workers in the primitively eusocial wasp Polistes chinensis antennalis.

Worker oviposition rate increased more rapidly with colony size than did queen oviposition. Allozyme and microsatellite markers revealed that the mean frequency of workers' sons among male adults in queen-right colonies was 0.39 ± 0.08 SE (n = 22). Genetic relatedness among female nestmates was high (0.654–0.796), showing that colonies usually had a single, monandrous queen. The mean sex allocation ratio (male investment/male and gyne investments) of 46 queen-right colonies was 0.47 ± 0.02, and for 25 orphaned colonies was 0.86 ± 0.04. The observed sex allocation ratio was likely to be under queen control. For queen-right colonies, the larger colonies invested more in males and produced reproductives protandrously and/or simultaneously, whereas the smaller colonies invested more in females and produced reproductives protogynously. Instead of positive relationships between colony size and worker oviposition rate, the frequency of workers' sons within queen-right colonies did not increase with colony size. These results suggest that queens control colony investment, even though they allow worker oviposition in queen-right colonies. Eggs laid by workers may be policed by the queen and/or fellow workers. Worker oviposition did not influence the outcome of sex allocation ratio as a straightforward function of colony size.

Using collections from the years 1892–1999, I determined maximum standard length within each of 1030 populations of riffle-inhabiting darters (Etheostoma spp.) representing five species from 788 sites in Illinois. Each site contained one to four riffle-inhabiting species of Etheostoma. Based on maximum-sized individuals in each collection, I calculated a ratio of standard lengths for all sympatric species pairs. Null models were developed using random pairings of body size measurements from sites with only one species of riffle Etheostoma to test whether body size ratios of sympatric species pairs are larger, indicating divergence, or smaller, indicating convergence, than ratios generated from allopatric populations. Results suggest that two of nine species pairs tend to converge in body size when sympatric. This suggests that convergence, in some cases, may facilitate the persistence of similar species in diverse communities. When variation in sympatric congener number is accounted for, the E. caeruleum:E. spectabile interaction results in significantly increasing size ratios as congener number increases. Etheostoma caeruleum and E. spectabile are the most ecologically similar and evolutionarily closely related species in the dataset. The divergence within this pair, associated with sympatric congener number, suggests that the degree of divergent character displacement may be indirectly enhanced by community structure.

Anolis lizards in the Greater Antilles partition the structural microhabitats available at a given site into four to six distinct categories. Most microhabitat specialists, or ecomorphs, have evolved only once on each island, yet closely related species of the same ecomorph occur in different geographic macrohabitats across the island. The extent to which closely related species of the same ecomorph have diverged to adapt to different geographic macrohabitats is largely undocumented. On the island of Hispaniola, members of the Anolis cybotes species group belong to the trunk-ground ecomorph category. Despite evolutionary stability of their trunk-ground microhabitat, populations of the A. cybotes group have undergone an evolutionary radiation associated with geographically distinct macrohabitats. A combined phylogeographic and morphometric study of this group reveals a strong association between macrohabitat type and morphology independent of phylogeny. This association results from long-term morphological evolutionary stasis in populations associated with mesic-forest environments (A. c. cybotes and A. marcanoi) and predictable morphometric changes associated with entry into new macrohabitat types (i.e., xeric forests, high-altitude pine forest, rock outcrops). Phylogeographic analysis of 73 new mitochondrial DNA sequences (1921 aligned sites) sampled from 68 geographic populations representing 12 recognized species and subspecies diagnoses 16 allopatric or parapatric groupings of populations differing from each other by 5–18% sequence divergence. At least some of these groupings appear to have attained species-level divergence from others. Evolutionary specialization to different macrohabitat types may be a major factor in the evolutionary diversification of Greater Antillean anoles.

If DNA replication is a major cause of mutation, then those life-history characters, which are expected to affect the number of male germline cell divisions, should also affect the male to female mutation bias (α_m). We tested this hypothesis by comparing several clades of bird species, which show variation both in suitable life-history characters (generation time as measured by age at first breeding and sexual selection as measured by frequency of extrapair paternity) and in α_m, which was estimated by comparing Z-linked and W-linked substitution rates in gametologous introns. α_m differences between clades were found to positively covary with both generation time and sexual selection, as expected if DNA replication causes mutation. The effects of extrapair paternity frequency on α_m suggests that increased levels of sexual selection cause higher mutation rates, which offers an interesting solution to the paradox of the loss of genetic variance associated with strong directional sexual selection. We also used relative rate tests to examine whether the observed differences in α_m between clades were due to differences in W-linked or Z-linked substitution rates. In one case, a significant difference in α_m between two clades was shown to be due to W-linked rates and not Z-linked rates, a result that suggests that mutation rates are not determined by replication alone.

Hybridizing California and Gambel's Quail (Callipepla californica and C. gambelii) are unlike many hybridizing avian species in that pairing primarily occurs within a flock, or covey, that is composed of several families. Coveys in the area of sympatry contain mostly hybrid individuals, relative to parental types. I tested whether individuals perceive covey mates as a single species and whether pairing within the covey causes inbreeding and a loss of reproductive success. Individuals discriminated between the parental species in captivity, but actual pairing in the wild was random with respect to species. Contrary to expectation, coveys were not more inbred than the local population. Results suggest that potential costs of inbreeding may be avoided through sex-biased dispersal and nonrandom pairing. Furthermore, breeding occurred earlier and with greater success in pairs formed within the resident covey, rather than outside it. These findings suggest that fitness benefits gained from pairing early within a mixed-species covey promote interspecific pairing. This study provides evidence that mating behaviors specific to local conditions maintain hybrid zones when genetic costs to interbreeding are small.

Despite their importance in evolutionary biology, heritability and the strength of natural selection have rarely been estimated in wild populations of iteroparous species or have usually been limited to one particular event during an organism's lifetime. Using an animal-model restricted maximum likelihood and phenotypic selection models, we estimated quantitative genetic parameters and the strength of lifetime selection on parturition date and litter size at birth in a natural population of North American red squirrels, Tamiasciurus hudsonicus. Litter size at birth and parturition date had low heritabilities (h² = 0.15 and 0.16, respectively). We considered potential effects of temporal environmental covariances between phenotypes and fitness and of spatial environmental heterogeneity in estimates of selection. Selection favored early breeders and females that produced litter sizes close to the population average. Stabilizing selection on litter size at birth may occur because of a trade-off between number of offspring produced per litter and offspring survival or a trade-off between a female's fecundity and her future reproductive success and survival.

Due to a trade-off between reproduction and life span, highly fertile individuals often live shorter lives than nonreproductive conspecifics. Perennial eusocial insects are exceptional in that reproductive queens live considerably longer than the nonreproductive workers. The two female castes may differ strongly in morphology, ontogeny, physiology, diet, behavior, and mating, and all these differences could be responsible for life span differences. In the ponerine ant Platythyrea punctata, morphological and ontogenetic caste differences do not exist. Instead, all workers are capable of producing diploid offspring through thelytokous parthenogenesis, and colonies are essentially clones. Here, we show that reproductives live significantly longer than nonreproductive workers. Reproductives stay in the nest during their whole life, whereas nonreproductives switch from intranidal tasks to foraging when they get older. Different work load and different hormone titers might proximately underlie the different life span of reproductives and nonreproductives in this ant.