Sarah K. Volkman

A genome-wide map of diversity in Plasmodium falciparum

Genetic variation allows the malaria parasite Plasmodium falciparum to overcome chemotherapeutic agents, vaccines and vector control strategies and remain a leading cause of global morbidity and mortality. Here we describe an initial survey of genetic variation across the P. falciparum genome. We performed extensive sequencing of 16 geographically diverse parasites and identified 46,937 SNPs, demonstrating rich diversity among P. falciparum parasites (π = 1.16 × 10−3) and strong correlation with gene function. We identified multiple regions with signatures of selective sweeps in drug-resistant parasites, including a previously unidentified 160-kb region with extremely low polymorphism in pyrimethamine-resistant parasites. We further characterized 54 worldwide isolates by genotyping SNPs across 20 genomic regions. These data begin to define population structure among African, Asian and American groups and illustrate the degree of linkage disequilibrium, which extends over relatively short distances in African parasites but over longer distances in Asian parasites. We provide an initial map of genetic diversity in P. falciparum and demonstrate its potential utility in identifying genes subject to recent natural selection and in understanding the population genetics of this parasite.

Amplification of pfmdr1 associated with mefloquine and halofantrine resistance in Plasmodium falciparum from Thailand

Drug resistance in Plasmodium falciparum is an expanding problem in most endemic areas. Recent studies have suggested the potential involvement of genes in the MDR gene family in resistance to quinoline-containing compounds in P. falciparum. In this study a molecular analysis of pfmdr 1 in recent isolates from Thailand was done (1) to further examine the role of pfmdr 1 in drug-resistant isolates and (2) to examine the reported association of pfmdr 1 intragenic alleles and chloroquine resistance. Most of the isolates (10 of 11) were resistant to all compounds tested. Analysis of pfmdr 1 revealed an apparent association between increased gene copy number and increased level of expression of pfmdr 1 and decreased susceptibility to mefloquine and halofantrine. Sequence analysis of pfmdr 1 in these isolates revealed no association of intragenic alleles with chloroquine resistance.

Use of high-density tiling microarrays to identify mutations globally and elucidate mechanisms of drug resistance in Plasmodium falciparum.

BackgroundThe identification of genetic changes that confer drug resistance or other phenotypic changes in pathogens can help optimize treatment strategies, support the development of new therapeutic agents, and provide information about the likely function of genes. Elucidating mechanisms of phenotypic drug resistance can also assist in identifying the mode of action of uncharacterized but potent antimalarial compounds identified in high-throughput chemical screening campaigns against Plasmodium falciparum.ResultsHere we show that tiling microarrays can detect de novo a large proportion of the genetic changes that differentiate one genome from another. We show that we detect most single nucleotide polymorphisms or small insertion deletion events and all known copy number variations that distinguish three laboratory isolates using readily accessible methods. We used the approach to discover mutations that occur during the selection process after transfection. We also elucidated a mechanism by which parasites acquire resistance to the antimalarial fosmidomycin, which targets the parasite isoprenoid synthesis pathway. Our microarray-based approach allowed us to attribute in vitro derived fosmidomycin resistance to a copy number variation event in the pfdxr gene, which enables the parasite to overcome fosmidomycin-mediated inhibition of isoprenoid biosynthesis.ConclusionsWe show that newly emerged single nucleotide polymorphisms can readily be detected and that malaria parasites can rapidly acquire gene amplifications in response to in vitro drug pressure. The ability to define comprehensively genetic variability in P. falciparum with a single overnight hybridization creates new opportunities to study parasite evolution and improve the treatment and control of malaria.

Genetic diversity and protective efficacy of the RTS,S/AS01 malaria vaccine

BACKGROUND The RTS,S/AS01 vaccine targets the circumsporozoite protein of Plasmodium falciparum and has partial protective efficacy against clinical and severe malaria disease in infants and children. We investigated whether the vaccine efficacy was specific to certain parasite genotypes at the circumsporozoite protein locus. METHODS We used polymerase chain reaction-based next-generation sequencing of DNA extracted from samples from 4985 participants to survey circumsporozoite protein polymorphisms. We evaluated the effect that polymorphic positions and haplotypic regions within the circumsporozoite protein had on vaccine efficacy against first episodes of clinical malaria within 1 year after vaccination. RESULTS In the per-protocol group of 4577 RTS,S/AS01-vaccinated participants and 2335 control-vaccinated participants who were 5 to 17 months of age, the 1-year cumulative vaccine efficacy was 50.3% (95% confidence interval [CI], 34.6 to 62.3) against clinical malaria in which parasites matched the vaccine in the entire circumsporozoite protein C-terminal (139 infections), as compared with 33.4% (95% CI, 29.3 to 37.2) against mismatched malaria (1951 infections) (P=0.04 for differential vaccine efficacy). The vaccine efficacy based on the hazard ratio was 62.7% (95% CI, 51.6 to 71.3) against matched infections versus 54.2% (95% CI, 49.9 to 58.1) against mismatched infections (P=0.06). In the group of infants 6 to 12 weeks of age, there was no evidence of differential allele-specific vaccine efficacy. CONCLUSIONS These results suggest that among children 5 to 17 months of age, the RTS,S vaccine has greater activity against malaria parasites with the matched circumsporozoite protein allele than against mismatched malaria. The overall vaccine efficacy in this age category will depend on the proportion of matched alleles in the local parasite population; in this trial, less than 10% of parasites had matched alleles. (Funded by the National Institutes of Health and others.).

A global transcriptional analysis of Plasmodium falciparum malaria reveals a novel family of telomere-associated lncRNAs

BackgroundMounting evidence suggests a major role for epigenetic feedback in Plasmodium falciparum transcriptional regulation. Long non-coding RNAs (lncRNAs) have recently emerged as a new paradigm in epigenetic remodeling. We therefore set out to investigate putative roles for lncRNAs in P. falciparum transcriptional regulation.ResultsWe used a high-resolution DNA tiling microarray to survey transcriptional activity across 22.6% of the P. falciparum strain 3D7 genome. We identified 872 protein-coding genes and 60 putative P. falciparum lncRNAs under developmental regulation during the parasites pathogenic human blood stage. Further characterization of lncRNA candidates led to the discovery of an intriguing family of lncRNA telomere-associated repetitive element transcripts, termed lncRNA-TARE. We have quantified lncRNA-TARE expression at 15 distinct chromosome ends and mapped putative transcriptional start and termination sites of lncRNA-TARE loci. Remarkably, we observed coordinated and stage-specific expression of lncRNA-TARE on all chromosome ends tested, and two dominant transcripts of approximately 1.5 kb and 3.1 kb transcribed towards the telomere.ConclusionsWe have characterized a family of 22 telomere-associated lncRNAs in P. falciparum. Homologous lncRNA-TARE loci are coordinately expressed after parasite DNA replication, and are poised to play an important role in P. falciparum telomere maintenance, virulence gene regulation, and potentially other processes of parasite chromosome end biology. Further study of lncRNA-TARE and other promising lncRNA candidates may provide mechanistic insight into P. falciparum transcriptional regulation.

Sequence-based association and selection scans identify drug resistance loci in the Plasmodium falciparum malaria parasite

Through rapid genetic adaptation and natural selection, the Plasmodium falciparum parasite—the deadliest of those that cause malaria—is able to develop resistance to antimalarial drugs, thwarting present efforts to control it. Genome-wide association studies (GWAS) provide a critical hypothesis-generating tool for understanding how this occurs. However, in P. falciparum, the limited amount of linkage disequilibrium hinders the power of traditional array-based GWAS. Here, we demonstrate the feasibility and power improvements gained by using whole-genome sequencing for association studies. We analyzed data from 45 Senegalese parasites and identified genetic changes associated with the parasites’ in vitro response to 12 different antimalarials. To further increase statistical power, we adapted a common test for natural selection, XP-EHH (cross-population extended haplotype homozygosity), and used it to identify genomic regions associated with resistance to drugs. Using this sequence-based approach and the combination of association and selection-based tests, we detected several loci associated with drug resistance. These loci included the previously known signals at pfcrt, dhfr, and pfmdr1, as well as many genes not previously implicated in drug-resistance roles, including genes in the ubiquitination pathway. Based on the success of the analysis presented in this study, and on the demonstrated shortcomings of array-based approaches, we argue for a complete transition to sequence-based GWAS for small, low linkage-disequilibrium genomes like that of P. falciparum.

Identification and functional validation of the novel antimalarial resistance locus PF10_0355 in Plasmodium falciparum.

The Plasmodium falciparum parasites ability to adapt to environmental pressures, such as the human immune system and antimalarial drugs, makes malaria an enduring burden to public health. Understanding the genetic basis of these adaptations is critical to intervening successfully against malaria. To that end, we created a high-density genotyping array that assays over 17,000 single nucleotide polymorphisms (∼1 SNP/kb), and applied it to 57 culture-adapted parasites from three continents. We characterized genome-wide genetic diversity within and between populations and identified numerous loci with signals of natural selection, suggesting their role in recent adaptation. In addition, we performed a genome-wide association study (GWAS), searching for loci correlated with resistance to thirteen antimalarials; we detected both known and novel resistance loci, including a new halofantrine resistance locus, PF10_0355. Through functional testing we demonstrated that PF10_0355 overexpression decreases sensitivity to halofantrine, mefloquine, and lumefantrine, but not to structurally unrelated antimalarials, and that increased gene copy number mediates resistance. Our GWAS and follow-on functional validation demonstrate the potential of genome-wide studies to elucidate functionally important loci in the malaria parasite genome.

Hybrid selection for sequencing pathogen genomes from clinical samples

We have adapted a solution hybrid selection protocol to enrich pathogen DNA in clinical samples dominated by human genetic material. Using mock mixtures of human and Plasmodium falciparum malaria parasite DNA as well as clinical samples from infected patients, we demonstrate an average of approximately 40-fold enrichment of parasite DNA after hybrid selection. This approach will enable efficient genome sequencing of pathogens from clinical samples, as well as sequencing of endosymbiotic organisms such as Wolbachia that live inside diverse metazoan phyla.

Modeling malaria genomics reveals transmission decline and rebound in Senegal.

Significance Traditional methods for estimating malaria transmission based on mosquito sampling are not standardized and are unavailable in many countries in sub-Saharan Africa. Such studies are especially difficult to implement when transmission is low, and low transmission is the goal of malaria elimination. Malaria-control efforts in Senegal have resulted in changes in population genomics evidenced by increased allele sharing among parasite genomes, often including genomic identity between independently sampled parasites. Fitting an epidemiological model to the observed data indicates falling transmission from 2006–2010 with a significant rebound in 2012–2013, an inference confirmed by incidence data. These results demonstrate that genomic approaches may help monitor transmission to assess initial and ongoing effectiveness of interventions to control malaria. To study the effects of malaria-control interventions on parasite population genomics, we examined a set of 1,007 samples of the malaria parasite Plasmodium falciparum collected in Thiès, Senegal between 2006 and 2013. The parasite samples were genotyped using a molecular barcode of 24 SNPs. About 35% of the samples grouped into subsets with identical barcodes, varying in size by year and sometimes persisting across years. The barcodes also formed networks of related groups. Analysis of 164 completely sequenced parasites revealed extensive sharing of genomic regions. In at least two cases we found first-generation recombinant offspring of parents whose genomes are similar or identical to genomes also present in the sample. An epidemiological model that tracks parasite genotypes can reproduce the observed pattern of barcode subsets. Quantification of likelihoods in the model strongly suggests a reduction of transmission from 2006–2010 with a significant rebound in 2012–2013. The reduced transmission and rebound were confirmed directly by incidence data from Thiès. These findings imply that intensive intervention to control malaria results in rapid and dramatic changes in parasite population genomics. The results also suggest that genomics combined with epidemiological modeling may afford prompt, continuous, and cost-effective tracking of progress toward malaria elimination.

Harnessing genomics and genome biology to understand malaria biology

Malaria is an important human disease and is the target of a global eradication campaign. New technological and informatics advancements in population genomics are being leveraged to identify genetic loci under selection in the malaria parasite and to find variants that are associated with key clinical phenotypes, such as drug resistance. This article provides a timely Review of how population-genetics-based strategies are being applied to Plasmodium falciparum both to identify genetic loci as key targets of interventions and to develop monitoring and surveillance tools that are crucial for the successful elimination and eradication of malaria.