A genomic and evolutionary approach reveals non-genetic drug resistance in malaria

Background: Drug resistance remains a major public health challenge for malaria treatment and eradication. Individual loci associated with drug resistance to many antimalarials have been identified, but their epistasis with other resistance mechanisms has not yet been elucidated. Results: We previously described two mutations in the cytoplasmic prolyl-tRNA synthetase (cPRS) gene that confer resistance to halofuginone. We describe here the evolutionary trajectory of halofuginone resistance of two independent drug resistance selections in Plasmodium falciparum. Using this novel methodology, we discover an unexpected non-genetic drug resistance mechanism that P. falciparum utilizes before genetic modification of the cPRS. P. falciparum first upregulates its proline amino acid homeostasis in response to halofuginone pressure. We show that this non-genetic adaptation to halofuginone is not likely mediated by differential RNA expression and precedes mutation or amplification of the cPRS gene. By tracking the evolution of the two drug resistance selections with whole genome sequencing, we further demonstrate that the cPRS locus accounts for the majority of genetic adaptation to halofuginone in P. falciparum. We further validate that copy-number variations at the cPRS locus also contribute to halofuginone resistance. Conclusions: We provide a three-step model for multi-locus evolution of halofuginone drug resistance in P. falciparum. Informed by genomic approaches, our results provide the first comprehensive view of the evolutionary trajectory malaria parasites take to achieve drug resistance. Our understanding of the multiple genetic and non-genetic mechanisms of drug resistance informs how we will design and pair future anti-malarials for clinical use. Electronic supplementary material The online version of this article (doi:10.1186/s13059-014-0511-2) contains supplementary material, which is available to authorized users

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

Background: Mounting 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. Results: We 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 parasite's 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. Conclusions: We 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.Organismic and Evolutionary BiologyStem Cell and Regenerative BiologyOther Research Uni

Release of Sequestered Malaria Parasites upon Injection of a Glycosaminoglycan

Severe human malaria is attributable to an excessive sequestration of Plasmodium falciparum–infected and uninfected erythrocytes in vital organs. Strains of P. falciparum that form rosettes and employ heparan sulfate as a host receptor are associated with development of severe forms of malaria. Heparin, which is similar to heparan sulfate in that it is composed of the same building blocks, was previously used in the treatment of severe malaria, but it was discontinued due to the occurrence of serious side effects such as intracranial bleedings. Here we report to have depolymerized heparin by periodate treatment to generate novel glycans (dGAG) that lack anticoagulant-activity. The dGAGs disrupt rosettes, inhibit merozoite invasion of erythrocytes and endothelial binding of P. falciparum–infected erythrocytes in vitro, and reduce sequestration in in vivo models of severe malaria. An intravenous injection of dGAGs blocks up to 80% of infected erythrocytes from binding in the micro-vasculature of the rat and releases already sequestered parasites into circulation. P. falciparum–infected human erythrocytes that sequester in the non-human primate Macaca fascicularis were similarly found to be released in to the circulation upon a single injection of 500 μg of dGAG. We suggest dGAGs to be promising candidates for adjunct therapy in severe malaria

Default Pathway of var2csa Switching and Translational Repression in Plasmodium falciparum

Antigenic variation is a subtle process of fundamental importance to the survival of a microbial pathogen. In Plasmodium falciparum malaria, PfEMP1 is the major variable antigen and adhesin expressed at the surface of the infected erythrocyte, which is encoded for by members of a family of 60 var-genes. Peri-nuclear repositioning and epigenetic mechanisms control their mono-allelic expression. The switching of PfEMP1 depends in part on variable transition rates and short-lived immune responses to shared minor epitopes. Here we show var-genes to switch to a common gene that is highly transcribed, but sparsely translated into PfEMP1 and not expressed at the erythrocyte surface. Highly clonal and adhesive P. falciparum, which expressed distinct var-genes and the corresponding PfEMP1s at onset, were propagated without enrichment or panning. The parasites successively and spontaneously switched to transcribe a shared var-gene (var2csa) matched by the loss of PfEMP1 surface expression and host cell-binding. The var2csa gene repositioned in the peri-nuclear area upon activation, away from the telomeric clusters and heterochromatin to transcribe spliced, full-length RNA. Despite abundant transcripts, the level of intracellular PfEMP1 was low suggesting post-transcriptional mechanisms to partake in protein expression. In vivo, off-switching and translational repression may constitute one pathway, among others, coordinating PfEMP1 expression

A Novel DBL-Domain of the P. falciparum 332 Molecule Possibly Involved in Erythrocyte Adhesion

Plasmodium falciparum malaria is brought about by the asexual stages of the parasite residing in human red blood cells (RBC). Contact between the erythrocyte surface and the merozoite is the first step for successful invasion and proliferation of the parasite. A number of different pathways utilised by the parasite to adhere and invade the host RBC have been characterized, but the complete biology of this process remains elusive. We here report the identification of an open reading frame (ORF) representing a hitherto unknown second exon of the Pf332 gene that encodes a cysteine-rich polypeptide with a high degree of similarity to the Duffy-binding-like (DBL) domain of the erythrocyte-binding-ligand (EBL) family. The sequence of this DBL-domain is conserved and expressed in all parasite clones/strains investigated. In addition, the expression level of Pf332 correlates with proliferation efficiency of the parasites in vitro. Antibodies raised against the DBL-domain are able to reduce the invasion efficiency of different parasite clones/strains. Analysis of the DBL-domain revealed its ability to bind to uninfected human RBC, and moreover demonstrated association with the iRBC surface. Thus, Pf332 is a molecule with a potential role to support merozoite invasion. Due to the high level of conservation in sequence, the novel DBL-domain of Pf332 is of possible importance for development of novel anti-malaria drugs and vaccines

Regulatory Hotspots in the Malaria Parasite Genome Dictate Transcriptional Variation

The determinants of transcriptional regulation in malaria parasites remain elusive. The presence of a well-characterized gene expression cascade shared by different Plasmodium falciparum strains could imply that transcriptional regulation and its natural variation do not contribute significantly to the evolution of parasite drug resistance. To clarify the role of transcriptional variation as a source of stain-specific diversity in the most deadly malaria species and to find genetic loci that dictate variations in gene expression, we examined genome-wide expression level polymorphisms (ELPs) in a genetic cross between phenotypically distinct parasite clones. Significant variation in gene expression is observed through direct co-hybridizations of RNA from different P. falciparum clones. Nearly 18% of genes were regulated by a significant expression quantitative trait locus. The genetic determinants of most of these ELPs resided in hotspots that are physically distant from their targets. The most prominent regulatory locus, influencing 269 transcripts, coincided with a Chromosome 5 amplification event carrying the drug resistance gene, pfmdr1, and 13 other genes. Drug selection pressure in the Dd2 parental clone lineage led not only to a copy number change in the pfmdr1 gene but also to an increased copy number of putative neighboring regulatory factors that, in turn, broadly influence the transcriptional network. Previously unrecognized transcriptional variation, controlled by polymorphic regulatory genes and possibly master regulators within large copy number variants, contributes to sweeping phenotypic evolution in drug-resistant malaria parasites

Antimalarial Exposure Delays Plasmodium falciparum Intra-Erythrocytic Cycle and Drives Drug Transporter Genes Expression

BACKGROUND: Multi-drug resistant Plasmodium falciparum is a major obstacle to malaria control and is emerging as a complex phenomenon. Mechanisms of drug evasion based on the intracellular extrusion of the drug and/or modification of target proteins have been described. However, cellular mechanisms related with metabolic activity have also been seen in eukaryotic systems, e.g. cancer cells. Recent observations suggest that such mechanism may occur in P. falciparum. METHODOLOGY/PRINCIPAL FINDINGS: We therefore investigated the effect of mefloquine exposure on the cell cycle of three P. falciparum clones (3D7, FCB, W2) with different drug susceptibilities, while investigating in parallel the expression of four genes coding for confirmed and putative drug transporters (pfcrt, pfmdr1, pfmrp1 and pfmrp2). Mefloquine induced a previously not described dose and clone dependent delay in the intra-erythrocytic cycle of the parasite. Drug impact on cell cycle progression and gene expression was then merged using a non-linear regression model to determine specific drug driven expression. This revealed a mild, but significant, mefloquine driven gene induction up to 1.5 fold. CONCLUSIONS/SIGNIFICANCE: Both cell cycle delay and induced gene expression represent potentially important mechanisms for parasites to escape the effect of the antimalarial drug

Adaptive Copy Number Evolution in Malaria Parasites

Copy number polymorphism (CNP) is ubiquitous in eukaryotic genomes, but the degree to which this reflects the action of positive selection is poorly understood. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive CNP. We provide compelling evidence that gch1 CNP is an adaptive consequence of selection by antifolate drugs, which target enzymes downstream in this pathway. (1) We compared gch1 CNP in parasites from Thailand (strong historical antifolate selection) with those from neighboring Laos (weak antifolate selection). Two percent of chromosomes had amplified copy number in Laos, while 72% carried multiple (2–11) copies in Thailand, and differentiation exceeded that observed at 73 synonymous SNPs. (2) We found five amplicon types containing one to greater than six genes and spanning 1 to >11 kb, consistent with parallel evolution and strong selection for this gene amplification. gch1 was the only gene occurring in all amplicons suggesting that this locus is the target of selection. (3) We observed reduced microsatellite variation and increased linkage disequilibrium (LD) in a 900-kb region flanking gch1 in parasites from Thailand, consistent with rapid recent spread of chromosomes carrying multiple copies of gch1. (4) We found that parasites bearing dhfr-164L, which causes high-level resistance to antifolate drugs, carry significantly (p = 0.00003) higher copy numbers of gch1 than parasites bearing 164I, indicating functional association between genes located on different chromosomes but linked in the same biochemical pathway. These results demonstrate that CNP at gch1 is adaptive and the associations with dhfr-164L strongly suggest a compensatory function. More generally, these data demonstrate how selection affects multiple enzymes in a single biochemical pathway, and suggest that investigation of structural variation may provide a fast-track to locating genes underlying adaptation

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

<p>Abstract</p> <p>Background</p> <p>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, <it>Plasmodium falciparum</it>, to identify and analyze the inheritance of 170 genome-wide CNVs.</p> <p>Results</p> <p>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 <it>de novo </it>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.</p> <p>Conclusions</p> <p>CNVs are a significant source of segregating and <it>de novo </it>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.</p

Comparative Transcriptional and Genomic Analysis of Plasmodium falciparum Field Isolates

Mechanisms for differential regulation of gene expression may underlie much of the phenotypic variation and adaptability of malaria parasites. Here we describe transcriptional variation among culture-adapted field isolates of Plasmodium falciparum, the species responsible for most malarial disease. It was found that genes coding for parasite protein export into the red cell cytosol and onto its surface, and genes coding for sexual stage proteins involved in parasite transmission are up-regulated in field isolates compared with long-term laboratory isolates. Much of this variability was associated with the loss of small or large chromosomal segments, or other forms of gene copy number variation that are prevalent in the P. falciparum genome (copy number variants, CNVs). Expression levels of genes inside these segments were correlated to that of genes outside and adjacent to the segment boundaries, and this association declined with distance from the CNV boundary. This observation could not be explained by copy number variation in these adjacent genes. This suggests a local-acting regulatory role for CNVs in transcription of neighboring genes and helps explain the chromosomal clustering that we observed here. Transcriptional co-regulation of physical clusters of adaptive genes may provide a way for the parasite to readily adapt to its highly heterogeneous and strongly selective environment