Skip to main content
  • Short communication
  • Open access
  • Published:

A missense mutation in the agouti signaling protein gene (ASIP) is associated with the no light points coat phenotype in donkeys

Abstract

Background

Seven donkey breeds are recognized by the French studbook and are characterized by a black, bay or grey coat colour including light cream-to-white points (LP). Occasionally, Normand bay donkeys give birth to dark foals that lack LP and display the no light points (NLP) pattern. This pattern is more frequent and officially recognized in American miniature donkeys. The LP (or pangare) phenotype resembles that of the light bellied agouti pattern in mouse, while the NLP pattern resembles that of the mammalian recessive black phenotype; both phenotypes are associated with the agouti signaling protein gene (ASIP).

Findings

We used a panel of 127 donkeys to identify a recessive missense c.349 T > C variant in ASIP that was shown to be in complete association with the NLP phenotype. This variant results in a cysteine to arginine substitution at position 117 in the ASIP protein. This cysteine is highly-conserved among vertebrate ASIP proteins and was previously shown by mutagenesis experiments to lie within a functional site. Altogether, our results strongly support that the identified mutation is causative of the NLP phenotype.

Conclusions

Thus, we propose to name the c.[349 T > C] allele in donkeys, the a nlp allele, which enlarges the panel of coat colour alleles in donkeys and ASIP recessive loss-of-function alleles in animals.

Background

Mutations in the gene ASIP (agouti signaling protein) result in various coat patterns in domestic mammals (http://omia.angis.org.au) including mouse (www.informatics.jax.org), dog [1], cat [2], rabbit [3], horse [4], sheep [5-8], rat [9] and alpaca [10]. Only a few coat colours, patterns and textures have been described in domestic donkeys (Equus asinus). In donkeys, the coat colour can be white or coloured, i.e. black, bay, grey and red with or without white spotting; hair texture is variable and includes the longhair phenotype, in addition to the common shorthair phenotype. Recently, we started to investigate the molecular aetiology of these phenotypes and identified three underlying loss-of-function mutations in the MC1R (melanocortin 1 receptor) and FGF5 (fibroblast growth factor) genes that are responsible respectively, for the red colour and longhair phenotype in donkeys [11,12]. Most coloured donkeys are born with a pangare or light points (LP) pattern that associates cream to grey-white hair on the belly, around the muzzle and around the eyes. In the American miniature donkey breed, all coat colours and patterns are admissible and foals with a no light points (NLP) coat are often obtained from LP breeding stock (Figure 1). This has led breeders to suspect a recessive inheritance pattern for the NLP pattern. For the seven French donkey breeds (Pyrenean, Berry Black, Poitou, Cotentin, Provence, Bourbonnais and Normand), the NLP pattern is not officially recognized. However, dark NLP donkeys are occasionally born to bay Normand parents (Figure 1).

Figure 1
figure 1

No light points phenotype in donkeys. Most coloured donkeys show a light cream to white coat on their belly and around their eyes and muzzle (Normand donkey, left). Bay Normand donkeys occasionally give birth to no light points (NLP) donkeys that are not officially recognized by the studbook (NLP donkeys, middle). The NLP phenotype is recognized in American miniature donkeys (NLP American miniature donkey, right).

Methods

Animals and ethics statement

One hundred and twenty seven donkeys from six breeds were included in the study. They were all sampled in France between September 2012 and October 2014 and included Normand (n = 35), Provence (n = 14), Poitou (n = 13), Pyrenean (n = 13) and Berry Black (n = 2) breeds and miniature donkeys (n = 50).

All donkeys were included at the owners’ request. Pictures and hair samples were sent directly by owners or collected by a veterinarian (MA). All animals were client-owned donkeys on which no harmful invasive procedure was performed; thus, according to the legal definitions in Europe (Subject 5f of Article 1, Chapter I of the Directive 2010/63/UE of the European Parliament and of the Council), no animal experiment was carried out.

DNA extraction

DNA was extracted from hair roots using a Maxwell® 16 Instrument (Promega Corporation, Madison, USA), according to the manufacturer’s protocol.

Sequencing of ASIP and genotyping

Reference genomic sequences were extracted from Ensembl (Equine ASIP gene, ENSECAG00000004241). PCR and sequencing primers were designed using Primer3 [13]. The three exons were amplified using three sets of primers [See Additional file 1: Table S1]. PCR amplicons were sequenced using Sanger dideoxy sequencing in both forward and reverse directions by GATC Biotech (GATC Biotech AG, Konstanz, Germany). Electropherograms were manually inspected with Chromas Lite (Technelysium Pty Ltd, South Brisbane, Australia). Multiple alignments were performed using Multalin ([14]; http://multalin.toulouse.inra.fr).

Protein sequence comparisons and impact of sequence variations

ASIP amino acid sequences from various species were collected from Ensembl (mouse: ENSMUST00000109697; horse: ENSECAT00000004772; cow: ENSBTAT00000048322; sheep: ENSOART00000010128; dog: ENSCAFT00000038625; cat: ENSFCAT00000011040; human: ENST00000568305; chicken: ENSGALT00000044768; zebrafish: ENSDART00000113083). Multiple alignments were performed using Multalin ([14]; http://multalin.toulouse.inra.fr). The putative impact of missense mutations was assessed using three different software, namely PolyPhen-2 (HumVar-trained PolyPhen-2 designed to distinguish mutations with drastic effects from other variations, including abundant mildly deleterious alleles, ([15]; http://genetics.bwh.harvard.edu/pph2/), SNAP ([16]; www.rostlab.org/services/snap/submit) and PROVEAN ([17]; http://provean.jcvi.org/seq_submit.php). The different domains of ASIP were schematized according to the previously published ASIP structure in mouse [18].

Accession numbers

Genomic coding sequences of ASIP from bay LP and NLP Normand donkeys were submitted to GeneBank. Accession numbers are KJ126712 for the LP allele and KP717040 for the NLP mutant allele.

Findings

Because whole-genome mapping tools are still lacking for donkey, we decided to screen directly for variants that affect ASIP function in two NLP and two LP control Normand donkeys that originated from a comprehensive panel of 127 donkeys from six breeds. The Ensembl ASIP equine genomic sequence was used to design three sets of intronic primers [See Additional file 1: Table S1] that allowed successful amplification of the three exons of the donkey ASIP gene. Then, we sequenced the three ASIP exonic amplicons and performed pair-wise base-to-base comparisons of the sequences between the LP and NLP donkeys and a horse reference sequence; we found that the coding sequences and the sequences that covered intron-exon boundaries were highly conserved between horse and donkey and detected only two variants between donkey sequences and the bay horse reference sequence [See Additional file 1: Table S2]. Only the c.349 T > C SNP (single nucleotide polymorphism) produced a substitution p.(Cys117Arg) and was consistent with the recessive mode of inheritance of the NLP pattern. Indeed both NLP donkeys were homozygous C/C for the mutant allele of c.349 T > C SNP, while one control donkey was heterozygous C/T, and the other control donkey and the bay horse were homozygous T/T for the reference allele [See Additional file 1: Table S2]. The second non-coding variant was located in the 3’UTR (untranslated region) region of the ASIP gene and was not associated with the NPL phenotype [See Additional file 1: Table S2].

PROVEAN, PolyPhen-2 and SNAP predicted that the p.(Cys117Arg) substitution was deleterious. Hence, the complete cohort of 127 donkeys was genotyped for SNP c.349 T > C. The nine NLP donkeys, including three NLP Normand and six NLP miniature donkeys, were all homozygous C/C, while the 118 LP donkeys were either homozygous T/T (n = 104) or heterozygous C/T (n = 14). The three NLP Normand donkeys were born to a single male mated with three females, which were all four heterozygous C/T. The complete concordance between the recessively-inherited NLP pattern and the c.349 T > C variant (Table 1) supported our hypothesis that this SNP is associated with the NLP trait (Chi square test p = 1.86 × 10−29).

Table 1 Genotypes for the c.349T > C variant in donkeys

To estimate the functional importance of the donkey ASIP cysteine 117 amino acid, we aligned the donkey ASIP protein sequence with the ASIP sequences of nine vertebrates and found that is was fully conserved (Figure 2). This result confirmed the 100% conservation previously reported for the 10 cysteine amino acids of the C-terminal Cys-rich domain of ASIP [18-21] the functional role of which was investigated by mutagenesis experiments. Perry and collaborators reported that in mouse, 13 mutated ASIP proteins displayed a partial (n = 4) or a total (n = 9) loss of activity [21]. In particular, they found that eight of the 10 cysteines located in the Cys-rich C-terminal tail of ASIP, including the murine cysteine 113 that corresponds to the donkey cysteine 117, were critical for protein activity [21]. Altogether these results strongly support that, in donkeys, the ASIP cysteine 117 has an essential role for ASIP function.

Figure 2
figure 2

Alignment between ASIP protein sequences from 10 vertebrate species. ASIP amino acid sequences from various species were collected from Ensembl (mouse: ENSMUST00000109697; horse: ENSECAT00000004772; cow: ENSBTAT00000048322; sheep: ENSOART00000010128; dog: ENSCAFT00000038625; cat: ENSFCAT00000011040; human: ENST00000568305; chicken: ENSGALT00000044768; zebrafish: ENSDART00000113083). ASIP sequences are identified with the name of the species on the left. The mouse ASIP sequence (reference sequence) is at the top of the alignment. Non-conserved residues in the 10 species analyzed are shown in grey. Conserved residues are indicated in black within the reference sequence and represented by black dots in other sequences. Dashes represent deletions. The conserved C117 residue in donkey ASIP that corresponds to the murine C113 residue is circled. The 10 functional cysteine residues that have been shown to be involved in the activity of ASIP are indicated with stars. Non-agouti mutations identified in other domestic mammalian species are reported above the corresponding position. Dog and alpaca mutations consist in missense mutations and a 19 amino acids (aa) in-frame deletion. Rabbit, rat, cat and horse mutations are frameshifts. Mouse and sheep regulatory mutations are not shown.

In conclusion, the complete association between the c.349 T > C mutation and the NLP phenotype and its inheritance pattern, on the one hand, and the high probability that the resulting substitution of the conserved cysteine 117 residue leads to loss of function of the mutated protein, on the other hand, support that this mutation is responsible for the NLP phenotype in donkeys. We thus propose to name the c.[349 T > C] allele, which can be easily detected with a DNA test, the a nlp allele in donkeys.

References

  1. Kerns JA, Newton J, Berryere TG, Rubin EM, Cheng JF, Schmutz SM, et al. Characterization of the dog Agouti gene and a nonagouti mutation in German Shepherd dogs. Mamm Genome. 2004;15:798–808.

    Article  CAS  PubMed  Google Scholar 

  2. Eizirik E, Yuhki N, Johnson WE, Menotti-Raymond M, Hannah SS, O’Brien SJ. Molecular genetics and evolution of melanism in the cat family. Curr Biol. 2003;13:448–53.

    Article  CAS  PubMed  Google Scholar 

  3. Fontanesi L, Forestier L, Allain D, Scotti E, Beretti F, Deretz-Picoulet S, et al. Characterization of the rabbit agouti signaling protein (ASIP) gene: transcripts and phylogenetic analyses and identification of the causative mutation of the nonagouti black coat colour. Genomics. 2010;95:166–75.

    Article  CAS  PubMed  Google Scholar 

  4. Rieder S, Taourit S, Mariat D, Langlois B, Guerin G. Mutations in the agouti (ASIP), the extension (MC1R), and the brown (TYRP1) loci and their association to coat color phenotypes in horses (Equus caballus). Mamm Genome. 2001;12:450–5.

    Article  CAS  PubMed  Google Scholar 

  5. Norris BJ, Whan VA. A gene duplication affecting expression of the ovine ASIP gene is responsible for white and black sheep. Genome Res. 2008;18:1282–93.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Fontanesi L, Dall’Olio S, Beretti F, Portolano B, Russo V. Coat colours in the Massese sheep breed are associated with mutations in the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes. Animal. 2011;5:8–17.

    Article  CAS  PubMed  Google Scholar 

  7. Royo LJ, Alvarez I, Arranz JJ, Fernandez I, Rodriguez A, Perez-Pardal L, et al. Differences in the expression of the ASIP gene are involved in the recessive black coat colour pattern in sheep: evidence from the rare Xalda sheep breed. Anim Genet. 2008;39:290–3.

    Article  CAS  PubMed  Google Scholar 

  8. Han JL, Yang M, Yue YJ, Guo TT, Liu JB, Niu CE, et al. Analysis of agouti signaling protein (ASIP) gene polymorphisms and association with coat color in Tibetan sheep (Ovis aries). Genet Mol Res. 2015;14:1200–9.

    Article  CAS  PubMed  Google Scholar 

  9. Kuramoto T, Nomoto T, Sugimura T, Ushijima T. Cloning of the rat agouti gene and identification of the rat nonagouti mutation. Mamm Genome. 2001;12:469–71.

    Article  CAS  PubMed  Google Scholar 

  10. Chandramohan B, Renieri C, La Manna V, La Terza A. The alpaca agouti gene: genomic locus, transcripts and causative mutations of eumelanic and pheomelanic coat color. Gene. 2013;521:303–10.

    Article  CAS  PubMed  Google Scholar 

  11. Legrand R, Tiret L, Abitbol M. Two recessive mutations in FGF5 are associated with the long-hair phenotype in donkeys. Genet Sel Evol. 2014;46:65.

    Article  PubMed Central  Google Scholar 

  12. Abitbol M, Legrand R, Tiret L. A missense mutation in melanocortin 1 receptor is associated with the red coat colour in donkeys. Anim Genet. 2014;45:878–80.

    Article  CAS  PubMed  Google Scholar 

  13. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3–new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Corpet F. Multiple sequence alignment with hierarchical clustering. Nucleic Acids Res. 1988;16:10881–90.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Bromberg Y, Rost B. SNAP: predict effect of non-synonymous polymorphisms on function. Nucleic Acids Res. 2007;35:3823–35.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Choi Y, Sims GE, Murphy S, Miller JR, Chan AP. Predicting the functional effect of amino acid substitutions and indels. PLoS One. 2012;7:e46688.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Miltenberger RJ, Wakamatsu K, Ito S, Woychik RP, Russell LB, Michaud EJ. Molecular and phenotypic analysis of 25 recessive, homozygous-viable alleles at the mouse agouti locus. Genetics. 2002;160:659–74.

    PubMed Central  CAS  PubMed  Google Scholar 

  19. Dinulescu DM, Cone RD. Agouti and agouti-related protein: analogies and contrasts. J Biol Chem. 2000;275:6695–8.

    Article  CAS  PubMed  Google Scholar 

  20. McNulty JC, Jackson PJ, Thompson DA, Chai B, Gantz I, Barsh GS, et al. Structures of the agouti signaling protein. J Mol Biol. 2005;346:1059–70.

    Article  CAS  PubMed  Google Scholar 

  21. Perry WL, Nakamura T, Swing DA, Secrest L, Eagleson B, Hustad CM, et al. Coupled site-directed mutagenesis/transgenesis identifies important functional domains of the mouse agouti protein. Genetics. 1996;144:255–64.

    PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We wish to thank owners and breeders for providing samples. Pictures for Normand donkey and miniature donkey were provided through the courtesy of INAM (Institut National Asin et Mulassier, France) and Carole De Koster (Francarolis Farm), respectively. This project was funded by the French National Donkey Institute, INAM.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Marie Abitbol.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MA designed the study and collected samples. RL and MA performed the experiments. MA and LT wrote the manuscript. All authors read and approved the final manuscript.

Additional file

Additional file 1: Table S1.

PCR and sequencing primers. Sequences and PCR temperatures from the intronic primers that were used to amplify and sequence the three ASIP coding exons. Table S2. Genomic variants in ASIP identified between donkey sequences and the horse reference sequence. Sequence variants in ASIP identified between two NLP donkeys, two control donkeys and the horse reference sequence (coding sequences and 5’end of the 3’UTR).

Rights and permissions

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Abitbol, M., Legrand, R. & Tiret, L. A missense mutation in the agouti signaling protein gene (ASIP) is associated with the no light points coat phenotype in donkeys. Genet Sel Evol 47, 28 (2015). https://doi.org/10.1186/s12711-015-0112-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12711-015-0112-x

Keywords