AN ANALYSIS OF LOCAL VARIABILITY OF FLOWER COLOR IN LINANTHUS PARRYAE

Genetic consequence of restricted habitat and population decline in endangered Isoetes sinensis (Isoetaceae)

BACKGROUND AND AIMS

Isoetes sinensis (Isoeteaceae) is a critically endangered aquatic quillwort in eastern China. Rapid decline of extant population size and local population extinction have occurred in recent years and have raised great concerns among conservationists.

METHODS

Amplified fragment length polymorphisms (AFLPs) were used to investigate the genetic variation and population structure of seven extant populations of the species.

KEY RESULTS

Eight primer combinations produced a total of 343 unambiguous bands of which 210 (61.2 %) were polymorphic. Isoetes sinensis exhibited a high level of intra-population genetic diversity (H(E) = 0.118; hs = 0.147; I = 0.192; P = 35.2 %). The genetic variation within each of the populations was not positively correlated with their size, suggesting recent population decline, which is well in accordance with field data of demographic surveys. Moreover, a high degree of genetic differentiation (F(ST) = 0.535; G(ST) = 0.608; theta(B) = 0.607) was detected among populations and no correlation was found between geographical and genetic distance, suggesting that populations were in disequilibrium of migration-drift. Genetic drift played a more important role than gene flow in the current population genetic structure of I. sinensis because migration of I. sinensis is predominantly water-mediated and habitat range was highly influenced by environment changes.

CONCLUSIONS

Genetic information obtained in the present study provides useful baseline data for formulating conservation strategies. Conservation management, including both reinforcement for in situ populations and ex situ conservation programmes should be carefully designed to avoid the potential risk of outbreeding depression by admixture of individuals from different regions. However, translocation within the same regional population should be considered as a measure of genetic enhancement to rehabilitate local populations. An ex situ conservation strategy for conserving all extant populations to maximize genomic representation of the species is also recommended.

Hierarchical analysis of population genetic structure in the monogynous ant Cataglyphis cursor using microsatellite and mitochondrial DNA markers

Despite having winged queens, female dispersal in the monogynous ant Cataglyphis cursor is likely to be restricted because colonies reproduce by fission. We investigated the pattern of population genetic structure of this species using eight microsatellite markers and a mitochondrial DNA (mtDNA) sequence, in order to examine the extent of female and nuclear gene flow in two types of habitat. Sampling was carried out at a large spatial scale (16 sites from 2.5 to 120 km apart) as well as at a fine spatial scale (two 4.5-km transects, one in each habitat type). The strong spatial clustering of mtDNA observed at the fine spatial scale strongly supported a restricted effective female dispersal. In agreement, patterns of the mtDNA haplotypes observed at large and fine spatial scales suggested that new sites are colonized by nearby sites. Isolation by distance and significant nuclear genetic structure have been detected at all the spatial scales investigated. The level of local genetic differentiation for mitochondrial marker was 15 times higher than for the nuclear markers, suggesting differences in dispersal pattern between the two sexes. However, male gene flow was not sufficient to prevent significant nuclear genetic differentiation even at short distances (500 m). Isolation-by-distance patterns differed between the two habitat types, with a linear decrease of genetic similarities with distance observed only in the more continuous of the two habitats. Finally, despite these low dispersal capacities and the potential use of parthenogenesis to produce new queens, no signs of reduction of nuclear genetic diversity was detected in C. cursor populations.

Geographical distance and physical barriers shape the genetic structure of Eurasian red squirrels (Sciurus vulgaris) in the Italian Alps

Red squirrels (Sciurus vulgaris) are widely distributed throughout Eurasia, occurring in many types of coniferous and mixed-deciduous forests. In fragmented landscapes, small and partly isolated populations with low immigration rates show reduced genetic diversity, but reforestation can increase gene flow and restore levels of genetic variation in a few decades. No studies have so far investigated the genetic structure of red squirrel in large, continuous forests. The Italian Alps are presently characterized by almost continuous, recently reconnected forest habitats, that were affected by deep landscape changes during last glaciations but remained mostly unchanged between 10 000 and 200 years bp, when forest cover was heavily reduced. In this study we analyse patterns of genetic variability of red squirrels in and between seven sites distributed over 250 km of Alpine habitat, using mitochondrial DNA (mtDNA) and microsatellites. We use isolation-by-distance (IBD) models to investigate the relative importance that past (Pleistocene glaciations) and recent (fragmentation, bottlenecks) events had on the present genetic situation. Both nuclear and mtDNA data indicate a significant differentiation among study sites and a significant correlation between genetic and geographical distance only over a large scale. No recent bottlenecks are recorded through microsatellites and demographic models strongly support equilibrium between gene flow and drift; however, mtDNA suggests that there may have been local demographic crashes, probably in correspondence with the 19th-century forest fragmentation. These findings indicate that local landscape factors other than geographical distance per se, such as barriers of unsuitable habitat, affect gene flow and determine differentiation.

Gene flow and inbreeding depression inferred from fine-scale genetic structure in an endangered heterostylous perennial, Primula sieboldii

We estimated the gene dispersal distance and the magnitude of inbreeding depression from the fine-scale genetic structure in the endangered heterostylous perennial Primula sieboldii. We indirectly estimated the neighbourhood size (Nb) and the standard deviation of gene dispersal distance (sigma(g)) from the detected genetic structure by using 10 microsatellite markers. We also estimated the fitness reduction in mating among neighbouring individuals caused by biparental inbreeding according to the genetic structure. We found clear fine-scale genetic structure (a significantly positive kinship coefficient within 42.3 m), and the indirect estimates of sigma(g) and Nb were 15.7 m and 50.9, respectively. These indirect estimates were similar to the direct estimates (18.4 m and 44.0). The slightly larger indirect estimate of Nb may reflect that inbreeding depression and genetic structure or rare long-distance dispersal that were overlooked in the direct estimate have elongated the long-term average of gene dispersal distance. P. sieboldii is also likely to suffer about 19% fitness reduction in progenies from mating among individuals 5 m apart. Our results suggest that biparental inbreeding and genetic structure can affect the range of gene dispersal and seed reproductive success in P. sieboldii.

Opposing natural selection from herbivores and pathogens may maintain floral-color variation in Claytonia virginica (Portulacaceae)

The maintenance of floral-color variation within natural populations is enigmatic because directional selection through pollinator preferences combined with random genetic drift should lead to the rapid loss of such variation. Fluctuating, balancing, and negative frequency-dependent selection mediated through pollinators have been identified as factors that may contribute to the maintenance of floral-color variation, and recently it has been suggested that indirect responses to selection on correlated characters through agents of selection other than pollinators may substantially shape the evolution of floral traits. Here, I provide empirical support for this latter view in Claytonia virginica (Portulacaceae) through a multiseason field study, a pollen supplementation study, and an artificial herbivory experiment. These studies show that most individuals fall into one of four discrete color classes, and suggest pollinator-mediated selection for increased floral redness in concurrent years. Floral color is also an indirect target of opposing directional selection via herbivores and pathogens that fluctuates through time. Taken together, these data suggest a novel mechanism by which floral-color variation may be maintained, and illustrate the importance of an inclusive, pluralistic view of selection when investigating the evolution of complex phenotypes.

Diversification before the most recent glaciation in Balanus glandula

A deep genetic cline between southern populations of the barnacle Balanus glandula (from about Monterey Bay southward) and northern populations (from northern California through Alaska) has recently been described. If this pattern is due to historical isolation and genetic drift, we expect it to have formed recently and represent a transient, nonequilibrium state. However, this cline appears to have formed well before the last glacial maximum. Our assays of sequence diversity at a region of mitochondrial cytochrome oxidase I, combined with coalescent estimators of the time of separation for these two regions, suggest that a late Pleistocene event more than 100 thousand years ago may be responsible for the initial separation. This suggests that either strong oceanographic mechanisms or natural selection have maintained the cline, because there has clearly been adequate time for this cline or polymorphism to resolve itself by genetic drift and migration. However, reliance on only a single mitochondrial marker for which the substitution rate has been estimated still limits the resolution of our analysis.

Population genetic structure and phylogeographical pattern of a relict tree fern, Alsophila spinulosa (Cyatheaceae), inferred from cpDNA atpB- rbcL intergenic spacers

Sequences of chloroplast DNA (cpDNA) atpB- rbcL intergenic spacers of individuals of a tree fern species, Alsophila spinulosa, collected from ten relict populations distributed in the Hainan and Guangdong provinces, and the Guangxi Zhuang region in southern China, were determined. Sequence length varied from 724 bp to 731 bp, showing length polymorphism, and base composition was with high A+T content between 63.17% and 63.95%. Sequences were neutral in terms of evolution (Tajima's criterion D=-1.01899, P>0.10 and Fu and Li's test D*=-1.39008, P>0.10; F*=-1.49775, P>0.10). A total of 19 haplotypes were identified based on nucleotide variation. High levels of haplotype diversity (h=0.744) and nucleotide diversity (Dij=0.01130) were detected in A. spinulosa, probably associated with its long evolutionary history, which has allowed the accumulation of genetic variation within lineages. Both the minimum spanning network and neighbor-joining trees generated for haplotypes demonstrated that current populations of A. spinulosa existing in Hainan, Guangdong, and Guangxi were subdivided into two geographical groups. An analysis of molecular variance indicated that most of the genetic variation (93.49%, P<0.001) was partitioned among regions. Wright's isolation by distance model was not supported across extant populations. Reduced gene flow by the Qiongzhou Strait and inbreeding may result in the geographical subdivision between the Hainan and Guangdong + Guangxi populations (FST=0.95, Nm=0.03). Within each region, the star-like pattern of phylogeography of haplotypes implied a population expansion process during evolutionary history. Gene genealogies together with coalescent theory provided significant information for uncovering phylogeography of A. spinulosa.

Partitioning nuclear and chloroplast variation at multiple spatial scales in the neotropical epiphytic orchid, Laelia rubescens

Insights into processes that lead to the distribution of genetic variation within plant species require recognition of the importance of both pollen and seed movement. Here we investigate the contributions of pollen and seed movement to overall gene flow in the Central American epiphytic orchid, Laelia rubescens. Genetic diversity and structure were examined at multiple spatial scales in the tropical dry forest of Costa Rica using nuclear (allozymes) and chloroplast restriction fragment length polymorphism (RFLP) markers, which were found to be diverse (allozymes, P = 73.3%; HE = 0.174; cpDNA, HE = 0.741). Nuclear genetic structure (FSTn) was low at every spatial scale (0.005-0.091). Chloroplast markers displayed more structure (0.073-0.254) but relatively similar patterns. Neither genome displayed significant isolation-by-distance. Pollen and seed dispersal rates did not differ significantly from one another (mp/ms = 1.40) at the broadest geographical scale, among sites throughout Costa Rica. However, relative contributions of pollen and seeds to gene flow were scale-dependent, with different mechanisms determining the dominant mode of gene flow at different spatial scales. Much seed dispersal is highly localized within the maternal population, while some seeds enter the air column and are dispersed over considerable distances. At the intermediate scale (10s to 100s of metres) pollinators are responsible for substantial pollen flow. This species appears capable of distributing its genes across the anthropogenically altered landscape that now characterizes its Costa Rican dry forest habitat.

Microsatellite variation reveals high levels of genetic variability and population structure in the gorgonian coral Pseudopterogorgia elisabethae across the Bahamas

The primary mechanism of gene flow in marine sessile invertebrates is larval dispersal. In Pseudopterogorgia elisabethae, a commercially important Caribbean gorgonian coral, a proportion of the larvae drop to the substratum within close proximity to the maternal colony, and most matings occur between individuals in close proximity to each other. Such limited dispersal of reproductive propagules suggests that gene flow is limited in this gorgonian. In this study, we characterized the population genetic structure of P. elisabethae across the Bahamas using six microsatellite loci. P. elisabethae was collected from 18 sites across the Bahamas. Significant deviations from Hardy-Weinberg equilibrium due to deficits of heterozygotes within populations were detected for all 18 populations in at least one of the six screened loci. Levels of genetic structure among populations of P. elisabethae were high and significant. A distance analysis placed populations within three groups, one formed by populations located within Exuma Sound, a semi-isolated basin, another consisting of populations located outside the basin and a third group comprising two populations from San Salvador Island. The patterns of genetic variation found in this study are concordant with the life-history traits of the species and in part with the geography of the Bahamas. Conservation and management plans developed for P. elisabethae should considered the high degree of genetic structure observed among populations of the species, as well as the high genetic diversity found in the San Salvador and the Exuma Sound populations.

Global Population Genetic Structure and Male-Mediated Gene Flow in the Green Sea Turtle (Chelonia mydas): Analysis of Microsatellite Loci

We assessed the degree of population subdivision among global populations of green sea turtles, Chelonia mydas, using four microsatellite loci. Previously, a single-copy nuclear DNA study indicated significant male-mediated gene flow among populations alternately fixed for different mitochondrial DNA haplotypes and that genetic divergence between populations in the Atlantic and Pacific Oceans was more common than subdivisions among populations within ocean basins. Even so, overall levels of variation at single-copy loci were low and inferences were limited. Here, the markedly more variable microsatellite loci confirm the presence of male-mediated gene flow among populations within ocean basins. This analysis generally confirms the genetic divergence between the Atlantic and Pacific. As with the previous study, phylogenetic analyses of genetic distances based on the microsatellite loci indicate a close genetic relationship among eastern Atlantic and Indian Ocean populations. Unlike the single-copy study, however, the results here cannot be attributed to an artifact of general low variability and likely represent recent or ongoing migration between ocean basins. Sequence analyses of regions flanking the microsatellite repeat reveal considerable amounts of cryptic variation and homoplasy and significantly aid in our understanding of population connectivity. Assessment of the allele frequency distributions indicates that at least some of the loci may not be evolving by the stepwise mutation model.

Genetic structure of the deep-sea coral Lophelia pertusa in the northeast Atlantic revealed by microsatellites and internal transcribed spacer sequences

The azooxanthellate scleractinian coral Lophelia pertusa has a near-cosmopolitan distribution, with a main depth distribution between 200 and 1000 m. In the northeast Atlantic it is the main framework-building species, forming deep-sea reefs in the bathyal zone on the continental margin, offshore banks and in Scandinavian fjords. Recent studies have shown that deep-sea reefs are associated with a highly diverse fauna. Such deep-sea communities are subject to increasing impact from deep-water fisheries, against a background of poor knowledge concerning these ecosystems, including the biology and population structure of L. pertusa. To resolve the population structure and to assess the dispersal potential of this deep-sea coral, specific microsatellites markers and ribosomal internal transcribed spacer (ITS) sequences ITS1 and ITS2 were used to investigate 10 different sampling sites, distributed along the European margin and in Scandinavian fjords. Both microsatellite and gene sequence data showed that L. pertusa should not be considered as one panmictic population in the northeast Atlantic but instead forms distinct, offshore and fjord populations. Results also suggest that, if some gene flow is occurring along the continental slope, the recruitment of sexually produced larvae is likely to be strongly local. The microsatellites showed significant levels of inbreeding and revealed that the level of genetic diversity and the contribution of asexual reproduction to the maintenance of the subpopulations were highly variable from site to site. These results are of major importance in the generation of a sustainable management strategy for these diversity-rich deep-sea ecosystems.

The rise and fall of isolation by distance in the anadromous brook charr (Salvelinus fontinalis Mitchill)

Geographic patterns of genetic diversity depend on a species' demographic properties in a given habitat, which may change over time. The rates at which patterns of diversity respond to changes in demographic properties and approach equilibrium are therefore pivotal in our understanding of spatial patterns of diversity. The brook charr Salvelinus fontinalis is a coastal fish exhibiting limited marine movements, such that a stable one-dimensional isolation-by-distance (IBD) pattern should be observed over the whole range. Its range, however, recently shifted northward such that northern populations may still be in the process of reaching equilibrium. We investigated variation in IBD patterns, genetic divergence, and allelic richness at six microsatellite markers in 2087 anadromous brook charr from 59 rivers along the most likely postglacial colonization route. We observed a decrease in allelic richness, together with an increase in differentiation and a decrease in IBD in the most recently colonized northern populations, as expected following recent colonization. Contrary to expectation, however, similar patterns were also observed at the southernmost part of the range, despite the fact that these populations are not considered to be newly colonized. We propose that the loss of dispersal capabilities associated with anadromy may have caused the southernmost populations to evolve relatively independently of one another. This study thus demonstrated that changes in a species' geographic range and dispersal capabilities may contribute to shaping geographic patterns of genetic diversity.

Host-plant-associated genetic differentiation in Northern French populations of the European corn borer

The phytophagous insects that damage crops are often polyphagous, feeding on several types of crop and on weeds. The refuges constituted by noncrop host plants may be useful in managing the evolution in pest species of resistance to the Bacillus thuringiensis toxins produced by transgenic crops. However, the benefits of these refuges may be limited because host-plant diversity may drive genetic divergence and possibly even host-plant-mediated sympatric speciation. The European corn borer, Ostrinia nubilalis Hübner (Lepidoptera: Crambidae), is the main pest of maize in Europe and North America, where it was introduced early in the 20th century. It has a wide host range but feeds principally on mugwort (Artemisia vulgaris L.) and maize (Zea mays L.). O. nubilalis is found on mugwort only in the northern part of France, whereas it is found on maize throughout France. The extent of genetic variation at allozyme markers was investigated in populations collected from the two host plants over the entire geographical distribution of the European corn borer on mugwort in France. Allelic differentiation between pairs of populations and hierarchical analyses of pools of samples from each host plant indicate that the group of populations feeding on maize differed from the group of populations feeding on mugwort. Our results suggest (1) host-plant-related divergent selection at the genomic region surrounding the Mpi locus and (2) limited gene flow between the populations feeding on mugwort and those infesting maize fields. These data indicate that adults emerging from mugwort would not be useful for managing the evolution of resistance to the B. thuringiensis toxins in European corn borer populations.

Genetic structure of age classes in Camellia japonica (Theaceae)

Camellia japonica L. (Theaceae), an insect- and bird-pollinated, broad-leaved evergreen tree, is widely distributed in Japan and the southern Korean peninsula. The species has a relatively even age distribution within populations, which may influence the spatial genetic structure of different age classes relative to species with typical L-shaped age distributions. To determine whether the internal spatial genetic structure found in seedlings and young individuals carries over into adults, we used allozyme loci, F-statistics, spatial autocorrelation statistics (Moran's I), and coancestry measures to examine changes in genetic structure among seven age classes in a population (60-m x 100-m area) in southern Korea. In seedlings, weak but significant positive values of Moran's I-statistics and coancestry measures were found for distances less than 14 m, which is consistent with a mechanism of limited seed dispersal combined with overlapping seed shadows. This spatial structure, however, dissipates in older age classes, and in adults genetic variation has an essentially random spatial distribution. Morisita's index of dispersion of individuals in each age class showed that seedlings and juveniles are more highly clustered than are older individuals. These results suggest that self-thinning changes the spatial relationships of individuals, and thus genotypes. A multilocus estimate of FST (0.008) shows a small but statistically significant difference in allele frequencies among age classes. In summary, intrapopulation genetic structure within and among age classes of C. japonica was significant but weak. Despite presumably limited seed dispersal, weak spatial genetic structure in juveniles suggests overlapping seed shadows followed by self-thinning during recruitment. The present study also demonstrates that studies of spatial genetic structure focusing on limited numbers of generations may not be sufficient to reveal the entire picture of genetic structure in populations with overlapping generations.

Contrasting patterns of divergence in quantitative traits and neutral DNA markers: analysis of clinal variation

Clinal variation in quantitative traits is often attributed to the effects of spatially varying selection. However, identical patterns can be produced by the interplay between purely stochastic processes (i.e. drift in combination with spatially restricted gene flow). One means of distinguishing between adaptive and nonadaptive causes of geographical variation is to compare relative levels of between-population divergence in quantitative traits and neutral DNA markers. Such comparisons can be used to test whether levels of trait divergence attributable to additive genetic effects (as measured by QST) exceed null expectations based on the level of divergence at neutral marker loci (as measured by FST). The purpose of this study was to use an approach based on 'QST vs. FST' contrasts to test for evidence of diversifying selection on body size of an Indian fruit bat, Cynopterus sphinx (Chiroptera: Pteropodidae). Specifically, relative levels of between-population divergence in body size and microsatellite DNA markers were compared to assess whether the observed pattern of clinal size variation could be explained by a neutral model of isolation by distance. QST for body size was calculated using unbiased estimators of within- and between-population variance of principal component scores. The association between body size variation and geographical/environmental distance was tested using pairwise and partial matrix correspondence tests (MCTs). Independent variables (representing causal hypotheses) were constructed as between-locality distance matrices. The effects of neutral genetic divergence were assessed by including a matrix of pairwise FST as an independent variable. Partial MCTs revealed highly significant associations between phenotypic divergence (QST) and both geographical and environmental distance, even when the effects of neutral genetic divergence (FST) were partialled out. Results of the tests confirmed that migration-drift equilibrium is not a sufficient explanation for the latitudinal pattern of clinal size variation in C. sphinx. The geographical patterning of pairwise QST is most likely attributable to spatially varying selection and/or the direct influence of latitudinally ordered environmental effects.

Evaluation of the Amazon River delta as a barrier to gene flow for the regional malaria vector, Anopheles aquasalis (Diptera: Culicidae) in northeastern Brazil

The Neotropical malaria vector Anopheles aquasalis Curry is distributed predominantly along the Atlantic and Pacific Coasts because of its tolerance for breeding in salt water. We tested the hypothesis that the freshwater Amazon River acts as a barrier to gene flow in northeastern Brazil, by examining variation at a 588-nucleotide fragment of the mitochondrial cytochrome oxidase Igene from five populations. We identified 15 haplotypes of which 5 were shared both (1) between sample localities and (2) across the Amazon River Delta. Sequence divergence ranged from 0.0017-0.0272 (average = 0.0137). Estimates of genetic subdivision based on the presence of the Amazon Riverwere greatest within localities (phi = 0.029) and among regions (phi = 0.018), followed by among localities (phi = 0.011), but none were significant. Parsimony, neighbor-joining, and Nested Clade Analyses were used to estimate relationships among populations and infer evolutionary processes. Two phylogenetically distinct clusters of populations were moderately supported by parsimony. Neighbor-joining trees were poorly resolved, thus providing no geographical resolution and no support for the Amazon River as a barrier to migration. Phylogeographic structure as detected by the Nested Clade Analysis was consistent with restricted gene flow coupled with isolation by distance. Taken together, these analyses suggest that the localities within this region of northeastern Brazil constitute a single large population of An. aquasalis that spans the Amazon Delta.

Invasion of vacant niches and subsequent sympatric speciation

An individual-based simulation study was conducted to examine the population dynamics of 'invasion of a vacant niche' and subsequent speciation (by reproductive isolation) when food resources are randomly distributed spatially within the habitat and the frequencies of different food types are bimodally distributed (i.e. smaller and larger sizes of food being most abundant). The initially vacant niche was that of unused larger sizes of food. When phenotypic variation for resource use (i.e. food sizes) was small in the initial population, and each female could choose a mate from anywhere in the habitat, the population could not invade the vacant niche. But when the dispersal distance of the offspring and the area within which a female could choose a mate were small (i.e. the genetic neighbourhood size was small), the population could, in most cases, evolve to use both smaller and larger food sizes and form sister species sympatrically, with each species utilizing one of the two niches (small and large sizes of food). When phenotypic variation in resource use in the initial population was large, the population could, in most cases, invade the vacant niche by evolving to use both smaller and larger sizes of food. The probability of speciation increased as the dispersal distance of offspring decreased. The results indicate that populations whose individuals have small Wright's genetic neighbourhoods may often exploit a vacant niche and diversify sympatrically in the process.

The influence of neighborhood size and habitat shape on the accumulation of deleterious mutations

To examine the impact of genetic neighborhood size and habitat shape on genetic load and the accumulation of deleterious mutation, individual-based simulations were performed in continuously distributed habitats. The risk of extinction increased as both the area of the habitat and the neighborhood size decreased. When the neighborhood area became smaller than the habitat area, habitat shape also began to influence the risk of extinction by mutation loads, expected time to extinction being shorter in longer and narrower habitats than in a square habitat. Both the number of homozygous deleterious loci per individual and the mutation load in the population increased as the neighborhood size and total population size decreased. Neighborhood size and total population size both independently affected the average number of homozygous deleterious loci per individual. In addition, as the ratio of the long to the short side of the rectangle of a habitat increased, the average number of homozygous deleterious loci increased. When the areas of the habitats were held constant, the average number of homozygous loci and the mutation loads were smallest for a regular square and largest for the longest, narrowest habitat. These results suggest that the spatial genetic structure of an individual is an important factor in the accumulation of deleterious mutations and the risk of extinction by mutation meltdown.

Stable two-allele polymorphisms maintained by fluctuating fitnesses and seed banks: protecting the blues in Linanthus parryae

Motivated by data demonstrating fluctuating relative and absolute fitnesses for white- versus blue-flowered morphs of the desert annual Linanthus parryae, we present conditions under which temporally fluctuating selection and fluctuating contributions to a persistent seed bank will maintain a stable single-locus polymorphism. In L. parryae, blue flower color is determined by a single dominant allele. To disentangle the underlying diversity-maintaining mechanism from the mathematical complications associated with departures from Hardy-Weinberg genotype frequencies and dominance, we successively analyze a haploid model, a diploid model with three distinguishable genotypes, and a diploid model with complete dominance. For each model, we present conditions for the maintenance of a stable polymorphism, then use a diffusion approximation to describe the long-term fluctuations associated with these polymorphisms. Our protected polymorphism analyses show that a genotype whose arithmetic and geometric mean relative fitnesses are both less than one can persist if its relative fitness exceeds one in years that produce the most offspring. This condition is met by data from a population of L. parryae whose white morph has higher fitness (seed set) only in years of relatively heavy rain fall. The data suggest that the observed polymorphism may be explained by fluctuating selection. However, the yearly variation in flower color frequencies cannot be fully explained by our simple models, which ignore age structure and possible selection in the seed bank. We address two additional questions--one mathematical, the other biological--concerning the applicability of diffusion approximations to intense selection and the applicability of long-term predictions to datasets spanning decades for populations with long-lived seed banks.

Perspective: Evolution of flower color in the desert annual Linanthus parryae: Wright revisited

Linanthus parryae, a diminutive desert annual with white or blue flowers, has been the focus of a long-standing debate among evolutionary biologists. At issue is whether the flower color polymorphism in this species is the product of random genetic drift, as Sewall Wright argued, or of natural selection, as proposed by Carl Epling and his colleagues. Our long-term studies of three polymorphic populations in the Mojave Desert demonstrate that flower color is subject to selection that varies in both time and space in its direction and magnitude. For all sites taken together, blue-flowered plants produced more seeds than white-flowered plants in years of relatively low seed production, whereas white-flowered plants had higher fitness in years of high seed production. Evidence of selection on flower color was found in two of the three study sites. Differences in fitness between the color morphs were sometimes large, with selection coefficients as high as 0.60 in some years. Our longest period of observations was at Pearblossom site 1, where plants reached appreciable densities in seven of the 11 years of study. Here we found significant differences in the seed production of the color morphs in six years, with three years of blue advantage and three years of white advantage. For all sites taken together, total spring precipitation (March and April) was positively correlated with both absolute and relative seed production of the color morphs. At Pearblossom site 1, blue-flowered plants typically had a fitness advantage in years of low spring precipitation, whereas white-flowered plants had a fitness advantage in years of high spring precipitation. This temporal variation in selection may contribute to the maintenance of the flower-color polymorphism at Pearblossom site 1, whereas gene flow from neighboring populations is proposed as the principal factor maintaining the polymorphism at the other study sites. We found no significant differences between the color morphs in pollinator visitation rate or in their carbon isotope ratio, a measure of water-use efficiency. Although the mechanism of selection remains elusive, our results refute Wright's conclusion that the flower color polymorphism in L. parryae is an example of isolation by distance, a key component of his shifting balance theory of evolution.

Population genetics and colony structure of the Argentine ant (Linepithema humile) in its native and introduced ranges

Introduced species often possess low levels of genetic diversity relative to source populations as a consequence of the small population sizes associated with founder events. Additionally, native and introduced populations of the same species can possess divergent genetic structuring at both large and small geographic scales. Thus, genetic systems that have evolved in the context of high diversity may function quite differently in genetically homogeneous introduced populations. Here we conduct a genetic analysis of native and introduced populations of the Argentine ant (Linepithema humile) in which we show that the population-level changes that have occurred during introduction have produced marked changes in the social structure of this species. Native populations of the Argentine ant are characterized by a pattern of genetic isolation by distance, whereas this pattern is absent in introduced populations. These differences appear to arise both from the effects of recent range expansion in the introduced range as well as from differences in gene flow within each range. Relatedness within nests and colonies is lower in the introduced range than in the native range as a consequence of the widespread genetic similarity that typifies introduced populations. In contrast, nestmates and colony-mates in the native range are more closely related, and local genetic differentiation is evident. Our results shed light on the problem posed for kin selection theory by the low levels of relatedness that are characteristic of many unicolonial species and suggest that the loss of genetic variation may be a common mechanism for the transition to a unicolonial colony structure.

Isolation by distance in equilibrium and nonequilibrium populations of four talitrid species in the Mediterranean Sea

Allozymic variation at 21-23 loci was studied in 28 populations of Talitrus saltator, 23 populations of Orchestia montagui, 13 populations of O. stephenseni, and five populations of Platorchestia platensis from the Mediterranean Basin. Different levels of gene flow (Nmtheta) were detected within each species at the scale of the whole Mediterranean: O. montagui and P. platensis had low population structure, with levels of Nmtheta > or = 1, whereas the T. saltator and 0. stephenseni populations have values of Nmtheta < 1. The relationship between Nmtheta and geographic distance was analyzed to test for the presence of an isolation by distance pattern in the spatial genetic variation within each species. A model of isolation by distance is useful to describe the pattern of genetic structuring of study species at the scale of the whole Mediterranean: geographic distance explained from 28% to 70% of the variation in gene flow. In the Aegean area all species showed an island model of genetic structuring regardless of the levels of gene flow.

Effects of local geographic barriers and latitude on population structure in Anopheles punctipennis (Diptera: Culicidae)

We sampled Anopheles punctipennis (Say) from 11 localities throughout Vermont to examine the effects of latitude and two local geographical boundaries, Lake Champlain and the Green Mountains, on the population genetic structure of this species. Thirty-five mitochondrial haplotypes were detected in 104 individuals using a variable region of the COI gene. When latitude was examined, we detected significant structure within localities and among localities within latitudinal regions. For geographic analysis, significant genetic structure was detected only within localities. Estimates of gene flow across geographic regions indicate that the Green Mountains, but not Lake Champlain, is a barrier to dispersal for this species. We found no correlation between genetic and geographic distances for An. punctipennis.

Extensive gene flow within sibling species in the deep-sea fish Beryx splendens

Molecular markers allow insights into the population biology and ecology of deep-sea organisms, which are usually hardly accessible to direct observation and poorly known. Such a study was undertaken here for the deep-sea fish Beryx splendens, a species of growing interest to fisheries. B. splendens populations were sampled on seamounts and continental margins in the southwestern Pacific (New Caledonia, New Zealand, southeastern Australia) and in the northeastern Atlantic. Two hundred and fifty individuals were characterised by their single-strand DNA conformation (SSCP) of a approximately 360-base-pair (bp) fragment of the mitochondrial cytochrome b gene, amplified by the polymerase chain reaction (PCR). Two major SSCP haplotypes were observed in New Caledonia, a and w, whose frequencies were negatively correlated along a north-to-south cline. All SSCP haplotypes in the total sample were sequenced on 273 bp. The phylogenetic tree of B. splendens haplotype sequences, rooted by two B. decadactylus sequences, showed that a and w belong to distinct mitochondrial clades, A and W, which are separated by approximately 4-6% nucleotide divergence. Thirty individuals from New Caledonia were characterised by their DNA fingerprint from arbitrary-primed PCR. The distribution of individual-pairwise similarity indices was strongly bimodal. The larger similarity values all corresponded to comparisons within a clade (A or W) while the lower values were all between clades. Therefore, there was a strict association between the mitochondrial type and the DNA (presumably, nuclear DNA) fingerprint of an individual. Altogether, these results point to the existence of two biological species (sp. A and sp. W) within the current taxon B. splendens. No within-species differentiation was detected at the regional scale (New Caledonia). A remarkable result is that the three cytochrome b haplotypes of northeastern Atlantic B. cf. splendens sp. A were also the three commonest in the southwestern Pacific populations of this species. Such a level of homogeneity in the distribution of haplotypes suggests there is, or recently has been, gene flow at the inter-oceanic scale.

Competition and dispersal in predator-prey waves

Dispersing predators and prey can exhibit complex spatio-temporal wave-like patterns if the interactions between them cause oscillatory dynamics. We study the effect of these predator-prey density waves on the competition between prey populations and between predator populations with different dispersal strategies. We first describe 1- and 2-dimensional simulations of both discrete and continuous predator-prey models. The results suggest that any population that diffuses faster, disperses farther, or is more likely to disperse will exclude slower diffusing, shorter dispersing, or less likely dispersing populations, everything else being equal. It also appears that it does not matter whether time, space, or state are discrete or continuous, nor what the exact interactions between the predators and prey are. So long as waves exist the competition between populations occurs in a similar fashion. We derive a theory that qualitatively explains the observed behaviour and calculate approximate analytical solutions that describe, to a reasonable extent, these behaviours. Predictions about the cost of dispersal are tested. If strong enough, cost can reverse the populations' relative competitive strengths or lead to coexistence because of the effect of spiral wave cores. The theory is also able to explain previous results of simulations of coexistence in host-parasitoid models (Comins, H. N., and Massell, M. P., 1996, J. Theor. Biol. 183, 19-28).

Gene flow among habitat patches on a fragmented landscape in the spider argiope trifasciata (Araneae: araneidae)

The banded garden spider (Argiope trifasciata) is a common orb weaver of old field habitats in the United States. In this study, we determined levels of genetic variability and gene flow among seven populations in central Pennsylvania, based on variation at eight allozyme loci. Mean heterozygosity (observed) per population was 7. 5% and mean polymorphism was 39.3%, consistent with levels of variability in other arthropods. Values of GST for the four polymorphic loci (mean GST=0.011) suggest that gene flow prevents the genetic differentiation of these populations. The average number of migrants per generation (Nm) among these populations is estimated to be 31.3. The lack of significant interpopulation genetic differentiation among these disjunct populations may result from spiderling aerial dispersal (ballooning), a more continuous distribution of suitable habitat in the past, and perhaps the use of roadside vegetation as gene flow corridors. On the other hand, the study populations did not exhibit isolation-by-distance, suggesting that suitable habitat in our study area is experienced as less than continuous by A. trifasciata. Thus, although A. trifasciata is an excellent ballooner, ballooning does not confer unlimited access among all populations, which suggests that ballooning may be a far less effective means of long-distance dispersal than previously thought.

Local genetic structure in red grouse (Lagopus lagopus scoticus): evidence from microsatellite DNA markers

Allelic variation at seven hypervariable tri- and tetranucleotide microsatellite loci was used to determine levels of population differentiation between 14 populations of red grouse (Lagopus lagopus scoticus) in northeast Scotland, UK. Despite the potential for long-distance dispersal in grouse, and a semicontinuous habitat, significant population divergence was observed (mean RST = 0.153; P < 0.01) and an isolation-by-distance effect detected (Mantel test: P < 0.001). Examination of the spatial trend in principal component scores derived from allele frequencies among populations highlighted a barrier to gene flow that was confounding a simple isolation-by-distance effect. This barrier corresponded to an area of unsuitable habitat for grouse associated with a river system that bisected the study area. Mean genetic relatedness was higher for males than for females in all but one of the study populations, suggesting that the territorial behaviour and natal philopatry displayed by cocks have a manifold effect in generating the observed spatial genetic structure. Lower female relatedness values suggest a higher level of female-mediated gene flow, which is sufficient to prevent the loss of genetic variation from within populations and the onset of inbreeding effects. The potential consequences of local subdivision for red grouse populations are discussed.

Local genetic structure within two rookeries of Chelonia mydas (the green turtle)

We used multilocus minisatellite DNA fingerprinting to examine the local genetic structure within nesting populations of green turtles (Chelonia mydas) in Tortuguero, Costa Rica and Melbourne, Florida, USA. In the Tortuguero population, there was a significant negative correlation between genetic similarity of pairs of nesting females and the distance between their nest sites both within years (r2 = 0.273; P < 0.001) and between years (r2 = 0.578; P < 0.001). Of the 122 female pairs scored for Tortuguero, 12.3 per cent had genetic similarity values resembling those of mother-offspring pairs. In the Melbourne population, however, no relationship between genetic similarity and distance was found (r2 = 0.017; P = 0.075). The distance-related genetic structure of the Tortuguero population indicates that these females exhibit low levels of dispersal from natal sites, and that nestmates return independently to nest near their natal sites. The lack of a similar structure in the Melbourne population suggests that females from this population may not return to natal sites with comparable precision. High levels of mortality among nests, hatchlings or maturing turtles produced in the Melbourne rookery may also be responsible for the absence of distance-related local genetic structure.

Estimates of gene flow in Drosophila pseudoobscura determined from nucleotide sequence analysis of the alcohol dehydrogenase region

The genetic structure of Drosophila pseudoobscura populations was inferred from a nucleotide sequence analysis of a 3.4-kb segment of the alcohol dehydrogenase (Adh) region. A total of 99 isochromosomal strains collected from 13 populations in North and South America were used to determine if any population departed from a neutral model and to estimate levels of gene flow between populations. This study also included the nucleotide sequences from two sibling species, D. persimilis and D. miranda. We estimated the neutral mutation parameter, 4N mu, in synonymous and noncoding sites for 17 subregions of Adh in each of nine populations with sample sizes greater than three. The nucleotide diversity data in the nine populations was tested for departures from an equilibrium neutral model with two statistical tests. The Tajima and the Hudson, Kreitman, Aguade tests showed that each population fails to reject a neutral model. Tests for genetic differentiation between populations fail to show any population substructure among the North American populations of D. pseudoobscura. The nucleotide diversity data is consistent with direct and indirect measures of gene flow that show extensive dispersal between populations of D. pseudoobscura.