Recent advances in malaria genomics and epigenomics

Parasite development and pathogenesis

Within their mammalian host and mosquito vector Plasmodium parasites complete a remarkable lifecycle, alternating between asexual and sexual replication (Fig. 2). Throughout the Plasmodium lifecycle, regulation of gene expression is orchestrated by a variety of mechanisms, including epigenetic, transcriptional, post-transcriptional, and translational control of gene expression. Owing to the absence of most canonical eukaryotic transcription factors in the Plasmodium genome [2], epigenetic control has long been recognized to play an important role in gene expression regulation.

Epigenetics lies at the very heart of gene expression, regulating access of the transcriptional machinery to chromatin [20] via (1) post-translational modifications (PTMs) of histones, (2) nucleosome occupancy, and (3) global chromatin architecture. In the past decade, various histone PTMs have been identified throughout the Plasmodium lifecycle (reviewed in [21]) and the existing catalog of modifications in Pf was recently extended to 232 distinct PTMs, 88 unique to Plasmodium [22]. The majority of detected PTMs show dynamic changes across the intra-erythrocytic development cycle (IDC), likely mirroring changes within chromatin organization linked to its transcriptional status. Methylation and acetylation of N-terminal histone tails are by far the most studied regulatory PTMs, linked either to a transcriptionally active chromatin structure (that is, euchromatin) or to transcriptionally inert heterochromatin. In Pf, various genes encoding putative epigenetic modulators (that is, proteins catalyzing either the addition or removal of histone PTM marks) have been identified [23], but only a few have been subjected to more detailed investigation [24, 25]. Many of the histone modifiers are essential for Plasmodium development, making them a promising target for antimalarial drugs [26]. In Pf, conditional knockdown of HDA2, a histone lysine deacetylase (HDAC) catalyzing the removal of acetyl groups from acetylated histone 3 lysine 9 (H3K9ac), resulted in elevated H3K9ac levels in previously defined heterochromatin regions [27]. H3K9ac is an epigenetic mark associated with transcriptionally active euchromatin [28] and HDA2 depletion resulted in the transcriptional activation of genes located in heterochromatin regions, leading to impaired asexual growth and an increased gametocyte conversion [27]. Interestingly, genes found to be dysregulated by HDA2 knockdown are also known to be associated with HP1, a key epigenetic player binding to tri-methylated H3K9 (H3K9me3), linked to transcriptionally repressed chromatin. Strikingly, conditional knockdown of PfHP1 recapitulated, to a much greater extent, the phenotype observed in HDA2-knockdown mutants [29]. HP1 is believed to act as a recruitment platform for histone lysine methyltransferases (HKMTs), required for maintenance and spreading of H3K9me3 marks [30], which is consistent with the reduction of H3K9me3 observed in HP1 knockdown cells [29]. In addition, bromodomain protein 1 (BDP1) was found to bind to H3K9ac and H3K14ac marks within transcription start sites (TSSs) in Pf, among them predominantly invasion-related genes (Fig. 3a), and BDP1-knockdown parasites consistently failed to invade new erythrocytes. BDP1 also appears to act as a recruitment platform for other effector proteins such as BDP2 and members of the apicomplexan AP2 (ApiAP2) transcription factor family [31].

In addition to histone PTMs, nucleosome organization plays a critical role in gene expression regulation in Plasmodium. In general, heterochromatin is substantially enriched in nucleosomes compared with euchromatin [32] and active promoters and intergenic regions in Pf show markedly reduced nucleosome occupancy [33]. In addition, common transcript features such as TSSs, transcription termination sites, and splice donor/acceptor sites show clearly distinguishable nucleosome positioning in Pf [34], but previously described dynamic changes in nucleosome positioning [32] appeared to be mostly restricted to TSSs during the IDC [34]. Uniquely in Plasmodium spp., canonical histones in intergenic regions are replaced by histone variant H2A.Z [28], which, in concert with the apicomplexan-specific H2B.Z, establishes a H2A.Z/H2B.Z double-variant nucleosome subtype enriched at AT-rich promoter regions and correlates with open chromatin and active gene transcription [35].

Within the confined space of the nucleus, chromosomes are tightly packed into a three-dimensional structure. This three-dimensional architecture allows interaction between otherwise distant chromatin regions possessing regulatory function and facilitates contacts with other nuclear sub-compartments such as the nucleolus and the nuclear envelope [36]. Until recently, knowledge of the chromosome architecture and chromatin interactions in Plasmodium was mostly restricted to single genomic loci based on fluorescence in situ hybridization experiments [37]. However, recent advances in deep-sequencing technologies [38] have for the first time enabled the genome-wide profiling of chromosome interactions at kilobase resolution in Plasmodium [37, 39]. In contrast to other eukaryotic organisms, the Pf nucleus appears to lack clearly defined chromosome territories and chromatin interactions are mainly restricted to intra-chromosomal contacts showing a clear distance-related dependency [37, 39]. Inter-chromosomal interactions are mostly absent in Pf and restricted to centromeres, telomeres, ribosomal DNA (rDNA) loci, and internal as well as subtelomeric-localized var genes (further discussed in the next section). This observed clustering appears to coincide with transcriptional activity of each cluster. Interestingly, using three-dimensional chromatin modeling, the highly transcribed rDNA genes were proposed to be localized to the nuclear periphery, which was previously mainly associated with transcriptionally silenced heterochromatin [40], indicative of perinuclear transcriptionally active compartments [37].

Transcription itself is initiated through binding of the transcriptional machinery to promoter regions in the nucleus, resulting in the synthesis of pre-mRNA molecules, which, following extensive processing and nuclear export, leads to the accumulation of mature mRNAs in the parasite cytosol [41]. A recent study found evidence for stage-specific transcription initiation from distinct TSSs of otherwise identical transcriptional units, giving rise to developmentally regulated mRNA isoforms [42]. While most canonical eukaryotic transcription factors are absent from the Plasmodium genome [2], the ApiAP2 family of DNA-binding proteins comprises by far the largest group of transcription factors in malaria parasites [43]. A collection of ApiAP2 proteins is expressed throughout all stages of the IDC [44], while other ApiAP2 proteins are expressed outside the IDC [45–47]. ApiAP2s appear to be among the main drivers of developmental progress throughout most Plasmodium lifecycle stages and their disruption abolishes or greatly reduces parasite development [45, 46]. They bind in a sequence-specific fashion to motifs generally distributed upstream of open reading frames (ORFs) and individual AP2s may have widespread influence; PfAP2-O has been shown to bind upstream of >500 genes (roughly 10 % of the parasite ORFs), potentially influencing a wide range of cellular activities [48].

Through forward genetic screens and comparative genomics, ApiAP2-G was discovered to function as a conserved master regulator of sexual commitment in Pf and Pb. ApiAP2-G binds to a conserved 6/8-mer nucleotide motif enriched upstream of gametocyte-specific genes and ap2-g itself, leading to an autoregulatory feedback loop [49, 50] (Fig. 3b). ApiAP2-G2, another ApiAP2 family member, acts downstream of ApiAP2-G during sexual development, functioning as a transcriptional repressor blocking expression of genes required for asexual development and influencing gametocyte sex ratios [50, 51]. During the asexual IDC, ap2-g displays characteristics of epigenetically silenced heterochromatin, such as H3K9me3 marks, binding to HP1 and localization to the nuclear periphery (reviewed in [52]) (Fig. 3b). However, the previously mentioned knockdowns of both PfHDA2 and HP1 resulted in increased gametocyte conversion, likely as a direct consequence of the loss of H3K9me3 marks and H3K9 hyperacetylation leading to ap2-g transcriptional activation [27, 29]. This opens the possibility of a bet-hedging mechanism for sexual commitment in Plasmodium, regulating stochastic, low-level activation of ap2-g sensitive to environmental stimuli, as has been shown for several blood-stage-expressed genes [52, 53]. PTMs such as lysine acetylation are not restricted to histones and a recent study has demonstrated that the “acetylome” impacts >1000 proteins and intriguingly is highly enriched in the ApiAP2 transcription factor family [54, 55], although the functional consequences of these PTMs have yet to be established.

Following their synthesis, eukaryotic mRNAs are processed and are finally translated by the ribosomal machinery. Translation has long been a focus of malaria research, not only because it represents a promising target for antimalarial drugs but also for its potential regulatory features [56]. The lack of correlation between transcript and protein level observed throughout the Plasmodium lifecycle has fueled researchers’ interest in post-transcriptional and translational control for decades [57]. Many features of post-transcriptional/translational control in malaria parasites are similar to the mechanisms found in other eukaryotes [41]. However, the advent of ribosome profiling [58] has enabled in-depth genome-wide analysis of the Plasmodium translatome. Throughout the IDC, transcription and translation are tightly coupled and only 8 % (approximately 300 transcripts) of the transcriptome was found to be translationally regulated [59]. These genes were found to be involved in merozoite egress and invasion, and while transcript level peaked during the late stages of the IDC, maximal translation was observed during the early ring stage. This observation resembles a general feature of gene expression in Plasmodium, whereby for a set of genes transcription and translation are uncoupled and mRNA translation occurs during a later developmental time point when compared with maximal transcriptional activity, most notably in female gametocytes [46, 60–64]. This is especially true for genes required for developmental progression and provides the parasite with the capability of rapid and timely protein synthesis without the need for preceding de novo mRNA synthesis. Recently, PfALBA1, a member of the DNA/RNA-binding Alba protein family, was postulated to act as master regulator during the Pf IDC, controlling translation of invasion-related transcripts (Fig. 3a) as well as regulating mRNA homeostasis of approximately 100 transcripts in blood-stage parasites [65]. In contrast to findings by Caro and colleagues [59], an earlier study using polysome profiling found a discrepancy between steady-state mRNA level and polysome-associated mRNAs among 30 % of genes (1280 transcripts) during the Pf IDC, indicative of translationally controlled gene expression [66]. Additionally, the results of this study, as well the findings of others, suggest upstream ORF translation and stop codon read-through in Pf [67–69], but the genome-wide extent of such mechanisms in Plasmodium spp. remains controversial [59]. Hence, expansion of these studies to other parasite life stages, such as the gametocyte, where translational control is firmly established, would surely give further insights into the extent of translation regulation in Plasmodium.

In addition to canonical protein-coding mRNAs, a vast number of genes encoding different ncRNAs have been identified within the Plasmodium genome in recent years, which are believed to exert a variety of regulatory functions (reviewed in [70]). Circular RNAs (circRNAs) are among the newest members of the still expanding catalogue of existing ncRNAs in Plasmodium [17]. Host microRNAs (miRNAs) have been shown to regulate parasite translation [71], and circRNAs therefore might act as sponge for host miRNAs, a mechanism described in other organisms [72]. Recent studies have especially increased our knowledge of the role of ncRNAs in var gene regulation (discussed in the next section) but, nevertheless, the biological role of the vast majority of these ncRNA species remains unclear.

Immune evasion

In their attempts to occupy a diverse range of host environments, protozoan parasites of the Plasmodium genus have evolved a plethora of molecular mechanisms to evade the host adaptive immune response. The host immune response to Plasmodium infection is dependent upon both host and parasite genomics and the developmental stage and phenotype of the invading parasite [73–75]. In the best-studied example in Plasmodium, virulence of Pf is attributed largely to monoallelic expression of just one of approximately 60 var genes that encode variant copies of the surface antigen, Pf erythrocyte membrane protein 1 (PfEMP1). The ability to switch expression from one var gene to another enables the invading parasite to alternate between phenotypes of variable cytoadherent and immunogenic properties [76–78]. PfEMP1 proteins are expressed at parasite-induced knobs on the infected erythrocyte surface, which are electron-dense features comprising many parasite proteins anchored to the erythrocyte cytoskeleton. Failure to present PfEMP1 in such knob structures greatly reduces the ability of the infected erythrocyte to bind to its specific host receptor [79].

Pf var gene regulation is complex and includes mechanisms of gene regulation such as chromosomal organization and subnuclear compartmentalization [80, 81], endogenous var gene clustering and var promoter–intron pairing [82, 83], transcriptional gene silencing via exoribonuclease-mediated RNA degradation [84], histone variant exchange at var promoters [85, 86], the effect of trans antisense long non-coding RNAs (lncRNAs) [87], and the presence or absence of histone modifications and their associated histone-modifying enzymes [27, 29, 40, 87–92] (Fig. 3c). Interest in delineating these mechanisms has continued, and even grown, as more research in the post-genomic area has highlighted the important differential role of the 5′ upstream promoter families into which the var genes can be subdivided into five classes (upsA to upsE), which correlate closely with the severity of malaria infection in the human host [93–98]. Pf var gene promoters are also essential components of the gene silencing mechanism and monoallelic expression. The upsC var promoter in particular is necessary to maintain chromosome-internal var genes in their silenced state and recently has been proposed to do so through the interaction of cis-acting MEE2-like sequence motifs and MEE2-interacting factors to reinforce var gene transcriptional repression [75, 83].

Monoallelic var gene transcription is also associated with the presence of H3K9me3 repressive marks at silent var gene loci (Fig. 3c). This histone modification is predicted but not proven to be imposed by the HKMT PfSET3 and is associated with perinuclear repressive centers and the binding of PfHP1, stimulating heterochromatin formation [40, 89, 90, 92]. Conditional disruption of one of these essential proteins, HP1, disrupts singular var gene expression and dysregulates antigenic variation [29]. In addition, conditional knockdown of PfHDA2 has been shown to result in a dramatic loss of monoallelic var gene expression [27]. This implicated PfHDA2 as an upstream regulator of HP1 binding as it facilitates the establishment of the H3K9me3 mark. The indispensable role of the dynamic histone lysine methylation of Plasmodium chromatin by histone lysine demethylases (HKDMs) and HKMTs in controlling the transcription of nearly all var genes has also been demonstrated. Knockout of the Pf hkmt gene encoding SET2/SETvs (vs, variant-silencing) resulted in reduced presence of the repressive H3K36me3 mark at the TSSs and intronic promoters of all var gene subtypes (Fig. 3c). Loss of this SETvs-dependent histone modification resulted in the loss of monoallelic var gene expression and expression of the entire var repertoire [98]. Furthermore, SETvs can directly interact with the C-terminal domain of RNA polymerase II, with SETvs disruption resulting in a loss of binding to RNA polymerase II and var gene switching [99].

Pf upsA-type var gene expression is also regulated by PfRNaseII, a chromatin-associated exoribonuclease. An inverse relationship exists between transcript levels of PfRNaseII and upsA-type var genes, with an increase in the latter corresponding to incidences of severe malaria in infected patients [84]. PfRNaseII is proposed to control upsA-type var gene transcription by marking TSSs and intronic promoter regions, degrading potential full-length transcripts to produce short-lived cryptic RNA molecules that are then further degraded by the exosome immediately upon expression (Fig. 3c). Disruption of the pfrnaseII gene resulted in loss of this degradation and the generation of full-length upsA var gene transcripts and intron-derived antisense lncRNA. These data illustrate the relationship between PfRNaseII and the control of monoallelic var gene transcription and suggest a correlation between lncRNA and var gene activation in Pf [84]. The role of lncRNAs in Pf var gene activation was again investigated in a study by Amit-Avraham and colleagues [87], which demonstrated dose-dependent transcriptional activation of var genes by overexpression of their individual antisense lncRNA transcripts. Disruption of antisense lncRNA expression by peptic nucleic acids resulted in downregulation of active var gene transcripts and induced var gene switching. The exact mechanism by which antisense lncRNAs act to promote the active transcription of a var gene is unknown. It has been postulated that antisense var transcripts may recruit chromatin-modifying enzymes that in turn would affect gene accessibility for the Pf transcriptional machinery. Antisense var gene lncRNAs would also contain a complementary sequence to var gene intronic insulator-like pairing elements that bind specific nuclear-binding proteins, therefore blocking the silencing activity of pairing elements by hybridization [87, 100].

The Plasmodium helical interspersed subtelomeric protein (PHIST) family of genes, which is unique to Pf, has also been implicated in the regulation of immune evasion as a result of its ability to bind to the intracellular acidic terminal segment of PfEMP1. Conditional knockdown of the essential PHIST protein PfE1605w reduced the capability of the infected host erythrocyte to adhere to the CD36 endothelial receptor, an important virulence feature of Pf. This study highlighted the importance not only of var genes and their controlled expression but also of other genes that are associated with anchoring PfEMP1 to the erythrocyte surface and creating the Plasmodium cytoadherence complex [101].

The list of regulatory mechanisms underlying var gene monoallelic expression is vast and much more may still be discovered in this area. However, immune evasion in the Plasmodium genus is not confined to Pf or var gene regulation. Indeed, var gene expression is exclusive to Pf, with much still to be gleaned in the areas of immune evasion in human malaria parasites such as P. vivax, P. knowlesi, Plasmodium ovale and Plasmodium malariae [13, 102–105]. In addition, PfEMP1 is just one of a number of variant surface antigens (VSAs) known to be expressed at the host erythrocyte surface upon infection with Pf, although it is the best characterized. Pf-infected erythrocytes also express VSAs of the multi-copy gene families of the proteins repetitive interspersed family (RIFIN), subtelomeric variable open reading frame (STEVOR), and Pf Maurer’s cleft 2 transmembrane (PfMC-2TM) [106]. The roles of these protein families in antigenic variation and pathology are generally poorly defined but are being elucidated; for example, RIFINs are implicated in the severity of Pf malaria in African children with blood group A. This tendency toward increased malaria pathogenicity is a result of their expression at the surface of infected host erythrocytes, from which they bind uninfected erythrocytes (preferentially, blood group A) to form rosette structures and mediate binding to the host microvasculature [107]. Thus, the combined roles of HP1 and HDA2 in regulating single var gene expression and the transcriptional regulator ApiAP2-G suggests that both processes share epigenetic regulatory mechanisms and that Plasmodium immune evasion and transmission to new hosts are inextricably linked [27, 29].

Immune evasion is not restricted to blood-stage Plasmodium; when the parasite passes through the mosquito it must also combat a sophisticated innate immune system that is very effective in reducing the parasite burden experienced by the vector. A forward genetic screen and WGS was used to identify the key parasite factor, the surface protein PfS47 (found on the surface of the ookinete as it penetrates the mosquito midgut), that appears to interact with and suppress the vector innate immune system [108]. PfS47 is thought to suppress signaling through the c-Jun N-terminal kinase (JNK) pathway that is critical to an effective immune response [109]. WGS demonstrated that PfS47 has a distinct population structure linked to global distribution. PfS47 is rapidly evolving and selected to achieve JNK suppression in diverse mosquito species, which becomes a key step in the adaptation of Pf to transmission in different vectors, thereby contributing to its broad global distribution [110].

Artemisinin resistance

The goals of MalariaGEN characterize a new approach to understanding parasite population biology. Through the generation and, these days more critically, the management and analysis of the colossal datasets that result from WGS of large numbers of samples, a well-organized study can draw meaningful conclusions. This was applied to perhaps the most serious threat to malaria control that has emerged in recent years—resistance to ART. Using these datasets in meta-analyses with clinical data describing the individual WGS-sequenced samples and outcomes of ART treatment allowed a path to be charted that associated SNPs with treatment features (such as delayed clearance) [111] and identified candidate genes [112]: in both studies a region of chromosome 13 was implicated (Fig. 3d). The precise gene encoding the protein KELCH13 was identified by a combination of “old-fashioned” selection of drug-resistant parasites in the laboratory followed by WGS and comparative genomics of the sensitive parental parasites and the progeny, as well as WGS of ART-resistant field isolates [113, 114]. The role of the kelch13 mutations in ART resistance was proven by direct genome engineering of kelch13 to generate resistant parasites [115, 116]. kelch13 SNPs have been used to map the alarmingly rapid spread of resistance throughout Southeast Asia [116] and it is clear that there is already significant but distinct kelch13 heterogeneity in African Pf strains, although there is no evidence of ART resistance [117–121]. However, in-depth analysis of Southeast Asian ART-resistant parasite genomes [122] revealed that a complex array of background mutations (Fig. 3d) in a variety of genes (encoding ferredoxin (FD), apicoplast ribosomal protein S10 (ARPS10), multidrug resistance protein 2 (MDR2), and Pf chloroquine resistance transporter (CRT)) that are not yet described in African parasites would explain why ART resistance is not (yet) a threat to the use of ART in that continent [121].

A further puzzle was the large number of independent SNPs that seemed capable of mediating ART resistance—typically drug resistance is generated by one or a small number of SNPs focused on either altering the target binding site for the drug or preventing drug access to a binding site buried in the target structure. KELCH proteins are propeller proteins with an iterated structural motif that serves as a platform for the assembly of multi-protein complexes. In addition, KELCH13 has a BTB/POZ domain that might be involved in homodimerization, E3 ubiquitin ligase binding, and transcriptional repression (reviewed in [123]). It has been suggested that ART-resistance-associated kelch13 SNPs might cause a degree of reduced binding of Pf phosphatidylinositol-3-kinase (PI3K), which in turn results in its reduced ubiquitination and consequent degradation of PI3K (Fig. 3d). Elevated levels of PI3K generate increased amounts of its lipid product phosphatidylinositol-3-phosphate (PI3P), which then changes the physiological state of the parasite cell through signaling in as yet unknown pathways [124] but through a mechanism predicated on the proposed abundance of PI3P in the lumen of the endoplasmic reticulum and its proposed role in protein export beyond the parasite vacuole within the host cell [125]. However, aspects of this view have been challenged [126] and further studies are clearly required to resolve the possible role of PI3K signaling in ART resistance. It will be of interest to see if PI3K signaling impacts upon the unfolded protein response implicated in ART resistance using population transcriptomics [127]. The WGS data and two proteomics studies [128, 129] that demonstrate the wide variety of proteins from different cellular compartments of the target parasite that interact with activated ART together suggest that ART resistance is a pleiotropic phenomenon [123]. Therefore, other interrogations, such as metabolomics (see next section), might also be needed to gain functional insights into ART’s mode of action.

Antimalarials

WGS has been instrumental in identifying the cellular target of novel Pf antimalarials as part of the drug discovery pipeline and in following the in vitro selection of resistant parasite lines and validation of observed genomic changes by reverse genetics as described for ART above. This approach has proved highly successful for spiralindolines [130], resulting in the identification of the target of NITD609 (also known as KAE609 or cipargamin) as the P-type ATPase PfATPase4. Furthermore, the translation elongation factor eEF2 has been identified as the target of the 2,6-disubstituted quinoline-4-carboxamide scaffold derivative DDD107498 [131]. WGS is not the only post-genome approach that is useful for determining modes of drug action; metabolomics has a similar potential for analyzing the metabolic changes produced in response to drug exposure and has been utilized in antibiotic [132] and anti-protozoan drug [133] investigations. A metabolomics-based approach also has the advantages that parasite lines resistant to the drug need not be generated and that the activity of pleiotropically acting drugs (such as ART) are directly observed rather than imputed from genomes of resistant parasites.

Vaccines

Post-genomic approaches have also identified promising new Pf vaccine candidates. For example, Pf reticulocyte-binding protein homologue 5 (RH5) binds to the human red cell surface receptor protein basigin, an interaction that is essential for erythrocyte invasion by Pf [134]. Recent WGS studies have shown that both the host and parasite proteins are highly conserved, that antibodies to RH5 block merozoite invasion of erythrocytes [135, 136], and that basigin itself is druggable by recombinant antibodies [137]. Although RH5–basigin interaction offers great promise, the challenges for vaccine development remain considerable and many promising candidates have fallen or will fall by the wayside due to an inability to formulate them to deliver effective vaccination, massive candidate gene sequence variability, and functional non-essentiality of the candidate. WGS will help identify non- or minimally variant candidates and should prove useful in monitoring the effect of vaccination and the analysis of “breakthrough” parasites (those developing in vaccinated individuals), as described in the next section. Effective subunit vaccines will be an invaluable additional approach to vaccination, supplementing other approaches such as the use of the promising but technologically challenging attenuated whole parasite, for example, sporozoite vaccine [138].

Surveillance

The identification of genome signatures of resistance through WGS in the laboratory and increasingly through large-scale genomic epidemiology provides a powerful tool to monitor the emergence of resistance in Plasmodium populations under selective pressure due to administration of both drugs and vaccines. In the case of drugs whose targets have been identified in the laboratory, specific, simple PCR-based assays can be devised. WGS of field parasites under drug pressure is still desirable, however, as alternative resistance mechanisms might emerge that would be missed by targeted assays and, with sufficient depth of sampling, new signatures of resistance would be identifiable from the sequence data. Similar surveillance of parasites that emerge post-vaccination might also be informative. An important analysis of the clinical trial of the RTS,S/AS01 malaria vaccine compared the strain-specific sequence of the gene encoding the circumsporozoite (CS) protein that comprised the vaccine with the CS gene sequences of strains in the infections actually encountered by immunized individuals (between 5 and 17 months of age) [139]. This study demonstrated that homologous protection was greater than protection against heterologous strains and that a cause of failure to protect was simply that the CS protein carried by the infecting parasites did not match that of the vaccine and so perhaps a protective effect was less likely [139]. Therefore, WGS has the power to guide vaccine design based upon the outcomes of trials.

Gene editing

A new era of genetic engineering has dawned with the discovery and development of the bacterial guide RNA template-targeted clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 recombinases as tools for the accurate editing of genomes. The technology has been successfully adapted to many species, including Plasmodium [140], Anopheles [141, 142], and humans (discussed in [143]). Currently, applications of CRISPR-Cas9 to Plasmodium manipulation are restricted to reverse genetic investigations of gene function. However, with the concepts of whole (pre-erythrocytic) parasite vaccines [144, 145], CRISPR-Cas9 offers an obvious route towards the generation of an immunogenic, non-pathogenic parasite that might be suitably safe to administer to humans as a vaccination strategy. Clearly, the engineering of human genomes at any stage of gestation is fraught with ethical considerations [146] and it is inconceivable that this will be applied to the improvement of human resistance to malaria in the foreseeable future. Conversely, although subject to similar ethical and ecological debate, significant conceptual advances towards the generation of CRISPR-Cas9-engineered Anopheles mosquitoes have been quickly achieved. Through harnessing the concept of gene drive, two independent teams have reported the generation of either engineered Anopheles stephensi (a major Indian vector of malaria) that is resistant to malaria [141] or sterile female Ag [142]. Again, owing to ecological considerations, it is unlikely that such engineered mosquitoes, although clearly feasible, will be released into the wild any time soon [147].