- Research article
- Open Access
Evolution of the hepcidin gene in primates
BMC Genomics volume 9, Article number: 120 (2008)
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.
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.
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.
Antimicrobial peptides (AMPs) are phylogenetically ancient and ubiquitous molecules that play an essential role in many aspects of host defence . 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 . 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 .
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 Liver Expressed) and its antimicrobial activity in vitro (-cidin, or Antimicrobial Peptide) [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 . 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 . 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.
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.
At the amino acid level, identity ranges from 98.8% in chimpanzee (PTR) to 83.3% in NWM Ateles (AFU) (see Table 1).
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  (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.
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 . 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  and residues Phe4, Pro5, Phe9, and Met21 that are involved in the oligomerization of hepcidin-25  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) , which may be functionally significant. While the conserved N-terminal region is recognised to be essential for interaction with ferroportin , the role of the C-terminal side is not well defined.
Shi et al  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.
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 . In this context, while some types of residue variations, such as those observed in mouse hepcidin 2 , 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.
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.
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).
Sequence alignment and Data Analysis
Multiple alignment of the nucleotide and amino acid sequences were performed using the program CLUSTAL W .
A phylogenetic tree was obtained by the UPGMA method, using the PHYLIP software package . Maximum Likelihood analysis was performed with the PAML software package  by comparing the log-likelihood ratios of the data using different NS-sites models . 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) .
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.
Zasloff M: Antimicrobial peptides of multicellular organisms. Nature. 415 (6870): 389-95. 10.1038/415389a. 2002 Jan 24
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. 10.1007/PL00000914.
Tossi A, Sandri L: Molecular diversity in gene-encoded, cationic antimicrobial polypeptides. Curr Pharm Des. 2002, 8 (9): 743-761. 10.2174/1381612023395475.
Park CH, Valore EV, Waring AJ, Ganz T: Hepcidin, a urinary antimicrobial peptide synthesized in the liver. J Biol Chem. 276 (11): 7806-7810. 10.1074/jbc.M008922200. 2001 Mar 16
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
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. 10.1073/pnas.151179498. 2001 Jul 17
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. 10.1152/ajpgi.00412.2005.
Ganz T, Nemeth E: Regulation of iron acquisition and iron distribution in mammals. Biochim Biophys Acta. 2006, 1763 (7): 690-9. 10.1016/j.bbamcr.2006.03.014.
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. 10.1074/jbc.M008923200. 2001 Mar 16
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.
Verga Falzacappa MV, Muckenthaler MU: Hepcidin: iron-hormone and anti-microbial peptide. Gene. 364: 37-44. 10.1016/j.gene.2005.07.020. 2005 Dec 30
Ganz T: Molecular control of iron transport. J Am Soc Nephrol. 2007, 18 (2): 394-400. 10.1681/ASN.2006070802.
Shi J, Camus AC: Hepcidins in amphibians and fishes: Antimicrobial peptides or iron-regulatory hormones?. Dev Comp Immunol. 2006, 30 (9): 746-55. 10.1016/j.dci.2005.10.009.
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. 10.1007/s11010-007-9543-1.
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. 10.1074/jbc.M511108200. 2006 Jul 21
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. 10.2174/1389203053027593.
Fry MM, Liggett JL, Baek SJ: Molecular cloning and expression of canine hepcidin. Vet Clin Pathol. 2004, 33 (4): 223-7. 10.1111/j.1939-165X.2004.tb00377.x.
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. 10.1182/blood-2005-05-2049. 2006 Jan 1
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. 10.1038/ng1053.
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. 10.1093/hmg/ddg225. 2003 Sep 1
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. 10.1182/blood-2004-01-0332. 2004 Oct 1
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. 10.1074/jbc.M205305200. 2002 Oct 4
Semple CA, Rolfe M, Dorin JR: Duplication and selection in the evolution of primate beta-defensin genes. Genome Biol. 2003, 4 (5): R31-10.1186/gb-2003-4-5-r31.
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
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. 10.1007/s00251-001-0412-x.
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
Aiyar A: The use of CLUSTAL W and CLUSTAL X for multiple sequence alignment. Methods Mol Biol. 2000, 132: 221-41.
Phylip software. [http://evolution.genetics.washington.edu/phylip.html]
Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13 (5): 555-6.
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.
Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15 (5): 568-73.
This study was supported by a grant 200502027001 from "Regione Friuli Venezia Giulia" and by a grant from IRCCS Burlo Garofolo (RC 2007).
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.
About this article
Cite this article
Segat, L., Pontillo, A., Milanese, M. et al. Evolution of the hepcidin gene in primates. BMC Genomics 9, 120 (2008). https://doi.org/10.1186/1471-2164-9-120
- Mature Peptide
- Porphyria Cutanea Tarda
- Host Defense Peptide
- Hepcidin Gene
- Identify Amino Acid Residue