Gene copy number and malaria biology

Prevalence of pvmrp1 Polymorphisms and Its Contribution to Antimalarial Response

As more sporadic cases of chloroquine resistance occur (CQR) in Plasmodium vivax (P. vivax) malaria, molecular markers have become an important tool to monitor the introduction and spread of drug resistance. P. vivax multidrug resistance-associated protein 1 (PvMRP1), as one of the members of the ATP-binding cassette (ABC) transporters, may modulate this phenotype. In this study, we investigated the gene mutations and copy number variations (CNVs) in the pvmrp1 in 102 P. vivax isolates from China, the Republic of Korea (ROK), Myanmar, Papua New Guinea (PNG), Pakistan, the Democratic People’s Republic of Korea (PRK), and Cambodia. And we also obtained 72 available global pvmrp1 sequences deposited in the PlasmoDB database to investigate the genetic diversity, haplotype diversity, natural selection, and population structure of pvmrp1. In total, 29 single nucleotide polymorphisms reflected in 23 non-synonymous, five synonymous mutations and one gene deletion were identified, and CNVs were found in 2.9% of the isolates. Combined with the antimalarial drug susceptibility observed in the previous in vitro assays, except the prevalence of S354N between the two CQ sensitivity categories revealed a significant difference, no genetic mutations or CNVs associated with drug sensitivity were found. The genetic polymorphism analysis of 166 isolates worldwide found that the overall nucleotide diversity (π) of pvmrp1 was 0.0011, with 46 haplotypes identified (Hd = 0.9290). The ratio of non-synonymous to synonymous mutations (dn/ds = 0.5536) and the neutrality tests statistic Fu and Li’s D* test (Fu and Li’s D* = −3.9871, p < 0.02) suggests that pvmrp1 had evolved under a purifying selection. Due to geographical differences, genetic differentiation levels of pvmrp1 in different regions were different to some extent. Overall, this study provides a new idea for finding CQR molecular monitoring of P. vivax and provides more sequences of pvmrp1 in Asia for subsequent research. However, further validation is still needed through laboratory and epidemiological field studies of P. vivax samples from more regions.

Extrachromosomal DNA amplicons in antimalarial-resistant Plasmodium falciparum

Extrachromosomal (ec) DNAs are genetic elements that exist separately from the genome. Since ecDNA can carry beneficial genes, they are a powerful adaptive mechanism in cancers and many pathogens. For the first time, we report ecDNA contributing to antimalarial resistance in Plasmodium falciparum, the most virulent human malaria parasite. Using pulse field gel electrophoresis combined with PCR-based copy number analysis, we detected two ecDNA elements that differ in migration and structure. Entrapment in the electrophoresis well and low susceptibility to exonucleases revealed that the biologically relevant ecDNA element is large and complex in structure. Using deep sequencing, we show that ecDNA originates from the chromosome and expansion of an ecDNA-specific sequence may improve its segregation or expression. We speculate that ecDNA is maintained using established mechanisms due to shared characteristics with the mitochondrial genome. Implications of ecDNA discovery in this organism are wide-reaching due to the potential for new strategies to target resistance development.

The contribution of extrachromosomal DNA to genome plasticity in malaria parasites

Malaria caused by the protozoan parasite Plasmodium falciparum continues to impose significant morbidity and mortality, despite substantial investment into drug and vaccine development and deployment. Underlying the resilience of this parasite is its remarkable ability to undergo genome modifications, thus, providing parasite populations with extensive genetic variability that accelerates selection of drug resistance and limits the efficacy of most vaccines. This genome plasticity is rooted in the mechanisms of DNA repair that parasites employ to maintain genome integrity, a process skewed toward homologous recombination through the evolutionary loss of classical nonhomologous end joining. Repair of DNA double-strand breaks have been shown to enable "shuffling" of antigen-encoding gene sequences to vastly increase antigen diversity and to enable copy number expansion of genes that contribute to drug resistance. The latter phenomenon has been proposed to be a major contributor to the rise of resistance to several classes of antimalarial drugs. In this issue of Molecular Microbiology, McDaniels and colleagues add yet another mechanism that malaria parasites use to reduce drug susceptibility by demonstrating that P. falciparum can maintain expanded arrays of drug resistance cassettes as stably replicating, circular, extrachromosomal DNAs, thus, expanding genome plasticity beyond the parasite's 14 nuclear chromosomes.

Transcriptional variation in malaria parasites: why and how

Transcriptional differences enable the generation of alternative phenotypes from the same genome. In malaria parasites, transcriptional plasticity plays a major role in the process of adaptation to fluctuations in the environment. Multiple studies with culture-adapted parasites and field isolates are starting to unravel the different transcriptional alternatives available to Plasmodium falciparum and the underlying molecular mechanisms. Here we discuss how epigenetic variation, directed transcriptional responses and also genetic changes that affect transcript levels can all contribute to transcriptional variation and, ultimately, parasite survival. Some transcriptional changes are driven by stochastic events. These changes can occur spontaneously, resulting in heterogeneity within parasite populations that provides the grounds for adaptation by dynamic natural selection. However, transcriptional changes can also occur in response to external cues. A better understanding of the mechanisms that the parasite has evolved to alter its transcriptome may ultimately contribute to the design of strategies to combat malaria to which the parasite cannot adapt.

Targeting Kinetoplastid and Apicomplexan Thymidylate Biosynthesis as an Antiprotozoal Strategy

Kinetoplastid and apicomplexan parasites comprise a group of protozoans responsible for human diseases, with a serious impact on human health and the socioeconomic growth of developing countries. Chemotherapy is the main option to control these pathogenic organisms and nucleotide metabolism is considered a promising area for the provision of antimicrobial therapeutic targets. Impairment of thymidylate (dTMP) biosynthesis severely diminishes the viability of parasitic protozoa and the absence of enzymatic activities specifically involved in the formation of dTMP (e.g. dUTPase, thymidylate synthase, dihydrofolate reductase or thymidine kinase) results in decreased deoxythymidine triphosphate (dTTP) levels and the so-called thymineless death. In this process, the ratio of deoxyuridine triphosphate (dUTP) versus dTTP in the cellular nucleotide pool has a crucial role. A high dUTP/dTTP ratio leads to uracil misincorporation into DNA, the activation of DNA repair pathways, DNA fragmentation and eventually cell death. The essential character of dTMP synthesis has stimulated interest in the identification and development of drugs that specifically block the biochemical steps involved in thymine nucleotide formation. Here, we review the available literature in relation to drug discovery studies targeting thymidylate biosynthesis in kinetoplastid (genera Trypanosoma and Leishmania) and apicomplexan (Plasmodium spp and Toxoplasma gondii) protozoans. The most relevant findings concerning novel inhibitory molecules with antiparasitic activity against these human pathogens are presented herein.

Complex DNA structures trigger copy number variation across the Plasmodium falciparum genome

Antimalarial resistance is a major obstacle in the eradication of the human malaria parasite, Plasmodium falciparum. Genome amplifications, a type of DNA copy number variation (CNV), facilitate overexpression of drug targets and contribute to parasite survival. Long monomeric A/T tracks are found at the breakpoints of many Plasmodium resistance-conferring CNVs. We hypothesize that other proximal sequence features, such as DNA hairpins, act with A/T tracks to trigger CNV formation. By adapting a sequence analysis pipeline to investigate previously reported CNVs, we identified breakpoints in 35 parasite clones with near single base-pair resolution. Using parental genome sequence, we predicted the formation of stable hairpins within close proximity to all future breakpoint locations. Especially stable hairpins were predicted to form near five shared breakpoints, establishing that the initiating event could have occurred at these sites. Further in-depth analyses defined characteristics of these 'trigger sites' across the genome and detected signatures of error-prone repair pathways at the breakpoints. We propose that these two genomic signals form the initial lesion (hairpins) and facilitate microhomology-mediated repair (A/T tracks) that lead to CNV formation across this highly repetitive genome. Targeting these repair pathways in P. falciparum may be used to block adaptation to antimalarial drugs.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Assessment of copy number variation in genes related to drug resistance in Plasmodium vivax and Plasmodium falciparum isolates from the Brazilian Amazon and a systematic review of the literature

Background

Parasite resistance to anti-malarials represents a great obstacle for malaria elimination. The majority of studies have investigated the association between single-nucleotide polymorphisms (SNPs) and drug resistance; however, it is becoming clear that the copy number variation (CNV) is also associated with this parasite phenotype. To provide a baseline for molecular surveillance of anti-malarial drug resistance in the Brazilian Amazon, the present study characterized the genetic profile of both markers in the most common genes associated with drug resistance in Plasmodium falciparum and Plasmodium vivax isolates. Additionally, these data were compared to data published elsewhere applying a systematic review of the literature published over a 20-year time period.

Methods

The genomic DNA of 67 patients infected by P. falciparum and P. vivax from three Brazilian States was obtained between 2002 and 2012. CNV in P. falciparum multidrug resistance gene-1 (pfmdr1), GTP cyclohydrolase 1 (pfgch1) and P. vivax multidrug resistance gene-1 (pvmdr1) were assessed by real-time PCR assays. SNPs in the pfmdr1 and pfcrt genes were assessed by PCR–RFLP. A literature search for studies that analysed CNP in the same genes of P. falciparum and P. vivax was conducted between May 2014 and March 2017 across four databases.

Results

All analysed samples of P. falciparum carried only one copy of pfmdr1 or pfgch1. Although the pfcrt K76T polymorphism, a determinant of CQ resistance, was present in all samples genotyped, the pfmdr1 N86Y was absent. For P. vivax isolates, an amplification rate of 20% was found for the pvmdr1 gene. The results of the study are in agreement with the low amplification rates for pfmdr1 gene evidenced in the Americas and Africa, while higher rates have been described in Southeast Asia. For P. vivax, very low rates of amplification for pvmdr1 have been described worldwide, with exceptions in French Guiana, Cambodia, Thailand and Brazil.

Conclusions

The present study was the first to evaluate gch1 CNV in P. falciparum isolates from Brazil, showing an absence of amplification of this gene more than 20 years after the withdrawal of the Brazilian antifolates therapeutic scheme. Furthermore, the rate of pvmdr1 amplification was significantly higher than that previously reported for isolates circulating in Northern Brazil.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-1806-z) contains supplementary material, which is available to authorized users.

Analysis of anti-malarial resistance markers in pfmdr1 and pfcrt across Southeast Asia in the Tracking Resistance to Artemisinin Collaboration

Background

Declining anti-malarial efficacy of artemisinin-based combination therapy, and reduced Plasmodium falciparum susceptibility to individual anti-malarials are being documented across an expanding area of Southeast Asia (SEA). Genotypic markers complement phenotypic studies in assessing the efficacy of individual anti-malarials.

Methods

The markers pfmdr1 and pfcrt were genotyped in parasite samples obtained in 2011–2014 at 14 TRAC (Tracking Resistance to Artemisinin Collaboration) sites in mainland Southeast Asia using a combination of PCR and next-generation sequencing methods.

Results

Pfmdr1 amplification, a marker of mefloquine and lumefantrine resistance, was highly prevalent at Mae Sot on the Thailand–Myanmar border (59.8% of isolates) and common (more than 10%) at sites in central Myanmar, eastern Thailand and western Cambodia; however, its prevalence was lower than previously documented in Pailin, western Cambodia. The pfmdr1 Y184F mutation was common, particularly in and around Cambodia, and the F1226Y mutation was found in about half of samples in Mae Sot. The functional significance of these two mutations remains unclear. Other previously documented pfmdr1 mutations were absent or very rare in the region. The pfcrt mutation K76T associated with chloroquine resistance was found in 98.2% of isolates. The CVIET haplotype made up 95% or more of isolates in western SEA while the CVIDT haplotype was common (30–40% of isolates) in north and northeastern Cambodia, southern Laos, and southern Vietnam.

Conclusions

These findings generate cause for concern regarding the mid-term efficacy of artemether–lumefantrine in Myanmar, while the absence of resistance-conferring pfmdr1 mutations and SVMNT pfcrt haplotypes suggests that amodiaquine could be an efficacious component of anti-malarial regimens in SEA.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1598-6) contains supplementary material, which is available to authorized users.

Dissecting malaria biology and epidemiology using population genetics and genomics

Molecular approaches have an increasingly recognized utility in surveillance of malaria parasite populations, not only in defining prevalence and incidence with higher sensitivity than traditional methods, but also in monitoring local and regional parasite transmission patterns. In this review, we provide an overview of population genetic and genomic studies of human-infecting Plasmodium species, highlighting recent advances in the field. In accordance with the renewed impetus for malaria eradication, many studies are now using genetic and genomic epidemiology to support local evidence-based intervention strategies. Microsatellite genotyping remains a popular approach for both Plasmodium falciparum and Plasmodium vivax. However, with the increasing availability of whole genome sequencing data enabling effective single nucleotide polymorphism-based panels tailored to a given study question and setting, this approach is gaining popularity. The availability of new reference genomes for Plasmodium malariae and Plasmodium ovale should see a surge in similar molecular studies on these currently neglected species. Genomic studies are revealing new insights into important adaptive mechanisms of the parasite including antimalarial drug resistance. The advent of new methodologies such as selective whole genome amplification for dealing with extensive human DNA in low density field isolates should see genome-wide approaches becoming routine for parasite surveillance once the economic costs outweigh the current cost benefits of targeted approaches.

Indels, structural variation, and recombination drive genomic diversity in Plasmodium falciparum

The malaria parasite Plasmodium falciparum has a great capacity for evolutionary adaptation to evade host immunity and develop drug resistance. Current understanding of parasite evolution is impeded by the fact that a large fraction of the genome is either highly repetitive or highly variable and thus difficult to analyze using short-read sequencing technologies. Here, we describe a resource of deep sequencing data on parents and progeny from genetic crosses, which has enabled us to perform the first genome-wide, integrated analysis of SNP, indel and complex polymorphisms, using Mendelian error rates as an indicator of genotypic accuracy. These data reveal that indels are exceptionally abundant, being more common than SNPs and thus the dominant mode of polymorphism within the core genome. We use the high density of SNP and indel markers to analyze patterns of meiotic recombination, confirming a high rate of crossover events and providing the first estimates for the rate of non-crossover events and the length of conversion tracts. We observe several instances of meiotic recombination within copy number variants associated with drug resistance, demonstrating a mechanism whereby fitness costs associated with resistance mutations could be compensated and greater phenotypic plasticity could be acquired.

Population Structure Shapes Copy Number Variation in Malaria Parasites

If copy number variants (CNVs) are predominantly deleterious, we would expect them to be more efficiently purged from populations with a large effective population size (Ne) than from populations with a small Ne. Malaria parasites (Plasmodium falciparum) provide an excellent organism to examine this prediction, because this protozoan shows a broad spectrum of population structures within a single species, with large, stable, outbred populations in Africa, small unstable inbred populations in South America and with intermediate population characteristics in South East Asia. We characterized 122 single-clone parasites, without prior laboratory culture, from malaria-infected patients in seven countries in Africa, South East Asia and South America using a high-density single-nucleotide polymorphism/CNV microarray. We scored 134 high-confidence CNVs across the parasite exome, including 33 deletions and 102 amplifications, which ranged in size from <500 bp to 59 kb, as well as 10,107 flanking, biallelic single-nucleotide polymorphisms. Overall, CNVs were rare, small, and skewed toward low frequency variants, consistent with the deleterious model. Relative to African and South East Asian populations, CNVs were significantly more common in South America, showed significantly less skew in allele frequencies, and were significantly larger. On this background of low frequency CNV, we also identified several high-frequency CNVs under putative positive selection using an FST outlier analysis. These included known adaptive CNVs containing rh2b and pfmdr1, and several other CNVs (e.g., DNA helicase and three conserved proteins) that require further investigation. Our data are consistent with a significant impact of genetic structure on CNV burden in an important human pathogen.

Zinc finger nuclease-based double-strand breaks attenuate malaria parasites and reveal rare microhomology-mediated end joining

Background

Genome editing of malaria parasites is key to the generation of live attenuated parasites used in experimental vaccination approaches. DNA repair in Plasmodium generally occurs only through homologous recombination. This has been used to generate transgenic parasites that lack one to three genes, leading to developmental arrest in the liver and allowing the host to launch a protective immune response. While effective in principle, this approach is not safe for use in humans as single surviving parasites can still cause disease. Here we use zinc-finger nucleases to generate attenuated parasite lines lacking an entire chromosome arm, by a timed induction of a double-strand break. Rare surviving parasites also allow the investigation of unconventional DNA repair mechanisms in a rodent malaria parasite.

Results

A single, zinc-finger nuclease-induced DNA double-strand break results in the generation of attenuated parasite lines that show varying degrees of developmental arrest, protection efficacy in an immunisation regime and safety, depending on the timing of zinc-finger nuclease expression within the life cycle. We also identify DNA repair by microhomology-mediated end joining with as little as four base pairs, resulting in surviving parasites and thus breakthrough infections.

Conclusions

Malaria parasites can repair DNA double-strand breaks with surprisingly small mini-homology domains located across the break point. Timely expression of zinc-finger nucleases could be used to generate a new generation of attenuated parasite lines lacking hundreds of genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0811-1) contains supplementary material, which is available to authorized users.

The molecular basis of antifolate resistance in Plasmodium falciparum: looking beyond point mutations

Drugs that target the folate-synthesis pathway have a long history of effectiveness against a variety of pathogens. As antimalarials, the antifolates were safe and well tolerated, but resistance emerged quickly and has persisted even with decreased drug pressure. The primary determinants of resistance in Plasmodium falciparum are well-described point mutations in the enzymes dihydropteroate synthase and dihydrofolate reductase targeted by the combination sulfadoxine-pyrimethamine. Recent work has highlighted the contributions of additional parasite adaptation to antifolate resistance. In fact, the evolution of antifolate-resistant parasites is multifaceted and complex. Gene amplification of the first enzyme in the parasite folate synthesis pathway, GTP-cyclohydrolase, is strongly associated with resistant parasites and potentially contributes to persistence of resistant parasites. Further understanding of how parasites adjust flux through the folate pathway is important to the further development of alternative agents targeting this crucial synthesis pathway.

A better resolution for integrating methods for monitoring Plasmodium falciparum resistance to antimalarial drugs

Effective chemotherapy is the mainstay of malaria control. However, resistance of falciparum malaria to antimalarial drugs compromised the efforts to eliminate the disease and led to the resurgence of malaria epidemics. Three main approaches are used to monitor antimalarial drug efficacy and drug resistance; namely, in vivo trials, in vitro/ex vivo assays and molecular markers of drug resistance. Each approach has its implications of use as well as its advantages and drawbacks. Therefore, there is a need to use an integrated approach that would give the utmost effect to detect resistance as early as its emergence and to track it once spread. Such integration becomes increasingly needed in the era of artemisinin-based combination therapy as a forward action to deter resistance. The existence of regional and global networks for the standardization of methodology, provision of high quality reagents for the assessment of antimalarial drug resistance and dissemination of open-access data would help in approaching an integrated resistance surveillance system on a global scale.

A survey of innovation through duplication in the reduced genomes of twelve parasites

We characterize the prevalence, distribution, divergence, and putative functions of detectable two-copy paralogs and segmental duplications in the Apicomplexa, a phylum of parasitic protists. Apicomplexans are mostly obligate intracellular parasites responsible for human and animal diseases (e.g. malaria and toxoplasmosis). Gene loss is a major force in the phylum. Genomes are small and protein-encoding gene repertoires are reduced. Despite this genomic streamlining, duplications and gene family amplifications are present. The potential for innovation introduced by duplications is of particular interest. We compared genomes of twelve apicomplexans across four lineages and used orthology and genome cartography to map distributions of duplications against genome architectures. Segmental duplications appear limited to five species. Where present, they correspond to regions enriched for multi-copy and species-specific genes, pointing toward roles in adaptation and innovation. We found a phylum-wide association of duplications with dynamic chromosome regions and syntenic breakpoints. Trends in the distribution of duplicated genes indicate that recent, species-specific duplicates are often tandem while most others have been dispersed by genome rearrangements. These trends show a relationship between genome architecture and gene duplication. Functional analysis reveals: proteases, which are vital to a parasitic lifecycle, to be prominent in putative recent duplications; a pair of paralogous genes in Toxoplasma gondii previously shown to produce the rate-limiting step in dopamine synthesis in mammalian cells, a possible link to the modification of host behavior; and phylum-wide differences in expression and subcellular localization, indicative of modes of divergence. We have uncovered trends in multiple modes of duplicate divergence including sequence, intron content, expression, subcellular localization, and functions of putative recent duplicates that highlight the role of duplications in the continuum of forces that have shaped these genomes.

Single-cell genomics for dissection of complex malaria infections

Most malaria infections contain complex mixtures of distinct parasite lineages. These multiple-genotype infections (MGIs) impact virulence evolution, drug resistance, intra-host dynamics, and recombination, but are poorly understood. To address this we have developed a single-cell genomics approach to dissect MGIs. By combining cell sorting and whole-genome amplification (WGA), we are able to generate high-quality material from parasite-infected red blood cells (RBCs) for genotyping and next-generation sequencing. We optimized our approach through analysis of >260 single-cell assays. To quantify accuracy, we decomposed mixtures of known parasite genotypes and obtained highly accurate (>99%) single-cell genotypes. We applied this validated approach directly to infections of two major malaria species, Plasmodium falciparum, for which long term culture is possible, and Plasmodium vivax, for which no long-term culture is feasible. We demonstrate that our single-cell genomics approach can be used to generate parasite genome sequences directly from patient blood in order to unravel the complexity of P. vivax and P. falciparum infections. These methods open the door for large-scale analysis of within-host variation of malaria infections, and reveal information on relatedness and drug resistance haplotypes that is inaccessible through conventional sequencing of infections.

Asexual populations of the human malaria parasite, Plasmodium falciparum, use a two-step genomic strategy to acquire accurate, beneficial DNA amplifications

Malaria drug resistance contributes to up to a million annual deaths. Judicious deployment of new antimalarials and vaccines could benefit from an understanding of early molecular events that promote the evolution of parasites. Continuous in vitro challenge of Plasmodium falciparum parasites with a novel dihydroorotate dehydrogenase (DHODH) inhibitor reproducibly selected for resistant parasites. Genome-wide analysis of independently-derived resistant clones revealed a two-step strategy to evolutionary success. Some haploid blood-stage parasites first survive antimalarial pressure through fortuitous DNA duplications that always included the DHODH gene. Independently-selected parasites had different sized amplification units but they were always flanked by distant A/T tracks. Higher level amplification and resistance was attained using a second, more efficient and more accurate, mechanism for head-to-tail expansion of the founder unit. This second homology-based process could faithfully tune DNA copy numbers in either direction, always retaining the unique DNA amplification sequence from the original A/T-mediated duplication for that parasite line. Pseudo-polyploidy at relevant genomic loci sets the stage for gaining additional mutations at the locus of interest. Overall, we reveal a population-based genomic strategy for mutagenesis that operates in human stages of P. falciparum to efficiently yield resistance-causing genetic changes at the correct locus in a successful parasite. Importantly, these founding events arise with precision; no other new amplifications are seen in the resistant haploid blood stage parasite. This minimizes the need for meiotic genetic cleansing that can only occur in sexual stage development of the parasite in mosquitoes.

Malaria parasites kill up to a million people around the world every year. Emergence of resistance to drugs remains a key obstacle against elimination of malaria. In the laboratory, parasites can efficiently acquire resistance to experimental antimalarials by changing DNA at the target locus. This happens efficiently even for an antimalarial that the parasite has never encountered in a clinical setting. In this study, we formally demonstrate how parasites achieve this feat: first, individual parasites in a population of millions randomly amplify large regions of DNA between short sequence repeats of adenines (A) or thymines (T) that are peppered throughout the malaria parasite genome. The rare lucky parasite that amplifies DNA coding for the target of the antimalarial, along with dozens of its neighboring genes, gains an evolutionary advantage and survives. In a second step, to withstand increasing drug pressure and to achieve higher levels of resistance, each parasite line makes additional copies of this region. This second expansion does not rely on the random A/T-based DNA rearrangements but, instead, a more precise amplification mechanism that retains the unique signature of co-amplified genes created earlier in each parasite. Generation of multiple copies of the target genes in the parasite genome may be the beginning of other beneficial changes for the parasite, including the future acquisition of mutations.

Direct evidence for the adaptive role of copy number variation on antifolate susceptibility in Plasmodium falciparum

Resistance to antimalarials targeting the folate pathway is widespread. GTP-cyclohydrolase (gch1), the first enzyme in this pathway, exhibits extensive copy number variation (CN) in parasite isolates from areas with a history of longstanding antifolate use. Increased CN of gch1 is associated with a greater number of point mutations in enzymes targeted by the antifolates, pyrimethamine and sulphadoxine. While these observations suggest that increases in gch1 CN are an adaptation to drug pressure, changes in CN have not been experimentally demonstrated to directly alter drug susceptibility. To determine if changes in gch1 expression alone modify pyrimethamine sensitivity, we manipulated gch1 CN in several parasite lines to test the effect on drug sensitivity. We report that increases in gch1 CN alter pyrimethamine resistance in most parasites lines. However we find evidence of a detrimental effect of very high levels of gch1 overexpression in parasite lines with high endogenous levels of gch1 expression, revealing the importance of maintaining balance in the folate pathway and implicating changes in gch1 expression in preserving proper metabolic flux. This work expands our understanding of parasite adaptation to drug pressure and provides a possible mechanism for how specific mutations become fixed within parasite populations.

Detection of copy number variation and single nucleotide polymorphisms in genes involved in drug resistance and other phenotypic traits in P. falciparum clinical isolates collected from Uganda

There is an increasing interest in mapping the genes of pathogens which underlie important phenotypic traits such as virulence and drug resistance. The Plasmodium falciparum genome exhibits sequence variation that contributes to the pathogenic mechanisms of the parasite. Determining the prevalence of resistance markers could provide a prediction about drug efficacy. Copy number polymorphism (CNP) of genes has been shown to influence important parasite phenotypes. In this work, CNPs within genes involved in drug resistance and other phenotypic traits namely P. falciparum multidrug resistance 1 (pfmdr-1), GTP cyclo hydrolase (gch1), Ring infected erythrocyte surface antigen precursor (resa) and a hypothetical protein coding gene were analyzed by quantitative real time-polymerase reaction (qRT-PCR) among clinical isolates collected from Uganda. The pfmdr-1 codons 86 and 1246 and P. falciparum chloroquine resistance (pfcrt) codon 76 were genotyped for single nucleotide polymorphisms (SNPs) by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP), and the proportion of resistance associated mutations were determined among mild and severe malaria cases using the chi-square test. Forty and 42 P. falciparum isolates collected from children with mild and severe malaria respectively were analyzed for CNPs. Seventy five and 81 P. falciparum isolates from children with mild or severe malaria were analyzed for SNPs. No pfmdr-1, gch1 or novel gene amplifications were identified among the P. falciparum clinical isolates. Although chloroquine was officially withdrawn from policy use since 7 years, all P. falciparum isolates presented the associated pfcrt K76T mutation, whatever the clinical status and no specific mutation in the pfmdr-1 gene was associated with disease type. In conclusion, this study provides baseline measures for continued surveillance for changes in copy number and SNP types among genes implicated in drug resistance and other important phenotypes that may have a potential role in parasite virulence mechanisms or drug treatment outcomes.

Outwitting evolution: fighting drug-resistant TB, malaria, and HIV

Although caused by vastly different pathogens, the world's three most serious infectious diseases, tuberculosis, malaria, and HIV-1 infection, share the common problem of drug resistance. The pace of drug development has been very slow for tuberculosis and malaria and rapid for HIV-1. But for each disease, resistance to most drugs has appeared quickly after the introduction of the drug. Learning how to manage and prevent resistance is a major medical challenge that requires an understanding of the evolutionary dynamics of each pathogen. This Review summarizes the similarities and differences in the evolution of drug resistance for these three pathogens.

The landscape of inherited and de novo copy number variants in a Plasmodium falciparum genetic cross

BACKGROUND

Copy number is a major source of genome variation with important evolutionary implications. Consequently, it is essential to determine copy number variant (CNV) behavior, distributions and frequencies across genomes to understand their origins in both evolutionary and generational time frames. We use comparative genomic hybridization (CGH) microarray and the resolution provided by a segregating population of cloned progeny lines of the malaria parasite, Plasmodium falciparum, to identify and analyze the inheritance of 170 genome-wide CNVs.

RESULTS

We describe CNVs in progeny clones derived from both Mendelian (i.e. inherited) and non-Mendelian mechanisms. Forty-five CNVs were present in the parent lines and segregated in the progeny population. Furthermore, extensive variation that did not conform to strict Mendelian inheritance patterns was observed. 124 CNVs were called in one or more progeny but in neither parent: we observed CNVs in more than one progeny clone that were not identified in either parent, located more frequently in the telomeric-subtelomeric regions of chromosomes and singleton de novo CNVs distributed evenly throughout the genome. Linkage analysis of CNVs revealed dynamic copy number fluctuations and suggested mechanisms that could have generated them. Five of 12 previously identified expression quantitative trait loci (eQTL) hotspots coincide with CNVs, demonstrating the potential for broad influence of CNV on the transcriptional program and phenotypic variation.

CONCLUSIONS

CNVs are a significant source of segregating and de novo genome variation involving hundreds of genes. Examination of progeny genome segments provides a framework to assess the extent and possible origins of CNVs. This segregating genetic system reveals the breadth, distribution and dynamics of CNVs in a surprisingly plastic parasite genome, providing a new perspective on the sources of diversity in parasite populations.

Recent Progress in Functional Genomic Research in Plasmodium falciparum

With the completion and near completion of many malaria parasite genome-sequencing projects, efforts are now being directed to a better understanding of gene functions and to the discovery of vaccine and drug targets. Inter- and intraspecies comparisons of the parasite genomes will provide invaluable insights into parasite evolution, virulence, drug resistance, and immune invasion. Genome-wide searches for loci under various selection pressures may lead to discovery of genes conferring drug resistance or encoding for protective antigens. In addition, the Plasmodium falciparum genome sequence provides the basis for the development of various microarrays to monitor gene expression and to detect nucleotide substitution and deletion/amplification. Genome-wide profiling of the parasite proteome, chromatin modification, and nucleosome position also depend on availability of the parasite genome. In this brief review, we will highlight some recent advances and studies in characterizing gene function and related phenotype in P. falciparum that were made possible by the genome sequence, particularly the development of a genome-wide diversity map and various high-throughput genotyping methods for genome-wide association studies (GWAS).

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

The Plasmodium falciparum parasite's 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.

Malaria infection with the human pathogen Plasmodium falciparum results in almost a million deaths each year, mostly in African children. Efforts to eliminate malaria are underway, but the parasite is adept at eluding both the human immune response and antimalarial treatments. Thus, it is important to understand how the parasite becomes resistant to drugs and to develop strategies to overcome resistance mechanisms. Toward this end, we used population genetic strategies to identify genetic loci that contribute to parasite adaptation and to identify candidate genes involved in drug resistance. We examined over 17,000 genetic variants across the parasite genome in over 50 strains in which we also measured responses to many known antimalarial compounds. We found a number of genetic loci showing signs of recent natural selection and a number of loci potentially involved in modulating the parasite's response to drugs. We further demonstrated that one of the novel candidate genes (PF10_0355) modulates resistance to the antimalarial compounds halofantrine, mefloquine, and lumefantrine. Overall, this study confirms that we can use genome-wide approaches to identify clinically relevant genes and demonstrates through functional testing that at least one of these candidate genes is indeed involved in antimalarial drug resistance.

How artemisinin-containing combination therapies slow the spread of antimalarial drug resistance

Antimalarial drug therapies containing artemisinins, 'ACTs', have become the mainstay for treating uncomplicated malaria in endemic countries. This is a major public health achievement requiring substantial political, financial and scientific input. The most compelling scientific argument for ACT deployment employed a very simple basic rationale that emphasised their role in slowing the origin of drug resistance while largely neglecting the additional role(s) of ACTs in slowing or preventing the spread of resistance once it has arisen. Recent reports suggest that early stages of resistance to artemisinins and/or its partner drugs could be occurring, thus it is timely to briefly review exactly how ACTs slow the origin and spread of resistance and to interpret the threat of resistance within this context.

Deamplification of pfmdr1-containing amplicon on chromosome 5 in Plasmodium falciparum is associated with reduced resistance to artelinic acid in vitro

Amplification of the Plasmodium falciparum multidrug resistance 1 gene (pfmdr1) has been implicated in multidrug resistance, including in vitro resistance to artelinic acid (AL). The stability and fitness of having multiple copies of pfmdr1 are important factors due to their potential effects on the resistance phenotype of parasites. These factors were investigated by using an AL-resistant line of P. falciparum (W2AL80) and clones originating from W2AL80. A rapid reduction in pfmdr1 copy number (CN) was observed in the uncloned W2AL80 line; 63% of this population reverted to a CN of <3 without exposure to the drug. Deamplification of the pfmdr1 amplicon was then determined in three clones, each initially containing three copies of pfmdr1. Interestingly, two outcomes were observed during 3 months without drug pressure. In one clone, parasites with fewer than 3 copies of pfmdr1 emerged rapidly. In two other clones, the reversion was significantly delayed. In all subclones, the reduction in pfmdr1 CN involved the deamplification of the entire amplicon (19 genes). Importantly, deamplification of the pfmdr1 amplicon resulted in partial reversal of resistance to AL and increased susceptibility to mefloquine. These results demonstrate that multiple copies of the pfmdr1-containing amplicon in AL-resistant parasites are unstable when drug pressure is withdrawn and have practical implications for the maintenance and spread of parasites resistant to artemisinin derivatives.

Genetic changes during laboratory propagation: copy number At the reticulocyte-binding protein 1 locus of Plasmodium falciparum

Comparative genomic hybridization studies have revealed elevated copy number (CN) at the reticulocyte-binding protein 1 gene (PfRh1) in fast growing lab-adapted parasites, while genetic manipulation demonstrates a causal link between cell invasion and PfRh1 CN. We therefore examined PfRh1 copy number variation (CNV) in 202 single clone parasite isolates from four countries to quantify the extent of CNV within natural populations. Surprisingly, we found that no natural parasite infections showed elevated CN. In contrast, 4/28 independent laboratory reference strains show elevated CN. One possibility is that amplification of PfRh1 (or neighboring loci) is selected during laboratory culture. In the case of FCR3 group of parasites, clone trees show that PfRh1 amplification arose in laboratory lines following establishment in culture. These data show that CNV at PfRh1 is rare or non-existent in natural populations, but can arise during laboratory propagation. We conclude that PfRh1 CNV is not an important determinant of gene expression, cell invasion or growth rate in natural parasite populations.