Progression of the canonical reference malaria parasite genome from 2002-2019

Annotation and visualization of parasite, fungi and arthropod genomes with Companion

As sequencing genomes has become increasingly popular, the need for annotation of the resulting assemblies is growing. Structural and functional annotation is still challenging as it includes finding the correct gene sequences, annotating other elements such as RNA and being able to submit those data to databases to share it with the community. Compared to de novo assembly where contiguous chromosomes are a sign of high quality, it is difficult to visualize and assess the quality of annotation. We developed the Companion web server to allow non-experts to annotate their genome using a reference-based method, enabling them to assess the output before submitting to public databases. In this update paper, we describe how we have included novel methods for gene finding and made the Companion server more efficient for annotation of genomes of up to 1 Gb in size. The reference set was increased to include genomes of interest for human and animal health from the fungi and arthropod kingdoms. We show that Companion outperforms existing comparable tools where closely related references are available.

Identification of domains in Plasmodium falciparum proteins of unknown function using DALI search on AlphaFold predictions

Plasmodium falciparum, the causative agent of malaria, poses a significant global health challenge, yet much of its biology remains elusive. A third of the genes in the P. falciparum genome lack annotations regarding their function, impeding our understanding of the parasite's biology. In this study, we employ structure predictions and the DALI search algorithm to analyse proteins encoded by uncharacterized genes in the reference strain 3D7 of P. falciparum. By comparing AlphaFold predictions to experimentally determined protein structures in the Protein Data Bank, we found similarities to known domains in 353 proteins of unknown function, shedding light on their potential functions. The lowest-scoring 5% of similarities were additionally validated using the size-independent TM-align algorithm, confirming the detected similarities in 88% of the cases. Notably, in over 70 P. falciparum proteins the presence of domains resembling heptatricopeptide repeats, which are typically involvement in RNA binding and processing, was detected. This suggests this family, which is important in transcription in mitochondria and apicoplasts, is much larger in Plasmodium parasites than previously thought. The results of this domain search provide a resource to the malaria research community that is expected to inform and enable experimental studies.

Natural immunity to malaria preferentially targets the endothelial protein C receptor-binding regions of PfEMP1s

Antibody responses to variant surface antigens (VSAs) produced by the malaria parasite Plasmodium falciparum may contribute to age-related natural immunity to severe malaria. One VSA family, P. falciparum erythrocyte membrane protein-1 (PfEMP1), includes a subset of proteins that binds endothelial protein C receptor (EPCR) in human hosts and potentially disrupts the regulation of inflammatory responses, which may lead to the development of severe malaria. We probed peptide microarrays containing segments spanning five PfEMP1 EPCR-binding domain variants with sera from 10 Malian adults and 10 children to determine the differences between adult and pediatric immune responses. We defined serorecognized peptides and amino acid residues as those that elicited a significantly higher antibody response than malaria-naïve controls. We aimed to identify regions consistently serorecognized among adults but not among children across PfEMP1 variants, potentially indicating regions that drive the development of immunity to severe malaria. Adult sera consistently demonstrated broader and more intense serologic responses to constitutive PfEMP1 peptides than pediatric sera, including peptides in EPCR-binding domains. Both adults and children serorecognized a significantly higher proportion of EPCR-binding peptides than peptides that do not directly participate in receptor binding, indicating a preferential development of serologic responses at functional residues. Over the course of a single malaria transmission season, pediatric serological responses increased between the start and the peak of the season, but waned as the transmission season ended. IMPORTANCE Severe malaria and death related to malaria disproportionately affect sub-Saharan children under 5 years of age, commonly manifesting as cerebral malaria and/or severe malarial anemia. In contrast, adults in malaria-endemic regions tend to experience asymptomatic or mild disease. Our findings indicate that natural immunity to malaria targets specific regions within the EPCR-binding domain, particularly peptides containing EPCR-binding residues. Epitopes containing these residues may be promising targets for vaccines or therapeutics directed against severe malaria. Our approach provides insight into the development of natural immunity to a binding target linked to severe malaria by characterizing an "adult-like" response as recognizing a proportion of epitopes within the PfEMP1 protein, particularly regions that mediate EPCR binding. This "adult-like" response likely requires multiple years of malaria exposure, as increases in pediatric serologic response over a single malaria transmission season do not appear significant.

The global transcriptome of Plasmodium falciparum mid-stage gametocytes (stages II-IV) appears largely conserved and gametocyte-specific gene expression patterns vary in clinical isolates

Our overall understanding of the developmental biology of malaria parasites has been greatly enhanced by recent advances in transcriptomic analysis. However, most of these investigations rely on laboratory strains (LS) that were adapted into in vitro culture many years ago, and the transcriptomes of clinical isolates (CI) circulating in human populations have not been assessed. In this study, RNA-seq was used to compare the global transcriptome of mid-stage gametocytes derived from three short-term cultured CI, with gametocytes derived from the NF54 reference laboratory strain. The core transcriptome appeared to be consistent between CI- and LS-derived gametocyte preparations, but some important differences were also observed. A majority of gametocyte-specific genes (43/53) appear to have relatively higher expression in CI-derived gametocytes than in LS-derived gametocytes, but a K-means clustering analysis showed that genes involved in flagellum- and microtubule-based processes (movement/motility) were more abundant in both groups, albeit with some differences between them. In addition, gametocytes from one CI described as CI group II gametocytes (CI:GGII) showed gene expression variation in the form of reduced gametocyte-specific gene expression compared to the other two CI-derived gametocytes (CI gametocyte group I, CI:GGI), although the mixed developmental stages used in our study is a potential confounder, only partially mitigated by the inclusion of multiple replicates for each CI. Overall, our study suggests that there may be subtle differences in the gene expression profiles of mid-stage gametocytes from CI relative to the NF54 reference strain of Plasmodium falciparum. Thus, it is necessary to deploy gametocyte-producing clinical parasite isolates to fully understand the diversity of gene expression strategies that may occur during the sequestered development of parasite sexual stages. IMPORTANCE Maturing gametocytes of Plasmodium falciparum are known to sequester away from peripheral circulation into the bone marrow until they are mature. Blocking gametocyte sequestration can prevent malaria transmission from humans to mosquitoes, but most studies aim to understand gametocyte development utilizing long-term adapted laboratory lines instead of clinical isolates. This is a particular issue for our understanding of the sexual stages, which are known to decrease rapidly during adaptation to long-term culture, meaning that many LS are unable to produce transmissible gametocytes. Using RNA-seq, we investigated the global transcriptome of mid-stage gametocytes derived from three clinical isolates and a reference strain (NF54). This identified important differences in gene expression profiles between immature gametocytes of CI and the NF54 reference strain of P. falciparum, suggesting increased investment in gametocytogenesis in clinical isolates. Our transcriptomic data highlight the use of clinical isolates in studying the morphological, cellular features and molecular biology of gametocytes.

A novel Modulator of Ring Stage Translation (MRST) gene alters artemisinin sensitivity in Plasmodium falciparum

The implementation of artemisinin (ART) combination therapies (ACTs) has greatly decreased deaths caused by Plasmodium falciparum malaria, but increasing ACT resistance in Southeast Asia and Africa could reverse this progress. Parasite population genetic studies have identified numerous genes, single-nucleotide polymorphisms (SNPs), and transcriptional signatures associated with altered artemisinin activity with SNPs in the Kelch13 (K13) gene being the most well-characterized artemisinin resistance marker. However, there is an increasing evidence that resistance to artemisinin in P. falciparum is not related only to K13 SNPs, prompting the need to characterize other novel genes that can alter ART responses in P. falciparum. In our previous analyses of P. falciparum piggyBac mutants, several genes of unknown function exhibited increased sensitivity to artemisinin that was similar to a mutant of K13. Further analysis of these genes and their gene co-expression networks indicated that the ART sensitivity cluster was functionally linked to DNA replication and repair, stress responses, and maintenance of homeostatic nuclear activity. In this study, we have characterized PF3D7_1136600, another member of the ART sensitivity cluster. Previously annotated as a conserved Plasmodium gene of unknown function, we now provide putative annotation of this gene as a Modulator of Ring Stage Translation (MRST). Our findings reveal that the mutagenesis of MRST affects gene expression of multiple translation-associated pathways during the early ring stage of asexual development via putative ribosome assembly and maturation activity, suggesting an essential role of MRST in protein biosynthesis and another novel mechanism of altering the parasite's ART drug response.IMPORTANCEPlasmodium falciparum malaria killed more than 600,000 people in 2021, though ACTs have been critical in reducing malaria mortality as a first-line treatment for infection. However, ACT resistance in Southeast Asia and emerging resistance in Africa are detrimental to this progress. Mutations to Kelch13 (K13) have been identified to confer increased artemisinin tolerance in field isolates, however, genes other than K13 are implicated in altering how the parasite responds to artemisinin prompts additional analysis. Therefore, in this study we have characterized a P. falciparum mutant clone with altered sensitivity to artemisinin and identified a novel gene (PF3D7_1136600) that is associated with alterations to parasite translational metabolism during critical timepoints for artemisinin drug response. Many genes of the P. falciparum genome remain unannotated, posing a challenge for drug-gene characterizations in the parasite. Therefore, through this study, we have putatively annotated PF3D7_1136600 as a novel MRST gene and have identified a potential link between MRST and parasite stress response mechanisms.

Why is Babesia not killed by artemisinin like Plasmodium?

Babesia spp. are intraerythrocytic apicomplexans that digest and utilize red blood cells in a similar way to intraerythrocytic Plasmodium spp., but unlike the latter, are not sensitive to artemisinin. A comparison of Babesia and Plasmodium genomes revealed that Babesia genomes, which are smaller than those of Plasmodium, lack numerous genes, and especially haem synthesis-related genes, that are found in the latter. Single-cell sequencing analysis showed that the different treatment groups of Babesia microti with expressed pentose phosphate pathway-related, DNA replication-related, antioxidation-related, glycolysis-related, and glutathione-related genes were not as sensitive to artemether as Plasmodium yoelii 17XNL. In particular, pentose phosphate pathway-related, DNA replication-related, and glutathione-related genes, which were actively expressed in P. yoelii 17XNL, were not actively expressed in B. microti. Supplying iron in vivo can promote the reproduction of B. microti. These results suggest that Babesia spp. lack a similar mechanism to that of malaria parasites through which the haem or iron in hemoglobin is utilized, and that this likely leads to their insensitivity to artemisinin.

Genomic variation during culture adaptation of genetically complex Plasmodium falciparum clinical isolates

Experimental studies on the biology of malaria parasites have mostly been based on laboratory-adapted lines, but there is limited understanding of how these may differ from parasites in natural infections. Loss-of-function mutants have previously been shown to emerge during culture of some Plasmodium falciparum clinical isolates in analyses focusing on single-genotype infections. The present study included a broader array of isolates, mostly representing multiple-genotype infections, which are more typical in areas where malaria is highly endemic. Genome sequence data from multiple time points over several months of culture adaptation of 28 West African isolates were analysed, including previously available sequences along with new genome sequences from additional isolates and time points. Some genetically complex isolates eventually became fixed over time to single surviving genotypes in culture, whereas others retained diversity, although proportions of genotypes varied over time. Drug resistance allele frequencies did not show overall directional changes, suggesting that resistance-associated costs are not the main causes of fitness differences among parasites in culture. Loss-of-function mutants emerged during culture in several of the multiple-genotype isolates, affecting genes (including AP2-HS, EPAC and SRPK1) for which loss-of-function mutants were previously seen to emerge in single-genotype isolates. Parasite clones were derived by limiting dilution from six of the isolates, and sequencing identified de novo variants not detected in the bulk isolate sequences. Interestingly, several of these were nonsense mutants and frameshifts disrupting the coding sequence of EPAC, the gene with the largest number of independent nonsense mutants previously identified in laboratory-adapted lines. Analysis of genomic identity by descent to explore relatedness among clones revealed co-occurring non-identical sibling parasites, illustrative of the natural genetic structure within endemic populations.

Genome-wide liver transcriptomic profiling of a malaria mouse model reveals disturbed immune and metabolic responses

BACKGROUND

The liver is responsible for a range of functions in vertebrates, such as metabolism and immunity. In malaria, the liver plays a crucial role in the interaction between the parasite and host. Although malarial hepatitis is a common clinical complication of severe malaria, other malaria-related liver changes have been overlooked during the blood stage of the parasite life-cycle, in contrast to the many studies that have focused on parasite invasion of and replication in the liver during the hepatic stage of the parasite.

METHODS

A rodent model of malaria was established using Plasmodium yoelii strain 17XL, a lethal strain of rodent malaria, for liver transcriptomic profiling.

RESULTS

Differentially expressed messenger RNAs were associated with innate and adaptive immune responses, while differentially expressed long noncoding RNAs were enriched in the regulation of metabolism-related pathways, such as lipid metabolism. The coexpression network showed that host genes were related to cellular transport and tissue remodeling. Hub gene analysis of P. yoelii indicated that ubiquitination genes that were coexpressed with the host were evolutionarily conserved.

CONCLUSIONS

Our analysis yielded evidence of activated immune responses, aberrant metabolic processes and tissue remodeling changes in the livers of mice with malaria during the blood stage of the parasite, which provided a systematic outline of liver responses during Plasmodium infection.

Protein KIC5 is a novel regulator of artemisinin stress response in the malaria parasite Plasmodium falciparum

Artemisinin combination therapies (ACTs) have led to a significant decrease in Plasmodium falciparum malaria mortality. This progress is now threatened by emerging artemisinin resistance (ART-R) linked originally in SE Asia to polymorphisms in the Kelch propeller protein (K13) and more recently to several other seemingly unrelated genetic mutations. To better understand the parasite response to ART, we are characterizing a P. falciparum mutant with altered sensitivity to ART that was created via piggyBac transposon mutagenesis. The transposon inserted near the putative transcription start site of a gene defined as a "Plasmodium-conserved gene of unknown function," now functionally linked to K13 as the Kelch13 Interacting Candidate 5 protein (KIC5). Phenotype analysis of the KIC5 mutant during intraerythrocytic asexual development identified transcriptional changes associated with DNA stress response and altered mitochondrial metabolism, linking dysregulation of the KIC5 gene to the parasite's ability to respond to ART exposure. Through characterization of the KIC5 transcriptome, we hypothesize that this gene may be essential under ART exposure to manage gene expression of the wild-type stress response at early ring stage, thereby providing a better understanding of the parasite's processes that can alter ART sensitivity.

Plasmodium falciparum Sexual Commitment Rate Variation among Clinical Isolates and Diverse Laboratory-Adapted Lines

Asexual blood-stage malaria parasites must produce sexual progeny to infect mosquitoes. It is important to understand the scope and causes of intraspecific variation in sexual commitment rates, particularly for the major human parasite P. falciparum. First, two alternative assay methods of measuring sexual commitment were compared to test a genetically modified P. falciparum line with elevated commitment rates inducible by overexpression of GDV1. The methods yielded correlated measurements with higher sensitivity and precision being achieved by one employing detection of the early gametocyte differentiation marker Pfs16. Thus, this was used to survey a diverse range of parasite lines and test each in multiple biological replicate assays in a serum-free medium supplemented with Albumax. There were differences among six recent clinical isolates from Ghana in their mean rates of sexual commitment per cycle, ranging from 3.3% to 12.2%. Among 13 diverse long-term laboratory-adapted lines, mean sexual commitment rates for most ranged from 4.7% to 13.4%, a few had lower rates with means from 0.3 to 1.6%, and one with a nonfunctional ap2-g gene always showed zero commitment. Among a subset of lines tested for the effects of exogenous choline to suppress commitment, there were significant differences. As expected, there was no effect in a line that had lost the gdv1 gene and that had generally low commitment, whereas the others showed quantitatively variable but significant responses to choline, suggesting potential trait variation. The results indicated the value of performing multiple replicate assays for understanding the variation of this key reproductive trait that likely affects transmission. IMPORTANCE Only sexual-stage malaria parasites are transmitted from human blood to mosquitoes. Thus, it is vital to understand variations in sexual commitment rates because these may be modifiable or susceptible to blocking. Two different methods of commitment rate measurement were first compared, demonstrating higher sensitivity and precision by the detection of an early differentiation marker, which was subsequently used to survey diverse lines. Clinical isolates from Ghana showed significant variation in mean per-cycle commitment rates and variation among biological replicates. Laboratory-adapted lines of diverse origins had a wider range with most being within the range observed for the clinical isolates, while a minority consistently had lower or zero rates. There was quantitative variation in the effects when adding choline to suppress commitment, indicating differing responsiveness of parasites to this environmental modification. Performing multiple assay replicates and comparisons of diverse isolates was important to understand this trait and its potential effects on transmission.

A manually curated annotation characterises genomic features of P. falciparum lncRNAs

BACKGROUND

Important regulation occurs at the level of transcription in Plasmodium falciparum and growing evidence suggests that these apicomplexan parasites have complex regulatory networks. Recent studies implicate long noncoding RNAs (lncRNAs) as transcriptional regulators in P. falciparum. However, due to limited research and the lack of necessary experimental tools, our understanding of their role in the malaria-causing parasite remains largely unelucidated. In this work, we address one of these limitations, the lack of an updated and improved lncRNA annotation in P. falciparum.

RESULTS

We generated long-read RNA sequencing data and integrated information extracted and curated from multiple sources to manually annotate lncRNAs. We identified 1119 novel lncRNAs and validated and refined 1250 existing annotations. Utilising the collated datasets, we generated evidence-based ranking scores for each annotation and characterised the distinct genomic contexts and features of P. falciparum lncRNAs. Certain features indicated subsets with potential biological significance such as 25 lncRNAs containing multiple introns, 335 lncRNAs lacking mutations in piggyBac mutagenic studies and lncRNAs associated with specific biologic processes including two new types of lncRNAs found proximal to var genes.

CONCLUSIONS

The insights and the annotation presented in this study will serve as valuable tools for researchers seeking to understand the role of lncRNAs in parasite biology through both bioinformatics and experimental approaches.

Cation ATPase (ATP4) Orthologue Replacement in the Malaria Parasite Plasmodium knowlesi Reveals Species-Specific Responses to ATP4-Targeting Drugs

Several unrelated classes of antimalarial compounds developed against Plasmodium falciparum target a parasite-specific P-type ATP-dependent Na+ pump, PfATP4. We have previously shown that other malaria parasite species infecting humans are less susceptible to these compounds. Here, we generated a series of transgenic Plasmodium knowlesi orthologue replacement (OR) lines in which the endogenous pkatp4 locus was replaced by a recodonized P. knowlesi atp4 (pkatp4) coding region or the orthologous coding region from P. falciparum, Plasmodium malariae, Plasmodium ovale subsp. curtisi, or Plasmodium vivax. Each OR transgenic line displayed a similar growth pattern to the parental P. knowlesi line. We found significant orthologue-specific differences in parasite susceptibility to three chemically unrelated ATP4 inhibitors, but not to comparator drugs, among the P. knowlesi OR lines. The PfATP4OR transgenic line of P. knowlesi was significantly more susceptible than our control PkATP4OR line to three ATP4 inhibitors: cipargamin, PA21A092, and SJ733. The PvATP4OR and PmATP4OR lines were similarly susceptible to the control PkATP4OR line, but the PocATP4OR line was significantly less susceptible to all ATP4 inhibitors than the PkATP4OR line. Cipargamin-induced inhibition of Na+ efflux was also significantly greater with the P. falciparum orthologue of ATP4. This confirms that species-specific susceptibility differences previously observed in ex vivo studies of human isolates are partly or wholly enshrined in the primary amino acid sequences of the respective ATP4 orthologues and highlights the need to monitor efficacy of investigational malaria drugs against multiple species. P. knowlesi is now established as an important in vitro model for studying drug susceptibility in non-falciparum malaria parasites. IMPORTANCE Effective drugs are vital to minimize the illness and death caused by malaria. Development of new drugs becomes ever more urgent as drug resistance emerges. Among promising compounds now being developed to treat malaria are several unrelated molecules that each inhibit the same protein in the malaria parasite-ATP4. Here, we exploited the genetic tractability of P. knowlesi to replace its own ATP4 genes with orthologues from five human-infective species to understand the drug susceptibility differences among these parasites. We previously estimated the susceptibility to ATP4-targeting drugs of each species using clinical samples from malaria patients. These estimates closely matched those of the corresponding "hybrid" P. knowlesi parasites carrying introduced ATP4 genes. Thus, species-specific ATP4 inhibitor efficacy is directly determined by the sequence of the gene. Our novel approach to understanding cross-species susceptibility/resistance can strongly support the effort to develop antimalarials that effectively target all human malaria parasite species.

Highly Variable Expression of Merozoite Surface Protein MSPDBL2 in Diverse Plasmodium falciparum Clinical Isolates and Transcriptome Scans for Correlating Genes

The merozoite surface protein MSPDBL2 of Plasmodium falciparum is under strong balancing selection and is a target of naturally acquired antibodies. Remarkably, MSPDBL2 is expressed in only a minority of mature schizonts of any cultured parasite line, and mspdbl2 gene transcription increases in response to overexpression of the gametocyte development inducer GDV1, so it is important to understand its natural expression. Here, MSPDBL2 in mature schizonts was analyzed in the first ex vivo culture cycle of 96 clinical isolates from 4 populations with various levels of infection endemicity in different West African countries, by immunofluorescence microscopy with antibodies against a conserved region of the protein. In most isolates, less than 1% of mature schizonts were positive for MSPDBL2, but the frequency distribution was highly skewed, as nine isolates had more than 3% schizonts positive and one had 73% positive. To investigate whether the expression of other gene loci correlated with MSPDBL2 expression, whole-transcriptome sequencing was performed on schizont-enriched material from 17 of the isolates with a wide range of proportions of schizonts positive. Transcripts of particular genes were highly significantly positively correlated with MSPDBL2 positivity in schizonts as well as with mspdbl2 gene transcript levels, showing overrepresentation of genes implicated previously as involved in gametocytogenesis but not including the gametocytogenesis master regulator ap2-g. Single-cell transcriptome analysis of a laboratory-adapted clone showed that most individual parasites expressing mspdbl2 did not express ap2-g, consistent with MSPDBL2 marking a developmental subpopulation that is distinct but likely to co-occur alongside sexual commitment. IMPORTANCE These findings contribute to understanding malaria parasite antigenic and developmental variation, focusing on the merozoite surface protein encoded by the single locus under strongest balancing selection. Analyzing the initial ex vivo generation of parasites grown from a wide sample of clinical infections, we show a unique and highly skewed pattern of natural expression frequencies of MSPDBL2, distinct from that of any other antigen. Bulk transcriptome analysis of a range of clinical isolates showed significant overrepresentation of sexual development genes among those positively correlated with MSPDBL2 protein and mspdbl2 gene expression, indicating the MSPDBL2-positive subpopulation to be often coincident with parasites developing sexually in preparation for transmission. Single-cell transcriptome data confirm the absence of a direct correlation with the ap2-g master regulator of sexual development, indicating that the MSPDBL2-positive subpopulation has a separate function in asexual survival and replication under conditions that promote terminal sexual differentiation.

Characterising genome architectures using genome decomposition analysis

Genome architecture describes how genes and other features are arranged in genomes. These arrangements reflect the evolutionary pressures on genomes and underlie biological processes such as chromosomal segregation and the regulation of gene expression. We present a new tool called Genome Decomposition Analysis (GDA) that characterises genome architectures and acts as an accessible approach for discovering hidden features of a genome assembly. With the imminent deluge of high-quality genome assemblies from projects such as the Darwin Tree of Life and the Earth BioGenome Project, GDA has been designed to facilitate their exploration and the discovery of novel genome biology. We highlight the effectiveness of our approach in characterising the genome architectures of single-celled eukaryotic parasites from the phylum Apicomplexa and show that it scales well to large genomes.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Identifying transcript 5' capped ends in Plasmodium falciparum

Background

The genome of the human malaria parasite Plasmodium falciparum is poorly annotated, in particular, the 5' capped ends of its mRNA transcripts. New approaches are needed to fully catalog P. falciparum transcripts for understanding gene function and regulation in this organism.

Methods

We developed a transcriptomic method based on next-generation sequencing of complementary DNA (cDNA) enriched for full-length fragments using eIF4E, a 5' cap-binding protein, and an unenriched control. DNA sequencing adapter was added after enrichment of full-length cDNA using two different ligation protocols. From the mapped sequence reads, enrichment scores were calculated for all transcribed nucleotides and used to calculate P-values of 5' capped nucleotide enrichment. Sensitivity and accuracy were increased by combining P-values from replicate experiments. Data were obtained for P. falciparum ring, trophozoite and schizont stages of intra-erythrocytic development.

Results

5' capped nucleotide signals were mapped to 17,961 non-overlapping P. falciparum genomic intervals. Analysis of the dominant 5' capped nucleotide in these genomic intervals revealed the presence of two groups with distinctive epigenetic features and sequence patterns. A total of 4,512 transcripts were annotated as 5' capped based on the correspondence of 5' end with 5' capped nucleotide annotated from full-length cDNA data.

Discussion

The presence of two groups of 5' capped nucleotides suggests that alternative mechanisms may exist for producing 5' capped transcript ends in P. falciparum. The 5' capped transcripts that are antisense, outside of, or partially overlapping coding regions may be important regulators of gene function in P. falciparum.

Chemoproteomics for Plasmodium Parasite Drug Target Discovery

Emerging Plasmodium parasite drug resistance is threatening progress towards malaria control and elimination. While recent efforts in cell-based, high-throughput drug screening have produced first-in-class drugs with promising activities against different Plasmodium life cycle stages, most of these antimalarial agents have elusive mechanisms of action. Though challenging to address, target identification can provide valuable information to facilitate lead optimization and preclinical drug prioritization. Recently, proteome-wide methods for direct assessment of drug-protein interactions have emerged as powerful tools in a number of systems, including Plasmodium. In this review, we will discuss current chemoproteomic strategies that have been adapted to antimalarial drug target discovery, including affinity- and activity-based protein profiling and the energetics-based techniques thermal proteome profiling and stability of proteins from rates of oxidation. The successful application of chemoproteomics to the Plasmodium blood stage highlights the potential of these methods to link inhibitors to their molecular targets in more elusive Plasmodium life stages and intracellular pathogens in the future.

Screening and characterization of hypothetical proteins of Plasmodium falciparum as novel vaccine candidates in the fight against malaria using reverse vaccinology

BACKGROUND

Plasmodium falciparum is the most deadly and leading cause of morbidity and mortality in Africa. About 90% of all malaria deaths in the world today occur in Sub-Saharan Africa especially in children aged < 5 years. In 2018, it was reported that there were 228 million malaria cases that resulted in 405,000 deaths from 91 countries. Currently, a fully effective and long-lasting preventive malaria vaccine is still elusive therefore more effort is needed to identify better effective vaccine candidates. The aim of this study was to identify and characterize hypothetical proteins as vaccine candidates derived from Plasmodium falciparum 3D7 genome by reverse vaccinology.

RESULTS

Of the 23 selected hypothetical proteins, 5 were predicted on the extracellular localization by WoLFPSORTv.2.0 program and all the 5 had less than 2 transmembrane regions that were predicted by TMHMMv2.0 and HMMTOP programs at default settings. Two out of the five proteins lacked secretory signal peptides as predicted by SignalP program. Among the 5 extracellular proteins, 3 were predicted to be antigenic by VaxiJen (score ≥ 0.5) and had negative GRAVY values ranging from - 1.156 to - 0.440. B cell epitope prediction by ABCpred and BCpred programs revealed a total of 15 antigenic epitopes. A total of 13 cytotoxic T cells were predicted from the 3 proteins using CTLPred online server. Only 2 out of the 13 CTL were antigenic, immunogenic, non-allergenic, and non-toxic using VaxiJen, IEDB, AllergenFp, and Toxinpred servers respectively in that order. Five HTL peptides from XP_001351030.1 protein are predicted inducers of all the three cytokines. STRING protein-protein network analysis of HPs revealed XP_001350955.1 closely interacts with nucleoside diphosphate kinase (PF13-0349) at 0.704, XP_001351030.1 interacts with male development protein1 (Mdv-1) at 0.645, and XP_001351047.1 with an uncharacterized protein (MAL8P1.53) at 0.400.

CONCLUSION

Reverse vaccinology is a promising strategy for the screening and identification of antigenic antigens with potential capacity to elicit cellular and humoral immune responses against P. falciparum infection. In this study, potential vaccine candidates of Plasmodium falciparum were identified and screened using standard bioinformatics tools. The vaccine candidates contained antigenic and immunogenic epitopes which could be considered for novel and effective vaccine targets. However, we strongly recommend in vivo and in vitro experiments to validate their immunogenicity and protective efficacy to completely decipher the vaccine targets against malaria.

Peculiarities of Plasmodium falciparum Gene Regulation and Chromatin Structure

The highly complex life cycle of the human malaria parasite, Plasmodium falciparum, is based on an orchestrated and tightly regulated gene expression program. In general, eukaryotic transcription regulation is determined by a combination of sequence-specific transcription factors binding to regulatory DNA elements and the packaging of DNA into chromatin as an additional layer. The accessibility of regulatory DNA elements is controlled by the nucleosome occupancy and changes of their positions by an active process called nucleosome remodeling. These epigenetic mechanisms are poorly explored in P. falciparum. The parasite genome is characterized by an extraordinarily high AT-content and the distinct architecture of functional elements, and chromatin-related proteins also exhibit high sequence divergence compared to other eukaryotes. Together with the distinct biochemical properties of nucleosomes, these features suggest substantial differences in chromatin-dependent regulation. Here, we highlight the peculiarities of epigenetic mechanisms in P. falciparum, addressing chromatin structure and dynamics with respect to their impact on transcriptional control. We focus on the specialized chromatin remodeling enzymes and discuss their essential function in P. falciparum gene regulation.

A PCR-Based Technique to Track the Geographic Origin of Plasmodium falciparum With 23-SNP Barcode Analysis

Increased population movement has increased the risk of reintroducing parasites to elimination areas and also dispersing drug-resistant parasites to new regions. Therefore, reliable and repeatable methods to trace back to the source of imported infections are essential. The recently developed 23-single-nucleotide polymorphism (SNP) barcode from organellar genomes of mitochondrion (mt) and apicoplast (apico) provides a valuable tool to locate the geographic origin of Plasmodium falciparum. This study aims to explore the feasibility of using the 23-SNP barcode for tracking P. falciparum by polymerase chain reaction and sequencing, while providing geographical haplotypes of isolates that originated from Central Africa. Based on 23-SNP barcode analysis, SNPs were found at seven loci; 27 isolates were confirmed to have originated in West Africa, and this study also showed four isolates from Central Africa (Equatorial Guinea, 3; Republic of Congo, 1) that originated in East Africa. This study provides the sequence data from Central Africa and fills 23-SNP barcode data gaps of sample origins.

Targeting Gametocytes of the Malaria Parasite Plasmodium falciparum in a Functional Genomics Era: Next Steps

Mosquito transmission of the deadly malaria parasite Plasmodium falciparum is mediated by mature sexual forms (gametocytes). Circulating in the vertebrate host, relatively few intraerythrocytic gametocytes are picked up during a bloodmeal to continue sexual development in the mosquito vector. Human-to-vector transmission thus represents an infection bottleneck in the parasite's life cycle for therapeutic interventions to prevent malaria. Even though recent progress has been made in the identification of genetic factors linked to gametocytogenesis, a plethora of genes essential for sexual-stage development are yet to be unraveled. In this review, we revisit P. falciparum transmission biology by discussing targetable features of gametocytes and provide a perspective on a forward-genetic approach for identification of novel transmission-blocking candidates in the future.

Mapping immune variation and var gene switching in naive hosts infected with Plasmodium falciparum

Falciparum malaria is clinically heterogeneous and the relative contribution of parasite and host in shaping disease severity remains unclear. We explored the interaction between inflammation and parasite variant surface antigen (VSA) expression, asking whether this relationship underpins the variation observed in controlled human malaria infection (CHMI). We uncovered marked heterogeneity in the host response to blood challenge; some volunteers remained quiescent, others triggered interferon-stimulated inflammation and some showed transcriptional evidence of myeloid cell suppression. Significantly, only inflammatory volunteers experienced hallmark symptoms of malaria. When we tracked temporal changes in parasite VSA expression to ask whether variants associated with severe disease rapidly expand in naive hosts, we found no transcriptional evidence to support this hypothesis. These data indicate that parasite variants that dominate severe malaria do not have an intrinsic growth or survival advantage; instead, they presumably rely upon infection-induced changes in their within-host environment for selection.

In silico identification of novel open reading frames in Plasmodium falciparum oocyte and salivary gland sporozoites using proteogenomics framework

Background

Plasmodium falciparum causes the deadliest form of malaria, which remains one of the most prevalent infectious diseases. Unfortunately, the only licensed vaccine showed limited protection and resistance to anti-malarial drug is increasing, which can be largely attributed to the biological complexity of the parasite’s life cycle. The progression from one developmental stage to another in P. falciparum involves drastic changes in gene expressions, where its infectivity to human hosts varies greatly depending on the stage. Approaches to identify candidate genes that are responsible for the development of infectivity to human hosts typically involve differential gene expression analysis between stages. However, the detection may be limited to annotated proteins and open reading frames (ORFs) predicted using restrictive criteria.

Methods

The above problem is particularly relevant for P. falciparum; whose genome annotation is relatively incomplete given its clinical significance. In this work, systems proteogenomics approach was used to address this challenge, as it allows computational detection of unannotated, novel Open Reading Frames (nORFs), which are neglected by conventional analyses. Two pairs of transcriptome/proteome were obtained from a previous study where one was collected in the mosquito-infectious oocyst sporozoite stage, and the other in the salivary gland sporozoite stage with human infectivity. They were then re-analysed using the proteogenomics framework to identify nORFs in each stage.

Results

Translational products of nORFs that map to antisense, intergenic, intronic, 3′ UTR and 5′ UTR regions, as well as alternative reading frames of canonical proteins were detected. Some of these nORFs also showed differential expression between the two life cycle stages studied. Their regulatory roles were explored through further bioinformatics analyses including the expression regulation on the parent reference genes, in silico structure prediction, and gene ontology term enrichment analysis.

Conclusion

The identification of nORFs in P. falciparum sporozoites highlights the biological complexity of the parasite. Although the analyses are solely computational, these results provide a starting point for further experimental validation of the existence and functional roles of these nORFs,

Essential Genes of the Parasitic Apicomplexa

Genome-scale mutagenesis screens for genes essential for apicomplexan parasite survival have been completed in three species: Plasmodium falciparum, the major human malaria parasite, Plasmodium berghei, a model rodent malaria parasite, and the more distantly related Toxoplasma gondii, the causative agent of toxoplasmosis. These three species share 2606 single-copy orthologs, 1500 of which have essentiality data in all three screens. In this review, we explore the overlap between these datasets to define the core essential genes of the phylum Apicomplexa. We further discuss the implications of these groundbreaking studies for understanding apicomplexan parasite biology, and we identify promising areas of focus for developing new pan-apicomplexan parasite interventions.

Structure of the Plasmodium-interspersed repeat proteins of the malaria parasite

The Plasmodium parasites that cause malaria replicate within blood cells of an infected host. These parasites send a small number of proteins to infected blood cell surfaces, allowing them to bind host molecules but also risking their detection by the host immune system. These proteins have diversified into large families, allowing the parasite to avoid detection by using antigenic variation. The most ubiquitous of these families is the Plasmodium-interspersed repeat (PIR) protein family. Here we present the structure of a PIR protein, revealing the architecture of its ectodomain and showing how it has diversified. Finally, we use structure-guided methods to understand which small variant surface antigen families are PIRs and to understand their evolution across malaria parasites.

The deadly symptoms of malaria occur as Plasmodium parasites replicate within blood cells. Members of several variant surface protein families are expressed on infected blood cell surfaces. Of these, the largest and most ubiquitous are the Plasmodium-interspersed repeat (PIR) proteins, with more than 1,000 variants in some genomes. Their functions are mysterious, but differential pir gene expression associates with acute or chronic infection in a mouse malaria model. The membership of the PIR superfamily, and whether the family includes Plasmodium falciparum variant surface proteins, such as RIFINs and STEVORs, is controversial. Here we reveal the structure of the extracellular domain of a PIR from Plasmodium chabaudi. We use structure-guided sequence analysis and molecular modeling to show that this fold is found across PIR proteins from mouse- and human-infective malaria parasites. Moreover, we show that RIFINs and STEVORs are not PIRs. This study provides a structure-guided definition of the PIRs and a molecular framework to understand their evolution.

CaaX-Like Protease of Cyanobacterial Origin Is Required for Complex Plastid Biogenesis in Malaria Parasites

Plasmodium parasites, which cause malaria, and related apicomplexans are important human and veterinary pathogens. These parasites represent a highly divergent and understudied branch of eukaryotes, and as such often defy the expectations set by model organisms. One striking example of unique apicomplexan biology is the apicoplast, an essential but nonphotosynthetic plastid derived from an unusual secondary (eukaryote-eukaryote) endosymbiosis. Endosymbioses are a major driver of cellular innovation, and apicoplast biogenesis pathways represent a hot spot for molecular evolution. We previously conducted an unbiased screen for apicoplast biogenesis genes in P. falciparum to uncover these essential and innovative pathways. Here, we validate a novel gene candidate from our screen and show that its role in apicoplast biogenesis does not match its functional annotation predicted by model eukaryotes. Our findings suggest that an uncharacterized chloroplast maintenance pathway has been reused for complex plastid biogenesis in this divergent branch of pathogens.

Plasmodium parasites and related apicomplexans contain an essential “complex plastid” organelle of secondary endosymbiotic origin, the apicoplast. Biogenesis of this complex plastid poses a unique challenge requiring evolution of new cellular machinery. We previously conducted a mutagenesis screen for essential apicoplast biogenesis genes to discover organellar pathways with evolutionary and biomedical significance. Here we validate and characterize a gene candidate from our screen, Pf3D7_0913500. Using a conditional knockdown strain, we show that Pf3D7_0913500 depletion causes growth inhibition that is rescued by the sole essential product of the apicoplast, isopentenyl pyrophosphate (IPP), and results in apicoplast loss. Because Pf3D7_0913500 had no previous functional annotation, we name it apicoplast-minus IPP-rescued 4 (AMR4). AMR4 has an annotated CaaX protease and bacteriocin processing (CPBP) domain, which in eukaryotes typically indicates a role in CaaX postprenylation processing. Indeed, AMR4 is the only putative CaaX-like protease in Plasmodium parasites which are known to require protein prenylation, and we confirm that the conserved catalytic residue of AMR4 (E352) is required for its apicoplast function. However, we unexpectedly find that AMR4 does not act in a CaaX postprenylation processing pathway in Plasmodium falciparum. Instead, we find that AMR4 is imported into the apicoplast and is derived from a cyanobacterial CPBP gene which was retained through both primary and secondary endosymbiosis. Our findings suggest that AMR4 is not a true CaaX protease, but instead it performs a conserved, uncharacterized chloroplast function that has been retained for complex plastid biogenesis.

Redesigned TetR-Aptamer System To Control Gene Expression in Plasmodium falciparum

Malaria elimination efforts have been repeatedly hindered by the evolution and spread of multidrug-resistant strains of Plasmodium falciparum. The absence of a commercially available vaccine emphasizes the need for a better understanding of Plasmodium biology in order to further translational research. This has been partly facilitated by targeted gene deletion strategies for the functional analysis of parasite genes. However, genes that are essential for parasite replication in erythrocytes are refractory to such methods, and require conditional knockdown or knockout approaches to dissect their function. One such approach is the TetR-DOZI system that employs multiple synthetic aptamers in the untranslated regions of target genes to control their expression in a tetracycline-dependent manner. Maintaining modified parasites with intact aptamer copies has been challenging since these repeats can be lost by recombination. By interspacing the aptamers with unique sequences, we created a stable genetic system that remains effective at controlling target gene expression.

One of the most powerful approaches to understanding gene function involves turning genes on and off at will and measuring the impact at the cellular or organismal level. This particularly applies to the cohort of essential genes where traditional gene knockouts are inviable. In Plasmodium falciparum, conditional control of gene expression has been achieved by using multicomponent systems in which individual modules interact with each other to regulate DNA recombination, transcription, or posttranscriptional processes. The recently devised TetR-DOZI aptamer system relies on the ligand-regulatable interaction of a protein module with synthetic RNA aptamers to control the translation of a target gene. This technique has been successfully employed to study essential genes in P. falciparum and involves the insertion of several aptamer copies into the 3′ untranslated regions (UTRs), which provide control over mRNA fate. However, aptamer repeats are prone to recombination and one or more copies can be lost from the system, resulting in a loss of control over target gene expression. We rectified this issue by redesigning the aptamer array to minimize recombination while preserving the control elements. As proof of concept, we compared the original and modified arrays for their ability to knock down the levels of a putative essential apicoplast protein (PF3D7_0815700) and demonstrated that the modified array is highly stable and efficient. This redesign will enhance the utility of a tool that is quickly becoming a favored strategy for genetic studies in P. falciparum.

IMPORTANCE Malaria elimination efforts have been repeatedly hindered by the evolution and spread of multidrug-resistant strains of Plasmodium falciparum. The absence of a commercially available vaccine emphasizes the need for a better understanding of Plasmodium biology in order to further translational research. This has been partly facilitated by targeted gene deletion strategies for the functional analysis of parasite genes. However, genes that are essential for parasite replication in erythrocytes are refractory to such methods, and require conditional knockdown or knockout approaches to dissect their function. One such approach is the TetR-DOZI system that employs multiple synthetic aptamers in the untranslated regions of target genes to control their expression in a tetracycline-dependent manner. Maintaining modified parasites with intact aptamer copies has been challenging since these repeats can be lost by recombination. By interspacing the aptamers with unique sequences, we created a stable genetic system that remains effective at controlling target gene expression.

Genomic and transcriptomic evidence for descent from Plasmodium and loss of blood schizogony in Hepatocystis parasites from naturally infected red colobus monkeys

Hepatocystis is a genus of single-celled parasites infecting, amongst other hosts, monkeys, bats and squirrels. Although thought to have descended from malaria parasites (Plasmodium spp.), Hepatocystis spp. are thought not to undergo replication in the blood–the part of the Plasmodium life cycle which causes the symptoms of malaria. Furthermore, Hepatocystis is transmitted by biting midges, not mosquitoes. Comparative genomics of Hepatocystis and Plasmodium species therefore presents an opportunity to better understand some of the most important aspects of malaria parasite biology. We were able to generate a draft genome for Hepatocystis sp. using DNA sequencing reads from the blood of a naturally infected red colobus monkey. We provide robust phylogenetic support for Hepatocystis sp. as a sister group to Plasmodium parasites infecting rodents. We show transcriptomic support for a lack of replication in the blood and genomic support for a complete loss of a family of genes involved in red blood cell invasion. Our analyses highlight the rapid evolution of genes involved in parasite vector stages, revealing genes that may be critical for interactions between malaria parasites and mosquitoes.

Hepatocystis parasites are single-celled organisms, closely related to the Plasmodium species which cause malaria. But Hepatocystis are distinct–unlike Plasmodium they are thought not to replicate in the blood and cause little or no disease in their mammalian hosts. They are transmitted from one host to the next, not by mosquitoes, but by biting midges. In this study we generated a genome sequence for Hepatocystis–the first time this data has ever been produced and analysed for this species. We compared genome sequences of Hepatocystis and Plasmodium, confirming that Hepatocystis is descended from Plasmodium. We strengthened support for the absence of replication in the blood and, in line with this finding, discovered that genes involved in interaction with red blood cells have been lost in Hepatocystis. Our analyses revealed rapid evolution of genes which are active when the parasite is in the insect vector, highlighting those which might be important for understanding interaction between malaria parasites and mosquitoes. Hepatocystis has a fascinating evolutionary story and is a powerful comparator for understanding malaria parasite biology.

Refining the transcriptome of the human malaria parasite Plasmodium falciparum using amplification-free RNA-seq

Background

Plasmodium parasites undergo several major developmental transitions during their complex lifecycle, which are enabled by precisely ordered gene expression programs. Transcriptomes from the 48-h blood stages of the major human malaria parasite Plasmodium falciparum have been described using cDNA microarrays and RNA-seq, but these assays have not always performed well within non-coding regions, where the AT-content is often 90–95%.

Results

We developed a directional, amplification-free RNA-seq protocol (DAFT-seq) to reduce bias against AT-rich cDNA, which we have applied to three strains of P. falciparum (3D7, HB3 and IT). While strain-specific differences were detected, overall there is strong conservation between the transcriptional profiles. For the 3D7 reference strain, transcription was detected from 89% of the genome, with over 78% of the genome transcribed into mRNAs. We also find that transcription from bidirectional promoters frequently results in non-coding, antisense transcripts. These datasets allowed us to refine the 5′ and 3′ untranslated regions (UTRs), which can be variable, long (> 1000 nt), and often overlap those of adjacent transcripts.

Conclusions

The approaches applied in this study allow a refined description of the transcriptional landscape of P. falciparum and demonstrate that very little of the densely packed P. falciparum genome is inactive or redundant. By capturing the 5′ and 3′ ends of mRNAs, we reveal both constant and dynamic use of transcriptional start sites across the intraerythrocytic developmental cycle that will be useful in guiding the definition of regulatory regions for use in future experimental gene expression studies.

The chromatin bound proteome of the human malaria parasite

Proteins interacting with DNA are fundamental for mediating processes such as gene expression, DNA replication and maintenance of genome integrity. Accumulating evidence suggests that the chromatin of apicomplexan parasites, such as Plasmodium falciparum, is highly organized, and this structure provides an epigenetic mechanism for transcriptional regulation. To investigate how parasite chromatin structure is being regulated, we undertook comparative genomics analysis using 12 distinct eukaryotic genomes. We identified conserved and parasite-specific chromatin-associated domains (CADs) and proteins (CAPs). We then used the chromatin enrichment for proteomics (ChEP) approach to experimentally capture CAPs in P. falciparum. A topological scoring analysis of the proteomics dataset revealed stage-specific enrichments of CADs and CAPs. Finally, we characterized, two candidate CAPs: a conserved homologue of the structural maintenance of chromosome 3 protein and a homologue of the crowded-like nuclei protein, a plant-like protein functionally analogous to animal nuclear lamina proteins. Collectively, our results provide a comprehensive overview of CAPs in apicomplexans, and contribute to our understanding of the complex molecular components regulating chromatin structure and genome architecture in these deadly parasites.

Strains used in whole organism Plasmodium falciparum vaccine trials differ in genome structure, sequence, and immunogenic potential

BACKGROUND

Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions.

METHODS

Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia.

RESULTS

While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8+ T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region.

CONCLUSIONS

These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.

Genomic analysis informs malaria evolution

Updated analysis of Plasmodium falciparum malaria reveals markers of antimalarial resistance

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Progression of the canonical reference malaria parasite genome from 2002-2019

Here we describe the ways in which the sequence and annotation of the Plasmodium falciparum reference genome has changed since its publication in 2002. As the malaria species responsible for the most deaths worldwide, the richness of annotation and accuracy of the sequence are important resources for the P. falciparum research community as well as the basis for interpreting the genomes of subsequently sequenced species. At the time of publication in 2002 over 60% of predicted genes had unknown functions. As of March 2019, this number has been significantly decreased to 33%. The reduction is due to the inclusion of genes that were subsequently characterised experimentally and genes with significant similarity to others with known functions. In addition, the structural annotation of genes has been significantly refined; 27% of gene structures have been changed since 2002, comprising changes in exon-intron boundaries, addition or deletion of exons and the addition or deletion of genes. The sequence has also undergone significant improvements. In addition to the correction of a large number of single-base and insertion or deletion errors, a major miss-assembly between the subtelomeres of chromosome 7 and 8 has been corrected. As the number of sequenced isolates continues to grow rapidly, a single reference genome will not be an adequate basis for interpreting intra-species sequence diversity. We therefore describe in this publication a population reference genome of P. falciparum, called Pfref1. This reference will enable the community to map to regions that are not present in the current assembly. P. falciparum 3D7 will continue to be maintained, with ongoing curation ensuring continual improvements in annotation quality.

Progression of the canonical reference malaria parasite genome from 2002-2019

Here we describe the ways in which the sequence and annotation of the Plasmodium falciparum reference genome has changed since its publication in 2002. As the malaria species responsible for the most deaths worldwide, the richness of annotation and accuracy of the sequence are important resources for the P. falciparum research community as well as the basis for interpreting the genomes of subsequently sequenced species. At the time of publication in 2002 over 60% of predicted genes had unknown functions. As of March 2019, this number has been significantly decreased to 33%. The reduction is due to the inclusion of genes that were subsequently characterised experimentally and genes with significant similarity to others with known functions. In addition, the structural annotation of genes has been significantly refined; 27% of gene structures have been changed since 2002, comprising changes in exon-intron boundaries, addition or deletion of exons and the addition or deletion of genes. The sequence has also undergone significant improvements. In addition to the correction of a large number of single-base and insertion or deletion errors, a major miss-assembly between the subtelomeres of chromosome 7 and 8 has been corrected. As the number of sequenced isolates continues to grow rapidly, a single reference genome will not be an adequate basis for interpreting intra-species sequence diversity. We therefore describe in this publication a population reference genome of P. falciparum, called Pfref1. This reference will enable the community to map to regions that are not present in the current assembly. P. falciparum 3D7 will continue to be maintained, with ongoing curation ensuring continual improvements in annotation quality.