Starvation induces changes in abundance and small RNA cargo of extracellular vesicles released from Plasmodium falciparum infected red blood cells.

The lethal malaria parasite Plasmodium falciparum needs to constantly respond and adapt to changes within the human host in order to survive and transmit. One such change is composed of nutritional limitation, which is augmented with increased parasite loads and intimately linked to severe disease development. Extracellular vesicles released from infected red blood cells have been proposed as important mediators of disease pathogenesis and intercellular communication but whether important for the parasite response to nutritional availability is unknown. Therefore, we investigated the abundance and small RNA cargo of extracellular vesicles released upon short-term nutritional starvation of P. falciparum in vitro cultures. We show that primarily ring-stage parasite cultures respond to glucose and amino acid deprivation with an increased release of extracellular vesicles. Small RNA sequencing of these extracellular vesicles further revealed human miRNAs and parasitic tRNA fragments as the main constituent biotypes. Short-term starvations led to alterations in the transcriptomic profile, most notably in terms of the over-represented biotypes. These data suggest a potential role for extracellular vesicles released from P. falciparum infected red blood cells in the response to nutritional perturbations, their potential as prognostic biomarkers and point towards an evolutionary conserved role among protozoan parasites.

PD-1 expression on mouse intratumoral NK cells and its effects on NK cell phenotype.

Although PD-1 was shown to be a hallmark of T cells exhaustion, controversial studies have been reported on the role of PD-1 on NK cells. Here, we found by flow cytometry and single cell RNA sequencing analysis that PD-1 can be expressed on MHC class I-deficient tumor-infiltrating NK cells in vivo. We also demonstrate distinct alterations in the phenotype of PD-1-deficient NK cells and a more mature phenotype which might reduce their capacity to migrate and kill in vivo. Tumor-infiltrating NK cells that express PD-1 were highly associated with the expression of CXCR6. Furthermore, our results demonstrate that PD-L1 molecules in membranes of PD-1-deficient NK cells migrate faster than in NK cells from wild-type mice, suggesting that PD-1 and PD-L1 form cis interactions with each other on NK cells. These data demonstrate that there may be a role for the PD-1/PD-L1 axis in tumor-infiltrating NK cells in vivo.

Characterization of Metronidazole-Resistant Giardia intestinalis Lines by Comparative Transcriptomics and Proteomics.

Metronidazole (MTZ) is a clinically important antimicrobial agent that is active against both bacterial and protozoan organisms. MTZ has been used extensively for more than 60 years and until now resistance has been rare. However, a recent and dramatic increase in the number of MTZ resistant bacteria and protozoa is of great concern since there are few alternative drugs with a similarly broad activity spectrum. To identify key factors and mechanisms underlying MTZ resistance, we utilized the protozoan parasite Giardia intestinalis, which is commonly treated with MTZ. We characterized two in vitro selected, metronidazole resistant parasite lines, as well as one revertant, by analyzing fitness aspects associated with increased drug resistance and transcriptomes and proteomes. We also conducted a meta-analysis using already existing data from additional resistant G. intestinalis isolates. The combined data suggest that in vitro generated MTZ resistance has a substantial fitness cost to the parasite, which may partly explain why resistance is not widespread despite decades of heavy use. Mechanistically, MTZ resistance in Giardia is multifactorial and associated with complex changes, yet a core set of pathways involving oxidoreductases, oxidative stress responses and DNA repair proteins, is central to MTZ resistance in both bacteria and protozoa.

A Detailed Gene Expression Map of Giardia Encystation.

Giardia intestinalis is an intestinal protozoan parasite that causes diarrheal infections worldwide. A key process to sustain its chain of transmission is the formation of infectious cysts in the encystation process. We combined deep RNAseq of a broad range of encystation timepoints to produce a high-resolution gene expression map of Giardia encystation. This detailed transcriptomic map of encystation confirmed a gradual change of gene expression along the time course of encystation, showing the most significant gene expression changes during late encystation. Few genes are differentially expressed early in encystation, but the major cyst wall proteins CWP-1 and -2 are highly up-regulated already after 3.5 h encystation. Several transcription factors are sequentially up-regulated throughout the process, but many up-regulated genes at 7, 10, and 14 h post-induction of encystation have binding sites in the upstream regions for the Myb2 transcription factor, suggesting that Myb2 is a master regulator of encystation. We observed major changes in gene expression of several meiotic-related genes from 10.5 h of encystation to the cyst stage, and at 17.5 h encystation, there are changes in many different metabolic pathways and protein synthesis. Late encystation, 21 h to cysts, show extensive gene expression changes, most of all in VSP and HCMP genes, which are involved in antigenic variation, and genes involved in chromatin modifications. This high-resolution gene expression map of Giardia encystation will be an important tool in further studies of this important differentiation process.

Red blood cell blood group A antigen level affects the ability of heparin and PfEMP1 antibodies to disrupt Plasmodium falciparum rosettes.

BACKGROUND: The histo-blood group ABO system has been associated with adverse outcomes in COVID-19, thromboembolic diseases and Plasmodium falciparum malaria. An integral part of the severe malaria pathogenesis is rosetting, the adherence of parasite infected red blood cells (RBCs) to uninfected RBCs. Rosetting is influenced by the host's ABO blood group (Bg) and rosettes formed in BgA have previously been shown to be more resilient to disruption by heparin and shield the parasite derived surface antigens from antibodies. However, data on rosetting in weak BgA subgroups is scarce and based on investigations of relatively few donors. METHODS: An improved high-throughput flow cytometric assay was employed to investigate rosetting characteristics in an extensive panel of RBC donor samples of all four major ABO Bgs, as well as low BgA expressing samples. RESULTS: All non-O Bgs shield the parasite surface antigens from strain-specific antibodies towards P. falciparum erythrocyte membrane protein 1 (PfEMP1). A positive correlation between A-antigen levels on RBCs and rosette tightness was observed, protecting the rosettes from heparin- and antibody-mediated disruption. CONCLUSIONS: These results provide new insights into how the ABO Bg system affects the disease outcome and cautions against interpreting the results from the heterogeneous BgA phenotype as a single group in epidemiological and experimental studies.

GPR43 regulates marginal zone B-cell responses to foreign and endogenous antigens.

Marginal zone (MZ) B cells are innate-like B cells that produce polyreactive antibodies with an affinity for microbial molecular patterns and carbohydrate ligands. MZ B cells have been shown to be important in mediating immunity to various bacteria including Streptococcus pneumoniae and are also implicated in inflammatory syndromes including lupus erythematosus. The intestinal microbiota is responsible for producing short-chain fatty acids, which can regulate immune cell function by several mechanisms including ligation of the G-protein-coupled receptor (GPR)43. Herein, we show that MZ B cells express Gpr43 messenger RNA and that the absence of this receptor impacts on MZ B-cell surface marker expression and antibody production. In T-cell-independent responses to the hapten 4-hydroxy-3-nitrophenylacetic acid (NP), mice deficient in GPR43 displayed higher serum titers of NP-specific antibodies. Moreover, in response to a pneumococcal polysaccharide vaccine, GPR43-deficient mice developed robust serum antibody responses and had markedly increased numbers of splenic antibody-secreting cells, compared with control mice. Finally, serum immunoglobulin M autoantibodies to double-stranded DNA and phosphatidylcholine were increased in resting 10-15-week-old mice lacking GPR43. Taken together, mice lacking GPR43 have heightened antibody responses to T-cell-independent antigens, which may be a result of impaired regulation of MZ B cells.

Levels of human proteins in plasma associated with acute paediatric malaria.

BACKGROUND: The intimate interaction between the pathophysiology of the human host and the biology of the Plasmodium falciparum parasite results in a wide spectrum of disease outcomes in malaria. Development of severe disease is associated with a progressively augmented imbalance in pro- and anti-inflammatory responses to high parasite loads and sequestration of parasitized erythrocytes. Although these phenomena collectively constitute common denominators for the wide variety of discrete severe malaria manifestations, the mechanistic rationales behind discrepancies in outcome are poorly understood. Exploration of the human pathophysiological response by variations in protein profiles in plasma presents an excellent opportunity to increase the understanding. This is ultimately required for better prediction, prevention and treatment of malaria, which is essential for ongoing elimination and eradication efforts. RESULTS: An affinity proteomics approach was used to analyse 541 paediatric plasma samples collected from community controls and patients with mild or severe malaria in Rwanda. Protein profiles were generated with an antibody-based suspension bead array containing 255 antibodies targetting 115 human proteins. Here, 57 proteins were identified with significantly altered levels (adjusted p-values < 0.001) in patients with malaria compared to controls. From these, the 27 most significant proteins (adjusted p-values < 10-14) were selected for a stringent analysis approach. Here, 24 proteins showed elevated levels in malaria patients and included proteins involved in acute inflammatory response as well as cell adhesion. The remaining three proteins, also implicated in immune regulation and cellular adhesivity, displayed lower abundance in malaria patients. In addition, 37 proteins (adjusted p-values < 0.05) were identified with increased levels in patients with severe compared to mild malaria. This set includes, proteins involved in tissue remodelling and erythrocyte membrane proteins. Collectively, this approach has been successfully used to identify proteins both with known and unknown association with different stages of malaria. CONCLUSION: In this study, a high-throughput affinity proteomics approach was used to find protein profiles in plasma linked to P. falciparum infection and malaria disease progression. The proteins presented herein are mainly involved in inflammatory response, cellular adhesion and as constituents of erythrocyte membrane. These findings have a great potential to provide increased conceptual understanding of host-parasite interaction and malaria pathogenesis.

Atrophy of skin-draining lymph nodes predisposes for impaired immune responses to secondary infection in mice with chronic intestinal nematode infection.

Intestinal nematodes suppress immune responses in the context of allergy, gut inflammation, secondary infection and vaccination. Several mechanisms have been proposed for this suppression including alterations in Th2 cell differentiation and increased Treg cell suppressive function. In this study, we show that chronic nematode infection leads to reduced peripheral responses to vaccination because of a generalized reduction in the available responsive lymphocyte pool. We found that superficial skin-draining lymph nodes (LNs) in mice that are chronically infected with the intestinal nematode Heligmosomides polygyrus, do not reach the same cellularity as worm-free mice upon subsequent BCG infection in the skin. B cells and T cells, all declined in skin-draining LN of H. polygyrus-infected mice, resulting in LNs atrophy and altered lymphocyte composition. Importantly, anti-helminthic treatment improved lymphocyte numbers in skin-draining LN, indicating that time after de-worming is critical to regain full-scale LN cellularity. De-worming, and time for the skin LN to recover cellularity, also mended responses to Bacille Calmette-Guerin (BCG) in the LN draining the footpad injection site. Thus, our findings show that chronic nematode infection leads to a paucity of lymphocytes in peripheral lymph nodes, which acts to reduce the efficacy of immune responses at these sites.

Coordinated Changes in Gene Expression Throughout Encystation of Giardia intestinalis.

Differentiation into infectious cysts through the process of encystation is crucial for transmission and survival of the intestinal protozoan parasite Giardia intestinalis. Hitherto the majority of studies have focused on the early events, leaving late encystation poorly defined. In order to further study encystation, focusing on the later events, we developed a new encystation protocol that generates a higher yield of mature cysts compared to standard methods. Transcriptome changes during the entire differentiation from trophozoites to cysts were thereafter studied using RNA sequencing (RNA-seq). A high level of periodicity was observed for up- and down-regulated genes, both at the level of the entire transcriptome and putative regulators. This suggests the trajectory of differentiation to be coordinated through developmentally linked gene regulatory activities. Our study identifies a core of 13 genes that are consistently up-regulated during initial encystation. Of these, two constitute previously uncharacterized proteins that we were able to localize to a new type of encystation-specific vesicles. Interestingly, the largest transcriptional changes were seen in the late phase of encystation with the majority of the highly up-regulated genes encoding hypothetical proteins. Several of these were epitope-tagged and localized to further characterize these previously unknown genetic components of encystation and possibly excystation. Finally, we also detected a switch of variant specific surface proteins (VSPs) in the late phase of encystation. This occurred at the same time as nuclear division and DNA replication, suggesting a potential link between the processes.

Strand-specific RNA sequencing in Plasmodium falciparum malaria identifies developmentally regulated long non-coding RNA and circular RNA.

BACKGROUND: The human malaria parasite Plasmodium falciparum has a complex and multi-stage life cycle that requires extensive and precise gene regulation to allow invasion and hijacking of host cells, transmission, and immune escape. To date, the regulatory elements orchestrating these critical parasite processes remain largely unknown. Yet it is becoming increasingly clear that long non-coding RNAs (lncRNAs) could represent a missing regulatory layer across a broad range of organisms. RESULTS: To investigate the regulatory capacity of lncRNA in P. falciparum, we harvested fifteen samples from two time-courses. Our sample set profiled 56 h of P. falciparum blood stage development. We then developed and validated strand-specific, non-polyA-selected RNA sequencing methods, and pursued the first assembly of P. falciparum strand-specific transcript structures from RNA sequencing data. This approach enabled the annotation of over one thousand lncRNA transcript models and their comprehensive global analysis: coding prediction, periodicity, stage-specificity, correlation, GC content, length, location relative to annotated transcripts, and splicing. We validated the complete splicing structure of three lncRNAs with compelling properties. Non-polyA-selected deep sequencing also enabled the prediction of hundreds of intriguing P. falciparum circular RNAs, six of which we validated experimentally. CONCLUSIONS: We found that a subset of lncRNAs, including all subtelomeric lncRNAs, strongly peaked in expression during invasion. By contrast, antisense transcript levels significantly dropped during invasion. As compared to neighboring mRNAs, the expression of antisense-sense pairs was significantly anti-correlated during blood stage development, indicating transcriptional interference. We also validated that P. falciparum produces circRNAs, which is notable given the lack of RNA interference in the organism, and discovered that a highly expressed, five-exon antisense RNA is poised to regulate P. falciparum gametocyte development 1 (PfGDV1), a gene required for early sexual commitment events.

The cytoplasmic prolyl-tRNA synthetase of the malaria parasite is a dual-stage target of febrifugine and its analogs.

The emergence of drug resistance is a major limitation of current antimalarials. The discovery of new druggable targets and pathways including those that are critical for multiple life cycle stages of the malaria parasite is a major goal for developing next-generation antimalarial drugs. Using an integrated chemogenomics approach that combined drug resistance selection, whole-genome sequencing, and an orthogonal yeast model, we demonstrate that the cytoplasmic prolyl-tRNA (transfer RNA) synthetase (PfcPRS) of the malaria parasite Plasmodium falciparum is a biochemical and functional target of febrifugine and its synthetic derivative halofuginone. Febrifugine is the active principle of a traditional Chinese herbal remedy for malaria. We show that treatment with febrifugine derivatives activated the amino acid starvation response in both P. falciparum and a transgenic yeast strain expressing PfcPRS. We further demonstrate in the Plasmodium berghei mouse model of malaria that halofuginol, a new halofuginone analog that we developed, is active against both liver and asexual blood stages of the malaria parasite. Halofuginol, unlike halofuginone and febrifugine, is well tolerated at efficacious doses and represents a promising lead for the development of dual-stage next-generation antimalarials.

Immune characterization of Plasmodium falciparum parasites with a shared genetic signature in a region of decreasing transmission.

As the intensity of malaria transmission has declined, Plasmodium falciparum parasite populations have displayed decreased clonal diversity resulting from the emergence of many parasites with common genetic signatures (CGS). We have monitored such CGS parasite clusters from 2006 to 2013 in Thiès, Senegal, using the molecular barcode. The first, and one of the largest observed clusters of CGS parasites, was present in 24% of clinical isolates in 2008, declined to 3.4% of clinical isolates in 2009, and then disappeared. To begin to explore the relationship between the immune responses of the population and the emergence and decline of specific parasite genotypes, we have determined whether antibodies to CGS parasites correlate with their prevalence. We measured (i) antibodies capable of inhibiting parasite growth in culture and (ii) antibodies recognizing the surfaces of infected erythrocytes (RBCs). IgG obtained from volunteers in 2009 showed increased reactivity to the surfaces of CGS-parasitized erythrocytes over IgG from 2008. Since P. falciparum EMP-1 (PfEMP-1) is a major variant surface antigen, we used var Ups quantitative reverse transcription-PCR (qRT-PCR) and sequencing with degenerate DBL1α domain primers to characterize the var genes expressed by CGS parasites after short-term in vitro culture. CGS parasites show upregulation of UpsA var genes and 2-cysteine-containing PfEMP-1 molecules and express the same dominant var transcript. Our work indicates that the CGS parasites in this cluster express similar var genes, more than would be expected by chance in the population, and that there is year-to-year variation in immune recognition of surface antigens on CGS parasite-infected erythrocytes. This study lays the groundwork for detailed investigations of the mechanisms driving the expansion or contraction of specific parasite clones in the population.

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.

Improved in vitro culture of Plasmodium falciparum permits establishment of clinical isolates with preserved multiplication, invasion and rosetting phenotypes.

To be able to robustly propagate P. falciparum at optimal conditions in vitro is of fundamental importance for genotypic and phenotypic studies of both established and fresh clinical isolates. Cryo-preserved P. falciparum isolates from Ugandan children with severe or uncomplicated malaria were investigated for parasite phenotypes under different in vitro growth conditions or studied directly from the peripheral blood. The parasite cultures showed a minimal loss of parasite-mass and preserved percentage of multiple infected pRBCs to that in peripheral blood, maintained adhesive phenotypes and good outgrowth and multiplication rates when grown in suspension and supplemented with gas. In contrast, abnormal and greatly fluctuating levels of multiple infections were observed during static growth conditions and outgrowth and multiplication rates were inferior. Serum, as compared to Albumax, was found necessary for optimal presentation of PfEMP1 at the pRBC surface and/or for binding of serum proteins (immunoglobulins). Optimal in vitro growth conditions of P. falciparum therefore include orbital shaking (50 rev/min), human serum (10%) and a fixed gas composition (5% O2, 5% CO2, 90% N2). We subsequently established 100% of 76 frozen patient isolates and found rosetting with schizont pRBCs in every isolate (>26% schizont rosetting rate). Rosetting during schizogony was often followed by invasion of the bound RBC as seen by regular and time-lapse microscopy as well as transmission electron microscopy. The peripheral parasitemia, the level of rosetting and the rate of multiplication correlated positively to one another for individual isolates. Rosetting was also more frequent with trophozoite and schizont pRBCs of children with severe versus uncomplicated malaria (p<0.002; p<0.004). The associations suggest that rosetting enhances the ability of the parasite to multiply within the human host.

Acquired antibodies to merozoite antigens in children from Uganda with uncomplicated or severe Plasmodium falciparum malaria.

Malaria can present itself as an uncomplicated or severe disease. We have here studied the quantity and quality of antibody responses against merozoite antigens, as well as multiplicity of infection (MOI), in children from Uganda. We found higher levels of IgG antibodies toward erythrocyte-binding antigen EBA181, MSP2 of Plasmodium falciparum 3D7 and FC27 (MSP2-3D7/FC27), and apical membrane antigen 1 (AMA1) in patients with uncomplicated malaria by enzyme-linked immunosorbent assay (ELISA) but no differences against EBA140, EBA175, MSP1, and reticulocyte-binding protein homologues Rh2 and Rh4 or for IgM against MSP2-3D7/FC27.Patients with uncomplicated malaria were also shown to have higher antibody affinities for AMA1 by surface plasmon resonance (SPR). Decreased invasion of two clinical P. falciparum isolates in the presence of patient plasma correlated with lower initial parasitemia in the patients, in contrast to comparisons of parasitemia to ELISA values or antibody affinities, which did not show any correlations. Analysis of the heterogeneity of the infections revealed a higher MOI in patients with uncomplicated disease, with the P. falciparum K1 MSP1 (MSP1-K1) and MSP2-3D7 being the most discriminative allelic markers. Higher MOIs also correlated positively with higher antibody levels in several of the ELISAs. In conclusion, certain antibody responses and MOIs were associated with differences between uncomplicated and severe malaria. When different assays were combined, some antibodies, like those against AMA1, seemed particularly discriminative. However, only decreased invasion correlated with initial parasitemia in the patient, signaling the importance of functional assays in understanding development of immunity against malaria and in evaluating vaccine candidates.

Genetic diversity of Plasmodium falciparum infections in mild and severe malaria of children from Kampala, Uganda.

Diversity in parasite virulence is one of the factors that contribute to the clinical outcome of malaria infections. The association between the severity of Plasmodium falciparum malaria and the number of distinct parasite populations infecting the host (multiplicity of infection) or polymorphism within any of the specific antigen genes was investigated. The study included 164 children presenting with mild and severe malaria from central Uganda where malaria is meso-endemic. The polymorphic regions of the circumsporozoite protein (csp), merozoite surface proteins 1 and 2 (msp1 and msp2), and glutamate-rich protein (glurp) were genotyped by polymerase chain reaction methods and fragment analysis by gel electrophoresis. In a subset of samples fragment analysis was also performed by fluorescent PCR genotyping followed by capillary electrophoresis. The multiplicity of infection (MOI), determined as the highest number of alleles detected within any of the four genetic loci, was significantly higher in severe than in mild malaria cases (mean 3.7 and 3.0, respectively, P=0.002). No particular genotype or allelic family of msp1 or msp2 was associated with severity of malaria, and nor did the genotyping method reveal any significant difference in MOI when only assessed by msp2 genotyping. Severity of malaria was not linked to the predominance of any particular msp1 or msp2 allelic types, independent of methods used for genotyping. Monitoring the dynamics of multiple clone infections in relation to disease outcome, host immune status and genetic factors will provide more insight into parasite virulence mechanisms.

Nuclear repositioning precedes promoter accessibility and is linked to the switching frequency of a Plasmodium falciparum invasion gene.

Variation of surface adhesins, such as the Plasmodium falciparum erythrocyte invasion ligand PfRh4, is critical for virulence and immune evasion in many microbes. While phenotypic switching is linked to transcriptional changes and chromatin function, the determinants of switching frequency remain poorly defined. By expressing a prokaryotic DNA methylase in P. falciparum, we directly assayed accessibility of transcriptionally active and silent chromatin at the PfRh4 locus. Parasites selected for in vivo PfRh4 activation show a reversible increase in promoter accessibility and exhibit perinuclear repositioning of the locus from a silent to a conserved activation domain. Forced activation of a proximal gene results in a similar repositioning of the PfRh4 locus, with a concomitant increase in PfRh4 activation in a subpopulation of parasites and promoter accessibility correlating with actively transcribed loci. Thus, nuclear repositioning is associated with increased P. falciparum switching frequency, while promoter accessibility is tightly linked to clonally active PfRh4 promoters.

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.

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.

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.