The Role of Melanocortin‐1 Receptor Polymorphism in Skin Cancer Risk Phenotypes

We have examined melanocortin‐1 receptor (MC1R) variant allele frequencies in the general population and in a collection of adolescent dizygotic and monozygotic twins to determine statistical associations of pigmentation phenotypes with increased skin cancer risk. This included hair and skin color, freckling, mole count and sun exposed skin reflectance. Nine variants were studied and designated as either strong R (OR = 63; 95% CI 32–140) or weak r (OR = 5; 95% CI 3–11) red hair alleles. Penetrance of each MC1R variant allele was consistent with an allelic model where effects were multiplicative for red hair but additive for skin reflectance. To assess the interaction of the brown eye color gene BEY2/OCA2 on the phenotypic effects of variant MC1R alleles we imputed OCA2 genotype in the twin collection. A modifying effect of OCA2 on MC1R variant alleles was seen on constitutive skin color, freckling and mole count. In order to study the individual effects of these variants on pigmentation phenotype we have established a series of human primary melanocyte strains genotyped for the MC1R receptor. These include strains which are MC1R wild‐type consensus, variant heterozygotes, and homozygotes for strong R alleles Arg151Cys and Arg160Trp. Ultrastructural analysis demonstrated that only consensus strains contained stage III and IV melanosomes in their terminal dendrites whereas Arg151Cys and Arg160Trp homozygous strains contained only immature stage I and II melanosomes. Such genetic association studies combined with the functional analysis of MC1R variant alleles in melanocytic cells should provide a link in understanding the association between pigmentary phototypes and skin cancer risk.     

Cutaneous pharmacologic cAMP induction induces melanization of the skin and improves recovery from ultraviolet injury in melanocortin 1 receptor‐intact or heterozygous skin

Homozygous loss of function of the melanocortin 1 receptor (MC1R) is associated with a pheomelanotic pigment phenotype and increased melanoma risk. MC1R heterozygosity is less well studied, although individuals inheriting one loss‐of‐function MC1R allele are also melanoma‐prone. Using the K14‐Scf C57BL/6J animal model whose skin is characterized by lifelong retention of interfollicular epidermal melanocytes like that of the human, we studied pigmentary, UV responses, and DNA repair capacity in the skin of variant Mc1r background. Topical application of forskolin, a skin‐permeable pharmacologic activator of cAMP induction to mimic native Mc1r signaling, increased epidermal eumelanin levels, increased the capacity of Mc1r‐heterozygous skin to resist UV‐mediated inflammation, and enhanced the skin's ability to clear UV photolesions from DNA. Interestingly, topical cAMP induction also promoted melanin accumulation, UV resistance, and accelerated clearance in Mc1r fully intact skin. Together, our findings suggest that heterozygous Mc1r loss is associated with an intermediately melanized and DNA repair‐proficient epidermal phenotype and that topical cAMP induction enhances UV resistance in Mc1r‐heterozygous or Mc1r‐wild‐type individuals by increasing eumelanin deposition and by improving nucleotide excision repair.     

Melanocortin‐1 receptor (MC1R) genotypes do not correlate with size in two cohorts of medium‐to‐giant congenital melanocytic nevi

Congenital melanocytic nevi (CMN) are cutaneous malformations whose prevalence is inversely correlated with projected adult size. CMN are caused by somatic mutations, but epidemiological studies suggest that germline genetic factors may influence CMN development. In CMN patients from the U.K., genetic variants in MC1R, such as p.V92M and loss‐of‐function variants, have been previously associated with larger CMN. We analyzed the association of MC1R variants with CMN characteristics in two distinct cohorts of medium‐to‐giant CMN patients from Spain (N = 113) and from France, Norway, Canada, and the United States (N = 53), similar at the clinical and phenotypical level except for the number of nevi per patient. We found that the p.V92M or loss‐of‐function MC1R variants either alone or in combination did not correlate with CMN size, in contrast to the U.K. CMN patients. An additional case–control analysis with 259 unaffected Spanish individuals showed a higher frequency of MC1R compound heterozygous or homozygous variant genotypes in Spanish CMN patients compared to the control population (15.9% vs. 9.3%; p = .075). Altogether, this study suggests that MC1R variants are not associated with CMN size in these non‐UK cohorts. Additional studies are required to define the potential role of MC1R as a risk factor in CMN development.     

Reverse genetic screen for loss‐of‐function mutations uncovers a frameshifting deletion in the melanophilin gene accountable for a distinctive coat color in Belgian Blue cattle

In the course of a reverse genetic screen in the Belgian Blue cattle breed, we uncovered a 10‐bp deletion (c.87_96del) in the first coding exon of the melanophilin gene (MLPH), which introduces a premature stop codon (p.Glu32Aspfs*1) in the same exon, truncating 94% of the protein. Recessive damaging mutations in the MLPH gene are well known to cause skin, hair, coat or plumage color dilution phenotypes in numerous species, including human, mice, dog, cat, mink, rabbit, chicken and quail. Large‐scale array genotyping undertaken to identify p.Glu32Aspfs*1 homozygous mutant animals revealed a mutation frequency of 5% in the breed and allowed for the identification of 10 homozygous mutants. As expression of a colored coat requires at least one wild‐type allele at the co‐dominant Roan locus encoded by the KIT ligand gene (KITLG), homozygous mutants for p.Ala227Asp corresponding with the missense mutation were excluded. The six remaining colored calves displayed a distinctive dilution phenotype as anticipated. This new coat color was named ‘cool gray’. It is the first damaging mutation in the MLPH gene described in cattle and extends the already long list of species with diluted color due to recessive mutations in MLPH and broadens the color palette of gray in this breed.     

Signatures of de‐domestication in autochthonous pig breeds and of domestication in wild boar populations from MC1R and NR6A1 allele distribution

Autochthonous pig breeds are usually reared in extensive or semi‐extensive production systems that might facilitate contact with wild boars and, thus, reciprocal genetic exchanges. In this study, we analysed variants in the melanocortin 1 receptor (MC1R) gene (which cause different coat colour phenotypes) and in the nuclear receptor subfamily 6 group A member 1 (NR6A1) gene (associated with increased vertebral number) in 712 pigs of 12 local pig breeds raised in Italy (Apulo‐Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano and Sarda) and south‐eastern European countries (Kr?kopolje from Slovenia, Black Slavonian and Turopolje from Croatia, Mangalitsa and Moravka from Serbia and East Balkan Swine from Bulgaria) and compared the data with the genetic variability at these loci investigated in 229 wild boars from populations spread in the same macro‐geographic areas. None of the autochthonous pig breeds or wild boar populations were fixed for one allele at both loci. Domestic and wild‐type alleles at these two genes were present in both domestic and wild populations. Findings of the distribution of MC1R alleles might be useful for tracing back the complex genetic history of autochthonous breeds. Altogether, these results indirectly demonstrate that bidirectional introgression of wild and domestic alleles is derived and affected by the human and naturally driven evolutionary forces that are shaping the Sus scrofa genome: autochthonous breeds are experiencing a sort of ‘de‐domestication’ process, and wild resources are challenged by a ‘domestication’ drift. Both need to be further investigated and managed.     

Loss‐of‐Function Mutations in MC4R Are Very Rare in the Greek Severely Obese Adult Population

Melanocortin‐4 receptor (MC4R) loss‐of‐function mutations are the commonest genetic cause of human monogenic obesity, so far. The contribution of MC4R coding mutations to severe obesity in the high‐obesity prone Greek population has not been investigated to date. We determined the MC4R coding sequence of 510 obese and 469 lean control subjects of Greek origin, and we estimated the prevalence and the penetrance on obesity of MC4R loss‐of‐function mutations. The functional impact of novel nonsynonymous variants detected was investigated in vitro. We found two novel synonymous mutations (L23L and I102I), four nonsynonymous mutations (T112M, S127L, N274S, and S295L), and two polymorphisms (V103I and I251L) previously described in literature. We also detected a novel mutation (L207V) in a severely obese 69‐year‐old female patient, although the mutation did not cosegregate with obesity in the corresponding pedigree and had no functional consequences on MC4R protein function. Loss‐of‐function mutations represented 75% of all nonsynonymous rare mutations identified among lean carriers and only 25% among obese subjects (P = 0.0001). The prevalence of loss‐of‐function mutations was lower in the obese group than in lean control subjects (0.20 vs. 0.64%) but this difference was not significant. Therefore, the estimated penetrance of deleterious MC4R mutations was very low (6.3%) in heterozygous Greek carriers of MC4R loss‐of‐function mutations. Our data suggest that MC4R loss‐of‐function mutations are rare in the Greek population. MC4R genetic deficiency is unlikely to explain the high propensity to develop severe obesity in this specific population.     

The miR‐379/miR‐410 cluster at the imprinted Dlk1‐Dio3 domain controls neonatal metabolic adaptation

In mammals, birth entails complex metabolic adjustments essential for neonatal survival. Using a mouse knockout model, we identify crucial biological roles for the miR‐379/miR‐410 cluster within the imprinted Dlk1‐Dio3 region during this metabolic transition. The miR‐379/miR‐410 locus, also named C14MC in humans, is the largest known placental mammal‐specific miRNA cluster, whose 39 miRNA genes are expressed only from the maternal allele. We found that heterozygote pups with a maternal—but not paternal—deletion of the miRNA cluster display partially penetrant neonatal lethality with defects in the maintenance of energy homeostasis. This maladaptive metabolic response is caused, at least in part, by profound changes in the activation of the neonatal hepatic gene expression program, pointing to as yet unidentified regulatory pathways that govern this crucial metabolic transition in the newborn's liver. Not only does our study highlight the physiological importance of miRNA genes that recently evolved in placental mammal lineages but it also unveils additional layers of RNA‐mediated gene regulation at the Dlk1‐Dio3 domain that impose parent‐of‐origin effects on metabolic control at birth and have likely contributed to mammal evolution.     

Effect of the MC1R gene on sexual dimorphism in melanin‐based colorations

Variants of the melanocortin‐1 receptor (MC1R) gene result in abrupt, naturally selected colour morphs. These genetic variants may differentially affect sexual dimorphism if one morph is naturally selected in the two sexes but another morph is naturally or sexually selected only in one of the two sexes (e.g. to confer camouflage in reproductive females or confer mating advantage in males). Therefore, the balance between natural and sexual selections can differ between MC1R variants, as suggest studies showing interspecific correlations between sexual dimorphism and the rate of nonsynonymous vs. synonymous amino acid substitutions at the MC1R. Surprisingly, how MC1R is related to within‐species sexual dimorphism, and thereby to sex‐specific selection, has not yet been investigated. We tackled this issue in the barn owl (Tyto alba), a species showing pronounced variation in the degree of reddish pheomelanin‐based coloration and in the number and size of black feather spots. We found that a valine (V)‐to‐isoleucine (I) substitution at position 126 explains up to 30% of the variation in the three melanin‐based colour traits and in feather melanin content. Interestingly, MC1R genotypes also differed in the degree of sexual colour dimorphism, with individuals homozygous for the II MC1R variant being 2 times redder and 2.5 times less sexually dimorphic than homozygous individuals for the VV MC1R variant. These findings support that MC1R interacts with the expression of sexual dimorphism and suggest that a gene with major phenotypic effects and weakly influenced by variation in body condition can participate in sex‐specific selection processes.     

Sex‐specific allelic transmission bias suggests sexual conflict at MC1R

Sexual conflict arises when selection in one sex causes the displacement of the other sex from its phenotypic optimum, leading to an inevitable tension within the genome – called intralocus sexual conflict. Although the autosomal melanocortin‐1‐receptor gene (MC1R) can generate colour variation in sexually dichromatic species, most previous studies have not considered the possibility that MC1R may be subject to sexual conflict. In the barn owl (Tyto alba), the allele MC1R_WHITE is associated with whitish plumage coloration, typical of males, and the allele MC1R_RUFOUS is associated with dark rufous coloration, typical of females, although each sex can express any phenotype. Because each colour variant is adapted to specific environmental conditions, the allele MC1R_WHITE may be more strongly selected in males and the allele MC1R_RUFOUS in females. We therefore investigated whether MC1R genotypes are in excess or deficit in male and female fledglings compared with the expected Hardy–Weinberg proportions. Our results show an overall deficit of 7.5% in the proportion of heterozygotes in males and of 12.9% in females. In males, interannual variation in assortative pairing with respect to MC1R explained the year‐specific deviations from Hardy–Weinberg proportions, whereas in females, the deficit was better explained by the interannual variation in the probability of inheriting the MC1R_WHITE or MC1R_RUFOUS allele. Additionally, we observed that sons inherit the MC1R_RUFOUS allele from their fathers on average slightly less often than expected under the first Mendelian law. Transmission ratio distortion may be adaptive in this sexually dichromatic species if males and females are, respectively, selected to display white and rufous plumages.     

Substantial differences in bias between single‐digest and double‐digest RAD‐seq libraries: A case study

The trade‐offs of using single‐digest vs. double‐digest restriction site‐associated DNA sequencing (RAD‐seq) protocols have been widely discussed. However, no direct empirical comparisons of the two methods have been conducted. Here, we sampled a single population of Gulf pipefish (Syngnathus scovelli) and genotyped 444 individuals using RAD‐seq. Sixty individuals were subjected to single‐digest RAD‐seq (sdRAD‐seq), and the remaining 384 individuals were genotyped using a double‐digest RAD‐seq (ddRAD‐seq) protocol. We analysed the resulting Illumina sequencing data and compared the two genotyping methods when reads were analysed either together or separately. Coverage statistics, observed heterozygosity, and allele frequencies differed significantly between the two protocols, as did the results of selection components analysis. We also performed an in silico digestion of the Gulf pipefish genome and modelled five major sources of bias: PCR duplicates, polymorphic restriction sites, shearing bias, asymmetric sampling (i.e., genotyping fewer individuals with sdRAD‐seq than with ddRAD‐seq) and higher major allele frequencies. This combination of approaches allowed us to determine that polymorphic restriction sites, an asymmetric sampling scheme, mean allele frequencies and to some extent PCR duplicates all contribute to different estimates of allele frequencies between samples genotyped using sdRAD‐seq versus ddRAD‐seq. Our finding that sdRAD‐seq and ddRAD‐seq can result in different allele frequencies has implications for comparisons across studies and techniques that endeavour to identify genomewide signatures of evolutionary processes in natural populations.     

Melanocortin‐1‐receptor (MC1R) variation is not associated with parasite burden in a neotropical bird,the bananaquit (Coereba flaveola)

It has been suggested that selection on melanocortin‐1‐receptor (MC1R) polymorphism, a common cause of melanic colour variation in vertebrates, results from pleiotropic effects of the gene in the immune system. Here we present the first test of whether MC1R variation is associated with differences in parasite abundance in a natural population. Bananaquits (Coereba flaveola) (Linnaeus, 1758) living on Grenada in the Caribbean exhibit a melanic plumage dimorphism as a result of a mutation in MC1R. The proportion of black individuals increases clinally towards the central, wetter parts of the island. We captured bananaquits through the cline and quantified parasite abundances. Avian malaria, feather mites, and mallophaga lice varied significantly in abundance across the cline; however, neither these infections, nor coccidia, nor arboviruses showed overall differences between the morphs. Feather mites tended to be more abundant on black individuals, in areas where the black morph was more common. This may result from differences in microhabitat use by the two morphs. These patterns do not support the idea that MC1R variation in itself results in differing susceptibility to parasites. © 2013 The Linnean Society of London     

Isolation and characterization of microsatellites in a marine food fish species,golden trevally Gnathanodon specious

The first set of six polymorphic microsatellites was isolated from a partial genomic DNA library of an important marine food fish Gnathanodon specious, and was characterized in 28 unrelated individuals. The allele number ranged from three to 12 with an average of seven/locus; average expected heterozygosity was 0.57 ranging from 0.10 to 0.87, whereas the observed heterozygosity ranged from 0.11 to 0.96 (average: 0.57). All six markers conform to Hardy–Weinberg and linkage equilibria. These markers will be useful for accessing genetic population structure of wild populations and for selective breeding program in fish farms.     

Cost‐effective development of highly polymorphic microsatellite in Japanese quail facilitated by next‐generation sequencing

Next‐generation sequencing technologies permit rapid and cost‐effective identification of numerous putative microsatellite loci. Here, from the genome sequences of Japanese quail, we developed microsatellite markers containing dinucleotide repeats and employed these for characterisation of genetic diversity and population structure. A total of 385 individuals from 12 experimental and one wild‐derived Japanese quail lines were genotyped with newly developed autosomal markers. The maximum number of alleles, expected heterozygosity and polymorphic information content (PIC) per locus were 10, 0.80 and 0.77 respectively. Approximately half of the markers were highly informative (PIC ≥ 0.50). The mean number of alleles per locus and observed heterozygosity within a line were in the range of 1.3–4.1 and 0.11–0.53 respectively. Compared with the wild‐derived line, genetic diversity levels were low in the experimental lines. Genetic differentiation (F_ST) between all pairs of the lines ranged from 0.13 to 0.83. Genetic clustering analyses based on multilocus genotypes of individuals showed that most individuals formed clearly defined clusters corresponding to the origins of the lines. These results suggest that Japanese quail experimental lines are highly structured. Microsatellite markers developed in this study may be effective for future genetic studies of Japanese quail.     

Abundant variations of <Emphasis Type="Italic">MC4R</Emphasis> gene revealed by Phylogenies of Yak (<Emphasis Type="Italic">Bos grunniens</Emphasis>) and other mammals

MC4R gene was proved to play important roles in body weight regulation in many mammals and exhibit higher homology among different species. The mutations MC4R significantly correlated to the restricted feeding weight, fat deposition and energy balance. In this work, ORF sequences of MC4R gene of Bos grunniens were cloned and phylogenetic relationships of yak and other mammals were analyzed on the basis of MC4R genes. Totally 290 variable sites were examined in 25 sequences from 22 different mammals, and 23 haplotypes were defined with a haplotype diversity of 0.9900. All the sequences were clustered into phylogenetic clades representing different orders or families. The individuals of Bos grunniens, Bos taurus and Ovis aries which belonged to the family of Bovidae were more divergent from the other orders or families and bovid animals may have branched out from the phylogenetic tree earlier than other mammals analyzed during 450 million years of vertebrate evolution. Amino acid sequences inferred from MC4R genes exhibited 54 variable sites, while high conservation of MC4R was observed within the same order or family. We concluded that coding region of MC4R gene displayed abundant variations among different mammal phylogenetic clades, whereas, the conservation of MC4R within order or family could be explained that MC4R gene may have been subjected to substantial constraints or strong purifying selection during several million years of mammal evolution.     

Non‐synonymous SNPs in MC1R gene are associated with the extended black variant in domestic ducks (Anas platyrhynchos)

In this study, we performed a sequence characterization of the duck melanocortin 1 receptor (alpha‐melanocyte stimulating hormone receptor) (MC1R) gene to analyze the relationship between MC1R polymorphism and the extended black variant in domestic ducks based on the extended black (E) and non‐extended black (e⁺) allele hypothesis of the duck MC1R gene. Both c.52G>A and c.376G>A substitutions are highly associated with the duck extended black variant (P < 0.01), but the novel c.52G>A substitution is more of a fit for the allele hypothesis of the duck MC1R gene.     

MC1R variation and melanoma risk in the Swedish population in relation to clinical and pathological parameters

The genetic background of cutaneous malignant melanoma (CMM) includes both germ line aberrations in high‐penetrance genes, like CDKN2A, and allelic variation in low‐penetrance genes like the melanocortin‐1 receptor gene, MC1R. Red‐hair colour associated MC1R alleles (RHC) have been associated with red hair, fair skin and risk of CMM. We investigated MC1R and CDKN2A variation in relation to phenotype, clinical factors and CMM risk in the Swedish population. The study cohort consisted of sporadic primary melanoma patients, familial melanoma patients and a control group. An allele‐dose dependent increase in melanoma risk for carriers of variant MC1R alleles (after adjusting for phenotype), with an elevated risk among familial CMM patients, was observed. This elevated risk was found to be significantly associated with an increased frequency of dysplastic nevi (DN) among familial patients compared to sporadic patients. MC1R variation was found to be less frequent among acral lentiginous melanomas (ALM) and dependent on tumour localisation. No association was found between CDKN2A gene variants and general melanoma risk. Two new variants in the POMC gene were identified in red haired individuals without RHC alleles.     

Beyond mean allelic effects: A locus at the major color gene MC1R associates also with differing levels of phenotypic and genetic (co)variance for coloration in barn owls

The mean phenotypic effects of a discovered variant help to predict major aspects of the evolution and inheritance of a phenotype. However, differences in the phenotypic variance associated to distinct genotypes are often overlooked despite being suggestive of processes that largely influence phenotypic evolution, such as interactions between the genotypes with the environment or the genetic background. We present empirical evidence for a mutation at the melanocortin‐1‐receptor gene, a major vertebrate coloration gene, affecting phenotypic variance in the barn owl, Tyto alba. The white MC1R allele, which associates with whiter plumage coloration, also associates with a pronounced phenotypic and additive genetic variance for distinct color traits. Contrarily, the rufous allele, associated with a rufous coloration, relates to a lower phenotypic and additive genetic variance, suggesting that this allele may be epistatic over other color loci. Variance differences between genotypes entailed differences in the strength of phenotypic and genetic associations between color traits, suggesting that differences in variance also alter the level of integration between traits. This study highlights that addressing variance differences of genotypes in wild populations provides interesting new insights into the evolutionary mechanisms and the genetic architecture underlying the phenotype.     

Association of the MC4R V103I Polymorphism With Obesity: A Chinese Case–control Study and Meta‐analysis in 55,195 Individuals

Recently, large‐scaled studies suggested a negative association of the infrequent allele of the melanocortin‐4 receptor (MC4R) V103I polymorphism with obesity. We conducted a Chinese case–control study, meta‐analysis in East Asians and in all populations, in order to further assess the association between the V103I polymorphism and reduced body weight and to explore whether the association varies among different ethnic groups. We conducted a case–control study to analyze this polymorphism in 2,012 children of two independent study groups from Beijing, China, no association was found between the V103I polymorphism and obesity or obesity‐related phenotypes (P > 0.10). In the meta‐analysis of 3,526 individuals from six East Asian studies, I103 carriers had a 31% lower risk for obesity (odds ratio (OR) = 0.69, 95% confidence interval (CI): 0.50–0.94, P = 0.02). Subsequently, we performed a large meta‐analysis of the six East Asian studies and 31 studies of other ethnic groups, involving 55,195 individuals with 19,822 obese cases and 35,373 nonobese controls. In total, the individuals with I103 allele had a 21% lower risk for obesity (OR = 0.79, 95% CI: 0.71–0.88, P < 0.0001). In conclusion, this study confirmed and extended the previous findings, suggesting the MC4R V103I polymorphism is negatively associated with human obesity. It provides significant evidence for the association in East Asian populations. Further large‐scaled studies in East Asian populations are needed to validate the association.