Xin-zhuan Su
Cui Zhang
Deirdre A. Joy2*- 1Laboratory of Malaria and Vector Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD, United States
- 2Parasitology and International Programs Branch, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD, United States
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.
Introduction
Malaria Parasites, Life Cycle, and Genome
Malaria is a parasitic disease caused by Plasmodium species (spp.), unicellular protozoan organisms in the phylum of Apicomplexa. The species that infect humans include Plasmodium falciparum, P. vivax, P. malariae, P. knowlesi, and P. ovale, with P. ovale recently recognized as two sub-species called Plasmodium ovale curtisi (classic type) and Plasmodium ovale wallikeri (Sutherland et al., 2010). Whereas, P. vivax is the most widespread species, P. falciparum is the most deadly to humans. There are several other species of malaria parasites that infect non-human primates including multiple species that infect African apes, and others that infect Asian and New World monkeys (Escalante et al., 2005; Liu et al., 2010; Boundenga et al., 2015; Martinelli and Culleton, 2018; Galinski, 2019). There are also malaria parasites of rodents such as P. berghei, P. chabaudi, and P. yoelii, often used as animal models of malaria, as well as parasites that infect birds (P. gallinaceum and P. relictum) (Duval and Ariey, 2012; Ramiro et al., 2012; Fecchio et al., 2020), lizards, bats, and ungulates (Schall, 1982; Schaer et al., 2013; Templeton et al., 2016; Galen et al., 2018). Human malaria cases are mostly reported in tropical and subtropical regions, particularly in Africa, south and southeast Asia, and south and central America. In 2018, an estimated 228 million cases of malaria occurred worldwide, with ~213 million (or 93%) in Africa (WHO, 2019).
Human malaria infections start with a bite from an infected female Anopheles mosquito, injecting sporozoites into the skin of a host (Figure 1). The sporozoites move to the liver through blood vessels and then traverse liver sinusoidal endothelial cells or Kupffer cells to infect hepatocytes (Baer et al., 2007; Tavares et al., 2013). Inside a hepatocyte, a parasite divides into thousands of merozoites through a process called schizogony. Mature merozoites enter the bloodstream after rupture of the infected hepatocyte and invade erythrocytes, starting a new cycles of schizogony within red blood cells (RBCs) which includes asexual replication of their haploid genome. Within RBCs parasites develop through ring, trophozoite, and schizont stages. The resulting mature schizont is segmented in appearance and contains 16–32 daughter merozoites. The infected RBC (iRBC) ruptures, releasing the daughter sporozoites to invade new RBCs. This cycle of erythrocytic development takes approximately 24 h for P. knowlesi, 48 h for P. falciparum, P. vivax, and P. ovale, and 72 h for P. malariae. In addition to merozoites, ruptured iRBCs also release various by-products of parasite metabolism such as hemozoin formed when the parasite digests hemoglobin. The release of parasite materials also triggers host responses leading to various clinical symptoms, including fever, chills, headache, dizziness, back pain, and myalgia (Ashley et al., 2018). Some patients may progress to severe malaria with coma (cerebral malaria), pulmonary edema, acute renal failure, jaundice, severe anemia, acidosis, hypoglycemia, and death (Trampuz et al., 2003; Sypniewska et al., 2017; Luzolo and Ngoyi, 2019). In response to changes in the host environment such as variation in metabolites and immune factors (Josling and Llinas, 2015; Brancucci et al., 2017), some of the parasites differentiate into sexual stages termed male and female gametocytes. When a mosquito takes a bloodmeal from an infected host, gametocytes will undergo sexual development in the midgut of the mosquito. Within minutes of entering the mosquito midgut, male and female gametocytes differentiate into male and female gametes that fertilize to form zygotes and then into motile ookinetes. Ookinetes penetrate through the mosquito midgut wall to develop into oocysts containing thousands of sporozoites. Mature sporozoites migrate to salivary glands and are injected into a new vertebrate host when the mosquito bites again, completing the life cycle.
Figure 1. Life cycle of malaria parasites. Malaria parasite infection starts with a mosquito bite injecting sporozoites into human body where the parasites travel to liver and invade hepatocytes. After replication in liver cells, thousands of mature merozoites are released into blood stream and invade red blood cells (RBCs). The parasites then replicate in RBCs that subsequently rupture, releasing more merozoites to invade new RBCs for another cycle. Malaria symptoms correlate with the rupture of infected RBCs. A small number of merozoites develop into male and female gametocytes after invading new RBCs. Gametocytes differentiate into male and female gametes in the mosquito midgut when another mosquito takes a blood meal. Male and female gametes fertilize to produce zygotes that develop into motile ookinetes and oocysts after penetrating the mosquito midgut. A large number of sporozoites from each oocysts moves to salivary glands and are injected into a new human host when the mosquito bites again.
The genomes of many malaria parasites infecting humans (P. falciparum, P. vivax, P. ovale, P. malariae, P. knowlesi), non-human primates (P. reichenowi, P. praefalciparum, P. blacklocki, P. adleri, P. billcollinsi, P. gaboni, P. cynomolgi, P. coatneyi, P. inui, and P. fragile), rodents (P. berghei, P. yoelii, P. chabaudi, and P. vinkei) and birds (P. relictum and P. gallinaceum) have been sequenced and analyzed (Shutler et al., 2005; Pain et al., 2008; Otto et al., 2014a,b, 2018; Lauron et al., 2015; Ansari et al., 2016; Auburn et al., 2016; Sundararaman et al., 2016; Pasini et al., 2017; Rutledge et al., 2017; Bohme et al., 2018; Su et al., 2019). Malaria parasites have relatively small haploid genomes, ranging from 20 to 35 megabases (Mb) that contain 14 chromosomes, a circular plastid genome of ~35 kb, and multiple copies of a 6-kb mitochondrial DNA (Gardner et al., 2002; Otto et al., 2014a,b; Auburn et al., 2016; Pasini et al., 2017; Bohme et al., 2018). The homologous genes in different species of malaria parasites are often found in syntenic blocks arranged in different order on the chromosomes (Carlton et al., 2005; Kooij et al., 2005). For more information on Plasmodium genome and genomics associated with parasite development, diseases, diagnosis, epidemiology, and vaccine development please consult (Kirchner et al., 2016; Bourgard et al., 2018; Escalante and Pacheco, 2019; Galinski, 2019; Garrido-Cardenas et al., 2019; Su et al., 2019; Videvall, 2019).
Evolution and Adaptation of Malaria Parasites in Great Apes
Origin of P. falciparum
Our understanding of the host-parasite interactions involved in the emergence of P. falciparum in humans has been greatly enabled by the discovery of multiple closely related malaria species in African apes beginning in the early 2000s. Three malaria species morphologically similar to P. falciparum were described in African apes during the last century (Reichenow, 1920; Blacklock and Adler, 1922; Adler, 1923), and of these only one, P. reichenowi, had been sequenced (Neafsey et al., 2005). This picture began to change with the discovery of a new malaria species identified in chimpanzees kept as pets in a village in Gabon (Ollomo et al., 2009). Based on the complete mitochondrial genome this parasite was found to fall within the Laverania subgenus that includes P. reichenowi and P. falciparum. A number of other studies quickly followed, leading to the discovery of additional novel malaria species, and a rethinking of the origin of P. falciparum in humans. DNA sequences of the cytochrome oxidase b (cyt b) gene in samples from a combination of both wild-living and captive chimpanzees and gorillas showed an unexpectedly diverse collection of Plasmodium lineages circulating in African apes, suggesting that P. falciparum likely originated from P. reichenowi via a single transfer event from chimpanzees to humans (Prugnolle et al., 2010). Sequencing of the Plasmodium mitochondrial genome and two nuclear genes, dihydrofolate reductase-thymidylate synthase (dhfr-ts) and the gene encoding the merozoite surface protein 2 (msp2), from blood samples of captive chimpanzees and bonobos in Uganda and the Democratic Republic of Congo identified two new Plasmodium species, and led the authors to conclude that P. falciparum arose in bonobos (Pan paniscus) (Krief et al., 2010). Yet another study of a small number of chimpanzees and gorillas from Cameroon previously kept as pets targeted two mitochondrial genes (cytb and cox1), one plastid gene (tufA), and one nuclear gene (ldh), and purportedly found P. falciparum in both chimpanzee and gorilla blood samples (Duval et al., 2010).
The contradictory findings on the origin of P. falciparum could be explained in part by sampling limitations, including small sample sizes and a general reliance on captive animals which introduces the possibility of contemporary infection by humans. These concerns were addressed in a large-scale study of non-invasive samples from wild-living chimpanzees, gorillas, and bonobos throughout central Africa (Liu et al., 2010). Mitochondrial, apicoplast, and nuclear genes from over 1,000 samples were sequenced, and most infected samples were found to contain multiple genetically diverse parasites, highlighting another potential problem with the previous studies. In samples with multiple genetically diverse parasites, conventional PCR methods risk generating recombinant artifacts that could mislead the data. To guard against this, samples were subjected to limiting dilution such that each positive PCR reaction contained no more than a single parasite. Using this approach, the number of host-specific parasite species within the Laverania subgenus was increased to seven: P. gaboni, P. billcollinsi, and P reichenowi in chimpanzees; P. praefalciparum, P. blackloci, and P. adleri in gorillas; and P. falciparum in humans (Liu et al., 2010). One additional parasite species, Plasmodium lomamiensis, was later identified in bonobos (Liu et al., 2017). P. praefalciparum in western gorillas was found to be genetically nearly identical to P. falciparum, which formed a monophyletic clade (Figure 2), strongly suggesting that the P. falciparum derived from a host switch from gorillas to humans. Genome sequences from P. reichenowi and two P. gaboni parasites would later allow for an estimation of the genetic diversity of these two species (Sundararaman et al., 2016). Both P. reichenowi and P. gaboni were estimated to be at least 10X more diverse than P. falciparum, strongly suggesting a relatively recent origin for P. falciparum. This was confirmed when additional Laverania genomes were used to estimate the divergence of P. falciparum and P. praefalciparum to be 40,000–60,000 years ago (Otto et al., 2018). Additionally, the successful jump to humans by P. falciparum appears to have involved multiple parasites over time rather than a single event (Otto et al., 2018). The relationships of Laverania species based on maximum likelihood analysis of a 3.4 kb region of the mitochondrial genome are presented in Figure 2.
Figure 2. Relationship among the Laverania species, adapted from Figure 1A in Plenderleith et al. (2019). Icons indicate the parasites usual host species: gorillas (Gorilla gorilla), chimpanzees (Pan troglodytes), bonobos (Pan paniscus), and humans. The tree was derived by maximum likelihood analysis of a 3.4 kb region of the mitochondrial genome (see Supplementary Figure 1A in Liu et al., 2010).
Molecular Basis of Host Trophism and Adaptation to Humans
The Laverania exhibit strong host trophism, with each parasite species infecting a single host species, although rare host switching events have occurred. Host trophism is not seen to the same extent in malaria parasites that infect Asian primates (P. cynomolgi and P. knowlesi). Barriers restricting parasite species to a single host species could exist at any of the host-parasite-vector interfaces, although considerable attention has focused on host-parasite interactions during the blood stage, specifically erythrocyte invasion (Martin et al., 2005; Rich et al., 2009; Rayner et al., 2011; Wanaguru et al., 2013; Galaway et al., 2019; Plenderleith et al., 2019; Proto et al., 2019).
The production of a greatly improved P. reichenowi reference genome, as well as a partial genome of P. gaboni, allowed for the use of comparative genomics to look for genetic loci involved in host restriction. The genomes of P. reichenowi, P. gaboni, and P. falciparum showed good conservation and synteny, but with striking differences in some erythrocyte invasion loci (Otto et al., 2014b). A large-scale genomics study of Laverania including multiple genomes from all known species found evidence for three interspecific gene transfers (Otto et al., 2018), one which would later be challenged by a reanalysis of the sequence reads (Plenderleith et al., 2019). Several genes involved in either invasion or pathogenesis were fixed in parasites infecting the different primate hosts (Otto et al., 2018). Both a gene transfer event (Sundararaman et al., 2016) and an ancient introgression involving some erythrocyte invasion and surface protein genes (Otto et al., 2018) may have laid the groundwork for P. falciparum's host switch to humans. There is strong evidence of a horizontal transfer of an 8 kb region of chromosome 4 from P. adleri to the ancestor of P. falciparum, although disagreements exist as to whether the 8 kb fragment containing Rh5 resulted from introgression or a horizontal transfer event (Sundararaman et al., 2016; Galaway et al., 2019). This region contained two erythrocyte invasion genes encoding the reticulocyte-binding-like homologous protein 5 (RH5) and the cysteine-rich protective antigen (CyRPA). RH5 differs from the other members of the RH family in that it is a secreted rather than a membrane-tethered protein and cannot be genetically deleted in any P. falciparum. A combination of ancestral sequence reconstruction and quantitative protein interaction assays was used to investigate the role of RH5 in the emergence of P. falciparum in humans (Galaway et al., 2019). Probabilistic-based approaches to deduce the likely sequence of the rh5 gene in the common ancestor of P. falciparum and P. praefalicarum converged on a single sequence, which the authors termed the ancestral introgressed Rh5 gene (IntRh5). P. falciparum RH5 and its orthologs were expressed as soluble recombinant proteins, and the binding affinities of RH5 orthologs for human, gorilla, and chimpanzee basigins were quantified using surface plasmon resonance (SPR) and confirmed with cell-based binding assays. Surprisingly, IntRH5 protein was found to bind to human and gorilla basigin with similar affinities. This promiscuous receptor-binding phenotype was also observed in all parasites belonging to clade A (P. adleri and P. gaboni) but not in parasites belonging to clade B (P. billcollinsi, P. reichenowi, P. preafalciparum). As P. falciparum is restricted to humans with no African apes having been found to be naturally infected with this parasite, the question remains as to how this restriction occurred if the ancestral character state was permissive. Others have shown that recombinant PfRH5 does not interact with gorilla basigin, and only weakly interacted with chimpanzee basigin (Wanaguru et al., 2013). Galaway et al. (2019) identified six mutations separating IntRh5 and PfRh5. A single mutation at residue 200 resulted in the complete loss of PfRH5 binding to gorilla basigin, suggesting that the host switch from gorillas to humans involved at least two steps, one in which erythrocytes from both hosts could be infected, followed by restriction in humans (Otto et al., 2018).
The erythrocyte binding-like (EBL) gene family also plays a critical role in erythrocyte invasion through the recognition of glycophorins, the predominant sialoglycoprotien on the erythrocyte surface. Different EBLs recognize different gylcophorins; for example, Pf EBA175 recognizes glycophorin A, and invasion of human erythrocytes by P. falciparum is heavily reliant on this EBL. The potential for sialic acid recognition to play a role in the exclusion of P. falciparum from chimpanzees and P. reichenowi from humans was recognized and investigated prior to the discovery of multiple Laverania species in African apes (Martin et al., 2005). The P. falciparum EBL165 gene (PfEbl165) is a pseudogene due to a frame-shift mutation in all known P. falciparum isolates studied so far (Proto et al., 2019). In contrast, all Ebl165 homologs in ape Laverania species encode an intact protein (Otto et al., 2014b, 2018; Sundararaman et al., 2016). Due to a mutation in the gene encoding cytidine monophospho-N-acetylneuraminic acid hydroxylase (CMAH), human erythrocytes produce only Neu5Ac sialic acids whereas ape erythrocytes produce both Neu5Ac and Neu5Gc. EBA165 binding is dependent on the specific sialic acid repertoire displayed, and the absence of ape-specific Neu5Gc sialic acid on human RBCs is sufficient to prevent PfEBL165 from binding to them (Proto et al., 2019). CRISPR-Cas9 mediated correction of the frameshift mutation in PfEbl165 led to the epigenetic silencing of a small number of genes in two regions of chromosome 4 and 11 that included PfEbl165 itself and PfRh4 (another invasion gene) (Proto et al., 2019). It appears as though a functional PfEBL165 is not compatible with robust growth of P. falciparum in humans, suggesting that the inactivation of this gene may have contributed to the parasites' establishment in this host. As an increasing number of candidate loci are studied, it is becoming clear that a single molecular event is not sufficient to explain the emergence of P. falciparum as a human pathogen. Accumulating evidence suggests that RH5, EBL165, other members of the RH5 complex, invasion proteins, and possibly exported proteins that remodel the infected erythrocyte surface all contribute to host specificity and tropism. Gene transfers between Plasmodium species have played an important role in the parasite evolution and host-parasite interactions (Plenderleith et al., 2019).
Host-Parasite Interaction and Evolution
Highly Polymorphic Immune Targets
Malaria parasites trigger an immune response the moment when they enter a host. To survive in this hostile environment, the parasite displays a range of strategies to evade host killing mechanisms, including variations in antigen epitopes targeted by host immune machinery and interference or suppression of specific arms of the host immune response (Recker et al., 2011; Yam and Preiser, 2017; Larsen et al., 2018; Xia et al., 2018). One of the consequences of these host-parasite interactions is increased genetic diversity at genes encoding proteins under host immune selection (Figure 3A), leading to genetic signatures of diversifying selection in the parasite genome. Some highly polymorphic genes include those encoding proteins such as the apical membrane antigen 1 (AMA1), merozoite surface protein 1 (MSP1), and circumsporozoite protein (CSP). These highly polymorphic genes as well as genome-wide polymorphisms such as single nucleotide polymorphisms (SNPs) and microsatellites (MSs) have been used for genotyping parasite strains, tracking parasite migration and disease outbreak, studying parasite molecular evolution, and evaluating host immune response and vaccine efficacy (Ranjit and Sharma, 1999; Volkman et al., 2001, 2007; Mu et al., 2002; Wootton et al., 2002; Joy et al., 2003; Jeffares et al., 2007; Ouattara et al., 2010; Drew et al., 2012; Tanabe et al., 2015; Dewasurendra et al., 2018; Stone et al., 2018). The high levels of genetic polymorphisms in immune targeted genes also provide signatures to identify potential vaccine candidates, particularly signatures of balancing selection within parasite populations (Mu et al., 2007; Tetteh et al., 2009; Ochola et al., 2010; Weedall and Conway, 2010; Amambua-Ngwa et al., 2012; Conway, 2015). For example, Tetteh et al. (2009) compared methods for prospectively identifying genes under balancing selection using 26 genes known or predicted to encode merozoite surface-exposed proteins and concluded that significant evidence of balancing selection could be detected using Hudson-Kreitman-Aguade (HKA) and Tajima's D (TjD) tests. In another study, high throughput sequencing of 65 clinical isolates from an endemic Gambian population identified 337 genes with at least three SNPs that had Tajima's D values > 0, including 25 genes with Tajima's D values > 1.5 (Amambua-Ngwa et al., 2012). However, the use of highly polymorphic antigens as vaccines for protection against malaria is challenging because antibody responses to these molecules are generally allele-specific or strain-specific (Graves et al., 2016; Laurens et al., 2017), which represents one of the major roadblocks for developing an effective malaria vaccine.
Figure 3. Potential impacts of host-parasite interactions and human intervention measures on parasite populations and evolution. (A) Host immunity will pressure parasites to diversify and/or express different sets of variant antigen genes, increasing diversity at immune targets. (B) A partially protective vaccine may selectively remove a specific parasite strain from the population. (C) Drug treatment will reduce population diversity by removing the drug sensitive parasites and selecting for one or a few drug resistant parasite clones that may then spread to many endemic regions. (D) Parasites passing through mosquitoes may generate new progeny through genetic recombination and chromosomal reassortment, reset gene expression profile particularly those related to parasite survival, and select for clones that are better adapted to specific mosquito species.
Diversity and Evolution of Gene Families of Malaria Parasites
The var Genes
In addition to highly polymorphic antigen genes, there are many polymorphic gene families in the Plasmodium genomes such as the P. falciparum var genes, SICA (Schizont Infected Cell Agglutination) var (SICAvar) genes in P. knowlesi, and the Plasmodium interspersed repeat (pir) multigene family, which include P. falciparum repetitive interspersed (rif) and subtelomeric variant open reading frame (stevor), found in all Plasmodium species (Su et al., 1995; Cunningham et al., 2010; Jemmely et al., 2010; Wahlgren et al., 2017; Galinski et al., 2018). Indeed, large multigene families are present in many Apicomplexa parasites (Reid, 2015). These gene families evolve at high rates and play critical roles in antigenic variation and immune evasion. The var genes encode P. falciparum erythrocyte membrane protein 1 (PfEMP1) proteins (Howard et al., 1988; Baruch et al., 1995; Su et al., 1995). Switches in expression among an estimated 60 var genes in every single parasite have been shown to correlate with variation in antigenic determinants that bind to different host receptors such as ICAM-1, CD36, EPCR, and other receptors (Smith et al., 1995; Wahlgren et al., 2017; Jensen et al., 2020). The binding of PfEMP1 proteins to host receptors is responsible for cytoadherence of infected erythrocytes (iRBCs) and the pathogenicity of severe malaria, particularly cerebral malaria (Smith et al., 2013; Jensen et al., 2020). The P. falciparum parasite evades host immunity via mutually exclusive expression (e.g., only one PfEMP1 is expressed on the surface of an iRBC at a time) of the highly diverse var family representing an almost unlimited gene pool at the parasite population level (Su et al., 1995; Chen et al., 1998; Hviid and Jensen, 2015; Deitsch and Dzikowski, 2017). Switches of var gene expression have been demonstrated in infections of human volunteers (Peters et al., 2002; Bachmann et al., 2019). Several mechanisms to explain the mutually exclusive var expression have been reported, including: (1) cis-acting DNA elements (promoters, silencers, activators) and anti-sense RNA transcripts acting to silence or activate the var genes; (2) epigenetic regulation through the presence or absence of specific epigenetic marks; (3) locations of the var genes in specific subnuclear compartments; (4) coordination between members of the var family to ensure activation of a new gene while simultaneously silencing of the previously active gene at the same time (Deitsch and Dzikowski, 2017). Similarly, there are also several potential mechanisms for the generation of new var genes and maintenance of the var gene repertoire. Although each parasite carries a set of approximately 60 var genes, a parasite strain usually has a unique set of var genes with diverse sequences that differ from those of other strains (Trimnell et al., 2006; Chen et al., 2011; Claessens et al., 2014; Day et al., 2017; Ruybal-Pesantez et al., 2017). In fact, it was reported that every child infected with P. falciparum had a unique var DBLα repertoire (Ruybal-Pesantez et al., 2017), and some specific types of homologous blocks or motifs have been associated with severe diseases (Avril et al., 2013; Rorick et al., 2013). Considering frequent mixed infection of multiple parasite strains in an individual patient in endemic regions (Conway et al., 1991; O'Brien et al., 2016; Zhu et al., 2019), the var gene repertoire in the parasite populations of an individual patient or in a village can be quite large.
What are the forces or mechanisms driving the high diversity of the var genes? Since the PfEMP1 proteins are expressed on the iRBC surface, they are expected to be under strong immune selection. In addition to switches in expression to different var genes, one obvious mechanism to evade immune recognition is substitutions of critical amino acids at antigenic epitopes. A second mechanism for generating var gene diversity is the production of hybrid copies of var genes through genetic recombination during meiosis in the mosquito host. If genetic recombination is the major driving force, then sequencing parasite strains in mixed infections from an individual patient would reveal shared haplotypes and breakpoints in some var genes among parasite strains from a mixed infection. Interestingly, transmission through a mosquito also resets the var gene switching process (Peters et al., 2002), suggesting that going through mosquitoes or a sexual cycle can change var gene expression profile and may also alter var gene composition in the genome at the same time. A third mechanism is generation of recombinant var genes through somatic recombination during asexual replication. Recombination events such as duplicative transposition of specific sequences were observed during parasite transformation in in vitro culture (Frank et al., 2008). By constructing large parasite clone trees and performing whole genome sequence analysis of var gene sequences in asexually replicating parasites, Claessens et al. (2014) found that the Var exon 1 recombined in a rate up to 0.2% of infected erythrocytes in vitro per life cycle, suggesting that var gene sequence polymorphism could be mainly generated during the asexual part of the life cycle (Claessens et al., 2014). Using targeted DNA double-strand breaks (DSBs) and long-read whole-genome sequencing, Zhang et al. (2019) showed that a cascade of recombination events could occur starting at a single DNA break within a subtelomeric region generating multiple new var genes (Zhang et al., 2019).
What is the mechanism mediating somatic recombination? One possible way of generating recombinant var genes during the asexual cycle is through microhomology-mediated end joining (MMEJ), involving alignment of microhomologous sequences flanking a DNA break and subsequent removal of the intermediate sequences (Grajcarek et al., 2019). Because of the absence of enzymes necessary for non-homologous end joining in Plasmodium, repair of DNA breaks generally required homologous templates (Lee et al., 2014). However, an end joining pathway for repair of DSBs without the presence of homologous sequence was observed in P. falciparum, with the repair junctions frequently containing short insertions in the surrounding sequences (Kirkman et al., 2014), suggesting that MMEJ may play a role in these events. Various enzymes for MMEJ are indeed present in P. falciparum parasites (Lee et al., 2014). Recently, MMEJ was shown to function in the P. yoelli genome and was able to generate a large number of deletion or insertion mutants of varying sizes in the central repeat region of P. yoelii CSP gene (Xu et al., 2019). Therefore, homologous regions among the different copies of the var genes may act as templates during DSB repair, leading to deletion, insertion, and translocation of repetitive DNA sequences and rapid evolution of new var genes.
The var genes have a GC content much higher than those of other protein coding genes in the P. falciparum genome (~34% vs. 20%), and how the genes with high GC content originated and are maintained in a genome with a high AT content remains largely unknown (Su et al., 1995; Gardner et al., 2002). One potential explanation is that the var gene came from the parasite hosts (humans and mosquitoes) through horizontal gene transfer; however, there is no evidence that the host genomes contain genes with structures and/or protein motifs observed in the PfEMP1. The var genes are not even present in the other human malaria parasites such as Plasmodium vivax, whose genome has high GC content in its genome (Carlton et al., 2008; Chan et al., 2015). Because var genes contain Duffy binding-like (DBL) domains that are also present in another gene family encoding erythrocyte invasion ligands (EBL proteins), it has been proposed that the var genes might have originated from the EBL genes (Smith et al., 2013). The var genes are also present in the Laverania parasites infecting African great apes (P. praefalciparum, P. blacklocki, P. adleri, P. billcollinsi, P. gaboni, and P. reichenowi) (Rask et al., 2010; Otto et al., 2014b, 2018; Laurens et al., 2017). Therefore, the var genes were already present when P. falciparum diverged from P. preafalciaprum to emerge as a human parasite 40,000–60,000 years ago (Otto et al., 2018). The question as how the var genes initiated in the ancestor remains unanswered. If var genes were derived from the genes containing DBL domains, why similar gene family expansion events did not occur in other malaria parasite species?
The pir Genes
Similar to the var genes, the pir multigene family is likely under host immune pressure and plays an important role in immune evasion (Cunningham et al., 2010). The pir gene family derives its name from the first initial of the parasite species name followed by “ir,” an abbreviation for interspersed repeat: yir for P. yoelii genes; bir for P. berghei genes; vir for P. vivax genes, and so on (Cunningham et al., 2010). Gene families of stevor, rif, and PfMC-2TM in P. falciparum were classified as pir gene families because of similar gene structures and sequence similarity in some introns (Janssen et al., 2002). The pir genes have a conserved three-exon structure, including a short first exon, a long second exon encoding conserved cysteine residues and a transmembrane domain, and a conserved third exon encoding cytoplasmic domain (Janssen et al., 2002). The pir genes are mostly distributed in the subtelomeric regions of chromosomes (except kir in P. knowlesi) with gene copies numbering from a few dozen to more than 800 (Cunningham et al., 2010). Distinct from the var genes, up to 40% of the cir gene repertoire are expressed throughout the intraerythrocytic cycle of development of the P. chaboudi parasite during infection (Lawton et al., 2012). The timing and level of transcription differs between cir genes, with ring, trophozoites, and schizonts expressing different sets of major cir genes (Cunningham et al., 2009). Additionally, individual CIR proteins were observed within differential localizations of iRBCs, including on the surface of iRBCs and merozoites, some of which could bind mouse erythrocytes and play a role in immune evasion (Yam et al., 2016; Yam and Preiser, 2017).
The mechanisms for generating pir gene diversity are expected to be similar to those of var genes, including amino acid substitutions at specific antigen epitopes in response to immune pressure and the generation of new recombinant genes through genetic recombination as well as MMEJ mediated somatic deletion, insertion and translocation. Genes with mosaic sequences or two genes with segments of 100% identity in the P. yoelii genome have been noted, suggesting crossovers between two copies of pir genes (Cunningham et al., 2010). However, rodent malaria parasites may have a lower degree of antigen diversity and undergo significantly less ectopic recombination than P. falciparum, partly due to the lack of components of the translesion (TLS) polymerases that are required for the recombination events that drive diversification of the multicopy gene families (Kirkman and Deitsch, 2020; Siao et al., 2020). Further discussions on functions of the pir genes can be found here (Cunningham et al., 2010; Su et al., 2019). More studies are necessary to better understand the mechanisms of gene expression regulation and the maintenance of genetic diversity of the pir genes.
The EBL Genes and Evolution of Virulence in Laboratories
One of the consequences of host immune pressure on parasites is the generation of parasites with altered virulence. Whereas it has been difficult to evaluate how host immune response or vaccination impacts parasite evolution and virulence in humans, changes in parasite growth and virulence can be measured in rodent malaria parasites in the laboratory. One example is the P. yoelii EBL protein (PyEBL) that is known to bind the Duffy antigen/receptor for chemokines (DARC), a transmembrane glycoprotein expressed on epithelial cells, endothelial cells, and erythrocytes (Swardson-Olver et al., 2002; Woolley et al., 2005; Culleton and Kaneko, 2010). Whereas several genes encoding proteins with DBL domains have been found in P. falciparum such as PfEBL-140, PfEBL-175, and PfEBL-181, there is only a single copy of the EBL gene in P. yoelii (Culleton and Kaneko, 2010). PyEBL contains a single DBL domain, unlike other EBLs which may have two DBL domains per molecule.
Immunity or host response to malaria infection is generally parasite strain specific (Fluck et al., 2004; Cheesman et al., 2006; Early et al., 2018). The principle of strain specific immunity has been applied to identify immune targets in rodent malaria parasites using a strategy called linkage group selection (LGS) (Culleton et al., 2005, 2011; Pattaradilokrat et al., 2009; Cheesman et al., 2010; Abkallo et al., 2017). In a LGS analysis, phenotype-specific selection pressure such as immunity is applied to uncloned progeny of a genetic cross between two malaria parasites that differ in phenotype (strain specific immunity in this case) (Culleton et al., 2005; Cheesman et al., 2010). PyEBL has also been identified by LGS experiments because it plays a role in parasite invasion of RBCs and is also the target of host immunity (Pattaradilokrat et al., 2009; Abkallo et al., 2017). Mutations in the PyEBL gene not only affect parasite growth and virulence, but also influence iRBC surface protein and host immune response (Otsuki et al., 2009; Peng et al., 2020). There are isogenic pairs of rodent malaria parasites that have different parasite growth rates in mice and produce very different diseases. P. yoelii yoelii 17X (or YM) infection of BALB/c or C57BL/6 mice causes host death within 7 days post infection (pi), whereas its isogenic strain P. y. yoelii 17XNL (or 17X) is cleared by the host around 21 days pi (Otsuki et al., 2009; Pattaradilokrat et al., 2009). Both parasites were derived from P. yoelii 17X during passages in separate laboratories (Pattaradilokrat et al., 2008). A replacement of Cys to Arg (at position 731 or C731R) at the second Cys position in region 6 of PyEBL was found to be the major determinant for the difference in virulence and parasite growth; however, this mutation cannot explain the entire parasite growth phenotype (Otsuki et al., 2009). A second example is the isogenic parasites of P. y. nigeriensis N67 (N67) and P. y. nigeriensis N67C (N67C) that also have very different disease phenotypes in C57BL/6 mice, including parasitemia, tissue pathology, and host mortality (Wu et al., 2014; Lacerda-Queiroz et al., 2017). Parasitemia in C57BL/6 infected with N67 (1 × 106 iRBCs) increases to ~40–50% day 5 pi, declines to below 10% on day 7 pi, and increases to ~60% at day 15 pi (Wu et al., 2014). The N67 infected mice die at ~day 20 pi. Mice infected with N67C have ~50% parasitemia day 6 pi and begin to die day 7 pi with declining parasitemia (Lacerda-Queiroz et al., 2017). Similarly, a substitution of Cys to Tyr (C741Y) in the protein trafficking domain of PyEBL between N67 and N67C parasites was found to influence disease phenotypes and host mortality (Peng et al., 2020). The substitution alters processing and trafficking of the PyEBL protein, reduces PyEBL binding to Band 3 of the RBC, increases phosphatidylserine (PS) surface exposure on iRBCs, and enhances iRBC osmotic fragility. Interestingly, the two studies have different conclusions on the mechanism of how the PyEBL substitutions affecting disease phenotypes, although the C731R and C741Y substitutions both appear to change PyEBL protein trafficking, from microneme localization to dense granules in the cytoplasm. For C731R substitution between 17X and 17XL (YM), parasite virulence was mostly attributed to altered erythrocyte invasion (Otsuki et al., 2009). In contrast, in the C741Y substitution in PyEBL between N67 and N67C was found to modulate host recognition of iRBCs and immune responses through interaction with RBC membrane molecules such as Band 3 and PS (Peng et al., 2020). A third example is the isogenic parasites P. berghei ANKA and P. berghei NK65: P. berghei ANKA infection is lethal with cerebral malaria symptoms and has been used as an experimental cerebral malaria (ECM) model, whereas P. berghei NK65 infection generally does not have neurological symptoms (Lacerda-Queiroz et al., 2011; Nacer et al., 2014). Even cloned lines of P. berghei ANKA can differ in their ability to induce ECM in the same host (Amani et al., 1998). However, the molecular mechanisms responsible for the differences in disease phenotype and virulence among these parasites are not clear. Additionally, the substitutions in PyEBL are likely not the only genetic changes contributing to the total disease phenotypes. Unknown differences including SNPs and variation in genes such as the pir multigene gene family may also contribute the differences in disease phenotypes and virulence (Cunningham et al., 2010). These parasite pairs are examples of the evolution of parasite genomes and changes in virulence following laboratory passages. There are additional parasite species and strains, including P. chabaudi and P. vinckei, that showed changes in parasite genome and virulence (Perkins et al., 2007). It can be expected that these types of phenotypic changes in parasite growth, disease virulence, and drug response are occurring daily in the huge parasite populations in endemic regions of the world. Additionally, mixed infections of multiple P. falciparum strains may lead to higher proportion of severe anemia, pulmonary complications, and multiple organ failure than those with single-strain P. falciparum infection (Kotepui et al., 2020). The dynamics of parasite genetic variation and gene expression are major roadblocks for developing effective vaccines and therapies.
Parasite Selection on the Host Genome
Whereas host immunity is one of the major factors driving parasite evolution, malaria parasites and many other pathogens also have tremendous impacts on the adaption and evolution of their hosts (Daub et al., 2013). The well-known Haldane malaria hypothesis suggests an evolutionary advantage of the thalassemic condition in protection against malaria (Akide-Ndunge et al., 2003; Canali, 2008). Many proteins expressed on the RBC surface such as sickle cell HbS variants, DARC, and glycophorin A (GYPA), and proteins that play a role in host immune response such as TNF-α and IFN-γ have been associated with protection from severe complications of malaria (Carter and Mendis, 2002; Kwiatkowski, 2005; Verra et al., 2009; Hedrick, 2012; Kariuki and Williams, 2020). The protective genes were mostly identified through association studies including genome-wide association studies (GWAS). For example, a recent GWAS study performed in a Tanzanian population revealed association of protective roles of sickle cell HbS variant, interleukin receptors (IL-23R and IL-12RBR2), and the kelch-like protein KLHL3 (Ravenhall et al., 2018). In another study, Ebel et al. analyzed mammalian adaptation in 490 Plasmodium- or Piroplasm-interacting proteins (PPIPs) and showed that blood parasites imposed specific selection on PPIPs, leading to a higher number of adaptive substitutions in the PPIPs than expected throughout mammalian evolution (Ebel et al., 2017). Interestingly, many immune PPIPs were also enriched for viral and bacterial interactions, with 48% of the PPIPs interacting with viruses and/or bacteria. Moreover, there was a 2.5X enrichment of adaptation when only PPIPs interacting with Plasmodium or Piroplasms were considered. Malaria parasites interact with a large number of host genes (Wu et al., 2015, 2020), which may enhance adaptive substitutions in these genes. There are many excellent reviews on the effects of parasite infections, particularly malaria parasite infections, on host genome evolution (Tishkoff and Williams, 2002; Kwiatkowski, 2005; Mangano and Modiano, 2014; Goheen et al., 2017).
Impact of Vaccines, Drugs, and Vectors on Parasite Populations and Disease
Types of Malaria Vaccines
Developing an effective vaccine has been one of the major goals of malaria research for more than 40 years, and we are still far away from achieving the goal. Malaria vaccines can be divided into several categories based on the parasite developmental stages such as pre-erythrocytic vaccines, blood stage vaccines including those against pregnancy-associated malaria, transmission-blocking vaccines, multistage or multiantigen vaccines, and whole organism based vaccines (Frimpong et al., 2018; Laurens, 2018; Beeson et al., 2019). Development of a cross-species vaccine that protects against multiple Plasmodium species or strains is an important goal. P. falciparum 3D7 strain has been the model strain for studies of genome structure and analyses of transcriptome and proteome. Subsequently, the majority of vaccines are based on sequences from the 3D7 parasite; however, the 3D7 strain is often not the most prevalent strain in malaria-endemic areas. Currently, the leading vaccine candidates are mostly based on polymorphic antigens with the vaccine target representing only a single variant strain. However, immunity to malaria is generally strain-specific, and the development of an effective vaccine is often hampered by strain-specific and/or partial immunity generated by vaccination.
Partially Protective Vaccines and Effects on Host Populations
The most promising and extensively tested malaria vaccine candidate against P. falciparum parasites is RTS,S/AS01 that consists of 18 copies of the central repeat and the C-terminal domain of PfCSP fused to hepatitis B virus surface antigen (HBsAg) (Draper et al., 2018). In a large-scale phase III trial involving 15,460 children, including 6,537 infants aged 6–12 weeks and 8,923 children aged 5–17 months, vaccine efficacy (VE) against clinical malaria was reported to be ~50–56%, and VE against severe malaria in the combined age children was ~35% in the per-protocol population at 11 month follow-up (Rts et al., 2011). However, in a 7 year follow-up analysis, the VEs of RTS,S/AS01 were estimated to be 35.9% in the first year after vaccination, 2.5% in the fourth year, and a negative VE in the fifth year (Olotu et al., 2016). More generalized convulsive seizures were observed in the RTS,S vaccinated group than the control group (comparator rabies vaccine group) in the older children (5–17 months) within 7 days after RTS,S/AS01 vaccination (Rts et al., 2011). Additionally, significantly more cases of meningitis were reported in the older children (Rts et al., 2012; Rts, 2015). Moreover, RTS,S vaccine was later found to be associated with significantly higher all-cause mortality in girls of both age groups (Klein et al., 2016). The higher rates of meningitis and mortality in girls have prompted the World Health Organization (WHO) to issue a position paper cautioning on malaria vaccines, with emphasis on the RTS,S/AS01 (World Health Organization, 2018) and discussions on the need of further studies on vaccine safety and negative non-specific effects (Aaby et al., 2015; Muller et al., 2015). Higher all-cause mortality was also observed in girls after vaccination with the high-titer measle vaccine (HTMV) (Aaby et al., 1993), which was later found to be likely associated with diphtheria-tetanus-pertussis (DTP) vaccine or interaction between HTMV and DTP (Aaby et al., 2003, 2007). Generally, women have a stronger cellular and humoral responses to infections and are more likely to have autoimmune diseases due to X-linked genes, hormones, and societal context (Fish, 2008; Billi et al., 2019). The observation of higher rate of severe symptoms in RTS,S vaccinated girls could be due to enhanced immune reactions from vaccine priming of the immune system.
A partial protective vaccine with enhanced adverse effects brings out some important issues in malaria vaccine development that has largely been ignored so far: Can a partially protective malaria vaccine prime for more severe disease symptoms and select for parasite population “resistant” to specific vaccines? The T-cell prime hypothesis proposed to explain clinical malaria and antibody dependent enhancement (ADE) observed in many viral infections may provide some answers to the questions (Porterfield, 1986; Riley, 1999). More than 20 years ago, Dr. Eleanor Riley proposed that T-cell priming is required for amplification of the inflammatory response to malaria, which can explain the patterns of clinical malaria in both endemic and non-endemic populations (Riley, 1999). According to the hypothesis, primary infections would induce low levels of IFN-γ and TNF-α; at the same time, antigen-specific T cells are primed by the infection. On reinfection, large amounts of cytokines such as IFN-γ and TNF-α are produced by the primed T cells, leading to an increased risk of cerebral malaria or other severe symptoms. She also cautioned that vaccinations with BCG or tetanus toxoid might prime T cells to respond to malaria antigens and suggested that having an appropriate balance between protective and pathogenic levels of Th1-derived cytokines will be crucial for a successful outcome of a vaccine (Riley, 1999). A partial protective vaccine may also prime the immune system in some individuals that can mount a strong inflammatory response when infected, leading to enhanced severe symptoms. Another related possibility is ADE that was observed in various viral infections, including Dengue virus (Waggoner et al., 2020), influenza (Winarski et al., 2019), Zika virus (Asad et al., 2019), coronavirus (Negro, 2020; Wan et al., 2020), and Acinetobacter baumannii bacteria (Wang-Lin et al., 2019). ADE is a well-known mechanism that viruses or bacteria may infect susceptible cells by interacting with antibodies or complement components (Fc or complement receptors), leading to enhanced viral entry into host cells and/or higher replication rate. Cross-reactive IgGs can opsonize viral particles by binding the virus through Fc receptors on host cell surface and may also mediate immune suppression to further increase viral load (Halstead et al., 2010). Monoclonal antibodies have been shown to destabilize the influenza HA stem domain, leading to faster kinetics of influenza virus fusion and enhanced respiratory disease (Winarski et al., 2019). Antibody enhancement of malaria transmission and complement activation in conjunction with antibodies in enhancing RBC invasion have been reported too (Biryukov et al., 2016; Stone et al., 2019). The RTS,S vaccine might prime T cells in some hosts and/or produce antibodies that can amplify immune response and enhance clinical symptoms. Addition of an arm to evaluate the potential of a vaccine to enhance adverse effects should be included before or during a vaccine trial.
Effects of Vaccination on Parasite Population and Gene Expression
Malaria parasite populations in endemic regions consist of a large number of strains that express different alleles of vaccine target proteins. Due to technical limitation, malaria vaccines are usually designed based on sequences from one to a limited number of alleles of a vaccine target. It is has been shown that anti-malaria immunity is mostly strain-specific (Cheesman et al., 2006; Graves et al., 2016; Laurens et al., 2017; Early et al., 2018), and a vaccine based on a limited number of alleles of a target may select for parasite populations with alternative alleles. Indeed, vaccination with pneumococcal conjugate vaccines (PCVs) changed the pneumococcal populations in children, resulting in statistically significant shift from vaccine-type population to non-vaccine-type populations (Quirk et al., 2018). Similarly, using the rodent malaria model P. chabaudi and recombinant AMA-1 antigen, it was shown that mono-allelic immunization increased the frequency of heterologous clones in mixed clone infections (Barclay et al., 2008). Moreover, analysis of parasite genotypes collected from the RTS,S phase III trial showed that the RTS,S vaccine had greater activity against malaria parasites with the matched PfCSP allele than against those with mismatched alleles (Neafsey et al., 2015). Therefore, vaccination with an allele of the target antigen will likely change target allele proportion in parasite populations, and large-scale vaccination may lead to vaccine mediated depletion of specific alleles targeted by the vaccine (Figure 3B). The impacts of this type of selection on parasite populations, including the possibility of selecting more virulent parasite strains, remain largely unknown and required further investigations. Interestingly, a vaccine may also alter parasite expression of variant antigen genes. Vaccination of 10 African volunteers with the PfSPZ vaccine showed parasites from individuals with intermediate antibody levels expressed only few var gene variants, whereas those with low antibody levels expressed a broad spectrum of multiple, predominantly subtelomeric var genes coding PfEMP1 binding to endothelial protein C receptor (EPCR) that is associated with severe childhood malaria (Bachmann et al., 2019). Therefore, vaccination certainly can impact parasite population dynamics, genetic diversity, and possibly virulence.
Drug Selection on Parasite Populations and Evolution
In addition to host immune response, antimalarial drugs have played a significant role in shaping parasite populations. Application of antimalarial drugs puts pressure on parasite populations and likely will select parasites with mutations that can confer resistance to the drugs, leading drug selective sweeps (Figure 3C). Some examples of drug selective sweep include mutations in P. falciparum chloroquine resistance transporter (PfCRT) and in dihydrofolate reductase-thymidylate synthase and dihydropteroate synthase (PfDHFR-TS and PfDHPS) conferring resistances to chloroquine and pyrimethamine/sulfadoxine (PS), respectively (Peterson et al., 1990; Fidock et al., 2000). Under drug pressure, a small number of parasites with resistant mutations will survive, whereas parasites without the mutations are killed by the drugs. With large-scale drug use, parasites with resistant mutations will spread, replacing parasites sensitive to the drugs. Indeed, worldwide chloroquine and PS selective sweeps have been reported (Wootton et al., 2002; Roper et al., 2004; McCollum et al., 2008), which can greatly reduce population diversity in many endemic areas generating population bottlenecks. Drug selective sweeps are important factors in shaping the current worldwide parasite populations and may have contributed to Malaria's Eve hypothesis debate (Rich et al., 1998; Su et al., 2003). The sweeps of parasites derived from one or a few mutants conferring resistance to a drug will result in a relatively homogeneous parasite population (Figure 3C), as was seen in the Brazilian Amazon after wide spread use of chloroquine (Wootton et al., 2002). A relatively homogeneous parasite population may be advantageous for vaccine development.
The Role of Mosquito Vectors in Altering Parasite Populations
Mosquito vectors can shape parasite populations in many ways (Figure 3D). First, genetic changes such as chromosome reassortment and crossover will occur inside a mosquito and generate new parasite strains if a patient is infected with two or more variant strains, which will increase parasite population diversity. A high multiplicity of infection is frequently observed in patients in regions with high transmission intensity (Kobbe et al., 2006; Felger et al., 2012). Second, parasite transmission through mosquitoes can also re-program gene expression profiles. For example, changes in composition and frequency of var gene transcripts were observed between cultured P. falciparum parasites used to infect mosquitoes and the parasites recovered from infected volunteers after mosquito bites, suggesting re-programing var gene expression profile (Peters et al., 2002). Similarly, attenuated parasite virulence was found to be associated with modified expression of the pir multi-gene family. Transmission of P. c. chabaudi through the mosquito changes gene expression of the pir multi-gene family in the erythrocytic cycle, leading to changes in parasite virulence and immune response in the mammalian hosts (Spence et al., 2013). Third, vector competency can also influence parasite transmission frequency. A study of P. vivax populations in southern Mexico showed that the distributions of parasite populations largely mirrored those of mosquito species, which was confirmed by feeding experiments in the laboratory (Joy et al., 2008). The results suggest a reciprocal selection between malaria parasites and mosquito vectors and local adaptation of the parasites to mosquito species in different environments (coastal or hill sites). P. falciparum isolates from Africa, Asia, or the Americas were shown to have low compatibility to malaria vectors from different continents, which is mainly mediated by a parasite protein called Pfs47 due to its interaction with mosquito immune system (Molina-Cruz et al., 2015). Mosquito phagocytes play an important role in mosquito vector competence and immunity against malaria infection (Kwon and Smith, 2019). Additionally, Serratia marcescens bacteria isolated from either laboratory-reared mosquitoes or wild populations in Burkina Faso shows great phenotypic variation that is directly correlated with its ability to inhibit Plasmodium development within the mosquito midgut (Bando et al., 2013).
Approaches for Studying Host-Parasite Interactions
Identification of Host-Parasite Interacting Genes Using Genetic Crosses
In addition to antibody response to antigens expressed on the surfaces of iRBCs, malaria infections also stimulate differential activation of immune cells and the production of cytokines and chemokines as part of the innate immune response. What are the parasite ligands that initially trigger host innate responses? What are the host receptors that recognize the parasite molecules? How are the host responses to malaria infections regulated? What are the impacts of the host-parasite interaction on the evolution of both the host and parasite genomes? As discussed above, host immunity drives diversity of parasite antigen genes, and parasite infections select for protective alleles in host response genes. Large numbers of parasite antigens and host protective genes have been identified (Kwiatkowski, 2005; Salinas et al., 2019; Vijayan and Chitnis, 2019; de Jong et al., 2020). However, the numbers of host and parasite molecules involved in host-parasite interactions are likely much larger than the numbers of genes so far investigated. For example, Wu et al. (2015) performed an immunogenetic screen using trans-species expression quantitative trait locus (ts-eQTL) analysis and significantly (LOD score ≥ 3.0) linked 1,054 host genes to a large number of parasite genetic markers/loci (Wu et al., 2015). In this study, microsatellite genotypes of the P. yoelii malaria parasite and phenotypes of host (mouse) gene expression levels day 4 pi were analyzed for significant linkages between parasite genetic markers and host gene expression levels. Additionally, host genes that were linked (LOD score ≥ 3.0 or 2.0) to parasite genetic loci were clustered based on genome-wide patterns of LOD scores (GPLSs), showing that genes functioning in related pathways generally have similar GPLSs. This study demonstrates that a large number of host genes respond to P. yoelii parasites as early as day 4 pi. Similarly, many parasite genetic loci were predicted to interact with host genes belonging to specific pathways. For example, a genetic locus on parasite chromosome 8 was significantly linked to a large number of host genes play a role in hematopoiesis, suggesting that a gene or genes in the locus may influence host hematopoiesis. The study represents the first QTL analysis involving two organisms, microsatellite genotypes from malaria parasites and gene expression as phenotypes from the host, with the advantage of identifying interacting genes from both the parasite and the host. This ts-QTL analysis is a novel approach for identifying interacting genes from both parasite and the host at the same time, taking advantage of genetic recombination and chromosome reassortment during the sexual stage of a pathogen in generating progeny with unique genomes inherited from the parents. In a follow-up study using mRNA isolated from the spleens of infected mice 24 h pi (earlier than day 4), many host receptors including Toll-like receptors (TLR), retinoic acid-inducible gene-I-like receptors (RLRs), G-protein coupled receptors (GPCRs), and olfactory receptors (ORs) were significantly linked to parasite genetic loci, providing important candidate genes for studying host molecules that recognize parasite ligands or pathogen-associated molecular patterns (PAMPs) (Wu et al., 2020). Ts-QTL holds the potential to identify more interacting genes using additional genetic crosses or host RNA samples collected at different time points throughout the parasite life cycle.
Transcriptome Analysis of Patient Samples
Another approach to identify host gene responses to malaria infections is to use samples collected from malaria patients, although a study using patient blood samples can be influenced by many factors such as variations in infection time, host genetic background and immune status, mixed parasite infections, nutrition status, antimalarial drug usages, etc. Mixed strain infections can change parasite growth dynamics, relapse time (for P. vivax), disease severity, and host mortality rate (Lover and Coker, 2015; Kotepui et al., 2020). Similarly, the genetic composition of each patient is unique, and variations in host genetic background will also affect host responses. Therefore, dissecting host-parasite genetic interactions using patient samples has proven to be much more complicated than using infections in inbred mice with cloned progeny from a genetic cross.
As discussed above, GWAS has been successful in identifying some host genes that contribute to protection against malaria. As an alternative, analysis of host gene expression in response to malaria infection can be useful too, particularly when transcripts from both the host and parasites are analyzed at the same time using RNAseq (Greenwood et al., 2016). Using blood samples of P. falciparum-infected West African children, Idaghdour et al. (2012) performed a joint genome-wide analysis of gene expression and genetic polymorphisms, and showed that malaria infection and host genotype both influenced immune gene expression throughout the genome (Idaghdour et al., 2012). The study implicated many genes in the complement system, antigen processing and presentation, and T-cell activation in response to malaria infections (Idaghdour et al., 2012). For example, the authors noted that the expression of five genetically regulated human leukocyte antigen (HLA) class II loci was negatively correlated with parasite load and molecules such as IL18R1, TLR4, TLR5, IFNGR1, and IFNGR2, suggesting impairment of antigen processing and immune signaling. Impairment of antigen processing and T cell activation, including suppression of various MHCII genes, were also observed in C57BL/6 mice infected with different P. yoelii strains and progeny of a P. yoelii genetic cross (Xia et al., 2018). In a dual transcriptome analyses of the host and parasite genes in samples from 46 malaria-infected Gambian children, Lee et al. (2018) performed multivariate analyses and showed association of disease severity with increased expression of granulopoiesis and IFN-γ related genes as well as inadequate suppression of IFN-I signaling (Lee et al., 2018). In another analysis by the same group, lower expression of cxcl10, HLA genes, and IFN-I associated genes, and higher expression of cathepsin G and matrix metallopeptidase 9 genes were associated with parasite growth inhibition (Georgiadou et al., 2019). Moreover, a blood transcriptional analysis in Senegalese patients with cerebral malaria (CM), severe non-cerebral malaria (NCM), or mild malaria (MM) identified host gene clusters that discriminated between CM and MM patients, showing that genes involved in BCR-, TCR-, TLR-, cytokine-, FcεRI-, and FCGR- signaling pathways played a role in cerebral malaria (Thiam et al., 2019). Associations between disease severity and PfEMP1 transcripts and specific domains were also observed in RNAseq data from 44 parasite isolates that caused SM and MM malaria in Papuan patients (Tonkin-Hill et al., 2018). These studies illustrate that there are signatures of host transcriptional responses to malaria infections and disease severity. Indeed, blood transcriptome analyses of pathogen-specific host immune responses have been proposed as a tool for improving the diagnosis and classification of patients with infectious diseases (Mejias et al., 2014). However, among the studies of host transcriptional response to malaria infections in animal models and in humans, the results are not always consistent, likely due to differences in study design, study objective, transmission intensity, sample collection time, as well as lack of replication. Caution is warranted when interpreting the results from field-based analyses.
There are other approaches for studying host-parasite interaction such as analysis of genome diversity and genetic association (Early et al., 2018; Damena et al., 2019; Malaria Genomic Epidemiology Network, 2019), computational prediction of host-parasite interactomes (Wuchty, 2011; Cuesta-Astroz et al., 2019), metabolite profiling to investigate transcriptomic changes during malaria infection (Tang et al., 2018), and systems biology/integrative approaches to study host immune responses (Mahanta et al., 2018). These are just few emerging and promising examples using -omics data to study host-parasite interactions.
Some Key Questions and Challenges
There are several area that require special attention: (1) The new Laverania species are identified through nucleic acid sequence analyses. Morphological and biological characterization of the parasites will be critical for establishing these parasites as solid independent species. (2) The molecular mechanism of host specificity of the Laverania species is interesting and requires further study. (3) The origin of the P. falciparum var genes remains unknown. How did this gene family with a high GC-content originated and how is it maintained in a AT-rich genome? (4) Similar to the var gene family, the origins and functions of the pir gene families require additional studies. The pir genes likely play important roles in virulence and immune evasion. (5) More studies are necessary to elucidate the mechanisms that generate and maintain diversity of the var and pir genes in a parasite cell. (6) The impacts of antimalarial drugs and vaccines on the diversity and evolution of malaria parasite populations are important issues related to disease control and management. (7) The long-term impacts of antimalarial drugs and vaccines on host health and immune system also require more attentions. ADE is a real possibility for malaria vaccines. (8) Insect vectors play important roles in the evolution, diversity and transmission of malaria parasites. This is an area that has not received sufficient attention. (9) The parasite molecules that trigger host innate immune responses remain largely uncharacterized. DNA and RNA have been implicated in stimulating type I interferon responses, and strain-specific responses have been observed. However, the specific sequences that can bind to host receptors and the mechanisms of the strain-specific innate immune responses require additional studies.
Conclusion
Host-parasite interactions are the main forces driving the evolution of both malaria parasites and their hosts. Genome diversity and gene function have co-evolved shaping the current genomes of malaria parasites and their hosts. Application of antimalarial drugs, deployment of vaccines, modulation of host immunity, and the development in mosquito vector control measures all have the potential to alter parasite populations including genome diversity and virulence. Proper monitoring of parasite populations in the context of drug, vaccine, and vector control programs is necessary to minimize unintended negative impacts on parasite population dynamics and disease severity.
Author Contributions
X-zS and DAJ conceive ideas and writing. CZ wrote and created Figure 1. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Division of Intramural Research, National Institute of Allergy and Infectious Diseases (NIAID), National Institutes of Health (NIH), USA.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank Bradley Otterson, NIH Library Writing Center, for manuscript editing assistance.
References
Aaby, P., Garly, M. L., Nielsen, J., Ravn, H., Martins, C., Bale, C., et al. (2007). Increased female-male mortality ratio associated with inactivated polio and diphtheria-tetanus-pertussis vaccines: observations from vaccination trials in Guinea-Bissau. Pediatr. Infect. Dis. J. 26, 247–252. doi: 10.1097/01.inf.0000256735.05098.01
Aaby, P., Jensen, H., Samb, B., Cisse, B., Sodemann, M., Jakobsen, M., et al. (2003). Differences in female-male mortality after high-titre measles vaccine and association with subsequent vaccination with diphtheria-tetanus-pertussis and inactivated poliovirus: reanalysis of West African studies. Lancet 361, 2183–2188. doi: 10.1016/S0140-6736(03)13771-3
Aaby, P., Knudsen, K., Whittle, H., Lisse, I. M., Thaarup, J., Poulsen, A., et al. (1993). Long-term survival after Edmonston-Zagreb measles vaccination in Guinea-Bissau: increased female mortality rate. J. Pediatr. 122, 904–908. doi: 10.1016/S0022-3476(09)90015-4
Aaby, P., Rodrigues, A., Kofoed, P. E., and Benn, C. S. (2015). RTS,S/AS01 malaria vaccine and child mortality. Lancet 386, 1735–1736. doi: 10.1016/S0140-6736(15)00693-5
Abkallo, H. M., Martinelli, A., Inoue, M., Ramaprasad, A., Xangsayarath, P., Gitaka, J., et al. (2017). Rapid identification of genes controlling virulence and immunity in malaria parasites. PLoS Pathog. 13:e1006447. doi: 10.1371/journal.ppat.1006447
Adler, S. (1923). Malaria in chinpanzees in Sierra Leone. Ann. Trop. Med. Parasitol. 17, 13–19. doi: 10.1080/00034983.1923.11684345
Akide-Ndunge, O. B., Ayi, K., and Arese, P. (2003). The Haldane malaria hypothesis: facts, artifacts, and a prophecy. Redox Rep. 8, 311–316. doi: 10.1179/135100003225002952
Amambua-Ngwa, A., Tetteh, K. K., Manske, M., Gomez-Escobar, N., Stewart, L. B., Deerhake, M. E., et al. (2012). Population genomic scan for candidate signatures of balancing selection to guide antigen characterization in malaria parasites. PLoS Genet. 8:e1002992. doi: 10.1371/journal.pgen.1002992
Amani, V., Boubou, M. I., Pied, S., Marussig, M., Walliker, D., Mazier, D., et al. (1998). Cloned lines of Plasmodium berghei ANKA differ in their abilities to induce experimental cerebral malaria. Infect. Immun. 66, 4093–4099. doi: 10.1128/IAI.66.9.4093-4099.1998
Ansari, H. R., Templeton, T. J., Subudhi, A. K., Ramaprasad, A., Tang, J., Lu, F., et al. (2016). Genome-scale comparison of expanded gene families in Plasmodium ovale wallikeri and Plasmodium ovale curtisi with Plasmodium malariae and with other Plasmodium species. Int. J. Parasitol. 46, 685–696. doi: 10.1016/j.ijpara.2016.05.009
Asad, S., Feitosa-Suntheimer, F., Gold, A., Londono-Renteria, B., and Colpitts, T. M. (2019). Quantification of antibody-dependent enhancement of the Zika virus in primary human cells. J. Vis. Exp. 143, 1–6. doi: 10.3791/58691
Ashley, E. A., Pyae Phyo, A., and Woodrow, C. J. (2018). Malaria. Lancet 391, 1608–1621. doi: 10.1016/S0140-6736(18)30324-6
Auburn, S., Bohme, U., Steinbiss, S., Trimarsanto, H., Hostetler, J., Sanders, M., et al. (2016). A new Plasmodium vivax reference sequence with improved assembly of the subtelomeres reveals an abundance of pir genes. Wellcome Open Res. 1:4. doi: 10.12688/wellcomeopenres.9876.1
Avril, M., Brazier, A. J., Melcher, M., Sampath, S., and Smith, J. D. (2013). DC8 and DC13 var genes associated with severe malaria bind avidly to diverse endothelial cells. PLoS Pathog. 9:e1003430. doi: 10.1371/journal.ppat.1003430
Bachmann, A., Bruske, E., Krumkamp, R., Turner, L., Wichers, J. S., Petter, M., et al. (2019). Controlled human malaria infection with Plasmodium falciparum demonstrates impact of naturally acquired immunity on virulence gene expression. PLoS Pathog. 15:e1007906. doi: 10.1371/journal.ppat.1007906
Baer, K., Roosevelt, M., Clarkson, A. B. Jr., van Rooijen, N., Schnieder, T., and Frevert, U. (2007). Kupffer cells are obligatory for Plasmodium yoelii sporozoite infection of the liver. Cell. Microbiol. 9, 397–412. doi: 10.1111/j.1462-5822.2006.00798.x
Bando, H., Okado, K., Guelbeogo, W. M., Badolo, A., Aonuma, H., Nelson, B., et al. (2013). Intra-specific diversity of Serratia marcescens in Anopheles mosquito midgut defines Plasmodium transmission capacity. Sci. Rep. 3:1641. doi: 10.1038/srep01641
Barclay, V. C., Chan, B. H., Anders, R. F., and Read, A. F. (2008). Mixed allele malaria vaccines: host protection and within-host selection. Vaccine 26, 6099–6107. doi: 10.1016/j.vaccine.2008.09.004
Baruch, D. I., Pasloske, B. L., Singh, H. B., Bi, X., Ma, X. C., Feldman, M., et al. (1995). Cloning the P. falciparum gene encoding PfEMP1, a malarial variant antigen and adherence receptor on the surface of parasitized human erythrocytes. Cell 82, 77–87. doi: 10.1016/0092-8674(95)90054-3
Beeson, J. G., Kurtovic, L., Dobano, C., Opi, D. H., Chan, J. A., Feng, G., et al. (2019). Challenges and strategies for developing efficacious and long-lasting malaria vaccines. Sci. Transl. Med. 11:eaau1458. doi: 10.1126/scitranslmed.aau1458
Billi, A. C., Kahlenberg, J. M., and Gudjonsson, J. E. (2019). Sex bias in autoimmunity. Curr. Opin. Rheumatol. 31, 53–61. doi: 10.1097/BOR.0000000000000564
Biryukov, S., Angov, E., Landmesser, M. E., Spring, M. D., Ockenhouse, C. F., and Stoute, J. A. (2016). Complement and antibody-mediated enhancement of red blood cell invasion and growth of Malaria parasites. EBioMedicine 9, 207–216. doi: 10.1016/j.ebiom.2016.05.015
Blacklock, E., and Adler, S. (1922). A parasite resembling Plasmodium falciparum in chimpanzee. Ann. Trop. Med. Prasitol. 16, 99–107. doi: 10.1080/00034983.1922.11684303
Bohme, U., Otto, T. D., Cotton, J. A., Steinbiss, S., Sanders, M., Oyola, S. O., et al. (2018). Complete avian malaria parasite genomes reveal features associated with lineage-specific evolution in birds and mammals. Genome Res. 28, 547–560. doi: 10.1101/gr.218123.116
Boundenga, L., Ollomo, B., Rougeron, V., Mouele, L. Y., Mve-Ondo, B., Delicat-Loembet, L. M., et al. (2015). Diversity of malaria parasites in great apes in Gabon. Malar. J. 14:111. doi: 10.1186/s12936-015-0622-6
Bourgard, C., Albrecht, L., Kayano, A., Sunnerhagen, P., and Costa, F. T. M. (2018). Plasmodium vivax biology: insights provided by genomics, transcriptomics and proteomics. Front. Cell. Infect. Microbiol. 8:34. doi: 10.3389/fcimb.2018.00034
Brancucci, N. M. B., Gerdt, J. P., Wang, C., de Niz, M., Philip, N., Adapa, S. R., et al. (2017). Lysophosphatidylcholine regulates sexual stage differentiation in the human malaria parasite Plasmodium falciparum. Cell 171, 1532–1544.e1515. doi: 10.1016/j.cell.2017.10.020
Canali, S. (2008). Researches on thalassemia and malaria in Italy and the origins of the “Haldane hypothesis”. Med. Secoli. 20, 827–846.
Carlton, J., Silva, J., and Hall, N. (2005). The genome of model malaria parasites, and comparative genomics. Curr. Issues Mol. Biol. 7, 23–37.
Carlton, J. M., Adams, J. H., Silva, J. C., Bidwell, S. L., Lorenzi, H., Caler, E., et al. (2008). Comparative genomics of the neglected human malaria parasite Plasmodium vivax. Nature 455, 757–763. doi: 10.1038/nature07327
Carter, R., and Mendis, K. N. (2002). Evolutionary and historical aspects of the burden of malaria. Clin. Microbiol. Rev. 15, 564–594. doi: 10.1128/CMR.15.4.564-594.2002
Chan, E. R., Barnwell, J. W., Zimmerman, P. A., and Serre, D. (2015). Comparative analysis of field-isolate and monkey-adapted Plasmodium vivax genomes. PLoS Negl. Trop. Dis. 9:e0003566. doi: 10.1371/journal.pntd.0003566
Cheesman, S., O'Mahony, E., Pattaradilokrat, S., Degnan, K., Knott, S., and Carter, R. (2010). A single parasite gene determines strain-specific protective immunity against malaria: the role of the merozoite surface protein I. Int. J. Parasitol. 40, 951–961. doi: 10.1016/j.ijpara.2010.02.003
Cheesman, S., Raza, A., and Carter, R. (2006). Mixed strain infections and strain-specific protective immunity in the rodent malaria parasite Plasmodium chabaudi chabaudi in mice. Infect. Immun. 74, 2996–3001. doi: 10.1128/IAI.74.5.2996-3001.2006
Chen, D. S., Barry, A. E., Leliwa-Sytek, A., Smith, T. A., Peterson, I., Brown, S. M., et al. (2011). A molecular epidemiological study of var gene diversity to characterize the reservoir of Plasmodium falciparum in humans in Africa. PLoS ONE 6:e16629. doi: 10.1371/journal.pone.0016629
Chen, Q., Fernandez, V., Sundstrom, A., Schlichtherle, M., Datta, S., Hagblom, P., et al. (1998). Developmental selection of var gene expression in Plasmodium falciparum. Nature 394, 392–395. doi: 10.1038/28660
Claessens, A., Hamilton, W. L., Kekre, M., Otto, T. D., Faizullabhoy, A., Rayner, J. C., et al. (2014). Generation of antigenic diversity in Plasmodium falciparum by structured rearrangement of Var genes during mitosis. PLoS Genet. 10:e1004812. doi: 10.1371/journal.pgen.1004812
Conway, D. J. (2015). Paths to a malaria vaccine illuminated by parasite genomics. Trends Genet. 31, 97–107. doi: 10.1016/j.tig.2014.12.005
Conway, D. J., Greenwood, B. M., and McBride, J. S. (1991). The epidemiology of multiple-clone Plasmodium falciparum infections in Gambian patients. Parasitology 103(Pt. 1), 1–6. doi: 10.1017/S0031182000059217
Cuesta-Astroz, Y., Santos, A., Oliveira, G., and Jensen, L. J. (2019). Analysis of predicted host-parasite interactomes reveals commonalities and specificities related to parasitic lifestyle and tissues tropism. Front. Immunol. 10:212. doi: 10.3389/fimmu.2019.00212
Culleton, R., and Kaneko, O. (2010). Erythrocyte binding ligands in malaria parasites: intracellular trafficking and parasite virulence. Acta Trop. 114, 131–137. doi: 10.1016/j.actatropica.2009.10.025
Culleton, R., Martinelli, A., Hunt, P., and Carter, R. (2005). Linkage group selection: rapid gene discovery in malaria parasites. Genome Res. 15, 92–97. doi: 10.1101/gr.2866205
Culleton, R. L., Inoue, M., Reece, S. E., Cheesman, S., and Carter, R. (2011). Strain-specific immunity induced by immunization with pre-erythrocytic stages of Plasmodium chabaudi. Parasite Immunol. 33, 73–78. doi: 10.1111/j.1365-3024.2010.01251.x
Cunningham, D., Fonager, J., Jarra, W., Carret, C., Preiser, P., and Langhorne, J. (2009). Rapid changes in transcription profiles of the Plasmodium yoelii yir multigene family in clonal populations: lack of epigenetic memory? PLoS ONE 4:e4285. doi: 10.1371/journal.pone.0004285
Cunningham, D., Lawton, J., Jarra, W., Preiser, P., and Langhorne, J. (2010). The pir multigene family of Plasmodium: antigenic variation and beyond. Mol. Biochem. Parasitol. 170, 65–73. doi: 10.1016/j.molbiopara.2009.12.010
Damena, D., Denis, A., Golassa, L., and Chimusa, E. R. (2019). Genome-wide association studies of severe P. falciparum malaria susceptibility: progress, pitfalls and prospects. BMC Med Genomics 12:120. doi: 10.1186/s12920-019-0564-x
Daub, J. T., Hofer, T., Cutivet, E., Dupanloup, I., Quintana-Murci, L., Robinson-Rechavi, M., et al. (2013). Evidence for polygenic adaptation to pathogens in the human genome. Mol. Biol. Evol. 30, 1544–1558. doi: 10.1093/molbev/mst080
Day, K. P., Artzy-Randrup, Y., Tiedje, K. E., Rougeron, V., Chen, D. S., Rask, T. S., et al. (2017). Evidence of strain structure in Plasmodium falciparum var gene repertoires in children from Gabon, West Africa. Proc. Natl. Acad. Sci. U.S.A. 114, E4103–E4111. doi: 10.1073/pnas.1613018114
de Jong, R. M., Tebeje, S. K., Meerstein-Kessel, L., Tadesse, F. G., Jore, M. M., Stone, W., et al. (2020). Immunity against sexual stage Plasmodium falciparum and Plasmodium vivax parasites. Immunol. Rev 293, 190–215. doi: 10.1111/imr.12828
Deitsch, K. W., and Dzikowski, R. (2017). Variant gene expression and antigenic variation by malaria parasites. Annu. Rev. Microbiol. 71, 625–641. doi: 10.1146/annurev-micro-090816-093841
Dewasurendra, R. L., Jeffreys, A., Gunawardena, S. A., Chandrasekharan, N. V., Rockett, K., Kwiatkowski, D., et al. (2018). Host genetic polymorphisms and serological response against malaria in a selected population in Sri Lanka. Malar. J. 17:473. doi: 10.1186/s12936-018-2622-9
Draper, S. J., Sack, B. K., King, C. R., Nielsen, C. M., Rayner, J. C., Higgins, M. K., et al. (2018). Malaria vaccines: recent advances and new horizons. Cell Host Microbe 24, 43–56. doi: 10.1016/j.chom.2018.06.008
Drew, D. R., Hodder, A. N., Wilson, D. W., Foley, M., Mueller, I., Siba, P. M., et al. (2012). Defining the antigenic diversity of Plasmodium falciparum apical membrane antigen 1 and the requirements for a multi-allele vaccine against malaria. PLoS ONE 7:e51023. doi: 10.1371/journal.pone.0051023
Duval, L., and Ariey, F. (2012). Ape Plasmodium parasites as a source of human outbreaks. Clin. Microbiol. Infect. 18, 528–532. doi: 10.1111/j.1469-0691.2012.03825.x
Duval, L., Fourment, M., Nerrienet, E., Rousset, D., Sadeuh, S. A., Goodman, S. M., et al. (2010). African apes as reservoirs of Plasmodium falciparum and the origin and diversification of the Laverania subgenus. Proc. Natl. Acad. Sci. U.S.A. 107, 10561–10566. doi: 10.1073/pnas.1005435107
Early, A. M., Lievens, M., MacInnis, B. L., Ockenhouse, C. F., Volkman, S. K., Adjei, S., et al. (2018). Host-mediated selection impacts the diversity of Plasmodium falciparum antigens within infections. Nat. Commun. 9:1381. doi: 10.1038/s41467-018-03807-7
Ebel, E. R., Telis, N., Venkataram, S., Petrov, D. A., and Enard, D. (2017). High rate of adaptation of mammalian proteins that interact with Plasmodium and related parasites. PLoS Genet. 13:e1007023. doi: 10.1371/journal.pgen.1007023
Escalante, A. A., Cornejo, O. E., Freeland, D. E., Poe, A. C., Durrego, E., Collins, W. E., et al. (2005). A monkey's tale: the origin of Plasmodium vivax as a human malaria parasite. Proc. Natl. Acad. Sci. U.S.A. 102, 1980–1985. doi: 10.1073/pnas.0409652102
Escalante, A. A., and Pacheco, M. A. (2019). Malaria molecular epidemiology: an evolutionary genetics perspective. Microbiol. Spectr. 7, 1–17. doi: 10.1128/microbiolspec.AME-0010-2019
Fecchio, A., Chagas, C. R. F., Bell, J. A., and Kirchgatter, K. (2020). Evolutionary ecology, taxonomy, and systematics of avian malaria and related parasites. Acta Trop. 204:105364. doi: 10.1016/j.actatropica.2020.105364
Felger, I., Maire, M., Bretscher, M. T., Falk, N., Tiaden, A., Sama, W., et al. (2012). The dynamics of natural Plasmodium falciparum infections. PLoS ONE 7:e45542. doi: 10.1371/journal.pone.0045542
Fidock, D. A., Nomura, T., Talley, A. K., Cooper, R. A., Dzekunov, S. M., Ferdig, M. T., et al. (2000). Mutations in the P. falciparum digestive vacuole transmembrane protein PfCRT and evidence for their role in chloroquine resistance. Mol. Cell 6, 861–871. doi: 10.1016/S1097-2765(05)00077-8
Fish, E. N. (2008). The X-files in immunity: sex-based differences predispose immune responses. Nat. Rev. Immunol. 8, 737–744. doi: 10.1038/nri2394
Fluck, C., Smith, T., Beck, H. P., Irion, A., Betuela, I., Alpers, M. P., et al. (2004). Strain-specific humoral response to a polymorphic malaria vaccine. Infect. Immun. 72, 6300–6305. doi: 10.1128/IAI.72.11.6300-6305.2004
Frank, M., Kirkman, L., Costantini, D., Sanyal, S., Lavazec, C., Templeton, T. J., et al. (2008). Frequent recombination events generate diversity within the multi-copy variant antigen gene families of Plasmodium falciparum. Int. J. Parasitol. 38, 1099–1109. doi: 10.1016/j.ijpara.2008.01.010
Frimpong, A., Kusi, K. A., Ofori, M. F., and Ndifon, W. (2018). Novel strategies for malaria vaccine design. Front. Immunol. 9:2769. doi: 10.3389/fimmu.2018.02769
Galaway, F., Yu, R., Constantinou, A., Prugnolle, F., and Wright, G. J. (2019). Resurrection of the ancestral RH5 invasion ligand provides a molecular explanation for the origin of P. falciparum malaria in humans. PLOS Biol. 17:e3000490. doi: 10.1371/journal.pbio.3000490
Galen, S. C., Borner, J., Martinsen, E. S., Schaer, J., Austin, C. C., West, C. J., et al. (2018). The polyphyly of Plasmodium: comprehensive phylogenetic analyses of the malaria parasites (order Haemosporida) reveal widespread taxonomic conflict. R Soc. Open Sci. 5:171780. doi: 10.1098/rsos.171780
Galinski, M. R. (2019). Functional genomics of simian malaria parasites and host-parasite interactions. Brief. Funct. Genomics 18, 270–280. doi: 10.1093/bfgp/elz013
Galinski, M. R., Lapp, S. A., Peterson, M. S., Ay, F., Joyner, C. J., KG, L. E. R., et al. (2018). Plasmodium knowlesi: a superb in vivo nonhuman primate model of antigenic variation in malaria. Parasitology 145, 85–100. doi: 10.1017/S0031182017001135
Gardner, M. J., Hall, N., Fung, E., White, O., Berriman, M., Hyman, R. W., et al. (2002). Genome sequence of the human malaria parasite Plasmodium falciparum. Nature 419, 498–511. doi: 10.1038/nature01097
Garrido-Cardenas, J. A., Gonzalez-Ceron, L., Manzano-Agugliaro, F., and Mesa-Valle, C. (2019). Plasmodium genomics: an approach for learning about and ending human malaria. Parasitol. Res. 118, 1–27. doi: 10.1007/s00436-018-6127-9
Georgiadou, A., Lee, H. J., Walther, M., van Beek, A. E., Fitriani, F., Wouters, D., et al. (2019). Modelling pathogen load dynamics to elucidate mechanistic determinants of host-Plasmodium falciparum interactions. Nat. Microbiol. 4, 1592–1602. doi: 10.1038/s41564-019-0474-x
Goheen, M. M., Campino, S., and Cerami, C. (2017). The role of the red blood cell in host defence against falciparum malaria: an expanding repertoire of evolutionary alterations. Br. J. Haematol. 179, 543–556. doi: 10.1111/bjh.14886
Grajcarek, J., Monlong, J., Nishinaka-Arai, Y., Nakamura, M., Nagai, M., Matsuo, S., et al. (2019). Genome-wide microhomologies enable precise template-free editing of biologically relevant deletion mutations. Nat. Commun. 10:4856. doi: 10.1038/s41467-019-12829-8
Graves, S. F., Kouriba, B., Diarra, I., Daou, M., Niangaly, A., Coulibaly, D., et al. (2016). Strain-specific Plasmodium falciparum multifunctional CD4(+) T cell cytokine expression in Malian children immunized with the FMP2.1/AS02A vaccine candidate. Vaccine 34, 2546–2555. doi: 10.1016/j.vaccine.2016.04.019
Greenwood, J. M., Ezquerra, A. L., Behrens, S., Branca, A., and Mallet, L. (2016). Current analysis of host-parasite interactions with a focus on next generation sequencing data. Zoology 119, 298–306. doi: 10.1016/j.zool.2016.06.010
Halstead, S. B., Mahalingam, S., Marovich, M. A., Ubol, S., and Mosser, D. M. (2010). Intrinsic antibody-dependent enhancement of microbial infection in macrophages: disease regulation by immune complexes. Lancet Infect. Dis. 10, 712–722. doi: 10.1016/S1473-3099(10)70166-3
Hedrick, P. W. (2012). Resistance to malaria in humans: the impact of strong, recent selection. Malar. J. 11:349. doi: 10.1186/1475-2875-11-349
Howard, R. J., Barnwell, J. W., Rock, E. P., Neequaye, J., Ofori-Adjei, D., Maloy, W. L., et al. (1988). Two approximately 300 kilodalton Plasmodium falciparum proteins at the surface membrane of infected erythrocytes. Mol. Biochem. Parasitol. 27, 207–223. doi: 10.1016/0166-6851(88)90040-0
Hviid, L., and Jensen, A. T. (2015). PfEMP1 - a parasite protein family of key importance in Plasmodium falciparum Malaria immunity and pathogenesis. Adv. Parasitol. 88, 51–84. doi: 10.1016/bs.apar.2015.02.004
Idaghdour, Y., Quinlan, J., Goulet, J. P., Berghout, J., Gbeha, E., Bruat, V., et al. (2012). Evidence for additive and interaction effects of host genotype and infection in malaria. Proc. Natl. Acad. Sci. U.S.A. 109, 16786–16793. doi: 10.1073/pnas.1204945109
Janssen, C. S., Barrett, M. P., Turner, C. M., and Phillips, R. S. (2002). A large gene family for putative variant antigens shared by human and rodent malaria parasites. Proc. Biol. Sci. 269, 431–436. doi: 10.1098/rspb.2001.1903
Jeffares, D. C., Pain, A., Berry, A., Cox, A. V., Stalker, J., Ingle, C. E., et al. (2007). Genome variation and evolution of the malaria parasite Plasmodium falciparum. Nat. Genet. 39, 120–125. doi: 10.1038/ng1931
Jemmely, N. Y., Niang, M., and Preiser, P. R. (2010). Small variant surface antigens and Plasmodium evasion of immunity. Fut. Microbiol. 5, 663–682. doi: 10.2217/fmb.10.21
Jensen, A. R., Adams, Y., and Hviid, L. (2020). Cerebral Plasmodium falciparum malaria: the role of PfEMP1 in its pathogenesis and immunity, and PfEMP1-based vaccines to prevent it. Immunol. Rev. 293, 230–252. doi: 10.1111/imr.12807
Josling, G. A., and Llinas, M. (2015). Sexual development in Plasmodium parasites: knowing when it's time to commit. Nat. Rev. Microbiol. 13, 573–587. doi: 10.1038/nrmicro3519
Joy, D. A., Feng, X., Mu, J., Furuya, T., Chotivanich, K., Krettli, A. U., et al. (2003). Early origin and recent expansion of Plasmodium falciparum. Science 300, 318–321. doi: 10.1126/science.1081449
Joy, D. A., Gonzalez-Ceron, L., Carlton, J. M., Gueye, A., Fay, M., McCutchan, T. F., et al. (2008). Local adaptation and vector-mediated population structure in Plasmodium vivax malaria. Mol. Biol. Evol. 25, 1245–1252. doi: 10.1093/molbev/msn073
Kariuki, S. N., and Williams, T. N. (2020). Human genetics and malaria resistance. Hum. Genet. 139, 801–811. doi: 10.1007/s00439-020-02142-6
Kirchner, S., Power, B. J., and Waters, A. P. (2016). Recent advances in malaria genomics and epigenomics. Genome Med. 8:92. doi: 10.1186/s13073-016-0343-7
Kirkman, L., and Deitsch, K. W. (2020). Vive la différence: exploiting the differences between rodent and human malarias. Trends Parasitol. 36, 504–511. doi: 10.1016/j.pt.2020.03.008
Kirkman, L. A., Lawrence, E. A., and Deitsch, K. W. (2014). Malaria parasites utilize both homologous recombination and alternative end joining pathways to maintain genome integrity. Nucleic Acids Res. 42, 370–379. doi: 10.1093/nar/gkt881
Klein, S. L., Shann, F., Moss, W. J., Benn, C. S., and Aaby, P. (2016). RTS,S malaria vaccine and increased mortality in girls. MBio 7:e00514–00516. doi: 10.1128/mBio.00514-16
Kobbe, R., Neuhoff, R., Marks, F., Adjei, S., Langefeld, I., von Reden, C., et al. (2006). Seasonal variation and high multiplicity of first Plasmodium falciparum infections in children from a holoendemic area in Ghana, West Africa. Trop. Med. Int. Health 11, 613–619. doi: 10.1111/j.1365-3156.2006.01618.x
Kooij, T. W., Carlton, J. M., Bidwell, S. L., Hall, N., Ramesar, J., Janse, C. J., et al. (2005). A Plasmodium whole-genome synteny map: indels and synteny breakpoints as foci for species-specific genes. PLoS Pathog. 1:e44. doi: 10.1371/journal.ppat.0010044
Kotepui, M., Kotepui, K. U., De Jesus Milanez, G., and Masangkay, F. R. (2020). Plasmodium spp. mixed infection leading to severe malaria: a systematic review and meta-analysis. Sci. Rep. 10:11068. doi: 10.1038/s41598-020-68082-3
Krief, S., Escalante, A. A., Pacheco, M. A., Mugisha, L., André, C., Halbwax, M., et al. (2010). On the diversity of malaria parasites in African apes and the origin of Plasmodium falciparum from bonobos. PLoS Pathog. 6:e1000765. doi: 10.1371/journal.ppat.1000765
Kwiatkowski, D. P. (2005). How malaria has affected the human genome and what human genetics can teach us about malaria. Am. J. Hum. Genet. 77, 171–192. doi: 10.1086/432519
Kwon, H., and Smith, R. C. (2019). Chemical depletion of phagocytic immune cells in Anopheles gambiae reveals dual roles of mosquito hemocytes in anti-Plasmodium immunity. Proc. Natl. Acad. Sci. U.S.A. 116, 14119–14128. doi: 10.1073/pnas.1900147116
Lacerda-Queiroz, N., Lima, O. C., Carneiro, C. M., Vilela, M. C., Teixeira, A. L., Teixeira-Carvalho, A., et al. (2011). Plasmodium berghei NK65 induces cerebral leukocyte recruitment in vivo: an intravital microscopic study. Acta Trop. 120, 31–39. doi: 10.1016/j.actatropica.2011.04.020
Lacerda-Queiroz, N., Riteau, N., Eastman, R. T., Bock, K. W., Orandle, M. S., Moore, I. N., et al. (2017). Mechanism of splenic cell death and host mortality in a Plasmodium yoelii malaria model. Sci. Rep. 7:10438. doi: 10.1038/s41598-017-10776-2
Larsen, M. D., Quintana, M. D. P., Ditlev, S. B., Bayarri-Olmos, R., Ofori, M. F., Hviid, L., et al. (2018). Evasion of classical complement pathway activation on Plasmodium falciparum-infected erythrocytes opsonized by PfEMP1-specific IgG. Front. Immunol. 9:3088. doi: 10.3389/fimmu.2018.03088
Laurens, M. B. (2018). The promise of a malaria vaccine-are we closer? Annu. Rev. Microbiol. 72, 273–292. doi: 10.1146/annurev-micro-090817-062427
Laurens, M. B., Kouriba, B., Bergmann-Leitner, E., Angov, E., Coulibaly, D., Diarra, I., et al. (2017). Strain-specific Plasmodium falciparum growth inhibition among Malian children immunized with a blood-stage malaria vaccine. PLoS ONE 12:e0173294. doi: 10.1371/journal.pone.0173294
Lauron, E. J., Aw Yeang, H. X., Taffner, S. M., and Sehgal, R. N. (2015). De novo assembly and transcriptome analysis of Plasmodium gallinaceum identifies the Rh5 interacting protein (ripr), and reveals a lack of EBL and RH gene family diversification. Malar. J. 14:296. doi: 10.1186/s12936-015-0814-0
Lawton, J., Brugat, T., Yan, Y. X., Reid, A. J., Bohme, U., Otto, T. D., et al. (2012). Characterization and gene expression analysis of the cir multi-gene family of Plasmodium chabaudi chabaudi (AS). BMC Genomics 13:125. doi: 10.1186/1471-2164-13-125
Lee, A. H., Symington, L. S., and Fidock, D. A. (2014). DNA repair mechanisms and their biological roles in the malaria parasite Plasmodium falciparum. Microbiol. Mol. Biol. Rev. 78, 469–486. doi: 10.1128/MMBR.00059-13
Lee, H. J., Georgiadou, A., Walther, M., Nwakanma, D., Stewart, L. B., Levin, M., et al. (2018). Integrated pathogen load and dual transcriptome analysis of systemic host-pathogen interactions in severe malaria. Sci. Transl. Med. 10:eaar3619. doi: 10.1126/scitranslmed.aar3619
Liu, W., Li, Y., Learn, G. H., Rudicell, R. S., Robertson, J. D., Keele, B. F., et al. (2010). Origin of the human malaria parasite Plasmodium falciparum in gorillas. Nature 467, 420–425. doi: 10.1038/nature09442
Liu, W., Sherrill-Mix, S., Learn, G. H., Scully, E. J., Li, Y., Avitto, A. N., et al. (2017). Wild bonobos host geographically restricted malaria parasites including a putative new Laverania species. Nat. Commun. 8:1635. doi: 10.1038/s41467-017-01798-5
Lover, A. A., and Coker, R. J. (2015). Do mixed infections matter? Assessing virulence of mixed-clone infections in experimental human and murine malaria. Infect. Genet. Evol. 36, 82–91. doi: 10.1016/j.meegid.2015.08.041
Luzolo, A. L., and Ngoyi, D. M. (2019). Cerebral malaria. Brain Res. Bull. 145, 53–58. doi: 10.1016/j.brainresbull.2019.01.010
Mahanta, A., Ganguli, P., Barah, P., Sarkar, R. R., Sarmah, N., Phukan, S., et al. (2018). Integrative approaches to understand the mastery in manipulation of host cytokine networks by protozoan parasites with emphasis on Plasmodium and Leishmania species. Front. Immunol. 9:296. doi: 10.3389/fimmu.2018.00296
Malaria Genomic Epidemiology Network. (2019). Insights into malaria susceptibility using genome-wide data on 17,000 individuals from Africa, Asia and Oceania. Nat. Commun. 10:5732. doi: 10.1038/s41467-019-13480-z
Mangano, V. D., and Modiano, D. (2014). An evolutionary perspective of how infection drives human genome diversity: the case of malaria. Curr. Opin. Immunol. 30, 39–47. doi: 10.1016/j.coi.2014.06.004
Martin, M. J., Rayner, J. C., Gagneux, P., Barnwell, J. W., and Varki, A. (2005). Evolution of human-chimpanzee differences in malaria susceptibility: relationship to human genetic loss of N-glycolylneuraminic acid. Proc. Natl. Acad. Sci. U. S.A. 102, 12819–12824. doi: 10.1073/pnas.0503819102
Martinelli, A., and Culleton, R. (2018). Non-human primate malaria parasites: out of the forest and into the laboratory. Parasitology 145, 41–54. doi: 10.1017/S0031182016001335
McCollum, A. M., Basco, L. K., Tahar, R., Udhayakumar, V., and Escalante, A. A. (2008). Hitchhiking and selective sweeps of Plasmodium falciparum sulfadoxine and pyrimethamine resistance alleles in a population from central Africa. Antimicrob. Agents Chemother. 52, 4089–4097. doi: 10.1128/AAC.00623-08
Mejias, A., Suarez, N. M., and Ramilo, O. (2014). Detecting specific infections in children through host responses: a paradigm shift. Curr. Opin. Infect. Dis. 27, 228–235. doi: 10.1097/QCO.0000000000000065
Molina-Cruz, A., Canepa, G. E., Kamath, N., Pavlovic, N. V., Mu, J., Ramphul, U. N., et al. (2015). Plasmodium evasion of mosquito immunity and global malaria transmission: The lock-and-key theory. Proc. Natl. Acad. Sci. U.S.A. 112, 15178–15183. doi: 10.1073/pnas.1520426112
Mu, J., Awadalla, P., Duan, J., McGee, K. M., Keebler, J., Seydel, K., et al. (2007). Genome-wide variation and identification of vaccine targets in the Plasmodium falciparum genome. Nat. Genet. 39, 126–130. doi: 10.1038/ng1924
Mu, J., Duan, J., Makova, K. D., Joy, D. A., Huynh, C. Q., Branch, O. H., et al. (2002). Chromosome-wide SNPs reveal an ancient origin for Plasmodium falciparum. Nature 418, 323–326. doi: 10.1038/nature00836
Muller, O., Tozan, Y., and Becher, H. (2015). RTS,S/AS01 malaria vaccine and child mortality. Lancet 386:1736. doi: 10.1016/S0140-6736(15)00694-7
Nacer, A., Movila, A., Sohet, F., Girgis, N. M., Gundra, U. M., Loke, P., et al. (2014). Experimental cerebral malaria pathogenesis–hemodynamics at the blood brain barrier. PLoS Pathog. 10:e1004528. doi: 10.1371/journal.ppat.1004528
Neafsey, D. E., Hartl, D. L., and Berriman, M. (2005). Evolution of noncoding and silent coding sites in the Plasmodium falciparum and Plasmodium reichenowi genomes. Mol. Biol. Evol. 22, 1621–1626. doi: 10.1093/molbev/msi154
Neafsey, D. E., Juraska, M., Bedford, T., Benkeser, D., Valim, C., Griggs, A., et al. (2015). Genetic diversity and protective efficacy of the RTS,S/AS01 malaria vaccine. N. Engl. J. Med. 373, 2025–2037. doi: 10.1056/NEJMoa1505819
Negro, F. (2020). Is antibody-dependent enhancement playing a role in COVID-19 pathogenesis? Swiss Med. Wkly. 150:w20249. doi: 10.4414/smw.2020.20249
O'Brien, J. D., Iqbal, Z., Wendler, J., and Amenga-Etego, L. (2016). Inferring strain mixture within clinical Plasmodium falciparum isolates from genomic sequence data. PLoS Comput. Biol. 12:e1004824. doi: 10.1371/journal.pcbi.1004824
Ochola, L. I., Tetteh, K. K., Stewart, L. B., Riitho, V., Marsh, K., and Conway, D. J. (2010). Allele frequency-based and polymorphism-versus-divergence indices of balancing selection in a new filtered set of polymorphic genes in Plasmodium falciparum. Mol. Biol. Evol. 27, 2344–2351. doi: 10.1093/molbev/msq119
Ollomo, B., Durand, P., Prugnolle, F., Douzery, E., Arnathau, C., Nkoghe, D., et al. (2009). A new malaria agent in African hominids. PLoS Pathog. 5:e1000446. doi: 10.1371/journal.ppat.1000446
Olotu, A., Fegan, G., Wambua, J., Nyangweso, G., Leach, A., Lievens, M., et al. (2016). Seven-year efficacy of RTS,S/AS01 malaria vaccine among young African children. N. Engl. J. Med. 374, 2519–2529. doi: 10.1056/NEJMoa1515257
Otsuki, H., Kaneko, O., Thongkukiatkul, A., Tachibana, M., Iriko, H., Takeo, S., et al. (2009). Single amino acid substitution in Plasmodium yoelii erythrocyte ligand determines its localization and controls parasite virulence. Proc. Natl. Acad. Sci. U.S.A. 106, 7167–7172. doi: 10.1073/pnas.0811313106
Otto, T. D., Bohme, U., Jackson, A. P., Hunt, M., Franke-Fayard, B., Hoeijmakers, W. A., et al. (2014a). A comprehensive evaluation of rodent malaria parasite genomes and gene expression. BMC Biol. 12:86. doi: 10.1186/s12915-014-0086-0
Otto, T. D., Gilabert, A., Crellen, T., Böhme, U., Arnathau, C., Sanders, M., et al. (2018). Genomes of all known members of a Plasmodium subgenus reveal paths to virulent human malaria. Nat. Microbiol. 3, 687–697. doi: 10.1038/s41564-018-0162-2
Otto, T. D., Rayner, J. C., Bohme, U., Pain, A., Spottiswoode, N., Sanders, M., et al. (2014b). Genome sequencing of chimpanzee malaria parasites reveals possible pathways of adaptation to human hosts. Nat. Commun. 5:4754. doi: 10.1038/ncomms5754
Ouattara, A., Mu, J., Takala-Harrison, S., Saye, R., Sagara, I., Dicko, A., et al. (2010). Lack of allele-specific efficacy of a bivalent AMA1 malaria vaccine. Malar. J. 9:175. doi: 10.1186/1475-2875-9-175
Pain, A., Bohme, U., Berry, A. E., Mungall, K., Finn, R. D., Jackson, A. P., et al. (2008). The genome of the simian and human malaria parasite Plasmodium knowlesi. Nature 455, 799–803. doi: 10.1038/nature07306
Pasini, E. M., Bohme, U., Rutledge, G. G., Voorberg-Van der Wel, A., Sanders, M., Berriman, M., et al. (2017). An improved Plasmodium cynomolgi genome assembly reveals an unexpected methyltransferase gene expansion. Wellcome Open Res. 2:42. doi: 10.12688/wellcomeopenres.11864.1
Pattaradilokrat, S., Cheesman, S. J., and Carter, R. (2008). Congenicity and genetic polymorphism in cloned lines derived from a single isolate of a rodent malaria parasite. Mol. Biochem. Parasitol. 157, 244–247. doi: 10.1016/j.molbiopara.2007.10.011
Pattaradilokrat, S., Culleton, R. L., Cheesman, S. J., and Carter, R. (2009). Gene encoding erythrocyte binding ligand linked to blood stage multiplication rate phenotype in Plasmodium yoelii yoelii. Proc. Natl. Acad. Sci. U.S.A. 106, 7161–7166. doi: 10.1073/pnas.0811430106
Peng, Y. C., Qi, Y., Zhang, C., Yao, X., Wu, J., Pattaradilokrat, S., et al. (2020). Plasmodium yoelii erythrocyte-binding-like protein modulates host cell membrane structure, immunity, and disease severity. MBio 11:e02995-19. doi: 10.1128/mBio.02995-19
Perkins, S. L., Sarkar, I. N., and Carter, R. (2007). The phylogeny of rodent malaria parasites: simultaneous analysis across three genomes. Infect. Genet. Evol. 7, 74–83. doi: 10.1016/j.meegid.2006.04.005
Peters, J., Fowler, E., Gatton, M., Chen, N., Saul, A., and Cheng, Q. (2002). High diversity and rapid changeover of expressed var genes during the acute phase of Plasmodium falciparum infections in human volunteers. Proc. Natl. Acad. Sci. U.S.A. 99, 10689–10694. doi: 10.1073/pnas.162349899
Peterson, D. S., Milhous, W. K., and Wellems, T. E. (1990). Molecular basis of differential resistance to cycloguanil and pyrimethamine in Plasmodium falciparum malaria. Proc. Natl. Acad. Sci. U.S.A. 87, 3018–3022. doi: 10.1073/pnas.87.8.3018
Plenderleith, L. J., Liu, W., Learn, G. H., Loy, D. E., Speede, S., Sanz, C. M., et al. (2019). Ancient introgression between two ape malaria parasite species. Genome Biol. Evol. 11, 3269–3274. doi: 10.1093/gbe/evz244
Porterfield, J. S. (1986). Antibody-dependent enhancement of viral infectivity. Adv. Virus Res. 31, 335–355. doi: 10.1016/S0065-3527(08)60268-7
Proto, W. R., Siegel, S. V., Dankwa, S., Liu, W., Kemp, A., Marsden, S., et al. (2019). Adaptation of Plasmodium falciparum to humans involved the loss of an ape-specific erythrocyte invasion ligand. Nat. Commun. 10:4512. doi: 10.1038/s41467-019-12294-3
Prugnolle, F., Durand, P., Neel, C., Ollomo, B., Ayala, F. J., Arnathau, C., et al. (2010). African great apes are natural hosts of multiple related malaria species, including Plasmodium falciparum. Proc. Natl. Acad. Sci. U.S.A. 107, 1458–1463. doi: 10.1073/pnas.0914440107
Quirk, S. J., Haraldsson, G., Erlendsdottir, H., Hjalmarsdottir, M. A., van Tonder, A. J., Hrafnkelsson, B., et al. (2018). Effect of vaccination on Pneumococci isolated from the nasopharynx of healthy children and the middle ear of children with otitis media in Iceland. J. Clin. Microbiol. 56:e01046-18. doi: 10.1128/JCM.01046-18
Ramiro, R. S., Reece, S. E., and Obbard, D. J. (2012). Molecular evolution and phylogenetics of rodent malaria parasites. BMC Evol. Biol. 12:219. doi: 10.1186/1471-2148-12-219
Ranjit, M. R., and Sharma, Y. D. (1999). Genetic polymorphism of falciparum malaria vaccine candidate antigen genes among field isolates in India. Am. J. Trop. Med. Hyg. 61, 103–108. doi: 10.4269/ajtmh.1999.61.103
Rask, T. S., Hansen, D. A., Theander, T. G., Gorm Pedersen, A., and Lavstsen, T. (2010). Plasmodium falciparum erythrocyte membrane protein 1 diversity in seven genomes–divide and conquer. PLoS Comput. Biol. 6:e1000933. doi: 10.1371/journal.pcbi.1000933
Ravenhall, M., Campino, S., Sepulveda, N., Manjurano, A., Nadjm, B., Mtove, G., et al. (2018). Novel genetic polymorphisms associated with severe malaria and under selective pressure in North-eastern Tanzania. PLoS Genet. 14:e1007172. doi: 10.1371/journal.pgen.1007172
Rayner, J. C., Liu, W., Peeters, M., Sharp, P. M., and Hahn, B. H. (2011). A plethora of Plasmodium species in wild apes: a source of human infection? Trends Parasitol. 27, 222–229. doi: 10.1016/j.pt.2011.01.006
Recker, M., Buckee, C. O., Serazin, A., Kyes, S., Pinches, R., Christodoulou, Z., et al. (2011). Antigenic variation in Plasmodium falciparum malaria involves a highly structured switching pattern. PLoS Pathog. 7:e1001306. doi: 10.1371/journal.ppat.1001306
Reichenow, E. (1920). Uber das vorkormmen der malaiaparasiten des menschen bei den Afrikanischen menschenaffen. Ventralbl. Bakteriol. Parasitenkd. Intekionskr. Hyg. Abt. 85, 207–216.
Reid, A. J. (2015). Large, rapidly evolving gene families are at the forefront of host-parasite interactions in Apicomplexa. Parasitology 142(Suppl. 1), S57–S70. doi: 10.1017/S0031182014001528
Rich, S. M., Leendertz, F. H., Xu, G., LeBreton, M., Djoko, C. F., Aminake, M. N., et al. (2009). The origin of malignant malaria. Proc. Natl. Acad. Sci. U.S.A. 106, 14902–14907. doi: 10.1073/pnas.0907740106
Rich, S. M., Licht, M. C., Hudson, R. R., and Ayala, F. J. (1998). Malaria's Eve: evidence of a recent population bottleneck throughout the world populations of Plasmodium falciparum. Proc. Natl. Acad. Sci. U.S.A. 95, 4425–4430. doi: 10.1073/pnas.95.8.4425
Riley, E. M. (1999). Is T-cell priming required for initiation of pathology in malaria infections? Immunol. Today 20, 228–233. doi: 10.1016/S0167-5699(99)01456-5
Roper, C., Pearce, R., Nair, S., Sharp, B., Nosten, F., and Anderson, T. (2004). Intercontinental spread of pyrimethamine-resistant malaria. Science 305:1124. doi: 10.1126/science.1098876
Rorick, M. M., Rask, T. S., Baskerville, E. B., Day, K. P., and Pascual, M. (2013). Homology blocks of Plasmodium falciparum var genes and clinically distinct forms of severe malaria in a local population. BMC Microbiol. 13:244. doi: 10.1186/1471-2180-13-244
Rts, S. C. T. P. (2015). Efficacy and safety of RTS,S/AS01 malaria vaccine with or without a booster dose in infants and children in Africa: final results of a phase 3, individually randomised, controlled trial. Lancet 386, 31–45. doi: 10.1016/S0140-6736(15)60721-8
Rts, S. C. T. P., Agnandji, S. T., Lell, B., Fernandes, J. F., Abossolo, B. P., Methogo, B. G., et al. (2012). A phase 3 trial of RTS,S/AS01 malaria vaccine in African infants. N. Engl. J. Med. 367, 2284–2295. doi: 10.1056/NEJMoa1208394
Rts, S. C. T. P., Agnandji, S. T., Lell, B., Soulanoudjingar, S. S., Fernandes, J. F., Abossolo, B. P., et al. (2011). First results of phase 3 trial of RTS,S/AS01 malaria vaccine in African children. N. Engl. J. Med. 365, 1863–1875. doi: 10.1056/NEJMoa1102287
Rutledge, G. G., Bohme, U., Sanders, M., Reid, A. J., Cotton, J. A., Maiga-Ascofare, O., et al. (2017). Plasmodium malariae and P. ovale genomes provide insights into malaria parasite evolution. Nature 542, 101–104. doi: 10.1038/nature21038
Ruybal-Pesantez, S., Tiedje, K. E., Tonkin-Hill, G., Rask, T. S., Kamya, M. R., Greenhouse, B., et al. (2017). Population genomics of virulence genes of Plasmodium falciparum in clinical isolates from Uganda. Sci. Rep. 7:11810. doi: 10.1038/s41598-017-11814-9
Salinas, N. D., Tang, W. K., and Tolia, N. H. (2019). Blood-stage malaria parasite antigens: structure, function, and vaccine potential. J. Mol. Biol. 431, 4259–4280. doi: 10.1016/j.jmb.2019.05.018
Schaer, J., Perkins, S. L., Decher, J., Leendertz, F. H., Fahr, J., Weber, N., et al. (2013). High diversity of West African bat malaria parasites and a tight link with rodent Plasmodium taxa. Proc. Natl. Acad. Sci. U.S.A. 110, 17415–17419. doi: 10.1073/pnas.1311016110
Schall, J. J. (1982). Lizards infected with malaria: physiological and behavioral consequences. Science 217, 1057–1059. doi: 10.1126/science.7112113
Shutler, D., Reece, S. E., Mullie, A., Billingsley, P. F., and Read, A. F. (2005). Rodent malaria parasites Plasmodium chabaudi and P. vinckei do not increase their rates of gametocytogenesis in response to mosquito probing. Proc. Biol. Sci. 272, 2397–2402. doi: 10.1098/rspb.2005.3232
Siao, M. C., Borner, J., Perkins, S. L., Deitsch, K. W., and Kirkman, L. A. (2020). Evolution of host specificity by malaria parasites through altered mechanisms controlling genome maintenance. mBio 11:e03272-19. doi: 10.1128/mBio.03272-19
Smith, J. D., Chitnis, C. E., Craig, A. G., Roberts, D. J., Hudson-Taylor, D. E., Peterson, D. S., et al. (1995). Switches in expression of Plasmodium falciparum var genes correlate with changes in antigenic and cytoadherent phenotypes of infected erythrocytes. Cell 82, 101–110. doi: 10.1016/0092-8674(95)90056-X
Smith, J. D., Rowe, J. A., Higgins, M. K., and Lavstsen, T. (2013). Malaria's deadly grip: cytoadhesion of Plasmodium falciparum-infected erythrocytes. Cell. Microbiol. 15, 1976–1983. doi: 10.1111/cmi.12183
Spence, P. J., Jarra, W., Levy, P., Reid, A. J., Chappell, L., Brugat, T., et al. (2013). Vector transmission regulates immune control of Plasmodium virulence. Nature 498, 228–231. doi: 10.1038/nature12231
Stone, W., Bousema, T., Sauerwein, R., and Drakeley, C. (2019). Two-faced immunity? The evidence for antibody enhancement of malaria transmission. Trends Parasitol. 35, 140–153. doi: 10.1016/j.pt.2018.11.003
Stone, W. J. R., Campo, J. J., Ouedraogo, A. L., Meerstein-Kessel, L., Morlais, I., Da, D., et al. (2018). Unravelling the immune signature of Plasmodium falciparum transmission-reducing immunity. Nat. Commun. 9:558. doi: 10.1038/s41467-017-02646-2
Su, X.-z., Mu, J., and Joy, D. A. (2003). The “Malaria's Eve” hypothesis and the debate concerning the origin of the human malaria parasite Plasmodium falciparum. Microbes Infect. 5, 891–896. doi: 10.1016/S1286-4579(03)00173-4
Su, X. Z., Heatwole, V. M., Wertheimer, S. P., Guinet, F., Herrfeldt, J. A., Peterson, D. S., et al. (1995). The large diverse gene family var encodes proteins involved in cytoadherence and antigenic variation of Plasmodium falciparum-infected erythrocytes. Cell 82, 89–100. doi: 10.1016/0092-8674(95)90055-1
Su, X. Z., Lane, K. D., Xia, L., Sa, J. M., and Wellems, T. E. (2019). Plasmodium genomics and genetics: new insights into malaria pathogenesis, drug resistance, epidemiology, and evolution. Clin. Microbiol. Rev. 32:e00019-19. doi: 10.1128/CMR.00019-19
Sundararaman, S. A., Plenderleith, L. J., Liu, W., Loy, D. E., Learn, G. H., Li, Y., et al. (2016). Genomes of cryptic chimpanzee Plasmodium species reveal key evolutionary events leading to human malaria. Nat. Commun. 7:11078. doi: 10.1038/ncomms11078
Sutherland, C. J., Tanomsing, N., Nolder, D., Oguike, M., Jennison, C., Pukrittayakamee, S., et al. (2010). Two nonrecombining sympatric forms of the human malaria parasite Plasmodium ovale occur globally. J. Infect. Dis. 201, 1544–1550. doi: 10.1086/652240
Swardson-Olver, C. J., Dawson, T. C., Burnett, R. C., Peiper, S. C., Maeda, N., and Avery, A. C. (2002). Plasmodium yoelii uses the murine Duffy antigen receptor for chemokines as a receptor for normocyte invasion and an alternative receptor for reticulocyte invasion. Blood 99, 2677–2684. doi: 10.1182/blood.V99.8.2677
Sypniewska, P., Duda, J. F., Locatelli, I., Althaus, C. R., Althaus, F., and Genton, B. (2017). Clinical and laboratory predictors of death in African children with features of severe malaria: a systematic review and meta-analysis. BMC Med. 15:147. doi: 10.1186/s12916-017-0906-5
Tanabe, K., Zollner, G., Vaughan, J. A., Sattabongkot, J., Khuntirat, B., Honma, H., et al. (2015). Plasmodium falciparum: genetic diversity and complexity of infections in an isolated village in western Thailand. Parasitol. Int. 64, 260–266. doi: 10.1016/j.parint.2013.09.011
Tang, Y., Gupta, A., Garimalla, S., Ma, H. C. E., Galinski, M. R., Styczynski, M. P., et al. (2018). Metabolic modeling helps interpret transcriptomic changes during malaria. Biochim. Biophys. Acta Mol. Basis Dis. 1864, 2329–2340. doi: 10.1016/j.bbadis.2017.10.023
Tavares, J., Formaglio, P., Thiberge, S., Mordelet, E., van Rooijen, N., Medvinsky, A., et al. (2013). Role of host cell traversal by the malaria sporozoite during liver infection. J. Exp. Med. 210, 905–915. doi: 10.1084/jem.20121130
Templeton, T. J., Asada, M., Jiratanh, M., Ishikawa, S. A., Tiawsirisup, S., Sivakumar, T., et al. (2016). Ungulate malaria parasites. Sci. Rep. 6:23230. doi: 10.1038/srep23230
Tetteh, K. K., Stewart, L. B., Ochola, L. I., Amambua-Ngwa, A., Thomas, A. W., Marsh, K., et al. (2009). Prospective identification of malaria parasite genes under balancing selection. PLoS ONE 4:e5568. doi: 10.1371/journal.pone.0005568
Thiam, A., Sanka, M., Ndiaye Diallo, R., Torres, M., Mbengue, B., Nunez, N. F., et al. (2019). Gene expression profiling in blood from cerebral malaria patients and mild malaria patients living in Senegal. BMC Med. Genomics 12:148. doi: 10.1186/s12920-019-0599-z
Tishkoff, S. A., and Williams, S. M. (2002). Genetic analysis of African populations: human evolution and complex disease. Nat. Rev. Genet. 3, 611–621. doi: 10.1038/nrg865
Tonkin-Hill, G. Q., Trianty, L., Noviyanti, R., Nguyen, H. H. T., Sebayang, B. F., Lampah, D. A., et al. (2018). The Plasmodium falciparum transcriptome in severe malaria reveals altered expression of genes involved in important processes including surface antigen-encoding var genes. PLoS Biol. 16:e2004328. doi: 10.1371/journal.pbio.2004328
Trampuz, A., Jereb, M., Muzlovic, I., and Prabhu, R. M. (2003). Clinical review: severe malaria. Crit Care 7, 315–323. doi: 10.1186/cc2183
Trimnell, A. R., Kraemer, S. M., Mukherjee, S., Phippard, D. J., Janes, J. H., Flamoe, E., et al. (2006). Global genetic diversity and evolution of var genes associated with placental and severe childhood malaria. Mol. Biochem. Parasitol. 148, 169–180. doi: 10.1016/j.molbiopara.2006.03.012
Verra, F., Mangano, V. D., and Modiano, D. (2009). Genetics of susceptibility to Plasmodium falciparum: from classical malaria resistance genes towards genome-wide association studies. Parasite Immunol. 31, 234–253. doi: 10.1111/j.1365-3024.2009.01106.x
Videvall, E. (2019). Genomic advances in avian malaria research. Trends Parasitol. 35, 254–266. doi: 10.1016/j.pt.2018.12.005
Vijayan, A., and Chitnis, C. E. (2019). Development of blood stage malaria vaccines. Methods Mol. Biol. 2013, 199–218. doi: 10.1007/978-1-4939-9550-9_15
Volkman, S. K., Barry, A. E., Lyons, E. J., Nielsen, K. M., Thomas, S. M., Choi, M., et al. (2001). Recent origin of Plasmodium falciparum from a single progenitor. Science 293, 482–484. doi: 10.1126/science.1059878
Volkman, S. K., Sabeti, P. C., deCaprio, D., Neafsey, D. E., Schaffner, S. F., Milner, D. A., et al. (2007). A genome-wide map of diversity in Plasmodium falciparum. Nat. Genet. 39, 113–119. doi: 10.1038/ng1930
Waggoner, J. J., Katzelnick, L. C., Burger-Calderon, R., Gallini, J., Moore, R. H., Kuan, G., et al. (2020). Antibody-dependent enhancement of severe disease is mediated by serum viral load in pediatric dengue virus infections. J. Infect. Dis. 6:29201. doi: 10.1093/infdis/jiz618
Wahlgren, M., Goel, S., and Akhouri, R. R. (2017). Variant surface antigens of Plasmodium falciparum and their roles in severe malaria. Nat. Rev. Microbiol 15, 479–491. doi: 10.1038/nrmicro.2017.47
Wan, Y., Shang, J., Sun, S., Tai, W., Chen, J., Geng, Q., et al. (2020). Molecular mechanism for antibody-dependent enhancement of coronavirus entry. J. Virol. 94:e02015-19. doi: 10.1128/JVI.02015-19
Wanaguru, M., Liu, W., Hahn, B. H., Rayner, J. C., and Wright, G. J. (2013). RH5-Basigin interaction plays a major role in the host tropism of Plasmodium falciparum. Proc. Natl. Acad. Sci. U.S.A. 110, 20735–20740. doi: 10.1073/pnas.1320771110
Wang-Lin, S. X., Olson, R., Beanan, J. M., MacDonald, U., Russo, T. A., and Balthasar, J. P. (2019). Antibody dependent enhancement of Acinetobacter baumannii infection in a mouse pneumonia model. J. Pharmacol. Exp. Ther. 368, 475–489. doi: 10.1124/jpet.118.253617
Weedall, G. D., and Conway, D. J. (2010). Detecting signatures of balancing selection to identify targets of anti-parasite immunity. Trends Parasitol. 26, 363–369. doi: 10.1016/j.pt.2010.04.002
WHO (2019). World Malaria Report 2019. Available online at: https://www.who.int/news-room/feature-stories/detail/world-malaria-report-2019.
Winarski, K. L., Tang, J., Klenow, L., Lee, J., Coyle, E. M., Manischewitz, J., et al. (2019). Antibody-dependent enhancement of influenza disease promoted by increase in hemagglutinin stem flexibility and virus fusion kinetics. Proc. Natl. Acad. Sci. U.S.A. 116, 15194–15199. doi: 10.1073/pnas.1821317116
Woolley, I. J., Wood, E. M., Sramkoski, R. M., Zimmerman, P. A., Miller, J. P., and Kazura, J. W. (2005). Expression of Duffy antigen receptor for chemokines during reticulocyte maturation: using a CD71 flow cytometric technique to identify reticulocytes. Immunohematology 21, 15–20.
Wootton, J. C., Feng, X., Ferdig, M. T., Cooper, R. A., Mu, J., Baruch, D. I., et al. (2002). Genetic diversity and chloroquine selective sweeps in Plasmodium falciparum. Nature 418, 320–323. doi: 10.1038/nature00813
World Health Organization (2018). Malaria vaccine: WHO position paper, January 2016 - recommendations. Vaccine 36, 3576–3577. doi: 10.1016/j.vaccine.2016.10.047
Wu, J., Cai, B., Sun, W., Huang, R., Liu, X., Lin, M., et al. (2015). Genome-wide analysis of host-Plasmodium yoelii interactions reveals regulators of the type i interferon response. Cell Rep. 12, 661–672. doi: 10.1016/j.celrep.2015.06.058
Wu, J., Tian, L., Yu, X., Pattaradilokrat, S., Li, J., Wang, M., et al. (2014). Strain-specific innate immune signaling pathways determine malaria parasitemia dynamics and host mortality. Proc. Natl. Acad. Sci. U.S.A. 111, E511–E520. doi: 10.1073/pnas.1316467111
Wu, J., Xia, L., Yao, X., Yu, X., Tumas, K. C., Sun, W., et al. (2020). The E3 ubiquitin ligase MARCH1 regulates antimalaria immunity through interferon signaling and T cell activation. Proc. Natl. Acad. Sci. U.S.A. 117, 16567–16578. doi: 10.1073/pnas.2004332117
Wuchty, S. (2011). Computational prediction of host-parasite protein interactions between P. falciparum and H. sapiens. PLoS ONE 6:e26960. doi: 10.1371/journal.pone.0026960
Xia, L., Wu, J., Pattaradilokrat, S., Tumas, K., He, X., Peng, Y. C., et al. (2018). Detection of host pathways universally inhibited after Plasmodium yoelii infection for immune intervention. Sci. Rep. 8:15280. doi: 10.1038/s41598-018-33599-1
Xu, R., Liu, Y., Fan, R., Liang, R., Yue, L., Liu, S., et al. (2019). Generation and functional characterisation of Plasmodium yoelii csp deletion mutants using a microhomology-based CRISPR/Cas9 method. Int. J. Parasitol. 49, 705–714. doi: 10.1016/j.ijpara.2019.04.003
Yam, X. Y., Brugat, T., Siau, A., Lawton, J., Wong, D. S., Farah, A., et al. (2016). Characterization of the Plasmodium interspersed repeats (PIR) proteins of Plasmodium chabaudi indicates functional diversity. Sci. Rep. 6:23449. doi: 10.1038/srep23449
Yam, X. Y., and Preiser, P. R. (2017). Host immune evasion strategies of malaria blood stage parasite. Mol. Biosyst. 13, 2498–2508. doi: 10.1039/C7MB00502D
Zhang, X., Alexander, N., Leonardi, I., Mason, C., Kirkman, L. A., and Deitsch, K. W. (2019). Rapid antigen diversification through mitotic recombination in the human malaria parasite Plasmodium falciparum. PLoS Biol. 17:e3000271. doi: 10.1371/journal.pbio.3000271
Zhu, S. J., Hendry, J. A., Almagro-Garcia, J., Pearson, R. D., Amato, R., Miles, A., et al. (2019). The origins and relatedness structure of mixed infections vary with local prevalence of P. falciparum malaria. Elife 8:e40845. doi: 10.7554/eLife.40845
Glossary
cyt b, cytochrome oxidase b gene; ADE, antibody dependent enhancement; AMA1, apical membrane antigen 1; CM, cerebral malaria; CSP, circumsporozoite protein; DBL, Duffy binding-like; CMAH, cytidine monophospho-N-acetylneuraminic acid hydroxylase; dhfr-ts, dihydrofolate reductase-thymidylate synthase gene; DTP, diphtheria-tetanus-pertussis; DARC, Duffy antigen/receptor for chemokines; EPCR, endothelial protein C receptor; EBL, erythrocyte binding-like; ECM, experimental cerebral malaria; GPCR, G-protein coupled receptor; GWAS, genome-wide association studies; GPLS, genome-wide patterns of LOD score; GYPA, glycophorin A; HBsAg, hepatitis B virus surface antigen; HTMV, high-titer measle vaccine; HKA, Hudson-Kreitman-Aguade; HLA, human leukocyte antigen; iRBC, infected erythrocyte; IntRh5, introgressed Rh5 gene; LGS, linkage group selection; MSP1, merozoite surface protein 1; MSP2, merozoite surface protein 2; MMEJ, microhomology-mediated end joining; MS, microsatellite; MM, mild malaria; OR, olfactory receptor; PfCRT, P. falciparum chloroquine resistance transporter; PfDHFR-TS, P. falciparum dihydrofolate reductase-thymidylate synthase; PfDHPS, P. falciparum dihydropteroate synthase; PfEbl165, P. falciparum EBL165 gene; PfEMP1, P. falciparum erythrocyte membrane protein 1; N67, Plasmodium y. nigeriensis N67; N67C, Plasmodium y. nigeriensis N67C; PyEBL, Plasmodium yoelii erythrocyte binding-like protein; PAMPs, pathogen-associated molecular patterns; PS, phosphatidylserine; PPIPs, Piroplasm-interacting proteins; pir, Plasmodium interspersed repeat; PCV, pneumococcal conjugate vaccine; RBC, red blood cell; rif, repetitive interspersed; RH5, reticulocyte-binding-like homologous protein 5; RLR, retinoic acid-inducible gene-I-like receptor; SICA, schizont infected cell agglutination; NCM, severe non-cerebral malaria; SNP, single nucleotide polymorphism; stevor, subtelomeric variant open reading frame; SPR, surface plasmon resonance; TjD, Tajima's D; CyRPA, cysteine-rich protective antigen; TLR, Toll-like receptors; Ts-eQTL, trans-species expression quantitative trait locus; VE, vaccine efficacy; WHO, World Health Organization.
Keywords: Plasmodium, genome diversity, population, vaccine, immunity, selection
Citation: Su X-z, Zhang C and Joy DA (2020) Host-Malaria Parasite Interactions and Impacts on Mutual Evolution. Front. Cell. Infect. Microbiol. 10:587933. doi: 10.3389/fcimb.2020.587933
Received: 11 August 2020; Accepted: 22 September 2020;
Published: 27 October 2020.
Edited by:
Yufeng Wang, University of Texas at San Antonio, United StatesReviewed by:
Fang Huang, National Institute of Parasitic Diseases, ChinaLiwang Cui, University of South Florida, United States
Copyright © 2020 Su, Zhang and Joy. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xin-zhuan Su, eHN1JiN4MDAwNDA7bmlhaWQubmloLmdvdg==; Deirdre A. Joy, ZGpveSYjeDAwMDQwO25pYWlkLm5paC5nb3Y=