The South American Plasmodium falciparum var gene repertoire is limited, highly shared and possibly lacks several antigenic types

Antigenic strain diversity predicts different biogeographic patterns of maintenance and decline of antimalarial drug resistance

The establishment and spread of antimalarial drug resistance vary drastically across different biogeographic regions. Though most infections occur in sub-Saharan Africa, resistant strains often emerge in low-transmission regions. Existing models on resistance evolution lack consensus on the relationship between transmission intensity and drug resistance, possibly due to overlooking the feedback between antigenic diversity, host immunity, and selection for resistance. To address this, we developed a novel compartmental model that tracks sensitive and resistant parasite strains, as well as the host dynamics of generalized and antigen-specific immunity. Our results show a negative correlation between parasite prevalence and resistance frequency, regardless of resistance cost or efficacy. Validation using chloroquine-resistant marker data supports this trend. Post discontinuation of drugs, resistance remains high in low-diversity, low-transmission regions, while it steadily decreases in high-diversity, high-transmission regions. Our study underscores the critical role of malaria strain diversity in the biogeographic patterns of resistance evolution.

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.

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.

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.

A systematic review on genetic diversity of var gene DBL1α domain from different geographical regions in Plasmodium falciparum isolates

Background The major variant surface antigen (VSA) in Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) encoded by var gene family has an important role in cytoadhesion/sequestration and rosetting by adhesion of uninfected erythrocytes to infected erythrocytes leading to disease severity. DBL1α domain in the PfEMP-1, protein is crucial in the cytoadhesion phenomena in P. falciparum infections and this review aims to analyse the genetic diversity of DBL1α domain sequences in PfEMP-1 from different geographical regions globally. Methods All available DBL1α sequence data was reviewed by using the electronic database PubMed, ResearchGate, Google, Google scholar, MEDLINE with the following Keywords-Plasmodium falciparum", "var gene", "DBL1α", "field isolate", "diversity", "polymorphism", "Africa", "America", "Asia" and "Caribbean" from different geographical regions across the world. Results A total of 240 studies were identified initially but only 20 studies qualified for this systematic review. The overall ratio of distinct sequences DBL1α domain was 24.62/1167 the highest in African region (33.59/766 isolates) and lowest in South America (5.6/215 isolates). In the 18 included studies, the presence of distinct DBL1α sequences was the highest in Oceania 55.32% (1186/2144) followed by Africa (38.43%), Asia (22.45%) and South America (16.48%), though the sample size in Oceania was comparatively smaller to that of Africa and South America. Conclusion This review highlights the ratio and percentage of distinct sequences of DBL1α domain of var gene in different geographical regions giving an idea of the existing diversity prevalent in this potential vaccine target gene which may contribute to designing the preventive measures towards disease severity.

Evolutionary analyses of the major variant surface antigen-encoding genes reveal population structure of Plasmodium falciparum within and between continents

Malaria remains a major public health problem in many countries. Unlike influenza and HIV, where diversity in immunodominant surface antigens is understood geographically to inform disease surveillance, relatively little is known about the global population structure of PfEMP1, the major variant surface antigen of the malaria parasite Plasmodium falciparum. The complexity of the var multigene family that encodes PfEMP1 and that diversifies by recombination, has so far precluded its use in malaria surveillance. Recent studies have demonstrated that cost-effective deep sequencing of the region of var genes encoding the PfEMP1 DBLα domain and subsequent classification of within host sequences at 96% identity to define unique DBLα types, can reveal structure and strain dynamics within countries. However, to date there has not been a comprehensive comparison of these DBLα types between countries. By leveraging a bioinformatic approach (jumping hidden Markov model) designed specifically for the analysis of recombination within var genes and applying it to a dataset of DBLα types from 10 countries, we are able to describe population structure of DBLα types at the global scale. The sensitivity of the approach allows for the comparison of the global dataset to ape samples of Plasmodium Laverania species. Our analyses show that the evolution of the parasite population emerging out of Africa underlies current patterns of DBLα type diversity. Most importantly, we can distinguish geographic population structure within Africa between Gabon and Ghana in West Africa and Uganda in East Africa. Our evolutionary findings have translational implications in the context of globalization. Firstly, DBLα type diversity can provide a simple diagnostic framework for geographic surveillance of the rapidly evolving transmission dynamics of P. falciparum. It can also inform efforts to understand the presence or absence of global, regional and local population immunity to major surface antigen variants. Additionally, we identify a number of highly conserved DBLα types that are present globally that may be of biological significance and warrant further characterization.

Bayes-optimal estimation of overlap between populations of fixed size

Measuring the overlap between two populations is, in principle, straightforward. Upon fully sampling both populations, the number of shared objects—species, taxonomical units, or gene variants, depending on the context—can be directly counted. In practice, however, only a fraction of each population’s objects are likely to be sampled due to stochastic data collection or sequencing techniques. Although methods exists for quantifying population overlap under subsampled conditions, their bias is well documented and the uncertainty of their estimates cannot be quantified. Here we derive and validate a method to rigorously estimate the population overlap from incomplete samples when the total number of objects, species, or genes in each population is known, a special case of the more general β-diversity problem that is particularly relevant in the ecology and genomic epidemiology of malaria. By solving a Bayesian inference problem, this method takes into account the rates of subsampling and produces unbiased and Bayes-optimal estimates of overlap. In addition, it provides a natural framework for computing the uncertainty of its estimates, and can be used prospectively in study planning by quantifying the tradeoff between sampling effort and uncertainty.

Understanding when two populations are composed of similar species is important for ecologists, epidemiologists, and population geneticists, and in principle it is easy: just sample the two populations, compare the sets of species identified in each, and count how many appear in both populations. In practice, however, this is difficult because sampling methods typically produce only a random subset of the total population, leaving current population overlap estimates biased. Knowing only the number of shared members between two of these partial population samples, this paper shows how we can nevertheless estimate the true overlap between the full populations, when those full populations’ sizes are known. Using Bayesian statistics, we can also compute credible intervals to produce error bars. We show that using this unbiased approach has a dramatic impact on the conclusions one might draw from previously published studies in the malaria literature, which used simple but biased methods. Because the method in this paper quantifies the tradeoff between sampling effort and uncertainty, we also show how to compute the number of samples required to ensure high-confidence results, which may be useful for planning future studies or budgeting lab reagents and time.

Malaria during pregnancy and newborn outcome in an unstable transmission area in Brazil: A population-based record linkage study

BACKGROUND

Malaria in pregnancy (MiP) is one of the major causes of mortality and morbidity in tropical regions, causing maternal anemia, intrauterine growth retardation, preterm birth, and low birth weight (LBW). The integration of the information systems on pregnancy and malaria could prove to be a useful method of improved decision making for better maternal-child health.

METHODS

A population-based observational study acquired information retrospectively from all live births that occurred between 2006 and 2014 in Cruzeiro do Sul (Acre, Brazil). Social and clinical data of the mother and newborn was extracted from the Information System of Live Births. Malaria episodes information was obtained from the Brazilian Epidemiological Surveillance Information System Malaria. A deterministic record linkage was performed to assess malaria impact on pregnancy.

RESULTS

The studied population presented a malaria incidence of 8.9% (1283 pregnant women infected), of which 63.9% infected by Plasmodium (P.) vivax. Reduction of newborn birth weight at term (small for gestational age (SGA) and LBW) has been found associated with P. vivax infection during pregnancy (SGA-OR 1.24, 95% CI 1.02-1.52, p = 0.035; term LBW-OR 1.39, 95% CI 1.03-1.88, p = 0.033). Additionally, P. falciparum infection during pregnancy has been found to be associated with preterm births (OR 1.54, 95% CI 1.09-2.18, p = 0.016), which is related with late preterm births (OR 1.59, 95% CI 1.11-2.27, p = 0.011).

CONCLUSIONS

Despite the decrease of malaria cases during the evaluation period and regardless of Plasmodium species, we present evidence of the deleterious effects of MiP in a low transmission area in the Amazonian region.

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.

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.

Evolutionary structure of Plasmodium falciparum major variant surface antigen genes in South America: Implications for epidemic transmission and surveillance

Strong founder effects resulting from human migration out of Africa have led to geographic variation in single nucleotide polymorphisms (SNPs) and microsatellites (MS) of the malaria parasite, Plasmodium falciparum. This is particularly striking in South America where two major founder populations of P. falciparum have been identified that are presumed to have arisen from the transatlantic slave trade. Given the importance of the major variant surface antigen of the blood stages of P. falciparum as both a virulence factor and target of immunity, we decided to investigate the population genetics of the genes encoding “Plasmodium falciparum Erythrocyte Membrane Protein 1” (PfEMP1) among several countries in South America, in order to evaluate the transmission patterns of malaria in this continent. Deep sequencing of the DBLα domain of var genes from 128 P. falciparum isolates from five locations in South America was completed using a 454 high throughput sequencing protocol. Striking geographic variation in var DBLα sequences, similar to that seen for SNPs and MS markers, was observed. Colombia and French Guiana had distinct var DBLα sequences, whereas Peru and Venezuela showed an admixture. The importance of such geographic variation to herd immunity and malaria vaccination is discussed.

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.

Phylogeography of var gene repertoires reveals fine-scale geospatial clustering of Plasmodium falciparum populations in a highly endemic area

Plasmodium falciparum malaria is a major global health problem that is being targeted for progressive elimination. Knowledge of local disease transmission patterns in endemic countries is critical to these elimination efforts. To investigate fine-scale patterns of malaria transmission, we have compared repertoires of rapidly evolving var genes in a highly endemic area. A total of 3680 high-quality DBLα-sequences were obtained from 68 P. falciparum isolates from ten villages spread over two distinct catchment areas on the north coast of Papua New Guinea (PNG). Modelling of the extent of var gene diversity in the two parasite populations predicts more than twice as many var gene alleles circulating within each catchment (Mugil = 906; Wosera = 1094) than previously recognized in PNG (Amele = 369). In addition, there were limited levels of var gene sharing between populations, consistent with local parasite population structure. Phylogeographic analyses demonstrate that while neutrally evolving microsatellite markers identified population structure only at the catchment level, var gene repertoires reveal further fine-scale geospatial clustering of parasite isolates. The clustering of parasite isolates by village in Mugil, but not in Wosera was consistent with the physical and cultural isolation of the human populations in the two catchments. The study highlights the microheterogeneity of P. falciparum transmission in highly endemic areas and demonstrates the potential of var genes as markers of local patterns of parasite population structure.

Clonal outbreak of Plasmodium falciparum infection in eastern Panama

Identifying the source of resurgent parasites is paramount to a strategic, successful intervention for malaria elimination. Although the malaria incidence in Panama is low, a recent outbreak resulted in a 6-fold increase in reported cases. We hypothesized that parasites sampled from this epidemic might be related and exhibit a clonal population structure. We tested the genetic relatedness of parasites, using informative single-nucleotide polymorphisms and drug resistance loci. We found that parasites were clustered into 3 clonal subpopulations and were related to parasites from Colombia. Two clusters of Panamanian parasites shared identical drug resistance haplotypes, and all clusters shared a chloroquine-resistance genotype matching the pfcrt haplotype of Colombian origin. Our findings suggest these resurgent parasite populations are highly clonal and that the high clonality likely resulted from epidemic expansion of imported or vestigial cases. Malaria outbreak investigations that use genetic tools can illuminate potential sources of epidemic malaria and guide strategies to prevent further resurgence in areas where malaria has been eliminated.

Antibody recognition of Plasmodium falciparum infected red blood cells by symptomatic and asymptomatic individuals in the Brazilian Amazon

In the Amazon Region, there is a virtual absence of severe malaria and few fatal cases of naturally occurring Plasmodium falciparum infections; this presents an intriguing and underexplored area of research. In addition to the rapid access of infected persons to effective treatment, one cause of this phenomenon might be the recognition of cytoadherent variant proteins on the infected red blood cell (IRBC) surface, including the var gene encoded P. falciparum erythrocyte membrane protein 1. In order to establish a link between cytoadherence, IRBC surface antibody recognition and the presence or absence of malaria symptoms, we phenotype-selected four Amazonian P. falciparum isolates and the laboratory strain 3D7 for their cytoadherence to CD36 and ICAM1 expressed on CHO cells. We then mapped the dominantly expressed var transcripts and tested whether antibodies from symptomatic or asymptomatic infections showed a differential recognition of the IRBC surface. As controls, the 3D7 lineages expressing severe disease-associated phenotypes were used. We showed that there was no profound difference between the frequency and intensity of antibody recognition of the IRBC-exposed P. falciparum proteins in symptomatic vs. asymptomatic infections. The 3D7 lineages, which expressed severe malaria-associated phenotypes, were strongly recognised by most, but not all plasmas, meaning that the recognition of these phenotypes is frequent in asymptomatic carriers, but is not necessarily a prerequisite to staying free of symptoms.

Natural antibody response to Plasmodium falciparum merozoite antigens MSP5, MSP9 and EBA175 is associated to clinical protection in the Brazilian Amazon

BACKGROUND

Antibodies have an essential role in the acquired immune response against blood stage P. falciparum infection. Although several antigens have been identified as important antibody targets, it is still elusive which antigens have to be recognized for clinical protection. Herein, we analyzed antibodies from plasmas from symptomatic or asymptomatic individuals living in the same geographic area in the Western Amazon, measuring their recognition of multiple merozoite antigens.

METHODS

Specific fragments of genes encoding merozoite proteins AMA1 and members of MSP and EBL families from circulating P. falciparum field isolates present in asymptomatic and symptomatic patients were amplified by PCR. After cloning and expression of different versions of the antigens as recombinant GST-fusion peptides, we tested the reactivity of patients' plasmas by ELISA and the presence of IgG subclasses in the most reactive plasmas.

RESULTS

11 out of 24 recombinant antigens were recognized by plasmas from either symptomatic or asymptomatic infections. Antibodies to MSP9 (X2(DF=1) = 9.26/p = 0.0047) and MSP5 (X2(DF=1) = 8.29/p = 0.0069) were more prevalent in asymptomatic individuals whereas the opposite was observed for MSP1 block 2-MAD20 (X2(DF=1) = 6.41/p = 0.0206, Fisher's exact test). Plasmas from asymptomatic individuals reacted more intensely against MSP4 (U = 210.5, p < 0.03), MSP5 (U = 212, p < 0.004), MSP9 (U = 189.5, p < 0.002) and EBA175 (U = 197, p < 0.014, Mann-Whitney's U test). IgG1 and IgG3 were predominant for all antigens, but some patients also presented with IgG2 and IgG4. The recognition of MSP5 (OR = 0.112, IC95% = 0.021-0.585) and MSP9 (OR = 0.125, IC95% = 0.030-0.529, cross tab analysis) predicted 8.9 and 8 times less chances, respectively, to present symptoms. Higher antibody levels against MSP5 and EBA175 were associated by odds ratios of 9.4 (IC95% = 1.29-69.25) and 5.7 (IC95% = 1.12-29.62, logistic regression), respectively, with an asymptomatic status.

CONCLUSIONS

Merozoite antigens were targets of cytophilic antibodies and antibodies against MSP5, MSP9 and EBA175 were independently associated with decreased symptoms.

Historical shifts in Brazilian P. falciparum population structure and drug resistance alleles

Previous work suggests that Brazilian Plasmodium falciparum has limited genetic diversity and a history of bottlenecks, multiple reintroductions due to human migration, and clonal expansions. We hypothesized that Brazilian P. falciparum would exhibit clonal structure. We examined isolates collected across two decades from Amapá, Rondônia, and Pará state (n = 190). By examining more microsatellites markers on more chromosomes than previous studies, we hoped to define the extent of low diversity, linkage disequilibrium, bottlenecks, population structure, and parasite migration within Brazil. We used retrospective genotyping of samples from the 1980s and 1990s to explore the population genetics of SP resistant dhfr and dhps alleles. We tested an existing hypothesis that the triple mutant dhfr mutations 50R/51I/108N and 51I/108N/164L developed in southern Amazon from a single origin of common or similar parasites. We found that Brazilian P. falciparum had limited genetic diversity and isolation by distance was rejected, which suggests it underwent bottlenecks followed by migration between sites. Unlike Peru, there appeared to be gene flow across the Brazilian Amazon basin. We were unable to divide parasite populations by clonal lineages and pairwise FST were common. Most parasite diversity was found within sites in the Brazilian Amazon, according to AMOVA. Our results challenge the hypothesis that triple mutant alleles arose from a single lineage in the Southern Amazon. SP resistance, at both the double and triple mutant stages, developed twice and potentially in different regions of the Brazilian Amazon. We would have required samples from before the 1980s to describe how SP resistance spread across the basin or describe the complex internal migration of Brazilian parasites after the colonization efforts of past decades. The Brazilian Amazon basin may have sufficient internal migration for drug resistance reported in any particular region to rapidly spread to other parts of basin under similar drug pressure.

Diversity of the var gene family of Indonesian Plasmodium falciparum isolates

Background

The large polymorphic protein PfEMP1 is encoded by the var gene family. PfEMP1 has been shown to play an important role as cytoadherence ligand on the surface of infected erythrocytes and thereby contributes to the distinct pathogenesis of malaria. The study explored the diversity of the DBL1α and DBL2β-C2 domains of the protein from Indonesian Plasmodium falciparum field isolates.

Methods

Samples of patients with severe and uncomplicated malaria from two different malaria-endemic areas in Indonesia were collected and DNA directly extracted. Dried blood on filter paper was prepared for RNA extraction. PCR amplicons were either cloned and subsequently sequenced or directly sequenced for analysis on nucleotide and amino acid level. Recently published as well as self-designed primers were used for amplification.

Results

Blood from eight patients was finally used for analysis. Seventy-one different sequences out of over 500 DBL1α sequenced clones were observed, resulting in an average of 8.9 different DBL1α sequences per isolate. The average DBL1α sequence similarity within isolates was similar to between isolates. Phylogenetic analysis demonstrated no clustering of sequences regarding strain or geographical origin. The DBL1α sequences were analysed by distribution of semi-conserved features (cysteine/PoLV1-4 grouping) and classified into six sequence groups. The DBL1α cys2 type was observed in all expressed sequences in vivo. Expression of certain DBL sequences implied potential involvement in the pathogenesis. As expected, the DBL2β-C2 domains showed high to moderate homology among each other.

Conclusion

The DBL1α domains of PfEMP1 from clinical Indonesian isolates showed high divergence among same isolates and some similarities with other Asia-Pacific strains. Further investigations of important var gene domains with a larger sample size are required to confirm with statistical significance observed associations with severe malaria in Indonesian samples.

Genetic diversity of expressed Plasmodium falciparum var genes from Tanzanian children with severe malaria

Background

Severe malaria has been attributed to the expression of a restricted subset of the var multi-gene family, which encodes for Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1). PfEMP1 mediates cytoadherence and sequestration of infected erythrocytes into the post-capillary venules of vital organs such as the brain, lung or placenta. var genes are highly diverse and can be classified in three major groups (ups A, B and C) and two intermediate groups (B/A and B/C) based on the genomic location, gene orientation and upstream sequences. The genetic diversity of expressed var genes in relation to severity of disease in Tanzanian children was analysed.

Methods

Children with defined severe (SM) and asymptomatic malaria (AM) were recruited. Full-length var mRNA was isolated and reversed transcribed into var cDNA. Subsequently, the DBL and N-terminal domains, and up-stream sequences were PCR amplified, cloned and sequenced. Sequences derived from SM and AM isolates were compared and analysed.

Results

The analysis confirmed that the var family is highly diverse in natural Plasmodium falciparum populations. Sequence diversity of amplified var DBL-1α and upstream regions showed minimal overlap among isolates, implying that the var gene repertoire is vast and most probably indefinite in endemic areas. var DBL-1α sequences from AM isolates were more diverse with more singletons found (p<0.05) than those from SM infections. Furthermore, few var DBL-1α sequences from SM patients were rare and restricted suggesting that certain PfEMP1 variants might induce severe disease.

Conclusions

The genetic sequence diversity of var genes of P. falciparum isolates from Tanzanian children is large and its relationship to disease severity has been studied. Observed differences suggest that different var genes might have fundamentally different roles in the host-parasite interaction. Further research is required to examine clear disease-associations of var gene subsets in different geographical settings. The importance of very strict clinical definitions and appropriate large control groups needs to be emphasized for future studies on disease associations of PfEMP1.

Analysis of var genes cloned from a Plasmodium falciparum isolate in China

OBJECTIVE

To analyse the var gene repertoire and characterise the chondroitin sulphate A (CSA)-binding activity of the Duffy-binding like (DBL) domains encoded by the var2csa gene of a Plasmodium falciparum (P. falciparum) isolate in Hainan Province, China.

METHODS

The sequences of var DBL1 regions were PCR-amplified, sequenced and the sequence characteristics was bioinformatically analysed. Recombinant proteins encoded by the var2csa genes were expressed and purified. The binding activities of the recombinant proteins to CSA receptor was detected by ELISA assays.

RESULTS

Fifty six unique DBL α sequences were obtained, and the sequences represented similar diversity to the var genes of the genome parasite 3D7. There are two var2csa genes in the P. falciparum isolated from Hainan Province. Unlike in other falciparum parasites such as HB3, the two var2csa genes are more diverged. The receptor-binding capacity of DBL-5ε and DBL-6ε domains of HN var2CSA was studied.

CONCLUSIONS

This work represented the diversity of var genes of a P. falciparum isolate in China.

South American Plasmodium falciparum after the malaria eradication era: clonal population expansion and survival of the fittest hybrids

Malaria has reemerged in many regions where once it was nearly eliminated. Yet the source of these parasites, the process of repopulation, their population structure, and dynamics are ill defined. Peru was one of malaria eradication's successes, where Plasmodium falciparum was nearly eliminated for two decades. It reemerged in the 1990s. In the new era of malaria elimination, Peruvian P. falciparum is a model of malaria reinvasion. We investigated its population structure and drug resistance profiles. We hypothesized that only populations adapted to local ecological niches could expand and repopulate and originated as vestigial populations or recent introductions. We investigated the genetic structure (using microsatellites) and drug resistant genotypes of 220 parasites collected from patients immediately after peak epidemic expansion (1999-2000) from seven sites across the country. The majority of parasites could be grouped into five clonal lineages by networks and AMOVA. The distribution of clonal lineages and their drug sensitivity profiles suggested geographic structure. In 2001, artesunate combination therapy was introduced in Peru. We tested 62 parasites collected in 2006-2007 for changes in genetic structure. Clonal lineages had recombined under selection for the fittest parasites. Our findings illustrate that local adaptations in the post-eradication era have contributed to clonal lineage expansion. Within the shifting confluence of drug policy and malaria incidence, populations continue to evolve through genetic outcrossing influenced by antimalarial selection pressure. Understanding the population substructure of P. falciparum has implications for vaccine, drug, and epidemiologic studies, including monitoring malaria during and after the elimination phase.

From in vivo to in vitro: dynamic analysis of Plasmodium falciparum var gene expression patterns of patient isolates during adaptation to culture

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1), encoded by the var gene family, plays a crucial role in disease virulence through its involvement in binding to various host cellular receptors during infection. Growing evidence suggests that differential expression of the various var subgroups may be involved in parasite virulence. To further explore this issue, we have collected isolates from symptomatic patients in south China-Myanmar border, and characterized their sequence diversity and transcription profiles over time of var gene family, and cytoadherence properties from the time of their initial collection and extending through a two month period of adaptation to culture. Initially, we established a highly diverse, DBLα (4 cysteines) subtype-enriched, but unique local repertoire of var-DBL1α sequences by cDNA cloning and sequencing. Next we observed a rapid transcriptional decline of upsA- and upsB-subtype var genes at ring stage through qRT-PCR assays, and a switching event from initial ICAM-I binding to the CD36-binding activity during the first week of adaptive cultivation in vitro. Moreover, predominant transcription of upsA var genes was observed to be correlated with those isolates that showed a higher parasitemia at the time of collection and the ICAM-1-binding phenotype in culture. Taken together, these data indicate that the initial stage of adaptive process in vitro significantly influences the transcription of virulence-related var subtypes and expression of PfEMP1 variants. Further, the specific upregulation of the upsA var genes is likely linked to the rapid propagation of the parasite during natural infection due to the A-type PfEMP1 variant-mediated growth advantages.

A molecular epidemiological study of var gene diversity to characterize the reservoir of Plasmodium falciparum in humans in Africa

BACKGROUND

The reservoir of Plasmodium infection in humans has traditionally been defined by blood slide positivity. This study was designed to characterize the local reservoir of infection in relation to the diverse var genes that encode the major surface antigen of Plasmodium falciparum blood stages and underlie the parasite's ability to establish chronic infection and transmit from human to mosquito.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated the molecular epidemiology of the var multigene family at local sites in Gabon, Senegal and Kenya which differ in parasite prevalence and transmission intensity. 1839 distinct var gene types were defined by sequencing DBLα domains in the three sites. Only 76 (4.1%) var types were found in more than one population indicating spatial heterogeneity in var types across the African continent. The majority of var types appeared only once in the population sample. Non-parametric statistical estimators predict in each population at minimum five to seven thousand distinct var types. Similar diversity of var types was seen in sites with different parasite prevalences.

CONCLUSIONS/SIGNIFICANCE

Var population genomics provides new insights into the epidemiology of P. falciparum in Africa where malaria has never been conquered. In particular, we have described the extensive reservoir of infection in local African sites and discovered a unique var population structure that can facilitate superinfection through minimal overlap in var repertoires among parasite genomes. Our findings show that var typing as a molecular surveillance system defines the extent of genetic complexity in the reservoir of infection to complement measures of malaria prevalence. The observed small scale spatial diversity of var genes suggests that var genetics could greatly inform current malaria mapping approaches and predict complex malaria population dynamics due to the import of var types to areas where no widespread pre-existing immunity in the population exists.

Interrupting malaria transmission: quantifying the impact of interventions in regions of low to moderate transmission

Malaria has been eliminated from over 40 countries with an additional 39 currently planning for, or committed to, elimination. Information on the likely impact of available interventions, and the required time, is urgently needed to help plan resource allocation. Mathematical modelling has been used to investigate the impact of various interventions; the strength of the conclusions is boosted when several models with differing formulation produce similar data. Here we predict by using an individual-based stochastic simulation model of seasonal Plasmodium falciparum transmission that transmission can be interrupted and parasite reintroductions controlled in villages of 1,000 individuals where the entomological inoculation rate is <7 infectious bites per person per year using chemotherapy and bed net strategies. Above this transmission intensity bed nets and symptomatic treatment alone were not sufficient to interrupt transmission and control the importation of malaria for at least 150 days. Our model results suggest that 1) stochastic events impact the likelihood of successfully interrupting transmission with large variability in the times required, 2) the relative reduction in morbidity caused by the interventions were age-group specific, changing over time, and 3) the post-intervention changes in morbidity were larger than the corresponding impact on transmission. These results generally agree with the conclusions from previously published models. However the model also predicted changes in parasite population structure as a result of improved treatment of symptomatic individuals; the survival probability of introduced parasites reduced leading to an increase in the prevalence of sub-patent infections in semi-immune individuals. This novel finding requires further investigation in the field because, if confirmed, such a change would have a negative impact on attempts to eliminate the disease from areas of moderate transmission.