Inactive alleles of cytochrome P450 2C19 may be positively selected in human evolution Genome evolution and evolutionary systems biology

© 2014 Janha et al.; licensee BioMed Central Ltd.Background: Cytochrome P450 CYP2C19 metabolizes a wide range of pharmacologically active substances and a relatively small number of naturally occurring environmental toxins. Poor activity alleles of CYP2C19 are very frequent worldwide, particularly in Asia, raising the possibility that reduced metabolism could be advantageous in some circumstances. The evolutionary selective forces acting on this gene have not previously been investigated. We analyzed CYP2C19 genetic markers from 127 Gambians and on 120 chromosomes from Yoruba, Europeans and Asians (Japanese + Han Chinese) in the Hapmap database. Haplotype breakdown was explored using bifurcation plots and relative extended haplotype homozygosity (REHH). Allele frequency differentiation across populations was estimated using the fixation index (FST) and haplotype diversity with coalescent models. Results: Bifurcation plots suggested conservation of alleles conferring slow metabolism (CYP2C19∗2 and ∗3). REHH was high around CYP2C19∗2 in Yoruba (REHH 8.3, at 133.3 kb from the core) and to a lesser extent in Europeans (3.5, at 37.7 kb) and Asians (2.8, at -29.7 kb). FST at the CYP2C19 locus was low overall (0.098). CYP2C19∗3 was an FST outlier in Asians (0.293), CYP2C19 haplotype diversity ST is low at the CYP2C19 locus, suggesting balancing selection overall. The biological factors responsible for these selective pressures are currently unknown. One possible explanation is that early humans were exposed to a ubiquitous novel toxin activated by CYP2C19. The genetic adaptation took place within the last 10,000 years which coincides with the development of systematic agricultural practices.This work was supported by the Medical Research Council Unit The Gambia and the European and Developing Countries Clinical Trials Partnership [grant number CG_ta_05_40204_018]

Malaria protection due to sickle haemoglobin depends on parasite genotype

Host genetic factors can confer resistance against malaria1, raising the question of whether this has led to evolutionary adaptation of parasite populations. Here we searched for association between candidate host and parasite genetic variants in 3,346 Gambian and Kenyan children with severe malaria caused by Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and three regions of the parasite genome, which is not explained by population structure or other covariates, and which is replicated in additional samples. The HbS-associated alleles include nonsynonymous variants in the gene for the acyl-CoA synthetase family member2-4 PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The alleles are in strong linkage disequilibrium and have frequencies that 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

Characterisation of the opposing effects of G6PD deficiency on cerebral malaria and severe malarial anaemia.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is believed to confer protection against Plasmodium falciparum malaria, but the precise nature of the protective effecthas proved difficult to define as G6PD deficiency has multiple allelic variants with different effects in males and females, and it has heterogeneous effects on the clinical outcome of P. falciparum infection. Here we report an analysis of multiple allelic forms of G6PD deficiency in a large multi-centre case-control study of severe malaria, using the WHO classification of G6PD mutations to estimate each individual's level of enzyme activity from their genotype. Aggregated across all genotypes, we find that increasing levels of G6PD deficiency are associated with decreasing risk of cerebral malaria, but with increased risk of severe malarial anaemia. Models of balancing selection based on these findings indicate that an evolutionary trade-off between different clinical outcomes of P. falciparum infection could have been a major cause of the high levels of G6PD polymorphism seen in human populations

Insights into malaria susceptibility using genome-wide data on 17,000 individuals from Africa, Asia and Oceania

The human genetic factors that affect resistance to infectious disease are poorly understood. Here we report a genome-wide association study in 17,000 severe malaria cases and population controls from 11 countries, informed by sequencing of family trios and by direct typing of candidate loci in an additional 15,000 samples. We identify five replicable associations with genome-wide levels of evidence including a newly implicated variant on chromosome 6. Jointly, these variants account for around one-tenth of the heritability of severe malaria, which we estimate as -23% using genome-wide genotypes. We interrogate available functional data and discover an erythroid-specific transcription start site underlying the known association in ATP2B4, but are unable to identify a likely causal mechanism at the chromosome 6 locus. Previously reported HLA associations do not replicate in these samples. This large dataset will provide a foundation for further research on thegenetic determinants of malaria resistance in diverse populations.Peer reviewe

A comparison of case-control and family-based association methods: the example of sickle-cell and malaria.

There has been much debate about the relative merits of population- and family-based strategies for testing genetic association, yet there is little empirical data that directly compare the two approaches. Here we compare case-control and transmission/disequilibrium test (TDT) study designs using a well-established genetic association, the protective effect of the sickle-cell trait against severe malaria. We find that the two methods give similar estimates of the level of protection (case-control odds ratio = 0.10, 95% confidence interval 0.03-0.23; family-based estimate of the odds ratio = 0.11, 95% confidence interval 0.04-0.25) and similar statistical significance of the result (case-control: chi2= 41.26, p= 10(-10), TDT: chi2= 39.06, p= 10(-10)) when 315 TDT cases are compared to 583 controls. We propose a family plus population control study design, which allows both case-control and TDT analysis of the cases. This combination is robust against the respective weaknesses of the case-control and TDT study designs, namely population structure and segregation distortion. The combined study design is especially cost-effective when cases are difficult to ascertain and, when the case-control and TDT results agree, offers greater confidence in the result

Investigation of malaria susceptibility determinants in the IFNG/IL26/IL22 genomic region.

Interferon-gamma, encoded by IFNG, is a key immunological mediator that is believed to play both a protective and a pathological role in malaria. Here, we investigate the relationship between IFNG variation and susceptibility to malaria. We began by analysing West African and European haplotype structure and patterns of linkage disequilibrium across a 100 kb genomic region encompassing IFNG and its immediate neighbours IL22 and IL26. A large case-control study of severe malaria in a West Africa population identified several weak associations with individual single-nucleotide polymorphisms in the IFNG and IL22 genes, and defined two IL22 haplotypes that are, respectively, associated with resistance and susceptibility. These data provide a starting point for functional and genetic analysis of the IFNG genomic region in malaria and other infectious and inflammatory conditions affecting African populations

Analysis of IL10 haplotypic associations with severe malaria

We investigated the association between severe malaria and genetic variation of IL10 in Gambian children, as several lines of evidence indicate that IL10 is protective against severe malaria and that IL10 production is genetically determined. We began by identifying five informative SNPs in the Gambian population that were genotyped in a combined case-control and intrafamilial study including 654 cases of severe malaria, 579 sets of parents and 459 ethnically matched controls. No significant associations were identified with individual SNPs. One haplotype of frequency 0.11 was strongly associated with protection against severe malaria in the case-control analysis (odds ratio 0.52, P=0.00002), but the transmission disequilibrium test in families showed no significant effect. These findings raise the question of whether IL10 associations with severe malaria might be confounded by foetal survival rates or other sources of transmission bias

Nucleotide diversity of the TNF gene region in an African village.

The wide variety of disease associations reported at the TNF locus raises the question of how much variation exists within a single population. To address this question, we sequenced the entire TNF gene in 72 chromosomes from healthy residents of a village in The Gambia, West Africa. We found 12 polymorphisms in 4393 nucleotides, of which five have not been previously described, giving an estimated nucleotide diversity (theta) of 5.6 x 10(-4). A significantly higher frequency of polymorphisms was found in the promoter region than in the coding region (8/1256 vs 0/882 nucleotides, P = 0.02). All polymorphisms with the exception of one rare allele were found to be present in Malawi, which is both geographically and genetically distant from The Gambia. Genotyping of 424 Gambian and 121 Malawian adults showed a significant frequency difference between the two populations for eight of the 12 polymorphisms, but the average fixation index across the variable sites was relatively low (F(ST) = 0.007). We conclude that, at the TNF locus, the nucleotide diversity found within a single African village is similar to the global value for human autosomal genes sampled across different continents