- Research article
- Open Access
Genetic study of the hepcidin gene (HAMP) promoter and functional analysis of the c.-582A > G variant
© Parajes et al; licensee BioMed Central Ltd. 2010
- Received: 12 July 2010
- Accepted: 10 December 2010
- Published: 10 December 2010
Hepcidin acts as the main regulator of iron homeostasis through regulation of intestinal absorption and macrophage release. Hepcidin deficiency causes iron overload whereas its overproduction is associated with anaemia of chronic diseases. The aims of the study were: to identify genetic variants in the hepcidin gene (HAMP) promoter, to asses the associations between the variants found and iron status parameters, and to functionally study the role on HAMP expression of the most frequent variant.
The sequencing of HAMP promoter from 103 healthy individuals revealed two genetic variants: The c.-153C > T with a frequency of 0.014 for allele T, which is known to reduce hepcidin expression and the c.-582A > G with a 0.218 frequency for allele G. In an additional group of 224 individuals, the c.-582A > G variant genotype showed no association with serum iron, transferrin or ferritin levels.
The c.-582G HAMP promoter variant decreased the transcriptional activity by 20% compared to c.-582A variant in cells from the human hepatoma cell line HepG2 when cotransfected with luciferase reporter constructs and plasmid expressing upstream stimulatory factor 1 (USF1) and by 12-14% when cotransfected with plasmid expressing upstream stimulatory factor 2 (USF2).
The c.-582A > G HAMP promoter variant is not associated with serum iron, transferrin or ferritin levels in the healthy population. The in vitro effect of the c.-582A > G variant resulted in a small reduction of the gene transactivation by allele G compared to allele A. Therefore the effect of the variant on the hepcidin levels in vivo would be likely negligible. Finally, the c.-153C > T variant showed a frequency high enough to be considered when a genetic analysis is done in iron overload patients.
- Iron Overload
- Serum Iron
- Serum Ferritin Level
- Transferrin Saturation
- Serum Transferrin
Hepcidin was first isolated from urine as an antimicrobial peptide  and thereafter, it has been shown to act as the main regulator of iron homeostasis through regulation of intestinal iron absorption and macrophage iron release [2–5]. The first evidence was provided by mice studies which showed increased hepcidin mRNA levels in the liver of iron overload mice . Further animal models and human disorders confirmed the relationship between this peptide and iron. Genetically modified mice showed that animals with reduced hepcidin expression developed severe iron overload , while those with increased expression had severe iron deficiency anaemia at birth . In humans, inactivating mutations of hepcidin result in a rare form of juvenile haemochromatosis , whereas hepcidin overexpression in inflammation causes anaemia of chronic diseases with features of iron restricted erythropoiesis . Besides, hepcidin has been described as a modifier peptide that exacerbates the phenotypic expression of patients with hereditary haemochromatosis (HFE) homozygous for the common HFE gene mutation, p.Cys282Tyr [11–13].
The gene encoding hepcidin, HAMP, is mainly expressed in the liver in situations of iron overload . It is upregulated by hemojuvelin , transferrin receptor 2  and HFE protein  mediated through the bone morphogenetic protein (BMP) signalling pathway . Hepcidin is also upregulated in inflammation by interleukin 6 (IL-6) and other cytokines through the Janus kinase/signal transducer and activator of transcription (JAK/STAT) signalling pathway [6, 10, 19]. Hepcidin is inhibited by hypoxia  and augmented erythropoiesis . More recently, matriptase-2 has been described as the first known negative regulator of hepcidin expression .
Lately, two genetic variants in the HAMP promoter have been described as modulators of iron overload. Firstly, the c.-582A > G variant has been associated with higher liver iron concentration and with higher serum ferritin levels in beta-thalassemic patients . Secondly, the c.-153C > T mutation, which is located within a BMP-Responsive Element, has been shown to contribute to a severe phenotype in HFE related haemochromatosis .
The aim of the study was to identify genetic variants in the HAMP promoter and their frequency distribution in the Galician population. Another objective was to assess possible associations between these variants and the serum iron, serum transferrin and serum ferritin levels in a random sample of Galician probands. Finally, we performed functional in vitro studies to determine the effect of variants on HAMP expression.
Identification of genetic variants in the hepcidin gene promoter
The sequencing of the HAMP promoter in the random sample of 103 healthy individuals allowed the identification of two genetic variants each with an allele frequency higher than 0.01. Therefore, both loci could be considered polymorphic. The first variant c.-582A > G was found in 49 promoter sequences (58 individuals were homozygous for allele A, 41 heterozygous AG and 4 homozygous GG), which resulted in an allele frequency of 0.762 for c.-582A and 0.218 for c.-582G. This variant corresponded to single nucleotide polymorphism (SNP) rs10421768 in the NCBI SNP database http://www.ncbi.nlm.nih.gov/snp/. The second variant c.-153C > T was found in only three heterozygous CT. Thus allele frequency of allele c.-153T was 1.4%. Furthermore, unique variants were found in two heterozygous individuals: c.-572C > T and c.-188C > T. No transcription factor binding sites were predicted for these unique variants http://www.biobase-international.com/pages/index.php?id=transfac.
c.-582A > G genotype and its relation with biochemical parameters
Serum iron, serum transferrin, serum ferritin and transferrin saturation according to HAMP c.-582A > G genotype in Galician population
c.-582A > G Genotype
A/G and G/G
SI mean ± SD (μg/dL)
55.6 ± 30.7
55.5 ± 31.2
ST mean ± SD (mg/dL)
265.7 ± 28.0
263.6 ± 30.8
SF mean ± SD (μg/dL)
95.9 ± 73.7
106.6 ± 74.5
TS mean ± SD (%)
15.0 ± 8.5
15.2 ± 8.8
Hepcidin is the main regulator of iron homeostasis and, therefore, a better understanding of its regulation is of great interest. In the present work, we sequenced, in an initial sample group, the promoter of the hepcidin gene, HAMP. This allowed the identification of two genetic polymorphic variants in our population. In a further group of individuals the c.-582A > G was genotyped and a possible relationship with biochemical iron status markers was studied. In addition, we performed functional-in-vitro analysis to characterize the effect of the c.-582A > G variant on the hepcidin gene transactivation.
In the present work, we first sequenced the promoter of the hepcidin gene HAMP in a group of 103 healthy individuals. This allowed the identification of two polymorphic sites. In another group of 224 individuals, the c.-582A > G variant was genotyped and its possible relationship with biochemical data as well as its effect on HAMP transactivation was studied. The frequency distribution of this SNP in the Galician population is similar to that of a cohort of Italian healthy donors  and furthermore to those of other European populations http://www.ncbi.nlm.nih.gov/snp/.
The c.-582A > G variant is located within an E-box, which is a responsive element for helix-loop-helix transcription factors (bHLH-ZIP). The canonical sequence of the E-box is CACGTG (consensus CANNTG). When the A is substituted by G, the E-box would be abolished and unable to bind the transcription factors. USF1 and USF2 are bHLH-ZIP transcription factors that have been shown to positively regulate hepcidin expression through E-boxes within the promoter . Therefore, we decided to study the in-vitro effect of the c.-582A > G variant on the HAMP transactivation in response to USF stimulation. Our results confirmed the ability of USF1 and USF2 to transactivate the HAMP promoter in HepG2 cells. Additionally, we found that the USF transactivation was maintained in the c.-582G promoter. However, the USF promoter transactivation was slightly reduced in the case of the c.-582G allele compared to c.-582A. The subtle decrease in transactivation might be caused by the fact that only one of the four E-boxes present in the promoter is affected. Nonetheless, it can not be ruled out that different endotoxin content or slightly incorrect DNA estimation in the plasmid preparations caused the small differences found between the two variants.
No significant association between the c.-582A > G genotype and serum iron, serum transferrin, transferrin saturation or ferritin levels were found, which might reflect no differences in liver iron concentration. These results are in agreement with Bruno and col. that did not find association either . Conversely, this association was found among a series of thalassemia major patients. The experiments performed in the HepG2 cells showed, at its maximum, a 20% transactivation reduction by the c.-582G allele, compared to c.-582A. In addition, one has to bear in mind that the in-vitro model used may not completely reflect what happens in-vivo. This data together with the fact that HAMP is regulated by many different signals [6, 10, 15–18, 20–22] suggested that the effect of the c.-582A > G polymorphism on hepcidin levels, and therefore on the iron status in-vivo is most likely negligible. We hypothesize that the c-.582A > G variant in the human HAMP promoter would have no effect on the hepcidin transcription in normal situations but it might have some effect in physiopathological situations where more hepcidin was needed. However, further studies would be necessary to confirm this hypothesis.
The variation at the second polymorphic site, c.-153C > T, in the promoter region of the hepcidin gene had been previously identified in a patient with massive iron overload homozygous for the p.Cys282Tyr mutation in HFE. Island and col. demonstrated that the c.-153C > T mutation, which is located in a BMP-Responsive Element, reduced basal hepcidin gene expression and impaired its response to BMPs and IL-6. This variation has been found at a very low frequency by other investigators: Island and col. did not found it in 200 chromosomes from healthy volunteers , Barton and col. found a frequency of 0.06% among 785 Haemochromatosis and Iron Overload Screening (HEIRS) Study participants . Both groups concluded that the variation in haemochromatosis patients is not worth to study. On the other hand, Aguilar-Martinez and col. found the variant with a 5% frequency in a group of haemochromatosis patients without the p.Cys282Tyr mutation in HFE, whilst in healthy individuals of the same population the allele frequency was only 0.4%. They suggested that this variant should be studied in a selected iron overload patients . Our finding of a 1.4% frequency of the c.-153T allele in the Galician population additionally supports this recommendation for genotyping this locus in iron overload Galician patients.
Due to the role of hepcidin as a key modulator of systemic iron levels, it is plausible that differences in the level of its expression may partly account for the phenotypic variant in iron metabolism between individuals. In this report, no association was found between serum iron, transferrin or ferritin levels and the c.-582A > G HAMP promoter variant in a healthy population. This result was supported by the functional in-vitro studies, which showed only a 10-20% reduction of the promoter transactivation by allele c.-582G compared to the c.-582A allele.
A total of 327 individuals were used for the present study. Initially, a sample of 103 healthy individuals was recruited from the university students and employees of the Hospital Clínico Universitario at Santiago de Compostela. Information about their parents' and grandparents' geographic origins was provided and all were from the Spanish Galician region, which is a relatively isolated European population at the westernmost continental edge. A second independent random sample of 224 individuals was randomly recruited from the general adult population of the Galician town La Estrada. Serum ferritin levels of this group were between 10-300 ng/mL. Detailed description has been reported elsewhere . Informed consent was obtained from each participant. All the procedures were conducted according to the principles expressed in the Declaration of Helsinki and to the ethical standards of the local Ethics Committee.
Serum iron and serum transferrin were assayed in an ADVIA-2400 Chemistry System (Siemens Medical Diagnostics, Germany) and serum ferritin was measured in an ADVIA-Centaur Immunoassay System (Siemens). Transferrin saturation (TS) was calculated with the following formula: TS (in %) = [serum iron (in μg/dL) × 71]/serum transferrin (in mg/dL). With these methods, normal concentrations in adults are as follows: serum iron, 45-160 μg/dL (males) and 40-150 μg/dL (females); serum transferrin, 215-365 mg/dL (males) and 250-380 mg/dL (females); transferrin saturation, 7-53%; and serum ferritin 10-295 ng/mL (males) and 25-325 ng/mL (females).
Genetic analysis of the HAMP promoter
Blood samples were collected and DNA was extracted from peripheral blood leucocytes following standard procedures. A fragment of 951bp that extended from nucleotide 1013 upstream the ATG translation initiation codon of HAMP to 9bp downstream of exon 1 was PCR amplified with the following pair of primers: forward 5'-GTACTCATCGGACTGTAGATGTTAGC-3' and reverse 5'-GTGACAGTCGCTTTTATGGGGCCTGC-3'. PCR products were purified and both strands were directly sequenced using the BigDye terminator kit and run in the 3730 }l DNA analyzer (Applied Biosystems, Foster City, CA, USA).
c.-582A > G variant genotype
The rs10421768 polymorphism located 582 bp upstream the ATG translation initiation codon was genotyped using the TaqMan® functionally tested assay c_2604942_10 in a 7300 Real-Time PCR system (Applied Biosystems).
The human hepatoma cell line HepG2 was obtained from the European Cell Culture Collection (ECACC, UK) and maintained in Dulbecco's modified media (DMEM) with 4500 mg/L glucose supplemented with 10% foetal bovine serum, 2 mM glutamine, 50 units/mL penicillin, 50 μg/mL streptomycin and 100 μg/mL neomycin.
Two 951 bp fragments, which included 942bp of the human HAMP promoter plus 9 bp of 5'UTR of exon 1 were PCR amplified from homozygous c.-582A and c.-582G genomic DNAs respectively (Forward 5'-GGCTCGAG GTACTCATCGGACTGTAGATGTTAGC-3', Reverse 5'-GGAAGCTT GTGACAGTCGCTTTTATGGGGCCTGC-3').
PCR products were double-digested with HindIII and XhoI restriction enzymes (Promega Corporation, Madison, WA, USA). VarA-Hep/Luc and VarG-Hep/Luc constructs were generated using T4 DNA ligase (Promega) by cloning digested products into the same restriction sites of pGL3-Control vector (Promega). The nucleotide sequence of the constructs was confirmed by DNA sequencing.
Plasmids pCMV-USF1 and pCMV-USF2, (kindly provided by Dr. M. Sawadogo University of Texas, MD USA and Dr. H K. Bayele University College London UK), were used for transactivation assays.
Approximately 1.4 × 105 HepG2 cells were plated in 12-well plates 24 h before transfection. Cells were transiently transfected with 1000 ng of DNA using FugeneHD transfection reagent at an 8:2 ratio according to manufacturer's protocol (Roche Applied Sciences, Indianapolis, USA). For promoter transactivation assays cells were transiently cotransfected with 350 ng of Hep/Luc constructs and 450 ng of pCMV-USF1, pCMV-USF2 or the combination of these latter two plasmids at a 1:1 ratio. In all cases, 200 ng of pRL-TK vector (Promega) was included as an internal control of transfection efficiency. When necessary, pcDNA3.1 empty vector was used to equalize the total amount of DNA. pGL3-Control vector was used in each experiment as an interassay control to normalize transfection efficiencies. Cells were incubated with transfection reagents for 48 h at 37°C in full DMEM media, followed by further 18 h incubation in fresh full media. Luciferase and Renilla activities were determined using Dual-Glo®Luciferase Assay System (Promega) according to manufacturer's protocol. Luminiscence was measured in a TRIAD LT luminometer (Dynex Technologies, Chantilly, VA, USA). All the experiments were independently performed in triplicate.
Transactivation of the HAMP promoter is expressed as fold induction of luciferase/renilla luminescence ratio of Hep/Luc expressing cells over the ratio obtained from cells transfected with pGL3-Control vector. Data are presented as mean values ± SEM.
Statistical analysis was performed using SPSS analysis software v.16.0 (SPSS Inc., Chicago, IL, USA). The differences found in the values of the biochemical indicators of iron status among the various genotypes at position c.-582 were tested by the Mann-Whitney and Kruskal-Wallis tests. The results from Luciferase assays were compared using one-way ANOVA for repeated measurements followed by Tukey's posthoc test. These analyses were performed using GraphPad Prism v.5.0 (GraphPad Software, San Diego, California). In all cases a p value less than 0.05 was considered significant.
The authors would like to thank Michelle Sawadogo, University of Texas and Henry K. Bayele, University College London for USF1 and USF2 expression plasmids.
This work was supported by a grant from the Fondo de Investigaciones Sanitarias del Instituto de Salud Carlos III (PI052249 to LL) and Xunta de Galicia (PGIDIT06PXIC9101136PN).
- Park CH, Valore EV, Waring AJ, Ganz T: Hepcidin, a urinary antimicrobial peptide synthesized in the liver. J Biol Chem. 2001, 276: 7806-7810. 10.1074/jbc.M008922200.View ArticlePubMedGoogle Scholar
- Fleming RE, Bacon BR: Orchestration of iron homeostasis. N Engl J Med. 2005, 352: 1741-1744. 10.1056/NEJMp048363.View ArticlePubMedGoogle Scholar
- Ganz T: Hepcidin-a regulator of intestinal iron absorption and iron recycling by macrophages. Best Pract Res Clin Haematol. 2005, 18: 171-182. 10.1016/j.beha.2004.08.020.View ArticlePubMedGoogle Scholar
- Leong WI, Lonnerdal B: Hepcidin, the recently identified peptide that appears to regulate iron absorption. J Nutr. 2004, 134: 1-4.PubMedGoogle Scholar
- Roy CN, Andrews NC: Anemia of inflammation: the hepcidin link. Curr Opin Hematol. 2005, 12: 107-111. 10.1097/00062752-200503000-00001.View ArticlePubMedGoogle Scholar
- 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. 2001, 7811-7819. 10.1074/jbc.M008923200. 16:Google Scholar
- 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. 2001, 17: 8780-8785. 10.1073/pnas.151179498.View ArticleGoogle Scholar
- Nicolas G, Bennoun M, Porteu A, Mativet S, Beaumont C, Grandchamp B, Sirito M, Sawadogo M, Kahn A, Vaulont S: Severe iron deficiency anemia in transgenic mice expressing liver hepcidin. Proc Natl Acad Sci USA. 2002, 99: 4596-4601. 10.1073/pnas.072632499.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. Nat Genet. 2003, 33: 21-22. 10.1038/ng1053.View ArticlePubMedGoogle Scholar
- Nemeth E, Rivera S, Gabayan V, Keller C, Taudorf S, Pedersen BK, Ganz T: IL-6 mediates hypoferremia of inflammation by inducing the synthesis of the iron regulatory hormone hepcidin. J Clin Invest. 2004, 113: 1271-1276.PubMed CentralView ArticlePubMedGoogle Scholar
- Island ML, Jouanolle AM, Mosser A, Deugnier Y, David V, Brissot P, Loréal O: A new mutation in the hepcidin promoter impairs its BMP response and contributes to a severe phenotype in HFE related hemochromatosis. Haematologica. 2009, 94: 720-724. 10.3324/haematol.2008.001784.PubMed CentralView ArticlePubMedGoogle Scholar
- Jacolot S, Le Gac G, Scotet V, Quere I, Mura C, Ferec C: HAMP as a modifier gene that increases the phenotypic expression of the HFE pC282Y homozygous genotype. Blood. 2004, 103: 2835-2840. 10.1182/blood-2003-10-3366.View ArticlePubMedGoogle Scholar
- Merryweather-Clarke AT, Cadet E, Bomford A, Capron D, Viprakasit V, Miller A, McHugh PJ, Chapman RW, Pointon JJ, Wimhurst VL, Livesey KJ, Tanphaichitr V, Rochette J, Robson KJ: Digenic inheritance of mutations in HAMP and HFE results in different types of haemochromatosis. Hum Mol Genet. 2003, 12: 2241-2247. 10.1093/hmg/ddg225.View ArticlePubMedGoogle Scholar
- Inamura J, Ikuta K, Jimbo J, Shindo M, Sato K, Torimoto Y, Kohgo Y: Upregulation of hepcidin by interleukin-1beta in human hepatoma cell lines. Hepatol Res. 2005, 33: 198-205. 10.1016/j.hepres.2005.08.005.View ArticlePubMedGoogle Scholar
- Babitt JL, Huang FW, Wrighting DM, Xia Y, Sidis Y, Samad TA, Campagna JA, Chung RT, Schneyer AL, Woolf CJ, Andrews NC, Lin HY: Bone morphogenetic protein signaling by hemojuvelin regulates hepcidin expression. Nat Genet. 2006, 38: 531-539. 10.1038/ng1777.View ArticlePubMedGoogle Scholar
- Fleming MD: The regulation of hepcidin and its effects on systemic and cellular iron metabolism. Hematology Am Soc Hematol Educ Program. 2008, 151-158.Google Scholar
- Schmidt PJ, Toran PT, Giannetti AM, Bjorkman PJ, Andrews NC: The transferrin receptor modulates Hfe-dependent regulation of hepcidin expression. Cell Metab. 2008, 7: 205-214. 10.1016/j.cmet.2007.11.016.PubMed CentralView ArticlePubMedGoogle Scholar
- Babitt JL, Huang FW, Xia Y, Sidis Y, Andrews NC, Lin HY: Modulation of bone morphogenetic protein signaling in vivo regulates systemic iron balance. J Clin Invest. 2007, 117: 1933-1939. 10.1172/JCI31342.PubMed CentralView ArticlePubMedGoogle Scholar
- Wrighting DM, Andrews NC: Interleukin-6 induces hepcidin expression through STAT3. Blood. 2006, 108: 3204-3209. 10.1182/blood-2006-06-027631.PubMed CentralView ArticlePubMedGoogle Scholar
- Peyssonnaux C, Zinkernagel AS, Schuepbach RA, Rankin E, Vaulont S, Haase VH, Nizet V, Johnson RS: Regulation of iron homeostasis by the hypoxia-inducible transcription factors (HIFs). J Clin Invest. 2007, 117: 1926-1932. 10.1172/JCI31370.PubMed CentralView ArticlePubMedGoogle Scholar
- Weizer-Stern O, Adamsky K, Amariglio N, Levin C, Koren A, Breuer W, Rachmilewitz E, Breda L, Rivella S, Cabantchik ZI, Rechavi G: Downregulation of hepcidin and haemojuvelin expression in the hepatocyte cell-line HepG2 induced by thalassaemic sera. Br J Haematol. 2006, 135: 129-138. 10.1111/j.1365-2141.2006.06258.x.View ArticlePubMedGoogle Scholar
- Du X, She E, Gelbart T, Truksa J, Lee P, Xia Y, Khovananth K, Mudd S, Mann N, Moresco EM, Beutler E, Beutler B: The serine protease TMPRSS6 is required to sense iron deficiency. Science. 2008, 320: 1088-1092. 10.1126/science.1157121.PubMed CentralView ArticlePubMedGoogle Scholar
- Andreani M, Radio FC, Testi M, De Bernardo C, Troiano M, Majore S, Bertucci P, Polchi P, Rosati R, Grammatico P: Association of hepcidin promoter c.-582 A > G variant and iron overload in thalassemia major. Haematologica. 2009, 94: 1293-1296. 10.3324/haematol.2009.006270.PubMed CentralView ArticlePubMedGoogle Scholar
- Bruno F, Bonalumi S, Camaschella C, Ferrari M, Cremonesi L: The -582A > G variant of the HAMP promoter is not associated with high serum ferritin levels in normal subjects. Haematologica. 2010, 95: 849-50. 10.3324/haematol.2009.018986.PubMed CentralView ArticlePubMedGoogle Scholar
- Bayele HK, McArdle H, Srai SK: Cis and trans regulation of hepcidin expression by upstream stimulatory factor. Blood. 2006, 108: 4237-4245. 10.1182/blood-2005-07-027037.View ArticlePubMedGoogle Scholar
- Barton JC, Leiendecker-Foster C, Li H, DelRio-LaFreniere S, Acton RT, Eckfeldt JH: HAMP promoter mutation nc.-153C > T in 785 HEIRS Study participants. Haematologica. 2009, 94: 1465-10.3324/haematol.2009.011486.PubMed CentralView ArticlePubMedGoogle Scholar
- Aguilar-Martinez P, Giansily-Blaizot M, Bismuth M, Cunat S, Igual H, Schved JF: HAMP promoter mutation nc.-153C > T in non p.C282Y homozygous patients with iron overload. Haematologica. 2010, 95: 687-688. 10.3324/haematol.2009.018028.PubMed CentralView ArticlePubMedGoogle Scholar
- Campos J, Gude F, Quinteiro C, Vidal C, Gonzalez-Quintela A: Gene by environment interaction: the -159C/T polymorphism in the promoter region of the CD14 gene modifies the effect of alcohol consumption on serum IgE levels. Alcohol Clin Exp Res. 2006, 30: 7-14. 10.1111/j.1530-0277.2006.00007.x.View ArticlePubMedGoogle Scholar
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.