Search

Summary Background Infection with Plasmodium falciparum leads to severe malaria and death in approximately 400 000 children each year in sub-Saharan Africa. Blood transfusion might benefit some patients with malaria but could potentially harm others. The aim of this study was to estimate the association between transfusion and death among children admitted to hospital with P falciparum malaria. Methods In this prospective, multicentre observational study, we analysed admissions to six tertiary care hospitals in The Gambia, Malawi, Gabon, Kenya, and Ghana that participated in the Severe Malaria in African Children network. Patients were enrolled if they were younger than 180 months and had a Giemsa-stained thick blood smear that was positive for P falciparum. Blood transfusion (whole blood at a target volume of 20 mL per kg) was administered at the discretion of the responsible physicians who were aware of local and international transfusion guidelines. The primary endpoint was death associated with transfusion, which was estimated using models adjusted for site and disease severity. We also aimed to identify factors associated with the decision to transfuse. The exploratory objective was to estimate optimal haemoglobin transfusion thresholds using generalised additive models. Findings Between Dec 19, 2000, and March 8, 2005, 26 106 patients were enrolled in the study, 25 893 of whom had their transfusion status recorded and were included in the primary analysis. 8513 (32·8%) patients received a blood transfusion. Patients were followed-up until discharge from hospital for a median of 2 days (IQR 1–4). 405 (4·8%) of 8513 patients who received a transfusion died compared with 689 (4·0%) of 17 380 patients who did not receive a transfusion. Transfusion was associated with decreased odds of death in site-adjusted analysis (odds ratio [OR] 0·82 [95% CI 0·71–0·94]) and after adjusting for the increased disease severity of patients who received a transfusion (0·50 [0·42–0·60]). Severe anaemia, elevated lactate concentration, respiratory distress, and parasite density were associated with greater odds of receiving a transfusion. Among all study participants, transfusion was associated with improved survival when the admission haemoglobin concentration was up to 77 g/L (95% CI 65–110). Among those with impaired consciousness (Blantyre Coma Score ≤4), transfusion was associated with improved survival at haemoglobin concentrations up to 105 g/L (95% CI 71–115). Among those with hyperlactataemia (blood lactate ≥5·0 mmol/L), transfusion was not significantly associated with harm at any haemoglobin concentration—ie, the OR of death comparing transfused versus not transfused was less than 1 at all haemoglobin concentrations (lower bound of the 95% CI for the haemoglobin concentration at which the OR of death equals 1: 90 g/L; no upper bound). Interpretation Our findings suggest that whole blood transfusion was associated with improved survival among children hospitalised with P falciparum malaria. Among those with impaired consciousness or hyperlactataemia, transfusion was associated with improved survival at haemoglobin concentrations above the currently recommended transfusion threshold. These findings highlight the need to do randomised controlled trials to test higher transfusion thresholds among African children with severe malaria complicated by these factors. Funding US National Institute of Allergy and Infectious Diseases.

More than 400,000 deaths from severe malaria (SM) are reported every year, mainly in African children. The diversity of clinical presentations associated with SM indicates important differences in disease pathogenesis that require specific treatment, and this clinical heterogeneity of SM remains poorly understood. Here, we apply tools from machine learning and model-based inference to harness large-scale data and dissect the heterogeneity in patterns of clinical features associated with SM in 2904 Gambian children admitted to hospital with malaria. This quantitative analysis reveals features predicting the severity of individual patient outcomes, and the dynamic pathways of SM progression, notably inferred without requiring longitudinal observations. Bayesian inference of these pathways allows us assign quantitative mortality risks to individual patients. By independently surveying expert practitioners, we show that this data-driven approach agrees with and expands the current state of knowledge on malaria progression, while simultaneously providing a data-supported framework for predicting clinical risk.

Host genetic factors can confer resistance against malaria, raising the question of whether this has led to evolutionary adaptation of parasite populations. In this study we investigated the correlation between host and parasite genetic variation in 4,171 Gambian and Kenya children ascertained with severe malaria due to Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and variation in three regions of the parasite genome, including nonsynonymous variants in the acyl-CoA synthetase family member PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The HbS-associated parasite alleles are in strong linkage disequilibrium and have frequencies which covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome, and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

The parasite Plasmodium falciparum is the main cause of severe malaria (SM). Despite treatment with antimalarial drugs, more than 450,000 SM deaths are reported every year, mainly in African children. The diversity of clinical presentations associated with SM indicate important differences in disease pathogenesis that often require specific treatment, and this clinical heterogeneity of SM is largely unresolved. In this study, we apply new machine learning and inference tools for large-scale data analysis to dissect the heterogeneity in patterns of clinical features associated with SM in 2,695 Gambian children admitted to hospital with Plasmodium falciparum malaria. This quantitative analysis, including the powerful HyperTraPS algorithm for inference of progressive processes, reveals pathways of SM symptom progression and features predicting the severity of individual patient outcomes. Notably, our approach allows the identification and dissection of disease progression pathways without the need for longitudinal observations. Learning these pathways and features from this rich dataset allows us to construct several quantitative measures of the mortality risk associated with a patient presenting with a given set of symptoms. By independently surveying expert practitioners, we show that this data-driven approach agrees with and expands the current state of knowledge on malaria progression, while simultaneously providing a data-supported framework for predicting clinical risk.

Additional file 1. List of differentially expressed proteins in fatal CM cases as measured by LC-MS/MS analysis.

INTRODUCTION Malaria parasites cause human disease by invading and replicating inside red blood cells. In the case of Plasmodium falciparum , this can lead to severe forms of malaria that are a major cause of childhood mortality in Africa. This species of parasite enters the red blood cell through interactions with surface proteins including the glycophorins GYPA and GYPB, which determine the polymorphic MNS blood group system. In a recent genome-wide association study, we identified alleles associated with protection against severe malaria near the cluster of genes encoding these invasion receptors. RATIONALE Investigation of genetic variants at this locus and their relation to severe malaria is challenging because of the high sequence similarity between the neighboring glycophorin genes and the relative lack of available sequence data capturing the genetic diversity of sub-Saharan Africa. To better assess whether variation in the glycophorin genes could explain the signal of association, we generated additional sequence data from sub-Saharan African populations and developed an analytical approach to characterize structural variation at this complex locus. RESULTS Using 765 newly sequenced human genomes from 10 African ethnic groups along with data from the 1000 Genomes Project, we generated a reference panel of haplotypes across the glycophorin region. In addition to single-nucleotide polymorphisms and short indels, we assayed large copy number variants (CNVs) using sequencing read depth and uncovered extensive structural diversity. By imputing from this reference panel into 4579 severe malaria cases and 5310 controls from three African populations, we found that a complex CNV, here called DUP4, is associated with resistance to severe malaria and fully explains the previously reported signal of association. In our sample, DUP4 is present only in east Africa, and this localization, as well as the extent of similarity between DUP4 haplotypes, suggests that it has recently increased in frequency, presumably under natural selection due to malaria. To evaluate the potential functional consequences of this structural variant, we analyzed high-coverage sequence-read data from multiple individuals to generate a model of the DUP4 chromosome structure. The DUP4 haplotype contains five glycophorin genes, including two hybrid genes that juxtapose the extracellular domain of GYPB with the transmembrane and intracellular domains of GYPA. Noting that these predicted hybrids are characteristic of the Dantu antigen in the MNS blood group system, we sequenced a Dantu positive individual and confirmed that DUP4 is the molecular basis of the Dantu NE blood group variant. CONCLUSION Although a role for GYPA and GYPB in parasite invasion is well known, a direct link between glycophorin polymorphisms and clinical susceptibility to malaria has been elusive. Here we have provided a systematic catalog of CNVs, describing structural diversity that may have functional importance at this locus. Our results identify a specific variant that encodes hybrid glycophorin proteins and is associated with protection against severe malaria. This discovery calls for further work to determine how this particular molecular rearrangement affects parasite invasion and the red blood cell response and may lead us toward new parasite vulnerabilities that can be utilized in future interventions against this deadly disease.

Interferon (IFN)-gamma is a critical mediator of immunity to malaria. This study explored the relationship between polymorphisms in the promoter region of the gene encoding IFN-gamma receptor 1 (IFNGR1) and susceptibility to malaria in African children. Four polymorphisms were found in the region between -1400 and +100 nt of the translational start site by sequencing, and analysis of 562 nuclear families revealed 6 haplotypes. Case-control analysis of 562 Gambian children with severe malaria and 569 umbilical cord blood samples (controls) showed that in Mandinka, the major Gambian ethnic group, heterozygotes for the IFNGR1-56 polymorphism were protected against cerebral malaria (odds ratio, 0.54; P=.016) and against death resulting from cerebral malaria (odds ratio, 0.22; P=.006). Analysis of a family study by transmission disequilibrium testing revealed a similar result. Further data are needed to validate this finding, but these results are reminiscent of those for other well-established heterozygote advantages, such as that associated with hemoglobin S.

Interferon Regulatory Factor 1 (IRF-1) is a member of the IRF family of transcription factors, which have key and diverse roles in the gene-regulatory networks of the immune system. IRF-1 has been described as a critical mediator of IFN-gamma signalling and as the major player in driving TH1 type responses. It is therefore likely to be crucial in both innate and adaptive responses against intracellular pathogens such as Plasmodium falciparum. Polymorphisms at the human IRF1 locus have been previously found to be associated with the ability to control P. falciparum infection in populations naturally exposed to malaria. In order to test whether genetic variation at the IRF1 locus also affects the risk of developing severe malaria, we performed a family-based test of association for 18 Single Nucleotide Polymorphisms (SNPs) across the gene in three African populations, using genotype data from 961 trios consisting of one affected child and his/her two parents (555 from The Gambia, 204 from Kenya and 202 from Malawi). No significant association with severe malaria or severe malaria subphenotypes (cerebral malaria and severe malaria anaemia) was observed for any of the SNPs/haplotypes tested in any of the study populations. Our results offer no evidence that the molecular pathways regulated by the transcription factor IRF-1 are involved in the immune-based pathogenesis of severe malaria.

Functional studies have demonstrated an interaction between the stimulatory G protein alpha subunit (G-alpha-s) and the malaria parasite at a cellular level. Obstruction of signal transduction via the erythrocyte G-alpha-s subunit reduced invasion by Plasmodium falciparum parasites. We sought to determine whether this signal pathway had an impact at the disease level by testing polymorphisms in the gene encoding G-alpha-s (GNAS) for association with severe malaria in a large multi-centre study encompassing family and case-control studies from The Gambia, Kenya and Malawi, and a case-control study from Ghana. We gained power to detect association using meta-analysis across the seven studies, with an overall sample size approximating 4,000 cases and 4,000 controls. Out of 12 SNPs investigated in the 19 kb GNAS region, four presented signals of association (P < 0.05) with severe malaria. The strongest single-locus association demonstrated an odds ratio of 1.13 (1.05-1.21), P = 0.001. Three of the loci presenting significant associations were clustered at the 5-prime end of the GNAS gene. Accordingly, haplotypes constructed from these loci demonstrated significant associations with severe malaria [OR = 0.88 (0.81-0.96), P = 0.005 and OR = 1.12 (1.03-1.20), P = 0.005]. The evidence presented here indicates that the influence of G-alpha-s on erythrocyte invasion efficacy may, indeed, alter individual susceptibility to disease.

Several lines of evidence link glucose-6-phosphate dehydrogenase (G6PD) deficiency to protection from severe malaria. Early reports suggested most G6PD deficiency in sub-Saharan Africa was because of the 202A/376G G6PD A− allele, and recent association studies of G6PD deficiency have employed genotyping as a convenient way to determine enzyme status. However, further work has suggested that other G6PD deficiency alleles are relatively common in some regions of West Africa. To investigate the consequences of unrecognized allelic heterogeneity on association studies, in particular studies of G6PD deficiency and malaria, we carried out a case–control analysis of 2488 Gambian children with severe malaria and 3875 controls. No significant association was found between severe malaria and the 202A/376G G6PD A− allele when analyzed alone, but pooling 202A/376G with other deficiency alleles revealed the signal of protection (male odds ratio (OR) 0.77, 95% CI 0.62–0.95, P=0.016; female OR 0.71, 95% CI 0.56–0.89, P=0.004). We have identified the 968C mutation as the most common G6PD A− allele in The Gambia. Our results highlight some of the consequences of allelic heterogeneity, particularly the increased type I error. They also suggest that G6PD-deficient male hemizygotes and female heterozygotes are protected from severe malaria.

Similarity between two individuals in the combination of genetic markers along their chromosomes indicates shared ancestry and can be used to identify historical connections between different population groups due to admixture. We use a genome-wide, haplotype-based, analysis to characterise the structure of genetic diversity and gene-flow in a collection of 48 sub-Saharan African groups. We show that coastal populations experienced an influx of Eurasian haplotypes over the last 7000 years, and that Eastern and Southern Niger-Congo speaking groups share ancestry with Central West Africans as a result of recent population expansions. In fact, most sub-Saharan populations share ancestry with groups from outside of their current geographic region as a result of gene-flow within the last 4000 years. Our in-depth analysis provides insight into haplotype sharing across different ethno-linguistic groups and the recent movement of alleles into new environments, both of which are relevant to studies of genetic epidemiology. DOI: http://dx.doi.org/10.7554/eLife.15266.001, eLife digest Our genomes contain a record of historical events. This is because when groups of people are separated for generations, the DNA sequence in the two groups’ genomes will change in different ways. Looking at the differences in the genomes of people from the same population can help researchers to understand and reconstruct the historical interactions that brought their ancestors together. The mixing of two populations that were previously separate is known as admixture. Africa as a continent has few written records of its history. This means that it is somewhat unknown which important movements of people in the past generated the populations found in modern-day Africa. Busby et al. have now attempted to use DNA to look into this and reconstruct the last 4000 years of genetic history in African populations. As has been shown in other regions of the world, the new analysis showed that all African populations are the result of historical admixture events. However, Busby et al. could characterize these events to unprecedented level of detail. For example, multiple ethnic groups from The Gambia and Mali all show signs of sharing the same set of ancestors from West Africa, Europe and Asia who mixed around 2000 years ago. Evidence of a migration of people from Central West Africa, known as the Bantu expansion, could also be detected, and was shown to carry genes to the south and east. An important next step will be to now look at the consequences of the observed gene-flow, and ask if it has contributed to spreading beneficial, or detrimental, mutations around Africa. DOI: http://dx.doi.org/10.7554/eLife.15266.002

Background Glucose-6-phosphate dehydrogenase (G6PD) deficiency exhibits considerable allelic heterogeneity which manifests with variable biochemical and clinical penetrance. It has long been thought that G6PD deficiency confers partial protection against severe malaria, however prior genetic association studies have disagreed with regard to the strength and specificity of a protective effect, which might reflect differences in the host genetic background, environmental influences, or in the specific clinical phenotypes considered. Methods A case-control association study of severe malaria was conducted in The Gambia, a region in West Africa where there is considerable allelic heterogeneity underlying expression of G6PD deficiency trait, evaluating the three major nonsynonymous polymorphisms known to be associated with enzyme deficiency (A968G, T542A, and C202T) in a cohort of 3836 controls and 2379 severe malaria cases. Results Each deficiency allele exhibited a similar trend toward protection against severe malaria overall (15–26 % reduced risk); however, in stratifying severe malaria to two of its constituent clinical subphenotypes, severe malarial anaemia (SMA) and cerebral malaria (CM), the three deficiency alleles exhibited trends of opposing effect, with risk conferred to SMA and protection with respect to CM. To assess the overall effect of G6PD deficiency trait, deficiency alleles found across all three loci were pooled. G6PD deficiency trait was found to be significantly associated with protection from severe malaria overall (OR 0.83 [0.75–0.92], \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P = 0.0006$$\end{document}P=0.0006), but this was limited to CM (OR 0.73 [0.61–0.87], \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P = 0.0005$$\end{document}P=0.0005), with a trend toward increased risk for SMA, especially in fully-deficient individuals (OR 1.43 [0.99–2.08], \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P = 0.056$$\end{document}P=0.056). Sex-stratified testing largely comported with these results, with evidence suggesting that protection by G6PD deficiency trait is conferred to both males and females, though susceptibility to SMA may be restricted to fully-deficient male hemizygotes. Conclusions In a part of Africa where multiple alleles contribute to expression of G6PD deficiency trait, these findings clarify and extend previous work done in populations where a single variant predominates, and taken together suggest a causal role for G6PD deficiency trait itself with respect to severe malaria, with opposing effects seen on two major clinical subphenotypes.

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5 million genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N ≤ 71,225 European ancestry, N ≤ 12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N = 29,136). We identified association between systolic or diastolic blood pressure and common variants in eight regions near the CYP17A1 (P = 7 × 10(-24)), CYP1A2 (P = 1 × 10(-23)), FGF5 (P = 1 × 10(-21)), SH2B3 (P = 3 × 10(-18)), MTHFR (P = 2 × 10(-13)), c10orf107 (P = 1 × 10(-9)), ZNF652 (P = 5 × 10(-9)) and PLCD3 (P = 1 × 10(-8)) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.

We report a genome-wide association (GWA) study of severe malaria in The Gambia. The initial GWA scan included 2,500 children genotyped on the Affymetrix 500K GeneChip, and a replication study included 3,400 children. We used this to examine the performance of GWA methods in Africa. We found considerable population stratification, and also that signals of association at known malaria resistance loci were greatly attenuated owing to weak linkage disequilibrium (LD). To investigate possible solutions to the problem of low LD, we focused on the HbS locus, sequencing this region of the genome in 62 Gambian individuals and then using these data to conduct multipoint imputation in the GWA samples. This increased the signal of association, from P = 4 × 10(-7) to P = 4 × 10(-14), with the peak of the signal located precisely at the HbS causal variant. Our findings provide proof of principle that fine-resolution multipoint imputation, based on population-specific sequencing data, can substantially boost authentic GWA signals and enable fine mapping of causal variants in African populations.

The tumor necrosis factor gene (TNF) and lymphotoxin-alpha gene (LTA) have long attracted attention as candidate genes for susceptibility traits for malaria, and several of their polymorphisms have been found to be associated with severe malaria (SM) phenotypes. In a large study involving >10,000 individuals and encompassing 3 African populations, we found evidence to support the reported associations between the TNF -238 polymorphism and SM in The Gambia. However, no TNF/LTA polymorphisms were found to be associated with SM in cohorts in Kenya and Malawi. It has been suggested that the causal polymorphisms regulating the TNF and LTA responses may be located some distance from the genes. Therefore, more-detailed mapping of variants across TNF/LTA genes and their flanking regions in the Gambian and allied populations may need to be undertaken to find any causal polymorphisms.

There is growing epidemiological and molecular evidence that ABO blood group affects host susceptibility to severe Plasmodium falciparum infection. The high frequency of common ABO alleles means that even modest differences in susceptibility could have a significant impact on the health of people living in malaria endemic regions. We performed an association study, the first to utilize key molecular genetic variation underlying the ABO system, genotyping >9000 individuals across three African populations. Using population- and family-based tests, we demonstrated that alleles producing functional ABO enzymes are associated with greater risk of severe malaria phenotypes (particularly malarial anemia) in comparison with the frameshift deletion underlying blood group O: case-control allelic odds ratio (OR), 1.2; 95% confidence interval (CI), 1.09-1.32; P = 0.0003; family-studies allelic OR, 1.19; 95% CI, 1.08-1.32; P = 0.001; pooled across all studies allelic OR, 1.18; 95% CI, 1.11-1.26; P = 2 x 10(-7). We found suggestive evidence of a parent-of-origin effect at the ABO locus by analyzing the family trios. Non-O haplotypes inherited from mothers, but not fathers, are significantly associated with severe malaria (likelihood ratio test of Weinberg, P = 0.046). Finally, we used HapMap data to demonstrate a region of low F(ST) (-0.001) between the three main HapMap population groups across the ABO locus, an outlier in the empirical distribution of F(ST) across chromosome 9 (approximately 99.5-99.9th centile). This low F(ST) region may be a signal of long-standing balancing selection at the ABO locus, caused by multiple infectious pathogens including P. falciparum.

The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations. Here we report completion of the project, having reconstructed the genomes of 2,504 individuals from 26 populations using a combination of low-coverage whole-genome sequencing, deep exome sequencing, and dense microarray genotyping. We characterized a broad spectrum of genetic variation, in total over 88 million variants (84.7 million single nucleotide polymorphisms (SNPs), 3.6 million short insertions/deletions (indels), and 60,000 structural variants), all phased onto high-quality haplotypes. This resource includes >99% of SNP variants with a frequency of >1% for a variety of ancestries. We describe the distribution of genetic variation across the global sample, and discuss the implications for common disease studies., Wellcome Trust (London, England) (Core Award 090532/Z/09/Z), Wellcome Trust (London, England) (Senior Investigator Award 095552/Z/11/Z ), Wellcome Trust (London, England) (WT095908), Wellcome Trust (London, England) (WT109497), Wellcome Trust (London, England) (WT098051), Wellcome Trust (London, England) (WT086084/Z/08/Z), Wellcome Trust (London, England) (WT100956/Z/13/Z ), Wellcome Trust (London, England) (WT097307), Wellcome Trust (London, England) (WT0855322/Z/08/Z ), Wellcome Trust (London, England) (WT090770/Z/09/Z ), Wellcome Trust (London, England) (Major Overseas program in Vietnam grant 089276/Z.09/Z), Medical Research Council (Great Britain) (grant G0801823), Biotechnology and Biological Sciences Research Council (Great Britain) (grant BB/I02593X/1), Biotechnology and Biological Sciences Research Council (Great Britain) (grant BB/I021213/1), Zhongguo ke xue ji shu qing bao yan jiu suo. Office of 863 Programme of China (2012AA02A201), National Basic Research Program of China (2011CB809201), National Basic Research Program of China (2011CB809202), National Basic Research Program of China (2011CB809203), National Natural Science Foundation of China (31161130357), Shenzhen Municipal Government of China (grant ZYC201105170397A), Canadian Institutes of Health Research (grant 136855), Quebec Ministry of Economic Development, Innovation, and Exports (PSR-SIIRI-195), Germany. Bundesministerium für Bildung und Forschung (0315428A), Germany. Bundesministerium für Bildung und Forschung (01GS08201), Germany. Bundesministerium für Bildung und Forschung (BMBF-EPITREAT grant 0316190A), Deutsche Forschungsgemeinschaft (Emmy Noether Grant KO4037/1-1), Beatriu de Pinos Program (2006 BP-A 10144), Beatriu de Pinos Program (2009 BP-B 00274), Spanish National Institute for Health (grant PRB2 IPT13/0001-ISCIII-SGEFI/FEDER), Japan Society for the Promotion of Science (fellowship number PE13075), Marie Curie Actions Career Integration (grant 303772), Fonds National Suisse del la Recherche, SNSF, Scientifique (31003A_130342), National Center for Biotechnology Information (U.S.) (U54HG3067), National Center for Biotechnology Information (U.S.) (U54HG3273), National Center for Biotechnology Information (U.S.) (U01HG5211), National Center for Biotechnology Information (U.S.) (U54HG3079), National Center for Biotechnology Information (U.S.) (R01HG2898), National Center for Biotechnology Information (U.S.) (R01HG2385), National Center for Biotechnology Information (U.S.) (RC2HG5552), National Center for Biotechnology Information (U.S.) (U01HG6513), National Center for Biotechnology Information (U.S.) (U01HG5214), National Center for Biotechnology Information (U.S.) (U01HG5715), National Center for Biotechnology Information (U.S.) (U01HG5718), National Center for Biotechnology Information (U.S.) (U01HG5728), National Center for Biotechnology Information (U.S.) (U41HG7635), National Center for Biotechnology Information (U.S.) (U41HG7497), National Center for Biotechnology Information (U.S.) (R01HG4960), National Center for Biotechnology Information (U.S.) (R01HG5701), National Center for Biotechnology Information (U.S.) (R01HG5214), National Center for Biotechnology Information (U.S.) (R01HG6855), National Center for Biotechnology Information (U.S.) (R01HG7068), National Center for Biotechnology Information (U.S.) (R01HG7644), National Center for Biotechnology Information (U.S.) (DP2OD6514), National Center for Biotechnology Information (U.S.) (DP5OD9154), National Center for Biotechnology Information (U.S.) (R01CA166661), National Center for Biotechnology Information (U.S.) (R01CA172652), National Center for Biotechnology Information (U.S.) (P01GM99568), National Center for Biotechnology Information (U.S.) (R01GM59290), National Center for Biotechnology Information (U.S.) (R01GM104390), National Center for Biotechnology Information (U.S.) (T32GM7790), National Center for Biotechnology Information (U.S.) (R01HL87699), National Center for Biotechnology Information (U.S.) (R01HL104608), National Center for Biotechnology Information (U.S.) (T32HL94284), National Center for Biotechnology Information (U.S.) (HHSN268201100040C), National Center for Biotechnology Information (U.S.) (HHSN272201000025C), Lundbeck Foundation (grant R170-2014-1039, Simons Foundation (SFARI award SF51), National Science Foundation (U.S.) (Research Fellowship DGE-1147470)

There is growing interest in the use of haplotype-based methods for detecting recent selection. Here, we describe a method that uses a sliding window to estimate similarity among the haplotypes associated with any given single-nucleotide polymorphism (SNP) allele. We used simulations of natural selection to provide estimates of the empirical power of the method to detect recently selected alleles and found it to be comparable in power to the popular long-range haplotype test and more powerful than methods based on nucleotide diversity. We then applied the method to a recently selected allele--the sickle mutation at the HBB locus--and found it to have a signal of selection that was significantly stronger than that of simulated models both with and without strong selection. Using this method, we also evaluated >4,000 SNPs on chromosome 20, indicating the applicability of the method to regional data sets.

Given the importance of Africa to studies of human origins and disease susceptibility, detailed characterization of African genetic diversity is needed. The African Genome Variation Project provides a resource with which to design, implement and interpret genomic studies in sub-Saharan Africa and worldwide. The African Genome Variation Project represents dense genotypes from 1,481 individuals and whole-genome sequences from 320 individuals across sub-Saharan Africa. Using this resource, we find novel evidence of complex, regionally distinct hunter-gatherer and Eurasian admixture across sub-Saharan Africa. We identify new loci under selection, including loci related to malaria susceptibility and hypertension. We show that modern imputation panels (sets of reference genotypes from which unobserved or missing genotypes in study sets can be inferred) can identify association signals at highly differentiated loci across populations in sub-Saharan Africa. Using whole-genome sequencing, we demonstrate further improvements in imputation accuracy, strengthening the case for large-scale sequencing efforts of diverse African haplotypes. Finally, we present an efficient genotype array design capturing common genetic variation in Africa.

Lipocalin 2 distinguishes severe and bacterial pneumonia from nonsevere and nonbacterial pneumonia with a high level of precision. The clinical impact of this biomarker requires large-scale clinical evaluation., Background. Pneumonia is the leading cause of death in children globally. Clinical algorithms remain suboptimal for distinguishing severe pneumonia from other causes of respiratory distress such as malaria or distinguishing bacterial pneumonia and pneumonia from others causes, such as viruses. Molecular tools could improve diagnosis and management. Methods. We conducted a mass spectrometry–based proteomic study to identify and validate markers of severity in 390 Gambian children with pneumonia (n = 204) and age-, sex-, and neighborhood-matched controls (n = 186). Independent validation was conducted in 293 Kenyan children with respiratory distress (238 with pneumonia, 41 with Plasmodium falciparum malaria, and 14 with both). Predictive value was estimated by the area under the receiver operating characteristic curve (AUC). Results. Lipocalin 2 (Lpc-2) was the best protein biomarker of severe pneumonia (AUC, 0.71 [95% confidence interval, .64–.79]) and highly predictive of bacteremia (78% [64%–92%]), pneumococcal bacteremia (84% [71%–98%]), and “probable bacterial etiology” (91% [84%–98%]). These results were validated in Kenyan children with severe malaria and respiratory distress who also met the World Health Organization definition of pneumonia. The combination of Lpc-2 and haptoglobin distinguished bacterial versus malaria origin of respiratory distress with high sensitivity and specificity in Gambian children (AUC, 99% [95% confidence interval, 99%–100%]) and Kenyan children (82% [74%–91%]). Conclusions. Lpc-2 and haptoglobin can help discriminate the etiology of clinically defined pneumonia and could be used to improve clinical management. These biomarkers should be further evaluated in prospective clinical studies.

The last author's name was spelled incorrectly. The correct name is: Malaria Genomic Epidemiology Network

Acquired immunity in vertebrates maintains polymorphisms in endemic pathogens, leading to identifiable signatures of balancing selection. To comprehensively survey for genes under such selection in the human malaria parasite Plasmodium falciparum, we generated paired-end short-read sequences of parasites in clinical isolates from an endemic Gambian population, which were mapped to the 3D7 strain reference genome to yield high-quality genome-wide coding sequence data for 65 isolates. A minority of genes did not map reliably, including the hypervariable var, rifin, and stevor families, but 5,056 genes (90.9% of all in the genome) had >70% sequence coverage with minimum read depth of 5 for at least 50 isolates, of which 2,853 genes contained 3 or more single nucleotide polymorphisms (SNPs) for analysis of polymorphic site frequency spectra. Against an overall background of negatively skewed frequencies, as expected from historical population expansion combined with purifying selection, the outlying minority of genes with signatures indicating exceptionally intermediate frequencies were identified. Comparing genes with different stage-specificity, such signatures were most common in those with peak expression at the merozoite stage that invades erythrocytes. Members of clag, PfMC-2TM, surfin, and msp3-like gene families were highly represented, the strongest signature being in the msp3-like gene PF10_0355. Analysis of msp3-like transcripts in 45 clinical and 11 laboratory adapted isolates grown to merozoite-containing schizont stages revealed surprisingly low expression of PF10_0355. In diverse clonal parasite lines the protein product was expressed in a minority of mature schizonts (, Author Summary The memory component of acquired immune responses selects for distinctive patterns of polymorphism in genes encoding important target antigens of pathogens. These are detectable by surveying for evidence of balancing selection, as previously illustrated in analyses of genes encoding malaria parasite antigens that are candidate targets of naturally acquired immunity. For a comprehensive screen to discover targets of immunity in the major human malaria parasite Plasmodium falciparum, an endemic population in West Africa was sampled and genome sequence data obtained from 65 clinical isolates, allowing analysis of polymorphism in almost all protein-coding genes. Antigen genes previously studied by capillary re-sequencing in independent population samples had highly concordant indices in the genome-wide analysis here, and this has identified other genes with stronger evidence of balancing selection, now prioritized for functional study and potential vaccine candidacy. The statistical signatures consistent with such selection were particularly common in genes with peak expression at the stage that invades erythrocytes, and members of several gene families were represented. The strongest signature was in the msp3-like gene PF10_0355, so we studied the transcript and protein product in parasites, revealing an unexpected pattern of phase variable expression. Variation in expression of polymorphic antigens under balancing selection may be more common than previously thought, requiring further study to assess vaccine candidacy.

The relevance of innate immune responses to Plasmodium falciparum infection, in particular the central role of natural killer (NK) cell-derived interferon gamma (IFN-γ), is becoming increasingly recognised. Recently, it has been shown that IFN-γ production in response to P. falciparum antigens is in part regulated by killer-cell immunoglobulin-like receptor (KIR) genes, and a study from malaria-exposed Melanesians suggested an association between KIR genotypes and susceptibility to infection. This prompted us to determine and compare the frequencies of 15 KIR genes in Gambian children presenting with either severe malaria (n = 133) or uncomplicated malaria (n = 188) and in cord-blood population control samples (n = 314) collected from the same area. While no significant differences were observed between severe and uncomplicated cases, proportions of individuals with KIR2DS2+C1 and KIR2DL2+C1 were significantly higher among malaria cases overall than in population control samples. In an exploratory analysis, activating KIR genes KIR2DS2, KIR3DS1 and KIR2DS5 were slightly higher in children in disease subgroups associated with the highest mortality. In addition, our data suggest that homozygosity for KIR genotype A might be associated with different malaria outcomes including protection from infection and higher blood parasitaemia levels in those that do get infected. These findings are consistent with a probable role of KIR genes in determining susceptibility to malaria, and further studies are warranted in different populations.

With the functional demonstration of a role in erythrocyte invasion by Plasmodium falciparum parasites, implications in the aetiology of common conditions that prevail in individuals of African origin, and a wealth of pharmacological knowledge, the stimulatory G protein (Gs) signal transduction pathway presents an exciting target for anti-malarial drug intervention. Having previously demonstrated a role for the G-alpha-s gene, GNAS, in severe malaria disease, we sought to identify other important components of the Gs pathway. Using meta-analysis across case-control and family trio (affected child and parental controls) studies of severe malaria from The Gambia and Malawi, we sought evidence of association in six Gs pathway candidate genes: adenosine receptor 2A (ADORA2A) and 2B (ADORA2B), beta-adrenergic receptor kinase 1 (ADRBK1), adenylyl cyclase 9 (ADCY9), G protein beta subunit 3 (GNB3), and regulator of G protein signalling 2 (RGS2). Our study amassed a total of 2278 cases and 2364 controls. Allele-based models of association were investigated in all genes, and genotype and haplotype-based models were investigated where significant allelic associations were identified. Although no significant associations were observed in the other genes, several were identified in ADORA2A. The most significant association was observed at the rs9624472 locus, where the G allele (approximately 20% frequency) appeared to confer enhanced risk to severe malaria [OR = 1.22 (1.09-1.37); P = 0.001]. Further investigation of the ADORA2A gene region is required to validate the associations identified here, and to identify and functionally characterize the responsible causal variant(s). Our results provide further evidence supporting a role of the Gs signal transduction pathway in the regulation of severe malaria, and request further exploration of this pathway in future studies.

The prevalence of CD36 deficiency in East Asian and African populations suggests that the causal variants are under selection by severe malaria. Previous analysis of data from the International HapMap Project indicated that a CD36 haplotype bearing a nonsense mutation (T1264G; rs3211938) had undergone recent positive selection in the Yoruba of Nigeria. To investigate the global distribution of this putative selection event, we genotyped T1264G in 3420 individuals from 66 populations. We confirmed the high frequency of 1264G in the Yoruba (26%). However, the 1264G allele is less common in other African populations and absent from all non-African populations without recent African admixture. Using long-range linkage disequilibrium, we studied two West African groups in depth. Evidence for recent positive selection at the locus was demonstrable in the Yoruba, although not in Gambians. We screened 70 variants from across CD36 for an association with severe malaria phenotypes, employing a case-control study of 1350 subjects and a family study of 1288 parent-offspring trios. No marker was significantly associated with severe malaria. We focused on T1264G, genotyping 10,922 samples from four African populations. The nonsense allele was not associated with severe malaria (pooled allelic odds ratio 1.0; 95% confidence interval 0.89-1.12; P = 0.98). These results suggest a range of possible explanations including the existence of alternative selection pressures on CD36, co-evolution between host and parasite or confounding caused by allelic heterogeneity of CD36 deficiency.

The tumour necrosis factor (TNF) gene and other genes flanking it in the major histocompatibility complex (MHC) class III region are potentially important mediators of both immunity and pathogenesis of malaria. We investigated the association of severe malaria with eleven haplotype tagging-polymorphisms for eleven MHC class III candidate genes, including TNF, lymphotoxin alpha (LTA), allograft inflammatory factor 1 (AIF1), and HLA-B associated transcript 2 (BAT2). An analysis of 2162 case-controls demonstrated the first evidence of association between a BAT2 polymorphism (rs1046089) and severe malaria.

Large-scale studies of genomic variation could assist efforts to eliminate malaria. But there are scientific, ethical and practical challenges to carrying out such studies in developing countries, where the burden of disease is greatest. The Malaria Genomic Epidemiology Network (MalariaGEN) is now working to overcome these obstacles, using a consortial approach that brings together researchers from 21 countries. © 2008 Macmillan Publishers Limited. All rights reserved.

Evidence from autopsy, mouse-model and in vitro binding studies suggests that adhesion of erythrocytes infected with Plasmodium falciparum to the human host intercellular adhesion molecule (ICAM)-1 receptor is important in the pathogenesis of severe malaria. Previous association studies between polymorphisms around the ICAM1 gene and susceptibility to severe malarial phenotypes have been inconclusive and often contradictory. We performed genetic association studies with 15 single-nucleotide-polymorphisms (SNPs) around the ICAM1 locus. All SNPs were screened in a family study of 1071 trios from Gambia, Malawi and Kenya. Two key non-synonymous SNPs with previously reported associations, rs5491 (K56M or ‘ICAM-1Kilifi’) and rs5498 (K469E), were tested in an additional 708 Gambian trios and a case-control study of 4058 individuals. None of the polymorphisms were associated with severe malaria phenotypes. Pooled results across our studies for ICAM-1Kilifi were, in severe malaria, odds ratio (OR) 1.02, 95% confidence interval (CI) 0.96 – 1.09, P=0.54, and cerebral malaria OR 1.07, CI 0.97 – 1.17, P=0.17. We assess the available epidemiological, population genetic and functional evidence which links ICAM-1Kilifi to severe malaria susceptibility.

There is presently much interest in utilizing patterns of linkage disequilibrium (LD) to further genetic association studies. This is particularly pertinent in the class III region of the human major histocompatibility complex (MHC), which has been extensively studied as a disease susceptibility locus in a number of ethnic groups. To date, however, few studies of LD in the MHC have considered non-Caucasian populations. With the advent of large-scale haplotyping of the human genome, the question of utilizing LD patterns across populations has come to the fore. We have previously used LD mapping to direct an MHC class III association study in a UK Caucasian population. As an extension of this, we sought to determine to what extent the pattern of LD observed in that study could be used to conduct a similar study in a West African Gambian population. We found that broad patterns of LD were similar in the two populations, resulting in similar candidate region delineations, but at a higher resolution, marker-specific patterns of LD and population-dependent allele frequencies confounded the choice of regional tagging SNPs. Our results have implications for the applicability of large-scale haplotype maps such as the HapMap to complex regions like the MHC.

To assess the hypothesis that nitric oxide is critical in the pathogenesis of cerebral malaria, we analysed genetic variation in the proximal promoter region of NOS2A, the gene encoding inducible nitric oxide synthase. Sequencing 72 Gambian chromosomes revealed 11 single nucleotide polymorphisms in 2.5 kB (theta=8.6 x 10(-4)). Genotyping 104 nuclear families identified six common haplotypes. A single haplotype, uniquely defined by the NOS2A-1659T allele, was associated with cerebral malaria by a transmission disequilibrium test of 334 affected children and their parents (P=0.02). An independent case-control study of 505 different children from the same population replicated the allelic association with cerebral malaria (odds ratio: 1.31, P=0.04). Taken together these data indicate a weak but significant association of the NOS2A locus with susceptibility to cerebral malaria. Despite high linkage disequilibrium across the region studied, this association would not have been detected without the initial construction of a dense marker set for haplotype tagging.

CD40 ligand (CD40L), a glycoprotein involved in B cell proliferation, antigen presenting cell activation, and Ig class switching, is important in the immune response to infection. Rare coding mutations in CD40L can lead to life-threatening immunodeficiency but the potential for common variants to alter disease susceptibility remains to be explored. To identify polymorphisms in CD40L, we sequenced 2.3 kb of the 5' flanking region and the first exon of the gene in DNA samples from 36 Gambian females and one chimpanzee. Diversity was lower than the average reported for other areas of the X chromosome, and only two polymorphisms were identified. The polymorphisms were genotyped in DNA samples from 957 Gambian individuals, cases and controls from a study of severe malaria. A significant reduction in risk for severe malaria (OR = 0.52, P = 0.002) was associated with males hemizygous for the CD40L-726C. Analysis by transmission disequilibrium test of 371 cases, for whom DNA from both parents was also available, confirmed the result was not due to stratification (P = 0.04). A similar but non-significant trend was found in females. This preliminary association of a common variant in CD40L with a malaria resistance phenotype encourages further genetic characterization of the role of CD40L in infectious disease.

Background During malaria infection the Toll-like receptor 9 (TLR9) is activated through induction with plasmodium DNA or another malaria motif not yet identified. Although TLR9 activation by malaria parasites is well reported, the implication to the susceptibility to severe malaria is not clear. The aim of this study was to assess the contribution of genetic variation at TLR9 to severe malaria. Methods This study explores the contribution of TLR9 genetic variants to severe malaria using two approaches. First, an association study of four common single nucleotide polymorphisms was performed on both family- and population-based studies from Malawian and Gambian populations (n>6000 individual). Subsequently, it was assessed whether TLR9 expression is affected by cis-acting variants and if these variants could be mapped. For this work, an allele specific expression (ASE) assay on a panel of HapMap cell lines was carried out. Results No convincing association was found with polymorphisms in TLR9 for malaria severity, in either Gambian or Malawian populations, using both case-control and family based study designs. Using an allele specific expression assay it was observed that TLR9 expression is affected by cis-acting variants, these results were replicated in a second experiment using biological replicates. Conclusion By using the largest cohorts analysed to date, as well as a standardized phenotype definition and study design, no association of TLR9 genetic variants with severe malaria was found. This analysis considered all common variants in the region, but it is remains possible that there are rare variants with association signals. This report also shows that TLR9 expression is potentially modulated through cis-regulatory variants, which may lead to differential inflammatory responses to infection between individuals.

Additional file 1. List of differentially expressed proteins in fatal CM cases as measured by LC-MS/MS analysis.

Additional file 2. Circulating proteasomes and PSMB9 in patients with cerebral malaria. Dot plots show the concentration of total circulating 20S proteasomes (left panels) and PSMB9 (right panels) in children with CM based on outcome (top), presence of seizures during admission (middle) and the correlation of these biomarkers with the number of seizures witnessed during admission (bottom).

Additional file 2. Circulating proteasomes and PSMB9 in patients with cerebral malaria. Dot plots show the concentration of total circulating 20S proteasomes (left panels) and PSMB9 (right panels) in children with CM based on outcome (top), presence of seizures during admission (middle) and the correlation of these biomarkers with the number of seizures witnessed during admission (bottom).

Background The sickle (βs) mutation in the beta-globin gene (HBB) occurs on five "classical" βs haplotype backgrounds in ethnic groups of African ancestry. Strong selection in favour of the βs allele – a consequence of protection from severe malarial infection afforded by heterozygotes – has been associated with a high degree of extended haplotype similarity. The relationship between classical βs haplotypes and long-range haplotype similarity may have both anthropological and clinical implications, but to date has not been explored. Here we evaluate the haplotype similarity of classical βs haplotypes over 400 kb in population samples from Jamaica, The Gambia, and among the Yoruba of Nigeria (Hapmap YRI). Results The most common βs sub-haplotype among Jamaicans and the Yoruba was the Benin haplotype, while in The Gambia the Senegal haplotype was observed most commonly. Both subtypes exhibited a high degree of long-range haplotype similarity extending across approximately 400 kb in all three populations. This long-range similarity was significantly greater than that seen for other haplotypes sampled in these populations (P < 0.001), and was independent of marker choice and marker density. Among the Yoruba, Benin haplotypes were highly conserved, with very strong linkage disequilibrium (LD) extending a megabase across the βs mutation. Conclusion Two different classical βs haplotypes, sampled from different populations, exhibit comparable and extensive long-range haplotype similarity and strong LD. This LD extends across the adjacent recombination hotspot, and is discernable at distances in excess of 400 kb. Although the multi-centric geographic distribution of βs haplotypes indicates strong subdivision among early Holocene sub-Saharan populations, we find no evidence that selective pressures imposed by falciparum malaria varied in intensity or timing between these subpopulations. Our observations also suggest that cis-acting loci, which may influence outcomes in sickle cell disease, could lie considerable distances away from β-globin.