Longitudinal Pooled Deep Sequencing of the Plasmodium vivax K12 Kelch Gene in Cambodia Reveals a Lack of Selection by Artemisinin

High frequency of artemisinin partial resistance mutations in the great lake region revealed through rapid pooled deep sequencing

In Africa, the first Plasmodium falciparum Kelch13 (K13) artemisinin partial resistance mutation 561H was first detected and validated in Rwanda. Surveillance to better define the extent of the emergence in Rwanda and neighboring countries as other mutations arise in East Africa is critical. We employ a novel scheme of liquid blood drop preservation combined with pooled sequencing to provide a cost-effective rapid assessment of resistance mutation frequencies at multiple collection sites across Rwanda and neighboring countries. Malaria-positive samples (n=5,465) were collected from 39 health facilities in Rwanda, Uganda, Tanzania, and the Democratic Republic of the Congo (DRC) between May 2022 and March 2023 and sequenced in 199 pools. In Rwanda, K13 561H and 675V were detected in 90% and 65% of sites with an average frequency of 19.0% (0-54.5%) and 5.0% (0-35.5%), respectively. In Tanzania, 561H had high frequency in multiple sites while it was absent from the DRC although 675V was seen at low frequency. Conceringly candidate mutations were observed: 441L, 449A, and 469F co-occurred with validated mutations suggesting they are arising under the same pressures. Other resistance markers associated with artemether-lumefantrine are common: P. falciparum multidrug resistance protein 1 N86 at 98.0% and 184F at 47.0% (0-94.3%) and P. falciparum chloroquine resistance transporter 76T at 14.7% (0-58.6%). Additionally, sulfadoxine-pyrimethamine-associated mutations show high frequencies. Overall, K13 mutations are rapidly expanding in the region further endangering control efforts with the potential of engendering partner drug resistance.

Drug resistance of Plasmodium falciparum and Plasmodium vivax isolates in Indonesia

This review article aims to investigate the genotypic profiles of Plasmodium falciparum and Plasmodium vivax isolates collected across a wide geographic region and their association with resistance to anti-malarial drugs used in Indonesia. A systematic review was conducted between 1991 and date. Search engines, such as PubMed, Science Direct, and Google Scholar, were used for articles published in English and Indonesian to search the literature. Of the 471 initially identified studies, 61 were selected for 4316 P. falciparum and 1950 P. vivax individual infections. The studies included 23 molecular studies and 38 therapeutic efficacy studies. K76T was the most common pfcrt mutation. K76N (2.1%) was associated with the haplotype CVMNN. By following dihydroartemisinin-piperaquine (DHA-PPQ) therapy, the mutant pfmdr1 alleles 86Y and 1034C were selected. Low prevalence of haplotype N86Y/Y184/D1246Y pfmdr1 reduces susceptibility to AS-AQ. SNP mutation pvmdr1 Y976F reached 96.1% in Papua and East Nusa Tenggara. Polymorphism analysis in the pfdhfr gene revealed 94/111 (84.7%) double mutants S108N/C59R or S108T/A16V in Central Java. The predominant pfdhfr haplotypes (based on alleles 16, 51, 59,108, 164) found in Indonesia were ANCNI, ANCSI, ANRNI, and ANRNL. Some isolates carried A437G (35.3%) or A437G/K540E SNPs (26.5%) in pfdhps. Two novel pfdhps mutant alleles, I588F/G and K540T, were associated with six pfdhps haplotypes. The highest prevalence of pvdhfr quadruple mutation (F57L/S58R/T61M/S117T) (61.8%) was detected in Papua. In pvdhps, the only polymorphism before and after 2008 was 383G mutation with 19% prevalence. There were no mutations in the pfk13 gene reported with validated and candidate or associated k13 mutation. An increased copy number of pfpm2, associated with piperaquine resistance, was found only in cases of reinfection. Meanwhile, mutation of pvk12 and pvpm4 I165V is unlikely associated with ART and PPQ drug resistance. DHA-PPQ is still effective in treating uncomplicated falciparum and vivax malaria. Serious consideration should be given to interrupt local malaria transmission and dynamic patterns of resistance to anti-malarial drugs to modify chemotherapeutic policy treatment strategies. The presence of several changes in pfk13 in the parasite population is of concern and highlights the importance of further evaluation of parasitic ART susceptibility in Indonesia.

Limited polymorphism in k13 gene of Plasmodium falciparum and k12 of Plasmodium vivax isolates imported from African and Asian countries between 2014 and 2019 in Hangzhou city, China

Background

Malaria causes major public health problems globally and drug resistance hinders its control and elimination. Molecular markers associated with drug resistance are considered as a beneficial tool to monitor the disease trends, evolution and distribution so as to help improve drug policy.

Methods

We collected 148 Plasmodium falciparum and 20 Plasmodium vivax isolates imported into Hangzhou city, China between 2014 and 2019. k13 gene of P. falciparum and k12 of P. vivax were sequenced. Polymorphisms and prevalence of k13 and k12 were analyzed.

Results

Most (98.65%, 146/148) P. falciparum infections were imported from Africa, and half P. vivax cases came from Africa and the other half from Asia. Nucleotide mutation prevalence was 2.03% (3/148) and the proportion of amino acid mutations was 0.68% (1/148). The amino acid mutation, A676S, was observed in an isolate from Nigeria. No mutation of k12 was observed from the parasites from African and Asian countries.

Conclusions

Limited polymorphism in k13 gene of P. falciparum isolates imported from African countries, but no evidence for the polymorphism of k12 in P. vivax samples from African and Asian countries was found. These results provide information for drug policy update in study region.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06579-6.

The molecular basis of antimalarial drug resistance in Plasmodium vivax

Plasmodium vivax is the most geographically widespread cause of human malaria and is responsible for the majority of cases outside of the African continent. While great progress has been made towards eliminating human malaria, drug resistant parasite strains pose a threat towards continued progress. Resistance has arisen to multiple antimalarials in P. vivax, including to chloroquine, which is currently the first line therapy for P. vivax in most regions. Despite its importance, an understanding of the molecular mechanisms of drug resistance in this species remains elusive, in large part due to the complex biology of P. vivax and the lack of in vitro culture. In this review, we will cover the extent and challenges of measuring clinical and in vitro drug resistance in P. vivax. We will consider the roles of candidate drug resistance genes. We will highlight the development of molecular approaches for studying P. vivax biology that provide the opportunity to validate the role of putative drug resistance mutations as well as identify novel mechanisms of drug resistance in this understudied parasite. Validated molecular determinants and markers of drug resistance are essential for the rapid and cost-effective monitoring of drug resistance in P. vivax, and will be useful for optimizing drug regimens and for informing drug policy in control and elimination settings.

•

Drug resistance is emerging in Plasmodium vivax, an important cause of malaria.

•

The complex biology of P. vivax and the limited range of research tools make it difficult to identify drug resistance.

•

The molecular mechanisms of drug resistance in P. vivax remain elusive.

•

This review highlights the extent of drug resistance, the putative mechanisms of resistance and new technologies for the study of P. vivax drug resistance.

Selection of Cytochrome b Mutants Is Rare among Plasmodium falciparum Patients Failing Treatment with Atovaquone-Proguanil in Cambodia

Atovaquone-proguanil remains effective against multidrug-resistant Plasmodium falciparum in Southeast Asia, but resistance is mediated by a single point mutation in cytochrome b (cytb) that can arise during treatment. Among 14 atovaquone-proguanil treatment failures in a clinical trial in Cambodia, only one recrudescence harbored the cytb mutation Y268C. Deep sequencing did not detect the mutation at baseline or in the first 3 days of treatment, suggesting that it arose de novo Further sequencing across cytb similarly found no low-frequency cytb mutations that were up-selected from baseline to recrudescence. Copy number amplification in dihydroorotate dehydrogenase (DHODH) and cytb as markers of atovaquone tolerance was also absent. Cytb mutation played a minor role in atovaquone-proguanil treatment failures in an active comparator clinical trial.

Polymorphisms in Plasmodium falciparum Kelch 13 and P. vivax Kelch 12 Genes in Parasites Collected from Three South Pacific Countries Prior to Extensive Exposure to Artemisinin Combination Therapies

The South Pacific countries Solomon Islands, Vanuatu, and Papua New Guinea (PNG) adopted artemisinin-based combination therapies (ACTs) in 2008. We examined Kelch 13 and Kelch 12 genes in parasites originating from these countries before or at ACT introduction. Four Kelch 13 and two Kelch 12 novel sequence polymorphisms, not associated with artemisinin resistance, were observed in parasites from Solomon Islands and Vanuatu. No polymorphisms were observed in PNG parasites. The findings provide useful baseline information.

Polymorphisms in Pvkelch12 and gene amplification of Pvplasmepsin4 in Plasmodium vivax from Thailand, Lao PDR and Cambodia

BACKGROUND

Mutations in Pfkelch13 and Pfplasmepsin2/3 gene amplification are well-established markers for artemisinin and piperaquine resistance in Plasmodium falciparum, a widespread problem in the Greater Mekong Subregion (GMS). The Plasmodium vivax parasite population has experienced varying drug pressure dependent on local drug policies. We investigated the correlation between drug pressure from artemisinins and piperaquine and mutations in the P. vivax orthologous genes Pvkelch12 and Pvplasmepsin4 (Pvpm4), as candidate resistance markers.

METHODS

Blood samples from 734 P. vivax patients were obtained from Thailand (n = 399), Lao PDR (n = 296) and Cambodia (n = 39) between 2007 and 2017. Pvkelch12 and Pvpm4 was amplified and sequenced to assess gene mutations. To assess PvPM4 gene amplification, a Taqman^® Real-Time PCR method was developed and validated. Selection of non-synonymous mutations was assessed by its ratio with synonymous mutations (Ka/Ks ratios). Mutation rates were compared to the estimated local drug pressure.

RESULTS

Polymorphisms in Pvkelch12 were rare. Pvkelch12 mutations V552I, K151Q and M124I were observed in 1.0% (7/734) of P. vivax samples. V552I was the most common mutation with a frequency of 0.7% (5/734), most of which (4/5) observed in Ubon Ratchathani, Thailand. Polymorphisms in Pvpm4 were more common, with a frequency of 40.3% (123/305) in 305 samples from Thailand, Lao PDR and Cambodia, but this was not related to the estimated piperaquine drug pressure in these areas (Pearson's χ² test, p = 0.50). Pvpm4 mutation V165I was most frequent in Tak, Thailand (40.2%, 43/107) followed by Pailin, Cambodia (43.5%, 37/85), Champasak, Lao PDR (40.4%, 23/57) and Ubon Ratchathani, Thailand (35.7%, 20/56). Pvpm4 amplification was not observed in 141 samples from Thailand and Cambodia. For both Pvkelch12 and Pvpm4, in all areas and at all time points, the Ka/Ks values were < 1, suggesting no purifying selection.

CONCLUSIONS

A novel real-time PCR-based method to assess P. vivax Pvpm4 gene amplification was developed. Drug pressure with artemisinins and piperaquine in the GMS was not clearly related to signatures of selection for mutations in the P. vivax orthologous resistance genes Pvkelch12 and Pvpm4 in areas under investigation. Current resistance of P. vivax to these drugs is unlikely and additional observations including analysis of associated clinical data from these regions could further clarify current findings.

Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics

During the last decade, the vast omics field has revolutionized biological research, especially the genomics, transcriptomics and proteomics branches, as technological tools become available to the field researcher and allow difficult question-driven studies to be addressed. Parasitology has greatly benefited from next generation sequencing (NGS) projects, which have resulted in a broadened comprehension of basic parasite molecular biology, ecology and epidemiology. Malariology is one example where application of this technology has greatly contributed to a better understanding of Plasmodium spp. biology and host-parasite interactions. Among the several parasite species that cause human malaria, the neglected Plasmodium vivax presents great research challenges, as in vitro culturing is not yet feasible and functional assays are heavily limited. Therefore, there are gaps in our P. vivax biology knowledge that affect decisions for control policies aiming to eradicate vivax malaria in the near future. In this review, we provide a snapshot of key discoveries already achieved in P. vivax sequencing projects, focusing on developments, hurdles, and limitations currently faced by the research community, as well as perspectives on future vivax malaria research.

A deep sequencing approach to estimate Plasmodium falciparum complexity of infection (COI) and explore apical membrane antigen 1 diversity

BACKGROUND

Humans living in regions with high falciparum malaria transmission intensity harbour multi-strain infections comprised of several genetically distinct malaria haplotypes. The number of distinct malaria parasite haplotypes identified from an infected human host at a given time is referred to as the complexity of infection (COI). In this study, an amplicon-based deep sequencing method targeting the Plasmodium falciparum apical membrane antigen 1 (pfama1) was utilized to (1) investigate the relationship between P. falciparum prevalence and COI, (2) to explore the population genetic structure of P. falciparum parasites from malaria asymptomatic individuals participating in the 2007 Demographic and Health Survey (DHS) in the Democratic Republic of Congo (DRC), and (3) to explore selection pressures on geospatially divergent parasite populations by comparing AMA1 amino acid frequencies in the DRC and Mali.

RESULTS

A total of 900 P. falciparum infections across 11 DRC provinces were examined. Deep sequencing of both individuals, for COI analysis, and pools of individuals, to examine population structure, identified 77 unique pfama1 haplotypes. The majority of individual infections (64.5%) contained polyclonal (COI > 1) malaria infections based on the presence of genetically distinct pfama1 haplotypes. A minimal correlation between COI and malaria prevalence as determined by sensitive real-time PCR was identified. Population genetic analyses revealed extensive haplotype diversity, the vast majority of which was shared across the sites. AMA1 amino acid frequencies were similar between parasite populations in the DRC and Mali.

CONCLUSIONS

Amplicon-based deep sequencing is a useful tool for the detection of multi-strain infections that can aid in the understanding of antigen heterogeneity of potential malaria vaccine candidates, population genetics of malaria parasites, and factors that influence complex, polyclonal malaria infections. While AMA1 and other diverse markers under balancing selection may perform well for understanding COI, they may offer little geographic or temporal discrimination between parasite populations.

Surveillance for sulfadoxine-pyrimethamine resistant malaria parasites in the Lake and Southern Zones, Tanzania, using pooling and next-generation sequencing

Background

Malaria in pregnancy (MiP) remains a major public health challenge in areas of high malaria transmission. Intermittent preventive treatment in pregnancy (IPTp) with sulfadoxine-pyrimethamine (SP) is recommended to prevent the adverse consequences of MiP. The effectiveness of SP for IPTp may be reduced in areas where the dhps581 mutation (a key marker of high level SP resistance) is found; this mutation was previously reported to be common in the Tanga Region of northern Tanzania, but there are limited data from other areas. The frequency of molecular markers of SP resistance was investigated in malaria parasites from febrile patients at health centres (HC) in seven regions comprising the Lake and Southern Zones of mainland Tanzania as part of the ongoing efforts to generate national-wide data of SP resistance.

Methods

A cross-sectional survey was conducted in the outpatient departments of 14 HCs in seven regions from April to June, 2015. 1750 dried blood spot (DBS) samples were collected (117 to 160 per facility) from consenting patients with positive rapid diagnostic tests for malaria, and no recent (within past 2 months) exposure to SP or related drugs. DNA was extracted from the DBS, pooled by HC, and underwent pooled targeted amplicon deep sequencing to yield estimates of mutated parasite allele frequency at each locus of interest.

Results

The dhps540 mutation was common across all 14 sites, ranging from 55 to 98.4% of sequences obtained. Frequency of the dhps581 mutation ranged from 0 to 2.4%, except at Kayanga HC (Kagera Region, Lake Zone) where 24.9% of sequences obtained were mutated. The dhfr164 mutation was detected only at Kanyanga HC (0.06%).

Conclusion

By pooling DNA extracts, the allele frequency of mutations in 14 sites could be directly determined on a single deep-sequencing run. The dhps540 mutant was very common at all locations. Surprisingly, the dhps581 was common at one health center, but rare in all the others, suggesting that there is geographic micro-heterogeneity in mutant distribution and that accurate surveillance requires inclusion of multiple sites. A better understanding of the effect of the dhps581 mutant on the efficacy of IPTp-SP is needed.

Divergent pattern of genomic variation in Plasmodium falciparum and P. vivax

The two main species causing malaria in humans, Plasmodium falciparum and P. vivax, differ significantly from each other in their evolutionary response to common drugs, but the reasons for this are not clear. Here we utilized the recently available large-scale genome sequencing data from these parasites and compared the pattern of single nucleotide polymorphisms, which may be related to these differences. We found that there was a five-fold higher preference for AT nucleotides compared to GC nucleotides at synonymous single nucleotide polymorphism sites in P. vivax. The preference for AT nucleotides was also present at non-synonymous sites, which lead to amino acid changes favouring those with codons of higher AT content. The substitution bias was also present at low and moderately conserved amino acid positions, but not at highly conserved positions. No marked bias was found at synonymous and non-synonymous sites in P. falciparum. The difference in the substitution bias between P. falciparum and P. vivax found in the present study may possibly contribute to their divergent evolutionary response to similar drug pressures.