Search

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed. Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

International audience; The last malaria parasite standing will display effective adaptations to selective forces. While substantial progress has been made in reducing malaria mortality, eradication will require elimination of all Plasmodium parasites, including those in asymptomatic infections. These typically chronic, low-density infections are difficult to detect, yet can persist for months. We argue that asymptomatic infection is the parasite’s best asset for survival but it can be exploited if studied as a new model for host–pathogen–vector interactions. Regular sampling from cohorts of asymptomatic individuals can provide a means to investigate continuous parasite development within its natural host. State-of-the-art techniques can now be applied to such infections. This approach may reveal key molecular drivers of chronic infections – a critical step for malaria eradication.

Background: Although thousands of clinical isolates of Plasmodium falciparum are being sequenced and analysed by short read technology, the data do not resolve the highly variable subtelomeric regions of the genomes that contain polymorphic gene families involved in immune evasion and pathogenesis. There is also no current standard definition of the boundaries of these variable subtelomeric regions. Methods: Using long-read sequence data (Pacific Biosciences SMRT technology), we assembled and annotated the genomes of 15 P. falciparum isolates, ten of which are newly cultured clinical isolates. We performed comparative analysis of the entire genome with particular emphasis on the subtelomeric regions and the internal var genes clusters. Results: The nearly complete sequence of these 15 isolates has enabled us to define a highly conserved core genome, to delineate the boundaries of the subtelomeric regions, and to compare these across isolates. We found highly structured variable regions in the genome. Some exported gene families purportedly involved in release of merozoites show copy number variation. As an example of ongoing genome evolution, we found a novel CLAG gene in six isolates. We also found a novel gene that was relatively enriched in the South East Asian isolates compared to those from Africa. Conclusions: These 15 manually curated new reference genome sequences with their nearly complete subtelomeric regions and fully assembled genes are an important new resource for the malaria research community. We report the overall conserved structure and pattern of important gene families and the more clearly defined subtelomeric regions.

Pathogen genome sequencing directly from clinical samples is quickly gaining importance in genetic and medical research studies. However, low DNA yield from blood-borne pathogens is often a limiting factor. The problem worsens in extremely base-biased genomes such as the AT-rich Plasmodium falciparum. We present a strategy for whole-genome amplification (WGA) of low-yield samples from P. falciparum prior to short-read sequencing. We have developed WGA conditions that incorporate tetramethylammonium chloride for improved amplification and coverage of AT-rich regions of the genome. We show that this method reduces amplification bias and chimera formation. Our data show that this method is suitable for as low as 10 pg input DNA, and offers the possibility of sequencing the parasite genome from small blood samples.

The most polymorphic gene family in P. falciparum is the ∼60 var genes distributed across parasite chromosomes, both in the subtelomeres and in internal regions. They encode hypervariable surface proteins known as P. falciparum erythrocyte membrane protein 1 (PfEMP1) that are critical for pathogenesis and immune evasion in Plasmodium falciparum. How var gene sequence diversity is generated is not currently completely understood. To address this, we constructed large clone trees and performed whole genome sequence analysis to study the generation of novel var gene sequences in asexually replicating parasites. While single nucleotide polymorphisms (SNPs) were scattered across the genome, structural variants (deletions, duplications, translocations) were focused in and around var genes, with considerable variation in frequency between strains. Analysis of more than 100 recombination events involving var exon 1 revealed that the average nucleotide sequence identity of two recombining exons was only 63% (range: 52.7–72.4%) yet the crossovers were error-free and occurred in such a way that the resulting sequence was in frame and domain architecture was preserved. Var exon 1, which encodes the immunologically exposed part of the protein, recombined in up to 0.2% of infected erythrocytes in vitro per life cycle. The high rate of var exon 1 recombination indicates that millions of new antigenic structures could potentially be generated each day in a single infected individual. We propose a model whereby var gene sequence polymorphism is mainly generated during the asexual part of the life cycle., Author Summary Malaria kills >600,000 people each year, with most deaths caused by Plasmodium falciparum. A family of proteins known as P. falciparum erythrocyte membrane protein 1, PfEMP1, is expressed on the surface of infected erythrocytes and plays an important role in pathogenesis. Each P. falciparum genome contains approximately 60 highly polymorphic var genes encoding the PfEMP1 proteins, and monoallelic expression with periodic switching results in immune evasion. Var gene polymorphism is thus critical to this survival strategy. We investigated how var gene diversity is generated by performing an in vitro evolution experiment, tracking var gene mutation in ‘real-time’ with whole genome sequencing. We found that genome structural variation is focused in and around var genes. These genetic rearrangements created new ‘chimeric’ var gene sequences during the mitotic part of the life cycle, and were consistent with processes of mitotic non-allelic homologous recombination. The recombinant var genes were always in frame and with conserved overall var gene architecture, and the recombination rate implies that many millions of rearranged var gene sequences are produced every 48-hour life cycle within infected individuals. In conclusion, we provide a detailed description of how new var gene sequences are continuously generated in the parasite genome, helping to explain long-term parasite survival within infected human hosts.

Sequence diversity in pathogen antigens is an obstacle to the development of interventions against many infectious diseases. In malaria caused by Plasmodium falciparum, the PfEMP1 family of variant surface antigens encoded by var genes are adhesion molecules that play a pivotal role in malaria pathogenesis and clinical disease. PfEMP1 is a major target of protective immunity, however, development of drugs or vaccines based on PfEMP1 is problematic due to extensive sequence diversity within the PfEMP1 family. Here we identified the PfEMP1 variants transcribed by P. falciparum strains selected for a virulence-associated adhesion phenotype (IgM-positive rosetting). The parasites transcribed a subset of Group A PfEMP1 variants characterised by an unusual PfEMP1 architecture and a distinct N-terminal domain (either DBLα1.5 or DBLα1.8 type). Antibodies raised in rabbits against the N-terminal domains showed functional activity (surface reactivity with live infected erythrocytes (IEs), rosette inhibition and induction of phagocytosis of IEs) down to low concentrations (, Author Summary Malaria remains one of the world's most deadly diseases. Life-threatening malaria is linked to a process called rosetting, in which malaria parasite-infected red blood cells bind to uninfected red cells to form aggregates that block blood flow in vital organs such as the brain. Current efforts to develop drugs or vaccines against rosetting are hindered by variation in the parasite rosette-mediating proteins, found on the surface of infected red cells. We studied these parasite-derived surface proteins and discovered that although they are variable, they share some common features. We raised antibodies against the rosette-mediating proteins, and found that they cross-reacted with multiple rosetting parasite strains from different countries around the world, including samples collected directly from African children with severe malaria. These findings provide new insights into malaria parasite interactions with human cells, and provide proof of principle that variable parasite molecules from virulent malaria parasites can induce strain-transcending antibodies. Hence, this work provides the foundation for the development of new therapies to treat or prevent life-threatening malaria.

Background The cytoadherence of Plasmodium falciparum is thought to be mediated by variant surface antigens (VSA), encoded by var, rif, stevor and pfmc-2tm genes. The last three families have rarely been studied in the context of cytoadherence. As most VSA genes are unique, the variability among sequences has impeded the functional study of VSA across different P. falciparum strains. However, many P. falciparum genomes have recently been sequenced, allowing the development of specific microarray probes for each VSA gene. Methods All VSA sequences from the HB3, Dd2 and IT/FCR3 genomes were extracted using HMMer software. Oligonucleotide probes were designed with OligoRankPick and added to the 3D7-based microarray chip. As a proof of concept, IT/R29 parasites were selected for and against rosette formation and the transcriptomes of isogenic rosetting and non-rosetting parasites were compared by microarray. Results From each parasite strain 50-56 var genes, 125-132 rif genes, 26-33 stevor genes and 3-8 pfmc-2tm genes were identified. Bioinformatic analysis of the new VSA sequences showed that 13 rif genes and five stevor genes were well-conserved across at least three strains (83-100% amino acid identity). The ability of the VSA-supplemented microarray chip to detect cytoadherence-related genes was assessed using P. falciparum clone IT/R29, in which rosetting is known to be mediated by PfEMP1 encoded by ITvar9. Whole transcriptome analysis showed that the most highly up-regulated gene in rosetting parasites was ITvar9 (19 to 429-fold up-regulated over six time points). Only one rif gene (IT4rifA_042) was up-regulated by more than four fold (five fold at 12 hours post-invasion), and no stevor or pfmc-2tm genes were up-regulated by more than two fold. 377 non-VSA genes were differentially expressed by three fold or more in rosetting parasites, although none was as markedly or consistently up-regulated as ITvar9. Conclusions Probes for the VSA of newly sequenced P. falciparum strains can be added to the 3D7-based microarray chip, allowing the analysis of the entire transcriptome of multiple strains. For the rosetting clone IT/R29, the striking transcriptional upregulation of ITvar9 was confirmed, and the data did not support the involvement of other VSA families in rosette formation.

Background Guanine-rich nucleic acid sequences are capable of folding into an intramolecular four-stranded structure called a G-quadruplex. When found in gene promoter regions, G-quadruplexes can downregulate gene expression, possibly by blocking the transcriptional machinery. Here we have used a genome-wide bioinformatic approach to identify Putative G-Quadruplex Sequences (PQS) in the Plasmodium falciparum genome, along with biophysical techniques to examine the physiological stability of P. falciparum PQS in vitro. Results We identified 63 PQS in the non-telomeric regions of the P. falciparum clone 3D7. Interestingly, 16 of these PQS occurred in the upstream region of a subset of the P. falciparum var genes (group B var genes). The var gene family encodes PfEMP1, the parasite's major variant antigen and adhesin expressed at the surface of infected erythrocytes, that plays a key role in malaria pathogenesis and immune evasion. The ability of the PQS found in the upstream regions of group B var genes (UpsB-Q) to form stable G-quadruplex structures in vitro was confirmed using 1H NMR, circular dichroism, UV spectroscopy, and thermal denaturation experiments. Moreover, the synthetic compound BOQ1 that shows a higher affinity for DNA forming quadruplex rather than duplex structures was found to bind with high affinity to the UpsB-Q. Conclusion This is the first demonstration of non-telomeric PQS in the genome of P. falciparum that form stable G-quadruplexes under physiological conditions in vitro. These results allow the generation of a novel hypothesis that the G-quadruplex sequences in the upstream regions of var genes have the potential to play a role in the transcriptional control of this major virulence-associated multi-gene family.

Severe malaria has a high mortality rate (15–20%) despite treatment with effective antimalarial drugs. Adjunctive therapies for severe malaria that target the underlying disease process are therefore urgently required. Adhesion of erythrocytes infected withPlasmodium falciparumto human cells has a key role in the pathogenesis of life-threatening malaria and could be targeted with antiadhesion therapy. Parasite adhesion interactions include binding to endothelial cells (cytoadherence), rosetting with uninfected erythrocytes and platelet-mediated clumping of infected erythrocytes. Recent research has started to define the molecular mechanisms of parasite adhesion, and antiadhesion therapies are being explored. However, many fundamental questions regarding the role of parasite adhesion in severe malaria remain unanswered. There is strong evidence that rosetting contributes to severe malaria in sub-Saharan Africa; however, the identity of other parasite adhesion phenotypes that are implicated in disease pathogenesis remains unclear. In addition, the possibility of geographic variation in adhesion phenotypes causing severe malaria, linked to differences in malaria transmission levels and host immunity, has been neglected. Further research is needed to realise the untapped potential of antiadhesion adjunctive therapies, which could revolutionise the treatment of severe malaria and reduce the high mortality rate of the disease.

International audience; Disrupted-In-Schizophrenia 1 (DISC1) is a risk factor for schizophrenia and other major mental illnesses. Its protein binding partners include the Nuclear Distribution Factor E Homologs (NDE1 and NDEL1), LIS1, and phosphodiesterases 4B and 4D (PDE4B and PDE4D). We demonstrate that NDE1, NDEL1 and LIS1, together with their binding partner dynein, associate with DISC1, PDE4B and PDE4D within the cell, and provide evidence that this complex is present at the centrosome. LIS1 and NDEL1 have been previously suggested to be synaptic, and we now demonstrate localisation of DISC1, NDE1, and PDE4B at syn-apses in cultured neurons. NDE1 is phosphorylated by cAMP-dependant Protein Kinase A (PKA), whose activity is, in turn, regulated by the cAMP hydrolysis activity of phosphodiesterases, including PDE4. We propose that DISC1 acts as an assembly scaffold for all of these proteins and that the NDE1/ NDEL1/LIS1/dynein complex is modulated by cAMP levels via PKA and PDE4.