Genome sequence of the human malaria parasite Plasmodium falciparum

DNA sequence and chromatin differentiate sequence-specific transcription factor binding in the human malaria parasite Plasmodium falciparum

Development of the malaria parasite, Plasmodium falciparum, is regulated by a limited number of sequence-specific transcription factors (TFs). However, the mechanisms by which these TFs recognize genome-wide binding sites is largely unknown. To address TF specificity, we investigated the binding of two TF subsets that either bind CACACA or GTGCAC DNA sequence motifs and further characterized two additional ApiAP2 TFs, PfAP2-G and PfAP2-EXP, which bind unique DNA motifs (GTAC and TGCATGCA). We also interrogated the impact of DNA sequence and chromatin context on P. falciparum TF binding by integrating high-throughput in vitro and in vivo binding assays, DNA shape predictions, epigenetic post-translational modifications, and chromatin accessibility. We found that DNA sequence context minimally impacts binding site selection for paralogous CACACA-binding TFs, while chromatin accessibility, epigenetic patterns, co-factor recruitment, and dimerization correlate with differential binding. In contrast, GTGCAC-binding TFs prefer different DNA sequence context in addition to chromatin dynamics. Finally, we determined that TFs that preferentially bind divergent DNA motifs may bind overlapping genomic regions due to low-affinity binding to other sequence motifs. Our results demonstrate that TF binding site selection relies on a combination of DNA sequence and chromatin features, thereby contributing to the complexity of P. falciparum gene regulatory mechanisms.

Expression and biochemical characterization of the putative insulinase enzyme PF11_0189 found in the Plasmodium falciparum genome

PF11_0189 is a putative insulin degrading enzyme present in Plasmodium falciparum genome. The catalytic domain of PF11_0189 is about 27 kDa. Substrate specificity study shows PF11_0189 acts upon different types of proteins. The substrate specificity is found to be highest when insulin is used as a substrate. Metal dependency study shows highest dependency of PF11_0189 towards zinc metal for its proteolytic activity. Chelation of zinc metal with EDTA shows complete absence of PF11_0189 activity. Peptide inhibitors, P-70 and P-121 from combinatorial peptide library prepared against PF11_0189 show inhibition with an IC50 value of 4.8 μM and 7.5 μM respectively. A proven natural anti-malarial peptide cyclosporin A shows complete inhibition against PF11_0189 with an IC50 value of 0.75 μM suggesting PF11_0189 as a potential target for peptide inhibitors. The study implicates that PF11_0189 is a zinc metalloprotease involved in catalysis of insulin. The study gives a preliminary insight into the mechanism of complications arising from glucose abnormalities during severe malaria.

Ancient Plasmodium genomes shed light on the history of human malaria

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

Activation loop phosphorylation and cGMP saturation of PKG regulate egress of malaria parasites

The cGMP-dependent protein kinase (PKG) is the sole cGMP sensor in malaria parasites, acting as an essential signalling hub to govern key developmental processes throughout the parasite life cycle. Despite the importance of PKG in the clinically relevant asexual blood stages, many aspects of malarial PKG regulation, including the importance of phosphorylation, remain poorly understood. Here we use genetic and biochemical approaches to show that reduced cGMP binding to cyclic nucleotide binding domain B does not affect in vitro kinase activity but prevents parasite egress. Similarly, we show that phosphorylation of a key threonine residue (T695) in the activation loop is dispensable for kinase activity in vitro but is essential for in vivo PKG function, with loss of T695 phosphorylation leading to aberrant phosphorylation events across the parasite proteome and changes to the substrate specificity of PKG. Our findings indicate that Plasmodium PKG is uniquely regulated to transduce signals crucial for malaria parasite development.

RNA polymerase III is involved in regulating Plasmodium falciparum virulence

While often undetected and untreated, persistent seasonal asymptomatic malaria infections remain a global public health problem. Despite the presence of parasites in the peripheral blood, no symptoms develop. Disease severity is correlated with the levels of infected red blood cells (iRBCs) adhering within blood vessels. Changes in iRBC adhesion capacity have been linked to seasonal asymptomatic malaria infections, however how this is occurring is still unknown. Here, we present evidence that RNA polymerase III (RNA Pol III) transcription in Plasmodium falciparum is downregulated in field isolates obtained from asymptomatic individuals during the dry season. Through experiments with in vitro cultured parasites, we have uncovered an RNA Pol III-dependent mechanism that controls pathogen proliferation and expression of a major virulence factor in response to external stimuli. Our findings establish a connection between P. falciparum cytoadhesion and a non-coding RNA family transcribed by Pol III. Additionally, we have identified P. falciparum Maf1 as a pivotal regulator of Pol III transcription, both for maintaining cellular homeostasis and for responding adaptively to external signals. These results introduce a novel perspective that contributes to our understanding of P. falciparum virulence. Furthermore, they establish a connection between this regulatory process and the occurrence of seasonal asymptomatic malaria infections.

Deciphering the Plasmodium falciparum perinuclear var gene expression site

The protozoan parasite Plasmodium falciparum, responsible for the deadliest form of human malaria, employs antigenic variation via monoallelic expression as a key survival strategy. The selective activation of one out of the 60-member var gene family is key to understanding the parasite's ability to cause severe disease and evade the host immune response. var gene activation is initiated by its relocation to a specialized expression site. While the perinuclear expression site (PES) plays a crucial role in enabling the expression of a single allele, the characteristics of this PES remain largely obscure. Recent breakthroughs in genome editing tools and the discovery of regulatory noncoding RNAs have shed light on this intriguing biological feature, offering significant insights into the mechanisms of pathogen virulence.

PfHDAC1 is an essential regulator of P. falciparum asexual proliferation and host cell invasion genes with a dynamic genomic occupancy responsive to artemisinin stress

Plasmodium falciparum, the deadly protozoan parasite responsible for malaria, has a tightly regulated gene expression profile closely linked to its intraerythrocytic development cycle. Epigenetic modifiers of the histone acetylation code have been identified as key regulators of the parasite's transcriptome but require further investigation. In this study, we map the genomic distribution of Plasmodium falciparum histone deacetylase 1 (PfHDAC1) across the erythrocytic asexual development cycle and find it has a dynamic occupancy over a wide array of developmentally relevant genes. Overexpression of PfHDAC1 results in a progressive increment in parasite load over consecutive rounds of the asexual infection cycle and is associated with enhanced gene expression of multiple families of host cell invasion factors (merozoite surface proteins, rhoptry proteins, etc.) and with increased merozoite invasion efficiency. With the use of class-specific inhibitors, we demonstrate that PfHDAC1 activity in parasites is crucial for timely intraerythrocytic development. Interestingly, overexpression of PfHDAC1 results in decreased sensitivity to frontline-drug dihydroartemisinin in parasites. Furthermore, we identify that artemisinin exposure can interfere with PfHDAC1 abundance and chromatin occupancy, resulting in enrichment over genes implicated in response/resistance to artemisinin. Finally, we identify that dihydroartemisinin exposure can interrupt the in vitro catalytic deacetylase activity and post-translational phosphorylation of PfHDAC1, aspects that are crucial for its genomic function. Collectively, our results demonstrate PfHDAC1 to be a regulator of critical functions in asexual parasite development and host invasion, which is responsive to artemisinin exposure stress and deterministic of resistance to it.

IMPORTANCE

Malaria is a major public health problem, with the parasite Plasmodium falciparum causing most of the malaria-associated mortality. It is spread by the bite of infected mosquitoes and results in symptoms such as cyclic fever, chills, and headache. However, if left untreated, it can quickly progress to a more severe and life-threatening form. The World Health Organization currently recommends the use of artemisinin combination therapy, and it has worked as a gold standard for many years. Unfortunately, certain countries in southeast Asia and Africa, burdened with a high prevalence of malaria, have reported cases of drug-resistant infections. One of the major problems in controlling malaria is the emergence of artemisinin resistance. Population genomic studies have identified mutations in the Kelch13 gene as a molecular marker for artemisinin resistance. However, several reports thereafter indicated that Kelch13 is not the main mediator but rather hinted at transcriptional deregulation as a major determinant of drug resistance. Earlier, we identified PfGCN5 as a global regulator of stress-responsive genes, which are known to play a central role in artemisinin resistance generation. In this study, we have identified PfHDAC1, a histone deacetylase as a cell cycle regulator, playing an important role in artemisinin resistance generation. Taken together, our study identified key transcriptional regulators that play an important role in artemisinin resistance generation.

ULTRA-Effective Labeling of Repetitive Genomic Sequence

In the age of long read sequencing, genomics researchers now have access to accurate repetitive DNA sequence (including satellites) that, due to the limitations of short read sequencing, could previously be observed only as unmappable fragments. Tools that annotate repetitive sequence are now more important than ever, so that we can better understand newly uncovered repetitive sequences, and also so that we can mitigate errors in bioinformatic software caused by those repetitive sequences. To that end, we introduce the 1.0 release of our tool for identifying and annotating locally-repetitive sequence, ULTRA (ULTRA Locates Tandemly Repetitive Areas). ULTRA is fast enough to use as part of an efficient annotation pipeline, produces state-of-the-art reliable coverage of repetitive regions containing many mutations, and provides interpretable statistics and labels for repetitive regions. It released under an open license, and available for download at https://github.com/TravisWheelerLab/ULTRA.

Exploring malaria parasite surface proteins to devise highly immunogenic multi-epitope subunit vaccine for Plasmodium falciparum

BACKGROUND

Malaria has remained a major health concern for decades among people living in tropical and sub-tropical countries. Plasmodium falciparum is one of the critical species that cause severe malaria and is responsible for major mortality. Moreover, the parasite has generated resistance against all WHO recommended drugs and therapies. Therefore, there is an urgent need for preventive measures in the form of reliable vaccines to achieve the target of a malaria-free world. Surface proteins are the preferable choice for subunit vaccine development because they are rapidly detected and engaged by host immune cells and vaccination-induced antibodies. Additionally, abundant surface or membrane proteins may contribute to the opsonization of pathogens by vaccine-induced antibodies.

RESULTS

In our study, we have listed all those surface proteins from the literature that could be functionally important and essential for infection and immune evasion of the malaria parasite. Eight Plasmodium surface and membrane proteins from the pre-erythrocyte and erythrocyte stages were shortlisted. Thirty-seven epitopes (B-cell, CTL, and HTL epitopes) from these proteins were predicted using immune-informatic tools and joined with suitable peptide linkers to design a vaccine construct. A TLR-4 agonist peptide adjuvant was added at the N-terminus of the multi-epitope series, followed by the PADRE sequence and EAAAK linker. The TLR-4 receptor was docked with the construct's anticipated model structure. The complex of vaccine and TLR-4, with the lowest energy -1514, was found to be stable under simulated physiological settings.

CONCLUSION

This study has provided a novel multi-epitope construct that may be exploited further for the development of an efficient vaccine for malaria.

tRNA modification reprogramming contributes to artemisinin resistance in Plasmodium falciparum

Plasmodium falciparum artemisinin (ART) resistance is driven by mutations in kelch-like protein 13 (PfK13). Quiescence, a key aspect of resistance, may also be regulated by a yet unidentified epigenetic pathway. Transfer RNA modification reprogramming and codon bias translation is a conserved epitranscriptomic translational control mechanism that allows cells to rapidly respond to stress. We report a role for this mechanism in ART-resistant parasites by combining tRNA modification, proteomic and codon usage analyses in ring-stage ART-sensitive and ART-resistant parasites in response to drug. Post-drug, ART-resistant parasites differentially hypomodify mcm5s2U on tRNA and possess a subset of proteins, including PfK13, that are regulated by Lys codon-biased translation. Conditional knockdown of the terminal s2U thiouridylase, PfMnmA, in an ART-sensitive parasite background led to increased ART survival, suggesting that hypomodification can alter the parasite ART response. This study describes an epitranscriptomic pathway via tRNA s2U reprogramming that ART-resistant parasites may employ to survive ART-induced stress.

A metagenomic analysis of the phase 2 Anopheles gambiae 1000 genomes dataset reveals a wide diversity of cobionts associated with field collected mosquitoes

The Anopheles gambiae 1000 Genomes (Ag1000G) Consortium previously utilized deep sequencing methods to catalogue genetic diversity across African An. gambiae populations. We analyzed the complete datasets of 1142 individually sequenced mosquitoes through Microsoft Premonition's Bayesian mixture model based (BMM) metagenomics pipeline. All specimens were confirmed as either An. gambiae sensu stricto (s.s.) or An. coluzzii with a high degree of confidence ( > 98% identity to reference). Homo sapiens DNA was identified in all specimens indicating contamination may have occurred either at the time of specimen collection, preparation and/or sequencing. We found evidence of vertebrate hosts in 162 specimens. 59 specimens contained validated Plasmodium falciparum reads. Human hepatitis B and primate erythroparvovirus-1 viral sequences were identified in fifteen and three mosquito specimens, respectively. 478 of the 1,142 specimens were found to contain bacterial reads and bacteriophage-related contigs were detected in 27 specimens. This analysis demonstrates the capacity of metagenomic approaches to elucidate important vector-host-pathogen interactions of epidemiological significance.

From multiplicity of infection to force of infection for sparsely sampled Plasmodium falciparum populations at high transmission

High multiplicity of infection or MOI, the number of genetically distinct parasite strains co-infecting a single human host, characterizes infectious diseases including falciparum malaria at high transmission. It accompanies high asymptomatic Plasmodium falciparum prevalence despite high exposure, creating a large transmission reservoir challenging intervention. High MOI and asymptomatic prevalence are enabled by immune evasion of the parasite achieved via vast antigenic diversity. Force of infection or FOI, the number of new infections acquired by an individual host over a given time interval, is the dynamic sister quantity of MOI, and a key epidemiological parameter for monitoring the impact of antimalarial interventions and assessing vaccine or drug efficacy in clinical trials. FOI remains difficult, expensive, and labor-intensive to accurately measure, especially in high-transmission regions, whether directly via cohort studies or indirectly via the fitting of epidemiological models to repeated cross-sectional surveys. We propose here the application of queuing theory to obtain FOI on the basis of MOI, in the form of either a two-moment approximation method or Little's law. We illustrate these methods with MOI estimates obtained under sparse sampling schemes with the recently proposed " v a r coding" method, based on sequences of the v a r multigene family encoding for the major variant surface antigen of the blood stage of malaria infection. The methods are evaluated with simulation output from a stochastic agent-based model, and are applied to an interrupted time-series study from Bongo District in northern Ghana before and immediately after a three-round transient indoor residual spraying (IRS) intervention. We incorporate into the sampling of the simulation output, limitations representative of those encountered in the collection of field data, including under-sampling of v a r genes, missing data, and usage of antimalarial drug treatment. We address these limitations in MOI estimates with a Bayesian framework and an imputation bootstrap approach. We demonstrate that both proposed methods give good and consistent FOI estimates across various simulated scenarios. Their application to the field surveys shows a pronounced reduction in annual FOI during intervention, of more than 70%. The proposed approach should be applicable to the many geographical locations where cohort or cross-sectional studies with regular and frequent sampling are lacking but single-time-point surveys under sparse sampling schemes are available, and for MOI estimates obtained in different ways. They should also be relevant to other pathogens of humans, wildlife and livestock whose immune evasion strategies are based on large antigenic variation resulting in high multiplicity of infection.

Size-dependent enhancement of gene expression by Plasmodium 5'UTR introns

BACKGROUND

Eukaryotic genes contain introns that are removed by the spliceosomal machinery during mRNA maturation. Introns impose a huge energetic burden on a cell; therefore, they must play an essential role in maintaining genome stability and/or regulating gene expression. Many genes (> 50%) in Plasmodium parasites contain predicted introns, including introns in 5' and 3' untranslated regions (UTR). However, the roles of UTR introns in the gene expression of malaria parasites remain unknown.

METHODS

In this study, an episomal dual-luciferase assay was developed to evaluate gene expression driven by promoters with or without a 5'UTR intron from four Plasmodium yoelii genes. To investigate the effect of the 5'UTR intron on endogenous gene expression, the pytctp gene was tagged with 3xHA at the N-terminal of the coding region, and parasites with or without the 5'UTR intron were generated using the CRISPR/Cas9 system.

RESULTS

We showed that promoters with 5'UTR introns had higher activities in driving gene expression than those without 5'UTR introns. The results were confirmed in recombinant parasites expressing an HA-tagged gene (pytctp) driven by promoter with or without 5'UTR intron. The enhancement of gene expression was intron size dependent, but not the DNA sequence, e.g. the longer the intron, the higher levels of expression. Similar results were observed when a promoter from one strain of P. yoelii was introduced into different parasite strains. Finally, the 5'UTR introns were alternatively spliced in different parasite development stages, suggesting an active mechanism employed by the parasites to regulate gene expression in various developmental stages.

CONCLUSIONS

Plasmodium 5'UTR introns enhance gene expression in a size-dependent manner; the presence of alternatively spliced mRNAs in different parasite developmental stages suggests that alternative slicing of 5'UTR introns is one of the key mechanisms in regulating parasite gene expression and differentiation.

Insights into the Mechanism of Catalytic Activity of Plasmodium Parasite Malate-Quinone Oxidoreductase

Plasmodium malate-quinone oxidoreductase (MQO) is a membrane flavoprotein catalyzing the oxidation of malate to oxaloacetate and the reduction of quinone to quinol. Recently, using a yeast expression system, we demonstrated that MQO, expressed in place of mitochondrial malate dehydrogenase (MDH), contributes to the TCA cycle and the electron transport chain in mitochondria, making MQO attractive as a promising drug target in Plasmodium malaria parasites, which lack mitochondrial MDH. However, there is little information on the structure of MQO and its catalytic mechanism, information that will be required to develop novel drugs. Here, we investigated the catalytic site of P. falciparum MQO (PfMQO) using our yeast expression system. We generated a model structure for PfMQO with the AI tool AlphaFold and used protein footprinting by acetylation with acetic anhydride to analyze the surface topology of the model, confirming the computational prediction to be reasonably accurate. Moreover, a putative catalytic site, which includes a possible flavin-binding site, was identified by this combination of protein footprinting and structural prediction model. This active site was analyzed by site-directed mutagenesis. By measuring enzyme activity and protein expression levels in the PfMQO mutants, we showed that several residues at the active site are essential for enzyme function. In addition, a single substitution mutation near the catalytic site resulted in enhanced sensitivity to ferulenol, an inhibitor of PfMQO that competes with malate for binding to the enzyme. This strongly supports the notion that the substrate binds to the proposed catalytic site. Then, the location of the catalytic site was demonstrated by structural comparison with a homologous enzyme. Finally, we used our results to propose a mechanism for the catalytic activity of MQO by reference to the mechanism of action of structurally or functionally homologous enzymes.

Diversity and selection analyses identify transmission-blocking antigens as the optimal vaccine candidates in Plasmodium falciparum

Background

A highly effective vaccine for malaria remains an elusive target, at least in part due to the under-appreciated natural parasite variation. This study aimed to investigate genetic and structural variation, and immune selection of leading malaria vaccine candidates across the Plasmodium falciparum's life cycle.

Methods

We analyzed 325 P. falciparum whole genome sequences from Zambia, in addition to 791 genomes from five other African countries available in the MalariaGEN Pf3k Rdatabase. Ten vaccine antigens spanning three life-history stages were examined for genetic and structural variations, using population genetics measures, haplotype network analysis, and 3D structure selection analysis.

Findings

Among the ten antigens analyzed, only three in the transmission-blocking vaccine category display P. falciparum 3D7 as the dominant haplotype. The antigens AMA1, CSP, MSP119 and CelTOS, are much more diverse than the other antigens, and their epitope regions are under moderate to strong balancing selection. In contrast, Rh5, a blood stage antigen, displays low diversity yet slightly stronger immune selection in the merozoite-blocking epitope region. Except for CelTOS, the transmission-blocking antigens Pfs25, Pfs48/45, Pfs230, Pfs47, and Pfs28 exhibit minimal diversity and no immune selection in epitopes that induce strain-transcending antibodies, suggesting potential effectiveness of 3D7-based vaccines in blocking transmission.

Interpretations

These findings offer valuable insights into the selection of optimal vaccine candidates against P. falciparum. Based on our results, we recommend prioritizing conserved merozoite antigens and transmission-blocking antigens. Combining these antigens in multi-stage approaches may be particularly promising for malaria vaccine development initiatives.

Funding

Purdue Department of Biological Sciences; Puskas Memorial Fellowship; National Institute of Allergy and Infectious Diseases (U19AI089680).

Charting new territory: The Plasmodium falciparum tRNA modification landscape

Ribonucleoside modifications comprising the epitranscriptome are present in all organisms and all forms of RNA, including mRNA, rRNA and tRNA, the three major RNA components of the translational machinery. Of these, tRNA is the most heavily modified and the tRNA epitranscriptome has the greatest diversity of modifications. In addition to their roles in tRNA biogenesis, quality control, structure, cleavage, and codon recognition, tRNA modifications have been shown to regulate gene expression post-transcriptionally in prokaryotes and eukaryotes, including humans. However, studies investigating the impact of tRNA modifications on gene expression in the malaria parasite Plasmodium falciparum are currently scarce. Current evidence shows that the parasite has a limited capacity for transcriptional control, which points to a heavier reliance on strategies for posttranscriptional regulation such as tRNA epitranscriptome reprogramming. This review addresses the known functions of tRNA modifications in the biology of P. falciparum while highlighting the potential therapeutic opportunities and the value of using P. falciparum as a model organism for addressing several open questions related to the tRNA epitranscriptome.

Benchmarking and Optimization of Methods for the Detection of Identity-By-Descent in Plasmodium falciparum

Genomic surveillance is crucial for identifying at-risk populations for targeted malaria control and elimination. Identity-by-descent (IBD) is being used in Plasmodium population genomics to estimate genetic relatedness, effective population size ( N e ) , population structure, and positive selection. However, a comprehensive evaluation of IBD segment detection tools is lacking for species with high rates of recombination. Here, we employ genetic simulations reflecting P. falciparum's high recombination rate and decreasing N e to benchmark IBD callers, including probabilistic (hmmIBD, isoRelate), identity-by-state-based (hap-IBD, phased IBD) and others (Refined IBD), using genealogy-based true IBD and downstream inference of population characteristics. Our findings reveal that low marker density per genetic unit, related to high recombination rates relative to mutation rates, significantly affects the quality of detected IBD segments. Most IBD callers suffer from high false negative rates, which can be improved with parameter optimization. Optimized parameters allow for more accurate capture of selection signals and population structure, but hmmIBD is unique in providing less biased estimates of N e . Empirical data subsampled from the MalariaGEN Pf7 database, representing different transmission settings, confirmed these patterns. We conclude that the detection of IBD in high-recombining species requires context-specific evaluation and parameter optimization and recommend that hmmIBD be used for quality-sensitive analysis, such as estimation of N e in these species.

Biomolecular interactions between Plasmodium and human host: A basis of targeted antimalarial therapy

Malaria is one of the serious health concerns worldwide as it remains a clinical challenge due to the complex life cycle of the malaria parasite and the morphological changes it undergoes during infection. The malaria parasite multiplies rapidly and spreads in the population by changing its alternative hosts. These various morphological stages of the parasite in the human host cause clinical symptoms (anemia, fever, and coma). These symptoms arise due to the preprogrammed biology of the parasite in response to the human pathophysiological response. Thus, complete elimination becomes one of the major health challenges. Although malaria vaccine(s) are available in the market, they still contain to cause high morbidity and mortality. Therefore, an approach for eradication is needed through the exploration of novel molecular targets by tracking the epidemiological changes the parasite adopts. This review focuses on the various novel molecular targets.

Positive-unlabeled learning identifies vaccine candidate antigens in the malaria parasite Plasmodium falciparum

Malaria vaccine development is hampered by extensive antigenic variation and complex life stages of Plasmodium species. Vaccine development has focused on a small number of antigens, many of which were identified without utilizing systematic genome-level approaches. In this study, we implement a machine learning-based reverse vaccinology approach to predict potential new malaria vaccine candidate antigens. We assemble and analyze P. falciparum proteomic, structural, functional, immunological, genomic, and transcriptomic data, and use positive-unlabeled learning to predict potential antigens based on the properties of known antigens and remaining proteins. We prioritize candidate antigens based on model performance on reference antigens with different genetic diversity and quantify the protein properties that contribute most to identifying top candidates. Candidate antigens are characterized by gene essentiality, gene ontology, and gene expression in different life stages to inform future vaccine development. This approach provides a framework for identifying and prioritizing candidate vaccine antigens for a broad range of pathogens.

Plasmodium falciparum Alba6 exhibits DNase activity and participates in stress response

Alba domain proteins, owing to their functional plasticity, play a significant role in organisms. Here, we report an intrinsic DNase activity of PfAlba6 from Plasmodium falciparum, an etiological agent responsible for human malignant malaria. We identified that tyrosine28 plays a critical role in the Mg2+ driven 5'-3' DNase activity of PfAlba6. PfAlba6 cleaves both dsDNA as well as ssDNA. We also characterized PfAlba6-DNA interaction and observed concentration-dependent oligomerization in the presence of DNA, which is evident from size exclusion chromatography and single molecule AFM-imaging. PfAlba6 mRNA expression level is up-regulated several folds following heat stress and treatment with artemisinin, indicating a possible role in stress response. PfAlba6 has no human orthologs and is expressed in all intra-erythrocytic stages; thus, this protein can potentially be a new anti-malarial drug target.

Genomics reveals heterogeneous Plasmodium falciparum transmission and selection signals in Zambia

BACKGROUND

Genomic surveillance is crucial for monitoring malaria transmission and understanding parasite adaptation to interventions. Zambia lacks prior nationwide efforts in malaria genomic surveillance among African countries.

METHODS

We conducted genomic surveillance of Plasmodium falciparum parasites from the 2018 Malaria Indicator Survey in Zambia, a nationally representative household survey of children under five years of age. We whole-genome sequenced and analyzed 241 P. falciparum genomes from regions with varying levels of malaria transmission across Zambia and estimated genetic metrics that are informative about transmission intensity, genetic relatedness between parasites, and selection.

RESULTS

We provide genomic evidence of widespread within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in Zambia. Our analysis reveals country-level clustering of parasites from Zambia and neighboring regions, with distinct separation in West Africa. Within Zambia, identity by descent (IBD) relatedness analysis uncovers local spatial clustering and rare cases of long-distance sharing of closely related parasite pairs. Genomic regions with large shared IBD segments and strong positive selection signatures implicate genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Furthermore, differences in selection signatures, including drug resistance loci, are observed between eastern and western Zambian parasite populations, suggesting variable transmission intensity and ongoing drug pressure.

CONCLUSIONS

Our findings enhance our understanding of nationwide P. falciparum transmission in Zambia, establishing a baseline for analyzing parasite genetic metrics as they vary over time and space. These insights highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Identifying the best PCR enzyme for library amplification in NGS

Background. PCR amplification is a necessary step in many next-generation sequencing (NGS) library preparation methods [1, 2]. Whilst many PCR enzymes are developed to amplify single targets efficiently, accurately and with specificity, few are developed to meet the challenges imposed by NGS PCR, namely unbiased amplification of a wide range of different sizes and GC content. As a result PCR amplification during NGS library prep often results in bias toward GC neutral and smaller fragments. As NGS has matured, optimized NGS library prep kits and polymerase formulations have emerged and in this study we have tested a wide selection of available enzymes for both short-read Illumina library preparation and long fragment amplification ahead of long-read sequencing.We tested over 20 different hi-fidelity PCR enzymes/NGS amplification mixes on a range of Illumina library templates of varying GC content and composition, and find that both yield and genome coverage uniformity characteristics of the commercially available enzymes varied dramatically. Three enzymes Quantabio RepliQa Hifi Toughmix, Watchmaker Library Amplification Hot Start Master Mix (2X) 'Equinox' and Takara Ex Premier were found to give a consistent performance, over all genomes, that mirrored closely that observed for PCR-free datasets. We also test a range of enzymes for long-read sequencing by amplifying size fractionated S. cerevisiae DNA of average size 21.6 and 13.4 kb, respectively.The enzymes of choice for short-read (Illumina) library fragment amplification are Quantabio RepliQa Hifi Toughmix, Watchmaker Library Amplification Hot Start Master Mix (2X) 'Equinox' and Takara Ex Premier, with RepliQa also being the best performing enzyme from the enzymes tested for long fragment amplification prior to long-read sequencing.

Two DOT1 enzymes cooperatively mediate efficient ubiquitin-independent histone H3 lysine 76 tri-methylation in kinetoplastids

In higher eukaryotes, a single DOT1 histone H3 lysine 79 (H3K79) methyltransferase processively produces H3K79me2/me3 through histone H2B mono-ubiquitin interaction, while the kinetoplastid Trypanosoma brucei di-methyltransferase DOT1A and tri-methyltransferase DOT1B efficiently methylate the homologous H3K76 without H2B mono-ubiquitination. Based on structural and biochemical analyses of DOT1A, we identify key residues in the methyltransferase motifs VI and X for efficient ubiquitin-independent H3K76 methylation in kinetoplastids. Substitution of a basic to an acidic residue within motif VI (Gx6K) is essential to stabilize the DOT1A enzyme-substrate complex, while substitution of the motif X sequence VYGE by CAKS renders a rigid active-site loop flexible, implying a distinct mechanism of substrate recognition. We further reveal distinct methylation kinetics and substrate preferences of DOT1A (H3K76me0) and DOT1B (DOT1A products H3K76me1/me2) in vitro, determined by a Ser and Ala residue within motif IV, respectively, enabling DOT1A and DOT1B to mediate efficient H3K76 tri-methylation non-processively but cooperatively, and suggesting why kinetoplastids have evolved two DOT1 enzymes.

Evaluating the performance of Plasmodium falciparum genetic metrics for inferring National Malaria Control Programme reported incidence in Senegal

BACKGROUND

Genetic surveillance of the Plasmodium falciparum parasite shows great promise for helping National Malaria Control Programmes (NMCPs) assess parasite transmission. Genetic metrics such as the frequency of polygenomic (multiple strain) infections, genetic clones, and the complexity of infection (COI, number of strains per infection) are correlated with transmission intensity. However, despite these correlations, it is unclear whether genetic metrics alone are sufficient to estimate clinical incidence.

METHODS

This study examined parasites from 3147 clinical infections sampled between the years 2012-2020 through passive case detection (PCD) across 16 clinic sites spread throughout Senegal. Samples were genotyped with a 24 single nucleotide polymorphism (SNP) molecular barcode that detects parasite strains, distinguishes polygenomic (multiple strain) from monogenomic (single strain) infections, and identifies clonal infections. To determine whether genetic signals can predict incidence, a series of Poisson generalized linear mixed-effects models were constructed to predict the incidence level at each clinical site from a set of genetic metrics designed to measure parasite clonality, superinfection, and co-transmission rates.

RESULTS

Model-predicted incidence was compared with the reported standard incidence data determined by the NMCP for each clinic and found that parasite genetic metrics generally correlated with reported incidence, with departures from expected values at very low annual incidence (< 10/1000/annual [‰]).

CONCLUSIONS

When transmission is greater than 10 cases per 1000 annual parasite incidence (annual incidence > 10‰), parasite genetics can be used to accurately infer incidence and is consistent with superinfection-based hypotheses of malaria transmission. When transmission was < 10‰, many of the correlations between parasite genetics and incidence were reversed, which may reflect the disproportionate impact of importation and focal transmission on parasite genetics when local transmission levels are low.

Expression and purification of active shikimate dehydrogenase from Plasmodium falciparum

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

Novel Ion Channel Genes in Malaria Parasites

Ion channels serve many cellular functions including ion homeostasis, volume regulation, signaling, nutrient acquisition, and developmental progression. Although the complex life cycles of malaria parasites necessitate ion and solute flux across membranes, the whole-genome sequencing of the human pathogen Plasmodium falciparum revealed remarkably few orthologs of known ion channel genes. Contrasting with this, biochemical studies have implicated the channel-mediated flux of ions and nutritive solutes across several membranes in infected erythrocytes. Here, I review advances in the cellular and molecular biology of ion channels in malaria parasites. These studies have implicated novel parasite genes in the formation of at least two ion channels, with additional ion channels likely present in various membranes and parasite stages. Computational approaches that rely on homology to known channel genes from higher organisms will not be very helpful in identifying the molecular determinants of these activities. Given their unusual properties, novel molecular and structural features, and essential roles in pathogen survival and development, parasite channels should be promising targets for therapy development.

Flexible and cost-effective genomic surveillance of P. falciparum malaria with targeted nanopore sequencing

Genomic surveillance of Plasmodium falciparum malaria can provide policy-relevant information about antimalarial drug resistance, diagnostic test failure, and the evolution of vaccine targets. Yet the large and low complexity genome of P. falciparum complicates the development of genomic methods, while resource constraints in malaria endemic regions can limit their deployment. Here, we demonstrate an approach for targeted nanopore sequencing of P. falciparum from dried blood spots (DBS) that enables cost-effective genomic surveillance of malaria in low-resource settings. We release software that facilitates flexible design of amplicon sequencing panels and use this software to design two target panels for P. falciparum. The panels generate 3-4 kbp reads for eight and sixteen targets respectively, covering key drug-resistance associated genes, diagnostic test antigens, polymorphic markers and the vaccine target csp. We validate our approach on mock and field samples, demonstrating robust sequencing coverage, accurate variant calls within coding sequences, the ability to explore P. falciparum within-sample diversity and to detect deletions underlying rapid diagnostic test failure.

Plasmodium falciparum topoisomerases: Emerging targets for anti-malarial therapy

Different metabolic pathways like DNA replication, transcription, and recombination generate topological constrains in the genome. These topological constraints are resolved by essential molecular machines known as topoisomerases. To bring changes in DNA topology, the topoisomerases create a single or double-stranded nick in the template DNA, hold the nicked ends to let the tangled DNA pass through, and finally re-ligate the breaks. The DNA nicking and re-ligation activities as well as ATPase activities (when present) in topoisomerases are subjected to inhibition by several anticancer and antibacterial drugs, thus establishing these enzymes as successful targets in anticancer and antibacterial therapies. The anti-topoisomerase drugs interfere with the functioning of these enzymes and result in the accumulation of DNA tangles or lethal genomic breaks, thereby promoting host cell (or organism) death. The potential of topoisomerases in the human malarial parasite, Plasmodium falciparum in antimalarial drug development has received little attention so far. Interestingly, the parasite genome encodes orthologs of topoisomerases found in eukaryotes, prokaryotes, and archaea, thus, providing an enormous opportunity for investigating these enzymes for antimalarial therapeutics. This review focuses on the features of Plasmodium falciparum topoisomerases (PfTopos) with respect to their closer counterparts in other organisms. We will discuss overall advances and basic challenges with topoisomerase research in Plasmodium falciparum and our attempts to understand the interaction of PfTopos with classical and new-generation topoisomerase inhibitors using in silico molecular docking approach. The recent episodes of parasite resistance against artemisinin, the only effective antimalarial drug at present, further highlight the significance of investigating new drug targets including topoisomerases in antimalarial therapeutics.

Using a mobile nanopore sequencing lab for end-to-end genomic surveillance of Plasmodium falciparum: A feasibility study

Genomic epidemiology holds promise for malaria control and elimination efforts, for example by informing on Plasmodium falciparum genetic diversity and prevalence of mutations conferring anti-malarial drug resistance. Limited sequencing infrastructure in many malaria-endemic areas prevents the rapid generation of genomic data. To address these issues, we developed and validated assays for P. falciparum nanopore sequencing in endemic sites using a mobile laboratory, targeting key antimalarial drug resistance markers and microhaplotypes. Using two multiplexed PCR reactions, we amplified six highly polymorphic microhaplotypes and ten drug resistance markers. We developed a bioinformatics workflow that allows genotyping of polyclonal malaria infections, including minority clones. We validated the panels on mock dried blood spot (DBS) and rapid diagnostic test (RDT) samples and archived DBS, demonstrating even, high read coverage across amplicons (range: 580x to 3,212x median coverage), high haplotype calling accuracy, and the ability to explore within-sample diversity of polyclonal infections. We field-tested the feasibility of rapid genotyping in Zanzibar in close collaboration with the local malaria elimination program using DBS and routinely collected RDTs as sample inputs. Our assay identified haplotypes known to confer resistance to known antimalarials in the dhfr, dhps and mdr1 genes, but no evidence of artemisinin partial resistance. Most infections (60%) were polyclonal, with high microhaplotype diversity (median HE = 0.94). In conclusion, our assays generated actionable data within a few days, and we identified current challenges for implementing nanopore sequencing in endemic countries to accelerate malaria control and elimination.

Plasmodium falciparum contains functional SCF and CRL4 ubiquitin E3 ligases, and CRL4 is critical for cell division and membrane integrity

Protein ubiquitination is essential for cellular homeostasis and regulation of several processes, including cell division and genome integrity. Ubiquitin E3 ligases determine substrate specificity for ubiquitination, and Cullin-RING E3 ubiquitin ligases (CRLs) make the largest group among the ubiquitin E3 ligases. Although conserved and most studied in model eukaryotes, CRLs remain underappreciated in Plasmodium and related parasites. To investigate the CRLs of human malaria parasite Plasmodium falciparum, we generated parasites expressing tagged P. falciparum cullin-1 (PfCullin-1), cullin-2 (PfCullin-2), Rbx1 (PfRbx1) and Skp1 (PfSkp1). PfCullin-1 and PfCullin-2 were predominantly expressed in erythrocytic trophozoite and schizont stages, with nucleocytoplasmic localization and chromatin association, suggesting their roles in different cellular compartments and DNA-associated processes. Immunoprecipitation, in vitro protein-protein interaction, and ubiquitination assay confirmed the presence of a functional Skp1-Cullin-1-Fbox (PfSCF) complex, comprising of PfCullin-1, PfRbx1, PfSkp1, PfFBXO1, and calcyclin binding protein. Immunoprecipitation, sequence analysis, and ubiquitination assay indicated that PfCullin-2 forms a functional human CRL4-like complex (PfCRL4), consisting of PfRbx1, cleavage and polyadenylation specificity factor subunit_A and WD40 repeat proteins. PfCullin-2 knock-down at the protein level, which would hinder PfCRL4 assembly, significantly decreased asexual and sexual erythrocytic stage development. The protein levels of several pathways, including protein translation and folding, lipid biosynthesis and transport, DNA replication, and protein degradation were significantly altered upon PfCullin-2 depletion, which likely reflects association of PfCRL4 with multiple pathways. PfCullin-2-depleted schizonts had poorly delimited merozoites and internal membraned structures, suggesting a role of PfCRL4 in maintaining membrane integrity. PfCullin-2-depleted parasites had a significantly lower number of nuclei/parasite than the normal parasites, indicating a crucial role of PfCRL4 in cell division. We demonstrate the presence of functional CRLs in P. falciparum, with crucial roles for PfCRL4 in cell division and maintaining membrane integrity.

Long-Read Sequencing and De Novo Genome Assembly Pipeline of Two Plasmodium falciparum Clones (Pf3D7, PfW2) Using Only the PromethION Sequencer from Oxford Nanopore Technologies without Whole-Genome Amplification

Antimalarial drug resistance has become a real public health problem despite WHO measures. New sequencing technologies make it possible to investigate genomic variations associated with resistant phenotypes at the genome-wide scale. Based on the use of hemisynthetic nanopores, the PromethION technology from Oxford Nanopore Technologies can produce long-read sequences, in contrast to previous short-read technologies used as the gold standard to sequence Plasmodium. Two clones of P. falciparum (Pf3D7 and PfW2) were sequenced in long-read using the PromethION sequencer from Oxford Nanopore Technologies without genomic amplification. This made it possible to create a processing analysis pipeline for human Plasmodium with ONT Fastq only. De novo assembly revealed N50 lengths of 18,488 kb and 17,502 kb for the Pf3D7 and PfW2, respectively. The genome size was estimated at 23,235,407 base pairs for the Pf3D7 clone and 21,712,038 base pairs for the PfW2 clone. The average genome coverage depth was estimated at 787X and 653X for the Pf3D7 and PfW2 clones, respectively. This study proposes an assembly processing pipeline for the human Plasmodium genome using software adapted to large ONT data and the high AT percentage of Plasmodium. This search provides all the parameters which were optimized for use with the software selected in the pipeline.

The genetic landscape of origins of replication in P. falciparum

Various origin mapping approaches have enabled genome-wide identification of origins of replication (ORI) in model organisms, but only a few studies have focused on divergent organisms. By employing three complementary approaches we provide a high-resolution map of ORIs in Plasmodium falciparum, the deadliest human malaria parasite. We profiled the distribution of origin of recognition complex (ORC) binding sites by ChIP-seq of two PfORC subunits and mapped active ORIs using NFS and SNS-seq. We show that ORIs lack sequence specificity but are not randomly distributed, and group in clusters. Licensing is biased towards regions of higher GC content and associated with G-quadruplex forming sequences (G4FS). While strong transcription likely enhances firing, active origins are depleted from transcription start sites. Instead, most accumulate in transcriptionally active gene bodies. Single molecule analysis of nanopore reads containing multiple initiation events, which could have only come from individual nuclei, showed a relationship between the replication fork pace and the distance to the nearest origin. While some similarities were drawn with the canonic eukaryote model, the distribution of ORIs in P. falciparum is likely shaped by unique genomic features such as extreme AT-richness-a product of evolutionary pressure imposed by the parasitic lifestyle.

CD8+ Trms against malaria liver-stage: prospects and challenges

Attenuated sporozoites provide a valuable model for exploring protective immunity against the malarial liver stage, guiding the design of highly efficient vaccines to prevent malaria infection. Liver tissue-resident CD8+ T cells (CD8+ Trm cells) are considered the host front-line defense against malaria and are crucial to developing prime-trap/target strategies for pre-erythrocytic stage vaccine immunization. However, the spatiotemporal regulatory mechanism of the generation of liver CD8+ Trm cells and their responses to sporozoite challenge, as well as the protective antigens they recognize remain largely unknown. Here, we discuss the knowledge gap regarding liver CD8+ Trm cell formation and the potential strategies to identify predominant protective antigens expressed in the exoerythrocytic stage, which is essential for high-efficacy malaria subunit pre-erythrocytic vaccine designation.

Large DNA fragment knock-in and sequential gene editing in Plasmodium falciparum: a preliminary study using suicide-rescue-based CRISPR/Cas9 system

CRISPR/Cas9 technology applied to Plasmodium falciparum offers the potential to greatly improve gene editing, but such expectations including large DNA fragment knock-ins and sequential gene editing have remained unfulfilled. Here, we achieved a major advance in addressing this challenge, especially for creating large DNA fragment knock-ins and sequential editing, by modifying our suicide-rescue-based system that has already been demonstrated to be highly efficient for conventional gene editing. This improved approach was confirmed to mediate efficient knock-ins of DNA fragments up to 6.3 kb, to produce "marker-free" genetically engineered parasites and to show potential for sequential gene editing. This represents an important advancement in establishing platforms for large-scale genome editing, which might gain a better understanding of gene function for the most lethal cause of malaria and contribute to adjusting synthetic biology strategies to live parasite malaria vaccine development. Site-directed knock-in of large DNA fragments is highly efficient using suicide-rescue-based CRISPR/Cas9 system, and sequential gene insertion is feasible but further confirmation is still needed.

In silico modeling revealed phytomolecules derived from Cymbopogon citratus (DC.) leaf extract as promising candidates for malaria therapy

The emergence of varying levels of resistance to currently available antimalarial drugs significantly threatens global health. This factor heightens the urgency to explore bioactive compounds from natural products with a view to discovering and developing newer antimalarial drugs with novel mode of actions. Therefore, we evaluated the inhibitory effects of sixteen phytocompounds from Cymbopogon citratus leaf extract against Plasmodium falciparum drug targets such as P. falciparum circumsporozoite protein (PfCSP), P. falciparum merozoite surface protein 1 (PfMSP1) and P. falciparum erythrocyte membrane protein 1 (PfEMP1). In silico approaches including molecular docking, pharmacophore modeling and 3D-QSAR were adopted to analyze the inhibitory activity of the compounds under consideration. The molecular docking results indicated that a compound swertiajaponin from C. citratus exhibited a higher binding affinity (-7.8 kcal/mol) to PfMSP1 as against the standard artesunate-amodiaquine (-6.6 kcal/mol). Swertiajaponin also formed strong hydrogen bond interactions with LYS29, CYS30, TYR34, ASN52, GLY55 and CYS28 amino acid residues. In addition, quercetin another compound from C. citratus exhibited significant binding energies -6.8 and -8.3 kcal/mol with PfCSP and PfEMP1, respectively but slightly lower than the standard artemether-lumefantrine with binding energies of -7.4 kcal/mol against PfCSP and -8.7 kcal/mol against PfEMP1. Overall, the present study provides evidence that swertiajaponin and other phytomolecules from C. citratus have modulatory properties toward P. falciparum drug targets and thus may warrant further exploration in early drug discovery efforts against malaria. Furthermore, these findings lend credence to the folkloric use of C. citratus for malaria treatment.Communicated by Ramaswamy H. Sarma.

Plasmodium, the Apicomplexa Outlier When It Comes to Protein Synthesis

Plasmodium is an obligate intracellular parasite that has numerous interactions with different hosts during its elaborate life cycle. This is also the case for the other parasites belonging to the same phylum Apicomplexa. In this study, we bioinformatically identified the components of the multi-synthetase complexes (MSCs) of several Apicomplexa parasites and modelled their assembly using AlphaFold2. It appears that none of these MSCs resemble the two MSCs that we have identified and characterized in Plasmodium. Indeed, tRip, the central protein involved in the association of the two Plasmodium MSCs is different from its homologues, suggesting also that the tRip-dependent import of exogenous tRNAs is not conserved in other apicomplexan parasites. Based on this observation, we searched for obvious differences that could explain the singularity of Plasmodium protein synthesis by comparing tRNA genes and amino acid usage in the different genomes. We noted a contradiction between the large number of asparagine residues used in Plasmodium proteomes and the single gene encoding the tRNA that inserts them into proteins. This observation remains true for all the Plasmodia strains studied, even those that do not contain long asparagine homorepeats.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

The protist Aurantiochytrium has universal subtelomeric rDNAs and is a host for mirusviruses

Viruses are the most abundant biological entities in the world's oceans, where they play important ecological and biogeochemical roles. Metagenomics is revealing new groups of eukaryotic viruses, although disconnected from known hosts. Among these are the recently described mirusviruses, which share some similarities with herpesviruses.1 50 years ago, "herpes-type" viral particles2 were found in a thraustochytrid member of the labyrinthulomycetes, a diverse group of abundant and ecologically important marine eukaryotes,3,4 but could not be further characterized by methods then available. Long-read sequencing has allowed us to connect the biology of mirusviruses and thraustochytrids. We sequenced the genome of the genetically tractable model thraustochytrid Aurantiochytrium limacinum ATCC MYA-1381 and found that its 26 linear chromosomes have an extraordinary configuration. Subtelomeric ribosomal DNAs (rDNAs) found at all chromosome ends are interspersed with long repeated sequence elements denoted as long repeated-telomere and rDNA spacers (LORE-TEARS). We identified two genomic elements that are related to mirusvirus genomes. The first is a ∼300-kbp episome (circular element 1 [CE1]) present at a high copy number. Strikingly, the second, distinct, mirusvirus-like element is integrated between two sets of rDNAs and LORE-TEARS at the left end of chromosome 15 (LE-Chr15). Similar to metagenomically derived mirusviruses, these putative A. limacinum mirusviruses have a virion module related to that of herpesviruses along with an informational module related to nucleocytoplasmic large DNA viruses (NCLDVs). CE1 and LE-Chr15 bear striking similarities to episomal and endogenous latent forms of herpesviruses, respectively, and open new avenues of research into marine virus-host interactions.

Effects of the G-quadruplex-binding drugs quarfloxin and CX-5461 on the malaria parasite Plasmodium falciparum

Plasmodium falciparum is the deadliest causative agent of human malaria. This parasite has historically developed resistance to most drugs, including the current frontline treatments, so new therapeutic targets are needed. Our previous work on guanine quadruplexes (G4s) in the parasite's DNA and RNA has highlighted their influence on parasite biology, and revealed G4 stabilising compounds as promising candidates for repositioning. In particular, quarfloxin, a former anticancer agent, kills blood-stage parasites at all developmental stages, with fast rates of kill and nanomolar potency. Here we explored the molecular mechanism of quarfloxin and its related derivative CX-5461. In vitro, both compounds bound to P. falciparum-encoded G4 sequences. In cellulo, quarfloxin was more potent than CX-5461, and could prevent establishment of blood-stage malaria in vivo in a murine model. CX-5461 showed clear DNA damaging activity, as reported in human cells, while quarfloxin caused weaker signatures of DNA damage. Both compounds caused transcriptional dysregulation in the parasite, but the affected genes were largely different, again suggesting different modes of action. Therefore, CX-5461 may act primarily as a DNA damaging agent in both Plasmodium parasites and mammalian cells, whereas the complete antimalarial mode of action of quarfloxin may be parasite-specific and remains somewhat elusive.

Lessons from the deep: mechanisms behind diversification of eukaryotic protein complexes

Genetic variation is the major mechanism behind adaptation and evolutionary change. As most proteins operate through interactions with other proteins, changes in protein complex composition and subunit sequence provide potentially new functions. Comparative genomics can reveal expansions, losses and sequence divergence within protein-coding genes, but in silico analysis cannot detect subunit substitutions or replacements of entire protein complexes. Insights into these fundamental evolutionary processes require broad and extensive comparative analyses, from both in silico and experimental evidence. Here, we combine data from both approaches and consider the gamut of possible protein complex compositional changes that arise during evolution, citing examples of complete conservation to partial and total replacement by functional analogues. We focus in part on complexes in trypanosomes as they represent one of the better studied non-animal/non-fungal lineages, but extend insights across the eukaryotes by extensive comparative genomic analysis. We argue that gene loss plays an important role in diversification of protein complexes and hence enhancement of eukaryotic diversity.

A paradoxical population structure of var DBLα types in Africa

The var multigene family encodes the P. falciparum erythrocyte membrane protein 1 (PfEMP1), which is important in host-parasite interaction as a virulence factor and major surface antigen of the blood stages of the parasite, responsible for maintaining chronic infection. Whilst important in the biology of P. falciparum, these genes (50 to 60 genes per parasite genome) are routinely excluded from whole genome analyses due to their hyper-diversity, achieved primarily through recombination. The PfEMP1 head structure almost always consists of a DBLα-CIDR tandem. Categorised into different groups (upsA, upsB, upsC), different head structures have been associated with different ligand-binding affinities and disease severities. We study how conserved individual DBLα types are at the country, regional, and local scales in Sub-Saharan Africa. Using publicly-available sequence datasets and a novel ups classification algorithm, cUps, we performed an in silico exploration of DBLα conservation through time and space in Africa. In all three ups groups, the population structure of DBLα types in Africa consists of variants occurring at rare, low, moderate, and high frequencies. Non-rare variants were found to be temporally stable in a local area in endemic Ghana. When inspected across different geographical scales, we report different levels of conservation; while some DBLα types were consistently found in high frequencies in multiple African countries, others were conserved only locally, signifying local preservation of specific types. Underlying this population pattern is the composition of DBLα types within each isolate DBLα repertoire, revealed to also consist of a mix of types found at rare, low, moderate, and high frequencies in the population. We further discuss the adaptive forces and balancing selection, including host genetic factors, potentially shaping the evolution and diversity of DBLα types in Africa.

The constraints of allotopic expression

Allotopic expression is the functional transfer of an organellar gene to the nucleus, followed by synthesis of the gene product in the cytosol and import into the appropriate organellar sub compartment. Here, we focus on mitochondrial genes encoding OXPHOS subunits that were naturally transferred to the nucleus, and critically review experimental evidence that claim their allotopic expression. We emphasize aspects that may have been overlooked before, i.e., when modifying a mitochondrial gene for allotopic expression━besides adapting the codon usage and including sequences encoding mitochondrial targeting signals━three additional constraints should be considered: (i) the average apparent free energy of membrane insertion (μΔGapp) of the transmembrane stretches (TMS) in proteins earmarked for the inner mitochondrial membrane, (ii) the final, functional topology attained by each membrane-bound OXPHOS subunit; and (iii) the defined mechanism by which the protein translocator TIM23 sorts cytosol-synthesized precursors. The mechanistic constraints imposed by TIM23 dictate the operation of two pathways through which alpha-helices in TMS are sorted, that eventually determine the final topology of membrane proteins. We used the biological hydrophobicity scale to assign an average apparent free energy of membrane insertion (μΔGapp) and a "traffic light" color code to all TMS of OXPHOS membrane proteins, thereby predicting which are more likely to be internalized into mitochondria if allotopically produced. We propose that the design of proteins for allotopic expression must make allowance for μΔGapp maximization of highly hydrophobic TMS in polypeptides whose corresponding genes have not been transferred to the nucleus in some organisms.

RNA Polymerase II evolution and adaptations: Insights from Plasmodium and other parasitic protists

The C-terminal domain (CTD) of RNA polymerase II plays a crucial role in regulating transcription dynamics in eukaryotes. The phosphorylation of serine residues within the CTD controls transcription initiation, elongation, and termination. While the CTD is highly conserved across eukaryotes, lower eukaryotes like protists, including Plasmodium, exhibit some differences. In this study, we performed a comparative analysis of CTD in eukaryotic systems to understand why the parasites evolved in this particular manner. The Plasmodium falciparum RPB1 is exceptionally large and feature a gap between the first and second heptad repeats, resulting in fifteen canonical heptad repeats excluding the initial repeat. Analysis of this intervening sequence revealed sub motifs of heptads where two serine residues occupy the first and fourth positions (S1X2X3S4). These motifs lie in the intrinsically disordered region of RPB1, a characteristic feature of the CTD. Interestingly, the S1X2X3S4 sub-motif was also observed in early-divergingeukaryotes like Leishmania major, which lack canonical heptad repeats. Furthermore, eukaryotes across the phylogenetic tree revealed a sigmoid pattern of increasing serine frequency in the CTD, indicating that serine enrichment is a significant step in the evolution of heptad-rich RPB1. Based on these observations and analysis, we proposed an evolutionary model for RNA Polymerase II CTD, encompassing organisms previously deemed exceptions, notably Plasmodium species. Thus, our study provides novel insights into the evolution of the CTD and will prompt further investigations into the differences exhibited by Plasmodium RNA Pol II and determine if they confer a survival advantage to the parasite.

Invasion-inhibitory peptides chosen by natural selection analysis as an antimalarial strategy

Plasmodium vivax's biological complexity has restricted in vitro culture development for characterising antigens involved in erythrocyte invasion and their immunological relevance. The murine model is proposed as a suitable alternative in the search for therapeutic candidates since Plasmodium yoelii uses homologous proteins for its invasion. The AMA-1 protein is essential for parasite invasion of erythrocytes as it is considered an important target for infection control. This study has focused on functional PyAMA-1 peptides involved in host-pathogen interaction; the protein is located in regions under negative selection as determined by bioinformatics analysis. It was found that pyama1 has two highly conserved regions amongst species (>70%) under negative selection. Fourteen synthetic peptides spanning both conserved regions were evaluated; 5 PyAMA-1 peptides having high specific binding (HABP) to murine erythrocytes were identified. The parasite's invasion inhibition capability was analysed through in vitro assays, suggesting that peptides 42681 (43-ENTERSIKLINPWDKYMEKY-62), 42903 (206-RYSSNDANNENQPFSFTPEK-225) and 42904 (221-FTPEKIENYKDLSYLTKNLR-240) had greater than 50% inhibition profile and restricted P. yoelii intra-erythrocyte development. This work proposes that the screening of conserved HABPs under negative selective pressure might be good candidates for developing a synthetic anti-malarial vaccine since they share functionally-relevant characteristics, such as interspecies conservation, specific RBC binding profile, invasion and parasite development inhibition capability, and the predicted B-epitopes within were recognised by sera obtained from experimentally-infected mice.

Conceptual risk assessment of mosquito population modification gene-drive systems to control malaria transmission: preliminary hazards list workshops

The field-testing and eventual adoption of genetically-engineered mosquitoes (GEMs) to control vector-borne pathogen transmission will require them meeting safety criteria specified by regulatory authorities in regions where the technology is being considered for use and other locales that might be impacted. Preliminary risk considerations by researchers and developers may be useful for planning the baseline data collection and field research used to address the anticipated safety concerns. Part of this process is to identify potential hazards (defined as the inherent ability of an entity to cause harm) and their harms, and then chart the pathways to harm and evaluate their probability as part of a risk assessment. The University of California Malaria Initiative (UCMI) participated in a series of workshops held to identify potential hazards specific to mosquito population modification strains carrying gene-drive systems coupled to anti-parasite effector genes and their use in a hypothetical island field trial. The hazards identified were placed within the broader context of previous efforts discussed in the scientific literature. Five risk areas were considered i) pathogens, infections and diseases, and the impacts of GEMs on human and animal health, ii) invasiveness and persistence of GEMs, and interactions of GEMs with target organisms, iii) interactions of GEMs with non-target organisms including horizontal gene transfer, iv) impacts of techniques used for the management of GEMs and v) evolutionary and stability considerations. A preliminary hazards list (PHL) was developed and is made available here. This PHL is useful for internal project risk evaluation and is available to regulators at prospective field sites. UCMI project scientists affirm that the subsequent processes associated with the comprehensive risk assessment for the application of this technology should be driven by the stakeholders at the proposed field site and areas that could be affected by this intervention strategy.

Extracellular vesicles could be a putative posttranscriptional regulatory mechanism that shapes intracellular RNA levels in Plasmodium falciparum

Plasmodium falciparum secretes extracellular vesicles (PfEVs) that contain parasite-derived RNA. However, the significance of the secreted RNA remains unexplored. Here, we compare secreted and intracellular RNA from asexual cultures of six P. falciparum lines. We find that secretion of RNA via extracellular vesicles is not only periodic throughout the asexual intraerythrocytic developmental cycle but is also highly conserved across P. falciparum isolates. We further demonstrate that the phases of RNA secreted via extracellular vesicles are discernibly shifted compared to those of the intracellular RNA within the secreting whole parasite. Finally, transcripts of genes with no known function during the asexual intraerythrocytic developmental cycle are enriched in PfEVs compared to the whole parasite. We conclude that the secretion of extracellular vesicles could be a putative posttranscriptional RNA regulation mechanism that is part of or synergise the classic RNA decay processes to maintain intracellular RNA levels in P. falciparum.

A new Plasmodium vivax reference genome for South American isolates

BACKGROUND

Plasmodium vivax is the second most important cause of human malaria worldwide, and accounts for the majority of malaria cases in South America. A high-quality reference genome exists for Papua Indonesia (PvP01) and Thailand (PvW1), but is lacking for South America. A reference genome specifically for South America would be beneficial though, as P. vivax is a genetically diverse parasite with geographical clustering.

RESULTS

This study presents a new high-quality assembly of a South American P. vivax isolate, referred to as PvPAM (P. vivax Peruvian AMazon). The genome was obtained from a low input patient sample from the Peruvian Amazon and sequenced using PacBio technology, resulting in a highly complete assembly with 6497 functional genes. Telomeric ends were present in 17 out of 28 chromosomal ends, and additional (sub)telomeric regions are present in 12 unassigned contigs. A comparison of multigene families between PvPAM and the PvP01 genome revealed remarkable variation in vir genes, and the presence of merozoite surface proteins (MSP) 3.6 and 3.7. Three dhfr and dhps drug resistance associated mutations are present in PvPAM, similar to those found in other Peruvian isolates. Mapping of publicly available South American whole genome sequencing (WGS) data to PvPAM resulted in significantly fewer variants and truncated reads compared to the use of PvP01 or PvW1 as reference genomes. To minimize the number of core genome variants in non-South American samples, PvW1 is most suited for Southeast Asian isolates, both PvPAM and PvW1 are suited for South Asian isolates, and PvPAM is recommended for African isolates. Interestingly, non-South American samples still contained the least subtelomeric variants when mapped to PvPAM, indicating high quality of the PvPAM subtelomeric regions.

CONCLUSIONS

Our findings show that the PvPAM reference genome more accurately represents South American P. vivax isolates in comparison to PvP01 and PvW1. In addition, PvPAM has a high level of completeness, and contains a similar number of annotated genes as PvP01 or PvW1. The PvPAM genome therefore will be a valuable resource to improve future genomic analyses on P. vivax isolates from the South American continent.

Cohesin contributes to transcriptional repression of stage-specific genes in the human malaria parasite

The complex life cycle of the human malaria parasite, Plasmodium falciparum, is driven by specific transcriptional programs, but it is unclear how most genes are activated or silenced at specific times. There is an association between transcription and spatial organization; however, the molecular mechanisms behind genome organization are unclear. While P. falciparum lacks key genome-organizing proteins found in metazoans, it has all core components of the cohesin complex. To investigate the role of cohesin in P. falciparum, we functionally characterize the cohesin subunit Structural Maintenance of Chromosomes protein 3 (SMC3). SMC3 knockdown during early stages of the intraerythrocytic developmental cycle (IDC) upregulates a subset of genes involved in erythrocyte egress and invasion, which are normally expressed at later stages. ChIP-seq analyses reveal that during the IDC, SMC3 enrichment at the promoter regions of these genes inversely correlates with gene expression and chromatin accessibility. These data suggest that SMC3 binding contributes to the repression of specific genes until their appropriate time of expression, revealing a new mode of stage-specific gene repression in P. falciparum.

Understanding the structure and function of Plasmodium aminopeptidases to facilitate drug discovery

Malaria continues to be the most widespread parasitic disease affecting humans globally. As parasites develop drug resistance at an alarming pace, it has become crucial to identify novel drug targets. Over the last decade, the metalloaminopeptidases have gained importance as potential targets for new antimalarials. These enzymes are responsible for removing the N-terminal amino acids from proteins and peptides, and their restricted specificities suggest that many perform unique and essential roles within the malaria parasite. This mini-review focuses on the recent progress in structure and functional data relating to the Plasmodium metalloaminopeptidases that have been validated or shown promise as new antimalarial drug targets.

Detection and quantification of Plasmodium falciparum in human blood by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry: a proof of concept study

BACKGROUND

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF) has revolutionized identification of bacteria and is becoming available in an increasing number of laboratories in malaria-endemic countries. The purpose of this proof-of-concept study was to explore the potential of MALDI-TOF as a diagnostic tool for direct detection and quantification of Plasmodium falciparum in human blood.

METHODS

Three different P. falciparum strains (3D7, HB3 and IT4) were cultured and synchronized following standard protocols. Ring-stages were diluted in fresh blood group 0 blood drawn in EDTA from healthy subjects to mimic clinical samples. Samples were treated with saponin and washed in PBS to concentrate protein material. Samples were analysed using a Microflex LT MALDI-TOF and resulting mass spectra were compared using FlexAnalysis software.

RESULTS

More than 10 peaks specific for P. falciparum were identified. The identified peaks were consistent among the three genetically unrelated strains. Identification was possible in clinically relevant concentrations of 0.1% infected red blood cells, and a close relationship between peak intensity and the percentage of infected red blood cells was seen.

CONCLUSION

The findings indicate that the method has the potential to detect and quantify P. falciparum at clinically relevant infection intensities and provides proof-of-concept for MALDI-TOF-based diagnosis of human malaria. Further research is needed to include other Plasmodium spp., wildtype parasite isolates and to increase sensitivity. MALDI-TOF may be a useful tool for mass-screening purposes and for diagnosis of malaria in settings where it is readily available.