Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

Evolution of the hepcidin gene in primates

Ludovica Segat1*, Alessandra Pontillo1, Michele Milanese1, Alessandro Tossi2 and Sergio Crovella13

Author Affiliations

1 Genetic Unit, IRCCS Burlo Garofolo and Department of Reproductive and Developmental Sciences, University of Trieste, Trieste, Italy

2 Department of Biochemistry, Biophysics and Macromolecular Chemistry, University of Trieste, Trieste, Italy

3 Department of Genetics, Federal University of Pernambuco, Recife, Brazil

For all author emails, please log on.

BMC Genomics 2008, 9:120  doi:10.1186/1471-2164-9-120


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/9/120


Received:28 August 2007
Accepted:5 March 2008
Published:5 March 2008

© 2008 Segat et al; licensee 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.

Abstract

Background

Hepcidin/LEAP-1 is an iron regulatory hormone originally identified as an antimicrobial peptide. As part of a systematic analysis of the evolution of host defense peptides in primates, we have sequenced the orthologous gene from 14 species of non-human primates.

Results

The sequence of the mature peptide is highly conserved amongst all the analyzed species, being identical to the human one in great apes and gibbons, with a single residue conservative variation in Old-World monkeys and with few substitutions in New-World monkeys.

Conclusion

Our analysis indicates that hepcidin's role as a regulatory hormone, which involves interaction with a conserved receptor (ferroportin), may result in conservation over most of its sequence, with the exception of the stretch between residues 15 and 18, which in New-World monkeys (as well as in other mammals) shows a significant variation, possibly indicating that this structural region is involved in other functions.

Background

Antimicrobial peptides (AMPs) are phylogenetically ancient and ubiquitous molecules that play an essential role in many aspects of host defence [1]. In plants and invertebrate animals they are central to protecting the host from infection, while in vertebrate animals that have developed an acquired immunity they not only provide a rapidly mobilised first line of defence as part of the innate immune system, but also act as signalling molecules to cellular components of both innate and adaptive immunity [2]. Being at the interface between host and commensal or pathogenic microbes, they generally undergo an accelerated evolution, which has resulted in many different AMP families and a remarkable sequence variation even within structural groups [3].

In some cases, molecules derived from AMPs have acquired functions beyond host defence, although they can maintain a role in the prevention of infection. The recently identified peptide hepcidin/LEAP-1 is a good example. It was independently isolated from human urine and plasma ultra filtrate, respectively, by two groups searching for novel AMPs. Both denominations thus reflect its hepatic expression (

    hep
- or
    L
iver
    E
xpressed) and its antimicrobial activity in vitro (-
    cidin
, or
    A
ntimicrobial
    P
eptide) [4,5]. Subsequently, the observation that disruption or mutation of the hepcidin gene in mouse and man results in severe forms of iron overload identified it also as a negative regulator of iron absorption, recycling and release from stores [6-8]. It functions by inducing the internalisation of the ubiquitous, vertebrate, iron export membrane protein ferroportin, at the major sites of iron flow (enterocytes, tissue macrophages and hepatocytes), causing a decrease in serum iron levels. It was thus proposed that hepcidin was the long sought iron regulatory hormone, whose synthesis is controlled by iron levels and erythroid demands [9,10]. Hepcidin maintains a role in host defence by acting as a bridge between immunity and iron-metabolism, as indicated by the fact that it is markedly induced by infection and inflammation [7,8,11]. This likely reflects an antimicrobial action via iron limitation rather than a direct antimicrobial activity [12,13].

In humans, hepcidin is encoded by three small exons in the HAMP gene (19q13.1), and translates to an 84aa pre-propeptide presenting a furin cleavage site. It is expressed and processed predominantly in liver for plasma delivery, and it is excreted by the kidney. In urine, the predominant form is of 25 amino acids (hep-25) with a charge of +2, but two shorter, inactive peptides (hep-22 and hep-20) were also detected [4]. It does not show significant sequence similarity to any other known AMP, but it has been associated with the defensin family of host defence peptides due to its salt-sensitive antimicrobial activity in vitro, cationic and amphipathic nature, and the presence of 4 disulfide bridges in its tertiary structure. However, unlike defensins, it assumes a hairpin conformation which is stabilised by four S-S bridges placed in a ladder-like manner, including an unusual and possibly reactive vicinal disulphide. Structurally, it is therefore more reminiscent of AMPs such as protegrin [8]. It is not clear whether the observed in vivo antimicrobial activity is biologically relevant, is simply a residue of a possible evolutionary origin as an AMP, or derives from the conserved structural characteristics imposed by its function as an iron-regulating hormone.

Hepcidin homologues and putative hepcidin precursors have been identified in several vertebrate animals, ranging from mammals and birds to amphibians and fish [11,13,14]. Homology amongst the mature peptide sequences is quite high, reflecting a high degree of structural conservation and pointing to a role as iron regulatory hormones also in non-mammalian animals. Some species however express more than one hepdicin gene, whose products may have functions other than iron regulation, and in fish they may indeed still act as AMPs [11,13]. It has in fact been hypothesised that the iron-regulatory hormones have evolved from ancient fish AMPs.

We have sequenced the hepcidin gene coding region (all three exons) in 14 non-human primate species and compared it with the human one (HSS), so as to gain information about its evolution in primates and determine whether it could usefully contribute to their phylogenetic analysis. This study was carried out in the context of a broader analysis of primate host defense peptide evolution, which has included the cathelicidin LL-37 and beta-defensins, and which has shown remarkably varied patterns of molecular evolution, ranging from conservation to positive selection for variation [15,16]. By investigating the degree of variation of the hepcidins in the closely related primate species we aimed to determine the extent of purifying selection during the evolution of this gene in primates and locate regions of variation that could provide useful information on the functional determinants of this multifunctional peptide hormone.

Results

The HAMP gene nucleotide sequences obtained from the analyzed primate species, and the deduced corresponding amino acid sequences, are reported in Figures 1 and 2. They clearly demonstrate a high identity to the human sequence at both the nucleotide and amino acid levels. The nucleotide identity for the entire coding and untranslated (UTR) regions ranges from 99.7% with chimpanzee (PTR) to 92.6% with the NWM Ateles fusciceps (AFU). An in-frame insertion of three nucleotides (namely GGA, leading to the insertion of a glycine between residue 38 and 39 of the pre-propeptide) occurs in Hylobates.

thumbnailFigure 1. Alignment of nucleotide sequences from coding and UTR regions of orthologous primate hepcidin genes. The UTR regions are underlined. Dots indicate nucleotide identity with the human sequence (GeneBank accession number: NM_021175). Non-synonymous substitutions are indicated in bold.

thumbnailFigure 2. Alignment of primate hepcidin deduced amino acid sequences. Dots indicate residue identity with the human sequence (GeneBank accession number: NP_066998). The signal peptide region is underlined and the mature peptide sequence is highlighted in grey.

At the amino acid level, identity ranges from 98.8% in chimpanzee (PTR) to 83.3% in NWM Ateles (AFU) (see Table 1).

Table 1. Percent identity of nucleotide (nt) and amino acid (aa) sequences calculated versus the human sequence. The highest and lowest values per each kind of identity (nt and aa) are underlined.

The HAMP gene coding sequences were then used to compute an UPGMA tree with all primate sequences, using the Canis lupus familiaris HAMP gene as the outgroup [17] (Figure 3), which reflects the expected phylogenetic relationship both among families and between families. A maximum-likelihood approach, using the PAML suite of programs, was then used to determine if positive selection acted on the primate HAMP genes, and globally found that it has not (ω<1). However, a site-analysis found that individual amino acid residues may be under positive selection (ω>1), including the Gly insertion between residue 38 and 39, and Phe39Gln and Gln40Leu substitutions in the Hylobates sequences (see Figure 2), and the His74Arg, Arg75Gln and Ser76Pro substitution specific to the NWM sequences. A subsequent branch-analysis also indicated that although the mean ω is <1 over the entire tree, positive selection may act on the Hylobates and NWM branches, confirming the site-analysis data.

thumbnailFigure 3. Rooted phylogenetic tree based on primate hepcidin nucleotide sequences, prepared using the UPGMA method and using the Canis lupus familiaris sequence as outgroup (GeneBank accession number NM_001007140)

Discussion

To assess whether the variations found at the amino acid level in the different primates species could have functional implications, these were analyzed in relation to the available data for the structure and activity of the human hepcidin, dividing the amino acid sequence into three parts, namely signal peptide, pro-region and mature peptide region (see Figure 2).

The entire sequence is quite conserved, but more so at the level of the mature peptide. In effect, the only significant variation in the pro-region corresponds to the insertion observed for the Hylobates species, whose effect is unpredictable as it falls outside the mature peptide. The putative mature peptide sequence is identical to the human one in great apes and gibbons, and shows a single conservative Lys83Arg variation in Old-World monkeys. Interestingly, a heterozygous Lys83Arg replacement leading to a peptide identical to that of the OWMs has been found in a patient with porphyria cutanea tarda and iron overload, but testing of the synthetic homologue revealed it to remain fully functional in iron regulation, so the effective role of this substitution is not yet clear [18]. Other mutations in the hepcidin gene previously reported to be associated with hemochromatosis in humans [19-21] were not found in the primate species examined (at least 5 individuals per species have been studied).

The pro-peptide convertase cleavage site, the functionally essential N-terminal region [18] and residues Phe4, Pro5, Phe9, and Met21 that are involved in the oligomerization of hepcidin-25 [22] are all conserved, as well as all cysteine residues that account for the formation of disulphide bond and confer hepcidin's characteristic structure. Only the New-World monkeys show significant variation in a stretch between residues 15 and 18 of the mature peptide sequence (74His, 75Arg and 76Ser, using the full pro-peptide sequence numbering) (see figure 3) which may be under positive selection. Alignment of the primate sequences with those of other mammals, as well as bird, amphibian and fish sequences [11,13,14], indicate that this region is in fact the most variable also with respect to other mammalian sequences (see figure 4a). This stretch is just C-terminal to the two Cys residues forming the vicinal disulphide at the peptide's hairpin turn (see figure 4b) [22], which may be functionally significant. While the conserved N-terminal region is recognised to be essential for interaction with ferroportin [18], the role of the C-terminal side is not well defined.

thumbnailFigure 4. A) Alignment of primate hepcidins with those of other mammals, amphibians and fish. Residues identical to the human sequence are highlighted in grey. Part of the pro-peptide corresponding to the consensus propeptide convertase (furin) cleavage site is shown in italics. Arrows indicate the residues that are involved in the oligomerization of human hepcidin-25. (GeneBank accession numbers for non-primate species: Pig NM_214117; Mouse-1 BC_02158; Mouse-2 AY_232841; Rat NM_053469; Pigeon see ref. [14]; Frog DN_020182; Zebra fish NM_205583; Salmon BQ_036900; Bass DQ_31605). B) Structure of human hepcidin-25, based on the coordinate file 1M4F.pdb. The structure schematically shows the ladder-like disulphide connectivity (numbered Cys residues are indicated) and location of the variable stretch (residues 15–18) in light grey.

Shi et al [13] have hypothesized that, during vertebrate evolution, hepcidin differentiated its functions towards an iron-regulatory hormone starting from a more ancient peptide characterized by antimicrobial activity. In this respect, it is interesting to note that substitutions in the variable region, concerning both primates and other animals, often concerns charged residues (see Figure 4); moreover, it has been shown for different types of AMPs that variation in charge are also favoured by positive selection [15,23]. Whereas for the latter peptides, charge modulation can be directly related to the efficiency of interaction with microbes, thus affecting antimicrobial activity, in the case of hepcidin, the effect of this variation on function/s is an open question.

Conclusion

Our results indicate that hepcidin is quite conserved in all primate species. This may derive from the fact that its role as an iron regulating hormone requires interaction with the conserved vertebrate iron transporter ferroportin [24]. In this context, while some types of residue variations, such as those observed in mouse hepcidin 2 [13], can result in a loss of the iron regulatory role, other variations, such as those observed in the stretch between residues 15 and 18 in NWM and other mammals, can more freely occur without apparently disrupting the hormonal function. It remains to be established whether these residues are involved in other possible functions of the hepcidin peptides, consistent with the possibility that they evolved under positive selection.

Methods

Samples

We have analyzed the HAMP gene in three species of great apes [Pan troglodytes (PTR), Gorilla gorilla (GGO), Pongo pygmaeus (PPY)], two species of Hylobatidae [Hylobates lar (HLA) and Hylobates concolor (HCO), now classified as Nomascus concolor], seven species of Cercopithecidae [Old World Monkeys (OWM) ; Presbytis obscurus (POB), Presbytis cristata (PCR), Presbytis melalophos (PME), Papio papio (PPA), Macaca fascicularis (MFA), Macaca fuscata (MFU), Cercopithecus aethiops (CAE)], and two Platyrrhine species [New World monkeys (NWM); Callithrix jacchus (CJA), Ateles fusciceps (AFU)].

We tried to minimize this risk of nucleotide intra-specific variations, by analyzing more than one individual for each species and determining a consensus sequence. This method allowed us to clear our data for single nucleotide polymorphisms variations in most species.

DNA extraction

Genomic DNA was extracted from hairs, liver, and muscle tissues of primates as described previously by Del Pero et al. [25] and Boniotto et al. [26], using the DNeasy Tissue kit (Qiagen).

PCR and sequencing

Genomic DNA amplification of the hepcidin three exons was performed using primers (see Table 2), designed on the basis of the human sequence available in GenBank (NM_021175) using the Primer Express 2.0 software (Applied Biosystems, Foster City, CA). The PCR reactions were carried out in a Thermal cycler 9700 (Applied Biosystems) using PCR Buffer 1×, 1 unit of Taq Gold, 0.4 mM dNTPs and variant concentration of MgCl2 (from 1.5 to 3 mM). The amplification conditions were 30 sec at 95°C, 30 sec at variable annealing temperature between the different templates (from 52°C up to 60°C), 30 sec at 72°C for 35 cycles. PCR products were observed, under UV light, in a 2% agarose gel, stained with ethidium bromide. DNA sequencing of PCR products was performed using the BigDye Terminator Cycle Sequencing Ready Reaction Kit v. 2.0 (Applied Biosystems). DNA sequences were detected and analyzed on an automated ABI Prism 3100 Genetic Analyser (Applied Biosystems).

Table 2. Primers used to perform the genomic DNA amplification. Additional couples of primers (e.g. 1N, 1NN and 2-3N) were used for each amplicon as it was not possible to amplify DNA from all the primates using a single set (e.g. 1 and 2–3).

Sequence alignment and Data Analysis

Multiple alignment of the nucleotide and amino acid sequences were performed using the program CLUSTAL W [27].

A phylogenetic tree was obtained by the UPGMA method, using the PHYLIP software package [28]. Maximum Likelihood analysis was performed with the PAML software package [29] by comparing the log-likelihood ratios of the data using different NS-sites models [30]. This type of analysis is able to identify amino acid residues with high posterior probabilities (greater than 0.95) of having evolved under positive selection. ω ratios for the tree were calculated using different models (free ratio, two- and three ratio) [31].

GeneBank accession numbers for non-human primate species

Pan troglodytes (PTR): EU076436; Gorilla gorilla (GGO): EU076444; Pongo pygmaeus (PPY): EU076437; Hylobates concolor (HCO): EU076448; Hylobates lar (HLA): EU076449; Presbytis obscurus (POB): EU076439; Presbytis cristata (PCR): EU076441; Presbytis melalophos (PME): EU076440; Papio papio (PPA): EU076438; Macaca fascicularis (MFA): EU076443; Macaca fuscata (MFU): EU076442; Cercopithecus aethiops (CAE): EU076446; Callithrix jacchus (CJA): EU076445; Ateles fusciceps (AFU): EU076447.

Authors' contributions

LS performed the molecular genetic studies and participated in data interpretation; MM carried out the molecular genetic studies; AP performed the sequence alignment and evolutionary analysis; AT participated in the study design and sequence analyses and wrote the manuscript. SC conceived the study, and participated in its design and coordination and revised the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This study was supported by a grant 200502027001 from "Regione Friuli Venezia Giulia" and by a grant from IRCCS Burlo Garofolo (RC 2007).

References

  1. Zasloff M: Antimicrobial peptides of multicellular organisms.

    Nature 415(6870):389-95.

    2002 Jan 24

    PubMed Abstract | Publisher Full Text OpenURL

  2. Yang D, Chertov O, Oppenheim JJ: The role of mammalian antimicrobial peptides and proteins in awakening of innate host defenses and adaptive immunity.

    Cell Mol Life Sci 2001, 58(7):978-89. PubMed Abstract | Publisher Full Text OpenURL

  3. Tossi A, Sandri L: Molecular diversity in gene-encoded, cationic antimicrobial polypeptides.

    Curr Pharm Des 2002, 8(9):743-761. PubMed Abstract | Publisher Full Text OpenURL

  4. Park CH, Valore EV, Waring AJ, Ganz T: Hepcidin, a urinary antimicrobial peptide synthesized in the liver.

    J Biol Chem 276(11):7806-7810.

    2001 Mar 16

    PubMed Abstract | Publisher Full Text OpenURL

  5. Krause A, Neitz S, Mägert HJ, Schulz A, Forssmann WG, Schulz-Knappe P, Adermann K: LEAP-1, a novel highly disulfide-bonded human peptide, exhibits antimicrobial activity.

    FEBS Lett 480(2–3):147-50.

    2000 Sep 1

    PubMed Abstract | Publisher Full Text OpenURL

  6. Nicolas G, Bennoun M, Devaux I, Beaumont C, Grandchamp B, Kahn A, Vaulont S: Lack of hepcidin gene expression and severe tissue iron overload in upstream stimulatory factor 2 (USF2) knockout mice.

    Proc Natl Acad Sci USA 98(15):8780-5.

    2001 Jul 17

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Ganz T, Nemeth E: Iron imports: IV. Hepcidin and regulation of body iron metabolism.

    Am J Physiol Gastrointest Liver Physiol 2006, 290(2):G199-203. PubMed Abstract | Publisher Full Text OpenURL

  8. Ganz T, Nemeth E: Regulation of iron acquisition and iron distribution in mammals.

    Biochim Biophys Acta 2006, 1763(7):690-9. PubMed Abstract | Publisher Full Text OpenURL

  9. Pigeon C, Ilyin G, Courselaud B, Leroyer P, Turlin B, Brissot P, Loréal O: A new mouse liver specific gene, encoding a protein homologous to human antimicrobial peptide hepcidin, is overexpressed during iron overload.

    J Biol Chem 276(11):7811-9.

    2001 Mar 16

    PubMed Abstract | Publisher Full Text OpenURL

  10. Nicolas G, Chauvet C, Viatte L, Danan JL, Bigard X, Devaux I, Beaumont C, Kahn A, Vaulont S: The gene encoding the iron regulatory peptide hepcidin is regulated by anemia, hypoxia, and inflammation.

    J Clin Invest 2002, 110(7):1037-44. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Verga Falzacappa MV, Muckenthaler MU: Hepcidin: iron-hormone and anti-microbial peptide.

    Gene 364:37-44.

    2005 Dec 30

    PubMed Abstract | Publisher Full Text OpenURL

  12. Ganz T: Molecular control of iron transport.

    J Am Soc Nephrol 2007, 18(2):394-400. PubMed Abstract | Publisher Full Text OpenURL

  13. Shi J, Camus AC: Hepcidins in amphibians and fishes: Antimicrobial peptides or iron-regulatory hormones?

    Dev Comp Immunol 2006, 30(9):746-55. PubMed Abstract | Publisher Full Text OpenURL

  14. Fu YM, Li SP, Wu YF, Chang YZ: Identification and expression analysis of hepcidin-like cDNAs from pigeon (Columba livia).

    Mol Cell Biochem 2007, 305(1–2):191-7. PubMed Abstract | Publisher Full Text OpenURL

  15. Zelezetsky I, Pontillo A, Puzzi L, Antcheva N, Segat L, Pacor S, Crovella S, Tossi A: Evolution of the primate cathelicidin. Correlation between structural variations and antimicrobial activity.

    J Biol Chem 281(29):19861-71.

    2006 Jul 21

    PubMed Abstract | Publisher Full Text OpenURL

  16. Crovella S, Antcheva N, Zelezetsky I, Boniotto M, Pacor S, Verga Falzacappa MV, Tossi A: Primate beta-defensins – structure, function and evolution.

    Curr Protein Pept Sci 2005, 6(1):7-21. PubMed Abstract | Publisher Full Text OpenURL

  17. Fry MM, Liggett JL, Baek SJ: Molecular cloning and expression of canine hepcidin.

    Vet Clin Pathol 2004, 33(4):223-7. PubMed Abstract | Publisher Full Text OpenURL

  18. Nemeth E, Preza GC, Jung CL, Kaplan J, Waring AJ, Ganz T: The N-terminus of hepcidin is essential for its interaction with ferroportin: structure-function study.

    Blood 107(1):328-33.

    2006 Jan 1

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Roetto A, Papanikolaou G, Politou M, Alberti F, Girelli D, Christakis J, Loukopoulos D, Camaschella C: Mutant antimicrobial peptide hepcidin is associated with severe juvenile hemochromatosis.

    Nature Genet 2003, 33(1):21-2. PubMed Abstract | Publisher Full Text OpenURL

  20. Merryweather-Clarke AT, Cadet E, Bomford A, Capron D, Viprakasi V, Miller A, McHugh PJ, Chapman RW, Pointon JJ, Wimhurst VLC, Livesey KJ, Tanphaichitr V, Rochette J, Robson KJH: Digenic inheritance of mutations in HAMP and HFE results in different types of haemochromatosis.

    Hum Molec Genet 12(17):2241-7.

    2003 Sep 1

    PubMed Abstract | Publisher Full Text OpenURL

  21. Matthes T, Aguilar-Martinez P, Pizzi-Bosman L, Darbellay R, Rubbia-Brandt L, Giostra E, Michel M, Ganz T, Beris P: Severe hemochromatosis in a Portuguese family associated with a new mutation in the 5-prime-UTR of the HAMP gene.

    Blood 104(7):2181-3.

    2004 Oct 1

    PubMed Abstract | Publisher Full Text OpenURL

  22. Hunter HN, Fulton DB, Ganz T, Vogel HJ: The solution structure of human hepcidin, a peptide hormone with antimicrobial activity that is involved in iron uptake and hereditary hemochromatosis.

    J Biol Chem 277(40):37597-603.

    2002 Oct 4

    PubMed Abstract | Publisher Full Text OpenURL

  23. Semple CA, Rolfe M, Dorin JR: Duplication and selection in the evolution of primate beta-defensin genes.

    Genome Biol 2003, 4(5):R31. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  24. Donovan A, Brownlie A, Zhou Y, Shepard J, Pratt SJ, Moynihan J, Paw BH, Drejer A, Barut B, Zapata A, Law TC, Brugnara C, Lux SE, Pinkus GS, Pinkus JL, Kingsley PD, Palis J, Fleming MD, Andrews NC, Zon LI: Positional cloning of zebrafish ferroportin1 identifies a conserved vertebrate iron exporter.

    Nature 403(6771):776-81.

    2000 Feb 17

    PubMed Abstract | Publisher Full Text OpenURL

  25. Del Pero M, Boniotto M, Zuccon D, Cervella P, Spano A, Amoroso A, Crovella S: Beta-defensin 1 gene variability among non-human primates.

    Immunogenetics 2002, 53(10–11):907-13. PubMed Abstract | Publisher Full Text OpenURL

  26. Boniotto M, Antcheva N, Zelezetsky I, Tossi A, Palumbo V, Verga Falzacappa MV, Sgubin S, Braida L, Amoroso A, Crovella S: A study of host defence peptide beta-defensin 3 in primates.

    Biochem J 374(Pt 3):707-14.

    2003 Sep 15

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Aiyar A: The use of CLUSTAL W and CLUSTAL X for multiple sequence alignment.

    Methods Mol Biol 2000, 132:221-41. PubMed Abstract OpenURL

  28. Phylip software [http://evolution.genetics.washington.edu/phylip.html] webcite

  29. Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood.

    Comput Appl Biosci 1997, 13(5):555-6. PubMed Abstract OpenURL

  30. Yang Z, Nielsen R, Goldman N, Pedersen AM: Codon-substitution models for heterogeneous selection pressure at amino acid sites.

    Genetics 2000, 155(1):431-49. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution.

    Mol Biol Evol 1998, 15(5):568-73. PubMed Abstract | Publisher Full Text OpenURL