Patterns of gene recombination shape var gene repertoires in Plasmodium falciparum: comparisons of geographically diverse isolates

A paradoxical population structure of var DBLα types in Africa

The var multigene family encodes the P. falciparum erythrocyte membrane protein 1 (PfEMP1), which is important in host-parasite interaction as a virulence factor and major surface antigen of the blood stages of the parasite, responsible for maintaining chronic infection. Whilst important in the biology of P. falciparum, these genes (50 to 60 genes per parasite genome) are routinely excluded from whole genome analyses due to their hyper-diversity, achieved primarily through recombination. The PfEMP1 head structure almost always consists of a DBLα-CIDR tandem. Categorised into different groups (upsA, upsB, upsC), different head structures have been associated with different ligand-binding affinities and disease severities. We study how conserved individual DBLα types are at the country, regional, and local scales in Sub-Saharan Africa. Using publicly-available sequence datasets and a novel ups classification algorithm, cUps, we performed an in silico exploration of DBLα conservation through time and space in Africa. In all three ups groups, the population structure of DBLα types in Africa consists of variants occurring at rare, low, moderate, and high frequencies. Non-rare variants were found to be temporally stable in a local area in endemic Ghana. When inspected across different geographical scales, we report different levels of conservation; while some DBLα types were consistently found in high frequencies in multiple African countries, others were conserved only locally, signifying local preservation of specific types. Underlying this population pattern is the composition of DBLα types within each isolate DBLα repertoire, revealed to also consist of a mix of types found at rare, low, moderate, and high frequencies in the population. We further discuss the adaptive forces and balancing selection, including host genetic factors, potentially shaping the evolution and diversity of DBLα types in Africa.

Measuring changes in Plasmodium falciparum census population size in response to sequential malaria control interventions

Here we introduce a new endpoint "census population size" to evaluate the epidemiology and control of Plasmodium falciparum infections, where the parasite, rather than the infected human host, is the unit of measurement. To calculate census population size, we rely on a definition of parasite variation known as multiplicity of infection (M O I v a r ), based on the hyper-diversity of the v a r multigene family. We present a Bayesian approach to estimate M O I v a r from sequencing and counting the number of unique DBLα tags (or DBLα types) of v a r genes, and derive from it census population size by summation of M O I v a r in the human population. We track changes in this parasite population size and structure through sequential malaria interventions by indoor residual spraying (IRS) and seasonal malaria chemoprevention (SMC) from 2012 to 2017 in an area of high-seasonal malaria transmission in northern Ghana. Following IRS, which reduced transmission intensity by > 90% and decreased parasite prevalence by ~40-50%, significant reductions in v a r diversity, M O I v a r , and population size were observed in ~2,000 humans across all ages. These changes, consistent with the loss of diverse parasite genomes, were short lived and 32-months after IRS was discontinued and SMC was introduced, v a r diversity and population size rebounded in all age groups except for the younger children (1-5 years) targeted by SMC. Despite major perturbations from IRS and SMC interventions, the parasite population remained very large and retained the v a r population genetic characteristics of a high-transmission system (high v a r diversity; low v a r repertoire similarity) demonstrating the resilience of P. falciparum to short-term interventions in high-burden countries of sub-Saharan Africa.

Decoding the impact of nuclear organization on antigenic variation in parasites

Antigenic variation as a strategy to evade the host adaptive immune response has evolved in divergent pathogens. Antigenic variation involves restricted, and often mutually exclusive, expression of dominant antigens and a periodic switch in antigen expression during infection. In eukaryotes, nuclear compartmentalization, including three-dimensional folding of the genome and physical separation of proteins in compartments or condensates, regulates mutually exclusive gene expression and chromosomal translocations. In this Review, we discuss the impact of nuclear organization on antigenic variation in the protozoan pathogens Trypanosoma brucei and Plasmodium falciparum. In particular, we highlight the relevance of nuclear organization in both mutually exclusive antigen expression and genome stability, which underlie antigenic variation.

Unravelling var complexity: Relationship between DBLα types and var genes in Plasmodium falciparum

The enormous diversity and complexity of var genes that diversify rapidly by recombination has led to the exclusion of assembly of these genes from major genome initiatives (e.g., Pf6). A scalable solution in epidemiological surveillance of var genes is to use a small ‘tag’ region encoding the immunogenic DBLα domain as a marker to estimate var diversity. As var genes diversify by recombination, it is not clear the extent to which the same tag can appear in multiple var genes. This relationship between marker and gene has not been investigated in natural populations. Analyses of in vitro recombination within and between var genes have suggested that this relationship would not be exclusive. Using a dataset of publicly-available assembled var sequences, we test this hypothesis by studying DBLα-var relationships for four study sites in four countries: Pursat (Cambodia) and Mae Sot (Thailand), representing low malaria transmission, and Navrongo (Ghana) and Chikwawa (Malawi), representing high malaria transmission. In all study sites, DBLα-var relationships were shown to be predominantly 1-to-1, followed by a second largest proportion of 1-to-2 DBLα-var relationships. This finding indicates that DBLα tags can be used to estimate not just DBLα diversity but var gene diversity when applied in a local endemic area. Epidemiological applications of this result are discussed.

Neutral vs. non-neutral genetic footprints of Plasmodium falciparum multiclonal infections

At a time when effective tools for monitoring malaria control and eradication efforts are crucial, the increasing availability of molecular data motivates their application to epidemiology. The multiplicity of infection (MOI), defined as the number of genetically distinct parasite strains co-infecting a host, is one key epidemiological parameter for evaluating malaria interventions. Estimating MOI remains a challenge for high-transmission settings where individuals typically carry multiple co-occurring infections. Several quantitative approaches have been developed to estimate MOI, including two cost-effective ones relying on molecular data: i) THE REAL McCOIL method is based on putatively neutral single nucleotide polymorphism loci, and ii) the varcoding method is a fingerprinting approach that relies on the diversity and limited repertoire overlap of the var multigene family encoding the major Plasmodium falciparum blood-stage antigen PfEMP1 and is therefore under selection. In this study, we assess the robustness of the MOI estimates generated with these two approaches by simulating P. falciparum malaria dynamics under three transmission conditions using an extension of a previously developed stochastic agent-based model. We demonstrate that these approaches are complementary and best considered across distinct transmission intensities. While varcoding can underestimate MOI, it allows robust estimation, especially under high transmission where repertoire overlap is extremely limited from frequency-dependent selection. In contrast, THE REAL McCOIL often considerably overestimates MOI, but still provides reasonable estimates for low and moderate transmission. Regardless of transmission intensity, results for THE REAL McCOIL indicate that an inaccurate tail at high MOI values is generated, and that at high transmission, an apparently reasonable estimated MOI distribution can arise from some degree of compensation between overestimation and underestimation. As many countries pursue malaria elimination targets, defining the most suitable approach to estimate MOI based on sample size and local transmission intensity is highly recommended for monitoring the impact of intervention programs.

CD36-A Host Receptor Necessary for Malaria Parasites to Establish and Maintain Infection

Plasmodium falciparum-infected erythrocytes (PfIEs) present P. falciparum erythrocyte membrane protein 1 proteins (PfEMP1s) on the cell surface, via which they cytoadhere to various endothelial cell receptors (ECRs) on the walls of human blood vessels. This prevents the parasite from passing through the spleen, which would lead to its elimination. Each P. falciparum isolate has about 60 different PfEMP1s acting as ligands, and at least 24 ECRs have been identified as interaction partners. Interestingly, in every parasite genome sequenced to date, at least 75% of the encoded PfEMP1s have a binding domain for the scavenger receptor CD36 widely distributed on host endothelial cells and many other cell types. Here, we discuss why the interaction between PfIEs and CD36 is optimal to maintain a finely regulated equilibrium that allows the parasite to multiply and spread while causing minimal harm to the host in most infections.

Bayesian estimation of community size and overlap from random subsamples

Counting the number of species, items, or genes that are shared between two groups, sets, or communities is a simple calculation when sampling is complete. However, when only partial samples are available, quantifying the overlap between two communities becomes an estimation problem. Furthermore, to calculate normalized measures of β-diversity, such as the Jaccard and Sorenson-Dice indices, one must also estimate the total sizes of the communities being compared. Previous efforts to address these problems have assumed knowledge of total community sizes and then used Bayesian methods to produce unbiased estimates with quantified uncertainty. Here, we address communities of unknown size and show that this produces systematically better estimates—both in terms of central estimates and quantification of uncertainty in those estimates. We further show how to use species, item, or gene count data to refine estimates of community size in a Bayesian joint model of community size and overlap.

Decoding the Host-Parasite Protein Interactions Involved in Cerebral Malaria Through Glares of Molecular Dynamics Simulations

Malaria causes millions of deaths every year. The malaria parasite spends a substantial part of its life cycle inside human erythrocytes. Inside erythrocytes, it synthesizes and displays various proteins onto the erythrocyte surface, such as Plasmodium falciparum erythrocytic membrane protein-1 (PfEMP1). This protein contains cysteine-rich interdomain region (CIDR) domains which have many subtypes based on sequence diversity and can cross-talk with host molecules. The CIDRα1.4 subtype can attach host endothelial protein C receptor (EPCR). This interaction facilitates infected erythrocyte adherence to brain endothelium and subsequent development of cerebral malaria. Through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method, we explored the mechanism of interaction in the CIDRα1-EPCR complex. We examined the structural behavior of two CIDRα1 molecules (encoded by HB3-isolate var03-gene and IT4-isolate var07-gene) with EPCR unbound and bound (complex) forms. ^HB3var03CIDRα1 in apo and complexed with EPCR was comparatively more stable than ^IT4var07CIDRα1. Both of the complexes adopted two distinct conformational energy states. The hydrophobic residues played a crucial role in the binding of both complexes. For ^HB3var03CIDRα1-EPCR, the dominant energetic components were total polar interactions, while in ^IT4var07CIDRα1-EPCR, the primary interaction was van der Waals and nonpolar solvation energy. The study also revealed details such as correlated conformational motions and secondary structure evolution. Further, it elucidated various hotspot residues involved in protein-protein recognition. Overall, our study provides additional information on the structural behavior of CIDR molecules in unbound and receptor-bound states, which will help to design potent inhibitors.

An accurate method for identifying recent recombinants from unaligned sequences

Motivation

Recombination is a fundamental process in molecular evolution, and the identification of recombinant sequences is thus of major interest. However, current methods for detecting recombinants are primarily designed for aligned sequences. Thus they struggle with analyses of highly diverse genes, such as the var genes of the malaria parasite Plasmodium falciparum, which are known to diversify primarily through recombination.

Results

We introduce an algorithm to detect recent recombinant sequences from a dataset without a full multiple alignment. Our algorithm can handle thousands of gene-length sequences without the need for a reference panel. We demonstrate the accuracy of our algorithm through extensive numerical simulations; in particular, it maintains its effectiveness in the presence of insertions and deletions. We apply our algorithm to a dataset of 17,335 DBLα types in var genes from Ghana, observing that sequences belonging to the same ups group or domain subclass recombine amongst themselves more frequently, and that non-recombinant DBLα types are more conserved than recombinant ones.

Availability

Source code is freely available at https://github.com/qianfeng2/detREC_program.

Supplementary information

Supplementary data are available at Bioinformatics online.

Analysis of var Gene Transcription Pattern Using DBLα Tags

The Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) antigens, which are encoded by a multigene family called var genes, are exported and inserted onto the surface of the infected erythrocytes. PfEMP1 plays a key role in the pathogenesis of severe malaria and are major targets of naturally acquired immunity. Studying the expression pattern of var genes in P. falciparum clinical isolates is crucial for understanding disease mechanism and immunity to malaria. However, var genes are highly variable, which makes it difficult to study their expression in clinical isolates obtained directly from malaria patients. In this chapter, we describe an approach for analysis of var gene expression that targets a region referred to as DBLα tag, which is relatively conserved in all var genes.

CRISPR-Cas9 Editing of the Plasmodium falciparum Genome: Special Applications

The virulence of Plasmodium falciparum has been attributed in large part to the expression on the surface of infected red blood cells of the variant surface antigen Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1). Different forms of this protein are encoded by individual members of the multicopy gene family called var. Two attributes of the var gene family are key to the pathogenesis of malaria caused by P. falciparum; the hyperrecombinogenic nature of the var gene family that continuously generates antigenic diversity within parasite populations, and the ability of parasites to express only a single var gene at a time and to switch which gene is expressed over the course of an infection. The unique attributes of CRISPR-Cas9 have been applied to help decipher the molecular mechanisms underlying these unusual properties of the var gene family, both as a source of the DNA double strand breaks that initiate var gene recombination and as a way to recruit molecular probes to specific regions of the genome. In this chapter, we describe these somewhat unusual applications of the CRISPR-Cas9 system.

Evolution of transcriptional control of antigenic variation and virulence in human and ape malaria parasites

Background

The most severe form of human malaria is caused by the protozoan parasite Plasmodium falciparum. This unicellular organism is a member of a subgenus of Plasmodium called the Laverania that infects apes, with P. falciparum being the only member that infects humans. The exceptional virulence of this species to humans can be largely attributed to a family of variant surface antigens placed by the parasites onto the surface of infected red blood cells that mediate adherence to the vascular endothelium. These proteins are encoded by a large, multicopy gene family called var, with each var gene encoding a different form of the protein. By changing which var gene is expressed, parasites avoid immune recognition, a process called antigenic variation that underlies the chronic nature of malaria infections.

Results

Here we show that the common ancestor of the branch of the Laverania lineage that includes the human parasite underwent a remarkable change in the organization and structure of elements linked to the complex transcriptional regulation displayed by the var gene family. Unlike the other members of the Laverania, the clade that gave rise to P. falciparum evolved distinct subsets of var genes distinguishable by different upstream transcriptional regulatory regions that have been associated with different expression profiles and virulence properties. In addition, two uniquely conserved var genes that have been proposed to play a role in coordinating transcriptional switching similarly arose uniquely within this clade. We hypothesize that these changes originated at a time of dramatic climatic change on the African continent that is predicted to have led to significant changes in transmission dynamics, thus selecting for patterns of antigenic variation that enabled lengthier, more chronic infections.

Conclusions

These observations suggest that changes in transmission dynamics selected for significant alterations in the transcriptional regulatory mechanisms that mediate antigenic variation in the parasite lineage that includes P. falciparum. These changes likely underlie the chronic nature of these infections as well as their exceptional virulence.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01872-z.

A suitable RNA preparation methodology for whole transcriptome shotgun sequencing harvested from Plasmodium vivax-infected patients

Plasmodium vivax is a world-threatening human malaria parasite, whose biology remains elusive. The unavailability of in vitro culture, and the difficulties in getting a high number of pure parasites makes RNA isolation in quantity and quality a challenge. Here, a methodological outline for RNA-seq from P. vivax isolates with low parasitemia is presented, combining parasite maturation and enrichment with efficient RNA extraction, yielding ~ 100 pg.µL⁻¹ of RNA, suitable for SMART-Seq Ultra-Low Input RNA library and Illumina sequencing. Unbiased coding transcriptome of ~ 4 M reads was achieved for four patient isolates with ~ 51% of transcripts mapped to the P. vivax P01 reference genome, presenting heterogeneous profiles of expression among individual isolates. Amongst the most transcribed genes in all isolates, a parasite-staged mixed repertoire of conserved parasite metabolic, membrane and exported proteins was observed. Still, a quarter of transcribed genes remain functionally uncharacterized. In parallel, a P. falciparum Brazilian isolate was also analyzed and 57% of its transcripts mapped against IT genome. Comparison of transcriptomes of the two species revealed a common trophozoite-staged expression profile, with several homologous genes being expressed. Collectively, these results will positively impact vivax research improving knowledge of P. vivax biology.

Febrile temperature causes transcriptional downregulation of Plasmodium falciparum Sirtuins through Hsp90-dependent epigenetic modification

Sirtuins (PfSIR2A and PfSIR2B) are implicated to play pivotal roles in the silencing of sub-telomeric genes and the maintenance of telomere length in P. falciparum 3D7 strain. Here, we identify the key factors that regulate the cellular abundance and activity of these two histone deacetylases. Our results demonstrate that PfSIR2A and PfSIR2B are transcriptionally downregulated at the mid-ring stage in response to febrile temperature. We found that the molecular chaperone PfHsp90 acts as a repressor of PfSIR2A & B transcription. By virtue of its presence in the PfSIR2A & B promoter proximal regions PfHsp90 helps recruiting H3K9me3, conferring heterochromatic state, and thereby leading to the downregulation of PfSIR2A & B transcription. Such transcriptional downregulation can be reversed by the addition of 17-(allylamino)-17-demethoxygeldanamycin or Radicicol, two potent inhibitors of PfHsp90. The reduced occupancy of PfSir2 at sub-telomeric var promoters leads to the de-repression of var genes. Thus, here we uncover how exposure to febrile temperature, a hallmark of malaria, enables the parasites to manipulate the expression of the two prominent epigenetic modifiers PfSir2A and PfSir2B.

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.

Sequence analysis of wheat subtelomeres reveals a high polymorphism among homoeologous chromosomes

Bread wheat, Triticum aestivum L., is one of the most important crops in the world. Understanding its genome organization (allohexaploid; AABBDD; 2n = 6x = 42) is essential for geneticists and plant breeders. Particularly, the knowledge of how homologous chromosomes (equivalent chromosomes from the same genome) specifically recognize each other to pair at the beginning of meiosis, the cellular process to generate gametes in sexually reproducing organisms, is fundamental for plant breeding and has a big influence on the fertility of wheat plants. Initial homologous chromosome interactions contribute to specific recognition and pairing between homologues at the onset of meiosis. Understanding the molecular basis of these critical processes can help to develop genetic tools in a breeding context to promote interspecific chromosome associations in hybrids or interspecific genetic crosses to facilitate the transfer of desirable agronomic traits from related species into a crop like wheat. The terminal regions of chromosomes, which include telomeres and subtelomeres, participate in chromosome recognition and pairing. We present a detailed molecular analysis of subtelomeres of wheat chromosome arms 1AS, 4AS, 7AS, 7BS and 7DS. Results showed a high polymorphism in the subtelomeric region among homoeologues (equivalent chromosomes from related genomes) for all the features analyzed, including genes, transposable elements, repeats, GC content, predicted CpG islands, recombination hotspots and targeted sequence motifs for relevant DNA-binding proteins. These polymorphisms might be the molecular basis for the specificity of homologous recognition and pairing in initial chromosome interactions at the beginning of meiosis in wheat.

Genome wide distribution of G-quadruplexes and their impact on gene expression in malaria parasites

Mechanisms of transcriptional control in malaria parasites are still not fully understood. The positioning patterns of G-quadruplex (G4) DNA motifs in the parasite's AT-rich genome, especially within the var gene family which encodes virulence factors, and in the vicinity of recombination hotspots, points towards a possible regulatory role of G4 in gene expression and genome stability. Here, we carried out the most comprehensive genome-wide survey, to date, of G4s in the Plasmodium falciparum genome using G4Hunter, which identifies G4 forming sequences (G4FS) considering their G-richness and G-skewness. We show an enrichment of G4FS in nucleosome-depleted regions and in the first exon of var genes, a pattern that is conserved within the closely related Laverania Plasmodium parasites. Under G4-stabilizing conditions, i.e., following treatment with pyridostatin (a high affinity G4 ligand), we show that a bona fide G4 found in the non-coding strand of var promoters modulates reporter gene expression. Furthermore, transcriptional profiling of pyridostatin-treated parasites, shows large scale perturbations, with deregulation affecting for instance the ApiAP2 family of transcription factors and genes involved in ribosome biogenesis. Overall, our study highlights G4s as important DNA secondary structures with a role in Plasmodium gene expression regulation, sub-telomeric recombination and var gene biology.

Evolution of Host Specificity by Malaria Parasites through Altered Mechanisms Controlling Genome Maintenance

Malaria remains one of the most prevalent and deadly infectious diseases of the developing world, causing approximately 228 million clinical cases and nearly half a million deaths annually. The disease is caused by protozoan parasites of the genus Plasmodium, and of the five species capable of infecting humans, infections with P. falciparum are the most severe. In addition to the parasites that infect people, there are hundreds of additional species that infect birds, reptiles, and other mammals, each exquisitely evolved to meet the specific challenges inherent to survival within their respective hosts. By comparing the unique strategies that each species has evolved, key insights into host-parasite interactions can be gained, including discoveries regarding the pathogenesis of human disease. Here, we describe the surprising observation that closely related parasites with different hosts have evolved remarkably different methods for repairing their genomes. This observation has important implications for the ability of parasites to maintain chronic infections and for the development of host immunity.

The protozoan parasites that cause malaria infect a wide variety of vertebrate hosts, including birds, reptiles, and mammals, and the evolutionary pressures inherent to the host-parasite relationship have profoundly shaped the genomes of both host and parasite. Here, we report that these selective pressures have resulted in unexpected alterations to one of the most basic aspects of eukaryotic biology, the maintenance of genome integrity through DNA repair. Malaria parasites that infect humans continuously generate genetic diversity within their antigen-encoding gene families through frequent ectopic recombination between gene family members, a process that is a crucial feature of the persistence of malaria globally. The continuous generation of antigen diversity ensures that different parasite isolates are antigenically distinct, thus preventing extensive cross-reactive immunity and enabling parasites to maintain stable transmission within human populations. However, the molecular basis of the recombination between gene family members is not well understood. Through computational analyses of the antigen-encoding, multicopy gene families of different Plasmodium species, we report the unexpected observation that malaria parasites that infect rodents do not display the same degree of antigen diversity as observed in Plasmodium falciparum and appear to undergo significantly less ectopic recombination. Using comparative genomics, we also identify key molecular components of the diversification process, thus shedding new light on how malaria parasites balance the maintenance of genome integrity with the requirement for continuous genetic diversification.

Protective Immunity against Severe Malaria in Children Is Associated with a Limited Repertoire of Antibodies to Conserved PfEMP1 Variants

Extreme diversity of the major Plasmodium falciparum antigen, PfEMP1, poses a barrier to identifying targets of immunity to malaria. Here, we used protein microarrays containing hundreds of variants of the DBLα domain of PfEMP1 to cover the diversity of Papua New Guinean (PNG) parasites. Probing the plasma of a longitudinal cohort of malaria-exposed PNG children showed that group 2 DBLα antibodies were moderately associated with a lower risk of uncomplicated malaria, whereas individual variants were only weakly associated with clinical immunity. In contrast, antibodies to 85 individual group 1 and 2 DBLα variants were associated with a 70%-100% reduction in severe malaria. Of these, 17 variants were strong predictors of severe malaria. Analysis of full-length PfEMP1 sequences from PNG samples shows that these 17 variants are linked to pathogenic CIDR domains. This suggests that whereas immunity to uncomplicated malaria requires a broad repertoire of antibodies, immunity to severe malaria targets a subset of conserved variants. These findings provide insights into antimalarial immunity and potential antibody biomarkers for disease risk.

Competition for hosts modulates vast antigenic diversity to generate persistent strain structure in Plasmodium falciparum

In their competition for hosts, parasites with antigens that are novel to the host immune system will be at a competitive advantage. The resulting frequency-dependent selection can structure parasite populations into strains of limited genetic overlap. For the causative agent of malaria, Plasmodium falciparum, the high recombination rates and associated vast diversity of its highly antigenic and multicopy var genes preclude such clear clustering in endemic regions. This undermines the definition of strains as specific, temporally persisting gene variant combinations. We use temporal multilayer networks to analyze the genetic similarity of parasites in both simulated data and in an extensively and longitudinally sampled population in Ghana. When viewed over time, populations are structured into modules (i.e., groups) of parasite genomes whose var gene combinations are more similar within than between the modules and whose persistence is much longer than that of the individual genomes that compose them. Comparison to neutral models that retain parasite population dynamics but lack competition reveals that the selection imposed by host immunity promotes the persistence of these modules. The modular structure is, in turn, associated with a slower acquisition of immunity by individual hosts. Modules thus represent dynamically generated niches in host immune space, which can be interpreted as strains. Negative frequency-dependent selection therefore shapes the organization of the var diversity into parasite genomes, leaving a persistence signature over ecological time scales. Multilayer networks extend the scope of phylodynamics analyses by allowing quantification of temporal genetic structure in organisms that generate variation via recombination or other non-bifurcating processes. A strain structure similar to the one described here should apply to other pathogens with large antigenic spaces that evolve via recombination. For malaria, the temporal modular structure should enable the formulation of tractable epidemiological models that account for parasite antigenic diversity and its influence on intervention outcomes.

A combination of computational modeling and empirical data reveals persistent strain structure despite vast antigenic diversity in the human malaria parasite Plasmodium falciparum, with potential consequences for the acquisition of immunity.

Rapid antigen diversification through mitotic recombination in the human malaria parasite Plasmodium falciparum

Malaria parasites possess the remarkable ability to maintain chronic infections that fail to elicit a protective immune response, characteristics that have stymied vaccine development and cause people living in endemic regions to remain at risk of malaria despite previous exposure to the disease. These traits stem from the tremendous antigenic diversity displayed by parasites circulating in the field. For Plasmodium falciparum, the most virulent of the human malaria parasites, this diversity is exemplified by the variant gene family called var, which encodes the major surface antigen displayed on infected red blood cells (RBCs). This gene family exhibits virtually limitless diversity when var gene repertoires from different parasite isolates are compared. Previous studies indicated that this remarkable genome plasticity results from extensive ectopic recombination between var genes during mitotic replication; however, the molecular mechanisms that direct this process to antigen-encoding loci while the rest of the genome remains relatively stable were not determined. Using targeted DNA double-strand breaks (DSBs) and long-read whole-genome sequencing, we show that a single break within an antigen-encoding region of the genome can result in a cascade of recombination events leading to the generation of multiple chimeric var genes, a process that can greatly accelerate the generation of diversity within this family. We also found that recombinations did not occur randomly, but rather high-probability, specific recombination products were observed repeatedly. These results provide a molecular basis for previously described structured rearrangements that drive diversification of this highly polymorphic gene family.

Comprehensive analysis of Fc-mediated IgM binding to the Plasmodium falciparum erythrocyte membrane protein 1 family in three parasite clones

PfEMP1 is a family of adhesive proteins expressed on the surface of Plasmodium falciparum-infected erythrocytes (IEs), where they mediate adhesion of IEs to a range of host receptors. Efficient PfEMP1-dependent IE sequestration often depends on soluble serum proteins, including IgM. Here, we report a comprehensive investigation of which of the about 60 var gene-encoded PfEMP1 variants per parasite genome can bind IgM via the Fc part of the antibody molecule, and which of the constituent domains of those PfEMP1 are involved. We erased the epigenetic memory of var gene expression in three distinct P. falciparum clones, 3D7, HB3, and IT4/FCR3 by promoter titration, and then isolated individual IEs binding IgM from malaria-unexposed individuals by fluorescence-activated single-cell sorting. The var gene transcription profiles of sub-clones measured by real-time qPCR were used to identify potential IgM-binding PfEMP1 variants. Recombinant DBL and CIDR domains corresponding to those variants were tested by ELISA and protein arrays to confirm their IgM-binding capacity. Selected DBL domains were used to raise specific rat anti-sera to select IEs with uniform expression of candidate PfEMP1 proteins. Our data document that IgM-binding PfEMP1 proteins are common in each of the three clones studied, and that the binding epitopes are mainly found in DBLε and DBLζ domains near the C-terminus.

Matched Placental and Circulating Plasmodium falciparum Parasites are Genetically Homologous at the var2csa ID1-DBL2X Locus by Deep Sequencing

In pregnancy-associated malaria, infected erythrocytes accumulate in the placenta. It is unclear if in polyclonal infections this results in distinct peripheral and placental parasite populations. We used long amplicon deep sequencing of Plasmodium falciparum var2csa ID1-DBL2X from 15 matched peripheral and placental samples collected at delivery from a high transmission area to determine genetic homology. Despite substantial sequence variation and detecting 23 haplotypes, the matched pairs mostly contained the same genetic variants, with 11 pairs sharing 100% of their variants, whereas others showed some heterogeneity. Thus, at delivery, peripheral and placental parasites appear to intermix and placental genotypes can be inferred through peripheral sampling.

Blood outgrowth endothelial cells (BOECs) as a novel tool for studying adhesion of Plasmodium falciparum-infected erythrocytes

The lack of suitable animal models for the study of cytoadhesion of P. falciparum-infected erythrocytes (IEs) has necessitated in vitro studies employing a range of cell lines of either human tumour origin (e.g., BeWo and C32 cells) or non-human origin (e.g., CHO cells). Of the human cells available, many were isolated from adults, or derived from a pool of donors (e.g., HBEC-5i). Here we demonstrate, for the first time, the successful isolation of blood outgrowth endothelial cells (BOECs) from frozen stabilates of peripheral blood mononuclear cells obtained from small-volume peripheral blood samples from paediatric malaria patients. BOECs are a sub-population of human endothelial cells, found within the peripheral blood. We demonstrate that these cells express receptors such as Intercellular Adhesion Molecule 1 (ICAM-1/CD54), Endothelial Protein C Receptor (EPCR/CD201), platelet/endothelial cell adhesion molecule 1 (PECAM-1/CD31), Thrombomodulin (CD141), and support adhesion of P. falciparum IEs.

Multiscale Immune Selection and the Transmission-Diversity Feedback in Antigenically Diverse Pathogen Systems

Antigenic diversity is commonly used by pathogens to enhance their transmission success. Within-host clonal antigenic variation helps to maintain long infectious periods, whereas high levels of allelic diversity at the population level significantly expand the pool of susceptible individuals. Diversity, however, is not necessarily a static property of a pathogen population but in many cases is generated by the very act of infection and transmission, and it is therefore expected to respond dynamically to changes in transmission and immune selection. We hypothesized that this coupling creates a positive feedback whereby infection and disease transmission promote the generation of diversity, which itself facilitates immune evasion and further infections. To investigate this link in more detail, we considered the human malaria parasite Plasmodium falciparum, one of the most important antigenically diverse pathogens. We developed an individual-based model in which antigenic diversity emerges as a dynamic property from the underlying transmission processes. Our results show that the balance between stochastic extinction and the generation of new antigenic variants is intrinsically linked to within-host and between-host immune selection. This in turn determines the level of diversity that can be maintained in a given population. Furthermore, the transmission-diversity feedback can lead to temporal lags in the response to natural or intervention-induced perturbations in transmission rates. Our results therefore have important implications for monitoring and assessing the effectiveness of disease control efforts.

RecQ helicases in the malaria parasite Plasmodium falciparum affect genome stability, gene expression patterns and DNA replication dynamics

The malaria parasite Plasmodium falciparum has evolved an unusual genome structure. The majority of the genome is relatively stable, with mutation rates similar to most eukaryotic species. However, some regions are very unstable with high recombination rates, driving the generation of new immune evasion-associated var genes. The molecular factors controlling the inconsistent stability of this genome are not known. Here we studied the roles of the two putative RecQ helicases in P. falciparum, PfBLM and PfWRN. When PfWRN was knocked down, recombination rates increased four-fold, generating chromosomal abnormalities, a high rate of chimeric var genes and many microindels, particularly in known ‘fragile sites’. This is the first identification of a gene involved in suppressing recombination and maintaining genome stability in Plasmodium. By contrast, no change in mutation rate appeared when the second RecQ helicase, PfBLM, was mutated. At the transcriptional level, however, both helicases evidently modulate the transcription of large cohorts of genes, with several hundred genes—including a large proportion of vars—showing deregulated expression in each RecQ mutant. Aberrant processing of stalled replication forks is a possible mechanism underlying elevated mutation rates and this was assessed by measuring DNA replication dynamics in the RecQ mutant lines. Replication forks moved slowly and stalled at elevated rates in both mutants, confirming that RecQ helicases are required for efficient DNA replication. Overall, this work identifies the Plasmodium RecQ helicases as major players in DNA replication, antigenic diversification and genome stability in the most lethal human malaria parasite, with important implications for genome evolution in this pathogen.

Human malaria is caused by Plasmodium parasites, with most of the mortality (almost half a million deaths each year) being caused by one species, Plasmodium falciparum. This parasite has an unusual genome: it is exceptionally biased towards A and T nucleotides rather than G and C, and it contains specific areas rich in hypervariable virulence-associated genes which evolve very rapidly to produce new variants. This evolution is probably vital for the parasite to evade the human immune system and maintain chronic infections, but how it is controlled at a molecular level remains unknown. We have identified a helicase in the parasite with a huge influence on genome stability and the rate of genome evolution. It appears to function by unwinding various unusual DNA structures, and if this fails then the genome becomes unstable. In addition, the transcription of many genes whose DNA tends to form secondary structures is affected, and DNA replication is impeded. If this helicase was expressed variably in different parasite strains infecting humans, it could influence the parasites’ ability to grow and replicate efficiently, and also, crucially, its ability to evolve and thus evade the human immune system.

Networks of genetic similarity reveal non-neutral processes shape strain structure in Plasmodium falciparum

Pathogens compete for hosts through patterns of cross-protection conferred by immune responses to antigens. In Plasmodium falciparum malaria, the var multigene family encoding for the major blood-stage antigen PfEMP1 has evolved enormous genetic diversity through ectopic recombination and mutation. With 50-60 var genes per genome, it is unclear whether immune selection can act as a dominant force in structuring var repertoires of local populations. The combinatorial complexity of the var system remains beyond the reach of existing strain theory and previous evidence for non-random structure cannot demonstrate immune selection without comparison with neutral models. We develop two neutral models that encompass malaria epidemiology but exclude competitive interactions between parasites. These models, combined with networks of genetic similarity, reveal non-neutral strain structure in both simulated systems and an extensively sampled population in Ghana. The unique population structure we identify underlies the large transmission reservoir characteristic of highly endemic regions in Africa.

Rethinking the extrinsic incubation period of malaria parasites

The time it takes for malaria parasites to develop within a mosquito, and become transmissible, is known as the extrinsic incubation period, or EIP. EIP is a key parameter influencing transmission intensity as it combines with mosquito mortality rate and competence to determine the number of mosquitoes that ultimately become infectious. In spite of its epidemiological significance, data on EIP are scant. Current approaches to estimate EIP are largely based on temperature-dependent models developed from data collected on parasite development within a single mosquito species in the 1930s. These models assume that the only factor affecting EIP is mean environmental temperature. Here, we review evidence to suggest that in addition to mean temperature, EIP is likely influenced by genetic diversity of the vector, diversity of the parasite, and variation in a range of biotic and abiotic factors that affect mosquito condition. We further demonstrate that the classic approach of measuring EIP as the time at which mosquitoes first become infectious likely misrepresents EIP for a mosquito population. We argue for a better understanding of EIP to improve models of transmission, refine predictions of the possible impacts of climate change, and determine the potential evolutionary responses of malaria parasites to current and future mosquito control tools.

The Plasmodium falciparum transcriptome in severe malaria reveals altered expression of genes involved in important processes including surface antigen-encoding var genes

Within the human host, the malaria parasite Plasmodium falciparum is exposed to multiple selection pressures. The host environment changes dramatically in severe malaria, but the extent to which the parasite responds to-or is selected by-this environment remains unclear. From previous studies, the parasites that cause severe malaria appear to increase expression of a restricted but poorly defined subset of the PfEMP1 variant, surface antigens. PfEMP1s are major targets of protective immunity. Here, we used RNA sequencing (RNAseq) to analyse gene expression in 44 parasite isolates that caused severe and uncomplicated malaria in Papuan patients. The transcriptomes of 19 parasite isolates associated with severe malaria indicated that these parasites had decreased glycolysis without activation of compensatory pathways; altered chromatin structure and probably transcriptional regulation through decreased histone methylation; reduced surface expression of PfEMP1; and down-regulated expression of multiple chaperone proteins. Our RNAseq also identified novel associations between disease severity and PfEMP1 transcripts, domains, and smaller sequence segments and also confirmed all previously reported associations between expressed PfEMP1 sequences and severe disease. These findings will inform efforts to identify vaccine targets for severe malaria and also indicate how parasites adapt to-or are selected by-the host environment in severe malaria.

Signatures of competition and strain structure within the major blood-stage antigen of Plasmodium falciparum in a local community in Ghana

The concept of niche partitioning has received considerable theoretical attention at the interface of ecology and evolution of infectious diseases. Strain theory postulates that pathogen populations can be structured into distinct nonoverlapping strains by frequency-dependent selection in response to intraspecific competition for host immune space. The malaria parasite Plasmodium falciparum presents an opportunity to investigate this phenomenon in nature, under conditions of high recombination rate and extensive antigenic diversity. The parasite's major blood-stage antigen, Pf EMP1, is encoded by the hyperdiverse var genes. With a dataset that includes thousands of var DBLα sequence types sampled from asymptomatic cases within an area of high endemicity in Ghana, we address how var diversity is distributed within isolates and compare this to the distribution of microsatellite allelic diversity within isolates to test whether antigenic and neutral regions of the genome are structured differently. With respect to var DBLα sequence types, we find that on average isolates exhibit significantly lower overlap than expected randomly, but that there also exists frequent pairs of isolates that are highly related. Furthermore, the linkage network of var DBLα sequence types reveals a pattern of nonrandom modularity unique to these antigenic genes, and we find that modules of highly linked DBLα types are not explainable by neutral forces related to var recombination constraints, microsatellite diversity, sampling location, host age, or multiplicity of infection. These findings of reduced overlap and modularity among the var antigenic genes are consistent with a role for immune selection as proposed by strain theory. Identifying the evolutionary and ecological dynamics that are responsible for the nonrandom structure in P. falciparum antigenic diversity is important for designing effective intervention in endemic areas.

In silico guided reconstruction and analysis of ICAM-1-binding var genes from Plasmodium falciparum

The Plasmodium falciparum variant surface antigen PfEMP1 expressed on the surface of infected erythrocytes is thought to play a major role in the pathology of severe malaria. As the sequence pool of the var genes encoding PfEMP1 expands there are opportunities, despite the high degree of sequence diversity demonstrated by this gene family, to reconstruct full-length var genes from small sequence tags generated from patient isolates. To test whether this is possible we have used a set of recently laboratory adapted ICAM-1-binding parasite isolates to generate sequence tags and, from these, to identify the full-length PfEMP1 being expressed by them. In a subset of the strains available we were able to produce validated, full-length var gene sequences and use these to conduct biophysical analyses of the ICAM-1 binding regions.

A re-assessment of gene-tag classification approaches for describing var gene expression patterns during human Plasmodium falciparum malaria parasite infections

PfEMP1 are variant parasite antigens that are inserted on the surface of Plasmodium falciparum infected erythrocytes (IE). Through interactions with various host molecules, PfEMP1 mediate IE sequestration in tissues and play a key role in the pathology of severe malaria. PfEMP1 is encoded by a diverse multi-gene family called var. Previous studies have shown that that expression of specific subsets of var genes are associated with low levels of host immunity and severe malaria. However, in most clinical studies to date, full-length var gene sequences were unavailable and various approaches have been used to make comparisons between var gene expression profiles in different parasite isolates using limited information. Several studies have relied on the classification of a 300 - 500 base-pair "DBLα tag" region in the DBLα domain located at the 5' end of most var genes. We assessed the relationship between various DBLα tag classification methods, and sequence features that are only fully assessable through full-length var gene sequences. We compared these different sequence features in full-length var gene from six fully sequenced laboratory isolates. These comparisons show that despite a long history of recombination,  DBLα sequence tag classification can provide functional information on important features of full-length var genes. Notably, a specific subset of DBLα tags previously defined as "group A-like" is associated with CIDRα1 domains proposed to bind to endothelial protein C receptor. This analysis helps to bring together different sources of data that have been used to assess var gene expression in clinical parasite isolates.

Population genomics of virulence genes of Plasmodium falciparum in clinical isolates from Uganda

Plasmodium falciparum causes a spectrum of malarial disease from asymptomatic to uncomplicated through to severe. Investigations of parasite virulence have associated the expression of distinct variants of the major surface antigen of the blood stages known as Pf EMP1 encoded by up to 60 var genes per genome. Looking at the population genomics of var genes in cases of uncomplicated malaria, we set out to determine if there was any evidence of a selective sweep of specific var genes or clonal epidemic structure related to the incidence of uncomplicated disease in children. By sequencing the conserved DBLα domain of var genes from six sentinel sites in Uganda we found that the parasites causing uncomplicated P. falciparum disease in children were highly diverse and that every child had a unique var DBLα repertoire. Despite extensive var DBLα diversity and minimal overlap between repertoires, specific DBLα types and groups were conserved at the population level across Uganda. This pattern was the same regardless of the geographic distance or malaria transmission intensity. These data lead us to propose that any parasite can cause uncomplicated malarial disease and that these diverse parasite repertoires are composed of both upsA and non-upsA var gene groups.

Characterisation of Plasmodium falciparum populations selected on the human endothelial receptors P-selectin, E-selectin, CD9 and CD151

The ability of the parasite Plasmodium falciparum to evade the immune system and be sequestered within human small blood vessels is responsible for severe forms of malaria. The sequestration depends on the interaction between human endothelial receptors and P. falciparum erythrocyte membrane protein 1 (PfEMP1) exposed on the surface of the infected erythrocytes (IEs). In this study, the transcriptomes of parasite populations enriched for parasites that bind to human P-selectin, E-selectin, CD9 and CD151 receptors were analysed. IT4_var02 and IT4_var07 were specifically expressed in IT4 parasite populations enriched for P-selectin-binding parasites; eight var genes (IT4_var02/07/09/13/17/41/44/64) were specifically expressed in isolate populations enriched for CD9-binding parasites. Interestingly, IT4 parasite populations enriched for E-selectin- and CD151-binding parasites showed identical expression profiles to those of a parasite population exposed to wild-type CHO-745 cells. The same phenomenon was observed for the 3D7 isolate population enriched for binding to P-selectin, E-selectin, CD9 and CD151. This implies that the corresponding ligands for these receptors have either weak binding capacity or do not exist on the IE surface. Conclusively, this work expanded our understanding of P. falciparum adhesive interactions, through the identification of var transcripts that are enriched within the selected parasite populations.

Pre-clinical and clinical development of the first placental malaria vaccine

INTRODUCTION

Malaria during pregnancy is a massive health problem in endemic areas. Placental malaria infections caused by Plasmodium falciparum are responsible for up to one million babies being born with a low birth weight every year. Significant efforts have been invested into preventing the condition. Areas covered: Pub Med was searched using the broad terms 'malaria parasite placenta' to identify studies of interactions between parasite and host, 'prevention of placental malaria' to identify current strategies to prevent placental malaria, and 'placental malaria vaccine' to identify pre-clinical vaccine development. However, all papers from these searches were not systematically included. Expert commentary: The first phase I clinical trials of vaccines are well underway. Trials testing efficacy are more complicated to carry out as only women that are exposed to parasites during pregnancy will contribute to endpoint measurements, further it may require extensive follow-up to establish protection. Future second generation vaccines may overcome the inherent challenges in making an effective placental malaria vaccine.

New var reconstruction algorithm exposes high var sequence diversity in a single geographic location in Mali

Background

Encoded by the var gene family, highly variable Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) proteins mediate tissue-specific cytoadherence of infected erythrocytes, resulting in immune evasion and severe malaria disease. Sequencing and assembling the 40–60 var gene complement for individual infections has been notoriously difficult, impeding molecular epidemiological studies and the assessment of particular var elements as subunit vaccine candidates.

Methods

We developed and validated a novel algorithm, Exon-Targeted Hybrid Assembly (ETHA), to perform targeted assembly of var gene sequences, based on a combination of Pacific Biosciences and Illumina data.

Results

Using ETHA, we characterized the repertoire of var genes in 12 samples from uncomplicated malaria infections in children from a single Malian village and showed them to be as genetically diverse as vars from isolates from around the globe. The gene var2csa, a member of the var family associated with placental malaria pathogenesis, was present in each genome, as were vars previously associated with severe malaria.

Conclusion

ETHA, a tool to discover novel var sequences from clinical samples, will aid the understanding of malaria pathogenesis and inform the design of malaria vaccines based on PfEMP1.

ETHA is available at: https://sourceforge.net/projects/etha/.

Electronic supplementary material

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

In Vitro Variant Surface Antigen Expression in Plasmodium falciparum Parasites from a Semi-Immune Individual Is Not Correlated with Var Gene Transcription

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) is considered to be the main variant surface antigen (VSA) of Plasmodium falciparum and is mainly localized on electron-dense knobs in the membrane of the infected erythrocyte. Switches in PfEMP1 expression provide the basis for antigenic variation and are thought to be critical for parasite persistence during chronic infections. Recently, strain transcending anti-PfEMP1 immunity has been shown to develop early in life, challenging the role of PfEMP1 in antigenic variation during chronic infections. In this work we investigate how P. falciparum achieves persistence during a chronic asymptomatic infection. The infected individual (MOA) was parasitemic for 42 days and multilocus var gene genotyping showed persistence of the same parasite population throughout the infection. Parasites from the beginning of the infection were adapted to tissue culture and cloned by limiting dilution. Flow cytometry using convalescent serum detected a variable surface recognition signal on isogenic clonal parasites. Quantitative real-time PCR with a field isolate specific var gene primer set showed that the surface recognition signal was not correlated with transcription of individual var genes. Strain transcending anti-PfEMP1 immunity of the convalescent serum was demonstrated with CD36 selected and PfEMP1 knock-down NF54 clones. In contrast, knock-down of PfEMP1 did not have an effect on the antibody recognition signal in MOA clones. Trypsinisation of the membrane surface proteins abolished the surface recognition signal and immune electron microscopy revealed that antibodies from the convalescent serum bound to membrane areas without knobs and with knobs. Together the data indicate that PfEMP1 is not the main variable surface antigen during a chronic infection and suggest a role for trypsin sensitive non-PfEMP1 VSAs for parasite persistence in chronic infections.

Recombination events among virulence genes in malaria parasites are associated with G-quadruplex-forming DNA motifs

Background

Malaria parasites of the genus Plasmodium possess large hyper-variable families of antigen-encoding genes. These are often variantly-expressed and are major virulence factors for immune evasion and the maintenance of chronic infections. Recombination and diversification of these gene families occurs readily, and may be promoted by G-quadruplex (G4) DNA motifs within and close to the variant genes. G4s have been shown to cause replication fork stalling, DNA breakage and recombination in model systems, but these motifs remain largely unstudied in Plasmodium.

Results

We examined the nature and distribution of putative G4-forming sequences in multiple Plasmodium genomes, finding that their co-distribution with variant gene families is conserved across different Plasmodium species that have different types of variant gene families. In P. falciparum, where a large set of recombination events that occurred over time in cultured parasites has been mapped, we found a strong spatial association between these recombination events and putative G4-forming sequences. Finally, we searched Plasmodium genomes for the three classes of helicase that can unwind G4s: Plasmodium spp. have no identifiable homologue of the highly efficient G4 helicase PIF1, but they do encode two putative RecQ helicases and one homologue of the RAD3-family helicase FANCJ.

Conclusions

Our analyses, conducted at the whole-genome level in multiple species of Plasmodium, support the concept that G4s are likely to be involved in recombination and diversification of antigen-encoding gene families in this important protozoan pathogen.

Electronic supplementary material

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

Mosquito Vectors and the Globalization of Plasmodium falciparum Malaria

Plasmodium falciparum malaria remains a devastating public health problem. Recent discoveries have shed light on the origin and evolution of Plasmodium parasites and their interactions with their vertebrate and mosquito hosts. P. falciparum malaria originated in Africa from a single horizontal transfer between an infected gorilla and a human, and became global as the result of human migration. Today, P. falciparum malaria is transmitted worldwide by more than 70 different anopheline mosquito species. Recent studies indicate that the mosquito immune system can be a barrier to malaria transmission and that the P. falciparum Pfs47 gene allows the parasite to evade mosquito immune detection. Here, we review the origin and globalization of P. falciparum and integrate this history with analysis of the biology, evolution, and dispersal of the main mosquito vectors. This new perspective broadens our understanding of P. falciparum population structure and the dispersal of important parasite genetic traits.

Maintenance of phenotypic diversity within a set of virulence encoding genes of the malaria parasite Plasmodium falciparum

Infection by the human malaria parasite Plasmodium falciparum results in a broad spectrum of clinical outcomes, ranging from severe and potentially life-threatening malaria to asymptomatic carriage. In a process of naturally acquired immunity, individuals living in malaria-endemic regions build up a level of clinical protection, which attenuates infection severity in an exposure-dependent manner. Underlying this shift in the immunoepidemiology as well as the observed range in malaria pathogenesis is the var multigene family and the phenotypic diversity embedded within. The var gene-encoded surface proteins Plasmodium falciparum erythrocyte membrane protein 1 mediate variant-specific binding of infected red blood cells to a diverse set of host receptors that has been linked to specific disease manifestations, including cerebral and pregnancy-associated malaria. Here, we show that cross-reactive immune responses, which minimize the within-host benefit of each additionally expressed gene during infection, can cause selection for maximum phenotypic diversity at the genome level. We further show that differential functional constraints on protein diversification stably maintain uneven ratios between phenotypic groups, in line with empirical observation. Our results thus suggest that the maintenance of phenotypic diversity within P. falciparum is driven by an evolutionary trade-off that optimizes between within-host parasite fitness and between-host selection pressure.

Phagocytosis-inducing antibodies to Plasmodium falciparum upon immunization with a recombinant PfEMP1 NTS-DBL1α domain

BACKGROUND

Individuals living in endemic areas gradually acquire natural immunity to clinical malaria, largely dependent on antibodies against parasite antigens. There are many studies indicating that the variant antigen PfEMP1 at the surface of the parasitized red blood cell (pRBC) is one of the major targets of the immune response. It is believed that antibodies against PfEMP1 confer protection by blocking sequestration (rosetting and cytoadherence), inducing antibody-dependent cellular-inhibitory effect and opsonizing pRBCs for phagocytosis.

METHODS

A recombinant NTS-DBL1α domain from a rosette-mediating PfEMP1 was expressed in Escherichia coli. The resulting protein was purified and used for immunization to generate polyclonal (goat) and monoclonal (mouse) antibodies. The antibodies' ability to opsonize and induce phagocytosis in vitro was tested and contrasted with the presence of opsonizing antibodies naturally acquired during Plasmodium falciparum infection.

RESULTS

All antibodies recognized the recombinant antigen and the surface of live pRBCs, however, their capacity to opsonize the pRBCs for phagocytosis varied. The monoclonal antibodies isotyped as IgG2b did not induce phagocytosis, while those isotyped as IgG2a were in general very effective, inducing phagocytosis with similar levels as those naturally acquired during P. falciparum infection. These monoclonal antibodies displayed different patterns, some of them showing a concentration-dependent activity while others showed a prozone-like effect. The goat polyclonal antibodies were not able to induce phagocytosis.

CONCLUSION

Immunization with an NTS-DBL1-α domain of PfEMP1 generates antibodies that not only have a biological role in rosette disruption but also effectively induce opsonization for phagocytosis of pRBCs with similar activity to naturally acquired antibodies from immune individuals living in a malaria endemic area. Some of the antibodies with high opsonizing activity were not able to disrupt rosettes, indicating that epitopes of the NTS-DBL1-α other than those involved in rosetting are exposed on the pRBC surface and are able to induce functional antibodies. The ability to induce phagocytosis largely depended on the antibody isotype and on the ability to recognize the surface of the pRBC regardless of the rosette-disrupting capacity.

Serological Conservation of Parasite-Infected Erythrocytes Predicts Plasmodium falciparum Erythrocyte Membrane Protein 1 Gene Expression but Not Severity of Childhood Malaria

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1), expressed on P. falciparum-infected erythrocytes, is a major family of clonally variant targets of naturally acquired immunity to malaria. Previous studies have demonstrated that in areas where malaria is endemic, antibodies to infected erythrocytes from children with severe malaria tend to be more seroprevalent than antibodies to infected erythrocytes from children with nonsevere malaria. These data have led to a working hypothesis that PfEMP1 variants associated with parasite virulence are relatively conserved in structure. However, the longevity of such serologically conserved variants in the parasite population is unknown. Here, using infected erythrocytes from recently sampled clinical P. falciparum samples, we measured serological conservation using pools of antibodies in sera that had been sampled 10 to 12 years earlier. The serological conservation of infected erythrocytes strongly correlated with the expression of specific PfEMP1 subsets previously found to be associated with severe malaria. However, we found no association between serological conservation per se and disease severity within these data. This contrasts with the simple hypothesis that P. falciparum isolates with a serologically conserved group of PfEMP1 variants cause severe malaria. The data are instead consistent with periodic turnover of the immunodominant epitopes of PfEMP1 associated with severe malaria.

Mosquito Passage Dramatically Changes var Gene Expression in Controlled Human Plasmodium falciparum Infections

Virulence of the most deadly malaria parasite Plasmodium falciparum is linked to the variant surface antigen PfEMP1, which is encoded by about 60 var genes per parasite genome. Although the expression of particular variants has been associated with different clinical outcomes, little is known about var gene expression at the onset of infection. By analyzing controlled human malaria infections via quantitative real-time PCR, we show that parasite populations from 18 volunteers expressed virtually identical transcript patterns that were dominated by the subtelomeric var gene group B and, to a lesser extent, group A. Furthermore, major changes in composition and frequency of var gene transcripts were detected between the parental parasite culture that was used to infect mosquitoes and Plasmodia recovered from infected volunteers, suggesting that P. falciparum resets its var gene expression during mosquito passage and starts with the broad expression of a specific subset of var genes when entering the human blood phase.

Recombination and Diversification of the Variant Antigen Encoding Genes in the Malaria Parasite Plasmodium falciparum

The most severe form of human malaria is caused by the protozoan parasite Plasmodium falciparum. These parasites invade and replicate within the circulating red blood cells of infected individuals leading to numerous disease manifestations, including severe anemia, altered circulation, and tissue inflammation. Malaria parasites are also known for their ability to maintain a chronic infection through antigenic variation, the ability to systematically alter the antigens displayed on the surface of infected cells and thereby avoid clearance by the host's antibody response. The genome of P. falciparum includes several large, multicopy gene families that encode highly variable forms of the surface proteins that are the targets of host immunity. Alterations in expression of genes within these families are responsible for antigenic variation. This process requires the continuous generation of new antigenic variants within these gene families, and studies have shown that new variants arise through extensive recombination and gene conversion events between family members. Malaria parasites possess an unusual complement of DNA repair pathways, thus the study of recombination between variant antigen encoding genes provides a unique view into the evolution of mobile DNA in an organism distantly related to the more closely studied model eukaryotes.

Global selection of Plasmodium falciparum virulence antigen expression by host antibodies

Parasite proteins called PfEMP1 that are inserted on the surface of infected erythrocytes, play a key role in the severe pathology associated with infection by the Plasmodium falciparum malaria parasite. These proteins mediate binding of infected cells to the endothelial lining of blood vessels as a strategy to avoid clearance by the spleen and are major targets of naturally acquired immunity. PfEMP1 is encoded by a large multi-gene family called var. Mutually-exclusive transcriptional switching between var genes allows parasites to escape host antibodies. This study examined in detail the patterns of expression of var in a well-characterized sample of parasites from Kenyan Children. Instead of observing clear inverse relationships between the expression of broad sub-classes of PfEMP1, we found that expression of different PfEMP1 groups vary relatively independently. Parasite adaptation to host antibodies also appears to involve a general reduction in detectable var gene expression. We suggest that parasites switch both between different PfEMP1 variants and between high and low expression states. Such a strategy could provide a means of avoiding immunological detection and promoting survival under high levels of host immunity.

The role of PfEMP1 as targets of naturally acquired immunity to childhood malaria: prospects for a vaccine

The Plasmodium falciparum erythrocyte membrane protein 1 antigens that are inserted onto the surface of P. falciparum infected erythrocytes play a key role both in the pathology of severe malaria and as targets of naturally acquired immunity. They might be considered unlikely vaccine targets because they are extremely diverse. However, several lines of evidence suggest that underneath this molecular diversity there are a restricted set of epitopes which may act as effective targets for a vaccine against severe malaria. Here we review some of the recent developments in this area of research, focusing on work that has assessed the potential of these molecules as possible vaccine targets.

Antigenic Variation in Plasmodium falciparum

Plasmodium falciparum is the protozoan parasite that causes most malaria-associated morbidity and mortality in humans with over 500,000 deaths annually. The disease symptoms are associated with repeated cycles of invasion and asexual multiplication inside red blood cells of the parasite. Partial, non-sterile immunity to P. falciparum malaria develops only after repeated infections and continuous exposure. The successful evasion of the human immune system relies on the large repertoire of antigenically diverse parasite proteins displayed on the red blood cell surface and on the merozoite membrane where they are exposed to the human immune system. Expression switching of these polymorphic proteins between asexual parasite generations provides an efficient mechanism to adapt to the changing environment in the host and to maintain chronic infection. This chapter discusses antigenic diversity and variation in the malaria parasite and our current understanding of the molecular mechanisms that direct the expression of these proteins.

Comparative study of Plasmodium falciparum erythrocyte membrane protein 1-DBLα domain variants with respect to antigenic variations and docking interaction analysis with glycosaminoglycans

The variant surface antigen PfEMP1 (Plasmodium falciparum erythrocyte membrane protein 1) encoded by the polymorphic multi-copy var gene family plays an important role in parasite biology and the host-parasite interactions. Sequestration and antigenic variation is an essential component in the survival and pathogenesis of Plasmodium falciparum and contributes to chronic infection. The DBLα domain of PfEMP1 is a potential target for immuno-epidemiological studies and has been visualized as a vaccine candidate against severe malaria. Specific host receptors like heparin, heparan sulphate, blood group A and complement receptor 1 have been reported to bind the DBLα domain. Although heparin has been experimentally shown to disrupt the parasite-host interaction and effectively disrupt rosetting, the binding sites for the DBLα domain and the mechanism behind heparin-mediated rosette inhibition have not been elucidated. In this study, 3D structures and epitopes of the DBLα domain in 3D7 and in two Indian isolates have been predicted and compared. We have carried out docking studies on DBLα domains with human GAG receptors (heparin and heparan sulphate) to predict the strength of association between the protein-ligand interactions. The DBLα domain structures showed extensive diversity and polymorphism in their binding sites. The docking results indicate that heparin binds more effectively with high affinity as compared to heparan sulphate with some common interacting residues. These common residues can play an important role in rosetting and will aid in the designing of inhibitors specific to the interactions between DBLα and heparin or heparan sulphate would be important in malaria treatment. Thus it may lead to the development of novel interference strategies to block red blood cell invasion and provide protection against malaria.

Plasmodium falciparum-infected erythrocyte knob density is linked to the PfEMP1 variant expressed

Members of the clonally variant Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) family mediate adhesion of infected erythrocytes (IEs) to vascular receptors. PfEMP1 expression is normally confined to nanoscale knob protrusions on the IE surface membrane. To investigate the relationship between the densities of these IE surface knobs and the PfEMP1 variant expressed, we used specific antibody panning to generate three sublines of the P. falciparum clone IT4, which expresses the PfEMP1 variants IT4VAR04, IT4VAR32b, and IT4VAR60. The knob density in each subline was then determined by atomic force microscopy (AFM) and scanning electron microscopy (SEM) and compared to PfEMP1 and knob-associated histidine-rich protein (KAHRP) expression. Selection for uniform expression of IT4VAR04 produced little change in knob density, compared to unselected IEs. In contrast, selection for IT4VAR32b expression increased knob density approximately 3-fold, whereas IEs selected for IT4VAR60 expression were essentially knobless. When IT4VAR60⁺ IEs were subsequently selected to express IT4VAR04 or IT4VAR32b, they again displayed low and high knob densities, respectively. All sublines expressed KAHRP regardless of the PfEMP1 expressed. Our study documents for the first time that knob density is related to the PfEMP1 variant expressed. This may reflect topological requirements to ensure optimal adhesive properties of the IEs.

Infections with Plasmodium falciparum malaria parasites are still responsible for many deaths, especially among children and pregnant women. New interventions are needed to reduce severe illness and deaths caused by this malaria parasite. Thus, a better understanding of the mechanisms behind the pathogenesis is essential. A main reason why Plasmodium falciparum malaria is more severe than disease caused by other malaria species is its ability to express variant antigens on the infected erythrocyte surface. These antigens are presented on membrane protrusions known as knobs. This study set out to investigate the interplay between different variant antigens on the surface of P. falciparum-infected erythrocytes and the density of the knobs on which the antigens are expressed. Such a direct analysis of this relationship has not been reported before but adds to the important understanding of the complexity of malaria antigen presentation.