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 malaria parasite Plasmodium falciparum kills on the order of a million African children each year 1 , and this is a small fraction of the number of infected individuals in the population [1][2][3] . In communities where everyone is repeatedly infected with P. falciparum, host genetic factors account for ~25% of the risk of severe malaria, that is, life-threatening forms of the disease 3 . The strongest known determinant of risk, hemoglobin S (HbS), accounts for 2% of the total variation, implying that only a small fraction of genetic resistance factors have so far been discovered 3 . Identifying the genetic basis of protective immunity against severe malaria may provide important insights for vaccine development.Here we examine the possibility of approaching this problem by genome-wide association (GWA) analysis. There are many unsolved methodological questions about how to conduct an effective GWA study in Africa 4 . High levels of ethnic diversity may result in false-positive associations owing to population structure. Variations in haplotype structure between different ethnic groups may reduce power to detect GWA signals, particularly when data are amalgamated across multiple study sites. Low LD implies the need for denser genotyping arrays than are currently available: a crude estimate is that an African GWA study with 1.5 million SNPs would have approximately the same statistical power as a European study with Jallow et al.Page 2Nat Genet. Author manuscript; available in PMC 2010 September 21. NIH-PA Author ManuscriptNIH-PA Author Manuscript NIH-PA Author Manuscript 0.6 million SNPs5, but this is based on HapMap data from a single ethnic group and a larger number of SNPs may be needed to achieve adequate power across different ethnic groups.We carried out an initial GWA study in Gambian children that explores these methodological questions. Genotyping of ~500,000 SNPs was conducted on 1,060 cases of severe malaria and 1...

Interactions of killer cell immunoglobulin-like receptors (KIRs) with major histocompatibility complex (MHC) class I ligands diversify natural killer cell responses to infection. By analyzing sequence variation in diverse human populations, we show that the KIR3DL1/S1 locus encodes two lineages of polymorphic inhibitory KIR3DL1 allotypes that recognize Bw4 epitopes of protein">HLA-A and HLA-B and one lineage of conserved activating KIR3DS1 allotypes, also implicated in Bw4 recognition. Balancing selection has maintained these three lineages for over 3 million years. Variation was selected at D1 and D2 domain residues that contact HLA class I and at two sites on D0, the domain that enhances the binding of KIR3D to HLA class I. HLA-B variants that gained Bw4 through interallelic microconversion are also products of selection. A worldwide comparison uncovers unusual KIR3DL1/S1 evolution in modern sub-Saharan Africans. Balancing selection is weak and confined to D0, KIR3DS1 is rare and KIR3DL1 allotypes with similar binding sites predominate. Natural killer cells express the dominant KIR3DL1 at a high frequency and with high surface density, providing strong responses to cells perturbed in Bw4 expression.

Many human genetic associations with resistance to malaria have been reported but few have been reliably replicated. We collected data on 11,890 cases of severe malaria due to Plasmodium falciparum and 17,441 controls from 12 locations in Africa, Asia and Oceania. There was strong evidence of association with the HBB, ABO, ATP2B4, G6PD and CD40LG loci but previously reported associations at 22 other loci did not replicate in the multi-centre analysis. The large sample size made it possible to identify authentic genetic effects that are heterogeneous across populations or phenotypes, a striking example being the main African form of G6PD deficiency, which reduced the risk of cerebral malaria but increased the risk of severe malarial anaemia. The finding that G6PD deficiency has opposing effects on different fatal complications of P. falciparum infection indicates that the evolutionary origins of this common human genetic disorder are more complex than previously supposed.

Interactions between HLA class I molecules and killer-cell immunoglobulin-like receptors (KIR) control natural killer cell (NK) functions in immunity and reproduction. Encoded by genes on different chromosomes, these polymorphic ligands and receptors correlate highly with disease resistance and susceptibility. Although studied at low-resolution in many populations, high-resolution analysis of combinatorial diversity of HLA class I and KIR is limited to Asian and Amerindian populations with low genetic diversity. At the other end of the spectrum is the West African population investigated here: we studied 235 individuals, including 104 mother-child pairs, from the Ga-Adangbe of Ghana. This population has a rich diversity of 175 KIR variants forming 208 KIR haplotypes, and 81 HLA-A, -B and -C variants forming 190 HLA class I haplotypes. Each individual we studied has a unique compound genotype of HLA class I and KIR, forming 1–14 functional ligand-receptor interactions. Maintaining this exceptionally high polymorphism is balancing selection. The centromeric region of the KIR locus, encoding HLA-C receptors, is highly diverse whereas the telomeric region encoding Bw4-specific KIR3DL1, lacks diversity in Africans. Present in the Ga-Adangbe are high frequencies of Bw4-bearing HLA-B*53:01 and Bw4-lacking HLA-B*35:01, which otherwise are identical. Balancing selection at key residues maintains numerous HLA-B allotypes having and lacking Bw4, and also those of stronger and weaker interaction with LILRB1, a KIR-related receptor. Correspondingly, there is a balance at key residues of KIR3DL1 that modulate its level of cell-surface expression. Thus, capacity to interact with NK cells synergizes with peptide binding diversity to drive HLA-B allele frequency distribution. These features of KIR and HLA are consistent with ongoing co-evolution and selection imposed by a pathogen endemic to West Africa. Because of the prevalence of malaria in the Ga-Adangbe and previous associations of cerebral malaria with HLA-B*53:01 and KIR, Plasmodium falciparum is a candidate pathogen.