1. Genetics and Genomics
  2. Microbiology and Infectious Disease

The Dantu blood group prevents parasite growth in vivo: Evidence from a controlled human malaria infection study

  1. Silvia N Kariuki  Is a corresponding author
  2. Alexander W Macharia
  3. Johnstone Makale
  4. Wilfred Nyamu
  5. Stephen L Hoffman
  6. Melissa C Kapulu
  7. Philip Bejon
  8. Julian C Rayner
  9. Thomas N Williams
  10. On behalf of for the CHMI-SIKA Study Team
  1. Centre for Geographic Medicine Research (Coast), KEMRI-Wellcome Trust Research Programme, Kenya
  2. Sanaria Incorporate, United States
  3. Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, United Kingdom
  4. Cambridge Institute for Medical Research, University of Cambridge, United Kingdom
  5. Institute for Global Health Innovation, Department of Surgery and Cancer, Imperial College London, United Kingdom
https://doi.org/10.7554/eLife.83874
Share this article
Cite this article
  1. Silvia N Kariuki
  2. Alexander W Macharia
  3. Johnstone Makale
  4. Wilfred Nyamu
  5. Stephen L Hoffman
  6. Melissa C Kapulu
  7. Philip Bejon
  8. Julian C Rayner
  9. Thomas N Williams
  10. On behalf of for the CHMI-SIKA Study Team
(2023)
The Dantu blood group prevents parasite growth in vivo: Evidence from a controlled human malaria infection study
eLife 12:e83874.
https://doi.org/10.7554/eLife.83874
  2. Download BibTeX
  3. Download .RIS

Abstract

Background:

The long co-evolution of Homo sapiens and Plasmodium falciparum has resulted in the selection of numerous human genetic variants that confer an advantage against severe malaria and death. One such variant is the Dantu blood group antigen, which is associated with 74% protection against severe and complicated P. falciparum malaria infections in homozygous individuals, similar to that provided by the sickle haemoglobin allele (HbS). Recent in vitro studies suggest that Dantu exerts this protection by increasing the surface tension of red blood cells, thereby impeding the ability of P. falciparum merozoites to invade them and reducing parasite multiplication. However, no studies have yet explored this hypothesis in vivo.

Methods:

We investigated the effect of Dantu on early phase P. falciparum (Pf) infections in a controlled human malaria infection (CHMI) study. 141 sickle-negative Kenyan adults were inoculated with 3.2 × 103 aseptic, purified, cryopreserved Pf sporozoites (PfSPZ Challenge) then monitored for blood-stage parasitaemia for 21 days by quantitative polymerase chain reaction (qPCR)analysis of the 18S ribosomal RNA P. falciparum gene. The primary endpoint was blood-stage P. falciparum parasitaemia of ≥500/μl while the secondary endpoint was the receipt of antimalarial treatment in the presence of parasitaemia of any density. On study completion, all participants were genotyped both for Dantu and for four other polymorphisms that are associated with protection against severe falciparum malaria: α+-thalassaemia, blood group O, G6PD deficiency, and the rs4951074 allele in the red cell calcium transporter ATP2B4.

Results:

The primary endpoint was reached in 25/111 (22.5%) non-Dantu subjects in comparison to 0/27 (0%) Dantu heterozygotes and 0/3 (0.0%) Dantu homozygotes (p=0.01). Similarly, 49/111 (44.1%) non-Dantu subjects reached the secondary endpoint in comparison to only 7/27 (25.9%) and 0/3 (0.0%) Dantu heterozygotes and homozygotes, respectively (p=0.021). No significant impacts on either outcome were seen for any of the other genetic variants under study.

Conclusions:

This study reveals, for the first time, that the Dantu blood group is associated with high-level protection against early, non-clinical, P. falciparum malaria infections in vivo. Learning more about the mechanisms involved could potentially lead to new approaches to the prevention or treatment of the disease. Our study illustrates the power of CHMI with PfSPZ Challenge for directly testing the protective impact of genotypes previously identified using other methods.

Funding:

The Kenya CHMI study was supported by an award from Wellcome (grant number 107499). SK was supported by a Training Fellowship (216444/Z/19/Z), TNW by a Senior Research Fellowship (202800/Z/16/Z), JCR by an Investigator Award (220266/Z/20/Z), and core support to the KEMRI-Wellcome Trust Research Programme in Kilifi, Kenya (203077), all from Wellcome. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. For the purpose of Open Access, the authors have applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.

Clinical trial number:

NCT02739763

Editor's evaluation

The large genetic association studies conducted in East Africa have shown that the Dantu blood group provides substantial protection against severe malaria since it increases the surface tension of red blood cells making it harder for malaria parasites to invade. In this important work, the authors show that parasite growth is indeed restricted in vivo by testing this hypothesis in adult Kenyan volunteers infected with P. falciparium under careful monitoring. They were able to show convincingly that indeed, parasite growth was reduced amongst Dantu adults.

https://doi.org/10.7554/eLife.83874.sa0

Introduction

Plasmodium falciparum malaria has been the pre-eminent cause of child morbidity and mortality in the tropics and sub-tropics for much of the last 5000 y. As a consequence, it has had a substantial impact on the human genome through the positive selection of multiple polymorphisms that confer a survival advantage against the disease (Williams, 2017). The best studied affect the biology of red blood cells (RBCs), which host malaria parasites for most of their life cycle in humans, the rs334 A>T βs sickle mutation in HBB (Williams, 2016), α-thalassaemia (Williams et al., 2005), and blood group O (Fry et al., 2008) all being important examples.

Recently, we identified a new variant which is associated with high-level protection against severe P. falciparum malaria to a degree that is close to that of sickle cell trait (HbAS), the strongest malaria-protective condition yet described (Malaria Genomic Epidemiology Network, 2014). The rare Dantu blood group antigen, which results from a genetic rearrangement within the glycophorin (GYP) cluster, was shown to confer 74% protection against severe malaria in homozygous individuals (Band et al., 2015; Ndila et al., 2018). Subsequent in vitro studies have suggested that this protection is explained by the resistance of Dantu RBCs to invasion by P. falciparum merozoites (Kariuki et al., 2020), thereby preventing infections from progressing to become severe or ultimately fatal. However, this hypothesis has not been tested directly in vivo to date.

In this study, we have investigated the impact of the Dantu blood group on in vivo P. falciparum parasite growth and clinical disease progression through a controlled human malaria infection (CHMI) study with aseptic, purified, cryopreserved P. falciparum sporozoites (PfSPZ Challenge) conducted in semi-immune Kenyan adults. To the best of our knowledge, this is the first time that CHMI has been used to directly explore the impact of Dantu on parasite growth in vivo.

Methods

Study design and population

The primary aim of the Kenya CHMI study was to investigate the impacts of naturally acquired immunity on early-phase malaria infections (Kapulu et al., 2018; Kapulu et al., 2021). Briefly, 161 healthy adult volunteers living in areas of varying malaria transmission were inoculated by direct venous inoculation (DVI) with 3.2 × 103 P. falciparum sporozoites (PfSPZ) of Sanaria PfSPZ Challenge (NF54) (Roestenberg et al., 2013; Mordmüller et al., 2015). With a sample size of 161 individuals, we had 80% power to detect a single variable with an effect size (r2) of 0.3 that accounts for 15% of the variability in parasite growth, as previously described (Kapulu et al., 2018; Kapulu et al., 2021). Because of its major impacts on both malaria susceptibility (Allison, 1954) and disease progression (Taylor et al., 2012), and in view of results from a previous CHMI with PfSPZ Challenge conducted in Gabon (Lell et al., 2018), recruitment was restricted to those who were negative for both sickle cell trait and disease. After inoculation, venous blood samples were collected twice daily from days 7–14, and then once every day from day 15 until the end of the experiment on day 21, and screened for parasitaemia by quantitative polymerase chain reaction (qPCR) analysis of the P. falciparum 18S ribosomal RNA gene. For each participant, the endpoint was considered met, and anti-malarial treatment administered, when a threshold of 500 P. falciparum parasites/µl was reached. Participants were treated earlier if signs and symptoms were observed in association with blood film positivity at any parasite density, and on day 21 post-inoculation regardless of outcome. While qPCR was used as the primary measurement of parasitaemia because it is much more sensitive than microscopy at low density (Bejon et al., 2006), thick film blood smears were also performed as an additional precaution, and participants were treated if they became blood film positive at any density, an approach that accords with that used in other CHMI studies (Sheehy et al., 2012; Murphy et al., 2012; Hodgson et al., 2014; Kamau et al., 2014; Hodgson et al., 2015; Seilie et al., 2019; Salkeld et al., 2022). Participants were recruited during 2016, 2017, and 2018 into three successive cohorts from three different malaria transmission zones: Kilifi North (no- to low-transmission) and Kilifi South (moderate transmission), both on the coast, and Ahero (moderate to high transmission) in Western Kenya. We tested for antimalarial drugs as described in a previous publication (Kapulu et al., 2021). Briefly, antimalarial drugs were measured retrospectively in all volunteers using samples collected on both the day before the challenge and 8 days after challenge. Plasma samples were tested by liquid chromatography-tandem mass spectrometry in two independent laboratories. Sulfadoxine, pyrimethamine, and chloroquine levels were measured at the Strathmore University in Nairobi, Kenya, while artemether and dihydroartemisinin concentrations were measured at the Mahidol Oxford Tropical Medicine Research Unit in Bangkok, Thailand. The study was conducted at the KEMRI-Wellcome Trust Research Programme in Kilifi, Kenya, and was registered on ClinicalTrials.gov (NCT02739763).

Genotyping for Dantu and other malaria-protective variants

Whole blood samples were collected at the point of recruitment into tubes containing EDTA and stored at –80°C pending batch processing at the end of the study. Genomic DNA was extracted following the manufacturer’s instructions using a QIAmp 96 DNA QIAcube HT kit on a QIAcube HT System (QIAGEN, Manchester, UK). Genotyping for Dantu was performed using ABI TaqMan SNP genotyping Assays-by-Design primers and probes on an ABI 7900HT PCR machine, as previously described (Kariuki et al., 2009). Dantu genotypes were inferred from the rs186873296 FREM3 allele, which is in strong linkage disequilibrium with the Dantu structural rearrangement (Band et al., 2015). The following primer sequence was used for the rs186873296 SNP: ATGTGAAGAAGCTGGGAACCCTGTC[A/G]TACAAGAAATGACAAAGAAAGCTT, with A being the reference allele. For comparative purposes, we also genotyped participants for a range of other polymorphisms that have been reproducibly associated with protection from severe and complicated malaria in other studies. We typed for the common African form of G6PD deficiency, caused by the G6PD c.202T mutation (Clarke et al., 2017), the blood group O mutation in ABO (Rowe et al., 2007) and the rs4951074 allele in ATP2B4 (Malaria Genomic Epidemiology Network, 2014) using TaqMan SNP genotyping assays, and for the -α3.7I deletional form of α+-thalassaemia by gap PCR (Wambua et al., 2006). To aid ease of comparison across all variants tested, the genotype groups are categorised as ‘Homozygous reference’ for individuals with two copies of the reference allele, ‘Heterozygous’ for individuals with one copy of the reference allele and one copy of the derived allele, and ‘Homozygous derived’ for individuals with two copies of the derived allele.

Statistical analysis

We considered three distinct outcomes, each capturing a different aspect of the parasitological and clinical progress of malaria infections: (a) whether infections progressed to reach the pre-defined treatment threshold of 500 parasites/µl; (b) whether or not participants received malaria treatment for either clinical or parasitological reasons; (c) the time from inoculation to treatment in participants who did receive treatment. We made (a) the primary endpoint for this analysis because the hypothesis we were testing concerned in vivo parasite growth rather than susceptibility to symptoms. We conducted between-genotype comparisons both by univariate analysis and by multivariate analysis with adjustment for other malaria-protective genes, anti-schizont antibody concentration, and location of residence, of which the latter two were significantly associated with the same outcomes in an earlier analysis of the same cohort (Kapulu et al., 2022). We used the Fisher’s exact test using the fmsb package (version 0.7.3) in R to test for any differences in the proportions of individuals that reached the pre-defined treatment threshold of 500 parasites/μl and to test for any differences in the proportions of individuals that required treatment across genotype groups. We also compared the treatment outcomes for the heterozygous- and homozygous-derived genotypes, individually, to the homozygous reference genotype by pairwise analysis.

To test for associations between each variant genotype and the treated or untreated categorical outcome, we used multivariate logistic regression using the glm function in the stats package (version 3.6.2) in R using additive models for each variant, where each variant genotype was coded as zero, one, or two copies of the derived allele. Each multivariable model adjusted for other malaria-protective variants, anti-schizont antibody concentration, and location of residence.

We used the Kruskal–Wallis test (stats package version 3.6.2) and the Dunn’s test (FSA package version 0.9.3) to investigate between-group differences in maximum parasitaemia. Finally, we compared time to treatment using Kaplan–Meier survival curves with univariate comparisons across genotype groups performed using the Log-Rank test, and Cox regression models for multivariate analyses, using the survival package (version 3.2.13) in R. All statistical analyses were performed using R V3.6.2 (R Development Core Team, 2017).

Results

Of 161 volunteers recruited to the study, 19 were excluded, either because non-CHMI parasite strains were detected in post-CHMI samples (suggesting the presence of coincidental natural infections) (n = 7), or because antimalarial drugs were detected in pre-CHMI samples (n = 12) (Figure 1) (Kapulu et al., 2021). After these exclusions, data from 142 individuals contributed to the current analysis. Genotyping revealed that 30 of these individuals were either heterozygous or homozygous for the Dantu allele. We did not get Dantu genotype data on one individual due to poor quality of the DNA sample; therefore, Dantu genotype data on 141 out of the 142 participants was used in the downstream analysis (Figure 1).

Figure 1
Download asset Open asset
Study design and participant recruitment.

The Dantu variant protects against P. falciparum growth in vivo

While infections progressed to the point of reaching the pre-determined treatment threshold of 500 parasites/µl in 25/142 (17.6%) of all volunteers, the proportions varied markedly by Dantu genotype. None of the thirty (0.0%) Dantu carriers reached this threshold in comparison to 25/111 (22.5%) of the non-Dantu individuals (p=0.01) (Figure 2, Table 1). The difference in the proportion of Dantu heterozygous (0/27; 0%) and non-Dantu individuals reaching the treatment threshold was strongly significant (p=0.004) but did not reach statistical significance in the case of Dantu homozygotes 0/3 (0.0%) because of the small number of individuals in this group.

Figure 2
Download asset Open asset
The impact of human genotype on parasite growth.

Following inoculation of volunteers with Pf sporozoites, parasitaemia was monitored by quantitative PCR (y-axis) over the full duration of the study (x-axis). The different panels indicate the specific variants that were studied. The red dots indicate only the individuals that exhibited febrile symptoms and met treatment criteria. The different genotype groups for all the variants are categorised as ‘Homozygous reference’ for individuals with two copies of the reference allele (red lines), ‘Heterozygous’ for individuals with one copy of the reference allele and one copy of the derived allele (green lines), and ‘Homozygous derived’ for individuals with two copies of the derived allele (blue lines). αα/αα, no α-thalassaemia; −α/αα, heterozygous α-thalassaemia; −α/−α, homozygous α-thalassaemia. * For G6PD, male and female were combined as C/CC = normal (wild type hemizygous males and homozygous females), CT = carrier females, and T/TT = G6PD-deficient hemizygous males and homozygous females. All volunteers were typed for all variants, and any one individual may carry a mixture of genotypes – the potential confounding effect of this was controlled for in multivariate analysis. The tables adjacent to each plot show the results from Fisher’s exact tests investigating differences in the proportion of participants that reached the pre-defined treatment threshold of 500 parasites/μl (n) compared to the total number within each genotype category (N).

Figure 2—source data 1

Related to Figure 2.

https://cdn.elifesciences.org/articles/83874/elife-83874-fig2-data1-v1.txt
Figure 2—source data 2

Related to Figure 2 table.

https://cdn.elifesciences.org/articles/83874/elife-83874-fig2-data2-v1.xlsx
Table 1
The proportion of participants reaching the pre-defined treatment parasitaemia threshold of 500 parasites/μl by genotype category.
VariantGenotypen/N%p value overallp value homozygous reference vs. heterozygousp value homozygous reference vs. homozygous derived
Dantu rs186873296Non-Dantu (AA)25/11122.50.010.0041
Heterozygous (AG)0/270.0
Dantu homozygous (GG)0/30.0
G6PD +202 rs1050828Homozygous reference (C/CC)20/11018.20.5050.7390.463
Heterozygous (CT)4/1723.5
Homozygous derived (T/TT)1/156.7
ABO rs8176719Non-O13/6420.30.5170.517-
O12/7516.0
α-thalassaemiaHomozygous reference (αα/αα)6/4812.50.4080.3280.29
Heterozygous (−α/αα)14/7020.0
Homozygous derived (−α/−α)5/2123.8
ATP2B4 rs4951074Homozygous reference (GG)11/6217.70.87610.749
Heterozygous (AG)11/5719.3
Homozygous derived (AA)3/2313.0

  1. n = the number of participants that reached the pre-defined treatment threshold of 500 parasites/μl and were treated; N = the total number within each genotype category; αα/αα, no α-thalassaemia; −α/αα, heterozygous α-thalassaemia; −α/−α, homozygous α-thalassaemia. * For G6PD, male and female were combined as C/CC = normal (wild type hemizygous males and homozygous females), CT = carrier females, and T/TT = G6PD -deficient hemizygous males and homozygous females. We used the Fisher’s exact test to investigate differences in the proportions of individuals that reached the pre-defined treatment threshold of 500 parasites/μl, both for global comparisons across genotype groups, and for separate comparisons between genotype pairs, where the proportion of individuals that reached the treatment threshold of 500 parasites/μl in the heterozygous- and homozygous -derived genotypes were compared to the homozygous reference genotype.

Table 1—source data 1

Related to Table 1.

https://cdn.elifesciences.org/articles/83874/elife-83874-table1-data1-v1.txt

Fewer Dantu carriers required malaria treatment

Although a threshold of 500 parasites/μl was considered the primary endpoint of the study and triggered the administration of antimalarial treatment per-protocol, some participants developed symptoms and were therefore treated before reaching this parasitaemia threshold. Only one quarter (7/27; 25.9%) of the Dantu heterozygotes and none (0/3; 0.0%) of the Dantu homozygotes (Figure 3, Table 2) received antimalarial treatment in comparison to 49/111 (44.1%) of the non-Dantu individuals. While this did not reach statistical significance on univariate analysis, we carried out multivariate regression analysis with the treated or untreated status as the dependent variable, and the Dantu variant as well as other genetic variants (as described below), anti-schizont antibody and location of residence as the independent variables. This analysis revealed that Dantu-carrying subjects overall were administered antimalarial treatment 83% less frequently (OR 0.17; 95% CI 0.04–0.55; p=0.007) than non-Dantu individuals (Table 3). In order to address the core issue of whether prior immunity was a confounder in our analysis, we used measurements of antibodies to whole schizont extract as a proxy indicator of transmission setting or ‘malaria exposure’ in our multivariate analyses. We compared anti-schizont antibody levels across Dantu genotype groups and found no differences (p=0.659) (Figure 3—figure supplement 1).

Figure 3 with 1 supplement see all
Download asset Open asset
The impact of each gene variant on the requirement for malaria treatment.

The proportion of individuals in each genotype category that required treatment over the course of the controlled human malaria infection (CHMI) study is shown on the y-axis. The number of treated individuals out of the total number in each genotype group is given in parenthesis above the bar graphs, while the p values from the Fisher’s exact tests comparing the differences in proportions of individuals that required treatment across genotype groups are also given above the bar graphs.

Figure 3—source data 1

Related to Figure 2.

https://cdn.elifesciences.org/articles/83874/elife-83874-fig3-data1-v1.txt
Table 2
The numbers and frequencies of individuals who received treatment before day 21 by genotypic category.
VariantGenotypen/N%Overall p valueHomozygous reference vs. heterozygousHomozygous reference vs. homozygous derived
Dantu rs186873296Non-Dantu (AA)49/11144.10.100.130.26
Heterozygous (AG)7/2725.9
Dantu homozygous (GG)0/30.0
G6PD +202 rs1050828*Homozygous reference (C/CC)46/11041.80.281.000.16
Heterozygous (CT)7/1741.2
Homozygous derived (T/TT)3/1520.0
ABO rs8176719Non-O28/6443.80.39-0.39
O27/7536.0
α-thalassaemiaHomozygous reference (αα/αα)19/4839.60.790.850.79
Heterozygous (−α/αα)30/7042.9
Homozygous derived (−α/−α)7/2133.3
ATP2B4 rs4951074Homozygous reference (GG)26/6241.90.611.000.45
Heterozygous (AG)23/5740.4
Homozygous derived (AA)7/2330.4
Table 2—source data 1

Related to Table 2.

https://cdn.elifesciences.org/articles/83874/elife-83874-table2-data1-v1.txt
Table 3
The association between each gene variant and the requirement for treatment after challenge.
Overall comparison across genotype groupsHomozygous reference vs. heterozygousHomozygous reference vs. homozygous derived
VariantOdds ratios95% CIp-valueOdds ratios95% CIp valueOdds ratios95% CIp value
Dantu rs1868732960.170.04–0.550.0070.200.04–0.830.0390NA – Inf0.990
G6PD +202 rs10508280.600.28–1.190.1571.710.41–6.590.4420.170.02–0.970.074
ABO rs81767190.400.15–1.030.0640.540.18–1.560.259---
α-thalassaemia0.910.46–1.730.7660.970.31–3.040.9570.870.18–3.220.758
ATP2B4 rs49510740.840.43–1.630.6190.530.17–1.580.2660.570.10–2.580.491
Table 3—source data 1

Related to Table 3.

https://cdn.elifesciences.org/articles/83874/elife-83874-table3-data1-v1.txt

Peak parasitaemias were lower in Dantu-carrying participants

The proportion of participants who became PCR-positive at any parasitaemia was similar at 86/111 (77.5%) in the non-Dantu and 20/27 (74.1%) and 2/3 (66.7%) among Dantu heterozygotes and homozygotes, respectively (p=0.745) (Figure 4). However, maximum parasitaemias were considerably higher in non-Dantu than in Dantu-carrying individuals. Peak parasitaemias reached 9694 parasites/μl in the non-Dantu group in comparison with 411 parasites/μl in the Dantu heterozygotes and only 3 parasites/μl among the Dantu homozygotes (non-Dantu vs. Dantu heterozygotes p=0.028; non-Dantu vs. Dantu homozygotes p=0.141; and non-Dantu vs. Dantu heterozygotes and homozygotes combined p=0.009) (Figure 4). Similarly, the median parasitaemia among those who did become PCR-positive was 112 parasites/μl in the non-Dantu, 13 parasites/μl in the Dantu heterozygous, and 2 parasites/μl in the Dantu homozygous groups, respectively, although these differences did not reach statistical significance (non-Dantu vs. Dantu heterozygotes p=0.108; non-Dantu vs. Dantu homozygotes p=0.256; and non-Dantu vs. combined Dantu heterozygotes and homozygotes p=0.068) (Figure 4).

Figure 4
Download asset Open asset
Peak parasitaemias were lower in Dantu variant carriers.

Maximum parasitaemia values for individuals across Dantu genotype groups, with dashed line indicating the treatment threshold of 500 parasites/μl. The table below the figure shows the numbers and frequencies of individuals in each genotype category that were PCR-positive over the course of the controlled human malaria infection (CHMI) study. n = the number of participants that were PCR-positive; N = the total number within the genotype category. Statistical comparisons of proportions of PCR-positive individuals across genotype groups and pairwise comparisons between genotype groups were performed using the Fisher’s exact test. Statistical comparisons of maximum and median parasitaemia between genotype groups were performed using the Kruskal–Wallis test, and post-hoc Dunn’s test for pairwise differences between the genotype groups.

Figure 4—source data 1

Related to Figure 4.

https://cdn.elifesciences.org/articles/83874/elife-83874-fig4-data1-v1.txt
Figure 4—source data 2

Related to Figure 4 table.

https://cdn.elifesciences.org/articles/83874/elife-83874-fig4-data2-v1.xlsx

Time to treatment was significantly longer in Dantu-carrying than non-Dantu individuals

Among the participants who did receive treatment, the time to treatment was significantly longer among Dantu-carrying than non-Dantu individuals. While the univariate comparisons across genotype groups, performed using the Log-Rank test in the Kaplan–Meier survival curves, did not reach statistical significance (Figure 5a), the multivariate Cox regression analysis with adjustments for other malaria-protective variants, anti-schizont antibody concentration, and location of residence showed that time to treatment was significantly longer among Dantu-carrying than non-Dantu individuals, the overall hazard ratio being 0.39 (CI 0.17–0.87; p=0.022) (Figure 5b). A dose-dependent effect of the Dantu genotype was also seen as none of the three Dantu homozygotes required treatment and the time to treatment was significantly longer in Dantu heterozygous than in non-Dantu individuals (HR = 0.41, p=0.042) (Figure 5b).

Figure 5
Download asset Open asset
Time to treatment was longer in Dantu variant carriers.

The impact of each gene variant on time to treatment was analysed by (a) Kaplan–Meier survival curves, with univariate comparisons across genotype groups performed using the Log-Rank test and (b) multivariate Cox regression models, with each variant genotype coded as zero, one, or two copies of the homozygous derived allele in an additive model, adjusting for the other four malaria-protective variants, anti-schizont antibody concentration, and location of residence. Pairwise analysis compared the time to treatment in the heterozygous- and homozygous-derived genotypes to the homozygous reference genotype.

Figure 5—source data 1

Related to Figure 5a.

https://cdn.elifesciences.org/articles/83874/elife-83874-fig5-data1-v1.txt
Figure 5—source data 2

Related to Figure 5b.

https://cdn.elifesciences.org/articles/83874/elife-83874-fig5-data2-v1.txt
Figure 5—source data 3

Related to Figure 5b table.

https://cdn.elifesciences.org/articles/83874/elife-83874-fig5-data3-v1.xlsx

No significant impacts were seen with regard to any of the outcomes under study for any of the remaining RBC polymorphisms

While our primary analysis focused on the Dantu genotype, individuals in malaria-endemic regions often carry more than one malaria protection-associated genotype (Ndila et al., 2018). We therefore typed individuals for α-thalassaemia, blood group O, G6PD deficiency and ATP2B4 alleles, and carried out the same set of analysis as above for each genotype. Each of these genotypic groups also included, by random chance, individuals of different genotypes for the other variants of interest. As the genotyping was conducted at the end of the study, and the study included only a relatively small number of individuals overall, it was not possible to limit each group to those who were reference homozygotes for all the other variants of interest. Instead, we used multivariate analyses to check that any differences seen on univariate analysis were not explained by the presence of other variants or confounders. No significant differences were seen in any of the study outcomes between the different genotype groups (Figures 2, 3 and 5, Tables 1–3). As noted above, the impact of Dantu on malaria treatment was independent of the presence or absence of these other genetic variants in multivariate analysis.

Discussion

Through the analysis of data from a CHMI study with PfSPZ Challenge injection conducted in Kenyan adults, we have shown that the Dantu blood group is associated with prevention of parasite growth in vivo. While more than 20% of Dantu-negative volunteers developed bloodstream malaria infections that reached a pre-defined threshold of 500 parasites/μl following controlled PfSPZ administration, this threshold was not reached by any of the 30 Dantu-positive volunteers. This is the first time that Dantu has been shown to protect against early-stage malaria infections. Our study suggests that Dantu protects against severe and complicated malaria (Ndila et al., 2018; Band et al., 2015) by preventing the disease from becoming established in its earliest phase. This is consistent with recent in vitro observations that have demonstrated a link between Dantu genotype and susceptibility to red blood cell invasion by P. falciparum merozoites, which we previously predicted would lead to reduced parasite growth in vivo (Kariuki et al., 2020). That mechanistic study revealed that the Dantu genotype protects red blood cells from invasion by increasing their surface tension, which reduces the ability of merozoites to deform their surface and hence productively invade. Critically, this tension difference was greater in Dantu homozygotes than heterozygotes, as was the reduction in invasion efficiency, supporting a dose-dependent protective effect. This dose-dependency was similarly reflected in our current in vivo study, as while 44.1% of the Dantu-negative volunteers developed symptoms that precipitated malaria treatment, this occurred in only 25.9% of Dantu heterozygous and in none of the three (0%) Dantu homozygous volunteers. Similarly, the maximum parasitaemia observed was considerably higher in the non-Dantu than in the Dantu heterozygous and homozygous participants.

In contrast to CHMI studies in malaria-naïve populations, where most participants typically develop clinical malaria (Church et al., 1997), less than one-fifth of the participants in this study reached a threshold of >500 parasites/μl. This is probably because most of the participants in our study were residents of malaria-endemic communities and will therefore have been partially immune to the disease. Indeed, in previous analyses of data from the same study we have shown that infection outcome was attributable to the degree of prior exposure as estimated by the titre of anti-schizont antibodies (Kapulu et al., 2022), as well as other measures of exposure including the antibody-dependent phagocytosis of both ring-infected and uninfected erythrocytes from parasite cultures (Musasia et al., 2022) and the breadth of antibodies to P. falciparum Variant Surface Antigens (Kimingi et al., 2022). As such, Dantu genotype was therefore clearly not the only factor at play in determining the clinical outcome in this study. However, our multivariate analysis adjusted for these factors (Kapulu et al., 2022) as well as other malaria-protective genetic variants, and the significant association in that analysis underscores the strong protective effect conferred by Dantu. There were also no differences observed in anti-schizont antibody levels across Dantu genotype groups, suggesting that differences in pre-existing anti-malaria immunity between Dantu and non-Dantu cannot explain the differences seen in this study.

The protective impact of Dantu against in vivo parasite growth was in stark contrast to that of the other genetic factors under study. Although consistent evidence has been found for protective effects against severe malaria by G6PD deficiency, blood group O, the rs4951074 allele in ATP2B4 and α+-thalassaemia in numerous previous studies (Malaria Genomic Epidemiology Network, 2014; Taylor et al., 2012), none of these polymorphisms had any significant impacts on any of the outcomes under investigation in this study. This is probably because, unlike Dantu, none of these conditions have a clearly established impact on merozoite red cell invasion, but instead influence the progress of malaria once the disease has become established. For example, recent studies have shown that α+-thalassaemia has no effect on either red blood cell invasion (Williams et al., 2002) or the development of uncomplicated malaria (Taylor et al., 2012; Wambua et al., 2006), but protects instead against the development of severe and complicated disease through mechanisms that include the reduced expression of red cell surface antigens that result in cytoadhesion (Krause et al., 2012; Opi et al., 2014). Similarly, it had no apparent impact on infectivity in a previous, smaller, CHMI study using PfSPZ Challenge, conducted in Tanzania (Shekalaghe et al., 2014). Because clinical guidelines mean that controlled human malaria challenge studies only ever reach relatively low-density parasitaemias, they may not adequately capture the impacts of genetic factors that influence the later, more severe stages of malaria disease. However, this study shows that they can be helpful in pointing towards the pathways leading to such outcomes for further study by other methods.

In conclusion, this study reveals the power of CHMI studies to deconvolute the malaria-protective effects of naturally occurring human genetic variants, establishes for the first time that the Dantu blood group provides strong protection against in vivo parasite growth, and emphasises the potential of Dantu-phenocopying interventions to limit P. falciparum growth in vivo.

Appendix 1

Members of the CHMI-SIKA Study Team

  • Centre for Geographic Medicine Research (Coast), Kenya Medical Research Institute-Wellcome Trust Research Programme, Kilifi, Kenya: Abdirahman I Abdi, Philip Bejon, Zaydah de Laurent, Mainga Hamaluba, Domtila Kimani, Rinter Kimathi, Kevin Marsh, Sam Kinyanjui, Khadija Said Mohammed, Moses Mosobo, Janet Musembi, Jennifer Musyoki, Michelle Muthui, Jedidah Mwacharo, Kennedy Mwai, Joyce M Ngoi, Omar Ngoto, Patricia Njuguna, Irene Nkumama, Francis Ndungu, Dennis Odera, Donwilliams Omuoyo, Faith Osier, Edward Otieno, Jimmy Shangala, James Tuju, Juliana Wambua, Thomas N Williams

  • Sanaria Inc, Rockville, MD, United States: Yonas Abebe, Peter F Billingsley, Stephen L Hoffman, Eric R James, Thomas L Richie, B Kim Lee Sim

  • Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University Oxford, Oxford, United Kingdom: Sam Kinyanjui, Kevin Marsh

  • Department of Pathology, University of Cambridge, Cambridge, United Kingdom: Peter C Bull

  • Pwani University, P. O. Box 195-80108, Kilifi, Kenya: Sam Kinyanjui, Cheryl Kivisi

  • Epidemiology and Biostatistics Division, School of Public Health, University of the Witwatersrand, Johannesburg, South Africa: Kennedy Mwai

  • Centre for Infectious Diseases, Heidelberg University Hospital, Heidelberg, Germany: Irene Nkumama, Dennis Odera, Faith Osier

  • Center for Research in Therapeutic Sciences, Strathmore University, Nairobi, Kenya: Bernhards Ogutu, Fredrick Olewe, John Ong’echa

  • Center for Research in Therapeutic Sciences, Strathmore University, Nairobi, Kenya: Bernhards Ogutu, Fredrick Olewe, John Ong’echa

Data availability

All data generated or analysed during this study are included in the manuscript and supporting file; Source Data files have been provided for Figures 2, 3, 4, and 5 and Tables 1, 2 and 3.

References

  24. Software
    1. R Development Core Team
    (2017) R: A language and environment for statistical computing
    R Foundation for Statistical Computing, Vienna, Austria.
    https://www.r-project.org
  36. Book
    1. Williams TN
    (2017)
    Host Genetics
    In: Gaur G, editors. Advances in Malaria Research. John Wiley & Sons, Inc. pp. 465–494.

Article and author information

Author details

  1. Silvia N Kariuki

    Centre for Geographic Medicine Research (Coast), KEMRI-Wellcome Trust Research Programme, Kilifi, Kenya
    Contribution
    Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing - original draft, Writing – review and editing
    For correspondence
    SNKariuki@kemri-wellcome.org
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0801-5285

  2. Alexander W Macharia

    Centre for Geographic Medicine Research (Coast), KEMRI-Wellcome Trust Research Programme, Kilifi, Kenya
    Contribution
    Investigation, Methodology
    Competing interests
    No competing interests declared

  3. Johnstone Makale

    Centre for Geographic Medicine Research (Coast), KEMRI-Wellcome Trust Research Programme, Kilifi, Kenya
    Contribution
    Investigation, Methodology
    Competing interests
    No competing interests declared

  4. Wilfred Nyamu

    Centre for Geographic Medicine Research (Coast), KEMRI-Wellcome Trust Research Programme, Kilifi, Kenya
    Contribution
    Investigation
    Competing interests
    No competing interests declared

  5. Stephen L Hoffman

    Sanaria Incorporate, Rockville, United States
    Contribution
    Conceptualization, Resources, Methodology, Writing – review and editing
    Competing interests
    SLH and members of the CHMI-SIKA Study Team (YA, PFB, ERJ, TLR, and BKLS; see Author Information) are salaried, full-time employees of Sanaria, the manufacturer of Sanaria PfSPZ Challenge. The authors have no additional financial interests

  6. Melissa C Kapulu

    1. Centre for Geographic Medicine Research (Coast), KEMRI-Wellcome Trust Research Programme, Kilifi, Kenya
    2. Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
    Contribution
    Conceptualization, Resources, Investigation, Methodology, Project administration, Writing – review and editing
    Competing interests
    Melissa C Kapulu has received grants from UKRI MRC UK and Wellcome Trust and has received travel expenses for the HIC-VAC Annual Meeting, Greenwood Africa Prize and WHO working group meeting on challenge studies. The author has no other competing interests to declare

  7. Philip Bejon

    1. Centre for Geographic Medicine Research (Coast), KEMRI-Wellcome Trust Research Programme, Kilifi, Kenya
    2. Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
    Contribution
    Conceptualization, Resources, Formal analysis, Supervision, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared

  8. Julian C Rayner

    Cambridge Institute for Medical Research, University of Cambridge, Cambridge, United Kingdom
    Contribution
    Conceptualization, Supervision, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9835-1014

  9. Thomas N Williams

    1. Centre for Geographic Medicine Research (Coast), KEMRI-Wellcome Trust Research Programme, Kilifi, Kenya
    2. Institute for Global Health Innovation, Department of Surgery and Cancer, Imperial College London, London, United Kingdom
    Contribution
    Conceptualization, Resources, Supervision, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4456-2382

Funding

Wellcome Trust (107499)

  • Melissa C Kapulu
  • Philip Bejon

Wellcome Trust (216444/Z/19/Z)

  • Silvia N Kariuki

Wellcome Trust (202800/Z/16/Z)

  • Thomas N Williams

Wellcome Trust (220266/Z/20/Z)

  • Julian C Rayner

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. For the purpose of Open Access, the authors have applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.

Ethics

Clinical trial registration: The study is registered on ClinicalTrials.gov (NCT02739763).

The study was approved by both the KEMRI Scientific and Ethics Review Unit (protocol KEMRI/SERU/CGMR-C/029/3190) in Kenya, and the University of Oxford Tropical Research Ethics Committee (OxTREC; protocol 2-16) in the UK. The study was conducted based on good clinical practice (GCP) following the principles of the Declaration of Helsinki.

Version history

  1. Preprint posted: September 22, 2022 (view preprint)
  2. Received: September 30, 2022
  3. Accepted: May 11, 2023
  4. Version of Record published: June 13, 2023 (version 1)

Copyright

© 2023, Kariuki et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 703
    views
  • 54
    downloads
  • 2
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Silvia N Kariuki
  2. Alexander W Macharia
  3. Johnstone Makale
  4. Wilfred Nyamu
  5. Stephen L Hoffman
  6. Melissa C Kapulu
  7. Philip Bejon
  8. Julian C Rayner
  9. Thomas N Williams
  10. On behalf of for the CHMI-SIKA Study Team
(2023)
The Dantu blood group prevents parasite growth in vivo: Evidence from a controlled human malaria infection study
eLife 12:e83874.
https://doi.org/10.7554/eLife.83874

Share this article

https://doi.org/10.7554/eLife.83874

Categories and tags

Research organisms

Further reading

    1. Genetics and Genomics

    Mono-methylated histones control PARP-1 in chromatin and transcription

    Gbolahan Bamgbose, Guillaume Bordet ... Alexei Tulin
    Research Article

    PARP-1 is central to transcriptional regulation under both normal and stress conditions, with the governing mechanisms yet to be fully understood. Our biochemical and ChIP-seq-based analyses showed that PARP-1 binds specifically to active histone marks, particularly H4K20me1. We found that H4K20me1 plays a critical role in facilitating PARP-1 binding and the regulation of PARP-1-dependent loci during both development and heat shock stress. Here, we report that the sole H4K20 mono-methylase, pr-set7, and parp-1 Drosophila mutants undergo developmental arrest. RNA-seq analysis showed an absolute correlation between PR-SET7- and PARP-1-dependent loci expression, confirming co-regulation during developmental phases. PARP-1 and PR-SET7 are both essential for activating hsp70 and other heat shock genes during heat stress, with a notable increase of H4K20me1 at their gene body. Mutating pr-set7 disrupts monomethylation of H4K20 along heat shock loci and abolish PARP-1 binding there. These data strongly suggest that H4 monomethylation is a key triggering point in PARP-1 dependent processes in chromatin.

    1. Cancer Biology
    2. Genetics and Genomics

    Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD)

    Ting Zhang, Alisa Ambrodji ... Steven M Offer
    Research Article

    Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

    1. Computational and Systems Biology
    2. Genetics and Genomics

    Imputation of 3D genome structure by genetic–epigenetic interaction modeling in mice

    Lauren Kuffler, Daniel A Skelly ... Gregory W Carter
    Research Article

    Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic–epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic–epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.