Characterization of Plasmodium falciparum structure in Nigeria with malaria SNPs barcode

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

Molecular markers for malaria genetic epidemiology: progress and pitfalls

Over recent years, progress in molecular markers for genotyping malaria parasites has enabled informative studies of epidemiology and transmission dynamics. Results have highlighted the value of these tools for surveillance to support malaria control and elimination strategies. There are many different types and panels of markers available for malaria parasite genotyping, and for end users, the nuances of these markers with respect to 'use case', resolution, and accuracy, are not well defined. This review clarifies issues surrounding different molecular markers and their application to malaria control and elimination. We describe available marker panels, use cases, implications for different transmission settings, limitations, access, cost, and data accuracy. The information provided can be used as a guide for molecular epidemiology and surveillance of malaria.

Baseline prevalence of molecular marker of sulfadoxine/pyrimethamine resistance in Ebonyi and Osun states, Nigeria: amplicon deep sequencing of dhps-540

BACKGROUND

Chemoprevention plays an important role in malaria control strategy. Perennial malaria chemoprevention (PMC) using sulfadoxine/pyrimethamine (SP) is a WHO-approved strategy to combat malaria in young children and may lead to drug pressure. Introducing SP-PMC may therefore be compromised due to the emergence of Plasmodium falciparum resistant to SP, particularly mutation at K540E of the dihydropteroate synthase (dhps) gene. Molecular surveillance of resistance markers can support assessment of antimalarial efficacy and effectiveness. High prevalence of 540E is associated with reduced effectiveness of SP, and areas with more than 50% prevalence are considered unsuitable for intermittent preventative treatment in pregnancy (IPTp) implementation. Assessing 540E prevalence is an important undertaking before implementation of SP-PMC.

METHODS

We conducted a rapid surveillance of dhps-540E to assess the suitability of SP as PMC in field studies from Ebonyi and Osun states in Nigeria. We used an in-house developed amplicon deep-sequencing method targeting part of the dhps gene.

RESULTS

Our data reveal that 18.56% of individuals evaluated carried the 540E mutation mixed with the WT K540. Mutant variant 540E alone was not found, and 80% of isolates harboured only WT (K540). Clonal analysis of the sequencing data shows a very low proportion of 540E circulating in both states.

CONCLUSIONS

Our data show that both states are suitable for SP-PMC implementation and, based on this finding, SP-PMC was implemented in Osun in 2022. Continuous monitoring of 540E will be required to ensure the chemoprevention effectiveness of SP in Nigeria.

CRISPR-Cas13 in malaria parasite: Diagnosis and prospective gene function identification

Malaria caused by Plasmodium is still a serious public health problem. Genomic editing is essential to understand parasite biology, elucidate mechanical pathways, uncover gene functions, identify novel therapeutic targets, and develop clinical diagnostic tools. Recent advances have seen the development of genomic diagnostic technologies and the emergence of genetic manipulation toolbox comprising a host of several systems for editing the genome of Plasmodium at the DNA, RNA, and protein level. Genomic manipulation at the RNA level is critical as it allows for the functional characterization of several transcripts. Of notice, some developed artificial RNA genome editing tools hinge on the endogenous RNA interference system of Plasmodium. However, Plasmodium lacks a robust RNAi machinery, hampering the progress of these editing tools. CRISPR-Cas13, which belongs to the VI type of the CRISPR system, can specifically bind and cut RNA under the guidance of crRNA, with no or minimal permanent genetic scar on genes. This review summarizes CRISPR-Cas13 system from its discovery, classification, principle of action, and diagnostic platforms. Further, it discusses the application prospects of Cas13-based systems in Plasmodium and highlights its advantages and drawbacks.

The many definitions of multiplicity of infection

The presence of multiple genetically different pathogenic variants within the same individual host is common in infectious diseases. Although this is neglected in some diseases, it is well recognized in others like malaria, where it is typically referred to as multiplicity of infection (MOI) or complexity of infection (COI). In malaria, with the advent of molecular surveillance, data is increasingly being available with enough resolution to capture MOI and integrate it into molecular surveillance strategies. The distribution of MOI on the population level scales with transmission intensities, while MOI on the individual level is a confounding factor when monitoring haplotypes of particular interests, e.g., those associated with drug-resistance. Particularly, in high-transmission areas, MOI leads to a discrepancy between the likelihood of a haplotype being observed in an infection (prevalence) and its abundance in the pathogen population (frequency). Despite its importance, MOI is not universally defined. Competing definitions vary from verbal ones to those based on concise statistical frameworks. Heuristic approaches to MOI are popular, although they do not mine the full potential of available data and are typically biased, potentially leading to misinferences. We introduce a formal statistical framework and suggest a concise definition of MOI and its distribution on the host-population level. We show how it relates to alternative definitions such as the number of distinct haplotypes within an infection or the maximum number of alleles detectable across a set of genetic markers. It is shown how alternatives can be derived from the general framework. Different statistical methods to estimate the distribution of MOI and pathogenic variants at the population level are discussed. The estimates can be used as plug-ins to reconstruct the most probable MOI of an infection and set of infecting haplotypes in individual infections. Furthermore, the relation between prevalence of pathogenic variants and their frequency (relative abundance) in the pathogen population in the context of MOI is clarified, with particular regard to seasonality in transmission intensities. The framework introduced here helps to guide the correct interpretation of results emerging from different definitions of MOI. Especially, it excels comparisons between studies based on different analytical methods.

Bias-corrected maximum-likelihood estimation of multiplicity of infection and lineage frequencies

Background

The UN’s Sustainable Development Goals are devoted to eradicate a range of infectious diseases to achieve global well-being. These efforts require monitoring disease transmission at a level that differentiates between pathogen variants at the genetic/molecular level. In fact, the advantages of genetic (molecular) measures like multiplicity of infection (MOI) over traditional metrics, e.g., R₀, are being increasingly recognized. MOI refers to the presence of multiple pathogen variants within an infection due to multiple infective contacts. Maximum-likelihood (ML) methods have been proposed to derive MOI and pathogen-lineage frequencies from molecular data. However, these methods are biased.

Methods and findings

Based on a single molecular marker, we derive a bias-corrected ML estimator for MOI and pathogen-lineage frequencies. We further improve these estimators by heuristical adjustments that compensate shortcomings in the derivation of the bias correction, which implicitly assumes that data lies in the interior of the observational space. The finite sample properties of the different variants of the bias-corrected estimators are investigated by a systematic simulation study. In particular, we investigate the performance of the estimator in terms of bias, variance, and robustness against model violations. The corrections successfully remove bias except for extreme parameters that likely yield uninformative data, which cannot sustain accurate parameter estimation. Heuristic adjustments further improve the bias correction, particularly for small sample sizes. The bias corrections also reduce the estimators’ variances, which coincide with the Cramér-Rao lower bound. The estimators are reasonably robust against model violations.

Conclusions

Applying bias corrections can substantially improve the quality of MOI estimates, particularly in areas of low as well as areas of high transmission—in both cases estimates tend to be biased. The bias-corrected estimators are (almost) unbiased and their variance coincides with the Cramér-Rao lower bound, suggesting that no further improvements are possible unless additional information is provided. Additional information can be obtained by combining data from several molecular markers, or by including information that allows stratifying the data into heterogeneous groups.

A genotyping assay to determine geographic origin and transmission potential of Plasmodium falciparum malaria cases

As countries work towards malaria elimination, it is important to monitor imported cases to prevent reestablishment of local transmission. The Plasmodium falciparum Pfs47 gene has strong geographic population structure, because only those parasites with Pfs47 haplotypes compatible with the mosquito vector species in a given continent are efficiently transmitted. Analysis of 4,971 world-wide Pfs47 sequences identified two SNPs (at 707 and 725 bp) as sufficient to establish the likely continent of origin of P. falciparum isolates. Pfs47 sequences from Africa, Asia, and the New World presented more that 99% frequency of distinct combinations of the SNPs 707 and 725 genotypes. Interestingly, Papua New Guinea Pfs47 sequences have the highest diversity in SNPs 707 and 725. Accurate and reproducible High-Resolution Melting (HRM) assays were developed to genotype Pfs47 SNPs 707 and 725 in laboratory and field samples, to assess the geographic origin and risk of local transmission of imported P. falciparum malaria cases.

Genetic diversity and population structure of Plasmodium falciparum in Nigeria: insights from microsatellite loci analysis

BACKGROUND

Malaria remains a public health burden especially in Nigeria. To develop new malaria control and elimination strategies or refine existing ones, understanding parasite population diversity and transmission patterns is crucial.

METHODS

In this study, characterization of the parasite diversity and structure of Plasmodium falciparum isolates from 633 dried blood spot samples in Nigeria was carried out using 12 microsatellite loci of P. falciparum. These microsatellite loci were amplified via semi-nested polymerase chain reaction (PCR) and fragments were analysed using population genetic tools.

RESULTS

Estimates of parasite genetic diversity, such as mean number of different alleles (13.52), effective alleles (7.13), allelic richness (11.15) and expected heterozygosity (0.804), were high. Overall linkage disequilibrium was weak (0.006, P < 0.001). Parasite population structure was low (Fst: 0.008-0.105, AMOVA: 0.039).

CONCLUSION

The high level of parasite genetic diversity and low population structuring in this study suggests that parasite populations circulating in Nigeria are homogenous. However, higher resolution methods, such as the 24 SNP barcode and whole genome sequencing, may capture more specific parasite genetic signatures circulating in the country. The results obtained can be used as a baseline for parasite genetic diversity and structure, aiding in the formulation of appropriate therapeutic and control strategies in Nigeria.

Polymorphisms in Plasmodium falciparum chloroquine resistance transporter (Pfcrt) and multidrug-resistant gene 1 (Pfmdr-1) in Nigerian children 10 years post-adoption of artemisinin-based combination treatments

The emergence and spread of Plasmodium falciparum parasites resistant to artemisinin derivatives and their partners in southeastern Asia threatens malaria control and elimination efforts, and heightens the need for an alternative therapy. We have explored the distribution of P. falciparum chloroquine resistance transporter (Pfcrt) and multidrug-resistant gene 1 (Pfmdr-1) haplotypes 10 years following adoption of artemisinin-based combination therapies in a bid to investigate the possible re-emergence of Chloroquine-sensitive parasites in Nigeria, and investigated the effect of these P. falciparum haplotypes on treatment outcomes of patients treated with artemisinin-based combination therapies. A total of 271 children aged <5 years with uncomplicated falciparum malaria were included in this study. Polymorphisms on codons 72-76 of the Pfcrt gene and codon 86 and 184 of Pfmdr-1 were determined using the high resolution melting assay. Of 240 (88.6%) samples successfully genotyped with HRM for Pfcrt, wildtype C72M74N75K76 (42.9%) and mutant C72I74E75T76 (53.8%) were observed. Also, wildtype N86Y184 (62.9%) and mutant N86F184 (21.1%), Y86Y184 (6.4%), and Y86F184 (0.4%) haplotypes of Pfmdr-1 were observed. Measures of responsiveness to ACTs were similar in children infected with P. falciparum crt haplotypes (C72I74E75T76 and C72M74N75K76) and major mdr-1 haplotypes (N86Y184, N86F184 and Y86Y184). Despite a 10 year gap since the malaria treatment policy changed to ACTs, over 50% of the P. falciparum parasites investigated in this study harboured the Chloroquine-resistant C72I74E75T76 haplotype, however this did not compromise the efficacy of artemisinin-based combination therapies. Should complete artemisinin resistance emerge from or spread to Nigeria, chloroquine might not be a good alternative therapy.

Polymorphisms in Plasmodium falciparum dihydropteroate synthetase and dihydrofolate reductase genes in Nigerian children with uncomplicated malaria using high-resolution melting technique

In 2005, the Nigerian Federal Ministry of Health revised the treatment policy for uncomplicated malaria with the introduction of artemisinin-based combination therapies (ACTs). This policy change discouraged the use of Sulphadoxine-pyrimethamine (SP) as the second-line treatment of uncomplicated falciparum malaria. However, SP is used as an intermittent preventive treatment of malaria in pregnancy (IPTp) and seasonal malaria chemoprevention (SMC) in children aged 3–59 months. There have been increasing reports of SP resistance especially in the non-pregnant population in Nigeria, thus, the need to continually monitor the efficacy of SP as IPTp and SMC by estimating polymorphisms in dihydropteroate synthetase (dhps) and dihydrofolate reductase (dhfr) genes associated with SP resistance. The high resolution-melting (HRM) assay was used to investigate polymorphisms in codons 51, 59, 108 and 164 of the dhfr gene and codons 437, 540, 581 and 613 of the dhps gene. DNA was extracted from 271 dried bloodspot filter paper samples obtained from children (< 5 years old) with uncomplicated malaria. The dhfr triple mutant I₅₁R₅₉N₁₀₈, dhps double mutant G₄₃₇G₅₈₁ and quadruple dhfr I₅₁R₅₉N₁₀₈ + dhps G₄₃₇ mutant haplotypes were observed in 80.8%, 13.7% and 52.8% parasites, respectively. Although the quintuple dhfr I₅₁R₅₉N₁₀₈ + dhps G₄₃₇E₅₄₀ and sextuple dhfr I₅₁R₅₉N₁₀₈ + dhps G₄₃₇E₅₄₀G₅₈₁ mutant haplotypes linked with in-vivo and in-vitro SP resistance were not detected, constant surveillance of these haplotypes should be done in the country to detect any change in prevalence.

Post-Genomic Approaches to Understanding Malaria Parasite Biology: Linking Genes to Biological Functions

Plasmodium species are evolutionarily distant from model eukaryotes, and as a consequence they exhibit many non-canonical cellular processes. In the post-genomic era, functional "omics" disciplines (transcriptomics, proteomics, and metabolomics) have accelerated our understanding of unique aspects of the biology of malaria parasites. Functional "omics" tools, in combination with genetic manipulations, have offered new opportunities to investigate the function of previously uncharacterized genes. Knowledge of basic parasite biology is fundamental to understanding drug modes of action, mechanisms of drug resistance, and relevance of vaccine candidates. This Perspective highlights recent "omics"-based discoveries in basic biology and gene function of the most virulent human malaria parasite, Plasmodium falciparum.