Sarah Auburn

Genetic loci associated with delayed clearance of Plasmodium falciparum following artemisinin treatment in Southeast Asia

The recent emergence of artemisinin-resistant Plasmodium falciparum malaria in western Cambodia could threaten prospects for malaria elimination. Identification of the genetic basis of resistance would provide tools for molecular surveillance, aiding efforts to contain resistance. Clinical trials of artesunate efficacy were conducted in Bangladesh, in northwestern Thailand near the Myanmar border, and at two sites in western Cambodia. Parasites collected from trial participants were genotyped at 8,079 single nucleotide polymorphisms (SNPs) using a P. falciparum-specific SNP array. Parasite genotypes were examined for signatures of recent positive selection and association with parasite clearance phenotypes to identify regions of the genome associated with artemisinin resistance. Four SNPs on chromosomes 10 (one), 13 (two), and 14 (one) were significantly associated with delayed parasite clearance. The two SNPs on chromosome 13 are in a region of the genome that appears to be under strong recent positive selection in Cambodia. The SNPs on chromosomes 10 and 13 lie in or near genes involved in postreplication repair, a DNA damage-tolerance pathway. Replication and validation studies are needed to refine the location of loci responsible for artemisinin resistance and to understand the mechanism behind it; however, two SNPs on chromosomes 10 and 13 may be useful markers of delayed parasite clearance in surveillance for artemisinin resistance in Southeast Asia.

Positive selection of a CD36 nonsense variant in sub-Saharan Africa, but no association with severe malaria phenotypes

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.

Tumor necrosis factor and lymphotoxin-α polymorphisms and severe malaria 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.

Validating discovered Cis-acting regulatory genetic variants: application of an allele specific expression approach to HapMap populations.

Background Localising regulatory variants that control gene expression is a challenge for genome research. Several studies have recently identified non-coding polymorphisms associated with inter-individual differences in gene expression. These approaches rely on the identification of signals of association against a background of variation due to other genetic and environmental factors. A complementary approach is to use an Allele-Specific Expression (ASE) assay, which is more robust to the effects of environmental variation and trans-acting genetic factors. Methodology/Principal Findings Here we apply an ASE method which utilises heterozygosity within an individual to compare expression of the two alleles of a gene in a single cell. We used individuals from three HapMap population groups and analysed the allelic expression of genes with cis-regulatory regions previously identified using total gene expression studies. We were able to replicate the results in five of the six genes tested, and refined the cis- associated regions to a small number of variants. We also showed that by using multi-populations it is possible to refine the associated cis-effect DNA regions. Conclusions/Significance We discuss the efficacy and drawbacks of both total gene expression and ASE approaches in the discovery of cis-acting variants. We show that the ASE approach has significant advantages as it is a cleaner representation of cis-acting effects. We also discuss the implication of using different populations to map cis-acting regions and the importance of finding regulatory variants which contribute to human phenotypic variation.

A genetic association study in the Gambia using tagging polymorphisms in the major histocompatibility complex class III region implicates a HLA-B associated transcript 2 polymorphism in severe malaria susceptibility

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 11 haplotype tagging-polymorphisms for 11 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 2,162 case-controls demonstrated the first evidence of association between a BAT2 polymorphism (rs1046089) and severe malaria.

Association of the GNAS Locus with 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 19xa0kb GNAS region, four presented signals of association (Pxa0<xa00.05) with severe malaria. The strongest single-locus association demonstrated an odds ratio of 1.13 (1.05–1.21), Pxa0=xa00.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 [ORxa0=xa00.88 (0.81–0.96), Pxa0=xa00.005 and ORxa0=xa01.12 (1.03–1.20), Pxa0=xa00.005]. The evidence presented here indicates that the influence of G-alpha-s on erythrocyte invasion efficacy may, indeed, alter individual susceptibility to disease.

SnoopCGH: software for visualizing comparative genomic hybridization data

Summary: Array-based comparative genomic hybridization (CGH) technology is used to discover and validate genomic structural variation, including copy number variants, insertions, deletions and other structural variants (SVs). The visualization and summarization of the array CGH data outputs, potentially across many samples, is an important process in the identification and analysis of SVs. We have developed a software tool for SV analysis using data from array CGH technologies, which is also amenable to short-read sequence data. Availability and implementation: SnoopCGH is written in java and is available from http://snoopcgh.sourceforge.net/ Contact: jg10@sanger.ac.uk; tc5@sanger.ac.uk

Genetic determinants of anti-malarial acquired immunity in a large multi-centre study

AbstractBackgroundMany studies report associations between human genetic factors and immunity to malaria but few have been reliably replicated. These studies are usually country-specific, use small sample sizes and are nnot directly comparable due to differences in methodologies. This study brings together samples and data collected from multiple sites across Africa and Asia to use standardized methods to look for consistent genetic effects on anti-malarial antibody levels.MethodsSera, DNA samples and clinical data were collected from 13,299 individuals from ten sites in Senegal, Mali, Burkina Faso, Sudan, Kenya, Tanzania, and Sri Lanka using standardized methods. DNA was extracted and typed for 202 Single Nucleotide Polymorphisms with known associations to malaria or antibody production, and antibody levels to four clinical grade malarial antigens [AMA1, MSP1, MSP2, and (NANP)4] plus total IgE were measured by ELISA techniques. Regression models were used to investigate the associations of clinical and genetic factors with antibody levels.ResultsMalaria infection increased levels of antibodies to malaria antigens and, as expected, stable predictors of anti-malarial antibody levels included age, seasonality, location, and ethnicity. Correlations between antibodies to blood-stage antigens AMA1, MSP1 and MSP2 were higher between themselves than with antibodies to the (NANP)4 epitope of the pre-erythrocytic circumsporozoite protein, while there was little or no correlation with total IgE levels. Individuals with sickle cell trait had significantly lower antibody levels to all blood-stage antigens, and recessive homozygotes for CD36 (rs321198) had significantly lower anti-malarial antibody levels to MSP2.ConclusionAlthough the most significant finding with a consistent effect across sites was for sickle cell trait, its effect is likely to be via reducing a microscopically positive parasitaemia rather than directly on antibody levels. However, this study does demonstrate a framework for the feasibility of combining data from sites with heterogeneous malaria transmission levels across Africa and Asia with which to explore genetic effects on anti-malarial immunity.

VarB: A variation browsing and analysis tool for variants derived from next-generation sequencing data

Summary: There is an immediate need for tools to both analyse and visualize in real-time single-nucleotide polymorphisms, insertions and deletions, and other structural variants from new sequence file formats. We have developed VarB software that can be used to visualize variant call format files in real time, as well as identify regions under balancing selection and informative markers to differentiate user-defined groups (e.g. populations). We demonstrate its utility using sequence data from 50 Plasmodium falciparum isolates comprising two different continents and confirm known signals from genomic regions that contain important antigenic and anti-malarial drug-resistance genes. Availability and implementation: The C++-based software VarB and user manual are available from www.pathogenseq.org/varb. Contact: taane.clark@lshtm.ac.uk

A Bayesian approach using covariance of single nucleotide polymorphism data to detect differences in linkage disequilibrium patterns between groups of individuals

MOTIVATIONnQuantifying differences in linkage disequilibrium (LD) between sub-groups can highlight genetic regions or sites under selection and/or associated with disease, and may have utility in trans-ethnic mapping studies.nnnRESULTSnWe present a novel pseudo Bayes factor (PBF) approach that assess differences in covariance of genotype frequencies from single nucleotide polymorphism (SNP) data from a genome-wide study. The magnitude of the PBF reflects the strength of evidence for a difference, while accounting for the sample size and number of SNPs, without the requirement for permutation testing to establish statistical significance. Application of the PBF to HapMap and Gambian malaria SNP data reveals regional LD differences, some known to be under selection.nnnAVAILABILITY AND IMPLEMENTATIONnThe PBF approach has been implemented in the BALD (Bayesian analysis of LD differences) C++ software, and is available from http://homepages.lshtm.ac.uk/tgclark/downloads.