Selective sweep and phylogenetic models for the emergence and spread of pyrimethamine resistance mutations in Plasmodium vivax

Population Genomics of Plasmodium malariae from Four African Countries

Plasmodium malariae is geographically widespread but neglected and may become more prevalent as P. falciparum declines. We completed the largest genomic study of African P. malariae to-date by performing hybrid capture and sequencing of 77 isolates from Cameroon (n=7), the Democratic Republic of the Congo (n=16), Nigeria (n=4), and Tanzania (n=50) collected between 2015 and 2021. There is no evidence of geographic population structure. Nucleotide diversity was significantly lower than in co-localized P. falciparum isolates, while linkage disequilibrium was significantly higher. Genome-wide selection scans identified no erythrocyte invasion ligands or antimalarial resistance orthologs as top hits; however, targeted analyses of these loci revealed evidence of selective sweeps around four erythrocyte invasion ligands and six antimalarial resistance orthologs. Demographic inference modeling suggests that African P. malariae is recovering from a bottleneck. Altogether, these results suggest that P. malariae is genomically atypical among human Plasmodium spp. and panmictic in Africa.

Establishment and evaluation of a qPCR method for the detection of pfmdr1 mutations in Plasmodium falciparum, the causal agent of fatal malaria

Drug resistance surveillance is a major integral part of malaria control programs. Molecular methods play a pivotal role in drug resistance detection and related molecular research. This study aimed to develop a rapid and accurate detection method for drug resistance of Plasmodium falciparum (P. falciparum). A quantitative real-time PCR (qPCR) assay has been developed that identifies the mutation at locus A256T in the P.falciparum multi-drug resistance(pfmdr1) gene producing amino acid change at position 86. The results of 198 samples detected by qPCR were consistent with nested PCR and sequencing, giving an accuracy of 94.3%. The sensitivity, specificity, positive and negative predictive value of qPCR were 85.7%, 97.6%, 90.0% and 96.4%, respectively. The results of qPCR are basically consistent with the nested PCR, which is expected to replace the nested PCR as a new molecular biological method for drug resistance detection, providing reliable technical support for global malaria prevention and control.

Antimalarial Drug Resistance: A Brief History of Its Spread in Indonesia

Malaria remains to be a national and global challenge and priority, as stated in the strategic plan of the Indonesian Ministry of Health and Sustainable Development Goals. In Indonesia, it is targeted that malaria elimination can be achieved by 2030. Unfortunately, the development and spread of antimalarial resistance inflicts a significant risk to the national malaria control programs which can lead to increased malaria morbidity and mortality. In Indonesia, resistance to widely used antimalarial drugs has been reported in two human species, Plasmodium falciparum and Plasmodium vivax. With the exception of artemisinin, resistance has surfaced towards all classes of antimalarial drugs. Initially, chloroquine, sulfadoxine-pyrimethamine, and primaquine were the most widely used antimalarial drugs. Regrettably, improper use has supported the robust spread of their resistance. Chloroquine resistance was first reported in 1974, while sulfadoxine-pyrimethamine emerged in 1979. Twenty years later, most provinces had declared treatment failures of both drugs. Molecular epidemiology suggested that variations in pfmdr1 and pfcrt genes were associated with chloroquine resistance, while dhfr and dhps genes were correlated with sulfadoxine-pyrimethamine resistance. Additionally, G453W, V454C and E455K of pfk13 genes appeared to be early warning sign to artemisinin resistance. Here, we reported mechanisms of antimalarial drugs and their development of resistance. This insight could provide awareness toward designing future treatment guidelines and control programs in Indonesia.

Genetic characterisation of the Theileria annulata cytochrome b locus and its impact on buparvaquone resistance in bovine

Control of tropical theileriosis, caused by the apicomplexan Theileria annulata, depends on the use of a single drug, buparvaquone, the efficacy of which is compromised by the emergence of resistance. The present study was undertaken to improve understanding of the role of mutations conferring buparvaquone resistance in T. annulata, and the effects of selection pressures on their emergence and spread. First, we investigated genetic characteristics of the cytochrome b locus associated with buparvaquone resistance in 10 susceptible and 7 resistant T. annulata isolates. The 129G (GGC) mutation was found in the Q₀₁ binding pocket and 253S (TCT) and 262S (TCA) mutations were identified within the Q₀₂ binding pocket. Next, we examined field isolates and identified cytochrome b mutations 129G (GGC), 253S (TCT) and 262S (TCA) in 21/75 buffalo-derived and 19/119 cattle-derived T. annulata isolates, providing evidence of positive selection pressure. Both hard and soft selective sweeps were identified, with striking differences between isolates. For example, 19 buffalo-derived and 7 cattle-derived isolates contained 129G (GGC) and 253S (TCT) resistance haplotypes at a high frequency, implying the emergence of resistance by a single mutation. Two buffalo-derived and 12 cattle-derived isolates contained equally high frequencies of 129G (GGC), 253S (TCT), 129G (GGC)/253S (TCT) and 262S (TCA) resistance haplotypes, implying the emergence of resistance by pre-existing or recurrent mutations. Phylogenetic analysis further revealed that 9 and 21 unique haplotypes in buffalo and cattle-derived isolates were present in a single lineage, suggesting a single origin. We propose that animal migration between farms is an important factor in the spread of buparvaquone resistance in endemic regions of Pakistan. The overall outcomes will be useful in understanding how drug resistance emerges and spreads, and this information will help design strategies to optimise the use and lifespan of the single most drug use to control tropical theileriosis.

Prevalence of antifolate drug resistance markers in Plasmodium vivax in China

The dihydrofolate reductase (dhfr) and dihydropteroate synthetase (dhps) genes of Plasmodium vivax, as antifolate resistance-associated genes were used for drug resistance surveillance. A total of 375 P. vivax isolates collected from different geographical locations in China in 2009-2019 were used to sequence Pvdhfr and Pvdhps. The majority of the isolates harbored a mutant type allele for Pvdhfr (94.5%) and Pvdhps (68.2%). The most predominant point mutations were S117T/N (77.7%) in Pvdhfr and A383G (66.8%) in Pvdhps. Amino acid changes were identified at nine residues in Pvdhfr. A quadruple-mutant haplotype at 57, 58, 61, and 117 was the most frequent (57.4%) among 16 distinct Pvdhfr haplotypes. Mutations in Pvdhps were detected at six codons, and the double-mutant A383G/A553G was the most prevalent (39.3%). Pvdhfr exhibited a higher mutation prevalence and greater diversity than Pvdhps in China. Most isolates from Yunnan carried multiple mutant haplotypes, while the majority of samples from temperate regions and Hainan Island harbored the wild type or single mutant type. This study indicated that the antifolate resistance levels of P. vivax parasites were different across China and molecular markers could be used to rapidly monitor drug resistance. Results provided evidence for updating national drug policy and treatment guidelines.

Current methods for the detection of antimalarial drug resistance in Plasmodium parasites infecting humans

Malaria is the world's deadliest parasitic disease. Great progress has been made in the fight against malaria over the past two decades, but this has recently begun to plateau, in part due to the global development of antimalarial drug resistance. The ability to track drug resistance is necessary to achieve progress in treatment, disease surveillance and epidemiology, which has prompted the development of advanced diagnostic methods. These new methods provide unprecedented access to information that can help to guide public health policies. Development of new technologies increases the potential for high throughput and reduced costs of diagnostic tests; improving the accessibility of tools to investigate the forces driving disease dynamics and, ultimately, clinical outcomes for malaria patients and public health. This literature review provides a summary of the methods currently available for the detection of antimalarial drug resistance from the examination of patients' blood samples. While no single method is perfect for every application, many of the newly developed methods give promise for more reliable and efficient characterisation of Plasmodium resistance in a range of settings. By exploiting the strengths of the tools available, we can develop a deeper understanding of the evolutionary and spatiotemporal dynamics of this disease. This will translate into more effective disease control, better-informed policy, and more timely and successful treatment for malaria patients.

Phylogenetic analysis suggests single and multiple origins of dihydrofolate reductase mutations in Plasmodium vivax

Pyrimethamine was first introduced for the treatment of malaria in Asia and Africa during the early 1980s, replacing chloroquine, and has become the first line of drugs in many countries. In recent years, development of pyrimethamine resistance in Plasmodium vivax has become a barrier to effective malaria control strategies. Here, we describe the use of meta-barcoded deep amplicon sequencing technology to assess the evolutionary origin of pyrimethamine resistance by analysing the flanking region of dihydrofolate reductase (dhfr) locus. The genetic modelling suggests that 58R and 173L single mutants and 58R/117N double mutants are present on a single lineage; suggesting a single origin of these mutations. The triple mutants (57L/58R/117N, 58R/61M/117N and 58R/117N/173L) share the lineage of 58R/117N, suggesting a common origin. In contrast, the 117N mutant is present on two separate lineages suggesting that there are multiple origins of this mutation. We characterised the allele frequency of the P. vivax dhfr locus. Our results support the view that the single mutation of 117N and double mutations of 58R/117N arise commonly, whereas the single mutation of 173L and triple mutations of 57L/58R/117N, 58R/61M/117N and 58R/117N/173L are less common. Our work will help to inform mitigation strategies for pyrimethamine resistance in P. vivax.

A novel metabarcoded 18S ribosomal DNA sequencing tool for the detection of Plasmodium species in malaria positive patients

Various PCR based methods have been described for the diagnosis of malaria, but most depend on the use of Plasmodium species-specific probes and primers; hence only the tested species are identified and there is limited available data on the true circulating species diversity. Sensitive diagnostic tools and platforms for their use are needed to detect Plasmodium species in both clinical cases and asymptomatic infections that contribute to disease transmission. We have recently developed for the first time a novel high throughput 'haemoprotobiome' metabarcoded DNA sequencing method and applied it for the quantification of haemoprotozoan parasites (Theleria and Babesia) of livestock. Here, we describe a novel, high throughput method using an Illumina MiSeq platform to demonstrate the proportions of Plasmodium species in metabarcoded DNA samples derived from human malaria patients. Plasmodium falciparum and Plasmodium vivax positive control gDNA was used to prepare mock DNA pools of parasites to evaluate the detection threshold of the assay for each of the two species. The different mock pools demonstrate the accurate detection ability and to show the proportions of each of the species being present. We then applied the assay to malaria-positive human samples to show the species composition of Plasmodium communities in the Punjab province of Pakistan and in the Afghanistan-Pakistan tribal areas. The diagnostic performance of the deep amplicon sequencing method was compared to an immunochromatographic assay that is widely used in the region. The deep amplicon sequencing showed that P. vivax was present in 69.8%, P. falciparum in 29.5% and mixed infection in 0.7% patients examined. The immunochromatographic assay showed that P. vivax was present in 65.6%, P. falciparum in 27.4%, mixed infection 0.7% patients and 6.32% malaria-positive cases were negative in immunochromatographic assay, but positive in the deep amplicon sequencing. Overall, metabarcoded DNA sequencing demonstrates better diagnostic performance, greatly increasing the estimated prevalence of Plasmodium infection. The next-generation sequencing method using metabarcoded DNA has potential applications in the diagnosis, surveillance, treatment, and control of Plasmodium infections, as well as to study the parasite biology.