Figures
Abstract
Background
Cerebral malaria (CM) represents a severe outcome of the Plasmodium falciparum infection. Recent genetic studies have correlated human genes with severe malaria susceptibility, but there is little data on genetic variants that increase the risk of developing specific malaria clinical complications. Nevertheless, susceptibility to experimental CM in the mouse has been linked to host genes including Transforming Growth Factor Beta 2 (TGFB2) and Heme oxygenase-1 (HMOX1). Here, we tested whether those genes were governing the risk of progressing to CM in patients with severe malaria syndromes.
Methodology/Principal Findings
We report that the clinical outcome of P. falciparum infection in a cohort of Angolan children (n = 430) correlated with nine TGFB2 SNPs that modify the risk of progression to CM as compared to other severe forms of malaria. This genetic effect was explained by two haplotypes harboring the CM-associated SNPs (Pcorrec. = 0.035 and 0.036). In addition, one HMOX1 haplotype composed of five CM-associated SNPs increased the risk of developing the CM syndrome (Pcorrec. = 0.002) and was under-transmitted to children with uncomplicated malaria (P = 0.036). Notably, the HMOX1-associated haplotype conferred increased HMOX1 mRNA expression in peripheral blood cells of CM patients (P = 0.012).
Conclusions/Significance
These results represent the first report on CM genetic risk factors in Angolan children and suggest the novel hypothesis that genetic variants of the TGFB2 and HMOX1 genes may contribute to confer a specific risk of developing the CM syndrome in patients with severe P. falciparum malaria. This work may provide motivation for future studies aiming to replicate our findings in larger populations and to confirm a role for these genes in determining the clinical course of malaria.
Citation: Sambo MR, Trovoada MJ, Benchimol C, Quinhentos V, Gonçalves L, Velosa R, et al. (2010) Transforming Growth Factor Beta 2 and Heme Oxygenase 1 Genes Are Risk Factors for the Cerebral Malaria Syndrome in Angolan Children. PLoS ONE 5(6): e11141. https://doi.org/10.1371/journal.pone.0011141
Editor: Mauricio Martins Rodrigues, Federal University of São Paulo, Brazil
Received: February 18, 2010; Accepted: May 21, 2010; Published: June 16, 2010
Copyright: © 2010 Sambo et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This project was supported by the Instituto Gulbenkian de Ciência. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Cerebral malaria (CM) is an acute life threatening complication of clinical malaria that afflicts a fraction of the P. falciparum-infected patients and it is perceived as a complex syndrome, rather than an isolated entity [1], [2], [3]. While potentially reversible, CM can develop into an unfettered neurovascular inflammatory syndrome leading to irreversible coma [3], [4]. Studies on malaria heritability indicate that about 25% of the risk to progress to severe malaria is determined by human genetic factors [5], [6]. The strong protection against severe malaria conferred by the HbS [7], G6PD [8], [9], [10] and the alpha-thalassemia mutations [11], [12], [13] illustrate the striking selection pressure of malaria on the human genome.
Genetic variation may in part determine the epidemiological patterns of susceptibility/resistance to malaria observed at the population level, whereby a fraction of the individuals infected with the malaria parasite would present parasitaemia and uncomplicated disease, a fraction of uncomplicated disease cases would progress to severe disease and a fraction of severe disease cases would develop CM syndrome [14]. Presumably, the genetic factors controlling malaria severity would affect the risk of disease progression from uncomplicated malaria, characterized by low levels of parasitaemia and fever, towards the different severe clinical complications and outcomes. Although many studies emphasize the role of genetic factors in controlling the severity of the malaria infection, in most instances the study design does not allow comparisons of CM with non-cerebral forms of severe malaria [15], [16], [17], [18] precluding the identification of genetic variants that specify the risk of developing CM.
We performed a genetic study to test the prior hypothesis that host factors influence specifically the risk of developing CM syndrome in the context of severe malaria. The study design aimed to identify candidate genetic variants associated to CM by comparing a group of CM patients with a group of severe malaria patients that did not exhibit CM. This hospital-based study was carried out in Luanda, a province of Angola that represents a mesoendemic stable region and where malaria is the first cause of morbidity and mortality [19]. Malaria transmission in Angola is perennial with enhancement during the rainy season and P. falciparum malaria is the most frequent infection (Angolan Ministry of Health, MINSA 2001).
The study of CM pathogenesis has benefited from investigations on experimental systems that use artificial laboratory infections to develop rodent models of experimental cerebral malaria (ECM). This report focused on four functional candidate genes (TGFB2, HMOX1, ICAM1 and CD36) that have been correlated with CM and other forms of severe malaria in mouse models. To our knowledge the TGFB2 gene has not been studied in the context of human malaria but it was suggested to confer resistance to ECM in a wild mouse strain. Heme oxygenase -1, encoded by Hmox1 gene, was shown to prevent the development of ECM but its role in human CM has not been addressed other than by an exploratory study in CM patients in Myanmar [20]. The ICAM 1 and CD36 genes have also been related to ECM in the mouse models [21], [22] but investigations on the role of these genes in human CM and severe malaria led to inconclusive and often contradictory results [15], [23], [24], [25], [26], [27], [28], [29].
This work represents an initial report on CM genetics in the Angolan population and though requiring further confirmation in larger samples from other populations, the results suggest that TGFB2 gene variants contribute to the risk of developing the CM syndrome when compared to other forms of severe malaria in Angolan children. The results also suggest that HMOX1 gene variants are associated to CM risk and control the HMOX1 mRNA expression in peripheral blood cells of children with CM.
Results
The sample population in this study consisted of 430 children diagnosed with malaria in the central pediatric hospital of Luanda (Angola) and the 319 uninfected apparently healthy controls, all attendees to the vaccination ward of the same hospital. A fraction of the CM patients (16.9%) showed severe malaria anemia and 11.5% also had hyperparasitaemia. Nevertheless, we did not find association between hyperparasitaemia and severe anemia (Chi-Square test, P value = 0.58) and 47.7% of the CM patients had low parasitaemia levels, suggesting that high levels of parasitaemia were not a requirement to develop CM. The ethnic group Kimbundu that predominates in Luanda was the most represented in our study (50.7%), while Umbundu was the second ethnic group (13.7%). Nevertheless, the ethnic composition did not differ among patient groups (Kruskal-Wallis Test, P-Value = 0.56) and there was no association between ethnic groups and malaria outcome (Chi-Square test, P-value = 0.74) suggesting that ethnic ancestry was not a major factor in stratifying the clinical malaria phenotypes.
All individuals were genotyped and analyzed for 54 single nucleotide polymorphisms (SNPs) that covered a region of 100Kb that encompassing the TGFB2 gene (Transforming growth factor, beta 2), 3.1 Kb in the CD36 gene (Thrombospondin receptor), 3.3 Kb in the HMOX1 (Heme oxygenase -1) gene region, 5.2 Kb in the ICAM 1 (Intercellular adhesion molecule 1, CD54) gene region as well as dbSNP rs334 which defines the HbS mutation. The complete list of SNPs that passed genotyping quality control and were subjected to association analysis is shown in Table S1. For the analyzed SNPs the Hardy-Weinberg (HWE) was preserved in all patient groups and uninfected apparently healthy controls, except for the HbS mutation that was in equilibrium in the healthy controls but not in the malaria patient groups.
To identify genetic polymorphisms specifically associated to CM, we used the strategy of comparing genotype frequencies between 130 patients with the CM syndrome and 158 patients with severe forms of malaria severe malaria in absence of CM (SnC). To better define CM-associated genetic factors that are exclusively detected when comparing CM to SnC patients, the study included a group of 142 patients with uncomplicated malaria (UM) and 319 uninfected healthy controls (UIF). The evaluation of the genetic association tests took into account that the sample size limited the resilience to multiple test corrections. Therefore, the initial single marker analysis was used to identify SNPs with suggestive effects on the CM phenotype and the subsequent analysis was pursued in the context of haplotype association using corrected P values (see methods section).
Single marker analysis
Genotype frequency analysis showed the frequency of the HbS mutation (rs334) did not differ significantly between CM and SnC patients but was notably higher in UM patients and UIF controls (Table 1). These results confirmed that the well-documented protective effect of the HbS mutation against severe malaria [7], [30] is measurable in the Angolan population and suggested that the sickle cell anemia trait is not specifically protecting from CM development.
To identify genetic polymorphisms that specifically modify the risk of developing CM, we searched for genotype differences between patients having CM and those having severe malaria in absence of CM (SnC). We found that the genotype frequencies in 9 out of the 26 SNPs analyzed in the TGFB2 gene and 3 out of 10 SNPs analyzed in HMOX1 gene showed suggestive association to CM (uncorrected P-value<0.05) when comparing CM patients with SnC patients but not in comparison with UM patients and UIF controls (Table 1). These results raised the possibility that variants of these genes could be involved in specifying the risk of developing CM in the context of severe malaria. Also, we noted that an additional SNP in the HMOX1 gene (rs2071748) showed significant association only when comparing CM with UM patients (Table 1). This finding suggested that the SNP rs2071748 could represent a risk factor for severe malaria, which was also detected when comparing SnC to UM patients (data not shown). This SNP maps very closely to the GT repeat polymorphism that has been associated to CM in a study in the Myanmar population [41]. We genotyped this microsatellite and found that the frequency of the GT-repeat short alleles tends to be higher in CM patients as compared to UIF controls (P = 0.0057) and to UM patients (P = 0.038) but not when compared to SnC patients (see Figure S1 and Table S2). This result suggests that the rs2071748 CM-association signal was also captured by the HMOX1-GT repeat. In contrast, two SNPs on the CD36 gene which was previously reported to be associated with disease severity [28], [31], [32], showed to be associated to CM when comparing to SnC patients but also with UM patients and UIF controls (Table 1). This result suggested that genetic variants of the CD36 gene play a role at various stages of disease progression but are not specifically involved in specifying the clinical outcome of the severe malaria. We also found that the SNPs of the ICAM 1 gene were not associated to CM (Table S3).
Linkage Disequilibrium analysis
We used linkage disequilibrium (LD) analysis to take into account the relative positions of the CM-associated SNPs in the TGFB2 and HMOX1 genes. Out of the 26 SNPs that were analyzed across 100 kb in the TGFB2 region, the nine SNPs specifically associated with CM clustered in a region spanning 48 Kb. LD analysis revealed a region of strong LD covering 27 kb of the structural TGFB2 gene that encompassed eight SNPs associated with CM. The remaining CM-associated SNP (rs1317681) maps 15 kb from this region and showed a weaker statistical association to CM (Table 1). The r2 among the eight SNPs varied from 0.6 to 1.0 suggesting that this region may represent a LD block in the TGFB2 gene (Figure 1). Similarly, we found that 5 out the 10 SNPs analyzed in the HMOX1 gene were mapping within a strong LD group (r2>0.9) that spanned approximately 8 kb. This region encompassed the three SNPs associated with CM, one SNP (rs9622194) marginally associated with CM (P-value = 0.06) (Table S3) and the SNP (rs2071748) associated with disease severity (Figure 2).
The diagram represents in the vertical axis, the log10 of the P value of association tests of SNPs listed in the y-axis, obtained from logistic regression analysis. The points in red, blue and green represent the results of the association with the SNPs between cases of cerebral malaria (CM) and controls with severe no cerebral malaria (SnC), uncomplicated malaria (UM) and uninfected (UIF), respectively. The upper diagram shows a scale representation of the structure of the TGFB2 gene: exons are represented by boxes and marked with its number. The bottom of the diagram represents the LD between a pair of SNPs measured by r2. The image, below and to the left, corresponds to the color scale representing the values of LD, from yellow (r2 = 0) to red (r2 = 1). This diagram is an adaptation of the figure produced by the snp. plotter R package.
The diagram represents in the vertical axis, the log10 of the P value of association tests of SNPs listed in the y-axis, obtained from logistic regression analysis. The points in red, blue and green represent the results of the association with the SNPs between cases of cerebral malaria (CM) and controls with severe no cerebral malaria (SnC), uncomplicated malaria (UM) and uninfected (UIF), respectively. The upper diagram shows a scale representation of the structure of the HMOX1 gene: exons are represented by boxes and marked with its number. The bottom of the diagram represents the LD between a pair of SNPs measured by r2. The image, below and to the left, corresponds to the color scale representing the values of LD, from yellow (r2 = 0) to red (r2 = 1). This diagram is an adaptation of the figure produced by the snp. plotter R package.
Haplotype association
The results of the LD analysis were used to reconstruct haplotypes (R package GAP) of the TGFB2 and HMOX1 genes and test whether the most frequent haplotypes (>5%) were associated to CM correcting statistical significance using permutation tests. Using the eight SNPs that showed exclusive CM-specific association within the TGFB2 LD block we reconstructed six haplotypes that were denoted TGFB2a to TGFB2e (Table 2). The analysis confirmed the overall association (Global Score Statistics) of the TGFB2 haplotypes when comparing CM with SnC patients (corrected P-value = 0.041). TGFB2c and TGFB2f haplotypes with distinct alleles in all seven CM-associated SNPs, showed opposite effects on the risk of developing CM (Table 2). The TGFB2c was significantly more frequent in CM patients than in those with other forms of severe malaria (corrected P-value = 0.035) and its genetic effect appears to be restricted to increase the risk of developing CM in the context of severe malaria. In contrast, the TGFB2f haplotype that was less frequent in the population significantly decreased the risk of progression to CM (corrected P value = 0.036). This protective effect is inferred from the reduction of the TGFB2f haplotype frequency among the CM patients that is still detectable when comparing CM patients to UIF controls (Table 2).
Haplotype reconstruction using the five SNPs in the LD block of the HOMX1 gene identified three haplotypes denoted HMOX1a, HMOX1b and HMOX1c (Table 3). Association analysis revealed that the HMOX1c haplotype frequency was significantly higher in CM patients as compared to SnC and UM patients as well as to uninfected controls (Table 3). These results suggest that the HMOX1c haplotype has a general role as a CM syndrome susceptibility factor particularly when compared to uncomplicated malaria patients (P = 0.002). On the other hand, the HMOX1b haplotype showed a CM protective effect with marginal significance when comparing CM with UM patients (corrected P-value = 0.056). In fact, the frequency of the HMOX1b haplotype was similar in CM patients and in patients with severe non-cerebral malaria. This finding suggested that this haplotype act at different stages of infection as a factor conferring protection from severe malaria (Table 3).
Transmission disequilibrium tests
We used transmission disequilibrium tests (TDT) in mother-affected child pairs to corroborate the haplotype associations observed in the case-control analysis. The TDT analysis was performed with the TRANSMIT software using 105 mother-affected child pairs with CM, 115 with SnC and 106 with UM. The TDT analysis detected the protective role of the HbS mutation as it was under-transmitted to CM children (P = 0.036) (Table 4). Moreover, we found that the TGFB2c haplotype (AGATGGTTA), deemed to be overrepresented in severe malaria patients that progressed to CM, was under-transmitted to SnC children that did not progress to cerebral malaria (P = 0.038) (Table 4). On the other hand, the HMOX1c haplotype (AGCGG) was found to be under-transmitted to UM children (P = 0.036), reinforcing the hypothesis that this haplotype represents a risk factor for severe malaria (Table 4). Overall, the TDT results corroborated the initial case-control association findings for haplotypes that showed relatively high frequency. Failure to confirm associations of less frequent haplotypes (e.g. TFB2f) by TDT is presumably due to reduced statistical power. Nevertheless, the results of the TDT analysis suggested that the associations of TGFB2 and HMOX1 genes to CM are not attributable to population stratification effects in groups of patients and controls.
HMOX1 SNPs control HMOX1 mRNA expression in children developing CM
To investigate the functional basis of HMOX1 association to malaria we examined whether the HMOX1 haplotypes associated to the clinical outcomes of malaria were also controlling the level of HMOX1 gene expression among children developing CM. Relative expression of HMOX1 mRNA in 42 children developing CM was quantified by Real-Time PCR in peripheral blood leukocytes. The basal levels were fairly high, as judged from the Real-Time PCR CTs with prominent curves. Quantitative trait locus (QTL) analysis detected suggestive association of the HMOX1c haplotype (AGCGG) to the level of HMOX1 mRNA expression (P≤0.01) (Table 5). This result revealed that the HMOX1c genetic variant that conferred increased risk to develop CM is also associated to higher expression of HMOX1 mRNA in the context of CM.
Discussion
In this work, we aimed to identify genetic variants specifying the emergence of CM in severe malaria patients by contrasting genotype frequencies in those patients presenting CM with those exhibiting signs of severe malaria but exempt of CM symptoms. Single marker analysis suggested that SNPs in the TGFB2 and HMOX1 genes showed specific association to the CM syndrome. In addition, haplotype association analysis showed that the TGFB2c haplotype specifically increased the risk of developing CM as compared to other severe malaria complications. This haplotype was also under-transmitted to patients with severe non-cerebral malaria. On the other hand, the HMOX1c haplotype seemed to play a general role in CM risk as it was overrepresented in CM patients as compared to SnC and UM patients as well as to uninfected controls. In accordance, this haplotype was under-transmitted to patients with uncomplicated malaria. Moreover, QTL analysis showed that the HMOX1c haplotype was associated to higher expression of HMOX1 mRNA in peripheral blood cells of CM patients.
We did not find that CD36 was a determinant of the risk of developing CM as compared to other severe malaria complications. In fact, the role of the CD36 gene in malaria is still debated and raises the question of whether different studies may capture different effects of the CD36 gene variants in severe and in cerebral malaria. We found that the SNP rs3211958 was not associated to CM although this SNP is known to be in complete LD with rs3211938 which showed contradictory results in different African populations [15], [26], [28], [33]. In our sample, the absence of CM-association with the ICAM 1 SNPs is in agreement with the results reported by a study in Gambia, Malawi and Kenya populations [26], though it has been reported that the heme oxygenase-1 has a down-regulating effect on the ICAM-1 expression levels [34], [35]. Such heterogeneity effects should be kept in mind when judging lack of replication of many reported SNP associations to disparate malaria clinical manifestations [36], [37], [38].
Though the pathogenesis and the symptoms of cerebral malaria in murine CM do not completely overlap the observations in human CM, the mouse models have brought critical insights to our understanding of CM pathogenesis [39], [40], [41]. Thus the TGFB2 gene was identified as a candidate gene to CM susceptibility in a genetic cross that used wild-derived mouse strains and where CM susceptibility segregated from hyperparasitaemia [42], suggesting that TGFB2 would be a gene conferring risk to CM in the context of severe malaria. Tgfb2−/− mice are not viable and studies with heterozygous mice in the context of malaria infection are warranted and may provide further evidence of involvement of the Tgfb2 gene in ECM pathogenesis.
Members of the TGF-beta cytokine family have been associated to the control of malaria infection and parasite growth [43], [44], [45] and TGFB isoforms were shown to accumulate in the brain of cerebral malaria patients [46]. The study of the TGFB2 haplotypic structure in the Angolan population led us to identify an LD block covering 27kb in the structural region of the gene. Although the CM-associated region in the TGFB2 gene covers part of the coding region, all the SNPs that we analyzed were intronic and have no predicted effects in folding, activity or function of the TGFB2 protein. Re-sequencing of relevant variants of the TGFB2 gene may provide further information on putative functional polymorphisms explanatory of the role of TGFB2 in CM pathogenesis. Studies of TGFB2 haplotypes in other populations and areas with different malaria endemicity would be required to determine whether the TGFB2 genetic variance represents a widespread CM risk factor. It should be noted that according to our observations the replication of TGFB2 specific association to CM in other populations is expected to be stronger when analyzed in the context of severe malaria.
The Heme oxygenase-1 activity has been shown to suppress ECM in mice [47] and one report in Myanmarese patients suggested that short alleles of a GT repeat in the HMOX1 promotor region are a risk factor for cerebral malaria [20]. Here, we show that one haplotype (HMOX1c) composed by SNPs covering a 8Kb LD block in the structural region of the HMOX1 gene confers susceptibility to CM which (Table 3). Strikingly, the occurrence of the allele A at position rs2071748 in the context of this haplotype generates the HMOX1b haplotype that marginally decreases the risk to progress from UM to CM. It is worth noting that the A at position rs2071748 is not able to mediate CM susceptibility in the context of the high frequency HMOX1a haplotype (Table 3). We also noted that the rs2071748 SNP was mapping in the HMOX1 promotor region very closely to the GT repeat polymorphism that showed that the short alleles were overrepresented in our CM patients and in Myanmar patients [20]. We also found that the HMOX1c haplotype governs the expression of HMOX1 mRNA in peripheral blood leucocytes (PBL) of CM patients. We speculate that HMOX1 gene expression measured in PBLs could reflect the ability of blood cells to express this gene in the brain. Thus, our results raise the possibility that HMOX1 variants may control the risk of CM through regulation of HMOX1 expression. Nevertheless, we cannot exclude that relative differences in PBLs lineages among CM patients could influence HMOX1 gene expression measurements. Overall, the HMOX1 association analysis suggests that the involvement of this gene in determining the outcome of malaria infection might be related to disease severity and not exclusively to the CM syndrome.
Our results highlight the value of a study design centered in a particular clinical manifestation of malaria, namely CM, and comparing that clinical manifestation with other frequent outcomes, including non-cerebral forms of malaria (SnC group) or uncomplicated malaria (UM group). This enabled the identification of genetic variants associated with CM in the context of severe malaria syndromes. Thus, the clinical stratification approach we used in this study could prove a useful tool to complement genome-wide association analysis of malaria [48], [49] allowing fine-resolution association analysis of candidate genes and taking into account the haplotypic context and the genetic variability within African populations [50], [51].
This work represents an initial study on genetics of CM in the Angolan population. The statistical evidence supporting the described genetic association findings was limited by the relatively small sample size, by the genetic heterogeneity of the Angolan population, by the presumed small genetic effects and by lack of knowledge of yet to describe causal genetic variants possibly located within the CM associated regions. Although these factors limit the extrapolation of the association results to other populations our findings suggest the novel hypothesis that in African populations the TGFB2 and HMOX1 genetic variation correlate with malaria genetic susceptibility.
This report may provide insights and motivation for future association and functional studies aiming to replicate our findings in larger populations and to confirm a role for TGFB2 and HMOX1 in determining the course of disease in the response to Plasmodium infections.
Materials and Methods
Subjects, Phenotypic and Clinical Criteria
A total of 749 children, living in Luanda and ranging from 6 months to 13 years of age were enrolled in the present study. Ethical permission for this study was granted by the Ethical Committee of the Hospital Pediátrico David Bernardino (HPDB) in Luanda that was appointed by the Angolan Ministry of Health. Written, informed consent was obtained from the parents or guardians of each child. Patients were selected among attendance to the HPDB. Uninfected controls were sampled from children randomly selected in the vaccination ward of the HPDB. The sample collection was carried out from February 2005 to May 2007 and comprised 130 CM children (cases), 158 patients with severe malaria but not CM (SnC), 142 patients with uncomplicated malaria (UM) and 319 uninfected controls (UIF). The mean age in months was 54.2 for CM cases, 45.9 for SnC patients, 50.3 for UM patients and 60.9 for UIF controls. Ethnicity was defined on a parental basis and the major ethnic groups were the Kimbundu and Umbundu. Malaria was diagnosed on the basis of a positive asexual parasitaemia detected on a Giemsa-stained thick smear. For parasitaemia quantification the number of parasites per 100 high-power microscopic fields was estimated and the parasite density was calculated from this value [52]. CM was defined according to the WHO criteria: a coma score <3 in Blantyre Scale for children <60 months or a coma score <7 in Glasgow Scale for children ≥60 months. Meningitis and encephalitis were ruled out by cerebrospinal fluid analysis after lumbar puncture. Exclusion criteria were a different known aetiology of encephalopathy and hypoglycaemia (glycaemia <40 mg/dl). A fraction of CM patients showed as part of the CM syndrome additional clinical complications such as severe malaria anaemia (SMA) and hyperparasitaemia (HPM). The severe non-CM group (SnC) included patients with SMA and/or HPM. Children with HPM had ≥100 red blood cells parasitized by one high-power microscopic field and SMA was defined by haemoglobin <5 g/dl or hematocrit <15%. Patients with consciousness disturbances or with other disease were excluded from this group. The uncomplicated malaria (UM) group represents patients with malaria diagnosis and febrile illness without any clinical finding suggestive of other causes of infection and with no manifestations of severe malaria. All the uncomplicated malaria patients were outpatients. Patient treatments followed the established hospital guidelines. Enrollment of uninfected controls excluded children with any clinical finding suggestive of illness and the uninfected status was confirmed by the absence of Plasmodium DNA in the peripheral blood as detected by PCR [53].
Sample Preparation, Candidate Genes and Genotyped SNPs
Genomic DNA was extracted from whole blood using the Chemagen Magnetic Bead Technology. DNA preparations were quantified using PicoGreen reagents according to the supplier instructions.
All the 749 children and their mothers were genotyped and analyzed for 54 SNPs in the following genes: TGFB2 (Transforming growth factor, beta 2), HBB (β-globin), CD36 (Thrombospondin receptor), HMOX1 (Heme oxygenase 1) and ICAM 1 (Intercellular adhesion molecule 1, CD54). A complete list of genotyped SNPs is available in Table S1.
When selecting SNPs to genotype we prioritized those that were previously associated with CM, preferentially confirmed by sequencing and those belonging to distinct LD blocks as judged from the HapMap database (http://www.hapmap.org). The SNP genotyping method used the Mass Array system to design multiplex reactions for PCR and iPlex primer extension (Sequenom) and the MALDI-TOF based Mass Array genotype platform (Sequenom). Genotyping quality control selected SNPs that yield correct genotypes according HapMap control samples, passed the Hardy-Weinberg equilibrium test (P>0.05) with a calling rate >90%.
HMOX1 gene specific expression by qRT-PCR
For gene-specific expression, total RNA was obtained from 500 µl of whole blood from 42 CM patients collected to RNAlater solution (to stabilize RNA) using the RiboPure-Blood Kit (Ambion). One microgram of total RNA was converted to cDNA (Transcriptor High Fidelity cDNA Synthesis Kit, Roche) and HMOX1 expression was quantified using TaqMan Gene Expression Assay Hs00157965_m1 from ABI with TaqMan Gene Expression master mix. The qRT-PCR was performed according to the manufacturer instructions on an ABI Prism 7900HT system. Relative quantification of HMOX1 specific mRNA was obtained after normalization for HPRT1 mRNA expression measured in the same PCR reaction.
Data analysis
Data analysis followed a step-wise strategy. Firstly, case-control association analysis was performed for individual SNPs. Genotypic association of SNPs with susceptibility to malaria was analyzed using logistic regression models implemented in the SNPassoc package [54] for the R statistical software (version 2.7.0). An uncorrected significance level of the likelihood ratio test (P-value <0.05) was considered as suggestive evidence for association. Secondly, linkage disequilibrium (LD) maps were constructed exclusively for genes that demonstrated multiple SNPs associated to CM only when comparing CM cases to SnC controls. This analysis only pursued for genes that showed multiple polymorphisms associated to CM, exclusively when comparing CM patients to SnC patients. Pairwise LD was estimated using the correlation coefficient (r2) and the LD map was visualized using the snp.plotter R package [55] and was used to identify LD blocks within gene regions. In a third step, we used the gene counting method available in the R package GAP to reconstruct haplotypes composed with multiple associated SNPs within an LD block [56]. We intended to determine whether particular genotypic combinations of the CM-associated SNPs could confer a single joint specific effect on the risk of developing CM. Association between haplotypes and susceptibility to malaria was also evaluated with logistic regression analysis, using the same package that calculates global and haplotype-specific score statistics. Significant global score statistics (GSS) indicate that haplotype-specific scores (hap-score) are valid. A significant positive hap-score means that the haplotype confers susceptibility to the trait and a negative hap-score relates to protection [57]. For each haplotype, empirical P-values were calculated by permutation tests using the Monte-Carlo method to generate 1000 data simulations.
Finally, a family-based test using the TRANSMIT software [58] was performed for HbS and for haplotypes significantly associated to CM. Due to difficulties in collecting both parents, we took the conservative option of restricting the analysis to mother-descendent pairs performing the TDT analysis in children and their mothers.
Haplotypic quantitative trait association analysis for HMOX1 gene expression levels was analyzed by the logistic regression model implemented in the Plink Package (version number 1.06) which calculates P values using empirical significance (http://pngu.mgh.harvard.edu/purcell/plink/) [59].
Supporting Information
Table S1.
Genotyped single nucleotide polymorphisms.
https://doi.org/10.1371/journal.pone.0011141.s001
(0.08 MB DOC)
Table S2.
Allelic distribution of the HMOX1 GT repeat in distinct patient groups and uninfected controls.
https://doi.org/10.1371/journal.pone.0011141.s002
(0.04 MB DOC)
Table S3.
Cerebral malaria association tests for single-nucleotide polymorphisms (SNPs) in the HMOX1, CD36 and ICAM1 genes.
https://doi.org/10.1371/journal.pone.0011141.s003
(0.07 MB DOC)
Figure S1.
Relative frequency distribution of the HMOX1 repeat alleles in distinct malaria phenotypes.
https://doi.org/10.1371/journal.pone.0011141.s004
(0.06 MB DOC)
Acknowledgments
We are grateful to the children that participated in this study. We thank the staff of the Hospital Pediátrico David Bernardino, the medical doctors and nurses of the Unidade de Cuidados II ward and specially Prof. Luís Bernardino, Prof. Victória E. Santo and nurse Domingas Augusto. We also acknowledge Drs. Filomena Silva, Emingarda Castelo Branco and Rosa Silva from the Instituto de Saúde Pública de Angola for the help on sample processing. We thank Drs. Rui Pinto and Luís Caetano from the Clínica Sagrada Esperança de Luanda for sample collection materials. We wish to thank Dr. Miguel Soares for insightful discussions and useful suggestions.
Author Contributions
Conceived and designed the experiments: MRS MJT CB VQ LAG SM OW AC CPG. Performed the experiments: MRS MJT CB VQ LAG MIM SM. Analyzed the data: MRS MJT RV MIM NS TGC CPG. Contributed reagents/materials/analysis tools: SM OW AC CPG. Wrote the paper: MRS RV NS TGC SM AC CPG.
References
- 1. Mackintosh CL, Beeson JG, Marsh K (2004) Clinical features and pathogenesis of severe malaria. Trends Parasitol 20: 597–603.
- 2. Issifou S, Kendjo E, Missinou MA, Matsiegui PB, Dzeing-Ella A, et al. (2007) Differences in presentation of severe malaria in urban and rural Gabon. Am J Trop Med Hyg 77: 1015–1019.
- 3. Newton CR, Krishna S (1998) Severe falciparum malaria in children: current understanding of pathophysiology and supportive treatment. Pharmacol Ther 79: 1–53.
- 4. Idro R, Jenkins NE, Newton CR (2005) Pathogenesis, clinical features, and neurological outcome of cerebral malaria. Lancet Neurol 4: 827–840.
- 5. Mackinnon MJ, Mwangi TW, Snow RW, Marsh K, Williams TN (2005) Heritability of malaria in Africa. PLoS Med 2: e340.
- 6. Rihet P, Traore Y, Abel L, Aucan C, Traore-Leroux T, et al. (1998) Malaria in humans: Plasmodium falciparum blood infection levels are linked to chromosome 5q31–q33. Am J Hum Genet 63: 498–505.
- 7. Williams TN, Mwangi TW, Wambua S, Alexander ND, Kortok M, et al. (2005) Sickle cell trait and the risk of Plasmodium falciparum malaria and other childhood diseases. J Infect Dis 192: 178–186.
- 8. Ruwende C, Hill A (1998) Glucose-6-phosphate dehydrogenase deficiency and malaria. J Mol Med 76: 581–588.
- 9. Guindo A, Fairhurst RM, Doumbo OK, Wellems TE, Diallo DA (2007) X-linked G6PD deficiency protects hemizygous males but not heterozygous females against severe malaria. PLoS Med 4(3): e66.
- 10. Ruwende C, Khoo SC, Snow RW, Yates SN, Kwiatkowski D, et al. (1995) Natural selection of hemi- and heterozygotes for G6PD deficiency in Africa by resistance to severe malaria. Nature 376: 246–249.
- 11. Williams TN, Wambua S, Uyoga S, Macharia A, Mwacharo JK, et al. (2005) Both heterozygous and homozygous alpha+ thalassemias protect against severe and fatal Plasmodium falciparum malaria on the coast of Kenya. Blood 106: 368–371.
- 12. Mockenhaupt FP, Ehrhardt S, Gellert S, Otchwemah RN, Dietz E, et al. (2004) Alpha(+)-thalassemia protects African children from severe malaria. Blood 104: 2003–2006.
- 13. Wambua S, Mwangi TW, Kortok M, Uyoga SM, Macharia AW, et al. (2006) The effect of alpha+-thalassaemia on the incidence of malaria and other diseases in children living on the coast of Kenya. PLoS Med 3: e158.
- 14. Kwiatkowski DP (2005) How malaria has affected the human genome and what human genetics can teach us about malaria. Am J Hum Genet 77: 171–192.
- 15. Aitman TJ, Cooper LD, Norsworthy PJ, Wahid FN, Gray JK, et al. (2000) Malaria susceptibility and CD36 mutation. Nature 405: 1015–1016.
- 16. Koch O, Awomoyi A, Usen S, Jallow M, Richardson A, et al. (2002) IFNGR1 gene promoter polymorphisms and susceptibility to cerebral malaria. J Infect Dis 185: 1684–1687.
- 17. Burgner D, Usen S, Rockett K, Jallow M, Ackerman H, et al. (2003) Nucleotide and haplotypic diversity of the NOS2A promoter region and its relationship to cerebral malaria. Hum Genet 112: 379–386.
- 18. Gyan BA, Goka B, Cvetkovic JT, Kurtzhals JL, Adabayeri V, et al. (2004) Allelic polymorphisms in the repeat and promoter regions of the interleukin-4 gene and malaria severity in Ghanaian children. Clin Exp Immunol 138: 145–150.
- 19. e Pinto EA, Alves JG (2008) The causes of death of hospitalized children in Angola. Trop Doct 38: 66–67.
- 20. Takeda M, Kikuchi M, Ubalee R, Na-Bangchang K, Ruangweerayut R, et al. (2005) Microsatellite polymorphism in the heme oxygenase-1 gene promoter is associated with susceptibility to cerebral malaria in Myanmar. Jpn J Infect Dis 58: 268–271.
- 21. Cunha-Rodrigues M, Portugal S, Febbraio M, Mota MM (2007) Bone marrow chimeric mice reveal a dual role for CD36 in Plasmodium berghei ANKA infection. Malar J 6: 32.
- 22. Bauer PR, Van Der Heyde HC, Sun G, Specian RD, Granger DN (2002) Regulation of endothelial cell adhesion molecule expression in an experimental model of cerebral malaria. Microcirculation 9: 463–470.
- 23. Kun JF, Klabunde J, Lell B, Luckner D, Alpers M, et al. (1999) Association of the ICAM-1Kilifi mutation with protection against severe malaria in Lambarene, Gabon. Am J Trop Med Hyg 61: 776–779.
- 24. Fernandez-Reyes D, Craig AG, Kyes SA, Peshu N, Snow RW, et al. (1997) A high frequency African coding polymorphism in the N-terminal domain of ICAM-1 predisposing to cerebral malaria in Kenya. Hum Mol Genet 6: 1357–1360.
- 25. Bellamy R, Kwiatkowski D, Hill AV (1998) Absence of an association between intercellular adhesion molecule 1, complement receptor 1 and interleukin 1 receptor antagonist gene polymorphisms and severe malaria in a West African population. Trans R Soc Trop Med Hyg 92: 312–316.
- 26. Fry AE, Auburn S, Diakite M, Green A, Richardson A, et al. (2008) Variation in the ICAM1 gene is not associated with severe malaria phenotypes. Genes Immun 9: 462–469.
- 27. Ndiaye R, Sakuntabhai A, Casademont I, Rogier C, Tall A, et al. (2005) Genetic study of ICAM1 in clinical malaria in Senegal. Tissue Antigens 65: 474–480.
- 28. Pain A, Urban BC, Kai O, Casals-Pascual C, Shafi J, et al. (2001) A non-sense mutation in Cd36 gene is associated with protection from severe malaria. Lancet 357: 1502–1503.
- 29. Fry AE, Ghansa A, Small KS, Palma A, Auburn S, et al. (2009) Positive selection of a CD36 nonsense variant in sub-Saharan Africa, but no association with severe malaria phenotypes. Hum Mol Genet 29: 29.
- 30. Williams TN, Mwangi TW, Roberts DJ, Alexander ND, Weatherall DJ, et al. (2005) An immune basis for malaria protection by the sickle cell trait. PLoS Med 2: e128.
- 31. Aitman TJ, Cooper LD, Norsworthy PJ, Wahid FN, Gray JK, et al. (2000) Malaria susceptibility and CD36 mutation. Nature 405: 1015–1016.
- 32. Omi K, Ohashi J, Patarapotikul J, Hananantachai H, Naka I, et al. (2003) CD36 polymorphism is associated with protection from cerebral malaria. Am J Hum Genet 72: 364–374.
- 33. Omi K, Ohashi J, Naka I, Patarapotikul J, Hananantachai H, et al. (2002) Polymorphisms of CD36 in Thai malaria patients. Southeast Asian J Trop Med Public Health 33: 1–4.
- 34. Soares MP, Seldon MP, Gregoire IP, Vassilevskaia T, Berberat PO, et al. (2004) Heme oxygenase-1 modulates the expression of adhesion molecules associated with endothelial cell activation. J Immunol 172: 3553–3563.
- 35. Seldon MP, Silva G, Pejanovic N, Larsen R, Gregoire IP, et al. (2007) Heme oxygenase-1 inhibits the expression of adhesion molecules associated with endothelial cell activation via inhibition of NF-kappaB RelA phosphorylation at serine 276. J Immunol 179: 7840–7851.
- 36. Ohashi J, Naka I, Patarapotikul J, Hananantachai H, Looareesuwan S, et al. (2001) Absence of association between the allele coding methionine at position 29 in the N-terminal domain of ICAM-1 (ICAM-1(Kilifi)) and severe malaria in the northwest of Thailand. Jpn J Infect Dis 54: 114–116.
- 37. Bienzle U, Eggelte TA, Adjei LA, Dietz E, Ehrhardt S, et al. (2005) Limited influence of haptoglobin genotypes on severe malaria in Ghanaian children. Trop Med Int Health 10: 668–671.
- 38. Aucan C, Walley AJ, Hennig BJ, Fitness J, Frodsham A, et al. (2003) Interferon-alpha receptor-1 (IFNAR1) variants are associated with protection against cerebral malaria in the Gambia. Genes Immun 4: 275–282.
- 39. de Souza JB, Riley EM (2002) Cerebral malaria: the contribution of studies in animal models to our understanding of immunopathogenesis. Microbes Infect 4: 291–300.
- 40. Engwerda C, Belnoue E, Gruner AC, Renia L (2005) Experimental models of cerebral malaria. Curr Top Microbiol Immunol 297: 103–143.
- 41. Stevenson MM, Gros P, Olivier M, Fortin A, Serghides LCerebral malaria: human versus mouse studies. Trends Parasitol 2010: 8.
- 42. Bagot S, Campino S, Penha-Goncalves C, Pied S, Cazenave PA, et al. (2002) Identification of two cerebral malaria resistance loci using an inbred wild-derived mouse strain. Proc Natl Acad Sci U S A 99: 9919–9923.
- 43. Omer FM, Kurtzhals JA, Riley EM (2000) Maintaining the immunological balance in parasitic infections: a role for TGF-beta? Parasitol Today 16: 18–23.
- 44. Omer FM, Riley EM (1998) Transforming growth factor beta production is inversely correlated with severity of murine malaria infection. J Exp Med 188: 39–48.
- 45. Walther M, Tongren JE, Andrews L, Korbel D, King E, et al. (2005) Upregulation of TGF-beta, FOXP3, and CD4+CD25+ regulatory T cells correlates with more rapid parasite growth in human malaria infection. Immunity 23: 287–296.
- 46. Deininger MH, Kremsner PG, Meyermann R, Schluesener HJ (2000) Differential cellular accumulation of transforming growth factor-beta1, -beta2, and -beta3 in brains of patients who died with cerebral malaria. J Infect Dis 181: 2111–2115.
- 47. Pamplona A, Ferreira A, Balla J, Jeney V, Balla G, et al. (2007) Heme oxygenase-1 and carbon monoxide suppress the pathogenesis of experimental cerebral malaria. Nat Med 13: 703–710.
- 48. Jallow M, Teo YY, Small KS, Rockett KA, Deloukas P, et al. (2009) Genome-wide and fine-resolution association analysis of malaria in West Africa. Nat Genet 24: 24.
- 49. The Malaria Genomic Epidemiology Network (2008) A global network for investigating the genomic epidemiology of malaria. Nature 456: 732–737.
- 50. Liu YJ, Papasian CJ, Liu JF, Hamilton J, Deng HW (2008) Is replication the gold standard for validating genome-wide association findings? PLoS ONE 3: e4037.
- 51. Chanock SJ, Manolio T, Boehnke M, Boerwinkle E, Hunter DJ, et al. (2007) Replicating genotype-phenotype associations. Nature 447: 655–660.
- 52. Greenwood BM, Armstrong JR (1991) Comparison of two simple methods for determining malaria parasite density. Trans R Soc Trop Med Hyg 85: 186–188.
- 53. Snounou G (1996) Detection and identification of the four malaria parasite species infecting humans by PCR amplification. Methods Mol Biol 50: 263–291.
- 54.
Gonzalez Fau - Armengol L Jr, Armengol L Fau - Sole X, Sole X Fau - Guino E, Guino E Fau - Mercader JM, Mercader Jm Fau - Estivill X, et al. SNPassoc: an R package to perform whole genome association studies.
- 55. Luna A, Nicodemus KK (2007) snp.plotter: an R-based SNP/haplotype association and linkage disequilibrium plotting package. Bioinformatics 23: 774–776.
- 56. Zhao JH, Lissarrague S, Essioux L, Sham PC (2002) GENECOUNTING: haplotype analysis with missing genotypes. Bioinformatics 18: 1694–1695.
- 57. Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA (2002) Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet 70: 425–434.
- 58. Clayton D (1999) A generalization of the transmission/disequilibrium test for uncertain-haplotype transmission. Am J Hum Genet 65: 1170–1177.
- 59. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.