Geraldine M. Clarke

A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits

To identify genetic factors influencing quantitative traits of biomedical importance, we conducted a genome-wide association study in 8,842 samples from population-based cohorts recruited in Korea. For height and body mass index, most variants detected overlapped those reported in European samples. For the other traits examined, replication of promising GWAS signals in 7,861 independent Korean samples identified six previously unknown loci. For pulse rate, signals reaching genome-wide significance mapped to chromosomes 1q32 (rs12731740, P = 2.9 × 10−9) and 6q22 (rs12110693, P = 1.6 × 10−9), with the latter ∼400 kb from the coding sequence of GJA1. For systolic blood pressure, the most compelling association involved chromosome 12q21 and variants near the ATP2B1 gene (rs17249754, P = 1.3 × 10−7). For waist-hip ratio, variants on chromosome 12q24 (rs2074356, P = 7.8 × 10−12) showed convincing associations, although no regional transcript has strong biological candidacy. Finally, we identified two loci influencing bone mineral density at multiple sites. On chromosome 7q31, rs7776725 (within the FAM3C gene) was associated with bone density at the radius (P = 1.0 × 10−11), tibia (P = 1.6 × 10−6) and heel (P = 1.9 × 10−10). On chromosome 7p14, rs1721400 (mapping close to SFRP4, a frizzled protein gene) showed consistent associations at the same three sites (P = 2.2 × 10−3, P = 1.4 × 10−7 and P = 6.0 × 10−4, respectively). This large-scale GWA analysis of well-characterized Korean population-based samples highlights previously unknown biological pathways.

Reappraisal of known malaria resistance loci in a large multicenter study

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. We tested 55 SNPs in 27 loci previously reported to associate with severe malaria. There was evidence of association at P < 1 × 10−4 with the HBB, ABO, ATP2B4, G6PD and CD40LG loci, but previously reported associations at 22 other loci did not replicate in the multicenter analysis. The large sample size made it possible to identify authentic genetic effects that are heterogeneous across populations or phenotypes, with 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 anemia. 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.

Basic statistical analysis in genetic case-control studies

This protocol describes how to perform basic statistical analysis in a population-based genetic association case-control study. The steps described involve the (i) appropriate selection of measures of association and relevance of disease models; (ii) appropriate selection of tests of association; (iii) visualization and interpretation of results; (iv) consideration of appropriate methods to control for multiple testing; and (v) replication strategies. Assuming no previous experience with software such as PLINK, R or Haploview, we describe how to use these popular tools for handling single-nucleotide polymorphism data in order to carry out tests of association and visualize and interpret results. This protocol assumes that data quality assessment and control has been performed, as described in a previous protocol, so that samples and markers deemed to have the potential to introduce bias to the study have been identified and removed. Study design, marker selection and quality control of case-control studies have also been discussed in earlier protocols. The protocol should take ∼1 h to complete.

Aspects of observing and claiming allele flips in association studies.

Significant allele flipping, where associations for the same disease occur at opposite alleles of the same bi‐allelic locus, is increasing. But when is a significant allele flip genuine? We address the statistical issues of claiming and observing genuine allele flips in actual samples. We show that unless an allele flip is genuine, the probability of observing a significant allele flip in samples ascertained similarly from a common population is negligible. We derive expressions for the expected values of commonly used measures of association, which confirm previous findings that the underlying mechanism of a genuine allele flip is variation in the haplotype frequencies and show further how this variation interacts with variation in the genetic effects to impact allele flipping. We show that for association testing at proxy SNPs, common in genome‐wide association studies, variation in haplotype frequencies must coincide with a reversal in the sign of linkage disequilibrium (LD) to trigger genuine allele flips. Using HapMap data and r, rather than r2, to highlight previously unobserved effects, we show that unless genetic effects are large, variation in LD is unlikely to cause genuine allele flips in samples drawn from the same population. However, as populations diverge, it is an increasingly viable cause of a genuine allele flip for sufficiently large genetic effect and/or sample sizes. We conclude that evidence of variation in local patterns of LD, ancestral composition of study samples, and environmental exposures between study populations can provide compelling practical evidence in defense of a genuine allele flip. Genet. Epidemiol. 34: 266–274, 2010.

Resistance to malaria through structural variation of red blood cell invasion receptors

Structural variants are mapped that are correlated with a reduced risk of severe malaria. Pathogens select for genomic variants Large-scale deletions and duplications of genes, referred to as structural variants (SVs), are common within the human genome and have been linked to disease. Examining a genomic region that appears to confer a selective benefit, Leffler et al. used fine mapping to identify a specific SV that reduces the risk of severe malaria by an estimated 40% (see the Perspective by Winzeler). Data from African individuals revealed that populations harbor different SVs in this region. Furthermore, by dissecting a highly complex genomic region, the authors identified the likely causal element. This element encodes hybrid genes that affect glycophorin proteins, which are used by the malarial parasite in infection and are associated with resistance to severe disease. Science, this issue p. eaam6393; see also p. 1122 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. A structural variant creating hybrid glycophorin genes is associated with protection from severe malaria. The reference haplotype carries three glycophorin genes, two of which (GYPA and GYPB) are expressed as proteins on the red blood cell surface. The malaria-protective haplotype carries five glycophorin genes, including two hybrid genes that encode the Dantu blood group antigen and are composed of a GYPB extracellular domain and GYPA intracellular domain. These glycophorins serve as receptors for malaria-parasite ligands during red blood cell invasion. The malaria parasite Plasmodium falciparum invades human red blood cells by a series of interactions between host and parasite surface proteins. By analyzing genome sequence data from human populations, including 1269 individuals from sub-Saharan Africa, we identify a diverse array of large copy-number variants affecting the host invasion receptor genes GYPA and GYPB. We find that a nearby association with severe malaria is explained by a complex structural rearrangement involving the loss of GYPB and gain of two GYPB-A hybrid genes, which encode a serologically distinct blood group antigen known as Dantu. This variant reduces the risk of severe malaria by 40% and has recently increased in frequency in parts of Kenya, yet it appears to be absent from west Africa. These findings link structural variation of red blood cell invasion receptors with natural resistance to severe malaria.

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 effect has 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. DOI: http://dx.doi.org/10.7554/eLife.15085.001

Glucose-6-phosphate dehydrogenase deficiency and the risk of malaria and other diseases in children in Kenya: a case-control and a cohort study

Summary Background The global prevalence of X-linked glucose-6-phosphate dehydrogenase (G6PD) deficiency is thought to be a result of selection by malaria, but epidemiological studies have yielded confusing results. We investigated the relationships between G6PD deficiency and both malaria and non-malarial illnesses among children in Kenya. Methods We did this study in Kilifi County, Kenya, where the G6PD c.202T allele is the only significant cause of G6PD deficiency. We tested the associations between G6PD deficiency and severe and complicated Plasmodium falciparum malaria through a case-control study of 2220 case and 3940 control children. Cases were children aged younger than 14 years, who visited the high dependency ward of Kilifi County Hospital with severe malaria between March 1, 1998, and Feb 28, 2010. Controls were children aged between 3–12 months who were born within the same study area between August 2006, and September 2010. We assessed the association between G6PD deficiency and both uncomplicated malaria and other common diseases of childhood in a cohort study of 752 children aged younger than 10 years. Participants of this study were recruited from a representative sample of households within the Ngerenya and Chonyi areas of Kilifi County between Aug 1, 1998, and July 31, 2001. The primary outcome measure for the case-control study was the odds ratio for hospital admission with severe malaria (computed by logistic regression) while for the cohort study it was the incidence rate ratio for uncomplicated malaria and non-malaria illnesses (computed by Poisson regression), by G6PD deficiency category. Findings 2863 (73%) children in the control group versus 1643 (74%) in the case group had the G6PD normal genotype, 639 (16%) versus 306 (14%) were girls heterozygous for G6PD c.202T, and 438 (11%) versus 271 (12%) children were either homozygous girls or hemizygous boys. Compared with boys and girls without G6PD deficiency, we found significant protection from severe malaria (odds ratio [OR] 0·82, 95% CI 0·70–0·97; p=0·020) among G6PD c.202T heterozygous girls but no evidence for protection among G6PD c.202T hemizygous boys and homozygous girls (OR 1·18, 0·99–1·40; p=0·056). Median follow-up for the mild disease cohort study was 2·24 years (IQR 2·22–2·85). G6PD c.202T had no effect on other common diseases of childhood in heterozygous girls (incidence rate ratio 0·98, 95% CI 0·86–1·11; p=0·82) or homozygous girls or hemizygous boys (0·93, 0·82–1·04; p=0·25), with the sole exception of a marginally significant increase in the incidence of helminth infections among heterozygous girls. Interpretation Heterozygous girls might be the driving force for the positive selection of G6PD deficiency alleles. Further studies are needed to definitively establish the mechanisms by which G6PD deficiency confers an advantage against malaria in heterozygous individuals. Such studies could lead to the development of new treatments. Funding Wellcome Trust, UK Medical Research Council, European Union, and Foundation for the National Institutes of Health (as part of the Bill & Melinda Gates Grand Challenges in Global Health Initiative).

The challenge of detecting epistasis (G x G interactions): Genetic Analysis Workshop 16.

Interest is increasing in epistasis as a possible source of the unexplained variance missed by genome‐wide association studies. The Genetic Analysis Workshop 16 Group 9 participants evaluated a wide variety of classical and novel analytical methods for detecting epistasis, in both the statistical and machine learning paradigms, applied to both real and simulated data. Because the magnitude of epistasis is clearly relative to scale of penetrance, and therefore to some extent, to the choice of model framework, it is not surprising that strong interactions under one model might be minimized or even disappear entirely under a different modeling framework. Genet. Epidemiol. 33 (Suppl. 1):S58–S67, 2009.

A Comparison of Sample Size and Power in Case-Only Association Studies of Gene-Environment Interaction

Assuming continuous, normally distributed environmental and categorical genotype variables, the authors compare 6 case-only designs for tests of association in gene-environment interaction. Novel tests modeling the environmental variable as either the response or the predictor and allowing a genetic variable with multiallelic variants are included. The authors show that tests imposing the same genotypic pattern of inheritance perform similarly regardless of whether genotype is the response variable or the predictor variable. The novel tests using the genetic variable as the response variable are advantageous because they are robust to non-normally distributed environmental exposures. Dominance deviance—deviation from additivity in the main or interaction effects—is key to test performance: When it is zero or modest, tests searching for a trend with increasing risk alleles are optimal; when it is large, tests for genotypic effects are optimal. However, the authors show that dominance deviance is attenuated when it is observed at a proxy locus, which is common in genome-wide association studies, so large dominance deviance is likely to be rare. The authors conclude that the trend test is the appropriate tool for large-scale association scans where the true gene-environment interaction model is unknown. The common practice of assuming a dominant pattern of inheritance can cause serious losses of power in the presence of any recessive, or modest dominant, effects.

A Flexible Approach for the Analysis of Rare Variants Allowing for a Mixture of Effects on Binary or Quantitative Traits

Multiple rare variants either within or across genes have been hypothesised to collectively influence complex human traits. The increasing availability of high throughput sequencing technologies offers the opportunity to study the effect of rare variants on these traits. However, appropriate and computationally efficient analytical methods are required to account for collections of rare variants that display a combination of protective, deleterious and null effects on the trait. We have developed a novel method for the analysis of rare genetic variation in a gene, region or pathway that, by simply aggregating summary statistics at each variant, can: (i) test for the presence of a mixture of effects on a trait; (ii) be applied to both binary and quantitative traits in population-based and family-based data; (iii) adjust for covariates to allow for non-genetic risk factors and; (iv) incorporate imputed genetic variation. In addition, for preliminary identification of promising genes, the method can be applied to association summary statistics, available from meta-analysis of published data, for example, without the need for individual level genotype data. Through simulation, we show that our method is immune to the presence of bi-directional effects, with no apparent loss in power across a range of different mixtures, and can achieve greater power than existing approaches as long as summary statistics at each variant are robust. We apply our method to investigate association of type-1 diabetes with imputed rare variants within genes in the major histocompatibility complex using genotype data from the Wellcome Trust Case Control Consortium.