Genetic variation at hair length candidate genes in elephants and the extinct woolly mammoth
BMC Evolutionary Biology volume 9, Article number: 232 (2009)
Like humans, the living elephants are unusual among mammals in being sparsely covered with hair. Relative to extant elephants, the extinct woolly mammoth, Mammuthus primigenius, had a dense hair cover and extremely long hair, which likely were adaptations to its subarctic habitat. The fibroblast growth factor 5 (FGF5) gene affects hair length in a diverse set of mammalian species. Mutations in FGF5 lead to recessive long hair phenotypes in mice, dogs, and cats; and the gene has been implicated in hair length variation in rabbits. Thus, FGF5 represents a leading candidate gene for the phenotypic differences in hair length notable between extant elephants and the woolly mammoth. We therefore sequenced the three exons (except for the 3' UTR) and a portion of the promoter of FGF5 from the living elephantid species (Asian, African savanna and African forest elephants) and, using protocols for ancient DNA, from a woolly mammoth.
Between the extant elephants and the mammoth, two single base substitutions were observed in FGF5, neither of which alters the amino acid sequence. Modeling of the protein structure suggests that the elephantid proteins fold similarly to the human FGF5 protein. Bioinformatics analyses and DNA sequencing of another locus that has been implicated in hair cover in humans, type I hair keratin pseudogene (KRTHAP1), also yielded negative results. Interestingly, KRTHAP1 is a pseudogene in elephantids as in humans (although fully functional in non-human primates).
The data suggest that the coding sequence of the FGF5 gene is not the critical determinant of hair length differences among elephantids. The results are discussed in the context of hairlessness among mammals and in terms of the potential impact of large body size, subarctic conditions, and an aquatic ancestor on hair cover in the Proboscidea.
Hair is a defining characteristic of mammals. The hair follicle is the only organ in mammals to undergo life-long cycles of growth, regression and quiescence . Hair development proceeds through a cycle of anagen in which hair follicles undergo rapid growth, catagen in which hair growth ceases due to apoptosis-driven regression, and telogen in which the hair follicle enters a period of relative quiescence . Longer hair can thus result from an increase in the length of time during which anagen proceeds. A loss-of-function mutation in the fibroblast growth factor 5 gene (designated Fgf5 in mice and rats, and FGF5 in other mammals) is responsible for the long hair phenotype present in angora mice, while a similar long hair phenotype occurs in mice homozygous for a null allele of Fgf5 produced by gene targeting . In mice, catagen does eventually occur even in the absence of functional FGF5, indicating that other factors are also involved in the cycle .
The FGF5 gene consists of three exons . In wild-type mice and other mammals, the FGF5 transcript is present in two isoforms, with the smaller transcript due to alternative splicing in which exon 2 is excluded from the mRNA [3, 4]. The shorter transcript antagonizes the activity of the longer transcript, suggesting that they function together in hair cycle regulation . In mice, the Fgf5 mutation causing the long-hair angora phenotype affects exon 1 . In dogs, sequencing of the FGF5 gene in 218 individuals from 14 breeds, including three dog breeds fixed for long hair and five breeds fixed for short hair, identified a missense mutation in exon 1 as responsible for the long haired phenotype . In domestic cats, the FGF5 gene was shown to be associated with hair length [6, 7], with four independent mutations in exon 1 or 3 considered to be functionally significant in controlling hair length in a survey of more than 380 individuals from 26 short- or long-haired breeds, non-breed cats and two pedigrees . In rabbits, mutations in FGF5 have been reported to have significant association with wool yield . Given that the species for which FGF5 is known to influence hair length belong to different superordinal placental clades that diverged ca. 97 million years ago (Mya) (Figure 1) [9, 10], it seems plausible to hypothesize that FGF5 may be a critical determinant of hair length across mammals.
The Elephantidae, a family of proboscideans comprising the living elephants, their extinct relatives, and the extinct mammoths , would constitute an important group for the comparative study of genes involved in regulating hair cover and growth. The woolly mammoth, Mammuthus primigenius, was covered with hair ranging in length from a few centimeters to over 90 cm, with coarse outer hairs, beneath which were shorter, thinner hairs forming densely packed underwool (maximum 2.5-8 cm long) that formed a thermal insulating layer [12, 13]. By contrast, extant elephants are sparsely covered with hair and have hair of short length . Given the role of FGF5 in determining hair length in a diverse set of mammalian taxa, we hypothesized that loss of function of this gene may play a role in the longer hair length of the woolly mammoth. We therefore generated and compared sequences of FGF5 from living elephants and the woolly mammoth. In addition, using a combination of bioinformatics analysis and DNA sequencing, another gene related to hair phenotype in humans, type I hair keratin pseudogene (KRTHAP1) , was examined, as were three other genes coding for hair keratin proteins .
The FGF5 gene was sequenced from a woolly mammoth and from two Asian elephants (Elephas maximus), two African savanna elephants (Loxodonta africana) and two African forest elephants (L. cyclotis) (Table 1) . For the woolly mammoth, protocols established for ancient DNA were used , with the complete FGF5 coding sequence obtained for Indigirka mammoth N2031, a ca. 11000-13000 year-old tooth from the Indigirka River basin, Federal Republic of Russia. PCR products from the Indigirka mammoth were amplified from multiple extracts at two different laboratories (Norfolk and Thunder Bay), with 2 or more PCRs performed per fragment, and PCR fragments cloned and sequenced. Among-clone variation was observed but no consistent differences among amplifications were detected. Partial sequence was also obtained from the Jarkov mammoth, a ca. 20,380 year-old sample from the Taimyr Peninsula, Federal Republic of Russia. (See Methods and Additional file 1 for further information on the mammoths and laboratory protocols and primers).
The complete 5' UTR, complete open reading frames (ORFs) of the three exons, and part of the promoter region of the FGF5 gene were sequenced in the elephants and the Indigirka woolly mammoth (Table 1). All of the elephants as well as the woolly mammoth had an uninterrupted open reading frame (without premature stop codons). In each of the elephantids, exon 1 was 593 bp in length (229 bp of 5' UTR, 364 bp of coding region); exon 2 was 104 bp; while the protein coding region of exon 3 was 348 bp, with 40 bp of the 3' UTR also sequenced. The four boundaries between exons and introns were identical across all elephantids sequenced, thus FGF5 does not vary at these splice sites among elephantids. Full sequencing of the two introns was not attempted due to their length: 7,729 bp and 11,312 bp for introns 1 and 2, respectively, in human, with even longer introns present in savanna elephant based on genomic traces (data not shown). Additionally, the mammoth genomic sequences  were found to have poor coverage of both introns in the mammoth (data not shown).
Only two mammoth-specific differences were found: one in the promoter and one in exon 3 (Table 1). The substitution in the promoter sequence (Table 1) did not alter any predicted transcription factor binding sites . The difference present in exon 3 of the mammoth was a guanine to adenine substitution at position 790 that was a silent mutation, i.e. did not alter the amino acid sequence (Table 1). This substitution also did not lead to a rare codon being present in the mammoth. The mammoth FGF5 amino acid sequence was identical to that of 2 of the 3 extant elephant species. This suggested that all elephantids including the long-haired mammoths had a functional FGF5 protein.
The nucleotide sequence of FGF5 varied among living elephantids (Table 1), although only one non-synonymous substitution was found among living or extinct elephantids. In both forest elephant individuals, exon 1 (nucleotide position 327) contained a codon for glycine at residue 33, whereas a codon for alanine was present in other elephantids at this position (Table 1, Figure 2). This is a physicochemically conservative change . Among other mammals for which FGF5 sequence is available, this G33A mutation was found to be present only in the in rodent FGF5 protein sequence (Figure 2). The mutation is located at the N-terminus of the protein (Figure 2). The N- and C-termini of other members of the FGF family play key roles in the specificity of interaction with the FGF-receptors (FGFRs). . However, this mutation resides within a mainly unstructured region, predicted to be an extended loop downstream of the signal peptide (Figure 3A and data not shown). Thus, the G33A mutation in forest elephants would be unlikely to affect the secondary structure of the FGF5 protein. Both glycine and alanine are nonpolar, neutral amino acids, with neighboring hydropathy indices [21, 22] in the hydrophobic range, which precludes prediction of their structural position, i.e. external or internal. The tertiary structure of this region could not be modeled, because the structure of this region has not been experimentally resolved in any member of the FGF family. Analysis using SIFT and POLYPHEN programs suggested that while the G33A mutation may increase the stability of the protein, the G33A mutation in forest elephants should not have serious consequences on FGF5 protein function (data not shown).
Human and mammoth FGF5 protein sequences were also compared [23, 24]. FGF5 is predicted to be secreted and has an almost identical signal peptide sequence for the two species (Figure 3A). The major difference between the two sequences is an insertion/deletion of three amino acids at the N-terminal region. This region is predicted to include many O-glucosylation sites, suggesting a putative difference in glycosylation between the two FGF5 proteins (Figure 3A). The two proteins are predicted to have a unique N-glycosylation site at position 110 of the human sequence.
The general features of the amino acids of mammalian FGF5 are shown in Additional file 1, as are the phylogenetic relationships among FGF5 amino acid sequences. The FGF5 proteins are very similar across elephantids but differ from those of other mammalian FGF5 proteins. To test whether any of these differences would be predicted to alter the three-dimensional structure of the elephantid FGF5 proteins, the tertiary structures of FGF5 in different species were predicted. This analysis predicted that all described mammalian FGF5 proteins fold similarly (data not shown). The structural analysis also revealed that the amino acid differences between human and mammoth FGF5 sequences (shown in green color in Figure 3B) do not correspond to residues known to interact with the FGF receptor (FGFR, yellow color) and heparin (blue color) (Figure 3B). The two amino acid differences between humans and mammoths that are included in the 3D model are predicted to be parts of loops and do not seem to affect the secondary or the tertiary conformation of the FGF5 molecule (Figure 3B).
To examine the quality of sequence traces for the recently published genome sequence of the woolly mammoth , the mammoth FGF5 DNA sequences generated for this study were compared using BLAST to homologous sequences generated by the Mammoth Genome Project http://mammoth.psu.edu . Five matching mammoth genome sequences were found, comprising sequence coverage of about 50% (714/1434 bp). Coverage of the mammoth genome varied by region for the FGF5 gene. Four genomic sequences matched the promoter and 5' UTR; these covered 99% (513/520 bp) of the corresponding region sequenced by the current study. There were 6 discrepancies among the traces covering this region. Only one genomic sequence was found that overlapped with the coding regions, and it covered only 22% (201/915 bp) of the sequence determined for the current study, with nine discrepancies found between genomic traces and our sequences. Overall, four of the five mammoth genomic sequences had discrepancies with the mammoth sequences generated for the current study, with a total of 15 nucleotide site discrepancies detected. The discrepancies relative to our sequence likely reflect damage present in the ancient DNA of the mammoths used to generate genomic sequences. Similar ancient DNA damage affected the mammoth sequences generated for the current study (although for the current study multiple clones from at least two independent PCRs per fragment were used to successfully generate a consensus sequence). For five PCR amplicons used in the current study to determine the sequence of the FGF5 promoter, the among-clone diversity was in the range 0-9 (see Additional file 1). Unlike PCR-based approaches where multiple PCRs can be performed and multiple clonal sequences per PCR determined to generate a consensus sequence, the mammoth genome does not currently have high enough coverage per base to be confident that observed differences among traces, individuals or species represent the true sequence rather than ancient DNA damage. Thus although the mammoth genome is extremely useful for designing mammoth-specific primers and for initial queries, our data suggest that PCR, cloning and sequencing would still be required to determine mammoth DNA sequences, to account both for ancient DNA damage and gaps in the low-coverage genome sequences.
Like FGF5, other loci have been identified that are associated with reduced hair cover. In humans the type I hair keratin pseudogene KRTHAP1 has a premature stop codon in the fourth exon, and protein is not detected in human hair follicles . In great apes, the orthologous gene has an intact ORF, with RNA expressed and protein translated in the hair follicles of chimpanzees (cHaA) and gorillas (gHaA) . Thus, while closely related primates with dense hair coverage express this gene, relatively hairless humans do not. Using the Loxodonta africana draft genome sequence, all of the homologous exons except for exon 7 for this gene were identified. Exon 1 displayed a predicted premature stop codon (Figure 4). Thus, as in humans, this gene appeared to be disrupted in the savanna African elephant. A 302 bp segment of exon 1 was therefore amplified and sequenced from the Indigirka mammoth to examine the region that contained the premature stop codon in the elephant. The stop codon was found to be present in the mammoth as well (Figure 4), suggesting that this mutation is not involved in hair phenotype differences among elephantids.
Among mammals, hair phenotype is affected by hair keratin genes . We therefore examined keratin genes reported as displaying either elephant or mammoth unique differences , using genomic sequences of savanna elephant or woolly mammoth [15, 26]. KRT25 was identified as having a unique alanine to serine change, but this was specific to only one of two woolly mammoths previously sequenced . Similarly, in the elephantids KRT27 and KRT83 were found to code for rare amino acid differences. However, only a methionine to valine change in KRT27 was found to be unique to mammoths, while the methionine present in elephant KRT27 was also present in the fully hair-covered hyrax. Thus the differences in KRT27 and KRT83 are unlikely to be associated with differences in hair cover.
To date, the FGF5 mutations found to produce long hair have been uncovered in phenotypic variants among laboratory or domestic mammals, including mice, rabbits, dogs and cats (Figure 1) [2, 5, 6, 8]. A role for FGF5 in inter- as opposed to intra-species differences in hair length has not been established. Nonetheless, the association of FGF5 mutations with long-hair phenotypes in a wide variety of distantly related mammals (Figure 1) suggested that FGF5 might be a determinant of hair length in mammals in general. To test this hypothesis we sequenced the open reading frames of all three exons, the 5' UTR and the promoter region of the FGF5 gene in the relatively hairless extant elephantids and in the woolly mammoth. Our data show that these regions of FGF5 are highly conserved among elephantids, including the woolly mammoth. Only one variant in the amino acid sequence was detected among elephantids, the G33A mutation in forest elephants coded by exon 1. However, our analysis suggested that this mutation would not greatly affect protein function, a conclusion also supported by the presence of same amino acid substitution in the wild-type sequence of FGF5 in murid rodents. While regulatory mechanisms may exist that would not be detectable by our study, and a role for FGF5 in the long hair of mammoths cannot be completely ruled out, the most parsimonious interpretation of our results suggests that FGF5 was not the major genetic determinant of long hair in mammoth.
Similarly, no differences were found among elephantids for partial sequences of several additional candidate genes such as KRTHAP1, KRT25, KRT27, and KRT83. Thus, none of the candidate genes examined thus far demonstrated a clear difference exclusive to mammoths, which would be necessary for establishing a role in their unique dense and long-haired phenotype relative to extant elephants. While a host of additional genes are known to influence hair development, many play other critical developmental roles and would likely be lethal if function were perturbed . Thus, future candidate genes will likely reside among the keratin and keratin-associated protein (KRTAP) genes, believed to play a role in the evolution of mammalian hair characteristics . Among mammals, KRTAP gene repertoires vary considerably, with homogenization within groups , although the genes have not been catalogued in elephants or other afrotheres. Once the savanna elephant genome is complete, keratin and KRTAP genes from this species may be identified as candidates for determining the hair differences among elephantids.
Among living mammals hairlessness is more pronounced among fully aquatic species of sirenians and cetaceans (Figure 1); thus the designation of humans and elephants as "hairless" is a relative term . In the case of elephants, hairlessness may be a thermoregulatory adaptation to large body size , which would be consistent with a gain of hair cover for the woolly mammoth [12, 29], since mammoths appear first in Africa before the lineage adapted to colder environments . In considering the evolution and genetics of hair cover in extant elephants, woolly mammoths and other proboscideans, a number of factors must be taken into account. First, it is difficult to determine based on outgroups whether hair cover was lost in extant elephantids or gained in the woolly mammoth. The presence of considerable hair cover in a distantly related outgroup to the elephantids, the American mastodon (Mammut americanum) [30, 31], does not necessarily suggest that hair cover is ancestral. Hair cover in the American mastodon may comprise a convergent adaptation to cold and/or aquatic habitats, rather than an ancestral state . Second, the proboscidean lineage that gave rise to both elephantids and mastodons is likely to have derived from aquatic or semi-aquatic ancestors [32, 33]. Although many semi-aquatic species are not hairless , proboscideans derive from a common ancestor with the fully aquatic and hairless sirenians ; while Proboscidea and Sirenia, along with Hyracoidea (hyraxes, which are not hairless), comprise the Paenungulata, one of the few unresolved trichotomies among extant mammalian orders . Additionally, some ancestral proboscideans were as large as living elephants , and if hairlessness is an adaptation to large body size in terrestrial mammals, it may have been the ancestral state in proboscideans [12, 28, 29].
A third consideration is that, both in the case of the hyrax-sirenian-proboscidean clade and the elephantid clade, the evidence suggests that divergence of the ancestral line into two and then three descendent lineages occurred in quick succession. Among the elephantids, nearly complete mtDNA sequences have been generated for all three genera including mammoths [34–36]; using the mastodon mito-genome as an outgroup suggests that Loxodonta diverged from the common ancestor of Elephas and Mammuthus ca. 7.6 Mya; followed by the divergence of the two latter genera ca. 6.7 Mya (Figure 1) . The order of divergence among the Paenungulata remains unresolved , suggesting, as in the case of the elephantids, that the two divergences that yielded the three mammalian orders occurred in rapid succession. The rapid divergence of lineages in both cases suggests that incongruent lineage sorting of alleles may have affected many loci, causing discrepancies between gene and species trees [37–39]. Interestingly, the mammoth and African elephant FGF5 sequences are identical at positions -112, -150 and -269 of the promoter (Table 1), while differing from the Asian elephant sequence even though the Asian elephant and mammoth are sister taxa . This suggests that this region of the genome may have been subject to incongruent lineage sorting in which the gene tree does not match the species tree , as has been reported for other gene segments . Thus both convergent evolution and incongruent lineage sorting may have affected genes involved in hair cover among the Proboscidea.
Although the gene for long hair in mammoths was not here identified, proboscideans remain an important group for understanding the evolution of hair cover. While most mammals have dense hair cover, humans and extant elephants are notable in being relatively hairless , and both are closely related to species with much greater hair cover (great apes and woolly mammoths, respectively). Both lineages are also noteworthy in being "genome-enabled"  for the study of genes affecting hair cover. The human and chimpanzee genomes have been sequenced , while the elephant genome is being sequenced , and substantial coverage for the mammoth genome is now available . Other than aquatic species, the number of other mammalian genera considered to be "hairless" is quite small . Thus, for a comparative approach to the evolution of hair cover, proboscideans comprise an important group for further research.
Modern elephant DNA was extracted from blood or tissue samples. Wild African savanna elephants Laf-KR0014 and Laf-KR0138 were from Kruger National Park, South Africa. Wild African forest elephants Lcy-LO3505 and Lcy-LO3508 were from Lopé National Park in Gabon. Asian elephants Ema-6 and Ema-10 were zoo animals at the Rosamond Gifford Zoo at Burnet Park, Syracuse, NY. Both Ema-6 (North American studbook number 27) and Ema-10 (North American studbook number 28) had been wild-caught, most likely in Thailand.
The mammoth tooth designated N2031, which is the focus of this project, is from the Indigirka River basin, Russian Federation. N2031 was found in 1965 on the Berelekh river (a tributary of the Indigirka river), in the Berelekh mammoth "cemetery" in situ. The approximate geological age is 11000 - 13000 years before present (BP; G. Boeskorov, personal communication). The sample was originally obtained from the Geological Museum, Geological Institute, Yakutsk. Partial sequences were also obtained from the Jarkov mammoth discovered in the Taimyr Pensinsula, Russian Federation and dated to ca. 20,380 years BP . In order to obtain material for DNA extraction, an electric drill with individual sterile drill bits were used at low speed to collect the bone powder and shavings.
DNA extraction, PCR, and sequencing
Extractions of mammoth samples in Norfolk were carried out in a room dedicated to ancient DNA work in a CleanSpot PCR hood (Coy Laboratory) following an established protocol . Likewise, all pre-amplification work in Thunder Bay was performed in a 'Clean Lab'. PCR amplifications were performed at least twice per primer pair. Primer sequences and details of ancient DNA extractions, PCR and sequencing are included in Additional file 1. All PCR products were cloned and sequenced since direct sequencing can lead to an erroneous sequence due to contamination and DNA damage in the extract. Cloning and sub-sampling individual representative amplified sequences provides a better representation of the original template amplified , therefore none of the mammoth consensus sequences generated in this study were determined from direct sequencing.
DNA from extant elephants (~50 ng) underwent amplification by PCR using 200 nM final concentration of each oligonucleotide primer in 1.5 mM MgCl2, with AmpliTaq Gold DNA Polymerase (Applied Biosystems Inc. [ABI]). Primers are listed in Additional file 1. For all primer pairs, PCR consisted of an initial 95°C for 9:45 min; with cycles of 20 sec at 94°C, followed by 30 sec at 60°C (3 cycles); 58°C, 56°C, 54°C, or 52°C (5 cycles each temperature); or 50°C (last 22 cycles), followed by 30 sec extension at 72°C; with a final extension of 3 min at 72°C. PCR products were enzyme-purified  and sequenced using the BigDye Terminator v3.1 Cycle Sequencing Kit (ABI). Extension products were purified with Sephadex G-50 (Amersham), and resolved on an ABI 3730 DNA Analyzer. The software Sequencher 4.5 (Gene Codes Corp.) was used to edit chromatograms and assemble contigs. Gene identity was established by homology to GenBank entries with BLAST . Direct sequences for elephants and consensus sequences for mammoths generated for FGF5 and KRTHAP1 have been deposited in GenBank [GenBank:FJ755444-FJ755451].
Protein sequence analysis and structural prediction
Sequences were collected from the NCBI and the ENSEMBL databases using both keyword and homology searches. Multiple protein sequence alignments were performed using MAFFT 6 (E-INS-i algorithm; scoring matrix: BLOSUM 62; gap opening penalty: 1.53; gap extension penalty: 0.00) . Pairwise alignments were performed using the Smith-Waterman algorithm . N- and O-glycosylation sites were predicted using the NetNGlyc 1.0 and NetOGlyc 3.1 webservers http://www.cbs.dtu.dk/services . The signal peptides were predicted using the SignalP 3.0 server http://www.cbs.dtu.dk/services/SignalP . The effects of mutations on protein function were predicted using the SIFT  and POLYPHEN programs . Tests on protein stability and secondary structure predictions were performed using the MuPro http://www.ics.uci.edu/~baldig/mutation.html and SSPro8 http://scratch.proteomics.ics.uci.edu webservers [53, 54]. The PDB database http://www.pdb.org was searched to check whether the FGF5 structure has been experimentally resolved, but with negative results. For this reason, homology modeling and fold recognition were performed using the SWISS-MODEL http://swissmodel.expasy.org  and PHYRE http://www.sbg.bio.ic.ac.uk/~phyre  web servers. Both programs identified the human FGF9 structure [PDB:1IHK] as the best candidate (most similar; E-value = 10-45) to build a structural model of FGF5. Therefore, the mammoth, elephant and human FGF5 proteins were modeled by using the human FGF9 as template. Pairwise structural alignments and model structural superimposition was performed using the SSAP http://cathdb.info/cgi-bin/SsapServer.pl [57, 58] and DaliLite http://www.ebi.ac.uk/Tools/dalilite  webservers. Tertiary structure figures were generated using PyMol (DeLano Scientific; http://pymol.org).
FGF5 sequences from therian (placental and marsupial) mammals were aligned using MAFFT 6 (G-INS-i algorithm with JTT200 scoring matrix; gap opening penalty: 1.53; gap extension penalty: 0.00), and examined for residue variation using the FINGERPRINT web server http://evol.mcmaster.ca/fingerprint . The phylogenetic relationships among FGF5 sequences were examined in a maximum likelihood framework in RAxML 7.0.4  using the best-fit JTT protein substitution matrix  with empirical residue frequencies and among-site rate heterogeneity modeled with Γ with four classes , after comparing the log-likelihood of all substitution models available in RAxML.
Genome project sequences
Elephant sequences of KRTHAP1, KRT25, KRT27 AND KRT83 were identified in the NCBI Loxodonta africana genome Trace Archives http://www.ncbi.nlm.nih.gov/Traces/home using MegaBlast . Human or chimpanzee KRTHAP1 exon sequences were used as queries and obtained from GenBank [GenBank:AJ401054 and Y16795] or from the UCSC Genome Browser http://genome.ucsc.edu  (Human March 2006 [hg18] assembly). Elephant trace files obtained by matches to primates were themselves used as queries against the elephant genomic trace files, to obtain additional elephant sequences, with the process repeated to obtain further upstream and downstream elephant traces and sequences. Mammoth sequences were obtained from the mammoth genome project BLAST server http://mammoth.psu.edu . Mammoth sequences with a score above 100 were used. The mammoth and elephant sequences were also verified using a BLAT search  against human sequences on the UCSC website to verify the identity of the locus.
Transcription factor binding site and rare codon analyses
Transcription factor binding sites of promoter regions were predicted using TFSEARCH http://mbs.cbrc.jp/research/db/TFSEARCH.html that uses the TRANSFAC database . The tRNA effect of the guanine-to-adenine mammoth-specific nucleotide substitution was examined using RARE CODON CALTOR http://www.doe-mbi.ucla.edu/~sumchan/caltor.html.
Krause K, Foitzik K: Biology of the hair follicle: the basics. Semin Cutan Med Surg. 2006, 25 (1): 2-10. 10.1016/j.sder.2006.01.002.
Hebert JM, Rosenquist T, Gotz J, Martin GR: FGF5 as a regulator of the hair growth cycle: evidence from targeted and spontaneous mutations. Cell. 1994, 78 (6): 1017-1025. 10.1016/0092-8674(94)90276-3.
Hattori Y, Yamasaki M, Itoh N: The rat FGF-5 mRNA variant generated by alternative splicing encodes a novel truncated form of FGF-5. Biochim Biophys Acta. 1996, 1306 (1): 31-33.
Suzuki S, Ota Y, Ozawa K, Imamura T: Dual-mode regulation of hair growth cycle by two Fgf-5 gene products. J Invest Dermatol. 2000, 114 (3): 456-463. 10.1046/j.1523-1747.2000.00912.x.
Housley DJ, Venta PJ: The long and the short of it: evidence that FGF5 is a major determinant of canine 'hair'-itability. Anim Genet. 2006, 37 (4): 309-315. 10.1111/j.1365-2052.2006.01448.x.
Kehler JS, David VA, Schaffer AA, Bajema K, Eizirik E, Ryugo DK, Hannah SS, O'Brien SJ, Menotti-Raymond M: Four independent mutations in the feline fibroblast growth factor 5 gene determine the long-haired phenotype in domestic cats. J Hered. 2007, 98 (6): 555-566. 10.1093/jhered/esm072.
Drogemuller C, Rufenacht S, Wichert B, Leeb T: Mutations within the FGF5 gene are associated with hair length in cats. Anim Genet. 2007, 38 (3): 218-221. 10.1111/j.1365-2052.2007.01590.x.
Li CX, Jiang MS, Chen SY, Lai SJ: [Correlation analysis between single nucleotide polymorphism of FGF5 gene and wool yield in rabbits]. Yi Chuan. 2008, 30 (7): 893-899.
Murphy WJ, Pringle TH, Crider TA, Springer MS, Miller W: Using genomic data to unravel the root of the placental mammal phylogeny. Genome Res. 2007, 17 (4): 413-421. 10.1101/gr.5918807.
Roca AL, Bar-Gal GK, Eizirik E, Helgen KM, Maria R, Springer MS, O'Brien SJ, Murphy WJ: Mesozoic origin for West Indian insectivores. Nature. 2004, 429 (6992): 649-651. 10.1038/nature02597.
Maglio VJ: Origin and evolution of the Elephantidae. Trans Am Phil Soc. 1973, 63 (3): 1-149. 10.2307/1006229.
Haynes G: Mammoths, mastodonts, and elephants: biology, behavior, and the fossil record. 1991, Cambridge: Cambridge University Press
Iacumin P, Davanzo S, Nikolaev V: Short-term climatic changes recorded by mammoth hair in the Arctic environment. Palaeogeogr Palaeoclimatol. 2005, 218 (3-4): 317-324. 10.1016/j.palaeo.2004.12.021.
Winter H, Langbein L, Krawczak M, Cooper DN, Jave-Suarez LF, Rogers MA, Praetzel S, Heidt PJ, Schweizer J: Human type I hair keratin pseudogene phihHaA has functional orthologs in the chimpanzee and gorilla: evidence for recent inactivation of the human gene after the Pan-Homo divergence. Hum Genet. 2001, 108 (1): 37-42. 10.1007/s004390000439.
Miller W, Drautz DI, Ratan A, Pusey B, Qi J, Lesk AM, Tomsho LP, Packard MD, Zhao F, Sher A, et al: Sequencing the nuclear genome of the extinct woolly mammoth. Nature. 2008, 456 (7220): 387-390. 10.1038/nature07446.
Roca AL, Georgiadis N, Pecon-Slattery J, O'Brien SJ: Genetic evidence for two species of elephant in Africa. Science. 2001, 293 (5534): 1473-1477. 10.1126/science.1059936.
Greenwood AD: Late Pleistocene DNA extraction and analysis. Techniques in Molecular Systematics and Evolution. Edited by: DeSalle R, Giribet G, Wheeler W. 2002, Basel: Birkhauser-Verlag
Heinemeyer T, Wingender E, Reuter I, Hermjakob H, Kel AE, Kel OV, Ignatieva EV, Ananko EA, Podkolodnaya OA, Kolpakov FA, et al: Databases on transcriptional regulation: TRANSFAC, TRRD and COMPEL. Nucleic Acids Res. 1998, 26 (1): 362-367. 10.1093/nar/26.1.362.
Livingstone CD, Barton GJ: Protein sequence alignments: a strategy for the hierarchical analysis of residue conservation. Comput Appl Biosci. 1993, 9 (6): 745-756.
Olsen SK, Ibrahimi OA, Raucci A, Zhang F, Eliseenkova AV, Yayon A, Basilico C, Linhardt RJ, Schlessinger J, Mohammadi M: Insights into the molecular basis for fibroblast growth factor receptor autoinhibition and ligand-binding promiscuity. Proc Natl Acad Sci USA. 2004, 101 (4): 935-940. 10.1073/pnas.0307287101.
Eisenberg D: Three-dimensional structure of membrane and surface proteins. Annu Rev Biochem. 1984, 53: 595-623. 10.1146/annurev.bi.53.070184.003115.
Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982, 157 (1): 105-132. 10.1016/0022-2836(82)90515-0.
Luo Y, Lu W, Mohamedali KA, Jang JH, Jones RB, Gabriel JL, Kan M, McKeehan WL: The glycine box: a determinant of specificity for fibroblast growth factor. Biochemistry. 1998, 37 (47): 16506-16515. 10.1021/bi9816599.
Hecht HJ, Adar R, Hofmann B, Bogin O, Weich H, Yayon A: Structure of fibroblast growth factor 9 shows a symmetric dimer with unique receptor- and heparin-binding interfaces. Acta Crystallogr D Biol Crystallogr. 2001, 57 (Pt 3): 378-384. 10.1107/S0907444900020813.
Millar SE: Molecular mechanisms regulating hair follicle development. J Invest Dermatol. 2002, 118 (2): 216-225. 10.1046/j.0022-202x.2001.01670.x.
Roca AL, O'Brien SJ: Genomic inferences from Afrotheria and the evolution of elephants. Curr Opin Genet Dev. 2005, 15 (6): 652-659. 10.1016/j.gde.2005.09.014.
Wu DD, Irwin DM, Zhang YP: Molecular evolution of the keratin associated protein gene family in mammals, role in the evolution of mammalian hair. BMC Evol Biol. 2008, 8: 241-10.1186/1471-2148-8-241.
Langdon JH: Parsimony of aquatic and terrestrial hypotheses: how many hypotheses do we need?. Water and Human Evolution Symposium Proceedings: 30 April, 1999; Ghent, Belgium. 1999, [http://users.ugent.be/~mvaneech/langdon.htm]
Ryder ML: Hair of the Mammoth. Nature. 1974, 249 (5453): 190-191. 10.1038/249190a0.
Shoshani J, Tassy P: The Proboscidea: Evolution and Palaeoecology of Elephants and their Relatives. 1996, New York: Oxford University Press
Shoshani J, Walter RC, Abraha M, Berhe S, Tassy P, Sanders WJ, Marchant GH, Libsekal Y, Ghirmai T, Zinner D: A proboscidean from the late Oligocene of Eritrea, a "missing link" between early Elephantiformes and Elephantimorpha, and biogeographic implications. Proc Natl Acad Sci USA. 2006, 103 (46): 17296-17301. 10.1073/pnas.0603689103.
Liu AGSC, Seiffert ER, Simons EL: Stable isotope evidence for an amphibious phase in early proboscidean evolution. Proc Natl Acad Sci USA. 2008, 105 (15): 5786-5791. 10.1073/pnas.0800884105.
Gaeth AP, Short RV, Renfree MB: The developing renal, reproductive, and respiratory systems of the African elephant suggest an aquatic ancestry. Proc Natl Acad Sci USA. 1999, 96 (10): 5555-5558. 10.1073/pnas.96.10.5555.
Krause J, Dear PH, Pollack JL, Slatkin M, Spriggs H, Barnes I, Lister AM, Ebersberger I, Paabo S, Hofreiter M: Multiplex amplification of the mammoth mitochondrial genome and the evolution of Elephantidae. Nature. 2006, 439 (7077): 724-727. 10.1038/nature04432.
Rogaev EI, Moliaka YK, Malyarchuk BA, Kondrashov FA, Derenko MV, Chumakov I, Grigorenko AP: Complete mitochondrial genome and phylogeny of Pleistocene mammoth Mammuthus primigenius. PLoS Biol. 2006, 4 (3): e73-10.1371/journal.pbio.0040073.
Rohland N, Malaspinas AS, Pollack JL, Slatkin M, Matheus P, Hofreiter M: Proboscidean mitogenomics: chronology and mode of elephant evolution using mastodon as outgroup. PLoS Biol. 2007, 5 (8): e207-10.1371/journal.pbio.0050207.
Roca AL: The mastodon mitochondrial genome: a mammoth accomplishment. Trends Genet. 2008, 24 (2): 49-52. 10.1016/j.tig.2007.11.005.
Capelli C, MacPhee RD, Roca AL, Brisighelli F, Georgiadis N, O'Brien SJ, Greenwood AD: A nuclear DNA phylogeny of the woolly mammoth (Mammuthus primigenius). Mol Phylogenet Evol. 2006, 40 (2): 620-627. 10.1016/j.ympev.2006.03.015.
Ebersberger I, Galgoczy P, Taudien S, Taenzer S, Platzer M, von Haeseler A: Mapping human genetic ancestry. Mol Biol Evol. 2007, 24 (10): 2266-2276. 10.1093/molbev/msm156.
Kohn MH, Murphy WJ, Ostrander EA, Wayne RK: Genomics and conservation genetics. Trends Ecol Evol. 2006, 21 (11): 629-637. 10.1016/j.tree.2006.08.001.
Chimpanzee Sequencing and Analysis Consortium: Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005, 437 (7055): 69-87. 10.1038/nature04072.
Mol D, Coppens Y, Tikhonov AN, Agenbroad LD, MacPhee RDE, Flemming C, Greenwood A, Buigues B, de Marliave C, van Geel B, et al: The Jarkov mammoth: 20,000-year-old carcass of a Siberian woolly mammoth Mammuthus primigenius (Blumenbach, 1799). Proceedings of the 1st International Congress "The World of Elephants": 16-20 October 2001; Rome, Italy. 2001, 305-309.
Calvignac S, Terme JM, Hensley SM, Jalinot P, Greenwood AD, Hänni C: Ancient DNA identification of early 20th century simian T-cell leukemia virus type 1. Mol Biol Evol. 2008, 25 (6): 1093-1098. 10.1093/molbev/msn054.
Cooper A, Poinar HN: Ancient DNA: do it right or not at all. Science. 2000, 289 (5482): 1139-10.1126/science.289.5482.1139b.
Hanke M, Wink M: Direct DNA sequencing of PCR-amplified vector inserts following enzymatic degradation of primer and dNTPs. Biotechniques. 1994, 17 (5): 858-860.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9 (4): 286-298. 10.1093/bib/bbn013.
Smith TF, Waterman MS: Identification of common molecular subsequences. J Mol Biol. 1981, 147 (1): 195-197. 10.1016/0022-2836(81)90087-5.
Julenius K, Molgaard A, Gupta R, Brunak S: Prediction, conservation analysis, and structural characterization of mammalian mucin-type O-glycosylation sites. Glycobiology. 2005, 15 (2): 153-164. 10.1093/glycob/cwh151.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340 (4): 783-795. 10.1016/j.jmb.2004.05.028.
Ng PC, Henikoff S: SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003, 31 (13): 3812-3814. 10.1093/nar/gkg509.
Ramensky V, Bork P, Sunyaev S: Human non-synonymous SNPs: server and survey. Nucleic Acids Res. 2002, 30 (17): 3894-3900. 10.1093/nar/gkf493.
Cheng J, Randall AZ, Sweredoski MJ, Baldi P: SCRATCH: a protein structure and structural feature prediction server. Nucleic Acids Res. 2005, W72-76. 10.1093/nar/gki396. 33 Web Server
Cheng J, Randall A, Baldi P: Prediction of protein stability changes for single-site mutations using support vector machines. Proteins. 2006, 62 (4): 1125-1132. 10.1002/prot.20810.
Arnold K, Bordoli L, Kopp J, Schwede T: The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006, 22 (2): 195-201. 10.1093/bioinformatics/bti770.
Bennett-Lovsey RM, Herbert AD, Sternberg MJ, Kelley LA: Exploring the extremes of sequence/structure space with ensemble fold recognition in the program Phyre. Proteins. 2008, 70 (3): 611-625. 10.1002/prot.21688.
Orengo CA, Taylor WR: SSAP: sequential structure alignment program for protein structure comparison. Methods Enzymol. 1996, 266: 617-635. full_text.
Taylor WR, Orengo CA: Protein structure alignment. J Mol Biol. 1989, 208 (1): 1-22. 10.1016/0022-2836(89)90084-3.
Holm L, Park J: DaliLite workbench for protein structure comparison. Bioinformatics. 2000, 16 (6): 566-567. 10.1093/bioinformatics/16.6.566.
Lou M, Golding GB: FINGERPRINT: visual depiction of variation in multiple sequence alignments. Mol Ecol Notes. 2007, 7 (6): 908-914. 10.1111/j.1471-8286.2007.01904.x.
Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22 (21): 2688-2690. 10.1093/bioinformatics/btl446.
Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8 (3): 275-282.
Yang Z: Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods. J Mol Evol. 1994, 39 (3): 306-314. 10.1007/BF00160154.
Zhang Z, Schwartz S, Wagner L, Miller W: A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000, 7 (1-2): 203-214. 10.1089/10665270050081478.
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res. 2002, 12 (6): 996-1006.
Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12 (4): 656-664.
The authors wish to thank R.D.E. MacPhee and D. Mol (CERPOLEX/MAMMUTHUS) for providing material from the Jarkov mammoth for analysis. We are grateful to the following colleagues for assistance with living elephant samples: N. Georgiadis (Mpala Research Centre, Laikipia, Kenya), R. Hanson and S.J. O'Brien (National Cancer Institute, NIH, Frederick, MD, USA), B. York and A. Baker (Rosamond Gifford Zoo at Burnet Park, Syracuse, NY, USA), and the governments of Gabon and South Africa. Samples were collected in full compliance with specific federal permits. We thank F. Hussain and D. Doyle for technical assistance, and J. Kehler for helpful advice. AR and YI thank M. Gadd and R. Ruggiero of the U.S. Fish and Wildlife Service African Elephant Conservation Fund for support. GB was supported by the Russian Foundation for Fundamental Research No. 09-04-98568-r_vostok_a. NN was supported by start-up funds from CSUF. SOK was supported by a DARPA postdoctoral fellowship.
SF, KS, SH, MT, and ADG performed ancient DNA extractions, PCR and sequencing experiments. ALR and YI performed all modern elephant DNA work. NN, SOK and YI performed the bioinformatic and phylogenetic analyses and contributed to the writing of the manuscript. NN performed the protein comparison and structural modeling analysis. GB provided mammoth samples and morphological information. ALR and ADG designed the study, contributed to the experimental work and analysis and wrote the manuscript (with contributions from the others).
Electronic supplementary material
Additional file 1: Supplementary information for laboratory and analytical procedures. Ancient DNA laboratory procedures, elephant and mammoth FGF5 primer sequences, FGF5 sequence properties and phylogenetic analysis, and mammoth FGF5 clone sequences. (PDF 484 KB)
About this article
Cite this article
Roca, A.L., Ishida, Y., Nikolaidis, N. et al. Genetic variation at hair length candidate genes in elephants and the extinct woolly mammoth. BMC Evol Biol 9, 232 (2009). https://doi.org/10.1186/1471-2148-9-232