Open Access

Characterization of hairless (Hr) and FGF5genes provides insights into the molecular basis of hair loss in cetaceans

BMC Evolutionary Biology201313:34

DOI: 10.1186/1471-2148-13-34

Received: 17 November 2012

Accepted: 30 January 2013

Published: 9 February 2013

Abstract

Background

Hair is one of the main distinguishing characteristics of mammals and it has many important biological functions. Cetaceans originated from terrestrial mammals and they have evolved a series of adaptations to aquatic environments, which are of evolutionary significance. However, the molecular mechanisms underlying their aquatic adaptations have not been well explored. This study provided insights into the evolution of hair loss during the transition from land to water by investigating and comparing two essential regulators of hair follicle development and hair follicle cycling, i.e., the Hairless (Hr) and FGF5 genes, in representative cetaceans and their terrestrial relatives.

Results

The full open reading frame sequences of the Hr and FGF5 genes were characterized in seven cetaceans. The sequence characteristics and evolutionary analyses suggested the functional loss of the Hr gene in cetaceans, which supports the loss of hair during their full adaptation to aquatic habitats. By contrast, positive selection for the FGF5 gene was found in cetaceans where a series of positively selected amino acid residues were identified.

Conclusions

This is the first study to investigate the molecular basis of the hair loss in cetaceans. Our investigation of Hr and FGF5, two indispensable regulators of the hair cycle, provide some new insights into the molecular basis of hair loss in cetaceans. The results suggest that positive selection for the FGF5 gene might have promoted the termination of hair growth and early entry into the catagen stage of hair follicle cycling. Consequently, the hair follicle cycle was disrupted and the hair was lost completely due to the loss of the Hr gene function in cetaceans. This suggests that cetaceans have evolved an effective and complex mechanism for hair loss.

Background

Hair is one of the distinguishing characteristics of mammals and it has many important biological functions [1]. Hair is produced by hair follicles (HFs), which are complex mini-organs in the skin that are formed during embryonic development (morphogenesis). New hair is generated continuously throughout life because the postnatal HFs experience cyclic phases of active growth (anagen stage), regression (catagen stage), and inactivity (telogen stage) [2]. Many genes and signaling pathways are involved in HF development [2, 3], including the hairless (Hr) gene and fibroblast growth factor 5.

The Hr gene is significantly expressed in skin and it encodes a putative zinc finger transcription factor of approximately 130 kDa [4]. Hr is a candidate gene that regulates basic HF functions [5]. More detailed biochemical analyses of the function of the encoded protein have shown that Hr is a transcriptional corepressor that interacts with nuclear receptors, including thyroid hormone receptor (TR), retinoic acid orphan receptor α (RORα) and vitamin D receptor (VDR), to regulate specific target genes involved with hair morphogenesis and HF cycling [6, 7]. It appears that Hr functions during the cellular transition to the first adult hair cycle because hair growth ceases completely in its absence, which results in a form of inherited total alopecia [8].

The fibroblast growth factor 5 (designated as Fgf5 in rats and mice, and FGF5 in other mammals), comprises three exons and it is an essential regulator of the HF development and cycling [9]. Recent studies have also suggested that the FGF5 gene is associated with hair length and it controls the cessation of the anagen stage [2, 1013].

Cetaceans (whales, dolphins and porpoises) are ecologically diverse and they inhabit waters that range from coastal to oceanic and from tropical to polar [14]. Numerous paleontological, morphological, embryological, and molecular studies have suggested that cetaceans evolved from terrestrial mammals [1519]. The transition from land to water and their subsequent adaptation to completely aquatic habitats make cetaceans remarkable and evolutionarily significant, although few studies have investigated the molecular basis of this process [2026].

Cetaceans generally lack a coat of hair, although some cetacean species retain a few hairs on their face while the fetus has whiskers in others. This is probably an adaptation that reduces friction and improves locomotion. However, the precise molecular mechanisms underlying hair loss are unclear. Given the important roles of the Hr and FGF5 genes during HF morphogenesis and HF cycling, the current study determined the full open reading frame (ORF) sequences of these two genes in seven representative cetacean species and compared them with orthologous sequences from terrestrial mammals. The goal was to determine whether evolutionary changes in these two genes were associated with the transition from land to water, and the hair loss of cetaceans during this adaptive process. To the best of our knowledge, this is the first study to investigate the molecular basis of hair loss in cetaceans.

Results

Hairless (Hr) and FGF5genes in cetaceans

Complete ORF sequences were determined for the Hr and FGF5 genes in seven cetaceans. As shown in Additional file 1, the cetacean Hr gene contained 18 exons and detailed information on each exon is provided in Table 1. They all shared the typical features of mammalian Hr genes and alignments of the deduced amino acid sequences of the cetacean Hr genes are shown in Additional file 2. No frame-shift mutations or premature stop codons were detected in cetaceans, although a series of apparent deletions and specific amino acid changes were found in important functional domains of the cetacean Hr genes (Figure 1, Additional files 2 and 3). In contrast to the toothed whale, the two baleen whales (Balaenoptera omurai and B. acutorostrata) lacked insertions/deletions (indels), so they had intact ORFs (Additional file 2).
Table 1

Exon organization of the seven Cetacean hairless ( Hr ) genes obtained in this study

Exon

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

Tursiops truncatus

612

784

151

188

165

90

116

82

164

219

166

70

131

120

116

165

129

63

Delphinus capensis

612

784

151

188

165

90

116

82

164

219

166

70

131

120

116

165

129

63

Neophocaena phocaenoides

612

784

151

188

165

90

116

82

164

213

166

70

131

120

116

165

129

63

Delphinapterus leucas

612

784

151

188

165

90

116

82

164

219

166

70

131

120

116

165

129

63

Lipotes vexillifer

612

784

151

188

165

90

116

82

164

237

166

70

131

120

116

165

129

63

Balaenoptera omurai

612

790

151

188

165

90

116

82

164

237

166

70

131

120

116

165

129

63

Balaenoptera acutorostrata

612

790

151

188

165

90

116

82

164

237

166

70

131

120

116

165

129

63

https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-13-34/MediaObjects/12862_2012_Article_2295_Fig1_HTML.jpg
Figure 1

Distribution of mutations in the three-dimensional structure of cetacean Hr . The RD3 and JmjC domains are colored yellow and green, respectively. Locations marked in red coincide with mutations found within cetaceans and toothed whales within the RD3 and JmjC regions.

The cetacean FGF5 also contained an uninterrupted ORF, with no premature stop codons. Three exons were identified based on an alignment of the cetacean FGF5 sequences with those from other mammals. Exons 2 and 3 were highly conserved in all of the mammals examined in this study (Additional files 1 and 4). As shown in Additional file 4, the FGF5 gene encoded a protein containing approximately 270 amino acid residues, including a signal peptide with 20 amino acid residues.

Phylogenetic analysis

The maximum likelihood and Bayesian analyses of all datasets yielded similar tree topologies (Additional files 5 and 6). This strongly supported the nesting of Cetacea within Artiodactyla and the monophyly of Cetartiodactyla. Overall, the relationships among the nine placental mammalian orders examined were in agreement with those reported in previous studies, except for Perssodactyla (with horse as a representative) and the sister relationship between Odontoceti (toothed whales) and Mysticeti (baleen whales) within Cetacea (e.g., [27, 28]).

Relaxation of selection for cetacean Hrgenes

A series of evolutionary models were performed with the likelihood framework to analyze the selective constraints on Hr, using the species tree shown in Figure 2 as the working topology. One-ratio model analyses of all mammals (dataset І: 17 sequences) showed that all of the branches in Figure 2 shared the same estimated ω of 0.2791 (Model A in Table 2), which indicated the existence of strong functional constraints on mammalian Hr genes. In the two-ratio model analyses (model B in Table 2), the ω value of the focal branch was 0.40432 and model B fitted the data significantly better than a one-ratio model, which assumed a single ω for all branches (P = 0.024, Table 2). A comparison of model B and model C (ω2 is fixed at 1, Table 2) also showed that ω2 was significantly less than 1, which suggested that the relaxation of the functional constraint on the Hr gene did not occur immediately after the common cetacean ancestor diverged from the terrestrial mammals.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-13-34/MediaObjects/12862_2012_Article_2295_Fig2_HTML.jpg
Figure 2

The ω values of FGF5 genes in distinct evolutionary lineages of cetaceans and other mammals using a phylogenetic tree derived from Chen et al. [29], Gatesy et al. [27], and Zhou et al. [19],[28]. Branches a-h relate to those shown in Table 3. The ω values of individual branches shown are based on the free-ratio model. In some cases, zero synonymous substitutions produced an ω value of infinity (n.a.). The estimated numbers of nonsynonymous and synonymous changes are shown in parentheses. The ω values of branches a, b, and c (marked with different colors) were estimated using the two-ratio models.

Table 2

Likelihood ratio tests using various models to determine the selective pressures on the cetacean Hr gene

Models

ω(d N /d S )

lnLa

npb

Models compared

2Δ(lnL)c

P-values

Dataset I: 17 sequences (16 terrestrial mammals plus the ancestral cetacean sequence)

      

A. All branches with the same ω

ω = 0.27908

−24914.624502

34

   

B. Branch a (ancestral to cetacean) with ω2, other branches with ω1

ω1 = 0.27502, ω2 = 0.40432

−24912.059374

35

A vs B

5.130256

0.024

C. Branch a with ω2 = 1, other branches with ω1

ω1 = 0.27427, ω2 = 1

−24925.508371

34

B vs C

26.89799

2.14E-07

Dataset II: 19 sequences (16 terrestrial mammals plus two baleen whales and the ancestral sequence of toothed whales)

      

D. All branches with the same ω

ω = 0.28065

−25530.203259

38

   

E. Branch b (ancestral to toothed whale) with ω2, other branches with ω1

ω1 = 0.27911, ω2 = 0.48639

−25528.780477

39

D vs E

2.845564

0.092

F. Branch b with ω2 = 1, other branches with ω1

ω1 = 0.27890, ω2 = 1

−25530.928574

38

E vs F

4.296194

0.038

Dataset III: 22 sequences (16 terrestrial mammals plus five toothed whales and the ancestral sequence of baleen whales)

      

G. All branches with the same ω

ω = 0.28417

−26239.723434

44

   

H. Branch c (ancestral to baleen whale) with ω2, other branches with ω1

ω1 = 0.28182, ω2 = 0.62531

−26236.569206

45

H vs G

6.308456

0.012

I. Branch c with ω2 = 1, other branches with ω1

ω1 = 0.28170, ω2 = 1

−26237.557075

44

H vs I

1.975738

0.160

Note: a, ln L is the log-likelihood score; b, number of parameters; c, twice the difference in ln L between the two models compared.

Interestingly, both the common ancestors of toothed whales and baleen whales had a higher ω compared with other mammals (model E and H; see Table 2). In addition, the likelihood ratio test (LRT) results detected no significant difference in the estimates using the two-ratio model where ω was not fixed for the branch of the common ancestor of baleen whales and that where it was fixed to 1 (P = 0.160; model H vs. I in Table 2), which suggested that the Hr gene was close to selective neutrality in baleen whales. Furthermore, a series of sites were polymorphic in cetaceans whereas they were highly conserved (monomorphic) in other mammals (Additional file 3). Overall, these results suggest that the cetacean Hr gene has experienced a significant relaxation of selection during its evolution from terrestrial mammals and subsequent lineage diversification.

Positive selection for the cetacean FGF5gene

In the branch-specific model analyses, the ω ratio calculated in the one-ratio model (M0) was 0.16615 (Table 3), which indicated a generally strong purifying selection for the mammalian FGF5 gene, whereas only the lineage leading to the common ancestor of toothed whales (branch b in Figure 2) contained significant evidence of positive selection (2ΔL=5.93932, P = 0.0148). In addition, the lineage leading to the common ancestor of baleen whales (branch c in Figure 2) had a higher ω value (ω1 = 0.93251) than the background (ω0 = 0.16491), although the two-ratio models did not fit significantly better than model M0 (Table 3). Significant LRT statistics and positively selected sites were detected in the lineage leading to the common ancestor of toothed whales (branch b: 2ΔL=9.50702, P = 0.0086) in the branch-site model, whereas most lineages outside the cetaceans showed no evidence of positive selection (Table 3 and Figure 2).
Table 3

Likelihood values and parameter estimates for the FGF5 gene

Models

ln La

Estimate of parameters

2ΔLb(P-value)

Positively selected sites

M0:one-ratio

−4450.77242

ω = 0.16615

  

Branch-specific models

    

Branch a (ancestral Cetacea)

    

Two-ratio

−4450.74677

ω0 = 0.16656, ω1 = 0.14336

0.0513

 

(P = 0.8208)

Branch b (ancestral Odontoceti)

    

Two-ratio

−4447.80276

ω0 = 0.16253, ω1 = 1.04812

5.93932

 

( P = 0.0148)

Branch c (ancestral Mysticeti)

    

Two-ratio

−4449.912614

ω0 = 0.16491, ω1= 0.93251

1.719612

 

(P = 0.189743)

Branch d (ancestral Cetartiodactyla)

    

Two-ratio

−4449.692636

ω0 = 0.16802, ω1 = 0.0001

2.159568

 

(P = 0.1416845)

Branch e (ancestral Carnivora)

    

Two-ratio

−4450.028289

ω0 = 0.16284, ω1 = 0.28977

1.488262

 

(P = 0.2224863)

Branch f (ancestral Chiroptera)

    

Two-ratio

−4450.762558

ω0 = 0.16637, ω1 = 0.14605

0.019724

 

(P = 0.8883106)

Branch g (ancestral Rodentia)

    

Two-ratio

−4450.367662

ω0 = 0.16377, ω1 = 0.33139

0.809516

 

(P = 0.3682634)

Branch h (ancestral Primates)

    

Two-ratio

−4450.765224

ω0 = 0.16585, ω1 = 0.17607

0.014396

 

(P = 0.9044964)

Branch-site models

    

Null

−4370.04122

p0 = 0.80258, p1 = 0.19742, ω0 = 0.05816, ω1 = 1

  

Branch a (ancestral Cetacea)

    

Alternative

−4370.04122

p0 = 0.80258, p1 = 0.19742, p2a = 0, p2b = 0, ω0 = 0.05816, ω1 = 1, ω2 = 1

0 (P = 1)

None

Branch b (ancestral Odontoceti)

    

Alternative

−4365.28771

p0 = 0.46391, p1 = 0.12028, p2a = 0.33019, p2b = 0.08561, ω0 = 0.05123, ω1 = 1, ω2 = 2.03039

9.50702

K24E(0.945);P34Q(0.914);S68P(0.535);S74N(0.947);S84T(0.94)

( P = 0.0086)

Branch c (ancestral Mysticeti)

    

Alternative

−4368.231030

p0 = 0, p1 = 0, p2a = 0.80784, p2b = 0.19216, ω0 = 0.05877, ω1 = 1, ω2 = 1

3.62038

M119T(0.777);V235A(0.937)

(P = 0.163623)

Branch d (ancestral Cetartiodactyla)

    

Alternative

−4370.041222

p0 = 0.80258, p1 = 0.19742, p2a = 0, p2b = 0,

0 (P = 1)

None

ω0 = 0.05816, ω1 = 1, ω2 = 1

Branch e (ancestral Carnivora)

    

Alternative

−4370.041222

p0 = 0.80258, p1 = 0.19742, p2a = 0, p2b = 0,

0 (P = 1)

186T(0.571)

ω0 = 0.05816, ω1 = 1, ω2 = 1

Branch f (ancestral Chiroptera)

    

Alternative

−4369.925398

p0 = 0.75669, p1 = 0.18750, p2a = 0.04473, p2b = 0.01108, ω0 = 0.05728, ω1 =1, ω2 =1

0.231644

84T(0.736)

(P = 0.8906337)

Branch g (ancestral Rodentia)

    

Alternative

−4370.041222

p0 = 0.80258, p1 = 0.19742, p2a =0, p2b =0,

0 (P = 1)

43G(0.621)

ω0 = 0.05816, ω1 = 1, ω2 = 1

Branch h (ancestral Primates)

    

Alternative

−4370.041222

p0 = 0.80258, p1 = 0.19742, p2a = 0, p2b = 0,

0 (P = 1)

None

  

ω0 = 0.05817, ω1 = 1, ω2 = 1

  

Note: a, ln L is the log-likelihood score; b, likelihood ratio test (LRT) to detect positive selection.

Five codons were shown to be under positive selection in the branch leading to the common ancestor of toothed whales (branch b in Figure 2) according to the branch-site model (Table 3), and four of these were shown to have undergone radical changes (Table 4). These positively selected amino acids did not correspond to residues known to interact with the FGF receptor (FGFR) and heparin (data not shown), but many of them were involved in or near a region rich in O-glycosylation and N-glycosylation sites, near the signal peptide region, or the FGFR-binding sites (Additional file 4).
Table 4

FGF5 candidate amino acid sites under positive selection identified in toothed whales

Amino acid position

PAML (branch-site model A)

Radical change

Conservative change

Change

24

0.945

Lys-Glu

Charge, polarity, and volume

34

0.914

Pro-Gln

Polarity

68

0.535

Ser-Pro

Polarity

74

0.947

Ser-Asn

Polarity and volume

84

0.94

Ser-Thr

 

Discussion

Relaxed selection for the cetacean Hrgene suggests its functional loss

The Hr gene is highly conserved and has traditionally been regarded as strongly functionally constrained during mammalian evolution because of its functional significance in HF cycling and the important roles of hair in mammals [1, 5] (Additional file 3 and Table 2). However, our data suggest that the cetacean Hr gene may have experienced a relaxation of selective pressure to become a pseudogene (Table 2). Pseudogenes that experience relaxed selective pressure are expected to have a higher ω ratio compared with functional genes, which are usually under purifying selection [22]. Our analyses of the ω ratio based on different datasets showed that the ω ratios for the Hr sequences of toothed whales, baleen whales, and all cetaceans were clearly higher than those of putative functional sequences in other mammals (Table 2). More importantly, the presence of a series of polymorphic (as opposed to conservative or monomorphic) sites in cetaceans may be further evidence for the relaxed selection of cetacean Hr genes. In addition, most of the mutations in cetacean Hr genes were found in or near important functional domains or conserved regions, which was also supported by the homology modeling analysis of of cetacean Hr genes (Figure 1, and Additional file 3). For example, H996Q and G1036A mutations in the JmjC domain probably have led to the loss of histone demethylase activity and constitutive methylation, which would have promoted the transcriptional repression of Hr-interacting signaling pathways genes and ultimately hair loss [7, 30]. This was corroborated by many natural mutants in the JmjC domain, which cause hair loss in mice and humans (Figure 3). Many mutations, including deletions, insertions, nonsense, missense and splice-sites, in the functional domains of the Hr genes from human patients or rodent models have also been reported to cause congenital atrichia (hair loss) (reviewed in [7]). Many, if not all, of the mutations (especially six base pair (bp) and 18-bp deletions in toothed whales) found in cetaceans are predicted to disrupt the local secondary structure and produce defects in HF regeneration, as found in humans and mouse models.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-13-34/MediaObjects/12862_2012_Article_2295_Fig3_HTML.jpg
Figure 3

Schematic representation of rat Hr structural and functional domains (reproduced from Thompson [7]; Hsieh et al. [31]) (A), and natural mutants of Hr that cause hair loss in mice and humans (B). The positions of the three repression domains (RD1, RD2, RD3), and JmjC domain are indicated yellow and red, respectively. Nuclear localization signal (NLS), and zinc-finger domain are boxed in black [7, 31]. Interacting domains 1 and 2 (ID1 and ID2) mediating the interaction between Hr and the retinoic acid receptor-related orphan receptor-alpha (RORα) [32] or with the thyroid hormone receptor (TR) [33] are boxed in green and blue, respectively. Specific cetacean, toothed whale, baleen whale, and Delphinoidea mutations are highlighted in red, blue, pink, and green, respectively [3441].

Furthermore, mutations leading to the absence of Hr function were found to differ to some extent in toothed and baleen whales, e.g., the 6-bp deletion in exon 2 was only present in toothed whales. Mutations in the cetacean Hr gene also exhibited a taxa-specific pattern. For example, the four delphinoids examined in this study all had an identical 18-bp deletion in exon 10. Another 6-bp deletion was identified in exon 10 of the finless porpoise. These different mutations in the Hr gene in different cetaceans may have resulted in a similar hair loss phenotype, which agrees with the fact that different mutations produce the same form of alopecia (hair loss) in humans and mice [7, 42].

In summary, this study suggests that the cetacean Hr gene has undergone various evolutionary changes that probably correspond to its loss of function. These findings are consistent with morphological evidence that adult whales have no body hair covering, whereas they have hair during their early development [14]. During the transition from land to water, the ancestor of cetaceans was faced with different habitats and survival challenges, and the presence of hair may have been a hindrance to locomotion. Therefore, hair was probably unnecessary for cetaceans based on the cetacean Hr gene.

Positive selection for the cetacean FGF5gene

In contrast to the relaxed selection for the Hr gene, the cetacean FGF5 gene was under strong positive selection, according to the significantly higher ω value on the branch leading to the toothed whale compared with the background and the large number of specific codons detected by the branch-site models (Tables 3). The ω value (ω1 = 0.93251) of the lineage leading to the common ancestor of baleen whales was not significantly higher than 1, but this value was higher than the background value (ω0 = 0.16491) (Tables 3). This elevated ω estimate for the FGF5 gene relative to other branches suggests the accelerated evolution of the FGF5 gene in baleen whales. However, no evidence of positive selection was detected in the common ancestor of Cetacea (Tables 3), which suggests that positive selection occurred after the Odontoceti-Mysticeti split. Furthermore, a series of potentially important adaptive amino acid changes were detected in toothed whales (Table 3 and Additional file 4) and most of these codon changes had radical effects on their physicochemical properties (charge, polarity, and volume) (Table 4). In general, more radical amino acid substitutions have greater functional effects during evolution [43]. In addition, these residues were located in or near a region rich in O-glycosylation and N-glycosylation sites, near the signal peptide region or near the FGFR-binding sites. Therefore, these amino acid changes may have affected the secondary or tertiary conformation of the FGF5 molecule and ultimately affected its function [44, 45].

Accumulating evidence suggests that the FGF5 gene controls the cessation of the anagen stage of the HF cycle and pathogenic mutations in this gene are known to produce long hair in phenotypic variants of many mammals [9, 10, 12, 13, 46]. It has been suggested that the epithelium and underlying mesenchyma interact in utero to form HFs (hair morphogenesis), before the HFs enter their three stages (anagen, catagen, and telogen) [3, 47], which is consistent with embryological evidence that cetaceans develop body hair in the womb [17, 48]. However, no extant whales retain any body hair after birth, with the exception of some snout hairs and hairs around the blowholes that act as sensory bristles in some baleen whales.

Humans and mice with mutations in the Hr gene develop apparently normal HFs but shed their hair completely soon after birth [4, 42]. However, cetaceans lose their body hair in the womb rather than after birth. In addition, no positive selection was detected in FGF5 genes in other branches, including the human and mouse genes. Overall, this study suggests that positive selection for FGF5 genes in toothed whales may have played an important role in terminating the hair growth cycle and accelerated entry into the catagen stage of hair growth. Interestingly, the current study showed that the Hr gene may have lost its function in cetaceans. This gene loss initiates a premature and abnormal catagen stage, which leads to the destruction of the normal HF architecture and abrogates the HF’s ability to cycle. This hypothesis may be tested by elucidating the molecular mechanism of hair loss in cetaceans and the differences in hair loss between cetaceans and Hr mutants in humans and mice. However, it is not easy to obtain cetacean biopsies for expression analyses, and so further research (e.g., expression experiment or histochemistry) will be required to test our hypotheses of Hr pseudogenisation and FGF5 positive selection in cetaceans. Further samples should also be obtained from cetaceans (especially baleen whales) and other fully aquatic marine mammals (e.g., sirenians: manatees and dugongs) to test this hypothesis and elucidate the molecular basis of hair loss.

Conclusions

The data presented in this study suggest that the cetacean Hr gene has undergone evolutionary changes related to its loss of function. By contrast, positive selection on the FGF5 gene was detected in cetaceans, including a series of positively selected amino acid residues. The evolutionary changes in these two genes may provide new insights into the molecular basis of significant hair loss in cetaceans during their transition from land to water. Many signaling pathways and factors are known to be involved in the regulation of HF morphogenesis and the HF cycle [2, 3]. Therefore, additional genes related to hair development should be investigated to improve our understanding of the molecular mechanism underlying hair loss in fully aquatic marine mammals.

Methods

Taxonomic coverage

Seven cetacean species (five odontocetes and two mysticetes) were sequenced during this study (Table 5). In addition, the full-length ORFs of FGF5 and Hr from 16 other mammals were searched and downloaded from the Ensemble Genome database (http://www.ensembl.org) and GenBank (http://www.ncbi.nlm.nih.gov) (Table 5). When multiple splice vari-ants of FGF5 and Hr were available, only the full-length form from a species was used in subsequent analyses.
Table 5

List of taxonomic samples and sequences used in this study

Order

Suborder

Species

Accession number

  

Family

Scientific name

Common name

Hr

FGF5

Cetacea

Odontoceti

Delphinidae

Tursiops truncatus

Common bottlenose dolphin

KC140206

KC140213

   

Delphinus capensis

Long-beaked common dolphin

KC140205

KC140212

  

Monodontidae

Delphinapterus leucas

Beluga or white whale

KC140208

KC140215

  

Phocoenidae

Neophocaena phocaenoides

Finless porpoise

KC140207

KC140214

  

Lipotidae

Lipotes vexillifer

Yangtze river dolphin or baiji

KC140209

KC140216

 

Mysticeti

Balaenopteridae

Balaenoptera acutorostrata

Common minke whale

KC140210

KC140217

   

Balaenoptera omurai

Omura's whale

KC140211

KC140218

Artiodactyla

 

Bovidae

Bos taurus

Cow

Genome sequence

Genome sequence

  

Suidae

Sus scrofa

Pig

Genome sequence

Genome sequence

Perissodactyla

 

Equidae

Equus caballus

Horse

Genome sequence

Genome sequence

Carnivora

 

Canidae

Canis familiaris

Dog

Genome sequence

Genome sequence

  

Ursidae

Ailuropoda melanoleuca

Giant panda

XM_002914796.1

XM_002912480.1

Chiroptera

 

Pteropodidae

Pteropus vampyrus

Megabat

Genome sequence

Genome sequence

  

Vespertilionidae

Myotis lucifugus

Microbat

Genome sequence

Genome sequence

Eulipotyphla

 

Erinaceidae

Erinaceus europaeus

Hedgehog

Genome sequence

Genome sequence

Rodentia

 

Muridae

Mus musculus

Mouse

Genome sequence

Genome sequence

   

Rattus norvegicus

Rat

Genome sequence

Genome sequence

  

Caviidae

Cavia porcellus

Guinea pig

Genome sequence

Genome sequence

Lagomorpha

 

Leporidae

Oryctolagus cuniculus

Rabbit

Genome sequence

Genome sequence

Primates

 

Hominidae

Homo sapiens

Human

Genome sequence

Genome sequence

  

Cercopithecidae

Macaca mulatta

Macaque

Genome sequence

Genome sequence

  

Callitrichidae

Callithrix jacchus

Marmoset

Genome sequence

Genome sequence

Proboscidea

 

Elephantidae

Loxodonta africana

Elephant

Genome sequence

Genome sequence

Amplification and sequencing of cetacean FGF5 and Hrgenes

Primers were designed for the conserved regions based on an alignment of genomic data from the cow Bos taurus (http://asia.ensembl.org/Bos_taurus/Info/Index) and bottlenose dolphin (Tursiops truncatus) (http://asia.ensembl.org/Tursiops_truncatus/Info/Index). The primer information is available upon request. Total genomic DNA was extracted from the muscle tissues using a standard phenol/chloroform procedure followed by ethanol precipitation [49]. For blood, we used the DNAeasy Blood Extraction Kit (Qiagen) in a separate laboratory facility. All PCR amplification were conducted using a BioRAD PTC-200 with 2×EasyTaq PCR SuperMix (TransGen Biotech) and the following profile: 34 cycles at 94°C for 5 min, 94°C for 30 s, 53°C −59°C for 30 s, and 72°C for 30 s, followed by a 10 min extension at 72°C. The amplified PCR products were purified and sequenced in both directions using an ABI 3730 automated genetic analyzer. Novel sequences were deposited in GenBank under accession numbers KC140205-KC140218.

Sequence alignment and statistical analyses

Nucleotide sequences from the coding regions of Hr and FGF5, and their deduced amino acid sequences, were aligned separately using the program CLUSTAL X [50] and manually adjusted with GeneDoc. The nucleotide sequence alignment was generated based on the protein sequence alignment. Ancestral sequences were reconstructed using Bayesian method [51], which was implemented in the BASEML program in PAML 4.5 [52]. A three-dimensional domain structure of the cetacean Hr was predicted using SWISS-MODEL (http://swissmodel.expasy.org) [5355].

Phylogenetic trees were reconstructed using maximum likelihood algorithms in MetaPIGA 2.0 [56] and Bayesian inference (BI) in MrBayes 3.1.2 [57] for each gene independently and for a combined dataset of FGF5 and Hr, using the African elephant (Loxodonta africana) as the outgroup. MRMODELTEST 2.3 [58] was used to select the optimal models for each partition based on the Akaike Information Criterion (AIC). Maximum Likelihood analyses were performed using MetaPIGA 2 [56] with 1000 replicate metaGA searches. The Bayesian analyses of the nucleotide matrix were performed using a codon model (general time irreversible (GTR) gamma invariant model) or mixed models (the GTR gamma invariant model for the first and second codon positions, and the GTR gamma model for the third codon position). Four Markov chains were run for 20 million generations in MrBayes 3.1.2, with sampling every 1000 generations. The stationarity of the likelihood scores of the sampled trees were checked using Tracer 1.4 [59]. The Bayesian posterior probabilities (PP) were obtained from the 50% majority rule consensus of the post burn-in trees sampled at stationarity, after removing the first 10% of trees as a “burn-in” stage. The aligned sequences and phylogenetic trees were deposited in TreeBase (http://purl.org/phylo/treebase/phylows/study/TB2:S13758).

Analysis of selective pressure

Analyses of selective pressure were carried out using a codon-based maximum likelihood method implemented in the CODEML program in the PAML 4.5 package [52]. A consensus tree that included all of the species employed in the present study was inferred from Chen et al. [29], Gatesy et al. [27], and Zhou et al. [19, 28], and used in the subsequent PAML analyses (Figure 2). In all PAML-based analyses, the alignment gaps were treated as ambiguous characters (setting: cleandata = 0). All models corrected the transition/transversion rate and codon usage biases (F3×4). Different starting ω values were also used to avoid local optima on the likelihood surface [60].

Three codon substitution models of maximum likelihood analysis were produced to detect selective pressure acting on the Hr and FGF5 genes: a site model, branch model, and branch-site model [61, 62]. The branch-specific models permitted variable ω ratios among branches but invariable ω ratios in the sites in the tree and they could be implemented to study changes in selective pressures in specific lineages [63]. In the branch-specific models, a ‘one-ratio’ model (M0) that assumed the same ω ratio for all branches [64] was compared with models where ω was allowed to differ in the background and a focal branch (two-ratio model). The branch-site models permitted the ω ratio to vary among sites and among lineages, which was useful for detecting positive selection that affected only a few sites in a few lineages [63, 65].

Significant for differences between two nested models were detected by using LRTs to calculate twice the log-likelihood (2ΔL) difference following a chi-square distribution, where the number of degrees of freedom was equal to the difference in the numbers of free parameters between models.

Declarations

Acknowledgments

Financial support was provided by the National Natural Science Foundation of China (NSFC) to GY (grant no. 30830016) and SXX (grant no. 31000953), the Priority Academic Program Development of Jiangsu Higher Education Institutions to GY and SXX, the Specialized Research Fund for the Doctoral Program of Higher Education (grant no. 20103207120010), the Ministry of Education of China to SXX, the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (grant no. 10KJB180002), and the Scientific Research Foundation of HNNU to ZC (grant no. 01046500138). We are grateful to Dr Anli Gao, Mr Xinrong Xu, and Dr Bingyao Chen for their contributions to sample collection.

Authors’ Affiliations

(1)
Jiangsu Key Laboratory for Biodiversity and Biotechnology, College of Life Sciences, Nanjing Normal University
(2)
College of Life Sciences, Henan Normal University

References

  1. Bergman J: Why mammal body hair is an evolutionary enigma. Creation Research Society Quarterly. 2004, 40: 240-243.Google Scholar
  2. Paus R, Cotsarelis G: The Biology of hair follicles. N Engl J Med. 1999, 341: 491-497. 10.1056/NEJM199908123410706.PubMedView ArticleGoogle Scholar
  3. Hardy MH: The secret life of the hair follicle. Trends Genet. 1992, 8 (2): 55-61. 10.1016/0168-9525(92)90350-D.PubMedView ArticleGoogle Scholar
  4. Cachon-Gonzalez MB, Fenner S, Coffin JM, Moran C, Best S, Stoye JP: Structure and expression of the hairless gene of mice. Proc Natl Acad Sci USA. 1994, 91 (16): 7717-7721. 10.1073/pnas.91.16.7717.PubMed CentralPubMedView ArticleGoogle Scholar
  5. Panteleyev AA, Botchkareva NV, Sundberg JP, Christiano AM, Paus R: The role of the hairless (hr) gene in the regulation of hair follicle catagen transformation. Am J Pathol. 1999, 155: 159-171. 10.1016/S0002-9440(10)65110-3.PubMed CentralPubMedView ArticleGoogle Scholar
  6. Hsieh JC, Sisk JM, Jurutka PW, Haussler CA, Slater SA, Haussler MR, Thompson CC: Physical and functional interaction between the vitamin D receptor and Hairless corepressor, two proteins required for hair cycling. J Biol Chem. 2003, 278: 38665-38674. 10.1074/jbc.M304886200.PubMedView ArticleGoogle Scholar
  7. Thompson CC: Hairless is a nuclear receptor corepressor essential for skin function. Nucl Recept Signal. 2009, 7: 1-11.Google Scholar
  8. Panteleyev AA, Botchkareva NV, van der Veen C, Christiano AM, Paus R: Pathobiology of the the hairless phenotype: dysregulation of hair follicle apoptosis and topobiology during the initiation of follicle cycling. J Invest Dermatol. 1998, 110: 577-Google Scholar
  9. Rosenquist TA, Martin GR: Fibroblast growth factor signalling in the hair growth cycle: expression of the fibroblast growth factor receptor and ligand genes in the murine hair follicle. Dev Dyn. 1996, 205: 379-86. 10.1002/(SICI)1097-0177(199604)205:4<379::AID-AJA2>3.0.CO;2-F.PubMedView ArticleGoogle Scholar
  10. 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.PubMedView ArticleGoogle Scholar
  11. 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.PubMedView ArticleGoogle Scholar
  12. 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.PubMed CentralPubMedView ArticleGoogle Scholar
  13. 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.PubMedView ArticleGoogle Scholar
  14. Jefferson TA, Webber MA, Pitman RL: Marine Mammals of the World: a Comprehensive Guide to their Identification. 2008, San Diego: Academic PressGoogle Scholar
  15. Gingerich P, Smith BH, Simons EL: Hind limbs of Eocene Basilosaurus: evidence of feet in whales. Science. 1990, 249: 154-7. 10.1126/science.249.4965.154.PubMedView ArticleGoogle Scholar
  16. Thewissen JG, Williams EM, Roe LJ, Hussain ST: Skeletons of terrestrial cetaceans and the relationship of whales to artiodactyls. Nature. 2001, 413: 277-281. 10.1038/35095005.PubMedView ArticleGoogle Scholar
  17. Thewissen JGM, Cooper LN, George JC, Bajpai S: From land to water: the origin of whales, dolphins and porpoises. Evol Edu Outreach. 2009, 2: 272-288. 10.1007/s12052-009-0135-2.View ArticleGoogle Scholar
  18. Xiong Y, Brandley MC, Xu S, Zhou K, Yang G: Seven new dolphin mitochondrial genomes and a time-calibrated phylogeny of whales. BMC Evol Biol. 2009, 9: 20-10.1186/1471-2148-9-20.PubMed CentralPubMedView ArticleGoogle Scholar
  19. Zhou X, Xu S, Yang Y, Zhou K, Yang G: Phylogenomic analyses and improved resolution of Cetartiodactyla. Mol Phylogenet Evol. 2011, 61: 255-264. 10.1016/j.ympev.2011.02.009.PubMedView ArticleGoogle Scholar
  20. McClellan DA, Palfreyman EJ, Smith MJ, Moss JL, Christensen RG, Sailsbery JK: Physicochemical evolution and molecular adaptation of the cetacean and artiodactyl cytochrome b proteins. Mol Biol Evol. 2005, 22: 437-455.PubMedView ArticleGoogle Scholar
  21. Wang Z, Yuan L, Rossiter SJ, Zuo X, Ru B, Zhong H, Han N, Jones G, Jepson PD, Zhang S: Adaptive evolution of 5’HoxD genes in the origin and diversification of the cetacean flipper. Mol Biol Evol. 2009, 26: 613-622.PubMedView ArticleGoogle Scholar
  22. Yu L, Jin W, Wang J, Zhang X, Chen M, Zhu Z, Lee H, Lee M, Zhang Y: Characterization of TRPC2, an essential genetic component of VNS chemoreception, provides insights into the evolution of pheromonal olfaction in secondary-adapted marine mammals. Mol Biol Evol. 2010, 27: 1467-1477. 10.1093/molbev/msq027.PubMedView ArticleGoogle Scholar
  23. Yu L, Jin W, Zhang X, Wang D, Zheng J, Yang G, Xu S, Cho S, Zhang Y: Evidence for positive selection on the leptin gene in Cetacea and Pinnipedia. PLoS One. 2011, 6 (10): e26579-10.1371/journal.pone.0026579.PubMed CentralPubMedView ArticleGoogle Scholar
  24. McGowen MR, Montgomery SH, Clark C, Gatesy J: Phylogeny and adaptive evolution of the brain-development gene microcephalin (MCPH1) in cetaceans. BMC Evol Biol. 2011, 11: 98-10.1186/1471-2148-11-98.PubMed CentralPubMedView ArticleGoogle Scholar
  25. Shen T, Xu S, Wang X, Yu W, Zhou K, Yang G: Adaptive evolution and functional constraint at TLR4 during the secondary aquatic adaptation and diversification of cetaceans. BMC Evol Biol. 2012, 12: 39-10.1186/1471-2148-12-39.PubMed CentralPubMedView ArticleGoogle Scholar
  26. Xu S, Chen Y, Cheng Y, Yang D, Zhou X, Xu J, Zhou K, Yang G: Positive selection at the ASPM gene coincides with brain size enlargements in cetaceans. Proc Roy Soc B-Biol Sci. 2012, 10.1098/rspb.2012.1729.Google Scholar
  27. Gatesy J, Milinkovitch MC, Waddelld V, Stanhope M: Stability of Cladistic Relationships between Cetacea and Higher Level Artiodactyl Taxa. Syst Biol. 1999, 48: 6-20. 10.1080/106351599260409.PubMedView ArticleGoogle Scholar
  28. Zhou X, Xu S, Xu J, Chen B, Zhou K, Yang G: Phylogenomic analysis resolves the interordinal relationships and rapid diversification of the Laurasiatherian mammals. Syt Biol. 2012, 61 (1): 150-164. 10.1093/sysbio/syr089.View ArticleGoogle Scholar
  29. Chen Z, Xu S, Zhou K, Yang G: Whale phylogeny and rapid radiation events revealed using novel retroposed elements and their flanking sequences. BMC Evol Biol. 2011, 11: 314-10.1186/1471-2148-11-314.PubMed CentralPubMedView ArticleGoogle Scholar
  30. Tsukada Y, Fang J, Erdjument-Bromage H, Warren ME, Borchers CH, Tempst P, Zhang Y: Histone demethylation by a family of JmjC domain-containing proteins. Nature. 2006, 439: 811-6. 10.1038/nature04433.PubMedView ArticleGoogle Scholar
  31. Hsieh JC, Slater SA, Whitfield GK, Dawson JL, Hsieh G, Sheedy C, Haussler CA, Haussler MR: Analysis of hairless corepressor mutants to characterize molecular cooperation with the vitamin D receptor in promoting the mammalian hair cycle. J Cell Biochem. 2010, 110: 671-686. 10.1002/jcb.22578.PubMed CentralPubMedView ArticleGoogle Scholar
  32. Moraitis AN, Giguere V, Thompson CC: Novel mechanism of nuclear receptor corepressor interaction dictated by activation function 2 helix determinants. Mol Cell Biol. 2002, 22: 6831-6841. 10.1128/MCB.22.19.6831-6841.2002.PubMed CentralPubMedView ArticleGoogle Scholar
  33. Potter GB, Zarach JM, Sisk JM, Thompson CC: The thyroid hormone-regulated corepressor hairless associates with histone deacetylases in neonatal rat brain. Mol Endocrinol. 2002, 16: 2547-60. 10.1210/me.2002-0115.PubMedView ArticleGoogle Scholar
  34. Ahmad W, Faiyaz ul Haque M, Brancolini V, Tsou HC, ul Haque S, Lam H, Aita VM, Owen J, deBlaquiere M, Frank J, Cserhalmi-Friedman PB, Leask A, McGrath JA, Peacocke M, Ahmad M, Ott J, Christiano AM: Alopecia universalis associated with a mutation in the human Hairless gene. Science. 1998, 279: 720-724. 10.1126/science.279.5351.720.PubMedView ArticleGoogle Scholar
  35. Cichon S, Anker M, Vogt IR, Rohleder H, Putzstuck M, Hillmer A, Farooq SA, Al-Dhafri KS, Ahmad M, Haque S, Rietschel M, Propping P, Kruse R, Nothen MM: Cloning, genomic organization, alternative transcripts and mutational analysis of the gene responsible for autosomal recessive universal congenital alopecia. Hum Mol Genet. 1998, 7: 1671-1679. 10.1093/hmg/7.11.1671.PubMedView ArticleGoogle Scholar
  36. Kim H, Wajid M, Kraemer L, Shimomura Y, Christiano AM: Nonsense mutations in the Hairless gene underlie APL in five families of Pakistani origin. J. Dermatol. Sci. 2007, 48: 207-211. 10.1016/j.jdermsci.2007.07.007.PubMed CentralPubMedView ArticleGoogle Scholar
  37. Klein I, Bergman R, Indelman M, Sprecher E: A novel missense mutation affecting the human hairless thyroid receptor interacting domain 2 causes congenital atrichia. J Invest Dermatol. 2002, 119: 920-922. 10.1046/j.1523-1747.2002.00268.x.PubMedView ArticleGoogle Scholar
  38. Kruse R, Cichon S, Anker M, Hillmer AM, Barros-Nunez P, Cantu JM, Leal E, Weinlich G, Schmuth M, Fritsch P, Ruzicka T, Propping P, Nothen MM: Novel hairless mutations in two kindreds with autosomal recessive popular atrichia. J Invest Dermatol. 1999, 113: 954-959. 10.1046/j.1523-1747.1999.00790.x.PubMedView ArticleGoogle Scholar
  39. Paradisi M, Chuang GS, Angelo C, Pedicelli C, Martinez-Mir A, Christiano AM: Atrichia with papular lesions resulting from a novel homozygous missense mutation in the hairless gene. Clin Exp Dermatol. 2003, 28: 535-538. 10.1046/j.1365-2230.2003.01333.x.PubMedView ArticleGoogle Scholar
  40. Zlotogorski A, Ahmad W, Christiano AM: Congenital atrichia in five Arab Palestinian families resulting from a deletion mutation in the human hairless gene. Hum Genet. 1998, 103: 400-404. 10.1007/s004390050840.PubMedView ArticleGoogle Scholar
  41. Zhang J, Fang S, Wang C: A Novel Nonsense Mutation and Polymorphisms in the Mouse Hairless Gene. J Invest Dermatol. 2005, 124: 1200-1205. 10.1111/j.0022-202X.2005.23744.x.PubMedView ArticleGoogle Scholar
  42. Panteleyev AA, Paus R, Ahmad W, Sundberg JP, Christiano AM: Molecular and functional aspects of the hairless (hr) gene in laboratory rodents and humans. Exp Dermatol. 1998, 7: 249-267.PubMedView ArticleGoogle Scholar
  43. Yampolsky LY, Stoltzfus A: The exchangeability of amino acids in proteins. Genetics. 2005, 170: 1459-1472. 10.1534/genetics.104.039107.PubMed CentralPubMedView ArticleGoogle Scholar
  44. Hecht HJ, Adar R, Hofmann B, Bogin O, Weich H, Yayon A: Structure of fibroblast growth factor 9 shows a symmetric dimmer with unique receptor- and heparin-binding interfaces. Acta Crystallogr D Biol Crystallogr. 2001, 57 (Pt3): 378-384.PubMedView ArticleGoogle Scholar
  45. Eswarakumar VP, Lax I, Schlessinger J: Cellular signaling by fibroblast growth factor receptors. Cytokine Growth Factor Rev. 2005, 16: 139-149. 10.1016/j.cytogfr.2005.01.001.PubMedView ArticleGoogle Scholar
  46. Li C, Jiang M, Chen S, Lai S: Correlation analysis between single nucleotide polymorphism of FGF5 gene and wool yield in rabbits. Yi Chuan. 2008, 30 (7): 893-899. 10.3724/SP.J.1005.2008.00893.PubMedView ArticleGoogle Scholar
  47. Chuong C-M: Molecular basis of epithelial appendage morphogenesis. Molecular biology intelligence unit 1. 1998, Austin, Tex: R.G. LandesGoogle Scholar
  48. Sutera R: The origin of whales and the power of independent evidence. Natl Cent Sci Educ. 2000, 20 (5): 33-41.Google Scholar
  49. Sambrook J, Russell DW: Molecular Cloning 4th edition. 2006, New York: Cold Spring Harbor Laboratory PressGoogle Scholar
  50. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 24: 4876-4882.View ArticleGoogle Scholar
  51. Yang Z, Kumar S, Nei M: A new method of inference of ancestral nucleotide and amino acid sequences. Genetics. 1995, 141: 1641-1650.PubMed CentralPubMedGoogle Scholar
  52. Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
  53. Arnold K, Bordoli L, Kopp J, Schwede T: The SWISS-MODEL Workspace: A web-based environment for protein structure homology modelling. Bioinformatics. 2006, 22: 195-201. 10.1093/bioinformatics/bti770.PubMedView ArticleGoogle Scholar
  54. Schwede T, Kopp J, Guex N, Peitsch MC: SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 2003, 31: 3381-3385. 10.1093/nar/gkg520.PubMed CentralPubMedView ArticleGoogle Scholar
  55. Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: An environment for comparative protein modelling. Electrophoresis. 1997, 18: 2714-2723. 10.1002/elps.1150181505.PubMedView ArticleGoogle Scholar
  56. Helaers R, Milinkovitch MC: MetaPIGA v2.0: maximum likelihood large phylogeny estimation using the metapopulation genetic algorithm and other stochastic heuristics. BMC Bioinformatics. 2010, 11: 379-10.1186/1471-2105-11-379.PubMed CentralPubMedView ArticleGoogle Scholar
  57. Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.PubMedView ArticleGoogle Scholar
  58. Nylander JAA: MrModeltest v2. Program distributed by the author. 2004, Sweden: Evolutionary Biology Centre, Uppsala University, Available from: http://www.abc.se/~nylander/mrmodeltest2/mrmodeltest2.htmlGoogle Scholar
  59. Rambaut A, Drummond AJ: Tracer v1.4. Distributed by the Authors. Available from: http://beast.bio.ed.ac.uk/Trace, 2007
  60. Suzuki Y, Nei M: Reliabilities of parsimony-based and likelihood-based methods for detecting positive selection at single amino acid sites. Mol Biol Evol. 2001, 18: 2179-2185. 10.1093/oxfordjournals.molbev.a003764.PubMedView ArticleGoogle Scholar
  61. Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics. 1998, 148: 929-936.PubMed CentralPubMedGoogle Scholar
  62. Yang Z, Wong WSW, Nielsen R: Bayes empirical Bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22: 1107-1118. 10.1093/molbev/msi097.PubMedView ArticleGoogle Scholar
  63. Yang Z, Nielsen R: Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. Mol Biol Evol. 2002, 19: 908-917. 10.1093/oxfordjournals.molbev.a004148.PubMedView ArticleGoogle Scholar
  64. Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15: 568-573. 10.1093/oxfordjournals.molbev.a025957.PubMedView ArticleGoogle Scholar
  65. 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: 2472-2479. 10.1093/molbev/msi237.PubMedView ArticleGoogle Scholar

Copyright

© Chen et al.; licensee BioMed Central Ltd. 2013

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