Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Highly Accessed Research article

Testing whether macroevolution follows microevolution: Are colour differences among swans (Cygnus) attributable to variation at the MC1R locus?

Marie A Pointer* and Nicholas I Mundy

Author Affiliations

Department of Zoology, University of Cambridge, Downing St, Cambridge, CB2 3EJ, UK

For all author emails, please log on.

BMC Evolutionary Biology 2008, 8:249  doi:10.1186/1471-2148-8-249

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/8/249


Received:6 May 2008
Accepted:12 September 2008
Published:12 September 2008

© 2008 Pointer and Mundy; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The MC1R (melanocortin-1 receptor) locus underlies intraspecific variation in melanin-based dark plumage coloration in several unrelated birds with plumage polymorphisms. There is far less evidence for functional variants of MC1R being involved in interspecific variation, in which spurious genotype-phenotype associations arising through population history are a far greater problem than in intraspecific studies. We investigated the relationship between MC1R variation and plumage coloration in swans (Cygnus), which show extreme variation in melanic plumage phenotypes among species (white to black).

Results

The two species with melanic plumage, C. atratus and C. melanocoryphus (black and black-necked swans respectively), both have amino acid changes at important functional sites in MC1R that are consistent with increased MC1R activity and melanism. Reconstruction of MC1R evolution over a newly generated independent molecular phylogeny of Cygnus and related genera shows that these putative melanizing mutations were independently derived in the two melanic lineages. However, interpretation is complicated by the fact that one of the outgroup genera, Coscoroba, also has a putative melanizing mutation at MC1R that has arisen independently but has nearly pure white plumage. Epistasis at other loci seems the most likely explanation for this discrepancy. Unexpectedly, the phylogeny shows that the genus Cygnus may not be monophyletic, with C. melanocoryphus placed as a sister group to true geese (Anser), but further data will be needed to confirm this.

Conclusion

Our study highlights the difficulty of extrapolating from intraspecific studies to understand the genetic basis of interspecific adaptive phenotypic evolution, even with a gene whose structure-function relationships are as well understood as MC1R as confounding variation make clear genotype/phenotype associations difficult at the macroevolutionary scale. However, the identification of substitutions in the black and black-necked swan that are known to be associated with melanic phenotypes, suggests Cygnus may be another example where there appears to be convergent evolution at MC1R. This study therefore provides a novel example where previously described intraspecific genotype/phenotype associations occur at the macroevolutionary level.

Background

A key question in evolutionary biology is whether the genetic basis of adaptations in macroevolution is the same as those of microevolution. This issue has been difficult to address since the different timescales involved lead to different methods of analysis, with population genetics being important for understanding microevolution and comparative phylogenetics being used for macroevolution. The growing understanding of the genetic basis of adaptive phenotypic traits, however, provides a new impetus to studies of this question.

Pigmentation is an excellent model for studying the genetic basis of phenotypic variation. In particular, melanization in vertebrates shows striking variability in pattern and colour that is frequently under strong natural or sexual selection, with the genetic control often well understood [1]. Many genes are involved in melanin production and deposition, however the melanocortin-1 receptor protein encoded at the MC1R locus plays a crucial role in controlling the type of melanin synthesized by melanocytes [2]. High activity of MC1R (either due to the presence of MSH agonist or high intrinsic MC1R activity) leads to synthesis of black/brown eumelanin, whereas low activity (due to an absence of MSH or presence of ASIP inhibitor) leads either to synthesis of reddish phaeomelanin, or an absence of melanin synthesis. Numerous studies have documented a strong association between colour variation within natural vertebrate populations and variations within the coding sequence of the MC1R locus [3-9]. One example is the lesser snow goose (Anser c. caerulescens), where the V85M substitution was found in one or both MC1R alleles in all melanic geese but was absent in white phase geese [6]. Interestingly, the same mutation was recently found to be associated with melanism in another bird, the red-footed booby (Sula sula) [10].

To date, almost all of the associations described between vertebrate coloration and MC1R have been intraspecific, and hence relating to microevolution. A previous study found no association between MC1R variation among species of Old World leaf warblers and the presence of small unmelanized patches [11], but it would have been surprising if such subtle plumage variation was affected by MC1R as polymorphisms known to associate with MC1R variation are more widespread across the pelage or plumage and not within small patches. A key question is whether MC1R variation plays a similar role in colour variation over longer evolutionary timescales as it is not yet clear whether similar or different molecular mechanisms underpin phenotypic variation at the macro- and microevolutionary scale. Methodologically, an absence of interbreeding means that the most direct means of demonstrating linkage/association between genotype and phenotype, i.e. in families or populations segregating for colour variation, is generally unavailable for interspecific studies. However, there is a good and growing understanding of MC1R structure-function relationships, which has arisen from in vitro studies of the biochemical function of MC1R [12] and studies of genotype-phenotype associations in domestic and wild species (for reviews see [13,1]. A crucial finding from these studies is that the same mutation leads to similar changes in MC1R function independent of other sequence variation in the protein. In addition, analyses of patterns of selection provide another route to identifying functionally important variation across phylogenies as demonstrated for MC1R and sexual dimorphism in galliformes [14].

The swans (Cygnus) provide a simple system of colour variation for investigating whether MC1R has a role in interspecific colour differences, and are in the same family as the lesser snow goose which shows intraspecific association between MC1R genotype and plumage melanism. All swan species are sexually monomorphic. The four species in the Northern hemisphere, Cygnus olor (mute swan), C. cygnus (whooper swan), C. buccinator (trumpeter swan) and C. columbianus (tundra swan) display pure white plumage. In the Southern hemisphere, C. atratus (black swan) has completely black plumage and C. melanocoryphus (black-necked swan) has white plumage on the torso and wings with black feathers on the head and neck (Figure 1).

thumbnailFigure 1. Reconstruction of MC1R evolution over an independent molecular phylogeny of swans and geese. A maximum likelihood tree of swans and geese using CytB and LDH B data. Amino acid substitutions at important functional sites in MC1R are shown in red, all other substitutions are shown in black. Numbers above branches are bootstrap support values from 100 replicates. Numbers below branches are the posterior support probabilities from Bayesian analysis. Anas platyrhynchos sequences were obtained from GenBank; Accession numbers EU009397 (CytB) and X68810 (LDH B).

In this study we address whether MC1R is associated with melanin-based plumage coloration in Cygnus. We determine whether variation at important functional sites in MC1R is associated with colour variation, and whether these variants are evolutionarily primitive or derived. In order to do this, we reconstruct MC1R evolution over a newly generated, independent molecular phylogeny of the genus Cygnus, since phylogenetic relationships both within the genus and among Cygnus and other genera of geese and swans (Anser, Cereopsis, Coscoroba) are poorly understood [15]

Results

Phylogenetic reconstruction of swans and geese using CytB and LDH B data

In total an 812 bp alignment (502 bp Cytochrome B, 310 bp Lactate dehydrogenase B) of all taxa was obtained. In MODELTEST, the TrN + I + G model [16] provided the best fit to the data. The maximum likelihood reconstruction obtained with this model in PAUP is shown in Figure 1. The identical topology was recovered in Bayesian analyses, and bootstrap support percentages and Bayesian support strengths are commensurate. Surprisingly, the genus Cygnus is not recovered as a monophyletic clade: four of the Cygnus species form a well-supported monophyletic group but C. melanocoryphus is placed as a sister group to Anser with moderate support. Elsewhere on the phylogeny, the coscoroba swan and the Cape Barren goose are strongly supported as sister taxa, which has previously been suggested based on the presence of a CR1 retrotransposon in LDH B [17].

MC1R variation and reconstruction in Cygnus

The full 945 bp coding sequence of MC1R was obtained for all species studied. There was substantial variation in the open reading frame among different species of Cygnus, with 19 synonymous and 33 non-synonymous variable positions found. All five species had unique predicted protein sequences. Sites of functional importance were identified from previous studies based on linkage/association between MC1R variants and coat colour, and in vitro studies of MC1R function, including site-directed mutagenesis. Three putatively important mutations (at sites 92, 100 and 230) were identified among Cygnus species: one of these (230H) occurs in C. atratus, two (92A, 100K) occur in C. melanocoryphus, and none occur in any of the white species. Reconstruction of MC1R evolution over the molecular phylogeny in PAML shows that these mutations were derived in terminal lineages (Figure 1).

All five individuals of C. atratus sampled have an R230H substitution. The same substitution is associated with melanism in arctic skuas (Stercorarius parasiticus) [6] and a histidine at this position is also associated with melanism in rock pocket mice (Chaetodipus intermedius) [18].

Cygnus melanocoryphus has two putative functional mutations that are both associated with a change of charge, an E92A mutation and an E100K mutation. A charge-changing mutation at the position homologous to 92 (E92K) is associated with melanism in mice, chicken, quail and bananaquit (Coereba flaveola)[19-21,8] and in vitro studies in mouse and chicken have shown that the mutation causes constitutive activation of MC1R [19,22]. The E100K mutation occurs at a strongly conserved site that when deleted is associated with melanism in jaguar and jaguarundi [4].

MC1R variation in other species, and tests of selection

Among other taxa studied, only Coscoroba had a mutation at an important functional site in MC1R. Coscoroba shares a derived R230H substitution with C. atratus, and maximum likelihood reconstruction over the phylogeny using PAML predicts an R residue at position 230 in the inferred last common ancestor between Coscoroba and Cygnus (posterior probability = 0.997) strongly suggesting that the R230H substitution occurred independently in the two lineages.

Lineage-specific analyses in PAML showed that the melanic C. atratus and C. melanocoryphus had a higher estimated dN/dS ratio (0.21) than other Cygnus lineages (0.12), but this was not significant. Site-specific analyses over the whole phylogeny did not provide any evidence for a category of sites under positive selection (not shown).

Discussion

MC1R variation and evolution of plumage coloration in swans and geese

Within the genus Cygnus, it is striking that the only variation at known sites of functional importance in MC1R was found in the two melanic species, C. atratus and C. melanocoryphus. These species both had independent derived mutations (R230H in C. atratus, E92A and E100K in C. melanocoryphus) consistent with an increase in melanism. Taken at face value, this suggests that MC1R is of functional importance for melanism among species of Cygnus, and also suggests that the ancestral swan may have been white.

However, the interpretation is clearly complicated by the finding of the R230H mutation in Coscoroba, which has close to pure white plumage. Either the R230H mutation has no influence on plumage coloration or it has a differential effect on phenotype in C. atratus and Coscoroba. The association of a 230H with melanism in two distantly related species (Arctic skua and rock pocket mouse) strongly suggests that this mutation does influence MC1R function. In addition, a growing body of evidence from in vitro and association studies shows that mutations in the same positions of MC1R generally have similar effects on MC1R function in very different lineages, i.e. independently of other variation at the locus. For example, an E92K mutation has been shown to cause constitutive activation (high activity in the absence of MSH) in natural MC1R variants of mouse, chicken and red-ruffed lemur (Varecia rubra) [19,22,12] as well as human MC1R subjected to in vitro mutagenesis [23] Another example is the R65C mutation, which leads to low MC1R activity in beach mice (Peromyscus polionotus) and mammoth [24,25]. Taken together, the presence of histidine at position 230 is therefore likely to affect MC1R activity in C. atratus and Coscoroba.

We therefore suggest that the most likely explanation for our results is that MC1R activity is increased in both C. atratus and Coscoroba, but whereas in C. atratus this leads to the expected melanism, in Coscoroba other genetic variation, e.g. at one or more loci downstream of MC1R, influences the phenotype. A precedent for the situation in Coscoroba is found in the red-ruffed lemur (Varecia rubra). This primate has a large amount of phaeomelanin in the coat but also an E94K mutation (homologous to position 92 in the swans), which is generally associated with large amounts of eumelanin [26]. Recent in vitro studies have shown that this mutation leads to the expected constitutive activation in MC1R, but that in this case the MC1R is still responsive to the inhibitor ASIP, so regulation of ASIP is almost certainly involved in the phenotype [12].

Reconstruction of Cygnus phylogeny

This is the first time to our knowledge that the relationships within the Cygnus genus have been investigated at the molecular level. Surprisingly, our results do not support the monophyly of Cygnus but suggest that C. melanocoryphus is more closely related to the emperor goose than to the remaining Cygnus species. These results conflict with a consensus parsimony tree obtained from 165 morphological traits in modern and fossil Anserinae [27] in which Cygnus was monophyletic, with C. melanocoryphus a sister taxon to C. atratus. Interestingly, C. melanocoryphus was first classified, not as Cygnus, but as the sole member of the genus Sthenelides [28]. However, the results are based on two loci, and resolution is not strong throughout the phylogeny, and further molecular data and denser sampling of other true geese (Anser and Branta) will be required to increase the resolution of the Cygnus phylogeny. However, it is important to note that the placement of C. melanocoryphus does not affect the conclusions of the MC1R analysis.

Elsewhere on the phylogeny, the coscoroba swan is removed from the Cygnus and emperor goose clade, and instead forms a strongly supported clade with the Cape Barren goose, a grouping that has received previous support in molecular datasets (CR1 transposon in LDH B [17], mitochondrial srRNA [29]). In contrast, in a recent analysis of the Anseriformes based on mtDNA control region sequences, a clade containing Cygnus, Coscoroba and Cereopsis to the exclusion of Anser was recovered, although statistical support for this grouping was rather weak [15].

Conclusion

We have likely uncovered an example where macroevolutionary change in plumage colour among species of Cygnus involves variation at the same locus that is involved in microevolution of melanism within a species of Anser, the snow goose A. caerulescens. However, the finding of a predicted melanizing mutation in MC1R in Coscoroba highlights the problems of providing a clear demonstration of functional involvement of a locus in a phenotype in cases where crosses/transgenics are not possible. A particular difficulty is that even though there is good information on the structure-function relationships of MC1R, and growing evidence that point mutations have effects on function that are largely independent of other sequence variation, the biochemical properties of MC1R do not necessarily predict the overall colour phenotype with epistatic interactions with other pigmentation loci also being key.

Methods

Sampling and DNA extraction

Samples comprised: muscle stored in ethanol from Cygnus olor (mute swan), C. cygnus (whooper swan), C. columbianus (tundra swan), Coscoroba coscoroba (coscoroba swan), Anser canagica (emperor goose) and Cereopsis novaehollandiae (cape barren goose) (National Museum of Scotland, Edinburgh); muscle stored in RNAlater (Ambion) from C. melanocoryphus (Natural History Museum, London); genomic DNA from five C. atratus individuals (Raoul Mulder, University of Melbourne). Genomic DNA was extracted using the QIAamp DNA Mini Kit (Qiagen) following the manufacturer's protocol.

Genomic DNA walking and sequencing of MC1R

The 3' UTR of MC1R from C. atratus was obtained using the BD GenomeWalker™ Universal Kit (BD Biosciences Clontech) according to manufacturer's instructions. Briefly, 2.5 μg of genomic DNA was digested with each of the four restriction enzymes; DraI, EcoRV, PvuII and StuI and ligated to the BD genomewalker adaptor. Primary PCRs were performed using the adaptor AP primer (kit) and 3' Swan O (5'-CCAACCCCTTCTGCAACTGCTTCTTCGG-3') obtained from preliminary partial MC1R sequences of C. atratus. PCR reactions were performed in 50 μl, with 2.5 U polymerase (Extensor Hi-Fidelity PCR Master Mix: Abgene), 0.4 μM of each primer and 1 μl adapted DNA. The cycling conditions were: 7 × 94°C 25 s, 72°C 3 min.; 32 × 94°C 25 s, 67°C 3 min.; 67°C 7 min. Secondary PCRs were performed using primers AP2 (kit) and 3'swan inner (5'-CAACCTCTTCCTCATCCTCATCATCTGC-3'). A 1.3 kb segment of MC1R including the entire coding region was amplified using a forward primer designed to the 5' UTR of the chicken MC1R gene (Chick 5' utr; 5'-GAGAAAGGGCCCTTTCTTC-3') and a reverse primer designed to the isolated 3' UTR of the C. atratus MC1R gene (Swan 3' utr; 5'-TAGCCTTTATTCGGTACCG-3') in a PCR containing 2.5U taq (Extensor Hi-Fidelity PCR Master Mix, Abgene), 2.25 mM MgCl2, 350 μM dNTPs, 0.2 μM of each primer and 50 ng DNA with cycling conditions of 94°C 3 min, 40 × 94°C 30 s, 62°C 45 s, 70°C 45 s; 70°C 5 min. For some samples, a new 5' UTR primer (Swan 5' utr; 5'-GCCCATGTCCTCTTGACC-3') was designed that bound downstream of the chicken 5' primer. PCR products were sequenced on both strands with BigDye v. 3.1 (Applied Biosystems) using PCR primers and the internal primers MSHR73 (5'-GGCGTAGAAGATGGTGATGTAGC-3') and MSHR74 (5'-GTGGACCGCTACATCACCAT-3').

Cytochrome B (CytB) and Lactate Dehydrogenase B (LDH B) sequencing

PCR primers Cyt B+ (5'-TGCCGACCTGCAAGTAGC-3'), Cyt B- (5'-GCGATTGAGGCGAGTTGG-3'), LdhF (5'-GGAAGACAAACTAAAAGGAGAAATGATGGA-3') and LdhR (5'-TTCCTCTGAAGCAGGTTGAGACGACTCTC-3') were designed to conserved regions identified in avian CytB and LDH B sequences available on GenBank. PCRs were performed with standard Taq and reagent concentrations. For Cytb, PCR conditions were: 94°C 3 min.; 35 × 94°C 30 s, 69°C 30 s, 70°C 45 s; 70°C 5 min. A product of approximately 1 kb was sequenced with the PCR primers and the two internal primers CytB 200-(5'-AGGCGAGTGAGGTGTCTG-3') and CytB 700-(5'-ATTTTGTCGCAGTCTGAYACG-3'). Sequences of 502 bp were obtained for all samples. For LDH B, PCR cycling conditions were: 94°C 2 min.; 35 × 94°C 1 min, 55°C 1 min, 72°C 2 min.; 72°C 5 min. The 600 bp products were sequenced with the PCR primers and the internal primers Ldh 390R (5'-AACTGCCACCCAAAAYAG-3'), Ldh 450F (5'-TAGTACTTTAATATAGAAG-3') and Ldh 750R (5'-TGACGAACACCTGCAGTTAC-3') providing sequences of 310 bp for all samples.

(CytB, LDH B and MC1R sequences have been deposited in GenBank, Accession numbers FJ170040FJ170063).

Phylogenetic analysis

Sequences of CytB and LDH B were concatenated and aligned with CLUSTAL-W [30] using MegAlign software (DNAstar). Positions containing gaps were removed before further analysis. Using MODELTEST [31], the model of sequence evolution with the highest maximum likelihood score as shown by the Akaike Information Criterion method was obtained and then implemented in maximum likelihood analyses in PAUP (v. 3.1, [32]). Branch support was examined by bootstrap resampling [33] with 100 replicates. Bayesian inference of phylogeny was performed in MrBayes (v. 3.1.2, [34]. The data were partitioned into CytB and LDH B and a GTR + I + Γ model of evolution applied to each partition, and 500,000 generations run, with sampling every 100 generations. Reconstruction of MC1R amino acid substitutions over the phylogeny obtained was performed in PAML (v.3.13, [35]).

Tests for Selection

For the MC1R dataset, dN/dS ratios were estimated in PAML using a codon-based substitution model. Two separate analyses were performed: 1). Lineage-specific models where the dN/dS was estimated for lineages with and without melanism, i.e. C. atratus and C. melanocoryphus to the rest of the species analysed. 2). Site specific models in which different codons are grouped into dN/dS categories across all lineages. Likelihood ratio tests were performed on nested sets of analyses.

Authors' contributions

MAP carried out all lab work and performed MODELTEST, PAUP and PAML analyses. NIM conceived of the study and performed the MrBayes analysis. Both authors were involved in the experimental design and wrote and approved the final manuscript.

Acknowledgements

We thank the National Museum of Scotland, Edinburgh and the Natural History Museum, London for providing muscle samples for gDNA extraction and Raoul Mulder, University of Melbourne for gDNA samples from Cygnus atratus. We also thank Adrian Friday and Brenda Bradley for helpful comments during manuscript preparation and George Speller for the swan and geese drawings used in Figure 1. This work was supported by NERC and the Leverhume Trust.

References

  1. Hoekstra HE: Genetics, development and evolution of adaptive pigmentation in vertebrates.

    Heredity 2006, 97:222-234. PubMed Abstract | Publisher Full Text OpenURL

  2. Jackson IJ: Homologous pigmentation mutations in human, mouse and other model organisms.

    Hum Mol Genet 1997, 6:1613-1624. PubMed Abstract | Publisher Full Text OpenURL

  3. Doucet SM, Shawkey MD, Rathburn MK, Mays HL Jr, Montgomerie R: Concordant evolution of plumage colour, feather microstructure, and a melanocortin receptor gene between mainland and island populations of a fairy-wren.

    Proc R Soc B 2004, 271:1663-1670. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. 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-453. PubMed Abstract | Publisher Full Text OpenURL

  5. Hoekstra HE, Nachman MW: Different genes underlie adaptive melanism in different populations of rock pocket mice.

    Mol Ecol 2003, 12(5):1185-1194. PubMed Abstract | Publisher Full Text OpenURL

  6. Mundy NI, Badcock NS, Hart T, Scribner K, Jansen K, Nadeau NJ: Conserved genetic basis of a quantitative trait involved in mate choice.

    Science 2004, 303:1870-1873. PubMed Abstract | Publisher Full Text OpenURL

  7. Rosenblum EB, Hoekstra HE, Nachman MW: Adaptive color variation and the evolution of the MC1R gene.

    Evolution Int J Org Evolution 2004, 58:1794-1808. OpenURL

  8. Theron E, Hawkins K, Bermingham E, Ricklefs RE, Mundy NI: The molecular basis of an avian plumage polymorphism in the wild: a point mutation in the melanocortin-1 receptor is perfectly associated with melanism in the bananaquit (Coereba flaveola).

    Curr Biol 2001, 11:550-557. PubMed Abstract | Publisher Full Text OpenURL

  9. Ritland K, Newton C, Marshall HD: Inheritance and population structure of the white-phased "Kermode" black bear.

    Curr Biol 2001, 11:1468-1472. PubMed Abstract | Publisher Full Text OpenURL

  10. Baião PC, Schreiber E, Parker PG: The genetic basis of the plumage polymorphism in red-footed boobies (Sula sula): a Melanocortin-1 receptor (MC1R) analysis.

    J Heredity 2007, 98:287-292. Publisher Full Text OpenURL

  11. Macdougall-Shackleton EA, Blanchard L, Igdoura SA, Gibbs HL: Unmelanized plumage patterns in Old World leaf warblers do not correspond to sequence variation at the melanocortin-1 receptor locus (MC1R).

    Mol Biol Evol 2003, 20:1675-1681. PubMed Abstract | Publisher Full Text OpenURL

  12. Haitina T, Ringholm A, Kelly J, Mundy NI, Schiöth HB: High diversity in functional properties of melanocortin -1 receptor (MC1R) in divergent primate species is more strongly associated with phylogeny than coat colour.

    Mol Biol Evol 2007, 24:2001-2008. PubMed Abstract | Publisher Full Text OpenURL

  13. Mundy NI: A window on the genetics of evolution: MC1R and plumage colouration in birds.

    Proc Biol Sci 2005, 272:1663-1640. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Nadeau NJ, Burke T, Mundy NI: Evolution of an avian pigmentation gene correlates with a measure of sexual selection.

    Proc R Soc B 2007, 274:1807-1813. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Donne-Goussé C, Laudet V, Hänni C: A molecular phylogeny of anseriformes based on mitochondrial DNA analysis.

    Mol Phylogenet Evol 2002, 23:339-356. PubMed Abstract | Publisher Full Text OpenURL

  16. Tateno Y, Takezaki N, Nei M: Relative efficiencies of the maximum-likelihood neighbour-joining, and maximum-parsimony methods when substitution rate varies with site.

    Mol Biol Evol 1994, 11:261-277. PubMed Abstract | Publisher Full Text OpenURL

  17. St John J, Cotter JP, Quinn TW: A recent chicken repeat 1 retrotransposition confirms the Coscoroba-Cape Barren goose clade.

    Mol Phylogenet Evol 2005, 37:83-90. PubMed Abstract | Publisher Full Text OpenURL

  18. Nachman MW, Hoekstra HE, D'Agostino SL: The genetic basis of adaptive melanism in pocket mice.

    Proc Natl Acad Sci USA 2003, 100:5268-5273. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Robbins LS, Nadeau JH, Johnson KR, Kelly MA, Roselli-Rehfuss L, Baack E, Mountjoy KG, Cone RD: Pigmentation phenotypes of variant extension locus alleles result from point mutations that alter MSH receptor function.

    Cell 1993, 72:827-834. PubMed Abstract | Publisher Full Text OpenURL

  20. Takeuchi S, Suzuki H, Hirose S, Yabuuchi M, Sato C, Yamamoto H, Takahashi S: Molecular cloning and sequence analysis of the chick melanocortin 1-receptor gene.

    Biochim Biophys Acta 1996, 1306:122-126. PubMed Abstract OpenURL

  21. Nadeau NJ, Minvielle F, Mundy NI: Association of a Glu92Lys substitution in MC1R with extended brown in Japanese quail (Coturnix japonica).

    Anim Genet 2006, 37:287-289. PubMed Abstract | Publisher Full Text OpenURL

  22. Ling MK, Lagerström MC, Fredriksson R, Okimoto R, Mundy NI, Takeuchi S, Schiöth HB: Association of feather colour with constitutively active melanocortin 1 receptors in chicken.

    Eur J Biochem 2003, 270:1441-1449. PubMed Abstract | Publisher Full Text OpenURL

  23. Schiöth HB, Phillips SR, Rudzish R, Birch-Machin MA, Wikberg JES, Rees JL: Loss of function mutations of the human melanocortin 1 receptor are common and are associated with red hair.

    Biochem Biophys Res Commun 1999, 260:488-491. PubMed Abstract | Publisher Full Text OpenURL

  24. Hoekstra HE, Hirshmann RJ, Bundey RA, Insel PA, Crossland JP: A single amino acid mutation contributes to adaptive beach mouse color pattern.

    Science 2006, 313:101-104. PubMed Abstract | Publisher Full Text OpenURL

  25. Römpler H, Rohland N, Lalueza-Fox C, Willerslev E, Kuznetsova T, Rabeder G, Bertranpetit J, Schöneberg T, Hofreiter M: Nuclear gene indicates coat-color polymorphism in mammoths.

    Science 2006, 313:62. PubMed Abstract | Publisher Full Text OpenURL

  26. Mundy NI, Kelly J: Evolution of a pigmentation gene, the melanocortin-1 receptor in primates.

    Am J Phys Anthropol 2003, 121:67-80. PubMed Abstract | Publisher Full Text OpenURL

  27. Livezey BC: A phylogenetic analysis of Geese and Swans (Anseriformes: Anserinae), including selected fossil species.

    Syst Biol 1996, 45:415-450. Publisher Full Text OpenURL

  28. Kingsley JS: The Standard Natural History. Boston S.E. Cassino & Co; 1885. OpenURL

  29. Zimmer R, Erdtmann B, Thomas WK, Quinn TW: Phylogenetic analysis of the Coscoroba coscoroba using mitochondrial srRNA gene sequences.

    Mol Phylogenet Evol 1994, 3:85-91. PubMed Abstract | Publisher Full Text OpenURL

  30. Thompson JD, Higgins DG, Gigson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting position-specific gap penalties and weight matrix choice.

    Nucleic Acid Res 1994, 22:4673-4680. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution.

    Bioinformatics 1998, 14:817-818. PubMed Abstract | Publisher Full Text OpenURL

  32. Swofford DL: Phylogenetic Analysis Using Parsimony. In Version 3.1 edition. Illinois Natural History Survey, Champaign Illinois; 1991. OpenURL

  33. Felsenstein J: Phylogenies and the comparative method.

    Am Nat 1985, 125:1-15. Publisher Full Text OpenURL

  34. Huelsenbeck JP, Rohnquist F: MrBayes: Bayesian inference of phylogenetic trees.

    Bioinformatics 2001, 17:754-755. PubMed Abstract | Publisher Full Text OpenURL

  35. Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood.

    Comput Appl Biosci 1997, 13(5):555-556. PubMed Abstract OpenURL