Background. Haptoglobin (Hp) genotype determines the efficiency of hemoglobin clearance after malaria-induced hemolysis and alters antioxidant and immune functions. The Hp2 allele is thought to have spread under strong selection pressure, but it is unclear whether this is due to protection from malaria or other diseases.
Methods. We monitored the incidence of febrile malaria and other childhood illnesses with regard to Hp genotype in a prospective cohort of 312 Kenyan children during 558.3 child-years of follow-up. We also conducted 7 cross-sectional surveys to determine the prevalence of Plasmodium falciparum parasitemia.
Results. The Hp2/2 genotype was associated with a 30% reduction in clinical malarial episodes (adjusted incidence rate ratio, 0.67; P = .008 for Hp2/2 vs. Hp1/1 and Hp2/1 combined). Protection increased with age; there was no protection in the first 2 years of life, 30% protection at ⩾2 years of age, and 50% protection from 4–10 years of age. Children with the Hp1/1 genotype had a significantly lower rate of nonmalarial fever (P = .001).
Conclusions. Balancing selection pressures may have influenced the spread of the Hp gene. Our observations suggest that the Hp2 allele may have spread as a result of protection from malaria, and the Hp1 allele may be sustained by protection from other infections.
Malaria is an important cause of childhood mortality [1] and has exerted selective pressure on the human genome. Haptoglobin (Hp) is an acute-phase plasma protein that binds rapidly and irreversibly with free hemoglobin (Hb) after the occurrence of malaria-induced hemolysis. Hp thus protects against Hb-iron-mediated oxidant damage [2] and modulates immune function and host response to inflammation. Binding of Hp-Hb complexes to CD163 macrophage receptors after hemolysis induces cytokine production [3]. Interestingly, Hp also alters leukocyte trafficking [4] and host response to bacterial lipopolysaccharides [5]. The functions of Hp are genotype dependent [6]. Two codominant alleles (Hp1 and Hp2) on chromosome 16q22 determine 3 phenotypes (the homodimer Hp1-1, the linear polymer Hp2-1, and the large circular polymer Hp2-2). The Hp2 allele is estimated to have arisen in India from an intragenic duplication of the Hp1 allele [7] and has since spread under strong selective pressure. Although the incidence and clinical expression of a number of diseases appear to be influenced by Hp type (reviewed in [6, 8]), it is unclear what selection pressures may have accounted for the spread of the Hp2 allele.
Clear functional differences exist between the Hp genotypes. Homozygosity for the Hp2 allele is associated with increased Hb-iron-induced oxidant damage [9], redox-active iron levels [10], and Hb-derived oxidation of lipids and proteins [11]. These differences are reflected in vivo. Serum vitamin C levels and serum ferroxidase levels in smokers are significantly reduced in individuals with the Hp2-2 phenotype (49.9 µmol/L and 52.9 U/L, respectively), but do not differ in individuals with the Hp1-1 (61.5 µmol/L and 63.4 U/L, respectively) or Hp2-1 (63.7 µmol/L and 65.2 U/L, respectively) phenotype [12, 13]. Oxidative damage to the erythrocyte membrane accelerates phagocytic removal of RBCs via aggregation of band 3 and deposition of immune complexes [14]. Interestingly, evidence suggests that this is a common mechanism for accelerated phagocytic removal of ring-parasitized erythrocytes in G6PD deficiency [15], sickle trait, and β-thalassaemia [16]. Moreover, Hp genotype appears to determine differences in immune function after hemolysis by altering the macrophage cytokine response [3] and modulating Th1 and Th2 balance [17].
Whether the spread of the Hp2 allele is due to protection from malaria or other diseases is controversial. Increased parasite densities were found in Hp knock-out mice [18], and in vitro Hp was toxic to malaria parasites (in order of toxicity, Hp1-1, followed by Hp2-2, followed by Hp2-1) [19]. On the basis of Hp phenotyping, a number of case-control and cross-sectional studies have suggested that Hp2 homozygosity has a protective effect and that there is increased susceptibility among individuals with the Hp1-1 phenotype, especially in relation to severe malaria and placental infection [20,21,22,23–24]. However, 2 additional studies did not find an association with Hp genotype [25, 26]. Previous studies have been cross-sectional or case-control in design. Here, we report on the effects of Hp genotype in a prospective cohort of children among whom we have measured the incidence of clinical malaria and other febrile diseases.
Study design. The study involved a cohort of 312 children (age, 0–10 years) living in the Ngerenya area of Kilifi District on the coast of Kenya and was conducted from August 1998 to August 2001. The study was nested within a larger study described in detail elsewhere [27, 28]. Children were mainly Giriama, a subgroup of the Mijikenda ethnic group. Figure 1 provides an overview of the study design. Children were recruited at the start of the study or at birth, if they were born into study households during the study period. The cohort was monitored by weekly active surveillance for clinical events, with a focus on malaria, and mothers could also bring their children to the dedicated research outpatient clinic at any time if the child was unwell. At each clinic visit, a clinician obtained a detailed clinical history, conducted an examination, and ascribed a primary and secondary diagnosis. During the study period, 7 cross-sectional surveys were conducted to assess the prevalence of Plasmodium falciparum parasitemia (figure 1). Children exited the study if informed consent was withdrawn or if they moved out of the study area for >2 months.
Construction of a study of the incidence of Plasmodium falciparum malaria among children on the coast of Kenya
Clinical definitions. For the purposes of this study, clinical malaria was defined according to 2 alternative definitions. “Malaria definition 1” was indicated either by a measured fever (axillary temperature, >37.5°C) at the time that the slides were prepared or by a clinical history of fever within the preceding 48 h, in conjunction with a slide positive for blood-stage asexual P. falciparum parasites at any density. “Malaria definition 2” was based on a measured fever in conjunction with the presence of P. falciparum parasites at any density (for children aged <1 year) and at a density of >2500 parasites/µL (for children aged ⩾1 year). Malaria definition 2, derived by means of multiple logistic regression methods, is associated with both a sensitivity and specificity of >80% for the diagnosis of clinical malaria in this area [28]. Asymptomatic malaria parasitemia was defined during cross-sectional surveys as a slide positive for P. falciparum in the absence of fever or other symptoms of clinical illness.
Nonmalarial fever was defined as a fever in a child whose blood slide was negative for P. falciparum and who had not received treatment with an antimalarial drug in the preceding 21 days; this definition included a range of nonmalarial febrile illnesses under other diagnostic categories. Upper respiratory tract infection was diagnosed in children whose principal symptoms were coryza or pharyngitis. Lower respiratory tract infection was diagnosed in children with cough and fast breathing and/or chest indrawing, in accordance with World Health Organization clinical criteria [29]. Gastroenteritis was defined as diarrhea (>3 watery stools/day), with or without vomiting. Fever of unknown cause was a diagnosis of exclusion and was only assigned if no etiology could be found for the fever. Helminth infection was diagnosed in children who presented to the clinic with a history of passing worms of any species, and skin infection was diagnosed in children who presented with a variety of dermatological conditions, including scabies, boils, and impetigo.
Laboratory procedures. Blood films were prepared with Giemsa stain and examined for malaria parasites in accordance with standard methods. Parasite densities (i.e., the number of parasites per microliter of whole blood) were recorded as a ratio of parasites to WBCs or, for “heavier” infections, as a ratio of parasites to RBCs. Parasite densities were calculated from complete hematological assessments, if available; if not available, parasite densities were calculated on the assumption of a WBC count of 8 ×; 103 cells/µL or a RBC count of 5 ×; 106 cells/µL. Hp genotype was determined by allele-specific PCR, as described previously [30]. Hb types (HbA and HbS) were characterized by electrophoresis using cellulose acetate gels (Helena), and the common African 3.7-kb α-globin deletion was typed by PCR reaction as described elsewhere [31].
Statistical analysis. Analyses were conducted using Stata software, version 8.0 (Stata Corp.). In all of our regression models, we accounted for the potential clustering of events within individual study children by using the sandwich estimator, which inflates CIs and adjusts P values upwards [32]. Variables that were not normally distributed were log-transformed before analysis. Multivariate regression models included the explanatory variables of age (as a continuous variable), sex, ethnic group, and season (defined in 3-month blocks); in addition, sickle cell trait and α+-thalassaemia genotypes were included in multivariate analyses of malaria or Hb levels.
We derived incidence rate ratios (IRRs) for malaria and other diseases using univariate and multivariate Poisson regression. Prior data suggested that the Hp2-2 phenotype has unique properties, compared with the Hp1-1 and Hp2-1 phenotypes [12, 30]. This study was designed to test the hypothesis that these functional differences would result in a reduced incidence of clinical malarial episodes in individuals with the Hp2/2 genotype. Thus, for malaria, the primary analysis compared individuals who had the Hp2/2 genotype with those who had the Hp2/1 or Hp1/1 genotype, and a secondary analysis treated each genotype separately. Because we did not have a prior hypothesis regarding the effect of Hp genotype on nonmalarial disease, each genotype was considered separately. Children were considered not to be at risk for malaria (and were omitted from both the numerator and denominator populations) for 21 days after receipt of treatment with an antimalarial drug. ORs for the prevalence of asymptomatic parasitemia over 7 cross-sectional surveys were derived by logistic regression analysis. Log-transformed parasite densities were analyzed using linear regression models. The likelihood ratio test was used to assess for interactions between explanatory variables, as appropriate; no significant interactions were found. We investigated the possibility that age may have acted as an effect modifier in the association between malaria and Hp genotype by comparing models that included or excluded interaction terms between Hp genotype and age, using the Wald test. The duration of follow-up for study children was accounted for in all analyses.
Ethics approval for the study was granted by the Kenya Medical Research Institute National Ethical Review Committee, and procedures complied with the Helsinki Declaration of 1975 (revised in 1983). Written informed consent was provided by children's parents or guardians before recruitment.
Characteristics of study population. A total of 4437 clinic visits were made by 312 children during 558.3 child-years of follow-up, and 7 cross-sectional surveys of parasite prevalence were conducted during this time (figure 1). Of these children, 87 (27.9%) carried the Hp1/1 genotype, 146 (46.8%) carried the Hp2/1 genotype, and 79 (25.3%) carried the Hp2/2 genotype. Hp genotypes were in Hardy-Weinberg equilibrium, which predicts allele and genotype frequencies within a nonevolving population. No association was found between Hp genotype and either sickle cell trait or α+-thalassaemia genotype. Children with the Hp1/1 genotype were slightly younger than children with the Hp2/1 or Hp2/2 genotype (P = .08); the effect of age was accounted for in all multivariate analyses. Hb levels were measured in 158 children at the start of the study (table 1). Hb levels were lower in children with the Hp1/1 genotype than in others (P = .08, by log likelihood ratio testing), probably because of their younger age. Table 1 summarizes the characteristics of the study population according to Hp genotype.
Hp genotype and malaria. In primary analysis, the incidence of clinical malaria was significantly lower in children with the Hp2/2 genotype, compared with children with the other genotypes; the degree of protection was ∼30% for each of the 2 definitions of clinical malaria considered (for definition 1 [i.e., any density of parasitemia], the adjusted IRR was 0.72 [95% CI, 0.54–0.95; P = .02]; for definition 2 [parasite density of >2500 parasites/µL for children aged ⩾1 year], the adjusted IRR was 0.67 [95% CI, 0.49–0.90; P = .008]) (table 2). A secondary analysis comparing the 3 genotypes also indicated a reduced incidence of clinical malaria for children with the Hp2/2 genotype (for definition 1, the adjusted IRR was 0.70 [95% CI, 0.50–0.96; P = .029]; for definition 2, the adjusted IRR was 0.70 [95% CI, 0.49–0.99; P = .049]). The incidence of malaria did not differ among children with the Hp1/1 genotype versus those with the Hp2/1 genotype (for definition 1, the adjusted IRR was 0.96 [95% CI, 0.75–1.22; P = .72]; for definition 2, the adjusted IRR was 1.08 [95% CI, 0.82–1.42; P = .58]).
For all age groups combined, the Hp2/2 genotype was associated with ∼30% protection against clinical malarial episodes. However, this rate of protection varied with age: there was no protection in the first 2 years of life, a 30% rate of protection at ⩾2 years of age, and a ∼50% rate of protection from 4–10 years (figure 2). Although a significant interaction was not found between age and protection over the full range of ages (χ24, 4.26; P = .37), the data support the impression of acquired protection with age.
Incidence rate ratios (IRRs) for clinical malaria in children with haptoglobin (Hp)2/2 genotype versus the Hp1/1 and Hp2/1 genotypes combined, by age. Because this was a longitudinal cohort study, some children contributed data to >1 age stratum. Adjusted IRRs for malaria in children with Hp2/2, compared with baseline (individuals with Hp1/1 and those with Hp2/1) were calculated individually for each age stratum using a Poisson regression model, which included each definition of malaria as a dependent variable and the explanatory variables hemoglobin type (HbAA or HbAS), α+-thalassaemia genotype, season (defined as 3-month blocks), sex, ethnic group, and age (as a continuous variable). Ninety-five percent CIs and significance values were adjusted to take account of potential within-subject clustering of events using the sandwich estimator [32].
The prevalence of asymptomatic parasitemia during the 7 cross-sectional surveys was marginally lower in the Hp2/2 genotype group than in the other groups (adjusted OR for Hp2/2 vs. Hp1/1 and Hp2/1 combined, 0.64 [95% CI, 0.42–0.99; P = .043]; adjusted OR for heterogeneity between the 3 genotypes, 0.62 [95% CI, 0.37–1.0; P = .065) (table 1). Parasite prevalence did not differ between the Hp1/1 genotype group and the Hp2/1 genotype group (adjusted OR, 0.95; 95% CI, 0.64–1.4; P = .79). The overall mean parasite density did not differ by Hp genotype, either during febrile malaria episodes (P = .24) or during asymptomatic infection (as detected by 7 cross-sectional surveys; P = .82) (figure 3). However, in subgroup analysis, individuals with the Hp2/2 genotype aged ⩾4 years had lower parasite densities during febrile malaria episodes, compared with children who carried the other genotypes (coefficient for Hp2/2 vs. Hp1/1 and Hp2/1 combined, -0.22; 95% CI -0.43 to -0.13; P = .038).
Geometric mean parasite densities, by haptoglobin (Hp) genotype. The data on asymptomatic parasitemia derive from 7 cross-sectional surveys of malaria parasite prevalence in the study cohort. The data on Hp1/1 and Hp2/1 reflect 291 slides that tested positive for malaria parasites (of 1197 blood slides) among 233 children, and the data for Hp2/2 reflect 85 slides that tested positive for malaria parasites (of 428 blood slides) among 79 children. The data on clinical malaria derive from episodes of malaria (defined as an axillary temperature of >37.5°C in association with a slide positive for blood-stage asexual Plasmodium falciparum at any parasite density) detected during active and passive surveillance over 558.3 child-years of follow-up. The data on Hp1/1 and Hp2/1 reflect 964 episodes of malaria among 233 children, and the data on Hp2/2 reflect 280 episodes of malaria among 79 children. Comparisons were made by linear regression with adjustment for the effects of age, season, sex, ethnic group, haemoglobin type (HbAA or HbAS), and α+-thalassaemia genotype. For asymptomatic parasitemia, the coefficient was 0.026 (95% CI, -0.20 to 0.25; P = .82), and for clinical malaria, the coefficient was -0.08 (95% CI, -0.22 to 0.06; P = .24). Bars, 95% CIs.
Hp genotype and nonmalarial diseases. We did not have a prior hypothesis regarding the effect of Hp genotype on nonmalarial childhood diseases. However, we found that the Hp1/1 genotype was associated with a 30% reduction in the incidence of nonmalarial febrile episodes (adjusted IRR, 0.70; 95% CI, 0.56–0.86; P = .001) (table 3). The incidence of nonmalarial febrile episodes did not differ in children with the Hp2/1 genotype versus those with the Hp2/2 genotype (adjusted IRR, 0.98; 95% CI, 0.79–1.21; P = .83). We did not find a significant association between Hp1/1 genotype and any particular primary nonmalarial clinical diagnosis. It is possible that the Hp1/1 genotype protects against a range of febrile illnesses, and our study may have been insufficiently powered to detect differences in individual diseases.
The association between Hp polymorphisms and malaria infection has been the subject of controversy and debate for >30 years. Some studies have suggested that the Hp2-2 phenotype is protective against malaria [20, 21, 23], and other studies have not found a clear association with Hp genotype [25, 26]. However, many of these studies were small, and all were of case-control or cross-sectional design. Here, we report a prospective cohort study that examined the association between Hp genotype and the incidence of clinical malaria and other childhood diseases. We found that the Hp2/2 genotype was associated with a 30% reduction in the incidence of clinical malaria (adjusted IRR for Hp2/2 vs. Hp1/1 and Hp2/1 combined, 0.67; 95% CI, 0.5–0.9; P = .008). Moreover, the degree of protection from malaria was not constant but increased with age, reaching a 50% rate of protection among children from 4–10 years of age. This suggests that, like the sickle cell trait [33], the Hp2/2 genotype may result in accelerated acquisition of immunity to malaria. Although our observations are in agreement with previous studies that have indicated a protective effect of Hp2-2 phenotype [20, 21, 23], we did not find a significant effect of the Hp2/1 genotype.
So how might the Hp2/2 genotype protect against malaria infection? We recently described a greater decrease in Hb level among children with the Hp2/2 genotype and speculated that this may have resulted from a failure to quench Hb-iron-mediated oxidant stress during malaria-induced hemolysis [30]. Studies suggest that reduced antioxidant protection may be specific to Hp2-2 and that there is little difference between the Hp2-1 and Hp1-1 phenotypes in this regard [12, 13]. Experimental evidence suggests that enhanced phagocytosis of ring-parasitized erythrocytes in sickle cell, G6PD deficiency, and β-thalassaemia is mediated by a mechanism essentially similar to the phagocytosis of oxidatively damaged RBCs [15, 16]. Increased reactive oxygen species (from the developing malaria parasite and from increased free Hb and iron levels in the Hp2/2 genotype) may similarly lead to membrane-bound hemichromes and iron, aggregation of band 3 proteins, and binding of autologous IgG and complement. It is possible, therefore, that enhanced immunity in the Hp2/2 genotype may be mediated by accelerated acquisition of antibodies to altered host antigens, such as the band 3 protein expressed on the ring-parasitized erythrocyte. Similar to other genetic polymorphisms associated with oxidant stress, such as G6PD deficiency and β-thalassaemia, Hp2/2 genotype (under conditions of oxidant stress) may protect against malaria at the expense of causing anemia [30].
Differences in host response to inflammation may represent an alternative explanation for the protective effect of the Hp2/2 genotype. The cytokine response produced by macrophages in response to hemolysis differs by Hp genotype. By binding to macrophage CD163 receptors, Hp and Hb complexes mimic antibody binding and trigger cytokine secretion [3]. This scavenging of Hp/Hb complexes induces an anti-inflammatory effect by strongly stimulating IL-10 secretion and heme oxygenase-1 synthesis [34]. Evidence suggests that the Hp2-2/Hb complex may promote a Th1 cytokine response by stimulating the release of less IL-10 and IL-6, compared with the Hp1-1/Hb complex [17]. Enhanced macrophage iron loading in individuals with the Hp2-2 phenotype after hemolysis [35] may also inhibit TNF-α and IFN-γ production [36]. The significance of these possible effects remains to be elucidated. However, it is possible that, by modulating the inflammatory response during hemolysis, individuals who carry the Hp2/2 genotype may have an altered clinical course of malaria infection.
Interestingly, we found that the Hp1/1 genotype was associated with significant protection against nonmalarial febrile illnesses. Mounting evidence suggests that Hp plays an important role in host defence against infection. Hp modulates leukocyte trafficking [4] and host response to bacterial lipopolysaccharides [5]. Hp (of the Hp1-1 type) is a ligand for the Mac-1 leukocyte integrin β2 (CD11b/CD18) receptor on monocytes, macrophages, and natural killer cells [4]. Moreover, Hp suppresses lipopolysaccharide-induced release of TNF-α, IL-10, and IL-12 in vitro and protects the host against the deleterious effect of lipopolysaccharides in vivo [5]. Furthermore, by binding to free Hb and reducing iron availability, Hp exerts a bacteriostatic effect on many organisms such as Escherichia coli [37]. Lower levels of available iron in the Hp1-1 phenotype compared to the Hp2-2 type [10] may thus limit bacterial growth. In clinical studies, Hp2-2 is associated with increased mortality among patients with tuberculosis [38] and HIV infection [39]. It is likely that many of the nonmalarial febrile illnesses suffered would be due to viral infections, and a limitation of the study is lack of microbiology data. However, it seems unlikely that this effect can be explained by protection from any organism individually, but rather that Hp1-1 may protect from a range of different febrile processes.
Our results suggest that homozygosity for the Hp2 allele protects against clinical malaria and that homozygosity for Hp1 protects against nonmalarial febrile illness. Both P. falciparum malaria and bacterial infection are major causes of childhood death in sub-Saharan Africa [1, 40] and have exerted selective pressure on the human genome. Although most febrile illnesses are likely to have been due to viral infections, it is intriguing to speculate that balancing selection pressures on the Hp alleles from malaria and bacteremia may have determined the spread and geographical distribution of the Hp gene. Establishment of the role of the Hp gene and mechanisms of protection may further our understanding of the pathophysiology of malaria and other childhood infections.
We thank the study participants and their families for taking part in the study. This paper is published with the permission of the Director of the Kenya Medical Research Institute.
Financial support. Wellcome Trust and the European Union via the BioMalpar Network 6 Initiative.
Potential conflicts of interest. All authors: no conflicts.
IDSA Members: For your free access to this journal, log in via the IDSA members area.
Open access options for authors visit Oxford Open
This journal enables compliance with the NIH Public Access Policy