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.     

Adaptive variation in beach mice produced by two interacting pigmentation genes

Little is known about the genetic basis of ecologically important morphological variation such as the diverse color patterns of mammals. Here we identify genetic changes contributing to an adaptive difference in color pattern between two subspecies of oldfield mice (Peromyscus polionotus). One mainland subspecies has a cryptic dark brown dorsal coat, while a younger beach-dwelling subspecies has a lighter coat produced by natural selection for camouflage on pale coastal sand dunes. Using genome-wide linkage mapping, we identified three chromosomal regions (two of major and one of minor effect) associated with differences in pigmentation traits. Two candidate genes, the melanocortin-1 receptor (Mc1r) and its antagonist, the Agouti signaling protein (Agouti), map to independent regions that together are responsible for most of the difference in pigmentation between subspecies. A derived mutation in the coding region of Mc1r, rather than change in its expression level, contributes to light pigmentation. Conversely, beach mice have a derived increase in Agouti mRNA expression but no changes in protein sequence. These two genes also interact epistatically: the phenotypic effects of Mc1r are visible only in genetic backgrounds containing the derived Agouti allele. These results demonstrate that cryptic coloration can be based largely on a few interacting genes of major effect.     

Exclusion of melanocortin-1 receptor (mc1r) and agouti as candidates for dominant black in dogs

The domestic dog exhibits a variety of coat colors that encompass a wide range of variation among different breeds. Very little is known about the molecular biology of dog pigmentation; current understanding is based mostly on traditional breeding experiments, which in some cases have suggested genetic interactions that are different from those reported in other mammals. We have examined the molecular genetics of dominant black, a uniform coat color characteristic of black Labrador retrievers or Newfoundlands that has been proposed to be caused by either variation in the melanocortin-1 receptor gene (Mc1r) or by variation in the Agouti gene (A). We identified several coding polymorphisms within Mc1r and several simple sequence repeat polymorphisms closely linked to A, and examined their inheritance in a Labrador retriever x greyhound cross that segregates dominant black. No single Mc1r allele was found consistently in animals carrying dominant black, and neither Mc1r nor A cosegregated with dominant black. These results refine our understanding of mammalian coat color inheritance and suggest that dominant black coat color in dogs is caused by a gene not previously implicated in pigment type switching.     

Sequence variation in the melanocortin-1 receptor (MC1R) pigmentation gene and its role in the cryptic coloration of two South American sand lizards

In reptiles, dorsal body darkness often varies with substrate color or temperature environment, and is generally presumed to be an adaptation for crypsis or thermoregulation. However, the genetic basis of pigmentation is poorly known in this group. In this study we analyzed the coding region of the melanocortin-1-receptor (MC1R) gene, and therefore its role underlying the dorsal color variation in two sympatric species of sand lizards (Liolaemus) that inhabit the southeastern coast of South America: L. occipitalis and L. arambarensis. The first is light-colored and occupies aeolic pale sand dunes, while the second is brownish and lives in a darker sandy habitat. We sequenced 630 base pairs of MC1R in both species. In total, 12 nucleotide polymorphisms were observed, and four amino acid replacement sites, but none of them could be associated with a color pattern. Comparative analysis indicated that these taxa are monomorphic for amino acid sites that were previously identified as functionally important in other reptiles. Thus, our results indicate that MC1R is not involved in the pigmentation pattern observed in Liolaemus lizards. Therefore, structural differences in other genes, such as ASIP, or variation in regulatory regions of MC1R may be responsible for this variation. Alternatively, the phenotypic differences observed might be a consequence of non-genetic factors, such as thermoregulatory mechanisms.     

Variation of the melanocortin 1 receptor gene in the macaques

Melanocortin 1 receptor (MC1R), a G-coupled seven-transmembrane receptor protein, plays a key role in the regulation of melanin synthesis in mammals. Sequence variation of the MC1R gene (MC1R) has been associated with pigmentation phenotypes in humans and in several animal species. The macaques (genus Macaca) are known to show a marked inter-specific variation in coat color although the causative genetic variation remains unclear. We investigated nucleotide sequences of the MC1R in 67 individuals of 18 macaque species with different coat color phenotypes including black and agouti. Twenty-eight amino acid replacements were identified in the macaques, but none of these amino acid replacements could explain the black coat color of Macaca silenus and the Sulawesi macaque species. Our molecular evolutionary analysis has revealed that nonsynonymous substitution/synonymous substitution (dN/dS) ratio of the MC1R has not been uniform in the macaque groups and, moreover, their coat color and dN/dS ratio were not related. These results suggest that the MC1R is unlikely to be responsible for the coat color variation of the macaques and functions of MC1R other than pigmentation might be associated with the different selective pressures on the MC1R in macaques.     

Sequence variation in the Mc1r gene for a group of polymorphic snakes

Studying the genetic factors underlying phenotypic traits can provide insight into dynamics of selection and molecular basis of adaptation, but this goal can be difficult for non-model organisms without extensive genomic resources. However, sequencing candidate genes for the trait of interest can facilitate the study of evolutionary genetics in natural populations. We sequenced the melanocortin-1 receptor (Mc1r) to study the genetic basis of color polymorphism in a group of snake species with variable black banding, the genera Sonora, Chilomeniscus, and Chionactis. Mc1r is an important gene in the melanin synthesis pathway and is associated with ecologically important variation in color pattern in birds, mammals, and other squamate reptiles. We found that Mc1r nucleotide sequence was variable and that within our focal Sonora species, there are both fixed and heterozygous nucleotide substitutions that result in an amino acid change and selection analyses indicated that Mc1r sequence was likely under purifying selection. However, we did not detect any statistical association with the presence or absence of black bands. Our results agree with other studies that have found no role for sequence variation in Mc1r and highlight the importance of comparative data for studying the phenotypic associations of candidate genes.     

Evolution of a pigmentation gene,the melanocortin-1 receptor,in primates

The melanocortin-1 receptor (MC1R) forms a critical switch in the production of orange/red pheomelanin and black/brown eumelanin pigments during hair development in mammals. The molecular evolution of the melanocortin-1 receptor gene was investigated in a broad range of primate species, including several groups with large differences in distribution of orange/red and black hairs. Primate MC1R has been subject to purifying selection throughout most of its evolution, with small changes in selective constraint being detected early in primate evolution. In contrast to the situation in humans and domestic mammals, many intraspecific and intrageneric differences in primate coat color cannot be attributed to changes in the MC1R coding sequence. Nevertheless, important changes in the biochemical function of MC1R are suggested by mutations in sites of known functional importance, particularly in New World monkeys and lemurs. The evolution of the MC1R in lion tamarins is anomalous, with a combination of a high nonsynonymous to synonymous substitution rate (dN/dS) ratio, deletions, and substitutions.     

Genetic, Comparative Genomic, and Expression Analyses of the Mc1r Locus in the Polychromatic Midas Cichlid Fish (Teleostei, Cichlidae Amphilophus sp.) Species Group

Natural populations of the Midas cichlid species in several different crater lakes in Nicaragua exhibit a conspicuous color polymorphism. Most individuals are dark and the remaining have a gold coloration. The color morphs mate assortatively and sympatric population differentiation has been shown based on neutral molecular data. We investigated the color polymorphism using segregation analysis and a candidate gene approach. The segregation patterns observed in a mapping cross between a gold and a dark individual were consistent with a single dominant gene as a cause of the gold phenotype. This suggests that a simple genetic architecture underlies some of the speciation events in the Midas cichlids. We compared the expression levels of several candidate color genes Mc1r, Ednrb1, Slc45a2, and Tfap1a between the color morphs. Mc1r was found to be up regulated in the gold morph. Given its widespread association in color evolution and role on melanin synthesis, the Mc1r locus was further investigated using sequences derived from a genomic library. Comparative analysis revealed conserved synteny in relation to the majority of teleosts and highlighted several previously unidentified conserved non-coding elements (CNEs) in the upstream and downstream regions in the vicinity of Mc1r. The identification of the CNEs regions allowed the comparison of sequences from gold and dark specimens of natural populations. No polymorphisms were found between in the population sample and Mc1r showed no linkage to the gold phenotype in the mapping cross, demonstrating that it is not causally related to the color polymorphism in the Midas cichlid.     

Sequence analysis of three pigmentation genes in the Newfoundland population of Canis latrans links the Golden Retriever Mc1r variant to white coat color in coyotes

Three genes, Mc1r, Agouti, and CBD103, interact in a type-switching process that controls much of the pigmentation variation observed in mammals. A deletion in the CBD103 gene is responsible for dominant black color in dogs, while the white-phased black bear (“spirit bear”) of British Columbia, Canada, is the lightest documented color variant caused by a mutation in Mc1r. Rare all-white animals have recently been discovered in a new northeastern population of the coyote in insular Newfoundland and Labrador, Canada. To investigate the causative gene and mutation of white coat in coyotes, we sequenced the three type-switching genes in white and dark-phased animals from Newfoundland. The only sequence variants unambiguously associated with white color were in Mc1r, and one of these variants causes the amino acid variant R306Ter, a premature stop codon also linked to coat color in Golden Retrievers and other dogs with yellow/red coats. The allele carrying R306Ter in coyotes matches that in the Golden Retriever at other variable amino acid sites and hence may have originated in these dogs. Coyotes experienced introgression with wolves and dogs as they colonized northeastern North America, and coyote/Golden Retriever interactions have been observed in Newfoundland. We speculate that natural selection, with or without a founder effect, may contribute to the observed frequency of white coyotes in Newfoundland, as it has contributed to the high frequency of white bears, and of a domestic dog-derived CBD allele in gray wolves.     

Loss‐of‐function mutations in the melanocortin 1 receptor cause disruption of dorso‐ventral countershading in teleost fish

The melanocortin 1 receptor (MC1R) is the central melanocortin receptor involved in vertebrate pigmentation. Mutations in this gene cause variations in coat coloration in amniotes. Additionally, in mammals MC1R is the main receptor for agouti‐signaling protein (ASIP), making it the critical receptor for the establishment of dorsal‐ventral countershading. In fish, Mc1r is also involved in pigmentation, but it has been almost exclusively studied in relation to melanosome dispersion activity and as a putative genetic factor involved in dark/light adaptation. However, its role as the crucial component for the Asip1‐dependent control of dorsal‐ventral pigmentation remains unexplored. Using CRISPR/Cas9, we created mc1r homozygous knockout zebrafish and found that loss‐of‐function of mc1r causes a reduction of countershading and a general paling of the animals. We find ectopic development of melanophores and xanthophores, accompanied by a decrease in iridophore numbers in the ventral region of mc1r mutants. We also reveal subtle differences in the role of mc1r in repressing pigment cell development between the skin and scale niches in ventral regions.     

Polymorphism of the Goat Agouti Signaling Protein Gene and Its Relationship with Coat Color in Italian and Spanish Breeds

Agouti signaling protein (ASIP) is one of the key players in the modulation of hair pigmentation in mammals. Binding to the melanocortin 1 receptor, ASIP induces the synthesis of phaeomelanin, associated with reddish brown, red, tan, and yellow coats. We have sequenced 2.8?kb of the goat ASIP gene in 48 individuals and identified two missense (Cys126Gly and Val128Gly) and two intronic polymorphisms. In silico analysis revealed that the Cys126Gly substitution may cause a structural change by disrupting a highly conserved disulfide bond. We studied its segregation in 12 Spanish and Italian goat breeds (N?=?360) with different pigmentation patterns and found striking differences in the frequency of the putative loss-of-function Gly(126) allele (Italian 0.43, Spanish Peninsular 0.08), but we did not observe a clear association with coat color. This suggests that the frequency of this putative loss-of-function allele has evolved under the influence of demographic rather than selection factors in goats from these two geographical areas.     

MC1R-dependent, melanin-based colour polymorphism is associated with cell-mediated response in the Eleonora's falcon

Colour polymorphism in vertebrates is usually under genetic control and may be associated with variation in physiological traits. The melanocortin 1 receptor (Mc1r) has been involved repeatedly in melanin-based pigmentation but it was thought to have few other physiological effects. However, recent pharmacological studies suggest that MC1R could regulate the aspects of immunity. We investigated whether variation at Mc1r underpins plumage colouration in the Eleonora's falcon. We also examined whether nestlings of the different morphs differed in their inflammatory response induced by phytohemagglutinin (PHA). Variation in colouration was due to a deletion of four amino acids at the Mc1r gene. Cellular immune response was morph specific. In males, but not in females, dark nestling mounted a lower PHA response than pale ones. Although correlative, our results raise the neglected possibility that MC1R has pleiotropic effects, suggesting a potential role of immune capacity and pathogen pressure on the maintenance of colour polymorphism in this species.     

Genetic and biochemical studies of the Agouti-attractin system

Pleiotropic effects of melanocortin signaling were first described nearly 100 years ago when mice carrying the lethal yellow (A(y)) allele of the Agouti coat color gene were recognized to develop increased growth and adiposity. Work from our laboratory and others over the last several years has demonstrated that the non-pigmentary effects of A(y) are caused by ectopic expression of Agouti protein, a paracrine signaling molecule whose normal function is to inhibit signaling through the melanocortin 1 receptor (Mc1r), but which can mimic the effects of Agouti-related protein (Agrp), a homologous neuropeptide produced in the medial portion of the arcuate nucleus that acts as a potent antagonist of the Mc3r and Mc4r. Recently we have used the genetics of pigmentation as an in vivo screening system to analyze other mutations in the Agouti-melanocortin pathway, leading to the identification of Attractin (Atrn), a widely expressed type I transmembrane protein that serves as an accessory receptor for Agouti protein. Surprisingly, homologs of Atrn are found in fruitflies and nematodes, even though Agouti and/or Agouti-related protein are found only in vertebrates. Insight into this apparent paradox now comes from studies of different Atrn alleles, in which we find hyperactivity, abnormal myelination, and widespread CNS vacuolation. We suggest that the neurodegenerative phenotype reflects the ancestral function of Atrn to facilitate and/or maintain cell-cell interactions in the nervous system. Expression in neurectodermal cells during vertebrate evolution may have allowed Atrn to be recruited by the Agouti-melanocortin system to control coat color.     

Adaptive evolution of fertilization proteins within a genus: variation in ZP2 and ZP3 in deer mice (Peromyscus)

Rapid evolution of reproductive proteins has been documented in a wide variety of taxa. In internally fertilized species, knowledge about the evolutionary dynamics of these proteins between closely related taxa is primarily limited to accessory gland proteins in the semen of Drosophila. Investigation of additional taxa and functional classes of proteins is necessary in order to determine if there is a general pattern of adaptive evolution of reproductive proteins between recently diverged species. We performed an evolutionary analysis of 2 egg coat proteins, ZP2 and ZP3, in 15 species of deer mice (genus Peromyscus). Both of these proteins are involved in egg-sperm binding, a critical step in maintaining species-specific fertilization. Here, we show that Zp2 and Zp3 gene trees are not consistent with trees based on nonreproductive genes, Mc1r and Lcat, where species formed monophyletic clades. In fact, for both of the reproductive genes, intraspecific amino acid variation was extensive and alleles were sometimes shared across species. We document positive selection acting on ZP2 and ZP3 and identify specific amino acid sites that are likely targets of selection using both maximum likelihood approaches and patterns of parallel amino acid change. In ZP3, positively selected sites are clustered in and around the region implicated in sperm binding in Mus, suggesting changes may impact egg-sperm binding and fertilization potential. Finally, we identify lineages with significantly elevated rates of amino acid substitution using a Bayesian mapping approach. These findings demonstrate that the pattern of adaptive reproductive protein evolution found at higher taxonomic levels can be documented between closely related mammalian species, where reproductive isolation has evolved recently.     

Characterization of the effect of Melanocortin 1 Receptor,a member of the hair color genetic locus,in alpaca (Lama pacos) fleece color differentiation

Little is known about the inheritance and influence of the fleece color gene Melanocortin 1 Receptor (MC1R). Melanocortin 1 Receptor (MC1R) is a well-known gene responsible for red versus black fleece pigmentation and is hypothesized to be a candidate gene for variation in alpaca coloration patterns. Inheritance of red versus black pigmentation in the context of genetic mutation is well understood in many domesticated mammals. We characterized the MC1R gene in a population of multi-colored alpacas in order to better understand its effect on coat color in the alpaca. Our characterization of the alpaca MC1R gene revealed 11 mutations. Of these one is a 4 bp deletion, four are silent mutations and six are single nucleotide polymorphisms (SNPs) that alter the amino acid sequence (T28V, M87V, S126G, T128I, S196F, R301C). No mutation correlated completely with fleece color in alpacas at the MC1R locus. This may be due to the epistatic relationship of MC1R with other coat color genes especially agouti signaling protein (ASIP).     

The genetics of adaptive coat color in gophers: coding variation at Mc1r is not responsible for dorsal color differences

The genetics of adaptation is a key problem in evolutionary biology. Pocket gophers of the species Thomomys bottae provide one of the most striking examples of coat color variation in mammals. Dorsal pelage color is strongly correlated with soil color across the range of the species, presumably reflecting the selective pressure exerted by predation. To investigate the genetic basis of coat color variation in T. bottae, we cloned and sequenced the melanocortin-1 receptor locus (Mc1r), a candidate pigmentation gene, in 5 dark and 5 light populations of the species. Our results show that, in contrast to many other species of mammals and other vertebrates, coding variation at Mc1r is not the main determinant of coat color variation in T. bottae. These results demonstrate that similar phenotypic variation may have a different genetic basis among different mammalian species.     

A Val85Met Mutation in Melanocortin-1 Receptor Is Associated with Reductions in Eumelanic Pigmentation and Cell Surface Expression in Domestic Rock Pigeons (Columba livia)

Variation in the melanocortin-1 receptor (Mc1r) is associated with pigmentation diversity in wild and domesticated populations of vertebrates, including several species of birds. Among domestic bird species, pigmentation variation in the rock pigeon ( Columba livia ) is particularly diverse. To determine the potential contribution of Mc1r variants to pigment diversity in pigeons, we sequenced Mc1r in a wide range of pigeon breeds and identified several single nucleotide polymorphisms, including a variant that codes for an amino acid substitution (Val85Met). In contrast to the association between Val85Met and eumelanism in other avian species, this change was associated with pheomelanism in pigeons. In vitro cAMP accumulation and protein expression assays revealed that Val85Met leads to decreased receptor function and reduced cell surface expression of the mutant protein. The reduced in vitro function is consistent with the observed association with reduced eumelanic pigmentation. Comparative genetic and cellular studies provide important insights about the range of mechanisms underlying diversity among vertebrates, including different phenotypic associations with similar mutations in different species.     

Agouti: from mouse to man, from skin to fat

The agouti protein regulates pigmentation in the mouse hair follicle producing a black hair with a subapical yellow band. Its effect on pigmentation is achieved by antagonizing the binding of alpha-melanocyte stimulating hormone (alpha-MSH) to melanocortin 1 receptor (Mc1r), switching melanin synthesis from eumelanin (black/brown) to phaeomelanin (red/yellow). Dominant mutations in the non-coding region of mouse agouti cause yellow coat colour and ectopic expression also results in obesity, type 11 diabetes, increased somatic growth and tumourigenesis. At least some of these pleiotropic effects can be explained by antagonism of other members of the melanocortin receptor family by agouti protein. The yellow coat colour is the result of agouti chronically antagonizing the binding of alpha-MSH to Mc1r and the obese phenotype results from agouti protein antagonizing the binding of alpha-MSH to Mc3r and/or Mc4r. Despite the existence of a highly homologous agouti protein in humans, agouti signal protein (ASIP), its role has yet to be defined. However it is known that human ASIP is expressed at highest levels in adipose tissue where it may antagonize one of the melanocortin receptors. The conserved nature of the agouti protein combined with the diverse phenotypic effects of agouti mutations in mouse and the different expression patterns of human and mouse agouti, suggest ASIP may play a role in human energy homeostasis and possibly human pigmentation.     

Discrete systems, cell-cell interactions and color pattern of animals. I. Conflicting dynamics and pattern formation

The morphogenesis of the color pattern of animals is modeled within the framework of a clonal model. The model takes into account the, sometimes conflicting, cell-cell and cell-background interactions. The color patterns of some reptiles and mammals are found to be consistent with the predictions of the model.