Mitotic evolution of Plasmodium falciparum shows a stable core genome but recombination in antigen families

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.

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.

Universal signatures of transposable element compartmentalization across eukaryotic genomes

The evolutionary mechanisms that drive the emergence of genome architecture remain poorly understood but can now be assessed with unprecedented power due to the massive accumulation of genome assemblies spanning phylogenetic diversity1,2. Transposable elements (TEs) are a rich source of large-effect mutations since they directly and indirectly drive genomic structural variation and changes in gene expression3. Here, we demonstrate universal patterns of TE compartmentalization across eukaryotic genomes spanning ~1.7 billion years of evolution, in which TEs colocalize with gene families under strong predicted selective pressure for dynamic evolution and involved in specific functions. For non-pathogenic species these genes represent families involved in defense, sensory perception and environmental interaction, whereas for pathogenic species, TE-compartmentalized genes are highly enriched for pathogenic functions. Many TE-compartmentalized gene families display signatures of positive selection at the molecular level. Furthermore, TE-compartmentalized genes exhibit an excess of high-frequency alleles for polymorphic TE insertions in fruit fly populations. We postulate that these patterns reflect selection for adaptive TE insertions as well as TE-associated structural variants. This process may drive the emergence of a shared TE-compartmentalized genome architecture across diverse eukaryotic lineages.

Strong positive selection biases identity-by-descent-based inferences of recent demography and population structure in Plasmodium falciparum

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD), yet strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we use simulations, a true IBD inference algorithm, and empirical data sets from different malaria transmission settings to investigate the extent of this bias and explore potential correction strategies. We analyze whole genome sequence data generated from 640 new and 3089 publicly available Plasmodium falciparum clinical isolates. We demonstrate that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discover that the removal of IBD peak regions partially restores the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness and extent of inbreeding. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Hyper-diverse antigenic variation and resilience to transmission-reducing intervention in falciparum malaria

Intervention against falciparum malaria in high transmission regions remains challenging, with relaxation of control efforts typically followed by rapid resurgence. Resilience to intervention co-occurs with incomplete immunity, whereby children eventually become protected from severe disease but not infection and a large transmission reservoir results from high asymptomatic prevalence across all ages. Incomplete immunity relates to the vast antigenic variation of the parasite, with the major surface antigen of the blood stage of infection encoded by the multigene family known as var. Recent deep sampling of var sequences from individual isolates in northern Ghana showed that parasite population structure exhibited persistent features of high-transmission regions despite the considerable decrease in prevalence during transient intervention with indoor residual spraying (IRS). We ask whether despite such apparent limited impact, the transmission system had been brought close to a transition in both prevalence and resurgence ability. With a stochastic agent-based model, we investigate the existence of such a transition to pre-elimination with intervention intensity, and of molecular indicators informative of its approach. We show that resurgence ability decreases sharply and nonlinearly across a narrow region of intervention intensities in model simulations, and identify informative molecular indicators based on var gene sequences. Their application to the survey data indicates that the transmission system in northern Ghana was brought close to transition by IRS. These results suggest that sustaining and intensifying intervention would have pushed malaria dynamics to a slow-rebound regime with an increased probability of local parasite extinction.

The problem of antimalarial resistance and its implications for drug discovery

INTRODUCTION

Malaria remains a devastating infectious disease with hundreds of thousands of casualties each year. Antimalarial drug resistance has been a threat to malaria control and elimination for many decades and is still of concern today. Despite the continued effectiveness of current first-line treatments, namely artemisinin-based combination therapies, the emergence of drug-resistant parasites in Southeast Asia and even more alarmingly the occurrence of resistance mutations in Africa is of great concern and requires immediate attention.

AREAS COVERED

A comprehensive overview of the mechanisms underlying the acquisition of drug resistance in Plasmodium falciparum is given. Understanding these processes provides valuable insights that can be harnessed for the development and selection of novel antimalarials with reduced resistance potential. Additionally, strategies to mitigate resistance to antimalarial compounds on the short term by using approved drugs are discussed.

EXPERT OPINION

While employing strategies that utilize already approved drugs may offer a prompt and cost-effective approach to counter antimalarial drug resistance, it is crucial to recognize that only continuous efforts into the development of novel antimalarial drugs can ensure the successful treatment of malaria in the future. Incorporating resistance propensity assessment during this developmental process will increase the likelihood of effective and enduring malaria treatments.

Direct long read visualization reveals metabolic interplay between two antimalarial drug targets

Increases in the copy number of large genomic regions, termed genome amplification, are an important adaptive strategy for malaria parasites. Numerous amplifications across the Plasmodium falciparum genome contribute directly to drug resistance or impact the fitness of this protozoan parasite. During the characterization of parasite lines with amplifications of the dihydroorotate dehydrogenase (DHODH) gene, we detected increased copies of an additional genomic region that encompassed 3 genes (~5 kb) including GTP cyclohydrolase I (GCH1 amplicon). While this gene is reported to increase the fitness of antifolate resistant parasites, GCH1 amplicons had not previously been implicated in any other antimalarial resistance context. Here, we further explored the association between GCH1 and DHODH copy number. Using long read sequencing and single read visualization, we directly observed a higher number of tandem GCH1 amplicons in parasites with increased DHODH copies (up to 9 amplicons) compared to parental parasites (3 amplicons). While all GCH1 amplicons shared a consistent structure, expansions arose in 2-unit steps (from 3 to 5 to 7, etc copies). Adaptive evolution of DHODH and GCH1 loci was further bolstered when we evaluated prior selection experiments; DHODH amplification was only successful in parasite lines with pre-existing GCH1 amplicons. These observations, combined with the direct connection between metabolic pathways that contain these enzymes, lead us to propose that the GCH1 locus is beneficial for the fitness of parasites exposed to DHODH inhibitors. This finding highlights the importance of studying variation within individual parasite genomes as well as biochemical connections of drug targets as novel antimalarials move towards clinical approval.

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.

Clinically relevant atovaquone-resistant human malaria parasites fail to transmit by mosquito

Long-acting injectable medications, such as atovaquone, offer the prospect of a "chemical vaccine" for malaria, combining drug efficacy with vaccine durability. However, selection and transmission of drug-resistant parasites is of concern. Laboratory studies have indicated that atovaquone resistance disadvantages parasites in mosquitoes, but lack of data on clinically relevant Plasmodium falciparum has hampered integration of these variable findings into drug development decisions. Here we generate atovaquone-resistant parasites that differ from wild type parent by only a Y268S mutation in cytochrome b, a modification associated with atovaquone treatment failure in humans. Relative to wild type, Y268S parasites evidence multiple defects, most marked in their development in mosquitoes, whether from Southeast Asia (Anopheles stephensi) or Africa (An. gambiae). Growth of asexual Y268S P. falciparum in human red cells is impaired, but parasite loss in the mosquito is progressive, from reduced gametocyte exflagellation, to smaller number and size of oocysts, and finally to absence of sporozoites. The Y268S mutant fails to transmit from mosquitoes to mice engrafted with human liver cells and erythrocytes. The severe-to-lethal fitness cost of clinically relevant atovaquone resistance to P. falciparum in the mosquito substantially lessens the likelihood of its transmission in the field.

Diverse evolutionary pathways challenge the use of collateral sensitivity as a strategy to suppress resistance

Drug resistance remains a major obstacle to malaria control and eradication efforts, necessitating the development of novel therapeutic strategies to treat this disease. Drug combinations based on collateral sensitivity, wherein resistance to one drug causes increased sensitivity to the partner drug, have been proposed as an evolutionary strategy to suppress the emergence of resistance in pathogen populations. In this study, we explore collateral sensitivity between compounds targeting the Plasmodium dihydroorotate dehydrogenase (DHODH). We profiled the cross-resistance and collateral sensitivity phenotypes of several DHODH mutant lines to a diverse panel of DHODH inhibitors. We focus on one compound, TCMDC-125334, which was active against all mutant lines tested, including the DHODH C276Y line, which arose in selections with the clinical candidate DSM265. In six selections with TCMDC-125334, the most common mechanism of resistance to this compound was copy number variation of the dhodh locus, although we did identify one mutation, DHODH I263S, which conferred resistance to TCMDC-125334 but not DSM265. We found that selection of the DHODH C276Y mutant with TCMDC-125334 yielded additional genetic changes in the dhodh locus. These double mutant parasites exhibited decreased sensitivity to TCMDC-125334 and were highly resistant to DSM265. Finally, we tested whether collateral sensitivity could be exploited to suppress the emergence of resistance in the context of combination treatment by exposing wildtype parasites to both DSM265 and TCMDC-125334 simultaneously. This selected for parasites with a DHODH V532A mutation which were cross-resistant to both compounds and were as fit as the wildtype parent in vitro. The emergence of these cross-resistant, evolutionarily fit parasites highlights the mutational flexibility of the DHODH enzyme.

Propensity of selecting mutant parasites for the antimalarial drug cabamiquine

We report an analysis of the propensity of the antimalarial agent cabamiquine, a Plasmodium-specific eukaryotic elongation factor 2 inhibitor, to select for resistant Plasmodium falciparum parasites. Through in vitro studies of laboratory strains and clinical isolates, a humanized mouse model, and volunteer infection studies, we identified resistance-associated mutations at 11 amino acid positions. Of these, six (55%) were present in more than one infection model, indicating translatability across models. Mathematical modelling suggested that resistant mutants were likely pre-existent at the time of drug exposure across studies. Here, we estimated a wide range of frequencies of resistant mutants across the different infection models, much of which can be attributed to stochastic differences resulting from experimental design choices. Structural modelling implicates binding of cabamiquine to a shallow mRNA binding site adjacent to two of the most frequently identified resistance mutations.

Uncovering the role of Rad51 in homologous recombination-mediated antigenic diversification in the human malaria parasite Plasmodium falciparum

The human malaria parasite Plasmodium falciparum maintains the chronicity of infections through antigenic variation, a well-coordinated immune evasion mechanism. The most prominent molecular determinant of antigenic variation in this parasite includes the members of the var multigene family. Homologous recombination (HR)-mediated genomic rearrangements have been implicated to play a major role in var gene diversification. However, the key molecular factors involved in the generation of diversity at var loci are less known. Here, we tested the hypothesis that PfRad51 could carry out recombination between var genes that are not homologous but homeologous in nature. We employed the whole-genome sequencing (WGS) approach to investigate recombination events among var sequences over 100 generations and compared the rate of sequence rearrangement at the var loci in both PfRad51-proficient and -deficient parasite lines. This brief report provides evidence that the loss of the key recombinase function renders the parasite with inefficient HR and results in fewer recombination events among the var sequences, thereby impacting the diversification of the var gene repertoire.

Decoding the impact of nuclear organization on antigenic variation in parasites

Antigenic variation as a strategy to evade the host adaptive immune response has evolved in divergent pathogens. Antigenic variation involves restricted, and often mutually exclusive, expression of dominant antigens and a periodic switch in antigen expression during infection. In eukaryotes, nuclear compartmentalization, including three-dimensional folding of the genome and physical separation of proteins in compartments or condensates, regulates mutually exclusive gene expression and chromosomal translocations. In this Review, we discuss the impact of nuclear organization on antigenic variation in the protozoan pathogens Trypanosoma brucei and Plasmodium falciparum. In particular, we highlight the relevance of nuclear organization in both mutually exclusive antigen expression and genome stability, which underlie antigenic variation.

Strong Positive Selection Biases Identity-By-Descent-Based Inferences of Recent Demography and Population Structure in Plasmodium falciparum

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD). Yet, strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we utilized simulations, a true IBD inference algorithm, and empirical datasets from different malaria transmission settings to investigate the extent of such bias and explore potential correction strategies. We analyzed whole genome sequence data generated from 640 new and 4,026 publicly available Plasmodium falciparum clinical isolates. Our findings demonstrated that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discovered that the removal of IBD peak regions partially restored the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Generation of a mutator parasite to drive resistome discovery in Plasmodium falciparum

In vitro evolution of drug resistance is a powerful approach for identifying antimalarial targets, however, key obstacles to eliciting resistance are the parasite inoculum size and mutation rate. Here we sought to increase parasite genetic diversity to potentiate resistance selections by editing catalytic residues of Plasmodium falciparum DNA polymerase δ. Mutation accumulation assays reveal a ~5-8 fold elevation in the mutation rate, with an increase of 13-28 fold in drug-pressured lines. Upon challenge with the spiroindolone PfATP4-inhibitor KAE609, high-level resistance is obtained more rapidly and at lower inocula than wild-type parasites. Selections also yield mutants with resistance to an "irresistible" compound, MMV665794 that failed to yield resistance with other strains. We validate mutations in a previously uncharacterised gene, PF3D7_1359900, which we term quinoxaline resistance protein (QRP1), as causal for resistance to MMV665794 and a panel of quinoxaline analogues. The increased genetic repertoire available to this "mutator" parasite can be leveraged to drive P. falciparum resistome discovery.

Analysis of genome instability and implications for the consequent phenotype in Plasmodium falciparum containing mutated MSH2-1 (P513T)

Malarial parasites exhibit extensive genomic plasticity, which induces the antigen diversification and the development of antimalarial drug resistance. Only a few studies have examined the genome maintenance mechanisms of parasites. The study aimed at elucidating the impact of a mutation in a DNA mismatch repair gene on genome stability by maintaining the mutant and wild-type parasites through serial in vitro cultures for approximately 400 days and analysing the subsequent spontaneous mutations. A P513T mutant of the DNA mismatch repair protein PfMSH2-1 from Plasmodium falciparum 3D7 was created. The mutation did not influence the base substitution rate but significantly increased the insertion/deletion (indel) mutation rate in short tandem repeats (STRs) and minisatellite loci. STR mutability was affected by allele size, genomic category and certain repeat motifs. In the mutants, significant telomere healing and homologous recombination at chromosomal ends caused extensive gene loss and generation of chimeric genes, resulting in large-scale chromosomal alteration. Additionally, the mutant showed increased tolerance to N-methyl-N'-nitro-N-nitrosoguanidine, suggesting that PfMSH2-1 was involved in recognizing DNA methylation damage. This work provides valuable insights into the role of PfMSH2-1 in genome stability and demonstrates that the genomic destabilization caused by its dysfunction may lead to antigen diversification.

Transmissibility of clinically relevant atovaquone-resistant Plasmodium falciparum by anopheline mosquitoes

Rising numbers of malaria cases and deaths underscore the need for new interventions. Long-acting injectable medications, such as those now in use for HIV prophylaxis, offer the prospect of a malaria "chemical vaccine", combining the efficacy of a drug (like atovaquone) with the durability of a biological vaccine. Of concern, however, is the possible selection and transmission of drug-resistant parasites. We addressed this question by generating clinically relevant, highly atovaquone-resistant, Plasmodium falciparum mutants competent to infect mosquitoes. Isogenic paired strains, that differ only by a single Y268S mutation in cytochrome b, were evaluated in parallel in southeast Asian (Anopheles stephensi) or African (Anopheles gambiae) mosquitoes, and thence in humanized mice. Fitness costs of the mutation were evident along the lifecycle, in asexual parasite growth in vitro and in a progressive loss of parasites in the mosquito. In numerous independent experiments, microscopic exam of salivary glands from hundreds of mosquitoes failed to detect even one Y268S sporozoite, a defect not rescued by coinfection with wild type parasites. Furthermore, despite uniformly successful transmission of wild type parasites from An. stephensi to FRG NOD huHep mice bearing human hepatocytes and erythrocytes, multiple attempts with Y268S-fed mosquitoes failed: there was no evidence of parasites in mouse tissues by microscopy, in vitro culture, or PCR. These studies confirm a severe-to-lethal fitness cost of clinically relevant atovaquone-resistant P. falciparum in the mosquito, and they significantly lessen the likelihood of their transmission in the field.

Very Low Rates of Spontaneous Gene Deletions and Gene Duplications in Dictyostelium discoideum

The study of spontaneous mutation rates has revealed a wide range of heritable point mutation rates across species, but there are comparatively few estimates for large-scale deletion and duplication rates. The handful of studies that have directly calculated spontaneous rates of deletion and duplication using mutation accumulation lines have estimated that genes are duplicated and deleted at orders of magnitude greater rates than the spontaneous point mutation rate. In our study, we tested whether spontaneous gene deletion and gene duplication rates are also high in Dictyostelium discoideum, a eukaryote with among the lowest point mutation rates (2.5 × 10-11 per site per generation) and an AT-rich genome (GC content of 22%). We calculated mutation rates of gene deletions and duplications using whole-genome sequencing data originating from a mutation accumulation experiment and determined the association between the copy number mutations and GC content. Overall, we estimated an average of 3.93 × 10-8 gene deletions and 1.18 × 10-8 gene duplications per gene per generation. While orders of magnitude greater than their point mutation rate, these rates are much lower compared to gene deletion and duplication rates estimated from mutation accumulation lines in other organisms (that are on the order of ~ 10-6 per gene/generation). The deletions and duplications were enriched in regions that were AT-rich even compared to the genomic background, in contrast to our expectations if low GC content was contributing to low mutation rates. The low deletion and duplication mutation rates in D. discoideum compared to other eukaryotes mirror their low point mutation rates, supporting previous work suggesting that this organism has high replication fidelity and effective molecular machinery to avoid the accumulation of mutations in their genome.

Unravelling var complexity: Relationship between DBLα types and var genes in Plasmodium falciparum

The enormous diversity and complexity of var genes that diversify rapidly by recombination has led to the exclusion of assembly of these genes from major genome initiatives (e.g., Pf6). A scalable solution in epidemiological surveillance of var genes is to use a small ‘tag’ region encoding the immunogenic DBLα domain as a marker to estimate var diversity. As var genes diversify by recombination, it is not clear the extent to which the same tag can appear in multiple var genes. This relationship between marker and gene has not been investigated in natural populations. Analyses of in vitro recombination within and between var genes have suggested that this relationship would not be exclusive. Using a dataset of publicly-available assembled var sequences, we test this hypothesis by studying DBLα-var relationships for four study sites in four countries: Pursat (Cambodia) and Mae Sot (Thailand), representing low malaria transmission, and Navrongo (Ghana) and Chikwawa (Malawi), representing high malaria transmission. In all study sites, DBLα-var relationships were shown to be predominantly 1-to-1, followed by a second largest proportion of 1-to-2 DBLα-var relationships. This finding indicates that DBLα tags can be used to estimate not just DBLα diversity but var gene diversity when applied in a local endemic area. Epidemiological applications of this result are discussed.

The mystery of persistent, asymptomatic Plasmodium falciparum infections

Plasmodium falciparum causes millions of malaria infections and hundreds of thousands of deaths annually. These parasites avoid the adaptive immune response by systematically cycling through a limited repertoire of variant surface antigens after which the number of circulating parasites drops to extremely low levels, coinciding with a loss of symptoms and eventual clearance of the infection. However, in regions with extended dry seasons or in individuals who no longer reside in endemic areas, asymptomatic infections have been observed to persist for many months or years, potentially serving as reservoirs for transmission. Recent work suggests the possibility that parasites can assume a state in which no variant surface antigens are expressed, thus rendering them virtually invisible to the immune system and enabling them to persist at low levels indefinitely.

Comparative overview of emerging RNA viruses: Epidemiology, pathogenesis, diagnosis and current treatment

From many decades, emerging infections have threatened humanity. The pandemics caused by different CoVs have already claimed and will continue to claim millions of lives. The SARS, Ebola, MERS epidemics and the most recent emergence of COVID-19 pandemic have threatened populations across borders. Since a highly pathogenic CoV has been evolved into the human population in the twenty-first century known as SARS, scientific advancements and innovative methods to tackle these viruses have increased in order to improve response preparedness towards the unpredictable threat posed by these rapidly emerging pathogens. Recently published review articles on SARS-CoV-2 have mainly focused on its pathogenesis, epidemiology and available treatments. However, in this review, we have done a systematic comparison of all three CoVs i.e., SARS, MERS and SARS-CoV-2 along with Ebola and Zika in terms of their epidemiology, virology, clinical features and current treatment strategies. This review focuses on important emerging RNA viruses starting from Zika, Ebola and the CoVs which include SARS, MERS and SARS-CoV-2. Each of these viruses has been elaborated on the basis of their epidemiology, virulence, transmission and treatment. However, special attention has been given to SARS-CoV-2 and the disease caused by it i.e., COVID-19 due to current havoc caused worldwide. At the end, insights into the current understanding of the lessons learned from previous epidemics to combat emerging CoVs have been described. The travel-related viral spread, the unprecedented nosocomial outbreaks and the high case-fatality rates associated with these highly transmissible and pathogenic viruses highlight the need for new prophylactic and therapeutic actions which include but are not limited to clinical indicators, contact tracing, and laboratory investigations as important factors that need to be taken into account in order to arrive at the final conclusion.

•

Recently published review articles on SARS-CoV-2 have mainly focused on its pathogenesis, epidemiology and available treatments.

•

The pandemics caused by different CoVs have already claimed and will continue to claim millions of lives.

•

This review focuses on important emerging RNA viruses starting from Zika, Ebola and the CoVs which include SARS, MERS and SARS-CoV-2.

•

Globally, numerous studies and researchers have recently started fighting this virus.

Using an antimalarial in mosquitoes overcomes Anopheles and Plasmodium resistance to malaria control strategies

The spread of insecticide resistance in Anopheles mosquitoes and drug resistance in Plasmodium parasites is contributing to a global resurgence of malaria, making the generation of control tools that can overcome these roadblocks an urgent public health priority. We recently showed that the transmission of Plasmodium falciparum parasites can be efficiently blocked when exposing Anopheles gambiae females to antimalarials deposited on a treated surface, with no negative consequences on major components of mosquito fitness. Here, we demonstrate this approach can overcome the hurdles of insecticide resistance in mosquitoes and drug resistant in parasites. We show that the transmission-blocking efficacy of mosquito-targeted antimalarials is maintained when field-derived, insecticide resistant Anopheles are exposed to the potent cytochrome b inhibitor atovaquone, demonstrating that this drug escapes insecticide resistance mechanisms that could potentially interfere with its function. Moreover, this approach prevents transmission of field-derived, artemisinin resistant P. falciparum parasites (Kelch13 C580Y mutant), proving that this strategy could be used to prevent the spread of parasite mutations that induce resistance to front-line antimalarials. Atovaquone is also highly effective at limiting parasite development when ingested by mosquitoes in sugar solutions, including in ongoing infections. These data support the use of mosquito-targeted antimalarials as a promising tool to complement and extend the efficacy of current malaria control interventions.

Effective control of malaria is hampered by resistance to vector-targeted insecticides and parasite-targeted drugs. This situation is exacerbated by a critical lack of chemical diversity in both interventions and, as such, new interventions are urgently needed. Recent laboratory studies have shown that an alternative approach based on treating Anopheles mosquitoes directly with antimalarial compounds can make mosquitoes incapable of transmitting the Plasmodium parasites that cause malaria. While promising, showing that mosquito-targeted antimalarials remain effective against wild parasites and mosquitoes, including drug- and insecticide-resistant populations in malaria-endemic countries, is crucial to the future viability of this approach. In this study, carried out in the US and Burkina Faso, we show that insecticide-resistance mechanisms found in highly resistant, natural Anopheles mosquito populations do not interfere with the transmission blocking activity of tarsal exposure to the antimalarial atovaquone, and that mosquito-targeted antimalarial exposure can block transmission of parasites resistant to the main therapeutic antimalarial drug artemisinin. By combining lab, and field-based studies in this way we have demonstrated that this novel approach can be effective in areas where conventional control measures are no longer as effective.

Stochastic expression of invasion genes in Plasmodium falciparum schizonts

Genetically identical cells are known to exhibit differential phenotypes in the same environmental conditions. These phenotypic variants are linked to transcriptional stochasticity and have been shown to contribute towards adaptive flexibility of a wide range of unicellular organisms. Here, we investigate transcriptional heterogeneity and stochastic gene expression in Plasmodium falciparum by performing the quasilinear multiple annealing and looping based amplification cycles (MALBAC) based amplification and single cell RNA sequencing of blood stage schizonts. Our data reveals significant transcriptional variations in the schizont stage with a distinct group of highly variable invasion gene transcripts being identified. Moreover, the data reflects several diversification processes including putative developmental “checkpoint”; transcriptomically distinct parasite sub-populations and transcriptional switches in variable gene families (var, rifin, phist). Most of these features of transcriptional variability are preserved in isogenic parasite cell populations (albeit with a lesser amplitude) suggesting a role of epigenetic factors in cell-to-cell transcriptional variations in human malaria parasites. Lastly, we apply quantitative RT-PCR and RNA-FISH approach and confirm stochastic expression of key invasion genes, such as, msp1, msp3, msp7, eba181 and ama1 which represent prime candidates for invasion-blocking vaccines.

Genetically identical cells can be phenotypically diverse to allow adaptive flexibility in a given environment. This phenotypic diversity is driven by epigenetic and transcriptional variability. Here, Tripathi et al. perform scRNA-seq of isogenic and non-isogenic Plasmodium falciparum schizont populations to explore transcriptional heterogeneity and stochastic gene expression during the course of development.

Current challenges and nanotechnology-based pharmaceutical strategies for the treatment and control of malaria

Malaria is one of the prevalent tropical diseases caused by the parasitic protozoan of the genus Plasmodium spp. With an estimated 228 million cases, it is a major public health concern with high incidence of morbidity and mortality worldwide. The emergence of drug-resistant parasites, inadequate vector control measures, and the non-availability of effective vaccine(s) against malaria pose a serious challenge to malaria eradication especially in underdeveloped and developing countries. Malaria treatment and control comprehensively relies on chemical compounds, which encompass various complications, including severe toxic effects, emergence of drug resistance, and high cost of therapy. To overcome the clinical failures of anti-malarial chemotherapy, a new drug development is of an immediate need. However, the drug discovery and development process is expensive and time consuming. In such a scenario, nanotechnological strategies may offer promising alternative approach for the treatment and control of malaria, with improved efficacy and safety. Nanotechnology based formulations of existing anti-malarial chemotherapeutic agents prove to exceed the limitations of existing therapies in relation to optimum therapeutic benefits, safety, and cost effectiveness, which indeed advances the patient's compliance in treatment. In this review, the shortcomings of malaria therapeutics and necessity of nanotechnological strategies for treating malaria were discussed.

Shared Mechanisms for Mutually Exclusive Expression and Antigenic Variation by Protozoan Parasites

Cellular decision-making at the level of gene expression is a key process in the development and evolution of every organism. Variations in gene expression can lead to phenotypic diversity and the development of subpopulations with adaptive advantages. A prime example is the mutually exclusive activation of a single gene from within a multicopy gene family. In mammals, this ranges from the activation of one of the two immunoglobulin (Ig) alleles to the choice in olfactory sensory neurons of a single odorant receptor (OR) gene from a family of more than 1,000. Similarly, in parasites like Trypanosoma brucei, Giardia lamblia or Plasmodium falciparum, the process of antigenic variation required to escape recognition by the host immune system involves the monoallelic expression of vsg, vsp or var genes, respectively. Despite the importance of this process, understanding how this choice is made remains an enigma. The development of powerful techniques such as single cell RNA-seq and Hi-C has provided new insights into the mechanisms these different systems employ to achieve monoallelic gene expression. Studies utilizing these techniques have shown how the complex interplay between nuclear architecture, physical interactions between chromosomes and different chromatin states lead to single allele expression. Additionally, in several instances it has been observed that high-level expression of a single gene is preceded by a transient state where multiple genes are expressed at a low level. In this review, we will describe and compare the different strategies that organisms have evolved to choose one gene from within a large family and how parasites employ this strategy to ensure survival within their hosts.

The Extracellular Milieu of Toxoplasma's Lytic Cycle Drives Lab Adaptation, Primarily by Transcriptional Reprogramming

Evolve and resequencing (E&R) was applied to lab adaptation of Toxoplasma gondii for over 1,500 generations with the goal of mapping host-independent in vitro virulence traits. Phenotypic assessments of steps across the lytic cycle revealed that only traits needed in the extracellular milieu evolved. Nonsynonymous single-nucleotide polymorphisms (SNPs) in only one gene, a P4 flippase, fixated across two different evolving populations, whereas dramatic changes in the transcriptional signature of extracellular parasites were identified. Newly developed computational tools correlated phenotypes evolving at different rates with specific transcriptomic changes. A set of 300 phenotype-associated genes was mapped, of which nearly 50% is annotated as hypothetical. Validation of a select number of genes by knockouts confirmed their role in lab adaptation and highlights novel mechanisms underlying in vitro virulence traits. Further analyses of differentially expressed genes revealed the development of a “pro-tachyzoite” profile as well as the upregulation of the fatty acid biosynthesis (FASII) pathway. The latter aligned with the P4 flippase SNP and aligned with a low abundance of medium-chain fatty acids at low passage, indicating this is a limiting factor in extracellular parasites. In addition, partial overlap with the bradyzoite differentiation transcriptome in extracellular parasites indicated that stress pathways are involved in both situations. This was reflected in the partial overlap between the assembled ApiAP2 and Myb transcription factor network underlying the adapting extracellular state with the bradyzoite differentiation program. Overall, E&R is a new genomic tool successfully applied to map the development of polygenic traits underlying in vitro virulence of T. gondii.

IMPORTANCE It has been well established that prolonged in vitro cultivation of Toxoplasma gondii augments progression of the lytic cycle. This lab adaptation results in increased capacities to divide, migrate, and survive outside a host cell, all of which are considered host-independent virulence factors. However, the mechanistic basis underlying these enhanced virulence features is unknown. Here, E&R was utilized to empirically characterize the phenotypic, genomic, and transcriptomic changes in the non-lab-adapted strain, GT1, during 2.5 years of lab adaptation. This identified the shutdown of stage differentiation and upregulation of lipid biosynthetic pathways as the key processes being modulated. Furthermore, lab adaptation was primarily driven by transcriptional reprogramming, which rejected the starting hypothesis that genetic mutations would drive lab adaptation. Overall, the work empirically shows that lab adaptation augments T. gondii’s in vitro virulence by transcriptional reprogramming and that E&R is a powerful new tool to map multigenic traits.

Plasmodium falciparum Multidrug Resistance Proteins (pfMRPs)

The capacity of the lethal Plasmodium falciparum parasite to develop resistance against anti-malarial drugs represents a central challenge in the global control and elimination of malaria. Historically, the action of drug transporters is known to play a pivotal role in the capacity of the parasite to evade drug action. MRPs (Multidrug Resistance Protein) are known in many phylogenetically diverse groups to be related to drug resistance by being able to handle a large range of substrates, including important endogenous substances as glutathione and its conjugates. P. falciparum MRPs are associated with in vivo and in vitro altered drug response, and might be important factors for the development of multi-drug resistance phenotypes, a latent possibility in the present, and future, combination therapy environment. Information on P. falciparum MRPs is scattered in the literature, with no specialized review available. We herein address this issue by reviewing the present state of knowledge.

Single-genome sequencing reveals within-host evolution of human malaria parasites

Population genomics of bulk malaria infections is unable to examine intrahost evolution; therefore, most work has focused on the role of recombination in generating genetic variation. We used single-cell sequencing protocol for low-parasitaemia infections to generate 406 near-complete single Plasmodium vivax genomes from 11 patients sampled during sequential febrile episodes. Parasite genomes contain hundreds of de novo mutations, showing strong signatures of selection, which are enriched in the ApiAP2 family of transcription factors, known targets of adaptation. Comparing 315 P. falciparum single-cell genomes from 15 patients with our P. vivax data, we find broad complementary patterns of de novo mutation at the gene and pathway level, revealing the importance of within-host evolution during malaria infections.

Single-cell genome sequencing of protozoan parasites

Despite considerable genetic variation within hosts, most parasite genome sequencing studies focus on bulk samples composed of millions of cells. Analysis of bulk samples is biased toward the dominant genotype, concealing cell-to-cell variation and rare variants. To tackle this, single-cell sequencing approaches have been developed and tailored to specific host-parasite systems. These are allowing the genetic diversity and kinship in complex parasite populations to be deciphered and for de novo genetic variation to be captured. Here, we outline the methodologies being used for single-cell sequencing of parasitic protozoans, such as Plasmodium and Leishmania spp., and how these tools are being applied to understand parasite biology.

Ape Origins of Human Malaria

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

The antimalarial MMV688533 provides potential for single-dose cures with a high barrier to Plasmodium falciparum parasite resistance

The emergence and spread of Plasmodium falciparum resistance to first-line antimalarials creates an imperative to identify and develop potent preclinical candidates with distinct modes of action. Here, we report the identification of MMV688533, an acylguanidine that was developed following a whole-cell screen with compounds known to hit high-value targets in human cells. MMV688533 displays fast parasite clearance in vitro and is not cross-resistant with known antimalarials. In a P. falciparum NSG mouse model, MMV688533 displays a long-lasting pharmacokinetic profile and excellent safety. Selection studies reveal a low propensity for resistance, with modest loss of potency mediated by point mutations in PfACG1 and PfEHD. These proteins are implicated in intracellular trafficking, lipid utilization, and endocytosis, suggesting interference with these pathways as a potential mode of action. This preclinical candidate may offer the potential for a single low-dose cure for malaria.

Evolution of transcriptional control of antigenic variation and virulence in human and ape malaria parasites

Background

The most severe form of human malaria is caused by the protozoan parasite Plasmodium falciparum. This unicellular organism is a member of a subgenus of Plasmodium called the Laverania that infects apes, with P. falciparum being the only member that infects humans. The exceptional virulence of this species to humans can be largely attributed to a family of variant surface antigens placed by the parasites onto the surface of infected red blood cells that mediate adherence to the vascular endothelium. These proteins are encoded by a large, multicopy gene family called var, with each var gene encoding a different form of the protein. By changing which var gene is expressed, parasites avoid immune recognition, a process called antigenic variation that underlies the chronic nature of malaria infections.

Results

Here we show that the common ancestor of the branch of the Laverania lineage that includes the human parasite underwent a remarkable change in the organization and structure of elements linked to the complex transcriptional regulation displayed by the var gene family. Unlike the other members of the Laverania, the clade that gave rise to P. falciparum evolved distinct subsets of var genes distinguishable by different upstream transcriptional regulatory regions that have been associated with different expression profiles and virulence properties. In addition, two uniquely conserved var genes that have been proposed to play a role in coordinating transcriptional switching similarly arose uniquely within this clade. We hypothesize that these changes originated at a time of dramatic climatic change on the African continent that is predicted to have led to significant changes in transmission dynamics, thus selecting for patterns of antigenic variation that enabled lengthier, more chronic infections.

Conclusions

These observations suggest that changes in transmission dynamics selected for significant alterations in the transcriptional regulatory mechanisms that mediate antigenic variation in the parasite lineage that includes P. falciparum. These changes likely underlie the chronic nature of these infections as well as their exceptional virulence.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01872-z.

Identification of a novel variant erythrocyte surface antigen-1 (VESA1) in Babesia orientalis

Babesia orientalis, belonging to the phylum Apicomplexa, is mainly accountable for water buffalo babesiosis, which adversely affected the livestock industry in China. Variant erythrocyte surface antigen-1 (VESA1), an antigen that helps infected erythrocytes to escape from host immune responses, was first reported in Babesia bovis. Various VESA1 proteins have also been characterized in other Babesia species. Nevertheless, there is no research on the identification and characterization of VESA1 proteins in Babesia orientalis. In this study, the BoVESA1 gene was amplified from both gDNA and cDNA. The results revealed that it is an intronless gene with a full length of 753 bp, encoding a protein of 250 amino acids with a predicted molecular weight of 28 kDa. The coding sequence (CDS) was cloned into the pGEX-6p-1 vector using a homologous recombination kit and expressed as a glutathione-S-transferase (GST)-fusion protein with a molecular weight of 53 kDa. The tertiary structure of BoVESA1 was predicted using the I-TASSER software. The recombinant protein was subjected to western blotting; the immunogenicity of recombinant BoVESA1 (rBoVESA1) was identified by incubating it with B. orientalis-positive serum. The native BoVESA1 was identified using the lysates of B. orientalis-infected water buffalo erythrocytes incubated with the anti-rBoVESA1 mouse serum. The results showed a band of ~ 28 kDa, which is similar to the predicted size. Immunofluorescence assay (IFA) using anti-rBoVESA1 serum probed indicated a strong signal in the infected RBCs, while the negative control showed no signal. In conclusion, the VESA1 protein was first identified in B. orientalis. This study facilitated further investigation of B. orientalis, and the results indicated that BoVESA1 may serve as a potential candidate antigen for diagnosis and detection of B. orientalis infection.

MalDA, Accelerating Malaria Drug Discovery

The Malaria Drug Accelerator (MalDA) is a consortium of 15 leading scientific laboratories. The aim of MalDA is to improve and accelerate the early antimalarial drug discovery process by identifying new, essential, druggable targets. In addition, it seeks to produce early lead inhibitors that may be advanced into drug candidates suitable for preclinical development and subsequent clinical testing in humans. By sharing resources, including expertise, knowledge, materials, and reagents, the consortium strives to eliminate the structural barriers often encountered in the drug discovery process. Here we discuss the mission of the consortium and its scientific achievements, including the identification of new chemically and biologically validated targets, as well as future scientific directions.

Assessing risks of Plasmodium falciparum resistance to select next-generation antimalarials

Strategies to counteract or prevent emerging drug resistance are crucial for the design of next-generation antimalarials. In the past, resistant parasites were generally identified following treatment failures in patients, and compounds would have to be abandoned late in development. An early understanding of how candidate therapeutics lose efficacy as parasites evolve resistance is important to facilitate drug design and improve resistance detection and monitoring up to the postregistration phase. We describe a new strategy to assess resistance to antimalarial compounds as early as possible in preclinical development by leveraging tools to define the Plasmodium falciparum resistome, predict potential resistance risks of clinical failure for candidate therapeutics, and inform decisions to guide antimalarial drug development.

How Many Is Enough? - Challenges of Multinucleated Cell Division in Malaria Parasites

Regulating the number of progeny generated by replicative cell cycles is critical for any organism to best adapt to its environment. Classically, the decision whether to divide further is made after cell division is completed by cytokinesis and can be triggered by intrinsic or extrinsic factors. Contrarily, cell cycles of some species, such as the malaria-causing parasites, go through multinucleated cell stages. Hence, their number of progeny is determined prior to the completion of cell division. This should fundamentally affect how the process is regulated and raises questions about advantages and challenges of multinucleation in eukaryotes. Throughout their life cycle Plasmodium spp. parasites undergo four phases of extensive proliferation, which differ over three orders of magnitude in the amount of daughter cells that are produced by a single progenitor. Even during the asexual blood stage proliferation parasites can produce very variable numbers of progeny within one replicative cycle. Here, we review the few factors that have been shown to affect those numbers. We further provide a comparative quantification of merozoite numbers in several P. knowlesi and P. falciparum parasite strains, and we discuss the general processes that may regulate progeny number in the context of host-parasite interactions. Finally, we provide a perspective of the critical knowledge gaps hindering our understanding of the molecular mechanisms underlying this exciting and atypical mode of parasite multiplication.

Single-cell sequencing of the small and AT-skewed genome of malaria parasites

Single-cell genomics is a rapidly advancing field; however, most techniques are designed for mammalian cells. We present a single-cell sequencing pipeline for an intracellular parasite, Plasmodium falciparum, with a small genome of extreme base content. Through optimization of a quasi-linear amplification method, we target the parasite genome over contaminants and generate coverage levels allowing detection of minor genetic variants. This work, as well as efforts that build on these findings, will enable detection of parasite heterogeneity contributing to P. falciparum adaptation. Furthermore, this study provides a framework for optimizing single-cell amplification and variant analysis in challenging genomes.

K13-Mediated Reduced Susceptibility to Artemisinin in Plasmodium falciparum Is Overlaid on a Trait of Enhanced DNA Damage Repair

Southeast Asia has been the hotbed for the development of drug-resistant malaria parasites, including those with resistance to artemisinin combination therapy. While mutations in the kelch propeller domain (K13 mutations) are associated with artemisinin resistance, a range of evidence suggests that other factors are critical for the establishment and subsequent transmission of resistance in the field. Here, we perform a quantitative analysis of DNA damage and repair in the malaria parasite Plasmodium falciparum and find a strong link between enhanced DNA damage repair and artemisinin resistance. This experimental observation is further supported when variations in seven known DNA repair genes are found in resistant parasites, with six of these mutations being associated with K13 mutations. Our data provide important insights on confounding factors that are important for the establishment and spread of artemisinin resistance and may explain why resistance has not yet arisen in Africa.

•

High-throughput MalariaCometChip to measure DNA damage level in P. falciparum

•

Subpopulation of Cambodian isolates possess enhanced DNA damage repair

•

Important link between enhanced DNA damage repair and artemisinin resistance

The ApiAP2 factor PfAP2-HC is an integral component of heterochromatin in the malaria parasite Plasmodium falciparum

Malaria parasites undergo a complex life cycle in the human host and the mosquito vector. The ApiAP2 family of DNA-binding proteins plays a dominant role in parasite development and life cycle progression. Most ApiAP2 factors studied to date act as transcription factors regulating stage-specific gene expression. Here, we characterized an ApiAP2 factor in Plasmodium falciparum that we termed PfAP2-HC. We demonstrate that PfAP2-HC specifically binds to heterochromatin throughout the genome. Intriguingly, PfAP2-HC does not bind DNA in vivo and recruitment of PfAP2-HC to heterochromatin is independent of its DNA-binding domain but strictly dependent on heterochromatin protein 1. Furthermore, our results suggest that PfAP2-HC functions neither in the regulation of gene expression nor in heterochromatin formation or maintenance. In summary, our findings reveal PfAP2-HC as a core component of heterochromatin in malaria parasites and identify unexpected properties and substantial functional divergence among the members of the ApiAP2 family of regulatory proteins.

•

The ApiAP2 factor AP2-HC is a core component of heterochromatin in malaria parasites

•

Binding of AP2-HC to heterochromatin strictly depends on heterochromatin protein 1

•

The AP2 DNA-binding domain of AP2-HC is dispensable for heterochromatin association

Cell wall protein variation, break-induced replication, and subtelomere dynamics in Candida glabrata

Candida glabrata is an opportunistic pathogen of humans, responsible for up to 30% of disseminated candidiasis. Adherence of C. glabrata to host cells is mediated by adhesin-like proteins (ALPs), about half of which are encoded in the subtelomeres. We performed a de novo assembly of two C. glabrata strains, BG2 and BG3993, using long single-molecule real-time (SMRT) reads, and constructed high-quality telomere-to-telomere assemblies of all 13 chromosomes to assess differences between C. glabrata strains. We documented variation between strains, and in agreement with earlier studies, found high (~0.5%-1%) frequencies of SNVs across the genome, including within subtelomeric regions. We documented changes in ALP gene structure and complement: there are large length differences in ALP genes in different strains, resulting from copy number variation in tandem repeats. We compared strains to characterize chromosome rearrangement events including within the poorly characterized subtelomeric regions. We show that rearrangements within the subtelomere regions all affect ALP-encoding genes, and 14/16 involve just the most terminal ALP gene. We present evidence that these rearrangements are mediated by break-induced replication. This study highlights the constrained nature of subtelomeric changes impacting ALP gene complement and subtelomere structure.

Genomic analysis reveals independent evolution of Plasmodium falciparum populations in Ethiopia

Background

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

Methods

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

Results

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

Conclusion

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

Supplementary Information

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

Evolutionary analyses of the major variant surface antigen-encoding genes reveal population structure of Plasmodium falciparum within and between continents

Malaria remains a major public health problem in many countries. Unlike influenza and HIV, where diversity in immunodominant surface antigens is understood geographically to inform disease surveillance, relatively little is known about the global population structure of PfEMP1, the major variant surface antigen of the malaria parasite Plasmodium falciparum. The complexity of the var multigene family that encodes PfEMP1 and that diversifies by recombination, has so far precluded its use in malaria surveillance. Recent studies have demonstrated that cost-effective deep sequencing of the region of var genes encoding the PfEMP1 DBLα domain and subsequent classification of within host sequences at 96% identity to define unique DBLα types, can reveal structure and strain dynamics within countries. However, to date there has not been a comprehensive comparison of these DBLα types between countries. By leveraging a bioinformatic approach (jumping hidden Markov model) designed specifically for the analysis of recombination within var genes and applying it to a dataset of DBLα types from 10 countries, we are able to describe population structure of DBLα types at the global scale. The sensitivity of the approach allows for the comparison of the global dataset to ape samples of Plasmodium Laverania species. Our analyses show that the evolution of the parasite population emerging out of Africa underlies current patterns of DBLα type diversity. Most importantly, we can distinguish geographic population structure within Africa between Gabon and Ghana in West Africa and Uganda in East Africa. Our evolutionary findings have translational implications in the context of globalization. Firstly, DBLα type diversity can provide a simple diagnostic framework for geographic surveillance of the rapidly evolving transmission dynamics of P. falciparum. It can also inform efforts to understand the presence or absence of global, regional and local population immunity to major surface antigen variants. Additionally, we identify a number of highly conserved DBLα types that are present globally that may be of biological significance and warrant further characterization.

Telomere length dynamics in response to DNA damage in malaria parasites

Malaria remains a major cause of morbidity and mortality in the developing world. Recent work has implicated chromosome end stability and the repair of DNA breaks through telomere healing as potent drivers of variant antigen diversification, thus associating basic mechanisms for maintaining genome integrity with aspects of host-parasite interactions. Here we applied long-read sequencing technology to precisely examine the dynamics of telomere addition and chromosome end stabilization in response to double-strand breaks within subtelomeric regions. We observed that the process of telomere healing induces the initial synthesis of telomere repeats well in excess of the minimal number required for end stability. However, once stabilized, these newly created telomeres appear to function normally, eventually returning to a length nearing that of intact chromosome ends. These results parallel recent observations in humans, suggesting an evolutionarily conserved mechanism for chromosome end repair.

Multistage and transmission-blocking targeted antimalarials discovered from the open-source MMV Pandemic Response Box

Chemical matter is needed to target the divergent biology associated with the different life cycle stages of Plasmodium. Here, we report the parallel de novo screening of the Medicines for Malaria Venture (MMV) Pandemic Response Box against Plasmodium asexual and liver stage parasites, stage IV/V gametocytes, gametes, oocysts and as endectocides. Unique chemotypes were identified with both multistage activity or stage-specific activity, including structurally diverse gametocyte-targeted compounds with potent transmission-blocking activity, such as the JmjC inhibitor ML324 and the antitubercular clinical candidate SQ109. Mechanistic investigations prove that ML324 prevents histone demethylation, resulting in aberrant gene expression and death in gametocytes. Moreover, the selection of parasites resistant to SQ109 implicates the druggable V-type H⁺-ATPase for the reduced sensitivity. Our data therefore provides an expansive dataset of compounds that could be redirected for antimalarial development and also point towards proteins that can be targeted in multiple parasite life cycle stages.

Here, Reader et al. screen the Medicines for Malaria Venture Pandemic Response Box in parallel against Plasmodiumasexual and liver stage parasites, stage IV/V gametocytes, gametes, oocysts and as endectocides. They identify two potent transmission-blocking drugs: a histone demethylase inhibitor ML324 and the antitubercular SQ109.

Proteome size reduction in Apicomplexans is linked with loss of DNA repair and host redundant pathways

Apicomplexans are alveolate parasites which include Plasmodium falciparum, the main cause of malaria, one of the world's biggest killers from infectious disease. Apicomplexans are characterized by a reduction in proteome size, which appears to result from metabolic and functional simplification, commensurate with their parasitic lifestyle. However, other factors may also help to explain gene loss such as population bottlenecks experienced during transmission, and the effect of reducing the overall genomic information content. The latter constitutes an 'informational constraint', which is proposed to exert a selective pressure to evolve and maintain genes involved in informational fidelity and error correction, proportional to the quantity of information in the genome (which approximates to proteome size). The dynamics of gene loss was examined in 41 Apicomplexan genomes using orthogroup analysis. We show that loss of genes involved in amino acid metabolism and steroid biosynthesis can be explained by metabolic redundancy with the host. We also show that there is a marked tendency to lose DNA repair genes as proteome size is reduced. This may be explained by a reduction in size of the informational constraint and can help to explain elevated mutation rates in pathogens with reduced genome size. Multiple Sequentially Markovian Coalescent (MSMC) analysis indicates a recent bottleneck, consistent with predictions generated using allele-based population genetics approaches, implying that relaxed selection pressure due to reduced population size might have contributed to gene loss. However, the non-randomness of pathways that are lost challenges this scenario. Lastly, we identify unique orthogroups in malaria-causing Plasmodium species that infect humans, with a high proportion of membrane associated proteins. Thus, orthogroup analysis appears useful for identifying novel candidate pathogenic factors in parasites, when there is a wide sample of genomes available.

Plasmodium falciparum multidrug resistance gene-1 polymorphisms in Northern Nigeria: implications for the continued use of artemether-lumefantrine in the region

Background

The analysis of single nucleotide polymorphism (SNPs) in drug-resistance associated genes is a commonly used strategy for the surveillance of anti-malarial drug resistance in populations of parasites. The present study was designed and performed to provide genetic epidemiological data of the prevalence of N86Y-Y184F-D1246Y SNPs in Plasmodium falciparum multidrug resistance 1 (pfmdr1) in the malaria hotspot of Northern Nigeria.

Methods

Plasmodium falciparum-positive blood samples on Whatman-3MM filter papers were collected from 750 symptomatic patients from four states (Kano, Kaduna, Yobe and Adamawa) in Northern Nigeria, and genotyped via BigDye (v3.1) terminator cycle sequencing for the presence of three SNPs in pfmdr1. SNPs in pfmdr1 were used to construct NYD, NYY, NFY, NFD, YYY, YYD, YFD and YFY haplotypes, and all data were analysed using Pearson Chi square and Fisher’s exact (FE) tests.

Results

The prevalence of the pfmdr1 86Y allele was highest in Kaduna (12.50%, ² = 10.50, P = 0.02), whilst the 184F allele was highest in Kano (73.10%, ² = 13.20, P = 0.00), and the pfmdr1 1246Y allele was highest in Yobe (5.26%, ² = 9.20, P = 0.03). The NFD haplotype had the highest prevalence of 69.81% in Kano (² = 36.10, P = 0.00), followed by NYD with a prevalence of 49.00% in Adamawa, then YFD with prevalence of 11.46% in Kaduna. The YYY haplotype was not observed in any of the studied states.

Conclusion

The present study suggests that strains of P. falciparum with reduced sensitivity to the lumefantrine component of AL exist in Northern Nigeria and predominate in the North-West region.

The contribution of extrachromosomal DNA to genome plasticity in malaria parasites

Malaria caused by the protozoan parasite Plasmodium falciparum continues to impose significant morbidity and mortality, despite substantial investment into drug and vaccine development and deployment. Underlying the resilience of this parasite is its remarkable ability to undergo genome modifications, thus, providing parasite populations with extensive genetic variability that accelerates selection of drug resistance and limits the efficacy of most vaccines. This genome plasticity is rooted in the mechanisms of DNA repair that parasites employ to maintain genome integrity, a process skewed toward homologous recombination through the evolutionary loss of classical nonhomologous end joining. Repair of DNA double-strand breaks have been shown to enable "shuffling" of antigen-encoding gene sequences to vastly increase antigen diversity and to enable copy number expansion of genes that contribute to drug resistance. The latter phenomenon has been proposed to be a major contributor to the rise of resistance to several classes of antimalarial drugs. In this issue of Molecular Microbiology, McDaniels and colleagues add yet another mechanism that malaria parasites use to reduce drug susceptibility by demonstrating that P. falciparum can maintain expanded arrays of drug resistance cassettes as stably replicating, circular, extrachromosomal DNAs, thus, expanding genome plasticity beyond the parasite's 14 nuclear chromosomes.

Torin 2 Derivative, NCATS-SM3710, Has Potent Multistage Antimalarial Activity through Inhibition of P. falciparum Phosphatidylinositol 4-Kinase (Pf PI4KIIIβ)

Drug resistance is a constant threat to malaria control efforts making it important to maintain a good pipeline of new drug candidates. Of particular need are compounds that also block transmission by targeting sexual stage parasites. Mature sexual stages are relatively resistant to all currently used antimalarials except the 8-aminoquinolines that are not commonly used due to potential side effects. Here, we synthesized a new Torin 2 derivative, NCATS-SM3710 with increased aqueous solubility and specificity for Plasmodium and demonstrate potent in vivo activity against all P. berghei life cycle stages. NCATS-SM3710 also has low nanomolar EC₅₀s against in vitro cultured asexual P. falciparum parasites (0.38 ± 0.04 nM) and late stage gametocytes (5.77 ± 1 nM). Two independent NCATS-SM3710/Torin 2 resistant P. falciparum parasite lines produced by growth in sublethal Torin 2 concentrations both had genetic changes in PF3D7_0509800, annotated as a phosphatidylinositol 4 kinase (Pf PI4KIIIβ). One line had a point mutation in the putative active site (V1357G), and the other line had a duplication of a locus containing Pf PI4KIIIβ. Both lines were also resistant to other Pf PI4K inhibitors. In addition NCATS-SM3710 inhibited purified Pf PI4KIIIβ with an IC₅₀ of 2.0 ± 0.30 nM. Together the results demonstrate that Pf PI4KIIIβ is the target of Torin 2 and NCATS-SM3710 and provide new options for potent multistage drug development.

Antimalarial pantothenamide metabolites target acetyl-coenzyme A biosynthesis in Plasmodium falciparum

Malaria eradication is critically dependent on new therapeutics that target resistant Plasmodium parasites and block transmission of the disease. Here, we report that pantothenamide bioisosteres were active against blood-stage Plasmodium falciparum parasites and also blocked transmission of sexual stages to the mosquito vector. These compounds were resistant to degradation by serum pantetheinases, showed favorable pharmacokinetic properties, and cleared parasites in a humanized mouse model of P. falciparum infection. Metabolomics revealed that coenzyme A biosynthetic enzymes converted pantothenamides into coenzyme A analogs that interfered with parasite acetyl-coenzyme A anabolism. Resistant parasites generated in vitro showed mutations in acetyl-coenzyme A synthetase and acyl-coenzyme A synthetase 11. Introduction and reversion of these mutations in P. falciparum using CRISPR-Cas9 gene editing confirmed the roles of these enzymes in the sensitivity of the malaria parasites to pantothenamides. These pantothenamide compounds with a new mode of action may have potential as drugs against malaria parasites.

Detection of simple and complex de novo mutations with multiple reference sequences

The characterization of de novo mutations in regions of high sequence and structural diversity from whole-genome sequencing data remains highly challenging. Complex structural variants tend to arise in regions of high repetitiveness and low complexity, challenging both de novo assembly, in which short reads do not capture the long-range context required for resolution, and mapping approaches, in which improper alignment of reads to a reference genome that is highly diverged from that of the sample can lead to false or partial calls. Long-read technologies can potentially solve such problems but are currently unfeasible to use at scale. Here we present Corticall, a graph-based method that combines the advantages of multiple technologies and prior data sources to detect arbitrary classes of genetic variant. We construct multisample, colored de Bruijn graphs from short-read data for all samples, align long-read–derived haplotypes and multiple reference data sources to restore graph connectivity information, and call variants using graph path-finding algorithms and a model for simultaneous alignment and recombination. We validate and evaluate the approach using extensive simulations and use it to characterize the rate and spectrum of de novo mutation events in 119 progeny from four Plasmodium falciparum experimental crosses, using long-read data on the parents to inform reconstructions of the progeny and to detect several known and novel nonallelic homologous recombination events.