Testing whether macroevolution follows microevolution: Are colour differences among swans (Cygnus) attributable to variation at the MC1R locus?
© Pointer and Mundy; licensee BioMed Central Ltd. 2008
Received: 06 May 2008
Accepted: 12 September 2008
Published: 12 September 2008
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).
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.
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.
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 . 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 . 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 . Interestingly, the same mutation was recently found to be associated with melanism in another bird, the red-footed booby (Sula sula) .
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 , 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  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 .
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 
Phylogenetic reconstruction of swans and geese using CytB and LDH Bdata
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  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 .
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)  and a histidine at this position is also associated with melanism in rock pocket mice (Chaetodipus intermedius) .
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 .
MC1Rvariation 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).
MC1Rvariation 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  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 . 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 .
Reconstruction of Cygnusphylogeny
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  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 . 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 , mitochondrial srRNA ). 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 .
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.
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 FJ170040–FJ170063).
Sequences of CytB and LDH B were concatenated and aligned with CLUSTAL-W  using MegAlign software (DNAstar). Positions containing gaps were removed before further analysis. Using MODELTEST , 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, ). Branch support was examined by bootstrap resampling  with 100 replicates. Bayesian inference of phylogeny was performed in MrBayes (v. 3.1.2, . 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, ).
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.
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.
- Hoekstra HE: Genetics, development and evolution of adaptive pigmentation in vertebrates. Heredity. 2006, 97: 222-234. 10.1038/sj.hdy.6800861.View ArticlePubMedGoogle Scholar
- Jackson IJ: Homologous pigmentation mutations in human, mouse and other model organisms. Hum Mol Genet. 1997, 6: 1613-1624. 10.1093/hmg/6.10.1613.View ArticlePubMedGoogle Scholar
- Doucet SM, Shawkey MD, Rathburn MK, Mays HL, 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. 10.1098/rspb.2004.2779.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1016/S0960-9822(03)00128-3.View ArticlePubMedGoogle Scholar
- Hoekstra HE, Nachman MW: Different genes underlie adaptive melanism in different populations of rock pocket mice. Mol Ecol. 2003, 12 (5): 1185-1194. 10.1046/j.1365-294X.2003.01788.x.View ArticlePubMedGoogle Scholar
- 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. 10.1126/science.1093834.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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. 10.1016/S0960-9822(01)00158-0.View ArticlePubMedGoogle Scholar
- Ritland K, Newton C, Marshall HD: Inheritance and population structure of the white-phased "Kermode" black bear. Curr Biol. 2001, 11: 1468-1472. 10.1016/S0960-9822(01)00448-1.View ArticlePubMedGoogle Scholar
- 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. 10.1093/jhered/esm030.View ArticleGoogle Scholar
- 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. 10.1093/molbev/msg186.View ArticlePubMedGoogle Scholar
- 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. 10.1093/molbev/msm134.View ArticlePubMedGoogle Scholar
- Mundy NI: A window on the genetics of evolution: MC1R and plumage colouration in birds. Proc Biol Sci. 2005, 272: 1663-1640. 10.1098/rspb.2005.3107.View ArticleGoogle Scholar
- 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. 10.1098/rspb.2007.0174.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1016/S1055-7903(02)00019-2.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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. 10.1016/j.ympev.2005.03.005.View ArticlePubMedGoogle Scholar
- 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. 10.1073/pnas.0431157100.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1016/0092-8674(93)90572-8.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1365-2052.2006.01442.x.View ArticlePubMedGoogle Scholar
- 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. 10.1046/j.1432-1033.2003.03506.x.View ArticlePubMedGoogle Scholar
- 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. 10.1006/bbrc.1999.0935.View ArticlePubMedGoogle Scholar
- 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. 10.1126/science.1126121.View ArticlePubMedGoogle Scholar
- 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-10.1126/science.1128994.View ArticlePubMedGoogle Scholar
- Mundy NI, Kelly J: Evolution of a pigmentation gene, the melanocortin-1 receptor in primates. Am J Phys Anthropol. 2003, 121: 67-80. 10.1002/ajpa.10169.View ArticlePubMedGoogle Scholar
- Livezey BC: A phylogenetic analysis of Geese and Swans (Anseriformes: Anserinae), including selected fossil species. Syst Biol. 1996, 45: 415-450. 10.2307/2413524.View ArticleGoogle Scholar
- Kingsley JS: The Standard Natural History. 1885, Boston S.E. Cassino & CoGoogle Scholar
- 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. 10.1006/mpev.1994.1011.View ArticlePubMedGoogle Scholar
- 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. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Swofford DL: Phylogenetic Analysis Using Parsimony. Version 3.1 edition. 1991, Illinois Natural History Survey, Champaign IllinoisGoogle Scholar
- Felsenstein J: Phylogenies and the comparative method. Am Nat. 1985, 125: 1-15. 10.1086/284325.View ArticleGoogle Scholar
- Huelsenbeck JP, Rohnquist F: MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13 (5): 555-556.PubMedGoogle Scholar
This article is published under license to 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.