- Research article
- Open Access
- Open Peer Review
Genetic polymorphisms located in genes related to immune and inflammatory processes are associated with end-stage renal disease: a preliminary study
BMC Medical Geneticsvolume 13, Article number: 58 (2012)
Chronic kidney disease progression has been linked to pro-inflammatory cytokines and markers of inflammation. These markers are also elevated in end-stage renal disease (ESRD), which constitutes a serious public health problem.
To investigate whether single nucleotide polymorphisms (SNPs) located in genes related to immune and inflammatory processes, could be associated with ESRD development.
Design and methods
A retrospective case-control study was carried out on 276 patients with ESRD and 288 control subjects. Forty-eight SNPs were genotyped via SNPlex platform. Logistic regression was used to assess the relationship between each sigle polymorphism and the development of ESRD.
Four polymorphisms showed association with ESRD: rs1801275 in the interleukin 4 receptor (IL4R) gene (OR: 0.66 (95%CI = 0.46-0.95); p = 0.025; overdominant model), rs4586 in chemokine (C-C motif) ligand 2 (CCL2) gene (OR: 0.70 (95%CI = 0.54-0.90); p = 0.005; additive model), rs301640 located in an intergenic binding site for signal transducer and activator of transcription 4 (STAT4) (OR: 1.82 (95%CI = 1.17-2.83); p = 0.006; additive model) and rs7830 in the nitric oxide synthase 3 (NOS3) gene (OR: 1.31 (95%CI = 1.01-1.71); p = 0.043; additive model). After adjusting for multiple testing, results lost significance.
Our preliminary data suggest that four genetic polymorphisms located in genes related to inflammation and immune processes could help to predict the risk of developing ESRD.
The prevalence of chronic kidney disease (CKD) and end-stage renal disease (ESRD) is growing worldwide, with the overall prevalence of CKD stages 3-5 being 6.83% in Spain . Therefore, it constitutes a serious public health problem, which causes substantial morbidity and mortality . Traditional risk factors for CKD progression include persistent proteinuria, dyslipidaemia, hypertension and smoking [3, 4]. However, it has been postulated that non-traditional risk factors, such as oxidative stress, inflammation and immune processes, may also be important contributors to the pathogenesis of cardiovascular disease as well as progression to ESRD . The inflammatory response involved in renal damage produces a release of proinflammatory cytokines and chemokines, which cause an increased inflow of leukocytes, intensification of interstitial nephritis and progressive fibrosis .
On the other hand, genetic susceptibility is also considered an important determining factor for the appearance and/or progression of ESRD and its complications [7, 8]. Therefore, the study of polymorphisms may help us to look further into the disease pathogenesis and its underlying causes, as well as to predict the predisposition to developing the disease in order to identify the at-risk population. The aim of this study was to investigate whether single nucleotide polymorphisms (SNPs) located in regions affecting immune and inflammatory processes could be associated with development of ESRD. Our preliminary data suggest that four polymorphisms showed an association with ESRD protection or development, which could help to predict the risk of developing ESRD.
A retrospective study was carried out on 276 kidney transplant recipients (case-group) and 288 normal subjects (control-group) from the Hospital Clínico Universitario of Valladolid (Spain). The study was conducted in accordance with the Declaration of Helsinki. All patients and controls gave their written consent and the Institutional Ethics Committee approved the study.
Subjects in the case-group were patients older than 18 with ESRD that received a cadaver renal graft between December 1995 and October 2008. Subjects from the control-group were patients that underwent routine analysis at the hospital and had no evidence of renal pathology or other immune or inflammatory diseases. Individuals were selected to have age and gender distribution similar to the case-group patients at the time of transplantation. In this retrospective study, the clinical data and DNA samples were collected in a transversal way for both the case-group and the control-group between June 2008 and December 2008. In order to ensure homogeneity, all patients were Caucasian.
DNA extraction, marker selection and genotyping
DNA was extracted from whole blood by the Chemagic Magnetic Separador Module1, CHEMAGEN® which uses a magnetic particles system to obtain DNA. The quantity of recovered DNA was quantified by using PicoGreen® dsDNA Quantitation Reagent (Molecular Probes, Inc., Eugene, Oregon, USA).
All SNPs were genotyped at the Spanish National Genotyping Centre (CeGen; http://www.cegen.org/) by the SNPlex genotyping system 48-plex (Applied Biosystems, Foster City, CA, USA) following the manufacturer’s recommended protocol (http://www3.appliedbiosystems.com/cms/groups/mcb_support/documents/generaldocuments/cms_042019.pdf). As quality control, two Centre d’Etude du Polymorphisme Humain (CEPH) samples (NA10860 and NA10861) from the HapMap database were included in all genotyping assays .
Published SNPs were selected by different criteria (see Additional file 1: Table S1 Content 1): SNPs associated with kidney disease or kidney complications in previous articles and/or SNPs located in genes involved in genetic pathways related to prevention or susceptibility to ESRD (inflammatory and immune processes),. With this purpose several databases were consulted: International HapMap Project (http://www.hapmap.org), dbSNP databases (http://www.ncbi.nlm.nih.gov/SNP/), as well as pathway databases such as GeneOntology (http://www.geneontology.org/) and KEGG (http://www.genome.jp/kegg/pathway.html). Besides, for each significant SNP, the biological implications were analyzed "in silico" via web-tools: SMART (http://smart.embl-heidelberg.de/) for identifying the protein domains, PSIPRED (http://bioinf.cs.ucl.ac.uk/psipred/) for analyzing protein secondary structure prediction, SIFT (http://blocks.fhcrc.org/sift/SIFT.html) for studying the amino acid change toleration, PATROCLES (http://www.patrocles.org/) for identifying putative microRNA binding sites, PROMO (http://alggen.lsi.upc.es/) for transcription factor binding site prediction and Human Splicing Finder v2.4.1 (http://www.umd.be/HSF/) for studying the pre-mRNA splicing sites.
In the control group, Hardy-Weinberg equilibrium (HWE) was assessed for each polymorphism using the Pearson Chi-square (χ2) statistic. Only SNPs that observed HWE were included in association analysis.
Unconditional logistic regression was used to assess the relationship between each single polymorphism and the development of ESRD. Five genetic models were fitted (dominant, co-dominant, over-dominant, recessive, and log-additive model). For each SNP, likelihood ratio test (LRT), and Akaike's information criteria (AIC) were used as measures of the goodness of fit between models and to choose the inheritance model that best fits the data. The next equation defines the logistic model: log(p/1-p) = α + βG + yZ; with p being the probability, G the categorical variable with the polymorphisms codified, Z the variables to adjust the model (α, β and γ must be estimated). Odds ratios (OR) and 95% confidence intervals (CI) were calculated to check the relative risk for association. This model was adjusted for age and gender by SNPStats (http://bioinfo.iconcologia.net/SNPStats_web)  (p-value < 0.05 was considered significant). Moreover, multiple testing corrections were performed by SFDR (Stratified False Discovery Rate) software version 1.6. (http://www.utstat.utoronto.ca/sun/Software/SFDR/index.html).
Table 1 shows the clinical characteristics of all cases and controls. The median age and male female ratio were similar between groups. The most predominant primary disease involving renal failure was glomerulonephritis (≈30%).
From 48 SNPs, the genotyping assay for six SNPs did not work (rs1570360, rs1808593, rs2070744, rs4762, rs5498, rs55634318); one SNP was monomorphic (rs175176); and five SNPs (rs1799969, rs1800471, rs1800871, rs4311, rs699) were discarded because they exceeded 10% of missing values. Therefore, we analyzed 36 SNPs, all of which fulfilled the minimum allele frequency (MAF) >0.05 for all samples and were in HWE in the control group.
Significant association with ESRD was found for 4 SNPs (Table 2): AG genotype of rs1801275 in interleukin 4 receptor (IL4R) showed reduced odds of ESRD assuming an overdominant model (adjusted OR = 0.66 (95%CI = 0.46-0.95), p = 0.025). In the case of rs4586 in chemokine (C-C motif) ligand 2 (CCL2), the presence of each additional copy of the minor allele was associated with reduced odds of ESRD (adjusted OR = 0.70 (95%CI = 0.54-0.90), p = 0.005) in a log-additive model. On the other hand, rs301640 in an intergenic binding site for signal transducer and activator of transcription 4 (STAT4) and rs7830 in nitric oxide synthase 3 (NOS3) were associated with elevated odds of ESRD assuming a log-additive model (adjusted OR = 1.82 (95%CI = 1.17-2.83), p = 0.006; and adjusted OR = 1.31 (95%CI = 1.01-1.71), p = 0.043, respectively). After applying the false discovery rate (FDR) for multiple test correction, we obtained that adjusted p values for each hypothesis testing were no significant (q-value > 0.05). Allelic and genotypic frequencies of significant SNPs are showed in Additional file 2: Table S2 Content 2.
It has been previously shown that biomarkers of inflammation are high even in the early stages of CKD [11, 12]. The increased inflammation could be caused by both genetic predisposition and environmental factors, being linked to the risk of CKD progression to ESRD. In our study, we have found strong suggestion of associations between ESRD and four SNPs located in IL4R (rs1801275), CCL2 (rs4586), NOS3 (rs7830), and an intergenic binding site for STAT4 (rs301640). All of these sites involve genes related to inflammation and immune response pathways.
On the one hand, two of the significant SNPs (rs1801275, and rs4586) reflected a certain protection against ESRD:
rs1801275 (A/G) is located at position g.54150A > G of the IL4R gene in chromosome 16. This position corresponds to exon 12, where it generates a missense change (p.Gln576Arg). IL-4 is a cytokine involved in Th2 immune response, and its effect depends on binding with IL4R. IL-4 has been found to be an important predictor of kidney injury. Modulation of the IL4 pathway during ESRD may be due to an intracellular regulation involving different pathways but also could be possible that polymorphisms within the IL4R gene could alter the signalling pathway of IL-4, leading to a progression or prevention of kidney damage . In our study, the presence of the AG genotype seems to indicate protection from ESRD development, which could be caused by a decrease in signalling of IL-4 through of its receptor. Although rs1801275 is a missense SNP which is located in a low complexity domain of the protein (SMART) and it has been predicted as tolerant (SIFT). To date, this SNP has shown associations with inflammatory and autoimmune diseases such as arthritis rheumatoid, asthma, allergic diseases and atopy [14, 15], but the entire role in CKD of IL-4 and its signalling pathway through IL4R still remains unknown and further studies are needed.
rs4586 (C/T) is a synonymous polymorphism located on exon 2 of the CCL2 gene in chromosome 17. The CCL2 gene encodes a key chemokine in recruiting mononuclear inflammatory cells to sites of inflammation. One example is the interstitium and glomerulus, where CCL2 causes renal interstitial and glomerular inflammation, leading to progressive renal injury . In our study, the presence of each additional copy of the C allele of rs4586 indicates protection against ESRD. Recent studies have suggested a beneficial effect of blocking the action of CCL2 on diabetic nephropathy and renal function through anti-fibrotic effects . By analyzing the CCL2 sequence via PATROCLES, we have found that the C allele of rs4586 generates a putative target site (TGCTGCTA) for six different microRNAs (hsa-miR-15a/15b/16/195/424/497), whereas the T allele disrupt this target site and consequently none of these microRNAs target this sequence. MicroRNAs are small RNA molecules (22 nucleotides) that have a great impact on posttranscriptional regulation and potentially large relevance to complex diseases. When a microRNA attaches to its target, it can silence expression via mRNA degradation or by preventing mRNA translation . Therefore we could hypothesize that C allele might exert its beneficial effect by binding a microRNA and blocking the transcription of CCL2 gene.
On the other hand, rs301640 and rs7830 polymorphisms indicated susceptibility to ESRD:
rs301640 (A/G) is a SNP located on chromosome 13, in an intergenic region between the eukaryotic translation initiation factor 4A1 pseudogene 6 (EIF4A1P6) and the gene encoding hsa-miR-3169. This SNP is located within a binding site for the transcription factor STAT4. In our study, the presence of each additional copy of A allele was associated with increased odds of ESRD. By using PROMO software [19, 20], we found that the A allele might disrupt the STAT4 binding site, which could potentially modify a distal enhancer. STAT4 plays an important role in Th1 differentiation by transmitting IL-12 signals to produce IFN-γ , which could induce pro-inflammatory cytokines leading to injury in target tissues.
rs7830 (G/T) is located in a region that belongs to two different genes. The sense strand corresponds to intron 26-27 of the NOS3 gene, whereas the antisense strand matches with the 3´-UTR of the autophagy related 9 homolog B (ATG9B) gene. On the one hand, NOS3 encodes for an enzyme that generates NO in endothelial cells and is involved in the regulation of vascular function . On the other hand, ATG9B encodes for a protein required for autophagy in several eukaryotic organisms, although its entire function is unknown . In our study, the T allele of rs7830 was associated with ESRD. We hypothesize that the effects of rs7830 on ESRD might be due to NOS3 rather than ATG9B. In fact, polymorphisms in NOS3 have been associated with atherosclerotic vascular diseases , renal dysfunction , and advanced diabetic nephropathy . Moreover, this SNP seems to be located within a silencer motif (TGGGGACT) , where the G to T allele change could disturb the splicing mechanism in NOS3 leading to different transcripts . Thus, this change might affect NOS3 expression and to be associated with the development of ESRD.
Although these findings have also been supported by previous studies where the SNPs described above have been associated with kidney disease or related (Additional file 1: Table S1 Content 1), we could not find any significant results after applying FDR correction. In regards to this, since p-value is depending on the sample size, it may be possible that we have not found any significant adjusted p-value because our sample size is not large enough to detect moderate effects. Thus only big effects would be detected in small populations. Moreover, it exists some controversy about adjusting the “p-value” after multiple tests on clinical-orientated studies [29, 30]. In addition, the weak association found for the studied polymorphisms could be due to an indirect involvement in ESRD. That is, it cannot be discarded that these SNPs might be in linkage disequilibrium with more powerful polymorphisms associated with ESRD. Another limitation of our study was that although we tried to select individuals for control-group with similar age to case-group, the comparison between means was significant (50.0 versus 52.8 years; p = 0.030, Table 1). We think that these differences are so slight that have a low clinical significance. However, logistic regression analysis was adjusted by age in order to avoid possible age interferences.
Therefore, this is a preliminary study which could be considered to generate hypothesis for future studies. In fact, the association of these SNPs with ESRD needs to be confirmed by replicating studies with a larger sample size, as well as functional studies should be performed in order to get further insights into ESRD susceptibility. Moreover, it will be interesting to include in the future some factors which may also influence ESRD development (diabetes, dyslipidemia, smoking, etc.) which could not be collected in this study. Besides, those patients with ESRD who are receiving dialysis and/or have not yet received kidney transplant should also be taken into account in further studies because they could carry on different genotypes.
In conclusion, our preliminary data suggest that four polymorphisms (rs1801275, rs301640, rs4586, rs7830) related to inflammatory and immunity processes showed an association with protection or development of ESRD. These results could help to predict the risk of developing the disease, and also to improve the understanding of the pathways involved in the disease pathogenesis.
Akaike’s information criteria
autophagy related 9 homolog B
chemokine (C-C motif) ligand 2
chronic kidney disease
eukaryotic translation initiation factor 4A1 pseudogene 6
End-stage renal disease
false discovery rate
interleukin 4 receptor
likelihood ratio test
nitric oxide synthase 3
single nucleotide polymorphism
signal transducer and activator of transcription 4.
Otero A, de Francisco A, Gayoso P, Garcia F: Prevalence of chronic renal disease in Spain: results of the EPIRCE study. Nefrologia. 2010, 30 (1): 78-86.
Kottgen A: Genome-wide association studies in nephrology research. Am J Kidney Dis. 2010, 56 (4): 743-758. 10.1053/j.ajkd.2010.05.018.
Gansevoort RT, Matsushita K, van der Velde M, Astor BC, Woodward M, Levey AS, Jong PE, Coresh J, de Jong PE, El-Nahas M, et al: Lower estimated GFR and higher albuminuria are associated with adverse kidney outcomes in both general and high-risk populations. A collaborative meta-analysis of general and high-risk population cohorts. Kidney Int. 2011, 80 (1): 93-104. 10.1038/ki.2010.531.
Reich HN, Gladman DD, Urowitz MB, Bargman JM, Hladunewich MA, Lou W, Fan SC, Su J, Herzenberg AM, Cattran DC, et al: Persistent proteinuria and dyslipidemia increase the risk of progressive chronic kidney disease in lupus erythematosus. Kidney Int. 2011, 9 (8): 914-920.
Rao M, Wong C, Kanetsky P, Girndt M, Stenvinkel P, Reilly M, Raj DS: Cytokine gene polymorphism and progression of renal and cardiovascular diseases. Kidney Int. 2007, 72 (5): 549-556. 10.1038/sj.ki.5002391.
Munshi R, Hsu C, Himmelfarb J: Advances in understanding ischemic acute kidney injury. BMC Med. 2011, 9 (1): 11-10.1186/1741-7015-9-11.
Kottgen A, Pattaro C, Boger CA, Fuchsberger C, Olden M, Glazer NL, Parsa A, Gao X, Yang Q, Smith AV, et al: New loci associated with kidney function and chronic kidney disease. Nat Genet. 2010, 42 (5): 376-384. 10.1038/ng.568.
Chambers JC, Zhang W, Lord GM, van der Harst P, Lawlor DA, Sehmi JS, Gale DP, Wass MN, Ahmadi KR, Bakker SJ, et al: Genetic loci influencing kidney function and chronic kidney disease. Nat Genet. 2010, 42 (5): 373-375. 10.1038/ng.566.
Ribases M, Ramos-Quiroga JA, Sanchez-Mora C, Bosch R, Richarte V, Palomar G, Gastaminza X, Bielsa A, Arcos-Burgos M, Muenke M, et al: Contribution of LPHN3 to the genetic susceptibility to ADHD in adulthood: a replication study. Genes Brain Behav. 2010, 10 (2): 149-157.
Sole X, Guino E, Valls J, Iniesta R, Moreno V: SNPStats: a web tool for the analysis of association studies. Bioinformatics. 2006, 22 (15): 1928-1929. 10.1093/bioinformatics/btl268.
Fried L, Solomon C, Shlipak M, Seliger S, Stehman-Breen C, Bleyer AJ, Chaves P, Furberg C, Kuller L, Newman A: Inflammatory and prothrombotic markers and the progression of renal disease in elderly individuals. J Am Soc Nephrol. 2004, 15 (12): 3184-3191. 10.1097/01.ASN.0000146422.45434.35.
Wolkow PP, Niewczas MA, Perkins B, Ficociello LH, Lipinski B, Warram JH, Krolewski AS: Association of urinary inflammatory markers and renal decline in microalbuminuric type 1 diabetics. J Am Soc Nephrol. 2008, 19 (4): 789-797. 10.1681/ASN.2007050556.
Nakamura E, Megumi Y, Kobayashi T, Kamoto T, Ishitoya S, Terachi T, Tachibana M, Matsushiro H, Habuchi T, Kakehi Y, et al: Genetic polymorphisms of the interleukin-4 receptor alpha gene are associated with an increasing risk and a poor prognosis of sporadic renal cell carcinoma in a Japanese population. Clin Cancer Res. 2002, 8 (8): 2620-2625.
Burgos PI, Causey ZL, Tamhane A, Kelley JM, Brown EE, Hughes LB, Danila MI, van Everdingen A, Conn DL, Jonas BL, et al: Association of IL4R single-nucleotide polymorphisms with rheumatoid nodules in African Americans with rheumatoid arthritis. Arthritis Res Ther. 2010, 12 (3): R75-10.1186/ar2994.
Tachdjian R, Mathias C, Al Khatib S, Bryce PJ, Kim HS, Blaeser F, O'Connor BD, Rzymkiewicz D, Chen A, Holtzman MJ, et al: Pathogenicity of a disease-associated human IL-4 receptor allele in experimental asthma. J Exp Med. 2009, 206 (10): 2191-2204. 10.1084/jem.20091480.
Zheng G, Wang Y, Xiang SH, Tay YC, Wu H, Watson D, Coombes J, Rangan GK, Alexander SI, Harris DC: DNA vaccination with CCL2 DNA modified by the addition of an adjuvant epitope protects against "nonimmune" toxic renal injury. J Am Soc Nephrol. 2006, 17 (2): 465-474. 10.1681/ASN.2005020164.
Kang YS, Lee MH, Song HK, Ko GJ, Kwon OS, Lim TK, Kim SH, Han SY, Han KH, Lee JE, et al: CCR2 antagonism improves insulin resistance, lipid metabolism, and diabetic nephropathy in type 2 diabetic mice. Kidney Int. 2010, 78 (9): 883-894. 10.1038/ki.2010.263.
Dai R, Ahmed SA: MicroRNA, a new paradigm for understanding immunoregulation, inflammation, and autoimmune diseases. Transl Res. 2011, 157 (4): 163-179. 10.1016/j.trsl.2011.01.007.
Messeguer X, Escudero R, Farre D, Nunez O, Martinez J, Alba MM: PROMO: detection of known transcription regulatory elements using species-tailored searches. Bioinformatics. 2002, 18 (2): 333-334. 10.1093/bioinformatics/18.2.333.
Farre D, Roset R, Huerta M, Adsuara JE, Rosello L, Alba MM, Messeguer X: Identification of patterns in biological sequences at the ALGGEN server: PROMO and MALGEN. Nucleic Acids Res. 2003, 31 (13): 3651-3653. 10.1093/nar/gkg605.
Wei L, Vahedi G, Sun HW, Watford WT, Takatori H, Ramos HL, Takahashi H, Liang J, Gutierrez-Cruz G, Zang C, et al: Discrete roles of STAT4 and STAT6 transcription factors in tuning epigenetic modifications and transcription during T helper cell differentiation. Immunity. 2010, 32 (6): 840-851. 10.1016/j.immuni.2010.06.003.
Nakayama T, Sato W, Kosugi T, Zhang L, Campbell-Thompson M, Yoshimura A, Croker BP, Johnson RJ, Nakagawa T: Endothelial injury due to eNOS deficiency accelerates the progression of chronic renal disease in the mouse. Am J Physiol Renal Physiol. 2009, 296 (2): F317-327.
Webber JL, Tooze SA: New insights into the function of Atg9. FEBS Lett. 2010, 584 (7): 1319-1326. 10.1016/j.febslet.2010.01.020.
Kullo IJ, Greene MT, Boerwinkle E, Chu J, Turner ST, Kardia SL: Association of polymorphisms in NOS3 with the ankle-brachial index in hypertensive adults. Atherosclerosis. 2008, 196 (2): 905-912. 10.1016/j.atherosclerosis.2007.02.008.
Popov AF, Hinz J, Schulz EG, Schmitto JD, Wiese CH, Quintel M, Seipelt R, Schoendube FA: The eNOS 786C/T polymorphism in cardiac surgical patients with cardiopulmonary bypass is associated with renal dysfunction. Eur J Cardiothorac Surg. 2009, 36 (4): 651-656. 10.1016/j.ejcts.2009.04.049.
Wang CH, Li F, Hiller S, Kim HS, Maeda N, Smithies O, Takahashi N: A modest decrease in endothelial NOS in mice comparable to that associated with human NOS3 variants exacerbates diabetic nephropathy. Proc Natl Acad Sci U S A. 2011, 108 (5): 2070-2075. 10.1073/pnas.1018766108.
Desmet FO, Hamroun D, Lalande M, Collod-Beroud G, Claustres M, Beroud C: Human Splicing Finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 2009, 37 (9): e67-10.1093/nar/gkp215.
Sironi M, Menozzi G, Riva L, Cagliani R, Comi GP, Bresolin N, Giorda R, Pozzoli U: Silencer elements as possible inhibitors of pseudoexon splicing. Nucleic Acids Res. 2004, 32 (5): 1783-1791. 10.1093/nar/gkh341.
Perneger TV: What's wrong with Bonferroni adjustments. BMJ. 1998, 316 (7139): 1236-1238. 10.1136/bmj.316.7139.1236.
Sterne JA, Davey Smith G: Sifting the evidence-what's wrong with significance tests?. BMJ. 2001, 322 (7280): 226-231. 10.1136/bmj.322.7280.226.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/13/58/prepub
We thank María Jesús Pérez (Hospital Clinico Universitario de Valladolid) who kindly collected the samples.
We thank Nick Weber who provided writing assistance for this manuscript.
This work was supported by grants from Instituto de Salud Carlos III (Ref: PI08/0738 and PI11/00245) to SR and Junta de Castilla y León (Ref: GRS 234/A/08) to ET. MAJS is supported by a grant from Instituto de Salud Carlos III (CM10/00105).
None of the authors has any potential financial conflict of interest related to this manuscript.
Study concept and design: ET, SR. Administrative, technical, or material support: EL, ET, MH, JIGH, JB, JFBM. Acquisition of data: EL, ET, MH, JIGH, JB. Statistical analysis and interpretation of data: MAJS, AFR, LSR, PF, JMGG, SR. Drafting of the manuscript: MAJS, AFR, SR. Critical revision of the manuscript for important intellectual content: PF, JFBM. Study supervision: SR. All authors read and approved the final manuscript.