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.     

EVIDENCE OF ADAPTATION FROM ANCESTRAL VARIATION IN YOUNG POPULATIONS OF BEACH MICE

To understand how organisms adapt to novel habitats, which involves both demographic and selective events, we require knowledge of the evolutionary history of populations and also selected alleles. There are still few cases in which the precise mutations (and hence, defined alleles) that contribute to adaptive change have been identified in nature; one exception is the genetic basis of camouflaging pigmentation of oldfield mice (Peromyscus polionotus) that have colonized the sandy dunes of Florida's Gulf Coast. To quantify the genomic impact of colonization as well as the signature of selection, we resequenced 5000 1.5‐kb noncoding loci as well as a 160‐kb genomic region surrounding the melanocortin‐1 receptor (Mc1r), a gene that contributes to pigmentation differences, in beach and mainland populations. Using a genome‐wide phylogenetic approach, we recovered a single monophyletic group comprised of beach mice, consistent with a single colonization event of the Gulf Coast. We also found evidence of a severe founder event, estimated to have occurred less than 3000 years ago. In this demographic context, we show that all beach subspecies share a single derived light Mc1r allele, which was likely selected from standing genetic variation that originated in the mainland. Surprisingly, we were unable to identify a clear signature of selection in the Mc1r region, despite independent evidence that this locus contributes to adaptive coloration. Nonetheless, these data allow us to reconstruct and compare the evolutionary history of populations and alleles to better understand how adaptive evolution, following the colonization of a novel habitat, proceeds in nature.     

Natural selection along an environmental gradient: a classic cline in mouse pigmentation

We revisited a classic study of morphological variation in the oldfield mouse ( Peromyscus polionotus ) to estimate the strength of selection acting on pigmentation patterns and to identify the underlying genes. We measured 215 specimens collected by Francis Sumner in the 1920s from eight populations across a 155-km, environmentally variable transect from the white sands of Florida's Gulf coast to the dark, loamy soil of southeastern Alabama. Like Sumner, we found significant variation among populations: mice inhabiting coastal sand dunes had larger feet, longer tails, and lighter pigmentation than inland populations. Most striking, all seven pigmentation traits examined showed a sharp decrease in reflectance about 55 km from the coast, with most of the phenotypic change occurring over less than 10 km. The largest change in soil reflectance occurred just south of this break in pigmentation. Geographic analysis of microsatellite markers shows little interpopulation differentiation, so the abrupt change in pigmentation is not associated with recent secondary contact or reduced gene flow between adjacent populations. Using these genetic data, we estimated that the strength of selection needed to maintain the observed distribution of pigment traits ranged from 0.0004 to 21%, depending on the trait and model used. We also examined changes in allele frequency of SNPs in two pigmentation genes, Mc1r and Agouti , and show that mutations in the cis -regulatory region of Agouti may contribute to this cline in pigmentation. The concordance between environmental variation and pigmentation in the face of high levels of interpopulation gene flow strongly implies that natural selection is maintaining a steep cline in pigmentation and the genes underlying it.     

An overview of the relationships between mimicry and crypsis

There have been many different and conflicting definitions of mimicry. Some of the definitions of mimicry include crypsis and others do not. Each definition includes different groups of phenomena and uses different criteria to distinguish mimetic from non-mimetic phenomena. The confusion is eliminated by a consideration of the criteria of all definitions. This shows that there are in fact three major criteria dividing six phenomona, rather than a single dichotomy between mimicry and crypsis (Table 2). The criteria are defined by the results of a mistake in discrimination between the model and mimìc: (a) the mistake does or does not depend upon relationship between mimic and background; (b) the mistake has or has no effect on the population dynamics or evolution of the model and (c) the mistake affects dynamics or evolution of one or of many models. The main reason for the contusion about mimicry and crypsis is that each author's definition includes differing and partially overlapping subsets of the six classes: crypsis; masquerade; Batesism; Müllerism; polymorphism and convergence.     

Genetic evidence for the convergent evolution of light skin in Europeans and East Asians

Human skin pigmentation shows a strong positive correlation with ultraviolet radiation intensity, suggesting that variation in skin color is, at least partially, due to adaptation via natural selection. We investigated the evolution of pigmentation variation by testing for the presence of positive directional selection in 6 pigmentation genes using an empirical F(ST) approach, through an examination of global diversity patterns of these genes in the Centre d'Etude du Polymorphisme Humain (CEPH)-Diversity Panel, and by exploring signatures of selection in data from the International HapMap project. Additionally, we demonstrated a role for MATP in determining normal skin pigmentation variation using admixture mapping methods. Taken together (with the results of previous admixture mapping studies), these results point to the importance of several genes in shaping the pigmentation phenotype and a complex evolutionary history involving strong selection. Polymorphisms in 2 genes, ASIP and OCA2, may play a shared role in shaping light and dark pigmentation across the globe, whereas SLC24A5, MATP, and TYR have a predominant role in the evolution of light skin in Europeans but not in East Asians. These findings support a case for the recent convergent evolution of a lighter pigmentation phenotype in Europeans and East Asians.     

An overview of the relationships between mimicry and crypsis

There have been many different and conflicting definitions of mimicry. Some of the definitions of mimicry include crypsis and others do not. Each definition includes different groups of phenomena and uses different criteria to distinguish mimetic from non-mimetic phenomena. The confusion is eliminated by a consideration of the criteria of all definitions. This shows that there are in fact three major criteria dividing six phenomona, rather than a single dichotomy between mimicry and crypsis (Table 2). The criteria are defined by the results of a mistake in discrimination between the model and mimìc: (a) the mistake does or does not depend upon relationship between mimic and background; (b) the mistake has or has no effect on the population dynamics or evolution of the model and (c) the mistake affects dynamics or evolution of one or of many models. The main reason for the contusion about mimicry and crypsis is that each author's definition includes differing and partially overlapping subsets of the six classes: crypsis; masquerade; Batesism; Müllerism; polymorphism and convergence.     

Gene flow and melanism in Lake Erie garter snake populations

Melanistic garter snakes ( Thamnophis sirtalis ) are unusually common near Lake Erie, apparently because selection for thermoregulatory ability in cool lake-shore habitats (which favours melanistic morphs) outweighs selection for crypsis (which favours striped morphs). However, morph frequencies are highly variable among sites, suggesting that random genetic drift also influences colour pattern. In an effort to better understand the evolutionary processes influencing garter snake colour patterns, we estimated Fx and Nm (the number of migrants per generation) among island and mainland populations from patterns of allozymic variation detected using electrophoresis. Estimates of Nm were high, ranging from 2.7 to 37.6 between pairs of study sites and making it unlikely that differences in morph frequencies among sites were solely the result of random genetic drift. Furthermore, differences in F _st estimates between colour pattern (a one-locus two-allele trait) and allozyme loci suggest that colour pattern alleles are not in Hardy-Weinberg equilibrium, most likely as a result of natural selection. Comparison of allozymic data from Lake Erie with those from more distant sites suggests that gene flow occurs over long distances in T. sirtalis.     

Testing the adaptive hypothesis of Batesian mimicry among hybridizing North American admiral butterflies

Batesian mimicry is characterized by phenotypic convergence between an unpalatable model and a palatable mimic. However, because convergent evolution may arise via alternative evolutionary mechanisms, putative examples of Batesian mimicry must be rigorously tested. Here, we used artificial butterfly facsimiles (N = 4000) to test the prediction that (1) palatable Limenitis lorquini butterflies should experience reduced predation when in sympatry with their putative model, Adelpha californica, (2) protection from predation on L. lorquini should erode outside of the geographical range of the model, and (3) mimetic color pattern traits are more variable in allopatry, consistent with relaxed selection for mimicry. We find support for these predictions, implying that this convergence is the result of selection for Batesian mimicry. Additionally, we conducted mark–recapture studies to examine the effect of mimicry and found that mimics survive significantly longer at sites where the model is abundant. Finally, in contrast to theoretical predictions, we found evidence that the Batesian model (A. californica) is protected from predation outside of its geographic range. We discuss these results considering the ongoing hybridization between L. lorquini and its sister species, L. weidemeyerii, and growing evidence that selection for mimicry predictably leads to a reduction in gene flow between nascent species.     

On transition bias in mitochondrial genes of pocket gophers

The relative contribution of mutation and purifying selection to transition bias has not been quantitatively assessed in mitochondrial protein genes. The observed transition/transversion (s/v) ratio is (μ _s P _s)/(μ _v P _v), where μ _s and μ _v denote mutation rate of transitions and transversions, respectively, andP _s andP _v denote fixation probabilities of transitions and transversions, respectively. Because selection against synonymous transitions can be assumed to be roughly equal to that against synonymous transversions,P _s/P_v ≈ 1 at fourfold degenerate sites, so that thes/v ratio at fourfold degenerate sites is approximately μ _s /μ _v , which is a measure of mutational contribution to transition bias. Similarly, thes/v ratio at nondegenerate sites is also an estimate of μ _s /μ _v if we assume that selection against nonsynonymous transitions is roughly equal to that against nonsynonymous transversions. In two mitochondrial genes, cytochrome oxidase subunit I (COI) and cytochromeb (cyt-b) in pocket gophers, thes/v ratio is about two at nondegenerate and fourfold degenerate sites for both the COI and the cyt-b genes. This implies that mutation contribution to transition bias is relatively small. In contrast, thes/v ratio is much greater at twofold degenerate sites, being 48 for COI and 40 for cyt-b. Given that the μ _s /μ _v ratio is about 2, theP _s/P_v ratio at twofold degenerate sites must be on the order of 20 or greater. This suggests a great effect of purifying selection on transition bias in mitochondrial protein genes because transitions are synonymous and transversions are nonsynonymous at twofold degenerate sites in mammalian mitochondrial genes. We also found that nonsynonymous mutations at twofold degenerate sites are more neutral than nonsynonymous mutations at nondegenerate sites, and that the COI gene is subject to stronger purifying selection than is the cyt-b gene. A model is presented to integrate the effect of purifying selection, codon bias, DNA repair and GC content ons/v ratio of protein-coding genes. Correspondence to: X. Xia     

Genetic diversity and differentiation of Kermode bear populations

The Kermode bear is a white phase of the North American black bear that occurs in low to moderate frequency on British Columbia's mid-coast. To investigate the genetic uniqueness of populations containing the white phase, and to ascertain levels of gene flow among populations, we surveyed 10 highly polymorphic microsatellite loci, assayed from trapped bear hairs. A total of 216 unique bear genotypes, 18 of which were white, was sampled among 12 localities. Island populations, where Kermodes are most frequent, show approximately 4% less diversity than mainland populations, and the island richest in white bears (Gribbell) exhibited substantial genetic isolation, with a mean pairwise FST of 0.14 with other localities. Among all localities, FST for the molecular variant underlying the coat-colour difference (A893G) was 0.223, which falls into the 95th percentile of the distribution of FST values among microsatellite alleles, suggestive of greater differentiation for coat colour than expected under neutrality. Control-region sequences confirm that Kermode bears are part of a coastal or western lineage of black bears whose existence predates the Wisconsin glaciation, but microsatellite variation gave no evidence of past population expansion. We conclude that Kermodism was established and is maintained in populations by a combination of genetic isolation and somewhat reduced population sizes in insular habitat, with the possible contribution of selective pressure and/or nonrandom mating.     

The genetic basis of adaptation: lessons from concealing coloration in pocket mice

Recent studies on the genetics of adaptive coat-color variation in pocket mice (Chaetodipus intermedius) are reviewed in the context of several on-going debates about the genetics of adaptation. Association mapping with candidate genes was used to identify mutations responsible for melanism in four different populations of C. intermedius. Here, I review four main results (i) a single gene, the melanocortin-1-receptor (Mc1r), appears to be responsible for most of the phenotypic variation in color in one population, the Pinacate site; (ii) four or fewer nucleotide changes at Mc1r appear to be responsible for the difference in receptor function; (iii) studies of migration-selection balance suggest that the selection coefficient associated with the dark Mc1r allele at the Pinacate site is large; and (iv) different (unknown) genes underlie the evolution of melanism on three other lava flows. These findings are discussed in light of the evolution of convergent phenotypes, the average size of phenotypic effects underlying adaptation, the evolution of dominance, and the distinction between adaptations caused by changes in gene dosage versus gene structure.     

The genetic consequences of selection in natural populations

The selection coefficient, s, quantifies the strength of selection acting on a genetic variant. Despite this parameter's central importance to population genetic models, until recently we have known relatively little about the value of s in natural populations. With the development of molecular genetic techniques in the late 20th century and the sequencing technologies that followed, biologists are now able to identify genetic variants and directly relate them to organismal fitness. We reviewed the literature for published estimates of natural selection acting at the genetic level and found over 3000 estimates of selection coefficients from 79 studies. Selection coefficients were roughly exponentially distributed, suggesting that the impact of selection at the genetic level is generally weak but can occasionally be quite strong. We used both nonparametric statistics and formal random‐effects meta‐analysis to determine how selection varies across biological and methodological categories. Selection was stronger when measured over shorter timescales, with the mean magnitude of s greatest for studies that measured selection within a single generation. Our analyses found conflicting trends when considering how selection varies with the genetic scale (e.g., SNPs or haplotypes) at which it is measured, suggesting a need for further research. Besides these quantitative conclusions, we highlight key issues in the calculation, interpretation, and reporting of selection coefficients and provide recommendations for future research.     

A decline of melanism in the peppered moth Biston betularia in The Netherlands

Melanic forms of the peppered moth Biston betularia were well established in The Netherlands by the end of the 19th century, indeed the first records of the black carbonaria form in 1867 are only about 20 years later than in England. Analysis of extensive sampling data collected by B. J. Lempke for a period of several years beginning in 1969 shows that carbonaria was at a frequency of about 60 to 70% in most of the country where epiphyte communities on trees were reduced due to the effects of air pollution. The pale typica and the three intermediate insularia forms were each at similar, low frequencies. Only in the extreme north and south-east of The Netherlands where epiphyte floras were richer was carbonaria at a lower frequency of less than 40%. Samples collected from seven localities in 1988 show that carbonaria has dramatically declined to a frequency of less than 10%. In contrast to England, the fully black form is being replaced not only by typica but also by the darkest of the insularia phenotypes. The decline in melanism coincides with a period of decreasing levels of sulphur dioxide and of increasing species diversity of lichens on trees.     

Skin layer‐specific transcriptional profiles in normal and recessive yellow (Mc1re/Mc1re) mice

The melanocortin 1 receptor (Mc1r) plays a central role in cutaneous biology, but is expressed at very low levels by a small fraction of cells in the skin. In humans, loss‐of‐function MC1R mutations cause fair skin, freckling, red hair, and increased predisposition to melanoma; in mice, Mc1r loss‐of‐function is responsible for the recessive yellow mutation, associated with pheomelanic hair and a decreased number of epidermal melanocytes. To better understand how Mc1r signaling affects different cutaneous phenotypes, we examined large‐scale patterns of gene expression in different skin components (whole epidermal sheets, basal epidermal cells and whole skins) of neonatal (P2.5) normal and recessive yellow mice, starting with a 26K mouse cDNA microarray. From c. 17 000 genes whose levels could be accurately measured in neonatal skin, we identified 883, 2097 and 552 genes that were uniquely expressed in the suprabasal epidermis, basal epidermis and dermis, respectively; specific biologic roles could be assigned for each class. Comparison of normal and recessive yellow mice revealed 69 differentially expressed genes, of which the majority had not been previously implicated in Mc1r signaling. Surprisingly, many of the Mc1r‐dependent genes are expressed in cells other than melanocytes, even though Mc1r expression in the skin is confined almost exclusively to epidermal melanocytes. These results reveal new targets for Mc1r signaling, and point to a previously unappreciated role for a Mc1r‐dependent paracrine effect of melanocytes on other components of the skin.     

Tests of hypotheses on predation as a factor maintaining polymorphic melanism in coastalplain fox squirrels (Sciurus niger L.)

Fox squirrels ( Sciurus niger ) in the south-eastern U.S. coastal plain differ from those in the rest of the species' range by having black heads with white nose and ears. Postcranially on the dorsum, they also show interindividual colour variation, ranging from all-light agouti to all-black non-agouti. I present results of experiments undertaken to test whether the evolution of these pigmentary features can be attributed to interactions with predators. Comparative static crypsis (when the squirrels are still) was tested by determining how well specimens of the squirrel morphs matched their backgrounds in terms of intensity (brightness) distributions and patch size. Comparative dynamic effects (when the squirrels are moving) were tested using captive red-tailed hawks ( Buteo jamaicensis ) and models of the fox squirrel morphs. Results indicated that morphs of coastal-plain fox squirrels with all-light backs are better static matches to unburned backgrounds than are the darker morphs with which they coexist. Against fire-blackened backgrounds, fox squirrel morphs with intermediate and all-black backs are better intensity matches than those with all-light backs, but the dark morphs are not better patch-length matches than light morphs. The superiority of dark morphs in intensity matching for a short time after a fire does not seem sufficient to account for their long-term maintenance. Hawks responded more slowly to moving fox squirrel models with intermediate amounts of black on the back than to all-light or all-dark morphs; this result suggests a possible factor that would favour retention of genes for dorsal blackness in the coastal-plain population. Patchy black and white heads appear to promote static crypsis of south-eastern fox squirrels, and hawks reacted more slowly to moving squirrel models with such colouration than to those with plain heads.     

Rapid evolution of fire melanism in replicated populations of pygmy grasshoppers

Evolutionary theory predicts an interactive process whereby spatiotemporal environmental heterogeneity will maintain genetic variation, while genetic and phenotypic diversity will buffer populations against stress and allow for fast adaptive evolution in rapidly changing environments. Here, we study color polymorphism patterns in pygmy grasshoppers (Tetrix subulata) and show that the frequency of the melanistic (black) color variant was higher in areas that had been ravaged by fires the previous year than in nonburned habitats, that, in burned areas, the frequency of melanistic grasshoppers dropped from ca. 50% one year after a fire to 30% after four years, and that the variation in frequencies of melanistic individuals among and within populations was genetically based on and represented evolutionary modifications. Dark coloration may confer a selective benefit mediated by enhanced camouflage in recently fire-ravaged areas characterized by blackened visual backgrounds before vegetation has recovered. These findings provide rare evidence for unusually large, extremely rapid adaptive contemporary evolution in replicated natural populations in response to divergent and fluctuating selection associated with spatiotemporal environmental changes.     

Population genomics of natural and experimental populations of guppies (Poecilia reticulata)

Convergent evolution represents one of the best lines of evidence for adaptation, but few cases of phenotypic convergence are understood at the genetic level. Guppies inhabiting the Northern Mountain Range of Trinidad provide a classic example of phenotypic convergent evolution, where adaptation to low or high predation environments has been found for a variety of traits. A major advantage of this system is the possibility of long‐term experimental studies in nature, including transplantation from high to low predation sites. We used genome scans of guppies from three natural high and low predation populations and from two experimentally established populations and their sources to examine whether phenotypic convergent evolution leaves footprints at the genome level. We used population‐genetic modelling approaches to reconstruct the demographic history and migration among sampled populations. Naturally colonized low predation populations had signatures of increased effective population size since colonization, while introduction populations had signatures of decreased effective population size. Only a small number of regions across the genome had signatures of selection in all natural populations. However, the two experimental populations shared many genomic regions under apparent selection, more than expected by chance. This overlap coupled with a population decrease since introduction provides evidence for convergent selection occurring in the two introduced populations. The lack of genetic convergence in the natural populations suggests that convergent evolution is lacking in these populations or that the effects of selection become difficult to detect after a long‐time period.     

Changes in C uptake in populations of Chlamydomonas reinhardtii selected at high CO2

Estimates of the effect of increased global atmospheric CO(2) levels on oceanic primary productivity depend on the physiological responses of contemporary phytoplankton populations. However, microalgal populations will possibly adapt to rising CO(2) levels in such a way that they become genetically different from contemporary populations. The unknown properties of these future populations introduce an undefined error into predictions of C pool dynamics, especially the presence and size of the biological C pump. To address the bias in predictions introduced by evolution, we measured the kinetics of CO(2) uptake in populations of Chlamydomonas reinhardtii that had been selected for growth at high CO(2) for 1000 generations. Following selection at high CO(2), the populations were unable to induce high-affinity CO(2) uptake, and one line had a lower rate of net CO(2) uptake. We attribute this to conditionally neutral mutations in genes affecting the C concentrating mechanism (CCM). Lower affinity CO(2) uptake, in addition to smaller population sizes, results in a significant reduction in net CO(2) uptake of about 38% relative to contemporary populations under the same conditions. This shows how predictions about the properties of communities in the future can be influenced by the effect of natural selection.     

Evolution of novel genes

Much progress in understanding the evolution of new genes has been accomplished in the past few years. Molecular mechanisms such as illegitimate recombination and LINE element mediated 3' transduction underlying exon shuffling, a major process for generating new genes, are better understood. The identification of young genes in invertebrates and vertebrates has revealed a significant role of adaptive evolution acting on initially rudimentary gene structures created as if by evolutionary tinkers. New genes in humans and our primate relatives add a new component to the understanding of genetic divergence between humans and non-humans.