- Research article
- Open Access
Molecular evolution of the keratin associated protein gene family in mammals, role in the evolution of mammalian hair
BMC Evolutionary Biology volume 8, Article number: 241 (2008)
The Erratum to this article has been published in BMC Evolutionary Biology 2009 9:213
Hair is unique to mammals. Keratin associated proteins (KRTAPs), which contain two major groups: high/ultrahigh cysteine and high glycine-tyrosine, are one of the major components of hair and play essential roles in the formation of rigid and resistant hair shafts.
The KRTAP family was identified as being unique to mammals, and near-complete KRTAP gene repertoires for eight mammalian genomes were characterized in this study. An expanded KRTAP gene repertoire was found in rodents. Surprisingly, humans have a similar number of genes as other primates despite the relative hairlessness of humans. We identified several new subfamilies not previously reported in the high/ultrahigh cysteine KRTAP genes. Genes in many subfamilies of the high/ultrahigh cysteine KRTAP genes have evolved by concerted evolution with frequent gene conversion events, yielding a higher GC base content for these gene sequences. In contrast, the high glycine-tyrosine KRTAP genes have evolved more dynamically, with fewer gene conversion events and thus have a lower GC base content, possibly due to positive selection.
Most of the subfamilies emerged early in the evolution of mammals, thus we propose that the mammalian ancestor should have a diverse KRTAP gene repertoire. We propose that hair content characteristics have evolved and diverged rapidly among mammals because of rapid divergent evolution of KRTAPs between species. In contrast, subfamilies of KRTAP genes have been homogenized within each species due to concerted evolution.
The availability of the rapidly increasing number of genome sequences provides opportunities for investigators to study evolutionary patterns that potentially account for morphological characteristics, and suggest the genetic basis for variation in phenotypes. In particular, gene families in which duplications, rate variation and pseudogenization occur frequently are likely involved in functional innovation and adaptation . Examples of such gene families are those involved in the perception systems, for example, the odorant receptors [2–4], the vomeronasal receptors [5–8], and the sweet/umami and bitter receptors [9–12]. Furthermore, these studies should facilitate our understanding of the general evolutionary trends in genomic complexity and lineage-specific adaptation . Here, we studied the evolutionary patterns of the keratin associated protein (KRTAP) gene family, whose encoded proteins are major components of hair, with the goal of revealing the underlying basis for unique mammalian hair and its phenotypic diversity.
Hair is a unique character found on all mammals, but not on other animals, where it plays a crucial role in the retention of heat within these homoiotherms and presumably contributed significantly to the rapid radiation of mammals and their rise to become the dominant terrestrial vertebrate . Other functions of hair include sexual dimorphism, attraction of mates, and protection of skin . An interesting event in hair evolution has been its loss in humans [15, 16], however; humans actually have a similar density of hair follicles to that seen in apes . Comparative studies have concluded that hair presents similar structure and modality of growth throughout mammals [17–19]. For example, the overall ultrastructure of hair and the distribution of keratins in monotremes are similar to that of marsupial and placental mammals , and the localization of acidic and basic keratins in marsupial hairs is similar to that in placentals . However, most studies have focused on keratins rather than the keratin associated proteins.
The major components of hair are alpha-keratins and keratin associated proteins, each of which are encoded by multigene families. The alpha-keratins include two multigene subfamilies, type I (acidic) and type II (basic) [20, 21], and form the intermediate filament cytoskeleton of all epithelia providing stability against stress . In humans, the alpha-keratin gene family has been extensively studied demonstrating that there are 54 functional genes that are clustered on chromosomes 12q13.13 and 17q21.2 that show differing expression patterns during hair development . Hair keratins form an intermediate filament (IF) network by co-polymerization of type I and type II members, in trichocytes, which are cells that populate the central hair-forming compartment of the anagen hair follicle [24, 25]. In the hair cortex, hair keratins IFs are embedded in an in interfilamentous matrix, which consists of hair keratin-associated proteins (KATAP, usually abbreviated as KAP) . KRTAP contains two major groups: high/ultrahigh cysteine (HS) and high glycine-tyrosine (HGT) that are considered to have originated independently, and are essential for the formation of rigid and resistant hair shafts through their extensive disulfide bond cross-linking with the abundant cysteine residues of hair keratins or hydrophobic interactions with keratins [23, 26–28]. The genes have been grouped into 27 subfamilies, termed KRTAP1 to KRTAP27, based upon phylogeny [23, 29]. In humans, about 100 KRTAP genes are identified, that are in five tandemly arranged clusters (chromosomal regions 11p15.5, 11q13.4, 17q21.2, 21q22.1, and 21q22.3) [30–34]. Previous research on KRTAP genes has focused on their function and expression, with little emphasis on the origin and evolution of this gene family. Here, we investigate the evolution of the KRTAP gene family, including phylogeny and classification, and the mechanisms involved, such as gene duplication, gene conversion, in our endeavor to resolve hair's evolutionary history and to explain the diversity observed in extant mammals. We find different repertoires among mammals which potentially explain the differing hair features of different lineages. An expanded KRTAP gene repertoire was found in rodents. Surprisingly, despite the lack of hair, human had a similar number of genes with other primates. Large-scale gene conversion events were detected in high cysteine KRTAP but fewer in high glycine/tyrosine, and the latter genes evolve more dynamically. Compared to the conserved structure and modality of keratins within mammals, the significant divergence of characteristics of hair among placental, marsupial and monotreme species is likely due to interspecific divergence of KRTAP sequences.
Inventory of KRTAP genes in mammals
The recent rapid increase in the availability of comparative genomic data is facilitating the illumination of evolutionary features of organisms. Particularly, mammals are well represented with data from: placental mammals including primates – human , chimpanzee  and rhesus macaque ; rodents – mouse  and rat ; carnivore – dog ; a marsupial, the opossum ; and a monotreme, the duck-billed platypus .
Some gene sequences in GenBank http://www.ncbi.nlm.nih.gov have been submitted independently by several researchers. KRTAP genes have been shown to have size polymorphism within populations [43–45], and sequences from different individuals at these loci are diverse and may not align well. For instance, we identified that a KRTAP gene cluster in Contig: NT_113931 actually corresponds to the KRTAP region on Chr17 of the human genome, except that it differs in sequence length for several KRTAP genes. Therefore, we removed the genes from Contig: NT_113931 from our analyses. Furthermore, since many of the draft genome sequences were generated by whole genome shotgun assembly, we took care to identify redundant KRTAP sequences and excluded these from our analysis.
Since KRTAP genes are clustered as tandem arrays at just a few chromosomal locations and the genome sequences have high coverage (except platypus), we should be able to obtain near-complete and non-redundant KRTAP gene inventories from these mammalian genomes after in-depth screening and examination (summarized in additional file 1). In contrast to the mammalian genomes, we were unable to identify any KRTAP-like sequences from the chicken, lizard, Xenopus tropicalis or zebrafish genomes. We also note that the keratin associated protein 10-4 annotated in the chicken genome in GenBank (gene id: 425968) is actually an oncogenic transcription factor (JAC), and is not KRTAP-like nor has homology with any mammalian KRTAP sequences. As an alternative strategy to identify KRTAP genes in the chicken, we used Mapviewer from NCBI http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/mapview/ to identify genes flanking the KRTAP gene clusters in mammals and searched the chicken genome using BLAST to identify orthologous genomic regions. While conserved synteny of the flanking genes was observed within the chicken genome, KRTAP-like genes could not be found. These results indicate that the KRTAP gene clusters are unique to mammals and have been inserted into in the ancestral mammalian genome.
We classified our identified KRTAP gene sequences into two categories, intact and pseudo-genes. Pseudogenes are defined as coding sequences that are disrupted by frameshifts and/or stop codons. In placental mammals, rodents appear to contain an expanded KRTAP gene repertoire with a lower level of pseudogenization relative to other mammals both in high (including ultrahigh) cysteine and high glycine KRTAPs (fig. 1). Surprisingly, although humans have the lowest number of intact genes (101) and the highest number of pseudogenes (21, 17.2% of total number), the complete repertoire (122) is still similar in number to that of other mammals. Potentially, changes in the levels of expression of KRTAP genes may account for the relative hairlessness of humans. Opossum has a slightly larger repertoire of genes, and the repertoire in platypus is probably much larger than that we reported, since a more complete genome sequence is required to accurately define gene number.
Classification of KRTAP gene family, and tandem cluster in the genomes
Without the influence of evolutionary constraint, pseudogenes evolve faster then functional genes. Accordingly, substitution bias could generate errors in constructing phylogenies; therefore all pseudogenes were excluded from our preliminary analysis on subfamily classification. The high levels of divergence between species and/or subgroups and the homogenization within them caused by frequent gene conversion (see following text), may also introduce a bias into the construction of a phylogeny of all of the high cysteine KRTAPs. To circumvent this problem, we first constructed a neighbor-joining tree of just the human high cysteine KRTAP protein sequences, and used this phylogeny to identify subfamilies (Additional file 2: figure 1). The majority of genes for each subfamily were found as subfamily-specific tandem cluster in the genome (Additional file 2: figure 1 and Additional file 1). Genes with close relationships in our phylogeny tended to be most closely linked in the genome. This positive correlation between phylogeny and chromosomal location has previously been observed in several other gene families [46, 47]. Thus we used chromosomal location in combination with phylogeny to refine our subfamily classification. In addition we also considered amino acid sequence composition to generate our final classification. Next we constructed phylogenies of all of the high cysteine KRTAPs for each species alone (Additional file 2: figures 2–8), or combined with the human sequences (Additional file 2: figures 9–15), as well as a phylogeny of the combined mouse and rat sequences (Additional file 2: figure 16). From these phylogenies, together with chromosomal location and amino acid composition, we could identify several new subfamilies not reported previously, which we name subfamilies 28–35. Additionally, we grouped the previously characterized subfamilies 14 and 15 into subfamily 13 based on our phylogeny. In a similar manner, we grouped the glycine-tyrosine rich KRTAPs into 6 subfamilies (subfamilies 6, 7, 8, 19, 20 and 21) based on the phylogenetic tree (fig. 2, and Additional file 2: figure 17). The previously defined human subfamily 22 was combined with subfamily 19. Our recommendations have been endorsed by the HUGO Gene Nomenclature Committee .
Our phylogeny of intact genes indicated that genes within subfamilies are clustered in the genomic location (fig. 3, and Additional file 1). When pseudogenes were included in the phylogenetic analysis they were also found to be most closely related to their genomic neighbors. To classify interrupted genes that were annotated in unmapped (i.e., unassembled) genomic locations, BLAST was performed against the non-redundant database in GenBank to identify the best hit KRTAP genes, and phylogenetic classification of these were used to help classify the pseudogenes.
It should be noted that there is typically only about one gene per species for subfamilies 16, 24–27, and 29 (table 1), thus, these subfamilies were not used for our subsequent analysis for gene conversion.
Chromosome distribution of KRTAP
KRTAP genes are distributed mainly at five genomic regions in placental and marsupial genomes: Cluster 1 contains genes from subfamilies KRTAP 1, 2, 3, 4, 9, 17, 16, and 29. Cluster 2 contains genes from subfamilies 13, 24–27 and all glycine-tyrosine rich KRTAPs. Cluster 3 possesses genes from subfamilies 10 and 12. Cluster 4 encodes genes of subfamily 28. Cluster 5 corresponds to genes of subfamily 5 (fig. 3). Some variation in gene distribution is observed in some species. In rodents, the new subfamilies 30 and 31 have been inserted into the genomic locations of subfamilies 4 and 9 respectively. Paralogous KRTAP gene clusters have been mapped to human chromosomes 11q13 and chr11p15, with genes residing in these two clusters intermingled in the phylogeny (fig. 3.). Thus suggests that the KRTAP gene cluster at 11q13 is derived from 11p15 potentially representing a segmental duplication. The new chromosome 11q13 gene cluster is unique to primates as an orthologous region is also found in the chimpanzee and rhesus macaque genomes but not in others mammals (additional file 1). In a similar manner the dog genome has generated a new cluster on chromosome 31 that includes three genes, one from subfamily 10 and two from subfamily 12. Similar mechanisms for the origin of new genes at new genomic locations have been observed for other gene families [49, 50].
Amino acid composition comparison of KRTAP subfamilies
Previous research has classified the keratin associated proteins by their amino acid composition into three major groups: high-sulfur (~16–30% cysteine), ultra-high sulfur (> 30% cysteine), and high-glycine/tyrosine . Subfamilies 1, 2, 3, 10, 12, 16, 29 and 31, belong to the high-sulfur group; subfamilies 4, 5, 9, 17, 28, 30, 32 and 33 are ultra-high sulfur (table 2). Many high cysteine genes also have a high content of serine. Interestingly, subfamilies 11, 13, 24–27, 29, 34 and 35 have high serine content but relative low cysteine (table 2). The newly identified subfamilies 28 and 30, for which there is no functional or expression data, have the highest cysteine content (39.5%, 50.1% respectively).
Concerted evolution yields a high GC (Guanine and Cytosine) content in high cysteine KRTAP gene family
The evolutionary patterns observed in a multigene family can be attributed to two traditional models: concerted evolution and the birth-and-death process . Members of a gene family under concerted evolution evolve in a concerted manner rather than independently, as a mutation occurring in one member will spread through the entire gene family by the repeated occurrence of unequal crossover and/or gene conversion . Concerted evolution therefore results in the distance between pairs of genes remaining low. Considering the abundant species-specific clusters in the KRTAP gene phylogeny, many should be candidates for concerted evolution. Gene conversion plays a parallel role to unequal crossing over, with their major difference being that the latter can change the copy number of a gene; however, it is difficult to distinguish between these two mechanisms.
We identified potential gene conversion events within each subfamily using the methods implemented in the GeneConv program  which identifies identical fragments shared between pairs of nucleotide sequences. We found a large number of KRTAP gene pairs for which gene conversion events are suggested in the high/ultrahigh cysteine subfamilies, but significantly fewer, only 4 pairs, in high glycine-tyrosine subfamilies (χ2 = 73.85, p << 10-10) (summarized in additional file 3). The distribution of gene conversion events differs between species, suggesting that different levels of gene conversion occur in each species. For example, 32 gene pairs in rat subfamily 10 suggested evidence for gene conversion, but only 7 pairs were identified for this subfamily in the mouse. The RDP2 program also detected a large number of recombination events within the high cysteine KRTAP genes but fewer in the high glycine/tyrosine genes (additional file 4), and typically suggested that different gene pairs were involved. The identification of different gene pairs by the two methods may reflect differences in how the two programs identify recombination events.
Evidence for concerted evolution was also found in changes in base composition of the genes. Gene conversion is a nonreciprocal recombination process in which a DNA segment of a recipient gene is copied from a donor gene and occurs during the repair of double strand breaks by recombination . Recently, it has been discovered that gene conversion introduces a GC nucleotide bias into sequences, the biased gene conversion (BCG) concept, resulting in the enrichment of GC content in DNA sequences undergoing concerted evolution [52, 53] resulting in a positive correlation between rate of gene conversion and GC content. Therefore we examined the GC content of KRTAP coding sequences, which we propose scales with the rate of gene conversion. To evaluate the contribution of concerted evolution, we also scaled the level of gene conversion by distance between paralogs in each subfamily. It has been reported that the frequency of gene conversion correlates negatively with the divergence between gene pairs . Synonymous nucleotide sites are expected to evolve neutrally and can be used to evaluate the relative evolutionary divergence between pairs of genes. Intriguingly, high cysteine KRTAP genes contain extremely high GC content, and the GC content is negatively correlated with the divergence within the high cysteine KRTAP gene subfamilies (fig. 4). Comparatively, high glycine/tyrosine KRTAP contain relatively lower GC content (~50%), and a higher synonymous substitution rate (fig. 4), consistent with the detection of fewer gene conversion events.
Phylogenetic relationships can also suggest sequence homogenization, and potentially uncover the level of gene conversion. In the tree of human and chimpanzee KRTAP protein sequences, large quantities of one-to-one orthologies were identified (additional file 2: figure 9), suggesting that intraspecific gene conversion in the past 5 MYRs ago has not obscured orthologous relationships of the human and chimpanzee genes. The number of these one-to-one orthologies is reduced as phylogenetic distances between species pairs increases (additional file 2: figures 9–15).
More dynamical evolution of high glycine/tyrosine KRTAPs
Higher GC content in high cysteine KRTAP compared to high glycine/tyrosine genes reveals that strong concerted evolution is occurring in former, and that they evolve in a more stably manner. Accordingly, we compared the evolutionary dynamics of the two kinds of KRTAPs by calculating the Pearson correlation coefficient of the number of genes within each subfamily between species. As expected, the coefficient value is significantly higher for high cysteine KRTAP than for high glycine/tyrosine genes (p < 10-4 by Wilcoxon signed ranks test) (fig. 5A), and the two groups of coefficient values are positively correlated (fig. 5B). In particular, subfamily 20 evolved with a very dramatic variation in gene numbers, e.g. mouse and opossum have 30 and 42 genes respectively, but humans only have 2. As the number of gene conversion events detected was close to zero, the false positive rate for detecting positive selection in subfamily 20 should be very low. Positive selection, which was detected on the rat and opossum lineages by likelihood ratio tests (table 3), thus is proposed to be the major force for the dynamic evolution of subfamily 20 and possibly for other high glycine/tyrosine KRTAP genes.
We have described the near-complete inventories of KRTAP genes from the genomes of human, chimpanzee, and rhesus macaque representing primates, mouse and rat representing rodents, dog representing carnivore, opossum representing marsupial, and platypus as a monotreme, and conducted comprehensive analyses of the evolutionary patterns for this gene family, including phylogenetic classification, and the detection of gene conversion.
Among the mammals that we studied, the genomes of each species have members of most of the subfamilies suggesting that the majority of KRTAP gene subfamilies originated and diverged before the mammalian radiation (fig. 1). Accordingly, the mammalian ancestor should also have had a high diversity of KRTAP genes and enjoyed a similar range and spectrum of hair characteristics that is seen in modern mammalian species. The genes for subfamilies that emerged in the early mammals account for about 74% total number of genes in present mammalian genomes. We advocate that the rapid emergence of the KRTAP gene family correlates with the evolution of mammalian hair and that the rapid emergence of plentiful hair contributed to the successful radiation of homothermal mammals by helping them to retain body heat since hair acts as an insulator . In contrast, Homo sapiens has recently lost body hair function, presumably because humans can obtain heat and can keep cold out by using clothing. Despite a lack of hair, humans actually have a similar density of hair follicles to apes , which possibly explains why human do not have a significant fewer number of KRTAP genes. Perhaps, the changes in human hair are due to the reduction in expression of the KRTAP genes. In contrast to humans, rodents have an expanded KRTAP gene family. Perhaps, mouse and rat need more hair as they are adapted to the nocturnal environment. Despite the proposed similar range and diversity of hair characteristics in the ancestor of mammals, the current traits, and content, of hair likely diverged significantly from their ancestral characteristics as KRTAP sequences have diverged dramatically among species due to gene conversion.
Both gene conversion and unequal crossing over have occurred during the evolution of the KRTAP. The evolution of several KRTAP gene subfamilies, such as subfamilies 1, 2 and 3 fit well with a model of punctuated equilibrium (table 1), where morphological divergence as well as speciation occurs in a burst-like manner with rapid evolutionary change followed by long periods of constancy . These KRTAP genes show evidence for fewer gene duplication or loss events, suggesting that a low level of unequal crossing over has occurred within these subfamilies, and that gene conversion may predominate. In addition, the syntenic arrangement of genes of these subfamilies has not been disturbed by post-duplication rearrangement events, suggesting a strong constraint has been exerted upon their distribution. In contrast, other KRTAP genes, such as subfamilies 4 and 5 (table 1) have evolved relatively dynamically likely including unequal crossing over as changes in the copy number of this subfamily is observed between species.
Changes in gene numbers within subfamilies between mammals may explain differences in the observed hair morphology between species. No genes of subfamily 9 were identified in opossum or platypus, suggesting that this subfamily originated after the divergence of placental and marsupial mammals. In a similar manner, genes of subfamily 12 emerged before the divergence of marsupials but after the divergence from monotremes (fig. 1). Intriguingly, subfamily 13 has been lost on the marsupial lineage, and the subfamilies 24, 26 and 27 which are adjacent to subfamily 13 in the genome, are also not found in the opossum genome. We used Mapviewer from NCBI to identify genes flanking this cluster in the human genome and searched the opossum genome using BLAST to identify an orthologous genomic region. An ortholog of the human gene that flanks the human KRTAP genes is only about 13 kbp away from the remainder of the KRTAP cluster in the opossum genome, which indicates the loss of this region in the opossum was due to a deletion event. Subfamily 30, which is a member of the ultra-high sulfur (~50%) group, along with subfamily 34 are unique to rodents, thus may be partially responsible for the unique characteristics of hair in rodents. Subfamily 35 is mouse-specific, while subfamilies 32 and 33 are platypus-specific, which could account for species-specific hair characteristics.
The high-glycine/tyrosine (HGT) gene repertoire evolve more dynamically within mammals with increased levels of pseudogenization (fig. 1), for example, mice possesses 56 genes, while primates have considerably fewer genes with only 27 and 25 genes in human and chimpanzee respectively. The gene number in platypus is likely to be underestimated because HGT cluster in platypus is shorter than that of other mammals (only about 100 kbp in platypus, but ~500 kbp in human, ~800 kbp in mouse, and ~350 kbp in opossum), and the incomplete nature of this genome sequence. However, the dynamics of HGT subfamilies does not appear to be due to unequal crossing over, as unequal crossing over should also generate a GC content bias that is not observed in these genes (personal communication, Gabriel Marais).
Tandemly arrayed paralogous genes with similar function can provide combinatorial complexity to biological diversity . This extraordinary evolutionary feature has been observed in many multigene families responsible for processes that face enormous external signals. For instance, the sensory system such as the olfactory receptor, vomeronasal receptor and sweet/umami receptor, bitter receptor genes, require a huge combination of diverse receptors for the diverse ligands they encounter in a tremendous range of molecular environments . Another classic example is the immunoglobulin and T-cell receptor protein superfamily which use recombination to generate large quantities of antigen recognition complexes to allow an immune responses to rapidly evolving pathogens . We therefore speculate that an analogous process occurs within the KRTAP gene family, where the KRTAP and keratin proteins have combined in unique combinations to generate the high diversity of hair phenotypes that are observed both between and within species, and even within individuals.
We have described the near-complete inventories of KRTAP genes in eight mammalian genomes. We found that the KRTAP family was unique to mammals, KRTAP gene repertoire was expanded in rodents, and surprisingly, humans had a similar number of genes as other primates, inconsistent with the hairlessness of humans. The high glycine-tyrosine KRTAP genes have evolved more dynamically, with fewer gene conversion events and thus have a lower GC content compared with high cysteine KRTAPs. We propose that the mammalian ancestor should have a diverse KRTAP gene repertoire, and that hair content characteristics have evolved and diverged rapidly among mammals because of rapid divergent evolution of KRTAPs between species caused by concerted evolution.
Retrieval of sequences
Some KRTAP genes have previously been annotated to possess introns, which generally complicates gene prediction and identification. All of the intron-containing KRTAP genes have short introns, and the sequences of these introns are similar to the repeated regions found within the exons, and many are predicted to be alternatively spiced. All of the introns could be included in a primary transcript that can be translated in-frame with the coding exons, they just have longer repetitive regions. Therefore, we hypothesize that all of the KRTAP genes can generate an mRNA sequence that is intron-less.
We identified KRTAP gene repertoires in the genome assemblies from human (Homo sapiens) (build36.2), chimpanzee (Pan troglodytes) (build2.1), rhesus macaque (Macaca mulatta) (build1.1), mouse (Mus musculus) (build36.1), rat (Rattus norvegicus) (RGSC v3.4), dog (Canis familiaris) (build2.1), opossum (Monodelphis domesticus) (MonDom5) and platypus (Ornithorhynchus anatinus) (build1.1). We used the BLASTn algorithm  to search these genomes using all known human KRTAP genes as queries. Each newly identified putative KRTAP gene was used as a query using BLAST http://blast.ncbi.nlm.nih.gov/Blast.cgi against the non-redundant GenBank database to check whether their best hit was a KRTAP gene. The chicken (Gallus gallus) (Build 2.1), zebrafish (Danio rerio) genomes (Zv6) in NCBI , lizard (Anolis carolinensis) genome at UCSC http://genome.ucsc.edu/ and the Xenopus tropicalis genome at JGI  genomes were also searched using BLAST for KRTAP-like sequences. Sequences that possessed an interrupting stop codons and/or frame-shifts caused by insertions or deletions were denoted as pseudogenes.
In order to classify the members of the KRTAP gene family, protein sequences were used to construct phylogenetic trees using the neighbor-joining method with p-distances with MEGA3.0 http://www.megasoftware.net after alignment with ClustalW http://www.ebi.ac.uk/Tools/clustalw/index.html. The reliability of the trees was evaluated by the interior branch tests with 1,000 replications.
Detection of recombination
We employed the GeneConv program http://www.math.wustl.edu/~sawyer/geneconv/ to identify potential gene conversion events in the KRTAP coding sequences. Gene conversion is a process where a segment of DNA from one allele of a gene is copied and replaces the sequence in another allele or gene. Accordingly, GeneConv extends a method previously described by  and detects this process by identifying shared fragments between pairs of sequences. Global Bonferroni corrected P values were calculated to evaluate the statistical significance of the observed fragment lengths and are compared to a simulated distribution (10,000 iterations) of the same number of sequences with similar variation. Lower P values suggest a greater probability that a gene conversion event has occurred. GeneConv has a significant limitation in that it is unable to distinguish between gene conversion and unequal crossing over events, but to date no other effective bioinformatic method have been described that can distinguish between these two types of events. We therefore, can not distinguish gene conversion and unequal crossing over.
The divergence among paralogs could also reflect the level of conversion which can homogenize paralogs. Nucleotide sequences were back-translated from protein sequences after alignment by ClustalW . Alignments were modified manually if necessary. Approximate synonymous substitution rate (Ks) values within subfamilies were calculated by the modified Nei-Gojobori (p-distance) method with a transition/transversion ratio of 2 .
Recombination events were also detected using the RDP2 software package http://darwin.uvigo.es/rdp/rdp.html. Evidence for recombination was detected by running RDP, BootScan, MaxChi and Chimaera with 1,000 permutations. Sequences were considered linear. The highest acceptable P value cut-off was set to 0.01. Bonferroni correction was employed.
Adaptive evolution analysis
The site-specific models M1a and M2a implemented in PAML http://abacus.gene.ucl.ac.uk/software/paml.html were used to detect potentially positively selected sites in the subfamily [66–68]. Considering the high false positive rate of likelihood ratio tests, particularly when there is frequent recombination [69, 70], we only detected positive selection in subfamily 20, a family in which gene conversion does not appear to occur and has expanded dramatically on the mouse, rat and opossum lineages.
We calculated the potential number of gene pairs where gene conversion could occur by ∑Nij(Nij-1)/2, where Nij is the intact gene numbers in subfamily i within species j, and N must be higher than 2 for gene conversion to be detected by the GeneConv program. The chi test was used to detect statistical significance in the difference in number of gene conversion events occurring between high cysteine and high glycine/tyrosine KRTAP. The Pearson correlation coefficient of the number of genes in each subfamily was determined for pairs of species (table 1). The evolutionary dynamics of the high cysteine and high glycine/tyrosine KRTAP genes was evaluated by comparing the values of the Pearson correlation coefficients.
Hughes AL: Adaptive evolution of genes and genomes. 1999, Oxford University Press US
Niimura Y, Nei M: Evolutionary dynamics of olfactory receptor genes in fishes and tetrapods. Proc Natl Acad Sci USA. 2005, 102 (17): 6039-6044. 10.1073/pnas.0501922102.
Niimura Y, Nei M: Evolutionary dynamics of olfactory and other chemosensory receptor genes in vertebrates. J Hum Genet. 2006, 51 (6): 505-517. 10.1007/s10038-006-0391-8.
Niimura Y, Nei M: Extensive gains and losses of olfactory receptor genes in mammalian evolution. PLoS ONE. 2007, 2 (8): e708-10.1371/journal.pone.0000708.
Grus WE, Shi P, Zhang J: Largest vertebrate vomeronasal type 1 receptor (V1R) gene repertoire in the semi-aquatic platypus. Mol Biol Evol. 2007, 24: 2153-2157. 10.1093/molbev/msm157.
Grus WE, Shi P, Zhang Y, Zhang J: Dramatic variation of the vomeronasal pheromone receptor gene repertoire among five orders of placental and marsupial mammals. Proc Natl Acad Sci USA. 2005, 102 (16): 5767-5772. 10.1073/pnas.0501589102.
Shi P, Zhang J: Comparative genomic analysis identifies an evolutionary shift of vomeronasal receptor gene repertoires in the vertebrate transition from water to land. Genome Res. 2007, 17 (2): 166-10.1101/gr.6040007.
Young JM, Trask BJ: V2R gene families degenerated in primates, dog and cow, but expanded in opossum. Trends Genet. 2007, 23 (5): 212-215. 10.1016/j.tig.2007.03.004.
Fischer A, Gilad Y, Man O, Pääbo S: Evolution of bitter taste receptors in humans and apes. Mol Biol Evol. 2005, 22 (3): 432-436. 10.1093/molbev/msi027.
Parry CM, Erkner A, le Coutre J: Divergence of T2R chemosensory receptor families in humans, bonobos, and chimpanzees. Proc Natl Acad Sci USA. 2004, 101 (41): 14830-14834. 10.1073/pnas.0404894101.
Shi P, Zhang J, Yang H, Zhang Y: Adaptive diversification of bitter taste receptor genes in mammalian evolution. Mol Biol Evol. 2003, 20 (5): 805-814. 10.1093/molbev/msg083.
Shi P, Zhang J: Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes. Mol Biol Evol. 2006, 23 (2): 292-300. 10.1093/molbev/msj028.
Maderson PFA: Mammalian skin evolution: a reevaluation. Exp Derm. 2003, 12 (3): 233-236. 10.1034/j.1600-0625.2003.00069.x.
Bergman J: Why mammal body hair is an evolutionary enigma. Creation Research Society Quarterly. 2004, 40: 240-243.
Schwartz GG, Rosenblum LA: Allometry of primate hair density and the evolution of human hairlessness. Am J phys Anthrop. 1981, 55: 9-12. 10.1002/ajpa.1330550103.
Wheeler PE: The evolution of bipedality and loss of functional body hair in hominids. J Hum Evol. 1984, 13 (1): 91-98. 10.1016/S0047-2484(84)80079-2.
Alibardi L: Fine structure and immunocytochemistry of monotreme hairs, with emphasis on the inner root sheath and trichohyalin-based cornification during hair evolution. J Morphol. 2004, 261 (3): 345-363. 10.1002/jmor.10254.
Alibardi L: Fine structure of marsupial hairs, with emphasis on trichohyalin and the structure of the inner root sheath. J Morphol. 2004, 261 (3): 390-402. 10.1002/jmor.10257.
Alibardi L: Comparative aspects of the inner root sheath in adult and developing hairs of mammals in relation to the evolution of hairs. J Anat. 2004, 205 (3): 179-200. 10.1111/j.0021-8782.2004.00324.x.
Rogers MA, Winter H, Wolf C, Heck M, Schweizer J: Characterization of a 190-kilobase pair domain of human type I hair keratin genes. J Biol Chem. 1998, 273 (41): 26683-26691. 10.1074/jbc.273.41.26683.
Rogers MA, Winter H, Langbein L, Wolf C, Schweizer J: Characterization of a 300 kbp region of human DNA containing the type II hair keratin gene domain. J Invest Dermatol. 2000, 114 (3): 464-472. 10.1046/j.1523-1747.2000.00910.x.
Hesse M, Zimek A, Weber K, Magin TM: Comprehensive analysis of keratin gene clusters in humans and rodents. Eur J Cell Biol. 2004, 83 (1): 19-26. 10.1078/0171-9335-00354.
Rogers MA, Langbein L, Praetzel-Wunder S, Winter H, Schweizer J: Human hair keratin-associated proteins (KAPs). Int Rev Cytol. 2006, 251: 209-263. 10.1016/S0074-7696(06)51006-X.
Langbein L, Rogers MA, Winter H, Praetzel S, Beckhaus U, Rackwitz H-R, Schweizer J: The catalog of human hair keratins. I. Expression of the nine type I members in the hair follicle. J Biol Chem. 1999, 274 (28): 19874-19884. 10.1074/jbc.274.28.19874.
Langbein L, Rogers MA, Winter H, Praetzel S, Schweizer J: The catalog of human hair keratins. II. Expression of the six type II members in the hair follicle and the combined catalog of human type I and II keratins. J Biol Chem. 2001, 276 (37): 35123-35132. 10.1074/jbc.M103305200.
Powell BC, Nesci A, Rogers GE: Regulation of keratin gene expression in hair follicle differentiation. Ann NY Acad Sci. 1991, 642: 1-20.
Powell BC, Arthur J, Nesci A: Characterization of a gene encoding a cysteine-rich keratin associated protein synthesized late in rabbit hair follicle differentiation. Differentiation. 1995, 58 (3): 227-232. 10.1046/j.1432-0436.1995.5830227.x.
Powell BC, Rogers GE: The role of keratin proteins and their genes in the growth, structure and properties of hair. Formation and Structure of Human Hair. Edited by: Jolles PZH, Höcker H. 1997, Basel:Birkhäuser Verlag, 59-148.
Rogers MA, Winter H, Langbein L, Wollschlager A, Praetzel-Wunder S, Jave-Suarez LF, Schweizer J: Characterization of human KAP24. 1, a cuticular hair keratin-associated protein with unusual amino-acid composition and repeat structure. J Invest Dermatol. 2007, 127: 1197-1204. 10.1038/sj.jid.5700702.
Rogers MA, Langbein L, Winter H, Ehmann C, Praetzel S, Korn B, Schweizer J: Characterization of a cluster of human high/ultrahigh sulfur keratin-associated protein genes embedded in the type I keratin gene domain on chromosome 17q12-21. J Biol Chem. 2001, 276 (22): 19440-19451. 10.1074/jbc.M100657200.
Rogers MA, Langbein L, Winter H, Ehmann C, Praetzel S, Schweizer J: Characterization of a first domain of human high glycine-tyrosine and high sulfur keratin-associated protein (KAP) genes on chromosome 21q22. J Biol Chem. 2002, 277 (50): 48993-49002. 10.1074/jbc.M206422200.
Rogers MA, Langbein L, Winter H, Beckmann I, Praetzel S, Schweizer J: Hair keratin associated proteins: characterization of a second high sulfur KAP gene domain on human chromosome 21. J Invest Dermatol. 2004, 122: 147-158. 10.1046/j.0022-202X.2003.22128.x.
Shibuya K, Obayashi I, Asakawa S, Minoshima S, Kudoh J, Shimizu N: A cluster of 21 keratin-associated protein genes within introns of another gene on human chromosome 21q22. 3. Genomics. 2004, 83 (4): 679-693. 10.1016/j.ygeno.2003.09.024.
Yahagi S, Shibuya K, Obayashi I, Masaki H, Kurata Y, Kudoh J, Shimizu N: Identification of two novel clusters of ultrahigh-sulfur keratin-associated protein genes on human chromosome 11. Biochem Biophys Res Commun. 2004, 318: 655-664. 10.1016/j.bbrc.2004.04.074.
International Human Genome Sequencing Consortium: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.
The Chimpanzee Sequencing and Analysis Consortium: Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005, 437: 69-87. 10.1038/nature04072.
Rhesus Macaque Genome Sequencing and Analysis: Evolutionary and biomedical insights from the rhesus macaque genome. Science. 2007, 316: 222-234. 10.1126/science.1139247.
Mouse Genome Sequencing Consortium: Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420: 520-562. 10.1038/nature01262.
Rat Genome Sequencing Project Consortium: Genome sequence of the Brown Norway rat yields insights into mammalian evolution. Nature. 2004, 428: 493-521. 10.1038/nature02426.
Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, Clamp M, Chang JL, Kulbokas EJ, Zody MC, Mauceli E, Xie X, Breen M, Wayne RK, Ostrander EA, Ponting CP, Galibert F, Smith DR, deJong PJ, Kirkness E, Alvarez P, Biagi T, Brockman W, Butler J, Chin C-W, Cook A, Cuff J, Daly MJ, DeCaprio D, Gnerre S, Grabherr M, Kellis M, Kleber M, Bardeleben C, Goodstadt L, Heger A, Hitte C, Kim L, Koepfli K-P, Parker HG, Pollinger JP, Searle SMJ, Sutter NB, Thomas R, Webber C, Lander ES: Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005, 438: 803-819. 10.1038/nature04338.
Mikkelsen TS, Wakefield MJ, Aken B, Amemiya CT, Chang JL, Duke S, Garber M, Gentles AJ, Goodstadt L, Heger A, Jurka J, Kamal M, Mauceli E, Searle SMJ, Sharpe T, Baker ML, Batzer MA, Benos PV, Belov K, Clamp M, Cook A, Cuff J, Das R, Davidow L, Deakin JE, Fazzari MJ, Glass JL, Grabherr M, Greally JM, Gu W, Hore TA, Huttley GA, Kleber M, Jirtle RL, Koina E, Lee JT, Mahony S, Marra MA, Miller RD, Nicholls RD, Oda M, Papenfuss AT, Parra ZE, Pollock DD, Ray DA, Schein JE, Speed TP, Thompson K, VandeBerg JL, Wade CM, Walker JA, Waters PD, Webber C, Weidman JR, Xie X, Zody MC, Graves JAM, Ponting CP, Breen M, Samollow PB, Lander ES, Lindblad-Toh K: Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences. Nature. 2007, 447: 167-177. 10.1038/nature05805.
Warren WC, Hillier LW, Graves JAM, Birney E, Ponting CP, Grützner F, Belov K, Miller W, Clarke L, Chinwalla AT, Yang S-P, Heger A, Locke DP, Miethke P, Waters PD, Veyrunes F, Fulton L, Fulton B, Graves T, Wallis J, Puente XS, López-Otín C, Ordóñez GR, Eichler EE, Chen L, Cheng Z, Deakin JE, Alsop A, Thompson K, Kirby P, Papenfuss AT, Wakefield MJ, Olender T, Lancet D, Huttley GA, Smit AFA, Pask A, Temple-Smith P, Batzer MA, Walker JA, Konkel MK, Harris RS, Whittington CM, Wong ESW, Gemmell NJ, Buschiazzo E, Jentzsch IMV, Merkel A, Schmitz J, Zemann A, Churakov G, Kriegs JO, Brosius J, Murchison EP, Sachidanandam R, Smith C, Hannon GJ, Tsend-Ayush E, McMillan D, Attenborough R, Rens W, Ferguson-Smith M, Lefèvre CM, Sharp JA, Nicholas KR, Ray DA, Kube M, Reinhardt R, Pringle TH, Taylor J, Jones RC, Nixon B, Dacheux J-L, Niwa H, Sekita Y, Huang X, Stark A, Kheradpour P, Kellis M, Flicek P, Chen Y, Webber C, Hardison R, Nelson J, Hallsworth-Pepin K, Delehaunty K, Markovic C, Minx P, Feng Y, Kremitzki C, Mitreva M, Glasscock J, Wylie T, Wohldmann P, Thiru P, Nhan MN, Pohl CS, Smith SM, Hou S, Renfree MB, Mardis ER, Wilson RK: Genome analysis of the platypus reveals unique signatures of evolution. Nature. 2008, 453 (7192): 175-183. 10.1038/nature06936.
Kariya N, Shimomura Y, Ito M: Size polymorphisms in the human ultrahigh sulfur hair keratin-associated protein 4, KAP4, gene family. J Invest Dermatol. 2005, 124: 1111-1118. 10.1111/j.0022-202X.2005.23662.x.
Rogers MA, Schweizer J: Human KAP genes, only the half of it? extensive size polymorphisms in hair keratin-associated protein genes. J Invest Dermatol. 2005, 124 (6): vii-ix. 10.1111/j.0022-202X.2005.23728.x.
Shimomura Y, Aoki N, Schweizer J, Langbein L, Rogers MA, Winter H, Ito M: Polymorphisms in the human high sulfur hair keratin-associated protein 1, KAP1, gene family. J Biol Chem. 2002, 277 (47): 45493-45501. 10.1074/jbc.M206398200.
Huntley S, Baggott DM, Hamilton AT, Tran-Gyamfi M, Yang S, Kim J, Gordon L, Branscomb E, Stubbs L: A comprehensive catalog of human KRAB-associated zinc finger genes: Insights into the evolutionary history of a large family of transcriptional repressors. Genome Res. 2006, 16 (5): 669-677. 10.1101/gr.4842106.
Nei M, Rooney AP: Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005, 39: 121-152. 10.1146/annurev.genet.39.073003.112240.
Hamilton AT, Huntley S, Tran-Gyamfi M, Baggott DM, Gordon L, Stubbs L: Evolutionary expansion and divergence in the ZNF91 subfamily of primate-specific zinc finger genes. Genome Res. 2006, 16 (5): 584-594. 10.1101/gr.4843906.
Niimura Y, Nei M: Evolution of olfactory receptor genes in the human genome. Proc Natl Acad Sci USA. 2003, 100 (21): 12235-12240. 10.1073/pnas.1635157100.
Sawyer S: Statistical tests for detecting gene conversion. Mol Biol Evol. 1989, 6 (5): 526-538.
Marais G: Biased gene conversion: implications for genome and sex evolution. Trends Genet. 2003, 19 (6): 330-338. 10.1016/S0168-9525(03)00116-1.
Galtier N: Gene conversion drives GC content evolution in mammalian histones. Trends Genet. 2003, 19 (2): 65-68. 10.1016/S0168-9525(02)00002-1.
Drouin G: Characterization of the gene conversions between the multigene family members of the yeast genome. J Mol Evol. 2002, 55 (1): 14-23. 10.1007/s00239-001-0085-y.
Gould SJ, Eldredge N: Punctuated equilibrium comes of age. Nature. 1993, 366 (6452): 223-227. 10.1038/366223a0.
Noonan JP, Grimwood J, Schmutz J, Dickson M, Myers RM: Gene conversion and the evolution of protocadherin gene cluster diversity. Genome Res. 2004, 14 (3): 354-366. 10.1101/gr.2133704.
Mombaerts P: Genes and ligands for odorant, vomeronasal and taste receptors. Nat Rev Neurosci. 2004, 5 (4): 263-278. 10.1038/nrn1365.
Flajnik MF: Comparative analyses of immunoglobulin genes: surprises and portents. Nat Rev Immunol. 2002, 2 (9): 688-698. 10.1038/nri889.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.
Kumar S, Tamura K, Nei M: MEGA3. 0: Integrated software tot molecular evolutionary genetics analysis and sequence alignment. Brief Bioinlbrm. 2000, 4 (5): 150-163.
Thompson JD, Higgins DG, Gibson TJ: CLUSTALW: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Rosenberg MS, Subramanian S, Kumar S: Patterns of transitional mutation biases within and among mammalian genomes. Mol Biol Evol. 2003, 20 (6): 988-993. 10.1093/molbev/msg113.
Martin DP, Williamson C, Posada D: RDP2: recombination detection and analysis from sequence alignments. Bioinformatics. 2005, 21 (2): 260-262. 10.1093/bioinformatics/bth490.
Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biol Sci. 1997, 13: 555-556.
Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15 (5): 568-573.
Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics. 1998, 148 (3): 929-936.
Anisimova M, Nielsen R, Yang Z: Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites. Genetics. 2003, 164 (3): 1229-1236.
Zhang J, Nielsen R, Yang Z: Evaluation of an improved branch-Site likelihood method for detecting positive selection at the molecular level. Mol Biol Evol. 2005, 22 (12): 2472-2479. 10.1093/molbev/msi237.
We are grateful to the invaluable suggestions and help with the classification of genes from Dr. Michael A. Rogers at the Section of Normal and Neoplastic Epidermal Differentiation, German Cancer Research Center, Heidelberg, Germany, and the kind help of Dr Matt Wright from the HUGO Gene Nomenclature Committee. We thank Dr. Peng Shi, Ms. Jing He, and Mr. Guo-Dong Wang for comments. This work was supported by grants from the National Basic Research Program of China (973 Program, 2007CB411600), the National Natural Science Foundation of China (30621092, 30430110), and Bureau of Science and Technology of Yunnan Province.
DDW, DMI and YPZ designed the research and outlined the manuscript together, and DDW drafted the manuscript. ALL authors have read and approved the final manuscript.
An erratum to this article is available at http://0-dx.doi.org.brum.beds.ac.uk/10.1186/1471-2148-9-213.
Electronic supplementary material
Additional file 1: table 1. KRTAP genes in the human, chimpanzee, rhesus macaque, dog, mouse, rat, opossum, and platypus genomes. (XLS 191 KB)
Additional file 2: figure 1–figure 17. Figure 1–Figure 16 are the phylogenetic trees of high/ultrahigh cysteine KAPs of human (Figure 1), chimpanzee (Figure 2), rhesus macaque(figure 3), dog (figure 4), mouse (figure 5), rat (figure 6), opossum (figure 7), platypus(figure 8), human and chimpanzee (figure 9), human and rhesus macaque (figure 10), human and dog (figure 11), human and mouse (figure 12), human and rat (figure 13), human and opossum (figure 14), human and platypus (figure 15), mouse and rat (figure 16). Figure 17 is the phylogenetic tree of high glycine/tyrosine KRTAPs. h represents human, c is chimpanzee, rh is rhesus macaque, d-dog, m-mouse, r is rat, o is opossum, and p is platypus. The values on the branches are reliabilities, which are evaluated by the interior branch tests with 1,000 replications. Only values higher than 50% are noted. (PDF 201 KB)
Additional file 3: table 2. Pairs of genes with significant statistical support for gene conversion. The significances are calculated by Bonferroni-corrected method. (DOC 203 KB)
Additional file 4: table 3. Results of recombination detection by RDP2 program with algorithms: RDP, BootScan, MaxChi and Chimaera with 1,000 permutations. Sequences were considered linear. The highest acceptable P value cut-off was set to 0.01, and the Bonferroni correction was employed. The numbers are the unique events (recombination signals). (DOC 52 KB)
Authors’ original submitted files for images
About this article
Cite this article
Wu, D., Irwin, D.M. & Zhang, Y. Molecular evolution of the keratin associated protein gene family in mammals, role in the evolution of mammalian hair. BMC Evol Biol 8, 241 (2008) doi:10.1186/1471-2148-8-241
- Gene Conversion
- Concerted Evolution
- Chicken Genome
- Gene Conversion Event
- Hair Keratin