Genetic evidence of regulatory gene variants of the STAT6, IL10R and FOXP3 locus as a susceptibility factor in uncomplicated malaria and parasitaemia in Congolese children
© Koukouikila-Koussounda et al.; licensee BioMed Central Ltd. 2013
Received: 26 July 2012
Accepted: 30 November 2012
Published: 8 January 2013
Regulatory T cells (Tregs) are a subset of T cells that play an important role in modulating T effector responses during infectious challenges. The aim of this study was to evaluate possible associations between regulatory gene polymorphisms and the risk of uncomplicated malaria and the control of Plasmodium falciparum parasite density levels.
Twelve regulatory single nucleotide polymorphisms (SNPs) in the promoter regions of FOXP3 (ss270137548, rs11091253), IL10RA (rs56356146, rs7925112), IL10RB (rs8178433, rs8178435, rs999788), STAT6 (rs3024941, rs3024943, rs3024944) and TNFRSF18 (ss2080581728, rs3753344) were genotyped in a cohort of Congolese children. Studied subjects were followed up (passively) during one year. The children who experienced one or several clinical episodes were genotyped as “uncomplicated malaria” group (n=179) and those children who did not experience any episode were genotyped as “asymptomatic children” group (n=138).
The prevalence of rs3024944CC genotype of STAT6 was significantly higher in the group of asymptomatic children compared to that of uncomplicated malaria (P=0.003). Similarly, the minor allele rs3024944C was more prevalent in the group of asymptomatic children (P=0.019). Two novel SNPs were observed including -163T/G (ss491228441) in IL10RA gene and -163C/T (ss491228440) in TNFRSF18 gene. The genotype ss491228441TT and the minor allele ss491228441G of the IL10RA were more frequent in the group of asymptomatic children (P=0.006 and P=0.007, respectively). The genotype rs11091253CT of the FOXP3 was associated with high parasite density levels. In addition, a new promoter IL10RA variant (ss491228441) contributes to shield against mild malaria.
The study indicated that the STAT6 promoter polymorphism rs3024944 was associated with uncomplicated malaria, whereas the FOXP3 promoter variant rs11091253 was associated with significant P. falciparum parasitaemia levels. These genetic data may contribute to the understanding of molecular mechanisms that regulate immune response to P. falciparum infections.
KeywordsPlasmodium falciparum Tregs FOXP3 IL10RA STAT6
Malaria is one of the major causes of morbidity and mortality in tropical and sub-tropical areas. Approximately 216 million clinical cases of malaria occurred worldwide in 2010, resulting in almost 655,000 deaths . Sub-Saharan Africa remains the region with the highest burden of the disease, accounting for nearly 91% of global malaria deaths. Malaria is caused by protozoan parasites of the genus Plasmodium and Plasmodium falciparum being the most virulent species. Plasmodium falciparum-infected individuals show different clinical phenotypes, which range from asymptomatic infections to severe forms of malaria . Several investigations demonstrated that immune gene polymorphisms in human host, such as those having an impact on expression of T regulatory cells (Treg), contribute to the variability of malaria phenotypes [2, 3]. The generation of Tregs has been demonstrated to modulate the immune response in malaria infection as well as in other infectious diseases [4–8]. A few of the Tregs are induced in response to infectious challenge and the others are considered as natural regulators . Expression of Tregs is well known to be regulated by number of human host genes, such as the forkhead box P3 (FOXP3), interleukin 10 receptors alpha and beta (IL10RA and IL10RB), the signal transducer and activator of transcription 6 (STAT6), and the tumour necrosis factor receptor superfamily member 18 (TNFRSF18) .
FOXP3 belongs to the family of transcription factors that play a role in various cellular processes and is described as an important regulator of natural Tregs development and function [10–12]. In humans, mutations in the FOXP3 gene are known to be associated with immune dysregulation polyendocrinophathy and enterophathy X-linked syndrome (IPEX) caused by lack of Tregs . IL10RA and IL10RB are two receptor chains forming the human interleukin 10-receptor (IL10R) which belongs to the class II cytokine receptors . The major role of IL10RA precursor is to mediate high affinity and ligand binding and subsequent signal transduction, whereas the IL10RB precursor is believed to contribute only in signalling process . Reports have demonstrated that mutations in the genes encoding IL10RA and IL10RB precursors of the IL10R annul the IL10 mediated immunomodulatory signalling and are strongly associated with inflammatory bowel disorders . Recently combined promoter haplotypes of the IL10RA and IL10RB genes have been shown to be associated with protection against severe malaria in Gabonese children . The STAT6 gene, which belongs to the STAT family of transcription factors, has been reported to mediate various cytokine-induced responses. Recent reports have shown that STAT6 gene polymorphisms play crucial roles in controlling helminth infection. More precisely, a promoter variant rs324013 was found to be associated with the control of Schistosoma haematobium in a Dogon population (where S. haematobium is endemic) in Mali . Similarly, the variant rs324015 was associated to the control of infection levels in Chinese children from an endemic area for Ascaris lumbricoides. TNFRSF18, also known as glucocorticoid-induced TNF receptor, is a member of TNF superfamily, which plays a key role in the co-stimulatory process for T cell activation for both conventional and CD25+CD4+ T cells. In addition, TNFRSF18 expressed on Tregs was demonstrated to play a crucial role in the maintenance of peripheral tolerance . Therefore, all these findings suggest that mutations in FOXP3, IL10RA, IL10RB, STAT6 and TNFRSF18 may lead to various immunopathological reactions in auto-immune disorders and/or parasitic infections.
In a previous investigation, regulatory SNPs in the promoter region of FOXP3, IL10RA, IL10RB, STAT6 and TNFRSF18 genes were identified in the Gabonese population in Central Africa. These identified SNPs were further validated for their allelic gene expression using transient transfection assays [14, 19–21]. Against this background, in the present study, the possible relationship between the reported SNP variants and risk of having uncomplicated P. falciparum malaria was investigated. For this purpose, this study utilized a case-control cohort of Congolese children, who were followed up (passively) during one year for surveillance of clinical malaria cases, comprising a group of children who experienced one or several uncomplicated malaria episodes and that of children who had no clinical episode during the follow up. As a second objective, the impact of these SNPs carriage on the P.falciparum parasite density levels was assessed.
Study area, study population and ethical consideration
The study was carried out in three districts (Kinsana, Mbouono and Ngoko) of Makélékélé health division in the southern area of Brazzaville. These districts are semi-urban with about 6,000 inhabitants and located along the Congo River. Malaria is transmitted throughout the year with P. falciparum being the predominant species and Anopheles gambiae the main mosquito vector . This study, as part of a cohort study for malaria surveillance in the study area (Ndounga et al, unpublished), is a community-based case-control study. The study population consisted of children aged one to nine years permanently living in the study area. All the study individuals belong to Bantu ethnicity and had equivalent risk for infection. Children were recruited between April and June 2010 and passively followed up during one year. At the time of enrolment 4 ml of whole blood were collected from all children and a thick and thin blood smears were performed. The inclusion criteria were: absence of clinical malaria in the previous two weeks and one week after enrolment (an axillary temperature of <37.5°C, no disease symptoms such as muscular pain and headaches), negative thick blood smear or positive blood smear with less than 5,000 parasites per microlitre of blood (p/μl), and not being homozygous for sickle cell trait. During the follow up, at any time, when a child was suspected to have fever or any malaria-related symptoms, parents or guardians were invited to present the child at the health centre for clinical examination. In case of malaria (positive thick and thin blood smears and axillary temperature ≥37.5°C), blood sample was collected, parasite density determined and children were treated with artesunate-amodiaquine or artemether-lumefantrine. No data were available about the present or past history of helminth infection.
The ethical approval was given by the Institutional Ethics Committee for Research on Health Sciences of the Republic of Congo. Written informed consent was obtained from parents or guardians of children.
Thick and thin blood films were stained with 10% Giemsa for 15 min and read by two independent, competent microscopists to determine the parasite density and Plasmodium species. Asexual parasites were counted against 200 leucocytes and expressed as the number of asexual parasites per μl, assuming the leucocytes count of 8,000/μl of blood.
Genomic DNA was extracted from 200 μl of peripheral whole blood samples obtained from patients and asymptomatic children using the commercially available QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions and stored at -20°C until use.
Primers used for genotyping and sequencing
Primer sequences (Forward/reverse)
94°C(5mn), [94°C(15S), 61°C(1mn) 72°C(1mn)]x35, 72°C(10mn)
94°C(2mn), [94°C(30S) 61°C(30S),72°C(30S)]x35, 72°C(2mn)
94°C(2mn), [94°C(30S) 61°C(30S),72°C(1mn)]x35, 72°C(2mn)
94°C(2mn), [94°C(30S) 61°C(30S),72°C(1mn)]x30, 72°C(2mn)
94°C(2mn), [94°C(30S) 61°C(30S),72°C(1mn)]x35, 72°C(2mn)
In this study, children without any reported malaria episode during the follow-up period were considered as asymptomatic and were classified in a group as asymptomatic children (as they could possibly be parasite carriers). Similarly, those who experienced at least one clinical malaria attack were classified in another group as uncomplicated malaria. Clinical and parasitological data for each child who had uncomplicated malaria were obtained from the longitudinal cohort study for malaria surveillance (Ndounga et al, unpublished). Those data included number of malaria attacks or episodes, parasite density and axillary temperature values.
Statistical analysis was performed using binary logistic regression model with Intercooled Stata v. 9.1 (STATA Corporation, Texas, USA) to determine the differences in allele frequencies and genotype distributions between uncomplicated malaria and asymptomatic children adjusted for age and gender. The level of significance was set at P<0.05. In addition, genetic models such as allelic, dominant and recessive were examined to evaluate the contributions of the genotype and allele frequencies using STATA. Kruskal-Wallis tests were employed to analyse the relationship between genotype distributions and the number of episodes as well as the risk of parasite density levels in the group of children with uncomplicated malaria. Genotype or haplotype frequencies were determined by simple gene counting and by using the expectation-maximum (EM) algorithm. The significance of deviation from Hardy-Weinberg equilibrium was tested using the random-permutation procedure as implemented in the Arlequin v. 18.104.22.168 software . Linkage disequilibrium (LD) analysis was performed using Haploview V. 3.2 program.
Characteristics of the study participants
Baseline characteristics of recruited children
Mild malaria (n=179)
Mean age (year)
4.73 ± 2.51
3.84 ± 2.56
Sex ratio (Male/Female)
Parasite carriers (%)
Hb AA carries (%)
Hb SS carriers (%)
Contribution of regulatory gene polymorphisms to uncomplicated clinical malaria
Distribution of gene variants in children from the group of uncomplicated malaria and that of asymptomatic children
Unc malar (n=179)
Asym (n=138 )
IL10RA ss491228441 T/G
STAT6 rs3024944 G/C
Distribution of STAT6 and IL10RA haplotypes
Unc. Malaria n=358(%)
STAT6 (rs3024941/ rs3024943 /rs3024944)
IL10RA (rs56356146 /rs7925112 /ss491228441)
Contribution of regulatory gene polymorphisms to clinical episodes and parasite density
Malaria has exerted selective pressure on the human genome that has led to changes, such as haemoglobin S which confers some resistance against severe malaria . Understanding the mechanisms leading to increased susceptibility or resistance against the disease would assist in the elaboration of better tools to fight the disease. In this search, immune responses that are on the front of the defence are particularly investigated. Tregs are suggested to be central in the control of P. falciparum infection. Therefore, investigation of human polymorphisms in genes associated with Treg activity may underlie both susceptibility to infection and level of Treg expression.
The present study evaluated the contribution of promoter SNPs of FOXP3, IL10RA, IL10RB, STAT6 and TNFRSF18 genes in P. falciparum infection outcome among Congolese children. Earlier studies have identified some regulatory SNPs on the promoter loci of these five genes from 40 unrelated Gabonese individuals who had earlier infection episodes with various parasites, including P. falciparum[14, 19–21]. These SNPs were then validated for their allelic gene expression using transient transfection assays. Some of the identified were already pre-described, while others were newly discovered. All the SNPs selected for the current investigation were found to alter expression of their respective gene [14, 19–21]. Therefore any possible correlations of these SNP variants with resistance or susceptibility to mild malaria attacks and/or the control of P. falciparum infection levels were evaluated.
The results obtained in this study demonstrated that the STAT6 rs3024944C/C genotype was significantly more prevalent in asymptomatic children. This suggests the contribution of this genotype to protection of Congolese children against mild clinical malaria. The SNP rs3024944 is located in a transcription factor-binding site TCF-1 which is a member of family of genes with homology to high mobility group I (HMG) proteins . TCF-1, in part with β-catenin, plays a pivotal role in T cell activation and is specifically expressed in T lymphocytes [25, 26]. And the fact that in the uncomplicated malaria group, children with the minor allele rs3024944C have low parasite densities compared to those having the major allele rs3024944G, that speculates that the rs3024944C/C variant directs activation of T cells in response to P. falciparum infection, thus contributing to protection against the disease through parasite elimination. In haplotype level, the major haplotype CGG was found to be correlated with the marginal risk of uncomplicated malaria. This emphasizes the role of the SNP rs3024944 on susceptibly/resistance to uncomplicated malaria among Congolese children. This is the first study that shows the correlation between a promoter SNP STAT6 rs3024944 and protection against uncomplicated malaria. However, in previous investigations, the promoter SNP rs324013 was found to be associated with the control of S. haematobium infection levels in a Dogon population in Mali . Similarly, the variant rs324015 was associated to the control of A. lumbricoides infection levels in Chinese children .
Two novel promoter SNPs were identified including, ss491228441 and ss491228440 of IL10RA and TNFRSF18 respectively. In both cases the frequency of the minor allele was low. Meanwhile the minor allele ss491228441G of IL10RA appeared to contribute to protection against uncomplicated malaria. Interestingly, the haplotype CTG involving this minor allele ss491228441G was significantly more present in asymptomatic children. This provides an additional evidence of the involvement of this allele in protection against uncomplicated clinical malaria. However, larger association studies in different malaria endemic areas including children and adults as well as functional studies of this ss491228441 SNP variant would be of great interest to verify the current observations.
The data also showed the association between the carriage of the FOXP3 rs11091253 SNP variant and P. falciparum parasite density levels. In particular, the rs11091253CT genotype was found to be associated with high infection levels, whereas the homozygous for the minor allele, the rs11091253TT genotype, was associated with lower infection levels. This suggests that FOXP3 contributes to the modulation of the immune response against P. falciparum infection. The SNP rs11091253 is located in a putative binding site for the transcription factors nuclear factor kappa B (NF-κB) and c-Rel, which is also a member of NF-κB family. These transcription factors have been demonstrated to be critical for regulation and transcription of FOXP3[27, 28]. Therefore the differential binding properties of the promoter carrying the rs11091253CT genotype and that of the promoter carrying the rs11091253TT genotype may explain the associations observed in this study. One possible explanation may be that the binding of NF-κB and/or c-Rel with the promoter carrying the rs11091253C/T genotype results in upregulation of FOXP3 which in turn leads to increased number of fully functional Tregs. Consequently, parasite multiplication is facilitated by suppression of either P. falciparum-responsive effector cells or pro-inflammatory cytokine response. This is in line with what was reported by previous study showing that upregulation of FOXP3, and CD4+CD25+ regulatory cells correlates with more rapid parasite growth in human malaria infection .
To summarize, selected regulatory SNPs in the promoter region of FOXP3, IL10RA, IL10RB, STAT6 and TNFRSF18 genes have been investigated for possible association with uncomplicated malaria and high P. falciparum parasite density among Congolese children. These findings indicate that the STAT6 promoter polymorphism rs3024944 is highly associated with protection against uncomplicated malaria, whereas the FOXP3 promoter variant rs11091253 is more likely to be involved in the control of P. falciparum parasite density levels. In addition, a new SNP detected in the promoter region of IL10RA (ss491228441) contributes to protection against clinical malaria. These findings deserve to be further investigated in a larger population of children in the Republic of Congo but also in other malaria endemic areas.
We thank the team of CERVE for its contribution in the field study. We are grateful to Prof Pascal Rihet for his helpful advices on results analysis. We also thank Silvelia Grummes and Viola Galinat for the excellent technical help for sequencing. FKK was a recipient of a training grant from CANTAM network. This work has been supported through the CANTAM project. CANTAM (Central Africa Network on Tuberculosis, HIV/AIDs and Malaria) is a network of excellence supported by EDCTP and NACCAP.
- World Health Organization: World malaria report 2011. 2011, Geneva: World Health OrganizationGoogle Scholar
- Maizels RM: Parasite immunomodulation and polymorphisms of the immune system. J Biol. 2009, 8: 62-10.1186/jbiol166.PubMed CentralView ArticlePubMedGoogle Scholar
- Weatherall DJ, Clegg JB: Genetic variability in response to infection: malaria and after. Genes Immun. 2002, 3: 331-337. 10.1038/sj.gene.6363878.View ArticlePubMedGoogle Scholar
- Belkaid Y, Rouse BT: Natural regulatory T cells in infectious disease. Nat Immunol. 2005, 6: 353-360.View ArticlePubMedGoogle Scholar
- Hisaeda H, Maekawa Y, Iwakawa D, Okada H, Himeno K, Kishihara K, Tsukumo S, Yasutomo K: Escape of malaria parasites from host immunity requires CD4+ CD25+ regulatory T cells. Nat Med. 2004, 10: 29-30. 10.1038/nm975.View ArticlePubMedGoogle Scholar
- Sakaguchi S, Sakaguchi N: Regulatory T cells in immunologic self-tolerance and autoimmune disease. Int Rev Immunol. 2005, 24: 211-226. 10.1080/08830180590934976.View ArticlePubMedGoogle Scholar
- Walther M, Tongren JE, Andrews L, Korbel D, King E, Fletcher H, Andersen RF, Bejon P, Thompson F, Dunachie SJ, Edele F, de Souza JB, Sinden RE, Gilbert SC, Riley EM, Hill AV: Upregulation of TGF-beta, FOXP3, and CD4+CD25+ regulatory T cells correlates with more rapid parasite growth in human malaria infection. Immunity. 2005, 23: 287-296. 10.1016/j.immuni.2005.08.006.View ArticlePubMedGoogle Scholar
- Velavan TP, Ojurongbe O: Regulatory T cells and parasites. J Biomed Biotechnol. 2011, 2011: 520940-PubMed CentralView ArticlePubMedGoogle Scholar
- Leoratti FM, Farias L, Alves FP, Suarez-Mutis MC, Coura JR, Kalil J, Camargo EP, Moraes SL, Ramasawmy R: Variants in the toll-like receptor signaling pathway and clinical outcomes of malaria. J Infect Dis. 2008, 198: 772-780. 10.1086/590440.View ArticlePubMedGoogle Scholar
- Gavin MA, Torgerson TR, Houston E, DeRoos P, Ho WY, Stray-Pedersen A, Ocheltree EL, Greenberg PD, Ochs HD, Rudensky AY: Single-cell analysis of normal and FOXP3-mutant human T cells: FOXP3 expression without regulatory T cell development. Proc Natl Acad Sci USA. 2006, 103: 6659-6664. 10.1073/pnas.0509484103.PubMed CentralView ArticlePubMedGoogle Scholar
- Sakaguchi S, Wing K, Miyara M: Regulatory T cells - a brief history and perspective. Eur J Immunol. 2007, 37 (Suppl 1): S116-S123.View ArticlePubMedGoogle Scholar
- Kim HP, Leonard WJ: CREB/ATF-dependent T cell receptor-induced FoxP3 gene expression: a role for DNA methylation. J Exp Med. 2007, 204: 1543-1551.PubMed CentralPubMedGoogle Scholar
- Sakaguchi S, Powrie F: Emerging challenges in regulatory T cell function and biology. Science. 2007, 317: 627-629. 10.1126/science.1142331.View ArticlePubMedGoogle Scholar
- Velavan TP, Buyukyazici B, Kremsner PG, Kun JF: Combined promoter haplotypes of the IL10R genes are associated with protection against severe malaria in Gabonese children. Immunogenetics. 2012, 64: 87-95. 10.1007/s00251-011-0563-3.View ArticlePubMedGoogle Scholar
- Gasche C, Grundtner P, Zwirn P, Reinisch W, Shaw SH, Zdanov A, Sarma U, Williams LM, Foxwell BM, Gangl A: Novel variants of the IL-10 receptor 1 affect inhibition of monocyte TNF-alpha production. J Immunol. 2003, 170: 5578-5582.View ArticlePubMedGoogle Scholar
- Glocker EO, Kotlarz D, Boztug K, Gertz EM, Schaffer AA, Noyan F, Perro M, Diestelhorst J, Allroth A, Murugan D, Hätscher N, Pfeifer D, Sykora KW, Sauer M, Kreipe H, Lacher M, Nustede R, Woellner C, Baumann U, Salzer U, Koletzko S, Shah N, Segal AW, Sauerbrey A, Buderus S, Snapper SB, Grimbacher B, Klein C: Inflammatory bowel disease and mutations affecting the interleukin-10 receptor. N Engl J Med. 2009, 361: 2033-2045. 10.1056/NEJMoa0907206.PubMed CentralView ArticlePubMedGoogle Scholar
- He H, Isnard A, Kouriba B, Cabantous S, Dessein A, Doumbo O, Chevillard C: A STAT6 gene polymorphism is associated with high infection levels in urinary schistosomiasis. Genes Immun. 2008, 9: 195-206. 10.1038/gene.2008.2.View ArticlePubMedGoogle Scholar
- Peisong G, Yamasaki A, Mao XQ, Enomoto T, Feng Z, Gloria-Bottini F, Bottini E, Shirakawa T, Sun D, Hopkin JM: An asthma-associated genetic variant of STAT6 predicts low burden of ascaris worm infestation. Genes Immun. 2004, 5: 58-62. 10.1038/sj.gene.6364030.View ArticlePubMedGoogle Scholar
- Hanel SA, Velavan TP, Kremsner PG, Kun JF: Novel and functional regulatory SNPs in the promoter region of FOXP3 gene in a Gabonese population. Immunogenetics. 2011, 63: 409-415. 10.1007/s00251-011-0524-x.View ArticlePubMedGoogle Scholar
- Velavan TP, Bechlars S, Huang X, Kremsner PG, Kun JF: Novel regulatory SNPs in the promoter region of the TNFRSF18 gene in a Gabonese population. Braz J Med Biol Res. 2011, 44: 418-420. 10.1590/S0100-879X2011007500036.View ArticlePubMedGoogle Scholar
- Velavan TP, Bechlars S, Tomiuk J, Kremsner PG, Kun JF: Molecular characterization of regulatory polymorphisms in the promoter region of the STAT6 gene in a Gabonese population. Mem Inst Oswaldo Cruz. 2011, 106: 65-69. 10.1590/S0074-02762011000100011.View ArticlePubMedGoogle Scholar
- Excoffier L, Lischer HE: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010, 10: 564-567. 10.1111/j.1755-0998.2010.02847.x.View ArticlePubMedGoogle Scholar
- Fortin A, Stevenson MM, Gros P: Susceptibility to malaria as a complex trait: big pressure from a tiny creature. Hum Mol Genet. 2002, 11: 2469-2478. 10.1093/hmg/11.20.2469.View ArticlePubMedGoogle Scholar
- Castrop J, van WD, Koomans-Bitter M, van de Wetering M, de WR, van DJ, Clevers H: The human TCF-1 gene encodes a nuclear DNA-binding protein uniquely expressed in normal and neoplastic T-lineage lymphocytes. Blood. 1995, 86: 3050-3059.PubMedGoogle Scholar
- Staal FJ, Burgering BM, van de Wetering M, Clevers HC: Tcf-1-mediated transcription in T lymphocytes: differential role for glycogen synthase kinase-3 in fibroblasts and T cells. Int Immunol. 1999, 11: 317-323. 10.1093/intimm/11.3.317.View ArticlePubMedGoogle Scholar
- Yu Q, Sharma A, Sen JM: TCF1 and beta-catenin regulate T cell development and function. Immunol Res. 2010, 47: 45-55. 10.1007/s12026-009-8137-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Hori S: c-Rel: a pioneer in directing regulatory T-cell lineage commitment?. Eur J Immunol. 2010, 40: 664-667. 10.1002/eji.201040372.View ArticlePubMedGoogle Scholar
- Long M, Park SG, Strickland I, Hayden MS, Ghosh S: Nuclear factor-kappaB modulates regulatory T cell development by directly regulating expression of Foxp3 transcription factor. Immunity. 2009, 31: 921-931. 10.1016/j.immuni.2009.09.022.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.