Reduced impact of pyrimethamine drug pressure on Plasmodium malariae dihydrofolate reductase gene

Geographical emergence of sulfadoxine-pyrimethamine drug resistance-associated P. falciparum and P. malariae alleles in co-existing Anopheles mosquito and asymptomatic human populations across Cameroon

Malaria molecular surveillance remains critical in detecting and tracking emerging parasite resistance to anti-malarial drugs. The current study employed molecular techniques to determine Plasmodium species prevalence and characterize the genetic diversity of Plasmodium falciparum and Plasmodium malariae molecular markers of sulfadoxine-pyrimethamine resistance in humans and wild Anopheles mosquito populations in Cameroon. Anopheles mosquito collections and parasitological survey were conducted in villages to determine Plasmodium species infection, and genomic phenotyping of anti-folate resistance was accomplished by sequencing the dihydrofolate-reductase (dhfr) and dihydropteroate-synthase (dhps) genes of naturally circulating P. falciparum and P. malariae isolates. The malaria prevalence in Elende was 73.5% with the 5-15 years age group harboring significant P. falciparum (27%) and P. falciparum + P. malariae (19%) infections. The polymorphism breadth of the pyrimethamine-associated Pfdhfr marker revealed a near fixation (94%) of the triple-mutant -A16I51R59N108I164. The Pfdhps backbone mediating sulfadoxine resistance reveals a high frequency of the V431A436G437K540A581A613 alleles (20.8%). Similarly, the Pmdhfr N50K55L57R58S59S114F168I170 haplotype (78.4%) was predominantly detected in the asexual blood stage. In contrast, the Pmdhps- S436A437occured at 37.2% frequency. The combined quadruple N50K55L57R58S59S114F168I170_ S436G437K540A581A613 (31.9%) was the major circulating haplotype with similar frequency in humans and mosquitoes. This study highlights the increasing frequency of the P. malariae parasite mostly common in asymptomatic individuals with apparent P. falciparum infection. Interventions directed at reducing malaria transmission such as the scaling-up of SP are favoring the emergence and spread of multiple drug-resistant alleles between the human and mosquito host systems.

Limited Polymorphism in the Dihydrofolate Reductase (dhfr) and dihydropteroate synthase genes (dhps) of Plasmodium knowlesi isolate from Thailand

BACKGROUND

The 2022 malaria WHO reported around 4000 P. knowlesi infections in the South-East Asia region. In the same period, 72 positive cases were reported by the Department of Disease Control in Thailand, suggesting a persistent infection. Little is known about dihydrofolate reductase (pkdhfr) and dihydropteroate synthase (pkdhps), putative antimalarial resistance markers for P. knowlesi. The relevant amplification and sequencing protocol are presently unavailable. In this study, we developed a protocol for amplifying and evaluating pkdhps mutations. The haplotype pattern of pkdhfr-pkdhps in Thai isolates was analyzed, and the effects of these pkdhps mutations were predicted by using a computer program.

METHODS

Pkdhps were amplified and sequenced from 28 P. knowlesi samples collected in 2008 and 2020 from nine provinces across Thailand. Combining pkdhfr sequencing data from previous work with pkdhps data to analyze polymorphisms of pkdhfr and pkdhps haplotype. Protein modeling and molecular docking were constructed using two inhibitors, sulfadoxine and sulfamethoxazole, and further details were obtained through analyses of protein-ligand interactions by using the Genetic Optimisation for Ligand Docking program. A phylogenetic tree cluster analysis was reconstructed to compare the P. knowlesi Malaysia isolates.

RESULTS

Five nonsynonymous mutations in the pkdhps were detected outside the equivalence of the binding pocket sites to sulfadoxine and sulfamethoxazole, which are at N391S, E421G, I425R, A449S, and N517S. Based on the modeling and molecular docking analyses, the N391S and N517S mutations located close to the enzyme-binding pocket demonstrated a different docking score and protein-ligand interaction in loop 2 of the enzyme. These findings indicated that it was less likely to induce drug resistance. Of the four haplotypes of pkdhfr-pkdhps, the most common one is the R34L pkdhfr mutation and the pkdhps quadruple mutation (GRSS) at E421G, I425R, A449S, and N517S, which were observed in P. knowlesi in southern Thailand (53.57%). Based on the results of neighbor-joining analysis for pkdhfr and pkdhps, the samples isolated from eastern Thailand displayed a close relationship with Cambodia isolates, while southern Thailand isolates showed a long branch separated from the Malaysian isolates.

CONCLUSIONS

A new PCR protocol amplification and evaluation of dihydropteroate synthase mutations in Knowlesi (pkdhps) has been developed. The most prevalent pkdhfr-pkdhps haplotypes (53.57%) in southern Thailand are R34L pkdhfr mutation and pkdhps quadruple mutation. Further investigation requires additional phenotypic data from clinical isolates, transgenic lines expressing mutant alleles, or recombinant proteins.

Plasmodium malariae structure and genetic diversity in sub-Saharan Africa determined from microsatellite variants and linked SNPs in orthologues of antimalarial resistance genes

Plasmodium malariae, a neglected human malaria parasite, contributes up to 10% of malaria infections in sub-Saharan Africa (sSA). Though P. malariae infection is considered clinically benign, it presents mostly as coinfections with the dominant P. falciparum. Completion of its reference genome has paved the way to further understand its biology and interactions with the human host, including responses to antimalarial interventions. We characterized 75 P. malariae isolates from seven endemic countries in sSA using highly divergent microsatellites. The P. malariae infections were highly diverse and five subpopulations from three ancestries (independent of origin of isolates) were determined. Sequences of 11 orthologous antimalarial resistance genes, identified low frequency single nucleotide polymorphisms (SNPs), strong linkage disequilibrium between loci that may be due to antimalarial drug selection. At least three sub-populations were detectable from a subset of denoised SNP data from mostly the mitochondrial cytochrome b coding region. This evidence of diversity and selection calls for including P. malariae in malaria genomic surveillance towards improved tools and strategies for malaria elimination.

Disparate selection of mutations in the dihydrofolate reductase gene (dhfr) of Plasmodium ovale curtisi and P. o. wallikeri in Africa

Plasmodium ovale curtisi and P. ovale wallikeri are both endemic in sub-Saharan Africa, the Middle East and Southeast Asia. Molecular surveillance data for drug resistance in P. ovale spp. is limited at present. We analysed polymorphisms in the podhfr, pocrt and pocytb genes of P. ovale spp. in 147 samples collected from travelers returning to China from Africa. Two podhfr mutations, S58R and S113N/T were detected in P. ovale curtisi with high/moderate frequencies of 52.17% and 17.39%, respectively. Evidence of positive selection (dN/dS = 2.41) was found for podhfr in P. ovale curtisi and decreased diversity (He) of microsatellite markers flanking the mutant alleles suggests that selective sweeps have occurred for both. Mutations E34G (1.50%) and L43V (1.50%) in pocrt of P. ovale curtisi, and E34G (3.70%), I102M (1.80%) and V111F (1.80%) of P. ovale wallikeri were found at low frequencies. Mutations R66K (6.20%), R75K (11.63%) and R95K (3.88%) of pocytb were found in both P. ovale curtisi and P. ovale wallikeri. These results suggest that the podhfr gene of P. ovale curtisi may be subject to drug selection in Africa, warranting further attention. We observed significant differences in the prevalence and distribution of podhfr mutations between the two P. ovale species, suggestive of fundamental biological differences between them.

The primate malaria parasites Plasmodium malariae, Plasmodium brasilianum and Plasmodium ovale spp.: genomic insights into distribution, dispersal and host transitions

During the twentieth century, there was an explosion in understanding of the malaria parasites infecting humans and wild primates. This was built on three main data sources: from detailed descriptive morphology, from observational histories of induced infections in captive primates, syphilis patients, prison inmates and volunteers, and from clinical and epidemiological studies in the field. All three were wholly dependent on parasitological information from blood-film microscopy, and The Primate Malarias” by Coatney and colleagues (1971) provides an overview of this knowledge available at that time. Here, 50 years on, a perspective from the third decade of the twenty-first century is presented on two pairs of primate malaria parasite species. Included is a near-exhaustive summary of the recent and current geographical distribution for each of these four species, and of the underlying molecular and genomic evidence for each. The important role of host transitions in the radiation of Plasmodium spp. is discussed, as are any implications for the desired elimination of all malaria species in human populations. Two important questions are posed, requiring further work on these often ignored taxa. Is Plasmodium brasilianum, circulating among wild simian hosts in the Americas, a distinct species from Plasmodium malariae? Can new insights into the genomic differences between Plasmodium ovale curtisi and Plasmodium ovale wallikeri be linked to any important differences in parasite morphology, cell biology or clinical and epidemiological features?

Plasmodium malariae, current knowledge and future research opportunities on a neglected malaria parasite species

Plasmodium malariae is often reported as a benign malaria parasite. There are limited data on its biology and disease burden in sub-Saharan Africa (sSA) possibly due to the unavailability of specific and affordable tools for routine diagnosis and large epidemiology studies. In addition, P. malariae occurs at low parasite densities and in co-infections with other species, predominately P. falciparum. The paucity of data on P. malariae infections limits the capacity to accurately determine its contribution to malaria and the effect of control interventions against P. falciparum on its prevalence. Here, we summarise the current knowledge on P. malariae epidemiology in sSA - overall prevalence ranging from 0-32%, as detected by different diagnostic methods; seroprevalence ranging from 0-56% in three countries (Mozambique, Benin and Zimbabwe), and explore the future application of next-generation sequencing technologies as a tool for enriching P. malariae genomic epidemiology. This will provide insights into important adaptive mechanisms of this neglected non-falciparum species, including antimalarial drug resistance, local and regional parasite transmission patterns and genomic signatures of selection. Improved diagnosis and genomic surveillance of non-falciparum malaria parasites in Africa would be helpful in evaluating progress towards elimination of all human Plasmodium species.

Genetic analysis of the orthologous crt and mdr1 genes in Plasmodium malariae from Thailand and Myanmar

Background

Plasmodium malariae is a widely spread but neglected human malaria parasite, which causes chronic infections. Studies on genetic polymorphisms of anti-malarial drug target genes in P. malariae are limited. Previous reports have shown polymorphisms in the P. malariae dihydrofolate reductase gene associated with pyrimethamine resistance and linked to pyrimethamine drug pressure. This study investigated polymorphisms of the P. malariae homologous genes, chloroquine resistant transporter and multidrug resistant 1, associated with chloroquine and mefloquine resistance in Plasmodium falciparum.

Methods

The orthologous P. malariae crt and mdr1 genes were studied in 95 patients with P. malariae infection between 2002 and 2016 from Thailand (N = 51) and Myanmar (N = 44). Gene sequences were analysed using BioEdit, MEGA7, and DnaSP programs. Mutations and gene amplifications were compared with P. falciparum and Plasmodium vivax orthologous genes. Protein topology models derived from the observed pmcrt and pmmdr1 haplotypes were constructed and analysed using Phyre2, SWISS MODEL and Discovery Studio Visualization V 17.2.

Results

Two non-synonymous mutations were observed in exon 2 (H53P, 40%) and exon 8 (E278D, 44%) of pmcrt. The topology model indicated that H53P and E278D were located outside of the transmembrane domain and were unlikely to affect protein function. Pmmdr1 was more diverse than pmcrt, with 10 non-synonymous and 3 synonymous mutations observed. Non-synonymous mutations were located in the parasite cytoplasmic site, transmembrane 11 and nucleotide binding domains 1 and 2. Polymorphisms conferring amino acid changes in the transmembrane and nucleotide binding domains were predicted to have some effect on PmMDR1 conformation, but were unlikely to affect protein function. All P. malariae parasites in this study contained a single copy of the mdr1 gene.

Conclusions

The observed polymorphisms in pmcrt and pmmdr1 genes are unlikely to affect protein function and unlikely related to chloroquine drug pressure. Similarly, the absence of pmmdr1 copy number variation suggests limited mefloquine drug pressure on the P. malariae parasite population, despite its long time use in Thailand for the treatment of falciparum malaria.

Indigenous Plasmodium malariae Infection in an Endemic Population at the Thai-Myanmar Border

Plasmodium malariae is a neglected malaria parasite. It has wide geographic distribution and, although often associated with mild malaria, is linked to a high burden of anemia and nephrotic syndromes. Here, we report a cohort study conducted in the Kanchanaburi Province of Thailand during May 2013-June 2014 in which P. malariae infection was detected. Of the 812 study participants, two were found to be infected with P. malariae. One had an infection that led to acute malaria, but the other was positive for P. malariae at multiple visits during the study and apparently had chronic asymptomatic infection. Such persistent infection may explain how P. malariae has been able to thrive at very low prevalence and represents a challenge for malaria elimination.

Plasmodium malariae and Plasmodium ovale infections in the China-Myanmar border area

Background

The Greater Mekong Subregion is aiming to achieve regional malaria elimination by 2030. Though a shift in malaria parasite species predominance by Plasmodium vivax has been recently documented, the transmission of the two minor Plasmodium species, Plasmodium malariae and Plasmodium ovale spp., is poorly characterized in the region. This study aims to determine the prevalence of these minor species in the China–Myanmar border area and their genetic diversity.

Methods

Epidemiology study was conducted during passive case detection in hospitals and clinics in Myanmar and four counties in China along the China–Myanmar border. Cross-sectional surveys were conducted in villages and camps for internally displaced persons to determine the prevalence of malaria infections. Malaria infections were diagnosed initially by microscopy and later in the laboratory using nested PCR for the SSU rRNA genes. Plasmodium malariae and P. ovale infections were confirmed by sequencing the PCR products. The P. ovale subtypes were determined by sequencing the Pocytb, Pocox1 and Pog3p genes. Parasite populations were evaluated by PCR amplification and sequencing of the MSP-1 genes. Antifolate sensitivity was assessed by sequencing the dhfr-ts and dhps genes from the P. malariae and P. ovale isolates.

Results

Analysis of 2701 blood samples collected from the China–Myanmar border by nested PCR targeting the parasite SSU rRNA genes identified 561 malaria cases, including 161 Plasmodium falciparum, 327 P. vivax, 66 P. falciparum/P. vivax mixed infections, 4 P. malariae and 3 P. ovale spp. P. vivax and P. falciparum accounted for >60 and ~30% of all malaria cases, respectively. In comparison, the prevalence of P. malariae and P. ovale spp. was very low and only made up ~1% of all PCR-positive cases. Nevertheless, these two species were often misidentified as P. vivax infections or completely missed by microscopy even among symptomatic patients. Phylogenetic analysis of the SSU rRNA, Pocytb, Pocox1 and Pog3p genes confirmed that the three P. ovale spp. isolates belonged to the subtype P. ovale curtisi. Low-level genetic diversity was detected in the MSP-1, dhfr and dhps genes of these minor parasite species, potentially stemming from the low prevalence of these parasites preventing their mixing. Whereas most of the dhfr and dhps positions equivalent to those conferring antifolate resistance in P. falciparum and P. vivax were wild type, a new mutation S113C corresponding to the S108 position in pfdhfr was identified in two P. ovale curtisi isolates.

Conclusions

The four human malaria parasite species all occurred sympatrically at the China–Myanmar border. While P. vivax has become the predominant species, the two minor parasite species also occurred at very low prevalence but were often misidentified or missed by conventional microscopy. These minor parasite species displayed low levels of polymorphisms in the msp-1, dhfr and dhps genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1605-y) contains supplementary material, which is available to authorized users.

Genetic variability of Plasmodium malariae dihydropteroate synthase (dhps) in four Asian countries

The dihydropteroate synthase (dhps) genes of 44 P. malariae strains from four Asian countries were isolated. Only a limited number of polymorphisms were observed. Comparison with homologous mutations in other Plasmodium species showed that these polymorphisms are unlikely to be associated with sulfadoxine resistance.

Accessible mutational trajectories for the evolution of pyrimethamine resistance in the malaria parasite Plasmodium vivax

Antifolate antimalarials, such as pyrimethamine, have experienced a dramatic reduction in therapeutic efficacy as resistance has evolved in multiple malaria species. We present evidence from one such species, Plasmodium vivax, which has experienced sustained selection for pyrimethamine resistance at the dihydrofolate reductase (DHFR) locus since the 1970s. Using a transgenic Saccharomyces cerevisiae model expressing the P. vivax DHFR enzyme, we assayed growth rate and resistance of all 16 combinations of four DHFR amino acid substitutions. These substitutions were selected based on their known association with drug resistance, both in natural isolates and in laboratory settings, in the related malaria species P. falciparum. We observed a strong correlation between the resistance phenotypes for these 16 P. vivax alleles and previously observed resistance data for P. falciparum, which was surprising since nucleotide diversity levels and common polymorphic variants of DHFR differ between the two species. Similar results were observed when we expressed the P. vivax alleles in a transgenic bacterial system. This suggests common constraints on enzyme evolution in the orthologous DHFR proteins. The interplay of negative trade-offs between the evolution of novel resistance and compromised endogenous function varies at different drug dosages, and so too do the major trajectories for DHFR evolution. In simulations, it is only at very high drug dosages that the most resistant quadruple mutant DHFR allele is favored by selection. This is in agreement with common polymorphic DHFR data in P. vivax, from which this quadruple mutant is missing. We propose that clinical dosages of pyrimethamine may have historically been too low to select for the most resistant allele, or that the fitness cost of the most resistant allele was untenable without a compensatory mutation elsewhere in the genome.