Abdominal pigmentation in Drosophila melanogaster females from natural Indian populations

Females of Drosophila melanogaster collected from five geographically distant populations in India were analysed for the intensity of pigmentation in the 5th, 6th and 7th segments of the abdomen. In all three segments, this intensity was found to vary among individuals of any given population, and, furthermore, different populations differ with respect to this phenotypic trait. Statistical analysis revealed significant intra- and interpopulational variation. A clinical pattern was detected: females from populations closer to the equator tended to have lighter cuticle, in which case differences between the three segments could not be detected and all three segments responded both independently and jointly to latitudinal variation, as indicated by a statistically significant positive correlation between latitude and pigmentation score. This is the first report on abdominal pigmentation analysis in natural populations of D. melanogaster that provides evidence that phenotypic flexibility reflects temperature differences, as a result of which abdominal pigmentation shows geographic differentiation.     

A Genome-Wide,Fine-Scale Map of Natural Pigmentation Variation in Drosophila melanogaster

Various approaches can be applied to uncover the genetic basis of natural phenotypic variation, each with their specific strengths and limitations. Here, we use a replicated genome-wide association approach (Pool-GWAS) to fine-scale map genomic regions contributing to natural variation in female abdominal pigmentation in Drosophila melanogaster, a trait that is highly variable in natural populations and highly heritable in the laboratory. We examined abdominal pigmentation phenotypes in approximately 8000 female European D. melanogaster, isolating 1000 individuals with extreme phenotypes. We then used whole-genome Illumina sequencing to identify single nucleotide polymorphisms (SNPs) segregating in our sample, and tested these for associations with pigmentation by contrasting allele frequencies between replicate pools of light and dark individuals. We identify two small regions near the pigmentation genes tan and bric-à-brac 1, both corresponding to known cis-regulatory regions, which contain SNPs showing significant associations with pigmentation variation. While the Pool-GWAS approach suffers some limitations, its cost advantage facilitates replication and it can be applied to any non-model system with an available reference genome.     

Skin and hair pigmentation variation in Island Melanesia

Skin and hair pigmentation are two of the most easily visible examples of human phenotypic variation. Selection-based explanations for pigmentation variation in humans have focused on the relationship between melanin and ultraviolet radiation, which is largely dependent on latitude. In this study, skin and hair pigmentation were measured as the melanin (M) index, using narrow-band reflectance spectroscopy for 1,135 individuals from Island Melanesia. Overall, the results show remarkable pigmentation variation, given the small geographic region surveyed. This variation is discussed in terms of differences between males and females, among islands, and among neighborhoods within those islands. The relationship of pigmentation to age, latitude, and longitude is also examined. We found that male skin pigmentation was significantly darker than females in 5 of 6 islands examined. Hair pigmentation showed a negative, but weak, correlation with age, while skin pigmentation showed a positive, but also weak, correlation with age. Skin and hair pigmentation varied significantly between islands as well as between neighborhoods within those islands. Bougainvilleans showed significantly darker skin than individuals from any other island considered, and are darker than a previously described African-American population. These findings are discussed in relation to prevailing hypotheses about the role of natural selection in shaping pigmentation variation in the human species, as well as the role of demographic processes such as admixture and drift in Island Melanesia.     

Relationships between environmental gradients and geographic variation in the intraspecific body size of three species of frogs (Anura)

The relationship between environmental gradients and patterns of geographic variation in body size has been a controversial topic for ectothermic organisms globally. To examine whether the patterns that generally hold in more temperate species also hold for tropical ones, we examined the intraspecific body size variation in three species of Neotropical frogs, Dendropsophus minutus, Hypsiboas faber and Physalaemus cuvieri, along different environmental gradients (e.g. temperature, precipitation and topography). We analysed four competing hypotheses: (i) the water availability hypothesis that predicts a negative relationship between body size and precipitation; (ii) the heat balance hypothesis that predicts a negative relationship between body size and temperature; (iii) the topography hypothesis that predicts a negative relationship between body size and altitude; and (iv) the mixed‐effect hypothesis that predicts that individuals occurring in wet and cold sites would be larger than individuals occurring in dry and warm sites. The spatial pattern of geographic variation in body size among populations of H. faber was associated with the mixed‐effect hypothesis. In localities with low precipitation seasonality and cold conditions, H. faber individuals were larger than in localities with high precipitation seasonality and warm conditions. Variation in the body size of D. minutus was the opposite of that predicted by the heat balance hypothesis. Individuals in localities with high temperatures were larger than in localities with low temperatures. On the other hand, variation in the body size of P. cuvieri was not associated with the variables used in this study. Our results suggest that intraspecific variation in anuran body size is more dependent on species‐specific response than on the region (i.e. temperate or tropical) where they occur.     

Spatial and environmental correlates of intraspecific morphological variation in three species of passerine birds from the Purus–Madeira interfluvium,Central Amazonia

Biogeographic studies in Amazonia typically describe biodiversity across interfluvia, rarely within them, where geographic variability in morphological traits might be observed. We tested for intraspecific phenotypic variation in three bird species within the Purus–Madeira interfluvium (Central Amazon) and whether phenotypes were correlated with environmental heterogeneity or geographic distance among sites. We compared coloration indexes derived from reflectance spectra and morphometrics of up to five adult individuals of each sex among 11 sites within the interfluvium and contrasted them with proxies for geographic distance and environmental variation (tree basal area and bird community). Environmental heterogeneity was minimally spatially autocorrelated, and there were no obvious geographical barriers to dispersal in the study region. The null hypothesis was that we would see either no phenotypic variation or random variation that was not explained by the tested variables. Half of the cases analyzed showed intraspecific morphological variation. Coloration varied more frequently than morphometrics, and color was better explained by environmental heterogeneity, particularly in males, whereas brightness also varied with geographic distance. Geographic distance explained the only case of variation in morphometrics. Our results indicate that coloration, particularly plumage brightness, is more labile than morphometric traits and that plumage color might be under stronger effects of local adaptation than brightness, which also seems to be under effects of neutral drift and gene flow among populations. Higher frequencies of association between male coloration and the environment suggest a role of non-arbitrary mechanisms of sexual selection on the expression of male phenotypes, whereas arbitrary intersexual selection might explain the randomly distributed variation that is not explained by environmental heterogeneity or geographic distance. We revealed intraspecific phenotypic variation in a spatial extent usually not considered in biogeographic studies in the Amazon and demonstrate that both local adaptation and neutral drift are important to explain intraspecific trait diversification at this geographical scale.     

Pleiotropic effects of regulatory variation in tan result in correlation of two pigmentation traits in Drosophila melanogaster

Traits with a common genetic basis frequently display correlated phenotypic responses to selection or environmental conditions. In Drosophila melanogaster, pigmentation of the abdomen and a trident‐shaped region on the thorax are genetically correlated. Here, we used a pooled replicated genomewide association approach (Pool‐GWAS) to identify the genetic basis of variation in thoracic trident pigmentation in two Drosophila melanogaster populations. We confirmed the previously reported large effect of ebony and the association of the cosmopolitan inversion In(3R)Payne. For the first time, we identified tan as another major locus contributing to variation in trident pigmentation. Intriguingly, the regulatory variants of tan that were most strongly associated with female abdominal pigmentation also showed a strong association with trident pigmentation. We validated this common genetic basis in transgenic assays and found qualitatively similar effects on trident and abdominal pigmentation. Further work is required to determine whether this genetic correlation is favoured by natural selection or reflects a neutral by‐product of a shared regulatory architecture.     

Extensive mtDNA variation within the yellow-pine chipmunk,Tamias amoenus (Rodentia: Sciuridae), and phylogeographic inferences for northwest North America

The yellow-pine chipmunk, Tamias amoenus, is common in xerophytic forests throughout much of northwest North America. We analyzed cytochrome b sequence variation from 155 individuals representing 57 localities across the distribution of T. amoenus including 10 additional species of Tamias. Maximum likelihood and parsimony tree estimation methods were used in conjunction with nested clade analysis to infer both deep and population-level processes. Our results indicate that two currently recognized subspecies of T. amoenus (T. a. canicaudus and T. a. cratericus) are not nested within other samples of T. amoenus. Maximum uncorrected levels of intraspecific sequence divergence within remaining samples of T. amoenus are >7%. Substantial geographic variation is characterized by 12 well-supported clades that correspond to distinct mountain ranges, but do not necessarily follow existing subspecific taxonomy. Significant association between geography and genealogy was detected within many of these clades and can be attributed to different population-level processes including past fragmentation, recent range expansion, and isolation by distance.     

Impact of body melanisation on desiccation resistance in montane populations of D. melanogaster: Analysis of seasonal variation

In the montane localities of subtropical regions, winter is the dry season and ectothermic drosophilids are expected to evolve desiccation resistance to cope with drier climatic conditions. An analysis of six montane populations (600–2226 m) of D. melanogaster showed variations for body melanisation (i.e. pigmentation) and desiccation resistance across seasons as well as along altitude. During winter season, plastic changes for melanisation of three posterior abdominal segments (5th + 6th + 7th) correspond with higher desiccation resistance. Thus, we analyzed genetic and plastic effects for these ecophysiological traits by comparing wild-caught and laboratory reared individuals of D. melanogaster for autumn as well as winter season. A ratio of slope values in wild vs. laboratory populations has shown a 1.64-fold plastic effect during autumn; and a two-fold effect during winter. For body melanisation and desiccation resistance, evolutionary response to altitudinal environmental gradient is similar to the phenotypic response across seasons. Thus, our observations are in agreement with the co-gradient hypothesis. Further, we tested the hypothesis whether a thicker cuticle (either due to melanisation or cuticular lipids) leads to lesser cuticular water loss and higher desiccation resistance across seasons as well as according to altitude. Based on within and between population analyses, body melanisation was found to be positively correlated with desiccation resistance but negatively with cuticular water loss. Interestingly, there were no changes in the amount of cuticular lipids per fly across seasons as well as along altitude; and therefore cuticular lipids did not account for desiccation resistance. Cuticular water loss exhibited negative correlation with body melanisation but not with cuticular lipids as well as with changes in body size across seasons. Thus, our data suggest that seasonal changes in body melanisation confer desiccation resistance in montane populations of D. melanogaster.     

Geographic variation in hybridization across a reinforcement contact zone of chorus frogs (Pseudacris)

Reinforcement contact zones, which are secondary contact zones where species are diverging in reproductive behaviors due to selection against hybridization, represent natural laboratories for studying speciation‐in‐action. Here, we examined replicate localities across the entire reinforcement contact zone between North American chorus frogs Pseudacris feriarum and P. nigrita to investigate geographic variation in hybridization frequencies and to assess whether reinforcement may have contributed to increased genetic divergence within species. Previous work indicated these species have undergone reproductive character displacement (RCD) in male acoustic signals and female preferences due to reinforcement. We also examined acoustic signal variation across the contact zone to assess whether signal characteristics reliably predict hybrid index and to elucidate whether the degree of RCD predicts hybridization rate. Using microsatellites, mitochondrial sequences, and acoustic signal information from >1,000 individuals across >50 localities and ten sympatric focal regions, we demonstrate: (1) hybridization occurs and (2) varies substantially across the geographic range of the contact zone, (3) hybridization is asymmetric and in the direction predicted from observed patterns of asymmetric RCD, (4) in one species, genetic distance is higher between conspecific localities where one or both have been reinforced than between nonreinforced localities, after controlling for geographic distance, (5) acoustic signal characters strongly predict hybrid index, and (6) the degree of RCD does not strongly predict admixture levels. By showing that hybridization occurs in all sympatric localities, this study provides the fifth and final line of evidence that reproductive character displacement is due to reinforcement in the chorus frog contact zone. Furthermore, this work suggests that the dual action of cascade reinforcement and partial geographic isolation is promoting genetic diversification within one of the reinforced species.     

Quantitative trait loci responsible for variation in sexually dimorphic traits in Drosophila melanogaster

To understand the mechanisms of morphological evolution and species divergence, it is essential to elucidate the genetic basis of variation in natural populations. Sexually dimorphic characters, which evolve rapidly both within and among species, present attractive models for addressing these questions. In this report, we map quantitative trait loci (QTL) responsible for variation in sexually dimorphic traits (abdominal pigmentation and the number of ventral abdominal bristles and sex comb teeth) in a natural population of Drosophila melanogaster. To capture the pattern of genetic variation present in the wild, a panel of recombinant inbred lines was created from two heterozygous flies taken directly from nature. High-resolution mapping was made possible by cytological markers at the average density of one per 2 cM. We have used a new Bayesian algorithm that allows QTL mapping based on all markers simultaneously. With this approach, we were able to detect small-effect QTL that were not evident in single-marker analyses. Our results show that at least for some sexually dimorphic traits, a small number of QTL account for the majority of genetic variation. The three strongest QTL account for >60% of variation in the number of ventral abdominal bristles. Strikingly, a single QTL accounts for almost 60% of variation in female abdominal pigmentation. This QTL maps to the chromosomal region that Robertson et al. have found to affect female abdominal pigmentation in other populations of D. melanogaster. Using quantitative complementation tests, we demonstrate that this QTL is allelic to the bric a brac gene, whose expression has previously been shown to correlate with interspecific differences in pigmentation. Multiple bab alleles that confer distinct phenotypes appear to segregate in natural populations at appreciable frequencies, suggesting that intraspecific and interspecific variation in abdominal pigmentation may share a similar genetic basis.     

Calls,colours, shape,and genes: a multi‐trait approach to the study of geographic variation in the Amazonian frog Allobates femoralis

Evolutionary divergence in behavioural traits related to mating may represent the initial stage of speciation. Direct selective forces are usually invoked to explain divergence in mate‐recognition traits, often neglecting a role for neutral processes or concomitant differentiation in ecological traits. We adopted a multi‐trait approach to obtain a deeper understanding of the mechanisms behind allopatric divergence in the Amazonian frog, Allobates femoralis. We tested the null hypothesis that geographic distance between populations correlates with genetic and phenotypic divergence, and compared divergence between mate‐recognition (acoustic) and ecological (coloration, body‐shape) traits. We quantified geographic variation in 39 phenotypic traits and a mitochondrial DNA marker among 125 individuals representing eight populations. Geographic variation in acoustic traits was pronounced and tracked the spatial genetic variation, which appeared to be neutral. Thus, the evolution of acoustic traits tracked the shared history of the populations, which is unexpected for pan‐Amazonian taxa or for mate‐recognition traits. Divergence in coloration appeared uncorrelated with genetic distance, and might be partly attributed to local selective pressures, and perhaps to Batesian mimicry. Divergence in body‐shape traits was low. The results obtained depict a complex evolutionary scenario and emphasize the importance of considering multiple traits when disentangling the forces behind allopatric divergence. ©2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 98 , 826–838.     

Geographic variation in advertisement calls in a tree frog species: gene flow and selection hypotheses

Background

In a species with a large distribution relative to its dispersal capacity, geographic variation in traits may be explained by gene flow, selection, or the combined effects of both. Studies of genetic diversity using neutral molecular markers show that patterns of isolation by distance (IBD) or barrier effect may be evident for geographic variation at the molecular level in amphibian species. However, selective factors such as habitat, predator, or interspecific interactions may be critical for geographic variation in sexual traits. We studied geographic variation in advertisement calls in the tree frog Hyla japonica to understand patterns of variation in these traits across Korea and provide clues about the underlying forces for variation.

Methodology

We recorded calls of H. japonica in three breeding seasons from 17 localities including localities in remote Jeju Island. Call characters analyzed were note repetition rate (NRR), note duration (ND), and dominant frequency (DF), along with snout-to-vent length.

Results

The findings of a barrier effect on DF and a longitudinal variation in NRR seemed to suggest that an open sea between the mainland and Jeju Island and mountain ranges dominated by the north-south Taebaek Mountains were related to geographic variation in call characters. Furthermore, there was a pattern of IBD in mitochondrial DNA sequences. However, no comparable pattern of IBD was found between geographic distance and call characters. We also failed to detect any effects of habitat or interspecific interaction on call characters.

Conclusions

Geographic variations in call characters as well as mitochondrial DNA sequences were largely stratified by geographic factors such as distance and barriers in Korean populations of H. japoinca. Although we did not detect effects of habitat or interspecific interaction, some other selective factors such as sexual selection might still be operating on call characters in conjunction with restricted gene flow.     

Evolution of abdominal pigmentation differences across species in the Drosophila dunni subgroup

The Drosophila dunni subgroup displays a nearly perfect latitudinal cline in abdominal pigmentation that likely resulted from selective forces acting in the habitat of each species during speciation. Here we characterize the nature of this clinal variation by developing a quantitative measure to assess variation in abdominal pigmentation within and between the D. dunni subgroup species. Using discriminant analysis, we confirm the existence of a cline and find that our quantitative measure of pigmentation distinguishes each of the species with singular efficacy. We then combine our quantitative phenotypic analysis of pigmentation with the phylogeny of the D. dunni subgroup species and map the species relationships into the three-dimensional morphological space defined by our pigmentation measures. In this manner, we can visualize how the species have traversed the morphological pigmentation space during the course of speciation. Our analysis reveals that natural selection has caused overall intensity of pigmentation among the northernmost species of the cline to converge. Along with this convergence in phenotype has been a relaxation in expression of sexual dimorphism in these species, indicating a possible shift in the relative intensity of natural and sexual selection. Our analysis indicates an accelerated rate of change in pigmentation for the darkest species in addition to this species evolving a novel abdominal pigmentation trait.     

GROWTH TEMPERATURE AND ADULT PIGMENTATION IN TWO DROSOPHILA SIBLING SPECIES: AN ADAPTIVE CONVERGENCE OF REACTION NORMS IN SYMPATRIC POPULATIONS?

Phenotypic plasticity of abdomen pigmentation was investigated in populations of the sibling species Drosophila melanogaster and D. simulans, living in sympatry in two French localities. Ten isofemale lines of each population and species were grown at different constant temperatures spanning their complete thermal range from 12 to 31°C. Genetic variability between isofemale lines was not affected by growth temperature, but was consistently less in D. simulans. For all traits, the dark pigmentation of the abdominal segments decreased according to growth temperature, in agreement with the thermal budget adaptive hypothesis. The shapes of the response curves were different between the abdominal segments, but for a given segment, quite similar in the two species. On average D. simulans was lighter than D. melanogaster, but the difference was mainly expressed at higher temperatures. An interesting result was the difference observed between the two localities: flies from the colder locality (Villeurbanne) were found to be darker than flies from the warmer locality (Bordeaux). Interestingly, this difference was expressed only at low temperatures, 21°C and below, that is, at temperatures encountered in natural conditions. This suggests an adaptive response resulting in a change of the shape of reaction norm and involving genotype-environment interactions. When comparing the genetic structure of geographic populations for quantitative traits, several laboratory environments should be preferred to a single one.     

Trait‐specific consequences of inbreeding on adaptive phenotypic plasticity

Environmental changes may stress organisms and stimulate an adaptive phenotypic response. Effects of inbreeding often interact with the environment and can decrease fitness of inbred individuals exposed to stress more so than that of outbred individuals. Such an interaction may stem from a reduced ability of inbred individuals to respond plastically to environmental stress; however, this hypothesis has rarely been tested. In this study, we mimicked the genetic constitution of natural inbred populations by rearing replicate Drosophila melanogaster populations for 25 generations at a reduced population size (10 individuals). The replicate inbred populations, as well as control populations reared at a population size of 500, were exposed to a benign developmental temperature and two developmental temperatures at the lower and upper margins of their viable range. Flies developed at the three temperatures were assessed for traits known to vary across temperatures, namely abdominal pigmentation, wing size, and wing shape. We found no significant difference in phenotypic plasticity in pigmentation or in wing size between inbred and control populations, but a significantly higher plasticity in wing shape across temperatures in inbred compared to control populations. Given that the norms of reaction for the noninbred control populations are adaptive, we conclude that a reduced ability to induce an adaptive phenotypic response to temperature changes is not a general consequence of inbreeding and thus not a general explanation of inbreeding–environment interaction effects on fitness components.     

云南引种印楝实生种群的表型变异

为了揭示印楝(Azadirachta indica)实生种群表型变异程度和变异规律, 以云南引种印楝人工林为研究对象, 基于9个种群90个单株14个表型性状严格细致的测量, 采用单因素方差分析、巢式方差分析、相关分析、协方差主成分分析(S法)和非加权配对算术平均法(UPGMA)聚类分析等数理方法, 分析了种群的表型变异。结果表明: 印楝种内表型性状在种群间和种群内均存在着较丰富的差异, 种群内的变异大于种群间的变异, 种群间的分化相对较小。对表型性状进行的变异系数多重比较和协方差主成分分析(S法)均显示, 结实和种子化学成分相关性状的变异是造成印楝表型变异的主要来源。利用种群间欧氏距离进行的UPGMA聚类分析结果进一步表明, 印楝9个种群可以分为4类, 表型性状并没有严格依地理距离而聚类。研究结果为印楝的遗传改良工作奠定了基础, 为制定育种策略和人工经营对策提供了科学依据。     

Biogeography of parasite abundance: latitudinal gradient and distance decay of similarity in the abundance of fleas and mites,parasitic on small mammals in the Palearctic,at three spatial scales

We tested whether biogeographic patterns characteristic for biological communities can also apply to populations and investigated geographic patterns of variation in abundance of ectoparasites (fleas and mites) collected from bodies of their small mammalian hosts (rodents and shrews) in the Palearctic at continental, regional and local scales. We asked whether (i) there is a relationship between latitude and abundance and (ii) similarity in abundance follows a distance decay pattern or it is better explained by variation in extrinsic biotic and abiotic factors. We analysed the effect of latitude on mean intraspecific abundance using general linear models including proportional abundance of its principal host as an additional predictor variable. Then, we examined the relative effect of geographic distance, biotic and abiotic dissimilarities among regions, subregions or localities on the intraspecific dissimilarity in abundance among regions, subregions or localities using Generalized Dissimilarity Modelling. We found no relationship between latitude and intraspecific flea or mite abundance. In both taxa, environmental dissimilarity explained the largest part of the deviance of spatial variation in abundance, whereas the effect of the dissimilarity in the principal host abundance was of secondary importance and the effect of geographic distance was minor. These patterns were generally consistent across the three spatial scales, although environmental variation and dissimilarity in principal host abundance were equally important at the local scale in fleas but not in mites. We conclude that biogeographic patterns related to latitude and geographic distance do not apply to spatial variation of ectoparasite abundance. Instead, the geographic distribution of abundance in arthropod ectoparasites depends on their responses, mainly to the off-host environment and to a lesser extent the abundance of their principal hosts.     

Genetic and developmental analysis of abdominal pigmentation differences across species in the Drosophila dunni subgroup

Abdominal pigmentation pattern varies dramatically among the species of the Drosophila dunni subgroup across the islands of the Caribbean. Previously, we developed a quantitative measure of abdominal pigmentation to assess phenotypic variation within and between species of this group. In this paper, we use this quantitative measure in an interspecific genetic analysis to decipher the underlying genetic basis of pigmentation differences between one of the lightest and the darkest species in the group. Our analysis shows that pigmentation expression in different areas of the abdomen is under separate genetic control. For these different abdominal regions, we detected a wide range of genetic effects, including X-linked, autosomal additive, near single-gene dominant, and sex-specific effects. Combining these genetic results with our earlier phenotypic and phylogenetic analyses, we present a simple conceptual model to explain how change in the control of expression of pigmentation has evolved throughout the D. dunni subgroup.     

Ecological genetics of adaptive color polymorphism in pocket mice: geographic variation in selected and neutral genes

Patterns of geographic variation in phenotype or genotype may provide evidence for natural selection. Here, we compare phenotypic variation in color, allele frequencies of a pigmentation gene (the melanocortin-1 receptor, Mc1r), and patterns of neutral mitochondrial DNA (mtDNA) variation in rock pocket mice (Chaetodipus intermedius) across a habitat gradient in southern Arizona. Pocket mice inhabiting volcanic lava have dark coats with unbanded, uniformly melanic hairs, whereas mice from nearby light-colored granitic rocks have light coats with banded hairs. This color polymorphism is a presumed adaptation to avoid predation. Previous work has demonstrated that two Mc1r alleles, D and d, differ by four amino acids, and are responsible for the color polymorphism: DD and Dd genotypes are melanic whereas dd genotypes are light colored. To determine the frequency of the two Mc1r allelic classes across the dark-colored lava and neighboring light-colored granite, we sequenced the Mc1r gene in 175 individuals from a 35-km transect in the Pinacate lava region. We also sequenced two neutral mtDNA genes, COIII and ND3, in the same individuals. We found a strong correlation between Mc1r allele frequency and habitat color and no correlation between mtDNA markers and habitat color. Using estimates of migration from mtDNA haplotypes between dark- and light-colored sampling sites and Mc1r allele frequencies at each site, we estimated selection coefficients against mismatched Mc1r alleles, assuming a simple model of migration-selection balance. Habitat-dependent selection appears strong but asymmetric: selection is stronger against light mice on dark rock than against melanic mice on light rock. Together these results suggest that natural selection acts to match pocket mouse coat color to substrate color, despite high levels of gene flow between light and melanic populations.