- Research article
- Open Access
- Open Peer Review
A polymorphic variant of the insulin-like growth factor 1 (IGF-1) receptor correlates with male longevity in the Italian population: a genetic study and evaluation of circulating IGF-1 from the "Treviso Longeva (TRELONG)" study
BMC Geriatricsvolume 9, Article number: 19 (2009)
An attenuation of the insulin-like growth factor 1 (IGF-1) signaling has been associated with elongation of the lifespan in simple metazoan organisms and in rodents. In humans, IGF-1 level has an age-related modulation with a lower concentration in the elderly, depending on hormonal and genetic factors affecting the IGF-1 receptor gene (IGF-1R).
In an elderly population from North-eastern Italy (n = 668 subjects, age range 70–106 years) we investigated the IGF-1R polymorphism G3174A (rs2229765) and the plasma concentration of free IGF-1. Frequency distributions were compared using χ2-test "Goodness of Fit" test, and means were compared by one-way analysis of variance (ANOVA); multiple regression analysis was performed using JMP7 for SAS software (SAS Institute, USA). The limit of significance for genetic and biochemical comparison was set at α = 0.05.
Males showed an age-related increase in the A-allele of rs2229765 and a change in the plasma level of IGF-1, which dropped significantly after 85 years of age (85+ group). In the male 85+ group, A/A homozygous subjects had the lowest plasma IGF-1 level. We found no clear correlation between rs2229765 genotype and IGF-1 in the females.
These findings confirm the importance of the rs2229765 minor allele as a genetic predisposing factor for longevity in Italy where a sex-specific pattern for IGF-1 attenuation with ageing was found.
Human longevity is the result of a complex interaction between genetics and environment (including lifestyle) and the contribution of each factor is difficult to quantify. Genetic components account for from 20 to 30% of lifespan, according to twin-based population studies [1, 2].
In model organisms like Caenorhabditis elegans, Drosophila and the mouse, several genes involved in different basic cell pathways (stress response, DNA repair, microbial immunity and inflammation, metabolism and calorie restriction) have been shown to be related to longevity [3–5]. Genes belonging to evolution-conserved pathways, like the insulin/insulin like-growth factor 1 (IGF-1) network, are of particular interest in view of their key roles in cell physiology [6, 7].
Insulin/IGF-1 signaling (IIS) triggers intracellular downstream transcription factors (for instance, belonging to the FOXO family) as the final step in a signal transduction pathway mediated by the homologous membrane receptors of insulin and IGF-1 . Down-regulation of IIS has been associated with increased longevity, probably as a consequence of the decrease in cell growth and reduction in oxidative metabolism [9–11]. In humans the level of circulating IGF-1 is low in the elderly, but is higher in childhood than adulthood, probably as a consequence of different hormonal regulation. IGF-1 expression is sensitive to the growth hormone (GH) concentration and is also influenced by sex hormones [12–14].
Genetic variations in genes belonging to the IIS pathway have been explored in relation to longevity, dementia, metabolic diseases and cancer [15–18]. In the Italian population, a synonymous polymorphism (rs2229765) on the IGF-1 receptor gene (IGF-1R) [GenBank:NM_000875] consisting of a G to A transition at nucleotide 3174 leading to the amino acid change Glu->Glu at position 1043 (E1043E) [GenBank:NP_000876] was evaluated and the A-allele showed a positive association with aging and a negative influence on circulating levels of IGF-1 . The same polymorphism has also been reported as predisposing factor for ischemic stroke in China .
We have collected a sample of elderly people (n = 668) from Treviso, Northeastern Italy, from 70 to 106 years old, homogeneous for ethnicity and geographical origin, and we have assessed their rs2229765 genotype and free IGF-1 plasma level.
Population sample recruitment
A comprehensive description of the study design has been reported . Briefly, starting from a list of about 14,000 Treviso inhabitants over 70 years of age, a population sample was built up comprising three classes (70–79, 80–89, 90 years and over). The list from which the sample has been drawn comes from the population register. The final sample size was of 668 independent subjects (311 males and 357 females). The participants were evaluated from the biologic, clinical and socio-economic point of view, with a blood sample collection and a structured interview. The study protocol presenting the inclusion criteria, collecting procedure and questionnaire to be administered to all people involved in the project, was submitted to and approved by the Ethical Committee of the National Institute on Research and Care of the Elderly (INRCA, Italy). This protocol includes an informed, written consent to be obtained from the participant, or from a legally responsible relative in case of mentally impaired subjects, for clinical and genetic studies.
Blood sampling and rs2229765 genotyping
Blood samples (about 30 ml) were collected by venipuncture; one part was used to isolate leukocytes for gDNA extraction, and the other was centrifuged at 2000 rpm for 10 minutes at 4°C (with sodium EDTA as anticoagulant) to separate the plasma fraction, which was then divided into aliquots and stored at -80°C until required. From the final population of 668 subjects, consent for blood collection was obtained from 590 Treviso inhabitants and, in all, 587 plasma samples were successfully prepared and stored for further analysis.
Genomic DNA was extracted from leukocytes using a semi-automated nucleic acid extractor (AB6100, Applera, USA) and stored at 4°C. To assess the rs2229765 genotype, a gDNA aliquot (about 50 ng) was amplified by polymerase chain reaction (PCR) using the following primers: forward 5'-tcttctccagtgtacgttcc-3', reverse 5'-ggaactttctctttaccacatg-3'. The resulting PCR product was digested by MnlI (New England Biolabs Inc, UK). The rs2229765 genotype was examined after loading the corresponding enzymatic digestions on a capillary electrophoresis unit with standard reference markers for instrumental lining-up (Agilent Technologies, USA).
From the available 590 gDNA samples, 510 samples gave unambiguous digestion result and were considered for subsequent analysis.
Plasma IGF-1 assay
Plasma IGF-1 was assayed by a specific sandwich-type enzyme-linked immunosorbent assay (ELISA) (Diagnostic System Laboratories, Inc, USA), according to the manufacturer's instructions. The kit sensitivity was 10 ng/ml and the intra-assay %CV was <10%. From the assayed 587 plasma samples, 564 gave a IGF-1 value above the kit sensitivity limit and were considered for further analysis.
Frequency distributions were compared using χ2-test "Goodness of Fit" test (Preacher, K. J. Calculation for the chi-square test: An interactive calculation tool for chi-square tests of goodness of fit and independence [Computer software]. http://www.quantpsy.org) and means were compared by one-way analysis of variance (ANOVA) followed by Tukey's or Dunnett's post-hoc test. Analyses were done using StatView program ver. 5.0. Multiple regression analysis was performed using JMP7 for SAS software (SAS Institute, USA). The limit of significance for genetic and biochemical comparison was set at α = 0.05.
Age-dependent distribution of rs2229765
The demographic and the key clinical data of the population recruited in the TRELONG study are summarized in [Table 1]. The sample was representative of the general elderly Treviso population composition, ranging from 70 to over 100 years, with six individuals more than 100 years. Male-to-female ratio was 0.87. We divided the sample by five-year intervals as shown or by comparing two groups: the (70–85) group and the 85+ one, in an attempt to find longevity-related traits. We chose this cut-off value as several data suggest that after 85 years peculiar biochemical and metabolic features can be highlighted in comparison to the general population and also according to Italian population-specific survival curves [19, 22, 23]. This division gave roughly two thirds of the sample under 85 and one third over 85. Males and females were equally represented in the (70–85) group, while in the 85+ group women were over-represented (male-to-female ratio 0.67).
We first assessed the genotypic distribution of rs2229765 in our sample in five-year groups. No significant difference came to light (data not shown).
Then, we compared the rs2229765 genotypes in the (70–85) and the 85+ groups [Table 2]. The genotypic distributions respected Hardy-Weinberg equilibrium in both groups (data not shown) but there were no differences for either genotypic or allelic frequency in the whole sample. However, we found a sex-specific pattern as rs2229765 genotypic distribution in females under and over 85 was different from (70–85) males. Remarkably, at allelic level (70–85) males had a reduced A frequency (34.4%) in comparison to 85+ males (43.7%).
Plasma IGF-1 concentration and relation between IGF-1 level and rs2229765 genotype
We measured plasma IGF-1 concentrations in our whole sample and plotted the results according to age but regardless of sex [Figure 1]. The IGF-1 concentration decreased with age, and simple linear regression showed a significant negative correlation between age and IGF-1 level (r = -0.18, p < 0.001, sample size n = 564). To control for confounding factors, we then performed a multiple regression analysis using IGF-1 level as dependent variable and age, sex, body mass index (BMI), smoking, hypertension and blood glucose level as covariates [Table 3]. The correlation was confirmed with an even stronger significance level (p < 0.0001).
The next step was to divide the population into five-year groups to check the longevity-related reduction in more detail [Table 4]. There was a significant decrease starting from 86 years in comparison to the 70–75 group, but the 91–95 group's mean IGF-1 level was not different from the first age bracket considered. However, the 96–100 group IGF-1 plasma level was significantly lower than the 70–75 group (a 24% reduction). This difference was even more evident in the extremely long-living subjects (40% less than the 70–75 group) but the small sample size did not give sufficient statistical power. This age-related IGF-1 reduction was evident by population stratification in ten-year groups, too [Table 4].
To verify whether the age-dependent effect on IGF-1 was influenced by sex, we compared IGF-1 levels by sex and five-year intervals [Table 5]. We found that in males above 85 years IGF-1 decreased significantly while the mean plasma IGF-1 was the same in females under and over this cut-off value. This sex-specific effect was present also by ten-year intervals [Table 5].
To check for any correlation between rs2229765 genotype and mean IGF-1 level we stratified the sample by rs2229765 genotype. No significant change in IGF-1 concentration came to light (IGF-1 mean ± standard deviation: G/G 186 ± 78, n = 174; G/A 175 ± 81, n = 258; A/A 179 ± 83, n = 78; p > 0.05, one-way ANOVA). When we did a similar analysis in the (70–85) and 85+ groups we found a gender-specific pattern [Figure 2A]. In fact, in the men over 85 the homozygous A class had the lowest IGF-1 levels (119 ± 50 ng/mL, n = 15). This reduction was also present in the heterozygous G/A group (160 ± 74 ng/mL, n = 48) in comparison to the homozygous G/G class (185 ± 74 ng/mL, n = 26). However, the above reported difference was not robust (p = 0.048 for 85+ males A/A in comparison to G/G group in the same age bracket) and a chance effect cannot be excluded. We have also plotted the individual male IGF-1 levels for the whole population according to rs2229765 genotype and we have calculated the linear regression using age as independent variable. We noticed a genotype effect as the faster age-related decrease correlated with the A/A genotype while the G/G genotype decreased more slowly [Figure 2B]. In details, slope, r value and significance for each genotype were the following: (G/G) slope -1.6; r = 0.15; p = 0.15; (G/A) slope -3.0; r = 0.27; p = 0.002; (A/A) slope -5.0; r = 0.54; p = 0.03. A similar regression analysis was performed for females and the regression parameters were: (G/G) slope 0.13; r = 0.003.; p = 0.97; (G/A) slope -2.8; r = 0.29; p = 0.0007; (A/A) slope -1.6; r = 0.19, p = 0.09.
Human longevity is influenced by genetic and environmental factors. The present study addressed the question whether the genetic variability due to rs2229765 polymorphism of the IGF-1R plays a role in human longevity. We also sought a correlation between these genetic data and the circulating level of IGF-1, an important molecule in the general cell metabolism, whose reduction was found to be beneficial for longevity in model organisms [5–7].
We collected a homogeneous sample of people from Treviso (Italy), aged from 70 to 106 years (TRELONG Study) . The sample was of particular interest to highlight the specific features of extremely long-living individuals (over 85 years) in comparison to those with a normal lifespan (between 70 and 85 years).
The polymorphism we assessed is a synonymous amino acidic substitution (E1043E) in the IGF-1R protein, whose functional significance is unclear. Genetic data about the rs2229765 did not indicate any major differences when the population was grouped in five-year age brackets or split in the two groups (70–85) and 85+, but there was an increase of the A-allele in the males over 85, suggesting a sex-specific effect on male longevity. The significance of this increase was not robust, (p = 0.04, χ2-test) probably because of the sample size of 85+ males. The (70–85) male group had a lower A-allele frequency than (70–85) females. This might be interpreted as a sex-specific effect but might also hide a chance or recruitment bias. A confirmatory study in an independent male population sample from Veneto Region is required to better clarify the above described limitations of the study. Our finding about rs2229765 A-allele was partially in agreement with a previous study in the Italian population by Bonafè et al., which reported an increase of the rs2229765 A-allele in individuals over 85 years of age regardless of sex. This aspect can be explained considering that our study design was different, based on an elderly population sample, while Bonafè et al. considered a wider sample ranging from 20 to 100+ years . Consequently, their increase of the A-allele in all the people over 85 years old might have depended on a larger difference between male A+ carriers in the two groups and does not exclude a sex-specific pattern. A gender-specific influence on longevity is a well-known aspect that can be partially explained at genetic level. Other polymorphisms show a gender-specific pattern in relation to longevity as well [24–26].
The second part of our study investigated plasma IGF-1 concentration, and a possible correlation between rs2229765 and IGF-1 circulating levels. We found a linear relationship between plasma IGF-1 and age using a simple and a multivariate regression model, confirming that IGF-1 continues to fall through life, however long. A substantial gender-specific effect was embedded in the general population analysis, as men and women contributed differently to the IGF-1 reduction: female IGF-1 in the (70–85) group was lower than in males but it did not drop further over time; in contrast, male IGF-1 was significantly lower in the 85+ group than the (70–85) group. Consequently, the drop in IGF-1 over time in our population was mainly due to a specific effect in males. This agrees with a reported sex-specific variation of IGF-1 level in model organisms, where it is influenced by sex hormones [27–29]. Bonafé et al. also reported a linear reduction of circulating free IGF-1 in the Italian population , in accordance with our data that magnified the effect for the elderly and found previously unreported sex-specific difference.
We found no difference in the entire sample as regards a correlation between IGF-1 plasma level and rs2229765 genotype. However, men showed a genotype-specific age-related decrease of plasma IGF-1 that declined faster when the subject was homozygous for the A-allele. In females this association was less evident, however the A-allele carriers had a faster plasma IGF-1 decline over age than the G/G subjects. This might suggest that in females rs2229765 impacts on circulating IGF-1 level, but this is masked by other genetic or environmental factors. This sex-specific relation was confirmed when we split our population at 85 years and divided men from women, as we observed that in the 85+ males the reduction of mean IGF-1 was unequally distributed among the genotypic classes G/G, G/A and A/A, the last group having the lowest IGF-1 concentration (119 ± 50 ng/mL), and the G/G genotype the highest (185 ± 74 ng/mL) that was comparable to the mean IGF-1 value of (70–85) men (205 ± 81 ng/mL). However, this correlation was statistically weak (p = 0.048, Tukey's test) and a chance effect should be taken into account. Recently, a genetic screening in Ashkenazi Jewish centenarians reported gender-specific modulations of circulating IGF-1 that was associated to IGF-1R genetic variability (including rs2229765), confirming that also this synonymous genetic alterations in the human IGF-1R might alter IGF signaling pathway and human longevity .
This genetic and biochemical study in the Italian population confirmed and better detailed a genetic contribution to longevity coming from the A-allele of the rs2229765 polymorphism in the IGF-1R gene, that correlated with male longevity.
Herskind A, McGue M, Holm N, Sørensen T, Harvald B, Vaupel J: The heritability of human longevity: a population-based study of 2872 Danish twin pairs born 1870–1900. Hum Genet. 1996, 97: 319-323. 10.1007/BF02185763.
vB Hjelmborg J, Iachine I, Skytthe A, Vaupel JW, McGue M, Koskenvuo M, Kaprio J, Pedersen NL, Christensen K: Genetic influence on human lifespan and longevity. Hum Genet. 2006, 119: 312-321. 10.1007/s00439-006-0144-y.
Antebi A: Genetics of aging in Caenorhabditis elegans. PLoS Genet. 2007, 3: 1565-1571. 10.1371/journal.pgen.0030129.
Herndon LA, Schmeissner PJ, Dudaronek JM, Brown PA, Listner KM, Sakano Y, Paupard MC, Hall DH, Driscoll M: Stochastic and genetic factors influence tissue-specific decline in ageing C. elegans. Nature. 2002, 419: 808-814. 10.1038/nature01135.
Kim SK: Common aging pathways in worms, flies, mice and humans. J Exp Biol. 2007, 210: 1607-1612. 10.1242/jeb.004887.
Samuelson AV, Carr CE, Ruvkun G: Gene activities that mediate increased life span of C. elegans insulin-like signaling mutants. Genes Dev. 2007, 21: 2976-2994. 10.1101/gad.1588907.
Selman C, Lingard S, Choudhury AI, Batterham RL, Claret M, Clements M, Ramadani F, Okkenhaug K, Schuster E, Blanc E, Piper MD, Al-Qassab H, Speakman JR, Carmignac D, Robinson IC, Thornton JM, Gems D, Partridge L, Withers DJ: Evidence for lifespan extension and delayed age-related biomarkers in insulin receptor substrate 1 null mice. FASEB J. 2008, 22: 807-818. 10.1096/fj.07-9261com.
Cameron AR, Anton S, Melville L, Houston NP, Dayal S, McDougall GJ, Stewart D, Rena G: Black tea polyphenols mimic insulin/insulin-like growth factor-1 signalling to the longevity factor FOXO1a. Aging Cell. 2008, 7: 69-77. 10.1111/j.1474-9726.2007.00353.x.
Murphy CT, McCarroll SA, Bargmann CI, Fraser A, Kamath RS, Ahringer J, Li H, Kenyon C: Genes that act downstream of DAF-16 to influence the lifespan of Caenorhabditis elegans. Nature. 2003, 424: 277-283. 10.1038/nature01789.
Pinkston-Gosse J, Kenyon C: DAF-16/FOXO targets genes that regulate tumor growth in Caenorhabditis elegans. Nat Genet. 2007, 39: 1403-1409. 10.1038/ng.2007.1.
van Heemst D, Beekman M, Mooijaart SP, Heijmans BT, Brandt BW, Zwaan BJ, Slagboom PE, Westendorp RG: Reduced insulin/IGF-1 signalling and human longevity. Aging Cell. 2005, 4: 79-85. 10.1111/j.1474-9728.2005.00148.x.
Hall K, Hilding A, Thorén M: Determinants of circulating insulin-like growth factor-I. J Endocrinol Invest. 1999, 22: 48-57.
Münzer T, Rosen CJ, Harman SM, Pabst KM, St Clair C, Sorkin JD, Blackman MR: Effects of GH and/or sex steroids on circulating IGF-I and IGFBPs in healthy, aged women and men. Am J Physiol Endocrinol Metab. 2006, 290: E1006-1013. 10.1152/ajpendo.00166.2005.
Ruiz-Torres A, Soares de Melo Kirzner M: Ageing and longevity are related to growth hormone/insulin-like growth factor-1 secretion. Gerontology. 2002, 48: 401-407. 10.1159/000065507.
Garcia J, Ahmadi A, Wonnacott A, Sutcliffe W, Nagga K, Soderkvist P, Marcusson J: Association of insulin-like growth factor-1 receptor polymorphism in dementia. Dement Geriatr Cogn Disord. 2006, 22: 439-444. 10.1159/000095803.
Kojima T, Kamei H, Aizu T, Arai Y, Takayama M, Nakazawa S, Ebihara Y, Inagaki H, Masui Y, Gondo Y, Sakaki Y, Hirose N: Association analysis between longevity in the Japanese population and polymorphic variants of genes involved in insulin and insulin-like growth factor 1 signaling pathways. Exp Gerontol. 2004, 39: 1595-1598. 10.1016/j.exger.2004.05.007.
Lunetta KL, D'Agostino RB, Karasik D, Benjamin EJ, Guo CY, Govindaraju R, Kiel DP, Kelly-Hayes M, Massaro JM, Pencina MJ, Seshadri S, Murabito JM: Genetic correlates of longevity and selected age-related phenotypes: a genome-wide association study in the Framingham Study. BMC Med Genet. 2007, 8 (Suppl 1): S13-10.1186/1471-2350-8-S1-S13.
Rodriguez S, Gaunt TR, Day IN: Molecular genetics of human growth hormone, insulin-like growth factors and their pathways in common disease. Hum Genet. 2007, 122: 1-21. 10.1007/s00439-007-0378-3.
Bonafè M, Barbieri M, Marchegiani F, Olivieri F, Ragno E, Giampieri C, Mugianesi E, Centurelli M, Franceschi C, Paolisso G: Polymorphic variants of insulin-like growth factor I (IGF-I) receptor and phosphoinositide 3-kinase genes affect IGF-I plasma levels and human longevity: cues for an evolutionarily conserved mechanism of life span control. J Clin Endocrinol Metab. 2003, 88: 3299-3304. 10.1210/jc.2002-021810.
Cheng J, Liu J, Li X, Peng J, Han S, Zhang R, Xu Y, Nie S: Insulin-like growth factor-1 receptor polymorphism and ischemic stroke: a case-control study in Chinese population. Acta Neurol Scand. 2008, 118: 333-8. 10.1111/j.1600-0404.2008.01040.x.
Gallucci M, Ongaro F, Bresolin F, Bernardi U, Salvato C, Minello A, Amici GP, Barasciutti E, Mazzuco S, Gajo GB, De Angeli S, Forloni GL, Albani D, Zanardo A, Regini C: The Treviso Longeva (Trelong) study: a biomedical, demographic, economic and social investigation on people 70 years and over in a typical town of North-East of Italy. Arch Gerontol Geriatr. 2007, 44: 173-192. 10.1016/j.archger.2007.01.025.
Yashin AI, De Benedictis G, Vaupel JW, Tan Q, Andreev KF, Iachine IA, Bonafe' M, De Luca M, Valensin S, Carotenuto L, Franceschi C: Genes, demography and life span: the contribution of demographic data and genetic studies on ageing and longevity. Am J Hum Genet. 1999, 65: 1178-1193. 10.1086/302572.
Barbieri M, Bonafè M, Rizzo MR, Ragno E, Olivieri F, Marchegiani F, Franceschi C, Paolisso G: Gender specific association of genetic variation in peroxisome proliferator-activated receptor (PPAR)gamma-2 with longevity. Exp Gerontol. 2004, 39: 1095-1100. 10.1016/j.exger.2004.03.034.
Cederholm T, Persson M, Andersson P, Stenvinkel P, Nordfors L, Madden J, Vedin I, Wretlind B, Grimble RF, Palmblad J: Polymorphisms in cytokine genes influence long-term survival differently in elderly male and female patients. J Intern Med. 2007, 262: 215-223. 10.1111/j.1365-2796.2007.01803.x.
Lio D, Scola L, Crivello A, Colonna-Romano G, Candore G, Bonafè M, Cavallone L, Franceschi C, Caruso C: Gender-specific association between -1082 IL-10 promoter polymorphism and longevity. Genes Immun. 2002, 3: 30-33. 10.1038/sj.gene.6363827.
Passarino G, Montesanto A, Dato S, Giordano S, Domma F, Mari V, Feraco E, De Benedictis G: Sex and age specificity of susceptibility genes modulating survival at old age. Hum Hered. 2006, 62: 213-220. 10.1159/000097305.
Harper JM, Wolf N, Galecki AT, Pinkosky SL, Miller RA: Hormone levels and cataract scores as sex-specific, mid-life predictors of longevity in genetically heterogeneous mice. Mech Ageing Dev. 2003, 124: 801-810. 10.1016/S0047-6374(03)00133-7.
McCulloch D, Gems D: Sex-specific effects of the DAF-12 steroid receptor on aging in Caenorhabditis elegans. Ann N Y Acad Sci. 2007, 1119: 253-259. 10.1196/annals.1404.018.
Shved N, Berishvili G, D'Cotta H, Baroiller JF, Segner H, Eppler E, Reinecke M: Ethinylestradiol differentially interferes with IGF-I in liver and extrahepatic sites during development of male and female bony fish. J Endocrinol. 2007, 195: 513-523. 10.1677/JOE-07-0295.
Suh Y, Atzmon G, Cho MO, Hwang D, Liu B, Leahy DJ, Barzilai N, Cohen P: Functionally significant insulin-like growth factor I receptor mutations in centenarians. Proc Natl Acad Sci USA. 2008, 105: 3438-3842. 10.1073/pnas.0705467105.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2318/9/19/prepub
This study was supported by grants from the Veneto Region, the Treviso Municipality, Treviso Province and Cassamarca Foundation. We are grateful to all people who kindly agreed to participate in this study and to J.D. Baggott for manuscript editing. LP is recipient of a fellowship from "Golgi Cenci Foundation", Abbiategrasso, Milan, Italy.
The authors declare that they have no competing interests.
DA drafted the paper and analyzed the results; MG drafted the paper, designed the study and performed the multivariate analysis; SB, LP, AV and MP performed the genotyping and ELISA assays; GBG, SDA, AZ collected the blood and plasma samples; GF critically read the manuscript.