Population size,habitat fragmentation,and the nature of adaptive variation in a stream fish

Whether and how habitat fragmentation and population size jointly affect adaptive genetic variation and adaptive population differentiation are largely unexplored. Owing to pronounced genetic drift, small, fragmented populations are thought to exhibit reduced adaptive genetic variation relative to large populations. Yet fragmentation is known to increase variability within and among habitats as population size decreases. Such variability might instead favour the maintenance of adaptive polymorphisms and/or generate more variability in adaptive differentiation at smaller population size. We investigated these alternative hypotheses by analysing coding-gene, single-nucleotide polymorphisms associated with different biological functions in fragmented brook trout populations of variable sizes. Putative adaptive differentiation was greater between small and large populations or among small populations than among large populations. These trends were stronger for genetic population size measures than demographic ones and were present despite pronounced drift in small populations. Our results suggest that fragmentation affects natural selection and that the changes elicited in the adaptive genetic composition and differentiation of fragmented populations vary with population size. By generating more variable evolutionary responses, the alteration of selective pressures during habitat fragmentation may affect future population persistence independently of, and perhaps long before, the effects of demographic and genetic stochasticity are manifest.     

Evolutionary consequences of habitat fragmentation: population size and density affect selection on inflorescence size in a perennial herb

Habitat fragmentation is considered to be one of the major threats to biological diversity worldwide. To date, however, its consequences have mainly been studied in an ecological context, while little is known about its effects on evolutionary processes. In this study we examined whether habitat fragmentation affects selection on plant phenotypic traits via changes in plant-pollinator interactions, using the self-incompatible perennial herb Phyteuma spicatum. Specifically, we hypothesized that limited pollination service in small or low-density populations leads to increased selection for traits that attract pollinators. We recorded mean seed production per capsule and per plant as a measure of pollination intensity and assessed selection gradients (i.e., trait-fitness relationships) in 16 natural populations of varying size and density over 2 years. Mean seed production was not related to population size or density, except for a marginal significant effect of density on the mean number of seeds per capsule in 1 year. Linear selection for flowering time and synchrony was consistent across populations; relative fitness was higher in earlier flowering plants and in plants flowering synchronously with others. Selection on inflorescence size, however, varied among populations, and linear selection gradients for inflorescence size were negatively related to plant population size and density in 1 year. Selection for increased inflorescence size decreased with increasing population size and density. Contrary to our expectation this appeared not to be related to changes in pollination intensity (mean seed production was not related to population size or density in this year), but was rather likely linked to differences in some other component of the abiotic or biotic environment. In summary, our results show that habitat fragmentation may influence selection on plant phenotypic traits, thereby highlighting potential evolutionary consequences of human-induced environmental change.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Population size is weakly related to quantitative genetic variation and trait differentiation in a stream fish

How population size influences quantitative genetic variation and differentiation among natural, fragmented populations remains unresolved. Small, isolated populations might occupy poor quality habitats and lose genetic variation more rapidly due to genetic drift than large populations. Genetic drift might furthermore overcome selection as population size decreases. Collectively, this might result in directional changes in additive genetic variation (V_A) and trait differentiation (Q_ST) from small to large population size. Alternatively, small populations might exhibit larger variation in V_A and Q_ST if habitat fragmentation increases variability in habitat types. We explored these alternatives by investigating V_A and Q_ST using nine fragmented populations of brook trout varying 50‐fold in census size N (179–8416) and 10‐fold in effective number of breeders, N_b (18–135). Across 15 traits, no evidence was found for consistent differences in V_A and Q_ST with population size and almost no evidence for increased variability of V_A or Q_ST estimates at small population size. This suggests that (i) small populations of some species may retain adaptive potential according to commonly adopted quantitative genetic measures and (ii) populations of varying sizes experience a variety of environmental conditions in nature, however extremely large studies are likely required before any firm conclusions can be made.     

Recent local adaptation of sockeye salmon to glacial spawning habitats

Salmonids spawn in highly diverse habitats, exhibit strong genetic population structuring, and can quickly colonize newly created habitats with few founders. Spawning traits often differ among populations, but it is largely unknown if these differences are adaptive or due to genetic drift. To test if sockeye salmon (Oncorhynchus nerka) populations are adapted to glacial, beach, and tributary spawning habitats, we examined variation in heritable phenotypic traits associated with spawning in 13 populations of wild sockeye salmon in Lake Clark, Alaska. These populations were commonly founded between 100 and 400 hundred sockeye salmon generations ago and exhibit low genetic divergence at 11 microsatellite loci (F _ST < 0.024) that is uncorrelated with spawning habitat type. We found that mean P _ST (phenotypic divergence among populations) exceeded neutral F _ST for most phenotypic traits measured, indicating that phenotypic differences among populations could not be explained by genetic drift alone. Phenotypic divergence among populations was associated with spawning habitat differences, but not with neutral genetic divergence. For example, female body color was lighter and egg color was darker in glacial than non-glacial habitats. This may be due to reduced sexual selection for red spawning color in glacial habitats and an apparent trade-off in carotenoid allocation to body and egg color in females. Phenotypic plasticity is an unlikely source of phenotypic differences because Lake Clark sockeye salmon spend nearly all their lives in a common environment. Our data suggest that Lake Clark sockeye salmon populations are adapted to spawning in glacial, beach and tributary habitats and provide the first evidence of a glacial spawning ecotype in salmonids. Glacial spawning habitats are often young (i.e., <200 years old) and ephemeral. Thus, local adaptation of sockeye salmon to glacial habitats appears to have occurred recently.     

The effects of selection,drift and genetic variation on life‐history trait divergence among insular populations of natterjack toad,Bufo calamita

Although loss of genetic variation is frequently assumed to be associated with loss of adaptive potential, only few studies have examined adaptation in populations with little genetic variation. On the Swedish west coast, the northern fringe populations of the natterjack toad Bufo calamita inhabit an atypical habitat consisting of offshore rock islands. There are strong among‐population differences in the amount of neutral genetic variation, making this system suitable for studies on mechanisms of trait divergence along a gradient of within‐population genetic variation. In this study, we examined the mechanisms of population divergence using Q_ST–F_ST comparisons and correlations between quantitative and neutral genetic variation. Our results suggest drift or weak stabilizing selection across the six populations included in this study, as indicated by low Q_ST–F_ST values, lack of significant population × temperature interactions and lack of significant differences among the islands in breeding pond size. The six populations included in this study differed in both neutral and quantitative genetic variation. Also, the correlations between neutral and quantitative genetic variation tended to be positive, however, the relatively small number of populations prevents any strong conclusions based on these correlations. Contrary to the majority of Q_ST–F_ST comparisons, our results suggest drift or weak stabilizing selection across the examined populations. Furthermore, the low heritability of fitness‐related traits may limit evolutionary responses in some of the populations.     

The role of environment and core‐margin effects on range‐wide phenotypic variation in a montane grasshopper

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range‐wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal‐related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120 000 years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among‐population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal‐related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range‐wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.     

Geographic variation in advertisement calls of a Microhylid frog – testing the role of drift and ecology

Acoustic signals for mating are important traits that could drive population differentiation and speciation. Ecology may play a role in acoustic divergence through direct selection (e.g., local adaptation to abiotic environment), constraint of correlated traits (e.g., acoustic traits linked to another trait under selection), and/or interspecific competition (e.g., character displacement). However, genetic drift alone can also drive acoustic divergence. It is not always easy to differentiate the role of ecology versus drift in acoustic divergence. In this study, we tested the role of ecology and drift in shaping geographic variation in the advertisement calls of Microhyla fissipes. We examined three predictions based on ecological processes: (1) the correlation between temperature and call properties across M. fissipes populations; (2) the correlation between call properties and body size across M. fissipes populations; and (3) reproductive character displacement (RCD) in call properties between M. fissipes populations that are sympatric with and allopatric to a congener M. heymonsi. To test genetic drift, we examined correlations among call divergence, geographic distance, and genetic distance across M. fissipes populations. We recorded the advertisement calls from 11 populations of M. fissipes in Taiwan, five of which are sympatrically distributed with M. heymonsi. We found geographic variation in both temporal and spectral properties of the advertisement calls of M. fissipes. However, the call properties were not correlated with local temperature or the callers' body size. Furthermore, we did not detect RCD. By contrast, call divergence, geographic distance, and genetic distance between M. fissipes populations were all positively correlated. The comparisons between phenotypic Q_st (P_st) and F_st values did not show significant differences, suggesting a role of drift. We concluded that genetic drift, rather than ecological processes, is the more likely driver for the geographic variation in the advertisement calls of M. fissipes.     

Anthropogenic fragmentation may not alter pre‐existing patterns of genetic diversity and differentiation in perennial shrubs

Many plant species have pollination and seed dispersal systems and evolutionary histories that have produced strong genetic structuring. These genetic patterns may be consistent with expectations following recent anthropogenic fragmentation, making it difficult to detect fragmentation effects if no prefragmentation genetic data are available. We used microsatellite markers to investigate whether severe habitat fragmentation may have affected the structure and diversity of populations of the endangered Australian bird‐pollinated shrub Grevillea caleyi R.Br., by comparing current patterns of genetic structure and diversity with those of the closely related G. longifolia R.Br. that has a similar life history but has not experienced anthropogenic fragmentation. Grevillea caleyi and G. longifolia showed similar and substantial population subdivision at all spatial levels (global F′_ST = 0.615 and 0.454; S_p = 0.039 and 0.066), marked isolation by distance and large heterozygous deficiencies. These characteristics suggest long‐term effects of inbreeding in self‐compatible species that have poor seed dispersal, limited connectivity via pollen flow and undergo population bottlenecks because of periodic fires. Highly structured allele size distributions, most notably in G. caleyi, imply historical processes of drift and mutation were important in isolated subpopulations. Genetic diversity did not vary with population size but was lower in more isolated populations for both species. Through this comparison, we reject the hypothesis that anthropogenic fragmentation has impacted substantially on the genetic composition or structure of G. caleyi populations. Our results suggest that highly self‐compatible species with limited dispersal may be relatively resilient to the genetic changes predicted to follow habitat fragmentation.     

Local adaptation and ecological differentiation under selection,migration, and drift in Arabidopsis lyrata*

How the balance between selection, migration, and drift influences the evolution of local adaptation has been under intense theoretical scrutiny. Yet, empirical studies that relate estimates of local adaptation to quantification of gene flow and effective population sizes have been rare. Here, we conducted a reciprocal transplant trial, a common garden trial, and a whole‐genome‐based demography analysis to examine these effects among Arabidopsis lyrata populations from two altitudinal gradients in Norway. Demography simulations indicated that populations within the two gradients are connected by gene flow (0.1 < 4N_em < 11) and have small effective population sizes (N_e < 6000), suggesting that both migration and drift can counteract local selection. However, the three‐year field experiments showed evidence of local adaptation at the level of hierarchical multiyear fitness, attesting to the strength of differential selection. In the lowland habitat, local superiority was associated with greater fecundity, while viability accounted for fitness differences in the alpine habitat. We also demonstrate that flowering time differentiation has contributed to adaptive divergence between these locally adapted populations. Our results show that despite the estimated potential of gene flow and drift to hinder differentiation, selection among these A. lyrata populations has resulted in local adaptation.     

Testing hypotheses driving genetic structure in the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla

Habitat fragmentation has often been implicated in the decline of many species. For habitat specialists and/or sedentary species, loss of habitat can result in population isolation and lead to negative genetic effects. However, factors other than fragmentation can often be important and also need to be considered when assessing the genetic structure of a species. We genotyped individuals from 13 populations of the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla in Florida to test three alternative hypotheses regarding the effects that habitat fragmentation might have on genetic structure. A map of potential habitat developed from recent satellite imagery suggested that Brown‐headed Nuthatch populations in southern Florida occupied smaller and more isolated habitat patches (i.e. were more fragmented) than populations in northern Florida. We also genotyped individuals from a small, isolated Brown‐headed Nuthatch population on Grand Bahama Island. We found that populations associated with more fragmented habitat in southern Florida had lower allelic richness than populations in northern Florida (P = 0.02), although there were no differences in heterozygosity. Although pairwise estimates of F_ST were low overall, values among southern populations were generally higher than northern populations. Population assignment tests identified K = 3 clusters corresponding to a northern cluster, a southern cluster and a unique population in southeast Florida; using sampling localities as prior information revealed K = 7 clusters, with greater structure only among southern Florida populations. The Bahamas population showed moderate to high differentiation compared with Florida populations. Overall, our results suggest that fragmentation could affect gene flow in Brown‐headed Nuthatch populations and is likely to become more pronounced over time.     

Population genetics and fitness in fragmented populations of the dioecious and endangered <Emphasis Type="Italic">Silene otites</Emphasis> (Caryophyllaceae)

Population fragmentation is often correlated with loss of genetic diversity and reduced fitness. Obligate out-crossing (dioecy) is expected to enhance genetic diversity, reduce genetic differentiation, and avoid inbreeding depression through frequent gene flow. However, in highly fragmented populations dioecy has only diminishing effects upon genetic structure as pollination limitations (e.g. flight distance of pollinators) most often restrict inter-population gene flow in insect pollinated species. In fragmented dry grasslands in northeastern Germany, we analysed genetic structure, fitness, and habitat quality of the endangered dioecious Silene otites (Caryophyllaceae). Using AFLP markers, a high level of differentiation among ten populations was found (F _st = 0.36), while the intra-population genetic diversities (H _E = 0.165–0.240) were similar as compared to hermaphroditic species. There was neither a correlation between geographic and genetic distance nor between genetic diversity and population size, which indicates reduced gene flow among populations and random genetic drift. Plant size was positively correlated with genetic diversity. Seed set and number of juveniles were positively related to population size. Higher total coverage resulted in reduced plant fitness, and the number of juveniles was negatively correlated to cryptogam cover. Additionally, we found a sex ratio bias towards more male plants in larger populations. Overall, our results indicate that on a regional geographic scale dioecy does not necessarily prevent genetic erosion in the case of habitat fragmentation, especially in the absence of long distance seed and pollen dispersal capacity.     

Gene flow does not prevent personality and morphological differentiation between two blue tit populations

Understanding the causes and consequences of population phenotypic divergence is a central goal in ecology and evolution. Phenotypic divergence among populations can result from genetic divergence, phenotypic plasticity or a combination of the two. However, few studies have deciphered these mechanisms for populations geographically close and connected by gene flow, especially in the case of personality traits. In this study, we used a common garden experiment to explore the genetic basis of the phenotypic divergence observed between two blue tit (Cyanistes caeruleus) populations inhabiting contrasting habitats separated by 25 km, for two personality traits (exploration speed and handling aggression), one physiological trait (heart rate during restraint) and two morphological traits (tarsus length and body mass). Blue tit nestlings were removed from their population and raised in a common garden for up to 5 years. We then compared adult phenotypes between the two populations, as well as trait‐specific Q_st and F_st. Our results revealed differences between populations similar to those found in the wild, suggesting a genetic divergence for all traits. Q_st–F_st comparisons revealed that the trait divergences likely result from dissimilar selection patterns rather than from genetic drift. Our study is one of the first to report a Q_st–F_st comparison for personality traits and adds to the growing body of evidence that population genetic divergence is possible at a small scale for a variety of traits including behavioural traits.     

Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

Genetic variation and effective population size in isolated populations of coastal cutthroat trout

Following glacial recession in southeast Alaska, waterfalls created by isostatic rebound have isolated numerous replicate populations of coastal cutthroat trout (Oncorhynchus clarkii clarkii) in short coastal streams. These replicate isolated populations offer an unusual opportunity to examine factors associated with the maintenance of genetic diversity. We used eight microsatellites to examine genetic variation within and differentiation among 12 population pairs sampled from above and below these natural migration barriers. Geological evidence indicated that the above-barrier populations have been isolated for 8,000–12,500 years. Genetic differentiation among below-barrier populations (F _ST = 0.10, 95% C.I. 0.08–0.12) was similar to a previous study of more southern populations of this species. Above-barrier populations were highly differentiated from adjacent below-barrier populations (mean pairwise F _ST = 0.28; SD 0.18) and multiple lines of evidence were consistent with asymmetric downstream gene flow that varied among streams. Each above-barrier population had reduced within-population genetic variation when compared to the adjacent below-barrier population. Within-population genetic diversity was significantly correlated with the amount of available habitat in above-barrier sites. Increased genetic differentiation of above-barrier populations with lower genetic diversity suggests that genetic drift has been the primary cause of genetic divergence. Long-term estimates of N _e based on loss of heterozygosity over the time since isolation were large (3,170; range 1,077–7,606) and established an upper limit for N _e if drift were the only evolutionary process responsible for loss of genetic diversity. However, it is likely that a combination of mutation, selection, and gene flow have also contributed to the genetic diversity of above-barrier populations. Contemporary above-barrier N _e estimates were much smaller than long-term N _e estimates, not correlated with within-population genetic diversity, and not consistent with the amount of genetic variation retained, given the approximate 10,000-year period of isolation. The populations isolated by waterfalls in this study that occur in larger stream networks have retained substantial genetic variation, which suggests that the amount of habitat in headwater streams is an important consideration for maintaining the evolutionary potential of isolated populations.     

Whole‐genome resequencing of extreme phenotypes in collared flycatchers highlights the difficulty of detecting quantitative trait loci in natural populations

Dissecting the genetic basis of phenotypic variation in natural populations is a long‐standing goal in evolutionary biology. One open question is whether quantitative traits are determined only by large numbers of genes with small effects, or whether variation also exists in large‐effect loci. We conducted genomewide association analyses of forehead patch size (a sexually selected trait) on 81 whole‐genome‐resequenced male collared flycatchers with extreme phenotypes, and on 415 males sampled independent of patch size and genotyped with a 50K SNP chip. No SNPs were genomewide statistically significantly associated with patch size. Simulation‐based power analyses suggest that the power to detect large‐effect loci responsible for 10% of phenotypic variance was <0.5 in the genome resequencing analysis, and <0.1 in the SNP chip analysis. Reducing the recombination by two‐thirds relative to collared flycatchers modestly increased power. Tripling sample size increased power to >0.8 for resequencing of extreme phenotypes (N = 243), but power remained <0.2 for the 50K SNP chip analysis (N = 1245). At least 1 million SNPs were necessary to achieve power >0.8 when analysing 415 randomly sampled phenotypes. However, power of the 50K SNP chip to detect large‐effect loci was nearly 0.8 in simulations with a small effective population size of 1500. These results suggest that reliably detecting large‐effect trait loci in large natural populations will often require thousands of individuals and near complete sampling of the genome. Encouragingly, far fewer individuals and loci will often be sufficient to reliably detect large‐effect loci in small populations with widespread strong linkage disequilibrium.     

Sexual selection and population divergence II. Divergence in different sexual traits and signal modalities in field crickets (Teleogryllus oceanicus)

Sexual selection can target many different types of traits. However, the relative influence of different sexually selected traits during evolutionary divergence is poorly understood. We used the field cricket Teleogryllus oceanicus to quantify and compare how five traits from each of three sexual signal modalities and components diverge among allopatric populations: male advertisement song, cuticular hydrocarbon (CHC) profiles and forewing morphology. Population divergence was unexpectedly consistent: we estimated the among‐population (genetic) variance‐covariance matrix, D , for all 15 traits, and D_max explained nearly two‐thirds of its variation. CHC and wing traits were most tightly integrated, whereas song varied more independently. We modeled the dependence of among‐population trait divergence on genetic distance estimated from neutral markers to test for signatures of selection versus neutral divergence. For all three sexual trait types, phenotypic variation among populations was largely explained by a neutral model of divergence. Our findings illustrate how phenotypic integration across different types of sexual traits might impose constraints on the evolution of mating isolation and divergence via sexual selection.     

Effect of migration and environmental heterogeneity on the maintenance of quantitative genetic variation: a simulation study

The paradox of high genetic variation observed in traits under stabilizing selection is a long‐standing problem in evolutionary theory, as mutation rates appear too low to explain observed levels of standing genetic variation under classic models of mutation–selection balance. Spatially or temporally heterogeneous environments can maintain more standing genetic variation within populations than homogeneous environments, but it is unclear whether such conditions can resolve the above discrepancy between theory and observation. Here, we use individual‐based simulations to explore the effect of various types of environmental heterogeneity on the maintenance of genetic variation (V_A) for a quantitative trait under stabilizing selection. We find that V_A is maximized at intermediate migration rates in spatially heterogeneous environments and that the observed patterns are robust to changes in population size. Spatial environmental heterogeneity increased variation by as much as 10‐fold over mutation–selection balance alone, whereas pure temporal environmental heterogeneity increased variance by only 45% at max. Our results show that some combinations of spatial heterogeneity and migration can maintain considerably more variation than mutation–selection balance, potentially reconciling the discrepancy between theoretical predictions and empirical observations. However, given the narrow regions of parameter space required for this effect, this is unlikely to provide a general explanation for the maintenance of variation. Nonetheless, our results suggest that habitat fragmentation may affect the maintenance of V_A and thereby reduce the adaptive capacity of populations.     

Toxic hydrogen sulfide and dark caves: phenotypic and genetic divergence across two abiotic environmental gradients in Poecilia mexicana

Divergent natural selection drives evolutionary diversification. It creates phenotypic diversity by favoring developmental plasticity within populations or genetic differentiation and local adaptation among populations. We investigated phenotypic and genetic divergence in the livebearing fish Poecilia mexicana along two abiotic environmental gradients. These fish typically inhabit nonsulfidic surface rivers, but also colonized sulfidic and cave habitats. We assessed phenotypic variation among a factorial combination of habitat types using geometric and traditional morphometrics, and genetic divergence using quantitative and molecular genetic analyses. Fish in caves (sulfidic or not) exhibited reduced eyes and slender bodies. Fish from sulfidic habitats (surface or cave) exhibited larger heads and longer gill filaments. Common-garden rearing suggested that these morphological differences are partly heritable. Population genetic analyses using microsatellites as well as cytochrome b gene sequences indicate high population differentiation over small spatial scale and very low rates of gene flow, especially among different habitat types. This suggests that divergent environmental conditions constitute barriers to gene flow. Strong molecular divergence over short distances as well as phenotypic and quantitative genetic divergence across habitats in directions classic to fish ecomorphology suggest that divergent selection is structuring phenotypic variation in this system.     

Divergent selection along climatic gradients in a rare central European endemic species,Saxifraga sponhemica

Background and Aims The effects of habitat fragmentation on quantitative genetic variation in plant populations are still poorly known. Saxifraga sponhemica is a rare endemic of Central Europe with a disjunct distribution, and a stable and specialized habitat of treeless screes and cliffs. This study therefore used S. sponhemica as a model species to compare quantitative and molecular variation in order to explore (1) the relative importance of drift and selection in shaping the distribution of quantitative genetic variation along climatic gradients; (2) the relationship between plant fitness, quantitative genetic variation, molecular genetic variation and population size; and (3) the relationship between the differentiation of a trait among populations and its evolvability.Methods Genetic variation within and among 22 populations from the whole distribution area of S. sponhemica was studied using RAPD (random amplified polymorphic DNA) markers, and climatic variables were obtained for each site. Seeds were collected from each population and germinated, and seedlings were transplanted into a common garden for determination of variation in plant traits.Key Results In contrast to previous results from rare plant species, strong evidence was found for divergent selection. Most population trait means of S. sponhemica were significantly related to climate gradients, indicating adaptation. Quantitative genetic differentiation increased with geographical distance, even when neutral molecular divergence was controlled for, and Q_ST exceeded F_ST for some traits. The evolvability of traits was negatively correlated with the degree of differentiation among populations (Q_ST), i.e. traits under strong selection showed little genetic variation within populations. The evolutionary potential of a population was not related to its size, the performance of the population or its neutral genetic diversity. However, performance in the common garden was lower for plants from populations with reduced molecular genetic variation, suggesting inbreeding depression due to genetic erosion.Conclusions The findings suggest that studies of molecular and quantitative genetic variation may provide complementary insights important for the conservation of rare species. The strong differentiation of quantitative traits among populations shows that selection can be an important force for structuring variation in evolutionarily important traits even for rare endemic species restricted to very specific habitats.