Skip to main content

Risk of pneumococcal bacteremia in Kenyan children with glucose-6-phosphate dehydrogenase deficiency



Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common enzyme deficiency state in humans. The clinical phenotype is variable and includes asymptomatic individuals, episodic hemolysis induced by oxidative stress, and chronic hemolysis. G6PD deficiency is common in malaria-endemic regions, an observation hypothesized to be due to balancing selection at the G6PD locus driven by malaria. G6PD deficiency increases risk of severe malarial anemia, a key determinant of invasive bacterial disease in malaria-endemic settings. The pneumococcus is a leading cause of invasive bacterial infection and death in African children. The effect of G6PD deficiency on risk of pneumococcal disease is undefined. We hypothesized that G6PD deficiency increases pneumococcal disease risk and that this effect is dependent upon malaria.


We performed a genetic case-control study of pneumococcal bacteremia in Kenyan children stratified across a period of falling malaria transmission between 1998 and 2010.


Four hundred twenty-nine Kenyan children with pneumococcal bacteremia and 2677 control children were included in the study. Among control children, G6PD deficiency, secondary to the rs1050828 G>A mutation, was common, with 11.2% (n = 301 of 2677) being hemi- or homozygotes and 33.3% (n = 442 of 1329) of girls being heterozygotes. We found that G6PD deficiency increased the risk of pneumococcal bacteremia, but only during a period of high malaria transmission (P = 0.014; OR 2.33, 95% CI 1.19–4.57). We estimate that the population attributable fraction of G6PD deficiency on risk of pneumococcal bacteremia in areas under high malaria transmission is 0.129.


Our data demonstrate that G6PD deficiency increases risk of pneumococcal bacteremia in a manner dependent on malaria. At the population level, the impact of G6PD deficiency on invasive pneumococcal disease risk in malaria-endemic regions is substantial. Our study highlights the infection-associated morbidity and mortality conferred by G6PD deficiency in malaria-endemic settings and adds to our understanding of the potential indirect health benefits of improved malaria control.

Peer Review reports


Glucose-6-phosphate dehydrogenase (G6PD) catalyzes the oxidation of glucose-6-phosphate, producing 6-phosphogluconolactone and dihydronicotinamide adenine dinucleotide phosphate (NADPH) [1]. NADPH is a reducing agent that is central to the mechanism by which cells mitigate oxidative stress. This is especially important for erythrocytes, which, not possessing mitochondria, lack alternative metabolic pathways for NADPH synthesis. In humans, G6PD is encoded by an X-linked locus: G6PD.

G6PD is highly polymorphic. Deleterious G6PD polymorphisms result in the most common human enzyme deficiency: G6PD deficiency [1]. The X-linked inheritance of such polymorphisms means that G6PD deficiency is limited to hemizygous males and homozygous females, while heterozygous females display an intermediate phenotype. G6PD deficiency can cause a wide spectrum of clinical disease, the severity of which is dependent on the level of residual enzyme activity. In individuals with symptomatic G6PD deficiency, disease commonly manifests as episodic hemolysis induced by oxidative stress (e.g., drugs, infection). More severe defects cause chronic hemolysis, characterized by chronic transfusion-dependent anemia, splenomegaly, and the formation of gallstones [2]. An immunodeficiency characterized by granulocyte dysfunction and a clinical phenotype similar to chronic granulomatous disease has also been described in severe enzymatic deficiency [3]. The major genetic determinant of G6PD deficiency in coastal Kenyan populations is the rs1050828 G>A mutation [4], which gives rise to a form of G6PD deficiency that is commonly referred to as the G6PD A− variant. G6PD A− is a WHO class 3 G6PD deficiency variant [5], predicted to result in only mild-to-moderate enzymatic deficiency (10–60% of normal). This genetic variant of G6PD deficiency is estimated to account for 85% of phenotypic variation in G6PD enzymatic activity in coastal Kenyan populations [4]. Among female heterozygous for the G6PD A− variant in coastal Kenya, G6PD enzymatic activity is determined by variation at a second SNP at the G6PD locus, rs1050829 [4]. rs1050829 is a WHO class 4 G6PD deficiency variant, predicted to result in a clinically asymptomatic decrease in G6PD enzymatic activity (< 40% reduction from wild type). Variation at rs1050829 does not affect residual enzymatic activity in individuals with G6PD deficiency secondary to the G6PD A− variant in African populations, as the G6PD A− variant is always inherited on a rs1050829:C background [4].

The geographic distribution of G6PD deficiency is strongly correlated with malaria transmission. This has led to the hypothesis that malaria has driven selection of G6PD deficiency variants in human populations. That hypothesis is supported by evidence of recent selection pressure at G6PD [6] and by studies defining the effect of G6PD variation on malaria risk [7, 8]. The exact nature of the selection pressure imposed by malaria at G6PD remains controversial. We have previously demonstrated that, in Kenyan children, G6PD status is not associated with uncomplicated malaria, but that G6PD deficiency secondary to the G6PD A− variant increases the risk of severe malarial anemia (SMA), whereas heterozygous females are protected against severe malaria [7]. In the same population, we found no evidence for an independent effect of C allele carriage at rs1050829 on risk of uncomplicated or severe malaria [7]. Furthermore, through a multi-population analysis of the effect of variation at G6PD on malaria risk, we have recently demonstrated that G6PD deficiency is associated with opposing additive effects on cerebral malaria and SMA [8]. In that analysis, G6PD deficiency was associated with increased susceptibility to SMA but a lower risk of cerebral malaria [8]. More recently, however, it has been suggested that this observation may represent an artifact of collider bias [9]: by excluding children with severe anemia from the case definition of cerebral malaria, the apparent protective association may be driven by an absence of severe anemia and not cerebral malaria per se. Notwithstanding on-going uncertainty regarding the exact nature of the balancing selection at the G6PD locus driven by malaria, there is a clear and reproducible association between G6PD deficiency and risk of SMA.

The effect of genetic variation at G6PD on infectious diseases other than malaria has not been widely studied. Specifically, the role for genetic variation at the G6PD locus as a determinant of invasive bacterial disease is undefined. Streptococcus pneumoniae is consistently among the most frequently isolated pathogens causing bacteremia in African children [10] and is a leading cause of mortality in children < 5 years [11]. Genetic variation at G6PD modifies risk of SMA. SMA is a major risk factor for community-acquired bacteremia in African children [12, 13], and the malaria-protective effects of HbAS carriage reduce the risk of some of the clinical syndromes that are associated with invasive pneumococcal diseases such as severe pneumonia and meningitis [14]. Moreover, severe anemia in African children, irrespective of its etiology, is associated with increased risk of community-acquired bacteremia [15]. We thus hypothesized that G6PD deficiency increases the risk of pneumococcal bacteremia in African children and that this effect is dependent on intense malaria transmission resulting in severe anemia among G6PD-deficient children. To test this hypothesis, we performed a genetic case-control analysis of the effect of G6PD deficiency on the risk of pneumococcal bacteremia in Kenyan children, across a period of falling malaria transmission.


Study population, genotyping, and quality control

Kenyan children under the age of 13 years, presenting to Kilifi County Hospital, Kenya, between 1 August 1998 and 30 October 2010, with community-acquired bacteremia were recruited to the study. Routine vaccination against S. pneumoniae was not available during the study period and was not introduced until January 2011. During the study period, a blood sample for bacterial culture was obtained from every child admitted to hospital, with the exception of elective surgical admissions and children admitted with trauma. Bacterial culture of blood was performed with a BACTEC 9050 instrument (Becton Dickinson, USA), and pneumococcal isolates identified by optochin susceptibility. Children with a mid-upper arm circumference < 11.5 cm (< 11 cm in children < 6 months of age) were considered to have severe malnutrition. HIV infection status was determined with two rapid antibody tests and with PCR for proviral DNA for children < 18 months. Malaria parasitemia was assessed using thick and thin Giemsa-stained blood films. Among children with all-cause bacteremia (n = 1816), 506 cases of pneumococcal bacteremia were recruited to the study. Healthy community controls (n = 3091) were recruited as part of a birth cohort study at 3–12 months of age from within the same population as cases between 1 May 2006 and 30 April 2008. All children have been subject to longitudinal follow-up.

Genome-wide genotyping was performed in all study samples as part of the Wellcome Trust Case Control Consortium 2 genome-wide association study (GWAS) of all-cause bacteremia [16]. Following genomic DNA extraction and quality control, DNA was whole-genome genotyped on Affymetrix SNP Chip 6.0 arrays. SNPs passing the following quality control (QC) metrics (MAF > 1%, genotype probability (info) > 0.975, plate effect P > 1 × 10−6, Hardy-Weinberg equilibrium P < 1 × 10−10, and SNP missingness > 2%) were included in the analysis. Following SNP QC, genotypes at 787,861 autosomal SNPs were taken forward for computation of principal components (EIGENSTRAT [17]) following linkage disequilibrium pruning, relatedness estimates (identity by descent, PLINK [18]) and heterozygosity (PLINK). Samples were then excluded from the association analysis as follows: relatedness (identity by descent > 0.4), population outliers, extreme heterozygosity, or low genotyping call rate [16].

Genotypes at rs1050828 were directly determined on the genotyping array. Genotyping quality at rs1050828 was good, passing SNP QC thresholds (as above) and with good cluster separation (visualized in Evoker). Genotypes at rs1050829 and the 4 malaria risk loci (ABO, rs8176719; Dantu blood group, rs192804806; sickle hemoglobin [HbS], rs334; ATP2B4, rs4951377) were not directly determined and were imputed with SHAPEIT [19] and IMPUTE2 [20], using 1000G phase 3 as a reference panel. Following quality control, 429 cases of S. pneumoniae bacteremia and 2677 healthy control samples were included in the association analysis.

Statistical analysis

Association analysis at G6PD

We used logistic regression to test for association between pneumococcal bacteremia and G6PD A− variant (rs1050828) genotype under additive (rs1050828 genotypes coded to reflect a monotonic change in G6PD biochemical activity), G6PD deficiency (G6PD deficiency c.f. normal G6PD activity), and G6PD heterozygous (rs1050828:T heterozygous girls c.f. normal G6PD activity) models. In a separate model, we used logistic regression to test for association between pneumococcal bacteremia and rs1050829 genotype among G6PD A− heterozygous females under an additive model. The first four principal components of genome-wide genotyping data were included in each model to account for population sub-structure. In addition, sex was included as a covariate in each model. Statistical analysis was performed in R.

Multinomial logistic regression

As control samples were recruited as a birth cohort towards the end of the study, we used multinomial logistic regression to estimate the effect of the G6PD variation on pneumococcal bacteremia risk in each of the three phases of malaria transmission (pre-decline, pre-2000; decline, 2000–2006; post-decline, post-2006). We used control status as the baseline stratum and cases of pneumococcal bacteremia during each of the specified time periods as strata. We tested for association under additive, G6PD deficiency, and G6PD heterozygous models (as defined above) at rs1050828, and the effect of rs1050829 genotype in rs1050828 heterozygous girls under an additive model. In each case, the first four principal components and sex were included as covariates in the model. To address the possibility of confounding secondary to other known malaria risk loci, we also included genotypes at each locus in the model, with genotypes coded to reflect the model of association observed in malaria (ABO, recessive; Dantu blood group, additive; HbS, heterozygote and recessive; ATP2B4, recessive).

To identify a subset of control samples well-matched to cases with respect to sex and ethnicity (as modeled by the first 4 principal components of genome-wide genotyping data), we used nearest neighbor matching as implemented in MatchIt [21]. The proportions of case children with HIV infection and malnutrition were compared between time periods using χ2 tests. Under the assumption that G6PD deficiency and HIV infection or malnutrition are independent, we performed case-only interaction analysis with logistic regression adjusted for sex and population structure. Statistical analysis was performed in R.

Bayesian comparison of models of association

We compared models of association at rs1050828 with pneumococcal bacteremia across three time periods, as estimated by multinomial logistic regression, using a Bayesian approach. We considered three models of effect across the specified time periods defined by the prior distributions on the effect size:

  • “Null”: effect size = 0, i.e., no association in any time period.

  • “Same”: effect size ~ N (0,a2) and fixed between time periods (ρ = 1).

  • “Pre-decline alone”: effect size ~ N (0,a2) in the pre-decline period only (with no effect in the other time periods).

For heterozygous and G6PD deficiency models, a = 0.5, and for additive models, a = 0.2. For each model, we calculated approximate Bayes factors [22] and posterior probabilities, assuming each model to be equally likely a priori. Statistical analysis was performed in R.

Population attributable fraction

Population attributable fractions were calculated as follows:

$$ \mathrm{PAF}=\frac{P\left(\mathrm{OR}-1\right)}{P\left(\mathrm{OR}-1\right)+1} $$

where PAF = population attributable fraction, OR = odds ratio, and P = population frequency of risk genotype(s).


Study population

Following quality control, 429 Kenyan children with pneumococcal bacteremia and 2677 control participants were included in our genetic association analysis. The first four principal components of genome-wide genotyping data capture self-reported ethnicity and confirm that the population structure of the control samples is representative of the cases (Fig. 1). The mean age among case samples was 1.9 years (range 0–13 years), and 38.7% were female. Among cases, at admission, 25.9% had malnutrition, 12.2% had malaria parasitemia, and 20.8% were HIV-infected (Table 1). Inpatient mortality was 24%. The use of a birth cohort as a healthy control population risks loss of study power through misclassification bias: the inclusion of children as healthy controls who have subsequently experienced an episode of pneumococcal bacteremia. Follow-up of the control children demonstrates that the effect of misclassification bias in our study is negligible. Following longitudinal follow-up to an average age of 5.1 years, there were 8 cases of all-cause bacteremia and 24 deaths among the control children included in the study.

Fig. 1
figure 1

Population structure of study samples. Plots of the first four principal components of genome-wide genotyping data. Individuals are color-coded according to self-reported ethnicity (a) and case-control status (b)

Table 1 Study sample demographics, case comorbidities, and inpatient mortality

Genetic variation causing G6PD deficiency

Among control samples, rs1050828 had a minor allele frequency (MAF) of 0.196. G6PD deficiency, defined by homozygous rs1050828:AA or hemizygous rs1050828:A genotypes, was common (n = 301 of 2677; 11.2%). 3.2% of girls (n = 42 of 1329) and 19.2% of boys (n = 259 of 1348) had homozygous and hemizygous A allele carriage at rs1050828, respectively, resulting in G6PD deficiency. One third (33.3%, n = 442 of 1329) of girls were heterozygous for rs1050828:A carriage. rs1050829 was well-imputed and common in this dataset (imputation info metric = 0.981; MAF = 0.396). In the study population, 27.9% (n = 1175 of 4213) were homozygous or hemizygous for the minor rs1050829:C allele. Among female rs1050828:T heterozygotes, 74.9% (n = 490 of 654) were rs1050829:TC heterozygotes and 25.1% (n = 164 of 654) were rs1050829:CC homozygotes. The 4 malaria risk loci previously described in this population [23] were common and well-imputed among the study samples (rs8176719, imputation info metric = 0.994, MAF = 0.259; rs334, imputation info metric = 0.932, MAF = 0.097; rs192804806, imputation info metric = 0.987, MAF = 0.092; rs4951377, imputation info metric = 0.995, MAF = 0.658).

G6PD status and risk of pneumococcal bacteremia

To investigate whether G6PD status is a determinant of pneumococcal bacteremia risk, we tested for association between the rs1050828 locus and pneumococcal bacteremia. Genotype at rs1050828 (Padditive = 0.261; OR = 1.08, 95% CI 0.94–1.25), G6PD deficiency (P = 0.080; OR = 1.30, 95% CI 0.96–1.74), and heterozygous carriage of the rs1050828:A allele (P = 0.326; OR = 0.83, 95% CI 0.56–1.19) were not associated with risk of pneumococcal bacteremia. In addition, among female rs1050829 G>A heterozygotes, rs1050829 carriage did not modify risk of pneumococcal bacteremia (Padditive = 0.452; OR = 0.75, 95% CI 0.36–1.58).

G6PD deficiency increases risk of pneumococcal bacteremia during an era of high malaria transmission

Over the course of case recruitment, between 1998 and 2008, there was a marked change in malaria transmission intensity in Kilifi [24]. Historically, malaria has been a major risk factor for community-acquired bacteremia in this population [13]. Declining malaria transmission may therefore confound any association between bacteremia and G6PD deficiency. To investigate this, we estimated the effect of G6PD status on pneumococcal bacteremia risk as malaria transmission declined over the period of the study. Dividing the study period into pre-decline (pre-2000), decline (2000–2006), and post-decline (post-2006) periods, we fitted multinomial regression models of pneumococcal bacteremia risk secondary to G6PD status, considering each time period as a stratum (Fig. 2). Under an additive model, the data provide support for a model in which decreasing G6PD activity (estimated by rs1050828 genotype) only increased the risk of pneumococcal bacteremia during the pre-decline period (ORPre-decline = 1.41, 95% CI 1.01–1.97, PPre-decline = 0.044). However, the data more strongly support a model in which G6PD deficiency, defined by homozygous rs1050828:AA or hemizygous rs1050828:A genotypes, increased risk of pneumococcal bacteremia, again in the pre-decline period alone (ORPre-decline = 2.33, 95% CI 1.19–4.57, PPre-decline = 0.014). We found no evidence for a significant effect of heterozygous carriage of rs1050828:A on pneumococcal bacteremia risk in any time period. We saw no evidence for a significant effect of rs1050829:C carriage on pneumococcal bacteremia risk in any time period.

Fig. 2
figure 2

Malaria transmission, G6PD deficiency, and risk of pneumococcal bacteremia. a Age-standardized, annual malaria parasite prevalence in Kilifi, Kenya, as estimated from parasite prevalence among trauma admissions [24, 25]. Ninety-five percent confidence intervals illustrated with red, dashed line. Pre-decline (pre-2000), decline (2000–2006), and post-decline (post-2006) periods used in the analysis are depicted. b Left panels: Log-transformed odds ratios and 95% confidence intervals of G6PD deficiency association with pneumococcal bacteremia risk in pre-decline, decline, and post-decline study periods. Right panels: Posterior probabilities of models of association with G6PD deficiency: “null,” no association with pneumococcal bacteremia in any time period; “same,” the same effect on bacteremia in all three time periods; and “pre-decline alone,” a non-zero effect on pneumococcal bacteremia in the pre-decline time period alone. Association statistics and model posterior probabilities are presented under additive (top), G6PD deficiency risk (middle), and G6PD heterozygous risk (bottom) models

To ensure that our model adequately controls for differences in sex and population structure, we used propensity score matching to identify a subset of control samples (n = 429) matched to cases for sex and population structure. In keeping with our primary analysis, there is again evidence for increased risk of pneumococcal disease among children with G6PD deficiency (ORPre-decline = 2.54, 95% CI 1.23–5.25, PPre-decline = 0.012). There is no evidence for effect of G6PD deficiency on pneumococcal disease risk in the decline and post-decline study periods (PDecline = 0.271, PPost-decline = 0.442).

G6PD deficiency increases risk of pneumococcal bacteremia independent of other risk factors for invasive infection

Both HIV infection and malnutrition are established risk factors for pneumococcal bacteremia in this population [26]. We therefore sought to understand whether changes in the epidemiology of malnutrition or HIV infection during the study could confound our findings. We compared the proportions of cases with HIV co-infection and malnutrition presenting during the pre-decline period and during later years. There was evidence for a temporal change in the proportion of cases with malnutrition (pre-2000, 12 of 47 children, 25.5%; post-2000, 86 of 370 children, 23.2%; P = 5.13 × 10−4), but not with HIV infection (pre-2000, 11 of 44 children, 25%; post-2000, 14 of 86 children, 16.3%; P = 0.233). There was no evidence for interaction between G6PD deficiency and malnutrition or HIV among cases of pneumococcal bacteremia in the study considered as a whole (PMalnutrition = 0.92; PHIV = 0.77) or during the pre-decline period alone (PMalnutrition = 0.85; PHIV = 0.48).

Genome-wide association studies have identified 4 genetic determinants of severe malaria risk in this population (ABO blood group, HbS, Dantu blood group, and ATP2B4) [23]. We investigated whether the observed association between G6PD deficiency and the risk of pneumococcal bacteremia was independent of genetic variation at these loci. Including genetic variation at these loci did not significantly alter the effect estimate of G6PD deficiency on pneumococcal bacteremia risk in the study considered as a whole (OR = 1.28, 95% CI 0.95–1.73, P = 0.101) or during the pre-decline period alone (OR = 2.32, 95% CI 1.18–4.56, P = 0.015).

The population attributable fraction of pneumococcal bacteremia

To better understand the impact of G6PD deficiency on pneumococcal disease risk, we calculated the population attributable fraction prior to the observed decline in malaria transmission at the study site. In the pre-decline period, the pneumococcal bacteremia population attributable fraction for G6PD deficiency was 0.129. For comparison, assuming an OR for pneumococcal bacteremia among children with sickle cell disease (HbSS) of 33.0 [27] and a mean HbSS frequency among children < 14 years in this population of 0.003 [27], the population attributable fraction of pneumococcal bacteremia for HbSS would be 0.088.


In Kenyan children, G6PD deficiency increases risk of pneumococcal bacteremia. This effect is only observed prior to a decline in malaria transmission at the study site. Heterozygous carriage of the rs1050828:A allele does not significantly affect risk of pneumococcal bacteremia in this population. These data demonstrate that G6PD deficiency increases risk of pneumococcal bacteremia in Kenyan children and that this effect is dependent on malaria.

Malaria is a well-established risk factor for community-acquired bacteremia [13]. Children with G6PD deficiency secondary to the G6PD A− variant have increased risk of SMA in this population [7]. Epidemiological [28] and immunological [29] data have defined a role for SMA in susceptibility to invasive bacterial disease, including pneumococcal bacteremia. In keeping with this, we observed a malaria-dependent increase in risk of pneumococcal bacteremia among children with G6PD deficiency secondary to the G6PD A− variant. By contrast, carriage of the rs1050829:C allele, a WHO class 4 deficiency variant, was not associated with increased risk of SMA in Kenyan children [7] and was not associated with increased risk of pneumococcal bacteremia in this study. This effect may therefore reflect an increased risk of SMA among G6PD-deficient children. It is important to note, however, that the increased risk of malaria-induced hemolysis conferred by G6PD deficiency will result in anemia and immunological deficits that will persist beyond the index malaria infection. Indeed, in this population, recent malaria infection is a stronger predictor of nontyphoidal Salmonella bacteremia than concurrent malaria [30].

We did not see any protective effect of heterozygosity at G6PD on pneumococcal bacteremia risk. Bacteremia complicating cerebral malaria is uncommon. If the malaria-protective effects of G6PD deficiency heterozygosity are restricted to cerebral malaria, these effects would not be predicted to modify invasive bacterial disease risk, consistent with our observations. That we do not observe a protective effect of G6PD deficiency heterozygosity on pneumococcal disease risk may suggest that G6PD deficiency alleles do indeed confer specific protection against cerebral malaria, but may equally represent a lack of study power to detect such an effect. Future studies are needed to clarify protective effects of G6PD deficiency heterozygosity on malaria and bacteremia risk.

Malaria control interventions have effects on mortality among African children in excess of those attributable to the observed reduction in severe malaria syndromes [31, 32]. In this study, G6PD deficiency was estimated to account for 13% of cases of pneumococcal bacteremia, but only during a period of high malaria transmission. This highlights the significant risk of morbidity and mortality secondary to invasive pneumococcal disease conferred by G6PD deficiency in malaria-endemic settings. These data also highlight that health outcomes for children with G6PD deficiency will be improved by better malaria control interventions and that these benefits will be both directly and indirectly attributable to improved malaria control.


Our data define a role for G6PD deficiency in susceptibility to pneumococcal bacteremia in Kenyan children. We propose a model in which G6PD deficiency increases risk of malaria-induced anemia in Kenyan children, thereby increasing risk of invasive pneumococcal disease. Our data highlight the importance of G6PD deficiency as a determinant of susceptibility to infection and further highlights the potential indirect benefits of improved malaria control in malaria-endemic populations.

Availability of data and materials

Genotype and phenotype data are available via the European Genotype Archive, with the accession code EGAD00010000950.



Glucose-6-phosphate dehydrogenase


Sickle cell hemoglobin


Dihydronicotinamide adenine dinucleotide phosphate


Population attributable fraction


Severe malarial anemia


  1. Beutler E. G6PD deficiency. Blood. 1994;84:3613–36.

    Article  CAS  Google Scholar 

  2. Cappellini MD, Fiorelli G. Glucose-6-phosphate dehydrogenase deficiency. Lancet. 2008;371:64–74.

    Article  CAS  Google Scholar 

  3. Siler U, Romao S, Tejera E, Pastukhov O, Kuzmenko E, Valencia RG, et al. Severe glucose-6-phosphate dehydrogenase deficiency leads to susceptibility to infection and absent NETosis. J Allergy Clin Immunol. 2017;139:212–3.

    Article  CAS  Google Scholar 

  4. Shah SS, Macharia A, Makale J, Uyoga S, Kivinen K, Craik R, et al. Genetic determinants of glucose-6-phosphate dehydrogenase activity in Kenya. BMC Med Genet. 2014;15:93.

    Article  Google Scholar 

  5. Yoshida A, Beutler E, Motulsky AG. Human glucose-6-phosphate dehydrogenase variants. Bull World Health Organ. 1971;45:243–53.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Tishkoff SA, Varkonyi R, Cahinhinan N, Abbes S, Argyropoulos G, Destro-Bisol G, et al. Haplotype diversity and linkage disequilibrium at human G6PD: recent origin of alleles that confer malarial resistance. Science. 2001;293:455–62.

    Article  CAS  Google Scholar 

  7. Uyoga S, Ndila CM, Macharia AW, Nyutu G, Shah S, Peshu N, et al. Glucose-6-phosphate dehydrogenase deficiency and the risk of malaria and other diseases in children in Kenya: a case-control and a cohort study. Lancet Haematol. 2015;2:e437–44.

    Article  Google Scholar 

  8. Clarke GM, Rockett K, Kivinen K, Hubbart C, Jeffreys AE, Rowlands K, et al. Characterisation of the opposing effects of G6PD deficiency on cerebral malaria and severe malarial anaemia. Elife. 2017;6:531.

    Article  Google Scholar 

  9. Watson JA, Leopold SJ, Simpson JA, Day NP, Dondorp AM, White NJ. Collider bias and the apparent protective effect of glucose-6-phosphate dehydrogenase deficiency on cerebral malaria. Elife. 2019;8:253.

    Google Scholar 

  10. Reddy EA, Shaw AV, Crump JA. Community-acquired bloodstream infections in Africa:a systematic review and meta-analysis. Lancet Infect Dis. 2010;10:417–32.

    Article  Google Scholar 

  11. O'Brien KL, Wolfson LJ, Watt JP, Henkle E, Deloria-Knoll M, McCall N, et al. Burden of disease caused by Streptococcus pneumoniae in children younger than 5 years: global estimates. Lancet. 2009;374:893–902.

    Article  Google Scholar 

  12. Church J, Maitland K. Invasive bacterial co-infection in African children with Plasmodium falciparum malaria: a systematic review. BMC Med. 2014;12:31–17.

    Article  Google Scholar 

  13. Scott JAG, Berkley JA, Mwangi I, Ochola L, Uyoga S, Macharia A, et al. Relation between falciparum malaria and bacteraemia in Kenyan children: a population-based, case-control study and a longitudinal study. Lancet. 2011;378:1316–23.

    Article  Google Scholar 

  14. Uyoga S, Macharia AW, Ndila CM, Nyutu G, Shebe M, Awuondo KO, et al. The indirect health effects of malaria estimated from health advantages of the sickle cell trait. Nat Commun. 2019;10:856–7.

    Article  Google Scholar 

  15. Calis JCJ, Phiri KS, Faragher EB, Brabin BJ, Bates I, Cuevas LE, et al. Severe anemia in Malawian children. N Engl J Med. 2008;358:888–99.

    Article  CAS  Google Scholar 

  16. Kenyan Bacteraemia Study Group, Wellcome Trust Case Control Consortium 2 (WTCCC2), Rautanen A, Pirinen M, Mills TC, Rockett KA, et al. Polymorphism in a lincRNA associates with a doubled risk of pneumococcal bacteremia in Kenyan children. Am J Hum Genet. 2016;98:1092–100.

    Article  Google Scholar 

  17. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38:904–9.

    Article  CAS  Google Scholar 

  18. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.

    Article  CAS  Google Scholar 

  19. Delaneau O, Marchini J, Zagury J-F. A linear complexity phasing method for thousands of genomes. Nat Methods. 2012;9:179–81.

    Article  CAS  Google Scholar 

  20. Howie BN, Donnelly P, Marchini J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 2009;5:e1000529 Schork NJ, editor.

    Article  Google Scholar 

  21. Ho DE, Imai K, King G, Stuart EA. MatchIt: nonparametric preprocessing for parametric causal inference. J Stat Softw. 2011;42:1–28.

    Article  Google Scholar 

  22. Wakefield J. Bayes factors for genome-wide association studies: comparison with P-values. Genet Epidemiol. 2009;33:79–86.

    Article  Google Scholar 

  23. Malaria Genomic Epidemiology Network. Insights into malaria susceptibility using genome-wide data on 17,000 individuals from Africa, Asia and Oceania. Nat Commun. 2019;10:5732–19.

    Article  CAS  Google Scholar 

  24. Mogeni P, Williams TN, Fegan G, Nyundo C, Bauni E, Mwai K, et al. Age, spatial, and temporal variations in hospital admissions with malaria in Kilifi County, Kenya: a 25-year longitudinal observational study. PLoS Med. 2016;13:e1002047 Grais RF, editor.

    Article  Google Scholar 

  25. O'Meara WP, Bejon P, Mwangi TW, Okiro EA, Peshu N, Snow RW, et al. Effect of a fall in malaria transmission on morbidity and mortality in Kilifi. Kenya Lancet. 2008;372:1555–62.

    Article  Google Scholar 

  26. Berkley JA, Lowe BS, Mwangi I, Williams T, Bauni E, Mwarumba S, et al. Bacteremia among children admitted to a rural hospital in Kenya. N Engl J Med. 2005;352:39–47.

    Article  CAS  Google Scholar 

  27. Williams TN, Uyoga S, Macharia A, Ndila C, McAuley CF, Opi DH, et al. Bacteraemia in Kenyan children with sickle-cell anaemia: a retrospective cohort and case-control study. Lancet. 2009;374:1364–70.

    Article  Google Scholar 

  28. Bronzan RN, Taylor TE, Mwenechanya J, Tembo M, Kayira K, Bwanaisa L, et al. Bacteremia in Malawian children with severe malaria: prevalence, etiology, HIV coinfection, and outcome. J Infect Dis. 2007;195:895–904.

    Article  Google Scholar 

  29. Cunnington AJ, de Souza JB, Walther M, Riley EM. Malaria impairs resistance to Salmonella through heme- and heme oxygenase-dependent dysfunctional granulocyte mobilization. Nat Med. 2012;18:120–7.

    Article  CAS  Google Scholar 

  30. Brent AJ, Oundo JO, Mwangi I, Ochola L, Lowe B, Berkley JA. Salmonella bacteremia in Kenyan children. Pediatr Infect Dis J. 2006;25:230–6.

    Article  Google Scholar 

  31. Snow RW, Korenromp EL, Gouws E. Pediatric mortality in Africa: Plasmodium falciparum malaria as a cause or risk? Am J Trop Med Hyg. 2004;71:16–24.

    Article  Google Scholar 

  32. Greenwood BM, Greenwood AM, Bradley AK, Snow RW, Byass P, Hayes RJ, et al. Comparison of two strategies for control of malaria within a primary health care programme in the Gambia. Lancet. 1988;1:1121–7.

    Article  CAS  Google Scholar 

Download references


We thank the Wellcome Trust Case Control Consortium 2 for allowing us access to genotype data.

The Kenyan Bacteraemia Study Group consists of the following:

Principal Investigators: Adrian V S Hill (Chair), Thomas N Williams, J Anthony G Scott, Stephen J Chapman.

Key Personnel: Anna Rautanen, Tara C Mills, Kirk Rockett, Anne W Ndungu, Vivek Naranbhai, Alex W Macharia, Sophie Uyoga, Carolyne Ndila, Neema Mturi, Patricia Njuguna, Shebe Mohammed, James A Berkley, Isaiah Mwangi, Salim Mwarumba, Barnes S Kitsao, Brett S Lowe, Susan C Morpeth, Iqbal Khandwalla.

The Kilifi DNA extraction Group: Alex W Macharia, Sophie Uyoga, Herbert Opi, Carolyne Ndila, Emily Nyatichi, Prophet Ingosi, Barnes Kitsao, Clement Lewa, Johnstone Makale, Adan Mohamed, Kenneth Magua, Mary Njoroge, Gideon Nyutu, Ruth Mwarabu, Metrine Tendwa, and Thomas N Williams.

The Kilifi Bacteraemia Surveillance Group: Ismail Ahmed, Samuel Akech, Alexander Balo Makazi, Mohammed Bakari Hajj, Andrew Brent, Charles Chesaro, Hiza Dayo, Richard Idro, Patrick Kosgei, Kathryn Maitland, Kevin Marsh, Laura Mwalekwa, Shalton Mwaringa, Charles Newton, Mwanajuma Ngama, Allan Pamba, Norbert Peshu, Anna Seale, Alison Talbert, and Thomas N Williams. This paper is published with the permission of the Director of KEMRI.


This study makes use of data generated by the Wellcome Trust Case Control Consortium 2 project (Grant Reference 085475/B/08/Z). JJG is funded by a National Institute for Health Research (NIHR) Clinical Lectureship. TNW and JAGS are supported by Senior Research Fellowships from the Wellcome Trust (202800 and 098532, respectively). AVSH is supported by a Wellcome Trust Senior Investigator Award (HCUZZ0) and by a European Research Council advanced grant (294557). The research was supported by the Wellcome Trust Core Award Grant Number 203141/Z/16/Z with additional support from the NIHR Oxford BRC. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR, or the Department of Health and Social Care. The funders had no role in the design of the study, the collection, analysis, and interpretation of data or in the writing of the manuscript.

Author information

Authors and Affiliations




JJG, MP, and AR performed the statistical and computational analysis. TNW, JAGS, SU, SM, PN, NM, and The Kenyan Bacteraemia Study Group recruited study subjects and compiled phenotypic data. JJG, SU, AVSH, JAGS, and TNW designed the study. JJG, AVSH, JAGS, and TNW wrote the manuscript. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to James J. Gilchrist or Thomas N. Williams.

Ethics declarations

Ethics approval and consent to participate

Ethical approval for study was obtained from the Kenya Medical Research Institute (KEMRI) National Scientific Steering and Research Committees and the Oxford Tropical Research Ethics Committee (OxTREC). Following an explanation of the study, written informed consent was obtained from the parent or guardian of each child included in the study.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gilchrist, J.J., Uyoga, S., Pirinen, M. et al. Risk of pneumococcal bacteremia in Kenyan children with glucose-6-phosphate dehydrogenase deficiency. BMC Med 18, 148 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: