We analyze the changes in the mean and variance components of a quantitative trait caused by changes in allele frequencies, concentrating on the effects of genetic drift. We use a general representation of epistasis and dominance that allows an arbitrary relation between genotype and phenotype for any number of diallelic loci. We assume initial and final Hardy-Weinberg and linkage equilibrium in our analyses of drift-induced changes. Random drift generates transient linkage disequilibria that cause correlations between allele frequency fluctuations at different loci. However, we show that these have negligible effects, at least for interactions among small numbers of loci. Our analyses are based on diffusion approximations that summarize the effects of drift in terms of F, the inbreeding coefficient, interpreted as the expected proportional decrease in heterozygosity at each locus. For haploids, the variance of the trait mean after a population bottleneck is var(Δz̄) = Σⁿ_k=1 F^kV_A(k), where n is the number of loci contributing to the trait variance, V_A(1) = V_A is the additive genetic variance, and V_A(k) is the kth-order additive epistatic variance. The expected additive genetic variance after the bottleneck, denoted ⟨V^*_A⟩, is closely related to var(Δz̄); ⟨V^*_A⟩ = (1 − F) Σⁿ_k=1 kF^k−1V_A(k). Thus, epistasis inflates the expected additive variance above V_A(1 − F), the expectation under additivity. For haploids (and diploids without dominance), the expected value of every variance component is inflated by the existence of higher order interactions (e.g., third-order epistasis inflates ⟨V^*_AA⟩). This is not true in general with diploidy, because dominance alone can reduce ⟨V^*_A⟩ below V_A(1 − F) (e.g., when dominant alleles are rare). Without dominance, diploidy produces simple expressions: var(Δz̄) = Σⁿ_k=1 (2F)^kV_A(k) and ⟨V^*_A⟩ = (1 − F) Σⁿ_k=1 k(2F)^k−1V_A(k). With dominance (and even without epistasis), var(Δz̄) and ⟨V^*_A⟩ no longer depend solely on the variance components in the base population. For small F, the expected additive variance simplifies to ⟨V^*_A⟩ ≃ (1 − F)V_A 4FV_AA 2FV_D 2FC_AD, where C_AD is a sum of two terms describing covariances between additive effects and dominance and additive × dominance interactions. Whether population bottlenecks lead to expected increases in additive variance depends primarily on the ratio of nonadditive to additive genetic variance in the base population, but dominance precludes simple predictions based solely on variance components. We illustrate these results using a model in which genotypic values are drawn at random

The fundamental equation in evolutionary quantitative genetics, the Lande equation, describes the response to directional selection as a product of the additive genetic variance and the selection gradient of trait value on relative fitness. Comparisons of both genetic variances and selection gradients across traits or populations require standardization, as both are scale dependent. The Lande equation can be standardized in two ways. Standardizing by the variance of the selected trait yields the response in units of standard deviation as the product of the heritability and the variance-standardized selection gradient. This standardization conflates selection and variation because the phenotypic variance is a function of the genetic variance. Alternatively, one can standardize the Lande equation using the trait mean, yielding the proportional response to selection as the product of the squared coefficient of additive genetic variance and the mean-standardized selection gradient. Mean-standardized selection gradients are particularly useful for summarizing the strength of selection because the mean-standardized gradient for fitness itself is one, a convenient benchmark for strong selection. We review published estimates of directional selection in natural populations using mean-standardized selection gradients. Only 38 published studies provided all the necessary information for calculation of mean-standardized gradients. The median absolute value of multivariate mean-standardized gradients shows that selection is on average 54% as strong as selection on fitness. Correcting for the upward bias introduced by taking absolute values lowers the median to 31%, still very strong selection. Such large estimates clearly cannot be representative of selection on all traits. Some possible sources of overestimation of the strength of selection include confounding environmental and genotypic effects on fitness, the use of fitness components as proxies for fitness, and biases in publication or choice of traits to study.

We examined barriers to gene flow in a hybrid zone of two subspecies of the song sparrow (Melospiza melodia). We focused on how mating signals and mate choice changed along an environmental gradient and gathered data on the morphology, genetics, ecology, and behavior across the zone. Melospiza m. heermanni of the Pacific slope of California and M. m. fallax of the Sonoran Desert, each distinct in plumage, meet across a steep environmental gradient in southeastern California. Although both subspecies occur in riparian habitat, their occupied habitat differs structurally, the former subspecies occurring in areas with denser understory and greater vertical heterogeneity. Song elements varied concomitantly, as predicted by the acoustic adaptation hypothesis, with heermanni having lower-pitched, more widely spaced elements. Females of both subspecies responded more strongly to homotypic than heterotypic song, and addition of subspecific plumage cues increased response if song was homotypic but not if heterotypic. Females thus assess multiple male traits, weighing song more heavily. Males of both subspecies showed significantly greater agonistic response to homotypic song. Microsatellite variation is correlated significantly with plumage variation across the zone and suggests limited gene flow between the taxa. The association of song and plumage with the environment and in turn with assortative mating suggests a means by which reproductive isolation may evolve or be maintained in hybrid zones.

The success of obligate endosymbiotic Wolbachia infections in insects is due in part to cytoplasmic incompatibility (CI), whereby Wolbachia bacteria manipulate host reproduction to promote their invasion and persistence within insect populations. The observed diversity of CI types raises the question of what the evolutionary pathways are by which a new CI type can evolve from an ancestral type. Prior evolutionary models assume that Wolbachia exists within a host individual as a clonal infection. While endosymbiotic theory predicts a general trend toward clonality, Wolbachia provides an exception in which there is selection to maintain diversity. Here, evolutionary trajectories are discussed that assume that a novel Wolbachia variant will co-exist with the original infection type within a host individual as a superinfection. Relative to prior models, this assumption relaxes requirements and allows additional pathways for the evolution of novel CI types. In addition to describing changes in the Wolbachia infection frequency associated with the hypothesized evolutionary events, the predicted impact of novel CI variants on the host population is also described. This impact, resulting from discordant evolutionary interests of symbiont and host, is discussed as a possible cause of Wolbachia loss from the host population or host population extinction. The latter is also discussed as the basis for an applied strategy for the suppression of insect pest populations. Model predictions are discussed relative to a recently published Wolbachia genome sequence and prior characterization of CI in naturally and artificially infected insects.

Cytoplasmically inherited symbiotic Wolbachia bacteria are known to induce a diversity of phenotypes on their numerous arthropod hosts including cytoplasmic incompatibility, male-killing, thelytokous parthenogenesis, and feminization. In the wasp Asobara tabida (Braconidae), in which all individuals harbor three genotypic Wolbachia strains (wAtab1, wAtab2 and wAtab3), the presence of Wolbachia is required for insect oogenesis. To elucidate the phenotype of each Wolbachia strain on host reproduction, especially on oogenesis, we established lines of A. tabida harboring different combinations of these three bacterial strains. We found that wAtab3 is essential for wasp oogenesis, whereas the two other strains, wAtab1 and wAtab2, seem incapable to act on this function. Furthermore, interline crosses showed that strains wAtab1 and wAtab2 induce partial (about 78%) cytoplasmic incompatibility of the female mortality type. These results support the idea that bacterial genotype is a major factor determining the phenotype induced by Wolbachia on A. tabida hosts. We discuss the implications of these findings for current hypotheses regarding the evolutionary mechanisms by which females of A. tabida have become dependent on Wolbachia for oogenesis.

The eastern Asian (EAS)–eastern North American (ENA) floristic disjunction is one of the best-known biogeographic patterns in the Northern Hemisphere. Recent paleontological and molecular analyses have illuminated the origins of the biogeographic pattern, but subsequent diversification and evolution of the disjunct floras in each of the two continents after isolation remains poorly understood. Although similar in climate and floristic composition, EAS has twice as many species as ENA in genera occurring in both regions. Explaining such differences in species diversity between regions with similar environmental conditions (diversity anomalies) is an important goal of the study of the global patterns of biodiversity. We used a phylogenetic approach to compare rates of net speciation and molecular evolution between the two regions. We first identified EAS–ENA disjunct sister clades from ten genera (Asarum, Buckleya, Carpinus, Carya, Cornus, Hamamelis, Illicium, Panax, Stewartia, and Styrax) that represent diverse angiosperm lineages using phylogenetic analyses of ITS (internal transcribed spacer of nuclear ribosomal DNA) sequence data. Species richness and substitution rate of ITS between sister clades were compared. The results revealed a pattern of greater species diversity in the EAS counterparts. A positive relationship between species diversity and ITS substitution rate was also documented. These results suggest greater net speciation and accelerated molecular evolution in EAS. The data support the idea that a regional difference in net speciation rate related to topographic heterogeneity contributes to the diversity anomaly between EAS and ENA. The close relationship between rates of ITS evolution and species richness further suggests that species production may be directly linked to rate of nucleotide substitution.

The biota of Hawaiian Islands is derived entirely from long distance dispersal, often followed by in situ speciation. Species descended from each colonist constitute monophyletic lineages that have diverged to varying degrees under similar spatial and temporal constraints. We partitioned the Hawaiian angiosperm flora into lineages and assessed morphological, ecological, and biogeographic characteristics to examine their relationships to variation in species number (S). Lineages with external bird dispersal (through adhesion) were significantly more species-rich than those with abiotic dispersal, but only weakly more species-rich than lineages with internal bird dispersal (involving fleshy fruits). Pollination mode and growth form (woody vs. herbaceous) had no significant effect on S, in contrast to studies of angiosperm families. S relates positively to the geographic and ecological range size of whole lineages, but negatively to local abundance and mean range sizes of constituent species. Species-rich lineages represent a large proportion of major adaptive shifts, although this appears to be an artifact of having more species. Examination of 52 sister species pairs in numerous lineages provides evidence for allopatric (including peripheral isolates) and parapatric (ecological) modes, with 15 cases of each. Although postspeciational dispersal may obscure these modes in many of the remaining cases, instances of sympatric and hybrid speciation are also discussed. Because speciation is both a consequence and a cause of ecological and biogeographic traits, speciation mode may be integral to relationships between traits. We discuss the role of speciation in shaping the regional species pool.

Species-specific obligate pollination mutualism between Glochidion trees (Euphorbiaceae) and Epicephala moths (Gracillariidae) involves a large number of interacting species and resembles the classically known fig–fig wasp and yucca–yucca moth associations. To assess the extent of parallel cladogenesis in Glochidion-Epicephala association, we reconstruct phylogenetic relationships of 18 species of Glochidion using nuclear ribosomal DNA sequences (internal and external transcribed spacers) and those of the corresponding 18 Epicephala species using mitochondrial (the cytochrome oxidase subunit I gene) and nuclear DNA sequences (the arginine kinase and elongation factor-1α genes). Based on the obtained phylogenies, we determine whether Glochidion and Epicephala have undergone parallel diversification using several different methods for investigating the level of cospeciation between phylogenies. These tests indicate that there is generally a greater degree of correlation between Glochidion and Epicephala phylogenies than expected in a random association, but the results are sensitive to selection of different phylogenetic hypotheses and analytical methods for evaluating cospeciation. Perfect congruence between phylogenies is not found in this association, which likely resulted from host shift by the moths. The observed significant discrepancy between Glochidion and Epicephala phylogenies implies that the one-to-one specificity between the plants and moths has been maintained through a complex speciation process or that there is an underestimated diversity of association between Glochidion trees and Epicephala moths.

This study aims at a better understanding of the evolutionary significance of viviparity in some freshwater gastropods. We use a phylogeny based on partial sequences of the mitochondrial 16S gene of representatives of the limnetic and pantropical Pachychilidae to infer the relationships within this particular group of cerithioideans and the evolution of reproductive strategies. The phylogeny presented herein implies a new systematization and suggests that viviparity has appeared three times among the Pachychilidae. This is supported by the finding of very distinct reproductive morphologies in different lineages of viviparous taxa that are exclusively found in Southeast Asia. Based on the observation that oviparity is the ancestral character state in this freshwater family, we conclude that viviparity has evolved subsequent to the exploration of freshwater. We present data showing that all Pachychilidae produce considerably larger but fewer egg capsules compared to most marine snails. In other studies on freshwater gastropods, this has been discussed as an adaptation to freshwater environments. In this context we hypothesize that the increased parental investment involved in the enlargement of eggs in concert with the reduction of clutch sizes was the driving factor that ultimately lead to the evolution of viviparity in the Asian taxa. Consequently, although not directly correlated with the colonization of the new adaptive zone, viviparity is strongly favored by other consequences of this step. Hence, we hypothesize that the production of large eggs, which is necessitated by the exploration of freshwater, represents a preadaptation existing in those ancestors from which viviparous pachychilid lineages eventually evolved in Southeast Asia.

A phylogenetic approach to the origin and maintenance of species diversity ideally requires the sampling of all species within a clade, confirmation that they are evolutionarily distinct entities, and knowledge of their geographical distributions. In the marine tropics such studies have mostly been of fish and reef-associated organisms, usually with high dispersal. In contrast, snails of the genus Echinolittorina (Littorinidae) are restricted to rocky shores, have a four-week pelagic development (and recorded dispersal up to 1400 km), and show different evolutionary patterns. We present a complete molecular phylogeny of Echinolittorina, derived from Bayesian analysis of sequences from nuclear 28S rRNA and mitochondrial 12S rRNA and COI genes (nodal support indicated by posterior probabilities, maximum likelihood, and neighbor-joining bootstrap). This consists of 59 evolutionarily significant units (ESUs), including all 50 known taxonomic species. The 26 ESUs found in the Indo-West Pacific region form a single clade, whereas the eastern Pacific and Atlantic species are basal. The earliest fossil occurred in the Tethys during the middle Eocene and we suggest that the Indo-West Pacific clade has been isolated since closure of the Tethyan seaway in the early Miocene. The geographical distributions of all species (based on more than 3700 locality records) appear to be circumscribed by barriers of low temperature, unsuitable sedimentary habitat, stretches of open water exceeding about 1400 km, and differences in oceanographic conditions on the continuum between oceanic and continental. The geographical ranges of sister species show little or no overlap, indicating that the speciation mode is predominantly allopatric. Furthermore, range expansion following speciation appears to have been limited, because a high degree of allopatry is maintained through three to five branching points of the phylogeny. This may be explained by infrequent long-distance colonization, habitat specialization on the oceanic/continental gradient, and perhaps by interspecific competition. In the eastern Pacific plus Atlantic we identify five cases of divergence on either side of the Isthmus of Panama, but our estimates of their ages pre-date the emergence of the Isthmus. There are three examples of sister relationships between species in the western Atlantic and eastern Atlantic, all resulting from dispersal to the east. Within the Indo-West Pacific, we find no geographical pattern of speciation events; narrowly endemic species of recent origin are present in both peripheral and central parts of the region. Evidence from estimated divergence times of sister species, and from a plot of the number of lineages over time, suggest that there has been no acceleration of diversification during the glacio-eustatic cycles of the Plio-Pleistocene. In comparison with reefal organisms, species of Echinolittorina on rocky shores may be less susceptible to extinction or isolation during sea-level fluctuations. The species richness of Echinolittorina in the classical biogeographic provinces conforms to the common pattern of highest diversity (11 species) in the central “East Indies Triangle” of the Indo-West Pacific, with a subsidiary focus in the eastern Pacific and western Atlantic, and lowest diversity in the eastern Atlantic. The diversity focus in the East Indies Triangle is produced by a mosaic of restricted allopatric species and overlap of a few widespread ones, and is the result of habitat specialization rather than historical vicariance. This study emphasizes the plurality of biogeographic histories and speciation patterns in the marine tropics.

Almost all of the more than 200 species of fungus-growing ants (Formicidae: Attini) cultivate litter-decomposing fungi in the family Lepiotaceae (Basidiomycota: Agaricales). The single exception to this rule is a subgroup of ant species within the lower attine genus Apterostigma, which cultivate pterulaceous fungi distantly related to the Lepiotaceae. Comparison of cultivar and ant phylogenies suggests that a switch from lepiotaceous to pterulaceous fungiculture occurred only once in the history of the fungus-growing ants. This unique switch occurred after the origin of the genus Apterostigma, such that the basal Apterostigma lineages retained the ancestral attine condition of lepiotaceous fungiculture, and none of the Apterostigma lineages in the monophyletic group of pterulaceous fungiculturists are known to have reverted back to lepiotaceous fungiculture. The origin of pterulaceous fungiculture in attine ants may have involved a unique transition from the ancestral cultivation of litter-decomposing lepiotaceous fungi to the cultivation of wood-decomposing pterulaceous fungi. Phylogenetic analyses further indicate that distantly related Apterostigma ant species sometimes cultivate the same cultivar lineage, indicating evolutionarily frequent, and possibly ongoing, exchanges of fungal cultivars between Apterostigma ant species. The pterulaceous cultivars form two sister clades, and different Apterostigma ant lineages are invariably associated with, and thus specialized on, only one of the two cultivar clades. However, within clades Apterostigma ant species are able to switch between fungi. This pattern of broad specialization by attine ants on defined cultivar clades, coupled with flexible switching between fungi within cultivar clades, is also found in other attine lineages and appears to be a general phenomenon of fungicultural evolution in all fungus-growing ants.

Mutualistic interactions can be exploited by cheaters that take the rewards offered by mutualists without providing services in return. The evolution of cheater species from mutualist ancestors is thought to be possible under particular ecological conditions. Here we provide a test of the first explicit model of the transition from mutualism to antagonism. We used the obligate pollination mutualism between yuccas and yucca moths to examine the origins of a nonpollinating cheater moth, Tegeticula intermedia, and its pollinating sister species, T. cassandra. Based on geographic distribution and ecological factors affecting the pollinators, previous research had indicated that the cheaters evolved in Florida as a result of sympatry of T. cassandra and another pollinator species. We used mitochondrial DNA (mtDNA) sequences and amplified fragment length polymorphism (AFLP) data to investigate the phylogeographic history of the pollinator-cheater sister pair and to test whether the cheaters arose in Florida. Contrary to predictions, phylogenetic and population genetic analyses suggested that the cheaters evolved in the western United States and subsequently spread eastward. Western populations of cheaters had the most ancestral haplotypes and the highest genetic diversity, and there was also significant genetic structure associated with a geographic split between eastern and western populations. In comparison, there was evidence for weak genetic structure between northern and southern pollinator populations, suggesting a long history in Florida. The western origin of the cheaters indicated that the pollinators have more recently become restricted to the southeastern United States. This was supported by AFLP analyses that indicated that the pollinators were more closely related to the western cheaters than they were to geographically proximate cheaters in the east. Shared mtDNA between pollinators and eastern cheaters suggested hybridization, possibly in a secondary contact zone. The results negate the out-of-Florida hypothesis and reveal instead a long, complex, and disparate history for the pollinator-cheater sister pair.

Immune system activation may benefit hosts by generating resistance to parasites. However, natural resources are usually limited, causing a trade-off between the investment in immunity and that in other life-history or sexually selected traits. Despite its importance for the evolution of host defense, state-dependent fitness costs of immunity received little attention under natural conditions. In a field experiment we manipulated the nutritional condition of male field crickets Gryllus campestris and subsequently investigated the effect of an induced immune response through inoculation of bacterial lipopolysaccharides. Immune system activation caused a condition-dependent reduction in body condition, which was proportional to the condition-gain during the preceding food-supplementation period. Independent of nutritional condition, the immune insult induced an enduring reduction in daily calling rate, whereas control-injected males fully regained their baseline level of sexual signaling following a temporary decline. Since daily calling rate affects female mate choice under natural conditions, this suggests a decline in male mating success as a cost of induced immunity. Food supplementation enhanced male life span, whereas the immune insult reduced longevity, independent of nutritional status. Thus, immune system activation ultimately curtails male fitness due to a combined decline in sexual display and life span. Our field study thus indicates a key role for fitness costs of induced immunity in the evolution of host defense. In particular, costs expressed in sexually selected traits might warrant the honest advertisement of male health status, thus representing an important mechanism in parasite-mediated sexual selection.

Quantitative genetics has been introduced to evolutionary biologists with the suggestion that microevolution could be directly linked to macroevolutionary patterns using, among other parameters, the additive genetic variance/ covariance matrix (G) which is a statistical representation of genetic constraints to evolution. However, little is known concerning the rate and pattern of evolution of G in nature, and it is uncertain whether the constraining effect of G is important over evolutionary time scales. To address these issues, seven species of field crickets from the genera Gryllus and Teleogryllus were reared in the laboratory, and quantitative genetic parameters for morphological traits were estimated from each of them using a nested full-sibling family design. We used three statistical approaches (T method, Flury hierarchy, and Mantel test) to compare G matrices or genetic correlation matrices in a phylogenetic framework. Results showed that G matrices were generally similar across species, with occasional differences between some species. We suggest that G has evolved at a low rate, a conclusion strengthened by the consideration that part of the observed across-species variation in G can be explained by the effect of a genotype by environment interaction. The observed pattern of G matrix variation between species could not be predicted by either morphological trait values or phylogeny. The constraint hypothesis was tested by comparing the multivariate orientation of the reconstructed ancestral G matrix to the orientation of the across-species divergence matrix (D matrix, based on mean trait values). The D matrix mainly revealed divergence in size and, to a much smaller extent, in a shape component related to the ovipositor length. This pattern of species divergence was found to be predictable from the ancestral G matrix in agreement with the expectation of the constraint hypothesis. Overall, these results suggest that the G matrix seems to have an influence on species divergence, and that macroevolution can be predicted, at least qualitatively, from quantitative genetic theory. Alternative explanations are discussed.

Predation is heterogeneously distributed across space and time, and is presumed to represent a major source of evolutionary diversification. In fishes, fast-starts—sudden, high-energy swimming bursts—are often important in avoiding capture during a predator strike. Thus, in the presence of predators, we might expect evolution of morphological features that facilitate increased fast-start speed. We tested this hypothesis using populations of western mosquitofish (Gambusia affinis) that differed in level of predation by piscivorous fish. Body morphology of G. affinis males, females, and juveniles diverged in a consistent manner between predatory environments. Fish collected from predator populations exhibited a larger caudal region, smaller head, more elongate body, and a posterior, ventral position of the eye relative to fish from predator-free populations. Divergence in body shape largely matched a priori predictions based on biomechanical principles, and was evident across space (multiple populations) and time (multiple years). We measured maximum burst-swimming speed for male mosquitofish and found that individuals from predator populations produced faster bursts than fish from predator-free populations (about 20% faster). Biomechanical models of fish swimming and intrapopulation morphology-speed correlations suggested that body shape differences were largely responsible for enhanced locomotor performance in fish from predator populations. Morphological differences also persisted in offspring raised in a common laboratory environment, suggesting a heritable component to the observed morphological divergence. Taken together, these results strongly support the hypothesis that divergent selection between predator regimes has produced the observed phenotypic differences among populations of G. affinis. Based on biomechanical principles and recent findings in other species, it appears that the general ecomorphological model described in this paper will apply for many aquatic taxa, and provide insight into the role of predators in shaping the body form of prey organisms.

How much of the variation in adaptive divergence can be explained by gene flow? The answer to this question should objectively reveal whether gene flow generally places a substantial constraint on evolutionary diversification. We studied multiple independent lake-stream population pairs of threespine stickleback (Gasterosteus aculeatus). For each pair, we quantified adaptive divergence based on morphological traits that have a genetic basis and are subject to divergent selection. We then estimated gene flow based on variation at five unlinked microsatellite loci. We found a consistent and significant pattern for morphological divergence to be positively correlated with genetic divergence and negatively correlated with gene flow. Statistical significance and the amount of variation explained varied within and among traits: 36.1–74.1% for body depth and 11.8–51.7% for gill raker number. Variation within each trait was the result of differences among methods for estimating genetic divergence and gene flow. Variation among traits likely reflects different strengths of divergent selection. We conclude that gene flow has a substantial effect on adaptive divergence in nature but that the magnitude of this effect varies among traits. An alternative explanation is that cause and effect are reversed: adaptive divergence is instead constraining gene flow. This effect seems relatively unimportant for our system because genetic divergence and gene flow were not correlated with ecologically relevant habitat features of lakes (surface area) or streams (width, depth, flow, canopy openness).

Because sharks possess an unusual suite of reproductive characteristics, including internal fertilization, sperm storage, relatively low fecundity, and reproductive modes that range from oviparity to viviparity, they can provide important insight into the evolution of mating systems and sexual selection. Yet, to date, few studies have characterized behavioral and genetic mating systems in natural populations of sharks or other elasmobranchs. In this study, highly polymorphic microsatellite loci were used to examine breeding biology of a large coastal shark, the lemon shark, Negaprion brevirostris, at a tropical lagoon nursery. Over six years, 910 lemon sharks were sampled and genotyped. Young were assigned into sibling groups that were then used to reconstruct genotypes of unsampled adults. We assigned 707 of 735 young sharks to one of 45 female genotypes (96.2%), and 485 (66.0%) were assigned to a male genotype. Adult female sharks consistently returned to Bimini on a biennial cycle to give birth. Over 86% of litters had multiple sires. Such high levels of polyandry raise the possibility that polyandry evolved in viviparous sharks to reduce genetic incompatibilities between mother and embryos. We did not find a relationship between relatedness of mates and the number of offspring produced, indicating that inbreeding avoidance was probably not driving pre- or postcopulatory mate choice. Adult male sharks rarely sired more than one litter at Bimini and may mate over a broader geographic area.

The relative importance of natural selection and genetic drift in determining patterns of phenotypic diversity observed in nature is still unclear. The natterjack toad (Bufo calamita) is one of a few amphibian species capable of breeding in saline ponds, even though water salinity represents a considerable stress for them. Results from two common-garden experiments showed a pattern of geographic variation in embryonic salinity tolerance among populations from either fresh or brackish environments, consistent with the hypothesis of local adaptation. Full-sib analysis showed increased variation in survival among sibships within population for all populations as osmotic stress was increased (broad-sense heritability increased as salinity raised). Nevertheless, toads native to the brackish water environment had the highest overall survival under brackish conditions. Levels of population genetic differentiation for salinity tolerance were higher than those of neutral genetic differentiation, the latter obtained through the analysis of eight microsatellite loci. Microsatellite markers also revealed little population differentiation, lack of an isolation-by-distance pattern, and moderate gene flow connecting the populations. Therefore, environmental stress tolerance appears to have evolved in absence of geographic isolation, and consequently we reject the null hypothesis of neutral differentiation.

Throughout their evolutionary histories, marsupial mammals have been taxonomically and morphologically less diverse than their sister taxa the placentals. Because of this, it has been proposed that the evolution of marsupials has been constrained by the functional requirements of their mode of reproduction. Marsupials give birth after short gestation times to immature neonates that immediately crawl, under the power of their precociously developed shoulder girdles, to the teat where they attach and complete their early development. Using a novel approach incorporating adult and embryological morphological data, this study is the first to both: (1) statistically support adult patterns of morphological divergence consistent with the constraint hypothesis, and (2) identify ontogenetic patterns of morphological change that demonstrate that the constraint was responsible, at least in part, for their formation. As predicted by the marsupial constraint, the shoulder girdles of adult marsupials are less diverse than those of adult placentals, and adult marsupial scapulae are less morphologically diverse than adult marsupial pelves. Furthermore, marsupials that complete an extensive crawl to the teat are restricted to a common pattern of ontogenetic scapular shape change, strongly supporting the hypothesis that the morphological development of the marsupial scapula has been limited evolutionarily by its obligate role in the crawl to the teat. Because this study establishes that ontogenetic and evolutionary morphological change is correlated within mammalian scapulae, it is probable that the marsupial constraint also restricted the morphological divergence of the scapula over evolutionary time by limiting ontogenetic change in the scapula. These findings, coupled with the importance of the shoulder girdle in mammalian locomotor specialization, support the conclusion that the low morphological diversity of marsupial forms over evolutionary time could be directly due to the constraint on marsupial morphological evolution caused by the functional requirements of the crawl to the teat.

The wider choice hypothesis suggests that hermaphroditic plants increase their female fitness by producing excess flowers, exposing more ovaries to the scrutiny of selective abortion, and thereby increasing the offspring quality obtained from the best, which escape abortion. Selective discrimination of fruit maturation has been well documented in many angiosperm species, but there has been no examination of how offspring quality varies in relation to the number of excess flowers. This relation must be positive over the entire range of excess flower number for the wider choice mechanism to be the selective force behind the evolution of very large excess floral displays. I examined offspring quality in relation to natural variation in nonfruiting flowers in 72 plants of Pultenaea gunnii, an Australian bush pea. Maximum seed mass and maximum seedling height at 30 days growth showed a nonlinear, saturating relation to excess flower number that is consistent with a theoretical model of wider choice. But the strongly diminishing marginal benefits of excess flowers in this relation make it unlikely that wider choice plays a major role in the evolution of floral display size in this species.

The immune system of invertebrates can mount different responses, including melanotic encapsulation and several antibacterial defense mechanisms. Variation of the efficacies of these responses is generally considered to be a product of the evolutionary pressure on each response due to infection by parasites. However, potential interactions and trade-offs among the different responses of the immune system could constrain the evolutionary potential of each response. In a natural population of the mosquito Anopheles gambiae, we measured the genetic association between the melanization response and an antibacterial response in two environmental qualities (well-fed and undernourished larvae). In both environments the two immune responses were positively genetically correlated: in full-sib families that were most likely to melanize a bead, injected bacteria were most likely to be cleared. Thus, our data do not support the idea of a trade-off among different outcomes of the invertebrate immune system, but rather that some families are overall immunologically superior to others.

Within the Diptera, two different selfish genetic elements are known to cause the production of female-biased sex ratios: maternally inherited bacteria that kill male zygotes (male-killers), and X chromosomes causing the degeneration of Y-bearing sperm in males (meiotic drive). We here develop a mathematical model for the dynamics of these two sex-ratio distorters where they co-occur. We show that X chromosome meiotic drive elements can be expected to substantially lower the equilibrium frequency of male-killers and can even lead to their extinction. Conversely, male-killers can also decrease the equilibrium frequency of X drivers and cause their extinction. Thus, we predict that there will be some complementarity in the incidence of X chromosome meiotic drive and male-killing in natural populations, with a lower than expected number of species bearing both elements.