How to combat emerging artemisinin resistance: Lessons from "The Three Little Pigs".

It is rare to come across an Aesop's fable in respectable journals. It might catch scientists outside the malaria field by surprise to learn that the famous story of "The Boy Who Cried Wolf" has been repeatedly compared to the threat from artemisinin-resistant malaria parasites, including the two latest reports on the rise of a specific haplotype in Cambodia and Thailand, sensationally dubbed "Super Malaria" by the media [1, 2]. The comparison to a children's tale should not negate the fact that malaria drug resistance is one of the most pressing threats to the global public health community. Here, the findings leading to this contentious discourse will be delineated in order to provide a perspective. Possible solutions will be presented to stimulate further research and discussion to solve one of the greatest public health challenges of our lifetime.

Plasmodium falciparum Falcipain-2a Polymorphisms in Southeast Asia and Their Association With Artemisinin Resistance.

Background: Falcipain-2a ([FP2a] PF3D7_1115700) is a Plasmodium falciparum cysteine protease and hemoglobinase. Functional FP2a is required for potent activity of artemisinin, and in vitro selection for artemisinin resistance selected for an FP2a nonsense mutation. Methods: To investigate associations between FP2a polymorphisms and artemisinin resistance and to characterize the diversity of the enzyme in parasites from the China-Myanmar border, we sequenced the full-length FP2a gene in 140 P falciparum isolates collected during 2004-2011. Results: The isolates were grouped into 8 different haplotype groups. Haplotype group I appeared in samples obtained after 2008, coinciding with implementation of artemisinin-based combination therapy in this region. In functional studies, compared with wild-type parasites, the FP2a haplotypes demonstrated increased ring survival, and all haplotype groups exhibited significantly reduced FP2a activity, with group I showing the slowest protease kinetics and reduced parasite fitness. Conclusions: These results suggest that altered hemoglobin digestion due to FP2a mutations may contribute to artemisinin resistance.

Combinatorial Genetic Modeling of pfcrt-Mediated Drug Resistance Evolution in Plasmodium falciparum.

The emergence of drug resistance continuously threatens global control of infectious diseases, including malaria caused by the protozoan parasite Plasmodium falciparum A critical parasite determinant is the P. falciparum chloroquine resistance transporter (PfCRT), the primary mediator of chloroquine (CQ) resistance (CQR), and a pleiotropic modulator of susceptibility to several first-line artemisinin-based combination therapy partner drugs. Aside from the validated CQR molecular marker K76T, P. falciparum parasites have acquired at least three additional pfcrt mutations, whose contributions to resistance and fitness have been heretofore unclear. Focusing on the quadruple-mutant Ecuadorian PfCRT haplotype Ecu1110 (K76T/A220S/N326D/I356L), we genetically modified the pfcrt locus of isogenic, asexual blood stage P. falciparum parasites using zinc-finger nucleases, producing all possible combinations of intermediate pfcrt alleles. Our analysis included the related quintuple-mutant PfCRT haplotype 7G8 (Ecu1110 + C72S) that is widespread throughout South America and the Western Pacific. Drug susceptibilities and in vitro growth profiles of our combinatorial pfcrt-modified parasites were used to simulate the mutational trajectories accessible to parasites as they evolved CQR. Our results uncover unique contributions to parasite drug resistance and growth for mutations beyond K76T and predict critical roles for the CQ metabolite monodesethyl-CQ and the related quinoline-type drug amodiaquine in driving mutant pfcrt evolution. Modeling outputs further highlight the influence of parasite proliferation rates alongside gains in drug resistance in dictating successful trajectories. Our findings suggest that P. falciparum parasites have navigated constrained pfcrt adaptive landscapes by means of probabilistically rare mutational bursts that led to the infrequent emergence of pfcrt alleles in the field.

Comparative genome-wide analysis and evolutionary history of haemoglobin-processing and haem detoxification enzymes in malarial parasites.

BACKGROUND: Malaria parasites have evolved a series of intricate mechanisms to survive and propagate within host red blood cells. Intra-erythrocytic parasitism requires these organisms to digest haemoglobin and detoxify iron-bound haem. These tasks are executed by haemoglobin-specific proteases and haem biocrystallization factors that are components of a large multi-subunit complex. Since haemoglobin processing machineries are functionally and genetically linked to the modes of action and resistance mechanisms of several anti-malarial drugs, an understanding of their evolutionary history is important for drug development and drug resistance prevention. METHODS: Maximum likelihood trees of genetic repertoires encoding haemoglobin processing machineries within Plasmodium species, and with the representatives of Apicomplexan species with various host tropisms, were created. Genetic variants were mapped onto existing three-dimensional structures. Genome-wide single nucleotide polymorphism data were used to analyse the selective pressure and the effect of these mutations at the structural level. RESULTS: Recent expansions in the falcipain and plasmepsin repertoires are unique to human malaria parasites especially in the Plasmodium falciparum and P. reichenowi lineage. Expansion of haemoglobin-specific plasmepsins occurred after the separation event of Plasmodium species, but the other members of the plasmepsin family were evolutionarily conserved with one copy for each sub-group in every Apicomplexan species. Haemoglobin-specific falcipains are separated from invasion-related falcipain, and their expansions within one specific locus arose independently in both P. falciparum and P. vivax lineages. Gene conversion between P. falciparum falcipain 2A and 2B was observed in artemisinin-resistant strains. Comparison between the numbers of non-synonymous and synonymous mutations suggests a strong selective pressure at falcipain and plasmepsin genes. The locations of amino acid changes from non-synonymous mutations mapped onto protein structures revealed clusters of amino acid residues in close proximity or near the active sites of proteases. CONCLUSION: A high degree of polymorphism at the haemoglobin processing genes implicates an imposition of selective pressure. The identification in recent years of functional redundancy of haemoglobin-specific proteases makes them less appealing as potential drug targets, but their expansions, especially in the human malaria parasite lineages, unequivocally point toward their functional significance during the independent and repetitive adaptation events in malaria parasite evolutionary history.

Origin of robustness in generating drug-resistant malaria parasites.

Biological robustness allows mutations to accumulate while maintaining functional phenotypes. Despite its crucial role in evolutionary processes, the mechanistic details of how robustness originates remain elusive. Using an evolutionary trajectory analysis approach, we demonstrate how robustness evolved in malaria parasites under selective pressure from an antimalarial drug inhibiting the folate synthesis pathway. A series of four nonsynonymous amino acid substitutions at the targeted enzyme, dihydrofolate reductase (DHFR), render the parasites highly resistant to the antifolate drug pyrimethamine. Nevertheless, the stepwise gain of these four dhfr mutations results in tradeoffs between pyrimethamine resistance and parasite fitness. Here, we report the epistatic interaction between dhfr mutations and amplification of the gene encoding the first upstream enzyme in the folate pathway, GTP cyclohydrolase I (GCH1). gch1 amplification confers low level pyrimethamine resistance and would thus be selected for by pyrimethamine treatment. Interestingly, the gch1 amplification can then be co-opted by the parasites because it reduces the cost of acquiring drug-resistant dhfr mutations downstream in the same metabolic pathway. The compensation of compromised fitness by extra GCH1 is an example of how robustness can evolve in a system and thus expand the accessibility of evolutionary trajectories leading toward highly resistant alleles. The evolution of robustness during the gain of drug-resistant mutations has broad implications for both the development of new drugs and molecular surveillance for resistance to existing drugs.

Biochemical and functional characterization of Plasmodium falciparum GTP cyclohydrolase I.

BACKGROUND: Antifolates are currently in clinical use for malaria preventive therapy and treatment. The drugs kill the parasites by targeting the enzymes in the de novo folate pathway. The use of antifolates has now been limited by the spread of drug-resistant mutations. GTP cyclohydrolase I (GCH1) is the first and the rate-limiting enzyme in the folate pathway. The amplification of the gch1 gene found in certain Plasmodium falciparum isolates can cause antifolate resistance and influence the course of antifolate resistance evolution. These findings showed the importance of P. falciparum GCH1 in drug resistance intervention. However, little is known about P. falciparum GCH1 in terms of kinetic parameters and functional assays, precluding the opportunity to obtain the key information on its catalytic reaction and to eventually develop this enzyme as a drug target. METHODS: Plasmodium falciparum GCH1 was cloned and expressed in bacteria. Enzymatic activity was determined by the measurement of fluorescent converted neopterin with assay validation by using mutant and GTP analogue. The genetic complementation study was performed in ∆folE bacteria to functionally identify the residues and domains of P. falciparum GCH1 required for its enzymatic activity. Plasmodial GCH1 sequences were aligned and structurally modeled to reveal conserved catalytic residues. RESULTS: Kinetic parameters and optimal conditions for enzymatic reactions were determined by the fluorescence-based assay. The inhibitor test against P. falciparum GCH1 is now possible as indicated by the inhibitory effect by 8-oxo-GTP. Genetic complementation was proven to be a convenient method to study the function of P. falciparum GCH1. A series of domain truncations revealed that the conserved core domain of GCH1 is responsible for its enzymatic activity. Homology modelling fits P. falciparum GCH1 into the classic Tunnelling-fold structure with well-conserved catalytic residues at the active site. CONCLUSIONS: Functional assays for P. falciparum GCH1 based on enzymatic activity and genetic complementation were successfully developed. The assays in combination with a homology model characterized the enzymatic activity of P. falciparum GCH1 and the importance of its key amino acid residues. The potential to use the assay for inhibitor screening was validated by 8-oxo-GTP, a known GTP analogue inhibitor.

Sideways: road to gene-by-gene functional screening in malaria parasites.

Genome-wide screening in apicomplexan species has transformed our understanding of these parasitic protozoa. Kimmel et al. report a 'knock sideways' system and provide a powerful use case for its feasibility in a gene-by-gene screening in Plasmodium falciparum. Carefully deployed, a novel toolkit helps to dissect the biological uniqueness of an important parasite.

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.

Pf7: an open dataset of Plasmodium falciparum genome variation in 20,000 worldwide samples.

We describe the MalariaGEN Pf7 data resource, the seventh release of Plasmodium falciparum genome variation data from the MalariaGEN network.  It comprises over 20,000 samples from 82 partner studies in 33 countries, including several malaria endemic regions that were previously underrepresented.  For the first time we include dried blood spot samples that were sequenced after selective whole genome amplification, necessitating new methods to genotype copy number variations.  We identify a large number of newly emerging crt mutations in parts of Southeast Asia, and show examples of heterogeneities in patterns of drug resistance within Africa and within the Indian subcontinent.  We describe the profile of variations in the C-terminal of the csp gene and relate this to the sequence used in the RTS,S and R21 malaria vaccines.  Pf7 provides high-quality data on genotype calls for 6 million SNPs and short indels, analysis of large deletions that cause failure of rapid diagnostic tests, and systematic characterisation of six major drug resistance loci, all of which can be freely downloaded from the MalariaGEN website.

Stepwise acquisition of pyrimethamine resistance in the malaria parasite.

The spread of high-level pyrimethamine resistance in Africa threatens to curtail the therapeutic lifetime of antifolate antimalarials. We studied the possible evolutionary pathways in the evolution of pyrimethamine resistance using an approach in which all possible mutational intermediates were created by site-directed mutagenesis and assayed for their level of drug resistance. The coding sequence for dihydrofolate reductase (DHFR) from the malaria parasite Plasmodium falciparum was mutagenized, and tests were carried out in Escherichia coli under conditions in which the endogenous bacterial enzyme was selectively inhibited. We studied 4 key amino acid replacements implicated in pyrimethamine resistance: N51I, C59R, S108N, and I164L. Using empirical estimates of the mutational spectrum in P. falciparum and probabilities of fixation based on the relative levels of resistance, we found that the predicted favored pathways of drug resistance are consistent with those reported in previous kinetic studies, as well as DHFR polymorphisms observed in natural populations. We found that 3 pathways account for nearly 90% of the simulated realizations of the evolution of pyrimethamine resistance. The most frequent pathway (S108N and then C59R, N51I, and I164L) accounts for more than half of the simulated realizations. Our results also suggest an explanation for why I164L is detected in Southeast Asia and South America, but not at significant frequencies in Africa.

SARS-CoV-2 RT-PCR positivity of individuals subsequent to completing quarantine upon entry into a country during a transmission-free period.

BACKGROUND: During the current coronavirus disease 2019 (COVID-19) pandemic, many countries require travellers to undergo a reverse transcription-polymerase chain reaction (RT-PCR) testing for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) before travelling across borders. However, in persons having recovered from COVID-19, RT-PCR positivity can persist for an extended period. MATERIALS AND METHODS: We describe three cases who sought fit-to-fly certificates in Thailand during the period free of local transmission but were tested positive for RT-PCR for SARS-CoV-2. All had returned from a country with an active outbreak of COVID-19. Their clinical courses are described; positive nasopharyngeal swab samples were processed for viral isolation and whole-genome sequencing (WGS); and serology as well as neutralizing antibody were assessed. The contact tracing was carried out for determining evidence of indigenous transmission among close contacts of those three cases. RESULTS: All three cases were completely asymptomatic. Chest computerized tomography was not compatible with COVID-19 pneumonia; cell cultures failed to rescue replication-competent virus; WGS revealed fragmented viral genetic material from nasopharyngeal swab samples; and serological tests demonstrated stable levels of antibodies, together with the presence of neutralizing antibody, suggesting past infection with negligible transmission risk. Contact tracing identified no transmission in high-risk close contact individuals. CONCLUSION: RT-PCR positivity for SARS-CoV-2 might detect fragmented viral genome. Issuance of a travel certificate in these circumstances is problematic. Serology tests can help to define past infection. A practical acceptable set of guidelines for issuance of a COVID-19 safety travel certification is a necessity.

Characterization of band 3-ankyrin-Protein 4.2 complex by biochemical and mass spectrometry approaches.

The elastic property of red blood cell is supported by interaction between red cell membrane and the intricate cytoskeleton network underlying the membrane bilayer cytoplasmic face. One of the major scaffold protein linkers is band 3-ankyrin complex. Defects occurring in this complex have been found in many inherited diseases, causing red blood cell abnormalities. Here we combined the power of mass spectrometry with conventional biochemical purification methods in order to study the native interactions among band 3, ankyrin and Protein 4.2. This approach provided in vivo evidence for the association between band 3 and N-terminal ankyrin purified directly from the cell membrane. The C-terminal regions of ankyrin were not found to be a stable partner of the band 3 complex. Protein 4.2 was shown here to be an integral part of the complex. Its association to the band 3-ankyrin complex could withstand harsh purification conditions. Our findings lend additional support to the interaction between band 3 and ankyrin N-terminal domain previously shown by in vitro binding assays and provide evidence for a band 3 core complex comprising of band 3, ankyrin and Protein 4.2.

Chromatin associated sense and antisense noncoding RNAs are transcribed from the var gene family of virulence genes of the malaria parasite Plasmodium falciparum.

Antigenic variation by the malaria parasite Plasmodium falciparum results from switches in expression between members of the multicopy var gene family. These genes encode the variant surface protein PfEMP-1, the primary determinant of the antigenic and cytoadherent properties of infected erythrocytes. Only a single var gene is expressed at a time while the remaining members of the family remain transcriptionally silent. How mutually exclusive var gene expression is regulated is poorly understood; however, it is generally thought to involve alterations in chromatin assembly and modification, resulting in a type of cellular memory. Recently, several aspects of the chromatin structure surrounding var genes have been described, in particular the histone modifications associated with the active and silent states of the genes as well as their subnuclear localization. Here, we demonstrate that this chromatin structure also includes the incorporation of long sense and antisense noncoding RNAs. These sterile transcripts initiate from a bidirectional promoter located within a conserved intron found in all var genes that was previously implicated in var gene silencing. Mapping of the 5' and 3' ends of the sterile transcripts indicates that they are nonpolyadenylated. RNA fluorescent in situ hybridization (RNA-FISH) analysis detects both the sense and antisense noncoding RNAs in distinct spots within the nucleus similar to the pattern described for the var genes themselves. Further, analysis by RNA chromatin immunoprecipitation (ChIP) indicates that the noncoding RNAs are physically associated with chromatin. These sterile transcripts therefore might act in a manner analogous to noncoding RNAs associated with silent, condensed chromatin found in other eukaryotic systems.

COVID-19 Transmission among Healthcare Workers at a Quarantine Facility in Thailand: Genomic and Outbreak Investigations.

During the COVID-19 pandemic, Thailand implemented a quarantine program at approved quarantine facilities for every international traveler. Here, we report an epidemiological and genomic investigation of a COVID-19 cluster consisting of seven healthcare workers (HCWs) at a quarantine facility and its partnered hospital in Thailand. Outbreak investigations were implemented to obtain contact tracing data and to establish chains of transmission. Genomic sequencing of SARS-CoV-2 with samples within the cohort was performed. Investigations of 951 HCWs and staff with quarantined travelers were implemented to determine the chain of transmission. Genomic and outbreak investigations identified the international travelers infected with the B.1.1.31 SARS-CoV-2 lineage as the source of this outbreak. The genomic data and the investigated timeline revealed a putative transmission chain among HCWs, pointing toward the transmission via the use of common living quarters at the investigated quarantine site. The evaluation of this cohort has led to a policy recommendation on quarantine facility management. International travel quarantine is an important strategy to contain importation of COVID-19 cases. However, a quarantine facility is likely to become a potential hotspot, requiring thorough preventive measures. Reducing the exposure risk by providing private living quarters and scheduling clinical duties at a quarantine site separated from the conventional healthcare workforce have been implemented.

Lumefantrine attenuates Plasmodium falciparum artemisinin resistance during the early ring stage.

Emerging artemisinin resistance in Plasmodium falciparum malaria has the potential to become a global public health crisis. In Southeast Asia, this phenomenon clinically manifests in the form of delayed parasite clearance following artemisinin treatment. Reduced artemisinin susceptibility is limited to the early ring stage window, which is sufficient to allow parasites to survive the short half-life of artemisinin exposure. A screen of known clinically-implemented antimalarial drugs was performed to identify a drug capable of enhancing the killing activity of artemisinins during this critical resistance window. As a result, lumefantrine was found to increase the killing activity of artemisinin against an artemisinin-resistant clinical isolate harboring the C580Y kelch13 mutation. Isobologram analysis revealed synergism during the early ring stage resistance window, when lumefantrine was combined with artemether, an artemisinin derivative clinically partnered with lumefantrine. These findings suggest that lumefantrine should be clinically explored as a partner drug in artemisinin-based combination therapies to control emerging artemisinin resistance.

Genomic surveillance of SARS-CoV-2 in Thailand reveals mixed imported populations, a local lineage expansion and a virus with truncated ORF7a.

Coronavirus Disease 2019 (COVID-19) is a global public health threat. Genomic surveillance of SARS-CoV-2 was implemented in March of 2020 at a major diagnostic hub in Bangkok, Thailand. Several virus lineages supposedly originated in many countries were found, and a Thai-specific lineage, designated A/Thai-1, has expanded to be predominant in Thailand. A virus sample in the SARS-CoV-2 A/Thai-1 lineage contains a frame-shift deletion at ORF7a, encoding a putative host antagonizing factor of the virus.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples.

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

A nuclear targeting system in Plasmodium falciparum.

BACKGROUND: The distinct differences in gene control mechanisms acting in the nucleus between Plasmodium falciparum and the human host could lead to new potential drug targets for anti-malarial development. New molecular toolkits are required for dissecting molecular machineries in the P. falciparum nucleus. One valuable tool commonly used in model organisms is protein targeting to specific sub-cellular locations. Targeting proteins to specified locations allows labeling of organelles for microscopy, or testing of how the protein of interest modulates organelle function. In recent years, this approach has been developed for various malaria organelles, such as the mitochondrion and the apicoplast. A tool for targeting a protein of choice to the P. falciparum nucleus using an exogenous nuclear localization sequence is reported here. METHODS: To develop a nuclear targeting system, a putative nuclear localization sequence was fused with green fluorescent protein (GFP). The nuclear localization sequence from the yeast transcription factor Gal4 was chosen because of its well-defined nuclear localization signal. A series of truncated Gal4 constructs was also created to narrow down the nuclear localization sequence necessary for P. falciparum nuclear import. Transfected parasites were analysed by fluorescent and laser-scanning confocal microscopy. RESULTS: The nuclear localization sequence of Gal4 is functional in P. falciparum. It effectively transported GFP into the nucleus, and the first 74 amino acid residues were sufficient for nuclear localization. CONCLUSIONS: The Gal4 fusion technique enables specific transport of a protein of choice into the P. falciparum nucleus, and thus provides a tool for labeling nuclei without using DNA-staining dyes. The finding also indicates similarities between the nuclear transport mechanisms of yeast and P. falciparum. Since the nuclear transport system has been thoroughly studied in yeast, this could give clues to research on the same mechanism in P. falciparum.

Epigenetic transcriptional gene silencing in Entamoeba histolytica: insight into histone and chromatin modifications.

We have previously discovered a unique mechanism of epigenetic transcriptional gene silencing in the Entamoeba histolytica trophozoites of strain HM-1:IMSS that resulted in the persistent downregulation of the amoebapore A (ap-a) gene, and that could be successfully applied to silence other virulence genes (cpA5, lgl1). In order to understand how the silencing is maintained throughout generations, we analysed whether modifications occurred at the chromatin level. Chromatin immunoprecipitation assays were done with antibodies specific to the methylated lysine 4 of E. histolytica histone H3. When the genes were in a transcriptionally silent state, the methylation levels of H3K4 in their coding region were significantly reduced. In contrast, the levels of core histone H3 were consistently higher in the silenced genes. Controlled chromatin digestion with micrococcal nuclease was used to assess changes in nucleosome compaction. We found a significant resistance to digestion in the promoter region of the silenced ap-a and cpA5 genes as compared to the parental strain that expresses those genes. Our data lend further support to the idea that histone modifications and heterochromatin formations are at the basis of the transcriptional silencing of genes in E. histolytica.