Genomes of cryptic chimpanzee Plasmodium species reveal key evolutionary events leading to human malaria

Ancient Plasmodium genomes shed light on the history of human malaria

Malaria-causing protozoa of the genus Plasmodium have exerted one of the strongest selective pressures on the human genome, and resistance alleles provide biomolecular footprints that outline the historical reach of these species1. Nevertheless, debate persists over when and how malaria parasites emerged as human pathogens and spread around the globe1,2. To address these questions, we generated high-coverage ancient mitochondrial and nuclear genome-wide data from P. falciparum, P. vivax and P. malariae from 16 countries spanning around 5,500 years of human history. We identified P. vivax and P. falciparum across geographically disparate regions of Eurasia from as early as the fourth and first millennia BCE, respectively; for P. vivax, this evidence pre-dates textual references by several millennia3. Genomic analysis supports distinct disease histories for P. falciparum and P. vivax in the Americas: similarities between now-eliminated European and peri-contact South American strains indicate that European colonizers were the source of American P. vivax, whereas the trans-Atlantic slave trade probably introduced P. falciparum into the Americas. Our data underscore the role of cross-cultural contacts in the dissemination of malaria, laying the biomolecular foundation for future palaeo-epidemiological research into the impact of Plasmodium parasites on human history. Finally, our unexpected discovery of P. falciparum in the high-altitude Himalayas provides a rare case study in which individual mobility can be inferred from infection status, adding to our knowledge of cross-cultural connectivity in the region nearly three millennia ago.

Plasmodium falciparum Glutamic Acid-Rich Protein-Independent Polyclonal Antibodies Inhibit Malaria Parasite Growth in Human Erythrocytes

Plasmodium falciparum glutamic acid-rich protein (PfGARP) is a recently characterized cell surface antigen encoded by Plasmodium falciparum, the causative agent of severe human malaria pathophysiology. Previously, we reported that the human erythrocyte band 3 (SLC4A1) serves as a host receptor for PfGARP. Antibodies against PfGARP did not affect parasite invasion and growth. We surmised that PfGARP may play a role in the rosetting and adhesion of malaria. Another study reported that antibodies targeting PfGARP exhibit potent inhibition of parasite growth. This inhibition occurred without the presence of any immune or complement components, suggesting the activation of an inherent density-dependent regulatory system. Here, we used polyclonal antibodies against PfGARP and a monoclonal antibody mAb7899 to demonstrate that anti-PfGARP polyclonal antibodies, but not mAb7899, exerted potent inhibition of parasite growth in infected erythrocytes independent of PfGARP. These findings suggest that an unknown malaria protein(s) is the target of growth arrest by polyclonal antibodies raised against PfGARP.

Role for gene conversion in the evolution of cell-surface antigens of the malaria parasite Plasmodium falciparum

While the malaria parasite Plasmodium falciparum has low average genome-wide diversity levels, likely due to its recent introduction from a gorilla-infecting ancestor (approximately 10,000 to 50,000 years ago), some genes display extremely high diversity levels. In particular, certain proteins expressed on the surface of human red blood cell-infecting merozoites (merozoite surface proteins (MSPs)) possess exactly 2 deeply diverged lineages that have seemingly not recombined. While of considerable interest, the evolutionary origin of this phenomenon remains unknown. In this study, we analysed the genetic diversity of 2 of the most variable MSPs, DBLMSP and DBLMSP2, which are paralogs (descended from an ancestral duplication). Despite thousands of available Illumina WGS datasets from malaria-endemic countries, diversity in these genes has been hard to characterise as reads containing highly diverged alleles completely fail to align to the reference genome. To solve this, we developed a pipeline leveraging genome graphs, enabling us to genotype them at high accuracy and completeness. Using our newly- resolved sequences, we found that both genes exhibit 2 deeply diverged lineages in a specific protein domain (DBL) and that one of the 2 lineages is shared across the genes. We identified clear evidence of nonallelic gene conversion between the 2 genes as the likely mechanism behind sharing, leading us to propose that gene conversion between diverged paralogs, and not recombination suppression, can generate this surprising genealogy; a model that is furthermore consistent with high diversity levels in these 2 genes despite the strong historical P. falciparum transmission bottleneck.

Expression and purification of active shikimate dehydrogenase from Plasmodium falciparum

Plasmodium falciparum is known to cause severe malaria, current treatment consists in artemisinin-based combination therapy, but resistance can lead to treatment failure. Knowledge concerning P. falciparum essential proteins can be used for searching new antimalarials, among these a potential candidate is shikimate dehydrogenase (SDH), an enzyme part of the shikimate pathway which is responsible for producing endogenous aromatic amino acids. SDH from P. falciparum (PfSDH) is unexplored by the scientific community, therefore, this study aims to establish the first protocol for active PfSDH expression. Putative PfSDH nucleotide sequence was used to construct an optimized expression vector pET28a+PfSDH inserted in E. coli BL21(DE3). As a result, optimal expression conditions were acquired by varying IPTG and temperature through time. Western Blot analysis was applied to verify appropriate PfSDH expression, solubilization and purification started with lysis followed by two-steps IMAC purification. Enzyme activity was measured spectrophotometrically by NADPH oxidation, optimal PfSDH expression occur at 0.1 mM IPTG for 48 hours growing at 37 °C and shaking at 200 rpm, recombinant PfSDH obtained after purification was soluble, pure and its physiological catalysis was confirmed. Thus, this study describes the first protocol for heterologous expression of PfSDH in soluble and active form.

Whole-genome surveillance identifies markers of Plasmodium falciparum drug resistance and novel genomic regions under selection in Mozambique

IMPORTANCE

Malaria is a devastating disease caused by Plasmodium parasites. The evolution of parasite drug resistance continues to hamper progress toward malaria elimination, and despite extensive efforts to control malaria, it remains a leading cause of death in Mozambique and other countries in the region. The development of successful vaccines and identification of molecular markers to track drug efficacy are essential for managing the disease burden. We present an analysis of the parasite genome in Mozambique, a country with one of the highest malaria burdens globally and limited available genomic data, revealing current selection pressure. We contribute additional evidence to limited prior studies supporting the effectiveness of SWGA in producing reliable genomic data from complex clinical samples. Our results provide the identity of genomic loci that may be associated with current antimalarial drug use, including artemisinin and lumefantrine, and reveal selection pressure predicted to compromise the efficacy of current vaccine candidates.

Neither the African-Centric S47 Nor P72 Variant of TP53 Is Associated With Reduced Risk of Febrile Malaria in a Malian Cohort Study

BACKGROUND

TP53 has been shown to play a role in inflammatory processes, including malaria. We previously found that p53 attenuates parasite-induced inflammation and predicts clinical protection to Plasmodium falciparum infection in Malian children. Here, we investigated whether p53 codon 47 and 72 polymorphisms are associated with differential risk of P. falciparum infection and uncomplicated malaria in a prospective cohort study of malaria immunity.

METHODS

p53 codon 47 and 72 polymorphisms were determined by sequencing TP53 exon 4 in 631 Malian children and adults enrolled in the Kalifabougou cohort study. The effects of these polymorphisms on the prospective risk of febrile malaria, incident parasitemia, and time to fever after incident parasitemia over 6 months of intense malaria transmission were assessed using Cox proportional hazards models.

RESULTS

Confounders of malaria risk, including age and hemoglobin S or C, were similar between individuals with or without p53 S47 and R72 polymorphisms. Relative to their respective common variants, neither S47 nor R72 was associated with differences in prospective risk of febrile malaria, incident parasitemia, or febrile malaria after parasitemia.

CONCLUSIONS

These findings indicate that p53 codon 47 and 72 polymorphisms are not associated with protection against incident P. falciparum parasitemia or uncomplicated febrile malaria.

A fast machine-learning-guided primer design pipeline for selective whole genome amplification

Addressing many of the major outstanding questions in the fields of microbial evolution and pathogenesis will require analyses of populations of microbial genomes. Although population genomic studies provide the analytical resolution to investigate evolutionary and mechanistic processes at fine spatial and temporal scales-precisely the scales at which these processes occur-microbial population genomic research is currently hindered by the practicalities of obtaining sufficient quantities of the relatively pure microbial genomic DNA necessary for next-generation sequencing. Here we present swga2.0, an optimized and parallelized pipeline to design selective whole genome amplification (SWGA) primer sets. Unlike previous methods, swga2.0 incorporates active and machine learning methods to evaluate the amplification efficacy of individual primers and primer sets. Additionally, swga2.0 optimizes primer set search and evaluation strategies, including parallelization at each stage of the pipeline, to dramatically decrease program runtime. Here we describe the swga2.0 pipeline, including the empirical data used to identify primer and primer set characteristics, that improve amplification performance. Additionally, we evaluate the novel swga2.0 pipeline by designing primer sets that successfully amplify Prevotella melaninogenica, an important component of the lung microbiome in cystic fibrosis patients, from samples dominated by human DNA.

Selective whole-genome amplification reveals population genetics of Leishmania braziliensis directly from patient skin biopsies

In Brazil, Leishmania braziliensis is the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) for Leishmania and show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multiple Leishmania species residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generate Leishmania genomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.

Role of Pfs47 in the dispersal of ancestral Plasmodium falciparum malaria through adaptation to different anopheline vectors

Plasmodium falciparum malaria originated when Plasmodium praefalciparum, a gorilla malaria parasite transmitted by African sylvan anopheline mosquitoes, adapted to humans. Pfs47, a protein on the parasite surface mediates P. falciparum evasion of the mosquito immune system by interacting with a midgut receptor and is critical for Plasmodium adaptation to different anopheline species. Genetic analysis of 4,971 Pfs47 gene sequences from different continents revealed that Asia and Papua New Guinea harbor Pfs47 haplotypes more similar to its ortholog in P. praefalciparum at sites that determine vector compatibility, suggesting that ancestral P. falciparum readily adapted to Asian vectors. Consistent with this observation, Pfs47-receptor gene sequences from African sylvan malaria vectors, such as Anopheles moucheti and An. marshallii, were found to share greater similarity with those of Asian vectors than those of vectors of the African An. gambiae complex. Furthermore, experimental infections provide direct evidence that transformed P. falciparum parasites carrying Pfs47 orthologs of P. praefalciparum or P. reichenowi were more effective at evading the immune system of the Asian malaria vector An. dirus than An. gambiae. We propose that high compatibility of ancestral P. falciparum Pfs47 with the receptors of Asian vectors facilitated the early dispersal of human malaria to the Asian continent, without having to first adapt to sub-Saharan vectors of the An. gambiae complex.

The ecology and epidemiology of malaria parasitism in wild chimpanzee reservoirs

Chimpanzees (Pan troglodytes) harbor rich assemblages of malaria parasites, including three species closely related to P. falciparum (sub-genus Laverania), the most malignant human malaria parasite. Here, we characterize the ecology and epidemiology of malaria infection in wild chimpanzee reservoirs. We used molecular assays to screen chimpanzee fecal samples, collected longitudinally and cross-sectionally from wild populations, for malaria parasite mitochondrial DNA. We found that chimpanzee malaria parasitism has an early age of onset and varies seasonally in prevalence. A subset of samples revealed Hepatocystis mitochondrial DNA, with phylogenetic analyses suggesting that Hepatocystis appears to cross species barriers more easily than Laverania. Longitudinal and cross-sectional sampling independently support the hypothesis that mean ambient temperature drives spatiotemporal variation in chimpanzee Laverania infection. Infection probability peaked at ~24.5 °C, consistent with the empirical transmission optimum of P. falciparum in humans. Forest cover was also positively correlated with spatial variation in Laverania prevalence, consistent with the observation that forest-dwelling Anophelines are the primary vectors. Extrapolating these relationships across equatorial Africa, we map spatiotemporal variation in the suitability of chimpanzee habitat for Laverania transmission, offering a hypothetical baseline indicator of human exposure risk.

Selective Whole-Genome Amplification as a Tool to Enrich Specimens with Low Treponema pallidum Genomic DNA Copies for Whole-Genome Sequencing

Syphilis is a sexually transmitted, disseminated acute and chronic infection caused by the bacterial pathogen Treponema pallidum subspecies pallidum . Primary syphilis typically presents as single or multiple mucocutaneous lesions and, if left untreated, can progress through multiple stages with various clinical manifestations.

Zoonotic origin of the human malaria parasite Plasmodium malariae from African apes

The human parasite Plasmodium malariae has relatives infecting African apes (Plasmodium rodhaini) and New World monkeys (Plasmodium brasilianum), but its origins remain unknown. Using a novel approach to characterise P. malariae-related sequences in wild and captive African apes, we found that this group comprises three distinct lineages, one of which represents a previously unknown, highly divergent species infecting chimpanzees, bonobos and gorillas across central Africa. A second ape-derived lineage is much more closely related to the third, human-infective lineage P. malariae, but exhibits little evidence of genetic exchange with it, and so likely represents a separate species. Moreover, the levels and nature of genetic polymorphisms in P. malariae indicate that it resulted from the zoonotic transmission of an African ape parasite, reminiscent of the origin of P. falciparum. In contrast, P. brasilianum falls within the radiation of human P. malariae, and thus reflects a recent anthroponosis.

Plasmodium malariae is a cause of malaria in humans and related species have been identified in non-human primates. Here, the authors use genomic analyses to establish that human P. malariae arose from a host switch of an ape parasite whilst a species infecting New World monkeys can be traced to a reverse zoonosis.

DNA in extracellular vesicles: from evolution to its current application in health and disease

Extracellular vesicle (EV) secretion is a highly conserved evolutionary trait in all organisms in the three domains of life. The packaging and release of EVs appears to be a bulk-flow process which takes place mainly under extreme conditions. EVs participate in horizontal gene transfer, which supports the survival of prokaryotic and eukaryotic microbes. In higher eukaryotes, almost all cells secrete a heterogeneous population of EVs loaded with various biomolecules. EV secretion is typically higher in cancer microenvironments, promoting tumor progression and metastasis. EVs are now recognized as additional mediators of autocrine and paracrine communication in health and disease. In this context, proteins and RNAs have been studied the most, but extracellular vesicle DNA (EV-DNA) has started to gain in importance in the last few years. In this review, we summarize new findings related to the loading mechanism(s), localization, and post-shedding function of EV-DNA. We also discuss the feasibility of using EV-DNA as a biomarker when performing a liquid biopsy, at the same time emphasizing the lack of data from clinical trials in this regard. Finally, we outline the potential of EV-DNA uptake and its interaction with the host genome as a promising tool for understanding the mechanisms of cancer evolution.

Protecting DNA in membrane vesicles seems to be a conserved phenomenon for the horizontal genetic flux between prokaryotes and lower eukaryotes.

Capturing and analyzing this vesicular DNA enables quick and non-invasive monitoring of natural ecosystems.

Cancer-derived extracellular vesicles containing DNA open up novel directions in cell-to-cell communication and therefore disease monitoring.

Complex and fluctuating conditions of the tumor microenvironment, mimicking natural ecosystems, could favor EV-DNA release, mediating tumor multi-clonal evolution and providing survival benefits.

Comparative analysis of the kinomes of Plasmodium falciparum, Plasmodium vivax and their host Homo sapiens

Background

Novel antimalarials should be effective across all species of malaria parasites that infect humans, especially the two species that bear the most impact, Plasmodium falciparum and Plasmodium vivax. Protein kinases encoded by pathogens, as well as host kinases required for survival of intracellular pathogens, carry considerable potential as targets for antimalarial intervention (Adderley et al. Trends Parasitol 37:508–524, 2021;  Wei et al. Cell Rep Med 2:100423, 2021). To date, no comprehensive P. vivax kinome assembly has been conducted; and the P. falciparum kinome, first assembled in 2004, requires an update. The present study, aimed to fill these gaps, utilises a recently published structurally-validated multiple sequence alignment (MSA) of the human kinome (Modi et al. Sci Rep 9:19790, 2019). This MSA is used as a scaffold to assist the alignment of all protein kinase sequences from P. falciparum and P. vivax, and (where possible) their assignment to specific kinase groups/families.

Results

We were able to assign six P. falciparum previously classified as OPK or ‘orphans’ (i.e. with no clear phylogenetic relation to any of the established ePK groups) to one of the aforementioned ePK groups. Direct phylogenetic comparison established that despite an overall high level of similarity between the P. falciparum and P. vivax kinomes, which will help in selecting targets for intervention, there are differences that may underlie the biological specificities of these species. Furthermore, we highlight a number of Plasmodium kinases that have a surprisingly high level of similarity with their human counterparts and therefore not well suited as targets for drug discovery.

Conclusions

Direct comparison of the kinomes of Homo sapiens, P. falciparum and P. vivax sheds additional light on the previously documented divergence of many P. falciparum and P. vivax kinases from those of their human host. We provide the first direct kinome comparison between the phylogenetically distinct species of P. falciparum and P. vivax, illustrating the key similarities and differences which must be considered in the context of kinase-directed antimalarial drug discovery, and discuss the divergences and similarities between the human and Plasmodium kinomes to inform future searches for selective antimalarial intervention.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08457-0.

Plasmodium development in Anopheles: a tale of shared resources

Interactions between the Anopheles mosquito vector and Plasmodium parasites shape how malaria is transmitted in endemic regions. The long association of these two organisms has led to evolutionary processes that minimize fitness costs of infection and benefit both players through shared nutrient resources, parasite immune suppression, and mosquito tolerance to infection. In this review we explore recent data describing how Plasmodium falciparum, the deadliest malaria parasite, associates with one of its most important natural mosquito hosts, Anopheles gambiae, and we discuss the implications of these findings for parasite transmission and vector control strategies currently in development.

Global diversity of the gene encoding the Pfs25 protein-a Plasmodium falciparum transmission-blocking vaccine candidate

Background

Vaccines against the sexual stages of the malarial parasite Plasmodium falciparum are indispensable for controlling malaria and abrogating the spread of drug-resistant parasites. Pfs25, a surface antigen of the sexual stage of P. falciparum, is a leading candidate for transmission-blocking vaccine development. While clinical trials have reported that Pfs25-based vaccines are safe and effective in inducing transmission-blocking antibodies, the extent of the genetic diversity of Pfs25 in malaria endemic populations has rarely been studied. Thus, this study aimed to investigate the global diversity of Pfs25 in P. falciparum populations.

Methods

A database of 307 Pfs25 sequences of P. falciparum was established. Population genetic analyses were performed to evaluate haplotype and nucleotide diversity, analyze haplotypic distribution patterns of Pfs25 in different geographical populations, and construct a haplotype network. Neutrality tests were conducted to determine evidence of natural selection. Homology models of the Pfs25 haplotypes were constructed, subjected to molecular dynamics (MD), and analyzed in terms of flexibility and percentages of secondary structures.

Results

The Pfs25 gene of P. falciparum was found to have 11 unique haplotypes. Of these, haplotype 1 (H1) and H2, the major haplotypes, represented 70% and 22% of the population, respectively, and were dominant in Asia, whereas only H1 was dominant in Africa, Central America, and South America. Other haplotypes were rare and region-specific, resulting in unique distribution patterns in different geographical populations. The diversity in Pfs25 originated from ten single-nucleotide polymorphism (SNP) loci located in the epidermal growth factor (EGF)-like domains and anchor domain. Of these, an SNP at position 392 (GGA/GCA), resulting in amino acid substitution 131 (Gly/Ala), defined the two major haplotypes. The MD results showed that the structures of H1 and H2 variants were relatively similar. Limited polymorphism in Pfs25 could likely be due to negative selection.

Conclusions

The study successfully established a Pfs25 sequence database that can become an essential tool for monitoring vaccine efficacy, designing assays for detecting malaria carriers, and conducting epidemiological studies of P. falciparum. The discovery of the two major haplotypes, H1 and H2, and their conserved structures suggests that the current Pfs25-based vaccines could be used globally for malaria control.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05078-6.

Population genomics and evidence of clonal replacement of Plasmodium falciparum in the Peruvian Amazon

Previous studies have shown that P. falciparum parasites in South America have undergone population bottlenecks resulting in clonal lineages that are differentially distributed and that have been responsible for several outbreaks different endemic regions. In this study, we explored the genomic profile of 18 P. falciparum samples collected in the Peruvian Amazon Basin (Loreto) and 6 from the Peruvian North Coast (Tumbes). Our results showed the presence of three subpopulations that matched previously typed lineages in Peru: Bv1 (n = 17), Clonet D (n = 4) and Acre-Loreto type (n = 3). Gene coverage analysis showed that none of the Bv1 samples presented coverage for pfhrp2 and pfhrp3. Genotyping of drug resistance markers showed a high prevalence of Chloroquine resistance mutations S1034C/N1042D/D1246Y in pfmdr1 (62.5%) and K45T in pfcrt (87.5%). Mutations associated with sulfadoxine and pyrimethamine treatment failure were found on 88.8% of the Bv1 samples which were triple mutants for pfdhfr (50R/51I/108N) and pfdhps (437G/540E/581G). Analysis of the pfS47 gene that allows P. falciparum to evade mosquito immune responses showed that the Bv1 lineage presented one pfS47 haplotype exclusive to Loreto and another haplotype that was present in both Loreto and Tumbes. Furthermore, a possible expansion of Bv1 was detected since 2011 in Loreto. This replacement could be a result of the high prevalence of CQ resistance polymorphisms in Bv1, which could have provided a selective advantage to the indirect selection pressures driven by the use of CQ for P. vivax treatment.

Common host variation drives malaria parasite fitness in healthy human red cells

The replication of Plasmodium falciparum parasites within red blood cells (RBCs) causes severe disease in humans, especially in Africa. Deleterious alleles like hemoglobin S are well-known to confer strong resistance to malaria, but the effects of common RBC variation are largely undetermined. Here, we collected fresh blood samples from 121 healthy donors, most with African ancestry, and performed exome sequencing, detailed RBC phenotyping, and parasite fitness assays. Over one-third of healthy donors unknowingly carried alleles for G6PD deficiency or hemoglobinopathies, which were associated with characteristic RBC phenotypes. Among non-carriers alone, variation in RBC hydration, membrane deformability, and volume was strongly associated with P. falciparum growth rate. Common genetic variants in PIEZO1, SPTA1/SPTB, and several P. falciparum invasion receptors were also associated with parasite growth rate. Interestingly, we observed little or negative evidence for divergent selection on non-pathogenic RBC variation between Africans and Europeans. These findings suggest a model in which globally widespread variation in a moderate number of genes and phenotypes modulates P. falciparum fitness in RBCs.

A population genetic perspective on the origin, spread and adaptation of the human malaria agents Plasmodium falciparum and Plasmodium vivax

Malaria is considered one of the most important scourges that humanity has faced during its history, being responsible every year for numerous deaths worldwide. The disease is caused by protozoan parasites, among which two species are responsible of the majority of the burden, Plasmodium falciparum and Plasmodium vivax. For these two parasite species, the questions of their origin (how and when they appeared in humans), of their spread throughout the world, as well as how they have adapted to humans have long been of interest to the scientific community. Here, we review the current knowledge that has accumulated on these different questions, thanks in particular to the analysis of the genetic and genomic variability of these parasites and comparison with related Plasmodium species infecting other host species (like non-human primates). In this paper we review the existing body of knowledge, including current research dealing with these questions, focusing particularly on genetic analysis and genomic variability of these parasites and comparison with related Plasmodium species infecting other species of host (such as non-human primates).

Unique Endomembrane Systems and Virulence in Pathogenic Protozoa

Virulence in pathogenic protozoa is often tied to secretory processes such as the expression of adhesins on parasite surfaces or the secretion of proteases to assisted in tissue invasion and other proteins to avoid the immune system. This review is a broad overview of the endomembrane systems of pathogenic protozoa with a focus on Giardia, Trichomonas, Entamoeba, kinetoplastids, and apicomplexans. The focus is on unique features of these protozoa and how these features relate to virulence. In general, the basic elements of the endocytic and exocytic pathways are present in all protozoa. Some of these elements, especially the endosomal compartments, have been repurposed by the various species and quite often the repurposing is associated with virulence. The Apicomplexa exhibit the most unique endomembrane systems. This includes unique secretory organelles that play a central role in interactions between parasite and host and are involved in the invasion of host cells. Furthermore, as intracellular parasites, the apicomplexans extensively modify their host cells through the secretion of proteins and other material into the host cell. This includes a unique targeting motif for proteins destined for the host cell. Most notable among the apicomplexans is the malaria parasite, which extensively modifies and exports numerous proteins into the host erythrocyte. These modifications of the host erythrocyte include the formation of unique membranes and structures in the host erythrocyte cytoplasm and on the erythrocyte membrane. The transport of parasite proteins to the host erythrocyte involves several unique mechanisms and components, as well as the generation of compartments within the erythrocyte that participate in extraparasite trafficking.

Ape Origins of Human Malaria

African apes harbor at least twelve Plasmodium species, some of which have been a source of human infection. It is now well established that Plasmodium falciparum emerged following the transmission of a gorilla parasite, perhaps within the last 10,000 years, while Plasmodium vivax emerged earlier from a parasite lineage that infected humans and apes in Africa before the Duffy-negative mutation eliminated the parasite from humans there. Compared to their ape relatives, both human parasites have greatly reduced genetic diversity and an excess of nonsynonymous mutations, consistent with severe genetic bottlenecks followed by rapid population expansion. A putative new Plasmodium species widespread in chimpanzees, gorillas, and bonobos places the origin of Plasmodium malariae in Africa. Here, we review what is known about the origins and evolutionary history of all human-infective Plasmodium species, the time and circumstances of their emergence, and the diversity, host specificity, and zoonotic potential of their ape counterparts.

Identification of a protein unique to the genus Plasmodium that contains a WD40 repeat domain and extensive low-complexity sequence

Proteins containing WD40 domains play important roles in the formation of multiprotein complexes. Little is known about WD40 proteins in the malaria parasite. This report contains the initial description of a WD40 protein that is unique to the genus Plasmodium and possibly closely related genera. The N-terminal portion of this protein consists of seven WD40 repeats that are highly conserved in all Plasmodium species. Following the N-terminal region is a central region that is conserved within the major Plasmodium clades, such as parasites of great apes, monkeys, rodents, and birds, but partially conserved across all Plasmodium species. This central region contains extensive low-complexity sequence and is predicted to have a disordered structure. Proteins with disordered structure generally function in molecular interactions. The C-terminal region is semi-conserved across all Plasmodium species and has no notable features. This WD40 repeat protein likely functions in some aspect of parasite biology that is unique to Plasmodium and this uniqueness makes the protein a possible target for therapeutic intervention.

Malaria in the 'Omics Era'

Genomics has revolutionised the study of the biology of parasitic diseases. The first Eukaryotic parasite to have its genome sequenced was the malaria parasite Plasmodium falciparum. Since then, Plasmodium genomics has continued to lead the way in the study of the genome biology of parasites, both in breadth—the number of Plasmodium species’ genomes sequenced—and in depth—massive-scale genome re-sequencing of several key species. Here, we review some of the insights into the biology, evolution and population genetics of Plasmodium gained from genome sequencing, and look at potential new avenues in the future genome-scale study of its biology.

Multi-omics approaches to improve malaria therapy

Malaria contributes to the most widespread infectious diseases worldwide. Even though current drugs are commercially available, the ever-increasing drug resistance problem by malaria parasites poses new challenges in malaria therapy. Hence, searching for efficient therapeutic strategies is of high priority in malaria control. In recent years, multi-omics technologies have been extensively applied to provide a more holistic view of functional principles and dynamics of biological mechanisms. We briefly review multi-omics technologies and focus on recent malaria progress conducted with the help of various omics methods. Then, we present up-to-date advances for multi-omics approaches in malaria. Next, we describe resistance phenomena to established antimalarial drugs and underlying mechanisms. Finally, we provide insight into novel multi-omics approaches, new drugs and vaccine developments and analyze current gaps in multi-omics research. Although multi-omics approaches have been successfully used in malaria studies, they are still limited. Many gaps need to be filled to bridge the gap between basic research and treatment of malaria patients. Multi-omics approaches will foster a better understanding of the molecular mechanisms of Plasmodium that are essential for the development of novel drugs and vaccines to fight this disastrous disease.

Genomic analysis reveals independent evolution of Plasmodium falciparum populations in Ethiopia

Background

Plasmodium falciparum parasite populations in Ethiopia have been experiencing local selective pressures from drugs and immunity, leading to evolutionary adaptation. However, there was a paucity of data on genomic characterization and evolutionary adaptations of P. falciparum isolates from the central area of Ethiopia.

Methods

Whole-genome analysis of 25 P. falciparum isolates from central Ethiopia, specifically from West Arsi, were studied to determine their genetic diversity, population structures, and signatures of selection in known drug resistance alleles against global isolates from Cambodia, Thailand, DR Congo, and Malawi.

Results

A total of 18,517 high-quality single-nucleotide polymorphisms (SNPs) were identified in Ethiopian P. falciparum isolates. About 84% of the Ethiopian P. falciparum isolates had a F_WS value > 0.95 showing a dominant single genotype infection in most isolates at the time of collection with little potential for out-crossing as expected in areas with low transmission intensity. Within-host diversity of Ethiopian infections was significantly different from East African (p < 0.001), but not Southeast Asian infections (P > 0.05). A significant population structure has been observed by PCA and population differentiation between Ethiopian parasites and East African (Fst ~ 10%) and Southeast Asian populations (Fst ~ 18%), suggesting limited gene flow and the independent evolution of the Ethiopian parasite population. Moreover, a total of 125 genes under balancing selection was found that include ama1, trap, eba175, and lsa3, previously identified as targets of human host immunity. Recent directional selection analysis using integrated standardized haplotype score (IHS) did not detect any selection signatures in the Pfcrt, Pfdhfr, Pfdhps, Pfmdr1, and PfK13 genes. However, known drug resistance-conferring mutations analysis showed that at least one SNP marker was fixed in these genes, but not in Pfdhps and PfK13.

Conclusion

Plasmodium falciparum populations in the central region of Ethiopia was structurally diverged from both Southeast Asian and other East African populations. Malaria infections in Ethiopia had low within-host diversity, and parasites carry fixed chloroquine resistance markers despite the withdrawal of this drug for the treatment of P. falciparum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03660-y.

Optimization of whole-genome sequencing of Plasmodium falciparum from low-density dried blood spot samples

Background

Whole-genome sequencing (WGS) is becoming increasingly useful to study the biology, epidemiology, and ecology of malaria parasites. Despite ease of sampling, DNA extracted from dried blood spots (DBS) has a high ratio of human DNA compared to parasite DNA, which poses a challenge for downstream genetic analyses. The effects of multiple methods for DNA extraction, digestion of methylated DNA, and amplification were evaluated on the quality and fidelity of WGS data recovered from DBS.

Methods

Low parasite density mock DBS samples were created, extracted either with Tween-Chelex or QIAamp, treated with or without McrBC, and amplified with one of three different amplification techniques (two sWGA primer sets and one rWGA). Extraction conditions were evaluated on performance of sequencing depth, percentiles of coverage, and expected SNP concordance.

Results

At 100 parasites/μL, Chelex-Tween-McrBC samples had higher coverage (5 × depth = 93% genome) than QIAamp extracted samples (5 × depth = 76% genome). The two evaluated sWGA primer sets showed minor differences in overall genome coverage and SNP concordance, with a newly proposed combination of 20 primers showing a modest improvement in coverage over those previously published.

Conclusions

Overall, Tween-Chelex extracted samples that were treated with McrBC digestion and are amplified using 6A10AD sWGA conditions had minimal dropout rate, higher percentages of coverage at higher depth, and more accurate SNP concordance than QiaAMP extracted samples. These findings extend the results of previously reported methods, making whole genome sequencing accessible to a larger number of low density samples that are commonly encountered in cross-sectional surveys.

Parasite and Host Erythrocyte Kinomics of Plasmodium Infection

Malaria remains a heavy public health and socioeconomic burden in tropical and subtropical regions. Increasing resistance against front-line treatments implies that novel targets for antimalarial intervention are urgently required. Protein kinases of both the parasites and their host cells possess strong potential in this respect. We present an overview of the updated kinome of Plasmodium falciparum, the species that is the largest contributor to malaria mortality, and of current knowledge pertaining to the function of parasite-encoded protein kinases during the parasite's life cycle. Furthermore, we detail recent advances in drug initiatives targeting Plasmodium kinases and outline the potential of protein kinases in the context of the growing field of host-directed therapies, which is currently being explored as a novel way to combat parasite drug resistance.

A Sarcina bacterium linked to lethal disease in sanctuary chimpanzees in Sierra Leone

Human and animal infections with bacteria of the genus Sarcina (family Clostridiaceae) are associated with gastric dilation and emphysematous gastritis. However, the potential roles of sarcinae as commensals or pathogens remain unclear. Here, we investigate a lethal disease of unknown etiology that affects sanctuary chimpanzees (Pan troglodytes verus) in Sierra Leone. The disease, which we have named "epizootic neurologic and gastroenteric syndrome" (ENGS), is characterized by neurologic and gastrointestinal signs and results in death of the animals, even after medical treatment. Using a case-control study design, we show that ENGS is strongly associated with Sarcina infection. The microorganism is distinct from Sarcina ventriculi and other known members of its genus, based on bacterial morphology and growth characteristics. Whole-genome sequencing confirms this distinction and reveals the presence of genetic features that may account for the unusual virulence of the bacterium. Therefore, we propose that this organism be considered the representative of a new species, named "Candidatus Sarcina troglodytae". Our results suggest that a heretofore unrecognized complex of related sarcinae likely exists, some of which may be highly virulent. However, the potential role of "Ca. S. troglodytae" in the etiology of ENGS, alone or in combination with other factors, remains a topic for future research.

Sub-Saharan African information potential to unveil adaptations to infectious disease

Sub-Saharan Africa is the most promising region of the world to conduct high-throughput studies to unveil adaptations to infectious diseases due to several reasons, namely, the longest evolving time-depth in the Homo sapiens phylogenetic tree (at least two-third older than any other worldwide region); the continuous burden of infectious diseases (still number one in health/life threat); and the coexistence of populations practising diverse subsistence modes (nomadic or seminomadic hunter-gatherers and agropastoralists, and sedentary agriculturalists, small urban and megacity groups). In this review, we will present the most up-to-date results that shed light on three main hypotheses related with this adaptation. One is the hypothesis of coevolution between host and pathogen, given enough time for the establishment of this highly dynamic relationship. The second hypothesis enunciates that the agricultural transition was responsible for the increase of the infectious disease burden, due to the huge expansion of the sedentary human population and the cohabitation with domesticates as main reservoirs of pathogens. The third hypothesis states that the boosting of our immune system against pathogens by past selection may have resulted in maladaptation of the developed hygienic societies, leading to an increase of allergic, inflammatory and autoimmune disorders. Further work will enlighten the biological mechanisms behind these main adaptations, which can be insightful for translation into diagnosis, prognosis and treatment interventions.

Zoonotic Transmissions and Host Switches of Malaria Parasites

Malaria is a deadly disease that affects the health of hundreds of millions of people annually. There are five Plasmodium parasite species that can naturally infect humans, including Plasmodium falciparum, Plasmodium vivax, Plasmodium malariae, Plasmodium ovale and Plasmodium knowlesi. Some of the parasites can also infect various non-human primates. Parasites mainly infecting monkeys such as Plasmodium cynomolgi (in fact P. knowlesi was considered as a parasite of monkeys for years) can also be transmitted to human hosts. Recently, many new Plasmodium species were discovered in African apes, and it is possible that some of the parasites can be transmitted to humans in the future. Here, we searched PubMed and the internet via Google and selected articles concerning zoonotic transmission and evolution of selected malaria parasite species. We reviewed the current advances in the relevant topics emphasizing on transmissions of malaria parasites between humans and non-human primates. We also briefly discuss the transmissions of some avian malaria parasites between wild birds and domestic fowls. Zoonotic malaria transmissions are widespread, which poses a threat to public health. More studies on parasite species identification in non-human primates, transmission, and evolution are needed to reduce or prevent transmission of malaria parasites from non-human primates to humans.

A mating-induced reproductive gene promotes Anopheles tolerance to Plasmodium falciparum infection

Anopheles mosquitoes have transmitted Plasmodium parasites for millions of years, yet it remains unclear whether they suffer fitness costs to infection. Here we report that the fecundity of virgin and mated females of two important vectors—Anopheles gambiae and Anopheles stephensi—is not affected by infection with Plasmodium falciparum, demonstrating that these human malaria parasites do not inflict this reproductive cost on their natural mosquito hosts. Additionally, parasite development is not impacted by mating status. However, in field studies using different P. falciparum isolates in Anopheles coluzzii, we find that Mating-Induced Stimulator of Oogenesis (MISO), a female reproductive gene strongly induced after mating by the sexual transfer of the steroid hormone 20-hydroxyecdysone (20E), protects females from incurring fecundity costs to infection. MISO-silenced females produce fewer eggs as they become increasingly infected with P. falciparum, while parasite development is not impacted by this gene silencing. Interestingly, previous work had shown that sexual transfer of 20E has specifically evolved in Cellia species of the Anopheles genus, driving the co-adaptation of MISO. Our data therefore suggest that evolution of male-female sexual interactions may have promoted Anopheles tolerance to P. falciparum infection in the Cellia subgenus, which comprises the most important malaria vectors.

Plasmodium falciparum, the deadliest form of human malaria, is transmitted when female Anopheles mosquitoes bite people and take a blood meal in order to develop eggs. To date, it is still poorly understood whether Anopheles mosquitoes that get infected with P. falciparum suffer fitness costs. Here, we find that the number of eggs produced by Anopheles gambiae and Anopheles stephensi females is not affected by P. falciparum infection, and that the mating status of the mosquitoes does not impact the parasite. However, in field experiments infecting a related species, Anopheles coluzzii, with P. falciparum using blood from donors in Burkina Faso, we find that interfering with the expression of a gene normally triggered by the sexual transfer of the steroid hormone 20-hydroxyecdysone induces increasing costs to egg development as females become more infected with P. falciparum, with no impacts on the parasite. The results of our study suggest that pathways triggered by mating may help Anopheles prevent reproductive costs associated with P. falciparum infection, providing new insights into evolutionary strategies adopted by anophelines in the face of a longstanding association with Plasmodium parasites.

Full Molecular Typing of Neisseria meningitidis Directly from Clinical Specimens for Outbreak Investigation

Neisseria meningitidis is a leading cause of bacterial meningitis and sepsis worldwide and an occasional cause of meningococcal urethritis. When isolates are unavailable for surveillance or outbreak investigations, molecular characterization of pathogens needs to be performed directly from clinical specimens, such as cerebrospinal fluid (CSF), blood, or urine. However, genome sequencing of specimens is challenging because of low bacterial and high human DNA abundances. We developed selective whole-genome amplification (SWGA), an isothermal multiple-displacement amplification-based method, to efficiently enrich, sequence, and de novo assemble N. meningitidis DNA from clinical specimens with low bacterial loads. SWGA was validated with 12 CSF specimens from invasive meningococcal disease cases and 12 urine specimens from meningococcal urethritis cases. SWGA increased the mean proportion of N. meningitidis reads by 2 to 3 orders of magnitude, enabling identification of at least 90% of the 1,605 N. meningitidis core genome loci for 50% of the specimens. The validated method was used to investigate two meningitis outbreaks recently reported in Togo and Burkina Faso. Twenty-seven specimens with low bacterial loads were processed by SWGA before sequencing, and 12 of 27 were successfully assembled to obtain the full molecular typing and vaccine antigen profile of the N. meningitidis pathogen, thus enabling thorough characterization of outbreaks. This method is particularly important for enhancing molecular surveillance in regions with low culture rates. SWGA produces enough reads for phylogenetic and allelic analysis at a low cost. More importantly, the procedure can be extended to enrich other important human bacterial pathogens.

Host-Malaria Parasite Interactions and Impacts on Mutual Evolution

Malaria is the most deadly parasitic disease, affecting hundreds of millions of people worldwide. Malaria parasites have been associated with their hosts for millions of years. During the long history of host-parasite co-evolution, both parasites and hosts have applied pressure on each other through complex host-parasite molecular interactions. Whereas the hosts activate various immune mechanisms to remove parasites during an infection, the parasites attempt to evade host immunity by diversifying their genome and switching expression of targets of the host immune system. Human intervention to control the disease such as antimalarial drugs and vaccination can greatly alter parasite population dynamics and evolution, particularly the massive applications of antimalarial drugs in recent human history. Vaccination is likely the best method to prevent the disease; however, a partially protective vaccine may have unwanted consequences that require further investigation. Studies of host-parasite interactions and co-evolution will provide important information for designing safe and effective vaccines and for preventing drug resistance. In this essay, we will discuss some interesting molecules involved in host-parasite interactions, including important parasite antigens. We also discuss subjects relevant to drug and vaccine development and some approaches for studying host-parasite interactions.

Genetic Markers of Adaptation of Plasmodium falciparum to Transmission by American Vectors Identified in the Genomes of Parasites from Haiti and South America

The malaria parasite, Plasmodium falciparum, was introduced into Hispaniola and other regions of the Americas through the slave trade spanning the 16th through the 19th centuries. During this period, more than 12 million Africans were brought across the Atlantic to the Caribbean and other regions of the Americas. Since malaria is holoendemic in West Africa, a substantial percentage of these individuals carried the parasite. St. Domingue on Hispaniola, now modern-day Haiti, was a major port of disembarkation, and malaria is still actively transmitted there. We undertook a detailed study of the phylogenetics of the Haitian parasites and those from Colombia and Peru utilizing whole-genome sequencing. Principal-component and phylogenetic analyses, based upon single nucleotide polymorphisms (SNPs) in protein coding regions, indicate that, despite the potential for millions of introductions from Africa, the Haitian parasites share an ancestral relationship within a well-supported monophyletic clade with parasites from South America, while belonging to a distinct lineage. This result, in stark contrast to the historical record of parasite introductions, is best explained by a severe population bottleneck experienced by the parasites introduced into the Americas. Here, evidence is presented for targeted selection of rare African alleles in genes which are expressed in the mosquito stages of the parasite's life cycle. These genetic markers support the hypothesis that the severe population bottleneck was caused by the required adaptation of the parasite to transmission by new definitive hosts among the Anopheles (Nyssorhynchus) spp. found in the Caribbean and South America.IMPORTANCE Historical data suggest that millions of P. falciparum parasite lineages were introduced into the Americas during the trans-Atlantic slave trade, which would suggest a paraphyletic origin of the extant isolates in the Western Hemisphere. Our analyses of whole-genome variants show that the American parasites belong to a well-supported monophyletic clade. We hypothesize that the required adaptation to American vectors created a severe bottleneck, reducing the effective introduction to a few lineages. In support of this hypothesis, we discovered genes expressed in the mosquito stages of the life cycle that have alleles with multiple, high-frequency or fixed, nonsynonymous mutations in the American populations which are rarely found in African isolates. These alleles appear to be in gene products critical for transmission through the anopheline vector. Thus, these results may inform efforts to develop novel transmission-blocking vaccines by identifying parasite proteins functionally interacting with the vector that are important for successful transmission. Further, to the best of our knowledge, these are the first whole-genome data available from Haitian P. falciparum isolates. Defining the genome of these parasites provides genetic markers useful for mapping parasite populations and monitoring parasite movements/introductions.

Humanized Mice and the Rebirth of Malaria Genetic Crosses

The first experimental crosses carried out with the human malaria parasite Plasmodium falciparum played a key role in determining the genetic loci responsible for drug resistance, virulence, invasion, growth rate, and transmission. These crosses relied on splenectomized chimpanzees to complete the liver stage of the parasite's life cycle and the subsequent transition to asexual blood stage culture followed by cloning of recombinant progeny in vitro. Crosses can now be routinely carried out using human-liver-chimeric mice infused with human erythrocytes to generate hundreds of unique recombinant progeny for genetic linkage mapping, bulk segregant analysis, and high-throughput 'omics readouts. The high number of recombinant progeny should allow for unprecedented power and efficiency in the execution of a systems genetics approach to study P. falciparum biology.

Selective whole genome amplification of Plasmodium malariae DNA from clinical samples reveals insights into population structure

The genomic diversity of Plasmodium malariae malaria parasites is understudied, partly because infected individuals tend to present with low parasite densities, leading to difficulties in obtaining sufficient parasite DNA for genome analysis. Selective whole genome amplification (SWGA) increases the relative levels of pathogen DNA in a clinical sample, but has not been adapted for P. malariae parasites. Here we design customized SWGA primers which successfully amplify P. malariae DNA extracted directly from unprocessed clinical blood samples obtained from patients with P. malariae-mono-infections from six countries, and further test the efficacy of SWGA on mixed infections with other Plasmodium spp. SWGA enables the successful whole genome sequencing of samples with low parasite density (i.e. one sample with a parasitaemia of 0.0064% resulted in 44% of the genome covered by ≥ 5 reads), leading to an average 14-fold increase in genome coverage when compared to unamplified samples. We identify a total of 868,476 genome-wide SNPs, of which 194,709 are unique across 18 high-quality isolates. After exclusion of the hypervariable subtelomeric regions, a high-quality core subset of 29,899 unique SNPs is defined. Population genetic analysis suggests that P. malariae parasites display clear geographical separation by continent. Further, SWGA successfully amplifies genetic regions of interest such as orthologs of P. falciparum drug resistance-associated loci (Pfdhfr, Pfdhps, Pfcrt, Pfk13 and Pfmdr1), and several non-synonymous SNPs were detected in these genes. In conclusion, we have established a robust SWGA approach that can assist whole genome sequencing of P. malariae, and thereby facilitate the implementation of much-needed large-scale multi-population genomic studies of this neglected malaria parasite. As demonstrated in other Plasmodia, such genetic diversity studies can provide insights into the biology underlying the disease and inform malaria surveillance and control measures.

An exported kinase family mediates species-specific erythrocyte remodelling and virulence in human malaria

The most severe form of human malaria is caused by Plasmodium falciparum. Its virulence is closely linked to the increase in rigidity of infected erythrocytes and their adhesion to endothelial receptors, obstructing blood flow to vital organs. Unlike other human-infecting Plasmodium species, P. falciparum exports a family of 18 FIKK serine/threonine kinases into the host cell, suggesting that phosphorylation may modulate erythrocyte modifications. We reveal substantial species-specific phosphorylation of erythrocyte proteins by P. falciparum but not by Plasmodium knowlesi, which does not export FIKK kinases. By conditionally deleting all FIKK kinases combined with large-scale quantitative phosphoproteomics we identified unique phosphorylation fingerprints for each kinase, including phosphosites on parasite virulence factors and host erythrocyte proteins. Despite their non-overlapping target sites, a network analysis revealed that some FIKKs may act in the same pathways. Only the deletion of the non-exported kinase FIKK8 resulted in reduced parasite growth, suggesting the exported FIKKs may instead support functions important for survival in the host. We show that one kinase, FIKK4.1, mediates both rigidification of the erythrocyte cytoskeleton and trafficking of the adhesin and key virulence factor PfEMP1 to the host cell surface. This establishes the FIKK family as important drivers of parasite evolution and malaria pathology.

The RH5-CyRPA-Ripr Complex as a Malaria Vaccine Target

Despite ongoing efforts, a highly effective vaccine against Plasmodium falciparum remains elusive. Vaccines targeting the pre-erythrocytic stages of the P. falciparum life cycle are the most advanced to date, affording moderate levels of efficacy in field trials. However, the discovery that the members of the merozoite PfRH5-PfCyRPA-PfRipr (RCR) complex are capable of inducing strain-transcendent neutralizing antibodies has renewed enthusiasm for the possibility of preventing disease by targeting the parasite during the blood stage of infection. With Phase I/II clinical trials now underway using first-generation vaccines against PfRH5, and more on the horizon for PfCyRPA and PfRipr, this review explores the rationale and future potential of the RCR complex as a P. falciparum vaccine target.

Evolutionary and population (epi)genetics of immunity to infection

Immune response is one of the functions that have been more strongly targeted by natural selection during human evolution. The evolutionary genetic dissection of the immune system has greatly helped to distinguish genes and functions that are essential, redundant or advantageous for human survival. It is also becoming increasingly clear that admixture between early Eurasians with now-extinct hominins such as Neanderthals or Denisovans, or admixture between modern human populations, can be beneficial for human adaptation to pathogen pressures. In this review, we discuss how the integration of population genetics with functional genomics in diverse human populations can inform about the changes in immune functions related to major lifestyle transitions (e.g., from hunting and gathering to farming), the action of natural selection to the evolution of the immune system, and the history of past epidemics. We also highlight the need of expanding the characterization of the immune system to a larger array of human populations-particularly neglected human groups historically exposed to different pathogen pressures-to fully capture the relative contribution of genetic, epigenetic, and environmental factors to immune response variation in humans.

Optimization of parasite DNA enrichment approaches to generate whole genome sequencing data for Plasmodium falciparum from low parasitaemia samples

Background

Owing to the large amount of host DNA in clinical samples, generation of high-quality Plasmodium falciparum whole genome sequencing (WGS) data requires enrichment for parasite DNA. Enrichment is often achieved by leukocyte depletion of infected blood prior to storage. However, leukocyte depletion is difficult in low-resource settings and limits analysis to prospectively-collected samples. As a result, approaches such as selective whole genome amplification (sWGA) are being used to enrich for parasite DNA. However, sWGA has had limited success in generating reliable sequencing data from low parasitaemia samples. In this study, enzymatic digestion with MspJI prior to sWGA and whole genome sequencing was evaluated to determine whether this approach improved genome coverage compared to sWGA alone. The potential of sWGA to cause amplification bias in polyclonal infections was also examined.

Methods

DNA extracted from laboratory-created dried blood spots was treated with a modification-dependent restriction endonuclease, MspJI, and filtered via vacuum filtration. Samples were then selectively amplified using a previously reported sWGA protocol and subjected to WGS. Genome coverage statistics were compared between the optimized sWGA approach and the previously reported sWGA approach performed in parallel. Differential amplification by sWGA was assessed by comparing WGS data generated from lab-created mixtures of parasite isolates, from the same geographical region, generated with or without sWGA.

Results

MspJI digestion did not enrich for parasite DNA. Samples that underwent vacuum filtration (without MspJI digestion) prior to sWGA had the highest parasite DNA concentration and displayed greater genome coverage compared to MspJI + sWGA and sWGA alone, particularly for low parasitaemia samples. The optimized sWGA (filtration + sWGA) approach was successfully used to generate WGS data from 218 non-leukocyte depleted field samples from Malawi. Sequences from lab-created mixtures of parasites did not show evidence of differential amplification of parasite strains compared to directly sequenced samples.

Conclusion

This optimized sWGA approach is a reliable method to obtain WGS data from non-leukocyte depleted, low parasitaemia samples. The absence of amplification bias in data generated from mixtures of isolates from the same geographic region suggests that this approach can be appropriately used for molecular epidemiological studies.

Ancient Introgression between Two Ape Malaria Parasite Species

The Laverania clade comprises the human malaria parasite Plasmodium falciparum as well as at least seven additional parasite species that infect wild African apes. A recent analysis of Laverania genome sequences (Otto TD, et al. 2018. Genomes of all known members of a Plasmodium subgenus reveal paths to virulent human malaria. Nat Microbiol. 3: 687–697) reported three instances of interspecies gene transfer, one of which had previously been described. Generating gene sequences from additional ape parasites and re-examining sequencing reads generated in the Otto et al. study, we identified one of the newly described gene transfers as an assembly artifact of sequences derived from a sample coinfected by two parasite species. The second gene transfer between ancestors of two divergent chimpanzee parasite lineages was confirmed, but involved a much larger number of genes than originally described, many of which encode exported proteins that remodel, or bind to, erythrocytes. Because successful hybridization between Laverania species is very rare, it will be important to determine to what extent these gene transfers have shaped their host interactions.

Extracellular Vesicles Could Carry an Evolutionary Footprint in Interkingdom Communication

Extracellular vesicles (EVs) are minute particles secreted by the cells of living organisms. Although the functional role of EVs is not yet clear, recent work has highlighted their role in intercellular communication. Here, we expand on this view by suggesting that EVs can also mediate communication among interacting organisms such as hosts, pathogens and vectors. This inter-kingdom communication via EVs is likely to have important evolutionary consequences ranging from adaptation of parasites to specialized niches in the host, to host resistance and evolution and maintenance of parasite virulence and transmissibility. A potential system to explore these consequences is the interaction among the human host, the mosquito vector and Plasmodium parasite involved in the malaria disease. Indeed, recent studies have found that EVs derived from Plasmodium infected red blood cells in humans are likely mediating the parasite's transition from the asexual to sexual stage, which might facilitate transmission to the mosquito vector. However, more work is needed to establish the adaptive consequences of this EV signaling among different taxa. We suggest that an integrative molecular approach, including a comparative phylogenetic analysis of the molecules (e.g., proteins and nucleic acids) derived from the EVs of interacting organisms (and their closely-related species) in the malaria system will prove useful for understanding interkingdom communication. Such analyses will also shed light on the evolution and persistence of host, parasite and vector interactions, with implications for the control of vector borne infectious diseases.

Falciparum malaria from coastal Tanzania and Zanzibar remains highly connected despite effective control efforts on the archipelago

BACKGROUND

Tanzania's Zanzibar archipelago has made significant gains in malaria control over the last decade and is a target for malaria elimination. Despite consistent implementation of effective tools since 2002, elimination has not been achieved. Importation of parasites from outside of the archipelago is thought to be an important cause of malaria's persistence, but this paradigm has not been studied using modern genetic tools.

METHODS

Whole-genome sequencing (WGS) was used to investigate the impact of importation, employing population genetic analyses of Plasmodium falciparum isolates from both the archipelago and mainland Tanzania. Ancestry, levels of genetic diversity and differentiation, patterns of relatedness, and patterns of selection between these two populations were assessed by leveraging recent advances in deconvolution of genomes from polyclonal malaria infections.

RESULTS

Significant decreases in the effective population sizes were inferred in both populations that coincide with a period of decreasing malaria transmission in Tanzania. Identity by descent analysis showed that parasites in the two populations shared long segments of their genomes, on the order of 5 cM, suggesting shared ancestry within the last 10 generations. Even with limited sampling, two of isolates between the mainland and Zanzibar were identified that are related at the expected level of half-siblings, consistent with recent importation.

CONCLUSIONS

These findings suggest that importation plays an important role for malaria incidence on Zanzibar and demonstrate the value of genomic approaches for identifying corridors of parasite movement to the island.

Proteomic Atomics Reveals a Distinctive Uracil-5-Methyltransferase

Carbon (C), hydrogen (H), nitrogen (N), oxygen (O), and sulfur (S) atoms intrigue as they are the foundation for amino acid (AA) composition and the folding and functions of proteins and thus define and control the survival of a cell, the smallest unit of life. Here, we calculated the proteomic atom distribution in >1500 randomly selected species across the entire current phylogenetic tree and identified uracil-5-methyltransferase (U5MTase) of the protozoan parasite Plasmodium falciparum (Pf, strain Pf3D7), with a distinct atom and AA distribution pattern. We determined its apicoplast location and in silico 3D protein structure to refocus attention exclusively on U5MTase with tremendous potential for therapeutic intervention in malaria. Around 300 million clinical cases of malaria occur each year in tropical and subtropical regions of the world, resulting in over one million deaths annually, placing malaria among the most serious infectious diseases. Genomic and proteomic research of the clades of parasites containing Pf is progressing slowly and the functions of most of the ∼5300 genes are still unknown. We applied a 'bottom-up' comparative proteomic atomics analysis across the phylogenetic tree to visualize a protein molecule on its actual basis - i. e., its atomic level. We identified a protruding Pf3D7-specific U5MTase, determined its 3D protein structure, and identified potential inhibitory drug molecules through in silico drug screening that might serve as possible remedies for the treatment of malaria. Besides, this atomic-based proteome map provides a unique approach for the identification of parasite-specific proteins that could be considered as novel therapeutic targets.

Genetic Evidence of Focal Plasmodium falciparum Transmission in a Pre-elimination Setting in Southern Province, Zambia

BACKGROUND

Southern Province, Zambia has experienced a dramatic decline in Plasmodium falciparum malaria transmission in the past decade and is targeted for elimination. Zambia's National Malaria Elimination Program recommends reactive case detection (RCD) within 140 m of index households to enhance surveillance and eliminate remaining transmission foci.

METHODS

To evaluate whether RCD captures local transmission, we genotyped 26 microsatellites from 106 samples collected from index (n = 27) and secondary (n = 79) cases detected through RCD in the Macha Hospital catchment area between January 2015 and April 2016.

RESULTS

Participants from the same RCD event harbored more genetically related parasites than those from different RCD events, suggesting that RCD captures, at least in part, infections related through local transmission. Related parasites clustered in space and time, up to at least 250 m from index households. Spatial analysis identified a putative focal transmission hotspot.

CONCLUSIONS

The current RCD strategy detects focal transmission events, although programmatic guidelines to screen within 140 m of index households may fail to capture all secondary cases. This study highlights the utility of parasite genetic data in assessing programmatic interventions, and similar approaches may be useful to malaria elimination programs seeking to tailor intervention strategies to the underlying transmission epidemiology.

VAR2CSA binding phenotype has ancient origin and arose before Plasmodium falciparum crossed to humans: implications in placental malaria vaccine design

VAR2CSA is a leading candidate for developing a placental malaria (PM) vaccine that would protect pregnant women living in malaria endemic areas against placental infections and improve birth outcomes. Two VAR2CSA-based PM vaccines are currently under clinical trials, but it is still unclear if the use of a single VAR2CSA variant will be sufficient to induce a broad enough humoral response in humans to cross-react with genetically diverse parasite populations. Additional immuno-focusing vaccine strategies may therefore be required to identify functionally conserved antibody epitopes in VAR2CSA. We explored the possibility that conserved epitopes could exist between VAR2CSA from the chimpanzee parasite Plasmodium reichenowi and Plasmodium falciparum sequences. Making use of VAR2CSA recombinant proteins originating from both species, we showed that VAR2CSA from P. reichenowi (Pr-VAR2CSA) binds to the placental receptor CSA with high specificity and affinity. Antibodies raised against Pr-VAR2CSA were able to recognize native VAR2CSA from different P. falciparum genotypes and to inhibit the interaction between CSA and P. falciparum-infected erythrocytes expressing different VAR2CSA variants. Our work revealed the existence of cross-species inhibitory epitopes in VAR2CSA and calls for pre-clinical studies assessing the efficacy of novel VAR2CSA-based cross-species boosting regimens.

Plasmodium infection in Kerteszia cruzii (Diptera: Culicidae) in the Atlantic tropical rain forest, southeastern Brazil

In Southeastern Brazil, Kerteszia cruzii (former Anopheles cruzii), a bromeliad mosquito species, is considered an efficient human Plasmodium spp. vector. In this region, recent studies showed asymptomatic or sub-patent Plasmodium falciparum infection. In areas of the Atlantic coast in Rio de Janeiro, Plasmodium simium infection was recently reported in both human and howler monkey. Considering that (1) few malaria cases are reported each year in areas across the tropical Atlantic rain forest in southeastern Brazil; (2) malaria elimination in Atlantic forest is challenged by circulation of P. falciparum and P. simium in humans; (3) the complexity of malaria epidemiology in this region; and (4) the public health importance of Kerteszia cruzii as a sylvatic vector; the major goal of this study is to evaluate Plasmodium infection in Ke. cruzii. Mosquito sampling collections were conducted in Esteiro do Morro and Sítio Itapuan, in Cananeia municipality, and Tapiraí municipality in Ribeira Valley, southeastern São Paulo state, Brazil. Influence of climate and landscape factors in Plasmodium infection in Ke. cruzii was addressed. Among the 1719 mosquitoes tested, 3 females collected in Sítio Itapuan and three from Tapiraí were found infected with either P. vivax or P. simium. Results of statistical analyses did not demonstrate association between Plasmodium infection in mosquito and the landscape. Mosquito infection was found in two landscape clusters, with Plasmodium detected in forest fringe mosquitoes. This finding shows that Ke. cruzii can facilitate transmission among human and non-human primates. Plasmodium falciparum was not identified in the samples analyzed. Spatiotemporal variation in local malaria incidence, low prevalence of Plasmodium, variations in humidity and temperature can explain the absence of mosquitoes infected with P. falciparum in the study.

Selective whole genome amplification and sequencing of Coxiella burnetii directly from environmental samples

Whole genome sequencing (WGS) is a widely available, inexpensive means of providing a wealth of information about an organism's diversity and evolution. However, WGS for many pathogenic bacteria remain limited because they are difficult, slow and/or dangerous to culture. To avoid culturing, metagenomic sequencing can be performed directly on samples, but the sequencing effort required to characterize low frequency organisms can be expensive. Recently developed methods for selective whole genome amplification (SWGA) can enrich target DNA to provide efficient sequencing. We amplified Coxiella burnetii (a bacterial select agent and human/livestock pathogen) from 3 three environmental samples that were overwhelmed with host DNA. The 68- to 147-fold enrichment of the bacterial sequences provided enough genome coverage for SNP analyses and phylogenetic placement. SWGA is a valuable tool for the study of difficult-to-culture organisms and has the potential to facilitate high-throughput population characterizations as well as targeted epidemiological or forensic investigations.

Resurrection of the ancestral RH5 invasion ligand provides a molecular explanation for the origin of P. falciparum malaria in humans

Many important infectious diseases are the result of zoonoses, in which pathogens that normally infect animals acquire mutations that enable the breaching of species barriers to permit the infection of humans. Our understanding of the molecular events that enable host switching are often limited, and yet this is a fundamentally important question. Plasmodium falciparum, the etiological agent of severe human malaria, evolved following a zoonotic transfer of parasites from gorillas. One gene—rh5—which encodes an essential ligand for the invasion of host erythrocytes, is suspected to have played a critical role in this host switch. Genome comparisons revealed an introgressed sequence in the ancestor of P. falciparum containing rh5, which likely allowed the ancestral parasites to infect both gorilla and human erythrocytes. To test this hypothesis, we resurrected the ancestral introgressed reticulocyte-binding protein homologue 5 (RH5) sequence and used quantitative protein interaction assays to demonstrate that this ancestral protein could bind the basigin receptor from both humans and gorillas. We also showed that this promiscuous receptor binding phenotype of RH5 was shared with the parasite clade that transferred its genome segment to the ancestor of P. falciparum, while the other lineages exhibit host-specific receptor binding, confirming the central importance of this introgression event for Plasmodium host switching. Finally, since its transfer to humans, P. falciparum, and also the RH5 ligand, have evolved a strong human specificity. We show that this subsequent restriction to humans can be attributed to a single amino acid mutation in the RH5 sequence. Our findings reveal a molecular pathway for the origin and evolution of human P. falciparum malaria and may inform molecular surveillance to predict future zoonoses.

This study reveals a molecular pathway by which Plasmodium falciparum malaria arose via zoonotic transfer from gorillas by comparing the host receptor binding properties of extant and “resurrected” ancestral versions of the parasite’s erythrocyte invasion ligand RH5.

Genetic mapping of fitness determinants across the malaria parasite Plasmodium falciparum life cycle

Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance.