Comparative genomics analyses of alpha-keratins reveal insights into evolutionary adaptation of marine mammals
Frontiers in Zoologyvolume 14, Article number: 41 (2017)
Diversity of hair in marine mammals was suggested as an evolutionary innovation to adapt aquatic environment, yet its genetic basis remained poorly explored. We scanned α-keratin genes, one major structural components of hair, in 16 genomes of mammalian species, including seven cetaceans, two pinnipeds, polar bear, manatee and five terrestrial species.
Extensive gene loss and high pseudogenization rate of α-keratin genes were identified in cetaceans when compared to terrestrial artiodactylans (average number of α-keratins 37.29 vs. 58.33; pseudogenization rate 29.89% vs. 8.00%), especially of hair follicle-specific keratin genes (average pseudogenization rate in cetaceans of 43.88% relative to 3.80% artiodactylian average). Compared to toothed whale, the much more number of intact functional α-keratin genes was examined in the baleen whale that had specific keratinized baleen. In contrast, the number of keratin genes in pinnipeds, polar bear and manatee were comparable to those of their respective terrestrial relatives. Additionally, four keratin genes (K39, K9, K42, and K74) were found to be pseudogenes or lost uniquely in cetaceans and manatees.
Species-specific evolution of α-keratin gene family identified in the marine mammals might be responsible for their different hair characteristics. Increased gene loss and pseudogenization rate identified in cetacean lineages was likely to contribute to hair-less phenotype to adaptation for complete aquatic environment. However, the fully aquatic manatee still remained the comparable number of intact genes to its terrestrial relative, probably due to its perioral bristles and bristle-like hairs on the oral disk. By contrast, similar evolution pattern of α-keratin gene repertoire in the pinnipeds, polar bear and their terrestrial relatives was likely due to abundant hair to keep warm when they went ashore. Interestingly, some keratin genes were exclusively lost in cetaceans and manatees, likely as a result of convergent hair-loss phenotype to inhabit completely aquatic environment in both groups.
Marine mammals are specific groups with evolutionary histories that their separate terrestrial relatives returned to the ocean on separate occasions and adapted to living all or part of their life in the aquatic environment. The living groups of marine mammals are generally including the following five groups: pinnipeds (seals, sea lions, fur seals, and walruses), cetaceans (whales, dolphins, and porpoises), sea otters, sirenians (dugongs and manatees), and polar bears . Despite their independent evolutionary origins, these clades of the marine mammals have developed a series of adaptations to full or part aquatic environments, including morphological (e.g. streamlined shape, paddle-like limbs and feet) and physiological features adaptations (such as superb diving skill, echolocation, the thickened blubber, and etc.) .
Hair emergence is one of the major innovations in the mammalian evolution. Hair plays a key role in protection from mechanical insults, facilitated homeothermy, sense the immediate surrounding, sexual dimorphism, attraction of mates and etc. [3, 4]. Similar to other mammals, marine mammals are characterized by the presence of hair coat but with different evolution pattern among them. In order to reduce drag in the water, cetaceans and sirenians lack of hair coat because both are completely aquatic though cetaceans have body hair temporarily at fetuses [5, 6]. In contrast, the walrus has lost much of their hair (fur) and are characterized by thick layers of blubber under the skin to keep warm due to spending considerable time on land . The Weddell seal grows a thin fur coat around its whole body except for small areas around the flippers . The polar bear’s fur consists of a dense layer of ‘under fur’ and an outer layer of guard hairs that appear white to tan but are actually transparent. The guard hair is 5–15 cm long over most of the bear’s body . The sea otter, however, has an exceptionally thick and dense coat of fur although almost time living in the sea . However, the genetic bases of hair evolution in these clades of marine mammals remain poorly explored.
Hair is a kind of strongly keratinized tissue and is mainly composed of alpha-keratins and keratin associated-proteins (KRTAPs) [9, 10], which are encoded by a large number ssof multigene families and arranged in clusters on chromosomes [9, 10]. Alpha-keratins and KRTAPs are fibre-reinforced structures consisting of intermediate filaments embedded in an amorphous protein matrix which can provide support for the stability and rigidity of epidermal cells and tissue morphology . The KRTAPs, as well known KAPs, are unique to mammals. According to the amino acid composition, KRTAP was divided into two major groups: high/ultrahigh cysteine (HS) and high glycine-tyrosine (HGT) that are essential for the formation of rigid and resistant hair shafts with α-keratins [9, 12]. By contrast, the α-keratin gene family was divided into two categories, i.e. epithelial keratin genes and hair keratin genes. A total of 54 functional keratin genes were identified in the human genome, which can be divided into 28 type I and 26 type II genes [10, 11]. The type I keratins consist of 17 epithelial and 11 hair keratins while the type II members comprise 20 epithelial and 6 hair keratins. In the 37 epithelial keratins, nine are specifically expressed in the hair follicle. Thus, the nine genes and 17 hair keratins are collectively referred to as hair follicle-specific keratins, which are reported to control the growth and variation of hair . Previous studies characterized KRTAP gene repertoire in 22 mammals and found that gene family repertoire expansion, contraction, and pseudogenization were related to hair diversities in mammals . In addition, an increase rate of hair keratin genes loss and pseudogenization were examined in two cetacean species when compared with fur terrestrial relatives, which was suggested to be associated with cetacean hairless phenotype . However, previous studies only determined in one or two cetacean species could not provide comprehensive insight into genetic basis of hair diversities in marine mammals. Thus, more marine mammals were added in the present study, including seven cetacean species, two pinniped species, polar bear as well as one manatee. We first scanned the α-keratin repertoire from the genomes of the 11 marine mammals, and compared with that of their terrestrial relatives to test if α-keratin gene diversities are related to their hair features. Moreover, we tested whether there was common genes loss along with pseudogenization in the full aquatic cetaceans and sirenians as response for convergent hair-loss phenotype. Specially, we also tested if baleen whales retained much more α-keratin genes than toothed whales due to its specific keratinized baleen. From all the analyses, we expected to address the underlying genetic basis of the different hair phenotype in marine mammals.
α-keratin and KRTAP genes scanning and identification
We first identified α-keratin gene repertoires in the cow (Bos taurus) (coverage 7×,Btau_4.6.1) from UCSC Genome Browser website (http://genome.ucsc.edu/) taking all the known α-keratin gene sequences of human  as queries using BLASTN and TBLASTN algorithm . All exon/intron junctions follow the cannonical AG/GT rule of splicing . Meanwhile, we combined the online websites GENEWISE prediction (http://www.ebi.ac.uk/Tools/psa/genewise/) to verifying the accuracy of the sequence. The putative α-keratin genes in the cow were then taken as queries to explore the α-keratin multigene family in the genome of 11 marine mammals, including bowhead whale (Balaena mysticetus, coverage 150×, http://www.bowhead-whale.org/downloads/), minke whale (B. acutorostrata, coverage 92×, BalAcu1.0, NCBI), sperm whale (Physeter macrocephalus, coverage 75×, Physeter_macrocephalus_2.0.2, NCBI), Yangtze River dolphin (Lipotes vexillifer, Lipotes_vexillifer_v1, coverage 115×, NCBI), killer whale (Orcinus orca, coverage 200×, Oorc_1.1, NCBI), bottlenose dolphin (Tursiops truncatus, coverage 30×, Ttru_1.4, UCSC), Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis, unpublished data), Weddell seal (Leptonychotes weddellii, coverage 82×, LepWed_1.0, NCBI), Florida manatee (Trichechus manatus latirostris, coverage 150×, TriManLat_1.0, NCBI), Polar bear (Ursus maritimus, coverage 101×, UrsMar_1.0, NCBI), Pacific walrus (Odobenus rosmarus, coverage 200×, Oros_1.0, NCBI, Additional file 1). The orthologous α-keratin genes were also scanned from the genome of their respective terrestrial relatives, such as sheep (Ovis aries, coverage 142×, Oar_v3.1, NCBI), alpaca (Vicugna pacos, coverage 22×, Vicugna_pacos-2.0.1, NCBI), African savanna elephant (Loxodonta africana, coverage 7×, Loxafr3.0, NCBI), giant panda (Ailuropoda melanoleuca, coverage 60×, AilMel_1.0, NCBI) (Additional file 1). Type I and type II keratins have their unique flanking sequences, which can ensures maximum searching these genes in mammalian genomes. In addition, in order to obtain as much as possible of the keratin repertoire, all the newly annotated gene sequences were taken as queries in blast searches against their own genomes. Finally, all genes were checked to test whether their best hit was a α-keratin gene from NCBI genome database http://www.ncbi.nlm.nih.gov/blast/Blast.cgi using BLASTN algorithm . All the identified α-keratin genes were separated into three categories: intact genes, incomplete genes (abbreviated to ic and pseudogenes (abbreviated to p), according to amino acid alignment and blast results. All the keratins were named according to the revised keratin nomenclature of human .
In addition, the published human  KRTAP genes were as queries to search unannotated genomic sequences of KRTAP genes in seal, walrus and manatee by using BLASTN algorithm . All KRTAP genes have single exon in mammals. The homology fragment more than 30 bp was retrieved in the present study in order to scan the KRTAP gene repertoire as fully as possible. Like α-keratin genes, all the newly annotated KRTAP gene sequences were used as queries in blast searches against the own genomes. In the end all KRTAP genes were checked whether their best hit was a KRTAP gene from NCBI genome database.
Phylogenetic inference for α-keratin genes classification
First, all nucleotide sequences were aligned using MUSCLE in MEGA6.05  and checked by eye. Second, to classify type I and II keratin genes to their respective family, we reconstructed the phylogenetic tree under the GTR + G+ I model using Maximum likelihood (ML) by RAxML  and Bayesian inference (BI) by MrBayes 3.2.3 . The ML tree evaluated the best tree for each cluster, and supported for the nodes obtained with 1000 replications. For the Bayesian analyses, two simultaneous independent runs were performed for 40 × 106 iterations of a Markov ChainMonte Carlo algorithm, with six simultaneous chains, sampling every 1000 generations. In case of misplacement in phylogenetic trees, we will recheck their sequence for correction. Based on the sequence homology and phylogenetic relationship, we divided α-keratin genes into different subfamilies. Finally, we mapped α-keratin gene organization of each species on their genomes.
Gene conversion and recombination detection
To detect whether the gene conversion and recombination presented among keratin gene family that grouped together in the same species rather than with the members of the same family, we used the RDP4 software  to detect gene conversion and recombination events using RDP, Geneconv, Bootscan, MaxChi and Chimaera with 1000 permutations and cut off p value of 0.001.
Likelihood analysis of gene gain and loss
To estimate the average gene gain/loss rate and to identify gene families that have undergone significant size changes, we applied CAFÉ v3.0 , a tool for the statistical analysis of the evolution of the gene family size. CAFÉ v3.0 can be estimated gene gain and loss rate among α-keratin gene family, calculated ancestral states of gene family sizes for each node in the phylogenetic tree and identified the average expansion or contraction gene family on each branch. For the CAFÉ analysis, we reconstructed the ultrametric tree as a starting tree for Bayesian inference of tree topologies and node ages using Markov chain Monte Carlo (MCMC) in BEAST v1.8.1 . We ran six independent MCMC, 3 × 108 steps long under individual gene models previously selected for reconstruction. We then checked for the convergence and stationarity using Tracer v1.6 . Finally, we extracted the maximum clade credibility tree for combined tree sets using TreeAnnotator v1.8.1 .
Results and discussion
α-keratin gene family related to hair features
The advent of mainstream sequencing technologies has allowed the whole genomes of many mammals to be assembled, which provides possibility to explore the genetic basis of adaptation to their specific lifestyles at genome level. In the present study, we scanned the α-keratin repertoire from the genomes of 11 marine mammals (including seven cetaceans, two pinnipeds, polar bear and manatee) and compared with their terrestrial relatives (Table 1) to explore the genetic basis of hair diversity in marine mammals as response for their adaptations to aquatic lifestyle with different degrees. The genomic regions containing keratin gene families are reported to be conservative with special flanking sequences of two clusters, which facilitate to obtain the complete keratin repertoire in mammalian genomes. For example, type I keratins are flanked by SMARCE1 and EIF1 genes while type II keratins are flanked by FAIM2 and EIF4B genes (Fig. 1). However, we have not retrieved the FAIM2 gene in the 5′ flanking regions of five species, including Weddell seal, giant panda, bottlenosed dolphin, sperm whale, and bowhead whale (Fig. 1). Actually, the FAIM2 gene was identified in such five genomes but located in the different chromosomes with α-keratins. Similar phenomenon was also found in the EIF1and EIF4B that was not found in the flanking region of Weddell seal and bottlenose dolphin, respectively (Fig. 1). A total of 780 α-keratin genes were identified in the 16 mammalian species examined in our study, of which 383 belong to type I keratins, 397 belong to type II keratins.
In contrast with the terrestrial mammals, the marine mammals have lost some α-keratins and are characterized with increased pseudogenization rate, likely related with their hair diversities in response to their living environment. Average of 44.27 α-keratin genes observed in the marine mammals has slightly differences to the terrestrial mammals which has averagely 58.8 α-keratins. Moreover, the pseudogenization rate of marine mammals (18.69%) is nearly three times of that of their terrestrial relatives (7.17%). In the case of the marine mammals, however, different evolution pattern of α-keratins among them might be responsible for their different hair characteristics. We found an average number of α-keratin genes in cetaceans of 37.29, which was significantly lower than that of artiodactylian species with an average number of 58.33. The number of α-keratin genes ranged from 32 in the sperm whale to 45 in the minke whale, while in artiodactylans it ranged from 54 in the alpaca to 61 in the cow (Table 1). The higher rate of pseudogenization rate was identified in cetacean species (average 29.89% compared to artiodactylans’ average of 8.00%). The cetacean species are nearly completely absent of hair coat with only a few bristle in order to reduce drag of swimming for adaptation to completely aquatic environment. There is evidence of an apparent reduction in the α-keratin gene repertoire and increased rate of pseudogenization in cetaceans, which may be associated with their hairless phenotype. For the fully aquatic manatee, body hairs are greatly reduced but with perioral bristles and bristle-like hairs on the oral disk that have specialized sensory and feeding function . Accordingly, the manatee still remained the comparable number (56 in manatee vs. 60 in elephant) and rate of pseudogenization (8.33% in manatee vs. 10.71% in elephant) to its terrestrial relative, i.e. elephant (Table 1). By contrast, gene repertoires in the pinnipeds (56) and polar bear (58) were similar to that of panda (58) (Table 1). In addition, we detected 49 and 45 intact genes in pinnipeds and polar bear, comparable to that of the panda with 50 ones. Similar evolution pattern of α-keratin gene repertoire in the pinnipeds, polar bear and their terrestrial relatives is likely due to abundant hair to keep warm when they went ashore to mate, give birth, molt or escape from predators.
Further, we found that the gene gain and loss of α-keratin was consistent with the hair diversification of marine mammals. we employed the software CAFÉ to estimate the gene gain and loss of α-keratin typeI and typeII at each ancestral node in the phylogenetic tree. It was noted that loss of 6.5 genes in the ∼53–56 million years since the cetacean ancestor split from the cow/sheep clade. By contrast, the gain of five genes in the common ancestor of the cow/sheep. However, average one gene gain or lose in other marine mammals and their terrestrial relatives.
Hair follicle-specific keratin genes consisted of 9 hair follicle-specific epithelial keratin genes (K71-K75, K25-K28) and 17 hair keratins (K31-K40, K81-K86) expressed in hair follicle and hair which could directly affect hair growth . A large number of diseases of hair were caused by mutations in hair follicle-specific keratin genes. For example, a mutation in K75 leads to loose anagen hair syndrome, which is characterized by easily pluckable hair . In addition, the K71 mutation will make almost all types of mouse hair abnormal and easily fall off . Subsequently, increased gene loss and pseudogenization rate were found to focus on hair follicle-specific keratin genes in cetacean lineages (Fig. 1). Only an average of 14 hair follicle-specific keratins were identified in cetaceans but with an average of 26.33 in artiodactylans. Most importantly, the averaged pseudogenized rate in cetaceans was 43.88%, 11 times higher than that in artiodactylans (3.80%). In summary, the unique gene loss and/or loss of their function in hair follicle-specific keratin genes might contribute to the hairless phenotype.
Species-specific evolution of α-keratin gene repertoire
The α-keratin gene family was further classified according to the phylogenetic relationship reconstructed using ML and Bayesian approaches. We found that almost all of the α-keratin genes of 16 species could be clustered into respective clades according to each subfamily, which can be regarded as orthologous genes in phylogenetic tree (Figs. 2 and 3). However, some genes within one specie grouped together rather than with the members of the same family in other species, such as among the genes of K31, K33A, K33B and K34 in type I, as well as the K6A, K6B, K6C, K81, K83, K87 and K86 in the type II (Figs. 2 and 3). This phenomenon suggested that such members within a repetitive family did not evolve independently of each other but under concerted evolution , which was also found in other mammals such as horse, dog, human, etc. . We further used the RDP4 to test whether the gene conversion and recombination resulted in these 19 gene groups supposed to be under concerted evolution. The result showed that only three gene group supposed to be under concerted evolution was due to significant level of gene conversion and recombination (Additional file 2).
Species-specific evolution of α-keratin gene family was identified in the marine mammals, including differences in the total number of genes, functional genes, as well as pseudogenes. Nine keratin subfamilies were completely absent in cetaceans when compared to its terrestrial relative of artiodactylans, such as six type I genes (K10, K25, K27, K28, K33A, K33B) and three type II genes (K1, K73, K77) (Fig. 1). Most importantly, the cetacean species were also exclusively lost of 10 keratin genes (type I: K10, K25, K27, K28, K33A, K33B, K37; type II: K1, K73, K77) when compared with other marine mammals (Fig. 1). Besides these lost genes, we also examined 12 pseudogenes in cetaceans (type I: K9, K26, K34, K38, K39, K42; type II: K2, K6B, K74, K75, K79, K82, Fig. 1). Previous study has reported that deficient for K6B (epithelial keratins) may contributed to the absence of hair and nail in mice . By contrast, mutations in K1 and K10 were shown to be associated with bullous congenital ichthyosiform erythroderma (BCIE) . Collectively, cetaceans were exclusively lost K6B, K1 and K10, which is likely to contribute to its hair-loss phenotype to adapt to the completely aquatic environment. For complete aquatic manatee, however, we found that only four keratin subfamilies, such as K9, K6C, K6B, K39, were uniquely lost when compared with its terrestrial relatives, the elephant (Fig. 1).
When only cetaceans were considered, we found the average number of α-keratin genes in baleen whales was 43.5, slightly higher than that in toothed whales (averaged 34.8). Specially, pseudogenization rate of hair follicle-specific keratin genes was higher in toothed whales (42.86%) than in baleen whales (35.29%). The much higher number of intact functional keratin genes in the baleen whales may be associated with the presence of keratinized baleens.
Sharing keratin genes loss and pseudogenes in cetaceans and manatee as response for convergent hair-less phenotype
Marine mammals are relatively independent in evolution but with similar phenotype changes and associated physiological features to adaptation to aquatic lifestyle. Remarkably, both cetaceans and sirenians are fully aquatic and therefore encompass convergent hair-less phenotype to adapt to aquatic environment, since the lack of fur improves their hydrodynamic and subaquatic movements. Cetaceans have apparent reduction of the α-keratin gene repertoire and increased rate of pseudogenization while the manatees have a comparable number of α-keratin and pseudogenization rate when compared to their respective terrestrial relatives. Importantly, both groups were found to share common gene loss in the four subfamilies, i.e. K39, K9, K42, K74 that play a key role in hair development, likely to be responsible for hair-less phenotype. For example, it has been reported that heterozygous mutations in K74 could cause autosomal dominant woolly hair, and/or hypotrichosis simplex  whereas mutations in K9 gene caused epidermolytic palmoplantar keratoderma, leading to manifests cytolysis and epidermal thickening . In addition, K39 as a member of the hair-keratins, was identified to play a key role in late hair differentiation . By contrast, when another components of hair, KRTAP, are considered, a significant gene number reduction (81 in manatee vs. 112 in elephant) and increased pseudogenization rate (41.98% in manatees vs. 20% in elephant) were detected in manatee. Similarly, it has been recently reported that dolphins only have 35 KRTAP genes (compared to 145 genes in the cow) but with highest pseudogenization rate (74% relative to the 19% of mammalian average) . However, we have detected average of 65 functional KRTAP genes in the pinnipeds that have abundant hair, which were comparable to that of terrestrial relatives, such as panda with 64 KRTAP genes . Collectively, a reduced number of functional KRTAP genes, high percentage of KRTAP pseudogenes, and common keratin gene loss in both cetaceans and manatees might be responsible for their convergent hair-less phenotype.
In this study, we provide a comprehensive characterization of α-keratin genes among marine mammals and shed light on the mechanisms involved in the evolution of this gene family. Our results show an apparent reduction in the α-keratin gene repertoire and increased rate of pseudogenization in cetaceans when compared to terrestrial counterparts, which may be associated with their hairless phenotype. In contrast, the evolution pattern of α-keratin genes in pinnipeds was comparable to that of their respective terrestrial relatives, which is well matched with its fur coat. Interestingly, we found a reduced number of functional KRTAP genes, high percentage of KRTAP pseudogenes, and common keratin gene loss in both cetaceans and manatees, which might be as a result of convergent hair-loss phenotype to inhabit completely aquatic environment. Species-specific evolution of α-keratin gene repertoire identified in marine mammals was likely to contribute to their hair diversities phenotype although their life or part inhabiting in aquatic environment. Our study first indicated the genetic basis of hair diversities of marine mammals.
Bullous congenital ichthyosiform erythroderma
Eukaryotic translation initiation factor 1
Eukaryotic translation initiation factor 4B
Fas apoptotic inhibitory molecule 2
SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily e, member 1
Uhen MD. Evolution of marine mammals: back to the sea after 300 million years. Anat Rec (Hoboken). 2007;290:514–22.
Berta A, Sumich JL, Kovacs KM. Marine mammals: evolutionary biology. Cambridge: Academic Press. 2005.
Bergman J. Why mammal body hair is an evolutionary enigma. CRS Quarterly. 2004;40:240–3.
Maderson PFA. Mammalian skin evolution: a reevaluation. Exp Dermatol. 2003;12:233–6.
Berta A. Pinnipedia, overview. J Zool. 2002;83:1525–31.
Drake SE, Crish SD, George JC, Stimmelmayr R, Thewissen JG. Sensory hairs in the bowhead whale, Balaena Mysticetus (Cetacea, Mammalia). Anat Rec (Hoboken). 2015;298:1327–35.
Yonenaga A. Future man-made fiber. High-performance and specialty fibers. Japan: Springer; 2016. p. 435–51.
Riedman M, Estes JA. The sea otter (Enhydra Lutris): behavior, ecology, and natural history. Biol Reprod. 1990;90:1–126.
Shimomura Y, Ito M. Human hair keratin-associated proteins. J Investig Dermatol Symp Proc. 2005;10:230–3.
Moll R, Divo M, Langbein L. The human keratins: biology and pathology. Histochem Cell Biol. 2008;129:705–33.
Schweizer J, Bowden PE, Coulombe PA, Langbein L, Lane EB, Magin TM, Maltais L, Omary MB, Parry DA, Rogers MA, Wright MW. New consensus nomenclature for mammalian keratins. J Cell Biol 2006;174:169-174.
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.
Khan I, Maldonado E, Vasconcelos V, O'Brien SJ, Johnson WE, Antunes A. Mammalian keratin associated proteins (KRTAPs) subgenomes: disentangling hair diversity and adaptation to terrestrial and aquatic environments. BMC Genomics. 2014;15:779.
Nery MF, Arroyo JI, Opazo JC. Increased rate of hair keratin gene loss in the cetacean lineage. BMC Genomics. 2014;15:869.
Altschul SF, Madden TL, Schäffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.
Cheng J, Liu C, Koopman WJ, Mountz JD. Characterization of human Fas gene. Exon/intron organization and promoter region. J Immunol. 1995;154:1239–45.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312.
Ronquist F, Huelsenbeck JP. MrBrayes 3: Bayesian phylogentetic inference under mixed models. Bioinformatics. 2003;19:1572–4.
Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015;1:vev003.
Han MV, Thomas GW, Lugo-Martinez J, Hahn MW. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol Biol Evol. 2013;30:1987–97.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
Reep RL, Stoll ML, Marshall CD, Homer BL, Samuelson DA. Microanatomy of facial vibrissae in the Florida manatee: the basis for specialized sensory function and Oripulation. Brain Behav Evol. 2001;58:1–14.
Schweizer J, Langbein L, Rogers MA, Winter H. Hair follicle-specific keratins and their diseases. Exp Cell Res. 2007;313:2010–20.
Chapalain V, Winter H, Langbein L, Le Roy JM, Labrèze C, Nikolic M, Schweizer J, Taïeb A. Is the loose anagen hair syndrome a keratin disorder? A clinical and molecular study. Arch Dermatol. 2002;138:501–6.
Peters T, Sedlmeier R, Büssow H, Runkel F, Lüers GH, Korthaus D, Fuchs H, Hrabé de Angelis M, Stumm G, Russ AP, Porter RM, Augustin M, Franz T. Alopecia in a novel mouse model RCO3 is caused by mK6irs1 deficiency. J Invest Dermatol. 2003;121:674–80.
Wojcik SM, Longley MA, Roop DR. Discovery of a novel Murine keratin 6 (K6) Isoform explains the absence of hair and nail defects in mice deficient for K6a and K6b. J Cell Biol. 2001;154:619–30.
McLean WH, Eady RA, Dopping-Hepenstal PJ, McMillan JR, Leigh IM, Navsaria HA, Higgins C, Harper JI, Paige DG, Morley SM. Mutations in the rod 1A domain of keratins 1 and 10 in bullous congenital ichthyosiform erythroderma (BCIE). J Invest Dermatol. 1994;102:4–30.
Raykova D, Klar J, Azhar A, Khan TN, Malik NA, Iqbal M, Tariq M, Baig SM, Dahl N. Autosomal recessive transmission of a rare KRT74 variant causes hair and nail ectodermal dysplasia: allelism with dominant woolly hair/hypotrichosis. PLoS One. 2014;9:e93607.
Torchard D, Blanchet-Bardon C, Serova O, Langbein L, Narod S, Janin N, Goguel AF, Bernheim A, Franke WW, Lenoir GM. Epidermolytic palmoplantar keratoderma cosegrates with a keratin 9 mutation in a pedigree with breast and ovarian cancer. Nat Genet. 1994;6:106–10.
Langbein L, Rogers MA, Praetzel-Wunder S, Bockler D, Schirmacher P, Schweizer J. Novel type I hair keratins K39 and K40 are the last to be expressed in differentiation of the hair: completion of the human hair keratin catalog. J Invest Dermatol. 2007;127:1532–5.
We thank members of the Jiangsu Key Laboratory for Biodiversity and Biotechnology, Nanjing Normal University, for their contributions to this paper. A special thank-you is also due from Zhengfei Wang and Professor Qi Wu for their technical support.
This work was supported by the National Natural Science Foundation of China (NSFC, grant number 31570379 to S X), the National Science Fund for Distinguished Young Scholars to GY (grant number 31325025), and the State Key Program of National Natural Science of China to GY (grant number 31630071); the National Key Programme of Research and Development, Ministry of Science and Technology (Grant number 2016YFC0503200 to GY and SX), the Priority Academic Program Development of Jiangsu Higher Education Institutions to GY. and S X, the Natural Science Foundation of Jiangsu Province of China (grant number BK20141449) to SX.
Availability of data and materials
Raw data are shown in Additional file 1.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.