- Open Access
Identification of a haplotype block in the 5q31 cytokine gene cluster associated with the susceptibility to severe malaria
Malaria Journalvolume 8, Article number: 232 (2009)
It has been previously demonstrated that a single nucleotide polymorphism (SNP) in the IL13 promoter region, IL13 -1055T>C (rs1800925), was associated with susceptibility to severe malaria in Thais. In the present study, fine association mapping for a cytokine gene cluster including IL4, IL5, and IL13 on chromosome 5q31 was conducted using the same malaria subjects to refine the region containing a primary variant or a haplotype susceptible to severe malaria.
A total of 82 SNPs spanning 522 kb of the 5q31 region were analysed in 368 patients with Plasmodium falciparum malaria (203 mild malaria and 165 severe malaria patients).
Only rs1881457 located in the promoter region of IL13, which is in linkage disequilibrium with rs1800925 (r2 = 0.73), showed a significant association with severe malaria after adjusting for multiple testing (P = 0.046 by permutation test). This SNP was in a haplotype block spanning 97 kb (from rs2069812 to rs2240032). The detected haplotype block contained the RAD50 gene and the promoter of IL13, but not the other genes.
A haplotype block in which a primary polymorphism associated with severe malaria is likely to be encoded was identified in Thai malaria patients.
Over the course of the last decade a number of studies have provided evidence for a linkage between the blood infection level of Plasmodium falciparum and the human chromosome 5q31 region in African populations [1–4]. In addition to malaria, the 5q31 region shows a linkage to the response against other infectious diseases such as schistosomiasis  and leishmaniasis . The 5q31-33 region contains genes encoding the T helper 2-type cytokines (the interleukin genes IL3, IL4, IL5, IL9, and IL13) and other immunologically active genes such as interferon regulatory factor-1 (IRF1). These genes are strong candidates for controlling the outcome of malaria infection.
In a previous study, three single nucleotide polymorphisms (SNPs) in the promoter regions of IL3, IL4, and IL13 were investigated. Of which, a SNP in the IL13 promoter region, IL13 -1055T>C (rs1800925), was found to be associated with susceptibility to severe malaria in Thais . However, a number of candidate genes or polymorphisms still remain to be analyzed. In addition, no other polymorphisms surrounding rs1800925 were analyzed and thus the possibility that the association of rs1800925 with severe malaria may have resulted from linkage disequilibrium (LD) from other polymorphisms could not be excluded. The aim of this study is to better define the genomic region showing the association with severe malaria on the 5q31 region.
A total of 368 adult patients with P. falciparum malaria (165 patients with severe P. falciparum malaria and 203 patients with mild malaria) living in northwest Thailand were enrolled in this study. All patients underwent treatment at the Hospital for Tropical Diseases, Faculty of Tropical Medicine, Mahidol University. Malarial infection by P. falciparum was confirmed by a positive blood smear for the asexual form of P. falciparum. Clinical manifestations of severe and mild malaria were classified according to the following definitions and criteria. A patient was classified as severe malaria when he/she has one or more of the following signs: high parasitaemia (>100,000 parasite/ml), hypoglycaemia (glucose <22 nmol/l), severe anaemia (haematocrit <20% or haemoglobin <7.0 g/dl), and a serum creatinine level of more than 3.0 mg/dl. In the present study, patients with cerebral malaria were not analyzed. Mild malaria was characterized by fever without other any underlying causes of infections and no manifestations of severe malaria as described above. All individuals were 13 years of age or older, and the mean ages of patients with mild malaria and those with severe malaria were 25.3 and 23.8, respectively. This study was approved by the institutional review board of the Faculty of Tropical Medicine, Mahidol University, and the Research Ethics Committees of the Faculty of Medicine, The University of Tokyo, and the Graduate School of Comprehensive Human Sciences, University of Tsukuba. Informed consent was obtained from all participants.
Genomic DNA was extracted from peripheral blood leukocytes using a QIAamp Blood Kit (Qiagen, Hilden, Germany). A total of 82 SNPs within a 522 kb region on human chromosome 5q31 were genotyped by using the DigiTag2 assay  or TaqMan assay (Table 1). These SNPs were selected to capture the LD structure on 5q31 in Asian populations .
The allele frequency at each SNP locus was compared between mild and severe malaria patients using the chi-square test and a permutation P value was calculated from 100000 permutations. A permutation P value of less than 0.05 was considered to be statistically significant. The pairwise LD coefficients (r2) between SNPs were calculated to evaluate the structure of LD on 5q31-33 in 368 Thai malaria patients. The frequencies of haplotypes consisting of rs2069812, rs2299015, rs2299014, rs2243677, rs2522414, rs2299013, rs2252775, rs2522394, rs2245460, rs2301713, rs3798135, rs2237060, rs2074369, rs2240032, rs1881457, and rs1800925 were estimated only for this haplotype block. All the statistical analyses were performed by using the Haploview software version 4.0 . The allelic state for each SNP (i.e., ancestral or derived) was inferred based on the genome sequence of Pan troglodytes (chimpanzee), obtained from the NCBI BLAST database (Table 1). When the genomic sequence of chimpanzee was not available, one of Macaca mulatta (rhesus macaque) was used.
Eighty-two SNPs including rs1800925 were analysed to evaluate the association of the 5q31 region with severity of malaria (Table 1). The permutation P value as well as the raw P value was calculated for each SNP to avoid any false positive findings due to multiple testing (Table 1 and Figure 1A). Only rs1881457 showed a significant association with severe malaria (raw P value = 0.002 and permutation P value = 0.046) and no SNPs in the other genes showed any such association (permutation P value > 0.05). When a derived allele is focused on in association test, rs1881457-C may be referred to as a protective allele against severe malaria.
In previous study, rs1800925 was found to be associated with severe malaria . Since rs1881457, showing the strongest association in the present study, is closely located to rs1800925, these SNPs may be in LD. In addition, a number of SNPs near rs1881457 and rs1800925 showed also raw P values of less than 0.05 (Table 1 and Figure 1), thus suggesting that some, if not all, of these SNPs are in the same haplotype block. To clarify the structure of the LD around rs1881457 and rs1800925, r2 values between the 82 SNPs were calculated. The LD analysis for the 5q31 region revealed that all the SNPs showing low P values were in a distinct haplotype block spanning 97 kb from rs2069812 to rs2240032 (Figure 1B). This block contained the RAD50 gene and the promoter of IL13, but none of the other candidate genes such as IL4, IL5, and IRF1.
Six frequent haplotypes were observed in the detected block and two of which, haplotypes 1 and 4, bore rs1881457-C (Table 2). Both haplotypes showed a decreased frequency in severe malaria patients in comparison to those with mild malaria, thus suggesting that the association of rs1881457-C with the protection against severe malaria was not caused by a specific haplotype.
In this study rs1881457 was found to be significantly associated with severe malaria, and this SNP was included in a haplotype block encompassing the whole RAD50 gene and the promoter of IL13 (Figure 1). Together with MRE11 and NBS1, RAD50 forms a conserved multiprotein complex, MRE11-RAD50-NBS1 (MRN), which plays an important role in double-strand break repair, cell cycle checkpoint control, meiotic recombination, and telomere maintenance [11–13]. In the immune system, the MRN complex is involved in B cell-specific immunoglobulin gene diversification (e.g., Ig class-switch recombination, somatic hypermutation, and gene conversion) [14, 15]. The polymorphisms of RAD50 may therefore influence the affinity and/or effector functions of antibodies. The IL13 gene encodes a immunoregulatory cytokine (Th2 cytokine) produced by activated Th2 cells. The Th2 cytokines down-regulate macrophage activity, and inhibit the production of pro-inflammatory cytokines such as TNF and IL1. It has been reported that increased concentrations of TNF and IL1β in serum are observed in severe malaria patients . The IL13 promoter polymorphisms may influence the expression of IL13. Thus, both RAD50 and IL13 seem to be plausible candidate genes associated with severe malaria.
The genes encoding the Th2 cytokines IL5, IL13, and IL4 are subject to coordinate regulation and are expressed in a cell lineage-specific manner [17, 18]. The expressions are regulated by a locus control region (LCR) located within a 25 kb region containing the 3' portion of RAD50 . Interestingly, the LCR is included in the haplotype block associated with severe malaria, raising a possibility that polymorphisms which influence the LCR activity could account for the observed association with the severity of malaria.
Since rs1881457 is located in the promoter region of IL13, the nucleotide change at this site may affect the binding ability of some transcription factor. The TFSEARCH (TFSEARCH: Searching Transcription Factor Binding Sites, http://www.rwcp.or.jp/papia/) program based on the TRANSFAC databases  was used to examine the possibility. The result indicated no possible binding site of transcription factor at rs1881457 regardless of alleles (rs1881457-A and rs1881457-C) with the default setting (threshold score = 0.85). Therefore, rs1881457 itself may not be a primary polymorphism associated with severe malaria, even though rs1881457 showed the strongest association observed in this study.
Among IL13 polymorphisms, rs1800925 in the IL13 promoter has been reported to be associated with various diseases [21–24]. This SNP is located within a putative primate-specific cis-regulatory element  and has been shown to affect the promoter activity of IL13  and IL13 production . In the present study rs1800925 and rs1881457 had a high r2 value (r2 = 0.73). Therefore, the possibility that rs1800925 is a primary SNP and the significant association of rs1881457 with severe malaria resulted from LD between these SNPs is not excluded. The future functional and association studies for rs1881457, rs1800925, and other polymorphisms, including those not analyzed in the present study, may thus help us to better understand the genetic susceptibility to severe malaria.
A haplotype block spanning 97 kb encompassing RAD50 gene and IL13 promoter region that was associated with severity of malaria was identified in a Thai population.
Flori L, Kumulungui B, Aucan C, Esnault C, Traore AS, Fumoux F, Rihet P: Linkage and association between Plasmodium falciparum blood infection levels and chromosome 5q31-q33. Genes Immun. 2003, 4: 265-268. 10.1038/sj.gene.6363960.
Garcia A, Marquet S, Bucheton B, Hillaire D, Cot M, Fievet N, Dessein AJ, Abel L: Linkage analysis of blood Plasmodium falciparum levels: interest of the 5q31-q33 chromosome region. Am J Trop Med Hyg. 1998, 58: 705-709.
Rihet P, Traore Y, Abel L, Aucan C, Traore-Leroux T, Fumoux F: Malaria in humans: Plasmodium falciparum blood infection levels are linked to chromosome 5q31-q33. Am J Hum Genet. 1998, 63: 498-505. 10.1086/301967.
Sakuntabhai A, Ndiaye R, Casademont I, Peerapittayamongkol C, Rogier C, Tortevoye P, Tall A, Paul R, Turbpaiboon C, Phimpraphi W, Trape JF, Spiegel A, Heath S, Mercereau-Puijalon O, Dieye A, Julier C: Genetic determination and linkage mapping of Plasmodium falciparum malaria related traits in Senegal. PLoS ONE. 2008, 3: e2000-10.1371/journal.pone.0002000.
Marquet S, Abel L, Hillaire D, Dessein H, Kalil J, Feingold J, Weissenbach J, Dessein AJ: Genetic localization of a locus controlling the intensity of infection by Schistosoma mansoni on chromosome 5q31-q33. Nat Genet. 1996, 14: 181-184. 10.1038/ng1096-181.
Jeronimo SM, Holst AK, Jamieson SE, Francis R, Martins DR, Bezerra FL, Ettinger NA, Nascimento ET, Monteiro GR, Lacerda HG, Miller EN, Cordell HJ, Duggal P, Beaty TH, Blackwell JM, Wilson ME: Genes at human chromosome 5q31.1 regulate delayed-type hypersensitivity responses associated with Leishmania chagasi infection. Genes Immun. 2007, 8: 539-551. 10.1038/sj.gene.6364422.
Ohashi J, Naka I, Patarapotikul J, Hananantachai H, Looareesuwan S, Tokunaga K: A single-nucleotide substitution from C to T at position -1055 in the IL-13 promoter is associated with protection from severe malaria in Thailand. Genes Immun. 2003, 4: 528-531. 10.1038/sj.gene.6364010.
Nishida N, Tanabe T, Takasu M, Suyama A, Tokunaga K: Further development of multiplex single nucleotide polymorphism typing method, the DigiTag2 assay. Anal Biochem. 2007, 364: 78-85. 10.1016/j.ab.2007.02.005.
Nuchnoi P, Ohashi J, Naka I, Nacapunchai D, Tokunaga K, Nishida N, Patarapotikul J: Linkage disequilibrium structure of the 5q31-33 region in a Thai population. J Hum Genet. 2008, 53: 850-856. 10.1007/s10038-008-0309-8.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.
D'Amours D, Jackson SP: The Mre11 complex: at the crossroads of dna repair and checkpoint signalling. Nat Rev Mol Cell Biol. 2002, 3: 317-327. 10.1038/nrm805.
Petrini JH, Stracker TH: The cellular response to DNA double-strand breaks: defining the sensors and mediators. Trends Cell Biol. 2003, 13: 458-462. 10.1016/S0962-8924(03)00170-3.
Bosch van den M, Bree RT, Lowndes NF: The MRN complex: coordinating and mediating the response to broken chromosomes. EMBO Rep. 2003, 4: 844-849. 10.1038/sj.embor.embor925.
Kracker S, Bergmann Y, Demuth I, Frappart PO, Hildebrand G, Christine R, Wang ZQ, Sperling K, Digweed M, Radbruch A: Nibrin functions in Ig class-switch recombination. Proc Natl Acad Sci USA. 2005, 102: 1584-1589. 10.1073/pnas.0409191102.
Yabuki M, Fujii MM, Maizels N: The MRE11-RAD50-NBS1 complex accelerates somatic hypermutation and gene conversion of immunoglobulin variable regions. Nat Immunol. 2005, 6: 730-736. 10.1038/ni1215.
Kwiatkowski D, Hill AV, Sambou I, Twumasi P, Castracane J, Manogue KR, Cerami A, Brewster DR, Greenwood BM: TNF concentration in fatal cerebral, non-fatal cerebral, and uncomplicated Plasmodium falciparum malaria. Lancet. 1990, 336: 1201-1204. 10.1016/0140-6736(90)92827-5.
Kelly BL, Locksley RM: Coordinate regulation of the IL-4, IL-13, and IL-5 cytokine cluster in Th2 clones revealed by allelic expression patterns. J Immunol. 2000, 165: 2982-2986.
Smale ST, Fisher AG: Chromatin structure and gene regulation in the immune system. Annu Rev Immunol. 2002, 20: 427-462. 10.1146/annurev.immunol.20.100301.064739.
Lee GR, Fields PE, Griffin TJ, Flavell RA: Regulation of the Th2 cytokine locus by a locus control region. Immunity. 2003, 19: 145-153. 10.1016/S1074-7613(03)00179-1.
Heinemeyer T, Wingender E, Reuter I, Hermjakob H, Kel AE, Kel OV, Ignatieva EV, Ananko EA, Podkolodnaya OA, Kolpakov FA, Podkolodny NL, Kolchanov NA: Databases on transcriptional regulation: TRANSFAC, TRRD and COMPEL. Nucleic Acids Res. 1998, 26: 362-367. 10.1093/nar/26.1.362.
Chang M, Li Y, Yan C, Callis-Duffin KP, Matsunami N, Garcia VE, Cargill M, Civello D, Bui N, Catanese JJ, Leppert MF, Krueger GG, Begovich AB, Schrodi SJ: Variants in the 5q31 cytokine gene cluster are associated with psoriasis. Genes Immun. 2008, 9: 176-181. 10.1038/sj.gene.6364451.
Black S, Teixeira AS, Loh AX, Vinall L, Holloway JW, Hardy R, Swallow DM: Contribution of functional variation in the IL13 gene to allergy, hay fever and asthma in the NSHD longitudinal 1946 birth cohort. Allergy. 2009, 64: 1172-1178. 10.1111/j.1398-9995.2009.01988.x.
Nedoszytko B, Niedoszytko M, Lange M, van Doormaal J, Glen J, Zablotna M, Renke J, Vales A, Buljubasic F, Jassem E, Roszkiewicz J, Valent P: Interleukin-13 promoter gene polymorphism -1112C/T is associated with the systemic form of mastocytosis. Allergy. 2009, 64: 287-294. 10.1111/j.1398-9995.2008.01827.x.
Llanes E, Quiralte J, Lopez E, Sastre B, Chacartegui M, del Pozo V, Palomino P, Lahoz C, Cardaba B: Analysis of polymorphisms in olive pollen allergy: IL13, IL4RA, IL5 and ADRB2 genes. Int Arch Allergy Immunol. 2009, 148: 228-238. 10.1159/000161583.
Cameron L, Webster RB, Strempel JM, Kiesler P, Kabesch M, Ramachandran H, Yu L, Stern DA, Graves PE, Lohman IC, Wright AL, Halonen M, Klimecki WT, Vercelli D: Th2 cell-selective enhancement of human IL13 transcription by IL13-1112C>T, a polymorphism associated with allergic inflammation. J Immunol. 2006, 177: 8633-8642.
van der Pouw Kraan TC, van Veen A, Boeije LC, van Tuyl SA, de Groot ER, Stapel SO, Bakker A, Verweij CL, Aarden LA, van der Zee JS: An IL-13 promoter polymorphism associated with increased risk of allergic asthma. Genes Immun. 1999, 1: 61-65. 10.1038/sj.gene.6363630.
The authors sincerely thank all patients who kindly participated in this study. This study was supported in part by research funds KAKENHI Grant-in-Aid for Scientific Research on Priority Areas "Comprehensive Genomics" and "Applied Genomics" from the Ministry of Education, Culture, Sports, Science and Technology of Japan.
The authors declare that they have no competing interests.
NI carried out the genotyping, helped to conduct statistical analyses, and wrote the manuscript. NN and KT helped to genotype the samples. JP, PN and HH collected blood samples, extracted DNA, and helped to genotype the samples. JP participated in the design of the study and coordination. NT was involved in the interpretation of the data and preparation of the manuscript. JO conceived of the study, and participated in its design, performed statistical analyses, and helped to write the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.